Untangling introductions and persistence in COVID-19 resurgence in Europe.

After the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting in late summer 2020 that was deadlier and more difficult to contain1. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave2. Here we build a phylogeographical model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the resurgence of COVID-19 in Europe. We inform this model using genomic, mobility and epidemiological data from 10 European countries and estimate that in many countries more than half of the lineages circulating in late summer resulted from new introductions since 15 June 2020. The success in onward transmission of newly introduced lineages was negatively associated with the local incidence of COVID-19 during this period. The pervasive spread of variants in summer 2020 highlights the threat of viral dissemination when restrictions are lifted, and this needs to be carefully considered in strategies to control the current spread of variants that are more transmissible and/or evade immunity. Our findings indicate that more effective and coordinated measures are required to contain the spread through cross-border travel even as vaccination is reducing disease burden.

Retrospective analysis of SARS-CoV-2 omicron invasion over delta in French regions in 2021-22: a status-based multi-variant model.

BACKGROUND: SARS-CoV-2 is a rapidly spreading disease affecting human life and the economy on a global scale. The disease has caused so far more then 5.5 million deaths. The omicron outbreak that emerged in Botswana in the south of Africa spread around the globe at further increased rates, and caused unprecedented SARS-CoV-2 infection incidences in several countries. At the start of December 2021 the first omicron cases were reported in France. METHODS: In this paper we investigate the spreading potential of this novel variant relatively to the delta variant that was also in circulation in France at that time. Using a dynamic multi-variant model accounting for cross-immunity through a status-based approach, we analyze screening data reported by Santé Publique France over 13 metropolitan French regions between 1st of December 2021 and the 30th of January 2022. During the investigated period, the delta variant was replaced by omicron in all metropolitan regions in approximately three weeks. The analysis conducted retrospectively allows us to consider the whole replacement time window and compare regions with different times of omicron introduction and baseline levels of variants' transmission potential. As large uncertainties regarding cross-immunity among variants persist, uncertainty analyses were carried out to assess its impact on our estimations. RESULTS: Assuming that 80% of the population was immunized against delta, a cross delta/omicron cross-immunity of 25% and an omicron generation time of 3.5 days, the relative strength of omicron to delta, expressed as the ratio of their respective reproduction rates, [Formula: see text], was found to range between 1.51 and 1.86 across regions. Uncertainty analysis on epidemiological parameters led to [Formula: see text] ranging from 1.57 to 2.34 on average over the metropolitan French regions, weighted by population size. CONCLUSIONS: Upon introduction, omicron spread rapidly through the French territory and showed a high fitness relative to delta. We documented considerable geographical heterogeneities on the spreading dynamics. The historical reconstruction of variant emergence dynamics provide valuable ground knowledge to face future variant emergence events.

Preparedness and vulnerability of African countries against importations of COVID-19: a modelling study.

BACKGROUND: The novel coronavirus disease 2019 (COVID-19) epidemic has spread from China to 25 countries. Local cycles of transmission have already occurred in 12 countries after case importation. In Africa, Egypt has so far confirmed one case. The management and control of COVID-19 importations heavily rely on a country's health capacity. Here we evaluate the preparedness and vulnerability of African countries against their risk of importation of COVID-19. METHODS: We used data on the volume of air travel departing from airports in the infected provinces in China and directed to Africa to estimate the risk of importation per country. We determined the country's capacity to detect and respond to cases with two indicators: preparedness, using the WHO International Health Regulations Monitoring and Evaluation Framework; and vulnerability, using the Infectious Disease Vulnerability Index. Countries were clustered according to the Chinese regions contributing most to their risk. FINDINGS: Countries with the highest importation risk (ie, Egypt, Algeria, and South Africa) have moderate to high capacity to respond to outbreaks. Countries at moderate risk (ie, Nigeria, Ethiopia, Sudan, Angola, Tanzania, Ghana, and Kenya) have variable capacity and high vulnerability. We identified three clusters of countries that share the same exposure to the risk originating from the provinces of Guangdong, Fujian, and the city of Beijing, respectively. INTERPRETATION: Many countries in Africa are stepping up their preparedness to detect and cope with COVID-19 importations. Resources, intensified surveillance, and capacity building should be urgently prioritised in countries with moderate risk that might be ill-prepared to detect imported cases and to limit onward transmission. FUNDING: EU Framework Programme for Research and Innovation Horizon 2020, Agence Nationale de la Recherche.

Tracing and analysis of 288 early SARS-CoV-2 infections outside China: A modeling study.

BACKGROUND: In the early months of 2020, a novel coronavirus disease (COVID-19) spread rapidly from China across multiple countries worldwide. As of March 17, 2020, COVID-19 was officially declared a pandemic by the World Health Organization. We collected data on COVID-19 cases outside China during the early phase of the pandemic and used them to predict trends in importations and quantify the proportion of undetected imported cases. METHODS AND FINDINGS: Two hundred and eighty-eight cases have been confirmed out of China from January 3 to February 13, 2020. We collected and synthesized all available information on these cases from official sources and media. We analyzed importations that were successfully isolated and those leading to onward transmission. We modeled their number over time, in relation to the origin of travel (Hubei province, other Chinese provinces, other countries) and interventions. We characterized the importation timeline to assess the rapidity of isolation and epidemiologically linked clusters to estimate the rate of detection. We found a rapid exponential growth of importations from Hubei, corresponding to a doubling time of 2.8 days, combined with a slower growth from the other areas. We predicted a rebound of importations from South East Asia in the successive weeks. Time from travel to detection has considerably decreased since first importation, from 14.5 ± 5.5 days on January 5, 2020, to 6 ± 3.5 days on February 1, 2020. However, we estimated 36% of detection of imported cases. This study is restricted to the early phase of the pandemic, when China was the only large epicenter and foreign countries had not discovered extensive local transmission yet. Missing information in case history was accounted for through modeling and imputation. CONCLUSIONS: Our findings indicate that travel bans and containment strategies adopted in China were effective in reducing the exportation growth rate. However, the risk of importation was estimated to increase again from other sources in South East Asia. Surveillance and management of traveling cases represented a priority in the early phase of the epidemic. With the majority of imported cases going undetected (6 out of 10), countries experienced several undetected clusters of chains of local transmissions, fueling silent epidemics in the community. These findings become again critical to prevent second waves, now that countries have reduced their epidemic activity and progressively phase out lockdown.

Host contact dynamics shapes richness and dominance of pathogen strains.

The interaction among multiple microbial strains affects the spread of infectious diseases and the efficacy of interventions. Genomic tools have made it increasingly easy to observe pathogenic strains diversity, but the best interpretation of such diversity has remained difficult because of relationships with host and environmental factors. Here, we focus on host-to-host contact behavior and study how it changes populations of pathogens in a minimal model of multi-strain interaction. We simulated a population of identical strains competing by mutual exclusion and spreading on a dynamical network of hosts according to a stochastic susceptible-infectious-susceptible model. We computed ecological indicators of diversity and dominance in strain populations for a collection of networks illustrating various properties found in real-world examples. Heterogeneities in the number of contacts among hosts were found to reduce diversity and increase dominance by making the repartition of strains among infected hosts more uneven, while strong community structure among hosts increased strain diversity. We found that the introduction of strains associated with hosts entering and leaving the system led to the highest pathogenic richness at intermediate turnover levels. These results were finally illustrated using the spread of Staphylococcus aureus in a long-term health-care facility where close proximity interactions and strain carriage were collected simultaneously. We found that network structural and temporal properties could account for a large part of the variability observed in strain diversity. These results show how stochasticity and network structure affect the population ecology of pathogens and warn against interpreting observations as unambiguous evidence of epidemiological differences between strains.

Novel coronavirus (2019-nCoV) early-stage importation risk to Europe, January 2020.

As at 27 January 2020, 42 novel coronavirus (2019-nCoV) cases were confirmed outside China. We estimate the risk of case importation to Europe from affected areas in China via air travel. We consider travel restrictions in place, three reported cases in France, one in Germany. Estimated risk in Europe remains high. The United Kingdom, Germany and France are at highest risk. Importation from Beijing and Shanghai would lead to higher and widespread risk for Europe.

Excess cases of influenza-like illnesses synchronous with coronavirus disease (COVID-19) epidemic, France, March 2020.

Several French regions where coronavirus disease (COVID-19) has been reported currently show a renewed increase in ILI cases in the general practice-based Sentinelles network. We computed the number of excess cases by region from 24 February to 8 March 2020 and found a correlation with the number of reported COVID-19 cases so far. The data suggest larger circulation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the French population than apparent from confirmed cases.

Epidemic Threshold in Continuous-Time Evolving Networks.

Current understanding of the critical outbreak condition on temporal networks relies on approximations (time scale separation, discretization) that may bias the results. We propose a theoretical framework to compute the epidemic threshold in continuous time through the infection propagator approach. We introduce the weak commutation condition allowing the interpretation of annealed networks, activity-driven networks, and time scale separation into one formalism. Our work provides a coherent connection between discrete and continuous time representations applicable to realistic scenarios.

The impact of regular school closure on seasonal influenza epidemics: a data-driven spatial transmission model for Belgium.

BACKGROUND: School closure is often considered as an option to mitigate influenza epidemics because of its potential to reduce transmission in children and then in the community. The policy is still however highly debated because of controversial evidence. Moreover, the specific mechanisms leading to mitigation are not clearly identified. METHODS: We introduced a stochastic spatial age-specific metapopulation model to assess the role of holiday-associated behavioral changes and how they affect seasonal influenza dynamics. The model is applied to Belgium, parameterized with country-specific data on social mixing and travel, and calibrated to the 2008/2009 influenza season. It includes behavioral changes occurring during weekend vs. weekday, and holiday vs. school-term. Several experimental scenarios are explored to identify the relevant social and behavioral mechanisms. RESULTS: Stochastic numerical simulations show that holidays considerably delay the peak of the season and mitigate its impact. Changes in mixing patterns are responsible for the observed effects, whereas changes in travel behavior do not alter the epidemic. Weekends are important in slowing down the season by periodically dampening transmission. Christmas holidays have the largest impact on the epidemic, however later school breaks may help in reducing the epidemic size, stressing the importance of considering the full calendar. An extension of the Christmas holiday of 1 week may further mitigate the epidemic. CONCLUSION: Changes in the way individuals establish contacts during holidays are the key ingredient explaining the mitigating effect of regular school closure. Our findings highlight the need to quantify these changes in different demographic and epidemic contexts in order to provide accurate and reliable evaluations of closure effectiveness. They also suggest strategic policies in the distribution of holiday periods to minimize the epidemic impact.

The potential value of crowdsourced surveillance systems in supplementing sentinel influenza networks: the case of France.

IntroductionParticipatory surveillance systems provide rich crowdsourced data, profiling individuals and their health status at a given time. We explored the usefulness of data from GrippeNet.fr, a participatory surveillance system, to estimate influenza-related illness incidence in France. Methods: GrippeNet.fr is an online cohort since 2012 averaging ca. 5,000 weekly participants reporting signs/symptoms suggestive of influenza. GrippeNet.fr has flexible criteria to define influenza-related illness. Different case definitions based on reported signs/symptoms and inclusions of criteria accounting for individuals' reporting and participation were used to produce influenza-related illness incidence estimates, which were compared to those from sentinel networks. We focused on the 2012/13 and 2013/14 seasons when two sentinel networks, monitoring influenza-like-illness (ILI) and acute respiratory infections (ARI) existed in France. Results: GrippeNet.fr incidence estimates agreed well with official temporal trends, with a higher accuracy for ARI than ILI. The influenza epidemic peak was often anticipated by one week, despite irregular participation of individuals. The European Centre for Disease Prevention and Control ILI definition, commonly used by participatory surveillance in Europe, performed better in tracking ARI than ILI when applied to GrippeNet.fr data. Conclusion: Evaluation of the epidemic intensity from crowdsourced data requires epidemic and intensity threshold estimations from several consecutive seasons. The study provides a standardised analytical framework for crowdsourced surveillance showing high sensitivity in detecting influenza-related changes in the population. It contributes to improve the comparability of epidemics across seasons and with sentinel systems. In France, GrippeNet.fr may supplement the ILI sentinel network after ARI surveillance discontinuation in 2014.

Predicting epidemic risk from past temporal contact data.

Understanding how epidemics spread in a system is a crucial step to prevent and control outbreaks, with broad implications on the system's functioning, health, and associated costs. This can be achieved by identifying the elements at higher risk of infection and implementing targeted surveillance and control measures. One important ingredient to consider is the pattern of disease-transmission contacts among the elements, however lack of data or delays in providing updated records may hinder its use, especially for time-varying patterns. Here we explore to what extent it is possible to use past temporal data of a system's pattern of contacts to predict the risk of infection of its elements during an emerging outbreak, in absence of updated data. We focus on two real-world temporal systems; a livestock displacements trade network among animal holdings, and a network of sexual encounters in high-end prostitution. We define the node's loyalty as a local measure of its tendency to maintain contacts with the same elements over time, and uncover important non-trivial correlations with the node's epidemic risk. We show that a risk assessment analysis incorporating this knowledge and based on past structural and temporal pattern properties provides accurate predictions for both systems. Its generalizability is tested by introducing a theoretical model for generating synthetic temporal networks. High accuracy of our predictions is recovered across different settings, while the amount of possible predictions is system-specific. The proposed method can provide crucial information for the setup of targeted intervention strategies.

Risk of MERS importation and onward transmission: a systematic review and analysis of cases reported to WHO.

BACKGROUND: The continuing circulation of MERS in the Middle East makes the international dissemination of the disease a permanent threat. To inform risk assessment, we investigated the spatiotemporal pattern of MERS global dissemination and looked for factors explaining the heterogeneity observed in transmission events following importation. METHODS: We reviewed imported MERS cases worldwide up to July 2015. We modelled importations in time based on air travel combined with incidence in Middle East. We used the detailed history of MERS case management after importation (time to hospitalization and isolation, number of hospitals visited,…) in logistic regression to identify risk factors for secondary transmission. We assessed changes in time to hospitalization and isolation in relation to collective and public health attention to the epidemic, measured by three indicators (Google Trends, ProMED-mail, Disease Outbreak News). RESULTS: Modelled importation events were found to reproduce both the temporal and geographical structure of those observed - the Pearson correlation coefficient between predicted and observed monthly time series was large (r = 0.78, p < 10(-4)). The risk of secondary transmission following importation increased with the time to case isolation or death (OR = 1.7 p = 0.04) and more precisely with the duration of hospitalization (OR = 1.7, p = 0.02). The average daily number of secondary cases was 0.02 [0.0,0.12] in the community and 0.20 [0.03,9.0] in the hospital. Time from hospitalisation to isolation decreased in periods of high public health attention (2.33 ± 0.34 vs. 6.44 ± 0.97 days during baseline attention). CONCLUSIONS: Countries at risk of importation should focus their resources on strict infection control measures for the management of potential cases in healthcare settings and on prompt MERS cases identification. Individual and collective awareness are key to substantially improve such preparedness.

Recalibrating disease parameters for increasing realism in modeling epidemics in closed settings.

BACKGROUND: The homogeneous mixing assumption is widely adopted in epidemic modelling for its parsimony and represents the building block of more complex approaches, including very detailed agent-based models. The latter assume homogeneous mixing within schools, workplaces and households, mostly for the lack of detailed information on human contact behaviour within these settings. The recent data availability on high-resolution face-to-face interactions makes it now possible to assess the goodness of this simplified scheme in reproducing relevant aspects of the infection dynamics. METHODS: We consider empirical contact networks gathered in different contexts, as well as synthetic data obtained through realistic models of contacts in structured populations. We perform stochastic spreading simulations on these contact networks and in populations of the same size under a homogeneous mixing hypothesis. We adjust the epidemiological parameters of the latter in order to fit the prevalence curve of the contact epidemic model. We quantify the agreement by comparing epidemic peak times, peak values, and epidemic sizes. RESULTS: Good approximations of the peak times and peak values are obtained with the homogeneous mixing approach, with a median relative difference smaller than 20 % in all cases investigated. Accuracy in reproducing the peak time depends on the setting under study, while for the peak value it is independent of the setting. Recalibration is found to be linear in the epidemic parameters used in the contact data simulations, showing changes across empirical settings but robustness across groups and population sizes. CONCLUSIONS: An adequate rescaling of the epidemiological parameters can yield a good agreement between the epidemic curves obtained with a real contact network and a homogeneous mixing approach in a population of the same size. The use of such recalibrated homogeneous mixing approximations would enhance the accuracy and realism of agent-based simulations and limit the intrinsic biases of the homogeneous mixing.

Host mobility drives pathogen competition in spatially structured populations.

Interactions among multiple infectious agents are increasingly recognized as a fundamental issue in the understanding of key questions in public health regarding pathogen emergence, maintenance, and evolution. The full description of host-multipathogen systems is, however, challenged by the multiplicity of factors affecting the interaction dynamics and the resulting competition that may occur at different scales, from the within-host scale to the spatial structure and mobility of the host population. Here we study the dynamics of two competing pathogens in a structured host population and assess the impact of the mobility pattern of hosts on the pathogen competition. We model the spatial structure of the host population in terms of a metapopulation network and focus on two strains imported locally in the system and having the same transmission potential but different infectious periods. We find different scenarios leading to competitive success of either one of the strain or to the codominance of both strains in the system. The dominance of the strain characterized by the shorter or longer infectious period depends exclusively on the structure of the population and on the the mobility of hosts across patches. The proposed modeling framework allows the integration of other relevant epidemiological, environmental and demographic factors, opening the path to further mathematical and computational studies of the dynamics of multipathogen systems.

Metapopulation epidemic models with heterogeneous mixing and travel behaviour.

BACKGROUND: Determining the pandemic potential of an emerging infectious disease and how it depends on the various epidemic and population aspects is critical for the preparation of an adequate response aimed at its control. The complex interplay between population movements in space and non-homogeneous mixing patterns have so far hindered the fundamental understanding of the conditions for spatial invasion through a general theoretical framework. To address this issue, we present an analytical modelling approach taking into account such interplay under general conditions of mobility and interactions, in the simplifying assumption of two population classes. METHODS: We describe a spatially structured population with non-homogeneous mixing and travel behaviour through a multi-host stochastic epidemic metapopulation model. Different population partitions, mixing patterns and mobility structures are considered, along with a specific application for the study of the role of age partition in the early spread of the 2009 H1N1 pandemic influenza. RESULTS: We provide a complete mathematical formulation of the model and derive a semi-analytical expression of the threshold condition for global invasion of an emerging infectious disease in the metapopulation system. A rich solution space is found that depends on the social partition of the population, the pattern of contacts across groups and their relative social activity, the travel attitude of each class, and the topological and traffic features of the mobility network. Reducing the activity of the less social group and reducing the cross-group mixing are predicted to be the most efficient strategies for controlling the pandemic potential in the case the less active group constitutes the majority of travellers. If instead traveling is dominated by the more social class, our model predicts the existence of an optimal across-groups mixing that maximises the pandemic potential of the disease, whereas the impact of variations in the activity of each group is less important. CONCLUSIONS: The proposed modelling approach introduces a theoretical framework for the study of infectious diseases spread in a population with two layers of heterogeneity relevant for the local transmission and the spatial propagation of the disease. It can be used for pandemic preparedness studies to identify adequate interventions and quantitatively estimate the corresponding required effort, as well as in an emerging epidemic situation to assess the pandemic potential of the pathogen from population and early outbreak data.

The representativeness of a European multi-center network for influenza-like-illness participatory surveillance.

BACKGROUND: The Internet is becoming more commonly used as a tool for disease surveillance. Similarly to other surveillance systems and to studies using online data collection, Internet-based surveillance will have biases in participation, affecting the generalizability of the results. Here we quantify the participation biases of Influenzanet, an ongoing European-wide network of Internet-based participatory surveillance systems for influenza-like-illness. METHODS: In 2011/2012 Influenzanet launched a standardized common framework for data collection applied to seven European countries. Influenzanet participants were compared to the general population of the participating countries to assess the representativeness of the sample in terms of a set of demographic, geographic, socio-economic and health indicators. RESULTS: More than 30,000 European residents registered to the system in the 2011/2012 season, and a subset of 25,481 participants were selected for this study. All age classes (10 years brackets) were represented in the cohort, including under 10 and over 70 years old. The Influenzanet population was not representative of the general population in terms of age distribution, underrepresenting the youngest and oldest age classes. The gender imbalance differed between countries. A counterbalance between gender-specific information-seeking behavior (more prominent in women) and Internet usage (with higher rates in male populations) may be at the origin of this difference. Once adjusted by demographic indicators, a similar propensity to commute was observed for each country, and the same top three transportation modes were used for six countries out of seven. Smokers were underrepresented in the majority of countries, as were individuals with diabetes; the representativeness of asthma prevalence and vaccination coverage for 65+ individuals in two successive seasons (2010/2011 and 2011/2012) varied between countries. CONCLUSIONS: Existing demographic and national datasets allowed the quantification of the participation biases of a large cohort for influenza-like-illness surveillance in the general population. Significant differences were found between Influenzanet participants and the general population. The quantified biases need to be taken into account in the analysis of Influenzanet epidemiological studies and provide indications on populations groups that should be targeted in recruitment efforts.

Integrating dynamical modeling and phylogeographic inference to characterize global influenza circulation.

Global seasonal influenza circulation involves a complex interplay between local (seasonality, demography, host immunity) and global factors (international mobility) shaping recurrent epidemic patterns. No studies so far have reconciled the two spatial levels, evaluating the coupling between national epidemics, considering heterogeneous coverage of epidemiological and virological data, integrating different data sources. We propose a novel combined approach based on a dynamical model of global influenza spread (GLEAM), integrating high-resolution demographic and mobility data, and a generalized linear model of phylogeographic diffusion that accounts for time-varying migration rates. Seasonal migration fluxes across global macro-regions simulated with GLEAM are tested as phylogeographic predictors to provide model validation and calibration based on genetic data. Seasonal fluxes obtained with a specific transmissibility peak time and recurrent travel outperformed the raw air-transportation predictor, previously considered as optimal indicator of global influenza migration. Influenza A subtypes supported autumn-winter reproductive number as high as 2.25 and an average immunity duration of 2 years. Similar dynamics were preferred by influenza B lineages, with a lower autumn-winter reproductive number. Comparing simulated epidemic profiles against FluNet data offered comparatively limited resolution power. The multiscale approach enables model selection yielding a novel computational framework for describing global influenza dynamics at different scales - local transmission and national epidemics vs. international coupling through mobility and imported cases. Our findings have important implications to improve preparedness against seasonal influenza epidemics. The approach can be generalized to other epidemic contexts, such as emerging disease outbreaks to improve the flexibility and predictive power of modeling.

Drivers and impact of the early silent invasion of SARS-CoV-2 Alpha.

SARS-CoV-2 variants of concern (VOCs) circulated cryptically before being identified as a threat, delaying interventions. Here we studied the drivers of such silent spread and its epidemic impact to inform future response planning. We focused on Alpha spread out of the UK. We integrated spatio-temporal records of international mobility, local epidemic growth and genomic surveillance into a Bayesian framework to reconstruct the first three months after Alpha emergence. We found that silent circulation lasted from days to months and decreased with the logarithm of sequencing coverage. Social restrictions in some countries likely delayed the establishment of local transmission, mitigating the negative consequences of late detection. Revisiting the initial spread of Alpha supports local mitigation at the destination in case of emerging events.

Dispersion patterns of SARS-CoV-2 variants Gamma, Lambda and Mu in Latin America and the Caribbean.

Latin America and Caribbean (LAC) regions were an important epicenter of the COVID-19 pandemic and SARS-CoV-2 evolution. Through the COVID-19 Genomic Surveillance Regional Network (COVIGEN), LAC countries produced an important number of genomic sequencing data that made possible an enhanced SARS-CoV-2 genomic surveillance capacity in the Americas, paving the way for characterization of emerging variants and helping to guide the public health response. In this study we analyzed approximately 300,000 SARS-CoV-2 sequences generated between February 2020 and March 2022 by multiple genomic surveillance efforts in LAC and reconstructed the diffusion patterns of the main variants of concern (VOCs) and of interest (VOIs) possibly originated in the Region. Our phylogenetic analysis revealed that the spread of variants Gamma, Lambda and Mu reflects human mobility patterns due to variations of international air passenger transportation and gradual lifting of social distance measures previously implemented in countries. Our results highlight the potential of genetic data to reconstruct viral spread and unveil preferential routes of viral migrations that are shaped by human mobility patterns.

The genomic evolutionary dynamics and global circulation patterns of respiratory syncytial virus.

Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infection in young children and the second leading cause of infant death worldwide. While global circulation has been extensively studied for respiratory viruses such as seasonal influenza, and more recently also in great detail for SARS-CoV-2, a lack of global multi-annual sampling of complete RSV genomes limits our understanding of RSV molecular epidemiology. Here, we capitalise on the genomic surveillance by the INFORM-RSV study and apply phylodynamic approaches to uncover how selection and neutral epidemiological processes shape RSV diversity. Using complete viral genome sequences, we show similar patterns of site-specific diversifying selection among RSVA and RSVB and recover the imprint of non-neutral epidemic processes on their genealogies. Using a phylogeographic approach, we provide evidence for air travel governing the global patterns of RSVA and RSVB spread, which results in a considerable degree of phylogenetic mixing across countries. Our findings highlight the potential of systematic global RSV genomic surveillance for transforming our understanding of global RSV spread.