Pathogen-mediated selection favours the maintenance of innate immunity gene polymorphism in a widespread wild ungulate

The characterization of toll-like receptor repertoire in Pinna nobilis after mass mortality events suggests adaptive introgression

The fan mussel Pinna nobilis is currently on the brink of extinction due to a multifactorial disease mainly caused to the highly pathogenic parasite Haplosporidium pinnae, meaning that the selection pressure outweighs the adaptive potential of the species. Hopefully, rare individuals have been observed somehow resistant to the parasite, stretching the need to identify the traits underlying this better fitness. Among the candidate to explore at first intention are fast-evolving immune genes, of which toll-like receptor (TLR). In this study, we examined the genetic diversity at 14 TLR loci across P. nobilis, Pinna rudis and P. nobilis × P. rudis hybrid genomes, collected at four physically distant regions, that were found to be either resistant or sensitive to the parasite H. pinnae. We report a high genetic diversity, mainly observed at cell surface TLRs compared with that of endosomal TLRs. However, the endosomal TLR-7 exhibited unexpected level of diversity and haplotype phylogeny. The lack of population structure, associated with a high genetic diversity and elevated dN/dS ratio, was interpreted as balancing selection, though both directional and purifying selection were detected. Interestingly, roughly 40% of the P. nobilis identified as resistant to H. pinnae were introgressed with P. rudis TLR. Specifically, they all carried a TLR-7 of P. rudis origin, whereas sensitive P. nobilis were not introgressed, at least at TLR loci. Small contributions of TLR-6 and TLR-4 single-nucleotide polymorphisms to the clustering of resistant and susceptible individuals could be detected, but their specific role in resistance remains highly speculative. This study provides new information on the diversity of TLR genes within the P. nobilis species after MME and additional insights into adaptation to H. pinnae that should contribute to the conservation of this Mediterranean endemic species.

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.

The Impact of Variation in the Toll-like Receptor 3 Gene on Epizootic Hemorrhagic Disease in Illinois Wild White-Tailed Deer (Odocoileus virginianus)

Epizootic hemorrhagic disease (EHD) leads to high mortality in white-tailed deer (Odocoileus virginianus) and is caused by a double-stranded RNA (dsRNA) virus. Toll-like receptor 3 (TLR3) plays a role in host immune detection and response to dsRNA viruses. We, therefore, examined the role of genetic variation within the TLR3 gene in EHD among 84 Illinois wild white-tailed deer (26 EHD-positive deer and 58 EHD-negative controls). The entire coding region of the TLR3 gene was sequenced: 2715 base pairs encoding 904 amino acids. We identified 85 haplotypes with 77 single nucleotide polymorphisms (SNPs), of which 45 were synonymous mutations and 32 were non-synonymous. Two non-synonymous SNPs differed significantly in frequency between EHD-positive and EHD-negative deer. In the EHD-positive deer, phenylalanine was relatively less likely to be encoded at codon positions 59 and 116, whereas leucine and serine (respectively) were detected less frequently in EHD-negative deer. Both amino acid substitutions were predicted to impact protein structure or function. Understanding associations between TLR3 polymorphisms and EHD provides insights into the role of host genetics in outbreaks of EHD in deer, which may allow wildlife agencies to better understand the severity of outbreaks.

Effect of population size and selection on Toll-like receptor diversity in populations of Galápagos mockingbirds

The interactions of evolutionary forces are difficult to analyse in free-living populations. However, when properly understood, they provide valuable insights into evolutionary biology and conservation genetics. This is particularly important for the interplay of genetic drift and natural selection in immune genes that confer resistance to disease. The Galápagos Islands are inhabited by four closely related species of mockingbirds (Mimus spp.). We used 12 different-sized populations of Galápagos mockingbirds and one population of their continental relative northern mockingbird (Mimus polyglottos) to study the effects of genetic drift on the molecular evolution of immune genes, the Toll-like receptors (TLRs: TLR1B, TLR4 and TLR15). We found that neutral genetic diversity was positively correlated with island size, indicating an important effect of genetic drift. However, for TLR1B and TLR4, there was little correlation between functional (e.g., protein) diversity and island size, and protein structural properties were largely conserved, indicating only a limited effect of genetic drift on molecular phenotype. By contrast, TLR15 was less conserved and even its putative functional polymorphism correlated with island size. The patterns observed for the three genes suggest that genetic drift does not necessarily dominate selection even in relatively small populations, but that the final outcome depends on the degree of selection constraint that is specific for each TLR locus.

Leveraging an existing whole-genome resequencing population data set to characterize toll-like receptor gene diversity in a threatened bird

Species recovery programs are increasingly using genomic data to measure neutral genetic diversity and calculate metrics like relatedness. While these measures can inform conservation management, determining the mechanisms underlying inbreeding depression requires information about functional genes associated with adaptive or maladaptive traits. Toll-like receptors (TLRs) are one family of functional genes, which play a crucial role in recognition of pathogens and activation of the immune system. Previously, these genes have been analysed using species-specific primers and PCR. Here, we leverage an existing short-read reference genome, whole-genome resequencing population data set, and bioinformatic tools to characterize TLR gene diversity in captive and wild tchūriwat'/tūturuatu/shore plover (Thinornis novaeseelandiae), a threatened bird endemic to Aotearoa New Zealand. Our results show that TLR gene diversity in tchūriwat'/tūturuatu is low, and forms two distinct captive and wild genetic clusters. The bioinformatic approach presented here has broad applicability to other threatened species with existing genomic resources in Aotearoa New Zealand and beyond.

Selection balancing at innate immune genes: adaptive polymorphism maintenance in Toll-like receptors

Balancing selection is a classic mechanism for maintaining variability in immune genes involved in host–pathogen interactions. However, it remains unclear how widespread the mechanism is across immune genes other than the major histocompatibility complex (MHC). Although occasional reports suggest that balancing selection (heterozygote advantage, negative frequency-dependent selection, and fluctuating selection) may act on other immune genes, the current understanding of the phenomenon in non-MHC immune genes is far from solid. In this review, we focus on Toll-like receptors (TLRs), innate immune genes directly involved in pathogen recognition and immune response activation, as there is a growing body of research testing the assumptions of balancing selection in these genes. After reviewing infection- and fitness-based evidence, along with evidence based on population allelic frequencies and heterozygosity levels, we conclude that balancing selection maintains variation in TLRs, though it tends to occur under specific conditions in certain evolutionary lineages rather than being universal and ubiquitous. Our review also identifies key gaps in current knowledge and proposes promising areas for future research. Improving our understanding of host–pathogen interactions and balancing selection in innate immune genes are increasingly important, particularly regarding threats from emerging zoonotic diseases.

Selection and demography drive range-wide patterns of MHC-DRB variation in mule deer

BACKGROUND

Standing genetic variation is important especially in immune response-related genes because of threats to wild populations like the emergence of novel pathogens. Genetic variation at the major histocompatibility complex (MHC), which is crucial in activating the adaptive immune response, is influenced by both natural selection and historical population demography, and their relative roles can be difficult to disentangle. To provide insight into the influences of natural selection and demography on MHC evolution in large populations, we analyzed geographic patterns of variation at the MHC class II DRB exon 2 locus in mule deer (Odocoileus hemionus) using sequence data collected across their entire broad range.

RESULTS

We identified 31 new MHC-DRB alleles which were phylogenetically similar to other cervid MHC alleles, and one allele that was shared with white-tailed deer (Odocoileus virginianus). We found evidence for selection on the MHC including high dN/dS ratios, positive neutrality tests, deviations from Hardy-Weinberg Equilibrium (HWE) and a stronger pattern of isolation-by-distance (IBD) than expected under neutrality. Historical demography also shaped variation at the MHC, as indicated by similar spatial patterns of variation between MHC and microsatellite loci and a lack of association between genetic variation at either locus type and environmental variables.

CONCLUSIONS

Our results show that both natural selection and historical demography are important drivers in the evolution of the MHC in mule deer and work together to shape functional variation and the evolution of the adaptive immune response in large, well-connected populations.

Neutral and Selective Processes Shape MHC Diversity in Roe Deer in Slovenia

Simple Summary

Disease prevention and appropriate wildlife management are among the major challenges in wildlife conservation. In the present study, we made a first assessment of the variability of major histocompatibility complex (MHC) genes in roe deer in Slovenia and evaluated local population adaptation by comparing MHC variability with neutral microsatellites. We discovered three new MHC DRB exon 2 alleles in addition to seven previously described in the literature. Moreover, we found evidence of historical positive selection, as selection analysis indicated that approx. 10% of the encoded amino acids were subjected to episodic positive selection. This study provides the basis for further research on immunogenetic variation in roe deer and highlights opportunities to incorporate genetic data into science-based population management.

Abstract

Disease control and containment in free-ranging populations is one of the greatest challenges in wildlife management. Despite the importance of major histocompatibility complex (MHC) genes for immune response, an assessment of the diversity and occurrence of these genes is still rare in European roe deer, the most abundant and widespread large mammal in Europe. Therefore, we examined immunogenetic variation in roe deer in Slovenia to identify species adaptation by comparing the genetic diversity of the MHC genes with the data on neutral microsatellites. We found ten MHC DRB alleles, three of which are novel. Evidence for historical positive selection on the MHC was found using the maximum likelihood codon method. Patterns of MHC allelic distribution were not congruent with neutral population genetic findings. The lack of population genetic differentiation in MHC genes compared to existing structure in neutral markers suggests that MHC polymorphism was influenced primarily by balancing selection and, to a lesser extent, by neutral processes such as genetic drift, with no clear evidence of local adaptation. Selection analyses indicated that approx. 10% of amino acids encoded under episodic positive selection. This study represents one of the first steps towards establishing an immunogenetic map of roe deer populations across Europe, aiming to better support science-based management of this important game species.