- Research article
- Open Access
Characterization of bovine MHC DRB3 diversity in global cattle breeds, with a focus on cattle in Myanmar
BMC Genetics volume 21, Article number: 95 (2020)
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.
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.
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 FST analysis showed a low level of genetic differentiation between the two Myanmar native breeds (FST = 0.003), and between these native breeds and the Holstein-Friesians (FST < 0.021). The average FST 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.
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.
Myanmar, located in tropical South East Asia, is divided into seven topographic regions: Northern Hills, Western Hills, Shan Plateau, Central Belt, Lower Myanmar Delta, Rakhine Coastal Region and Tanintharyi Coastal Strip. The country has a tropical monsoon climate with three seasons: hot (mid-February to mid-May), rainy (mid-May to mid-October) and cool season (mid-October to mid-February).
Myanmar’s livestock population consists of about 18.7 million head of cattle which are distributed throughout the country, but most concentrated in the central region . Currently, the Myanmar cattle population can be separated into the two native breeds (Shwe Ni and Pyer Sein) and the Holstein-Friesian crossbreed. These populations are predominantly made up of animals with a low degree of artificial selection but which are adapted to the harsh tropical environment, resistant to tropical diseases and external parasites, and able to thrive on low-quality roughages and grasses . Myanmar native cattle are zebu type (Bos indicus) cattle and are characterized by their geographical location and coat color. The indigenous breeds Shwe Ni and Pyer Sein are mainly used for draught and to a lesser extent dairy production . Only older animals who are no longer fit for work are slaughtered for beef consumption.
Pyer Sein breed was obtained by crossing Indian and native cattle breeds. Indian breeds included Red Shindi, Hariarna and Thari imported during British colonial times and after Myanmar independence to improve the local milk production. Shwe Ni, the true indigenous Myanmar cattle are mostly red or pied and moderately small, compact and with fine bones . In addition, European dairy breeds including Friesians, Jersey, Guernsey and Norwegian Reds were brought into the country in the second half of the twentieth century. Nowadays, the Holstein-Friesian crossbreed is the most popular dairy cow in Myanmar.
The major histocompatibility complex (MHC) is a major component of the adaptive immune system; with MHC genes encoding the cell-surface glycoproteins that bind small peptide fragments derived from host- and pathogen-expressed proteins via proteolysis . The bovine MHC, known as the bovine leukocyte antigen (BoLA), includes three copies of the DRB class II gene with the BoLA-DRB3 variant being the most highly expressed and polymorphic . This locus’ polymorphisms influence both the magnitude and epitope specificity of the antigen-specific T cell response to infectious diseases . Polymorphic sites within the BoLA-DRB3 gene are mainly found at the peptide-binding sites that comprise the α1 and β1 domains of the protein, with most maintained by balancing or overdominance selection [7, 8]. The BoLA-DRB3 gene is associated with resistance/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) [9, 10]. For these reasons, the study of the genetic variability in bovine MHC genes is of interest for evolutionary biologists and animal sciences researchers, as well as, veterinarians, breeders and cattle owners. The analysis of the genetic diversity of the BoLA-DRB3 gene began more than 25 years ago with pioneering studies using serotype analysis [11, 12]. Later, different genotyping techniques were applied (sequencing of cloned genomic DNA, cDNA, cloned PCR products or PCR-RFLP) to several different bovine breeds (e.g., [13,14,15]). Currently, polymerase chain reaction-sequence based typing (PCR-SBT) [16,17,18,19,20,21,22,23,24,25,26,27] and target next generation sequencing (Target-NGS)  are the most commonly employed and powerful tools used in this type of analysis. However, these techniques have been applied in studies which have included only a handful of the over 800 recognized cattle breeds worldwide [16,17,18,19,20,21,22,23,24,25,26,27,28]. Consequently, there are still a number of breeds that remain uncharacterized, and this number increases when local native bovine breeds are considered [17, 18, 26, 29]. As a result of these studies there are 144 described BoLA-DRB3 alleles for bovines, and 303 subtypes listed in the Immuno Polymorphism Database (IPD-MHC) [30, 31]. In addition, previous studies [16,17,18,19,20,21,22,23,24,25,26,27] have shown an even distribution of BoLA-DRB3 polymorphisms between the major bovine types (Bos taurus and B. indicus) and breeds. This is likely as a result of various factors including breed origin, artificial or natural selection and geographical distribution.
The aim of this work was to assess, at both the allele and molecular levels, the inherent genetic diversity and structure of the BoLA-DRB3 gene in Myanmar cattle.
Distribution of BoLA-DRB3 alleles between selected Myanmar cattle breeds
PCR-SBT genotyping allowed us to identify 72 BoLA-DRB3 alleles (69 previously reported variants and three new alleles; Table 1) for the breeds selected in this study. The number of alleles (na) was 57 in Pyer Sein cattle (54 previously reported and three new), 43 in Shwe Ni (41 previously reported and two new), and 33 in the Holstein-Friesian crossbreed (32 previously reported and one new) (Tables 1 and 2). Nucleotide and predicted amino acid sequences of the three new allele variants are shown in Fig. 1.
The three new variants were assigned allele names by IPD-MHC, namely, BoLA-DRB3*002:02, which differs from BoLA-DRB3*002:01 at 147 positions; BoLA-DRB3*079:01, which differs from DRB3*001:01 at seven positions (76, 108, 129, 206, 218, 220 and 254); and BoLA-DRB3*080:01, which differs from BoLA-DRB3*003:01 at five positions (108, 157, 173, 199 and 207). All three new BoLA-DRB3 allele variants shared about 90–94% and 83.7–90.7% nucleotide and amino acid similarity with the BoLA-DRB3 cDNA clone NR1 (correspond to GenBank accession number D45357 and the allele BoLA-DRB3*016:01), respectively .
A Venn diagram was constructed using data obtained in this study and from previous reports [18, 23, 29] which include 102 BoLA-DRB3 alleles. Data were grouped in terms of the breed’s geographical origin as follows: Myanmar native breeds (Pyer Sein and Shwe Ni) and Holstein-Friesian crossbreed; Asian native (Philippine native and Japanese Black); Zebu (Bolivian Nellore, Bolivian Gir, Peruvian Nellore-Brahman and Philippine Brahman); and European (Chilean Hereford, Chilean Black Angus, Chilean Red Angus, Japanese Jersey, Japanese Shorthorn and Japanese Holstein) breeds. Six alleles were not present in any of these breed groups. This analysis revealed that out of the 96 alleles identified in the five cattle groups, only nine were detected in the Myanmar native breeds (Fig. 2a), three of which exhibited gene frequencies that were higher than 0.5% (Fig. 2b). Two other variants were only present in Myanmar native breeds and the Holstein-Friesian crossbreed. Together, these 11 alleles represent about 15% of the 73 alleles detected in the Myanmar native cattle and Holstein-Friesian crossbreed. Twenty-six other alleles were only found in Myanmar cattle populations and Asian native or Zebu breeds, or a combination of these groups. In addition, the BoLA-DRB3 NJ tree including all the previously reported alleles and the three new variants ones showed that the variants detected in Myanmar cattle populations were interspersed among the various clusters (Fig. 3).
Three, six and five alleles appeared with frequencies of > 5% in Pyer Sein, Shwe Ni and Holstein-Friesians, respectively. Three of these high-frequency (> 5%) alleles (BoLA-DRB3*002:01, *017:03 and *022:01) were common in at least two out of three Myanmar populations (Table 1). These common alleles accounted for a low proportion of the cumulative gene frequencies (19.02, 34.78 and 43.13% in Pyer Sein, Shwe Ni and Holstein-Friesians, respectively), revealing an even gene frequency distribution (Fig. S1).
Nucleotide and amino acid diversity in the BoLA-DRB3 alleles found in Myanmar cattle populations
The results of genetic diversity at the DNA and amino acid levels obtained for Myanmar cattle breeds and for breeds previously reported are shown in Table 3. The π values within Pyer Sein, Shwe Ni and Holstein-Friesian crossbreed were 0.090, 0.080 and 0.080, respectively, while the mean number of pairwise differences was 20.96, 17.89 and 20.09, respectively. These nucleotide diversity values all fall within the upper end of the range reported (πrange = 0.068–0.083; NPDrange = 16.31–20.04) for other bovine breeds when using PCR-SBT genotyping methods [18, 22, 23, 26, 29]. The average dN and dS substitutions in Myanmar cattle breeds was calculated across BoLA-DRB3 exon 2 and the antigen-binding site (ABS). As expected, the dN/dS ratio was higher when only the ABS was analyzed. As shown in Table 3, the values obtained in Myanmar cattle were similar to those estimated for other cattle breeds (dN/dS total = 0.054–0.067; dN/dS ABS = 0.247–0.282).
Gene diversity, HWE, and neutrality testing of BoLA-DRB3 variants found in Myanmar cattle populations
As mentioned above, na ranged from 33 in the Holstein-Friesian crossbreed to 58 in the Pyer Sein native breed, while he and ho were both higher than 0.86 in all three populations (Table 2). These indexes are evidence of high diversity values for Myanmar cattle populations, which is similar to the results reported for other bovine breeds which have been evaluated by PCR-SBT, and characteristic of MHC class II DR genes [18, 22, 23, 26, 29]. When we evaluated the populations using the HWE test, two of the three Myanmar populations were in equilibrium, while native breed Shwe Ni significantly deviated from its theoretical values (Table 2), probably because of the significant proportions of homozygous animals found in this study (FIS = 0.1197, p < 0.0001). As demonstrated in Table 2, other breeds have also been seen to be in disequilibrium, as a result of an excess or deficit in the proportion of homozygous animals within the population.
It is widely accepted that the genetic diversity of MHC class II genes can be maintained by balancing selection. Thus, we performed a Slatkin’s exact neutrality test (Table 2) to evaluate this phenomenon in our populations. The BoLA-DRB3 gene frequency profile in both the Pyer Sein and Shwe Ni cattle showed an even distribution (p = 0.006 and 0.010, respectively), consistent with the theoretical proportion expected under balancing selection pressures as opposed to positive or neutral selection (p > 0.025). A similarly even BoLA-DRB3 gene frequency was observed in other cattle breeds, including Japanese Black, Yacumeño Creole and Bolivian Gir. Conversely, BoLA-DRB3 gene frequency distributions in the Holstein-Friesian crossbreed were more compatible with neutral selection, which is similar to the results obtained for the majority of the cattle breeds analyzed to date (Table 2).
BoLA-DRB3 genetic structure and levels of population differentiation in Myanmar cattle
The average FST analysis showed a low level of genetic differentiation between Myanmar native breeds (FST = 0.003), similar to those estimated in Holstein populations (0–0.0067) . FST values between the native breeds and the Holstein-Friesian crossbreed varied from 0.019 to 0.021. These values were within the range estimated for differences within Taurine or Zebu breeds and lower than those obtained when comparing breeds from different groups (Fig. 4 and Table S2).
The average FST values across all Myanmar Holstein-Friesian and native breed populations were 0.0136 and 0.0121, respectively (p < 0.001). Significant differences were observed in nine out of the fifteen native and one out of six Holstein-Friesian crossbreed populations (p < 0.05). In addition, FST values for comparisons between Myanmar native breeds ranged between 0.003 and 0.024 and between 0 and 0.031 for Myanmar Holstein-Friesian crossbreed populations (Table S3a and b). As mentioned above, similar genetic distance values were observed among Holstein populations from different countries .
Genetic differentiation of BoLA-DRB3 alleles in Myanmar breeds. Comparison with zebu and taurine breeds
First, BoLA-DRB3 allele frequencies from Myanmar and previously reported breeds included in our dataset were used to generate Nei’s DA and DS genetic distance matrices. Then, dendrograms were constructed from these distance matrices using UPGMA and NJ algorithms. All trees revealed congruent topologies, which were consistent with the historical and geographical origin of the breeds. As expected, this tree revealed two main clusters which included most of the Taurine and Zebuine breeds (Fig. 5a), with Japanese Jersey and Chilean Hereford located outside of these clusters. The Holstein-Friesian crossbreed fell into the Taurine cluster, while the Myanmar native breeds were located in a sub-cluster within the Zebuine cluster close to Philippine populations.
Second, we used BoLA-DRB3 allele frequencies to perform three PCA analyses among breeds. In these PCA, the first two components accounted for 41.11% of the data variability. The first PC accounted for 27.12% of the total variance and, as shown in a previous study , first component clearly exhibited a differentiation pattern between Zebu (negative values) and Taurine (positive values) breeds, while native breeds from Myanmar and the Philippines were located near the origin of the plot (Fig. 5b). The first PC was primarily determined by differences in the frequency of the same alleles reported by Takeshima et al. . The second PC explained 13.99% of the total variation and showed a gradient among Taurine breeds, with Japanese Black and Japanese Jersey located at opposite ends. Furthermore, this component discriminated between Myanmar and Philippine native breeds. The second PC was identical to PC1 reported in the study mentioned above . Finally, the third PC accounted for 13.64% of the variance and allowed the differentiation of Chilean Hereford cattle from other Taurine breeds. As shown in Fig. 5b and c, the Myanmar Holstein-Friesian crossbreed was located within the Taurine cloud but in an intermediate position between the Japanese Holstein and Myanmar native breeds, supporting the presence of the same level of gene introgression in Myanmar Holstein populations, which is also supported by the presence of unique BoLA-DRB3 alleles within these populations. These PCA results agree with the overall clustering observed after NJ or UPGMA tree construction.
In addition, we analyzed protein pockets (pocket 1, pocket 4, pocket 6, pocket 7, and pocket 9) involved in the antigen-binding function of the MHC complex using PCA. As shown in Fig. S2 a-e, only the distribution pattern of pocket 4 was similar to the PCAs created using the allelic frequency data, while PCA of the remaining pockets did not exhibit a spatial distribution related to the geographical or historical origin of the breeds. The position of the Myanmar native breeds in pocket 4 was the result of positive PC1 and PC2 values for the presence of amino acid motifs GFDQKEV, SYDRENY, SFDREYY, SFDDEAY, KFDRAAY, and GYDREYY (amino acid positions 13, 26, 28, 70, 71, 74 and 78). Also, pocket 1 showed quite similar distribution pattern, but for PC2.
Finally, PCA was performed at the Myanmar population level, to evaluate the impact of the ten sampling sites (six for native breeds and four for Holstein-Friesians) on our results. This analysis showed that Myanmar native populations grouped in a narrow cluster that diverged clearly from the Myanmar Holstein-Friesian crossbred populations (Fig. 6), in agreement with the FST analysis described above. Furthermore, PCA showed that some Myanmar native populations (Bago, Mandalay and Yangon) seemed to be closer to Zebu breeds (Nellore, Gir and Brahman), while others (Kayin, Magway and Sagaing) were more closely related to the Philippine native breeds. However, PCA results at the Myanmar population level did not show a clear correlation between the genetic relationship of BoLA-DRB3 alleles and geographical distribution. By contrast, Myanmar Holstein-Friesian populations showed a more dispersed distribution when compared to the compact cloud reported for Holstein populations from other countries , which may be the result of the differences in the degree of admixture between these populations.
In this study, we carried out the first genetic characterization of the BoLA-DRB3 gene in two Myanmar native breeds and the Myanmar Holstein-Friesian crossbreed population using PCR-SBT. This analysis allowed us to detect 71 alleles, including three new variants. Analyses of D-loop mitochondrial haplotypes revealed a large percentage of novel haplotypes in Myanmar native breeds, suggesting that these indigenous populations may possess novel polymorphisms throughout the genome of these native breeds [32, 33]. The BoLA-DRB3 allele NJ tree, generated using the nucleotide sequence of the β1 domain of all the reported alleles and the three new ones, evidenced that the variants detected in Myanmar cattle populations were not grouped in specific clusters of the dendrogram, but were interspersed along the tree. Similar results were reported in South American native cattle breeds  and agree with the ancient origin of the MHC alleles proposed by the trans-specific theory of MHC alleles . According to this theory, these alleles originated before cattle domestication and speciation, meaning that many variants will be common across species and geographies. This is supported by the fact that BoLA-DRB3*1303 (gi|1,261,491,414|gb|KY094633.1|) which was detected in this work was previously described by a study focused on Indonesian Bos javanicus.
In order to understand the distribution of the BoLA-DRB3 alleles across breeds, a Venn diagram was constructed. This analysis demonstrated that 11 BoLA-DRB3 alleles were only detected in Myanmar cattle populations, at least for the breeds included in this analysis. Three of these alleles were entirely unique variants, described in this study (BoLA-DRB3*002:02, BoLA-DRB3*079:01 and BoLA-DRB3*080:01). It is worth noting that a review of the IPD–MHC database showed that some variants had been previously detected in other Asian breeds (BoLA-DRB3*031:03, *058:01 and *064:02 in Chinese yellow cattle, BoLA-DRB3*013:03 in Korean Hanwoo cattle and BoLA-DRB3*049:01 in the Indian Red Sindhi); African Zebu breeds (BoLA-DRB3*038:01 in Ethiopian Arsi and BoLA-DRB3*013:03 in Sudanese Baggara) and other Asian Bos species, including BoLA-DRB3*013:03 in Indonesian Bos javanicus. In previous studies privative alleles were also identified in native breeds from Asia and South America [21, 26, 28]. In addition, some alleles detected in Creole cattle have also been detected in African breeds (Boran, N’Dama, Ethiopian Arsi, and Gudali) and Bos indicus (Brahman). Together, these results show that BoLA-DRB3 variants could be classified using the following categories: worldwide geographical distribution, present only in one major bovine type (Taurus or Indicus), located in a geographical region, and detected in only one or a few related breeds. The current geographical allele distribution could be the result of several factors, including breed origin, founder effect, natural or artificial selection and recent or historical gene introgression.
The large number of BoLA-DRB3 variants indicate a high level of genetic diversity in Myanmar cattle populations. High genetic variability was not unexpected for the BoLA-DRB3 locus, with the values recorded in this study falling in line with observations made in Philippine native and Hanwoo cattle breeds [21, 28]. A high level of diversity was also observed in DNA sequences when nucleotide diversity and mean number of pairwise differences were estimated. In addition, a higher number of dn than ds changes (especially in antigen peptide binding sites) were found in this population which is also similar to the findings reported for other cattle breeds. Different mechanisms for the maintenance of the high diversity of MHC loci have been proposed, including balancing selection (i.e., overdominance and negative frequency dependence) [35,36,37].
The neutrality test showed that both Myanmar indigenous breeds had an even allele frequency distribution, while the most common allele (> 0.5) accounted for only 19.02% of alleles in Pyer Sein and 34.78% of alleles in Shwe Ni. These results are compatible with the hypothesis of balancing selection rather than positive selection against one or more variants or neutral selection strategies. This finding has also been reported in Yacumeño creole cattle as well as in other cattle breeds (Bolivian Gir, Japanese Jersey and Japanese Black) .
Takeshima et al.  reported a significant deviation from HWE in class II BoLA-DQA1 genes in Japanese Holstein cows with mastitis caused by Escherichia or Streptococcus bacteria. Similar results were observed in human, model and non-model species . This evidence was interpreted using the overdominance selection theory, which proposes that heterozygous individuals recognize a broader spectrum of foreign antigens than homozygous ones, thereby increasing their fitness and selective advantage against pathogens [39, 40]. In this study, only the Shwe Ni native breed showed a significant deviation from HWE theorical proportions resulting from an excessive proportion of homozygotic animals within the population. However, none of the Myanmar cattle populations evidenced a significant excess of heterozygotes. These results agree with data observed in other cattle populations exhibiting HWE or excessive numbers of homozygotes, and could be the consequence of low values for overdominance selection coefficients in respect of the stochastic forces (e.g., genetic drift, inbreeding), population structure (Whalund effect) and the low-resolution power of HWE methods .
The genetic structure and breed relationships of BoLA-DRB3 polymorphisms were studied using an FST index, PCA and dendrograms. Together, these analyses showed that Myanmar native breeds, which are mainly characterized by their geographical location and coat color, exhibited a low level of genetic divergence and tended to be closely related in both PCA and dendrograms. These results were reinforced when PCA was performed at the population level, and was able to demonstrate that there was a low level of genetic differentiation among Myanmar cattle populations (FST = 0.0121). These breeds grouped in a narrow cluster with a low level of genetic distance between them. Furthermore, no correlation between genetic and geographical distance for the BoLA-DRB3 alleles was observed among these populations. In previous studies, assignation to the Chinese yellow cattle populations using microsatellites and SNP markers did not allow researchers to distinguish unique genetic differences for each of the Chinese breeds, and individuals could not be correctly assigned to their breed based on their genotype, despite two specific yellow cattle genetic profiles being described [41, 42]. For these reasons native cattle breeds from Eastern Asia are usually described based on their geographical region, without breeding associations or herd books that create breeding barriers, mainly characteristic of the European breeds. These types of analysis that allow to identify genetic entities are important in the development and design of conservation programs created to protect these native cattle breeds.
In addition, PCA and dendrograms showed that Myanmar native breeds were located close to Philippine populations and fell into a sub-cluster within the Zebuine cluster. Furthermore, PCA showed that some Myanmar native populations (Bago, Mandalay, and Yangon) seemed to be closer to some Zebu breeds (Gir, and Brahman), while others (Kayin, Magway, and Sagaing) were more closely related to the Philippine native breeds. This is in agreement with Myanmar’s location in South East Asia, where it shares borders with Bangladesh and India to the Northwest, China to the Northeast, Laos to the East and Thailand to the Southeast, and the fact that Myanmar native breeds are phenotypically humped cattle. Studies based on blood proteins and SRY gene polymorphisms support our conclusion that these indigenous cattle are zebu type animals [43, 44]. In addition, Chen et al.  and Lwin et al.  showed that Myanmar native cattle had zebu mitochondrial haplotypes, confirming our hypothesis.
When breed relationships were studied using the gene frequencies described for the amino acid motifs found in the ABS pockets implicated in antigen-binding function, only the PCA from pocket 4 demonstrated a similar distribution pattern to the PCA pattern developed using allelic frequency data. By contrast, the remaining PCAs did not exhibit a spatial distribution compatible with the geographical or historical origin of the analyzed breeds. The position of the Myanmar native breeds in the pocket 4 PCA was as a result of their positive PC1 and PC2 values for amino acid motifs GFDQKEV, SYDRENY, SFDREYY, SFDDEAY, KFDRAAY, and GYDREYY. Previous studies showed that pocket 4 has been found to be important for the binding of peptides may be due to this pocket being located in the center of the PBC . In addition, it has been reported in cattle that immune responses and disease resistance is significantly related to differences in the pocket 4 motif [45, 46]. Takeshima et al.  reported that PCAs based on allelic frequencies in the amino acid motif of the ABS pockets showed a clear differentiation between Taurine and Zebu breeds, because of an enrichment for particular amino acid motifs in specific pockets. Further studies involving structural modeling and molecular simulation of the BoLA-DRβ protein are needed to elucidate whether these differences play a role in its function.
The Myanmar Holstein-Friesian crossbreed was located within the Taurine cloud but in an intermediate position between Japanese Holstein and Myanmar native breeds, probably as a consequence of gene introgression in Myanmar Holstein-Friesian cattle, which is also supported by the presence of BoLA-DRB3 alleles unique to the native and Holstein-Friesian breeds from Myanmar. The average level of genetic differentiation observed among the Myanmar Holstein-Friesian crossbreed populations (FST = 0.0136) was higher than those previously reported in Holstein populations from other countries (FST = 0.009) . Its distribution was more dispersed in the PCA when compared with the more compact cloud reported for other Holstein populations . These differences are likely the result of the different degree of admixture between this crossbreed and the Myanmar native cattle.
Holstein cattle are mainly raised in temperate and cold regions because they are sensitive to heat stress and susceptible to tropical disease. However, this breed has been introduced to tropical and subtropical countries several times with varying degrees of success, and it is possible to find Holstein crossbreed populations in these environments. In 1978, pregnant Holstein cows from New Zealand and Australia were imported to Myanmar to improve the local dairy cows . Later, semen was introduced from North America and Europe [47, 48]. The present BoLA-DRB3 gene results would indicate that foreign Holstein cattle could be crossed with local cattle to create an adapted Holstein crossbreed population. To support this, four BoLA-DRB3 alleles (BoLA-DRB3*022:01, BoLA-DRB3*026:01, BoLA-DRB3*031:01 and BoLA-DRB3*079:01), including the novel BoLA-DRB3*079:01 variant, were present in the Myanmar Holstein-Friesian population. Notably, these alleles were common in the Myanmar native cattle population, but they had a very low frequency or were absent in Holstein populations from temperate or cold regions. This is consistent with the results reported by Takeshima et al. , which demonstrated that Bolivian Holsteins, raised in the humid subtropical plains of Santa Cruz (Bolivia), diverged from cattle bred in temperate and cold geographical regions. In addition, the Bolivian Holstein population was found between Yacumeño Creole cattle and the compact cloud of Holstein populations when evaluated as part of a larger data set and shared common alleles with this Creole breed. In this sense, the introgression of genes from locally adapted native cattle in Holstein-like populations following natural selection may have contributed to increased fitness in tropical and subtropical regions.
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, using the PCR-SBT assay. These results revealed the presence of three new alleles and demonstrated a high degree of genetic diversity in this gene, which could contribute to the adaptation of these breeds to the harsh subtropical environmental of Myanmar. Furthermore, dendrograms and PCA showed that Myanmar native breeds were closely related to one another and to the Philippine native breeds at this locus. Therefore, this study increases our knowledge of the worldwide variability of cattle BoLA-DRB3 genes, an important locus for immune response and protection against pathogens.
Sample populations and genomic DNA extraction
Blood samples were obtained from 294 cattle belonging to the Myanmar Pyer Sein (N = 163), Shwe Ni (N = 69) native breeds and the Holstein-Friesian crossbreed (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) (Table 1; Fig. S3). Blood samples were randomly collected from adult animals according to the availability and permission of local farmers in the respective regions. A database that included 2428 BoLA-DRB3 genotypes from European, Zebuine and Native breeds was also included (Table S1) [18, 23, 26, 29]. These genotypes data were used for recalculating the genetic diversity indexes at alleles and sequence level in breeds details in Table S1 for further comparison with Myanmar results. The genomic DNA was extracted from whole blood using a PureLink Genomic DNA Mini Kit (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s instructions.
BoLA-DRB3 alleles were genotyped using PCR-SBT. Briefly, BoLA-DRB3 exon 2 was amplified using PCR according to the method described by Takeshima et al. . 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). The raw sequence data were analyzed using Assign 400ATF ver. 126.96.36.199 software (Conexio Genomics, Fremantle, Australia).
Measures of genetic variability
Allele frequencies and the observed number of alleles (na) were obtained by direct counting. The distribution of alleles across breeds was represented by a Venn plot created using the R package ‘VennDiagram’ . The observed (ho) and unbiased expected (he) heterozygosity of the BoLA-DRB3 locus were estimated according to the values and assumptions described in Nei  using the Arlequin 3.5 software for population genetic analyses . Potential deviations from Hardy-Weinberg equilibrium (HWE) were estimated using FIS statistics  for each breed calculated using the exact test included in Genepop 4.7 software . The Ewens–Watterson–Slatkin exact test of neutrality was estimated using the method described by Slatkin  and implemented in the Arlequin 3.5 program.
Genetic structure and population differentiation
Genetic structure and genetic differentiation among the breeds were assessed using Wright’s FST statistics, calculated using the variance-based method of Weir and Cockerham . This parameter was estimated using Arlequin 3.5 and Genepop 4.7 software. The FST values were represented graphically using the pairFstMatrix.r function implemented in the Arlequin 3.5 software.
Analysis of the genetic relationship between breeds
To condense the genetic variation at the BoLA-DRB3 locus, allele frequencies were used to perform a Principal Component Analysis (PCA) according to the Cavalli-Sforza  method, implemented in Past software . Nei’s standard genetic distances Ds  and DA  were calculated from allele frequencies to perform a cluster analysis using the unweighted pair-group method with arithmetic mean (UPGMA)  and the neighbor-joining algorithm (NJ) . Confidence intervals for the groupings were estimated by bootstrap re-sampling of the data using 1000 replicates. Genetic distances and trees were computed using the Populations 1.2.28 software . The trees were then visualized using TreeView .
Genetic diversity at the sequence level
Nucleotide diversity (π) and pairwise comparisons of nucleotide substitutions between alleles (defined as the average number of differences between pairs of DNA sequences) were calculated using Arlequin 3.5. The mean number of nonsynonymous (dN) and synonymous (dS) nucleotide substitutions per site calculated as an average over all sequence pairs was estimated within each group using the Nei-Gojobori model  and Jukes–Cantor’s formula. These parameters were estimated using Arlequin 3.5 and MEGA X . The BoLA-DRB3 allele tree was constructed using a distance matrix based on the NJ method. To test the significance of each of the branches, 1000 bootstrap replicate calculations were performed. The allele tree was constructed using MEGA X software.
Availability of data and materials
The data sets supporting the results of this article are included in this manuscript and its additional information files. In addition, the new BoLA-DRB3 variants generated during the current study are available in the IPD-MHC database (accession numbers BoLA09869, BoLA09870, BoLA09871; https://www.ebi.ac.uk/ipd/mhc/) and in the European Nucleotide Archive (ENA) repository (accession numbers LC466585, LC466586, LC466588; https://www.ebi.ac.uk/ena/).
Bovine leukocyte antigen
Major histocompatibility complex
Antigen recognition site
Next generation sequencing
Polymerase chain reaction
- PCR- SBT:
Polymerase chain reaction-sequence-based typing
Highland Bolivian Creole cattle
Mean number pairwise difference
- dn :
- ds :
- na :
Number of alleles
- ho :
- he :
Principal component analysis
- FST :
FAOSTAT database; http://www.fao.org/faostat/en/#data/QA. Accessed 20 Mar 2020.
Polat M, Moe HH, Shimogiri T, Moe KK, Takeshima SN, Aida Y. The molecular epidemiological study of bovine leukemia virus infection in Myanmar cattle. Arch Virol. 2017;162(2):425–37.
Daing T. Animal genetic resources in Myanmar. Report Soc Res Native Livest. 2004;21:9–18.
Klein J. Natural History of the Major Histocompatability Complex. 99th ed. New York: Wiley; 1986.
Aida Y. Characterization and expression of bovine MHC class II genes. Bull Soc Fr Jpn Sci Vet. 1995;6:17–24.
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.
Hughes AL, Nei M. Nucleotide substitution at major histocompatibility com plex class II loci: evidence for overdominant selection. Proc Natl Acad Sci U S A. 1989;86:958–62.
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.
Takeshima SN, Aida Y. Structure, function and disease susceptibility of the bovine major histocompatibility complex. Anim Sci J. 2006;77:138–50.
Nayeri S, Stothard P. Tissues, metabolic pathways and genes of key importance in lactating dairy cattle. Springer Sci Rev. 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.
Sigurdardottir S, Borsch C, Gustafsson K, Andersson L. Cloning and sequence analysis of 14 DRB alleles of the bovine major histocompatibility complex by using the polymerase chain reaction. Anim Genet. 1991;22:199–209.
Mikko S, Andersson L. Extensive MHC class II DRB3 diversity in African and European cattle. Immunogenetics. 1995;42:408–13.
Giovambattista G, Golijow CD, Dulout FN, Lojo MM. Gene frequencies of DRB3.2 locus of argentine creole cattle. Anim Genet. 1996;27(1):55–6.
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.
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.
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.
Bohórquez MD, Ordoñez D, Suárez CF, Vicente B, Vieira C, López-Abán J, et al. Major histocompatibility complex class II (DRB3) genetic diversity in Spanish Morucha and Colombian Normande cattle compared to taurine and zebu populations. Front Genet. 2020;10:1293.
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.
Takeshima SN, Corbi-Botto C, Giovambattista G, Aida Y. Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations. BMC Genet. 2018;19(1):33.
IPD–MHC database, https://www.ebi.ac.uk/ipd/mhc/group/BoLA. Accessed 20 Mar 2020.
Maccari G, Robinson J, Ballingall K, Guethlein LA, Grimholt U, Kaufman J, Ho CS, De Groot NG, Flicek P, Bontrop RE, Hammond JA, Marsh SGE. IPD-MHC 2.0: an improved inter-species database for the study of the major histocompatibility complex. Nucleic Acids Res. 2017;45:D860–4.
Chen S, Lin BZ, Baig M, Mitra B, Lopes RJ, Santos AM, Magee DA, Azevedo M, Tarroso P, Sasazaki S, Ostrowski S, Mahgoub O, Chaudhuri TK, Zhang YP, Costa V, Royo LJ, Goyache F, Luikart G, Boivin N, Fuller DQ, Mannen H, Bradley DG, Beja-Pereira A. Zebu cattle are an exclusive legacy of the South Asia Neolithic. Mol Biol Evol. 2010;27:1–6.
Lwin M, Mon SLY, Nagano Y, Kawabe K, Mannen H, Kamoto S, Shimogiri T. Genetic diversity of Myanmar cattle breeds using complete mitochondrial D-loop sequence. J Anim Genet. 2018;46:57–67.
Klein J, Satta Y, Takahata N, O'hUigin C. Trans-species MHC polymorphisms and the origin of species in primates. J Med Primatol. 1993;22:57–64.
Clarke B, Kirby DR. Maintenance of histocompatibility polymorphisms. Nature. 1966;211(5052):999–1000.
Doherty P, Zinkernagel R. Enhanced immunological surveillance in mice heterozygous at the H-2 gene complex. Nature. 1975;256:50–2.
Piertney SB, Oliver MK. The evolutionary ecology of the major histocompatibility complex. Heredity. 2006;96(1):7–21.
Takeshima SN, Matsumoto Y, Chen J, Yoshida T, Mukoyama H, Aida Y. Evidence for cattle major histocompatibility complex (BoLA) class II DQA1 gene hetero-zygote advantage against clinical mastitis caused by streptococci and Escherichia species. Tissue Antigens. 2008;72:525–31.
Hedrick PW, Thomson G, et al. Genetics. 1983;104(3):449–56 52.
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.
Rogberg-Muñoz A, Wei S, Ripoli MV, Guo BL, Carino MH, Castillo N, Villegas Castagnaso EE, Lirón JP, Morales Durand HF, Melucci L, Villarreal E, Peral-García P, Wei YM, Giovambattista G. Foreign meat identification by DNA breed assignment for the Chinese market. Meat Sci. 2014;98(4):822–7.
Rogberg-Muñoz A, Wei S, Ripoli MV, Guo BL, Carino MH, Lirón JP, Prando AJ, Vaca RJ, Peral-García P, Wei YM, Giovambattista G. Effectiveness of a 95 SNPpanel for the screening of breed label fraud in the Chinese meat market. Meat Sci. 2016;111:47–52.
Nomura K, Takahashi Y, Amano T, Yamagata T, Mannen H, Nishibori M, Yamamoto Y, Kurosawa Y, Tanaka K, Okabayashi H, Tsunoda K, Nozawa K, Nyunt MM, Daing T, Hla T, Win N. Composition of genes controlling blood protein types of native cattle and mithan in Myanmar and their phylogenetic study. Report Soc Res Native Livest. 2004;21:77–89.
Tanaka K, Kurita J, Mannen H, Kurosawa Y, Nozawa K, Nishibori M, Yamamoto Y, Okabayashi H, Tsunoda K, Yamagata T, Suzuki Y, Kinosita K, Maeda Y, Nyunt M, Daing T, Hla T, Win N, Tur T, Aung P, Cho A. An assay for paternal gene flow between the taurus and indicus cattle in Myanmar using polymorphisms in SRY gene. Report Soc Res Native Livest. 2004;21:91–9.
Sharif S, Mallard BA, Sargeant JM. Presence of glutamine at position 74 of pocket 4 in the BoLA-DR antigen binding groove is associated with occurrence of clinical mastitis caused by Staphylococcus species. Vet Immunol Immunopathol. 2000;6:231–8.
Baxter R, Craigmile SC, Haley C, Douglas AJ, Williams JL, Glass EJ. BoLA-DR peptide binding pockets are fundamental for foot-and-mouth disease virus vaccine design in cattle. Vaccine. 2009;28(1):28–37.
National Consultative Committee, Myanmar. 2002. National Report on Animal Genetic Resources, the Union of Myanmar; [Cited February 2018]. Available from http://www.fao.org/tempref/docrep/fao/010/a1250e/annexes/CountryReports/Myanmar.pdf.
Than H, Aung M, Daw SK, Ye HW. Improvement of cattle production in Myanmar through the use of progesterone RIA to increase efficiency and quality of artificial insemination services. In: Radioimmunoassay and related techniques to improve artificial insemination programmes for cattle reared under tropical and sub-tropical conditions, Proceedings of a final Research Co-ordination Meeting; 1999. p. 45–50.
Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12(1):35.
Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89: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 32.
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, vol. 35. 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 36.
Nei M. Genetic distance between populations. Am Nat. 1972;106(949):283 37.
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 38.
Sneath PHA, Sokal RR. Numerical taxonomy: The principles and practice of numerical classification. Hardcover edn. United States: Freeman & Company Limited, W H; 1973. 39.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25 40.
Langella O. Populations: a population genetic software. In: vol. CNRS UPR9034., 1.2.28 edn; 1999. http://www.mybiosoftware.com/populations-12-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 42.
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.
We thank LBVD staff, UVS staff, local veterinarians and farm owners for their kind assistance in blood sampling at the cattle farms in Myanmar, and the Support Unit at the Biomaterial Analysis, RIKEN BSI Research Resource Center for help with the sequence analysis. We also thank all members of the Viral Infectious Diseases Field, RIKEN, for technical assistance, help, and suggestions.
The study was supported by grants from the Project of the NARO Bio-oriented Technology Research Advancement Institution [the Special Scheme Project on Regional Developing Strategy (Grant No. 16817983) and the Special Scheme Project on Vitalizing Management Entities of Agriculture, Forestry and Fisheries (Grant No. 16930548)], and was supported a JSPS International Research Fellow [FY2018 Invitational Fellowship for Research in Japan (Long-term)]. These funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All animals were handled by RIKEN, Japan in strict accordance with good animal practice following the guidelines of RIKEN. The study was approved by the RIKEN Animal Experiments Committee (approval number H29–2-104). The written informed consent to use the animals in the present study were obtained from the owners of the cattle.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Cumulative gene frequency plot of BoLA-DRB3 alleles in Pyer Sein (red), Shwe Ni (violet) and Holstein-Friesian crossbreed (orange) populations.
Principal component analysis of BoLa-DRB3 gene pocket amino acid motifs frequencies in 18 populations: a. Pocket 1, b. Pocket 4, c. Pocket 6, d. Pocket 7, and e. Pocket 9. BW = Pyer Sein, GR = Shwe Ni, HoMy = Myanmar Holstein-Friesian crossbreed, WaJa = Japanese Black, HoJa = Holstein, ShJa = Japanese Shorthorn, JeJa = Jersey, HeCh = Chilean Hereford, BACh = Chilean Black Angus, RACh = Chilean Red Angus, ONCh = Chilean Overo Negro, OCCh = Chilean Overo Colorado, NaPh = Philippine Native, NaxBrPh = Native x Brahman Philippine crossbreed, BrPh = Philippine Brahman, NeBo = Bolivian Nellore, GirBo = Bolivian Gir, and BrxNe = Peruvian Brahman × Nellore crossbreed.
Geographical distribution of sampling sites of Myanmar cattle that are indicated by a dot.
Detailed information about the populations analyzed.
Genetic distance between pairs of populations estimated by FST in (a) six Myanmar native (KN = Kayin, BN = Bago, SN = Sagaing, MdN = Mandalay, MgN = Magway, and YN = Yangon) and (b) four Myanmar Holstein-Friesian crossbreed (KF = Kayin, BF = Bago, SF = Sagaing, and YF = Yangon) populations.
Genetic distance between pairs of populations estimated by FST in (a) six Myanmar native (KN = Kayin, BN = Bago, SN = Sagaing, MdN = Mandalay, MgN = Magway, and YN = Yangon), and (b) four Myanmar Holstein-Friesian crossbreed (KF = Kayin, BF = Bago, SF = Sagaing, and YF = Yangon) populations.
About this article
Cite this article
Giovambattista, G., Moe, K.K., Polat, M. et al. Characterization of bovine MHC DRB3 diversity in global cattle breeds, with a focus on cattle in Myanmar. BMC Genet 21, 95 (2020). https://doi.org/10.1186/s12863-020-00905-8
- Myanmar cattle
- Breed diversity
- Breed/population identification