- 1Department of Infectious Diseases, Kobe Institute of Health, Kobe City, Japan
- 2Planning Division, Sewage Works Department, Public Construction Projects Bureau, Kobe City, Japan
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of the current coronavirus disease 2019 (COVID-19) pandemic and associated respiratory infections, has been detected in the feces of patients. Therefore, determining SARS-CoV-2 RNA levels in sewage may help to predict the number of infected people within the area. In this study, we quantified SARS-CoV-2 RNA copy number using reverse transcription quantitative real-time PCR with primers and probes targeting the N gene, which allows the detection of both wild-type and variant strain of SARS-CoV-2 in sewage samples from two wastewater treatment plants (WWTPs) in Kobe City, Japan, during the fourth and fifth pandemic waves of COVID-19 between February 2021 and October 2021. The wastewater samples were concentrated via centrifugation, yielding a pelleted solid fraction and a supernatant, which was subjected to polyethylene glycol (PEG) precipitation. The SARS-CoV-2 RNA was significantly and frequently detected in the solid fraction than in the PEG-precipitated fraction. In addition, the copy number in the solid fraction was highly correlated with the number of COVID-19 cases in the WWTP basin (WWTP-A: r = 0.8205, p < 0.001; WWTP-B: r = 0.8482, p < 0.001). The limit of capturing COVID-19 cases per 100,000 people was 0.75 cases in WWTP-A and 1.20 cases in WWTP-B, respectively. Quantitative studies of RNA in sewage can be useful for administrative purposes related to public health, including issuing warnings and implementing preventive measures within sewage basins.
Introduction
The coronavirus disease 2019 (COVID-19) outbreak, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), started in December 2019 and then spread worldwide in 2020 (Huang et al., 2020). In Japan, the COVID-19 pandemic can be divided into five major waves as of December 2021 (Worldometer, 2020). COVID-19 infections occurred mainly due to alpha and delta variants in the fourth (March to June 2021) and fifth (August to October 2021) waves, respectively (Hodcroft, 2021). SARS-CoV-2 mainly causes symptoms of upper respiratory tract infections, but it may also cause severe pneumonia including acute respiratory distress syndrome (ARDS; Wu et al., 2020). Because SARS-CoV-2 can also infect the digestive organs, high levels of the virus can be detected in the feces of infected individuals (Lescure et al., 2020). Stool samples from 48.1% of patients tested positive for SARS-CoV-2 RNA and 70.3% of these patients had stool viral RNA that remained positive even when respiratory specimens were negative, as shown in a meta-analysis study (Cheung et al., 2020). Thus, SARS-CoV-2 is not only a case of respiratory distress, but is also one of the most important causative agents of human gastroenteritis.
Since municipal wastewater contains microorganisms derived from human feces, the concentration of pathogens in sewage is affected by infectious disease epidemics, mainly gastroenteritis, in the watershed population. For example, a previous study showed that the RNA copy number of human gastroenteric norovirus in sewage was significantly related to the number of gastroenteritis cases in the wastewater treatment plant (WWTP) basin (Kazama et al., 2017). In addition, hepatitis E virus was detectable in raw sewage when 1%–4% of residents in a WWTP basin were infected (Miura et al., 2016). Monitoring the prevalence of SARS-CoV-2 in wastewater is also considered to be an effective approach for predicting the COVID-19 epidemic following the detection of SARS-CoV-2 RNA in sewage in several countries (Ahmed et al., 2020a; La Rosa et al., 2020b; Wang et al., 2020; Wurtzer et al., 2020, 2022; Albastaki et al., 2021; Wehrendt et al., 2021; Kevill et al., 2022). In some of these studies, the RNA copy number of SARS-CoV-2 in the sewage was correlated with the number of COVID-19 clinical cases (Medema et al., 2020a; Carrillo-Reyes et al., 2021; Nagarkar et al., 2021; Street et al., 2021; Wurtz et al., 2021; Monteiro et al., 2022). In Japan, SARS-CoV-2 RNA was first detected in secondary treated wastewater before chlorination at a WWTP in Yamanashi Prefecture in April 2020 (Haramoto et al., 2020). Although experiments to detect SARS-CoV-2 RNA have been carried out on wastewater in other regions in Japan for 1–3 months (Hata et al., 2021; Kitamura et al., 2021; Nagashima et al., 2021; Torii et al., 2021), few studies have reported the relationship between the number of COVID-19 cases and the amount of RNA detected with long-term monitoring.
Owing to the low concentration of pathogenic viruses in wastewater, a method for concentrating and detecting these viruses is necessary (Haramoto et al., 2018). Various methods aimed at concentrating RNA in environmental samples to improve detection have been evaluated (Ahmed et al., 2020b; Weidhaas et al., 2021). Among them, electronegative membrane adsorption, polyethylene glycol (PEG) precipitation, and ultrafiltration are frequently used to detect non-enveloped viruses such as poliovirus and norovirus (World Health Organization, 2003; Kazama et al., 2017). Several studies have employed these methods to enrich SARS-CoV-2 RNA (Foladori et al., 2020; La Rosa et al., 2020a; Sangkham, 2021); however, since SARS-CoV-2 is an enveloped virus, the concentration efficiency of the virus differs from that of non-enveloped viruses. Recently, comparative studies of extraction methods showed that SARS-CoV-2 RNA was more abundantly detected in the solid fraction of sewage samples, i.e., the pellet obtained by centrifugation of sewage samples (Kitamura et al., 2021; Westhaus et al., 2021). On the other hand, another study reported that approximately 90% of the SARS-CoV-2 RNA was present in the liquid phase of the influent wastewater (Weidhaas et al., 2021). The assessment of RNA concentration from sewage solids is important, even though solid residues are largely removed in studies on wastewater treatment.
Evaluating the efficiency of the process after wastewater concentration requires the use of a control virus. Pepper mild mottle virus (PMMoV) is the most abundant virus in human feces (Zhang et al., 2006), owing to which it can be easily quantified without spiking in a wastewater sample. High concentrations of PMMoV have been detected in water environment (Rosario et al., 2009; Haramoto et al., 2013; Hughes et al., 2017), and this virus has been used as an internal control for virus detection in wastewater in several studies (D’Aoust et al., 2021; Gerrity et al., 2021; Rosiles-González et al., 2021).
In the present study, we examined the pelleted solid fraction and the product of PEG precipitation of the supernatant fraction of wastewater samples for SARS-CoV-2 RNA. Wastewater was collected once a week from two WWTPs in Kobe, Japan, during the fourth and fifth pandemic waves of COVID-19, and the relationship between the SARS-CoV-2 RNA copy number in the two sample types was correlated with the reported number of COVID-19 cases in the corresponding sewage basin.
Materials and Methods
Sample Collection
Influent wastewater samples were collected once a week from 24 February to 27 October 2021 at WWTP-A (n = 36) and WWTP-B (n = 36) in Kobe City, Japan. The samples were grabbed from the influent, which comprised wastewater before treatment at the WWTPs. All sampling was performed at a fixed time every Wednesday, except on May 6 (WWTP-A and WWTP-B) and August 12 (WWTP-B), in which samples were collected on a Thursday. The samples were collected in sterile plastic bottles and kept frozen at −20°C until analysis. As of December 2021, the city had 1,515,907 inhabitants, of which 98.7% were covered by six WWTPs. WWTP-A and WWTP-B covered 51.5% of the population, received 51.3% of the total wastewater, and treated a total flow of 364,100 m3 per day. The amounts of rainfall (mm/day) and influent flow (m3/day) were measured as routine work at each WWTP.
RNA Extraction
Viral RNA was extracted from each sewage sample after centrifugation, to produce a solid fraction, and after PEG precipitation of the supernatant, to produce a PEG-precipitated fraction, following the procedures of previous studies with minor modifications (Jones and Johns, 2009; Kitamura et al., 2021). Specifically, 160 ml of each sample was divided equally into four aliquots (40 ml each) held in 50 ml tubes and centrifuged at 10,000 × g for 30 min. RNA was extracted from the resulting pellet (solid fraction sample) using the NucleoBond RNA Soil kit (Macherey-Nagel, Düren, Germany) following the manufacturer’s instructions. Meanwhile, the entire supernatant was precipitated using PEG 8000 (final concentration 10%; Promega, Madison, WI, United States) and NaCl (final concentration 1 M; Wako, Tokyo, Japan) by incubating at 4°C overnight with gentle rotation. After centrifugation at 10,000 × g for 60 min, the precipitate was resuspended in 500 μl of phosphate buffer solution (pH 7.0, 0.067 mol/L; Nacalai Tesque, Kyoto, Japan). RNA was extracted from 140 μl of the PEG-precipitated suspension using a QIAamp Viral RNA Kit (Qiagen, Hilden, Germany). RNA was also extracted from 140 μl of raw unconcentrated sewage samples using a QIAamp Viral RNA Kit (Qiagen).
Reverse Transcription-Quantitative PCR
To quantify viral RNA in the samples, reverse transcription-quantitative PCR (RT-qPCR) was performed using the Thermal Cycler Dice Real Time System III (Takara Bio, Shiga, Japan). SARS-CoV-2 RNA was quantified in the solid fraction, the PEG-precipitated sample, and unconcentrated sewage samples using the TaqMan Fast Virus 1-Step Master Mix (Applied Biosystems, Foster City, CA, United States) with combination of CDC 2019-nCoV_N1 and CDC 2019-nCoV_N2 primers and probes, which can be used to detect both the wild type and variant strains. The primer sequences used are described in Supplementary Table S1. Thermal cycling conditions included an initial incubation at 50°C for 5 min and initial denaturation at 95°C for 20 s, followed by 45 cycles of denaturation at 95°C for 3 s and annealing and extension at 60°C for 30 s, as per the manufacturer’s instructions. PMMoV RNA was also quantified in the same samples using the One Step PrimeScript III RT-qPCR Mix (Takara Bio). Thermal cycling conditions for PMMoV included an initial incubation at 52°C for 5 min and initial denaturation at 95°C for 10 s, followed by 45 cycles of denaturation at 95°C for 5 s and annealing and extension at 60°C for 30 s, as per the manufacturer’s instructions. All RT-qPCR analyses included both positive (standard DNA/RNA) and negative (water) controls. The analysis of SARS-CoV-2 samples was performed in duplicate; samples in which only one of the reactions showed a positive amplification were considered as negative overall. To obtain a standard curve for each assay, 10-fold dilution series of a standard plasmid DNA (PMMoV; 5 × 103–5 × 106; Haramoto et al., 2013) or RNA (2.5 × 100 and 5 × 100–5 × 104; SARS-CoV-2 RNA positive control; Takara Bio) solutions were prepared for each assay. RNA copy numbers were calculated from the Ct values using the standard curves. The limit of quantification for SARS-CoV-2 was set at 2.5 copies/reaction (Supplementary Figure S1).
Calculation of Viral RNA Copy Number and Recovery Rate
The copy number of the viral RNA calculated using RT-qPCR was corrected to copy/L as previously described (Qiu et al., 2022), as follows:
where Vextracted RNA is the total volume of the extracted RNA, VRNA in each PCR reaction is the volume of RNA assayed in a RT-PCR reaction, Vwastewater concentrate is the sample volume after concentration, Vwastewater concentrate for RNA extractions is the volume of wastewater concentrate used for RNA extraction, and Vinitial wastewater is the volume of initial wastewater sample processed.
Meanwhile, the recovery rate (%) was calculated using PMMoV quantitative value as follows:
Statistical Analysis
The daily newly reported number of COVID-19 cases was obtained from the Coronavirus Infection Status Report, which is a public database from Kobe City (2021). This database includes the symptomatic cases reported by hospitals and private COVID-19 test centers, as well as asymptomatic cases tested for contact tracing conducted by Public Health Management Center, Kobe City. The number of COVID-19 cases in each investigated basin of the WWTPs were provided by the Public Health Division, Public Health Management Center, Kobe City. Statistical analyses were performed using GraphPad Prism 8 software (GraphPad Prism Software, San Diego, CA, United States). The slope, intercept, and coefficient of determination (R2) values between standard RNA and Ct value, and case numbers and RNA copy numbers were calculated using linear regression; in the latter case, regression through the origin was used. The detection frequency of SARS-CoV-2 between the solid and liquid fractions was assessed using Fisher’s extract test. The Mann–Whitney U test was used to compare RNA copy numbers of the different sewage treatment samples. To compare case numbers and RNA copy numbers, and amount of rainfall/influent flow and RNA copy numbers, the correlation coefficient (r) was calculated using Spearman’s correlation coefficient.
Results
Comparison of the Solid Fraction and PEG-Precipitated Sewage Samples
The PMMoV was used as a control for the RNA extraction process. The PMMoV RNA copy numbers extracted from the non-enriched (raw), solid, and PEG-precipitated samples were in the range 1.2 × 109–1.6 × 1010, 1.5 × 107–2.0 × 108, and 3.1 × 107–5.5 × 108 copies/L, respectively (Figure 1A). The recovery rates calculated from the PMMoV RNA copy numbers of solid and PEG-precipitated samples were in the range 0.5%–3.8% and 1.7%–20% (Figure 1B). The RNA copy number and recovery rate of PMMoV RNA in PEG-precipitated fraction were significantly higher than in solid fraction. The detection frequency of SARS-CoV-2 RNA in solid samples was significantly higher than that in PEG-precipitated samples in both WWTP-A and WWTP-B (Figure 2A). SARS-CoV-2 RNA copy numbers in the solid fraction and PEG-precipitated fraction were in the range 3.1 × 102–3.8 × 104 and 7.6 × 102–2.4 × 104, respectively (Figure 2B). While SARS-CoV-2 RNA copy numbers in the solid fraction samples were significantly higher than in the PEG-precipitated fraction from WWTP-A, no significant difference was observed between both the fraction from WWTP-B.
Figure 1. Quantification of pepper mild mottle virus (PMMoV) RNA. (A) RNA samples from wastewater treatment plant (WWTP)-A (orange circles) and WWTP-B (blue squares) were quantified for PMMoV RNA copy number via RT-qPCR. (B) Recovery rate (%) of concentrated samples calculated from PMMoV RNA copy number. Raw: unconcentrated raw wastewater; solid: solid fraction; PEG: PEG-precipitated fraction. Bars indicate the mean ± SD (n = 36). Statistical significance was calculated using the Mann–Whitney U test. ***p < 0.001.
Figure 2. Detection of SARS-CoV-2 RNA in the solid and PEG-precipitated fractions. (A) Positivity rate of SARS-CoV-2 RNA in the solid and PEG-precipitated fractions. Values of p were calculated using Fisher’s extract test. (B) The SARS-CoV-2 RNA copy number of positive samples in solid and PEG-precipitated fractions from WWTP-A (orange circles) and WWTP-B (blue squares) were quantified using RT-qPCR. Bars indicate the mean ± SD. Statistical significance was calculated using the Mann–Whitney U test. *p < 0.05, NS, not significant.
Association Between RNA Copy Number and Infected Case Number
The SARS-CoV-2 RNA copy number in the solid fraction was highly correlated with the number of COVID-19 cases reported between 24 February 2021 and 27 October 2021 (WWTP-A: r = 0.8205, p < 0.001; WWTP-B: r = 0.8482, p < 0.001; Figure 3). While SARS-CoV-2 RNA in the PEG-precipitated fraction and unconcentrated raw samples was significantly correlated with the number of COVID-19 cases between 24 February 2021 and 27 October 2021 (PEG: WWTP-A: r = 0.6237, p < 0.001; WWTP-B: r = 0.7803, p < 0.001, and Raw: WWTP-A: r = 0.6285, p < 0.001; WWTP-B: r = 0.4517, p = 0.0057; Supplementary Figure S2), the correlation between RNA copy number and COVID-19 cases was lower than that in the solid fraction. The relationships between SARS-CoV-2 RNA and COVID-19 cases in WWTPs basin were evaluated using linear regression analysis (Figure 4). In this study, detection limit of SARS-CoV-2 RNA by RT-PCR was 2.5 copy/reaction, which is calculated to be 625 copy/L in the solid fraction. The limit of capturing COVID-19 cases per 100,000 people calculated using slope were 0.75 cases in WWTP-A and 1.20 cases in WWTP-B, respectively. When the effect of rain on viral RNA concentrations was evaluated, no inverse correlation was found between the amount of rainfall and SARS-CoV-2 RNA in solid fraction (WWTP-A: r = 0.1973, p = 0.2047; WWTP-B: r = 0.1539, p = 0.3701; Supplementary Figures S3A,B). Likewise, no inverse correlation between the amount of influent flow and RNA concentration was observed (WWTP-A: r = 0.4088, p = 0.0133; WWTP-B: r = 0.1652, p = 0.3356; Supplementary Figures S3C,D).
Figure 3. SARS-CoV-2 RNA copy number in the solid fraction and COVID-19 case numbers. The SARS-CoV-2 copy numbers in the solid fraction from (A) WWTP-A (orange circle) and (B) WWTP-B (blue squares) are plotted. The number of new COVID-19 cases per day in WWTP basin is indicated by the gray bars, and the seven-day moving average is indicated by the red line.
Figure 4. Correlation of COVID-19 cases and SARS-CoV-2 copy number. The daily reported COVID-19 cases per 100,000 people and SARS-CoV-2 RNA copy number in WWTP-A (A) and WWTP-B (B) were analyzed using linear regression. Lines indicate the linear regression and dotted lines indicate 95% CI.
Discussion
In this study, SARS-CoV-2 RNA was detected at higher levels in the solid fractions than in the PEG-precipitated fractions, and the RNA copy numbers reflected the infection trend numbers during the fourth and fifth pandemic waves in Kobe, Japan.
The fraction in which higher RNA copy numbers were measured differed between non-enveloped PMMoV and enveloped SARS-CoV-2; PMMoV RNA was higher in the PEG-precipitated liquid fraction, whereas SARS-CoV-2 RNA was higher in the solid fraction. This difference may be related to the structure of the viruses. In one study, up to 26% of enveloped viruses, such as mouse hepatitis virus and bacteriophage φ6, were bound to the solid fraction, whereas only 6% of non-enveloped viruses, such as bacteriophages MS2 and T3, were similarly bound in wastewater samples (Ye et al., 2016). The results of our study are consistent with those of a previous study that reported the detection efficiency of PMMoV in the PEG-precipitated liquid fraction significantly higher than in the solid fraction (Graham et al., 2021). Similar to our findings, RNA extraction from the solid fraction was better than that from liquid fraction for detecting SARS-CoV-2 in previous reports comparing enrichment methods, including pelleted solid fraction, PEG precipitation, electronegative membrane adsorption, and ultrafiltration (Kitamura et al., 2021; Westhaus et al., 2021). In contrast, another study reported that approximately 90% of the SARS-CoV-2 RNA was present in the liquid phase of the influent wastewater compared to the RNA sorbed on the influent solids (Weidhaas et al., 2021). Tomasino et al. (2021) reported that no significant differences were observed in Ct values of SARS-CoV-2 RNA between the liquid and solid phases. These differences are thought to be due to the centrifuge conditions for solid collection. In this study, since centrifugal condition described in previous studies (Kitamura et al., 2021; Westhaus et al., 2021) did not completely precipitate the solid, we set a strong centrifugal condition for ease of work and efficient solid recovery. The different solid removal strategies may result in a high or low representation of the virus in the solid fraction. To minimize errors in RT-PCR detection and quantification, it is recommended that SARS-CoV-2 RNA should be concentrated from both liquid and solid phases of wastewater (Ahmed et al., 2022). Furthermore, in this study, SARS-CoV-2 RNA could be detected using a non-enrichment method as a result of the large number of COVID-19 cases. A previous study reported that the RNA copy number of unenriched wastewater correlated with the number of cases in Marseille, France (Wurtz et al., 2021). These results suggest that the wastewater enrichment methods are not always necessary in areas with high number of COVID-19 cases.
Recovery rate using the PMMoV copy number have been conducted for concentration efficiency of wastewater. Previous studies that compared the liquid and solid fraction showed that the detection efficiency of PMMoV ranged 8.0%–30% in the PEG-precipitated liquid fraction and 6.0%–17% in the solid fraction (Graham et al., 2021), and 12%–102% from liquid fractions and 9.4%–62% from solid fractions (Alamin et al., 2022). Our reported recovery values tended to be lower than the range reported in previous studies. The copy numbers of PMMoV RNA were in the range 8.2 × 106–3.1 × 108 copy/L in the liquid fraction and 1.6 × 102–1.0 × 107 copy/L in the solid fraction (Kitamura et al., 2021). Hasing et al. (2021) reported that the copy numbers of PMMoV were median values of 8.98 × 106 (interquartile range, 6.38 × 106–1.20 × 107) copies per 100 ml in the liquid fraction, and 1.71 × 106 (interquartile range, 1.52 × 106–2.58 × 106) copies per 100 ml in the solids. PMMoV copy numbers in our study did not deviate significantly from the ranges in previous studies, suggesting that the process of wastewater concentration had been performed properly.
In our study, collection of wastewater samples was performed by spot grab sampling, which can affect the SARS-CoV-2 RNA concentration because wastewater flow in a WWTP is increased by the rain. Rainfall was expected to have little effect on RNA concentrations as no inverse correlation between amount of rainfall/influent flow and SARS-CoV-2 RNA in solid was observed. Previous studies showed that composite samples, which were collected by flow-weighting for 24 h, were more detectable than grab samples for SARS-CoV-2 RNA in wastewater (Gerrity et al., 2021; Monteiro et al., 2022). However, grab sampling had been performed for SARS-CoV-2 RNA detection (Randazzo et al., 2020; Carrillo-Reyes et al., 2021; Kitamura et al., 2021; Street et al., 2021; Wehrendt et al., 2021), and correlated with COVID-19 cases (Kitamura et al., 2021; Street et al., 2021). While grab samples have the limitation of low sensitivity, they could be a useful sampling method because they have the advantage of being collected quickly, do not require automated equipment, and were able to reflect COVID-19 cases in our study.
Our data showed that SARS-CoV-2 RNA was detected at higher concentrations as the number of COVID-19 cases increased. This result is consistent with that of a previous study in Tokyo, Japan, which compared SARS-CoV-2 RNA levels in the solid fraction of wastewater with the number of COVID-19 cases from June 2020 to August 2020 (Kitamura et al., 2021). The number of SARS-CoV-2 RNA in primary settled solids collected from primary clarifier was correlated with COVID-19 cases in a study conducted in California, United States (Graham et al., 2021; Wolfe et al., 2021). The number of positive COVID-19 cases has also been correlated with SARS-CoV-2 RNA in wastewater in other countries (Medema et al., 2020a; Carrillo-Reyes et al., 2021; Wurtz et al., 2021); however, the results differed on whether the detection of SARS-CoV-2 RNA increased before or coincident with the number of COVID-19 cases. The viral load in wastewater preceded clinical data by 4 days to 2 weeks in some studies (Medema et al., 2020a; Randazzo et al., 2020; Trottier et al., 2020; Claro et al., 2021; Wu et al., 2022), whereas no time difference was reported in other studies (Peccia et al., 2020; Kitamura et al., 2021). Our data also showed no time difference when correlating the number of COVID-19 cases and SARS-CoV-2 RNA levels. In addition, a study performed daily composite collection of wastewater and reported that although the trend in SARS-CoV-2 RNA levels preceded the number of cases during the first infection wave in France, both measures followed a similar curve in the second infection wave (Wurtz et al., 2021). The difference in the trends of the two waves was probably due to differences in the duration of recognition of the number of cases, that is, in the early stages of a pandemic it is difficult to determine the number of cases as reporting is relatively late, which suggests that the detection of viral RNA precedes case load. In Kobe City, because of active investigation of close contacts and efforts to ascertain the number of infected people, the time difference in reporting the number of cases may be reduced; thus, the viral RNA level and the number of cases correlate without an apparent time difference. In addition, the detection of SARS-CoV-2 RNA from WWTP-B occurred earlier than in WWTP-A. In fact, the number of COVID-19 cases in the WWTP-B basin tended to peak earlier than in WWTP-A.
Linear regression analysis between SARS-CoV-2 RNA copy number and COVID-19 cases showed that the limit of capturing COVID-19 cases per 100,000 people was 0.75 cases in WWTP-A and 1.20 cases in WWTP-B, respectively. SARS-CoV-2 RNA in wastewater was quantifiable in some WWTP basins with daily positive test rates of less than 1 per 10,000 people (Wilder et al., 2021). To detect of SARS-CoV-2, approximately 0.12% and 0.09% of the total population in the WWTP basin area were required to be assessed (Chavarria-Miró et al., 2021). SARS-CoV-2 in wastewater samples collected from five WWTPs in Japan was more likely to be detected when there were more than 10 confirmed cases of COVID-19 per 100,000 people in the basin area, but it was detectable in wastewater even before the number of cases reached 1 per 100,000 people (Hata et al., 2021). The current study found that the capturing COVID-19 cases was equal to or higher than in previous studies. In this study, RNA concentration was detectable in the range 7.6 × 102–2.4 × 104 copy/L when COVID-19 cases per 100,000 people were in the range 0.77–34.4. Previous studies showed that SARS-CoV-2 RNA was detected in the following ranges at the following COVID-19 case rates (per 100,000 people); 1.7 × 103–3.8 × 105 copy/L at 4.8–57.3 cases (D’Aoust et al., 2021), 1.2 × 101–2.2 × 103 copy/L at 0.1–100 cases (Medema et al., 2020a), and 3.0 × 103–2.0 × 104 copy/L at 30–174 cases (Westhaus et al., 2021). Medema et al. (2020b) reported a simulation model of the number of COVID-19 infected people in the population and concentration of SARS-CoV-2 RNA in sewage, and estimated that RNA copy number was approximately 102–105 copy/L at 10–100 cases per 100,000 people. Our results are consistent with these previous studies, suggesting that RNA concentration reflects COVID-19 case numbers.
In Japan, the fourth and fifth waves of SARS-CoV-2 RNA were useful predictions of manifesting COVID-19 cases. The present data indicate that the detection of SARS-CoV-2 RNA in sewage can be used to monitor and predict trends in SARS-CoV-2 infections. This monitoring may provide valuable data even when the number of patients diagnosed with COVID-19 at clinical sites becomes low owing to mass vaccination. In the mass vaccination era, the number of asymptomatic cases is expected to increase, making it more difficult to determine the actual number of cases in the community. The usefulness of wastewater-based epidemiology, which can determine the number of both infected and asymptomatic persons in a community, will increase in the future. These results show the potential of using sewage monitoring, such as RNA levels, in public health, including responding to and the issuing of health warnings within sewerage basins.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding authors.
Author Contributions
YT, TM, and TI: conceptualization. YT and EI: methodology and writing—original draft preparation. NO and TM: resources. YT, EI, SM, AM, RN, and NN: investigation. YT: formal analysis and funding acquisition. RN and TI: supervision. YT, EI, SM, AM, RN, NN, NO, TM, and TI: writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This study was partially supported by a grant from the Daido Life Welfare Foundation for Regional Health and Welfare Research from 2021 to YT.
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.
Acknowledgments
The authors thank Kanna Kodama at the Planning Division, Sewage Works Department, Public Construction Projects Bureau, Kobe City; Satsuki Shirai at the Facility Management Division, Central Sewage Management and Water Environment Control Center, Public Construction Projects Bureau, Kobe City; and Tomohiro Mori at the Facility Management Division, Eastern Sewage Management and Water Environment Control Center, Public Construction Projects Bureau, Kobe City, for their cooperation in wastewater sampling. The authors also thank Kengo Mukai at the Kobe Institute of Health for their useful advice.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.892447/full#supplementary-material
Supplementary Figure S1 | Detection limit of SARS-CoV-2 RNA using RT-qPCR. Each copy of the standard SARS-CoV-2 RNA (Ct value) was confirmed using RT-qPCR measurement. The slope, intercept, and coefficient of determination (R2) values were calculated using linear regression.
Supplementary Figure S2 | SARS-CoV-2 RNA copy number in the PEG-precipitated or unconcentrated sample and COVID-19 case numbers. The SARS-CoV-2 copy numbers in the liquid fraction concentrated using PEG precipitation (A,B) or unconcentrated raw sample (C,D) from WWTP-A (A,C; orange circles) and WWTP-B (B,D; blue squares) are plotted. The number of new COVID-19 cases per day in WWTP basin is indicated by the gray bars, and the seven-day moving average is indicated by the red line.
Supplementary Figure S3 | The amount of rainfall and SARS-CoV-2 RNA copy number in the solid fraction. The SARS-CoV-2 RNA copy number in the solid fraction from WWTP-A (A,C; orange circles) and WWTP-B (B,D; blue squares) are plotted. The amount of rainfall in each WWTPs is indicated by the purple bars (A,B) and the amount of influent flow is indicated by the black line (C,D).
References
Ahmed, W., Angel, N., Edson, J., Bibby, K., Bivins, A., O’Brien, J. W., et al. (2020a). 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. Sci. Total Environ. 728:138764. doi: 10.1016/J.SCITOTENV.2020.138764
Ahmed, W., Bivins, A., Bertsch, P. M., Bibby, K., Choi, P. M., Farkas, K., et al. (2020b). Surveillance of SARS-CoV-2 RNA in wastewater: methods optimisation and quality control are crucial for generating reliable public health information. Curr. Opin. Environ. Sci. Health. doi: 10.1016/j.coesh.2020.09.003 [Epub ahead of print]
Ahmed, W., Simpson, S. L., Bertsch, P. M., Bibby, K., Bivins, A., Blackall, L. L., et al. (2022). Minimizing errors in RT-PCR detection and quantification of SARS-CoV-2 RNA for wastewater surveillance. Sci. Total Environ. 805:149877. doi: 10.1016/j.scitotenv.2021.149877
Alamin, M., Tsuji, S., Hata, A., Hara-Yamamura, H., and Honda, R. (2022). Selection of surrogate viruses for process control in detection of SARS-CoV-2 in wastewater. Sci. Total Environ. 823:153737. doi: 10.1016/j.scitotenv.2022.153737
Albastaki, A., Naji, M., Lootah, R., Almeheiri, R., Almulla, H., Almarri, I., et al. (2021). First confirmed detection of SARS-COV-2 in untreated municipal and aircraft wastewater in Dubai, UAE: The use of wastewater based epidemiology as an early warning tool to monitor the prevalence of COVID-19. Sci. Total Environ. 760:143350. doi: 10.1016/J.SCITOTENV.2020.143350
Carrillo-Reyes, J., Barragán-Trinidad, M., and Buitrón, G. (2021). Surveillance of SARS-CoV-2 in sewage and wastewater treatment plants in Mexico. J. Water Process Eng. 40:101815. doi: 10.1016/j.jwpe.2020.101815,101815
Chavarria-Miró, G., Anfruns-Estrada, E., Martínez-Velázquez, A., Vázquez-Portero, M., Guix, S., Paraira, M., et al. (2021). Time evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in wastewater during the first pandemic wave of COVID-19 in the metropolitan area of Barcelona, Spain. Appl. Environ. Microbiol. 87, e02750–e02720. doi: 10.1128/AEM.02750-20
Cheung, K. S., Hung, I. F. N., Chan, P. P. Y., Lung, K. C., Tso, E., Liu, R., et al. (2020). Gastrointestinal manifestations of SARS-CoV-2 infection and virus load in fecal samples from a Hong Kong cohort: systematic review and meta-analysis. Gastroenterology 159, 81–95. doi: 10.1053/j.gastro.2020.03.065
Claro, I. C. M., Cabral, A. D., Augusto, M. R., Duran, A. F. A., Graciosa, M. C. P., Fonseca, F. L. A., et al. (2021). Long-term monitoring of SARS-COV-2 RNA in wastewater in Brazil: a more responsive and economical approach. Water Res. 203:117534. doi: 10.1016/j.watres.2021.117534
D’Aoust, P. M., Mercier, E., Montpetit, D., Jia, J. J., Alexandrov, I., Neault, N., et al. (2021). Quantitative analysis of SARS-CoV-2 RNA from wastewater solids in communities with low COVID-19 incidence and prevalence. Water Res. 188:116560. doi: 10.1016/j.watres.2020.116560
Foladori, P., Cutrupi, F., Segata, N., Manara, S., Pinto, F., Malpei, F., et al. (2020). SARS-CoV-2 from faeces to wastewater treatment: what do we know? A review. Sci. Total Environ. 743:140444. doi: 10.1016/J.SCITOTENV.2020.140444
Gerrity, D., Papp, K., Stoker, M., Sims, A., and Frehner, W. (2021). Early-pandemic wastewater surveillance of SARS-CoV-2 in southern Nevada: methodology, occurrence, and incidence/prevalence considerations. Water Res X. 10:100086. doi: 10.1016/j.wroa.2020.100086
Graham, K. E., Loeb, S. K., Wolfe, M. K., Catoe, D., Sinnott-Armstrong, N., Kim, S., et al. (2021). SARS-CoV-2 RNA in wastewater settled solids is associated with COVID-19 cases in a large urban sewer shed. Environ. Sci. Technol. 55, 488–498. doi: 10.1021/acs.est.0c06191
Haramoto, E., Kitajima, M., Hata, A., Torrey, J. R., Masago, Y., Sano, D., et al. (2018). A review on recent progress in the detection methods and prevalence of human enteric viruses in water. Water Res. 135, 168–186. doi: 10.1016/j.watres.2018.02.004
Haramoto, E., Kitajima, M., Kishida, N., Konno, Y., Katayama, H., Asami, M., et al. (2013). Occurrence of pepper mild mottle virus in drinking water sources in Japan. Appl. Environ. Microbiol. 79, 7413–7418. doi: 10.1128/AEM.02354-13
Haramoto, E., Malla, B., Thakali, O., and Kitajima, M. (2020). First environmental surveillance for the presence of SARS-CoV-2 RNA in wastewater and river water in Japan. Sci. Total Environ. 737:140405. doi: 10.1016/J.SCITOTENV.2020.140405
Hasing, M., Yu, J., Qiu, Y., Maal-Bared, R., Bhavanam, S., Lee, B., et al. (2021). Comparison of detecting and quantitating SARS-CoV-2 in wastewater using moderate-speed centrifuged solids versus an ultrafiltration method. Water 13:2166. doi: 10.3390/w13162166
Hata, A., Hara-Yamamura, H., Meuchi, Y., Imai, S., and Honda, R. (2021). Detection of SARS-CoV-2 in wastewater in Japan during a COVID-19 outbreak. Sci. Total Environ. 758:143578. doi: 10.1016/J.SCITOTENV.2020.143578
Hodcroft, E. (2021). CoVariants. Available at: https://covariants.org/ (Accessed 17 December 2021).
Huang, C., Wang, Y., Li, X., Ren, L., Zhao, J., Hu, Y., et al. (2020). Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395, 497–506. doi: 10.1016/S0140-6736(20)30183-5
Hughes, B., Beale, D. J., Dennis, P. G., Cook, S., and Ahmed, W. (2017). Cross-comparison of human wastewater-associated molecular markers in relation to fecal indicator bacteria and enteric viruses in recreational beach waters. Appl. Environ. Microbiol. 83, e00028–e000217. doi: 10.1128/AEM.00028-17
Jones, T. H., and Johns, M. W. (2009). Improved detection of F-specific RNA coliphages in fecal material by extraction and polyethylene glycol precipitation. Appl. Environ. Microbiol. 75, 6142–6146. doi: 10.1128/AEM.00436-09
Kazama, S., Miura, T., Masago, Y., Konta, Y., Tohma, K., Manaka, T., et al. (2017). Environmental surveillance of Norovirus genogroups I and II for sensitive detection of epidemic variants. Appl. Environ. Microbiol. 83, e03406–e034016. doi: 10.1128/AEM.03406-16
Kevill, J. L., Pellett, C., Farkas, K., Brown, M. R., Bassano, I., Denise, H., et al. (2022). A comparison of precipitation and filtration-based SARS-CoV-2 recovery methods and the influence of temperature, turbidity, and surfactant load in urban wastewater. Sci. Total Environ. 808:151916. doi: 10.1016/J.SCITOTENV.2021.151916
Kitamura, K., Sadamasu, K., Muramatsu, M., and Yoshida, H. (2021). Efficient detection of SARS-CoV-2 RNA in the solid fraction of wastewater. Sci. Total Environ. 763:144587. doi: 10.1016/j.scitotenv.2020.144587
Kobe City (2021). Coronavirus Infection Status Report (in Japanese). Available at: https://www.city.kobe.lg.jp/a73576/kenko/health/infection/protection/covid_19.html (Accessed December 10, 2021).
La Rosa, G., Bonadonna, L., Lucentini, L., Kenmoe, S., and Suffredini, E. (2020a). Coronavirus in water environments: occurrence, persistence and concentration methods—a scoping review. Water Res. 179:115899. doi: 10.1016/J.WATRES.2020.115899
La Rosa, G., Iaconelli, M., Mancini, P., Bonanno Ferraro, G., Veneri, C., Bonadonna, L., et al. (2020b). First detection of SARS-CoV-2 in untreated wastewaters in Italy. Sci. Total Environ. 736:139652. doi: 10.1016/J.SCITOTENV.2020.139652
Lescure, F. X., Bouadma, L., Nguyen, D., Parisey, M., Wicky, P. H., Behillil, S., et al. (2020). Clinical and virological data of the first cases of COVID-19 in Europe: a case series. Lancet Infect. Dis. 20, 697–706. doi: 10.1016/S1473-3099(20)30200-0
Medema, G., Been, F., Heijnen, L., and Petterson, S. (2020b). Implementation of environmental surveillance for SARS-CoV-2 virus to support public health decisions: opportunities and challenges. Curr. Opin. Environ. Sci Health. 17, 49–71. doi: 10.1016/j.coesh.2020.09.006
Medema, G., Heijnen, L., Elsinga, G., Italiaander, R., and Brouwer, A. (2020a). 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. Environ. Sci. Technol. Lett. 7, 511–516. doi: 10.1021/acs.estlett.0c00357
Miura, T., Lhomme, S., Le Saux, J. C., Le Mehaute, P., Guillois, Y., Couturier, E., et al. (2016). Detection of hepatitis E virus in sewage after an outbreak on a French Island. Food Environ. Virol. 8, 194–199. doi: 10.1007/s12560-016-9241-9
Monteiro, S., Rente, D., Cunha, M. V., Gomes, M. C., Marques, T. A., Lourenço, A. B., et al. (2022). A wastewater-based epidemiology tool for COVID-19 surveillance in Portugal. Sci. Total Environ. 804:150264. doi: 10.1016/J.SCITOTENV.2021.150264
Nagarkar, M., Keely, S. P., Jahne, M., Wheaton, E., Hart, C., Smith, B., et al. (2021). SARS-CoV-2 monitoring at three sewersheds of different scales and complexity demonstrates distinctive relationships between wastewater measurements and COVID-19 case data. Sci. Total Environ. 816:151534. doi: 10.1016/J.SCITOTENV.2021.151534
Nagashima, M., Kawakami, M., Hayashi, M., Kumagai, R., Kasuya, F., Yoshida, I., et al. (2021). RNA Detection by RT-qPCR and non-isolation of SARS-CoV-2 in concentrated wastewater (June–August 2020, Tokyo). Jpn. J. Infect. Dis. doi: 10.7883/yoken.JJID.2021.055 [Epub ahead of print]
Peccia, J., Zulli, A., Brackney, D. E., Grubaugh, N. D., Kaplan, E. H., Casanovas-Massana, A., et al. (2020). Measurement of SARS-CoV-2 RNA in wastewater tracks community infection dynamics. Nat. Biotechnol. 38, 1164–1167. doi: 10.1038/S41587-020-0684-Z
Qiu, Y., Yu, J., Pabbaraju, K., Lee, B. E., Gao, T., Ashbolt, N. J., et al. (2022). Validating and optimizing the method for molecular detection and quantification of SARS-CoV-2 in wastewater. Sci. Total Environ. 812:151434. doi: 10.1016/j.scitotenv.2021.151434
Randazzo, W., Truchado, P., Cuevas-Ferrando, E., Simón, P., Allende, A., and Sánchez, G. (2020). SARS-CoV-2 RNA in wastewater anticipated COVID-19 occurrence in a low prevalence area. Water Res. 181:115942. doi: 10.1016/J.WATRES.2020.115942
Rosario, K., Symonds, E. M., Sinigalliano, C., Stewart, J., and Breitbart, M. (2009). Pepper mild mottle virus as an indicator of fecal pollution. Appl. Environ. Microbiol. 75, 7261–7267. doi: 10.1128/AEM.00410-09
Rosiles-González, G., Carrillo-Jovel, V. H., Alzate-Gaviria, L., Betancourt, W. Q., Gerba, C. P., Moreno-Valenzuela, O. A., et al. (2021). Environmental surveillance of SARS-CoV-2 RNA in wastewater and groundwater in Quintana Roo, Mexico. Food Environ. Virol. 13, 457–469. doi: 10.1007/s12560-021-09492-y
Sangkham, S. (2021). A review on detection of SARS-CoV-2 RNA in wastewater in light of the current knowledge of treatment process for removal of viral fragments. J. Environ. Manag. 299:113563. doi: 10.1016/J.JENVMAN.2021.113563
Street, R., Mathee, A., Mangwana, N., Dias, S., Sharma, J. R., Ramharack, P., et al. (2021). Spatial and temporal trends of SARS-CoV-2 RNA from wastewater treatment plants over 6 weeks in Cape Town, South Africa. Int. J. Environ. Res. Public Health 18:12085. doi: 10.3390/IJERPH182212085
Tomasino, M. P., Semedo, M., Vieira, E., Moreira, P., Ferraz, E., Rocha, A., et al. (2021). SARS-CoV-2 RNA detected in urban wastewater from Porto, Portugal: method optimization and continuous 25-week monitoring. Sci. Total Environ. 792:148467. doi: 10.1016/j.scitotenv.2021.148467
Torii, S., Furumai, H., and Katayama, H. (2021). Applicability of polyethylene glycol precipitation followed by acid guanidinium thiocyanate-phenol-chloroform extraction for the detection of SARS-CoV-2 RNA from municipal wastewater. Sci. Total Environ. 756:143067. doi: 10.1016/J.SCITOTENV.2020.143067
Trottier, J., Darques, R., Ait Mouheb, N., Partiot, E., Bakhache, W., Deffieu, M. S., et al. (2020). Post-lockdown detection of SARS-CoV-2 RNA in the wastewater of Montpellier, France. One Heal 10:100157. doi: 10.1016/J.ONEHLT.2020.100157
Wang, S., Green, H. C., Wilder, M. L., Du, Q., Kmush, B. L., Collins, M. B., et al. (2020). High-throughput wastewater analysis for substance use assessment in Central New York during the COVID-19 pandemic. Environ. Sci. Process Impacts 22, 2147–2161. doi: 10.1039/D0EM00377H
Wehrendt, D. P., Massó, M. G., Gonzales Machuca, A., Vargas, C. V., Barrios, M. E., Campos, J., et al. (2021). A rapid and simple protocol for concentration of SARS-CoV-2 from sewage. J. Virol. Methods 297:114272. doi: 10.1016/J.JVIROMET.2021.114272
Weidhaas, J., Aanderud, Z. T., Roper, D. K., VanDerslice, J., Gaddis, E. B., Ostermiller, J., et al. (2021). Correlation of SARS-CoV-2 RNA in wastewater with COVID-19 disease burden in sewersheds. Sci. Total Environ. 775:145790. doi: 10.1016/j.scitotenv.2021.145790
Westhaus, S., Weber, F. A., Schiwy, S., Linnemann, V., Brinkmann, M., Widera, M., et al. (2021). Detection of SARS-CoV-2 in raw and treated wastewater in Germany – suitability for COVID-19 surveillance and potential transmission risks. Sci. Total Environ. 751:141750. doi: 10.1016/J.SCITOTENV.2020.141750
Wilder, M. L., Middleton, F., Larsen, D. A., Du, O., Fenty, A., Zeng, T., et al. (2021). Co-quantification of crAssphage increases confidence in wastewater-based epidemiology for SARS-CoV-2 in low prevalence areas. Water Res. X. 11:100100. doi: 10.1016/j.wroa.2021.100100
Wolfe, M. K., Topol, A., Knudson, A., Simpson, A., White, B., Vugia, D. J., et al. (2021). High-frequency, high-throughput quantification of SARS-CoV-2 RNA in wastewater settled solids at eight publicly owned treatment works in Northern California shows strong association with COVID-19 incidence. mSystems 6:e0082921. doi: 10.1128/mSystems.00829-21
World Health Organization (2003). Guidelines for environmental surveillance of poliovirus circulation. Available at: https://apps.who.int/iris/handle/10665/67854 (Accessed December 21, 2021).
Worldometer (2020). COVID-19 coronavirus pandemic. Available at: https://www.worldometers.info/coronavirus/japan (Accessed December 20, 2021).
Wu, C., Chen, X., Cai, Y., Xia, J., Zhou, X., Xu, S., et al. (2020). Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease 2019 pneumonia in Wuhan, China. JAMA Intern. Med. 180, 934–943. doi: 10.1001/jamainternmed.2020.0994
Wu, F., Xiao, A., Zhang, J., Moniz, K., Endo, N., Armas, F., et al. (2022). SARS-CoV-2 RNA concentrations in wastewater foreshadow dynamics and clinical presentation of new COVID-19 cases. Sci. Total Environ. 805:150121. doi: 10.1016/J.SCITOTENV.2021.150121
Wurtz, N., Lacoste, A., Jardot, P., Delache, A., Fontaine, X., Verlande, M., et al. (2021). Viral RNA in city wastewater as a key indicator of COVID-19 recrudescence and containment measures effectiveness. Front. Microbiol. 12:664477. doi: 10.3389/FMICB.2021.664477
Wurtzer, S., Marechal, V., Mouchel, J. M., Maday, Y., Teyssou, R., Richard, E., et al. (2020). 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. Euro Surveill. 25:2000776. doi: 10.2807/1560-7917.ES.2020.25.50.2000776
Wurtzer, S., Waldman, P., Levert, M., Cluzel, N., Almayrac, J. L., Charpentier, C., et al. (2022). SARS-CoV-2 genome quantification in wastewaters at regional and city scale allows precise monitoring of the whole outbreaks dynamics and variants spreading in the population. Sci. Total Environ. 810:152213. doi: 10.1016/J.SCITOTENV.2021.152213
Ye, Y., Ellenberg, R. M., Graham, K. E., and Wigginton, K. R. (2016). Survivability, partitioning, and recovery of enveloped viruses in untreated municipal wastewater. Environ. Sci. Technol. 50, 5077–5085. doi: 10.1021/acs.est.6b00876
Keywords: COVID-19, SARS-CoV-2, sewage, wastewater, environmental surveillance
Citation: Tanimoto Y, Ito E, Miyamoto S, Mori A, Nomoto R, Nakanishi N, Oka N, Morimoto T and Iwamoto T (2022) SARS-CoV-2 RNA in Wastewater Was Highly Correlated With the Number of COVID-19 Cases During the Fourth and Fifth Pandemic Wave in Kobe City, Japan. Front. Microbiol. 13:892447. doi: 10.3389/fmicb.2022.892447
Edited by:
Ziyong Sun, Huazhong University of Science and Technology, ChinaReviewed by:
Serena Manara, University of Trento, ItalyMojtaba Pourakbar, Maragheh University of Medical Sciences, Iran
Rashed Noor, Independent University, Bangladesh, Bangladesh
Marta Rusiñol Arantegui, University of Barcelona, Spain
Patrick Monfort, Centre National de la Recherche Scientifique (CNRS), France
Copyright © 2022 Tanimoto, Ito, Miyamoto, Mori, Nomoto, Nakanishi, Oka, Morimoto and Iwamoto. 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: Yoshihiko Tanimoto, eW9zaGloaWtvMDIxOHRhbmltb3RvQGdtYWlsLmNvbQ==; Tomotada Iwamoto, a3gydC1pd210QGFzYWhpLW5ldC5vci5qcA==
†These authors have contributed equally to this work