SEARCH · PubMed Health

Results for “phylogenetic comparative analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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 large mitochondrial genome of Syndiclis anlungensis (Lauraceae): Genome structure, comparative analysis, and phylogenetic relationships among Syndiclis species.

The complete mitochondrial genome (mitogenome) of Syndiclis anlungensis, a critically endangered tropical tree, was determined in this study. The mitogenome spans 2,368,454 bp across four contigs and harbors 41 protein-coding genes, 22 tRNA genes, and three rRNA genes. Potential mutation regions, including 1317 repeat sequences and 698 simple sequence repeats (SSRs), were accurately located in the S. anlungensis mitogenome. Sixty-five transferred fragments of the repeats were found between its mitochondrial and chloroplast genomes. When compared to three other Laurales mitogenomes, extensive gene order shuffling is evident, leaving only five conserved gene clusters intact. Codon usage analysis reveals a pronounced A/T bias in both mitochondrial and chloroplast genes, and three mitochondrial genes (atp9, rps19, and sdh3) stand out for their high divergence across eleven Syndiclis taxa. Selection analyses indicate strong purifying pressure on rpl2, rpl16, and sdh3 (Ka/Ks < 1), with no positive selection detected. Using 41 mitochondrial protein-coding gene sequences from sixteen and three individuals of Syndiclis and Beilschmiedia species, respectively, our phylogenetic tree recovers Syndiclis as monophyletic, with two well-supported clades: one includes S. anlungensis, S. chinensis, S. lotungensis, S. marlipoensis, and a putative new Syndiclis species from Yunnan; the other contains S. furfuracea, S. hongkongensis, S. kwangsiensis, and three putative new Syndiclis species from Guangdong and Vietnam.

Genome, Mitochondrial↗

PubMed2026

Chloroplast genome comparative analysis and phylogenetic relationships of 15 Syringa species (Oleaceae).

Syringa is a crucial shrub genus in the family Oleaceae, which has significant ornamental, economic, and medicinal value. However, research on the chloroplast genome (CPG) phylogeny and lineage diversification of this genus remains limited. In this study, all 15 Syringa CPGs exhibited a characteristic quadripartite structure, with genome lengths ranging from 154,019-158,020 bp. These CPGs were highly conserved and moderately differentiated, each containing 130-132 genes. Analysis of inverted repeat (IR) boundaries indicated structural conservation, with six genes: rps19, rpl2, ycf1, trnN, ndhF, and trnH present at the IR/single-copy (SC) junctions. The small single copy (SSC) region displayed greater sequence variability than the IR regions. ycf1, ndhH, trnL-rpl32, ndhF-ycf1, and rbcL-accD were identified as potential molecular markers and rps11, ycf2, and ycf4 may have contributed to the adaptive evolution of Syringa. Phylogenetic reconstruction based on whole CPG data supported the monophyly of the 15 species, which were divided into three distinct subclades. Molecular dating estimated that Syringa diverged from its sister genus approximately 58 million years ago, with most Syringa species diversifying further approximately 47.49 million years ago during the Eocene. Our findings will hopefully stimulate further studies on this genus that may enhance biodiversity knowledge.

Journal Article↗

PubMed2025

Molecular Evolution and Zoonotic Potential of Muju Virus (Orthohantavirus puumalaense) in Craseomys regulus, Republic of Korea.

Orthohantavirus puumalaense causes hemorrhagic fever with renal syndrome in Europe, with Puumala virus (PUUV) as its primary representative. Muju virus (MUJV), harbored by Craseomys regulus, an Arvicolinae rodent species endemic to the Republic of Korea (ROK), is also a genotype of O. puumalaense. However, their genomic diversity and zoonotic potential remain largely unknown. To investigate their prevalence, 185 voles were collected from 23 regions of the ROK between 2012 and 2023. Serological assays detected anti-PUUV immunoglobulin G antibodies in five samples (3.1%), whereas reverse-transcription polymerase chain reaction confirmed MUJV RNA in identical specimens (2.7%). Amplicon-based nanopore sequencing facilitates near-complete genome recovery, enabling high-resolution comparative analysis. Phylogenetic analysis revealed distinct genetic lineages in Gangwon and Jeollabuk Provinces. Evolutionary rate estimates indicated greater sequence divergence in the S and L segments than in the M segment. A zoonotic risk assessment revealed that most MUJV variants exhibited moderate-to-high spillover potential. The molecular detection of MUJV in Cheorwon, Gangwon Province, expands its known geographic range and provides the first molecular evidence of MUJV circulation in this region. These findings highlight the need for continued surveillance and seroprevalence studies of MUJV to assess its potential for human exposure and public health relevance in the ROK.

Animals↗

PubMed2026

Genomic and functional characterization of novel therapeutic lytic bacteriophages targeting multidrug-resistant Enterobacter cloacae.

The alarming rates at which extensively drug-resistant (XDR) and pandrug-resistant (PDR) Enterobacter cloacae in hospitals are increasing has begun to severely limit treatment options, and thus the urgency for alternative interventions, including bacteriophage therapy. The purpose of the study was to isolate and molecularly characterize phages that can infect E. cloacae, and, furthermore, to assess the antimicrobial efficacy of the four novel lytic bacteriophages (MMRP1, MMRP2, MMRP3, and MMRP4) against antimicrobial-resistant E. cloacae isolates and to evaluate their potential as alternative therapeutic strategies. These novel phages were characterized by plaque morphology, transmission electron microscopy (TEM), host range testing, thermal and chloroform stability assays, bacterial reduction assays, and whole-genome sequencing (WGS). Among 27 clinical isolates, MDR, XDR, and PDR phenotypes were observed in 20 (74.1%), six (22.2%), and one (3.7%) isolates, respectively. All four phages produced clear lytic plaques (0.5-3.0 mm) with titers reaching up to 6 × 1010 PFU/mL, and the phage cocktail lysed 81.4% (22 of 27 isolates) of clinical isolates with high host specificity. TEM revealed that all four E. cloacae-infecting phages (MMRP1, MMRP2, MMRP3, and MMRP4) belong to the class Caudoviricetes, exhibiting icosahedral capsids, tailed morphology, and double-stranded DNA genomes, consistent with current ICTV classification criteria. Whole genome sequencing and comparative phylogenetic analysis further resolved the taxonomic placement of these phages at the family level, positioning MMRP1 within the family Demerecviridae and MMRP4 within the family Straboviridae. All phages were stable from -20 to 40 °C and were unaffected by exposure to chloroform. Phage cocktail reduced bacterial OD₆₀₀ to ≤ 0.3 within 4 h in the bacterial reduction test. WGS revealed large circular dsDNA genomes of ~132 kbp (MMRP1) and ~149 kbp (MMRP4), GC content of 38%, and modular architectures encoding structural, lytic, and replication gene modules. The most striking and highlighted suggestion that in vitro evaluation of MMRP1 and MMRP4 are highly recommended to more deeper future experimental studies to combat MDR E. cloacae nosocomial infections supported by genomic foundation and eventually, the possibility to be suitable for phage-engineering applications in clinical settings.

Enterobacter cloacae↗

PubMed2026

Comparative genomics and phylogenetic analysis of three Malvaceae species on the basis of chloroplast genomes.

INTRODUCTION: The Malvaceae family shows rich species diversity and has substantial economic and medicinal value. However, the frequent interspecific hybridization among members of this family has resulted in confused phylogenetic relationships among the groups, limiting the usefulness of traditional classification methods. METHODS: This study aimed to investigate the phylogenetic relationships among selected taxa of Malvaceae by evaluating 23 chloroplast (CP) genomes, including three newly assembled CP genomes. Among these three genomes, the CP genome of Hibiscus schizopetalus L. was reported for the first time, while the CP genomes of Alcea rosea L. and Hibiscus grewiifolius L., which have been deposited in NCBI, were re-analyzed here alongside newly generated data for comparative purposes. In addition, 20 downloaded CP genomes encompassing 13 genera were analyzed using SNPs in whole CP genomes data. RESULTS: The results showed that the genomes ranged from 160,403 to 161,978 base pairs in length and consisted of small single copies (SSCs) and large single copies (LSCs) separated by two inverted repeat sequences (IRs), forming a typical quadripartite circular structure. The entire genome sequence showed relative conservation across species in terms of structure, GC content, codon usage, and gene composition. The mutation sites were mainly located in the LSC and SSC regions, and the variability in the non-coding regions was higher than that in the coding regions. The nucleotide polymorphism (Pi) analysis identified the non-coding regions such as ndhF-rpl32 and psbZ-trnG as high variable hotspots. A maximum likelihood phylogenetic tree was constructed based on SNPs in whole CP genomes data. The phylogenetic analysis divided these 23 species into five highly supported clades. It also revealed a close sister-group relationship between Abelmoschus and Hibiscus species, suggesting that Hibiscus may have a separate lineage from okra species. DISCUSSION: In conclusion, the increasing availability of CP genome resources will enhance our understanding of the classification and evolutionary patterns of the Malvaceae family. The development of molecular markers will provide important molecular evidence for precise identification and classification revision of plants in this family.

Malvaceae↗

PubMed2026

Transcriptome mining and comparative genomics reveal 36 putative novel marafivirus species and conserved evolution of the marafibox regulatory element.

BACKGROUND: Marafiviruses are plant-infecting RNA viruses associated with several economically important crops, but their genomic diversity remains incompletely characterized. OBJECTIVE: This study aimed to identify previously unrecognized marafivirus genomes and investigate their genomic features and evolutionary relationships. METHODS: Publicly available plant transcriptome datasets were systematically mined to detect marafivirus-like sequences. Recovered genomes were analyzed using comparative sequence analysis, phylogenetic reconstruction, and genome organization characterization. RESULTS: A total of 62 marafivirus-like genomes were recovered from 33 independent sources representing diverse plant hosts. Polyprotein-based comparative and phylogenetic analyses grouped these genomes into 36 lineages likely representing novel species. All newly identified viruses clustered within the Marafivirus clade. Genome organization analysis revealed conserved polyprotein architecture and widespread presence of the marafibox promoter element. Conservation of additional open reading frames among closely related isolates aided identification of potentially functional genes. CONCLUSION: These findings substantially expand the known diversity of marafiviruses and demonstrate the effectiveness of transcriptome mining for discovering previously unrecognized plant viruses.

Phylogeny↗

PubMed2026

Di-, tri-, and tetranucleotide frequencies covary with lifespan and genome size across protostome invertebrates.

Animal lifespans span orders of magnitude, yet how genome sequence covaries with lifespan remains poorly characterized outside vertebrates. Although promoter CpG density has been linked to vertebrate longevity due to its gene-regulatory function through DNA methylation, it is unclear whether such patterns are promoter- and CpG-specific, or if they reflect broader sequence evolution. We curated maximum lifespan estimates for 466 protostome species spanning eight phyla with available genome assemblies and quantified mono-, di-, tri-, and tetranucleotide composition across whole genomes, intergenic regions, and six gene-associated regions (two upstream regions, exons, introns, and two downstream regions) defined using Benchmarking Universal Single-Copy Orthologs. Dinucleotide observed/expected ratios showed significant associations with lifespan and genome size in different ways. Lifespan-associated motifs were most pronounced in gene-associated non-coding regions, especially in introns and downstream regions, whereas genome-size effects were strongest in whole-genome and intergenic sequence. Tri- and tetranucleotide observed/expected ratios broadly recapitulated this regional organization. In contrast, GC content was not associated with lifespan across regions, indicating that the observed signals are not explained by mononucleotide composition but instead by how those nucleotides are arranged into short sequence motifs. These results suggest that lifespan and genome size show distinct but overlapping associations with regional sequence composition across invertebrate species and that lifespan-associated motif evolution extends beyond vertebrate promoter methylation architectures.

CpG density↗

PubMed2024

Molecular epidemiology of enteroviruses from Guatemalan wastewater isolated from human lung fibroblasts.

The Global Specialized Polio Laboratory at CDC supports the Global Poliovirus Laboratory Network with environmental surveillance (ES) to detect the presence of vaccine strain polioviruses, vaccine-derived polioviruses, and wild polioviruses in high-risk countries. Environmental sampling provides valuable supplementary information, particularly in areas with gaps in surveillance of acute flaccid paralysis (AFP) mainly in children less than 15 years. In collaboration with Guatemala's National Health Laboratory (Laboratorio Nacional de Salud Guatemala), monthly sewage collections allowed screening enterovirus (EV) presence without incurring additional costs for sample collection, transport, or concentration. Murine recombinant fibroblast L-cells (L20B) and human rhabdomyosarcoma (RD) cells are used for the isolation of polioviruses following a standard detection algorithm. Though non-polio-Enteroviruses (NPEV) can be isolated, the algorithm is optimized for the detection of polioviruses. To explore if other EV's are present in sewage not found through standard methods, five additional cell lines were piloted in a small-scale experiment, and next-generation sequencing (NGS) was used for the identification of any EV types. Human lung fibroblast cells (HLF) were selected based on their ability to isolate EV-A genus. Sewage concentrates collected between 2020-2021 were isolated in HLF cells and any cytopathic effect positive isolates used for NGS. A large variety of EVs, including echoviruses 1, 3, 6, 7, 11, 13, 18, 19, 25, 29; coxsackievirus A13, B2, and B5, EV-C99, EVB, and polioviruses (Sabin 1 and 3) were identified through genomic typing in NGS. When the EV genotypes were compared by phylogenetic analysis, it showed many EV's were genomically like viruses previously isolated from ES collected in Haiti. Enterovirus occurrence did not follow a seasonality, but more diverse EV types were found in ES collection sites with lower populations. Using the additional cell line in the existing poliovirus ES algorithm may add value by providing data about EV circulation, without additional sample collection or processing. Next-generation sequencing closed gaps in knowledge providing molecular epidemiological information on multiple EV types and full genome sequences of EVs present in wastewater in Guatemala.

Humans↗

PubMed2026

Comparative Analysis of Chloroplast Genomes Reveals Molecular Evolution and Phylogenetic Relationships in Fraxinus (Fraxinus mandshurica).

Fraxinus mandshurica (Manchurian ash) is an ecologically and economically valuable hardwood tree native to Northeast Asia, yet its genomic resources remain limited. We assembled its complete chloroplast (cp) genome (155,559 bp) using hybrid PacBio and Illumina sequencing and performed comparative, phylogenetic, and evolutionary analyses. The cp genome exhibits a typical quadripartite structure encoding 132 gene copies, comprising 114 unique genes (80 protein-coding, 30 tRNA, and 4 rRNA genes), with 18 genes duplicated in the inverted repeat (IR) regions. Simple sequence repeat analysis revealed dominance of mononucleotide A/T repeats. Phylogenetic analysis of 53 complete cp genomes strongly supported the monophyly of Oleaceae and resolved F. mandshurica as sister to the North American F. nigra, consistent with previously proposed Miocene intercontinental dispersal scenarios between East Asia and North America. Most protein-coding genes were under strong purifying selection (Ka/Ks << 1), whereas petB, rpl2, and several ndh genes showed elevated Ka/Ks values that are suggestive of altered selective constraint but are based on very few substitutions and are therefore not, on their own, evidence of positive selection. Nucleotide diversity (Pi) analysis identified 15 hypervariable intergenic spacers (mean Pi = 0.067), among which trnM-CAU-rps14, ndhJ-ndhK, and petL-petG represent promising candidate barcode regions requiring further validation. This study provides a high-quality, fully annotated cp genome of F. mandshurica and a valuable genomic resource for future phylogenetic, population genetic, and conservation studies of this important genus.

Fraxinus↗

PubMed2026

Molecular epidemiology and genomic characteristics of clinical Acinetobacter baumannii isolates from patients with hospital-acquired pneumonia in China, 2019-2020: a multicentre retrospective study.

BACKGROUND: Carbapenem-resistant Acinetobacter baumannii (CRAB) is a leading cause of hospital-acquired pneumonia (HAP) with high mortality. However, large-scale nationwide data of HAP-causing CRAB in China remain limited. METHODS: Here, we performed a nationwide multicentre retrospective study to characterise the molecular epidemiology and genomic features of 802 A. baumannii isolates from patients with HAP across 33 tertiary hospitals in China during 2019-2020. Antimicrobial susceptibility testing (AST), whole-genome sequencing (WGS), phylogenetic and comparative genomic analysis were used to investigate molecular epidemiology of HAP-causing CRAB strain. Clinical comparative analyses were carried out on data from 500 patients with HAP stratified by distinct antimicrobial susceptibility and genomic profiles, and a Galleria mellonella infection model was utilised for in vivo virulence assessment. FINDINGS: The overall carbapenem resistance rate of A. baumannii was 82.0% (658/802), with marked regional variations. CRAB exhibited high resistance to conventional agents but remained largely susceptible to polymyxin, tigecycline, cefiderocol and sulbactam-durlobactam. Among enrolled patients, CRAB infection was linked to substantially higher mortality (39.0% vs. 17.5%), and multivariate analysis confirmed ICU admission and advanced age as independent risk factors for patients with CRAB infection. Molecular typing revealed STPas2 (96.2%) as the absolutely predominant type; STOxf208, STOxf195, STOxf540, and STOxf369 were the most prevalent Oxford sequence types with obvious geographic stratification and divergent comorbidity profiles among corresponding patients. A total of 654 CRAB isolates harboured carbapenemase genes, with blaOXA-23 dominating at 98.8%. Genomic analysis revealed lineage-specific features: STOxf208 carried more virulence genes, while STOxf540 harboured a broader antimicrobial resistance genes (ARGs). The STOxf208 clone mainly belonged to KL2 (62.1%) and KL7 (36.8%) serotypes, with KL2 strains possessing richer ARGs and virulence factors, and in vivo virulence assays further validated that KL2 strains possessed higher pathogenicity than KL7 strains. INTERPRETATION: This study demonstrates the extremely high prevalence and clonal dominance of CRAB in Chinese patients with HAP, providing critical evidence for clinical treatment, antimicrobial stewardship, and targeted infection control. FUNDING: National Key Research and Development Program of China (2024YFE0106200), National Natural Science Foundation of China (U22A20338, 82502763, W2621007), Zhejiang Provincial Natural Science Foundation of China (LQN25H190006), Zhejiang Provincial Postdoctoral Science Foundation (ZJ2025058).

Acinetobacter baumannii↗

PubMed2024

Comparative analysis of chloroplast genomes in ten holly (Ilex) species: insights into phylogenetics and genome evolution.

In order to clarify the chloroplast genomes and structural features of ten Ilex species and provide insights into the phylogeny and genome evolution of the genus Ilex, we conducted a comparative analysis of chloroplast genomes using bioinformatics methods. The chloroplast genomes of ten Ilex species were obtained, and their structural features and variations were compared. The results indicated that all chloroplast genomes in the genus Ilex exhibit a double-stranded circular structure, with sizes ranging from 157,356 to 158,018 bp, showing minimal differences in size. The chloroplast genomes of the ten Ilex species have a relatively conservative gene count, with a total of 134 to 135 genes, including 88 or 89 protein-coding genes, and a conserved number of 8 rRNA genes. Each chloroplast genome contains 3 to 123 SSR (Simple Sequence Repeat) sites, predominantly composed of mononucleotide and trinucleotide repeats, with no detection of pentanucleotide or hexanucleotide repeats. The variation in dispersed repeat sequences among Ilex species is minimal, with a total repeat sequence number ranging from 1 to 14, concentrated in the length range of 30 to 42 base pairs. The expansion and contraction of chloroplast genome boundaries among Ilex species are relatively stable, with only minor variations observed in individual species. Variations in non-coding regions are more pronounced than those in coding regions, with the variability in the Large Single Copy region (LSC) being the highest, while the variability in the Inverted Repeat region A (IRa) is the lowest. The divergence time among Ilex species was estimated using the MCMC-tree module, revealing the evolutionary relationships among these species, their common ancestors, and their differentiation throughout the evolutionary process. The research findings provide a valuable reference for the systematic study and molecular marker development of Ilex plants.

Genome, Chloroplast↗

PubMed2026

Chromosome-Level Genome Assembly and Annotation of the Chinese Lizard Gudgeon (Saurogobio dabryi).

The Chinese lizard gudgeon (Saurogobio dabryi) is an economically important freshwater species within the Cyprinidae family, abundant in the middle and lower reaches of the Yangtze River and its adjacent basins. As a promising species suitable for aquaculture in China, the lack of genomic resources has rendered the genetic breeding and conservation research. Here, we present the first chromosome-level genome assembly of S. dabryi using PacBio HiFi long reads, short reads, and Hi-C sequencing data. The final assembly reaches a total size of 1.09 Gb and Hi-C scaffolding anchors 99.55% of the assembled contigs onto 25 chromosomes, with a scaffold N50 reaching 43.15 Mb. The final genome assembly shows a BUSCO completeness of 98.39%. We annotated 659.55 Mb repetitive sequences and 26,036 protein-coding genes, 99.47% of which are functionally annotated. Comparative phylogenomic analysis clarifies the phylogenetic position of Saurogobio within Gobioninae. This high-quality genome provides a critical genetic basis for exploring cyprinid phylogeny, benthic adaptive evolution, genetic improvement, and conservation efforts of S. dabryi.

Saurogobio dabryi↗

PubMed2026

Genomic and phenotypic characterization of mcr-9 carrying Enterobacter oligotrophicus recovered from bovine mastitis.

Bovine mastitis (BM) is a leading cause of economic loss in the dairy industry, driven by decreased milk yields, involuntary culling, and substantial veterinary costs. A single Gram-negative isolate recovered from BM was characterized in this study. For this, antimicrobial susceptibility testing (AST) was performed using the Neg-Urine-Combo 98 panel. Whole-genome sequencing (WGS) was employed to identify antimicrobial resistance genes (ARG), virulence factors (VF) genes, plasmid replicons and prophage sequences. Comparative genomic analysis was performed through phylogenetic analysis. The BM isolate was phenotypically identified as Citrobacter rodentium, however, WGS analysis reclassified the isolate as Enterobacter oligotrophicus. Phenotypic AST revealed a resistance profile of 12%, with the isolate exhibiting resistance to β-lactams antimicrobials, specifically amoxicillin-clavulanate, ampicillin and the cephalosporin-cefoxitin. Conversely, susceptibility was demonstrated for the remaining tested antimicrobials. Genomic profiling identified 31 ARG, 10 VF genes and 6 prophage sequences within the E. oligotrophicus genome. This study provides the first evidence of E. oligotrophicus as a causative agent of BM, expanding the known spectrum of pathogens affecting the dairy industry by delivering the second complete genome of E. oligotrophicus available globally. The identification of 31 ARG, 10 VF, and 6 prophage sequences underscore the potential pathogenic risk and environmental resilience of this isolate. These findings highlight the critical role of WGS-based surveillance in identifying non-conventional mastitis pathogens and underscore the need for targeted mitigation strategies in veterinary medicine.

Animals↗

PubMed2026

CAKR: commutative algebra k-mer representations for genomics.

Despite the availability of various sequence analysis models, comparative genomic analysis remains a challenge in genomics, genetics, and phylogenetics. Commutative algebra, a fundamental tool in algebraic geometry and number theory, has rarely been used in data and biological sciences. In this study, we introduce commutative algebra k-mer representations as a nonlinear algebraic framework for analyzing genomic sequences. This representation bridges commutative algebra, algebraic topology, combinatorics, and machine learning to establish a mathematical framework for comparative genomic analysis. We evaluate its effectiveness on three tasks including genetic variant classification, phylogenetic tree reconstruction, and viral classification, typically requiring alignment-based, alignment-free, and machine-learning approaches, respectively. In this work, we show that commutative algebra k-mer representations outperform five state-of-the-art sequence analysis methods across twelve primary datasets, with two additional supplementary fragment-placement benchmarks, especially in viral classification, and maintain relatively stable predictive accuracy as dataset size increases, underscoring scalability and robustness.

Genomics↗

PubMed2025

CAKL: Commutative algebra k-mer learning of genomics.

Despite the availability of various sequence analysis models, comparative genomic analysis remains a challenge in genomics, genetics, and phylogenetics. Commutative algebra, a fundamental tool in algebraic geometry and number theory, has rarely been used in data and biological sciences. In this study, we introduce commutative algebra k-mer learning (CAKL) as the first-ever nonlinear algebraic framework for analyzing genomic sequences. CAKL bridges between commutative algebra, algebraic topology, combinatorics, and machine learning to establish a new mathematical paradigm for comparative genomic analysis. We evaluate its effectiveness on three tasks-genetic variant identification, phylogenetic tree analysis, and viral genome classification-typically requiring alignment-based, alignment-free, and machine-learning approaches, respectively. Across eleven datasets, CAKL outperforms five state-of-the-art sequence analysis methods, particularly in viral classification, and maintains stable predictive accuracy as dataset size increases, underscoring its scalability and robustness. This work ushers in a new era in commutative algebraic data analysis and learning.

Journal Article↗

PubMed2025

Genomic analysis of the liverpool epidemic strain of pseudomonas aeruginosa infecting persons with cystic fibrosis reveals likely Canadian origins.

INTRODUCTION: The Liverpool Epidemic Strain (LES) of Pseudomonas aeruginosa is one of several known strains to be transmissible between persons with cystic fibrosis (CF) (pwCF) and the only known strain to have infected large proportions of CF populations on two continents. Despite its prevalence, efforts to understand its spread have proven elusive. METHODS: We leveraged a prospective collection of P. aeruginosa isolates from pwCF attending the Southern Alberta Adult CF clinic from 1986 to 2020 to identify all individuals with LES infection. LES isolates collected every 1-2 years from each pwCF were sequenced and compared with 171 published LES genomes by phylogenetic analysis. RESULTS: Of 395 pwCF screened, ten pwCF infected with the LES were identified, from whom 46 LES isolates were sequenced. The earliest LES isolate was recovered in 1986, ∼2 years earlier than the previously oldest published LES isolate recovered in the UK. Phylogenetic analysis identified a diverse set of isolates at the root of the LES phylogeny that formed four clades, one of which gave rise to a "classic LES" clade. Canadian isolates formed a paraphyletic group that included the root of this clade and out of which the UK LES clade emerged. We estimated the date of the most recent common ancestor (MRCA) of the UK LES clade as 1977. CONCLUSIONS: Our study provides genomic evidence in support of a silent epidemic of LES infection occurring in the late 1970s among pwCF first originating in Canada and being spread to the UK, where transmission markedly accelerated.

Humans↗

PubMed2026

Emergence of cefiderocol resistance in carbapenem-resistant Escherichia coli ST167 prior to clinical use: A multifactored resistance landscape.

OBJECTIVES: Cefiderocol is a novel siderophore cephalosporin with potent activity against multidrug-resistant Gram-negative bacteria. Here, we reported the prevalence and mechanisms of cefiderocol resistance in carbapenem-resistant Escherichia coli (CREC) in China before its clinical use. METHODS: A total of 443 non-duplicate CREC isolates collected from 67 hospitals in China (2013-2021) underwent antimicrobial susceptibility testing according to CLSI guidelines. Whole-genome sequencing, transcriptomic analysis, siderophore quantification, and targeted genetic manipulation were performed to investigate the underlying resistance mechanisms. RESULTS: Among the 443 CREC isolates, 102 (23.0%) were resistant to cefiderocol, and 34 (7.6%) showed intermediate susceptibility. Multivariable logistic regression identified ST167 lineage (OR, 3.05; 95% CI, 1.12-8.29; P = 0.028), blaNDM-5 carriage (OR, 9.04; 95% CI, 2.96-27.57; P < 0.001), and cirA truncation (OR, 49.56; 95% CI, 20.33-120.79; P < 0.001) as independent factors associated with cefiderocol resistance. Among ST167 isolates, cefiderocol-resistant isolates showed increased yersiniabactin carriage and siderophore production but comparable TonB-dependent transporter expression profiles. Phylogenetic analysis revealed that cefiderocol-resistant ST167 isolates clustered into a distinct subclade enriched with resistance-associated determinants, including a recurrent FhuA P50S substitution detected in 59/64 (92.2%) resistant isolates. Functional assays showed that the P50S substitution increased cefiderocol minimum inhibitory concentration (0.032-0.125 µg/mL), particularly in an NDM-5-producing background (0.032-0.5 µg/mL). CONCLUSIONS: Cefiderocol resistance is highly prevalent among high-risk ST167 CREC isolates before the clinical introduction of cefiderocol in China, highlighting the need for continued surveillance of this epidemic lineage. Cefiderocol resistance is mediated by multiple resistance determinants, and we identify the recurrent FhuA P50S substitution as a novel contributor to reduced cefiderocol susceptibility.

Antimicrobial resistance↗

PubMed2024

Expansion of Oropouche virus in non-endemic Brazilian regions: analysis of genomic characterisation and ecological drivers.

BACKGROUND: Oropouche virus (OROV) is an arbovirus endemic in the Amazon region that closely resembles other arboviruses in terms of human disease, leading to potential misdiagnoses. The virus ecology has mostly restricted its occurrence to the Amazon biome; however, after a large 2023-24 OROV epidemic in the Brazilian Amazon region, outbreaks are being reported across Brazil and in other countries in Latin America. Here, we investigate the OROV spread outside Amazonia. METHODS: In this genomic and epidemiological study, OROV cases from January, 2023, to July, 2024, provided by the General Coordination of Public Health Laboratories of Brazil on Aug 1, 2024, were compared by geographical location (Amazon vs non-Amazon) and municipal population size, and a linear mixed model was employed to assess the relationship between agricultural area size and cases. OROV-positive samples from central laboratories of five non-Amazonian Brazilian states were sequenced using an amplicon-based approach. Bayesian phylogeographical analysis was performed with near full-length viral genomes, incorporating individual travel histories when relevant. The estimated dates of viral introductions in each sampled location were then contextualised with public epidemiological data. FINDINGS: Epidemic data show that outside the Amazon region, OROV cases frequency was 3·9-times higher in small municipalities than in large municipalities. The planted areas of some agricultural products, such as banana plantations, were positively correlated (r=0·39, p<0·0001) with OROV cases. The linear mixed model revealed that, besides banana, cassava also has larger (p<0·05) planted areas in municipalities with OROV cases when compared with those with no cases. The phylogenetic analysis of 32 new OROV genomes reconstructed multiple exportation events of the newly identified reassortant lineage from the Amazon to other Brazilian regions between January and March, 2024. At least three of the previously described OROV phylogenetic clades circulating in the Amazon were the source of viral introductions. Molecular clock analysis estimated that viral introductions happened from 50 days to 100 days before detecting the outbreaks in each state. INTERPRETATION: Our results confirm that the novel OROV reassortant lineage spread from the Amazon to other regions in early 2024, successfully establishing local transmission. The fact that outbreaks were observed in small municipalities, instead of large urban centres, suggests that local ecological conditions that are ideal for OROV vector occurrence, such as the banana plantation environment, might be important factors driving its spread in Brazil. FUNDING: DECIT, CNPq, FAPEAM, and Inova-Fiocruz. TRANSLATION: For the Portuguese translation of the abstract see Supplementary Materials section.

Brazil↗