RESUMO
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.
Assuntos
COVID-19/transmissão , COVID-19/virologia , SARS-CoV-2/isolamento & purificação , COVID-19/epidemiologia , COVID-19/prevenção & controle , Europa (Continente)/epidemiologia , Genoma Viral/genética , Humanos , Incidência , Locomoção , Filogenia , Filogeografia , SARS-CoV-2/classificação , SARS-CoV-2/genética , SARS-CoV-2/patogenicidade , Fatores de Tempo , Viagem/estatística & dados numéricosRESUMO
Since the latter part of 2020, SARS-CoV-2 evolution has been characterised by the emergence of viral variants associated with distinct biological characteristics. While the main research focus has centred on the ability of new variants to increase in frequency and impact the effective reproductive number of the virus, less attention has been placed on their relative ability to establish transmission chains and to spread through a geographic area. Here, we describe a phylogeographic approach to estimate and compare the introduction and dispersal dynamics of the main SARS-CoV-2 variants - Alpha, Iota, Delta, and Omicron - that circulated in the New York City area between 2020 and 2022. Notably, our results indicate that Delta had a lower ability to establish sustained transmission chains in the NYC area and that Omicron (BA.1) was the variant fastest to disseminate across the study area. The analytical approach presented here complements non-spatially-explicit analytical approaches that seek a better understanding of the epidemiological differences that exist among successive SARS-CoV-2 variants of concern.
Assuntos
COVID-19 , Humanos , COVID-19/epidemiologia , COVID-19/virologia , Cidade de Nova Iorque/epidemiologia , SARS-CoV-2/genéticaRESUMO
The dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron-BA.1 variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the United States became increasingly significant. The number of detected introductions varied from 96 and 101 for Alpha and Delta to 39 for Omicron-BA.1. Most of these introductions left a low number of descendants (<10), suggesting a limited impact on the evolution of the pandemic in Galicia. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.
Assuntos
COVID-19 , SARS-CoV-2 , Espanha/epidemiologia , COVID-19/epidemiologia , COVID-19/transmissão , COVID-19/virologia , Humanos , SARS-CoV-2/genética , Genoma Viral , Filogenia , PandemiasRESUMO
During the first phase of the COVID-19 epidemic, New York City rapidly became the epicenter of the pandemic in the United States. While molecular phylogenetic analyses have previously highlighted multiple introductions and a period of cryptic community transmission within New York City, little is known about the circulation of SARS-CoV-2 within and among its boroughs. We here perform phylogeographic investigations to gain insights into the circulation of viral lineages during the first months of the New York City outbreak. Our analyses describe the dispersal dynamics of viral lineages at the state and city levels, illustrating that peripheral samples likely correspond to distinct dispersal events originating from the main metropolitan city areas. In line with the high prevalence recorded in this area, our results highlight the relatively important role of the borough of Queens as a transmission hub associated with higher local circulation and dispersal of viral lineages toward the surrounding boroughs.
Assuntos
COVID-19/epidemiologia , COVID-19/transmissão , SARS-CoV-2/genética , Genoma Viral/genética , Humanos , Cidade de Nova Iorque/epidemiologia , Filogenia , Filogeografia , Prevalência , SARS-CoV-2/classificação , SARS-CoV-2/isolamento & purificaçãoRESUMO
BackgroundThe earliest recognised infections by the SARS-CoV-2 Omicron variant (Pango lineage B.1.1.529) in Belgium and Switzerland suggested a connection to an international water polo tournament, held 12-14 November 2021 in Brno, Czechia.AimTo study the arrival and subsequent spread of the Omicron variant in Belgium and Switzerland, and understand the overall importance of this international sporting event on the number of infections in the two countries.MethodsWe performed intensive forward and backward contact tracing in both countries, supplemented by phylogenetic investigations using virus sequences of the suspected infection chain archived in public databases.ResultsThrough contact tracing, we identified two and one infected athletes of the Belgian and Swiss water polo teams, respectively, and subsequently also three athletes from Germany. In Belgium and Switzerland, four and three secondary infections, and three and one confirmed tertiary infections were identified. Phylogenetic investigation demonstrated that this sporting event played a role as the source of infection, but without a direct link with infections from South Africa and not as a superspreading event; the virus was found to already be circulating at that time in the countries involved.ConclusionThe SARS-CoV-2 Omicron variant started to circulate in Europe several weeks before its identification in South Africa on 24 November 2021. Accordingly, it can be assumed that travel restrictions are usually implemented too late to prevent the spread of newly detected SARS-CoV-2 variants to other regions. Phylogenetic analysis may modify the perception of an apparently clear result of intensive contact tracing.
Assuntos
COVID-19 , Esportes Aquáticos , Humanos , SARS-CoV-2/genética , Bélgica/epidemiologia , Suíça/epidemiologia , República Tcheca , Filogenia , COVID-19/epidemiologia , AlemanhaRESUMO
Since the start of the COVID-19 pandemic, an unprecedented number of genomic sequences of SARS-CoV-2 have been generated and shared with the scientific community. The unparalleled volume of available genetic data presents a unique opportunity to gain real-time insights into the virus transmission during the pandemic, but also a daunting computational hurdle if analyzed with gold-standard phylogeographic approaches. To tackle this practical limitation, we here describe and apply a rapid analytical pipeline to analyze the spatiotemporal dispersal history and dynamics of SARS-CoV-2 lineages. As a proof of concept, we focus on the Belgian epidemic, which has had one of the highest spatial densities of available SARS-CoV-2 genomes. Our pipeline has the potential to be quickly applied to other countries or regions, with key benefits in complementing epidemiological analyses in assessing the impact of intervention measures or their progressive easement.
Assuntos
COVID-19/transmissão , COVID-19/virologia , Genoma Viral , Filogeografia , SARS-CoV-2/genética , Bélgica , COVID-19/epidemiologia , Evolução Molecular , Genômica , Humanos , Funções Verossimilhança , Mutação , Isolamento de Pacientes , Filogenia , Distanciamento Físico , Análise Espaço-Temporal , Fluxo de TrabalhoRESUMO
Multidrug-resistant Salmonella enterica subspecies enterica 4,[5],12:i:- sequence type 34 represents a worldwide public health risk. To determine its origin in the United States, we reconstructed a time-scaled phylogeny with a discrete trait geospatial model. The clone in the United States was introduced from Europe on multiple occasions in the early 2000s.
Assuntos
Salmonelose Animal , Salmonella enterica , Animais , Antibacterianos/farmacologia , Farmacorresistência Bacteriana Múltipla/genética , Europa (Continente)/epidemiologia , Testes de Sensibilidade Microbiana , Salmonella enterica/genética , Estados Unidos/epidemiologiaRESUMO
Nontyphoidal Salmonella enterica (NTS) poses a major public health risk worldwide that is amplified by the existence of antimicrobial-resistant strains, especially those resistant to quinolones and extended-spectrum cephalosporins (ESC). Little is known on the dissemination of plasmids harboring the acquired genetic determinants that confer resistance to these antimicrobials across NTS serotypes from livestock in the United States. NTS isolates (n = 183) from U.S. swine clinical cases retrieved during 2014 to 2016 were selected for sequencing based on their phenotypic resistance to enrofloxacin (quinolone) or ceftiofur (3rd-generation cephalosporin). De novo assemblies were used to identify chromosomal mutations and acquired antimicrobial resistance genes (AARGs). In addition, plasmids harboring AARGs were identified using short-read assemblies and characterized using a multistep approach that was validated by long-read sequencing. AARGs to quinolones [qnrB15, qnrB19, qnrB2, qnrD, qnrS1, qnrS2, and aac(6')Ib-cr] and ESC (blaCMY-2, blaCTX-M-1, blaCTX-M-27, and blaSHV-12) were distributed across serotypes and were harbored by several plasmids. In addition, chromosomal mutations associated with resistance to quinolones were identified in the target enzyme and efflux pump regulation genes. The predominant plasmid harboring the prevalent qnrB19 gene was distributed across serotypes. It was identical to a plasmid previously reported in S. enterica serovar Anatum from swine in the United States (GenBank accession number KY991369.1) and similar to Escherichia coli plasmids from humans in South America (GenBank accession numbers GQ374157.1 and JN979787.1). Our findings suggest that plasmids harboring AARGs encoding mechanisms of resistance to critically important antimicrobials are present in multiple NTS serotypes circulating in swine in the United States and can contribute to resistance expansion through horizontal transmission.
Assuntos
Resistência às Cefalosporinas/genética , Cefalosporinas/farmacologia , Plasmídeos/genética , Quinolonas/farmacologia , Salmonella enterica/genética , Animais , Antibacterianos/farmacologia , Proteínas de Bactérias/genética , Farmacorresistência Bacteriana Múltipla/genética , Enrofloxacina/farmacologia , Escherichia coli/efeitos dos fármacos , Escherichia coli/genética , Testes de Sensibilidade Microbiana/métodos , Salmonella enterica/efeitos dos fármacos , Sorogrupo , América do Sul , Suínos , Estados UnidosRESUMO
Phylogeographic analyses are able to exploit the location data associated with sampled molecular sequences to reconstruct the spatio-temporal dispersal history of a pathogen. Visualisation software is commonly used to facilitate the interpretation of the accompanying estimation results, as these are not always easily interpretable. spread.gl is a powerful, open-source and feature-rich browser application that enables smooth, intuitive and user-friendly visualisation of both discrete and continuous phylogeographic inference results, enabling the animation of pathogen geographic dispersal through time. spread.gl can render and combine the visualisation of several data layers, including a geographic layer (e.g., a world map), multiple layers that contain information extracted from the input phylogeny, and different types of layers that represent environmental data. As such, users can explore which environmental data may have shaped pathogen dispersal patterns, that can subsequently be formally tested through more principled statistical analyses. We showcase the visualisation features of spread.gl on several representative pathogen dispersal examples, including the smooth animation of a phylogeny encompassing over 17,000 genomic sequences resulting from a large-scale SARS-CoV-2 analysis.
RESUMO
Different factors influence the spread of SARS-CoV-2, from the inherent transmission capabilities of the different variants to the control measurements put in place. Here we studied the introduction of the Alpha, Delta, and Omicron-BA.1 variants of concern (VOCs) into Spain. For this, we collected genomic data from the GISAID database and combined it with connectivity data from different countries with Spain to perform a phylodynamic Bayesian analysis of the introductions. Our findings reveal that the introductions of these VOCs predominantly originated from France, especially in the case of Alpha. As travel restrictions were eased during the Delta and Omicron-BA.1 waves, the number of introductions from distinct countries increased, with the United Kingdom and Germany becoming significant sources of the virus. The largest number of introductions detected corresponded to the Delta wave, which was associated with fewer restrictions and the summer period, when Spain receives a considerable number of tourists. This research underscores the importance of monitoring international travel patterns and implementing targeted public health measures to manage the spread of SARS-CoV-2.
RESUMO
The dynamics of SARS-CoV-2 transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the USA became increasingly significant. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.
RESUMO
The emergence of recombinant viruses is a threat to public health, as recombination may integrate variant-specific features that together result in escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown. We identified a Delta-Omicron (AY.45-BA.1) recombinant in an immunosuppressed transplant recipient treated with monoclonal antibody Sotrovimab. The single recombination breakpoint is located in the spike N-terminal domain adjacent to the Sotrovimab binding site. While Delta and BA.1 are sensitive to Sotrovimab neutralization, the Delta-Omicron recombinant is highly resistant. To our knowledge, this is the first described instance of recombination between circulating SARS-CoV-2 variants as a functional mechanism of resistance to treatment and immune escape.
RESUMO
Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused millions of infections and fatalities globally since its emergence in late 2019. The virus was first detected in Finland in January 2020, after which it rapidly spread among the populace in spring. However, compared to other European nations, Finland has had a low incidence of SARS-CoV-2. To gain insight into the origins and turnover of SARS-CoV-2 lineages circulating in Finland in 2020, we investigated the phylogeographic and -dynamic history of the virus. Methods: The origins of SARS-CoV-2 introductions were inferred via Travel-aware Bayesian time-measured phylogeographic analyses. Sequences for the analyses included virus genomes belonging to the B.1 lineage and with the D614G mutation from countries of likely origin, which were determined utilizing Google mobility data. We collected all available sequences from spring and fall peaks to study lineage dynamics. Results: We observed rapid turnover among Finnish lineages during this period. Clade 20C became the most prevalent among sequenced cases and was replaced by other strains in fall 2020. Bayesian phylogeographic reconstructions suggested 42 independent introductions into Finland during spring 2020, mainly from Italy, Austria, and Spain. Conclusions: A single introduction from Spain might have seeded one-third of cases in Finland during spring in 2020. The investigations of the original introductions of SARS-CoV-2 to Finland during the early stages of the pandemic and of the subsequent lineage dynamics could be utilized to assess the role of transboundary movements and the effects of early intervention and public health measures.
RESUMO
Background: The emergence of recombinant viruses is a threat to public health. Recombination of viral variants may combine variant-specific features that together catalyze viral escape from treatment or immunity. The selective advantages of recombinant SARS-CoV-2 isolates over their parental lineages remain unknown. Methods: Multi-method amplicon and metagenomic sequencing of a clinical swab and the in vitro grown virus allowed for high-confidence detection of a novel recombinant variant. Mutational, phylogeographic, and structural analyses determined features of the recombinant genome and spike protein. Neutralization assays using infectious as well as pseudotyped viruses and point mutants thereof defined the recombinant's sensitivity to a panel of monoclonal antibodies and sera from vaccinated and/or convalescent individuals. Results: A novel Delta-Omicron SARS-CoV-2 recombinant was identified in an unvaccinated, immunosuppressed kidney transplant recipient treated with monoclonal antibody Sotrovimab. The recombination breakpoint is located in the spike N-terminal domain, adjacent to the Sotrovimab quaternary binding site, and results in a 5'-Delta AY.45 and a 3'-Omicron BA.1 mosaic spike protein. Delta and BA.1 are sensitive to Sotrovimab neutralization, whereas the Delta-Omicron recombinant is highly resistant to Sotrovimab, both with and without the RBD resistance mutation E340D. Conclusions: Recombination between circulating SARS-CoV-2 variants can functionally contribute to immune escape. It is critical to validate phenotypes of mosaic viruses and monitor immunosuppressed COVID-19 patients treated with monoclonal antibodies for the selection of recombinant and immune escape variants. (Funded by NYU, the National Institutes of Health, and others).
RESUMO
Healthcare workers (HCWs) are known to be at higher risk of developing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections although whether these risks are equal across all occupational roles is uncertain. Identifying these risk factors and understand SARS-CoV-2 transmission pathways in healthcare settings are of high importance to achieve optimal protection measures. We aimed to investigate the implementation of a voluntary screening program for SARS-CoV-2 infections among hospital HCWs and to elucidate potential transmission pathways though phylogenetic analysis before the vaccination era. HCWs of the University Hospital of Liège, Belgium, were invited to participate in voluntary reverse transcriptase-polymerase chain reaction (RT-PCR) assays performed every week from April to December 2020. Phylogenetic analysis of SARS-CoV-2 genomes were performed for a subgroup of 45 HCWs. 5095 samples were collected from 703 HCWs. 212 test results were positive, 15 were indeterminate, and 4868 returned negative. 156 HCWs (22.2%) tested positive at least once during the study period. All SARS-CoV-2 test results returned negative for 547 HCWs (77.8%). Nurses (p < 0.05), paramedics (p < 0.05), and laboratory staff handling respiratory samples (p < 0.01) were at higher risk for being infected compared to the control non-patient facing group. Our phylogenetic analysis revealed that most positive samples corresponded to independent introduction events into the hospital. Our findings add to the growing evidence of differential risks of being infected among HCWs and support the need to implement appropriate protection measures based on each individual's risk profile to guarantee the protection of both HCWs and patients. Furthermore, our phylogenetic investigations highlight that most positive samples correspond to distinct introduction events into the hospital.
Assuntos
COVID-19 , Bélgica/epidemiologia , COVID-19/diagnóstico , COVID-19/epidemiologia , Atenção à Saúde , Pessoal de Saúde , Hospitais Universitários , Humanos , Recursos Humanos em Hospital , Filogenia , SARS-CoV-2/genéticaRESUMO
An adequate SARS-CoV-2 genomic surveillance strategy has proven to be essential for countries to obtain a thorough understanding of the variants and lineages being imported and successfully established within their borders. During 2020, genomic surveillance in Belgium was not structurally implemented but performed by individual research laboratories that had to acquire the necessary funds themselves to perform this important task. At the start of 2021, a nationwide genomic surveillance consortium was established in Belgium to markedly increase the country's genomic sequencing efforts (both in terms of intensity and representativeness), to perform quality control among participating laboratories, and to enable coordination and collaboration of research projects and publications. We here discuss the genomic surveillance efforts in Belgium before and after the establishment of its genomic sequencing consortium, provide an overview of the specifics of the consortium, and explore more details regarding the scientific studies that have been published as a result of the increased number of Belgian SARS-CoV-2 genomes that have become available.
Assuntos
COVID-19 , Pandemias , Humanos , Bélgica/epidemiologia , COVID-19/epidemiologia , Genoma Viral , Genômica , SARS-CoV-2/genética , Sequenciamento de Nucleotídeos em Larga EscalaRESUMO
As with other pathogens, competitive interactions between Bordetella pertussis strains drive infection risk. Vaccines are thought to perturb strain diversity through shifts in immune pressures; however, this has rarely been measured because of inadequate data and analytical tools. We used 3344 sequences from 23 countries to show that, on average, there are 28.1 transmission chains circulating within a subnational region, with the number of chains strongly associated with host population size. It took 5 to 10 years for B. pertussis to be homogeneously distributed throughout Europe, with the same time frame required for the United States. Increased fitness of pertactin-deficient strains after implementation of acellular vaccines, but reduced fitness otherwise, can explain long-term genotype dynamics. These findings highlight the role of vaccine policy in shifting local diversity of a pathogen that is responsible for 160,000 deaths annually.
Assuntos
Bordetella pertussis , Coqueluche , Bordetella pertussis/genética , Europa (Continente) , Genótipo , Humanos , Vacina contra Coqueluche , Coqueluche/epidemiologia , Coqueluche/prevenção & controleRESUMO
Advances in sequencing technologies have tremendously reduced the time and costs associated with sequence generation, making genomic data an important asset for routine public health practices. Within this context, phylogenetic and phylogeographic inference has become a popular method to study disease transmission. In a Bayesian context, these approaches have the benefit of accommodating phylogenetic uncertainty, and popular implementations provide the possibility to parameterize the transition rates between locations as a function of epidemiological and ecological data to reconstruct spatial spread while simultaneously identifying the main factors impacting the spatial spread dynamics. Recent developments enable researchers to make use of travel history data of infected individuals in the reconstruction of pathogen spread, offering increased inference accuracy and mitigating sampling bias. Here, we describe a detailed workflow to reconstruct the spatial spread of a pathogen through Bayesian phylogeographic analysis in discrete space using these novel approaches, implemented in BEAST. The individual protocols focus on how to incorporate molecular data, covariates of spread, and individual travel history data into the analysis. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Creating a SARS-CoV-2 MSA using sequences from GISAID Basic Protocol 2: Setting up a discrete trait phylogeographic reconstruction in BEAUti Basic Protocol 3: Phylogeographic reconstruction incorporating travel history information Basic Protocol 4: Visualizing ancestral spatial trajectories for specific taxa.
Assuntos
COVID-19/epidemiologia , COVID-19/virologia , SARS-CoV-2/genética , Viagem/estatística & dados numéricos , Teorema de Bayes , COVID-19/genética , COVID-19/transmissão , Biologia Computacional/métodos , Bases de Dados de Ácidos Nucleicos , Humanos , Filogenia , Filogeografia/métodos , SARS-CoV-2/classificação , SARS-CoV-2/isolamento & purificação , Análise de Sequência de DNA/métodos , Software , Estados Unidos/epidemiologiaRESUMO
At the end of 2020, several new variants of SARS-CoV-2-designated variants of concern-were detected and quickly suspected to be associated with a higher transmissibility and possible escape of vaccine-induced immunity. In Belgium, this discovery has motivated the initiation of a more ambitious genomic surveillance program, which is drastically increasing the number of SARS-CoV-2 genomes to analyse for monitoring the circulation of viral lineages and variants of concern. In order to efficiently analyse the massive collection of genomic data that are the result of such increased sequencing efforts, streamlined analytical strategies are crucial. In this study, we illustrate how to efficiently map the spatio-temporal dispersal of target mutations at a regional level. As a proof of concept, we focus on the Belgian province of Liège that has been consistently sampled throughout 2020, but was also one of the main epicenters of the second European epidemic wave. Specifically, we employ a recently developed phylogeographic workflow to infer the regional dispersal history of viral lineages associated with three specific mutations on the spike protein (S98F, A222V and S477N) and to quantify their relative importance through time. Our analytical pipeline enables analysing large data sets and has the potential to be quickly applied and updated to track target mutations in space and time throughout the course of an epidemic.
Assuntos
Genoma Viral , Mutação , SARS-CoV-2/genética , Glicoproteína da Espícula de Coronavírus/genética , Bélgica , Monitoramento Epidemiológico , HumanosRESUMO
Following the first wave of SARS-CoV-2 infections in spring 2020, Europe experienced a resurgence of the virus starting late summer that was deadlier and more difficult to contain. Relaxed intervention measures and summer travel have been implicated as drivers of the second wave. Here, we build a phylogeographic model to evaluate how newly introduced lineages, as opposed to the rekindling of persistent lineages, contributed to the COVID-19 resurgence in Europe. We inform this model using genomic, mobility and epidemiological data from 10 West European countries and estimate that in many countries more than 50% of the lineages circulating in late summer resulted from new introductions since June 15th. The success in onwards transmission of these lineages is predicted by SARS-CoV-2 incidence during this period. Relatively early introductions from Spain into the United Kingdom contributed to the successful spread of the 20A.EU1/B.1.177 variant. The pervasive spread of variants that have not been associated with an advantage in transmissibility highlights the threat of novel variants of concern that emerged more recently and have been disseminated by holiday travel. Our findings indicate that more effective and coordinated measures are required to contain spread through cross-border travel.