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

Tracking Seasonal Influenza Trends in South Tyrol During 2022/2023 Using Genomic Surveillance Data

BACKGROUND

Influenza seasons are characterized by a complex interplay of co-circulating strains with high spatial and temporal heterogeneities. Surveillance is crucial for monitoring the spread and evolution of the virus and design effective public health response strategies.

AIM

We combined epidemiological, virological, and genomic surveillance data to provide a comprehensive analysis of influenza subtypes circulating in the South Tyrol region (Italy) during season 2022/2023, leveraging phylogenetic and phylodynamic approaches.

METHODS

Clinical samples were collected from patients exhibiting influenza-like symptoms and evaluated by molecular diagnostics. Whole genome sequencing was conducted, and the hemagglutinin (HA) gene sequences were used for phylogenetic analysis. A birth-death skyline model was applied to estimate strain-specific effective reproduction numbers (Re) and attack rates.

RESULTS

Out of 4891 samples tested, 862 tested positive for influenza, of which 224 genomes were sequenced. Phylogenetic analysis of HA gene revealed A(H3N2) strains predominantly clustering in clade 3C.2a1b.2a.2b, followed by 3C.2a1b.2a.1b. A(H1N1pdm09) strains predominantly clustered in clade 6B.1A.5a.2a. Exclusive circulation of B (Victoria) subtype strains aligned with the global trend, all falling within clade V1A.3a.2. Phylogenetic analyses indicate that the strains isolated in the South Tyrol region closely resembled those circulating in the rest of Italy and Austria. Re peaked at 1.16-1.35 (95%CI) for A(H3N2), 1.06-1.34 for A(H1N1pdm09) and 1.02-1.29 for B (Victoria). 95%CI of attack rates were 6.3%-33.5% for A(H3N2), 0.6%-5.0% for A(H1N1pdm09), and 0.8%-6.5% for B (Victoria).

CONCLUSION

Epidemiological estimates from traditional surveillance data can be corroborated by those derived from genomic sequencing, providing more robust assessments of viral transmissibility and attack rates with limited additional effort.

Bayesian mixture models for phylogenetic source attribution from consensus sequences and time since infection estimates

In stopping the spread of infectious diseases, pathogen genomic data can be used to reconstruct transmission events and characterize population-level sources of infection. Most approaches for identifying transmission pairs do not account for the time passing since the divergence of pathogen variants in individuals, which is problematic in viruses with high within-host evolutionary rates. This prompted us to consider possible transmission pairs in terms of phylogenetic data and additional estimates of time since infection derived from clinical biomarkers. We develop Bayesian mixture models with an evolutionary clock as a signal component and additional mixed effects or covariate random functions describing the mixing weights to classify potential pairs into likely and unlikely transmission pairs. We demonstrate that although sources cannot be identified at the individual level with certainty, even with the additional data on time elapsed, inferences into the population-level sources of transmission are possible, and more accurate than using only phylogenetic data without time since infection estimates. We apply the proposed approach to estimate age-specific sources of HIV infection in Amsterdam tranamission networks among men who have sex with men between 2010 and 2021. This study demonstrates that infection time estimates provide informative data to characterize transmission sources, and shows how phylogenetic source attribution can then be done with multi-dimensional mixture models.

Phylogenetic Study of Local Patterns Influenza A(H3N2) Virus Transmission in a Semi-Isolated Population in a Remote Island in Japan Between 2011 and 2013

BACKGROUND

Influenza A outbreak risk is impacted by the potential for importation and local transmission. Reconstructing transmission history with phylogenetic analysis of genetic sequences can help assess outbreak risk but relies on regular collection of genetic sequences. Few influenza genetic sequences are collected in Japan, which makes phylogenetic analysis challenging, especially in rural, remote settings. We generated influenza A genetic sequences from nasopharyngeal swabs (NPS) samples collected using rapid influenza diagnostic tests and used them to analyze the transmission dynamics of influenza in a remote island in Japan.

METHODS

We generated 229 whole genome sequences of influenza A/H3N2 collected during 2011/12 and 2012/13 influenza seasons in Kamigoto Island, Japan, of which 178 sequences passed the quality check. We built time-resolved phylogenetic trees from hemagglutinin sequences to classify the circulating clades by comparing the Kamigoto sequences to global sequences. Spatiotemporal transmission patterns were then analyzed for the largest local clusters.

RESULTS

Using a time-resolved phylogenetic tree, we showed that the sequences clustered in six independent transmission groups (1 in 2011/12, 5 in 2012/13). Sequences were closely related to strains from mainland Japan. All 2011/12 strains were identified as clade 3C.2 (n = 29), while 2012/13 strains fell into two clades: clade 3C.2 (n = 129) and 3C.3a (n = 20). Clusters reported in 2012/13 circulated simultaneously in the same regions. The spatiotemporal analysis of the largest cluster revealed that while the first sequences were reported in the busiest district of Kamigoto, the later sequences were scattered across the island.

CONCLUSION

Kamigoto Island was exposed to repeated importations of Influenza A(H3N2), mostly from mainland Japan, sometimes leading to local transmission and ultimately outbreaks. As independent groups of sequences overlapped in time and space, cases may be wrongly allocated to the same transmission group in the absence of genomic surveillance, thereby underestimating the risk of importations. Our analysis highlights how NPS could be used to better understand influenza transmission patterns in little-studied settings and improve influenza surveillance in Japan.

Local patterns of spread of influenza A H3N2 virus in coastal Kenya over a 1-year period revealed through virus sequence data

The patterns of spread of influenza A viruses in local populations in tropical and sub-tropical regions are unclear due to sparsity of representative spatiotemporal sequence data. We sequenced and analyzed 58 influenza A(H3N2) virus genomes sampled between December 2015 and December 2016 from nine health facilities within the Kilifi Health and Demographic Surveillance System (KHDSS), a predominantly rural region, covering approximately 891 km2 along the Kenyan coastline. The genomes were compared with 1571 contemporaneous global sequences from 75 countries. We observed at least five independent introductions of A(H3N2) viruses into the region during the one-year period, with the importations originating from Africa, Europe, and North America. We also inferred 23 virus location transition events between the nine facilities included in the study. International virus imports into the study area were captured at the facilities of Chasimba, Matsangoni, Mtondia, and Mavueni, while all four exports from the region were captured from the Chasimba facility, all occurring to Africa destinations. A strong spatial clustering of virus strains at all locations was observed associated with local evolution. Our study shows that influenza A(H3N2) virus epidemics in local populations appear to be characterized by limited introductions followed by significant local spread and evolution. Knowledge of the viral lineages that circulate within specific populations in understudied tropical and subtropical regions is required to understand the full diversity and global ecology of influenza viruses and to inform vaccination strategies within these populations.

An Intense Out-of-Season Rebound of Influenza Activity After the Relaxation of Coronavirus Disease 2019 Restrictions in Beijing, China

BACKGROUND

The aim of this study was to investigate the changes of epidemic characteristics of influenza activity pre- and post-coronavirus disease 2019 (COVID-19) in Beijing, China.

METHODS

Epidemiologic data were collected from the influenza surveillance system in Beijing. We compared epidemic intensity, epidemic onset and duration, and influenza transmissibility during the 2022-2023 season with pre-COVID-19 seasons from 2014 to 2020.

RESULTS

The overall incidence rate of influenza in the 2022-2023 season was significantly higher than that of the pre-COVID-19 period, with the record-high level of epidemic intensity in Beijing. The onset and duration of the influenza epidemic period in 2022-2023 season was notably later and shorter than that of the 2014-2020 seasons. Maximum daily instantaneous reproduction number (Rt) of the 2022-2023 season (Rt = 2.31) was much higher than that of the pre-COVID-19 period (Rt = 1.49). The incidence of influenza A(H1N1) and A(H3N2) were the highest among children aged 0-4 years and 5-14 years, respectively, in the 2022-2023 season.

CONCLUSIONS

A late, intense, and short-term peak influenza activity was observed in the 2022-2023 season in Beijing. Children <15 years old were impacted the most by the interruption of influenza circulation during the COVID-19 pandemic. Maintaining continuous surveillance and developing targeted public health strategies of influenza is necessary.

Application of phylodynamics to identify spread of antimicrobial-resistant Escherichia coli between humans and canines in an urban environment

The transmission of antimicrobial resistant bacteria in the urban environment is poorly understood. We utilized genomic sequencing and phylogenetics to characterize the transmission dynamics of antimicrobial resistant Escherichia coli (AMR-Ec) cultured from putative canine (caninep) and human feces present on urban sidewalks in San Francisco, California. We isolated a total of fifty-six AMR-Ec isolates from human (n = 20) and caninep (n = 36) fecal samples from the Tenderloin and South of Market (SoMa) neighborhoods of San Francisco. We then analyzed phenotypic and genotypic antimicrobial resistance (AMR) of the isolates, as well as clonal relationships based on cgMLST and single nucleotide polymorphisms (SNPs) of the core genomes. Using Bayesian inference, we reconstructed the transmission dynamics between humans and caninesp from multiple local outbreak clusters using the marginal structured coalescent approximation (MASCOT). Our results provide evidence for multiple sharing events of AMR-Ec between humans and caninesp. In particular, we found one instance of likely transmission from caninesp to humans as well as an additional local outbreak cluster consisting of one caninep and one human sample. Based on this analysis, it appears that non-human feces act as an important reservoir of clinically relevant AMR-Ec within the urban environment for this study population. This work showcases the utility of genomic epidemiology to reconstruct potential pathways by which antimicrobial resistance spreads.

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.

Reconstructing relative transmission rates in Bayesian phylodynamics: Two-fold transmission advantage of Omicron in Berlin, Germany during December 2021

Phylodynamic methods have lately played a key role in understanding the spread of infectious diseases. During the coronavirus disease (COVID-19) pandemic, large scale genomic surveillance has further increased the potential of dynamic inference from viral genomes. With the continual emergence of novel severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) variants, explicitly allowing transmission rate differences between simultaneously circulating variants in phylodynamic inference is crucial. In this study, we present and empirically validate an extension to the BEAST2 package birth-death skyline model (BDSKY), BDSKY[Formula: see text], which introduces a scaling factor for the transmission rate between independent, jointly inferred trees. In an extensive simulation study, we show that BDSKY[Formula: see text] robustly infers the relative transmission rates under different epidemic scenarios. Using publicly available genome data of SARS-CoV-2, we apply BDSKY[Formula: see text] to quantify the transmission advantage of the Omicron over the Delta variant in Berlin, Germany. We find the overall transmission rate of Omicron to be scaled by a factor of two with pronounced variation between the individual clusters of each variant. These results quantify the transmission advantage of Omicron over the previously circulating Delta variant, in a crucial period of pre-established non-pharmaceutical interventions. By inferring variant- as well as cluster-specific transmission rate scaling factors, we show the differences in transmission dynamics for each variant. This highlights the importance of incorporating lineage-specific transmission differences in phylodynamic inference.

Phylodynamics Uncovers the Transmission of Antibiotic-Resistant Escherichia coli between Canines and Humans in an Urban Environment

The role of canines in transmitting antibiotic resistant bacteria to humans in the urban environment is poorly understood. To elucidate this role, we utilized genomic sequencing and phylogenetics to characterize the burden and transmission dynamics of antibiotic resistant Escherichia coli (ABR-Ec) cultured from canine and human feces present on urban sidewalks in San Francisco, California. We collected a total of fifty-nine ABR-Ec from human (n=12) and canine (n=47) fecal samples from the Tenderloin and South of Market (SoMa) neighborhoods of San Francisco. We then analyzed phenotypic and genotypic antibiotic resistance (ABR) of the isolates, as well as clonal relationships based on cgMLST and single nucleotide polymorphisms (SNPs) of the core genomes. Using Bayesian inference, we reconstructed the transmission dynamics between humans and canines from multiple local outbreak clusters using the marginal structured coalescent approximation (MASCOT). Overall, we found human and canine samples to carry similar amounts and profiles of ABR genes. Our results provide evidence for multiple transmission events of ABR-Ec between humans and canines. In particular, we found one instance of likely transmission from canines to humans as well as an additional local outbreak cluster consisting of one canine and one human sample. Based on this analysis, it appears that canine feces act as an important reservoir of clinically relevant ABR-Ec within the urban environment. Our findings support that public health measures should continue to emphasize proper canine feces disposal practices, access to public toilets and sidewalk and street cleaning. Importance: Antibiotic resistance in E. coli is a growing public health concern with global attributable deaths projected to reach millions annually. Current research has focused heavily on clinical routes of antibiotic resistance transmission to design interventions while the role of alternative reservoirs such as domesticated animals remain less well understood. Our results suggest canines are part of the transmission network that disseminates high-risk multidrug resistance in E. coli within the urban San Francisco community. As such, this study highlights the need to consider canines, and potentially domesticated animals more broadly, when designing interventions to reduce the prevalence of antibiotic resistance in the community. Additionally, it showcases the utility of genomic epidemiology to reconstruct the pathways by which antimicrobial resistance spreads.

Transmission dynamics of seasonal influenza in a remote island population

Seasonal influenza outbreaks remain an important public health concern, causing large numbers of hospitalizations and deaths among high-risk groups. Understanding the dynamics of individual transmission is crucial to design effective control measures and ultimately reduce the burden caused by influenza outbreaks. In this study, we analyzed surveillance data from Kamigoto Island, Japan, a semi-isolated island population, to identify the drivers of influenza transmission during outbreaks. We used rapid influenza diagnostic test (RDT)-confirmed surveillance data from Kamigoto island, Japan and estimated age-specific influenza relative illness ratios (RIRs) over eight epidemic seasons (2010/11 to 2017/18). We reconstructed the probabilistic transmission trees (i.e., a network of who-infected-whom) using Bayesian inference with Markov-chain Monte Carlo method and then performed a negative binomial regression on the inferred transmission trees to identify the factors associated with onwards transmission risk. Pre-school and school-aged children were most at risk of getting infected with influenza, with RIRs values consistently above one. The maximal RIR values were 5.99 (95% CI 5.23, 6.78) in the 7-12 aged-group and 5.68 (95%CI 4.59, 6.99) in the 4-6 aged-group in 2011/12. The transmission tree reconstruction suggested that the number of imported cases were consistently higher in the most populated and busy districts (Tainoura-go and Arikawa-go) ranged from 10-20 to 30-36 imported cases per season. The number of secondary cases generated by each case were also higher in these districts, which had the highest individual reproduction number (R_eff: 1.2-1.7) across the seasons. Across all inferred transmission trees, the regression analysis showed that cases reported in districts with lower local vaccination coverage (incidence rate ratio IRR = 1.45 (95% CI 1.02, 2.05)) or higher number of inhabitants (IRR = 2.00 (95% CI 1.89, 2.12)) caused more secondary transmissions. Being younger than 18 years old (IRR = 1.38 (95%CI 1.21, 1.57) among 4-6 years old and 1.45 (95% CI 1.33, 1.59) 7-12 years old) and infection with influenza type A (type B IRR = 0.83 (95% CI 0.77, 0.90)) were also associated with higher numbers of onwards transmissions. However, conditional on being infected, we did not find any association between individual vaccination status and onwards transmissibility. Our study showed the importance of focusing public health efforts on achieving high vaccine coverage throughout the island, especially in more populated districts. The strong association between local vaccine coverage (including neighboring regions), and the risk of transmission indicate the importance of achieving homogeneously high vaccine coverage. The individual vaccine status may not prevent onwards transmission, though it may reduce the severity of infection.

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.

Cov2clusters: genomic clustering of SARS-CoV-2 sequences

Background

The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada.

Results

We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data.

Conclusions

Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08936-4.

Characterizing the Countrywide Epidemic Spread of Influenza A(H1N1)pdm09 Virus in Kenya between 2009 and 2018

The spatiotemporal patterns of spread of influenza A(H1N1)pdm09 viruses on a countrywide scale are unclear in many tropical/subtropical regions mainly because spatiotemporally representative sequence data are lacking. We isolated, sequenced, and analyzed 383 A(H1N1)pdm09 viral genomes from hospitalized patients between 2009 and 2018 from seven locations across Kenya. Using these genomes and contemporaneously sampled global sequences, we characterized the spread of the virus in Kenya over several seasons using phylodynamic methods. The transmission dynamics of A(H1N1)pdm09 virus in Kenya were characterized by (i) multiple virus introductions into Kenya over the study period, although only a few of those introductions instigated local seasonal epidemics that then established local transmission clusters, (ii) persistence of transmission clusters over several epidemic seasons across the country, (iii) seasonal fluctuations in effective reproduction number (Re) associated with lower number of infections and seasonal fluctuations in relative genetic diversity after an initial rapid increase during the early pandemic phase, which broadly corresponded to epidemic peaks in the northern and southern hemispheres, (iv) high virus genetic diversity with greater frequency of seasonal fluctuations in 2009-2011 and 2018 and low virus genetic diversity with relatively weaker seasonal fluctuations in 2012-2017, and (v) virus spread across Kenya. Considerable influenza virus diversity circulated within Kenya, including persistent viral lineages that were unique to the country, which may have been capable of dissemination to other continents through a globally migrating virus population. Further knowledge of the viral lineages that circulate within understudied low-to-middle-income tropical and subtropical regions is required to understand the full diversity and global ecology of influenza viruses in humans and to inform vaccination strategies within these regions.

Fundamental Identifiability Limits in Molecular Epidemiology

Viral phylogenies provide crucial information on the spread of infectious diseases, and many studies fit mathematical models to phylogenetic data to estimate epidemiological parameters such as the effective reproduction ratio (Re) over time. Such phylodynamic inferences often complement or even substitute for conventional surveillance data, particularly when sampling is poor or delayed. It remains generally unknown, however, how robust phylodynamic epidemiological inferences are, especially when there is uncertainty regarding pathogen prevalence and sampling intensity. Here, we use recently developed mathematical techniques to fully characterize the information that can possibly be extracted from serially collected viral phylogenetic data, in the context of the commonly used birth-death-sampling model. We show that for any candidate epidemiological scenario, there exists a myriad of alternative, markedly different, and yet plausible "congruent" scenarios that cannot be distinguished using phylogenetic data alone, no matter how large the data set. In the absence of strong constraints or rate priors across the entire study period, neither maximum-likelihood fitting nor Bayesian inference can reliably reconstruct the true epidemiological dynamics from phylogenetic data alone; rather, estimators can only converge to the "congruence class" of the true dynamics. We propose concrete and feasible strategies for making more robust epidemiological inferences from viral phylogenetic data.

Transmissibility and transmission of respiratory viruses

Human respiratory virus infections lead to a spectrum of respiratory symptoms and disease severity, contributing to substantial morbidity, mortality and economic losses worldwide, as seen in the COVID-19 pandemic. Belonging to diverse families, respiratory viruses differ in how easy they spread (transmissibility) and the mechanism (modes) of transmission. Transmissibility as estimated by the basic reproduction number (R0) or secondary attack rate is heterogeneous for the same virus. Respiratory viruses can be transmitted via four major modes of transmission: direct (physical) contact, indirect contact (fomite), (large) droplets and (fine) aerosols. We know little about the relative contribution of each mode to the transmission of a particular virus in different settings, and how its variation affects transmissibility and transmission dynamics. Discussion on the particle size threshold between droplets and aerosols and the importance of aerosol transmission for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza virus is ongoing. Mechanistic evidence supports the efficacies of non-pharmaceutical interventions with regard to virus reduction; however, more data are needed on their effectiveness in reducing transmission. Understanding the relative contribution of different modes to transmission is crucial to inform the effectiveness of non-pharmaceutical interventions in the population. Intervening against multiple modes of transmission should be more effective than acting on a single mode.

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.