Experimental Challenges for Reduced Genomes: The Cell Model Escherichia coli

Bacterial genome reduction for optimal chassis of synthetic biology: a review

Bacteria with streamlined genomes, that harbor full functional genes for essential metabolic networks, are able to synthesize the desired products more effectively and thus have advantages as production platforms in industrial applications. To obtain streamlined chassis genomes, a large amount of effort has been made to reduce existing bacterial genomes. This work falls into two categories: rational and random reduction. The identification of essential gene sets and the emergence of various genome-deletion techniques have greatly promoted genome reduction in many bacteria over the past few decades. Some of the constructed genomes possessed desirable properties for industrial applications, such as: increased genome stability, transformation capacity, cell growth, and biomaterial productivity. The decreased growth and perturbations in physiological phenotype of some genome-reduced strains may limit their applications as optimized cell factories. This review presents an assessment of the advancements made to date in bacterial genome reduction to construct optimal chassis for synthetic biology, including: the identification of essential gene sets, the genome-deletion techniques, the properties and industrial applications of artificially streamlined genomes, the obstacles encountered in constructing reduced genomes, and the future perspectives.

In pursuit of a minimal CHO genome: Establishment of large-scale genome deletions

Chinese hamster ovary (CHO) cells are the most commonly used mammalian cell line for the production of complex therapeutic glycoproteins. As CHO cells have evolved as part of a multicellular organism, they harbor many cellular functions irrelevant for their application as production hosts in industrial bioprocesses. Consequently, CHO cells have been the target for numerous genetic engineering efforts in the past, but a tailored host cell chassis holistically optimized for its specific task in a bioreactor is still missing. While the concept of genome reduction has already been successfully applied to bacterial production cells, attempts to create higher eukaryotic production hosts exhibiting reduced genomes have not been reported yet. Here, we present the establishment and application of a large-scale genome deletion strategy for targeted excision of large genomic regions in CHO cells. We demonstrate the feasibility of genome reduction in CHO cells using optimized CRISPR/Cas9 based experimental protocols targeting large non-essential genomic regions with high efficiency. Achieved genome deletions of non-essential genetic regions did not introduce negative effects on bioprocess relevant parameters, although we conducted the largest reported genomic excision of 864 kilobase pairs in CHO cells so far. The concept presented serves as a directive to accelerate the development of a significantly genome-reduced CHO host cell chassis paving the way for a next generation of CHO cell factories.

Reduction-to-synthesis: the dominant approach to genome-scale synthetic biology

Advances in systems and synthetic biology have propelled the construction of reduced bacterial genomes. Genome reduction was initially focused on exploring properties of minimal genomes, but more recently it has been deployed as an engineering strategy to enhance strain performance. This review provides the latest updates on reduced genomes, focusing on dual-track approaches of top-down reduction and bottom-up synthesis for their construction. Using cases from studies that are based on established industrial workhorse strains, we discuss the construction of a series of synthetic phenotypes that are candidates for biotechnological applications. Finally, we address the possible uses of reduced genomes for biotechnological applications and the needed future research directions that may ultimately lead to the total synthesis of rationally designed genomes.

Comparison of wild-type KT2440 and genome-reduced EM42 Pseudomonas putida strains for muconate production from aromatic compounds and glucose

Pseudomonas putida KT2440 is a robust, aromatic catabolic bacterium that has been widely engineered to convert bio-based and waste-based feedstocks to target products. Towards industrial domestication of P. putida KT2440, rational genome reduction has been previously conducted, resulting in P. putida strain EM42, which exhibited characteristics that could be advantageous for production strains. Here, we compared P. putida KT2440- and EM42-derived strains for cis,cis-muconic acid production from an aromatic compound, p-coumarate, and in separate strains, from glucose. To our surprise, the EM42-derived strains did not outperform the KT2440-derived strains in muconate production from either substrate. In bioreactor cultivations, KT2440- and EM42-derived strains produced muconate from p-coumarate at titers of 45 g/L and 37 g/L, respectively, and from glucose at 20 g/L and 13 g/L, respectively. To provide additional insights about the differences in the parent strains, we analyzed growth profiles of KT2440 and EM42 on aromatic compounds as the sole carbon and energy sources. In general, the EM42 strain exhibited reduced growth rates but shorter growth lags than KT2440. We also observed that EM42-derived strains resulted in higher growth rates on glucose compared to KT2440-derived strains, but only at the lowest glucose concentrations tested. Transcriptomics revealed that genome reduction in EM42 had global effects on transcript levels and showed that the EM42-derived strains that produce muconate from glucose exhibit reduced modulation of gene expression in response to changes in glucose concentrations. Overall, our results highlight that additional studies are warranted to understand the effects of genome reduction on microbial metabolism and physiology, especially when intended for use in production strains.

Complete sequence of the genome-reduced Escherichia coli DGF-298

We report the complete genome sequence and annotation of Escherichia coli DGF-298, a genome-reduced E. coli strain with interesting properties for systems and synthetic biology. DGF-298 has a single circular chromosome of 2,991,126 bp and 2,831 genes, including 2,691 coding sequences, with a mean G + C content of ~51%.

The genetics of aerotolerant growth in an alphaproteobacterium with a naturally reduced genome

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living alphaproteobacterium Zymomonas mobilis has a naturally reduced genome containing fewer than 2,000 protein-coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model alphaproteobacterium. Here, we use genome-scale CRISPRi-seq (clustered regularly interspaced short palindromic repeats interference sequencing) to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most "generally essential" genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as a reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rhodobacter nitrogen fixation (Rnf) respiratory complex are required for the anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts. IMPORTANCE The inherent complexity of biological systems is a major barrier to our understanding of cellular physiology. Bacteria with markedly fewer genes than their close relatives, or reduced genome bacteria, are promising biological models with less complexity. Reduced genome bacteria can also have superior properties for industrial use, provided the reduction does not overly restrict strain robustness. Naturally reduced genome bacteria, such as the alphaproteobacterium Zymomonas mobilis, have fewer genes but remain environmentally robust. In this study, we show that Z. mobilis is a simplified genetic model for Alphaproteobacteria, a class with important impacts on the environment, human health, and industry. We also identify genes that are only required in the absence of atmospheric oxygen, uncovering players that maintain and utilize the cellular energy state. Our findings have broad implications for the genetics of Alphaproteobacteria and industrial use of Z. mobilis to create biofuels and bioproducts.

Efforts to Minimise the Bacterial Genome as a Free-Living Growing System

Exploring the minimal genetic requirements for cells to maintain free living is an exciting topic in biology. Multiple approaches are employed to address the question of the minimal genome. In addition to constructing the synthetic genome in the test tube, reducing the size of the wild-type genome is a practical approach for obtaining the essential genomic sequence for living cells. The well-studied Escherichia coli has been used as a model organism for genome reduction owing to its fast growth and easy manipulation. Extensive studies have reported how to reduce the bacterial genome and the collections of genomic disturbed strains acquired, which were sufficiently reviewed previously. However, the common issue of growth decrease caused by genetic disturbance remains largely unaddressed. This mini-review discusses the considerable efforts made to improve growth fitness, which was decreased due to genome reduction. The proposal and perspective are clarified for further accumulated genetic deletion to minimise the Escherichia coli genome in terms of genome reduction, experimental evolution, medium optimization, and machine learning.

Mutations identified in engineered Escherichia coli with a reduced genome

Characterizing genes that regulate cell growth and survival in model organisms is important for understanding higher organisms. Construction of strains harboring large deletions in the genome can provide insights into the genetic basis of cell growth compared with only studying wild-type strains. We have constructed a series of genome-reduced strains with deletions spanning approximately 38.9% of the E. coli chromosome. Strains were constructed by combining large deletions in chromosomal regions encoding nonessential gene groups. We also isolated strains Δ33b and Δ37c, whose growth was partially restored by adaptive laboratory evolution (ALE). Genome sequencing of nine strains, including those selected following ALE, identified the presence of several Single Nucleotide Variants (SNVs), insertions, deletions, and inversions. In addition to multiple SNVs, two insertions were identified in ALE strain Δ33b. The first was an insertion at the promoter region of pntA, which increased cognate gene expression. The second was an insertion sequence (IS) present in sibE, encoding the antitoxin in a toxin-antitoxin system, which decreased expression of sibE. 5 strains of Δ37c independently isolated following ALE harboring multiple SNVs and genetic rearrangements. Interestingly, a SNV was identified in the promoter region of hcaT in all five strains, which increased hcaT expression and, we predict, rescued the attenuated Δ37b growth. Experiments using defined deletion mutants suggested that hcaT encodes a 3-phenylpropionate transporter protein and is involved in survival during stationary phase under oxidative stress. This study is the first to document accumulation of mutations during construction of genome-reduced strains. Furthermore, isolation and analysis of strains derived from ALE in which the growth defect mediated by large chromosomal deletions was rescued identified novel genes involved in cell survival.

Systematic Large Fragment Deletions in the Genome of Synechococcus elongatus and the Consequent Changes in Transcriptomic Profiles

Genome streamlining, as a natural process in the evolution of microbes, has become a common approach for generating ideal chassis cells for synthetic biology studies and industrial applications. However, systematic genome reduction remains a bottleneck in the generation of such chassis cells with cyanobacteria, due to very time-consuming genetic manipulations. Synechococcus elongatus PCC 7942, a unicellular cyanobacterium, is a candidate for systematic genome reduction, as its essential and nonessential genes have been experimentally identified. Here, we report that at least 20 of the 23 over 10 kb nonessential gene regions could be deleted and that stepwise deletions of these regions could be achieved. A septuple-deletion mutant (genome reduced by 3.8%) was generated, and the effects of genome reduction on the growth and genome-wide transcription were investigated. In the ancestral triple to sextuple mutants (b, c, d, e1), an increasingly large number of genes (up to 998) were upregulated relative to the wild type, while slightly fewer genes (831) were upregulated in the septuple mutant (f). In a different sextuple mutant (e2) derived from the quintuple mutant d, much fewer genes (232) were upregulated. Under the standard conditions in this study, the mutant e2 showed a higher growth rate than the wild type, e1 and f. Our results indicate that it is feasible to extensively reduce the genomes of cyanobacteria for generation of chassis cells and for experimental evolutionary studies.

Growth rate-associated transcriptome reorganization in response to genomic, environmental, and evolutionary interruptions

The genomic, environmental, and evolutionary interruptions caused the changes in bacterial growth, which were stringently associated with changes in gene expression. The growth and gene expression changes remained unclear in response to these interruptions that occurred combinative. As a pilot study, whether and how bacterial growth was affected by the individual and dual interruptions of genome reduction, environmental stress, and adaptive evolution were investigated. Growth assay showed that the presence of the environmental stressors, i.e., threonine and chloramphenicol, significantly decreased the growth rate of the wild-type Escherichia coli, whereas not that of the reduced genome. It indicated a canceling effect in bacterial growth due to the dual interruption of the genomic and environmental changes. Experimental evolution of the reduced genome released the canceling effect by improving growth fitness. Intriguingly, the transcriptome architecture maintained a homeostatic chromosomal periodicity regardless of the genomic, environmental, and evolutionary interruptions. Negative epistasis in transcriptome reorganization was commonly observed in response to the dual interruptions, which might contribute to the canceling effect. It was supported by the changes in the numbers of differentially expressed genes (DEGs) and the enriched regulons and functions. Gene network analysis newly constructed 11 gene modules, one out of which was correlated to the growth rate. Enrichment of DEGs in these modules successfully categorized them into three types, i.e., conserved, responsive, and epistatic. Taken together, homeostasis in transcriptome architecture was essential to being alive, and it might be attributed to the negative epistasis in transcriptome reorganization and the functional differentiation in gene modules. The present study directly connected bacterial growth fitness with transcriptome reorganization and provided a global view of how microorganisms responded to genomic, environmental, and evolutionary interruptions for survival from wild nature.

Biocircuits in plants and eukaryotic algae

As one of synthetic biology's foundations, biocircuits are a strategy of genetic parts assembling to recognize a signal and to produce a desirable output to interfere with a biological function. In this review, we revisited the progress in the biocircuits technology basis and its mandatory elements, such as the characterization and assembly of functional parts. Furthermore, for a successful implementation, the transcriptional control systems are a relevant point, and the computational tools help to predict the best combinations among the biological parts planned to be used to achieve the desirable phenotype. However, many challenges are involved in delivering and stabilizing the synthetic structures. Some research experiences, such as the golden crops, biosensors, and artificial photosynthetic structures, can indicate the positive and limiting aspects of the practice. Finally, we envision that the modulatory structural feature and the possibility of finer gene regulation through biocircuits can contribute to the complex design of synthetic chromosomes aiming to develop plants and algae with new or improved functions.

Global coordination of the mutation and growth rates across the genetic and nutritional variety in Escherichia coli

Fitness and mutability are the primary traits of living organisms for adaptation and evolution. However, their quantitative linkage remained largely deficient. Whether there is any general relationship between the two features and how genetic and environmental variables influence them remained unclear and were addressed here. The mutation and growth rates of an assortment of Escherichia coli strain collections, including the wild-type strains and the genetically disturbed strains of either reduced genomes or deletion of the genes involved in the DNA replication fidelity, were evaluated in various media. The contribution of media to the mutation and growth rates was differentiated depending on the types of genetic disturbance. Nevertheless, the negative correlation between the mutation and growth rates was observed across the genotypes and was common in all media. It indicated the comprehensive association of the correlated mutation and growth rates with the genetic and medium variation. Multiple linear regression and support vector machine successfully predicted the mutation and growth rates and the categories of genotypes and media, respectively. Taken together, the study provided a quantitative dataset linking the mutation and growth rates, genotype, and medium and presented a simple and successful example of predicting bacterial growth and mutability by data-driven approaches.

Random genome reduction coupled with polyhydroxybutyrate biosynthesis to facilitate its accumulation in Escherichia coli

Genome reduction has been emerged as a powerful tool to construct ideal chassis for synthetic biology. Random genome reduction couple genomic deletion with growth and has the potential to construct optimum genome for a given environment. Recently, we developed a transposon-mediated random deletion (TMRD) method that allows the random and continuous reduction of Escherichia coli genome. Here, to prove its ability in constructing optimal cell factories, we coupled polyhydroxybutyrate (PHB) accumulation with random genome reduction and proceeded to reduce the E. coli genome. Five mutants showed high biomass and PHB yields were selected from 18 candidates after ten rounds of genome reduction. And eight or nine genomic fragments (totally 230.1–270.0 Kb) were deleted in their genomes, encompassing 4.95%–5.82% of the parental MG1655 genome. Most mutants displayed better growth, glucose utilization, protein expression, and significant increase of electroporation efficiency compared with MG1655. The PHB content and concentration enhanced up to 13.3%–37.2% and 60.2%–102.9% when batch fermentation was performed in M9-glucose medium using the five mutants. Particularly, in mutant H16, lacking 5.28% of its genome, the increase of biomass and PHB concentration were more than 50% and 100% compared with MG1655, respectively. This work expands the strategy for creating streamlined chassis to improve the production of high value-added products.

Data Incompleteness May form a Hard-to-Overcome Barrier to Decoding Life’s Mechanism

Simple Summary

The influence of data incompleteness on the correctness of conclusions about the structure and functions of the objects under study is widely discussed in the literature. It was noted that even a small percentage of missing data can lead to incorrect conclusions and imperfect knowledge. In particular, incompleteness can lead to critical errors in the qualitative and quantitative assessments of interactions in biological systems and a distorted understanding of the functioning mechanisms of living systems. In this brief review, we attempt to demonstrate the extent of this incompleteness in functional information about living systems using the best-studied examples. We suggest that this incompleteness may form seemingly insurmountable barriers in deciphering the mechanisms of the functioning of complex systems with unpredictable properties arising from the interaction of the system components.

Abstract

In this brief review, we attempt to demonstrate that the incompleteness of data, as well as the intrinsic heterogeneity of biological systems, may form very strong and possibly insurmountable barriers for researchers trying to decipher the mechanisms of the functioning of live systems. We illustrate this challenge using the two most studied organisms: E. coli, with 34.6% genes lacking experimental evidence of function, and C. elegans, with identified proteins for approximately 50% of its genes. Another striking example is an artificial unicellular entity named JCVI-syn3.0, with a minimal set of genes. A total of 31.5% of the genes of JCVI-syn3.0 cannot be ascribed a specific biological function. The human interactome mapping project identified only 5–10% of all protein interactions in humans. In addition, most of the available data are static snapshots, and it is barely possible to generate realistic models of the dynamic processes within cells. Moreover, the existing interactomes reflect the de facto interaction but not its functional result, which is an unpredictable emerging property. Perhaps the completeness of molecular data on any living organism is beyond our reach and represents an unsolvable problem in biology.

Reduction of the Bacterial Genome by Transposon-Mediated Random Deletion

Genome reduction is an important strategy in synthetic biology for constructing functional chassis cells or minimal genomes. However, the limited knowledge of complex gene functions and interactions makes genome reduction by rational design encounter a bottleneck. Here, we present an iterative and random genome reduction method for Escherichia coli, named "transposon-mediated random deletion (TMRD)". TMRD generates random double-strand breaks (DSBs) in the genome by combining Tn5 transposition with the CRISPR/Cas9 system and allows genomic deletions of various sizes at random positions during DSB repair through the intracellular alternative end-joining mechanism. Using E. coli MG1655 as the original strain, a pool of cells with multiple random genomic deletions were obtained after five reduction cycles. The growth rates of the obtained cells were comparable to that of MG1655, while the electroporation efficiency increased by at least 2 magnitudes. TMRD can generate a small E. coli library carrying multiple and random genomic deletions while enriching the cells with environmental fitness in the population. TMRD has the potential to be widely applied in the construction of minimal genomes or chassis cells for metabolic engineering.

Improving the Intensity of Integrated Expression for Microbial Production

Chromosomal integration of exogenous genes is preferred for industrially related fermentation, as plasmid-mediated fermentation leads to extra metabolic burden and genetic instability. Moreover, with the development and advancement of genome engineering and gene editing technologies, inserting genes into chromosomes has become more convenient; integration expression is extensively utilized in microorganisms for industrial bioproduction and expected to become the trend of recombinant protein expression. However, in actual research and application, it is important to enhance the expression of heterologous genes at the host genome level. Herein, we summarized the basic principles and characteristics of genomic integration; furthermore, we highlighted strategies to improve the expression of chromosomal integration of genes and pathways in host strains from three aspects, including chassis cell optimization, regulation of expression elements in gene expression cassettes, optimization of gene dose level and integration sites on chromosomes. Moreover, we reviewed and summarized the relevant studies on the application of integrated expression in the exploration of gene function and the various types of industrial microorganism production. Consequently, this review would serve as a reference for the better application of integrated expression.

Correlated chromosomal periodicities according to the growth rate and gene expression

Linking genetic information to population fitness is crucial to understanding living organisms. Despite the abundant knowledge of the genetic contribution to growth, the overall patterns/features connecting genes, their expression, and growth remain unclear. To reveal the quantitative and direct connections, systematic growth assays of single-gene knockout Escherichia coli strains under both rich and poor nutritional conditions were performed; subsequently, the resultant growth rates were associated with the original expression levels of the knockout genes in the parental genome. Comparative analysis of growth and the transcriptome identified not only the nutritionally differentiated fitness cost genes but also a significant correlation between the growth rates of the single-gene knockout strains and the original expression levels of these knockout genes in the parental strain, regardless of the nutritional variation. In addition, the coordinated chromosomal periodicities of the wild-type transcriptome and the growth rates of the strains lacking the corresponding genes were observed. The common six-period periodicity was somehow attributed to the essential genes, although the underlying mechanism remains to be addressed. The correlated chromosomal periodicities associated with the gene expression-growth dataset were highly valuable for bacterial growth prediction and discovering the working principles governing minimal genetic information.

Plasmid DNA Production in Proteome-Reduced Escherichia coli

The design of optimal cell factories requires engineering resource allocation for maximizing product synthesis. A recently developed method to maximize the saving in cell resources released 0.5% of the proteome of Escherichia coli by deleting only three transcription factors. We assessed the capacity for plasmid DNA (pDNA) production in the proteome-reduced strain in a mineral medium, lysogeny, and terrific broths. In all three cases, the pDNA yield from biomass was between 33 and 53% higher in the proteome-reduced than in its wild type strain. When cultured in fed-batch mode in shake-flask, the proteome-reduced strain produced 74.8 mg L^-1 pDNA, which was four times greater than its wild-type strain. Nevertheless, the pDNA supercoiled fraction was less than 60% in all cases. Deletion of recA increased the pDNA yields in the wild type, but not in the proteome-reduced strain. Furthermore, recA mutants produced a higher fraction of supercoiled pDNA, compared to their parents. These results show that the novel proteome reduction approach is a promising starting point for the design of improved pDNA production hosts.

Clustering of Bacterial Growth Dynamics in Response to Growth Media by Dynamic Time Warping

Bacterial growth curves, representing population dynamics, are still poorly understood. The growth curves are commonly analyzed by model-based theoretical fitting, which is limited to typical S-shape fittings and does not elucidate the dynamics in their entirety. Thus, whether a certain growth condition results in any particular pattern of growth curve remains unclear. To address this question, up-to-date data mining techniques were applied to bacterial growth analysis for the first time. Dynamic time warping (DTW) and derivative DTW (DDTW) were used to compare the similarity among 1015 growth curves of 28 Escherichia coli strains growing in three different media. In the similarity evaluation, agglomerative hierarchical clustering, assessed with four statistic benchmarks, successfully categorized the growth curves into three clusters, roughly corresponding to the three media. Furthermore, a simple benchmark was newly proposed, providing a highly improved accuracy (~99%) in clustering the growth curves corresponding to the growth media. The biologically reasonable categorization of growth curves suggested that DTW and DDTW are applicable for bacterial growth analysis. The bottom-up clustering results indicate that the growth media determine some specific patterns of population dynamics, regardless of genomic variation, and thus have a higher priority of shaping the growth curves than the genomes do.