This invention relates, e.g., to the identification of non-essential genes of bacteria, and of a minimal set of genes required to support viability of a free-living organism.
One consequence of progress in the new field of synthetic biology is an emerging view of cells as assemblages of parts that can be put together to produce an organism with a desired phenotype (1). That perspective begs the question: “How few parts would it take to construct a cell?” In an environment that is free from stress and provides all necessary nutrients, what would comprise the simplest free-living organism? This problem has been approached theoretically and experimentally in our laboratory and elsewhere.
In a comparison of the first two bacterial genomes sequenced, Mushegian and Koonin projected that the 256 orthologous genes shared by the Gram negative Haemophilus influenzae and the Gram positive M. genitalium genomes are a close approximation of a minimal gene set for bacterial life (2). More recently Gil et al. proposed a 206 protein-coding gene core of a minimal bacterial gene set based on analysis of several free-living and endosymbiotic bacterial genomes (3).
In 1999 some of the present inventors reported the first use of global transposon mutagenesis to experimentally determine the genes not essential for laboratory growth of M. genitalium (4). Since then there have been numerous other experimental determinations of bacterial essential gene sets using our approach and other methods such as site directed gene knockouts and antisense RNA (5-12). Most of these studies were done with human pathogens, often with the aim of identifying essential genes that might be used as antibiotic targets. Almost all of these organisms contain relatively large genomes that include many paralogous gene families. Disruption or deletion of such genes shows they are non-essential but does not determine if their products perform essential biological functions. It is only through gene essentiality studies of bacteria that have near minimal genomes that we bring empirical verification to the compositions of hypothetical minimal gene sets.
The Mollicutes, generically known as the mycoplasmas, are an excellent experimental platform for experimentally defining a minimal gene set. These wall-less bacteria evolved from more conventional progenitors in the Firmicutes taxon by a process of massive genome reduction. Mycoplasmas are obligate parasites that live in relatively unchanging niches requiring little adaptive capability. M. genitalium, a human urogenital pathogen, is the extreme manifestation of this genomic parsimony, having only 482 protein-coding genes and the smallest genome at ˜580 kb of any known free-living organism capable of being grown in pure culture (13). The bacteria can grow independently on an agar plate free of other living cells. While more conventional bacteria with larger genomes used in gene essentiality studies have on average 26% of their genes in paralogous gene families, M. genitalium has only 6% (Table 1). Thus, with its lack of genomic redundancy and contingencies for different environmental conditions, M. genitalium is already close to being a minimal bacterial cell.
The 1999 report by some of the present inventors on the essential microbial gene for M. genitalium and its closest relative, Mycoplasma pneumoniae, mapped 2200 transposon insertion sites in these two species, and identified 130 putatively non-essential M. genitalium protein-coding genes or M. pneumoniae orthologs of M. genitalium genes. In that report (Hutchison et al. (1999) Science 286, 2165-9), those authors estimated that 265 to 350 of the protein-coding genes of M. genitalium are essential under laboratory growth conditions (4). However proof of gene dispensability requires isolation and characterization of pure clonal populations, which they did not do. In that report, the authors grew Tn4001 transformed cells in mixed pools for several weeks, and then isolated genomic DNA from those mixtures of mutants. They sequenced amplicons from inverse PCRs using that DNA as a template to identify the transposon insertion sites in the mycoplasma genomes. Most of the genes containing transposon insertions encoded either hypothetical proteins or other proteins not expected to be essential. Nonetheless, some of the putatively disrupted genes, such as isoleucyl and tyrosyl-tRNA synthetases (MG345 & MG455), DNA replication gene dnaA (MG469), and DNA polymerase III, subunit alpha (MG261) are thought to perform essential functions. They hypothesized how genes generally thought to be essential might be disrupted: a gene may be tolerant of the transposon insertion and not actually disrupted, cells could contain two copies of a gene, or the gene product may be supplied by other cells in the same mixed pool of mutants.
Disclosed herein is an expanded study in which we have isolated and characterized M. genitalium Tn4001 insertion mutants that were present in individual colonies picked from agar plates. This analysis has provided a new, more thorough, estimate of the number of essential genes in this minimalist bacterium.
a-2i show global transposon mutagenesis of M. genitalium. The locations of transposon insertions from the current study are noted by a Δ below the insertion site on the map. The letters over the Gene Loci (MG###) refer to the functional category of the gene product as listed.
The inventors have identified 100 protein-coding genes that are non-essential for sustaining the growth of an organism, such as a bacterium, in a rich bacterial culture medium, e.g. SP4. Such a culture medium contains all of the salts, growth factors, nutrients etc. required for bacterial growth under laboratory conditions. A minimal set of genes required for sustaining the viability of a free-living organism under laboratory conditions is extrapolated from the identification of these non-essential genes. By a “minimal gene set” is meant the minimal set of genes whose expression allows the viability (e.g., survival, growth, replication, proliferation, etc.) of a free-living organism in a particular rich bacterial medium as discussed above.
The 100 protein-coding genes of M. genitalium that were disrupted in the bacteria and nevertheless retained viability, and are thus dispensable (non-essential) for growth, are listed in Table 2, where they are grouped by their functional roles. The 382 genes that were not disrupted are summarized in Table 3, where they are also grouped by functional roles. These genes form part of a minimal essential gene set. Other genes may also be part of a minimal gene set. At minimum, these other genes include protein-coding genes for ABC transporters for phosphate and/or phosphonate, and certain lipoproteins and/or glycerophosphoryl diester phosphodiesterases; and RNA-encoding genes.
As noted above, the some of the present inventors published a preliminary study in 1999 that reported putative sets of genes that appeared to be either essential or disposable for viability. Table 4 lists genes identified in the present study as being dispensable, but which were not so identified in the 1999 paper. Table 5 lists genes identified in the present study as being required for growth, but which were not so identified in the 1999 paper.
One aspect of the invention is a set of protein-coding genes that provides the information required for replication of a free-living organism under axenic conditions in a rich bacterial culture medium, such as SP4, (e.g., a minimal set of protein-coding genes),
wherein the gene set lacks at least 40 of the 101 protein-coding genes listed in Table 2 (the “lacking genes”), or functional equivalents thereof, wherein at least one of the genes in Table 4 is among the lacking genes;
wherein the set comprises between 350 and 382 of the 382 protein-coding genes listed in Table 3, or functional equivalents thereof, including at least one of the genes in Table 5; and
wherein the set comprises no more than 450 protein-coding genes.
A set of genes that “provides the information” required for replication of a free-living organism can be in any form that can be transcribed (e.g. into mRNA, rRNA or tRNA) and, in the case of protein-encoding sequences, translated into protein, wherein the transcription/translation products provide functions that allow the free-living organism to function.
This set of protein-coding genes is smaller than the complete complement of genes found in M. genitalium (482 genes), the smallest known set of naturally occurring genes in a free-living organism.
A set of protein-coding genes of the invention can lack at least about 55 (e.g. at least about, 70, 80 or 90) of the genes listed in Table 2), and/or it can comprise at least about 360 (e.g. at least about 370 or 380) of the genes listed in Table 3.
A set of the invention can further comprise:
genes encoding an ABC transporter for phosphate import, selected from the group consisting of (a) MG410, MG411 and MG412, and (b) MG289, MG290 and MG291, and functional equivalents thereof; and/or
a lipoprotein-encoding gene selected from the group consisting of MG185 and MG260, and functional equivalents thereof; and/or
a glycerophosphoryl diester phosphodiesterase gene selected from the group consisting of MG293 and MG385, and functional equivalents thereof.
Furthermore, a set of the invention can further comprise the 43 RNA-coding genes of Mycoplasma genitalium, or functional equivalents thereof.
The genes in a set of the invention may constitute a chromosome; and/or may be from M. genitalium.
Another aspect of the invention is a free-living organism that can grow and replicate under axenic conditions in a rich bacterial culture medium (such as SP4), whose set of genes consists of a set of the invention, e.g. a set that comprises at least one gene involved in hydrogen or ethanol production.
Another aspect of the invention is a method for determining the function of a gene, comprising inserting, mutating or removing the gene into/in/from such a free-living organism, and measuring a property of the organism.
Another aspect of the invention is a method of hydrogen or ethanol production, comprising growing a free-living organism of that invention that comprises at least one gene involved in hydrogen or ethanol production, in a suitable medium such that hydrogen or ethanol is produced.
Another aspect of the invention is an effective subset of a set as noted above. An “effective subset,” as used herein, refers to a subset that provides the information required for replication of a free-living organism in a rich bacterial culture medium, such as SP4.
A minimal gene set of the invention has a variety of applications. For example, a minimal gene set of the invention can be introduced into cells of a microorganism, such as a bacterium, which lack a genome or a functional genome (e.g. ghost cells) and used experimentally to investigate requirements for cell growth, protein synthesis, replication or other bacterial functions under varying conditions. One or more of the minimal genes in the ghost cells can be modified or substituted with orthologous genes or genes or substituted with non-orthologous genes that express proteins which perform the same function(s), to allow structure/function studies of those genes. Cells comprising a minimal gene set of the invention can be modified to further comprise one or more expressible heterologous genes, either integrated into the genome or replicating on one or more independent plasmids. These cells can be used, e.g., to study properties or activities of the heterologous genes (e.g., structure/function studies), or to produce useful amounts of the heterologous proteins (e.g. biologic drugs, vaccines, catalytic enzymes, energy sources, etc).
As noted, a minimal gene set is one that provides the information required for replication of a free-living organism in a rich bacterial culture medium. The minimal gene set described herein was identified based on genes that were shown to be non-essential for bacterial growth in the medium SP4 (whose composition is described in reference #17), in the presence of tetracycline selection (the tetM tetracycline resistance gene is present in the transposon used to inactivate the genes which were shown to be non-essential). The set of non-essential genes may be different for organisms grown under different conditions (e.g. in different bacterial medium, under different selection conditions, etc). In general, a culture medium that supports growth and proliferation of a minimal organism (containing a gene set as discussed herein), with as few environmental stresses as possible, contains energy sources such as glucose, arginine or urea; protein or peptides; all amino acids; nucleotides; vitamins; cofactors; fatty acids and other membrane components such as cholesterol; enzyme cofactors; salts; minerals and buffers.
Such a medium is SP4 (Spiroplasma medium), which is a highly nutritious mixture of beef heart infusion, peptone supplemented with yeast extract, CMRL 1066 Medium and 17% fetal bovine serum. The yeast extract provides diphosphopyridine nucleotides and the serum provides cholesterol and a source of protein. (See, e.g., Tully et al. (1979) J. Infect. Dis 139, 478-82.) In particular, SP4 medium contains the following components:
Mix
Add Aseptically
The term “gene,” as used herein, refers to a polynucleotide comprising a protein-coding or RNA-coding sequence, in an expressible form, e.g. operably linked to an expression control sequence. The “coding sequences” of the gene generally do not include expression control sequences, unless they are embedded within the coding sequence. In different embodiments of the invention, the coding sequences of the genes listed in Tables 2 to 5 can be under the control of the naturally occurring expression control sequences or they can be under the control of heterologous expression control sequences, or combinations thereof.
An “expression control sequence,” as used herein, refers to a polynucleotide sequence that regulates expression of a polypeptide coded for by a polynucleotide to which it is functionally (“operably”) linked. Expression can be regulated at the level of the mRNA or polypeptide. Thus, the term expression control sequence includes mRNA-related elements and protein-related elements. Such elements include promoters, domains within promoters, ribosome binding sequences, transcriptional terminators, etc. An expression control sequence is operably linked to a nucleotide sequence when the expression control sequence is positioned in such a manner to effect or achieve expression of the coding sequence. For example, when a promoter is operably linked 5′ to a coding sequence, expression of the coding sequence is driven by the promoter.
The minimal gene set suggested in the Examples herein is composed of genes or sequences from Mycoplasma genitalium (M. genitalium) G37 (ATCC 33530). The complete genome of this bacterium is provided as Genbank accession number L43967. The individual genes are annotated in the Genbank listing as MG001, MG002 through MG470. The sequences of the genes were published on the TIGR web site in early October, 2005.
However, any of a variety of other protein- or RNA-coding genes or sequences can be substituted in a minimal gene set for the exemplified protein- or RNA-coding gene or sequences, provided that the protein or RNA encoded by the substituting gene can be expressed and that it provides a sufficient amount of the activity, function and/or structure to substitute for the M. genitalium gene or sequence in a minimal gene set. Such substitutes are sometimes referred to herein as “functional equivalents” of the exemplified genes or coding sequences.
Suitable genes or coding sequences that can be substituted include, for example, an active mutant, variant, polymorph etc. of a M. genitalium gene; or a corresponding (orthologous) gene from another bacterium, such as a different Mycoplasma species (e.g., M. capricolum). Furthermore, genes or sequences from the minimal gene set can be substituted with orthologous genes from an evolutionarily more diverse organism, such as an archaebacterium or a eukaryotic organism. Genes from eukaryotic organisms which must be post-translationally modified in order to function by a mechanism unavailable in a bacterial host cannot, of course, be used. Similarly, expression control sequences from eukaryotic genes can be used only if they can function in the background of a bacterial cell.
In one embodiment of the invention, genes from the minimal gene set are replaced by non-orthologous gene displacement (by a different set of genes providing an equivalent function or activity). For example, genes from the glycolytic pathway of M. genitalium as shown in the Examples can be substituted with genes from a different organism that utilizes a different source for generating energy (such as hydrolysis of urea, fermentation of arginine, etc.).
For example, M. genitalium generates energy via glycolysis. One can substitute a different energy generation system from another organism that would make most of the genes that express the enzymes of the glycolytic pathway superfluous. For instance energy generation in Ureaplasma parvum, a bacterium closely related to M. genitalium is based on the hydrolysis of urea. That system includes 8 genes that encode the urease enzyme complex, two ammonium transporters, and as yet unidentified nickel ion transporter (presumably one of several U. parvum cation transporters), and possibly a urea transporter (no transporter has been identified, and the very small urea molecule may enter the cell by diffusion). We expect that substitution of these 11-12 U. parvum genes for 15-20 M. genitalium genes encoding glycolytic enzymes and carbohydrate transporters would produce an organism with fewer genes capable more robust growth as is seen with U. parvum.
As used herein, the term “polynucleotide” includes a single stranded DNA corresponding to the single strand provided in the Genbank® listing, or to the complete complement thereto, or to the double stranded form of the molecule. Also included are RNA and DNA-like or RNA-like materials, such as branched DNAs, peptide nucleic acids (PNA) or locked nucleic acids (LNA).
Functional equivalents of genes can also include a variety of variant polynucleotides, provided that the variant polynucleotide can provide at least a measurable amount of the function of the original polynucleotide from which it varies. Preferably, the variant can provide at least about 50%, 75%, 90% or 95% of the function of the original polynucleotide. For example, a functional variant of a polynucleotide as described herein includes a polynucleotide that includes degenerate codons; or that is an active fragment of the original polynucleotide; or that exhibits at least about 90% identity (e.g. at least about 95% or 98% identity) with the original polynucleotide; or that can hybridize specifically to the original polynucleotide under conditions of high stringency.
Unless otherwise indicated, the term “about,” as used herein, refers to plus or minus 10%. Thus, about 90%, as used above, includes 81% to 99%. As used herein, the end points of a range are included with the range.
Functional variant polynucleotides may take a variety of forms, including, e.g., naturally or non-naturally occurring polymorphisms, including single nucleotide polymorphisms (SNPs), allelic variants, and mutants. They may comprise, e.g., one or more additions, insertions, deletions, substitutions, transitions, transversions, inversions, chromosomal translocations, variants resulting from alternative splicing events, or the like, or any combinations thereof.
The degree of sequence identity can be obtained by conventional algorithms, such as those described by Lipman and Pearson (Proc. Natl. Acad. Sci. 80:726-730, 1983) or Martinez/Needleman-Wunsch (Nucl Acid Research 11:4629-4634, 1983).
A polynucleotide that hybridizes specifically to a second polynucleotide under conditions of high stringency hybridizes preferentially to that polynucleotide. Conditions of “high stringency,” as used herein, means, for example, incubating a blot or other hybridization reaction overnight (e.g., at least 12 hours) with a long polynucleotide probe in a hybridization solution containing, e.g., about 5×SSC, 0.5% SDS, 100 μg/ml denatured salmon sperm DNA and 50% formamide, at 42° C. Blots can be washed at high stringency conditions that allow, e.g., for less than 5% by mismatch (e.g., wash twice in 0.1×SSC and 0.1% SDS for 30 min at 65° C.), thereby selecting sequences having, e.g., 95% or greater sequence identity. Other non-limiting examples of high stringency conditions include a final wash at 65° C. in aqueous buffer containing 30 mM NaCl and 0.5% SDS. Another example of high stringent conditions is hybridization in 7% SDS, 0.5 M NaPO4, pH7, 1 mM EDTA at 50° C., e.g., overnight, followed by one or more washes with a 1% SDS solution at 42° C. Whereas high stringency washes can allow for less than 5% mismatch, reduced or low stringency conditions can permit up to 20% nucleotide mismatch. Hybridization at low stringency can be accomplished as above, but using lower formamide conditions, lower temperatures and/or lower salt concentrations, as well as longer periods of incubation time.
The minimal gene set suggested herein has been derived by taking into account some of the following factors. Furthermore, the minimal gene set may be modified, e.g. for growth under other culture conditions, taking into account some of the following factors:
Although the noted protein-coding genes appear to be essential for growth under the conditions of the experiments described herein, additional protein-coding genes may be required under other conditions. For example, we isolated mutants in DNA metabolism genes that were expendable for the duration of our experiment, but might be necessary for the long-term survival of the organism. These were six genes involved in recombination and DNA repair: recA (MG339), recU (MG352), Holliday junction DNA helicases ruvA (MG358) and ruvB (MG359), formamidopyrimidine-DNA glycosylase mutM (MG262.1), which excises oxidized purines from DNA, and a likely DNA damage inducible protein gene (MG360). Perhaps because of an accumulation of cell damage over time, mutants in chromosome segregation protein SMC (MG298) and hypothetical gene MG115, which is similar to the cinA gene of Streptococcus pneumoniae competence-inducible (cin) operon, grew more poorly after repeated passage.
Even with its near minimal gene set M. genitalium has apparent enzymatic redundancy. We disrupted two complete ABC transporter gene cassettes for phosphate (MG410, MG411, MG412) and putatively phosphonate (MG289, MG290, MG291) import. The PhoU regulatory protein gene (MG409) was not disrupted, suggesting it is needed for both cassettes. Phosphate is an essential metabolite that must be imported. Either phosphate might be imported by both transporters as a result of relaxed substrate specificity by the phosphonate system, or there is a metabolic capacity to interconvert phosphate and phosphonate. Although we disrupted both of these three gene cassettes, cells presumably need at least one phosphate transporter. Therefore, a minimal gene set preferably contains three ABC transporter genes for phosphate importation. Relaxed substrate specificity is a recurring theme proposed and shown for several M. genitalium enzymes as a mechanism by which this bacterium meets its metabolic needs with fewer genes (21, 22).
M. genitalium generates ATP through glycolysis, and although none of the genes encoding enzymes involved in the initial glycolytic reactions were disrupted, mutations in two energy generation genes suggested there may be still more unexpected genomic redundancy in this essential pathway. We identified viable insertion mutants in genes encoding lactate/malate dehydrogenase (MG460) and the dihydrolipoamide dehydrogenase subunit of the pyruvate dehydrogenase complex (MG271). Mutations in either of these dehydrogenases would be expected to have glycolytic ATP production, and unbalanced NAD+ and NADH levels, which are the primary oxidizing and reducing agents in glycolysis. These mutations should have greatly reduced growth rate and accelerated acidification of the growth medium While the MG271 mutants grew about 20% slower than wild type cells, inexplicably, the lactate dehydrogenase mutants grow ˜20% faster. We also isolated a mutant in glycerol-3-phospate dehydrogenase (MG039), a phospholipid biosynthesis enzyme. The loss of functions in these mutants could have been compensated for by other M. genitalium dehydrogenases or reductases. This could be another case of mycoplasma enzymes having a relaxed substrate specificity as has been reported for lactate/malate dehydrogenase (21) and nucleotide kinases (22).
Under our laboratory conditions we identified 100 non-essential protein-coding genes. It appears that the remaining 382 M. genitalium protein-coding genes, plus three phosphate transporter genes, and 43 RNA-coding genes comprise the essential genes set for this minimal cell (Table 3). We disrupted genes in only 5 of the 12 M. genitalium paralogous gene families. Only for the two families comprised of lipoproteins MG185 and MG260 and glycerophosphoryl diester phosphodiesterases MG293 and MG385 did we disrupt all members. Accordingly, these families' functions may be essential, and we expanded our projection of the essential gene set to 387 genes to include them (one each of MG185 or MG260, and one each of MG293 and MG385). This is a significantly greater number of essential genes than the 265-350 predicted in the inventors' previous study of M. genitalium (4), or in the gene knockout/disruption study that identified 279 essential genes in B. subtilis, which is a more conventional bacterium from the same Firmicutes taxon as M. genitalium (6). Similarly, our finding of 387 essential protein-coding genes greatly exceeds theoretical projections of how many genes comprise a minimal genome such as Mushegian and Koonin's 256 genes shared by both H. influenzae and M. genitalium (2), and the 206 gene core of a minimal bacterial gene set proposed by Gil et al (3). One of the surprises about the present essential gene set is its inclusion of 108 hypothetical proteins and proteins of unknown function.
These data suggest that a genome constructed to encode the 387 protein-coding and 43 structural RNA genes could sustain a viable synthetic cell, which has been referred to hypothetically as a Mycoplasma laboratorium (24). A variety of mechanisms can be used for preparing such a viable synthetic cell. For example, the minimal gene set can be introduced into a ghost cell, from which the resident genome has been removed or disabled. In one embodiment, ribosomes, membranes and other cellular components important for gene regulation, transcription, translation, post-transcriptional modification, secretion, uptake of nutrients or other substances, etc., are present in the ghost cell. In another embodiment, one or more of these components is prepared synthetically.
In one embodiment of the invention, the genes in the minimal gene set, or a subset of those genes, are cloned into conventional vectors, e.g. to form a library. The DNA to be cloned can be obtained from any suitable source, including naturally occurring genes, genes previously cloned into a different vector, or artificially synthesized genes. The genes may be cloned by in vitro, synthetic procedures, such as those disclosed in co-pending PCT application PCT/US06/16349, filed 1 May 2006, “Amplification and Cloning of Single DNA Molecules Using Rolling Circle Amplification,” incorporated by reference herein in its entirety. For example, synthetically prepared genes of the gene set may be amplified and assembled to font′ a synthetic gene or genome. This can be performed by diluting DNA molecules, such that each sample of diluted DNA contains, on average, one molecule of DNA, in fragments of about 5 kb, for example, and then converting to single stranded DNA circles, and then amplifying the DNA circles using Φ29 polymerase.
As a library, the gene sets of the invention can be arranged in any font′, in single or multiple copies, and can be arranged in individual oligonucleotides each having a section of one of the genes, one of the genes, or more than one of the genes. These oligonucleotides can be arranged as cassettes. The cassettes can be joined up to form larger gene assemblies, including a minimal genome comprising or consisting of all the genes of the gene set of the invention. The genes can be assembled by a method such as that described in PCT International Patent Application No. PCT/US06/31214, filed 11 Aug. 2006, “Method For In Vitro Recombination Employing a 3′ Exonuclease Activity,” incorporated by reference herein in its entirety. PCT/US06/31214 describes methods of joining cassettes of genes into larger assemblies, and can be used to produce a single DNA molecule comprising the gene set of the invention. In particular, that application describes an in vitro method, using isolated proteins, for joining two or more double-stranded (ds) DNA molecules of interest, wherein the distal region of the first DNA molecule and the proximal region of the second DNA molecule of each pair share a region of sequence identity, comprising (a) treating the DNA molecules with an enzyme having an exonuclease activity, under conditions effective to yield single-stranded overhanging portions of each DNA molecule which contain a sufficient length of the region of sequence homology to hybridize specifically to the region of sequence homology of its pair; (b) incubating the treated DNA molecules of (a) under conditions effective to achieve specific annealing of the single-stranded overhanging portions; and (c) treating the incubated DNA molecules in (b) under conditions effective to fill in remaining single-stranded gaps and to seal the nicks thus formed, wherein the region of sequence identity comprises at least 20 non-palindromic nucleotides (nt).
The DNA molecules of the library may have a size of any practical length. The lower size limit for a dsDNA to circularize is about 200 base pairs. Therefore, the total length of the joined fragments (including, in some cases, the length of the vector) is preferably at least about 200 bp in length. The DNAs can take the form of either a circle or a linear molecule. The library may include from two to a very large number of DNA molecules, which can be joined together. In general, at least about 10 fragments can be joined.
More particularly, the number of DNA molecules or cassettes that may be joined to produce an end product, in one or several assembly stages, may be at least or no greater than about 2, 3, 4, 6, 8, 10, 15, 20, 25, 50, 100, 200, 500, 1000, 5000, or 10,000 DNA molecules, for example in the range of about 4 to about 100 molecules. The DNA molecules or cassettes in a library of the invention may each have a starting size in a range of at least or no greater than about 80 bs, 100 bs, 500 bs, 1 kb, 3 kb, 5 kb, 6 kb, 10 kb, 18 kb, 20 kb, 25 kb, 32 kb, 50 kb, 65 kb, 75 kb, 150 kb, 300 kb, 500 kb, 600 kb, or larger, for example in the range of about 3 kb to about 100 kb. According to the invention, methods may be used for assembly of about 100 cassettes of about 6 kb each, into a DNA molecule of about 600 kb.
One embodiment of the invention is to join cassettes, such as 5-6 kb DNA molecules representing adjacent regions of a gene or genome included in a gene set of the invention, to create combinatorial assemblies. For example, it may be of interest to modify a bacterial genome, such as a putative minimal genome or a minimal genome, so that one or more of the genes is eliminated or mutated, and/or one or more additional genes is added. Such modifications can be carried out by dividing the genome into suitable cassettes, e.g. of about 5-6 kb, and assembling a modified genome by substituting a cassette containing the desired modification for the original cassette. Furthermore, if it is desirable to introduce a variety of changes simultaneously (e.g. a variety of modifications of a gene of interest, the addition of a variety of alternative genes, the elimination of one or more genes, etc.), one can assemble a large number of genomes simultaneously, using a variety of cassettes corresponding to the various modifications, in combinatorial assemblies. After the large number of modified sequences is assembled, preferably in a high throughput manner, the properties of each of the modified genomes can be tested to determine which modifications confer desirable properties on the genome (or an organism comprising the genome). This “mix and match” procedure produces a variety of test genomes or organisms whose properties can be compared. The entire procedure can be repeated as desired in a recursive fashion.
Methods of cloning, as well as many of the other molecular biological methods used in conjunction with the present invention, are discussed, e.g., in Sambrook, et al. (1989), Molecular Cloning, a Laboratory Manual, Cold Harbor Laboratory Press, Cold Spring Harbor, N.Y.; Ausubel et al. (1995). Current Protocols in Molecular Biology, N.Y., John Wiley & Sons; Davis et al. (1986), Basic Methods in Molecular Biology, Elsevier Sciences Publishing, Inc., New York; Hames et al. (1985), Nucleic Acid Hybridization, IL Press; Dracopoli et al. Current Protocols in Human Genetics, John Wiley & Sons, Inc.; and Coligan et al. Current Protocols in Protein Science, John Wiley & Sons, Inc.
Another aspect of the invention is a set of genes or polynucleotides on the invention which are in a free-living organism. The organism may be in a dormant or resting state (e.g., lyophilized, stored in a suitable solution, such as glycerol, or stored in culture medium), or it may growing and/or replicating, for example in a rich culture medium, such as SP4.
Another aspect of the invention is a set of polypeptides encoded by a set of genes or polynucleotides of the invention. The polypeptides may be, e.g., in a free-living organism.
Another aspect of the invention is a set of genes or polynucleotides of the invention that are recorded on computer readable media. As used herein, “computer readable media” refers to any medium that can be read and accessed directly by a computer. Such media include, but are not limited to: magnetic storage media, such as floppy discs, hard disc storage medium, and magnetic tape; optical storage media such as CD-ROM; electrical storage media such as RAM and ROM; and hybrids of these categories such as magnetic/optical storage media. The skilled artisan will readily appreciate how any of the presently known computer readable media can be used to create a manufacture comprising computer readable medium having recorded thereon a polynucleotide or amino acid sequence of the present invention.
As used herein, “recorded” refers to a process for storing information on computer readable medium. The skilled artisan can readily adopt any of the presently known methods for recording information on computer readable medium to generate manufactures comprising the nucleotide or amino acid sequence information of the present invention.
A variety of data storage structures are available to a skilled artisan for creating a computer readable medium having recorded thereon a set of nucleotide or amino acid sequences of the present invention. The choice of the data storage structure will generally be based on the means chosen to access the stored information. In addition, a variety of data processor programs and formats can be used to store the nucleotide sequence information of the present invention on computer readable medium. The sequence information can be represented in a word processing text file, formatted in commercially-available software such as WordPerfect® and Microsoft® Word, or represented in the form of an ASCII file, stored in a database application, such as DB2®, Sybase®, Oracle®, or the like. The skilled artisan can readily adapt any number of data processor structuring formats (e.g., text file or database) in order to obtain computer readable medium having recorded thereon the nucleotide sequence information of the present invention.
By providing a set of nucleotide or amino acid sequences of the invention in computer readable form, the skilled artisan can routinely access the sequence information for a variety of purposes. For example, one skilled in the art can use the nucleotide or amino acid sequences of the invention in computer readable form to compare the sequences with orthologous sequences that can be substituted for the present sequences in an alternative version of the minimal genome. Computer software is publicly available which allows a skilled artisan to access sequence information provided in a computer readable medium for analysis and comparison to other sequences. A variety of known algorithms are disclosed publicly and a variety of commercially available software for conducting search means are and can be used in the computer-based systems of the present invention. Examples of such software include, but are not limited to, MacPattern (EMBL), BLASTN and BLASTX (NCBIA).
For example, software which implements the BLAST (Altschul et al. (1990) J. Mol. Biol. 215:403-410) and BLAZE (Brutlag et al. (1993) Comp. Chem. 17:203-207) search algorithms on a Sybase system can be used to identify open reading frames (ORFs) of the sequences of the invention which contain homology to ORFs or proteins from other libraries. Such ORFs are protein encoding fragments and are useful in producing commercially important proteins such as enzymes used in various reactions and in the production of commercially useful metabolites.
In the foregoing and in the following example, all temperatures are set forth in uncorrected degrees Celsius; and, unless otherwise indicated, all parts and percentages are by weight.
A. Cells and plasmids. We obtained wild type M. genitalium 037 (ATCC® Number: 33530™) from the American Type Culture Collection (Manassas, Va.). As part of this project we re-sequenced and re-annotated the genome of this bacterium. The new M. genitalium 037 sequence (Genbank accession number L43967) differed from the previous M. genitalium (13) genome sequence at 34 sites. Several genes previously listed as having frameshifts were merged including MG016, MG017, and MG018 (DEAD helicase) and MG419 and MG420 (DNA polymerase III gamma/tau subunit). Our transposon mutagenesis vector was the plasmid pIVT-1, which contains the Tn4001 transposon with a tetracycline resistance gene (tetM)(15), and was a gift from Dr. Kevin Dybvig at the University of Alabama at Birmingham.
B. Transformation of M. genitalium with Tn4001 by electroporation. Confluent flasks of M. genitalium cells were harvested by scraping into electroporation buffer (EB) comprised of 8 mM HEPES+272 mM sucrose at pH 7.4. We washed and then resuspended the cells in a total volume of 200-300 μl EB. On ice, 100 μl cells were mixed with 30 μg pIVT-1 plasmid DNA and transferred to a 2 mm chilled electroporation cuvette (BioRad®, Hercules, Calif.). We electroporated using 2500 V, 25 μF, and 100Ω. After electroporation we resuspended the cells in 1 ml of 37° C. SP4 medium and allowed the cells to recover for 2 hours at 37° C. with 5% CO2. Aliquots of 200 μl of cells were spread onto SP4 agar plates containing 2 mg/l tetracycline hydrochloride (VWR®, Bridgeport, N.J.). The plates were incubated for 3-4 weeks at 37° C. with 5% CO2 until colonies were visible. When colonies were 3-4 weeks old, we transferred individual M. genitalium colonies into SP4 medium +7 mg/L tetracycline in 96 well plates. We incubated the plates at 37° C. with 5% CO2 until the SP4 in most of the wells began to turn acidic and became yellow or orange (˜4 days). We froze those mutant stock cells at −80° C.
C. Amplification of isolated colonies for DNA extraction. We inoculated 4 ml SP4 containing 7 μg/ml tetracycline in 6 well plates with 20 μl transposon mutant stock cells and incubated the plates at 37° C. with 5% CO2 until the cells reached 100% confluence. To extract genomic DNA from confluent cells, we scraped the cells and then transferred the cell suspension to a tube for pelleting by centrifugation. Thus any non-adherent cells were not lost. We washed the cells in PBS (MediaTech®, Herndon, Va.) and then resuspended them in a mixture of 100 μl PBS and 100 μl of the chaotropic MTL buffer from a Qiagen® MagAttract® DNA Mini M48 Kit (QIAGEN, Valencia, Calif.). Tubes were stored at −20° C. until the genomic DNA could be extracted using a Qiagen® BioRobot® M48 workstation (Qiagen®).
D. Location of Tn4001 tet insertion sites by DNA sequencing from M. genitalium genomic templates. Our 20 μl sequencing reactions contained ˜0.5 μg of genomic DNA, 6.4 pmol of the 30 base oligonucleotide GTACTCAATGAATTAGGTGGAAGACCGAGG (SEQ ID NO:1) (Integrated DNA Technologies®, Coralville, Iowa). The primer binds in the tetM gene 103 basepairs from one of the transposon/genome junctions. Using BLAST we located the insertion site on the M. genitalium genome.
E. Quantitative PCR to determine colony homogeneity and genes duplication. We designed quantitative PCR primers (Integrated DNA Technologies®) flanking transposon insertion sites using the default conditions for the primer design software Primer Express 1.5 (Applied Biosystems®). Using quantitative PCR done on an Applied Biosystems® 7700 Sequence Detection System, we determined the amounts of the target genes lacking a Tn4001 insertion in genomic DNA prepared from mutant colonies relative to a the amount of the those genes in wild type M. genitalium. Reactions were done in Eurogentec qPCR Mastermix® Plus SYBR Green (San Diego, Calif.). Genomic DNA concentrations were normalized after determining their relative amounts using a TagMan® quantitative PCR specific for the 16S rRNA gene that was done in Eurogentec qPCR Mastermix® Plus. We calculated the amounts of target genes lacking the transposon in mutant genomic DNA preparations relative to the amounts in wild type using the delta-delta Ct method (16).
We sequenced across the transposon-genome junctions of our mutants using a primer specific for Tn4001tet. Presence of a transposon in the central region of a gene of a viable bacterium indicated that gene was disrupted and therefore non-essential (dispensable). We considered transposon insertions disruptive only if they were after the first three codons and before the 3′-most 20% of the coding sequence of a gene. Thus, non-disruptive mutations resulting from transposon mediated duplication of short sequences at the insertion site (18, 19), and potentially inconsequential COOH-terminal insertions do not result in erroneous determination of gene expendability. Without wishing to be bound by any particular theory, it is suggested that these disruptions actually occurred, even though theoretically, some genes might tolerate transposon insertions, and we did not confirm the absence of the gene products. To exclude the possibility that gene disruptions were the result of a transposon insertion in one copy of a duplicated gene, we used PCR to detect genes lacking the insertion. This showed us that almost all of our colonies contained both disrupted and wild type versions of the genes identified as having the Tn4001. Further analysis using quantitative PCR showed most colonies were mixtures of two or more mutants, thus we operationally refer to them and any DNA isolated from them as colonies rather than clones. This cell clumping led us to isolate individual mutants using filter cloning. To do this we forced cells through 0.22 μm filters before plating to break up clumps of cells possibly containing multiple different mutants. We used these cells to produce subcolonies which we both sequenced and analyzed using quantitative PCR. For each disrupted gene we subcloned at least one primary colony.
In total we analyzed 3,152 M. genitalium transposon insertion mutant primary colonies, and subcolonies to determine the locations of Tn4001tet inserts. For 75% of these we generated sequence data that enabled us to map the transposon insertion sites. Colonies containing multiple Tn4001tet insertions cannot be characterized using this approach. Only 62% of primary colonies generated useful sequence. This was likely because of the tendency of mycoplasma cells to form persistent cell aggregates leading to colonies containing mixtures of multiple mutants that proved refractory to sequencing. For subcolonies the success rate was 82%. Of the successfully sequenced subcolonies in 59% the transposon insert was at a different site than in the parental primary colony. The rate at which we identified mutants with previously unhit insertion sites on the genome was higher for the primary colonies than the subcolonies. However the rate of accumulation of new insertion sites dropped after our first 600 colonies, indicating we were approaching saturation mutagenesis of all non-lethal insertion sites (
We mapped a total of 2293 different transposon insertion sites on the genome (
We wanted to determine if any of our disrupted genes were in cells bearing two copies of the gene. Unexpectedly, PCRs using primers flanking the transposon insertion sites produced amplicons of the size expected for wild type templates from all 5 colonies initially tested. End-stage analysis of PCRs could not tell us if the wild type sequences we amplified were the result of a low level of transposon jumping out of the target gene, or if there was a gene duplication. To address this, for at least one colony or subcolony for each disrupted gene we used quantitative PCR to measure how many copies of contaminating wild type versions of that gene there were in the sequenced DNA preps.
Analysis of the quantitative PCR results showed most colonies were mixtures of multiple mutants. This was likely a consequence of our high transformation efficiency and the tendency of mycoplasma cells to aggregate. The direct genomic sequencing identified only the plurality member of the population. To address this issue we adapted our mutant isolation protocol to include one or two rounds of filter cloning. Existing colonies of interest were filter subcloned. We isolated 10 subcolonies and the sites of their Tn4001 insertions were determined. We took both rapidly growing colonies and M. genitalium colonies that were delayed in their appearance. Often only a minority of the subcolonies had inserts in the same location as found with the parental colony. After filter cloning we still found that almost every subcolony had some low level of a wild type copy of the disrupted gene. This is likely the result of Tn4001 jumping(20). After subcloning we were able to isolate gene disruption mutant colonies for 99 of our 100 different disrupted M. genitalium genes that had less than 1% wild type sequence.
Several mutants manifested remarkable phenotypes. While many of the mutants grew slowly, mutants in lactate/malate dehydrogenase (MG460), and conserved hypothetical proteins MG414 and MG415 mutants had doubling times up to 20% faster than wild type M. genitalium (data not shown). Cells with transposon insertions in the transketolase gene (MG066), which encodes a membrane protein and pentose phosphate pathway enzyme, grew in chains of clumped cells rather than in the monolayers characteristic of wild type M. genitalium. Other mutant cells grew in suspension rather than adhering to plastic. Some cells would lyse when washed with PBS, and thus had to be processed in either SP4 medium or 100% serum.
We isolated mutants with transposon insertions at some sites much more frequently than others (
As noted above, at least 387 protein-coding genes and all of the RNA genes are essential and could form a minimal set. However, it seems unlikely that all of those “one-at-a time” dispensable genes could be eliminated simultaneously. To determine a subset that can be simultaneously deleted, a wild type chromosome is constructed synthetically. The synthetic genome is constructed hierarchically from chemically synthesized oligonucleotides. Subsets of the dispensable genes are then removed. The synthetic natural chromosome and the reduced genome are tested for viability by transplantation into cells from which the resident chromosome has been removed. Rapid advances in gene synthesis technology and efforts at developing genome transplantation methods allow the confirmation that the M. genitalium essential gene set described above is a true minimal gene set, or provide a basis to modify that gene set.
Mycoplasma genitalium
genitalium
Bacillus subtilis
Escherichia coli (K-12)
Haemophilus influenzae
Helicobacter pylori
Mycobacterium bovis
Pseudomonas aeruginosa
Staphylococcus aureus
We used a common definition for members of paralogous gene families requiring they have 30% identity over 60% of the length of the longer protein sequence (a single linkage clustering then defines the families).
Mycoplasma genitalium genes with Tn4001tet insertions that are disrupted. Genes
All information is based on the M. genitalium genome sequence and annotation reported herein. Genes are grouped by main biological roles. The columns are as follows:
M. genitalium gene locus
Gene symbol
Gene common name
Mycoplasma genitalium protein coding genes that were not disrupted in this study.
All information is based on the M. genitalium genome sequence and annotation reported herein. Genes are grouped by main biological roles. The columns for the protein coding genes are as follows:
M. genitalium gene locus
Gene symbol
Gene common name
Mycoplasma genitalium genes with Tn4001tet insertions that were not reported as
All information is based on the new M. genitalium genome sequence and annotation reported here. Genes are grouped by main biological roles. The columns are as follows:
M. genitalium gene locus
Gene symbol
Gene common name
Mycoplasma genitalium genes with Tn4001tet insertions that were not reported
All information is based on the M. genitalium genome sequence and annotation reported herein. Genes are grouped by main biological roles. The columns for these protein coding genes are as follows:
M. genitalium gene locus
Gene symbol
Gene common name
From the foregoing description, one skilled in the art can easily ascertain the essential characteristics of this invention, and without departing from the spirit and scope thereof, can make changes and modifications of the invention to adapt it to various usage and conditions and to utilize the present invention to its fullest extent. The preceding specific embodiments are to be construed as merely illustrative, and not limiting of the scope of the invention in any way whatsoever. The entire disclosure of all applications, patents, publications (including U.S. provisional application 60/725,295, filed Oct. 12, 2005) cited above and in the figures, are hereby incorporated in their entirety by reference.
This application is a divisional of Ser. No. 11/546,364, filed Oct. 12, 2006, which claims the benefit of U.S. provisional application 60/725,295, filed Oct. 12, 2005, each of which is incorporated by reference herein in their entireties, including all Tables, Figures, and Claims.
Aspects of this invention were made with government support (DOE grant number DE-FG02-02ER63453). The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
60725295 | Oct 2005 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11546364 | Oct 2006 | US |
Child | 14733743 | US |