Local-Scale phylodynamics reveal differential community impact of SARS-CoV-2 in metropolitan US county

SARS-CoV-2 transmission is largely driven by heterogeneous dynamics at a local scale, leaving local health departments to design interventions with limited information. We analyzed SARS-CoV-2 genomes sampled between February 2020 and March 2022 jointly with epidemiological and cell phone mobility data to investigate fine scale spatiotemporal SARS-CoV-2 transmission dynamics in King County, Washington, a diverse, metropolitan US county. We applied an approximate structured coalescent approach to model transmission within and between North King County and South King County alongside the rate of outside introductions into the county. Our phylodynamic analyses reveal that following stay-at-home orders, the epidemic trajectories of North and South King County began to diverge. We find that South King County consistently had more reported and estimated cases, COVID-19 hospitalizations, and longer persistence of local viral transmission when compared to North King County, where viral importations from outside drove a larger proportion of new cases. Using mobility and demographic data, we also find that South King County experienced a more modest and less sustained reduction in mobility following stay-at-home orders than North King County, while also bearing more socioeconomic inequities that might contribute to a disproportionate burden of SARS-CoV-2 transmission. Overall, our findings suggest a role for local-scale phylodynamics in understanding the heterogeneous transmission landscape. One Sentence SummaryAnalysis of SARS-CoV-2 genomes in King County, Washington show that diverse areas in the same metropolitan region can have different epidemic dynamics.

Trends in risk factors and symptoms associated with SARS-CoV-2 and Rhinovirus test positivity in King County, Washington: A Test-Negative Design Study of the Greater Seattle Coronavirus Assessment Network

ImportanceFew US studies have reexamined risk factors for SARS-CoV-2 positivity in the context of widespread vaccination and new variants or considered risk factors for co-circulating endemic viruses, such as rhinovirus. ObjectiveTo understand how risk factors and symptoms associated with SARS-CoV-2 test positivity changed over the course of the pandemic and to compare these to the factors associated with rhinovirus test positivity. DesignThis test-negative design study used multivariable logistic regression to assess associations between SARS-CoV-2 and rhinovirus test positivity and self-reported demographic and symptom variables over a 22-month period. SettingKing County, Washington, June 2020-April 2022 Participants23,278 symptomatic individuals of all ages enrolled in a cross-sectional community surveillance study. ExposuresSelf-reported data for 15 demographic and health behavior variables and 16 symptoms. Main Outcome(s) and Measure(s)RT-PCR confirmed SARS-CoV-2 or rhinovirus infection. ResultsClose contact with a SARS-CoV-2 case (adjusted odds ratio, aOR 4.3, 95% CI 3.7-5.0) and loss of smell/taste (aOR 3.7, 95% CI 3.0-4.5) were the variables most associated with SARS-CoV-2 test positivity, but both attenuated during the Omicron period. Contact with a vaccinated case (aOR 2.4, 95% CI 1.7-3.3) was associated with a lower odds of test positivity than contact with an unvaccinated case (aOR 4.4, 95% CI 2.7-7.3). Sore throat was associated with Omicron infection (aOR 2.3, 95% CI 1.6-3.2) but not Delta. Vaccine effectiveness for participants fully vaccinated with a booster dose was 43% (95% CI 11-63%) for Omicron and 92% (95% CI 61-100%) for Delta. Variables associated with rhinovirus test positivity included age <12 years (aOR 4.0, 95% CI 3.5-4.6) and reporting a runny or stuffy nose (aOR 4.6, 95% CI 4.1-5.2). Race, region, and household crowding were significantly associated with both SARS-CoV-2 and rhinovirus test positivity. Conclusions and RelevanceEstimated risk factors and symptoms associated with SARS-CoV-2 infection have changed over time. There was a shift in reported symptoms between the Delta and Omicron variants as well as reductions in the protection provided by vaccines. Racial and socioeconomic disparities persisted in the third year of SARS-CoV-2 circulation and were also present in rhinovirus infection, although the causal pathways remain unclear. Trends in testing behavior and availability may influence these results. Key Points QuestionWhat are the characteristics associated with SARS-CoV-2 and rhinovirus infection? FindingsIn this test-negative design study of 23,278 participants, reporting close contact with a SARS-CoV-2 case was the strongest risk factor associated with test positivity. Loss of smell and taste was associated with the Delta variant, but not the Omicron variant. Vaccination and prior infection provided greater protection against Delta infection than Omicron Infection. Young age was the strongest predictor of rhinovirus positivity. Sociodemographic disparities were present for both SARS-CoV-2 and rhinovirus. MeaningMonitoring factors associated with respiratory pathogen test positivity remains important to identify at-risk populations in the post-SARS-CoV-2 pandemic period.

Mapping the emergence of SARS-CoV-2 Omicron variants on a university campus

Novel variants continue to emerge in the SARS-CoV-2 pandemic. University testing programs may provide timely epidemiologic and genomic surveillance data to inform public health responses. We conducted testing from September 2021 to February 2022 in a university population under vaccination and indoor mask mandates. A total of 3,048 of 24,393 individuals tested positive for SARS-CoV-2 by RT-PCR; whole genome sequencing identified 209 Delta and 1,730 Omicron genomes of the 1,939 total sequenced. Compared to Delta, Omicron had a shorter median serial interval between genetically identical, symptomatic infections within households (2 versus 6 days, P=0.021). Omicron also demonstrated a greater peak reproductive number (2.4 versus 1.8) and a 1.07 (95% confidence interval: 0.58, 1.57; P<0.0001) higher mean cycle threshold value. Despite near universal vaccination and stringent mitigation measures, Omicron rapidly displaced the Delta variant to become the predominant viral strain and led to a surge in cases in a university population.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance.

Investment in Africa over the past year with regards to SARS-CoV-2 genotyping has led to a massive increase in the number of sequences, exceeding 100,000 genomes generated to track the pandemic on the continent. Our results show an increase in the number of African countries able to sequence within their own borders, coupled with a decrease in sequencing turnaround time. Findings from this genomic surveillance underscores the heterogeneous nature of the pandemic but we observe repeated dissemination of SARS-CoV-2 variants within the continent. Sustained investment for genomic surveillance in Africa is needed as the virus continues to evolve, particularly in the low vaccination landscape. These investments are very crucial for preparedness and response for future pathogen outbreaks. One-Sentence SummaryExpanding Africa SARS-CoV-2 sequencing capacity in a fast evolving pandemic.

SARS-CoV-2 variant dynamics across US states show consistent differences in effective reproduction numbers

Accurately estimating relative transmission rates of SARS-CoV-2 variants remains a scientific and public health priority. Recent studies have used the sample proportions of different variants from genetic sequence data to describe variant frequency dynamics and relative transmission rates, but frequencies alone cannot capture the rich epidemiological behavior of SARS-CoV-2. Here, we extend methods for inferring the effective reproduction number of an epidemic using confirmed case data to jointly estimate variant-specific effective reproduction numbers and frequencies of cocirculating variants using cases and sequences across states in the US from January 2021 to March 2022. Our method can be used to infer structured relationships between effective reproduction numbers across time series allowing us to estimate fixed variant-specific growth advantages. We use this model to estimate the effective reproduction number of SARS-CoV-2 Variants of Concern and Variants of Interest in the United States and estimate consistent growth advantages of particular variants across different locations.

Associations between SARS-CoV-2 variants and risk of COVID-19 hospitalization among confirmed cases in Washington State: a retrospective cohort study

BackgroundThe COVID-19 pandemic is dominated by variant viruses; the resulting impact on disease severity remains unclear. Using a retrospective cohort study, we assessed the hospitalization risk following infection with seven SARS-CoV-2 variants. MethodsOur study includes individuals with positive SARS-CoV-2 RT-PCR in the Washington Disease Reporting System with available viral genome data, from December 1, 2020 to January 14, 2022. The analysis was restricted to cases with specimens collected through sentinel surveillance. Using a Cox proportional hazards model with mixed effects, we estimated hazard ratios (HR) for hospitalization risk following infection with a variant, adjusting for age, sex, calendar week, and vaccination. Findings58,848 cases were sequenced through sentinel surveillance, of which 1705 (2.9%) were hospitalized due to COVID-19. Higher hospitalization risk was found for infections with Gamma (HR 3.20, 95%CI 2.40-4.26), Beta (HR 2.85, 95%CI 1.56-5.23), Delta (HR 2.28 95%CI 1.56-3.34) or Alpha (HR 1.64, 95%CI 1.29-2.07) compared to infections with ancestral lineages; Omicron (HR 0.92, 95%CI 0.56-1.52) showed no significant difference in risk. Following Alpha, Gamma, or Delta infection, unvaccinated patients show higher hospitalization risk, while vaccinated patients show no significant difference in risk, both compared to unvaccinated, ancestral lineage cases. Hospitalization risk following Omicron infection is lower with vaccination. ConclusionInfection with Alpha, Gamma, or Delta results in a higher hospitalization risk, with vaccination attenuating that risk. Our findings support hospital preparedness, vaccination, and genomic surveillance. SummaryHospitalization risk following infection with SARS-CoV-2 variant remains unclear. We find a higher hospitalization risk in cases infected with Alpha, Beta, Gamma, and Delta, but not Omicron, with vaccination lowering risk. Our findings support hospital preparedness, vaccination, and genomic surveillance.

Rapid and parallel adaptive mutations in spike S1 drive clade success in SARS-CoV-2

Given the importance of variant SARS-CoV-2 viruses with altered receptor-binding or antigenic phenotypes, we sought to quantify the degree to which adaptive evolution is driving accumulation of mutations in the SARS-CoV-2 genome. Here we assessed adaptive evolution across genes in the SARS-CoV-2 genome by correlating clade growth with mutation accumulation as well as by comparing rates of nonsynonymous to synonymous divergence, clustering of mutations across the SARS-CoV-2 phylogeny and degree of convergent evolution of individual mutations. We find that spike S1 is the focus of adaptive evolution, but also identify positively-selected mutations in other genes that are sculpting the evolutionary trajectory of SARS-CoV-2. Adaptive changes in S1 accumulated rapidly, resulting in a remarkably high ratio of nonsynonymous to synonymous divergence that is 2.5X greater than that observed in HA1 at the beginning of the 2009 H1N1 pandemic.

Recombination patterns in coronaviruses

As shown during the SARS-CoV-2 pandemic, phylogenetic and phylodynamic methods are essential tools to study the spread and evolution of pathogens. One of the central assumptions of these methods is that the shared history of pathogens isolated from different hosts can be described by a branching phylogenetic tree. Recombination breaks this assumption. This makes it problematic to apply phylogenetic methods to study recombining pathogens, including, for example, coronaviruses. Here, we introduce a Markov chain Monte Carlo approach that allows inference of recombination networks from genetic sequence data under a template switching model of recombination. Using this method, we first show that recombination is extremely common in the evolutionary history of SARS-like coronaviruses. We then show how recombination rates across the genome of the human seasonal coronaviruses 229E, OC43 and NL63 vary with rates of adaptation. This suggests that recombination could be beneficial to fitness of human seasonal coronaviruses. Additionally, this work sets the stage for Bayesian phylogenetic tracking of the spread and evolution of SARS-CoV-2 in the future, even as recombinant viruses become prevalent.

SARS-CoV-2 Epidemiology on a Public University Campus in Washington State

BackgroundTesting programs have been utilized as part of SARS-CoV-2 mitigation strategies on university campuses, and it is not known which strategies successfully identify cases and contain outbreaks. ObjectiveEvaluation of a testing program to control SARS-CoV-2 transmission at a large university. DesignProspective longitudinal study using remote contactless enrollment, daily mobile symptom and exposure tracking, and self-swab sample collection. Individuals were tested if the participant was (1) exposed to a known case, developed new symptoms, or reported high-risk behavior, (2) a member of a group experiencing an outbreak, or (3) at baseline upon enrollment. SettingAn urban, public university during Autumn quarter of 2020 ParticipantsStudents, staff, and faculty. MeasurementsSARS-CoV-2 PCR testing was conducted, and viral genome sequencing was performed. ResultsWe enrolled 16,476 individuals, performed 29,783 SARS-CoV-2 tests, and detected 236 infections. Greek community affiliation was the strongest risk factor for testing positive. 75.0% of positive cases reported at least one of the following: symptoms (60.8%), exposure (34.7%), or high-risk behaviors (21.5%). 88.1% of viral genomes (52/59) sequenced from Greek-affiliated students were genetically identical to at least one other genome detected, indicative of rapid SARS-CoV-2 spread within this group, compared to 37.9% (11/29) of genomes from non-Greek students and employees. LimitationsObservational study. ConclusionIn a setting of limited resources during a pandemic, we prioritized testing of individuals with symptoms and high-risk exposure during outbreaks. Rapid spread of SARS- CoV-2 occurred within outbreaks without evidence of further spread to the surrounding community. A testing program focused on high-risk populations may be effective as part of a comprehensive university-wide mitigation strategy to control the SARS-CoV-2 pandemic.

Specific allelic discrimination of N501Y and other SARS-CoV-2 mutations by ddPCR detects B.1.1.7 lineage in Washington State

Real-time epidemiological tracking of variants of interest can help limit the spread of more contagious forms of SARS-CoV-2, such as those containing the N501Y mutation. Typically, genetic sequencing is required to be able to track variants of interest in real-time. However, sequencing can take time and may not be accessible in all laboratories. Genotyping by RT-ddPCR offers an alternative to sequencing to rapidly detect variants of concern through discrimination of specific mutations such as N501Y that is associated with increased transmissibility. Here we describe the first cases of the B.1.1.7 lineage of SARS-CoV-2 detected in Washington State by using a combination of RT-PCR, RT-ddPCR, and next-generation sequencing. We screened 1,035 samples positive for SARS-CoV-2 by our CDC-based laboratory developed assay using ThermoFishers multiplex RT-PCR COVID-19 assay over four weeks from late December 2020 to early January 2021. S gene dropout candidates were subsequently assayed by RT-ddPCR to confirm four mutations within the S gene associated with the B.1.1.7 lineage: a deletion at amino acid (AA) 69-70 (ACATGT), deletion at AA 145, (TTA), N501Y mutation (TAT), and S982A mutation (GCA). All four targets were detected in two specimens, and follow-up sequencing revealed a total of 10 mutations in the S gene and phylogenetic clustering within the B.1.1.7 lineage. As variants of concern become increasingly prevalent, molecular diagnostic tools like RT-ddPCR can be utilized to quickly, accurately, and sensitively distinguish more contagious lineages of SARS-CoV-2.

A Novel SARS-CoV-2 Variant of Concern, B.1.526, Identified in New York

Recent months have seen surges of SARS-CoV-2 infection across the globe with considerable viral evolution1-3. Extensive mutations in the spike protein may threaten efficacy of vaccines and therapeutic monoclonal antibodies4. Two signature mutations of concern are E484K, which plays a crucial role in the loss of neutralizing activity of antibodies, and N501Y, a driver of rapid worldwide transmission of the B.1.1.7 lineage. Here, we report the emergence of variant lineage B.1.526 that contains E484K and its alarming rise to dominance in New York City in early 2021. This variant is partially or completely resistant to two therapeutic monoclonal antibodies in clinical use and less susceptible to neutralization by convalescent plasma or vaccinee sera, posing a modest antigenic challenge. The B.1.526 lineage has now been reported from all 50 states in the US and numerous other countries. B.1.526 rapidly replaced earlier lineages in New York upon its emergence, with an estimated transmission advantage of 35%. Such transmission dynamics, together with the relative antibody resistance of its E484K sub-lineage, likely contributed to the sharp rise and rapid spread of B.1.526. Although SARS-CoV-2 B.1.526 initially outpaced B.1.1.7 in the region, its growth subsequently slowed concurrent with the rise of B.1.1.7 and ensuing variants.

Viral genome sequencing places White House COVID-19 outbreak into phylogenetic context

In October 2020, an outbreak of at least 50 COVID-19 cases was reported surrounding individuals employed at or visiting the White House. Here, we applied genomic epidemiology to investigate the origins of this outbreak. We enrolled two individuals with exposures linked to the White House COVID-19 outbreak into an IRB-approved research study and sequenced their SARS-CoV-2 infections. We find these viral sequences are identical to each other, but are distinct from over 190,000 publicly available SARS-CoV-2 genomes. These genomes fall as part of a lineage circulating in the USA since April or May 2020 and detected in Virginia and Michigan. Looking forwards, sequencing of additional community SARS-CoV-2 infections collected in the USA prior to October 2020 may shed further light on its geographic ancestry. In sequencing of SARS-CoV-2 infections collected after October 2020, it may be possible to identify infections that likely descend from the White House COVID-19 outbreak.

Evidence for adaptive evolution in the receptor-binding domain of seasonal coronaviruses

Seasonal coronaviruses (OC43, 229E, NL63 and HKU1) are endemic to the human population, regularly infecting and reinfecting humans while typically causing asymptomatic to mild respiratory infections. It is not known to what extent reinfection by these viruses is due to waning immune memory or antigenic drift of the viruses. Here, we address the influence of antigenic drift on immune evasion of seasonal coronaviruses. We provide evidence that at least two of these viruses, OC43 and 229E, are undergoing adaptive evolution in regions of the viral spike protein that are exposed to human humoral immunity. This suggests that reinfection may be due, in part, to positively-selected genetic changes in these viruses that enable them to escape recognition by the immune system. It is possible that, as with seasonal influenza, these adaptive changes in antigenic regions of the virus would necessitate continual reformulation of a vaccine made against them.

Viral genomes reveal patterns of the SARS-CoV-2 outbreak in Washington State

The rapid spread of SARS-CoV-2 has gravely impacted societies around the world. Outbreaks in different parts of the globe are shaped by repeated introductions of new lineages and subsequent local transmission of those lineages. Here, we sequenced 3940 SARS-CoV-2 viral genomes from Washington State to characterize how the spread of SARS-CoV-2 in Washington State (USA) was shaped by differences in timing of mitigation strategies across counties, as well as by repeated introductions of viral lineages into the state. Additionally, we show that the increase in frequency of a potentially more transmissible viral variant (614G) over time can potentially be explained by regional mobility differences and multiple introductions of 614G, but not the other variant (614D) into the state. At an individual level, we see evidence of higher viral loads in patients infected with the 614G variant. However, using clinical records data, we do not find any evidence that the 614G variant impacts clinical severity or patient outcomes. Overall, this suggests that at least to date, the behavior of individuals has been more important in shaping the course of the pandemic than changes in the virus. One Sentence SummaryLocal outbreak dynamics of SARS-CoV-2 in Washington State (USA) were driven by regionally different mitigation measures and repeated introductions of unique viral variants with different viral loads.

Cryptic transmission of SARS-CoV-2 in Washington State

Following its emergence in Wuhan, China, in late November or early December 2019, the SARS-CoV-2 virus has rapidly spread throughout the world. Genome sequencing of SARS-CoV-2 strains allows for the reconstruction of transmission history connecting these infections. Here, we analyze 346 SARS-CoV-2 genomes from samples collected between 20 February and 15 March 2020 from infected patients in Washington State, USA. We found that the large majority of SARS-CoV-2 infections sampled during this time frame appeared to have derived from a single introduction event into the state in late January or early February 2020 and subsequent local spread, indicating cryptic spread of COVID-19 before active community surveillance was implemented. We estimate a common ancestor of this outbreak clade as occurring between 18 January and 9 February 2020. From genomic data, we estimate an exponential doubling between 2.4 and 5.1 days. These results highlight the need for large-scale community surveillance for SARS-CoV-2 and the power of pathogen genomics to inform epidemiological understanding.

Preliminary support for a 'dry swab, extraction free' protocol for SARS-CoV-2 testing via RT-qPCR

Structured AbstractO_ST_ABSBackgroundC_ST_ABSThe urgent need for massively scaled clinical testing for SARS-CoV-2, along with global shortages of critical reagents and supplies, has necessitated development of streamlined laboratory testing protocols. Conventional nucleic acid testing for SARS-CoV-2 involves collection of a clinical specimen with a nasopharyngeal swab in transport medium, nucleic acid extraction, and quantitative reverse transcription PCR (RT-qPCR) (1). As testing has scaled across the world, the global supply chain has buckled, rendering testing reagents and materials scarce (2). To address shortages, we developed SwabExpress, an end-to-end protocol developed to employ mass produced anterior nares swabs and bypass the requirement for transport media and nucleic acid extraction. MethodsWe evaluated anterior nares swabs, transported dry and eluted in low-TE buffer as a direct-to-RT-qPCR alternative to extraction-dependent viral transport media. We validated our protocol of using heat treatment for viral activation and added a proteinase K digestion step to reduce amplification interference. We tested this protocol across archived and prospectively collected swab specimens to fine-tune test performance. ResultsAfter optimization, SwabExpress has a low limit of detection at 2-4 molecules/uL, 100% sensitivity, and 99.4% specificity when compared side-by-side with a traditional RT-qPCR protocol employing extraction. On real-world specimens, SwabExpress outperforms an automated extraction system while simultaneously reducing cost and hands-on time. ConclusionSwabExpress is a simplified workflow that facilitates scaled testing for COVID-19 without sacrificing test performance. It may serve as a template for the simplification of PCR-based clinical laboratory tests, particularly in times of critical shortages during pandemics.

A Genomic Survey of SARS-CoV-2 Reveals Multiple Introductions into Northern California without a Predominant Lineage

The COVID-19 pandemic caused by the novel coronavirus SARS-CoV-2 has spread globally, resulting in >300,000 reported cases worldwide as of March 21st, 2020. Here we investigate the genetic diversity and genomic epidemiology of SARS-CoV-2 in Northern California using samples from returning travelers, cruise ship passengers, and cases of community transmission with unclear infection sources. Virus genomes were sampled from 29 patients diagnosed with COVID-19 infection from Feb 3rd through Mar 15th. Phylogenetic analyses revealed at least 8 different SARS-CoV-2 lineages, suggesting multiple independent introductions of the virus into the state. Virus genomes from passengers on two consecutive excursions of the Grand Princess cruise ship clustered with those from an established epidemic in Washington State, including the WA1 genome representing the first reported case in the United States on January 19th. We also detected evidence for presumptive transmission of SARS-CoV-2 lineages from one community to another. These findings suggest that cryptic transmission of SARS-CoV-2 in Northern California to date is characterized by multiple transmission chains that originate via distinct introductions from international and interstate travel, rather than widespread community transmission of a single predominant lineage. Rapid testing and contact tracing, social distancing, and travel restrictions are measures that will help to slow SARS-CoV-2 spread in California and other regions of the USA.

Effects of weather-related social distancing on city-scale transmission of respiratory viruses

BACKGROUNDUnusually high snowfall in western Washington State in February 2019 led to widespread school and workplace closures. We assessed the impact of social distancing caused by this extreme weather event on the transmission of respiratory viruses. METHODSResidual specimens from patients evaluated for acute respiratory illness at hospitals in the Seattle metropolitan area were screened for a panel of respiratory viruses. Transmission models were fit to each virus, with disruption of contact rates and care-seeking informed by data on local traffic volumes and hospital admissions. RESULTSDisruption in contact patterns reduced effective contact rates during the intervention period by 16% to 95%, and cumulative disease incidence through the remainder of the season by 3% to 9%. Incidence reductions were greatest for viruses that were peaking when the disruption occurred and least for viruses in early epidemic phase. CONCLUSIONHigh-intensity, short-duration social distancing measures may substantially reduce total incidence in a respiratory virus epidemic if implemented near the epidemic peak. One sentence summaryDisruptions of school and work due to heavy snowfall in the Seattle metro area reduced the total size of respiratory virus epidemics by up to 9%.