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.

SARS-CoV-2 evolution during prolonged infection in immunocompromised patients.

Prolonged infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in immunocompromised patients provides an opportunity for viral evolution, potentially leading to the generation of new pathogenic variants. To investigate the pathways of viral evolution, we carried out a study on five patients experiencing prolonged SARS-CoV-2 infection (quantitative polymerase chain reaction-positive for 79-203 days) who were immunocompromised due to treatment for lymphoma or solid organ transplantation. For each timepoint analyzed, we generated at least two independent viral genome sequences to assess the heterogeneity and control for sequencing error. Four of the five patients likely had prolonged infection; the fifth apparently experienced a reinfection. The rates of accumulation of substitutions in the viral genome per day were higher in hospitalized patients with prolonged infection than those estimated for the community background. The spike coding region accumulated a significantly greater number of unique mutations than other viral coding regions, and the mutation density was higher. Two patients were treated with monoclonal antibodies (bebtelovimab and sotrovimab); by the next sampled timepoint, each virus population showed substitutions associated with monoclonal antibody resistance as the dominant forms (spike K444N and spike E340D). All patients received remdesivir, but remdesivir-resistant substitutions were not detected. These data thus help elucidate the trends of emergence, evolution, and selection of mutational variants within long-term infected immunocompromised individuals. IMPORTANCE: SARS-CoV-2 is responsible for a global pandemic, driven in part by the emergence of new viral variants. Where do these new variants come from? One model is that long-term viral persistence in infected individuals allows for viral evolution in response to host pressures, resulting in viruses more likely to replicate efficiently in humans. In this study, we characterize replication in several hospitalized and long-term infected individuals, documenting efficient pathways of viral evolution.

Avian H7N9 influenza viruses are evolutionarily constrained by stochastic processes during replication and transmission in mammals.

H7N9 avian influenza viruses (AIVs) have caused over 1,500 documented human infections since emerging in 2013. Although wild-type H7N9 AIVs can be transmitted by respiratory droplets in ferrets, they have not yet caused widespread outbreaks in humans. Previous studies have revealed molecular determinants of H7N9 AIV host switching, but little is known about potential evolutionary constraints on this process. Here, we compare patterns of sequence evolution for H7N9 AIV and mammalian H1N1 viruses during replication and transmission in ferrets. We show that three main factors-purifying selection, stochasticity, and very narrow transmission bottlenecks-combine to severely constrain the ability of H7N9 AIV to effectively adapt to mammalian hosts in isolated, acute spillover events. We find rare evidence of natural selection favoring new, potentially mammal-adapting mutations within ferrets but no evidence of natural selection acting during transmission. We conclude that human-adapted H7N9 viruses are unlikely to emerge during typical spillover infections. Our findings are instead consistent with a model in which the emergence of a human-transmissible virus would be a rare and unpredictable, though highly consequential, 'jackpot' event. Strategies to control the total number of spillover infections will limit opportunities for the virus to win this evolutionary lottery.

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.

A multivalent nucleoside-modified mRNA vaccine against all known influenza virus subtypes.

Seasonal influenza vaccines offer little protection against pandemic influenza virus strains. It is difficult to create effective prepandemic vaccines because it is uncertain which influenza virus subtype will cause the next pandemic. In this work, we developed a nucleoside-modified messenger RNA (mRNA)-lipid nanoparticle vaccine encoding hemagglutinin antigens from all 20 known influenza A virus subtypes and influenza B virus lineages. This multivalent vaccine elicited high levels of cross-reactive and subtype-specific antibodies in mice and ferrets that reacted to all 20 encoded antigens. Vaccination protected mice and ferrets challenged with matched and mismatched viral strains, and this protection was at least partially dependent on antibodies. Our studies indicate that mRNA vaccines can provide protection against antigenically variable viruses by simultaneously inducing antibodies against multiple antigens.

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.

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.

Repeated introductions and intensive community transmission fueled a mumps virus outbreak in Washington State.

In 2016/2017, Washington State experienced a mumps outbreak despite high childhood vaccination rates, with cases more frequently detected among school-aged children and members of the Marshallese community. We sequenced 166 mumps virus genomes collected in Washington and other US states, and traced mumps introductions and transmission within Washington. We uncover that mumps was introduced into Washington approximately 13 times, primarily from Arkansas, sparking multiple co-circulating transmission chains. Although age and vaccination status may have impacted transmission, our data set could not quantify their precise effects. Instead, the outbreak in Washington was overwhelmingly sustained by transmission within the Marshallese community. Our findings underscore the utility of genomic data to clarify epidemiologic factors driving transmission and pinpoint contact networks as critical for mumps transmission. These results imply that contact structures and historic disparities may leave populations at increased risk for respiratory virus disease even when a vaccine is effective and widely used.

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.

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.

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.

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.

Genomic epidemiology supports multiple introductions and cryptic transmission of Zika virus in Colombia.

BACKGROUND: Colombia was the second most affected country during the American Zika virus (ZIKV) epidemic, with over 109,000 reported cases. Despite the scale of the outbreak, limited genomic sequence data were available from Colombia. We sought to sequence additional samples and use genomic epidemiology to describe ZIKV dynamics in Colombia. METHODS: We sequenced ZIKV genomes directly from clinical diagnostic specimens and infected Aedes aegypti samples selected to cover the temporal and geographic breadth of the Colombian outbreak. We performed phylogeographic analysis of these genomes, along with other publicly-available ZIKV genomes from the Americas, to estimate the frequency and timing of ZIKV introductions to Colombia. RESULTS: We attempted PCR amplification on 184 samples; 19 samples amplified sufficiently to perform sequencing. Of these, 8 samples yielded sequences with at least 50% coverage. Our phylogeographic reconstruction indicates two separate introductions of ZIKV to Colombia, one of which was previously unrecognized. We find that ZIKV was first introduced to Colombia in February 2015 (95%CI: Jan 2015 - Apr 2015), corresponding to 5 to 8 months of cryptic ZIKV transmission prior to confirmation in September 2015. Despite the presence of multiple introductions, we find that the majority of Colombian ZIKV diversity descends from a single introduction. We find evidence for movement of ZIKV from Colombia into bordering countries, including Peru, Ecuador, Panama, and Venezuela. CONCLUSIONS: Similarly to genomic epidemiological studies of ZIKV dynamics in other countries, we find that ZIKV circulated cryptically in Colombia. More accurately dating when ZIKV was circulating refines our definition of the population at risk. Additionally, our finding that the majority of ZIKV transmission within Colombia was attributable to transmission between individuals, rather than repeated travel-related importations, indicates that improved detection and control might have succeeded in limiting the scale of the outbreak within Colombia.

Comprehensive mapping of adaptation of the avian influenza polymerase protein PB2 to humans.

Viruses like influenza are infamous for their ability to adapt to new hosts. Retrospective studies of natural zoonoses and passaging in the lab have identified a modest number of host-adaptive mutations. However, it is unclear if these mutations represent all ways that influenza can adapt to a new host. Here we take a prospective approach to this question by completely mapping amino-acid mutations to the avian influenza virus polymerase protein PB2 that enhance growth in human cells. We identify numerous previously uncharacterized human-adaptive mutations. These mutations cluster on PB2's surface, highlighting potential interfaces with host factors. Some previously uncharacterized adaptive mutations occur in avian-to-human transmission of H7N9 influenza, showing their importance for natural virus evolution. But other adaptive mutations do not occur in nature because they are inaccessible via single-nucleotide mutations. Overall, our work shows how selection at key molecular surfaces combines with evolutionary accessibility to shape viral host adaptation.

Tracking dengue virus type 1 genetic diversity during lineage replacement in an hyperendemic area in Colombia.

Dengue virus (DENV) is a flavivirus responsible for the most common and burdensome arthropod-borne viral disease of humans[1]. DENV evolution has been extensively studied on broad geographic and time scales, using sequences from a single gene[2,3]. It is believed that DENV evolution in humans is dominated primarily by purifying selection due to the constraint of maintaining fitness in both humans and mosquitoes[4,5]. Few studies have explored DENV evolutionary dynamics using whole genome sequences, nor have they explored changes in viral diversity that occur during intra-epidemic periods. We used deep sequencing of the viral coding region to characterize DENV-1 evolution in a Colombian population sampled during two high-prevalence dengue seasons in which serotype dominance shifted. Our data demonstrate patterns of strain extinction and replacement within DENV-1 as its prevalence waned and DENV-3 became established. A comparison of whole-genome versus single-gene-based phylogenetic analyses highlights an important difference in evolutionary patterns. We report a trend of higher nonsynonymous to synonymous diversity ratios among non-structural (NS) genes, and statistically significantly higher values among these ratios in the NS1 gene after DENV-1 strain replacement. These results suggest that positive selection could be driving DENV evolution within individual communities. Signals of positive selection coming from distinct samples may be drowned out when combining multiple regions with differing patterns of endemic transmission as commonly done by large-scale geo-temporal assessments. Here, we frame our findings within a small, local transmission history which aids significance. Moreover, these data suggest that the NS1 gene, rather than the E gene, may be a target of positive selection, although not mutually exclusive, and potentially useful sentinel of adaptive changes at the population level.

Diversity of Influenza A(H5N1) Viruses in Infected Humans, Northern Vietnam, 2004-2010.

Influenza viruses exist in each host as a collection of genetically diverse variants, which might enhance their adaptive potential. To assess the genetic and functional diversity of highly pathogenic avian influenza A(H5N1) viruses within infected humans, we used deep-sequencing methods to characterize samples obtained from infected patients in northern Vietnam during 2004-2010 on different days after infection, from different anatomic sites, or both. We detected changes in virus genes that affected receptor binding, polymerase activity, or interferon antagonism, suggesting that these factors could play roles in influenza virus adaptation to humans. However, the frequency of most of these mutations remained low in the samples tested, implying that they were not efficiently selected within these hosts. Our data suggest that adaptation of influenza A(H5N1) viruses is probably stepwise and depends on accumulating combinations of mutations that alter function while maintaining fitness.

Within-Host Evolution of Human Influenza Virus.

The rapid global evolution of influenza virus begins with mutations that arise de novo in individual infections, but little is known about how evolution occurs within hosts. We review recent progress in understanding how and why influenza viruses evolve within human hosts. Advances in deep sequencing make it possible to measure within-host genetic diversity in both acute and chronic influenza infections. Factors like antigenic selection, antiviral treatment, tissue specificity, spatial structure, and multiplicity of infection may affect how influenza viruses evolve within human hosts. Studies of within-host evolution can contribute to our understanding of the evolutionary and epidemiological factors that shape influenza virus's global evolution.

Emergence of Oseltamivir-Resistant H7N9 Influenza Viruses in Immunosuppressed Cynomolgus Macaques.

Antiviral compounds (eg, the neuraminidase inhibitor oseltamivir) are invaluable for the treatment of individuals infected with influenza A viruses of the H7N9 subtype (A[H7N9]), which have infected and killed hundreds of persons. However, oseltamivir treatment often leads to the emergence of resistant viruses in immunocompromised individuals. To better understand the emergence and properties of oseltamivir-resistant A(H7N9) viruses in immunosuppressed individuals, we infected immunosuppressed cynomolgus macaques with an A(H7N9) virus and treated them with oseltamivir. Disease severity and mortality were higher in immunosuppressed than in immunocompetent animals. Oseltamivir treatment at 2 different doses reduced A(H7N9) viral titers in infected animals, but even high-dose oseltamivir did not block viral replication sufficiently to suppress the emergence of resistant variants. Some resistant variants were not appreciably attenuated in cultured cells, but an oseltamivir-resistant A(H7N9) virus did not transmit among ferrets. These findings are useful for the control of A(H7N9) virus infections in clinical settings.

Influenza Evolution: New Insights into an Old Foe.

Influenza viruses steadily evolve to escape detection by antibodies, necessitating vaccine updates. A new study uses a massively parallel approach, deep mutational scanning, to catalogue antibody escape mutations. This approach defines potential pathways of viral evolution, beyond those already observed in natural infections, and may help predict its future directions.

Within-Host Evolution of Simian Arteriviruses in Crab-Eating Macaques.

Simian arteriviruses are a diverse clade of viruses infecting captive and wild nonhuman primates. We recently reported that Kibale red colobus virus 1 (KRCV-1) causes a mild and self-limiting disease in experimentally infected crab-eating macaques, while simian hemorrhagic fever virus (SHFV) causes lethal viral hemorrhagic fever. Here we characterize how these viruses evolved during replication in cell culture and in experimentally infected macaques. During passage in cell culture, 68 substitutions that were localized in open reading frames (ORFs) likely associated with host cell entry and exit became fixed in the KRCV-1 genome. However, we did not detect any strong signatures of selection during replication in macaques. We uncovered patterns of evolution that were distinct from those observed in surveys of wild red colobus monkeys, suggesting that these species may exert different adaptive challenges for KRCV-1. During SHFV infection, we detected signatures of selection on ORF 5a and on a small subset of sites in the genome. Overall, our data suggest that patterns of evolution differ markedly among simian arteriviruses and among host species. IMPORTANCE: Certain RNA viruses can cross species barriers and cause disease in new hosts. Simian arteriviruses are a diverse group of related viruses that infect captive and wild nonhuman primates, with associated disease severity ranging from apparently asymptomatic infections to severe, viral hemorrhagic fevers. We infected nonhuman primate cell cultures and then crab-eating macaques with either simian hemorrhagic fever virus (SHFV) or Kibale red colobus virus 1 (KRCV-1) and assessed within-host viral evolution. We found that KRCV-1 quickly acquired a large number of substitutions in its genome during replication in cell culture but that evolution in macaques was limited. In contrast, we detected selection focused on SHFV ORFs 5a and 5, which encode putative membrane proteins. These patterns suggest that in addition to diverse pathogenic phenotypes, these viruses may also exhibit distinct patterns of within-host evolution both in vitro and in vivo.