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

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.

Sarcoptic mange severity is associated with reduced genomic variation and evidence of selection in Yellowstone National Park wolves (Canis lupus)

Population genetic theory posits that molecular variation buffers against disease risk. Although this "monoculture effect" is well supported in agricultural settings, its applicability to wildlife populations remains in question. In the present study, we examined the genomics underlying individual-level disease severity and population-level consequences of sarcoptic mange infection in a wild population of canids. Using gray wolves (Canis lupus) reintroduced to Yellowstone National Park (YNP) as our focal system, we leveraged 25 years of observational data and biobanked blood and tissue to genotype 76,859 loci in over 400 wolves. At the individual level, we reported an inverse relationship between host genomic variation and infection severity. We additionally identified 410 loci significantly associated with mange severity, with annotations related to inflammation, immunity, and skin barrier integrity and disorders. We contextualized results within environmental, demographic, and behavioral variables, and confirmed that genetic variation was predictive of infection severity. At the population level, we reported decreased genome-wide variation since the initial gray wolf reintroduction event and identified evidence of selection acting against alleles associated with mange infection severity. We concluded that genomic variation plays an important role in disease severity in YNP wolves. This role scales from individual to population levels, and includes patterns of genome-wide variation in support of the monoculture effect and specific loci associated with the complex mange phenotype. Results yielded system-specific insights, while also highlighting the relevance of genomic analyses to wildlife disease ecology, evolution, and conservation.

Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence

The prevalence and intensity of parasites in wild hosts varies across space and is a key determinant of infection risk in humans, domestic animals and threatened wildlife. Because the immune system serves as the primary barrier to infection, replication and transmission following exposure, we here consider the environmental drivers of immunity. Spatial variation in parasite pressure, abiotic and biotic conditions, and anthropogenic factors can all shape immunity across spatial scales. Identifying the most important spatial drivers of immunity could help pre-empt infectious disease risks, especially in the context of how large-scale factors such as urbanization affect defence by changing environmental conditions. We provide a synthesis of how to apply macroecological approaches to the study of ecoimmunology (i.e. macroimmunology). We first review spatial factors that could generate spatial variation in defence, highlighting the need for large-scale studies that can differentiate competing environmental predictors of immunity and detailing contexts where this approach might be favoured over small-scale experimental studies. We next conduct a systematic review of the literature to assess the frequency of spatial studies and to classify them according to taxa, immune measures, spatial replication and extent, and statistical methods. We review 210 ecoimmunology studies sampling multiple host populations. We show that whereas spatial approaches are relatively common, spatial replication is generally low and unlikely to provide sufficient environmental variation or power to differentiate competing spatial hypotheses. We also highlight statistical biases in macroimmunology, in that few studies characterize and account for spatial dependence statistically, potentially affecting inferences for the relationships between environmental conditions and immune defence. We use these findings to describe tools from geostatistics and spatial modelling that can improve inference about the associations between environmental and immunological variation. In particular, we emphasize exploratory tools that can guide spatial sampling and highlight the need for greater use of mixed-effects models that account for spatial variability while also allowing researchers to account for both individual- and habitat-level covariates. We finally discuss future research priorities for macroimmunology, including focusing on latitudinal gradients, range expansions and urbanization as being especially amenable to large-scale spatial approaches. Methodologically, we highlight critical opportunities posed by assessing spatial variation in host tolerance, using metagenomics to quantify spatial variation in parasite pressure, coupling large-scale field studies with small-scale field experiments and longitudinal approaches, and applying statistical tools from macroecology and meta-analysis to identify generalizable spatial patterns. Such work will facilitate scaling ecoimmunology from individual- to habitat-level insights about the drivers of immune defence and help predict where environmental change may most alter infectious disease risk.

Sequence diversity of major histo-compatibility complex class II DQA1 in Indian Tharparkar cattle: novel alleles and in-silico analysis

Exon 2 of MHC class II gene codes for the first domain of the molecule that forms the peptide-binding groove and its polymorphism partly explains functional MHC diversity. A 850 bp DQA1 gene fragment spanning from intron I to exon III was typed by sequencing of 40 Tharparkar cattle of various agro-climatic zones of northern India along with 10 Tharparkar crossbreds. On analysis of nucleotide sequences, a total of 30 polymorphic sites (1 insertion and 29 SNPs) were identified in 14 MHC alleles leading to amino acid changes in 5 places in 249 bp (exon 2). Five new BoLa DQA1 alleles were identified and reported. The within group mean distance was highest in Tharparkar herd of Bikaner (0.045) and lowest (0.020) in that of Surathgarh (breeding tract) whereas, between groups mean distance was highest in Bikaner Tharparkar-Suratgarh Tharparkar pair. There was excess of nonsynonymous over synonymous nucleotide substitutions in the present study. The effects of these substitutions were predicted using I-Mutant and Panther online resources. The mean ratio of dN/dS was found to be >1.0 at 12 codons with two mutation hotspots at 13th codon (P = 0.002) and 64th codon (P = 0.01). The phylo-geographic analysis revealed that alleles 5, 7 and 13 formed a different cluster with alleles 7 and 13 grouped by the most frequent allele (BoLa-DQA*1401).

Characterization of the peripheral blood transcriptome and adaptive evolution of the MHC I and TLR gene families in the wolf (Canis lupus)

Background

The wolf (Canis lupus) is one of the most widely distributed terrestrial mammals, because it is well adapted to various ecological niches and their corresponding pathogen environments. Immunological competence is a crucial factor involved in adapting to a changing environment and fighting pathogen infection in animals. In this study, the peripheral blood transcriptome of wolves was generated via RNA-seq to advance understanding of the wolf immunome, with a special focus on the major histocompatibility complex class I (MHC I) and toll-like receptor (TLR) gene families, which are involved in pathogen recognition and defense.

Results

The blood transcriptomic libraries of eight wolves originating from Tibet and Inner Mongolia were sequenced, and approximately 383 million reads were generated. Using a genome-guided assembly strategy, we obtained 123,851 unigenes, with a mean length of 845 bp and an N50 length of 1121 bp. On the basis of BLAST searches against the NCBI non-redundant protein database (Nr), a total of 36,192 (29.22%) unigenes were annotated. For functional classification, 24,663 unigenes were assigned to 13,016 Gene Ontology (GO) terms belonging to 51 sub-categories of the three main GO categories. Additionally, 7682 unigenes were classified into 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, in which the most represented functional sub-categories were signal transduction and the immune system, and 16,238 unigenes were functionally classified into 25 Eukaryotic Orthologous Groups (KOG) categories. We observed an overall higher ω (d_N/d_S) value at antigen-binding sites (ABSs) than at non-ABS regions as well as clear evidence of intergenic/intragenic recombination events at wolf MHC I loci. Additionally, our analysis revealed that carnivorous TLRs were dominated by purifying selection, with mean ω values at each TLR locus ranging from 0.173 to 0.527. However, we also found significant instances of positive selection that acted on several codons in pathogen recognition domains and were linked to species-specific differences in pathogen recognition.

Conclusions

This study represents the first attempt to characterize the blood transcriptome of the wolf and to highlight the value of investigating the immune system. Balancing selection and recombination have contributed to the historical evolution of wolf MHC I genes. Moreover, TLRs in carnivores have undergone adaptive evolution against the background of purifying selection, and a high level of adaptive evolution was detected in the wolf TLR system.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3983-0) contains supplementary material, which is available to authorized users.

Choosy Wolves? Heterozygote Advantage But No Evidence of MHC-Based Disassortative Mating

A variety of nonrandom mate choice strategies, including disassortative mating, are used by vertebrate species to avoid inbreeding, maintain heterozygosity and increase fitness. Disassortative mating may be mediated by the major histocompatibility complex (MHC), an important gene cluster controlling immune responses to pathogens. We investigated the patterns of mate choice in 26 wild-living breeding pairs of gray wolf (Canis lupus) that were identified through noninvasive genetic methods and genotyped at 3 MHC class II and 12 autosomal microsatellite (STR) loci. We tested for deviations from random mating and evaluated the covariance of genetic variables at functional and STR markers with fitness proxies deduced from pedigree reconstructions. Results did not show evidences of MHC-based disassortative mating. Rather we found a higher peptide similarity between mates at MHC loci as compared with random expectations. Fitness values were positively correlated with heterozygosity of the breeders at both MHC and STR loci, whereas they decreased with relatedness at STRs. These findings may indicate fitness advantages for breeders that, while avoiding highly related mates, are more similar at the MHC and have high levels of heterozygosity overall. Such a pattern of MHC-assortative mating may reflect local coadaptation of the breeders, while a reduction in genetic diversity may be balanced by heterozygote advantages.

Multilocus detection of wolf x dog hybridization in italy, and guidelines for marker selection

Hybridization and introgression can impact the evolution of natural populations. Several wild canid species hybridize in nature, sometimes originating new taxa. However, hybridization with free-ranging dogs is threatening the genetic integrity of grey wolf populations (Canis lupus), or even the survival of endangered species (e.g., the Ethiopian wolf C. simensis). Efficient molecular tools to assess hybridization rates are essential in wolf conservation strategies. We evaluated the power of biparental and uniparental markers (39 autosomal and 4 Y-linked microsatellites, a melanistic deletion at the β-defensin CBD103 gene, the hypervariable domain of the mtDNA control-region) to identify the multilocus admixture patterns in wolf x dog hybrids. We used empirical data from 2 hybrid groups with different histories: 30 presumptive natural hybrids from Italy and 73 Czechoslovakian wolfdogs of known hybrid origin, as well as simulated data. We assessed the efficiency of various marker combinations and reference samples in admixture analyses using 69 dogs of different breeds and 99 wolves from Italy, Balkans and Carpathian Mountains. Results confirmed the occurrence of hybrids in Italy, some of them showing anomalous phenotypic traits and exogenous mtDNA or Y-chromosome introgression. Hybridization was mostly attributable to village dogs and not strictly patrilineal. The melanistic β-defensin deletion was found only in Italian dogs and in putative hybrids. The 24 most divergent microsatellites (largest wolf-dog FST values) were equally or more informative than the entire panel of 39 loci. A smaller panel of 12 microsatellites increased risks to identify false admixed individuals. The frequency of F1 and F2 was lower than backcrosses or introgressed individuals, suggesting hybridization already occurred some generations in the past, during early phases of wolf expansion from their historical core areas. Empirical and simulated data indicated the identification of the past generation backcrosses is always uncertain, and a larger number of ancestry-informative markers is needed.

Sexual selection explains more functional variation in the mammalian major histocompatibility complex than parasitism

Understanding drivers of genetic diversity at the major histocompatibility complex (MHC) is vitally important for predicting how vertebrate immune defence might respond to future selection pressures and for preserving immunogenetic diversity in declining populations. Parasite-mediated selection is believed to be the major selective force generating MHC polymorphism, and while MHC-based mating preferences also exist for multiple species including humans, the general importance of mate choice is debated. To investigate the contributions of parasitism and sexual selection in explaining among-species variation in MHC diversity, we applied comparative methods and meta-analysis across 112 mammal species, including carnivores, bats, primates, rodents and ungulates. We tested whether MHC diversity increased with parasite richness and relative testes size (as an indicator of the potential for mate choice), while controlling for phylogenetic autocorrelation, neutral mutation rate and confounding ecological variables. We found that MHC nucleotide diversity increased with parasite richness for bats and ungulates but decreased with parasite richness for carnivores. By contrast, nucleotide diversity increased with relative testes size for all taxa. This study provides support for both parasite-mediated and sexual selection in shaping functional MHC polymorphism across mammals, and importantly, suggests that sexual selection could have a more general role than previously thought.

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.

MHC class II association study in eight breeds of dog with hypoadrenocorticism

Canine hypoadrenocorticism is an endocrine disorder characterised by inadequate secretion of steroid hormones from the adrenal glands. Pathology results from immune-mediated destruction of the adrenal cortex, which is similar to that seen in the human Addison's disease. Both the canine and human diseases have similar clinical presentation, with the diagnosis based on performing a dynamic adrenocorticotropic hormone stimulation test. MHC class II has previously been associated with the human and canine diseases. In the current study, we conducted an MHC class II association study in eight breeds of dog with diagnoses of hypoadrenocorticism. We demonstrated significant differences in dog leukocyte antigen (DLA) haplotype frequencies in six of these breeds: Cocker spaniel, Springer spaniel, Labrador, West Highland white terrier (WHWT), Bearded collie, and Standard poodle. In the Springer spaniel, the DLA-DRB1*015:01--DQA1*006:01--DQB1*023:01 haplotype was significantly associated with disease risk (p = 0.014, odds ratio (OR) = 5.14) and showed a similar trend in the Cocker spaniel. This haplotype is related to one associated with hypoadrenocorticism in the Nova Scotia duck tolling retriever. Similar haplotypes shared between breeds were demonstrated, with DLA-DRB1*001:01--DQA1*001:01--DQB1*002:01 more prevalent in both affected Labrador (p = 0.0002, OR = 3.06) and WHWT (p = 0.01, OR = 2.11). Other haplotypes that have not previously been associated with the disease were identified. The inter-breed differences in DLA haplotypes associated with susceptibility to canine hypoadrenocorticism could represent divergent aetiologies. This could have implications for clinical diagnosis and future comparative studies. Alternatively, it may suggest that the gene of interest is closely linked to the MHC.

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.

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.

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.