Influenza A virus (IAV) genetic exchange through reassortment has the potential to accelerate viral evolution and has played a critical role in the generation of multiple pandemic strains. For reassortment to occur, distinct viruses must co-infect the same cell. The spatio-temporal dynamics of viral dissemination within an infected host therefore define opportunity for reassortment. Here, we used wild type and synonymously barcoded variant viruses of a pandemic H1N1 strain to examine the within-host viral dynamics that govern reassortment in guinea pigs, ferrets and swine. The first two species are well-established models of human influenza, while swine are a natural host and a frequent conduit for cross-species transmission and reassortment. Our results show reassortment to be pervasive in all three hosts but less frequent in swine than in ferrets and guinea pigs. In ferrets, tissue-specific differences in the opportunity for reassortment are also evident, with more reassortants detected in the nasal tract than the lower respiratory tract. While temporal trends in viral diversity are limited, spatial patterns are clear, with heterogeneity in the viral genotypes detected at distinct anatomical sites revealing extensive compartmentalization of reassortment and replication. Our data indicate that the dynamics of viral replication in mammals allow diversification through reassortment but that the spatial compartmentalization of variants likely shapes their evolution and onward transmission.

Periodic resurgences of COVID-19 in the coming years can be expected, while public health interventions may be able to reduce their intensity. We used a transmission model to assess how the use of booster doses and non-pharmaceutical interventions (NPIs) amid ongoing pathogen evolution might influence future transmission waves. We find that incidence is likely to increase as NPIs relax, with a second seasonally driven surge expected in autumn 2022. However, booster doses can greatly reduce the intensity of both waves and reduce cumulative deaths by 20% between 7 January 2022 and 7 January 2023. Reintroducing NPIs during the autumn as incidence begins to increase again could also be impactful. Combining boosters and NPIs results in a 30% decrease in cumulative deaths, with potential for greater impacts if variant-adapted boosters are used. Reintroducing these NPIs in autumn 2022 as transmission rates increase provides similar benefits to sustaining NPIs indefinitely (307 000 deaths with indefinite NPIs and boosters compared with 304 000 deaths with transient NPIs and boosters). If novel variants with increased transmissibility or immune escape emerge, deaths will be higher, but vaccination and NPIs are expected to remain effective tools to decrease both cumulative and peak health system burden, providing proportionally similar relative impacts.

Epidemiological models are commonly fit to case and pathogen sequence data to estimate parameters and to infer unobserved disease dynamics. Here, we present an inference approach based on sequence data that is well suited for model fitting early on during the expansion of a viral lineage. Our approach relies on a trajectory of segregating sites to infer epidemiological parameters within a Sequential Monte Carlo framework. Using simulated data, we first show that our approach accurately recovers key epidemiological quantities under a single-introduction scenario. We then apply our approach to SARS-CoV-2 sequence data from France, estimating a basic reproduction number of approximately 2.3-2.7 under an epidemiological model that allows for multiple introductions. Our approach presented here indicates that inference approaches that rely on simple population genetic summary statistics can be informative of epidemiological parameters and can be used for reconstructing infectious disease dynamics during the early expansion of a viral lineage.

Segmentation of viral genomes provides the potential for genetic exchange within coinfected cells. However, for this potential to be realized, coinfecting genomes must mix during the viral life cycle. The efficiency of reassortment, in turn, dictates its potential to drive evolution. The opportunity for mixing within coinfected cells may vary greatly across virus families, such that the evolutionary implications of genome segmentation differ as a result of core features of the viral life cycle. To investigate the relationship between viral replication compartments and genetic exchange, we quantified reassortment in mammalian orthoreovirus (reovirus). Reoviruses carry a 10-segmented, double-stranded RNA genome, which is replicated within proteinaceous structures termed inclusion bodies. We hypothesized that inclusions impose a barrier to reassortment. We quantified reassortment between wild-type (wt) and variant (var) reoviruses that differ by one nucleotide per segment. Studies of wt/var systems in both T1L and T3D backgrounds revealed frequent reassortment without bias toward particular genotypes. However, reassortment was more efficient in the T3D serotype. Since T1L and T3D viruses exhibit different inclusion body morphologies, we tested the impact of this phenotype on reassortment. In both serotypes, reassortment levels did not differ by inclusion morphology. Reasoning that the merging of viral inclusions may be critical for genome mixing, we then tested the effect of blocking merging. Reassortment proceeded efficiently even under these conditions. Our findings indicate that reovirus reassortment is highly efficient despite the localization of many viral processes to inclusion bodies, and that the robustness of this genetic exchange is independent of inclusion body structure and fusion. IMPORTANCE Quantification of reassortment in diverse viral systems is critical to elucidate the implications of genome segmentation for viral evolution. In principle, genome segmentation offers a facile means of genetic exchange between coinfecting viruses. In practice, there may be physical barriers within the cell that limit the mixing of viral genomes. Here, we tested the hypothesis that localization of the various stages of the mammalian orthoreovirus life cycle within cytoplasmic inclusion bodies compartmentalizes viral replication and limits genetic exchange. Contrary to this hypothesis, our data indicate that reovirus reassortment occurs readily within coinfected cells and is not strongly affected by the structure or dynamics of viral inclusion bodies. We conclude that the potential for reassortment to contribute to reovirus evolution is high.

A reanalysis of SARS-CoV-2 deep sequencing data from donor-recipient pairs indicates that transmission bottlenecks are very narrow (one to three virions).

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.

Animal models are frequently used to characterize the within-host dynamics of emerging zoonotic viruses. More recent studies have also deep-sequenced longitudinal viral samples originating from experimental challenges to gain a better understanding of how these viruses may evolve in vivo and between transmission events. These studies have often identified nucleotide variants that can replicate more efficiently within hosts and also transmit more effectively between hosts. Quantifying the degree to which a mutation impacts viral fitness within a host can improve identification of variants that are of particular epidemiological concern and our ability to anticipate viral adaptation at the population level. While methods have been developed to quantify the fitness effects of mutations using observed changes in allele frequencies over the course of a host’s infection, none of the existing methods account for the possibility of cellular coinfection. Here, we develop mathematical models to project variant allele frequency changes in the context of cellular coinfection and, further, integrate these models with statistical inference approaches to demonstrate how variant fitness can be estimated alongside cellular multiplicity of infection. We apply our approaches to empirical longitudinally sampled H5N1 sequence data from ferrets. Our results indicate that previous studies may have significantly underestimated the within-host fitness advantage of viral variants. These findings underscore the importance of considering the process of cellular coinfection when studying within-host viral evolutionary dynamics.

How a history of influenza virus infections contributes to protection is not fully understood, but such protection might explain the contrasting age distributions of cases of the two lineages of influenza B, B/Victoria and B/Yamagata. Fitting a statistical model to those distributions using surveillance data from New Zealand, we found they could be explained by historical changes in lineage frequencies combined with cross-protection between strains of the same lineage. We found additional protection against B/Yamagata in people for whom it was their first influenza B infection, similar to the immune imprinting observed in influenza A. While the data were not informative about B/Victoria imprinting, B/Yamagata imprinting could explain the fewer B/Yamagata than B/Victoria cases in cohorts born in the 1990s and the bimodal age distribution of B/Yamagata cases. Longitudinal studies can test if these forms of protection inferred from historical data extend to more recent strains and other populations.

Transmission of influenza A viruses (IAV) between hosts is subject to numerous physical and biological barriers that impose genetic bottlenecks, constraining viral diversity and adaptation. The bottlenecks within hosts and their potential impacts on evolutionary pathways taken during infection are poorly understood. To address this, we created highly diverse IAV libraries bearing molecular barcodes on two gene segments, enabling high-resolution tracking and quantification of unique virus lineages within hosts. Here we show that IAV infection in lungs is characterized by multiple within-host bottlenecks that result in “islands” of infection in lung lobes, each with genetically distinct populations. We perform site-specific inoculation of barcoded IAV in the upper respiratory tract of ferrets and track viral diversity as infection spreads to the trachea and lungs. We detect extensive compartmentalization of discrete populations within lung lobes. Bottleneck events and localized replication stochastically sample individual viruses from the upper respiratory tract or the trachea that become the dominant genotype in a particular lobe. These populations are shaped strongly by founder effects, with limited evidence for positive selection. The segregated sites of replication highlight the jackpot-style events that contribute to within-host influenza virus evolution and may account for low rates of intrahost adaptation.

The evolutionary mechanisms by which SARS-CoV-2 viruses adapt to mammalian hosts and, potentially, undergo antigenic evolution depend on the ways genetic variation is generated and selected within and between individual hosts. Using domestic cats as a model, we show that SARS-CoV-2 consensus sequences remain largely unchanged over time within hosts, while dynamic sub-consensus diversity reveals processes of genetic drift and weak purifying selection. We further identify a notable variant at amino acid position 655 in Spike (H655Y), which was previously shown to confer escape from human monoclonal antibodies. This variant arises rapidly and persists at intermediate frequencies in index cats. It also becomes fixed following transmission in two of three pairs. These dynamics suggest this site may be under positive selection in this system and illustrate how a variant can quickly arise and become fixed in parallel across multiple transmission pairs. Transmission of SARS-CoV-2 in cats involved a narrow bottleneck, with new infections founded by fewer than ten viruses. In RNA virus evolution, stochastic processes like narrow transmission bottlenecks and genetic drift typically act to constrain the overall pace of adaptive evolution. Our data suggest that here, positive selection in index cats followed by a narrow transmission bottleneck may have instead accelerated the fixation of S H655Y, a potentially beneficial SARS-CoV-2 variant. Overall, our study suggests species- and context-specific adaptations are likely to continue to emerge. This underscores the importance of continued genomic surveillance for new SARS-CoV-2 variants as well as heightened scrutiny for signatures of SARS-CoV-2 positive selection in humans and mammalian model systems.