Genetic drift and purifying selection shape within-host influenza A virus populations during natural swine infections

Patterns of within-host influenza A virus (IAV) diversity and evolution have been described in natural human infections, but these patterns remain poorly characterized in non-human hosts. Elucidating these dynamics is important to better understand IAV biology and the evolutionary processes that govern spillover into humans. Here, we sampled an IAV outbreak in pigs during a week-long county fair to characterize viral diversity and evolution in this important reservoir host. Nasal wipes were collected on a daily basis from all pigs present at the fair, yielding up to 421 samples per day. Subtyping of PCR-positive samples revealed the co-circulation of H1N1 and H3N2 subtype swine IAVs. PCR-positive samples with robust Ct values were deep-sequenced, yielding 506 sequenced samples from a total of 253 pigs. Based on higher-depth re-sequenced data from a subset of these initially sequenced samples (260 samples from 168 pigs), we characterized patterns of within-host IAV genetic diversity and evolution. We find that IAV genetic diversity in single-subtype infected pigs is low, with the majority of intrahost Single Nucleotide Variants (iSNVs) present at frequencies of <10%. The ratio of the number of nonsynonymous to the number of synonymous iSNVs is significantly lower than under the neutral expectation, indicating that purifying selection shapes patterns of within-host viral diversity in swine. The dynamic turnover of iSNVs and their pronounced frequency changes further indicate that genetic drift also plays an important role in shaping IAV populations within pigs. Taken together, our results highlight similarities in patterns of IAV genetic diversity and evolution between humans and swine, including the role of stochastic processes in shaping within-host IAV dynamics.

Prevention of respiratory virus transmission by resident memory CD8+ T cells

An ideal vaccine both attenuates virus growth and disease in infected individuals and reduces the spread of infections in the population, thereby generating herd immunity. Although this strategy has proved successful by generating humoral immunity to measles, yellow fever and polio, many respiratory viruses evolve to evade pre-existing antibodies1. One approach for improving the breadth of antiviral immunity against escape variants is through the generation of memory T cells in the respiratory tract, which are positioned to respond rapidly to respiratory virus infections2-6. However, it is unknown whether memory T cells alone can effectively surveil the respiratory tract to the extent that they eliminate or greatly reduce viral transmission following exposure of an individual to infection. Here we use a mouse model of natural parainfluenza virus transmission to quantify the extent to which memory CD8+ T cells resident in the respiratory tract can provide herd immunity by reducing both the susceptibility of acquiring infection and the extent of transmission, even in the absence of virus-specific antibodies. We demonstrate that protection by resident memory CD8+ T cells requires the antiviral cytokine interferon-γ (IFNγ) and leads to altered transcriptional programming of epithelial cells within the respiratory tract. These results suggest that tissue-resident CD8+ T cells in the respiratory tract can have important roles in protecting the host against viral disease and limiting viral spread throughout the population.

Within-host evolutionary dynamics and tissue compartmentalization during acute SARS-CoV-2 infection

The global evolution of SARS-CoV-2 depends in part upon the evolutionary dynamics within individual hosts with varying immune histories. To characterize the within-host evolution of acute SARS-CoV-2 infection, we sequenced saliva and nasal samples collected daily from vaccinated and unvaccinated individuals early during infection. We show that longitudinal sampling facilitates high-confidence genetic variant detection and reveals evolutionary dynamics missed by less-frequent sampling strategies. Within-host dynamics in both unvaccinated and vaccinated individuals appeared largely stochastic; however, in rare cases, minor genetic variants emerged to frequencies sufficient for forward transmission. Finally, we detected significant genetic compartmentalization of viral variants between saliva and nasal swab sample sites in many individuals. Altogether, these data provide a high-resolution profile of within-host SARS-CoV-2 evolutionary dynamics.IMPORTANCEWe detail the within-host evolutionary dynamics of SARS-CoV-2 during acute infection in 31 individuals using daily longitudinal sampling. We characterized patterns of mutational accumulation for unvaccinated and vaccinated individuals, and observed that temporal variant dynamics in both groups were largely stochastic. Comparison of paired nasal and saliva samples also revealed significant genetic compartmentalization between tissue environments in multiple individuals. Our results demonstrate how selection, genetic drift, and spatial compartmentalization all play important roles in shaping the within-host evolution of SARS-CoV-2 populations during acute infection.

Transmission Bottleneck Size Estimation from De Novo Viral Genetic Variation

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, these approaches have the potential to substantially underestimate true transmission bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arise de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these 2 respiratory viruses.

Genetic drift and purifying selection shape within-host influenza A virus populations during natural swine infections

Patterns of within-host influenza A virus (IAV) diversity and evolution have been described in natural human infections, but these patterns remain poorly characterized in non-human hosts. Elucidating these dynamics is important to better understand IAV biology and the evolutionary processes that govern spillover into humans. Here, we sampled an IAV outbreak in pigs during a week-long county fair to characterize viral diversity and evolution in this important reservoir host. Nasal wipes were collected on a daily basis from all pigs present at the fair, yielding up to 421 samples per day. Subtyping of PCR-positive samples revealed the co-circulation of H1N1 and H3N2 subtype IAVs. PCR-positive samples with robust Ct values were deep-sequenced, yielding 506 sequenced samples from a total of 253 pigs. Based on higher-depth re-sequenced data from a subset of these initially sequenced samples (260 samples from 168 pigs), we characterized patterns of within-host IAV genetic diversity and evolution. We find that IAV genetic diversity in single-subtype infected pigs is low, with the majority of intra-host single nucleotide variants (iSNVs) present at frequencies of <10%. The ratio of the number of nonsynonymous to the number of synonymous iSNVs is significantly lower than under the neutral expectation, indicating that purifying selection shapes patterns of within-host viral diversity in swine. The dynamic turnover of iSNVs and their pronounced frequency changes further indicate that genetic drift also plays an important role in shaping IAV populations within pigs. Taken together, our results highlight similarities in patterns of IAV genetic diversity and evolution between humans and swine, including the role of stochastic processes in shaping within-host IAV dynamics.

Transmission bottleneck size estimation from de novo viral genetic variation

Sequencing of viral infections has become increasingly common over the last decade. Deep sequencing data in particular have proven useful in characterizing the roles that genetic drift and natural selection play in shaping within-host viral populations. They have also been used to estimate transmission bottleneck sizes from identified donor-recipient pairs. These bottleneck sizes quantify the number of viral particles that establish genetic lineages in the recipient host and are important to estimate due to their impact on viral evolution. Current approaches for estimating bottleneck sizes exclusively consider the subset of viral sites that are observed as polymorphic in the donor individual. However, allele frequencies can change dramatically over the course of an individual's infection, such that sites that are polymorphic in the donor at the time of transmission may not be polymorphic in the donor at the time of sampling and allele frequencies at donor-polymorphic sites may change dramatically over the course of a recipient's infection. Because of this, transmission bottleneck sizes estimated using allele frequencies observed at a donor's polymorphic sites may be considerable underestimates of true bottleneck sizes. Here, we present a new statistical approach for instead estimating bottleneck sizes using patterns of viral genetic variation that arose de novo within a recipient individual. Specifically, our approach makes use of the number of clonal viral variants observed in a transmission pair, defined as the number of viral sites that are monomorphic in both the donor and the recipient but carry different alleles. We first test our approach on a simulated dataset and then apply it to both influenza A virus sequence data and SARS-CoV-2 sequence data from identified transmission pairs. Our results confirm the existence of extremely tight transmission bottlenecks for these two respiratory viruses, using an approach that does not tend to underestimate transmission bottleneck sizes.

Effect of Norovirus Inoculum Dose on Virus Kinetics, Shedding, and Symptoms

The effect of norovirus dose on outcomes such as virus shedding and symptoms after initial infection is not well understood. We performed a secondary analysis of a human challenge study by using Bayesian mixed-effects models. As the dose increased from 4.8 to 4,800 reverse transcription PCR units, the total amount of shed virus in feces increased from 4.5 × 1011 to 3.4 × 1012 genomic equivalent copies; in vomit, virus increased from 6.4 × 105 to 3.0 × 107 genomic equivalent copies. Onset time of viral shedding in feces decreased from 1.4 to 0.8 days, and time of peak viral shedding decreased from 2.3 to 1.5 days. Time to symptom onset decreased from 1.5 to 0.8 days. One type of symptom score increased. An increase in norovirus dose was associated with more rapid shedding and symptom onset and possibly increased severity. However, the effect on virus load and shedding was inconclusive.

Epidemiological inference for emerging viruses using segregating sites

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.

Semi-infectious particles contribute substantially to influenza virus within-host dynamics when infection is dominated by spatial structure

Influenza is an ribonucleic acid virus with a genome that comprises eight segments. Experiments show that the vast majority of virions fail to express one or more gene segments and thus cannot cause a productive infection on their own. These particles, called semi-infectious particles (SIPs), can induce virion production through complementation when multiple SIPs are present in an infected cell. Previous within-host influenza models did not explicitly consider SIPs and largely ignore the potential effects of coinfection during virus infection. Here, we constructed and analyzed two distinct models explicitly keeping track of SIPs and coinfection: one without spatial structure and the other implicitly considering spatial structure. While the model without spatial structure fails to reproduce key aspects of within-host influenza virus dynamics, we found that the model implicitly considering the spatial structure of the infection process makes predictions that are consistent with biological observations, highlighting the crucial role that spatial structure plays during an influenza infection. This model predicts two phases of viral growth prior to the viral peak: a first phase driven by fully infectious particles at the initiation of infection followed by a second phase largely driven by coinfections of fully infectious particles and SIPs. Fitting this model to two sets of data, we show that SIPs can contribute substantially to viral load during infection. Overall, the model provides a new interpretation of the in vivo exponential viral growth observed in experiments and a mechanistic explanation for why the production of large numbers of SIPs does not strongly impede viral growth. Being simple and predictive, our model framework serves as a useful tool to understand coinfection dynamics in spatially structured acute viral infections.

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.

Influenza A virus reassortment in mammals gives rise to genetically distinct within-host subpopulations

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.

The role of booster vaccination and ongoing viral evolution in seasonal circulation of SARS-CoV-2

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.

Heterogeneity in viral populations increases the rate of deleterious mutation accumulation

RNA viruses have high mutation rates, with the majority of mutations being deleterious. We examine patterns of deleterious mutation accumulation over multiple rounds of viral replication, with a focus on how cellular coinfection and heterogeneity in viral output affect these patterns. Specifically, using agent-based intercellular simulations we find, in agreement with previous studies, that coinfection of cells by viruses relaxes the strength of purifying selection, and thereby increases the rate of deleterious mutation accumulation. We further find that cellular heterogeneity in viral output exacerbates the rate of deleterious mutation accumulation, regardless of whether this heterogeneity in viral output is stochastic or is due to variation in cellular multiplicity of infection. These results highlight the need to consider the unique life histories of viruses and their population structure to better understand observed patterns of viral evolution.

Influenza A virus undergoes compartmentalized replication in vivo dominated by stochastic bottlenecks

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.

Transmission of influenza A viruses (IAV) between hosts and replication within host impose genetic bottlenecks, constraining viral diversity and adaptation. Here, Amato et al. perform site-specific inoculation of barcoded IAV of ferrets and track viral diversity as infection spreads to the lower respiratory tract and conclude that narrow population bottlenecks are an important feature of the within-host infection dynamics.

The changing epidemiology of SARS-CoV-2

We have come a long way since the start of the COVID-19 pandemic-from hoarding toilet paper and wiping down groceries to sending our children back to school and vaccinating billions. Over this period, the global community of epidemiologists and evolutionary biologists has also come a long way in understanding the complex and changing dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19. In this Review, we retrace our steps through the questions that this community faced as the pandemic unfolded. We focus on the key roles that mathematical modeling and quantitative analyses of empirical data have played in allowing us to address these questions and ultimately to better understand and control the pandemic.

Unrecognized introductions of SARS-CoV-2 into the US state of Georgia shaped the early epidemic

In early 2020, as diagnostic and surveillance responses for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ramped up, attention focused primarily on returning international travelers. Here, we build on existing studies characterizing early patterns of SARS-CoV-2 spread within the USA by analyzing detailed clinical, molecular, and viral genomic data from the state of Georgia through March 2020. We find evidence for multiple early introductions into Georgia, despite relatively sparse sampling. Most sampled sequences likely stemmed from a single or small number of introductions from Asia three weeks prior to the state's first detected infection. Our analysis of sequences from domestic travelers demonstrates widespread circulation of closely related viruses in multiple US states by the end of March 2020. Our findings indicate that the exclusive focus on identifying SARS-CoV-2 in returning international travelers early in the pandemic may have led to a failure to recognize locally circulating infections for several weeks and point toward a critical need for implementing rapid, broadly targeted surveillance efforts for future pandemics.

Comment on "Genomic epidemiology of superspreading events in Austria reveals mutational dynamics and transmission properties of SARS-CoV-2"

[Figure: see text].

Recombinant SARS-CoV-2 genomes circulated at low levels over the first year of the pandemic

Viral recombination can generate novel genotypes with unique phenotypic characteristics, including transmissibility and virulence. Although the capacity for recombination among betacoronaviruses is well documented, recombination between strains of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has not been characterized in detail. Here, we present a lightweight approach for detecting genomes that are potentially recombinant. This approach relies on identifying the mutations that primarily determine SARS-CoV-2 clade structure and then screening genomes for ones that contain multiple mutational markers from distinct clades. Among the over 537,000 genomes queried that were deposited on GISAID.org prior to 16 February 2021, we detected 1,175 potential recombinant sequences. Using a highly conservative criteria to exclude sequences that may have originated through de novo mutation, we find that at least 30 per cent (n = 358) are likely of recombinant origin. An analysis of deep-sequencing data for these putative recombinants, where available, indicated that the majority are high quality. Additional phylogenetic analysis and the observed co-circulation of predicted parent clades in the geographic regions of exposure further support the feasibility of recombination in this subset of potential recombinants. An analysis of these genomes did not reveal evidence for recombination hotspots in the SARS-CoV-2 genome. While most of the putative recombinant sequences we detected were genetic singletons, a small number of genetically identical or highly similar recombinant sequences were identified in the same geographic region, indicative of locally circulating lineages. Recombinant genomes were also found to have originated from parental lineages with substitutions of concern, including D614G, N501Y, E484K, and L452R. Adjusting for an unequal probability of detecting recombinants derived from different parent clades and for geographic variation in clade abundance, we estimate that at most 0.2-2.5 per cent of circulating viruses in the USA and UK are recombinant. Our identification of a small number of putative recombinants within the first year of SARS-CoV-2 circulation underscores the need to sustain efforts to monitor the emergence of new genotypes generated through recombination.

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.

Lineage-specific protection and immune imprinting shape the age distributions of influenza B cases

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.

Fitness Estimation for Viral Variants in the Context of Cellular Coinfection

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.

Characterizing Norovirus Transmission from Outbreak Data, United States

Norovirus is the leading cause of acute gastroenteritis outbreaks in the United States. We estimated the basic (R₀) and effective (R_e) reproduction numbers for 7,094 norovirus outbreaks reported to the National Outbreak Reporting System (NORS) during 2009–2017 and used regression models to assess whether transmission varied by outbreak setting. The median R₀ was 2.75 (interquartile range [IQR] 2.38–3.65), and median R_e was 1.29 (IQR 1.12–1.74). Long-term care and assisted living facilities had an R₀ of 3.35 (95% CI 3.26–3.45), but R₀ did not differ substantially for outbreaks in other settings, except for outbreaks in schools, colleges, and universities, which had an R₀ of 2.92 (95% CI 2.82–3.03). Seasonally, R₀ was lowest (3.11 [95% CI 2.97–3.25]) in summer and peaked in fall and winter. Overall, we saw little variability in transmission across different outbreaks settings in the United States.

Avian Influenza A Viruses Reassort and Diversify Differently in Mallards and Mammals

Reassortment among co-infecting influenza A viruses (IAVs) is an important source of viral diversity and can facilitate expansion into novel host species. Indeed, reassortment played a key role in the evolution of the last three pandemic IAVs. Observed patterns of reassortment within a coinfected host are likely to be shaped by several factors, including viral load, the extent of viral mixing within the host and the stringency of selection. These factors in turn are expected to vary among the diverse host species that IAV infects. To investigate host differences in IAV reassortment, here we examined reassortment of two distinct avian IAVs within their natural host (mallards) and a mammalian model system (guinea pigs). Animals were co-inoculated with A/wildbird/California/187718-36/2008 (H3N8) and A/mallard/Colorado/P66F1-5/2008 (H4N6) viruses. Longitudinal samples were collected from the cloaca of mallards or the nasal tract of guinea pigs and viral genetic exchange was monitored by genotyping clonal isolates from these samples. Relative to those in guinea pigs, viral populations in mallards showed higher frequencies of reassortant genotypes and were characterized by higher genotype richness and diversity. In line with these observations, analysis of pairwise segment combinations revealed lower linkage disequilibrium in mallards as compared to guinea pigs. No clear longitudinal patterns in richness, diversity or linkage disequilibrium were present in either host. Our results reveal mallards to be a highly permissive host for IAV reassortment and suggest that reduced viral mixing limits avian IAV reassortment in a mammalian host.

Recombinant SARS-CoV-2 genomes are currently circulating at low levels

Viral recombination can generate novel genotypes with unique phenotypic characteristics, including transmissibility and virulence. Although the capacity for recombination among betacoronaviruses is well documented, there is limited evidence of recombination between SARS-CoV-2 strains. By identifying the mutations that primarily determine SARS-CoV-2 clade structure, we developed a lightweight approach for detecting recombinant genomes. Among the over 537,000 genomes queried, we detect 1175 putative recombinants that contain multiple mutational markers from distinct clades. Additional phylogenetic analysis and the observed co-circulation of predicted parent clades in the geographic regions of exposure further support the feasibility of recombination in these detected cases. An analysis of these detected cases did not reveal any evidence for recombination hotspots in the SARS-CoV-2 genome. Although most recombinant genotypes were detected a limited number of times, at least two recombinants are now widely transmitted. Recombinant genomes were also found to contain substitutions of concern for elevated transmissibility and lower vaccine efficacy, including D614G, N501Y, E484K, and L452R. Adjusting for an unequal probability of detecting recombinants derived from different parent clades, and for geographic variation in clade abundance, we estimate that at most 5% of circulating viruses in the USA and UK are recombinant. While the phenotypic characterization of detected recombinants was beyond the scope of our analysis, the identification of transmitted recombinants involving substitutions of concern underscores the need to sustain efforts to monitor the emergence of new genotypes generated through recombination.

The Effect of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Mitigation Strategies on Seasonal Respiratory Viruses: A Tale of 2 Large Metropolitan Centers in the United States

To assess the impact of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic on seasonal respiratory viruses, absolute case counts and viral reproductive rates from 2019-2020 were compared against previous seasons. Our findings suggest that the public health measures implemented to reduce SARS-CoV-2 transmission significantly reduced the transmission of other respiratory viruses.

Reanalysis of deep-sequencing data from Austria points towards a small SARS-COV-2 transmission bottleneck on the order of one to three virions.. [PMID: 0 DOI: 10.1101/2021.02.22.432096]

An early analysis of SARS-CoV-2 deep-sequencing data that combined epidemiological and genetic data to characterize the transmission dynamics of the virus in and beyond Austria concluded that the size of the virus’s transmission bottleneck was large – on the order of 1000 virions. We performed new computational analyses using these deep-sequenced samples from Austria. Our analyses included characterization of transmission bottleneck sizes across a range of variant calling thresholds and examination of patterns of shared low-frequency variants between transmission pairs in cases where de novo genetic variation was present in the recipient. From these analyses, among others, we found that SARS-CoV-2 transmission bottlenecks are instead likely to be very tight, on the order of 1-3 virions. These findings have important consequences for understanding how SARS-CoV-2 evolves between hosts and the processes shaping genetic variation observed at the population level.

Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck

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.

Asymptomatic Infection and Transmission of Pertussis in Households: A Systematic Review

We conducted a systematic review to describe the frequency of mild, atypical, and asymptomatic infection among household contacts of pertussis cases and to explore the published literature for evidence of asymptomatic transmission. We included studies that obtained and tested laboratory specimens from household contacts regardless of symptom presentation and reported the proportion of cases with typical, mild/atypical, or asymptomatic infection. After screening 6789 articles, we included 26 studies. Fourteen studies reported household contacts with mild/atypical pertussis. These comprised up to 46.2% of all contacts tested. Twenty-four studies reported asymptomatic contacts with laboratory-confirmed pertussis, comprising up to 55.6% of those tested. Seven studies presented evidence consistent with asymptomatic pertussis transmission between household contacts. Our results demonstrate a high prevalence of subclinical infection in household contacts of pertussis cases, which may play a substantial role in the ongoing transmission of disease. Our review reveals a gap in our understanding of pertussis transmission.

Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck

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.

484. The Impact of SARS-CoV-2 on Reproduction Rates of Seasonal Influenza and RSV

Background

The COVID-19 pandemic caused by SARS-CoV-2 has precipitated a global health crisis. In an effort to decrease person-to-person transmission, societal-level non-pharmacologic interventions (NPIs) to maintain social distancing have been enacted. As SARS-CoV-2 shares similar routes of transmission with other respiratory viruses, implementation of these NPIs may have decreased transmission for multiple viral pathogens. We compared influenza and respiratory syncytial (RSV) rates in prior seasons to rates during the 2019 - 2020 season at two large academic centers in Atlanta and Boston.

Methods

The clinical records were queried for adults with respiratory virus testing conducted at the Emory Healthcare system and associated clinics in Atlanta and the Mass General Brigham (MGB) Healthcare System in Boston. Total cases for influenza A and B, RSV and SARS-CoV-2 were analyzed for each week of the past 5 seasons (07/01/2015-05/30/2020) for the Atlanta and Boston sites. Systematic changes in viral infection rates were calculated using viral reproduction rates, R(t), between consecutive weeks. R(t) is the ratio of the number of positive cases in one week to the number of positive cases in the previous week. We used statistical bootstrapping to determine whether R(t) for influenza and RSV were lower in 2019–2020 following the introduction of SARS-CoV-2. Analyses were conducted using R (v 4.0.0).

Absolute respiratory virus activity by season, Boston (panel A) v. Atlanta (panel B)

Results

For the 2019–2020 Atlanta season, R(t) < 1 (which reflects steady decline in infection rates) occurred at week 28 for influenza A, week 33 for influenza B, and week 35 for RSV, which corresponded with the increase of SARS-Cov-2 cases. The R(t) of these viruses stayed at or near 1 during weeks 33–35 in prior seasons, and R(t) was greater than 1 up to week 47. Data from MGB sites showed similar trends with a sudden decline in R(t) to < 1 at the start of the SARS-CoV-2 pandemic.

Conclusion

We note decreased transmission of influenza and RSV during a time window where widespread movement restrictions and social distancing were imposed to control COVID-19. This trend was most pronounced for influenza A in Atlanta and influenza B in Boston. These data suggest that NPIs can have important effects across multiple pathogens.

Disclosures

Kraft Colleen, MD, MSc, Rebiotix (Advisor or Review Panel member)

Revealing fine-scale spatiotemporal differences in SARS-CoV-2 introduction and spread

Evidence-based public health approaches that minimize the introduction and spread of new SARS-CoV-2 transmission clusters are urgently needed in the United States and other countries struggling with expanding epidemics. Here we analyze 247 full-genome SARS-CoV-2 sequences from two nearby communities in Wisconsin, USA, and find surprisingly distinct patterns of viral spread. Dane County had the 12th known introduction of SARS-CoV-2 in the United States, but this did not lead to descendant community spread. Instead, the Dane County outbreak was seeded by multiple later introductions, followed by limited community spread. In contrast, relatively few introductions in Milwaukee County led to extensive community spread. We present evidence for reduced viral spread in both counties following the statewide "Safer at Home" order, which went into effect 25 March 2020. Our results suggest patterns of SARS-CoV-2 transmission may vary substantially even in nearby communities. Understanding these local patterns will enable better targeting of public health interventions.

Full genome viral sequences inform patterns of SARS-CoV-2 spread into and within Israel

Full genome sequences are increasingly used to track the geographic spread and transmission dynamics of viral pathogens. Here, with a focus on Israel, we sequence 212 SARS-CoV-2 sequences and use them to perform a comprehensive analysis to trace the origins and spread of the virus. We find that travelers returning from the United States of America significantly contributed to viral spread in Israel, more than their proportion in incoming infected travelers. Using phylodynamic analysis, we estimate that the basic reproduction number of the virus was initially around 2.5, dropping by more than two-thirds following the implementation of social distancing measures. We further report high levels of transmission heterogeneity in SARS-CoV-2 spread, with between 2-10% of infected individuals resulting in 80% of secondary infections. Overall, our findings demonstrate the effectiveness of social distancing measures for reducing viral spread.

Full genome viral sequences inform patterns of SARS-CoV-2 spread into and within Israel

Full genome sequences are increasingly used to track the geographic spread and transmission dynamics of viral pathogens. Here, with a focus on Israel, we sequence 212 SARS-CoV-2 sequences and use them to perform a comprehensive analysis to trace the origins and spread of the virus. We find that travelers returning from the United States of America significantly contributed to viral spread in Israel, more than their proportion in incoming infected travelers. Using phylodynamic analysis, we estimate that the basic reproduction number of the virus was initially around 2.5, dropping by more than two-thirds following the implementation of social distancing measures. We further report high levels of transmission heterogeneity in SARS-CoV-2 spread, with between 2-10% of infected individuals resulting in 80% of secondary infections. Overall, our findings demonstrate the effectiveness of social distancing measures for reducing viral spread.

Cellular co-infection can modulate the efficiency of influenza A virus production and shape the interferon response

During viral infection, the numbers of virions infecting individual cells can vary significantly over time and space. The functional consequences of this variation in cellular multiplicity of infection (MOI) remain poorly understood. Here, we rigorously quantify the phenotypic consequences of cellular MOI during influenza A virus (IAV) infection over a single round of replication in terms of cell death rates, viral output kinetics, interferon and antiviral effector gene transcription, and superinfection potential. By statistically fitting mathematical models to our data, we precisely define specific functional forms that quantitatively describe the modulation of these phenotypes by MOI at the single cell level. To determine the generality of these functional forms, we compare two distinct cell lines (MDCK cells and A549 cells), both infected with the H1N1 strain A/Puerto Rico/8/1934 (PR8). We find that a model assuming that infected cell death rates are independent of cellular MOI best fits the experimental data in both cell lines. We further observe that a model in which the rate and efficiency of virus production increase with cellular co-infection best fits our observations in MDCK cells, but not in A549 cells. In A549 cells, we also find that induction of type III interferon, but not type I interferon, is highly dependent on cellular MOI, especially at early timepoints. This finding identifies a role for cellular co-infection in shaping the innate immune response to IAV infection. Finally, we show that higher cellular MOI is associated with more potent superinfection exclusion, thus limiting the total number of virions capable of infecting a cell. Overall, this study suggests that the extent of cellular co-infection by influenza viruses may be a critical determinant of both viral production kinetics and cellular infection outcomes in a host cell type-dependent manner.

During influenza A virus (IAV) infection, the number of virions to enter individual cells can be highly variable. Cellular co-infection appears to be common and plays an essential role in facilitating reassortment for IAV, yet little is known about how cellular co-infection influences infection outcomes at the cellular level. Here, we combine quantitative in vitro infection experiments with statistical model fitting to precisely define the phenotypic consequences of cellular co-infection in two cell lines. We reveal that cellular co-infection can increase and accelerate the efficiency of IAV production in a cell line-dependent fashion, identifying it as a potential determinant of viral replication kinetics. We also show that induction of type III, but not type I, interferon is highly dependent upon the number of virions that infect a given cell, implicating cellular co-infection as an important determinant of the host innate immune response to infection. Altogether, our findings show that cellular co-infection plays a crucial role in determining infection outcome. The integration of experimental and statistical modeling approaches detailed here represents a significant advance in the quantitative study of influenza virus infection and should aid ongoing efforts focused on the construction of mathematical models of IAV infection.

Considering indirect benefits is critical when evaluating SARS-CoV-2 vaccine candidates

Significant progress has already been made in development and testing of SARS-CoV-2 vaccines, and Phase III clinical trials have begun for 6 novel vaccine candidates to date. These Phase III trials seek to demonstrate direct benefits of a vaccine on vaccine recipients. However, vaccination is also known to bring about indirect benefits to a population through the reduction of virus circulation. The indirect effects of SARS-CoV-2 vaccination can play a key role in reducing case counts and COVID-19 deaths. To illustrate this point, we show through simulation that a vaccine with strong indirect effects has the potential to reduce SARS-CoV-2 circulation and COVID-19 deaths to a greater extent than an alternative vaccine with stronger direct effects but weaker indirect effects. Protection via indirect effects may be of particular importance in the context of this virus, because elderly individuals are at an elevated risk of death but are also less likely to be directly protected by vaccination due to immune senescence. We therefore encourage ongoing data collection and model development aimed at evaluating the indirect effects of forthcoming SARS-CoV-2 vaccines.

Distinct patterns of SARS-CoV-2 transmission in two nearby communities in Wisconsin, USA

Evidence-based public health approaches that minimize the introduction and spread of new SARS-CoV-2 transmission clusters are urgently needed in the United States and other countries struggling with expanding epidemics. Here we analyze 247 full-genome SARS-CoV-2 sequences from two nearby communities in Wisconsin, USA, and find surprisingly distinct patterns of viral spread. Dane County had the 12th known introduction of SARS-CoV-2 in the United States, but this did not lead to descendant community spread. Instead, the Dane County outbreak was seeded by multiple later introductions, followed by limited community spread. In contrast, relatively few introductions in Milwaukee County led to extensive community spread. We present evidence for reduced viral spread in both counties, and limited viral transmission between counties, following the statewide Safer-at-Home public health order, which went into effect 25 March 2020. Our results suggest that early containment efforts suppressed the spread of SARS-CoV-2 within Wisconsin.

CpG-creating mutations are costly in many human viruses

Mutations can occur throughout the virus genome and may be beneficial, neutral or deleterious. We are interested in mutations that yield a C next to a G, producing CpG sites. CpG sites are rare in eukaryotic and viral genomes. For the eukaryotes, it is thought that CpG sites are rare because they are prone to mutation when methylated. In viruses, we know less about why CpG sites are rare. A previous study in HIV suggested that CpG-creating transition mutations are more costly than similar non-CpG-creating mutations. To determine if this is the case in other viruses, we analyzed the allele frequencies of CpG-creating and non-CpG-creating mutations across various strains, subtypes, and genes of viruses using existing data obtained from Genbank, HIV Databases, and Virus Pathogen Resource. Our results suggest that CpG sites are indeed costly for most viruses. By understanding the cost of CpG sites, we can obtain further insights into the evolution and adaptation of viruses.

Intrahost cytomegalovirus population genetics following antibody pretreatment in a monkey model of congenital transmission

Human cytomegalovirus (HCMV) infection is the leading non-genetic cause of congenital birth defects worldwide. While several studies have addressed the genetic composition of viral populations in newborns diagnosed with HCMV, little is known regarding mother-to-child viral transmission dynamics and how therapeutic interventions may impact within-host viral populations. Here, we investigate how preexisting CMV-specific antibodies shape the maternal viral population and intrauterine virus transmission. Specifically, we characterize the genetic composition of CMV populations in a monkey model of congenital CMV infection to examine the effects of passively-infused hyperimmune globulin (HIG) on viral population genetics in both maternal and fetal compartments. In this study, 11 seronegative, pregnant monkeys were challenged with rhesus CMV (RhCMV), including a group pretreated with a standard potency HIG preparation (n = 3), a group pretreated with a high-neutralizing potency HIG preparation (n = 3), and an untreated control group (n = 5). Targeted amplicon deep sequencing of RhCMV glycoprotein B and L genes revealed that one of the three strains present in the viral inoculum (UCD52) dominated maternal and fetal viral populations. We identified minor haplotypes of this strain and characterized their dynamics. Many of the identified haplotypes were consistently detected at multiple timepoints within sampled maternal tissues, as well as across tissue compartments, indicating haplotype persistence over time and transmission between maternal compartments. However, haplotype numbers and diversity levels were not appreciably different between control, standard-potency, and high-potency pretreatment groups. We found that while the presence of maternal antibodies reduced viral load and congenital infection, it had no apparent impact on intrahost viral genetic diversity at the investigated loci. Interestingly, some minor haplotypes present in fetal and maternal-fetal interface tissues were also identified as minor haplotypes in corresponding maternal tissues, providing evidence for a loose RhCMV mother-to-fetus transmission bottleneck even in the presence of preexisting antibodies.

Within-host infectious disease models accommodating cellular coinfection, with an application to influenza

Within-host models are useful tools for understanding the processes regulating viral load dynamics. While existing models have considered a wide range of within-host processes, at their core these models have shown remarkable structural similarity. Specifically, the structure of these models generally consider target cells to be either uninfected or infected, with the possibility of accommodating further resolution (e.g. cells that are in an eclipse phase). Recent findings, however, indicate that cellular coinfection is the norm rather than the exception for many viral infectious diseases, and that cells with high multiplicity of infection are present over at least some duration of an infection. The reality of these cellular coinfection dynamics is not accommodated in current within-host models although it may be critical for understanding within-host dynamics. This is particularly the case if multiplicity of infection impacts infected cell phenotypes such as their death rate and their viral production rates. Here, we present a new class of within-host disease models that allow for cellular coinfection in a scalable manner by retaining the low-dimensionality that is a desirable feature of many current within-host models. The models we propose adopt the general structure of epidemiological ‘macroparasite’ models that allow hosts to be variably infected by parasites such as nematodes and host phenotypes to flexibly depend on parasite burden. Specifically, our within-host models consider target cells as ‘hosts’ and viral particles as ‘macroparasites’, and allow viral output and infected cell lifespans, among other phenotypes, to depend on a cell’s multiplicity of infection. We show with an application to influenza that these models can be statistically fit to viral load and other within-host data, and demonstrate using model selection approaches that they have the ability to outperform traditional within-host viral dynamic models. Important in vivo quantities such as the mean multiplicity of cellular infection and time-evolving reassortant frequencies can also be quantified in a straightforward manner once these macroparasite models have been parameterized. The within-host model structure we develop here provides a mathematical way forward to address questions related to the roles of cellular coinfection, collective viral interactions, and viral complementation in within-host viral dynamics and evolution.

Causes and Consequences of Spatial Within-Host Viral Spread

The spread of viral pathogens both between and within hosts is inherently a spatial process. While the spatial aspects of viral spread at the epidemiological level have been increasingly well characterized, the spatial aspects of viral spread within infected hosts are still understudied. Here, with a focus on influenza A viruses (IAVs), we first review experimental studies that have shed light on the mechanisms and spatial dynamics of viral spread within hosts. These studies provide strong empirical evidence for highly localized IAV spread within hosts. Since mathematical and computational within-host models have been increasingly used to gain a quantitative understanding of observed viral dynamic patterns, we then review the (relatively few) computational modeling studies that have shed light on possible factors that structure the dynamics of spatial within-host IAV spread. These factors include the dispersal distance of virions, the localization of the immune response, and heterogeneity in host cell phenotypes across the respiratory tract. While informative, we find in these studies a striking absence of theoretical expectations of how spatial dynamics may impact the dynamics of viral populations. To mitigate this, we turn to the extensive ecological and evolutionary literature on range expansions to provide informed theoretical expectations. We find that factors such as the type of density dependence, the frequency of long-distance dispersal, specific life history characteristics, and the extent of spatial heterogeneity are critical factors affecting the speed of population spread and the genetic composition of spatially expanding populations. For each factor that we identified in the theoretical literature, we draw parallels to its analog in viral populations. We end by discussing current knowledge gaps related to the spatial component of within-host IAV spread and the potential for within-host spatial considerations to inform the development of disease control strategies.

Selection on non-antigenic gene segments of seasonal influenza A virus and its impact on adaptive evolution

Most studies on seasonal influenza A/H3N2 virus adaptation have focused on the main antigenic gene, hemagglutinin. However, there is increasing evidence that the genome-wide genetic background of novel antigenic variants can influence these variants’ emergence probabilities and impact their patterns of dominance in the population. This suggests that non-antigenic genes may be important in shaping the viral evolutionary dynamics. To better understand the role of selection on non-antigenic genes in the adaptive evolution of seasonal influenza viruses, we have developed a simple population genetic model that considers a virus with one antigenic and one non-antigenic gene segment. By simulating this model under different regimes of selection and reassortment, we find that the empirical patterns of lineage turnover for the antigenic and non-antigenic gene segments are best captured when there is both limited viral coinfection and selection operating on both gene segments. In contrast, under a scenario of only neutral evolution in the non-antigenic gene segment, we see persistence of multiple lineages for long periods of time in that segment, which is not compatible with observed molecular evolutionary patterns. Further, we find that reassortment, occurring in coinfected individuals, can increase the speed of viral adaptive evolution by primarily reducing selective interference and genetic linkage effects. Together, these findings suggest that, for influenza, with six internal or non-antigenic gene segments, the evolutionary dynamics of novel antigenic variants are likely to be influenced by the genome-wide genetic background as a result of linked selection among both beneficial and deleterious mutations.

Transmission Bottleneck Size Estimation from Pathogen Deep-Sequencing Data, with an Application to Human Influenza A Virus

The bottleneck governing infectious disease transmission describes the size of the pathogen population transferred from the donor to the recipient host. Accurate quantification of the bottleneck size is particularly important for rapidly evolving pathogens such as influenza virus, as narrow bottlenecks reduce the amount of transferred viral genetic diversity and, thus, may decrease the rate of viral adaptation. Previous studies have estimated bottleneck sizes governing viral transmission by using statistical analyses of variants identified in pathogen sequencing data. These analyses, however, did not account for variant calling thresholds and stochastic viral replication dynamics within recipient hosts. Because these factors can skew bottleneck size estimates, we introduce a new method for inferring bottleneck sizes that accounts for these factors. Through the use of a simulated data set, we first show that our method, based on beta-binomial sampling, accurately recovers transmission bottleneck sizes, whereas other methods fail to do so. We then apply our method to a data set of influenza A virus (IAV) infections for which viral deep-sequencing data from transmission pairs are available. We find that the IAV transmission bottleneck size estimates in this study are highly variable across transmission pairs, while the mean bottleneck size of 196 virions is consistent with a previous estimate for this data set. Furthermore, regression analysis shows a positive association between estimated bottleneck size and donor infection severity, as measured by temperature. These results support findings from experimental transmission studies showing that bottleneck sizes across transmission events can be variable and influenced in part by epidemiological factors.IMPORTANCE The transmission bottleneck size describes the size of the pathogen population transferred from the donor to the recipient host and may affect the rate of pathogen adaptation within host populations. Recent advances in sequencing technology have enabled bottleneck size estimation from pathogen genetic data, although there is not yet a consistency in the statistical methods used. Here, we introduce a new approach to infer the bottleneck size that accounts for variant identification protocols and noise during pathogen replication. We show that failing to account for these factors leads to an underestimation of bottleneck sizes. We apply this method to an existing data set of human influenza virus infections, showing that transmission is governed by a loose, but highly variable, transmission bottleneck whose size is positively associated with the severity of infection of the donor. Beyond advancing our understanding of influenza virus transmission, we hope that this work will provide a standardized statistical approach for bottleneck size estimation for viral pathogens.

Preexisting antibodies can protect against congenital cytomegalovirus infection in monkeys

Human cytomegalovirus (HCMV) is the most common congenital infection and a known cause of microcephaly, sensorineural hearing loss, and cognitive impairment among newborns worldwide. Natural maternal HCMV immunity reduces the incidence of congenital infection, but does not prevent the disease altogether. We employed a nonhuman primate model of congenital CMV infection to investigate the ability of preexisting antibodies to protect against placental CMV transmission in the setting of primary maternal infection and subsequent viremia, which is required for placental virus exposure. Pregnant, CD4+ T cell-depleted, rhesus CMV-seronegative (RhCMV-seronegative) rhesus monkeys were treated with either standardly produced hyperimmune globulin (HIG) from RhCMV-seropositive macaques or dose-optimized, potently RhCMV-neutralizing HIG prior to intravenous challenge with an RhCMV mixture. HIG passive infusion provided complete protection against fetal loss in both groups. The dose-optimized, RhCMV-neutralizing HIG additionally inhibited placental transmission of RhCMV and reduced viral replication and diversity. Our findings suggest that the presence of durable and potently neutralizing antibodies at the time of primary infection can prevent transmission of systemically replicating maternal RhCMV to the developing fetus, and therefore should be a primary target of vaccines to eliminate this neonatal infection.

The effective rate of influenza reassortment is limited during human infection

We characterise the evolutionary dynamics of influenza infection described by viral sequence data collected from two challenge studies conducted in human hosts. Viral sequence data were collected at regular intervals from infected hosts. Changes in the sequence data observed across time show that the within-host evolution of the virus was driven by the reversion of variants acquired during previous passaging of the virus. Treatment of some patients with oseltamivir on the first day of infection did not lead to the emergence of drug resistance variants in patients. Using an evolutionary model, we inferred the effective rate of reassortment between viral segments, measuring the extent to which randomly chosen viruses within the host exchange genetic material. We find strong evidence that the rate of effective reassortment is low, such that genetic associations between polymorphic loci in different segments are preserved during the course of an infection in a manner not compatible with epistasis. Combining our evidence with that of previous studies we suggest that spatial heterogeneity in the viral population may reduce the extent to which reassortment is observed. Our results do not contradict previous findings of high rates of viral reassortment in vitro and in small animal studies, but indicate that in human hosts the effective rate of reassortment may be substantially more limited.

The influenza virus is an important cause of disease in the human population. During the course of an infection the virus can evolve rapidly. An important mechanism of viral evolution is reassortment, whereby different segments of the influenza genome are shuffled with other segments, producing new viral combinations. Here we study natural selection and reassortment during the course of infections occurring in human hosts. Examining viral genome sequence data from these infections, we note that genetic variants that were acquired during the growth of viruses in culture are selected against in the human host. In addition, we find evidence that the effective rate of reassortment is low. We suggest that the spatial separation between viruses in different parts of the host airway may limit the extent to which genetically distinct segments reassort with one another. Within the global population of influenza viruses, reassortment remains an important factor. However, reassortment is not so rapid as to exclude the possibility of interactions between genome segments affecting the course of influenza evolution during a single infection.

The effects of a deleterious mutation load on patterns of influenza A/H3N2's antigenic evolution in humans

Recent phylogenetic analyses indicate that RNA virus populations carry a significant deleterious mutation load. This mutation load has the potential to shape patterns of adaptive evolution via genetic linkage to beneficial mutations. Here, we examine the effect of deleterious mutations on patterns of influenza A subtype H3N2's antigenic evolution in humans. By first analyzing simple models of influenza that incorporate a mutation load, we show that deleterious mutations, as expected, act to slow the virus's rate of antigenic evolution, while making it more punctuated in nature. These models further predict three distinct molecular pathways by which antigenic cluster transitions occur, and we find phylogenetic patterns consistent with each of these pathways in influenza virus sequences. Simulations of a more complex phylodynamic model further indicate that antigenic mutations act in concert with deleterious mutations to reproduce influenza's spindly hemagglutinin phylogeny, co-circulation of antigenic variants, and high annual attack rates.

DOI:http://dx.doi.org/10.7554/eLife.07361.001

Each year, up to 15% of the world's population experience symptoms of an influenza infection, also commonly known as flu. The most common culprit is a strain of the virus called influenza type A subtype H3N2. One reason that so many people become infected each year is that this virus evolves rapidly. Within a few years, proteins on the surface of the virus known as antigens become less recognizable to the immune system of a person who has been previously infected. This means that the person can become ill with the virus again because their immune system cannot mount an effective response to the evolved virus strain.

Influenza virus strains evolve rapidly because their genetic material accumulates mutations quickly. Although some of these mutations are beneficial to the virus, other mutations are harmful and reduce the ability of the virus to spread. Sometimes beneficial mutations may occur alongside harmful ones, but it is not known how the harmful mutations affect the evolution of the virus.

Here, Koelle and Rasmussen used computer models of H3N2 influenza to examine the effect of harmful mutations on the evolution of this virus population. The models show that harmful mutations limit how quickly the antigens can evolve. Also, the presence of these harmful mutations effectively acts as a sieve: they allow only large changes in the antigens to establish in the virus population.

The models suggest that there are three routes by which large changes in the antigens on H3N2 viruses may occur. The first is by a single mutation that has a big effect on the antigens in viruses that only carry a few harmful mutations, but these large mutations would not happen very often. Another route may be through more common mutations that have only a small or moderate benefit, which would allow the virus to become more common in the population before it acquires a beneficial mutation with a much greater effect. The third possibility is that a large beneficial mutation may arise in viruses that have many harmful mutations. These harmful mutations may initially limit the ability of the virus to spread, but over time, some of these harmful mutations may then be lost.

Koelle and Rasmussen found that the computer models could recreate the patterns of virus evolution that have been observed in real strains of H3N2. Researchers use predictions of influenza evolution to help them decide which virus strains should be included in flu vaccines each year. Koelle and Rasmussen findings indicate that harmful mutations should be considered when making these predictions.

DOI:http://dx.doi.org/10.7554/eLife.07361.002

Minimal within-host dengue models highlight the specific roles of the immune response in primary and secondary dengue infections

In recent years, the within-host viral dynamics of dengue infections have been increasingly characterized, and the relationship between aspects of these dynamics and the manifestation of severe disease has been increasingly probed. Despite this progress, there are few mathematical models of within-host dengue dynamics, and the ones that exist focus primarily on the general role of immune cells in the clearance of infected cells, while neglecting other components of the immune response in limiting viraemia. Here, by considering a suite of mathematical within-host dengue models of increasing complexity, we aim to isolate the critical components of the innate and the adaptive immune response that suffice in the reproduction of several well-characterized features of primary and secondary dengue infections. By building up from a simple target cell limited model, we show that only the innate immune response is needed to recover the characteristic features of a primary symptomatic dengue infection, while a higher rate of viral infectivity (indicative of antibody-dependent enhancement) and infected cell clearance by T cells are further needed to recover the characteristic features of a secondary dengue infection. We show that these minimal models can reproduce the increased risk of disease associated with secondary heterologous infections that arises as a result of a cytokine storm, and, further, that they are consistent with virological indicators that predict the onset of severe disease, such as the magnitude of peak viraemia, time to peak viral load, and viral clearance rate. Finally, we show that the effectiveness of these virological indicators to predict the onset of severe disease depends on the contribution of T cells in fuelling the cytokine storm.