Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA

Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA | medRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-P4HH5NV'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA View ORCID Profile Elana M. G. Chan , View ORCID Profile Alessandro Zulli , View ORCID Profile Alexandria B. Boehm doi: https://doi.org/10.1101/2025.07.21.25331943 Elana M. G. Chan 1 Department of Civil and Environmental Engineering, Stanford University , 473 Via Ortega, Stanford, California 94305, United States Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Elana M. G. Chan Alessandro Zulli 1 Department of Civil and Environmental Engineering, Stanford University , 473 Via Ortega, Stanford, California 94305, United States Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Alessandro Zulli Alexandria B. Boehm 1 Department of Civil and Environmental Engineering, Stanford University , 473 Via Ortega, Stanford, California 94305, United States Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Alexandria B. Boehm For correspondence: aboehm{at}stanford.edu Abstract Full Text Info/History Metrics Supplementary material Data/Code Preview PDF Abstract Surveillance of rotavirus infections remains critical because vaccines are underutilized in the USA. Using wastewater solids measurements of rotavirus RNA collected over one year from 185 wastewater treatment plants (WWTPs) in the USA, we inferred spatiotemporal occurrence patterns of rotavirus infections and compared occurrence patterns to clinical metrics of infections and markers of vaccination coverage. We also estimated infection prevalence from wastewater measurements using available data on rotavirus RNA shedding in feces. Nationally, wastewater measurements of rotavirus RNA were correlated with clinical metrics of infection and exhibited elevated winter-spring concentrations beginning in the South. WWTP service areas characterized by markers of high vaccination coverage generally experienced a shorter duration of elevated rotavirus concentrations compared to areas characterized by markers of low vaccination coverage. Rotavirus infection prevalence estimates were highly uncertain and sensitive to shedding parameters. Wastewater monitoring of vaccine-preventable diseases is valuable for informing where vaccination campaigns should be targeted. Introduction Rotavirus infection remains the leading cause of diarrheal deaths worldwide, particularly affecting individuals under 5 and over 70 years, despite available vaccines. 1 Globally, rotavirus deaths declined between 1990 and 2019 but there was a significant uptrend in Canada and the United States—despite North America having some of the highest levels of vaccination coverage in the world. 1 , 2 Still, vaccines have alleviated the burden of rotavirus disease. About 2.7 million cases were reported each year in the United States, but vaccines have since annually averted 280,000 clinic visits; 62,000 emergency department visits; and 45,000 hospitalizations. 3 Continued surveillance of infections is necessary to understand the epidemiology of rotavirus in the post-vaccination era and to identify populations inadequately covered by vaccines. However, rotavirus disease is not nationally notifiable in the United States, 4 so comprehensive clinical surveillance data remains limited. Two live-attenuated, oral rotavirus vaccines are currently licensed in the United States. RotaTeq (Merck and Company, Whitehouse Station, New Jersey) was licensed in 2006 and is a three-dose series. 5 Rotarix (GlaxoSmithKline Biologicals, Rixensart, Belgium) was licensed in 2008 and is a two-dose series. 5 Both vaccines are administered before 8 months of age and highly effective at preventing severe infection during the first two years of life. 5 , 6 Whether vaccines confer lifelong immunity is unknown, but waning effectiveness has been reported in some children. 6 , 7 Older children and adults born prior to licensure of vaccines may experience repeat infections throughout their lives, although with reduced disease severity, as immunity from natural rotavirus infection is incomplete. 8 – 10 Nonetheless, studies have shown that rotavirus vaccines provide some indirect protection to unvaccinated populations. 11 , 12 Despite the important direct and indirect benefits of rotavirus vaccines, vaccines are underutilized in the United States. Coverage for recommended childhood vaccines decreased in the United States during the COVID-19 pandemic. 13 Estimated rotavirus vaccination coverage was 75.1% in the 2020–2021 birth cohort, which falls below the Healthy People 2020 target of 80% and lags behind other vaccines for the same birth cohort (e.g., poliovirus: 91.1%; measles, mumps, and rubella: 90.3%; varicella: 89.9%). 13 , 14 Several geographic, sociodemographic, and healthcare characteristics are documented to influence rotavirus vaccination coverage rates. 14 – 16 Still, vaccines have impacted rotavirus seasonality in the United States according to available clinical testing data. Prior to vaccines, rotavirus followed an annual, winter–spring occurrence pattern. 4 , 17 Following vaccine introduction, seasonality shifted to a biennial occurrence pattern, with low incidence during even years and elevated incidence during odd years. 4 , 17 Symptomatically and asymptomatically infected individuals shed rotavirus RNA in feces and other excreta; 18 studies have demonstrated that community-level monitoring of rotavirus infections can be conducted via wastewater. 19 – 24 Monitoring rotavirus in wastewater is complementary to clinical surveillance because clinical testing has declined since vaccine introduction, is biased towards those with severe disease and access to healthcare, and is reported on a voluntary basis by healthcare providers. 4 In contrast, a wastewater sample represents an unbiased composite biological sample of an entire contributing population and measurements can be available within 24 hours. Wastewater measurements are inclusive of individuals with self-limiting symptoms who are less likely to seek clinical care, such as adults, and may better represent population-level disease occurrence. 8 , 9 In this study, we aimed to characterize spatiotemporal occurrence patterns in rotavirus disease across the United States using wastewater measurements of rotavirus RNA. Samples were collected from 185 wastewater treatment plants (WWTPs) multiple times a week over approximately one year. We compared occurrence patterns of rotavirus disease inferred from wastewater to clinical metrics of rotavirus infections, estimates of rotavirus vaccination coverage, and population characteristics associated with rotavirus vaccination coverage. As a secondary aim, we estimated rotavirus infection prevalence using wastewater rotavirus RNA concentrations and available data on shedding of rotavirus RNA in human feces. Although previous studies have linked rotavirus RNA concentrations in wastewater to clinical disease occurrence (as measured by clinical test positivity rates and reported case numbers), 19 – 24 those studies were conducted over small geographic scales and none of them compared rotavirus RNA concentrations to markers of vaccination coverage to assess the utility of wastewater monitoring for identifying vaccination coverage disparities. Results Wastewater samples were typically collected three times per week from 185 WWTPs across 40 states and the District of Columbia between 3 September 2023 and 12 December 2024. Across a total of 33,471 wastewater samples, measured rotavirus RNA concentrations ranged from non-detect to 8.4 × 10 7 gene copies per gram (gc/g), with a median concentration of 5.1 × 10 4 gc/g and interquartile range (IQR) of 8.5 × 10 3 –2.4 × 10 5 gc/g. We detected pepper mild mottle virus (PMMoV) RNA in all wastewater samples, suggesting adequate viral recovery throughout the analysis period ( Figure S1 ). Occurrence of rotavirus wastewater events across the United States We examined patterns in rotavirus RNA concentrations in wastewater at the WWTP, state, and national scales to identify periods when concentrations were elevated, suggesting higher than background levels of rotavirus infections; these periods are referred to as rotavirus “wastewater events”. Some locations had more than one wastewater event, so the longest wastewater event is the one described below. Nationally aggregated, PMMoV-normalized wastewater measurements of rotavirus RNA were variable over the study period, with relatively high values observed during the single national wastewater event occurring between 1 January 2024 and 23 June 2024 (duration: 173 days) ( Figure 1a ). An example rotavirus wastewater event at the WWTP scale is shown in Figure 1b for San Jose, California (duration: 316 days). PMMoV-normalized rotavirus RNA wastewater measurements and events for all other WWTPs and states are included in the Supporting Information (SI) ( Figures S2–S3 ). At the WWTP scale, wastewater event durations ranged from 20–444 days (median: 182 days); two WWTPs did not experience any wastewater event. At the state scale, the maximum wastewater event duration occurred in Utah (306 days) and the minimum duration occurred in Vermont (14 days); the median duration was 150 days ( Figure S4 ). Download figure Open in new tab Figure 1. Wastewater rotavirus RNA concentrations and events. ( a ) Nationally aggregated, smoothed, PMMoV-normalized wastewater rotavirus RNA concentrations with the national rotavirus wastewater event shaded. ( b ) Smoothed, PMMoV-normalized wastewater rotavirus RNA concentrations with the rotavirus wastewater event shaded for an individual WWTP (San Jose, California). ( c ) Count of states in a wastewater event each day, grouped by census region. ( d ) Smoothed, PMMoV-normalized wastewater rotavirus RNA concentrations across all samples collected between September–December in 2023 versus 2024. The center line represents the median, the box limits represent the 25th and 75th percentiles, the whiskers represent 1.5 x IQR from the box limits, and points represent values beyond the whiskers. Concentrations are displayed on a log 10 scale, and the horizontal bar denotes a significant difference. Abbreviations: IQR = interquartile range, PMMoV = pepper mild mottle virus, WWTP = wastewater treatment plant. Figure 1c displays the number of states experiencing a wastewater event each day, stratified by census region. Wastewater events appeared to onset first in the South followed by the Midwest, West, and Northeast. At least half of states in the region were in onset for 21 consecutive days beginning 10 December 2023 in the South, 1 January 2024 in the Midwest, 25 January 2024 in the West, and 7 February 2024 in the Northeast. Wastewater events appeared to offset first in the South followed by the Northeast, Midwest, and West. Less than half of states in the region were in onset for 21 consecutive days following onset beginning 24 May 2024 in the South, 27 June 2024 in the Northeast, 1 July 2024 in the Midwest, and 10 July 2024 in the West. Across all wastewater samples collected between September–December, PMMoV-normalized wastewater measurements of rotavirus RNA were higher in 2023 compared to 2024 (Wilcoxon test W = 39,152,179; p < 0.0001; n = 8,929 in 2023; n = 6,278 in 2024) ( Figure 2d ). Download figure Open in new tab Figure 2. Rotavirus RNA wastewater concentrations and occurrence of rotavirus infections across the United States. Weekly smoothed, PMMoV-normalized rotavirus RNA concentration in wastewater (left y-axis; blue) and weekly clinical surveillance metric of rotavirus infections (right y-axis; red) for the United States. ( A ) Percent of rotavirus antigen tests reported as positive each week from the National Respiratory and Enteric Virus Surveillance System (NREVSS) is used as the clinical metric. ( B ) Patients per 100 million with a rotaviral enteritis encounter diagnosis each week from Epic Cosmos is used as the clinical metric. Wastewater measurements of rotavirus RNA compared to clinical metrics of rotavirus infection Weekly PMMoV-normalized rotavirus RNA wastewater measurements were correlated with weekly rotavirus clinical metrics at the national scale ( Figure 2 ). The correlation coefficient was larger between wastewater measurements and the proportion of patients with a rotaviral enteritis encounter diagnosis from Epic Cosmos (Kendall’s tau = 0.84 [95% CI: 0.80, 0.88], p < 0.0001; n = 67) versus rotavirus test positivity rate from the National Respiratory and Enteric Virus Surveillance System (NREVSS) (Kendall’s tau = 0.49 [95% CI: 0.36, 0.63], p < 0.0001; n = 67). Rotavirus wastewater event durations compared to estimates of rotavirus vaccination coverage Rotavirus wastewater event durations significantly differed across tertiles of vaccination coverage at the sewershed scale using vaccination estimates from Epic Cosmos (χ 2 = 7.04, p = 0.030); the effect size was small (η 2 = 0.028) with a significant difference between the middle and highest tertile ( Figure 3 ). At the state scale, rotavirus wastewater event durations did not significantly differ across tertiles of vaccination coverage using rotavirus vaccination estimates from both the US Centers for Disease Control (CDC) National Center for Immunization and Respiratory Diseases (NCIRD) (χ 2 = 0.67, p = 0.72) and Epic Cosmos (χ 2 = 2.14, p = 0.34) ( Figure S5 ). The state-level rotavirus vaccination estimates from NCIRD and Epic Cosmos were not correlated ( Figure S6 ). Download figure Open in new tab Figure 3. Rotavirus wastewater event durations among sewersheds by tertile of rotavirus vaccination coverage. Distribution of rotavirus wastewater event durations among sewersheds by tertile of rotavirus vaccination coverage in the predominant county of the sewershed. The center line represents the median, the box limits represent the 25th and 75th percentiles, and the whiskers represent 1.5 x IQR from the box limits. Jittered points represent individual data points. Rotavirus vaccination coverage estimates from Epic Cosmos ranged from 4.69–84.6% among sewersheds. Horizontal bar denotes significant pairwise differences. Abbreviations: IQR = interquartile range. Rotavirus wastewater event durations compared to population characteristics associated with rotavirus vaccination coverage Ghaswalla et al. 14 identified several population characteristics—pertaining to social, geographical, family, and healthcare factors—associated with vaccination coverage in the United States. The distribution of wastewater event durations among sewersheds by population characteristic tertile or census region is shown in Figure 4 ; refer to the SI for tertile cutoffs ( Figures S7–S10 ). Horizontal bars denote significant pairwise differences following Kruskal Wallis tests with p < 0.05. Refer to Table 1 for chi square test statistics, p values, and eta-squared effect size estimates for each Kruskal-Wallis test. We did not observe a significant difference across tertiles for Hispanic or Latino children or pediatrician-to-child ratio. We observed a significant difference across tertiles for publicly insured children; however, the effect size was small (η 2 = 0.024) and no significant pairwise differences were identified. We observed a significant difference across groups for all other characteristics with a small effect size for foreign-born children (η 2 = 0.044), children below poverty (η 2 = 0.044), average family size (η 2 = 0.034), and urbanicity (η 2 = 0.053); a moderate effect size for Black or African American children (η 2 = 0.073), uninsured children (η 2 = 0.074), and children born to younger mothers (η 2 = 0.098); and a large effect size for census region (η 2 =0.14). Download figure Open in new tab Figure 4. Rotavirus wastewater event durations by population characteristics. Distribution of rotavirus wastewater event durations among sewersheds by population characteristic tertiles ( a–j ) or census region ( k ). The center line represents the median, the box limits represent the 25th and 75th percentiles, and the whiskers represent 1.5 x IQR from the box limits. Jittered points represent individual data points. Horizontal bars denote significant pairwise differences. Refer to the Supporting Information for tertile cutoffs of population characteristics in a–j. Abbreviations: IQR = interquartile range. View this table: View inline View popup Download powerpoint Table 1. Association between rotavirus wastewater event duration and population characteristics at the sewershed scale Estimated population shedding rotavirus RNA from measurements of rotavirus RNA in wastewater solids We adopted a virus-independent mass balance model established by Wolfe et al. 25 to estimate the population fraction shedding rotavirus RNA into wastewater each day (F shed ). The model uses PMMoV-normalized measurements of rotavirus RNA in wastewater solids (C rota_ww /C PMMoV_ww ) and fecal shedding concentrations for PMMoV RNA (C PMMoV_feces ) and rotavirus RNA (C rota_feces ) as inputs. Arts et al. 26 and Zheng et al. 27 provide large datasets on fecal shedding which were used to characterize the distributions of C PMMoV_feces and C rota_feces , respectively. The distributions of C PMMoV_feces and C rota_feces were reasonably modeled as log 10 -normal distributions upon examination of quantile-quantile plots. The mean and standard deviation of log 10 -transformed concentrations of PMMoV RNA measured in feces in gc/g from Arts et al. 26 is 7.94 and 1.57 log 10 gc/g, respectively. The mean and standard deviation of log 10 -transformed concentrations of rotavirus RNA measured in feces in gc/g from Zheng et al. 27 is 7.01 and 2.16 log 10 gc/g, respectively. Figure S11 shows F shed expressed as a percentage modeled from C rota_ww /C PMMoV_ww . Daily estimates of the percentage of people shedding rotavirus RNA for an individual WWTP (San Jose, California) and across all WWTPs for the United States are shown in Figure 5 . At the San Jose, California WWTP, C rota_ww /C PMMoV_ww ranged from 0.000011 to 0.0018, corresponding to a percentage of people shedding rotavirus RNA of 0.0094% (25th percentile = 0.00015%, 75th percentile = 0.61%) and 1.5% (25th percentile = 0.026%, 75th percentile = 110%). Nationally aggregated C rota_ww /C PMMoV_ww across all WWTPs in the study ranged from 0.000053 to 0.0031, corresponding to a percentage of people shedding rotavirus RNA of 0.040% (25th percentile = 0.00059%, 75th percentile = 2.5%) and 3.0% (25th percentile = 0.046%, 75th percentile = 184%). Download figure Open in new tab Figure 5. Daily estimates of the percentage of the population shedding rotavirus RNA. ( a ) Estimates for an individual WWTP (San Jose, California). ( b ) National estimates across all WWTPs. Solid line represents the median estimate; dotted lines represent the 25th and 75th percentile estimates. Abbreviation: WWTP = wastewater treatment plant. The sensitivity analysis suggests the model is highly sensitive to all input parameters. We calculated F shed for each parameter at its 75th percentile (F shed,75 ) and 25th percentile (F shed,25 ) holding all other parameters at their 50th percentile. The F shed,75 :F shed,25 ratio was 0.0013 for C rota_feces , 987 for C PMMoV_feces , and 100 for C rota_ww /C PMMoV_ww . The fraction of the population shedding rotavirus RNA decreases as the concentration of rotavirus in feces increases, increases as the concentration of PMMoV in feces increases, and increases as the PMMoV-normalized rotavirus RNA concentration in wastewater increases. Discussion Rotavirus disease is not nationally notifiable in the United States, yet surveillance of rotavirus infections is critical for identifying gaps in rotavirus vaccination coverage and understanding rotavirus epidemiology in the post-vaccination era. We used wastewater RNA concentrations of rotavirus collected over more than one year across the United States to infer rotavirus disease occurrence in the United States and compare occurrence patterns to clinical metrics of rotavirus infection and markers of rotavirus vaccination coverage. Rotavirus RNA was detected in wastewater samples throughout the study period in varying concentrations. Prior to our study, concentrations of rotavirus RNA in wastewater solids have only been reported in California. 19 Rotavirus has been detected in wastewater sludge in several countries, 28 but our study is the first to report concentrations in wastewater solids across the United States. Boehm et al. 19 retrospectively measured rotavirus RNA concentrations in wastewater solids at two WWTPs in California between February 2021–April 2023 and observed median concentrations on the order of 10 3 gc/g. Herein, we observed median concentrations of rotavirus RNA in wastewater solids on the order of 10 4 gc/g based on samples from 185 WWTPs across 40 states and the District of Columbia between September 2023– December 2024. Studies that measured rotavirus RNA concentrations in raw sewage observed seasonal occurrence patterns, as we observed in our study, with concentrations ranging on the order of 10 2 –10 8 gc/L and 1–10 4 fluorescent foci/L. 29 – 34 Nationally, a rotavirus wastewater event occurred 1 January to 23 June which aligned with the winter–spring occurrence pattern of rotavirus infections according to clinical testing data. 4 , 17 Wastewater measurements of rotavirus RNA were higher in fall 2023 (leading up to an even-numbered season) than in fall 2024 (leading up to an odd-numbered season) across the United States. We did not have longitudinal wastewater measurements from two consecutive winter– spring seasons, but our observations suggest a deviation from the typical biennial occurrence pattern of rotavirus. According to clinical testing data, odd-numbered seasons are generally associated with elevated incidence compared to even-numbered seasons in the post-vaccination era. 4 , 17 The COVID-19 pandemic impacted the biennial pattern of rotavirus occurrence. 35 Mathematical modeling suggests rotavirus transmission will not return to pre-pandemic equilibrium before 2030 which may explain why wastewater rotavirus RNA measurements seemingly did not align with the biennial occurrence pattern. 35 Continued wastewater monitoring of rotavirus RNA can be a useful tool for empirically understanding how temporal patterns in rotavirus occurrence continue to develop. The spatial progression of rotavirus occurrence in wastewater lacked national coherency. Elevated occurrence of rotavirus in wastewater began first in the South followed by the Midwest, West, and Northeast. Prior to vaccines, annual rotavirus occurrence also lacked national coherency, beginning in the Southwest and ending in the Northeast. 36 , 37 In the post-vaccination era, Burnett et al. 38 suggest areas with high birth rate and low vaccination coverage have the most rapid accumulation of susceptible individuals and consequently drive annual rotavirus occurrence. Indeed, rotavirus occurrence began in Oklahoma, Arkansas, and the western Gulf coast—all areas in the South with high birth rate and low vaccination coverage—according to clinical testing data across five seasons in the post-vaccination era. 38 Wastewater monitoring of rotavirus RNA over consecutive seasons can be used alongside clinical testing data to characterize patterns in the spatial progression of rotavirus in the post-vaccination era. Weekly wastewater measurements of rotavirus RNA were positively correlated with clinical metrics of rotavirus infections obtained from both NREVSS and Epic Cosmos at the national scale. Although the two clinical metrics themselves were positively correlated at the national scale ( Figure S12 ), the two datasets are distinct. Epic Cosmos is a research dataset that aggregates electronic health record data from 295 million patients. 39 , 40 In contrast, NREVSS relies on voluntary reporting of rotavirus tests; fewer than 100 rotavirus test results were submitted to NREVSS each week during our analysis period. As shown in Figure 3 , test positivity from NREVSS was visually prone to more week-to-week variability compared to patient diagnoses from Epic Cosmos. Moreover, data can be accessed at subnational scales from Epic Cosmos unlike from NREVSS. 40 For example, we observed wastewater measurements of rotavirus RNA were still positively correlated with clinical metrics of rotavirus infections at the state scale using data from Epic Cosmos ( Figure S13 , Table S1 ). Nonetheless, electronic health record datasets like Epic Cosmos still suffer biases not present in wastewater monitoring data, such as bias towards individuals with symptomatic disease and access to healthcare. Rotavirus encounter diagnoses in Epic Cosmos mostly occur in the under 5 year age group ( Figure S14 ), suggesting infections in older populations are subclinical. Likewise, other studies observed a weak or lack of a positive correlation between rotavirus RNA concentrations in wastewater and traditional measures of disease surveillance, suggesting the presence of asymptomatic or subclinical infections in the contributing community. 19 – 24 Wastewater monitoring represents a less-biased approach for inferring population-level disease occurrence compared to clinical or syndromic measures of disease surveillance. Rotavirus wastewater event durations were generally shorter in sewersheds with markers of high rotavirus vaccination coverage, suggesting the utility of wastewater monitoring for identifying gaps in vaccination coverage. Exceptions include wastewater event durations not being associated with the proportion of Hispanic or Latino children or the pediatrician-to-child ratio and not being highest in the Northeast. 14 However, not all indirect markers of vaccination coverage identified by Ghaswalla et al. 14 were directly available from publicly available sources. For example, we obtained public data on pediatrician-to-child ratio which describes access to pediatricians but not whether children actually seek care from pediatricians over family physicians. We also only conducted bivariate analyses between wastewater event durations and each indirect marker of rotavirus vaccination coverage. Future work might consider combining all markers of rotavirus vaccination coverage into a multivariable model with wastewater event duration as the outcome variable. Additionally, other metrics (e.g., peak wastewater measurement) may be important to consider as the outcome variable. Herein, we were interested in the effect of vaccination on prolonged rotavirus occurrence rather than an isolated outbreak given that vaccines have been documented to shorten rotavirus season. 41 A limitation with our selected dependent variable, however, is that the wastewater event duration may be underestimated for some WWTPs and consequently some states. The wastewater sampling start and end dates differed across WWTPs ( Table S2 ), so sampling may have ended in the midst of a wastewater event for some WWTPs. Nonetheless, the national wastewater event ended before 1 July 2024, which is when a large number of WWTPs ended sampling (mainly in California); it is likely durations were underestimated for only a few WWTPs. Prior to vaccines, approximately 2.7 million cases of rotavirus were reported each year in the United States, 3 corresponding to about 1% of the population based on the population estimate from the 2000 census. 42 Annual rotavirus incidence estimates in the United States are not available for the post-vaccination era but are expected to be substantially lower given the effectiveness of vaccines at preventing severe infection. 5 , 6 Infected individuals may shed rotavirus RNA for consecutive days; 27 therefore, modeled F shed values are likely more representative of disease prevalence rather than incidence. Thus, it is challenging to evaluate the accuracy of our model because rotavirus is monitored in terms of incidence (i.e., reported cases) rather than prevalence. Nonetheless, we expect our model estimates higher rotavirus prevalence than could be determined from case reporting. The age-adjusted prevalence of asymptomatic rotavirus infection is 11% and rotavirus may be detected in about 20–30% of individuals without diarrhea. 43 , 44 Rotavirus infection in adults also usually results in self-limiting symptoms, 8 so adult infections are likely not captured in case reporting. As indicated by our sensitivity analysis, the model is highly sensitive to all parameters. Fecal shedding dynamics are highly individualized, 27 and our model demonstrates that precisely estimating disease occurrence from wastewater measurements is extremely challenging. Infected individuals who shed rotavirus RNA in feces in high concentrations or for long durations (“super shedders”) likely contribute the most rotavirus RNA to wastewater. Future work is needed to better characterize fecal shedding dynamics and refine models so disease occurrence may be estimated from wastewater measurements with less uncertainty. Additionally, rotavirus RNA in wastewater solids was measured using droplet digital reverse transcription-polymerase chain reaction (ddRT-PCR) without first subjecting samples to a “heat snap” (i.e., incubating samples at 99°C for 5 minutes). 45 Rotavirus is a double-stranded RNA (dsRNA), 2 and the additional denaturing step may be important for quantifying RNA of dsRNA viruses using ddRT-PCR. 46 , 47 Including a heat snap step may increase measured rotavirus RNA concentrations in wastewater which could affect modeled F shed reported herein. Our study is the first to assess spatiotemporal occurrence patterns of rotavirus RNA in wastewater across the United States and compare occurrence patterns to markers of rotavirus vaccination coverage. Wastewater monitoring has been used previously to support polio eradication efforts. 48 Wastewater monitoring of other vaccine-preventable diseases may similarly be valuable for identifying disparities in vaccination coverage and informing where vaccination campaigns should be targeted. 49 , 50 Although this study focused on the United States, the global burden of rotavirus is disproportionately highest in Africa, Oceania, and South Asia and rotavirus vaccines have lower performance in countries with high child mortality. 1 , 6 Therefore, wastewater monitoring for rotavirus could be especially impactful in low- and middle-income countries. Continued work on adapting wastewater monitoring methods to decentralized or non-sewered settings is needed to scale wastewater monitoring of rotavirus globally. 51 – 55 Methods Wastewater monitoring Samples were collected from 185 WWTPs across 40 states and the District of Columbia ( Figure 6 ). The number of WWTPs in states ranged from 1–54 (median: 2 WWTPs). The reported population serviced by the WWTPs ranged from 3,000–4,000,000 (median: 89,847). Wastewater samples were typically collected three times per week from WWTPs between 3 September 2023 and 12 December 2024, resulting in 33,471 total samples ( Table S2 ). Following environmental molecular biology best practices, 56 concentrations of rotavirus and pepper mild mottle virus (PMMoV) RNA in the solid fraction of wastewater samples were measured in gene copies per gram (gc/g) using ddRT-PCR. Non-detect rotavirus RNA concentrations were set to 500 gc/g (approximately half the assay detection limit). PMMoV is an indigenous wastewater virus of dietary origin. 57 , 58 Full methodological details and data collected between 1 January 2022 to 30 June 2024 are available in a data descriptor by Boehm et al. 45 Data between 1 July 2024 and 12 December 2024, which were measured using the exact methods reported by Boehm et al. 45 , are uniquely presented in this paper and are available through the Stanford Digital Repository ( https://doi.org/10.25740/mv671wb2446 ). 59 Based on our in silico evaluation described in the SI, it is improbable the rotavirus assay used herein amplifies either the RotaTeq or Rotarix human vaccine strain, animal rotavirus strains, or the bovine rota-coronavirus vaccine spiked into samples as an exogenous control due to the presence of base pair mismatches with the assay primers and/or probe ( Text S1 ). Download figure Open in new tab Figure 6. Map of national wastewater monitoring locations. Locations of WWTPs are symbolized by circles, and states are colored by census region. Created in ArcGIS Pro (version 3.1.1) using 2023 state cartographic boundaries (20-meter resolution) from the US Census Bureau (accessed 5 Aug 2024). 60 Abbreviation: WWTP = wastewater treatment plant. Clinical surveillance We obtained clinical surveillance metrics of rotavirus infections from the National Respiratory and Enteric Virus Surveillance System (NREVSS) and Epic Cosmos (Epic Systems Corporation, Verona, Wisconsin, United States). NREVSS is a laboratory-based sentinel surveillance system that monitors respiratory and enteric virus disease in the United States. Participating laboratories voluntarily report to the Centers for Disease Control (CDC) the number of fecal specimens tested and the number of those specimens that tested positive for rotavirus on a weekly basis. 4 We obtained the publicly available, centered, three-week moving average rotavirus test positivity rate each week (Sunday–Saturday) for the United States for weeks ending 9 September 2023 through 14 December 2024 from the public-facing NREVSS dashboard. 61 The number of specimens tested for rotavirus across the entire NREVSS network during these weeks ranged from 16–80 tests (median: 41 tests). Epic Cosmos is an electronic health record dataset obtained from participating organizations that use Epic Systems. 39 , 40 Cosmos combines data from 295 million patients across over 37,000 clinics and 1,500 hospitals; the dataset is representative of patients who seek healthcare in the United States. 39 We used Cosmos to identify diagnoses of rotaviral enteritis over the study period. First, we geographically filtered the dataset for patients in the United States and temporally filtered the dataset to include patients alive between 1 September 2023 and 14 December 2024. From this subset, we determined the number of patients with a rotaviral enteritis encounter diagnosis each week (Sunday–Saturday). We used these data to calculate the weekly proportion of all patients alive between 1 September 2023 and 14 December 2024 with a rotavirus encounter diagnosis in the United States. Diagnoses were identified using the International Classification of Diseases (ICD-10) code A08.0 for rotaviral enteritis. 62 Vaccination coverage We obtained rotavirus vaccination coverage estimates from the CDC’s National Center for Immunization and Respiratory Diseases (NCIRD) and Epic Cosmos. NCIRD conducts the National Immunization Surveys (NIS) which are phone surveys used to estimate vaccination coverage for vaccinations recommended by the Advisory Committee on Immunization Practices (ACIP). 63 National, regional, state, and selected local area vaccination coverage estimates by birth year and 2- and 4-year birth cohorts are available for rotavirus vaccination from NCIRD. 64 For this analysis, we used state-level rotavirus vaccination coverage estimates from the latest available 2-year birth cohort (2020–2021). To obtain rotavirus vaccination estimates from Epic Cosmos, we first filtered the dataset for patients alive and ≤ 5 years of age in the United States between 1 January 2019 and 12 December 2024. We selected children ≤ 5 years because we aimed to estimate patterns in rotavirus vaccination coverage in the past 5 years. From this subset, we determined the number of patients with any rotavirus immunization encounter grouped by (1) state and (2) county. We used these data to calculate the proportion of all patients alive and ≤ 5 years between 1 January 2019 and 12 December 2024 with a rotavirus immunization encounter at the state and county level. Population characteristics A systematic review on rotavirus vaccination coverage in the United States by Ghaswalla et al. 14 identified social, geographical, family, and healthcare characteristics associated with vaccination coverage. 14 Using publicly available data, we determined the value of 11 variables pertaining to these characteristics for each WWTP sewershed (i.e., the geographic area serviced by a WWTP) ( Table 2 ). 65 – 68 Most WWTPs (n = 113) provided a sewershed boundary. For WWTPs that did not provide a sewershed boundary (n = 72), we approximated a boundary based on the reported zip code(s) serviced by the WWTP and 2023 United States ZIP Code Boundaries published by Esri (Source: TomTom, US Postal Service, Esri). 69 However, sewershed boundaries do not align with legal boundaries (e.g., countries or census tracts) at which publicly available data are often reported. Refer to the SI for details about how we derived sewershed-level estimates for each characteristic using data from each source ( Equations S1–S2 , Table S3 ). View this table: View inline View popup Download powerpoint Table 2. Characteristics associated with rotavirus vaccination coverage in the United States as reported by Ghaswalla et al. 14 Data analysis Unless otherwise stated, we normalized rotavirus RNA concentrations by PMMoV RNA concentrations to correct for differences in virus recovery and wastewater fecal strength and then smoothed the data to reduce effects of outliers. 25 , 45 We smoothed measured PMMoV-normalized concentrations at each WWTP by calculating the five-sample, centered, trimmed, moving average. 45 Trimming removes the highest and lowest value in the moving average window prior to calculating the average. A shrinking window was used to calculate moving average values at the start and end of a time series (e.g., a three-sample window was used to calculate the moving average value for the first sample). Herein, “wastewater measurement” refers to the smoothed, PMMoV-normalized concentration unless otherwise stated. We spatially aggregated wastewater measurements across WWTP sewersheds using a population-weighted averaging approach which has been previously described. 70 – 73 We used WWTP service populations ( Table S2 ) to weight wastewater measurements across WWTPs in each state and 2023 state populations to weight wastewater measurements across states for the United States ( Equations S3–S4 ). 74 Sewershed boundaries do not align with jurisdictional boundaries at which public health data are reported, so spatially aggregating wastewater measurements to state and national scales allowed us to compare wastewater measurements to clinical or vaccination data. For modeling (described later) and data visualization purposes only, we generated daily wastewater measurements for each WWTP prior to spatial aggregation. We used linear interpolation between adjacent wastewater measurements to generate a wastewater measurement each day that was either interpolated or based on measured values; no interpolation was performed if there was a gap in sampling > 10 days. Year comparison Rotavirus disease is known to follow a biennial occurrence pattern in the post-vaccination era, 4 , 17 and wastewater monitoring data were available in both 2023 and 2024 during fall (September–December). Log 10 -transformed wastewater measurements of rotavirus RNA were not normally distributed in 2023 and 2024 upon examination of quantile-quantile plots. Therefore, we used the two-sided Wilcoxon rank sum test to evaluate the null hypothesis that log 10 -transformed wastewater measurements of rotavirus RNA collected between September–December did not differ in 2023 compared to 2024. Wastewater events We inferred the onset and offset date of rotavirus epidemics for WWTPs using wastewater rotavirus RNA concentrations as measured (i.e., prior to smoothing and PMMoV normalization). Here, a wastewater event refers to the period between onset and offset dates, and the duration of a wastewater event is the number of days between corresponding onset and offset dates. We defined onset as the first date that wastewater rotavirus RNA concentrations in all wastewater samples collected in the past 14 days were ≥ 25,000 gc/g (approximately half the median concentration detected across all samples in this study). We defined offset as the first date after onset that wastewater rotavirus RNA concentrations in ≤ 50% of wastewater samples collected in the past 14 days were ≥ 25,000 gc/g. WWTPs did not have an onset or offset designation for the first 14 days of data availability. To identify wastewater events at the state scale, we defined onset as the first date that ≥ 50% of WWTPs in the state were in onset and offset as the first date after onset that < 50% of WWTPs in the state were in onset. 70 – 72 Likewise, we determined wastewater events at the national scale by defining onset as the first date that ≥ 50% of states in the United States were in onset and offset as the first date after onset that < 50% of states in the United States were in onset. 70 – 72 For all spatial scales (WWTP, state, national), we identified all occurrences of wastewater events. Sometimes more than one wastewater event occurred, so the longest wastewater event is referred to as the “wastewater event” below. Correlation with clinical metrics We evaluated the correlation between wastewater measurements of rotavirus RNA and clinical metrics of rotavirus infections. We assessed the correlation on a weekly basis at the national scale using two clinical rotavirus metrics: (1) rotavirus test positivity rate from NREVSS and (2) proportion of patients with a rotaviral enteritis encounter diagnosis from Epic Cosmos. Clinical metrics were obtained on a weekly basis; we averaged wastewater measurements across all samples for each WWTP on a weekly basis to match the reporting frequency of the clinical metrics prior to spatially aggregating measurements to the national scale. Nationally aggregated weekly wastewater measurements were not normally distributed (Shapiro-Wilk test, p < 0.0001); therefore, we used a two-sided Kendall’s tau correlation to test the null hypothesis that weekly wastewater measurements of rotavirus RNA are not correlated with weekly clinical metrics of rotavirus infections for the United States. We used the KendallTauB function in the DescTools library to determine the 95% confidence interval (CI) of Kendall’s tau estimates. 75 Because we conducted two hypothesis tests (i.e., one test for each clinical metric), we used an alpha value of 0.05 / 2 = 0.025 to account for multiple hypothesis testing when assessing statistical significance. Correlation with vaccination coverage We evaluated the correlation between the rotavirus wastewater event duration and rotavirus vaccination coverage at the state and sewershed scales. At the state scale, rotavirus vaccination coverage estimates were available from (1) NCIRD and (2) Epic Cosmos. Rotavirus vaccination coverage is not directly estimated at the sewershed scale, so we used the vaccination coverage estimate of the county that a sewershed predominantly occupies ( Table S2 ) to represent the vaccination coverage of sewersheds ( Figure S15 ). County-level rotavirus vaccination coverage estimates are only available from Epic Cosmos; therefore, we only used one estimate of rotavirus vaccination coverage to evaluate this correlation at the sewershed scale. For both spatial scales, we grouped states or sewersheds into vaccination coverage tertiles (lowest, middle, and highest levels of vaccination coverage). Wastewater event durations were not normally distributed among states or sewersheds within all vaccination tertiles (Shapiro-Wilk, p < 0.05); therefore, we used Kruskal-Wallis test to assess the null hypothesis that the wastewater event duration among states or sewersheds is equal across vaccination coverage tertiles. Because we conducted two hypothesis tests at the state scale (i.e., one test for each vaccination estimate), we used an adjusted p value of 0.05 / 2 = 0.025 to account for multiple hypothesis testing when assessing statistical significance. We used a p value of 0.05 to assess statistical significance at the sewershed scale because we only conducted one hypothesis test at this scale. When the null hypothesis was rejected, we used the kruskall_effsize function in the rstatix library to determine the eta-squared effect size measure and the Conover-Inman post-hoc test with a Bonferroni correction to determine significant pairwise differences. 76 , 77 Correlation with population characteristics We evaluated the correlation between rotavirus wastewater event duration and each population characteristic associated with rotavirus vaccination coverage at the sewershed scale. For each population characteristic except census region, we grouped sewersheds into tertiles (lowest, middle, and highest levels of the population characteristic); census region was already a categorical variable. Wastewater event durations were not normally distributed among sewersheds within all groupings for each population characteristic (Shapiro-Wilk, p < 0.05); therefore, we used Kruskal-Wallis test to assess the null hypothesis that the wastewater event duration among sewersheds is equal across groups for each population characteristic. We used a p value of 0.05 to assess statistical significance because we did not conduct exploratory multiple hypothesis testing. 14 When the null hypothesis was rejected, we used the kruskall_effsize function in the rstatix library to determine the eta-squared effect size measure and the Conover-Inman post-hoc test with a Bonferroni correction to determine significant pairwise differences. 76 , 77 Estimation of shedding population We estimated the fraction of the population shedding rotavirus RNA into wastewater from rotavirus RNA measurements in wastewater solids. We used a mass balance model established by Wolfe et al. 25 for SAR-CoV-2 RNA. The model is not virus specific and was adopted by Boehm et al. 78 for human norovirus GII RNA. Herein, we adopted the model for rotavirus RNA. Briefly, the model assumes fecal shedding is the predominant source of viral RNA into the waste stream and that the fraction of wastewater solids of fecal origin can be approximated from the concentration of PMMoV RNA in wastewater solids ( Equation S5 ). Similar to Boehm et al. 78 , we neglect viral RNA decay given the time sewage spends in the system, including the primary clarifier, is short (< 1 day) and that the first order decay rate constant for rotavirus RNA and PMMoV RNA can reasonably be assumed to be both small and similar in magnitude ( Equation S6 ). 25 , 79 – 82 Lastly, we assume the solid-liquid partitioning coefficient for rotavirus RNA and PMMoV RNA are equivalent like Boehm et al. 78 , leading to the simplified model shown in Equation 1 . F shed is the population fraction shedding rotavirus RNA in feces; C rota_ww and C PMMoV_ww are the concentration of rotavirus RNA and PMMoV RNA in wastewater solids, respectively; and C PMMoV_feces and C rota_feces are the concentration of PMMoV RNA and rotavirus RNA shed in feces, respectively. Arts et al. 26 provide a large dataset on PMMoV RNA shedding in feces which we used to characterize the distribution of C PMMoV_feces . Zheng et al. 27 compiled data on rotavirus RNA shedding in feces from infected individuals which we used to characterize the distribution of C rota_feces . Using these distributions in a Monte Carlo simulation, we modeled F shed as a function of the PMMoV-normalized concentration of rotavirus RNA in wastewater solids (C rota_ww /C PMMoV_ww ). We varied C rota_ww /C PMMoV_ww as log 10 (C rota_ww /C PMMoV_ww ) from -5 to -2 in 0.1 increments. For each C rota_ww /C PMMoV_ww , we conducted ten thousand simulations, randomly drawing C PMMoV_feces and C rota_feces from the distributions given by Arts et al. 26 and Zheng et al. 27 , to estimate the median and IQR of F shed . We also conducted the Monte Carlo simulation to model F shed each day using longitudinal measurements of C rota_ww /C PMMoV_ww from an individual WWTP and spatially aggregated across all WWTPs to represent the United States. We conducted a sensitivity analysis to evaluate the uncertainty of each model parameter. Following the method outlined by Xue et al. 83 , we determined the 25th, 50th, and 75th percentiles for each model input parameter (C rota_ww /C PMMoV_ww , C PMMoV_feces , C rota_feces ). For C rota_ww /C PMMoV_ww , we determined percentile values based on nationally aggregated wastewater measurements. Then, holding all other parameters at their 50th percentile value, we calculated F shed for each parameter twice: once at its 25th percentile value (F shed,25 ) and once at its 75th percentile value (F shed,75 ). We calculated the F shed,75 :F shed,25 ratio for each input parameter. A ratio less than 1 indicates F shed decreases as the input parameter increases, a ratio equal to 1 indicates F shed is not affected by input parameter increases, and a ratio greater than 1 indicates F shed increases as the input parameter increases (add other refs - see Sarah’s QMRA paper) All calculations, analyses, and visualizations were conducted in R (version 4.1.3). This study was reviewed by the Stanford University Committee for the Protection of Human Subjects and determined to be exempt from oversight. Author contributions Elana M. G. Chan: conceptualization, data curation, formal analysis, investigation, methodology, software, validation, visualization, writing - original draft, writing - review & editing; Alessandro Zulli: conceptualization, data curation, methodology, writing - review & editing; Alexandria B Boehm: conceptualization, data curation, funding acquisition, investigation, methodology, project administration, resources, supervision, validation, writing - original draft, writing - review & editing. Competing interests All authors declare no financial or non-financial competing interests. Data availability The datasets generated and analyzed for this study are available through the Stanford Digital Repository ( https://doi.org/10.25740/mv671wb2446 ). Code availability The underlying code for this study is available through the Stanford Digital Repository ( https://doi.org/10.25740/mv671wb2446 ). Supplementary information Rotavirus_SI.pdf Acknowledgments We thank the participating wastewater treatment plant staff for collecting samples for this project. Data used in this study came from Epic Cosmos, a dataset created in collaboration with a community of Epic health systems representing more than 284 million patient records from over 1,500 hospitals and 36,000 clinics from all 50 states, D.C., Lebanon, and Saudi Arabia. We thank Dr. Rebecca Linfield for assistance with understanding the Epic Cosmos dataset. This work was supported by a gift to ABB from the Sergey Brin Family Foundation. The funder played no role in study design, data collection, analysis and interpretation of data, or the writing of this manuscript. References 1. ↵ Du , Y. et al. Global burden and trends of rotavirus infection-associated deaths from 1990 to 2019: an observational trend study . Virol. J . 19 , 166 ( 2022 ). OpenUrl PubMed 2. ↵ Muhsen , K. & Omar , M . Chapter 114 - Rotavirus . in Molecular Medical Microbiology (Third Edition) (eds. Tang , Y.-W. et al. ) 2321 – 2338 ( Academic Press , 2024 ). doi: 10.1016/B978-0-12-818619-0.00052-6 . OpenUrl CrossRef 3. ↵ Cortese , M. M. & Haber , P. Chapter 19: Rotavirus . in Epidemiology and Prevention of Vaccine-Preventable Diseases ( 2024 ). 4. ↵ Baker , J. M. , Tate , J. E. , Gautam , R. & Parashar , U. D. Chapter 13: Rotavirus . in Manual for the Surveillance of Vaccine-Preventable Diseases ( 2025 ). 5. ↵ Cortese , M. M. , Parashar , U. D. & U.S. Centers for Disease Control and Prevention. Prevention of Rotavirus Gastroenteritis Among Infants and Children Recommendations of the Advisory Committee on Immunization Practices (ACIP) . MMWR Recomm. Rep . 58 , 1 – 25 ( 2009 ). OpenUrl PubMed 6. ↵ Burnett , E. , Parashar , U. D. & Tate , J. E. Real-world effectiveness of rotavirus vaccines, 2006–19: a literature review and meta-analysis . Lancet Glob. Health 8 , e1195 – e1202 ( 2020 ). OpenUrl 7. ↵ European Centre for Disease Prevention and Control (ECDC) . Expert Opinion on Rotavirus Vaccination in Infancy . doi 10.2900/362947 ( 2017 ). OpenUrl CrossRef 8. ↵ Anderson , E. J. & Weber , S. G . Rotavirus infection in adults . Lancet Infect. Dis . 4 , 91 – 99 ( 2004 ). OpenUrl CrossRef PubMed Web of Science 9. ↵ Arakaki , L. , Tollefson , D. , Kharono , B. & Drain , P. K . Prevalence of rotavirus among older children and adults with diarrhea: A systematic review and meta -analysis . Vaccine 39 , 4577 – 4590 ( 2021 ). OpenUrl PubMed 10. ↵ Crawford , S. E. et al. Rotavirus infection . Nat. Rev. Dis. Primer 3 , 1 – 16 ( 2017 ). OpenUrl CrossRef 11. ↵ Chavers , T. , Cates , J. , Burnett , E. , Parashar , U. D. & Tate , J. E . Indirect protection from rotavirus vaccines: a systematic review . Expert Rev. Vaccines 23 , 789 – 795 ( 2024 ). OpenUrl PubMed 12. ↵ Anderson , E. J. et al. Indirect Protection of Adults From Rotavirus by Pediatric Rotavirus Vaccination . Clin. Infect. Dis . 56 , 755 – 760 ( 2013 ). OpenUrl CrossRef PubMed 13. ↵ Hill , H. A. et al. Decline in Vaccination Coverage by Age 24 Months and Vaccination Inequities Among Children Born in 2020 and 2021 — National Immunization Survey-Child, United States, 2021–2023 . MMWR Morb. Mortal. Wkly. Rep . 73 , ( 2024 ). 14. ↵ Ghaswalla , P. K. , D’Angelo , J. & Abu-Elyazeed , R . Rotavirus vaccination in the US: a systematic review of vaccination coverage and completion . Hum. Vaccines Immunother . 17 , 872 – 879 ( 2021 ). OpenUrl 15. Valier , M. R . Vital Signs: Trends and Disparities in Childhood Vaccination Coverage by Vaccines for Children Program Eligibility — National Immunization Survey-Child, United States, 2012–2022 . MMWR Morb. Mortal. Wkly. Rep. 73 , ( 2024 ). 16. ↵ Rogers , M. A. M. , Kim , C. & Hofstetter , A. M . Geospatial Variation in Rotavirus Vaccination in Infants, United States, 2010–2017 . Emerg. Infect. Dis. 25 , 1993 – 1995 ( 2019 ). OpenUrl PubMed 17. ↵ Pindyck , T. , Tate , J. E. & Parashar , U. D . A decade of experience with rotavirus vaccination in the United States – vaccine uptake, effectiveness, and impact . Expert Rev. Vaccines 17 , 593 – 606 ( 2018 ). OpenUrl PubMed 18. ↵ Mukhopadhya , I. et al. Rotavirus Shedding in Symptomatic and Asymptomatic Children using Reverse Transcription-Quantitative PCR . J. Med. Virol . 85 , 1661 – 1668 ( 2013 ). OpenUrl CrossRef PubMed 19. ↵ Boehm , A. B. et al. A retrospective longitudinal study of adenovirus group F, norovirus GI and GII, rotavirus, and enterovirus nucleic acids in wastewater solids at two wastewater treatment plants: solid-liquid partitioning and relation to clinical testing data . mSphere 9 , e00736 – 23 ( 2024 ). OpenUrl CrossRef PubMed 20. Wolken , M. et al. PreK-12 school and citywide wastewater monitoring of the enteric viruses astrovirus, rotavirus, and sapovirus . Sci. Total Environ . 931 , 172683 ( 2024 ). OpenUrl PubMed 21. Zheng , X. et al. Tracking diarrhea viruses and mpox virus using the wastewater surveillance network in Hong Kong . Water Res . 255 , 121513 ( 2024 ). OpenUrl 22. Valdivia-Carrera , C. A. et al. Surveillance of SARS-CoV-2, rotavirus, norovirus genogroup II, and human adenovirus in wastewater as an epidemiological tool to anticipate outbreaks of COVID-19 and acute gastroenteritis in a city without a wastewater treatment plant in the Peruvian Highlands . Sci. Total Environ . 905 , 167161 ( 2023 ). OpenUrl PubMed 23. Wang , H. et al. Measures against COVID-19 affected the spread of human enteric viruses in a Swedish community, as found when monitoring wastewater . Sci. Total Environ . 895 , 165012 ( 2023 ). OpenUrl CrossRef PubMed 24. ↵ Ito , E. et al. Detection of Rotavirus Vaccine Strains in Oysters and Sewage and Their Relationship with the Gastroenteritis Epidemic . Appl. Environ. Microbiol . 87 , e02547 – 20 ( 2021 ). OpenUrl PubMed 25. ↵ Wolfe , M. K. et al. Scaling of SARS-CoV-2 RNA in Settled Solids from Multiple Wastewater Treatment Plants to Compare Incidence Rates of Laboratory-Confirmed COVID-19 in Their Sewersheds . Environ. Sci. Technol. Lett . 8 , 398 – 404 ( 2021 ). OpenUrl PubMed 26. ↵ Arts , P. J. et al. Longitudinal and quantitative fecal shedding dynamics of SARS-CoV-2, pepper mild mottle virus, and crAssphage . mSphere 8 , e00132 – 23 ( 2023 ). OpenUrl PubMed 27. ↵ Zheng , G. , Chan , E. M. G. & Boehm , A. Systematic review and meta-analysis of enteric virus shedding in human excretions . 2025.03.01.25323153 Preprint at doi: 10.1101/2025.03.01.25323153 ( 2025 ). OpenUrl Abstract / FREE Full Text 28. ↵ Gholipour , S. et al. Occurrence of viruses in sewage sludge: A systematic review . Sci. Total Environ . 824 , 153886 ( 2022 ). OpenUrl PubMed 29. ↵ Silva-Sales , M. , Martínez-Puchol , S. , Gonzales-Gustavson , E. , Hundesa , A. & Gironès , R . High Prevalence of Rotavirus A in Raw Sewage Samples from Northeast Spain . Viruses 12 , 318 ( 2020 ). OpenUrl PubMed 30. Zhou , N. et al. Continuous detection and genetic diversity of human rotavirus A in sewage in eastern China, 2013–2014 . Virol. J. 13 , 153 ( 2016 ). OpenUrl PubMed 31. Barril , P. A. et al. Rotavirus seasonality in urban sewage from Argentina: Effect of meteorological variables on the viral load and the genetic diversity . Environ. Res . 138 , 409 – 415 ( 2015 ). OpenUrl CrossRef 32. Bosch , A. , Pinto , R. M. , Blanch , A. R. & Jofre , J. T . Detection of human rotavirus in sewage through two concentration procedures . Water Res . 22 , 343 – 348 ( 1988 ). OpenUrl CrossRef 33. Hejkal , T. W. , Smith , E. M. & Gerba , C. P . Seasonal occurrence of rotavirus in sewage . Appl. Environ. Microbiol . 47 , 588 – 590 ( 1984 ). OpenUrl Abstract / FREE Full Text 34. ↵ Bucardo , F. , Lindgren , P.-E. , Svensson , L. & Nordgren , J . Low Prevalence of Rotavirus and High Prevalence of Norovirus in Hospital and Community Wastewater after Introduction of Rotavirus Vaccine in Nicaragua . PLOS ONE 6 , e25962 ( 2011 ). OpenUrl CrossRef PubMed 35. ↵ Lappe , B. L. et al. Predicting norovirus and rotavirus resurgence in the United States following the COVID-19 pandemic: a mathematical modelling study . BMC Infect. Dis . 23 , 254 ( 2023 ). OpenUrl PubMed 36. ↵ Turcios , R. M. et al. TEMPORAL AND GEOGRAPHIC TRENDS OF ROTAVIRUS ACTIVITY IN THE UNITED STATES, 1997–2004 . Pediatr. Infect. Dis. J. 25 , 451 ( 2006 ). OpenUrl CrossRef PubMed Web of Science 37. ↵ Török , T. J. et al. Visualizing geographic and temporal trends in rotavirus…L : The Pediatric Infectious Disease Journal . 38. ↵ Burnett , E. , Parashar , U. D. , Winn , A. , Curns , A. T. & Tate , J. E . Major Changes in Spatiotemporal Trends of US Rotavirus Laboratory Detections After Rotavirus Vaccine Introduction—2009–2021 . Pediatr. Infect. Dis. J . 41 , 759 ( 2022 ). OpenUrl PubMed 39. ↵ Epic Systems Corporation . About Epic Cosmos . https://cosmos.epic.com/about/ ( 2025 ). 40. ↵ Tarabichi , Y. , et al. The Cosmos Collaborative: A Vendor-Facilitated Electronic Health Record Data Aggregation Platform . ACI Open 5 , e36 – e46 ( 2021 ). OpenUrl CrossRef 41. ↵ Hallowell , B. D. , Parashar , U. D. , Curns , A. , DeGroote , N. P. & Tate , J. E. Trends in the Laboratory Detection of Rotavirus Before and After Implementation of Routine Rotavirus Vaccination — United States, 2000–2018 . MMWR Morb. Mortal. Wkly. Rep. 68 , ( 2019 ). 42. ↵ U.S. Census Bureau . Historical Population Change Data (1910-2020) . https://www.census.gov/data/tables/time-series/dec/popchange-data-text.html ( 2021 ). 43. ↵ Phillips , G. , Lopman , B. , Rodrigues , L. C. & Tam , C. C . Asymptomatic Rotavirus Infections in England: Prevalence, Characteristics, and Risk Factors . Am. J. Epidemiol . 171 , 1023 – 1030 ( 2010 ). OpenUrl CrossRef PubMed Web of Science 44. ↵ Barrón-Romero , B. L. , Barreda-González , J. , Doval-Ugalde , R. , Zermeño-Eguia Liz , J. & Huerta-Peña , M . Asymptomatic rotavirus infections in day care centers . J. Clin. Microbiol . 22 , 116 – 118 ( 1985 ). OpenUrl Abstract / FREE Full Text 45. ↵ Boehm , A. B. et al. Human pathogen nucleic acids in wastewater solids from 191 wastewater treatment plants in the United States . Sci. Data 11 , 1141 ( 2024 ). OpenUrl PubMed 46. ↵ Ye , Y. , Chang , P. H. , Hartert , J. & Wigginton , K. R . Reactivity of Enveloped Virus Genome, Proteins, and Lipids with Free Chlorine and UV254 . Environ. Sci. Technol . 52 , 7698 – 7708 ( 2018 ). OpenUrl CrossRef PubMed 47. ↵ Anderson , C. E. , Tong , J. , Zambrana , W. , Boehm , A. B. & Wolfe , M. K . Investigating the Efficacy of Various Handwashing Methods against Enveloped and Non-Enveloped Viruses . ( 2023 ) doi: 10.4269/ajtmh.22-0287 . OpenUrl CrossRef 48. ↵ Asghar , H. et al. Environmental surveillance for polioviruses in the Global Polio Eradication Initiative . J. Infect. Dis . 210 Suppl 1 , S294 – 303 ( 2014 ). OpenUrl CrossRef PubMed 49. ↵ Smith , T. , Cassell , G. & Bhatnagar , A . Wastewater Surveillance Can Have a Second Act in COVID-19 Vaccine Distribution . JAMA Health Forum 2 , e201616 ( 2021 ). OpenUrl 50. ↵ McCarthy , D. K. et al. Does wastewater surveillance have a role in vaccine-preventable disease control? Comparison of clinical and wastewater surveillance data for hepatitis A, E, measles, rubella and influenza viruses NICD, South Africa, 2021-2024 . Int. J. Infect. Dis. 152 , 107380 ( 2025 ). OpenUrl 51. ↵ Capone , D. , et al. Impact of sampling depth on pathogen detection in pit latrines . PLoS Negl. Trop. Dis. 15 , e0009176 ( 2021 ). OpenUrl CrossRef PubMed 52. Yang , Z . Low-cost and rapid sensors for wastewater surveillance at low-resource settings . Nat. Water 1 , 405 – 407 ( 2023 ). OpenUrl 53. Grijalva , L. M. et al. Persistence of respiratory, enteric, and fecal indicator viruses in fecal sludge from on-site sanitation in Dakar, Senegal . 2024.08.07.24302194 Preprint at doi: 10.1101/2024.08.07.24302194 ( 2024 ). OpenUrl Abstract / FREE Full Text 54. Bivins , A. et al. Passive sampling to scale wastewater surveillance of infectious disease: Lessons learned from COVID-19 . Sci. Total Environ . 835 , 155347 ( 2022 ). OpenUrl CrossRef PubMed 55. ↵ Hamilton , K. A. , Wade , M. J. , Barnes , K. G. , Street , R. A. & Paterson , S . Wastewater-based epidemiology as a public health resource in low- and middle-income settings . Environ. Pollut . 351 , 124045 ( 2024 ). OpenUrl CrossRef PubMed 56. ↵ Borchardt , M. A. et al. The Environmental Microbiology Minimum Information (EMMI) Guidelines: qPCR and dPCR Quality and Reporting for Environmental Microbiology . Environ. Sci. Technol . 55 , 10210 – 10223 ( 2021 ). OpenUrl CrossRef PubMed 57. ↵ Kitajima , M. , Sassi , H. P. & Torrey , J. R . Pepper mild mottle virus as a water quality indicator . Npj Clean Water 1 , 1 – 9 ( 2018 ). OpenUrl 58. ↵ Symonds , E. M. , Nguyen , K. H. , Harwood , V. J. & Breitbart , M . Pepper mild mottle virus: A plant pathogen with a greater purpose in (waste)water treatment development and public health management . Water Res . 144 , 1 – 12 ( 2018 ). OpenUrl CrossRef PubMed 59. ↵ Chan , E. , Zulli , A. & Boehm , A . Data and R code for ‘Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA’ . Stanford Digital Repository doi: 10.25740/mv671wb2446 ( 2025 ). OpenUrl CrossRef 60. ↵ U.S. Census Bureau. Cartographic Boundary Files . https://www.census.gov/geographies/mapping-files/time-series/geo/cartographic-boundary.html . 61. ↵ U.S. Centers for Disease Control and Prevention . The National Respiratory and Enteric Virus Surveillance System (NREVSS): Interactive Dashboard . https://app.powerbigov.us/view?r=eyJrIjoiMWQ4MGRkZDYtOWZhNS00YzZhLWE2MTYtMDNjMzY5MjQ4ZjYyIiwidCI6IjljZTcwODY5LTYwZGItNDRmZC1hYmU4LWQyNzY3MDc3ZmM4ZiJ9 ( 2024 ). 62. ↵ World Health Organization . International Statistical Classification of Diseases and Related Health Problems 10th Revision . https://icd.who.int/browse10/2019/en ( 2019 ). 63. ↵ U.S. Centers for Disease Control and Prevention. About the National Immunization Surveys (NIS) . https://www.cdc.gov/nis/about/index.html ( 2025 ). 64. ↵ National Center for Immunization and Respiratory Diseases. Vaccination Coverage among Young Children (0-35 Months) . https://data.cdc.gov/Child-Vaccinations/Vaccination-Coverage-among-Young-Children-0-35-Mon/fhky-rtsk/about_data ( 2024 ). 65. ↵ U.S. Census Bureau . 2018-2022 American Community Survey 5-Year Data . https://www.census.gov/programs-surveys/acs/data.html . 66. The American Board of Pediatrics . 2024 Pediatric State and County Distributions Data . Chapel Hill, NC : The American Board of Pediatrics https://www.abp.org/dashboards/general-pediatricians-us-state-and-county-maps ( 2024 ). 67. U.S. Census Bureau . 2020 Census Urban-Rural Classification . https://www.census.gov/programs-surveys/geography/guidance/geo-areas/urban-rural.html ( 2024 ). 68. ↵ U.S. Census Bureau . Census Regions and Divisions of the United States . ( 1984 ). 69. ↵ Esri Data and Maps. USA ZIP Code Boundaries . https://www.arcgis.com/home/item.html?id=d6f7ee6129e241cc9b6f75978e47128b ( 2023 ). 70. ↵ Chan , E. M. G. & Boehm , A. B. Respiratory Virus Season Surveillance in the United States Using Wastewater Metrics, 2023–2024 . ACS EST Water 5 , 985 – 992 ( 2025 ). OpenUrl 71. Zulli , A. , Varkila , M. R. J. , Parsonnet , J. , Wolfe , M. K. & Boehm , A. B . Observations of Respiratory Syncytial Virus (RSV) Nucleic Acids in Wastewater Solids Across the United States in the 2022–2023 Season: Relationships with RSV Infection Positivity and Hospitalization Rates . ACS EST Water 4 , 1657 – 1667 ( 2024 ). OpenUrl 72. ↵ Schoen , M. E. , Bidwell , A. L. , Wolfe , M. K. & Boehm , A. B. United States Influenza 2022– 2023 Season Characteristics as Inferred from Wastewater Solids, Influenza Hospitalization, and Syndromic Data . Environ. Sci. Technol . 57 , 20542 – 20550 ( 2023 ). OpenUrl CrossRef PubMed 73. ↵ Schenk , H. , Rauch , W. , Zulli , A. & Boehm , A. B . SARS-CoV-2 surveillance in US wastewater: Leading indicators and data variability analysis in 2023–2024 . PLOS ONE 19 , e0313927 ( 2024 ). OpenUrl CrossRef PubMed 74. ↵ U.S. Census Bureau, Population Division. Annual Estimates of the Resident Population for the United States, Regions, States, District of Columbia, and Puerto Rico: April 1, 2020 to July 1, 2023 (NST-EST2023-POP) . ( 2023 ). 75. ↵ Signorell , A. et al. DescTools: Tools for Descriptive Statistics . ( 2023 ). 76. ↵ Kassambara , A. rstatix: Pipe-Friendly Framework for Basic Statistical Tests . ( 2023 ). 77. ↵ Dinno , A. conover.test: Conover-Iman Test of Multiple Comparisons Using Rank Sums . ( 2024 ). 78. ↵ Boehm , A. B. et al. Human norovirus (HuNoV) GII RNA in wastewater solids at 145 United States wastewater treatment plants: comparison to positivity rates of clinical specimens and modeled estimates of HuNoV GII shedders . J. Expo. Sci. Environ. Epidemiol . 34 , 440 – 447 ( 2024 ). OpenUrl PubMed 79. ↵ Kraay , A. N. M. , et al. Modeling environmentally mediated rotavirus transmission: The role of temperature and hydrologic factors . Proc. Natl. Acad. Sci. 115 , E2782 – E2790 ( 2018 ). OpenUrl Abstract / FREE Full Text 80. Boehm , A. B. , Silverman , A. I. , Schriewer , A. & Goodwin , K . Systematic review and meta-analysis of decay rates of waterborne mammalian viruses and coliphages in surface waters . Water Res . 164 , 114898 ( 2019 ). OpenUrl CrossRef 81. Roldan-Hernandez , L. , Graham , K. E. , Duong , D. & Boehm , A. B . Persistence of Endogenous SARS-CoV-2 and Pepper Mild Mottle Virus RNA in Wastewater-Settled Solids . ACS EST Water 2 , 1944 – 1952 ( 2022 ). OpenUrl 82. ↵ Burnet , J.-B. , Cauchie , H.-M. , Walczak , C. , Goeders , N. & Ogorzaly , L . Persistence of endogenous RNA biomarkers of SARS-CoV-2 and PMMoV in raw wastewater: Impact of temperature and implications for wastewater-based epidemiology . Sci. Total Environ . 857 , 159401 ( 2023 ). OpenUrl CrossRef PubMed 83. ↵ Xue , J. et al. A Probabilistic Arsenic Exposure Assessment for Children Who Contact Chromated Copper Arsenate (CCA)-Treated Playsets and Decks, Part 2: Sensitivity and Uncertainty Analyses . Risk Anal. 26 , 533 – 541 ( 2006 ). OpenUrl PubMed View the discussion thread. Back to top Previous Next Posted July 22, 2025. Download PDF Supplementary Material Data/Code Email Thank you for your interest in spreading the word about medRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA Message Subject (Your Name) has forwarded a page to you from medRxiv Message Body (Your Name) thought you would like to see this page from the medRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA Elana M. G. Chan , Alessandro Zulli , Alexandria B. Boehm medRxiv 2025.07.21.25331943; doi: https://doi.org/10.1101/2025.07.21.25331943 Share This Article: Copy Citation Tools Wastewater concentrations of rotavirus RNA are associated with infection and vaccination metrics in the USA Elana M. G. Chan , Alessandro Zulli , Alexandria B. Boehm medRxiv 2025.07.21.25331943; doi: https://doi.org/10.1101/2025.07.21.25331943 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Epidemiology Subject Areas All Articles Addiction Medicine (568) Allergy and Immunology (863) Anesthesia (297) Cardiovascular Medicine (4421) Dentistry and Oral Medicine (443) Dermatology (381) Emergency Medicine (606) Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1507) Epidemiology (15212) Forensic Medicine (30) Gastroenterology (1121) Genetic and Genomic Medicine (6581) Geriatric Medicine (667) Health Economics (996) Health Informatics (4520) Health Policy (1366) Health Systems and Quality Improvement (1611) Hematology (539) HIV/AIDS (1264) Infectious Diseases (except HIV/AIDS) (15906) Intensive Care and Critical Care Medicine (1103) Medical Education (620) Medical Ethics (144) Nephrology (667) Neurology (6580) Nursing (345) Nutrition (998) Obstetrics and Gynecology (1141) Occupational and Environmental Health (956) Oncology (3324) Ophthalmology (970) Orthopedics (369) Otolaryngology (420) Pain Medicine (435) Palliative Medicine (129) Pathology (663) Pediatrics (1689) Pharmacology and Therapeutics (691) Primary Care Research (710) Psychiatry and Clinical Psychology (5431) Public and Global Health (9212) Radiology and Imaging (2193) Rehabilitation Medicine and Physical Therapy (1368) Respiratory Medicine (1194) Rheumatology (593) Sexual and Reproductive Health (709) Sports Medicine (529) Surgery (709) Toxicology (99) Transplantation (288) Urology (265) (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'9ff0890d58764807',t:'MTc3OTMzNDU3MQ=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();