Single nucleotide polymorphisms in the bovine MHC region of Japanese Black cattle are associated with bovine leukemia virus proviral load
© The Author(s) 2017
Received: 10 January 2017
Accepted: 21 March 2017
Published: 4 April 2017
Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leukosis, a malignant B cell lymphoma that has spread worldwide and causes serious problems for the cattle industry. The BLV proviral load, which represents the BLV genome integrated into host genome, is a useful index for estimating disease progression and transmission risk. Here, we conducted a genome-wide association study to identify single nucleotide polymorphisms (SNPs) associated with BLV proviral load in Japanese Black cattle. The study examined 93 cattle with a high proviral load and 266 with a low proviral load. Three SNPs showed a significant association with proviral load. One SNP was detected in the CNTN3 gene on chromosome 22, and two (which were not in linkage disequilibrium) were detected in the bovine major histocompatibility complex region on chromosome 23. These results suggest that polymorphisms in the major histocompatibility complex region affect proviral load. This is the first report to detect SNPs associated with BLV proviral load in Japanese Black cattle using whole genome association study, and understanding host factors may provide important clues for controlling the spread of BLV in Japanese Black cattle.
KeywordsBovine leukemia virus Whole genome association study Major histocompatibility complex
Bovine leukemia virus (BLV), which infects cattle worldwide [1–6], belongs to the family Retroviridae (genus Deltaretrovirus), together with human T cell leukemia virus types 1 and 2 (HTLV-1 and -2) . Historically, the economic losses caused by BLV infection were thought to be related only to the onset of bovine leukosis, which occurs in only 1–5% of BLV-infected cows within 5 years post-infection . However, recent reports show that BLV infection also reduces milk production [6, 8–11] and causes a high incidence of infectious disease  and reproductive inefficiency, resulting in high culling rates ; thus BLV eradication is of utmost importance.
Previous studies show that the proviral load is an important index for estimating the stage of BLV infection because it is associated with disease progression [14–16], lymphocyte count , viral biokinetics , and virus shedding into saliva and nasal secretions . Indeed, one study shows that cattle with a low proviral load are not a source of BLV transmission . Therefore, determining host factors associated with an increased proviral load is important if we are to develop eradication programs for BLV.
Studies of BLV-associated host factors identified polymorphisms within the bovine major histocompatibility complex (MHC) (BoLA) [21–29]. Recently, Miyasaka et al. revealed that polymorphisms within BoLA class II haplotypes were strongly associated with BLV proviral load in Japanese Black cattle, the main breed of beef cattle in Japan, but less so with that in European breeds . However, no group has undertaken a genome-wide association study (GWAS) to identify such host factors.
SNPs showing a significant association with BLV proviral load
Reference cluster IDb
Minor alleled (case)
Minor alleled (control)
1.91 × 10−7
1.91 × 10−7
5.37 × 10−7
Genes within or near these regions were then analyzed using the UMD3.1 genome assembly tool. Hapmap57616-rs29026690 (27,421,348 bp on BTA23) was located between ENSBTAG00000000580 and ABHD16A (ENSBTAG00000000578) (Additional file 1: Table S1), whereas ARS-BFGL-NGS-113235 (28,223,274 bp on BTA23) was located between the 4th and 5th exons of PRR3 (ENSBTAG00000006914). These SNPs reside within the BoLA class III and class I regions, respectively (Figs. 2b, 3; Additional file 1: Table S1). Therefore, the gene density was much higher than that in other areas of the genome, and a number of candidate genes that could be used to estimate proviral load were present around the detected SNPs . Hapmap33580-BTA-136506 was located on the centromeric side of BTA22, at a distance of 6.5 kb from the CONTACTIN3 (CNTN3) gene (Table 1; Additional file 1: Table S1, Additional file 2: Fig. S1).
To the best of our knowledge, this is the first report to detect SNPs associated with BLV proviral load in Japanese Black cattle using GWAS. Two of the identified SNPs were located in the BoLA region. We found it interesting that these two SNPs were located within the class III and class I regions because a previous study reported involvement of only class II genes . The genome reference sequences for the BoLA region have many gaps, mainly because class I genes were difficult to genotype, making associations with class I genes difficult to determine. Target resequencing of high density SNPs across the MHC region using a next generation sequencer should be undertaken to confirm which genes are truly responsible for regulating the proviral load. Our result showed that the MHC polymorphism is important factor for proviral load. The reason why MHC polymorphisms were associated with proviral load is the polymorphism of classical MHC directly associate with antigen presentation and the difference of antigen presentation in each allele leads to the immunological difference in each host.
Taken together, the results described herein show that MHC genotyping of class III and class I alleles can identify cows with a low proviral load. In the farm with high infection rate, eliminating high proviral load cow is an effective way for eradicating BLV because proviral load is major risk factor for transmitting BLV to other host . Therefore, farmer should frequently check the proviral load because the proviral load is variable, although it is not cost-effective. Taken together with the information of our finding 3 SNPs and our previously report about resistant BoLA class II allele , we can identify the BLV resistant cow. It will be helpful to develop a low cost method of eradicating BLV from farms because we can reduce the frequently measurement of proviral load.
Study conception and design: YA. Data acquisition, analysis, and interpretation: SS, MP and ST. Contribution of reagents/materials/analysis tools: YA and YS. Drafting and revising the manuscript: ST, YA, and SS. Approval of the final version: ST, SS, MP, YS, and YA. All authors read and approved the final manuscript.
We thank the Support Unit at the Biomaterial Analysis, RIKEN BSI Research Resources Center, for help with sequence analysis.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets analyzed herein are available from the corresponding author upon request.
Ethics approval and consent to participate
All protocols for the collection of blood samples were reviewed and approved by the Shirakawa Institute of Animal Genetics Committee on Animal Research (H21-2).
The study was supported by Grants-in-Aid for Scientific Research (A and C) from the Japan Society for the Promotion of Science (JSPS) (Grant Nos. 18255013, 25450405), by a grant from the Integration Research for Agriculture and Interdisciplinary Fields in Japan (Grant No. 14538311). This research was also supported by grants from the Project of the NARO Bio-oriented Technology Research Advancement Institution (Grant No. 16817983) (the special scheme project on regional developing strategy).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Erskine RJ, Bartlett PC, Byrem TM, Render CL, Febvay C, Houseman JT. Using a herd profile to determine age-specific prevalence of bovine leukemia virus in Michigan dairy herds. Vet Med Int. 2012;2012:350374.View ArticlePubMedPubMed CentralGoogle Scholar
- Polat M, Ohno A, Takeshima SN, Kim J, Kikuya M, Matsumoto Y, Mingala CN, Onuma M, Aida Y. Detection and molecular characterization of bovine leukemia virus in Philippine cattle. Arch Virol. 2015;160:285–96.View ArticlePubMedGoogle Scholar
- Polat M, Takeshima SN, Hosomichi K, Kim J, Miyasaka T, Yamada K, Arainga M, Murakami T, Matsumoto Y, de la Barra Diaz V, et al. A new genotype of bovine leukemia virus in South America identified by NGS-based whole genome sequencing and molecular evolutionary genetic analysis. Retrovirology. 2016;13:4.View ArticlePubMedPubMed CentralGoogle Scholar
- Polat M, Moe HH, Shimogiri T, Moe KK, Takeshima SN, Aida Y. The molecular epidemiological study of bovine leukemia virus infection in Myanmar cattle. Arch Virol. 2017;162:425–37.Google Scholar
- Zhang R, Jiang J, Sun W, Zhang J, Huang K, Gu X, Yang Y, Xu X, Shi Y, Wang C. Lack of association between bovine leukemia virus and breast cancer in Chinese patients. Breast Cancer Res. 2016;18:101.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Y, Fan W, Mao Y, Yang Z, Lu G, Zhang R, Zhang H, Szeto C, Wang C. Bovine leukemia virus infection in cattle of China: association with reduced milk production and increased somatic cell score. J Dairy Sci. 2016;99:3688–97.View ArticlePubMedGoogle Scholar
- Aida Y, Murakami H, Takahashi M, Takeshima SN. Mechanisms of pathogenesis induced by bovine leukemia virus as a model for human T-cell leukemia virus. Front Microbiol. 2013;4:328.View ArticlePubMedPubMed CentralGoogle Scholar
- Emanuelson U, Scherling K, Pettersson H. Relationships between herd bovine leukemia-virus infection status and reproduction, disease incidence, and productivity in Swedish dairy herds. Prev Vet Med. 1992;12:121–31.View ArticleGoogle Scholar
- Da Y, Shanks RD, Stewart JA, Lewin HA. Milk and fat yields decline in bovine leukemia virus-infected Holstein cattle with persistent lymphocytosis. Proc Natl Acad Sci USA. 1993;90:6538–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Norby B, Bartlett PC, Byrem TM, Erskine RJ. Effect of infection with bovine leukemia virus on milk production in Michigan dairy cows. J Dairy Sci. 2016;99:2043–52.View ArticlePubMedGoogle Scholar
- Nekouei O, VanLeeuwen J, Stryhn H, Kelton D, Keefe G. Lifetime effects of infection with bovine leukemia virus on longevity and milk production of dairy cows. Prev Vet Med. 2016;133:1–9.View ArticlePubMedGoogle Scholar
- Sandev N, Koleva M, Binev R, Ilieva D. Influence of enzootic bovine leukosis virus upon the incidence of subclinical mastitis in cows at a different stage of infection. Vet Arh. 2004;74:411–6.Google Scholar
- Bartlett PC, Norby B, Byrem TM, Parmelee A, Ledergerber JT, Erskine RJ. Bovine leukemia virus and cow longevity in Michigan dairy herds. J Dairy Sci. 2013;96:1591–7.View ArticlePubMedGoogle Scholar
- Jimba M, Takeshima SN, Matoba K, Endoh D, Aida Y. BLV-CoCoMo-qPCR: quantitation of bovine leukemia virus proviral load using the CoCoMo algorithm. Retrovirology. 2010;7:91.View ArticlePubMedPubMed CentralGoogle Scholar
- Jimba M, Takeshima SN, Murakami H, Kohara J, Kobayashi N, Matsuhashi T, Ohmori T, Nunoya T, Aida Y. BLV-CoCoMo-qPCR: a useful tool for evaluating bovine leukemia virus infection status. BMC Vet Res. 2012;8:167.View ArticlePubMedPubMed CentralGoogle Scholar
- Somura Y, Sugiyama E, Fujikawa H, Murakami K. Comparison of the copy numbers of bovine leukemia virus in the lymph nodes of cattle with enzootic bovine leukosis and cattle with latent infection. Arch Virol. 2014;159:2693–7.View ArticlePubMedGoogle Scholar
- Ohno A, Takeshima SN, Matsumoto Y, Aida Y. Risk factors associated with increased bovine leukemia virus proviral load in infected cattle in Japan from 2012 to 2014. Virus Res. 2015;210:283–90.View ArticlePubMedGoogle Scholar
- Panei CJ, Takeshima SN, Omori T, Nunoya T, Davis WC, Ishizaki H, Matoba K, Aida Y. Estimation of bovine leukemia virus (BLV) proviral load harbored by lymphocyte subpopulations in BLV-infected cattle at the subclinical stage of enzootic bovine leucosis using BLV-CoCoMo-qPCR. BMC Vet Res. 2013;9:95.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuan Y, Kitamura-Muramatsu Y, Saito S, Ishizaki H, Nakano M, Haga S, Matoba K, Ohno A, Murakami H, Takeshima SN, Aida Y. Detection of the BLV provirus from nasal secretion and saliva samples using BLV-CoCoMo-qPCR-2: comparison with blood samples from the same cattle. Virus Res. 2015;210:248–54.View ArticlePubMedGoogle Scholar
- Juliarena MA, Barrios CN, Ceriani MC, Esteban EN. Hot topic: bovine leukemia virus (BLV)-infected cows with low proviral load are not a source of infection for BLV-free cattle. J Dairy Sci. 2016;99:4586–9.View ArticlePubMedGoogle Scholar
- Nikbakht Brujeni G, Ghorbanpour R, Esmailnejad A. Association of BoLA-DRB3.2 alleles with BLV infection profiles (persistent lymphocytosis/lymphosarcoma) and lymphocyte subsets in Iranian Holstein cattle. Biochem Genet. 2016;54:194–207.View ArticlePubMedGoogle Scholar
- Miyasaka T, Takeshima SN, Jimba M, Matsumoto Y, Kobayashi N, Matsuhashi T, Sentsui H, Aida Y. Identification of bovine leukocyte antigen class II haplotypes associated with variations in bovine leukemia virus proviral load in Japanese Black cattle. Tissue Antigens. 2013;81:72–82.View ArticlePubMedGoogle Scholar
- Juliarena MA, Poli M, Sala L, Ceriani C, Gutierrez S, Dolcini G, Rodriguez EM, Marino B, Rodriguez-Dubra C, Esteban EN. Association of BLV infection profiles with alleles of the BoLA-DRB3.2 gene. Anim Genet. 2008;39:432–8.View ArticlePubMedGoogle Scholar
- Aida Y, Nishino Y, Amanuma H, Murakami K, Okada K, Ikawa Y. The role of tumor-associated antigen in bovine leukemia virus-induced lymphosarcoma. Leukemia. 1997;11(Suppl 3):216–8.PubMedGoogle Scholar
- Zanotti M, Poli G, Ponti W, Polli M, Rocchi M, Bolzani E, Longeri M, Russo S, Lewin HA, van Eijk MJ. Association of BoLA class II haplotypes with subclinical progression of bovine leukaemia virus infection in Holstein–Friesian cattle. Anim Genet. 1996;27:337–41.PubMedGoogle Scholar
- Xu A, van Eijk MJ, Park C, Lewin HA. Polymorphism in BoLA-DRB3 exon 2 correlates with resistance to persistent lymphocytosis caused by bovine leukemia virus. J Immunol. 1993;151:6977–85.PubMedGoogle Scholar
- Stear MJ, Dimmock CK, Newman MJ, Nicholas FW. BoLA antigens are associated with increased frequency of persistent lymphocytosis in bovine leukaemia virus infected cattle and with increased incidence of antibodies to bovine leukaemia virus. Anim Genet. 1988;19:151–8.View ArticlePubMedGoogle Scholar
- Lewin HA, Wu MC, Stewart JA, Nolan TJ. Association between BoLA and subclinical bovine leukemia virus infection in a herd of Holstein–Friesian cows. Immunogenetics. 1988;27:338–44.View ArticlePubMedGoogle Scholar
- Lewin HA, Bernoco D. Evidence for BoLA-linked resistance and susceptibility to subclinical progression of bovine leukaemia virus infection. Anim Genet. 1986;17:197–207.View ArticlePubMedGoogle Scholar
- Sasaki S, Ibi T, Watanabe T, Matsuhashi T, Ikeda S, Sugimoto Y. Variants in the 3′ UTR of general transcription factor IIF, polypeptide 2 affect female calving efficiency in Japanese Black cattle. BMC Genet. 2013;14:41.View ArticlePubMedPubMed CentralGoogle Scholar
- Takeshima SN, Kitamura-Muramatsu Y, Yuan Y, Polat M, Saito S, Aida Y. BLV-CoCoMo-qPCR-2: improvements to the BLV-CoCoMo-qPCR assay for bovine leukemia virus by reducing primer degeneracy and constructing an optimal standard curve. Arch Virol. 2015;160:1325–32.View ArticlePubMedGoogle Scholar
- Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012;44:821-U136.Google Scholar
- Consortium TMs. Complete sequence and gene map of a human major histocompatibility complex. The MHC sequencing consortium. Nature. 1999;401:921–3.View ArticleGoogle Scholar