Large-effect Mutations Research Articles

Energy landscapes of peptide-MHC binding.

Molecules of the Major Histocompatibility Complex (MHC) present short protein fragments on the cell surface, an important step in T cell immune recognition. MHC-I molecules process peptides from intracellular proteins; MHC-II molecules act in antigen-presenting cells and present peptides derived from extracellular proteins. Here we show that the sequence-dependent energy landscapes of MHC-peptide binding encode class-specific nonlinearities (epistasis). MHC-I has a smooth landscape with global epistasis; the binding energy is a simple deformation of an underlying linear trait. This form of epistasis enhances the discrimination between strong-binding peptides. In contrast, MHC-II has a rugged landscape with idiosyncratic epistasis: binding depends on detailed amino acid combinations at multiple positions of the peptide sequence. The form of epistasis affects the learning of energy landscapes from training data. For MHC-I, a low-complexity problem, we derive a simple matrix model of binding energies that outperforms current models trained by machine learning. For MHC-II, higher complexity prevents learning by simple regression methods. Epistasis also affects the energy and fitness effects of mutations in antigen-derived peptides (epitopes). In MHC-I, large-effect mutations occur predominantly in anchor positions of strong-binding epitopes. In MHC-II, large effects depend on the background epitope sequence but are broadly distributed over the epitope, generating a bigger target for escape mutations due to loss of presentation. Together, our analysis shows how an energy landscape of protein-protein binding constrains the target of escape mutations from T cell immunity, linking the complexity of the molecular interactions to the dynamics of adaptive immune response.

Environmental memory alters the fitness effects of adaptive mutations in fluctuating environments.

Evolution in a static laboratory environment often proceeds via large-effect beneficial mutations that may become maladaptive in other environments. Conversely, natural settings require populations to endure environmental fluctuations. A sensible assumption is that the fitness of a lineage in a fluctuating environment is the time average of its fitness over the sequence of static conditions it encounters. However, transitions between conditions may pose entirely new challenges, which could cause deviations from this time average. To test this, we tracked hundreds of thousands of barcoded yeast lineages evolving in static and fluctuating conditions and subsequently isolated 900 mutants for pooled fitness assays in 15 environments. Here we find that fitness in fluctuating environments indeed often deviates from the time average, leading to fitness non-additivity. Moreover, closer examination reveals that fitness in one component of a fluctuating environment is often strongly influenced by the previous component. We show that this environmental memory is especially common for mutants with high variance in fitness across tested environments. We use a simple mathematical model and whole-genome sequencing to propose mechanisms underlying this effect, including lag time evolution and sensing mutations. Our results show that environmental fluctuations impact fitness and suggest that variance in static environments can explain these impacts.

Impacts of pleiotropy and migration on repeated genetic adaptation.

Observations of genetically repeated evolution (repeatability) in complex organisms are incongruent with the Fisher-Orr model, which implies that repeated use of the same gene should be rare when mutations are pleiotropic (i.e. affect multiple traits). When spatially divergent selection occurs in the presence of migration, mutations of large effect are more strongly favored, and hence, repeatability is more likely, but it is unclear whether this observation is limited by pleiotropy. Here, we explore this question using individual-based simulations of a two-patch model incorporating multiple quantitative traits governed by mutations with pleiotropic effects. We explore the relationship between fitness trade-offs and repeatability by varying the alignment between mutation effect and spatial variation in trait optima. While repeatability decreases with increasing trait dimensionality, trade-offs in mutation effects on traits do not strongly limit the contribution of a locus of large effect to repeated adaptation, particularly under increased migration. These results suggest that repeatability will be more pronounced for local rather than global adaptation. Whereas pleiotropy limits repeatability in a single-population model, when there is local adaptation with gene flow, repeatability can occur if some loci are able to produce alleles of large effect, even when there are pleiotropic trade-offs.

The loci of environmental adaptation in a model eukaryote

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Evolutionary rescue in resistance to pesticides.

Evolutionary rescue occurs when the genetic evolution of adaptation saves a population from decline or extinction after environmental change. The evolution of resistance to pesticides is a special scenario of abrupt environmental change, where rescue occurs under (very) strong selection for one or a few de novo resistance mutations of large effect. Here, a population genetic model of evolutionary rescue with density-dependent population change is developed, with a focus on deriving results that are important to resistance management. Massive stochastic simulations are used to generate observations, which are accurately predicted using analytical approximations. Key results include the probability density function for the time to resistance and the probability of population extinction. The distribution of resistance times shows a lag period, a narrow peak and a long tail. Surprisingly, the mean time to resistance can increase with the strength of selection because, if a mutation does not occur early on, then its emergence is delayed by the pesticide reducing the population size. The probability of population extinction shows a sharp transition, in that when extinction is possible, it is also highly likely. Consequently, population suppression and (local) eradication can be theoretically achievable goals, as novel strategies to delay resistance evolution.

The Effects of De Novo Mutation on Gene Expression and the Consequences for Fitness in Chlamydomonas reinhardtii.

Mutation is the ultimate source of genetic variation, the bedrock of evolution. Yet, predicting the consequences of new mutations remains a challenge in biology. Gene expression provides a potential link between a genotype and its phenotype. But the variation in gene expression created by de novo mutation and the fitness consequences of mutational changes to expression remain relatively unexplored. Here, we investigate the effects of >2,600 de novo mutations on gene expression across the transcriptome of 28 mutation accumulation lines derived from 2 independent wild-type genotypes of the green algae Chlamydomonas reinhardtii. We observed that the amount of genetic variance in gene expression created by mutation (Vm) was similar to the variance that mutation generates in typical polygenic phenotypic traits and approximately 15-fold the variance seen in the limited species where Vm in gene expression has been estimated. Despite the clear effect of mutation on expression, we did not observe a simple additive effect of mutation on expression change, with no linear correlation between the total expression change and mutation count of individual MA lines. We therefore inferred the distribution of expression effects of new mutations to connect the number of mutations to the number of differentially expressed genes (DEGs). Our inferred DEE is highly L-shaped with 95% of mutations causing 0-1 DEG while the remaining 5% are spread over a long tail of large effect mutations that cause multiple genes to change expression. The distribution is consistent with many cis-acting mutation targets that affect the expression of only 1 gene and a large target of trans-acting targets that have the potential to affect tens or hundreds of genes. Further evidence for cis-acting mutations can be seen in the overabundance of mutations in or near differentially expressed genes. Supporting evidence for trans-acting mutations comes from a 15:1 ratio of DEGs to mutations and the clusters of DEGs in the co-expression network, indicative of shared regulatory architecture. Lastly, we show that there is a negative correlation with the extent of expression divergence from the ancestor and fitness, providing direct evidence of the deleterious effects of perturbing gene expression.

Genome-wide association in Drosophila identifies a role for Piezo and Proc-R in sleep latency.

Sleep latency, the amount of time that it takes an individual to fall asleep, is a key indicator of sleep need. Sleep latency varies considerably both among and within species and is heritable, but lacks a comprehensive description of its underlying genetic network. Here we conduct a genome-wide association study of sleep latency. Using previously collected sleep and activity data on a wild-derived population of flies, we calculate sleep latency, confirming significant, heritable genetic variation for this complex trait. We identify 520 polymorphisms in 248 genes contributing to variability in sleep latency. Tests of mutations in 23 candidate genes and additional putative pan-neuronal knockdown of 9 of them implicated CG44153, Piezo, Proc-R and Rbp6 in sleep latency. Two large-effect mutations in the genes Proc-R and Piezo were further confirmed via genetic rescue. This work greatly enhances our understanding of the genetic factors that influence variation in sleep latency.

Schrödinger’s yeast: the challenge of using transformation to compare fitness among Saccharomyces cerevisiae that differ in ploidy or zygosity

How the number of genome copies modifies the effect of random mutations remains poorly known. In yeast, researchers have investigated these effects for knock-out or other large-effect mutations, but have not accounted for differences at the mating-type locus. We set out to compare fitness differences among strains that differ in ploidy and/or zygosity using a panel of spontaneously arising mutations acquired in haploid yeast from a previous study. To ensure no genetic differences, even at the mating-type locus, we embarked on a series of transformations, which first sterilized and then temporarily introduced plasmid-borne mating types. Despite these attempts to equalize the haplotypes, fitness variation introduced during transformation swamped the differences among the original mutation-accumulation lines. While colony size looked normal, we observed a bi-modality in the maximum growth rate of our transformed yeast and determined that many of the slow growing lines were respiratory deficient (“petite”). Not previously reported, we found that yeast that were TID1/RDH54 knockouts were less likely to become petite. Even for lines with the same petite status, however, we found no correlation in fitness between the two replicate transformations performed. These results pose a challenge for any study using transformation to measure the fitness effect of genetic differences among strains. By attempting to hold haplotypes constant, we introduced more mutations that overwhelmed our ability to measure fitness differences between the genetic states. In this study, we transformed over one hundred different lines of yeast, using two independent transformations, and found that this common laboratory procedure can cause large changes to the microbe studied. Our study provides a cautionary tale of the need to use multiple transformants in fitness assays.

Forces driving transposable element load variation during Arabidopsis range expansion.

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Environmental dependence of mutational (co)variances of adaptive traits.

Standing genetic variation, and capacity to adapt to environment change, will ultimately depend on the fitness effects of mutations across the range of environments experienced by contemporary, panmictic, populations. We investigated how mild perturbations in diet and temperature affect mutational (co)variances of traits that evolve under climatic adaptation, and contribute to individual fitness in Drosophila serrata. We assessed egg-to-adult viability, development time and wing size of 64 lines that had diverged from one another via spontaneous mutation over 30 generations of brother-sister mating. Our results suggested most mutations have directionally concordant (i.e., synergistic) effects in all environments and both sexes. However, elevated mutational variance under reduced macronutrient conditions suggested environment-dependent variation in mutational effect sizes for development time. We also observed evidence for antagonistic effects under standard versus reduced macronutrient conditions, where these effects were further contingent on temperature (for development time) or sex (for size). Diet also influenced the magnitude and sign of mutational correlations between traits, although this result was largely due to a single genotype (line), which may reflect a rare, large effect mutation. Overall, our results suggest environmental heterogeneity and environment-dependency of mutational effects could contribute to the maintenance of genetic variance.

Host genomic variation shapes gut microbiome diversity in threespine stickleback fish.

A major focus of host-microbe research is to understand how genetic differences, of various magnitudes, among hosts translate to differences in their microbiomes. This has been challenging for animal hosts, including humans, because it is difficult to control environmental variables tightly enough to isolate direct genetic effects on the microbiome. Our work in stickleback fish is a significant contribution because our experimental approach allowed strict control over environmental factors, including standardization of the microbiome from the earliest stage of development and unrestricted co-housing of fish in a truly common environment. Furthermore, we measured host genetic variation over 2,000 regions of the stickleback genome, comparing this information and microbiome composition data among fish from very similar and very different genetic backgrounds. Our findings highlight how differences in the host genome influence microbiome diversity and make a case for future manipulative microbiome experiments that use host systems with naturally occurring genetic variation.

Effects of population size change on the genetics of adaptation following an abrupt change in environment.

Since the rediscovery of Mendelian genetics over a century ago, there has been much debate about the evolutionary importance of mutations with large phenotypic effects. While population genetic models predict that large-effect mutations will typically contribute to adaptation following an abrupt change in environment, the prediction applies to populations of stable size and overlooks the effects of population size change on adaptation (e.g., population decline following habitat loss; growth during range expansion). We evaluate the phenotypic and fitness effects of mutations contributing to adaptation immediately following an abrupt environmental shift that alters both selection and population size dynamics. We show that large-effect mutations are likely to contribute to adaptation in populations declining to a new carrying capacity, somewhat smaller-effect mutations contribute to evolutionary rescue, and small-effect mutations predominate in growing populations. We also show that the relative contributions of positively selected and overdominant mutations to adaptation depend on interactions between the phenotypic effect size distribution for new mutations and the specific form of population size change during adaptation (i.e., growth, decline, or evolutionary rescue). Our results illustrate how population size dynamics can shape the genetic basis of adaptation, which should motivate empirical comparisons of populations adapting in different demographic contexts.

High-Density Linkage Map Constructed from a Skim Sequenced Diploid Potato Population Reveals Transmission Distortion and QTLs for Tuber Yield and Pollen Shed

The reinvention of potato, from a tetraploid clonal crop into a diploid seed-based hybrid crop, requires insight in the mutational load, recombination landscape, and the genetic basis of fertility. Genomics-based breeding and QTL discovery rely on efficient genotyping strategies such as skim sequencing, to gather genotypic information. The application of skim sequencing to full-sib population of non-inbred parents remains challenging. Here, we report on an R implementation of the OutcrossSeq pipeline for diploids. We applied this pipeline to a large diploid skim sequenced potato population. We used the resulting bin-markers for the construction of high-density parent specific linkage maps, highlighting variation in parental recombination rate and structural variations. We subsequently explored transmission ratio distortion and non-independent assortment of alleles, indicative of large-effect deleterious mutations. Finally, we identified QTLs for seedling tuber yield in pots and pollen shed. This study showcases the range of genetic analyses, from marker inference, identification of transmission ratio distortion, and linkage map construction to QTL mapping, resulting in new insights that contribute to breeding diploid potato.

The landscape of antibody binding affinity in SARS-CoV-2 Omicron BA.1 evolution.

The Omicron BA.1 variant of SARS-CoV-2 escapes convalescent sera and monoclonal antibodies that are effective against earlier strains of the virus. This immune evasion is largely a consequence of mutations in the BA.1 receptor binding domain (RBD), the major antigenic target of SARS-CoV-2. Previous studies have identified several key RBD mutations leading to escape from most antibodies. However, little is known about how these escape mutations interact with each other and with other mutations in the RBD. Here, we systematically map these interactions by measuring the binding affinity of all possible combinations of these 15 RBD mutations (215=32,768 genotypes) to 4 monoclonal antibodies (LY-CoV016, LY-CoV555, REGN10987, and S309) with distinct epitopes. We find that BA.1 can lose affinity to diverse antibodies by acquiring a few large-effect mutations and can reduce affinity to others through several small-effect mutations. However, our results also reveal alternative pathways to antibody escape that does not include every large-effect mutation. Moreover, epistatic interactions are shown to constrain affinity decline in S309 but only modestly shape the affinity landscapes of other antibodies. Together with previous work on the ACE2 affinity landscape, our results suggest that the escape of each antibody is mediated by distinct groups of mutations, whose deleterious effects on ACE2 affinity are compensated by another distinct group of mutations (most notably Q498R and N501Y).

Breeding Diploid F1 Hybrid Potatoes for Propagation from Botanical Seed (TPS): Comparisons with Theory and Other Crops.

This paper reviews the progress and the way ahead in diploid F1 hybrid potato breeding by comparisons with expectations from the theory of inbreeding and crossbreeding, and experiences from other diploid outbreeding crops. Diploid potatoes can be converted from an outbreeding species, in which self-pollination is prevented by a gametophytic self-incompatibility system, into one where self-pollination is possible, either through a dominant self-incompatibility inhibitor gene (Sli) or knockout mutations in the incompatibility locus. As a result, diploid F1 hybrid breeding can be used to produce genetically uniform potato cultivars for propagation from true potato seeds by crossing two near-homozygous inbred lines, derived from a number of generations of self-pollination despite inbreeding depression. Molecular markers can be used to detect and remove deleterious recessive mutations of large effect, including those in tight repulsion linkage. Improvements to the inbred lines can be made by introducing and stacking genes and chromosome segments of large desirable effect from wild relatives by backcrossing. Improvements in quantitative traits require a number of cycles of inbreeding and crossbreeding. Seed production can be achieved by hand pollinations. F1 hybrid planting material can be delivered to farmers as true seeds or young plants, and mini-tubers derived from true seeds.

Novel Protein-Protein Interactions Highlighting the Crosstalk between Hypoplastic Left Heart Syndrome, Ciliopathies and Neurodevelopmental Delays.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease (CHD) affecting 1 in 5000 newborns. We constructed the interactome of 74 HLHS-associated genes identified from a large-scale mouse mutagenesis screen, augmenting it with 408 novel protein–protein interactions (PPIs) using our High-Precision Protein–Protein Interaction Prediction (HiPPIP) model. The interactome is available on a webserver with advanced search capabilities. A total of 364 genes including 73 novel interactors were differentially regulated in tissue/iPSC-derived cardiomyocytes of HLHS patients. Novel PPIs facilitated the identification of TOR signaling and endoplasmic reticulum stress modules. We found that 60.5% of the interactome consisted of housekeeping genes that may harbor large-effect mutations and drive HLHS etiology but show limited transmission. Network proximity of diabetes, Alzheimer’s disease, and liver carcinoma-associated genes to HLHS genes suggested a mechanistic basis for their comorbidity with HLHS. Interactome genes showed tissue-specificity for sites of extracardiac anomalies (placenta, liver and brain). The HLHS interactome shared significant overlaps with the interactomes of ciliopathy- and microcephaly-associated genes, with the shared genes enriched for genes involved in intellectual disability and/or developmental delay, and neuronal death pathways, respectively. This supported the increased burden of ciliopathy variants and prevalence of neurological abnormalities observed among HLHS patients with developmental delay and microcephaly, respectively.

Mutators Enhance Adaptive Micro-Evolution in Pathogenic Microbes.

Adaptation to the changing environmental conditions experienced within a host requires genetic diversity within a microbial population. Genetic diversity arises from mutations which occur due to DNA damage from exposure to exogenous environmental stresses or generated endogenously through respiration or DNA replication errors. As mutations can be deleterious, a delicate balance must be obtained between generating enough mutations for micro-evolution to occur while maintaining fitness and genomic integrity. Pathogenic microorganisms can actively modify their mutation rate to enhance adaptive micro-evolution by increasing expression of error-prone DNA polymerases or by mutating or decreasing expression of genes required for DNA repair. Strains which exhibit an elevated mutation rate are termed mutators. Mutators are found in varying prevalence in clinical populations where large-effect beneficial mutations enhance survival and are predominately caused by defects in the DNA mismatch repair (MMR) pathway. Mutators can facilitate the emergence of antibiotic resistance, allow phenotypic modifications to prevent recognition and destruction by the host immune system and enable switching to metabolic and cellular morphologies better able to survive in the given environment. This review will focus on recent advances in understanding the phenotypic and genotypic changes occurring in MMR mutators in both prokaryotic and eukaryotic pathogens.

Domestication reshaped the genetic basis of inbreeding depression in a maize landrace compared to its wild relative, teosinte.

Inbreeding depression is the reduction in fitness and vigor resulting from mating of close relatives observed in many plant and animal species. The extent to which the genetic load of mutations contributing to inbreeding depression is due to large-effect mutations versus variants with very small individual effects is unknown and may be affected by population history. We compared the effects of outcrossing and self-fertilization on 18 traits in a landrace population of maize, which underwent a population bottleneck during domestication, and a neighboring population of its wild relative teosinte. Inbreeding depression was greater in maize than teosinte for 15 of 18 traits, congruent with the greater segregating genetic load in the maize population that we predicted from sequence data. Parental breeding values were highly consistent between outcross and selfed offspring, indicating that additive effects determine most of the genetic value even in the presence of strong inbreeding depression. We developed a novel linkage scan to identify quantitative trait loci (QTL) representing large-effect rare variants carried by only a single parent, which were more important in teosinte than maize. Teosinte also carried more putative juvenile-acting lethal variants identified by segregation distortion. These results suggest a mixture of mostly polygenic, small-effect partially recessive effects in linkage disequilibrium underlying inbreeding depression, with an additional contribution from rare larger-effect variants that was more important in teosinte but depleted in maize following the domestication bottleneck. Purging associated with the maize domestication bottleneck may have selected against some large effect variants, but polygenic load is harder to purge and overall segregating mutational burden increased in maize compared to teosinte.

Allen Orr and the genetics of adaptation.

Over most of the 20th century, evolutionary biologists predominantly subscribed to a strong form of "micro-mutationism," in which adaptive phenotypic divergence arises from allele frequency changes at many loci, each with a small effect on the phenotype. To be sure, there were well-known examples of large-effect alleles contributing to adaptation, yet such cases were generally regarded as atypical and unrepresentative of evolutionary change in general. In 1998, Allen Orr published a landmark theoretical paper in Evolution, which showed that both small- and large-effect mutations are likely to contribute to "adaptive walks" of a population to an optimum. Coupled with a growing set of empirical examples of large-effect alleles contributing to divergence (e.g., from QTL studies), Orr's paper provided a mathematical formalism that converted many evolutionary biologists from micro-mutationism to a more pluralistic perspective on the genetic basis of evolutionary change. We revisit the theoretical insights emerging from Orr's paper within the historical context leading up to 1998, and track the influence of this paper on the field of evolutionary biology through an examination of its citations over the last two decades and an analysis of the extensive body of theoretical and empirical research that Orr's pioneering paper inspired.

The magic of genomics in creating hybrid potato

The magic of genomics in creating hybrid potato