The Dynamics of Adaptation to Stress from Standing Genetic Variation and de novo Mutations.

Adaptation from standing genetic variation is an important process underlying evolution in natural populations, but we rarely get the opportunity to observe the dynamics of fitness and genomic changes in real time. Here, we used experimental evolution and Pool-Seq to track the phenotypic and genomic changes of genetically diverse asexual populations of the yeast Saccharomyces cerevisiae in four environments with different fitness costs. We found that populations rapidly and in parallel increased in fitness in stressful environments. In contrast, allele frequencies showed a range of trajectories, with some populations fixing all their ancestral variation in <30 generations and others maintaining diversity across hundreds of generations. We detected parallelism at the genomic level (involving genes, pathways, and aneuploidies) within and between environments, with idiosyncratic changes recurring in the environments with higher stress. In particular, we observed a tendency of becoming haploid-like in one environment, whereas the populations of another environment showed low overall parallelism driven by standing genetic variation despite high selective pressure. This work highlights the interplay between standing genetic variation and the influx of de novo mutations in populations adapting to a range of selective pressures with different underlying trait architectures, advancing our understanding of the constraints and drivers of adaptation.

Shepherd: accurate clustering for correcting DNA barcode errors.

MOTIVATION: DNA barcodes are short, random nucleotide sequences introduced into cell populations to track the relative counts of hundreds of thousands of individual lineages over time. Lineage tracking is widely applied, e.g. to understand evolutionary dynamics in microbial populations and the progression of breast cancer in humans. Barcode sequences are unknown upon insertion and must be identified using next-generation sequencing technology, which is error prone. In this study, we frame the barcode error correction task as a clustering problem with the aim to identify true barcode sequences from noisy sequencing data. We present Shepherd, a novel clustering method that is based on an indexing system of barcode sequences using k-mers, and a Bayesian statistical test incorporating a substitution error rate to distinguish true from error sequences. RESULTS: When benchmarking with synthetic data, Shepherd provides barcode count estimates that are significantly more accurate than state-of-the-art methods, producing 10-150 times fewer spurious lineages. For empirical data, Shepherd produces results that are consistent with the improvements seen on synthetic data. These improvements enable higher resolution lineage tracking and more accurate estimates of biologically relevant quantities, e.g. the detection of small effect mutations. AVAILABILITY AND IMPLEMENTATION: A Python implementation of Shepherd is freely available at: https://www.github.com/Nik-Tavakolian/Shepherd. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Forest and woodland replacement patterns following drought-related mortality.

Forest vulnerability to drought is expected to increase under anthropogenic climate change, and drought-induced mortality and community dynamics following drought have major ecological and societal impacts. Here, we show that tree mortality concomitant with drought has led to short-term (mean 5 y, range 1 to 23 y after mortality) vegetation-type conversion in multiple biomes across the world (131 sites). Self-replacement of the dominant tree species was only prevalent in 21% of the examined cases and forests and woodlands shifted to nonwoody vegetation in 10% of them. The ultimate temporal persistence of such changes remains unknown but, given the key role of biological legacies in long-term ecological succession, this emerging picture of postdrought ecological trajectories highlights the potential for major ecosystem reorganization in the coming decades. Community changes were less pronounced under wetter postmortality conditions. Replacement was also influenced by management intensity, and postdrought shrub dominance was higher when pathogens acted as codrivers of tree mortality. Early change in community composition indicates that forests dominated by mesic species generally shifted toward more xeric communities, with replacing tree and shrub species exhibiting drier bioclimatic optima and distribution ranges. However, shifts toward more mesic communities also occurred and multiple pathways of forest replacement were observed for some species. Drought characteristics, species-specific environmental preferences, plant traits, and ecosystem legacies govern postdrought species turnover and subsequent ecological trajectories, with potential far-reaching implications for forest biodiversity and ecosystem services.

Experimental Resurrection of Ancestral Mammalian CPEB3 Ribozymes Reveals Deep Functional Conservation.

Self-cleaving ribozymes are genetic elements found in all domains of life, but their evolution remains poorly understood. A ribozyme located in the second intron of the cytoplasmic polyadenylation binding protein 3 gene (CPEB3) shows high sequence conservation in mammals, but little is known about the functional conservation of self-cleaving ribozyme activity across the mammalian tree of life or during the course of mammalian evolution. Here, we use a phylogenetic approach to design a mutational library and a deep sequencing assay to evaluate the in vitro self-cleavage activity of numerous extant and resurrected CPEB3 ribozymes that span over 100 My of mammalian evolution. We found that the predicted sequence at the divergence of placentals and marsupials is highly active, and this activity has been conserved in most lineages. A reduction in ribozyme activity appears to have occurred multiple different times throughout the mammalian tree of life. The in vitro activity data allow an evaluation of the predicted mutational pathways leading to extant ribozyme as well as the mutational landscape surrounding these ribozymes. The results demonstrate that in addition to sequence conservation, the self-cleavage activity of the CPEB3 ribozyme has persisted over millions of years of mammalian evolution.

Phased nucleotide inserts for sequencing low-diversity RNA samples from in vitro selection experiments.

In vitro selection combined with high-throughput sequencing is a powerful experimental approach with broad application in the engineering and characterization of RNA molecules. Diverse pools of starting sequences used for selection are often flanked by fixed sequences used as primer binding sites. These low diversity regions often lead to data loss from complications with Illumina image processing algorithms. A common method to alleviate this problem is the addition of fragmented bacteriophage PhiX genome, which improves sequence quality but sacrifices a portion of usable sequencing reads. An alternative approach is to insert nucleotides of variable length and composition ("phased inserts") at the beginning of each molecule when adding sequencing adaptors. This approach preserves read depth but reduces the length of each read. Here, we test the ability of phased inserts to replace PhiX in a low-diversity sample generated for a high-throughput sequencing based ribozyme activity screen. We designed a pool of 4096 RNA sequence variants of the self-cleaving twister ribozyme from Oryza sativa For each unique sequence, we determined the fraction of ribozyme cleaved during in vitro transcription via deep sequencing on an Illumina MiSeq. We found that libraries with the phased inserts produced high-quality sequence data without the addition of PhiX. We found good agreement between previously published data on twister ribozyme variants and our data produced with phased inserts even when PhiX was omitted. We conclude that phased inserts can be implemented following in vitro selection experiments to reduce or eliminate the use of PhiX and maximize read depth.

Saccharomyces yeast hybrids on the rise.

Saccharomyces hybrid yeasts are receiving increasing attention as a powerful model system to understand adaptation to environmental stress and speciation mechanisms, using experimental evolution and omics techniques. We compiled all genomic resources available from public repositories of the eight recognized Saccharomyces species and their interspecific hybrids. We present the newest numbers on genomes sequenced, assemblies, annotations, and sequencing runs, and an updated species phylogeny using orthogroup inference. While genomic resources are highly skewed towards Saccharomyces cerevisiae, there is a noticeable movement to use wild, recently discovered yeast species in recent years. To illustrate the degree and potential causes of reproductive isolation, we reanalyzed published data on hybrid spore viabilities across the entire genus and tested for the role of genetic, geographic, and ecological divergence within and between species (28 cross types and 371 independent crosses). Hybrid viability generally decreased with parental genetic distance likely due to antirecombination and negative epistasis, but notable exceptions emphasize the importance of strain-specific structural variation and ploidy differences. Surprisingly, the viability of crosses within species varied widely, from near reproductive isolation to near-perfect viability. Geographic and ecological origins of the parents predicted cross viability to an extent, but with certain caveats. Finally, we highlight publication trends in the field and point out areas of special interest, where hybrid yeasts are particularly promising for innovation through research and development, and experimental evolution and fermentation.

Genotype network intersections promote evolutionary innovation.

Evolutionary innovations are qualitatively novel traits that emerge through evolution and increase biodiversity. The genetic mechanisms of innovation remain poorly understood. A systems view of innovation requires the analysis of genotype networks-the vast networks of genetic variants that produce the same phenotype. Innovations can occur at the intersection of two different genotype networks. However, the experimental characterization of genotype networks has been hindered by the vast number of genetic variants that need to be functionally analyzed. Here, we use high-throughput sequencing to study the fitness landscape at the intersection of the genotype networks of two catalytic RNA molecules (ribozymes). We determined the ability of numerous neighboring RNA sequences to catalyze two different chemical reactions, and we use these data as a proxy for a genotype to fitness map where two functions come in close proximity. We find extensive functional overlap, and numerous genotypes can catalyze both functions. We demonstrate through evolutionary simulations that these numerous points of intersection facilitate the discovery of a new function. However, the rate of adaptation of the new function depends upon the local ruggedness around the starting location in the genotype network. As a consequence, one direction of adaptation is more rapid than the other. We find that periods of neutral evolution increase rates of adaptation to the new function by allowing populations to spread out in their genotype network. Our study reveals the properties of a fitness landscape where genotype networks intersect and the consequences for evolutionary innovations. Our results suggest that historic innovations in natural systems may have been facilitated by overlapping genotype networks.

Recombining Your Way Out of Trouble: The Genetic Architecture of Hybrid Fitness under Environmental Stress.

Hybridization between species can either promote or impede adaptation. But we know very little about the genetic basis of hybrid fitness, especially in nondomesticated organisms, and when populations are facing environmental stress. We made genetically variable F2 hybrid populations from two divergent Saccharomyces yeast species. We exposed populations to ten toxins and sequenced the most resilient hybrids on low coverage using ddRADseq to investigate four aspects of their genomes: 1) hybridity, 2) interspecific heterozygosity, 3) epistasis (positive or negative associations between nonhomologous chromosomes), and 4) ploidy. We used linear mixed-effect models and simulations to measure to which extent hybrid genome composition was contingent on the environment. Genomes grown in different environments varied in every aspect of hybridness measured, revealing strong genotype-environment interactions. We also found selection against heterozygosity or directional selection for one of the parental alleles, with larger fitness of genomes carrying more homozygous allelic combinations in an otherwise hybrid genomic background. In addition, individual chromosomes and chromosomal interactions showed significant species biases and pervasive aneuploidies. Against our expectations, we observed multiple beneficial, opposite-species chromosome associations, confirmed by epistasis- and selection-free computer simulations, which is surprising given the large divergence of parental genomes (∼15%). Together, these results suggest that successful, stress-resilient hybrid genomes can be assembled from the best features of both parents without paying high costs of negative epistasis. This illustrates the importance of measuring genetic trait architecture in an environmental context when determining the evolutionary potential of genetically diverse hybrid populations.

Hybridization Outcomes Have Strong Genomic and Environmental Contingencies.

AbstractExtreme F2 phenotypes known as transgressive segregants can cause increased or decreased fitness in hybrids beyond the ranges seen in parental populations. Despite the usefulness of transgression for plant and animal breeding and its potential role in hybrid speciation, the genetic mechanisms and predictors of transgressive segregation remain largely untested. We generated seven hybrid crosses between five widely divergent Saccharomyces yeast species and measured the fitness of the parents and their viable F1 and F2 hybrids in seven stressful environments. We found that on average 16.6% of all replicate F2 hybrids had higher fitness than both parents. Against our predictions, transgression frequency was not a function of parental genetic and phenotypic distances across test environments. Within environments, some relationships were significant, but not in the predicted direction; for example, genetic distance was negatively related to transgression in ethanol and hydrogen peroxide. Significant effects of hybrid cross, test environment, and cross × environment interactions suggest that the amount of transgression produced in a hybrid cross is highly context specific and that outcomes of hybridization differ even among crosses made from the same two parents. If the goal is to reliably predict hybrid fitness and forecast the evolutionary potential of admixed populations, we need more efforts to identify patterns beyond the idiosyncrasies caused by specific genomic or environmental contexts.

Intramolecular phenotypic capacitance in a modular RNA molecule.

Phenotypic capacitance refers to the ability of a genome to accumulate mutations that are conditionally hidden and only reveal phenotype-altering effects after certain environmental or genetic changes. Capacitance has important implications for the evolution of novel forms and functions, but experimentally studied mechanisms behind capacitance are mostly limited to complex, multicomponent systems often involving several interacting protein molecules. Here we demonstrate phenotypic capacitance within a much simpler system, an individual RNA molecule with catalytic activity (ribozyme). This naturally occurring RNA molecule has a modular structure, where a scaffold module acts as an intramolecular chaperone that facilitates folding of a second catalytic module. Previous studies have shown that the scaffold module is not absolutely required for activity, but dramatically decreases the concentration of magnesium ions required for the formation of an active site. Here, we use an experimental perturbation of magnesium ion concentration that disrupts the folding of certain genetic variants of this ribozyme and use in vitro selection followed by deep sequencing to identify genotypes with altered phenotypes (catalytic activity). We identify multiple conditional mutations that alter the wild-type ribozyme phenotype under a stressful environmental condition of low magnesium ion concentration, but preserve the phenotype under more relaxed conditions. This conditional buffering is confined to the scaffold module, but controls the catalytic phenotype, demonstrating how modularity can enable phenotypic capacitance within a single macromolecule. RNA's ancient role in life suggests that phenotypic capacitance may have influenced evolution since life's origins.

Negative Epistasis in Experimental RNA Fitness Landscapes.

Mutations and their effects on fitness are a fundamental component of evolution. The effects of some mutations change in the presence of other mutations, and this is referred to as epistasis. Epistasis can occur between mutations in different genes or within the same gene. A systematic study of epistasis requires the analysis of numerous mutations and their combinations, which has recently become feasible with advancements in DNA synthesis and sequencing. Here we review the mutational effects and epistatic interactions within RNA molecules revealed by several recent high-throughput mutational studies involving two ribozymes studied in vitro, as well as a tRNA and a snoRNA studied in yeast. The data allow an analysis of the distribution of fitness effects of individual mutations as well as combinations of two or more mutations. Two different approaches to measuring epistasis in the data both reveal a predominance of negative epistasis, such that higher combinations of two or more mutations are typically lower in fitness than expected from the effect of each individual mutation. These data are in contrast to past studies of epistasis that used computationally predicted secondary structures of RNA that revealed a predominance of positive epistasis. The RNA data reviewed here are more similar to that found from mutational experiments on individual protein enzymes, suggesting that a common thermodynamic framework may explain negative epistasis between mutations within macromolecules.

RNA sequence to structure analysis from comprehensive pairwise mutagenesis of multiple self-cleaving ribozymes.

Self-cleaving ribozymes are RNA molecules that catalyze the cleavage of their own phosphodiester backbones. These ribozymes are found in all domains of life and are also a tool for biotechnical and synthetic biology applications. Self-cleaving ribozymes are also an important model of sequence-to-function relationships for RNA because their small size simplifies synthesis of genetic variants and self-cleaving activity is an accessible readout of the functional consequence of the mutation. Here, we used a high-throughput experimental approach to determine the relative activity for every possible single and double mutant of five self-cleaving ribozymes. From this data, we comprehensively identified non-additive effects between pairs of mutations (epistasis) for all five ribozymes. We analyzed how changes in activity and trends in epistasis map to the ribozyme structures. The variety of structures studied provided opportunities to observe several examples of common structural elements, and the data was collected under identical experimental conditions to enable direct comparison. Heatmap-based visualization of the data revealed patterns indicating structural features of the ribozymes including paired regions, unpaired loops, non-canonical structures, and tertiary structural contacts. The data also revealed signatures of functionally critical nucleotides involved in catalysis. The results demonstrate that the data sets provide structural information similar to chemical or enzymatic probing experiments, but with additional quantitative functional information. The large-scale data sets can be used for models predicting structure and function and for efforts to engineer self-cleaving ribozymes.

The evolutionary and ecological potential of yeast hybrids.

Recent findings in yeast genetics and genomics have advanced our understanding of the evolutionary potential unlocked by hybridization, especially in the genus Saccharomyces. We now have a clearer picture of the prevalence of yeast hybrids in the environment, their ecological and evolutionary history, and the genetic mechanisms driving (and constraining) their adaptation. Here, we describe how the instability of hybrid genomes determines fitness across large evolutionary scales, highlight new hybrid strain engineering techniques, and review tools for comparative hybrid genome analysis. The recent push to take yeast research back 'into the wild' has resulted in new genomic and ecological resources. These provide an arena for quantitative genetics and allow us to investigate the architecture of complex traits and mechanisms of adaptation to rapidly changing environments. The vast genetic diversity of hybrid populations can yield insights beyond those possible with isogenic lines. Hybrids offer a limitless supply of genetic variation that can be tapped for industrial strain improvement but also, combined with experimental evolution, can be used to predict population responses to future climate change - a fundamental task for biologists.

Patterns of Genomic Instability in Interspecific Yeast Hybrids With Diverse Ancestries.

The genomes of hybrids often show substantial deviations from the features of the parent genomes, including genomic instabilities characterized by chromosomal rearrangements, gains, and losses. This plastic genomic architecture generates phenotypic diversity, potentially giving hybrids access to new ecological niches. It is however unclear if there are any generalizable patterns and predictability in the type and prevalence of genomic variation and instability across hybrids with different genetic and ecological backgrounds. Here, we analyzed the genomic architecture of 204 interspecific Saccharomyces yeast hybrids isolated from natural, industrial fermentation, clinical, and laboratory environments. Synchronous mapping to all eight putative parental species showed significant variation in read depth indicating frequent aneuploidy, affecting 44% of all hybrid genomes and particularly smaller chromosomes. Early generation hybrids with largely equal genomic content from both parent species were more likely to contain aneuploidies than introgressed genomes with an older hybridization history, which presumably stabilized the genome. Shared k-mer analysis showed that the degree of genomic diversity and variability varied among hybrids with different parent species. Interestingly, more genetically distant crosses produced more similar hybrid genomes, which may be a result of stronger negative epistasis at larger genomic divergence, putting constraints on hybridization outcomes. Mitochondrial genomes were typically inherited from the species also contributing the majority nuclear genome, but there were clear exceptions to this rule. Together, we find reliable genomic predictors of instability in hybrids, but also report interesting cross- and environment-specific idiosyncrasies. Our results are an important step in understanding the factors shaping divergent hybrid genomes and their role in adaptive evolution.

Genomic Evidence of an Ancient East Asian Divergence Event in Wild Saccharomyces cerevisiae.

Comparative genome analyses have suggested East Asia to be the cradle of the domesticated microbe Brewer's yeast (Saccharomyces cerevisiae), used in the food and biotechnology industry worldwide. Here, we provide seven new, high-quality long-read genomes of nondomesticated yeast strains isolated from primeval forests and other natural environments in China and Taiwan. In a comprehensive analysis of our new genome assemblies, along with other long-read Saccharomycetes genomes available, we show that the newly sequenced East Asian strains are among the closest living relatives of the ancestors of the global diversity of Brewer's yeast, confirming predictions made from short-read genomic data. Three of these strains (termed the East Asian Clade IX Complex here) share a recent ancestry and evolutionary history suggesting an early divergence from other S. cerevisiae strains before the larger radiation of the species, and prior to its domestication. Our genomic analyses reveal that the wild East Asian strains contain elevated levels of structural variations. The new genomic resources provided here contribute to our understanding of the natural diversity of S. cerevisiae, expand the intraspecific genetic variation found in this heavily domesticated microbe, and provide a foundation for understanding its origin and global colonization history.

Primers to highly conserved elements optimized for qPCR-based telomere length measurement in vertebrates.

Telomere length dynamics are an established biomarker of health and ageing in animals. The study of telomeres in numerous species has been facilitated by methods to measure telomere length by real-time quantitative PCR (qPCR). In this method, telomere length is determined by quantifying the amount of telomeric DNA repeats in a sample and normalizing this to the total amount of genomic DNA. This normalization requires the development of genomic reference primers suitable for qPCR, which remains challenging in nonmodel organism with genomes that have not been sequenced. Here we report reference primers that can be used in qPCR to measure telomere lengths in any vertebrate species. We designed primer pairs to amplify genetic elements that are highly conserved between evolutionarily distant taxa and tested them in species that span the vertebrate tree of life. We report five primer pairs that meet the specificity and reproducibility standards of qPCR. In addition, we demonstrate an approach to choose the best primers for a given species by testing the primers on multiple individuals within a species and then applying an established computational tool. These reference primers can facilitate qPCR-based telomere length measurements in any vertebrate species of ecological or economic interest.