Summary

Quantification and Whole Genome Characterization of SARS-CoV-2 RNA in Wastewater and Air Samples

Published: June 30, 2023
doi:

Summary

This protocol aims to quantify SARS-CoV-2 RNA in wastewater and air samples to be used for wastewater-based epidemiology studies and to assess the exposure risk to SARS-CoV-2 in indoor and outdoor aerosols. This protocol also describes a tiled amplicon long-template sequencing approach for SARS-CoV-2 whole genome characterization.

Abstract

Wastewater-based epidemiology has emerged as a promising and efficacious surveillance system for SARS-CoV-2 and other infectious diseases in many nations. The process typically involves wastewater concentration, nucleic acid extraction, amplification of selected genomic segments, and detection and quantification of the amplified genomic segment. This methodology can similarly be leveraged to detect and quantify infectious agents, such as SARS-CoV-2, in air samples. Initially, SARS-CoV-2 was presumed to spread primarily through close personal contact with droplets generated by an infected individual while speaking, sneezing, coughing, singing, or breathing. However, a growing number of studies have reported the presence of SARS-CoV-2 RNA in the air of healthcare facilities, establishing airborne transmission as a viable route for the virus. This study presents a composite of established protocols to facilitate environmental detection, quantification, and sequencing of viruses from both wastewater and air samples.

Introduction

In December 2019, a novel disease called COVID-19 emerged, caused by a previously unknown coronavirus, SARS-CoV-21. The resulting global pandemic has presented a significant challenge to clinical and public health laboratories worldwide, as a large number of individuals require testing to accurately assess virus transmission and prevalence in the community. However, in many regions, achieving the necessary level of testing in a timely and spatially comprehensive manner is economically unfeasible2,3. Current surveillance systems based on individual clinical diagnostics rely heavily on symptom severity and individual reporting, as well as the extent to which these symptoms overlap with existing diseases circulating in the population4,5,6,7,8,9,10. Consequently, a high number of asymptomatic cases contributes to a significant underestimation of disease burden7,11.

Due to these challenges, wastewater-based epidemiology (WBE) for COVID-19 surveillance was proposed as a complementary surveillance strategy. WBE was first described in 200112, and was initially used to trace cocaine and other illegal drugs13. This approach relies on the assumption that it is possible to calculate the initial concentration of any substance that is stable in wastewater and excreted by humans8,12. WBE has been successfully implemented in many countries as a complementary and efficient surveillance system for SARS-CoV-23,8,14,15,16. The majority of methods to detect human viruses in aquatic environments follow these steps: concentration, nucleic acid extraction, amplification of the genomic segment (or segments) chosen, and detection/quantification of the amplified genomic segment3.

Another important environment for the detection and quantification of SARS-CoV-2 is in air samples. Initially, SARS-CoV-2 was thought to be transmitted mainly through close personal contact with respiratory droplets from aerosols generated by an infected person while speaking, sneezing, coughing, singing, or breathing17. However, several studies began to report the presence of SARS-CoV-2 RNA in the air, especially in healthcare facilities and other enclosed spaces18,19,20,21. Evidence of SARS-CoV-2 viability in air samples taken indoors in hospitals and other enclosed spaces has been found when the virus concentration was sufficiently high22,23,24. Outdoor studies have generally found no evidence of SARS-CoV-2, except in crowded outdoor spaces21,25,26,27,28,29. As of now, airborne transmission of SARS-CoV-2 has been recognized as a mode of transmission30,31. A recent review study shows the differences between outdoors, where risks of airborne transmission are minimal outside of crowded areas, and indoors, where larger risks could be present in poorly ventilated environments in which strong sources (i.e., number of infected people) could be present. A recent comprehensive review study has highlighted the substantial differences between the risks of airborne transmission in outdoor versus indoor environments, particularly in crowded areas with poor ventilation. The study indicates that the risk of airborne transmission is minimal in outdoor environments, where there is a larger volume of air available for the dilution and dispersion of virus particles32. These findings have important implications for public health policies and guidelines related to COVID-19. By recognizing the significant differences in transmission risks between indoor and outdoor environments, policymakers can develop more effective strategies to mitigate the spread of the virus and protect public health.

There are a variety of methods and protocols for the detection, quantification, and sequencing of SARS-CoV-2 from different environmental samples. This method article aims to present a combination of well-established protocols that allow laboratories with different capacity levels to perform environmental detection, quantification, and sequencing of viruses from wastewater and air samples.

Protocol

All methods described here have been published elsewhere and contain small modifications from the original methods.

1. Wastewater collection and sample pre-processing

NOTE: Due to the low concentrations of SARS-CoV-2 RNA in environmental samples, the implementation of a concentration step is crucial for a successful detection33,34,35. Described here is the first reported method for the detection of SARS-CoV-2 in wastewater36.

  1. Collection and concentration of wastewater sample(s)
    1. Collect 1 L of 24 h composite wastewater sample. Store the sample at 4 °C and proceed with the protocol within 24 h.
    2. Homogenize the sample by gently shaking the bottle. Collect 70 mL into a 100 mL centrifugal tube or divide the sample into two centrifugal tubes of 50 mL.
    3. Spike 20 µL of mengovirus (3.2 x 103 copies/µL) to 70 mL of each sample as an internal control of the concentration process. Alternatively, spike 10 µL of mengovirus (3.2 x 103 copies/µL) into each 50 mL centrifuge tube.
    4. Homogenize the sample and centrifuge at 700 x g for 10 min to remove large particles and organisms (pellet).
    5. Use the resulting supernatant for concentration with centrifugal ultrafiltration devices with a cutoff of 10 KDa by centrifuging at 4,000 x g for 40 min at 4 °C. Elute the concentrate by centrifuging at 700 x g for 40 min and inverting the position of the ultrafiltration device.
    6. Measure the volume of the resulting concentrate. The volume will change depending on the quantity of solids in the sample that clog the centrifugal filters. In this study, the concentrate volumes obtained were between 200-1,200 µL.
    7. Use the resulting concentrate for RNA extraction using a commercial kit or an in-house method. For the samples used in this study, a specific commercial kit for viral RNA extraction with an elution volume of 50 µL is used.
    8. Measure the concentration of the extracted RNA with a fluorometric RNA quantification kit. Divide the extracted RNA into aliquots and freeze at -80 °C until further analysis.
      NOTE: For each RNA isolation procedure set, a negative control of isolation (NCI) should be included, containing only buffers to detect possible contamination during the extraction.
  2. Collection of air samples using a Coriolis compact air sampler
    1. Choose a sampling strategy according to the research goal by defining the sampling frequency and sampling locations.
    2. Set the flow rate on the air sampler (here, 50 L/min for 30 min). Please note that a flow rate of more than 200 L/min can degrade viral RNA, so it is recommended to use a flow rate below 200 L/min.
    3. Place a sterile cone in the air sampler before starting the collection by pulsing the start. Once the sampling is finished, add 5-15 mL of sterile phosphate-buffered saline (PBS) into the cone. Gently vortex or shake the sample by hand for 15 s. Store the samples for up to 24 h at 4 °C or at -80 °C in cryotubes.
    4. An optional concentration step can be performed. Clean and decontaminate the cones after each experiment using bleach.
    5. Extract RNA using a commercial kit or an in-house method. For the samples in this study, a specific commercial kit for viral RNA extraction with an elution volume of 50 µL is used.
    6. Measure the concentration of the extracted RNA with a fluorometric RNA quantification kit. Divide the extracted RNA into aliquots and freeze at -80 °C until further analysis.
      ​NOTE: For each RNA isolation procedure set, a negative control of isolation (NCI) should be included, containing only buffers to detect possible contamination during the extraction.

2. Quantification of SARS-CoV-2 RNA by real-time-quantitative polymerase chain reaction (RT-qPCR)

NOTE: The below protocol is according to the CDC 2019-Novel Coronavirus (2019-nCoV) RT-PCR diagnostic panel37. Divide the primer/probe mix into several aliquots to avoid freezing and thawing cycles.

  1. Prepare the reaction mastermix for each target (mengovirus; SARS-CoV-2 [N1 and N2 genes]) as in Table 1, in a clean hood in the reagent setup room. Mix the primers/probe mix and enzyme by inversion five times. Alternatively, light pulse-vortex the primers/probe mix and enzymefive times.
    NOTE: Keep all reagents cold during preparation and use. In order to minimize freeze-thaw cycles, it is highly recommended to freeze in aliquots.
  2. Spin down (1,000 x g for 15-30 s) the tubes to collect the contents at the bottom and place the tubes in a cold rack or on ice.
  3. Label 1.5 mL microcentrifuge tubes for each target. Add to each microcentrifuge tube the amount of each reagent needed (volume per reaction times the number of reactions including the needed controls). Mix by pipetting up and down and centrifuge for 5 s to collect the contents at the bottom.
  4. Keep the microcentrifuge tubes in a cold rack and dispense 15 µL into strip PCR tubes or a 96-well plate in a cooling rack. Cover the plate and move it to the nucleic acid handling area (keep it in a cooling rack).
  5. Thaw an aliquot of extracted RNA and gently vortex for 5 s. Pipette 5 µL of RNA (in addition to 5 µL of non-template control [NTC], negative control of isolation [NCI], and positive control [PC]) to each reaction well or tube containing the previously prepared mastermix. Change gloves often as needed to avoid cross-contamination.
  6. Cover the entire reaction plate or tubes and gently vortex. Centrifuge (1,000 x g for 15-30 s) the plate or PCR tubes.
  7. Start the 25 µL RT-qPCR with the following cycling conditions: reverse transcriptase at 45 °C for 10 min; polymerase activation at 95 °C for 10 min; 45 cycles of denaturation at 95 °C for 15 s; and annealing/extension at 60 °C for 45s.
  8. Calculate the recovery efficiency of mengovirus controls as reported by Conte et al.38:
    Equation 1 × 100

3. Sequencing variants in wastewater and data analysis

NOTE: The described protocol is a modified protocol created by Quick et al.39,40. It uses two sets of primers for SARS-CoV-2 genome amplification by PCR tiling methodology-ARTIC primers and VarSkip primers. A combination of primers is used to guarantee the best genome coverage and to minimize the possibility of novel mutations causing primers of one type to fail. In general, the protocol is divided into three parts: reverse transcription (RT) and amplicon generation, sequencing library preparation, and sequencing and data analysis.

  1. Reverse transcription and amplicon generation
    1. Ensure the input samples have a known qPCR cycling threshold (Ct) value to correctly dilute them in PCR-grade water. The Ct value is obtained during the quantification with qPCR. The dilutions are as follows: if the Ct value is 12-15, dilute the samples 1:100; if the Ct value is 15-18, dilute the samples 1:10; if the Ct value is 18-35, do not dilute the sample. Samples above a Ct value of 35 have a high chance of not working.
    2. In a PCR plate, pipette 16 µL of each sample to its position on the plate. Add 4 µL of reverse transcription (RT) mastermix (5x). Cover the plate and start the PCR reaction with the following conditions: 25 °C for 2 min, 55 °C for 20 min, and 95°C for 1 min.
    3. Prepare 10 µM dilutions of each primer set (ARTIC pool A and pool B, and VarSkip pool A and pool B). Prepare the mastermix for each set of primers (ARTIC set A and B, and VarSkip set A and B) as in Table 2. Four mastermixes will result from this.
    4. On a new PCR plate, pipette 20 µL of the mastermix to each corresponding well. Add 5 µL of the RT sample to each mix.
      NOTE: Each sample will have four reaction mixtures-ARTIC pool A and B, and VarSkip pool A and B.
    5. Cover the plate and start the PCR reaction as follows: polymerase activation at 98 °C for 30 s; 35 cycles of denaturation at 98 °C for 15 s; and annealing/extension at 65 °C for 4 min. Spin down (1,000 x g for 15-30 s) the plate. Combine the ARTIC reaction pools and separately combine the VarSkip reaction pools. Each sample will have two 50 µL amplicon reactions.
    6. Add 50 µL of bead-based reagent for PCR purification to each well and mix well by pipetting. Incubate at room temperature for 5 min.
      NOTE: Before every use, mix the PCR purification reagent vigorously to resuspend the beads.
    7. Place the plate on a magnetic separation stand and wait for the beads to form a pellet and the liquid to clear (~5 min). Pipette and discard the supernatant. Keep the PCR plate on the magnetic stand.
    8. Maintaining the PCR plate on the magnetic stand, add 200 µL of 80% ethanol to each sample without touching the bead pellet. Remove the ethanol.
    9. Repeat the previous step. Leave the plate on the magnetic stand uncovered for 30 s to allow the ethanol to evaporate.
    10. Remove the plate from the magnetic stand and add 15 µL of PCR-grade water to each sample. Resuspend the pellet by pipetting. Incubate at room temperature for 2 min.
    11. Place the plate back on the magnetic stand and allow the beads to pellet (2 min). Carefully pipette 15 µL of the supernatant from each sample to a new PCR plate.
    12. Measure the concentration of each sample with a fluorometric RNA quantification kit . Take approximately 50 ng of each sample in 12.5 µL to the next step. If required, dilute the sample in PCR-grade water.
  2. Library preparation
    1. Add 1.75 µL of reaction bufferfrom the Illumina DNA library preparation kit and 0.75 µL of enzyme mix to each sample. Cover the plate, vortex briefly, and spin down (1,000 x g for 15-30 s). Incubate at 21°C for 5 min and at 65 °C for 5 min.
    2. In a new plate, pipette 3 µL of PCR-grade water for each sample to a new PCR plate. Add 0.75 µL of the prepared samples, 1.25 µL of barcodes for sequencing library preparation, and 5 µL of T4 DNA ligase mastermix. Mix well by pipetting, briefly spin down (1,000 x g for 15-30 s) in a centrifuge and incubate at 21 °C for 20 min and at 65 °C for 10 min.
      NOTE: If using less than 25 samples, double all the volumes in this step.
    3. Pool all the samples together by adding 10 µL of each sample to the same low-binding tube. Take 480 µL to the next step.
    4. Add 192 µL of bead-based reagent for PCR purification to the pool and mix well by pipetting. Incubate at room temperature for 10 min.
    5. Place the tube on a magnetic stand and wait for the supernatant to clear and a bead pellet to form. Remove the supernatant. Remove the tube from the magnetic stand.
    6. Add 700 µL of the short fragment buffer (SFB) and mix by pipetting. Place back on the magnetic stand and wait for the pellet to form and the liquid to clear (~5 min). Pipette out the supernatant and discard it.
    7. Repeat the previous step. Leave the tube on the magnetic stand. Add 100 µL of 80% ethanol without touching the bead. Pipette out the ethanol and allow the bead to dry for 30 s.
      NOTE: Do not allow the bead to over-dry.
    8. Remove the tube from the magnetic stand and add 35 µL of PCR-grade water. Mix by pipetting and incubate at room temperature for 2 min. Place the tube on the magnetic stand and allow the bead to pellet and the liquid to clear. Pipette 35 µL of the supernatant to a new tube.
    9. Measure the RNA concentration of the pooled library. Pipette the volume needed to reach an amount of 30-50 ng of RNA and fill it with PCR-grade water to reach a final volume of 30 µL.
    10. Prepare the adapter ligation reaction mix: 30 µL of the pooled library, 5 µL of Adapter Mix II, 10 µL of the ligation reaction buffer (5x), and 5 µL of T4 DNA ligase. Mix by pipetting and spin-down (1,000 x g for 15-30 s). Incubate at room temperature for 10 min.
    11. Add 20 µL of the bead-based reagent for PCR purification and mix by pipetting. Incubate for 10 min at room temperature.
    12. Place the tube on the magnetic stand and wait for the bead pellet to form and the supernatant to become colorless. Carefully pipette out and discard the supernatant.
    13. Remove from the magnetic stand and add 125 µL of SFB. Mix by pipetting and place back on the magnetic stand to separate the beads. When the liquid is clear, pipette and discard the supernatant.
    14. Repeat the previous step. Leave the tube open on the magnetic stand for 30 s for some of the leftover liquid to evaporate.
    15. Remove the tube from the magnetic stand and add 15 µL of the elution buffer. Mix well by pipetting and briefly spin-down (1,000 x g for 15-30 s). Incubate at room temperature for 5 min. Place on the magnetic stand for 2 min.
    16. Pipette 15 µL of the supernatant into a new low-binding tube. This is the final library. Measure the concentration with a fluorometric DNA quantification kit.
      NOTE: In the next step, 12 µL is needed. If the concentration is very high and less than 12 µL of the library is needed, fill up to 12 µL with the elution buffer reagent.
  3. Flow cell loading and sequencing
    1. Insert the flow cell (R9.4.1) into the real-time DNA and RNA sequencing device.
    2. Prepare the priming mix by adding 30 µL of flush tether (FLT) reagent to a tube of flush buffer (FB) reagent. Vortex to mix and spin down (1,000 x g for 15-30 s).
    3. Open the priming port cover. Using a 1,000 µL pipette, set to 200 µL and insert the tip into the priming port. Turn the volume setting wheel and increase the volume until the liquid in the tip is seen. Discard the tip.
      NOTE: Only turn the wheel until a few µL of liquid in the tip is seen. Pulling larger amounts could damage the flow cell.
    4. Slowly add 800 µL of the priming mix to the priming port, taking care not to introduce bubbles. Wait for 5 min.
    5. In the meantime, prepare the libraries for loading. Mix 37.5 µL of the sequencing buffer, 25.5 µL of the loading buffer, and 12 µL of the library.
      ​NOTE: Mix the loading buffer before pipetting. The beads in this reagent settle very quickly.
    6. Open the port cover. Pipette 200 µL of the priming mix into the priming port.
      NOTE: While pipetting into the priming port with the cover port open, one will notice small drops coming from the sample port. Pipette slow enough so the sample port can pull the drops in without spilling over.
    7. Before loading the prepared library, mix well by pipetting. Using a 100 µL pipette set to 75 µL, take the library and pipette drop by drop into the sample port. Do not touch the port and wait for each drop to be absorbed into the port before loading the next drop to avoid overflowing of the liquid.
    8. Close the cover port and priming port. Open the software and start the sequencing experiment. Select basecalling On. For a library with 96 multiplexed samples, 10-12 h of sequencing is usually enough to generate sufficient data.
      NOTE: Using the software, one can insert the sample sheet in .csv format. The sample sheet requires the following data: flow_cell_id, position_id, sample_id, experiment_id, flow_cell_product_code, and kit. The flow cell ID is needed (listed on the side of the flow cell) together with the kit used. For the suggested protocol, the LSK109 kit should be selected.
  4. Sequencing data analysis
    1. Insert the compressed fastq reads into the wf-artic workflow41. Run the workflow separately with two different primer schemes (ARTIC/V4.1 and NEB-VarSkip/v1a). Two ". primertrimmed.rg.sorted.bam" files are created for each sample.
    2. Merge the BAM files into a single file with the "samtools merge" command42. Insert the merged BAM file into Freyja tool, leaving default settings, to recover relative lineage abundances43.
    3. To create a FASTA file, insert the merged BAM file into the "samtools mpileup" and "ivar" tools with default settings44,45. For mutation calling, load the FASTA file into the Nextclade tool46,47.

Representative Results

The results summarized in Table 3 show examples of the detection and quantification of SARS-CoV-2 RNA in wastewater and air samples using the method described in this article. Wastewater samples were collected from wastewater treatment plants in Spain and Slovenia and were considered positive if the Ct was less than 40 in at least two of the three replicates, with quantification considered valid if the Ct had a variation of less than 5%. In Spain and Portugal, indoor and outdoor air samples were collected, and the same rules were applied. Duplicates were used for the air samples, not triplicates as for the wastewater samples, as the aim of the study was to detect SARS-CoV-2 RNA and not to quantify it. Furthermore, an RT-qPCR run was only considered valid if all controls used (as described in this protocol) behaved as expected and if no significant deviations were observed in the mengovirus RNA control. For these samples, the mengovirus recovery efficiency varied between 13.7% and 19.7%. The inhibition was assessed by diluting tenfold and 100-fold extracted RNA and by comparing the resulting RT-qPCR results, as well as using a commercial kit. The instructions for each inhibition control kit should be followed as described by the manufacturers.

Wastewater samples had a standard deviation from 1.91% to 13.98% among the triplicates and ranged from 3.05 x 103 to 2.83 x 108 gene copies/L. Air samples ranged from 6.17 x 103 to 5.48 x 109 with a standard deviation between the duplicates between 0.54% and 10.95%. This is in accordance with previous studies6,48,49,50.

As described in this protocol, the samples were sequenced and an example of three sequenced sample results is summarized in Table 4 and Figure 1. In all three samples, the lineage BA.5.x was assigned as the most prevalent by the Freyja tool (95.3%, 95.4%, and 99.8%).

Figure 1
Figure 1: SARS-CoV-2 lineage prevalence in selected wastewater samples. Lineages were assigned using the Freyja tool. The summary of prevalence numbers does not necessarily reach 100%, as some of the data is not of sufficient quality. Please click here to view a larger version of this figure.

Table 1: Reagents and volumes to be added of each to prepare the mastermix for SARS-CoV2 quantification by RT-qPCR. Please click here to download this Table.

Table 2: Reagents and volumes to be added of each to prepare the mastermix to amplify the complete SARS-CoV-2 genomes in the samples. Please click here to download this Table.

Table 3: Summary of RT-qPCR quantification results of SARS-CoV-2 RNA from wastewater and air samples. The N1 gene was used for quantification and the N2 was used for confirmation (positive or negative). The results are expressed as copies/L. Please click here to download this Table.

Table 4: Detected mutations and lineage prevalence. Lineages were assigned using the Freyja tool and according to the mutations of interest found using Nextclade. Genome coverage for each sample is shown. The summary of prevalence numbers does not necessarily reach 100%, as some of the data is not of sufficient quality. Please click here to download this Table.

Discussion

Microbial and viral detection and quantification using (RT-)qPCR methods have garnered widespread acceptance due to their remarkable sensitivity. However, these techniques face numerous challenges when analyzing environmental samples. Wastewater samples contain an abundance of inhibitory substances that can skew measurements and generate misleading results. To tackle these limitations and enhance precision, a complex protocol was conceived, designed, and implemented. This protocol was tailored by combining protocols from the scientific literature and to specifically address the restrictions inherent in qPCR methods by incorporating a series of stringent controls at each stage of the process. By integrating these rigorous quality control measures, this protocol offers a robust solution to the challenges faced by qPCR-based detection and quantification methods in complex environmental samples51,52. The negative control of isolation (NCI) and non-template control (NTC) are essential tools to assess the likelihood of contamination during the RNA extraction and RT-qPCR reaction. Moreover, the integration of mengovirus control RNA provides a crucial means to evaluate inter-sample variability, as well as the efficiency of the concentration step and the inhibitory effects of agents present in complex environmental samples. In every RT-qPCR reaction, positive controls must be included to confirm the effectiveness of the amplification process. It is critical to carefully monitor each step of the protocol, including the storage and processing time of the samples. The degradation of the samples can be influenced by temperature and the complex composition of wastewater, making it essential to store the samples at a temperature of 4 °C and process them within 24 h. By adhering to these stringent quality control measures, researchers and public health workers can confidently and accurately analyze environmental samples to uncover vital insights into microbial and viral populations and better understand the impact of environmental factors on human health48.

The selection of RNA targets for SARS-CoV-2 detection is a critical factor that can impact the sensitivity of RT-qPCR assays. In this study, the authors evaluated the RdRp, E, N1, and N2 targets and found that the RdRp and E assays had lower sensitivity than the N1 and N2 assays, consistent with previous research53. The N1 and N2 assays use double quencher probes, which have been shown to enhance their sensitivity in (RT-)qPCR assays54. Those reasons lead to the decision of using N1 for quantification and N2 as a confirmation of the RNA presence, in order to eliminate doubts when very low concentrations are found. The observed discrepancies among RT-qPCR targets was previously reported36.

Several studies have employed similar protocols and have reported the successful detection and quantification of SARS-CoV-2 RNA in wastewater during the ongoing COVID-19 pandemic6,48,49,55,56. The reported SARS-CoV-2 RNA concentrations are in accordance with the ones reported in this study and protocol48,50,55. Some of these studies have explored the relationship between active cases and concentrations of RT-qPCR target genes, particularly N1 and N26,36,49,50,55,57,58,59. To reduce variability in the obtained gene concentrations, it has been proposed to multiply the experimental concentrations by the flow rate at the time of sampling to obtain the initial viral load36. This approach has been incorporated into the VATar COVID-19 project in Spain, a national monitoring system for COVID-19 in wastewater60.

As the next phase of COVID-19 wastewater-based epidemiology (WBE) studies, researchers have attempted to detect SARS-CoV-2 mutations in wastewater to estimate the circulation of variants of concern within communities. Previous studies have shown a strong correlation between the variants identified in wastewater and those present in the population at the same time. These findings suggest that wastewater-based monitoring could serve as a valuable tool in tracking the spread of SARS-CoV-2 variants61,62. WBE has shown great promise in monitoring SARS-CoV-2 and can be used for the developing future surveillance strategies for pandemics. This approach is particularly useful in the early stages of a pandemic, when there are limited individual tests and many cases may be asymptomatic, and has a complementary strategy to clinical surveillance. Therefore, WBE can be especially useful in places with low economic capacity that cannot afford other more expensive surveillance strategies. The described method can be easily adapted to other viruses that might be relevant to monitor in the future.

The results of air sampling in indoor and outdoor environments revealed that SARS-CoV-2 RNA is present in both environments, although in lower concentrations in outdoor environments38,63,64. These findings suggest that activities and exposures in which mask-wearing and social distancing are difficult to maintain, such as eating, drinking, or attending music festivals, could pose a higher risk for COVID-19 transmission. To mitigate such risks, it is critical to consider appropriate ventilation measures and the use of masks, especially with the emergence of new variants of concern that are more transmissible, such as the Delta (B.1.617.2) and Omicron (B.1.1.529) variants. In light of the lifting of COVID-19 restrictions in many areas, maintaining the use of masks is recommended to keep community transmission of SARS-CoV-2 at a minimum and prevent future outbreaks of the disease5,21,26,65.

Disclosures

The authors have nothing to disclose.

Acknowledgements

This work was performed with financial support from the Regional Government of Castilla y Leon and the FEDER program (projects CLU 2017-09, UIC315 and VA266P20).

Materials

Adapter+A25+A2:D19+A2:D20+A2+A2:D19 Oxford Nanopore EXP-AMII001 Sequencing
AllPrep PowerViral DNA/RNA Kit Qiagen 28000-50 RNA extraction kit
AMPure XP Beckman Coulter A63880 PCR Purification, NGS Clean-up, PCR clean-up
ARTIC SARS-CoV-2 Amplicon Panel IDT 10011442 SARS-CoV-2 genome amplification
Blunt/TA Ligase Master Mix NEB M0367S Library preparation
CENTRICON PLUS­70 10KDA. Fisher Scientific 10296062 Concentration filters
CORIOLIS COMPACT AIR SAMPLER Bertin Technologies 083-DU001 Air sampler
Duran laboratory bottles Merck Z305200-10EA Sampling Bottles
Flow Cell (R9.4.1) Oxford Nanopore FLO-MIN106D Sequencing
General labarotory consumables (tubes, qPCR plates, etc)
Ligation Sequencing Kit Oxford Nanopore SQK-LSK109 Sequencing
LunaScript RT SuperMix Kit NEB E3010  cDNA synthesis
Mengovirus extraction control Kit Biomérieux KMG Concentration control
Nalgene General Long-Term Storage Cryogenic Tubes Thermofisher 5011-0012 Sample storage
Native Barcoding Expansion 1-12 (PCR-free Oxford Nanopore EXP-NBD104 Barcoding
NEBNext Ultra II End Repair/dA-Tailing Module NEB E7595 DNA repair
NEBNext VarSkip Short SARS-CoV-2 Primer Mixes NEB E7658 SARS-CoV-2 genome amplification
NEBNext Quick Ligation Reaction Buffer NEB B6058S Sequencing 
Phosphate buffered saline Merck P4474 Collection buffer
Phosphate-buffered saline (PBS, 1X), sterile-filtered Thermofisher J61196.AP Elution of air samples
Q5 Hot Start High-Fidelity 2X Master Mix NEB M0494S hot start DNA polymerase
Qubit RNA HS Assay Kit Thermofisher Q32852 RNA quantitation
SARS-CoV-2 RUO qPCR Primer & Probe Kit IDT 10006713 Primer-Probe mix and qPCR positive control
TaqPath 1-Step RT-qPCR Master Mix Thermofisher A15299 RT-qPCR kit

References

  1. Naming the Coronavirus Disease (COVID-19) and the Virus that Causes it. World Health Organization Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease (2020)
  2. Lab Workplace Safety. Centers for Disease Control and Prevention Available from: https://www.cdc.gov/coronavirus/2019-ncov/lab/lab-safety-practices.html (2020)
  3. Gonçalves, J., et al. Centralized and decentralized wastewater-based epidemiology to infer COVID-19 transmission – A brief review. One Health. 15, 100405 (2022).
  4. Dawood, F. S., et al. Estimated global mortality associated with the first 12 months of 2009 pandemic influenza A H1N1 virus circulation: a modelling study. The Lancet Infectious Diseases. 12 (9), 687-695 (2012).
  5. Gonçalves, J., Koritnik, T., Paragi, M. Assessment of weather and atmospheric pollution as a co-factor in the spread of SARS-CoV-2. Acta Bio-Medica: Atenei Parmensis. 92 (3), e2021094 (2021).
  6. Gonçalves, J., et al. Detection of SARS-CoV-2 RNA in hospital wastewater from a low COVID-19 disease prevalence area. The Science of The Total Environment. 755, 143226 (2021).
  7. Mizumoto, K., Kagaya, K., Zarebski, A., Chowell, G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance. 25 (10), 2000180 (2020).
  8. Polo, D., et al. Making waves: Wastewater-based epidemiology for COVID-19 – approaches and challenges for surveillance and prediction. Water Research. 186, 116404 (2020).
  9. Shmueli, G., Burkom, H. Statistical challenges facing early outbreak detection in biosurveillance. Technometrics. 52 (1), 39-51 (2010).
  10. Simonsen, L., et al. Global mortality estimates for the 2009 influenza pandemic from the GLaMOR project: A modeling study. PLoS Medicine. 10 (11), e1001558 (2013).
  11. Oran, D. P., Topol, E. J. Prevalence of asymptomatic SARS-CoV-2 infection: a narrative reivew. Annals of Internal Medicine. 173, 362-367 (2020).
  12. Daughton, C., Jones-Lepp, T. Pharmaceuticals and Personal Care Products in the Environment: Scientific and Regulatory Issues. ACS Symposium Series. , (2001).
  13. Zuccato, E., et al. Cocaine in surface waters: a new evidence-based tool to monitor community drug abuse. Environmental Health. 4, 14 (2005).
  14. Aguiar-Oliveira, M. d. e. L., et al. Wastewater-based epidemiology (WBE) and viral detection in polluted surface water: A valuable tool for COVID-19 surveillance-a brief review. International Journal of Environmental Research and Public Health. 17, 9251 (2020).
  15. García-Encina, P. A. Wastewater-based epidemiology (WBE). Water and Environment Journal. 35 (4), 1162-1163 (2021).
  16. Mao, K., Zhang, H., Pan, Y., Yang, Z. Biosensors for wastewater-based epidemiology for monitoring public health. Water Research. 191, 116787 (2021).
  17. Shereen, M. A., Khan, S., Kazmi, A., Bashir, N., Siddique, R. COVID-19 infection: Origin, transmission, and characteristics of human coronaviruses. Journal of Advanced Research. 24, 91-98 (2020).
  18. Chia, P. Y., et al. Detection of air and surface contamination by SARS-CoV-2 in hospital rooms of infected patients. Nature Communications. 11 (1), 2800 (2020).
  19. Lei, H., et al. SARS-CoV-2 environmental contamination associated with persistently infected COVID-19 patients. Influenza and Other Respiratory Viruses. 14 (6), 688-699 (2020).
  20. Razzini, K., et al. SARS-CoV-2 RNA detection in the air and on surfaces in the COVID-19 ward of a hospital in Milan, Italy. The Science of The Total Environment. 742, 140540 (2020).
  21. da Silva, P. G., Gonçalves, J., Nascimento, M. S. J., Sousa, S. I. V., Mesquita, J. R. Detection of SARS-CoV-2 in the indoor and outdoor areas of urban public transport systems of three major cities of Portugal in 2021. International Journal of Environmental Research and Public Health. 19 (10), 5955 (2022).
  22. Barbieri, P., et al. Molecular detection of SARS-CoV-2 from indoor air samples in environmental monitoring needs adequate temporal coverage and infectivity assessment. Environmental Research. 198, 111200 (2021).
  23. Lednicky, J., et al. Earliest detection to date of SARS-CoV-2 in Florida: Identification together with influenza virus on the main entry door of a university building, February 2020. PLoS One. 16 (1), 0245352 (2021).
  24. Santarpia, J. L., et al. Aerosol and surface contamination of SARS-CoV-2 observed in quarantine and isolation care. Scientific Reports. 10 (1), 12732 (2020).
  25. Chirizzi, D., et al. SARS-CoV-2 concentrations and virus-laden aerosol size distributions in outdoor air in north and south of Italy. Environment International. 146, 106255 (2021).
  26. Hadei, M., et al. Presence of SARS-CoV-2 in the air of public places and transportation. Atmospheric Pollution Research. 12 (3), 302-306 (2021).
  27. Moreno, T., et al. Tracing surface and airborne SARS-CoV-2 RNA inside public buses and subway trains. Environment International. 147, 106326 (2021).
  28. Mouchtouri, V. A., et al. Environmental contamination of SARS-CoV-2 on surfaces, air-conditioner and ventilation systems. International Journal of Hygiene and Environmental Health. 230, 113599 (2020).
  29. Setti, L., et al. Airborne transmission route of COVID-19: why 2 meters/6 feet of inter-personal distance could not be enough. International Journal of Environmental Research and Public Health. 17, 2932 (2020).
  30. SARS-CoV-2 Transmission. Centers for Disease Control and Prevention (CDC) Available from: https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/sars-cov-2-transmission.html (2021)
  31. Coronavirus Disease (COVID-19): How is it Transmitted. World Health Organization Available from: https://www.who.int/news-room/q-a-detail/coronavirus-disease-covid-19-how-is-it-transmitted (2021)
  32. Dinoi, A., et al. A review on measurements of SARS-CoV-2 genetic material in air in outdoor and indoor environments: Implication for airborne transmission. The Science of the Total Environment. 809, 151137 (2022).
  33. Bosch, A., et al. Analytical methods for virus detection in water and food. Food Analytical Methods. 4, 4-12 (2011).
  34. Gonçalves, J., et al. Surveillance of human enteric viruses in coastal waters using concentration with methacrylate monolithic supports prior to detection by RT-qPCR. Marine Pollution Bulletin. 128, 307-317 (2018).
  35. La Rosa, G., Muscillo, M. Molecular detection of viruses in water and sewage. Viruses in Food and Water. , 97-125 (2013).
  36. Medema, G., Heijnen, L., Elsinga, G., Italiaander, R., Brouwer, A. Presence of SARS-Coronavirus-2 RNA in sewage and correlation with reported COVID-19 prevalence in the early stage of the epidemic in the Netherlands. Environmental Science & Technology Letters. 7 (7), 511-516 (2020).
  37. CDC – 2019-nCoV Real-Time RT-PCR Diagnostic Panel Fact Sheet for Healthcare Providers. Centers for Disease Control and Prevention Available from: https://stacks.cdc.gov/view/cdc/85028 (2020)
  38. Conte, M. Airborne concentrations of SARS-CoV-2 in indoor community environments in Italy. Environmental Science and Pollution Research International. 29 (10), 13905-13916 (2022).
  39. nCoV-2019 sequencing protocol v3 (LoCost). protocols.io Available from: https://www.protocols.io/view/ncov-2019-sequencing-protocol-v3-locost-bh42j8ye (2020)
  40. Tyson, J. R. . Improvements to the ARTIC multiplex PCR method for SARS-CoV-2 genome sequencing using nanopore. , (2020).
  41. . ARTIC SARS-CoV-2 Workflow Available from: https://github.com/epi2me-labs/wf-artic (2022)
  42. Li, H., et al. The sequence alignment/map format and SAMtools. Bioinformatics. 25 (16), 2078-2079 (2009).
  43. . Freyja Available from: https://github.com/andersen-lab/Freyja (2022)
  44. Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 27 (21), 2987-2993 (2011).
  45. Grubaugh, N. D., et al. An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar. Genome Biology. 20 (1), 8 (2019).
  46. Hadfield, J., et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. 34 (23), 4121-4123 (2018).
  47. Aksamentov, I., Roemer, C., Hodcroft, E. B., Neher, R. A. Nextclade: clade assignment, mutation calling and quality control for viral genomes. Journal of Open Source Software. 6 (67), 3773 (2021).
  48. Markt, R., et al. Detection and stability of SARS-CoV-2 fragments in wastewater: impact of storage temperature. Pathogens. 10 (9), 1215 (2021).
  49. Kocamemi, B. A., et al. First Data-Set on SARS-CoV-2 Detection for Istanbul Wastewaters in Turkey. MedRxiv. , (2020).
  50. Randazzo, W., et al. SARS-CoV-2 RNA in wastewater anticipated COVID-19 occurrence in a low prevalence area. Water Research. 181, 115942 (2020).
  51. Hoorfar, J., et al. Practical considerations in design of internal amplification controls for diagnostic PCR assays. Journal of Clinical Microbiology. 42 (5), 1863-1868 (2004).
  52. Parshionikar, S. U., Cashdollar, J., Shay Fout, G. Development of homologous viral internal controls for use in RT-PCR assays of waterborne enteric viruses. Journal of Virological Methods. 121, 39-48 (2004).
  53. Nalla, A. K. Comparative performance of SARS-CoV-2 detection assays using seven different primer-probe sets and one assay kit. Journal of Clinical Microbiology. 58 (6), 00557 (2020).
  54. Hirotsu, Y., Mochizuki, H., Omata, M. Double-quencher probes improve detection sensitivity toward Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) in a reverse-transcription polymerase chain reaction (RT-PCR) assay. Journal of Virological Methods. 284, 113926 (2020).
  55. Ahmed, W. First confirmed detection of SARS-CoV-2 in untreated wastewater in Australia: A proof of concept for the wastewater surveillance of COVID-19 in the community. The Science of The Total Environment. 728, 138764 (2020).
  56. Bar-Or, I., et al. Detection of SARS-CoV-2 variants by genomic analysis of wastewater samples in Israel. The Science of the Total Environment. 789, 148002 (2021).
  57. La Rosa, G., Bonadonna, L., Lucentini, L., Kenmoe, S., Suffredini, E. Coronavirus in water environments: Occurrence, persistence and concentration methods – A scoping review. Water Research. 179, 115899 (2020).
  58. Wu, F., et al. SARS-CoV-2 titers in wastewater are higher than expected from clinically confirmed cases. mSystems. 5, 00614 (2020).
  59. Wurtzer, S., et al. Evaluation of lockdown effect on SARS-CoV-2 dynamics through viral genome quantification in waste water, Greater Paris, France, 5 March to 23 April 2020. European Communicable Disease Bulletin. 25 (50), 2000776 (2020).
  60. . VATar COVID-19 | Caso de Exito – Ministerio para la Transición Ecologica y el Reto Demografico Available from: https://esri.es/es-es/descubre-los-gis/casos-de-exito/administracion-/vatar-covod19-miteco-cs (2022)
  61. Nemudryi, A., et al. Temporal detection and phylogenetic assessment of SARS-CoV-2 in municipal wastewater. Cell Reports. Medicine. 1 (6), 100098 (2020).
  62. Rios, G., et al. Monitoring SARS-CoV-2 variants alterations in Nice neighborhoods by wastewater nanopore sequencing. The Lancet Regional Health. Europe. 10, 100202 (2021).
  63. Gomes da Silva, P. Environmental dissemination of SARS-CoV-2 in a University Hospital during the COVID-19 5th wave Delta variant peak in Castile-León, Spain. International Journal of Environmental Research and Public Health. 20, 1574 (2023).
  64. Gonçalves, J., et al. . Exposure assessment of SARS-CoV-2 and Nov GII/GII in aerosols generated by a municipal wastewater treatment plant. , (2022).
  65. Lednicky, J. A., et al. Isolation of SARS-CoV-2 from the air in a car driven by a COVID patient with mild illness. International Journal of Infectious Diseases. 108, 212-216 (2021).

Play Video

Cite This Article
Gonçalves, J., Gomes da Silva, P., Koritnik, T., Bosilj, M., Torres-Franco, A., Diaz, I., Rodriguéz, E., Marcos, E., Mesquita, J. R., García-Encina, P. Quantification and Whole Genome Characterization of SARS-CoV-2 RNA in Wastewater and Air Samples. J. Vis. Exp. (196), e65053, doi:10.3791/65053 (2023).

View Video