Corrigendum: The selective advantage of the lac operon for Escherichia coli is conditional on diet and microbiota composition.

[This corrects the article DOI: 10.3389/fmicb.2021.709259.].

Shared Evolutionary Path in Social Microbiomes.

Social networks can influence the ecology of gut bacteria, shaping the species composition of the gut microbiome in humans and other animals. Gut commensals evolve and can adapt at a rapid pace when colonizing healthy hosts. Here, we aimed at assessing the impact of host-to-host bacterial transmission on Escherichia coli evolution in the mammalian gut. Using an in vivo experimental evolution approach in mice, we found a transmission rate of 7% (±3% 2× standard error [2SE]) of E. coli cells per day between hosts inhabiting the same household. Consistent with the predictions of a simple population genetics model of mutation-selection-migration, the level of shared events resulting from within host evolution is greatly enhanced in cohoused mice, showing that hosts undergoing the same diet and habit are not only expected to have similar microbiome species compositions but also similar microbiome evolutionary dynamics. Furthermore, we estimated the rate of mutation accumulation of E. coli to be 3.0 × 10-3 (±0.8 × 10-3 2SE) mutations/genome/generation, irrespective of the social context of the regime. Our results reveal the impact of bacterial migration across hosts in shaping the adaptive evolution of new strains colonizing gut microbiomes.

Ecotype formation and prophage domestication during gut bacterial evolution.

How much bacterial evolution occurs in our intestines and which factors control it are currently burning questions. The formation of new ecotypes, some of which capable of coexisting for long periods of time, is highly likely in our guts. Horizontal gene transfer driven by temperate phages that can perform lysogeny is also widespread in mammalian intestines. Yet, the roles of mutation and especially lysogeny as key drivers of gut bacterial adaptation remain poorly understood. The mammalian gut contains hundreds of bacterial species, each with many strains and ecotypes, whose abundance varies along the lifetime of a host. A continuous high input of mutations and horizontal gene transfer events mediated by temperate phages drives that diversity. Future experiments to study the interaction between mutations that cause adaptation in microbiomes and lysogenic events with different costs and benefits will be key to understand the dynamic microbiomes of mammals. Also see the video abstract here: https://youtu.be/Zjqsiyb5Pk0.

Intrahost evolution of the gut microbiota.

A massive number of microorganisms, belonging to different species, continuously divide inside the guts of animals and humans. The large size of these communities and their rapid division times imply that we should be able to watch microbial evolution in the gut in real time, in a similar manner to what has been done in vitro. Here, we review recent findings on how natural selection shapes intrahost evolution (also known as within-host evolution), with a focus on the intestines of mice and humans. The microbiota of a healthy host is not as static as initially thought from the information measured at only one genomic marker. Rather, the genomes of each gut-colonizing species can be highly dynamic, and such dynamism seems to be related to the microbiota species diversity. Genetic and bioinformatic tools, and analysis of time series data, allow quantification of the selection strength on emerging mutations and horizontal transfer events in gut ecosystems. The drivers and functional consequences of gut evolution can now begin to be grasped. The rules of this intrahost microbiota evolution, and how they depend on the biology of each species, need to be understood for more effective development of microbiota therapies to help maintain or restore host health.

Horizontal gene transfer among host-associated microbes.

Horizontal gene transfer is an important evolutionary force, facilitating bacterial diversity. It is thought to be pervasive in host-associated microbiomes, where bacterial densities are high and mobile elements are frequent. These genetic exchanges are also key for the rapid dissemination of antibiotic resistance. Here, we review recent studies that have greatly extended our knowledge of the mechanisms underlying horizontal gene transfer, the ecological complexities of a network of interactions involving bacteria and their mobile elements, and the effect of host physiology on the rates of genetic exchanges. Furthermore, we discuss other, fundamental challenges in detecting and quantifying genetic exchanges in vivo, and how studies have contributed to start overcoming these challenges. We highlight the importance of integrating novel computational approaches and theoretical models with experimental methods where multiple strains and transfer elements are studied, both in vivo and in controlled conditions that mimic the intricacies of host-associated environments.

Epistasis decreases with increasing antibiotic pressure but not temperature.

Predicting mutational effects is essential for the control of antibiotic resistance (ABR). Predictions are difficult when there are strong genotype-by-environment (G × E), gene-by-gene (G × G or epistatic) or gene-by-gene-by-environment (G × G × E) interactions. We quantified G × G × E effects in Escherichia coli across environmental gradients. We created intergenic fitness landscapes using gene knock-outs and single-nucleotide ABR mutations previously identified to vary in the extent of G × E effects in our environments of interest. Then, we measured competitive fitness across a complete combinatorial set of temperature and antibiotic dosage gradients. In this way, we assessed the predictability of 15 fitness landscapes across 12 different but related environments. We found G × G interactions and rugged fitness landscapes in the absence of antibiotic, but as antibiotic concentration increased, the fitness effects of ABR genotypes quickly overshadowed those of gene knock-outs, and the landscapes became smoother. Our work reiterates that some single mutants, like those conferring resistance or susceptibility to antibiotics, have consistent effects across genetic backgrounds in stressful environments. Thus, although epistasis may reduce the predictability of evolution in benign environments, evolution may be more predictable in adverse environments. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Fluctuating selection on bacterial iron regulation in the mammalian gut.

Iron is critical in host-microbe interactions, and its availability is tightly regulated in the mammalian gut. Antibiotics and inflammation can perturb iron availability in the gut, which could alter host-microbe interactions. Here, we show that an adaptive allele of iscR, a major regulator of iron homeostasis of Escherichia coli, is under fluctuating selection in the mouse gut. In vivo competitions in immune-competent, immune-compromised, and germ-free mice reveal that the selective pressure on an iscR mutant E. coli is modulated by the presence of antibiotics, the microbiota, and the immune system. In vitro assays show that iron availability is an important mediator of the iscR allele fitness benefits or costs. We identify Lipocalin-2, a host's immune protein that prevents bacterial iron acquisition, as a major host mechanism underlying fluctuating selection of iscR. Our results provide a remarkable example of strong fluctuating selection acting on bacterial iron regulation in the mammalian gut.

Host lung microbiota promotes malaria-associated acute respiratory distress syndrome.

Severe malaria can manifest itself with a variety of well-recognized clinical phenotypes that are highly predictive of death - severe anaemia, coma (cerebral malaria), multiple organ failure, and respiratory distress. The reasons why an infected individual develops one pathology rather than another remain poorly understood. Here we use distinct rodent models of infection to show that the host microbiota is a contributing factor for the development of respiratory distress syndrome and host mortality in the context of malaria infections (malaria-associated acute respiratory distress syndrome, MA-ARDS). We show that parasite sequestration in the lung results in sustained immune activation. Subsequent production of the anti-inflammatory cytokine IL-10 by T cells compromises microbial control, leading to severe lung disease. Notably, bacterial clearance with linezolid, an antibiotic commonly used in the clinical setting to control lung-associated bacterial infections, prevents MA-ARDS-associated lethality. Thus, we propose that the host's anti-inflammatory response to limit tissue damage can result in loss of microbial control, which promotes MA-ARDS. This must be considered when intervening against life-threatening respiratory complications.

Adaptive Potential of Epigenetic Switching During Adaptation to Fluctuating Environments.

Epigenetic regulation of gene expression allows for the emergence of distinct phenotypic states within the clonal population. Due to the instability of epigenetic inheritance, these phenotypes can intergenerationally switch between states in a stochastic manner. Theoretical studies of evolutionary dynamics predict that the phenotypic heterogeneity enabled by this rapid epigenetic switching between gene expression states would be favored under fluctuating environmental conditions, whereas genetic mutations, as a form of stable inheritance system, would be favored under a stable environment. To test this prediction, we engineered switcher and non-switcher yeast strains, in which the uracil biosynthesis gene URA3 is either continually expressed or switched on and off at two different rates (slow and fast switchers). Competitions between clones with an epigenetically controlled URA3 and clones without switching ability (SIR3 knockout) show that the switchers are favored in fluctuating environments. This occurs in conditions where the environments fluctuate at similar rates to the rate of switching. However, in stable environments, but also in environments with fluctuation frequency higher than the rate of switching, we observed that genetic changes dominated. Remarkably, epigenetic clones with a high, but not with a low, rate of switching can coexist with non-switchers even in a constant environment. Our study offers an experimental proof of concept that helps defining conditions of environmental fluctuation under which epigenetic switching provides an advantage.

SARS-CoV-2 introductions and early dynamics of the epidemic in Portugal.

Background: Genomic surveillance of SARS-CoV-2 in Portugal was rapidly implemented by the National Institute of Health in the early stages of the COVID-19 epidemic, in collaboration with more than 50 laboratories distributed nationwide. Methods: By applying recent phylodynamic models that allow integration of individual-based travel history, we reconstructed and characterized the spatio-temporal dynamics of SARS-CoV-2 introductions and early dissemination in Portugal. Results: We detected at least 277 independent SARS-CoV-2 introductions, mostly from European countries (namely the United Kingdom, Spain, France, Italy, and Switzerland), which were consistent with the countries with the highest connectivity with Portugal. Although most introductions were estimated to have occurred during early March 2020, it is likely that SARS-CoV-2 was silently circulating in Portugal throughout February, before the first cases were confirmed. Conclusions: Here we conclude that the earlier implementation of measures could have minimized the number of introductions and subsequent virus expansion in Portugal. This study lays the foundation for genomic epidemiology of SARS-CoV-2 in Portugal, and highlights the need for systematic and geographically-representative genomic surveillance.

Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution.

Background and objectives: To understand how organisms evolve, it is fundamental to study how mutations emerge and establish. Here, we estimated the rate of mutation accumulation of SARS-CoV-2 in vitro and investigated the repeatability of its evolution when facing a new cell type but no immune or drug pressures. Methodology: We performed experimental evolution with two strains of SARS-CoV-2, one carrying the originally described spike protein (CoV-2-D) and another carrying the D614G mutation that has spread worldwide (CoV-2-G). After 15 passages in Vero cells and whole genome sequencing, we characterized the spectrum and rate of the emerging mutations and looked for evidences of selection across the genomes of both strains. Results: From the frequencies of the mutations accumulated, and excluding the genes with signals of selection, we estimate a spontaneous mutation rate of 1.3 × 10 -6 ± 0.2 × 10-6 per-base per-infection cycle (mean across both lineages of SARS-CoV-2 ± 2SEM). We further show that mutation accumulation is larger in the CoV-2-D lineage and heterogeneous along the genome, consistent with the action of positive selection on the spike protein, which accumulated five times more mutations than the corresponding genomic average. We also observe the emergence of mutators in the CoV-2-G background, likely linked to mutations in the RNA-dependent RNA polymerase and/or in the error-correcting exonuclease protein. Conclusions and implications: These results provide valuable information on how spontaneous mutations emerge in SARS-CoV-2 and on how selection can shape its genome toward adaptation to new environments. Lay Summary: Each time a virus replicates inside a cell, errors (mutations) occur. Here, via laboratory propagation in cells originally isolated from the kidney epithelium of African green monkeys, we estimated the rate at which the SARS-CoV-2 virus mutates-an important parameter for understanding how it can evolve within and across humans. We also confirm the potential of its Spike protein to adapt to a new environment and report the emergence of mutators-viral populations where mutations occur at a significantly faster rate.

Diet leaves a genetic signature in a keystone member of the gut microbiota.

Switching from a low-fat and high-fiber diet to a Western-style high-fat and high-sugar diet causes microbiota imbalances that underlay many pathological conditions (i.e., dysbiosis). Although the effects of dietary changes on microbiota composition and functions are well documented, their impact in gut bacterial evolution remains unexplored. We followed the emergence of mutations in Bacteroides thetaiotaomicron, a prevalent fiber-degrading microbiota member, upon colonization of the murine gut under different dietary regimens. B. thetaiotaomicron evolved rapidly in the gut and Western-style diet selected for mutations that promote degradation of mucin-derived glycans. Periodic dietary changes caused fluctuations in the frequency of such mutations and were associated with metabolic shifts, resulting in the maintenance of higher intraspecies genetic diversity compared to constant dietary regimens. These results show that dietary changes leave a genetic signature in microbiome members and suggest that B. thetaiotaomicron genetic diversity could be a biomarker for dietary differences among individuals.

Competition dynamics in long-term propagations of Schizosaccharomyces pombe strain communities.

Experimental evolution studies with microorganisms such as bacteria and yeast have been an increasingly important and powerful tool to draw long-term inferences of how microbes interact. However, while several strains of the same species often exist in natural environments, many ecology and evolution studies in microbes are typically performed with isogenic populations of bacteria or yeast. In the present study, we firstly perform a genotypic and phenotypic characterization of two laboratory and eight natural strains of the yeast Schizosaccharomyces pombe. We then propagated, in a rich resource environment, yeast communities of 2, 3, 4, and 5 strains for hundreds of generations and asked which fitness-related phenotypes-maximum growth rate or relative competitive fitness-would better predict the outcome of a focal strain during the propagations. While the strain's growth rates would wrongly predict long-term coexistence, pairwise competitive fitness with a focal strain qualitatively predicted the success or extinction of the focal strain by a simple multigenotype population genetics model, given the initial community composition. Interestingly, we have also measured the competitive fitness of the ancestral and evolved communities by the end of the experiment (≈370 generations) and observed frequent maladaptation to the abiotic environment in communities with more than three members. Overall, our results aid establishing pairwise competitive fitness as good qualitative measurement of long-term community composition but also reveal a complex adaptive scenario when trying to predict the evolutionary outcome of those communities.

Molecular signatures of resource competition: Clonal interference favors ecological diversification and can lead to incipient speciation.

Microbial ecosystems harbor an astonishing diversity that can persist for long times. To understand how such diversity is structured and maintained, ecological and evolutionary processes need to be integrated at similar timescales. Here, we study a model of resource competition that allows for evolution via de novo mutation, and focus on rapidly adapting asexual populations with large mutational inputs, as typical of many bacteria species. We characterize the adaptation and diversification of an initially maladapted population and show how the eco-evolutionary dynamics are shaped by the interaction between simultaneously emerging lineages - clonal interference. We find that in large populations, more intense clonal interference can foster diversification under sympatry, increasing the probability that phenotypically and genetically distinct clusters coexist. In smaller populations, the accumulation of deleterious and compensatory mutations can push further the diversification process and kick-start speciation. Our findings have implications beyond microbial populations, providing novel insights about the interplay between ecology and evolution in clonal populations.

The Selective Advantage of the lac Operon for Escherichia coli Is Conditional on Diet and Microbiota Composition.

The lac operon is one of the best known gene regulatory circuits and constitutes a landmark example of how bacteria tune their metabolism to nutritional conditions. It is nearly ubiquitous in Escherichia coli strains justifying the use of its phenotype, the ability to consume lactose, for species identification. Lactose is the primary sugar found in milk, which is abundant in mammals during the first weeks of life. However, lactose is virtually non-existent after the weaning period, with humans being an exception as many consume dairy products throughout their lives. The absence of lactose during adulthood in most mammals and the rarity of lactose in the environment, means that the selective pressure for maintaining the lac operon could be weak for long periods of time. Despite the ability to metabolize lactose being a hallmark of E. coli's success when colonizing its primary habitat, the mammalian intestine, the selective value of this trait remains unknown in this ecosystem during adulthood. Here we determine the competitive advantage conferred by the lac operon to a commensal strain of E. coli when colonizing the mouse gut. We find that its benefit, which can be as high as 11%, is contingent on the presence of lactose in the diet and on the presence of other microbiota members in the gut, but the operon is never deleterious. These results help explaining the pervasiveness of the lac operon in E. coli, but also its polymorphism, as lac-negative E. coli strains albeit rare can naturally occur in the gut.

DNA Breaks-Mediated Fitness Cost Reveals RNase HI as a New Target for Selectively Eliminating Antibiotic-Resistant Bacteria.

Antibiotic resistance often generates defects in bacterial growth called fitness cost. Understanding the causes of this cost is of paramount importance, as it is one of the main determinants of the prevalence of resistances upon reducing antibiotics use. Here we show that the fitness costs of antibiotic resistance mutations that affect transcription and translation in Escherichia coli strongly correlate with DNA breaks, which are generated via transcription-translation uncoupling, increased formation of RNA-DNA hybrids (R-loops), and elevated replication-transcription conflicts. We also demonstrated that the mechanisms generating DNA breaks are repeatedly targeted by compensatory evolution, and that DNA breaks and the cost of resistance can be increased by targeting the RNase HI, which specifically degrades R-loops. We further show that the DNA damage and thus the fitness cost caused by lack of RNase HI function drive resistant clones to extinction in populations with high initial frequency of resistance, both in laboratory conditions and in a mouse model of gut colonization. Thus, RNase HI provides a target specific against resistant bacteria, which we validate using a repurposed drug. In summary, we revealed key mechanisms underlying the fitness cost of antibiotic resistance mutations that can be exploited to specifically eliminate resistant bacteria.

Rampant loss of social traits during domestication of a Bacillus subtilis natural isolate.

Most model bacteria have been domesticated in laboratory conditions. Yet, the tempo with which a natural isolate diverges from its ancestral phenotype under domestication to a novel laboratory environment is poorly understood. Such knowledge, however is essential to understanding the rate of evolution, the time scale over which a natural isolate can be propagated without loss of its natural adaptive traits, and the reliability of experimental results across labs. Using experimental evolution, phenotypic assays, and whole-genome sequencing, we show that within a week of propagation in a common laboratory environment, a natural isolate of Bacillus subtilis acquires mutations that cause changes in a multitude of traits. A single adaptive mutational step in the gene coding for the transcriptional regulator DegU impairs a DegU-dependent positive autoregulatory loop and leads to loss of robust biofilm architecture, impaired swarming motility, reduced secretion of exoproteases, and to changes in the dynamics of sporulation across environments. Importantly, domestication also resulted in improved survival when the bacteria face pressure from cells of the innate immune system. These results show that degU is a target for mutations during domestication and underscores the importance of performing careful and extremely short-term propagations of natural isolates to conserve the traits encoded in their original genomes.

Massive dissemination of a SARS-CoV-2 Spike Y839 variant in Portugal.

Genomic surveillance of SARS-CoV-2 was rapidly implemented in Portugal by the National Institute of Health in collaboration with a nationwide consortium of >50 hospitals/laboratories. Here, we track the geotemporal spread of a SARS-CoV-2 variant with a mutation (D839Y) in a potential host-interacting region involving the Spike fusion peptide, which is a target motif of anti-viral drugs that plays a key role in SARS-CoV-2 infectivity. The Spike Y839 variant was most likely imported from Italy in mid-late February and massively disseminated in Portugal during the early epidemic, becoming prevalent in the Northern and Central regions of Portugal where it represented 22% and 59% of the sampled genomes, respectively, by 30 April. Based on our high sequencing sampling during the early epidemics [15.5% (1275/8251) and 6.0% (1500/24987) of all confirmed cases until the end of March and April, respectively], we estimate that, between 14 March and 9 April (covering the epidemic exponential phase) the relative frequency of the Spike Y839 variant increased at a rate of 12.1% (6.1%-18.2%, CI 95%) every three days, being potentially associated with 24.8% (20.8-29.7%, CI 95%; 3177-4542 cases, CI 95%) of all COVID-19 cases in Portugal during this period. Our data supports population/epidemiological (founder) effects contributing to the Y839 variant superspread. The potential existence of selective advantage is also discussed, although experimental validation is required. Despite huge differences in genome sampling worldwide, SARS-CoV-2 Spike D839Y has been detected in 13 countries in four continents, supporting the need for close surveillance and functional assays of Spike variants.