Abstract
This is the first study focused on the presence of SARS-CoV-2 in different freshwater environments of an urban setting. Groundwater and surface water reservoirs for drinking water as well as water from receiving rivers of Monterrey Metropolitan Area were sampled repeatedly during a peak phase between October 2020 and January 2021, and the virus RNA was measured by quantitative reverse transcription polymerase chain reaction. Forty-four percent of groundwater samples demonstrated detectable virus loads between 2.9 and 5.6 copies/ml. A significant correlation with sucralose concentration in groundwater reaffirmed the hypothesis of leaching and infiltration of effluents from surface and/or failing sewage pipes, and emphases the importance of water disinfection. Thirteen percent of surface water dam samples tested positive, with values varying between 3.3 and 3.8 copies/ml. Finally, 21% of river samples marked positive for viral RNA, with concentrations ranging from 0.3 to 7.0 copies/ml. Raw wastewater samples taken in the same period show viral loads of up to 3535 copies/ml, demonstrating a dilution effect and/or efficiency of wastewater facilities. Variations of viral loads over time and at sub-metropolitan level in groundwater and surface water in general reflects the reported variation of infections for Monterrey.
1. Introduction
Since the outbreak of the COVID-19 disease, various routes of transmission of SARS-CoV-2 have been verified, and others hypothesized. The current knowledge is that the main transmission between people occurs though respiratory droplets (diameter >5-10 µm) produced by infected individuals when coughing or sneezing (Chan, Yuan et al. 2020; Li et al. 2020; WHO, 2020b). Regarding the presence and persistence of SARS-CoV-2 in wastewater, there is sufficient evidence that indicates that wastewaters may contain both RNA fragments and viable particles of SARS-CoV-2 (Langone et al., 2020; Ihsanullah et al., 2021). A study in stool samples suggests that the SARS-CoV-2 may replicate for 11 days in the gastrointestinal tract of patients even after samples from the respiratory tract are becoming negative (Wu et al. 2020a). According to another experiment, SARS-CoV-2 remained viable for 2 to 6 h in adult’s feces and up to 2 days in children’s feces (Liu 2020b) Several studies have reported the new coronavirus in raw wastewater in Netherlands (Medema et al., 2020), Italy (La Rosa et al., 2020a; Rimoldi et al., 2020); USA (Wu et al., 2020b); France (Wurtzer et al., 2020); Australia (Ahmed et al., 2020a); Spain (Randazzo et al., 2020). Wastewater treatment plants (WWTPs) with tertiary disinfection have been found negative for SARS-CoV-2 (Rimoldi et al., 2020), while effluents from secondary treatments has been found positive (Randazzo et al., 2020; Rimoldi et al., 2020). The presence of SARS-CoV-2 RNA in sewage sludge was reported in a 10-week monitoring study in New Haven, Connecticut, USA (Peccia et al., 2020). These reports show the utility of viral RNA monitoring in municipal wastewater for SARS-CoV-2 infection surveillance at a population-wide level.
Although several authors have hypothesized potential routes in water environments, so far there exists little evidence of the presence of SARS-CoV-2 virus in freshwaters (La Rosa et al. 2020b; Ihsanullah et al. (2021). Water safety starts with the protection of water resources in the catchment, therefore it is mandatory to prevent surface and groundwaters from coming into contact with fecal material. It is expected that groundwater benefits from the pathogens removal due to soil filtration, adsorption on sediment grains and progressive inactivation, and viruses in surface waters are exposed to several potentially inactivating stressors, including sunlight, oxidants, and predation by microorganisms (Langone et al., 2021). Rimoldi et al. (2020) detected viral RNA in three receiving rivers in the Milan area indicating the partial efficiency of the sewage system of the metropolitan area. Haramoto et al. (2020) collected three river samples between March and May 2020 in Japan and reported that none tested positive for SARS-CoV-2 RNA. Guerrero-Latorre et al. (2021) reported virus loads during a peak of the outbreak from three different sites of a river receiving untreated sewage from Quito, Ecuador. To our knowledge, so far there is no evidence of the presence of the virus in surface water reservoirs and aquifers reported.
Considering the rigorous biosafety requirements necessary for working with infectious coronavirus, surrogate markers are useful for assessing method performance and for quality assurance and control during monitoring campaigns (Casanova et al., 2009; Ye et al., 2016; Ahmed et al., 2020b). Similarly, markers of human input in sewer systems could be used as an alternative to trace the presence of SARS-CoV-2 virus. Sucralose is one of the most popular artificial sweeteners and serves as tracer of human wastewater (Kokotou et al., 2012; Voss et al., 2019), and its concentration is correlated with people connected to the sewage system. It is stable at a broad pH range, is heat stable, non-volatile, highly polar and chiral. It is also strongly recalcitrant, only degrading under strongly oxidizing conditions, and is not metabolized by animals or microbes (Soh et al. 2011).
In the present study we evaluated the presence of genetic material from SARS-CoV-2 RNA in different freshwater environments of the Monterrey Metropolitan Area (5.3 million inhabitants) in Northern Mexico. The aim of the study was to perform a survey of viral dispersion and potential implications in the aquatic environment during a peak phase of the epidemic. To that end, we collected groundwater, river water and water from drinking water reservoirs repeatedly between October 2020 and January 2021, and measured SARS-CoV-2 RNA by quantitative reverse transcription polymerase chain reaction (qRT–PCR). Regarding groundwater, we also measured the concentration of artificial sweetener sucralose.
2. Material and Methods
2.1 Study area
The Monterrey Metropolitan Area (MMA) is the second important city in Mexico in terms of population and economy (INEGI, 2021). It comprises 12 municipalities with a total population of approximately 5.3 million inhabitants. The climate is semi-arid with a mean annual temperature and rainfall of 22.3 °C and 622 mm, respectively, with a dry (November-April) and rain season (May-October). Maximum rainfall values occur at mountain flanks while minimum occur at the valley and lower elevations. The urban area is bordered at west and south by mountain ranges varying in composition from clastic marine to carbonate sedimentary rocks reaching elevations up to 2100 m above sea level (masl) (Fig. 1). These mountain ranges, the Sierra Madre Oriental, represents a Mesozoic sedimentary belt that expose sequences of limestones, lutites, with minor sandstones of Late Jurassic to late Cretaceous periods (Santiago Carrasco et al., 2008).
(a) Location of study area, (b) regional view showing surface water reservoirs with sampling points, (c) urban area with main features and groundwater/river sampling points.
The MMA sits in a valley at 580 masl on Quaternary alluvial deposits eroded from the surrounding mountain ranges. The valley is mostly composed of fluvial and alluvial sedimentary deposits as terraces that occurred during accumulation-erosion cycles on the early Quaternary (Martinez and Werner, 1997). Most of these recent deposits can be found as riverbeds in La Silla, Santa Catarina and Pesquería rivers, which flow eastward through the city towards the Gulf Coastal Plain.
Water for MMA is supplied from surface water (58%) and groundwater (42%) reservoirs (SADM, 2021). Surface water is extracted from El Cuchillo dam (4.69 m3/s), Cerro Prieto dam (2.83 m3/s), and La Boca dam (0.45 m3/s). Water from El Cuchillo dam and Cerro Prieto dam is conveyed 108 km and 133 km to MMA, respectively, while La Boca dam connects to the Cerro Prieto aqueduct (Fig. 1b). Raw water from all three dams is purified in the San Roque and La Boca treatment plants before distribution in the city through two water supply pipelines with over 70 km length each (Fig. 1c).
Groundwater is extracted from several aquifer units and wellfields and disinfected locally before introducing into the supply network (Torres-Martinez et al. 2020) (Fig 1a): Buenos Aires (BA) well field (2.11 m3/s) located in a side valley close to the city consists of La Huasteca horizontal filtrating gallery and 23 deep wells with water table depths between 20 and 120 m below ground, extracting water from Early Cretaceous limestone formations; the Santiago (SA) Groundwater System (1.27 m3/s) consists of La Estanzuela spring and three horizontal filtrating gallery’s: San Francisco, Cola de Caballo I and II; the Monterrey Metropolitan Zone (ZM) aquifer (1.08 m3/s) includes wells throughout the metropolitan area, providing water from an unconfined aquifer which consists of altered lutites, conglomerates, gravel, sand and clay, with an average depth to groundwater of 20 m; finally the Mina well field (1.20 m3/s) located about 35 km northwest from MMA.
Used water is treated to over 90% by public wastewater facilities that includes primary and secondary stages in the treatment process. The most important wastewater treatment plants (WWTPs) are Dulces Nombres (7.5 m3/s), Norte (4.0 m3/s), Noreste (1.9 m3/s) and Cadereyta (0.25 m3/s) (Fig 1c). All the mentioned WWTPs discharge the treated water directly or indirectly to Pesquería river, with the exception of Cadereyta WWTP that discharges to the Santa Catarina river. Both rivers are tributaries of San Juan river which in turn flows into the Rio Bravo/Rio Grande. Pesquería and Santa Catarina rivers had a discharge decreasing from 5.5 to 4.4 m3/s at Pesquería hydrometric station, and from 5.6 to 2.5 m3/s at Cadereyta hydrometric station between October 2020 and December 2020, respectively (SMN, 2021).
2.2 Field methods
Groundwater and surface water grab samples were collected in different sites and occasions between October 2020 and January 2021. For groundwater, 42 sites corresponding to production wells of supplying aquifer units of Monterrey (BA well field, ST system and ZM aquifer) were sampled initially between October 29 and November 3, 2020. Of these, 37 wells were public drinking water supply wells and five wells were for industrial purposes. A subset of wells (n=10) was re-sampled two more times in cycles of approximately one month to observe changes over time.
Similarly, samples were taken from three sites of three surface water reservoirs supplying Monterrey (El Cuchillo, Cerro Prieto and La Boca) in October 22-23, 2020, and repeated two more times. Finally, a total of 12 river water grab samples were taken along the three urban rivers Pesquería, Santa Catarina and La Silla in December 10-11,2020, and repeated in January 5-6, 2021. The river sites were selected strategically upstream and downstream of WWTPs discharges to the rivers. For reference, 24-h composite samples of influent of a Dulces Nombres WWTP were taken weekly during the same period.
All samples were collected in sterile 125 ml-HDPE bottles and stored at 4 °C and analyzed within 48 hours. It is known that SARS-CoV-2 is highly stable at 4 °C (Chin et al.; 2020). Groundwater included samples for analysis of sucralose, using 125 ml-HDPE bottles.
2.3 Laboratory methods
2.3.1 RNA and DNA Extraction – QIAamp® Viral RNA Mini
We followed standard procedures to extract and purify nucleic acids from the water samples. Briefly, after viral thermal inactivation (95 °C; 5 minutes), a volume of 500 µl of the water sample were centrifuged for 10 minutes at 1500 G. Then, a volume of 140 µl of the supernatant was added to a mix containing 0.56 µL of Buffer AVL solution (Qiagen, USA) and 5.6 µL of carrier RNA-AVE solution (Qiagen, USA) in a 1.5 mL microcentrifuge tube. This mix was vortexed for 15 seconds and incubated at room temperature (15–25°C) for 10 min and briefly centrifuged to remove drops from the interior surface of the lid. A volume of 560 µl ethanol (96–100%) was added to the sample, and mixed by pulse-vortexing for 15 s.
After mixing, the tube was briefly centrifuged to remove drops from interior surface of the lid. Then, this solution (∼630 µL) was filtered through a QIAamp Mini column (Qiagen, USA) to retain the nucleic acids originally present in the sample. The retained was repeatedly washed with different buffer solutions to elute contaminants and purify the nucleic acids. To this aim, the solution was loaded into the column contained in a 2 mL collection tube, the cap of the tube was closed, and the tube with the column was centrifuged at 6000 x g (8000 rpm) for 1 min.
After centrifugation, the QIAamp Mini column was placed into a clean 2 ml collection tube, and the filtrate was discarded. In a first raising step, 500 µL of 96% ethanol were loaded into the column contained in the 2 mL collection tube, the cap of the tube was closed, and the tube with the column was centrifuged at 6000 x g (8000 rpm) for 1 min. Following these two centrifugation stages, 500 µl of Buffer AW1 (Qiagen, USA) were added into the QIAamp Mini column, the cap of the container tube was closed, and the tube with the column was centrifuged at 6000 x g (8000 rpm) for 1 min. As before, the QIAamp Mini column was placed into a clean 2 ml collection tube, and the filtrate was discarded. In a fourth centrifugation cycle, the QIAamp Mini column was added with 500 µl Buffer AW2 (Qiagen, USA), the cap of the container tube was closed, and the tube with the column was centrifuged at high speed (20,000 x g; 14,000 rpm) for 3 min.
Afterwards, the QIAamp Mini column was placed in a clean 1.5 ml microcentrifuge tube and the filtrate was discarded. In a fifth centrifugation cycle, the QIAamp Mini column was added with 60 µl Buffer AVE (Qiagen, USA) equilibrated to room temperature, the cap of the container tube was closed, and the tube with the column was centrifuged at high speed (6,000 x g; 8,000 rpm) for 1 min.
For DNA extraction, 500 µl of water sample were centrifuged for 5 minutes at 5000 G; 400 µl of the centrifuge supernatant were discarded. The remaining 100 µl were added with 20 µl of proteinase K solution and 80 µl of Buffer ATL (Qiagen, USA), vortexed, and incubated at 56°C for at least 1 hour. The rest of the extraction protocol is analogous to that previously described.
2.3.2 RNA and DNA amplification
We amplified RNA segments of SARS-CoV-2 using two sets of primers (commonly referred as N1 and N2) in each amplification reaction. Both of these primers directed to sequences that encode the N protein of SARS-CoV-2. These primer sets have been recommended and extensively used for the diagnostic of COVID-19 in human samples (Gonzalez-Gonzalez et al., 2020; Nalia et al., 2020) and wastewaters (Medema et al., 2020; Wu et al., 2020; Nemudryi et al., 2020; Randazzo et al., 2020; Haramoto et al., 2020; Scherchan et al., 2020; Peccia et al., 2020).
Similarly, we used two sets of primers to amplify the LAC and LAM regions of the genome of Escherichia coli in the same reaction. E. coli is used as biological indicator of the presence of fecal content in water (Bej et al. 1990; Mo et al., 2002; Reza et al. 2014). The sequences of both the forward and reverse primers used are shown in Table S1.
The quantitative amplification was conducted in a quantitative PCR thermal cycle (Rotor gene Q 5plex, Qiagen, Germany). For the amplification of RNA of SARS-CoV-2 sequences, the amplification mix (final volume of 20 µL) consisted in 10 µL of 2X QuantiNova Syber Green RT-Master Mix, 0.2 µL of QN SYBR Green RT-Mix, 1µL of 10x primer mix (0.5 µM final concentration), and 8.8 µL of RNA extract. For the amplification of DNA sequences of E. coli, the amplification mix (final volume of 20 µL) consisted in 10 µL of 2X QuantiNova Syber Green RT-Master Mix, 1µL of 10x primer mix (0.5 µM final concentration), and 9.0 µL of DNA extract. The amplification cycle consisted in 10 minutes of reverse transcription at 50 °C and 2 minutes of amplification activation at 95 C, followed by 40 iterative cycles of denaturation for 5 seconds at 95 °C and combined annealing and extension for 10 s at 60 °C.
A calibration curve was built to establish the conversion between CT values and equivalent gene copies per milliliter (GC/ml). For this purpose, we used commercial synthetic genetic material that contained the complete N gene from SARS-CoV-2 (Integrated DNA Technologies, IA, USA). Samples containing different concentrations of synthetic nucleic acids of SARS-CoV-2 (in the range of 10 to 100,000 copies mL-1) were prepared by successive dilutions from stocks. This plasmid has been used before as a positive control in amplification assays of SARS-CoV-2 genetic material (González-Gonzalez et al., 2021).
2.3.3 Sucralose quantification
Sucralose is used as artificial sweetener and useful tracer for evidencing the presence of human wastewater in groundwater (Kokotou et al., 2012; Voss et al., 2019). Sucralose was determined using high performance liquid chromatography and mass spectrometric detection (HPLC-MS/MS) after solid-phase extraction (SOE). Isotope-labelled internal standards and an external calibration in tap water were used for quantification. Details of the analytical method are given in the Table S2. The analysis was performed at DVGW-Technologiezentrum Wasser, Karlsruhe, Germany.
2.4 Monitoring of COVID-19 Cases in Monterrey Metropolitan Area
To get an indication of the sensitivity of the monitoring of urban water cycle, a proxy for the period prevalence of COVID-19 in the MMA was created using the reported number of COVID-19 cases per day and the normalized cumulative number of reported COVID-19 cases per day for the year 2020. The normalization was performed by dividing the cumulative number of reported cases by the population size.
3. Results
3.1 Reported cases
The number of reported COVID-19 cases in each of the 12 municipalities and MMA shows that the pandemics evolved at different rates in each of the municipalities as the epidemic spread during 2020 (Fig. 2a). The first infection was reported on March 10, and the number of cases remained relatively low until mid-May, when another increase occurred, and from June 10, the infection maintains a constant increase in the MMA, with exception of November, when the number of cases dropped. Santiago and Monterrey municipalities reported the most cases, followed by Santa Catarina, Guadalupe and San Nicolas.
Reported cases for MMA and its 12 municipalities: (a) Reported daily cases of infection; and (b) normalized cumulative cases. Note: Date obtained from CONACyT (2021). The vertical blue shades indicate the sampling periods.
However, it is worth to mention that these numbers are not directly comparable to other countries or regions because the collection method is not necessarily standardized, and probably the sampling efforts are unequal and asynchrony respect to the real infection dates (Sims and Kasprzyk-Hordern, 2020). As can be observed, the freshwater sampling for this report was performed during the second peak of the outbreak of the epidemic: end of October, end of November 2020, mid-December 2020 and beginning January 2021 (Fig. 2ab).
3.2 Groundwater
Two field campaigns were performed for groundwater. Regarding the first campaign, the qRT–PCR concentration threshold (Ct) average values for SARS-CoV-2 ranged from 30.2 to over 40 (Table 1). Interestingly, near half of them (19 of 40) were positive, and 38% of the positives were below the value of 33. In this study, a sample was arbitrarily defined “positive” when a Ct value was detected in at least two of three replicates. Two of seven cases in BA well field were detected positive with Ct values of 30.2 and 32.4. Galeria 4 is a well at the entrance to the Huasteca highway, with a high urban development in the area prior to the entrance, while Pozo 39 is in the lower area of the Sierra Madre close to ranches and houses. Five out of eight cases in SA system were reported positive with a Ct value between 32.5 and 36.3. Estanzuela, constructed in 1910 as one of the first groundwater sources for Monterrey water supply system, is in a woodland-rural environment, while Cola de Caballo Tunnels and San Francisco Tunnel represent horizontal galleries in piedmont shrubland. Finally, Margarita is a well located in an environment of urban development. Thirteen out of 26 production wells in ZM aquifer indicate positive Ct’s, with values between 30.3 and 34.2. These sites are dispersed in the urbanized MMA. A trend shows a higher proportion of sites affected in the downstream area in the northeastern portion (Apodaca) and no positive sites in the southeastern portion (Contry) of ZM aquifer.
Summary of the results of determination of SARS-CoV-2 RNA and sucralose for groundwater of Monterrey. Note: The Ct value represents the average of triplicate analysis for each sample (Table S3 and S4), ‘n.d.’ indicates not detected and ‘-’ indicates not measured.
Sucralose was detected in 22 out of 40 samples (55%) (Table 1), and its concentrations varied between 0.07 and 2.9 µg/l. In the BA well field, which represents dessert and piedmont shrubland with low population density, none of the samples showed detectable levels. In the SA system, one site (Andares) showed concentrations of sucralose close to the detection limit and one site (Margaritas) one of the highest concentrations. These sites represent residential areas. In the ZM aquifer, 20 out of 25 well sites (80%) had detectable concentrations of sucralose, whose values ranged between 0.1 and 2.7 µg/l. These results are in general consistent with the land use distribution, all except one site with detected sucralose lie in urbanized or industrial plots. In addition, we found a significant correlation between sucralose and Ct values (r2=0.62, n=0.043), but no correlation between Ct values and depth to groundwater.
Those cases which were positive in the first sampling campaign and not close to each other were repeated for a second campaign (Table 1). In the second sampling campaign only 3 out of 10 sites tested positive for SARS-CoV-2. This suggests that groundwater was less affected in the second date, and only three sites are consistently positive in both dates, namely California 2, Lincoln 2 and Puentes 1 in the Monterrey municipality. It is notable that the depth-to-water table of these sites was less than 22 m.
3.3 Surface water
Two sampling campaigns were performed in surface water reservoirs between the end of October and mid-December 2020 (Table 2). For the first period in October 2020 none of the samples were detected positive. For the second sampling period two sites were found positive, one in La Boca dam (33.8) and another in Cerro Prieto dam (33.6). It was not possible to analyze the correlation between Ct values of SARS-CoV-2 and E. coli because only two pairs have quantitative data.
Results of determination of SARS-CoV-2 RNA and E. coli for surface water reservoirs. Note: The Ct value represents the average of triplicate analysis for each sample (Table S5), ‘n.d.’ indicates not detected, and ‘-’ indicates not measured.
With respect to rivers waters, two sampling campaigns were performed in December 2020 and in January 2021 respectively. In December, three out of twelve samples were tested positive, with Ct values ranging from 32.7 to 34.2. The sites with positive values are Pesqueria river downstream of WWTP Norte, Santa Catarina river upstream of WWTP Cadereyta, and La Silla river upstream of Tolteca Park. For the second sampling period, two out of twelve samples were positive, namely Pesquería river upstream WWTP Norte and La Silla river at upstream of Tolteca Park (Table 3). The result for La Silla river is notable because this river receives no treated wastewaters of domestic origin. The Ct values of SARS-CoV-2 correlated with those of E. coli (r2=0.75, n=0.088), however the correlation is weak due the low number of pairs.
Results of the determination of SARS-CoV-2 RNA for rivers in MMA. Note: The Ct value represents the average of triplicate analysis for each sample (Table S6) and ‘n.d.’ means not detected.
3.4 Wastewater
For reference, raw wastewater from influent of Dulces Nombres WWTP was measured for SARS-CoV-2. Between October 25, 2020, and December 13, 2020, 3 out of 8 samples (38%) were positive. The Ct value of positive samples ranged from 23.5 to 31.2 (Table S7).
4 Discussion
4.1 Contextualization of findings in freshwater environments
This is the first study that quantifies the presence of SARS-CoV-2 in different freshwater environments of an urban setting. Previous studies aimed to detect the virus in freshwater focused on receiving rivers (Table 4). For example, Rimoldi et al. (2020) collected grab samples in three sites of receptor rivers of the Milan area on April 14 and 22, 2020. In a first sampling round, all three samples were positive, while in a second round only one out of three was positive. A quantitative analysis was not performed. Similarly, Haramoto et al. (2020) collected grab water samples in a river in Yamanashi Prefecture, Japan, on three different occasions between April 22 and May 7, 2020; they reported that none tested positive for SARS-CoV-2 RNA.
Selected studies on municipal wastewater/sludge and receiving river waters. Note: ‘NA’ means not applied and ‘WW’ indicates wastewater.
Guerrero-Latorre et al. (2020) reported virus loads during a peak of the outbreak (June 5, 2020) from three different sites of a river receiving untreated sewage from Quito city. The authors used RT-qPCR for these determinations and used two different primer sets, namely N1 and N2. All samples were found positive, and the values ranged from 284 to 3190 GC/ml and from 207 to 2230 GC/ml in assays using the N1 and N2 target region, respectively. These values could be related clearly to COVID-19 cases reported in the contributing areas.
4.2 Explanation of viral loads in receiving waters
In the present study, 21% of all river water samples (5 out of 24) were positive regarding viral RNA, and the positive tested samples varied between 2.5 to 7.0 GC/ml, and between 0.3 to 0.8 GC/ml for the first (10-11 December 2020) and second campaign (5-6 January 2021), respectively. Importantly, during this period no significant rainfall was recorded in the Monterrey area that could have had an impact on virus concentration in waters (Tables S8 and S9). These loads are two to three orders of magnitude lower than those reported by Guerrero-Latorre et al. (2020) for Quito’s river. This could be because Monterrey treats more than 90% of its municipal wastewater, while the urban rivers of Quito are impacted by the direct discharge of sewage water from the city (3 million inhabitants). Similarly, the negative results derived from the analysis of river water samples from Yamanashi Prefecture, Japan (Haramoto et al. 2020) and Milan, Italy (Rimoldi et al. 2020) could be possibly attributed to the fact that both studies collected water from rivers receiving treated wastewater.
It is expected that wastewaters from WWTPs with a complete treatment were tested as negative. Thus, the occurrence of SARS-CoV-2 RNA in a few samples in La Silla and Pesquería river water could stem from different sources coexisting in the same basin. For example, aliquots of non-treated sewage can be present because of illicit discharges, malfunction of sewerage systems, and their increased relative contribution during dry periods (Mosley et al. 2015). The lack of separation of the urban runoff waters from the domestic effluents, which causes combined sewer overflows (CSOs), could also be a reason (Rimoldi et al. 2020). CSOs occur usually during high rainfall events. However, the accumulated rainfall between December 2020 and January 2021 in Monterrey was only in the order of 3 mm.
Another reason of high aliquots of untreated sewage in river water could be the organization of local football derbies whose high loads in short time periods may be overburden the capacity of WWTPs releasing untreated wastewater to Pesquería river (SADM, 2020). The case of La Silla river is notable because it receives no relevant treated municipal wastewater due to the sanitary drainage to the other two rivers, therefore illicit discharges or sewerage system malfunction is a plausible explanation for the presence of SARS-CoV-2 genetic material in this water course.
Regarding dam water, only 12% of the samples (2 out of 16) were tested positive for SARS-CoV-2 RNA, with no positive result in the first campaign (22-23 October 2020). The positive samples (which contain 3.3 and 3.8 viral copies/ml) occurred during the second campaign (14-15 December 2020) and only in one site of La Boca and in one site of Cerro Prieto dam, respectively. In both cases urbanization is observed nearby, which suggests that the presence of the virus might be due to failure in the local sewerage system. The observed values are comparable to the range of the urban rivers in Monterrey. The lack of virus loads in the first campaign and the presence in two of the nine sites may reflect the increasing trend of reported cases of infection in the corresponding municipalities during the same period (Fig 1a).
4.3 Viral loads in groundwater reaffirming human sewage impact
The share of groundwater samples containing detectable SARS-CoV-2 RNA was surprisingly high. Twenty-two out of 50 samples (44%) showed viral loads between 2.9 and 5.6 GC/ml. This finding suggests that a fraction of untreated sewage has entered the groundwater system. The origin of the untreated sewerage may be from the surface or from a leaky sewerage system. Torres-Hernández et al. (2020) used isotopic and chemical evidence to determine that nitrate pollution in groundwater from Monterrey was mainly derived from sewage leaks in the urban area. It is evident that organic and viral loads could have been entered to the groundwater system using the same pathway. The significant correlation between SARS-CoV-2 concentrations and sucralose at 0.05 level is another remarkable confirmation of the contribution of raw wastewater to the groundwater and reaffirms possible leaching and infiltrations of effluents from health care facilities, sewage, solid landfill, and drainage water as well as failing sewage pipes in the MMA.
The observed viral load may not be higher than detected in rivers and dams of Monterrey, but it represents certain concern regarding the vulnerability of local groundwater systems to the coronavirus. From the three aquifer units used for water supply, the SA system (63%) is most affected, followed by the ZM aquifer (54%), and the BA well field (22%). Nevertheless, the viral loads observed in wells of the first sampling campaign (29 October 2020 - 4 November 2020) could only partly be reproduced one month later (26-30 November 2020), indicating a decrease in the viral load. This demonstrates how dynamic is the groundwater system in relation to the presence of the coronavirus; this decrease in viral load in groundwater appears to follow the decreasing trend of reported cases of infection during the month of November 2020 (Fig 2a).
From the sampled municipalities of MMA during the first campaign, Apodaca leads with 63% positive samples, followed by Monterrey (50%), and San Nicolas (50%). Coincidently, these are the most affected municipalities considering the officially reported daily cases of infection in Fig. 2a. Guadalupe is also among the most affected municipalities; however, it only is represented by one sampled well. Santiago, the southernmost municipality is the exception as it shows a relatively lower number of cases of infection, but a high incidence of positive cases (63%). This could indicate a different dynamic. Indistinctly, the high level of positive samples in municipalities with highest COVID infections suggests that groundwater samples approximately mirror the infection situation at municipality level.
4.4 Environmental and public health implications
This study is the first evidence that SARS-CoV-2 may find its ways to groundwater through possible leaching and infiltrations of effluents from health care facilities, sewage, solid landfill and drainage water, as well as leakages from sewage pipes, as hypothesized by Ihsanullah et al. (2021). This finding derives in two different lectures. First, groundwater in MMA is currently disinfected by gas chlorination removing pathogenic viruses and bacteria. Since coronaviruses are sensitive to oxidants like chlorine (La Rosa et al., 2020b), it is important to continue strengthening and advancing treatment processes of groundwater, especially in wells located in shallow aquifers in combination with places where sewage effluents from health care facilities, sewage, solid landfill and drainage water are not treated or treated inefficiently (Guerrero-Latorre et al. 2020) and expected to infiltrate, or where sewage pipes could be leaky (Torres-Hernandez et al., 2020). Second, the correlation of the Covid-19 cases with SARS-CoV-2 RT-qPCR results in groundwater (and surface water) has quite potential for verifying the reported number of Covid-19 cases with the real situation.
The concentration of SARS-CoV-2 RNA in untreated wastewater from selected studies worldwide was in the range of not detected to 5600 GC/ml (Table 4). In our study, the monitoring of the influent of Dulces Nombres WWTP showed that between October 25, 2020, and December 13, 2020, 3 out of 8 samples (38%) were positive to SARS-CoV-2, and that the maximum load was 3535 GC/ml (Table S7). This number is quite comparable to other studies of raw WW during outbreaks (Nemudryi et al., 2020; Randazzo et al. 2020; Torttier et al, 2020; Wurtzer et al., 2020; Medema et al., 2020). This shows that the concentration of SARS-CoV-2 in surface (<5.6 GC/ml) and groundwater (<7 GC/ml) in MMA is approximately three orders of magnitude lower than in raw wastewater. This means that the viral load could not be eradicated completed as observed in Haramoto et al. (2020), however the result is similar to Rimoldi et al. (2020).
Importantly, the presence of genetic material of SARS-CoV-2 in freshwater should not be considered as direct indicator of the presence of infective viral particles. Establishing infectivity requires the successful isolation and culture of actual viral particles from water samples (Giacobbo et al., 2021). However, the presence of detectable amounts of genetic material from SARS-CoV-2 in freshwater should not be taken lightly.
4.5 Future directions
This study shows the importance of monitoring programs to determine the fate of SARS-CoV-2 in the urban water cycle. Until now, there is no evidence and enough data to confirm if the water or wastewater containing SARS-CoV-2 could be the potential source of its transmission. Some studies predicted a low risk of SARS-CoV-2 transmission via wastewater (Chin et al., 2020; Rimoldi et al., 2020; Ihsanullah et al., 2021), but the investigations in this topic are still in an early phase. Future research should be oriented towards the understanding of the fate of the virus in different water environments, factors which are favorable and unfavorable for its development, and the potential exposure risk of water contaminated with this virus in its different variants. Experimental research must go hand in hand with the development of conceptual and computational analysis.
A look on studies performed so far shows that there is a lack of standardized protocols for the sampling, detection and quantification of SARS-CoV-2 in water and wastewater (Table 4). For example, in some studies grab samples were taken, while in others 24-hour composite samples were collected. There are significant differences not only in the sample collection, but also in sample storage and treatment, and use (or not) of genetic or chemical traces (i.e., chemical agents indicating human activity or viral tracers used for normalization purposes) among others. This may lead to discrepancies in the results. Currently, the RT-qPCR has been employed widely for detection of SARS-CoV-2 in water samples, however it is imperative to develop standard sampling procedure for accurately extraction, isolation, detection and quantification of the virus.
4 Conclusions
This study evaluated the presence of SARS-CoV-2 RNA in the urban water cycle of a metropolitan area. The aim of the study was to perform a survey of viral dispersion and potential implications in the aquatic environment during a peak phase of the epidemic between October 2020 and January 2021.
Twenty-one percent of all river water samples were positive regarding viral RNA, and the positive tested samples varied between 2.5 to 7.0 GC/ml, and between 0.3 to 0.8 GC/ml for the first and second campaign. These low viral loads demonstrate the dilution effect and/or efficiency of wastewater treatment facilities in MMA. Raw wastewater measured during the same period shows viral loads of up to 3535 GC/ml. Some samples that tested positive in the river stem possibly from aliquots of non-treated sewage due to illicit discharges, malfunction of sewerage system and their relative contribution during the dry period. A further indication of water aliquots of fecal origin of water is a weak correlation with E. coli bacteria. These findings highlight the importance of a good coverage of wastewater treatment facilities such as in MMA. The low degree of wastewater treatment coverage might be a factor of increased risk for COVID-19 pandemic.
Two out of 16 samples from 3 different surface water reservoirs tested positive to SARS-CoV-2, with differences between the first and second campaign. The lack of the presence of the virus in the first campaign is related to a decrease in the reported cases of infection during November due to lockdown measures. With values of 3.3 and 3.8 GC/ml, the viral load is a bit lower than in urban rivers.
Twenty-two out of 50 groundwater samples exhibited detectable viral loads that varied between 2.9 and 5.6 GC/ml. This share is unexpected high and can be attributed to a contribution from untreated sewage and/or leaky sewage pipes to the groundwater system. The significant correlation between SARS-CoV-2 loads and sucralose concentration in groundwater reaffirms the possibility of leaching and infiltrations of effluents from surface and/or failing sewage pipes. These findings highlight the importance of water disinfection processes in drinking water facilities. Even if it is not probable that freshwater sources may be highly viable for transmission of SARS-CoV-2, the possibility via fecal-oral transmission cannot be discarded at this point. Moreover, differences of viral loads in groundwater between the municipalities inside the metropolitan area are observed; they appear to be consistent with the trends of reported daily cases of infections in each municipality. This shows that groundwater has potential for epidemic surveillance and should be considered as an alternative.
Future investigations should be oriented to the understanding of the fate of the virus in different water environments, factors which are favorable and unfavorable for its development, and the potential exposure risk of water contaminated with this virus in its different variants. Experimental research must go hand in hand with the development of conceptual and computational analysis. It is imperative to develop standard sampling, extraction, isolation, detection, and quantification procedures of the virus in aquatic environments to avoid discrepancies in the results.
Data Availability
NA
Acknowledgements
We would like to acknowledge the financial support received from Consejo Nacional de Ciencia y Tecnología (CONACyT) through the Fondo Desarrollo Tecnológico e Innovación COVID-19 grant No. 312558. Thanks to A. Torres-Hernandez and Christian F. Narváez-Montoya for assisting in the field work. Korine M. S. assisted with the graphical design.