Direct inference of the distribution of fitness effects of spontaneous mutations from recombinant inbred C. elegans mutation accumulation lines.

The distribution of fitness effects (DFE) of new mutations plays a central role in evolutionary biology. Estimates of the DFE from experimental Mutation Accumulation (MA) lines are compromised by the complete linkage disequilibrium (LD) between mutations in different lines. To reduce LD, we constructed two sets of recombinant inbred lines from a cross of two C. elegans MA lines. One set of lines ("RIAILs") was intercrossed for ten generations prior to ten generations of selfing; the second set of lines ("RILs") omitted the intercrossing. Residual LD in the RIAILs is much less than in the RILs, which affects the inferred DFE when the sets of lines are analyzed separately. The best-fit model estimated from all lines (RIAILs + RILs) infers a large fraction of mutations with positive effects (∼40%); models that constrain mutations to have negative effects fit much worse. The conclusion is the same using only the RILs. For the RIAILs, however, models that constrain mutations to have negative effects fit nearly as well as models that allow positive effects. When mutations in high LD are pooled into haplotypes, the inferred DFE becomes increasingly negative-skewed and leptokurtic. We conclude that the conventional wisdom - most mutations have effects near zero, a handful of mutations have effects that are substantially negative and mutations with positive effects are very rare - is likely correct, and that unless it can be shown otherwise, estimates of the DFE that infer a substantial fraction of mutations with positive effects are likely confounded by LD.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans.

Albendazole and ivermectin are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channel (GluCl) genes, but it is unknown whether these genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required loss of two GluCl genes (avr-14 and avr-15) or loss of three GluCl genes (avr-14, avr-15, and glc-1). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance.

Competitive fitness and homologous recombination of SARS-CoV-2 variants of concern.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to emerge and cocirculate in humans and wild animals. The factors driving the emergence and replacement of novel variants and recombinants remain incompletely understood. Herein, we comprehensively characterized the competitive fitness of SARS-CoV-2 wild type (WT) and three variants of concern (VOCs), Alpha, Beta and Delta, by coinfection and serial passaging assays in different susceptible cells. Deep sequencing analyses revealed cell-specific competitive fitness: the Beta variant showed enhanced replication fitness during serial passage in Caco-2 cells, whereas the WT and Alpha variant showed elevated fitness in Vero E6 cells. Interestingly, a high level of neutralizing antibody sped up competition and completely reshaped the fitness advantages of different variants. More importantly, single clone purification identified a significant proportion of homologous recombinants that emerged during the passage history, and immune pressure reduced the frequency of recombination. Interestingly, a recombination hot region located between nucleotide sites 22,995 and 28,866 of the viral genomes could be identified in most of the detected recombinants. Our study not only profiled the variable competitive fitness of SARS-CoV-2 under different conditions, but also provided direct experimental evidence of homologous recombination between SARS-CoV-2 viruses, as well as a model for investigating SARS-CoV-2 recombination.

Non-additive effects between genotypes: Implications for competitive fitness assays.

Competitive fitness assays are widely used in evolutionary biology and typically rely on a reference strain to compare different focal genotypes. This approach implicitly relies on the absence of interaction between the competing genotypes. In other words, the performance of the reference strain must not depend on the competitor. This report scrutinized this assumption by competing diverged Drosophila simulans populations against a common reference strain. We detected strong evidence for interaction between the competing genotypes: (1) Frequency-dependent selection was common with opposite effects in genetically diverged populations. (2) Temporal heterogeneity of fitness estimates, which can be partially attributed to a competitor-specific delay in the eclosion of the reference strain. We propose that this inconsistent behavior of the reference strain can be considered a specific case of a genotype × environment interaction. Focal populations could modify the environment of the reference strain, either indirectly by altering the microbiome composition and food availability or directly by genotype-specific cannibalism. Our results provide new insights into the interaction of diverged genotypes and have important implications for the interpretation of competitive fitness assays.

Inter-species competition of surface bacterial flora of pomegranate and their role in spoilage.

The surface of fruits is heterogenous in term of its microenvironment hence dictate the kind of microflora that develops during storage. A better understanding of spoilage organisms would lead to better preservation methods. The pomegranate was chosen, since its sturdy and spoils slow at room temperature and is ideal for studying fruit spoilage in-situ. In the current study we isolated organisms from fruit surface and study the spoilage and competition amongst microbial species. Total 17 unique bacterial isolates from pomegranate were identified. The 16S rRNA gene identification placed them in 8 major genera (Acinetobacter, Micrococcus, Pantoea, Microbacterium, Strenotrophomonas, Bacillus, Staphylococcus and Exiguobacterium). Competition assay among isolate suggested that Exiguobacterium is dominant species followed by Micrococcus, Pantoea and Bacillus. The consortium of 3 different combinations (5 bacteria each) of isolated bacteria showed the spoilage phenotype on pomegranate. Except for 3 bacterial isolates, the rest of the isolates produced any one or multiple enzymes associated with the food spoilage (cellulase, amylase, lactase, pectinase and protease). The isolates were checked for the presence of genes associated with antibiotic resistance and 78.9% of the tested micro-organisms were blaTEM positive. Aminoglycoside resistance genes were present in 10% of the tested microbes. This study demonstrated interspecies competition amongst spoilage organisms. This understanding of surface flora of fruit would give better insights to preserve fruits.

The Cell-Cell Communication Signal Indole Controls the Physiology and Interspecies Communication of Acinetobacter baumannii.

Many bacteria utilize quorum sensing (QS) to control group behavior in a cell density-dependent manner. Previous studies have demonstrated that Acinetobacter baumannii employs an N-acyl-L-homoserine lactone (AHL)-based QS system to control biological functions and virulence. Here, we report that indole controls biological functions, virulence and AHL signal production in A. baumannii. The biosynthesis of indole is performed by A1S_3160 (AbiS, Acinetobacter baumannii indole synthase), which is a novel indole synthase annotated as an alpha/beta hydrolase in A. baumannii. Heterologous expression of AbiS in an Escherichia coli indole-deficient mutant also rescued the production of indole by using a distinct biosynthetic pathway from the tryptophanase TnaA, which produces indole directly from tryptophan in E. coli. Moreover, we revealed that indole from A. baumannii reduced the competitive fitness of Pseudomonas aeruginosa by inhibiting its QS systems and type III secretion system (T3SS). As A. baumannii and P. aeruginosa usually coexist in human lungs, our results suggest the crucial roles of indole in both the bacterial physiology and interspecies communication. IMPORTANCE Acinetobacter baumannii is an important human opportunistic pathogen that usually causes high morbidity and mortality. It employs the N-acyl-L-homoserine lactone (AHL)-type quorum sensing (QS) system, AbaI/AbaR, to regulate biological functions and virulence. In this study, we found that A. baumannii utilizes another QS signal, indole, to modulate biological functions and virulence. It was further revealed that indole positively controls the production of AHL signals by regulating abaI expression at the transcriptional levels. Furthermore, indole represses the QS systems and type III secretion system (T3SS) of P. aeruginosa and enhances the competitive ability of A. baumannii. Together, our work describes a QS signaling network where a pathogen uses to control the bacterial physiology and pathogenesis, and the competitive ability in microbial community.

Absence of Ribosome Modulation Factor Alters Growth and Competitive Fitness of Escherichia coli.

During stationary phase in Escherichia coli, the expression of the ribosome modulation factor (RMF) protein participates in the dimerization of two 70S ribosomes, ultimately creating a 100S particle. 100S ribosomes are commonly thought to function to preserve ribosomes as growth ceases and cells begin to catabolize intracellular components, including proteins, during their transition into stationary phase. Here, we show that the rates of stationary-phase ribosomal degradation are increased in an rmf mutant strain that cannot produce 100S ribosomes, resulting in deficiencies in outgrowth upon reinoculation into fresh medium. Upon coinoculation in LB medium, the mutant exhibits a delay in entry into log phase, differences in growth rates, and an overall reduction in relative fitness during competition. Unexpectedly, the rmf mutant exhibited shorter generation times than wild-type cells during log phase, both in monoculture and during competition. These doubling times of ∼13 min suggest that failure to maintain ribosomal balance affects the control of cell division. Though the timing of entry into and exit from log phase is altered, 100S ribosomes are not essential for long-term viability of the rmf mutant when grown in monoculture. IMPORTANCE Ribosomes are the sole source in any cell for new protein synthesis that is vital to maintain life. While ribosomes are frequently consumed as sources of nutrients under low-nutrient conditions, some ribosomes appear to be preserved for later use. The failure to maintain the availability of these ribosomes can lead to a dire consequence upon the influx of new nutrients, as cells are unable to efficiently replenish their metabolic machinery. It is important to study the repercussions, consequences, and mechanisms of survival in cells that cannot properly maintain the availability of their ribosomes in order to better understand their mechanisms of survival during competition under nutrient-depleted conditions.

Meiotic cDNA libraries reveal gene truncations and mitochondrial proteins important for competitive fitness in Saccharomyces cerevisiae.

Gametogenesis is an evolutionarily conserved developmental program whereby a diploid progenitor cell undergoes meiosis and cellular remodeling to differentiate into haploid gametes, the precursors for sexual reproduction. Even in the simple eukaryotic organism Saccharomyces cerevisiae, the meiotic transcriptome is very rich and complex, thereby necessitating new tools for functional studies. Here, we report the construction of 5 stage-specific, inducible complementary DNA libraries from meiotic cells that represent over 84% of the genes found in the budding yeast genome. We employed computational strategies to detect endogenous meiotic transcript isoforms as well as library-specific gene truncations. Furthermore, we developed a robust screening pipeline to test the effect of each complementary DNA on competitive fitness. Our multiday proof-of-principle time course revealed 877 complementary DNAs that were detrimental for competitive fitness when overexpressed. The list included mitochondrial proteins that cause dose-dependent disruption of cellular respiration as well as library-specific gene truncations that expose a dominant negative effect on competitive growth. Together, these high-quality complementary DNA libraries provide an important tool for systematically identifying meiotic genes, transcript isoforms, and protein domains that are important for a specific biological function.

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.

Newly identified parasitic nematode beta-tubulin alleles confer resistance to benzimidazoles.

Infections by parasitic nematodes cause large health and economic burdens worldwide. We use anthelmintic drugs to reduce these infections. However, resistance to anthelmintic drugs is extremely common and increasing worldwide. It is essential to understand the mechanisms of resistance to slow its spread. Recently, four new parasitic nematode beta-tubulin alleles have been identified in benzimidazole (BZ) resistant parasite populations: E198I, E198K, E198T, and E198stop. These alleles have not been tested for the ability to confer resistance or for any effects that they might have on organismal fitness. We introduced these four new alleles into the sensitive C. elegans laboratory-adapted N2 strain and exposed these genome-edited strains to both albendazole and fenbendazole. We found that all four alleles conferred resistance to both BZ drugs. Additionally, we tested for fitness consequences in both control and albendazole conditions over seven generations in competitive fitness assays. We found that none of the edited alleles had deleterious effects on fitness in control conditions and that all four alleles conferred strong and equivalent fitness benefits in BZ drug conditions. Because it is unknown if previously validated alleles confer a dominant or recessive BZ resistance phenotype, we tested the phenotypes caused by five of these alleles and found that none of them conferred a dominant BZ resistance phenotype. Accurate measurements of resistance, fitness effects, and dominance caused by the resistance alleles allow for the generation of better models of population dynamics and facilitate control practices that maximize the efficacy of this critical anthelmintic drug class.

The Disordered C-Terminus of the Chaperone DnaK Increases the Competitive Fitness of Pseudomonas putida and Facilitates the Toxicity of GraT.

Chaperone proteins are crucial for proper protein folding and quality control, especially when cells encounter stress caused by non-optimal temperatures. DnaK is one of such essential chaperones in bacteria. Although DnaK has been well characterized, the function of its intrinsically disordered C-terminus has remained enigmatic as the deletion of this region has been shown to either enhance or reduce its protein folding ability. We have shown previously that DnaK interacts with toxin GraT of the GraTA toxin-antitoxin system in Pseudomonas putida. Interestingly, the C-terminal truncation of DnaK was shown to alleviate GraT-caused growth defects. Here, we aim to clarify the importance of DnaK in GraT activity. We show that DnaK increases GraT toxicity, and particularly important is the negatively charged motif in the DnaK C-terminus. Given that GraT has an intrinsically disordered N-terminus, the assistance of DnaK is probably needed for re-modelling the toxin structure. We also demonstrate that the DnaK C-terminal negatively charged motif contributes to the competitive fitness of P. putida at both high and optimal growth temperatures. Thus, our data suggest that the disordered C-terminal end of DnaK enhances the chaperone functionality.

Respiratory supercomplexes enhance electron transport by decreasing cytochrome c diffusion distance.

Respiratory chains are crucial for cellular energy conversion and consist of multi-subunit complexes that can assemble into supercomplexes. These structures have been intensively characterized in various organisms, but their physiological roles remain unclear. Here, we elucidate their function by leveraging a high-resolution structural model of yeast respiratory supercomplexes that allowed us to inhibit supercomplex formation by mutation of key residues in the interaction interface. Analyses of a mutant defective in supercomplex formation, which still contains fully functional individual complexes, show that the lack of supercomplex assembly delays the diffusion of cytochrome c between the separated complexes, thus reducing electron transfer efficiency. Consequently, competitive cellular fitness is severely reduced in the absence of supercomplex formation and can be restored by overexpression of cytochrome c. In sum, our results establish how respiratory supercomplexes increase the efficiency of cellular energy conversion, thereby providing an evolutionary advantage for aerobic organisms.

PhenoMIP: High-Throughput Phenotyping of Diverse Caenorhabditis elegans Populations via Molecular Inversion Probes.

Whether generated within a lab setting or isolated from the wild, variant alleles continue to be an important resource for decoding gene function in model organisms such as Caenorhabditis elegans With advances in massively parallel sequencing, multiple whole-genome sequenced (WGS) strain collections are now available to the research community. The Million Mutation Project (MMP) for instance, analyzed 2007 N2-derived, mutagenized strains. Individually, each strain averages ∼400 single nucleotide variants amounting to ∼80 protein-coding variants. The effects of these variants, however, remain largely uncharacterized and querying the breadth of these strains for phenotypic changes requires a method amenable to rapid and sensitive high-throughput analysis. Here we present a pooled competitive fitness approach to quantitatively phenotype subpopulations of sequenced collections via molecular inversion probes (PhenoMIP). We phenotyped the relative fitness of 217 mutant strains on multiple food sources and classified these into five categories. We also demonstrate on a subset of these strains, that their fitness defects can be genetically mapped. Overall, our results suggest that approximately 80% of MMP mutant strains may have a decreased fitness relative to the lab reference, N2 The costs of generating this form of analysis through WGS methods would be prohibitive while PhenoMIP analysis in this manner is accomplished at less than one-tenth of projected WGS costs. We propose methods for applying PhenoMIP to a broad range of population selection experiments in a cost-efficient manner that would be useful to the community at large.

Quantitative benzimidazole resistance and fitness effects of parasitic nematode beta-tubulin alleles.

Infections by parasitic nematodes inflict a huge burden on the health of humans and livestock throughout the world. Anthelmintic drugs are the first line of defense against these infections. Unfortunately, resistance to these drugs is rampant and continues to spread. To improve treatment strategies, we must understand the genetics and molecular mechanisms that underlie resistance. Studies of the fungus Aspergillus nidulans and the free-living nematode Caenorhabditis elegans discovered that a beta-tubulin gene is mutated in benzimidazole (BZ) resistant strains. In parasitic nematode populations, three beta-tubulin alleles, F167Y, E198A, and F200Y, have long been correlated with resistance. Additionally, improvements in sequencing technologies have identified new alleles - E198V, E198L, E198K, E198I, and E198Stop - also correlated with BZ resistance. However, none of these alleles have been proven to cause resistance. To empirically demonstrate this point, we independently introduced the F167Y, E198A, and F200Y alleles as well as two of the newly identified alleles, E198V and E198L, into the BZ susceptible C. elegans N2 genetic background using the CRISPR-Cas9 system. These genome-edited strains were exposed to both albendazole and fenbendazole to quantitatively measure animal responses to BZs. We used a range of concentrations for each BZ compound to define response curves and found that all five of the alleles conferred resistance to BZ compounds equal to a loss of the entire beta-tubulin gene. These results prove that the parasite beta-tubulin alleles cause resistance. The E198V allele is found at low frequencies along with the E198L allele in natural parasite populations, suggesting that it could affect fitness. We performed competitive fitness assays and demonstrated that the E198V allele reduces animal health, supporting the hypothesis that this allele might be less fit in field populations. Overall, we present a powerful platform to quantitatively assess anthelmintic resistance and effects of specific resistance alleles on organismal fitness in the presence or absence of the drug.

Three biofilm types produced by a model pseudomonad are differentiated by structural characteristics and fitness advantage.

Model bacterial biofilm systems suggest that bacteria produce one type of biofilm, which is then modified by environmental and physiological factors, although the diversification of developing populations might result in the appearance of adaptive mutants producing altered structures with improved fitness advantage. Here we compare the air-liquid (A-L) interface viscous mass (VM) biofilm produced by Pseudomonas fluorescens SBW25 and the wrinkly spreader (WS) and complementary biofilm-forming strain (CBFS) biofilm types produced by adaptive SBW25 mutants in order to better understand the link between these physical structures and the fitness advantage they provide in experimental microcosms. WS, CBFS and VM biofilms can be differentiated by strength, attachment levels and rheology, as well as by strain characteristics associated with biofilm formation. Competitive fitness assays demonstrate that they provide similar advantages under static growth conditions but respond differently to increasing levels of physical disturbance. Pairwise competitions between biofilms suggest that these strains must be competing for at least two growth-limiting resources at the A-L interface, most probably O2 and nutrients, although VM and CBFS cells located lower down in the liquid column might provide an additional fitness advantage through the colonization of a less competitive zone below the biofilm. Our comparison of different SBW25 biofilm types illustrates more generally how varied biofilm characteristics and fitness advantage could become among adaptive mutants arising from an ancestral biofilm-forming strain and raises the question of how significant these changes might be in a range of medical, biotechnological and industrial contexts where diversification and change may be problematic.

An Overexpression Experiment Does Not Support the Hypothesis That Avoidance of Toxicity Determines the Rate of Protein Evolution.

The misfolding avoidance hypothesis postulates that sequence mutations render proteins cytotoxic and therefore the higher the gene expression, the stronger the operation of selection against substitutions. This translates into prediction that relative toxicity of extant proteins is higher for those evolving faster. In the present experiment, we selected pairs of yeast genes which were paralogous but evolving at different rates. We expressed them artificially to high levels. We expected that toxicity would be higher for ones bearing more mutations, especially that overcrowding should rather exacerbate than reverse the already existing differences in misfolding rates. We did find that the applied mode of overexpression caused a considerable decrease in fitness and that the decrease was proportional to the amount of excessive protein. However, it was not higher for proteins which are normally expressed at lower levels (and have less conserved sequence). This result was obtained consistently, regardless whether the rate of growth or ability to compete in common cultures was used as a proxy for fitness. In additional experiments, we applied factors that reduce accuracy of translation or enhance structural instability of proteins. It did not change a consistent pattern of independence between the fitness cost caused by overexpression of a protein and the rate of its sequence evolution.

Marked Neurospora crassa Strains for Competition Experiments and Bayesian Methods for Fitness Estimates.

The filamentous fungus Neurospora crassa, a model microbial eukaryote, has a life cycle with many features that make it suitable for studying experimental evolution. However, it has lacked a general tool for estimating relative fitness of different strains in competition experiments. To remedy this need, we constructed N. crassa strains that contain a modified csr-1 locus and developed an assay for detecting the proportion of the marked strain using a post PCR high resolution melting assay. DNA extraction from spore samples can be performed on 96-well plates, followed by a PCR step, which allows many samples to be processed with ease. Furthermore, we suggest a Bayesian approach for estimating relative fitness from competition experiments that takes into account the uncertainty in measured strain proportions. We show that there is a fitness effect of the mating type locus, as mating type mat a has a higher competitive fitness than mat A The csr-1* marker also has a small fitness effect, but is still a suitable marker for competition experiments. As a proof of concept, we estimate the fitness effect of the qde-2 mutation, a gene in the RNA interference pathway, and show that its competitive fitness is lower than what would be expected from its mycelial growth rate alone.

Competitive fitness analysis using Convolutional Neural Network.

We developed a procedure for estimating competitive fitness by using Caenorhabditis elegans as a model organism and a Convolutional Neural Network (CNN) as a tool. Competitive fitness is usually the most informative fitness measure, and competitive fitness assays often rely on green fluorescent protein (GFP) marker strains. CNNs are a class of deep learning neural networks, which are well suited for image analysis and object classification. Our model analyses involved image classification of nematodes as wild-type vs. GFP-expressing, and counted both categories. The performance was analyzed with (i) precision and recall parameters, and (ii) comparison of the wild-type frequency calculated from the model against that obtained by visual scoring of the same images. The average precision and recall varied from 0.79 to 0.87 and from 0.84 to 0.92, respectively, depending on worm density in the images. Compared with manual counting, the model decreased counting time at least 20-fold while preventing human errors. Given the rapid development in the field of CNN, the model, which is fully available on GitHub, can be further optimized and adapted for other image-based uses.

Fitness Barriers Limit Reversion of a Proofreading-Deficient Coronavirus.

The 3'-to-5' exoribonuclease in coronavirus (CoV) nonstructural protein 14 (nsp14-ExoN) mediates RNA proofreading during genome replication. ExoN catalytic residues are arranged in three motifs: I (DE), II (E), and III (D). Alanine replacement of the motif I residues (AA-E-D; four nucleotide substitutions) in murine hepatitis virus (MHV) and severe acute respiratory syndrome (SARS)-CoV yields viable mutants with impaired replication and fitness, increased mutation rates, and attenuated virulence in vivo Despite these impairments, MHV- and SARS-CoV ExoN motif I AA mutants (ExoN-AA) have not reverted at motif I in diverse in vitro and in vivo environments, suggesting that profound fitness barriers prevent motif I reversion. To test this hypothesis, we engineered MHV-ExoN-AA with 1, 2, or 3 nucleotide mutations along genetic pathways to AA-to-DE reversion. We show that engineered intermediate revertants were viable but had no increased replication or competitive fitness compared to that of MHV-ExoN-AA. In contrast, a low-passage-number (passage 10 [P10]) MHV-ExoN-AA showed increased replication and competitive fitness without reversion of ExoN-AA. Finally, engineered reversion of ExoN-AA to ExoN-DE in the presence of ExoN-AA passage-adaptive mutations resulted in significant fitness loss. These results demonstrate that while reversion is possible, at least one alternative adaptive pathway is more rapidly advantageous than intermediate revertants and may alter the genetic background to render reversion detrimental to fitness. Our results provide an evolutionary rationale for lack of ExoN-AA reversion, illuminate potential multiprotein replicase interactions and coevolution, and support future studies aimed at stabilizing attenuated CoV ExoN-AA mutants.IMPORTANCE Coronaviruses encode an exoribonuclease (ExoN) that is important for viral replication, fitness, and virulence, yet coronaviruses with a defective ExoN (ExoN-AA) have not reverted under diverse experimental conditions. In this study, we identify multiple impediments to MHV-ExoN-AA reversion. We show that ExoN-AA reversion is possible but evolutionarily unfavorable. Instead, compensatory mutations outside ExoN-AA motif I are more accessible and beneficial than partial reversion. We also show that coevolution between replicase proteins over long-term passage partially compensates for ExoN-AA motif I but renders the virus inhospitable to a reverted ExoN. Our results reveal the evolutionary basis for the genetic stability of ExoN-inactivating mutations, illuminate complex functional and evolutionary relationships between coronavirus replicase proteins, and identify potential mechanisms for stabilization of ExoN-AA coronavirus mutants.

TCF1 and LEF1 Control Treg Competitive Survival and Tfr Development to Prevent Autoimmune Diseases.

CD4+ Foxp3+ T regulatory (Treg) cells are key players in preventing lethal autoimmunity. Tregs undertake differentiation processes and acquire diverse functional properties. However, how Treg's differentiation and functional specification are regulated remains incompletely understood. Here, we report that gradient expression of TCF1 and LEF1 distinguishes Tregs into three distinct subpopulations, particularly highlighting a subset of activated Treg (aTreg) cells. Treg-specific ablation of TCF1 and LEF1 renders the mice susceptible to systemic autoimmunity. TCF1 and LEF1 are dispensable for Treg's suppressive capacity but essential for maintaining a normal aTreg pool and promoting Treg's competitive survival. As a consequence, the development of T follicular regulatory (Tfr) cells, which are a subset of aTreg, is abolished in TCF1/LEF1-conditional knockout mice, leading to unrestrained T follicular helper (Tfh) and germinal center B cell responses. Thus, TCF1 and LEF1 act redundantly to control the maintenance and functional specification of Treg subsets to prevent autoimmunity.