MHC class II DR classification based on antigen-binding groove natural selection

Transcriptional response to an alternative diet on liver, muscle, and rumen of beef cattle

Feed cost represents a major economic determinant within cattle production, amounting to an estimated 75% of the total variable costs. Consequently, comprehensive approaches such as optimizing feed utilization through alternative feed sources, alongside the selection of feed-efficient animals, are of great significance. Here, we investigate the effect of two diets, traditional corn-grain fed and alternative by-product based, on 14 phenotypes related to feed, methane emission and production efficiency and on multi-tissue transcriptomics data from liver, muscle, and rumen wall, derived from 52 Nellore bulls, 26 on each diet. To this end, diets were contrasted at the level of phenotype, gene expression, and gene-phenotype network connectivity. As regards the phenotypic level, at a P value < 0.05, significant differences were found in favour of the alternative diet for average daily weight gain at finishing, dry matter intake at finishing, methane emission, carcass yield and subcutaneous fat thickness at the rib-eye muscle area. In terms of the transcriptional level of the 14,776 genes expressed across the examined tissues, we found 487, 484, and 499 genes differentially expressed due to diet in liver, muscle, and rumen, respectively (P value < 0.01). To explore differentially connected phenotypes across both diet-based networks, we focused on the phenotypes with the largest change in average number of connections within diets and tissues, namely methane emission and carcass yield, highlighting, in particular, gene expression changes involving SREBF2, and revealing the largest differential connectivity in rumen and muscle, respectively. Similarly, from examination of differentially connected genes across diets, the top-ranked most differentially connected regulators within each tissue were MEOX1, PTTG1, and BASP1 in liver, muscle, and rumen, respectively. Changes in gene co-expression patterns suggest activation or suppression of specific biological processes and pathways in response to dietary interventions, consequently impacting the phenotype. The identification of genes that respond differently to diets and their associated phenotypic effects serves as a crucial stepping stone for further investigations, aiming to build upon our discoveries. Ultimately, such advancements hold the promise of improving animal welfare, productivity, and sustainability in livestock farming.

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

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

Cross-population genetic analysis revealed genetic variation and selection in the Ovar-DRB1 gene of Indian sheep breeds

In sheep, MHC variability is studied widely to explore disease association. The aim of the current study was to explore the genetic diversity of Ovar-DRB diversity across sheep breeds of India. Here, Ovar-DRB1 locus was studied across 20 sheep breeds. DRB1 was amplified (301 bp) and sequenced using a PCR-sequence-based typing approach. Results revealed a high degree of heterozygosity across breeds (mean: 73.99%). Overall mean distance for DRB1 was highest in Sangamneri (0.18) and lowest in Madgyal sheep (0.10). There was a higher rate of transition, across breeds. Further, 39 alleles were isolated in different breeds, out of which 10 were new. To allow easy access and use of the immune-polymorphic database, an online database management system was launched (http://www.mhcdbms.in/). Nucleotide content across breeds for the DRB1 region revealed the richness of GC content (59.26%). Wu-Kabat index revealed vast genetic variation across peptide binding sites (PBS) of DRB1. Residues 6, 66, 69, 52, and 81, were polymorphic showing utility for antigen presentation. All breeds were under positive selection for DRB1 locus (dN > dS). Study revealed the importance of DRB locus diversity for beta chain specifically at PBS across sheep breeds of the Indian subcontinent and presented evidence of positive selection for DRB owing to its evolutionary significance.

Diversity of the BoLA-DRB3 gene in cattle breeds from tropical and subtropical regions of Argentina

Bovine leukocyte antigens (BoLA) have been widely studied because of their primary function in the recognition of pathogens by the immune system. To date, however, the characterization of the BoLA-DRB3 gene in Latin American Zebu and mixed zebuine breeds is scarce. By a sequence-based typing method, here we sequenced exon 2 of BoLA class II DRB3 gene in 264 animals from the five most commonly used breeds in northern Argentina (Creole, Brahman, Braford, Brangus, and Nellore).The Bos taurus, Bos indicus, and mixed breeds analyzed here contained 61 previously reported alleles. Genetic diversity was high at both allelic and nucleotide sequence levels, particularly in the mixed breeds Braford and Brangus. In contrast to previous reports on DRB3 diversity, no evidence of balancing selection was found in our data. Differentiation among breeds was highly significant, as shown by FST (F_ST = 0.052, P < 0.001) and cluster analyses. In accordance with historical origin of the breeds, UPGMA trees and metric multidimensional scaling (MDS) analyses showed that Creole is distantly related to the other zebuine breeds. Among them, Brahman, Braford, and Brangus exhibited the closest affiliations. Despite the overall differentiation of the breeds, analysis of the peptide binding regions at the aminoacid level revealed that the key aminoacids involved in peptide recognition are greatly conserved suggesting little influence of domestication and breeding in functional MHC variability. In sum, this is the first report of BoLA-DRB3 diversity in pure and mixed Bos indicus cattle breeds from Argentina. Knowledge of BoLA-DRB3 variability in breeds adapted to tropical and subtropical environments contributes not only to the characterization of MHC diversity but also to the design of peptide-based vaccines.

BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment.

Association Between BoLA-DRB3.2 Polymorphism and Bovine Papillomavirus Infection for Bladder Tumor Risk in Podolica Cattle

Blood samples from 260 unrelated cattle (132 animals affected by papillomavirus-associated bladder tumors and 128 healthy) were genotyped using the classic polymerase chain reaction/restriction fragment length polymorphism method to screen MHC class II bovine leukocyte antigen-DRB3. 2 polymorphism. The DRB3^*22 allele was significantly (p ≤ 0.01) detected in healthy cattle, thus appearing to have a negative association (protective effect) with virus infection of the urinary bladder known to represent a bladder tumor risk for cattle living free at pasture. Considering the two sequence alleles identified in animals carrying DRB3^*22, DRB3^*011:01 allele from samples of animals harboring the unexpressed bovine papillomaviruses (BPV)-2 E5 gene was characterized by amino acid residues believed to have a protective effect against BPV infection such as arginine at position 71 (R⁷¹) in pocket 4, histidine at position 11 (H¹¹) in pocket 6, and both glutamine at position 9 (Q⁹) and serine at position 57 (S⁵⁷) in pocket 9 of the antigen-binding groove. The DRB3^*011:02v allele from affected animals was characterized by amino acids believed to be susceptibility residues such as lysine (K⁷¹), tyrosine (Y¹¹), glutamic acid (E⁹), and aspartic acid (D⁵⁷) in these pockets. These results suggest that animals harboring the DRB3^*011:01 allele may have a lower risk of BPV infection and, consequently, a reduced risk of bladder tumors.

Association of Bovine Leukemia Virus-Induced Lymphoma with BoLA-DRB3 Polymorphisms at DNA, Amino Acid, and Binding Pocket Property Levels

Bovine leukemia virus (BLV) causes enzootic bovine leucosis, a malignant B-cell lymphoma in cattle. The DNA sequence polymorphisms of bovine leukocyte antigen (BoLA)-DRB3 have exhibited a correlation with BLV-induced lymphoma in Holstein cows. However, the association may vary between different cattle breeds. Furthermore, little is known about the relationship between BLV-induced lymphoma and DRB3 at the amino acid and structural diversity levels. Here, we comprehensively analyzed the correlation between BLV-induced lymphoma and DRB3 at DNA, amino acid, and binding pocket property levels, using 106 BLV-infected asymptomatic and 227 BLV-induced lymphoma Japanese black cattle samples. DRB3*011:01 was identified as a resistance allele, whereas DRB3*005:02 and DRB3*016:01 were susceptibility alleles. Amino acid association studies showed that positions 9, 11, 13, 26, 30, 47, 57, 70, 71, 74, 78, and 86 were associated with lymphoma susceptibility. Structure and electrostatic charge modeling further indicated that binding pocket 9 of resistance DRB3 was positively charged. In contrast, alleles susceptible to lymphoma were neutrally charged. Altogether, this is the first association study of BoLA-DRB3 polymorphisms with BLV-induced lymphoma in Japanese black cattle. In addition, our results further contribute to understanding the mechanisms regarding how BoLA-DRB3 polymorphisms mediate susceptibility to BLV-induced lymphoma.

Polymorphism of duck MHC class molecules

Major histocompatibility complex class I (MHC I) molecules are critically involved in defense against pathogens, and their high polymorphism is advantageous to a range of immune responses, especially in duck displaying biased expression of one MHC I gene. Here, we examined MHC I polymorphism in two duck (Anas platyrhynchos) breeds from China: Shaoxing (SX) and Jinding (JD). Twenty-seven unique UAA alleles identified from the MHC I genes of these breeds were analyzed concerning amino acid composition, homology, and phylogenetic relationships. Based on amino acid sequence homology, allelic groups of Anas platyrhynchos MHC I (Anpl-MHC I) were established and their distribution was analyzed. Then, highly variable sites (HVSs) in peptide-binding domains (PBD) were estimated and located in the three-dimensional structure of Anpl-MHC I. The UAA alleles identified showed high polymorphism, based on full-length sequence homology. By adding the alleles found here to known Anpl-MHC I genes from domestic ducks, they could be divided into 17 groups and four novel groups were revealed for SX and JD ducks. The UAA alleles of the two breeds were not divergent from the MHC I of other duck breeds, and HVSs were mostly located in the peptide-binding groove (PBG), suggesting that they might determine peptide-binding characteristics and subsequently influence peptide presentation and recognition. The results from the present study enrich Anpl-MHC I polymorphism data and clarify the distribution of alleles with different peptide-binding specificities, which might also accelerate effective vaccine development and help control various infections in ducks.

Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations

Background

Bovine leukocyte antigens (BoLAs) are used extensively as markers of disease and immunological traits in cattle. However, until now, characterization of BoLA gene polymorphisms in Zebu breeds using high resolution typing methods has been poor. Here, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 421 cattle (116 Bolivian Nellore, 110 Bolivian Gir, and 195 Peruvian Nellore-Brahman). Data from 1416 Taurine and Zebu samples were also included in the analysis.

Results

We identified 46 previously reported alleles and no novel variants. Of note, 1/3 of the alleles were detected only in Zebu cattle. Comparison of the degree of genetic variability at the population and sequence levels with genetic distance in the three above mentioned breeds and nine previously reported breeds revealed that Zebu breeds had a gene diversity score higher than 0.86, a nucleotide diversity score higher than 0.06, and a mean number of pairwise differences greater than 16, being similar to those estimated for other cattle breeds. A neutrality test revealed that only Nellore-Brahman cattle showed the even gene frequency distribution expected under a balanced selection scenario. The F_ST index and the exact G test showed significant differences across all cattle populations (F_ST = 0.057; p <  0.001). Neighbor-joining trees and principal component analysis identified two major clusters: one comprising mainly European Taurine breeds and a second comprising Zebu breeds. This is consistent with the historical and geographical origin of these breeds. Some of these differences may be explained by variation of amino acid motifs at antigen-binding sites.

Conclusions

The results presented herein show that the historical divergence between Taurine and Zebu cattle breeds is a result of origin, selection, and adaptation events, which would explain the observed differences in BoLA-DRB3 gene diversity between the two major bovine types. This allelic information will be important for investigating the relationship between the major histocompatibility complex and disease, and contribute to an ongoing effort to catalog bovine MHC allele frequencies according to breed and location.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0618-7) contains supplementary material, which is available to authorized users.

The great diversity of major histocompatibility complex class II genes in Philippine native cattle

Bovine leukocyte antigens (BoLA) are extensively used as markers for bovine disease and immunological traits. However, none of the BoLA genes in Southeast Asian breeds have been characterized by polymerase chain reaction (PCR)-sequence-based typing (SBT). Therefore, we sequenced exon 2 of the BoLA class II DRB3 gene from 1120 individual cows belonging to the Holstein, Sahiwal, Simbrah, Jersey, Brahman, and Philippine native breeds using PCR-SBT. Several cross-breeds were also examined. BoLA-DRB3 PCR-SBT identified 78 previously reported alleles and five novel alleles. The number of BoLA-DRB3 alleles identified in each breed from the Philippines was higher (71 in Philippine native cattle, 58 in Brahman, 46 in Holstein × Sahiwal, and 57 in Philippine native × Brahman) than that identified in breeds from other countries (e.g., 23 alleles in Japanese Black and 35 in Bolivian Yacumeño cattle). A phylogenetic tree based on the D_A distance calculated from the BoLA-DRB3 allele frequency showed that Philippine native cattle from different Philippine islands are closely related, and all of them are closely similar to Philippine Brahman cattle but not to native Japanese and Latin American breeds. Furthermore, the BoLA-DRB3 allele frequency in Philippine native cattle from Luzon Island, located in the Northern Philippines was different from that in cattle from Iloilo, Bohol, and Leyte Islands, which are located in the Southern Philippines. Therefore, we conclude that Philippine native cattle can be divided into two populations, North and South areas. Moreover, a neutrality test revealed that Philippine native cattle from Leyte showed significantly greater genetic diversity, which may be maintained by balancing selection. This study shows that Asian breeds have high levels of BoLA-DRB3 polymorphism. This finding, especially the identification of five novel BoLA-DRB3 alleles, will be helpful for future SBT studies of BoLA-DRB3 alleles in East Asian cattle.

Haplotype determination of the upstream regulatory region and the second exon of the BoLA-DRB3 gene in Holstein cattle

Polymorphisms of the BoLA-DRB3 gene are located primarily in the second exon [antigen binding site (ABS)] and, to a lesser extent, in the upstream regulatory region (URR). It can be hypothesised that exon 2 and the URR are under different types of natural selection. The aim of this work was to determine the URR-exon 2 haplotypes; 34 Holstein samples were genotyped by direct sequencing. A total of 7 URR alleles and 23 exon 2 alleles were detected, and 3 of the URR alleles were novel. Our results may suggest that no relationship exists between the URR and exon 2 of the BoLA-DRB3 gene (linkage disequilibrium P value > 0.05), most likely due to recombination over time. Our results also suggest that both regions of class II genes may be included in the development of new genotyping methods based on next-generation DNA sequencing technologies.

Association of BoLA DRB3 alleles with variability in immune response among the crossbred cattle vaccinated for foot-and-mouth disease (FMD)

Polymorphism of bovine leukocyte antigen (BoLA) DRB3 gene is being intensively investigated for potential association with economically important diseases of cattle. Accordingly, we investigated the association of DRB3 Exon 2 polymorphism as evidenced by the variation in the binding pockets with variability in immune response to inactivated trivalent (O, A and Asia1) foot and mouth disease virus (FMDV) vaccine in a closed population of crossbred cattle. Antibody titer of ≥ 1.8 was set as the cut off value to distinguish the protected (≥ 1.8) and unprotected (<1.8) animals. Eleven different alleles of over 3% frequency were detected in the population. We found that DRB3 alleles 0201, 0801 and 1501 always ranked high for protective immune response whereas alleles 0701, 1103 and 1101 consistently ranked low for unprotected immune response for all the three serotypes. Rank correlation of DRB3 alleles among the three serotypes was positive, high in magnitude and statistically significant (P<0.05). Logistic regression analysis revealed that odds of protection from the vaccine were highest for all the three serotypes if allele (∗)1501 was present and strengthened the results of allele ranking. Predicted amino acid substitution in the peptide binding pockets revealed that all the important sites had high Wu-Kabat index. Similarly, specific residues in pockets were crucial for immune response to FMD vaccine. There were specific substitutions in un-protected alleles such as absence of acidic amino acids substituted by basic amino acid at β71, presence of non-polar cysteine or basic histidine at β30 and presence of polar tyrosine at β37. From the observations, we hypothesize that the substitutions lead to unique conformational changes in the protein products of the studied alleles that would associate with the protective or unprotective antibody response to FMDV vaccine. The knowledge has potential implications in future selection programs if integrated with the complete BoLA haplotype details and production traits of the herd.

Minor histocompatibility antigen UTY as target for graft-versus-leukemia and graft-versus-haematopoiesis in the canine model

Male patients with female-stem-cell donors have better prognosis compared to female-to-male combinations due to Y-encoded minor histocompatibility antigens recognized by female-alloimmune-effector lymphocytes in the context of a graft-versus-leukemia (GvL) effect. We provide data in a dog-model that the minor histocompatibility antigen UTY might be a promising target to further improve GvL-immune reactions after allogeneic-stem-cell transplantations. Female-canine-UTY-specific T cells (CTLs) were stimulated in vitro using autologous-DCs loaded with three HLA-A2-restricted-UTY-derived peptides (3-fold-expansion), and specific T cell responses were determined in 3/6 female dogs. CTLs specifically recognized/lysed autologous-female-peptide-loaded DCs, but not naïve-autologous-female DCs and monocytes. They mainly recognized bone-marrow (BM) and to a lower extent DCs, monocytes, PBMCs and B-cells from DLA-identical-male littermates and peptide-loaded T2-cells in an MHC-I-restricted manner. A UTY-/male-specific reactivity was also obtained in vivo after stimulation of a female dog with DLA-identical-male PBMCs. In summary, we demonstrated natural UTY processing and presentation in dogs. We showed that female-dog CTLs were specifically stimulated by HLA-A2-restricted-UTY peptides, thereby enabling recognition of DLA-identical-male cells, mainly BM cells. These observations suggest UTY as a promising candidate-antigen to improve GvL-reactions in the course of immunotherapy.

Character of chicken polymorphic major histocompatibility complex class II alleles of 3 Chinese local breeds

To better understand the major histocompatibility complex (MHC) genetic character of domestic birds, we sequenced and analyzed chicken MHC II (B-L) genes of 3 local chicken breeds, derived from 3 separate areas in China. We amplified cDNA sequences from 105 individuals, accounting for 35 alleles. Some of the same B-LB alleles with a high frequency were found in all samples. The putative B-L α-chain had few polymorphic sites, whereas the B-L β-chain had several polymorphic sites. Most of the mutation positions were located in the B-LB β1 domain encoded by exon 2, especially in the peptide-binding region. This indicated that the highly polymorphic peptide-binding region could potentiate binding diverse antigen epitopes. The comparison of 3-D molecule structures of chicken B-L and human HLA-DR1 revealed a distinctly structural similarity, but the chicken B-L molecule had more polymorphic sites than the human HLA-DR1 molecule, which presumably might be a mechanism to compensate for responding to a wider array of pathogens due to fewer loci for chicken. Moreover, some conserved sites in human and chicken MHC class II molecules reflected their common ancestry and similar functions. These results suggest that the chicken B-L gene showed more polymorphic sites and distinctly dominant trans-breed alleles, potentially to adapt to pathogens.

Identification and diversity of bovine major histocompatibility complex class II haplotypes in Japanese Black and Holstein cattle in Japan

Bovine leukocyte antigen (BoLA), the major histocompatibility complex of cattle, is one of the most polymorphic gene clusters. We genotyped a population of 109 Japanese Black and 39 Holstein cattle to analyze their BoLA class II haplotypes, BoLA-DRB3 locus, 5 BoLA-DQA loci, and 5 BoLA-DQB loci. We identified 26 previously reported DRB3 alleles, 22 previously reported and 3 new DQA alleles, and 24 previously reported and 6 new DQB alleles. A dendrogram was constructed based on the predicted amino acid sequences of the α1 or β1 domains encoded by BoLA-DQA or -DQB alleles, which revealed that DQA alleles were clustered into 5 loci, whereas DQB alleles could not be clearly assigned to specific DQB loci. The BoLA-DRB3-DQA-DQB haplotypes were sorted by sequential analytical processes, and 42 distinct haplotypes, including 11 previously published haplotypes and 31 novel haplotypes, were defined. Strong linkage disequilibrium was present in the BoLA genes. We also compared DRB3-DQA1 haplotype frequencies between 507 Japanese Black and 143 Holstein cattle. Thirty-nine DRB3-DQA1 haplotypes were identified, including 29 haplotypes from Japanese Black and 22 haplotypes from Holstein cattle. The majority of the haplotypes could be identified in both breeds, although several haplotypes were identified in only a single breed. This is the first report presenting a detailed study of the BoLA class II haplotype in Japanese Black and Holstein cattle in Japan.