Genomic analysis of Ovis aries (Ovar) MHC class IIa loci

Genome-wide association with footrot in hair and wool sheep

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum. Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH. We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1. These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete.

Placental Transcriptome Adaptations to Maternal Nutrient Restriction in Sheep

Placental development is modified in response to maternal nutrient restriction (NR), resulting in a spectrum of fetal growth rates. Pregnant sheep carrying singleton fetuses and fed either 100% (n = 8) or 50% (NR; n = 28) of their National Research Council (NRC) recommended intake from days 35–135 of pregnancy were used to elucidate placentome transcriptome alterations at both day 70 and day 135. NR fetuses were further designated into upper (NR NonSGA; n = 7) and lower quartiles (NR SGA; n = 7) based on day 135 fetal weight. At day 70 of pregnancy, there were 22 genes dysregulated between NR SGA and 100% NRC placentomes, 27 genes between NR NonSGA and 100% NRC placentomes, and 22 genes between NR SGA and NR NonSGA placentomes. These genes mediated molecular functions such as MHC class II protein binding, signaling receptor binding, and cytokine activity. Gene set enrichment analysis (GSEA) revealed significant overrepresentation of genes for natural-killer-cell-mediated cytotoxicity in NR SGA compared to 100% NRC placentomes, and alterations in nutrient utilization pathways between NR SGA and NR NonSGA placentomes at day 70. Results identify novel factors associated with impaired function in SGA placentomes and potential for placentomes from NR NonSGA pregnancies to adapt to nutritional hardship.

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.

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.

Characterisation of major histocompatibility complex class IIa haplotypes in an island sheep population

The ovine MHC class IIa is known to consist of six to eight loci located in close proximity on chromosome 20, forming haplotypes that are typically inherited without recombination. Here, we characterise the class IIa haplotypes within the Soay sheep (Ovis aries) on St. Kilda to assess the diversity present within this unmanaged island population. We used a stepwise sequence-based genotyping strategy to identify alleles at seven polymorphic MHC class IIa loci in a sample of 118 Soay sheep from four cohorts spanning 15 years of the long-term study on St. Kilda. DRB1, the most polymorphic MHC class II locus, was characterised first in all 118 sheep and identified six alleles. Using DRB1 homozygous animals, the DQA (DQA1, DQA2 and DQA2-like) and DQB (DQB1, DQB2 and DQB2-like) loci were sequenced, revealing eight haplotypes. Both DQ1/DQ2 and DQ2/DQ2-like haplotype configurations were identified and a single haplotype carrying three DQB alleles. A test sample of 94 further individuals typed at the DRB1 and DQA loci found no exceptions to the eight identified haplotypes and a haplotype homozygosity of 21.3%. We found evidence of historic positive selection at DRB1, DQA and DQB. The limited variation at MHC class IIa loci in Soay sheep enabled haplotype characterisation but showed that no single locus could capture the full extent of the expressed variation in the region.

Allelic nomenclature for the duplicated MHC class II DQ genes in sheep

The principal MHC class II molecules involved in the presentation of peptides to the antigen specific receptors on CD4⁺ T cells genes in sheep are derived from DR and DQ genes. Allelic nomenclature systems for the DRB1 and its partner DRA loci are available for Ovid's; however, no official nomenclature is available for the DQ genes which creates ambiguity within the research community. Ovine MHC haplotypes include at least two pairs of DQA and DQB genes, termed DQA1, DQB1 and DQA2, DQB2 and both sets are polymorphic and both seem to be functional. In a number of haplotypes, the DQA1 locus appears to be absent (DQA1-null) and is replaced by a second locus termed DQA2-like. Here, we identify families of alleles based on sequence similarity and phylogenetic clustering which correspond to each of the DQA and DQB genes identified in previous genomic and transcript analyses of homozygous animals. Using such criteria to cluster sequences, we have named 82 full-length and partial cDNA transcripts derived from domestic sheep (Ovis aries) which correspond to alleles at the Ovar-DQA1, DQA2, DQA2-like, DQB1, DQB2 and DQB2-like genes and provide associated sequence resources available to the research community through the IPD-MHC Database. This sets the basis for naming and annotation of DQ genes within the ovine MHC and may be used as a template for DQ genes in other ruminant species which will ultimately support research in livestock infectious disease.

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.

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.

Conserved haplotype blocks within the sheep MHC and low SNP heterozygosity in the Class IIa subregion

This report describes single-nucleotide polymorphisms (SNPs) in the sheep major histocompatibility complex (MHC) class II and class III regions and provides insights into the internal structure of this important genomic complex. MHC haplotypes were deduced from sheep family trios based on genotypes from 20 novel SNPs representative of the class II region and 10 previously described SNPs spanning the class III region. All 30 SNPs exhibited Hardy-Weinberg proportions in the sheep population studied. Recombination within an extended sire haplotype was observed within the class II region for 4 of 20 sheep chromosomes, thereby supporting the presence of separated IIa and IIb subregions similar to those present in cattle. SNP heterozygosity varied across the class II and III regions. One segment of the class IIa subregion manifested very low heterozygosity for several SNPs spanning approximately 120 Kbp. This feature corresponds to a subregion within the human MHC class II region previously described as a 'SNP desert' because of its paucity of SNPs. Linkage disequilibrium (LD) was reduced at the junction separating the putative class IIb and IIa subregions and also between the class IIa and the class III subregions. The latter observation is consistent with either an unmapped physical separation at this location or more likely a boundary characterized by more frequent recombination between two conserved subregions, each manifesting high within-block LD. These results identify internal blocks of loci in the sheep MHC, within which recombination is relatively rare.

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 novel HURRAH protocol reveals high numbers of monomorphic MHC class II loci and two asymmetric multi-locus haplotypes in the Père David's deer

The Père David's deer is a highly inbred, but recovered, species, making it interesting to consider their adaptive molecular evolution from an immunological perspective. Prior to this study, genomic sequencing was the only method for isolating all functional MHC genes within a certain species. Here, we report a novel protocol for isolating MHC class II loci from a species, and its use to investigate the adaptive evolution of this endangered deer at the level of multi-locus haplotypes. This protocol was designated “HURRAH” based on its various steps and used to estimate the total number of MHC class II loci. We confirmed the validity of this novel protocol in the giant panda and then used it to examine the Père David's deer. Our results revealed that the Père David's deer possesses nine MHC class II loci and therefore has more functional MHC class II loci than the eight genome-sequenced mammals for which full MHC data are currently available. This could potentially account at least in part for the strong survival ability of this species in the face of severe bottlenecking. The results from the HURRAH protocol also revealed that: (1) All of the identified MHC class II loci were monomorphic at their antigen-binding regions, although DRA was dimorphic at its cytoplasmic tail; and (2) these genes constituted two asymmetric functional MHC class II multi-locus haplotypes: DRA1*01 ∼ DRB1 ∼ DRB3 ∼ DQA1 ∼ DQB2 (H1) and DRA1*02 ∼ DRB2 ∼ DRB4 ∼ DQA2 ∼ DQB1 (H2). The latter finding indicates that the current members of the deer species have lost the powerful ancestral MHC class II haplotypes of nine or more loci, and have instead fixed two relatively weak haplotypes containing five genes. As a result, the Père David's deer are currently at risk for increased susceptibility to infectious pathogens.

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.

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.

Common promoter deletion is associated with 3.9-fold differential transcription of ovine CCR5 and reduced proviral level of ovine progressive pneumonia virus

Chemokine (C-C motif) Receptor 5 (CCR5) is a chemokine receptor that regulates immune cell recruitment in inflammation and serves as a coreceptor for human immunodeficiency virus (HIV). A human CCR5 coding deletion (termed delta-32) results in strong resistance to HIV infection, and sequence variants in CCR5 regulatory regions have been implicated in delayed progression to acquired immune deficiency syndrome. Both ovine progressive pneumonia virus (OPPV), also known as maedi-visna, and HIV are macrophage-tropic lentiviruses, have similar genomic structures, and cause lifelong persistent host infection, suggesting CCR5 may have a role in regulating OPPV provirus levels. Therefore, the ovine CCR5 genomic sequence was determined, and sequence variants were obtained from the open reading frame and surrounding regulatory sites. One CCR5 variant contained a 4-base deletion within a binding site for octamer transcription factors in the promoter region. A test for differential transcription from each allele in heterozygous animals showed a 3.9-fold transcription difference (P < 0.0001). OPPV proviral levels were also measured in 351 naturally exposed Rambouillet, Polypay and Columbia sheep. Deletion homozygotes showed reduced OPPV proviral levels among these animals (P < 0.01). The association of this CCR5 promoter deletion with OPPV levels will need to be validated in additional populations before the deletion can be recommended for widespread use in marker-assisted selection. However, because of the large impact on transcription and because CCR5 has roles in inflammation, recruitment of effector cells, and cell-mediated immunity, this deletion may play a role in the control of infections of many diverse pathogens of sheep.