Studying Natural Selection in the Era of Ubiquitous Genomes

Repeated global adaptation across plant species

Global adaptation occurs when all populations of a species undergo selection toward a common optimum. This can occur by a hard selective sweep with the emergence of a new globally advantageous allele that spreads throughout a species' natural range until reaching fixation. This evolutionary process leaves a temporary trace in the region affected, which is detectable using population genomic methods. While selective sweeps have been identified in many species, there have been few comparative and systematic studies of the genes involved in global adaptation. Building upon recent findings showing repeated genetic basis of local adaptation across independent populations and species, we asked whether certain genes play a more significant role in driving global adaptation across plant species. To address this question, we scanned the genomes of 17 plant species to identify signals of repeated global selective sweeps. Despite the substantial evolutionary distance between the species analyzed, we identified several gene families with strong evidence of repeated positive selection. These gene families tend to be enriched for reduced pleiotropy, consistent with predictions from Fisher's evolutionary model and the cost of complexity hypothesis. We also found that genes with repeated sweeps exhibit elevated levels of gene duplication. Our findings contrast with recent observations of increased pleiotropy in genes driving local adaptation, consistent with predictions based on the theory of migration-selection balance.

Minus the Error: Estimating d N/d S and Testing for Natural Selection in the Presence of Residual Alignment Errors

Errors in multiple sequence alignments (MSAs) are known to bias many comparative evolutionary methods. In the context of natural selection analyses, specifically codon evolutionary models, excessive rates of false positives result. A characteristic signature of error-driven findings is unrealistically high estimates of dN/dS (e.g., >100), affecting only a small fraction (e.g., ~0.1%) of the alignment. Despite the widespread use of codon models to assess alignment quality, their potential for error correction remains unexplored. We present BUSTED-E: a novel method designed to detect positive selection while concurrently identifying alignment errors. This method is a straightforward adaptation of the BUSTED flexible branch-site random effects model used to fit distributions of dN/dS, with an important modification: it integrates an "error-sink" component representing an abiological evolutionary regime (dN/dS > 100), and provides the option for masking errors in the MSA for downstream analyses. Statistical performance of BUSTED-E on data simulated without errors shows that there is a small loss of power, which can be mitigated by model averaged inference. Using four published empirical datasets, we show BUSTED-E reduces unrealistic rates of positive selection detection, often by an order of magnitude, and improves biological interpretability of results. BUSTED-E also detects errors that are largely distinct from other popular alignment cleaning tools (HMMCleaner and BMGE). Overall, BUSTED-E is a robust and scalable solution for improving the accuracy of evolutionary analyses in the presence of residual alignment errors, contributing to a more nuanced understanding of natural selection and adaptive evolution.

Somatic genome architecture and molecular evolution are decoupled in "young" linage-specific gene families in ciliates

The evolution of lineage-specific gene families remains poorly studied across the eukaryotic tree of life, with most analyses focusing on the recent evolution of de novo genes in model species. Here we explore the origins of lineage-specific genes in ciliates, a ~1 billion year old clade of microeukaryotes that are defined by their division of somatic and germline functions into distinct nuclei. Previous analyses on conserved gene families have shown the effect of ciliates' unusual genome architecture on gene family evolution: extensive genome processing-the generation of thousands of gene-sized somatic chromosomes from canonical germline chromosomes-is associated with larger and more diverse gene families. To further study the relationship between ciliate genome architecture and gene family evolution, we analyzed lineage specific gene families from a set of 46 transcriptomes and 12 genomes representing x species from eight ciliate classes. We assess how the evolution lineage-specific gene families occurs among four groups of ciliates: extensive fragmenters with gene-size somatic chromosomes, non-extensive fragmenters with "large'' multi-gene somatic chromosomes, Heterotrichea with highly polyploid somatic genomes and Karyorelictea with 'paradiploid' somatic genomes. Our analyses demonstrate that: 1) most lineage-specific gene families are found at shallow taxonomic scales; 2) extensive genome processing (i.e., gene unscrambling) during development likely influences the size and number of young lineage-specific gene families; and 3) the influence of somatic genome architecture on molecular evolution is increasingly apparent in older gene families. Altogether, these data highlight the influences of genome architecture on the evolution of lineage-specific gene families in eukaryotes.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Examining Signatures of Natural Selection in Antifungal Resistance Genes Across Aspergillus Fungi

Certain Aspergillus fungi cause aspergillosis, a set of diseases that typically affect immunocompromised individuals. Most cases of aspergillosis are caused by Aspergillus fumigatus, which infects millions of people annually. Some closely related so-called cryptic species, such as Aspergillus lentulus, can also cause aspergillosis, albeit at lower frequencies, and they are also clinically relevant. Few antifungal drugs are currently available for treating aspergillosis and there is increasing worldwide concern about the presence of antifungal drug resistance in Aspergillus species. Furthermore, isolates from both A. fumigatus and other Aspergillus pathogens exhibit substantial heterogeneity in their antifungal drug resistance profiles. To gain insights into the evolution of antifungal drug resistance genes in Aspergillus, we investigated signatures of positive selection in 41 genes known to be involved in drug resistance across 42 susceptible and resistant isolates from 12 Aspergillus section Fumigati species. Using codon-based site models of sequence evolution, we identified ten genes that contain 43 sites with signatures of ancient positive selection across our set of species. None of the sites that have experienced positive selection overlap with sites previously reported to be involved in drug resistance. These results identify sites that likely experienced ancient positive selection in Aspergillus genes involved in resistance to antifungal drugs and suggest that historical selective pressures on these genes likely differ from any current selective pressures imposed by antifungal drugs.

Global climate change, diet, and the complex relationship between human host and microbiome: Towards an integrated picture

Dietary changes can alter the human microbiome with potential detrimental consequences for health. Given that environment, health, and evolution are interconnected, we ask: Could diet-driven microbiome perturbations have consequences that extend beyond their immediate impact on human health? We address this question in the context of the urgent health challenges posed by global climate change. Drawing on recent studies, we propose that not only can diet-driven microbiome changes lead to dysbiosis, they can also shape life-history traits and fuel human evolution. We posit that dietary shifts prompt mismatched microbiome-host genetics configurations that modulate human longevity and reproductive success. These mismatches can also induce a heritable intra-holobiont stress response, which encourages the holobiont to re-establish equilibrium within the changed nutritional environment. Thus, while mismatches between climate change-related genetic and epigenetic configurations within the holobiont increase the risk and severity of diseases, they may also affect life-history traits and facilitate adaptive responses. These propositions form a framework that can help systematize and address climate-related dietary challenges for policy and health interventions.