Underdetected dispersal and extensive local transmission drove the 2022 mpox epidemic.

The World Health Organization declared mpox a public health emergency of international concern in July 2022. To investigate global mpox transmission and population-level changes associated with controlling spread, we built phylogeographic and phylodynamic models to analyze MPXV genomes from five global regions together with air traffic and epidemiological data. Our models reveal community transmission prior to detection, changes in case reporting throughout the epidemic, and a large degree of transmission heterogeneity. We find that viral introductions played a limited role in prolonging spread after initial dissemination, suggesting that travel bans would have had only a minor impact. We find that mpox transmission in North America began declining before more than 10% of high-risk individuals in the USA had vaccine-induced immunity. Our findings highlight the importance of broader routine specimen screening surveillance for emerging infectious diseases and of joint integration of genomic and epidemiological information for early outbreak control.

Did a Single Amino Acid Change Make Ebola Virus More Virulent?

A mutation in the Ebola virus glycoprotein arose early during the 2013-2016 epidemic and dominated the viral population. Two studies by Diehl et al. and Urbanowicz et al. now reveal that this mutation is associated with higher infectivity to human cells, representing the clearest example of Ebola's functional adaptation to human hosts.

Spike deep mutational scanning helps predict success of SARS-CoV-2 clades.

SARS-CoV-2 variants acquire mutations in the spike protein that promote immune evasion1 and affect other properties that contribute to viral fitness, such as ACE2 receptor binding and cell entry2,3. Knowledge of how mutations affect these spike phenotypes can provide insight into the current and potential future evolution of the virus. Here we use pseudovirus deep mutational scanning4 to measure how more than 9,000 mutations across the full XBB.1.5 and BA.2 spikes affect ACE2 binding, cell entry or escape from human sera. We find that mutations outside the receptor-binding domain (RBD) have meaningfully affected ACE2 binding during SARS-CoV-2 evolution. We also measure how mutations to the XBB.1.5 spike affect neutralization by serum from individuals who recently had SARS-CoV-2 infections. The strongest serum escape mutations are in the RBD at sites 357, 420, 440, 456 and 473; however, the antigenic effects of these mutations vary across individuals. We also identify strong escape mutations outside the RBD; however, many of them decrease ACE2 binding, suggesting they act by modulating RBD conformation. Notably, the growth rates of human SARS-CoV-2 clades can be explained in substantial part by the measured effects of mutations on spike phenotypes, suggesting our data could enable better prediction of viral evolution.

Emergence and expansion of SARS-CoV-2 B.1.526 after identification in New York.

SARS-CoV-2 infections have surged across the globe in recent months, concomitant with considerable viral evolution1-3. Extensive mutations in the spike protein may threaten the efficacy of vaccines and therapeutic monoclonal antibodies4. Two signature spike mutations of concern are E484K, which has 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 the variant lineage B.1.526 (also known as the Iota variant5), which contains E484K, and its rise to dominance in New York City in early 2021. This variant is partially or completely resistant to two therapeutic monoclonal antibodies that are in clinical use and is less susceptible to neutralization by plasma from individuals who had recovered from SARS-CoV-2 infection or serum from vaccinated individuals, posing a modest antigenic challenge. The presence of the B.1.526 lineage has now been reported in all 50 states in the United States and in many other countries. B.1.526 rapidly replaced earlier lineages in New York, with an estimated transmission advantage of 35%. These transmission dynamics, together with the relative antibody resistance of its E484K sub-lineage, are likely to have 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 concurrently with the rise of B.1.1.7 and ensuing variants.

Estimating the reproduction number and transmission heterogeneity from the size distribution of clusters of identical pathogen sequences.

Quantifying transmission intensity and heterogeneity is crucial to ascertain the threat posed by infectious diseases and inform the design of interventions. Methods that jointly estimate the reproduction number R and the dispersion parameter k have however mainly remained limited to the analysis of epidemiological clusters or contact tracing data, whose collection often proves difficult. Here, we show that clusters of identical sequences are imprinted by the pathogen offspring distribution, and we derive an analytical formula for the distribution of the size of these clusters. We develop and evaluate an inference framework to jointly estimate the reproduction number and the dispersion parameter from the size distribution of clusters of identical sequences. We then illustrate its application across a range of epidemiological situations. Finally, we develop a hypothesis testing framework relying on clusters of identical sequences to determine whether a given pathogen genetic subpopulation is associated with increased or reduced transmissibility. Our work provides tools to estimate the reproduction number and transmission heterogeneity from pathogen sequences without building a phylogenetic tree, thus making it easily scalable to large pathogen genome datasets.

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.

Urgent need for a non-discriminatory and non-stigmatizing nomenclature for monkeypox virus.

We propose a novel, non-discriminatory classification of monkeypox virus diversity. Together with the World Health Organization, we named three clades (I, IIa and IIb) in order of detection. Within IIb, the cause of the current global outbreak, we identified multiple lineages (A.1, A.2, A.1.1 and B.1) to support real-time genomic surveillance.

Comparative diagnostic utility of SARS-CoV-2 rapid antigen and molecular testing in a community setting.

BACKGROUND: SARS-CoV-2 antigen-detection rapid diagnostic tests (Ag-RDTs) have become widely utilized but longitudinal characterization of their community-based performance remains incompletely understood. METHODS: This prospective longitudinal study at a large public university in Seattle, WA utilized remote enrollment, online surveys, and self-collected nasal swab specimens to evaluate Ag-RDT performance against real-time reverse transcription polymerase chain reaction (rRT-PCR) in the context of SARS-CoV-2 Omicron. Ag-RDT sensitivity and specificity within 1 day of rRT-PCR were evaluated by symptom status throughout the illness episode and Orf1b cycle threshold (Ct). RESULTS: From February to December 2022, 5,757 participants reported 17,572 Ag-RDT results and completed 12,674 rRT-PCR tests, of which 995 (7.9%) were rRT-PCR-positive. Overall sensitivity and specificity were 53.0% (95% CI: 49.6-56.4%) and 98.8% (98.5-99.0%), respectively. Sensitivity was comparatively higher for Ag-RDTs used 1 day after rRT-PCR (69.0%), 4 to 7 days post-symptom onset (70.1%), and Orf1b Ct ≤20 (82.7%). Serial Ag-RDT sensitivity increased with repeat testing ≥2 (68.5%) and ≥4 (75.8%) days after an initial Ag-RDT-negative result. CONCLUSION: Ag-RDT performance varied by clinical characteristics and temporal testing patterns. Our findings support recommendations for serial testing following an initial Ag-RDT-negative result, especially among recently symptomatic persons or those at high-risk for SARS-CoV-2 infection.

Ebola Virus Transmission Initiated by Relapse of Systemic Ebola Virus Disease.

During the 2018-2020 Ebola virus disease (EVD) outbreak in North Kivu province in the Democratic Republic of Congo, EVD was diagnosed in a patient who had received the recombinant vesicular stomatitis virus-based vaccine expressing a ZEBOV glycoprotein (rVSV-ZEBOV) (Merck). His treatment included an Ebola virus (EBOV)-specific monoclonal antibody (mAb114), and he recovered within 14 days. However, 6 months later, he presented again with severe EVD-like illness and EBOV viremia, and he died. We initiated epidemiologic and genomic investigations that showed that the patient had had a relapse of acute EVD that led to a transmission chain resulting in 91 cases across six health zones over 4 months. (Funded by the Bill and Melinda Gates Foundation and others.).

Changing genomic epidemiology of COVID-19 in long-term care facilities during the 2020-2022 pandemic, Washington State.

BACKGROUND: Long-term care facilities (LTCFs) are vulnerable to disease outbreaks. Here, we jointly analyze SARS-CoV-2 genomic and paired epidemiologic data from LTCFs and surrounding communities in Washington state (WA) to assess transmission patterns during 2020-2022, in a setting of changing policy. We describe sequencing efforts and genomic epidemiologic findings across LTCFs and perform in-depth analysis in a single county. METHODS: We assessed genomic data representativeness, built phylogenetic trees, and conducted discrete trait analysis to estimate introduction sizes over time, and explored selected outbreaks to further characterize transmission events. RESULTS: We found that transmission dynamics among cases associated with LTCFs in WA changed over the course of the COVID-19 pandemic, with variable introduction rates into LTCFs, but decreasing amplification within LTCFs. SARS-CoV-2 lineages circulating in LTCFs were similar to those circulating in communities at the same time. Transmission between staff and residents was bi-directional. CONCLUSIONS: Understanding transmission dynamics within and between LTCFs using genomic epidemiology on a broad scale can assist in targeting policies and prevention efforts. Tracking facility-level outbreaks can help differentiate intra-facility outbreaks from high community transmission with repeated introduction events. Based on our study findings, methods for routine tree building and overlay of epidemiologic data for hypothesis generation by public health practitioners are recommended. Discrete trait analysis added valuable insight and can be considered when representative sequencing is performed. Cluster detection tools, especially those that rely on distance thresholds, may be of more limited use given current data capture and timeliness. Importantly, we noted a decrease in data capture from LTCFs over time. Depending on goals for use of genomic data, sentinel surveillance should be increased or targeted surveillance implemented to ensure available data for analysis.

Sentinel Surveillance System Implementation and Evaluation for SARS-CoV-2 Genomic Data, Washington, USA, 2020-2021.

Genomic data provides useful information for public health practice, particularly when combined with epidemiologic data. However, sampling bias is a concern because inferences from nonrandom data can be misleading. In March 2021, the Washington State Department of Health, USA, partnered with submitting and sequencing laboratories to establish sentinel surveillance for SARS-CoV-2 genomic data. We analyzed available genomic and epidemiologic data during presentinel and sentinel periods to assess representativeness and timeliness of availability. Genomic data during the presentinel period was largely unrepresentative of all COVID-19 cases. Data available during the sentinel period improved representativeness for age, death from COVID-19, outbreak association, long-term care facility-affiliated status, and geographic coverage; timeliness of data availability and captured viral diversity also improved. Hospitalized cases were underrepresented, indicating a need to increase inpatient sampling. Our analysis emphasizes the need to understand and quantify sampling bias in phylogenetic studies and continue evaluation and improvement of public health surveillance systems.

Acute SARS-CoV-2 infections harbor limited within-host diversity and transmit via tight transmission bottlenecks.

The emergence of divergent SARS-CoV-2 lineages has raised concern that novel variants eliciting immune escape or the ability to displace circulating lineages could emerge within individual hosts. Though growing evidence suggests that novel variants arise during prolonged infections, most infections are acute. Understanding how efficiently variants emerge and transmit among acutely-infected hosts is therefore critical for predicting the pace of long-term SARS-CoV-2 evolution. To characterize how within-host diversity is generated and propagated, we combine extensive laboratory and bioinformatic controls with metrics of within- and between-host diversity to 133 SARS-CoV-2 genomes from acutely-infected individuals. We find that within-host diversity is low and transmission bottlenecks are narrow, with very few viruses founding most infections. Within-host variants are rarely transmitted, even among individuals within the same household, and are rarely detected along phylogenetically linked infections in the broader community. These findings suggest that most variation generated within-host is lost during transmission.

Virus genomes reveal factors that spread and sustained the Ebola epidemic.

The 2013-2016 West African epidemic caused by the Ebola virus was of unprecedented magnitude, duration and impact. Here we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region by analysing 1,610 Ebola virus genomes, which represent over 5% of the known cases. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic 'gravity' model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already sown the seeds for an international epidemic, rendering these measures ineffective at curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing that these countries were susceptible to substantial outbreaks but at lower risk of introductions. Finally, we reveal that this large epidemic was a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help to inform interventions in future epidemics.

Associations Between Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Variants and Risk of Coronavirus Disease 2019 (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.

Limited Predictability of Amino Acid Substitutions in Seasonal Influenza Viruses.

Seasonal influenza viruses repeatedly infect humans in part because they rapidly change their antigenic properties and evade host immune responses, necessitating frequent updates of the vaccine composition. Accurate predictions of strains circulating in the future could therefore improve the vaccine match. Here, we studied the predictability of frequency dynamics and fixation of amino acid substitutions. Current frequency was the strongest predictor of eventual fixation, as expected in neutral evolution. Other properties, such as occurrence in previously characterized epitopes or high Local Branching Index (LBI) had little predictive power. Parallel evolution was found to be moderately predictive of fixation. Although the LBI had little power to predict frequency dynamics, it was still successful at picking strains representative of future populations. The latter is due to a tendency of the LBI to be high for consensus-like sequences that are closer to the future than the average sequence. Simulations of models of adapting populations, in contrast, show clear signals of predictability. This indicates that the evolution of influenza HA and NA, while driven by strong selection pressure to change, is poorly described by common models of directional selection such as traveling fitness waves.

Helicobacter pylori diversification during chronic infection within a single host generates sub-populations with distinct phenotypes.

Helicobacter pylori chronically infects the stomach of approximately half of the world's population. Manifestation of clinical diseases associated with H. pylori infection, including cancer, is driven by strain properties and host responses; and as chronic infection persists, both are subject to change. Previous studies have documented frequent and extensive within-host bacterial genetic variation. To define how within-host diversity contributes to phenotypes related to H. pylori pathogenesis, this project leverages a collection of 39 clinical isolates acquired prospectively from a single subject at two time points and from multiple gastric sites. During the six years separating collection of these isolates, this individual, initially harboring a duodenal ulcer, progressed to gastric atrophy and concomitant loss of acid secretion. Whole genome sequence analysis identified 1,767 unique single nucleotide polymorphisms (SNPs) across isolates and a nucleotide substitution rate of 1.3x10-4 substitutions/site/year. Gene ontology analysis identified cell envelope genes among the genes with excess accumulation of nonsynonymous SNPs (nSNPs). A maximum likelihood tree based on genetic similarity clusters isolates from each time point separately. Within time points, there is segregation of subgroups with phenotypic differences in bacterial morphology, ability to induce inflammatory cytokines, and mouse colonization. Higher inflammatory cytokine induction in recent isolates maps to shared polymorphisms in the Cag PAI protein, CagY, while rod morphology in a subgroup of recent isolates mapped to eight mutations in three distinct helical cell shape determining (csd) genes. The presence of subgroups with unique genetic and phenotypic properties suggest complex selective forces and multiple niches within the stomach during chronic infection.

Quantifying within-host diversity of H5N1 influenza viruses in humans and poultry in Cambodia.

Avian influenza viruses (AIVs) periodically cross species barriers and infect humans. The likelihood that an AIV will evolve mammalian transmissibility depends on acquiring and selecting mutations during spillover, but data from natural infection is limited. We analyze deep sequencing data from infected humans and domestic ducks in Cambodia to examine how H5N1 viruses evolve during spillover. Overall, viral populations in both species are predominated by low-frequency (<10%) variation shaped by purifying selection and genetic drift, and half of the variants detected within-host are never detected on the H5N1 virus phylogeny. However, we do detect a subset of mutations linked to human receptor binding and replication (PB2 E627K, HA A150V, and HA Q238L) that arose in multiple, independent humans. PB2 E627K and HA A150V were also enriched along phylogenetic branches leading to human infections, suggesting that they are likely human-adaptive. Our data show that H5N1 viruses generate putative human-adapting mutations during natural spillover infection, many of which are detected at >5% frequency within-host. However, short infection times, genetic drift, and purifying selection likely restrict their ability to evolve extensively during a single infection. Applying evolutionary methods to sequence data, we reveal a detailed view of H5N1 virus adaptive potential, and develop a foundation for studying host-adaptation in other zoonotic viruses.

Characterising the epidemic spread of influenza A/H3N2 within a city through phylogenetics.

Infecting large portions of the global population, seasonal influenza is a major burden on societies around the globe. While the global source sink dynamics of the different seasonal influenza viruses have been studied intensively, its local spread remains less clear. In order to improve our understanding of how influenza is transmitted on a city scale, we collected an extremely densely sampled set of influenza sequences alongside patient metadata. To do so, we sequenced influenza viruses isolated from patients of two different hospitals, as well as private practitioners in Basel, Switzerland during the 2016/2017 influenza season. The genetic sequences reveal that repeated introductions into the city drove the influenza season. We then reconstruct how the effective reproduction number changed over the course of the season. While we did not find that transmission dynamics in Basel correlate with humidity or school closures, we did find some evidence that it may positively correlated with temperature. Alongside the genetic sequence data that allows us to see how individual cases are connected, we gathered patient information, such as the age or household status. Zooming into the local transmission outbreaks suggests that the elderly were to a large extent infected within their own transmission network. In the remaining transmission network, our analyses suggest that school-aged children likely play a more central role than pre-school aged children. These patterns will be valuable to plan interventions combating the spread of respiratory diseases within cities given that similar patterns are observed for other influenza seasons and cities.

Persistent HIV-1 replication maintains the tissue reservoir during therapy.

Lymphoid tissue is a key reservoir established by HIV-1 during acute infection. It is a site associated with viral production, storage of viral particles in immune complexes, and viral persistence. Although combinations of antiretroviral drugs usually suppress viral replication and reduce viral RNA to undetectable levels in blood, it is unclear whether treatment fully suppresses viral replication in lymphoid tissue reservoirs. Here we show that virus evolution and trafficking between tissue compartments continues in patients with undetectable levels of virus in their bloodstream. We present a spatial and dynamic model of persistent viral replication and spread that indicates why the development of drug resistance is not a foregone conclusion under conditions in which drug concentrations are insufficient to completely block virus replication. These data provide new insights into the evolutionary and infection dynamics of the virus population within the host, revealing that HIV-1 can continue to replicate and replenish the viral reservoir despite potent antiretroviral therapy.

Characteristics of COVID-19 in Homeless Shelters : A Community-Based Surveillance Study.

BACKGROUND: Homeless shelters are a high-risk setting for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission because of crowding and shared hygiene facilities. OBJECTIVE: To investigate SARS-CoV-2 case counts across several adult and family homeless shelters in a major metropolitan area. DESIGN: Cross-sectional, community-based surveillance study. (ClinicalTrials.gov: NCT04141917). SETTING: 14 homeless shelters in King County, Washington. PARTICIPANTS: A total of 1434 study encounters were done in shelter residents and staff, regardless of symptoms. INTERVENTION: 2 strategies were used for SARS-CoV-2 testing: routine surveillance and contact tracing ("surge testing") events. MEASUREMENTS: The primary outcome measure was test positivity rate of SARS-CoV-2 infection at shelters, determined by dividing the number of positive cases by the total number of participant encounters, regardless of symptoms. Sociodemographic, clinical, and virologic variables were assessed as correlates of viral positivity. RESULTS: Among 1434 encounters, 29 (2% [95% CI, 1.4% to 2.9%]) cases of SARS-CoV-2 infection were detected across 5 shelters. Most (n = 21 [72.4%]) were detected during surge testing events rather than routine surveillance, and most (n = 21 [72.4% {CI, 52.8% to 87.3%}]) were asymptomatic at the time of sample collection. Persons who were positive for SARS-CoV-2 were more frequently aged 60 years or older than those without SARS-CoV-2 (44.8% vs. 15.9%). Eighty-six percent of persons with positive test results slept in a communal space rather than in a private or shared room. LIMITATION: Selection bias due to voluntary participation and a relatively small case count. CONCLUSION: Active surveillance and surge testing were used to detect multiple cases of asymptomatic and symptomatic SARS-CoV-2 infection in homeless shelters. The findings suggest an unmet need for routine viral testing outside of clinical settings for homeless populations. PRIMARY FUNDING SOURCE: Gates Ventures.