- Open Access
Sequencing and phylogenetic analysis of the gp51 gene from Korean bovine leukemia virus isolates
Virology Journalvolume 12, Article number: 64 (2015)
Bovine Leukemia virus (BLV) infection of cattle has been reported in Korea for more than three decades. However, to date, there have been few studies regarding Korean BLV since 1980s. Thus, the purpose of this study is to perform a diagnosis and molecular characterization of BLV strains circulating in Korea and to estimate genetic diversity of different genotypes of BLV.
To investigate the distribution of BLV variants in the world and assess the evolutionary history of Korean BLV isolates, a comprehensive molecular analysis of the BLV env gp51 gene was conducted using recent worldwide BLV isolates. The isolates included 50 samples obtained from two cattle farms in southeastern Korea in 2014.
Sequence and phylogenetic analyses of partial 444-nt fragment sequences and complete gp51 sequences of BLV revealed eight distinct genotypes of BLV showing geographic distribution of the world. Most Korean BLV isolates were found to belong to genotype 1 which is a major genotype prevailed throughout the world, and only four isolates from one farm were classified as genotype 3 related to the US and Japan isolates. Analysis of amino acids of Korean BLV isolates showed several sequence substitutions in the leader peptide, conformational epitope, and neutralizing domain regions. The observations suggest the possibility of affecting on viral infectivity and formation.
Korean BLV isolates showed the close relationship to genotype 1 and 3. Further study to identify the diversity of BLV circulating in Korea is necessary with samples collected nationwide because this study is the first report of BLV genotype 3 being in circulation in Korea.
Bovine Leukemia virus (BLV) belongs to the genus Deltaretrovirus of the Retroviridae family and is the causative agent of enzootic bovine leukosis (EBL), a disease that results in economic loss for the cattle industry . BLV infected cattle are asymptomatic; 30 ~ 70% of infected animals develop persistent lymphocytosis (PL) and 0.1 ~ 10% develop lymphoid tumors [2-4]. BLV infection can be found worldwide, but many European countries have successfully eradicated EBL in recent years [1,4-14].
BLV has an envelope (Env) protein complex that is composed of surface glycoprotein (SU) subunits, which are anchored to virions by their association with transmembrane (TM) protein subunits . The envelope glycoproteins of BVL play a crucial role in the virus life cycle; envelope proteins are responsible for cellular tropism because they contain the recognition site for the cell surface receptor required for virus entry and are the natural target for neutralizing antibodies .
Recent phylogenetic analysis of BLV env genes showed that this virus could be classified into 8 or more genotypes [17,18]. EBL was first reported by the serological method in Korea in the 1980s , and recently, its prevalence in dairy cattle was found to be over 50% [20-23]. To date, there has been no reports of the BLV genotype circulating in Korea since partial env gene sequences that belonged to genotype 1 were deposited in 2009 . It is known that some of the genotype groupings correlate with the geographical origin of the strain [16,25-34], but with the increased movement of cattle and people around the world, a variety of genotypes have been reported in places where only a specific genotype was dominant in the past [17,33]. In Japanese studies, it was reported that many genotypes of BLV (genotypes 1 ~ 5) co-circulated nationwide  with genotype 1 being predominant and genotypes 3 and 5 being widely distributed .
This study was performed to identify the genetic diversity of BLV in Korean cattle. The gp51 sequence in 50 BLV isolates was analyzed in our study.
A total of 185 whole blood samples were collected from dairy cattle on two farms (119 from farm A and 66 from farm B) in southeast Korea. Farm A is located in Gyeongsangbuk-do and farm B in Gyeongsangnam-do. The samples were collected between April and July 2014.
Viral DNA extraction
Peripheral blood lymphocytes (PBLs) were isolated from whole blood samples and stored at −70°C until further use. PBLs isolation was performed using Luecosep tube (Greiner Bio-One Gmbh, Frickenhausen, Germany) according to the manufacturer’s protocol. Isolated PBLs were resuspended in 1 ml of PBS. Viral DNA was extracted from 200 μL of PBL solution with QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany) following the manufacturer’s recommendations. DNA samples were stored at −20°C until further use.
PCR amplification of the env gene fragment
BLV gp51 gene amplification for diagnosis was carried out using primers sets previously described by Fechner et al. . The PCR reaction was performed using a Maxime PCR PreMix Kit (iNtRON Biotechnongy, Gyeonggi-do, Korea). The first PCR was performed using a 20-μL reaction volume. For one reaction, the assay was optimized to contain 1× PCR buffer, 1 μL of DNA (~1 μg), 1 μL of each of the env-specific primers BLV-env-1 (5′-TCTGTGCCAAGTCTCCCAGATA-3′) and BLV-env-2 (5′-AACAACAACCTCTGGGAAGGG-3′) (20 pmol/μL), 2.5 mM of each dNTP, 2.5 units of i-Taq polymerase. The reaction parameters were 2 minutes at 94°C, followed by 30 cycles of 30 seconds at 95°C, 30 seconds at 58°C and 60 seconds at 72°C and 4 minutes at 72°C. The nested PCR was performed using a 20-μL reaction volume. For one reaction, the assay was optimized to contain 0.5 μL of PCR product of the first PCR, 1× PCR buffer, 0.5 μL of each of the env-specific primers BLV-env-3 (5′-CCCACAAGGGCGGCGCCGGTTT-3′) and BLV-env-4 (5′-GCGAGGCCGGGTCCAGAGCTGG-3′) (20 pmol/μL), 2.5 mM of each dNTP, 2.5 units of i-Taq polymerase. The PCR reactions were performed using a T3000 thermocycler 48 (Biometra, Rundolf wissell-str, Germany) with parameters of 2 minutes at 94°C, followed by 30 cycles of 30 seconds at 95°C, 30 seconds at 58°C and 60 seconds at 72°C and 4 minutes at 72°C. The PCR products were analyzed by electrophoresis on a 1.5% agarose gel containing ethidium bromide (1 μg/mL).
The Gp51 full gene (903 bp) for sequencing was amplified using HotStar Hifidelity DNA Polymerase (QIAGEN, Hilgen, Germany) and a primer set (BLV-F-4807; CTGGCGTTTGYTGAAGGCCTT, BLV-gp51-R; ACTACGTCTGACCCGGGTA) following the manufacturer’s recommendation. The reaction mixture contained 3 μL of template, 1 μL of HotStar Hifidelity DNA Polymerase, 10 μL of 5X PCR buffer, 1 μL of each primer (10 pmol/μL) and 34 μL of distilled H2O. The PCR reaction was performed on a T3000 thermocycler 48 (Biometra, Rundolf wissell-str, Germany) with parameters of 5 min at 95°C, followed by 40 cycles of 15 s at 94°C, 60 s at 58°C, and 60 s at 72°C and 10 min at 72°C. The PCR products were analyzed by electrophoresis on a 1.5% agarose gel containing ethidium bromide (1 μg/mL).
Sequencing of the gp51
The gp51 PCR products were purified using a PCR purification kit (QIAGEN, Hilgen, Germany). The purified amplicons were cloned into the pDrive vector (QIAGEN, Hilgen, Germany) and sequenced by a local company (Macrogen Inc. Seoul, Korea). The sequences were analyzed using the DNAstar program (DNASTAR Inc., Ver. 5.0).
To elucidate the genetic relationships between Korean and worldwide BLV isolates, currently available BLV env sequences were downloaded from the GenBank database and combined with 17 unique Korean BLV env fragment sequences. After online alignment using MAFFT v. 7 (http://mafft.cbrc.jp/alignment/server/) , sequences with too many missing characters were manually excluded. To consider all of the possible genotypes of the worldwide BLV isolates, two different datasets were generated for phylogenetic analyses: partial gp51 (444 bp), ranging from 93 to 240 amino acids, and complete gp51 (903 bp), ranging from 1 to 301 amino acids . To estimate the genetic divergence of the gp51 gene, the Kimura-2-parameter (K2P)  of the genetic distance was calculated using MEGA v.5  for both intra- and inter-genotype distances. To compare the gp51 amino acid sequences from Korean to other isolates from other countries, protein sequences based on a BLV reference sequence (i.e., K02120 and AY995174) were obtained by MEGA.
For accurate and robust phylogenetic analyses of the BLV gp51 sequence, the Maximum likelihood (ML) and Bayesian Inference (BI) approaches were performed. To determine the most appropriate nucleotide substitution models, jModelTest v.2.1.4  was used among 88 candidates on a BioNJ tree, under the Akaike Information Criterion (AIC), corrected AIC(c), and Bayesian Information Criterion (BIC). The closest model available in the corresponding ML and Bayesian analyses were chosen for each dataset.
For the ML analyses, GARLI 2.01 , based on genetic algorithms for the ML search, and RAxML 7.4.2 , implemented in raxmlGUI v.1.3 , were used to compare the phylogenetic results. Because RAxML only allows the most complex GTR (General Time Reversible) substitution model but GARLI is available for the standard nested submodels of the GTR, the following models were selected after model estimation for each dataset: GTR + I (proportion of the invariable sites) (903 bp) and GTR + G (Gamma distribution of the rates among sites) (444 bp) for RAxML and TPM3uf + I (903 bp) and TPM2 + G (444 bp) for GARLI. Bootstrap supports were calculated using 1,000 bootstrap replicates for the four datasets, and the bootstrap trees were summarized with a 50% majority rule consensus tree by the SumTrees script (v.3.3.1) implemented in DendroPy v.3.10.1 .
For the Bayesian analyses, two programs were used. The first was MrBayes v.3.2.2 , but such analyses have been reported to produce biased posterior probabilities in certain cases (e.g., the “star-tree paradox” ). To address this issue, Phycas v.1.2.0  was used to as a second approach for the polytomy prior . As the best-fit models, HKY (Hasegawa, Kishino and Yano) + I (903 bp) and HKY + G (444 bp) were specified for each dataset in MrBayes. Two runs with four chains were carried out simultaneously for 10 million generations, and the trees were sampled every 100 generations. The convergence diagnostics were evaluated when the average standard deviation of the split frequencies of the two runs was less than 0.01. In Phycas, the same substitution models of MrBayes were used for each dataset, but with a 1 million generation. Tree samples prior to reaching a stationary posterior distribution were discarded, and the remaining samples were summarized with SumTrees.
Diagnosis and identification of Korean BLV isolates
A total of 185 samples (119 samples from farm A and 66 samples from farm B) were analyzed, and 78 samples (42.2%) were positive for BLV env DNA, as detected by nested PCR for a 444-bp fragment (37 samples out of 119 in farm A and 41 out of 66 in farm B) (Table 1). Another PCR reaction to amplify the full gp51 gene was also performed. Nineteen out of 37 positive samples from farm A and 31 out of 41 positive samples from farm B were successfully amplified. As a result of the phylogenic analysis based on the Korean isolates, the genotypes from farm A was divided into two groups, genotype 1 and genotype 3 (Figures 1 and 2). Out of the 19 samples analyzed from farm A, 15 were genotype 1 (78.9%) and 4 were genotype 3 (21.1%). We found that all 31 isolates from farm B belong to genotype 1. Of all of the cloned gp51 sequences, 46 samples were genotype 1 (92%) and 4 samples were genotype 3 (8%). All of the sequences of BLV gp51 obtained in this study can be found in Additional file 1: Table S1 and are deposited in the GenBank databases with accession numbers KP20146-KP201482.
Genotyping of the BLV gp51 gene based on global phylogenetic analyses
ML and BI trees based on the complete gp51 gene showed congruent topologies. The existence of eight genotypes of BLV was highly supported by high bootstrap and posterior probability values of more than 95 and 0.99, respectively (Figure 2). Phylogenetic analyses based on the partial 444 bp of gp51 also supported the distinct eight genotypes of BLV, which was previously described [17,18]. However, there were some discrepancies in the genotypic phylogenetic relationships between the ML and Bayesian analyses (Figure 1, Additional file 1: Figure S1). The use of different substitution models had little effect on the overall analyses of the genotype separation, regardless of the datasets. Supporting values for each node were relatively lower in ML than BI (Figures 1 and 2, Additional file 1: Figure S1).
Evolutionary relationships among the genotypes showed some correlation between the genotypes and geographic affinities. For instance, genotypes 7 and 4, mostly from European countries, have a close relationship with a node receiving less than 80% support, except for MrBayes (74/67/0.9/68: data not shown on the tree). Compared to the G4-G7 relationship, genotypes 2 and 3 are more closely related to each other, as has been previously described [17,18,33,36,51,52], with 88/85/0.92/67. However, the phylogenetic relationships between these groups were not fully resolved in the partial gp51 analyses.
Worldwide distribution of BLV genotypes based on partial gp51 sequences is summarized in Table 2. Genotype 1 is found in 10 countries (i.e., Korea, Japan, USA, Costa Rica, Argentina, Uruguay, Brazil, Iran, Australia, and Germany) and is the most prevalent genotype worldwide, whereas genotypes 4, 7, and 8 are mainly distributed in Europe and Russia. However, the second most widely distributed genotype, G4, is found more frequently in Europe and America. The prevalent BLV genotypes in Central and South America appear to have more variability, including G1, G2, G4, G5, G6, and G7. Among these, genotypes 5 and 6 were exclusively found in Argentinean, Brazilian, and Costa Rican isolates. Interestingly, genotype 3 isolates from the USA are closely related to those from Korea and Japan.
Genetic variation of the complete gp51 gene
The complete gp51 sequences in this study include 87 isolates with 903 nucleotides, of which, 142 nucleotides were parsimoniously informative. The average evolutionary divergence of the sequence pairs is 2.3%; similarly, 2.5% divergence was observed in the partial 444 bp from 109 sequences. Specifically, the minimum to maximum divergence range was 0.11% to 1.91% for genotype 1 and 0.33% to 1.58% for genotype 3. The estimated average divergence over sequence pairs within and between genotypes is described in Table 3.
More divergent isolate sequences were observed in genotype 2 (from Argentina, Brazil, and Japan), with an average at 1.42% for the complete gp51 sequence, whereas genotype 4 included diverse isolates at 1.38% for the partial 444 bp of gp51 (Table 3, diagonal line). Genotype 4 was found in 9 out of 20 countries, including Argentina, Chile, Belarus, Russia, Ukraine, Poland, Belgium, France, and Germany, and the nucleotide variation within the genotype varied from 0% to 2.79% (Figure 1, Tables 2 and 3). The average genetic distance between genotypes did not exceed 5% in either datasets, except for G5 to G3 (5.13%) and G5 to G2 (5.31%) for the complete gp51 sequence (Table 3). In general, genetic variation within genotype 1 appeared to be minimal compared to that of other genotypes.
Amino acid changes in the complete gp51 sequence
Eighty-seven deduced full-length amino acid sequences of BLV gp51 were aligned with the annotation of epitopes and functional domains as previously described [39,53-57] (Table 4). Specifically, amino acid sequences from the Korean BLV isolates were compared to 87 BLV isolates found throughout the world. Various single amino acid changes caused by a point mutation were discovered over the full length of gp51. A total of 85 amino acid variations were observed among the 301 amino acid sites, and 40 of these (47%) were observed in multiple isolates with at least two types of amino acids. Table 4 represents the distribution of the amino acid substitutions based on the parsimony-informative sites, according to genotype. Significant amino acid changes were not observed in the majority of Korean BLV isolates from genotype 1 compared to other isolates from different countries. However, it is interesting that only two isolates, GNCN-10 and 11, have 3 amino acid substitutions at 12 (Q- > P), 23 (T- > A), and 28 (C- > S), located in the leader peptide region, and 2 changes at 54 (K- > Q) and 69 (S- > L), corresponding to the conformational epitope region. In particular, one substitution at 140 (V- > A) in GNCN-1, GNCN-2, and GNCN-3 from Gyeongsangnam-do, occurs in the region of the second neutralizing domain (ND2), which is a binding zinc peptide that interacts with zinc ions and viral fusion proteins . In addition, two isolates, GBGS-6 and 7, in Gyeongsangbuk-do had another amino acid change at 134 (D- > N) in the same region. Additionally, the Korean genotype 3 isolates were identical with those from Japan and the USA, which are based on amino acid substitutions. Interestingly, amino acid substitutions were more highly variable in the regions containing the conformational epitope region (Table 4). These changes were manly found in genotypes 4, 5, 6, and 7, encompassing European, Russian, and South American isolates. Additionally, some of the amino acid changes had a trend towards geographic clustering; for instance, substitutions found at 73 (A- > P) and 121 (R- > H) were only found in Russian and European isolates and substitutions at 15 (I- > V) and 37 (S- > T) were only found in Costa Rican isolates.
The phylogenetic analysis presented in this study supports previously established data indicating that BLV has eight genotype groups. We also discovered unknown phylogenetic relationships between newly analyzed sequences (i.e., more recently unpublished sequences from 2014), including isolates from Korea. Due to the different number of isolates for each genotype and two different gp51 regions sequenced for the two datasets (i.e., 903 bp and 444 bp), the phylogenetic results and genetic divergence values for the isolates are not the same in all analyses. Nevertheless, the eight genotypes of BLV were completely identified through comprehensive phylogenetic analyses. Our results are consistent with previous studies [17,18].
Maximum likelihood and Bayesian approaches showed that a complex of several BLV genotypes has worldwide distribution, and phylogenetic separation of the eight major genotypes was well resolved by all of the analyses. Additionally, we found that the phylogenetic relationships among the genotypes are correlated to the geographical origin of the isolates (e.g., G4-G7 in Europe/Russia, G5-G6 in South America) based on the complete gp51 sequence. As shown in Figures 1 and 2, genotype 1 is the most dominant BLV type in the world and is found in Australia, Iran, USA, Argentina, Brazil, Uruguay, Costa Rica, Japan, and Germany. Genotype 3 included only isolates from three geographically distinct countries, Japan, the USA, and Korea. Similarly, genotypes 5 and 6 were exclusively found in three countries, Costa Rica, Argentina, and Brazil in Central and South America. In addition, isolates of genotypes 8 and 7 mostly originated from Eurasian areas, excluding Korea and Japan. Noticeably, genotypes 1 and 4 covered large geographic areas from Europe to America, suggesting the possibility of extensive cattle trading between countries (Table 2).
Phylogenetic analyses of the complete gp51 gene revealed that 17Korean BLV isolates belonged to genotypes 1 and 3 (Additional file 1: Figure S1 and S2, Additional file 1: Figure S1). Seven genotype 1 isolates were found in Gyeongsangbuk-do (GBGS) and 8 isolates were found in Gyeongsangnam-do (GNCN). Indeed, genotype 1 was reported in Korea in 2009  and was expected to be the dominant form found in our study. Two isolates from the Gyeongsangbuk-do region were assigned to genotype 3, a result that was also seen in a recent Japanese study . Considering the geographic affinity of the two countries, these results are not unusual.
As shown in Table 3, the average distance between genotype 1 isolates (0.74% for the complete gp51 sequence and 0.83% for the partial gp51 sequence) was relatively lower than those observed within G2-G8. The low level of genetic diversity did not reflect geographic barriers, but perhaps indicates a single origin of BLV infection or international importation of cattle across the border from several countries.
Analyses of the BLV gp51 amino acid sequences in isolates collected from multiple geographic locations showed specific sequence conservation depending on the genotype (Table 4). Considering this result, a few nucleotide mutations or amino acid substitutions may affect the ability of BLV to survive and may infect hosts. The gp51 region characterized by amino acid changes of several functional sites remained unchanged over a long period in this study (e.g., 216/301 conserved region (72%)) (Table 4). Based on Korean genotype 1 BLV isolates, several amino acid changes were observed at positions 12, 23, and 28 in the leader peptide region and at positions 54 and 69 in the conformational epitope region in only two isolates from the GNCN locality of Gyeongsangnam-do. In addition, unique antigenic variations were identified in the second neutralizing domain, position 140 of GNCN-1, 2, and 3 along with position 134 of GBGS-6 and 7. The biological significance of these changes in the Korean isolates needs to be investigated by functional studies in vivo, but the changes in the conformational epitope overlapping with the second neutralizing domain region indicate a possibility of them affecting viral infectivity and formation.
This study shows that two genotypes of BLV are circulating throughout Korea. It is notable that most of these isolates belong to genotype 1 (46 out of 50 isolates), which is referred to as the dominant BLV genotype around the world, whereas only 4 isolates belong to genotype 3. Further study to identify the diversity of BLV circulating in Korea is necessary with samples collected nationwide because we identified genotype 3 in Korean isolates. Two distinct genotypes imply that there were two independent introduction events in Korea. At this point, the exact geographic origin of the Korean isolates remains uncertain, even though they are highly similar (99%) to Costa Rican (EF065640) and Japanese (LC007985) genotype 1 isolates and highly similar (99%) to Japanese (EF065650) and USA (EF065648) genotype 3 isolates.
Burny A, Cleuter Y, Kettmann R, Mammerickx M, Marbaix G, Portetelle D, et al. Bovine leukaemia: facts and hypotheses derived from the study of an infectious cancer. Vet Microbiol. 1988;17(3):197–218.
Kettmann R, Cleuter Y, Mammerickx M, Meunier-Rotival M, Bernardi G, Burny A, et al. Genomic integration of bovine leukemia provirus: comparison of persistent lymphocytosis with lymph node tumor form of enzootic. Proc Natl Acad Sci U S A. 1980;77(5):2577–81.
Kettmann R, Portetelle D, Mammerickx M, Cleuter Y, Dekegel D, Galoux M, et al. Bovine leukemia virus: an exogenous RNA oncogenic virus. Proc Natl Acad Sci U S A. 1976;73(4):1014–8.
Mirsky ML, Olmstead CA, Da Y, Lewin HA. The prevalence of proviral bovine leukemia virus in peripheral blood mononuclear cells at two subclinical stages of infection. J Virol. 1996;70(4):2178–83.
Acaite J, Tamosiunas V, Lukauskas K, Milius J, Pieskus J. The eradication experience of enzootic bovine leukosis from Lithuania. Prev Vet Med. 2007;82(1–2):83–9.
Alfonso R, Almansa JE, Barrera JC. [Serological prevalence and evaluation of the risk factors of bovine enzootic leukosis in the Bogota savannah and the Ubate and Chiquinquira Valleys, Colombia]. Rev Sci Tech. 1998;17(3):723–32.
D’Angelino JL, Garcia M, Birgel EH. Epidemiological study of enzootic bovine leukosis in Brazil. Trop Anim Health Prod. 1998;30(1):13–5.
Jacobs RM, Pollari FL, McNab WB, Jefferson B. A serological survey of bovine syncytial virus in Ontario: associations with bovine leukemia and immunodeficiency-like viruses, production records, and management practices. Can J Vet Res Revue Canadienne de Recherche Veterinaire. 1995;59(4):271–8.
Marin C, de Lopez NM, Alvarez L, Lozano O, Espana W, Castanos H, et al. Epidemiology of bovine leukemia in Venezuela. Annales de Recherches Veterinaires Ann Vet Res. 1978;9(4):743–6.
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.
Nuotio L, Rusanen H, Sihvonen L, Neuvonen E. Eradication of enzootic bovine leukosis from Finland. Prev Vet Med. 2003;59(1–2):43–9.
Trono KG, Perez-Filgueira DM, Duffy S, Borca MV, Carrillo C. Seroprevalence of bovine leukemia virus in dairy cattle in Argentina: comparison of sensitivity and specificity of different detection methods. Vet Microbiol. 2001;83(3):235–48.
VanLeeuwen JA, Keefe GP, Tremblay R, Power C, Wichtel JJ. Seroprevalence of infection with Mycobacterium avium subspecies paratuberculosis, bovine leukemia virus, and bovine viral diarrhea virus in maritime Canada dairy cattle. Can Vet J La Revue Veterinaire Canadienne. 2001;42(3):193–8.
Wang CT. Bovine leukemia virus infection in Taiwan: epidemiological study. J Vet Med Sci Jpn Soc Vet Sci. 1991;53(3):395–8.
Johnston ER, Albritton LM, Radke K. Envelope proteins containing single amino acid substitutions support a structural model of the receptor-binding domain of bovine leukemia virus surface protein. J Virol. 2002;76(21):10861–72.
Mamoun R, Morisson M, Rebeyrotte N, Busetta B, Couez D, Kettmann R, et al. Sequence variability of bovine leukemia virus env gene and its relevance to the structure and antigenicity of the glycoproteins. J Virol. 1990;64(9):4180–8.
Balic D, Lojkic I, Periskic M, Bedekovic T, Jungic A, Lemo N, et al. Identification of a new genotype of bovine leukemia virus. Arch Virol. 2012;157(7):1281–90.
Rola-Luszczak M, Pluta A, Olech M, Donnik I, Petropavlovskiy M, Gerilovych A, et al. The molecular characterization of bovine leukaemia virus isolates from Eastern Europe and Siberia and its impact on phylogeny. PLoS One. 2013;8(3):e58705.
Choi W. Survey for antibodies to bovine leukemia virus in dairy and Korean native cattle. Korean J Vet Res. 1982;22:23–6.
Suh G, Lee C, Lee C, Hur T, Kang S, Son D, et al. Prevalence of anti-bovine leukemia virus antibodies in dairy and Korean native cattle. J Vet Clin. 2003;20(2):172–6.
Jung K, Shim H, Baek J. Investigation of bovine leukemia virus infection in dairy farms of northern Gyeonggi province, Korea. Korean J Vet Serv. 2012;35(4):333–7.
Suh G, Lee J, Lee C, Hur T, Son D, Ahn B, et al. Establishment of a bovine leukemia virus-free dairy herd in Korea. J Vet Sci. 2005;6(3):227–30.
Chu K, Hyong S, Im J, Seo L. Seroprevalence of infection with Neospora caninum, Mycobacterium paratuberculosis, bovine leukosis and Brucella abortus of dairy cattle in Jeonbuk-Iksan area. Korean J Vet Serv. 2007;30(1):95–102.
Lim S, Jeong W, Tark D, Yang D, Kweon C. Agar gel immunodiffusion analysis using baculovirus-expressed recombinant bovine leukemia virus envelope glycoprotein (gp51/gp30T-). J Vet Sci. 2009;10(4):331–6.
Coulston J, Naif H, Brandon R, Kumar S, Khan S, Daniel RC, et al. Molecular cloning and sequencing of an Australian isolate of proviral bovine leukaemia virus DNA: comparison with other isolates. General Virol. 1990;71(Pt 8):1737–46.
Felmer R, Munoz G, Zuniga J, Recabal M. Molecular analysis of a 444 bp fragment of the bovine leukaemia virus gp51 env gene reveals a high frequency of non-silent point mutations and suggests the presence of two subgroups of BLV in Chile. Vet Microbiol. 2005;108(1–2):39–47.
Hemmatzadeh F. Sequencing and phylogenetic analysis of gp51 gene of bovine leukaemia virus in Iranian isolates. Vet Res Commun. 2007;31(6):783–9.
Molteni E, Agresti A, Meneveri R, Marozzi A, Malcovati M, Bonizzi L, et al. Molecular characterization of a variant of proviral bovine leukaemia virus (BLV). Zentralblatt fur Veterinarmedizin Reihe B J Vet Med Series B. 1996;43(4):201–11.
Monti G, Schrijver R, Beier D. Genetic diversity and spread of Bovine leukaemia virus isolates in Argentine dairy cattle. Arch Virol. 2005;150(3):443–58.
Zhao X, Buehring GC. Natural genetic variations in bovine leukemia virus envelope gene: possible effects of selection and escape. Virology. 2007;366(1):150–65.
Camargos MF, Pereda A, Stancek D, Rocha MA, dos Reis JK, Greiser-Wilke I, et al. Molecular characterization of the env gene from Brazilian field isolates of Bovine leukemia virus. Virus Genes. 2007;34(3):343–50.
Camargos MF, Stancek D, Rocha MA, Lessa LM, Reis JK, Leite RC. Partial sequencing of env gene of bovine leukaemia virus from Brazilian samples and phylogenetic analysis. Vet Med B Infect Dis Vet Publ Health. 2002;49(7):325–31.
Rodriguez SM, Golemba MD, Campos RH, Trono K, Jones LR. Bovine leukemia virus can be classified into seven genotypes: evidence for the existence of two novel clades. J General Virol. 2009;90(Pt 11):2788–97.
Licursi M, Inoshima Y, Wu D, Yokoyama T, González ET, Sentsui H. Provirus variants of bovine leukemia virus in naturally infected cattle from Argentina and Japan. Vet Microbiol. 2003;96(1):17–23.
Asfaw Y, Tsuduku S, Konishi M, Murakami K, Tsuboi T, Wu D, et al. Distribution and superinfection of bovine leukemia virus genotypes in Japan. Arch Virol. 2005;150(3):493–505.
Matsumura K, Inoue E, Osawa Y, Okazaki K. Molecular epidemiology of bovine leukemia virus associated with enzootic bovine leukosis in Japan. Virus Res. 2011;155(1):343–8.
Fechner H, Blankenstein P, Looman AC, Elwert J, Geue L, Albrecht C, et al. Provirus variants of the bovine leukemia virus and their relation to the serological status of naturally infected cattle. Virology. 1997;237(2):261–9.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
Moratorio G, Fischer S, Bianchi S, Tomé L, Rama G, Obal G, et al. A detailed molecular analysis of complete bovine leukemia virus genomes isolated from B-cell lymphosarcomas. Vet Res. 2013;44:19.
Kimura M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.
Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.
Zwickl DJ. Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion [Dissertation]. Austin: The University of Texas at Austin; 2006.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.
Silvestro D, Michalak I. raxmlGUI: a graphical front-end for RAxML. Org Divers Evol. 2011;12(4):335–7.
Sukumaran J, Holder MT. DendroPy: a Python library for phylogenetic computing. Bioinformatics. 2010;26(12):1569–71.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42.
Yang Z, Rannala B. Branch-length prior influences Bayesian posterior probability of phylogeny. Syst Biol. 2005;54(3):455–70.
Phycas 1.2.0 [http://www.phycas.org]
Lewis PO, Holder MT, Holsinger KE. Polytomies and Bayesian phylogenetic inference. Syst Biol. 2005;54(2):241–53.
Moratorio G, Obal G, Dubra A, Correa A, Bianchi S, Buschiazzo A, et al. Phylogenetic analysis of bovine leukemia viruses isolated in South America reveals diversification in seven distinct genotypes. Arch Virol. 2010;155(4):481–9.
D’Angelino RH, Pituco EM, Villalobos EM, Harakava R, Gregori F, Del Fava C. Detection of bovine leukemia virus in brains of cattle with a neurological syndrome: pathological and molecular studies. BioMed Res Int. 2013;2013:425646.
Dube S, Abbott L, Dube DK, Dolcini G, Gutierrez S, Ceriani C, et al. The complete genomic sequence of an in vivo low replicating BLV strain. Virol J. 2009;6:120.
Dube S, Dolcini G, Abbott L, Mehta S, Dube D, Gutierrez S, et al. The complete genomic sequence of a BLV strain from a Holstein cow from Argentina. Virology. 2000;277(2):379–86.
Rice NR, Stephens RM, Couez D, Deschamps J, Kettmann R, Burny A, et al. The nucleotide sequence of the env gene and post-env region of bovine leukemia virus. Virology. 1984;138(1):82–93.
Sagata N, Yasunaga T, Tsuzuku-Kawamura J, Ohishi K, Ogawa Y, Ikawa Y. Complete nucleotide sequence of the genome of bovine leukemia virus: its evolutionary relationship to other retroviruses. Proc Natl Acad Sci U S A. 1985;82(3):677–81.
Willems L, Thienpont E, Kerkhofs P, Burny A, Mammerickx M, Kettmann R. Bovine leukemia virus, an animal model for the study of intrastrain variability. J Virol. 1993;67(2):1086–9.
Gatot JS, Callebaut I, Van Lint C, Demonte D, Kerkhofs P, Portetelle D, et al. Bovine leukemia virus SU protein interacts with zinc, and mutations within two interacting regions differently affect viral fusion and infectivity in vivo. J Virol. 2002;76(16):7956–67.
Derse D, Diniak AJ, Casey JW, Deininger PL. Nucleotide sequence and structure of integrated bovine leukemia virus long terminal repeats. Virology. 1985;141:162–6.
Stancek D, Camargos MF, Reis JKP, Rocha MA, Lessa LM, Leite FSC, et al. Higher frequency of non-silent point mutations in some bovine leukemia virus (BLV) strains isolated in Brazil in comparison with the strains circulating in several other countries. J Vet Med B Infect Dis Vet Public Health. 2002. In press 2002.
This study was funded by Ministry of Agriculture, Food and Rural Affairs, Republic of Korea to Animal and Plant Quarantine Agency (I-1543083-2013-13-01).
The authors declare that they have no competing interests.
Conceived and designed the experiments: YKS EJK EJL. Analyzed the data: EJL HKJ Contributed reagents/materials/analysis tools: YKS EJK. Wrote the paper: EJL YKS. Contributed sampling: YKS EJK. Performed the laboratory experiments: HKJ BHK. Reviewed manuscript draft and provided editorial comments: JYS ISC KKL. Supervised the study: YKS. All authors read and approved the final manuscript.
Figure S1. A Bayesian inference phylogenetic tree based on the partial gp51 (444 bp) sequences from different geographic regions. Korean isolates are shown in blue bold-italic names. The remaining isolates in the tree are denoted by country of origin, author with published date (or direct submitted date), and accession number. Genotypes 1 through 8 are indicated by vertical lines with the symbol ‘G’. Two numbers at the branches indicate BI posterior probabilities values (MrBayes(0–1)/Phycas(0–100)). An asterisk indicates unpublished or direct submission to GenBank sequences. + denotes that the selected sequences were used for the complete gp51 analyses. *Means that the sequences have not yet been investigated by a published phylogenetic analysis. The tree is rooted on genotype 5. Table S1. Accession number and geographic information regarding the fifty Korean BLV isolate sequences used in this study.