Hidden MHC genetic diversity in the Iberian ibex (Capra pyrenaica)

Diversity and evolution of the MHC class II DRB gene in the Capra sibirica experienced a demographic fluctuation in China

The major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates, and their proteins play a critical role in adaptive immunity for defense against a variety of pathogens. MHC diversity was lost in many species after experiencing a decline in size. To understand the variation and evolution of MHC genes in the Siberian ibex, Capra sibirica, which has undergone a population decline, we analyzed the variation of the second exon of MHC class II DRB genes in samples collected from five geographic localities in Xinjiang, China, that belong to three diverged mitochondrial clades. Consequently, we identified a total of 26 putative functional alleles (PFAs) with 260 bp in length from 43 individuals, and found one (for 27 individuals) to three (for 5 individuals) PFAs per individual, indicating the presence of one or two DRB loci per haploid genome. The Casi-DRB1*16 was the most frequently occurring PFA, Casi-DRB1*22 was found in only seven individuals, 14 PFAs occurred once, 7 PFAs twice, implying high frequency of rare PFAs. Interestingly, more than half (15) of the PFAs were specific to clade I, only two and three PFAs were specific to clades II and III, respectively. So, we assume that the polygamy and sexual segregation nature of this species likely contributed to the allelic diversity of DRB genes. Genetic diversity indices showed that PFAs of clade II were lower in nucleotide, amino acid, and supertype diversity compared to those of the other two clades. The pattern of allele sharing and FST values between the three clades was to some extent in agreement with the pattern observed in mitochondrial DNA divergence. In addition, recombination analyses revealed no evidence for significant signatures of recombination events. Alleles shared by clades III and the other two clades diverged 6 million years ago, and systematic neighbor grids showed Trans-species polymorphism. Together with the PAML and MEME analyses, the results indicated that the DRB gene in C. sibirica evolved under balancing and positive selection. However, by comparison, it can be clearly seen that different populations were under different selective pressures. Our results are valuable in understanding the diversity and evolution of the DRB gene in a mountain living C. sibirica and in making decisions on future long-term protection strategies.

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

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

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.

Assessing the levels of intraspecific admixture and interspecific hybridization in Iberian wild goats (Capra pyrenaica)

Iberian wild goats (Capra pyrenaica, also known as Iberian ibex, Spanish ibex, and Spanish wild goat) underwent strong genetic bottlenecks during the 19th and 20th centuries due to overhunting and habitat destruction. From the 1970s to 1990s, augmentation translocations were frequently carried out to restock Iberian wild goat populations (very often with hunting purposes), but they were not systematically planned or recorded. On the other hand, recent data suggest the occurrence of hybridization events between Iberian wild goats and domestic goats (Capra hircus). Augmentation translocations and interspecific hybridization might have contributed to increase the diversity of Iberian wild goats. With the aim of investigating this issue, we have genotyped 118 Iberian wild goats from Tortosa-Beceite, Sierra Nevada, Muela de Cortes, Gredos, Batuecas, and Ordesa and Monte Perdido by using the Goat SNP50 BeadChip (Illumina). The analysis of genotypic data indicated that Iberian wild goat populations are strongly differentiated and display low diversity. Only three Iberian wild goats out from 118 show genomic signatures of mixed ancestry, a result consistent with a scenario in which past augmentation translocations have had a limited impact on the diversity of Iberian wild goats. Besides, we have detected eight Iberian wild goats from Tortosa-Beceite with signs of domestic goat introgression. Although rare, hybridization with domestic goats could become a potential threat to the genetic integrity of Iberian wild goats; hence, measures should be taken to avoid the presence of uncontrolled herds of domestic or feral goats in mountainous areas inhabited by this iconic wild ungulate.

A robust sequencing assay of a thousand amplicons for the high-throughput population monitoring of Alpine ibex immunogenetics

Polymorphism for immune functions can explain significant variation in health and reproductive success within species. Drastic loss in genetic diversity at such loci constitutes an extinction risk and should be monitored in species of conservation concern. However, effective implementations of genome‐wide immune polymorphism sets into high‐throughput genotyping assays are scarce. Here, we report the design and validation of a microfluidics‐based amplicon sequencing assay to comprehensively capture genetic variation in Alpine ibex (Capra ibex). This species represents one of the most successful large mammal restorations recovering from a severely depressed census size and a massive loss in diversity at the major histocompatibility complex (MHC). We analysed 65 whole‐genome sequencing sets of the Alpine ibex and related species to select the most representative markers and to prevent primer binding failures. In total, we designed ~1,000 amplicons densely covering the MHC, further immunity‐related genes as well as randomly selected genome‐wide markers for the assessment of neutral population structure. Our analysis of 158 individuals shows that the genome‐wide markers perform equally well at resolving population structure as RAD‐sequencing or low‐coverage genome sequencing data sets. Immunity‐related loci show unexpectedly high degrees of genetic differentiation within the species. Such information can now be used to define highly targeted individual translocations. Our design strategy can be realistically implemented into genetic surveys of a large range of species. In conclusion, leveraging whole‐genome sequencing data sets to design targeted amplicon assays allows the simultaneous monitoring of multiple genetic risk factors and can be translated into species conservation recommendations.

Genetic epidemiology of the Alpine ibex reservoir of persistent and virulent brucellosis outbreak

While it is now broadly accepted that inter-individual variation in the outcomes of host-pathogen interactions is at least partially genetically controlled, host immunogenetic characteristics are rarely investigated in wildlife epidemiological studies. Furthermore, most immunogenetic studies in the wild focused solely on the major histocompatibility complex (MHC) diversity despite it accounts for only a fraction of the genetic variation in pathogen resistance. Here, we investigated immunogenetic diversity of the Alpine ibex (Capra ibex) population of the Bargy massif, reservoir of a virulent outbreak of brucellosis. We analysed the polymorphism and associations with disease resistance of the MHC Class II Drb gene and several non-MHC genes (Toll-like receptor genes, Slc11A1) involved in the innate immune response to Brucella in domestic ungulates. We found a very low neutral genetic diversity and a unique MHC Drb haplotype in this population founded few decades ago from a small number of individuals. By contrast, other immunity-related genes have maintained polymorphism and some showed significant associations with the brucellosis infection status hence suggesting a predominant role of pathogen-mediated selection in their recent evolutionary trajectory. Our results highlight the need to monitor immunogenetic variation in wildlife epidemiological studies and to look beyond the MHC.