ORIGINAL RESEARCH article

Front. Microbiol., 03 July 2025

Sec. Infectious Agents and Disease

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1635733

Droplet digital RT-PCR method for SARS-CoV-2 variants detection in clinical and wastewater samples

  • Zhejiang Key Laboratory of Public Health Detection and Pathogenesis Research, Department of Microbiology, Zhejiang Provincial Center for Disease Control and Prevention (Zhejiang CDC), Hangzhou, China

Article metrics

View details

2

Citations

2,5k

Views

623

Downloads

Abstract

Objective:

To establish a sensitive, specific, and precise quantitative detection method for SARS-CoV-2 variants using droplet digital RT-PCR (RT-ddPCR).

Methods:

Dual primer-probe sets targeting the SARS-CoV-2 nucleocapsid (N) and spike (S) genes were designed. The annealing temperature for RT-ddPCR was optimized using a gradient PCR system. The sensitivity, defined as the limit of detection (LOD), was determined by serially diluting SARS-CoV-2 RNA. The specificity of the RT-ddPCR assay was evaluated using SARS-CoV-2 variants and common respiratory viruses. Precision and repeatability were assessed by quantitatively repeating the detection on serial dilutions of SARS-CoV-2 RNA. Additionally, the results of RT-ddPCR for clinical and environmental wastewater samples were compared with those from RT-qPCR.

Results:

The optimal annealing temperature was 53.5°C. The LOD for the N and S genes of the original SARS-CoV-2 strain was 4.26 (95% CI: 3.12–9.89) and 3.87 (95% CI: 2.77–7.75) copies/reaction. The Delta strain exhibited LODs of 4.65 (N gene, 95% CI: 3.28–9.64) and 6.12 (S gene, 95% CI: 4.33–15.59) copies/reaction. The Omicron showed 4.07 (N gene, 95% CI: 3.11–6.26) and 4.58 (S gene, 95% CI: 3.43–7.40) copies/reaction. Importantly, the RT-ddPCR assay was repeatable with a coefficient of variation of less than 10% when RNA concentrations of SARS-CoV-2 were between 73.50 and 7,500 copies/reaction. The high specificity of the RT-ddPCR assay was demonstrated by its ability to correctly detect the thirty SARS-CoV-2 variants, while not other common respiratory viruses. For 148 clinical pharyngeal swab specimens, the positive rate for both RT-ddPCR and RT-qPCR was 86.49%, and a coincidence rate of 98.65% and a Kappa value of 0.94. Quantitative comparison of RT-ddPCR and RT-qPCR in 50 wastewater samples with low viral load, RT-ddPCR assay detected 50 positives for dual gene targets (N and S genes), whereas RT-qPCR assay only 21 exhibited concurrent positivity for dual gene targets, while 25 showed S gene detection, and 4 were negative for dual gene targets, suggesting our RT-ddPCR assay enabled absolute quantification of SARS-CoV-2 variants with low viral load.

Conclusion:

The RT-ddPCR assay developed in this study can be used for SARS-CoV-2 variants detection and quantitative analysis of clinical and environmental samples.

1 Introduction

As of May 2025, more than 777 million confirmed cases of coronavirus disease 2019 (COVID-19), caused by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), were reported (WHO, 2025), thus substantially affecting human health and global economic growth. SARS-CoV-2, a member of the β-coronavirus family, has an envelope and a non-segmented, positive-sense single-stranded RNA genome. It shares 79 and 50% genome sequence similarity with SARS-CoV and MERS-CoV, respectively (Lu et al., 2020; Rabaan et al., 2020). According to the World Health Organization (WHO), real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) is the most frequently used method for diagnosing SARS-CoV-2 infection. However, the RT-qPCR assay has limited sensitivity, is vulnerable to variables such as standard curves, and is prone to false negatives in samples with low viral loads, thus hindering the prevention and management of outbreaks (Safiabadi Tali et al., 2021).

Droplet digital RT-PCR (RT-ddPCR) is based on the principle of limiting dilution PCR, with the PCR reaction mixture is uniformly partitioned into tens of thousands of independent micro-droplets. Consequently, some micro-droplets contain one or more template copies, whereas others lack the template. Each micro-droplet is subjected to PCR independently. A significant increase in fluorescence signal is detected for droplets containing template nucleic acid, whereas droplets lacking template nucleic acid maintain the background fluorescence intensity. Micro-droplets were dichotomously classified as positive or negative based on fluorescence signal thresholds (Xu et al., 2023). Finally, on the basis of the Poisson distribution, the number of positive micro-droplets is converted into a nucleic acid copy number, thus enabling absolute quantification of the target nucleic acid (Kojabad et al., 2021). Compared with RT-qPCR, RT-ddPCR exhibits higher sensitivity for samples with low viral load and stronger detection specificity. The RT-ddPCR method has been applied in the quantitative analysis of viruses such as hepatitis B virus, human immunodeficiency virus, Zika virus, enterovirus, parechovirus, and herpes simplex virus types 1 and 2 (Urso et al., 2016; Hui et al., 2018; Hayashi et al., 2022; Zhu et al., 2022).

The objective of this study was to establish an RT-ddPCR method for SARS-CoV-2 variants detection with high efficiency, specificity, and sensitivity. This method aims to improve detection accuracy in clinical and wastewater settings, decrease the occurrence of false-negative results, mitigate potential transmission risks, and contribute to more effective diagnostic strategies for SARS-CoV-2 variants infections (Ishak et al., 2021).

2 Materials and methods

2.1 Strains and samples

The original strain of SARS-CoV-2 and variants used for sensitivity evaluation were obtained from the laboratory of the Zhejiang Provincial Center for Disease Control and Prevention (Zhejiang CDC), China. The original strain, SARS-CoV-2/E6/WGF/2020/ZJ8, had a titer of 3.76 × 106 TCID50/mL; the Delta strain, SARS-CoV-2/Vero/LXG/2021/ZJ28, had a titer of 3.73 × 105 TCID50/mL; and the Omicron strain, SARS-CoV-2/E6/Gabriol/2022/ZJ60, had a titer of 3.16 × 105 TCID50/mL.

Seven common respiratory viruses and SARS-CoV-2 pseudovirus quantification reference material (Fantasiabio, Zhejiang, China; RFKSS001) were used to evaluate the specificity of RT-ddPCR. Seven common respiratory viruses were obtained from Zhejiang CDC, including influenza A (H1N1)pdm 09, Victoria lineage of influenza B virus, respiratory syncytial virus subtype A, human parainfluenza virus type III, adenovirus type 7, human coronavirus OC43, and human coronavirus 229E. The original strain and thirty SARS-CoV-2 variants, including the three variants of concern (VOCs) (Alpha, Delta, and Omicron) were used for specificity evaluation. The Omicron variants comprised 28 sublineages: BA.1.1, BA.2.12.1, BA.2.3, XBB.1.5.4, XBB.1.9.2, EG.5.1, EG.5.1.1, HK.3, XBB.1.16, FU.1, XBB.1.22, BA.4.1, BA.5.2, BF.7, BA.5.2.48, DY.2, BQ.1.1, BF.7.14, JN.1, JN.1.4.5, LB.1.2, JN.1.16, KP.2, KP.3.1.1, JN.1.67.1, XDV.1, XDV.1.5, and NB.1.

For quantitative comparison between RT-ddPCR and RT-qPCR, a total of 148 nasopharyngeal samples collected from fever clinics between June 2022 and December 2024 at three hospitals: Yiwu Central Hospital, Hangzhou First People’s Hospital, and Children’s Hospital, Zhejiang University School of Medicine, China. All samples were collected in 3 mL universal viral transport medium and transported to the laboratory immediately or storage at −80°C until use. Additionally, 50 environmental wastewater samples were obtained from the inlet of the Qianjiang sewage-treatment plant in Xiaoshan District, Hangzhou City, Zhejiang Province, China. The initial volume of each wastewater sample was 5 L.

2.2 Wastewater preparation

The samples were transported to the laboratory at 4°C and concentrated through hyperfiltration. Firstly, 400 mL of the 5 L untreated wastewater was centrifuged at 2,500 × g for 20 min at 4°C after blending. The liquid supernatant was then placed in an aseptic bottle and concentrated with a tubular ultrafiltration membrane connected to an ultrafiltration device. Next, the viruses retained on the ultrafiltration membrane were eluted with 2 mL of 3% beef extract solution. The eluates were stored at −80°C until use.

2.3 RNA extraction

RNeasy Mini Kit (Qiagen, Hilden, Germany; 74,104) was used to extract RNA according to the manufacturer’s instructions. A total of 50 μL viral RNA was extracted from 200 μL clinical and concentrated wastewater samples. The RNA was stored at −80°C until use.

2.4 Primer and probe sets

SARS-CoV-2 specific primer-probe sets targeting the N and S genes were developed for RT-ddPCR and RT-qPCR are provided in Table 1. The primer-probe sets were synthesized by Shanghai Sangon Company (Shanghai, China). All primers and probes were prepared at a concentration of 20 μmol/L, and primer-probe mix of the N and S genes were prepared in working solutions with a volume of 2:2:1.

Table 1

Target geneTypeSequence (5′-3′)
N geneForward primerACATTGGCACCCGCAATCC
Reverse primerGCTTGACTGCCGCCTCTGCT
ProbeFAM-5’-CGTGCTACAACTTCCTCAAGGAACA-3’-BHQ1
S geneForward primerTTGATCACAGGCAGACTTCAAAGT
Reverse primerAGCTCTGATTTCTGCAGCTCTAATT
ProbeVIC-5’-TGCAGACATATGTGACTCA-3’-BHQ1

SARS-CoV-2 primer-probes used in this study.

2.5 RT-ddPCR for SARS-CoV-2 quantification

Absolute quantification of SARS-CoV-2 RNA was performed with the One-Step RT-ddPCR Advanced Kit for Probes (Bio-Rad, Hercules, USA; 1,864,021). The reaction mixture volume of 20 μL comprised 5 μL Supermix, 2 μL primer-probe mix, 2 μL reverse transcriptase, 1 μL 300 mM dithiothreitol (DTT), 8 μL nuclease-free water, and 2 μL RNA template. According to the manufacturer’s guidelines, the reaction mixture was used to produce droplets with a Bio-Rad Auto Droplet Generator (Bio-Rad, USA). Thermal cycling for all RT-ddPCR assays was performed with a T100 Thermal Cycler (Bio-Rad, USA) under the following amplification conditions: 45°C for 60 min; 95°C for 10 min; 40 cycles of 95°C for 30 s and 52°C to 60°C for 1 min; 98°C for 10 min; and storage at 4°C. Droplet signals were read in different channels using the QX200 Droplet Reader (Bio-Rad, USA), FAM channel for the N gene, and VIC channel for the S gene. Data were considered valid if the total number of droplets in each tube was ≥ 10,000. A sample was considered positive if the number of droplets exceeded three and negative if the number of droplets was three or fewer. Results are expressed as copies per reaction (copies/reaction).

2.6 Optimization of annealing temperature and DTT concentration for RT-ddPCR

We used the original strain to optimize the RT-ddPCR annealing temperature and DTT concentration. The annealing temperature of RT-ddPCR was optimized by testing eight temperatures (60.0, 59.4, 58.4, 56.9, 55.1, 53.5, 52.5, and 52.0°C) with the described method. The optimal annealing temperature was determined on the basis of the signal discrimination and nucleic acid copy number of each reaction.

In PCR systems, DTT is frequently employed as a protein reducing agent to preserve the sulfhydryl group of cysteine in proteins in a reduced state and to safeguard PCR reaction enzymes. Subsequently, with the optimal annealing temperature, a 300 mM DTT concentration in the RT-ddPCR reaction system was optimized. Two concentrations, 0.5 μL and 1 μL of 300 mM DTT, were tested in four repeated experiments, and the results were analyzed to determine the optimal reaction conditions.

2.7 RT-qPCR reaction system

An AgPath-ID™ One-step RT-PCR Kit (Thermo Fisher, Carlsbad, CA, USA; AM1005) was used to prepare a 25 μL RT-qPCR reaction system comprising 12.5 μL 2 × RT-PCR buffer, 1.5 μL primer-probe mix, 1 μL reverse transcriptase, 8 μL nuclease-free water, and 2 μL RNA template. The RT-qPCR reaction conditions were as follows: 45°C for 10 min; 95°C for 10 min; and 40 cycles of 94°C for 15 s and 53.5°C for 35 s. The cycle threshold (CT) value was derived from the amplification of the N gene and S gene of SARS-CoV-2 variants using the ABI 7500 Real Time PCR System (Applied Biosystems, USA).

2.8 Evaluation of RT-ddPCR sensitivity, specificity, and repeatability

We assessed the sensitivity of the RT-ddPCR method by determining the LOD for three strains of SARS-CoV-2 (original, Delta, and Omicron). The initial nucleic acid concentrations of these strains were quantified with identical primer-probe systems in a preliminary experiment. Samples were initially diluted by a factor of 10 and subsequently 2-fold serially diluted for samples with low nucleic acid concentrations. The diluted nucleic acid samples underwent RT-ddPCR detection to determine the positive detection rate. LOD was defined as the concentration corresponding to the 95% confidence interval of nucleic acid copies. A lower LOD indicated higher detection sensitivity.

Additionally, influenza A (H1N1)pdm 09, Victoria lineage of influenza B virus, respiratory syncytial virus subtype A, parainfluenza virus type III, adenovirus type 7, human coronavirus OC43, human coronavirus 229E, and SARS-CoV-2 pseudovirus quantification reference material were used to evaluate the specificity of the RT-ddPCR method developed in this study.

SARS-CoV-2 pseudovirus nucleic acid standards were utilized to extract nucleic acids and evaluate the repeatability of RT-ddPCR. The concentrations of positive quantitative reference materials ranged from 1 × 102 to 1 × 106 copies/mL, accompanied by a negative control group. Each concentration group underwent 16 tests under optimal RT-ddPCR conditions. The resulting data were used to calculate the mean, standard deviation, and coefficient of variation to assess the repeatability of the experiment. A coefficient of variation less than 10% indicated good repeatability.

2.9 Comparison of RT-ddPCR and RT-qPCR for quantitative detection of SARS-CoV-2 variants

The original strain and thirty variants of SARS-CoV-2 were selected: Alpha (B.1.1.7), Delta (B.1.617.2), BA.1.1, BA.2.12.1, BA.2.3, XBB.1.5.4, XBB.1.9.2, EG.5.1, EG.5.1.1, HK.3, XBB.1.16, FU.1, XBB.1.22, BA.4.1, BA.5.2, BF.7, BA.5.2.48, DY.2, BQ.1.1, BF.7.14, JN.1, JN.1.4.5, LB.1.2, JN.1.16, KP.2, KP.3.1.1, JN.1.67.1, XDV.1, XDV.1.5, and NB.1. The RT-ddPCR and RT-qPCR methods were used for quantitative detection, and reaction systems were as previously described.

2.10 Comparison of RT-ddPCR and RT-qPCR for quantitative detection in clinical and wastewater samples

A total of 148 clinical specimens from three hospitals (Yiwu Central Hospital, Hangzhou First People’s Hospital, and Children’s Hospital, Zhejiang University School of Medicine, China) were analyzed with the RT-ddPCR method described previously. Additionally, 50 wastewater samples from a sewage treatment plant in Hangzhou underwent the same extraction and amplification methods. The results obtained via RT-ddPCR subsequently compared with those obtained via RT-qPCR.

2.11 Statistical analysis

RT-ddPCR data were analyzed using QuantSoft 1.7.4.0917 software (Bio-Rad). The LOD, the coincidence rate, and Kappa value were calculated through SPSS 26.0 software, with a threshold for positive detection set at 95%. Statistical significance was set at p < 0.05.

3 Results

3.1 Determination of the optimal annealing temperature and DTT concentration for RT-ddPCR

Eight gradient annealing temperatures (60.0, 59.4, 58.4, 56.9, 55.1, 53.5, 52.5, and 52.0°C) were tested, revealing that as the temperature decreases, the RT-ddPCR amplification efficiency gradually increases (see Table 2). Improved separation between positive and negative droplets was notably observed within the range of 53.5–52.0°C (see Figure 1). To ensure PCR specificity and optimize amplification efficiency, we identified 53.5°C as the optimal annealing temperature.

Table 2

Temperature (°C)N gene (copies/reaction)S gene (copies/reaction)
60.055211
59.450424
58.450268
56.9534410
55.1520576
53.5530534
52.5504570
52.0670746

Copy numbers of SARS-CoV-2 N and S genes at various annealing temperatures, detected by RT-ddPCR.

Figure 1

The results comparing the optimal volume of 300 mM DTT in the reaction system, between 0.5 μL and 1 μL, were replicated four times each (see Figure 2). Using of 0.5 μL 300 mM DTT in a 20 μL reaction system enabled clear differentiation between positive and negative droplets. Consequently, 0.5 μL is considered the optimal concentration of 300 mM DTT. Therefore, the final optimized RT-ddPCR reaction system consisted of 20 μL comprising 5 μL SuperMix, 2 μL primer-probe mix, 2 μL reverse transcriptase, 0.5 μL DTT (300 mM), 8.5 μL nuclease-free water, and 2 μL RNA template.

Figure 2

3.2 Evaluation of RT-ddPCR sensitivity, specificity, and repeatability

Three strains of SARS-CoV-2, original, Delta, and Omicron, were initially diluted to concentrations of 105, 104, and 104 copies/reaction, respectively, and followed by 2-fold serial dilution for RT-ddPCR analysis. Each dilution underwent 16 repetitions to achieve a 95% positive detection rate. Probit regression analysis was performed on the average copy numbers and positive detection rates of the N and S genes across various dilutions for each strain. The fitted curves demonstrated that the p values from Probit analysis for all three strains were <0.05, indicating statistical significance. The LOD for the original strain was 4.26 copies/reaction (95% CI: 3.12–9.89) for the N gene and 3.87 copies/reaction (95% CI: 2.77–7.75) for the S gene. For the Delta strain, the LOD was 4.65 copies/reaction (95% CI: 3.28–9.64) for the N gene and 6.12 copies/reaction (95% CI: 4.33–15.59) for the S gene. For the Omicron strain, the LOD was 4.07 copies/reaction (95% CI: 3.11–6.26) for the N gene and 4.58 copies/reaction (95% CI: 3.43–7.40) for the S gene. Detailed results can be found in Table 3.

Table 3

StrainDilutionTarget geneAverage copy number (copies/reaction)Positive rate (%)LOD (copies/reaction)95% CI (copies/reaction)
original10–5-16N gene0.5837.504.263.12–9.89
10–5-82.2362.50
10–5-42.9987.50
10–5-28.70100.00
10−528.75100.00
10–5-16S gene0.9150.003.872.77–7.75
10–5-83.7993.75
10–5-45.79100.00
10–5-212.26100.00
10−550.50100.00
Delta10–4-32N gene0.4431.254.653.28–9.64
10–4-163.6887.50
10–4-89.99100.00
10–4-425.63100.00
10–4-276.50100.00
10−4189.50100.00
10–4-32S gene1.6150.006.124.33–15.59
10–4-165.8393.75
10–4-812.05100.00
10–4-426.38100.00
10–4-288.88100.00
10−4215.88100.00
Omicron10–4-16N gene0.166.254.073.11–6.26
10–4-80.2112.50
10–4-41.2862.50
10–4-23.8087.50
10−411.84100.00
10–4-16S gene0.156.254.583.43–7.40
10–4-80.5431.25
10–4-41.4043.75
10–4-24.5993.75
10−413.83100.00

Sensitivity results for three strains of SARS-CoV-2 (original, Delta, Omicron).

The specificity of RT-ddPCR was evaluated with influenza A (H1N1)pdm 09, Victoria lineage of influenza B virus, respiratory syncytial virus subtype A, parainfluenza virus type III, adenovirus type 7, human coronavirus OC43, human coronavirus 229E, and SARS-CoV-2 pseudovirus quantification reference material (see Figure 3). Positive droplets were observed for only the SARS-CoV-2 pseudovirus quantification reference material, whereas no positive detection was observed for the other viruses. This funding indicates that the RT-ddPCR method exhibits specificity for SARS-CoV-2 detection without cross-reacting with other respiratory viruses.

Figure 3

The repeatability evaluation results of SARS-CoV-2 pseudovirus quantitative reference materials using the RT-ddPCR method were presented in Table 4. These results were compared with the target values specified in the instructions: L1: 1 × 102 copies/mL, L2: 1 × 103 copies/mL, L3: 1 × 104 copies/mL, L4: 1 × 105 copies/mL, and L5: 1 × 106 copies/mL. The errors ranged from 1.38 to 11.38%, thus demonstrating consistency between the quantitative results and the target values of the reference materials. For nucleic acid copy numbers >73.50 copies/reaction, the coefficients of variation, calculated from the mean and standard deviation, were all below 10%, thereby indicating that the RT-ddPCR had excellent repeatability within the range of 73.50–7,500 copies/reaction, thus ensuring the reliability experimental outcomes.

Table 4

NumberTarget geneRepetitionsAverage copy number (copies/reaction)Standard deviationCoefficient of variation (%)
L1N gene161.760.7542.69
S gene1.260.7660.51
L2N gene168.212.8434.61
S gene8.002.3529.39
L3N gene1673.507.4710.16
S gene76.887.559.82
L4N gene16769.2544.525.79
S gene78944.775.67
L5N gene167,535461.876.13
S gene7,090434.576.13

RT-ddPCR detection results of SARS-CoV-2 pseudovirus quantitative reference material.

3.3 Comparison of RT-ddPCR and RT-qPCR detection among SARS-CoV-2 variants

The RT-ddPCR results for the original strain and 30 variants of SARS-CoV-2 revealed copy numbers ranging from 3 to 46,010 copies/reaction for the N gene, and from 1.4 to 54,300 copies/reaction for the S gene. Correspondingly, the N gene CT values in RT-qPCR ranged from 37.69 to 23.55, whereas the S gene CT values ranged from 37.47 to 22.51. These findings are detailed in Table 5 and illustrated in Figure 4. The consistency between the results of the two methods across the SARS-CoV-2 variants suggested that RT-ddPCR was well-suited to detection of SARS-CoV-2 variants.

Table 5

LineageSublineageN geneS gene
CT valueCopy number (copies/reaction)CT valueCopy number (copies/reaction)
Original24.2746,01023.4752,700
AlphaB.1.1.730.9423230.11322
DeltaB.1.617.228.712,33527.902,870
OmicronBA.1BA.1.131.2823129.91448
BA.2BA.2.12.135.6319.934.7450
BA.2.334.9318833.18403
JN.128.6855527.81600
JN.1.4.534.191734.356.8
LB.1.231.196830.5570
KP.235.304.633.767.2
KP.3.1.131.845530.9661
JN.1.1629.2251829.61210
JN.1.67.132.841334.701.4
BA.4BA.4.126.056,56025.716,830
BA.5BA.5.228.161,45628.261,102
BA.5.2.4831.329630.6596
BF.729.9522329.43234
BF.7.1432.545731.8744
DY.229.9923329.50234
BQ.1.131.0512630.64101
XBBXBB.1.5.424.1823,09023.4934,620
XBB.1.9.237.69337.476.8
EG.5.128.244,77026.6410,110
EG.5.1.133.0222231.98330
HK.323.5539,52022.5154,300
XBB.1.1628.9671327.731,165
FU.129.5249428.231,130
XBB.1.2231.031,05429.871780
XDVXDV.132.642331.9330
XDV.1.528.6260329.00281
NB.129.6218628.76280

Detection results of SARS-CoV-2 variants by RT-ddPCR and RT-qPCR.

Figure 4

3.4 Quantitative detection results of clinical and wastewater samples by RT-ddPCR

We compared the established RT-ddPCR method with RT-qPCR in the detection of 148 clinical specimens and 50 wastewater samples. The RT-ddPCR analysis consistently detected more than 10,000 effective droplets per sample, thus ensuring data reliability. Among the 148 clinical specimens, both RT-ddPCR and RT-qPCR identified 128 positives for dual gene targets and 20 negatives, and the positive rate were 86.49% between two methods. There were 146 samples with consistent results, with a concordance rate of 98.65% and a Kappa value of 0.94. At 95% confidence interval, the sensitivity was 99.22% (95.02–99.96%) and the specificity was 95.00% (73.06–99.74%). RT-ddPCR quantification results ranged from 1.8 to 61,600 copies/reaction for the N gene and from 1.4 to 60,600 copies/reaction for the S gene. Corresponding RT-qPCR CT values ranged from 18.90 to 38.72 for the N gene and 18.51 to 38.86 for the S gene (see Supplementary Table 1).

In the analysis of low-concentration wastewater samples, RT-ddPCR detected 50 positives, whereas RT-qPCR identified 21 were positive for both the N and S genes, 25 were positive for the S gene, and 4 were negative. RT-ddPCR quantification results ranged from 0.6 to 163 copies/reaction for the N gene, and from 0.6 to 136 copies/reaction for the S gene. In contrast, in RT-qPCR, the CT values ranged from 37.92 to 32.01 for the N gene, and from 38.16 to 30.33 for the S gene. Overall, RT-ddPCR had slightly superior quantification performance to that of RT-qPCR in detecting low viral loads (see Table 6).

Table 6

SampleTarget geneCT valueCopy number (copies/reaction)SampleTarget geneCT valueCopy number (copies/reaction)
1N gene35.368.326N gene-1.2
S gene34.4711.3S gene38.010.6
2N gene36.325.327N gene-3.6
S gene35.560.6S gene34.044.7
3N gene37.921.928N gene-3.7
S gene37.861.2S gene35.026.2
4N gene36.674.829N gene-5.9
S gene35.561.8S gene34.606.3
5N gene33.345730N gene-2.3
S gene32.6153S gene36.551.7
6N gene36.542.431N gene-3.5
S gene35.503.6S gene34.104.6
7N gene36.443.532N gene-4.6
S gene35.847.5S gene34.284.6
8N gene36.852233N gene-5.8
S gene34.1712.8S gene33.965.8
9N gene35.727.434N gene-3.6
S gene35.224.3S gene33.136
10N gene34.294235N gene-4.3
S gene32.4342S gene34.818.7
11N gene34.047036N gene-0.6
S gene32.2754S gene35.970.6
12N gene32.0816337N gene-5
S gene30.50136S gene34.258.2
13N gene35.3112238N gene-3.2
S gene32.63105S gene34.761.3
14N gene34.4710439N gene-1.9
S gene32.3897S gene35.461.9
15N gene33.689940N gene-3.4
S gene31.9776S gene35.824
16N gene37.115.841N gene-8.1
S gene36.012.2S gene34.002.9
17N gene34.4618.242N gene-8
S gene33.2132S gene34.478
18N gene37.146.243N gene-3.2
S gene34.763.6S gene34.533.2
19N gene36.421.644N gene-1.3
S gene37.131.2S gene36.240.6
20N gene32.822645N gene-1.2
S gene31.8954S gene34.091.2
21N gene32.019046N gene-4.3
S gene30.33126S gene34.373.7
22N gene-1.247N gene-0.7
S gene38.162.6S gene-2
23N gene-4.148N gene-0.6
S gene32.587.1S gene-3.9
24N gene-1.349N gene-0.6
S gene35.722.6S gene-1.2
25N gene-350N gene-3.9
S gene36.023S gene-2.8

Detection of SARS-CoV-2 in wastewater samples by RT-ddPCR and RT-qPCR.

4 Discussion

Since the onset of the COVID-19 pandemic, SARS-CoV-2 has been mutating and spreading globally, thus posing a considerable public health threat. Traditional RT-qPCR methods often fail to accurately quantify viral copies, and may potentially fail to detect low concentrations of nucleic acids. Developing a more precise, efficient, and quantifiable detection method is imperative for effective epidemic prevention and control.

We developed an RT-ddPCR method that enhances accuracy by quantifying SARS-CoV-2 variants with dual primers and probes targeting the N and S genes. SARS-CoV-2 N protein has evolutionary conservation, which was used for great diagnostic marker (Eltayeb et al., 2024). Notably, the SARS-CoV-2 S glycoprotein is composed of two subunits, S1 and S2. In the prefusion state, the S1 subunit mediates binding to the host cell receptor angiotensin converting enzyme 2 (ACE2), while the S2 subunit drives viral envelope fusion with the host membrane (Wrapp et al., 2020). Mutations in SARS-CoV-2 variants predominantly localize to the receptor-binding domain (RBD) of the S1 subunit (Walls et al., 2020). As the primary target for neutralizing antibodies, the S protein elicits potent humoral immunity, with the RBD harboring the dominant neutralizing epitopes responsible for over 90% of neutralizing activity (Jackson et al., 2022). In this study, our forward primer (5′-3′, 2986–3009), reverse primer (5′-3′, 3036–3060), and probe (5′-3′, 3011–3029) of S gene were targeting the central helices (CH) in relatively high conservation region, demonstrating the universal detection efficacy across SARS-CoV-2 variants. Highly conserved primer-probe sets are essential for the detection of SARS-CoV-2 variants.

Unlike traditional methods, RT-ddPCR does not rely on standard curves, thereby enabling direct quantification of nucleic acid concentrations and demonstrating superior sensitivity (Huggett et al., 2015; Park et al., 2021). Tao has highlighted that ddPCR has greater sensitivity than RT-qPCR in detecting low viral loads, and has benefits of requiring minimal nucleic acid amounts, without a need for repeated sampling or extensive reagents (Suo et al., 2020). Herein, we conducted sensitivity experiments on three SARS-CoV-2 strains and achieved a lowest detection limit of <6.12 copies/reaction with a 95% positive detection rate, thus validating our method’s high sensitivity in the detection of samples with low viral load. In terms of specificity, our RT-ddPCR method showed no cross-reactivity with common respiratory viruses such as influenza A and B viruses, respiratory syncytial virus, parainfluenza virus, and other coronaviruses. Additionally, repeatability studies using SARS-CoV-2 pseudovirus RNA revealed a coefficient of variation <10% at nucleic acid concentrations ranging from 73.50 to 7,500 copies/reaction, in agreement with the highly reproducible results observed in various laboratory tests by Whale et al. (2017). Furthermore, RT-ddPCR accurately quantified SARS-CoV-2 variants, thereby underscoring its excellent specificity. Compared with RT-qPCR, RT-ddPCR provided advantages in precise quantification and accuracy in detecting clinical and environmental samples.

Despite its advantages, RT-ddPCR had several limitations. Quantification accuracy may be compromised at high target concentrations (≥105 copies/reaction) due to saturation effects, reducing the confidence in detection (Quan et al., 2018). In such cases, pre-dilution experiments should be conducted before detection. Additionally, RT-ddPCR necessitated higher standards for instruments, equipment, and experimental personnel, thus contributing to its lower adoption than other methods (Kojabad et al., 2021).

In conclusion, this study established an RT-ddPCR detection method that accurately quantified low concentrations of SARS-CoV-2 variants, exhibiting robust specificity, high sensitivity, and excellent repeatability. This method was well-suited to early clinical detection of SARS-CoV-2 infections and tracing viral presence in the environment. We believe that this method may be applied to provide valuable insights in clinical diagnosis and treatment, thus ultimately mitigating the risk and effects of viral transmission.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Ethics Committee of Chinese Academy of Medical Sciences & Peking Union Medical College, Chinese Academy of Medical Sciences & Peking Union Medical College (CAMS&PUMC-IEC-2023-001). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

FW: Writing – review & editing, Writing – original draft, Investigation, Formal analysis, Data curation. YS: Data curation, Writing – review & editing, Formal analysis. LG: Writing – review & editing, Investigation. LS: Writing – review & editing, Investigation. BZ: Writing – review & editing, Formal analysis, Data curation. XL: Writing – review & editing, Investigation. YC: Writing – review & editing, Conceptualization. WS: Writing – review & editing, Formal analysis, Data curation. HM: Conceptualization, Writing – review & editing. YZ: Conceptualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences (2022-I2M-CoV19-006); National Key Research and Development Program of China (2021YFC2301200); Zhejiang Science and Technology Plan for Disease Prevention and Control (2025JK002).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1635733/full#supplementary-material

References

  • 1

    EltayebA.Al-SarrajF.AlharbiM.AlbiheyriR.MattarE.Abu ZeidI. M.et al. (2024). Overview of the SARS-CoV-2 nucleocapsid protein. Int. J. Biol. Macromol.260:129523. doi: 10.1016/j.ijbiomac.2024.129523

  • 2

    HayashiS.IsogawaM.KawashimaK.ItoK.ChuaypenN.MorineY.et al. (2022). Droplet digital PCR assay provides intrahepatic HBV cccDNA quantification tool for clinical application. Sci. Rep.12:2133. doi: 10.1038/s41598-022-05882-9

  • 3

    HuggettJ. F.CowenS.FoyC. A. (2015). Considerations for digital PCR as an accurate molecular diagnostic tool. Clin. Chem.61, 79–88. doi: 10.1373/clinchem.2014.221366

  • 4

    HuiY.WuZ.QinZ.ZhuL.LiangJ.LiX.et al. (2018). Micro-droplet digital polymerase chain reaction and real-time quantitative polymerase chain reaction technologies provide highly sensitive and accurate detection of Zika virus. Virol. Sin.33, 270–277. doi: 10.1007/s12250-018-0037-y

  • 5

    IshakA.AlRawashdehM. M.EsagianS. M.NikasI. P. (2021). Diagnostic, prognostic, and therapeutic value of droplet digital PCR (ddPCR) in COVID-19 patients: a systematic review. J. Clin. Med.10:5712. doi: 10.3390/jcm10235712

  • 6

    JacksonC. B.FarzanM.ChenB.ChoeH. (2022). Mechanisms of SARS-CoV-2 entry into cells. Nat. Rev. Mol. Cell Biol.23, 3–20. doi: 10.1038/s41580-021-00418-x

  • 7

    KojabadA. A.FarzanehpourM.GalehH. E. G.DorostkarR.JafarpourA.BolandianM.et al. (2021). Droplet digital PCR of viral DNA/RNA, current progress, challenges, and future perspectives. J. Med. Virol.93, 4182–4197. doi: 10.1002/jmv.26846

  • 8

    LuR. J.ZhaoX.LiJ.NiuP.YangB.WuH.et al. (2020). Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet395, 565–574. doi: 10.1016/S0140-6736(20)30251-8

  • 9

    ParkC.LeeJ.HassanZ. U.KuK. B.KimS.-J.KimH. G.et al. (2021). Comparison of digital PCR and quantitative PCR with various SARS-CoV-2 primer-probe sets. J. Microbiol. Biotechnol.31, 358–367. doi: 10.4014/jmb.2009.09006

  • 10

    QuanP.-L.SauzadeM.BrouzesE. (2018). dPCR: a technology review. Sensors18:1271. doi: 10.3390/s18041271

  • 11

    RabaanA. A.Al-AhmedS. H.HaqueS.SahR.TiwariR.MalikY. S.et al. (2020). SARS-CoV-2, SARS-CoV, and MERS-COV: a comparative overview. Infez. Med.28, 174–184

  • 12

    Safiabadi TaliS. H.LeBlancJ. J.SadiqZ.OyewunmiO. D.CamargoC.NikpourB.et al. (2021). Tools and techniques for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)/COVID-19 detection. Clin. Microbiol. Rev.34:e00228-20. doi: 10.1128/CMR.00228-20

  • 13

    SuoT.LiuX.FengJ.GuoM.HuW.GuoD.et al. (2020). ddPCR: a more accurate tool for SARS-CoV-2 detection in low viral load specimens. Emerg Microbes Infect9, 1259–1268. doi: 10.1080/22221751.2020.1772678

  • 14

    UrsoC.PierucciF.SollaiM.ArviaR.MassiD.ZakrzewskaK. (2016). Detection of Merkel cell polyomavirus and human papillomavirus DNA in porocarcinoma. J. Clin. Virol.78, 71–73. doi: 10.1016/j.jcv.2016.03.008

  • 15

    WallsA. C.ParkY.-J.TortoriciM. A.WallA.McGuireA. T.VeeslerD. (2020). Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell181, 281–292.e6. doi: 10.1016/j.cell.2020.02.058

  • 16

    WhaleA. S.DevonshireA. S.Karlin-NeumannG.ReganJ.JavierL.CowenS.et al. (2017). International interlaboratory digital PCR study demonstrating high reproducibility for the measurement of a rare sequence variant. Anal. Chem.89, 1724–1733. doi: 10.1021/acs.analchem.6b03980

  • 17

    WHO COVID-19 dashboard (2025). Available online at: https://data.who.int/dashboards/covid19/cases (Accessed May 8, 2025).

  • 18

    WrappD.WangN.CorbettK. S.GoldsmithJ. A.HsiehC.-L.AbionaO.et al. (2020). Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science367, 1260–1263. doi: 10.1126/science.abb2507

  • 19

    XuD.ZhangW.LiH.LiN.LinJ.-M. (2023). Advances in droplet digital polymerase chain reaction on microfluidic chips. Lab Chip23, 1258–1278. doi: 10.1039/d2lc00814a

  • 20

    ZhuX.LiuP.LuL.ZhongH.XuM.JiaR.et al. (2022). Development of a multiplex droplet digital PCR assay for detection of enterovirus, parechovirus, herpes simplex virus 1 and 2 simultaneously for diagnosis of viral CNS infections. Virol. J.19:70. doi: 10.1186/s12985-022-01798-y

Summary

Keywords

droplet digital RT-PCR, SARS-CoV-2 variants, clinical, wastewater, sensitivity, specificity, repeatability

Citation

Wang F, Sun Y, Gong L, Su L, Zhou B, Lou X, Chen Y, Shi W, Mao H and Zhang Y (2025) Droplet digital RT-PCR method for SARS-CoV-2 variants detection in clinical and wastewater samples. Front. Microbiol. 16:1635733. doi: 10.3389/fmicb.2025.1635733

Received

27 May 2025

Accepted

18 June 2025

Published

03 July 2025

Volume

16 - 2025

Edited by

Xia Cai, Fudan University, China

Reviewed by

Li Guo, Institute of Pathogen Biology, China

Zheng Shen, Zhejiang University, China

Updates

Copyright

*Correspondence: Haiyan Mao, ; Yanjun Zhang,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics