Genetic diversity of the class II major histocompatibility DRA locus in European, Asiatic and African domestic donkeys

Newly identified variability of the antigen binding site coding sequences of the equine major histocompatibility complex class I and class II genes

The major histocompatibility complex (MHC) with its class I and II genes plays a crucial role in the immune response to pathogens by presenting oligopeptide antigens to various immune response effector cells. In order to counteract the vast variability of infectious agents, MHC class I and II genes usually retain high levels of SNPs mainly concentrated in the exons encoding the antigen binding sites. The aim of the study was to reveal new variability of selected MHC genes with a special focus on MHC class I physical haplotypes. Long-range NGS to was used to identify exon 2-exon 3 alleles in three genetically distinct horse breeds. A total of 116 allelic variants were found in the MHC class I genes Eqca-1, Eqca-2, Eqca-7 and Eqca-Ψ, 112 of which were novel. The MHC class II DRA locus was confirmed to comprise five exon 2 alleles, and no new sequences were observed. Additional variability in terms of 15 novel exon 2 alleles was identified in the DQA1 locus. Extensive overall variability across the entire MHC region was confirmed by an analysis of MHC-linked microsatellite loci. Both diversifying and purifying selection were detected within the MHC class I and II loci analyzed.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae

Background

The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family.

Results

Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3.

Conclusions

The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Genetics, Evolution, and Physiology of Donkeys and Mules

The genus Equus is made up of donkeys, horses, and zebras. Despite significant variation in chromosome number across these species, interspecies breeding results in healthy, although infertile, hybrid offspring. Most notable among these are the horse-donkey hybrids, the mule and hinny. Donkeys presently are used for everything from companion animals to beasts of burden. Although closely related from an evolutionary standpoint, differences in anatomy and physiology preclude the assumption that they can be treated identically to the domestic horse. Veterinarians should be aware of these differences and adjust their practice accordingly.

Microsatellite markers for evaluating the diversity of the natural killer complex and major histocompatibility complex genomic regions in domestic horses

Genotyping microsatellite markers represents a standard, relatively easy, and inexpensive method of assessing genetic diversity of complex genomic regions in various animal species, such as the major histocompatibility complex (MHC) and/or natural killer cell receptor (NKR) genes. MHC-linked microsatellite markers have been identified and some of them were used for characterizing MHC polymorphism in various species, including horses. However, most of those were MHC class II markers, while MHC class I and III sub-regions were less well covered. No tools for studying genetic diversity of NKR complex genomic regions are available in horses. Therefore, the aims of this work were to establish a panel of markers suitable for analyzing genetic diversity of the natural killer complex (NKC), and to develop additional microsatellite markers of the MHC class I and class III genomic sub-regions in horses. Nine polymorphic microsatellite loci were newly identified in the equine NKC. Along with two previously reported microsatellites flanking this region, they constituted a panel of 11 loci allowing to characterize genetic variation in this functionally important part of the horse genome. Four newly described MHC class I/III-linked markers were added to 11 known microsatellites to establish a panel of 15 MHC markers with a better coverage of the class I and class III sub-regions. Major characteristics of the two panels produced on a group of 65 horses of 13 breeds and on five Przewalski's horses showed that they do reflect genetic variation within the horse species.

Genetic monitoring of horses in the Czech Republic: A large-scale study with a focus on the Czech autochthonous breeds

We propose the first comprehensive in-depth study monitoring horses in the Czech Republic. We scanned 9,289 animals from 44 populations for 17 equine STRs. Other equids analysed involved Equus przewalskii and Equus asinus. The total of 228 different alleles were detected, with the mean number of 13.4 per locus. The highest allelic richness (AR) was found in the Welsh Part Bred (6.01), followed by the Camargue (5.93) and Czech Sport Pony (5.91), whereas the Friesian exhibited the lowest AR (3.06). Interpopulation differences explained approximately nine per cent of the total genetic diversity. Reynold's genetic distance ranged from 0.003 between the Czech Warmblood and the Slovak Warmblood to 0.404 between the Friesian and donkeys. Close genetic proximity between the Silesian Noriker and Noriker was revealed. The Moravian Warmblood was better differentiated and more distant from the Czech Warmblood than the Kinsky Horse and retained the original genes of the old Austro-Hungarian tribes. A high gene flow level and a lack of genetic structure were found in the seven studied populations. Despite the historical bottlenecks and previous inbreeding, the Czech-Moravian Belgian Horse, Hucul, Old Kladruber Horse and Silesian Noriker did not suffer a serious loss of genetic diversity due to genetic drift/low effective population size. A NeighborNet dendrogram revealed breeds not classified in their groups according to the nomenclature (the Friesian, Hafling and Merens).

MHC haplotype diversity in Persian Arabian horses determined using polymorphic microsatellites

Previous research on the equine major histocompatibility complex (MHC) demonstrated strong correlations between haplotypes defined by polymorphic intra-MHC microsatellites and haplotypes defined using classical serology. Here, we estimated MHC diversity in a sample of 124 Arabian horses from an endangered strain native to Iran (Persian Asil Arabians), using a validated 10-marker microsatellite panel. In a group of 66 horses related as parent-offspring pairs or half-sibling groups, we defined 51 MHC haplotypes, 49 of which were new. In 47 of the remaining 58 unrelated horses, we could assign one previously identified MHC haplotype, and by default, we gave provisional haplotype status to the remaining constellation of microsatellite alleles. In these horses, we found 21 haplotypes that we had previously defined and 31 provisional haplotypes, two of which had been identified in an earlier study. This gave a total of 78 new MHC haplotypes. The final 11 horses were MHC heterozygotes that we could not phase using information from any of the previously validated or provisional haplotypes. However, we could determine that these horses carried a total of 22 different undefined haplotypes. In the overall population sample, we detected three homozygous horses and one maternally inherited recombinant from 21 informative segregations. Virtually all of the horses tested were MHC heterozygotes, and most unrelated horses (98%) were heterozygous for rare microsatellite-defined haplotypes found less than three times in the sampled horses. This is evidence for a very high level of MHC haplotype variation in the Persian Asil Arabian horse.

Polymorphism at expressed DQ and DR loci in five common equine MHC haplotypes

The polymorphism of major histocompatibility complex (MHC) class II DQ and DR genes in five common equine leukocyte antigen (ELA) haplotypes was determined through sequencing of mRNA transcripts isolated from lymphocytes of eight ELA homozygous horses. Ten expressed MHC class II genes were detected in horses of the ELA-A3 haplotype carried by the donor horses of the equine bacterial artificial chromosome (BAC) library and the reference genome sequence: four DR genes and six DQ genes. The other four ELA haplotypes contained at least eight expressed polymorphic MHC class II loci. Next generation sequencing (NGS) of genomic DNA of these four MHC haplotypes revealed stop codons in the DQA3 gene in the ELA-A2, ELA-A5, and ELA-A9 haplotypes. Few NGS reads were obtained for the other MHC class II genes that were not amplified in these horses. The amino acid sequences across haplotypes contained locus-specific residues, and the locus clusters produced by phylogenetic analysis were well supported. The MHC class II alleles within the five tested haplotypes were largely non-overlapping between haplotypes. The complement of equine MHC class II DQ and DR genes appears to be well conserved between haplotypes, in contrast to the recently described variation in class I gene loci between equine MHC haplotypes. The identification of allelic series of equine MHC class II loci will aid comparative studies of mammalian MHC conservation and evolution and may also help to interpret associations between the equine MHC class II region and diseases of the horse.

A field survey for the seroprevalence of Theileria equi and Babesia caballi in donkeys from Nuu Division, Kenya

Equine piroplasmosis is one of the most significant tick-borne disease of equids. The prevalence of this disease in donkeys of semi-arid Kenya remains largely unexplored. The primary objective of this study was to demonstrate the extent to which donkeys in Nuu division, Kenya have been exposed to the haemoprotozoans Babesia caballi and Theileria equi, the causative agents of equine piroplasmosis. The study also assessed the effect of age and sex on seroprevalence. A stratified sampling approach was used and three hundred and fourteen donkeys were sampled across nine sub-locations in Nuu division, Mwingi district. Serodiagnosis was via competitive inhibition enzyme linked immunosorbent assays (cELISA). The seroprevalence of T. equi was 81.2% (95% CI: 76.4-85.4). There was no significant difference in sub-location seropositivity, gender seropositivity or age related seropositivity. Antibodies against B. caballi were not detected (95% CI: 0-1.2). Findings from this study suggest that T. equi infection is endemic in Nuu division, Mwingi where it exists in a state of endemic stability. Existence of the infection should be communicated to animal health practitioners and donkey owning communities in the area.

Isolation and cDNA characteristics of MHC-DRA genes from gayal (Bos frontalis) and gaytle (Bos frontalis × Bos taurus)

The mammalian major histocompatibility complex (MHC) plays important roles in pathogen recognition and disease resistance. In the present study, the coding sequence and the 5′- and 3′-untranslated regions of MHC class II DR alpha chain (the DRA gene) from rare gayal and gaytle were cloned and analyzed to dissect structural and functional variations. The nucleotide and amino acid sequences for the DRA genes in gayal (Bofr-DRA) and gaytle (Bofr × BoLA-DRA) were almost identical to those for cattle and yak (99%). Compared to yak, two amino acids substitutions in the signal peptide (SP) domain for gayal were found within all Bos animals. Except for only one replacement in the amino acid within the α2 domain of the DRA protein in gayal, the additional residues were highly conserved across the species investigated. The 20 peptide-binding sites (PBS) of Bofr-DRA and Bofr × BoLA-DRA were essentially reserved in the α1 domain among all species investigated. The lesser degree of substitution in Bofr-DRA is concordant with the concept that the DRA gene is highly conserved among all mammals. The very high degree of conservativity of the DRA gene among ruminants, including gayal, suggests its recent evolutionary separation.

Molecular characteristics of the MHC-DRA genes from yak (Bos grunniens) and Chinese yakow (Bos grunniens × Bos taurus)

Two full-length cDNAs (762 bp) of the DRA gene from yak and Chinese yakow were isolated and analysed to identify structural and functional variations. The sequences for DRA in yak (Bogr-DRA) and Chinese yakow (Bogr × BoLA-DRA) were essentially identical to those for cattle (99%) and buffalo (97%). Except for two substitutions in the amino acids comprising the domain for signal peptide (SP) in yak, the additional residues were highly conserved across the species investigated. Peptide-binding site (PBS) of Bogr-DRA and Bogr × BoLA-DRA was highly reserved in the α1 domain among all species investigated. The lack of mutation in Bogr-DRA is consistent with the conception that the gene is highly conserved among all mammalian species. The very high conservation of the DRA gene among ruminants, including yak, may be due to its recent evolutionary detachment.

Unraveling the effects of selection and demography on immune gene variation in free-ranging plains zebra (Equus quagga) populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.