Plastid genome evolution and phylogenomics with broad taxon sampling: insights into intrafamilial classification of Hamamelidaceae

Abstract

Hamamelidaceae, within the order Saxifragales, comprises 27 genera and approximately 120 species. The family has a pantropical and temperate distribution across the Americas, Asia, Africa, and Australia. Previous molecular investigations, constrained by limited taxon sampling and inadequate genetic markers, supported a five-subfamily classification system. However, these studies predominantly focused on Asian taxa, resulting in poor resolution of the evolutionary relationships among American, African, and Australian genera. To address these sampling gaps, we employed near-complete generic sampling (26 of 27 genera) to investigate plastome architecture, structural variation, and phylogenetic relationships. We newly sequenced and assembled 15 plastid genomes representing geographically and taxonomically underrepresented genera and analyzed them alongside 59 publicly available plastomes retrieved from GenBank. Plastid genomes exhibited conserved quadripartite architecture with sizes ranging from 158, 076 bp to 160, 814 bp, minimal structural variation, consistent GC content (37.7-38.2%), and identical gene order. Inverted repeat (IR) regions had limited size variation (26, 211-26, 429 bp). Simple sequence repeat (SSR) distribution (2, 219 loci) showed no clear correlation with the genus-level phylogenetic relationships. We identified ten hypervariable regions, including coding sequences ( accD, ycf1, clpP, ndhF , and rpl22 ) and intergenic spacers ( rpl33-rps18 , the trnG-UCC intron, trnH-GUG-psbA, accD-psaI , and petA-psbJ ), as promising candidate regions for future applications in species delimitation and phylogenetic studies. Phylogenetic analyses revealed largely congruent topologies across datasets and methods, providing improved resolution and strong support for most subfamilial and tribal relationships compared with previous studies. This study highlights the utility of plastid genome data for resolving deep-level phylogenetic relationships within Hamamelidaceae. The genome architecture reflects the high conservation of plastid genomes, while the identified mutation hotspots represent potential resources for future taxonomic and phylogenetic studies. Our results support the existing subfamily classification while improving geographical coverage and generic representation, providing a robust framework for future taxonomic and evolutionary studies of this globally distributed and taxonomically complex family.