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

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.

Advances in understanding the genetic architecture of antibody response to paratuberculosis in sheep by heritability estimate and LDLA mapping analyses and investigation of candidate regions using sequence-based data

BACKGROUND

Paratuberculosis is a contagious and incurable disease that is caused by Mycobacterium avium subsp. paratuberculosis (MAP) with significant negative effects on animal welfare and farm profitability. Based on a large naturally infected flock over 12 years, we analyzed repeated enzyme-linked immunosorbent assay tests (ELISA), OvineSNP50 BeadChip genotypes and whole-genome sequences imputed from 56 influential animals. The main goals were to estimate the genetic parameters of proxy traits for resistance to MAP, identify genomic regions associated with the host's immune response against MAP and search for candidate genes and causative mutations through association and functional annotation analyses of polymorphisms identified by sequencing.

RESULTS

Two variables were derived from ELISA tests. The first, a binary variable, assessed the infection status of each animal over the entire productive life, while the second considered the level of antibody recorded over time. Very similar results were obtained for both variables. Heritability estimates of about 0.20 were found and a significant region capturing 18% and 13% of the genetic variance was detected on ovine chromosome 20 by linkage disequilibrium and linkage analysis on OvineSNP50 positions. Functional annotation and association analyses on the imputed sequence polymorphisms that were identified in this region were carried out. No significant variants showed a functional effect on the genes that mapped to this region, most of which belong to the major histocompatibility complex class II (MHC II). However, the conditional analysis led to the identification of two significant polymorphisms that can explain the genetic variance associated with the investigated genomic region.

CONCLUSIONS

Our results confirm the involvement of the host's genetics in susceptibility to MAP in sheep and suggest that selective breeding may be an option to limit the infection. The estimated heritability is moderate with a relevant portion being due to a highly significant region on ovine chromosome 20. The results of the combined use of sequence-based data and functional analyses suggest several genes belonging to the MHC II as the most likely candidates, although no mutations in their coding regions showed a significant association. Nevertheless, information from genotypes of two highly significant polymorphisms in the region can enhance the efficiency of selective breeding programs.

Population-wide genetic analysis of Ovar-DQA1 and DQA2 loci across sheep breeds in India revealed their evolutionary importance and fitness of sheep in a tropical climate

Genetic variability at the major histocompatibility complex (MHC) is important in any species due to significant role played by MHC for antigen presentation. DQA locus has not been studied for its genetic variability across sheep population in India. In the present study, MHC of sheep at DQA1 and DQA2 loci were evaluated across 17 Indian sheep breeds. Results revealed high degree of heterozygosity (10.34% to 100% for DQA1 and 37.39 to 100% for DQA2). 18 DQA1 alleles and 22 DQA2 alleles were isolated in different breeds. Nucleotide content for DQA region revealed richness of AT content (54.85% for DQA1 and 53.89% for DQA2). DQA1 and DQA2 sequences clustered independently. We could see evidence of divergence of DQA as DQA1 and DQA2 across sheep breeds. Wu-Kabat variability index revealed vast genetic variation across DQA1 and DQA2, specifically at peptide binding sites (PBS) that consisted 21 residues for DQA1 and 17 residues for DQA2. Evolutionary analysis revealed the presence of positive and balancing selection for DQA1 locus, however DQA2 was under purifying selection across sheep breeds. Higher heterozygosity and large diversity at both loci especially at PBS indicated the fitness of the sheep population for evading pathogens and adapt to the harsh tropical climate.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

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.

Identifications of immune-responsive genes for adaptative traits by comparative transcriptome analysis of spleen tissue from Kazakh and Suffolk sheep

Aridity and heat are significant environmental stressors that affect sheep adaptation and adaptability, thus influencing immunity, growth, reproduction, production performance, and profitability. The aim of this study was to profile mRNA expression levels in the spleen of indigenous Kazakh sheep breed for comparative analysis with the exotic Suffolk breed. Spleen histomorphology was observed in indigenous Kazakh sheep and exotic Suffolk sheep raised in Xinjiang China. Transcriptome sequencing of spleen tissue from the two breeds were performed via Illumina high-throughput sequencing technology and validated by RT-qPCR. Blood cytokine and IgG levels differed between the two breeds and IgG and IL-1β were significantly higher in Kazakh sheep than in Suffolk sheep (p < 0.05), though spleen tissue morphology was the same. A total of 52.04 Gb clean reads were obtained and the clean reads were assembled into 67,271 unigenes using bioinformatics analysis. Profiling analysis of differential gene expression showed that 1158 differentially expressed genes were found when comparing Suffolk with Kazakh sheep, including 246 up-regulated genes and 912 down-regulated genes. Utilizing gene ontology annotation and pathway analysis, 21 immune- responsive genes were identified as spleen-specific genes associated with adaptive traits and were significantly enriched in hematopoietic cell lineage, natural killer cell-mediated cytotoxicity, complement and coagulation cascades, and in the intestinal immune network for IgA production. Four pathways and up-regulated genes associated with immune responses in indigenous sheep played indispensable and promoting roles in arid and hot environments. Overall, this study provides valuable transcriptome data on the immunological mechanisms related to adaptive traits in indigenous and exotic sheep and offers a foundation for research into adaptive evolution.

A novel technique for retrospective genetic analysis of the response to vaccination or infection using cell-free DNA from archived sheep serum and plasma

Genetic variation is associated with differences in disease resistance and susceptibility among individuals within a population. To date, molecular genetic analyses of host responses have relied on extraction of genomic DNA from whole blood or tissue samples. However, such samples are not routinely collected during large-scale field studies. We demonstrate that cell-free genomic DNA (cfDNA) may be extracted and amplified from archived plasma samples, allowing retrospective analysis of host genetic diversity. This technique was also applicable to archived serum samples up to 35 years old and to different ruminant species. As proof of concept, we used this cfDNA approach to genotype the major histocompatibility complex (MHC) class II DRB1 locus of 224 Merino sheep which had participated in field trials of a commercial Haemonchus contortus vaccine, Barbervax^®, in Australia. This identified a total of 51 different DRB1 alleles and their relative frequencies. This is the first study to examine host MHC diversity using DNA extracted from archived plasma samples, an approach that may be applied to retrospective analyses of genetic diversity and responses to vaccination or infection across different species and populations.

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.

MHC haplotype diversity in Icelandic horses determined by polymorphic microsatellites

The Icelandic horse has been maintained as a closed population in its eponymous homeland for many generations, with no recorded introductions of new horses of any breed since the year 1000 CE. Here we determined the diversity of major histocompatibility complex (MHC) haplotypes in 156 Icelandic horses from two groups, based on a panel of 12 polymorphic intra-MHC microsatellites tested in families of various composition. We identified a total of 79 MHC haplotypes in these two groups, including one documented intra-MHC recombination event from a total of 147 observed meioses. None of these MHC haplotypes have been previously described in any other horse breed. Only one MHC homozygote was found in the entire population studied. These results indicate a very high level of MHC heterozygosity and haplotype diversity in the Icelandic horse. The environment in Iceland is remarkable for its lack of common agents of equine infectious disease, including equine herpesvirus type 1, influenza virus, and streptococcus equi. The driving forces for maintenance of MHC heterozygosity in Icelandic horses must thus be sought outside of these major horse pathogens. Based on our results, we propose that intra-MHC recombination may play a major role in the generation of novel haplotypes.

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.

Comparative MHC nomenclature: report from the ISAG/IUIS-VIC committee 2018

Significant progress has been made over the last decade in defining major histocompatibility complex (MHC) diversity at the nucleotide, allele, haplotype, diplotype, and population levels in many non-human species. Much of this progress has been driven by the increased availability and reduced costs associated with nucleotide sequencing technologies. This report provides an update on the activities of the comparative MHC nomenclature committee which is a standing committee of both the International Society for Animal Genetics (ISAG) and the International Union of Immunological Societies (IUIS) where it operates under the umbrella of the Veterinary Immunology Committee (VIC). A previous report from this committee in 2006 defined the role of the committee in providing guidance in the development of a standardized nomenclature for genes and alleles at MHC loci in non-human species. It described the establishment of the Immuno Polymorphism Database, IPD-MHC, which continues to provide public access to high quality MHC sequence data across a range of species. In this report, guidelines for the continued development of a universal MHC nomenclature framework are described, summarizing the continued development of each species section within the IPD-MHC project.