- Research article
- Open Access
Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations
BMC Geneticsvolume 19, Article number: 33 (2018)
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.
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 FST index and the exact G test showed significant differences across all cattle populations (FST = 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.
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.
Domestic bovines comprise two major breed types: humpless (Bos taurus) and humped (Bos indicus or Bos taurus indicus). Based on archeological and molecular evidence, it is estimated that the Taurine and Zebu lineages divergence occurred from a common ancestor between 330,000 and 2 million years ago, respectively [1,2,3,4]. Since then, the two types of cattle have accumulated many genetic variations resulting in a high genetic distance; this has contributed to highly differentiated phenotypes and adaptation. Currently, Zebu breeds are found worldwide in tropical and subtropical regions; indeed, beef production in South America’s warm regions is based mainly on Zebu cattle (Nellore, Brahman, and Gir breeds).
The bovine MHC, named bovine leukocyte antigen (BoLA), comprises three DRB class II copies, with the BoLA-DRB3 gene being the most highly expressed and polymorphic . Polymorphic sites within the BoLA-DRB3 gene are located mainly at peptide-binding sites; these polymorphisms are mostly maintained by balancing selection . Several lines of evidence show that polymorphisms influence both the magnitude and epitope specificity of antigen-specific T cell responses to infectious diseases . In this sense, the BoLA-DRB3 gene is associated with differences in susceptibility to infectious disease (e.g., bovine leukosis virus-induced lymphocytosis, mastitis, and dermatophilosis), different immunological and production traits (e.g., milk yield), and different vaccine responses (e.g., to foot-and-mouth disease and Theileria parva.) [8, 9].
Since the first pioneering studies based on serotype analysis [10, 11], others have used different genotyping techniques to examine the genetic diversity of the BoLA-DRB3 gene in a large number of bovine breeds. Currently, polymerase chain reaction (PCR)-sequence-based typing (SBT) [12,13,14,15,16,17,18,19,20,21,22] and next generation sequencing (NGS)  are the most powerful methods. More than 25 years’ work has identified 107 BoLA-DRB3 alleles in cattle; these are listed in the Immuno Polymorphism Database (IPD)–MHC database (https://www.ebi.ac.uk/ipd/mhc/group/BoLA). Forty-two of them were reported in Bos indicus.
SBT has better discriminatory power and precision than PCR-RFLP to identify diversity of BoLA-DRB3 alleles in different cattle breeds . However, BoLA-DRB3 diversity has been characterized at the DNA sequence level only in 12 out of over 800 breeds of cattle recognized worldwide [12,13,14,15,16,17,18,19,20,21,22,23]. Studied breeds include European taurine (e.g., Hereford, Angus, Shorthorn, Holstein, Jersey, Overo Colorado, and Overo Negro), Asian taurine (Japanese Black, Philippine Native, and Korean Native), American Creole (e.g., Yacumeño and Hartón del Valle), and different types of cross-breed (e.g., Holstein × Sahiwal, Philippine Native × Brahman, and Philippine Native × Holstein).
Until now, several private Zebu BoLA-DRB3 alleles have been reported by authors using indirect techniques such as PCR-RFLP, followed by cloning and sequencing [24,25,26,27]. These studies focused mainly on screening and analysis of a few animals from some Zebu breeds (e.g., Brahman, Dairy Gir, Boran, and Fulani); however, no study has used SBT or NGS assays to examine genetic diversity of BoLA-DRB3 in large samples. With this in mind, the aim of the present study was to examine the genetic diversity of the BoLA-DRB3 gene in three Zebu breeds from Bolivia and Peru, and to compare the results with those reported for Taurine breeds. The results may contribute to an ongoing effort to catalog bovine MHC allele frequency according to breed and location.
Sample populations and genomic DNA extraction
Blood samples were obtained from 421 Zebu cattle belonging to Bolivian Nellore (N = 116), Bolivian Gir (N = 110), and Peruvian Nellore-Brahman (N = 195) breeds. These breeds were selected because they are bred most widely in South America. In addition, previously reported data from 1416 animals from Taurine and Zebu breeds from Japan, the Philippines, Colombia, Chile, and Bolivia were included for comparison [12,13,14,15,16,17,18,19,20,21,22]. Detailed information about the breeds/populations analyzed is provided in Table 1. Genomic DNA was extracted from 40 μl of whole blood spotted on FTA Elute cards (Whatman, Tokyo, Japan), according to the manufacturer’s instructions.
Cattle were handled by veterinarians from local farms, RIKEN, and Universidad Autónoma Gabriel René Moreno, in strict accordance with good animal practice following the Universidad Austral de Chile Institutional guidelines. Veterinarian Collaborator explained about the sample collection to farmer in verbal way, and their future use in the research project. This study was approved by Committees on the Ethics of Animals for Research at the National from University of LA PLATA (Argentine: Certificate date May 26th, 2014), Universidad Autónoma Gabriel Moreno (Bolivia, and Universidad Austral de Chile (Certificate No. 153–2014).
BoLA-DRB3 typing using PCR-sequence based typing (SBT)
BoLA-DRB3 alleles were genotyped using PCR-SBT. Briefly, exon 2 of BoLA-DRB3 was amplified by PCR as described by Takeshima et al. , using primers DRB3FRW and DRB3REV . The PCR fragments were purified using an ExoSAP-IT PCR Product Purification Kit (USB Corp., Cleveland, OH) and sequenced using the ABI PRISM BigDye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Foster City, CA). Raw sequence data were analyzed using Assign 400ATF ver. 184.108.40.206 software (Conexio Genomics, Fremantle, Australia).
Measures of genetic variability
Allele frequencies and the number of alleles (Na) were obtained by direct counting, and 95% confidence intervals (CI) for allele frequencies were calculated using the binomial distribution implemented in R with the binom.confint function (http://cran.r-project.org/web/packages/binom/). The distribution of alleles across major groups of breeds was represented by a Venn plot using the R package ‘VennDiagram’. The observed (Ho) and expected unbiased (He) heterozygosity for the BoLA-DRB3 locus were estimated as described previously  using ARLEQUIN 3.5 software for population genetic analyses . Potential deviations from Hardy-Weinberg equation (HWE) were estimated for each breed by FIS statistics  using the exact test included in GENEPOP 4.0 . The Ewens–Watterson–Slatkin exact test for neutrality was estimated as described by Slatkin  and implemented in ARLEQUIN 3.5 software.
Genetic structure and differentiation
Genetic structure and genetic differentiation among breeds were assessed using Wright’s F statistics (FST) and the previously described variance-based method . These parameters were estimated using ARLEQUIN 3.5 and GENEPOP 4.0. FST values were represented graphically using the pairFstMatrix.r function in the R statistical environment.
Principal component analysis (PCA)
To condense the genetic variation at the BoLA-DRB3 locus, allele frequencies were used to perform principal component analysis (PCA) as described previously  using PAST software . Nei’s standard genetic distance (Ds)  and Nei’s unbiased genetic distance (D A )  were calculated from allele frequencies, and cluster analysis was performed using the UPGMA  and the neighbor-joining (NJ)  algorithms. For these analyses, only pure breeds and a Nellore-Brahman mixed population were included due to limitations inherent in the dendrograms. CIs for the groupings were estimated by bootstrap re-sampling of data using 1000 replications. Genetic distances and trees were calculated using POPULATIONS 1.2.28 software , and trees were constructed using TREEVIEW .
Genetic diversity at the nucleotide and amino acid sequence levels
Nucleotide diversity (π) and pairwise comparison of nucleotide substitutions between alleles were calculated using ARLEQUIN 3.5. The mean number of nonsynonymous (dN) and synonymous (dS) nucleotide substitutions per site was estimated for each pair using the Jukes–Cantor’s formula as described by Nei and Gojobori . These parameters were estimated using MEGA 5 .
Distribution of BoLA-DRB3 alleles among analyzed zebu cattle breeds
PCR-SBT genotyping allowed to identify 46 previously reported BoLA-DRB3 alleles, while no new variants were observed (Table 2) and all genotyping result were shown in additional table (Additional file 1). The Na ranged from 19 BoLA-DRB3 alleles in Gir to 33 variants in Nellore-Brahman (Table 2). Of these, six alleles reported in the IPD-MHC public database (https://www.ebi.ac.uk/ipd/mhc/) were from Bos indicus and 26 alleles were from B. taurus; 14 alleles were shared among these two major bovine types. Furthermore, a Venn diagram was constructed using data obtained in this study and from previous reports [12,13,14,15,16,17,18,19,20,21,22]; these data, that included 14 breeds and cross-breeds, were then used to evaluate the distribution of the BoLA-DRB3 allele (Fig. 1.). To achieve this, data were grouped in terms of geographical origin, i.e., British, European Continental, American Creole, and Zebu cattle breeds. The plot clearly shows that the Zebu group, that comprised Bolivian Gir, Bolivian Nellore, Peruvian Nellore-Brahman and Philippine Brahman breeds, harbors 21 out of the 90 alleles identified in the four cattle groups (Fig. 1); these 21 alleles represent about 30% of the 66 alleles detected in the Zebu group.
The Na within each Zebu population with frequencies > 5% varied; from five in Bolivian Nellore and Peruvian Nellore-Brahman to ten in Bolivian Gir. Some of the high-frequency (> 5%) alleles were present in at least two out of three Zebu populations; these alleles were BoLA-DRB3*2801, *2201, *3101, and *3601 (Table 2). These common alleles account for 65.95% of all alleles sampled in Bolivian Nellore, 67.18% in Peruvian Nellore-Brahman, and 85.00% in Bolivian Gir (Fig. 2).
Nucleotide and amino acid diversity of BoLA-DRB3 alleles in analyzed zebu breeds
Genetic diversity was examined at the DNA level using two parameters: π and the number of pairwise differences (NPD). The results obtained were then compared with those previously reported for ten other breeds (Table 3). The π ranged from 0.068 in the Nellore × Brahman population to 0.078 in the Bolivian Gir population, while the mean NPD varied from 16.91 in the Nellore × Brahman to 19.43 in Bolivian Gir. The 13 cattle breeds analyzed by PCR-SBT varied in terms of the number and composition of BoLA-DRB3 alleles. However, with the exception of Shorthorn and Native Philippine cattle, the values observed for the π and NPD indices were similar to those estimated for other cattle breeds (Table 3).
We also calculated the average number for dn and ds in Bolivian Nellore, Peruvian Nellore × Brahman, and Bolivian Gir. As shown in Table 4, the values obtained for these populations (0.037 to 0.042 for ds and 0.103 to 0.108 for dn) when analyzing the whole of exon 2 suggested a dn/ds difference range for this genome region of 0.058 to 0.066. The dn/ds differences were more prominent when only the antigen-binding sites (ABS) were analyzed (dn/ds > 0.132; Table 3).
Gene diversity, HWE, and neutrality testing of BoLA-DRB3 in the analyzed zebu populations
Genetic diversity within the three analyzed Zebu populations was estimated according to the Na and the observed and expected heterozygosities. As mentioned above, the Na ranged from 19 in Bolivian Gir to 33 in Peruvian Nellore-Brahman, while the values of expected heterozygosity ranged from 0.86 in Nellore-Brahman to 0.92 in Bolivian Gir (Table 4). The observed heterozygosity in Peruvian Nellore-Brahman and Bolivian Gir varied from 0.76 to 0.88, respectively (Table 4). Regarding HWE, two out of three populations were in equilibrium, while Nellore × Brahman showed a significant deviation from the theoretical values (Table 4). This deviation can be explained by a significant excess of homozygous animals (FIS = 0.113, p value = 0.033). To assess balancing selection, Slatkin’s exact neutrality test was performed (Table 4). The BoLA-DRB3 gene frequency profile in Gir cattle showed an even distribution (p = 0.008) despite the finding that its genotype proportions did not deviate from the theoretical proportions. This even gene frequency distribution is consistent with the theoretical proportion expected under recent balancing selection as opposed to positive or neutral selection (p > 0.025). Conversely, the Nellore and Nellore × Brahman gene frequency distributions were compatible with neutral selection (Table 4).
Genetic structure and differentiation of BoLA-DRB3 in zebu breeds
To examine population genetic structure and differentiation levels among cattle breeds, we calculated two parameters: the FST index and the exact G test. There were significant differences in the FST parameter across all cattle populations (FST = 0.057; p value ≤0.0001); pairwise comparison ranged from 0.022 (between Philippine Brahman and Philippine Native) to 0.133 (between Peruvian Nellore × Brahman and Chilean Hereford) (Fig. 3). However, there was around 5% differentiation among Zebu populations. The exact G test for differentiation indicated that, on average, differences in allele frequency distribution among populations were highly significant (exact p value ≤0.0001). Furthermore, when all pairs of breeds were tested, significant results (p value ≤0.0001) were also obtained for all comparisons, indicating a significant level of genetic structure for BoLA-DRB3 among breeds. In conclusion, when Zebu breeds were included in the analysis, the average FST value was 1 to 5% higher than that obtained when analyzing taurine breeds alone  or a population of the same breeds , respectively.
Genetic differentiation of BoLA-DRB3 alleles in zebu and taurine breeds
To assess the genetic differentiation of BoLA-DRB3 alleles between Zebu cattle and nine breeds previously studied by SBT, we performed two types of analysis: dendrograms and PCA. First, BoLA-DRB3 allele frequencies were used to generate Nei’s DA and DS genetic distances matrices; data from Bolivian Nellore, Bolivian Gir, Peruvian Brahman × Nellore, and nine previously reported breeds (Japanese Jersey, Japanese Shorthorn, Japanese Holstein, Japanese Black, Chilean Hereford, Philippine Native, Philippine Brahman, Yacumeño, and Hartón del Valle; [12,13,14,15,16,17,18,19,20,21,22]) were used. Dendrograms based on distance were constructed using UPGMA and NJ algorithms. NJ cluster analysis based on Nei’ DA or DS genetic distances showed two major clusters: one comprising mainly European taurine breeds and Yacumeño Creole cattle and one comprising Zebu breeds, Hartón del Valle, and Philippine Native breeds (Fig. 4a). Japanese Jersey or Chilean Hereford was located outside of these clusters. UPGMA cluster analysis showed a similar trend, with Jersey and Hereford breeds located in separated clades (data not shown).
Second, we used allele frequencies of BoLA-DRB3 from the 12 breeds mentioned above to perform PCA. The results are reported in Fig. 4b, which shows the first and second principal components (PCs) for BoLA-DRB3 gene frequency. The first two components accounted cumulatively for 44.68% of the variability in the data. The first PC accounted for 28.34% of the total variance, and showed a clear differentiation pattern between Zebu (positive values) and Taurine breeds (negative values), whereas native breeds were located near the origin of the plot. The first PC was determined mainly by 26 alleles, nine with a positive value (BoLA-DRB3*0902, *2201, *2710, *2801, *3001, *3101, *3501, *3601, and *5702) and 17 with a negative value (BoLA-DRB3*0101, *0201, *0301, *0701, *0801, *1001, *1101, *1201, *1501 *1601, *1801, *2006, *2502, *3202, *4501, *14011, and *20012). It is noteworthy that almost all alleles with higher negative values for PC1 were previously reported in B. Taurus in the MHC diversity database, whereas the variants that account for PC1 positive values were reported in both B. Taurus and/or B. indicus. The second PC explained 16.34% of the total variation and shows a gradient among taurine breeds, with Hereford and Jersey located at opposite ends. The second PC was determined mainly by differences in the frequency of 20 alleles; seven with a positive value (BoLA-DRB3*1002, *1501, *1601, *1801, *3202, *14011, and *20012) and 13 with a negative value (BoLA-DRB3*0101, *0201, *0301, *0501, *0701, *0801, *1101, *1201, *2006, *2502, *3701, *4401, and *4501). Finally, the third PC accounted for 14.95% of the variance, similar to PC2; this component allowed differentiation of the B. taurus breed, but with Holstein and Jersey at opposite ends. These PCA results are in agreement with the overall clustering observed after NJ or UPGMA tree construction. Native breeds occupy an intermediate position in both analyses.
BoLA class II molecule binds various peptides derived from antigens via five antigen binding pockets named pocket 1, pocket 4, pocket 6, pocket 7 and pocket 9 . To assess whether observed differences in allelic frequency between Taurine and Zebu populations is reflected within amino acid motifs in each pocket, we analyzed protein pockets implicated in antigen-binding function by PCA. With respect to taurine breeds, Zebu breeds were clearly differentiated in component 2 of pocket 1 due to a higher frequency of GVFT and VGFT and a lower frequency of VVFT and GMFT amino acid motifs (Fig. 5). Component 1 of pockets 4 and 7 allowed differentiation of Zebu (positive values) from Taurine breeds (negative values). GLDEREY and GLDRREV amino acid motifs in pocket 4, and DYWIR, DFWFR, and EYWIR amino acid motifs in pocket 7, provided the major contribution and explain the distribution observed in PCA plots. With the exception of Philippine Brahman, Zebu breeds showed a similar trend in pocket 9, which was accounted for by QYS, QFA, and END amino acid motifs. Finally, component 1 of pocket 6 only discriminated Nellore from Nellore crossbred by means of the HH motif.
Domestic cattle were derived from the now-extinct wild aurochs (Bos primigenius) . Zooarcheological evidence revealed that this wild species was widespread throughout Eurasia and North Africa, even though, events of bovine domestication occurred only few times in the human history and in precise geographic regions (called domestication center). In this sense, different local wild auroch subspecies populations would have been domesticated in five domestication centers: three for B. taurus (southwest Asia, north Africa, and northeastern China) at about 10,000 Years B.P., and two for B. indicus (Indus Valley and southern India) at around 8000 Years B.P. Mitochondrial DNA (mtDNA) has been extensively used to reconstruct the ancestry and domestic history of cattle. Analysis of modern and ancient mtDNAs showed that auroch diversity had an uneven distribution with multiple mtDNA clades and a complex branching structure [46, 47]. Since domestication process occurred in geographically restricted areas, only a fraction of the whole wild genetic variation was retained in each domestication center (sampling effect). Furthermore, domestic cattle mtDNA phylogenies exhibited two major distinct clades, including taurine T1, T2, T3 and T4, and Zebu I1 and I2 haplogroups .
Regarding MHC Class II genes, trans-species theory supports the hypothesis that MHC alleles preceded the origin of related species and survived through the speciation process . Phylogenetic analysis of Bovidae DRB sequences also supports this theory [49, 50]. Bearing this in mind, it is expected that domestic cattle BoLA-DRB3 diversity was present in the auroch population before the divergence of major bovine types. Just as mtDNA haplogroups, BoLA-DRB3 alleles could have shown an uneven distribution in the B. primigenius population at the time of domestication, keeping only a fraction of the total diversity in each cattle lineage. Likewise, a later bottleneck having occurred during the more recent breed formation may have increased the genetic differentiation at BoLA-DRB3 alleles. On the other hand, genetic structure could have been diluted due to migration and gene introgression events occurred since domestication. Most of the detected alleles in Zebu were reported in the IPD-MHC database as Taurine, while a slightly smaller number of variants was shared among these two major bovine types, and only six were described as B. indicus. In contrast, Venn diagram, used to evaluate the distribution of the BoLA-DRB3 allele among group of breeds, clearly shows that the Zebu group, comprised of Bolivian Gir, Bolivian Nellore, Peruvian Nellore-Brahman and Philippine Brahman breeds, harbors 30% of Zebu privative alleles. This data further supports the hypothesis that BoLA-DRB3 diversity has an uneven distribution.
MHC class II gene diversity can be maintained by some form of balancing or overdominant selection [51, 52]. To assess this hypothesis, allele diversity, gene frequency profile, HWE and Slatkin’s exact neutrality tests were performed. In this sense, between 19 and 33 alleles were detected within the analyzed Zebu breeds, resulting in values of expected heterozygosity higher than 0.85. Furthermore, common alleles accounted for 65.95% of all alleles sampled in Bolivian Nellore, 67.18% in Peruvian Nellore-Brahman, and 85.00% in Bolivian Gir. Despite the fact Zebu breeds exhibit higher resistance to tropical pathogens; these results demonstrate that analyzed populations showed allele diversity and gene frequency distribution within the range reported for Taurine breeds [12, 13, 22].
Regarding HWE, none of the analyzed populations showed a significant excess of heterozygotes animals that could reflect overdominant selection, while only Gir gene frequency profile showed an even distribution consistent with the theoretical proportion expected under recent balancing selection. Even though this evolutionary force could not be detected in most cattle populations studied until now [12, 22], this even gene distribution was observed in Yacumeño and Japanese Black breeds. Furthermore, several of the analyzed breeds showed neutrality test p values ranging from 0.025 to 0.1 which, though not significant, could reflect the same bias towards balancing selection. Conversely, Nellore and Nellore × Brahman gene frequency distributions were compatible with neutral selection.
Observed diversity at allele level was also present when analyzing the nucleotide and amino acid variability. The 13 cattle breeds analyzed by PCR-SBT varied in terms of the number and composition of BoLA-DRB3 alleles. However, with the exception of Shorthorn and Native Philippine cattle, the values observed in Zebu breeds for the π and NPD indices were similar to those estimated for other cattle breeds [12, 20,21,22]. As expected, the dn/ds differences were more prominent when only the ABS were analyzed, which is in concordance with the biological function of these sites. Of note, the average number for dn and ds across the entire exon 2, as well as the ABS, in Zebu breeds was slightly smaller than the average reported for previously analyzed breeds [12, 20,21,22].
Genetic structure and differentiation of BoLA-DRB3 in Zebu breeds were assessed by calculating FST index and the exact G test. Both estimations revealed significant differences across all cattle populations. Noteworthy, when Zebu breeds were included in the analysis, the average FST value was 1 to 5% higher than that obtained when analyzing taurine breeds alone  or populations of the same Taurine breed , respectively. Moreover, dendrograms based on Nei’ DA or DS genetic distance were constructed using UPGMA and NJ algorithms. These trees showed similar relationship between breeds: two major clusters revealing topologies consistent with the historical and geographical origin. In addition, PCA results are in agreement with the overall clustering observed after NJ or UPGMA tree construction, which showed a clear divergence between Taurine and Zebu populations. Also, native breeds occupy an intermediate position in both analyses, probably due to some level of Zebu gene introgression [53, 54]. Finally, both dendrograms and PCA showed that Chilean Hereford and Japanese Jersey were the most divergent populations among B. taurus breeds. This topology could be partially explained by the ancient genetic structure present among and within B. primigenius sub-species; and the sampling effect during the domestication process. Furthermore, subsequent selection, local adaptation and breed origin could also have contributed to the observed genetic differentiation at BoLA-DRB3 gene.
PCA based on allelic frequencies of amino acid motif of pockets implicated in antigen-binding function showed a clear differentiation between Taurine and Zebu breeds, as consequence of an enrichment of particular amino acid motif in specific pocket. These differences could be consistent with the widely reported Zebu cattle high disease tolerance/susceptibility to tropical infectious disease, such as ticks and intestinal parasites [47, 55]. Further studies involving structural modeling and molecular simulation of the BoLA-DRβ protein are needed to elucidate whether these differences play a role in the function of the resultant protein.
The present study is the first to use PCR-SBT to examine the allelic distribution of the BoLA-DRB3 gene in Zebu populations. The results show that the historical divergence between Taurine and Zebu cattle breeds is due to origin, selection, and adaptation events, which explains the observed differences in BoLA-DRB3 gene diversity between these two major bovine types.
Bovine leukocyte antigens
- DA :
Nei’s unbiased genetic distance
Nei’s standard genetic distance
Wright’s F statistics
expected unbiased heterozygosity
Immuno Polymorphism Database
number of alleles
next generation sequencing
number of pairwise differences
principal component analysis
polymerase chain reaction
Polymerase chain reaction sequence-based typing
Achilli A, Olivieri A, Pellecchia M, Uboldi C, Colli L, Al-Zahery N, Accetturo M, Pala M, Kashani BH, Perego UA, et al. Mitochondrial genomes of extinct aurochs survive in domestic cattle. Curr Biol. 2008;18(4):R157–8.
Bradley D, Loftus R, Cunningham P, MacHugh D. Genetics and domestic cattle origins. Evol Anthropol. 1998;6(3):79–86.
Hiendleder S, Lewalski H, Janke A. Complete mitochondrial genomes of Bos taurus and Bos indicus provide new insights into intraspecies variation, taxonomy and domestication. Cytogenet Genome Res. 2008;120(1–2):150–6.
Loftus RT, MacHugh DE, Bradley DG, Sharp PM, Cunningham P. Evidence for two independent domestications of cattle. Proc Natl Acad Sci U S A. 1994;91(7):2757–61.
Aida Y. Characterization and expression of bovine MHC class II genes. Bull Soc Fr Jpn Sci Vet. 1995;6:17–24.
Satta Y, Ohuigin C, Takahata N, Klein J. Intensity of natural-selection at the major histocompatibility complex loci. Proc Natl Acad Sci U S A. 1994;91(15):7184–8.
Xu A, van Eijk MJ, Park C, Lewin HA. Polymorphism in BoLA-DRB3 exon 2 correlates with resistance to persistent lymphocytosis caused by bovine leukemia virus. J Immunol. 1993;151(12):6977–85.
Takeshima SN, Aida Y. Structure, function and disease susceptibility of the bovine major histocompatibility complex. Anim Sci J. 2006;77(2):138–50.
Nayeri S, Stothard P. Tissues, metabolic pathways and genes of key importance in lactating dairy cattle. Springer Science Reviews. 2016;4(2):49–77.
Davies CJ, Andersson L, Joosten I, Mariani P, Gasbarre LC, Hensen EJ. Characterization of bovine MHC class II polymorphism using three typing methods: serology, RFLP and IEF. Eur J Immunogenet. 1992;19(5):253–62.
Davies CJ, Joosten I, Andersson L, Arriens MA, Bernoco D, Bissumbhar B, Byrns G, van Eijk MJ, Kristensen B, Lewin HA, et al. Polymorphism of bovine MHC class II genes. Joint report of the fifth international bovine lymphocyte antigen (BoLA) workshop, Interlaken, Switzerland, 1 august 1992. Eur J Immunogenet. 1994;21(4):259–89.
Giovambattista G, Takeshima SN, Ripoli MV, Matsumoto Y, Franco LA, Saito H, Onuma M, Aida Y. Characterization of bovine MHC DRB3 diversity in Latin American creole cattle breeds. Gene. 2013;519(1):150–8.
Miyasaka T, Takeshima SN, Matsumoto Y, Kobayashi N, Matsuhashi T, Miyazaki Y, Tanabe Y, Ishibashi K, Sentsui H, Aida Y. The diversity of bovine MHC class II DRB3 and DQA1 alleles in different herds of Japanese black and Holstein cattle in Japan. Gene. 2011;472(1–2):42–9.
Miyasaka T, Takeshima SN, Sentsui H, Aida Y. Identification and diversity of bovine major histocompatibility complex class II haplotypes in Japanese black and Holstein cattle in Japan. J Dairy Sci. 2012;95(1):420–31.
Takeshima S, Ikegami M, Morita M, Nakai Y, Aida Y. Identification of new cattle BoLA-DRB3 alleles by sequence-based typing. Immunogenetics. 2001;53(1):74–81.
Takeshima S, Nakai Y, Ohta M, Aida Y. Short communication: characterization of DRB3 alleles in the MHC of Japanese shorthorn cattle by polymerase chain reaction-sequence-based typing. J Dairy Sci. 2002;85(6):1630–2.
Takeshima S, Saitou N, Morita M, Inoko H, Aida Y. The diversity of bovine MHC class II DRB3 genes in Japanese black, Japanese shorthorn, Jersey and Holstein cattle in Japan. Gene. 2003;316:111–8.
Takeshima SN, Matsumoto Y, Aida Y. Short communication: establishment of a new polymerase chain reaction-sequence-based typing method for genotyping cattle major histocompatibility complex class II DRB3. J Dairy Sci. 2009;92(6):2965–70.
Takeshima SN, Matsumoto Y, Miyasaka T, Arainga-Ramirez M, Saito H, Onuma M, Aida Y. A new method for typing bovine major histocompatibility complex class II DRB3 alleles by combining two established PCR sequence-based techniques. Tissue Antigens. 2011;78(3):208–13.
Takeshima SN, Miyasaka T, Polat M, Kikuya M, Matsumoto Y, Mingala CN, Villanueva MA, Salces AJ, Onuma M, Aida Y. The great diversity of major histocompatibility complex class II genes in Philippine native cattle. Meta Gene. 2014;2(0):176–90.
Takeshima SN, Miyasaka T, Matsumoto Y, Xue G, Diaz Vde L, Rogberg-Munoz A, Giovambattista G, Ortiz M, Oltra J, Kanemaki M, et al. Assessment of biodiversity in Chilean cattle using the distribution of major histocompatibility complex class II BoLA-DRB3 allele. Tissue Antigens. 2015;85(1):35–44.
Takeshima SN, Giovambattista G, Okimoto N, Matsumoto Y, Rogberg-Munoz A, Acosta TJ, Onuma M, Aida Y. Characterization of bovine MHC class II DRB3 diversity in South American Holstein cattle populations. Tissue Antigens. 2015;86(6):419–30.
Lee BY, Hur TY, Jung YH, Kim H. Identification of BoLA-DRB3.2 alleles in Korean native cattle (Hanwoo) and Holstein populations using a next generation sequencer. Anim Genet. 2012;43(4):438–41.
Gelhaus A, Schnittger L, Mehlitz D, Horstmann RD, Meyer CG. Sequence and PCR-RFLP analysis of 14 novel BoLA-DRB3 alleles. Anim Genet. 1995;26(3):147–53.
Mikko S, Andersson L. Extensive MHC class II DRB3 diversity in African and European cattle. Immunogenetics. 1995;42(5):408–13.
Maillard JC, Renard C, Chardon P, Chantal I, Bensaid A. Characterization of 18 new BoLA-DRB3 alleles. Anim Genet. 1999;30(3):200–3.
da Mota AF, Martinez ML, Coutinho LL. Genotyping BoLA-DRB3 alleles in Brazilian dairy Gir cattle (Bos indicus) by temperature-gradient gel electrophoresis (TGGE) and direct sequencing. Eur J Immunogenet. 2004;31(1):31–5.
Baxter R, Hastings N, Law A, Glass EJ. A rapid and robust sequence-based genotyping method for BoLA-DRB3 alleles in large numbers of heterozygous cattle. Anim Genet. 2008;39(5):561–3.
Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89(3):583–90.
Schneider S, Roessli D, Excoffier L. Arlequin version 2.000: a software for population genetics data analysis. In: 2.000 edn: genetics and biometry laboratory. Switzerland: University of Geneva; 2000.
Weir B, Cockerham C. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.
Rousset F. genepop'007: a complete re-implementation of the genepop software for windows and Linux. Mol Ecol Resour. 2008;8(1):103–6.
Slatkin M. A correction to the exact test based on the Ewens sampling distribution. Genet Res. 1996;68(3):259–60.
Cavalli-Sforza LL, Menozzi P, Piazza A. The history and geography of human genes. United States of America: Princeton University Press; 1994.
Hammer Ø, Harper DAT, Ryan PD. PAST: paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001;4(1):9.
Nei M. Genetic distance between populations. Am Nat. 1972;106(949):283.
Nei M, Tajima F, Tateno Y. Accuracy of estimated phylogenetic trees from molecular data. II. Gene frequency data. J Mol Evol. 1983;19(2):153–70.
Sneath PHA, Sokal RR. Numerical taxonomy: The principles and practice of numerical classification. Hardcover edn. United States: Freeman & Company Limited, W H; 1973.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.
Langella O. Populations: a population genetic software. In: vol. CNRS UPR9034., 1.2.28 edn; 1999. http://www.mybiosoftware.com/populations-1-2-32-population-genetic-software.html.
Page RD. TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996;12(4):357–8.
Nei M, Gojobori T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986;3(5):418–26.
Kumar S, Nei M, Dudley J, Tamura K. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008;9(4):299–306.
Takeshima SN, Sarai A, Saitou N, Aida Y. MHC class II DR classification based on antigen-binding groove natural selection. Biochem Biophys Res Commun. 2009;385(2):137–42.
Clutton-Brock J. A natural history of domesticated mammals. 2nd edition edn. United Kingdom: Cambridge university press; 1999.
Bradley DG, Magee DA. Genetics and domestic cattle origins. In: Documenting domestication: new genetic and archaeological paradigms. Zeder MA, Bradley DG, Emshwiller E, Smith BD. editors. 1 edn. U. S. A: University of California Press; 2006: 317–328.
Magee DA, MacHugh DE, Edwards CJ. Interrogation of modern and ancient genomes reveals the complex domestic history of cattle. Animal Frontiers. 2014;4(3):7–22.
Klein J. Origin of major histocompatibility complex polymorphism: the trans-species hypothesis. Hum Immunol. 1987;19(3):155–62.
Sena L, Schneider MP, Brenig B, Honeycutt RL, Womack JE, Skow LC. Polymorphisms in MHC-DRA and -DRB alleles of water buffalo (Bubalus bubalis) reveal different features from cattle DR alleles. Anim Genet. 2003;34(1):1–10.
Mishra SK, Niranjan SK, Banerjee B, Dubey PK, Gonge DS, Mishra BP, Kataria RS. High genetic diversity and distribution of Bubu-DQA alleles in swamp buffaloes (Bubalus bubalis carabanesis): identification of new Bubu-DQA loci and haplotypes. Immunogenetics. 2016;68(6–7):439–47.
Hedrick PW, Thomson G. Evidence for balancing selection at Hla. Genetics. 1983;104(3):449–56.
Hughes AL, Nei M. Nucleotide substitution at major histocompatibility complex class II loci: evidence for overdominant selection. Proc Natl Acad Sci U S A. 1989;86(3):958–62.
Liron JP, Peral-Garcia P, Giovambattista G. Genetic characterization of argentine and Bolivian creole cattle breeds assessed through microsatellites. J Hered. 2006;97(4):331–9.
Liron JP, Bravi CM, Mirol PM, Peral-Garcia P, Giovambattista G. African matrilineages in American creole cattle: evidence of two independent continental sources. Anim Genet. 2006;37(4):379–82.
Porto Neto LR, Jonsson NN, D'Occhio MJ, Barendse W. Molecular genetic approaches for identifying the basis of variation in resistance to tick infestation in cattle. Vet Parasitol. 2011;180(3–4):165–72.
Mikko S, Spencer M, Morris B, Stabile S, Basu T, Stormont C, Andersson L. A comparative analysis of Mhc DRB3 polymorphism in the American bison (Bison bison). J Hered. 1997;88(6):499–503.
Russell GC, Marello KL, Gallagher A, McKeever DJ, Spooner RL. Amplification and sequencing of expressed DRB second exons from Bos indicus. Immunogenetics. 1994;39(6):432–36.
Sitte K, East IJ, Lavin MF, Jazwinska EC. Identification and characterization of new BoLA-DRB3 alleles by heteroduplex analysis and direct sequencing. Anim Genet. 1995;26(6):413–17.
Wang K, Sun D, Zhang Y. Sequencing of 15 new BoLA-DRB3 alleles. Int J Immunogenet. 2008;35(4–5):331–32.
Das DN, Sri Hari VG, Hatkar DN, Rengarajan K, Saravanan R, Suryanarayana VV, Murthy LK. Genetic diversity and population genetic analysis of bovine MHC class II DRB 3.2 locus in three Bos indicus cattle breeds of Southern India. Int J Immunogenet. 2012;39(6):508–19.
We thank Drs. Guangai Xue, Veronica de la Barra Diaz, Manuel Ortiz, Jorge Oltra, Misao Kanemaki, Luz Angela Alvarez Franco, and Hideki Saito, the veterinarians in Chile, Argentina, Paraguay, Bolivia, and Peru, and members of the National University of LA PLATA and Universidad Austral de Chile for obtaining samples from many farms. We are grateful to the Support Unit at the Bio-material Analysis, RIKEN BSI Research Resources Center, for help with sequence analysis.
This work was supported by Grants-in-Aid for Scientific Research: A(18255013) (supported the design of the study and the trip cost and sampling cost), A(16H02590) (supported the trip cost and sampling cost), C(16 K08039) (supported the genotyping cost and analyzing cost, and interpretation of data), and C(25450405)) (supported the genotyping cost and analyzing cost, interpretation of data, and in writing the manuscript), from the Japan Society for the Promotion of Science (JSPS).
Availability of data and materials
The datasets used and/or analyzed in the current study are available from the corresponding author upon reasonable request.
All animals were handled by veterinarians from RIKEN, Universidad Austral de Chile and LAVET, in strict accordance with good animal practice following the Universidad Austral de Chile Institutional guidelines. This study was approved by the Committee on the Ethics of Animals for Research at the National University of LA PLATA (Certificate date May 26th, 2014) and by the Committee on the Ethics of Animals for Research at Universidad Austral de Chile (Certificate No. 153–2014). The animal has derived from private owners and the consent was taken in verbal form, because each farmer commonly contacts with veterinarian in same way.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Information and genotyping data used in this study. (XLSX 173 kb)