Genomic surveillance reveals dynamic shifts in the connectivity of COVID-19 epidemics.

The maturation of genomic surveillance in the past decade has enabled tracking of the emergence and spread of epidemics at an unprecedented level. During the COVID-19 pandemic, for example, genomic data revealed that local epidemics varied considerably in the frequency of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage importation and persistence, likely due to a combination of COVID-19 restrictions and changing connectivity. Here, we show that local COVID-19 epidemics are driven by regional transmission, including across international boundaries, but can become increasingly connected to distant locations following the relaxation of public health interventions. By integrating genomic, mobility, and epidemiological data, we find abundant transmission occurring between both adjacent and distant locations, supported by dynamic mobility patterns. We find that changing connectivity significantly influences local COVID-19 incidence. Our findings demonstrate a complex meaning of "local" when investigating connected epidemics and emphasize the importance of collaborative interventions for pandemic prevention and mitigation.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic.

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations.

Emergence of an early SARS-CoV-2 epidemic in the United States.

The emergence of the COVID-19 epidemic in the United States (U.S.) went largely undetected due to inadequate testing. New Orleans experienced one of the earliest and fastest accelerating outbreaks, coinciding with Mardi Gras. To gain insight into the emergence of SARS-CoV-2 in the U.S. and how large-scale events accelerate transmission, we sequenced SARS-CoV-2 genomes during the first wave of the COVID-19 epidemic in Louisiana. We show that SARS-CoV-2 in Louisiana had limited diversity compared to other U.S. states and that one introduction of SARS-CoV-2 led to almost all of the early transmission in Louisiana. By analyzing mobility and genomic data, we show that SARS-CoV-2 was already present in New Orleans before Mardi Gras, and the festival dramatically accelerated transmission. Our study provides an understanding of how superspreading during large-scale events played a key role during the early outbreak in the U.S. and can greatly accelerate epidemics.

Emergence and rapid transmission of SARS-CoV-2 B.1.1.7 in the United States.

The highly transmissible B.1.1.7 variant of SARS-CoV-2, first identified in the United Kingdom, has gained a foothold across the world. Using S gene target failure (SGTF) and SARS-CoV-2 genomic sequencing, we investigated the prevalence and dynamics of this variant in the United States (US), tracking it back to its early emergence. We found that, while the fraction of B.1.1.7 varied by state, the variant increased at a logistic rate with a roughly weekly doubling rate and an increased transmission of 40%-50%. We revealed several independent introductions of B.1.1.7 into the US as early as late November 2020, with community transmission spreading it to most states within months. We show that the US is on a similar trajectory as other countries where B.1.1.7 became dominant, requiring immediate and decisive action to minimize COVID-19 morbidity and mortality.

Transmission, infectivity, and neutralization of a spike L452R SARS-CoV-2 variant.

We identified an emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California, a state in the western United States. Named B.1.427/B.1.429 to denote its two lineages, the variant emerged in May 2020 and increased from 0% to >50% of sequenced cases from September 2020 to January 2021, showing 18.6%-24% increased transmissibility relative to wild-type circulating strains. The variant carries three mutations in the spike protein, including an L452R substitution. We found 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation common to variants B.1.1.7, B.1.351, and P.1. Antibody neutralization assays revealed 4.0- to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California exhibiting decreased antibody neutralization warrants further investigation.

Genomic Epidemiology of SARS-CoV-2 in Guangdong Province, China.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2 infection and was first reported in central China in December 2019. Extensive molecular surveillance in Guangdong, China's most populous province, during early 2020 resulted in 1,388 reported RNA-positive cases from 1.6 million tests. In order to understand the molecular epidemiology and genetic diversity of SARS-CoV-2 in China, we generated 53 genomes from infected individuals in Guangdong using a combination of metagenomic sequencing and tiling amplicon approaches. Combined epidemiological and phylogenetic analyses indicate multiple independent introductions to Guangdong, although phylogenetic clustering is uncertain because of low virus genetic variation early in the pandemic. Our results illustrate how the timing, size, and duration of putative local transmission chains were constrained by national travel restrictions and by the province's large-scale intensive surveillance and intervention measures. Despite these successes, COVID-19 surveillance in Guangdong is still required, because the number of cases imported from other countries has increased.

Coast-to-Coast Spread of SARS-CoV-2 during the Early Epidemic in the United States.

The novel coronavirus SARS-CoV-2 was first detected in the Pacific Northwest region of the United States in January 2020, with subsequent COVID-19 outbreaks detected in all 50 states by early March. To uncover the sources of SARS-CoV-2 introductions and patterns of spread within the United States, we sequenced nine viral genomes from early reported COVID-19 patients in Connecticut. Our phylogenetic analysis places the majority of these genomes with viruses sequenced from Washington state. By coupling our genomic data with domestic and international travel patterns, we show that early SARS-CoV-2 transmission in Connecticut was likely driven by domestic introductions. Moreover, the risk of domestic importation to Connecticut exceeded that of international importation by mid-March regardless of our estimated effects of federal travel restrictions. This study provides evidence of widespread sustained transmission of SARS-CoV-2 within the United States and highlights the critical need for local surveillance.

Travel Surveillance and Genomics Uncover a Hidden Zika Outbreak during the Waning Epidemic.

The Zika epidemic in the Americas has challenged surveillance and control. As the epidemic appears to be waning, it is unclear whether transmission is still ongoing, which is exacerbated by discrepancies in reporting. To uncover locations with lingering outbreaks, we investigated travel-associated Zika cases to identify transmission not captured by reporting. We uncovered an unreported outbreak in Cuba during 2017, a year after peak transmission in neighboring islands. By sequencing Zika virus, we show that the establishment of the virus was delayed by a year and that the ensuing outbreak was sparked by long-lived lineages of Zika virus from other Caribbean islands. Our data suggest that, although mosquito control in Cuba may initially have been effective at mitigating Zika virus transmission, such measures need to be maintained to be effective. Our study highlights how Zika virus may still be "silently" spreading and provides a framework for understanding outbreak dynamics. VIDEO ABSTRACT.

Inference of Infectious Disease Transmission through a Relaxed Bottleneck Using Multiple Genomes Per Host.

In recent times, pathogen genome sequencing has become increasingly used to investigate infectious disease outbreaks. When genomic data is sampled densely enough amongst infected individuals, it can help resolve who infected whom. However, transmission analysis cannot rely solely on a phylogeny of the genomes but must account for the within-host evolution of the pathogen, which blurs the relationship between phylogenetic and transmission trees. When only a single genome is sampled for each host, the uncertainty about who infected whom can be quite high. Consequently, transmission analysis based on multiple genomes of the same pathogen per host has a clear potential for delivering more precise results, even though it is more laborious to achieve. Here, we present a new methodology that can use any number of genomes sampled from a set of individuals to reconstruct their transmission network. Furthermore, we remove the need for the assumption of a complete transmission bottleneck. We use simulated data to show that our method becomes more accurate as more genomes per host are provided, and that it can infer key infectious disease parameters such as the size of the transmission bottleneck, within-host growth rate, basic reproduction number, and sampling fraction. We demonstrate the usefulness of our method in applications to real datasets from an outbreak of Pseudomonas aeruginosa amongst cystic fibrosis patients and a nosocomial outbreak of Klebsiella pneumoniae.

Global spread characteristics of CTX-M-type extended-spectrum ß-lactamases: A genomic epidemiology analysis.

BACKGROUND: Extended-spectrum ß-lactamases (ESBLs) producing bacteria have spread worldwide and become a global public health concern. Plasmid-mediated transfer of ESBLs is an important route for resistance acquisition. METHODS: We collected 1345 complete sequences of plasmids containing CTX-Ms from public database. The global transmission pattern of plasmids and evolutionary dynamics of CTX-Ms have been inferred. We applied the pan-genome clustering based on plasmid genomes and evolution analysis to demonstrate the transmission events. FINDINGS: Totally, 48 CTX-Ms genotypes and 186 incompatible types of plasmids were identified. The geographical distribution of CTX-Ms showed significant differences across countries and continents. CTX-M-14 and CTX-M-55 were found to be the dominant genotypes in Asia, while CTX-M-1 played a leading role in Europe. The plasmids can be divided into 12 lineages, some of which forming distinct geographical clusters in Asia and Europe, while others forming hybrid populations. The Inc types of plasmids are lineage-specific, with the CTX-M-1_IncI1-I (Alpha) and CTX-M-65_IncFII (pHN7A8)/R being the dominant patterns of cross-host and cross-regional transmission. The IncI-I (Alpha) plasmids with the highest number, were presumed to form communication groups in Europe-Asia and Asia-America-Oceania, showing the transmission model as global dissemination and regional microevolution. Meanwhile, the main kinetic elements of blaCTX-Ms showed genotypic preferences. ISEcpl and IS26 were most frequently involved in the transfer of CTX-M-14 and CTX-M-65, respectively. IS15 has become a crucial participant in mediating the dissemination of blaCTX-Ms. Interestingly, blaTEM and blaCTX-Ms often coexisted in the same transposable unit. Furthermore, antibiotic resistance genes associated with aminoglycosides, sulfonamides and cephalosporins showed a relatively high frequency of synergistic effects with CTX-Ms. CONCLUSIONS: We recognized the dominant blaCTX-Ms and mainstream plasmids of different continents. The results of this study provide support for a more effective response to the risks associated with the evolution of blaCTX-Ms-bearing plasmids, and lay the foundation for genotype-specific epidemiological surveillance of resistance, which are of important public health implications.

Spatial scale of tuberculosis transmission in Lima, Peru.

Spatially targeted interventions may be effective alternatives to individual or population-based prevention strategies against tuberculosis (TB). However, their efficacy may depend on the mechanisms that lead to geographically constrained hotspots. Local TB incidence may reflect high levels of local transmission; conversely, they may point to frequent travel of community members to high-risk areas. We used whole-genome sequencing to explore patterns of TB incidence and transmission in Lima, Peru. Between 2009 and 2012, we recruited incident pulmonary TB patients and their household contacts, whom we followed for the occurrence of TB disease. We used whole-genome sequences of 2,712 Mycobacterial tuberculosis isolates from 2,440 patients to estimate pariwise genomic distances and compared these to the spatial distance between patients' residences. Genomic distances increased rapidly as spatial distances increased and remained high beyond 2 km of separation. Next, we divided the study catchment area into 1 × 1 km grid-cell surface units and used household spatial coordinates to locate each TB patient to a specific cell. We estimated cell-specific transmission by calculating the proportion of patients in each cell with a pairwise genomic distance of 10 or fewer single-nucleotide polymorphisms. We found that cell-specific TB incidence and local transmission varied widely but that cell-specific TB incidence did not correlate closely with our estimates of local transmission (Cohen's k = 0.27). These findings indicate that an understanding of the spatial heterogeneity in the relative proportion of TB due to local transmission may help guide the implementation of spatially targeted interventions.

Genomic Epidemiology of Treponema pallidum and Circulation of Strains With Diminished tprK Antigen Variation Capability in Seattle, 2021-2022.

BACKGROUND: The incidence of syphilis continues to increase in the United States, yet little is known about Treponema pallidum genomic epidemiology within American metropolitan areas. METHODS: We performed whole-genome sequencing and tprK deep sequencing of 28 T. pallidum-containing specimens, collected mostly from remnant Aptima swab specimens from 24 individuals from Seattle Sexual Health Clinic during 2021-2022. RESULTS: All 12 individuals infected with Nichols-lineage strains were men who have sex with men, while a specific SS14 cluster (mean, 0.33 single-nucleotide variant) included 1 man who has sex with women and 5 women. All T. pallidum strains sequenced were azithromycin resistant via 23S ribosomal RNA A2058G mutation. Identical T. pallidum genomic sequences were found in pharyngeal and rectal swab specimens taken concurrently from the same individuals. The tprK sequences were less variable between patient-matched specimens and between epidemiologically linked clusters. We detected a 528-base pair deletion in the tprK donor site locus, eliminating 9 donor sites, in T. pallidum genomes of 3 individuals with secondary syphilis, associated with diminution of TprK diversity. CONCLUSIONS: We developed an end-to-end workflow for public health genomic surveillance of T. pallidum from remnant Aptima swab specimens. tprK sequencing may assist in linking cases beyond routine T. pallidum genome sequencing. T. pallidum strains with deletions in tprK donor sites currently circulate and are associated with diminished TprK antigenic diversity.

Genomic epidemiology and evolution of rhinovirus in western Washington State, 2021-22.

BACKGROUND: Human rhinoviruses (RV) primarily cause the common cold, but infection outcomes vary from subclinical to severe cases, including asthma exacerbations and fatal pneumonia in immunocompromised individuals. To date, therapeutic strategies have been hindered by the high diversity of serotypes. Global surveillance efforts have traditionally focused on sequencing VP1 or VP2/VP4 genetic regions, leaving gaps in our understanding of RV genomic diversity. METHODS: We sequenced 1,078 RV genomes from nasal swabs of symptomatic and asymptomatic individuals to explore viral evolution during two epidemiologically distinct periods in Washington State: when the COVID-19 pandemic affected the circulation of other seasonal respiratory viruses except for RV (February - July 2021), and when the seasonal viruses reemerged with the severe RSV and influenza outbreak (November-December 2022). We constructed maximum likelihood and BEAST-phylodynamic trees to characterize intra-genotype evolution. RESULTS: We detected 99 of 168 known genotypes and observed inter-genotypic recombination and genotype cluster swapping from 2021 to 2022. We found a significant association between the presence of symptoms and viral load, but not with RV species or genotype. Phylodynamic trees, polyprotein selection pressure, and Shannon entropy revealed co-circulation of divergent clades within genotypes with high amino acid constraints throughout polyprotein. DISCUSSION: Our study underscores the dynamic nature of RV genomic epidemiology within a localized geographic region, as more than 20% of existing genotypes within each RV species co-circulated each studied month. Our findings also emphasize the importance of investigating correlations between rhinovirus genotypes and serotypes to understand long-term immunity and cross-protection.

Analysis of a large SARS-CoV-2 (Alpha) outbreak in a Catalan prison using conventional and genomic epidemiology.

Enforcing strict protocols that prevent transmission of airborne infections in prisons is challenging. We examine a large SARS-CoV-2 outbreak in a Catalan penitentiary center from February-March 2021, prior to vaccination deployment. The aim was to describe the evolution of the outbreak using classical and genomic epidemiology and the containment strategy applied. The outbreak was initially detected in one module but spread to four, infecting 7 staff members and 140 incarcerated individuals, 6 of whom were hospitalized (4.4%). Genomic analysis confirmed a single origin (B.1.1.7). Contact tracing identified transmission vectors between modules and prevented further viral spread. In future similar scenarios, the control strategy described here may help limiting transmission of airborne infections in correctional settings.

Genomic Epidemiology of Large Blastomycosis Outbreak, Ontario, Canada, 2021.

Using phylogenomic analysis, we provide genomic epidemiology analysis of a large blastomycosis outbreak in Ontario, Canada, caused by Blastomyces gilchristii. The outbreak occurred in a locale where blastomycosis is rarely diagnosed, signaling a possible shift in geographically associated incidence patterns. Results elucidated fungal population genetic structure, enhancing understanding of the outbreak.

Recurrence, Microevolution, and Spatiotemporal Dynamics of Legionella pneumophila Sequence Type 1905, Portugal, 2014-2022.

We investigated molecular evolution and spatiotemporal dynamics of atypical Legionella pneumophila serogroup 1 sequence type 1905 and determined its long-term persistence and linkage to human disease in dispersed locations, far beyond the large 2014 outbreak epicenter in Portugal. Our finding highlights the need for public health interventions to prevent further disease spread.

Emergence of Poultry-Associated Human Salmonella enterica Serovar Abortusovis Infections, New South Wales, Australia.

Salmonella enterica serovar Abortusovis is a ovine-adapted pathogen that causes spontaneous abortion. Salmonella Abortusovis was reported in poultry in 2009 and has since been reported in human infections in New South Wales, Australia. Phylogenomic analysis revealed a clade of 51 closely related isolates from Australia originating in 2004. That clade was genetically distinct from ovine-associated isolates. The clade was widespread in New South Wales poultry production facilities but was only responsible for sporadic human infections. Some known virulence factors associated with human infections were only found in the poultry-associated clade, some of which were acquired through prophages and plasmids. Furthermore, the ovine-associated clade showed signs of genome decay, but the poultry-associated clade did not. Those genomic changes most likely led to differences in host range and disease type. Surveillance using the newly identified genetic markers will be vital for tracking Salmonella Abortusovis transmission in animals and to humans and preventing future outbreaks.

Coronavirus GenBrowser for monitoring the transmission and evolution of SARS-CoV-2.

Genomic epidemiology is important to study the COVID-19 pandemic, and more than two million severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic sequences were deposited into public databases. However, the exponential increase of sequences invokes unprecedented bioinformatic challenges. Here, we present the Coronavirus GenBrowser (CGB) based on a highly efficient analysis framework and a node-picking rendering strategy. In total, 1,002,739 high-quality genomic sequences with the transmission-related metadata were analyzed and visualized. The size of the core data file is only 12.20 MB, highly efficient for clean data sharing. Quick visualization modules and rich interactive operations are provided to explore the annotated SARS-CoV-2 evolutionary tree. CGB binary nomenclature is proposed to name each internal lineage. The pre-analyzed data can be filtered out according to the user-defined criteria to explore the transmission of SARS-CoV-2. Different evolutionary analyses can also be easily performed, such as the detection of accelerated evolution and ongoing positive selection. Moreover, the 75 genomic spots conserved in SARS-CoV-2 but non-conserved in other coronaviruses were identified, which may indicate the functional elements specifically important for SARS-CoV-2. The CGB was written in Java and JavaScript. It not only enables users who have no programming skills to analyze millions of genomic sequences, but also offers a panoramic vision of the transmission and evolution of SARS-CoV-2.

Genomic epidemiology and molecular characteristics of blaNDM-1-positive carbapenem-resistant Pseudomonas aeruginosa belonging to international high-risk clone ST773 in the Gauteng region, South Africa.

PURPOSE: The emergence of carbapenem-resistant P. aeruginosa (CRPA) harbouring acquired carbapenemase genes (blaVIM, blaIMP and blaNDM) has become a global public health threat. Three CRPA isolates included in the study had an extensively drug-resistant phenotype with susceptibility to colistin only and were positive for the blaNDM-1 gene. The current study aimed to investigate the genomic epidemiology and molecular characteristics of the blaNDM-1-positive CRPA isolates collected from the Gauteng region, South Africa. METHODS: Short read whole genome sequencing (WGS) was performed to determine sequence types (STs), genetic relatedness, resistome, virulome and the genetic environment of the blaNDM-1 gene. RESULTS: The WGS and phylogenetic analyses revealed that the study isolates belonged to an international high-risk clone ST773 and belonged to the same clade with eight blaNDM-1-positive ST773 isolates from Hungary, India, Nigeria, South Korea and USA. The study isolates harboured a wide repertoire of intrinsic and acquired antibiotic resistance genes (ARGs) related with mobile genetic elements, porins and efflux pumps, as well as virulence factor genes. The clade-specific ARGs (blaNDM-1, floR2/cmlA9, rmtB4, tetG) were found in a putative integrative and conjugative element (ICE) region similar to ICE6660-like. CONCLUSION: As ICE carrying the blaNDM-1 gene can easily spread to other P. aeruginosa isolates and other Gram-negative bacteria, the findings in this study highlight the need for appropriate management strategies and active surveillance of CRPA isolates in the Gauteng region, South Africa.

Investigating the cause of a 2021 winter wave of COVID-19 in a border region in eastern Germany: a mixed-methods study, August to November 2021.

It is so far unclear how the COVID-19 winter waves started and what should be done to prevent possible future waves. In this study, we deciphered the dynamic course of a winter wave in 2021 in Saxony, a state in Eastern Germany neighbouring the Czech Republic and Poland. The study was carried out through the integration of multiple virus genomic epidemiology approaches to track transmission chains, identify emerging variants and investigate dynamic changes in transmission clusters. For identified local variants of interest, functional evaluations were performed. Multiple long-lasting community transmission clusters have been identified acting as driving force for the winter wave 2021. Analysis of the dynamic courses of two representative clusters indicated a similar transmission pattern. However, the transmission cluster caused by a locally occurring new Delta variant AY.36.1 showed a distinct transmission pattern, and functional analyses revealed a replication advantage of it. This study indicated that long-lasting community transmission clusters starting since early autumn caused by imported or locally occurring variants all contributed to the development of the 2021 winter wave. The information we achieved might help future pandemic prevention.