Skip to main content

Development of a duplex real-time recombinase aided amplification assay for the simultaneous and rapid detection of PCV3 and PCV4

Abstract

Background

Porcine circoviruses 3 (PCV3) and 4 (PCV4) are emerging pathogens with global implications for swine industry, disturbing the diagnosis of PCVs associated diseases due to a range of similar clinical symptoms and increasingly coinfections. A rapid and accurate method for detection of PCV3 and PCV4 is critical for controlling the transmission of associated disease.

Methods

We developed a duplex real-time recombinase aided amplification (RAA) assay for detection of both PCV3 and PCV4 simultaneously. The assay was completed within 20 min at 39℃ with the designed optimal primers and probes.

Results

The established assay was more convenient and simpler operation compared with conventional molecular biological assays. The assay achieved a detection limit of 73.67 copies/reaction for each circovirus (at 95% probability by probit regression analysis) and showed high specificity and no cross-reactivity with other important porcine viruses (including PCV2). The intra- and inter-group coefficients of variation (CV) were ranged from 2.08 to 4.97%, indicating high stability and reliability. Comparative analysis with PCV3 and PCV4 qPCR on 60 clinical samples and artificially spiked samples indicated high congruence (the kappa value was 0.966 and 1, respectively, with p < 0.001), with only minor discrepancies, validating effectiveness of the duplex RAA assay in detecting co-infections and its suitability for preliminary clinical diagnosis of PCV3 and PCV4.

Conclusions

This study provides a robust basis for multiplex detection of veterinary pathogens using RAA technique, enhancing the field’s capacity to control PCV3 and PCV4, and supporting reliable aid for epidemiological understanding of emerging circoviruses.

Introduction

Circoviruses are notable members among eukaryotic circular Rep-encoding ssDNA (CRESS DNA) viruses, belonging to the genus Circovirus of the Circoviridae family [1]. As ubiquitous pathogens in swine industry, porcine circoviruses (PCVs) are small, icosahedral, and nonenveloped viruses containing a covalently closed, circular single-stranded DNA. Currently, four PCV species have been characterized, including porcine circovirus type 1 (PCV1), porcine circovirus type 2 (PCV2), porcine circovirus type 3 (PCV3), and porcine circovirus type 4 (PCV4) [2]. PCV1 was initially isolated and identified from pig kidney cell lines as a contaminant in cell culture, and its non-pathogenicity was subsequently indicated by in vivo experiments [2, 3]. PCV2 is an important pathogen causing porcine circovirus associated disease (PCVAD), which is associated with multiple clinical disease syndrome in pigs, and causes great economic losses to the global animal husbandry industry [4]. With the development of metagenomic next-generation sequencing in the field of pathogenic microorganism detection and diagnosis, PCV3 and PCV4 were discovered successively in the United States and China [5, 6]. As newly emerging pathogens, PCV3 and PCV4 had been detected in domestic pig and wild boar herds in many countries and regions worldwide, and occasionally found in other healthy animals [7,8,9,10,11]. Current clinical symptoms associated with PCV3 and PCV4 infection mainly included porcine dermatitis and nephropathy syndrome (PDNS), reproductive failure, respiratory disorders and diarrhea [2, 7, 8, 12].

Under field or experimental conditions, PCV3 or PCV4 infection alone seemed like mild and may even not present subclinical signs, while co-infection or secondary infection with other PCVs and swine pathogens can exacerbate clinical manifestations [1, 7, 8]. Due to above factors and similar clinical symptoms between them, differential diagnostic assay for rapid and sensitive detection of PCV3 and PCV4 is critical for controlling the transmission of associated disease. Currently, the most commonly used diagnostic methods for PCV3 and PCV4 infections are singleplex or multiplex and SYBR green-based or TaqMan probe-based real-time qPCR [7, 8]. Other novel molecular biological techniques, including droplet digital PCR (ddPCR), loop-mediated isothermal amplification (LAMP), lateral flow strip (LFD), CRISPR-Cas12a/Cas13a technology, recombinase aided amplification (RAA), also had been subsequently applied for rapidly detecting single infection of PCV3 or PCV4 [7, 13,14,15,16]. qPCR offers comparable sensitivity for the detection of PCV3 and PCV4 but is constrained by its duration, procedural complexity, and the need for specialized equipment and expertise, rendering it impractical for rapid diagnostic purposes. Furthermore, existing rapid detection approaches for PCV3 and PCV4 primarily focus on individual pathogens, failing to provide simultaneous identification of both. This limitation necessitates multiple screenings and results in a lack of diagnostic comprehensiveness. Therefore, it is necessary to develop a suitable duplex rapid diagnostic method for these two emerging circoviruses.

Recombinase aided amplification (RAA) is a novel isothermal nucleic acid amplification technique that stands out for its simple procedure and speed. It operates at a constant temperature of 37–42 °C and can complete the amplification process within 15–30 min. The RAA system hinges on three key enzymes: a recombinase that pairs specific primers to the DNA template, a single-stranded DNA-binding protein (SSB) that facilitates the opening of the DNA structure, and a DNA polymerase responsible for amplification and extension [17, 18]. Additionally, its amplification result can visually output via multiple ways, containing gel electrophoresis, direct blue illumination, specific probes-based real-time fluorescence monitoring, lateral flow dipsticks and CRISPR-Cas12a/Cas13a system [19,20,21,22,23]. Those advantages make singleplex RAA widely applicable for rapid detection of various pathogenic microorganisms in public and veterinary health.

The duplex or multiplex real-time RAA, which employs specific primers and probes to amplify and detect multiple target DNA sequences concurrently in a single reaction, offers a rapid, sensitive, specific and convenient method for simultaneous pathogen detection through real-time fluorescence monitoring. However, only few studies have reported the application of multiplex RAA technique in the detection of pathogens at present [24,25,26].

In this study, we developed a duplex real-time RAA assay for the rapid, sensitive, and simultaneous detection of PCV3 and PCV4, two emerging porcine circoviruses. This duplex RAA-based system is not only sensitive but also more rapid, convenient and simpler to operate than conventional microbiological, PCR, or qPCR-based techniques, offering a robust research foundation for the multiplex detection of veterinary pathogens using the RAA technique, providing reliable technical assist to PCV3 and PCV4 clinical diagnosis, limiting geographic distribution and timely preventing and controlling transmission of associated disease.

Materials and methods

Standard plasmid for duplex real-time RAA

The complete genome of PCV3 (2, 000 bp) and PCV4 (1, 770 bp) were synthesized according to the PCV3 and PCV4 genomic sequence (GenBank accession no. KT869077 and MT311852) and cloned into the pUC57 vector by Sangon Biotech Co., Ltd. (Sangon Biotech Co., Ltd., Shanghai, China) to construct standard plasmids for duplex real-time RAA assay, designated as pUC57-PCV3 and pUC57-PCV4. The concentration of the standard plasmids was measured using NanoDrop and Qubit 2.0 fluorometer (Thermo Fisher Scientific, MA, USA). The DNA copy number was calculated using the following formula: DNA copy number/µL = [plasmid concentration (ng/µL)×10− 9 × 6.02 × 1023]/[DNA length (nt)×660]. The gene sequencing was processed by Sangon Biotech Co., Ltd. (Sangon Biotech Co., Ltd., Shanghai, China). The standard plasmids pUC57-PCV3 and pUC57-PCV4 were stored at -20℃ until use.

Design and synthesis of primers and probes

The whole genomic sequences of 30 strains (including genotype 3a, 3b and 3c) of PCV3 and 29 strains (including genotype 4a, 4b and 4c) of PCV4 from GenBank database were aligned respectively using the Lasergene version 7.1 software (DNASTAR Inc., WI, USA) to identify the highly conserved regions. Primer Premier version 5.0 software (PREMIER Biosoft, CA, USA) was used to design RAA primers and probes according to the criteria suggested in the TwistAmp™ amplification guidelines (TwistDx Ltd., Cambridge, UK). An in silico check for the design of the primer pair and probe was conducted utilizing the NCBI Primer-BLAST web tool (https://www.ncbi.nlm.nih.gov/tools/primer-blast). For the development of a duplex real-time RAA assay detecting PCV3 and PCV4, two primer pairs and two RAA probes were designed to target the ORF2 gene of PCV3 (Fig. 1). Concurrently, three primer pairs and two RAA probes were designed to target the ORF1 or ORF2 gene of PCV4 (Figs. 2 and 3). Details of the primers and probes for duplex real-time RAA in this study are shown in Table 1. All the primers and probes were synthesized by Sangon Biotech Co., Ltd. (Sangon Biotech Co., Ltd., Shanghai, China).

Fig. 1
figure 1

Positions of the duplex real-time RAA primers and probes in the ORF2 (Cap gene) sequences of different PCV3 strains in the GenBank database. Dots represent nucleotide residues that match the majority. The forward primers are marked with red bars, the reverse primers are marked with blue bars, and probes are marked with green bars

Fig. 2
figure 2

Positions of the duplex real-time RAA primers and probes in the ORF1 (Rep gene) sequences of different PCV4 strains in the GenBank database. Dots represent nucleotide residues that match the majority. The forward primers are marked with red bars, the reverse primers are marked with blue bars, and probes are marked with green bars

Fig. 3
figure 3

Positions of the duplex real-time RAA primers and probes in the ORF2 (Cap gene) sequences of different PCV4 strains in the GenBank database. Dots represent nucleotide residues that match the majority. The forward primers are marked with red bars, the reverse primers are marked with blue bars, and probes are marked with green bars

Table 1 Primers and probes used for the duplex RAA

Duplex real-time RAA assay

The reaction mixture of duplex real-time RAA assay was prepared using the commercial probe-based RAA Kit (Hangzhou Zhongce Biotechnology Co., Ltd., Hangzhou, China). The reaction mixture (50 µL reaction volume) contained Buffer A (25 µL), 2 µL each of PCV3 and PCV4 forward and reverse primers (10 µM), 0.6 µL each of PCV3 and PCV4 probes (10 µM), template DNA (5 µL), nuclease-free water (8.3 µL) and Buffer B (2.5 µL). Briefly, components except template DNA and Buffer B were premixed and transferred into the reaction tube containing RAA enzyme dry powder and overturned to mix thoroughly. Add 5 µL of template DNA to the reaction tube and pipetting 2.5 µL of Buffer B into the tube lid. After a brief vortex mixing and centrifugation, the reaction tube containing the reaction mixture was promptly placed in a ViiA™ 7 real-time PCR instrument (Thermo Fisher Scientific, MA, USA) and incubated at 39 °C for approximately 20 min (1 cycle at 39 °C for 40 s and 40 cycles at 39 °C for 30 s). The fluorescence signal was monitored in real-time, which collected every 30 s. Samples that gave rise to an exponential amplification curve above the negative control threshold within 20 min were judged as positive.

Optimization of reaction system

Optimization experiments were conducted to determine the optimal duplex real-time RAA system. In the optimization experiments, 2.5 µL standard plasmid pUC57-PCV3 (20.0 ng/µL) and 2.5 µL standard plasmid pUC57-PCV4 (20.0 ng/µL) were mixed to serve as templates in the reaction system. 6 sets of primer-probe combinations (PCV3 F1R1 P1/PCV4 F1R1 P1, PCV3 F1R1 P1/PCV4 F2R2 P1, PCV3 F1R1 P1/ PCV4 F3R3 P2, PCV3 F2R2 P2/PCV4 F1R1 P1, PCV3 F2R2 P2/PCV4 F2R2 P1, PCV3 F2R2 P2/PCV4 F3R3 P2) were screened using commercial probe-based RAA Kit under reaction conditions recommended by the kit manual to determine the optimal probe and primer combination.

The concentrations of the best performing primer-probe combination were furtherly optimized by setting the concentration of a single variable component while all other conditions were fixed as recommended by the kit manual. Primer pairs concentration was diluted to 200, 400, and 600 nM, and combinations of concentration were set in Table 2. Probe concentration gradient was set to 60, 90, 120,150 and 180 nM.

Table 2 PCV3 and PCV4 primer pairs concentration ratio

Nucleic acid extraction

Viral nucleic acids of clinical samples and artificially spiked samples were extracted from 200 µL of sample lysate with the MagNA Pure LC Total Nucleic Acid Isolation Kit (Roche Diagnostics, Branchburg, NJ, USA), following the manufacturer’s instructions. The extracted nucleic acids were eluted with 50 µL of nuclease-free water and stored at -80℃ until used for duplex real-time RAA assay and qPCR assay. Standard plasmids of PCV3 and PCV4 for duplex real-time RAA were extracted using a Plasmid Mini Kit (Sangon Biotech Co., Ltd., Shanghai, China), according to the manufacturer’s instructions.

Analytical specificity

The specificity of the duplex real-time RAA assay for PCV3 and PCV4 detection was evaluated by using viral DNA or cDNA from other important swine pathogens and standard plasmids (pUC57-PCV3 and pUC57-PCV4) under the optimal reaction system and conditions. The swine pathogens included FMDV (Foot and mouth disease virus), CSFV (Classical swine fever virus), PRV (Pseudorabies virus), PPV (Porcine parvovirus), PCV2 (Porcine circovirus 2), PEDV (Porcine epidemic diarrhea virus), TGEV (Transmissible gastroenteritis virus), PRRSV (Porcine reproductive and respiratory syndrome virus), and PDCoV (Porcine deltacoronavirus). Viral nucleic acids of the above 9 viral pathogens were extracted and maintained in our laboratory. Viral cDNA was obtained using the PrimeScript™ cDNA Synthesis Kit (Takara Biomedical Technology Co., Ltd., Dalian, China), according to the manufacturer’s instructions.

Analytical sensitivity

To evaluate the sensitivity of the duplex real-time RAA assay for PCV3 and PCV4, standard plasmids (pUC57-PCV3 and pUC57-PCV4) was subjected to 10-fold serial dilution ranging from 107 to 100 copies per 2.5 µL. Each dilution was used as a template to carry out the duplex real-time RAA assay according to the optimized system. For a more accurate analysis of the detection limit, 8 independent runs were performed using the above dilution series, and the data were used for the calculation of the 95% limit of detection (LOD) by probit regression analysis using SPSS 24.0 software (SPSS Inc., Chicago, IL, USA).

Analysis of repeatability and reproducibility

Three different dilutions (104, 103, and 102 copies per 2.5 µL) of standard plasmids (pUC57-PCV3 and pUC57-PCV4) were used as templates for the standard plasmids (pUC57-PCV3 and pUC57-PCV4) to assess the intra-assay repeatability and inter-assay reproducibility. In each reaction system, the template consisted of 2.5 µL of diluted pUC57-PCV3 and 2.5 µL of diluted pUC57-PCV4, and both at the same concentration. Each dilution reaction was tested in triplicate in one run or in three independent runs on separate days. The coefficients of variation (CVs) were obtained by calculating the cycle threshold (CT).

Detection of clinical samples and artificially spiked samples

To evaluate the clinical performance of the duplex real-time RAA assay for PCV3 and PCV4, 60 samples (55 clinical samples and 5 artificially spiked samples) were prepared for detection by the developed assay. 55 clinical samples were screened using the PCV3 qPCR assay (according to local standards in China) and PCV4 qPCR assay reported before respectively [6], and the results showed no qPCR positive samples of PCV4. Therefore, standard plasmid of PCV4 (pUC57-PCV4) was mixed with 3 qPCR positive samples of PCV3 and 5 artificially spiked samples respectively. Specifically, 50 µL of pUC57-PCV4 (205.8 ng/µL) was individually mixed with 1 g, 2 g, and 2.6 g of lymph node samples that tested qPCR positive for PCV3. For the artificially spiked samples, 50 µL of pUC57-PCV4 (205.8 ng/µL) was also individually mixed with 2 g, 3.5 g, and 5 g of pork, as well as 1.5 g and 2 g of lymph node samples. Viral nucleic acid was extracted from each of these mixtures and stored at -20 °C before use. 60 samples were also detected in parallel using the PCV3 and PCV4 qPCR assays respectively for comparison.

Statistical analysis

Kappa statistics were applied to compare the coincidence rates between the duplex real-time RAA assay and qPCR assays using SPSS 24.0 software (SPSS Inc., Chicago, IL, USA). Statistical analyses and data plotting were performed using the GraphPad Prism software (La Jolla, CA, USA).

Results

Screening of the optimal primers and probe for the duplex real-time RAA assay

To determine the best diagnostic target for duplex real-time RAA of PCV3 and PCV4, the whole genomic sequences of 30 different strains (including genotype 3a, 3b and 3c) of PCV3 and 29 PCV4 strains (including genotype 4a, 4b and 4c) were aligned (dots indicate nucleotide residues that matched the majority). According to the alignment results, a region of PCV3 ORF2 gene exhibited high conservation was selected as the amplification candidate for duplex real-time RAA assay, and two pairs of primers and two probes of PCV3 were designed based on it. Three pairs of primers and two probes of PCV4 were designed based on conserved region of ORF1 and ORF2 gene (Table 1). As shown in the alignment results of PCV3 or PCV4 genomic sequences, forward primers are marked with red bars, reverse primers are marked with blue bars, and probes marked with green bars (Figs. 1, 2 and 3). The two thymine (T) residues in designed probes of PCV3 and PCV4 labeled with a fluorophore (HEX or FAM) and quencher (BHQ1) are fully conserved among the representative strains (Table 1). Based on above designed primers and probes, 6 set of primer-probe combinations were generated and selected as candidates to determine the optimal probe and primer combination. The screening experiments showed that PCV3 F2R2 P2 and PCV4 F2R2 P1 respectively produced better amplification performance in targeting their own conserved region with more fast speed and higher fluorescence value, compared with other preliminarily designed candidates (Fig. 4). Therefore, we confirmed the PCV3 F2R2 P2/PCV4 F2R2 P1 set to be the optimal primer-probe combination and used in the subsequent duplex real-time RAA assay for PCV3 and PCV4.

Fig. 4
figure 4

Screening the optimal primer-probe combination for the duplex real-time RAA for PCV3 and PCV4. Representative amplification results of the duplex real-time RAA assay using 6 primer-probe combinations. Amplification curves in different colors represent different primer-probe combinations. (A) The amplification result of the FAM channel, corresponding to the amplification curves of PCV4; (B) The amplification result for the HEX channel, corresponding to the amplification curves of PCV3; (C) The amplification results of both channels

Reaction optimization for the duplex real-time RAA assay

To achieve an optimized duplex RAA reaction system, an evaluation of the primer pairs and probes concentrations were undertaken through a series of dilution gradient assays, respectively. PCV3 and PCV4 primer pairs concentration was diluted to 200, 400, and 600 nM, and volume ratio of 6 groups were set in Table 2. The primer concentration grading experiment showed that group 5 (400 nM of PCV3 primer pair and 400 nM of PCV4 primer pair) had better performance in both PCV3 and PCV4 amplification curve among other experimental groups (Fig. 5A and C). Given that the optimal concentration of the PCV3 and PCV4 primer pairs was consistent, probes of PCV3 and PCV4 were simultaneously diluted to a series of identical concentrations (60, 90, 120,150 and 180 nM) to select the optimal working concentration. As shown in Figs. 5D and F and 150 nM of PCV3 and PCV4 probes generated the best amplification results with stronger fluorescence signal and the shorter reaction time. Consequently, for subsequent duplex real-time RAA assay, the optimal primer pairs and probe concentration were respectively determined at 400 nM and 150 nM.

Fig. 5
figure 5

Optimization of primer pairs and probes concentrations for the duplex real-time RAA assay of PCV3 and PCV4. Amplification curves in different colors represent different concentration combinations of primer pairs or different concentrations of probes. The amplification result of the FAM channel, corresponding to the amplification curves of PCV4. The amplification result for the HEX channel, corresponding to the amplification curves of PCV3. (A-C) Representative amplification results of the duplex real-time RAA assay using 6 concentration combinations of primer pairs; (D-F) Representative amplification results of the duplex real-time RAA assay using 5 concentrations of probes

Analytical specificity of the duplex real-time RAA assay

The specificity of PCV3 and PCV4 duplex real-time RAA assay was determined using the nucleic acid from 9 clinically important swine pathogens and standard plasmids (pUC57-PCV3 and pUC57-PCV4) under the optimal reaction system. As shown in Fig. 6A, pUC57-PCV3, pUC57-PCV4 and mixture of both were positively amplified, and the remaining other swine pathogens (including FMDV、CSFV、PRV、PPV、PCV2、PEDV、TGEV、PRRSV、PDCoV) were negative. These results demonstrated that duplex real-time RAA assay had no cross-reaction with other common swine viruses, exhibiting a good specificity for PCV3 and PCV4 detection.

Fig. 6
figure 6

Specificity and sensitivity analysis of the duplex real-time RAA assay for PCV3 or PCV4. (A) Analytical specificity test of the duplex real-time RAA assay. Amplification curves in blue corresponding to the FAM channel, and those in green corresponding to the HEX channel. Curves 1: pUC57-PCV3 and pUC57-PCV4, 2: pUC57-PCV3, 3: pUC57-PCV4, 4–13: nucleic acids of FMDV, CSFV, PRV, PPV, PCV2, PEDV, TGEV, PRRSV, PDCoV, and the negative control, respectively; (B) Analytical sensitivity test of the duplex real-time RAA assay. Amplification curves in blue corresponding to the FAM channel, and those in green corresponding to the HEX channel. Curves 1–9, correspond to 107-100 copies and the negative control, respectively; (C) Probit regression analysis of the duplex real-time RAA assay results, the detection limit of the duplex real-time RAA assay at 95% reliability (73.67 copies/reaction) is marked with a red number

Analytical sensitivity of the duplex real-time RAA assay

Tenfold serial dilutions of PCV3 and PCV4 standard plasmids (107 to 100 copies per reaction) were used as amplification templates under the optimized system to determine the sensitivity of the duplex real-time RAA assay. Both PCV3 and PCV4 amplification curves showed a concentration dependent effect in fluorescence intensity across the tested concentration range of templates. Compared with the negative control, a discernible and consistent amplification signal could still be detected when the PCV3 and PCV4 standard plasmids concentration reached at 102 copies/reaction (Fig. 6B). To ascertain a more accurate limit of detection (LOD), probit regression analysis was conducted on the duplex real-time RAA assay (Fig. 6C). And results further revealed an LOD of 73.67 copies/reaction for PCV3 and PCV4 with 95% confidence interval, demonstrating the reliability of sensitivity of the duplex real-time RAA assay.

Analysis of repeatability and reproducibility

The repeatability and reproducibility of the PCV3 and PCV4 duplex real-time RAA assay were assessed by determining the CVs of the intra- and inter- batch assays for three different concentrations (104, 103, and 102 copies per reaction) of standard plasmids (pUC57-PCV3 and pUC57-PCV4). The intra-batch assay was tested in triplicate within one run and the inter-batch assay was determined in three independent runs performed within 3 separate days. The results indicated that the intra-batch assay CVs of PCV3 and PCV4 ranged from 2.08 to 3.91%, while the inter-batch assay CVs ranged from 3.33 to 4.97% (Table 3). The CV values for all measurements was below 5%, demonstrating that the duplex real-time RAA assay had dependable repeatability and reproducibility for PCV3 and PCV4 detection.

Table 3 Repeatability and reproducibility of duplex real-time RAA assay for PCV3 and PCV4 detection

Detection of clinical samples and artificially spiked samples

55 clinical samples (including blood, kidney, lymph node, muscle, environmental samples) preserved in our laboratory were respectively detected by PCV3 qPCR assay (according to local standards in China) and PCV4 qPCR assay reported before [6]. However, all samples displayed negative results by PCV4 qPCR. The diluted standard plasmid (different concentrations of pUC57-PCV4) was added into 3 PCV3-positive samples and 5 artificially spiked samples respectively, thus total 60 samples were used to evaluate the clinical performance of the duplex real-time RAA assay. All 8 samples contained pUC57-PCV4 were tested positively by both duplex real-time RAA and PCV4 qPCR assay. PCV3 positive/negative samples of duplex real-time RAA and PCV4 qPCR assay were 26/34 and 25/35, respectively. Additionally, 3 samples were detected positive for both PCV3 and PCV4 by duplex real-time RAA assay, which was consistent with the expected results. The comparison results between duplex real-time RAA and qPCR assays are shown in Table 4 Compared with the results of PCV3 and PCV4 qPCR assay respectively, the kappa value of duplex real-time RAA was 0.966 (p < 0.001) and 1(p < 0.001), and the overall agreement was 98.33% (59/60) and 100% (60/60). These results indicated that duplex real-time RAA assay for PCV3 and PCV4 detection could be applied in clinical practice.

Table 4 Comparison of the duplex real-time RAA with the qPCR assays on clinical samples

Discussion

Since their discovery, PCV3 and PCV4 have been confirmed in both domestic pigs and wild boar populations across various global regions. Accumulating evidence indicates that these emerging circoviruses have gradually shown detrimental effects on sows and piglets, and with similar clinical symptoms, making it difficult to distinguish between them clinically. In addition, co-infections with PCVs have been documented in a substantial number of case reports involving PCV3 and PCV4. For PCV2, there are commercial vaccines and numerous detection methods currently. While, the unclear pathogenesis of PCV3 and PCV4 poses significant challenges to prevent and control diseases associated with these viruses. Thus, the development of an accurate, rapid and simultaneous detection method for PCV3 and PCV4 is imperative for advancing research, diagnostics studies and epidemiological investigation of these emerging circoviruses. The RAA technique has garnered significant attention in the molecular diagnostics research community due to its distinctive advantages, including speed, isothermal reactivity and simple procedure [17, 18]. The establishment of duplex RAA approach in the detection of foodborne pathogenic bacteria [26], mycoplasma, chlamydia [27] and virus [25, 28] had been increasingly reported. However, a duplex real-time RAA approach for simultaneous detection of circoviruses is yet to be reported. Based on the above, our study is designed to develop a duplex real-time RAA assay capable of concurrently detecting PCV3 and PCV4 in clinical samples, addressing the need for a rapid and comprehensive diagnostic approach.

Establishing a rapid and accurate duplex RAA assay relies on the specific primer pairs and probes. The design of these elements necessitates rigorous screening and experimental validation. Unlike conventional qPCR, RAA primers and probes adhere to distinct design principles, such as a longer base length, which limited the use of existing primer and probe design software. For multiple RAA reactions, it is necessary to design multiple primer pairs and probes [19, 22, 26]. This process must account for the increased potential for binding between primers or probes and must circumvent the formation of undesirable structures, such as primer dimers, to ensure assay specificity and efficiency. In this study, we respectively conducted comparative analysis by aligning the genomic sequences of 30 strains (including genotype 3a, 3b and 3c) of PCV3 [29, 30] and 29 strains (including genotype 4a, 4b and 4c) of PCV4 [31] with low homology to identify the highly conserved regions. Utilizing alignment results and with the assist of software and BLAST tool, we developed 5 primer pairs (2 for PCV3 and 3 for PCV4) and 4 RAA probes (2 for PCV3 and 2 for PCV4) according to the principles (Figs. 1, 2 and 3; Table 1), and generated 6 sets of primer-probe combinations. Through the screening experiment, we identified the optimal primer-probe combination that demonstrated the best amplification performances and ensured no cross-interference among primers or between probes (Fig. 4). Subsequent research could aim to develop and refine specialized software for RAA primers and probes, thereby enhancing the efficiency of RAA establishment. Next, the concentration of optimal primers and probes were also variable factors in the duplex RAA reaction system, which might complicate reaction system. To enhance amplification efficiency and prevent mutual suppression of reactions due to concentration disparities, the concentrations of RAA primer pairs and probes were optimized (Fig. 5; Table 2). This optimization ensures a more reliable and consistent reaction outcome, which is crucial for our duplex real-time RAA assay.

In the context of co-infections with PCV2, PCV3 and PCV4 in pigs, several studies have reported the development of duplex and multiplex qPCR assays [32,33,34]. As a predominant diagnostic technique for simultaneous amplification and detection of multiple pathogen sequences, qPCR typically requires a minimum of 1 h to complete the amplification process. In contrast, our duplex real-time RAA assay achieves the same within 20 min at 39 °C, significantly reducing the reaction duration and improving the detection efficiency. Recently, the singleplex RAA technique has shown progress in the rapid detection of PCV3 and PCV4, enabling minimum detection limits ranging from 1 to 50 copies/µL of standard plasmid DNA [13, 16, 35, 36]. Here, our duplex real-time RAA assay exhibited both PCV3 and PCV4 detection limits of 73.67 copies per reaction at 95% probability, equivalent to 29.47 copies/µL of standard plasmid DNA (Fig. 6B and C). This data is within the established range of above singleplex RAA assays, thereby confirming the comparable sensitivity for the detection of two emerging circoviruses. Specificity of our assay was also performed good, and no cross reactivity with common swine pathogens including PCV2 (Fig. 6A). The coefficients of variation (CV) of the intra- and inter-group for this method ranged from 2.08 to 3.91% and 3.33–4.97%, providing good repeatability and reproducibility (Table 3). The above results suggested that newly developed duplex RAA assay could provide a relevant tool for more rapid and accurate diagnosis of PCV3 and PCV4. Additionally, our duplex real-time RAA assay was further applied for a pilot study of clinical samples (including blood, kidney, lymph node, muscle, environmental samples) and artificially spiked samples (due to undetected PCV4 in Zhejiang province). The comparative analysis of detection results between duplex real-time RAA assay and qPCR assay of PCV3 and PCV4 indicated that only one sample tested PCV3 positive by qPCR but negative in duplex real-time RAA assay, and 3 samples were detected positive for both PCV3 and PCV4 by duplex RAA assay as expected (Table 4). These findings demonstrated that our assay not only effectively detected co-infections but also produced detection results that are comparable to those of qPCR, rendering it suitable for preliminary clinical detecting of PCV3 and PCV4.

Several studies have shown the suitability of combining RAA technology with visualization techniques for diagnostics for resource-limited laboratories and on-site screening [26, 37, 38]. This combination yields rapid results without the need for sophisticated laboratory equipment and is cost-effective. The integration of PCV3 and PCV4 duplex real-time RAA assay with complementary product analysis technologies is essential for overcoming current technical challenges, representing an ongoing research direction of us. In this vein, the development of portable, user-friendly device equipped with thermostatic control, real-time fluorescence detection and multiple fluorescent channels is a priority. This device is essential for the field application of duplex real-time RAA assays, ensuring robust and reliable diagnostic capabilities. Furthermore, lateral flow test strip technology, recognized for its rapid and portable diagnostic potential [20, 28, 38], is another avenue we are exploring for integration with duplex RAA assay of PCV3 and PCV4. This integration aims to facilitate the visualization of diagnostic outcomes, enhancing the accessibility and immediacy of detection results for on-site applications.

Conclusion

A duplex real-time RAA assay capable of simultaneous and rapid detection of PCV3 and PCV4 was successfully established in this study. The detection operation is simple and could efficiently complete at 39 °C for a duration of 20 min. This method exhibits good specificity, sensitivity, repeatability, and reproducibility, ensuring it as a reliable tool for accurate diagnosis and epidemiological studies for PCV3 and PCV4. Furthermore, clinical performance of the developed assay demonstrates significant potential in differential diagnosis of coinfections and preliminary clinical detection.

Data availability

No datasets were generated or analysed during the current study.

Abbreviations

PCV3:

Porcine circoviruses 3

PCV4:

Porcine circoviruses 4

RAA:

Recombinase Aided Amplification

CV:

Coefficients of variation

CRESS:

DNA Eukaryotic circular Rep-encoding ssDNA

PCVs:

Porcine circoviruses

PCVAD:

Porcine circovirus associated disease

PDNS:

Porcine dermatitis and nephropathy syndrome

PCR:

Polymerase chain reaction

qPCR:

Quantitative polymerase chain reaction

ddPCR:

Droplet digital PCR

LAMP:

Loop-mediated isothermal amplification

LFD:

Lateral flow strip

CRISPR:

Clustered regularly interspaced short palindromic repeats

SSB:

Single-stranded DNA-binding protein

FMDV:

Foot and mouth disease virus

CSFV:

Classical swine fever virus

PRV:

Pseudorabies virus

PPV:

Porcine parvovirus

PCV2:

Porcine circoviruses 2

PEDV:

Porcine epidemic diarrhea virus

TGEV:

Transmissible Gastroenteritis Virus

PRRSV:

Porcine reproductive and respiratory syndrome virus

PDCoV:

Porcine deltacoronavirus

LOD:

Limit of detection

CT:

Cycle threshold

References

  1. Opriessnig T, Karuppannan AK, Castro A, Xiao CT. Porcine circoviruses: current status, knowledge gaps and challenges. Virus Res. 2020;286:198044. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.virusres.2020.198044.

    Article  CAS  PubMed  Google Scholar 

  2. Niu G, Chen S, Li X, Zhang L, Ren L. Advances in crosstalk between porcine circoviruses and host. Viruses. 2022;14(7):1419. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/v14071419.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Tischer I, Mields W, Wolff D, Vagt M, Griem W. Studies on epidemiology and pathogenicity of porcine circovirus. Arch Virol. 1986;91(3–4):271–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/BF01314286.

    Article  CAS  PubMed  Google Scholar 

  4. Meng XJ. Porcine circovirus type 2 (PCV2): pathogenesis and interaction with the immune system. Annu Rev Anim Biosci. 2013;1:43–64. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev-animal-031412-103720.

    Article  CAS  PubMed  Google Scholar 

  5. Palinski R, Pineyro P, Shang P, Yuan F, Guo R, Fang Y, et al. A Novel Porcine Circovirus distantly related to known circoviruses is Associated with Porcine Dermatitis and Nephropathy Syndrome and Reproductive failure. J Virol. 2017;91(1):e01879–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/JVI.01879-16.

    Article  CAS  PubMed  Google Scholar 

  6. Zhang HH, Hu WQ, Li JY, Liu TN, Zhou JY, Opriessnig T, et al. Novel circovirus species identified in farmed pigs designated as porcine circovirus 4, Hunan province, China. Transbound Emerg Dis. 2020;67(3):1057–61. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/tbed.13446.

    Article  CAS  PubMed  Google Scholar 

  7. Kroeger M, Temeeyasen G, Pineyro PE. Five years of porcine circovirus 3: what have we learned about the clinical disease, immune pathogenesis, and diagnosis. Virus Res. 2022;314:198764. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.virusres.2022.198764.

    Article  CAS  PubMed  Google Scholar 

  8. Wang D, Mai J, Yang Y, Xiao CT, Wang N. Current knowledge on epidemiology and evolution of novel porcine circovirus 4. Vet Res. 2022;53(1):38. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13567-022-01053-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Xu T, Chen XM, Fu Y, Ai Y, Wang DM, Wei ZY, et al. Cross-species transmission of an emerging porcine circovirus (PCV4): first molecular detection and retrospective investigation in dairy cows. Vet Microbiol. 2022;273:109528. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.vetmic.2022.109528.

    Article  CAS  PubMed  Google Scholar 

  10. Zhang LH, Wang TX, Fu PF, Zhao YY, Li HX, Wang DM, et al. First Molecular Detection and Genetic Analysis of a Novel Porcine Circovirus (Porcine Circovirus 4) in dogs in the World. Microbiol Spectr. 2023;11(2):e0433322. https://doiorg.publicaciones.saludcastillayleon.es/10.1128/spectrum.04333-22.

    Article  CAS  PubMed  Google Scholar 

  11. Xu T, Deng LS, Jian ZJ, Xu L, Li FQ, Lai SY, et al. First report on identification and genomic analysis of a novel porcine circovirus (porcine circovirus 4) in cats. Front Microbiol. 2023;14:1258484. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2023.1258484.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Sirisereewan C, Thanawongnuwech R, Kedkovid R. Current understanding of the pathogenesis of Porcine Circovirus 3. Pathogens. 2022;11(1):64. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/pathogens11010064.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Li Y, Zhao Y, Li C, Yang K, Li Z, Shang W, et al. Rapid detection of porcine circovirus type 4 via multienzyme isothermal rapid amplification. Front Vet Sci. 2022;9:949172. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fvets.2022.949172.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kim DY, Kim HR, Park JH, Kwon NY, Kim JM, Kim JK, et al. Detection of a novel porcine circovirus 4 in Korean pig herds using a loop-mediated isothermal amplification assay. J Virol Methods. 2022;299:114350. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jviromet.2021.114350.

    Article  CAS  PubMed  Google Scholar 

  15. Liu Y, Zhang X, Han X, Liu J, Yao L. Development of a droplet digital PCR method for detection of porcine circovirus 4. BMC Vet Res. 2023;19(1):129. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12917-023-03690-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wang J, Zhu X, Yin D, Cai C, Liu H, Yang Y, et al. Rapid and Easy-Read Porcine Circovirus Type 4 detection with CRISPR-Cas13a-Based lateral Flow Strip. Microorganisms. 2023;11(2):354. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/microorganisms11020354.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Deng H, Gao Z. Bioanalytical applications of isothermal nucleic acid amplification techniques. Anal Chim Acta. 2015;853:30–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.aca.2014.09.037.

    Article  CAS  PubMed  Google Scholar 

  18. Munawar MA. Critical insight into recombinase polymerase amplification technology. Expert Rev Mol Diagn. 2022;22(7):725–37. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/14737159.2022.2109964.

    Article  CAS  PubMed  Google Scholar 

  19. Tu F, Yang X, Xu S, Chen D, Zhou L, Ge X, et al. Development of a fluorescent probe-based real-time reverse transcription recombinase-aided amplification assay for the rapid detection of classical swine fever virus. Transbound Emerg Dis. 2021;68(4):2017–27. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/tbed.13849.

    Article  CAS  PubMed  Google Scholar 

  20. Zeng J, Wang W, Zhou L, Ge X, Han J, Guo X, et al. A nucleic acid detection assay combining reverse transcription recombinase-aided amplification with a lateral flow dipstick for the rapid visual detection of porcine deltacoronavirus. Virulence. 2022;13(1):1471–85. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/21505594.2022.2116157.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Guo Y, Wang Y, Fan Z, Zhao X, Bergmann SM, Dong H, et al. Establishment and evaluation of qPCR and real-time recombinase-aided amplification assays for detection of largemouth bass ranavirus. J Fish Dis. 2022;45(7):1033–43. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/jfd.13627.

    Article  CAS  PubMed  Google Scholar 

  22. Cui H, Tu F, Zhang C, Zhang C, Zhao K, Liu J, et al. Real-time reverse transcription recombinase-aided amplification assay for Rapid amplification of the N gene of SARS-CoV-2. Int J Mol Sci. 2022;23(23):15269. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms232315269.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Li X, Zhu S, Zhang X, Ren Y, He J, Zhou J, et al. Advances in the application of recombinase-aided amplification combined with CRISPR-Cas technology in quick detection of pathogenic microbes. Front Bioeng Biotechnol. 2023;11:1215466. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fbioe.2023.1215466.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Ma B, Li J, Chen K, Yu X, Sun C, Zhang M. Multiplex recombinase polymerase amplification assay for the simultaneous detection of three foodborne pathogens in Seafood. Foods. 2020;9(3):278. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/foods9030278.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Li G, Wu M, Li J, Cai W, Xie Y, Si G, et al. Rapid detection of porcine deltacoronavirus and porcine epidemic diarrhea virus using the duplex recombinase polymerase amplification method. J Virol Methods. 2021;292:114096. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jviromet.2021.114096.

    Article  CAS  PubMed  Google Scholar 

  26. Xiang S, Zhang H, Cha X, Lin Y, Shang Y. A new duplex recombinase polymerase amplification (D-RPA) method for the simultaneous and Rapid Detection of Shigella and Bacillus cereus in Food. Foods. 2023;12(9):1889. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/foods12091889.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Nie MZ, Zhang RQ, Zhao MC, Tan H, Hu YX, Fan GH, et al. Development of a duplex recombinase-aided amplification assay for direct detection of Mycoplasma pneumoniae and Chlamydia trachomatis in clinical samples. J Microbiol Methods. 2022;198:106504. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.mimet.2022.106504.

    Article  CAS  PubMed  Google Scholar 

  28. Li XN, Shen XX, Li MH, Qi JJ, Wang RH, Duan QX, et al. Applicability of duplex real time and lateral flow strip reverse-transcription recombinase aided amplification assays for the detection of Enterovirus 71 and Coxsackievirus A16. Virol J. 2019;16(1):166. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12985-019-1264-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Arruda B, Pineyro P, Derscheid R, Hause B, Byers E, Dion K, et al. PCV3-associated disease in the United States swine herd. Emerg Microbes Infect. 2019;8(1):684–98. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/22221751.2019.1613176.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Pan Y, Qiu S, Chen R, Zhang T, Liang L, Wang M, et al. Molecular detection and phylogenetic analysis of porcine circovirus type 3 in tibetan pigs on the Qinghai-Tibet Plateau of China. Virol J. 2022;19(1):64. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12985-022-01792-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Xu T, You D, Wu F, Zhu L, Sun XG, Lai SY, et al. First molecular detection and genetic analysis of porcine circovirus 4 in the Southwest of China during 2021–2022. Front Microbiol. 2022;13:1052533. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2022.1052533.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Hou CY, Xu T, Zhang LH, Cui JT, Zhang YH, Li XS, et al. Simultaneous detection and differentiation of porcine circovirus 3 and 4 using a SYBR Green capital I, Ukrainian-based duplex quantitative PCR assay. J Virol Methods. 2021;293:114152. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jviromet.2021.114152.

    Article  CAS  PubMed  Google Scholar 

  33. Zou J, Liu H, Chen J, Zhang J, Li X, Long Y, et al. Development of a TaqMan-Probe-based Multiplex Real-Time PCR for the Simultaneous Detection of Porcine Circovirus 2, 3, and 4 in East China from 2020 to 2022. Vet Sci. 2022;10(1):29. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/vetsci10010029.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Mi Y, Huang D, Zhuo Y, Li M, Yue J, Zhong H, et al. Development and preliminary application of a quadruplex real-time PCR assay for differential detection of porcine circovirus 1–4 in Chengdu, China. Front Vet Sci. 2024;11:1337461. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fvets.2024.1337461.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Li Y, Yu Z, Jiao S, Liu Y, Ni H, Wang Y. Development of a recombinase-aided amplification assay for rapid and sensitive detection of porcine circovirus 3. J Virol Methods. 2020;282:113904. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jviromet.2020.113904.

    Article  CAS  PubMed  Google Scholar 

  36. Jiang G, Yang X, Li Z, Mao J, Zeng P, Wang D, et al. Recombinant polymerase amplification coupled with CRISPR/Cas12a Detection System for Rapid Visual Detection of Porcine Circovirus 3. Anim (Basel). 2024;14(17):2527. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ani14172527.

    Article  Google Scholar 

  37. Ngoc LTN, Lee YC. Current trends in RNA virus detection via Nucleic Acid Isothermal amplification-based platforms. Biosens (Basel). 2024;14(2):97. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/bios14020097.

    Article  CAS  Google Scholar 

  38. Chen K, Zhang J, Wang S, Yi Z, Fu Y. Duplex recombinase aided amplification-lateral flow dipstick assay for rapid distinction of Mycobacterium tuberculosis and Mycobacterium avium complex. Front Cell Infect Microbiol. 2024;14:1454096. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fcimb.2024.1454096.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

This research was supported by the Zhejiang Provincial Natural Science Foundation of China (No. LQ23C18001), National Natural Science Foundation of China (No. 32172819), Zhejiang Province Science and Technology Cooperation Project of “Three Rural and Nine Parties” (No. 2023SNJF059) and the High-level Talents Special Support Plan of Zhejiang Province (No. 2021R52041).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: R.S, R.X and L.Z; Methodology: R.S, S.S and H.L; Software: R.S and Y.S; Validation: Y.W, Y.Y and R.S; Formal Analysis: R.S and Y.Y; Investigation, R.S and Y.Y; Resources: R.X, W.F and X.L; Data Curation: R.X and L.Z; Writing - Original Draft Preparation: R.S; Writing - Review & Editing: Y.Y, R.X and L.Z; Visualization: R.S and H.L; Supervision: Y.Y, R.X and L.Z; Project Administration: Y.Y, R.X and L.Z; Funding Acquisition: R.S, R.X, W.F and X.L. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yongle Yang, Ronghui Xie or Lingyan Zhao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, R., Liu, H., Sun, S. et al. Development of a duplex real-time recombinase aided amplification assay for the simultaneous and rapid detection of PCV3 and PCV4. Virol J 22, 23 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12985-025-02625-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12985-025-02625-w

Keywords