Genomic organisation and allelic diversity within coding and non-coding regions of the Ovar-DRB1 locus

Exploring genetic diversity and variation of Ovar-DRB1 gene in Sudan Desert Sheep using targeted next-generation sequencing

INTRODUCTION

The Ovar-DRB1 gene, a crucial element of the Major Histocompatibility Complex (MHC) Class II region, initiates adaptive immunity by presenting antigens to T-cells. Genetic diversity in sheep, particularly in MHC Class II genes like Ovar-DRB1, directly influences the specturm of presented antigens impacting immune responses and disease susceptability. Understanding the allelic diversity of Ovar-DRB1 gene in Sudan Desert Sheep (SDS) is essential for uncovering the genetic basis of immune responses and disease resistance, given the the breeds significance in Sudan's unique environment.

METHODS

Utilizing Targeted Next-Generation Sequencing (NGS) we explore allelic diversity in Ovar-DRB1 gene within SDS. Successfully ampliying and and sequencing the second exon of this gene in 288 SDS samples representing six breeds provided a comprehensive allelic profile, enabling a detalied examination of the gene's genetic makeup.

RESULTS

We identifed forty-six alleles, including four previously unreported, enrichness the genetic diversity of SDS breeds. These alleles exhibiting non-uniform distribution, varying frequencies across breeds, indicating a breed-specific genetic landscape. Certain alleles, known and novel, show higher frequencies in specific populations, suggesting potential associations with adaptive immune responses. Identifying these alleles sets the stage for investigating their functional roles and implications for disease resistance. Genetic differentiation among SDS breeds, as indicated by FST values and clustering analyses, highlights a unique genetic makeup shaped by geographic and historical factors. These differentiation patterns among SDS breeds have broader implications for breed conservation and targeted breeding to enhance disease resistance in specific populations.

CONCLUSION

This study unveils Ovar-DRB1 gene allelic diversity in SDS breeds through targeted NGS and genetic analyses, revealing new alleles that underscore the breeds' unique genetic profile. Insights into the genetic factors governing immune responses and disease resistance emerge, promising for optimization of breeding strategies for enhanced livestock health in Sudan's unique environment.

Association of Ovar-DRB1 alleles with innate immune responses in sheep

Background

Major histocompatibility complex (MHC) is the best characterised genetic region associated with adaptive immune responses, including humoral and cell‐mediated immunities.

Objectives

In this study, the association of MHC class II alleles with inflammatory cytokines and acute‐phase proteins was evaluated in sheep population.

Methods

Allelic diversity of second exon of ovine DRB1 locus (Ovar‐DRB1.2) was determined in 100 indigenous Iranian Lori‐Bakhtiari fat‐tailed sheep using restriction fragment length polymorphism and direct sequencing methods. The association of DRB1.2 alleles with inflammatory cytokines (interleukin‐1β, IL‐1β; IL‐6 and tumour necrosis factor‐α) and acute‐phase proteins (serum amyloid A, alpha‐1‐acid glycoprotein and haptoglobin) was examined using generalised linear model and multivariate regression analysis.

Results

Seven distinct RsaI restriction patterns and fourteen alleles were identified in this population. Allele DRB1*2101 showed a negative influence on the IL‐6 response and was associated with lower serum level of IL‐6. DRB1.2 heterozygous individuals also showed higher haptoglobin concentration than homozygotes.

Conclusions

These results provide additional support for the association between Ovar‐DRB1 alleles and regulation of immune responses in sheep population. Description of MHC polymorphism and its role in the controlling of immune responses will increase our understanding of host–pathogen interactions, and ultimately facilitate the selection of disease‐resistant flocks in genetic breeding programs.

Documentation of extensive genetic diversity in the Ovar-DRB1 gene in native Turkish sheep

Indigenous breeds have a high level of genetic diversity that might contribute to develop animal breeds with desired traits such as disease resistance and high productivity. Major histocompatibility complex (MHC) is a key component of adaptive immune system and consists of highly polymorphic genes that take part in adaptive immune response and disease resistance. Exploring and understanding the effect of polymorphisms in MHC could be beneficial to future animal breeding strategies. In this study, we sequenced the highly polymorphic Exon2 of the ovine DRB1 gene using Sanger sequencing to explore the diversity of this gene in six indigenous Turkish sheep breeds and two crossbreeds. In total, 894 haplotypes from 447 sheep were investigated, and 69 different haplotypes including 27 novel ones were identified. Among the identified haplotypes there were common and breed specific haplotypes. There was a relatively high diversity of the alleles within indigenous breeds. Allelic diversity patterns were mostly associated with geographical differences. The results of this study highlight the genetic variation within indigenous breeds which has important implications for biodiversity and the adaptability of breeds to specific environments. There is value to further studies which include other genomic regions and traits, and these could guide breeding strategies.

A chromosome-scale reference genome and genome-wide genetic variations elucidate adaptation in yak

Yak is an important livestock animal for the people indigenous to the harsh, oxygen‐limited Qinghai‐Tibetan Plateau and Hindu Kush ranges of the Himalayas. The yak genome was sequenced in 2012, but its assembly was fragmented because of the inherent limitations of the Illumina sequencing technology used to analyse it. An accurate and complete reference genome is essential for the study of genetic variations in this species. Long‐read sequences are more complete than their short‐read counterparts and have been successfully applied towards high‐quality genome assembly for various species. In this study, we present a high‐quality chromosome‐scale yak genome assembly (BosGru_PB_v1.0) constructed with long‐read sequencing and chromatin interaction technologies. Compared to an existing yak genome assembly (BosGru_v2.0), BosGru_PB_v1.0 shows substantially improved chromosome sequence continuity, reduced repetitive structure ambiguity, and gene model completeness. To characterize genetic variation in yak, we generated de novo genome assemblies based on Illumina short reads for seven recognized domestic yak breeds in Tibet and Sichuan and one wild yak from Hoh Xil. We compared these eight assemblies to the BosGru_PB_v1.0 genome, obtained a comprehensive map of yak genetic diversity at the whole‐genome level, and identified several protein‐coding genes absent from the BosGru_PB_v1.0 assembly. Despite the genetic bottleneck experienced by wild yak, their diversity was nonetheless higher than that of domestic yak. Here, we identified breed‐specific sequences and genes by whole‐genome alignment, which may facilitate yak breed identification.

Identification of the amino acids in the Major Histocompatibility Complex class II region of Scottish Blackface sheep that are associated with resistance to nematode infection

Lambs with the Major Histocompatibility Complex DRB1*1101 allele have been shown to produce fewer nematode eggs following natural and deliberate infection. These sheep also possess fewer adult Teladorsagia circumcincta than sheep with alternative alleles at the DRB1 locus. However, it is unclear if this allele is responsible for the reduced egg counts or merely acts as a marker for a linked gene. This study defined the MHC haplotypes in a population of naturally infected Scottish Blackface sheep by PCR amplification and sequencing, and examined the associations between MHC haplotypes and faecal egg counts by generalised linear mixed modelling. The DRB1*1101 allele occurred predominately on one haplotype and a comparison of haplotypes indicated that the causal mutation or mutations occurred in or around this locus. Additional comparisons with another resistant haplotype indicated that mutations in or around the DQB2*GU191460 allele were also responsible for resistance to nematode infections. Further analyses identified six amino acid substitutions in the antigen binding site of DRB1*1101 that were significantly associated with reductions in the numbers of adult T. circumcincta.

Molecular Cloning and Bioinformatics Analysis of DQA Gene from Mink (Neovison vison)

In the present study, we cloned, sequenced, and explored the structural and functional characteristics of the major histocompatibility complex (MHC)-DQA gene from mink (Neovison vison) for the first time. The full-length sequence of DQA gene was 1147-bp-long, contained a coding region of 768-bp, which was predicted to encoding 255 amino acid residues. The comparison between DQA from mink (Neovison vison) and other MHC-DQA molecules from different animal species showed that nucleotide and encoded amino acid sequences of the mink DQA gene exhibited high similarity with the ferret (Mustela pulourius furo). Phylogenetic analysis revealed that mink (Neovison vison) DQA is grouped with that of ferret (Mustela pulourius furo). The cloned sequence contained a 23-amino acid NH₂-terminal signal sequence with the signal peptide cutting site located in amino acids 23–24, and had three Asn-Xaa-Ser/Thr sequons. Three cysteine residues were also identified (Cys-85, Cys-121, and Cys-138). The 218 to 240 amino acids were predicted to be the transmembrane domains. The prediction of the secondary structure revealed three α-helixes and fourteen β-sheets in Neovison vison DQA protein, while random coil was a major pattern. In this study, the whole CDS sequence of Neovison vison DQA gene was successfully cloned, which was valuable for exploring the function and antiviral molecular mechanisms underlying the molecule. The findings of the present study have laid the foundation for the disease resistance and breeding of mink.

Genetic characterization of some Turkish sheep breeds based on the sequencing of the Ovar-DRB1 gene in the major histocompatibility complex (MHC) gene region

In this research, Ovar-DRB1 gene in the major histocompatibility complex (MHC) gene region was surveyed by DNA sequencing in some of the native sheep breeds that are reared in Turkey. A total of 80 samples were collected from eight different Turkish native sheep breeds, and these samples were used for DNA sequencing.

The exon 2 region of Ovar-DRB1 in the MHC gene region was polymerase chain reaction (PCR) amplified and sequenced. A total of 25 new alleles were revealed in the Ovar-DRB1 gene in Turkish native sheep breeds with 24 variable sites; only 13 sites were parsimony informative. The average pairwise genetic distance was 0.029 % for the Ovar-DRB1 gene exon 2 region. The sequence variations at eight different positions (7026, 7036, 7040, 7053, 7059, 7069, 7131 and 7214) are found in all of the studied samples. G → C transversion at position 7081 is only seen in Akkaraman sheep breed, whereas T → C transition at position 7097 is only seen in one sample from the Akkaraman breed. Overall, two main groups were detected among the 25 alleles from Turkish native sheep breeds. All Daǧliç and Kivircik alleles and one allele from Karayaka, Malya and Sakiz are grouped together while all the other breeds are grouped in the other branch.

A novel in silico minigene vaccine based on CD4+ T-helper and B-cell epitopes of EG95 isolates for vaccination against cystic echinococcosis

EG95 oncospheral antigen plays a crucial role in Echinococcus granulosus pathogenicity. Considering the diversity of antigen among different EG95 isolates, it seems to be an ideal antigen for designing a universal multivalent minigene vaccine, so-called multi-epitope vaccine. This is the first in silico study to design a construct for the development of global EG95-based hydatid vaccine against E. granulosus in intermediate hosts. After antigen sequence selection, the three-dimensional structure of EG95 was modeled and multilaterally validated. The preliminary parameters for B-cell epitope prediction were implemented such as the possible transmembrane helix, signal peptide, post-translational modifications and allergenicity. The high ranked linear and conformational B-cell epitopes derived from several online web-servers (e.g., ElliPro, BepiPred v1.0, BcePred, ABCpred, SVMTrip, IEDB algorithms, SEPPA v2.0 and Discotope v2.0) were utilized for multiple sequence alignment and then for engineering the vaccine construct. T-helper based epitopes were predicted by molecular docking between the high frequent ovar class II allele (Ovar-DRB1*1202) and hexadecamer fragments of the EG95 protein. Having used the immune-informatics tools, we formulated the first EG95-based minigene vaccine based on T-helper epitope with high-binding affinity to the ovar MHC allele. This designed construct was analyzed for different physicochemical properties. It was also codon-optimized for high-level expression in Escherichia coli k12. Taken all, we propose the present in silico vaccine constructs as a promising platform for the generation of broadly protective vaccines for species and genus-specific immunization of the natural hosts of the parasite.

Structural and functional diversity arising from intra- and inter-haplotype combinations of duplicated DQA and B loci within the ovine MHC

In sheep, the A and B loci encoding the α and β chains of the classical class II MHC molecules are DRA and DRB and DQA and DQB. Previous analyses described the duplication of the DQA and DQB genes. The majority of haplotypes include DQA1 and DQA2 loci, however, in a number of haplotypes, DQA1 appears absent and these haplotypes have been described as DQA1 null. In these haplotypes, the DQA2 locus is found in combination with a second locus which appeared more closely related to DQA2 than DQA1, hence the description of this locus as DQA2-like. Here we combine our previous analysis of the DQA transcripts with an analysis of the associated DQB transcripts in ten haplotypes from MHC homozygous animals. This allows the potential for surface expression of different haplotype combinations of DQA and B genes and the functional significance of DQA2-like and its predicted DQB partner to be determined. Atypical DQB transcripts (DQB2-like) were identified in haplotypes classified as DQA1-null and conserved DQB2-like orthologues were identified in other Bovidae indicating trans-species conservation of the allelic lineage. Functional combinations detected by co-transfection of DQ1, DQ2 and DQ2-like genes demonstrates the potential for a wide range of DQ molecules derived from both intra- and inter-haplotype as well as inter-locus combinations. We provide evidence that DQA2-like and B2-like genes form an evolutionary conserved pair which generates structurally distinct class II molecules that are likely to present a distinct range of peptides to CD4⁺ T cells.

A quick and robust MHC typing method for free-ranging and captive primate species

Gene products of the major histocompatibility complex (MHC) of human and non-human primates play a crucial role in adaptive immunity, and most of the relevant genes not only show a high degree of variability (polymorphism) but also copy number variation (CNV) is observed. Due to this diversity, MHC proteins influence the capability of individuals to cope with various pathogens. MHC and/or MHC-linked gene products such as odorant receptor genes are thought to influence mate choice and reproductive success. Therefore, MHC typing of wild and captive primate populations is considered to be useful in conservation biology, which is, however, often hampered by the need of invasive and time-consuming methods. All intact Mhc-DRB genes in primates appear to possess a complex and highly divergent microsatellite, DRB-STR. A panel of 154 pedigreed olive baboons (Papio anubis) was examined for their DRB content by DRB-STR analysis of genomic DNA. Using the same methodology on DNA of feces samples, DRB variability of a silvery gibbon population (Hylobates moloch) (N = 24), an endangered species, could successfully be studied. In both species, length determination of the DRB-STR resulted in the definition of unique genotyping patterns that appeared to be specific for a certain chromosome. Moreover, the different STR lengths were shown to segregate with the allelic variation of the respective gene. The results obtained expand data gained previously on DRB-STR typing in macaques, great apes, and humans and strengthen the conclusion that this protocol is applicable in molecular ecology, conservation biology, and colony management, especially of endangered primate species.

The genetic architecture of the MHC class II region in British Texel sheep

Understanding the structure of the major histocompatibility complex, especially the number and frequency of alleles, loci and haplotypes, is crucial for efficient investigation of the way in which the MHC influences susceptibility to disease. Nematode infection is one of the most important diseases suffered by sheep, and the class II region has been repeatedly associated with differences in susceptibility and resistance to infection. Texel sheep are widely used in many different countries and are relatively resistant to infection. This study determined the number and frequency of MHC class II genes in a small flock of Texel sheep. There were 18 alleles at DRB1, 9 alleles at DQA1, 13 alleles at DQB1, 8 alleles at DQA2 and 16 alleles at DQB2. Several haplotypes had no detectable gene products at DQA1, DQB1 or DQB2, and these were defined as null alleles. Despite the large numbers of alleles, there were only 21 distinct haplotypes in the population. The relatively small number of observed haplotypes will simplify finding disease associations because common haplotypes provide more statistical power but complicate the discrimination of causative mutations from linked marker loci.

Effect of swine leukocyte antigen-DQA gene variation on diarrhea in Large White, Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA-DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA-DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA-DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.

The diversity of major histocompatibility complex class II DRB1 gene in sheep breeds from Xinjiang, China

Exon 2 of the ovine leukocyte antigen OLA-DRB1 locus was examined in sheep from the Xinjiang Karakul Ram and Bashibai populations, and three generations of hybrids were derived from a cross between Bashibai and Altai Argali wild sheep. This identified 12 novel alleles and 30 previously reported alleles. A neighbor-joining tree of the amino acid sequences of these 42 alleles revealed allelic clusters shared across the study populations. There were significant differences in allelic frequency between Karakul Ram and Bashibai sheep. DRB1*K18cC was the most frequent allele in Kararul Ram with a frequency of 21.2%, while DRB1*2F10c8 (13.2%) and DRB1*0803 (13.2%) were the most frequent alleles found in Bashibai sheep; the alleles DRB1*2F16c2, DRB1*1601, and DRB1*0803 occurred most frequently in F1, F2, and F3 populations, with frequencies of 17.6%, 14.3%, and 20%, respectively. Although many alleles were shared by Bashibai and hybrid sheep, some alleles differed between them, especially in the F1 generation of the Bashibai × Altai Argali cross. The hybrid-specific alleles indicated the introgression of Altai Argali alleles into hybrid flocks. A population tree based on the OLA-DRB1 allelic frequency in each population indicated that the Bashibai sheep and three hybrid populations were similar, with Karakul Ram being genetically distinct.

MHC-DRB1/DQB1 Gene Polymorphism and Its Association with Resistance/Susceptibility to Cystic Echinococcosis in Chinese Merino Sheep

The aim of this study was to analyze the relationship between polymorphism of the MHC-DRB1/DQB1 gene and its resistance to Cystic Echinococcosis (C.E), as well as to screen out the molecular genetic marker of antiechinococcosis in Chinese Merino sheep. The MHCII-DRB1/DQB1 exon 2 was amplified by polymerase chain reaction (PCR) from DNA samples of healthy and hydatidosis sheep. PCR products were characterized by restriction fragment length polymorphism (RFLP) technique. Five restriction enzymes (Mval, HaeIII, SacI, SacII, and Hin1I) were employed to cut DRB1, while seven restriction enzymes (MroxI, ScaI, SacII, NciI, TaqI, Mval, and HaeIII) were employed to cut DQB1.Results showed that frequencies of patterns Mvalbb (P < 0.01), SacIab in DRB1 exon 2 (P < 0.05), and TaqIaa, HaeIIInn (P < 0.01) in DQB1 exon 2 were significantly higher in the healthy group compared with the C.E individuals, which implied that there was a strong association between these genotypes and hydatidosis resistance or susceptibility. Chi-square test showed that individuals with the genic haplotype DRB1-SacIab/DRB1-Mvalbb/DQB1-TaqIaa/DQB1-HaeIIInn (P < 0.01) were relatively resistant to C.E, while individuals with the genic haplotypes DRB1-Mvalbc/DQB1-Mvalyy/DQB1-TaqIab/DQB1-HaeIIImn (P < 0.01) and DRB1-Mvalbb/DQB1-Mvalcc/DQB1-TaqIab/DQB1-HaeIIImn (P < 0.01) were more susceptible to C.E. In addition, to confirm these results, a fielding experiment was performed with Chinese Merino sheep which were artificially infected with E.g. The result was in accordance with the results of the first study. In conclusion, MHC-DRB1/DQB1 exon 2 plays an important role as resistant to C.E in Chinese Merino sheep. In addition, the molecular genetic marker of antiechinococcosis (DRB1-SacIab/DRB1-Mvalbb/DQB1-TaqIaa/DQB1-HaeIIInn) was screened out in Chinese Merino sheep.

A physical map of a BAC clone contig covering the entire autosome insertion between ovine MHC Class IIa and IIb

Background

The ovine Major Histocompatibility Complex (MHC) harbors genes involved in overall resistance/susceptibility of the host to infectious diseases. Compared to human and mouse, the ovine MHC is interrupted by a large piece of autosome insertion via a hypothetical chromosome inversion that constitutes ~25% of ovine chromosome 20. The evolutionary consequence of such an inversion and an insertion (inversion/insertion) in relation to MHC function remains unknown. We previously constructed a BAC clone physical map for the ovine MHC exclusive of the insertion region. Here we report the construction of a high-density physical map covering the autosome insertion in order to address the question of what the inversion/insertion had to do with ruminants during the MHC evolution.

Results

A total of 119 pairs of comparative bovine oligo primers were utilized to screen an ovine BAC library for positive clones and the orders and overlapping relationships of the identified clones were determined by DNA fingerprinting, BAC-end sequencing, and sequence-specific PCR. A total of 368 positive BAC clones were identified and 108 of the effective clones were ordered into an overlapping BAC contig to cover the consensus region between ovine MHC class IIa and IIb. Therefore, a continuous physical map covering the entire ovine autosome inversion/insertion region was successfully constructed. The map confirmed the bovine sequence assembly for the same homologous region. The DNA sequences of 185 BAC-ends have been deposited into NCBI database with the access numbers HR309252 through HR309068, corresponding to dbGSS ID 30164010 through 30163826.

Conclusions

We have constructed a high-density BAC clone physical map for the ovine autosome inversion/insertion between the MHC class IIa and IIb. The entire ovine MHC region is now fully covered by a continuous BAC clone contig. The physical map we generated will facilitate MHC functional studies in the ovine, as well as the comparative MHC evolution in ruminants.

A single nomenclature and associated database for alleles at the major histocompatibility complex class II DRB1 locus of sheep

The development of standardised nomenclatures with associated databases containing reference sequences for alleles at polymorphic loci within the major histocompatibility complex (MHC) has been facilitated by the development of the immuno polymorphism database (IPD). Recently, included within IPD-MHC is information on allelic diversity within sheep species (IPD-MHC-OLA). Here, we present the first report of progress in populating the sheep IPD-MHC database with alleles at the class II MHC DRB1 locus. The sequence of 63 Ovar-DRB1 alleles within 24 allelic families is now held within the database, each meeting the minimum requirement of a complete second exon. These sequences are derived from a combination of genomic and cDNA-based approaches and represent the most extensive collection of validated alleles at the sheep DRB1 locus yet described. Although these 63 alleles probably represent only a fraction of the DRB1 allelic diversity in sheep species worldwide, we encourage the research community to use the official allelic nomenclature and to contribute allelic sequences to the database via its web-based submission tool. In time, the IPD-MHC-OLA resource will underpin population-based MHC genotyping studies and help to simplify meta-analyses of multi-source data from wild and domestic sheep populations.

A complete DNA sequence map of the ovine major histocompatibility complex

Background

The ovine Major Histocompatibility Complex (MHC) harbors clusters of genes involved in overall resistance/susceptibility of an animal to infectious pathogens. However, only a limited number of ovine MHC genes have been identified and no adequate sequence information is available, as compared to those of swine and bovine. We previously constructed a BAC clone-based physical map that covers entire class I, class II and class III region of ovine MHC. Here we describe the assembling of a complete DNA sequence map for the ovine MHC by shotgun sequencing of 26 overlapping BAC clones.

Results

DNA shotgun sequencing generated approximately 8-fold genome equivalent data that were successfully assembled into a finished sequence map of the ovine MHC. The sequence map spans approximately 2,434,000 nucleotides in length, covering almost all of the MHC loci currently known in the sheep and cattle. Gene annotation resulted in the identification of 177 protein-coding genes/ORFs, among which 145 were not previously reported in the sheep, and 10 were ovine species specific, absent in cattle or other mammals. A comparative sequence analyses among human, sheep and cattle revealed a high conservation in the MHC structure and loci order except for the class II, which were divided into IIa and IIb subregions in the sheep and cattle, separated by a large piece of non-MHC autosome of approximately 18.5 Mb. In addition, a total of 18 non-protein-coding microRNAs were predicted in the ovine MHC region for the first time.

Conclusion

An ovine MHC DNA sequence map was successfully assembled by shotgun sequencing of 26 overlapping BAC clone. This makes the sheep the second ruminant species for which the complete MHC sequence information is available for evolution and functional studies, following that of the bovine. The results of the comparative analysis support a hypothesis that an inversion of the ancestral chromosome containing the MHC has shaped the MHC structures of ruminants, as we currently observed in the sheep and cattle. Identification of relative large numbers of microRNAs in the ovine MHC region helps to provide evidence that microRNAs are actively involved in the regulation of MHC gene expression and function.

Trans-species polymorphism and selection in the MHC class II DRA genes of domestic sheep

Highly polymorphic genes with central roles in lymphocyte mediated immune surveillance are grouped together in the major histocompatibility complex (MHC) in higher vertebrates. Generally, across vertebrate species the class II MHC DRA gene is highly conserved with only limited allelic variation. Here however, we provide evidence of trans-species polymorphism at the DRA locus in domestic sheep (Ovis aries). We describe variation at the Ovar-DRA locus that is far in excess of anything described in other vertebrate species. The divergent DRA allele (Ovar-DRA*0201) differs from the sheep reference sequences by 20 nucleotides, 12 of which appear non-synonymous. Furthermore, DRA*0201 is paired with an equally divergent DRB1 allele (Ovar-DRB1*0901), which is consistent with an independent evolutionary history for the DR sub-region within this MHC haplotype. No recombination was observed between the divergent DRA and B genes in a range of breeds and typical levels of MHC class II DR protein expression were detected at the surface of leukocyte populations obtained from animals homozygous for the DRA*0201, DRB1*0901 haplotype. Bayesian phylogenetic analysis groups Ovar-DRA*0201 with DRA sequences derived from species within the Oryx and Alcelaphus genera rather than clustering with other ovine and caprine DRA alleles. Tests for Darwinian selection identified 10 positively selected sites on the branch leading to Ovar-DRA*0201, three of which are predicted to be associated with the binding of peptide antigen. As the Ovis, Oryx and Alcelaphus genera have not shared a common ancestor for over 30 million years, the DRA*0201 and DRB1*0901 allelic pair is likely to be of ancient origin and present in the founding population from which all contemporary domestic sheep breeds are derived. The conservation of the integrity of this unusual DR allelic pair suggests some selective advantage which is likely to be associated with the presentation of pathogen antigen to T-cells and the induction of protective immunity.

Sequence-based genotyping of the sheep MHC class II DRB1 locus

The immunopolymorphism database (IPD) provides a single nomenclature for alleles at the major histocompatibility complex (MHC) loci for a range of different species. The minimum requirements for inclusion of a sheep class II DRB1 sequence is a submission that includes all polymorphic sites within the second exon from at least two independent polymerase chain reactions (PCR). In order to meet these requirements, we have developed a DNA-based genotyping method for the rapid analysis of allelic diversity at the DRB1 locus in domestic sheep, Ovis aries. Using a series of primers located within introns flanking exon 2 and genomic DNA from a cohort of 214 sheep representing 15 different breeds and crossbreeds, the complete exon 2 sequences of 38 Ovar-DRB1 alleles were obtained. This sequence resource allowed the development of a generic set of locus-specific primers which amplify a fragment that includes all polymorphic sites within the second exon. Bidirectional sequence analysis of the PCR product provides a composite sequence where each polymorphic site is represented by the corresponding International Union of Biochemistry nucleotide code. A Basic Local Alignment Search Tool search of alleles held within the IPD or National Center for Biotechnology Information databases allows individual allele sequences to be identified. Low levels of homozygosity (7.48%) within the cohort and verification of previously genotyped samples confirmed the broad allelic specificity of this method. It improves on currently available methods and is broadly applicable to the analysis of MHC diversity in studies investigating linkages with resistance or susceptibility to disease.