David
Lisboa
a,
Devrim
Kaya
ab,
Michael
Harry
a,
Casey
Kanalos
a,
Gabriel
Davis
a,
Oumaima
Hachimi
a,
Shana
Jaaf
a,
David
Mickle
a,
Dana
Alegre
c,
Katherine
Carter
c,
Steven
Carrell
c,
Mark
Dasenko
c,
Nathan
Davidson
c,
Justin
Elser
c,
Matthew
Geniza
c,
Anne-Marie
Girard
c,
Brent
Kronmiller
c,
Matthew
Peterson
c,
Elizabeth
Zepeda
c,
Christine
Kelly
a and
Tyler S.
Radniecki
*a
aSchool of Chemical, Biological and Environmental Engineering, Oregon State University, Corvallis, OR, USA. E-mail: tyler.radniecki@oregonstate.edu
bSchool of Public Health and SDSU-IV, San Diego State University, San Diego, CA, USA
cCenter for Quantitative Life Sciences, Oregon State University, Corvallis, OR, USA
First published on 5th July 2024
The evaluation of COVID-19 policy effectiveness on university campuses, particularly in mitigating spread to neighboring cities (i.e., “campus spill-over”), is challenging due to asymptomatic transmission, biases in case reporting, and spatial case reporting limitations. Wastewater surveillance offers a less biased and more spatially precise alternative to conventional clinical surveillance, thus providing reliable data for university COVID-19 policy evaluation. Wastewater surveillance data spanning the academic terms from Fall 2020 through Spring 2022 was used to evaluate the impact of university COVID-19 policies. During the campus closure to external visitors (09/21/2020–9/15/2021), campus viral concentrations and variant compositions were dissimilar from those of the host and neighboring cities (MAPE = 0.25 ± 0.14; Bray–Curtis = 0.68 ± 0.1, respectively), indicating relative isolation of the campus from its surroundings. Upon the campus reopening to visitors (9/15/2021–2/27/2022), the viral concentrations and variant compositions matched more closely with the host and neighboring cities (MAPE = 0.21 ± 0.1; Bray–Curtis = 0.14 ± 0.08, respectively). Furthermore, post-lifting of campus and state mask mandates (2/27/2022–6/12/2022), the campus, host and neighboring city viral concentrations and variant compositions became indistinguishable (MAPE = 0.06 ± 0.02; Bray–Curtis = 0.07 ± 0.05, respectively). This data suggests that university COVID-19 policies effectively prevented campus-spill over, with no significant contribution to COVID-19 spread into the surrounding communities. Conversely, it was the surrounding communities that led to the spread of COVID-19 onto the campus. Therefore, wastewater surveillance proves instrumental in monitoring COVID-19 trends in surrounding areas, aiding in predicting the impact of easing campus restrictions on campus health.
Water impactThis study utilized wastewater surveillance to assess the efficacy of university COVID-19 policies to protect campus and community health. The findings demonstrated that university COVID-19 policies were effective at preventing campus spill-over to neighboring communities. This research highlights wastewater surveillance's potential as a vital tool in evaluating the effectiveness of university public health policies. |
Traditional clinical surveillance methods, such as PCR testing, contact tracing, and isolating, have been essential in both identifying and isolating cases to prevent COVID-19 spread.5 However, there are significant limitations within this. These include asymptomatic transmission, underreporting, test availability, and time delays due to infection-to-symptom onset.6 Additionally, early COVID-19 testing were shown to disproportionately represent both Black and Hispanic populations due to geographic disparities of testing sites.7–9 Although clinical surveillance is useful in monitoring county and state-wide data, it often fails to capture the whole picture of COVID-19 viral transmission, especially within dynamic populations such as university campuses.
In contrast, wastewater surveillance offers a less biased and more spatially precise alternative.10–12 This method uses PCR-based approaches to quantify pathogens, like the SARS-CoV-2 virus, that are shed in an infected individual's stool and enters the waste stream after flushing. Additionally, pathogen variant composition can be identified via amplicon sequencing.13 This tool has been used to provide an accurate representation of the pathogenic activity within a neighborhood or community, providing an additional source of data to assist public health officials in making informed decisions on how to respond to the pandemic dynamics.14–18 Wastewater surveillance as a whole is cost-effective, non-invasive, and can be performed in areas where clinical testing is constrained.19,20
Despite its advantages, wastewater surveillance is also linked to potential biases and constraints such as variations in water use, rainwater infiltration and inundation (I and I), and sampling frequency. This study addresses these challenges through regular maintenance of autosamplers and excluding data from significant rainfall events. Additionally, SARS-CoV-2 quantification underwent rigorous quality control measures to mitigate the risk of inaccurate data.
Wastewater surveillance has been utilized to identify COVID-19 hotspots at the community scale12 and to provide initial screening efforts at the building scale.21 At the university level, wastewater surveillance has been successful in monitoring the viral burden of a campus community.20,22 Furthermore, wastewater surveillance has been effective as a tool to identify locations of positive individuals and provide intervention to contain outbreaks.21,23–25 However, to date, wastewater surveillance has not been used to evaluate the effectiveness of university COVID-19 policies.
A common concern of many living in university host cities was the increased risk of COVID-19 spreading through the community and neighboring communities due to the influx of out-of-area students.26–29 Universities responded to this concern, in addition to concern for the well-being of their students, faculty and staff, with a variety of COVID-19 policies. These policies included campus closures, isolation policies in residential halls for infected students, and mask mandates. However, quantifying the impact of college students entering into a community, as well as the effectiveness of university COVID-19 policies to reduce SARS-CoV-2 transmission, can be challenging using traditional clinical surveillance metrics due to the small-geographic scale of a university campus,26,30 as well as the clinical testing biases mentioned above.
This study utilized wastewater surveillance to investigate the effectiveness of university COVID-19 policies at Oregon State University (OSU) in reducing the spread of COVID-19 on campus. Additionally, it evaluated the effects of rainwater infiltration and inundation (I and I) as well as autosampler maintenance times to make recommendations towards more accurate data collection. Furthermore, wastewater surveillance data of the university host city (Corvallis, OR) and that of a neighboring city (Albany, OR) were compared to the campus community to help determine the effectiveness of university COVID-19 policies in preventing the spread of COVID-19 from the campus into its surrounding communities, (i.e. campus spill-over).
Group | Location | Number of samples | Sampling period |
---|---|---|---|
Main line | East Campus | 127 | 09/21/20–06/04/22 |
International Living-Learning Center (ILLC) | 133 | 09/27/20–06/18/22 | |
Arnold Dining | 142 | 10/04/20–06/11/22 | |
Building clusters | Poling/NW Weatherford | 114 | 09/27/20–06/11/22 |
Goss Stadium | 139 | 10/04/20–06/04/22 | |
Halsell | 109 | 09/27/20–06/11/22 | |
Beth Ray Center | 135 | 10/07/20–06/04/22 | |
Isolated buildings | Sackett | 115 | 09/21/20–06/11/22 |
West | 112 | 09/30/20–06/08/22 | |
Finley | 149 | 09/21/20–06/18/22 | |
Callahan | 123 | 09/21/20–06/08/22 | |
Gem | 118 | 09/27/20–06/11/22 | |
Hawley/Buxton | 118 | 09/21/20–06/08/22 | |
Wilson/McNary/Tebeau | 118 | 09/21/20–06/11/22 | |
SE Weatherford | 112 | 09/30/20–06/11/22 | |
WW treatment plants | Corvallis | 374 | 04/21/20–06/05/22 |
Albany | 151 | 09/23/20–08/12/22 |
Campus sampling sites were divided into three categories: main line, building clusters, and isolated buildings. Main line locations were those in which most of the wastewater effluents from buildings across campus intersected and flowed. Building clusters contained wastewater effluent from at least ten campus buildings which included a mix of residential halls and academic buildings. Finally, isolated buildings were those in which wastewater effluent contributing to the sample was only from residential halls. Isolated buildings were primarily a single residential hall but could also include residential hall groups (e.g. Hawley/Buxton and Wilson/McNary/Tebeau). The International Living-Learning Center (ILLC) and East Campus locations were used to represent the entire campus' wastewater data as the majority of the university's effluent flowed through these points before exiting to the city. The East Campus location was used as the primary source for representing the entire campus while the ILLC location was used when sampling failures occurred at the East Campus location.
SARS-CoV-2 targets (N1 and N2) and Human RNaseP (internal control, RP) were quantified via RT-ddPCR using BioRad's 2019-nCoV CDC ddPCR Triplex Probe Assay and the One-Step RT-ddPCR Advanced Kit for Probes (Bio-Rad Laboratories, Hercules, CA) on a QX-200 ddPCR system with an automated droplet generator and droplet reader, per manufacturer's instructions (Bio-Rad Laboratories). All assay conditions were performed as specified in the Bio-Rad assay protocol with a template concentration of 5.5 μL of RNA per reaction.31 One-step thermal cycling conditions were as follows: reverse transcription at 50 °C for 60 min, enzyme activation at 95 °C for 10 min; 40 cycles of denaturation at 94 °C for 30s and annealing/extension at 55 °C for 60 s; enzyme inactivation at 98 °C for 10 min; droplet stabilization at 4 °C for 30 min to a maximum of overnight.
A minimum of three positive droplets was required for a sample to be identified as positive for SARS-CoV-2. Samples negative for SARS-CoV-2 detection were assigned a concentration of ½ the limit of detection. All samples were run in duplicate and the N1 and N2 concentrations were averaged together for a final concentration per sample. As determined by Layton et al., the limit of blank for N1 and N2 were 2.0 and 4.2 copies per reaction and the estimated limit of detection for N1 and N2 were 4 and 12 copies per reaction, respectively.16 Additionally, the process recovery efficiency (using bovine coronavirus as a surrogate) was determined to be 57%.
Each ddPCR plate contains the following controls in duplicate; field blanks (FBs), extraction blanks (EBs), negative control reactions (containing human genomic DNA), positive control reactions (containing synthetic RNA of N1 and N2 and synthetic DNA of RP), and no template controls (NTCs). For a ddPCR run to be considered valid, a minimum of 6000 droplets per reaction must have been produced, all negative controls (FBs, EBs, NTCs and negative control reactions) must have less than 3 positive droplets, and all positive controls (synthetic RNA of N1 and N2 and synthetic DNA of RP) must have passed the three-positive droplet threshold. Per manufacturer's instructions, for an individual ddPCR reaction to be considered valid, at least one of the targets (N1, N2 or RP) must have passed the three-positive droplet threshold. Reactions that did not pass this threshold were excluded from further analyses.
Sequence reads were demultiplexed, trimmed to the reference sequence (Wuhan-Hu1, GenBank accession no. NC_045512.2), and coordinate-sorted with SAMtools version 1.10 (Genome Research Limited, http://www.sanger.ac.uk). The primer sequences were trimmed and GATK version 4.2.0.0 (Broad Institute, https://www.broadinstitute.org) was used to identify mutations compared to the reference genome sequence.
A multilocus sequence typing approach was then used to identify SARS-CoV-2 variants by matching amplicon mutations to known SARS-CoV-2 variant mutation sequences. For a positive identification of a variant, a lower limit of 5% of sequence reads (with a minimum of six total reads) was required to span the mutation site. Additionally, a minimum of two different mutations unique to each variant was required.
To determine the effects of rainwater inundation and infiltration (I and I), spearman correlations were conducted between 48 h rainfall totals and wastewater conductivity. Additionally, regression analyses were performed comparing East Campus wastewater SARS-CoV-2 concentrations to reported on-campus COVID-19 cases under two conditions: 1) using all samples regardless of rainfall and 2) removing samples if the 48 h rainfall total was greater than 15 mm. This cutoff was chosen using the IQR method of outlier detection.
During the 2021–2022 academic year, the OSU Corvallis Campus was opened to outside visitors, non-essential employees, and non-resident hall students, and policies followed guidelines provided by the CDC. These requirements included: receiving the first set of vaccinations before arriving on campus, weekly testing for those exempt from vaccination, wearing face coverings at all times in shared spaces, and maintaining a distance of 6 feet at all times.34 During Spring 2022, the mask mandate was lifted, and all other policies remained in place. During the 2021–2022 academic year, the OSU Corvallis Campus was opened to outside visitors, non-essential employees, and non-resident hall students. Classes were held in person and the resident hall population increased to ∼4700 students across the 13 residential halls. A state-wide mask mandate remained in place until March 12, 2022, requiring masks to be worn inside.
![]() | (1) |
In the comparison of rainfall and wastewater conductivity, the interquartile range (IQR) method of outlier detection was used to identify and filter out outliers that would skew the data. This method defines an outlier as a data point in which its value is 1.5 times greater than or 1.5 times less than the IQR.40
Additionally, the Fisher's Z-transformation was utilized when comparing between different correlations. This test transforms the samples to become normally distributed and then returns a “z-score” that allows for the testing of the significance of the difference between two correlation coefficients.41 The Fisher's Z-transformation was used to assess whether rainfall events resulted in significant variations in the wastewater data.
![]() | (2) |
y i = log-transformed SARS-CoV-2 concentration of another community
n = number of samples collected
The Bray–Curtis dissimilarity was used to quantify the differences between communities based on the variant composition of the wastewater samples (eqn (3)). Bray–Curtis returns a value in the range of 0–1 where 1 indicates the populations have complete dissimilarity from each other.
![]() | (3) |
S i = sum of total variant percentage counted in the community i
S j = sum of total variance percentage counted in community j
C ij = the sum of only the lesser counts for each variant in both sites
![]() | (4) |
σ = standard deviation
= sample mean
These failure events impacted the data collected by creating gaps in the data, potentially leading to an underrepresentation of the viral load during those times. To mitigate this impact, data from the days with failure events (7 out of 144 samples) were excluded from this study to ensure that the analysis accounted for these missing points. Thus, autosamplers used in campus wastewater surveillance efforts should undergo preventative maintenance approximately every 4–5 weeks to minimize risk of these disruptions and biases.
![]() | ||
Fig. 1 Correlation between conductivity and rainfall at different sampling points across campus. * = p < 0.05, ** = p < 0.01, *** = p < 0.001. |
The removal of samples with 48 h rainfall totals greater than 15 mm increased the correlation from rs(71) = 0.67 to rs(63) = 0.71 (Fig. 2A and B). However, a Fisher's r to z transformation identified that these values were not significantly different from each other (p > 0.05). Additionally, 71% of sampling events occurred during periods of little-to-no precipitation, and more than 90% of the rainfalls were less than 15 mm over 48 hours. Thus, it was concluded that rainfall events were not of significant consequence to this study, and no further transformation of the data was necessary.
It should also be noted that while the influence of I and I on wastewater surveillance results is often ignored, it is also highly site specific. For instance, in new housing developments with new separated sewer systems, the influence of I and I on wastewater surveillance would be negligible. However, when conducting wastewater surveillance in established communities with combined sewer systems or older conveyance lines, the influence of I and I can be substantial. Especially if a sizable area drains to that location. Thus, we recommend that the influence of I and I be investigated at all sampling locations through identifying potential sources of I and I as well as through measuring wastewater parameters, such as conductivity, that may indicate if I and I is occurring during the sampling period.
Additionally, during the campus closure, the viral wastewater concentrations at the East Campus location were weakly correlated with those in Corvallis and Albany, showing r values of 0.31 and −0.18 (Fig. 3A and D) and MAPE values of 0.87 and 1.16, respectively (Fig. 4B). Notably, except for the period from January 17, 2021, to March 14, 2021, the East Campus wastewater viral concentrations were generally lower than those observed in Corvallis or Albany. Meanwhile, during the campus closure, the correlation of wastewater viral concentrations between Corvallis and Albany remained moderate, with a rs value of 0.50 and a MAPE of 0.39 (Fig. 3G and 4B), indicating some interaction between these communities. Finally, during the campus closure, the wastewater SARS-CoV-2 variant composition remained distinct for the East Campus location compared with the Corvallis and Albany wastewater SARS-CoV-2 variant compositions with a Bray Curtis dissimilarity index of 0.63 and 0.73, respectively (Fig. 5C).
The high correlation between East Campus SARS-CoV-2 wastewater concentrations and the reported campus COVID-19 cases suggests that the wastewater accurately captured the campus disease burden. Additionally, the moderate correlation between the campus and Corvallis wastewater SARS-CoV-2 concentrations in combination with the low correlation between the wastewater SARS-CoV-2 concentrations and variant composition with the cities of Corvallis and Albany suggests that the campus was insular from the outside communities and did not strongly influence the adjacent communities. This is unlike Corvallis and Albany themselves, which have moderate correlations between wastewater SARS-CoV-2 concentrations and variant compositions suggesting that there was an influence of these two communities on one another.
These results are consistent with other non-wastewater studies that estimated the effects of opening the campus to the surrounding community. Arnold et al.,26 conducted a clinical study that quantified the presence of anti-SARS-CoV-2 antibodies of Pennsylvania State University residents and those residing in the surrounding neighborhoods during the time students were allowed to return to campus. There was little evidence of a significant increase in COVID-19 cases in the surrounding community following the return of students to campus. Similarly, Valesano et al.,28 conducted a genomic analysis of SARS-CoV-2 in positive individuals both on-campus and within the community. This study found that most of the variants present in the outside community were not linked to variants within the campus population. Thus, these findings collectively demonstrate minimal evidence of a campus spillover into the surrounding communities.
The weakened correlation between reported campus COVID-19 cases and East Campus wastewater SARS-CoV-2 concentrations was indicative of the campus community becoming less insular from the surrounding communities. In particular, the shallower slope between wastewater SARS-CoV-2 concentrations and reported campus COVID-19 cases (i.e. fewer reported COVID-19 cases per wastewater SARS-CoV-2 concentration) indicated that infected individuals arrived on campus, contributed to SARS-CoV-2 signals in the wastewater but were not reported as campus-based COVID-19 infections, since their residences were located off-campus.
Similarly, the increased correlation between East Campus and Corvallis wastewater SARS-CoV-2 concentrations suggests a greater interaction between these two communities while the improved, but still weak, correlation between East Campus and Albany SARS-CoV-2 wastewater concentrations suggests that these two communities are still largely insular from one another. Additionally, the increased similarity of the wastewater SARS-CoV-2 variant composition also indicates increased interactions between the three communities (Fig. 5). However, the rapid rise of the Delta variant being dominant across the state at this time weakens the usefulness of this comparison.16 Nonetheless, various Omicron variants were detected in the Corvallis wastewater before the East Campus wastewater, indicating that these variants did not enter the city communities through the campus community members (Fig. 5B).
Taken as a whole, the increased correlation between East Campus and Corvallis indicates that either East Campus is influencing Corvallis or vice versa. However, the increased correlation between Corvallis and Albany combined with the weakened correlation between East Campus and Albany suggests that there was not a campus spill-over effect into either Corvallis or Albany. Instead, this data indicates the exact opposite with Corvallis spilling over onto campus. These findings coincide with a study performed by Platt et al.,4 where in-class instruction with a mask mandate did not increase viral transmission on campus, but students who had become infected had potential exposure points from interactions with the outer community.
This theory of the surrounding communities increasing the COVID-19 burden on campus is also supported by the transmission rates and vaccination coverage reported in the Corvallis (Benton County) and Albany (Linn County) communities during this time. The non-detections of SARS-CoV-2 in the East Campus wastewater at the start of the campus opening period indicated a low level of COVID-19 transmission on campus (Fig. 4A). This was likely due to the near 100% vaccination level of university students, faculty, and staff at that time, as required by university policy. In contrast, at this time the COVID-19 transmission rates in Corvallis (Benton County) and Albany (Linn County) were classified as “high” to “substantial” and the single dose vaccination rates in Corvallis and Albany, were much lower at rates of 68% and 49%, respectively (Fig. S3†).35,37 Thus, the introduction of unvaccinated individuals from communities with relatively high transmission rates onto campus may have likely contributed to the onto-campus spill-over effect.
The further increase in correlation between East Campus wastewater SARS-CoV-2 concentrations with those from both Corvallis and Albany suggests an even greater interaction between the three communities after the mask mandate was lifted (Fig. 4). At this time, transmission rates in both communities remained classified as “high” to “substantial” and the single dose vaccination rates in Corvallis and Albany were still lower than the university campus at 80% and 60%, respectively (Fig. S3†). Similar to before the mask mandate was lifted, the stronger correlation between Corvallis and Albany compared to East Campus indicates that campus spillover was not occurring but rather the surrounding communities were spilling over into the campus. Finally, the wastewater SARS-CoV-2 variant composition of East Campus, Corvallis, and Albany also increased in similarity after the mask mandate was lifted (Fig. 5). However, just as with Delta, a single variant, Omicron, was dominant across the state at this time weakening the usefulness of this comparison.16
This study also demonstrated that the enactment of campus closures and mask mandates were highly effective in reducing SARS-CoV-2 viral burdens and dampened the influence of campus residents on the surrounding communities and vice versa. Additionally, the findings of this study indicate that the campus spill-over hypothesis can be dispelled, showing that the increase in the SARS-CoV-2 viral load on the campus primarily originated from the surrounding community rather than from the campus spreading the virus to the community.
As the campus COVID-19 restrictions on campus were adjusted, the viral dynamics on the campus began to mirror those of the host city. This suggests that university policymakers can use the viral activity patterns of the host city as a benchmark to estimate the consequences of lifting COVID-19 restrictions on campus. Additionally, through wastewater surveillance's ability to track viral trends down to the neighborhood and building scale, this method could be utilized to create targeted interventions that alleviate transmission potential from building to building. On a university-scale, this could include enforcing testing and isolation of that area, without having to spend resources to test the entire campus. Overall, wastewater surveillance in conjunction with traditional clinical surveillance methods, can serve as a valuable tool towards refining policy decisions based on accurate, real-time data.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ew00168k |
This journal is © The Royal Society of Chemistry 2025 |