# Inter-laboratory comparison of eleven quantitative or digital PCR assays for detection of proviral bovine leukemia virus in blood samples

Aneta Pluta, Juan Pablo Jaworski, Casey Droscha, Sophie Vanderweele, Tasia

M. Taxis, Stephen Valas, Dragan Brnić, Andreja Jungić, María José Ruano, Azucena Sánchez, et al.

► To cite this version:

Aneta Pluta, Juan Pablo Jaworski, Casey Droscha, Sophie Vanderweele, Tasia M. Taxis, et al. Inter-laboratory comparison of eleven quantitative or digital PCR assays for detection of proviral bovine leukemia virus in blood samples. BMC Veterinary Research, 2024, 20 (1), pp.381. 10.1186/s12917-024-04228-z. https://doi.org/10.1186/s12917-024-04228-z

HAL Id: anses-04687241

https://anses.hal.science/anses-04687241v1

Submitted on 4 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Inter-laboratory comparison of eleven quantitative or digital PCR assays for detection of proviral bovine leukemia virus in blood samples

Aneta Pluta¹,¹³*, Juan Pablo Jaworski², Casey Droscha³, Sophie VanderWeele³, Tasia M. Taxis⁴, Stephen Valas⁵, Dragan Brnić⁶, Andreja Jungić⁶, María José Ruano⁷, Azucena Sánchez⁷, Kenji Murakami⁸, Kurumi Nakamura⁸, Rodrigo Puentes⁹, MLaureana De Brun⁹, Vanesa Ruiz², Marla Eliana Ladera Gómez¹⁰, Pamela Lendez¹⁰, Guillermina Dolcini¹⁰, Marcelo Fernandes Camargos¹¹, Antônio Fonseca¹¹, Subarna Barua¹², Chengming Wang¹², Aleksandra Giza¹³ and Jacek Kuźmak¹

**Abstract**

Bovine leukemia virus (BLV) is the etiological agent of enzootic bovine leukosis and causes a persistent infection that can leave cattle with no symptoms. Many countries have been able to successfully eradicate BLV through improved detection and management methods. However, with the increasing novel molecular detection methods there have been few efforts to standardize these results at global scale. This study aimed to determine the interlaboratory accuracy and agreement of 11 molecular tests in detecting BLV. Each qPCR/ddPCR method varied by target gene, primer design, DNA input and chemistries. DNA samples were extracted from blood of BLV-seropositive cattle and lyophilized to grant a better preservation during shipping to all participants around the globe. Twenty nine out of 44 samples were correctly identified by the 11 labs and all methods exhibited a diagnostic sensitivity between 74 and 100%. Agreement amongst different assays was linked to BLV copy numbers present in samples and the characteristics of each assay (i.e., BLV target sequence). Finally, the mean correlation value for all assays was within the range of strong correlation. This study highlights the importance of continuous need for standardization and harmonization amongst assays and the different participants. The results underscore the need of an international calibrator to estimate the efficiency (standard curve) of the different assays and improve quantitation accuracy. Additionally, this will inform future participants about the variability associated with emerging chemistries, methods, and technologies used to study BLV. Altogether, by improving tests performance worldwide it will positively aid in the eradication efforts.

**Keywords** Bovine leukemia virus (BLV), Quantitative real-time PCR (qPCR), Proviral DNA, DdPCR, BLV international network, Update on the efforts in harmonization qPCR

<footnote>*Correspondence:
Aneta Pluta
aneta.pluta@piwet.pulawy.pl
Full list of author information is available at the end of the article</footnote>
## Introduction

Bovine leukemia virus (BLV) is a deltaretrovirus from the *Orthoretrovirinae* subfamily of the *Retroviridae* family. An essential step in the BLV replication cycle is the integration of DNA copy of its RNA genome into the DNA of a host cell [1]. Once integrated, the proviral DNA is replicated along with the host's DNA during cellular divisions, as for any cellular gene. The BLV is the etiologic agent of enzootic bovine leukemia (EBL). BLV causes a persistent infection in cattle, and in most cases this infection is asymptomatic [2]. In one-third of infected animals the infection progresses to a state of persistent lymphocytosis, and in 1 to 10% of infected cattle it develops into lymphosarcoma [2]. BLV induces high economic losses due to trade restrictions, replacement cost, reduced milk production, immunosuppression, and increased susceptibility to pneumonia, diarrhea, mastitis, and so on [3–6]. BLV is globally distributed with a high prevalence, except for Western Europe and Oceania, where the virus has been successfully eradicated through detection and elimination of BLV-infected animals [7, 8]. The agar gel immunodiffusion and ELISA for the detection of BLV-specific antibodies in sera and milk are the World Organization for Animal Health (WOAH, founded as OIE) prescribed tests for serological diagnosis but ELISA, due to its high sensitivity and ability to test many samples at a very low cost, is highly recommended [9]. Despite the advantages of serologic testing, there are some scenarios in which direct detection of the BLV genomic fragment was important to improve BLV detection. The most frequent cases is the screening of calves with maternal antibodies, acute infection, animals without persistent antibody response and animal subproducts (i.e., semen). In this regard, nucleic acid amplification tests such as real-time quantitative PCR (qPCR) allows for a rapid and highly sensitive detection of BLV proviral DNA (BLV DNA) that can be used to test infected and asymptomatic animals, before the elicitation of anti-BLV specific antibodies and when proviral load (PVL) are still low [10]. Furthermore, qPCR assays can serve as confirmatory tests for the clarification of inconclusive and discordant serological test results usually associated with these cases [11]. For these reasons, the inclusion of qPCR in combination with other screening tests might increase control programs efficiency. Additionally, qPCR allows the estimation of BLV PVL which is important for studying the dynamics of BLV infection (i.e., basic research). Further, considering that BLV PVL correlates with the risk of BLV transmission, this feature of qPCR can be exploited for developing rational segregation programs [12, 13]. The results of Kobayashi et al. suggest that high PVL is also a significant risk factor for progression to EBL and should therefore be used as a parameter to identify cattle for culling from

the herd well before EBL progression [14]. Several qPCRs have been developed globally for the quantitation of BLV DNA. Although most assays have been properly validated by each developer, a proper standardization and harmonization of such tests is currently lacking. Considering that standardization and harmonization of qPCR methods and results are essential for comparisons of data from BLV laboratories around the world, this could directly impact international surveillance programs and collaborative research. We built a global collaborative network of BLV reference laboratories to evaluate the interlaboratory variability of different qPCRs and sponsored a harmonization of assays to hopefully impact international surveillance programs and research going forward.

In 2018 we conducted the first global trial of this kind to assess the interlaboratory variability of six qPCRs for the detection of BLV DNA [15]. Since this complex process is a continuous rather than a one-time effort, we now started a second study of this type. In this follow up study, we built a more comprehensive sample panel, accounting for a broader geographical diversification. Additionally, we increased the number of participants to ten collaborating laboratories plus one WOAH reference lab and tested novel methodologies including digital PCR (ddPCR) and FRET-qPCR. Finally, we established the next steps towards the international standardization of molecular assays for the detection of BLV DNA.

## Materials and methods

### Participants

The eleven laboratories that took part in the study were:(i) the Auburn University College of Veterinary Medicine (Auburn, Alabama, United States); (ii) Antel-Bio, a division of CentralStar Cooperative (Michigan, United States); (iii) Laboratórios Federais de Defesa Agropecuária de Minas Gerais (LFDA-MG, Pedro Leopoldo, Brasil); (iv) Centro de Investigación Veterinaria de Tandil (CIVETAN, Buenos Aires, Argentina); (v) the Faculty of Agriculture Iwate University (Iwate, Japan); (vi) Universidad de la República de Uruguay (UdelaR, Montevideo, Uruguay); (vii) the Croatian Veterinary Institute (Zagreb, Croatia); (viii) Instituto Nacional de Tecnología Agropecuaria (INTA, Buenos Aires, Argentina); (ix) Laboratório Central de Veterinaria (LCV, Madrid, Spain); (x) the National Veterinary Research Institute (NVRI, Puławy, Poland) and (xi) the French Agency for Food, Environmental and Occupational Health and Safety (Anses, Niort, France). All European laboratories participating in this study are acting as national reference laboratories for EBL, NVRI acts as WOAH reference laboratory for EBL, while the remaining laboratories are nationally renowned entities for BLV diagnostics. The eleven participating
methods are referred to below as qPCR1 – qPCR5, ddPCR6, qPCR7 – qPCR11, respectively.

## Sample collection and DNA extraction

A total of 42 DNA samples obtained from blood of naturally BLV-infected dairy cattle from Poland, Moldova, Pakistan, Ukraine, Canada and United States were used for this study. Thirty-six of them were archival DNA samples obtained between 2012–2018 as described in our previous studies on samples from Poland (n=21) [16, 17], Moldova (n=4) [18], Pakistan (n=5) [19] and Ukraine (n=6) [15, 20]. Between 2020–2021 6 peripheral blood and serum samples from naturally BLV-infected cattle were obtained from three dairy farms of Alberta, Canada and two dairy farms of Michigan, US. Serological testing and sample processing were conducted by the laboratories from which the samples originated. The genomic DNA from Canadian and US samples was extracted from whole blood using a Quick DNA Miniprep Plus kit (Zymo Research) and a DNeasy Blood & Tissue Kit (Qiagen), respectively in University of Calgary and Michigan State University and sent to the NVRI in the form of DNA solutions. Additionally, one plasmid DNA sample (pBLV344) was kindly supplied by Luc Willems (University of Liège, Belgium) and DNA extracted from FLK-BLV cells were included as positive controls. Finally, DNA extracted from PBL of a serologically negative cattle was included as negative control. At the NVRI, the DNA concentration in all samples was estimated by spectrophotometry using a NanoPhotometer (Implen). Each sample was divided into eleven identical aliquots containing between 800 and 4,000 ng of lyophilised genomic DNA. Eleven identical sets of these samples were lyophilized (Alpha 1–4 LSC basic, Martin Christ Gefriertrocknungsanlagen GmbH) and distributed to participating laboratories. At the NVRI, all samples were coded (identification [21] run numbers 1 to 44) to perform a blinded testing. The samples, together with instructions for their preparation (Additional file 1), were shipped by air at room temperature (RT).

## Examination of DNA quality/stability

Since different extraction methods and lyophilization process were employed for the preparation of the DNA samples, it was necessary to test the quality of the DNA at the NVRI laboratory. For that purpose, one complete set of samples (n=44) was tested by Fragment Analyzer (Agilent Technologies), before and after freeze-drying, to assess DNA quality by calculating a Genomic Quality Number (GQN) for every sample. Low GQN value (<2.5) represents sheared or degraded DNA. A high GQN (>9) represents undegraded DNA. In addition, quality of DNA was assessed by determination of copy number of

the histone H3 family 3A (*H3F3A*) housekeeping gene using quantitative real-time PCR (qPCR) [22]. The qPCR results were expressed as the number of *H3F3A* gene copies per 300 ng of DNA in each sample. Grubbs' test was performed to determine outliers. To test the stability of DNA, samples were stored for 20 days at RT (10 days) and at +4 °C (10 days) and were retested by Fragment Analyzer and qPCR 21 days later. A Mann–Whitney U-test was used to compare the median values between fresh and stored samples (time 0 and time 1), respectively.

## Description of BLV qPCR protocols used by participating laboratories

All participating laboratories performed their qPCR or ddPCR using a variety of different equipment, reagents, and reaction conditions, which had been set up, validated, and evaluated previously and are currently used as working protocols. The specific features of each of these protocols are described below and summarized in Table 1.

All laboratories applied standard procedures for avoiding false-positive results indicative of DNA contamination, such as the use of separate rooms for preparing reaction mixtures, adding the samples, and performing the amplification reaction. One of the ten BLV qPCRs used LTR region and the remaining nine qPCRs used the *pol* gene as the target sequence for amplification, while the ddPCR amplified the *env* gene.

### Method qPCR1

The BLV qPCR amplifying a 187-bp *pol* gene was performed according to a previously published methods [23, 24]. A real-time fluorescence resonance energy transfer (FRET) PCR was carried out in a 20-µl PCR mixture containing 10 µl handmade reaction master mix and 10 µl genomic DNA. The PCR buffer was 4.5 mM MgCl2, 50 mM KCl, 20 mM Tris–HCl, pH 8.4, supplemented with 0.05% each Tween20 and Non-idet P-40, and 0.03% acetylated BSA (Roche Applied Science). For each 20 µl total reaction volume, the nucleotides were used at 0.2 mM each and 1.5 U Platinum Taq DNA polymerase (Invitrogen, Carlsbad, CA, USA) was used. Primers were used at 1 µM, LCRed640 probe was used at 0.2 µM, and 6-FAM probe was used at 0.1 µM. Amplification was performed in the Roche Light Cycler 480 II (Roche Molecular Biochemicals) using 10 min denaturation step at 95 °C, followed by 18 high-stringency step-down thermal cycles and 30 low-stringency fluorescence acquisition cycles.

A plasmid containing the BLV-PCR amplicon region was diluted ten-fold from 1×10⁵ copies to 10 copies per 10 µl and was used as a standard to measure the BLV copy numbers.
<table><caption>**Table 1** Primers and probe for qPCR and ddPCR used in this study</caption><thead><tr><th>PCR assay</th><th>Institution/ country</th><th>Detection method</th><th>Target</th><th>Primer/probe position</th><th>Sequence (5'-3')</th><th>Amplicon size (bp)</th><th>Reference</th></tr></thead><tbody><tr><td>qPCR1</td><td>AUCVM/US</td><td>FRET qPCR</td><td>pol</td><td>Fw, 2321-2345<br>Rv, 2484-2507<br>Anchor probe,<br>2407-2430<br>Reporter probe,<br>2431-2454</td><td>CCTCAATTCCCTTTA<br>AACTAGAACG<br>ATGGGCTTTGTAAGA<br>GCATTTGTA<br>GACGGGCCAGGCAAT<br>AATCCAGT-6-FAM<br>LCRed640-TTCCCGGTA<br>CGGAAACCAAATGG-<br>phosphate</td><td>187</td><td>[23, 24]</td></tr><tr><td>qPCR2</td><td>CentralStar/US</td><td>Multiplex TaqMan qPCR</td><td>pol,<br>β-actin,<br>spike-in control</td><td>N/A</td><td>N/A</td><td>N/A</td><td>[25]</td></tr><tr><td>qPCR3</td><td>LFDA-MG/ Brasil</td><td>TaqMan qPCR</td><td>LTR</td><td>BLV.LTR.114.Fw, 37-54<br>BLV.LTR.114.Rv, 131-150<br>BLV.LTR.114.Probe,<br>108-130</td><td>GCCCCGTAAACCAGA<br>CAG<br>CTCAGGGTGTGGATT<br>TCTCG<br>FAM-TACCTCCCC/ZEN/<br>AACTTCCCCTTTCC-IABkF</td><td>114</td><td>[26]</td></tr><tr><td>qPCR4</td><td>CIVETAN/Argentina</td><td>SYBRGreen qPCR</td><td>pol</td><td>Fw, 4618-4636<br>Rv, 4657-4676</td><td>CACCATTCACCCCA<br>CTTG<br>TCAGAGCCCTTGGGT<br>GTTTC</td><td>59</td><td>[27]</td></tr><tr><td>qPCR5</td><td>FA-IU/Japan</td><td>Duplex TaqMan qPCR</td><td>pol,<br>RPPH1</td><td>N/A</td><td>N/A</td><td>N/A</td><td>[28]</td></tr><tr><td>ddPCR6</td><td>UdelaR/Urugway</td><td>TaqMan ddPCR</td><td>env</td><td>Fw, 5384-5406<br>Rv, 5525-5546<br>Probe, 5477-5497</td><td>CAGTGACTGGGTTCC<br>CTCTGTC<br>AGGGCGAGRCCGGGT<br>CCAGAG<br>HEX-CCCTCCCTGGGC<br>TCCCGARA-BHQ1</td><td>162</td><td>[29]</td></tr><tr><td>qPCR7</td><td>VEINST/Croatia</td><td>TaqMan qPCR</td><td>pol</td><td>MRBLVL Fw, 2321-2340<br>MRBLVR Rv, 2421-2440<br>MRBLV Probe, 2341-2360</td><td>CCTCAATTCCCTTTA<br>AACTA<br>GTACCGGGAAGACTG<br>GATTA<br>6FAM-GAACGCCTCCAG<br>GCCCTTCA-BHQ1</td><td>120</td><td>[30]</td></tr><tr><td>qPCR8</td><td>INTA/Argentina</td><td>SYBRGreen qPCR</td><td>pol</td><td>BLVpol5f Fw, 2321-2340<br>BLVpol3r Rv, 2421-2440</td><td>CCTCAATTCCCTTTA<br>AACTA<br>GTACCGGGAAGACTG<br>GATTA</td><td>120</td><td>[30, 31]</td></tr><tr><td>qPCR9</td><td>LCV/Spain</td><td>Duplex TaqMan qPCR</td><td>pol<br>β-actin</td><td>MRBLVL Fw, 2321-2340<br>MRBLVR Rv, 2421-2440<br>Probe, 2341-2360<br>ACT-1005-F Fw,<br>1005-1029<br>ACT-1135-R Rv,<br>1135-1114<br>ACT-1081 Probe,<br>1081-1105</td><td>CCTCAATTCCCTTTA<br>AACTA<br>GTACCGGGAAGACTG<br>GATTA<br>Cy5-GAACGCCTCCAG<br>GCCCTTCA-BHQ1<br>CAGCACAATGAAGAT<br>CAAGATCATC<br>CGGACTCATCGTACTCCT<br>GCTT<br>HEX-TCGCTGTCCACC<br>TTCCAGCAGATGT-BHQ1</td><td>120<br>130</td><td>[30, 32]</td></tr><tr><td>qPCR10</td><td>NVRI/Poland</td><td>TaqMan qPCR</td><td>pol</td><td>MRBLVL Fw, 2321-2340<br>MRBLVR Rv, 2421-2440<br>MRBLV Probe, 2341-2360</td><td>CCTCAATTCCCTTTA<br>AACTA<br>GTACCGGGAAGACTG<br>GATTA<br>6FAM-GAACGCCTCCAG<br>GCCCTTCA-BHQ1</td><td>120</td><td>[30]</td></tr><tr><td>qPCR11</td><td>Anses/France</td><td>TaqMan qPCR</td><td>pol</td><td>MRBLVL Fw, 2321-2340<br>MRBLVR Rv, 2421-2440<br>MRBLV Probe, 2341-2360</td><td>CCTCAATTCCCTTTA<br>AACTA<br>GTACCGGGAAGACTG<br>GATTA<br>6FAM-GAACGCCTCCAG<br>GCCCTTCA-BHQ1</td><td>120</td><td>[30]</td></tr></tbody></table>

N/A no data available
### Method qPCR2

A BLV proviral load qPCR assay developed by AntelBio, a division of CentralStar Cooperative Inc. on Applied Biosystems 7500 Real-Time PCR system [25, 33]. This multiplex assay amplifies the BLV *pol* gene along with the bovine *β-actin* gene and an internal amplification control, “Spike”. A quantitative TaqMan PCR was carried out in a 25-µl PCR mixture containing 12.5 µl of 2X Inhibi-Taq Multiplex HotStart qPCR MasterMix (Empirical Bioscience), 16 nM each BLV primer, 16 nM each *β-actin* primer, 8 nM each spike primer, 8 nM BLV FAM-probe, 8 nM *β-actin* Cy5-probe, 4 nM spike JOE-probe, 1 µl of an internal spike-in control (10,000 copies per µl), 7.25 µl of nuclease-free water and 4 µl of DNA sample for each qPCR reaction. The thermal PCR protocol was as follows: 95 °C for 10 min, 40×(95 °C for 15 s, 60 °C for 1 min). Copy numbers of both the BLV *pol* gene and bovine *β-Actin* were derived using a plasmid containing target sequences, quantified by ddPCR, diluted 1×10⁶ copies per µl to 10 copies per µl in tenfold dilutions. DNA concentrations of each sample were measured using a Qubit 4 Fluorometer and used in combination with the qPCR copy numbers to calculate BLV copies per 100 ng.

### Method qPCR3

The qPCR assays for the BLV LTR gene were performed according to a previously published methods [26]. Genomic DNA was amplified by TaqMan PCR with 10 µl of GoTaq Probe qPCR Master Mix×2 (Promega), 0.6 pmol/µl each primer, 0.3 pmol/µl double-quenched probe and 100 ng genomic DNA. Amplification was performed in the CFX96 cycler (BioRad) according to the protocol: 5 min denaturation at 95°C followed by 45 cycles (60 s at 94°C and 60 s at 60°C). The efficiency of each reaction was calculated from the serial dilution of DNA extracted from BLV persistently infected fetal lamb kidney (FLK) cells, starting at a concentration of 100 ng/µl [21]. The detection limit was tested using a plasmid containing the target of the qPCRs, starting at 10³ ng/µl.

### Method qPCR4

The quantitative real-time PCR was done with the primers for the BLV *pol* gene as previously described [34]. The qPCR reaction mix contained 1×PCR Master Mix with SYBR Green (FastStart Universal SYBR Green Master Rox, Roche), 0.3 µM each primer and 30 ng of extracted genomic DNA. Amplification was performed in Quant-Studio 5 Real-Time PCR System (Applied Biosystems) under the following conditions: 2 min at 50 °C, 10 min at 95 °C, 40 cycles of 15 s at 95 °C and 60 s at 60 °C. A standard curve of six tenfold serial dilutions of pBLV, containing 1×10⁶ to 10 BLV copies, was built and run 3 times for

validation of the method. The number of provirus copies per reaction (100 ng) was calculated.

### Method qPCR5

BLV PVLs were determined by using qPCR kit, RC202 (Takara Bio, Shiga, Japan) [28, 35]. This qPCR assay amplifies the BLV *pol* gene along with the bovine *RPPH1* gene as an internal control. Briefly, 100 ng genomic DNA was amplified by TaqMan PCR with four primers for *pol* gene and *RPPH1* gene according to the manufacturer's instructions: 30 s denaturation at 95 °C followed by 45 cycles (5 s at 95 °C and 30 s at 60 °C). The qPCR was performed on a QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific K.K., Tokyo, Japan). Standard curve was generated by creating tenfold serial dilutions of the standard plasmid included in the kit. The standards for calibration ranged from 1 to 10⁶ copies/reaction and were run in duplicate. The number of provirus copies per 100 ng was calculated.

### Method ddPCR6

The digital droplet PCR (ddPCR) assay for the *env* gene of the BLV was performed using the protocol previously described by [28, 29]. An absolute quantification by TaqMan ddPCR was performed in a typical 20-µl assay, 1 µl of DNA sample was mixed with 1 µl of each primer (10 µM), 0.5 µl of probe (10 µM), and 2×Supermix emulsified with oil (Bio-Rad). The droplets were transferred to a 96-well plate (Eppendorf). The PCR assay was performed in a thermocycler (C1000 touch cycler; Bio-Rad) with the following parameters: initial denaturation of 10 min at 95 °C, then 40 cycles of 30 s at 94 °C, and 1 min at 58 °C, with final deactivation of the enzyme for 10 min at 98 °C. The presence of fluorescent droplets determined the number of resulting positive events that were analyzed in the software (QuantaSoft v.1.7.4; Bio-Rad), using dot charts. The number of provirus copies per 100 ng were calculated. Each sample was run in duplicate, and results were averaged.

### Method qPCR7

This qPCR method for the BLV *pol* gene is a modified option of widely available quantitative TaqMan qPCR described by Rola-Łuszczak et al. [11], using the same primers and standards. A quantitative TaqMan PCR was performed in a 20 µl PCR mix containing 10 µl of 2×ORA qPCR Probe ROX L Mix (highQu, Kraichtal, Germany), 2 µl primer/probe mix (final concentration 400 nM of each of the primers, 200 nM of BLV probe), and 3 µl extracted genomic DNA. Amplification was performed in the Rotor-Gene Q system (Qiagen) with an initial denaturation step and polymerase activation at 95 °C
for 3 min, followed by 45 cycles of 95 °C for 5 s and 60 °C for 30 s. As a standard, plasmid pBLV1 (NVRI, Pulawy, PL) containing a BLV *pol* fragment was used. Tenfold dilutions of plasmid DNA were made from 1×10¹⁰ copies to 1×10¹ copies per reaction and used to generate the standard curve and estimate BLV copy number per 100 ng.

### Method qPCR8

Proviral load quantification was assessed by SYBR Green real-time quantitative PCR (qPCR) using the *pol* gene as the target sequence [36]. Briefly, 12-µl PCR mixture contained Fast Start Universal SYBR Green Master Mix (Roche), 800 nM each BLV *pol* primers and 1 µl DNA as template. The reactions were incubated at 50 °C for 2 min and 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s, 55 °C for 15 s and 60 °C for 1 min. All samples were tested in duplicate on a StepOne Plus machine (Applied Biosystems). A positive and negative control, as well as a no-template control, were included in each plate. After the reaction was completed, the specificity of the ampli-cons was checked by analyzing the individual dissociation curves. As a standard, plasmid pBLV1 (NVRI, Pulawy, PL) containing a BLV *pol* fragment was used. Tenfold dilutions of plasmid DNA were made from 1×10⁶ to 10 copies per µl and used to generate the standard curve and estimate BLV copy number per 100 ng.

### Method qPCR9

This qPCR method is a modified option of widely available quantitative TaqMan qPCR described by Rola-Łuszczak et al. [11], using the same primers and standards. The detection of BLV genome was combined with an endogenous control system (Toussaint 2007) in a duplex assay. Briefly, 20-µl qPCR reaction contained AhPath ID™ One-Step RT-PCR Reagents with ROX (Applied Biosystems, CA, USA) – 10 µl of 2×RT-PCR buffer and 0.8 µl of 25×RT-PCR enzyme mix, 400 nM each primer for *pol* gene, 100 nM BLV specific probe, 40 nM each β-actin primer, 40 nM β-actin specific probe and 2 µl DNA sample. All samples were tested in ABI7500 Real-Time PCR System (Applied Biosystems) according to the following protocol: 10 min at 48 °C (reverse transcription), 10 min at 95 °C (inactivation reverse transcriptase / activation Taq polymerase) followed by 45 cycles (15 s at 95 °C and 60 s at 60 °C). As a standard, plasmid pBLV1 (NVRI, Pulawy, PL) containing a BLV *pol* fragment was used. Tenfold dilutions of plasmid DNA were made from 1×10⁴ copies to 0.1 copies per µl and used to generate the standard curve and estimate BLV copy number per 100 ng.

### Method qPCR10

The BLV qPCR was performed as published previously [11]. A quantitative TaqMan PCR was carried out in a 25-µl PCR mixture containing 12.5 µl of 2×QuantiTect Multiplex PCR NoROX master mix (Qiagen), 0.4 µM each primer, 0.2 µM specific BLV probe, and 500 ng of extracted genomic DNA. Amplification was performed in the Rotor-Gene Q system (Qiagen) using an initial denaturation step and polymerase activation at 95 °C for 15 min, followed by 50 cycles of 94 °C for 60 s and 60 °C for 60 s. All samples were amplified in duplicate. As a standard, the pBLV1 plasmid (NVRI, Pulawy, PL), containing a 120-bp BLV *pol* fragment, was used. Tenfold dilutions of this standard were made from 1×10⁶ copies per µl to 100 copies per µl and were used to estimate the BLV copy numbers per 100 ng.

### Method qPCR11

This qPCR method for the BLV *pol* gene is a modified option of widely available quantitative TaqMan qPCR described by Rola-Łuszczak et al. [11], using the same primers and standards. The reaction mixture contained 400 nM of each primer, 200 nM of probe, 10 µl of 2×Sso-Fast probes supermix (Bio-Rad), 5 µl of DNA sample and H₂O up to 20 µl of the final volume. PCR assays were carried out on a CFX96 thermocycler (Bio-Rad) under the following amplification profile: 98 °C for 3 min, followed by 45 cycles of 95 °C for 5 s and 60 °C for 30 s. As a standard, plasmid pBLV1 (NVRI, Pulawy, PL) containing a BLV *pol* fragment was used. Tenfold dilutions of plasmid DNA were used to generate the standard curve and estimate BLV copy number per 100 ng.

## Analysis of BLV *pol*, *env* and LTR sequences targeted by particular qPCR/ddPCR assays

In order to assess full-length *pol*, *env* and LTR sequence variability among BLV genotypes, all BLV sequences (n=2191) available on 30 September 2023 in GenBank (https://www.ncbi.nlm.nih.gov/GenBank/) repository were retrieved. From the collected sequences, 100 *pol*, *env* and LTR sequences, which were characterized by the highest level of sequence variability and divergence, were selected for the further analysis. A *pol*-based, *env*-based and LTR-based maximum likelihood (ML) phylogenetic trees (see Additional file 6) was constructed to assign genotypes to the unassigned BLV genomes [37–39]. For all genes and LTR region the Tamura-Nei model and Bootstrap replications (1,000) were applied. In this analysis, *pol* sequences were assigned to 7 BLV genotypes (G1, G2, G3, G4, G6, G9, and G10), while *env* and LTR sequences were assigned to 10 BLV genotypes (G1, G2, G3, G4, G5,
G6, G7, G8, G9, and G10). Phylogeny of the same isolates assigned to particular genotypes by ML method was confirmed by Mr. Bayes analysis [40–42] (data not shown). From this analysis, a total of 100 full-length *pol*, *env* and LTR sequences were used for multiple-sequence alignment (MSA) using ClustalW algorithm, implemented in MEGA X. For all sequences, nucleotide diversity (π), defined as the average number of nucleotide differences per site between two DNA sequences in all possible pairs in the sample population, was estimated using MEGA X. To measure the relative variation in different positions of aligned genes and LTR region the Shannon's entropy (a quantitative measure of diversity in the alignment, where H=0 indicates complete conservation) was estimated using BioEdit v. 7.2.5 software 64. The statistical analyses were performed using DATAtab e.U. Graz, Austria and GraphPad Software by Dotmatics, Boston.

## Results

### Examination of the quality and stability of DNA samples

To test the quality of DNA samples, the *H3F3A* copy number of each individual sample was assessed by qPCR at the NVRI. Copy numbers were normalized to DNA mass input and results were expressed as copy numbers per 300 ng of total DNA. The respective values were tested by Grubbs' test. The results for 43 DNA samples (sample ID: 42 with BLV genome plasmid was excluded) followed a normal distribution (Shapiro-Wilk 0.97; *P*=0.286), with a mean value of 35,626 copies (95% confidence interval [43] 33,843 to 37,408 copies), a minimum value of 19,848 copies and a maximum value of 46,951 copies (see Additional file 2). Despite a low value for sample ID: 40 no significant outlier was detected in the dataset (*P*>0.05). Therefore, it can be assumed that the DNA quality was acceptable for all samples present in the panel. Next, DNA stability was assessed by retesting the *H3F3A* copy numbers in each sample (n=43) after a combined storage consisting in 10 days at RT and 10 days at +4°C. A Mann-Whitney U-test was used to compare the median values between fresh and stored samples (time 0 and time 1, respectively), and no significant difference was observed at the 5% level (*P*=0.187) (Fig. 1A).

In addition, the quality of DNA samples after lyophilization was analyzed. DNA from individual samples (n=43) was assessed with the genomic DNA quality number on the Fragment Analyzer system. The GQN from all lyophilized samples ranged from 4.0 to 9.7—that represented undegraded DNA. There was no significant difference in GQN values between `before freeze-drying` and `after freeze-drying` groups with respect to the corresponding DNA samples (*P*=0.236) (Fig. 1B). Altogether, these results suggested that sample storage, lyophilization and shipping has a minimal impact in DNA

stability and further testing during the interlaboratory trial.

### Detection of BLV proviral DNA by different qPCR assays

A total of 44 DNA samples, including two positive (ID: 42 and 43) and one negative (ID: 32) controls, were blinded and independently tested by eleven laboratories using their own qPCR methods (Table 2). All laboratories measured the concentration of DNA in samples (Additional file 3). BLV provirus copy number was normalized to DNA concentration and expressed per 100 ng of genomic DNA for each test.

Except for the positive (pBLV344 and FLK cell line) and the negative controls, all samples had previously shown detectable levels of BLV-specific antibodies (BLV-Abs) by enzyme-linked immunosorbent assays (ELISA). During the current interlaboratory study, both the positive and negative controls were assessed adequately by all eleven PCR tests. Of all 43 positive samples, 43, 35, 37, 36, 40, 32, 40, 42, 42, 42 and 41 samples were detected as positive by the qPCR1, qPCR2, qPCR3, qPCR4, qPCR5, ddPCR6, qPCR7, qPCR8, qPCR9, qPCR10 and qPCR11 methods, respectively. Based on these observations, the most sensitive method was the qPCR1, and the method with the lowest sensitivity was the ddPCR6. Twenty-nine out of 44 samples were identified correctly by all qPCRs. The remaining 15 samples gave discordant results. Comparison of qualitative results (positive versus negative) from all eleven methods revealed 87.33% overall agreement and a kappa value of 0.396 (Cohen's kappa method adapted by Fleiss) [44, 45]. The levels of agreement among the results from the eleven methods are represented in Table 3. The maximum agreement was seen between two methods (qPCR9 and qPCR10 [100% agreement and a Cohen's kappa value of 1.000]) that used similar protocols and targeted the same region of BLV *pol*.

### Analysis of BLV *pol*, *env* and LTR sequences targeted by particular PCR assays

Due to differences in performance observed among the *pol*-based qPCR assays (the qPCR1, qPCR2, qPCR4, qPCR5 and qPCR7- qPCR11 methods), and considering that the *env*-based ddPCR6 and LTR-based qPCR3 assay showed the lowest sensitivity and the poorest agreement with the other assays, the degree of sequence variability between the *pol*, *env* and LTR genes was addressed. From the MSAs for *pol*, *env* and LTR, the nucleotide diversity (π) was calculated. The π value for *pol* gene was lower than that for LTR and *env* gene (π<sub>pol</sub>, 0.023 [standard deviation {SD}, 0.018]; π<sub>LTR</sub>, 0.024 [SD, 0.011]; π<sub>env</sub>, 0.037 [SD, 0.013]). From this analysis, *pol* sequences appeared to be less variable than *env* and LTR sequences. In addition, we performed a Shannon
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Assessment of the stability of DNA samples. A Shown are copy numbers of the H3F3A housekeeping gene in 43 DNA samples that were stored in 10 days at RT and 10 days at +4°C and tested twice with a 21-day interval. A Mann-Whitney U-test was used to compare the median values between two groups (*P* = 0.187); B Shown are GQN values (n = 43) tested twice with a 21-day interval: 'before freeze-drying' and 'after freeze-drying'. A Mann-Whitney U-test results between two groups (*P* = 0.236)</figcaption></figure>

entropy-based per-site variability profile of the *pol*, *env* and LTR sequences used in this study (Fig. 2A-C).

The all-observed entropy plots were homogeneous along the whole sequences. Considering the three regions of *pol* gene, the highest entropy (4.67) occurred in the region targeted by the qPCR1 primers, whereas the entropy for qPCR7-qPCR11 and qPCR4 primers were 1.57 and 0.38, respectively. For the LTR region targeted by qPCR3 primers and for *env* gene targeted by ddPCR6, the total entropy was equal to 4.46 and 7.85, respectively. This analysis showed a marked region of variability for LTR and *env* fragments. Interestingly, we noted that the qPCR7-qPCR11 targeted the most conserved regions of reverse transcriptase and qPCR4

primers targeted the most-conserved region of virus integrase (Fig. 2A-C; see also Additional file 7).

## Quantitation of BLV proviral DNA by different qPCR/ddPCR assays

To analyze whether the range of copy numbers detected by each qPCR was comparable to those of the others, Kruskal-Wallis one-way analysis of variance (ANOVA) was used. The violin plots were used to visualize the ANOVA results (Fig. 3A-B).

The grouping variable revealed significant differences among the distributions of proviral DNA copy numbers with the various qPCRs (*P* < 0.001). These results showed
<table><caption>**Table 2** Origin of the samples included in the study (n=44) and the respective results. Copy numbers were then normalized to DNA concentration and expressed as per 100 ng of genomic DNA</caption><thead><tr><th>Sample No.</th><th>Origin</th><th>Genotype</th><th>ELISA</th><th>qPCR result (no. of BLV DNA copies per 100 ng of total DNA)<sup>a</sup></th><th>qPCR1 (pol)</th><th>qPCR2 (pol)</th><th>qPCR3 (LTR)</th><th>qPCR4 (pol)</th><th>qPCR5 (pol)</th><th>ddPCR6 (env)</th><th>qPCR7 (pol)</th><th>qPCR8 (pol)</th><th>qPCR9 (pol)</th><th>qPCR10 (pol)</th><th>qPCR11 (pol)</th></tr></thead><tbody><tr><td>1</td><td>Poland</td><td>G4</td><td>+</td><td>17253.9</td><td>7820.2</td><td>962541.0</td><td>285583.3</td><td>3490.2</td><td>4480.0</td><td>36700.0</td><td>10820.0</td><td>24700.0</td><td>11613.8</td><td>84463.7</td></tr><tr><td>2</td><td>Poland</td><td>G4</td><td>+</td><td>2042.6</td><td>465.6</td><td>48555.9</td><td>19950.0</td><td>947.4</td><td>442.0</td><td>159.0</td><td>489.0</td><td>209.0</td><td>835.8</td><td>20127.2</td></tr><tr><td>3</td><td>Ukraine</td><td>G4</td><td>+</td><td>30365.9</td><td>8329.2</td><td>184258.8</td><td>398733.3</td><td>19476.0</td><td>8320.0</td><td>50800.0</td><td>28699.0</td><td>16400.0</td><td>13441.1</td><td>165171.6</td></tr><tr><td>4</td><td>Ukraine</td><td>G8</td><td>+</td><td>6994.0</td><td>4818.0</td><td>14457.6</td><td>120130.0</td><td>7180.4</td><td>4.4</td><td>5710.0</td><td>5459.0</td><td>3900.0</td><td>3384.3</td><td>52440.6</td></tr><tr><td>5</td><td>Ukraine</td><td>G8</td><td>+</td><td>6653.5</td><td>6354.8</td><td>64721.9</td><td>109896.7</td><td>4753.6</td><td>1734.0</td><td>7110.0</td><td>9199.0</td><td>1660.0</td><td>46.3</td><td>249756.1</td></tr><tr><td>6</td><td>Ukraine</td><td><sup>b</sup>Ns</td><td>+</td><td>11586.5</td><td>10479.8</td><td>191854.7</td><td>217700.0</td><td>11272.9</td><td>1501.2</td><td>4970.0</td><td>12250.0</td><td>15500.0</td><td>9476.4</td><td>48900.5</td></tr><tr><td>7</td><td>Ukraine</td><td>Ns</td><td>+</td><td>89.1</td><td>52.3</td><td>43.3</td><td>5613.3</td><td>7.7</td><td>10.4</td><td>12.6</td><td>51.0</td><td>1.7</td><td>27.0</td><td>502.7</td></tr><tr><td>8</td><td>Poland</td><td>G8</td><td>+</td><td>1146.9</td><td>381.9</td><td>119.5</td><td>8993.3</td><td>479.7</td><td>NEG</td><td>173.0</td><td>267.0</td><td>118.0</td><td>388.9</td><td>2738.9</td></tr><tr><td>9</td><td>Poland</td><td>G7</td><td>+</td><td>0.3</td><td>NEG</td><td>NEG</td><td>NEG</td><td>1.0</td><td>NEG</td><td>12.3</td><td>17.0</td><td>NEG</td><td>NEG</td><td>6.1</td></tr><tr><td>10</td><td>Moldova</td><td>G7</td><td>+</td><td>130.5</td><td>1345.1</td><td>2533.8</td><td>7090.0</td><td>21.1</td><td>76.0</td><td>188.0</td><td>224.0</td><td>3.4</td><td>207.3</td><td>2935.0</td></tr><tr><td>11</td><td>Ukraine</td><td>Ns</td><td>+</td><td>1618.1</td><td>1274.6</td><td>308.8</td><td>18996.7</td><td>934.5</td><td>NEG</td><td>1480.0</td><td>948.0</td><td>209.0</td><td>816.1</td><td>8483.0</td></tr><tr><td>12</td><td>Moldova</td><td>G7</td><td>+</td><td>5964.5</td><td>6473.7</td><td>127613.4</td><td>235810.0</td><td>1673.2</td><td>5320.0</td><td>21600.0</td><td>11541.0</td><td>5950.0</td><td>16235.6</td><td>171017.9</td></tr><tr><td>13</td><td>Poland</td><td>G4</td><td>+</td><td>41743.1</td><td>2816.6</td><td>580303.9</td><td>378866.7</td><td>15555.0</td><td>6100.0</td><td>21300.0</td><td>16300.0</td><td>18400.0</td><td>28410.3</td><td>505675.6</td></tr><tr><td>14</td><td>Poland</td><td>Ns</td><td>+</td><td>2.0</td><td>NEG</td><td>NEG</td><td>NEG</td><td>2.5</td><td>NEG</td><td>11.0</td><td>4.0</td><td>0.1</td><td>2.1</td><td>97.9</td></tr><tr><td>15</td><td>Poland</td><td>G8</td><td>+</td><td>13.0</td><td>NEG</td><td>NEG</td><td>NEG</td><td>6.4</td><td>NEG</td><td>NEG</td><td>12.0</td><td>0.1</td><td>3.9</td><td>170.7</td></tr><tr><td>16</td><td>Moldova</td><td>Ns</td><td>+</td><td>6.8</td><td>NEG</td><td>14.1</td><td>NEG</td><td>4.0</td><td>4.6</td><td>40.8</td><td>6.0</td><td>0.1</td><td>1.8</td><td>97.8</td></tr><tr><td>17</td><td>Poland</td><td>G4</td><td>+</td><td>40.2</td><td>69.6</td><td>1002.7</td><td>736.7</td><td>11.3</td><td>12.0</td><td>18.8</td><td>18.0</td><td>1.2</td><td>19.8</td><td>368.3</td></tr><tr><td>18</td><td>Poland</td><td>G7</td><td>+</td><td>1.5</td><td>NEG</td><td>306.5</td><td>536.7</td><td>NEG</td><td>3.8</td><td>23.7</td><td>NEG</td><td>0.1</td><td>2.5</td><td>35.0</td></tr><tr><td>19</td><td>Poland</td><td>G7</td><td>+</td><td>18092.8</td><td>807.5</td><td>246292.1</td><td>432900.0</td><td>3215.4</td><td>7760.0</td><td>47700.0</td><td>16849.0</td><td>19800.0</td><td>18718.6</td><td>218610.3</td></tr><tr><td>20</td><td>Moldova</td><td>G7</td><td>+</td><td>2188.4</td><td>19195.8</td><td>268180.6</td><td>360833.3</td><td>2957.0</td><td>5700.0</td><td>29100.0</td><td>3863.0</td><td>985.0</td><td>19494.0</td><td>129483.2</td></tr><tr><td>21</td><td>Poland</td><td>G4</td><td>+</td><td>423.9</td><td>1593.3</td><td>295.7</td><td>8503.3</td><td>202.7</td><td>214.0</td><td>460.0</td><td>267.0</td><td>7.4</td><td>277.0</td><td>5211.4</td></tr><tr><td>22</td><td>Poland</td><td>G4</td><td>+</td><td>33005.6</td><td>14517.9</td><td>335322.0</td><td>364700.0</td><td>13565.7</td><td>12040.0</td><td>15400.0</td><td>4539.0</td><td>12200.0</td><td>28424.2</td><td>152574.9</td></tr><tr><td>23</td><td>Poland</td><td>G4</td><td>+</td><td>4776.8</td><td>4750.2</td><td>197.2</td><td>80110.0</td><td>3385.1</td><td>1800.0</td><td>13200.0</td><td>2050.0</td><td>1280.0</td><td>2269.0</td><td>71669.1</td></tr><tr><td>24</td><td>Poland</td><td>G7</td><td>+</td><td>4.9</td><td>NEG</td><td>81.6</td><td>NEG</td><td>1.9</td><td>8.2</td><td>34.1</td><td>5.0</td><td>0.1</td><td>1.4</td><td>99.2</td></tr><tr><td>25</td><td>Poland</td><td>G8</td><td>+</td><td>38.1</td><td>1152.2</td><td>8.0</td><td>533.3</td><td>30.0</td><td>NEG</td><td>144.0</td><td>26.0</td><td>0.2</td><td>33.8</td><td>524.7</td></tr><tr><td>26</td><td>Poland</td><td>G8</td><td>+</td><td>1017.4</td><td>921.7</td><td>233.9</td><td>19566.7</td><td>828.3</td><td>NEG</td><td>1850.0</td><td>850.0</td><td>74.4</td><td>781.9</td><td>9819.5</td></tr><tr><td>27</td><td>Poland</td><td>G4</td><td>+</td><td>41388.9</td><td>9513.7</td><td>NEG</td><td>218583.3</td><td>10410.8</td><td>4000.0</td><td>30800.0</td><td>6033.0</td><td>10900.0</td><td>18878.0</td><td>131461.4</td></tr><tr><td>28</td><td>Poland</td><td>G7</td><td>+</td><td>1.0</td><td>NEG</td><td>480.1</td><td>NEG</td><td>NEG</td><td>NEG</td><td>17.5</td><td>9.0</td><td>0.3</td><td>0.7</td><td>NEG</td></tr><tr><td>29</td><td>Poland</td><td>G4</td><td>+</td><td>23913.0</td><td>1886.5</td><td>173311.8</td><td>334533.3</td><td>16739.7</td><td>7570.0</td><td>12700.0</td><td>2476.0</td><td>7590.0</td><td>27163.2</td><td>170384.4</td></tr><tr><td>30</td><td>Poland</td><td>G8</td><td>+</td><td>7.8</td><td>113.5</td><td>NEG</td><td>326.7</td><td>4.8</td><td>NEG</td><td>6.3</td><td>12.0</td><td>1.0</td><td>3.8</td><td>41.5</td></tr><tr><td>31</td><td>Pakistan</td><td>G6</td><td>+</td><td>8920.5</td><td>3293.8</td><td>13751.5</td><td>57813.3</td><td>1633.2</td><td>860.0</td><td>3570.0</td><td>2386.0</td><td>825.0</td><td>3197.8</td><td>49053.1</td></tr></tbody></table>
<table><caption>**Table 2 (continued)**</caption><thead><tr><th rowspan="2">Sample No.</th><th rowspan="2">Origin</th><th rowspan="2">Genotype</th><th rowspan="2">ELISA</th><th colspan="8">qPCR result (no. of BLV DNA copies per 100 ng of total DNA)<sup>a</sup></th></tr><tr><th>qPCR1<sup>a</sup><br>(pol)</th><th>qPCR2<br>(pol)</th><th>qPCR3<br>(LTR)</th><th>qPCR4<br>(pol)</th><th>qPCR5<br>(pol)</th><th>ddPCR6<br>(env)</th><th>qPCR7<sup>d</sup><br>(pol)</th><th>qPCR8<sup>d</sup><br>(pol)</th><th>qPCR9<sup>d</sup><br>(pol)</th><th>qPCR10<sup>d</sup><br>(pol)</th><th>qPCR11<sup>d</sup><br>(pol)</th></tr></thead><tbody><tr><td>32</td><td>Poland</td><td><sup>c</sup>Na</td><td>-</td><td><sup>d</sup>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td></tr><tr><td>33</td><td>Pakistan</td><td>G6</td><td>+</td><td>1308.8</td><td>417.5</td><td>29221.6</td><td>15846.7</td><td>387.0</td><td>351.8</td><td>1490.0</td><td>800.0</td><td>103.0</td><td>856.7</td><td>18479.7</td></tr><tr><td>34</td><td>Pakistan</td><td>G1</td><td>+</td><td>1.6</td><td>NEG</td><td>259.8</td><td>NEG</td><td>NEG</td><td>NEG</td><td>NEG</td><td>8.0</td><td>0.1</td><td>0.7</td><td>NEG</td></tr><tr><td>35</td><td>Pakistan</td><td>G1</td><td>+</td><td>7000.0</td><td>1082.7</td><td>104381.2</td><td>76063.3</td><td>5485.9</td><td>1700.0</td><td>9480.0</td><td>4156.0</td><td>1200.0</td><td>3869.7</td><td>23367.2</td></tr><tr><td>36</td><td>Pakistan</td><td>G6</td><td>+</td><td>326.3</td><td>52.8</td><td>2742.6</td><td>5933.3</td><td>240.3</td><td>176.0</td><td>187.0</td><td>209.0</td><td>16.2</td><td>188.2</td><td>1499.4</td></tr><tr><td>37</td><td>Canada</td><td>G1</td><td>+</td><td>385.7</td><td>3083.5</td><td>437.2</td><td>9630.0</td><td>245.0</td><td>132.0</td><td>430.0</td><td>326.0</td><td>26.5</td><td>169.9</td><td>4185.5</td></tr><tr><td>38</td><td>Canada</td><td>G1</td><td>+</td><td>1452.8</td><td>84.4</td><td>78032.9</td><td>24483.3</td><td>911.6</td><td>338.0</td><td>2000.0</td><td>949.0</td><td>133.0</td><td>591.4</td><td>16415.2</td></tr><tr><td>39</td><td>Canada</td><td>G1</td><td>+</td><td>1245.8</td><td>871.8</td><td>6963.9</td><td>19863.3</td><td>667.6</td><td>268.0</td><td>2140.0</td><td>989.0</td><td>111.0</td><td>410.5</td><td>15095.0</td></tr><tr><td>40</td><td>USA</td><td>G1</td><td>+</td><td>30652.2</td><td>4620.7</td><td>406854.9</td><td>269223.3</td><td>10184.5</td><td>4040.0</td><td>25900.0</td><td>29479.0</td><td>2660.0</td><td>9719.0</td><td>70397.3</td></tr><tr><td>41</td><td>USA</td><td>G1</td><td>+</td><td>51885.2</td><td>23907.4</td><td>607640.5</td><td>179676.7</td><td>22906.4</td><td>5280.0</td><td>32400.0</td><td>51006.0</td><td>6980.0</td><td>16097.0</td><td>458976.1</td></tr><tr><td>42</td><td>Belgium</td><td>G4</td><td>pBLV344<sup>e</sup></td><td>4.07E+10</td><td>4.15E+09</td><td>4E+09</td><td>7E+06</td><td>1.45E+09</td><td>2.9E+07</td><td>7.48E+09</td><td>2.27E+08</td><td>2.55E+11</td><td>2.59E+09</td><td>1.22E+11</td></tr><tr><td>43</td><td>USA</td><td>G1</td><td>BLV-FLK<sup>f</sup></td><td>307692.3</td><td>162180.6</td><td>863440.0</td><td>3250666.7</td><td>75773.5</td><td>28480.0</td><td>168000.0</td><td>244769.0</td><td>482000.0</td><td>138382.8</td><td>4286730.3</td></tr><tr><td>44</td><td>Poland</td><td>G8</td><td>+</td><td>1.7</td><td>41.8</td><td>NEG</td><td>763.3</td><td>10.4</td><td>NEG</td><td>NEG</td><td>5.0</td><td>0.2</td><td>2.0</td><td>5.9</td></tr><tr><td>Total</td><td></td><td>POSITIVE</td><td>43</td><td>43</td><td>35</td><td>37</td><td>36</td><td>40</td><td>32</td><td>40</td><td>42</td><td>42</td><td>42</td><td>41</td></tr><tr><td></td><td></td><td>NEGATIVE</td><td>1</td><td>1</td><td>9</td><td>7</td><td>8</td><td>4</td><td>12</td><td>4</td><td>2</td><td>2</td><td>2</td><td>3</td></tr></tbody></table>

<footnote><sup>a</sup> NEG negative result, <sup>b</sup>Na not sequenced, <sup>c</sup>Na not applicable, <sup>d</sup>in the qPCR two primer sets for pol gene MRBLVL and MRBLVR were used, <sup>e</sup>pure plasmid DNA BLV344 (gift from Luc Willems), <sup>f</sup>DNA isolated from Fetal Lamb Kidney cell line persistently infected with BLV (FLK-BLV), qPCR1 - Auburn University College of Veterinary Medicine, US, qPCR2 CentralStar Cooperative, Michigan, US, qPCR3 - Laboratório Federal de Defesa Agropecuária de Minas Gerais (LFDA-MG), Brasil; qPCR4 - Centro de Investigación Veterinaria de Tandil (CIVETAN), Argentina; qPCR5 - Faculty of Agriculture Iwate University, Japan; ddPCR6 - Universidad de la República de Uruguay (UdelaR), Montevideo, Uruguay; qPCR7 - Croatian Veterinary Institute, Zagreb, Croatia; qPCR8 - Instituto Nacional de Tecnología Agropecuaria (INTA), Argentina; qPCR9 - Laboratorio Central de Veterinaria LCV Algete, Spain; qPCR10 - National Veterinary Institute (NVRI), Pulawy, Poland; qPCR11 - French Agency for Food, Environmental and Occupational Health and Safety (Anses), Niort, France</footnote>
<table><caption>**Table 3** Assessment of qPCR results by pairwise comparison</caption><thead><tr><th rowspan="2">PCR assay</th><th colspan="10">Kappa value or agreement<sup>a</sup> with the results of the following PCR assay</th></tr><tr><th>qPCR1</th><th>qPCR2</th><th>qPCR3</th><th>qPCR4</th><th>qPCR5</th><th>ddPCR6</th><th>qPCR7</th><th>qPCR8</th><th>qPCR9</th><th>qPCR10</th><th>qPCR11</th></tr></thead><tbody><tr><td>qPCR1</td><td>-</td><td>0.166</td><td>0.219</td><td>0.189</td><td>0.377</td><td>0.117</td><td>0.377</td><td>0.656</td><td>0.656</td><td>0.656</td><td>0.482</td></tr><tr><td>qPCR2</td><td>81.8%</td><td>-</td><td>0.391</td><td>0.927</td><td>0.560</td><td>0.441</td><td>0.384</td><td>0.313</td><td>0.313</td><td>0.312</td><td>0.443</td></tr><tr><td>qPCR3</td><td>86.4%</td><td>81.8%</td><td>-</td><td>0.438</td><td>0.075</td><td>0.539</td><td>0.486</td><td>0.163</td><td>0.402</td><td>0.402</td><td>0.115</td></tr><tr><td>qPCR4</td><td>84.1%</td><td>97.7%</td><td>84.1%</td><td>-</td><td>0.431</td><td>0.488</td><td>0.431</td><td>0.137</td><td>0.353</td><td>0.352</td><td>0.495</td></tr><tr><td>qPCR5</td><td>93.2%</td><td>88.6%</td><td>79.5%</td><td>86.4%</td><td>-</td><td>0.276</td><td>0.450</td><td>0.645</td><td>0.290</td><td>0.290</td><td>0.845</td></tr><tr><td>ddPCR6</td><td>75.0%</td><td>79.6%</td><td>84.1%</td><td>81.8%</td><td>77.3%</td><td>-</td><td>0.421</td><td>0.070</td><td>0.225</td><td>0.225</td><td>0.326</td></tr><tr><td>qPCR7</td><td>93.2%</td><td>84.1%</td><td>88.6%</td><td>86.4%</td><td>90.9%</td><td>81.8%</td><td>-</td><td>0.290</td><td>0.290</td><td>0.290</td><td>0.535</td></tr><tr><td>qPCR8</td><td>97.7%</td><td>84.1%</td><td>84.1%</td><td>66.8%</td><td>95.5%</td><td>72.7%</td><td>90.9%</td><td>-</td><td>0.476</td><td>0.476</td><td>0.365</td></tr><tr><td>qPCR9</td><td>97.7%</td><td>84.1%</td><td>88.6%</td><td>86.4%</td><td>90.9%</td><td>77.3%</td><td>90.9%</td><td>95.5%</td><td>-</td><td>1.00</td><td>0.365</td></tr><tr><td>qPCR10</td><td>97.7%</td><td>84.1%</td><td>88.6%</td><td>86.4%</td><td>90.9%</td><td>77.3%</td><td>90.9%</td><td>95.5%</td><td>100%</td><td>-</td><td>0.365</td></tr><tr><td>qPCR11</td><td>95.5%</td><td>86.4%</td><td>81.8%</td><td>88.6%</td><td>97.7%</td><td>79.6%</td><td>93.2%</td><td>93.2%</td><td>93.2%</td><td>93.2%</td><td>-</td></tr></tbody></table>

that the abilities of qPCRs/ddPCR to determine the proviral DNA copy number differed. A Dunn-Bonferroni test was used to compare the groups in pairs to find out which was significantly different. The Dunn-Bonferroni test revealed that the pairwise group comparisons of qPCR2-qPCR4, qPCR3-ddPCR6, qPCR4-qPCR5, qPCR4-ddPCR6, qPCR4-qPCR9, qPCR4-qPCR10, qPCR5-qPCR11, ddPCR6-qPCR11 and qPCR9-qPCR11 have an adjusted P value less than 0.05 and thus, it can be assumed that these groups were significantly different in each pair (see Additional file 4). The Pareto chart was used to show the average copy number values of all methods in descending order. These Pareto charts were prepared based on 80-20 rule, which states that 80% of effects come from 20% of the various causes [46]. The methods that generated the highest copy numbers was qPCR3 and qPCR4, on the other hand the lowest copy numbers and/or highest negative results were generated by ddPCR6 (Fig. 4).

The correlations between copy numbers detected by different qPCRs and ddPCR assays were calculated. The Kendall's Tau correlation coefficient measured between each pair of the assays was shown in the Additional file 5 and in Fig. 5 as a correlation heatmap. The average correlation for all qPCRs and ddPCR assays was strong (Kendall's tau = 0.748; *P* < 0.001).

Since the differences between PCR tests may be influenced by the number of BLV proviral copies present in each sample, we compared the average number of BLV copies between a group of genomic DNA samples that gave concordant results (group I [n=28]) and a group that gave discordant results (group II [n=15]). The mean number of copies was 73,907 (minimum, 0; maximum, 4,286,730) in group I, and 3,479 (minimum, 0; maximum,

218,583) in group II, and this difference was statistically significant (*P* < 0.001 by a Mann-Whitney U-test) (Fig. 6).

The results show that the concordant results group had considerably higher copy numbers (median, 5,549.0) than the discordant results group (median, 6.3).

## Discussion

BLV control and eradication programs consist of correct identification and subsequent segregation/elimination of BLV-infected animals [47]. Detection of BLV-infected cows by testing for BLV-specific antibodies in serum by agar gel immunodiffusion and ELISA is the key step and standard to be implemented of EBL eradication programs according to WOAH (https://www.woah.org/en/disease/enzootic-bovine-leukosis/) [9]. Despite the low cost and high throughput of serological tests, there are several scenarios where highly specific and sensitive molecular assays for the detection of BLV DNA might improve detection and program efficiency.

In this perspective, qPCR assays can detect small quantities of proviral DNA during acute infection, in which animals show very low levels of anti-BLV antibodies [43, 48–50]. qPCR methods can also work as confirmatory tests to clarify ambiguous and inconsistent serological test results [11]. Such quantitative features of qPCRs are crucial when eradication programs progress and prevalence decreases. Moreover, qPCR allows not only the detection of BLV infection but also estimation of the BLV PVL, which directly correlates with the risk of disease transmission [51, 52]. This feature of qPCR allows for a rational segregation of animals based on the stratified risk of transmission. These considerations allow for greater precision in the management of BLV within large

<footnote><sup>a</sup> Values above the blank cells are kappa values; values below the blank cells are percentages of agreement</footnote>
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 2** Sequence variability measured as per-site entropy. A Multiple alignment of the *pol* gene showing the locations of qPCR fragments in regions of the *pol* gene for the qPCR1 (highlighted in pink), qPCR4 (highlighted in yellow) and for the qPCR7, qPCR8, qPCR9, qPCR10 and qPCR11 assays (highlighted in orange). B Multiple alignment of the *env* gene targeted by ddPCR6 (highlighted by blue rectangle). C Multiple alignment of the LTR region by qPCR3 (highlighted in mint)</figcaption></figure>
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 3** Comparison of detection of BLV proviral DNA copy numbers by eleven testing methods. Shown is a box plot of data from Kruskal-Wallis ANOVA, a rank test. The DNA copy numbers for 41 samples, determined independently by each of the 11 qPCRs, were used for the variance analysis. In this analysis, the positive controls (sample ID 42 and ID 43) and negative control (sample ID 32) were excluded. A Violin plot for graphical presentation of the ANOVA of proviral copy number values. B Violin plot for ANOVA analysis of variance, copy number values are presented on a logarithmic scale (Log1.2) for better illustration of copy number differences between PCR methods</figcaption></figure>

herds with a high prevalence of BLV ELISA-positive animals to effectively reduce herd prevalence [13, 53]. BLV is a global burden and the lack of technical standardization of molecular detection systems remains a huge obstacle to compare surveillance data globally based on the first interlaboratory trial performed in 2018 [15]. In the 2018 study we observed an adjusted level of agreement of 70% comparing qualitative qPCR results; however,

inconsistencies amongst methods were larger when low number of copies of BLV DNA were compared. Samples with low copies of BLV DNA (<20 copies per 100 ng) accounted for the higher variability and discrepancies amongst tests. We concluded from the first interlaboratory trial that standardizing protocols to improve sensitivity of assays with lower detection rates was necessary.
<figure><img src="image_9.png" /><figcaption>**Fig. 4** A Pareto chart with the proviral BLV copy mean values for eleven PCR assay arranged in descending order. Pareto charts were prepared based on 80–20 rule, which states that 80% of effects come from 20% of the various causes</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 5** The heatmap of Kendall's tau correlation coefficients between copy numbers detected by ten qPCRs and one ddPCR. Statistically significant differences in the distribution of copy numbers, a moderate, strong and very strong correlation between particular qPCRs/ddPCR was observed. The strength of the association, for absolute values of r, 0–0.19 is regarded as very weak, 0.2–0.39 as weak, 0.40–0.59 as moderate, 0.6–0.79 as strong and 0.8–1 as very strong correlation</figcaption></figure>

In this follow up study, we re-tested the TaqMan BLV qPCR developed and validated by NVRI (acting as reference WOAH laboratory) and the one adapted from this original protocol to be used with SYBR Green dye, allowing a significant reduction in costs [11]. Another 3 laboratories also performed NVRI's qPCR with slight modifications (i.e., Spain performed a multiplex assay for internal normalization). The remaining 6 labs introduced novel methodologies to the trial including one ddPCR (UY).

To compare different qPCR methods, a more comprehensive sample panel, accounting for a more geographical diversification was used in this trial. The amounts of BLV DNA in these samples were representative of the different BLV proviral loads found in field samples (from 1 to > 10,000 copies of BLV proviral DNA). Of note, 34% of reference samples had less than 100 copies of BLV DNA per 100 ng; samples were lyophilized to grant better preservation and reduced variability during distribution to participants around the globe.
<figure><img src="image_11.png" /><figcaption>**Fig. 6** Impact of BLV proviral copy numbers on the level of agreement. Violin plot for graphical presentation of Mann-Whitney U test. The test was performed to compare BLV provirus copy number in two groups of samples: 28 samples with fully concordant results from all eleven qPCR/ddPCR assays (left) and 15 samples with discordant results from different qPCR/ddPCR assays (right) (*P* < 0.001). Sample ID 42 was excluded from the statistical analysis</figcaption></figure>

The panel included a single negative control and two positive controls. Diagnostic sensitivity (DxSn) was estimated for each qPCR. Considering the 43 positive samples, the DxSn for the different qPCRs were: qPCR1=100%, qPCR2=82%, qPCR3=86%, qPCR4=84%, qPCR5=93%, ddPCR6=74%, qPCR7=93%, qPCR8=98%, qPCR9=98%, qPCR10=98% and qPCR11=95%. The most sensitive method was the qPCR1, and the method with the lowest sensitivity was the ddPCR6 method. Twenty-nine out of 44 samples were identified correctly by all qPCRs. The remaining 15 samples gave discordant results. The comparison of qualitative qPCR results among all raters revealed an overall observed agreement of 87%, indicating strong interrater reliability (Cohen's kappa = 0.396) [54, 55].

There are several factors that contribute to variability in qPCR results (i.e., number of copies of target input, sample acquisition, processing, storage and shipping, DNA purification, target selection, assay design, calibrator, data analysis, etc.). For that reason and as expected, the level of agreement among sister qPCRs (qPCR7, qPCR9-11) sharing similar protocols was higher compared to the rest of assays; this was also true for qPCR8 which targets the same region of BLV *pol* gene (shares same primers) but has a particular set-up to be used with SYBR Green chemistry. Oppositely, lower sensitivity and larger discrepancy against other tests was observed for the ddPCR6 and qPCR2-4.

Based on these observations we investigated which factors might have accounted for larger assessment variability amongst tests. In the first place, we observed that the use of different chemistries was not detrimental for the sensitivity and agreement among tests; similar DxSn and comparable level of agreement were obtained comparing TaqMan (qPCR7, 10, 11) vs SYBR Green (qPCR8) chemistries while targeting identical BLV sequence and using same standards. Also, when a multiplex qPCR (TaqMan) targeting the same BLV sequence and using the same standard was compared to previous ones, agreement was kept high, indicating that the lower sensitivity described for some multiplex qPCRs did not take place in this comparison. The use of an international calibrator and the efficiency estimation (standard curve) might inform variability associated with different chemistries. In contrast, another multiplex assay targeting another region of BLV *pol* (qPCR2) showed much lower sensitivity and agreement. As qPCR2 is performed as service by private company and oligonucleotide sequences were not available, we were not able to investigate in which proportion each of these two variables contributed to the lower performance of this assay, but we note the addition of 4 µl genomic DNA to this assay that would have an impact the DxSn. In this regard, there is substantial evidence showing that the variability of target sequence among strains from different geographical areas, might affect the sensitivity of BLV qPCRs. Previous studies comparing the *pol*, *gag*, *tax* and *env* genes reported that the *pol*
gene was the most suitable region to target for diagnostic purposes, since it provided the most-sensitive assays [11, 15, 56–59]. This might be due in part to higher sequence conservation of *pol* among strains from different geographical areas. Supporting this observation, it is noticeable how JPN qPCR improved their performance in the current trial, by targeting *pol* in place of *tax*, as it did in the previous interlaboratory trial. Since it is a commercial test, we cannot exclude other factors contributing for the performance upgrade observed for this qPCR. In the current study, qPCR3 and ddPCR6 targeting LTR and *env* sequences, showed lower performances than other assays. Standardization of DNA input into each qPCR would have likely resulted in higher concordance in results. For instance, qPCR1 added 10 µl of genomic DNA per reaction and ddPCR6 added 1 µl of genomic DNA, impacting the resulting sensitivity differences.

Since the sensitivity of each assay and, consequently, the level of agreement among assays might also be influenced by the number of BLV DNA copies present in each sample [48], we compared the average number of BLV DNA copies between a group of genomic DNA samples that gave concordant results and a group that gave discordant results, and observed that samples that gave discordant results had significantly lower numbers of BLV DNA copies than samples that gave concordant results. Related to this point, the degradation of target DNA during lyophilization, shipment and resuspension, could have been more significant in low-copy compared to high-copy samples. Consequently, the degradation of target DNA in samples with low copies of BLV DNA might have accounted for the greater level of discrepancy within this subset of samples. The rational of adding a large proportion of such samples (34% samples with less than 100 BLV copies per 100 ng of total DNA) was to mimic what is frequently observed in surveillance programs (i.e., hyperacute infection, chronic asymptomatic infection, etc.).

Quantitative methods for the detection of BLV DNA copies are important for segregation programs based on animal level of BLV PVL, as well as for scientific research and the study of BLV dynamics. When the numbers of copies of BLV DNA detected by different assays were compared, in the present study, we observed that although the ability to quantify BLV DNA differed among qPCRs/ddPCR and there were statistically significant differences in the distribution of copy numbers among assays, a strong average correlation was found for the eleven qPCRs/ddPCR. In this regard, the lack of an international calibrator (standard curve) could be a major contributor to the increment of quantitative variation amongst laboratories. For that reason, plasmid pBLV1 containing *pol* 120 bp sequence was originally

constructed for use as standard for quantification and shared with some collaborators (i.e., qPCR7, qPCR8, qPCR 9, qPCR10 and qPCR11). Remarkably, the laboratories used pBLV1 standard in the current trial obtained the most comparable results, indicating that the use of an international standard may have significant impact on the convergence of results; such standard reference material should be prepared under identical conditions. To avoid further variability a detailed protocol for lyophilized DNA sample resuspension, quantitation and template input into each qPCR should be shared with all participants.

## Conclusions

BLV DNA was detected with different level of sensitivity in serologically positive samples from different origin and classified into different BLV genotypes. Overall agreement was high; however, we found significant differences in results for the samples with low BLV DNA copy numbers. This second interlaboratory study demonstrated that differences in target sequence, DNA input and calibration curve standards can increase interlaboratory variability considerably. Next steps should focus on (i) standard unification (international gold standard) to estimate individual test efficiency and improve quantitative accuracy amongst tests; (ii) building a new panel of samples with low BLV DNA copy numbers to re-evaluate sensitivity and quantitation of molecular methods. Since no variation was observed in samples from different genotypes, all samples will be collected in Poland to standardize the collection, purification, lyophilization and shipping steps with precise instructions for suspension and constant input volume for the PCR reaction. Finally, we believe that following this standardization approach we will be able to improve overall agreement amongst tests, improving the diagnostic of BLV around the world.

## Abbreviations

<table><tr><td>ANOVA</td><td>One-way analysis of variance</td></tr><tr><td>BLV</td><td>Bovine leukemia virus</td></tr><tr><td>BLV-Abs</td><td>BLV-specific antibodies</td></tr><tr><td>ddPCR</td><td>Digital PCR</td></tr><tr><td>DxSn</td><td>Diagnostic sensitivity</td></tr><tr><td>EBL</td><td>Enzootic bovine leukemia</td></tr><tr><td>ELISA</td><td>Enzyme-linked immunosorbent assays</td></tr><tr><td>FRET</td><td>Real-time fluorescence resonance energy transfer PCR</td></tr><tr><td>GQN</td><td>Genomic quality number</td></tr><tr><td>H3F3A</td><td>Histone H3 family 3A housekeeping gene</td></tr><tr><td>ML</td><td>Maximum likelihood phylogenetic tree</td></tr><tr><td>MSA</td><td>Multiple-sequence alignment</td></tr><tr><td>PBLs</td><td>Peripheral blood leukocytes</td></tr><tr><td>PBS</td><td>Phosphate-buffered saline</td></tr><tr><td>PVL</td><td>Proviral load</td></tr><tr><td>qPCR</td><td>Quantitative real-time PCR</td></tr><tr><td>RT</td><td>Room temperature</td></tr><tr><td>WOAH</td><td>World Organisation for Animal Health</td></tr></table>
**Supplementary Information**

The online version contains supplementary material available at https://doi.org/10.1186/s12917-024-04228-z.

Additional file 1. Copy of the instruction included with the panel of 44 DNA samples sent to participating laboratories for dilution of the lyophilisates

Additional file 2. Detection of the H3F3A gene copy number in 43 DNA samples; no outlier was found for any samples (*P* <0.05) (two-sided).

Additional file 3. Concentration values of 44 DNA samples measured by the 11 participating laboratories (given in ng per µl)

Additional file 4. Post hoc - Dunn-Bonferroni-Tests. The Dunn-Bonferroni test revealed that the pairwise group comparisons of qPCR2 - qPCR4, qPCR3 - ddPCR6, qPCR4 - qPCR5, qPCR4 - ddPCR6, qPCR4 - qPCR9, qPCR4 - qPCR10, qPCR5 - qPCR11, ddPCR6 - qPCR11 and qPCR9 - qPCR11 have an adjusted *p*-value less than 0,05

Additional file 5. Kendall's Tau correlation coefficient values measured between each pair of assays. The numbers 1 to 11 in the first column and last row of the table indicate the names of the assays qPCR1-qPCR5, ddPCR6, qPCR7-qPCR11 respectively

Additional file 6. Maximum-likelihood phylogenetic analysis of full-length BLV-pol gene sequences representing 7 BLV genotypes (G1, G2, G3, G4, G6, G9, and G10) (A); (B) env-based sequences assigned to 10 BLV genotypes (G1, G2, G3, G4, G5, G6, G7, G8, G9, and G10); (C) LTR-based sequences representing 10 BLV genotypes (G1-G10). For all genes and LTR region the Tamura-Nei model and Bootstrap replications (1,000) were applied in MEGA X

Additional file 7. Multiple sequence alignment of reverse transcriptase, integrase, envelope and LTR sequences in the context of the specific primers used by different qPCR assays. (A) Multiple sequence alignment of reverse transcriptase (pol gene) sequences in the context of qPCR7, qPCR8, qPCR9, qPCR10 and qPCR11 assay primers. (B) Multiple sequence alignment of integrase (pol gene) sequences in the context of qPCR4 assay primers. (C) Multiple sequence alignment of env gene sequences in the context of ddPCR6. (D) Sequence alignment of LTR region sequences in the context of qPCR3 method primers

**Acknowledgements**

The authors thank Luc Willems (University of Liège, Belgium) for plasmid DNA sample pBLV344; Marlena Smagacz and Eliza Czarnecka (National Veterinary Research Institute, Poland) for lyophilizing DNA samples and DNA analysis, respectively; Ali Sakhawat (Animal Quarantine Department, Pakistan), Vitaliy Bolotin (National Scientific Center IECVM, Ukraine), Frank van der Meer and Sulav Shrestha (University of Calgary, Canada) for sharing material.

**Authors' contributions**

Proposed the conception and design of the study, A.P.; data curation, A.P., J.P.J., C.D., S.V., D.B., A.S., K.M., R.P., G.D., M.F.C. and CH.W.; investigation, A.P., VR., S.VW., S.V., A.J., M.J.R., K.N., M.L.B., M.L.G., PL., A.F., A.G. and S.B., formal analysis, A.P.; statistical analysis, A.P.; database analysis, A.P., visualization of the results, A.P.; resources, A.P., T.M.T. and J.K.; writing—original draft preparation, A.P., J.P.J.; writing—review and editing, A.P., J.P.J., C.D., S.VW., T.M.T. and J.K.; project administration, A.P. All authors read and approved the submitted version.

**Funding**

The APC was funded by the National Veterinary Research Institute, Pulawy, Poland.

**Availability of data and materials**

Not applicable.

**Data availability**

No datasets were generated or analysed during the current study.

**Declarations**

**Ethics approval and consent to participate**

The study was approved by the Veterinary Sciences Animal Care Committee No. AC21-0210, Canada; the Institutional Animal Care and Use Committee No. PROTO202000096 from 4/13/2020 to 4/14/2023, Michigan State University, United States and the Ethics Review Board, COMSATS Institute of Information Technology, Islamabad, Pakistan, no. CIIT/Bio/ERB/17/26. Blood samples from Polish, Moldovan and Ukrainian cattle, naturally infected with BLV, were selected from collections at local diagnostic laboratories as part of the Enzootic bovine leukosis (EBL) monitoring program between 2012 and 2018 and sent to the National Veterinary Research Institute (NVRI) in Pulawy for confirmation study. The approval for collection of these samples from ethics committee was not required according to Polish regulation ("Act on the Protection of Animals Used for Scientific or Educational Purposes", Journal of Laws of 2015). All methods were carried out in accordance with relevant guidelines and regulations. The owners of the cattle herds from which the DNA samples originated, the district veterinarians caring for these farms and the ministries of agriculture were informed and consented to the collection of blood from the animals for scientific purposes and the sending of samples to NVRI.

**Consent for publication**

Not applicable.

**Competing interests**

The authors declare no competing interests.

**Author details**

¹Department of Biochemistry, National Veterinary Research Institute, Pulawy 24-100, Poland. ²Instituto de Virologia E Innovaciones Tecnológicas (IVIT), Centro de Investigaciones en Ciencias Veterinarias y Agronómicas (CICVya), Instituto Nacional de Tecnología Agropecuaria (INTA) - CONICET, Buenos Aires, Argentina. ³CentralStar Cooperative, 4200 Forest Rd, Lansing, MI 48910, USA. ⁴Department of Animal Science, College of Agriculture and Natural Resources, Michigan State University, East Lansing, Michigan 48824, USA. ⁵Niotr Laboratory, Unit Pathology and Welfare of Ruminants, French Agency for Food, Environmental and Occupational Health and Safety (Anses), Ploufragan-Plouzané, Niot, France. ⁶Croatian Veterinary Institute, Savska Cesta 143, Zagreb 10000, Croatia. ⁷Laboratorio Central de Veterinaria (LCV), Ministry of Agriculture, Fisheries and Food, Carretera M-106 (Km 1,4), Madrid, Algete 28110, Spain. ⁸Department of Veterinary Sciences, Faculty of Agriculture, Iwate University, 3-18-8 Ueda, Morioka 020-8550, Japan. ⁹Departamento de Patobiología, Facultad de Veterinaria, Unidad de Microbiología, Universidad de La República, Ruta 8, Km 18, Montevideo 13000, Uruguay. ¹⁰Laboratorio de Virología, Departamento SAMP, Centro de Investigación Veterinaria de Tandil-CIVETAN (CONICET/UNCPBA/CICPBA), Buenos Aires, Argentina. ¹¹Laboratório Federal de Defesa Agropecuária de Minas Gerais, Pedro Leopoldo, Brazil. ¹²Department of Pathobiology, College of Veterinary Medicine, Auburn University, Auburn, AL 36849-5519, USA. ¹³Department of Omics Analyses, National Veterinary Research Institute, 24-100, Pulawy, Poland.

Received: 24 November 2023 Accepted: 9 August 2024

Published online: 26 August 2024

**References**

1. Coffin JM, Hughes SH, Varmus HE. (Eds.). 1650–1655 (1997). Retroviruses. Cold Spring Harbor Laboratory Press.

2. Ghysdael J, Bruck C, Kettmann R, Burny A. Bovine leukemia virus. Curr Top Microbiol Immunol. 1984;112:1–19.

3. Ott SL, Johnson R, Wells SJ. Association between bovine-leukosis virus seroprevalence and herd-level productivity on US dairy farms. Prev Vet Med. 2003;61:249–62.

4. Bartlett PC, et al. Options for the control of bovine leukemia virus in dairy cattle. J Am Vet Med Assoc. 2014;244:914–22.

5. Kuczewski A, et al. Economic evaluation of 4 bovine leukemia virus control strategies for Alberta dairy farms. J Dairy Sci. 2019;102:2578–92.
6. Frie MC, Coussens PM. Bovine leukemia virus: a major silent threat to proper immune responses in cattle. Vet Immunol Immunopathol. 2015;163:103–14.

7. Panel, E.A. Scientific opinion on enzootic bovine leukaosis. EFSA J. 2015;13:4188.

8. OIE. World Animal Health Information Database - Version: 1.4. World Animal Health Information Database. Paris, France: World Organisation for Animal Health; 2009. Available from: http://www.oie.int. Accessed 16 Aug 2024.

9. Health, W.O.f.A. Manual of diagnostic tests and vaccines for terrestrial animals. Infect Bursal Dis. 2012;12:549–65.

10. Hutchinson HC, et al. Bovine leukemia virus detection and dynamics following experimental inoculation. Res Vet Sci. 2020;133:269–75.

11. Rola-Luszczak M, Finnegan C, Olech M, Choudhury B, Kuzmak J. Development of an improved real time PCR for the detection of bovine leukaemia provirus nucleic acid and its use in the clarification of inconclusive serological test results. J Virol Methods. 2013;189:258–64.

12. Nakada S, Kohara J, Makita K. Estimation of circulating bovine leukemia virus levels using conventional blood cell counts. J Dairy Sci. 2018;101:11229–36.

13. Ruggiero VJ, Bartlett PC. Control of Bovine Leukemia Virus in Three US Dairy Herds by Culling ELISA-Positive Cows. Vet Med Int. 2019;19:3202184.

14. Kobayashi T, et al. Increasing Bovine leukemia virus (BLV) proviral load is a risk factor for progression of Enzootic bovine leucosis: A prospective study in Japan. Prev Vet Med. 2020;178: 104680.

15. Jaworski JP, Pluta A, Rola-Luszczak M, McGowan SL, Finnegan C, Heenemann K, Carignano HA, Alvarez I, Murakami K, Willems L, Vahlenkamp TW, Trono KG, Choudhury, B, Kužmak J. Interlaboratory Comparison of Six Real-Time PCR Assays for Detection of Bovine Leukemia Virus Proviral DNA. J Clin Microbiol. 2018;56(7):e00304-18. https://doi.org/10.1128/JCM.00304-18.

16. Pluta A, Rola-Luszczak M, Douville RN, Kužmak J. Bovine leukemia virus long terminal repeat variability: identification of single nucleotide polymorphisms in regulatory sequences. Virol J. 2018;15:165.

17. Pluta A, Willems L, Douville RN, Kužmak J. Effects of Naturally Occurring Mutations in Bovine Leukemia Virus 5'-LTR and Tax Gene on Viral Transcriptional Activity. Pathog. 2020;9(10):836. https://doi.org/10.3390/pathogens9100836.

18. Pluta A, et al. Molecular characterization of bovine leukemia virus from Moldovan dairy cattle. Arch Virol. 2017;162:1563–76.

19. Rola-Luszczak M, Sakhawat A, Pluta A, Rylo A, Bomba A, Bibi N, Kužmak J. Molecular Characterization of the env Gene of Bovine Leukemia Virus in Cattle from Pakistan with NGS-Based Evidence of Virus Heterogeneity. Pathogens (Basel, Switzerland). 2021;10(7):910. https://doi.org/10.3390/pathogens10070910.

20. Rola-Luszczak M, 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.

21. Pinheiro de Oliveira TF, et al. Detection of contaminants in cell cultures, sera and trypsin. Biologicals. 2013;41:407–14.

22. Pluta A, Blazhko NV, Ngirande C, Joris T, Willems L, Kužmak J. Analysis of Nucleotide Sequence of Tax, miRNA and LTR of Bovine Leukemia Virus in Cattle with Different Levels of Persistent Lymphocytosis in Russia. Pathogens. 2021;10(2):246. https://doi.org/10.3390/pathogens10020246.

23. Yang Y, et al. Bovine leukemia virus infection in cattle of China: Association with reduced milk production and increased somatic cell score. J Dairy Sci. 2016;99:3688–97.

24. DeGraves FJ, Gao D, Kaltenboeck B. High-sensitivity quantitative PCR platform. Biotechniques. 2003;34(106–110):112–105.

25. Fonseca Junior AA, et al. Evaluation of three different genomic regions for detection of bovine leukemia virus by real-time PCR. Braz J Microbiol. 2021;52:2483–8.

26. Farias MVN, et al. Toll-like receptors, IFN-gamma and IL-12 expression in bovine leukemia virus-infected animals with low or high proviral load. Res Vet Sci. 2016;107:190–5.

27. Holland PM, Abramson RD, Watson R, Gelfand DH. Detection of specific polymerase chain reaction product by utilizing the 5'---3' exonuclease activity of Thermus aquaticus DNA polymerase. Proc Natl Acad Sci U S A. 1991;88:7276–80.

28. De Brun ML, et al. Development of a droplet digital PCR assay for quantification of the proviral load of bovine leukemia virus. J Vet Diagn Invest. 2022;34:439–47.

29. Rola-Luszczak M, Finnegan C, Olech M, Choudhury B, Kužmak J. Development of an improved real time PCR for the detection of bovine leukaemia provirus nucleic acid and its use in the clarification of inconclusive serological test results. J Virol Methods. 2013;189:258–64.

30. Petersen MI, Alvarez I, Trono KG, Jaworski JP. Quantification of bovine leukemia virus proviral DNA using a low-cost real-time polymerase chain reaction. J Dairy Sci. 2018;101:6366–74.

31. Toussaint JF, Sailleau C, Breard E, Zientara S, De Clercq K. Bluetongue virus detection by two real-time RT-qPCRs targeting two different genomic segments. J Virol Methods. 2007;140:115–23.

32. John EE, et al. Development of a predictive model for bovine leukemia virus proviral load. J Vet Intern Med. 2022;36:1827–36.

33. Farias MVN, et al. Toll-like receptors, IFN-γ and IL-12 expression in bovine leukemia virus-infected animals with low or high proviral load. Res Vet Sci. 2016;107:190–5.

34. Yoneyama S, et al. Comparative Evaluation of Three Commercial Quantitative Real-Time PCRs Used in Japan for Bovine Leukemia Virus. Viruses. 2022;14:1182.

35. Polat M, Takeshima SN, Aida Y. Epidemiology and genetic diversity of bovine leukemia virus. Virol J. 2017;14:209.

36. Lee E, et al. Molecular epidemiological and serological studies of bovine leukemia virus (BLV) infection in Thailand cattle. Infect Genet Evol. 2016;41:245–54.

37. Duran-Yelken S, Alkan F. Molecular analysis of the env, LTR, and pX regions of bovine leukemia virus in dairy cattle of Türkiye. Virus Genes. 2024;60:173–85.

38. Lv G, Wang J, Lian S, Wang H, Wu R. The Global Epidemiology of Bovine Leukemia Virus: Current Trends and Future Implications. Animals. 2024;14(2):297. https://doi.org/10.3390/ani14020297.

39. Úsuga-Monroy C, Díaz FJ, Echeverri-Zuluaga JJ, González-Herrera LG, López-Herrera A. Presence of bovine leukemia virus genotypes 1 and 3 in Antioquia, Colombia. Revista UDCA Actualidad & Divulgación Científica. 2018;21:119–26.

40. Úsuga-Monroy C, Díaz FJ, González-Herrera LG, Echeverry-Zuluaga JJ, López-Herrera A. Phylogenetic analysis of the partial sequences of the env and tax BLV genes reveals the presence of genotypes 1 and 3 in dairy herds of Antioquia. Colombia VirusDisease. 2023;34:483–97.

41. Martin D, et al. Comparative study of PCR as a direct assay and ELISA and AGID as indirect assays for the detection of bovine leukaemia virus. J Vet Med B Infect Dis Vet Public Health. 2001;48:97–106.

42. Cohen J. A Coefficient of Agreement for Nominal Scales. Educ Psychol Measur. 1960;20:37–46.

43. Feinstein AR, Cicchetti DV. High agreement but low kappa: I. The problems of two paradoxes. J Clin Epidemiol. 1990;43:543–9.

44. Woolhouse MEJ, et al. Heterogeneities in the transmission of infectious agents: Implications for the design of control programs. Proc Natl Acad Sci. 1997;94:338–42.

45. Ohshima K, Okada K, Numakunai S, Kayano H, Goto T. An eradication program without economic loss in a herd infected with bovine leukemia virus (BLV). Nihon Juigaku Zasshi. 1988;50:1074–8.

46. Juliarena MA, Gutierrez SE, Ceriani C. Determination of proviral load in bovine leukemia virus-infected cattle with and without lymphocytosis. Am J Vet Res. 2007;68:1220–5.

47. Mirsky ML, Olmstead CA, Da Y, Lewin HA. The prevalence of proviral bovine leukemia virus in peripheral blood mononuclear cells at two subclinical stages of infection. J Virol. 1996;70:2178–83.

48. Eaves FW, Molloy JB, Dimmock CK, Eaves LE. A field evaluation of the polymerase chain reaction procedure for the detection of bovine leukaemia virus proviral DNA in cattle. Vet Microbiol. 1994;39:313–21.

49. Juliarena MA, Barrios CN, Ceriani MC, Esteban EN. Hot topic: Bovine leukemia virus (BLV)-infected cows with low proviral load are not a source of infection for BLV-free cattle. J Dairy Sci. 2016;99:4586–9.

50. Yuan Y, 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–54.

51. Taxis TM, et al. Reducing bovine leukemia virus prevalence on a large midwestern dairy farm by using lymphocyte counts, ELISA antibody testing, and proviral load. The Bovine Practitioner. 2020;54:136–44.
52. McHugh ML. Interrater reliability: the kappa statistic. Biochem Med (Zagreb). 2012;22:276-82.

53. Cicchetti DV, Feinstein AR. High agreement but low kappa: II. Resolving the paradoxes J Clin Epidemiol. 1990;43:551-8.

54. Heenemann K, et al. Development of a Bovine leukemia virus polymerase gene-based real-time polymerase chain reaction and comparison with an envelope gene-based assay. J Vet Diagn Invest. 2012;24:649-55.

55. Lew AE, et al. Sensitive and specific detection of proviral bovine leukemia virus by 5' Taq nuclease PCR using a 3' minor groove binder fluorogenic probe. J Virol Methods. 2004;115:167-75.

56. Takeshima SN, Sasaki S, Meripet P, Sugimoto Y, Aida Y. Single nucleotide polymorphisms in the bovine MHC region of Japanese Black cattle are associated with bovine leukemia virus proviral load. Retrovirology. 2017;14:24.

57. Debacq C, et al. Reduced proviral loads during primo-infection of sheep by Bovine Leukemia virus attenuated mutants. Retrovirology. 2004;1:31.

58. Kuckleburg CJ, et al. Detection of bovine leukemia virus in blood and milk by nested and real-time polymerase chain reactions. J Vet Diagn Invest. 2003;15:72-6.

59. Dube S, et al. Degenerate and specific PCR assays for the detection of bovine leukaemia virus and primate T cell leukaemia/lymphoma virus pol DNA and RNA: phylogenetic comparisons of amplified sequences from cattle and primates from around the world. J Gen Virol. 1997;78(Pt 6):1389-98.

**Publisher's Note**

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