Serological and molecular diversity in the cattle MHC class I region

A new method for ecologists to estimate heterozygote excess and deficit for multi-locus gene families

The fixation index, F IS, has been a staple measure to detect selection, or departures from random mating in populations. However, current Next Generation Sequencing (NGS) cannot easily estimate F IS, in multi-locus gene families that contain multiple loci having similar or identical arrays of variant sequences of ≥1 kilobase (kb), which differ at multiple positions. In these families, high-quality short-read NGS data typically identify variants, but not the genomic location, which is required to calculate F IS (based on locus-specific observed and expected heterozygosity). Thus, to assess assortative mating, or selection on heterozygotes, from NGS of multi-locus gene families, we need a method that does not require knowledge of which variants are alleles at which locus in the genome. We developed such a method. Like F IS, our novel measure, 1 H IS, is based on the principle that positive assortative mating, or selection against heterozygotes, and some other processes reduce within-individual variability relative to the population. We demonstrate high accuracy of 1 H IS on a wide range of simulated scenarios and two datasets from natural populations of penguins and dolphins. 1 H IS is important because multi-locus gene families are often involved in assortative mating or selection on heterozygotes. 1 H IS is particularly useful for multi-locus gene families, such as toll-like receptors, the major histocompatibility complex in animals, homeobox genes in fungi and self-incompatibility genes in plants.

Diversity and antigenic potentials of Mycoplasmopsis bovis secreted and outer membrane proteins within a core genome of strains isolated from North American bison and cattle

Mycoplasmopsis bovis is a worldwide economically important pathogen of cattle that can cause or indirectly contribute to bovine respiratory disease. M. bovis is also a primary etiological agent of respiratory disease in bison with high mortality rates. A major challenge in the development of an efficacious M. bovis vaccine is the design of antigens that contain both MHC-1 and MHC-2 T-cell epitopes, and that account for population level diversity within the species. Publicly available genomes and sequence read archive libraries of 381 M. bovis strains isolated from cattle (n = 202) and bison (n = 179) in North America were used to identify a core genome of 575 genes, including 38 that encode either known or predicted secreted or outer membrane proteins. The antigenic potentials of the proteins were characterized by the presence and strength of their T-cell epitopes, and their protein variant diversity at the population-level. The proteins had surprisingly low diversity and varying predictive levels of T-cell antigenicity. These results provide a reference for the selection or design of antigens for vaccine testing against strains infecting North American cattle and bison.

CD4 T Cell Responses to Theileria parva in Immune Cattle Recognize a Diverse Set of Parasite Antigens Presented on the Surface of Infected Lymphoblasts

Parasite-specific CD8 T cell responses play a key role in mediating immunity against Theileria parva in cattle (Bos taurus), and there is evidence that efficient induction of these responses requires CD4 T cell responses. However, information on the antigenic specificity of the CD4 T cell response is lacking. The current study used a high-throughput system for Ag identification using CD4 T cells from immune animals to screen a library of ∼40,000 synthetic peptides representing 499 T. parva gene products. Use of CD4 T cells from 12 immune cattle, representing 12 MHC class II types, identified 26 Ags. Unlike CD8 T cell responses, which are focused on a few dominant Ags, multiple Ags were recognized by CD4 T cell responses of individual animals. The Ags had diverse properties, but included proteins encoded by two multimember gene families: five haloacid dehalogenases and five subtelomere-encoded variable secreted proteins. Most Ags had predicted signal peptides and/or were encoded by abundantly transcribed genes, but neither parameter on their own was reliable for predicting antigenicity. Mapping of the epitopes confirmed presentation by DR or DQ class II alleles and comparison of available T. parva genome sequences demonstrated that they included both conserved and polymorphic epitopes. Immunization of animals with vaccine vectors expressing two of the Ags demonstrated induction of CD4 T cell responses capable of recognizing parasitized cells. The results of this study provide detailed insight into the CD4 T cell responses induced by T. parva and identify Ags suitable for use in vaccine development.

Granzyme B Is an Essential Mediator in CD8+ T Cell Killing of Theileria parva-Infected Cells

There is established evidence that cytotoxic CD8⁺ T cells are important mediators of immunity against the bovine intracellular protozoan parasite Theileria parva However, the mechanism by which the specific CD8⁺ T cells kill parasitized cells is not understood. Although the predominant pathway used by human and murine CD8⁺ T cells to kill pathogen-infected cells is granule exocytosis, involving the release of perforin and granzyme B, there is to date a lack of published information on the biological activities of bovine granzyme B. The present study set out to define the functional activities of bovine granzyme B and determine its role in mediating the killing of T. parva-parasitized cells. DNA constructs encoding functional and nonfunctional forms of bovine granzyme B were produced, and the proteins expressed in Cos-7 cells were used to establish an enzymatic assay to detect and quantify the expression of functional granzyme B protein. Using this assay, the levels of killing of different T. parva-specific CD8⁺ T cell clones were found to be significantly correlated with the levels of granzyme B protein but not the levels of mRNA transcript expression. Experiments using inhibitors specific for perforin and granzyme B confirmed that CD8⁺ T cell killing of parasitized cells is dependent on granule exocytosis and, specifically, granzyme B. Further studies showed that the granzyme B-mediated death of parasitized cells is independent of caspases and that granzyme B activates the proapoptotic molecule Bid.

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.

Identification of immediate early gene products of bovine herpes virus 1 (BHV-1) as dominant antigens recognized by CD8 T cells in immune cattle

In common with other herpes viruses, bovine herpes virus 1 (BHV-1) induces strong virus-specific CD8 T-cell responses. However, there is a paucity of information on the antigenic specificity of the responding T-cells. The development of a system to generate virus-specific CD8 T-cell lines from BHV-1-immune cattle, employing Theileria-transformed cell lines for antigen presentation, has enabled us to address this issue. Use of this system allowed the study to screen for CD8 T-cell antigens that are efficiently presented on the surface of virus-infected cells. Screening of a panel of 16 candidate viral gene products with CD8 T-cell lines from 3 BHV-1-immune cattle of defined MHC genotypes identified 4 antigens, including 3 immediate early (IE) gene products (ICP4, ICP22 and Circ) and a tegument protein (UL49). Identification of the MHC restriction specificities revealed that the antigens were presented by two or three class I MHC alleles in each animal. Six CD8 T-cell epitopes were identified in the three IE proteins by screening of synthetic peptides. Use of an algorithm (NetMHCpan) that predicts the peptide-binding characteristics of restricting MHC alleles confirmed and, in some cases refined, the identity of the epitopes. Analyses of the epitope specificity of the CD8 T-cell lines showed that a large component of the response is directed against these IE epitopes. The results indicate that these IE gene products are dominant targets of the CD8 T-cell response in BHV-I-immune cattle and hence are prime-candidate antigens for the generation of a subunit vaccine.

CD8 T-cell responses against the immunodominant Theileria parva peptide Tp249-59 are composed of two distinct populations specific for overlapping 11-mer and 10-mer epitopes

Immunity against Theileria parva is associated with CD8 T-cell responses that exhibit immunodominance, focusing the response against limited numbers of epitopes. As candidates for inclusion in vaccines, characterization of responses against immunodominant epitopes is a key component in novel vaccine development. We have previously demonstrated that the Tp249-59 and Tp1214-224 epitopes dominate CD8 T-cell responses in BoLA-A10 and BoLA-18 MHC I homozygous animals, respectively. In this study, peptide-MHC I tetramers for these epitopes, and a subdominant BoLA-A10-restricted epitope (Tp298-106 ), were generated to facilitate accurate and rapid enumeration of epitope-specific CD8 T cells. During validation of these tetramers a substantial proportion of Tp249-59 -reactive T cells failed to bind the tetramer, suggesting that this population was heterogeneous with respect to the recognized epitope. We demonstrate that Tp250-59 represents a distinct epitope and that tetramers produced with Tp50-59 and Tp49-59 show no cross-reactivity. The Tp249-59 and Tp250-59 epitopes use different serine residues as the N-terminal anchor for binding to the presenting MHC I molecule. Molecular dynamic modelling predicts that the two peptide-MHC I complexes adopt structurally different conformations and Tcell receptor β sequence analysis showed that Tp249-59 and Tp250-59 are recognized by non-overlapping T-cell receptor repertoires. Together these data demonstrate that although differing by only a single residue, Tp249-59 and Tp250-59 epitopes form distinct ligands for T-cell receptor recognition. Tetramer analysis of T. parva-specific CD8 T-cell lines confirmed the immunodominance of Tp1214-224 in BoLA-A18 animals and showed in BoLA-A10 animals that the Tp249-59 epitope response was generally more dominant than the Tp250-59 response and confirmed that the Tp298-106 response was subdominant.

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.

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.

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.

Reproduction of bovine neonatal pancytopenia (BNP) by feeding pooled colostrum reveals variable alloantibody damage to different haematopoietic lineages

Bovine neonatal pancytopenia (BNP) is a recently described haemorrhagic disease of calves characterised by thrombocytopenia, leucopenia and bone marrow depletion. Feeding colostrum from cows that have previously produced a BNP affected calf has been shown to induce the disease in some calves, leading to the hypothesis that alloantibodies in colostrum from dams of affected calves mediate destruction of blood and bone marrow cells in the recipient calves. The aims of the current experimental study were first to confirm the role of colostrum-derived antibody in mediating the disease and second to investigate the haematopoietic cell lineages and maturation stages depleted by the causative antibodies. Clinical, haematological and pathological changes were examined in 5 calves given a standardised pool of colostrum from known BNP dams, and 5 control calves given an equivalent pool of colostrum from non-BNP dams. All calves fed challenge colostrum showed progressive depletion of bone marrow haematopoietic cells and haematological changes consistent with the development of BNP. Administration of a standardised dose of the same colostrum pool to each calf resulted in a consistent response within the groups, allowing detailed interpretation of the cellular changes not previously described. Analyses of blood and serial bone marrow changes revealed evidence of differential effects on different blood cell lineages. Peripheral blood cell depletion was confined to leucocytes and platelets, while bone marrow damage occurred to the primitive precursors and lineage committed cells of the thrombocyte, lymphocyte and monocyte lineages, but only to the more primitive precursors in the neutrophil, erythrocyte and eosinophil lineages. Such differences between lineages may reflect cell type-dependent differences in levels of expression or conformational nature of the target antigens.

Theileria annulata-transformed cell lines are efficient antigen-presenting cells for in vitro analysis of CD8 T cell responses to bovine herpesvirus-1

Continuously growing cell lines infected with the protozoan parasite Theileria annulata can readily be established by in vitro infection of leukocytes with the sporozoite stage of the parasite. The aim of the current study was to determine whether such transformed cell lines could be used as antigen presenting cells to analyse the antigenic specificity of bovine CD8 T cell responses to viral infections. Bovine herpes virus 1 (BHV-1), which is known to induce CD8 T cell responses, was used as a model. T. annulata- transformed cells were shown to express high levels of CD40 and CD80 and were susceptible to infection with BHV-1, vaccinia and canarypox viruses. The capacity of the cells to generate antigen-specific CD8 T cell lines was initially validated using a recombinant canarypox virus expressing a defined immunodominant T. parva antigen (Tp1). Autologous T. annulata-transformed cells infected with BHV-1 were then used successfully to generate specific CD8 T cell lines and clones from memory T cell populations of BHV-1-immune animals. These lines were BHV-1-specific and class I MHC-restricted. In contrast to previous studies, which reported recognition of the glycoproteins gB and gD, the CD8 T cell lines generated in this study did not recognise these glycoproteins. Given the ease with which T. annulata-transformed cell lines can be established and maintained in vitro and their susceptibility to infection with poxvirus vectors, these cell lines offer a convenient and efficient in vitro system to analyse the fine specificity of virus-specific CD8 T cell responses in cattle.

Interspecific hybridization increases MHC class II diversity in two sister species of newts

Our understanding of the evolutionary mechanisms generating variation within the highly polymorphic major histocompatibility complex (MHC) genes remains incomplete. Assessing MHC variation across multiple populations, of recent and ancient divergence, may facilitate understanding of geographical and temporal aspects of variation. Here, we applied 454 sequencing to perform a large-scale, comprehensive analysis of MHC class II in the closely related, hybridizing newts, Lissotriton vulgaris (Lv) and Lissotriton montandoni (Lm). Our study revealed an extensive (299 alleles) geographically structured polymorphism. Populations at the southern margin of the Lv distribution, inhabited by old and distinct lineages (southern Lv), exhibited moderate MHC variation and strong population structure, indicating little gene flow or extensive local adaptation. Lissotriton vulgaris in central Europe and the northern Balkans (northern Lv) and almost all Lm populations had a high MHC variation. A much higher proportion of MHC alleles was shared between Lm and northern Lv than between Lm and southern Lv. Strikingly, the average pairwise F(ST) between northern Lv and Lm was significantly lower than between northern and southern Lv for MHC, but not for microsatellites. Thus, high MHC variation in Lm and northern Lv may result from gene flow between species. We hypothesize that the interspecific exchange of MHC genes may be facilitated by frequency-dependent selection. A marginally significant correlation between the MHC and microsatellite allelic richness indicates that demographic factors may have contributed to the present-day pattern of MHC variation, but unequivocal signatures of adaptive evolution in MHC class II sequences emphasize the role of selection on a longer timescale.

Escape from CD8+ T cell response by natural variants of an immunodominant epitope from Theileria parva is predominantly due to loss of TCR recognition

Polymorphism of immunodominant CD8(+) T cell epitopes can facilitate escape from immune recognition of pathogens, leading to strain-specific immunity. In this study, we examined the TCR β-chain (TRB) diversity of the CD8(+) T cell responses of cattle against two immunodominant epitopes from Theileria parva (Tp1(214-224) and Tp2(49-59)) and investigated the role of TCR recognition and MHC binding in determining differential recognition of a series of natural variants of the highly polymorphic Tp2(49-59) epitope by CD8(+) T cell clones of defined TRB genotype. Our results show that both Tp1(214-224) and Tp2(49-59) elicited CD8(+) T cell responses using diverse TRB repertoires that showed a high level of stability following repeated pathogenic challenge over a 3-y period. Analysis of single-alanine substituted versions of the Tp2(49-59) peptide demonstrated that Tp2(49-59)-specific clonotypes had a broad range of fine specificities for the epitope. Despite this diversity, all natural variants exhibited partial or total escape from immune recognition, which was predominantly due to abrogation of TCR recognition, with mutation resulting in loss of the lysine residue at P8, playing a particularly dominant role in escape. The levels of heterozygosity in individual Tp2(49-59) residues correlated closely with loss of immune recognition, suggesting that immune selection has contributed to epitope polymorphism.

Demonstration of strain-specific CD8 T cell responses to Theileria annulata

The present study set out to examine the nature and specificity of the bovine CD8 T cell response at the clonal level in a group of eight animals immunized with a cloned population of Theileria annulata. The results demonstrated that immunized animals generated parasite-specific CD8 T cells that produced IFNγ in response to parasite stimulation but had highly variable levels of cytotoxicity for parasitized cells. The study also demonstrated that these parasite-specific CD8 T cells could be propagated and cloned in vitro from the memory T cell pool of cattle immunized with live T. annulata parasites. Within the small group of animals studied, there was evidence that responses were preferentially directed to antigens presented by an A10⁺ class I major histocompatibility complex (MHC) haplotype, suggesting that responses restricted by products of this haplotype may be dominant. The A10-restricted responses showed differential recognition of different parasite isolates and clones. By using a cloned population of parasites both for immunization of the animals and for in vitro analyses of the responses, we obtained unambiguous evidence that at least a proportion of CD8 T cells restricted by one MHC haplotype were parasite strain restricted.

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.

454 sequencing reveals extreme complexity of the class II Major Histocompatibility Complex in the collared flycatcher

Background

Because of their functional significance, the Major Histocompatibility Complex (MHC) class I and II genes have been the subject of continuous interest in the fields of ecology, evolution and conservation. In some vertebrate groups MHC consists of multiple loci with similar alleles; therefore, the multiple loci must be genotyped simultaneously. In such complex systems, understanding of the evolutionary patterns and their causes has been limited due to challenges posed by genotyping.

Results

Here we used 454 amplicon sequencing to characterize MHC class IIB exon 2 variation in the collared flycatcher, an important organism in evolutionary and immuno-ecological studies. On the basis of over 152,000 sequencing reads we identified 194 putative alleles in 237 individuals. We found an extreme complexity of the MHC class IIB in the collared flycatchers, with our estimates pointing to the presence of at least nine expressed loci and a large, though difficult to estimate precisely, number of pseudogene loci. Many similar alleles occurred in the pseudogenes indicating either a series of recent duplications or extensive concerted evolution. The expressed alleles showed unambiguous signals of historical selection and the occurrence of apparent interlocus exchange of alleles. Placing the collared flycatcher's MHC sequences in the context of passerine diversity revealed transspecific MHC class II evolution within the Muscicapidae family.

Conclusions

454 amplicon sequencing is an effective tool for advancing our understanding of the MHC class II structure and evolutionary patterns in Passeriformes. We found a highly dynamic pattern of evolution of MHC class IIB genes with strong signals of selection and pronounced sequence divergence in expressed genes, in contrast to the apparent sequence homogenization in pseudogenes. We show that next generation sequencing offers a universal, affordable method for the characterization and, in perspective, genotyping of MHC systems of virtually any complexity.

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.

CD8+ T-cell responses to Theileria parva are preferentially directed to a single dominant antigen: Implications for parasite strain-specific immunity

Although immunodominance of CD8⁺ T-cell responses is a well-recognised feature of viral infections, its role in responses to more antigenically complex pathogens is less clear. In previous studies we have observed that CD8⁺ T-cell responses to Theileria parva exhibit different patterns of parasite strain specificity in cattle of different MHC genotypes. In the current study, we demonstrated that animals homozygous for the A10 and A18 MHC haplotypes have detectable responses to only one of 5 T. parva antigens. Over 60% of the responding T cells from the A18⁺ and A10⁺ animals recognised defined epitopes in the Tp1 and Tp2 antigens, respectively. Comparison of T-cell receptor β chain expression profiles of CD8⁺ T-cell lines and CD8⁺ T cells harvested ex vivo confirmed that the composition of the T-cell lines was representative of the in vivo memory CD8⁺ T-cell populations. Analysis of the Tp1 and Tp2 antigens revealed sequence polymorphism, which was reflected by differential recognition by T-cell lines. In conclusion, we have demonstrated a profound immunodominance in the CD8⁺ T-cell response to T. parva, which we propose is a major determinant of the parasite strain specificity of the response and hence immune protection.

Dissection of the clonal composition of bovine alphabeta T cell responses using T cell receptor Vbeta subfamily-specific PCR and heteroduplex analysis

Although techniques that permit analysis of the clonal composition of T cell populations have been used extensively to provide a better understanding of the mechanisms that influence efficacy of T cell responses in humans and mice, such methods are lacking for other animal species. In this paper we report the establishment and validation of a panel of Vbeta subfamily-specific semi-nested PCR assays, and a CDR3beta heteroduplex technique for analysing the clonal diversity of bovine alphabeta T cell responses. Development of these methods was based on available sequence data for 48 functional Vbeta genes classified within 17 subfamilies. These techniques were used to determine the clonal composition of parasite-reactive CD8(+) T cells obtained from two animals immunised with the protozoan parasite Theileria parva. Analyses of uncloned T cell lines as well as large panels of cloned T cells derived from each of these lines confirmed the specificity and sensitivity of the assays. Specific PCR products were obtained from 96% of the T cell clones examined, indicating that the currently identified Vbeta genes represent most of the functional Vbeta subfamilies in cattle. Heteroduplex analyses, coupled with sequencing of PCR products, identified over 20 clonal expansions within each of the T cell lines, distributed over a large number of Vbeta subfamilies, although a limited number of clonotypes numerically dominated the response in both animals. The development and validation of these methods provides for the first time a generic set of molecular tools that can be used to perform detailed analysis of the TCR diversity and clonal composition of bovine T cell responses.

Characterization of the fine specificity of bovine CD8 T-cell responses to defined antigens from the protozoan parasite Theileria parva

Immunity against the bovine intracellular protozoan parasite Theileria parva has been shown to be mediated by CD8 T cells. Six antigens targeted by CD8 T cells from T. parva-immune cattle of different major histocompatibility complex (MHC) genotypes have been identified, raising the prospect of developing a subunit vaccine. To facilitate further dissection of the specificity of protective CD8 T-cell responses and to assist in the assessment of responses to vaccination, we set out to identify the epitopes recognized in these T. parva antigens and their MHC restriction elements. Nine epitopes in six T. parva antigens, together with their respective MHC restriction elements, were successfully identified. Five of the cytotoxic-T-lymphocyte epitopes were found to be restricted by products of previously described alleles, and four were restricted by four novel restriction elements. Analyses of CD8 T-cell responses to five of the epitopes in groups of cattle carrying the defined restriction elements and immunized with live parasites demonstrated that, with one exception, the epitopes were consistently recognized by animals of the respective genotypes. The analysis of responses was extended to animals immunized with multiple antigens delivered in separate vaccine constructs. Specific CD8 T-cell responses were detected in 19 of 24 immunized cattle. All responder cattle mounted responses specific for antigens for which they carried an identified restriction element. By contrast, only 8 of 19 responder cattle displayed a response to antigens for which they did not carry an identified restriction element. These data demonstrate that the identified antigens are inherently dominant in animals with the corresponding MHC genotypes.

Identification of a novel foot-and-mouth disease virus specific T-cell epitope with immunodominant characteristics in cattle with MHC serotype A31

To identify foot-and-mouth disease virus (FMDV) specific T-cell epitopes within the entire polyprotein sequence of the virus, 442 overlapping pentadecapeptides were tested in proliferation assays using lymphocytes from cattle experimentally infected with FMDV. Four months post-infection cells from all investigated animals (n = 4) responded by proliferation and interferon-gamma production to a peptide located on the structural protein 1D (VP1), amino acid residues 66-80. Major histocompatibility complex (MHC) serotyping of the investigated cattle indicated that all animals shared the MHC serotype A31 which comprises the class II allele DRB3 0701. This may explain the common recognition of this newly discovered epitope. Responses to other peptides could only be observed in one animal and rapidly declined during the time course of the study. These observations point to an immunodominant role of this epitope located on the protein 1D in cattle with MHC serotype A31.

The origins of polypeptide domains

Three decades ago Gilbert posited that novel proteins arise by re-shuffling genomic sequences encoding polypeptide domains. Today, with numerous genomes and countless genes sequenced, it is well established that recombination of sequences encoding polypeptide domains plays a major role in protein evolution. There is, however, less evidence to suggest how the novel polypeptide domains, themselves, arise. Recent comparisons of genomes from closely related species have revealed numerous species-specific exons, supporting models of domain origin based on "exonization" of intron sequences. Also, a mechanism for the origin of novel polypeptide domains has been proposed based on analyses of insertion-based polymorphisms between orthologous genes across broad phylogenetic spectra and between allelic variants of genes within species. This review discusses these processes and how each might participate in the evolutionary emergence of novel polypeptide domains.

BoLA class I allele diversity and polymorphism in a herd of cattle

Major histocompatibility complex class I genes are among the most polymorphic genes characterized. The high level of polymorphism is essential for generating host immune responses. In humans, three distinct genomic loci encode human leukocyte antigen (HLA) class I genes, allowing individuals to express up to six different HLA class I molecules. In cattle, the number of distinct genomic loci are currently at least six, and the number of different bovine leukocyte antigens (BoLA) class I molecules that are expressed in individual animals are variable. The extent of allele variation within the cattle population is unknown. In this study, the number and variety of BoLA class I sequences expressed by 36 individuals were determined from full-length BoLA class I cDNA clones. Twenty distinct BoLA class I alleles were identified, with only four being previously reported. The number of expressed BoLA class I alleles in individual animals ranged between one and four, with none of the animals having an identical complement of BoLA class I molecules. Variation existed in the number of BoLA class I alleles expressed as well as the composition of expressed alleles, however, several BoLA class I alleles were found in multiple individual animals. Polymorphic amino acid sites were analyzed for positive and negative selection using the ADAPTSITE program. In the antigen recognition sites (ARS), there were eight positions that were predicted to be under positive selection and three positions that were predicted to be under negative selection from 62 positions. In contrast, for non-antigen recognition sites (non-ARS), there were three positions that were predicted to be under positive selection and 20 that were predicted to be under negative selection from 278, indicating that positive selection of amino acids occurs at a greater frequency within the antigen recognition sites.

The E5 oncoprotein of BPV-4 does not interfere with the biosynthetic pathway of non-classical MHC class I

The major histocompatibility complex (MHC) class I region in mammals contains both classical and non-classical MHC class I genes. Classical MHC class I molecules present antigenic peptides to cytotoxic T lymphocytes, whereas non-classical MHC class I molecules have a variety of functions. Both classical and non-classical MHC molecules interact with natural killer cell receptors and may under some circumstances prevent cell death by natural killer cytotoxicity. The E5 oncoprotein of BPV-4 down-regulates the expression of classical MHC class I on the cell surface and retains the complex in the Golgi apparatus. The inhibition of classical MHC class I to the cell surface results from both the impaired acidification of the Golgi, due to the interaction of E5 with subunit c of the H+ V-ATPase, and to the physical binding of E5 to the heavy chain of MHC class I. Despite the profound effect of E5 on classical MHC class I, E5 does not retain a non-classical MHC class I in the Golgi, does not inhibit its transport to the cell surface and does not bind its heavy chain. We conclude that, as is the case for HPV-16 E5, BPV-4 E5 does not down-regulate certain non-classical MHC class I, potentially providing a mechanism for the escape of the infected cell from attack by both cytotoxic T lymphocytes and NK cells.

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

Major histocompatibility complex (MHC) class I genes play a crucial role in the immune defence against intracellular pathogens. An important evolutionary strategy is to generate and maintain a high level of diversity in these genes. Humans express three highly polymorphic classical MHC class I genes (HLA-A, HLA-B and HLA-C). In contrast, some species, for example rat and rhesus macaque, maintain diversity by generation of haplotypes that vary considerably with regard to the number and combination of transcribed genes. Cattle appear to use both strategies. We show that various combinations of six apparently classical genes, three of which are highly polymorphic, are transcribed on different haplotypes. Although additional sequences were identified in both cDNA and gDNA, it was not possible to assign them to any of these defined genes. Most were highly divergent or were non-classical class I genes. Thus, we found little evidence for frequent duplication and deletion of classical class I genes as reported in some other species. However, the maintenance of class I diversity in cattle may involve limited gene shuffling and deletion, possibly as a result of unequal crossing-over within the class I region.