Human Immunology through the Lens of Evolutionary Genetics.

Pathogen-imposed selection pressures have been paramount during human evolution. Detecting such selection signatures in ancient and modern human genomes can thus help us to identify genes of temporal and spatial immunological relevance. Admixture with ancient hominins and between human populations has been a source of genetic diversity open to selection by infections. Furthermore, cultural transitions, such as the advent of agriculture, have exposed humans to new microbial threats, with impacts on host defense mechanisms. The integration of population genetics and systems immunology holds great promise for the increased understanding of the factors driving immune response variation between individuals and populations.

Evidence that RNA Viruses Drove Adaptive Introgression between Neanderthals and Modern Humans.

Neanderthals and modern humans interbred at least twice in the past 100,000 years. While there is evidence that most introgressed DNA segments from Neanderthals to modern humans were removed by purifying selection, less is known about the adaptive nature of introgressed sequences that were retained. We hypothesized that interbreeding between Neanderthals and modern humans led to (1) the exposure of each species to novel viruses and (2) the exchange of adaptive alleles that provided resistance against these viruses. Here, we find that long, frequent-and more likely adaptive-segments of Neanderthal ancestry in modern humans are enriched for proteins that interact with viruses (VIPs). We found that VIPs that interact specifically with RNA viruses were more likely to belong to introgressed segments in modern Europeans. Our results show that retained segments of Neanderthal ancestry can be used to detect ancient epidemics.

Ancient familial Mediterranean fever mutations in human pyrin and resistance to Yersinia pestis.

Familial Mediterranean fever (FMF) is an autoinflammatory disease caused by homozygous or compound heterozygous gain-of-function mutations in MEFV, which encodes pyrin, an inflammasome protein. Heterozygous carrier frequencies for multiple MEFV mutations are high in several Mediterranean populations, suggesting that they confer selective advantage. Among 2,313 Turkish people, we found extended haplotype homozygosity flanking FMF-associated mutations, indicating evolutionarily recent positive selection of FMF-associated mutations. Two pathogenic pyrin variants independently arose >1,800 years ago. Mutant pyrin interacts less avidly with Yersinia pestis virulence factor YopM than with wild-type human pyrin, thereby attenuating YopM-induced interleukin (IL)-1ß suppression. Relative to healthy controls, leukocytes from patients with FMF harboring homozygous or compound heterozygous mutations and from asymptomatic heterozygous carriers released heightened IL-1ß specifically in response to Y. pestis. Y. pestis-infected MefvM680I/M680I FMF knock-in mice exhibited IL-1-dependent increased survival relative to wild-type knock-in mice. Thus, FMF mutations that were positively selected in Mediterranean populations confer heightened resistance to Y. pestis.

150 years of Darwin's theory of intercellular flow of hereditary information.

Charles Darwin's Pangenesis theory, which proposed an intercellular mechanism for the flow of hereditary information, is gaining new ground.

Continuously fluctuating selection reveals fine granularity of adaptation.

Temporally fluctuating environmental conditions are a ubiquitous feature of natural habitats. Yet, how finely natural populations adaptively track fluctuating selection pressures via shifts in standing genetic variation is unknown1,2. Here we generated genome-wide allele frequency data every 1-2 generations from a genetically diverse population of Drosophila melanogaster in extensively replicated field mesocosms from late June to mid-December (a period of approximately 12 total generations). Adaptation throughout the fundamental ecological phases of population expansion, peak density and collapse was underpinned by extremely rapid, parallel changes in genomic variation across replicates. Yet, the dominant direction of selection fluctuated repeatedly, even within each of these ecological phases. Comparing patterns of change in allele frequency to an independent dataset procured from the same experimental system demonstrated that the targets of selection are predictable across years. In concert, our results reveal a fitness relevance of standing variation that is likely to be masked by inference approaches based on static population sampling or insufficiently resolved time-series data. We propose that such fine-scaled, temporally fluctuating selection may be an important force contributing to the maintenance of functional genetic variation in natural populations and an important stochastic force impacting genome-wide patterns of diversity at linked neutral sites, akin to genetic draft.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

What Has a Century of Quantitative Genetics Taught Us About Nature's Genetic Tool Kit?

The complexity of heredity has been appreciated for decades: Many traits are controlled not by a single genetic locus but instead by polymorphisms throughout the genome. The importance of complex traits in biology and medicine has motivated diverse approaches to understanding their detailed genetic bases. Here, we focus on recent systematic studies, many in budding yeast, which have revealed that large numbers of all kinds of molecular variation, from noncoding to synonymous variants, can make significant contributions to phenotype. Variants can affect different traits in opposing directions, and their contributions can be modified by both the environment and the epigenetic state of the cell. The integration of prospective (synthesizing and analyzing variants) and retrospective (examining standing variation) approaches promises to reveal how natural selection shapes quantitative traits. Only by comprehensively understanding nature's genetic tool kit can we predict how phenotypes arise from the complex ensembles of genetic variants in living organisms.

Mutation bias reflects natural selection in Arabidopsis thaliana.

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences1. Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana. In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome-mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions (r = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias2 reduces the occurrence of deleterious mutations in Arabidopsis, challenging the prevailing paradigm that mutation is a directionless force in evolution.

Emergence of methicillin resistance predates the clinical use of antibiotics.

The discovery of antibiotics more than 80 years ago has led to considerable improvements in human and animal health. Although antibiotic resistance in environmental bacteria is ancient, resistance in human pathogens is thought to be a modern phenomenon that is driven by the clinical use of antibiotics1. Here we show that particular lineages of methicillin-resistant Staphylococcus aureus-a notorious human pathogen-appeared in European hedgehogs in the pre-antibiotic era. Subsequently, these lineages spread within the local hedgehog populations and between hedgehogs and secondary hosts, including livestock and humans. We also demonstrate that the hedgehog dermatophyte Trichophyton erinacei produces two ß-lactam antibiotics that provide a natural selective environment in which methicillin-resistant S. aureus isolates have an advantage over susceptible isolates. Together, these results suggest that methicillin resistance emerged in the pre-antibiotic era as a co-evolutionary adaptation of S. aureus to the colonization of dermatophyte-infected hedgehogs. The evolution of clinically relevant antibiotic-resistance genes in wild animals and the connectivity of natural, agricultural and human ecosystems demonstrate that the use of a One Health approach is critical for our understanding and management of antibiotic resistance, which is one of the biggest threats to global health, food security and development.

Pandemic-scale phylogenomics reveals the SARS-CoV-2 recombination landscape.

Accurate and timely detection of recombinant lineages is crucial for interpreting genetic variation, reconstructing epidemic spread, identifying selection and variants of interest, and accurately performing phylogenetic analyses1-4. During the SARS-CoV-2 pandemic, genomic data generation has exceeded the capacities of existing analysis platforms, thereby crippling real-time analysis of viral evolution5. Here, we use a new phylogenomic method to search a nearly comprehensive SARS-CoV-2 phylogeny for recombinant lineages. In a 1.6 million sample tree from May 2021, we identify 589 recombination events, which indicate that around 2.7% of sequenced SARS-CoV-2 genomes have detectable recombinant ancestry. Recombination breakpoints are inferred to occur disproportionately in the 3' portion of the genome that contains the spike protein. Our results highlight the need for timely analyses of recombination for pinpointing the emergence of recombinant lineages with the potential to increase transmissibility or virulence of the virus. We anticipate that this approach will empower comprehensive real-time tracking of viral recombination during the SARS-CoV-2 pandemic and beyond.

Genetic mechanisms of animal camouflage: an interdisciplinary perspective.

Camouflage is a classic example of a trait wherein animals respond to natural selection to avoid predation or attract prey. This unique phenomenon has attracted significant recent attention and the rapid development of integrative research methods is facilitating advances in our understanding of the in-depth genetic mechanisms of camouflage. In this review article, we revisit camouflage definitions and strategies and then we examine the underlying mechanisms of the two most common forms of camouflage, crypsis and masquerade, that have recently been elucidated using multiple approaches. We also discuss unresolved questions related to camouflage. Ultimately, we highlight the implications of camouflage for informing various key issues in ecology and evolution.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

Evolving a favorable distribution for mutation effects.

Tandem-repeat DNA sequences appear to be singularly capable of yielding abundant repeat-number mutations with a potentially advantageous distribution of fitness effects. Although knowing the rates and relative proportions of deleterious, neutral and beneficial mutations is fundamental for understanding evolvability, analysis of adaptation routinely overlooks small-effect mutations arising in tandem repeats.

Molecular and evolutionary processes generating variation in gene expression.

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Constraining models of dominance for nonsynonymous mutations in the human genome.

Dominance is a fundamental parameter in genetics, determining the dynamics of natural selection on deleterious and beneficial mutations, the patterns of genetic variation in natural populations, and the severity of inbreeding depression in a population. Despite this importance, dominance parameters remain poorly known, particularly in humans or other non-model organisms. A key reason for this lack of information about dominance is that it is extremely challenging to disentangle the selection coefficient (s) of a mutation from its dominance coefficient (h). Here, we explore dominance and selection parameters in humans by fitting models to the site frequency spectrum (SFS) for nonsynonymous mutations. When assuming a single dominance coefficient for all nonsynonymous mutations, we find that numerous h values can fit the data, so long as h is greater than ~0.15. Moreover, we also observe that theoretically-predicted models with a negative relationship between h and s can also fit the data well, including models with h = 0.05 for strongly deleterious mutations. Finally, we use our estimated dominance and selection parameters to inform simulations revisiting the question of whether the out-of-Africa bottleneck has led to differences in genetic load between African and non-African human populations. These simulations suggest that the relative burden of genetic load in non-African populations depends on the dominance model assumed, with slight increases for more weakly recessive models and slight decreases shown for more strongly recessive models. Moreover, these results also demonstrate that models of partially recessive nonsynonymous mutations can explain the observed severity of inbreeding depression in humans, bridging the gap between molecular population genetics and direct measures of fitness in humans. Our work represents a comprehensive assessment of dominance and deleterious variation in humans, with implications for parameterizing models of deleterious variation in humans and other mammalian species.

The distribution of fitness effects during adaptive walks using a simple genetic network.

The tempo and mode of adaptation depends on the availability of beneficial alleles. Genetic interactions arising from gene networks can restrict this availability. However, the extent to which networks affect adaptation remains largely unknown. Current models of evolution consider additive genotype-phenotype relationships while often ignoring the contribution of gene interactions to phenotypic variance. In this study, we model a quantitative trait as the product of a simple gene regulatory network, the negative autoregulation motif. Using forward-time genetic simulations, we measure adaptive walks towards a phenotypic optimum in both additive and network models. A key expectation from adaptive walk theory is that the distribution of fitness effects of new beneficial mutations is exponential. We found that both models instead harbored distributions with fewer large-effect beneficial alleles than expected. The network model also had a complex and bimodal distribution of fitness effects among all mutations, with a considerable density at deleterious selection coefficients. This behavior is reminiscent of the cost of complexity, where correlations among traits constrain adaptation. Our results suggest that the interactions emerging from genetic networks can generate complex and multimodal distributions of fitness effects.

Polygenic selection to a changing optimum under self-fertilisation.

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

The strength and pattern of natural selection on gene expression in rice.

Levels of gene expression underpin organismal phenotypes1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown1,2. Here we assayed gene expression in rice (Oryza sativa)3, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic6 (0.04%), and transcripts that display lower levels of expression and stochastic noise7-9 and higher levels of plasticity9 are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity9. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes4,5, but that the efficacy of selection is genetically constrained under drought conditions10. Drought selected for earlier flowering11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering13-marking this gene as a drought-escape gene11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Evolutionary and functional analyses of LRP5 in archaic and extant modern humans.

BACKGROUND: The human lineage has undergone a postcranial skeleton gracilization (i.e. lower bone mass and strength relative to body size) compared to other primates and archaic populations such as the Neanderthals. This gracilization has been traditionally explained by differences in the mechanical load that our ancestors exercised. However, there is growing evidence that gracilization could also be genetically influenced. RESULTS: We have analyzed the LRP5 gene, which is known to be associated with high bone mineral density conditions, from an evolutionary and functional point of view. Taking advantage of the published genomes of archaic Homo populations, our results suggest that this gene has a complex evolutionary history both between archaic and living humans and within living human populations. In particular, we identified the presence of different selective pressures in archaics and extant modern humans, as well as evidence of positive selection in the African and South East Asian populations from the 1000 Genomes Project. Furthermore, we observed a very limited evidence of archaic introgression in this gene (only at three haplotypes of East Asian ancestry out of the 1000 Genomes), compatible with a general erasing of the fingerprint of archaic introgression due to functional differences in archaics compared to extant modern humans. In agreement with this hypothesis, we observed private mutations in the archaic genomes that we experimentally validated as putatively increasing bone mineral density. In particular, four of five archaic missense mutations affecting the first ß-propeller of LRP5 displayed enhanced Wnt pathway activation, of which two also displayed reduced negative regulation. CONCLUSIONS: In summary, these data suggest a genetic component contributing to the understanding of skeletal differences between extant modern humans and archaic Homo populations.