The DY genes of the cattle MHC: expression and comparative analysis of an unusual class II MHC gene pair

The DRB3 gene of the bovine major histocompatibility complex: discovery, diversity and distribution of alleles in commercial breeds of cattle and applications for development of vaccines

The bovine Major Histocompatibility Complex (MHC), also known as the Bovine Leucocyte Antigen (BoLA) complex, is the genomic region that encodes the most important molecules for antigen presentation to initiate immune responses. The first evidence of MHC in bovines pointed to a locus containing 2 antigens, one detected by cytotoxic antiserum (MHC class I) and another studied by mixed lymphocyte culture tests (MHC class II). The most studied gene in the BoLA region is the highly polymorphic BoLA-DRB3, which encodes a β chain with a peptide groove domain involved in antigen presentation for T cells that will develop and co-stimulate cellular and humoral effector responses. BoLA-DRB3 alleles have been associated with outcomes in infectious diseases such as mastitis, trypanosomiasis, and tick loads, and with production traits. To catalog these alleles, 2 nomenclature methods were proposed, and the current use of both systems makes it difficult to list, comprehend and apply these data effectively. In this review we have organized the knowledge available in all of the reports on the frequencies of BoLA-DRB3 alleles. It covers information from studies made in at least 26 countries on more than 30 breeds; studies are lacking in countries that are important producers of cattle livestock. We highlight practical applications of BoLA studies for identification of markers associated with resistance to infectious and parasitic diseases, increased production traits and T cell epitope mapping, in addition to genetic diversity and conservation studies of commercial and creole and locally adapted breeds. Finally, we provide support for the need of studies to discover new BoLA alleles and uncover unknown roles of this locus in production traits.

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.

Organization of the Addax Major Histocompatibility Complex Provides Insights Into Ruminant Evolution

Ruminants are critical as prey in transferring solar energy fixed by plants into carnivorous species, yet the genetic signature of the driving forces leading to the evolutionary success of the huge number of ruminant species remains largely unknown. Here we report a complete DNA map of the major histocompatibility complex (MHC) of the addax (Addax nasomaculatus) genome by sequencing a total of 47 overlapping BAC clones previously mapped to cover the MHC region. The addax MHC is composed of 3,224,151 nucleotides, harboring a total of 150 coding genes, 50 tRNA genes, and 14 non-coding RNA genes. The organization of addax MHC was found to be highly conserved to those of sheep and cattle, highlighted by a large piece of chromosome inversion that divided the MHC class II into IIa and IIb subregions. It is now highly possible that all of the ruminant species in the family of Bovidae carry the same chromosome inversion in the MHC region, inherited from a common ancestor of ruminants. Phylogenetic analysis indicated that DY, a ruminant-specific gene located at the boundary of the inversion and highly expressed in dendritic cells, was possibly evolved from DQ, with an estimated divergence time ~140 million years ago. Homology modeling showed that the overall predicted structure of addax DY was similar to that of HLA-DQ2. However, the pocket properties of P1, P4, P6, and P9, which were critical for antigen binding in the addax DY, showed certain distinctive features. Structural analysis suggested that the populations of peptide antigens presented by addax DY and HLA-DQ2 were quite diverse, which in theory could serve to promote microbial regulation in the rumen by ruminant species, contributing to enhanced grass utilization ability. In summary, the results of our study helped to enhance our understanding of the MHC evolution and provided additional supportive evidence to our previous hypothesis that an ancient chromosome inversion in the MHC region of the last common ancestor of ruminants may have contributed to the evolutionary success of current ruminants on our planet.

A high-density BAC physical map covering the entire MHC region of addax antelope genome

BACKGROUND

The mammalian major histocompatibility complex (MHC) harbours clusters of genes associated with the immunological defence of animals against infectious pathogens. At present, no complete MHC physical map is available for any of the wild ruminant species in the world.

RESULTS

The high-density physical map is composed of two contigs of 47 overlapping bacterial artificial chromosome (BAC) clones, with an average of 115 Kb for each BAC, covering the entire addax MHC genome. The first contig has 40 overlapping BAC clones covering an approximately 2.9 Mb region of MHC class I, class III, and class IIa, and the second contig has 7 BAC clones covering an approximately 500 Kb genomic region that harbours MHC class IIb. The relative position of each BAC corresponding to the MHC sequence was determined by comparative mapping using PCR screening of the BAC library of 192,000 clones, and the order of BACs was determined by DNA fingerprinting. The overlaps of neighboring BACs were cross-verified by both BAC-end sequencing and co-amplification of identical PCR fragments within the overlapped region, with their identities further confirmed by DNA sequencing.

CONCLUSIONS

We report here the successful construction of a high-quality physical map for the addax MHC region using BACs and comparative mapping. The addax MHC physical map we constructed showed one gap of approximately 18 Mb formed by an ancient autosomal inversion that divided the MHC class II into IIa and IIb. The autosomal inversion provides compelling evidence that the MHC organizations in all of the ruminant species are relatively conserved.

Fine-mapping host genetic variation underlying outcomes to Mycobacterium bovis infection in dairy cows

BACKGROUND

Susceptibility to Mycobacterium bovis infection in cattle is governed in part by host genetics. However, cattle diagnosed as infected with M. bovis display varying signs of pathology. The variation in host response to infection could represent a continuum since time of exposure or distinct outcomes due to differing pathogen handling. The relationships between host genetics and variation in host response and pathological sequelae following M. bovis infection were explored by genotyping 1966 Holstein-Friesian dairy cows at 538,231 SNPs with three distinct phenotypes. These were: single intradermal cervical comparative tuberculin (SICCT) test positives with visible lesions (VLs), SICCT-positives with undetected visible lesions (NVLs) and matched controls SICCT-negative on multiple occasions.

RESULTS

Regional heritability mapping identified three loci associated with the NVL phenotype on chromosomes 17, 22 and 23, distinct to the region on chromosome 13 associated with the VL phenotype. The region on chromosome 23 was at genome-wide significance and candidate genes overlapping the mapped window included members of the bovine leukocyte antigen class IIb region, a complex known for its role in immunity and disease resistance. Chromosome heritability analysis attributed variance to six and thirteen chromosomes for the VL and NVL phenotypes, respectively, and four of these chromosomes were found to explain a proportion of the phenotypic variation for both the VL and NVL phenotype. By grouping the M. bovis outcomes (VLs and NVLs) variance was attributed to nine chromosomes. When contrasting the two M. bovis infection outcomes (VLs vs NVLs) nine chromosomes were found to harbour heritable variation. Regardless of the case phenotype under investigation, chromosome heritability did not exceed 8% indicating that the genetic control of bTB resistance consists of variants of small to moderate effect situated across many chromosomes of the bovine genome.

CONCLUSIONS

These findings suggest the host genetics of M. bovis infection outcomes is governed by distinct and overlapping genetic variants. Thus, variation in the pathology of M. bovis infected cattle may be partly genetically determined and indicative of different host responses or pathogen handling. There may be at least three distinct outcomes following M. bovis exposure in dairy cattle: resistance to infection, infection resulting in pathology or no detectable pathology.

What chickens might tell us about the MHC class II system

Almost all knowledge about the structure and function of MHC class II molecules outside of mammals comes from work with chickens. Most of the genes implicated in the class II system are present in chickens, so it is likely that the machinery of antigen processing and peptide-loading is similar to mammals. However, there is only one isotype (lineage) of classical class II genes, with one monomorphic DR-like BLA gene and two polymorphic BLB genes, located near one DMA and two DMB genes. The DMB2 and BLB2 genes are widely expressed at high levels, whereas the DMB1 and BLB1 genes are only expressed at highest levels in spleen and intestine, suggesting the possibility of two class II systems in chickens.

Selection and trans-species polymorphism of major histocompatibility complex class II genes in the order Crocodylia

Major Histocompatibility Complex (MHC) class II genes encode for molecules that aid in the presentation of antigens to helper T cells. MHC characterisation within and between major vertebrate taxa has shed light on the evolutionary mechanisms shaping the diversity within this genomic region, though little characterisation has been performed within the Order Crocodylia. Here we investigate the extent and effect of selective pressures and trans-species polymorphism on MHC class II α and β evolution among 20 extant species of Crocodylia. Selection detection analyses showed that diversifying selection influenced MHC class II β diversity, whilst diversity within MHC class II α is the result of strong purifying selection. Comparison of translated sequences between species revealed the presence of twelve trans-species polymorphisms, some of which appear to be specific to the genera Crocodylus and Caiman. Phylogenetic reconstruction clustered MHC class II α sequences into two major clades representing the families Crocodilidae and Alligatoridae. However, no further subdivision within these clades was evident and, based on the observation that most MHC class II α sequences shared the same trans-species polymorphisms, it is possible that they correspond to the same gene lineage across species. In contrast, phylogenetic analyses of MHC class II β sequences showed a mixture of subclades containing sequences from Crocodilidae and/or Alligatoridae, illustrating orthologous relationships among those genes. Interestingly, two of the subclades containing sequences from both Crocodilidae and Alligatoridae shared specific trans-species polymorphisms, suggesting that they may belong to ancient lineages pre-dating the divergence of these two families from the common ancestor 85-90 million years ago. The results presented herein provide an immunogenetic resource that may be used to further assess MHC diversity and functionality in Crocodylia.

The major histocompatibility complex in bovines: a review

Productivity in dairy cattle and buffaloes depends on the genetic factors governing the production of milk and milk constituents as well as genetic factors controlling disease resistance or susceptibility. The immune system is the adaptive defense system that has evolved in vertebrates to protect them from invading pathogens and also carcinomas. It is remarkable in the sense that it is able to generate an enormous variety of cells and biomolecules which interact with each other in numerous ways to form a complex network that helps to recognize, counteract, and eliminate the apparently limitless number of foreign invading pathogens/molecules. The major histocompatibility complex which is found to occur in all mammalian species plays a central role in the development of the immune system. It is an important candidate gene involved in susceptibility/resistance to various diseases. It is associated with intercellular recognition and with self/nonself discrimination. It plays major role in determining whether transplanted tissue will be accepted as self or rejected as foreign.

Allelic Diversity of Major Histocompatibility Complex Class II DRB Gene in Indian Cattle and Buffalo

The present study was conducted to study the diversity of MHC-DRB3 alleles in Indian cattle and buffalo breeds. Previously reported BoLA-DRB exon 2 alleles of Indian Zebu cattle, Bos taurus cattle, buffalo, sheep, and goats were analyzed for the identities and divergence among various allele sequences. Comparison of predicted amino acid residues of DRB3 exon 2 alleles with similar alleles from other ruminants revealed considerable congruence in amino acid substitution pattern. These alleles showed a high degree of nucleotide and amino acid polymorphism at positions forming peptide-binding regions. A higher rate of nonsynonymous substitution was detected at the peptide-binding regions, indicating that BoLA-DRB3 allelic sequence evolution was driven by positive selection.

Giant panda genomic data provide insight into the birth-and-death process of mammalian major histocompatibility complex class II genes

To gain an understanding of the genomic structure and evolutionary history of the giant panda major histocompatibility complex (MHC) genes, we determined a 636,503-bp nucleotide sequence spanning the MHC class II region. Analysis revealed that the MHC class II region from this rare species contained 26 loci (17 predicted to be expressed), of which 10 are classical class II genes (1 DRA, 2 DRB, 2 DQA, 3 DQB, 1 DYB, 1 DPA, and 2 DPB) and 4 are non-classical class II genes (1 DOA, 1 DOB, 1 DMA, and 1 DMB). The presence of DYB, a gene specific to ruminants, prompted a comparison of the giant panda class II sequence with those of humans, cats, dogs, cattle, pigs, and mice. The results indicated that birth and death events within the DQ and DRB-DY regions led to major lineage differences, with absence of these regions in the cat and in humans and mice respectively. The phylogenetic trees constructed using all expressed alpha and beta genes from marsupials and placental mammals showed that: (1) because marsupials carry loci corresponding to DR, DP, DO and DM genes, those subregions most likely developed before the divergence of marsupials and placental mammals, approximately 150 million years ago (MYA); (2) conversely, the DQ and DY regions must have evolved later, but before the radiation of placental mammals (100 MYA). As a result, the typical genomic structure of MHC class II genes for the giant panda is similar to that of the other placental mammals and corresponds to BTNL2∼DR1∼DQ∼DR2∼DY∼DO_box∼DP∼COL11A2. Over the past 100 million years, there has been birth and death of mammalian DR, DQ, DY, and DP genes, an evolutionary process that has brought about the current species-specific genomic structure of the MHC class II region. Furthermore, facing certain similar pathogens, mammals have adopted intra-subregion (DR and DQ) and inter-subregion (between DQ and DP) convergent evolutionary strategies for their alpha and beta genes, respectively.

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

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

Comparative analysis of the bovine MHC class IIb sequence identifies inversion breakpoints and three unexpected genes

The bovine major histocompatibility complex (MHC) or BoLA is organized differently from typical mammalian MHCs in that a large portion of the class II region, called class IIb, has been transposed to a position near the centromere on bovine chromosome 23. Gene mapping indicated that the rearrangement resulted from a single inversion, but the boundaries and gene content of the inverted segment have not been fully determined. Here, we report the genomic sequence of BoLA IIb. Comparative sequence analysis with the human MHC revealed that the proximal inversion breakpoint occurred approximately 2.5 kb from the 3' end of the glutamate-cysteine ligase, catalytic subunit (GCLC) locus and that the distal breakpoint occurred about 2 kb from the 5' end from a divergent class IIDRbeta-like sequence designated DSB. Gene content, order and orientation of BoLA IIb are consistent with the single inversion hypothesis when compared with the corresponding region of the human class II MHC (HLA class II). Differences with HLA include the presence of a single histone H2B gene located between the proteasome subunit, beta type, 9 (PSMB9) and DMB loci and a duplicated TAP2 with a variant splice site. BoLA IIb spans approximately 450 kb DNA, with 20 apparently intact genes and no obvious pseudogenes. The region contains 227 simple sequence repeats (SSRs) and approximately 167 kb of retroviral-related repetitive DNA. Nineteen of the 20 genes identified in silico are supported by bovine EST data indicating that the functional gene content of BoLA IIb has not been diminished because it has been transposed from the remainder of BoLA genes.

The genomic sequence and analysis of the swine major histocompatibility complex

We describe the generation and analysis of an integrated sequence map of a 2.4-Mb region of pig chromosome 7, comprising the classical class I region, the extended and classical class II regions, and the class III region of the major histocompatibility complex (MHC), also known as swine leukocyte antigen (SLA) complex. We have identified and manually annotated 151 loci, of which 121 are known genes (predicted to be functional), 18 are pseudogenes, 8 are novel CDS loci, 3 are novel transcripts, and 1 is a putative gene. Nearly all of these loci have homologues in other mammalian genomes but orthologues could be identified with confidence for only 123 genes. The 28 genes (including all the SLA class I genes) for which unambiguous orthology to genes within the human reference MHC could not be established are of particular interest with respect to porcine-specific MHC function and evolution. We have compared the porcine MHC to other mammalian MHC regions and identified the differences between them. In comparison to the human MHC, the main differences include the absence of HLA-A and other class I-like loci, the absence of HLA-DP-like loci, and the separation of the extended and classical class II regions from the rest of the MHC by insertion of the centromere. We show that the centromere insertion has occurred within a cluster of BTNL genes located at the boundary of the class II and III regions, which might have resulted in the loss of an orthologue to human C6orf10 from this region.

Nomenclature for factors of the swine leukocyte antigen class II system, 2005

A systematic nomenclature for the genes and alleles of the swine major histocompatibility complex (MHC) is essential to the development and communication of research in swine immunology. The Swine Leukocyte Antigen (SLA) Nomenclature Committee of the International Society for Animal Genetics (ISAG) has reviewed all of the DNA-sequence information for MHC class II genes, available in GenBank/EMBL/DDBJ databases, and the associated published reports to develop such a systematic nomenclature. This article summarizes the proposed nomenclature, which parallels the World Health Organization's nomenclature for factors of the human MHC. The SLA class II genes expressed on the cell membrane will be noted as SLA-DRA, SLA-DRB1, SLA-DQA, and SLA-DQB1. Nomenclature assignments for all SLA class II GenBank sequences are now noted. The committee will add new SLA class II allele designations, as they are discovered, and will maintain a publicly available list of all recognized genes and alleles using the Immuno Polymorphism Database (IPD). The sequences will be available from the IPD-MHC section of the database which contains non-human MHC sequences (http://www.ebi.ac.uk/ipd/mhc/sla/).

Conservation of promoter, coding and intronic regions of the non-classical MHC class II DYA gene suggests evolution under functional constraints

The major histocompatibility complex (MHC) in ruminants contains a unique pair of class II genes (DYA and DYB) of unknown function. As functional genes show higher levels of nucleotide conservation than pseudogenes we compared the DYA genes from sheep and cattle, two species which diverged from a common ancestor approximately 20 million years ago. Comparative analysis identified levels of nucleotide conservation in immediate promoter (97%), coding (94%) and intronic regions (91%) comparable with functional MHC genes. The Ovar-DYA transcript revealed an open reading frame encoding a 288 amino acid protein compared with a 253 amino acid protein associated with the BoLA-DYA transcript. A dinucleotide deletion in exon 4 of the Ovar-DYA transcript combined with alternative exon 5 splice sites introduces unusual diversity to the cytoplasmic domain of the Ovar-DYalpha polypeptide. The degree of conservation between these class II MHC genes is consistent with evolution under purifying selection suggesting that these genes retain a unique function in ruminants.

Comparative genomics of major histocompatibility complexes

The major histocompatibility complex (MHC) is a gene dense region found in all jawed vertebrates examined to date. The MHC contains a high percentage of immune genes, in particular genes involved in antigen presentation, which are generally highly polymorphic. The region plays an important role in disease resistance. The clustering of MHC genes could be advantageous for co-evolution or regulation, and its study in many species is desirable. Even though some linkage of MHC genes is apparent in all gnathostomes, the genomic organization can differ greatly by species, suggesting rapid evolution of MHC genes after divergence from a common ancestor. Previous reviews of comparative MHC organization have been written when relatively fragmentary sequence and mapping data were available on many species. This review compares maps of MHC gene orders in commonly studied species, where extensive sequencing has been performed.

The cattle major histocompatibility complex: is it unique?

Major histocompatibility complex (MHC) class I genes encode highly polymorphic molecules that are expressed on virtually every cell type, and have been identified in all but the most primitive vertebrates. They play a number of crucial roles in the immune response to infectious disease. Most information regarding MHC genes has been generated from humans and mice but, because of the great variability found in the MHC system, it is not always possible to extrapolate from these to other species. Many strategies have evolved to maximise the ability of the MHC to protect individuals and populations against pathogens. Cattle MHC class I genes exhibit a number of unusual features. Evidence from mapping studies, haplotype and phylogenetic analyses suggests the presence of six classical class I loci, in contrast to the more usual two or three, and these are expressed in various combinations of one, two or three on different haplotypes. Although it remains difficult to assign alleles to loci, it appears that none of the loci are expressed on all haplotypes. There is currently limited information relating to polymorphism, but various approaches suggest diversity is high, and may vary between breeds/populations. Functional consequences of variable MHC haplotype composition are discussed. Identifying unique features of the MHC in cattle will lead to new insights into evolution of the immune system.