MHC class II genes in European wolves: a comparison with dogs

Novel insights on aquatic mammal MHC evolution: Evidence from manatee DQB diversity

The low diversity in marine mammal major histocompatibility complex (MHC) appears to support the hypothesis of reduced pathogen selective pressure in aquatic systems compared to terrestrial environments. However, the lack of characterization of the aquatic and evolutionarily distant Sirenia precludes drawing more generalized conclusions. Therefore, we aimed to characterize the MHC DQB diversity of two manatee species and compare it with those reported for marine mammals. Our results identified 12 and 6 alleles in T. inunguis and T. manatus, respectively. Alleles show high rates of nonsynonymous substitutions, suggesting loci are evolving under positive selection. Among aquatic mammals, Pinnipeda DQB had smaller numbers of alleles, higher synonymous substitution rate, and a dN/dS ratio closer to 1, suggesting it may be evolving under more relaxed selection compared to fully aquatic mammals. This contradicts one of the predictions of the hypothesis that aquatic environments impose reduced pathogen pressure to mammalian immune system. These results suggest that the unique evolutionary trajectories of mammalian MHC may impose challenges in drawing ecoevolutionary conclusions from comparisons across distant vertebrate lineages.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Whole-genome resequencing of temporally stratified samples reveals substantial loss of haplotype diversity in the highly inbred Scandinavian wolf population

Genetic drift can dramatically change allele frequencies in small populations and lead to reduced levels of genetic diversity, including loss of segregating variants. However, there is a shortage of quantitative studies of how genetic diversity changes over time in natural populations, especially on genome-wide scales. Here, we analyzed whole-genome sequences from 76 wolves of a highly inbred Scandinavian population, founded by only one female and two males, sampled over a period of 30 yr. We obtained chromosome-level haplotypes of all three founders and found that 10%–24% of their diploid genomes had become lost after about 20 yr of inbreeding (which approximately corresponds to five generations). Lost haplotypes spanned large genomic regions, as expected from the amount of recombination during this limited time period. Altogether, 160,000 SNP alleles became lost from the population, which may include adaptive variants as well as wild-type alleles masking recessively deleterious alleles. Although not sampled, we could indirectly infer that the two male founders had megabase-sized runs of homozygosity and that all three founders showed significant haplotype sharing, meaning that there were on average only 4.2 unique haplotypes in the six copies of each autosome that the founders brought into the population. This violates the assumption of unrelated founder haplotypes often made in conservation and management of endangered species. Our study provides a novel view of how whole-genome resequencing of temporally stratified samples can be used to visualize and directly quantify the consequences of genetic drift in a small inbred population.

Haplotype structures and polymorphisms of dog leukocyte antigen (DLA) class I loci shaped by intralocus and interlocus recombination events

The dog leukocyte antigen (DLA) class I genomic region is located on chromosome 12, and the class I genomic region is composed of at least two distinct haplotypic gene structures, DLA-88-DLA-12 and DLA-88-DLA-88L. However, detailed information of the genomic differences among DLA-88, DLA-12, and DLA-88L are still lacking at the full-length gene level, and therefore, DLA allelic sequences classified for each of these loci are limited in number so far. In this study, we determined the DNA sequence of a 95-kb DLA class I genomic region including DLA-88, DLA-12/88L, and DLA-64 with three DLA homozygous dogs and of 37 full-length allelic gene sequences for DLA-88 and DLA-12/88L loci in 26 DLA class I homozygous dogs. Nucleotide diversity profiles of the 95-kb regions and sequence identity scores of the allelic sequences suggested that DLA-88L is a hybrid gene generated by interlocus and/or intralocus gene conversion between DLA-88 and DLA-12. The putative minimum conversion tract was estimated to be at least an 850-bp segment in length located from the 5´flanking untranslated region to the end of intron 2. In addition, at least one DLA-12 allele (DLA-12*004:01) was newly generated by interlocus gene conversion. In conclusion, the analysis for the occurrence of gene conversion within the dog DLA class I region revealed intralocus gene conversion tracts in 17 of 27 DLA-88 alleles and two of 10 DLA-12 alleles, suggesting that intralocus gene conversion has played an important role in expanding DLA allelic variations.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

Appeasing Pheromones for the Management of Stress and Aggression during Conservation of Wild Canids: Could the Solution Be Right under Our Nose?

Simple Summary

Many canid species are declining globally. It is important to conserve these species that often serve as important predators within ecosystems. Continued human expansion and the resulting habitat fragmentation necessitate conservation interventions, such as translocation, artificial pack formation, and captive breeding programs. However, chronic stress often occurs during these actions, and can result in aggression, and the physiological suppression of immunity and reproduction. Limited options are currently available for stress and aggression management in wild canids. Pheromones provide a promising natural alternative for stress management; an appeasing pheromone has been identified for multiple domestic species and may reduce stress and aggression behaviours. Many pheromones are species-specific, and the appeasing pheromone has been found to have slight compositional changes across species. In this review, the benefits of a dog appeasing pheromone and the need to investigate species-specific derivatives to produce more pronounced and beneficial behavioural and physiological modulation in target species as a conservation tool are examined.

Abstract

Thirty-six species of canid exist globally, two are classified as critically endangered, three as endangered, and five as near threatened. Human expansion and the coinciding habitat fragmentation necessitate conservation interventions to mitigate concurrent population deterioration. The current conservation management of wild canids includes animal translocation and artificial pack formation. These actions often cause chronic stress, leading to increased aggression and the suppression of the immune and reproductive systems. Castration and pharmaceutical treatments are currently used to reduce stress and aggression in domestic and captive canids. The undesirable side effects make such treatments inadvisable during conservation management of wild canids. Pheromones are naturally occurring chemical messages that modulate behaviour between conspecifics; as such, they offer a natural alternative for behaviour modification. Animals are able to distinguish between pheromones of closely related species through small compositional differences but are more likely to have greater responses to pheromones from individuals of the same species. Appeasing pheromones have been found to reduce stress- and aggression-related behaviours in domestic species, including dogs. Preliminary evidence suggests that dog appeasing pheromones (DAP) may be effective in wild canids. However, the identification and testing of species-specific derivatives could produce more pronounced and beneficial behavioural and physiological changes in target species. In turn, this could provide a valuable tool to improve the conservation management of many endangered wild canids.

Characterization of two MHC II genes (DOB, DRB) in white-tailed deer (Odocoileus virginianus)

BACKGROUND

The major histocompatibility complex (MHC) is responsible for detecting and addressing foreign pathogens inside the body. While the general structure of MHC genes is relatively well conserved among mammalian species, it is notably different among ruminants due to a chromosomal inversion that splits MHC type II genes into two subregions (IIa, IIb). Recombination rates are reportedly high between these subregions, and a lack of linkage has been documented in domestic ruminants. However, no study has yet examined the degree of linkage between these subregions in a wild ruminant. The white-tailed deer (Odocoileus virginianus), a popular ruminant of the Cervidae family, is habitually plagued by pathogens in its natural environment (e.g. Haemonchus contortus, Elaeophora). Due to the association between MHC haplotypes and disease susceptibility, a deeper understanding of MHC polymorphism and linkage between MHC genes can further aid in this species' successful management. We sequenced MHC-DRB exon 2 (IIa) and MHC-DOB exon 2 (IIb) on the MiSeq platform from an enclosed white-tailed deer population located in Alabama.

RESULTS

We identified 12 new MHC-DRB alleles, and resampled 7 alleles, which along with other published alleles brings the total number of documented alleles in white-tailed deer to 30 for MHC-DRB exon 2. The first examination of MHC-DOB in white-tailed deer found significantly less polymorphism (11 alleles), as was expected of a non-classical MHC gene. While MHC-DRB was found to be under positive, diversifying selection, MHC-DOB was found to be under purifying selection for white-tailed deer. We found no significant linkage disequilibrium between MHC-DRB and MHC-DOB, suggesting that these loci are unlikely to be closely linked.

CONCLUSIONS

Overall, this study identified 12 new MHC-DRB exon 2 alleles and characterized a new, non-classical, MHC II gene (MHC-DOB) for white-tailed deer. We also found a lack of significant linkage between these two loci, which supports previous findings of a chromosomal inversion within the MHC type II gene region in ruminants, and suggests that white-tailed deer may have a recombination hotspot between these MHC regions similar to that found for Bos taurus.

Molecular evolution of MHC class II DRB exon 2 in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Raccoon dogs, Nyctereutes procyonoides, are native to East Asia, but have been introduced into western Russia and eastern Europe. To determine allelic diversity and elucidate the evolution of major histocompatibility complex (MHC) class II genes in the raccoon dog, we analysed a 237-bp region of DRB exon 2 from 36 individuals of native and introduced populations from Japan and Russia. We detected 23 DRB alleles (Nypr-DRBs), 22 of which were novel. Some alleles were found across the species’ range, while others were geographically restricted. For both native and introduced populations, the ratio of non-synonymous to synonymous substitution rates for codons at predicted antigen-binding sites was significantly greater than 2, indicating that Nypr-DRBs have evolved under positive selection. Mixed effect model evolution analysis and an algorithm to detect recombination showed five positively selected codons and one recombination breakpoint, respectively. Overall, our results suggest that the diversity of MHC class II DRB in N. procyonoides was influenced and maintained by recombination, pathogen-driven positive selection, geographical barriers and the founder effect. A Bayesian phylogenetic tree revealed no evidence of trans-species polymorphism (TSP), but instead showed monophyly for the Nypr-DRB alleles within a larger clade of canid sequences. The lack of TSP may have been due to long-term divergence of raccoon dogs from other canids, or to their having encountered different sets of pathogens due to occupying a different ecological niche.

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia-one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty-three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population-for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human-induced changes to the coastal ecosystem it inhabits.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore

Background

Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci.

Results

Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection.

Conclusion

Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Electronic supplementary material

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

Diversified Anchoring Features the Peptide Presentation of DLA-88*50801: First Structural Insight into Domestic Dog MHC Class I

Canines represent a crucial animal model for studying human diseases and organ transplantation, as well as the evolution of domestic animals. MHCs, with a central role in cellular immunity, are commonly used in the study of dog population genetics and genome evolution. However, the molecular basis for the peptide presentation of dog MHC remains largely unknown. In this study, peptide presentation by canine MHC class I DLA-88*50801 was structurally determined, revealing diversified anchoring modes of the binding peptides. Flexible and large pockets composed of both hydrophobic and hydrophilic residues can accommodate pathogen-derived peptides with diverse anchor residues, as confirmed by thermostability measurements. Furthermore, DLA-88*50801 contains an unusual α2 helix with a large coil in the TCR contact region. These results further our understanding of canine T cell immunity through peptide presentation of MHC class I and shed light on the molecular basis for vaccine development for canine infectious diseases, for example, canine distemper virus.

Spatial pattern of adaptive and neutral genetic diversity across different biomes in the lesser anteater (Tamandua tetradactyla)

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and activation of the adaptive immune response and are thought to be regulated by natural selection, especially due to pathogen‐driven selective pressure. In this study, we investigated the spatial distribution of MHC class IIDRB exon 2 gene diversity of the lesser anteater (Tamandua tetradactyla) across five Brazilian biomes using next‐generation sequencing and compared the MHC pattern with that of neutral markers (microsatellites). We found a noticeable high level of diversity in DRB (60 amino acid alleles in 65 individuals) and clear signatures of historical positive selection acting on this gene. Higher allelic richness and proportion of private alleles were found in rain forest biomes, especially Amazon forest, a megadiverse biome, possibly harboring greater pathogen richness as well. Neutral markers, however, showed a similar pattern to DRB, demonstrating the strength of demography as an additional force to pathogen‐driven selection in shaping MHC diversity and structure. This is the first characterization and description of diversity of a MHC gene for any member of the magna‐order Xenarthra, one of the basal lineages of placental mammals.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.

Major histocompatibility complex class II polymorphisms in forest musk deer (Moschus berezovskii) and their probable association with purulent disease

Genes of the major histocompatibility complex (MHC) family are crucial in immune responses because they present pathogenic peptides to T cells. In this study, we analysed the genetic variation in forest musk deer (Moschus berezovskii) MHC II genes and its potential association with musk deer purulent disease. In total, 53 purulent disease-susceptible and 46 purulent disease-resistant individuals were selected for MHC II exon 2 fragment analysis. Among them, 16 DQ alleles and four additional DR alleles were identified, with DQ exon 2 fragments displaying a low level of polymorphism. The nonsynonymous substitutions exceeded the synonymous substitutions in the peptide-binding sites of DQA2, DQB1 and DQB2. Then, 28 MHC II alleles were used to analyse the distribution patterns of purulent disease between the susceptible and resistant groups. Among them, three alleles (DQA1*01, DQA1*02 and DQA2*04) were found to be resistant, and five alleles (DRB3*07, DQA1*03, DQA1*04, DQA2*05 and DQA2*06) were found to increase susceptibility. Additionally, three haplotypes were found to be putatively associated with musk deer purulent disease. However, these three haplotypes were only found in the resistant or susceptible group, and their frequencies were low. The results from our study support a contributory role of MHC II polymorphisms in the development of purulent disease in forest musk deer.

MHC variability in an isolated wolf population in Italy

Small, isolated populations may experience increased extinction risk due to reduced genetic variability at important functional genes, thus decreasing the population's adaptive potential. The major histocompatibility complex (MHC), a key immunological gene cluster, usually shows high variability maintained by positive or balancing selection in response to challenges by pathogens. Here we investigated for the first time, the variability of 3 MHC class II genes (DRB1, DQA1, and DQB1) in 94 samples collected from Italian wolves. The Italian wolf population has been long isolated south of the Alps and is presently recovering from a recent bottleneck that decreased the population to less than 100 individuals. Despite the bottleneck, Italian wolves show remarkable MHC variability with 6-9 alleles per locus, including 2 recently described alleles at DRB1. MHC sequences show signatures of historical selective pressures (high d N/d S ratio, ω > 1.74) but no evidence of ongoing selection. Variation at the MHC genes and 12 background microsatellite loci were not apparently affected by the recent bottleneck. Although MHC alleles of domestic dog origin were detected in 8 genetically admixed individuals, these alleles were rare or absent in nonadmixed wolves. Thus, despite known hybridization events between domestic dogs and Italian wolves, the Italian wolf population does not appear affected by deep introgression of domestic dog MHC alleles.

Evolution of MHC class I genes in the endangered loggerhead sea turtle (Caretta caretta) revealed by 454 amplicon sequencing

Background

In evolutionary and conservation biology, parasitism is often highlighted as a major selective pressure. To fight against parasites and pathogens, genetic diversity of the immune genes of the major histocompatibility complex (MHC) are particularly important. However, the extensive degree of polymorphism observed in these genes makes it difficult to conduct thorough population screenings.

Methods

We utilized a genotyping protocol that uses 454 amplicon sequencing to characterize the MHC class I in the endangered loggerhead sea turtle (Caretta caretta) and to investigate their evolution at multiple relevant levels of organization.

Results

MHC class I genes revealed signatures of trans-species polymorphism across several reptile species. In the studied loggerhead turtle individuals, it results in the maintenance of two ancient allelic lineages. We also found that individuals carrying an intermediate number of MHC class I alleles are larger than those with either a low or high number of alleles.

Conclusions

Multiple modes of evolution seem to maintain MHC diversity in the loggerhead turtles, with relatively high polymorphism for an endangered species.

Domestication does not narrow MHC diversity in Sus scrofa

The Major Histocompatibility Complex (MHC) is a multigene family of outstanding polymorphism. MHC molecules bind antigenic peptides in the peptide-binding region (PBR) that consists of five binding pockets (P). In this study, we compared the genetic diversity of domestic pigs to that of the modern representatives of their wild ancestors, the wild boar, in two MHC loci, the oligomorphic DQA and the polymorphic DRB1. MHC nucleotide polymorphism was compared with the actual functional polymorphism in the PBR and the binding pockets P1, P4, P6, P7, and P9. The analysis of approximately 200 wild boars collected throughout Europe and 120 domestic pigs from four breeds (three pureblood, Pietrain, Leicoma, and Landrace, and one mixed Danbred) revealed that wild boars and domestic pigs share the same levels of nucleotide and amino acid polymorphism, allelic richness, and heterozygosity. Domestication did not appear to act as a bottleneck that would narrow MHC diversity. Although the pattern of polymorphism was uniform between the two loci, the magnitude of polymorphism was different. For both loci, most of the polymorphism was located in the PBR region and the presence of positive selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous substitutions in the PBR. P4 and P6 were the most polymorphic binding pockets. Functional polymorphism, i.e., the number and the distribution of pocket variants within and among populations, was significantly narrower than genetic polymorphism, indicative of a hierarchical action of selection pressures on MHC loci.

Extensive polymorphism and evidence of selection pressure on major histocompatibility complex DLA-DRB1, DQA1 and DQB1 class II genes in Croatian grey wolves

The genes of the major histocompatibility complex (MHC) are a key component of the mammalian immune system and have become important molecular markers for measuring fitness-related genetic variation in wildlife populations. Because of human persecution and habitat fragmentation, the grey wolf has become extinct from a large part of Western and Central Europe, and remaining populations have become isolated. In Croatia, the grey wolf population, part of the Dinaric-Balkan population, shrank nearly to extinction during the 20th century, and is now legally protected. Using the cloning-sequencing method, we investigated the genetic diversity and evolutionary history of exon 2 of MHC class II DLA-DRB1, DQA1 and DQB1 genes in 77 individuals. We identified 13 DRB1, 7 DQA1 and 11 DQB1 highly divergent alleles, and 13 DLA-DRB1/DQA1/DQB1 haplotypes. Selection analysis comparing the relative rates of non-synonymous to synonymous mutations (d(N)/d(S)) showed evidence of positive selection pressure acting on all three loci. Trans-species polymorphism was found, suggesting the existence of balancing selection. Evolutionary codon models detected considerable difference between alpha and beta chain gene selection patterns: DRB1 and DQB1 appeared to be under stronger selection pressure, while DQA1 showed signs of moderate selection. Our results suggest that, despite the recent contraction of the Croatian wolf population, genetic variability in selectively maintained immune genes has been preserved.

MHC variability supports dog domestication from a large number of wolves: high diversity in Asia

The process of dog domestication is still somewhat unresolved. Earlier studies indicate that domestic dogs from all over the world have a common origin in Asia. So far, major histocompatibility complex (MHC) diversity has not been studied in detail in Asian dogs, although high levels of genetic diversity are expected at the domestication locality. We sequenced the second exon of the canine MHC gene DLA-DRB1 from 128 Asian dogs and compared our data with a previously published large data set of MHC alleles, mostly from European dogs. Our results show that Asian dogs have a higher MHC diversity than European dogs. We also estimated that there is only a small probability that new alleles have arisen by mutation since domestication. Based on the assumption that all of the currently known 102 DLA-DRB1 alleles come from the founding wolf population, we simulated the number of founding wolf individuals. Our simulations indicate an effective population size of at least 500 founding wolves, suggesting that the founding wolf population was large or that backcrossing has taken place.

Impact of historical founder effects and a recent bottleneck on MHC variability in Commander Arctic foxes (Vulpes lagopus)

Populations of Arctic foxes (Vulpes lagopus) have been isolated on two of the Commander Islands (Bering and Mednyi) from the circumpolar distributed mainland population since the Pleistocene. In 1970-1980, an epizootic outbreak of mange caused a severe population decline on Mednyi Island. Genes of the major histocompatibility complex (MHC) play a primary role in infectious disease resistance. The main objectives of our study were to compare contemporary variation of MHC class II in mainland and island Arctic foxes, and to document the effects of the isolation and the recent bottleneck on MHC polymorphism by analyzing samples from historical and contemporary Arctic foxes. In 184 individuals, we found 25 unique MHC class II DRB and DQB alleles, and identified evidence of balancing selection maintaining allelic lineages over time at both loci. Twenty different MHC alleles were observed in mainland foxes and eight in Bering Island foxes. The historical Mednyi population contained five alleles and all contemporary individuals were monomorphic at both DRB and DQB. Our data indicate that despite positive and diversifying selection leading to elevated rates of amino acid replacement in functionally important antigen-binding sites, below a certain population size, balancing selection may not be strong enough to maintain genetic diversity in functionally important genes. This may have important fitness consequences and might explain the high pathogen susceptibility in some island populations. This is the first study that compares MHC diversity before and after a bottleneck in a wild canid population using DNA from museum samples.

Major Histocompatibility Complex (MHC) markers in conservation biology

Human impacts through habitat destruction, introduction of invasive species and climate change are increasing the number of species threatened with extinction. Decreases in population size simultaneously lead to reductions in genetic diversity, ultimately reducing the ability of populations to adapt to a changing environment. In this way, loss of genetic polymorphism is linked with extinction risk. Recent advances in sequencing technologies mean that obtaining measures of genetic diversity at functionally important genes is within reach for conservation programs. A key region of the genome that should be targeted for population genetic studies is the Major Histocompatibility Complex (MHC). MHC genes, found in all jawed vertebrates, are the most polymorphic genes in vertebrate genomes. They play key roles in immune function via immune-recognition and -surveillance and host-parasite interaction. Therefore, measuring levels of polymorphism at these genes can provide indirect measures of the immunological fitness of populations. The MHC has also been linked with mate-choice and pregnancy outcomes and has application for improving mating success in captive breeding programs. The recent discovery that genetic diversity at MHC genes may protect against the spread of contagious cancers provides an added impetus for managing and protecting MHC diversity in wild populations. Here we review the field and focus on the successful applications of MHC-typing for conservation management. We emphasize the importance of using MHC markers when planning and executing wildlife rescue and conservation programs but stress that this should not be done to the detriment of genome-wide diversity.

Adaptive molecular evolution of the Major Histocompatibility Complex genes, DRA and DQA, in the genus Equus

BACKGROUND

Major Histocompatibility Complex (MHC) genes are central to vertebrate immune response and are believed to be under balancing selection by pathogens. This hypothesis has been supported by observations of extremely high polymorphism, elevated nonsynonymous to synonymous base pair substitution rates and trans-species polymorphisms at these loci. In equids, the organization and variability of this gene family has been described, however the full extent of diversity and selection is unknown. As selection is not expected to act uniformly on a functional gene, maximum likelihood codon-based models of selection that allow heterogeneity in selection across codon positions can be valuable for examining MHC gene evolution and the molecular basis for species adaptations.

RESULTS

We investigated the evolution of two class II MHC genes of the Equine Lymphocyte Antigen (ELA), DRA and DQA, in the genus Equus with the addition of novel alleles identified in plains zebra (E. quagga, formerly E. burchelli). We found that both genes exhibited a high degree of polymorphism and inter-specific sharing of allele lineages. To our knowledge, DRA allelic diversity was discovered to be higher than has ever been observed in vertebrates. Evidence was also found to support a duplication of the DQA locus. Selection analyses, evaluated in terms of relative rates of nonsynonymous to synonymous mutations (dN/dS) averaged over the gene region, indicated that the majority of codon sites were conserved and under purifying selection (dN <dS). However, the most likely evolutionary codon models allowed for variable rates of selection across codon sites at both loci and, at the DQA, supported the hypothesis of positive selection acting on specific sites.

CONCLUSIONS

Observations of elevated genetic diversity and trans-species polymorphisms supported the conclusion that balancing selection may be acting on these loci. Furthermore, at the DQA, positive selection was occurring at antigen binding sites, suggesting that a few selected residues may play a significant role in equid immune function. Future studies in natural equid populations will be valuable for understanding the functional significance of the uniquely diverse DRA locus and for elucidating the mechanism maintaining diversity at these MHC loci.

Major histocompatibility complex diversity in the endangered Ethiopian wolf (Canis simensis)

The major histocompatibility complex (MHC) influences immune response to infection and vaccination. In most species, MHC genes are highly polymorphic, but few wild canid populations have been investigated. In Ethiopian wolves, we identified four DLA (dog leucocyte antigen)-DRB1, two DLA-DQA1 and five DQB1 alleles. Ethiopian wolves, the world's rarest canids with fewer than 500 animals worldwide, are further endangered and threatened by rabies. Major rabies outbreaks in the Bale Mountains of southern Ethiopia (where over half of the Ethiopian wolf population is located) have killed over 75% of wolves in the affected sub-populations. In 2004, following a rabies outbreak, 77 wolves were vaccinated, and 19 were subsequently recaptured to monitor the effectiveness of the intervention. Pre- and post-vaccination rabies antibody titres were available for 18 animals, and all of the animals sero-converted after vaccination. We compared the haplotype frequencies of this group of 18 with the post-vaccination antibody titre, and showed that one haplotype was associated with a lower response (uncorrected P < 0.03). In general, Ethiopian wolves probably have an adequate amount of MHC variation to ensure the survival of the species. However, we sampled only the largest Ethiopian wolf population in Bale, and did not take the smaller populations further north into consideration.

Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae

The allelic diversity of the DRB locus in major histocompatibility complex (MHC) genes was analyzed in the brown bear (Ursus arctos) from the Hokkaido Island of Japan, Siberia, and Kodiak of Alaska. Nineteen alleles of the DRB exon 2 were identified from a total of 38 individuals of U. arctos and were highly polymorphic. Comparisons of non-synonymous and synonymous substitutions in the antigen-binding sites of deduced amino acid sequences indicated evidence for balancing selection on the bear DRB locus. The phylogenetic analysis of the DRB alleles among three genera (Ursus, Tremarctos, and Ailuropoda) in the family Ursidae revealed that DRB allelic lineages were not separated according to species. This strongly shows trans-species persistence of DRB alleles within the Ursidae.

Increased genetic risk or protection for canine autoimmune lymphocytic thyroiditis in Giant Schnauzers depends on DLA class II genotype

Dogs represent an excellent comparative model for autoimmune thyroiditis as several dog breeds develop canine lymphocytic thyroiditis (CLT), which is clinically similar to Hashimoto's thyroiditis in human. We obtained evidence that dog leukocyte antigen (DLA) class II genotype function as either genetic risk factor that predisposes for CLT or as protective factor against the disease. Genetic diversity at their DLA-DRB1, -DQA1, and -DQB1 loci were defined and potential association to major histocompatibility complex II haplotypes and alleles was analyzed. Giant Schnauzers carrying the DLA-DRB1*01201/DQA1*00101/DQB1*00201 haplotype showed an increased risk (odds ratio of 6.5) for developing CLT. The same risk haplotype has, to date, been observed in three different breeds affected by this disease, Giant Schnauzer, Dobermann, and Labrador Retriever, indicating that it is a common genetic risk factor in a variety of breeds affected by this disease. Importantly, protection for development of the disease was found in dogs carrying the DLA-DRB1*01301/DQA1*00301/DQB1*00501 haplotype (odds ratio of 0.3).

Major histocompatibility complex variation at class II DQA locus in the brown hare (Lepus europaeus)

The major histocompatability complex (MHC) is a multigene family of receptors that bind and present antigenic peptides to T-cells. Genes of the MHC are characterized by an outstanding genetic polymorphism, which is considered to be maintained by positive selection. Sites involved in peptide binding form binding pockets (P) that are collectively termed the peptide-binding region (PBR). In this study, we examined the level of MHC genetic diversity within and among natural populations of brown hare (Lepus europaeus) from Europe and Anatolia choosing for analysis of the second exon of the DQA locus, one of the most polymorphic class II loci. We aimed at an integrated population genetic analysis of L. europeaus by (i) correlating MHC polymorphism to genetic variability and phylogenetic status estimated previously from maternally (mtDNA) and biparentally (allozymes, microsatellites) inherited loci; and (ii) comparing full-length exon amino acid polymorphism with functional polymorphism in the PBR and the binding pockets P1, P6 and P9. A substantial level of DQA exon 2 polymorphism was detected with two completely different set of alleles between the Anatolian and European populations. However, the phylogeny of full-length exon 2 Leeu-DQA alleles did not show a strong phylogeographic signal. The presence of balancing selection was supported by a statistically significant excess of nonsynonymous substitutions over synonymous in the PBR and a trans-species pattern of evolution detected after phylogenetic reconstruction. The differentiating patterns detected between genetic and functional polymorphism, i.e. the number and the distribution of pocket variants within and among populations, indicated a hierarchical action of selection pressures.

Polymorphism and selection in the major histocompatibility complex DRA and DQA genes in the family Equidae

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.

Selection, diversity and evolutionary patterns of the MHC class II DAB in free-ranging Neotropical marsupials

BACKGROUND

Research on the genetic architecture and diversity of the MHC has focused mainly on eutherian mammals, birds and fish. So far, studies on model marsupials used in laboratory investigations indicated very little or even no variation in MHC class II genes. However, natural levels of diversity and selection are unknown in marsupials as studies on wild populations are virtually absent. We used two endemic South American mouse opossums, Gracilinanus microtarsus and Marmosops incanus, to investigate characteristic features of MHC selection. This study is the first investigation of MHC selection in free-ranging Neotropical marsupials. In addition, the evolutionary history of MHC lineages within the group of marsupials was examined.

RESULTS

G. microtarsus showed extensive levels of MHC diversity within and among individuals as 47 MHC-DAB alleles and high levels of sequence divergence were detected at a minimum of four loci. Positively selected codon sites were identified, of which most were congruent with human antigen binding sites. The diversity in M. incanus was rather low with only eight observed alleles at presumably two loci. However, these alleles also revealed high sequence divergence. Again, positive selection was identified on specific codon sites, all congruent with human ABS and with positively selected sites observed in G. microtarsus. In a phylogenetic comparison alleles of M. incanus interspersed widely within alleles of G. microtarsus with four alleles being present in both species.

CONCLUSION

Our investigations revealed extensive MHC class II polymorphism in a natural marsupial population, contrary to previous assumptions. Furthermore, our study confirms for the first time in marsupials the presence of three characteristic features common at MHC loci of eutherian mammals, birds and fish: large allelic sequence divergence, positive selection on specific sites and trans-specific polymorphism.

DLA-DRB1, DQA1, and DQB1 alleles and haplotypes in North American Gray Wolves

The canine major histocompatibility complex contains highly polymorphic genes, many of which are critical in regulating immune response. Since domestic dogs evolved from Gray Wolves (Canis lupus), common DLA class II alleles should exist. Sequencing was used to characterize 175 Gray Wolves for DLA class II alleles, and data from 1856 dogs, covering 85 different breeds of mostly European origin, were available for comparison. Within wolves, 28 new alleles were identified, all occurring in at least 2 individuals. Three DLA-DRB1, 8 DLA-DQA1, and 6 DLA-DQB1 alleles also identified in dogs were present. Twenty-eight haplotypes were identified, of which 2 three-locus haplotypes, and many DLA-DQA1/DQB1 haplotypes, are also found in dogs. The wolves studied had relatively few dog DLA alleles and may therefore represent a remnant population descended from Asian wolves. The single European wolf included carried a haplotype found in both these North American wolves and in many dog breeds. Furthermore, one wolf DQB1 allele has been found in Shih Tzu, a breed of Asian origin. These data suggest that the wolf ancestors of Asian and European dogs may have had different gene pools, currently reflected in the DLA alleles present in dog breeds.

Dog leucocyte antigen class II diversity and relationships among indigenous dogs of the island nations of Indonesia (Bali), Australia and New Guinea

The genetic polymorphism at the dog leucocyte antigen (DLA) class II loci DQA1, DQB1 and DRB1 was studied in a large genetically diverse population of feral and wild-type dogs from the large island nations of Indonesia (Bali), Australia and New Guinea (Bali street dog, dingo and New Guinea singing dog, respectively). Sequence-based typing (SBT) of the hypervariable region of DLA-DRB1, -DQA1 and -DQB1 alleles was used to determine genetic diversity. No new DQA1 alleles were recognized among the three dog populations, but five novel DLA-DRB1 and 2 novel DLA-DQB1 allele sequences were detected. Additional unknown alleles were postulated to exist in Bali street dogs, as indicated by the large percentage of individuals (15%-33%) that had indeterminate DRB1, DQA1 and DQB1 alleles by SBT. All three groups of dogs possessed alleles that were relatively uncommon in conventional purebreds. The New Guinea singing dog and dingo shared alleles that were not present in the Bali street dogs. These findings suggested that the dingo was more closely related to indigenous dogs from New Guinea. Feral dog populations, in particular large ones such as that of Bali, show genetic diversity that existed prior to phenotypic selection for breeds originating from their respective regions. This diversity needs to be identified and maintained in the face of progressive Westernization. These populations deserve further study as potential model populations for the evolution of major histocompatibility complex alleles, for the study of canine genetic diversity, for the development of dog breeds and for studies on the comigration of ancestral human and dog populations.

Clonal origin and evolution of a transmissible cancer

The transmissible agent causing canine transmissible venereal tumor (CTVT) is thought to be the tumor cell itself. To test this hypothesis, we analyzed genetic markers including major histocompatibility (MHC) genes, microsatellites, and mitochondrial DNA (mtDNA) in naturally occurring tumors and matched blood samples. In each case, the tumor is genetically distinct from its host. Moreover, tumors collected from 40 dogs in 5 continents are derived from a single neoplastic clone that has diverged into two subclades. Phylogenetic analyses indicate that CTVT most likely originated from a wolf or an East Asian breed of dog between 200 and 2500 years ago. Although CTVT is highly aneuploid, it has a remarkably stable genotype. During progressive growth, CTVT downmodulates MHC antigen expression. Our findings have implications for understanding genome instability in cancer, natural transplantation of allografts, and the capacity of a somatic cell to evolve into a transmissible parasite.

Diversity and duplication of DQB and DRB-like genes of the MHC in baleen whales (suborder: Mysticeti)

The molecular diversity and phylogenetic relationships of two class II genes of the baleen whale major histocompatibility complex were investigated and compared to toothed whales and out-groups. Amplification of the DQB exon 2 provided sequences showing high within-species and between-species nucleotide diversity and uninterrupted reading frames consistent with functional class II loci found in related mammals (e.g., ruminants). Cloning of amplified products indicated gene duplication in the humpback whale and triplication in the southern right whale, with average nucleotide diversity of 5.9 and 6.3%, respectively, for alleles of each species. Significantly higher nonsynonymous divergence at sites coding for peptide binding (32% for humpback and 40% for southern right) suggested that these loci were subject to positive (overdominant) selection. A population survey of humpback whales detected 23 alleles, differing by up to 21% of their inferred amino acid sequences. Amplification of the DRB exon 2 resulted in two groups of sequences. One was most similar to the DRB3 of the cow and present in all whales screened to date, including toothed whales. The second was most similar to the DRB2 of the cow and was found only in the bowhead and right whales. Both loci showed low diversity among species and apparent loss of function or altered function including interruption of reading frames. Finally, comparison of inferred protein sequence of the DRB3-like locus suggested convergence with the DQB, perhaps resulting from intergenic conversion or recombination.

Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Major histocompatibility complex (MHC) genes play important role in host-parasite interactions and parasites are crucial factors influencing the population dynamics of hosts. We described the structure and diversity of exon 2 of the MHC class II DQA gene in three species of voles (Arvicolinae) exhibiting regular multi-annual fluctuations of population density and analysed the processes leading to the observed MHC polymorphism. By using cloning-sequencing methodology and capillary electrophoresis-single strand conformation polymorphism, we described seven sequences in the water, eight in the common, and seven in the bank voles coming from an area of 70 km(2) around the Nozeroy canton in the Jura Mountains (Franche Comté, France). All exon 2 sequences translate to give unique amino acid sequences and positive selection was found to act very intensively on antigen binding sites. We documented the presence of recombination at vole DQA region but its importance in generating allelic polymorphism seems to be relatively limited. For the first time within rodents, we documented the duplication of the DQA gene in all three species with both copies being transcriptionally active. Phylogenetic analysis of allelic sequences revealed extensive trans-species polymorphism within the subfamily although no alleles were shared between species in our data set. We discuss possible role of parasites in forming the recent polymorphism pattern of the DQA locus in voles.

MHC promoter polymorphism in grey wolves and domestic dogs

A functional immune system requires a tight control over major histocompatibility complex (MHC) gene transcription, as the abnormal MHC expression patterns of severe immunodeficiency and autoimmune diseases demonstrate. Although the regulation of MHC expression has been well documented in humans and mice, little is known in other species. In this study, we detail the level of polymorphism in wolf and dog MHC gene promoters. The promoter regions of the DRB, DQA and DQB locus were sequenced in 90 wolves and 90 dogs. The level of polymorphism was high in the DQB promoters, with variation found within functionally relevant regions, including binding sites for transcription factors. Clear associations between DQB promoters and exon 2 alleles were noted in wolves, indicating strong linkage disequilibrium in this region. Low levels of polymorphism were found within the DRB and DQA promoter regions. However, a variable site was identified within the T box, a TNF-alpha response element, of the DQA promoter. Furthermore, we identified a previously unrecognised 18-base-pair deletion within exon 1 of the DQB locus.

Patterns of variability at the major histocompatibility class II alpha locus in Atlantic salmon contrast with those at the class I locus

In order to investigate the mechanisms creating and maintaining variability at the major histocompatibility (MH) class II alpha (DAA) locus we examined patterns of polymorphism in two isolated Atlantic salmon populations which share a common post-glacial origin. As expected from their common origin, but contrary to the observation at the MH class I locus, these populations shared the majority of DAA alleles: out of 17 sequences observed, 11 were common to both populations. Recombination seems to play a more important role in the origin of new alleles at the class II alpha locus than at the class I locus. A greater than expected proportion of sites inferred to be positively selected (potentially peptide binding residues, PBRs) were found to be involved in recombination events, suggesting a mechanism for increasing MH variability through an interaction between recombination and natural selection. Thus it appears that although selection and recombination are important mechanisms for the evolution of both class II alpha and class I loci in the Atlantic salmon, the pattern of variability differs markedly between these classes of MH loci.

MHC class II DRB variability and parasite load in the striped mouse (Rhabdomys pumilio) in the Southern Kalahari

Major histocompatibility complex (MHC) variability is believed to be maintained by pathogen-driven selection, mediated either through heterozygous advantage or frequency-dependent selection. However, empirical support for these hypotheses under natural conditions is rare. In this study, we investigated the genetic constitution of the functionally important MHC class II gene (DRB exon 2) and the parasite load in a population of the striped mouse (Rhabdomys pumilio) in the Southern Kalahari. Fifty-eight individuals were genetically examined and the endoparasite load was quantified by counting fecal helminth eggs by using a modified McMaster technique. Thirty-four animals (58.6%) were infected. We identified 20 different MHC alleles with high levels of sequence divergence between alleles. Particularly, the antigen-binding sites revealed a significant higher rate of nonsynonymous substitutions (d(N)) than synonymous substitutions (d(S)), giving strong evidence of balancing selection. Heterozygosity did influence the infection status (being infected or not) and the individual fecal egg count (FEC) value with significantly higher values observed in homozygous individuals. Furthermore, a positive relationship was found between specific alleles and parasite load. The allele Rhpu-DRB*1 significantly occurred more frequently in infected individuals and in individuals with high FEC values (high parasite load). Individuals with the allele Rhpu-DRB*1 had a 1.5-fold higher chance of being infected than individuals without this allele (odds ratio test, P < 0.05). Contrarily, the allele Rhpu-DRB*8 significantly occurred more frequent in individuals with low FEC values. Our results support the hypotheses that MHC polymorphism in R. pumilio is maintained through pathogen-driven selection acting by both heterozygosity advantage and frequency-dependent selection.