Abstract
The genus Stachys, a large and taxonomically intricate group within the Lamiaceae family, includes numerous species valued in traditional medicine, yet their evolutionary relationships are often obscured by morphological complexity. To resolve these persistent difficulties, we sequenced and assembled the complete chloroplast genomes of two medicinally important Iranian species, Stachys persica and Stachys germanica, and compared them against seven congeneric genomes. Despite an overall picture of strong structural conservation across the genus, we discovered that evolutionary variation concentrates in specific genomic hotspots. Several non-coding intergenic spacers, most notably trnL-trnF, along with the matK and ycf1 genes emerged as exceptionally variable regions, making them promising candidates for species-level DNA barcoding. In a more surprising finding, we detected clear signatures of positive selection acting on the photosynthetic gene petB, suggesting that even core energy pathways have experienced adaptive refinement within this lineage. Phylogenetically, whole-plastome data placed S. germanica as sister to S. byzantina, with S. persica as their closest relative, a relationship that individual barcode loci failed to recover reliably. Taken together, this study provides a robust evolutionary framework for Stachys, identifies practical molecular tools for authenticating medicinal materials, and demonstrates the power of whole-plastome sequencing to untangle taxonomically intricate plant groups.