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.

Cerebral small vessel disease burden and longitudinal cognitive decline from age 73 to 82: the Lothian Birth Cohort 1936.

Slowed processing speed is considered a hallmark feature of cognitive decline in cerebral small vessel disease (SVD); however, it is unclear whether SVD's association with slowed processing might be due to its association with overall declining general cognitive ability. We quantified the total MRI-visible SVD burden of 540 members of the Lothian Birth Cohort 1936 (age: 72.6 ± 0.7 years; 47% female). Using latent growth curve modelling, we tested associations between total SVD burden at mean age 73 and changes in general cognitive ability, processing speed, verbal memory and visuospatial ability, measured at age 73, 76, 79 and 82. Covariates included age, sex, vascular risk and childhood cognitive ability. In the fully adjusted models, greater SVD burden was associated with greater declines in general cognitive ability (standardised ß: -0.201; 95% CI: [-0.36, -0.04]; pFDR = 0.022) and processing speed (-0.222; [-0.40, -0.04]; pFDR = 0.022). SVD burden accounted for between 4 and 5% of variance in declines of general cognitive ability and processing speed. After accounting for the covariance between tests of processing speed and general cognitive ability, only SVD's association with greater decline in general cognitive ability remained significant, prior to FDR correction (-0.222; [-0.39, -0.06]; p = 0.008; pFDR = 0.085). Our findings do not support the notion that SVD has a specific association with declining processing speed, independent of decline in general cognitive ability (which captures the variance shared across domains of cognitive ability). The association between SVD burden and declining general cognitive ability supports the notion of SVD as a diffuse, whole-brain disease and suggests that trials monitoring SVD-related cognitive changes should consider domain-specific changes in the context of overall, general cognitive decline.

Gross Chromosomal Rearrangements in Kluyveromyces marxianus Revealed by Illumina and Oxford Nanopore Sequencing.

Kluyveromyces marxianus (K. marxianus) is an increasingly popular industrially relevant yeast. It is known to possess a highly efficient non-homologous end joining (NHEJ) pathway that promotes random integration of non-homologous DNA fragments into its genome. The nature of the integration events was traditionally analyzed by Southern blot hybridization. However, the precise DNA sequence at the insertion sites were not fully explored. We transformed a PCR product of the Saccharomyces cerevisiae URA3 gene (ScURA3) into an uracil auxotroph K. marxianus otherwise wildtype strain and picked 24 stable Ura+ transformants for sequencing analysis. We took advantage of rapid advances in DNA sequencing technologies and developed a method using a combination of Illumina MiSeq and Oxford Nanopore sequencing. This approach enables us to uncover the gross chromosomal rearrangements (GCRs) that are associated with the ScURA3 random integration. Moreover, it will shine a light on understanding DNA repair mechanisms in eukaryotes, which could potentially provide insights for cancer research.

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.

Using a Robust and Sensitive GFP-Based cGMP Sensor for Real-Time Imaging in Intact Caenorhabditis elegans.

cGMP plays a role in sensory signaling and plasticity by regulating ion channels, phosphodiesterases, and kinases. Studies that primarily used genetic and biochemical tools suggest that cGMP is spatiotemporally regulated in multiple sensory modalities. FRET- and GFP-based cGMP sensors were developed to visualize cGMP in primary cell culture and Caenorhabditis elegans to corroborate these findings. While a FRET-based sensor has been used in an intact animal to visualize cGMP, the requirement of a multiple emission system limits its ability to be used on its own as well as with other fluorophores. Here, we demonstrate that a C. elegans codon-optimized version of the cpEGFP-based cGMP sensor FlincG3 can be used to visualize rapidly changing cGMP levels in living, behaving C. elegans We coexpressed FlincG3 with the blue-light-activated guanylyl cyclases BeCyclOp and bPGC in body wall muscles, and found that the rate of change in FlincG3 fluorescence correlated with the rate of cGMP production by each cyclase. Furthermore, we show that FlincG3 responds to cultivation temperature, NaCl concentration changes, and sodium dodecyl sulfate in the sensory neurons AFD, ASEL/R, and PHB, respectively. Intriguingly, FlincG3 fluorescence in ASEL and ASER decreased in response to a NaCl concentration upstep and downstep, respectively, which is opposite in sign to the coexpressed calcium sensor jRGECO1a and previously published calcium recordings. These results illustrate that FlincG3 can be used to report rapidly changing cGMP levels in an intact animal, and that the reporter can potentially reveal unexpected spatiotemporal landscapes of cGMP in response to stimuli.

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.

Direct Transfer of a Mycoplasma mycoides Genome to Yeast Is Enhanced by Removal of the Mycoides Glycerol Uptake Factor Gene glpF.

We previously discovered that intact bacterial chromosomes can be directly transferred to a yeast host cell where they can propagate as centromeric plasmids by fusing bacterial cells with S accharomyces cerevisiae spheroplasts. Inside the host any desired number of genetic changes can be introduced into the yeast centromeric plasmid to produce designer genomes that can be brought to life using a genome transplantation protocol. Earlier research demonstrated that the removal of restriction-systems from donor bacteria, such as Mycoplasma mycoides, Mycoplasma capricolum, or Haemophilus influenzae increased successful genome transfers. These findings suggested that other genetic factors might also impact the bacteria-to-yeast genome transfer process. In this study, we demonstrated that the removal of a particular genetic factor, the glycerol uptake facilitator protein gene glpF from M. mycoides, significantly increased direct genome transfer by up to 21-fold. Additionally, we showed that intact bacterial cells were endocytosed by yeast spheroplasts producing organelle-like structures within these yeast cells. These might lead to the possibility of creating novel synthetic organelles.

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.

One step engineering of the small-subunit ribosomal RNA using CRISPR/Cas9.

Bacteria are indispensable for the study of fundamental molecular biology processes due to their relatively simple gene and genome architecture. The ability to engineer bacterial chromosomes is quintessential for understanding gene functions. Here we demonstrate the engineering of the small-ribosomal subunit (16S) RNA of Mycoplasma mycoides, by combining the CRISPR/Cas9 system and the yeast recombination machinery. We cloned the entire genome of M. mycoides in yeast and used constitutively expressed Cas9 together with in vitro transcribed guide-RNAs to introduce engineered 16S rRNA genes. By testing the function of the engineered 16S rRNA genes through genome transplantation, we observed surprising resilience of this gene to addition of genetic elements or helix substitutions with phylogenetically-distant bacteria. While this system could be further used to study the function of the 16S rRNA, one could envision the "simple" M. mycoides genome being used in this setting to study other genetic structures and functions to answer fundamental questions of life.

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.

Elevated IL-33 expression is associated with pediatric eosinophilic esophagitis, and exogenous IL-33 promotes eosinophilic esophagitis development in mice.

We tested whether the T helper (Th) type 2 (Th2) cell agonist and allergenic ligand IL-33 was associated with eosinophilic esophagitis (EoE) development in a pediatric cohort and whether IL-33 protein could induce disease symptoms in mice. Biopsies from EoE patients or controls were used to measure IL-33 mRNA and protein expression. Increased expression of IL-33 mRNA was found in the esophageal mucosa in EoE. IL-33 protein was detected in cells negative for CD45, mast cells, and epithelial cell markers near blood vessels. Circulating levels of IL-33 were not increased. The time course for IL-33 gene expression was quantified in an established Aspergillus fumigatus allergen mouse model of EoE. Because IL-33 induction was transient in this model and chronicity of IL-33 expression has been demonstrated in humans, naive mice were treated with recombinant IL-33 for 1 wk and esophageal pathology was evaluated. IL-33 application produced changes consistent with phenotypically early EoE, including transmural eosinophilia, mucosal hyperproliferation, and upregulation of eosinophilic genes and chemokines. Th2 cytokines, including IL-13, along with innate lymphoid cell group 2, Th1/17, and M2 macrophage marker genes, were increased after IL-33 application. IL-33-induced eosinophilia was ablated in IL-13 null mice. In addition, IL-33 induced a profound inhibition of the regulatory T cell gene signature. We conclude that IL-33 gene expression is associated with pediatric EoE development and that application of recombinant protein in mice phenocopies the early clinical phase of the human disease in an IL-13-dependent manner. IL-33 inhibition of esophageal regulatory T cell function may induce loss of antigenic tolerance, thereby providing a mechanistic rationale for EoE development.

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.

Designed surface residue substitutions in [NiFe] hydrogenase that improve electron transfer characteristics.

Photobiological hydrogen production is an attractive, carbon-neutral means to convert solar energy to hydrogen. We build on previous research improving the Alteromonas macleodii "Deep Ecotype" [NiFe] hydrogenase, and report progress towards creating an artificial electron transfer pathway to supply the hydrogenase with electrons necessary for hydrogen production. Ferredoxin is the first soluble electron transfer mediator to receive high-energy electrons from photosystem I, and bears an electron with sufficient potential to efficiently reduce protons. Thus, we engineered a hydrogenase-ferredoxin fusion that also contained several other modifications. In addition to the C-terminal ferredoxin fusion, we truncated the C-terminus of the hydrogenase small subunit, identified as the available terminus closer to the electron transfer region. We also neutralized an anionic patch surrounding the interface Fe-S cluster to improve transfer kinetics with the negatively charged ferredoxin. Initial screening showed the enzyme tolerated both truncation and charge neutralization on the small subunit ferredoxin-binding face. While the enzyme activity was relatively unchanged using the substrate methyl viologen, we observed a marked improvement from both the ferredoxin fusion and surface modification using only dithionite as an electron donor. Combining ferredoxin fusion and surface charge modification showed progressively improved activity in an in vitro assay with purified enzyme.

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.

A broad survey reveals substitution tolerance of residues ligating FeS clusters in [NiFe] hydrogenase.

BACKGROUND: In order to understand the effects of FeS cluster attachment in [NiFe] hydrogenase, we undertook a study to substitute all 12 amino acid positions normally ligating the three FeS clusters in the hydrogenase small subunit. Using the hydrogenase from Alteromonas macleodii "deep ecotype" as a model, we substituted one of four amino acids (Asp, His, Asn, Gln) at each of the 12 ligating positions because these amino acids are alternative coordinating residues in otherwise conserved-cysteine positions found in a broad survey of NiFe hydrogenase sequences. We also hoped to discover an enzyme with elevated hydrogen evolution activity relative to a previously reported "G1" (H230C/P285C) improved enzyme in which the medial FeS cluster Pro and the distal FeS cluster His were each substituted for Cys. RESULTS: Among all the substitutions screened, aspartic acid substitutions were generally well-tolerated, and examination suggests that the observed deficiency in enzyme activity may be largely due to misprocessing of the small subunit of the enzyme. Alignment of hydrogenase sequences from sequence databases revealed many rare substitutions; the five substitutions present in databases that we tested all exhibited measurable hydrogen evolution activity. Select substitutions were purified and tested, supporting the results of the screening assay. Analysis of these results confirms the importance of small subunit processing. Normalizing activity to quantity of mature small subunit, indicative of total enzyme maturation, weakly suggests an improvement over the "G1" enzyme. CONCLUSIONS: We have comprehensively screened 48 amino acid substitutions of the hydrogenase from A. macleodii "deep ecotype", to understand non-canonical ligations of amino acids to FeS clusters and to improve hydrogen evolution activity of this class of hydrogenase. Our studies show that non-canonical ligations can be functional and also suggests a new limiting factor in the production of active enzyme.