Visualizing Neighborhood COVID-19 Levels, Trends, and Inequities in Wastewater: An Equity-Centered Approach and Comparison to CDC Methods.

CONTEXT: Monitoring neighborhood-level SARS-CoV-2 wastewater concentrations can help guide public health interventions and provide early warning ahead of lagging COVID-19 clinical indicators. To date, however, U.S. Centers for Disease Control and Prevention's (CDC) National Wastewater Surveillance System (NWSS) has provided methodology solely for communicating national and state-level "wastewater viral activity levels." PROGRAM: In October 2022, the Boston Public Health Commission (BPHC) began routinely sampling wastewater at 11 neighborhood sites to better understand COVID-19 epidemiology and inequities across neighborhoods, which vary widely in sociodemographic and socioeconomic characteristics. We developed equity-centered methods to routinely report interpretable and actionable descriptions of COVID-19 wastewater levels, trends, and neighborhood-level inequities. APPROACH AND IMPLEMENTATION: To produce these data visualizations, spanning October 2022 to December 2023, we followed four general steps: (1) smoothing raw values; (2) classifying current COVID-19 wastewater levels; (3) classifying current trends; and (4) reporting and visualizing results. EVALUATION: COVID-19 wastewater levels corresponded well with lagged COVID-19 hospitalizations and deaths over time, with "Very High" wastewater levels coinciding with winter surges. When citywide COVID-19 levels were at the highest and lowest points, levels and trends tended to be consistent across sites. In contrast, when citywide levels were moderate, neighborhood levels and trends were more variable, revealing inequities across neighborhoods, emphasizing the importance of neighborhood-level results. Applying CDC/NWSS state-level methodology to neighborhood sites resulted in vastly different neighborhood-specific wastewater cut points for "High" or "Low," obscured inequities between neighborhoods, and systematically underestimated COVID-19 levels during surge periods in neighborhoods with the highest COVID-19 morbidity and mortality. DISCUSSION: Our methods offer an approach that other local jurisdictions can use for routinely monitoring, comparing, and communicating neighborhood-level wastewater levels, trends, and inequities. Applying CDC/NWSS methodology at the neighborhood-level can obscure and perpetuate COVID-19 inequities. We recommend jurisdictions adopt equity-focused approaches in neighborhood-level wastewater surveillance for valid community comparisons.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting.

After introducing pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in Streptococcus pneumoniae. Predicting which pneumococcal strains will become common in carriage after vaccination can enhance vaccine design, public health interventions, and understanding of pneumococcal evolution. Invasive pneumococcal isolates were collected during 1998-2018 by the Active Bacterial Core surveillance (ABCs). Carriage data from Massachusetts (MA) and Southwest United States were used to calculate weights. Using pre-vaccine data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in carriage to the proportion in invasive data. Genomic data were processed under bioinformatic pipelines to define genetically similar sequence clusters (i.e., strains), and accessory genes (COGs) present in 5-95% of isolates. Weights were applied to adjust observed strain proportions and COG frequencies. The negative frequency-dependent selection (NFDS) model predicted strain proportions by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine COG frequencies. Inverse-invasiveness weighting increased the correlation of COG frequencies between invasive and carriage data in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13; and between different epochs in the invasive data. Weighting the invasive data significantly improved the NFDS model's accuracy in predicting strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R2 increasing from 0.254 before weighting to 0.545 after weighting. The weighting system adjusted invasive disease data to better represent the pneumococcal carriage population, allowing the NFDS mechanism to predict strain proportions in carriage in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance.IMPORTANCEStreptococcus pneumoniae, a common colonizer in the human nasopharynx, can cause invasive diseases including pneumonia, bacteremia, and meningitis mostly in children under 5 years or older adults. The PCV7 was introduced in 2000 in the United States within the pediatric population to prevent disease and reduce deaths, followed by PCV13 in 2010, PCV15 in 2022, and PCV20 in 2023. After the removal of vaccine serotypes, the prevalence of carriage remained stable as the vacated pediatric ecological niche was filled with certain non-vaccine serotypes. Predicting which pneumococcal clones, and which serotypes, will be most successful in colonization after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. While carriage data, which are collected from the pneumococcal population that is competing to colonize and transmit, are most directly relevant to evolutionary studies, invasive disease data are often more plentiful. Previously, evolutionary models based on negative frequency-dependent selection (NFDS) on the accessory genome were shown to predict which non-vaccine strains and serotypes were most successful in colonization following the introduction of PCV7. Here, we show that an inverse-invasiveness weighting system applied to invasive disease surveillance data allows the NFDS model to predict strain proportions in the projected carriage population in the post-PCV13/pre-PCV15 and pre-PCV20 epoch. The significance of our research lies in using a sample of invasive disease surveillance data to extend the use of NFDS as an evolutionary mechanism to predict post-PCV13 population dynamics. This has shown that we can correct for biased sampling that arises from differences in virulence and can enrich the value of genomic data from disease surveillance and advance our understanding of how NFDS impacts carriage population dynamics after both PCV7 and PCV13 vaccination.

Successful treatment of SARS-CoV-2 in an immunocompromised patient with persistent infection for 245 days: A case report.

Background: Immunocompromised patients receiving B-cell-depleting therapies are at increased risk of persistent SARS-CoV-2 infection, with many experiencing fatal outcomes. We report a successful outcome in a patient with rheumatoid arthritis (RA) on rituximab diagnosed with COVID-19 in July 2020 with persistent infection for over 245 days. Results: The patient received numerous treatment courses for persistent COVID-19 infection, including remdesivir, baricitinib, immunoglobulin and high doses of corticosteroids followed by a prolonged taper due to persistent respiratory symptoms and cryptogenic organizing pneumonia. Her clinical course was complicated by Pseudomonas aeruginosa sinusitis with secondary bacteremia, and cytomegalovirus (CMV) viremia and pneumonitis. SARS-CoV-2 positive RNA samples were extracted from two nasopharyngeal swabs and sequenced using targeted amplicon Next-Generation Sequencing which were analyzed for virus evolution over time. Viral sequencing indicated lineage B.1.585.3 SARS-CoV-2 accumulated Spike protein mutations associated with immune evasion and resistance to therapeutics. Upon slowly decreasing the patient's steroids, she had resolution of her symptoms and had a negative nasopharyngeal SARS-CoV-2 PCR and serum CMV PCR in March 2021. Conclusion: A patient with RA on B-cell depleting therapy developed persistent SARS-CoV-2 infection allowing for virus evolution and had numerous complications, including viral and bacterial co-infections with opportunistic pathogens. Despite intra-host evolution with a more immune evasive SARS-CoV-2 lineage, it was cleared after 245 days with reconstitution of the patient's immune system.

Sink-traps are a major source for carbapenemase-producing Enterobacteriaceae transmission.

OBJECTIVE: We studied the extent of carbapenemase-producing Enterobacteriaceae (CPE) sink contamination and transmission to patients in a nonoutbreak setting. METHODS: During 2017-2019, 592 patient-room sinks were sampled in 34 departments. Patient weekly rectal swab CPE surveillance was universally performed. Repeated sink sampling was conducted in 9 departments. Isolates from patients and sinks were characterized using pulsed-field gel electrophoresis (PFGE), and pairs of high resemblance were sequenced by Oxford Nanopore and Illumina. Hybrid assembly was used to fully assemble plasmids, which are shared between paired isolates. RESULTS: In total, 144 (24%) of 592 CPE-contaminated sinks were detected in 25 of 34 departments. Repeated sampling (n = 7,123) revealed that 52%-100% were contaminated at least once during the sampling period. Persistent contamination for >1 year by a dominant strain was common. During the study period, 318 patients acquired CPE. The most common species were Klebsiella pneumoniae, Escherichia coli, and Enterobacter spp. In 127 (40%) patients, a contaminated sink was the suspected source of CPE acquisition. For 20 cases with an identical sink-patient strain, temporal relation suggested sink-to-patient transmission. Hybrid assembly of specific sink-patient isolates revealed that shared plasmids were structurally identical, and SNP differences between shared pairs, along with signatures for potential recombination events, suggests recent sharing of the plasmids. CONCLUSIONS: CPE-contaminated sinks are an important source of transmission to patients. Although traditionally person-to-person transmission has been considered the main route of CPE transmission, these data suggest a change in paradigm that may influence strategies of preventing CPE dissemination.

Quantitatively assessing early detection strategies for mitigating COVID-19 and future pandemics.

Researchers and policymakers have proposed systems to detect novel pathogens earlier than existing surveillance systems by monitoring samples from hospital patients, wastewater, and air travel, in order to mitigate future pandemics. How much benefit would such systems offer? We developed, empirically validated, and mathematically characterized a quantitative model that simulates disease spread and detection time for any given disease and detection system. We find that hospital monitoring could have detected COVID-19 in Wuhan 0.4 weeks earlier than it was actually discovered, at 2,300 cases (standard error: 76 cases) compared to 3,400 (standard error: 161 cases). Wastewater monitoring would not have accelerated COVID-19 detection in Wuhan, but provides benefit in smaller catchments and for asymptomatic or long-incubation diseases like polio or HIV/AIDS. Air travel monitoring does not accelerate outbreak detection in most scenarios we evaluated. In sum, early detection systems can substantially mitigate some future pandemics, but would not have changed the course of COVID-19.

Inferring bacterial transmission dynamics using deep sequencing genomic surveillance data.

Identifying and interrupting transmission chains is important for controlling infectious diseases. One way to identify transmission pairs - two hosts in which infection was transmitted from one to the other - is using the variation of the pathogen within each single host (within-host variation). However, the role of such variation in transmission is understudied due to a lack of experimental and clinical datasets that capture pathogen diversity in both donor and recipient hosts. In this work, we assess the utility of deep-sequenced genomic surveillance (where genomic regions are sequenced hundreds to thousands of times) using a mouse transmission model involving controlled spread of the pathogenic bacterium Citrobacter rodentium from infected to naïve female animals. We observe that within-host single nucleotide variants (iSNVs) are maintained over multiple transmission steps and present a model for inferring the likelihood that a given pair of sequenced samples are linked by transmission. In this work we show that, beyond the presence and absence of within-host variants, differences arising in the relative abundance of iSNVs (allelic frequency) can infer transmission pairs more precisely. Our approach further highlights the critical role bottlenecks play in reserving the within-host diversity during transmission.

Quantitatively assessing early detection strategies for mitigating COVID-19 and future pandemics.

Researchers and policymakers have proposed systems to detect novel pathogens earlier than existing surveillance systems by monitoring samples from hospital patients, wastewater, and air travel, in order to mitigate future pandemics. How much benefit would such systems offer? We developed, empirically validated, and mathematically characterized a quantitative model that simulates disease spread and detection time for any given disease and detection system. We find that hospital monitoring could have detected COVID-19 in Wuhan 0.4 weeks earlier than it was actually discovered, at 2,300 cases (standard error: 76 cases) compared to 3,400 (standard error: 161 cases). Wastewater monitoring would not have accelerated COVID-19 detection in Wuhan, but provides benefit in smaller catchments and for asymptomatic or long-incubation diseases like polio or HIV/AIDS. Monitoring of air travel provides little benefit in most scenarios we evaluated. In sum, early detection systems can substantially mitigate some future pandemics, but would not have changed the course of COVID-19.

Estimated preventable COVID-19-associated deaths due to non-vaccination in the United States.

While some studies have previously estimated lives saved by COVID-19 vaccination, we estimate how many deaths could have been averted by vaccination in the US but were not because of a failure to vaccinate. We used a simple method based on a nationally representative dataset to estimate the preventable deaths among unvaccinated individuals in the US from May 30, 2021 to September 3, 2022 adjusted for the effects of age and time. We estimated that at least 232,000 deaths could have been prevented among unvaccinated adults during the 15 months had they been vaccinated with at least a primary series. While uncertainties exist regarding the exact number of preventable deaths and more granular data are needed on other factors causing differences in death rates between the vaccinated and unvaccinated groups to inform these estimates, this method is a rapid assessment on vaccine-preventable deaths due to SARS-CoV-2 that has crucial public health implications. The same rapid method can be used for future public health emergencies.

Addressing spatial misalignment in population health research: a case study of US congressional district political metrics and county health data.

Areal spatial misalignment, which occurs when data on multiple variables are collected using mismatched boundary definitions, is a ubiquitous obstacle to data analysis in public health and social science research. As one example, the emerging sub-field studying the links between political context and health in the United States faces significant spatial misalignment-related challenges, as the congressional districts (CDs) over which political metrics are measured and administrative units, e.g., counties, for which health data are typically released, have a complex misalignment structure. Standard population-weighted data realignment procedures can induce measurement error and invalidate inference, which has prompted the development of fully model-based approaches for analyzing spatially misaligned data. One such approach, atom-based regression models (ABRM), holds particular promise but has scarcely been used in practice due to the lack of appropriate software or examples of implementation. ABRM use "atoms", the areas created by intersecting all sets of units on which variables of interest are measured, as the units of analysis and build models for the atom-level data, treating the atom-level variables (generally unmeasured) as latent variables. In this paper, we demonstrate the feasibility and strengths of the ABRM in a case study of the association between political representatives' voting behavior (CD-level) and COVID-19 mortality rates (county-level) in a post-vaccine period. The adjusted ABRM results suggest that more conservative voting record is associated with an increase in COVID-19 mortality rates, with estimated associations smaller in magnitude but consistent in direction with those of standard realignment methods. The results also indicate that ABRM may enable more robust confounding adjustment and more realistic uncertainty estimates, properly representing the uncertainties arising from all analytic procedures. We also implement the ABRM in modern optimized Bayesian computing programs and make our code publicly available, which may enable these methods to be more widely adopted.

Early Introduction and Rise of the Omicron Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variant in Highly Vaccinated University Populations.

BACKGROUND: The Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly transmissible in vaccinated and unvaccinated populations. The dynamics that govern its establishment and propensity toward fixation (reaching 100% frequency in the SARS-CoV-2 population) in communities remain unknown. Here, we describe the dynamics of Omicron at 3 institutions of higher education (IHEs) in the greater Boston area. METHODS: We use diagnostic and variant-specifying molecular assays and epidemiological analytical approaches to describe the rapid dominance of Omicron following its introduction into 3 IHEs with asymptomatic surveillance programs. RESULTS: We show that the establishment of Omicron at IHEs precedes that of the state and region and that the time to fixation is shorter at IHEs (9.5-12.5 days) than in the state (14.8 days) or region. We show that the trajectory of Omicron fixation among university employees resembles that of students, with a 2- to 3-day delay. Finally, we compare cycle threshold values in Omicron vs Delta variant cases on college campuses and identify lower viral loads among college affiliates who harbor Omicron infections. CONCLUSIONS: We document the rapid takeover of the Omicron variant at IHEs, reaching near-fixation within the span of 9.5-12.5 days despite lower viral loads, on average, than the previously dominant Delta variant. These findings highlight the transmissibility of Omicron, its propensity to rapidly dominate small populations, and the ability of robust asymptomatic surveillance programs to offer early insights into the dynamics of pathogen arrival and spread.

Prediction of post-PCV13 pneumococcal evolution using invasive disease data enhanced by inverse-invasiveness weighting.

Background: After introduction of pneumococcal conjugate vaccines (PCVs), serotype replacement occurred in the population of Streptococcus pneumoniae. Predicting which pneumococcal clones and serotypes will become more common in carriage after vaccination can enhance vaccine design and public health interventions, while also improving our understanding of pneumococcal evolution. We sought to use invasive disease data to assess how well negative frequency-dependent selection (NFDS) models could explain pneumococcal carriage population evolution in the post-PCV13 epoch by weighting invasive data to approximate strain proportions in the carriage population. Methods: Invasive pneumococcal isolates were collected and sequenced during 1998-2018 by the Active Bacterial Core surveillance (ABCs) from the Centers for Disease Control and Prevention (CDC). To predict the post-PCV13 population dynamics in the carriage population using a NFDS model, all genomic data were processed under a bioinformatic pipeline of assembly, annotation, and pangenome analysis to define genetically similar sequence clusters (i.e., strains) and a set of accessory genes present in 5% to 95% of the isolates. The NFDS model predicted the strain proportion by calculating the post-vaccine strain composition in the weighted invasive disease population that would best match pre-vaccine accessory gene frequencies. To overcome the biases of invasive disease data, serotype-specific inverse-invasiveness weights were defined as the ratio of the proportion of the serotype in the carriage data to the proportion in the invasive data, using data from 1998-2001 in the United States, before conjugate vaccine introduction. The weights were applied to adjust both the observed strain proportion and the accessory gene frequencies. Results: Inverse-invasiveness weighting increased the correlation of accessory gene frequencies between invasive and carriage data with reduced residuals in linear or logit scale for pre-vaccine, post-PCV7, and post-PCV13. Similarly, weighting increased the correlation of accessory gene frequencies between different time periods in the invasive data. By weighting the invasive data, we were able to use the NFDS model to predict strain proportions in the carriage population in the post-PCV13 epoch, with the adjusted R-squared between predicted and observed strain proportions increasing from 0.176 to 0.544 after weighting. Conclusions: The weighting system adjusted the invasive disease surveillance data to better represent the carriage population of S. pneumoniae. The NFDS mechanism predicted the strain proportions in the projected carriage population as estimated from the weighted invasive disease frequencies in the post-PCV13 epoch. Our methods enrich the value of genomic sequences from invasive disease surveillance, which is readily available, easy to collect, and of direct interest to public health.

Monkeypox Virus Infections in Humans.

Human monkeypox is a viral zoonosis endemic to West and Central Africa that has recently generated increased interest and concern on a global scale as an emerging infectious disease threat in the midst of the slowly relenting COVID-2019 disease pandemic. The hallmark of infection is the development of a flu-like prodrome followed by the appearance of a smallpox-like exanthem. Precipitous person-to-person transmission of the virus among residents of 100 countries where it is nonendemic has motivated the immediate and widespread implementation of public health countermeasures. In this review, we discuss the origins and virology of monkeypox virus, its link with smallpox eradication, its record of causing outbreaks of human disease in regions where it is endemic in wildlife, its association with outbreaks in areas where it is nonendemic, the clinical manifestations of disease, laboratory diagnostic methods, case management, public health interventions, and future directions.

Relationship of political ideology of US federal and state elected officials and key COVID pandemic outcomes following vaccine rollout to adults: April 2021-March 2022.

Background: Scant research, including in the United States, has quantified relationships between the political ideologies of elected representatives and COVID-19 outcomes among their constituents. Methods: We analyzed observational cross-sectional data on COVID-19 mortality rates (age-standardized) and stress on hospital intensive care unit (ICU) capacity for all 435 US Congressional Districts (CDs) in a period of adult vaccine availability (April 2021-March 2022). Political metrics comprised: (1) ideological scores based on each US Representative's and Senator's concurrent overall voting record and their specific COVID-19 votes, and (2) state trifectas (Governor, State House, and State Senate under the same political party control). Analyses controlled for CD social metrics, population density, vaccination rates, the prevalence of diabetes and obesity, and voter political lean. Findings: During the study period, the higher the exposure to conservatism across several political metrics, the higher the COVID-19 age-standardized mortality rates, even after taking into account the CD's social characteristics; similar patterns occurred for stress on hospital ICU capacity for Republican trifectas and US Senator political ideology scores. For example, in models mutually adjusting for CD political and social metrics and vaccination rates, Republican trifecta and conservative voter political lean independently remained significantly associated with an 11%-26% higher COVID-19 mortality rate. Interpretation: Associations between the political ideologies of US federal elected officials and state concentrations of political party power with population health warrant greater consideration in public health analyses and monitoring dashboards. Funding: This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Linking contact tracing with genomic surveillance to deconvolute SARS-CoV-2 transmission on a university campus.

Contact tracing and genomic data, approaches often used separately, have both been important tools in understanding the nature of SARS-CoV-2 transmission. Linked analysis of contact tracing and sequence relatedness of SARS-CoV-2 genomes from a regularly sampled university environment were used to build a multilevel transmission tracing and confirmation system to monitor and understand transmission on campus. Our investigation of an 18-person cluster stemming from an athletic team highlighted the importance of linking contact tracing and genomic analysis. Through these findings, it is suggestive that certain safety protocols in the athletic practice setting reduced transmission. The linking of traditional contact tracing with rapid-return genomic information is an effective approach for differentiating between multiple plausible transmission scenarios and informing subsequent public health protocols to limit disease spread in a university environment.

Causal discovery for the microbiome.

Measurement and manipulation of the microbiome is generally considered to have great potential for understanding the causes of complex diseases in humans, developing new therapies, and finding preventive measures. Many studies have found significant associations between the microbiome and various diseases; however, Koch's classical postulates remind us about the importance of causative reasoning when considering the relationship between microbes and a disease manifestation. Although causal discovery in observational microbiome data faces many challenges, methodological advances in causal structure learning have improved the potential of data-driven prediction of causal effects in large-scale biological systems. In this Personal View, we show the capability of existing methods for inferring causal effects from metagenomic data, and we highlight ways in which the introduction of causal structures that are more flexible than existing structures offers new opportunities for causal reasoning. Our observations suggest that microbiome research can further benefit from tools developed in the past 5 years in causal discovery and learn from their applications elsewhere.

Understanding Early Pandemic Severe Acute Respiratory Syndrome Coronavirus 2 Transmission in a Medical Center by Incorporating Public Sequencing Databases to Mitigate Bias.

BACKGROUND: Throughout the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, healthcare workers (HCWs) have faced risk of infection from within the workplace via patients and staff as well as from the outside community, complicating our ability to resolve transmission chains in order to inform hospital infection control policy. Here we show how the incorporation of sequences from public genomic databases aided genomic surveillance early in the pandemic when circulating viral diversity was limited. METHODS: We sequenced a subset of discarded, diagnostic SARS-CoV-2 isolates between March and May 2020 from Boston Medical Center HCWs and combined this data set with publicly available sequences from the surrounding community deposited in GISAID with the goal of inferring specific transmission routes. RESULTS: Contextualizing our data with publicly available sequences reveals that 73% (95% confidence interval, 63%-84%) of coronavirus disease 2019 cases in HCWs are likely novel introductions rather than nosocomial spread. CONCLUSIONS: We argue that introductions of SARS-CoV-2 into the hospital environment are frequent and that expanding public genomic surveillance can better aid infection control when determining routes of transmission.

Analyzing the impact of a real-life outbreak simulator on pandemic mitigation: An epidemiological modeling study.

An app-based educational outbreak simulator, Operation Outbreak (OO), seeks to engage and educate participants to better respond to outbreaks. Here, we examine the utility of OO for understanding epidemiological dynamics. The OO app enables experience-based learning about outbreaks, spreading a virtual pathogen via Bluetooth among participating smartphones. Deployed at many colleges and in other settings, OO collects anonymized spatiotemporal data, including the time and duration of the contacts among participants of the simulation. We report the distribution, timing, duration, and connectedness of student social contacts at two university deployments and uncover cryptic transmission pathways through individuals' second-degree contacts. We then construct epidemiological models based on the OO-generated contact networks to predict the transmission pathways of hypothetical pathogens with varying reproductive numbers. Finally, we demonstrate that the granularity of OO data enables institutions to mitigate outbreaks by proactively and strategically testing and/or vaccinating individuals based on individual social interaction levels.