Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The F_ST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

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.

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.

A novel workflow to improve genotyping of multigene families in wildlife species: An experimental set-up with a known model system

Genotyping complex multigene families in novel systems is particularly challenging. Target primers frequently amplify simultaneously multiple loci leading to high PCR and sequencing artefacts such as chimeras and allele amplification bias. Most genotyping pipelines have been validated in nonmodel systems whereby the real genotype is unknown and the generation of artefacts may be highly repeatable. Further hindering accurate genotyping, the relationship between artefacts and genotype complexity (i.e. number of alleles per genotype) within a PCR remains poorly described. Here, we investigated the latter by experimentally combining multiple known major histocompatibility complex (MHC) haplotypes of a model organism (chicken, Gallus gallus, 43 artificial genotypes with 2-13 alleles per amplicon). In addition to well-defined 'optimal' primers, we simulated a nonmodel species situation by designing 'cross-species' primers based on sequence data from closely related Galliform species. We applied a novel open-source genotyping pipeline (ACACIA; https://gitlab.com/psc_santos/ACACIA), and compared its performance with another, previously published pipeline (AmpliSAS). Allele calling accuracy was higher when using ACACIA (98.5% versus 97% and 77.8% versus 75% for the 'optimal' and 'cross-species' data sets, respectively). Systematic allele dropout of three alleles owing to primer mismatch in the 'cross-species' data set explained high allele calling repeatability (100% when using ACACIA) despite low accuracy, demonstrating that repeatability can be misleading when evaluating genotyping workflows. Genotype complexity was positively associated with nonchimeric artefacts, chimeric artefacts (nonlinearly by levelling when amplifying more than 4-6 alleles) and allele amplification bias. Our study exemplifies and demonstrates pitfalls researchers should avoid to reliably genotype complex multigene families.

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.

Allelic diversity and selection at the MHC class I and class II in a bottlenecked bird of prey, the White-tailed Eagle

Background

Genes of the Major Histocompatibility Complex (MHC) are essential for adaptive immune response in vertebrates, as they encode receptors that recognize peptides derived from the processing of intracellular (MHC class I) and extracellular (MHC class II) pathogens. High MHC diversity in natural populations is primarily generated and maintained by pathogen-mediated diversifying and balancing selection. It is, however, debated whether selection at the MHC can counterbalance the effects of drift in bottlenecked populations. The aim of this study was to assess allelic diversity of MHC genes in a recently bottlenecked bird of prey, the White-tailed Eagle Haliaeetus albicilla, as well as to compare mechanisms that shaped the evolution of MHC class I and class II in this species.

Results

We showed that significant levels of MHC diversity were retained in the core Central European (Polish) population of White-tailed Eagles. Ten MHC class I and 17 MHC class II alleles were recovered in total and individual birds showed high average MHC diversity (3.80 and 6.48 MHC class I and class II alleles per individual, respectively). Distribution of alleles within individuals provided evidence for the presence of at least three class I and five class II loci the White-tailed Eagle, which suggests recent duplication events. MHC class II showed greater sequence polymorphism than MHC class I and there was much stronger signature of diversifying selection acting on MHC class II than class I. Phylogenetic analysis provided evidence for trans-species similarity of class II, but not class I, sequences, which is likely consistent with stronger balancing selection at MHC class II.

Conclusions

Relatively high MHC diversity retained in the White-tailed Eagles from northern Poland reinforces high conservation value of local eagle populations. At the same time, our study is the first to demonstrate contrasting patterns of allelic diversity and selection at MHC class I and class II in an accipitrid species, supporting the hypothesis that different mechanisms can shape evolutionary trajectories of MHC class I and class II genes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1338-3) contains supplementary material, which is available to authorized users.

Selection on MHC in a Context of Historical Demographic Change in 2 Closely Distributed Species of Tuco-tucos (Ctenomys australis and C. talarum)

Selection necessarily acts within the same current and historical demographic framework as neutral evolutionary processes, and the outcome of the interplay between these forces may vary according to their relative strength. In this study, we compare the variation at a major histocompatibility complex (MHC) locus (DRB exon 2), typically subject to strong diversifying selection, and mitochondrial diversity (D-loop) across populations encompassing the entire distribution of 2 species of South American subterranean rodents: Ctenomys australis and C. talarum (tuco-tucos). Although these species are parapatric along most of their distribution, historically they have followed distinct demographic trajectories associated with sea level changes during the Quaternary. We surveyed 8 populations of C. australis and 15 of C. talarum, from which we analyzed 70 and 212 D-loop haplotypes and 91 and 346 DRB genotypes, respectively. Both species have gone through a recent demographic expansion; however, the signal of this process only encompasses the entire distribution of one of the species: C. australis. While balancing selection on MHC in C. talarum-enhanced DRB diversity at the local level compared to D-loop, although not promoting divergence among populations, in C. australis local diversifying selection may have driven higher population differentiation at DRB than at D-loop. Our findings reinforce the idea that the relative strength of selection acting on MHC genes varies spatially and temporally within and among species, even between species using the same macrohabitat and exposed to similar immune challenges.

Historical gene flow and profound spatial genetic structure among golden pheasant populations suggested by multi-locus analysis

Major histocompatibility complex (MHC) is a good marker system for geographical genetics since they are functional genes in the immune system that are likely to affect the fitness of the individual, and the survival and evolutionary potential of a population in a changing environment. Golden pheasant (Chrysolophus pictus) is a wild Phasianidae distributed in central and north China. In this study, we used a locus-specific genotyping technique for MHC IIB genes of golden pheasant. Combining with microsatellites (simple sequence repeat, SSR) and mitochondrial DNA (mtDNA) D-loop region, we investigated the demographic history and illuminate genetic structure of this bird in detail. SYR (south of Yangtze river) - NYR (north of Yangtze river) lineages, separated by Yangtze River, were defined in genetic structure of MHC IIB. NYR was supposed as refuge during glacial period, suggested by diversity parameters and more ancient alleles in this region. Based on this hypothesis, there was gene flow from NYR to SYR, which was proved by three pieces of evidence: (1) distinct demographic histories of SYR (kept stable) and NYR (experienced expansion); (2) specific affiliation of LC in genetic structure of SSR and MHC genes; (3) significant gene flow from NYR to SYR. Moreover, we also found balancing selection by combination of three Grouping A2's regions (SC, QL and North) into one in Grouping B4 (NYR) and no pattern of isolation by distance (IBD) found in MHC IIB, whereas for SSR we found a relatively strong and significant IBD. Several mechanisms in the evolution of MHC IIB genes, including recombination, historically positive selection, trans-species evolution and concerted evolution, were shown by molecular and phylogenetic analysis. Overall these results suggest the Yangtze River was inferred to be a geological barrier for this avian and NYR might experience population expansion, which invaded into a neighboring region. This study contributes to the understanding of the effects of geographic features on contemporary patterns of genetic variation in the golden pheasant in China, and helps us to define the adaptive unite (AU) for this avian.

Coloniality and migration are related to selection on MHC genes in birds

The major histocompatibility complex (MHC) plays a key role in pathogen recognition as a part of the vertebrate adaptive immune system. The great diversity of MHC genes in natural populations is maintained by different forms of balancing selection and its strength should correlate with the diversity of pathogens to which a population is exposed and the rate of exposure. Despite this prediction, little is known about how life-history characteristics affect selection at the MHC. Here, we examined whether the strength of balancing selection on MHC class II genes in birds (as measured with nonsynonymous nucleotide substitutions, dN) was related to their social or migratory behavior, two life-history characteristics correlated with pathogen exposure. Our comparative analysis indicated that the rate of nonsynonymous substitutions was higher in colonial and migratory species than solitary and resident species, suggesting that the strength of balancing selection increases with coloniality and migratory status. These patterns could be attributed to: (1) elevated transmission rates of pathogens in species that breed in dense aggregations, or (2) exposure to a more diverse fauna of pathogens and parasites in migratory species. Our study suggests that differences in social structure and basic ecological traits influence MHC diversity in natural vertebrate populations.

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.

Development of novel associations between MHC alleles and susceptibility to parasitic infections in an isolated population of an endangered mammal

The role of pathogens in dynamics of endangered species is not fully understood, and the effect of infection often interacts with other processes affecting those species, such as fragmentation and isolation or loss of genetic variation. Small, isolated populations are prone to losing functional alleles due to demographic processes and genetic drift, which may diminish their ability to resist infection if immune genes are affected. Demographic processes may also alter the selective pressure exerted by a parasite, as they influence the rate of parasite transmission between individuals. In the present paper we studied changes in parasite infection levels and genetic variability in an isolated population of spotted suslik (Spermophillus suslicus). Over a three-year period (approx. three generations), when the population size remained relatively stable, we observed a considerable increase in parasite prevalence and infection intensity, followed by the development of novel associations between MHC DRB alleles and parasite burden. Contrary to expectations, the change in MHC allele frequency over time was not consistent with the effect of the allele - for instance, Spsu-DRB*07, associated with higher intensity of infection with a nematode Capillaria sp., increased in frequency from 11.8 to 20.2%. Yet, we found no signatures of selection in the studied loci. Our results show that an isolated, stable population may experience a sudden increase in parasitic infections, resulting in a development of novel associations between MHC alleles and parasite susceptibility/resistance, even though no signatures of selection can be found.

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake, the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long-term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.

Genetic variation at the MHC DRB1 locus is similar across Gunnison's prairie dog (Cynomys gunnisoni) colonies regardless of plague history

Yersinia pestis was introduced to North America around 1900 and leads to nearly 100% mortality in prairie dog (Cynomys spp.) colonies during epizootic events, which suggests this pathogen may exert a strong selective force. We characterized genetic diversity at an MHC class II locus (DRB1) in Gunnison's prairie dog (C. gunnisoni) and quantified population genetic structure at the DRB1 versus 12 microsatellite loci in three large Arizona colonies. Two colonies, Seligman (SE) and Espee Ranch (ES), have experienced multiple plague‐related die‐offs in recent years, whereas plague has never been documented at Aubrey Valley (AV). We found fairly low allelic diversity at the DRB1 locus, with one allele (DRB1*01) at high frequency (0.67–0.87) in all colonies. Two other DRB1 alleles appear to be trans‐species polymorphisms shared with the black‐tailed prairie dog (C. ludovicianus), indicating that these alleles have been maintained across evolutionary time frames. Estimates of genetic differentiation were generally lower at the MHC locus (F_ST = 0.033) than at microsatellite markers (F_ST = 0.098). The reduced differentiation at DRB1 may indicate that selection has been important for shaping variation at MHC loci, regardless of the presence or absence of plague in recent decades. However, genetic drift has probably also influenced the DRB1 locus because its level of differentiation was not different from that of microsatellites in an F_ST outlier analysis. We then compared specific MHC alleles to plague survivorship in 60 C. gunnisoni that had been experimentally infected with Y. pestis. We found that survival was greater in individuals that carried at least one copy of the most common allele (DRB1*01) compared to those that did not (60% vs. 20%). Although the sample sizes of these two groups were unbalanced, this result suggests the possibility that this MHC class II locus, or a nearby linked gene, could play a role in plague survival.

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.

Spatial and temporal variation at major histocompatibility complex class IIB genes in the endangered Blakiston's fish owl

INTRODUCTION

Quantifying intraspecific genetic variation in functionally important genes, such as those of the major histocompatibility complex (MHC), is important in the establishment of conservation plans for endangered species. The MHC genes play a crucial role in the vertebrate immune system and generally show high levels of diversity, which is likely due to pathogen-driven balancing selection. The endangered Blakiston's fish owl (Bubo blakistoni) has suffered marked population declines on Hokkaido Island, Japan, during the past several decades due to human-induced habitat loss and fragmentation. We investigated the spatial and temporal patterns of genetic diversity in MHC class IIβ genes in Blakiston's fish owl, using massively parallel pyrosequencing.

RESULTS

We found that the Blakiston's fish owl genome contains at least eight MHC class IIβ loci, indicating recent gene duplications. An analysis of sequence polymorphism provided evidence that balancing selection acted in the past. The level of MHC variation, however, was low in the current fish owl populations in Hokkaido: only 19 alleles were identified from 174 individuals. We detected considerable spatial differences in MHC diversity among the geographically isolated populations. We also detected a decline of MHC diversity in some local populations during the past decades.

CONCLUSIONS

Our study demonstrated that the current spatial patterns of MHC variation in Blakiston's fish owl populations have been shaped by loss of variation due to the decline and fragmentation of populations, and that the short-term effects of genetic drift have counteracted the long-term effects of balancing selection.

Characterization and genotyping of the DRB1 gene of the major histocompatibility complex (MHC) in the Marmota monax, animal model of hepatitis B

The major histocompatibility complex (MHC)-containing genes are among the most polymorphic in vertebrates. MHC genes code for proteins that are critical in the immune system response. In this study, the polymorphism of the second exon of the MHC class II DRB gene was characterized in the Eastern woodchuck (Marmota monax). Woodchucks chronically infected with the woodchuck hepatitis virus (WHV) represent the best available animal model for the study of chronic hepatitis B infection in humans. In the genotyped animals we found fifteen alleles, which were expressed in two independent loci and that were named DRB1A and DRB1B in this work. The 15 alleles investigated showed an elevated divergence. A significant excess of non-synonymous substitutions was detected, which could indicate that a historical positive selection is acting in the woodchuck DRB1 genes. This hypothesis was confirmed in our study by the high variability in or near the antigen binding sites (ABS) and by the results obtained in sequence variability analyses. This analysis identified the presence of a microsatellite sequence that is located at the start of the second intron, which could further allow the development of a fast and cheap semiautomatic sequencing method.

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

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre-eminent system for the study of selective pressures that arise from host-pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population-genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.

Spatial patterns of neutral and functional genetic variations reveal patterns of local adaptation in raccoon (Procyon lotor) populations exposed to raccoon rabies

Local adaptation is necessary for population survival and depends on the interplay between responses to selective forces and demographic processes that introduce or retain adaptive and maladaptive attributes. Host-parasite systems are dynamic, varying in space and time, where both host and parasites must adapt to their ever-changing environment in order to survive. We investigated patterns of local adaptation in raccoon populations with varying temporal exposure to the raccoon rabies virus (RRV). RRV infects approximately 85% of the population when epizootic and has been presumed to be completely lethal once contracted; however, disease challenge experiments and varying spatial patterns of RRV spread suggest some level of immunity may exist. We first assessed patterns of local adaptation in raccoon populations along the eastern seaboard of North America by contrasting spatial patterns of neutral (microsatellite loci) and functional, major histocompatibility complex (MHC) genetic diversity and structure. We explored variation of MHC allele frequencies in the light of temporal population exposure to RRV (0-60 years) and specific RRV strains in infected raccoons. Our results revealed high levels of MHC variation (66 DRB exon 2 alleles) and pronounced genetic structure relative to neutral microsatellite loci, indicative of local adaptation. We found a positive association linking MHC genetic diversity and temporal RRV exposure, but no association with susceptibility and resistance to RRV strains. These results have implications for landscape epidemiology studies seeking to predict the spread of RRV and present an example of how population demographics influence the degree to which populations adapt to local selective pressures.

Genomic Porosity between Invasive Chondrostoma nasus and Endangered Endemic Parachondrostoma toxostoma (Cyprinidae): The Evolution of MHC IIB Genes

Two cyprinid species, Parachondrostoma toxostoma, an endemic threatened species, and Chondrostoma nasus, an invasive species, live in sympatry in southern France and form two sympatric zones where the presence of intergeneric hybrids is reported. To estimate the potential threat to endemic species linked to the introduction of invasive species, we focused on the DAB genes (functional MHC IIB genes) because of their adaptive significance and role in parasite resistance. More specifically, we investigated (1) the variability of MHC IIB genes, (2) the selection pattern shaping MHC polymorphism, and (3) the extent to which trans-species evolution and intergeneric hybridization affect MHC polymorphism. In sympatric areas, the native species has more diversified MHC IIB genes when compared to the invasive species, probably resulting from the different origins and dispersal of both species. A similar level of MHC polymorphism was found at population level in both species, suggesting similar mechanisms generating MHC diversity. In contrast, a higher number of DAB-like alleles per specimen were found in invasive species. Invasive species tended to express the alleles of two DAB lineages, whilst native species tended to express the alleles of only the DAB3 lineage. Hybrids have a pattern of MHC expression intermediate between both species. Whilst positive selection acting on peptide binding sites (PBS) was demonstrated in both species, a slightly higher number of positively selected sites were identified in C. nasus, which could result from parasite-mediated selection. Bayesian clustering analysis revealed a similar pattern of structuring for the genetic variation when using microsatellites or the MHC approach. We confirmed the importance of trans-species evolution for MHC polymorphism. In addition, we demonstrated bidirectional gene flow for MHC IIB genes in sympatric areas. The positive significant correlation between MHC and microsatellites suggests that demographic factors may contribute to MHC variation on a short time scale.

Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana)

BACKGROUND

Small, isolated populations often experience loss of genetic variation due to random genetic drift. Unlike neutral or nearly neutral markers (such as mitochondrial genes or microsatellites), major histocompatibility complex (MHC) genes in these populations may retain high levels of polymorphism due to balancing selection. The relative roles of balancing selection and genetic drift in either small isolated or bottlenecked populations remain controversial. In this study, we examined the mechanisms maintaining polymorphisms of MHC genes in small isolated populations of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) by comparing genetic variation found in MHC and microsatellite loci. There are few studies of this kind conducted on highly endangered primate species.

RESULTS

Two MHC genes were sequenced and sixteen microsatellite loci were genotyped from samples representing three isolated populations. We isolated nine DQA1 alleles and sixteen DQB1 alleles and validated expression of the alleles. Lowest genetic variation for both MHC and microsatellites was found in the Shennongjia (SNJ) population. Historical balancing selection was revealed at both the DQA1 and DQB1 loci, as revealed by excess non-synonymous substitutions at antigen binding sites (ABS) and maximum-likelihood-based random-site models. Patterns of microsatellite variation revealed population structure. FST outlier analysis showed that population differentiation at the two MHC loci was similar to the microsatellite loci.

CONCLUSIONS

MHC genes and microsatellite loci showed the same allelic richness pattern with the lowest genetic variation occurring in SNJ, suggesting that genetic drift played a prominent role in these isolated populations. As MHC genes are subject to selective pressures, the maintenance of genetic variation is of particular interest in small, long-isolated populations. The results of this study may contribute to captive breeding and translocation programs for endangered species.

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.

Geographic variation of the major histocompatibility complex in Eastern Atlantic grey seals (Halichoerus grypus)

Pathogen-driven balancing selection maintains high genetic diversity in many vertebrates, particularly in the major histocompatibility complex (MHC) immune system gene family, which is often associated with disease susceptibility. In large natural populations where subpopulations face different pathogen pressures, the MHC should show greater genetic differentiation within a species than neutral markers. We examined genetic diversity at the MHC-DQB locus and nine putatively neutral microsatellite markers in grey seals (Halichoerus grypus) from eight United Kingdom (UK) colonies, the Faeroe Islands and Sable Island, Canada. Five DQB alleles were identified in grey seals, which varied in prevalence across the grey seal range. Among the seal colonies, significant differences in DQB allele and haplotype frequencies and in average DQB heterozygosity were observed. Additionally, the DQB gene exhibited greater differentiation among colonies compared with neutral markers, yet a weaker pattern of isolation by distance (IBD). After correcting for the underlying IBD pattern, subpopulations breeding in similar habitats were more similar to one another in DQB allele frequencies than populations breeding in different habitats, but the same did not hold true for microsatellites, suggesting that habitat-specific pathogen pressure influences MHC evolution. Overall, the data are consistent with selection at MHC-DQB loci in grey seals with both varying selective pressures and geographic population structure appearing to influence the DQB genetic composition of breeding colonies.

MHC heterozygosity and survival in red junglefowl

Genes of the major histocompatibility complex (MHC) form a vital part of the vertebrate immune system and play a major role in pathogen resistance. The extremely high levels of polymorphism observed at the MHC are hypothesised to be driven by pathogen-mediated selection. Although the exact nature of selection remains unclear, three main hypotheses have been put forward; heterozygote advantage, negative frequency-dependence and fluctuating selection. Here, we report the effects of MHC genotype on survival in a cohort of semi-natural red junglefowl (Gallus gallus) that suffered severe mortality as a result of an outbreak of the disease coccidiosis. The cohort was followed from hatching until 250 days of age, approximately the age of sexual maturity in this species, during which time over 80% of the birds died. We show that on average birds with MHC heterozygote genotypes survived infection longer than homozygotes and that this effect was independent of genome-wide heterozygosity, estimated across microsatellite loci. This MHC effect appeared to be caused by a single susceptible haplotype (CD_c) the effect of which was masked in all heterozygote genotypes by other dominant haplotypes. The CD_c homozygous genotype had lower survival than all other genotypes, but CD_c heterozygous genotypes had survival probabilities equal to the most resistant homozygote genotype. Importantly, no heterozygotes conferred greater resistance than the most resistant homozygote genotype, indicating that the observed survival advantage of MHC heterozygotes was the product of dominant, rather than overdominant processes. This pattern and effect of MHC diversity in our population could reflect the processes ongoing in similarly small, fragmented natural populations.

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well-known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.

Characterization of major histocompatibility complex (MHC) DRB exon 2 and DRA exon 3 fragments in a primary terrestrial rabies vector (Procyon lotor)

The major histocompatibility complex (MHC) presents a unique system to explore links between genetic diversity and pathogens, as diversity within MHC is maintained in part by pathogen driven selection. While the majority of wildlife MHC studies have investigated species that are of conservation concern, here we characterize MHC variation in a common and broadly distributed species, the North American raccoon (Procyon lotor). Raccoons host an array of broadly distributed wildlife diseases (e.g., canine distemper, parvovirus and raccoon rabies virus) and present important human health risks as they persist in high densities and in close proximity to humans and livestock. To further explore how genetic variation influences the spread and maintenance of disease in raccoons we characterized a fragment of MHC class II DRA exon 3 (250bp) and DRB exon 2 (228 bp). MHC DRA was found to be functionally monomorphic in the 32 individuals screened; whereas DRB exon 2 revealed 66 unique alleles among the 246 individuals screened. Between two and four alleles were observed in each individual suggesting we were amplifying a duplicated DRB locus. Nucleotide differences between DRB alleles ranged from 1 to 36 bp (0.4–15.8% divergence) and translated into 1 to 21 (1.3–27.6% divergence) amino acid differences. We detected a significant excess of nonsynonymous substitutions at the peptide binding region (P = 0.005), indicating that DRB exon 2 in raccoons has been influenced by positive selection. These data will form the basis of continued analyses into the spatial and temporal relationship of the raccoon rabies virus and the immunogenetic response in its primary host.

Does reduced MHC diversity decrease viability of vertebrate populations?

Loss of genetic variation may render populations more vulnerable to pathogens due to inbreeding depression and depletion of variation in genes responsible for immunity against parasites. Here we review the evidence for the significance of variation in genes of the Major Histocompatibility Complex (MHC) for conservation efforts. MHC molecules present pathogen-derived antigens to the effector cells of the immune system and thus trigger the adaptive immune response. Some MHC genes are the most variable functional genes in the vertebrate genome. Their variation is clearly of adaptive significance and there is considerable evidence that its maintenance is mainly due to balancing selection imposed by pathogens. However, while the evidence for selection shaping MHC variation on the historical timescale is compelling, a correlation between levels of MHC variation and variation at neutral loci is often observed, indicating that on a shorter timescale drift also substantially affects MHC, leading to depletion of MHC diversity. The evidence that the loss of MHC variation negatively affects population survival is so far equivocal and difficult to separate from effects of general inbreeding. Some species with depleted MHC variation seem to be particularly susceptible to infection, but other species thrive and expand following severe bottlenecks that have drastically limited their MHC variation. However, while the latter demonstrate that MHC variation is not always critical for population survival, these species may in fact represent rare examples of survival despite of the loss of MHC variation. There is clearly a compelling need for data that would disclose the possible consequences of MHC diversity for population viability. In particular, we need more data on the impact of MHC allelic richness on the abundance of parasites or prevalence of disease in populations, while controlling for the role of general inbreeding. Before such evidence accumulates, captive breeding programs and other conservation measures aimed at inbreeding avoidance should be favoured over those protecting only MHC variation, especially since inbreeding avoidance programs would usually conserve both types of genetic diversity simultaneously.

How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings

Major histocompatibility complex (MHC) genes have been put forward as a model for studying how genetic diversity is maintained in wild populations. Pathogen-mediated selection (PMS) is believed to generate the extraordinary levels of MHC diversity observed. However, establishing the relative importance of the three proposed mechanisms of PMS (heterozygote advantage, rare-allele advantage and fluctuating selection) has proved extremely difficult. Studies have attempted to differentiate between mechanisms of PMS using two approaches: (i) comparing MHC diversity with that expected under neutrality and (ii) relating MHC diversity to pathogen regime. Here, we show that in many cases the same predictions arise from the different mechanisms under these approaches, and that most studies that have inferred one mechanism of selection have not fully considered the alternative explanations. We argue that, while it may be possible to demonstrate that particular mechanisms of PMS are occurring, resolving their relative importance within a system is probably impossible. A more realistic target is to continue to demonstrate when and where the different mechanisms of PMS occur, with the aim of determining their relative importance across systems. We put forward what we believe to be the most promising approaches that will allow us to progress towards achieving this.