Skip to main content
  • Research article
  • Open access
  • Published:

Genetic variation in Japanese Holstein cattle for EBL development



Infection with bovine leukemia virus (BLV), the causative agent for enzootic bovine leukosis (EBL), is increasing in dairy farms of Japan. The tendency of tumor development following BLV infection in certain cow families and bull lines has previously been described. We therefore hypothesized the existence of a genetic component which differentiates cattle susceptibility to the disease.


We analyzed routinely collected large-scale data including postmortem inspection data, which were combined with pedigree information and epidemiological data of BLV infection. A total of 6,022 postmortem inspection records of Holstein cattle, raised on 226 farms served by a regional abattoir over 10 years from 2004 to 2015, were analyzed for associations between sire information and EBL development. We then identified statistically the relative susceptibility to EBL development for the progeny of specific sires and paternal grandsires (PGSs). The heritability of EBL development was calculated as 0.19. Similarly, proviral loads (PVLs) of progeny from identified sires and PGSs were analyzed, but no significant differences were found.


These observations suggest that because EBL development in our Holstein population is, at least in part, influenced by genetic factors independent of PVL levels, genetic improvement for lower incidence of EBL development in cattle notwithstanding BLV infection is possible.


Among the many environmental and managemental factors affecting dairy cattle production, enzootic bovine leukosis (EBL), caused by the bovine leukemia virus (BLV) [1], an oncogenic retrovirus, undoubtedly has a negative impact on dairy productivity. It is now believed that 79.1% of dairy farms and 34.7% of dairy cattle in Japan are infected with BLV [2, 3]. BLV is transmitted mainly by transferring infected blood through contaminated needles, rectal palpation sleeves, dehorners or bloodsucking insects. Because infected cattle are lifelong carriers and remain a reservoir of infection for other cattle, reducing herd BLV prevalence through segregation or gradual culling of BLV-positive animals is a lengthy process. The high prevalence of BLV may affect cattle farms by decreasing milk yield and necessitating carcass condemnation at slaughterhouses [4]. The number of reported cases of Bovine Leukosis (BL), which includes EBL and Sporadic Bovine Leukosis (SBL), has been increasing every year since BL was designated as a notifiable disease in Japan in 1998.

The possibility of genetic contributions to susceptibility to BL has been reported in previous studies [5,6,7]. Although about 1 to 5% of BLV infected cattle develop lymphosarcoma at ages 4 to 5, more than 90% of them remain asymptomatic throughout their productive life [7]. A recent study indicated that although increasing age and high proviral load (PVL) increase risks of EBL development, not all aged cattle with higher PVL develop neoplastic tumors. Rather, only 14% of the high risk group with age older than 4.5 years and PVL higher than 52,000 copies/105 cells were diagnosed with EBL [8]. Genetic components that promote tumor development from both paternal and maternal pedigrees have been described in several previous reports published in the middle of the 20th century, before the identification of BLV as the causative agent of EBL. These previous reports have shown evidence of a family of cattle highly susceptible to lymphoma development that all descended from one bull in East Prussia [9], a significant difference in predisposition between two bulls in the same herd [10], and a significant positive correlation between daughters and dams with regard to the development of leukotic tumors [11]. These reports indicate that the susceptibility to the development of neoplastic tumors, now known to be EBL, could be heritable in cattle populations. However, as these data are based on only a few herds with a small data set, it is difficult to statistically evaluate the genetic variability associated with EBL development among pedigrees. Moreover, such a small data set makes it difficult to eliminate the influence of environmental factors on EBL development.

Since the year 2003, the Japanese dairy industry has participated in the International Bull Evaluation Service (INTERBULL) program, and male stock information has been accumulated in the Domestic Animal Improvement (DAI) databank of Japan. In 2003, a national program under the Cattle Traceability Law was implemented to trace each animal from birth to slaughter. All cattle were mandatorily registered with an individual identification (ID) number in the database established at the National Livestock Breeding Center (Fukushima, Japan). Each cattle ID entry retains extensive information, including dam ID, breed, sex and all translocation history from birth to slaughter, and is linked with the DAI databank. In 1998, bovine leukosis was added to the list of reportable animal diseases in the Domestic Animal Infectious Diseases Control Law. Subsequently, cattle diagnosed as EBL were subjected to compulsory destruction, and diagnostic information at the time of postmortem inspection was stored by the veterinary officers of Municipal Meat Inspection Stations.

The objective of this study was to use these data to assess the information supporting a genetic variation regarding EBL development as well as BLV infection in Holstein cattle, and to gain insight into the mode of inheritance. These objectives were examined through a case-control study of the Holstein cattle registered in the DAI databank for the period of 2004 to 2015.


Study 1-Analysis of postmortem inspection data

Postmortem inspection records of 6,022 Holstein cattle raised on 226 farms were analyzed for associations between pedigree and EBL development (Fig. 1). Among this study population, 69 cattle were diagnosed as EBL at the local abattoir, while the remaining cattle exhibited no symptoms of development of lymphoid tumors (Table 1). In this study area, none of the cattle were reported as EBL during the study period and thus all EBL diagnosis was made at the abattoir. From cattle IDs, pedigree information including sire, paternal grand sire (PGS) and maternal grand sire (MGS) of each animal were analyzed. The 6,022 Holstein cattle analyzed in Study1 were sired by 1,070 bulls (defined in this study as Sire). These Sires were bred by 288 bulls (defined as PGS). Regarding maternal lineage, although the lineage information was not available for dams of 51 cattle and grandsires for 291 cattle among 6,022 cattle, 5,971 cattle were born from 4,921 dams. These dams were sired by more than 958 sires (defined as MGS). For the EBL case population, 52 sires, 39 PGSs and more than 53 MGSs were identified. All sires with at least 50 daughters or granddaughters were included in the case/control file, from which tables were built and Fisher’s exact test was performed to compare cell contents and compute the odds ratio (OR) for each sire (Table 1). In this analysis, each sire/grandsire was compared against all other sires and grandsires. ORs are presented as 95% confidence intervals. Among all sires of the study population, 13 sires had daughters numbering at least 50 heads. Sire9 and Sire13 present significantly high frequencies of EBL cases (OR = 3.99 (1.04–12.24), p = 0.046 and OR = 6.04 (1.95–16.97), p = 0.006). PGS1 and PGS3 present significantly low frequencies of EBL cases (p = 0.033 and p = 0.049), while PGS6 and PGS12 present significantly high frequencies of EBL cases (OR = 2.52 (1.05–5.95), p = 0.042 and OR = 3.07 (1.00–8.24), p = 0.049) in the granddaughter population. None of the MGS exhibited a significant difference in the development of EBL in their progenies (Table 1). There are no significant differences in longevity between cattle that developed lymphoma and cattle that did not develop lymphoma (5.64 years vs. 5.96 years, p = 0.154, Mann-Whitney U Test).

Table 1 No. of daughters/granddaughters slaughtered and examined for EBL at postmortem inspection listed by sire/grandsire
Fig. 1
figure 1

Schematic overview of the analysis on cattle data

The data from four data sources were combined according to cattle ID, and the cattle data of which ID was registered to DAI data bank was further extracted. KSLHS: Kanagawa Shonan Livestock Hygiene Station, MISK: Meat Inspection Station, Kanagawa Prefectural Government, DAI: Domestic Animal Improvement, PGS: Paternal grandsire, MGS: Maternal grand sire

Study 2-Heritability estimate of BLV infection and EBL susceptibility

Based on the data of Study 1, we identified the sires and PGSs that exhibited significant differences in the development of EBL in their daughter and granddaughter populations. These data indicated that some lineages either more susceptible or less susceptible to EBL tumor development could exist. To estimate the heritability of BLV infection and EBL susceptibility among Holstein cows, the dataset from IDs was linked to the lineage information available at the DAI databank. In this analysis, the heritability was calculated by using the BLUP method; the diagnostic information regarding EBL and the presence or absence of BLV antibodies were treated as random effects, while age at analysis, the year of analysis, on-farm BLV infection and EBL development were treated as fixed effects. The heritability of EBL development was 0.19, where the additive genetic variance (\({\sigma }_{g}^{2})\) was 0.0508 and the residual variance (\({\sigma }_{r}^{2})\) was 0.222.

Study 3- PVL of progenies of detected sires and grandsires at farm level

Previously, it was found that cattle exhibiting disease-related tumors following BLV infection had higher PVLs [8, 12]. The question then arises whether the difference in EBL frequencies among pedigrees is associated with PVL at the individual animal level. Therefore, we analyzed the data of PVL levels of progenies of identified sires and PGSs. The PVL of these cattle ranged from zero to 374,025 copies/105 cells, the median was 42,072 copies/105 cells (Table 2). There was no significant difference in PVL levels between the progenies of the identified sires and PGSs and the progenies of other sires and grandsires.

Table 2 Seroprevalence and PVL per sires with significant odds ratio


Most BLV-infected animals remain asymptomatic, although a small fraction of infected animals exhibit tumor development during the course of disease progression. The tendency for tumor development after BLV infection in certain cow families and bull lines has previously been described [9,10,11]. However, there has been no statistical evidence to support the existence of genetic components that influence the degree of susceptibility to the disease. In this study, we analyzed routinely collected large-scale data including postmortem inspection data, which were combined with EBL diagnostic data and pedigree information. We then identified statistically the sires and PGSs that were more susceptible or less susceptible to the development of EBL (Table 1). These findings were made possible because the information related to BLV epidemiological data at farms and EBL diagnostic data at the local abattoir were all linked through each cattle ID. Based on this information, it was possible to estimate the heritability of susceptibility to the BLV infection and EBL development as 0.08 and 0.19, respectively. The former estimate agrees with the heritability calculated from BLV prevalence in milk [5]. These observations suggest that development of EBL in our Holstein population is, at least in part, influenced by genetic factors.

Previous studies have shown that the specific alleles of the major histocompatibility complex, bovine leukocyte antigen (BoLA) genes, are significantly associated with PVL in BLV-infected animals. BoLA-DRB3*009:02, BoLA-DRB3*002:01 and DRB3*014:01:01 are reported as the alleles associated with low PVL in several studies [13,14,15]. In the population examined in Study 3, we analyzed the BoLA DRB3 haplotype of granddaughters of identified PGSs. However, paternal granddaughters with these particular alleles were found in only 23.3% (BoLA-DRB3*009:02: 6.7%, BoLA-DRB3*002:01: 0% and DRB3*014:01:01: 20%) of those examined (data not shown). In addition, as 99% of Holstein cows were asymptomatic, these results indicate that genetic factors that affect disease development may be not only BoLA but also others yet to be identified.

A limitation of our study is that in the analysis of Study 1, we assessed the frequency of EBL development at the time of slaughtering. It is possible that the difference of the frequency of EBL development could be biased by cattle that were slaughtered before the age of tumor development. In the population of Study 1, the median longevities of the progenies of PGS1 which display no EBL development was slightly shorter than that of other PGSs (5.0 years vs. 5.7 years, p < 0.001, Mann-Whitney U test). This may have biased the results and should be noted when interpreting the data. However, the differences in longevity of the progenies of PGS1 are less than a year, and do not fully explain the absence of EBL cattle because the age of occurrence of EBL ranges from 1.6 years to 9.8 years in this study.

In our study, Holstein cattle were from a small region in Japan, but ID numbers and male stock information accumulated in the DAI data bank enabled us to trace definitive information on individual animals, sire information and EBL development. It should be noted that BLV seroprevalence rates in individual cows and farms were similar, if not the same, to those reported for the entire population of Holstein cattle in Japan [3]. The seropositive rates were also similar to those of the US [16, 17] and Canada [4] due at least in part to Japan’s participation in INTERBULL since 2003. These data indicate that Holstein sire usage is similar in these countries, exhibiting genetic homology and homogeneity. It was reported that genetic correlation is increasing in dairy countries, resulting from the use of genetically superior sires [18]. VanRaden et al. [19] also reported that genetic makeups of Holstein sires in major dairy countries including Japan are similar to those of the US. Identification of the sires and PGSs with reduced frequency of EBL development in their progenies suggests that sires potentially identified through the use of methods tested in this study could be applied to reduce economic losses related to EBL in major dairy countries.

It is still an open question as to whether focusing on prevention of EBL development through genetic selection rather than reduction in BLV infection or control of PVL levels once infected is most beneficial to dairy production and/or population. In this study, all EBL cases were found at a local abattoir and therefore most of the dairy cattle infected with BLV had survived through their productive lives without exhibiting clinical symptoms or fatal lymphoma. Our observations indicate that sire selection from the grandsires less susceptible to developing EBL in their progenies could serve as one method in reducing EBL-incidence in dairy producing countries.


In this study, we demonstrated that among the many sires used in Japan, there exist sires and PGSs exhibiting increased or decreased frequency of tumor development in their progenies, although PVL did not differ among the progenies of sires and PGSs examined. The observed level of heritability indicates that genetic contributions are of importance in EBL development in Holstein cattle. Although the sources of genetic variability or patterns of inheritance have not been elucidated, the use of PGSs that are less susceptible to tumor development could reduce the frequency of EBL in the cattle population. It is possible that more such PGSs could be identified if analyses similar to this study are carried out elsewhere in the world.


The investigation was divided into three components, as follows:

Study 1 - Analysis of postmortem inspection data

Study area and data collection

In this study, we used the data from a local abattoir within jurisdiction of the Meat Inspection Station, Kanagawa Prefectural Government (MISK), as collected from 2004 to 2015. This abattoir receives 4,000 to 5,000 cattle annually for slaughtering, mainly from Kanagawa prefecture and nearby Shizuoka prefecture, both of which are located in the central part of Japan near the capital of Tokyo. The diagnostic data of the postmortem inspection of all slaughtered cattle were recorded with each cattle ID by the veterinary officers. A total of 67,189 inspection records were included in the data set before editing. Data were further edited to include only the records of Holstein cattle received at the local abattoir from farms within the jurisdiction of the Kanagawa Shonan Livestock Hygiene Service Center (KSLHS), so that information of the cattle diagnosed with EBL in these farms was also available. Additionally, data were further edited to include only Holstein cattle registered to the DAI databank so that pedigree information could be examined. A total of 6,022 postmortem records were included in the final data set (Fig. 1). These data constituted the case group (all cattle diagnosed with EBL) and the control group (all remaining cattle), and tabulated by pedigree information (sire, PGS, and MGS).

EBL diagnosis

For all animals slaughtered at the local abattoir, when neoplastic lesions are found at inspection, the carcass is separately stored and further laboratory examinations are performed: lymphoma or lymphosarcoma is diagnosed by histologic examination, followed by additional immunohistochemical analysis to determine T-cell or B-cell origin. Subsequently, antibodies against BLV in serum are detected with a passive hemagglutination test (Nisseiken Co., Ltd., Ome, Japan).

Study 2-Heritability estimate of BLV infection and EBL development

Study area

During the course of data collection in Study 1, a serological survey of 4,114 randomly selected cattle was performed at the KSLHS. The results of serological testing were cross-referenced with the final data set of Study 1 to extract cattle ID found in both data sets. These IDs were tabulated with results of the BLV serological test and diagnostic information regarding development of EBL, from which their pedigrees were further evaluated (Fig. 1).

Heritability estimate

A total of 1,940 cattle confirmed through ID to be present in all three relevant databases, namely the serological test, postmortem inspection data and DAI data banks, were selected for this study. Based on this information, the heritability for EBL development was calculated using a threshold model, presented in matrix notation as y = Xβ + Zα + e, where y is a vector of observations, β is a vector with the fixed effects of year of birth (17 levels), age of testing (13 levels) and 98 farms (3 prevalence levels), α is a vector of additive genetic effects (i.e., breeding values), e is a vector of residual effects, and X and Z are prevalence matrices that link individual susceptibility to the BLV infection and development of lymphoid tumors to their respective fixed and random effects.

Using the THRGIBBS1F90 software program [20], Gibbs sampling programs for categorical traits and Bayesian analysis were performed to estimate the variance components in the models. Single chains of 100,000 cycles were run with the first 50,000 cycles used as the burn-in period, followed by post-Gibbs analysis using the POSTGIBBSF90 software program [21]. These repeated analyses enabled the estimation of (co)variance, G-matrix and R-matrix, resulting in the calculation of the heritability for BLV infection and development of EBL. In this analysis, there was no significant difference between sire usage in the study area and that of Japan as a whole (data not shown).

Study 3-PVL of progenies of sires on farm

Study area and sample collection

To analyze whether or not genetic variability is associated with susceptibility to BLV infection and/or PVL levels, the epidemiological data of daughters/granddaughters of identified sires/PGSs were obtained from the survey conducted within the same study area in 2015. Sampling and testing procedures were performed as previously described in detail [8]. In addition to the previous data from only the year 2015, blood samples were further collected using the same methods from a neighboring city in 2016. Serological tests were conducted and PVL levels of seropositive cattle were measured. The results of the epidemiological surveys for the two cities were analyzed together. A total of 2,179 cattle from 58 farms were included in the final data set. Among these data, cattle raised in farms without seropositive animals were excluded. Additionally, data were further edited to include only Holstein cattle registered to the DAI databank. A total of 1,069 cattle data from 49 farms were included in the final data set (Fig. 1). These records were tabulated by the results of serological test, PVL levels for seropositive animals, and distinct pedigrees.

Statistical analysis

Odds ratio (OR) was estimated by Fisher’s exact test. Longevity for each cow was calculated from the date of birth to the date of culling or date of death. Longevities and BLV PVLs (copy numbers/105 cells) were used as a continuous variable. Continuous data was presented as median (range) values and compared using a Mann-Whitney U test. All statistical analyses were performed with SPSS, R version 3.2.1 (The R Foundation for Statistical Computing, Vienna, Austria) and EZR version 1.36 (Saitama Medical Center, Jichi Medical University, Saitama, Japan), which is a graphical user interface for R.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request and with permission of Kanagawa Prefectural Government.



Enzootic Bovine Leukosis


Bovine Leukemia Virus


Bovine Leukosis


Sporadic Bovine Leukosis


Proviral load


Domestic Animal Improvement


Individual identification


Meat Inspection Station, Kanagawa Prefectural Government


Kanagawa Shonan Livestock Hygiene Service Center


Paternal grand sire


Maternal grand sire


Odds ratio


Bovine leukocyte antigen


  1. Gillet N, Florins A, Boxus M, Burteau C, Nigro A, Vandermeers F, Balon H, Bouzar AB, Defoiche J, Burny A, et al. Mechanisms of leukemogenesis induced by bovine leukemia virus: prospects for novel anti-retroviral therapies in human. Retrovirology. 2007;4:18.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Murakami K, Kobayashi S, Konishi M, Kameyama K, Yamamoto T, Tsutsui T. The recent prevalence of bovine leukemia virus (BLV) infection among Japanese cattle. Vet Microbiol. 2011;148(1):84–8.

  3. Murakami K, Kobayashi S, Konishi M, Kameyama K, Tsutsui T. Nationwide survey of bovine leukemia virus infection among dairy and beef breeding cattle in Japan from 2009–2011. The Journal of Veterinary Medical science. 2013;75(8):1123–6.

  4. 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.

    Article  PubMed  Google Scholar 

  5. Abdalla EA, Rosa GJ, Weigel KA, Byrem T. Genetic analysis of leukosis incidence in United States Holstein and Jersey populations. J Dairy Sci. 2013;96(9):6022–9.

    Article  CAS  PubMed  Google Scholar 

  6. Abdalla EA, Weigel KA, Byrem TM, Rosa GJM. Short communication: Genetic correlation of bovine leukosis incidence with somatic cell score and milk yield in a US Holstein population. J Dairy Sci. 2016;99(3):2005–9.

    Article  CAS  PubMed  Google Scholar 

  7. Rodriguez SM, Florins A, Gillet N, de Brogniez A, Sanchez-Alcaraz MT, Boxus M, Boulanger F, Gutierrez G, Trono K, Alvarez I, et al. Preventive and therapeutic strategies for bovine leukemia virus: lessons for HTLV. Viruses. 2011;3(7):1210–48.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Kobayashi T, Inagaki Y, Ohnuki N, Sato R, Murakami S, Imakawa K. Increasing Bovine leukemia virus (BLV) proviral load is a risk factor for progression of Enzootic bovine leucosis: A prospective study in Japan. Prev Vet Med. 2019;S0167-5877(18):30795–5.

  9. Czymoch O. Contribution to the etiology of bovine leukosis in East Prussia. Z Hyg Infektionskrankh. 1937;52(189).

  10. Fortner J. The leukosis of cattle. MhTierheilk. 1953;5.

  11. Henricson B, Olson H. Statistical studies on bovine leukosis. Acta Vet Scand. 1961;2(supple.2):55–62.

  12. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 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(5):338–44.

    Article  CAS  PubMed  Google Scholar 

  14. Nagaoka Y, Kabeya H, Onuma M, Kasai N, Okada K, Aida Y. Ovine MHC class II DRB1 alleles associated with resistance or susceptibility to development of bovine leukemia virus-induced ovine lymphoma. Cancer Res. 1999;59(4):975–81.

    CAS  PubMed  Google Scholar 

  15. Takeshima SN, Matsumoto Y, Aida Y. Short communication: Establishment of a new polymerase chain reaction-sequence-based typing method for genotyping cattle major histocompatibility complex class II DRB3. J Dairy Sci. 2009;92(6):2965–70.

    Article  CAS  PubMed  Google Scholar 

  16. LaDronka RM, Ainsworth S, Wilkins MJ, Norby B, Byrem TM, Bartlett PC. Prevalence of Bovine Leukemia Virus Antibodies in US Dairy Cattle. Vet Med Int. 2018;2018:5831278.

    Article  PubMed  PubMed Central  Google Scholar 

  17. White TL, Moore DA. Reasons for whole carcass condemnations of cattle in the United States and implications for producer education and veterinary intervention. J Am Vet Med Assoc. 2009;235(8):937–41.

    Article  PubMed  Google Scholar 

  18. Battagin M, Forabosco F, Jakobsen JH, Penasa M, Lawlor TJ, Cassandro M. International genetic evaluation of Holstein bulls for overall type traits and body condition score. J Dairy Sci. 2012;95(8):4721–31.

    Article  CAS  PubMed  Google Scholar 

  19. VanRaden PM, Sanders AH, Tooker ME, Miller RH, Norman HD, Kuhn MT, Wiggans GR. Development of a national genetic evaluation for cow fertility. J Dairy Sci. 2004;87(7):2285–92.

    Article  CAS  PubMed  Google Scholar 

  20. Misztal I, Legarra A, Aguilar I. Computing procedures for genetic evaluation including phenotypic, full pedigree, and genomic information. J Dairy Sci. 2009;92(9):4648–55.

    Article  CAS  PubMed  Google Scholar 

  21. Misztal I, Tsuruta S, Strabel T, Auvray B, Druet T, Lee DH: BLUPF90 AND RELATED PROGRAMS (BGF90). 7th World Congress on Genetics Applied to Livestock Production, August 19–23, 2002, Montpellier, France 2002, 28 – 07.

Download references


We would like to express our gratitude to the officers of the Holstein Cattle Association of Japan for their helpful suggestions on the lineage evaluation of sires, the staff of the KSLHS for blood sample collection, the staff of MISK for providing information of EBL cattle, and Ms. Chiaki Mashima for supporting data editing. We also thank Mr. Robert Moriarty for his English editing.


This work was supported by Livestock Promotional Funds of Japan Racing Association (JRA) (to K.I), and by Science Research Promotion Fund of the Japan Private School Promotion Funding, Japan (to T.K).

Author information

Authors and Affiliations



Y.I. AND K.I. conceived and designed the study; Y.I., T.K. AND Y.S. performed the analysis; K.K. contributed in analysis support and results discussion: Y.I., T.K., Y.S. AND K.I. wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kazuhiko Imakawa.

Ethics declarations

Ethics approval and consent to participate

The research was conducted by using data obtained through the information publicly available from Kanagawa Prefectural Government reports. All data or information excludes personally identifiable information such as farmers’ names. The pedigree information of the Holstein cows registered in Kanagawa was obtained from the DAI databank through web-search using their IDs. The dairy farm leaders gave verbal consent for the blood sample collection, serological test of BLV antibodies and PVL measurements. The serum samples and EDTA-treated blood samples used in this study were collected by KSLHS from 2015 to 2016 national whole herd test surveillance targeting Mycobacterium avium subsp. paratuberculosis (MAP) infection and Johne’s disease.

Consent for publication

Not applicable.

Competing interests

Non-financial competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Inagaki, Y., Kobayashi, T., Suda, Y. et al. Genetic variation in Japanese Holstein cattle for EBL development. BMC Vet Res 16, 407 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: