Major Histocompatibility Complex Class II (DRB3) Genetic Diversity in Spanish Morucha and Colombian Normande Cattle Compared to Taurine and Zebu Populations

First characterization of major histocompatibility complex class II DRB3 diversity in cattle breeds raised in Egypt

Genes encoding bovine leukocyte antigen (BoLA) enable the immune system to identify pathogens. Therefore, these genes have been used as genetic markers for infectious and autoimmune diseases as well as for immunological traits in cattle. Although BoLA polymorphisms have been reported in various cattle breeds worldwide, they have not been studied in cattle populations in Egypt. In this study, we characterized BoLA-DRB3 in two local Egyptian populations and one foreign population using polymerase chain reaction-sequence-based typing (PCR-SBT) method. Fifty-four previously reported BoLA-DRB3 alleles and eight new alleles (BoLA-DRB3*005:08, *015:07, *016:03, *017:04, *020:02:02, *021:03, *164:01, and *165:01) were identified. Alignment analysis of the eight new alleles revealed 90.7-98.9 %, and 83.1-97.8 % nucleotide and amino acid identities, respectively, with the BoLA-DRB3 cDNA clone NR-1. Interestingly, BoLA-DRB3 in Egyptian cattle showed a high degree of allelic diversity in native (na = 28, hE > 0.95), mixed (na = 61, hE > 0.96), and Holstein (na = 18, hE > 0.88) populations. BoLA-DRB3*002:01 (14.3 %), BoLA-DRB3*001:01 (8.5 %), and BoLA-DRB3*015:01 (20.2 %) were the most frequent alleles in native, mixed, and Holstein populations, respectively, indicating that the genetic profiles differed in each population. Based on the allele frequencies of BoLA-DRB3, genetic variation among Egyptian, Asian, African, and American breeds was examined using Nei's distances and principal component analysis. The results suggested that native and mixed cattle populations were most closely associated with African breeds in terms of their gene pool, whereas Holstein cattle were more distinct from the other breeds and were closely related to Holstein cattle populations from other countries.

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

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

Identifying major histocompatibility complex class II-DR molecules in bovine and swine peripheral blood monocyte-derived macrophages using mAb-L243

Major histocompatibility complex class II (MHC-II) molecules are involved in immune responses against pathogens and vaccine candidates' immunogenicity. Immunopeptidomics for identifying cancer and infection-related antigens and epitopes have benefited from advances in immunopurification methods and mass spectrometry analysis. The mouse anti-MHC-II-DR monoclonal antibody L243 (mAb-L243) has been effective in recognising MHC-II-DR in both human and non-human primates. It has also been shown to cross-react with other animal species, although it has not been tested in livestock. This study used mAb-L243 to identify Staphylococcus aureus and Salmonella enterica serovar Typhimurium peptides binding to cattle and swine macrophage MHC-II-DR molecules using flow cytometry, mass spectrometry and two immunopurification techniques. Antibody cross-reactivity led to identifying expressed MHC-II-DR molecules, together with 10 Staphylococcus aureus peptides in cattle and 13 S. enterica serovar Typhimurium peptides in swine. Such data demonstrates that MHC-II-DR expression and immunocapture approaches using L243 mAb represents a viable strategy for flow cytometry and immunopeptidomics analysis of bovine and swine antigen-presenting cells.

Kompetitive Allele Specific PCR Genotyping of 89 SNPs in Romanian Spotted and Romanian Brown Cattle Breeds and Their Association with Clinical Mastitis

Mastitis is the most common production disease in the dairy sector worldwide, its incidence being associated with both cows' exposure to bacteria and the cows' genetic make-up for resistance to pathogens. The objective of our study was to analyse 89 missense SNPs belonging to six genes (CXCR2, CXCL8, TLR4, BRCA1, LTF, BOLA-DRB3), which were found to be associated with genetic resistance or susceptibility to mastitis. A total of 298 cattle (250 Romanian Spotted and 48 Romanian Brown) were genotyped by Kompetitive Allele Specific PCR (KASP) and a chi-squared test was used for genetic association studies with clinical mastitis. A total of 35 SNPs (39.3%) among the selected 89 SNPs were successfully genotyped, of which 31 markers were monomorphic. The polymorphic markers were found in two genes: TLR4 (rs460053411) and BOLA-DRB3 (rs42309897, rs208816121, rs110124025). The polymorphic SNPs with MAF > 5% and call rates > 95% were used for the association study. The results showed that rs110124025 in the BOLA-DRB3 gene was significantly associated with mastitis prevalence (p ≤ 0.05) in both investigated breeds. Current results show that the SNP rs110124025 in the BOLA-DRB3 gene can be used as a candidate genetic marker in selection for mastitis resistance in Romanian dairy cattle.

Immunoinformatics analysis of candidate proteins for controlling bovine paratuberculosis

BACKGROUND

Paratuberculosis is debilitating chronic enteritis usually characterized by diarrhea, decreased milk production, and progressive cachexia. Mycobacterium avium subspecies paratuberculosis (MAP) causes significant economic losses by affecting dairy herds globally. Development of protective vaccines is considered as one of the most effective controlling measures for MAP infections. In the current study, hydrophilic parts of MAP2191 and FAP-P proteins as two vaccine candidates were analyzed using immunoinformatics approaches.

METHODS

After selecting the most hydrophilic parts of MAP2191 and FAP-P, helper and cytotoxic T-cell epitopes of ht-MAP2191 and ht-FAP-P were identified. The immunogenic, toxicity and physicochemical properties were assessed. Secondary structures of these proteins were predicted, and their tertiary structures were modeled, refined, and validated. Linear and conformational epitopes of corresponding B-cells were recognized. Then ht-MAP2191 and ht-FAP-P epitopes were employed for molecular docking simulations.

RESULTS

The results indicated that ht-MAP2191 and ht-FAP-P were immunogenic, non-allergenic, and non-toxic and possess potent T-cell and B-cell epitopes. Eventually, these protein constructs were docked favorably against TLR4.

CONCLUSION

According to the findings, ht-MAP2191 and ht-FAP-P could be effective protein-based vaccine candidates for paratuberculosis. It should be noted that to examine their efficacy, further in vitro and in vivo experiments are underway.

Comparing Class II MHC DRB3 Diversity in Colombian Simmental and Simbrah Cattle Across Worldwide Bovine Populations

The major histocompatibility complex (MHC) exerts great influence on responses to infectious diseases and vaccination due to its fundamental role in the adaptive immune system. Knowledge about MHC polymorphism distribution among breeds can provide insights into cattle evolution and diversification as well as population-based immune response variability, thus guiding further studies. Colombian Simmental and Simbrah cattle’s BoLA-DRB3 genetic diversity was compared to that of taurine and zebuine breeds worldwide to estimate functional diversity. High allele richness was observed for Simmental and Simbrah cattle; nevertheless, high homozygosity was associated with individual low sequence variability in both the β1 domain and the peptide binding region (PBR), thereby implying reduced MHC-presented peptide repertoire size. There were strong signals of positive selection acting on BoLA-DRB3 in all populations, some of which were poorly structured and displayed common alleles accounting for their high genetic similarity. PBR sequence correlation analysis suggested that, except for a few populations exhibiting some divergence at PBR, global diversity regarding potential MHC-presented peptide repertoire could be similar for the cattle populations analyzed here, which points to the retention of functional diversity in spite of the selective pressures imposed by breeding.

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

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

A pooled testing system to rapidly identify cattle carrying the elite controller BoLA-DRB3*009:02 haplotype against bovine leukemia virus infection

As genetically resistant individuals, the “elite controllers” (ECs) of human immunodeficiency virus infection have been focused on as the keys to developing further functional treatments in medicine. In the livestock production field, identifying the ECs of bovine leukemia virus (BLV) infection in cattle is desired to stop BLV transmission chains on farms. Cattle carrying the bovine leukocyte antigen (BoLA)‐DRB3*009:02 allele (DRB3*009:02) have a strong possibility of being BLV ECs. Most of cattle carrying this allele maintain undetectable BLV proviral loads and do not shed virus even when infected. BLV ECs can act as transmission barriers when placed between uninfected and infected cattle in a barn. To identify cattle carrying DRB3*009:02 in large populations more easily, we developed a pooled testing system. It employs a highly sensitive, specific real‐time PCR assay and TaqMan MGB probes (DRB3*009:02‐TaqMan assay). Using this system, we determined the percentage of DRB3*009:02‐carrying cattle on Kyushu Island, Japan. Our pooled testing system detected cattle carrying the DRB3*009:02 allele from a DNA pool containing one DRB3*009:02‐positive animal and 29 cattle with other alleles. Its capacity is sufficient for herd‐level screening for DRB3*009:02‐carrying cattle. The DRB3*009:02‐TaqMan assay showed high‐discriminative sensitivity and specificity toward DRB3*009:02, making it suitable for identifying DRB3*009:02‐carrying cattle in post‐screening tests on individuals. We determined that the percentage of DRB3*009:02‐carrying cattle in Kyushu Island was 10.56%. With its ease of use and reliable detection, this new method strengthens the laboratory typing for DRB3*009:02‐carrying cattle. Thus, our findings support the use of BLV ECs in the field.

Nanovaccines against Animal Pathogens: The Latest Findings

Nowadays, safe and efficacious vaccines represent powerful and cost-effective tools for global health and economic growth. In the veterinary field, these are undoubtedly key tools for improving productivity and fighting zoonoses. However, cases of persistent infections, rapidly evolving pathogens having high variability or emerging/re-emerging pathogens for which no effective vaccines have been developed point out the continuing need for new vaccine alternatives to control outbreaks. Most licensed vaccines have been successfully used for many years now; however, they have intrinsic limitations, such as variable efficacy, adverse effects, and some shortcomings. More effective adjuvants and novel delivery systems may foster real vaccine effectiveness and timely implementation. Emerging vaccine technologies involving nanoparticles such as self-assembling proteins, virus-like particles, liposomes, virosomes, and polymeric nanoparticles offer novel, safe, and high-potential approaches to address many vaccine development-related challenges. Nanotechnology is accelerating the evolution of vaccines because nanomaterials having encapsulation ability and very advantageous properties due to their size and surface area serve as effective vehicles for antigen delivery and immunostimulatory agents. This review discusses the requirements for an effective, broad-coverage-elicited immune response, the main nanoplatforms for producing it, and the latest nanovaccine applications for fighting animal pathogens.

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

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

Immunoinformatics Approach to Design Multi-Epitope- Subunit Vaccine against Bovine Ephemeral Fever Disease

Bovine ephemeral fever virus (BEFV) is an overlooked pathogen, recently gaining widespread attention owing to its associated enormous economic impacts affecting the global livestock industries. High endemicity with rapid spread and morbidity greatly impacts bovine species, demanding adequate attention towards BEFV prophylaxis. Currently, a few suboptimum vaccines are prevailing, but were confined to local strains with limited protection. Therefore, we designed a highly efficacious multi-epitope vaccine candidate targeted against the geographically distributed BEFV population. By utilizing immunoinformatics technology, all structural proteins were targeted for B- and T-cell epitope prediction against the entire allele population of BoLA molecules. Prioritized epitopes were adjoined by linkers and adjuvants to effectively induce both cellular and humoral immune responses in bovine. Subsequently, the in silico construct was characterized for its physicochemical parameters, high immunogenicity, least allergenicity, and non-toxicity. The 3D modeling, refinement, and validation of ligand (vaccine construct) and receptor (bovine TLR7) then followed molecular docking and molecular dynamic simulation to validate their stable interactions. Moreover, in silico cloning of codon-optimized vaccine construct in the prokaryotic expression vector (pET28a) was explored. This is the first time HTL epitopes have been predicted using bovine datasets. We anticipate that the designed construct could be an effective prophylactic remedy for the BEF disease that may pave the way for future laboratory experiments.

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

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

Genetic assessment of BoLA-DRB3 polymorphisms by comparing Bangladesh, Ethiopian, and Korean cattle

Attributable to their major function in pathogen recognition, the use of bovine leukocyte antigens (BoLA) as disease markers in immunological traits in cattle is well established. However, limited report exists on polymorphism of the BoLA gene in zebu cattle breeds by high resolution typing methods. Thus, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 100 animals (Boran, n = 13; Sheko, n = 20; Fogera, n = 16; Horro, n = 19), Hanwoo cattle (n = 18) and Bangladesh Red Chittagong zebu (n = 14). Out of the 59 detected alleles, 43 were already deposited under the Immuno Polymorphism Database for major histocompatibility complex (IPD-MHC) while 16 were unique to this study. Assessment of the level of genetic variability at the population and sequence levels with genetic distance in the breeds considered in this study showed that Zebu breeds had a gene diversity score greater than 0.752, nucleotide diversity score greater than 0.152, and mean number of pairwise differences higher than 14, being very comparable to those investigated for other cattle breeds. Regarding neutrality tests analyzed, we investigated that all the breeds except Hanwoo had an excess number of alleles and could be expected from a recent population expansion or genetic hitchhiking. Howbeit, the observed heterozygosity was not significantly (p < 0.05) higher than the expected heterozygosity. The Hardy Weinberg equilibrium (HWE) analysis revealed non-significant excess of heterozygote animals, indicative of plausible over-dominant selection. The pairwise FST values suggested a low genetic variation among all the breeds (FST = 0.056; p < 0.05), besides the rooting from the evolutionary or domestication history of the cattle. No detached clade was observed in the evolutionary divergence study of the BoLA-DRB3 gene, inferred from the phylogenetic tree based on the maximum likelihood model. The investigation herein indicated the clear differences in BoLA-DRB3 gene variability between African and Asian cattle breeds.

A comparative analysis of SLA-DRB1 genetic diversity in Colombian (creoles and commercial line) and worldwide swine populations

Analysing pig class II mayor histocompatibility complex (MHC) molecules is mainly related to antigen presentation. Identifying frequently-occurring alleles in pig populations is an important aspect to be considered when developing peptide-based vaccines. Colombian creole pig populations have had to adapt to local conditions since entering Colombia; a recent census has shown low amounts of pigs which is why they are considered protected by the Colombian government. Commercial hybrids are more attractive regarding production. This research has been aimed at describing the allele distribution of Colombian pigs from diverse genetic backgrounds and comparing Colombian SLA-DRB1 locus diversity to that of internationally reported populations. Twenty SLA-DRB1 alleles were identified in the six populations analysed here using sequence-based typing. The amount of alleles ranged from six (Manta and Casco Mula) to nine (San Pedreño). Only one allele (01:02) having > 5% frequency was shared by all three commercial line populations. Allele 02:01:01 was shared by five populations (around > 5% frequency). Global F_ST indicated that pig populations were clearly structured, as 20.6% of total allele frequency variation was explained by differences between populations (F_ST = 0.206). This study’s results confirmed that the greatest diversity occurred in wild boars, thereby contrasting with low diversity in domestic pig populations.

Lipopolysaccharide triggers different transcriptional signatures in taurine and indicine cattle macrophages: Reactive oxygen species and potential outcomes to the development of immune response to infections

Macrophages are classified upon activation as classical activated M1 and M2 anti-inflammatory regulatory populations. This macrophage polarization is well characterized in humans and mice, but M1/M2 profile in cattle has been far less explored. Bos primigenius taurus (taurine) and Bos primigenius indicus (indicine) cattle display contrasting levels of resistance to infection and parasitic diseases such as C57BL/6J and Balb/c murine experimental models of parasite infection outcomes based on genetic background. Thus, we investigated the differential gene expression profile of unstimulated and LPS stimulated monocyte-derived macrophages (MDMs) from Holstein (taurine) and Gir (indicine) breeds using RNA sequencing methodology. For unstimulated MDMs, the contrast between Holstein and Gir breeds identified 163 Differentially Expressed Genes (DEGs) highlighting the higher expression of C-C chemokine receptor type five (CCR5) and BOLA-DQ genes in Gir animals. LPS-stimulated MDMs from Gir and Holstein animals displayed 1,257 DEGs enriched for cell adhesion and inflammatory responses. Gir MDMs cells displayed a higher expression of M1 related genes like Nitric Oxide Synthase 2 (NOS2), Toll like receptor 4 (TLR4), Nuclear factor NF-kappa-B 2 (NFKB2) in addition to higher levels of transcripts for proinflammatory cytokines, chemokines, complement factors and the acute phase protein Serum Amyloid A (SAA). We also showed that gene expression of inflammatory M1 population markers, complement and SAA genes was higher in Gir in buffy coat peripheral cells in addition to nitric oxide concentration in MDMs supernatant and animal serum. Co-expression analyses revealed that Holstein and Gir animals showed different transcriptional signatures in the MDMs response to LPS that impact on cell cycle regulation, leukocyte migration and extracellular matrix organization biological processes. Overall, the results suggest that Gir animals show a natural propensity to generate a more pronounced M1 inflammatory response than Holstein, which might account for a faster immune response favouring resistance to many infection diseases.

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

Background

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

Methods

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

Results

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

Conclusion

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