Characterization of bovine MHC class II DRB3 diversity in South American Holstein cattle 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.

BoLA-DRB3 Polymorphism Associated with Bovine Leukemia Virus Infection and Proviral Load in Holstein Cattle in Egypt

Bovine leukemia virus (BLV) is the etiological agent of enzootic bovine leukosis, the most prevalent neoplastic disease of cattle worldwide. The immune response to BLV and disease susceptibility and resistance in cattle are strongly correlated with the bovine leukocyte antigen (BoLA)-DRB3 allelic polymorphism. BLV infection continues to spread in Egypt, in part because the relationships between BLV infection, proviral load in Egypt, and BoLA-DRB3 polymorphism are unknown. Here, we identified 18 previously reported alleles in 121 Holstein cows using a polymerase chain reaction sequence-based typing method. Furthermore, BoLA-DRB3 gene polymorphisms in these animals were investigated for their influence on viral infection. BoLA-DRB3*015:01 and BoLA-DRB3*010:01 were identified as susceptible and resistant alleles, respectively, for BLV infection in the tested Holsteins. In addition, BoLA-DRB3*012:01 was associated with low PVL in previous reports but high PVL in Holstein cattle in Egypt. This study is the first to demonstrate that the BoLA-DRB3 polymorphism confers resistance and susceptibility to PVL and infections of BLV in Holstein cattle in Egypt. Our results can be useful for the disease control and eradication of BLV through genetic selection.

Identifying Pathogen and Allele Type Simultaneously in a Single Well Using Droplet Digital PCR

In the transmission control of chronic and untreatable livestock diseases such as bovine leukemia virus (BLV) infection, the removal of viral superspreaders is a fundamental approach. On the other hand, selective breeding of cattle with BLV-resistant capacity is also critical for reducing the viral damage to productivity by keeping infected cattle. To provide a way of measuring BLV proviral load (PVL) and identifying susceptible/resistant cattle simply and rapidly, we developed a fourplex droplet digital PCR method targeting the BLV pol gene, BLV-susceptible bovine major histocompatibility complex (BoLA)-DRB3*016:01 allele, resistant DRB3*009:02 allele, and housekeeping RPP30 gene (IPATS-BLV). IPATS-BLV successfully measured the percentage of BLV-infected cells and determined allele types precisely. Furthermore, it discriminated homozygous from heterozygous carriers. Using this method to determine the impact of carrying these alleles on the BLV PVL, we found DRB3*009:02-carrying cattle could suppress the PVL to a low or undetectable level, even with the presence of a susceptible heterozygous allele. Although the population of DRB3*016:01-carrying cattle showed significantly higher PVLs compared with cattle carrying other alleles, their individual PVLs were highly variable. Because of the simplicity and speed of this single-well assay, our method has the potential of being a suitable platform for the combined diagnosis of pathogen level and host biomarkers in other infectious diseases satisfying the two following characteristics of disease outcomes: (i) pathogen level acts as a critical maker of disease progression; and (ii) impactful disease-related host genetic biomarkers are already identified. IMPORTANCE While pathogen-level quantification is an important diagnostic of disease severity and transmissibility, disease-related host biomarkers are also useful in predicting outcomes in infectious diseases. In this study, we demonstrate that combined proviral load (PVL) and host biomarker diagnostics can be used to detect bovine leukemia virus (BLV) infection, which has a negative economic impact on the cattle industry. We developed a fourplex droplet digital PCR assay for PVL of BLV and susceptible and resistant host genes named IPATS-BLV. IPATS-BLV has inherent merits in measuring PVL and identifying susceptible and resistant cattle with superior simplicity and speed because of a single-well assay. Our new laboratory technique contributes to strengthening risk-based herd management used to control within-herd BLV transmission. Furthermore, this assay design potentially improves the diagnostics of other infectious diseases by combining the pathogen level and disease-related host genetic biomarker to predict disease outcomes.

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.

A novel missense mutation (rs464874590) within BoLA-DOB gene associated with the heat-resistance in Chinese cattle

Bovine leukocyte antigen, class II, DO beta (BoLA-DOB) is related to antigen presentation, which can triggered by multicul factors. And the condition of immune function determines how much cattle load to heat stress. To evaluate the relationship between heat-resistance and single nucleotide polymorphisms (SNPs) in BoLA-DOB gene, our study has taken further analysis in Chinese indigenous cattle for the first time. A missense single nucleotide polymorphism (rs464874590) was detected in BoLA-DOB gene. We directly sequenced rs464874590 (NM_001013600.1 g.7122762 A > G) in BoLA-DOB gene of 522 individuals of 26 cattle breeds. The frequency of allele G gradually decreases from south to north with distinct climatic distribution characteristics. Further association analysis was carried out between different genotypes and environmental parameters, including annual mean temperature (T), relative humidity (RH), and temperature-humidity index (THI). The result showed that three genotypes were significantly correlated with T, H, and THI (P < 0.01), indicating that GG genotype was distributed in areas with hot and moist conditions. Therefore, our results suggested that the rs464874590 could be applied as a genetic marker to detect the heat-resistance of Chinese indigenous cattle.

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.

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

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

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

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

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

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

Genetic Variation and Population Differentiation in the Bovine Lymphocyte Antigen DRB3.2 Locus of South African Nguni Crossbred Cattle

Simple Summary

Indigenous cattle breeds are important to their native environment as they confer significant and often unique adaptability traits. The Nguni is one such breeds that is indigenous to the Southern African region. This breed confers resistance to thermal stress and diseases, amongst other factors. The bovine major histocompatibility complex is an important region, which codes for alleles that have been associated with a plethora of diseases. In the current study, the genetic diversity within this region was assessed in Nguni crossbred cattle. This was done to detect the gene pool of the Nguni breed, and to identify genes that might be important within this breed. The populations displayed a high degree of genetic diversity, and some alleles were common throughout the populations and accounted for a significant portion of the total alleles. This high genetic diversity could account for the great adaptability of the Nguni breed to Southern Africa.

Abstract

The bovine lymphocyte antigen (BoLA-DRB3) gene is an important region that codes for glycoproteins responsible for the initiation of an immune response. BoLA-DRB3 alleles have been demonstrated to be associated with disease resistance/tolerance. Therefore, great genetic diversity is correlated with better adaptation, fitness, and robustness. The current study was conducted to assess the population genetic structure of the BoLA-DRB3 gene in Nguni crossbred cattle using polymerase chain reaction-sequence based typing (PCR-SBT). High genetic diversity was detected, with 30 alleles, 11 of which are novel to the study. Alleles DRB3*0201, DRB3*0701, DRB*0901, and DRB*1601 were present in all populations and accounted for nearly around 50% of all observed alleles. A mean genetic diversity (H_E) of 0.93 was detected. The high overall genetic diversity is possibly associated with pathogen-assisted selection and heterozygote advantage. Such high diversity might explain the hardiness of the Nguni crossbred cattle to the Southern African region. Low population genetic structure was identified (F_ST = 0.01), suggesting possible gene flow between populations and retention of similar alleles. The study was undertaken to bridge the dearth of such studies in South African breeds and it is imperative for effective sustainability of indigenous breeds and the implementation of effective breeding strategies.

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

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

Genetic variability and biodiversity of Ukrainian Gray cattle by the BoLA-DRB3 gene

At the current stage of genetic studies of cattle, more and more attention is being drawn to autochthonous breeds. Native cattle have a number of prominent phenotypic traits and have preserved unique genes and their combinations lost by modern commercial breeds, which would be valuable to use in selective programs. We surveyed polymorphism of the Ukrainian autochthonous Gray breed according to alleles of exon 2 of the BoLA-DRB3 gene. The uniqueness of the gene lies in the broad variability of its allele variants. Significant informativeness at DNA level is quite important for genetic studies. We surveyed allele polymorphism using the PCR-RLFP method on DNA isolated from 88 samples of blood of cows and 5 samples of sperm. We identified 28 alleles, of which 23 variants were nomenclature ones and 5 (jba, *jab, *jbb, *nad and *nda) were “without established nomenclature”, their share accounting for 8.9%. Four alleles *06, *12, *16 and *jba had a frequency above 5% and occupied 69.9% of the breed’s allele fund overall. The commonest allele was BoLA-DRB3.2*16 (44.1%). In total, we found 40 genotypes. Considering the significant dominance of variant *16, as expected, 5 genotypes with its inclusion occurred: *16/*16, *12/*16, *06/*16, *16/*24 and *jba/*16. It was present in the genotype of two out three studied animals. Parameters of heterozygosity, effective number of alleles, Shannon and Pielou indices indicate that Ukrainian Gray cattle are characterized by lowest level of genetic variability and biodiversity according to the BoLA-DRB3 gene compared with other breeds. Due to significant dominance of allele *16, the breed has no inbred motifs. We noted deviation toward increase in homozygosity without deviations from the norm of the distribution according to Hardy-Weinberg equilibrium. The obtained results will be used for genetic-populational programs with the purpose of improving the genetic potential of cattle breeds in terms of economically beneficial traits and diseases of cattle.

BoLA-DRB3 genetic diversity in Highland Creole cattle from Bolivia

The genetic diversity of the BoLA-DRB3 gene has been reported in different cattle breeds owing to its central role in the immune response. However, it is still unknown in hundreds of cattle breeds, especially native populations. Here, we studied BoLA-DRB3 genetic diversity in Highland Creole cattle (CrAl) from Western Bolivia, raised at altitudes between 3800 and 4200 m. DNAs from 48 CrAl cattle were genotyped for BoLA-DRB3 exon 2 alleles using polymerase chain reaction-sequence-based typing (PCR-SBT). The results were compared with 1341 previously reported data from Tropical Creole cattle and other breeds raised in the region. Twenty-three BoLA-DRB3 alleles were identified in CrAl, including the BoLA-DRB3*029:02 variant previously detected in other Creole cattle. Observed and expected heterozygosity were 0.87 and 0.93, respectively. Nucleotide diversity and the number of pairwise difference values were 0.078 and 19.46, respectively. The average number of nonsynonymous and synonymous substitutions were 0.037 and 0.097 for the entire BoLA-DRB3 exon 2, and 0.129 and 0.388 for the antigen-binding site, respectively. Venn analysis and the review of the IPD-MHC database and the literature showed that 2 of 64 alleles were only detected in CrAl, including BoLA-DRB3*029:01 previously reported in African cattle and *048:01 detected in Philippine cattle. Two additional alleles, BoLA-DRB3*007:02 and *029:02, were only present in CrAl and Lowland Creole cattle. Principal Component Analysis (PCA) showed that Bolivian Creole cattle breeds were closely located but they were distant from the Colombian Hartón del Valle Creole. F_ST analysis showed a low degree of genetic differentiation between Highland and Lowland Bolivian Creole cattle (F_ST = 0.015). The present results contribute to increasing our knowledge of BoLA-DRB3 genetic diversity in cattle breeds.

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.

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

Bovine leukocyte antigens (BoLA) have been used as disease markers and immunological traits in cattle due to their primary role in pathogen recognition by the immune system. A higher MHC allele diversity in a population will allow presenting a broader peptide repertoire. However, loss of overall diversity due to domestication process can decrease a population's peptide repertoire. Within the context of zebu and taurine cattle populations, BoLA-DRB3 genetic diversity in Spanish Morucha and Colombian Normande cattle was analyzed and an approach to estimate functional diversity was performed. Sequence-based typing was used for identifying 29, 23, 27, and 28 alleles in Spanish Morucha, Nariño-, Boyacá-, and Cundinamarca-Normande cattle, respectively. These breeds had remarkably low heterozygosity levels and the Hardy-Weinberg principle revealed significant heterozygote deficiency. F_ST and D_A genetic distance showed that Colombian Normande populations had greater variability than other phenotypically homogeneous breeds, such as Holstein. It was also found that Spanish Morucha cattle were strongly differentiated from other cattle breeds. Spanish Morucha had greater divergence in the peptide-binding region regarding other cattle breeds. However, peptide-binding region covariation indicated that the potential peptide repertoire seemed equivalent among cattle breeds. Despite the genetic divergence observed, the extent of the potential peptide repertoire in the cattle populations studied appears to be similar and thus their pathogen recognition potential should be equivalent, suggesting that functional diversity might persist in the face of bottlenecks imposed by domestication and breeding.

Genetic Diversity of Bovine Major Histocompatibility Complex Class II DRB3 locus in cattle breeds from Asia compared to those from Africa and America

Genetic polymorphisms and diversity of BoLA-DRB3.2 are essential because of DRB3 gene's function in innate immunity and its association with infectious diseases resistance or tolerance in cattle. The present study was aimed at assessing the level of genetic diversity of DRB3 in the exon 2 (BoLA-DRB3.2) region in African, American and Asian cattle breeds. Amplification of exon 2 in 174 cattle revealed 15 haplotypes. The breeds with the highest number of haplotypes were Brangus (10), Sokoto Gudali (10) and Dajal (9), while the lowest number of haplotypes were found in Holstein and Sahiwal with 4 haplotypes each. Medium Joining network obtained from haplotypic data showed that all haplotypes condensed around a centric area and each sequence (except in H-3, H-51 and H-106) representing almost a specific haplotype. The BoLA-DRB3.2 sequence analyses revealed a non-significant higher rate of non-synonymous (dN) compared to synonymous substitutions (dS). The ratio of dN/dS substitution across the breeds were observed to be greater than one suggesting that variation at the antigen-binding sites is under positive selection; thus increasing the chances of these breeds to respond to wide array of pathogenic attacks. An analysis of molecular variance revealed that 94.01 and 5.99% of the genetic variation was attributable to differences within and among populations, respectively. Generally, results obtained suggest that within breed genetic variation across breeds is higher than between breeds. This genetic information will be important for investigating the relationship between BoLADRB3.2 and diseases in various cattle breeds studied with attendant implication on designing breeding programs that will aim at selecting individual cattle that carry resistant alleles.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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