Identification of novel bovine leukocyte antigen alleles and association of BoLA-DRB3.2*020:02:01 with resistance to Theileria orientalis infection in crossbred Kedah-Kelantan cattle: a pilot study

The bovine leukocyte antigen (BoLA) gene is a significant genetic part of the immune system and has been used as a disease marker in cattle. In this study, we detected Theileria orientalis, T. sinensis, Anaplasma marginale, Anaplasma platys, Candidatus Mycoplasma haemobos and Trypanosoma evansi by PCR amplification and sequencing of the amplicons. The allelic association of the BoLA-DRB3.2 gene with blood pathogen disease resistance and susceptibility in 87 Kedah-Kelantan x Brahman (KKB) and 38 Bali cattle was determined by Fisher's exact test and Cochran Mantel Haenszel (CMH) correction test. Sequence-based typing of the BoLA-DRB3.2 gene identified 43 alleles (27 previously reported alleles and 16 novel alleles) across the two cattle breeds. Alignment analysis of the 16 novel alleles revealed 90.7-95.8% and 85-92% nucleotide and amino acid identities, with the reference allele, BoLA-DRB3*016:01 cDNA clone NR-1. BoLA-DRB3*009:02 (25.6%) and BoLA-DRB3*036:01 (36%) were the most frequent alleles in KKB and Bali cattle, respectively. In KKB cattle, BoLA-DRB3*020:02:01 was significantly associated with resistance to T. orientalis whereas *007:01 and *009:02 were significantly associated with resistance to C. Mycoplasma haemobos. Also, DRB3*017:01 was associated with susceptibility to T. orientalis in KKB cattle. In the Bali cattle, BoLA-DRB3*015:01 was found to be a genetic marker of susceptibility to C. Mycoplasma haemobos infection. Therefore, this study identified BoLA-DRB3.2 alleles associated with resistance and susceptibility to T. orientalis infection in KKB cattle and susceptibility to C. Mycoplasma haemobos infection in Bali cattle for the first time. Therefore, this study suggests that these BoLA-DRB3 resistance alleles could be used as candidate markers for selection, whereas susceptibility alleles could be used as candidate markers for culling in the beef industry.

Toll like receptor 4 (TLR4) gene polymorphism and its association with somatic cell score and milk production traits in Indian dromedary camels

Our study aimed to explore the genetic variation in the Toll-like receptor 4 (TLR4) gene and establish its association with somatic cell score (SCS) and milk production traits in four Indian camel breeds namely Bikaneri, Kachchhi, Jaisalmeri and Mewari. TLR4 gene fragment of 573 bp spanning 5' UTR, exon-1 and partial intron-1 region was amplified and genotyped using the PCR-sequence based typing method. Only one SNP located at position C472T was identified. Genotyping revealed two alleles (C and T) and three genotypes: CC, CT and TT. The genotype frequencies for CC, CT and TT were 0.116, 0.326 and 0.558 and allele frequencies for C and T alleles were 0.279 and 0.721, respectively. Association study inferred that the effect of genotype on SCS, lactation yield (LY) and peak yield (PY) was non-significant however heterozygote (CT) genotypes recorded lower SCS and higher LY and PY. It can be concluded that the TLR4 gene possesses limited genetic variation, depicting polymorphism at a single locus in Indian camel breeds with a predominance of the TT genotype. The association study indicated that heterozygote animals possess better udder health and production performance, the statistical significance of which needs to be established using a large data set.

Association between BoLA-DRB3 polymorphism and bovine leukemia virus proviral load in Vietnamese Holstein Friesian cattle

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis. Polymorphism in bovine leukocyte antigen (BoLA)-DRB3 allele can influence the host immune response to pathogens, including BLV. However, association between specific BoLA-DRB3 alleles and BLV proviral load (PVL), which is a useful index for estimating disease progression and transmission risk, in Vietnamese cattle are unknown. Here, association study of BoLA-DRB3 allele frequency between cattle with high or low PVL demonstrated BoLA-DRB3*12:01 associates with high PVL in Vietnamese Holstein-Friesian (HF) crossbred cattle. This is the first study to demonstrate that BoLA-DRB3 polymorphism confers susceptibility to BLV high PVL in HF crossbred kept in Vietnam. Our results may be useful in disease control and eradiation for BLV through genetic selection. This article is protected by copyright. All rights reserved.

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.

High throughput analysis of MHC-I and MHC-DR diversity of Brazilian cattle populations

The major histocompatibility complex (MHC) contains many genes that play key roles in initiating and regulating immune responses. This includes the polymorphic MHCI and MHCII genes that present epitopes to CD8+ and CD4+ T-cells, respectively. Consequently, the characterisation of the repertoire of MHC genes is an important component of improving our understanding of the genetic variation that determines the outcomes of immune responses. In cattle, MHC (BoLA) research has predominantly focused on Holstein-Friesian animals (as the most economically important breed globally), although the development of high-throughput approaches has allowed the BoLA-DRB3 repertoire to be studied in a greater variety of breeds. In a previous study we reported on the development of a MiSeq-based method to enable high-throughput and high-resolution analysis of bovine MHCI repertoires. Herein, we report on the expansion of this methodology to incorporate analysis of the BoLA-DRB3 and its application to analyse MHC diversity in a large cohort of cattle from Brazil (>500 animals), including representatives from the three major Bos indicus breeds present in Brazil - Guzerat, Gir and Nelore. This large-scale description of paired MHCI-DRB3 repertoires in Bos indicus cattle has identified a small number of novel DRB3 alleles, a large number of novel MHCI alleles and haplotypes, and provided novel insights into MHCI-MHCII association - further expanding our knowledge of bovine MHC diversity.

Bovine major histocompatibility complex (BoLA) heterozygote advantage against the outcome of bovine leukemia virus infection

Bovine leukemia virus (BLV) causes enzootic bovine leucosis. Host genetic heterozygosity at the major histocompatibility complex can enhance the ability to combat infectious diseases. However, heterozygote advantage is loci specific and depends on disease type. Bovine leukocyte antigen (BoLA)-DRB3 polymorphisms are related with BLV-infection outcome; however, whether BoLA-DRB3 heterozygotes have an advantage against BLV-induced lymphoma and proviral load (PVL) remains unclear. By analyzing 1567 BLV-infected individuals, we found that BoLA-DRB3 heterozygous status was significantly associated with lymphoma resistance irrespective of cattle breeds (p < 0.0001). Similarly, decreased PVL was observed in BoLA-DRB3 heterozygotes (p = 0.0407 for Holstein cows; p = 0.0889 for Japanese Black cattle). Our report provides first evidence of BoLA-DRB3 heterozygote advantage against BLV infection outcome.

Integral Use of Immunopeptidomics and Immunoinformatics for the Characterization of Antigen Presentation and Rational Identification of BoLA-DR-Presented Peptides and Epitopes

MHC peptide binding and presentation is the most selective event defining the landscape of T cell epitopes. Consequently, understanding the diversity of MHC alleles in a given population and the parameters that define the set of ligands that can be bound and presented by each of these alleles (the immunopeptidome) has an enormous impact on our capacity to predict and manipulate the potential of protein Ags to elicit functional T cell responses. Liquid chromatography-mass spectrometry analysis of MHC-eluted ligand data has proven to be a powerful technique for identifying such peptidomes, and methods integrating such data for prediction of Ag presentation have reached a high level of accuracy for both MHC class I and class II. In this study, we demonstrate how these techniques and prediction methods can be readily extended to the bovine leukocyte Ag class II DR locus (BoLA-DR). BoLA-DR binding motifs were characterized by eluted ligand data derived from bovine cell lines expressing a range of DRB3 alleles prevalent in Holstein-Friesian populations. The model generated (NetBoLAIIpan, available as a Web server at www.cbs.dtu.dk/services/NetBoLAIIpan) was shown to have unprecedented predictive power to identify known BoLA-DR-restricted CD4 epitopes. In summary, the results demonstrate the power of an integrated approach combining advanced mass spectrometry peptidomics with immunoinformatics for characterization of the BoLA-DR Ag presentation system and provide a prediction tool that can be used to assist in rational evaluation and selection of bovine CD4 T cell epitopes.

Genetic variation in Japanese Holstein cattle for EBL development

BACKGROUND

Infection with bovine leukemia virus (BLV), the causative agent for enzootic bovine leukosis (EBL), is increasing in dairy farms of Japan. The tendency of tumor development following BLV infection in certain cow families and bull lines has previously been described. We therefore hypothesized the existence of a genetic component which differentiates cattle susceptibility to the disease.

RESULTS

We analyzed routinely collected large-scale data including postmortem inspection data, which were combined with pedigree information and epidemiological data of BLV infection. A total of 6,022 postmortem inspection records of Holstein cattle, raised on 226 farms served by a regional abattoir over 10 years from 2004 to 2015, were analyzed for associations between sire information and EBL development. We then identified statistically the relative susceptibility to EBL development for the progeny of specific sires and paternal grandsires (PGSs). The heritability of EBL development was calculated as 0.19. Similarly, proviral loads (PVLs) of progeny from identified sires and PGSs were analyzed, but no significant differences were found.

CONCLUSIONS

These observations suggest that because EBL development in our Holstein population is, at least in part, influenced by genetic factors independent of PVL levels, genetic improvement for lower incidence of EBL development in cattle notwithstanding BLV infection is possible.

Characterization of bovine MHC DRB3 diversity in global cattle breeds, with a focus on cattle in Myanmar

Background

Myanmar cattle populations predominantly consist of native cattle breeds (Pyer Sein and Shwe), characterized by their geographical location and coat color, and the Holstein-Friesian crossbreed, which is highly adapted to the harsh tropical climates of this region. Here, we analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus that has been linked to the immune response, in Myanmar cattle populations.

Methods

Blood samples (n = 294) were taken from two native breeds (Pyer Sein, n = 163 and Shwe Ni, n = 69) and a cattle crossbreed (Holstein-Friesian, n = 62) distributed across six regions of Myanmar (Bago, n = 38; Sagaing, n = 77; Mandalay, n = 46; Magway, n = 46; Kayin, n = 43; Yangon, n = 44). In addition, a database that included 2428 BoLA-DRB3 genotypes from European (Angus, Hereford, Holstein, Shorthorn, Overo Negro, Overo Colorado, and Jersey), Zebuine (Nellore, Brahman and Gir), Asian Native from Japan and Philippine and Latin-American Creole breeds was also included. Furthermore, the information from the IPD–MHC database was also used in the present analysis. DNA was genotyped using the sequence-based typing method. DNA electropherograms were analyzed using the Assign 400ATF software.

Results

We detected 71 distinct alleles, including three new variants for the BoLA-DRB3 gene. Venn analysis showed that 11 of these alleles were only detected in Myanmar native breeds and 26 were only shared with Asian native and/or Zebu groups. The number of alleles ranged from 33 in Holstein-Friesians to 58 in Pyer Seins, and the observed versus unbiased expected heterozygosity were higher than 0.84 in all the three the populations analyzed. The F_ST analysis showed a low level of genetic differentiation between the two Myanmar native breeds (F_ST = 0.003), and between these native breeds and the Holstein-Friesians (F_ST <  0.021). The average F_ST value for all the Myanmar Holstein-Friesian crossbred and Myanmar native populations was 0.0136 and 0.0121, respectively. Principal component analysis (PCA) and tree analysis showed that Myanmar native populations grouped in a narrow cluster that diverged clearly from the Holstein-Friesian populations. Furthermore, the BoLA-DRB3 allele frequencies suggested that while some Myanmar native populations from Bago, Mandalay and Yangon regions were more closely related to Zebu breeds (Gir and Brahman), populations from Kayin, Magway and Sagaing regions were more related to the Philippines native breeds. On the contrary, PCA showed that the Holstein-Friesian populations demonstrated a high degree of dispersion, which is likely the result of the different degrees of native admixture in these populations.

Conclusion

This study is the first to report the genetic diversity of the BoLA-DRB3 gene in two native breeds and one exotic cattle crossbreed from Myanmar. The results obtained contribute to our understanding of the genetic diversity and distribution of BoLA-DRB3 gene alleles in Myanmar, and increases our knowledge of the worldwide variability of cattle BoLA-DRB3 genes, an important locus for immune response and protection against pathogens.

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.

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.

A new method for typing bovine major histocompatibility complex class II DRB3 alleles by combining two established PCR sequence-based techniques

Recently, two polymerase chain reaction sequence-based typing (PCR-SBT) methods were reported for the genotyping of the bovine leukocyte antigen (BoLA)-DRB3. One technique is a single PCR-SBT (sPCR-SBT) method that generates heterozygous sequences that are subsequently analyzed by the haplofinder program, while the other technique is a nested PCR-SBT (nPCR-SBT) method that allows the analysis of heterozygous sequences using the assign 400ATF software. In this study, these techniques were compared and then integrated to produce an improved genotyping method. The primer set used for sPCR-SBT was more accurate than those used for nPCR-SBT. Combining sPCR-SBT with the assign 400ATF software previously reported for nPCR-SBT enables rapid and accurate genotyping of a large number of DNA samples.

Distributions of HLA-A and -B alleles and haplotypes in the Yi ethnic minority of Yunnan, China: relationship to other populations

Objective: To investigate the distributions of human leukocyte antigen (HLA)-A and -B alleles and HLA-A-B haplotypes in the Yi ethnic minority of the Yunnan Province, situated in southwestern China. Methods: DNA typing for HLA-A and -B loci was performed using the polymerase chain reaction-sequence-based typing (PCR-SBT) method on 114 randomly selected healthy individuals of the Yi population. The allelic frequencies of HLA-A and -B loci were calculated by direct counting and HLA-A-B haplotypes were estimated using the expectation maximization algorithm. Results: A total of 17 HLA-A and 38 HLA-B alleles were found in the Yi population. The most frequent alleles were A*2402 (32.46%), A*1101 (26.32%), and A*0203 (10.09%) at the HLA-A locus and B*4601 (12.28%), B*1525 (10.09%), B*4001 (8.77%), and B*3802 (7.89%) at the HLA-B locus. The predominant HLA-A-B haplotypes were A*2402-B*1525 (7.86%) and A*0203-B*3802 (5.64%), followed by A*1101-B*4001 (4.69%). Phylogenetic analysis indicates that the Yi population in the Honghe, Yunnan Province of China basically belongs to groups of southeastern Asian origin, but shares some characteristics with northeastern Asian groups. Conclusion: The present study may add to the understanding of HLA polymorphism in the Yi ethnic group that was poorly defined previously, and provide useful information for bone marrow transplantation, anthropological research, and forensic sciences as well as for disease-association studies.