Skip to main content

ORIGINAL RESEARCH article

Front. Public Health, 12 January 2024
Sec. Infectious Diseases: Epidemiology and Prevention
This article is part of the Research Topic Women in Science: Infectious Diseases: Epidemiology and Prevention 2023 View all 15 articles

Ready for new waves: optimizing SARS-CoV-2 variants monitoring in pooled samples with droplet digital PCR

Antonella Pacini,&#x;Antonella Pacini1,2Franco Paredes,&#x;Franco Paredes1,3Sofia Heckel,,Sofia Heckel1,2,3Guadalupe Ibarra,Guadalupe Ibarra1,3Maria Victoria Petreli,Maria Victoria Petreli1,3Marilina PerezMarilina Perez1Yanina Agnella,Yanina Agnella1,4Laura PiskulicLaura Piskulic5Maria Belen AllasiaMaria Belen Allasia5Luis CaprileLuis Caprile1Alejandro ColaneriAlejandro Colaneri6Juliana Sesma,,
Juliana Sesma1,2,7*
  • 1Molecular Biology Department, Hospital Provincial de Rosario, Rosario, Argentina
  • 2Instituto de Inmunología Clínica y Experimental de Rosario, CONICET, Rosario, Argentina
  • 3Facultad de Ciencias Bioquímicas y Farmacéuticas de Rosario, Universidad Nacional de Rosario, Rosario, Argentina
  • 4Facultad de Ciencias Veterinarias, Universidad Nacional de Rosario, Casilda, Argentina
  • 5Área Estadística y Procesamiento de Datos, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Rosario, Argentina
  • 6Consejo Nacional de Investigaciones Científicas y Técnicas, Rosario, Argentina
  • 7Facultad de Ciencias Médicas, Universidad Nacional de Rosario, Rosario, Argentina

Introduction: The declaration of the end of the Public Health Emergency for COVID-19 on May 11th, 2023, has shifted the global focus led by WHO and CDC towards monitoring the evolution of SARS-CoV-2. Augmenting these international endeavors with local initiatives becomes crucial to not only track the emergence of new variants but also to understand their spread. We present a cost-effective digital PCR-based pooled sample testing methodology tailored for early variant surveillance.

Methods: Using 1200 retrospective SARS-CoV-2 positive samples, either negative or positive for Delta or Omicron, we assessed the sensitivity and specificity of our detection strategy employing commercial TaqMan variant probes in a 1:9 ratio of variant-positive to variant-negative samples.

Results: The study achieved 100% sensitivity and 99% specificity in 10-sample pools, with an Area Under the Curve (AUC) exceeding 0.998 in ROC curves, using distinct commercial TaqMan variant probes.

Discussion: The employment of two separate TaqMan probes for both Delta and Omicron establishes dual validation routes, emphasizing the method’s robustness. Although we used known samples to model realistic emergence scenarios of the Delta and Omicron variants, our main objective is to demonstrate the versatility of this strategy to identify future variant appearances. The utilization of two divergent variants and distinct probes for each confirms the method’s independence from specific variants and probes. This flexibility ensures it can be tailored to recognize any subsequent variant emergence, given the availability of its sequence and a specific probe. Consequently, our approach stands as a robust tool for tracking and managing any new variant outbreak, reinforcing our global readiness against possible future SARS-CoV-2 waves.

Introduction

The dynamic nature of COVID-19 pandemic is characterized by the emergence and succession of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants, categorized as variants of concern (VOCs), variants of interest (VOIs), and, more recently, variants being monitored (VBM) (15).

As of this writing, the World Health Organization (WHO) recognizes Omicron (B.1.1.529 and descendant lineages) as the only VOC, while the Alpha (B1.1.7), Beta (B.1.351), Gamma (P1), and Delta (B.1.617.2) are considered VBM (1). The declaration of the end of the Public Health Emergency for COVID-19 on May 11th, 2023, has shifted the global focus led by WHO and CDC towards monitoring the evolution of SARS-CoV-2.

Traditional methods for tracking variants of SARS-CoV-2, such as Next-Generation Sequencing (NGS), provide definitive identification of mutations. NGS would be the method of choice when one seeks to monitor the emergence of new variants with unknown sequences. However, in the event of a wave of a new VOC, the timeliness of NGS falls short of the immediate demands for treatment decisions and for the surveillance required in the face of a quickly spreading virus (6, 7). Furthermore, the cost and infrastructure required for NGS are substantial, limiting its reach to countries with the necessary resources. In this context, more accessible and rapid methods would be needed. One such alternative is RT-PCR genotyping panels, which have enabled the quick identification of circulating SARS-CoV-2 variants (8, 9). However, these strategies primarily rely on individual testing, thereby escalating the time, labor, and costs involved.

Droplet digital PCR (ddPCR), an increasingly popular method for nucleic acid detection and quantification, offers an alternative approach (1012). This technique compartmentalizes DNA molecules into thousands of droplets, amplifying them simultaneously, leading to several advantages, including decreased inhibitor concentration, increased relative sample concentration, and improved accuracy. The benefits of ddPCR have led to its widespread adoption in clinical research for the diagnosis of various diseases (1316). Within our ddPCR process, two TaqMan probes with different fluorophores are utilized for amplification. These probes only differ in the sequence corresponding to the mutation of the variant under study. The probe with the HEX-fluorophore typically hybridizes specifically to the wild-type (WT) sequence, while the probe with the FAM fluorophore is specific for the sequence carrying the point mutation of interest. Our current digital PCR technology, featuring two channels of detection, supports the use of multiple probes, as evidenced in several studies (17, 18). Effective multiplexing is achieved by employing the same fluorophore for various probes, with distinct signal differentiation accomplished by adjusting the concentration of each probe. This allows for the separation of signal clusters in a two-dimensional plot, enhancing the versatility of digital PCR. Our system is designed with primers and probes that amplify a specific fragment size, typically between 60 to 100 base pairs. Thus, the size of the target marker does not affect the intensity of fluorescence in our assays. Rather, the fluorescence intensity within a positive droplet is predominantly influenced by the concentration of primers and probes used in the reaction. A higher concentration of these components generally results in stronger fluorescence signals, which is a crucial aspect for multiplexing applications. Therefore, by varying the concentrations of different probes, we can effectively multiplex without compromising the accuracy and efficiency of the assay.

In keeping with advancements in the field, newer versions of ddPCR technology now allow up to six channels to be read, significantly enhancing the analytical capabilities of ddPCR. This advancement enables the quantification of up to 12 targets in a single well, representing a substantial increase in throughput and efficiency for complex analyses, further broadening its application in research and diagnostics.

Following amplification, the fluorescence of each droplet is read, and the concentrations of mutant and WT molecules are determined using a Poisson distribution (19).

The utilization of pooling assays in diagnostic PCR offers significant advantages, including reduced costs, increased testing capacity, and shorter turnaround times, as multiple samples can be simultaneously analyzed in a single assay, providing an efficient and cost-effective approach to detect and monitor pathogens in diagnostics assays (2022).

In this paper, we present a systematic variant detection approach that we evaluated on a total of 1,200 retrospective samples collected between August 2021 and May 2022 using commercially TaqMan variant probes. We analyzed the Delta and Omicron variants in pools of 10 SARS-CoV-2 positive samples, achieving a sensitivity of 100% and a specificity of 99% compared to individual samples.

During the initial phase of the Delta variant’s emergence, we evaluated our platform using a double-entry chart combinatorial pool-testing strategy on 41 SARS-CoV-2 positive samples. This approach efficiently reduced the number of samples needing ddPCR testing from 41 to 13. Upon identifying positive intersections of rows and columns, we performed individual RT-PCR tests on the intersecting samples and were able to detect all three that carried the Delta variant. This led to notable savings in both reagents and time.

While our study focused on known, retrospective samples and targeted SARS-CoV-2 variants that may no longer be relevant for current testing, this was a necessary step to fine-tune the method. The true intent behind this work is to ensure that, in the future, this methodology can be applied to monitor the circulation of any newly emerging variant. The ddPCR pooling strategy is especially effective for detecting these SARS-CoV-2 variants when they first begin to circulate and their incidence is very low, allowing for pooled testing. In this landscape, both NGS and ddPCR techniques can complement each other. While NGS stands out in the identification of new mutations, ddPCR excels in tracking their distribution, making it an invaluable tool for rapid response to any rapidly evolving pathogens with pandemic potential.

Materials and methods

Sample collection

Nasopharyngeal samples were collected between August 2021 and May 2022, mixed with 2 mL saline solution and stored at 4°C until extraction. In Biosafety Level 2 containment, the virus was inactivated and then extracted automatically using the Applied Biosystems MagMAX kit and the Thermo Scientific KingFisher purification system.

Individual RT-qPCR test

RT-qPCR was performed in the clinical laboratory to detect the presence of SARS-CoV-2 RNA with the PerkinElmer SARS-CoV-2 RT-qPCR Reagent kit®, following manufacturer instructions. We received de-identified positive samples with a code. SARS-CoV-2 positive samples were genotyped by RT-PCR to detect the variants of concern (Delta and Omicron). For that, probes P681H.CCT.CAT, L452R.CTG.CGG, Q954H.CAA.CAT and P681R.CCT.CGT (ThermoFisher Scientific®) and CoV_B1429_L452R (Biorad®) were used. At the time of assay, the spike protein mutation P681R was used to identify and distinguish the Delta variant and L452R was specific for Delta and Epsilon. On the other hand, Q954H was specific for Omicron while P681H was specific for Omicron and Alpha variants.

Pooling preparation

To confirm that the samples employed for pool creation were devoid of the variants under study, we utilized samples collected prior to the emergence of these variants. We designated these samples, which tested positive for SARS-CoV-2 but negative for the specific variant in study, as “WT-SARS-CoV-2 samples.” The SARS-CoV-2 sample carrying the mutation in study was designated VARIANT-SARS-CoV-2 (Va-SARS-CoV-2) sample. For the construction of sample pools, we combined one Va-SARS-CoV-2 sample with nine WT-SARS-CoV-2 samples.

Primer/probe thermal gradient optimization

Optimus pool size and dilution were used as a template to assess the optimal annealing temperature of the individual primer/probes. The standard ddPCR cycling program was modified by replacing the annealing temperature step with a thermal gradient from 55°C to 60°C for 1 min extension time. These experiments were performed with P681H and Q954H from Thermo Fisher Scientific® and L452R from Thermo Fisher Scientific® and Biorad®.

Pooled ddPCR test

One-Step RT-ddPCR Advanced Kit for Probes (BioRad) was used for RT-ddPCR according to the manufacturer’s recommendations. 5,5 μL of pooled samples were added to 16,5 ul of master mix with variants/WT primers and probes according to the manufacturer’s recommendations. VOCs probes (P681H, L452R, Q954H, and P681R) were labeled with carboxyfluorescein (FAM), WT probes were labeled with HEX fluorophore. RT-ddPCR reactions were set up in 96-well ddPCR plates. Plates were covered with a pierceable foil heat seal, sealed using a PX1 PCR Plate Sealer (BioRad), mixed by vortexing, and centrifuged before droplet generation. Droplet generation was done in the QX200 AUTO DG (BioRad), which dispensed droplets into a new 96-well plate. Then, plates were heat-sealed with a pierceable foil and thermocycled on a C1000 touch thermal cycler (BioRad) as follows: 50°C for 50 min; 95°C for 10 min; 5 cycles of 94°C for 15 s and 55°C for 50s, 35 cycles of 94°C for 10 s and 55°C for 45 s, 98°C for 10 min, followed by a 25 min hold at 4°C for droplet stabilization. Droplets were then read on a QX200 Droplet Reader (BioRad) set up to read FAM and HEX channels and analyzed using the QuantaSoft Analysis Pro 1.0.596 software. The software performed the automated droplet count for each case, which was then manually reviewed to determine whether positive droplets fell outside gating parameters. Besides the cut-off determined and according to manufacturer instructions, a minimum number of accepted droplets is required to ensure optimal detection in the analysis: positive assays must have at least 6,000 accepted droplets, negative ones must have 10,000 or more (23).

Statistical analysis

The overall performance of the new technique was evaluated by calculating the area under the Receiver Operator Characteristics (ROC) curve and its 95% confidence interval using the method described by Hanley and McNeil (24). The best cut-off point was selected using the Youden Index (25), with a priority for the test to have a 100% sensitivity. The estimates of sensitivity and specificity were determined using a 95% confidence interval for previously established cut-off values.

Data were processed using R. Among other packages, pROC and OptimalCutpoints were utilized.

Analysis of the ddPCR data was performed with QuantaSoft analysis software v.1.7.4.0917 (Bio-Rad) to calculate the concentration of the targets.

Ethical approval

The study protocol was approved on December 10th, 2020 by the Bioethics Committee of Universidad Nacional de Rosario, Facultad de Ciencias Médicas (resolution: 3733/2022). Signed informed consent was not necessary as de-identified samples of SARS-CoV-2 tests by HOSPITAL PROVINCIAL DE ROSARIO were used. This study was conducted in accordance with the principles of the 1964 Declaration of Helsinki.

Results

Assay optimization: pool size and thermal gradient

While our previous work demonstrated the ability to detect a SARS-CoV-2 positive sample within a 34-sample pool (26), identifying specific Va-SARS-CoV-2 sequences among SARS-CoV-2 positive samples proves to be more complex. This complexity arises because, in the former case, we were seeking a specific sequence within a vast array of DNA without similar sequences. In contrast, now we face the added challenge of differentiating between the target sequence and the wild-type sequence, which differ by only a single nucleotide.

In our initial tests, we purposefully mixed one Delta sample (with a CT value of 27) with nine WT-SARS-CoV-2 samples. We intentionally chose a sample with a high CT value, which holds clinical relevance, to challenge the system’s capability in detecting variants at a low concentration (27). Assay conditions were then modified for optimization. As depicted in Figures 1AD the dilution of the pool with water (1/2, 1/4, and 1/30) and the reduction of the annealing temperature expanded the separation between the WT and mutant droplet clouds (Figures 2A,B). This increase in the separation of the clouds in a graph was already described by the manufacturer (28). The optimal condition was found to be a 1/30 pool dilution in water combined with an annealing temperature of 55°C.

Figure 1
www.frontiersin.org

Figure 1. Assay optimization: influence of sample dilution on assay performance. A Delta sample (CT = 27) was combined with nine WT-SARS-CoV-2 samples (CT values ranging from 14 to 20). The combined pool was then analyzed by ddPCR with the L452R probe (Catalogue # CVAAAAD). We examined the distribution of droplet populations across various dilution conditions with water: (A) undiluted, (B) 1/2 dilution, (C) 1/4 dilution, and (D) 1/30 dilution. The diagrams illustrate the variation in droplet distribution under each dilution condition.

Figure 2
www.frontiersin.org

Figure 2. Assay optimization: impact of thermal gradient on annealing temperature. A Delta sample (CT = 24) was combined with nine WT-SARS-CoV-2 samples (CT values ranging from 14 to 20). This pooled sample was then analyzed by ddPCR with the L452R probe (CoV_B1429_L452R Biorad®). Assays were conducted at three different annealing temperatures: 55°C, 58°C, and 60°C. (A) A 2-D plot displaying the distribution of the droplet populations at each of the three annealing temperatures. The smooth line represents 60°C, the dashed line represents 58°C, and the purple dotted line represents 55°C. (B) 1-D plots of the ddPCR reactions at different temperatures, with blue representing positive droplets for the variant, green for the WT, and gray for negative droplets.

Specificity and sensitivity assessment

Recognizing the competitive amplification environment where the mutant sequence is outnumbered by the predominance of WT sequence, it was crucial to establish conditions enhancing the method’s specificity. We assessed the specificity of our method in accordance with the guidelines of the European Pharmacopoeia (29). We obtained 1,000 WT SARS-CoV-2 samples that predate the Delta variants’ emergence. These samples were divided into 100 pools for analysis, and each pool was examined by ddPCR with the 4 probes named in Table 1. Each probe had to be run independently to ascertain their distinct cutoff points. The cut-off points were the maximum number of positive drops tolerated in a negative sample.

Table 1
www.frontiersin.org

Table 1. Details of sensitivity, specificity, and area under the curve (AUC) for all four probes used for SARS-CoV-2 variant detection.

Evaluating the sensitivity of our platform presented a unique challenge due to the need to detect a specific variant sequence amidst a majority of WT-SARS-CoV-2 samples. To accurately measure this sensitivity, we designed our tests according to the guidelines provided by the European Pharmacopoeia (27). Our sensitivity analysis was initiated with a Delta SARS-CoV-2 sample (CT = 27) combined with nine WT SARS-CoV-2 samples (with CT values ranging from 14 to 20). This high-CT Delta sample was purposefully chosen to test the system’s capacity to detect a low-concentration variant amid a higher-concentration of WT SARS-CoV-2 samples. To ensure the robustness of our findings, each measurement was repeated six times across four independent experiments, resulting in a total of 24 repetitions for each probe. This process was conducted separately for each of the two Delta variant-specific probes. A similar study design was applied for the two probes specific to the Omicron variant. Table 1 presents the results. The sensitivity at the best cut-off value was 100% for all probes, while the specificity was 100% for the Q954H and L452R probes, and 99% for the P681R and P681H probes.

The Receiver Operating Characteristic (ROC) curves were used to analyze the specificity and sensitivity of the four different probes (Figures 3AD). Table 1 shows that the areas under the curve (AUC) were 1 for the Q954H and L452R probes, and above 0.998 for the P681R and P681H probes.

Figure 3
www.frontiersin.org

Figure 3. Sensitivity and specificity assessment. Receiver operating characteristic (ROC) curves are illustrated for the following probes: (A) L452R, (B) Q954H, (C) P681R, and (D) P681H. These curves were derived from the analysis of pooled samples consisting of one Va-SARS-CoV-2 sample (CT = 27) and nine WT-SARS-CoV-2 samples (CT values ranging from 14 to 25).

In regards to the P681R and P681H probes, we prioritized having 100% sensitivity over specificity as we did not want to miss any positive sample. Conversely, if a pool tested false positive, this error could be eliminated by de-pooling and individually testing each sample via RT-PCR.

Clinical case study

We had the opportunity to evaluate our platform during the emergence of the Delta variant wave from September 17th to 22nd, 2021. At that time, the incidence of Delta cases in our country was still quite low, representing about less than 2% of total cases (30). We selected 41 leftover SARS-CoV-2 positive samples from that week and applied them to a double-entry chart combinatorial pool-testing strategy, also known as matrix or grid testing (20). This approach entails testing each sample twice in a grid format, which facilitates the cross-matching of results to accurately identify positive samples. The methodology is visually represented in Table 2. As the incidence was 2%, pools were assembled with either six or seven samples each.

Table 2
www.frontiersin.org

Table 2. Double-entry chart with analyzed samples for the real case.

Table 2 shows that pools B, E, 5, 6, and 7 returned positive results for the Delta variant (light gray). The cells located at the intersection of these columns and rows, indicated in dark gray, were designated as “putative positive samples.” We proceeded to “deconstruct” these five pools and individually re-tested the six dark gray samples using RT-PCR. This process resulted in the confirmation of samples 19–509, 19–510, and 22–626 as Delta variant positives.

It is important to highlight that this approach effectively reduced the number of samples requiring ddPCR testing from 41 down to 13 (plus an additional 6 samples individually tested by RT-PCR). Furthermore, this method ensured each sample was tested twice: once in column pools and once in row pools.

We employed this technology briefly. Within just one week, the prevalence of the delta variant in our country jumped from 2% to over 10%. Given this rapid increase and the subsequent high prevalence, pooling samples became impractical, as the majority were expected to test positive for the variant (30).

Conclusion

Despite the declaration of the end of the Public Health Emergency for COVID-19 on May 11th, 2023, national efforts to monitor the emergence and spread of new variants continue. While these strategies may vary from country to country, they generally involve genomic surveillance of a percentage of patients, travelers, and wastewater samples (1).

Although next-generation sequencing remains the primary method for detecting and tracing new SARS-CoV-2 variants, there are scenarios where its inherent challenges – such as cost, extended turnaround times, and the need for specialized personnel – could benefit from complementary methodologies (6, 7). Techniques like Sanger sequencing and RT-PCR genotyping using single-nucleotide polymorphism (SNP) can serve as supplementary tools, providing faster and more cost-effective solutions in certain contexts (8, 31, 32). Digital PCR is a sophistication of the traditional RT-PCR that has higher sensitivity to detect a low copy number mutation (33). Additionally, pooling samples with ddPCR offers the significant benefit of reducing reagent costs. The extent of these savings is influenced by several factors, including the positivity rate, which in turn guides the optimal group size, as well as the inherent cost of ddPCR testing. When comparing the costs, if ddPCR is approximately twice as expensive as RT-PCR, pooling 1,000 samples in groups of 10 with a positivity rate of 2% can result in reagent cost reductions of at least 50%.

Our study presents a robust and flexible platform for the detection of specific SARS-CoV-2 variants amidst a backdrop of 10 WT-SARS-CoV-2 samples. We have successfully adapted and optimized our earlier SARS-CoV-2 detection methodology (26) leveraging a ddPCR pooling strategy, to identify two different variants, Delta and Omicron.

A pivotal strength of our platform lies in its inherent adaptability. It is designed to be modifiable with minimal alterations needed to detect future emerging variants. Through the strategic selection of two distinct TaqMan probes for each variant, we established two independent methods for each variant detection, enhancing the reliability and versatility of our platform.

Our platform’s sensitivity and specificity were carefully evaluated following the European Pharmacopoeia. We demonstrated 100% sensitivity at the best cut-off value for all probes, with specificity ranging from 99% to 100% across the four probes used. The real-world efficacy of our platform was validated during the emergence of the Delta variant, accurately identifying Delta-positive samples within the tested pools.

In conclusion, we have developed a robust, sensitive, and flexible platform capable of accurately identifying specific SARS-CoV-2 variants within predominantly WT-SARS-CoV-2 sample pools. Moreover, the underlying principles and techniques of this platform hold significant potential for broader applications. It is not only adaptable for detecting variants from other viruses beyond SARS-CoV-2 but also versatile enough to work effectively across various sample matrices. Digital PCR’s ability to detect different types of viruses, including respiratory pathogens and others like Ebola and Dengue, in blood samples, has been validated (34, 35). For instance, our work on the detection of SARS-CoV-2 in saliva (26), and the effectiveness of digital PCR in other matrices such as urine and feces for SARS-CoV-2 (36) and bacteria in blood (34), underlines the flexibility of this technology in a variety of diagnostic scenarios.

The adaptability and precision of this platform may thus contribute profoundly not only to the current fight against COVID-19, but also to future endeavors in viral epidemiology and public health. By simply modifying the targeted genetic sequences, our methodology could serve as a versatile foundation for rapid response to new viral threats, enhancing our capabilities in global health surveillance and response.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the following link: https://doi.org/10.6084/m9.figshare.24915639.

Ethics statement

The studies involving humans were approved by COMITÉ DE BIOSEGURIDAD, Facultad de Ciencias Medicas, Universidad Nacional de Rosario. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from remnants from SARS-CoV-2 tests conducted at the Hospital Provincial de Rosario. Signed informed consent was waived as patient information was anonymized and de-identified prior to analysis. This study was carried out in accordance with the principles of the 1964 Declaration of Helsinki. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

AP: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. FP: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. SH: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. GI: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. MVP: Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft. MaP: Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. YA: Investigation, Methodology, Writing – review & editing. LP: Formal analysis, Writing – review & editing. MA: Formal analysis, Writing – review & editing. LC: Investigation, Writing – original draft. AC: Conceptualization, Funding acquisition, Investigation, Resources, Writing – review & editing. JS: Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the MINISTERIO DE CIENCIA, TECNOLOGIA E INNOVACION (grant number SF 06-COVID FEDERAL: IF-2020-37418385-APN-SSFCTEI#MCT, 2020). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

We would like to express our sincere gratitude to the Hospital Provincial de Rosario and its dedicated healthcare personnel for their invaluable support and cooperation throughout this study. Their tireless efforts in patient care and sample collection have been essential to the success of our research. We would also like to extend our appreciation to Dr. Oscar Bottasso and Dr. Javier Girardini, for their valuable contributions in reviewing and editing the manuscript. Their insightful feedback and expertise have greatly improved the clarity and quality of this publication.

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.

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.

References

1. WHO coronavirus (COVID-19) dashboard | WHO coronavirus (COVID-19) dashboard with vaccination data. Available at: https://covid19.who.int/ (Accessed Jun 7, 2023).

Google Scholar

2. Tracking SARS-CoV-2 variants. Available at: https://www.who.int/activities/tracking-SARS-CoV-2-variants. (Accessed Jun 7, 2023).

Google Scholar

3. Nabaes Jodar, MS, Torres, C, Mojsiejczuk, L, Acuña, D, Valinotto, LE, Goya, S, et al. The lambda variant in Argentina: analyzing the evolution and spread of SARS-CoV-2 lineage C.37. Viruses. (2023) 15:1382. doi: 10.3390/v15061382

PubMed Abstract | Crossref Full Text | Google Scholar

4. Torres, C, Nabaes Jodar, M, Acuña, D, Montaño, RMZ, Culasso, ACA, Amadio, AF, et al. Omicron waves in Argentina: dynamics of SARS-CoV-2 lineages BA.1, BA.2 and the emerging BA.2.12.1 and BA.4/BA.5. Viruses. (2023) 15:312. doi: 10.3390/v15020312

PubMed Abstract | Crossref Full Text | Google Scholar

5. Hijano, DR, Ferrolino, JA, Gu, Z, Brazelton, JN, Zhu, H, Suganda, S, et al. Digital PCR to measure SARS-CoV-2 RNA, variants, and outcomes in youth. J Pediatric Infect Dis Soc. 12:618–26. doi: 10.1093/jpids/piad101/7416854

Crossref Full Text | Google Scholar

6. Greninger, AL, Bard, JD, Colgrove, RC, Graf, EH, Hanson, KE, Hayden, MK, et al. Clinical and infection prevention applications of severe acute respiratory syndrome coronavirus 2 genotyping: an Infectious Diseases Society of America/American Society for Microbiology consensus review document. J Clin Microbiol. (2022) 60:e0165921–502. doi: 10.1128/JCM.01659-21

PubMed Abstract | Crossref Full Text | Google Scholar

7. Rosenthal, SH, Gerasimova, A, Ruiz-Vega, R, Livingston, K, Kagan, RM, Liu, Y, et al. Development and validation of a high throughput SARS-CoV-2 whole genome sequencing workflow in a clinical laboratory. Sci Rep. (2022) 12:2054. doi: 10.1038/s41598-022-06091-0

PubMed Abstract | Crossref Full Text | Google Scholar

8. Lai, E, Kennedy, EB, Lozach, J, Hayashibara, K, Davis-Turak, J, Becker, D, et al. A method for variant agnostic detection of SARS-CoV-2, rapid monitoring of circulating variants, and early detection of emergent variants such as omicron. J Clin Microbiol. (2022) 60:e0034222. doi: 10.1128/jcm.00342-22

PubMed Abstract | Crossref Full Text | Google Scholar

9. Perchetti, GA, Zhu, H, Mills, MG, Shrestha, L, Wagner, C, Bakhash, SM, et al. Specific allelic discrimination of N501Y and other SARS-CoV-2 mutations by ddPCR detects B.1.1.7 lineage in Washington state. J Med Virol. (2021) 93:5931–41. doi: 10.1002/jmv.27155

PubMed Abstract | Crossref Full Text | Google Scholar

10. Pernet, O, Weisenhaus, M, Stafylis, C, Williams, C, Campan, M, Pettersson, J, et al. SARS-CoV-2 viral variants can rapidly be identified for clinical decision making and population surveillance using a high-throughput digital droplet PCR assay. Sci Rep. (2023) 13:7612. doi: 10.1038/s41598-023-34188-7

PubMed Abstract | Crossref Full Text | Google Scholar

11. Stafylis, C, Pernet, O, Hernandez-Tamayo, C, Kovacs, A, Emerson, J, Ward, PM, et al. Surveillance of severe acute respiratory syndrome coronavirus 2 and variants using Digital Droplet polymerase chain reaction at a large university and healthcare system in California. Open Forum Infect Dis. (2023) 10:ofad147. doi: 10.1093/ofid/ofad147

PubMed Abstract | Crossref Full Text | Google Scholar

12. Hijano, DR, Gu, Z, Brazelton, J, Zhu, H, Suganda, S, Glasgow, HL, et al. Prediction of symptomatic severe acute respiratory syndrome coronavirus 2 infection by quantitative Digital polymerase chain reaction normalized to international units. Open Forum Infect Dis. (2022) 9:ofac490. doi: 10.1093/ofid/ofac490

PubMed Abstract | Crossref Full Text | Google Scholar

13. Li, H, Bai, R, Zhao, Z, Tao, L, Ma, M, Ji, Z, et al. Application of droplet digital PCR to detect the pathogens of infectious diseases. Biosci Rep. (2018) 38:BSR20181170. doi: 10.1042/BSR20181170

PubMed Abstract | Crossref Full Text | Google Scholar

14. Ruelle, J, Yfantis, V, Duquenne, A, and Goubau, P. Validation of an ultrasensitive digital droplet PCR assay for HIV-2 plasma RNA quantification. J Int AIDS Soc. (2014) 17:19675. doi: 10.7448/IAS.17.4.19675

PubMed Abstract | Crossref Full Text | Google Scholar

15. Sedlak, RH, Cook, L, Cheng, A, Magaret, A, and Jerome, KR. Clinical utility of droplet digital PCR for human cytomegalovirus. J Clin Microbiol. (2014) 52:2844–8. doi: 10.1128/JCM.00803-14

PubMed Abstract | Crossref Full Text | Google Scholar

16. Hijano, DR, De Cardenas, JB, Maron, G, Garner, CD, Ferrolino, JA, Dallas, RH, et al. Clinical correlation of influenza and respiratory syncytial virus load measured by digital PCR. PLoS One. (2019) 14:e0220908. doi: 10.1371/journal.pone.0220908

PubMed Abstract | Crossref Full Text | Google Scholar

17. Du, Y, Yan, Z, Song, K, Jin, J, Xiao, L, Sun, Z, et al. Development and evaluation of a multiplex droplet digital polymerase chain reaction method for simultaneous detection of five biothreat pathogens. Front Microbiol. (2022) 13:13. doi: 10.3389/fmicb.2022.970973

PubMed Abstract | Crossref Full Text | Google Scholar

18. Whale, AS, Huggett, JF, and Tzonev, S. Fundamentals of multiplexing with digital PCR. Biomol Detect Quantif. (2016) 10:15–23. doi: 10.1016/j.bdq.2016.05.002

PubMed Abstract | Crossref Full Text | Google Scholar

19. Suo, T, Liu, X, Feng, J, Guo, M, Hu, W, Guo, D, et al. ddPCR: a more accurate tool for SARS-CoV-2 detection in low viral load specimens. Emerg Microbes Infect. (2020) 9:1259–68. doi: 10.1080/22221751.2020.1772678

PubMed Abstract | Crossref Full Text | Google Scholar

20. Grobe, N, Cherif, A, Wang, X, Dong, Z, and Kotanko, P. Sample pooling: Burden or solution? Clin Microbiol Infect. (2021) 27:1212–20. doi: 10.1016/j.cmi.2021.04.007

PubMed Abstract | Crossref Full Text | Google Scholar

21. Cherif, A, Grobe, N, Wang, X, and Kotanko, P. Simulation of Pool testing to identify patients with coronavirus disease 2019 under conditions of limited test availability. JAMA Netw Open. 3:e2013075. doi: 10.1001/jamanetworkopen.2020.13075

PubMed Abstract | Crossref Full Text | Google Scholar

22. Zhuang, X, Zhao, Z, Feng, X, Lui, GCY, Chan, D, Lee, SS, et al. Integrating magnetic-bead-based sample extraction and molecular barcoding for the one-step pooled RT-qPCR assay of viral pathogens without retesting. Anal Chem. (2023) 95:6182–90. doi: 10.1021/acs.analchem.3c00885

PubMed Abstract | Crossref Full Text | Google Scholar

23. Xu, J, Kirtek, T, Xu, Y, Zheng, H, Yao, H, Ostman, E, et al. Digital Droplet PCR for SARS-CoV-2 resolves borderline cases. Am J Clin Pathol. (2021) 155:815–22. doi: 10.1093/ajcp/aqab041

PubMed Abstract | Crossref Full Text | Google Scholar

24. Hanley, JA, and McNeil, BJ. The meaning and use of the area under a receiver operating characteristic (ROC) Curve1. Radiology. (1982) 143:29–36. doi: 10.1148/radiology.143.1.7063747

Crossref Full Text | Google Scholar

25. Youden, WJ . Index for rating diagnostic tests. Cancer. (1950) 3:32–5. doi: 10.1002/1097-0142(1950)3:1<32::aid-cncr2820030106>3.0.co;2-3

Crossref Full Text | Google Scholar

26. Heckel, S, Pacini, A, Paredes, F, Petreli, MV, Perez, M, Adriani, N, et al. Practical considerations to establish a validated platform for pooled detection of SARS-CoV-2 by droplet digital PCR. PLoS One. (2022) 17:e0271860. doi: 10.1371/journal.pone.0271860

PubMed Abstract | Crossref Full Text | Google Scholar

27. Goswami, C, Sheldon, M, Bixby, C, Keddache, M, Bogdanowicz, A, Wang, Y, et al. Identification of SARS-CoV-2 variants using viral sequencing for the Centers for Disease Control and Prevention genomic surveillance program. BMC Infect Dis. (2022) 22:404. doi: 10.1186/s12879-022-07374-7

PubMed Abstract | Crossref Full Text | Google Scholar

28. Bio-rad Laboratories.Droplet Digital PCR Droplet Digital PCR Applications Guide. Bulletin. 6407.

Google Scholar

29. European pharmacopoeia (Ph. Eur.) 11th edition—European Directorate for the Quality of medicines & HealthCare. https://www.edqm.eu/en/web/edqm/european-pharmacopoeia-ph.-eur.-11th-edition (Accessed Jun 7, 2023).

Google Scholar

30. De, M, De, S, Nación, LA, Vizzotti, DC, and De Gabinete, J. Covid-19 Situación de nuevas variantes SARS-CoV-2 en Argentina, S35/2022. Informe Técnico Ministerio de Salud, Argentina Agosto (2022).

Google Scholar

31. Scott, L, Hsiao, NY, Moyo, S, Singh, L, Tegally, H, Dor, G, et al. Track Omicron's spread with molecular data. Science. (2021) 374:1454–5. Epub 2021 Dec 9. doi: 10.1126/science.abn4543

PubMed Abstract | Crossref Full Text | Google Scholar

32. Murillo, E, Palacio-Rua, K, Afanador-Ayala, C, García-Correa, JF, and Zuluaga, AF. Validation of a new strategy for the identification of SARS-CoV-2 variants by sequencing the spike gene by sanger. Enferm Infecc Microbiol Clin. (2023) 41:284–9. doi: 10.1016/j.eimce.2022.10.003

PubMed Abstract | Crossref Full Text | Google Scholar

33. Mills, MG, Hajian, P, Bakhash, SM, Xie, H, Mantzke, D, Zhu, H, et al. Rapid and accurate identification of SARS-CoV-2 omicron variants using droplet digital PCR (RT-ddPCR). J Clin Virol. (2022) 154:105218. doi: 10.1016/j.jcv.2022.105218

PubMed Abstract | Crossref Full Text | Google Scholar

34. Wouters, Y, Dalloyaux, D, Christenhusz, A, Roelofs, HMJ, Wertheim, HF, Bleeker-Rovers, CP, et al. Droplet digital polymerase chain reaction for rapid broad-spectrum detection of bloodstream infections. Microb Biotechnol. (2020) 13:657–68. doi: 10.1111/1751-7915.13491

PubMed Abstract | Crossref Full Text | Google Scholar

35. Mairiang, D, Songjaeng, A, Hansuealueang, P, Malila, Y, Lertsethtakarn, P, Silapong, S, et al. Application of one-step reverse transcription droplet digital pcr for dengue virus detection and quantification in clinical specimens. Diagnostics. (2021) 11:639. doi: 10.3390/diagnostics11040639

PubMed Abstract | Crossref Full Text | Google Scholar

36. Shinde, M, Lavania, M, Rawal, J, Chavan, N, and Shinde, P. Evaluation of droplet digital qRT-PCR (dd qRT-PCR) for quantification of SARS CoV-2 RNA in stool and urine specimens of COVID-19 patients. Front Med (Lausanne). (2023) 10:10. doi: 10.3389/fmed.2023.1148688

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Pool testing, droplet digital PCR, COVID-19 pandemic, new variant outbreak, SARS-CoV-2 wave surveillance, Delta/Omicron tracking

Citation: Pacini A, Paredes F, Heckel S, Ibarra G, Petreli MV, Perez M, Agnella Y, Piskulic L, Allasia MB, Caprile L, Colaneri A and Sesma J (2024) Ready for new waves: optimizing SARS-CoV-2 variants monitoring in pooled samples with droplet digital PCR. Front. Public Health. 11:1340420. doi: 10.3389/fpubh.2023.1340420

Received: 17 November 2023; Accepted: 21 December 2023;
Published: 12 January 2024.

Edited by:

Marina Pekmezovic, CureVac AG, Germany

Reviewed by:

Tesfaye Gelanew, Armauer Hansen Research Institute (AHRI), Ethiopia
Danijela Miljanovic, University of Belgrade, Serbia

Copyright © 2024 Pacini, Paredes, Heckel, Ibarra, Petreli, Perez, Agnella, Piskulic, Allasia, Caprile, Colaneri and Sesma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Juliana Sesma, sesma@rosario-conicet.gov.ar

These authors have contributed equally to this work and share first authorship

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.