Minimum information and guidelines for reporting a multiplexed assay of variant effect.

Multiplexed assays of variant effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines have led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Minimum information and guidelines for reporting a Multiplexed Assay of Variant Effect.

Multiplexed Assays of Variant Effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines has led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Microbial experimental evolution in a massively multiplexed and high-throughput era.

Experimental evolution with microbial model systems has transformed our understanding of the basic rules underlying ecology and evolution. Experiments leveraging evolution as a central feature put evolutionary theories to the test, and modern sequencing and engineering tools then characterized the molecular basis of adaptation. As theory and experimentations refined our understanding of evolution, a need to increase throughput and experimental complexity has emerged. Here, we summarize recent technologies that have made high-throughput experiments practical and highlight studies that have capitalized on these tools, defining an exciting new era in microbial experimental evolution. Multiple research directions previously limited by experimental scale are now accessible for study and we believe applying evolutionary lessons from in vitro studies onto these applied settings has the potential for major innovations and discoveries across ecology and medicine.

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.

Bulk-Fitness Measurements Using Barcode Sequencing Analysis in Yeast.

The use of DNA barcodes for determining changes in genotype frequencies has been instrumental to increase the scale at which we can phenotype strain libraries by using next-generation sequencing technologies. Here, we describe the determination of strain fitness for thousands of yeast strains simultaneously in a single assay using recent innovations that increase the precision of these measurements, such as the inclusion of unique-molecular identifiers (UMIs) and purification by solid-phase reverse immobilization (SPRI) beads.

Barcoded bulk QTL mapping reveals highly polygenic and epistatic architecture of complex traits in yeast.

Mapping the genetic basis of complex traits is critical to uncovering the biological mechanisms that underlie disease and other phenotypes. Genome-wide association studies (GWAS) in humans and quantitative trait locus (QTL) mapping in model organisms can now explain much of the observed heritability in many traits, allowing us to predict phenotype from genotype. However, constraints on power due to statistical confounders in large GWAS and smaller sample sizes in QTL studies still limit our ability to resolve numerous small-effect variants, map them to causal genes, identify pleiotropic effects across multiple traits, and infer non-additive interactions between loci (epistasis). Here, we introduce barcoded bulk quantitative trait locus (BB-QTL) mapping, which allows us to construct, genotype, and phenotype 100,000 offspring of a budding yeast cross, two orders of magnitude larger than the previous state of the art. We use this panel to map the genetic basis of eighteen complex traits, finding that the genetic architecture of these traits involves hundreds of small-effect loci densely spaced throughout the genome, many with widespread pleiotropic effects across multiple traits. Epistasis plays a central role, with thousands of interactions that provide insight into genetic networks. By dramatically increasing sample size, BB-QTL mapping demonstrates the potential of natural variants in high-powered QTL studies to reveal the highly polygenic, pleiotropic, and epistatic architecture of complex traits.

The genetic basis of differential autodiploidization in evolving yeast populations.

Spontaneous whole-genome duplication, or autodiploidization, is a common route to adaptation in experimental evolution of haploid budding yeast populations. The rate at which autodiploids fix in these populations appears to vary across strain backgrounds, but the genetic basis of these differences remains poorly characterized. Here, we show that the frequency of autodiploidization differs dramatically between two closely related laboratory strains of Saccharomyces cerevisiae, BY4741 and W303. To investigate the genetic basis of this difference, we crossed these strains to generate hundreds of unique F1 segregants and tested the tendency of each segregant to autodiplodize across hundreds of generations of laboratory evolution. We find that variants in the SSD1 gene are the primary genetic determinant of differences in autodiploidization. We then used multiple laboratory and wild strains of S. cerevisiae to show that clonal populations of strains with a functional copy of SSD1 autodiploidize more frequently in evolution experiments, while knocking out this gene or replacing it with the W303 allele reduces autodiploidization propensity across all genetic backgrounds tested. These results suggest a potential strategy for modifying rates of spontaneous whole-genome duplications in laboratory evolution experiments in haploid budding yeast. They may also have relevance to other settings in which eukaryotic genome stability plays an important role, such as biomanufacturing and the treatment of pathogenic fungal diseases and cancers.

Dynamics and variability in the pleiotropic effects of adaptation in laboratory budding yeast populations.

Evolutionary adaptation to a constant environment is driven by the accumulation of mutations which can have a range of unrealized pleiotropic effects in other environments. These pleiotropic consequences of adaptation can influence the emergence of specialists or generalists, and are critical for evolution in temporally or spatially fluctuating environments. While many experiments have examined the pleiotropic effects of adaptation at a snapshot in time, very few have observed the dynamics by which these effects emerge and evolve. Here, we propagated hundreds of diploid and haploid laboratory budding yeast populations in each of three environments, and then assayed their fitness in multiple environments over 1000 generations of evolution. We find that replicate populations evolved in the same condition share common patterns of pleiotropic effects across other environments, which emerge within the first several hundred generations of evolution. However, we also find dynamic and environment-specific variability within these trends: variability in pleiotropic effects tends to increase over time, with the extent of variability depending on the evolution environment. These results suggest shifting and overlapping contributions of chance and contingency to the pleiotropic effects of adaptation, which could influence evolutionary trajectories in complex environments that fluctuate across space and time.

Phenotypic and molecular evolution across 10,000 generations in laboratory budding yeast populations.

Laboratory experimental evolution provides a window into the details of the evolutionary process. To investigate the consequences of long-term adaptation, we evolved 205 Saccharomyces cerevisiae populations (124 haploid and 81 diploid) for ~10,000,000 generations in three environments. We measured the dynamics of fitness changes over time, finding repeatable patterns of declining adaptability. Sequencing revealed that this phenotypic adaptation is coupled with a steady accumulation of mutations, widespread genetic parallelism, and historical contingency. In contrast to long-term evolution in E. coli, we do not observe long-term coexistence or populations with highly elevated mutation rates. We find that evolution in diploid populations involves both fixation of heterozygous mutations and frequent loss-of-heterozygosity events. Together, these results help distinguish aspects of evolutionary dynamics that are likely to be general features of adaptation across many systems from those that are specific to individual organisms and environmental conditions.

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis.

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

Microwave-Generated Steam Decontamination of N95 Respirators Utilizing Universally Accessible Materials.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has caused a severe, international shortage of N95 respirators, which are essential to protect health care providers from infection. Given the contemporary limitations of the supply chain, it is imperative to identify effective means of decontaminating, reusing, and thereby conserving N95 respirator stockpiles. To be effective, decontamination must result in sterilization of the N95 respirator without impairment of respirator filtration or user fit. Although numerous methods of N95 decontamination exist, none are universally accessible. In this work, we describe a microwave-generated steam decontamination protocol for N95 respirators for use in health care systems of all sizes, geographies, and means. Using widely available glass containers, mesh from commercial produce bags, a rubber band, and a 1,100-W commercially available microwave, we constructed an effective, standardized, and reproducible means of decontaminating N95 respirators. Employing this methodology against MS2 phage, a highly conservative surrogate for SARS-CoV-2 contamination, we report an average 6-log10 plaque-forming unit (PFU) (99.9999%) and a minimum 5-log10 PFU (99.999%) reduction after a single 3-min microwave treatment. Notably, quantified respirator fit and function were preserved, even after 20 sequential cycles of microwave steam decontamination. This method provides a valuable means of effective decontamination and reuse of N95 respirators by frontline providers facing urgent need.IMPORTANCE Due to the rapid spread of coronavirus disease 2019 (COVID-19), there is an increasing shortage of protective gear necessary to keep health care providers safe from infection. As of 9 April 2020, the CDC reported 9,282 cumulative cases of COVID-19 among U.S. health care workers (CDC COVID-19 Response Team, MMWR Morb Mortal Wkly Rep 69:477-481, 2020, https://doi.org/10.15585/mmwr.mm6915e6). N95 respirators are recommended by the CDC as the ideal method of protection from COVID-19. Although N95 respirators are traditionally single use, the shortages have necessitated the need for reuse. Effective methods of N95 decontamination that do not affect the fit or filtration ability of N95 respirators are essential. Numerous methods of N95 decontamination exist; however, none are universally accessible. In this study, we describe an effective, standardized, and reproducible means of decontaminating N95 respirators using widely available materials. The N95 decontamination method described in this work will provide a valuable resource for hospitals, health care centers, and outpatient practices that are experiencing increasing shortages of N95 respirators due to the COVID-19 pandemic.

Chance and necessity in the pleiotropic consequences of adaptation for budding yeast.

Mutations that a population accumulates during evolution in one 'home' environment may cause fitness gains or losses in other environments. Such pleiotropic fitness effects determine the evolutionary fate of the population in variable environments and can lead to ecological specialization. It is unclear how the pleiotropic outcomes of evolution are shaped by the intrinsic randomness of the evolutionary process and by the deterministic variation in selection pressures across environments. Here, to address this question, we evolved 20 replicate populations of the yeast Saccharomyces cerevisiae in 11 laboratory environments and measured their fitness across multiple conditions. We found that evolution led to diverse pleiotropic fitness gains and losses, driven by multiple types of mutations. Approximately 60% of this variation is explained by the home environment of a clone and the most common parallel genetic changes, whereas about 40% is attributed to the stochastic accumulation of mutations whose pleiotropic effects are unpredictable. Although populations are typically specialized to their home environment, generalists also evolved in almost all of the conditions. Our results suggest that the mutations that accumulate during evolution incur a variety of pleiotropic costs and benefits with different probabilities. Thus, whether a population evolves towards a specialist or a generalist phenotype is heavily influenced by chance.

Experimental evolution for niche breadth in bacteriophage T4 highlights the importance of structural genes.

Ecologists have long studied the evolution of niche breadth, including how variability in environments can drive the evolution of specialism and generalism. This concept is of particular interest in viruses, where niche breadth evolution may explain viral disease emergence, or underlie the potential for therapeutic measures like phage therapy. Despite the significance and potential applications of virus-host interactions, the genetic determinants of niche breadth evolution remain underexplored in many bacteriophages. In this study, we present the results of an evolution experiment with a model bacteriophage system, Escherichia virus T4, in several host environments: exposure to Escherichia coli C, exposure to E. coli K-12, and exposure to both E. coli C and E. coli K-12. This experimental framework allowed us to investigate the phenotypic and molecular manifestations of niche breadth evolution. First, we show that selection on different hosts led to measurable changes in phage productivity in all experimental populations. Second, whole-genome sequencing of experimental populations revealed signatures of selection. Finally, clear and consistent patterns emerged across the host environments, especially the presence of new mutations in phage structural genes-genes encoding proteins that provide morphological and biophysical integrity to a virus. A comparison of mutations found across functional gene categories revealed that structural genes acquired significantly more mutations than other categories. Our findings suggest that structural genes are central determinants in bacteriophage niche breadth.

High-resolution lineage tracking reveals travelling wave of adaptation in laboratory yeast.

In rapidly adapting asexual populations, including many microbial pathogens and viruses, numerous mutant lineages often compete for dominance within the population1-5. These complex evolutionary dynamics determine the outcomes of adaptation, but have been difficult to observe directly. Previous studies have used whole-genome sequencing to follow molecular adaptation6-10; however, these methods have limited resolution in microbial populations. Here we introduce a renewable barcoding system to observe evolutionary dynamics at high resolution in laboratory budding yeast. We find nested patterns of interference and hitchhiking even at low frequencies. These events are driven by the continuous appearance of new mutations that modify the fates of existing lineages before they reach substantial frequencies. We observe how the distribution of fitness within the population changes over time, and find a travelling wave of adaptation that has been predicted by theory11-17. We show that clonal competition creates a dynamical 'rich-get-richer' effect: fitness advantages that are acquired early in evolution drive clonal expansions, which increase the chances of acquiring future mutations. However, less-fit lineages also routinely leapfrog over strains of higher fitness. Our results demonstrate that this combination of factors, which is not accounted for in existing models of evolutionary dynamics, is critical in determining the rate, predictability and molecular basis of adaptation.

Proteome-wide signatures of function in highly diverged intrinsically disordered regions.

Intrinsically disordered regions make up a large part of the proteome, but the sequence-to-function relationship in these regions is poorly understood, in part because the primary amino acid sequences of these regions are poorly conserved in alignments. Here we use an evolutionary approach to detect molecular features that are preserved in the amino acid sequences of orthologous intrinsically disordered regions. We find that most disordered regions contain multiple molecular features that are preserved, and we define these as 'evolutionary signatures' of disordered regions. We demonstrate that intrinsically disordered regions with similar evolutionary signatures can rescue function in vivo, and that groups of intrinsically disordered regions with similar evolutionary signatures are strongly enriched for functional annotations and phenotypes. We propose that evolutionary signatures can be used to predict function for many disordered regions from their amino acid sequences.

Evolthon: A community endeavor to evolve lab evolution.

In experimental evolution, scientists evolve organisms in the lab, typically by challenging them to new environmental conditions. How best to evolve a desired trait? Should the challenge be applied abruptly, gradually, periodically, sporadically? Should one apply chemical mutagenesis, and do strains with high innate mutation rate evolve faster? What are ideal population sizes of evolving populations? There are endless strategies, beyond those that can be exposed by individual labs. We therefore arranged a community challenge, Evolthon, in which students and scientists from different labs were asked to evolve Escherichia coli or Saccharomyces cerevisiae for an abiotic stress-low temperature. About 30 participants from around the world explored diverse environmental and genetic regimes of evolution. After a period of evolution in each lab, all strains of each species were competed with one another. In yeast, the most successful strategies were those that used mating, underscoring the importance of sex in evolution. In bacteria, the fittest strain used a strategy based on exploration of different mutation rates. Different strategies displayed variable levels of performance and stability across additional challenges and conditions. This study therefore uncovers principles of effective experimental evolutionary regimens and might prove useful also for biotechnological developments of new strains and for understanding natural strategies in evolutionary arms races between species. Evolthon constitutes a model for community-based scientific exploration that encourages creativity and cooperation.

Modular epistasis and the compensatory evolution of gene deletion mutants.

Screens for epistatic interactions have long been used to characterize functional relationships corresponding to protein complexes, metabolic pathways, and other functional modules. Although epistasis between adaptive mutations is also common in laboratory evolution experiments, the functional basis for these interactions is less well characterized. Here, we quantify the extent to which gene function (as determined by a genome-wide screen for epistasis among deletion mutants) influences the rate and genetic basis of compensatory adaptation in a set of 37 gene deletion mutants nested within 16 functional modules. We find that functional module has predictive power: mutants with deletions in the same module tend to adapt more similarly, on average, than those with deletions in different modules. At the same time, initial fitness also plays a role: independent of the specific functional modules involved, adaptive mutations tend to be systematically more beneficial in less-fit genetic backgrounds, consistent with a general pattern of diminishing returns epistasis. We measured epistatic interactions between initial gene deletion mutations and the mutations that accumulate during compensatory adaptation and found a general trend towards positive epistasis (i.e. mutations tend to be more beneficial in the background in which they arose). In two functional modules, epistatic interactions between the initial gene deletions and the mutations in their descendant lines caused evolutionary entrenchment, indicating an intimate functional relationship. Our results suggest that genotypes with similar epistatic interactions with gene deletion mutations will also have similar epistatic interactions with adaptive mutations, meaning that genome scale maps of epistasis between gene deletion mutations can be predictive of evolutionary dynamics.

Experimental Studies of Evolutionary Dynamics in Microbes.

Evolutionary dynamics in laboratory microbial evolution experiments can be surprisingly complex. In the past two decades, observations of these dynamics have challenged simple models of adaptation and have shown that clonal interference, hitchhiking, ecological diversification, and contingency are widespread. In recent years, advances in high-throughput strain maintenance and phenotypic assays, the dramatically reduced cost of genome sequencing, and emerging methods for lineage barcoding have made it possible to observe evolutionary dynamics at unprecedented resolution. These new methods can now begin to provide detailed measurements of key aspects of fitness landscapes and of evolutionary outcomes across a range of systems. These measurements can highlight challenges to existing theoretical models and guide new theoretical work towards the complications that are most widely important.

Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments.

Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.

Selection maintains signaling function of a highly diverged intrinsically disordered region.

Intrinsically disordered regions (IDRs) are characterized by their lack of stable secondary or tertiary structure and comprise a large part of the eukaryotic proteome. Although these regions play a variety of signaling and regulatory roles, they appear to be rapidly evolving at the primary sequence level. To understand the functional implications of this rapid evolution, we focused on a highly diverged IDR in Saccharomyces cerevisiae that is involved in regulating multiple conserved MAPK pathways. We hypothesized that under stabilizing selection, the functional output of orthologous IDRs could be maintained, such that diverse genotypes could lead to similar function and fitness. Consistent with the stabilizing selection hypothesis, we find that diverged, orthologous IDRs can mostly recapitulate wild-type function and fitness in S. cerevisiae We also find that the electrostatic charge of the IDR is correlated with signaling output and, using phylogenetic comparative methods, find evidence for selection maintaining this quantitative molecular trait despite underlying genotypic divergence.