Toward the Complete Functional Characterization of a Minimal Bacterial Proteome.

Recently, we presented a whole-cell kinetic model of the genetically minimal bacterium JCVI-syn3A that described the coupled metabolic and genetic information processes and predicted behaviors emerging from the interactions among these networks. JCVI-syn3A is a genetically reduced bacterial cell that has the fewest number and smallest fraction of genes of unclear function, with approximately 90 of its 452 protein-coding genes (that is less than 20%) unannotated. Further characterization of unclear JCVI-syn3A genes strengthens the robustness and predictive power of cell modeling efforts and can lead to a deeper understanding of biophysical processes and pathways at the cell scale. Here, we apply computational analyses to elucidate the functions of the products of several essential but previously uncharacterized genes involved in integral cellular processes, particularly those directly affecting cell growth, division, and morphology. We also suggest directed wet-lab experiments informed by our analyses to further understand these "missing puzzle pieces" that are an essential part of the mosaic of biological interactions present in JCVI-syn3A. Our workflow leverages evolutionary sequence analysis, protein structure prediction, interactomics, and genome architecture to determine upgraded annotations. Additionally, we apply the structure prediction analysis component of our work to all 452 protein coding genes in JCVI-syn3A to expedite future functional annotation studies as well as the inverse mapping of the cell state to more physical models requiring all-atom or coarse-grained representations for all JCVI-syn3A proteins.

Synthetic chromosomes, genomes, viruses, and cells.

Synthetic genomics is the construction of viruses, bacteria, and eukaryotic cells with synthetic genomes. It involves two basic processes: synthesis of complete genomes or chromosomes and booting up of those synthetic nucleic acids to make viruses or living cells. The first synthetic genomics efforts resulted in the construction of viruses. This led to a revolution in viral reverse genetics and improvements in vaccine design and manufacture. The first bacterium with a synthetic genome led to construction of a minimal bacterial cell and recoded Escherichia coli strains able to incorporate multiple non-standard amino acids in proteins and resistant to phage infection. Further advances led to a yeast strain with a synthetic genome and new approaches for animal and plant artificial chromosomes. On the horizon there are dramatic advances in DNA synthesis that will enable extraordinary new opportunities in medicine, industry, agriculture, and research.

Fundamental behaviors emerge from simulations of a living minimal cell.

We present a whole-cell fully dynamical kinetic model (WCM) of JCVI-syn3A, a minimal cell with a reduced genome of 493 genes that has retained few regulatory proteins or small RNAs. Cryo-electron tomograms provide the cell geometry and ribosome distributions. Time-dependent behaviors of concentrations and reaction fluxes from stochastic-deterministic simulations over a cell cycle reveal how the cell balances demands of its metabolism, genetic information processes, and growth, and offer insight into the principles of life for this minimal cell. The energy economy of each process including active transport of amino acids, nucleosides, and ions is analyzed. WCM reveals how emergent imbalances lead to slowdowns in the rates of transcription and translation. Integration of experimental data is critical in building a kinetic model from which emerges a genome-wide distribution of mRNA half-lives, multiple DNA replication events that can be compared to qPCR results, and the experimentally observed doubling behavior.

Kinetic Modeling of the Genetic Information Processes in a Minimal Cell.

JCVI-syn3A is a minimal bacterial cell with a 543 kbp genome consisting of 493 genes. For this slow growing minimal cell with a 105 min doubling time, we recently established the essential metabolism including the transport of required nutrients from the environment, the gene map, and genome-wide proteomics. Of the 452 protein-coding genes, 143 are assigned to metabolism and 212 are assigned to genetic information processing. Using genome-wide proteomics and experimentally measured kinetic parameters from the literature we present here kinetic models for the genetic information processes of DNA replication, replication initiation, transcription, and translation which are solved stochastically and averaged over 1,000 replicates/cells. The model predicts the time required for replication initiation and DNA replication to be 8 and 50 min on average respectively and the number of proteins and ribosomal components to be approximately doubled in a cell cycle. The model of genetic information processing when combined with the essential metabolic and cell growth networks will provide a powerful platform for studying the fundamental principles of life.

Polar Effects of Transposon Insertion into a Minimal Bacterial Genome.

Global transposon mutagenesis is a valuable tool for identifying genes required for cell viability. Here we present a global analysis of the orientation of viable Tn5-Puror (Tn5-puromycin resistance) insertions into the near-minimal bacterial genome of JCVI-syn2.0. Sixteen of the 478 protein-coding genes show a noticeable asymmetry in the orientation of disrupting insertions of Tn5-Puror Ten of these are located in operons, upstream of essential or quasi-essential genes. Inserts transcribed in the same direction as the downstream gene are favored, permitting read-through transcription of the essential or quasi-essential gene. Some of these genes were classified as quasi-essential solely because of polar effects on the expression of downstream genes. Three genes showing asymmetry in Tn5-Puror insertion orientation prefer the orientation that avoids collisions between read-through transcription of Tn5-Puror and transcription of an adjacent gene. One gene (JCVISYN2_0132 [abbreviated here as "_0132"]) shows a strong preference for Tn5-Puror insertions transcribed upstream, away from the downstream nonessential gene _0133. This suggested that expression of _0133 due to read-through from Tn5-Puror is lethal when _0132 function is disrupted by transposon insertion. This led to the identification of genes _0133 and _0132 as a toxin-antitoxin pair. The three remaining genes show read-through transcription of Tn5-Puror directed downstream and away from sizable upstream intergenic regions (199 bp to 363 bp), for unknown reasons. In summary, polar effects of transposon insertion can, in a few cases, affect the classification of genes as essential, quasi-essential, or nonessential and sometimes can give clues to gene function.IMPORTANCE In studies of the minimal genetic requirements for life, we used global transposon mutagenesis to identify genes needed for a minimal bacterial genome. Transposon insertion can disrupt the function of a gene but can also have polar effects on the expression of adjacent genes. In the Tn5-Puror construct used in our studies, read-through transcription from Tn5-Puror can drive expression of downstream genes. This results in a preference for Tn5-Puror insertions transcribed toward a downstream essential or quasi-essential gene within the same operon. Such polar effects can have an impact on the classification of genes as essential, quasi-essential, or nonessential, but this has been observed in only a few cases. Also, polar effects of Tn5-Puror insertion can sometimes give clues to gene function.

Essential metabolism for a minimal cell.

JCVI-syn3A, a robust minimal cell with a 543 kbp genome and 493 genes, provides a versatile platform to study the basics of life. Using the vast amount of experimental information available on its precursor, Mycoplasma mycoides capri, we assembled a near-complete metabolic network with 98% of enzymatic reactions supported by annotation or experiment. The model agrees well with genome-scale in vivo transposon mutagenesis experiments, showing a Matthews correlation coefficient of 0.59. The genes in the reconstruction have a high in vivo essentiality or quasi-essentiality of 92% (68% essential), compared to 79% in silico essentiality. This coherent model of the minimal metabolism in JCVI-syn3A at the same time also points toward specific open questions regarding the minimal genome of JCVI-syn3A, which still contains many genes of generic or completely unclear function. In particular, the model, its comparison to in vivo essentiality and proteomics data yield specific hypotheses on gene functions and metabolic capabilities; and provide suggestions for several further gene removals. In this way, the model and its accompanying data guide future investigations of the minimal cell. Finally, the identification of 30 essential genes with unclear function will motivate the search for new biological mechanisms beyond metabolism.

One way that researchers can test whether they understand a biological system is to see if they can accurately recreate it as a computer model. The more they learn about living things, the more the researchers can improve their models and the closer the models become to simulating the original. In this approach, it is best to start by trying to model a simple system. Biologists have previously succeeded in creating 'minimal bacterial cells'. These synthetic cells contain fewer genes than almost all other living things and they are believed to be among the simplest possible forms of life that can grow on their own. The minimal cells can produce all the chemicals that they need to survive ­ in other words, they have a metabolism. Accurately recreating one of these cells in a computer is a key first step towards simulating a complete living system. Breuer et al. have developed a computer model to simulate the network of the biochemical reactions going on inside a minimal cell with just 493 genes. By altering the parameters of their model and comparing the results to experimental data, Breuer et al. explored the accuracy of their model. Overall, the model reproduces experimental results, but it is not yet perfect. The differences between the model and the experiments suggest new questions and tests that could advance our understanding of biology. In particular, Breuer et al. identified 30 genes that are essential for life in these cells but that currently have no known purpose. Continuing to develop and expand models like these to reproduce more complex living systems provides a tool to test current knowledge of biology. These models may become so advanced that they could predict how living things will respond to changing situations. This would allow scientists to test ideas sooner and make much faster progress in understanding life on Earth. Ultimately, these models could one day help to accelerate medical and industrial processes to save lives and enhance productivity.

Tuning Gene Activity by Inducible and Targeted Regulation of Gene Expression in Minimal Bacterial Cells.

Functional genomics studies in minimal mycoplasma cells enable unobstructed access to some of the most fundamental processes in biology. Conventional transposon bombardment and gene knockout approaches often fail to reveal functions of genes that are essential for viability, where lethality precludes phenotypic characterization. Conditional inactivation of genes is effective for characterizing functions central to cell growth and division, but tools are limited for this purpose in mycoplasmas. Here we demonstrate systems for inducible repression of gene expression based on clustered regularly interspaced short palindromic repeats-mediated interference (CRISPRi) in Mycoplasma pneumoniae and synthetic Mycoplasma mycoides, two organisms with reduced genomes actively used in systems biology studies. In the synthetic cell, we also demonstrate inducible gene expression for the first time. Time-course data suggest rapid kinetics and reversible engagement of CRISPRi. Targeting of six selected endogenous genes with this system results in lowered transcript levels or reduced growth rates that agree with lack or shortage of data in previous transposon bombardment studies, and now produces actual cells to analyze. The ksgA gene encodes a methylase that modifies 16S rRNA, rendering it vulnerable to inhibition by the antibiotic kasugamycin. Targeting the ksgA gene with CRISPRi removes the lethal effect of kasugamycin and enables cell growth, thereby establishing specific and effective gene modulation with our system. The facile methods for conditional gene activation and inactivation in mycoplasmas open the door to systematic dissection of genetic programs at the core of cellular life.

Minimal Cells-Real and Imagined.

A minimal cell is one whose genome only encodes the minimal set of genes necessary for the cell to survive. Scientific reductionism postulates the best way to learn the first principles of cellular biology would be to use a minimal cell in which the functions of all genes and components are understood. The genes in a minimal cell are, by definition, essential. In 2016, synthesis of a genome comprised of only the set of essential and quasi-essential genes encoded by the bacterium Mycoplasma mycoides created a near-minimal bacterial cell. This organism performs the cellular functions common to all organisms. It replicates DNA, transcribes RNA, translates proteins, undergoes cell division, and little else. In this review, we examine this organism and contrast it with other bacteria that have been used as surrogates for a minimal cell.

Design and synthesis of a minimal bacterial genome.

We used whole-genome design and complete chemical synthesis to minimize the 1079-kilobase pair synthetic genome of Mycoplasma mycoides JCVI-syn1.0. An initial design, based on collective knowledge of molecular biology combined with limited transposon mutagenesis data, failed to produce a viable cell. Improved transposon mutagenesis methods revealed a class of quasi-essential genes that are needed for robust growth, explaining the failure of our initial design. Three cycles of design, synthesis, and testing, with retention of quasi-essential genes, produced JCVI-syn3.0 (531 kilobase pairs, 473 genes), which has a genome smaller than that of any autonomously replicating cell found in nature. JCVI-syn3.0 retains almost all genes involved in the synthesis and processing of macromolecules. Unexpectedly, it also contains 149 genes with unknown biological functions. JCVI-syn3.0 is a versatile platform for investigating the core functions of life and for exploring whole-genome design.

Designer diatom episomes delivered by bacterial conjugation.

Eukaryotic microalgae hold great promise for the bioproduction of fuels and higher value chemicals. However, compared with model genetic organisms such as Escherichia coli and Saccharomyces cerevisiae, characterization of the complex biology and biochemistry of algae and strain improvement has been hampered by the inefficient genetic tools. To date, many algal species are transformable only via particle bombardment, and the introduced DNA is integrated randomly into the nuclear genome. Here we describe the first nuclear episomal vector for diatoms and a plasmid delivery method via conjugation from Escherichia coli to the diatoms Phaeodactylum tricornutum and Thalassiosira pseudonana. We identify a yeast-derived sequence that enables stable episome replication in these diatoms even in the absence of antibiotic selection and show that episomes are maintained as closed circles at copy number equivalent to native chromosomes. This highly efficient genetic system facilitates high-throughput functional characterization of algal genes and accelerates molecular phytoplankton research.

Bacterial genome reduction using the progressive clustering of deletions via yeast sexual cycling.

The availability of genetically tractable organisms with simple genomes is critical for the rapid, systems-level understanding of basic biological processes. Mycoplasma bacteria, with the smallest known genomes among free-living cellular organisms, are ideal models for this purpose, but the natural versions of these cells have genome complexities still too great to offer a comprehensive view of a fundamental life form. Here we describe an efficient method for reducing genomes from these organisms by identifying individually deletable regions using transposon mutagenesis and progressively clustering deleted genomic segments using meiotic recombination between the bacterial genomes harbored in yeast. Mycoplasmal genomes subjected to this process and transplanted into recipient cells yielded two mycoplasma strains. The first simultaneously lacked eight singly deletable regions of the genome, representing a total of 91 genes and ∼ 10% of the original genome. The second strain lacked seven of the eight regions, representing 84 genes. Growth assay data revealed an absence of genetic interactions among the 91 genes under tested conditions. Despite predicted effects of the deletions on sugar metabolism and the proteome, growth rates were unaffected by the gene deletions in the seven-deletion strain. These results support the feasibility of using single-gene disruption data to design and construct viable genomes lacking multiple genes, paving the way toward genome minimization. The progressive clustering method is expected to be effective for the reorganization of any mega-sized DNA molecules cloned in yeast, facilitating the construction of designer genomes in microbes as well as genomic fragments for genetic engineering of higher eukaryotes.

Rescue of mutant fitness defects using in vitro reconstituted designer transposons in Mycoplasma mycoides.

With only hundreds of genes contained within their genomes, mycoplasmas have become model organisms for precise understanding of cellular processes, as well as platform organisms for predictable engineering of microbial functions for mission-critical applications. Despite the availability of "whole genome writing" in Mycoplasma mycoides, some traditional methods for genetic engineering are underdeveloped in mycoplasmas. Here we demonstrate two facile transposon-mediated approaches for introducing genes into the synthetic cell based on M. mycoides. The marker-less approach involves preparing a fragment containing only a small genomic region of interest with flanking transposase-binding sites, followed by in vitro transposase loading and introduction into the cells. The marker-driven approach involves cloning an open reading frame (ORF) of interest into a vector containing a marker for mycoplasma transformation, as well as sites for transposase loading and random genomic integration. An innovative feature of this construct is to use a single promoter to express the transformation marker and the introduced ORF. The marker-driven approach can be conveniently applied to any exogenous or synthetic gene without any information on the effect of the gene on the strain, whereas the marker-less approach requires that the fragment has a recognizable effect. Using the marker-less method, we found that a region containing the nusG gene rescues a slow growth phenotype of a strain containing a larger deletion encompassing this gene. Using the marker-driven approach, we better defined this finding, thereby establishing that nusG is required for a normal growth rate in synthetic M. mycoides. These methods are suitable for complementation tests to identify genes responsible for assorted functions lacking in deletion mutants. These approaches are also expected to facilitate rapid testing of various natural and engineered genes or gene clusters from numerous sources in M. mycoides.

Transferring whole genomes from bacteria to yeast spheroplasts using entire bacterial cells to reduce DNA shearing.

Direct cell-to-cell transfer of genomes from bacteria to yeast facilitates genome engineering for bacteria that are not amenable to genetic manipulation by allowing instead for the utilization of the powerful yeast genetic tools. Here we describe a protocol for transferring whole genomes from bacterial cells to yeast spheroplasts without any DNA purification process. The method is dependent on the treatment of the bacterial and yeast cellular mixture with PEG, which induces cell fusion, engulfment, aggregation or lysis. Over 80% of the bacterial genomes transferred in this way are complete, on the basis of structural and functional tests. Excluding the time required for preparing starting cultures and for incubating cells to form final colonies, the protocol can be completed in 3 h.

Assembly of eukaryotic algal chromosomes in yeast.

BACKGROUND: Synthetic genomic approaches offer unique opportunities to use powerful yeast and Escherichia coli genetic systems to assemble and modify chromosome-sized molecules before returning the modified DNA to the target host. For example, the entire 1 Mb Mycoplasma mycoides chromosome can be stably maintained and manipulated in yeast before being transplanted back into recipient cells. We have previously demonstrated that cloning in yeast of large (> ~ 150 kb), high G + C (55%) prokaryotic DNA fragments was improved by addition of yeast replication origins every ~100 kb. Conversely, low G + C DNA is stable (up to at least 1.8 Mb) without adding supplemental yeast origins. It has not been previously tested whether addition of yeast replication origins similarly improves the yeast-based cloning of large (>150 kb) eukaryotic DNA with moderate G + C content. The model diatom Phaeodactylum tricornutum has an average G + C content of 48% and a 27.4 Mb genome sequence that has been assembled into chromosome-sized scaffolds making it an ideal test case for assembly and maintenance of eukaryotic chromosomes in yeast. RESULTS: We present a modified chromosome assembly technique in which eukaryotic chromosomes as large as ~500 kb can be assembled from cloned ~100 kb fragments. We used this technique to clone fragments spanning P. tricornutum chromosomes 25 and 26 and to assemble these fragments into single, chromosome-sized molecules. We found that addition of yeast replication origins improved the cloning, assembly, and maintenance of the large chromosomes in yeast. Furthermore, purification of the fragments to be assembled by electroelution greatly increased assembly efficiency. CONCLUSIONS: Entire eukaryotic chromosomes can be successfully cloned, maintained, and manipulated in yeast. These results highlight the improvement in assembly and maintenance afforded by including yeast replication origins in eukaryotic DNA with moderate G + C content (48%). They also highlight the increased efficiency of assembly that can be achieved by purifying fragments before assembly.

Direct transfer of whole genomes from bacteria to yeast.

Transfer of genomes into yeast facilitates genome engineering for genetically intractable organisms, but this process has been hampered by the need for cumbersome isolation of intact genomes before transfer. Here we demonstrate direct cell-to-cell transfer of bacterial genomes as large as 1.8 megabases (Mb) into yeast under conditions that promote cell fusion. Moreover, we discovered that removal of restriction endonucleases from donor bacteria resulted in the enhancement of genome transfer.

Sequence analysis of a complete 1.66 Mb Prochlorococcus marinus MED4 genome cloned in yeast.

Marine cyanobacteria of the genus Prochlorococcus represent numerically dominant photoautotrophs residing throughout the euphotic zones in the open oceans and are major contributors to the global carbon cycle. Prochlorococcus has remained a genetically intractable bacterium due to slow growth rates and low transformation efficiencies using standard techniques. Our recent successes in cloning and genetically engineering the AT-rich, 1.1 Mb Mycoplasma mycoides genome in yeast encouraged us to explore similar methods with Prochlorococcus. Prochlorococcus MED4 has an AT-rich genome, with a GC content of 30.8%, similar to that of Saccharomyces cerevisiae (38%), and contains abundant yeast replication origin consensus sites (ACS) evenly distributed around its 1.66 Mb genome. Unlike Mycoplasma cells, which use the UGA codon for tryptophane, Prochlorococcus uses the standard genetic code. Despite this, we observed no toxic effects of several partial and 15 whole Prochlorococcus MED4 genome clones in S. cerevisiae. Sequencing of a Prochlorococcus genome purified from yeast identified 14 single base pair missense mutations, one frameshift, one single base substitution to a stop codon and one dinucleotide transversion compared to the donor genomic DNA. We thus provide evidence of transformation, replication and maintenance of this 1.66 Mb intact bacterial genome in S. cerevisiae.

Cloning the Acholeplasma laidlawii PG-8A genome in Saccharomyces cerevisiae as a yeast centromeric plasmid.

Cloning of whole genomes of the genus Mycoplasma in yeast has been an essential step for the creation of the first synthetic cell. The genome of the synthetic cell is based on Mycoplasma mycoides, which deviates from the universal genetic code by encoding tryptophan rather than the UGA stop codon. The feature was thought to be important because bacterial genes might be toxic to the host yeast cell if driven by a cryptic promoter active in yeast. As we move to expand the range of bacterial genomes cloned in yeast, we extended this technology to bacteria that use the universal genetic code. Here we report cloning of the Acholeplasma laidlawii PG-8A genome, which uses the universal genetic code. We discovered that only one A. laidlawii gene, a surface anchored extracellular endonuclease, was toxic when cloned in yeast. This gene was inactivated in order to clone and stably maintain the A. laidlawii genome as a centromeric plasmid in the yeast cell.

Assembly of large, high G+C bacterial DNA fragments in yeast.

The ability to assemble large pieces of prokaryotic DNA by yeast recombination has great application in synthetic biology, but cloning large pieces of high G+C prokaryotic DNA in yeast can be challenging. Additional considerations in cloning large pieces of high G+C DNA in yeast may be related to toxic genes, to the size of the DNA, or to the absence of yeast origins of replication within the sequence. As an example of our ability to clone high G+C DNA in yeast, we chose to work with Synechococcus elongatus PCC 7942, which has an average G+C content of 55%. We determined that no regions of the chromosome are toxic to yeast and that S. elongatus DNA fragments over ~200 kb are not stably maintained. DNA constructs with a total size under 200 kb could be readily assembled, even with 62 kb of overlapping sequence between pieces. Addition of yeast origins of replication throughout allowed us to increase the total size of DNA that could be assembled to at least 454 kb. Thus, cloning strategies utilizing yeast recombination with large, high G+C prokaryotic sequences should include yeast origins of replication as a part of the design process.

Chemical synthesis of the mouse mitochondrial genome.

We describe a one-step, isothermal assembly method for synthesizing DNA molecules from overlapping oligonucleotides. The method cycles between in vitro recombination and amplification until the desired length is reached. As a demonstration of its simplicity and robustness, we synthesized the entire 16.3-kilobase mouse mitochondrial genome from 600 overlapping 60-mers.

Creation of a bacterial cell controlled by a chemically synthesized genome.

We report the design, synthesis, and assembly of the 1.08-mega-base pair Mycoplasma mycoides JCVI-syn1.0 genome starting from digitized genome sequence information and its transplantation into a M. capricolum recipient cell to create new M. mycoides cells that are controlled only by the synthetic chromosome. The only DNA in the cells is the designed synthetic DNA sequence, including "watermark" sequences and other designed gene deletions and polymorphisms, and mutations acquired during the building process. The new cells have expected phenotypic properties and are capable of continuous self-replication.