Genetic and proteomic analysis of the MHC class I repertoire from four ovine haplotypes

Large scale transcriptional analysis of MHC class I haplotype diversity in sheep

Domestic sheep (Ovis aries) have been an important component of livestock agricultural production for thousands of years. Preserving genetic diversity within livestock populations maintains a capacity to respond to changing environments and rapidly evolving pathogens. MHC genetic diversity can influence immune functionality at individual and population levels. Here, we focus on defining functional MHC class I haplotype diversity in a large cohort of Scottish Blackface sheep pre-selected for high levels of MHC class II DRB1 diversity. Using high-throughput amplicon sequencing with three independent sets of barcoded primers we identified 134 MHC class I transcripts within 38 haplotypes. Haplotypes were identified with between two and six MHC class I genes, plus variable numbers of conserved sequences with very low read frequencies. One or two highly transcribed transcripts dominate each haplotype indicative of two highly polymorphic, classical MHC class I genes. Additional clusters of medium, low, and very low expressed transcripts are described, indicative of lower transcribed classical, non-classical and genes whose function remains to be determined.

Uncovering novel MHC alleles from RNA-Seq data: expanding the spectrum of MHC class I alleles in sheep

BACKGROUND

Major histocompatibility complex (MHC) class I glycoproteins present selected peptides or antigens to CD8 + T cells that control the cytotoxic immune response. The MHC class I genes are among the most polymorphic loci in the vertebrate genome, with more than twenty thousand alleles known in humans. In sheep, only a very small number of alleles have been described to date, making the development of genotyping systems or functional studies difficult. A cost-effective way to identify new alleles could be to use already available RNA-Seq data from sheep. Current strategies for aligning RNA-Seq reads against annotated genome sequences or transcriptomes fail to detect the majority of class I alleles. Here, I combine the alignment of RNA-Seq reads against a specific reference database with de novo assembly to identify alleles. The method allows the comprehensive discovery of novel MHC class I alleles from RNA-Seq data (DinoMfRS).

RESULTS

Using DinoMfRS, virtually all expressed MHC class I alleles could be determined. From 18 animals 75 MHC class I alleles were identified, of which 69 were novel. In addition, it was shown that DinoMfRS can be used to improve the annotation of MHC genes in the sheep genome sequence.

CONCLUSIONS

DinoMfRS allows for the first time the annotation of unknown, more divergent MHC alleles from RNA-Seq data. Successful application to RNA-Seq data from 16 animals has approximately doubled the number of known alleles in sheep. By using existing data, alleles can now be determined very inexpensively for populations that have not been well studied. In addition, MHC expression studies or evolutionary studies, for example, can be greatly improved in this way, and the method should be applicable to a broader spectrum of other multigene families or highly polymorphic genes.

Sequence diversity and differences at the highly duplicated MHC-I gene reflect viral susceptibility in sympatric pinniped species

Differences in disease susceptibility among species can result from rapid host-pathogen coevolution and differences in host species ecology that affect the strength and direction of natural selection. Among two sympatric pinniped species that differ in sociality and putative disease exposure, we investigate observed differences in susceptibility through an analysis of a highly variable, duplicated gene family involved in the vertebrate immune response. Using high-throughput amplicon sequencing, we characterize diversity at the two exons that encode the peptide binding region of the major histocompatibility complex class I (MHC-I) gene in harbor (N = 60) and gray (N = 90) seal populations from the Northwest Atlantic. Across species, we identified 106 full-length exon 2 and 103 exon 3 sequence variants and a minimum of 11 duplicated MHC-I loci. The sequence variants clustered in 15 supertypes defined by the physiochemical properties of the peptide binding region, including a putatively novel Northwest Atlantic MHC-I diversity sublineage. Trans-species polymorphisms, dN/dS ratios, and evidence of gene conversion among supertypes are consistent with balancing selection acting on this gene. High functional redundancy suggests particularly strong selection among gray seals at the novel Northwest Atlantic MHC-I diversity sublineage. At exon 2, harbor seals had a significantly greater number of variants per individual than gray seals, but fewer supertypes. Supertype richness and private supertypes are hypothesized to contribute to observed differences in disease resistance between species, as consistently, across the North Atlantic and many disease outbreaks, gray seals appear to be more resistant to respiratory viruses than harbor seals.

Immunological Homeostasis at the Ovine Placenta May Reflect the Degree of Maternal Fetal Interaction

Successful mammalian pregnancies are a result of complex physiological, endocrinological, and immunological processes that combine to create an environment where the mother is tolerant to the semi-allogeneic fetus. Our knowledge of the mechanisms that contribute to maternal tolerance is derived mainly from human and murine studies of haemochorial placentation. However, as this is the most invasive type of placentation it cannot be assumed that identical mechanisms apply to the less invasive epitheliochorial placentation found in other species such as ruminants. Here, we examine three features associated with reproductive immune regulation in a transformed ovine trophoblast cell line and ex-vivo ovine reproductive tissues collected at term, namely: major histocompatibility complex (MHC) expression, Indoleamine 2,3 dioxygenase-1 (IDO-1) expression, and Natural Killer (NK) cell infiltration. High levels of MHC class I protein expression were detected at the surface of the trophoblast cell line using a pan-MHC class I specific monoclonal antibody. The majority of MHC class I transcripts isolated from the cell line clustered with classical MHC alleles. Transcriptional analysis of placental tissues identified only classical MHC class I transcripts. We found no evidence of constitutive transcription of IDO-1 in either the trophoblast cell line or placental tissues. Ex-vivo tissues collected from the materno-fetal interface were negative for cells expressing NKp46/NCR1. Collectively, these observations suggest that the relatively non-invasive synepitheliochorial placentation found in sheep has a more limited requirement for local immunoregulation compared to the more invasive haemochorial placentation of primates and rodents.

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.

Sequence-Based Genotyping of Expressed Swine Leukocyte Antigen Class I Alleles by Next-Generation Sequencing Reveal Novel Swine Leukocyte Antigen Class I Haplotypes and Alleles in Belgian, Danish, and Kenyan Fattening Pigs and Göttingen Minipigs

The need for typing of the swine leukocyte antigen (SLA) is increasing with the expanded use of pigs as models for human diseases and organ-transplantation experiments, their use in infection studies, and for design of veterinary vaccines. Knowledge of SLA sequences is furthermore a prerequisite for the prediction of epitope binding in pigs. The low number of known SLA class I alleles and the limited knowledge of their prevalence in different pig breeds emphasizes the need for efficient SLA typing methods. This study utilizes an SLA class I-typing method based on next-generation sequencing of barcoded PCR amplicons. The amplicons were generated with universal primers and predicted to resolve 68-88% of all known SLA class I alleles dependent on amplicon size. We analyzed the SLA profiles of 72 pigs from four different pig populations; Göttingen minipigs and Belgian, Kenyan, and Danish fattening pigs. We identified 67 alleles, nine previously described haplotypes and 15 novel haplotypes. The highest variation in SLA class I profiles was observed in the Danish pigs and the lowest among the Göttingen minipig population, which also have the highest percentage of homozygote individuals. Highlighting the fact that there are still numerous unknown SLA class I alleles to be discovered, a total of 12 novel SLA class I alleles were identified. Overall, we present new information about known and novel alleles and haplotypes and their prevalence in the tested pig populations.

A comprehensive mapping of the structure and gene organisation in the sheep MHC class I region

Background

The major histocompatibility complex (MHC) is a chromosomal region that regulates immune responsiveness in vertebrates. This region is one of the most important for disease resistance because it has been associated with resistance or susceptibility to a wide variety of diseases and because the MHC often accounts for more of the variance than other loci. Selective breeding for disease resistance is becoming increasingly common in livestock industries, and it is important to determine how this will influence MHC polymorphism and resistance to diseases that are not targeted for selection. However, in sheep the order and sequence of the protein coding genes is controversial. Yet this information is needed to determine precisely how the MHC influences resistance and susceptibility to disease.

Methods

CHORI bacterial artificial chromosomes (BACs) known to contain sequences from the sheep MHC class I region were sub-cloned, and the clones partially sequenced. The resulting sequences were analysed and re-assembled to identify gene content and organisation within each BAC. The low resolution MHC class I physical map was then compared to the cattle reference genome, the Chinese Merino sheep MHC map published by Gao, et al. (2010) and the recently available sheep reference genome.

Results

Immune related class I genes are clustered into 3 blocks; beta, kappa and a novel block not previously identified in other organisms. The revised map is more similar to Bovidae maps than the previous sheep maps and also includes several genes previously not annotated in the Chinese Merino BAC assembly and others not currently annotated in the sheep reference chromosome 20. In particular, the organisation of nonclassical MHC class I genes is similar to that present in the cattle MHC. Sequence analysis and prediction of amino acid sequences of MHC class I classical and nonclassical genes was performed and it was observed that the map contained one classical and eight nonclassical genes together with three possible pseudogenes.

Conclusions

The comprehensive physical map of the sheep MHC class I region enhances our understanding of the genetic architecture of the class I MHC region in sheep and will facilitate future studies of MHC function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1992-4) contains supplementary material, which is available to authorized users.

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.

The functional significance of cattle major histocompatibility complex class I genetic diversity

Current concerns about food security highlight the importance of maintaining productive and disease-resistant livestock populations. Major histocompatibility complex (MHC) class I genes have a central role in immunity. A high level of diversity in these genes allows populations to survive despite exposure to rapidly evolving pathogens. This review aims to describe the key features of MHC class I genetic diversity in cattle and to discuss their role in disease resistance. Discussion centers on data derived from the cattle genome sequence and studies addressing MHC class I gene expression and function. The impact of intensive selection on MHC diversity is also considered. A high level of complexity in MHC class I genes and functionally related gene families is revealed. This highlights the need for increased efforts to determine key genetic components that govern cattle immune responses to disease, which is increasingly important in the face of changing human and environmental demands.

Natural selection on marine carnivores elaborated a diverse family of classical MHC class I genes exhibiting haplotypic gene content variation and allelic polymorphism

Pinnipeds, marine carnivores, diverged from terrestrial carnivores ~45 million years ago, before their adaptation to marine environments. This lifestyle change exposed pinnipeds to different microbiota and pathogens, with probable impact on their MHC class I genes. Investigating this question, genomic sequences were determined for 71 MHC class I variants: 27 from harbor seal and 44 from gray seal. These variants form three MHC class I gene lineages, one comprising a pseudogene. The second, a candidate nonclassical MHC class I gene, comprises a nonpolymorphic transcribed gene related to dog DLA-79 and giant panda Aime-1906. The third is the diversity lineage, which includes 62 of the 71 seal MHC class I variants. All are transcribed, and they minimally represent six harbor and 12 gray seal MHC class I genes. Besides species-specific differences in gene number, seal MHC class I haplotypes exhibit gene content variation and allelic polymorphism. Patterns of sequence variation, and of positions for positively selected sites, indicate the diversity lineage genes are the seals' classical MHC class I genes. Evidence that expansion of diversity lineage genes began before gray and harbor seals diverged is the presence in both species of two distinctive sublineages of diversity lineage genes. Pointing to further expansion following the divergence are the presence of species-specific genes and greater MHC class I diversity in gray seals than harbor seals. The elaboration of a complex variable family of classical MHC class I genes in pinnipeds contrasts with the single, highly polymorphic classical MHC class I gene of dog and giant panda, terrestrial carnivores.

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.

Difference in number of loci of swine leukocyte antigen classical class I genes among haplotypes

The structure of the entire genomic region of swine leukocyte antigen (SLA)-the porcine major histocompatibility complex--was recently elucidated in a particular haplotype named Hp-1.0 (H01). However, it has been suggested that there are differences in the number of loci of SLA genes, particularly classical class I genes, among haplotypes. To clarify the between-haplotype copy number variance in genes of the SLA region, we sequenced the genomic region carrying SLA classical class I genes on two different haplotypes, revealing increments of up to six in the number of classical class I genes in a single haplotype. All of the SLA-1(-like) (SLA-1 and newly designated SLA-12) and SLA-3 genes detected in the haplotypes thus analyzed were transcribed in the individual. The process by which duplication of SLA classical class I genes was likely to have occurred was interpreted from an analysis of repetitive sequences adjacent to the duplicated class I genes.