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

Diversity of major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes and their interactions in domestic horses

The immunogenome is the part of the genome that underlies immune mechanisms and evolves under various selective pressures. Two complex regions of the immunogenome, major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes, play an important role in the response to selective pressures of pathogens. Their importance is expressed by their genetic polymorphism at the molecular level, and their diversity associated with different types of diseases at the population level. Findings of associations between specific combinations of MHC/NKR haplotypes with different diseases in model species suggest that these gene complexes did not evolve independently. No such associations have been described in horses so far. The aim of the study was to detect associations between MHC and NKR gene/microsatellite haplotypes in three horse breed groups (Camargue, African, and Romanian) by statistical methods; chi-square test, Fisher's exact test, Pearson's goodness-of-fit test and logistic regression. Associations were detected for both MHC/NKR genes and microsatellites; the most significant associations were found between the most variable KLRA3 gene and the EQCA-1 or EQCA-2 genes. This finding supports the assumption that the KLRA3 is an important receptor for MHC I and that interactions of these molecules play important roles in the horse immunity and reproduction. Despite some limitations of the study such as low numbers of horses or lack of knowledge of the selected genes functions, the results were consistent across different statistical methods and remained significant even after overconservative Bonferroni corrections. We therefore consider them biologically plausible.

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.

Microsatellite markers of the major histocompatibility complex genomic region of domestic camels

We identified and characterized 11 polymorphic microsatellite markers suitable for routine testing (three in the MHC class I sub-region, four in MHC class II and four in the MHC class III sub-region) of dromedaries and Bactrian camels. In total, 38 dromedaries and 33 Bactrian camels were genotyped, and interspecific differences were observed in the numbers of alleles and in allelic frequencies, as well as in the observed heterozygosity. These loci may be used as markers to study the adaptive genetic diversity of the MHC region in Old World camels.

Inter- and intrabreed diversity of the major histocompatibility complex (MHC) in primitive and draft horse breeds

Background

Polymorphism of major histocompatibility complex (MHC) genes ensures effective immune responses against a wide array of pathogens. However, artificial selection, as performed in the case of domestic animals, may influence MHC diversity. Here, we investigate and compare the MHC diversity of three populations of horses, for which different breeding policies were applied, to evaluate the impact of artificial selection and the environment on MHC polymorphism.

Methods

Samples of DNA were taken from 100 Polish draft horses, 38 stabled Konik Polski horses and 32 semiferal Konik Polski horses. MHC alleles and haplotype diversity within and between these populations of horses was estimated from 11 MHC microsatellite loci.

Results

MHC diversity measured based on allelic richness, observed heterozygosity, expected heterozygosity and polymorphism content was similar across the MHC microsatellite loci in all three populations. The highest expected heterozygosity was detected in semiferal primitive horses (He = 0.74), while the lowest was calculated for draft horses (He = 0.65). In total, 203 haplotypes were determined (111 in Polish draft horses, 43 in semiferal Konik Polski horses and 49 in stabled Konik Polski horses), and four haplotypes were shared between the two populations of Koniks. None of these haplotypes were present in any of the previously investigated horse breeds. Intra-MHC recombination events were detected in all three populations. However, the population of semiferal Konik horses showed the highest recombination frequency among the three horse populations. In addition, three recombination events were detected.

Conclusions

These results showed that despite the different breeding policies, the MHC allele and haplotype diversity was similarly high in all three horse populations. Nevertheless, the proportion of new haplotypes in the offspring was the highest in semiferal Konik Polski horses, which indicates the influence of the environment on MHC diversity in horses. Thus, we speculate that the genetic makeup of the domestic horse MHC might be more strongly influenced by the environment than by artificial selection. Moreover, intra-MHC conversion, insertion, and deletion and intra-MHC recombination may be proposed as mechanisms underlying the generation of new MHC haplotypes in horses.

Associations between the presence of specific antibodies to the West Nile Virus infection and candidate genes in Romanian horses from the Danube delta

The West Nile virus (WNV) is a mosquito-borne flavivirus causing meningoencephalitis in humans and animals. Due to their particular susceptibility to WNV infection, horses serve as a sentinel species. In a population of Romanian semi-feral horses living in the Danube delta region, we have analyzed the distribution of candidate polymorphic genetic markers between anti WNV-IgG seropositive and seronegative horses. Thirty-six SNPs located in 28 immunity-related genes and 26 microsatellites located in the MHC and LY49 complex genomic regions were genotyped in 57 seropositive and 32 seronegative horses. The most significant association (p_corr < 0.0002) was found for genotypes composed of markers of the SLC11A1 and TLR4 genes. Markers of five other candidate genes (ADAM17, CXCR3, IL12A, MAVS, TNFA), along with 5 MHC class I and LY49-linked microsatellites were also associated with the WNV antibody status in this model horse population. The OAS1 gene, previously associated with WNV-induced clinical disease, was not associated with the presence of anti-WNV antibodies.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

Genetic susceptibility to West Nile virus infection in Camargue horses

West Nile virus (WNV) is a mosquito-borne zoonotic neurotropic virus capable to cause lethal meningoencephalitis (WNE) in infected hosts such as birds, horses, and humans. Due to their sensitivity, horses serve as sentinel species in areas at risk. We studied a population of Camargue horses living in Southern France in two zones with endemic WNV circulation where WNV outbreaks were recorded in 2000 and 2003-4. Two sets of microsatellite markers located in MHC and Ly49 genomic regions were genotyped as well as multiple SNPs in ten immunity-related candidate gene regions. Associations between genetic polymorphisms and resistance/susceptibility to WNE were tested. While single marker associations were weak, compound two-gene genotypes of SNPs located within the MAVS, NCR2 and IL-10 genes and microsatellites HMS082 and CZM013 were associated with susceptibility to WNE. Combinations of microsatellite markers CZM009, ABGe17402 and ABGe9019 were associated with simple seroconversion without clinical signs of WNE (resistance). In addition, a distribution of polymorphic markers between WNV-IgG seropositive horses and a control group of WNV-IgG seronegative horses was tested. One SNP in the OAS1 gene (NC_009151.3:g.21961328A>G) was significantly associated with the seropositive phenotype (p_corr = 0.023; OR = 40.5 CI (4.28; 383.26); RR = 8.18 CI (1.27; 52.89) in the Camargue breed. In compound genotypes, SNP markers for SLC11A1, MAVS, OAS1, TLR4, ADAM17 and NCR2 genes and ten microsatellites showed non-random distribution between seropositive and seronegative groups of horses. Further analysis of associated markers could contribute to our understanding of anti-WNV defense mechanisms in horses.

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).