Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing

Improving the genetic stability of bacterial growth control for long-term bioproduction

Using microorganisms for bioproduction requires the reorientation of metabolic fluxes from biomass synthesis to the production of compounds of interest. We previously engineered a synthetic growth switch in Escherichia coli based on inducible expression of the β- and β'-subunits of RNA polymerase. Depending on the level of induction, the cells stop growing or grow at a rate close to that of the wild-type strain. This strategy has been successful in transforming growth-arrested bacteria into biofactories with a high production yield, releasing cellular resources from growth towards biosynthesis. However, high selection pressure is placed on a growth-arrested population, favoring mutations that allow cells to escape from growth control. Accordingly, we made the design of the growth switch more robust by building in genetic redundancy. More specifically, we added the rpoA gene, encoding for the α-subunit of RNA polymerase, under the control of a copy of the same inducible promoter used for expression control of ββ'. The improved growth switch is much more stable (escape frequency <10-9), while preserving the capacity to improve production yields. Moreover, after a long period of growth inhibition the population can be regenerated within a few generations. This opens up the possibility to alternate biomass accumulation and product synthesis over a longer period of time and is an additional step towards the dynamical control of bioproduction.

Clinically relevant mutations in regulatory regions of metabolic genes facilitate early adaptation to ciprofloxacin in Escherichia coli

The genomic landscape associated with early adaptation to ciprofloxacin is poorly understood. Although the interplay between core metabolism and antimicrobial resistance is being increasingly recognized, mutations in metabolic genes and their biological role remain elusive. Here, we exposed Escherichia coli to increasing gradients of ciprofloxacin with intermittent transfer-bottlenecking and identified mutations in three non-canonical targets linked to metabolism including a deletion (tRNA-ArgΔ414-bp) and point mutations in the regulatory regions of argI (ARG box) and narU. Our findings suggest that these mutations modulate arginine and carbohydrate metabolism, facilitate anaerobiosis and increased ATP production during ciprofloxacin stress. Furthermore, mutations in the regulatory regions of argI and narU were detected in over 70% of sequences from clinical E. coli isolates and were overrepresented among ciprofloxacin-resistant isolates. In sum, we have identified clinically relevant mutations in the regulatory regions of metabolic genes as a central theme that drives physiological changes necessary for adaptation to ciprofloxacin stress.

An additional proofreader contributes to DNA replication fidelity in mycobacteria

The removal of mis-incorporated nucleotides by proofreading activity ensures DNA replication fidelity. Whereas the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, it has been shown that proofreading in a majority of bacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase, despite the presence of a DnaQ homolog that is structurally and functionally distinct from E. coli DnaQ. However, the biological functions of this type of noncanonical DnaQ remain unclear. Here, we provide independent evidence that noncanonical DnaQ functions as an additional proofreader for mycobacteria. Using the mutation accumulation assay in combination with whole-genome sequencing, we showed that depletion of DnaQ in Mycolicibacterium smegmatis leads to an increased mutation rate, resulting in AT-biased mutagenesis and increased insertions/deletions in the homopolymer tract. Our results showed that mycobacterial DnaQ binds to the β clamp and functions synergistically with the PHP domain proofreader to correct replication errors. Furthermore, the loss of dnaQ results in replication fork dysfunction, leading to attenuated growth and increased mutagenesis on subinhibitory fluoroquinolones potentially due to increased vulnerability to fork collapse. By analyzing the sequence polymorphism of dnaQ in clinical isolates of Mycobacterium tuberculosis (Mtb), we demonstrated that a naturally evolved DnaQ variant prevalent in Mtb lineage 4.3 may enable hypermutability and is associated with drug resistance. These results establish a coproofreading model and suggest a division of labor between DnaQ and PHP domain proofreader. This study also provides real-world evidence that a mutator-driven evolutionary pathway may exist during the adaptation of Mtb.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

Variation in the Spectrum of New Mutations among Inbred Strains of Mice

The mouse serves as a mammalian model for understanding the nature of variation from new mutations, a question that has both evolutionary and medical significance. Previous studies suggest that the rate of single-nucleotide mutations (SNMs) in mice is ∼50% of that in humans. However, information largely comes from studies involving the C57BL/6 strain, and there is little information from other mouse strains. Here, we study the mutations that accumulated in 59 mouse lines derived from four inbred strains that are commonly used in genetics and clinical research (BALB/cAnNRj, C57BL/6JRj, C3H/HeNRj, and FVB/NRj), maintained for eight to nine generations by brother-sister mating. By analyzing Illumina whole-genome sequencing data, we estimate that the average rate of new SNMs in mice is ∼μ = 6.7 × 10-9. However, there is substantial variation in the spectrum of SNMs among strains, so the burden from new mutations also varies among strains. For example, the FVB strain has a spectrum that is markedly skewed toward C→A transversions and is likely to experience a higher deleterious load than other strains, due to an increased frequency of nonsense mutations in glutamic acid codons. Finally, we observe substantial variation in the rate of new SNMs among DNA sequence contexts, CpG sites, and their adjacent nucleotides playing an important role.

Measuring the burden of hundreds of BioBricks defines an evolutionary limit on constructability in synthetic biology

Engineered DNA will slow the growth of a host cell if it redirects limiting resources or otherwise interferes with homeostasis. Escape mutants that alleviate this burden can rapidly evolve and take over cell populations, making genetic engineering less reliable and predictable. Synthetic biologists often use genetic parts encoded on plasmids, but their burden is rarely characterized. We measured how 301 BioBrick plasmids affected Escherichia coli growth and found that 59 (19.6%) were burdensome, primarily because they depleted the limited gene expression resources of host cells. Overall, no BioBricks reduced the growth rate of E. coli by >45%, which agreed with a population genetic model that predicts such plasmids should be unclonable. We made this model available online for education ( https://barricklab.org/burden-model ) and added our burden measurements to the iGEM Registry. Our results establish a fundamental limit on what DNA constructs and genetic modifications can be successfully engineered into cells.

The DNA damage response of Escherichia coli, revisited: Differential gene expression after replication inhibition

In 1967, in this journal, Evelyn Witkin proposed the existence of a coordinated DNA damage response in Escherichia coli, which later came to be called the "SOS response." We revisited this response using the replication inhibitor azidothymidine (AZT) and RNA-Seq analysis and identified several features. We confirm the induction of classic Save our ship (SOS) loci and identify several genes, including many of the pyrimidine pathway, that have not been previously demonstrated to be DNA damage-inducible. Despite a strong dependence on LexA, these genes lack LexA boxes and their regulation by LexA is likely to be indirect via unknown factors. We show that the transcription factor "stringent starvation protein" SspA is as important as LexA in the regulation of AZT-induced genes and that the genes activated by SspA change dramatically after AZT exposure. Our experiments identify additional LexA-independent DNA damage inducible genes, including 22 small RNA genes, some of which appear to activated by SspA. Motility and chemotaxis genes are strongly down-regulated by AZT, possibly as a result of one of more of the small RNAs or other transcription factors such as AppY and GadE, whose expression is elevated by AZT. Genes controlling the iron siderophore, enterobactin, and iron homeostasis are also strongly induced, independent of LexA. We confirm that IraD antiadaptor protein is induced independent of LexA and that a second antiadaptor, IraM is likewise strongly AZT-inducible, independent of LexA, suggesting that RpoS stabilization via these antiadaptor proteins is an integral part of replication stress tolerance.

Escherichia coli DNA replication: the old model organism still holds many surprises

Research on Escherichia coli DNA replication paved the groundwork for many breakthrough discoveries with important implications for our understanding of human molecular biology, due to the high level of conservation of key molecular processes involved. To this day, it attracts a lot of attention, partially by virtue of being an important model organism, but also because the understanding of factors influencing replication fidelity might be important for studies on the emergence of antibiotic resistance. Importantly, the wide access to high-resolution single-molecule and live-cell imaging, whole genome sequencing, and cryo-electron microscopy techniques, which were greatly popularized in the last decade, allows us to revisit certain assumptions about the replisomes and offers very detailed insight into how they work. For many parts of the replisome, step-by-step mechanisms have been reconstituted, and some new players identified. This review summarizes the latest developments in the area, focusing on (a) the structure of the replisome and mechanisms of action of its components, (b) organization of replisome transactions and repair, (c) replisome dynamics, and (d) factors influencing the base and sugar fidelity of DNA synthesis.

Synthetic genomes unveil the effects of synonymous recoding

Engineering the genetic code of an organism provides the basis for (i) making any organism safely resistant to natural viruses and (ii) preventing genetic information flow into and out of genetically modified organisms while (iii) allowing the biosynthesis of genetically encoded unnatural polymers1-4. Achieving these three goals requires the reassignment of multiple of the 64 codons nature uses to encode proteins. However, synonymous codon replacement-recoding-is frequently lethal, and how recoding impacts fitness remains poorly explored. Here, we explore these effects using whole-genome synthesis, multiplexed directed evolution, and genome-transcriptome-translatome-proteome co-profiling on multiple recoded genomes. Using this information, we assemble a synthetic Escherichia coli genome in seven sections using only 57 codons to encode proteins. By discovering the rules responsible for the lethality of synonymous recoding and developing a data-driven multi-omics-based genome construction workflow that troubleshoots synthetic genomes, we overcome the lethal effects of 62,007 synonymous codon swaps and 11,108 additional genomic edits. We show that synonymous recoding induces transcriptional noise including new antisense RNAs, leading to drastic transcriptome and proteome perturbation. As the elimination of select codons from an organism's genetic code results in the widespread appearance of cryptic promoters, we show that synonymous codon choice may naturally evolve to minimize transcriptional noise. Our work provides the first genome-scale description of how synonymous codon changes influence organismal fitness and paves the way for the construction of functional genomes that provide genetic firewalls from natural ecosystems and safely produce biopolymers, drugs, and enzymes with an expanded chemistry.

Manufacturing DNA in E. coli yields higher-fidelity DNA than in vitro enzymatic synthesis

Biotechnologies such as gene therapy have brought DNA vectors to the forefront of pharmaceuticals. The quality of starting material plays a pivotal role in determining final product quality. Here, we examined the fidelity of DNA replication using enzymatic methods (in vitro) compared to plasmid DNA produced in vivo in E. coli. Next-generation sequencing approaches rely on in vitro polymerases, which have inherent limitations in sensitivity. To address this challenge, we introduce a novel assay based on loss-of-function (LOF) mutations in the conditionally toxic sacB gene. Our findings show that DNA production in E. coli results in significantly fewer LOF mutations (80- to 3,000-fold less) compared to enzymatic DNA replication methods such as polymerase chain reaction (PCR) and rolling circle amplification (RCA). These results suggest that using DNA produced by PCR or RCA may introduce a substantial number of mutation impurities, potentially affecting the quality and yield of final pharmaceutical products. Our study underscores that DNA synthesized in vitro has a significantly higher mutation rate than DNA produced traditionally in E. coli. Therefore, utilizing in vitro enzymatically produced DNA in biotechnology and biomanufacturing may entail considerable fidelity-related risks, while using DNA starting material derived from E. coli substantially mitigates this risk.

Correction of non-random mutational biases along a linear bacterial chromosome by the mismatch repair endonuclease NucS

The linear chromosome of Streptomyces exhibits a highly compartmentalized structure with a conserved central region flanked by variable arms. As double strand break (DSB) repair mechanisms play a crucial role in shaping the genome plasticity of Streptomyces, we investigated the role of EndoMS/NucS, a recently characterized endonuclease involved in a non-canonical mismatch repair (MMR) mechanism in archaea and actinobacteria, that singularly corrects mismatches by creating a DSB. We showed that Streptomyces mutants lacking NucS display a marked colonial phenotype and a drastic increase in spontaneous mutation rate. In vitro biochemical assays revealed that NucS cooperates with the replication clamp to efficiently cleave G/T, G/G and T/T mismatched DNA by producing DSBs. These findings are consistent with the transition-shifted mutational spectrum observed in the mutant strains and reveal that NucS-dependent MMR specific task is to eliminate G/T mismatches generated by the DNA polymerase during replication. Interestingly, our data unveil a crescent-shaped distribution of the transition frequency from the replication origin towards the chromosomal ends, shedding light on a possible link between NucS-mediated DSBs and Streptomyces genome evolution.

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

Characterization of transfusion-relevant bacteria reference strains in a lyophilized format

BACKGROUND AND OBJECTIVES

Blood safety measures used by blood establishments to increase blood component safety can be validated using Transfusion-Relevant Bacterial Reference Strains (TRBRS). Ultra-cold storage conditions and manual preparation of the current TRBRS may restrict their practical use. To address this issue, the ISBT Transfusion-Transmitted Infectious Diseases Working Party's Bacterial Subgroup organized an international study to validate TRBRS in a user-friendly, lyophilised format.

MATERIALS AND METHODS

Two bacterial strains Klebsiella pneumoniae PEI-B-P-08 and Staphylococcus aureus PEI-B-P-63 were manufactured as lyophilised material. The lyophilised bacteria were distributed to 11 different labs worldwide to assess the robustness for enumeration, identification and determination of growth kinetics in platelet concentrates (PCs).

RESULTS

Production of lyophilised TRBRS had no impact on the growth properties compared with the traditional format. The new format allows a direct low-quantity spiking of approximately 30 bacteria in PCs for transfusion-relevant experiments. In addition, the lyophilised bacteria exhibit long-term stability across a broad temperature range and can even be directly rehydrated in PCs without losing viability. Interlaboratory comparative study demonstrated the robustness of the new format as 100% of spiked PC exhibited growth.

CONCLUSION

Lyophilised TRBRS provide a user-friendly material for transfusion-related studies. TRBRS in the new format have improved features that may lead to a more frequent use in the quality control of transfusion-related safety measures in the future.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

Spontaneous mutations and mutational responses to penicillin treatment in the bacterial pathogen Streptococcus pneumoniae D39

Bacteria with functional DNA repair systems are expected to have low mutation rates due to strong natural selection for genomic stability. However, our study of the wild-type Streptococcus pneumoniae D39, a pathogen responsible for many common diseases, revealed a high spontaneous mutation rate of 0.02 per genome per cell division in mutation-accumulation (MA) lines. This rate is orders of magnitude higher than that of other non-mutator bacteria and is characterized by a high mutation bias in the A/T direction. The high mutation rate may have resulted from a reduction in the overall efficiency of selection, conferred by the tiny effective population size in nature. In line with this, S. pneumoniae D39 also exhibited the lowest DNA mismatch-repair (MMR) efficiency among bacteria. Treatment with the antibiotic penicillin did not elevate the mutation rate, as penicillin did not induce DNA damage and S. pneumoniae lacks a stress response pathway. Our findings suggested that the MA results are applicable to within-host scenarios and provide insights into pathogen evolution.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00220-6.

Origin, evolution, and maintenance of gene-strand bias in bacteria

Gene-strand bias is a characteristic feature of bacterial genome organization wherein genes are preferentially encoded on the leading strand of replication, promoting co-orientation of replication and transcription. This co-orientation bias has evolved to protect gene essentiality, expression, and genomic stability from the harmful effects of head-on replication-transcription collisions. However, the origin, variation, and maintenance of gene-strand bias remain elusive. Here, we reveal that the frequency of inversions that alter gene orientation exhibits large variation across bacterial populations and negatively correlates with gene-strand bias. The density, distance, and distribution of inverted repeats show a similar negative relationship with gene-strand bias explaining the heterogeneity in inversions. Importantly, these observations are broadly evident across the entire bacterial kingdom uncovering inversions and inverted repeats as primary factors underlying the variation in gene-strand bias and its maintenance. The distinct catalytic subunits of replicative DNA polymerase have co-evolved with gene-strand bias, suggesting a close link between replication and the origin of gene-strand bias. Congruently, inversion frequencies and inverted repeats vary among bacteria with different DNA polymerases. In summary, we propose that the nature of replication determines the fitness cost of replication-transcription collisions, establishing a selection gradient on gene-strand bias by fine-tuning DNA sequence repeats and, thereby, gene inversions.

Trade-offs, trade-ups, and high mutational parallelism underlie microbial adaptation during extreme cycles of feast and famine

Microbes are evolutionarily robust organisms capable of rapid adaptation to complex stress, which enables them to colonize harsh environments. In nature, microbes are regularly challenged by starvation, which is a particularly complex stress because resource limitation often co-occurs with changes in pH, osmolarity, and toxin accumulation created by metabolic waste. Often overlooked are the additional complications introduced by eventual resource replenishment, as successful microbes must withstand rapid environmental shifts before swiftly capitalizing on replenished resources to avoid invasion by competing species. To understand how microbes navigate trade-offs between growth and survival, ultimately adapting to thrive in environments with extreme fluctuations, we experimentally evolved 16 Escherichia coli populations for 900 days in repeated feast/famine conditions with cycles of 100-day starvation before resource replenishment. Using longitudinal population-genomic analysis, we found that evolution in response to extreme feast/famine is characterized by narrow adaptive trajectories with high mutational parallelism and notable mutational order. Genetic reconstructions reveal that early mutations result in trade-offs for biofilm and motility but trade-ups for growth and survival, as these mutations conferred positively correlated advantages during both short-term and long-term culture. Our results demonstrate how microbes can navigate the adaptive landscapes of regularly fluctuating conditions and ultimately follow mutational trajectories that confer benefits across diverse environments.

Measuring the burden of hundreds of BioBricks defines an evolutionary limit on constructability in synthetic biology

Engineered DNA will slow the growth of a host cell if it redirects limiting resources or otherwise interferes with homeostasis. Populations of engineered cells can rapidly become dominated by "escape mutants" that evolve to alleviate this burden by inactivating the intended function. Synthetic biologists working with bacteria rely on genetic parts and devices encoded on plasmids, but the burden of different engineered DNA sequences is rarely characterized. We measured how 301 BioBricks on high-copy plasmids affected the growth rate of Escherichia coli. Of these, 59 (19.6%) negatively impacted growth. The burden imposed by engineered DNA is commonly associated with diverting ribosomes or other gene expression factors away from producing endogenous genes that are essential for cellular replication. In line with this expectation, BioBricks exhibiting burden were more likely to contain highly active constitutive promoters and strong ribosome binding sites. By monitoring how much each BioBrick reduced expression of a chromosomal GFP reporter, we found that the burden of most, but not all, BioBricks could be wholly explained by diversion of gene expression resources. Overall, no BioBricks reduced the growth rate of E. coli by >45%, which agreed with a population genetic model that predicts such plasmids should be "unclonable" because escape mutants will take over during growth of a bacterial colony or small laboratory culture from a transformed cell. We made this model available as an interactive web tool for synthetic biology education and added our burden measurements to the iGEM Registry descriptions of each BioBrick.

Mutational Signatures in Wild Type Escherichia coli Strains Reveal Predominance of DNA Polymerase Errors

While mutational processes operating in the Escherichia coli genome have been revealed by multiple laboratory experiments, the contribution of these processes to accumulation of bacterial polymorphism and evolution in natural environments is unknown. To address this question, we reconstruct signatures of distinct mutational processes from experimental data on E. coli hypermutators, and ask how these processes contribute to differences between naturally occurring E. coli strains. We show that both mutations accumulated in the course of evolution of wild-type strains in nature and in the lab-grown nonmutator laboratory strains are explained predominantly by the low fidelity of DNA polymerases II and III. By contrast, contributions specific to disruption of DNA repair systems cannot be detected, suggesting that temporary accelerations of mutagenesis associated with such disruptions are unimportant for within-species evolution. These observations demonstrate that accumulation of diversity in bacterial strains in nature is predominantly associated with errors of DNA polymerases.

Environmental and genetic influence on the rate and spectrum of spontaneous mutations in Escherichia coli

Spontaneous mutations are the ultimate source of novel genetic variation on which evolution operates. Although mutation rate is often discussed as a single parameter in evolution, it comprises multiple distinct types of changes at the level of DNA. Moreover, the rates of these distinct changes can be independently influenced by genomic background and environmental conditions. Using fluctuation tests, we characterized the spectrum of spontaneous mutations in Escherichia coli grown in low and high glucose environments. These conditions are known to affect the rate of spontaneous mutation in wild-type MG1655, but not in a ΔluxS deletant strain - a gene with roles in both quorum sensing and the recycling of methylation products used in E. coli's DNA repair process. We find an increase in AT>GC transitions in the low glucose environment, suggesting that processes relating to the production or repair of this mutation could drive the response of overall mutation rate to glucose concentration. Interestingly, this increase in AT>GC transitions is maintained by the glucose non-responsive ΔluxS deletant. Instead, an elevated rate of GC>TA transversions, more common in a high glucose environment, leads to a net non-responsiveness of overall mutation rate for this strain. Our results show how relatively subtle changes, such as the concentration of a carbon substrate or loss of a regulatory gene, can substantially influence the amount and nature of genetic variation available to selection.

CRISPR-Cas9-assisted genome editing in E. coli elevates the frequency of unintended mutations

Cas-assisted lambda Red recombineering techniques have rapidly become a mainstay of bacterial genome editing. Such techniques have been used to construct both individual mutants and massive libraries to assess the effects of genomic changes. We have found that a commonly used Cas9-assisted editing method results in unintended mutations elsewhere in the genome in 26% of edited clones. The unintended mutations are frequently found over 200 kb from the intended edit site and even over 10 kb from potential off-target sites. We attribute the high frequency of unintended mutations to error-prone polymerases expressed in response to dsDNA breaks introduced at the edit site. Most unintended mutations occur in regulatory or coding regions and thus may have phenotypic effects. Our findings highlight the risks associated with genome editing techniques involving dsDNA breaks in E. coli and likely other bacteria and emphasize the importance of sequencing the genomes of edited cells to ensure the absence of unintended mutations.

Impressive pan-genomic diversity of E. coli from a wild animal community near urban development reflects human impacts

Human and domesticated animal waste infiltrates global freshwater, terrestrial, and marine environments, widely disseminating fecal microbes, antibiotics, and other chemical pollutants. Emerging evidence suggests that guts of wild animals are being invaded by our microbes, including Escherichia coli, which face anthropogenic selective pressures to gain antimicrobial resistance (AMR) and increase virulence. However, wild animal sources remain starkly under-represented among genomic sequence repositories. We sequenced whole genomes of 145 E. coli isolates from 55 wild and 13 domestic animal fecal samples, averaging 2 (ranging 1-7) isolates per sample, on a preserve imbedded in a human-dominated landscape in California Bay Area, USA, to assess AMR, virulence, and pan-genomic diversity. With single nucleotide polymorphism analyses we predict potential transmission routes. We illustrate the usefulness of E. coli to aid our understanding of and ability to surveil the emergence of zoonotic pathogens created by the mixing of human and wild bacteria in the environment.

Evolution of pH-sensitive transcription termination during adaptation to repeated long-term starvation

Fluctuating environments that consist of regular cycles of co-occurring stress are a common challenge faced by cellular populations. For a population to thrive in constantly changing conditions, an ability to coordinate a rapid cellular response is essential. Here, we identify a mutation conferring an arginine-to-histidine (Arg to His) substitution in the transcription terminator Rho. The rho R109H mutation frequently arose in E. coli populations experimentally evolved under repeated long-term starvation conditions, during which feast and famine result in drastic environmental pH fluctuations. Metagenomic sequencing revealed that populations containing the rho mutation also possess putative loss-of-function mutations in ydcI, which encodes a recently characterized transcription factor associated with pH homeostasis. Genetic reconstructions of these mutations show that the rho allele confers a plastic alkaline-induced reduction of Rho function that, when found in tandem with a ΔydcI allele, leads to intracellular alkalinization and genetic assimilation of Rho mutant function. We further identify Arg to His substitutions at analogous sites in rho alleles from species originating from fluctuating alkaline environments. Our results suggest that Arg to His substitutions in global regulators of gene expression can serve to rapidly coordinate complex responses through pH sensing and shed light on how cellular populations across the tree of life use environmental cues to coordinate rapid responses to complex, fluctuating environments.

Prolonged carriage of OXA-244-carbapenemase-producing Escherichia coli complicates epidemiological investigations

The rapid increase of OXA-244-producing Escherichia coli, predominantly driven by genetically clustered isolates of sequence type (ST)38, has been observed in at least nine European countries, including Germany. However, the reasons for the spread of OXA-244-producing E. coli remain unclear. Here, we aim to evaluate the possibility of prolonged carriage. We identified a total of six different patients with repeated detection of OXA-244-producing E. coli isolates, which were subjected to both short and long-read whole-genome sequencing (WGS). Besides allelic differences using core genome multilocus sequence typing (cgMLST) analyses, we obtained numbers of single-nucleotide polymorphisms (SNPs) to calculate individual base-pair substitution (BPS) rates. To assess possible re-exposure and risk factors for prolonged carriage, case interviews were conducted. The time between detections ranged from eleven months to more than three years. Initial isolates originated in three+ out of six cases from clinical samples, whereas remaining samples were from screening, mostly in the inpatient setting. As expected, cgMLST analyses showed low numbers of allelic differences between isolates of each case ranging from 1 to 4, whereas numbers of SNPs were between 2 and 99 (mean = 36), thus clearly highlighting the discrepancy between these different bacterial typing approaches. For five out of six cases, observed BPS rates suggest that patients can be colonized with OXA-244-producing E. coli, including ST38 cluster isolates, for extensively long times. Thus, we may have previously missed the epidemiological link between cases because exposure to OXA-244-producing E. coli could have occurred in a time frame, which has not been evaluated in previous investigations. Our results may help to guide future epidemiological investigations as well as to support the interpretation of genetic diversity of OXA-244-producing E. coli, particularly among ST38 cluster isolates.

Genetic mutation in Escherichia coli genome during adaptation to the murine intestine is optimized for the host diet

Mammalian gut microbes colonize the intestinal tract of their host and adapt to establish a microbial ecosystem. The host diet changes the nutrient profile of the intestine and has a high impact on microbiota composition. Genetic mutations in Escherichia coli, a prevalent species in the human gut, allow for adaptation to the mammalian intestine, as reported in previous studies. However, the extent of colonization fitness in the intestine elevated by genetic mutation and the effects of diet change on these mutations in E. coli are still poorly known. Here, we show that notable mutations in sugar metabolism-related genes (gatC, araC, and malI) were detected in the E. coli K-12 genome just 2 weeks after colonization in the germ-free mouse intestine. In addition to elevated fitness by deletion of gatC, as previously reported, deletion of araC and malI also elevated E. coli fitness in the murine intestine in a host diet-dependent manner. In vitro cultures of medium containing nutrients abundant in the intestine (e.g., galactose, N-acetylglucosamine, and asparagine) also showed increased E. coli fitness after deletion of the genes-of-interest associated with their metabolism. Furthermore, the host diet was found to influence the developmental trajectory of gene mutations in E. coli. Taken together, we suggest that genetic mutations in E. coli are selected in response to the intestinal environment, which facilitates efficient utilization of nutrients abundant in the intestine under laboratory conditions. Our study offers some insight into the possible adaptation mechanisms of gut microbes.IMPORTANCEThe gut microbiota is closely associated with human health and is greatly impacted by the host diet. Bacteria such as Escherichia coli live in the gut all throughout the life of a human host and adapt to the intestinal environment. Adaptive mutations in E. coli are reported to enhance fitness in the mammalian intestine, but to what extent is still poorly known. It is also unknown whether the host diet affects what genes are mutated and to what extent fitness is affected. This study suggests that genetic mutations in the E. coli K-12 strain are selected in response to the intestinal environment and facilitate efficient utilization of abundant nutrients in the germ-free mouse intestine. Our study provides a better understanding of these intestinal adaptation mechanisms of gut microbes.

Glycolate as alternative carbon source for Escherichia coli

The physiology of different Escherichia coli stains was analyzed for growth with glycolate as a potentially promising sustainable sole source of carbon and energy. Different E. coli strains showed large differences regarding lag phases after provision of glycolate. Whereas E. coli W showed fast adaptation, E. coli BW25113, JM101, and BL21 (DE3) needed extensive time for adaption (up to 30 generations) until the attainable µmax was reached, which, at 30 °C, amounted to 0.20-0.25 h-1 for all strains. The overexpression of genes encoding glycolate degradation did neither overcome the need for adaptation of E. coli BL21 (DE3) nor improve growth of E. coli W. Rather, high level expression of proteins involved in uptake and initial degradation steps had an adverse effect on growth. Overall, the results show a promising capacity of E. coli strains for growth on glycolate.

Sub-inhibitory antibiotic treatment selects for enhanced metabolic efficiency

Bacterial growth and metabolic rates are often closely related. However, under antibiotic selection, a paradox in this relationship arises: antibiotic efficacy decreases when bacteria are metabolically dormant, yet antibiotics select for resistant cells that grow fastest during treatment. That is, antibiotic selection counterintuitively favors bacteria with fast growth but slow metabolism. Despite this apparent contradiction, antibiotic resistant cells have historically been characterized primarily in the context of growth, whereas the extent of analogous changes in metabolism is comparatively unknown. Here, we observed that previously evolved antibiotic-resistant strains exhibited a unique relationship between growth and metabolism whereby nutrient utilization became more efficient, regardless of the growth rate. To better understand this unexpected phenomenon, we used a simplified model to simulate bacterial populations adapting to sub-inhibitory antibiotic selection through successive bottlenecking events. Simulations predicted that sub-inhibitory bactericidal antibiotic concentrations could select for enhanced metabolic efficiency, defined based on nutrient utilization: drug-adapted cells are able to achieve the same biomass while utilizing less substrate, even in the absence of treatment. Moreover, simulations predicted that restoring metabolic efficiency would re-sensitize resistant bacteria exhibiting metabolic-dependent resistance; we confirmed this result using adaptive laboratory evolutions of Escherichia coli under carbenicillin treatment. Overall, these results indicate that metabolic efficiency is under direct selective pressure during antibiotic treatment and that differences in evolutionary context may determine both the efficacy of different antibiotics and corresponding re-sensitization approaches.IMPORTANCEThe sustained emergence of antibiotic-resistant pathogens combined with the stalled drug discovery pipelines highlights the critical need to better understand the underlying evolution mechanisms of antibiotic resistance. To this end, bacterial growth and metabolic rates are often closely related, and resistant cells have historically been characterized exclusively in the context of growth. However, under antibiotic selection, antibiotics counterintuitively favor cells with fast growth, and slow metabolism. Through an integrated approach of mathematical modeling and experiments, this study thereby addresses the significant knowledge gap of whether antibiotic selection drives changes in metabolism that complement, and/or act independently, of antibiotic resistance phenotypes.

Bacterial Mutation During Seasonal Epidemics

Rapidly evolving bacterial pathogens pose a unique challenge for long-term plant disease management. In this study, we investigated the types and rate of mutations in bacterial populations during seasonal disease epidemics. Two phylogenetically distinct strains of the bacterial spot pathogen, Xanthomonas perforans, were marked, released in tomato fields, and recaptured at several time points during the growing season. Genomic variations in recaptured isolates were identified by comparative analysis of their whole-genome sequences. In total, 180 unique variations (116 substitutions, 57 insertions/deletions, and 7 structural variations) were identified from 300 genomes, resulting in the overall host-associated mutation rate of ∼0.3 to 0.9/genome/week. This result serves as a benchmark for bacterial mutation during epidemics in similar pathosystems. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Establishing a synthetic orthogonal replication system enables accelerated evolution in E. coli

The evolution of new function in living organisms is slow and fundamentally limited by their critical mutation rate. Here, we established a stable orthogonal replication system in Escherichia coli. The orthogonal replicon can carry diverse cargos of at least 16.5 kilobases and is not copied by host polymerases but is selectively copied by an orthogonal DNA polymerase (O-DNAP), which does not copy the genome. We designed mutant O-DNAPs that selectively increase the mutation rate of the orthogonal replicon by two to four orders of magnitude. We demonstrate the utility of our system for accelerated continuous evolution by evolving a 150-fold increase in resistance to tigecycline in 12 days. And, starting from a GFP variant, we evolved a 1000-fold increase in cellular fluorescence in 5 days.

Ultrahigh-throughput screening-assisted in vivo directed evolution for enzyme engineering

BACKGROUND

Classical directed evolution is a powerful approach for engineering biomolecules with improved or novel functions. However, it traditionally relies on labour- and time-intensive iterative cycles, due in part to the need for multiple molecular biology steps, including DNA transformation, and limited screening throughput.

RESULTS

In this study, we present an ultrahigh throughput in vivo continuous directed evolution system with thermosensitive inducible tunability, which is based on error-prone DNA polymerase expression modulated by engineered thermal-responsive repressor cI857, and genomic MutS mutant with temperature-sensitive defect for fixation of mutations in Escherichia coli. We demonstrated the success of the in vivo evolution platform with β-lactamase as a model, with an approximately 600-fold increase in the targeted mutation rate. Furthermore, the platform was combined with ultrahigh-throughput screening methods and employed to evolve α-amylase and the resveratrol biosynthetic pathway. After iterative rounds of enrichment, a mutant with a 48.3% improvement in α-amylase activity was identified via microfluidic droplet screening. In addition, when coupled with an in vivo biosensor in the resveratrol biosynthetic pathway, a variant with 1.7-fold higher resveratrol production was selected by fluorescence-activated cell sorting.

CONCLUSIONS

In this study, thermal-responsive targeted mutagenesis coupled with ultrahigh-throughput screening was developed for the rapid evolution of enzymes and biosynthetic pathways.

First single-strain enrichments of Electrothrix cable bacteria, description of E. aestuarii sp. nov. and E. rattekaaiensis sp. nov., and proposal of a cable bacteria taxonomy following the rules of the SeqCode

Cable bacteria are electrically conductive, filamentous Desulfobulbaceae, which are morphologically, functionally, and phylogenetically distinct from the other members of this family. Cable bacteria have not been obtained in pure culture and were therefore previously described as candidate genera, Candidatus Electrothrix and Ca. Electronema; a representative of the latter is available as single-strain sediment enrichment. Here we present an improved workflow to obtain the first single-strain enrichments of Ca. Electrothrix and report their metagenome-assembled genomes (MAGs) and morphology. Based on these results and on previously published high-quality MAGs and morphological data of cable bacteria from both candidate genera, we propose to adopt the genus names Electrothrix and Electronema following the rules of the Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode), with Electrothrix communis RBTS and Electronema aureum GSTS, respectively, as the nomenclatural types of the genera. Furthermore, based on average nucleotide identity (ANI) values < 95 % with any described species, we propose two of our three single-strain enrichment cultures as novel species of the genus Electrothrix, with the names E. aestuarii sp. nov. and E. rattekaaiensis sp. nov., according to the SeqCode.

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

Trade-offs constrain adaptive pathways to the type VI secretion system survival

The Type VI Secretion System (T6SS) is a nano-harpoon used by many bacteria to inject toxins into neighboring cells. While much is understood about mechanisms of T6SS-mediated toxicity, less is known about the ways that competitors can defend themselves against this attack, especially in the absence of their own T6SS. Here we subjected eight replicate populations of Escherichia coli to T6SS attack by Vibrio cholerae. Over ∼500 generations of competition, isolates of the E. coli populations evolved to survive T6SS attack an average of 27-fold better, through two convergently evolved pathways: apaH was mutated in six of the eight replicate populations, while the other two populations each had mutations in both yejM and yjeP. However, the mutations we identified are pleiotropic, reducing cellular growth rates, and increasing susceptibility to antibiotics and elevated pH. These trade-offs help us understand how the T6SS shapes the evolution of bacterial interactions.

Competition and evolutionary selection among core regulatory motifs in gene expression control

Gene products that are beneficial in one environment may become burdensome in another, prompting the emergence of diverse regulatory schemes that carry their own bioenergetic cost. By ensuring that regulators are only expressed when needed, we demonstrate that autoregulation generally offers an advantage in an environment combining mutation and time-varying selection. Whether positive or negative feedback emerges as dominant depends primarily on the demand for the target gene product, typically to ensure that the detrimental impact of inevitable mutations is minimized. While self-repression of the regulator curbs the spread of these loss-of-function mutations, self-activation instead facilitates their propagation. By analyzing the transcription network of multiple model organisms, we reveal that reduced bioenergetic cost may contribute to the preferential selection of autoregulation among transcription factors. Our results not only uncover how seemingly equivalent regulatory motifs have fundamentally different impact on population structure, growth dynamics, and evolutionary outcomes, but they can also be leveraged to promote the design of evolutionarily robust synthetic gene circuits.

Ultra-high static magnetic field induces a change in the spectrum but not frequency of DNA spontaneous mutations in Arabidopsis thaliana

Biological effects of magnetic fields have been extensively studied in plants, microorganisms and animals, and applications of magnetic fields in regulation of plant growth and phytoprotection is a promising field in sustainable agriculture. However, the effect of magnetic fields especially ultra-high static magnetic field (UHSMF) on genomic stability is largely unclear. Here, we investigated the mutagenicity of 24.5, 30.5 and 33.0 T UHSMFs with the gradient of 150, 95 and 0 T/m, respectively, via whole genome sequencing. Our results showed that 1 h exposure of Arabidopsis dried seeds to UHSMFs has no significant effect on the average rate of DNA mutations including single nucleotide variations and InDels (insertions and deletions) in comparison with the control, but 33.0 T and 24.5 T treatments lead to a significant change in the rate of nucleotide transitions and InDels longer than 3 bp, respectively, suggesting that both strength and gradient of UHSMF impact molecular spectrum of DNA mutations. We also found that the decreased transition rate in UHSMF groups is correlated with the upstream flanking sequences of G and C mutation sites. Furthermore, the germination rate of seeds exposed to 24.5 T SMF with -150 T/m gradient showed a significant decrease at 24 hours after sowing. Overall, our data lay a basis for precisely assessing the potential risk of UHSMF on DNA stability, and for elucidating molecular mechanism underlying gradient SMF-regulated biological processes in the future.

A rugged yet easily navigable fitness landscape

Fitness landscape theory predicts that rugged landscapes with multiple peaks impair Darwinian evolution, but experimental evidence is limited. In this study, we used genome editing to map the fitness of >260,000 genotypes of the key metabolic enzyme dihydrofolate reductase in the presence of the antibiotic trimethoprim, which targets this enzyme. The resulting landscape is highly rugged and harbors 514 fitness peaks. However, its highest peaks are accessible to evolving populations via abundant fitness-increasing paths. Different peaks share large basins of attraction that render the outcome of adaptive evolution highly contingent on chance events. Our work shows that ruggedness need not be an obstacle to Darwinian evolution but can reduce its predictability. If true in general, the complexity of optimization problems on realistic landscapes may require reappraisal.

Gene birth in a model of non-genic adaptation

BACKGROUND

Over evolutionary timescales, genomic loci can switch between functional and non-functional states through processes such as pseudogenization and de novo gene birth. Particularly, de novo gene birth is a widespread process, and many examples continue to be discovered across diverse evolutionary lineages. However, the general mechanisms that lead to functionalization are poorly understood, and estimated rates of de novo gene birth remain contentious. Here, we address this problem within a model that takes into account mutations and structural variation, allowing us to estimate the likelihood of emergence of new functions at non-functional loci.

RESULTS

Assuming biologically reasonable mutation rates and mutational effects, we find that functionalization of non-genic loci requires the realization of strict conditions. This is in line with the observation that most de novo genes are localized to the vicinity of established genes. Our model also provides an explanation for the empirical observation that emerging proto-genes are often lost despite showing signs of adaptation.

CONCLUSIONS

Our work elucidates the properties of non-genic loci that make them fertile for adaptation, and our results offer mechanistic insights into the process of de novo gene birth.

The geologic history of primary productivity

The rate of primary productivity is a keystone variable in driving biogeochemical cycles today and has been throughout Earth's past.1 For example, it plays a critical role in determining nutrient stoichiometry in the oceans,2 the amount of global biomass,3 and the composition of Earth's atmosphere.4 Modern estimates suggest that terrestrial and marine realms contribute near-equal amounts to global gross primary productivity (GPP).5 However, this productivity balance has shifted significantly in both recent times6 and through deep time.7,8 Combining the marine and terrestrial components, modern GPP fixes ≈250 billion tonnes of carbon per year (Gt C year-1).5,9,10,11 A grand challenge in the study of the history of life on Earth has been to constrain the trajectory that connects present-day productivity to the origin of life. Here, we address this gap by piecing together estimates of primary productivity from the origin of life to the present day. We estimate that ∼1011-1012 Gt C has cumulatively been fixed through GPP (≈100 times greater than Earth's entire carbon stock). We further estimate that 1039-1040 cells have occupied the Earth to date, that more autotrophs than heterotrophs have ever existed, and that cyanobacteria likely account for a larger proportion than any other group in terms of the number of cells. We discuss implications for evolutionary trajectories and highlight the early Proterozoic, which encompasses the Great Oxidation Event (GOE), as the time where most uncertainty exists regarding the quantitative census presented here.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

Impact of transient acquired hypermutability on the inter- and intra-species competitiveness of Pseudomonas aeruginosa

Once acquired, hypermutation is unrelenting, and in the long-term, leads to impaired fitness due to its cumulative impact on the genome. This raises the question of why hypermutators arise so frequently in microbial ecosystems. In this work, we explore this problem by examining how the transient acquisition of hypermutability affects inter- and intra-species competitiveness, and the response to environmental insults such as antibiotic challenge. We do this by engineering Pseudomonas aeruginosa to allow the expression of an important mismatch repair gene, mutS, to be experimentally controlled over a wide dynamic range. We show that high levels of mutS expression induce genomic stasis (hypomutation), whereas lower levels of induction lead to progressively higher rates of mutation. Whole-genome sequence analyses confirmed that the mutational spectrum of the inducible hypermutator is similar to the distinctive profile associated with mutS mutants obtained from the airways of people with cystic fibrosis (CF). The acquisition of hypermutability conferred a distinct temporal fitness advantage over the wild-type P. aeruginosa progenitor strain, in both the presence and the absence of an antibiotic selection pressure. However, over a similar time-scale, acquisition of hypermutability had little impact on the population dynamics of P. aeruginosa when grown in the presence of a competing species (Staphylococcus aureus). These data indicate that in the short term, acquired hypermutability primarily confers a competitive intra-species fitness advantage.

Comparative genomics hints at dispensability of multiple essential genes in two Escherichia coli L-form strains

Despite the critical role of bacterial cell walls in maintaining cell shapes, certain environmental stressors can induce the transition of many bacterial species into a wall-deficient state called L-form. Long-term induced Escherichia coli L-forms lose their rod shape and usually hold significant mutations that affect cell division and growth. Besides this, the genetic background of L-form bacteria is still poorly understood. In the present study, the genomes of two stable L-form strains of E. coli (NC-7 and LWF+) were sequenced and their gene mutation status was determined and compared with their parental strains. Comparative genomic analysis between two L-forms reveals both unique adaptions and common mutated genes, many of which belong to essential gene categories not involved in cell wall biosynthesis, indicating that L-form genetic adaptation impacts crucial metabolic pathways. Missense variants from L-forms and Lenski's long-term evolution experiment (LTEE) were analyzed in parallel using an optimized DeepSequence pipeline to investigate predicted mutation effects (α) on protein functions. We report that the two L-form strains analyzed display a frequency of 6-10% (0% for LTEE) in mutated essential genes where the missense variants have substantial impact on protein functions (α<0.5). This indicates the emergence of different survival strategies in L-forms through changes in essential genes during adaptions to cell wall deficiency. Collectively, our results shed light on the detailed genetic background of two E. coli L-forms and pave the way for further investigations of the gene functions in L-form bacterial models.

Differential adenine methylation analysis reveals increased variability in 6mA in the absence of methyl-directed mismatch repair

IMPORTANCE

Methylation greatly influences the bacterial genome by guiding DNA repair and regulating pathogenic and stress-response phenotypes. But, the rate of epigenetic changes and their consequences on molecular phenotypes are underexplored. Through a detailed characterization of genome-wide adenine methylation in a commonly used laboratory strain of Escherichia coli, we reveal that mismatch repair deficient populations experience an increase in epimutations resulting in a genome-wide reduction of 6mA methylation in a manner consistent with genetic drift. Our findings highlight how methylation patterns evolve and the constraints on epigenetic evolution due to post-replicative DNA repair, contributing to a deeper understanding of bacterial genome evolution and how epimutations may introduce semi-permanent variation that can influence adaptation.

Mycobacterial DnaQ is an Alternative Proofreader Ensuring DNA Replication Fidelity

Remove of mis-incorporated nucleotides ensures replicative fidelity. Although the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, proofreading in mycobacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase despite the presence of an alternative DnaQ homolog. Here, we show that depletion of DnaQ in Mycolicibacterium smegmatis results in increased mutation rate, leading to AT-biased mutagenesis and elevated insertions/deletions in homopolymer tract. We demonstrated that mycobacterial DnaQ binds to the b-clamp and functions synergistically with the PHP domain to correct replication errors. Further, we found that the mycobacterial DnaQ sustains replicative fidelity upon chromosome topological stress. Intriguingly, we showed that a naturally evolved DnaQ variant prevalent in clinical Mycobacterium tuberculosis isolates enables hypermutability and is associated with extensive drug resistance. These results collectively establish that the alternative DnaQ functions in proofreading, and thus reveal that mycobacteria deploy two proofreaders to maintain replicative fidelity.

Shifts in Mutation Bias Promote Mutators by Altering the Distribution of Fitness Effects

AbstractRecent experimental evidence demonstrates that shifts in mutational biases-for example, increases in transversion frequency-can change the distribution of fitness effects of mutations (DFE). In particular, reducing or reversing a prevailing bias can increase the probability that a de novo mutation is beneficial. It has also been shown that mutator bacteria are more likely to emerge if the beneficial mutations they generate have a larger effect size than observed in the wild type. Here, we connect these two results, demonstrating that mutator strains that reduce or reverse a prevailing bias have a positively shifted DFE, which in turn can dramatically increase their emergence probability. Since changes in mutation rate and bias are often coupled through the gain and loss of DNA repair enzymes, our results predict that the invasion of mutator strains will be facilitated by shifts in mutation bias that offer improved access to previously undersampled beneficial mutations.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

Model of neural development by differentiating human induced pluripotent stem cells into neural progenitor cells to study the neurodevelopmental toxicity of lead

Lead (Pb) exposure causes immeasurable damage to multiple human systems, particularly the central nervous system (CNS). In this study, human induced pluripotent stem cells (hiPSCs) were differentiated into neural progenitor cells (NPCs) to investigate the neurotoxic effects of Pb. The hiPSCs were treated with 0, 0.5, 1.0, 2.5, 5.0 and 10.0 μmol/L Pb for 7 days, whereas embryoid bodies (EBs) and NPCs were treated with 0, 0.1, 0.5, and 1.0 μmol/L Pb for 7 days. Pb exposure disrupted the cell cycle and caused apoptosis in hiPSCs, EBs, and NPCs. Besides, Pb inhibited the differentiation of NPCs and EBs. Whole exome sequencing revealed 2509, 2413, and 1984 single nucleotide variants (SNVs) caused by Pb in hiPSCs, EBs, and NPCs, respectively. The common mutation sites in the exon region were mostly nonsynonymous mutations. We identified 18, 19, and 18 common deleterious mutations in hiPSCs, EBs, and NPCs, respectively. Additionally, Online Mendelian Inheritance in Man database analysis revealed 30, 20, and 13 genes related to CNS disorders in hiPSCs, EBs, and NPCs, respectively. Our findings suggest that this in vitro model may supplement animal models and be applied to the study of neurodevelopmental toxicity in the future.

Phenotype switching of the mutation rate facilitates adaptive evolution

The mutation rate plays an important role in adaptive evolution. It can be modified by mutator and anti-mutator alleles. Recent empirical evidence hints that the mutation rate may vary among genetically identical individuals: evidence from bacteria suggests that the mutation rate can be affected by expression noise of a DNA repair protein and potentially also by translation errors in various proteins. Importantly, this non-genetic variation may be heritable via a transgenerational epigenetic mode of inheritance, giving rise to a mutator phenotype that is independent from mutator alleles. Here, we investigate mathematically how the rate of adaptive evolution is affected by the rate of mutation rate phenotype switching. We model an asexual population with two mutation rate phenotypes, non-mutator and mutator. An offspring may switch from its parental phenotype to the other phenotype. We find that switching rates that correspond to so-far empirically described non-genetic systems of inheritance of the mutation rate lead to higher rates of adaptation on both artificial and natural fitness landscapes. These switching rates can maintain within the same individuals both a mutator phenotype and intermediary mutations, a combination that facilitates adaptation. Moreover, non-genetic inheritance increases the proportion of mutators in the population, which in turn increases the probability of hitchhiking of the mutator phenotype with adaptive mutations. This in turns facilitates the acquisition of additional adaptive mutations. Our results rationalize recently observed noise in the expression of proteins that affect the mutation rate and suggest that non-genetic inheritance of this phenotype may facilitate evolutionary adaptive processes.

Retron-mediated multiplex genome editing and continuous evolution in Escherichia coli

While there are several genome editing techniques available, few are suitable for dynamic and simultaneous mutagenesis of arbitrary targeted sequences in prokaryotes. Here, to address these limitations, we present a versatile and multiplex retron-mediated genome editing system (REGES). First, through systematic optimization of REGES, we achieve efficiency of ∼100%, 85 ± 3%, 69 ± 14% and 25 ± 14% for single-, double-, triple- and quadruple-locus genome editing, respectively. In addition, we employ REGES to generate pooled and barcoded variant libraries with degenerate RBS sequences to fine-tune the expression level of endogenous and exogenous genes, such as transcriptional factors to improve ethanol tolerance and biotin biosynthesis. Finally, we demonstrate REGES-mediated continuous in vivo protein evolution, by combining retron, polymerase-mediated base editing and error-prone transcription. By these case studies, we demonstrate REGES as a powerful multiplex genome editing and continuous evolution tool with broad applications in synthetic biology and metabolic engineering.

Selective Pressure by Rifampicin Modulates Mutation Rates and Evolutionary Trajectories of Mycobacterial Genomes

Resistance to the frontline antibiotic rifampicin constitutes a challenge to the treatment and control of tuberculosis. Here, we analyzed the mutational landscape of Mycobacterium smegmatis during long-term evolution with increasing concentrations of rifampicin, using a mutation accumulation assay combined with whole-genome sequencing. Antibiotic treatment enhanced the acquisition of mutations, doubling the genome-wide mutation rate of the wild-type cells. While antibiotic exposure led to extinction of almost all wild-type lines, the hypermutable phenotype of the ΔnucS mutant strain (noncanonical mismatch repair deficient) provided an efficient response to the antibiotic, leading to high rates of survival. This adaptative advantage resulted in the emergence of higher levels of rifampicin resistance, an accelerated acquisition of drug resistance mutations in rpoB (β RNA polymerase), and a wider diversity of evolutionary pathways that led to drug resistance. Finally, this approach revealed a subset of adaptive genes under positive selection with rifampicin that could be associated with the development of antibiotic resistance. IMPORTANCE Rifampicin is the most important first-line antibiotic against mycobacterial infections, including tuberculosis, one of the top causes of death worldwide. Acquisition of rifampicin resistance constitutes a major global public health problem that makes the control of the disease challenging. Here, we performed an experimental evolution assay under antibiotic selection to analyze the response and adaptation of mycobacteria, leading to the acquisition of rifampicin resistance. This approach explored the total number of mutations that arose in the mycobacterial genomes under long-term rifampicin exposure, using whole-genome sequencing. Our results revealed the effect of rifampicin at a genomic level, identifying different mechanisms and multiple pathways leading to rifampicin resistance in mycobacteria. Moreover, this study detected that an increase in the rate of mutations led to enhanced levels of drug resistance and survival. In summary, all of these results could be useful to understand and prevent the emergence of drug-resistant isolates in mycobacterial infections.