Search PubMedSearch

SEARCH · Search PubMed

Results for “phylogeny”

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.

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

88 recordsLinked to original sources

Uce-based phylogeny and classification of Megachilini.

The generic-level classification of the bee tribe Megachilini (Megachilidae) has remained controversial due to poor phylogenetic resolution at the base of the group, particularly among the brood parasitic genera and the numerous dauber ("Chalicodoma s. l.") lineages. We present a phylogenomic analysis of Megachilini based on ultraconserved elements (UCEs), sampling 52 ingroup taxa with emphasis on the dauber lineages. We also present a combined UCE + six-gene analysis to improve taxon coverage, resulting in a dataset with 127 ingroup taxa. Maximum likelihood, coalescent, and Bayesian analyses of multiple UCE matrices recover largely congruent topologies with substantially improved support relative to previous studies. Our results strongly support the monophyly of Megachilini, the early divergence of Noteriades and Gronoceras, and a single origin of brood parasitism. All remaining non-parasitic Megachilini form a moderately supported clade sister to the brood parasitic lineage. The leafcutter bees are monophyletic and nested within dauber lineages. Several major dauber clades are consistently recovered, including an exclusively Australian clade corresponding to the Hackeriapis group of subgenera, while several recognized subgenera are paraphyletic. The lineage known as Morphella, previously placed in synonymy with the subgenus Callomegachile, was not closely related to that subgenus and is here treated as a valid subgenus. Divergence-time analyses place the crown age of Megachilini in the late Eocene to early Oligocene, with major extant lineages diversifying during the Miocene. Limited morphological diagnosability of several clades indicates that splitting non-parasitic lineages into numerous genera would result in an impractical classification that would widen the gap between taxonomists and non-specialists and exacerbate the taxonomic impediment in bees. We therefore advocate retaining a single genus Megachile for non-parasitic Megachilini (excluding Noteriades and Gronoceras), as the classification best supported by phylogenomic evidence and most robust to future taxon sampling.

Animals

Viral surveillance beyond detection: JMTV and the need for ensemble approaches in emerging virus discovery.

The recent report by T. Murillo, L. E. Enrique Chaves-González, S. Temmam, S. Bermúdez, et al. (Microbiol Spectr 14:e04078-25, 2026, https://doi.org/10.1128/spectrum.04078-25) expands the known geographic and ecological range of Jingmen tick virus (JMTV) by detecting the virus in Amblyomma mixtum ticks collected from horses in Costa Rica. This is an important finding because A. mixtum can feed on wildlife, domestic animals, and humans, creating a possible interface for virus movement across various hosts. The study also places the Costa Rican virus in a wider phylogenetic context, linking it to JMTV diversity reported from other regions. However, the detection of viral RNA in ticks should not be interpreted as proof of local disease, human infection, or active transmission, especially in the absence of supporting results. Instead, it reflects an important signal for careful viral surveillance. Here, I discuss how JMTV illustrates the need for ensemble approaches that combine field sampling, phylogeny, segment-level genome analysis, serology, experimental validation, and data-driven virus discovery tools.

emerging viruses

First insights into the role of evolutionary history in shaping venom composition of Vipera ammodytes.

Understanding intraspecific venom variation requires distinguishing the contributions of neutral population history from natural selection. This study aims to determine whether venom variation in Vipera ammodytes species complex is structured across eight phylogenetic groups. Despite a complex evolutionary history, venom composition did not differ among phylogenetic groups within the analytical framework used, suggesting that shared ancestry alone does not explain venom variation. Whether local adaptation to environmental conditions explains the observed variation remains an open question for future studies.

Animals

Pervasive hybridization and introgression in Diervilleae (Caprifoliaceae).

Diervilleae (Caprifoliaceae) is a horticulturally important lineage with striking floral diversity and a long history of interspecific crossing, suggesting reticulate evolution. We integrated nuclear SNPs and whole plastome data to reconstruct a phylogenomic backbone for the tribe and to identify hybrids, cultivated accessions, and introgression among lineages. Nuclear and plastid phylogenies consistently recover Weigela and Diervilla as reciprocally monophyletic and resolve four major lineages within Weigela, providing a reproducible framework for revising sectional limits and species boundaries. Cultivated accessions form a well supported clade sister to W. florida and show predominantly W. florida ancestry while retaining contributions from multiple wild lineages, consistent with recurrent crossing, backcrossing, and selection. Analyses of wild populations reveal recurrent hybrids and enable plausible parental combinations to be inferred. Tests across the genome further indicate strong evidence for historical introgression across Diervilleae, with the strongest signals involving W. middendorffiana, W. maximowiczii, and Diervilla. Fossil evidence, divergence time estimation, and paleodistribution modelling together suggest range expansion during the Miocene and Pliocene followed by climate driven contraction, providing a spatiotemporal context for episodic contact, introgression, and the East Asia-North America disjunction.

Hybridization, Genetic

Integrating genomic distance analyses in the description of a new family, genus, and species of sponge-associated antipatharians (black corals).

Antipatharians (black corals) are among the least studied coral groups, with much of their diversity still undescribed. Here, we present an integrative morphological, phylogenomic and genomic distance study of deep-sea antipatharians sampled in high seas areas of the North Pacific Ocean and from New Zealand's Exclusive Economic Zone. These corals grow on hexactinellid sponges - a unique characteristic in the order Antipatharia. Using a dataset of ultra-conserved elements and exons, combined with morphological analyses, we reconstruct phylogenomic relationships and formally describe a new family (Eidikopathidae fam. nov.), a new genus (Eidikopathesgen. nov.), and two new species (E. korallispongiasp. nov., E. zealandkoralliasp. nov.). Morphologically, the new family is distinguished by a corallum consisting of a network of loose branches that fuse with the sponge skeletal framework. Phylogenomic analyses recovered consistent topologies with strong nodal support, corroborating the distinct evolutionary placement of this sponge-associated lineage. Pairwise genomic distances estimated using the Tamura-Nei model were concordant with patristic genomic distances, identifying Pteridopathidae as the genetically closest family to Eidikopathidae fam. nov., followed by Myriopathidae and Stylopathidae, which were recovered as sister families in the phylogeny. This pattern shows that genomic distance complements, rather than simply mirrors, tree topology by quantifying accumulated sequence divergence among lineages. Together, these results provide the first genomic distance framework for Antipatharia, offering a baseline for future systematic, evolutionary, and biodiversity studies on this fundamental shallow, mesophotic and deep-sea coral group.

Animals

Genome-wide SNP data support species boundaries in sympatric Polylepis Ruiz & Pav. (Rosaceae) species from Bolivia and Ecuador.

Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.

Bolivia

Algae-to-host horizontal gene transfer in Paramecium bursaria is associated with host adaptation during endosymbiosis.

Paramecium bursaria maintains a stable endosymbiosis with green algae, yet the evolutionary consequences of this association remain unclear. Here, we screened the host genome for algal-derived horizontally transferred genes (HTGs) using a lineage-aware workflow designed to detect horizontal gene transfer (HGT) between two defined lineages. We identified 16 candidate HTGs, including four putative newly transferred genes and 12 homologous transferred genes, most of which were functionally associated with redox homeostasis and metabolism. Five HTGs showed symbiosis-dependent expression. RNAi knockdown of GH32s and SATs reduced host proliferation, total cell area, and motility, while GH32s knockdown also reduced endosymbiont load. Duplication patterns suggest that most transfers may have occurred after the P. bursaria lineage diverged from the sampled Paramecium species but before its lineage-specific whole-genome duplication (WGD). The HTGs also showed host-associated shifts in GC content and gene length, while representative HTGs retained conserved domains and functional motifs. Together, our results support algae-to-host HGT in P. bursaria and suggest that some transferred genes may contribute to metabolic integration during endosymbiosis.

Gene Transfer, Horizontal

Longitudinal whole-genome analysis of bluetongue virus identifies conserved serotype-specific genomes and distinct genomic constellations within a Colorado sheep flock (2021-2023).

Bluetongue virus (BTV) is a segmented double-stranded RNA virus of ruminants transmitted by Culicoides spp. biting midges. Although the genome consists of ten segments, classification into serotypes is primarily based on genome segment 2. However, reassortment among genomic segments is a major driver of BTV evolution and diversity. This study used longitudinal whole-genome sequencing to characterize BTV genomes collected from 2021 to 2023 within a single sheep flock in Colorado, where multiple serotypes co-circulate. Whole-genome sequences were generated from fourteen blood samples representing four serotypes: BTV-6, -11, -13, and -17. Longitudinal sampling identified multiple BTV serotypes within individual sheep across consecutive years. Tanglegram analysis comparing segment phylogenies to the segment 2 tree demonstrated incongruent topologies across all genomic segments, suggestive of reassortment or the circulation of distinct genomic constellations. Nucleotide-level comparisons revealed high sequence homology among same-serotype samples from the same year, while the greatest genetic divergence was observed among BTV-17 genomes collected in different years. Additionally, all BTV-13 genomes contained a previously undescribed nonsynonymous substitution in segment 10 predicted to extend the encoded protein by three amino acids. Together, these findings demonstrate that highly conserved BTV genomes and distinct genomic constellations can be detected at the flock level across multiple years. This longitudinal whole-genome approach reveals the genetic complexity of endemic BTV populations, including novel variants and genomic patterns consistent with reassortment that are lost with conventional serotyped-based approaches, highlighting the need to integrate whole-genome characterization into endemic BTV monitoring programs.

Animals

Cytonuclear conflict and reticulate evolution in the Morelloid clade (Solanum, Solanaceae): Insights from genome skimming and network Phylogenomics.

The Morelloid clade (black nightshades) is one of the most strongly supported clades within the megadiverse Solanum genus. It comprises 76 globally distributed, non-spiny herbaceous and suffrutescent species. While often erroneously considered poisonous weeds, several species are economically important as orphan crops. The clade is closely related to tomato and potato but, due to a lack of focused breeding efforts, remains a putative reservoir of genetic diversity for crop improvement. Despite this potential, we lack fundamental knowledge on the evolution of the Morelloid clade. The group includes polyploid species with unknown parental origins-likely reflecting reticulate processes such as hybridization, introgression, and associated backcrossing events. Prior analyses have been unable to disentangle these processes, leaving the mechanisms underlying reticulate evolution in the Morelloid clade poorly understood. Here, we use genome skimming to produce a well-supported maximum likelihood plastid phylogeny from complete circularized plastomes and a coalescent-based species tree from combined Angiosperms353 and conserved ortholog set nuclear markers. Our dataset, composed of previously published data and deep genome skimming from herbarium samples, spans 26 Morelloid species. To investigate phylogenetic discordance, we used a nuclear phylogenetic network, multispecies coalescent simulations, a fused rooted nuclear chloroplast tree, and quantification of nuclear gene tree concordance. We show that incongruence between nuclear and plastid trees is pervasive and cannot be explained by incomplete lineage sorting alone. Instead, our results demonstrate that events consistent with repeated chloroplast capture have shaped the reticulate evolutionary history of the clade, especially among African polyploid and Pan-American diploid lineages.

Phylogeny

Molecular adaptation of caspase genes to salinity stress in the tropical sea cucumber Stichopus monotuberculatus: A comparative analysis across echinoderms.

Apoptosis is an essential physiological process that plays a critical role in development and tissue homeostasis. Caspases, as central regulators of apoptosis, are crucial in controlling inflammation and cell death. In this study, we investigated the caspase gene family in Stichopus monotuberculatus to explore their potential roles in salinity stress adaptation. Five caspase genes were identified from the genome of S. monotuberculatus, including Smcaspase3, Smcaspase6, Smcaspase8a, Smcaspase8b, and Smcaspase8c. Phylogenetic analysis revealed that these Smcaspase genes clustered into distinct caspase subfamilies and showed high conservation with homologs from other echinoderms and representative vertebrates. Conserved motif and gene structure analyses showed relatively similar structural patterns within each clade, whereas divergence was observed among different subfamilies. Promoter analysis identified numerous cis-acting elements related to gene regulation, immune response, and growth and development. Expression profiling under salinity stress showed that Smcaspase8a was significantly upregulated, particularly under prolonged stress, whereas the other genes exhibited limited transcriptional responses. Our findings highlight caspase function in salinity stress and provide the foundation of molecular salinity adaptation mechanisms in S. monotuberculatus.

Animals

Assembly and Characterization of the First Complete Mitochondrial Genome of Tussilago farfara L.: Insights into Biological Functions and Phylogenetic Relationships within the Asteraceae Family.

Tussilago farfara L., a member of the Asteraceae family, is an economically valuable species due to its edible and medicinal properties. To elucidate the structural characteristics, genetic mechanisms, and evolutionary pathways of the organelle genomes of T. farfara, we sequenced, assembled, and annotated its mitochondrial genome for the first time. The complete mitochondrial genome of T. farfara spans 306,024 bp and contains 33 mitochondrial protein-coding genes (PCGs), 3 rRNAs, and 22 tRNAs. Analysis of the nucleotide substitution rate and genetic diversity revealed that most mitochondrial genome genes may have undergone purifying selection, indicating a slow evolutionary rate and a relatively conserved genomic structure. We further identified 13 fragments of chloroplast-derived DNA integrated into the mitochondrial genome, evidencing intracellular gene transfer. Collinearity analysis showed that Arctium lappa shares the most extensive mitochondrial homologous sequences and the highest sequence similarity with T. farfara. Phylogenetic analysis based on the mitochondrial genome helped to clarify the evolutionary and taxonomic position of T. farfara within the Asteraceae family. The mitochondrial genome sequence of T. farfara provides a valuable genomic resource for species identification and for evolutionary studies within the Asteraceae family.

Genome, Mitochondrial

Evolutionary architecture and lineage-specific diversification of Forkhead box transcription factors in Perna viridis.

The Forkhead box (Fox) transcription factors are evolutionarily conserved regulators of development, cell cycle, and apoptosis across metazoans. This study provides the first comprehensive genome-wide analysis of the Fox gene family in the Asian green mussel (Perna viridis). We identified 28 Fox genes distributed across 10 chromosomes. Comparative analysis reveals the absence of the FoxI, FoxQ1, FoxR and FoxS subfamily, consistent with other bivalves and indicative of lineage-specific gene loss during molluscan evolution. Notably, gene duplications in the FoxAB, FoxD, FoxH, FoxN1-4, FoxQ2 and FoxQD subfamilies may reflect functional diversification associated with environmental adaptation. Exon-intron structural variability, including intron loss in several paralogues, suggests structural diversification and potential regulatory variation. Phylogenetic reconstruction confirmed the monophyly of core Fox classes while highlighting divergent expansion patterns in lophotrochozoans. Selection analyses showed strong purifying selection across duplicated Fox paralogs, supporting functional conservation after lineage-specific expansion. Gene Ontology enrichment linked Fox genes to stress response, apoptosis, and transcriptional regulation. By integrating phylogenetic, structural, and transcriptomic analyses, this study provides a genomic framework for understanding Fox gene organisation, evolution, and tissue-associated expression patterns in Perna viridis and establishes a comparative resource for future functional studies in bivalves.

Animals

LitCTL1: A novel C-type lectin involved in the mucosal and cellular immunity of the common periwinkle Littorinalittorea.

C-type lectins (CTLs) are vital pattern-recognition receptors (PRRs) that mediate innate immune responses in mollusks, yet their characterization in Caenogastropoda, the largest gastropod group, remains limited. This study characterizes LitCTL1, a novel secreted single-domain C-type lectin from the common periwinkle, Littorina littorea. The 199-amino acid polypeptide contains a conserved carbohydrate recognition domain with canonical QPD and WND motifs and is predicted to form a homodimer. Uniquely, LitCTL1 was localized in both circulating hemocytes and mucus-secreting epithelial cells of the foot, mantle, and hypobranchial gland - the first report of such dual localization for a molluscan lectin, linking systemic and mucosal defense. Expression analysis revealed that LitCTL1 is constitutively expressed in hemocytes. Functional assays with recombinant LitCTL1 demonstrated its role as a potent opsonin with hemagglutinating activity, significantly enhancing hemocyte spreading and the phagocytosis of zymosan. Genomic analysis reveals that LitCTL1 belongs to a rapidly diversifying, genus-specific expansion distinct from conserved perlucin-like lineages. These results identify LitCTL1 as a key effector molecule in both systemic and mucosal innate immunity, likely reflecting an evolutionary adaptation to the microbial challenges of the intertidal environment.

Animals

Genetic diversity and recombination of NA-PRRSV field strains in Vietnam: Implications for vaccine efficacy.

Porcine reproductive and respiratory syndrome (PRRS) causes severe reproductive losses in pregnant sows and piglets, resulting in substantial economic impact on the swine industry worldwide. However, due to the significant genetic diversity and rapid evolutionary changes of the pathogen, continuous surveillance and detailed genetic analysis of circulating strains are essential. The current study aimed to evaluate the genetic diversity of the hypervariable (HV) region of non-structural protein 2 (nsp2) among North American PRRSV strains isolated from swine farms in Vietnam. Phylogenetic analysis and multiple sequence alignment were conducted to determine subtype classification and assess genetic variability. A total of 48 field isolates were obtained, of which 12.5% belonged to classical NA-PRRSV, 16.6% to NADC30-like and 70.9% to HP-PRRSV, primarily distributed across sublineages 1.4, 5.1, 8.7 and 8.9. Amino acid comparisons found multiple insertions, deletions and substitutions at various positions within the hypervariable region of nsp2. The study revealed substantial genetic variation in the HV region of nsp2 among NA-PRRSV field strains, largely associated with recombination and immune escape. These findings highlight epidemiological risks to vaccine efficacy and underscore the need for continuous molecular surveillance to support effective PRRSV control in Vietnam.

PRRSV

Ramu stunt virus genome reveals previously unreported segments and nucleocapsid domain duplication in Mechlorovirus.

Ramu stunt virus (RmSV), a member of the genus Mechlorovirus within the family Phenuiviridae, was previously described as a six-segmented RNA virus infecting sugarcane. In this study, we re-examined type material and additional isolates using high-throughput sequencing and RT-PCR validation, revealing that RmSV possesses a nine-segmented genome, making it the largest reported in the Phenuiviridae. This expanded architecture includes duplicated RNA segments (RNA 2a and RNA 2b) encoding nucleocapsid-like proteins and two novel segments (RNA 7 and RNA 8). Comparative analysis showed that RNA 2a and 2b share about 84% amino acid identity, while RNA 5 encodes a third nucleocapsid homolog, indicating unprecedented domain redundancy. Structural modeling confirmed that all three nucleocapsid proteins maintain a conserved fold despite low sequence identity, with electrostatic mapping suggesting differential RNA-binding potential. Additionally, RNA 6 encodes a hypothetical protein structurally similar to the rice stripe virus disease-specific S-protein, implicating a role in symptom development. Transcript abundance analysis revealed RNA 6 as the most highly expressed segment across isolates. These findings revise the genomic composition of RmSV, highlight mechanisms of genome plasticity and adaptive evolution in plant-infecting bunyaviruses, and underscore practical implications for diagnostic assay design, resistance breeding, and biosecurity surveillance.

Genome, Viral

Genome-Wide Characterization of β-Glucosidase (TaBGLU) Genes in Bread Wheat and Their Expression Under Drought, Cold, and Combined Stress.

Glycoside hydrolase 1 (GH1) β-glucosidases were known to activate hormone conjugates and defense metabolites, yet their genomic organization and stress-response dynamics in wheat remained incompletely defined. We therefore performed an integrated characterization of TaBGLUs spanning phylogeny, gene structure and conserved motifs, subcellular localization, promoter cis-elements, Gene Ontology enrichment, protein-protein interaction networks, and targeted expression profiling. Wheat TaBGLUs partitioned into well-supported clades that shared canonical GH1 catalytic residues and a largely conserved motif scaffold. Subcellular localization predictions indicated predominant nuclear and chloroplast targeting, with a smaller cohort directed to secretory or endomembrane compartments. Promoters were enriched for light-responsive, hormone-related (ABA, JA/SA, auxin, GA) and stress-associated (MYB/WRKY, heat, low temperature) cis-elements, and functional annotations were consistent with roles in carbohydrate and cell-wall metabolism, hormone homeostasis, and defense. Network analysis revealed a densely connected TaBGLU submodule embedded within broader carbohydrate and defense interaction networks, suggesting coordinated or cooperative functions. Expression profiling under cold, drought, and combined drought and cold demonstrated broad stress inducibility, with early activation detected by 6 h, cold-responsive maxima typically at 12 h, drought-responsive peaks predominating at 24 h, and combined stress eliciting both earlier and more sustained expression maxima between 12-24 h. Representative strongly responsive genes included TaBGLU20, TaBGLU44, TaBGLU6, and TaBGLU23, which showed pronounced late induction under combined stress, TaBGLU30, which exhibited an earlier combined-stress peak, and TaBGLU12, which displayed a marked late drought-specific response. Taken together, this integrated genomic, regulatory, and expression atlas refined the wheat BGLU repertoire relative to previous gene model inventories, highlighted candidate TaBGLUs with central network positions and strong stress inducibility, and provided concrete entry points for functional validation and breeding for improved stress resilience.

Triticum

Molecular characterization of 16 MAPK genes in silver carp (Hypophthalmichthys molitrix) and the differences of their mRNA expression between Qiandao Lake and Taihu Lake.

Mitogen-activated protein kinase (MAPK), a serine-threonine protein kinase, is involved in a variety of stress-induced responses and also plays an important regulatory role in cell metabolism. In the study the open reading frames (ORFs) of 16 MAPK genes in silver carp (Hypophthalmichthys molitrix) were obtained and verified, with the evaluations of their taxonomy, structures, conserved motifs, and evolutionary linkages. And the expression patterns of these genes in the silver carp from Qiandao Lake and Taihu Lake were explored for better understanding the response of MAPK genes to different water environment. MAPK genes of silver carp were divided into three subfamilies, including extracellular signal-regulated kinase (ERK) subfamily, p38 subfamily and C-Jun N-terminal kinase (JNK) subfamily. All these genes possessed similar structures and conserved motifs of MAPK family. Realtime qPCR revealed that the expression patterns of 10 MAPK genes (ScMAPK1, ScMAPK3, ScMAPK4, ScMAPK7, ScMAPK15, ScMAPK8a, ScMAPK8b, ScMAPK9, ScMAPK10 and ScMAPK11) in head kidney, spleen and gill of silver carp in Taihu Lake and Qiandao Lake were different. These findings provide a basis for further research on the function of MAPK in silver carp.

Animals

Genomic and Phenotypic Characterization of Two Novel Enterobacter Phages With EDTA-Enhanced Antibiofilm Activity.

Multidrug-resistant members of the Enterobacter cloacae complex (ECC) are increasingly linked to difficult-to-treat infections and biofilm-mediated antimicrobial tolerance. Here, two lytic phages, vB_EhoIP_HHH and vB_EluM_RZH, displaying podovirus-like and myovirus-like morphology, respectively, were isolated from the River Chelt. HHH has a 39,582 bp genome (51.2% GC, 63 ORFs), while RZH has a 174,197 bp genome (39.4% GC, 314 ORFs), with neither genome carrying antimicrobial resistance, virulence or lysogeny-associated genes. VIRIDIC and VICTOR analyses placed HHH within Kayfunavirus and RZH within Karamvirus, supporting their classification as distinct species. Both phages demonstrated rapid adsorption, short latent periods and stability across physiological pH and temperature ranges. A phage cocktail targeting MDR ECC strain was evaluated with EDTA against established biofilms. Crystal violet assays showed the greatest biomass reduction at MOI 10 with 0.5-0.75 mM EDTA. Bliss independence analysis revealed localized synergy within this window but significant overall antagonism at higher EDTA concentrations. CFU enumeration confirmed greater activity against 24 h than 48 h biofilms. The optimized combination also reduced recoverable bacteria in a fibroblast infection model while maintaining low LDH release. These findings identify two novel lytic Enterobacter phages and support a narrow EDTA concentration window for enhanced phage-mediated antibiofilm activity.

Biofilms