Short communication: characterization of DRB3 alleles in the MHC of Japanese shorthorn cattle by polymerase chain reaction-sequence-based typing

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

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

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.

Genome-Wide Association Studies for Milk Somatic Cell Score in Romanian Dairy Cattle

Mastitis is one of the most frequently encountered diseases in dairy cattle, negatively affecting animal welfare and milk production. For this reason, contributions to understanding its genomic architecture are of great interest. Genome-wide association studies (GWAS) have identified multiple loci associated with somatic cell score (SCS) and mastitis in cattle. However, most of the studies have been conducted in different parts of the world on various breeds, and none of the investigations have studied the genetic architecture of mastitis in Romanian dairy cattle breeds up to this point in time. In this study, we report the first GWAS for SCS in dairy cattle breeds from Romania. For GWAS, we used an Axiom Bovine v3 SNP-chip (>63,000 Single Nucleotide Polymorphism -SNPs) and 33,330 records from 690 cows belonging to Romanian Spotted (RS) and Romanian Brown (RB) cattle. The results found one SNP significantly associated with SCS in the RS breed and 40 suggestive SNPs with -log₁₀ (p) from 4 to 4.9 for RS and from 4 to 5.4 in RB. From these, 14 markers were located near 12 known genes (AKAP8, CLHC1, MEGF10, SATB2, GATA6, SPATA6, COL12A1, EPS8, LUZP2, RAMAC, IL12A and ANKRD55) in RB cattle, 3 markers were close to ZDHHC19, DAPK1 and MMP7 genes, while one SNP overlapped the HERC3 gene in RS cattle. Four genes (HERC3, LUZP2, AKAP8 and MEGF10) associated with SCS in this study were previously reported in different studies. The most significant SNP (rs110749552) associated with SCS was located within the HERC3 gene. In both breeds, the SNPs and position of association signals were distinct among the three parities, denoting that mastitis is controlled by different genes that are dependent according to parity. The current results contribute to an expansion in the body of knowledge regarding the proportion of genetic variability explained by SNPs for SCS in dairy cattle.

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.

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.

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.

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.

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.

BoLA-DRB3 Polymorphism is Associated with Differential Susceptibility to Bovine Leukemia Virus-Induced Lymphoma and Proviral Load

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leucosis. However, less than 5% of BLV-infected cattle will develop lymphoma, suggesting that, in addition to viral infection, host genetic polymorphisms might play a role in disease susceptibility. Bovine leukocyte antigen (BoLA)-DRB3 is a highly polymorphic gene associated with BLV proviral load (PVL) susceptibility. Due to the fact that PVL is positively associated with disease progression, it is believed that controlling PVL can prevent lymphoma development. Thus, many studies have focused on the relationship between PVL and BoLA-DRB3. Despite this, there is little information regarding the relationship between lymphoma and BoLA-DRB3. Furthermore, whether or not PVL-associated BoLA-DRB3 is linked to lymphoma-associated BoLA-DRB3 has not been clarified. Here, we investigated whether or not lymphoma-associated BoLA-DRB3 is correlated with PVL-associated BoLA-DRB3. We demonstrate that two BoLA-DRB3 alleles were specifically associated with lymphoma resistance (*010:01 and *011:01), but no lymphoma-specific susceptibility alleles were found; furthermore, two other alleles, *002:01 and *012:01, were associated with PVL resistance and susceptibility, respectively. In contrast, lymphoma and PVL shared two resistance-associated (DRB3*014:01:01 and *009:02) BoLA-DRB3 alleles. Interestingly, we found that PVL associated alleles, but not lymphoma associated alleles, are related with the anti-BLV gp51 antibody production level in cows. Overall, our study is the first to demonstrate that the BoLA-DRB3 polymorphism confers differential susceptibility to BLV-induced lymphoma and PVL.

Comparison of cattle BoLA-DRB3 typing by PCR-RFLP, direct sequencing, and high-resolution DNA melting curve analysis

Major histocompatibility complex (MHC) represents an important genetic marker for manipulation to improve the health and productivity of cattle. It is closely associated with numerous disease susceptibilities and immune responses. Bovine MHC, also called bovine leukocyte antigen (BoLA), is considered as a suitable marker for genetic diversity studies. In cattle, most of the polymorphisms are located in exon 2 of BoLA-DRB3, which encodes the peptide-binding cleft. In this study, the polymorphism of the BoLA-DRB3.2 gene in Holstein's calves was studied using high resolution melting curve analysis (HRM). Observed HRM results were compared to PCR-RFLP and direct sequencing techniques. Eight different HRM and seven different RFLP profiles were identified among the population studied. By comparing to sequencing data, HRM could completely discriminate all genotypes (eight profiles), while the RFLP failed to distinguish between the genotypes *1101/*1001 and *1104/*1501. According to the results, the HRM analysis method gave more accurate results than RFLP by differentiating between the BoLA-DRB3.2 genotypes. Due to the Co-dominant nature of the MHC alleles, HRM technique could be used for investigating the polymorphisms of genotypes and their associations with immune responses.

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.

Invited review: Microbiota of the bovine udder: Contributing factors and potential implications for udder health and mastitis susceptibility

Various body sites of vertebrates provide stable and nutrient-rich ecosystems for a diverse range of commensal, opportunistic, and pathogenic microorganisms to thrive. The collective genomes of these microbial symbionts (the microbiome) provide host animals with several advantages, including metabolism of indigestible carbohydrates, biosynthesis of vitamins, and modulation of innate and adaptive immune systems. In the context of the bovine udder, however, the relationship between cow and microbes has been traditionally viewed strictly from the perspective of host-pathogen interactions, with intramammary infections by mastitis pathogens triggering inflammatory responses (i.e., mastitis) that are often detrimental to mammary tissues and cow physiology. This traditional view has been challenged by recent metagenomic studies indicating that mammary secretions of clinically healthy quarters can harbor genomic markers of diverse bacterial groups, the vast majority of which have not been associated with mastitis. These observations have given rise to the concept of "commensal mammary microbiota," the ecological properties of which can have important implications for understanding the pathogenesis of mastitis and offer opportunities for development of novel prophylactic or therapeutic products (or both) as alternatives to antimicrobials. Studies conducted to date have suggested that an optimum diversity of mammary microbiota is associated with immune homeostasis, whereas the microbiota of mastitic quarters, or those with a history of mastitis, are considerably less diverse. Whether disruption of the diversity of udder microbiota (dysbiosis) has a role in determining mastitis susceptibility remains unknown. Moreover, little is known about contributions of various biotic and abiotic factors in shaping overall diversity of udder microbiota. This review summarizes current understanding of the microbiota within various niches of the udder and highlights the need to view the microbiota of the teat apex, teat canal, and mammary secretions as interconnected niches of a highly dynamic microbial ecosystem. In addition, host-associated factors, including physiological and anatomical parameters, as well as genetic traits that may affect the udder microbiota are briefly discussed. Finally, current understanding of the effect of antimicrobials on the composition of intramammary microbiota is discussed, highlighting the resilience of udder microbiota to exogenous perturbants.

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.

Characterization of bovine MHC class II DRB3 diversity in South American Holstein cattle populations

Holstein cattle dominate the global milk production industry because of their outstanding milk production, however, this breed is susceptible to tropical endemic pathogens and suffers from heat stress and thus fewer Holstein populations are raised in tropical areas. The bovine major histocompatibility complex (BoLA)-DRB3 class II gene is used as a marker for disease and immunological traits, and its polymorphism has been studied extensively in Holstein cattle from temperate and cold regions. We studied the genetic diversity of the BoLA-DRB3 gene in South American Holstein populations to determine whether tropical populations have diverged from those bred in temperate and cold regions by selection and/or crossbreeding with local native breeds. We specifically studied Exon 2 of this gene from 855 South American Holstein individuals by a polymerase chain reaction (PCR) sequence-based typing method. We found a high degree of gene diversity at the allelic (Na > 20 and He > 0.87) and molecular (π > 0.080) levels, but a low degree of population structure (FST = 0.009215). A principal components analysis and tree showed that the Bolivian subtropical population had the largest genetic divergence compared with Holsteins bred in temperate or cold regions, and that this population was closely related to Bolivian Creole cattle. Our results suggest that Holstein genetic divergence can be explained by selection and/or gene introgression from local germplasms. This is the first examination of BoLA-DRB3 in Holsteins adapted to tropical environments, and contributes to an ongoing effort to catalog bovine MHC allele frequencies by breed and location.

Genetic structure and relationships of 16 Asian and European cattle populations using DigiTag2 assay

In this study, we genotyped 117 autosomal single nucleotide polymorphisms using a DigiTag2 assay to assess the genetic diversity, structure and relationships of 16 Eurasian cattle populations, including nine cattle breeds and seven native cattle. Phylogenetic and principal component analyses showed that Bos taurus and Bos indicus populations were clearly distinguished, whereas Japanese Shorthorn and Japanese Polled clustered with European populations. Furthermore, STRUCTURE analysis demonstrated the distinct separation between Bos taurus and Bos indicus (K=2), and between European and Asian populations (K=3). In addition, Japanese Holstein exhibited an admixture pattern with Asian and European cattle (K=3‐5). Mongolian (K=13‐16) and Japanese Black (K=14‐16) populations exhibited admixture patterns with different ancestries. Bos indicus populations exhibited a uniform genetic structure at K=2‐11, thereby suggesting that there are close genetic relationships among Bos indicus populations. However, the Bhutan and Bangladesh populations formed a cluster distinct from the other Bos indicus populations at K=12‐16. In conclusion, our study could sufficiently explain the genetic construction of Asian cattle populations, including: (i) the close genetic relationships among Bos indicus populations; (ii) the genetic influences of European breeds on Japanese breeds; (iii) the genetic admixture in Japanese Holstein, Mongolian and Japanese Black cattle; and (iv) the genetic subpopulations in Southeast Asia.

Assessment of biodiversity in Chilean cattle using the distribution of major histocompatibility complex class II BoLA-DRB3 allele

Bovine leukocyte antigens (BoLAs) are used extensively as markers for bovine disease and immunological traits. In this study, we estimated BoLA-DRB3 allele frequencies using 888 cattle from 10 groups, including seven cattle breeds and three crossbreeds: 99 Red Angus, 100 Black Angus, 81 Chilean Wagyu, 49 Hereford, 95 Hereford × Angus, 71 Hereford × Jersey, 20 Hereford × Overo Colorado, 113 Holstein, 136 Overo Colorado, and 124 Overo Negro cattle. Forty-six BoLA-DRB3 alleles were identified, and each group had between 12 and 29 different BoLA-DRB3 alleles. Overo Negro had the highest number of alleles (29); this breed is considered in Chile to be an 'Old type' European Holstein Friesian descendant. By contrast, we detected 21 alleles in Holstein cattle, which are considered to be a 'Present type' Holstein Friesian cattle. Chilean cattle groups and four Japanese breeds were compared by neighbor-joining trees and a principal component analysis (PCA). The phylogenetic tree showed that Red Angus and Black Angus cattle were in the same clade, crossbreeds were closely related to their parent breeds, and Holstein cattle from Chile were closely related to Holstein cattle in Japan. Overall, the tree provided a thorough description of breed history. It also showed that the Overo Negro breed was closely related to the Holstein breed, consistent with historical data indicating that Overo Negro is an 'Old type' Holstein Friesian cattle. This allelic information will be important for investigating the relationship between major histocompatibility complex (MHC) and disease.

The great diversity of major histocompatibility complex class II genes in Philippine native cattle

Bovine leukocyte antigens (BoLA) are extensively used as markers for bovine disease and immunological traits. However, none of the BoLA genes in Southeast Asian breeds have been characterized by polymerase chain reaction (PCR)-sequence-based typing (SBT). Therefore, we sequenced exon 2 of the BoLA class II DRB3 gene from 1120 individual cows belonging to the Holstein, Sahiwal, Simbrah, Jersey, Brahman, and Philippine native breeds using PCR-SBT. Several cross-breeds were also examined. BoLA-DRB3 PCR-SBT identified 78 previously reported alleles and five novel alleles. The number of BoLA-DRB3 alleles identified in each breed from the Philippines was higher (71 in Philippine native cattle, 58 in Brahman, 46 in Holstein × Sahiwal, and 57 in Philippine native × Brahman) than that identified in breeds from other countries (e.g., 23 alleles in Japanese Black and 35 in Bolivian Yacumeño cattle). A phylogenetic tree based on the D_A distance calculated from the BoLA-DRB3 allele frequency showed that Philippine native cattle from different Philippine islands are closely related, and all of them are closely similar to Philippine Brahman cattle but not to native Japanese and Latin American breeds. Furthermore, the BoLA-DRB3 allele frequency in Philippine native cattle from Luzon Island, located in the Northern Philippines was different from that in cattle from Iloilo, Bohol, and Leyte Islands, which are located in the Southern Philippines. Therefore, we conclude that Philippine native cattle can be divided into two populations, North and South areas. Moreover, a neutrality test revealed that Philippine native cattle from Leyte showed significantly greater genetic diversity, which may be maintained by balancing selection. This study shows that Asian breeds have high levels of BoLA-DRB3 polymorphism. This finding, especially the identification of five novel BoLA-DRB3 alleles, will be helpful for future SBT studies of BoLA-DRB3 alleles in East Asian cattle.

The major histocompatibility complex in bovines: a review

Productivity in dairy cattle and buffaloes depends on the genetic factors governing the production of milk and milk constituents as well as genetic factors controlling disease resistance or susceptibility. The immune system is the adaptive defense system that has evolved in vertebrates to protect them from invading pathogens and also carcinomas. It is remarkable in the sense that it is able to generate an enormous variety of cells and biomolecules which interact with each other in numerous ways to form a complex network that helps to recognize, counteract, and eliminate the apparently limitless number of foreign invading pathogens/molecules. The major histocompatibility complex which is found to occur in all mammalian species plays a central role in the development of the immune system. It is an important candidate gene involved in susceptibility/resistance to various diseases. It is associated with intercellular recognition and with self/nonself discrimination. It plays major role in determining whether transplanted tissue will be accepted as self or rejected as foreign.

Identification and diversity of bovine major histocompatibility complex class II haplotypes in Japanese Black and Holstein cattle in Japan

Bovine leukocyte antigen (BoLA), the major histocompatibility complex of cattle, is one of the most polymorphic gene clusters. We genotyped a population of 109 Japanese Black and 39 Holstein cattle to analyze their BoLA class II haplotypes, BoLA-DRB3 locus, 5 BoLA-DQA loci, and 5 BoLA-DQB loci. We identified 26 previously reported DRB3 alleles, 22 previously reported and 3 new DQA alleles, and 24 previously reported and 6 new DQB alleles. A dendrogram was constructed based on the predicted amino acid sequences of the α1 or β1 domains encoded by BoLA-DQA or -DQB alleles, which revealed that DQA alleles were clustered into 5 loci, whereas DQB alleles could not be clearly assigned to specific DQB loci. The BoLA-DRB3-DQA-DQB haplotypes were sorted by sequential analytical processes, and 42 distinct haplotypes, including 11 previously published haplotypes and 31 novel haplotypes, were defined. Strong linkage disequilibrium was present in the BoLA genes. We also compared DRB3-DQA1 haplotype frequencies between 507 Japanese Black and 143 Holstein cattle. Thirty-nine DRB3-DQA1 haplotypes were identified, including 29 haplotypes from Japanese Black and 22 haplotypes from Holstein cattle. The majority of the haplotypes could be identified in both breeds, although several haplotypes were identified in only a single breed. This is the first report presenting a detailed study of the BoLA class II haplotype in Japanese Black and Holstein cattle in Japan.

Distribution of BoLA-DRB3 allelic frequencies and identification of two new alleles in Iranian buffalo breed

The role of the major histocompatibility complex (MHC) in the immune response makes it an attractive candidate gene for associations with disease resistance and susceptibility. This study describes genetic variability in the BoLA-DRB3 in Iranian buffaloes. Heminested PCR-RFLP method was used to identify the frequency of BoLA-DRB3 alleles. The BoLA-DRB3 locus is highly polymorphic in the study herd (12 alleles). Almost 63.50% of the alleles were accounted for by four alleles (BoLA-DRB3.2 ∗48, ∗20, ∗21, and obe) in Iranian buffalo. The DRB3.2 ∗48 allele frequency (24.20%) was higher than the others. The frequencies of the DRB3.2 ∗20 and DRB3.2 ∗21 are 14.52 and 14.00, respectively, and obe and gbb have a new pattern. Significant distinctions have been found between Iranian buffalo and other cattle breed studied. In the Iranian buffaloes studied alleles associated with resistance to various diseases are found.

Polymorphism of exon 2 of BoLA-DRB3 gene and its relationship with somatic cell score in Beijing Holstein cows

In the present study, the exon 2 and 3' end sequence of intron 1 of BoLA-DRB3 gene of 211 Beijing Holstein cows was amplified and a uniform fragment of 284 bp was obtained. The genetic polymorphism was investigated by PCR-RFLP using two restriction endonucleases HaeIII and BstYI. Seven genotypes were detected by digesting the PCR products with HaeIII. The frequency of AA, AB, AC, AD, BB, BC and BF genotypes was 0.4638, 0.0193, 0.0193, 0.3478, 0.0290, 0.0386 and 0.0822, respectively. Three genotypes were found by digesting the PCR products with BstYI. The frequency of AA, AB and BB genotypes was 0.0569, 0.2844 and 0.6587, respectively. The relationship between the polymorphisms in exon 2 of BoLA-DRB3 gene and somatic cell score (SCS) in Beijing Holstein cows was analyzed by least squares linear model. No significant difference was detected among least squares means of SCS for seven HaeIII-RFLP genotypes (P > 0.05). As for BstYI-RFLP analysis, least squares mean of SCS for AA was significantly lower than that for AB (P < 0.05) or BB (P < 0.05), and no significant difference was detected between AB and BB genotypes (P > 0.05). BstYI AA was the most favorable genotype and BstYI BB was the most unfavorable genotype for mastitis resistance. The information found in the present study is very important for improving mastitis resistance in dairy cattle by marker assisted selection.

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

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

Analysis of relationship between bovine lymphocyte antigen DRB3.2 alleles, somatic cell count and milk traits in Iranian Holstein population

The major histocompatibility complex (MHC) is a gene complex closely linked to the vertebrate immune system due to its importance in antigen recognition and immune response to pathogens. To improve our understanding of the MHC and disease resistance in dairy cattle, we gathered 5119 test day records of somatic cell count (SCC) and performance traits of 262 Holstein dairy cows to determine whether the DRB region of the MHC contains alleles that are associated with elevated SCC, milk yield, protein and fat percent of milk. To this purpose, genotyping of animals for DRB3 gene was investigated by polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) assay. A two-step PCR was carried out so as to amplify a 284 base-pair fragment of exon 2 of the target gene. Second PCR products were treated with three restriction endonuclease enzymes RsaI, BstYI and HaeIII. Twenty-eight BoLA-DRB3 alleles were identified including one novel allele (*40). The results in general are in good accordance with allele frequencies of Holstein cattle populations reported by previous studies. Analyses of associations were modeled based on repeated measurement anova and generalized logistic linear methods for production traits and SCC data, respectively. The results of this study showed a significant relationship between the elevated SCC reflecting an increased probability of occurrence to subclinical mastitis and DRB3.2 allele *8 (p < 0.03). The results also revealed significant positive relationships of alleles*22 (p < 0.01) and allele*11 (p < 0.05) with milk fat percent as well as of alleles*24 (p < 0.03) and *22 (p < 0.05) with protein percent. The present study failed to find any association between milk yield and tested alleles. Because of the lack of consistency among results of similar studies, we suggest further investigations to determine the precise nature of these associations with the high polymorphic bovine MHC region to be performed based on haplotypes.

Characterization of genetic polymorphism of the bovine lymphocyte antigen DRB3.2 locus in Kankrej cattle (Bos indicus)

Bovine lymphocyte antigen DRB 3.2 (BoLA-DRB3.2) gene encodes for the beta chain of the major histocompatibility complex (MHC) class II molecule in cattle, which is a glycoprotein present on the surface of antigen-presenting cells. This locus shows extensive polymorphism in it. The objective of the present study was to genotype the BoLA-DRB3.2 locus in Kankrej cattle (n = 50) by PCR-RFLP. Bovine DNA was isolated from aliquots of whole blood. Primers specific for exon 2 of the bovine lymphocyte antigen (BoLA)-DRB3 gene were used to amplify the region. The 304-bp amplified product of the DRB3 gene was separately digested with restriction endonucleases RsaI, BstYI, and Hae III. Twenty-four BoLA-DRB 3.2 alleles were identified with frequencies ranging from 1 to 22.0%. Twenty-one alleles of the total 24 alleles were similar to those reported earlier; 3 alleles were new and had not been reported previously. The allele BoLA-DRB3.2*34 occurred at the highest frequency of 22% (approx.) in the Kankrej animals studied. Six alleles (BoLA-DRB3.2 *34, *15, *06, *20, *37, and *20) accounted for almost 71% of the total alleles observed to be present in the Kankrej animals. All the new alleles observed were present at frequencies of 1%. The results obtained in the present study demonstrated that the BoLA DRB3.2 locus is highly polymorphic in the Kankrej cattle.

Establishment of a sequence-based typing system for BoLA-DQA1 exon 2

In cattle, bovine leukocyte antigens (BoLAs) have been extensively used as markers for bovine diseases and immunological traits. Here, we developed a rapid, high-resolution sequence-based typing (SBT) system for BoLA-DQA1. We amplified 355 bp of BoLA-DQA1 by fully nested polymerase chain reaction (PCR) using the newly constructed primers and then performed direct sequencing of each product. Using this method, we investigated the locus in 51 animals whose BoLA haplotypes had been characterized at the Fifth International BoLA Workshop. We identified 15 distinct DQA1 alleles, and there is no conflict between the typing result of PCR-SBT and restriction fragment length polymorphism analysis. Together with the previously developed method for typing BoLA-DRB3, the PCR-SBT for BoLA-DQA1 clearly provides a useful tool for detailed class IIa haplotype analysis.

PCR based RFLP genotyping of bovine lymphocyte antigen DRB3.2 in Iranian Holstein population

Major Histocompatibility Complex (MHC) class II locus DRB3 was investigated by PCR based restriction fragment length polymorphism (PCR-RFLP) assay. A total of 262 Holstein cows participating in the national recording system were sampled from 10 herds. A two-step polymerase chain reaction was carried out in order to amplify a 284 base-pair fragment of exon 2 of the target gene. Second PCR products were treated with three restriction endonucleas enzymes RsaI, BstYI and HaeIII. Digested fragments were analyzed by polyacrylamid gel electrophoresis. Twenty-eight BoLA-DRB3 alleles were identified. Identified alleles are: BoLA-DRB3.2 *3, *6, *7, *8, *9, *10, *11, *12, *13, *14, *15, *16,20, *21, *22, *23, *24, *25, *26, *27, *28, *32, *36, *37, *40, *51, *iaa and *ibb. The BoLA-DRB3.2*40 allele that was observed in this study has not been reported previously. The calculated frequencies were as follows: 2.29, 1.34, 0.19, 14.5, 0.38, 3.05, 12.21, 1.34, 2.29, 1.34, 2.48, 9.16, 0.95, 0.77, 6.68, 9.16, 17.94, 1.15, 0.57, 1.15, 0.95, 0.57, 0.38, 1.91, 0.38, 5.73, 0.19 and 0.95% respectively. The six most frequently observed alleles (BoLA-DRB3.2 *8, *11, *16, *22, *23 and *24) accounted for 69.65% of the alleles in these 10 herds. The results of this study confirm the allelic distribution of six most frequent alleles in Holstein population's worldwide.

The diversity of bovine MHC class II DRB3 genes in Japanese Black, Japanese Shorthorn, Jersey and Holstein cattle in Japan

We sequenced exon 2 of the major histocompatibility complex (MHC) class II DRB3 gene from 471 individuals in four different Japanese populations of cattle (201 Japanese Black, 101 Holstein, 100 Japanese Shorthorn, and 69 Jersey cattle) using a new method for sequence-based typing (SBT). We identified the 34 previously reported alleles and four novel alleles. These alleles were 80.0-100.0% identical at the nucleotide level and 77.9-100.0% identical at the amino acid level to the bovine MHC (BoLA)-DRB3 cDNA clone NR1. Among the 38 alleles, eight alleles were found in only one breed in this study. However, these alleles did not form specific clusters on a phylogenetic tree of 236-base pairs (bp) nucleotide sequences. Furthermore, these breeds exhibited similar variations with respect to average frequencies of nucleotides and amino acids, as well as synonymous and non-synonymous substitutions, in all pairwise comparisons of the alleles found in this study. By contrast, analysis of the frequencies of the various BoLA-DRB3 alleles in each breed indicated that DRB3*1101 was the most frequent allele in Holstein cattle (16.8%), DRB3*4501 was the most frequent allele in Jersey cattle (18.1%), DRB3*1201 was the most frequent allele in Japanese Shorthorn cattle (16.0%) and DRB3*1001 was the most frequent allele in Japanese Black cattle (17.4%), indicating that the frequencies of alleles were differed in each breed. In addition, a population tree based on the frequency of BoLA-DRB3 alleles in each breed suggested that Holstein and Japanese Black cattle were the most closely related, and that Jersey cattle were more different from both these breeds than Japanese Shorthorns.

Technical note: DNA typing for ovine MHC DRB1 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP)

The ovine major histocompatibilty complex (Ovar) class II DRB1 second exon was amplified by polymerase chain reaction (PCR) from DNA samples of 52 Suffolk sheep. Polymerase chain reaction products were characterized by the restriction fragment length polymorphism (RFLP) technique using nine restriction enzymes, RsaI, HaeIII, SacI, SacII, DdeI, NciI, Hin1I, EcoRI, and BstNI, yielding 13 types. Sequencing of cloned PCR products identified 16 Ovar-DRB1 alleles. Collectively, all PCR-RFLP patterns exactly matched those predicted from DNA sequences. These findings strongly indicate that the PCR-RFLP method using a combination of nine restriction endonucleases is a very powerful tool in Ovar typing.

Sequences and diversity of 17 new Ovar-DRB1 alleles from three breeds of sheep

To investigate the genetic diversity of the sheep MHC (Ovar) class II DRB1 locus, we amplified exon 2 of Ovar-DRB1 alleles by polymerase chain reaction (PCR) and determined the nucleotide sequences of both resultant strands after cloning. In our study of a total of 97 sheep of three breeds, namely, Suffolk, Cheviot and Corriedale, we identified 18 previously published alleles and 17 new alleles. These alleles were 83.4 to 94.1% identical at the nucleotide level and 71.4 to 90.9% identical at the amino acid level to Ovar-DRB1*0101. We identified six new alleles in Cheviot sheep and 11 new alleles in Suffolk sheep. Furthermore, we identified 15, 6 and 1 allele in Suffolk, Cheviot and Corriedale sheep, respectively, that have only been found in these breeds to date. Analysis of the frequencies of the various Ovar-DRB1 alleles in each breed indicated that Ovar-DRB1*0702 was the most frequent allele in Suffolk sheep (23.9%), Ovar-DRB1*0203 was the most frequent allele in Cheviot sheep (27.5%) and Ovar-DRB1*0201 was the most frequent allele in Corriedale sheep (25.0%). A comparative analysis of the positions of polymorphic residues in the first extracellular domain of the DRB genes of sheep, humans and mice revealed an extraordinary similarity amongst the positions of polymorphic residues that are associated with the antigen recognition site (ARS). Moreover, the extent of polymorphism seems to be similar in sheep, humans and mice.

The influence of ovine MHC class II DRB1 alleles on immune response in bovine leukemia virus infection

We have reported previously that the alleles of the ovine leukocyte antigen (OLA)-DRB1 gene that encode the Arg-Lys (RK) motif and the Ser-Arg (SR) motif at positions beta70/71 of the OLA-DRbeta1 domain are associated with resistance and susceptibility, respectively, to development of bovine leukemia virus (BLV)-induced ovine lymphoma. Here, to investigate the different immune response in sheep that carried alleles associated with resistance and susceptible for 30 weeks after infection with BLV, we selected sheep that had the RK/RK or SR/SR genotype among the 52 sheep analyzed by polymerase chain reaction-restriction fragment length polymorphism and DNA sequencing of PCR product for the OLA-DRB1 exon 2 and infected them with BLV. Although the number of BLV-infected cells and virus titer had been maintaining low levels throughout the experimental period, the sheep with the RK/RK genotype could induce expansion of CD5- B-cells and rapid production of neutralizing antibody in the early phase of infection. The level of incorporation of [3H]thymidine by peripheral blood mononuclear cells from the sheep with RK/RK genotype gave a strong response to BLV virion antigen and synthetic antigenic peptides that corresponded to T-helper epitope of the BLV envelope glycoprotein gp51. In contrast, the sheep with SR/SR genotype showed a strong response to BLV virion antigen and synthetic antigenic peptides that corresponded to T-cytotoxic and B-cell epitopes. In such cases, the animals with the RK/RK strongly expressed IFN-gamma, the animals with SR/SR genotype strongly expressed IL-2. To determine the proliferating cells, we tried a blocking assay with monoclonal antibodies such as anti-CD4, -CD8 and -DR molecule. We found that these proliferating cells were MHC-restricted CD4+ T-cells.