The genomic landscape of compensatory evolution

The metabolic domestication syndrome of budding yeast

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

Alterations in sperm DNA methylation may as a mediator of paternal air pollution exposure and offspring birth outcomes: Insight from a birth cohort study

Paternal exposure to environmental risk factors influences the offspring health. This study aimed to evaluate the association between paternal air pollution exposure mediated by sperm DNA methylation and adverse birth outcomes in offspring. We recruited 1607 fertile men and their partners from 2014 to 2016 and collected semen samples to detect sperm DNA methylation. Multivariate linear regression and weighted quantile sum regression models were used to assess the associations between paternal air pollution exposure and offspring birth outcomes. A critical exposure window was identified. Reduced representation bisulfite sequencing was used to detect sperm DNA methylation. The results demonstrated that high paternal exposure to PM2.5 (β = -211.31, 95% CI: (-386.37, -36.24)), PM10 (β = -178.20, 95% CI: (-277.13, -79.27)), and NO2 (β = -84.22, 95% CI: (-165.86, -2.57)) was negatively associated with offspring's birthweight, especially in boys. Additionally, an early exposure window of 15-69 days before fertilization was recognized to be the key exposure window, which increased the risk of low birth weight and small for gestational age. Furthermore, paternal co-exposure to six air pollutants contributed to lower birthweight (β = -51.91, 95% CI: (-92.72, -11.10)) and shorter gestational age (β = -1.72, 95% CI: (-3.26, -0.17)) and PM2.5 was the most weighted pollutant. Paternal air pollution exposure resulted in 10,328 differentially methylated regions and the IGF2R gene was the key gene involved in the epigenetic process. These differentially methylated genes were predominantly associated with protein binding, transcriptional regulation, and DNA templating. These findings indicate that spermatogenesis is a susceptible window during which paternal exposure to air pollution affects sperm DNA methylation and the birth outcomes of offspring.

Interaction with a phage gene underlie costs of a β-lactamase

The fitness cost of an antibiotic resistance gene (ARG) can differ across host strains, creating refuges that allow the maintenance of an ARG in the absence of direct selection for its resistance phenotype. Despite the importance of such ARG-host interactions for predicting ARG dynamics, the basis of ARG fitness costs and their variability between hosts are not well understood. We determined the genetic basis of a host-dependent cost of a β-lactamase, blaTEM-116*, that conferred a significant cost in one Escherichia coli strain but was close to neutral in 11 other Escherichia spp. strains. Selection of a blaTEM-116*-encoding plasmid in the strain in which it initially had a high cost resulted in rapid and parallel compensation for that cost through mutations in a P1-like phage gene, relAP1. When the wild-type relAP1 gene was added to a strain in which it was not present and in which blaTEM-116* was neutral, it caused the ARG to become costly. Thus, relAP1 is both necessary and sufficient to explain blaTEM-116* costs in at least some host backgrounds. To our knowledge, these findings represent the first demonstrated case of the cost of an ARG being influenced by a genetic interaction with a phage gene. The interaction between a phage gene and a plasmid-borne ARG highlights the complexity of selective forces determining the maintenance and spread of ARGs and, by extension, encoding phage and plasmids in natural bacterial communities.IMPORTANCEAntibiotic resistance genes (ARGs) play a major role in the increasing problem of antibiotic resistance in clinically relevant bacteria. Selection of these genes occurs in the presence of antibiotics, but their eventual success also depends on the sometimes substantial costs they impose on host bacteria in antibiotic-free environments. We evolved an ARG that confers resistance to penicillin-type antibiotics in one host in which it did confer a cost and in one host in which it did not. We found that costs were rapidly and consistently reduced through parallel genetic changes in a gene encoded by a phage that was infecting the costly host. The unmutated version of this gene was sufficient to cause the ARG to confer a cost in a host in which it was originally neutral, demonstrating an antagonism between the two genetic elements and underlining the range and complexity of pressures determining ARG dynamics in natural populations.

Deep Conservation and Unexpected Evolutionary History of Neighboring lncRNAs MALAT1 and NEAT1

Long non-coding RNAs (lncRNAs) have begun to receive overdue attention for their regulatory roles in gene expression and other cellular processes. Although most lncRNAs are lowly expressed and tissue-specific, notable exceptions include MALAT1 and its genomic neighbor NEAT1, two highly and ubiquitously expressed oncogenes with roles in transcriptional regulation and RNA splicing. Previous studies have suggested that NEAT1 is found only in mammals, while MALAT1 is present in all gnathostomes (jawed vertebrates) except birds. Here we show that these assertions are incomplete, likely due to the challenges associated with properly identifying these two lncRNAs. Using phylogenetic analysis and structure-aware annotation of publicly available genomic and RNA-seq coverage data, we show that NEAT1 is a common feature of tetrapod genomes except birds and squamates. Conversely, we identify MALAT1 in representative species of all major gnathostome clades, including birds. Our in-depth examination of MALAT1, NEAT1, and their genomic context in a wide range of vertebrate species allows us to reconstruct the series of events that led to the formation of the locus containing these genes in taxa from cartilaginous fish to mammals. This evolutionary history includes the independent loss of NEAT1 in birds and squamates, since NEAT1 is found in the closest living relatives of both clades (crocodilians and tuataras, respectively). These data clarify the origins and relationships of MALAT1 and NEAT1 and highlight an opportunity to study the change and continuity in lncRNA structure and function over deep evolutionary time.

Evolutionary origin of germline pathogenic variants in human DNA mismatch repair genes

BACKGROUND

Mismatch repair (MMR) system is evolutionarily conserved for genome stability maintenance. Germline pathogenic variants (PVs) in MMR genes that lead to MMR functional deficiency are associated with high cancer risk. Knowing the evolutionary origin of germline PVs in human MMR genes will facilitate understanding the biological base of MMR deficiency in cancer. However, systematic knowledge is lacking to address the issue. In this study, we performed a comprehensive analysis to know the evolutionary origin of human MMR PVs.

METHODS

We retrieved MMR gene variants from the ClinVar database. The genomes of 100 vertebrates were collected from the UCSC genome browser and ancient human sequencing data were obtained through comprehensive data mining. Cross-species conservation analysis was performed based on the phylogenetic relationship among 100 vertebrates. Rescaled ancient sequencing data were used to perform variant calling for archeological analysis.

RESULTS

Using the phylogenetic approach, we traced the 3369 MMR PVs identified in modern humans in 99 non-human vertebrate genomes but found no evidence for cross-species conservation as the source for human MMR PVs. Using the archeological approach, we searched the human MMR PVs in over 5000 ancient human genomes dated from 45,045 to 100 years before present and identified a group of MMR PVs shared between modern and ancient humans mostly within 10,000 years with similar quantitative patterns.

CONCLUSION

Our study reveals that MMR PVs in modern humans were arisen within the recent human evolutionary history.

Disruption of Leishmania flagellum attachment zone architecture causes flagellum loss

Leishmania are flagellated eukaryotic parasites that cause leishmaniasis and are closely related to the other kinetoplastid parasites such as Trypanosoma brucei. In all these parasites there is a cell membrane invagination at the base of the flagellum called the flagellar pocket, which is tightly associated with and sculpted by cytoskeletal structures including the flagellum attachment zone (FAZ). The FAZ is a complex interconnected structure linking the flagellum to the cell body and has critical roles in cell morphogenesis, function and pathogenicity. However, this structure varies dramatically in size and organisation between these different parasites, suggesting changes in protein localisation and function. Here, we screened the localisation and function of the Leishmania orthologues of T. brucei FAZ proteins identified in the genome-wide protein tagging project TrypTag. We identified 27 FAZ proteins and our deletion analysis showed that deletion of two FAZ proteins in the flagellum, FAZ27 and FAZ34 resulted in a reduction in cell body size, and flagellum loss in some cells. Furthermore, after null mutant generation, we observed distinct and reproducible changes to cell shape, demonstrating the ability of the parasite to adapt to morphological perturbations resulting from gene deletion. This process of adaptation has important implications for the study of Leishmania mutants.

The Pseudomonas aeruginosa Resistome: Permanent and Transient Antibiotic Resistance, an Overview

One of the most concerning characteristics of Pseudomonas aeruginosa is its low susceptibility to several antibiotics of common use in clinics, as well as its facility to acquire increased resistance levels. Consequently, the study of the antibiotic resistance mechanisms of this bacterium is of relevance for human health. For such a study, different types of resistance should be distinguished. The intrinsic resistome is composed of a set of genes, present in the core genome of P. aeruginosa, which contributes to its characteristic, species-specific, phenotype of susceptibility to antibiotics. Acquired resistance refers to those genetic events, such as the acquisition of mutations or antibiotic resistance genes that reduce antibiotic susceptibility. Finally, antibiotic resistance can be transiently acquired in the presence of specific compounds or under some growing conditions. The current article provides information on methods currently used to analyze intrinsic, mutation-driven, and transient antibiotic resistance in P. aeruginosa.

An FGFR2 mutation as the potential cause of a new phenotype including early-onset osteoporosis and bone fractures: a case report

Osteoporosis is a systemic, multifactorial disorder of bone mineralization. Many factors contributing to the development of osteoporosis have been identified so far, including gender, age, nutrition, lifestyle, exercise, drug use, as well as a range of comorbidities. In addition to environmental and lifestyle factors, molecular genetic factors account for 50-85% of osteoporosis cases. For example, the vitamin D receptor (VDR), collagen type I (COL1), estrogen receptor (ER), apolypoprotein Е (ApoE), bone morphogenetic protein (BMP), and Low-density lipoprotein receptor-related protein 5 (LRP5) are all involved in the pathogenesis of osteoporosis. Among the candidate genes, the pathogenic variants in which are involved in the pathogenesis of osteoporosis is FGFR2. Additionally, FGFs/FGFRs-dependent signaling has been shown to regulate skeletal development and has been linked to a plethora of heritable disorders of the musculoskeletal system. In this study we present the clinical, biochemical and radiological findings, as well as results of molecular genetic testing of a 13-year-old male proband with heritable osteoporosis, arthralgia and multiple fractures and a family history of abnormal bone mineralization and fractures. Whole exome sequencing found a heterozygous previously undescribed variant in the FGFR2 gene (NM_000141.5) (GRCh37.p13 ENSG00000066468.16: g.123298133dup; ENST00000358487.5:c.722dup; ENSP00000351276.5:p.Asn241LysfsTer43). The same variant was found in two affected relatives. These data lead us to believe that the variant in FGFR2 found in our proband and his relatives could be related to their phenotype. Therefore, modern methods of molecular genetic testing can allow us to differentiate between osteogenesis imperfecta and other bone mineralization disorders.

Fitness Costs of Antibiotic Resistance Impede the Evolution of Resistance to Other Antibiotics

Antibiotic resistance is a major threat to global health, claiming the lives of millions every year. With a nearly dry antibiotic development pipeline, novel strategies are urgently needed to combat resistant pathogens. One emerging strategy is the use of sequential antibiotic therapy, postulated to reduce the rate at which antibiotic resistance evolves. Here, we use the soft agar gradient evolution (SAGE) system to carry out high-throughput in vitro bacterial evolution against antibiotic pressure. We find that evolution of resistance to the antibiotic chloramphenicol (CHL) severely affects bacterial fitness, slowing the rate at which resistance to the antibiotics nitrofurantoin and streptomycin emerges. In vitro acquisition of compensatory mutations in the CHL-resistant cells markedly improves fitness and nitrofurantoin adaptation rates but fails to restore rates to wild-type levels against streptomycin. Genome sequencing reveals distinct evolutionary paths to resistance in fitness-impaired populations, suggesting resistance trade-offs in favor of mitigation of fitness costs. We show that the speed of bacterial fronts in SAGE plates is a reliable indicator of adaptation rates and evolutionary trajectories to resistance. Identification of antibiotics whose mutational resistance mechanisms confer stable impairments may help clinicians prescribe sequential antibiotic therapies that are less prone to resistance evolution.

Global analysis of suppressor mutations that rescue human genetic defects

BACKGROUND

Genetic suppression occurs when the deleterious effects of a primary "query" mutation, such as a disease-causing mutation, are rescued by a suppressor mutation elsewhere in the genome.

METHODS

To capture existing knowledge on suppression relationships between human genes, we examined 2,400 published papers for potential interactions identified through either genetic modification of cultured human cells or through association studies in patients.

RESULTS

The resulting network encompassed 476 unique suppression interactions covering a wide spectrum of diseases and biological functions. The interactions frequently linked genes that operate in the same biological process. Suppressors were strongly enriched for genes with a role in stress response or signaling, suggesting that deleterious mutations can often be buffered by modulating signaling cascades or immune responses. Suppressor mutations tended to be deleterious when they occurred in absence of the query mutation, in apparent contrast with their protective role in the presence of the query. We formulated and quantified mechanisms of genetic suppression that could explain 71% of interactions and provided mechanistic insight into disease pathology. Finally, we used these observations to predict suppressor genes in the human genome.

CONCLUSIONS

The global suppression network allowed us to define principles of genetic suppression that were conserved across diseases, model systems, and species. The emerging frequency of suppression interactions among human genes and range of underlying mechanisms, together with the prevalence of suppression in model organisms, suggest that compensatory mutations may exist for most genetic diseases.

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.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.

Overexpression profiling reveals cellular requirements in the context of genetic backgrounds and environments

Overexpression can help life adapt to stressful environments, making an examination of overexpressed genes valuable for understanding stress tolerance mechanisms. However, a systematic study of genes whose overexpression is functionally adaptive (GOFAs) under stress has yet to be conducted. We developed a new overexpression profiling method and systematically identified GOFAs in Saccharomyces cerevisiae under stress (heat, salt, and oxidative). Our results show that adaptive overexpression compensates for deficiencies and increases fitness under stress, like calcium under salt stress. We also investigated the impact of different genetic backgrounds on GOFAs, which varied among three S. cerevisiae strains reflecting differing calcium and potassium requirements for salt stress tolerance. Our study of a knockout collection also suggested that calcium prevents mitochondrial outbursts under salt stress. Mitochondria-enhancing GOFAs were only adaptive when adequate calcium was available and non-adaptive when calcium was deficient, supporting this idea. Our findings indicate that adaptive overexpression meets the cell's needs for maximizing the organism's adaptive capacity in the given environment and genetic context.

Multiple phosphorylation of the Cdc48/p97 cofactor protein Shp1/p47 occurs upon cell stress in budding yeast

The homohexameric p97 complex, composed of Cdc48 subunits in yeast, is a crucial component of protein quality control pathways including ER-associated degradation. The complex acts to segregate protein complexes in an ATP-dependent manner, requiring the engagement of cofactor proteins that determine substrate specificity. The function of different Cdc48 cofactors and how they are regulated remains relatively poorly understood. In this study, we assess the phosphorylation of Cdc48 adaptor proteins, revealing a unique and distinctive phosphorylation pattern of Shp1/p47 that changed in response to TORC1 inhibition. Site-directed mutagenesis confirmed that this pattern corresponded to phosphorylation at residues S108 and S315 of Shp1, with the double-phosphorylated form becoming predominant upon TORC1 inhibition, ER-stress, and oxidative stress. Finally, we assessed candidate kinases and phosphatases responsible for Shp1 phosphorylation and identified two regulators. We found that cells lacking the kinase Mpk1/Slt2 show reduced Shp1 phosphorylation, whereas impaired PP1 phosphatase catalytic subunit (Glc7) activity resulted in increased Shp1 phosphorylation. Overall, these findings identify a phosphoregulation of Shp1 at multiple sites by Mpk1 kinase and PP1 phosphatase upon various stresses.

Experimental approaches to study evolutionary cell biology using yeasts

The past century has witnessed tremendous advances in understanding how cells function. Nevertheless, how cellular processes have evolved is still poorly understood. Many studies have highlighted surprising molecular diversity in how cells from diverse species execute the same processes, and advances in comparative genomics are likely to reveal much more molecular diversity than was believed possible until recently. Extant cells remain therefore the product of an evolutionary history that we vastly ignore. Evolutionary cell biology has emerged as a discipline aiming to address this knowledge gap by combining evolutionary, molecular, and cellular biology thinking. Recent studies have shown how even essential molecular processes, such as DNA replication, can undergo fast adaptive evolution under certain laboratory conditions. These developments open new lines of research where the evolution of cellular processes can be investigated experimentally. Yeasts naturally find themselves at the forefront of this research line. Not only do they allow the observation of fast evolutionary adaptation, but they also provide numerous genomic, synthetic, and cellular biology tools already developed by a large community. Here we propose that yeasts can serve as an "evolutionary cell lab" to test hypotheses, principles, and ideas in evolutionary cell biology. We discuss various experimental approaches available for this purpose, and how biology at large can benefit from them.

Complexity and self-organization in the evolution of cell polarization

Cellular life exhibits order and complexity, which typically increase over the course of evolution. Cell polarization is a well-studied example of an ordering process that breaks the internal symmetry of a cell by establishing a preferential axis. Like many cellular processes, polarization is driven by self-organization, meaning that the macroscopic pattern emerges as a consequence of microscopic molecular interactions at the biophysical level. However, the role of self-organization in the evolution of complex protein networks remains obscure. In this Review, we provide an overview of the evolution of polarization as a self-organizing process, focusing on the model species Saccharomyces cerevisiae and its fungal relatives. Moreover, we use this model system to discuss how self-organization might relate to evolutionary change, offering a shift in perspective on evolution at the microscopic scale.

Fitness of evolving bacterial populations is contingent on deep and shallow history but only shallow history creates predictable patterns

Long-term evolution experiments have tested the importance of genetic and environmental factors in influencing evolutionary outcomes. Differences in phylogenetic history, recent adaptation to distinct environments and chance events, all influence the fitness of a population. However, the interplay of these factors on a population's evolutionary potential remains relatively unexplored. We tracked the outcome of 2000 generations of evolution of four natural isolates of Escherichia coli bacteria that were engineered to also create differences in shallow history by adding previously identified mutations selected in a separate long-term experiment. Replicate populations started from each progenitor evolved in four environments. We found that deep and shallow phylogenetic histories both contributed significantly to differences in evolved fitness, though by different amounts in different selection environments. With one exception, chance effects were not significant. Whereas the effect of deep history did not follow any detectable pattern, effects of shallow history followed a pattern of diminishing returns whereby fitter ancestors had smaller fitness increases. These results are consistent with adaptive evolution being contingent on the interaction of several evolutionary forces but demonstrate that the nature of these interactions is not fixed and may not be predictable even when the role of chance is small.

Evolution of Habitat-Dependent Antibiotic Resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that usually causes difficult-to-treat infections due to its low intrinsic antibiotic susceptibility and outstanding capacity for becoming resistant to antibiotics. In addition, it has a remarkable metabolic versatility, being able to grow in different habitats, from natural niches to different and changing inpatient environments. Study of the environmental conditions that shape genetic and phenotypic changes of P. aeruginosa toward antibiotic resistance supposes a novelty, since experimental evolution assays are usually performed with well-defined antibiotics in regular laboratory growth media. Therefore, in this work we address the extent to which the nutrients’ availability may constrain the evolution of antibiotic resistance. We determined that P. aeruginosa genetic trajectories toward resistance to tobramycin, ceftazidime, and ceftazidime-avibactam are different when evolving in laboratory rich medium, urine, or synthetic sputum. Furthermore, our study, linking genotype with phenotype, showed a clear impact of each analyzed environment on both the fitness and resistance level associated with particular resistance mutations. This indicates that the phenotype associated with specific resistance mutations is variable and dependent on the bacterial metabolic state in each particular habitat. Our results support that the design of evolution-based strategies to tackle P. aeruginosa infections should be based on robust patterns of evolution identified within each particular infection and body location.

IMPORTANCE Predicting evolution toward antibiotic resistance (AR) and its associated trade-offs, such as collateral sensitivity, is important to design evolution-based strategies to tackle AR. However, the effect of nutrients' availability on such evolution, particularly those that can be found under in vivo infection conditions, has been barely addressed. We analyzed the evolutionary patterns of P. aeruginosa in the presence of antibiotics in different media, including urine and synthetic sputum, whose compositions are similar to the ones in infections, finding that AR evolution differs, depending on growth conditions. Furthermore, the representative mutants isolated under each condition tested render different AR levels and fitness costs, depending on nutrients’ availability, supporting the idea that environmental constraints shape the phenotypes associated with specific AR mutations. Consequently, the selection of AR mutations that render similar phenotypes is environment dependent. The analysis of evolution patterns toward AR requires studying growth conditions mimicking those that bacteria face during in vivo evolution.

Gene loss and compensatory evolution promotes the emergence of morphological novelties in budding yeast

Deleterious mutations are generally considered to be irrelevant for morphological evolution. However, they could be compensated by conditionally beneficial mutations, thereby providing access to new adaptive paths. Here we use high-dimensional phenotyping of laboratory-evolved budding yeast lineages to demonstrate that new cellular morphologies emerge exceptionally rapidly as a by-product of gene loss and subsequent compensatory evolution. Unexpectedly, the capacities for invasive growth, multicellular aggregation and biofilm formation also spontaneously evolve in response to gene loss. These multicellular phenotypes can be achieved by diverse mutational routes and without reactivating the canonical regulatory pathways. These ecologically and clinically relevant traits originate as pleiotropic side effects of compensatory evolution and have no obvious utility in the laboratory environment. The extent of morphological diversity in the evolved lineages is comparable to that of natural yeast isolates with diverse genetic backgrounds and lifestyles. Finally, we show that both the initial gene loss and subsequent compensatory mutations contribute to new morphologies, with their synergistic effects underlying specific morphological changes. We conclude that compensatory evolution is a previously unrecognized source of morphological diversity and phenotypic novelties.

Idiosyncratic epistasis leads to global fitness-correlated trends

Epistasis can markedly affect evolutionary trajectories. In recent decades, protein-level fitness landscapes have revealed extensive idiosyncratic epistasis among specific mutations. By contrast, other work has found ubiquitous and apparently nonspecific patterns of global diminishing-returns and increasing-costs epistasis among mutations across the genome. Here, we used a hierarchical CRISPR gene drive system to construct all combinations of 10 missense mutations from across the genome in budding yeast and measured their fitness in six environments. We show that the resulting fitness landscapes exhibit global fitness-correlated trends but that these trends emerge from specific idiosyncratic interactions. We thus provide experimental validation of recent theoretical work arguing that fitness-correlated trends can emerge as the generic consequence of idiosyncratic epistasis.

Rapid decline of ceftazidime resistance in antibiotic-free and sub-lethal 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.

The relative fitness of the de novo variants in general Lithuanian population vs. in individuals with intellectual disability

The effect of a variant on an organism is always multifaceted and can be considered from multiple perspectives-biochemical, medical, or evolutionary. However, the relationship between the effects of amino acid substitution on protein activity, human health, and an individual's evolutionary fitness is not trivial. We uncover that the general Lithuanian population is characterized by a "mirror reflection" of the de novo variant fitness effect, confirming the theory of neutrality. Meanwhile, in the group of individuals with intellectual disability, compared with the reference exome de novo variants significantly changed the composition of the amino acid. Therefore, it predicts that, both in terms of the number of amino acids and changes in their relative fitness, the structure of the proteins encoded by the studied amino acids undergo significant changes following the de novo variant, leading to possible changes in protein function associated with phenotypic traits. These results suggest that the analysis of relative fitness of exome sequences with de novo variants can predict the future phenotype. Therefore even in those cases, then only a few of all functional prediction analysis tools predict a variant as damaging, the negative relative fitness or even adaptability of the genome variant should be carefully evaluated considering both its direct function and the global background of the possible disease-associated mechanism regardless of the phenotype being studied.

Gene losses may contribute to subterranean adaptations in naked mole-rat and blind mole-rat

Background

Naked mole-rats (Heterocephalus glaber, NMRs) and blind mole-rats (Spalax galili, BMRs) are representative subterranean rodents that have evolved many extraordinary traits, including hypoxia tolerance, longevity, and cancer resistance. Although multiple candidate loci responsible for these traits have been uncovered by genomic studies, many of them are limited to functional changes to amino acid sequence and little is known about the contributions of other genetic events. To address this issue, we focused on gene losses (unitary pseudogenes) and systematically analyzed gene losses in NMRs and BMRs, aiming to elucidate the potential roles of pseudogenes in their adaptation to subterranean lifestyle.

Results

We obtained the pseudogene repertoires in NMRs and BMRs, as well as their respective aboveground relatives, guinea pigs and rats, on a genome-wide scale. As a result, 167, 139, 341, and 112 pseudogenes were identified in NMRs, BMRs, guinea pigs, and rats, respectively. Functional enrichment analysis identified 4 shared and 2 species-specific enriched functional groups (EFGs) in subterranean lineages. Notably, the pseudogenes in these EFGs might be associated with either regressive (e.g., visual system) or adaptive (e.g., altered DNA damage response) traits. In addition, several pseudogenes including TNNI3K and PDE5A might be associated with specific cardiac features observed in subterranean lineages. Interestingly, we observed 20 convergent gene losses in NMRs and BMRs. Given that the functional investigations of these genes are generally scarce, we provided functional evidence that independent loss of TRIM17 in NMRs and BMRs might be beneficial for neuronal survival under hypoxia, supporting the positive role of eliminating TRIM17 function in hypoxia adaptation. Our results also suggested that pseudogenes, together with positively selected genes, reinforced subterranean adaptations cooperatively.

Conclusions

Our study provides new insights into the molecular underpinnings of subterranean adaptations and highlights the importance of gene losses in mammalian evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01243-0.

Experimental Evolution Expands the Breadth of Adaptation to an Environmental Gradient Correlated With Genome Reduction

Whether and how adaptive evolution adjusts the breadth of adaptation in coordination with the genome are essential issues for connecting evolution with ecology. To address these questions, experimental evolution in five Escherichia coli strains carrying either the wild-type genome or a reduced genome was performed in a defined minimal medium (C0). The ancestral and evolved populations were subsequently subjected to fitness and chemical niche analyses across an environmental gradient with 29 combinations of eight chemical components of the minimal medium. The results showed that adaptation was achieved not only specific to the evolutionary condition (C0), but also generally, to the environmental gradient; that is, the breadth of adaptation to the eight chemical niches was expanded. The magnitudes of the adaptive improvement and the breadth increase were both correlated with genome reduction and were highly significant in two out of eight niches (i.e., glucose and sulfate). The direct adaptation-induced correlated adaptation to the environmental gradient was determined by only a few genome mutations. An additive increase in fitness associated with the stepwise fixation of mutations was consistently observed in the reduced genomes. In summary, this preliminary survey demonstrated that evolution finely tuned the breadth of adaptation correlated with genome reduction.

Mutation bias shapes the spectrum of adaptive substitutions

How do mutational biases influence the process of adaptation? A common assumption is that selection alone determines the course of adaptation from abundant preexisting variation. Yet, theoretical work shows broad conditions under which the mutation rate to a given type of variant strongly influences its probability of contributing to adaptation. Here we introduce a statistical approach to analyzing how mutation shapes protein sequence adaptation. Using large datasets from three different species, we show that the mutation spectrum has a proportional influence on the types of changes fixed in adaptation. We also show via computer simulations that a variety of factors can influence how closely the spectrum of adaptive substitutions reflects the spectrum of variants introduced by mutation.

Evolutionary adaptation often occurs by the fixation of beneficial mutations. This mode of adaptation can be characterized quantitatively by a spectrum of adaptive substitutions, i.e., a distribution for types of changes fixed in adaptation. Recent work establishes that the changes involved in adaptation reflect common types of mutations, raising the question of how strongly the mutation spectrum shapes the spectrum of adaptive substitutions. We address this question with a codon-based model for the spectrum of adaptive amino acid substitutions, applied to three large datasets covering thousands of amino acid changes identified in natural and experimental adaptation in Saccharomyces cerevisiae, Escherichia coli, and Mycobacterium tuberculosis. Using species-specific mutation spectra based on prior knowledge, we find that the mutation spectrum has a proportional influence on the spectrum of adaptive substitutions in all three species. Indeed, we find that by inferring the mutation rates that best explain the spectrum of adaptive substitutions, we can accurately recover the species-specific mutation spectra. However, we also find that the predictive power of the model differs substantially between the three species. To better understand these differences, we use population simulations to explore the factors that influence how closely the spectrum of adaptive substitutions mirrors the mutation spectrum. The results show that the influence of the mutation spectrum decreases with increasing mutational supply (Nμ) and that predictive power is strongly affected by the number and diversity of beneficial mutations.

Molecular species selectivity of lipid transport creates a mitochondrial sink for di-unsaturated phospholipids

Mitochondria depend on the import of phospholipid precursors for the biosynthesis of phosphatidylethanolamine (PE) and cardiolipin, yet the mechanism of their transport remains elusive. A dynamic lipidomics approach revealed that mitochondria preferentially import di-unsaturated phosphatidylserine (PS) for subsequent conversion to PE by the mitochondrial PS decarboxylase Psd1p. Several protein complexes tethering mitochondria to the endomembrane system have been implicated in lipid transport in yeast, including the endoplasmic reticulum (ER)-mitochondrial encounter structure (ERMES), ER-membrane complex (EMC), and the vacuole and mitochondria patch (vCLAMP). By limiting the availability of unsaturated phospholipids, we created conditions to investigate the mechanism of lipid transfer and the contributions of the tethering complexes in vivo. Under these conditions, inactivation of ERMES components or of the vCLAMP component Vps39p exacerbated accumulation of saturated lipid acyl chains, indicating that ERMES and Vps39p contribute to the mitochondrial sink for unsaturated acyl chains by mediating transfer of di-unsaturated phospholipids. These results support the concept that intermembrane lipid flow is rate-limited by molecular species-dependent lipid efflux from the donor membrane and driven by the lipid species' concentration gradient between donor and acceptor membrane.

Epistasis, aneuploidy, and functional mutations underlie evolution of resistance to induced microtubule depolymerization

Cells with blocked microtubule polymerization are delayed in mitosis, but eventually manage to proliferate despite substantial chromosome missegregation. While several studies have analyzed the first cell division after microtubule depolymerization, we have asked how cells cope long-term with microtubule impairment. We allowed 24 clonal populations of yeast cells with beta-tubulin mutations preventing proper microtubule polymerization, to evolve for ˜150 generations. At the end of the laboratory evolution experiment, cells had regained the ability to form microtubules and were less sensitive to microtubule-depolymerizing drugs. Whole-genome sequencing identified recurrently mutated genes, in particular for tubulins and kinesins, as well as pervasive duplication of chromosome VIII. Recreating these mutations and chromosome VIII disomy prior to evolution confirmed that they allow cells to compensate for the original mutation in beta-tubulin. Most of the identified mutations did not abolish function, but rather restored microtubule functionality. Analysis of the temporal order of resistance development in independent populations repeatedly revealed the same series of events: disomy of chromosome VIII followed by a single additional adaptive mutation in either tubulins or kinesins. Since tubulins are highly conserved among eukaryotes, our results have implications for understanding resistance to microtubule-targeting drugs widely used in cancer therapy.

Non-genetic and genetic rewiring underlie adaptation to hypomorphic alleles of an essential gene

Adaptive evolution to cellular stress is a process implicated in a wide range of biological and clinical phenomena. Two major routes of adaptation have been identified: non-genetic changes, which allow expression of different phenotypes in novel environments, and genetic variation achieved by selection of fitter phenotypes. While these processes are broadly accepted, their temporal and epistatic features in the context of cellular evolution and emerging drug resistance are contentious. In this manuscript, we generated hypomorphic alleles of the essential nuclear pore complex (NPC) gene NUP58. By dissecting early and long-term mechanisms of adaptation in independent clones, we observed that early physiological adaptation correlated with transcriptome rewiring and upregulation of genes known to interact with the NPC; long-term adaptation and fitness recovery instead occurred via focal amplification of NUP58 and restoration of mutant protein expression. These data support the concept that early phenotypic plasticity allows later acquisition of genetic adaptations to a specific impairment. We propose this approach as a genetic model to mimic targeted drug therapy in human cells and to dissect mechanisms of adaptation.

Ploidy and recombination proficiency shape the evolutionary adaptation to constitutive DNA replication stress

In haploid budding yeast, evolutionary adaptation to constitutive DNA replication stress alters three genome maintenance modules: DNA replication, the DNA damage checkpoint, and sister chromatid cohesion. We asked how these trajectories depend on genomic features by comparing the adaptation in three strains: haploids, diploids, and recombination deficient haploids. In all three, adaptation happens within 1000 generations at rates that are correlated with the initial fitness defect of the ancestors. Mutations in individual genes are selected at different frequencies in populations with different genomic features, but the benefits these mutations confer are similar in the three strains, and combinations of these mutations reproduce the fitness gains of evolved populations. Despite the differences in the selected mutations, adaptation targets the same three functional modules in strains with different genomic features, revealing a common evolutionary response to constitutive DNA replication stress.

Shortening of membrane lipid acyl chains compensates for phosphatidylcholine deficiency in choline-auxotroph yeast

Phosphatidylcholine (PC) is an abundant membrane lipid component in most eukaryotes, including yeast, and has been assigned multiple functions in addition to acting as building block of the lipid bilayer. Here, by isolating S. cerevisiae suppressor mutants that exhibit robust growth in the absence of PC, we show that PC essentiality is subject to cellular evolvability in yeast. The requirement for PC is suppressed by monosomy of chromosome XV or by a point mutation in the ACC1 gene encoding acetyl-CoA carboxylase. Although these two genetic adaptations rewire lipid biosynthesis in different ways, both decrease Acc1 activity, thereby reducing average acyl chain length. Consistently, soraphen A, a specific inhibitor of Acc1, rescues a yeast mutant with deficient PC synthesis. In the aneuploid suppressor, feedback inhibition of Acc1 through acyl-CoA produced by fatty acid synthase (FAS) results from upregulation of lipid synthesis. The results show that budding yeast regulates acyl chain length by fine-tuning the activities of Acc1 and FAS and indicate that PC evolved by benefitting the maintenance of membrane fluidity.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

De Novo Mutation and Rapid Protein (Co-)evolution during Meiotic Adaptation in Arabidopsis arenosa

A sudden shift in environment or cellular context necessitates rapid adaptation. A dramatic example is genome duplication, which leads to polyploidy. In such situations, the waiting time for new mutations might be prohibitive; theoretical and empirical studies suggest that rapid adaptation will largely rely on standing variation already present in source populations. Here, we investigate the evolution of meiosis proteins in Arabidopsis arenosa, some of which were previously implicated in adaptation to polyploidy, and in a diploid, habitat. A striking and unexplained feature of prior results was the large number of amino acid changes in multiple interacting proteins, especially in the relatively young tetraploid. Here, we investigate whether selection on meiosis genes is found in other lineages, how the polyploid may have accumulated so many differences, and whether derived variants were selected from standing variation. We use a range-wide sample of 145 resequenced genomes of diploid and tetraploid A. arenosa, with new genome assemblies. We confirmed signals of positive selection in the polyploid and diploid lineages they were previously reported in and find additional meiosis genes with evidence of selection. We show that the polyploid lineage stands out both qualitatively and quantitatively. Compared with diploids, meiosis proteins in the polyploid have more amino acid changes and a higher proportion affecting more strongly conserved sites. We find evidence that in tetraploids, positive selection may have commonly acted on de novo mutations. Several tests provide hints that coevolution, and in some cases, multinucleotide mutations, might contribute to rapid accumulation of changes in meiotic proteins.

Mutation E48K in PB1 Polymerase Subunit Improves Stability of a Candidate Live Attenuated Influenza B Virus Vaccine

Influenza B virus (IBV) is a major respiratory pathogen of humans, particularly in the elderly and children, and vaccines are the most effective way to control it. In previous work, incorporation of two mutations (E580G, S660A) along with the addition of an HA epitope tag in the PB1 segment of B/Brisbane/60/2008 (B/Bris) resulted in an attenuated strain that was safe and effective as a live attenuated vaccine. A third attempted mutation (K391E) in PB1 was not always stable. Interestingly, viruses that maintained the K391E mutation were associated with the mutation E48K. To explore the contribution of the E48K mutation to stability of the K391E mutation, a vaccine candidate was generated by inserting both mutations, along with attenuating mutations E580G and S660A, in PB1 of B/Bris (B/Bris PB1att 4M). Serial passages of the B/Bris PB1att 4M vaccine candidate in eggs and MDCK indicated high stability. In silico structural analysis revealed a potential interaction between amino acids at positions 48 and 391. In mice, B/Bris PB1att 4M was safe and provided complete protection against homologous challenge. These results confirm the compensatory effect of mutation E48K to stabilize the K391E mutation, resulting in a safer, yet still protective, IBV LAIV vaccine.

Rampant False Detection of Adaptive Phenotypic Optimization by ParTI-Based Pareto Front Inference

Organisms face tradeoffs in performing multiple tasks. Identifying the optimal phenotypes maximizing the organismal fitness (or Pareto front) and inferring the relevant tasks allow testing phenotypic adaptations and help delineate evolutionary constraints, tradeoffs, and critical fitness components, so are of broad interest. It has been proposed that Pareto fronts can be identified from high-dimensional phenotypic data, including molecular phenotypes such as gene expression levels, by fitting polytopes (lines, triangles, tetrahedrons, and so on), and a program named ParTI was recently introduced for this purpose. ParTI has identified Pareto fronts and inferred phenotypes best for individual tasks (or archetypes) from numerous data sets such as the beak morphologies of Darwin’s finches and mRNA concentrations in human tumors, implying evolutionary optimizations of the involved traits. Nevertheless, the reliabilities of these findings are unknown. Using real and simulated data that lack evolutionary optimization, we here report extremely high false-positive rates of ParTI. The errors arise from phylogenetic relationships or population structures of the organisms analyzed and the flexibility of data analysis in ParTI that is equivalent to p-hacking. Because these problems are virtually universal, our findings cast doubt on almost all ParTI-based results and suggest that reliably identifying Pareto fronts and archetypes from high-dimensional phenotypic data are currently generally difficult.

Model-guided development of an evolutionarily stable yeast chassis

First-principle metabolic modelling holds potential for designing microbial chassis that are resilient against phenotype reversal due to adaptive mutations. Yet, the theory of model-based chassis design has rarely been put to rigorous experimental test. Here, we report the development of Saccharomyces cerevisiae chassis strains for dicarboxylic acid production using genome-scale metabolic modelling. The chassis strains, albeit geared for higher flux towards succinate, fumarate and malate, do not appreciably secrete these metabolites. As predicted by the model, introducing product-specific TCA cycle disruptions resulted in the secretion of the corresponding acid. Adaptive laboratory evolution further improved production of succinate and fumarate, demonstrating the evolutionary robustness of the engineered cells. In the case of malate, multi-omics analysis revealed a flux bypass at peroxisomal malate dehydrogenase that was missing in the yeast metabolic model. In all three cases, flux balance analysis integrating transcriptomics, proteomics and metabolomics data confirmed the flux re-routing predicted by the model. Taken together, our modelling and experimental results have implications for the computer-aided design of microbial cell factories.

Genetic interaction network has a very limited impact on the evolutionary trajectories in continuous culture-grown populations of yeast

Background

The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes.

Results

Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation.

Conclusions

Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01830-9.

Dynamics of bacterial adaptation

Determining pattern in the dynamics of population evolution is a long-standing focus of evolutionary biology. Complementing the study of natural populations, microbial laboratory evolution experiments have become an important tool for addressing these dynamics because they allow detailed and replicated analysis of evolution in response to controlled environmental and genetic conditions. Key findings include a tendency for smoothly declining rates of adaptation during selection in constant environments, at least in part a reflection of antagonism between accumulating beneficial mutations, and a large number of beneficial mutations available to replicate populations leading to significant, but relatively low genetic parallelism, even as phenotypic characteristics show high similarity. Together, there is a picture of adaptation as a process with a varied and largely unpredictable genetic basis leading to much more similar phenotypic outcomes. Increasing sophistication of sequencing and genetic tools will allow insight into mechanisms behind these and other patterns.

Gene Loss Predictably Drives Evolutionary Adaptation

Loss of gene function is common throughout evolution, even though it often leads to reduced fitness. In this study, we systematically evaluated how an organism adapts after deleting genes that are important for growth under oxidative stress. By evolving, sequencing, and phenotyping over 200 yeast lineages, we found that gene loss can enhance an organism’s capacity to evolve and adapt. Although gene loss often led to an immediate decrease in fitness, many mutants rapidly acquired suppressor mutations that restored fitness. Depending on the strain’s genotype, some ultimately even attained higher fitness levels than similarly adapted wild-type cells. Further, cells with deletions in different modules of the genetic network followed distinct and predictable mutational trajectories. Finally, losing highly connected genes increased evolvability by facilitating the emergence of a more diverse array of phenotypes after adaptation. Together, our findings show that loss of specific parts of a genetic network can facilitate adaptation by opening alternative evolutionary paths.

Historical Contingency Causes Divergence in Adaptive Expression of the lac Operon

Populations of Escherichia coli selected in constant and fluctuating environments containing lactose often adapt by substituting mutations in the lacI repressor that cause constitutive expression of the lac operon. These mutations occur at a high rate and provide a significant benefit. Despite this, eight of 24 populations evolved for 8,000 generations in environments containing lactose contained no detectable repressor mutations. We report here on the basis of this observation. We find that, given relevant mutation rates, repressor mutations are expected to have fixed in all evolved populations if they had maintained the same fitness effect they confer when introduced to the ancestor. In fact, reconstruction experiments demonstrate that repressor mutations have become neutral or deleterious in those populations in which they were not detectable. Populations not fixing repressor mutations nevertheless reached the same fitness as those that did fix them, indicating that they followed an alternative evolutionary path that made redundant the potential benefit of the repressor mutation, but involved unique mutations of equivalent benefit. We identify a mutation occurring in the promoter region of the uspB gene as a candidate for influencing the selective choice between these paths. Our results detail an example of historical contingency leading to divergent evolutionary outcomes.

Suboptimal Global Transcriptional Response Increases the Harmful Effects of Loss-of-Function Mutations

The fitness impact of loss-of-function mutations is generally assumed to reflect the loss of specific molecular functions associated with the perturbed gene. Here, we propose that rewiring of the transcriptome upon deleterious gene inactivation is frequently nonspecific and mimics stereotypic responses to external environmental change. Consequently, transcriptional response to gene deletion could be suboptimal and incur an extra fitness cost. Analysis of the transcriptomes of ∼1,500 single-gene deletion Saccharomyces cerevisiae strains supported this scenario. First, most transcriptomic changes are not specific to the deleted gene but are rather triggered by perturbations in functionally diverse genes. Second, gene deletions that alter the expression of dosage-sensitive genes are especially harmful. Third, by elevating the expression level of downregulated genes, we could experimentally mitigate the fitness defect of gene deletions. Our work shows that rewiring of genomic expression upon gene inactivation shapes the harmful effects of mutations.

Evidence for Horizontal and Vertical Transmission of Mtr-Mediated Extracellular Electron Transfer among the Bacteria

Some bacteria and archaea have evolved the means to use extracellular electron donors and acceptors for energy metabolism, a phenomenon broadly known as extracellular electron transfer (EET). One such EET mechanism is the transmembrane electron conduit MtrCAB, which has been shown to transfer electrons derived from metabolic substrates to electron acceptors, like Fe(III) and Mn(IV) oxides, outside the cell. Although most studies of MtrCAB-mediated EET have been conducted in Shewanella oneidensis MR-1, recent investigations in Vibrio and Aeromonas species have revealed that the electron-donating proteins that support MtrCAB in Shewanella are not as representative as previously thought. This begs the question of how widespread the capacity for MtrCAB-mediated EET is, the changes it has accrued in different lineages, and where these lineages persist today. Here, we employed a phylogenetic and comparative genomics approach to identify the MtrCAB system across all domains of life. We found mtrCAB in the genomes of numerous diverse Bacteria from a wide range of environments, and the patterns therein strongly suggest that mtrCAB was distributed through both horizontal and subsequent vertical transmission, and with some cases indicating downstream modular diversification of both its core and accessory components. Our data point to an emerging evolutionary story about metal-oxidizing and -reducing metabolism, demonstrates that this capacity for EET has broad relevance to a diversity of taxa and the biogeochemical cycles they drive, and lays the foundation for further studies to shed light on how this mechanism may have coevolved with Earth's redox landscape. IMPORTANCE While many metabolisms make use of soluble, cell-permeable substrates like oxygen or hydrogen, there are other substrates, like iron or manganese, that cannot be brought into the cell. Some bacteria and archaea have evolved the means to directly "plug in" to such environmental electron reservoirs in a process known as extracellular electron transfer (EET), making them powerful agents of biogeochemical change and promising vehicles for bioremediation and alternative energy. Yet the diversity, distribution, and evolution of EET mechanisms are poorly constrained. Here, we present findings showing that the genes encoding one such EET system (mtrCAB) are present in a broad diversity of bacteria found in a wide range of environments, emphasizing the ubiquity and potential impact of EET in our biosphere. Our results suggest that these genes have been disseminated largely through horizontal transfer, and the changes they have accrued in these lineages potentially reflect adaptations to changing environments.

High-throughput synthetic rescue for exhaustive characterization of suppressor mutations in human genes

Inherited or acquired mutations can lead to pathological outcomes. However, in a process defined as synthetic rescue, phenotypic outcome created by primary mutation is alleviated by suppressor mutations. An exhaustive characterization of these mutations in humans is extremely valuable to better comprehend why patients carrying the same detrimental mutation exhibit different pathological outcomes or different responses to treatment. Here, we first review all known suppressor mutations' mechanisms characterized by genetic screens on model species like yeast or flies. However, human suppressor mutations are scarce, despite some being discovered based on orthologue genes. Because of recent advances in high-throughput screening, developing an inventory of human suppressor mutations for pathological processes seems achievable. In addition, we review several screening methods for suppressor mutations in cultured human cells through knock-out, knock-down or random mutagenesis screens on large scale. We provide examples of studies published over the past years that opened new therapeutic avenues, particularly in oncology.

The Whi2p-Psr1p/Psr2p complex regulates interference competition and expansion of cells with competitive advantage in yeast colonies

Yeast form complex highly organized colonies in which cells undergo spatiotemporal phenotypic differentiation in response to local gradients of nutrients, metabolites, and specific signaling molecules. Colony fitness depends on cell interactions, cooperation, and the division of labor between differentiated cell subpopulations. Here, we describe the regulation and dynamics of the expansion of papillae that arise during colony aging, which consist of cells that overcome colony regulatory rules and disrupt the synchronized colony structure. We show that papillae specifically expand within the U cell subpopulation in differentiated colonies. Papillae emerge more frequently in some strains than in others. Genomic analyses further revealed that the Whi2p-Psr1p/Psr2p complex (WPPC) plays a key role in papillae expansion. We show that cells lacking a functional WPPC have a sizable interaction-specific fitness advantage attributable to production of and resistance to a diffusible compound that inhibits growth of other cells. Competitive superiority and high relative fitness of whi2 and psr1psr2 strains are particularly pronounced in dense spatially structured colonies and are independent of TORC1 and Msn2p/Msn4p regulators previously associated with the WPPC function. The WPPC function, described here, might be a regulatory mechanism that balances cell competition and cooperation in dense yeast populations and, thus, contributes to cell synchronization, pattern formation, and the expansion of cells with a competitive fitness advantage.

Evolutionary Repair Experiments as a Window to the Molecular Diversity of Life

Comparative genomics reveals an unexpected diversity in the molecular mechanisms underlying conserved cellular functions, such as DNA replication and cytokinesis. However, the genetic bases and evolutionary processes underlying this 'molecular diversity' remain to be explained. Here, we review a tool to generate alternative mechanisms for conserved cellular functions and test hypotheses concerning the generation of molecular diversity - evolutionary repair experiments, in which laboratory microbial populations adapt in response to a genetic perturbation. We summarize the insights gained from evolutionary repair experiments, the spectrum and dynamics of compensatory mutations, and the alternative molecular mechanisms used to repair perturbed cellular functions. We relate these experiments to the modifications of conserved functions that have occurred outside the laboratory. We end by proposing strategies to improve evolutionary repair experiments as a tool to explore the molecular diversity of life.

Predicting Evolution Using Regulatory Architecture

The limits of evolution have long fascinated biologists. However, the causes of evolutionary constraint have remained elusive due to a poor mechanistic understanding of studied phenotypes. Recently, a range of innovative approaches have leveraged mechanistic information on regulatory networks and cellular biology. These methods combine systems biology models with population and single-cell quantification and with new genetic tools, and they have been applied to a range of complex cellular functions and engineered networks. In this article, we review these developments, which are revealing the mechanistic causes of epistasis at different levels of biological organization-in molecular recognition, within a single regulatory network, and between different networks-providing first indications of predictable features of evolutionary constraint.

The evolutionary plasticity of chromosome metabolism allows adaptation to constitutive DNA replication stress

Many biological features are conserved and thus considered to be resistant to evolutionary change. While rapid genetic adaptation following the removal of conserved genes has been observed, we often lack a mechanistic understanding of how adaptation happens. We used the budding yeast, Saccharomyces cerevisiae, to investigate the evolutionary plasticity of chromosome metabolism, a network of evolutionary conserved modules. We experimentally evolved cells constitutively experiencing DNA replication stress caused by the absence of Ctf4, a protein that coordinates the enzymatic activities at replication forks. Parallel populations adapted to replication stress, over 1000 generations, by acquiring multiple, concerted mutations. These mutations altered conserved features of two chromosome metabolism modules, DNA replication and sister chromatid cohesion, and inactivated a third, the DNA damage checkpoint. The selected mutations define a functionally reproducible evolutionary trajectory. We suggest that the evolutionary plasticity of chromosome metabolism has implications for genome evolution in natural populations and cancer.

WHI-2 Regulates Intercellular Communication via a MAP Kinase Signaling Complex

The formation of the fungal mycelial network is facilitated by somatic cell fusion of germinating asexual spores (or germlings). Neurospora crassa germlings in close proximity display chemotropic growth that is dependent upon an intracellular network of mitogen-activated protein kinase (MAPK) signaling cascades. Approximately 80 genes involved in intercellular communication and fusion have been identified, including three mutants with similar morphological phenotypes: Δwhi-2, Δcsp-6, and Δamph-1. Here we show that WHI-2 localizes to the cell periphery and regulates endocytosis, mitochondrial organization, sporulation, and cell fusion. WHI-2 was required to transduce signals through a conserved MAPK pathway (NRC-1/MEK-2/MAK-2) and target transcription factors (PP-1/ADV-1). The amph-1 locus encodes a Bin/Amphiphysin/Rvs domain-containing protein and mis-expression of whi-2 compensated for the cell fusion and endocytosis deficiencies of a Δamph-1 mutant. The csp-6 locus encodes a haloacid dehalogenase phosphatase whose activity was essential for cell fusion. Although fusion-deficient with themselves, cells that lacked whi-2, csp-6, or amph-1 showed a low frequency of chemotropic interactions with wild type cells. We hypothesize that WHI-2 could be important for signal perception during chemotropic interactions via a role in endocytosis.