SEARCH · Search PubMed

Results for “genomic epidemiology”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

The SARS-CoV-2 Integrated Genomic Epidemiology Database (IGED): Linking viral genomes with patient-level metadata to advance statewide genomic surveillance in California.

In July 2021, the California Code of Regulations Title 17 required all laboratories performing SARS‑CoV‑2 whole genome sequencing (WGS) to report their sequencing results to the California Department of Public Health (CDPH). These viral genomic data and patient metadata were compiled into the Integrated Genomic Epidemiology Database (IGED). Linking anonymized viral sequences with patient‑level information enabled monitoring of infectiousness, pathogenicity, transmission dynamics, evolution, and vaccine evasion among emerging SARS‑CoV‑2 lineages. Laboratories performing SARS-CoV-2 WGS transmitted sequencing results to CDPH through Electronic Laboratory Reporting (ELR) and non-ELR pathways. CDPH applied uniform reporting requirements but allowed flexibility in specific data formats to accommodate diverse data systems. To preserve data quality and interoperability across heterogeneous sources, CDPH implemented standardization, validation, and deduplication protocols. Snowflake, a cloud‑based data storage and analytics platform, and Posit Connect, a cloud deployment and automation platform, supported the management, processing, and integration of data within the IGED. The IGED established links between SARS‑CoV‑2 WGS data and epidemiologic metadata for 801,418 sequences, representing 81.7% of all sequences reported in California. Lineages reported to the IGED showed strong concordance with lineage proportions in GISAID. Sequences reported to the IGED had average turnaround times longer than one month, and the majority of sequencing was performed in Southern California and Los Angeles. The IGED enhanced genomic surveillance through predictive modeling and monitoring concerning evolutionary trends such as recombination and saltations in persistent infections. Development of the IGED highlighted the need for standardized data requirements, sustained funding for sequencing, incentives for data submission, and interdisciplinary collaboration to build an effective genomic surveillance system. This framework for linking genomic and epidemiologic data has not only generated critical insights for SARS‑CoV‑2 but also provided the foundation for CDPH and other public health organizations to develop similar IGED‑like systems for other priority pathogens as genomic surveillance expands.

Journal Article↗

PubMed2025

Invasive Streptococcus dysgalactiae subspecies equisimilis compared with Streptococcus pyogenes in Australia, 2011-23, and the emergence of a multi-continent stG62647 lineage: a retrospective clinical and genomic epidemiology study.

BACKGROUND: Streptococcus dysgalactiae subspecies equisimilis (SDSE) is closely related to Streptococcus pyogenes, with overlapping disease manifestations. We compared the clinical and genomic epidemiology of invasive SDSE with invasive S pyogenes across different settings in Australia and phylogenetically contextualised the SDSE sequences within a global cohort of genomes. METHODS: In this retrospective clinical and genomic epidemiology study, cases of invasive SDSE isolated from normally sterile sites were identified and whole-genome sequenced across five hospital networks in temperate southeast Australia (Melbourne and Sydney) and the tropical Top End of the Northern Territory. SDSE disease incidence, case demographics, clinical outcomes, and longitudinal lineage dynamics were compared between southeast Australia and the Top End and to co-collected invasive S pyogenes cases in each region. SDSE genomes and lineages were also contextualised within 1166 global SDSE sequences. Genomic transmission clusters (not necessarily direct transmission) were inferred between isolates from different individuals by single-linkage clustering at a single nucleotide polymorphism threshold of less than or equal to seven for SDSE and less than or equal to five for S pyogenes based on previous transmission analyses. FINDINGS: Between Jan 1, 2011, and Feb 28, 2023, there were 693 invasive SDSE cases and 995 invasive S pyogenes cases. Invasive SDSE occurred almost exclusively in adults. The overall invasive SDSE incidence in southeast Australia was similar to invasive S pyogenes (incidence rate ratio [IRR] 1·15, 95% CI 0·91-1·46; p=0·26) and increased over the study period (IRR 1·06 per year, 95% CI 1·05-1·08; p<0·0001) from 1·30 cases per 10 000 admissions in 2011 to 3·72 cases per 10 000 admissions in the first 2 months of 2023 (95% CI 2·13-6·07). In southeast Australia, where stringent COVID-19 non-pharmaceutical interventions (NPIs) were implemented between 2020 and 2021, the SDSE incidence plateaued during 2020-21 but did not significantly decline (IRR 1·09 compared with 2017-19, 95% CI 0·88-1·35; p=0·47). By contrast, S pyogenes incidence substantially declined in 2020-21 in southeast Australia (IRR 0·35 compared to 2017-19, 95% CI 0·22-0·52; p=0·017). In the Top End, SDSE incidence was lower than S pyogenes (IRR 0·24, 95% CI 0·19-0·31; p<0·0001). However, crude incidence remained higher than southeast Australia (crude IRR 1·24, 95% CI 1·07-1·42; p=0·0037) and disproportionately affected First Nations Australians in the Top End compared with non-First Nations individuals (IRR 3·36, 95% CI 2·33-4·85; p<0·0001). Comparing 2020-21 with 2017-19, there was no decline in SDSE (IRR 1·27, 95% CI 0·73-2·24; p=0·45) or S pyogenes (IRR 0·97, 95% CI 0·80-1·18; p=0·81) incidence in the Top End, which did not implement prolonged stringent COVID-19 NPIs. Analysing the available genomes of invasive cases and in lineages for which more than or equal to five invasive cases occurred, only 24 (6%) of 384 SDSE cases were assigned to genomic transmission clusters, compared with 271 (52%) of 524 S pyogenes cases. An stG62647 lineage encompassed 113 (26%) of 436 sequenced SDSE genomes. Analysis of available SDSE sequences from Australia, western Europe, and North America inferred concurrent international expansion of the stG62647 lineage in all three regions between 1990 and 2005. INTERPRETATION: We identified a substantial burden of invasive SDSE, dominated by the emergent stG62647 lineage. The contrasting epidemiology between species in the different Australian regions, during COVID-19 NPIs, and genomic infection patterns indicates transmission dynamic, pathogen population, and host-pathogen interaction differences between SDSE and S pyogenes and indicates implications for disease control measures. FUNDING: Australian National Health and Medical Research Council.

Humans↗

PubMed2025

Global diversity and evolution of Salmonella enterica serovar Panama: a genomic epidemiology study.

BACKGROUND: Non-typhoidal Salmonella is a globally important bacterial pathogen, typically associated with foodborne gastrointestinal infection. Some non-typhoidal Salmonella serovars can also colonise typically sterile sites in people to cause invasive non-typhoidal Salmonella disease. Salmonella enterica serovar Panama is responsible for a substantial number of cases of human bloodstream infection, but despite its global dissemination, numerous outbreaks, and a reported association with invasive non-typhoidal Salmonella disease, S enterica serovar Panama (S Panama) is understudied. We aimed to describe the genomic epidemiology and evolutionary history of S Panama to provide a vital baseline of understanding for this globally important serovar. METHODS: In this genomic epidemiology study, we analysed S Panama genomes derived from historical collections, national surveillance datasets, and publicly available epidemiological and whole-genome sequencing data which span the years 1931-2019. Maximum likelihood and Bayesian phylodynamic approaches were used to investigate population structure and evolutionary history and to infer geotemporal dissemination. A combination of different bioinformatic approaches with short-read and long-read data were used to characterise geographical and clade-specific trends in antimicrobial resistance (AMR) and genetic markers for invasiveness. FINDINGS: We analysed 836 S Panama genomes, of which 559 (67%) were sequenced as part of this study. The collection represents all inhabited continents and includes isolates collected between 1931 and 2019. We identified the presence of four geographically linked S Panama clades (C1 [ie, the Latin America and the Caribbean clade; n=338], C2 [ie, the European clade; n=124], C3 [ie, the Martinique clade; n=131], and C4 [ie, the Asia and Oceania clade; n=104]) and regional trends in AMR profiles. Most isolates (715 [86%] of 836) were pan-susceptible to antibiotics and belonged to clades circulating in Latin America and the Caribbean (64%, n=458). Most antibiotic-resistant isolates in our collection (113 [93%] of 121) fell within clades C4 (ie, the Asia and Oceania clade) and C2 (ie, the European clade), the latter of which had the highest invasiveness index values based on the conservation of 196 extraintestinal predictor genes. INTERPRETATION: This first large-scale phylogenetic analysis of S Panama has revealed important information about the population structure, AMR, global ecology, and genetic markers of invasiveness of the identified genomic subtypes. Our findings provide an important baseline for understanding S Panama infection. The presence of multidrug-resistant clades with elevated invasiveness index values should be monitored through ongoing surveillance, as such clades could pose an increased public health risk. FUNDING: UK Research and Innovation Global Challenges Research Fund and Biotechnology and Biological Sciences Research Council, UK Medical Research Council, Wellcome Trust, John Lennon Memorial Scholarship, Institut Pasteur, Santé publique France, Fondation Le Roch-Les Mousquetaires, Investissement d'Avenir Programme, and Australian National Health and Medical Research Council.

Humans↗

PubMed2026

Emerging terbinafine-resistant Trichophyton indotineae between 2018 and 2023: a multinational genomic epidemiology study.

BACKGROUND: Trichophyton species cause the greatest burden of dermatophytosis worldwide, with the Trichophyton mentagrophytes species complex being particularly associated with the emergence of new aggressive infections. One emerging species, Trichophyton indotineae is notable for its clinical resistance to terbinafine antifungal treatment and rapid global spread. In this study we aim to characterise the epidemiology of this emerging pathogen using genomics. METHODS: In this retrospective genomic epidemiology study, to better understand the epidemiology of this disease, we sourced isolates collected from patients with severe cases of dermatophytosis (identified either by internal transcribed spacer sequencing or phenotypic characterisation) in the UK, Ireland, France, Canada, and India for the period 2014-23, including the T indotineae type strain from Japan. We used whole-genome sequencing to confirm 90 isolates were T indotineae, and antifungal susceptibility testing to assess susceptibility to terbinafine. FINDINGS: 103 cases of severe dermatophytosis caused by Trichophyton species collected between 2018 and 2023 in the UK, France, Canada, Ireland, and India were included in this study. Susceptibility testing indicated that 63 (70%) of 90 T indotineae isolates were resistant to terbinafine (minimum inhibitory concentration [MIC] ≥0·5 mg/L). Pairwise genetic distances showed very high identity with only 147 (range 1-414) single-nucleotide polymorphisms (SNPs) separating isolates that were nested within a monophyletic phylogeny, supporting a single evolutionary origin of T indotineae. Genome-wide analyses identified multiple non-synonymous SNPs in SQLE (ERG1), the squalene epoxidase target of terbinafine, that were associated with terbinafine in vitro resistance of 0·5 mg/L or higher. However, six isolates exhibited high MIC values without SQLE mutations, suggesting the presence of alternative resistance mechanisms. INTERPRETATION: That no clear geographical clustering of isolates was observed confirms the rapid transcontinental spread of T indotineae from its likely centre of diversity in Asia. Our findings highlight the importance of better genomic surveillance to understand and manage this severe and rapidly emerging terbinafine-resistant dermatophyte. FUNDING: None.

Terbinafine↗

PubMed2026

Effect of repeated mass drug administration on the transmission of yaws: a retrospective genomic epidemiology study.

BACKGROUND: Yaws, a neglected tropical disease caused by Treponema pallidum subspecies pertenue (T p pertenue), has evaded eradication, in part due to a high proportion of asymptomatic cases. Repeated mass drug administration (MDA), whereby an entire population is repeatedly treated irrespective of disease, could provide a solution. Here, we aimed to investigate the effect of MDA on the genomic epidemiology of T p pertenue. METHODS: We conducted a retrospective genomic epidemiology study on samples collected during a cluster-randomised trial of mass administration of azithromycin for yaws eradication in the Namatanai District of Papua New Guinea. Participants were in 38 wards (administrative units encompassing several villages) in three local-level government areas (LLGs). The experimental group received an initial round of MDA followed by two further rounds 6 months and 12 months after the first round. The control group received one round of MDA followed by two rounds of treatment targeting clinical cases and contacts only, on the same schedule as the MDA in the experimental group. A follow-up survey on both groups was done 18 months after the first MDA round. Swab samples were collected at each round from ulcerative and nodular skin lesions, and blood was collected by finger-prick for serological testing at 18 months. Metadata on ulcer size (cm) and duration (days) were recorded at each round, and treponemal and non-treponemal antibodies were recorded at 18 months. Samples from swabs positive for T p pertenue underwent library preparation and whole-genome sequencing. We examined the phylogenetic relationships between genomes, linking them with geospatial and patient metadata to understand the impact of MDA on T p pertenue diversity and transmission. FINDINGS: Swabs collected from 297 individuals with active yaws from April 30, 2018, to Nov 2, 2019, yielded 222 good-quality Tp pertenue genomes. We identified 20 sublineages of T p pertenue in the control group and 21 in the experimental group at the beginning of the study. At the end of the study, there were 13 sublineages in the control group and three in the experimental group, of which two persisted in both groups. Three sublineages not detected at baseline were observed in the control group after commencing MDA. The two sublineages that persisted in both groups had non-synonymous mutations in penicillin-binding proteins. One of these sublineages evolved macrolide resistance in three individuals and was associated with lowered treponemal antibody (p=0·0036) and longer ulcer duration (p=0·015). Despite the study taking place within a small island, sublineages were geographically clustered, with pairs of samples from the same ward (odds ratio 7·1, 95% CI 5·7-8·8; p<0·0001) or neighbouring wards (4·3, 3·3-5·4; p<0·0001) more likely to share the same sublineages compared with pairs from different LLGs. Additionally, older individuals were more likely to share sublineages than were younger individuals (1·5, 1·2-1·9; p<0·0001). INTERPRETATION: Repeated MDA was successful in reducing and maintaining the genetic diversity of T p pertenue at a low level but was associated with the development of macrolide resistance. Yaws re-emergence after MDA was attributed to multiple sublineages, of which the majority were detected in the population before MDA. Participants within the same ward were more likely to share sublineages than those that were more widely geographically separated, suggesting that re-emergence was driven by local transmission. These findings could inform future yaws elimination strategies. FUNDING: European Research Council, EU, Provincial Deputation of Barcelona, Barberà Solidària Foundation, Wellcome, and Fundació "la Caixa".

Adolescent↗

PubMed2025

Genomic Epidemiology of the Main SARS-CoV-2 Variants Circulating in Italy During the Omicron Era.

Since early 2022 the Omicron variant has rapidly spread worldwide, becoming the dominant variant to date. The study aimed to investigate the clinical and epidemiological characteristics of COVID-19 patients and reconstruct the genomic epidemiology of main SARS-CoV-2 Omicron sublineages in Italy in 2022. A total of 8970 SARS-CoV-2 samples were studied, and phylogenetic analyses were focused on BA.1, BA.2, and BA.5 subvariants. More than half of subjects received three doses of vaccine and experienced a reinfection. A significant larger proportion of unvaccinated subjects presented reinfection compared with vaccinated. Clusters presented a tMRCA between September-November 2021 (BA.1), November 2021-January 2022 (BA.2), and October 2021-May 2022 (BA.5). Re values showed the highest level between September-October, January-February 2022, and May 2022 for BA.1, BA.2 and BA.5, respectively. Limited number of studied variant sequences are included in clusters. The spread rate of the studied variant exceeded its evolutionary rate. No single sublineage had sufficient time to differentiate into large clusters, but only into small and fragmented groups sharing the same recent ancestor. These analyses dissect the epidemiological dynamics of Omicron sublineages in Italy over a period of great epidemiological changes in the COVID-19 epidemic.

Humans↗

PubMed2026

Clinical Outcomes and Genomic Epidemiology of Multidrug-Resistant Methicillin-Resistant Staphylococcus aureus Keratitis.

PURPOSE: To characterize the clinical features, management, antimicrobial resistance patterns, and genomic epidemiology of methicillin-resistant Staphylococcus aureus (MRSA) keratitis at two North American centers. DESIGN: Retrospective interventional case series combined with laboratory investigation PARTICIPANTS: Seventy eyes of 67 patients presenting laboratory-confirmed MRSA keratitis were included METHODS: We performed a multicenter retrospective case series of patients with culture-proven MRSA keratitis treated between 2005 and 2022. Demographic and clinical data were collected. Antimicrobial susceptibility testing was conducted, and multidrug resistance (MDR) was defined as resistance to ≥3 antibiotic classes. A subset of isolates underwent whole-genome sequencing with core genome multilocus sequence typing. Vancomycin susceptibility, heteroresistance screening, and tolerance testing were performed on available isolates. MAIN OUTCOME MEASURES: Antimicrobial susceptibility and multidrug resistance rates, vancomycin phenotypic profiles, MRSA genotypic distribution, and final best-corrected visual acuity RESULTS: Median age was 63.5 years, and 61.4% were female. Ocular surface disease (67.7%) and prior ocular surgery (65.2%) were common. Only 25.4% had significant healthcare exposure in the preceding year. Most isolates (85.7%) were MDR. Fluoroquinolone susceptibility was low (moxifloxacin 19.7%). All isolates were susceptible to vancomycin (MIC₉₀ 2 µg/mL), and no vancomycin-intermediate, heteroresistant, or tolerant phenotypes were identified. Whole genome sequencing (n = 41) demonstrated predominance of clonal complexes 5 (68.3%) and 8 (29.2%). Visual outcomes were poor, with most patients (85.2%) having a final visual acuity worse than 20/60 among those with follow-up. CONCLUSIONS: MRSA keratitis is associated with high rates of multidrug resistance and poor visual outcomes despite guideline-based therapy. Infections were predominantly caused by CC5 MDR strains despite limited recent healthcare exposure. These findings highlight the persistence of highly resistant MRSA lineages in community-associated corneal infection and underscore the need for ongoing antimicrobial surveillance and optimized treatment strategies.

Humans↗

PubMed2026

Temporal shifts in gyrA mutation types and sublineage replacement in ST11 Salmonella enterica serovar Enteritidis over a decade (2014-2023): A genomic epidemiological study in Guangxi, China.

The overuse or abuse of antibiotics drives the global health threat of antimicrobial resistance. Although bans on certain veterinary antibiotics, such as colistin, have proven effective, the impact of fluoroquinolone stewardship on the evolution of the foodborne pathogen Salmonella enterica serovar Enteritidis (S. Enteritidis) remains unclear. Here, we conducted a decade-long (2014-2023) retrospective longitudinal genomic epidemiological analysis of 441 ST11 S. Enteritidis isolates from Guangxi, China, alongside a global reference dataset of 4297 genomes. Our aim was to elucidate the effect of real-world antibiotic stewardship on the shift of gyrA point mutations and lineage distribution. Surveillance identified three global epidemic clade sublineages (GEC-L2, L3, L4), with the multidrug-resistant GEC-L4 (i.e., GC-c or MMC2), characterized by the gyrA mutation with amino acid substitution D87Y, being domestically dominant (70.07%, 309/441). Following China's 2016 ban on the veterinary use of critical fluoroquinolones, the proportion of the highly resistant GEC-L4 sublineage decreased continuously (from 86.84% in 2017 to 56.00% in 2023), while the less resistant GEC-L3 sublineage (i.e., GC-b or MMC1), mainly characterized by gyrA D87G, increased simultaneously (from 13.16% to 44.00%). This phenomenon might be attributed to the fact that the GEC-L4 sublineage exhibited a higher fitness cost compared with the GEC-L3 sublineage, as confirmed by the competition assay. A Random Forest Model validated that the gyrA mutation with amino acid substitution D87Y was the paramount feature for these sublineages' identification. In contrast, global data showed a continuous increase in gyrA mutations (from 8.63% in 2006 to 68.85% in 2024), primarily D87Y (from 1.44% to 31.15%) and D87N (from 4.32% to 22.95%), correlating with rising average fluoroquinolone consumption. This study provides direct genomic evidence that national-level antibiotic stewardship can drive the replacement of highly resistant sublineages with moderately resistant ones. These findings offer crucial scientific evidence for evaluating the impact of antibiotic management policies and inform strategies for the rational use of antimicrobials.

China↗

PubMed2026

Comparison of phylogenetic metrics of transmission between symptomatic and asymptomatic tuberculosis in individuals who were incarcerated in Brazil in 2008-24: a retrospective genomic epidemiology study.

BACKGROUND: Tuberculosis control efforts have traditionally targeted symptomatic individuals; however, the role of asymptomatic cases in sustaining transmission is increasingly recognised. We aimed to quantify the contribution of asymptomatic tuberculosis to recent transmission using genomic and epidemiological data from a high-transmission setting. METHODS: We conducted a retrospective genomic epidemiology study of Mycobacterium tuberculosis isolates collected in Mato Grosso do Sul, Brazil, between Aug 25, 2008, and March 19, 2024. Available isolates underwent whole-genome sequencing. Demographic, clinical, incarceration history, and laboratory metadata were obtained from surveillance records. From Jan 1, 2017, to March 19, 2024, active case finding was conducted in the state's three largest prisons (all male-only facilities), during which sputum samples were collected from individuals irrespective of symptoms and tested using GeneXpert and culture. Comparisons of transmission between individuals with and without symptoms were restricted to individuals who were incarcerated and were identified through active case finding and for whom high-quality, M tuberculosis lineage 4 genomes were available. Metrics of recent transmission included phylogenetic clustering, time-scaled haplotype density (THD), local branching index (LBI), and transmission probabilities inferred using Bayesian Reconstruction and Evolutionary Analysis of Transmission Histories. FINDINGS: 4448 tuberculosis cases were notified in Mato Grosso do Sul in 2008-24. After excluding cases for which M tuberculosis isolates were not available or had low sequencing quality, who had contaminated cultures or mixed infection, or who were infected with non-lineage 4 M tuberculosis, we included 2362 lineage 4 M tuberculosis isolates with high-quality genome sequences. 1849 (78·3%) of 2362 isolates were part of a genomic cluster. Among 2362 individuals with tuberculosis, 1137 (48·1%) were incarcerated at diagnosis. Of these individuals, 505 were identified through active case finding in three male-only prisons. The median age was 30 years (IQR 25-37); 304 (60·2%) had mixed ethnicity, 90 (17·8%) were White, 56 (11·1%) were Black, 13 (2·6%) were Indigenous, and six (1·2%) were Asian. 277 (54·9%) had symptomatic disease and 228 (45·1%) had asymptomatic tuberculosis. There were no significant differences between symptomatic and asymptomatic individuals in phylogenetic clustering (213 [76·9%] of 277 vs 195 [85·5%] of 228; p=0·37), THD (median 0·39 [IQR 0·06-0·62] vs 0·50 [0·09-0·65]; p=0·12), or LBI (0·00863 [0·00810-0·00988] vs 0·00871 [0·00829-0·01020]; p=0·088). Bayesian transmission trees showed no significant difference in the number of secondary infections inferred from symptomatic compared with asymptomatic individuals (p=0·56). These findings were consistent across genomic clusters and robust to model assumptions. INTERPRETATION: We identified no differences in transmission between individuals who were symptomatic and those who were asymptomatic using multiple genomic measures. In this high-transmission setting, where systematic screening is implemented, our findings indicate that asymptomatic tuberculosis substantially contributes to tuberculosis transmission at the population level. These results suggest that symptom-based case detection alone is likely to be insufficient to interrupt transmission and highlight the importance of expanded screening strategies in high-risk populations. FUNDING: US National Institutes of Health and the Brazilian National Research Council (CNPq).

Humans↗

PubMed2026

Genomic epidemiology of dengue virus 2 and 3 reveals repeated introductions and exportations of several lineages in Colombia.

Dengue fever, a major mosquito-borne viral disease, is transmitted by Aedes mosquitoes and poses a significant global health burden. Despite extensive research, the spatiotemporal dynamics of dengue virus (DENV) lineages in Colombia remain understudied. Here, we analyze 11,443 complete genome sequences from Colombia and the Americas to map the genomic epidemiology of DENV-2 and DENV-3. Phylogeographic reconstruction revealed multiple independent introductions and exportations of the DENV-2 II and III lineages, as well as the DENV-3 lineage III_C.2, underscoring Colombia's critical role as both a source and a sink of viral traffic within the Americas. Antigenic profiling demonstrated distinct clustering of emergent lineages in antigenic space, consistent with immune-escape-driven turnover. These results highlight the necessity of sustained, high-resolution genomic surveillance to guide targeted public health interventions and mitigate dengue transmission across the region.

Dengue Virus↗

PubMed2025

Genomic epidemiology of clade Ia monkeypox viruses circulating in the Central African Republic in 2022-24: a retrospective cross-sectional study.

BACKGROUND: The spread of monkeypox virus (Orthopoxvirus monkeypox) clade Ib from the Democratic Republic of the Congo to neighbouring countries has raised global concerns, leading to WHO declaring mpox a public health emergency on Aug 14, 2024. We applied genomic epidemiology to investigate the causes of recurrent mpox outbreaks in the Central African Republic. We aimed to determine whether frequent zoonotic spillovers or increased human-to-human transmissions are driving mpox epidemiology. METHODS: We performed a retrospective cross-sectional study of monkeypox virus genomic sequences among PCR-confirmed mpox cases detected in the Central African Republic between Feb 17, 2022, and Sept 17, 2024. We used hybridisation capture coupled to high throughput sequencing to analyse 46 samples from mpox outbreaks that occurred in eight of the 20 prefectures (14 of 35 health districts). Near-complete genomes were used for phylogenomic analyses. FINDINGS: Between Jan 10, 2022, and Sept 15, 2024, 89 mpox cases were confirmed, including 53 cases in the first 9 months of 2024. We generated 41 near-complete genomes from this period, including 33 from 2024. All new and already published monkeypox virus genomes from the Central African Republic belonged to clade Ia. These genomes spanned the phylogenetic diversity of clade Ia viruses, and most likely represented several dozen independent transmission events to humans. The monkeypox virus phylogenetic diversity was geographically structured within the country. Plausibly linked cases often showed indistinguishable genomes. Conversely, we detected identical genomes in cases that epidemiological information would suggest were independent outbreaks. Finally, we found that three distinct viruses caused cases in the capital city of Bangui in July, 2024, with all three detected on the same day (July 24, 2024). We did not detect substantial enrichment of APOBEC3 editing, suggesting limited human-to-human transmission. INTERPRETATION: The data indicate that mpox epidemiology in the Central African Republic is primarily driven by short-lived outbreaks resulting from many independent zoonotic spillover events, particularly in rural areas. Although evidence remains limited, in Bangui additional factors such as movement of people and importation of bushmeat from other regions might be introducing the virus into urban settings. Similar spillover patterns have been observed in the Democratic Republic of the Congo. The poorly understood nature of monkeypox virus reservoirs in both countries is a regional concern, as frequent spillovers increase the risk of outbreaks leading to sustained human transmission. Beyond strengthening surveillance and developing countermeasures, it is important to better understand the reservoirs and focus on reducing transmission opportunities to prevent further outbreaks. FUNDING: Pasteur Institute of Bangui, Africa CDC, AFROSCREEN, WHO, the Helmholtz Institute for One Health, and the Deutsche Forschungsgemeinschaft.

Humans↗

PubMed2026

First nationwide full-genome characterisation of human-derived Andes virus in Chile: a retrospective genomic epidemiology study.

BACKGROUND: Andes virus (ANDV) is the only hantavirus known to transmit between humans and causes hantavirus cardiopulmonary syndrome in Chile and Argentina. In Chile, ANDV genomic diversity remains incompletely characterised. This study aimed to characterise the genetic diversity, geographical structure, and molecular signatures of ANDV using human clinical samples collected over a 13-year period (2011-24). METHODS: We conducted a retrospective genomic epidemiology study of ANDV infections in Chile. Clinical samples from patients with confirmed ANDV, collected between March 9, 2011, and June 27, 2024, were analysed and sequenced. Clinical and epidemiological data were obtained from diagnostic laboratories and surveillance programmes. Consensus sequences for the S, M, and L segments were generated, and genetic clustering and divergence were assessed using phylogenetic inference and variant calling. FINDINGS: We analysed clinical samples from 58 infected individuals and identified two major genomic variants of ANDV with distinct geographical distributions, defined by regionally structured patterns of nucleotide and amino acid substitutions across the S, M, and L segments: ANDV Chi-North (central Chile) and ANDV-South (southern Chile). No consistent clustering by clinical severity was observed, and no recurrent non-synonymous substitutions were uniquely associated with severe disease. Substitutions previously associated with person-to-person transmission in outbreaks in Argentina were not consistently observed in Chilean sequences, including in four person-to-person transmission cases. Although some substitutions described in ANDV-like viruses were present in the Chi-North lineage, this lineage remained phylogenetically distinct and geographically restricted to central Chile. INTERPRETATION: To our knowledge, this study provides the first nationwide genomic characterisation of human-derived ANDV in Chile. The identification of geographically structured variants indicates that ANDV diversity in Chile is driven by regional diversification rather than clinical outcome. The absence of consistent amino acid signatures associated with disease severity or person-to-person transmission suggests that these phenotypes are unlikely to be explained by viral genetic variation alone. These findings refine current understanding of ANDV evolution and highlight the need for continued integrated genomic surveillance in endemic regions. FUNDING: Agencia Nacional de Investigación y Desarrollo de Chile and National Institutes of Health.

Humans↗

PubMed2026

Characterisation of Bordetella pertussis virulence and macrolide resistance in Australia by targeted culture-independent sequencing: a genomic epidemiology study.

BACKGROUND: Bordetella pertussis continues to circulate globally despite widespread vaccination, with a notable epidemic in 2024. Its resurgence is confounded by the emergence of pertactin-deficient, macrolide-resistant B pertussis strains in Asia and Europe, which are under-recognised by conventional diagnostics. We aimed to apply targeted culture-independent next-generation sequencing (tNGS) of respiratory specimens to improve global B pertussis diagnostic capability and genomic surveillance. METHODS: We did a nationwide genomic epidemiology study of B pertussis RT-PCR-positive respiratory specimens that were retrospectively and prospectively collected by diagnostic and public health laboratories in six of seven states and territories of Australia. Specimens underwent tNGS and macrolide-resistant B pertussis-specific PCR, and an opportunistic subset from New South Wales and Queensland were cultured for confirmatory susceptibility testing and whole-genome sequencing. Sequencing data were analysed for genome recovery, virulence profiles, and macrolide resistance mutations, and were compared with international macrolide-resistant B pertussis genomes and ancestral Australian genomes. The performance of the tNGS approach was assessed with logistic regression relative to RT-PCR cycle threshold values, and sensitivity and specificity values were calculated. FINDINGS: 255 respiratory specimens positive for B pertussis were included in the study. 64 (25%) were retrospectively collected between Jan 12, 2012, and Dec 31, 2023, and 191 (75%) were prospectively collected between Jan 1 and Oct 28, 2024. Of these 255 specimens, 148 (58%) yielded near-complete B pertussis genomes through tNGS. Seven co-circulating lineages of B pertussis were documented, including two associated with macrolide-resistance. Eight epidemiologically unrelated and geographically dispersed cases of macrolide-resistant B pertussis with a 23S rRNA 2037A→G mutation were identified by tNGS and confirmed by whole-genome sequencing. Three of these were further validated by phenotypic testing. The estimated prevalence of macrolide resistance among Australian cases positive for B pertussis was 4% (eight of 188). INTERPRETATION: tNGS can recover near-complete B pertussis genomes directly from clinical specimens, enabling identification of macrolide resistance mutations and high-resolution phylogenetic analysis. These findings show that tNGS complements PCR-based surveillance by providing genome-wide assessment of resistance, virulence, and genomic diversity in a single workflow. FUNDING: NSW Health Prevention Research Support Program.

Macrolides↗

PubMed2025

Global emergence and transmission dynamics of carbapenemase-producing Citrobacter freundii sequence type 22 high-risk international clone: a retrospective, genomic, epidemiological study.

BACKGROUND: Carbapenemase-producing Citrobacter (CPC) species have recently been recognised as emerging pathogens associated with nosocomial infections in humans. The increased rate of Citrobacter freundii infections is a public health concern and there is a paucity of genomic data regarding its global transmission dynamics. We aimed to characterise the genetic features of CPC species, and their associated carbapenemase-encoding plasmids, obtained from hospitalised patients in China and from publicly available global data, with a particular focus on high-risk clones. METHODS: This was a retrospective, genomic epidemiological study of CPC species obtained from a tertiary hospital in Zhejiang Province, China, from March 5, 2013, to March 5, 2023. We used antimicrobial susceptibility testing, short-read and long-read whole-genome sequencing, phylogenomic analysis, and plasmid structure analysis. A global dataset of complete plasmid sequences encoding blaKPC, blaNDM, and blaIMP was constructed from the National Center for Biotechnology Information (NCBI) RefSeq database to provide insights into their diversity and distribution. All carbapenemase-producing Citrobacter freundii genomes from the NCBI GenBank database were incorporated in the comparative genomic analyses. Bayesian phylogeographical analysis and growth rate assays were carried out to characterise the high-risk C freundii sequence type (ST) 22 clone. FINDINGS: 1724 Citrobacter species isolates were collected from diverse clinical specimens, with 48 identified as CPC species. Citrobacter koseri (22 [46%] of 48) and C freundii (20 [42%]) were the predominant CPC species. Comparative analysis found C freundii carried significantly higher median numbers of plasmid replicons (5·0 [IQR 3·3-6·0] vs 2·0 [2·0-3·0]; p<0·0001) and acquired antimicrobial resistance genes (12·0 [7·3-15·8] vs 3·0 [3·0-5·3]; p<0·0001) than did C koseri. Molecular characterisation identified Inc-type plasmids, In823::Kl.pn.I3/In1589-like/In837-like integrons, Tn6296/Tn125/Tn5060 transposons, and insertion sequences (eg, IS26, IS3000, IS5, ISAba125, ISCR1), collectively facilitating the dissemination of carbapenemase genes. Global analysis of 3126 carbapenemase-encoding plasmids found epidemic plasmids with broad host ranges and global diversity. Phylogenetic investigation of predominant carbapenemase-encoding plasmids showed their persistence across geographical regions, temporal spans, and Enterobacterales species, exhibiting high genetic similarity to our clinical plasmids. A phylogenetic tree of 726 global carbapenemase-producing C freundii genomes showed that ST22 (227 [31·3%]) represents the predominant multidrug-resistant clone across community, health-care, and environmental niches. Transmission across continents contributes to the global predominance of the ST22 clone, which carries a high load of resistance genes (median 15·0 [IQR 11·0-17·0] vs 12·0 [3·0-16·0]; p<0·0001) and enhanced plasmid maintenance capacity (median replicons 5·0 [IQR 4·0-7·0] vs 4·0 [3·0-6·0]; p<0·0001) relative to non-ST22 clones. INTERPRETATION: Our study provides evidence to suggest that Citrobacter species are emerging carriers of carbapenem-resistance genes. These findings provide insight into the population structure of CPC species and highlight C freundii ST22 as a prominent high-risk international clone. FUNDING: National Natural Science Foundation of China, National Health Commission Scientific Research Fund-Zhejiang Provincial Major Health Science and Technology Plan Project, Zhejiang Province Natural Science Foundation Project, Outstanding Youth Foundation of Jiangsu Province of China, the Priority Academic Program Development of Jiangsu Higher Education Institutions, and Postgraduate Research and Practice Innovation Program of Jiangsu Province.

Citrobacter freundii↗

PubMed2025

Temporal and geographical lineage dynamics of invasive Streptococcus pyogenes in Australia from 2011 to 2023: a retrospective, multicentre, clinical and genomic epidemiology study.

BACKGROUND: Defining the temporal dynamics of invasive Streptococcus pyogenes (group A Streptococcus) and differences between hyperendemic and lower-incidence regions provides crucial insights into pathogen evolution and, in turn, informs preventive measures. We aimed to examine the clinical and temporal lineage dynamics of S pyogenes across different disease settings in Australia to improve understanding of drivers of pathogen diversity. METHODS: In this retrospective, multicentre, clinical and genomic epidemiology study, we identified cases of invasive S pyogenes infection from normally sterile sites between Jan 1, 2011, and Feb 28, 2023. Data were collected from five hospital networks across low-incidence regions in temperate southeast Australia and the hyperendemic, tropical, and largely remote Top End of the Northern Territory of Australia. The crude incidence rate ratio (IRR) of bloodstream S pyogenes infection comparing the Top End and southeast Australia and in First Nations people compared with non-First Nations people was estimated by quasi-Poisson regression. We estimated odds ratios (ORs) of intensive care unit (ICU) admission, in-hospital mortality, and 30-day mortality for the Top End versus southeast Australia using logistic regression. Retrieved and successfully sequenced isolates were assigned lineages at whole-genome resolution. Temporal trends in the composition of co-circulating lineages were compared between the two regions. We used an S pyogenes-specific multistrain simulated transmission model to examine the relationship between host population-specific parameters and observed pathogen lineage dynamics. The prevalence of accessory genes (those present in 5-95% of all genomes) was compared across geographies and temporal periods to investigate genomic drivers of diversity. FINDINGS: We identified 500 cases of invasive S pyogenes infection in patients in the Top End and 495 cases in patients in southeast Australia. The crude IRR of bloodstream infection for the Top End compared with southeast Australia was 5·97 (95% CI 4·61-7·73) across the entire study period; in the Top End, infection disproportionately affected First Nations people compared with non-First Nations people (5·41, 4·28-6·89). The odds of in-hospital mortality (OR 0·43, 95% CI 0·26-0·70), 30-day mortality (0·38, 0·23-0·63), and ICU admission (0·42, 0·30-0·59) were lower in the Top End than in southeast Australia. Longitudinal lineage analysis of 642 S pyogenes genomes identified waves of replacement with distinct lineages in the Top End, whereas southeast Australia had a small number of dominant lineages that persisted and cycled in frequency. The transmission model qualitatively reproduced a similar pattern of replacement with distinct lineages when using a high transmission rate, small population size, and high levels of human movement-characteristics similar to those of communities in the hyperendemic Top End. Using a lower transmission rate, larger population size, and lower levels of migration similar to those of communities in urbanised southeast Australia, the transmission model qualitatively reproduced a pattern of dominant lineages that cycled in frequency. Despite distinct circulating lineages, the prevalence of accessory genes in the bacterial population was maintained across geographies and temporal periods. INTERPRETATION: In a hyperendemic setting, the replacement of distinct S pyogenes lineages occurred in waves, which could be linked to the disproportionate burden of disease and sparse human population in this setting. The maintenance of bacterial gene frequency could be consistent with multilocus selection. These findings suggest that lineage-specific interventions-such as vaccines under development-should consider disease setting and, without broad cross-protection, might lead to lineage replacement. FUNDING: National Health and Medical Research Council, and Leducq Foundation.

Humans↗

PubMed2026

Genomic epidemiology of clinically critical antibiotic resistance in Salmonella enterica causing bloodstream infections across six Chinese provinces, 1994-2023.

Clinically critical antibiotic-resistant Salmonella enterica (S. enterica) causing bloodstream infections remains a public health challenge. Here, we aim to reveal the emergence and trends of clinically important antibiotic resistance in S. enterica causing bloodstream infections using 833 isolates from six Chinese provincial-level administrative areas during 1994-2023. We identified 48 serovars and 64 sequence types (STs). Overall, 8.52% of 833 isolates were resistant or had decreased susceptibility to ciprofloxacin, 4.32% and 6.84% reported resistance or decreased susceptibility to third- and fourth-generation cephalosporins (3GCs and 4GCs), 1.80% reported resistance to fosfomycin, and 2.16% reported resistance to azithromycin. Across these six regions, azithromycin and fosfomycin resistance is increasing, as is decreased susceptibility or resistance to ciprofloxacin, 3GCs, and 4GCs, especially among younger children and elderly people. Clinically prioritized antibiotic resistance also varies by region, serovar, and age group. S. Paratyphi A genotype 2.3.3 strains are mainly divided into 2 lineages distributed in Guangxi and Shanghai. Within the scope of this passive surveillance dataset, S. Typhi genotype 4.3.1.2.1 was identified as the earliest documented case among the collected isolates. Our retrospective and longitudinal genomic epidemiology study provides critical data for the formulation of treatment guidelines and policies for bloodstream infections and for the monitoring and control of antimicrobial resistance.

Humans↗

PubMed2026

Genomic epidemiology and antimicrobial resistance profile of Shigella isolated from diarrhoea diseases in under-five children in Blantyre, Malawi.

Antimicrobial resistance (AMR) in Shigella is rising globally, complicating shigellosis management. Whole-genome sequence analysis (WGSA) has advanced our understanding of AMR and transmission dynamics, yet contemporary whole-genome sequencing data from Shigella in sub-Saharan Africa remain scarce. In this study, we applied WGSA to 27 Shigella isolates collected from children presenting with diarrhoea at Ndirande Health Centre in Malawi (2022-2023), as part of the Enterics for Global Health Shigella surveillance study. Serotyping, AMR profiling and phylogenetic analysis revealed Shigella sonnei as the dominant serogroup, with distinct genetic clustering relative to global reference strains among S. sonnei, Shigella flexneri and Shigella boydii. We identified 16 AMR genes linked to ten antimicrobial classes with qnrS1 and qnrB19 genes conferring resistance to fluoroquinolone, alongside IncFIB(K) and IncFII plasmid replicon markers. Importantly, no azithromycin resistance determinants were detected both genotypically and phenotypically, providing baseline evidence that warrants continued surveillance of current first-line treatment. However, the detection of fluoroquinolone resistance genes with plasmid replicon markers in the absence of phenotypic resistance might indicate a silent reservoir with epidemic potential. This is the first contemporary Shigella data from a large-scale diarrhoea disease surveillance study in Malawi, providing essential baseline information for guiding antibiotic treatment and future vaccine development efforts, contributing to the efforts to combat shigellosis in Malawi and other similar regions.

Humans↗

PubMed2026

Genomic Epidemiology and Clinical Characteristics of Mpox Lineage C.1 Outbreak in Thailand, 2023-2024.

Since 2022, human monkeypox virus (hMPXV) has emerged in non-endemic regions, including Thailand. However, the genomic dynamics and clinical correlates of local transmission remain incompletely defined. Whole-genome sequencing was performed on hMPXV from 16 patients in Thailand (2023-2024) using targeted amplicon NGS. Phylogenetic analyses integrated global reference sequences. Mutational profiles, specifically non-synonymous substitutions and APOBEC3-associated signatures, were analyzed in relation to clinical data. Phylogenetic reconstruction identified three temporal phases. Early 2022 cases (clade IIb lineages A and B) were interspersed with global sequences, consistent with multiple introductions. In contrast, 2023-2024 cases were dominated by lineage C.1. All 16 genomes belonged to C.1 (one C.1.1), and formed a distinct mid-2023 cluster, designated C.1/Thai/Cluster, supporting sustained local transmission. APOBEC3-associated mutations were pervasive across the C.1 lineage overall, including within C.1/Thai/Cluster, without evidence of significant enrichment specific to this cluster. The cohort comprised exclusively male patients (81% HIV-positive, MSM), with predominantly genital painful lesions and a median recovery time of 23 days. No significant associations were detected between viral genetic variation and clinical outcomes. Mpox transmission in Thailand evolved from multiple introductions to sustained C.1-dominated local spread, underscoring the importance of continued genomic surveillance.

Humans↗