SARS-CoV-2 Variants Detection Strategies in Wastewater Samples Collected in the Bangkok Metropolitan Region

Wastewater surveillance is considered a promising approach for COVID-19 surveillance in communities. In this study, we collected wastewater samples between November 2020 and February 2022 from twenty-three sites in the Bangkok Metropolitan Region to detect the presence of SARS-CoV-2 and its variants for comparison to standard clinical sampling. A total of 215 wastewater samples were collected and tested for SARS-CoV-2 RNA by real-time PCR with three targeted genes (N, E, and ORF1ab); 102 samples were positive (42.5%). The SARS-CoV-2 variants were determined by a multiplex PCR MassARRAY assay to distinguish four SARS-CoV-2 variants, including Alpha, Beta, Delta, and Omicron. Multiple variants of Alpha–Delta and Delta–Omicron were detected in the wastewater samples in July 2021 and January 2022, respectively. These wastewater variant results mirrored the country data from clinical specimens deposited in GISAID. Our results demonstrated that wastewater surveillance using multiple signature mutation sites for SARS-CoV-2 variant detection is an appropriate strategy to monitor the presence of SARS-CoV-2 variants in the community at a low cost and with rapid turn-around time. However, it is essential to note that sequencing surveillance of wastewater samples should be considered complementary to whole genome sequencing of clinical samples to detect novel variants.


Wastewater Enrichment and RNA Extraction
The samples were thoroughly mixed by inversion. Then, 40 mL of mixed wastewater from 200-250 mL of the collected samples, maintained at 4 • C, was transferred to a fresh 50 mL centrifuge tube with glass beads (HiMedia, West Chester, PA, USA), then homogenized by a sample disruption instrument, FastPrep-24 5G (MP Biomedicals, Solon, OH, USA), for 60 s, at 6 m/s. The homogenized sample was centrifuged at 3000× g for 10 min at 4 • C to remove the precipitate. The supernatant was carefully transferred to a new tube. A ZR Urine RNA Isolation Kit (Zymo Research, Irvine, CA, USA) was used to isolate the total RNA from the samples, following the manufacturer's protocol. The enrichment process was conducted prior to the RNA extraction. Briefly, 30-35 mL of wastewater supernatant was filtered through a 1.6 µm pore size glass fiber filtration membrane (ZRC GF Filter). Next, 700 µL of Urine RNA Buffer was applied to the filter, and the flow-through was collected in a nuclease-free tube. Following the manufacturer's protocol, the RNA was purified with a Zymo-Spin IC Column and eluted with 50 µL of DNase/RNase-Free Water. The RNA was immediately tested for a SARS-CoV-2 real-time PCR assay and was kept at −80 • C for further analysis.

SARS-CoV-2 Real-Time PCR Detection
The presence of SARS-CoV-2 RNA was determined using the Fosun COVID-19 RT-PCR Detection Kit (Fosun, Shanghai, China) by detecting three different genes: the E gene, ORF1ab gene, and N gene. The real-time RT-PCR was performed following the manufacturer's protocol. Briefly, 10 µL of the extracted RNA was added to a mixture of 14 µL of SARS-CoV-2 reaction reagent and 6 µL of RT-PCR enzyme. The exogenous internal positive control (IPC) RNA was included in the PCR reaction reagent to monitor the appearance of potential PCR inhibitors. Thermal cycling was performed at 50 • C for 15 min for reverse transcription, followed by 95 • C for 3 min, 5 cycles of 95 • C for 5 s to 60 • C for 40 s, and then 40 cycles of 95 • C for 5 s to 60 • C for 40 s using a Bio-Rad CFX 96 PCR instrument (Bio-Rad, Hercules, CA, USA). According to the manufacturer, samples with a PCR Ct value below 40 were considered positive. ogy, Armed Forces Research Institute of Medical Sciences (AFRIMS), Bangkok, Thailand, was used to evaluate the wastewater enrichment extraction protocol. Eight virus concentrations (200 µL each) from 10 7 to 1 copies/mL were spiked into 50 mL of SARS-CoV-2 negative wastewater. The RNA from the SARS-CoV-2-spiked wastewater samples was extracted as indicated above. The RNA from the virus isolate was extracted using the ZR Urine RNA Isolation Kit (Zymo Research, USA) without the above enrichment step. The extraction was conducted in triplicate for each virus concentration. The SARS-CoV-2 RNA was detected by real-time PCR, as indicated above. The PCR Ct value from the SARS-CoV-2-spiked wastewater was compared with the RNA from the direct extraction from the virus isolate in the VTM at the same viral concentration ( Table 2). The limit of detection (LOD) was determined to be the last concentration where the PCR was positive.  --10  ------1 ------* Mean of triplicate assays.

SARS-CoV-2 Variants of Concern (VOC) Detection by Multiplex PCR MassARRAY (PMA)
The PCR-positive RNA samples with a Ct value less than 34 were further characterized for the SARS-CoV-2 variants using the multiplex PCR-nucleotide mass spectrometry technology (MassARRAY) or PMA developed in our laboratory [18]. Briefly, cDNA was generated using 8 µL of RNA extract via reverse transcription using the Superscript III First-Strand Synthesis System (Thermo Fisher Scientific, Waltham, MA, USA), following the manufacturer's protocol. Ten primer pairs and probes were designed to detect eight specific locations on the spike gene to differentiate the Alpha (S982A, DEL69/70), Beta (DEL241/243, A701V), Delta (DEL157/158, L452R), and Omicron (E484A, S477N) variants. In addition, other targets were included as a control, including D614G for VOCs and the N gene as an optimistic gene target.
The multiplex PCR consists of three steps. First, the amplification of 10 targets: 3 µL of cDNA was added into PCR reagents comprising 0.5 µL of PCR buffer (10×), 0.4 µL of MgCl 2 (25 mM), 0.1 µL of dNTPs (25 mM), 0.2 µL of PCR enzyme (5 U/µL), 1 µL of amplification primer mix, and HPLC-grade H 2 O to a total volume of 6 µL. The thermocycling conditions were 95 • C for 2 min, followed by 45 cycles at 95 • C for 30 s, 56 • C for 30 s, and 72 • C for 1 min, with a final incubation at 72 • C for 5 min.
Second, to eliminate excess dNTPs from the previous step, 0.17 µL of SAP buffer (10×), 0.3 µL of SAP enzyme, and 1.53 µL of HPLC-grade H 2 O were added to the first-step PCR products to a total volume of 8 µL and incubated in thermal cycles at 37 • C for 40 min, followed by 85 • C for 5 min.
Third, the single-base extension reaction was performed with an iPLEX Pro Reagent Kit (Agena Bioscience, San Diego, CA, USA) to hybridize and elongate the specific extension primer for each target to distinguish the mutation nucleotide from the wild type. The reaction was performed with terminator nucleotides (ddNTPs), 0.2 µL of iPLEX buffer (10×), 0.2 µL of iPLEX Terminator Mix (10×), 0.04 µL of iPLEX Pro enzyme (33 U/µL), and 0.94 µL of the extension primer mix, mixed with the second-step products, and H 2 O was added to a total volume of 10 µL. The reaction was performed at 95 • C for 30 s, followed by Viruses 2023, 15, 876 5 of 14 95 • C for 5 s, 5 cycles of 52 • C for 5 s to 80 • C for 5 s, and then 40 cycles of 95 • C for 5 s, 52 • C for 5 s to 80 • C for 5 s, with a final extension at 72 • C for 3 min.
Next, 29 µL of HPLC-grade H 2 O and 13 µL of clean resin were added to the extension products for desalination. Afterward, the supernatant was dispensed onto the SpectroCHIP Array through the MassARRAY Nanodispenser and scanned using a MassARRAY Analyzer. The results were analyzed using MassARRAY Typer Software (v.4.1.8.3). The mutations were distinguished with nucleotide matrix-assisted laser desorption/ionization-time of flight mass spectrophotometry (MALDI-TOF MS) based on different masses. The peaks in the mass spectrum were identified as mutations. The reported mass was generated via the MassARRAY Typer Software. The mass of nucleotides detected by the MassARRAY differentiated the wild type from the mutations. The samples with more than one variant in the wastewater specimen, wild type and mutant nucleotides were considered heterozygous.

Whole-Genome Sequencing by Next-Generation Sequencing (NGS)
The RNA from seven SARS-CoV-2-positive wastewater samples with a PCR Ct value lower than 26 were characterized for whole genome sequences using NGS technology. The libraries were performed using ARTIC protocol v3 and v4 primers for the viral RNA amplification [19]. The DNA library preparation was performed using an Illumina DNA Prep kit (Illumina, San Diego, CA, USA). The sequencing was achieved using a MiSeq Reagent Kit v2 (2 × 250 nucleotides) and run on an Illumina MiSeq platform (Illumina, USA). The sequencing reads were assembled by mapping to a reference sequence (SARS-CoV-2 isolate Wuhan-Hu-1) using BWA v.0.7.17 [20].

Data Analysis
We retrieved the sequence number of the SARS-CoV-2 variants of concern (Alpha, Beta, Gamma, Delta, and Omicron) in Thailand from GISAID on 7th May 2022 (Supplement Table S1). This data package contains information about the identified variants, sample locations, and collection dates [21]. This data provided information regarding each VOC's prevalence over time and was used as a reference to compare whether the PMA results in this study corresponded with the VOC emergent period or trends in Thailand determined from the whole genome sequencing data of clinical samples.

Country's Reported Case Number
From November 2020 to February 2022, Thailand reported 2,888,147 cases of COVID-19. During this time, Bangkok reported nearly 600,000 cases, accounting for approximately 20% of the national cases ( Figure 1).

Country's Reported Case Number
From November 2020 to February 2022, Thailand reported 2,888,147 cases of COVID-19. During this time, Bangkok reported nearly 600,000 cases, accounting for approximately 20% of the national cases ( Figure 1).

Evaluation of RNA Extraction Protocol
The efficacy of the RNA extraction from wastewater samples used in this study was assessed by comparing the sensitivity of direct viral extraction and the wastewater enrichment extraction protocols. The detection limit of the direct extraction method was 10 2 copies/mL on all the tested genes (PCR Ct values = 32.13, 31.15, and 29.88 for the E, ORF1ab, and N, respectively), whereas the limit of detection of the wastewater extraction was 10 3 copies/mL on all the tested genes (PCR Ct values = 32.09, 30.59, and 29.68 for the E,

Evaluation of RNA Extraction Protocol
The efficacy of the RNA extraction from wastewater samples used in this study was assessed by comparing the sensitivity of direct viral extraction and the wastewater enrichment extraction protocols. The detection limit of the direct extraction method was 10 2 copies/mL on all the tested genes (PCR Ct values = 32.13, 31.15, and 29.88 for the E, ORF1ab, and N, respectively), whereas the limit of detection of the wastewater extraction was 10 3 copies/mL on all the tested genes (PCR Ct values = 32.09, 30.59, and 29.68 for the E, ORF1ab, and N, respectively). The wastewater extraction protocol showed a 10 times higher detection limit than direct extraction (Table 2).

Limit of Detection of the SARS-CoV-2 RNA Assay from a Wastewater Sample
The limit of detection (LOD) of the SARS-CoV-2 RNA PCR from a wastewater sample was further calculated from the previous experiments. The LOD unit (copies per liter of wastewater) was calculated from the total copies of the spiked virus isolates at each dilution. Six concentrations of the virus isolate from 2 × 10 6 to 20 copies were spiked into 50 mL negative wastewater specimens for final concentrations ranging from 40,000 to 0.4 copies/mL of wastewater, and the limit of detection of SARS-CoV-2 Delta in wastewater was 4 copies/mL (Table 3).

SARS-CoV-2 PCR Results in Wastewater Samples from the Community
A total of 215 samples were collected from 14 sites in Bangkok from November 2020 to February 2022 ( Table 1). The wastewater samples from two hotels (A and B) that were used as quarantine hotels for travelers were SARS-CoV-2 PCR-positive at the initiation of sampling and continued to be positive for 13 and 9 weeks, respectively ( Table 1) Table S2).  Table S2).
Twenty lavatory wastewater samples were collected from the aircraft from abroad that landed at Suvarnabhumi International Airport on 29 November and 21 December 2021. Seven of the twenty samples were positive with PCR Ct values of 29.14-33.33, 30.09-34.4, and 30.0-33.8 for the N, E, and ORF1ab genes, respectively, demonstrating that at least one COVID-19 patient traveled within the aircraft.  (Table 4).

Screening of SARS-CoV-2 Variants in Wastewater
Fifty-eight wastewater samples (n = 58) tested positive by PCR and, with a Ct value less than 34, were further tested for viral variants by a multiplex PCR MassARRAY assay (PMA). Ten PCR targets were simultaneously amplified in the same assay; eight targets were the mutation points for differentiating SARS-CoV-2 variants: Alpha (S982A, DEL69/70), Beta (DEL241/243, A701V), Delta (DEL157/158, L452R), and Omicron (E484A, S477N), and two targets were included as a control, including D614G for VOCs and the N gene as an optimistic gene target. The SARS-CoV-2 Beta variant using DEL241/243 and A701V as markers was not detected from any of the tested wastewater samples (Figure 2). The Alpha variant (S982A and DEL69/70 markers) was first detected in wastewater collected in January 2021, the same period as the first clinical sample deposited in the GISAID. The first reported case of the Delta variant in Thailand was in January 2021, and cases started to increase rapidly in May 2021, as deposited in GISAID (Supplement Table S1). Unfortunately, our study did not collect wastewater samples between March and June 2021. The Delta variant (DEL157/158 and L452R mutants) was first detected in wastewater from our study in July 2021. Mixed detection of Alpha and Delta variants was found from July to August 2021 (Figure 2). In January 2022, our study first detected the Omicron variant (DEL69/70, E484A, and S477N) in the wastewater, as well as mixed with Delta variants (Figure 2). The Omicron variant was first reported in GISAID Thailand in November 2021, and the prevalence sharply increased in January 2022 (Supplement Table S1).   The specificity of the PMA was validated with the whole genome sequence from nextgeneration sequencing of the same wastewater specimens (Figure 3). All seven whole genome sequences were deposited into GISAID with the following accession numbers: EPI_ISL_16616378, EPI_ISL_16616380, EPI_ISL_16637171, EPI_ISL_16616381, EPI_ISL_16616382, EPI_ISL_16616383, and EPI_ISL_16616384. In addition, the target point mutations were identified from the Sequence Alignment Map file using Unipro UGENE [22] to extract the nucleotide coverage per base. The suspected SARS-CoV-2 variants of all seven samples concluded from the eight specific point mutations reported by the PMA and the NGS are the same (Figure 3). The specificity of the PMA was validated with the whole genome sequence from nextgeneration sequencing of the same wastewater specimens (Figure 3). All seven whole genome sequences were deposited into GISAID with the following accession numbers: EPI_ISL_16616378, EPI_ISL_16616380, EPI_ISL_16637171, EPI_ISL_16616381, EPI_ISL_16616382, EPI_ISL_16616383, and EPI_ISL_16616384. In addition, the target point mutations were identified from the Sequence Alignment Map file using Unipro UGENE [22] to extract the nucleotide coverage per base. The suspected SARS-CoV-2 variants of all seven samples concluded from the eight specific point mutations reported by the PMA and the NGS are the same (Figure 3).

Discussion
Wastewater surveillance is a cost-effective tool to survey outbreaks. Wastewater surveillance is used in many countries to track COVID-19 circulation and SARS-CoV-2 variants at the community level [12][13][14]23]. In combination with other indicators, wastewater surveillance is beneficial in targeting the appropriate COVID-19 responses and interven-

Discussion
Wastewater surveillance is a cost-effective tool to survey outbreaks. Wastewater surveillance is used in many countries to track COVID-19 circulation and SARS-CoV-2 variants at the community level [12][13][14]23]. In combination with other indicators, wastewater surveillance is beneficial in targeting the appropriate COVID-19 responses and interventions by providing an early indication (4-7 days) of changes in the incidence and levels of virus circulation [24][25][26].
Our study demonstrates a simple and easy method for detecting SARS-CoV-2 RNA and its variants from wastewater specimens. There have been many surveillance studies to detect viruses from wastewater since 2020 that have used different methods of samplecollection, enrichment, and purification of viral RNA [27]. The studies that require special equipment, for example, using the ultracentrifugation at the enrichment step to increase the viral RNA copy, raise limitations in conducting wastewater surveillance in resourcelimited countries. Before extraction, this study used commercially available extraction kits (Zemo Research, Irvine, CA, USA) to enrich the viral RNA on a filter. Our study's LOD of wastewater SARS-CoV-2 detection protocol is four copies/mL (Table 3). Similar results on the LOD of wastewater PCR detection have been reported by Thongpradit et al.; the LOD for three gene targets, N, ORF1ab, and S, was 1.63, 1.20, and 1.51 copies/mL, respectively. The study demonstrates the ability to detect SARS-CoV-2 RNA in wastewater, even with only one patient in the COVID-19 isolation ward [28].
To demonstrate the practicality of our method, eight different wastewater sampling locations (fourteen sites) were included in this study: a hotel, market, factory, construction camp, hospital, field hospital, condominium, and aircraft. Positive samples were found at all the sites, including the sewage lavatory from the aircraft, where access to clinical samples is extremely limited. The sampling frequency varied between 1 and 16 times at each sampling site. The viral copy was not determined in this study, but we used the PCR Ct value to estimate the amount of the virus. The PCR Ct value varied from 23.31-37.76, 22.19-36.44, and 20.64-37.35 when detected in N, E, and ORF1ab, respectively. Of the 102 samples, 86 were PCR-positive for all three genes, 15 were positive for two genes, and one was positive for only the ORF1ab gene (PCR Ct value = 33.8). This result is comparable to detecting SARS-CoV-2 RNA in clinical samples, in which the assay's sensitivity is increased when multiple gene targets are used for detection [29]. The absence of SARS-CoV-2 RNA in the wastewater samples from the market in the second week demonstrates that interventions such as closing and cleaning help to control and eliminate the virus from the system. The continuing SARS-CoV-2 RNA-positive results from the wastewater collected from the condominium confirm that there were COVID-19 patients in the building during the sampling period, even though there were no reports of COVID-19 cases.
Aircraft wastewater surveillance can form part of a broader strategy for public health surveillance at national borders [30]. Aircraft wastewater samples were collected from 20 flights in our study, and seven samples were positive. The PCR Ct values of the N, E, ORF1ab range were 29.14-33.33, 30.08-34.4, and 30.0-33.8, respectively. These PCR Ct value levels were similar to the wastewater tested from the hospital when four patients were admitted (N = 31.94, ORF1ab = 33.89, and S = 32.01), as reported previously [28]. This finding revealed that COVID-19-infected travelers entered the country despite strict travel regulations from November to December 2021. SARS-CoV-2 RNA was not detected in the wastewater collected from clean areas, such as the non-COVID-19 IPD and OR wards. The wastewater study at the hospital site demonstrated that the COVID-19 prevention measures used for the clean areas, such as the non-COVID-19 IPD and OR wards, are proper for screening and excluding COVID-19-infected patients. The hospital infection control measures included collecting the contact history for the seven days before admission, the temperature measurement, and a SARS-CoV-2 RT-PCR test from each patient before hospital admission. The patient caregiver was also tested for SARS-CoV-2 by RT-PCR on the first day and every subsequent seven days. SARS-CoV-2 detection in wastewater can help monitor the effectiveness of the control measures.
The study also highlights that genomic surveillance of SARS-CoV-2 in wastewater could track clinically relevant mutations, which might be underrepresented in clinical sequencing data. The wastewater samples contained low viral loads of SARS-CoV-2 because they were diluted by the massive volume of water and low viral RNA quality, due to degradation from the treatment process [15,16]. Testing of the SARS-CoV-2 viral variants by the whole genome was therefore limited; the limit of detection for successful whole genome sequencing is a PCR Ct value less than 26 [18]. The PCR Ct values from the wastewater samples in this study and others [28,31] were mostly higher than 30. An alternative approach for SARS-CoV-2 variant detection is applying target amplification and detecting short genome sequences characteristic of genome variants [17,[32][33][34] using digital PCR technology. Our study used multiplex PCR to amplify and then used MassARRAY technology (PMA) to detect the short spike gene fragments in wastewater. This technology allows the detection of wild type and mutant nucleotides for the SARS-CoV-2 variant in the same reaction. It has been successfully used to simultaneously identify SARS-CoV-2 variants from biological specimens on nine mutation sites on a sample with a PCR Ct lower than 27 [35]. Recently, PMA was developed to detect 24 mutation targets simultaneously in one reaction [18]. The sensitivity of each target varied from 10 to 1000 copies/mL.
In this study, primers and probes were designed to detect four significant VOCs from January 2021 to February 2022, including Alpha, Beta, Delta, and Omicron variants, using PMA. These included eight unique mutation points and two control sites (the N gene as the SARS-CoV-2 positive control and D614G as the VOC site). The overall sensitivity (call rate) of 58 wastewater samples was 100% (n = 26), 90% (n = 25), 80% (n = 6), and 70% (n = 1) (Figure 2). A mixture of genomic material (wild type and mutant) of multiple SARS-CoV-2 variants present in wastewater samples can be distinguished by the spectrum of a mass detected from MassARRAY ( Figure 4). MassARRAY technology allows the detection of heterozygous mutation points. It is an assay to detect coinfection and multiple infections from the same sample in one reaction. The process takes 8 h, and the reagent cost is cheaper than the whole genome sequencing assay [18]. PMA is the method of choice for variant identification in wastewater samples, in which the pooled sample can contain multiple virus strains (variants). In contrast, the presence of multiple SARS-CoV-2 variants could not be determined from the consensus whole genome sequence obtained from the NGS technology in this study, as only the major mutation was called for each position. Therefore, all the point mutations found in each position were extracted directly from the variant calling file (Figure 3).
COVID-19 cases infected with the Beta variant were rare in Thailand (112 reported in the GISAID), so it is not unexpected for the wastewater testing to reveal no infection. The association between the detection of the SARS-CoV-2 Alpha variant in wastewater and the occurrence of the variant in clinical specimens in January 2021 from our study was consistent with the results from the USA, Canada, Israel, and Europe [12,13,[33][34][35]. The SARS-CoV-2 variant results from wastewater collected from January 2021 to February 2022 indicated the replacement of consecutive variants in circulation over time (Figure 2). The Delta variant was first detected in wastewater in July 2021, while the sharp rise of the Delta variant in the clinical samples deposited in GISAID was found in May 2021; unfortunately, there were no samples collected and tested during April and June 2021 in our study. The SARS-CoV-2 Omicron variant was first reported in a patient in Thailand in November 2021, but it was first found in wastewater in January 2022. Two wastewater samples from the community were collected in November 2021. The disagreement between the variant's presence in the clinical and wastewater samples from our study was due to the low number of samples (the variant data were received from 58 of 104 positive samples) and inconsistency of the sampling time and period.
The limitations of this study are the inconsistency of the sampling number and the frequency of the sampling time. In addition, the comparable data of the SARS-CoV-2 variant detection in wastewater samples and GISAID data from many laboratories in Thailand reveal an alternative method for monitoring the relative prevalence of variants in each community, especially in a location where randomized sequencing from clinical specimens might be a challenge. Indeed, there are limitations to detecting new pathogens using the PMA method. However, it is essential to note that sequencing surveillance of wastewater samples should be considered complementary to whole genome sequencing of clinical samples. Therefore, it is necessary to monitor the whole genome sequence from patients or wastewater, and then update any new mutation points in the PCR and PMA panels. The results will be further implemented as the basis for COVID-19 outbreak and variant screening in communities that will work adequately in Thailand. rate) of 58 wastewater samples was 100% (n = 26), 90% (n = 25), 80% (n = 6), and 70% (n = 1) (Figure 2). A mixture of genomic material (wild type and mutant) of multiple SARS-CoV-2 variants present in wastewater samples can be distinguished by the spectrum of a mass detected from MassARRAY ( Figure 4). MassARRAY technology allows the detection of heterozygous mutation points. It is an assay to detect coinfection and multiple infections from the same sample in one reaction. The process takes 8 h, and the reagent cost is cheaper than the whole genome sequencing assay [18]. PMA is the method of choice for variant identification in wastewater samples, in which the pooled sample can contain multiple virus strains (variants). In contrast, the presence of multiple SARS-CoV-2 variants could not be determined from the consensus whole genome sequence obtained from the NGS technology in this study, as only the major mutation was called for each position. Therefore, all the point mutations found in each position were extracted directly from the variant calling file (Figure 3). COVID-19 cases infected with the Beta variant were rare in Thailand (112 reported in the GISAID), so it is not unexpected for the wastewater testing to reveal no infection. The association between the detection of the SARS-CoV-2 Alpha variant in wastewater and the occurrence of the variant in clinical specimens in January 2021 from our study was consistent with the results from the USA, Canada, Israel, and Europe [12,13,[33][34][35]. The SARS-CoV-2 variant results from wastewater collected from January 2021 to February 2022 indicated the replacement of consecutive variants in circulation over time (Figure 2). The Delta variant was first detected in wastewater in July 2021, while the sharp rise of the Delta variant in the clinical samples deposited in GISAID was found in May 2021; unfortunately, there were no samples collected and tested during April and June 2021 in our study. The SARS-CoV-2 Omicron variant was first reported in a patient in Thailand in November 2021, but it was first found in wastewater in January 2022. Two wastewater samples from the community were collected in November 2021. The disagreement between the variant's presence in the clinical and wastewater samples from our study was due to the low number of samples (the variant data were received from 58 of 104 positive samples) and inconsistency of the sampling time and period.

Conclusions
Our findings demonstrate a simple method for SARS-CoV-2 detection from wastewater collected from various settings. It can be used to monitor the circulation of the SARS-CoV-2 virus in the community and assess the infectious control measures in a COVID-19-free environment, such as an operating room, non-COVID-19 ward, or aircraft. In addition, it can be expanded to other contagious infectious diseases such as gastrointestinal tract infection or Mpox. The samples for wastewater surveillance should be collected continuously (at least two specimens per site) and consistently. In addition, the primers for multiplex PCR should be updated frequently to avoid missed SARS-CoV-2 variant identification.