Generation and maintenance of diversity in the cattle MHC class I region

Fully Phased Population-Prevalent East African Cattle BoLA-I Alleles Determined Using PacBio HiFi Long-Read Sequencing Represent Five Novel Specificities With Distinctive Peptide Binding Potential

Due to factors such as lower biosecurity, greater wildlife/farm animal interfaces, and environmental challenges, cattle in sub-Saharan Africa are exposed to more diverse and intensive bacterial, viral and protozoan pathogen challenges than cattle in Europe and other high-income regions of the world. Classical class I genes of the major histocompatibility complex (MHC) contribute to protection from diseases caused by these pathogens by refining a huge pool of potential pathogen-derived peptide ligands into a smaller ensemble for presentation to CD8+ T cells. Knowledge of population-prevalent MHC alleles is therefore critical for evidence-based approaches to vaccine design and improved understanding of pathogen resistance. Whereas variation in MHC molecules is understood in most detail for European Bos taurus, the alleles expressed by Africa's cattle remain poorly defined. We have leveraged recent improvements in the accuracy of PacBio high-fidelity (HiFi) circular consensus sequencing (CCS) and adapted stringent sequence filtering algorithms to identify hundreds of as yet uncharacterised fully phased BoLA-I alleles from multiple populations of African taurine (Ankole) and indicine (Zebu) cattle in East Africa. The analysis highlights a convergence of population-prevalent class I MHC allelic repertoires in taurine and indicine cattle, likely due to the similar pathogen-driven selective pressures. Our analysis of the anchor residue accommodating pockets of these prevalent alleles revealed extremely high levels of polymorphism, which contrast with Holstein alleles that exhibit a more limited repertoire of MHC specificity-determining pocket residues, potentially constraining the breadth of peptide presentation. However, in the context of considerable sequence and physicochemical variation in the pocket-forming residues, it was possible to discern overlaps in the predicted peptide binding spectrum. Interrogation of potential differences in peptide binding specificities with European B. taurus alleles revealed that the fully phased African cattle class I MHC alleles represent five novel specificities. We envisage that this novel finding will find broad application in assessing potentially achievable vaccination coverages of future pathogen-encoded vaccine candidates against important intracellular pathogens. One aim of future research should be to leverage recent improvements in the sensitivity of mass spectrometry combined with immunoprecipitation of peptides bound to African cattle MHC to search directly for T-cell epitopes in the context of the inferred 'supertype' diversity.

Characterization and transcription of non-classical class I major histocompatibility complex (MHC) genes in buffaloes

The objective of the present study was to characterize non-classical class I major histocompatibility complex (MHC) genes in buffaloes and evaluate the expression of these genes in different tissue components of the placenta of buffaloes during pregnancy and in trophoblastic cells after stimulation using lipopolysaccharide (LPS). To do this, DNA was extracted from the blood of buffaloes and was subjected to PCR testing and sequencing of the genes NC3 and MICB. The RNA extracted from the placentome and intercotyledonary region of buffaloes in their first (n = 6), second (n = 6) and third (n = 6) trimesters of gestation was subjected to real-time PCR. Explants were created using the chorioallantoic membrane and two experimental groups were established: control and stimulated with LPS for four hours to evaluate the gene expression profile. Analysis on the sequences obtained showed that the genes NC3 and MICB of buffaloes were homologous with those of cattle, with high similarity in the analysis on the sequence variation pattern. The gene expression analysis showed that the genes assessed were transcribed at stages and in placental tissue that differed from what was seen in cattle. The transcription of these genes varied in the tissues studied, with greater transcription of MICB in the intercotyledonary region over the first third of gestation, while the genes studied in the placentome presented low rates of transcription. The trophoblastic cells of the chorioallantoic membrane stimulated with LPS for six hours did not present non-classic MFC-I transcription alterations. The present study therefore provides additional knowledge regarding the immune regulation of placental tissues of buffaloes.

The antigen recognition portion of African buffalo class I MHC is highly polymorphic, consistent with a complex pathogen challenge environment, and the 3' region suggests distinct haplotype configurations

African buffalo (Syncerus caffer) have been distinct from the Auroch lineage leading to domestic cattle for 5 million years, and are reservoirs of multiple pathogens, that affect introduced domestic cattle. To date, there has been no analysis of the class I MHC locus in African buffalo. We present the first data on African buffalo class I MHC, which demonstrates that gene and predicted protein coding sequences are approximately 86-87% similar to that of African domestic cattle in the peptide binding region. The study also shows concordance in the distribution of codons with elevated posterior probabilities of positive selection in the buffalo class I MHC and known antigen binding sites in cattle. Overall, the diversity in buffalo class I sequences appears greater than that in cattle, perhaps related to a more complex pathogen challenge environment in Africa. However, application of NetMHCpan suggested broad clustering of peptide binding specificities between buffalo and cattle. Furthermore, in the case of at least 20 alleles, critical peptide-binding residues appear to be conserved with those of cattle, including at secondary anchor residues. Alleles with six different length transmembrane regions were detected. This preliminary analysis suggests that like cattle, but unlike most other mammals, African buffalo appears to exhibit configuration (haplotype) variation in which the loci are expressed in distinct combinations.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.

High-resolution characterization of the structural features and genetic variation of six feline leukocyte antigen class I loci via single molecule, real-time (SMRT) sequencing

Of the 12 full-length feline leukocyte antigen class I (FLAI) loci, 3 are presumed to be classical: FLAI-E, FLAI-H, and FLAI-K. As diversity is a class Ia hallmark, multi-allelism is an important surrogate supporting a classical designation, in the absence of direct demonstration of T-cell restriction. Conversely, limited polymorphism at an expressed locus suggests regulation of immune effectors with invariant receptors, and non-classical status. FLAI-A, FLAI-J, FLAI-L, and FLAI-O are putative class Ib genes in cats. For both classes, identifying prevalent variants across outbred populations can illuminate specific genotypes to be prioritized for immune studies, as shared alleles direct shared responses. Since variation is concentrated in exons 2 and 3, which encode the antigen-binding domains, partial-length cloning/sequencing can be used for allele discovery, but is laborious and occasionally ambiguous. Here we develop a targeted approach to FLAI genotyping, using the single-molecule real-time (SMRT) platform, which allows full-length (3.4-kb) reads without assembly. Consensus sequences matched full-length Sanger references. Thirty-one new class Ia genes were found in 17 cats. Alleles segregated strongly by loci, and the origins of formerly difficult-to-assign sequences were resolved. Although not targeted, FLAI-L and FLAI-J, and the pseudogene FLAI-F, were also returned. Eighteen class Ib alleles were identified. Diversity was restricted and outside hypervariable regions. Both class Ib genes were transcriptionally active. Novel alternative splicing of FLAI-L was observed. SMRT sequencing of FLAI amplicons is useful for full-length genotyping at feline class Ia loci. High-throughput sequencing could allow highly accurate allele surveys in large cat cohorts.

Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep

Background

Because very large numbers of single nucleotide polymorphisms (SNPs) are now available throughout the genome, they are particularly suitable for the detection of genomic regions where a reduction in heterozygosity has occurred and they offer new opportunities to improve the accuracy of inbreeding (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F) estimates. Runs of homozygosity (ROH) are contiguous lengths of homozygous segments of the genome where the two haplotypes inherited from the parents are identical. Here, we investigated the occurrence and distribution of ROH using a medium-dense SNP panel to characterize autozygosity in 516 Valle del Belice sheep and to identify the genomic regions with high ROH frequencies.

Results

We identified 11,629 ROH and all individuals displayed at least one ROH longer than 1 Mb. The mean value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F$$\end{document}F estimated from ROH longer than1 Mb was 0.084 ± 0.061. ROH that were shorter than 10 Mb predominated. The highest and lowest coverages of Ovis aries chromosomes (OAR) by ROH were on OAR24 and OAR1, respectively. The number of ROH per chromosome length displayed a specific pattern, with higher values for the first three chromosomes. Both number of ROH and length of the genome covered by ROH varied considerably between animals. Two hundred and thirty-nine SNPs were considered as candidate markers that may be under directional selection and we identified 107 potential candidate genes. Six genomic regions located on six chromosomes, corresponding to ROH islands, are presented as hotspots of autozygosity, which frequently coincided with regions of medium recombination rate. According to the KEGG database, most of these genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. A genome scan revealed the presence of ROH islands in genomic regions that harbor candidate genes for selection in response to environmental stress and which underlie local adaptation.

Conclusions

These results suggest that natural selection has, at least partially, a role in shaping the genome of Valle del Belice sheep and that ROH in the ovine genome may help to detect genomic regions involved in the determinism of traits under selection.

Electronic supplementary material

The online version of this article (10.1186/s12711-017-0360-z) contains supplementary material, which is available to authorized users.

Improved Prediction of Bovine Leucocyte Antigens (BoLA) Presented Ligands by Use of Mass-Spectrometry-Determined Ligand and in Vitro Binding Data

Peptide binding to MHC class I molecules is the single most selective step in antigen presentation and the strongest single correlate to peptide cellular immunogenicity. The cost of experimentally characterizing the rules of peptide presentation for a given MHC-I molecule is extensive, and predictors of peptide–MHC interactions constitute an attractive alternative. Recently, an increasing amount of MHC presented peptides identified by mass spectrometry (MS ligands) has been published. Handling and interpretation of MS ligand data is, in general, challenging due to the polyspecificity nature of the data. We here outline a general pipeline for dealing with this challenge and accurately annotate ligands to the relevant MHC-I molecule they were eluted from by use of GibbsClustering and binding motif information inferred from in silico models. We illustrate the approach here in the context of MHC-I molecules (BoLA) of cattle. Next, we demonstrate how such annotated BoLA MS ligand data can readily be integrated with in vitro binding affinity data in a prediction model with very high and unprecedented performance for identification of BoLA-I restricted T-cell epitopes. The prediction model is freely available at http://www.cbs.dtu.dk/services/NetMHCpan/NetBoLApan. The approach has here been applied to the BoLA-I system, but the pipeline is readily applicable to MHC systems in other species.

Rapid identification of bovine MHCI haplotypes in genetically divergent cattle populations using next-generation sequencing

The major histocompatibility complex (MHC) region contains many genes that are key regulators of both innate and adaptive immunity including the polymorphic MHCI and MHCII genes. Consequently, the characterisation of the repertoire of MHC genes is critical to understanding the variation that determines the nature of immune responses. Our current knowledge of the bovine MHCI repertoire is limited with only the Holstein-Friesian breed having been studied in any depth. Traditional methods of MHCI genotyping are of low resolution and laborious and this has been a major impediment to a more comprehensive analysis of the MHCI repertoire of other cattle breeds. Next-generation sequencing (NGS) technologies have been used to enable high throughput and much higher resolution MHCI typing in a number of species. In this study we have developed a MiSeq platform approach and requisite bioinformatics pipeline to facilitate typing of bovine MHCI repertoires. The method was validated initially on a cohort of Holstein-Friesian animals and then demonstrated to enable characterisation of MHCI repertoires in African cattle breeds, for which there was limited or no available data. During the course of these studies we identified >140 novel classical MHCI genes and defined 62 novel MHCI haplotypes, dramatically expanding the known bovine MHCI repertoire.

Expression of bovine non-classical major histocompatibility complex class I proteins in mouse P815 and human K562 cells

Major histocompatibility complex class I (MHC-I) proteins can be expressed as cell surface or secreted proteins. To investigate whether bovine non-classical MHC-I proteins are expressed as cell surface or secreted proteins, and to assess the reactivity pattern of monoclonal antibodies with non-classical MHC-I isoforms, we expressed the MHC proteins in murine P815 and human K562 (MHC-I deficient) cells. Following antibiotic selection, stably transfected cell lines were stained with H1A or W6/32 antibodies to detect expression of the MHC-I proteins by flow cytometry. Two non-classical proteins (BoLA-NC1*00501 and BoLA-NC3*00101) were expressed on the cell surface in both cell lines. Surprisingly, the BoLA-NC4*00201 protein was expressed on the cell membrane of human K562 but not mouse P815 cells. Two non-classical proteins (BoLA-NC1*00401, which lacks a transmembrane domain, and BoLA-NC2*00102) did not exhibit cell surface expression. Nevertheless, Western blot analyses demonstrated expression of the MHC-I heavy chain in all transfected cell lines. Ammonium-sulfate precipitation of proteins from culture supernatants showed that BoLA-NC1*00401 was secreted and that all surface expressed proteins where shed from the cell membrane by the transfected cells. Interestingly, the surface expressed MHC-I proteins were present in culture supernatants at a much higher concentration than BoLA-NC1*00401. This comprehensive study shows that bovine non-classical MHC-I proteins BoLA-NC1*00501, BoLA-NC3*00101, and BoLA-NC4*00201 are expressed as surface isoforms with the latter reaching the cell membrane only in K562 cells. Furthermore, it demonstrated that BoLA-NC1*00401 is a secreted isoform and that significant quantities of membrane associated MHC-I proteins can be shed from the cell membrane.

Sequence diversity between class I MHC loci of African native and introduced Bos taurus cattle in Theileria parva endemic regions: in silico peptide binding prediction identifies distinct functional clusters

There is strong evidence that the immunity induced by live vaccination for control of the protozoan parasite Theileria parva is mediated by class I MHC-restricted CD8(+) T cells directed against the schizont stage of the parasite that infects bovine lymphocytes. The functional competency of class I MHC genes is dependent on the presence of codons specifying certain critical amino acid residues that line the peptide binding groove. Compared with European Bos taurus in which class I MHC allelic polymorphisms have been examined extensively, published data on class I MHC transcripts in African taurines in T. parva endemic areas is very limited. We utilized the multiplexing capabilities of 454 pyrosequencing to make an initial assessment of class I MHC allelic diversity in a population of Ankole cattle. We also typed a population of exotic Holstein cattle from an African ranch for class I MHC and investigated the extent, if any, that their peptide-binding motifs overlapped with those of Ankole cattle. We report the identification of 18 novel allelic sequences in Ankole cattle and provide evidence of positive selection for sequence diversity, including in residues that predominantly interact with peptides. In silico functional analysis resulted in peptide binding specificities that were largely distinct between the two breeds. We also demonstrate that CD8(+) T cells derived from Ankole cattle that are seropositive for T. parva do not recognize vaccine candidate antigens originally identified in Holstein and Boran (Bos indicus) cattle breeds.

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.

Pathogenicity of Bovine Neonatal Pancytopenia-associated vaccine-induced alloantibodies correlates with Major Histocompatibility Complex class I expression

Bovine Neonatal Pancytopenia (BNP), a fatal bleeding syndrome of neonatal calves, is caused by maternal alloantibodies absorbed from colostrum and is characterized by lymphocytopenia, thrombocytopenia and bone marrow hypoplasia. An inactivated viral vaccine is the likely source of alloantigens inducing BNP-associated alloantibodies in the dam. In this study the specificity of BNP alloantibodies was assessed and was linked to the pathology of BNP. We demonstrated that Major Histocompatibility Complex class I (MHC I) and Very Late Antigen-3, an integrin α3/β1 heterodimer, were the major targets of BNP alloantibodies. However, alloantibody binding to various bovine cell types correlated with MHC I expression, rather than integrin β1 or α3 expression. Likewise, alloantibody-dependent complement-mediated cell lysis correlated strongly with MHC I expression. Examination of several tissues of third trimester bovine foetuses revealed that cells, shown to be affected in calves with BNP, were characterized by high MHC class I expression and high levels of alloantibody binding. We conclude that in spite of the heterogeneous specificity of BNP associated maternal alloantibodies, MHC I-specific antibodies mediate the pathogenicity of BNP in the calf and that cells with high MHC I expression were preferentially affected in BNP.

The assembly and characterisation of two structurally distinct cattle MHC class I haplotypes point to the mechanisms driving diversity

In cattle, there are six classical MHC class I genes that are variably present between different haplotypes. Almost all known haplotypes contain between one and three genes, with an allele of Gene 2 present on the vast majority. However, very little is known about the sequence and therefore structure and evolutionary history of this genomic region. To address this, we have refined the MHC class I region in the Hereford cattle genome assembly and sequenced a complete A14 haplotype from a homozygous Holstein. Comparison of the two haplotypes revealed extensive variation within the MHC class Ia region, but not within the flanking regions, with each gene contained within a conserved 63- to 68-kb sequence block. This variable region appears to have undergone block gene duplication and likely deletion at regular breakpoints, suggestive of a site-specific mechanism. Phylogenetic analysis using complete gene sequences provided evidence of allelic diversification via gene conversion, with breakpoints between each of the extracellular domains that were associated with high guanine-cytosine (GC) content. Advancing our knowledge of cattle MHC class I evolution will help inform investigations of cattle genetic diversity and disease resistance.

Transcription of non-classic major histocompatibility complex (MHC) class I in the bovine placenta throughout gestation and after Brucella abortus infection

Transcription of non-classical major histocompatibility complex class I (MHC-I) was assessed in the bovine placenta throughout gestation. Additionally, the effect of Brucella abortus infection on expression of non-classical MHC-I was also evaluated using a chorioallantoic membrane explant model of infection. The non-classical MHC-I genes MICB and NC3 had higher levels of transcription in the intercotyledonary region when compared to the placentome, which had higher levels of transcription at the second trimester of gestation. NC1 and classical MHC-I had very low levels of transcription throughout gestation. Trophoblastic cells of B. abortus-infected chorioallantoic membrane explants had an increase in transcription of non-classical MHC-I at 4h post infection. Therefore, this study provides an analysis of non-classical MHC-I transcription at different stages of gestation and different placental tissues, and during B. abortus infection. These findings provide additional knowledge on immune regulation in placental tissues, a known immune-privileged site.

Haematopoietic depletion in vaccine-induced neonatal pancytopenia depends on both the titre and specificity of alloantibody and levels of MHC I expression

Bovine Neonatal Pancytopenia (BNP) is a disease of calves characterised by haematopoietic depletion, mediated by ingestion of alloantibodies in colostrum. It has been linked epidemiologically to vaccination of the dams of affected calves with a particular vaccine (Pregsure) containing a novel adjuvant. Evidence suggests that BNP-alloantibodies are directed against MHC I molecules, induced by contaminant bovine cellular material from Madin-Darby Bovine Kidney (MDBK) cells used in the vaccine's production. We aimed to investigate the specificity of BNP-alloantibody for bovine MHC I alleles, particularly those expressed by MDBK cells, and whether depletion of particular cell types is due to differential MHC I expression levels. A complement-mediated cytotoxicity assay was used to assess functional serum alloantibody titres in BNP-dams, Pregsure-vaccinated dams with healthy calves, cows vaccinated with an alternative product and unvaccinated controls. Alloantibody specificity was investigated using transfected mouse lines expressing the individual MHC I alleles identified from MDBK cells and MHC I-defined bovine leukocyte lines. All BNP-dams and 50% of Pregsure-vaccinated cows were shown to have MDBK-MHC I specific alloantibodies, which cross-reacted to varying degrees with other MHC I genotypes. MHC I expression levels on different blood cell types, assessed by flow cytometry, were found to correlate with levels of alloantibody-mediated damage in vitro and in vivo. Alloantibody-killed bone marrow cells were shown to express higher levels of MHC I than undamaged cells. The results provide evidence that MHC I-specific alloantibodies play a dominant role in the pathogenesis of BNP. Haematopoietic depletion was shown to be dependent on the titre and specificity of alloantibody produced by individual cows and the density of surface MHC I expression by different cell types. Collectively, the results support the hypothesis that MHC I molecules originating from MDBK cells used in vaccine production, coupled with a powerful adjuvant, are responsible for the generation of pathogenic alloantibodies.

Characterisation of major histocompatibility complex class I genes at the fetal-maternal interface of marsupials

Major histocompatibility complex class I molecules (MHC-I) are expressed at the cell surface and are responsible for the presentation of self and non-self antigen repertoires to the immune system. Eutherian mammals express both classical and non-classical MHC-I molecules in the placenta, the latter of which are thought to modulate the maternal immune response during pregnancy. Marsupials last shared a common ancestor with eutherian mammals such as humans and mice over 160 million years ago. Since, like eutherians, they have an intra-uterine development dependent on a placenta, albeit a short-lived and less invasive one, they provide an opportunity to investigate the evolution of MHC-I expression at the fetal-maternal interface. We have characterised MHC-I mRNA expression in reproductive tissues of the tammar wallaby (Macropus eugenii) from the time of placental attachment to day 25 of the 26.5 day pregnancy. Putative classical MHC-I genes were expressed in the choriovitelline placenta, fetus, and gravid endometrium throughout the whole of this period. The MHC-I classical sequences were phylogenetically most similar to the Maeu-UC (50/100 clones) and Maeu-UA genes (7/100 clones). Expression of three non-classical MHC-I genes (Maeu-UD, Maeu-UK and Maeu-UM) were also present in placental samples. The results suggest that expression of classical and non-classical MHC-I genes in extant marsupial and eutherian mammals may have been necessary for the evolution of the ancestral therian placenta and survival of the mammalian fetus at the maternal-fetal interface.

Bovine Neonatal Pancytopenia is a heritable trait of the dam rather than the calf and correlates with the magnitude of vaccine induced maternal alloantibodies not the MHC haplotype

Bovine Neonatal Pancytopenia (BNP), a bleeding syndrome of neonatal calves, is caused by alloantibodies absorbed from the colostrum of particular cows. A commercial BVD vaccine is the likely source of alloantigens eliciting BNP associated alloantibodies. We hypothesized that the rare occurrence of BNP in calves born to vaccinated dams could be associated with genetic differences within dams and calves. We found that the development of BNP within calves was a heritable trait for dams, not for calves and had a high heritability of 19%. To elucidate which genes play a role in the development of BNP we sequenced candidate genes and characterized BNP alloantibodies. Alloantigens present in the vaccine have to be presented to the dam’s immune system via MHC class II, however sequencing of DRB3 showed no differences in MHC class II haplotype between BNP and non-BNP dams. MHC class I, a highly polymorphic alloantigen, is an important target of BNP alloantibodies. Using a novel sequence based MHC class I typing method, we found no association of BNP with MHC class I haplotype distribution in dams or calves. Alloantibodies were detected in both vaccinated BNP and non-BNP dams and we found no differences in alloantibody characteristics between these groups, but alloantibody levels were significantly higher in BNP dams. We concluded that the development of BNP in calves is a heritable trait of the dam rather than the calf and genetic differences between BNP and non-BNP dams are likely due to genes controlling the quantitative alloantibody response following vaccination.

Characterization of binding specificities of bovine leucocyte class I molecules: impacts for rational epitope discovery

The binding of peptides to classical major histocompatibility complex (MHC) class I proteins is the single most selective step in antigen presentation. However, the peptide-binding specificity of cattle MHC (bovine leucocyte antigen, BoLA) class I (BoLA-I) molecules remains poorly characterized. Here, we demonstrate how a combination of high-throughput assays using positional scanning combinatorial peptide libraries, peptide dissociation, and peptide-binding affinity binding measurements can be combined with bioinformatics to effectively characterize the functionality of BoLA-I molecules. Using this strategy, we characterized eight BoLA-I molecules, and found the peptide specificity to resemble that of human MHC-I molecules with primary anchors most often at P2 and P9, and occasional auxiliary P1/P3/P5/P6 anchors. We analyzed nine reported CTL epitopes from Theileria parva, and in eight cases, stable and high affinity binding was confirmed. A set of peptides were tested for binding affinity to the eight BoLA proteins and used to refine the predictors of peptide-MHC binding NetMHC and NetMHCpan. The inclusion of BoLA-specific peptide-binding data led to a significant improvement in prediction accuracy for reported T. parva CTL epitopes. For reported CTL epitopes with weak or no predicted binding, these refined prediction methods suggested presence of nested minimal epitopes with high-predicted binding affinity. The enhanced affinity of the alternative peptides was in all cases confirmed experimentally. This study demonstrates how biochemical high-throughput assays combined with immunoinformatics can be used to characterize the peptide-binding motifs of BoLA-I molecules, boosting performance of MHC peptide-binding prediction methods, and empowering rational epitope discovery in cattle.

Bovine neonatal pancytopenia (BNP): novel insights into the incidence, vaccination-associated epidemiological factors and a potential genetic predisposition for clinical and subclinical cases

Bovine neonatal pancytopenia (BNP) is a haemorrhagic disease of newborn calves elicited by colostrum from specific cows. Two studies have indicated that BNP-inducing colostrum might be associated with alloantibodies directed against MHC class I in response to vaccination with a distinct inactivated viral vaccine. However, the proportion of alloantibody-producing individuals by far exceeds the proportion of clinical BNP cases in the vaccinated population. This raises the question about the incidence of subclinical, unrecognised cases and also suggests further factors involved in BNP pathogenesis, e.g., genetic predisposition. Our results on neonatal calves from a closely monitored resource population confirmed the hypothesis of a genetic predisposition for clinical BNP and suggest that the predisposition is also involved in subclinical BNP-cases. No indication was obtained for a higher frequency of subclinical BNP-cases compared with clinical cases. Neither time point nor frequency of vaccination was a relevant factor for BNP in our resource population.

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.

Polymorphisms and tissue expression of the feline leukocyte antigen class I loci FLAI-E, FLAI-H, and FLAI-K

Cytotoxic CD8+ T-cell immunosurveillance for intracellular pathogens, such as viruses, is controlled by classical major histocompatibility complex (MHC) class Ia molecules, and ideally, these antiviral T-cell populations are defined by the specific peptide and restricting MHC allele. Surprisingly, despite the utility of the cat in modeling human viral immunity, little is known about the feline leukocyte antigen class I complex (FLAI). Only a few coding sequences with uncertain locus origin and expression patterns have been reported. Of 19 class I genes, three loci--FLAI-E, FLAI-H, and FLAI-K--are predicted to encode classical molecules, and our objective was to evaluate their status by analyzing polymorphisms and tissue expression. Using locus-specific, PCR-based genotyping, we amplified 33 FLAI-E, FLAI-H, and FLAI-K alleles from 12 cats of various breeds, identifying, for the first time, alleles across three distinct loci in a feline species. Alleles shared the expected polymorphic and invariant sites in the α1/α2 domains, and full-length cDNA clones possessed all characteristic class Ia exons. Alleles could be assigned to a specific locus with reasonable confidence, although there was evidence of potentially confounding interlocus recombination between FLAI-E and FLAI-K. Only FLAI-E, FLAI-H, and FLAI-K origin alleles were amplified from cDNAs of multiple tissue types. We also defined hypervariable regions across these genes, which permitted the assignment of names to both novel and established alleles. As predicted, FLAI-E, FLAI-H, and FLAI-K fulfill the major criteria of class Ia genes. These data represent a necessary prerequisite for studying epitope-specific antiviral CD8+ T-cell responses in cats.

Cattle MHC nomenclature: is it possible to assign sequences to discrete class I genes?

The cattle major histocompatibility complex (MHC) region contains a variable number of classical class I genes encoding polymorphic molecules involved in antigen presentation. Six classical class I genes have been described, but assigning sequences to these genes has proved problematic. We propose a refinement of the existing nomenclature, which currently names the 97 known classical class I sequences in a single series. Phylogenetic analysis of the 3' portion of the coding region allows segregation of these into six groups; thus, we have prefixed existing names with the appropriate number. Although it is clear that some of these groups correspond to discrete genes, it is currently not possible to state definitively that all do. However, the main groupings are consistent, and in conjunction with other evidence, we feel it is now appropriate to rename the sequences accordingly. Segregation of sequences into groups in this way will facilitate ongoing research and future use of the cattle MHC section of the Immuno Polymorphism Database.

Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA

BACKGROUND

Artificial selection has caused rapid evolution in domesticated species. The identification of selection footprints across domesticated genomes can contribute to uncover the genetic basis of phenotypic diversity.

METHODOLOGY/MAIN FINDINGS

Genome wide footprints of pig domestication and selection were identified using massive parallel sequencing of pooled reduced representation libraries (RRL) representing ∼2% of the genome from wild boar and four domestic pig breeds (Large White, Landrace, Duroc and Pietrain) which have been under strong selection for muscle development, growth, behavior and coat color. Using specifically developed statistical methods that account for DNA pooling, low mean sequencing depth, and sequencing errors, we provide genome-wide estimates of nucleotide diversity and genetic differentiation in pig. Widespread signals suggestive of positive and balancing selection were found and the strongest signals were observed in Pietrain, one of the breeds most intensively selected for muscle development. Most signals were population-specific but affected genomic regions which harbored genes for common biological categories including coat color, brain development, muscle development, growth, metabolism, olfaction and immunity. Genetic differentiation in regions harboring genes related to muscle development and growth was higher between breeds than between a given breed and the wild boar.

CONCLUSIONS/SIGNIFICANCE

These results, suggest that although domesticated breeds have experienced similar selective pressures, selection has acted upon different genes. This might reflect the multiple domestication events of European breeds or could be the result of subsequent introgression of Asian alleles. Overall, it was estimated that approximately 7% of the porcine genome has been affected by selection events. This study illustrates that the massive parallel sequencing of genomic pools is a cost-effective approach to identify footprints of selection.

Extensive polymorphism and evidence of immune selection in a highly dominant antigen recognized by bovine CD8 T cells specific for Theileria annulata

Although parasite strain-restricted CD8 T cell responses have been described for several protozoa, the precise role of antigenic variability in immunity is poorly understood. The tick-borne protozoan parasite Theileria annulata infects leukocytes and causes an acute, often fatal lymphoproliferative disease in cattle. Building on previous evidence of strain-restricted CD8 T cell responses to T. annulata, this study set out to identify and characterize the variability of the target antigens. Three antigens were identified by screening expressed parasite cDNAs with specific CD8 T cell lines. In cattle expressing the A10 class I major histocompatibility complex haplotype, A10-restricted CD8 T cell responses were shown to be focused entirely on a single dominant epitope in one of these antigens (Ta9). Sequencing of the Ta9 gene from field isolates of T. annulata demonstrated extensive sequence divergence, resulting in amino acid polymorphism within the A10-restricted epitope and a second A14-restricted epitope. Statistical analysis of the allelic sequences revealed evidence of positive selection for amino acid substitutions within the region encoding the CD8 T cell epitopes. Sequence differences in the A10-restricted epitope were shown to result in differential recognition by individual CD8 T cell clones, while clones also differed in their ability to recognize different alleles. Moreover, the representation of these clonal specificities within the responding CD8 T cell populations differed between animals. As well as providing an explanation for incomplete protection observed after heterologous parasite challenge of vaccinated cattle, these results have important implications for the choice of antigens for the development of novel subunit vaccines.

A first comparative map of copy number variations in the sheep genome

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.

MHC class I bound to an immunodominant Theileria parva epitope demonstrates unconventional presentation to T cell receptors

T cell receptor (TCR) recognition of peptide-MHC class I (pMHC) complexes is a crucial event in the adaptive immune response to pathogens. Peptide epitopes often display a strong dominance hierarchy, resulting in focusing of the response on a limited number of the most dominant epitopes. Such T cell responses may be additionally restricted by particular MHC alleles in preference to others. We have studied this poorly understood phenomenon using Theileria parva, a protozoan parasite that causes an often fatal lymphoproliferative disease in cattle. Despite its antigenic complexity, CD8+ T cell responses induced by infection with the parasite show profound immunodominance, as exemplified by the Tp1(214-224) epitope presented by the common and functionally important MHC class I allele N*01301. We present a high-resolution crystal structure of this pMHC complex, demonstrating that the peptide is presented in a distinctive raised conformation. Functional studies using CD8+ T cell clones show that this impacts significantly on TCR recognition. The unconventional structure is generated by a hydrophobic ridge within the MHC peptide binding groove, found in a set of cattle MHC alleles. Extremely rare in all other species, this feature is seen in a small group of mouse MHC class I molecules. The data generated in this analysis contribute to our understanding of the structural basis for T cell-dependent immune responses, providing insight into what determines a highly immunogenic p-MHC complex, and hence can be of value in prediction of antigenic epitopes and vaccine design.

Induction of a cross-reactive CD8(+) T cell response following foot-and-mouth disease virus vaccination

Foot-and-mouth disease virus (FMDV) causes a highly contagious infection in cloven-hoofed animals. Current inactivated FMDV vaccines generate short-term, serotype-specific protection, mainly through neutralizing antibody. An improved understanding of the mechanisms of protective immunity would aid design of more effective vaccines. We have previously reported the presence of virus-specific CD8(+) T cells in FMDV-vaccinated and -infected cattle. In the current study, we aimed to identify CD8(+) T cell epitopes in FMDV recognized by cattle vaccinated with inactivated FMDV serotype O. Analysis of gamma interferon (IFN-γ)-producing CD8(+) T cells responding to stimulation with FMDV-derived peptides revealed one putative CD8(+) T cell epitope present within the structural protein P1D, comprising residues 795 to 803 of FMDV serotype O UKG/2001. The restricting major histocompatibility complex (MHC) class I allele was N*02201, expressed by the A31 haplotype. This epitope induced IFN-γ release, proliferation, and target cell killing by αβ CD8(+) T cells, but not CD4(+) T cells. A protein alignment of representative samples from each of the 7 FMDV serotypes showed that the putative epitope is highly conserved. CD8(+) T cells from FMDV serotype O-vaccinated A31(+) cattle recognized antigen-presenting cells (APCs) loaded with peptides derived from all 7 FMDV serotypes, suggesting that CD8(+) T cells recognizing the defined epitope are cross-reactive to equivalent peptides derived from all of the other FMDV serotypes.

MHC gene copy number variation in Tasmanian devils: implications for the spread of a contagious cancer

Tasmanian devils face extinction owing to the emergence of a contagious cancer. Devil facial tumour disease (DFTD) is a clonal cancer spread owing to a lack of major histocompatibility complex (MHC) barriers in Tasmanian devil populations. We present a comprehensive screen of MHC diversity in devils and identify 25 MHC types and 53 novel sequences, but conclude that overall levels of MHC diversity at the sequence level are low. The majority of MHC Class I variation can be explained by allelic copy number variation with two to seven sequence variants identified per individual. MHC sequences are divided into two distinct groups based on sequence similarity. DFTD cells and most devils have sequences from both groups. Twenty per cent of individuals have a restricted MHC repertoire and contain only group I or only group II sequences. Counterintuitively, we postulate that the immune system of individuals with a restricted MHC repertoire may recognize foreign MHC antigens on the surface of the DFTD cell. The implication of these results for management of DFTD and this endangered species are discussed.

Evolutionary analysis of two classical MHC class I loci of the medaka fish, Oryzias latipes: haplotype-specific genomic diversity, locus-specific polymorphisms, and interlocus homogenization

The major histocompatibility complex (MHC) region of the teleost medaka (Oryzias latipes) contains two classical class I loci, UAA and UBA, whereas most lower vertebrates possess or express a single locus. To elucidate the allelic diversification and evolutionary relationships of these loci, we compared the BAC-based complete genomic sequences of the MHC class I region of three medaka strains and the PCR-based cDNA sequences of two more strains and two wild individuals, representing nine haplotypes. These were derived from two geographically distinct medaka populations isolated for four to five million years. Comparison of the genomic sequences showed a marked diversity in the region encompassing UAA and UBA even between the strains derived from the same population, and also showed an ancient divergence of these loci. cDNA analysis indicated that the peptide-binding domains of both UAA and UBA are highly polymorphic and that most of the polymorphisms were established in a locus-specific manner before the divergence of the two populations. Interallelic recombination between exons 2 and 3 encoding these domains was observed. The second intron of the UAA genes contains a highly conserved region with a palindromic sequence, suggesting that this region contributed to the recombination events. In contrast, the alpha3 domain is extremely homogenized not only within each locus but also between UAA and UBA regardless of populations. Two lineages of the transmembrane and cytoplasmic regions are also shared by UAA and UBA, suggesting that these two loci evolved with intimate genetic interaction through gene conversion or unequal crossing over.

A snapshot of the Mamu-B genes and their allelic repertoire in rhesus macaques of Chinese origin

The major histocompatibility complex class I gene repertoire was investigated in a large panel of rhesus macaques of Chinese origin. As observed in Indian animals, subjects of Chinese derivation display Mamu-B gene copy number variation, and the sum of expressed genes varies among haplotypes. In addition, these genes display differential transcription levels. The majority of the Mamu-B alleles discovered during this investigation appear to be unique for the population studied. Only one particular Mamu-B haplotype is shared between Indian and Chinese animals, and it must have been present in the progenitor stock. Hence, the data highlight the fact that most allelic polymorphism, and most of the Mamu-B haplotypes themselves, are of relatively recent origin and were most likely generated after the separation of the Indian and Chinese rhesus macaque populations.

Genomic location and characterisation of MIC genes in cattle

Major histocompatibility complex (MHC) class I chain-related (MIC) genes have been previously identified and characterised in human. They encode polymorphic class I-like molecules that are stress-inducible, and constitute one of the ligands of the activating natural killer cell receptor NKG2D. We have identified three MIC genes within the cattle genome, located close to three non-classical MHC class I genes. The genomic position relative to other genes is very similar to the arrangement reported in the pig MHC region. Analysis of MIC cDNA sequences derived from a range of cattle cell lines suggest there may be four MIC genes in total. We have investigated the presence of the genes in distinct and well-defined MHC haplotypes, and show that one gene is consistently present, while configuration of the other three genes appears variable.

Comparison of predicted binders in Rhipicephalus (Boophilus) microplus intestine protein variants BM86 Campo Grande strain, BM86 and BM95

This paper reports the sequence analysis of Bm86 Campo Grande strain comparing it with Bm86 and Bm95 antigens from the preparations TickGardPLUS and GavacTM, respectively. The PCR product was cloned into pMOSBlue and sequenced. The secondary structure prediction tool PSIPRED was used to calculate alpha helices and beta strand contents of the predicted polypeptide. The hydrophobicity profile was calculated using the algorithms from the Hopp and Woods method, in addition to identification of potential MHC class-I binding regions in the antigens. Pair-wise alignment revealed that the similarity between Bm86 Campo Grande strain and Bm86 is 0.2% higher than that between Bm86 Campo Grande strain and Bm95 antigens. The identities were 96.5% and 96.3% respectively. Major suggestive differences in hydrophobicity were predicted among the sequences in two specific regions.

An MHC-restricted CD8+ T-cell response is induced in cattle by foot-and-mouth disease virus (FMDV) infection and also following vaccination with inactivated FMDV

Foot-and-mouth disease virus (FMDV) causes a highly contagious disease of cloven-hooved animals that carries enormous economic consequences. CD8(+) cytotoxic T lymphocytes play an important role in protection and disease outcome in viral infections but, to date, the role of the CD8(+) T-cell immune response to FMDV remains unclear. This study aimed to investigate major histocompatibility complex (MHC) class I-restricted CD8(+) T-cell responses to FMDV in vaccinated and in infected cattle. An in vitro assay was used to detect antigen-specific gamma interferon release by CD8(+) T cells in FMDV-infected cattle of known MHC class I genotypes. A significant MHC class I-restricted CD8(+) T-cell response was detected to both FMDV strain O1 BFS and a recombinant fowlpox virus expressing the structural proteins of FMDV. Antigen-specific MHC class I-restricted CD8(+) T-cell responses were also detected in cattle vaccinated with inactivated FMDV. These responses were shown to be directed, at least in part, to epitopes within the structural proteins (P12A region) of the virus. By using mouse cells expressing single cattle MHC class I alleles, it was possible to identify the restriction elements in each case. Identification of these epitopes will facilitate the quantitative and qualitative analysis of FMDV-specific memory CD8(+) T cells in cattle and help to ensure that potential vaccines induce a qualitatively appropriate CD8(+) T-cell response.

Placental expression of major histocompatibility complex class I in bovine somatic clones

Abnormally increased placental expression of major histocompatibility complex class I (MHC-I) molecules at the trophoblastic surface has been suggested previously to be the cause of early fetal loss in nuclear transfer (NT) bovine pregnancies. Here, we report the lack of expression of MHC-I at the trophoblastic surface at D30 and D60 and in placentomes from D60 to term in placentas obtained by NT from three different genotypes and by artificial insemination, whatever the outcome of the pregnancy. MHC-I expression was assessed by immunohistochemistry using four different antibodies, including a novel beta2-microglobulin antibody. The MHC-I type of the clones was established using reference strand-mediated conformation analysis (RSCA); however, since it proved problematic to type the recipient animals in the same way, outcome of pregnancy could not be related to MHC compatibility. In conclusion, the present study provides no evidence to support abnormal expression of MHC-I on the trophoblastic surface in clones as a major cause of fetal loss during pregnancy after NT.

Complexity in cattle KIR genes: transcription and genome analysis

Killer immunoglobulin (Ig)-like receptors (KIRs) are the major functional natural killer (NK) cell receptors in human. The presence of KIR genes has only recently been demonstrated in other (non-primate) species, and their expression, genomic arrangement, and function in these species have yet to be investigated. In this study, we describe the KIR gene family in cattle. KIR sequences were amplified from cDNA derived from four animals. Seventeen new sequences were identified in total. Some are alleles of two previously described genes, and the remainder are representative of at least four additional genes. These cDNA data, together with analysis of the cattle genome sequence, confirm that, as in humans, cattle have multiple inhibitory and activating KIR genes, with variable haplotype composition, and putative framework genes. In contrast to human, the majority of the cattle KIR genes encode three Ig-domain KIRs; most of the inhibitory genes encode only one immunoreceptor tyrosine-based inhibitory motif (ITIM), and the activating genes encode molecules with arginine rather than the more usual lysine in the transmembrane domain. A divergent gene, 2DL1, encodes a two Ig-domain KIR with an unusual D0-D2 structure, and a distinct signaling domain with two ITIMs. Similarity to pig and human two Ig-domain (D0-D2) KIRs suggest these may be more related to an ancestral gene than the other cattle KIR genes. Cattle have multiple NKG2A-related genes and at least one Ly49 gene; thus, the data presented here suggest that they have the potential to express more major histocompatibility complex-binding NK receptors than other species.

Complexity in the cattle CD94/NKG2 gene families

Natural killer cell responses are controlled to a large extent by the interaction of an array of inhibitory and activating receptors with their ligands. The mostly nonpolymorphic CD94/NKG2 receptors in both humans and mice were shown to recognize a single nonclassical MHC class I molecule in each case. In this paper, we describe the CD94/NKG2 gene family in cattle. NKG2 and CD94 sequences were amplified from cDNA derived from four animals. Four CD94 sequences, ten NKG2A, and three NKG2C sequences were identified in total. In contrast to human, we show that cattle have multiple distinct NKG2A genes, some of which show minor allelic variation. All of the sequences designated NKG2A have two tyrosine-based inhibitory motifs in the cytoplasmic domain and one putative gene has, in addition, a charged residue in the transmembrane domain. NKG2C appears to be essentially monomorphic in cattle. All of the NKG2A sequences are similar apart from NKG2A-01, which, in contrast, shares the majority of its carbohydrate recognition domain with NKG2-C. Most of the genes appear to generate multiple alternatively spliced forms. These findings suggest that the CD94/NKG2A heterodimers in cattle, in contrast to other species, are binding several different ligands. Because NKG2C is not polymorphic, this raises questions as to the combined functional capacity of the CD94/NKG2 gene families in cattle.