CC BY
84
find
Foods and Raw Materials, 2018, vol. 6, no. 2
ISSN 2308-4057 (Print) ISSN 2310-9599 (Online)
Research Article DOI: http://doi .org/10.21603/2308-4057-2018-2-314-324
Open Access Available online at http:jfrm.ru
Genetic identification of bovine leukaemia virus
Irina M. Donnika , Ramil R. Vafinb'* , Aram G. Galstyanb , Anna S. Krivonogova0 , Aigul Y. Shaevad , Khamid Kh. Gilmanovd , Rufiya G. Karimovad , Sergey V. Tyulkine , and Jacek Kuzmakf
a Russian Academy of Sciences, Leninsky Ave. 14, Moscow 119991, Russian Federation
bAll-Russian Research Institute of Brewing, Non-Alcoholic and Wine Industry, Rossolimo Str. 7, Moscow 119021, Russian Federation
c Ural Federal Agrarian Research Centre of the Ural branch of the Russian Academy of Science, Belinskogo Str. 112A, Ekaterinburg 620142, Russian Federation
d N.E. Bauman Kazan State Academy of Veterinary Medicine, Sibirsky Tract Str. 35, Kazan 420029, Russian Federation
e Kazan State Agrarian University, K. Marx Str. 65, Kazan 420015, Russian Federation
f National Veterinary Research Institute, Partyzantow Ave. 57, Pulawy, Poland
* e-mail: [email protected]
Received August 27, 2018; Accepted in revised form October 10, 2018; Published December 20, 2018
Abstract: Molecular genetic research methods make it possible to evaluate the genetic diversity of bovine leukemia virus (BLV) and are the most informative approaches to its genetic identification. Molecular genetic research methods work well for the phylogenetic analysis of sequenced nucleotide DNA sequences of the provirus, as well as for the polymerase chain reaction-restriction fragment length polymorphism analysis (PCR-RFLP) according to the phylogenetic classification of the pathogen. The purpose of the research was to study the scientific and methodological approaches to the genetic identification of bovine leukemia virus, integrated into the molecular monitoring of infection of cattle with BLV genotypes. The authors used PCR-RFLP-genotyping and comparative phylogenetic analysis of aligned nucleotide sequences of the env gene fragment of the BLV provirus isolates to detect the genotypic affiliation of the cattle from twenty-one livestock farms of the Republic of Tatarstan. As a result, isolates of four out of ten BLV genotypes were found in the Tatarstani cattle, namely genotypes 1, 4, 7, and 8. The research involved a comparative analysis of 505 nucleotide sequences of a fragment of the BLV env gene, including those deposited in GenBank NCBI. The analysis confirms the inconsistency of several earlier PCR-RFLP typing strategies with the current approach in assessing the genotypic diversity by phylogenetic analysis. The improved strategy of PCR-RFLP genotyping of BLV corresponds with its modern phylogenetic classification. The strategy makes it possible to identify all the known genotypes of the viral pathogen. Its validity has been proved by in silico modelling of restrictogrammes and a phylogenetic analysis of the env gene fragment of 57 reference isolates of ten BLV genotypes that generate 57 genotype-associated combinations of diagnostically significant PCR-RFLP profiles.
Keywords: Bovine leukaemia virus, BLV, cattle, gene, genotype, genetic identification, PCR, RFLP, sequencing
Please cite this article in press as: Donnik I.M., Vafin R.R., Galstyan A.G., et al. Genetic identification of bovine leukaemia virus. Foods and Raw Materials, 2018, vol. 6, no. 2, pp. 314-324. DOI: http://doi.org/10.21603/2308-4057-2018-2-314-324.
INTRODUCTION
Enzootic Bovine Leukosis (EBL) is a chronic infectious disease of a tumorous nature. It causes significant economic damage to the dairy and beef cattle industry due to poor production, low quality, cattle
mortality, and expensive epidemic prevention measures [1, 2].
Foodstuffs of infected animals can be dangerous to humans due to the harmful metabolites it contains. The causative agent affects all kinds of raw material (milk,
Copyright © 2018, Donnik et al. This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material for any purpose, even commercially, provided the original work is properly cited and states its license.
meat, by-products) and products, remaining a potential source of human infection [3-5].
Pasteurization of milk inactivates the virus but does not degrade its genome. The genetic material of the provirus maintains its integrity in canned meat [6]. Moreover, there are dairy products with partial pasteurization regime, e.g. classic cheeses, granulated cottage cheese, powdered milk with a low heating temperature, etc. The temperature processing parameters used in accordance with the regulatory and technical documentation cannot destroy harmful metabolites and, in some cases, do not kill the virus [7].
According to some researches, DNA of BLV provirus was found in epithelial cells of human mammary glands, including those of breast cancer patients. The hypothesis states that BLV may destabilize the host genome, thus leading to cancerous degeneration of cells [8-12].
Obtaining high-quality raw materials of animal origin is the most important challenge for meat and dairy industry. The challenge includes the development of functional and gerodietic foods [13-15].
According to the requirements of the Technical Regulations of the Customs Union "On safety of milk and dairy products" (TR CU 033/2013), BLV preventive measures and eradication activities are extremely important, given the significant prevalence of this incurable disease in the Russian Federation [16, 17].
An early genetic diagnosing of the pathogen is part of the system of anti-epizootic measures, followed by the removal of infected animals from the herd. Molecular genetic research methods make it possible to assess the genetic diversity of BLV [18]. This is the most informative approach to the gene identification of the virus. Molecular genetic research methods work well for the phylogenetic analysis of nucleotide DNA sequences of the provirus, as well as for PCR-RFLP analysis according to the phylogenetic classification of the pathogen [19].
The current phylogenetic classification of BLV includes ten genotypes. The first seven genotypes were described by Argentinean scientists in 2009 [20], while genotype 8 was described by researchers from Russia [21-23], Croatia [24], and a European team of scientists [25] in 2011-2013. Genotype 9 was investigated by a team of Argentinean, Chilean, and Japanese scientists in 2016 [26]. Genotype 10 was described by a team of researchers from Thailand and South Korea [27] in 2016.
The objective of the current research was to study the scientific and methodological approaches to the genetic identification of BLV integrated into the molecular monitoring of infection of cattle herds with BLV genotypes. The following tasks were set:
- to establish the genotypes of BLV isolates in Tearstain cattle;
- to define the types of BLV isolates with deciphered nucleotide sequences of the env gene fragment, depending on the chosen gene identification strategy;
- to improve the strategy of PCR-RFLP-genotyping of BLV and make it consistent with the modern phylogenetic BLV classification.
STAGY OBYECTS AND METHODS
The research involved a total of 179 blood samples from AGID-positive cows. The samples were provided by agricultural enterprises from 21 districts of the Republic of Tatarstan. The samples were genetically examined for BLV. The examination included a phylogenetic analysis of sequenced env gene fragment of the pathogen and a PCR-RFLP-genotyping consistent with the phylogenetic classification of the infectious agent.
To extract DNA from the whole conserved blood, we used a commercial PCR diagnostic kit, 'DNA-Sorb B', produced by the Central Research Institute of Epidemiology of the Federal Supervisory Service for Consumer Rights and Human Welfare, Ministry of Health of the Russian Federation.
Nested PCR with extracted samples of BLV proviral DNA was performed with external (env5032 and env5608) and internal (env5099 and env5521) primers initiating the generation of env gene 444 bp fragment of causative agent at the final stage of reaction [28].
Restriction endonucleases used in PCR-RFLP-genotyping of BLV were consistent with its phylogenetic classification: BstYI (isoshizomer BstX2I), HphI (isoshizomer AsuHPI), Haelll, PvuII, SspI. The NEBcutter v.2.0 web resource was used for PCR-RFLP modelling.
For the detection of the obtained results of PCR and PCR-RFLP analysis, 2.5% agarose gel horizontal electrophoresis was applied with a TBE buffer (pH 8.0) containing ethidium bromide. The electrophoregrammes were examined in a UV-transilluminator (X = 310 nm). The sizes of the generated DNA fragments were compared with standard DNA molecular weight markers (SibEnzyme Ltd, Russia).
Sequencing of the PCR amplification products of the env gene fragment of detected BLV provirus isolates was performed on ABI PRISM 3100 Genetic Analyser (Applied Biosystems, USA) in the laboratory of Scientific and Technical Complex Sintol (Russia). Internal oligonucleotide primers env5099 and env5521 were used as sequencing. The sequenced fragments of the env gene of BLV provirus isolates were aligned with the corresponding nucleotide sequences of the reference BLV isolates from GenBank with the help of BLAST and MEGA-4 programmes. The last stage included a phylogenetic analysis.
RESULTS AND DISCUSSION
The study featured a PCR-RFLP-genotyping and a comparative phylogenetic analysis of the aligned sequences of the env gene fragment of BLV provirus isolates from 21 districts of the Republic of Tatarstan.
As a result, out of 179 identified isolates, ten isolates belonged to genotype 1; 106 isolates belonged to the cluster of genotype 4; 55 were characterized as genotype 7, and the remaining eight provirus isolates belonged to genotype 8 (Table 1).
According to the results obtained by PCR-RFLP-genotyping and phylogenetic analysis of sequenced env gene fragment, there are four out of ten currently known BLV genotypes in Tatarstan: 1, 4, 7, and 8.
Fig. 1 shows the genotypes of BLV isolated with the help of phylogenetic analysis of nucleotide sequences of env gene fragment.
An additional assessment of the heterogeneity of the reference BLV representatives for the env gene included an analysis of the intra- and intergenotypic heterogeneity of genotypes. The data in Table 2 indicate that it is impossible to use the 'heterogeneous' criterion for assessing the genetic diversity of BLV.
As part of the next task, BLV isolates with the decoded nucleotide sequences of the env gene fragment were identified according to the chosen genetic identification strategy. The degree of consistency of genotypic approaches was assessed by comparing the data of the in silico PCR-RFLP and the phylogenetic analyses.
A comparative analysis of 505 nucleotide sequences of the BLV env gene locus, including those deposited with GenBank NCBI, confirms the inconsistency of a number of earlier PCR-RFLP typing strategies [28-30] with the current approach in assessing the genotypic diversity by means of phylogenetic analysis.
Thus, the BLV isolates that belong to the Belgian subgroup according to D. Beier et al. (2001) [28], belong to genotype 4 according to the phylogenetic classification; Australian subgroup can be referred to genotypes 1, 3, 6, 8, or 9; Japanese subgroup - to genotypes 1, 6, or 7. In addition, the genotyping strategy [28] includes 11 additional unique combinations of PCR-RFLP profiles, conditionally identical to 11 unclassifiable BLV subgroups (Table 3).
Besides, BLV isolates that belong to genotype 1 and 7 according to the phylogenetic analysis, can be referred to Australian and Japanese subgroups, as well as to three unclassifiable subgroups, according to the strategy of D. Beier et al. (2001) [28]; genotypes 2, 5, and 10 belong to two unclassifiable subgroups;
genotypes 3 and 8 - to Australian subgroup; genotype 4 - to the Belgian subgroup and four unclassifiable subgroups; the genotype 5 - to two unclassifiable subgroups; genotype 6 - to the Australian, Japanese and two unclassifiable subgroups; genotype 9 - to the Australian and one unclassifiable subset (Table 3).
BLV isolates that were genotyped according to M. Licursi et al. (2002) [29] as genotype 1 may belong to genotypes 1, 4, 6, or 7 according to the phylogenetic classification; genotype 3 - to genotypes 1, 6, or 7; genotype 5 - to genotypes 1, 3, 6, 7, or 9; genotype 6 -to genotypes 2, 4, 5, or 7 (Table 4).
For the genotyping strategy described in [29], there are 19 unique combinations of PCR-RFLP profiles that are conditionally identical to 19 unclassifiable BLV genotypes (Table 4).
Besides, BLV isolates that are genotyped according to phylogenetic analysis as genotype 1 may refer to 1, 3, 5, and three unclassifiable BLV genotypes according to the strategy of M. Licursi et al. (2002) [29]; genotype 2 belongs to genotype 6 and two unclassifiable genotypes; the genotype 3 - to genotype 5 and one unclassifiable genotype; genotype 4 - to genotypes 1 and 6 and five unclassifiable genotypes; genotype 5 - to genotype 6 and two unclassifiable genotypes; genotype 6 - to genotypes 1, 3, and 5 and three unclassifiable genotypes; genotype 7 - to genotypes 1, 3, and 6 and four unclassifiable genotypes; genotype 8 - to one unclassifiable genotype; genotype 9 - to genotype 5; genotype 10 - to three unclassifiable genotypes (Table 4).
It should be mentioned that, when analyzing in silico PCR-RFLP data from 505 BLV representatives, we found not a single nucleotide sequence of the env gene fragment that would belong to genotypes 2 and 4 according to M. Licursi et al. (2002) (Table 4). This fact did not make it possible to prove the actual existence of PCR-RFLP profiles indicated for these two BLV genotypes.
Table 1. Distribution of 179 genotyped samples of BLV provirus DNA according to 21 districts of the Republic of Tatarstan, Russian Federation
Districts Number of analysed samples 12 34 BLV genotypes 5 6 7 8 9 10
1 Aznakaevsky 10 - - -9 -- 1 - --
2 Al'keyevsky 13 - - -5 -- 8 - --
3 Arsky 7 -- -6 -- 1 - --
4 Buinsky 7 -- -2 -- 3 2 --
5 Vysogorsky 4 -- -4 -- - - --
6 Drozhanovsky 12 -- -4 -- 7 1 - -
7 Zainsky 8 -- -7 -- - 1 - -
8 Klaibitsky 7 -- -7 -- - - --
9 Laishevsky 13 -- - 12 -- 1 - --
10 Leninogorsky 19 -- - 15 -- 4 - --
11 Mamadyshksy 10 10 -
12 Menzelinsky 6 -- -6 -- - - --
13 Musl'umovsky 4 -- -- -- 4 - --
14 Nizhnekamensky 14 -- -8 -- 5 1 - -
15 Pestrechinsky 1 -- -1 -- - - --
16 Rybnoslobodsky 8 -- -8 -- - - --
17 Sarmanovsky 2 -- -2 -- - - --
18 Spassky 9 -- -3 -- 6 - --
19 Tukaevsky 9 -- -4 -- 3 2 --
20 T'ulyachinsky 6 -- -- -- 6 - --
21 Chistopolsky 10 -- -3 -- 6 1 - -
Total number of samples 179 10 - - 106 - 55 8 - -
- 1S-c9/JQ353640/genotype 4 - 1S-c10/JQ35365Û/genotype 4
— 1 S-c 1 Q/JQ353652/genoty pe 4
- ♦ N084/KC886623/genotype 4
- ♦ N018/KC867145/genotype 4
- NK11/JQ686117/geriotype 4 ♦ NÛ25/KC867150/genotype 4
I-1 BY/HQ902258/genotype 4
♦ N029/KC888608/genotype 4
^-♦ N030/KC886609/genotype 4
I ♦ N050/KC866613/genoty pe 4
I I-♦ N054/KC886614/genotype 4
I— ♦ N094/KC886627/genotype 4
-♦ N022/KC867148/genotype 4
♦ N015/KC867143/genotype 4
1 ♦ N062/KC886615/gsnotype 4 BG/EF065638/geno1ype 4
.-3/U87872/genotype 4
1-♦ N10/HM 102355/genotype 4
-♦ N001/KC867136/genotype 4
-♦ N074/KC886620/genotype 4
- ♦ N031/KC886610/genotype 4
— ♦ N019/KC867146/genotype 4 ♦ N097/KC886628/genotype 4
— ♦NI 22/KC8S6632/genotype 4 I-♦ N023/KC867149/genotype 4
_♦ N069/KC886619/genotype 4
L + N092/KC866626/genotype 4
~ |-♦ N100/KC886630/genotype 4
_♦ N090/KC886625/geno1ype 4
U— ♦ N123/KCB86633/genotype 4
I-♦ N027/KC886S07/genotype 4
—-♦ N034/KC886611/genotype 4
' ♦ N078/KC886621/genotype 4
I-♦ N002/KC867137/genotype 4
I ♦ N007/KC867141/genotype 4
-♦ N082/KC886622/genotype 4
♦ N72/JF683619/genotype 4
I-CRLC-1/EF065655/genotype 5
j,— CFiAS-1/EF065635/genntype 5
1--— CRGC/EF065639/genotype 5
GS3/MF574055/genotype 6
-QH1/MF574057/genotype 6
I— Pucallpa-7/LC075552/genotype 6
'-Paraguay-96/LC075556/genotype S
SC2/MF574060/genotype 6 PL-1238/FJ808582/genotype 6
151/AY185360/genotype 6
Pa51-A3/KU233547/genotype 10
-ML45-B3/KU233540/genotype 10
-L1/LC154066/genotype 10
1S-c1/JQ353649/genotype 7
l2/S83530/genotype 7
-30/DÛ059417/genotype 7
♦ N003/KC867138/genotype 7 -♦ N017/KC867144/geno1ype 7
♦ N005/KC867139/genotype 7
♦ N067/KC886618/genotype 7
- 3S/JF720351/genotype 7
♦ N021/KC867147/genotype 7
— ♦ N099/KC888629/genotype 7 4T-c19/JQ353655/genotype 7 -4T-c11/JQ353656/genotype 7
1 S-c4/JQ353651/genotype 7 176/AY515276/genotype 7
— NK17/JQ686120/genotype 7
♦ N035/KC886612/genotype 7 I-♦ N28/HM1O2350/genotype 7
N013/KCS67142/genotype 7
-14/ AY515274/g e n oty pe 7
♦ N066/KC886617genoty pe 7
-4S/JF720352/geno1ype 7
-1S-c6/JQ353633/genotype 7
I-AL-1453/FJS0S577/genotype 2
AL-164/FJ808574/genotype 2
-PL-4960/FJ808590/geno1ype 2
I-ARGSF8/AF485773/genotype 2
-Monetro-1/LC075563/genotype 9
- Portachello-20/LC075567/genotype 9
- JPFU/EF065650/genotype 3
I-USCA-1/EF065647/genotype 3
I USCA-2/EF065648/genotype 3
_N121/KC886631/genotype 8
I <-♦ N142/KC886634/genotype 8
>-♦ N174/JF713455/genotype 8
♦ N174/JF713455/genotype 8
.-MKC2137/J0675759/genotype 8
1 ♦ N089/KC886624/genotype 8 ♦ N006/KC867140/genotype 8
-4-6/HM563764/genotype 8
— ♦ N063/KC886616igenotype 8
M1/ELG Cro/08/GU724606/genotype 8
-ELG Cro/VRA/09/JN990072/genotype 8
Kurtiistan/EU266062/genotype 1
_J AL-2106/FJ808578/genotype 1
~1-UruC06ll/FM955558/genotype 1
VdM/M35239/genotype 1 Cow 527/AF007764/genotype 1
1-AL-63/FJ808571/genotype 1
23/U87873/genotype 1
Fig. 1. Dendrogramme of 99 isolates of 10 BLV genotypes, based on a phylogenetic analysis of the env gene fragment [MEGA-4: algorithm NJ, 400 nt, 99 seq.] Legend: black diamond marks GenBank NCBI nucleotide sequences of the env gene fragment of BLV provirus isolates in the Republic of Tatarstan.
Donnik I.M. et al. Foods and Raw Materials, 2018, vol. 6, no. 2, pp. 314-324 Table 2. Intra- and intergenotypic heterogeneity of reference BLV representatives according to env gene (% ratio)
PHYLOGENETIC CLASSIFICATION OF BLV
GENOTYPE 1 2 3 4 5 6 7 8 9 10
1 0-5 3-6 3-7 3-7 3-7 3-7 3-8 2-6 3-6 4-7
2 3-6 0-1 3-4 2-4 4-5 3-5 3-5 2-4 3 4-6
3 3-7 3-4 0-2 3-5 4-5 3-5 3-6 3-4 3 4-6
4 3-7 2-4 3-5 0-3 3-5 2-5 2-6 2-4 2-4 3-5
5 3-7 4-5 4-5 3-5 0-2 4-6 3-5 4-5 4-5 5-6
6 3-7 3-5 3-5 2-5 4-6 0-4 3-6 2-5 3-5 3-5
7 3-8 3-5 3-6 2-6 3-5 3-6 0-4 2-5 3-5 4-6
8 2-6 2-4 3-4 2-4 4-5 2-5 2-5 0-2 2-3 3-5
9 3-6 3 3 2-4 4-5 3-5 3-5 2-3 0-1 4-5
10 4-7 4-6 4-6 3-5 5-6 3-5 4-6 3-5 4-5 0-2
Table 3. Comparison of in silico data for PCR-RFLP (typification according to D. Beier et al., 2001) and the phylogenetic analysis of the BLV env gene fragment
PCR-RFLP
PCR product
RFLP fragments (bp)
BLV genotypes
N
ping (bp) PvuII BamHI BclI 1 2 3 4 5 6 7 8 9 10
Belgian 444 280/164 444 225/219 - - - 142 - - - - - - 142
Australian 444 444 316/128 225/219 57 - 4 - - 28 70 21 19 - 199
Japanese 444 444 316/128 219/121/104 8 - - - - 6 1 - - - 15
? 444 444 444 225/219 43 - - 1 - - 2 - 3 17 66
? 444 444 316/128 444 1 - - - - 14 2 - - - 17
? 444 444 316/128 225/191/28 1 1
? 444 280/164 316/128 225/219 - 36 - - 10 - 3 - - - 49
? 444 280/164 316/128 444 - 1 1
? 444 280/164 316/128 219/189/36 - - - - 1 - - - - - 1
? 444 280/164 444 444 - - - 1 - - - - - - 1
? 444 208/164 253/191 225/219 - - - 1 - - - - - - 1
? 444 280/164 444 219/121/104 - - - 4 - - - - - - 4
? 444 444 242/128/74 225/219 - - - - - 1 - - - - 1
? 444 444 444 444 7 7
Legend. N is the number of analysed BLV isolates with an established PCR-RFLP profile; "?" - an unclassifiable BLV taxon.
Table 4. Comparison of in silico data of PCR-RFLP (typification according to M. Licursi et al., 2002) and phylogenetic analysis of a fragment of the BLV env gene
PCR-RFLP PCR product RFLP fragments (bp)
BLV genotypes
N
'ping (bp) BclI HaeIII PvuII
1 444 225/219 198/94/87/32/27/6 444
2 444 219/121/104 312/94/32/6 444
3 444 219/121/104 285/94/32/27/6 444
4 444 219/121/104 198/94/87/32/27/6 444
5 444 225/219 285/94/32/27/6 444
6 444 225/219 198/94/87/32/27/6 280/164
? 444 444 198/94/87/32/27/6 444
? 444 225/191/28 198/119/94/27/6 444
? 444 225/219 312/94/32/6 444
? 444 444 198/94/87/32/27/6 280.164
? 444 219/189/36 198/94/87/32/27/6 280/164
? 444 225/219 285/94/32/27/6 280/164
? 444 444 198/87/49/45/32/27/6 444
? 444 225/219 225/94/87/32/6 444
? 444 225/219 285/94/32/21/6/6 444
? 444 225/219 198/121/87/32/6 280/164
? 444 225/219 198/119/94/27/6 280/164
? 444 219/121/104 285/94/32/27/6 280/164
? 444 225/219 198/87/49/45/32/27/6 444
? 444 225/219 198/119/94/27/6 444
? 444 444 198/121/87/32/6 444
? 444 225/219 198/87/49/45/32/27/6 280/164
? 444 444 198/94/81/32/21/6/6 444
? 444 225/219 198/94/81/32/27/6/6 444
? 444 225/219 279/94/32/27/6/6 444
12 3 4 56789 10
98 - - 1 - 28 68 - - - 195
8 - - -- 61-- - 15
- 3 - - 1 1 - 22 -35 - 139 9 - 3 - - -
- - - - 12 - - - 7
1 - 1 ------ - 2 1 - - - -
1
- - - - - --21 -- -- -- -- -- -- 21 -- --
- - - - - 1 1 - -
6 11
28 186 20 1 1 2 1
3 1 21 1
1 1
4 2 1 2 1 1 6 11
Legend. N is the number of analysed BLV isolates with an established PCR-RFLP profile; "?" - an unclassifiable BLV taxon.
1
1
Table 5 compares the in silico data of PCR-RFLP according to the strategy by H. Fechner et al. (1997) [30] and the phylogenetic analysis of a fragment of the BLV env gene.
Thus, the BLV isolates identified according to H. Fechner et al. (1997) [30] as variant group A belong to genotype 4 according to the phylogenetic classification; variant group B - to genotypes 1, 6, or 7; variant group C - to genotypes 1, 3, 6, 7, or 9; variant group D - to genotypes 1, 4, or 7; variant group F - to genotypes 2, 5, or 7; variant group G - to genotype 1 (Table 5).
For this typing strategy [30], there are 17 unique combinations of PCR-RFLP profiles that are conditionally identical to 17 unclassifiable variant BLV groups (Table 5). Besides, BLV isolates genotyped by phylogenetic analysis as genotype 1 are characterized as variant groups B, C, D, G and three unclassifiable variant groups of BLV according to the strategy of H. Fechner et al. (1997) [30]; the genotype 2 belongs to variant group F and one unclassifiable variant group; genotype 3 - to variant group C; genotype 4 - to
variant groups A, D and three unclassifiable variant groups; genotype 5 - to variant group F and two unclassifiable variant groups; genotype 6 - to variant groups B and C and four unclassifiable variant groups; genotype 7 - to variant groups B, C, D, F and two unclassifiable variant groups; genotype 8 - to variant group E; genotype 9 - to variant group C and one unclassifiable variant group; genotype 10 - to three unclassifiable variant groups of BLV (Table 5).
The priority task of the research was to improve the strategy of PCR-RFLP genotyping of BLV by making it consistent with the phylogenetic classification and taking into account the update information on the genetic diversity of the ten known BLV genotypes.
505 BLV isolates were generated during the analysis of restriction mappings of the env gene locus according to 5 restriction enzymes. The interpretation of their env-PCR-RFLP profiles actually reflects the strategy of PCR-RFLP genotyping of BLV in accordance with its phylogenetic classification. The data are represented in Table 6.
Table 5. Comparison of in silico data of PCR-RFLP (genotyping according to H. Fechner et al., 1997) and the phylogenetic analysis of a fragment of the BLV env gene
PCR- PCR RFLP product genotyping (bp)
RFLP fragments (bp)
BLV genotypes
N
BamHI Bcll
Bgll HaelII
Pvull 12 3 4 56789 10
up o
>
A 444
B 444 C 444
D 444 E 444 F 444
G 444
?
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?
444 444 444 444 444 444 444 444 444 444 444 444 444 444 444 444 444
444
316/128 316/128
225/219
219/121/104 225/219
444
316/128 316/128 316/128 316/128 316/128 444
316/128 316/128 316/128 316/128 444
253/191 444
316/128
316/128
242/128/74
316/128
444
444
444
225/219 225/219 225/219 225/219 444
225/191/28
225/219
444
219/189/36 225/219 444 444
225/219 219/121/104 444 444
225/219 444
225/219 444
225/219
328/116 444
328/116 444
328/116 444
328/116 328/116 328/116 444
328/116 328/116 444
328/116 328/116 444
328/116 328/116 328/116 328/116 328/116 444
328/116 444
328/116 328/116 444
198/94/87/32/27/6 198/121/87/32/6 198/119/94/27/6 285/94/32/27/6 285/94/32/27/6 444
198/94/87/32/27/6 444 285/94/32/21/6/6 285/94/32/27/6 198/87/49/45/32/27/6 198/119/94/27/6 198/94/87/32/27/6 225/94/8732/6 198/94/87/32/27/6 312/94/32/6 198/94/87/32/27/6 198/119/94/27/6 285/94/32/27/6 198/94/87/32/27/6 198/94/87/32/27/6 285/94/32/27/6 198/87/49/45/32/27/6 444
280/164 - -
142 ■
8 -
56 - 4
444 42 -
444 - -
280/164 - 36
444 1 -
444 1 -
444 1 -
444 1 -
280/164 - 1
280/164 - -
280/164 - -
6 1
28 70
- - 142
- - 15
19 - 177
- - 2 - - -- - - 21 - -9 - 3 - - -
45 21 48
2 -
3 -
198/94/87/32/27/6 198/94/87/32/27/6 285/94/32/27/6 198/121/87/32/6 198/94/87/32/27/6 198/94/87/32/27/6 198/121/87/32/6 198/94/81/32/27/6/6 444 198/94/81/32/27/6/6 444 279/94/32/27/6/6 444
280/164 - -
280/164 - -
280/164 - -
444 - -
444 - -
444 - -
444 - -
11
1
1
11 11
Legend. N is the number of analysed BLV isolates with an established PCR-RFLP profile; "?" - an unclassifiable BLV taxon.
1
1
1
1
Donnik I.M. et al. Foods and Raw Materials, 2018, vol. 6, no. 2, pp. 314-324 Table 6. An improved strategy for PCR-RFLP-genotyping of BLV, consistent with the phylogenetic classification
G BLV isolate GenBank PCR _RFLP fragments (bp)_C N
A/N product PvuII SspI HphI HaeIII BstYI
(bp)
1 AL-63 FJ808571 444 444 399/45 224/220 198/94/87/32/27/6 198/128/118 1 56
1 Cow 527 AF007764 444 444 399/45 224/220 285/94/32/27/6 198/128/118 2 8
1 23 U87873 444 444 399/45 224/220 312/94/32/6 198/128/118 3 1
1 AL-2106 FJ808578 444 444 399/45 224/220 198/94/87/32/27/6 246/198 4 42
1 UruC06II FM955558 444 444 399/45 224/220 285/94/32/27/6 246/198 5 1
1 VdM M35239 444 444 399/45 224/181/39 198/94/87/32/27/6 316/128 6 1
1 Kurdistan EU266062 444 444 399/45 220/196/28 198/119/94/27/6 198/128/118 7 1
2 AL-164 FJ808574 444 280/164 399/45 224/220 198/94/87/32/27/6 198/128/118 8 34
2 PL-4960 FJ808590 444 280/164 399/45 224/220 198/87/49/45/32/27/6 198/128/118 9 1
2 ARGSF8 AF485773 444 280/164 399/45 444 198/94/87/32/27/6 198/128/118 10 1
2 AL-1453 FJ808577 444 280/164 444 224/220 198/94/87/32/27/6 198/128/118 11 1
3 USCA-1 EF065647 444 444 399/45 444 285/94/32/21/6/6 198/128/96/22 12 1
3 USCA-2 EF065648 444 444 399/45 444 285/94/32/27/6 198/128/96/22 13 2
3 JPFU EF065650 444 444 399/45 444 285/94/32/27/6 198/128/118 14 1
4 BG EF065638 444 280/164 399/45 224/220 198/94/87/32/27/6 444 15 115
4 3 U87872 444 444 399/45 224/220 198/94/87/32/27/6 444 16 1
4 1S-c16 JQ353652 444 280/164 399/45 444 198/94/87/32/27/6 444 17 16
4 N023 KC867149 444 280.164 399.45 224.220 198/94/87/32/27/6 253/191 18 1
4 1 BY HQ902258 444 280/164 444 224/220 198/94/87/32/27/6 444 19 7
4 N034 KC886611 444 280/164 399/45 224/220 198/121/87/32/6 444 20 1
4 1S-c9 JQ353640 444 280/164 399/45 224/220 198/119/94/27/6 444 21 1
4 NK11 JQ686117 444 280/164 399/45 224/220 285/94/32/27/6 444 22 6
4 1S-c10 JQ353650 444 280/164 399/45 220/145/79 198/94/87/32/27/6 444 23 1
5 CRAS-1 EF065635 444 280/164 399/45 224/181/39 198/94/87/32/27/6 316/128 24 8
5 CRGC EF065639 444 280/164 399/45 224/181/39 285/94/32/27/6 316/128 25 1
5 CRLC-1 EF065655 444 280/164 444 224/181/39 198/94/87/32/27/6 316/128 26 2
6 PL-1238 FJ808582 444 444 399/45 224/220 285/94/32/27/6 316/128 27 7
6 151 AY185360 444 444 399/45 224/220 198/94/87/32/27/6 316/128 28 27
6 GS3 MF574055 444 444 399/45 444 198/94/87/32/27/6 316/128 29 11
6 SC2 MF574060 444 444 399/45 224/220 198/94/87/32/27/6 242/128/74 30 1
6 QH1 MF574057 444 444 213/186/45 444 198/94/81/32/21/6/6 316/128 31 1
6 Pucallpa-7 LC075552 444 444 399/45 444 198/94/87/32/27/6 316/79/49 32 1
6 Paraguay-96 LC075556 444 444 399/45 444 198/121/87/32/6 316/128 33 1
7 N28 HM102356 444 444 444 224/137/83 198/94/87/32/27/6 294/128/22 34 7
7 176 AY515276 444 444 444 224/137/83 198/94/87/32/27/6 316/128 35 53
7 I2 S83530 444 444 444 224/220 285/94/32/27/6 316/128 36 1
7 14 AY515274 444 444 444 145/137/83/79 198/94/87/32/27/6 316/128 37 1
7 30 DQ059417 444 444 444 444 198/87/49/45/32/27/6 316/128 38 1
7 3S JF720351 444 280/164 444 224/137/83 198/94/87/32/27/6 316/128 39 3
7 4T-c19 JQ353655 444 444 399/45 224/137/83 198/94/87/32/27/6 316/128 40 3
7 1S-c4 JQ353651 444 444 444 224/137/83 198/94/87/32/27/6 316/79/49 41 1
7 NK17 JQ686120 444 444 444 224/137/83 198/87/49/45/32/27/6 316/128 42 2
7 4S JF720352 444 444 444 224/137/83 198/119/94/27/6 316/128 43 1
7 1S-c6 JQ353633 444 444 444 224/137/83 198/121/87/32/6 316/128 44 1
7 4T-c11 JQ353656 444 444 444 224/137/83 285/94/32/27/6 316/128 45 1
7 N067 KC886618 444 444 444 224/137/44/39 198/94/87/32/27/6 316/128 46 1
7 1S-c1 JQ353649 444 444 444 224/220 198/94/87/32/27/6 444 47 2
8 M1/ELG Cro/08 GU724606 444 444 399/45 224/220 225/94/87/32/6 198/128/118 48 13
8 N174 JF713455 444 444 399/45 224/220 225/94/87/32/6 316/128 49 4
8 ELG Cro/VRA/09 JN990072 444 444 444 224/220 225/94/87/32/6 198/128/118 50 2
8 4-6 HM563764 444 444 399/45 224/137/83 225/94/87/32/6 198/128/118 51 1
8 MKC2137 JQ675759 444 444 399/45 444 225/94/87/32/6 198/128/118 52 1
9 Monetro-1 LC075563 444 444 399/45 224/171/49 285/94/32/27/6 198/128/118 53 19
9 Portachello-20 LC075567 444 444 399/45 224/171/49 285/94/32/27/6 246/198 54 3
10 Pa51-A3 KU233547 444 444 399/45 224/220 198/94/81/32/27/6/6 444 55 12
10 ML45-B3 KU233540 444 444 399/45 224/220 279/94/32/27/6/6 444 56 11
10 L1 LC154066 444 444 444 224/220 198/94/81/32/27/6/6 444 57 1
Legend. G - genotype; C - combination; N - the number of analysed BLV isolates with the established combination of PCR-RFLP profiles.
Fig. 2. Electrophoregramme of combinations of PCR-RFLP profiles (C1 and C4) of genotype 1 (an improved BLV genotyping strategy).
Legend. 1, 7) DNA markers 100 bp + 50 bp (SibEnzyme). 2-6) PCR-RFLP profile of BLV N-1 provirus isolate (C1, genotype 1): 2) PvuII-RFLP (444 bp); 3) SspI-RFLP (399/45 bp); 4) HphI-RFLP (224/220 bp); 5) HaeIII-RFLP (198/94/87/32/27/6 bp); 6) BstYI-RFLP (198/128/118 bp). 8-12) PCR-RFLP profile of BLV N-4 provirus isolate (C4, genotype 1): 8) PvuII-RFLP (444 bp); 9) SspI-RFLP (399/45 bp); 10) HphI-RFLP (224/220 bp); 11) HaeIII-RFLP (198/94/87/32 / 27.6 bp); 12) BstYI-RFLP (246/198 bp).
Fig. 3. Electrophoregramme of combinations of PCR-RFLP profiles (C17 and C18) of genotype
4 (improved BLV genotyping strategy).
Legend. 1, 7) DNA markers 100 bp + 50 bp (SibEnzyme). 2-6) PCR-RFLP-profile of the BLV N015 provirus isolate (C17, genotype 4): 2) PvwII-RFLP (280/164 bp); 3) SspI-RFLP (399/45 bp); 4) HphI-RFLP (444 bp); 5) HaeIII-RFLP (198/94/87/32/27/6 bp); 6) BstYI-RFLP (444 bp); 8-12) PCR-RFLP-profile of the BLV N023 provirus isolate (C18, genotype 8) PvwII-RFLP (280/164 bp); 9) SspI-RFLP (399/45 bp); 10) HphI-RFLP (224/220 bp); 11) HaeIII-RFLP (198/94/87/32 / 27.6 bp); 12) BstYI-RFLP (253/191 bp).
The PCR-RFLP-genotyping strategy of BLV, which we have improved, is consistent with its phylogenetic classification. The new classification makes it possible to identify all the ten currently known genotypes of the viral pathogen (Table 6).
It should be noted that genotype 1 is characterized by seven combinations of env-PCR-RFLP profiles (C 1-7), genotype 2 - by four combinations (C 8-11); genotype 3 - by three combinations (C 12-14); genotype 4 - by nine combinations (C15-23); genotype
5 - by three combinations (C 24-26); genotype 6 - by seven combinations (C27-33), genotype 7 - by fourteen combinations (C 34-47); genotype 8 - by five combinations (C 48-52); genotype 9 - by two combinations (C 53-54); genotype 10 - by 3 three combinations (C 55-57) (Table 6).
Fig. 4. Electrophoregramme of combinations of PCR-RFLP profiles (C46 and C48) of genotypes 7 and 8 (an improved BLV genotyping strategy).
Legend. 1, 7) DNA markers 100 bp + 50 bp (SibEnzyme). 2-6) PCR-RFLP-profile of the BLV N067 provirus isolate (C46, genotype 7): 2) PvuII-RFLP (444 bp); 3) SspI-RFLP (444 bp); 4) HphI-RFLP (224/137/44/39 bp); 5) HaeIII-RFLP (198/94/87/32/27/6 bp); 6) BstYI-RFLP (316/128 bp). 8-12) PCR-RFLP profile of the BLV N006 provirus isolate (C48, eighth genotype 8) PvuII-RFLP (444 bp); 9) SspI-RFLP (399/45 bp); 10) HphI-RFLP (224/220 bp); 11) HaeIII-RFLP (225/94/87/32/6 bp); 12) BstYI-RFLP (198/128/118 bp).
It should be emphasized that genotypes 8 and 9 can be easily identified even with the use of one restrictase, HaeIII, generating RFLP fragments (225/94/87/32/6 bp) that are characteristic of genotype 8; HphI -generating RFLP fragments (224/171/49 bp) that are characteristic of genotype 9. Representatives of genotypes 2 (BstYI and PvuII), 3 (HaeIII and HphI), and 5 (HphI and PvuII) can be identified with two restriction enzymes (Table 6).
Figs. 2-4 show illustrative examples of the implementation of the strategy of BLV PCR-RFLP-genotyping in accordance with its phylogenetic classification.
As one can see from the electrophoregramme in Fig. 2 (tracks 2-6), the PCR-RFLP profile of the BLV N-1 provirus isolate is identified as combination 1 (C1) of the env-PCR-RFLP profile of genotype 1, which includes at least 56 isolates deposited in the GenBank NCBI (Table 6).
The PCR-RFLP profile of the BLV N-4 provirus isolate (Fig. 2, tracks 8-12) characterizes combination 4 (C4) of the env-PCR-RFLP profile of genotype 1, with at least 42 identified representatives (Table 6).
The PCR-RFLP profile of the BLV N015 provirus isolate (Fig. 3, tracks 2-6) is identified as combination
17 (C17) of the env-PCR-RFLP profile of genotype 4. It includes at least 16 representatives (Table 6), two of which affect cattle in Tatarstan. According to GenBank NCBI, these nucleotide sequences of the env gene fragment are isolate N015 (GenBank A/N: KC867143) and isolate N062 (GenBank A/N: KC886615).
The PCR-RFLP profile of the BLV N023 provirus isolate (Fig. 3, tracks 8-12) is identified as combination
18 (C18) of the env-PCR-RFLP profile of genotype 4 (Table 6). Its nucleotide sequence of the env gene fragment from the GenBank NCBI is the only variant for the given combination (isolate N023, GenBank A/N: KC867149).
- GS3/MF574055/genotype 6
-QH1/MF574057/genotype 6
Pucallpa-7/LC075552/genotype 6 Paraguay-96/LC075556/genotype 6
-SC2/MF574060/genotype 6
- PL-1238/FJ808582/genotype 6
- 151/AY185360/genotype 6
- Pa51-A3/KU233547/genotype 10
- ML45-B3/KU233540/genotype 10 -L1/LC154066/genotype 10
■CE
Ct
-1Si16/JQ353652/genotype4
-1 S-c9/JQ353640/genotype 4
-NK11/JQ686117/genotyps 4
-1S-c10/JQ353650/genotype 4
-N023/KC8S7149/genotype 4
-N034/KC88S611/genotype 4
-1 BY/HO902258/genotype 4
BG/EF065638/genotype 4
-3/U87872/genotype 4
j-CRAS-1/EF06S635/genotype 5
- CRGC/EF065639/genotype 5
-CRLC-1/EF065655/genotype 5
-l2/S83530/genotype 7
1S-c1/JQ353649/genotype 7
- 30/DQ059417/genotype 7
N067/KC886618/genotype 7
J-4T-c19/JQ353655/genotype 7
"I-4T-C11/JQ353S56/genotype 7
— NK17/JQ686120/genotype 7
-1S-c4/JQ353651/genotype 7
-3S/JF720351/genotype 7
-N28/HM102356/genotype 7
-14/AY515274/genotype 7
-1S-c6/JQ353633/genotype 7
176/AY515276/genotype 7 -4S/JF720352/genotype 7
-PL-4960/FJ808590/genotype 2
-AL-1453/FJ808577/genotype 2
AL-164/FJ808574/genotype 2 -ARGSF8/AF485773/genotype 2
4-6/HM 563764/genotype 8
N174/JF713455/genotype 8 MKC2137/JQ675759/genotype 8
M1/ELG Cro/08/GU724506/genotype 8
-ELG Cro/VRA/09/JN990072/genotype 8
I-JPFU/EF065650/genotype 3
-|i-USCA-1/EF065647/genotype 3
' USCA-2/EF065648/geno1ype 3
_I-MonetrD-1/LC075563/genotype 9
I— Portachello-20/LC075567/genotype 9
-Kurdi3tan/EU266062/gBnotypD 1
_[AL-2106/FJ808578/genotype 1
I-UruC06ll/FM955558/genotype 1
-VdM/M35239/genotype 1
-Cow 527/AF007764/genotype 1
- AL-63/FJ808571/genotype 1 -23/U87873/genotype 1
Fig. 5. Phylogramme of 57 reference isolates of the 10 currently known BLV genotypes, built on the basis of a phylogenetic analysis of the env gene fragment [MEGA-4: NJ algorithm, 400 nt, 57 seq.].
The PCR-RFLP profile of the BLV N067 provirus isolate (Fig. 4, tracks 2-6) characterizes itself as combination 46 (C46) of the env-PCR-RFLP profile of genotype 7 (Table 6). Its nucleotide sequence of the env gene fragment from GenBank NCBI is the only one for this combination (isolate N067, GenBank A/N: KC886618).
The PCR-RFLP profile of the BLV N006 provirus isolate (Fig. 4, tracks 8-12) belongs to combination 48 (C48) of the env-PCR-RFLP profile of genotype 8. It includes at least 13 representatives (Table 6), three of which affect cattle populations in Tatarstan. According to GenBank NCBI, their nucleotide sequences of the env gene fragment are isolate N063 (GenBank A/N: KC886616), isolate N006 (GenBank A/N: KC867140), and isolate N089 (GenBank A/N: KC886624).
The improved strategy of PCR-RFLP-genotyping corresponds with the modern phylogenetic classification of BLV and makes it possible to identify all its known genotypes. Its accuracy is based upon in silico modelling of restrictogrammes and the phylogenetic analysis of the env gene fragment of 57 reference isolates of the ten known BLV genes (Fig. 5). They produce 57 diagnostically significant genotype-associated combinations of PCR-RFLP profiles.
CONCLUSION
To determine the genotypes of BLV isolates obtained from Tatarstani cattle, we performed a phylogenetic analysis of the env gene fragment sequences and a PCR-RFLP analysis that corresponded with the phylogenetic classification of the infectious agent. The genotypic composition of 179 identified BLV isolates detected in cattle from livestock farms of 21 districts of the Republic of Tatarstan was represented by genotypes 1 (10 isolates), 4 (106 isolates), 7 (55 isolates), and 8 (8 isolates). Thus, we state the fact that four out of ten currently known BLV genotypes circulate on the territory of the Republic of Tatarstan, namely genotypes 1, 4, 7, and 8.
After that, we classified the BLV isolates with decoded nucleotide sequences of the env gene locus according to the chosen genetic identification strategy. Subsequently, we assessed the degree of consistency of genotypic approaches by comparing in silico PCR-RFLP data and the results of the phylogenetic analysis. We used 505 corresponding sequences, including those deposited in GenBank NCBI. As a result, we managed to prove that a number of previously used PCR-RFLP typing strategies were inconsistent with the current approach in assessing the genotypic diversity of BLV with the help of the phylogenetic analysis. The inconsistency of the three PCR-RFLP strategies for BLV typing with the modern phylogenetic classification is associated, among other things, with the on-going knowledge acquisition in the sphere of the genetic diversity of the ten known BLV genotypes.
During the final stage of the research, we improved the strategy of PCR-RFLP-genotyping of BLV to make it consistent with the phylogenetic classification. The new version takes into account the new data about the genetic diversity of BLV. It also includes an interpretation of the PCR-RFLP profiles of 505 BLV isolates. The interpretation resulted from a restriction mapping of the env gene fragment according to 5 restriction endonucleases. The improved strategy of PCR-RFLP-genotyping allows one to identify all the currently known BLV genotypes. The improved strategy owes its accuracy to in silico modelling of restrictogrammes and the phylogenetic analysis of the env gene of 57 reference isolates of ten BLV genes that generate 57 diagnostically significant genotype-associated combinations of PCR-RFLP profiles
CONFLICT OF INTEREST
The authors declare that there are no conflict of interest related to this article.
The research was supported by the section of storage and processing of agricultural products of the
FUNDING
Department of Agricultural Sciences (Russian Academy of Sciences) and V.M. Gorbatov Federal Research Center for Food Systems.
REFERENCES
1. Mishchenko V.A., Petrova O.N., Karaulov A.K., and Mishchenko A.V. Problema leykoza krupnogo rogatogo skota [The issues of bovine leukaemia]. Vladimir: FGBI ARRIAH Publ., 2018. 38 p. (In Russ.).
2. Donnik I.M., Petropavlovsky M.V., Krivonogova A.S., et al. Revisiting the issue of the molecular-genetic structure of the causative agent of the bovine leukemia virus in the Russian Federation. Indian Journal of Science and Technology, 2016, vol. 9, no. 42. DOI: https://doi.org/10.17485/ijst/2016/v9i42/104253.
3. Olaya-Galán N.N., Corredor-Figueroa A.P., Guzmán-Garzón T.C., et al. Bovine leukaemia virus DNA in fresh milk and raw beef for human consumption. Epidemiology & Infection, 2017, vol. 145, no. 15, pp. 3125-3130. DOI: https://doi.org/10.1017/S0950268817002229.
4. Villalobos-Cortes A. Enzootic bovine leukosis and the risk to human health. African Journal of Biotechnology, 2017, vol. 16, no. 15, pp. 763-770. DOI: https://doi.org/10.5897/AJB2016.15736.
5. Alabbody H.H.K. Bovine leukosis and the possibility to cause cancer in humans: A scientific review. Iraqi Journal of Veterinary Medicine, 2018, vol. 42, no. 1, pp. 52-60.
6. Reichert M., Grundboeck-Jusko J., Rulka J., Stec J., and Kozaczynski W. Influence of selected technological treatments on the BLV provirus DNA occurring in the tissues of leukaemic cattle. Bulletin of the Veterinary Institute in Pulawy, 1994, vol. 38, no. 2, pp. 52-60.
7. Sviridenko G.M. and Semova E.G. Leucose in the cattle and safety of milk products. Dairy industry, 2003, no. 7, pp. 8-10. (In Russ.).
8. Mesa G., Ulloa J.C., Uribe A.M., and Gutierrez M.F. Bovine Leukemia Virus Gene Segment Detected in Human Breast Tissue. Open Journal of Medical Microbiology, 2013, vol. 3, no. 1, pp. 84-90. DOI: https://doi.org/10.4236/ojmm.2013.31013.
9. Buehring G.C., Shen H.M., Jensen H.M., et al. Bovine Leukemia Virus DNA in Human Breast Tissue. Emerging Infectious Diseases, 2014, vol. 20, no. 5, pp. 772-782. DOI: https://doi.org/10.3201/eid2005.131298.
10. Buehring G.C., Shen H.M., Jensen H.M., et al. Exposure to Bovine Leukemia Virus Is Associated with Breast Cancer: A Case-Control Study. PLoS One, 2015, vol. 10, no. 9. DOI: https://doi.org/10.1371/journal.pone.0134304.
11. Buehring G.C., Shen H., Schwartz D.A., and Lawson J.S. Bovine leukemia virus linked to breast cancer in Australian women and identified before breast cancer development. PLoS One, 2017, vol. 12, no. 6. DOI: https://doi.org/10.1371/journal.pone.0179367.
12. Khazipov N.Z., Vafin R.R., Shaeva A.Yu., and Zaynullin L.I. Transformation of cells under the influence of bovine leukemia virus is a real risk of developing the oncological diseases in humans. Modern problems of science and education, 2013, no. 6. Available at: http://www.science-education.ru/ru/article/view?id=11792. (accessed 27 July 2018). (In Russ.).
13. Petrov A.N., Galstyan A.G., Radaeva I.A., et al. Indicators of quality of canned milk: Russian and international priorities. Foods and Raw Materials, 2017, vol. 5, no. 2, pp. 151-161. DOI: https://doi.org/10.21179/2308-4057-2017-2-151-161.
14. Petrov A.N., Khanferyan R.A., and Galstyan A.G. Current aspects of counteraction of foodstuff's falsification. Problems of Nutrition, 2016, vol. 85, no. 5, pp. 86-92. (In Russ.).
15. Galstyan A.G., Radaeva I.A., Chervetsov V.V., et al. Improvement of the canned milk products quality due to the application of the pasteurized raw milk. Dairy industry, 2015, no. 5, pp. 42-44. (In Russ.).
16. Kuznetsova T.V., Kuznetsov A.A., and Kirillova S.V. Tekhnicheskiy reglamennt TS "O bezopasnosti moloka i molochnoy produktsii" i ehkonomicheskie aspekty ego realizatsii molochnymi tovaroproizvoditelyami Rossii [On safety of milk and dairy products' and the economic aspects of its implementation by Russian dairy producers].
Agro-food policy in Russia, 2015, vol. 48, no. 12, pp. 35-38. (In Russ.).
17. Lysov A., Petropavlovskiy M., Donnik I., and Krivonogova A. The system of individual veterinary and zootechnie activities for healthy improvement from BLV of cattle on the example of Tyumen region. Issues of Legal Regulation in Veterinary Medicine, 2017, no. 3, pp. 40-43. (In Russ.).
18. Vinogradova I.V., Gladyr' E.A., Kovalyuk N.V., et al. Izuchenie genotipicheskogo raznoobraziya, molekulyarno-geneticheskoy struktury i filogeneticheskiy analiz vozbuditelya leykoza krupnogo rogatogo skota, tsirkuliruyushchego v populyatsii zhivotnykh iz raznykh regionov Rossii [A study of the genotypic diversity, molecular-genetic structure and a phylogenetic analysis of the causative agent of bovine leukaemia in cattle from different regions of Russia]. Ekaterinburg: Ural State Agrarian University Publ., 2014. pp. 3-18. (In Russ.).
19. Vafin R.R., Khazipov N.Z., Shaeva A.Y., et al. Genotypic identification of the bovine leukemia virus.
Molecular Genetics, Microbiology and Virology, 2014, vol. 29, no. 4, pp. 195-203. DOI: https://doi.org/10.3103/S0891416814040120.
20. Rodriguez S.M., Golemba M.D., Campos R.H., Trono K., and Jones L.R. Bovine leukemia virus can be classified into seven genotypes: evidence for the existence of two novel clades. Journal of General Virology, 2009, vol. 90, no. 11, pp. 2788-2797. DOI: https://doi.org/10.1099/vir0.011791-0.
21. Shaeva A.Y., Vafin R.R., Khazipov N.Z., et al. Genotypic identification of BLV isolates detected in livestock farms of the republic of Tatarstan. Uchenye zapiski Kazanskoy gosudarstvennoy akademii veterinarnoy meditsiny im. N.EH. Baumana [Scientific notes of the Kazan State Academy of Veterinary Medicine. N.E. Bauman], 2011, vol. 208, pp. 330-337. (In Russ.).
22. Shaeva A.Y., Garaeva Z.R., Vafin R.R., Khazipov N.Z., and Alimov A.M. Identification of a new BLV genotype.
Uchenye zapiski Kazanskoy gosudarstvennoy akademii veterinarnoy meditsiny im. N.EH. Baumana [Scientific notes of the Kazan State Academy of Veterinary Medicine. N.E. Bauman], 2012, vol. 211, pp. 192-197. (In Russ.).
23. Vafin R.R., Khazipov N.Z., Shaeva A.Y., et al. Finding a new (8th) genotype of BLV in various regions of the world. Fundamental research, 2013, no. 10, part 7, pp. 1467-1471. (In Russ.).
24. Balic D., Lojkic I., Periskic M., et al. Identification of a new genotype of bovine leukemia virus. Archives of Virology, 2012, vol. 157, no. 7, pp. 1281-1290. DOI: https://doi.org/10.1007/s00705-012-1300-4.
25. Rola-Luszczak M., Pluta A., Olech M., et al. The molecular characterization of bovine leukaemia virus isolates from Eastern Europe and Siberia and its impact on phylogeny. PLoS One, 2013, vol. 8, no. 3. DOI: https://doi.org/10.1371/journal.pone.0058705.
26. Polat M., Takeshima S.N., Hosomichi K., et al. A new genotype of bovine leukemia virus in South America identified by NGS-based whole genome sequencing and molecular evolutionary genetic analysis. Retrovirology. 2016, vol. 13, no. 4. DOI: https://doi.org/10.1186/s12977-016-0239-z.
27. Lee E., Kim E.J., Ratthanophart J., et al. Molecular epidemiological and serological studies of bovine leukemia virus (BLV) infection in Thailand cattle. Infection, Genetics and Evolution, 2016, vol. 41, no. 245-254. DOI: https://doi.org/10.1016Zj.meegid.2016.04.010.
28. Beier D., Blankenstein P., Marquardt O., and Kuzmak J. Identification of different BLV provirus isolates by PCR, RFLPA and DNA sequencing. Berliner und Munchener Tierarztliche Wochenschrift, 2001, vol. 114, no. 7-8, pp. 252-256.
29. Licursi M., Inoshima Y., Wu D., et al. Genetic heterogeneity among bovine leukemia virus genotypes and its relation to humoral responses in hosts. Virus Research, 2002, vol. 86, no. 1-2, pp. 101-110.
30. Fechner H., Blankenstein P., Looman A.C., et al. Provirus Variants of The Bovine Leukemia Virus and Their Relation to the Serological Status of Naturally Infected Cattle. Virology, 1997, vol. 237, no. 2, pp. 261-269. DOI: https://doi.org/10.1006/viro.1997.8784.
ORCID IDs
Irina M. Donnik https://orcid.org/0000-0001-5611^1427 RamilR. Vafin https://orcid.org/0000-0003-0914-0053 Aram G. Galstyan https://orcid.org/0000-0002-0786-2055 Anna S. Krivonogova https://orcid.org/0000-0001-9112-0830 Aigul Y. Shaeva https://orcid.org/0000-0003-3623-0791 KhamidKh. Gilmanov https://orcid.org/0000-0001-7053-6925 Rufiya G. Karimova https://orcid.org/0000-0002-7022-9498 Sergey V. Tyulkin https://orcid.org/0000-0001-5379-237X Jacek Kuzmak https://orcid.org/0000-0003-1021-8276
