New evidence of bovine leukemia virus circulating in Myanmar cattle through epidemiological and molecular characterization

Kyaw Kyaw Moe; Meripet Polat; Liushiqi Borjigin; Ryosuke Matsuura; Si Thu Hein; Hla Hla Moe; Yoko Aida

doi:10.1371/journal.pone.0229126

Abstract

Bovine leukemia virus (BLV) is the etiological agent of enzootic bovine leukosis, which is the most common neoplastic disease of cattle. BLV infects cattle worldwide and causes serious problems for the cattle industry. In this study, we examined the prevalence of BLV infection and the distribution of BLV genotypes in cattle in the northern, central, and southern parts of Myanmar. The prevalence of BLV infection among Myanmar cattle (37.04%) in this study was markedly higher than the prevalence (9.1%) observed in our earlier study in which BLV was detected from the limited number of cattle only from a small area of Myanmar. Phylogenetic analysis of partial env-gp51 sequence of the isolated BLV strains revealed that there are at least three BLV genotypes (genotype-1, genotype-6, and genotype-10) in Myanmar, which have also been detected in the neighboring countries. We performed this study to estimate the BLV proviral load, which is a major diagnosis index for determining the virus transmission risk. The cattle of the three test regions with warm, wet, and humid climatic conditions (upper Sagaing, Yangon, and Kayin) exhibited a high mean proviral load, while cattle of three other regions with low annual rainfall and very high temperature (Mandalay, Magway, and upper Bago) exhibited a low mean proviral load. Further, the level of proviral load and the prevalence of BLV infection in Myanmar native cattle (N = 235) were lower than that in the hybrid cattle (Holstein Friesian × Myanmar native) (N = 62). We also observed that the cattle with high risk for BLV transmission, which have high proviral load, may enhance the BLV infection rate. Hence, to control BLV transmission, it is necessary to eliminate these cattle with high-risk for BLV transmission and to diagnose BLV provirus in cattle in the remaining regions/states of Myanmar sharing a boundary with neighboring countries.

Citation: Moe KK, Polat M, Borjigin L, Matsuura R, Hein ST, Moe HH, et al. (2020) New evidence of bovine leukemia virus circulating in Myanmar cattle through epidemiological and molecular characterization. PLoS ONE 15(2): e0229126. https://doi.org/10.1371/journal.pone.0229126

Editor: William M. Switzer, Centers for Disease Control and Prevention, UNITED STATES

Received: April 18, 2019; Accepted: January 30, 2020; Published: February 21, 2020

Copyright: © 2020 Moe et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: The study was supported by the Grants-in-Aid for Scientific Research (A) from the Japan Society for the Promotion of Science (JSPS) (Grant No. 16H02590). This research was also supported by grants from the Project of the NARO Bio-oriented Technology Research Advancement Institution [the Special Scheme Project on Regional Developing Strategy (Grant No. 16817983) and the Special Scheme Project on Vitalizing Management Entities of Agriculture, Forestry and Fisheries (Grant No. 16930548)]. Kyaw Kyaw Moe is a JSPS International Research Fellow [FY2018 Invitational Fellowship for Research in Japan (Long-term)].

Competing interests: The authors have declared that no competing interests exist.

Introduction

Bovine leukemia virus (BLV) infects cattle worldwide and is the etiologic agent for enzootic bovine leukosis (EBL). BLV belongs to the genus Deltaretrovirus of the Retroviridae family [1, 2]. BLV is closely related to the human T-cell leukemia virus type 1 and 2 (HTLV-1 and -2) [1]. BLV has also been detected in the human tissues, such as human breast tissue, which suggests a risk for the infection and proliferation of this virus in humans [3].

Majority of the BLV-infected cattle are asymptomatic. However, approximately one-third of the infected animals develop persistent lymphocytosis (PL), which is characterized by the polyclonal proliferation of CD5+ B lymphocytes [1]. Only 1–5% of the BLV-infected cattle exhibit malignant monoclonal B-cell lymphosarcoma. The clinical symptoms of lymphoma, caused by BLV infection, depend on the site of the tumor, but include digestive disturbance, loss of appetite, weight loss, weakness, reduction in milk production or general debility, and various neurological manifestations [4]. In addition, recent studies have reported that BLV infection reduces milk production [5, 6] and is associated with high incidence of infectious disease [7] and reproductive inefficiency, resulting in high culling rates of the BLV-infected but healthy cattle [8]. Hence, BLV eradication is of utmost importance. BLV is transmitted through blood containing infected lymphocytes and through the colostrum and milk. Additionally, farm practices such as tattooing, dehorning, rectal palpation, and injection have also been associated with BLV transmission [1]. Blood-sucking flies are one of the important factors for BLV transmission [9, 10]. BLV infection rates tend to be high in older cattle populations [11].

The BLV genome comprises the essential genes gag, pro, pol, and env, which encode structural proteins and enzymes, the regulatory genes tax and rex, and the accessory genes R3 and G4, and two identical long terminal repeats (LTRs) [1, 12]. Similar to the other members of the Retroviridae family, the BLV env gene encodes the envelope glycoprotein precursor pr72env, which is important for viral infection and syncytium formation. BLV env gene also encodes a mature surface glycoprotein (gp51) and a transmembrane protein (gp30) [13–15]. The glycoprotein gp51 plays a crucial role in the viral life cycle and influences the capacity of BLV to enter cells and has been identified as a target of specific neutralizing antibodies [16–18]. Therefore, the gp51 region has been widely used for BLV genotyping and phylogenetic studies for the identification of viral strains isolated from different geographical regions. Globally, at least 11 different genotypes have been identified [19–36].

Previous studies have demonstrated that the proviral load, which represents the BLV genome integrated into the host genome, is an important index for estimating the BLV infection stage. Proviral load is associated with disease progression [37, 38, 39], lymphocyte count [40], viral biokinetics [41], viral infectivity [42], and virus shedding into the saliva and nasal secretions [43]. Among the cows with high risk for BLV transmission, the cows in which BLV provirus was detected in the nasal secretion and saliva had a proviral load of over 14,000 copies/10⁵ cells and 18,000 copies/10⁵ cells in their blood sample, respectively [43]. It has been suggested that these cows are strongly responsible for BLV transmission in healthy cows via direct contact.

Myanmar, situated in the Southeast Asian mainland, is an agro-based country. Myanmar native cattle are zebu type (Bos indicus) and is the main draught animal used in cultivation. They are also used for transportation in the rural area and for meat and milk production [44]. Local breeds of cattle, such as Pyar Sein and Shwe Ni, are widely reared throughout the country. These breeds are adapted to the harsh native environment and are resistant to tropical diseases and external parasites. For the development of local draught cattle, zebu type Shindi and Thari breeds were imported into Myanmar in the 1950s and were crossed with the local breeds. In the 1970s, the local breeds in Myanmar were crossed with the imported dairy cattle breeds such as Holstein Frisian and Jersey to improve their dairy performance. Currently, dairy cattle play a secondary role in ruminant production owing to demand for milk and milk products among the urban populations and improved nutritional standards and quality of the life of the rural populations. In 2017, Myanmar owned 17.14 million heads of cattle (about 1.1% of the world cattle population), based on the FAOSTAT database (http://www.fao.org/faostat/en/#data/QA).

Recently, we have detected BLV infection in different cattle farms of the Nay Pyi Taw Union Territory, Myanmar, and examined the genetic variability in the BLV strains using blood samples of 66 cattle [33]. About 9.1% of the samples were positive for BLV provirus and a phylogenetic tree, constructed using partial and complete env-gp51 sequences, indicated that Myanmar strains were genotype-10 among the 11 globally detected genotypes. However, only a limited number of samples were tested in the previous study, which were collected from a small and definite area. If samples are collected from different parts of the country, more genotypes might be detected. Moreover, there is a need to examine the diversity of BLV in Myanmar, as three BLV genotypes (genotype-1, genotype-6, and genotype-10) have been recorded in the countries neighboring Myanmar; genotype-1, genotype-6, and genotype-10 in Thailand [32], genotype-6 and genotype-10 in China [45, 46], and genotype-6 in India [36]. In this study, we examined the detailed epidemiological information on the prevalence of BLV in Myanmar using BLV-CoCoMo-quantitative polymerase chain reaction (qPCR) assay [39, 47]. We extracted DNA from 297 bovine blood samples, which were collected from different locations throughout the country. Furthermore, we demonstrated that the level of BLV proviral load was associated with different regional variations and cattle breeds in Myanmar. Additionally, we investigated the genetic diversity of BLV strains detected in Myanmar by phylogenetic analysis of BLV env-gp51 gene sequence. We detected the presence of at least three BLV genotypes (genotype-1, genotype-6, and genotype-10) in Myanmar.

Materials and methods

Animals and extraction of genomic DNA

Blood samples were collected from 297 cattle in 42 different private farms of 11 townships located in 5 regions and 1 state in the northern, central, and southern parts of Myanmar between February and March 2018, as shown in Fig 1. The cattle breeds included native Pyar Sein and Shwe Ni breeds (N = 235), and crossbred Holstein-Friesian breed (N = 62). The age of the cattle ranged from 6-months to 18-years. Most of the farms were small holding farms having less than 10 cattle per farm. The permission from owner was required to access farms and cattle. The genomic DNA was extracted from 200 μL of whole bovine blood using PureLink Genomic DNA Mini Kit (Invitrogen, Carlsbad, CA), following manufacturer’s instructions.

Download:

Fig 1. Map of Myanmar showing the number of cattle included in this study and bovine leukemia virus (BLV) genotypes detected in the tested area.

Blood sampling was performed in eleven townships (indicated by black color circles) located in one state and five regions of Myanmar. The gray color indicates the blood sampling area in previous study. The figure is created by modifying "map of Myanmar", copyright belongs to (© Keisuke Inoue (Licensed under CC BY 4.0 https://www.freemap.jp/item/asia/myanmar.html).

https://doi.org/10.1371/journal.pone.0229126.g001

Evaluation of BLV provirus by BLV-CoCoMo-qPCR-2.

BLV-CoCoMo-qPCR-2 (RIKEN Genesis, Kanagawa, Japan) using THUNDERBIRD Probe qPCR Mix (Toyobo, Tokyo, Japan) was used to determine the presence of BLV and to measure proviral load as described previously [47]. A 183-bp sequence of the BLV LTR was amplified using the degenerate primer pair, CoCoMo-FRW and CoCoMo-REV, and a 15 bp 6-carboxyfluorescein (FAM)-labeled LTR probe. A 151-bp sequence of BoLA-DRA (internal control) was amplified using the primer pair DRA-FW and DRA-RW and a FAM-labeled DRA probe as previously described [47]. The proviral load was calculated as follows: (number of BLV LTR copies/number of BoLA-DRA copies) x 10⁵ [47].

PCR amplification of BLV env-gp51 gene fragments and nucleotide sequences

The partial sequence of BLV env-gp51 gene of the BLV provirus-positive samples was amplified by nested PCR using PrimeSTAR GXL DNA Polymerase (Takara Bio Inc., Otsu, Japan), as described previously [28, 48, 49]. The BLV env-gp51-positive PCR products were purified using Exo-SAP IT (USB Corp., Cleveland, OH) and were sequenced on an ABI 3730xl DNA Analyzer using the BigDye Terminator v 3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA). The sequence included a 502 bp region of the env-gp51 gene, corresponding to the 5078 to 5579 nucleotide position on the BLV whole genome sequence (National Center for Biotechnology Information (NCBI) GenBank accession number K02120) [50]. The nucleotide sequence was edited, aligned, and identified using MEGA 7 software [51].

Construction of the phylogenetic tree

The partial sequence of BLV env-gp51 from 66 BLV-positive DNA samples, which were successfully amplified by nested PCR, was aligned with the 44 BLV env-gp51 sequences from GenBank (representative of the eleven known BLV genotypes) using MEGA 7 software. Phylogenetic analysis of a partial (502 bp) sequence of the env-gp51 gene was conducted using MEGA 7. For robust and accurate phylogenetic analysis of the BLV env-gp51 sequence, the ‘‘find best DNA/Protein models” tool of MEGA 7 software was used to choose the best fit model. The Kimura-2 parameter model with gamma distribution (K2+G) was chosen as the model with the best fit to analyze the BLV env-gp51 sequence with the smallest Akaike information criterion (AIC) values. Two phylogenetic trees were constructed using the maximum likelihood (ML) algorithm with the K2+G model of nucleotide substitution in MEGA 7. The reliability of the phylogenetic relationships was evaluated by nonparametric bootstrap analysis with 1000 replicates. The deduction of protein sequence through translation of nucleotide to amino acid sequence was performed using MEGA 7 [51].

Statistical analysis

Brunner-Munzel test was used to calculate the significance of the differences of the means of proviral loads between two breeds. Tukey's test after the analysis of variance was used to determine the significance of the means of proviral loads between the different geographical regions. P<0.05 was considered significant.

Animal handling and research ethics

All animals were handled by RIKEN, Japan in strict accordance with good animal practice following the guidelines of RIKEN. The study was approved by the RIKEN Animal Experiments Committee (approval number H29-2-104).

Results

Infection rate prevalence of BLV

Previously, we had analyzed BLV infection among the cattle of different farms in the Nay Pyi Taw Union Territory of central Myanmar as highlighted in gray in Fig 1 and we observed an infection rate of 9.1% [33]. In this study, to investigate the prevalence of BLV infection in Myanmar, we collected 297 blood samples of cattle in different farms from 11 townships located in 5 regions and 1 state, which are flanked by Nay Pyi Taw Union Territory (Fig 1). The genomic DNA isolated from these samples was first screened for BLV infection by BLV-CoCoMo-qPCR-2 targeting the BLV LTRs. Surprisingly, BLV infection was detected in all the regions and state tested in this study, in addition to the Nay Pyi Taw Union Territory. Out of a total 297 cattle DNA samples from 42 farms throughout Myanmar, 110 samples (37.04% with a 95% confidence interval extends from 31.74% to 42.66%) from 31 farms were positive for BLV provirus, as determined by BLV-CoCoMo-qPCR-2 (Table 1 and Fig 2A). This indicated that the BLV infection rate was higher than that previously observed in the Nay Pyi Taw Union Territory.

Download:

Fig 2.

Proviral load and prevalence of bovine leukemia virus (BLV) infection among the cattle of six test areas in Myanmar (A) and among the hybrid (Holstein-Friesian crossbred with native breed) cattle and among the native cattle breed (B). The mean proviral load is indicated for each region/state and each cattle breed. Significant differences between two cattle breeds and between the different geographical regions were calculated by Brunner-Munzel test and Tukey's test, respectively. P<0.05 was considered significant.

https://doi.org/10.1371/journal.pone.0229126.g002

Download:

Table 1. Summary of bovine leukemia virus (BLV) detection in Myanmar cattle as determined by CoCoMo-quantitative PCR (qPCR).

https://doi.org/10.1371/journal.pone.0229126.t001

Estimation of BLV proviral load

In BLV-CoCoMo-qPCR assay, we amplified a single-copy host gene, BoLA-DRA gene to calculate cell number, in parallel with the viral genomic DNA, and calculated the proviral load (expressed as the number of copies of provirus per 100,000 peripheral blood mononuclear cells [PBMCs]). Then we summarized the level of proviral load with BLV infection rate in each tested region/state of Myanmar, as shown in Fig 2A. The BLV-positive farms in two regions (Sagaing and Yangon) and a state (Kayin) exhibited a high mean proviral load (ranging from 4,663 to 9,665 copies), while that in the other three regions (Mandalay, Magway, and Bago) exhibited a low mean proviral load (ranging from 125 to 1,958 copies). Interestingly, Yangon region, where the highest mean proviral load (mean 9665 copies) was observed, exhibited the highest prevalence of BLV infection (58.70%) among the 6 test areas. Particularly, the highest mean proviral load (mean 22,576 copies) was observed among the dairy cattle in farm YG7 located in the Hmawbi Township of Yangon region. The copy number in this farm ranged from 195 to 57,870 copies per 10⁵ cells. Additionally, this farm exhibited the highest infection rate (88.89%) among the 42 farms (Table 1). Seven out of the eight farms tested (87.5%) were positive for BLV provirus in the Yangon region (Table 1). The second highest mean proviral load (9,082 copies) among the 6 test areas was observed in Kayin state. Out of the three cattle farms in Myawaddy Township of Kayin state, the farm KY2 tested positive for BLV provirus (45.0%), with the copy number ranging from 72 to 29,490 copies per 10⁵ cells. The third highest mean proviral load (4,663 copies) among the 6 test areas was observed in Sagaing region. All the six farms in this region tested positive for BLV provirus (100%). The highest mean proviral load (9,000 copies) and the highest prevalence of BLV infection (55.0%) in Sagaing region was observed among the dairy cattle of farm SG1 in Tamu Township with the proviral load ranging from 263 to 44,260 copies per 10⁵ cells (Table 1). In contrast to the cattle of Sagaing, Yangon and Kayin, the rate of BLV-positive farm in Mandalay, Magway, and Bago regions, which had low mean proviral load, was 83.3%, 63.3%, and 66.7%, respectively.

We then compared the prevalence of BLV infection between the native cattle of Myanmar (Pyar Sein and Shwe Ni) (N = 62) and the hybrid cattle (Holstein Friesian × Myanmar native) (N = 235) (Fig 2B). The BLV infection prevalence among the hybrid cattle was 54.84% (34 out of 62 heads), with the copy number ranging from 46 to 57,870 copies per 10⁵ cells (mean 10,420 copies). The prevalence of BLV infection among Myanmar native cattle was 32.34% (76 out of 235 heads), with the copy number ranging from 21 to 36,467 copies per 10⁵ cells (mean 2,367 copies) (Fig 2B). These results demonstrated that the prevalence of BLV infection and the level of proviral load among Myanmar native cattle tend to be significantly lower than those observed among the hybrid cattle. Nonetheless, the difference in sample size is quite large (62 vs 235), and the hybrid cattle is unevenly distributed among 6 different geographical regions, additional studies with more animals each geographical region are required to further confirm these findings.

Phylogenetic analysis of the partial sequence of BLV env-gp51

A previous phylogenetic analysis based on partial and complete env-gp51 sequence of the BLV strains isolated from only the Nay Pyi Taw Union Territory revealed that the presence of BLV genotype-10 in Myanmar [33]. In this study, to gain further insight into the degree of genetic variability among the BLV strains in Myanmar, phylogenetic characterization was carried out after sequencing the env-gp51 gene. We selected 83 out of the 110 BLV-positive DNA samples as determined by BLV-CoCoMo-qPCR-2, 66 samples were successfully amplified for the BLV env-gp51 gene using nested PCR (Table 1). This result showed the same tendency as discussed in our previous publications [28, 39, 52], indicating that BLV-CoCoMo-qPCR-2 has better sensitivity as compared with nested PCR. Indeed, Table 2 shows that proviral load of positive cows (N = 66), as determined by nested PCR for BLV env-gp51 gene, were detected at a significantly higher level as compared with that of cows negative for provirus (N = 17) as determined by with nested PCR. Then the 66 DNA samples were sequenced and subjected to phylogenetic analysis. The 502 bp sequence corresponding to the nucleotide positions 5,078 to 5,579 of the full-length BLV genome was aligned with the 44 nucleotide sequences from known BLV strains representing all existing BLV genotypes (genotype-1 to genotype-11). The similarity of the 502 bp env-gp51 nucleotide sequence ranged from 95.2 to 100% among the 66 Myanmar BLV strains used in this study (S1 Table). The degree of similarity of the nucleotide sequence with the 44 nucleotide sequences of 11 known genotypes deposited in GenBank ranged from 92.9 to 100% (S2 Table).

Download:

Table 2. Association between proviral load as determined by BLV-CoCoMo-qPCR-2 and BLV infection rate as determined by nested PCR for BLV env-gp51 gene [p Value = 0.000454, odds ratio (95%CI) = 7.125 (2.18–23.3)].

https://doi.org/10.1371/journal.pone.0229126.t002

We constructed a ML phylogenetic tree using the K2+G model of nucleotide substitution [51], as shown in Fig 3. The results of phylogenetic analysis were similar to those published in the previous studies and BLV strains were classified into 10 genotypes. The 66 Myanmar BLV strains isolated in this study showed high homology that 26 distinctive sequences of Myanmar BLV strains were used in the phylogenetic analysis. Among these distinctive sequences 4 distinctive sequences from 23 strains clustered with genotype-1. The degree of similarity of these sequences ranged from 99 to 100% when compared with the 9 nucleotide sequences of the BLV genotype-1 (S3 Table). Further, the 10 distinctive sequences of the 16 BLV strains clustered with that of genotype-6 and were 95.7 to 100% similar to the 7 nucleotide sequences of the genotype-6 (S4 Table). The 12 distinctive sequences of the remaining 27 strains clustered with that of the genotype-10 and were 97.7 to 99.8% similar to 6 nucleotide sequence of the genotype-10 (S5 Table).

Download:

Fig 3. Maximum likelihood phylogenetic tree based on 502 bp nucleotide sequence of the env-gp51 gene from 66 Myanmar bovine leukemia virus (BLV) strains and 44 known BLV strains isolated elsewhere.

Myanmar BLV strains are indicated by the sample ID along with the region/state name and country name. The remaining isolates in the phylogenetic tree are indicated by accession number and country of origin. Distinctive sequences of Myanmar BLV strains are indicated by filled (λ) circles. The 26 typical BLV isolates (accession numbers LC466589–LC466614) were aligned as shown in Fig 5A and indicated by filled (λ) circles. Myanmar BLV strains isolated in our previous study are indicated by filled triangles (▲). The genotypes are indicated by the number towards the right of the figure.

https://doi.org/10.1371/journal.pone.0229126.g003

Phylogenetic analysis indicated four separate subgroups within the genotype-6: G-6a, G-6b, G-6c and G-6d (Fig 3). To confirm the presence of 4 subgroups of genotype-6, a separate ML phylogenetic tree was constructed with the sequences of Myanmar genotype-6 and genotype-10 strains together with other reference sequences of genotype-6, and genotype-10 worldwide (Fig 4). As shown in Fig 4 the sequence of Genotype-6 clustered into 6 subgroups. This subgrouping is supported by bootstrap value ranging from 76% to 99% (Fig 4) and the inter-subgenotype pairwise distance value ranging from 2.3% to 4.3% (S6 Table).

Download:

Fig 4. Maximum likelihood phylogenetic tree based on 502 bp nucleotide sequence of the env-gp51 gene of Myanmar bovine leukemia virus (BLV) strains and other strains representing genotype-6 and genotype-10.

Myanmar BLV strains isolated in this study are indicated by the sample ID along with the region/state name and country name. The remaining isolates in the phylogenetic tree are indicated by the accession number and country of origin. Myanmar BLV strains in genotype-6 and genotype-10 are indicated by filled dark circles (●) and open circles (○). The genotypes are indicated by the number towards the right of the figure. The bar at the bottom of the figure denotes the distance.

https://doi.org/10.1371/journal.pone.0229126.g004

BLV genotype distribution in Myanmar

Phylogenetic analysis revealed that there are 3 BLV genotypes in Myanmar (Figs 3 and 4), as summarized in Fig 1. BLV genotype-1 was the most prevalent genotype and was detected in all of the six tested areas of Myanmar. The genotype-10 was detected in 5 tested areas except Mandalay region. The genotype-6 was detected in the 3 regions of southern Myanmar (G6a in Yangon and Kayin, G6b in Kayin and, G6c in Bago and Yangon).

Nucleotide and amino acid substitution in the BLV env gp51 gene of BLV strains isolated in Myanmar

One typical sequence among those that showed 100% identity was submitted to GenBank (accession numbers LC466589–LC466614). Twenty six selected nucleotide sequences from the 66 Myanmar BLV strains were aligned with that of the Japanese K02120 strain (GenBank accession number K02120), a reference sequence (Fig 5A). All of these are indicated by filled circles in the ML phylogenetic trees (Figs 3 and 4).

Download:

Fig 5. Alignment of the nucleotide and deduced amino acid sequences of the partial bovine leukemia virus (BLV) env-gp51 gene.

Nucleotide sequence (A) and deduced amino acid sequence (B) of the twenty six Myanmar BLV strains were aligned with the reference sequence of genotype-1 (Japanese K02120 strain: GenBank accession number K02120). The Myanmar BLV strains identified in this study are indicated by sample ID. The genotypes are indicated by the vertical lines towards the far left of the figure. The number on the top of the nucleotide sequences indicate the first and last nucleotide positions, and the position of nucleotide substitutions. The number above the deduced amino acid sequence is the amino acid residue number that indicate the start and end of each domain and the position of amino acid substitution. The first, second, and third neutralizing domains (ND) as well as other epitopes are shown at the top of the alignment. Dots indicate identity with the reference BLV strain K02120.

https://doi.org/10.1371/journal.pone.0229126.g005

The nucleotide sequence of all BLV strains that clustered in genotype-1 are mostly similar to that of the reference strain K02120. Myanmar BLV genotype-1, genotype-6 and genotype-10 exhibited two common silent substitutions at nucleotide (nt) 363 and nt 399, corresponding to amino acid residues 121 and 133, respectively. In contrast, Myanmar BLV strains in the three subgenotypes of genotype-6 exhibited diversity in their nucleotide sequence (Fig 5A). Two strains, MW 18 and HB9 in the BLV subtype 6a, exhibited six unique silent substitutions of nt 390, 507, 516, 618, 651, 243, corresponding to amino acid residues 130, 169, 172, 206, 217 and 243, respectively. The remained strains of subtype 6c exhibited seven unique silent substitutions at nt 378, 471, 567, 585, 663, 702 and 750, corresponding to amino acid residues 126, 157, 189, 195, 221, 236 and 250.

All or most of the Myanmar BLV strains belonging to genotype-6 and genotype-10 exhibited six common silent substitutions in the third base of the residues 135, 161, 184, 185, 205, and 236, which were not affected by the variations in nt 405, 483, 552, 555, 615 and 717, respectively. All of the BLV strains that clustered in the genotype-10 exhibited a unique silent substitution in the third base of the residue 207 (nucleotide 619) (Fig 5A).

Most of the silent nucleotide substitution of Myanmar BLV strains were located within the epitopes, including the viral G epitope (nucleotide 363), the second neutralizing domain (ND) (nucleotides 390, 399, 405, and 447), the CD8+ T-cell epitope (nucleotides 471, 483, 507 and 516), the E epitope (nucleotides 552, 555 and 567), the third neutralizing domain (ND) (nucleotides 651 and 663), and the B epitope (nucleotide 708) as shown in Fig 5B.

Twenty-six deduced amino acid sequences representing the sixty-six different Myanmar BLV strains were aligned with the predicted amino acid sequence of K02120. Fig 5B shows the distribution of amino acid changes within the region of env-gp51 (amino acid positions 87 to 253), which includes a portion of the first ND (residues 97–106), second ND (residues 131–150) and third ND (residues 210–225), a portion of the CD4+ T-cell epitope (residues 87–113) and CD8+ T-cell epitope (residues 154–182), and the viral F (residues 95), G (residues 121), E (residues 175–194), B (residues 228–238) and D (residues 251–253) epitopes [26, 53]. Myanmar BLV strains exhibited highly conserved regions in the predicted amino acid sequences of the partial env-gp51 gene. The deduced amino acid sequence of the Myanmar genotype-1 strains exhibited high homology to that of the Japanese K02120 strain. Six amino acid variations were detected among the Myanmar BLV strains. Three amino acid substitutions were specifically detected in all or most of the genotype-10 strains. They were substitution of V106A within the CD4+ T-cell epitope and the first ND, substitution of S137Q and substitution of Q143R within the second ND. Furthermore, the genotype-6 strains and genotype-10 strains exhibited a common amino acid substitution of I144T within the second ND and the genotype-6c strains and genotype-10 strains exhibited amino acid substitution of P177S within a portion of CD8+T-cell epitope and the E epitope. All of the Myanmar strains of the three genotypes exhibited substitution of S252T within the D epitope (Fig 5B).

Discussion

This study focused on the epidemiologic and molecular characterization of BLV strains in Myanmar and leads to three major conclusions. Firstly, we used BLV-CoCoMo-qPCR targeting BLV LTRs to determine the BLV infection rate among cattle in Myanmar. The data clearly indicated that the prevalence of BLV infection among Myanmar cattle was 37.04% with a 95% confidence interval extends from 31.74% to 42.66%, which is markedly higher than the prevalence of BLV infection (9.1%) reported in our earlier study. The low prevalence detected in our earlier study may be because BLV was detected from a limited number of DNA samples collected only from a small area of Myanmar [33]. Secondly, we, for the first time in Myanmar, successfully distinguished the cattle with a high proviral load from those with a low proviral load. This indicated that there was a correlation between the level of BLV proviral load and the BLV infection rate, geographic region, and cattle breeds. Third, the phylogenetic analysis of the partial env-gp51 gene sequence of the BLV strains isolated in Myanmar revealed that there were three genotypes, genotype-1, genotype-6 and genotype-10 of BLV in the country. The results of this study are more accurate than that from our earlier study on BLV strains isolated from Nay Pyi Taw Union Territory, which detected only genotype-10 in Myanmar [33]. This study also demonstrated that the three BLV genotypes, genotype-1, genotype-6 and genotype-10, were similar to those identified in the neighboring countries of Myanmar [32, 36, 45, 46]. Interestingly, the BLV strains belonging to genotype-6 clustered into 3 subgenotypes and the BLV strains belong to the subtype-6a had distinct nucleotide and amino acid substitutions.

Previously, we had shown an infection rate of 9.1% using cattle from several farms in the Nay Pyi Taw Union Territory of central Myanmar [33]. In this study, we investigated the prevalence of BLV in 11 townships in five regions and 1 state located in the northern, central and southern parts of Myanmar. We demonstrated that the BLV infection prevalence was 37.04% (110 out of 297 head). When we compared the BLV infection rate among the 6 test areas, we observed that the highest infection rate was in Yangon region (58.70%) and the lowest was in Kayin State (24.44%). The BLV infection rate among the cattle of Myanmar reported in this study was higher than that reported in some Asian countries: Philippines (4.8–9.7%) [28] and India (27.9%) [36]. However, the BLV infection rate in this study was close to or lower than that in Korea (42.16%) [29], Thailand (58.74%) [32], and China (49.11%) [46].

BLV proviral load is a major diagnosis index for determining the virus transmission risk. To the best of our knowledge, this is the first study to estimate the proviral load among Myanmar cattle. We observed that the levels of BLV proviral load vary by geographical region within Myanmar. BLV-positive cattle in the three test areas (Sagaing, Yangon, and Kayin) exhibited high mean proviral load, while 3 other test areas (Mandalay, Magway and Bago) exhibited low mean proviral load. Furthermore, we demonstrated that the level of mean proviral load in the BLV-positive dairy farms correlated with the BLV infection rate in the dairy farms. Yangon region, which exhibited the highest mean proviral load also exhibited the highest prevalence of BLV infection among the 6 test areas. Additionally, dairy cattle in the farm YG7 located in Hmawbi township of Yangon region exhibited the highest mean proviral load and the highest infection rate (88.89%) among the 42 farms in the region. Our previous report [43] demonstrated that the cows with a proviral load of over 14,000 copies/10⁵ cells secreted BLV provirus into the nasal secretion. It has been suggested that the cows with the proviral load over 10,000/10⁵ cells have a higher risk of BLV transmission, while those with the proviral load under 10,000/10⁵ cells have a lower risk of BLV transmission. In our study we found that the three regions, which had high mean proviral load, were also associated with the high risk for BLV transmission (6, 6, and 4 for heads for Sagaing, Yangon and Kayin, respectively) as compared with that in the three regions which showed low mean proviral load (0, 0, and 1 heads for Mandalay, Magway and Bago, respectively) (Fig 2A). Similarly, the presence of cows with high risk for BLV transmission among the hybrids (Holstein Friesian × Myanmar native) (9 heads) was higher compared to that among the Myanmar native herds (8 heads) (Fig 2B). Thus, our observation indicates that depending on the geographical region and the population of cows with high risk for BLV transmission (proviral load over 10,000/10⁵ cells) may increase the BLV infection rate.

The climate of Yangon region and Kayin state, situated in southern Myanmar, is warm and is associated with high annual rainfall and humidity. Two sampling townships (Tamu and Kale) in the northern part of Sagaing region are geographically hilly and receive high annual rainfall. High rainfall cause floods in these area with long lasting wet condition. The climatic conditions may be one possible factor for the difference in BLV infection prevalence. Warm, wet, and humid climatic conditions in Sagaing, Yangon and Kayin may favor the blood sucking flies to transmit BLV from the infected animal to the healthy animal in the cattle farms, which result in higher proviral load and infection rate. Blood-sucking flies are regarded as one of the important mediators of BLV transmission [9, 10]. Interestingly, Mandalay and Magway regions, which are located in the central dry zone of Myanmar, and the sampling area in northern part of Bago region, which is closer to the central dry zone, have low annual rainfall and very high temperature (above 40°C in summer). Hence, these regions may be inhospitable for the survival of blood sucking insects.

The phylogenetic analysis based on partial and complete env-gp51 gene sequence in our previous study indicated that the isolated Myanmar strains were of genotype-10 [33]. The limited number of samples tested in the previous study, which was collected only from a small area of Myanmar, may be the reason for isolating only a single genotype. Hence, we hypothesized that if the samples are collected from different regions of the country, more genotypes of BLV may be identified. Indeed, three BLV genotypes (genotype-1, genotype-6 and genotype-10) were detected in our study based on the phylogenetic analysis of the partial env-gp51 nucleotide sequence of the BLV strains isolated from cattle reared in different regions of the country. Genotype-1 was the most prevalent genotype and was detected in the six test areas, while genotype-10 was detected in five test areas and genotype-6 was detected in only three test areas of Myanmar, as shown in Fig 1. Interestingly, the presence of two genotypes in the same herd was detected in some cattle farms of Sagaing, Bago, Yangon regions and Kayin state.

It is interesting to determine the origin of 3 different BLV genotypes discovered in this study. Myanmar has a history of importing exotic cattle breeds. Holstein Frisian and Jersey from Australia and New Zealand were imported in the 1970s for the artificial insemination project [44] and BLV genotype-1 seems to be existing for many years ago together with foreign cattle importation. The BLV genotype-6 has been detected in Brazil [23], Argentina [54], Peru, Paraguay and, Bolivia [31] in South America and Philippines [28], Thailand [32], India [36], and China [46] in Asia. Shindi and Thari breeds were imported to Myanmar many decades ago for the development of local draught cattle [44], but these cattle breeds came from Pakistan and not from India. The recently reported BLV infected regions in India [36] and China [46] are very far from the southern part of Myanmar where BLV genotype-6 was detected. Although the Philippines is situated in Southeast Asia, it is also very far from Myanmar and there is no cattle trade between the two countries. Based on these facts, it is unlikely that the origin of Myanmar BLV genotype-6 is India, China, or the Philippines. BLV genotype-1, genotype-6, and genotype-10 have also been reported in Thailand [32]. Furthermore, our phylogenetic analysis also indicated that the Myanmar BLV genotype-6 subtypes, G6a and G6b strains, exhibited homology to Thailand genotype-6 strains but not to Brazilian or Indian genotype-6 isolates. Myanmar genotype-10 strains isolated in our previous study and in this study were closely related to the genotype-10 strains of Thailand (Figs 3 and 4). Although there is no official information on live cattle trade between Myanmar and Thailand, there is a report by World Food and Agriculture Organization (FAO) that Thailand exported frozen semen of Brahman cattle and dairy cattle to its neighboring countries including Myanmar [55]. Recently, dairy and beef cattle have been imported into Myanmar through border trade. Therefore, we could speculate that the BLV strains, particularly genotype-6 and genotype-10, may have been transmitted to Myanmar through live cattle movement or cattle imported from Thailand.

Several studies have reported the existence of subgroups within BLV genotype-6, based on partial and/or complete env-gp51 sequence [28, 31, 32, 36]. We also reported an important finding of additional subgenotypes within BLV genotype-6 in this study. The topology of the phylogenetic tree in Fig 4 clearly elucidates the clustering of genotype-6 and its 3 subgroups. This finding was supported not only by the phylogenetic analysis of env-gp51 partial gene sequence with significant branch support (bootstrap values) and pairwise distance values, but also by the unique nucleotide substitutions. Interestingly, most of the Myanmar BLV strains in G6 clustered as subgroup 6c separately. The distinct features of the G6c strain were unique silent nucleotide substitutions at seven positions: 378, 471, 567, 585, 663, 708, and 750, and an amino acid substitution namely proline to serine at residue 177 within a portion of CD8+T-cell epitope and the E epitope, which were not observed in the other G6 strains isolated elsewhere (Fig 5A and 5B).

Most of the nucleotide and amino acid substitutions of Myanmar BLV genotype-6 and genotype-10 strains were located within the epitopes, including CD4+ T-cell epitope, the first ND, the viral G epitope, the second ND, the CD8+ T-cell epitope, the E epitope, the third ND and the B epitope (Fig 5A and 5B). These results were consistent with previous reports on the BLV genotype-6 and genotype-10 strains isolated elsewhere [28, 31, 32, 36, 54].

The present study indicated that the BLV prevalence (37.04%) is relatively higher and three BLV genotypes (genotype-1, genotype-6 and genotype-10) are commonly detected in Myanmar cattle than those observed in our previous study (9.1% prevalence and only genotype-10). BLV may have been transmitted to Myanmar through live cattle movement or cattle import, as these 3 BLV genotypes have also been detected in the neighboring countries. The introduction of infected cattle such as dairy breeds is considered as a major driver for BLV infection. To establish an appropriate cattle management policy for the control of BLV infection in Myanmar, it is necessary to issue a regulation to eliminate the cattle with high risk for BLV transmission, which exhibit a proviral load of over 10,000/10⁵ cells, and to allow rearing of only the cattle with low risk for BLV transmission, which exhibit a proviral load of under 10,000/10⁵ cells. Hence, there is a need to identify and characterize BLV strains among the cattle in the 7 out of 8 remaining regions/states of Myanmar which have boundaries with the neighboring countries to prevent BLV transmission into the country.

Supporting information

S1 Table. Evolutionary divergence among the nucleotide sequence of 66 Myanmar bovine leukemia virus (BLV) strains isolated in this study.

https://doi.org/10.1371/journal.pone.0229126.s001

(XLSX)

S2 Table. Evolutionary divergence between the nucleotide sequence of 66 Myanmar bovine leukemia virus (BLV) strains and 40 known strains of 10 genotypes used for the phylogenetic analysis.

https://doi.org/10.1371/journal.pone.0229126.s002

(XLS)

S3 Table. Evolutionary divergence of nucleotide sequence between the 23 Myanmar bovine leukemia virus (BLV) genotype-1 strains and 10 known strains of genotype-1 used for the phylogenetic analysis.

https://doi.org/10.1371/journal.pone.0229126.s003

(XLS)

S4 Table. Evolutionary divergence of nucleotide sequence between the 16 Myanmar bovine leukemia virus (BLV) genotype-6 strains and 7 known strains of genotype-6 used in the phylogenetic analysis.

https://doi.org/10.1371/journal.pone.0229126.s004

(XLS)

S5 Table. Evolutionary divergence of nucleotide sequence between the 26 Myanmar bovine leukemia virus (BLV) genotype-10 strains and 6 known strains of genotype-10 used in the phylogenetic analysis.

https://doi.org/10.1371/journal.pone.0229126.s005

(XLS)

S6 Table. Inter-subgenotype pairwise evolutionary divergence among the 3 subgenotypes of Myanmar genotype-6 strain.

https://doi.org/10.1371/journal.pone.0229126.s006

(XLS)

Acknowledgments

We thank LBVD staff, UVS staff, local veterinarians and farm owners for their kind assistance in blood sampling at the cattle farms in Myanmar, and the Support Unit at the Biomaterial Analysis, RIKEN BSI Research Resource Center for help with the sequence analysis. We also thank all members of the Viral Infectious Diseases Field, RIKEN, for technical assistance, help, and suggestions.

References

1. 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. pmid:24265629
- View Article
- PubMed/NCBI
- Google Scholar
2. Polat M, Takeshima SN, Aida Y. Epidemiology and genetic diversity of bovine leukemia virus. Virol J. 2017; 14: 209. pmid:29096657
- View Article
- PubMed/NCBI
- Google Scholar
3. Buehring GC, Shen HM, Jensen HM, Choi KY, Sun DJ, Nuovo G. Bovine leukemia virus DNA in human breast tissue. Emerg Infect Dis. 2014; 20: 772–782. pmid:24750974
- View Article
- PubMed/NCBI
- Google Scholar
4. OIE (World Organization for Animal Health). Enzootic bovine leucosis. In: Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. 2018; available on http://www.oie.int/fileadmin/Home/eng/Heath_standards/tahm/2.04.10_EBL.pdf
5. 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–6541. pmid:8341665
- View Article
- PubMed/NCBI
- Google Scholar
6. 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–2052. pmid:26723124
- View Article
- PubMed/NCBI
- Google Scholar
7. Sandev N, Koleva M, Binev R, Ilieva D. Influence of enzootic bovine leucosis virus upon the incidence of subclinical mastitis in cows at a different stage of infection. Vet Arh. 2004; 74: 411–4166.
- View Article
- Google Scholar
8. 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–1597. pmid:23332856
- View Article
- PubMed/NCBI
- Google Scholar
9. Kobayashi S, Tsutsui T, Yamamoto T, Hayama Y, Kameyama K, Konishi M, et al. Risk factors associated with within-herd transmission of bovine leukemia virus on dairy farms in Japan. BMC Vet Res. 2010; 6: 1. pmid:20055982
- View Article
- PubMed/NCBI
- Google Scholar
10. Ooshiro M, Konnai S, Katagiri Y, Afuso M, Arakaki N, Tsuha O, et al. Horizontal transmission of bovine leukemia virus from lymphocytotic cattle, and beneficial effects of insect vector control. Vet Rec. 2013; 173: 527.
- View Article
- Google Scholar
11. Mousavi S, Haghparast A, Mohammadi G, Tabatabaeizadeh SE. Prevalence of bovine leukemia virus (BLV) infection in the northeast of Iran. Veterinary research forum: an international quarterly journal. 2014; 5: 135–139.
- View Article
- Google Scholar
12. Gillet N, Florins A, Boxus M, Burteau C, Nigro A, Vandermeers F, et al. Mechanisms of leukemogenesis induced by bovine leukemia virus: prospects for novel antiretroviral therapies in human. Retrovirology. 2007; 4:18. pmid:17362524
- View Article
- PubMed/NCBI
- Google Scholar
13. Mamoun R, Astier T, Guillemain B, Duplan J. Bovine lymphosarcoma: expression of BLV-related proteins in cultured cells. J Gen Virol. 1983; 64: 1895–1905. pmid:6310032
- View Article
- PubMed/NCBI
- Google Scholar
14. Zarkik S, Decroly E, Wattiez R, Seidah NG, Burny A, Ruysschaert JM. Comparative processing of bovine leukemia virus envelope glycoprotein gp72 by subtilisin/kexin-like mammalian convertases. FEBS Lett. 1997; 406: 205–210. pmid:9109419
- View Article
- PubMed/NCBI
- Google Scholar
15. Bai L, Otsuki H, Sato H, Kohara J, Isogai E, Takeshima S-n, et al. Identification and characterization of common B cell epitope in bovine leukemia virus via high-throughput peptide screening system in infected cattle. Retrovirology. 2015; 12:106. pmid:26715158
- View Article
- PubMed/NCBI
- Google Scholar
16. Callebaut I, Voneche V, Mager A, Fumiere O, Krchnak V, Merza M, et al. Mapping of B-neutralizing and T-helper cell epitopes on the bovine leukemia virus external glycoprotein gp51. J Virol. 1993; 67: 5321–5327. pmid:7688821
- View Article
- PubMed/NCBI
- Google Scholar
17. Johnston ER, Radke K. The SU and TM envelope protein subunits of bovine leukemia virus are linked by disulfide bonds, both in cells and in virions. J Virol. 2000; 74: 2930–2935. pmid:10684314
- View Article
- PubMed/NCBI
- Google Scholar
18. 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: 4180–4188. pmid:1696635
- View Article
- PubMed/NCBI
- Google Scholar
19. 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. J Vet Med B Infect Dis Vet Public Health. 2002; 49: 325–331. pmid:12420867
- View Article
- PubMed/NCBI
- Google Scholar
20. 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. J Gen Virol. 1990; 71: 1737–1746. pmid:2167927
- View Article
- PubMed/NCBI
- Google Scholar
21. Beier D, Blankenstein P, Marquardt O, Kuzmak J. Identification of different BLV provirus isolates by PCR, RFLPA and DNA sequencing. Berl Munch Tierarztl Wochenschr. 2001; 114: 252–256. pmid:11505797
- View Article
- PubMed/NCBI
- Google Scholar
22. 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: 39–47. pmid:15885933
- View Article
- PubMed/NCBI
- Google Scholar
23. 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: 343–350. pmid:16917740
- View Article
- PubMed/NCBI
- Google Scholar
24. Monti G, Schrijver R, Beier D. Genetic diversity and spread of bovine leukaemia virus isolates in Argentine dairy cattle. Arch Virol. 2005 150: 443–458. pmid:15592892
- View Article
- PubMed/NCBI
- Google Scholar
25. 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: 343–348. pmid:21095208
- View Article
- PubMed/NCBI
- Google Scholar
26. 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:1281–1290. pmid:22488472
- View Article
- PubMed/NCBI
- Google Scholar
27. 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: e58705. pmid:23527009
- View Article
- PubMed/NCBI
- Google Scholar
28. Polat M, Ohno A, Takeshima SN, Kim J, Kikuya M, Matsumoto Y, et al. Detection and molecular characterization of bovine leukemia virus in Philippine cattle. Arch Virol. 2015; 160: 285–296. pmid:25399399
- View Article
- PubMed/NCBI
- Google Scholar
29. Lee E, Kim EJ, Joung HK, Kim BH, Song JY, Cho IS, et al. Sequencing and phylogenetic analysis of the gp51 gene from Korean bovine leukemia virus isolates. Virol J. 2015; 12: 64. pmid:25879943
- View Article
- PubMed/NCBI
- Google Scholar
30. Ochirkhuu N, Konnai S, Odbileg R, Nishimori A, Okagawa T, Murata S, et al. Detection of bovine leukemia virus and identification of its genotype in Mongolian cattle. Arch Virol. 2016; 161: 985–991. pmid:26711456
- View Article
- PubMed/NCBI
- Google Scholar
31. Polat M, Takeshima SN, Hosomichi K, Kim J, Miyasaka T, Yamada K, et al. (2016) A new genotype of bovine leukemia virus in South America identified by NGS-based whole genome sequencing and molecular evolutionary genetic analysis. Retrovirology 13: 4. pmid:26754835
- View Article
- PubMed/NCBI
- Google Scholar
32. Lee E, Kim EJ, Ratthanophart J, Vitoonpong R, Kim BH, Cho IS, et al. Molecular epidemiological and serological studies of bovine leukemia virus (BLV) infection in Thailand cattle. Infect Genet Evol. 2016; 41: 245–254. pmid:27090024
- View Article
- PubMed/NCBI
- Google Scholar
33. 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–437. pmid:27771791
- View Article
- PubMed/NCBI
- Google Scholar
34. Pandey GS, Simulundu E, Mwiinga D, Samui KL, Mweene AS, Kajihara M, et al. Clinical and subclinical bovine leukemia virus infection in a dairy cattle herd in Zambia. Arch Virol. 2017; 162: 1051–1056. pmid:28025710
- View Article
- PubMed/NCBI
- Google Scholar
35. Pluta A, Rola-Łuszczak M, Kubis P, Balov S, Moskalik R, Choudhury B, et al. Molecular characterization of bovine leukemia virus from Moldovan dairy cattle. Arch Virol. 2017; 162: 1563–1576. pmid:28213870
- View Article
- PubMed/NCBI
- Google Scholar
36. Gautam S, Mishra N, Kalaiyarasu S, Jhade SK, Sood R. Molecular characterization of bovine leukaemia virus (BLV) strains reveals existence of genotype 6 in cattle in India with evidence of a new subgenotype. Transbound Emerg Dis. 2018; 1–11.
- View Article
- Google Scholar
37. 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. Veterinarski Arhiv. 2004; 74: 411–416.
- View Article
- Google Scholar
38. 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–1597. pmid:23332856
- View Article
- PubMed/NCBI
- Google Scholar
39. 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. pmid:21044304
- View Article
- PubMed/NCBI
- Google Scholar
40. 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–290. pmid:26321160
- View Article
- PubMed/NCBI
- Google Scholar
41. Panei CJ, Takeshima SN, Omori T, Nunoya T, Davis WC, Ishizaki H, et al. 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. pmid:23641811
- View Article
- PubMed/NCBI
- Google Scholar
42. Sato H, Watanuki S, Murakami H, Sato R, Ishizaki H, Aida Y. Development of a luminescence syncytium induction assay (LuSIA) for easily detecting and quantitatively measuring bovine leukemia virus infection. Arch Virol. 2018; 163: 1519–1530. pmid:29455325
- View Article
- PubMed/NCBI
- Google Scholar
43. Yuan Y, Kitamura-Muramatsu Y, Saito S, Ishizaki H, Nakano M, Haga S, et al. 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–254. pmid:26298004
- View Article
- PubMed/NCBI
- Google Scholar
44. National Consultative Committee, Myanmar. National report on animal genetic resources, The Union of Myanmar. 2013.
45. Wang M, Wang Y, Baloch AR, Pan Y, Xu F, Tian L, et al. Molecular epidemiology and characterization of bovine leukemia virus in domestic yaks (Bos grunniens) on the Qinghai‑Tibet Plateau, China. Arch Virol. 2018; 163: 659–670. pmid:29224130
- View Article
- PubMed/NCBI
- Google Scholar
46. Yang Y, Chu S,.Shang S, Wang C. Genotyping and single nucleotide polymorphism analysis of bovine leukemia virus in Chinese dairy cattle. J Dairy Sci. 2019; 102: 3469–3473. pmid:30712932
- View Article
- PubMed/NCBI
- Google Scholar
47. 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–1332. pmid:25731158
- View Article
- PubMed/NCBI
- Google Scholar
48. 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: 493–505. pmid:15592894
- View Article
- PubMed/NCBI
- Google Scholar
49. 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: 481–489. pmid:20169372
- View Article
- PubMed/NCBI
- Google Scholar
50. 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 USA. 1985; 82: 677–681. pmid:2983308
- View Article
- PubMed/NCBI
- Google Scholar
51. Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular Biology and Evolution. 2016; 33: 1870–1874. pmid:27004904
- View Article
- PubMed/NCBI
- Google Scholar
52. Jimba M, Takeshima SN, Murakami H, Kohara J, Kobayashi N, et al. (2012) BLV-CoCoMo-qPCR: a useful tool for evaluating bovine leukemia virus infection status. BMC Vet Res 8: 167. pmid:22995575
- View Article
- PubMed/NCBI
- Google Scholar
53. Zhao X, Buehring GC. Natural genetic variations in bovine leukemia virus envelope gene: possible effects of selection and escape. Virology. 2007; 366: 150–165. pmid:17498765
- View Article
- PubMed/NCBI
- Google Scholar
54. 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 Gen Virol. 2009; 90: 2788–2797. pmid:19587134
- View Article
- PubMed/NCBI
- Google Scholar
55. FAO (World Food and Agriculture Organization). The second report on the state of the world’s animal genetic resources for food and agriculture, edited by Scherf BD & Pilling D. FAO Commission on Genetic Resources for Food and Agriculture Assessments. Rome. 2015. Available from: http://www.fao.org/3/a-i4787e/index.html

[ref1] 1. 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. pmid:24265629
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Polat M, Takeshima SN, Aida Y. Epidemiology and genetic diversity of bovine leukemia virus. Virol J. 2017; 14: 209. pmid:29096657
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Buehring GC, Shen HM, Jensen HM, Choi KY, Sun DJ, Nuovo G. Bovine leukemia virus DNA in human breast tissue. Emerg Infect Dis. 2014; 20: 772–782. pmid:24750974
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. OIE (World Organization for Animal Health). Enzootic bovine leucosis. In: Manual of Diagnostic Tests and Vaccines for Terrestrial Animals. 2018; available on http://www.oie.int/fileadmin/Home/eng/Heath_standards/tahm/2.04.10_EBL.pdf

[ref5] 5. 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–6541. pmid:8341665
View Article
PubMed/NCBI
Google Scholar

[15] View Article

[16] PubMed/NCBI

[17] Google Scholar

[ref6] 6. 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–2052. pmid:26723124
View Article
PubMed/NCBI
Google Scholar

[19] View Article

[20] PubMed/NCBI

[21] Google Scholar

[ref7] 7. Sandev N, Koleva M, Binev R, Ilieva D. Influence of enzootic bovine leucosis virus upon the incidence of subclinical mastitis in cows at a different stage of infection. Vet Arh. 2004; 74: 411–4166.
View Article
Google Scholar

[23] View Article

[24] Google Scholar

[ref8] 8. 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–1597. pmid:23332856
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref9] 9. Kobayashi S, Tsutsui T, Yamamoto T, Hayama Y, Kameyama K, Konishi M, et al. Risk factors associated with within-herd transmission of bovine leukemia virus on dairy farms in Japan. BMC Vet Res. 2010; 6: 1. pmid:20055982
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref10] 10. Ooshiro M, Konnai S, Katagiri Y, Afuso M, Arakaki N, Tsuha O, et al. Horizontal transmission of bovine leukemia virus from lymphocytotic cattle, and beneficial effects of insect vector control. Vet Rec. 2013; 173: 527.
View Article
Google Scholar

[34] View Article

[35] Google Scholar

[ref11] 11. Mousavi S, Haghparast A, Mohammadi G, Tabatabaeizadeh SE. Prevalence of bovine leukemia virus (BLV) infection in the northeast of Iran. Veterinary research forum: an international quarterly journal. 2014; 5: 135–139.
View Article
Google Scholar

[37] View Article

[38] Google Scholar

[ref12] 12. Gillet N, Florins A, Boxus M, Burteau C, Nigro A, Vandermeers F, et al. Mechanisms of leukemogenesis induced by bovine leukemia virus: prospects for novel antiretroviral therapies in human. Retrovirology. 2007; 4:18. pmid:17362524
View Article
PubMed/NCBI
Google Scholar

[40] View Article

[41] PubMed/NCBI

[42] Google Scholar

[ref13] 13. Mamoun R, Astier T, Guillemain B, Duplan J. Bovine lymphosarcoma: expression of BLV-related proteins in cultured cells. J Gen Virol. 1983; 64: 1895–1905. pmid:6310032
View Article
PubMed/NCBI
Google Scholar

[44] View Article

[45] PubMed/NCBI

[46] Google Scholar

[ref14] 14. Zarkik S, Decroly E, Wattiez R, Seidah NG, Burny A, Ruysschaert JM. Comparative processing of bovine leukemia virus envelope glycoprotein gp72 by subtilisin/kexin-like mammalian convertases. FEBS Lett. 1997; 406: 205–210. pmid:9109419
View Article
PubMed/NCBI
Google Scholar

[48] View Article

[49] PubMed/NCBI

[50] Google Scholar

[ref15] 15. Bai L, Otsuki H, Sato H, Kohara J, Isogai E, Takeshima S-n, et al. Identification and characterization of common B cell epitope in bovine leukemia virus via high-throughput peptide screening system in infected cattle. Retrovirology. 2015; 12:106. pmid:26715158
View Article
PubMed/NCBI
Google Scholar

[52] View Article

[53] PubMed/NCBI

[54] Google Scholar

[ref16] 16. Callebaut I, Voneche V, Mager A, Fumiere O, Krchnak V, Merza M, et al. Mapping of B-neutralizing and T-helper cell epitopes on the bovine leukemia virus external glycoprotein gp51. J Virol. 1993; 67: 5321–5327. pmid:7688821
View Article
PubMed/NCBI
Google Scholar

[56] View Article

[57] PubMed/NCBI

[58] Google Scholar

[ref17] 17. Johnston ER, Radke K. The SU and TM envelope protein subunits of bovine leukemia virus are linked by disulfide bonds, both in cells and in virions. J Virol. 2000; 74: 2930–2935. pmid:10684314
View Article
PubMed/NCBI
Google Scholar

[60] View Article

[61] PubMed/NCBI

[62] Google Scholar

[ref18] 18. 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: 4180–4188. pmid:1696635
View Article
PubMed/NCBI
Google Scholar

[64] View Article

[65] PubMed/NCBI

[66] Google Scholar

[ref19] 19. 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. J Vet Med B Infect Dis Vet Public Health. 2002; 49: 325–331. pmid:12420867
View Article
PubMed/NCBI
Google Scholar

[68] View Article

[69] PubMed/NCBI

[70] Google Scholar

[ref20] 20. 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. J Gen Virol. 1990; 71: 1737–1746. pmid:2167927
View Article
PubMed/NCBI
Google Scholar

[72] View Article

[73] PubMed/NCBI

[74] Google Scholar

[ref21] 21. Beier D, Blankenstein P, Marquardt O, Kuzmak J. Identification of different BLV provirus isolates by PCR, RFLPA and DNA sequencing. Berl Munch Tierarztl Wochenschr. 2001; 114: 252–256. pmid:11505797
View Article
PubMed/NCBI
Google Scholar

[76] View Article

[77] PubMed/NCBI

[78] Google Scholar

[ref22] 22. 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: 39–47. pmid:15885933
View Article
PubMed/NCBI
Google Scholar

[80] View Article

[81] PubMed/NCBI

[82] Google Scholar

[ref23] 23. 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: 343–350. pmid:16917740
View Article
PubMed/NCBI
Google Scholar

[84] View Article

[85] PubMed/NCBI

[86] Google Scholar

[ref24] 24. Monti G, Schrijver R, Beier D. Genetic diversity and spread of bovine leukaemia virus isolates in Argentine dairy cattle. Arch Virol. 2005 150: 443–458. pmid:15592892
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref25] 25. 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: 343–348. pmid:21095208
View Article
PubMed/NCBI
Google Scholar

[92] View Article

[93] PubMed/NCBI

[94] Google Scholar

[ref26] 26. 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:1281–1290. pmid:22488472
View Article
PubMed/NCBI
Google Scholar

[96] View Article

[97] PubMed/NCBI

[98] Google Scholar

[ref27] 27. 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: e58705. pmid:23527009
View Article
PubMed/NCBI
Google Scholar

[100] View Article

[101] PubMed/NCBI

[102] Google Scholar

[ref28] 28. Polat M, Ohno A, Takeshima SN, Kim J, Kikuya M, Matsumoto Y, et al. Detection and molecular characterization of bovine leukemia virus in Philippine cattle. Arch Virol. 2015; 160: 285–296. pmid:25399399
View Article
PubMed/NCBI
Google Scholar

[104] View Article

[105] PubMed/NCBI

[106] Google Scholar

[ref29] 29. Lee E, Kim EJ, Joung HK, Kim BH, Song JY, Cho IS, et al. Sequencing and phylogenetic analysis of the gp51 gene from Korean bovine leukemia virus isolates. Virol J. 2015; 12: 64. pmid:25879943
View Article
PubMed/NCBI
Google Scholar

[108] View Article

[109] PubMed/NCBI

[110] Google Scholar

[ref30] 30. Ochirkhuu N, Konnai S, Odbileg R, Nishimori A, Okagawa T, Murata S, et al. Detection of bovine leukemia virus and identification of its genotype in Mongolian cattle. Arch Virol. 2016; 161: 985–991. pmid:26711456
View Article
PubMed/NCBI
Google Scholar

[112] View Article

[113] PubMed/NCBI

[114] Google Scholar

[ref31] 31. Polat M, Takeshima SN, Hosomichi K, Kim J, Miyasaka T, Yamada K, et al. (2016) A new genotype of bovine leukemia virus in South America identified by NGS-based whole genome sequencing and molecular evolutionary genetic analysis. Retrovirology 13: 4. pmid:26754835
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref32] 32. Lee E, Kim EJ, Ratthanophart J, Vitoonpong R, Kim BH, Cho IS, et al. Molecular epidemiological and serological studies of bovine leukemia virus (BLV) infection in Thailand cattle. Infect Genet Evol. 2016; 41: 245–254. pmid:27090024
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref33] 33. 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–437. pmid:27771791
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref34] 34. Pandey GS, Simulundu E, Mwiinga D, Samui KL, Mweene AS, Kajihara M, et al. Clinical and subclinical bovine leukemia virus infection in a dairy cattle herd in Zambia. Arch Virol. 2017; 162: 1051–1056. pmid:28025710
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref35] 35. Pluta A, Rola-Łuszczak M, Kubis P, Balov S, Moskalik R, Choudhury B, et al. Molecular characterization of bovine leukemia virus from Moldovan dairy cattle. Arch Virol. 2017; 162: 1563–1576. pmid:28213870
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref36] 36. Gautam S, Mishra N, Kalaiyarasu S, Jhade SK, Sood R. Molecular characterization of bovine leukaemia virus (BLV) strains reveals existence of genotype 6 in cattle in India with evidence of a new subgenotype. Transbound Emerg Dis. 2018; 1–11.
View Article
Google Scholar

[136] View Article

[137] Google Scholar

[ref37] 37. 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. Veterinarski Arhiv. 2004; 74: 411–416.
View Article
Google Scholar

[139] View Article

[140] Google Scholar

[ref38] 38. 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–1597. pmid:23332856
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref39] 39. 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. pmid:21044304
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref40] 40. 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–290. pmid:26321160
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref41] 41. Panei CJ, Takeshima SN, Omori T, Nunoya T, Davis WC, Ishizaki H, et al. 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. pmid:23641811
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref42] 42. Sato H, Watanuki S, Murakami H, Sato R, Ishizaki H, Aida Y. Development of a luminescence syncytium induction assay (LuSIA) for easily detecting and quantitatively measuring bovine leukemia virus infection. Arch Virol. 2018; 163: 1519–1530. pmid:29455325
View Article
PubMed/NCBI
Google Scholar

[158] View Article

[159] PubMed/NCBI

[160] Google Scholar

[ref43] 43. Yuan Y, Kitamura-Muramatsu Y, Saito S, Ishizaki H, Nakano M, Haga S, et al. 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–254. pmid:26298004
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref44] 44. National Consultative Committee, Myanmar. National report on animal genetic resources, The Union of Myanmar. 2013.

[ref45] 45. Wang M, Wang Y, Baloch AR, Pan Y, Xu F, Tian L, et al. Molecular epidemiology and characterization of bovine leukemia virus in domestic yaks (Bos grunniens) on the Qinghai‑Tibet Plateau, China. Arch Virol. 2018; 163: 659–670. pmid:29224130
View Article
PubMed/NCBI
Google Scholar

[167] View Article

[168] PubMed/NCBI

[169] Google Scholar

[ref46] 46. Yang Y, Chu S,.Shang S, Wang C. Genotyping and single nucleotide polymorphism analysis of bovine leukemia virus in Chinese dairy cattle. J Dairy Sci. 2019; 102: 3469–3473. pmid:30712932
View Article
PubMed/NCBI
Google Scholar

[171] View Article

[172] PubMed/NCBI

[173] Google Scholar

[ref47] 47. 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–1332. pmid:25731158
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

[ref48] 48. 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: 493–505. pmid:15592894
View Article
PubMed/NCBI
Google Scholar

[179] View Article

[180] PubMed/NCBI

[181] Google Scholar

[ref49] 49. 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: 481–489. pmid:20169372
View Article
PubMed/NCBI
Google Scholar

[183] View Article

[184] PubMed/NCBI

[185] Google Scholar

[ref50] 50. 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 USA. 1985; 82: 677–681. pmid:2983308
View Article
PubMed/NCBI
Google Scholar

[187] View Article

[188] PubMed/NCBI

[189] Google Scholar

[ref51] 51. Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular Biology and Evolution. 2016; 33: 1870–1874. pmid:27004904
View Article
PubMed/NCBI
Google Scholar

[191] View Article

[192] PubMed/NCBI

[193] Google Scholar

[ref52] 52. Jimba M, Takeshima SN, Murakami H, Kohara J, Kobayashi N, et al. (2012) BLV-CoCoMo-qPCR: a useful tool for evaluating bovine leukemia virus infection status. BMC Vet Res 8: 167. pmid:22995575
View Article
PubMed/NCBI
Google Scholar

[195] View Article

[196] PubMed/NCBI

[197] Google Scholar

[ref53] 53. Zhao X, Buehring GC. Natural genetic variations in bovine leukemia virus envelope gene: possible effects of selection and escape. Virology. 2007; 366: 150–165. pmid:17498765
View Article
PubMed/NCBI
Google Scholar

[199] View Article

[200] PubMed/NCBI

[201] Google Scholar

[ref54] 54. 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 Gen Virol. 2009; 90: 2788–2797. pmid:19587134
View Article
PubMed/NCBI
Google Scholar

[203] View Article

[204] PubMed/NCBI

[205] Google Scholar

[ref55] 55. FAO (World Food and Agriculture Organization). The second report on the state of the world’s animal genetic resources for food and agriculture, edited by Scherf BD & Pilling D. FAO Commission on Genetic Resources for Food and Agriculture Assessments. Rome. 2015. Available from: http://www.fao.org/3/a-i4787e/index.html

Figures

Abstract

Introduction

Materials and methods

Animals and extraction of genomic DNA

Evaluation of BLV provirus by BLV-CoCoMo-qPCR-2.

PCR amplification of BLV env-gp51 gene fragments and nucleotide sequences

Construction of the phylogenetic tree

Statistical analysis

Animal handling and research ethics

Results

Infection rate prevalence of BLV

Estimation of BLV proviral load

Phylogenetic analysis of the partial sequence of BLV env-gp51

BLV genotype distribution in Myanmar

Nucleotide and amino acid substitution in the BLV env gp51 gene of BLV strains isolated in Myanmar

Discussion

Supporting information

S1 Table. Evolutionary divergence among the nucleotide sequence of 66 Myanmar bovine leukemia virus (BLV) strains isolated in this study.

S2 Table. Evolutionary divergence between the nucleotide sequence of 66 Myanmar bovine leukemia virus (BLV) strains and 40 known strains of 10 genotypes used for the phylogenetic analysis.

S3 Table. Evolutionary divergence of nucleotide sequence between the 23 Myanmar bovine leukemia virus (BLV) genotype-1 strains and 10 known strains of genotype-1 used for the phylogenetic analysis.

S4 Table. Evolutionary divergence of nucleotide sequence between the 16 Myanmar bovine leukemia virus (BLV) genotype-6 strains and 7 known strains of genotype-6 used in the phylogenetic analysis.

S5 Table. Evolutionary divergence of nucleotide sequence between the 26 Myanmar bovine leukemia virus (BLV) genotype-10 strains and 6 known strains of genotype-10 used in the phylogenetic analysis.

S6 Table. Inter-subgenotype pairwise evolutionary divergence among the 3 subgenotypes of Myanmar genotype-6 strain.

Acknowledgments

References