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.

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.

Major histocompatibility complex and other allergy-related candidate genes associated with insect bite hypersensitivity in Icelandic horses.

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses caused by bites of insects. IBH is a multifactorial disease with contribution of genetic and environmental factors. Candidate gene association analysis of IBH was performed in a group of 89 Icelandic horses all born in Iceland and imported to Europe. Horses were classified in IBH-affected and non-affected based on clinical signs and history of recurrent dermatitis, and on the results of an in vitro sulfidoleukotriene (sLT)-release assay with Culicoides nubeculosus and Simulium vittatum extract. Different genetic markers were tested for association with IBH by the Fisher's exact test. The effect of the major histocompatibility complex (MHC) gene region was studied by genotyping five microsatellites spanning the MHC region (COR112, COR113, COR114, UM011 and UMN-JH34-2), and exon 2 polymorphisms of the class II Eqca-DRA gene. Associations with Eqca-DRA and COR113 were identified (p < 0.05). In addition, a panel of 20 single nucleotide polymorphisms (SNPs) in 17 candidate allergy-related genes was tested. During the initial screen, no marker from the panel was significantly (p < 0.05) associated with IBH. Five SNPs associated with IBH at p < 0.10 were therefore used for analysis of combined genotypes. Out of them, SNPs located in the genes coding for the CD14 receptor (CD14), interleukin 23 receptor (IL23R), thymic stromal lymphopoietin (TSLP) and transforming growth factor beta 3 (TGFB3) molecules were associated with IBH as parts of complex genotypes. These results are supported by similar associations and by expression data from different horse populations and from human studies.

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 and conservation in a small endangered horse population.

The Old Kladruber horses arose in the 17th century as a breed used for ceremonial purposes. Currently, grey and black coat colour varieties exist as two sub-populations with different recent breeding history. As the population underwent historical bottlenecks and intensive inbreeding, loss of genetic variation is considered as the major threat. Therefore, genetic diversity in neutral and non-neutral molecular markers was examined in the current nucleus population. Fifty microsatellites, 13 single nucleotide polymorphisms (SNPs) in immunity-related genes, three mutations in coat colour genes and one major histocompatibility (MHC-DRA) gene were studied for assessing genetic diversity after 15 years of conservation. The results were compared to values obtained in a similar study 13 years ago. The extent of genetic diversity of the current population was comparable to other breeds, despite its small size and isolation. The comparison between 1997 and 2010 did not show differences in the extent of genetic diversity and no loss of allele richness and/or heterozygosity was observed. Genetic differences identified between the black and grey sub-populations observed 13 years ago persisted. Deviations from the Hardy-Weinberg equilibrium found in 19 microsatellite loci and in five SNP loci are probably due to selective breeding. No differences between neutral and immunity-related markers were found. No changes in the frequencies of markers associated with two diseases, melanoma and insect bite hypersensitivity, were observed, due probably to the short interval of time between comparisons. It, thus, seems that, despite its small size, previous bottlenecks and inbreeding, the molecular variation of Old Kladruber horses is comparable to other horse breeds and that the current breeding policy does not compromise genetic variation of this endangered population.

Immunogenomic analysis of insect bite hypersensitivity in a model horse population.

Equine insect bite hypersensitivity (IBH) is a seasonal IgE-mediated dermatosis caused by bites of insects of the genus Culicoides. A familial predisposition for the disease has been shown but, except for the MHC, the genes involved have not been identified so far. An immunogenomic analysis of IBH was performed in a model population of Old Kladruby horses, all living in the same environment. Clinical signs of IBH were used as phenotypic manifestation of IBH. Furthermore, total serum IgE levels were determined in the sera of these horses and used as an independent phenotypic marker for the immunogenetic analysis. Single nucleotide polymorphisms (SNPs) in candidate immunity-related genes were used for association analyses. Genotypes composed of two to five genes encoding interferon gamma -IFNG, transforming growth factor beta 1 -TGFB1, Janus kinase 2 -JAK2, thymic stromal lymphopoietin -TSLP, and involucrin -IVL were associated with IBH, indicating a role of the genes in the pathogenesis of IBH. These findings were supported by analysis of gene expression in skin biopsies of 15 affected and 15 unaffected horses. Two markers associated with IBH, IFNG and TGFB1, showed differences in mRNA expression in skin biopsies from IBH-affected and non-affected horses (p<0.05). Expression of the gene coding for the CD14 receptor molecule -CD14 was different in skin biopsies at p<0.06. When total IgE levels were treated as binary traits, genotypes of IGHE, ELA-DRA, and IL10/b were associated with this trait. When treated as a continuous trait, total IgE levels were associated with genes IGHE, FCER1A, IL4, IL4R, IL10, IL1RA, and JAK2. This first report on non-MHC genes associated with IBH in horses is thus supported by differences in expression of genes known to play a role in allergy and immunity.