Local-scale phylodynamics reveal differential community impact of SARS-CoV-2 in a 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.

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.

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

Background

The coronavirus disease 2019 (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 7 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants.

Methods

Our study includes individuals with positive SARS-CoV-2 reverse transcription polymerase chain reaction (RT-PCR) in the Washington Disease Reporting System with available viral genome data, from 1 December 2020 to 14 January 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.

Results

In total, 58 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% confidence interval [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 .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.

Conclusions

Infection 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.

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

The rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has gravely affected societies around the world. Outbreaks in different parts of the globe have been shaped by repeated introductions of new viral lineages and subsequent local transmission of those lineages. Here, we sequenced 3940 SARS-CoV-2 viral genomes from Washington State (USA) to characterize how the spread of SARS-CoV-2 in Washington State in early 2020 was shaped by differences in timing of mitigation strategies across counties and by repeated introductions of viral lineages into the state. In addition, 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 observed evidence of higher viral loads in patients infected with the 614G variant. However, using clinical records data, we did not find any evidence that the 614G variant affects clinical severity or patient outcomes. Overall, this suggests that with regard to D614G, the behavior of individuals has been more important in shaping the course of the pandemic in Washington State than this variant of the virus.

Cryptic transmission of SARS-CoV-2 in Washington state

After its emergence in Wuhan, China, in late November or early December 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus rapidly spread globally. Genome sequencing of SARS-CoV-2 allows the reconstruction of its transmission history, although this is contingent on sampling. We analyzed 453 SARS-CoV-2 genomes collected between 20 February and 15 March 2020 from infected patients in Washington state in the United States. We find that most SARS-CoV-2 infections sampled during this time derive from a single introduction in late January or early February 2020, which subsequently spread locally before active community surveillance was implemented.

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.

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. On March 11, 2020, the World Health Organization declared Coronavirus Disease 2019 (COVID-19) a pandemic. 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, strongly suggesting cryptic spread of COVID-19 during the months of January and February 2020, 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 introductions and spread and the power of pathogen genomics to inform epidemiological understanding.

Single suture craniosynostosis: Identification of rare variants in genes associated with syndromic forms

We report RNA-Sequencing results on a cohort of patients with single suture craniosynostosis and demonstrate significant enrichment of heterozygous, rare, and damaging variants among key craniosynostosis-related genes. Genetic burden analysis identified a significant increase in damaging variants in ATR, EFNA4, ERF, MEGF8, SCARF2, and TGFBR2. Of 391 participants, 15% were found to have damaging and potentially causal variants in 29 genes. We observed transmission in 96% of the affected individuals, and thus penetrance, epigenetics, and oligogenic factors need to be considered when recommending genetic testing in patients with nonsyndromic craniosynostosis.

Characterization of the zinc-containing metalloprotease encoded by zpx and development of a species-specific detection method for Enterobacter sakazakii

Enterobacter sakazakii causes a severe form of neonatal meningitis that occurs as sporadic cases as well as outbreaks. The disease has been epidemiologically associated with consumption of reconstituted, dried infant formulas. Very little information is available regarding pathogenicity of the organism and production of virulence factors. Clinical and environmental strains were screened for production of factors which have activity against Chinese hamster ovary (CHO) cells in tissue culture. Polymyxin B lysate and sonicate preparations but not culture supernatants from the strains caused "rounding" of CHO cells. Subsequent studies showed that the CHO cell-rounding factor is a proteolytic enzyme that has activity against azocasein. The cell-bound protease was isolated by using a combination of polymyxin B lysis, followed by sonication of cells harvested from tryptone broth. The protease was purified to homogeneity by sequential ammonium sulfate precipitation, gel filtration chromatography with Sephadex G-100, hydrophobic interaction chromatography with phenyl-Sepharose CL-4B, and a second gel filtration with Sephadex G-100. In addition to activity against azocasein, the purified protease also exhibits activity against azocoll and insoluble casein but not elastin. The protease has a molecular weight of 38,000 and an isoelectric point of 4.4. It is heat labile and for maximal activity against azocasein has an optimum temperature of 37 degrees C and a pH range of 5 to 7. Proteolytic activity is inhibited by ortho-phenanthroline and Zincov but is not affected by phenylmethylsulfonyl fluoride, N-ethylmaleimide, and trypsin inhibitors, which demonstrates that the protease is a zinc-containing metalloprotease. The metalloprotease does not hemagglutinate chicken or sheep erythrocytes. Twenty-three to 27 of the first 42 N-terminal amino acid residues of the metalloprotease are identical to proteases produced by Serratia proteamaculans, Pectobacterium carotovorum, and Anabaena sp. PCR analysis using primers designed from a consensus nucleotide sequence showed that 135 E. sakazakii strains possessed the metalloprotease gene, zpx, and 25 non-E. sakazakii strains did not. The cloned zpx gene of strain 29544 consists of 1,026 nucleotides, and the deduced amino acid sequence of the metalloprotease has 341 amino acid residues, which corresponds to a theoretical protein size of 37,782 with a theoretical pI of 5.23. The sequence possesses three well-characterized zinc-binding and active-site motifs present in other bacterial zinc metalloproteases.