Phylodynamic Methods Research Articles

The resurgence of influenza A/H3N2 virus in Australia after the relaxation of COVID-19 restrictions during the 2022 season.

This study retrospectively analyzed the genetic characteristics of influenza A H3N2 (A/H3N2) viruses circulating in New South Wales (NSW), the Australian state with the highest number of influenza cases in 2022, and explored the phylodynamics of A/H3N2 transmission within Australia during this period. Sequencing was performed on 217 archived specimens, and A/H3N2 evolution and spread within Australia were analyzed using phylogenetic and phylodynamic methods. Hemagglutinin genes of all analyzed NSW viruses belonged to subclade 3C.2a1b.2a.2 and clustered together with the 2022 vaccine strain. Complete genome analysis of NSW viruses revealed highly frequent interclade reassortments between subclades 3C.2a1b.2a.2 and 3C.2a1b.1a. The estimated earliest introduction time of the dominant subgroup 3C.2a1b.2a.2a.1 in Australia was February 22, 2022 (95% highest posterior density: December 19, 2021-March 13, 2022), following the easing of Australian travel restrictions, suggesting a possible international source. Phylogeographic analysis revealed that Victoria drove the transmission of A/H3N2 viruses across the country during this season, while NSW did not have a dominant role in viral dissemination to other regions. This study highlights the importance of continuous surveillance and genomic characterization of influenza viruses in the postpandemic era, which can inform public health decision-making and enable early detection of novel strains with pandemic potential.

Ecological Changes Exacerbating the Spread of Invasive Ticks has Driven the Dispersal of Severe Fever with Thrombocytopenia Syndrome Virus Throughout Southeast Asia.

Severe fever with thrombocytopenia syndrome virus (SFTSV) is a tick-borne virus recognized by the World Health Organization as an emerging infectious disease of growing concern. Utilizing phylodynamic and phylogeographic methods, we have reconstructed the origin and transmission patterns of SFTSV lineages and the roles demographic, ecological, and climatic factors have played in shaping its emergence and spread throughout Asia. Environmental changes and fluctuations in tick populations, exacerbated by the widespread use of pesticides, have contributed significantly to its geographic expansion. The increased adaptability of Lineage L2 strains to the Haemaphysalis longicornis vector has facilitated the dispersal of SFTSV through Southeast Asia. Increased surveillance and proactive measures are needed to prevent further spread to Australia, Indonesia, and North America.

Inferring enterovirus D68 transmission dynamics from the genomic data of two 2022 North American outbreaks

Enterovirus D68 (EV-D68) has emerged as a significant cause of acute respiratory illness in children globally, notably following its extensive outbreak in North America in 2014. A recent outbreak of EV-D68 was observed in Ontario, Canada, from August to October 2022. Our phylogenetic analysis revealed a notable genetic similarity between the Ontario outbreak and a concurrent outbreak in Maryland, USA. Utilizing Bayesian phylodynamic modeling on whole genome sequences (WGS) from both outbreaks, we determined the median peak time-varying reproduction number (Rt) to be 2.70, 95% HPD (1.76, 4.08) in Ontario and 2.10, 95% HPD (1.41, 3.17) in Maryland. The Rt trends in Ontario closely matched those derived via EpiEstim using reported case numbers. Our study also provides new insights into the median infection duration of EV-D68, estimated at 7.94 days, 95% HPD (4.55, 12.8) in Ontario and 10.8 days, 95% HPD (5.85, 18.6) in Maryland, addressing the gap in the existing literature surrounding EV-D68’s infection period. We observed that the estimated Time since the Most Recent Common Ancestor (TMRCA) and the epidemic’s origin coincided with the easing of COVID-19 related social contact restrictions in both areas. This suggests that the relaxation of non-pharmaceutical interventions, initially implemented to control COVID-19, may have inadvertently facilitated the spread of EV-D68. These findings underscore the effectiveness of phylodynamic methods in public health, demonstrating their broad application from local to global scales and underscoring the critical role of pathogen genomic data in enhancing public health surveillance and outbreak characterization.

Inferring the evolutionary history of the Sino-Himalayan biodiversity hotspot using a Bayesian birth–death skyline model

The current status of the Sino-Himalayan region as a biodiversity hotspot, particularly for flora, has often been linked to the uplift of the Sino-Tibetan Plateau and Himalayan and Hengduan Mountains. However, the relationship between the topological development of the region and the onset of diversification is yet to be confirmed. Here, we apply Bayesian phylodynamic methods to a large phylogeny of angiosperm species from the Sino-Himalayas to infer changes in their rates of evolution through time. We find strong evidence for high diversification rates in the Paleocene, late Miocene and Pliocene, and for negative diversification rates in the Quaternary driven by an increase in extinction rates. Our analyses suggest that changes in global palaeotemperatures are unlikely to be a driving force for these rate shifts. Instead, the collision of the Indian continent with Eurasia and coeval topographic change in the Sino-Himalayas, the Miocene Grassland Expansion and the impact of Pleistocene glaciations on this altitudinally variable region may have driven these rates. We also demonstrate the strong influence of the choice of change times on the shape of inferred piecewise-constant trajectories in Bayesian phylodynamics, and advocate for the use of prior information when making this decision.

Characterising HIV-1 transmission in Victoria, Australia: a molecular epidemiological study

In Australia the incidence of HIV has declined steadily, yet sustained reduction of HIV transmission in this setting requires improved public health responses. As enhanced public health responses and prioritisation of resources may be guided by molecular epidemiological data, here we aimed to assess the applicability of these approaches in Victoria, Australia. A comprehensive collection of HIV-1 pol sequences from individuals diagnosed with HIV in Victoria, Australia, between January 1st 2000 and December 31st 2020 were deidentified and used as the basis of our assessment. These sequences were subtyped and surveillance drug resistance mutations (SDRMs) identified, before definition of transmission groups was performed using HIV-TRACE (0.4.4). Phylodynamic methods were applied using BEAST (2.6.6), assessing effective reproductive numbers for large groups, and additional demographic data were integrated to provide a high resolution view of HIV transmission in Victoria on a decadal time scale. Based on standard settings for HIV-TRACE, 70% (2438/3507) of analysed HIV-1 pol sequences were readily assigned to a transmission group. Individuals in transmission groups were more commonly males (aOR 1.50), those born in Australia (aOR 2.13), those with probable place of acquisition as Victoria (aOR 6.73), and/or those reporting injectable drug use (aOR 2.13). SDRMs were identified in 375 patients (10.7%), with sustained transmission of these limited to a subset of smaller groups. Informative patterns of epidemic growth, stabilisation, and decline were observed; many transmission groups showed effective reproductive numbers (R e ) values reaching greater than 4.0, representing considerable epidemic growth, while others maintained low R e values. This study provides a high resolution view of HIV transmission in Victoria, Australia, and highlights the potential of molecular epidemiology to guide and enhance public health responses in this setting. This informs ongoing discussions with community groups on the acceptability and place of molecular epidemiological approaches in Australia. National Health and Medical Research Council, Australian Research Council.

The emergence of highly resistant and hypervirulent Klebsiella pneumoniae CC14 clone in a tertiary hospital over 8 years

BackgroundKlebsiella pneumoniae is a major bacterial and opportunistic human pathogen, increasingly recognized as a healthcare burden globally. The convergence of resistance and virulence in K. pneumoniae strains has led to the formation of hypervirulent and multidrug-resistant strains with dual risk, limiting treatment options. K. pneumoniae clones are known to emerge locally and spread globally. Therefore, an understanding of the dynamics and evolution of the emerging strains in hospitals is warranted to prevent future outbreaks.MethodsIn this study, we conducted an in-depth genomic analysis on a large-scale collection of 328 multidrug-resistant (MDR) K. pneumoniae strains recovered from 239 patients from a single major hospital in the western coastal city of Jeddah in Saudi Arabia from 2014 through 2022. We employed a broad range of phylogenetic and phylodynamic methods to understand the evolution of the predominant clones on epidemiological time scales, virulence and resistance determinants, and their dynamics. We also integrated the genomic data with detailed electronic health record (EHR) data for the patients to understand the clinical implications of the resistance and virulence of different strains.ResultsWe discovered a diverse population underlying the infections, with most strains belonging to Clonal Complex 14 (CC14) exhibiting dominance. Specifically, we observed the emergence and continuous expansion of strains belonging to the dominant ST2096 in the CC14 clade across hospital wards in recent years. These strains acquired resistance mutations against colistin and extended spectrum β-lactamase (ESBL) and carbapenemase genes, namely blaOXA-48 and blaOXA-232, located on three distinct plasmids, on epidemiological time scales. Strains of ST2096 exhibited a high virulence level with the presence of the siderophore aerobactin (iuc) locus situated on the same mosaic plasmid as the ESBL gene. Integration of ST2096 with EHR data confirmed the significant link between colonization by ST2096 and the diagnosis of sepsis and elevated in-hospital mortality (p-value < 0.05).ConclusionsOverall, these results demonstrate the clinical significance of ST2096 clones and illustrate the rapid evolution of an emerging hypervirulent and MDR K. pneumoniae in a clinical setting.

Phylodynamic analysis of foot-and-mouth disease virus evolution in Mar Chiquita, Argentina.

Foot-and-mouth disease is a highly contagious disease affecting cloven-hoofed animals, resulting in considerable economic losses. Its causal agent is foot-and-mouth disease virus (FMDV), a picornavirus. Due to its error-prone replication and rapid evolution, the transmission and evolutionary dynamics of FMDV can be studied using genomic epidemiological approaches. To analyze FMDV evolution and identify possible transmission routes in an Argentinean region, field samples that tested positive for FMDV by PCR were obtained from 21 farms located in the Mar Chiquita district. Whole FMDV genome sequences were obtained by PCR amplification in seven fragments and sequencing using the Sanger technique. The genome sequences obtained from these samples were then analyzed using phylogenetic, phylogeographic, and evolutionary approaches. Three local transmission clusters were detected among the sampled viruses. The dataset was analyzed using Bayesian phylodynamic methods with appropriate coalescent and relaxed molecular clock models. The estimated mean viral evolutionary rate was 1.17 × 10- 2 substitutions/site/year. No significant differences in the rate of viral evolution were observed between farms with vaccinated animals and those with unvaccinated animals. The most recent common ancestor of the sampled sequences was dated to approximately one month before the first reported case in the outbreak. Virus transmission started in the south of the district and later dispersed to the west, and finally arrived in the east. Different transmission routes among the studied herds, such as non-replicating vectors and close contact contagion (i.e., aerosols), may be responsible for viral spread.

Congruity of genomic and epidemiological data in modelling of local cholera outbreaks.

Cholera continues to be a global health threat. Understanding how cholera spreads between locations is fundamental to the rational, evidence-based design of intervention and control efforts. Traditionally, cholera transmission models have used cholera case-count data. More recently, whole-genome sequence data have qualitatively described cholera transmission. Integrating these data streams may provide much more accurate models of cholera spread; however, no systematic analyses have been performed so far to compare traditional case-count models to the phylodynamic models from genomic data for cholera transmission. Here, we use high-fidelity case-count and whole-genome sequencing data from the 1991 to 1998 cholera epidemic in Argentina to directly compare the epidemiological model parameters estimated from these two data sources. We find that phylodynamic methods applied to cholera genomics data provide comparable estimates that are in line with established methods. Our methodology represents a critical step in building a framework for integrating case-count and genomic data sources for cholera epidemiology and other bacterial pathogens.

Reconstructing relative transmission rates inBayesian phylodynamics: Two-fold transmission advantage ofOmicron inBerlin, 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.

Population genomics of diarrheagenic Escherichia coli uncovers high connectivity between urban and rural communities in Ecuador.

Human movement may be an important driver of transmission dynamics for enteric pathogens but has largely been underappreciated except for international 'travelers' diarrhea or cholera. Phylodynamic methods, which combine genomic and epidemiological data, are used to examine rates and dynamics of disease matching underlying evolutionary history and biogeographic distributions, but these methods often are not applied to enteric bacterial pathogens. We used phylodynamics to explore the phylogeographic and evolutionary patterns of diarrheagenic E. coli in northern Ecuador to investigate the role of human travel in the geographic distribution of strains across the country. Using whole genome sequences of diarrheagenic E. coli isolates, we built a core genome phylogeny, reconstructed discrete ancestral states across urban and rural sites, and estimated migration rates between E. coli populations. We found minimal structuring based on site locations, urban vs. rural locality, pathotype, or clinical status. Ancestral states of phylogenomic nodes and tips were inferred to have 51% urban ancestry and 49% rural ancestry. Lack of structuring by location or pathotype E. coli isolates imply highly connected communities and extensive sharing of genomic characteristics across isolates. Using an approximate structured coalescent model, we estimated rates of migration among circulating isolates were 6.7 times larger for urban towards rural populations compared to rural towards urban populations. This suggests increased inferred migration rates of diarrheagenic E. coli from urban populations towards rural populations. Our results indicate that investments in water and sanitation prevention in urban areas could limit the spread of enteric bacterial pathogens among rural populations.

Lessons learned from community engagement regarding phylodynamic research with molecular HIV surveillance data.

The widespread implementation of molecular HIV surveillance (MHS) has resulted in an increased discussion about the ethical, human rights and public health implications of MHS. We narrate our process of pausing our research that uses data collected through MHS in response to these growing concerns and summarize the key lessons we learned through conversations with community members. The original study aimed to describe HIV transmission patterns by age and race/ethnicity among men who have sex with men in King County, Washington, by applying probabilistic phylodynamic modelling methods to HIV-1 pol gene sequences collected through MHS. In September 2020, we paused the publication of this research to conduct community engagement: we held two public-facing online presentations, met with a national community coalition that included representatives of networks of people living with HIV, and invited two members of this coalition to provide feedback on our manuscript. During each of these meetings, we shared a brief presentation of our methods and findings and explicitly solicited feedback on the perceived public health benefit and potential harm of our analyses and results. Some community concerns about MHS in public health practice also apply to research using MHS data, namely those related to informed consent, inference of transmission directionality and criminalization. Other critiques were specific to our research study and included feedback about the use of phylogenetic analyses to study assortativity by race/ethnicity and the importance of considering the broader context of stigma and structural racism. We ultimately decided the potential harms of publishing our study-perpetuating racialized stigma about men who have sex with men and eroding the trust between phylogenetics researchers and communities of people living with HIV-outweighed the potential benefits. HIV phylogenetics research using data collected through MHS data is a powerful scientific technology with the potential to benefit and harm communities of people living with HIV. Addressing criminalization and including people living with HIV in decision-making processes have the potential to meaningfully address community concerns and strengthen the ethical justification for using MHS data in both research and public health practice. We close with specific opportunities for action and advocacy by researchers.

Genetic determination of regional connectivity in modelling the spread of COVID-19 outbreak for more efficient mitigation strategies

For the COVID-19 pandemic, viral transmission has been documented in many historical and geographical contexts. Nevertheless, few studies have explicitly modeled the spatiotemporal flow based on genetic sequences, to develop mitigation strategies. Additionally, thousands of SARS-CoV-2 genomes have been sequenced with associated records, potentially providing a rich source for such spatiotemporal analysis, an unprecedented amount during a single outbreak. Here, in a case study of seven states, we model the first wave of the outbreak by determining regional connectivity from phylogenetic sequence information (i.e. “genetic connectivity”), in addition to traditional epidemiologic and demographic parameters. Our study shows nearly all of the initial outbreak can be traced to a few lineages, rather than disconnected outbreaks, indicative of a mostly continuous initial viral flow. While the geographic distance from hotspots is initially important in the modeling, genetic connectivity becomes increasingly significant later in the first wave. Moreover, our model predicts that isolated local strategies (e.g. relying on herd immunity) can negatively impact neighboring regions, suggesting more efficient mitigation is possible with unified, cross-border interventions. Finally, our results suggest that a few targeted interventions based on connectivity can have an effect similar to that of an overall lockdown. They also suggest that while successful lockdowns are very effective in mitigating an outbreak, less disciplined lockdowns quickly decrease in effectiveness. Our study provides a framework for combining phylodynamic and computational methods to identify targeted interventions.

Phylodynamics of SARS-CoV-2 in France, Europe, and the world in 2020

Although France was one of the most affected European countries by the COVID-19 pandemic in 2020, the dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) movement within France, but also involving France in Europe and in the world, remain only partially characterized in this timeframe. Here, we analyzed GISAID deposited sequences from January 1 to December 31, 2020 (n = 638,706 sequences at the time of writing). To tackle the challenging number of sequences without the bias of analyzing a single subsample of sequences, we produced 100 subsamples of sequences and related phylogenetic trees from the whole dataset for different geographic scales (worldwide, European countries, and French administrative regions) and time periods (from January 1 to July 25, 2020, and from July 26 to December 31, 2020). We applied a maximum likelihood discrete trait phylogeographic method to date exchange events (i.e., a transition from one location to another one), to estimate the geographic spread of SARS-CoV-2 transmissions and lineages into, from and within France, Europe, and the world. The results unraveled two different patterns of exchange events between the first and second half of 2020. Throughout the year, Europe was systematically associated with most of the intercontinental exchanges. SARS-CoV-2 was mainly introduced into France from North America and Europe (mostly by Italy, Spain, the United Kingdom, Belgium, and Germany) during the first European epidemic wave. During the second wave, exchange events were limited to neighboring countries without strong intercontinental movement, but Russia widely exported the virus into Europe during the summer of 2020. France mostly exported B.1 and B.1.160 lineages, respectively, during the first and second European epidemic waves. At the level of French administrative regions, the Paris area was the main exporter during the first wave. But, for the second epidemic wave, it equally contributed to virus spread with Lyon area, the second most populated urban area after Paris in France. The main circulating lineages were similarly distributed among the French regions. To conclude, by enabling the inclusion of tens of thousands of viral sequences, this original phylodynamic method enabled us to robustly describe SARS-CoV-2 geographic spread through France, Europe, and worldwide in 2020.

The Omicron Lineages BA.1 and BA.2 (Betacoronavirus SARS-CoV-2) Have Repeatedly Entered Brazil through a Single Dispersal Hub.

Brazil currently ranks second in absolute deaths by COVID-19, even though most of its population has completed the vaccination protocol. With the introduction of Omicron in late 2021, the number of COVID-19 cases soared once again in the country. We investigated in this work how lineages BA.1 and BA.2 entered and spread in the country by sequencing 2173 new SARS-CoV-2 genomes collected between October 2021 and April 2022 and analyzing them in addition to more than 18,000 publicly available sequences with phylodynamic methods. We registered that Omicron was present in Brazil as early as 16 November 2021 and by January 2022 was already more than 99% of samples. More importantly, we detected that Omicron has been mostly imported through the state of São Paulo, which in turn dispersed the lineages to other states and regions of Brazil. This knowledge can be used to implement more efficient non-pharmaceutical interventions against the introduction of new SARS-CoV variants focused on surveillance of airports and ground transportation.

Model misspecification misleads inference of the spatial dynamics of disease outbreaks

Epidemiology has been transformed by the advent of Bayesian phylodynamic models that allow researchers to infer the geographic history of pathogen dispersal over a set of discrete geographic areas [1, 2]. These models provide powerful tools for understanding the spatial dynamics of disease outbreaks, but contain many parameters that are inferred from minimal geographic information (i.e., the single area in which each pathogen was sampled). Consequently, inferences under these models are inherently sensitive to our prior assumptions about the model parameters. Here, we demonstrate that the default priors used in empirical phylodynamic studies make strong and biologically unrealistic assumptions about the underlying geographic process. We provide empirical evidence that these unrealistic priors strongly (and adversely) impact commonly reported aspects of epidemiological studies, including: 1) the relative rates of dispersal between areas; 2) the importance of dispersal routes for the spread of pathogens among areas; 3) the number of dispersal events between areas, and; 4) the ancestral area in which a given outbreak originated. We offer strategies to avoid these problems, and develop tools to help researchers specify more biologically reasonable prior models that will realize the full potential of phylodynamic methods to elucidate pathogen biology and, ultimately, inform surveillance and monitoring policies to mitigate the impacts of disease outbreaks.

State-dependent evolutionary models reveal modes of solid tumour growth

Spatial properties of tumour growth have profound implications for cancer progression, therapeutic resistance and metastasis. Yet, how spatial position governs tumour cell division remains difficult to evaluate in clinical tumours. Here, we demonstrate that faster division on the tumour periphery leaves characteristic genetic patterns, which become evident when a phylogenetic tree is reconstructed from spatially sampled cells. Namely, rapidly dividing peripheral lineages branch more extensively and acquire more mutations than slower-dividing centre lineages. We develop a Bayesian state-dependent evolutionary phylodynamic model (SDevo) that quantifies these patterns to infer the differential division rates between peripheral and central cells. We demonstrate that this approach accurately infers spatially varying birth rates of simulated tumours across a range of growth conditions and sampling strategies. We then show that SDevo outperforms state-of-the-art, non-cancer multi-state phylodynamic methods that ignore differential sequence evolution. Finally, we apply SDevo to single-time-point, multi-region sequencing data from clinical hepatocellular carcinomas and find evidence of a three- to six-times-higher division rate on the tumour edge. With the increasing availability of high-resolution, multi-region sequencing, we anticipate that SDevo will be useful in interrogating spatial growth restrictions and could be extended to model non-spatial factors that influence tumour progression.

Fitting stochastic epidemic models to gene genealogies using linear noise approximation.

Phylodynamics is a set of population genetics tools that aim at reconstructing demographic history of a population based on molecular sequences of individuals sampled from the population of interest. One important task in phylodynamics is to estimate changes in (effective) population size. When applied to infectious disease sequences such estimation of population size trajectories can provide information about changes in the number of infections. To model changes in the number of infected individuals, current phylodynamic methods use non-parametric approaches (e.g., Bayesian curve-fitting based on change-point models or Gaussian process priors), parametric approaches (e.g., based on differential equations), and stochastic modeling in conjunction with likelihood-free Bayesian methods. The first class of methods yields results that are hard to interpret epidemiologically. The second class of methods provides estimates of important epidemiological parameters, such as infection and removal/recovery rates, but ignores variation in the dynamics of infectious disease spread. The third class of methods is the most advantageous statistically, but relies on computationally intensive particle filtering techniques that limits its applications. We propose a Bayesian model that combines phylodynamic inference and stochastic epidemic models, and achieves computational tractability by using a linear noise approximation (LNA) - a technique that allows us to approximate probability densities of stochastic epidemic model trajectories. LNA opens the door for using modern Markov chain Monte Carlo tools to approximate the joint posterior distribution of the disease transmission parameters and of high dimensional vectors describing unobserved changes in the stochastic epidemic model compartment sizes (e.g., numbers of infectious and susceptible individuals). In a simulation study, we show that our method can successfully recover parameters of stochastic epidemic models. We apply our estimation technique to Ebola genealogies estimated using viral genetic data from the 2014 epidemic in Sierra Leone and Liberia.

Searching for Fast-Growing HIV Variants Using Scalable Phylodynamic Methods

Searching for Fast-Growing HIV Variants Using Scalable Phylodynamic Methods

Phylodynamic signatures in the emergence of community-associated MRSA

Community-associated, methicillin-resistant Staphylococcus aureus (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (Re) and sustained transmission (Re > 1) coincided with spread of progenitor methicillin-susceptible S. aureus (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated S. aureus lineages from America, Asia, Australasia, and Europe. Surges in Re were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers.

Phylodynamic analysis of SARS-CoV-2 spread in Rio de Janeiro, Brazil, highlights how metropolitan areas act as dispersal hubs for new variants.

During the first semester of 2021, all of Brazil has suffered an intense wave of COVID-19 associated with the Gamma variant. In July, the first cases of Delta variant were detected in the state of Rio de Janeiro. In this work, we have employed phylodynamic methods to analyse more than 1 600 genomic sequences of Delta variant collected until September in Rio de Janeiro to reconstruct how this variant has surpassed Gamma and dispersed throughout the state. After the introduction of Delta, it has initially spread mostly in the homonymous city of Rio de Janeiro, the most populous of the state. In a second stage, dispersal occurred to mid- and long-range cities, which acted as new close-range hubs for spread. We observed that the substitution of Gamma by Delta was possibly caused by its higher viral load, a proxy for transmissibility. This variant turnover prompted a new surge in cases, but with lower lethality than was observed during the peak caused by Gamma. We reason that high vaccination rates in the state of Rio de Janeiro were possibly what prevented a higher number of deaths.