Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Moderate Genetic Diversity of MHC Genes in an Isolated Small Population of Black-and-White Snub-Nosed Monkeys (Rhinopithecus bieti)

Genetic diversity is an essential indicator that echoes the natural selection and environmental adaptation of a species. Isolated small populations are vulnerable to genetic drift, inbreeding, and limited gene flow; thus, assessing their genetic diversity is critical in conservation. In this study, we studied the genetic diversity of black-and-white snub-nosed monkeys (Rhinopithecus bieti) using neutral microsatellites and five adaptive major histocompatibility complex (MHC) genes. Two DQA1 alleles, two DQB1 alleles, two DRB1 alleles, two DRB5 alleles, and three DPB1 alleles were isolated from a population. The results indicate that neutral microsatellites demonstrate a high degree of heterozygosity and polymorphism, while adaptive MHC genes display a high degree of heterozygosity and moderate polymorphism. The results also show that balancing selection has prominently influenced the MHC diversity of the species during evolution: (1) significant positive selection is identified at several amino acid sites (primarily at and near antigen-binding sites) of the DRB1, DRB5, and DQB1 genes; (2) phylogenetic analyses display the patterns of trans-species evolution for all MHC loci. This study provides valuable genetic diversity insights into black-and-white snub-nosed monkeys, which dwell at the highest altitude and have experienced the harshest environmental selection of all primates globally since the Pleistocene. Such results provide valuable scientific evidence and a reference for making or amending conservation strategies for this endangered primate species.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

The quest for good genes: Epigamic traits, fitness, MHC and multilocus heterozygosity in the guppy

The 'good genes' hypothesis for the evolution of male secondary sexual traits poses that female preferences for such traits are driven by indirect genetic benefits. However, support for the hypothesis remains ambiguous, and, in particular, the genetic basis for the benefits has rarely been investigated. Here, we use seminatural populations of Trinidadian guppies to investigate whether sexually selected traits (orange, black and iridescent colouration, gonopodium length and body size) predict fitness measured as the number of grandoffspring, a metric that integrates across fitness components and sexes. Furthermore, we tested whether two potential sources of genetic benefits-major histocompatibility complex (MHC) genotypes and multilocus heterozygosity (MLH)-are significant predictors of fitness and of the size of sexually selected traits. We found a significant, nonlinear effect of the area of black pigmentation and male body size on the number of grandoffspring, suggesting stabilizing selection on black area, and nonlinear selection favouring small body size. MLH was heritable (h2  = 0.14) and significantly predicted the number of grandoffspring, indicating the potential for genetic benefits based on heterozygosity. We also found support for local heterozygosity effects, which may reflect a noneven distribution of genetic load across the genome. MHC genotype was not significantly associated with any tested fitness component, or with the load of Gyrodactylus parasites. Neither MHC nor MLH was significant predictor of sexually selected traits. Overall, our results highlight the role of heterozygosity in determining fitness, but do not provide support for male sexually selected traits being indicators of genetic quality.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Diversifying selection and climatic effects on major histocompatibility complex class II gene diversity in the greater horseshoe bat

Heterogeneous pathogenic stress can shape major histocompatibility complex (MHC) diversity by influencing the functional plasticity of the immune response. Therefore, MHC diversity could reflect environmental stress, demonstrating its importance in uncovering the mechanisms of adaptive genetic variation. In this study, we combined neutral microsatellite loci, an immune-related MHC II-DRB locus, and climatic factors to unravel the mechanisms affecting the diversity and genetic differentiation of MHC genes in the greater horseshoe bat (Rhinolophus ferrumequinum), a species with a wide geographical distribution that has three distinct genetic lineages in China. First, increased genetic differentiation at the MHC locus among populations compared using microsatellites indicated diversifying selection. Second, the genetic differentiation of MHC and microsatellites were significantly correlated, suggesting that demographic processes exist. However, MHC genetic differentiation was significantly correlated with geographical distance among populations, even after controlling for the neutral markers, suggesting a major effect of selection. Third, although the MHC genetic differentiation was larger than that for microsatellites, there was no significant difference in the genetic differentiation between the two markers among genetic lineages, indicating the effect of balancing selection. Fourth, combined with climatic factors, MHC diversity and supertypes showed significant correlations with temperature and precipitation, but not with the phylogeographic structure of R. ferrumequinum, suggesting an effect of local adaptation driven by climate on MHC diversity. Moreover, the number of MHC supertypes varied between populations and lineages, suggesting regional characteristics and support for local adaptation. Taken together, the results of our study provide insights into the adaptive evolutionary driving forces at different geographic scales in R. ferrumequinum. In addition, climate factors may have played a vital role in driving adaptive evolution in this species.

Recombination, selection and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.

Urbanization processes drive divergence at the major histocompatibility complex in a common waterbird

Urban sprawl is one of the most common landscape alterations occurring worldwide, and there is a growing list of species that are recognised to have adapted to urban life. To be successful, processes of urban colonization by wildlife require a broad spectrum of phenotypic (e.g., behavioural or physiological) adjustments, but evidence for genetic adaptations is much scarcer. One hypothesis proposes that different pathogen-driven selective pressures between urban and non-urban landscapes leads to adaptations in host immune genes. Here, we examined urbanization-related differentiation at the key pathogen-recognition genes of vertebrate adaptive immunity-the major histocompatibility complex (MHC)-in a common waterbird, the Eurasian coot (Fulica atra). Samples were collected from an old urban population (established before the 1950s), a new urban population (established in the 2000s), and two rural populations from central Poland. We found strong significant divergence (as measured with Jost's D) at the MHC class II between the old urban population and the remaining (new urban and rural) populations. Also, there was a moderate, but significant divergence at the MHC between the new urban population and two rural populations, while no divergence was found between the two rural populations. The total number of MHC alleles and the number of private (population-specific) MHC alleles was lower in old urban populations, as compared to the rural ones. These patterns of differentiation at the MHC were not consistent with patterns found for neutral genetic markers (microsatellites), which showed few differences between the populations. Our results indicate that MHC allele composition depended on the level of anthropogenic disturbance and the time which passed since urban colonization, possibly due to the processes of genotype sorting and local adaptation. As such, our study contributes to the understanding of genetic mechanisms associated with urbanization processes in wildlife.

Depletion of MHC supertype during domestication can compromise immunocompetence

The major histocompatibility complex (MHC) plays an important role in infectious disease resistance. The presence of certain MHC alleles and functionally similar groups of MHC alleles (i.e., supertypes) has been associated with resistance to particular parasite species. Farmed and domesticated fish stocks are often depleted in their MHC alleles and supertype diversity, possibly as a consequence of artificial selection for desirable traits, inbreeding (loss of heterozygosity), genetic drift (loss of allelic diversity) and/or reduced parasite biodiversity. Here we quantify the effects of depletion of MHC class II genotype and supertype variation on resistance to the parasite Gyrodactylus turnbulli in guppies (Poecilia reticulata). Compared to the descendants of wild‐caught guppies, ornamental fish had a significantly reduced MHC variation (i.e., the numbers of MHC alleles and supertypes per individual, and per population). In addition, ornamental fish were significantly more susceptible to G. turnbulli infections, accumulating peak intensity 10 times higher than that of their wildtype counterparts. Four out of 13 supertypes were associated with a significantly reduced parasite load, and the presence of some supertypes had a dramatic effect on the intensity of infection. Remarkably, the ornamental and wildtype fish differed in the supertypes that were associated with parasite resistance. Analysis with a genetic algorithm showed that resistance‐conferring supertypes of the wildtype and ornamental fish shared two unique amino acids in the peptide‐binding region of the MHC that were not found in any other alleles. These data show that the supertype demarcation captures some, but not all, of the variation in the immune function of the alleles. This study highlights the importance of managing functional MHC diversity in livestock, and suggests there might be some immunological redundancy among MHC supertypes.

Balancing selection versus allele and supertype turnover in MHC class II genes in guppies

Selection pressure from parasites is thought to be a major force shaping the extreme polymorphism of the major histocompatibility complex (MHC) genes, but the modes and consequences of selection remain unclear. Here, we analyse MHC class II and microsatellite diversity in 16 guppy (Poecilia reticulata) populations from two islands (Trinidad and Tobago) that have been separated for at least 10 ky. Within-population MHC diversity was high, but allele sharing was limited within islands and even lower between islands, suggesting relatively fast turnover of alleles. Allelic lineages strongly supported in phylogenetic analyses tended to be island-specific, suggesting rapid lineage sorting, and an expansion of an allelic lineage private to Tobago was observed. New alleles appear to be generated locally at a detectably high frequency. We did not detect a consistent signature of local adaptation, but F_ST outlier analysis suggested that balancing selection may be the more general process behind spatial variation in MHC allele frequencies in this system, particularly within Trinidad. We found no evidence for divergent allele advantage within populations, or for decreased genetic structuring of MHC supertypes compared to MHC alleles. The dynamic and complex nature of MHC evolution we observed in guppies, coupled with some evidence for balancing selection shaping MHC allele frequencies, are consistent with Red Queen processes of host-parasite coevolution.

Identification of Novel Alleles and Structural Haplotypes of Major Histocompatibility Complex Class I and DRB Genes in Domestic Cat (Felis catus) by a Newly Developed NGS-Based Genotyping Method

The major histocompatibility complex (MHC) is a highly polymorphic and duplicated genomic region that encodes transplantation and immune regulatory molecules. Although it is well-known that particular MHC allelic polymorphisms and haplotypes are genetically relate to immune-mediated diseases detailed information of the cat MHC (Feline Leukocyte Antigen; FLA) genetic and haplotypic structure and diversity is limited in comparison to humans and many other species. In this study, to better understand the degree and types of allele and allelic haplotype diversity of FLA-class I (FLA-I) and FLA-DRB loci in domestic cats, we identified six expressible FLA-I loci in peripheral white blood cells by in silico estimation of the coding exons and NGS-based amplicon sequencing using five unrelated cats. We then used a newly developed NGS-based genotyping method to genotype and annotate 32 FLA-I and 16 FLA-DRB sequences in two families of 20 domestic cats. A total of 14 FLA-I and seven FLA-DRB were identified as novel polymorphic sequences. Phylogenetic analyses grouped the sequences into six FLA-I (FLA-E/H/K, FLA-A, FLA-J, FLA-L, FLA-O and a tentatively named FLA-E/H/K_Rec) and four FLA-DRB (FLA-DRB1, FLA-DRB3, FLA-DRB4, and FLA-DRB5) lineages. Pedigree analysis of two cat families revealed eight distinct FLA structural haplotypes (Class I - DRB) with five to eight FLA-I and two to three FLA-DRB transcribed loci per haplotype. It is evident that the eight FLA haplotypes were generated by gene duplications and deletions, and rearrangements by genetic recombination with the accumulation and/or inheritance of novel polymorphisms. These findings are useful for further genetic diversity analysis and disease association studies among cat breeds and in veterinary medicine.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

MHC structuring and divergent allele advantage in a urodele amphibian: a hierarchical multi-scale approach

Proteins encoded by extraordinarily polymorphic major histocompatibility complex (MHC) genes are involved in the adaptive immune response. Balancing selection is believed to maintain MHC polymorphism in the long term, although neutral processes also play a role in shaping MHC diversity. However, the relative contribution of these processes is poorly understood. Here we characterized MHC class II variation of a low-dispersal, pond-breeding newt (Triturus carnifex) over a restricted, geographically structured area. We aimed to (1) evaluate the contribution of selection and neutral processes to shaping MHC diversity at two geographic scales, and (2) test for signatures of divergent allele advantage (DAA), which is a potentially important mechanism of balancing selection. The dominant role of selection in shaping MHC variation was suggested by the lack of correlation between MHC and neutral (microsatellite) variation. Although most variation occurred within populations for both types of markers, they differed in the extent of structuring at the two spatial scales. MHC structuring was more pronounced at local scales, suggesting the role of local selection, while structuring was not detectable at a larger scale, possibly due to the effect of balancing selection. Microsatellites showed the opposite pattern. As expected under DAA, the observed genotypes combined more sequence diversity than expected under a random association of alleles. Thus, DAA may contribute to maintaining MHC polymorphism, which is ancient, as supported by signatures of historical positive selection and trans-species polymorphism. Our results point to the importance of a multi-scale approach in studying MHC variation, especially in low-dispersal taxa, which are genetically structured at fine spatial scales.

A genome-wide scan for genes under balancing selection in the plant pathogen Ralstonia solanacearum

Background

Plant pathogens are under significant selective pressure by the plant host. Consequently, they are expected to have adapted to this condition or contribute to evading plant defenses. In order to acquire long-term fitness, plant bacterial pathogens are usually forced to maintain advantageous genetic diversity in populations. This strategy ensures that different alleles in the pathogen’s gene pool are maintained in a population at frequencies larger than expected under neutral evolution. This selective process, known as balancing selection, is the subject of this work in the context of a common bacterial phytopathogen. We performed a genome-wide scan of Ralstonia solanacearum species complex, an aggressive plant bacterial pathogen that shows broad host range and causes a devastating disease called ‘bacterial wilt’.

Results

Using a sliding window approach, we analyzed 57 genomes from three phylotypes of the R. solanacearum species complex to detect signatures of balancing selection. A total of 161 windows showed extreme values in three summary statistics of population genetics: Tajima’s D, θ_w and Fu & Li’s D*. We discarded any confounding effects due to demographic events by means of coalescent simulations of genetic data. The prospective windows correspond to 78 genes with known function that map in any of the two main replicons (1.7% of total number of genes). The candidate genes under balancing selection are related to primary metabolism and other basal activities (51.3%) or directly associated to virulence (48.7%), the latter being involved in key functions targeted to dismantle plant defenses or to participate in critical stages in the pathogenic process.

Conclusions

We identified various genes under balancing selection that play a significant role in basic metabolism as well as in virulence of the R. solanacearum species complex. These genes are useful to understand and monitor the evolution of bacterial pathogen populations and emerge as potential candidates for future treatments to induce specific plant immune responses.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1456-6) contains supplementary material, which is available to authorized users.

Selection, drift, and introgression shape MHC polymorphism in lizards

The major histocompatibility complex (MHC) has long served as a model for the evolution of adaptive genetic diversity in wild populations. Pathogen-mediated selection is thought to be a main driver of MHC diversity, but it remains elusive to what degree selection shapes MHC diversity in complex biogeographical scenarios where other evolutionary processes (e.g. genetic drift and introgression) may also be acting. Here we focus on two closely related green lizard species, Lacerta trilineata and L. viridis, to address the evolutionary forces acting on MHC diversity in populations with different biogeographic structure. We characterized MHC class I exon 2 and exon 3, and neutral diversity (microsatellites), to study the relative importance of selection, drift, and introgression in shaping MHC diversity. As expected, positive selection was a significant force shaping the high diversity of MHC genes in both species. Moreover, introgression significantly increased MHC diversity in mainland populations, with a primary direction of gene flow from L. viridis to L. trilineata. Finally, we found significantly fewer MHC alleles in island populations, but maintained MHC sequence and functional diversity, suggesting that positive selection counteracted the effect of drift. Overall, our data support that different evolutionary processes govern MHC diversity in different biogeographical scenarios: positive selection occurs broadly while introgression acts in sympatry and drift when the population sizes decrease.

Balancing selection and introgression of newt immune-response genes

The importance of interspecific introgression as a source of adaptive variation is increasingly recognized. Theory predicts that beneficial genetic variants cross species boundaries easily even when interspecific hybridization is rare and gene flow is strongly constrained throughout the genome. However, it remains unclear whether certain classes of genes are particularly prone to adaptive introgression. Genes affected by balancing selection (BS) may constitute such a class, because forms of BS that favour novel, initially rare alleles, should facilitate introgression. We tested this hypothesis in hybridizing newts by comparing 13 genes with signatures of BS, in particular an excess of common non-synonymous polymorphisms, to the genomic background (154 genes). Parapatric hybridizing taxa were less differentiated in BS candidate genes than more closely related allopatric lineages, while the opposite was observed in the control genes. Coalescent and forward simulations that explored neutral and BS scenarios under isolation and migration showed that processes other than differential gene flow are unlikely to account for this pattern. We conclude that BS, probably involving a form of novel allele advantage, promotes introgression. This mechanism may be a source of adaptively relevant variation in hybridizing species over prolonged periods.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.

Spatial patterns of immunogenetic and neutral variation underscore the conservation value of small, isolated American badger populations

Small and isolated populations often exhibit low genetic diversity due to drift and inbreeding, but may simultaneously harbour adaptive variation. We investigate spatial distributions of immunogenetic variation in American badger subspecies (Taxidea taxus), as a proxy for evaluating their evolutionary potential across the northern extent of the species' range. We compared genetic structure of 20 microsatellites and the major histocompatibility complex (MHC DRB exon 2) to evaluate whether small, isolated populations show low adaptive polymorphism relative to large and well-connected populations. Our results suggest that gene flow plays a prominent role in shaping MHC polymorphism across large spatial scales, while the interplay between gene flow and selection was stronger towards the northern peripheries. The similarity of MHC alleles within subspecies relative to their neutral genetic differentiation suggests that adaptive divergence among subspecies can be maintained despite ongoing gene flow along subspecies boundaries. Neutral genetic diversity was low in small relative to large populations, but MHC diversity within individuals was high in small populations. Despite reduced neutral genetic variation, small and isolated populations harbour functional variation that likely contribute to the species evolutionary potential at the northern range. Our findings suggest that conservation approaches should focus on managing adaptive variation across the species range rather than protecting subspecies per se.

Egernia stokesii (gidgee skink) MHC I positively selected sites lack concordance with HLA peptide binding regions

Genes of the major histocompatibility complex (MHC) play an important role in vertebrate disease resistance, kin recognition and mate choice. Mammalian MHC is the most widely characterised of all vertebrates, and attention is often given to the peptide binding regions of the MHC because they are presumed to be under stronger selection than non-peptide binding regions. For vertebrates where the MHC is less well understood, researchers commonly use the amino acid positions of the peptide binding regions of the human leukocyte antigen (HLA) to infer the peptide binding regions within the MHC sequences of their taxon of interest. However, positively selected sites within MHC have been reported to lack correspondence with the HLA in fish, frogs, birds and reptiles including squamates. Despite squamate diversity, the MHC has been characterised in few snakes and lizards. The Egernia group of scincid lizards is appropriate for investigating mechanisms generating MHC variation, as their inclusion will add a new lineage (i.e. Scincidae) to studies of selection on the MHC. We aimed to identify positively selected sites within the MHC of Egernia stokesii and then determine if these sites corresponded with the peptide binding regions of the HLA. Six positively selected sites were identified within E. stokesii MHC I, only two were homologous with the HLA. E. stokesii positively selected sites corresponded more closely to non-lizard than other lizard taxa. The characterisation of the MHC of more intermediate taxa within the squamate order is necessary to understand the evolution of the MHC across all vertebrates.

PCR Strategies for Complete Allele Calling in Multigene Families Using High-Throughput Sequencing Approaches

The characterization of multigene families with high copy number variation is often approached through PCR amplification with highly degenerate primers to account for all expected variants flanking the region of interest. Such an approach often introduces PCR biases that result in an unbalanced representation of targets in high-throughput sequencing libraries that eventually results in incomplete detection of the targeted alleles. Here we confirm this result and propose two different amplification strategies to alleviate this problem. The first strategy (called pooled-PCRs) targets different subsets of alleles in multiple independent PCRs using different moderately degenerate primer pairs, whereas the second approach (called pooled-primers) uses a custom-made pool of non-degenerate primers in a single PCR. We compare their performance to the common use of a single PCR with highly degenerate primers using the MHC class I of the Iberian lynx as a model. We found both novel approaches to work similarly well and better than the conventional approach. They significantly scored more alleles per individual (11.33 ± 1.38 and 11.72 ± 0.89 vs 7.94 ± 1.95), yielded more complete allelic profiles (96.28 ± 8.46 and 99.50 ± 2.12 vs 63.76 ± 15.43), and revealed more alleles at a population level (13 vs 12). Finally, we could link each allele's amplification efficiency with the primer-mismatches in its flanking sequences and show that ultra-deep coverage offered by high-throughput technologies does not fully compensate for such biases, especially as real alleles may reach lower coverage than artefacts. Adopting either of the proposed amplification methods provides the opportunity to attain more complete allelic profiles at lower coverages, improving confidence over the downstream analyses and subsequent applications.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.

MHC, parasites and antler development in red deer: no support for the Hamilton & Zuk hypothesis

The Hamilton-Zuk hypothesis proposes that the genetic benefits of preferences for elaborated secondary sexual traits have their origins in the arms race between hosts and parasites, which maintains genetic variance in parasite resistance. Infection, in turn, can be reflected in the expression of costly sexual ornaments. However, the link between immune genes, infection and the expression of secondary sexual traits has rarely been investigated. Here, we explored whether the presence and identity of functional variants (supertypes) of the highly polymorphic major histocompatibility complex (MHC), which is responsible for the recognition of parasites, predict the load of lung and gut parasites and antler development in the red deer (Cervus elaphus). While we found MHC supertypes to be associated with infection by a number of parasite species, including debilitating lung nematodes, we did not find support for the Hamilton-Zuk hypothesis. On the contrary, we found that lung nematode load was positively associated with antler development. We also found that the supertypes that were associated with resistance to certain parasites at the same time cause susceptibility to others. Such trade-offs may undermine the potential genetic benefits of mate choice for resistant partners.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Contrasting patterns of selection and drift between two categories of immune genes in prairie-chickens

Immune-receptor genes of the adaptive immune system, such as the major histocompatibility complex (MHC), are involved in recognizing specific pathogens and are known to have high rates of adaptive evolution, presumably as a consequence of rapid co-evolution between hosts and pathogens. In contrast, many 'mediating' genes of the immune system do not interact directly with specific pathogens and are involved in signalling (e.g. cytokines) or controlling immune cell growth. As a consequence, we might expect stronger selection at immune-receptor than mediating genes, but these two types of genes have not been compared directly in wild populations. Here, we tested the hypothesis that selection differs between MHC (class I and II) and mediating genes by comparing levels of population differentiation across the range of greater prairie-chickens (Tympanuchus cupido). As predicted, there was stronger population differentiation and isolation by distance at immune receptor (MHC) than at either mediating genes or neutral microsatellites, suggesting a stronger role of local adaptation at the MHC. In contrast, mediating genes displayed weaker differentiation between populations than neutral microsatellites, consistent with selection favouring similar alleles across populations for mediating genes. In addition to selection, drift also had a stronger effect on immune receptor (MHC) than mediating genes as indicated by the stronger decline of MHC variation in relation to population size. This is the first study in the wild to show that the effects of selection and drift on immune genes vary across populations depending on their functional role.

Lack of Spatial Immunogenetic Structure among Wolverine (Gulo gulo) Populations Suggestive of Broad Scale Balancing Selection

Elucidating the adaptive genetic potential of wildlife populations to environmental selective pressures is fundamental for species conservation. Genes of the major histocompatibility complex (MHC) are highly polymorphic, and play a key role in the adaptive immune response against pathogens. MHC polymorphism has been linked to balancing selection or heterogeneous selection promoting local adaptation. However, spatial patterns of MHC polymorphism are also influenced by gene flow and drift. Wolverines are highly vagile, inhabiting varied ecoregions that include boreal forest, taiga, tundra, and high alpine ecosystems. Here, we investigated the immunogenetic variation of wolverines in Canada as a surrogate for identifying local adaptation by contrasting the genetic structure at MHC relative to the structure at 11 neutral microsatellites to account for gene flow and drift. Evidence of historical positive selection was detected at MHC using maximum likelihood codon-based methods. Bayesian and multivariate cluster analyses revealed weaker population genetic differentiation at MHC relative to the increasing microsatellite genetic structure towards the eastern wolverine distribution. Mantel correlations of MHC against geographical distances showed no pattern of isolation by distance (IBD: r = -0.03, p = 0.9), whereas for microsatellites we found a relatively strong and significant IBD (r = 0.54, p = 0.01). Moreover, we found a significant correlation between microsatellite allelic richness and the mean number of MHC alleles, but we did not observe low MHC diversity in small populations. Overall these results suggest that MHC polymorphism has been influenced primarily by balancing selection and to a lesser extent by neutral processes such as genetic drift, with no clear evidence for local adaptation. This study contributes to our understanding of how vulnerable populations of wolverines may respond to selective pressures across their range.

Major histocompatibility complex similarity and sexual selection: different does not always mean attractive

Females that mate multiply have the possibility to exert postcopulatory choice and select more compatible sperm to fertilize eggs. Prior work suggests that dissimilarity in major histocompatibility complex (MHC) plays an important role in determining genetic compatibility between partners. Favouring a partner with dissimilar MHC alleles would result in offspring with high MHC diversity and therefore with enhanced survival thanks to increased resistance to pathogens and parasites. The high variability of MHC genes may further allow discrimination against the sperm from related males, reducing offspring homozygosity and inbreeding risk. Despite the large body of work conducted at precopulatory level, the role of MHC similarity between partners at postcopulatory level has been rarely investigated. We used an internal fertilizing fish with high level of multiple matings (Poecilia reticulata) to study whether MHC similarity plays a role in determining the outcome of fertilization when sperm from two males compete for the same set of eggs. We also controlled for genomewide similarity by determining similarity at 10 microsatellite loci. Contrary to prediction, we found that the more MHC-similar male sired more offspring while similarity at the microsatellite loci did not predict the outcome of sperm competition. Our results suggest that MHC discrimination may be involved in avoidance of hybridization or outbreeding rather than inbreeding avoidance. This, coupled with similar findings in salmon, suggests that the preference for MHC-dissimilar mates is far from being unanimous and that pre- and postcopulatory episodes of sexual selection can indeed act in opposite directions.

Large-scale genotyping of highly polymorphic loci by next-generation sequencing: how to overcome the challenges to reliably genotype individuals?

Studying the different roles of adaptive genes is still a challenge in evolutionary ecology and requires reliable genotyping of large numbers of individuals. Next-generation sequencing (NGS) techniques enable such large-scale sequencing, but stringent data processing is required. Here, we develop an easy to use methodology to process amplicon-based NGS data and we apply this methodology to reliably genotype four major histocompatibility complex (MHC) loci belonging to MHC class I and II of Alpine marmots (Marmota marmota). Our post-processing methodology allowed us to increase the number of retained reads. The quality of genotype assignment was further assessed using three independent validation procedures. A total of 3069 high-quality MHC genotypes were obtained at four MHC loci for 863 Alpine marmots with a genotype assignment error rate estimated as 0.21%. The proposed methodology could be applied to any genetic system and any organism, except when extensive copy-number variation occurs (that is, genes with a variable number of copies in the genotype of an individual). Our results highlight the potential of amplicon-based NGS techniques combined with adequate post-processing to obtain the large-scale highly reliable genotypes needed to understand the evolution of highly polymorphic functional genes.