Compensatory relationship between low-complexity regions and gene paralogy in the evolution of prokaryotes.

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life.

Rapid evolutionary repair by secondary perturbation of a primary disrupted transcriptional network.

The discrete steps of transcriptional rewiring have been proposed to occur neutrally to ensure steady gene expression under stabilizing selection. A conflict-free switch of a regulon between regulators may require an immediate compensatory evolution to minimize deleterious effects. Here, we perform an evolutionary repair experiment on the Lachancea kluyveri yeast sef1Δ mutant using a suppressor development strategy. Complete loss of SEF1 forces cells to initiate a compensatory process for the pleiotropic defects arising from misexpression of TCA cycle genes. Using different selective conditions, we identify two adaptive loss-of-function mutations of IRA1 and AZF1. Subsequent analyses show that Azf1 is a weak transcriptional activator regulated by the Ras1-PKA pathway. Azf1 loss-of-function triggers extensive gene expression changes responsible for compensatory, beneficial, and trade-off phenotypes. The trade-offs can be alleviated by higher cell density. Our results not only indicate that secondary transcriptional perturbation provides rapid and adaptive mechanisms potentially stabilizing the initial stage of transcriptional rewiring but also suggest how genetic polymorphisms of pleiotropic mutations could be maintained in the population.

Sexually antagonistic coevolution between the sex chromosomes of Drosophila melanogaster.

Antagonistic interactions between the sexes are important drivers of evolutionary divergence. Interlocus sexual conflict is generally described as a conflict between alleles at two interacting loci whose identity and genomic location are arbitrary, but with opposite fitness effects in each sex. We build on previous theory by suggesting that when loci under interlocus sexual conflict are located on the sex chromosomes it can lead to cycles of antagonistic coevolution between them and therefore between the sexes. We tested this hypothesis by performing experimental crosses using Drosophila melanogaster where we reciprocally exchanged the sex chromosomes between five allopatric wild-type populations in a round-robin design. Disrupting putatively coevolved sex chromosome pairs resulted in increased male reproductive success in 16 of 20 experimental populations (10 of which were individually significant), but also resulted in lower offspring egg-to-adult viability that affected both male and female fitness. After 25 generations of experimental evolution these sexually antagonistic fitness effects appeared to be resolved. To formalize our hypothesis, we developed population genetic models of antagonistic coevolution using fitness expressions based on our empirical results. Our model predictions support the conclusion that antagonistic coevolution between the sex chromosomes is plausible under the fitness effects observed in our experiments. Together, our results lend both empirical and theoretical support to the idea that cycles of antagonistic coevolution can occur between sex chromosomes and illustrate how this process, in combination with autosomal coadaptation, may drive genetic and phenotypic divergence between populations.

The Effect of the Stringent Response and Oxidative Stress Response on Fitness Costs of De Novo Acquisition of Antibiotic Resistance.

Resistance evolution during exposure to non-lethal levels of antibiotics is influenced by various stress responses of bacteria which are known to affect growth rate. Here, we aim to disentangle how the interplay between resistance development and associated fitness costs is affected by stress responses. We performed de novo resistance evolution of wild-type strains and single-gene knockout strains in stress response pathways using four different antibiotics. Throughout resistance development, the increase in minimum inhibitory concentration (MIC) is accompanied by a gradual decrease in growth rate, most pronounced in amoxicillin or kanamycin. By measuring biomass yield on glucose and whole-genome sequences at intermediate and final time points, we identified two patterns of how the stress responses affect the correlation between MIC and growth rate. First, single-gene knockout E. coli strains associated with reactive oxygen species (ROS) acquire resistance faster, and mutations related to antibiotic permeability and pumping out occur earlier. This increases the metabolic burden of resistant bacteria. Second, the ΔrelA knockout strain, which has reduced (p)ppGpp synthesis, is restricted in its stringent response, leading to diminished growth rates. The ROS-related mutagenesis and the stringent response increase metabolic burdens during resistance development, causing lower growth rates and higher fitness costs.

Pairs of Mutually Compensatory Frameshifting Mutations Contribute to Protein Evolution.

Insertions and deletions of lengths not divisible by 3 in protein-coding sequences cause frameshifts that usually induce premature stop codons and may carry a high fitness cost. However, this cost can be partially offset by a second compensatory indel restoring the reading frame. The role of such pairs of compensatory frameshifting mutations (pCFMs) in evolution has not been studied systematically. Here, we use whole-genome alignments of protein-coding genes of 100 vertebrate species, and of 122 insect species, studying the prevalence of pCFMs in their divergence. We detect a total of 624 candidate pCFM genes; six of them pass stringent quality filtering, including three human genes: RAB36, ARHGAP6, and NCR3LG1. In some instances, amino acid substitutions closely predating or following pCFMs restored the biochemical similarity of the frameshifted segment to the ancestral amino acid sequence, possibly reducing or negating the fitness cost of the pCFM. Typically, however, the biochemical similarity of the frameshifted sequence to the ancestral one was not higher than the similarity of a random sequence of a protein-coding gene to its frameshifted version, indicating that pCFMs can uncover radically novel regions of protein space. In total, pCFMs represent an appreciable and previously overlooked source of novel variation in amino acid sequences.

Rapid Decline of Ceftazidime Resistance in Antibiotic-Free and Sublethal Environments Is Contingent on Genetic Background.

Trade-offs of antibiotic resistance evolution, such as fitness cost and collateral sensitivity (CS), could be exploited to drive evolution toward antibiotic susceptibility. Decline of resistance may occur when resistance to other drug leads to CS to the first one and when compensatory mutations, or genetic reversion of the original ones, reduce fitness cost. Here we describe the impact of antibiotic-free and sublethal environments on declining ceftazidime resistance in different Pseudomonas aeruginosa resistant mutants. We determined that decline of ceftazidime resistance occurs within 450 generations, which is caused by newly acquired mutations and not by reversion of the original ones, and that the original CS of these mutants is preserved. In addition, we observed that the frequency and degree of this decline is contingent on genetic background. Our results are relevant to implement evolution-based therapeutic approaches, as well as to redefine global policies of antibiotic use, such as drug cycling.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves.

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the "nuclear compensation hypothesis," a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with "doubly uniparental" mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Epistasis-Driven Evolution of the SARS-CoV-2 Secondary Structure.

Epistasis is an evolutionary phenomenon whereby the fitness effect of a mutation depends on the genetic background in which it arises. A key source of epistasis in an RNA molecule is its secondary structure, which contains functionally important topological motifs held together by hydrogen bonds between Watson-Crick (WC) base pairs. Here we study epistasis in the secondary structure of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by examining properties of derived alleles arising from substitution mutations at ancestral WC base-paired and unpaired (UP) sites in 15 conserved topological motifs across the genome. We uncover fewer derived alleles and lower derived allele frequencies at WC than at UP sites, supporting the hypothesis that modifications to the secondary structure are often deleterious. At WC sites, we also find lower derived allele frequencies for mutations that abolish base pairing than for those that yield G·U "wobbles," illustrating that weak base pairing can partially preserve the integrity of the secondary structure. Last, we show that WC sites under the strongest epistatic constraint reside in a three-stemmed pseudoknot motif that plays an essential role in programmed ribosomal frameshifting, whereas those under the weakest epistatic constraint are located in 3' UTR motifs that regulate viral replication and pathogenicity. Our findings demonstrate the importance of epistasis in the evolution of the SARS-CoV-2 secondary structure, as well as highlight putative structural and functional targets of different forms of natural selection.

Natural variants suppress mutations in hundreds of essential genes.

The consequence of a mutation can be influenced by the context in which it operates. For example, loss of gene function may be tolerated in one genetic background, and lethal in another. The extent to which mutant phenotypes are malleable, the architecture of modifiers and the identities of causal genes remain largely unknown. Here, we measure the fitness effects of ~ 1,100 temperature-sensitive alleles of yeast essential genes in the context of variation from ten different natural genetic backgrounds and map the modifiers for 19 combinations. Altogether, fitness defects for 149 of the 580 tested genes (26%) could be suppressed by genetic variation in at least one yeast strain. Suppression was generally driven by gain-of-function of a single, strong modifier gene, and involved both genes encoding complex or pathway partners suppressing specific temperature-sensitive alleles, as well as general modifiers altering the effect of many alleles. The emerging frequency of suppression and range of possible mechanisms suggest that a substantial fraction of monogenic diseases could be managed by modulating other gene products.

Systematic analysis of bypass suppression of essential genes.

Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole-genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane-associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line-specific essentiality in the Cancer Dependency Map (DepMap) project.

Evolutionary Rate Correlation between Mitochondrial-Encoded and Mitochondria-Associated Nuclear-Encoded Proteins in Insects.

The mitochondrion is a pivotal organelle for energy production, and includes components encoded by both the mitochondrial and nuclear genomes. Functional and evolutionary interactions are expected between the nuclear- and mitochondrial-encoded components. The topic is of broad interest in biology, with implications to genetics, evolution, and medicine. Here, we compare the evolutionary rates of mitochondrial proteins and ribosomal RNAs to rates of mitochondria-associated nuclear-encoded proteins, across the major orders of holometabolous insects. There are significant evolutionary rate correlations (ERCs) between mitochondrial-encoded and mitochondria-associated nuclear-encoded proteins, which are likely driven by different rates of mitochondrial sequence evolution and correlated changes in the interacting nuclear-encoded proteins. The pattern holds after correction for phylogenetic relationships and considering protein conservation levels. Correlations are stronger for both nuclear-encoded OXPHOS proteins that are in contact with mitochondrial OXPHOS proteins and for nuclear-encoded mitochondrial ribosomal amino acids directly contacting the mitochondrial rRNAs. We find that ERC between mitochondrial- and nuclear-encoded proteins is a strong predictor of nuclear-encoded proteins known to interact with mitochondria, and ERC shows promise for identifying new candidate proteins with mitochondrial function. Twenty-three additional candidate nuclear-encoded proteins warrant further study for mitochondrial function based on this approach, including proteins in the minichromosome maintenance helicase complex.

Extremely fast amelioration of plasmid fitness costs by multiple functionally diverse pathways.

The acquisition of plasmids is often accompanied by fitness costs such that compensatory evolution is required to allow plasmid survival, but it is unclear whether compensatory evolution can be extensive or rapid enough to maintain plasmids when they are very costly. The mercury-resistance plasmid pQBR55 drastically reduced the growth of its host, Pseudomonas fluorescens SBW25, immediately after acquisition, causing a small colony phenotype. However, within 48 h of growth on agar plates we observed restoration of the ancestral large colony morphology, suggesting that compensatory mutations had occurred. Relative fitness of these evolved strains, in lab media and in soil microcosms, varied between replicates, indicating different mutational mechanisms. Using genome sequencing we identified that restoration was associated with chromosomal mutations in either a hypothetical DNA-binding protein PFLU4242, RNA polymerase or the GacA/S two-component system. Targeted deletions in PFLU4242, gacA or gacS recapitulated the ameliorated phenotype upon plasmid acquisition, indicating three distinct mutational pathways to compensation. Our data shows that plasmid compensatory evolution is fast enough to allow survival of a plasmid despite it imposing very high fitness costs upon its host, and indeed may regularly occur during the process of isolating and selecting individual plasmid-containing clones.

Mechanisms of suppression: The wiring of genetic resilience.

Recent analysis of genome sequences has identified individuals that are healthy despite carrying severe disease-associated mutations. A possible explanation is that these individuals carry a second genomic perturbation that can compensate for the detrimental effects of the disease allele, a phenomenon referred to as suppression. In model organisms, suppression interactions are generally divided into two classes: genomic suppressors which are secondary mutations in the genome that bypass a mutant phenotype, and dosage suppression interactions in which overexpression of a suppressor gene rescues a mutant phenotype. Here, we describe the general properties of genomic and dosage suppression, with an emphasis on the budding yeast. We propose that suppression interactions between genetic variants are likely relevant for determining the penetrance of human traits. Consequently, an understanding of suppression mechanisms may guide the discovery of protective variants in healthy individuals that carry disease alleles, which could direct the rational design of new therapeutics.

Coadaptation of mitochondrial and nuclear genes, and the cost of mother's curse.

Strict maternal inheritance renders the mitochondrial genome susceptible to accumulating mutations that harm males, but are otherwise benign or beneficial for females. This 'mother's curse' effect can degrade male survival and fertility if unopposed by counteracting evolutionary processes. Coadaptation between nuclear and mitochondrial genomes-with nuclear genes evolving to compensate for male-harming mitochondrial substitutions-may ultimately resolve mother's curse. However, males are still expected to incur a transient fitness cost during mito-nuclear coevolution, and it remains unclear how severe such costs should be. We present a population genetic analysis of mito-nuclear coadaptation to resolve mother's curse effects, and show that the magnitude of the 'male mitochondrial load'-the negative impact of mitochondrial substitutions on male fitness components-may be large, even when genetic variation for compensatory evolution is abundant. We also find that the male load is surprisingly sensitive to population size: male fitness costs of mito-nuclear coevolution are particularly pronounced in both small and large populations, and minimized in populations of intermediate size. Our results reveal complex interactions between demography and genetic constraints during the resolution of mother's curse, suggesting potentially widespread species differences in susceptibility to mother's curse effects.

Migration promotes plasmid stability under spatially heterogeneous positive selection.

Bacteria-plasmid associations can be mutualistic or antagonistic depending on the strength of positive selection for plasmid-encoded genes, with contrasting outcomes for plasmid stability. In mutualistic environments, plasmids are swept to high frequency by positive selection, increasing the likelihood of compensatory evolution to ameliorate the plasmid cost, which promotes long-term stability. In antagonistic environments, plasmids are purged by negative selection, reducing the probability of compensatory evolution and driving their extinction. Here we show, using experimental evolution of Pseudomonas fluorescens and the mercury-resistance plasmid, pQBR103, that migration promotes plasmid stability in spatially heterogeneous selection environments. Specifically, migration from mutualistic environments, by increasing both the frequency of the plasmid and the supply of compensatory mutations, stabilized plasmids in antagonistic environments where, without migration, they approached extinction. These data suggest that spatially heterogeneous positive selection, which is common in natural environments, coupled with migration helps to explain the stability of plasmids and the ecologically important genes that they encode.

Plasmid stability is enhanced by higher-frequency pulses of positive selection.

Plasmids accelerate bacterial adaptation by sharing ecologically important traits between lineages. However, explaining plasmid stability in bacterial populations is challenging owing to their associated costs. Previous theoretical and experimental studies suggest that pulsed positive selection may explain plasmid stability by favouring gene mobility and promoting compensatory evolution to ameliorate plasmid cost. Here we test how the frequency of pulsed positive selection affected the dynamics of a mercury-resistance plasmid, pQBR103, in experimental populations of Pseudomonas fluorescens SBW25. Plasmid dynamics varied according to the frequency of Hg2+ positive selection: in the absence of Hg2+ plasmids declined to low frequency, whereas pulses of Hg2+ selection allowed plasmids to sweep to high prevalence. Compensatory evolution to ameliorate the cost of plasmid carriage was widespread across the entire range of Hg2+ selection regimes, including both constant and pulsed Hg2+ selection. Consistent with theoretical predictions, gene mobility via conjugation appeared to play a greater role in promoting plasmid stability under low-frequency pulses of Hg2+ selection. However, upon removal of Hg2+ selection, plasmids which had evolved under low-frequency pulse selective regimes declined over time. Our findings suggest that temporally variable selection environments, such as those created during antibiotic treatments, may help to explain the stability of mobile plasmid-encoded resistance.

Plasmid persistence: costs, benefits, and the plasmid paradox.

Plasmids are extrachromosomal DNA elements that can be found throughout bacteria, as well as in other domains of life. Nonetheless, the evolutionary processes underlying the persistence of plasmids are incompletely understood. Bacterial plasmids may encode genes for traits that are sometimes beneficial to their hosts, such as antimicrobial resistance, virulence, heavy metal tolerance, and the catabolism of unique nutrient sources. In the absence of selection for these traits, however, plasmids generally impose a fitness cost on their hosts. As such, plasmid persistence presents a conundrum: models predict that costly plasmids will be lost over time or that beneficial plasmid genes will be integrated into the host genome. However, laboratory and comparative studies have shown that plasmids can persist for long periods, even in the absence of positive selection. Several hypotheses have been proposed to explain plasmid persistence, including host-plasmid co-adaptation, plasmid hitchhiking, cross-ecotype transfer, and high plasmid transfer rates, but there is no clear evidence that any one model adequately resolves the plasmid paradox.

The Roles of Compensatory Evolution and Constraint in Aminoacyl tRNA Synthetase Evolution.

Mitochondrial protein translation requires interactions between transfer RNAs encoded by the mitochondrial genome (mt-tRNAs) and mitochondrial aminoacyl tRNA synthetase proteins (mt-aaRS) encoded by the nuclear genome. It has been argued that animal mt-tRNAs have higher deleterious substitution rates relative to their nuclear-encoded counterparts, the cytoplasmic tRNAs (cyt-tRNAs). This dynamic predicts elevated rates of compensatory evolution of mt-aaRS that interact with mt-tRNAs, relative to aaRS that interact with cyt-tRNAs (cyt-aaRS). We find that mt-aaRS do evolve at significantly higher rates (exemplified by higher dN and dN/dS) relative to cyt-aaRS, across mammals, birds, and Drosophila. While this pattern supports a model of compensatory evolution, the level at which a gene is expressed is a more general predictor of protein evolutionary rate. We find that gene expression level explains 10-56% of the variance in aaRS dN/dS, and that cyt-aaRS are more highly expressed in addition to having lower dN/dS values relative to mt-aaRS, consistent with more highly expressed genes being more evolutionarily constrained. Furthermore, we find no evidence of positive selection acting on either class of aaRS protein, as would be expected under a model of compensatory evolution. Nevertheless, the signature of faster mt-aaRS evolution persists in mammalian, but not bird or Drosophila, lineages after controlling for gene expression, suggesting some additional effect of compensatory evolution for mammalian mt-aaRS. We conclude that gene expression is the strongest factor governing differential amino acid substitution rates in proteins interacting with mitochondrial versus cytoplasmic factors, with important differences in mt-aaRS molecular evolution among taxonomic groups.

Amelioration of the Fitness Costs of Antibiotic Resistance Due To Reduced Outer Membrane Permeability by Upregulation of Alternative Porins.

The fitness cost of antibiotic resistance is a key parameter in determining the evolutionary success of resistant bacteria. Studies of the effect of antibiotic resistance on bacterial fitness are heavily biased toward target alterations. Here we investigated how the costs in the form of a severely impaired growth rate associated with resistance due to absence of two major outer membrane porins can be genetically compensated. We performed an evolution experiment with 16 lineages of a double mutant of Escherichia coli with the ompCF genes deleted, and reduced fitness and increased resistance to different classes of antibiotics, including the carbapenems ertapenem and meropenem. After serial passage for only 250 generations, the relative growth rate increased from 0.85 to 0.99 (susceptible wild type set to 1.0). Compensation of the costs followed two different adaptive pathways where upregulation of expression of alternative porins bypassed the need for functional OmpCF porins. The first compensatory mechanism involved mutations in the phoR and pstS genes, causing constitutive high-level expression of the PhoE porin. The second mechanism involved mutations in the hfq and chiX genes that disrupted Hfq-dependent small RNA regulation, causing overexpression of the ChiP porin. Although susceptibility was restored in compensated mutants with PhoE overexpression, evolved mutants with high ChiP expression maintained the resistance phenotype. Our findings may explain why porin composition is often altered in resistant clinical isolates and provide new insights into how bypass mechanisms may allow genetic adaptation to a common multidrug resistance mechanism.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes.

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.