Evolutionary genome engineering using a restriction-modification system

Recent progress in adaptive laboratory evolution of industrial microorganisms

Adaptive laboratory evolution (ALE) is a technique for the selection of strains with better phenotypes by long-term culture under a specific selection pressure or growth environment. Because ALE does not require detailed knowledge of a variety of complex and interactive metabolic networks, and only needs to simulate natural environmental conditions in the laboratory to design a selection pressure, it has the advantages of broad adaptability, strong practicability, and more convenient transformation of strains. In addition, ALE provides a powerful method for studying the evolutionary forces that change the phenotype, performance, and stability of strains, resulting in more productive industrial strains with beneficial mutations. In recent years, ALE has been widely used in the activation of specific microbial metabolic pathways and phenotypic optimization, the efficient utilization of specific substrates, the optimization of tolerance to toxic substance, and the biosynthesis of target products, which is more conducive to the production of industrial strains with excellent phenotypic characteristics. In this paper, typical examples of ALE applications in the development of industrial strains and the research progress of this technology are reviewed, followed by a discussion of its development prospects.

Genomic Features and Pervasive Negative Selection in Rhodanobacter Strains Isolated from Nitrate and Heavy Metal Contaminated Aquifer

Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation.

IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

Genetic Diversity for Accelerating Microbial Adaptive Laboratory Evolution

Adaptive laboratory evolution (ALE) is a widely used and highly effective tool for improving microbial phenotypes and investigating the evolutionary roots of biological phenomena. Serving as the raw materials of evolution, mutations have been extensively utilized to increase the chances of engineering molecules or microbes with tailor-made functions. The generation of genetic diversity is therefore a core technology for accelerating ALE, and a high-quality mutant library is crucial to its success. Because of its importance, technologies for generating genetic diversity have undergone rapid development in recent years. Here, we review the existing techniques for the construction of mutant libraries, briefly introduce their mechanisms and applications, discuss ongoing and emerging efforts to apply engineering technologies in the construction of mutant libraries, and suggest future perspectives for library construction.

MamA essentiality in Mycobacterium smegmatis is explained by the presence of an apparent cognate restriction endonuclease

OBJECTIVE

Restriction-Modification (R-M) systems are ubiquitous in bacteria and were considered for years as rudimentary immune systems that protect bacterial cells from foreign DNA. Currently, these R-M systems are recognized as important players in global gene expression and other cellular processes such us virulence and evolution of genomes. Here, we report the role of the unique DNA methyltransferase in Mycobacterium smegmatis, which shows a moderate degree of sequence similarity to MamA, a previously characterized methyltransferase that affects gene expression in Mycobacterium tuberculosis and is important for survival under hypoxic conditions.

RESULTS

We found that depletion of mamA levels impairs growth and produces elongated cell bodies. Microscopy revealed irregular septation and unevenly distributed DNA, with large areas devoid of DNA and small DNA-free cells. Deletion of MSMEG_3214, a predicted endonuclease-encoding gene co-transcribed with mamA, restored the WT growth phenotype in a mamA-depleted background. Our results suggest that the mamA-depletion phenotype can be explained by DNA cleavage by the apparent cognate restriction endonuclease MSMEG_3214. In addition, in silico analysis predicts that both MamA methyltransferase and MSMEG_3214 endonuclease recognize the same palindromic DNA sequence. We propose that MamA and MSMEG_3214 constitute a previously undescribed R-M system in M. smegmatis.

Functional Traits Co-Occurring with Mobile Genetic Elements in the Microbiome of the Atacama Desert

Mobile genetic elements (MGEs) play an essential role in bacterial adaptation and evolution. These elements are enriched within bacterial communities from extreme environments. However, very little is known if specific genes co-occur with MGEs in extreme environments and, if so, what their function is. We used shotgun-sequencing to analyse the metagenomes of 12 soil samples and characterized the composition of MGEs and the genes co-occurring with them. The samples ranged from less arid coastal sites to the inland hyperarid core of the Atacama Desert, as well as from sediments below boulders, protected from UV-irradiation. MGEs were enriched at the hyperarid sites compared with sediments from below boulders and less arid sites. MGEs were mostly co-occurring with genes belonging to the Cluster Orthologous Group (COG) categories “replication, recombination and repair,” “transcription” and “signal transduction mechanisms.” In general, genes coding for transcriptional regulators and histidine kinases were the most abundant genes proximal to MGEs. Genes involved in energy production were significantly enriched close to MGEs at the hyperarid sites. For example, dehydrogenases, reductases, hydrolases and chlorite dismutase and other enzymes linked to nitrogen metabolism such as nitrite- and nitro-reductase. Stress response genes, including genes involved in antimicrobial and heavy metal resistance genes, were rarely found near MGEs. The present study suggests that MGEs could play an essential role in the adaptation of the soil microbiome in hyperarid desert soils by the modulation of housekeeping genes such as those involved in energy production.

Semi-supervised prediction of gene regulatory networks using machine learning algorithms

Use of computational methods to predict gene regulatory networks (GRNs) from gene expression data is a challenging task. Many studies have been conducted using unsupervised methods to fulfill the task; however, such methods usually yield low prediction accuracies due to the lack of training data. In this article, we propose semi-supervised methods for GRN prediction by utilizing two machine learning algorithms, namely, support vector machines (SVM) and random forests (RF). The semi-supervised methods make use of unlabelled data for training. We investigated inductive and transductive learning approaches, both of which adopt an iterative procedure to obtain reliable negative training data from the unlabelled data. We then applied our semi-supervised methods to gene expression data of Escherichia coli and Saccharomyces cerevisiae, and evaluated the performance of our methods using the expression data. Our analysis indicated that the transductive learning approach outperformed the inductive learning approach for both organisms. However, there was no conclusive difference identified in the performance of SVM and RF. Experimental results also showed that the proposed semi-supervised methods performed better than existing supervised methods for both organisms.

Role of Restriction-Modification Systems in Prokaryotic Evolution and Ecology

Restriction-modification (R-M) systems are able to methylate or cleave DNA depending on methylation status of their recognition site. It allows them to protect bacterial cells from invasion by foreign DNA. Comparative analysis of a large number of available bacterial genomes and methylomes clearly demonstrates that the role of R-M systems in bacteria is wider than only defense. R-M systems maintain heterogeneity of a bacterial population and are involved in adaptation of bacteria to change in their environmental conditions. R-M systems can be essential for host colonization by pathogenic bacteria. Phase variation and intragenomic recombinations are sources of the fast evolution of the specificity of R-M systems. This review focuses on the influence of R-M systems on evolution and ecology of prokaryotes.

Restriction-Modification Systems as Mobile Genetic Elements in the Evolution of an Intracellular Symbiont

Long-term vertical transmission of intracellular bacteria causes massive genomic erosion and results in extremely small genomes, particularly in ancient symbionts. Genome reduction is typically preceded by the accumulation of pseudogenes and proliferation of mobile genetic elements, which are responsible for chromosome rearrangements during the initial stage of endosymbiosis. We compared the genomes of an endosymbiont of termite gut flagellates, “Candidatus Endomicrobium trichonymphae,” and its free-living relative Endomicrobium proavitum and discovered many remnants of restriction-modification (R-M) systems that are consistently associated with genome rearrangements in the endosymbiont genome. The rearrangements include apparent insertions, transpositions, and the duplication of a genomic region; there was no evidence of transposon structures or other mobile elements. Our study reveals a so far unrecognized mechanism for genome rearrangements in intracellular symbionts and sheds new light on the general role of R-M systems in genome evolution.

To be or not to be: regulation of restriction-modification systems and other toxin-antitoxin systems

One of the simplest classes of genes involved in programmed death is that containing the toxin–antitoxin (TA) systems of prokaryotes. These systems are composed of an intracellular toxin and an antitoxin that neutralizes its effect. These systems, now classified into five types, were initially discovered because some of them allow the stable maintenance of mobile genetic elements in a microbial population through postsegregational killing or the death of cells that have lost these systems. Here, we demonstrate parallels between some TA systems and restriction–modification systems (RM systems). RM systems are composed of a restriction enzyme (toxin) and a modification enzyme (antitoxin) and limit the genetic flux between lineages with different epigenetic identities, as defined by sequence-specific DNA methylation. The similarities between these systems include their postsegregational killing and their effects on global gene expression. Both require the finely regulated expression of a toxin and antitoxin. The antitoxin (modification enzyme) or linked protein may act as a transcriptional regulator. A regulatory antisense RNA recently identified in an RM system can be compared with those RNAs in TA systems. This review is intended to generalize the concept of TA systems in studies of stress responses, programmed death, genetic conflict and epigenetics.

Diverse functions of restriction-modification systems in addition to cellular defense

Restriction-modification (R-M) systems are ubiquitous and are often considered primitive immune systems in bacteria. Their diversity and prevalence across the prokaryotic kingdom are an indication of their success as a defense mechanism against invading genomes. However, their cellular defense function does not adequately explain the basis for their immaculate specificity in sequence recognition and nonuniform distribution, ranging from none to too many, in diverse species. The present review deals with new developments which provide insights into the roles of these enzymes in other aspects of cellular function. In this review, emphasis is placed on novel hypotheses and various findings that have not yet been dealt with in a critical review. Emerging studies indicate their role in various cellular processes other than host defense, virulence, and even controlling the rate of evolution of the organism. We also discuss how R-M systems could have successfully evolved and be involved in additional cellular portfolios, thereby increasing the relative fitness of their hosts in the population.

Conservation of Dcm-mediated cytosine DNA methylation in Escherichia coli

In Escherichia coli, cytosine DNA methylation is catalyzed by the DNA cytosine methyltransferase (Dcm) protein and occurs at the second cytosine in the sequence 5'CCWGG3'. Although the presence of cytosine DNA methylation was reported over 35 years ago, the biological role of 5-methylcytosine in E. coli remains unclear. To gain insight into the role of cytosine DNA methylation in E. coli, we (1) screened the 72 strains of the ECOR collection and 90 recently isolated environmental samples for the presence of the full-length dcm gene using the polymerase chain reaction; (2) examined the same strains for the presence of 5-methylcytosine at 5'CCWGG3' sites using a restriction enzyme isoschizomer digestion assay; and (3) quantified the levels of 5-methyl-2'-deoxycytidine in selected strains using liquid chromatography tandem mass spectrometry. Dcm-mediated cytosine DNA methylation is conserved in all 162 strains examined, and the level of 5-methylcytosine ranges from 0.86% to 1.30% of the cytosines. We also demonstrate that Dcm reduces the expression of ribosomal protein genes during stationary phase, and this may explain the highly conserved nature of this DNA modification pathway.