Closing a gap in the physical map of the ovine major histocompatibility complex

MHC class IIa haplotypes derived by high-throughput SNP screening in an isolated sheep population

Investigating the current evolutionary processes acting on a highly polymorphic gene region, such as the major histocompatibility complex (MHC), requires extensive population data for both genotypes and phenotypes. The MHC consists of several tightly linked loci with both allelic and gene content variation, making it challenging to genotype. Eight class IIa haplotypes have previously been identified in the Soay sheep (Ovis aries) of St. Kilda using Sanger sequencing and cloning, but no single locus is representative of all haplotypes. Here, we exploit the closed nature of the island population of Soay sheep and its limited haplotypic variation to identify a panel of SNPs that enable imputation of MHC haplotypes. We compared MHC class IIa haplotypes determined by Sanger sequence-based genotyping of 135 individuals to their SNP profiles generated using the Ovine Infinium HD BeadChip. A panel of 11 SNPs could reliably determine MHC diplotypes, and two additional SNPs within the DQA1 gene enabled detection of a recombinant haplotype affecting only the SNPs downstream of the expressed genes. The panel of 13 SNPs was genotyped in 5951 Soay sheep, of which 5349 passed quality control. Using the Soay sheep pedigree, we were able to trace the origin and inheritance of the recombinant SNP haplotype. This SNP-based method has enabled the rapid generation of locus-specific MHC genotypes for large numbers of Soay sheep. This volume of high-quality genotypes in a well-characterized population of free-living sheep will be valuable for investigating the mechanisms maintaining diversity at the MHC.

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.

[Screening of tissues pooled cDNA library using probes by restricted fragments of BAC positive clones of ovine MHC]

Under the premise what we have known bacterial artificial chromosome(BAC)clone sequence information and gene annotation predicted in the Chinese Merino sheep major histocompatibility complex (MHC) region, the digested fragments from 6 BAC clones that were located in the MHC region of the Chinese Merino sheep genome BAC library, which were used to screen the cDNA library using plaque in situ hybridization as probes. The full length of positive cDNA clones (sequences) isolated were completely sequenced, and the sequences obtained were aligned with the corresponding known sequence information and the BAC clones with gene annotation. Meanwhile, the sequence similarity was searched in NCBI Blastn database. This work aimed at verification of accuracy of the gene annotation results and initial analysis of gene (sequence) function. At last, 27 positive cDNA clones (sequences) in total were screened through two runs of hybridization. It was also found that these sequences could be positioned in the corresponding BAC clones, and 25 sequences were located in exon area of the annotated gene. It was verified that 23 sequences had the highest sequence similarity with those in the Bos taurus by searching against the NCBI Blastn database; moreover, the function of these sequences were closely relate to immunology.

[Confirmation and expression analysis of three predicted genes in sheep MHC region]

Previous DNA sequencing of BAC clones covering entire ovine MHC (OLA) region resulted in identification of approximately 130 functional genes in the region, of which 8 were predicted by computer software to be exclusively existed in sheep, but not in any other species known to date. In the present study, we successfully identified and cloned cDNA sequence of OaN2, OaN5, and OaN6 from representative sheep tissues, confirmed their existence in reality. The sequences obtained experimentally exactly identical to those predicted previously. The length of cDNA fragments for OaN2, OaN5, and OaN6 was 270 bp, 309 bp, and 205 bp, respectively, with GenBank accession number assigned as JF330782 (OaN2), JF330783 (OaN5), and JF330784 (OaN6). Northern analyses indicated that the mRNA transcripts of OaN2 were mainly seen in ovine mesenteric lymph nodes and spleen, while OaN5 was observed in only in mesenteric lymph nodes. In contrast, OaN6 transcripts were detected in all tissues except for liver and heart. Western blot showed that OaN2 protein expression level was detected in mesenteric lymph nodes, spleen, and liver, essentially consistent with that of mRNA transcripts. Immunohistochemistry analysis showed that OaN2 protein was highly expressed in ovine mesenteric lymph nodes, moderately expressed in, and not expressed in heart, liver, and pancreas, consistent with the results of Western blotting. The cloning and expression analysis of 3 novel genes provide a basis for revealing their specificities and would be helpful to further study of their expression profile and their potential functions.

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