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

The genome of the grey wolf, one of the most widely distributed land mammal species, has been subjected to both stochastic factors, including biogeographical subdivision and population fragmentation, and strong selection during the domestication of the dog. To explore the effects of drift and selection on the partitioning of MHC variation in the diversification of species, we present nine DQA, 10 DQB, and 17 DRB1 sequences of the second exon for European wolves and compare them with sequences of North American wolves and dogs. The relatively large number of class II alleles present in both European and North American wolves attests to their large historical population sizes, yet there are few alleles shared between these regions at DQB and DRB1. Similarly, the dog has an extensive array of class II MHC alleles, a consequence of a genetically diverse origin, but allelic overlap with wolves only at DQA. Although we might expect a progression from shared alleles to shared allelic lineages during differentiation, the partitioning of diversity between wolves and dogs at DQB and DRB1 differs from that at DQA. Furthermore, an extensive region of nucleotide sequence shared between DRB1 and DQB alleles and a shared motif suggests intergenic recombination may have contributed to MHC diversity in the Canidae.     

A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution

The major histocompatibility complex (MHC) has an integral role in the immune system, and hence diversity at its genes may be of particular importance for the health of populations. In large populations, balancing selection maintains diversity in MHC genes, but theoretical expectations indicate that this form of selection is absent or inefficient in small populations. We examine the level of diversity at three MHC class II loci in the wolf population of Scandinavia, a population naturally recolonized with a genetic contribution from as few as three founders, and in four neighbouring wolf populations. In the Scandinavian wolf population, two alleles were found for each locus and the distribution of alleles is compatible with their linkage into two haplotypes. Changes in the level of heterozygosity over time since recolonization demonstrate the effects of the proposed arrival of an immigrant wolf. The maintenance of diversity is shown to be compatible with a neutral, random allocation of alleles, in conjunction with crossing between packs. A total of 15 DRB1, seven DQA and 10 DQB1 alleles are found in four neighbouring wolf populations, with substantial sharing across populations. Even in these larger populations, bottlenecks and fragmentation with consequent genetic drift are likely to have resulted in few indicators for balancing selection and significant differentiation of populations.     

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.     

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.     

Genes of domestic mammals augmented by backcrossing with wild ancestors

Both archaeological data and the presence of few mitochondrial DNA lineages suggest that most widespread domestic mammals (cattle, sheep, goats, pigs and dogs) derive from only a handful of domestication events. However, each of these species shows a high level of diversity at the nuclear genes of the major histocompatibility complex (MHC). Through simulations incorporating various degrees of population subdivision, growth rate and selection, we demonstrate that the numerous MHC DRB alleles that are present in modern domestic mammals implies that substantial backcrossing with wild ancestors, either accidental or intentional, has been important in shaping the genetic diversity of our domesticates. These results support the view that, contrary to common assumption, domestic and wild lineages might not have been clearly separated throughout their history.     

Genome Sequencing Highlights the Dynamic Early History of Dogs

To identify genetic changes underlying dog domestication and reconstruct their early evolutionary history, we generated high-quality genome sequences from three gray wolves, one from each of the three putative centers of dog domestication, two basal dog lineages (Basenji and Dingo) and a golden jackal as an outgroup. Analysis of these sequences supports a demographic model in which dogs and wolves diverged through a dynamic process involving population bottlenecks in both lineages and post-divergence gene flow. In dogs, the domestication bottleneck involved at least a 16-fold reduction in population size, a much more severe bottleneck than estimated previously. A sharp bottleneck in wolves occurred soon after their divergence from dogs, implying that the pool of diversity from which dogs arose was substantially larger than represented by modern wolf populations. We narrow the plausible range for the date of initial dog domestication to an interval spanning 11–16 thousand years ago, predating the rise of agriculture. In light of this finding, we expand upon previous work regarding the increase in copy number of the amylase gene (AMY2B) in dogs, which is believed to have aided digestion of starch in agricultural refuse. We find standing variation for amylase copy number variation in wolves and little or no copy number increase in the Dingo and Husky lineages. In conjunction with the estimated timing of dog origins, these results provide additional support to archaeological finds, suggesting the earliest dogs arose alongside hunter-gathers rather than agriculturists. Regarding the geographic origin of dogs, we find that, surprisingly, none of the extant wolf lineages from putative domestication centers is more closely related to dogs, and, instead, the sampled wolves form a sister monophyletic clade. This result, in combination with dog-wolf admixture during the process of domestication, suggests that a re-evaluation of past hypotheses regarding dog origins is necessary.     

Diversification of porcine MHC class II genes: evidence for selective advantage

The major histocompatibility complex (MHC) is an immunological gene-dense region of high diversity in mammalian species. Sus scrofa was domesticated by at least six independent events over Eurasia during the Holocene period. It has been hypothesized that the level and distribution of MHC variation in pig populations reflect genetic selection and environmental influences. In an effort to define the complexity of MHC polymorphisms and the role of selection in the generation of class II gene diversity (DQB, DRB1, and pseudogene ΨDRB3), DNA from globally distributed unrelated domestic pigs of European and Asian origins and a Suidae out-group was analyzed. The number of pseudogene alleles identified (ΨDRB3 33) was greater than those found in the expressed genes (DQB 20 and DRB1 23) but the level of observed heterozygosity (ΨDRB3 0.452, DQB 0.732, and DRB1 0.767) and sequence diversity (ΨDRB3 0.029, DQB 0.062, and DRB1 0.074) were significantly lower in the pseudogene, respectively. The substitution ratios reflected an excess of d _N (DQB 1.476, DRB1 1.724, and ΨDRB3 0.508) and the persistence of expressed gene alleles suggesting the influence of balancing selection, while the pseudogene was undergoing purifying selection. The lack of a clear MHC phylogeographic tree, coupled with close genetic distances observed between the European and Asian populations (DQB 0.047 and DRB1 0.063) suggested that unlike observations using mtDNA, the MHC diversity lacks phylogeographic structure and appears to be globally uniform. Taken together, these results suggest that, despite regional differences in selective breeding and environments, no skewing of MHC diversity has occurred. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species'' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species'' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.     

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber)

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.     

On the origin of mongrels: evolutionary history of free-breeding dogs in Eurasia

Although a large part of the global domestic dog population is free-ranging and free-breeding, knowledge of genetic diversity in these free-breeding dogs (FBDs) and their ancestry relations to pure-breed dogs is limited, and the indigenous status of FBDs in Asia is still uncertain. We analyse genome-wide SNP variability of FBDs across Eurasia, and show that they display weak genetic structure and are genetically distinct from pure-breed dogs rather than constituting an admixture of breeds. Our results suggest that modern European breeds originated locally from European FBDs. East Asian and Arctic breeds show closest affinity to East Asian FBDs, and they both represent the earliest branching lineages in the phylogeny of extant Eurasian dogs. Our biogeographic reconstruction of ancestral distributions indicates a gradual westward expansion of East Asian indigenous dogs to the Middle East and Europe through Central and West Asia, providing evidence for a major expansion that shaped the patterns of genetic differentiation in modern dogs. This expansion was probably secondary and could have led to the replacement of earlier resident populations in Western Eurasia. This could explain why earlier studies based on modern DNA suggest East Asia as the region of dog origin, while ancient DNA and archaeological data point to Western Eurasia.     

Balancing selection and heterozygote advantage in major histocompatibility complex loci of the bottlenecked Finnish wolf population

Maintaining effective immune response is an essential factor in the survival of small populations. One of the most important immune gene regions is the highly polymorphic major histocompatibility complex (MHC). We investigated how a population bottleneck and recovery have influenced the diversity and selection in three MHC class II loci, DLA‐DRB1, DLA‐DQA1 and DLA‐DQB1, in the Finnish wolf population. We studied the larger Russian Karelian wolf population for comparison and used 17 microsatellite markers as reference loci. The Finnish and Karelian wolf populations did not differ substantially in their MHC diversities ( = 0.047, P = 0.377), but differed in neutral microsatellite diversities ( = 0.148, P = 0.008). MHC allele frequency distributions in the Finnish population were more even than expected under neutrality, implying balancing selection. In addition, an excess of nonsynonymous compared to synonymous polymorphisms indicated historical balancing selection. We also studied association between helminth (Trichinella spp. and Echinococcus canadensis) prevalence and MHC diversity at allele and SNP level. MHC‐heterozygous wolves were less often infected by Trichinella spp. and carriers of specific MHC alleles, SNP haplotypes and SNP alleles had less helminth infections. The associated SNP haplotypes and alleles were shared by different MHC alleles, which emphasizes the necessity of single‐nucleotide‐level association studies also in MHC. Here, we show that strong balancing selection has had similar effect on MHC diversities in the Finnish and Russian Karelian wolf populations despite significant genetic differentiation at neutral markers and small population size in the Finnish population.     

Genetic wealth,population health: Major histocompatibility complex variation in captive and wild ring‐tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring‐tailed lemur (Lemur catta) as a model, we compared two populations under long‐term study: a relatively “open,” wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003–2013) and a “closed,” captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980–2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC‐DRB diversity and, across populations, we compared the number of unique MHC‐DRB alleles and their distributions. Wild individuals possessed more MHC‐DRB alleles than did captive individuals, and overall, the wild population had more unique MHC‐DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive‐breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.     

Low major histocompatibility complex diversity in the Tasmanian devil predates European settlement and may explain susceptibility to disease epidemics

The Tasmanian devil (Sarcophilus harrisii) is at risk of extinction owing to the emergence of a contagious cancer known as devil facial tumour disease (DFTD). The emergence and spread of DFTD has been linked to low genetic diversity in the major histocompatibility complex (MHC). We examined MHC diversity in historical and ancient devils to determine whether loss of diversity is recent or predates European settlement in Australia. Our results reveal no additional diversity in historical Tasmanian samples. Mainland devils had common modern variants plus six new variants that are highly similar to existing alleles. We conclude that low MHC diversity has been a feature of devil populations since at least the Mid-Holocene and could explain their tumultuous history of population crashes.     

MHC-DRB3 variation in a free-living population of the European bison, Bison bonasus

MHC genes play a crucial role in pathogen recognition and are the most polymorphic genes in vertebrates. Loss of variation in these genes in bottlenecked species is thought to put their survival at risk. We examined variation at the MHC II DRB3 locus in the European bison, Bison bonasus, a species that has undergone an extreme bottleneck: the current population originated from only 12 founders. We also tested for the association of DRB3 genes with the incidence of posthitis, a disease affecting the reproductive organs of bulls and posing a new threat to the survival of the species. We found very limited MHC diversity, with only four alleles segregating in a sample of 172 individuals from a free‐ranging Białowieża population. The alleles were highly divergent and revealed the hallmark of positive selection acting on them in the past, that is, a significant excess of nonsynonymous substitutions. This excess was concentrated in putative antigen‐binding sites, suggesting that selection was driven by pathogens. However, we did not observe departures from Hardy–Weinberg equilibrium, an indicator of strong ongoing selection. Neither have we found a significant association between DRB3 alleles or genotypes and susceptibility to posthitis. Alleles conferring resistance to males may have been lost during the extreme bottleneck the species had undergone.     

Genetic assessment of the Iberian wolf Canis lupus signatus captive breeding program

The main goal of ex situ conservation programs is to improve the chances of long term survival of natural populations by founding and managing captive colonies that can serve as a source of individuals for future reintroductions or to reinforce existing populations. The degree in which a captive breeding program has captured the genetic diversity existing in the source wild population has seldom been evaluated. In this study we evaluate the genetic diversity in wild and captive populations of the Iberian wolf, Canis lupus signatus, in order to assess how much genetic diversity is being preserved in the ongoing ex situ conservation program for this subspecies. A sample of domestic dogs was also included in the analysis for comparison. Seventy-four wolves and 135 dogs were genotyped at 13 unlinked microsatellite loci. The results show that genetic diversity in Iberian wolves is comparable in magnitude to that of other wild populations of gray wolf. Both the wild and the captive Iberian wolf populations have a similarly high genetic diversity indicating that no substantial loss of diversity has occurred in the captive-breeding program. The effective number of founders of the program was estimated as ∼ ∼16, suggesting that all founders in the studbook pedigree were genetically independent. Our results emphasize also the genetic divergence between wolves and domestic dogs and indicate that our set of 13 microsatellite loci provide a powerful diagnostic test to distinguish wolves, dogs and their hybrids.     

Major histocompatibility complex alleles associated with parasite susceptibility in wild giant pandas

Major histocompatibility complex (MHC) polymorphism is thought to be driven by antagonistic coevolution between pathogens and hosts, mediated through either overdominance or frequency-dependent selection. However, investigations under natural conditions are still rare for endangered mammals which often exhibit depleted variation, and the mechanism of selection underlying the maintenance of characteristics remains a considerable debate. In this study, 87 wild giant pandas were used to investigate MHC variation associated with parasite load. With the knowledge of the MHC profile provided by the genomic data of the giant panda, seven DRB1, seven DQA1 and eight DQA2 alleles were identified at each single locus. Positive selection evidenced by a significantly higher number of non-synonymous substitutions per non-synonymous codon site relative to synonymous substitutions per synonymous codon site could only be detected at the DRB1 locus, which leads to the speculation that DRB1 may have a more important role in dealing with parasite infection for pandas. Coprological analyses revealed that 55.17% of individuals exhibited infection with 1–2 helminthes and 95.3% of infected pandas carried Baylisascaris shroederi. Using a generalized linear model, we found that Aime-DRB1*10 was significantly associated with parasite infection, but no resistant alleles could be detected. MHC heterozygosity of the pandas was found to be uncorrelated with the infection status or the infection intensity. These results suggested that the possible selection mechanisms in extant wild pandas may be frequency dependent rather than being determined by overdominance selection. Our findings could guide the candidate selection for the ongoing reintroduction or translocation of pandas.     

Substantial functional diversity accompanies limited major histocompatibility complex class II variability in golden jackal (Canis aureus): A comparison between two wild Canis species in Croatia

The golden jackal (Canis aureus) is one of the less studied carnivores and research on its major histocompatibility complex (MHC) variability is just at its early stages. MHC genes encode cell-surface receptors that serve to bind and present antigens to T cells, which is essential to initiating specific immunological responses in vertebrates. In this paper we present for the first time patterns of genetic diversity and natural selection on MHC class II DLA-DRB1, DQA1 and DQB1 loci in the golden jackal using samples from two geographically distinct regions in Croatia and further compare them to the values found in its congener grey wolf (Canis lupus). Diversity of golden jackals at all three loci was markedly lower than that of grey wolves (allelic richness values were 4, 2 and 3 in jackal versus 11.9, 6.6 and 10.2 in wolves for DRB1, DQA1 and DQB1, respectively) and can be attributed to a genetic drift rather than to the lack of historical positive selection. The finding of high evolutionary distances (16.3% for DRB1 and 8.5% for DQB1) and a substantial number of codons predicted to be under the influence of positive selection (11 for DRB1 and 9 for DQB1) suggests that the investigated golden jackal population still contains considerable functional diversity necessary for the presentation of varied foreign peptides. In contrast to neutral genetic variation, our results suggest that the Dalmatian population has a higher MHC diversity than the Slavonian population, casting doubt on its supposed isolation and calling for a more extensive investigation of the MHC variability of southern Balkan jackal populations.     

Testing the myth: tolerant dogs and aggressive wolves

Cooperation is thought to be highly dependent on tolerance. For example, it has been suggested that dog–human cooperation has been enabled by selecting dogs for increased tolerance and reduced aggression during the course of domestication (‘emotional reactivity hypothesis’). However, based on observations of social interactions among members of captive packs, a few dog–wolf comparisons found contradictory results. In this study, we compared intraspecies aggression and tolerance of dogs and wolves raised and kept under identical conditions by investigating their agonistic behaviours and cofeeding during pair-wise food competition tests, a situation that has been directly linked to cooperation. We found that in wolves, dominant and subordinate members of the dyads monopolized the food and showed agonistic behaviours to a similar extent, whereas in dogs these behaviours were privileges of the high-ranking individuals. The fact that subordinate dogs rarely challenged their higher-ranking partners suggests a steeper dominance hierarchy in dogs than in wolves. Finally, wolves as well as dogs showed only rare and weak aggression towards each other. Therefore, we suggest that wolves are sufficiently tolerant to enable wolf–wolf cooperation, which in turn might have been the basis for the evolution of dog–human cooperation (canine cooperation hypothesis).