3BIO: Journal of Biological Science. Technology and Management 8. : 28-43 http://journals. id/index. php/3bio e-ISSN: 2655-8777 In Silico Phylogenetic Analysis of Lamiaceae Based on ITS, matK, and rbcL DNA Barcodes Chayra Endlessa1. Topik Hidayat1*. Siti Sriyati1 Biology Study Program. Faculty of Mathematics and Science Education. Universitas Pendidikan Indonesia *) Corresponding author. e-mail: topikhidayat@upi. Received: 2024-10-22 Accepted for publication: 2025-10-27 Abstract Lamiaceae, widely used as herbal medicine, is increasingly vulnerable to adulteration driven by market demand, compromising product safety and efficacy. Prevention is challenging due to morphological similarities. DNA-based phylogenetics offer an alternative for accurate species authentication. However. Lamiaceae phylogenetics remain complicated by inconsistencies between morphological and DNA data. This study reconstructed Lamiaceae phylogeny using partial ITS, matK, and rbcL barcodes to evaluate their potential application in species authentication and adulteration prevention. Sequences for 52 species across 11 genera (Spathodea campanulata: outgrou. were obtained from NCBI GenBank, aligned, and trimmed. Four maximum parsimony (MP) trees were constructed in MEGA 11 . hree single-barcode, one concatenate. The concatenated dataset was also analyzed by maximum likelihood (ML). Tree robustness was evaluated with bootstrapping, consistency index (CI), and retention index (RI). matK had the longest mean sequence . 6 b. , rbcL the highest homology . 5%), and ITS the most parsimony-informative sites . 3%). MP trees exhibited moderate homoplasy . ean CI = 0. but strong synapomorphic signal . ean RI = 0. Individual barcodes produced similar genus groupings, yet misplaced several species. Concatenation corrected these positions across MP and ML trees, resolving six robust monophyletic clades . ootstrap >70%), broadly consistent with earlier phylogenies: Callicarpa. Scutellaria. Clerodendrum. Lamium, and Stachys. Salvia. Thymus. Origanum, and Mentha. Orthosiphon and Ocimum. Topological discrepancies with prior studies likely reflect differences in barcode choice and taxon sampling. Concatenated barcodes improved phylogenetic resolution in Lamiaceae, producing clades that identify potential adulterants and guide DNA marker development for species authentication and adulterant detection. Keywords: DNA barcodes. Lamiaceae. molecular phylogenetics. species authentication Introduction Lamiaceae, comprising approximately 7000 species across 236 genera . , is the largest family in the order Lamiales . Owing to the presence of therapeutic essential oils, 13% of its species are considered medicinal plants . The genus Salvia is among the oldest medicinal plants with antibacterial, antioxidant, and antiinflammatory properties . Additionally, notable species like Lavandula angustifolia and Origanum vulgare are extensively used in cosmetics and culinary applications . Increasing research on the benefits of herbal medicine over conventional treatments has driven global demand . , resulting in the adulteration of Lamiaceae species. For instance, peppermint oil (Mentha piperit. is often adulterated with menthol oil from Mentha arvensis . , while Origanum vulgare has been substituted with Origanum majorana . In some cases, the high market demand has driven some species to the brink of extinction due to exploitation . The unique biochemical profiles of Lamiaceae . suggest that adulteration compromises both the effectiveness and safety of herbal medicines. However, prevention efforts remain challenging due to morphological similarities among species, particularly in powdered form . For example, products containing Scutellaria may pose health risks if adulterated with the hepatotoxic genus Teucrium . Given the limitations of morphological identification, a DNA-based phylogenetic approach offers a reliable alternative for more accurate species authentication, as DNA sequences are less susceptible DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 to environmental factors. Using DNA barcodesAishort, standardized DNA sequencesAiin molecular phylogenetics has proven effective in estimating evolutionary relationships . However. Lamiaceae phylogenetics complex and unresolved . A study on Thymus struggled to delineate species boundaries due to high morphological similarity . Morphological convergence has also occurred across different lineages, as seen in staminal characters: the four stamens of the tribe Mentheae . ubfamily Nepetoidea. independently reduced to two in the subtribes Salviinae and Menthinae . Hybridization and polyploidy have also led to genomic alterations, creating complex phylogenetic networks and blurring taxonomic boundaries, as observed in Mentha. Stachys, and Ocimum . A study on Salvia further highlighted discrepancies between morphologyand DNA-based phylogenies, with the former supporting its monophyly and the latter suggesting paraphyly . Many phylogenetic studies on Lamiaceae have relied on a single genomic source, either nuclear or chloroplast DNA, limiting phylogenetic resolution. For example, researchers have studied Scutellaria using only chloroplast sequences . , and the most comprehensive Lamiaceae phylogeny, covering 78% of genera, also relied on chloroplast markers . This underscores the need to combine DNA barcodes to improve phylogenetic accuracy . Standard plant barcodes include the chloroplast regions maturase-K . atK) and ribulose1,5-bisphosphate carboxylase/oxygenase . bcL), with the nuclear internal transcribed spacer (ITS) recommended as a supplementary marker . , . To address this need, the present in silico study incorporated Lamiaceae species beyond those found in Indonesia by retrieving sequences from GenBank. This study aimed to obtain phylogenetic insights based on partial ITS, matK, and rbcL barcodes to support species authentication and prevent adulteration in Lamiaceae herbal products. Methodology Collection of Lamiaceae ITS, matK, and rbcL Sequences Partial DNA sequences from 52 Lamiaceae species across 11 genera were obtained, with Spathodea campanulata (Bignoniacea. as the outgroup and sister taxon . (Table The chosen species were relevant to previous studies addressing phylogenetic challenges and adulteration in Lamiaceae. Partial sequences were retrieved from the National Center for Biotechnology Information (NCBI) GenBank . ttps://w. gov/genbank/) using species names and barcode regions. The ITS dataset included partial ITS1 and ITS2 regions and the complete 5. ITS was included to complement plastid markers for its strong discriminatory ability at the species level and Endlessa et al. its recommendation as a core plant barcode . The fastevolving matK and highly conserved rbcL markers were selected as the standard barcode pair . to resolve recent and older phylogenetic relationships. All sequences were saved in FASTA format to create the DNA database. Bioinformatics Analysis DNA sequences were aligned separately for each barcode using the ClustalW . algorithm implemented in Molecular Evolutionary Genetics Analysis (MEGA) v11. under default parameters. Sequences were trimmed at the 5Ao and 3Ao ends so that all taxa within each marker had equal length, using the shortest sequence in the dataset as a reference. This ensured length uniformity within barcodes while maintaining length differences among ITS, matK, and rbcL. Phylogenetic analyses were conducted using the maximum parsimony (MP) model to evaluate the relative performance of ITS, matK, and rbcL, both individually and in combination. Four MP trees were constructed in MEGA under default settings: one from each barcode and the concatenated dataset. MP was selected as a computationally efficient method that infers topology based on the parsimony principle, or the minimum number of evolutionary changes . This model is appropriate for relatively small datasets of closely related sequences . , as in this study. To address the long-branch attraction associated with MP and increase phylogenetic resolution, the concatenated dataset was also used to construct a maximum likelihood (ML) tree in MEGA using the TamuraNei substitution model . Node support for both models was evaluated with 1000 bootstrap replicates, with higher percentages indicating stronger branch support . The robustness of the MP trees was further assessed using the consistency index (CI) and retention index (RI), which were calculated by MEGA after the MP analysis. These indices quantify phylogenetic signals by indicating how closely character similarities among taxa reflect their evolutionary relationships, with the signal reduced by homoplasy . CI values approaching 1 indicate low homoplasy . , whereas RI values close to 1 reflect a higher proportion of synapomorphies . DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Table 1. Partial Lamiaceae ITS, matK, and rbcL sequences from the NCBI GenBank. No. Genus Accession Number Species ITS Orthosiphon aristatus FJ593403 LC456391 MW789616 Orthosiphon stamineus AY506663 KM658969 MH069809 Thymus serpyllum KR150171 MF350183 MK105914 Thymus caespititius GU381457 HM850802 HM850398 Thymus vulgaris AY506646 OP243225 MN972464 Thymus quinquecostatus EU556524 LC618903 LC618880 Thymus mongolicus MH808603 MN433407 MN185199 Ocimum basilicum MW150025 MF694868 ON755091 Ocimum tenuiflorum MW150027 MF468149 JN114828 Ocimum gratissimum MW150026 MH552359 MW150006 Mentha piperita KY072944 KX783716 JQ230988 Mentha spicata GU381394 KC571807 KU499887 Mentha arvensis KY072946 MG224998 HQ590183 Mentha aquatica KR611529 KP172053 KC584892 Mentha suaveolens GU381395 KP172057 MG223550 Mentha longifolia KR611531 HQ902745 ON755102 Mentha canadensis KY072951 MT929800 JN407303 Salvia splendens MF622186 KX783777 ON755108 Salvia przewalskii MH808595 MN433404 JQ934026 Salvia officinalis KJ584196 HE967482 ON755112 Salvia miltiorrhiza MT039859 FJ513168 JQ934009 Salvia rosmarinus OQ165223 MH552339 MT931624 Salvia fruticosa KJ584194 HQ902726 HM590078 Salvia plebeia KU563788 MH660151 JQ934021 Clerodendrum cyrtophyllum KP092826 KJ888428 KJ939237 Clerodendrum japonicum KP092847 MK551817 GQ436521 Clerodendrum bungei EU591963 MH659049 JQ618463 Clerodendrum colebrookianum KX079329 MK551754 MK241954 Orthosiphon Thymus Ocimum Mentha Salvia Clerodendrum DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 No. Genus Endlessa et al. Accession Number Species ITS Callicarpa dichotoma KP092811 LC680459 LC694383 Callicarpa americana ON820115 MF350069 KY626890 Callicarpa macrophylla KP092818 OP032135 KF443315 Callicarpa kochiana KP092816 OP032127 KJ688019 Callicarpa giraldii FJ593347 OP032121 MH657300 Lamium album JX893229 MN311840 FJ395588 Lamium amplexicaule MN718246 MN433402 OL434812 Lamium galeobdolon KF529538 ON286905 JN891020 Origanum majorana JX162957 KX783725 JQ230991 Origanum vulgare AY506647 MF694869 ON755119 Origanum onites JX163054 HQ902752 HQ902807 Origanum dictamnus EU252137 FR719089 FR720564 Stachys sylvatica KF529644 JN895511 ON755118 Stachys palustris KF529624 JN894812 HE574636 Stachys floridana KF529590 OL434945 HQ644074 Stachys recta KF529631 KJ204541 KJ746271 Stachys arvensis KF529568 HM850806 MG224452 Stachys cretica KF529583 HQ902708 HQ902776 Scutellaria baicalensis MH711530 MH660079 KT280158 Scutellaria lateriflora MK356052 MG225186 HQ590266 Scutellaria indica MH808599 FJ513171 MN167869 Scutellaria barbata MF193539 FJ513170 FJ513144 Scutellaria viscidula MF193526 HQ676587 HQ676583 Scutellaria rehderiana JX893232 HQ676588 FJ513147 Spathodea campanulata MF616608 MF476853 MT933895 Callicarpa Lamium Origanum Stachys Scutellaria Outgroup: Spathodea (Bignoniacea. DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Results and Discussion Characteristics of Lamiaceae Nuclear and Chloroplast DNA Barcodes Table 2 summarizes the characteristics of the partial DNA barcodes used for reconstructing Lamiaceae phylogeny across a global dataset. Among the barcodes, matK had the longest average sequence length . 6 base pair. , whereas rbcL showed the highest homology . 5%), indicating a slower evolutionary rate. Homology refers to traits in organisms that arise from common ancestry . and are often conserved due to consistent inheritance, which explains the reduced variability observed in rbcL. Halmschlag et al. also found rbcL to be less variable than ITS and matK in their phylogeny of 89 Lamiaceae species from converted land in Sumatra. Indonesia. The conserved nature of rbcL likely reflects its essential role in encoding Rubisco, a key enzyme in photosynthesis and plant adaptation . The highest proportion of parsimony-informative (P. sites was observed in ITS . 3%), while rbcL showed the lowest . 6%). Comparable values were found by Halmschlag et al. , with 34% In ITS and 10% in rbcL. This consistent trend may be explained by differences in substitution rates, as the nuclear genome evolves approximately ten times faster than the chloroplast genome, where homologous recombination maintains genomic integrity . Pi sites represent inherited mutations that infer shared ancestry and close relationships among organisms . In multiple alignments, they highlight unique sequences within specific taxa, thereby increasing the resolution of phylogenetic The low number of Pi sites in rbcL further supports its conserved nature within Lamiaceae. Overall, the dataset exhibited higher homology than Pi sites, indicating low levels of nucleotide substitution in Lamiaceae. Endlessa et al. Lamiaceae Phylogenetics Based on Partial Nuclear and Chloroplast DNA Barcodes Four maximum parsimony (MP) phylogenetic trees were constructed from ITS (Figure . , matK (Figure . , rbcL (Figure . , and the concatenated dataset (Figure . maximum likelihood (ML) tree was also inferred from the combined dataset to improve phylogenetic resolution and mitigate long-branch attraction associated with MP (Figure Bootstrap values are presented on the branches of each tree, while the consistency and retention indices (CI and RI) are summarized in Table 3. MP tree topologies are assumed to reflect the principle of parsimony, which favors the minimum evolutionary changes . As a quantitative indicator of phylogenetic signals, the MP trees exhibited an average CI of 0. 63 and an average RI of 0. The relatively low CI may indicate the presence of homoplasy, which reduces phylogenetic signal through convergent evolution . Long-branch attraction is also assumed to promote homoplasy, as the MP algorithm may incorrectly group distant taxa that have accumulated more mutations . High RI values, however, indicate synapomorphic signals in the DNA sequences, providing evolutionary information for MP analysis . In this study. RI values were close to 1, with 0. 75 in ITS as the lowest, suggesting that many synapomorphies were retained in the Lamiaceae DNA sequences. Therefore, although the MP trees showed relatively low CI values, the consistently higher RI values supported the robustness of the MP model in recovering reliable evolutionary relationships for this dataset. Bootstrap support was assessed with 1000 replicates, where values of Ou70% are associated with a 95% probability of the true evolutionary relationships . The MP trees based on individual barcodes showed lower support, with 44% of nodes in the rbcL tree, 68% in the matK tree, and 76% in the ITS tree exceeding 70%. In contrast, 88% of nodes in the combined Table 2. Characteristics of the Lamiaceae partial DNA barcodes. DNA Barcodes Sample Size Length Range . Average Length A SD Post-trimming Length . Homology (%) Parsimony -informative (P. Sites (%) ITS 478 - 699 5 A 39. 679 - 866 6 A 35. 523 - 907 5 A 71. ITS matK rbcL 478 - 907 2 A 94. *Note: bp = base pair DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Figure 1. Maximum parsimony phylogenetic tree based on the partial ITS barcode. The numbers on the branches represent bootstrap values from 1000 replicates. DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Figure 2. Maximum parsimony phylogenetic tree based on the partial matK barcode. The numbers on the branches represent bootstrap values from 1000 replicates. DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Figure 3. Maximum parsimony phylogenetic tree based on the partial rbcL barcode. The numbers on the branches represent bootstrap values from 1000 replicates. DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 tree were over 70%. Similarly, 92% of nodes in the ML tree were above 70%, displaying congruence between both models. Individual barcodes demonstrated similar groupings for closely related genera. Consistent clusterings were observed between Orthosiphon and Ocimum. Thymus. Origanum, and Mentha. and Lamium with Stachys. However, three grouping patterns were not recovered in the rbcL tree: Salvia as the sister group to Thymus. Origanum, and Mentha. Clerodendrum as sister to Lamium and Stachys. Clerodendrum. Lamium, and Stachys as sister to Scutellaria. These results suggest conserved sequences in the dataset, as reflected by the higher homology relative to Pi sites. Despite these similarities, individual barcodes did not consistently place every species with its correct genus. In contrast, combining the three barcodes resolved these inconsistencies, grouping Lamiaceae into six robust monophyletic clades in both the MP and ML trees (Figures 4 and . Minor differences remained in the placement of one to two Origanum. Scutellaria, and Callicarpa species. Nevertheless, higher bootstrap values and consistent overall topologies across both trees indicated that the improvement resulted from characters, which increased branch support and phylogenetic resolution . Therefore, only the combined trees are discussed further. Five clades were classified into four monophyletic subfamiliesAiNepetoideae. Ajugoideae. Lamioideae, and ScutellarioideaeAiconsistent with one of the most extensive phylogenetic studies of Lamiaceae . Nepetoideae is sister to the lineage comprising Ajugoideae. Scutellarioideae, and Lamioideae . Callicarpa was separated from the rest of Lamiaceae in Clade I, and together with the four subfamilies, form the phylogenetic backbone of Lamiaceae . 1 Callicarpa: Incertae sedis Callicarpa (Clade I) is considered an incertae sedis taxon due to its unresolved phylogenetic placement within Lamiaceae . Previous research based on five chloroplast DNA identified Callicarpa as the sister group to the subfamily Prostantheroideae . In contrast, analysis incorporating chloroplast and mitochondrial DNA placed it as a sister to the rest of the family . The present study also maintains its incertae sedis status due to incomplete sampling and limited genomic representation, which have contributed to inconsistencies in its phylogenetic position. In the MP tree. Callicarpa diverged from the main branches leading to the remaining clades, while in the ML tree, it diverged from Clades II and i. Morphologically. Callicarpa is characterized by its peltate or capitate stigma and drupes containing four stony pyrenes . 2 Subfamily Nepetoideae Nepetoideae encompassed several genera distributed across clades IV Ae VI. As the largest subfamily within Endlessa et al. Lamiaceae, it comprises approximately 118 genera divided into three tribes: Mentheae. Ocimeae, and Elsholtizeae . , . Clades IV and V formed a monophyletic clade within Mentheae, while clade VI belonged to Ocimeae. Clade IV contained the paraphyletic genus Salvia . ubtribe Salviina. , which is characterized by two fertile stamens and a distinctive staminal lever mechanism in pollination . Salvia was initially classified as monophyletic based on its staminal morphological traits . However, a DNA-based phylogenetic analysis has redefined it as part of a larger clade that includes multiple genera, revealing that some Salvia species are more closely related to other genera than each other . Consequently. Salvia and its related genera are now recognized as a paraphyletic group. Clade V was one of the largest clades, comprising Thymus. Origanum, and Mentha of the subtribe Menthinae. Members of this group are distinguished by reticulate pollen grains and a circular abscission scar . In the study. Thymus and Origanum each formed monophyletic groups clustered as sister clades, with Mentha recovered as their closest relative. However, earlier chloroplast DNA analysis reported Thymus as paraphyletic to Origanum, pointing to possible introgression rather than strict shared ancestry . Furthermore. Mentha was recovered here as paraphyletic, contrary to previous work identifying it as monophyletic using ITS, trnK, and trnL-trnF barcodes . Clade VI included the monophyletic genera Orthosiphon and Ocimum, members of the tribe Ocimeae identified by dorsifixed anthers . This close relationship has also been supported by previous studies . Although Orthosiphon has been extensively investigated for its therapeutic essential oils, phylogenetic research remains relatively scarce compared to other genera, and available DNA sequence data in GenBank are limited. Thus, this study did not resolve additional evolutionary relationships within Orthosiphon, yielding results consistent with Sudarmono et al. Meanwhile. Ocimum was classified as monophyletic, despite earlier research suggesting polyphyly . , likely due to the absence of clear synapomorphies and high morphological similarities among species . , . Overall, discrepancies between findings of the current and prior studies may reflect differences in DNA barcodes and the scope of taxon sampling. 3 Subfamily Ajugoideae and Lamioideae Clade i comprised the genera Clerodendrum. Lamium, and Stachys. Clerodendrum was recovered as monophyletic, forming a sister relationship with the latter genera. It is a member of the subfamily Ajugoideae, the third largest in Lamiaceae, consisting of 23 genera . Previous phylogenetic analysis of Ajugoideae using four chloroplast markers classified it into four clades . , later revised into tribal status, with Clerodendrum placed within Clerodendreae . Although previously considered polyphyletic based on DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 plastid data . , nuclear DNA has redefined Clerodendrum as monophyletic, except for two Australian and Indian Ocean species more closely related to the genus Volkameria . The same study also identified three major lineages, including an Asian-Australian clade . Although the present sampling included only a subset of the Asian clade, the combined ITS, matK, and rbcL sequences recovered an evolutionary pattern consistent with Satthaphorn et al. These barcodes thus provide valuable information for reconstructing phylogenetic relationships within Clerodendrum. Synapomorphies of this genus include brightly colored accrescent calyces and protandrous reproductive organs . Lamium and Stachys are categorized within Lamioideae, in the tribes Lamieae and Stachydeae, respectively . Members of this subfamily are characterized by tricolpate pollen grains and spatulate embryos . This analysis resolved Lamium as non-monophyletic, with L. album separated from a clade containing L. amplexicaule and L. Bendiksby et . , using nuclear and chloroplast DNA from 79 Lamium samples, generally supported the monophyly of the genus, though L. galeobdolon demonstrated variable placements, and L. amplexicaule appeared polyphyletic in the chloroplast Lamium is divided into three subgenera, with L. album and L. amplexicaule belonging to subgenus Lamium, which is further subdivided into different sections (L. in Lamiotypus and L. amplexicaule in Amplexicaul. The non-monophyly of Lamium observed here may reflect these section-level divergences within subgenus Lamium and the uncertain position of L. galeobdolon, which has been assigned within the genus or to separate genera . Stachys was found to form a monophyletic group, with members of Stachydeae often characterized by campanulated calyces and strongly two-lipped corollas . However, it is recognized as a paraphyletic genus including several smaller genera nestled within . As a result, this taxonomic complexity complicates the phylogenetic placement of Stachys, as its species are distributed across Stachydeae. The present study grouped Lamium and Stachys as sisters to Clerodendrum in Clade i, with Lamioideae and Ajugoideae forming a sister group to Scutellarioideae. This topology differs from Scheen et al. , who reported Lamioideae as more closely related to Scutellarioideae. Another large-scale analysis using four chloroplast markers . revealed similar findings. Such topological incongruence among studies is likely caused by variations in DNA barcodes and sampling. Expanding taxon representation remains one of the most effective strategies to improve the accuracy of Lamiaceae phylogeny, particularly at the subfamily level, as demonstrated in this case . Endlessa et al. 4 Subfamily Scutellarioideae Clade II consisted of a single genus. Scutellaria, within the subfamily Scutellarioideae. It formed a monophyletic group, consistent with the findings of a previous study using three chloroplast DNA regions . Scutellaria is distinguished by a unique morphological structure called the scutellum, a projecting appendage on the upper lip of the bilabiate calyx . Based on inflorescence and bract characteristics. Scutellaria has been classified into two subgenera, subg. Scutellaria and Apeltanthus . In Salimov et al. , these subgenera were divided into three clades: subg. Apeltanthus formed a clade with several species of subg. Scutellaria, while the other two clades consisted solely of subg. Scutellaria. Thus, the latter subgenus was inferred to be paraphyletic to subg. Apeltanthus. However, this paraphyly was not detected in the present study, as sampling only included species from subg. Scutellaria. This study identified Lamioideae and Ajugoideae as the sister group to Scutellarioideae, in contrast to several earlier DNA-based analyses mentioned in the previous subsection. Morphological evidence provides additional support for a closer relationship between Lamioideae and Scutellarioideae. In particular, calyx structure and xylem tissue characteristics indicate that species in both subfamilies generally possess a higher density of fibre cells, tracheids, and vessels in the calyx tube compared to those in other Lamiaceae subfamilies, especially Nepetoideae . These findings are consistent with the Lamiaceae phylogenetic backbone proposed by Li et al. , supporting also a closer evolutionary relationship between Lamioideae and Scutellarioideae than with Ajugoideae. Lamiaceae Phylogenetics: Species Authentication and Adulteration Prevention The use of plants in the treatment of various diseases has long been integral to human health. Approximately 80% of the global population relies on plants as their primary source of traditional healthcare, largely due to their perceived lower side effects compared to conventional medicines . , . However, the increasing demand for plant-based remedies has also led to rising cases of adulteration, including within the Lamiaceae The phylogenetic reconstruction of Lamiaceae using the concatenated dataset presented in this study provides a valuable framework for addressing these issues, particularly in species authentication to prevent adulteration. The tree topologies derived from the combined sequence revealed relationships among taxa that support accurate species authentication. These phylogenies can highlight closely related species that are potential adulterants, thereby guiding the development of DNA markers for DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Figure 4. Maximum parsimony phylogenetic tree based on the concatenated partial ITS, matK, and rbcL barcodes. The numbers on the branches represent bootstrap values from 1000 replicates. The bolded words on the far right indicate subfamilies and an incertae sedis taxon. DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Figure 5. Maximum likelihood phylogenetic tree based on the concatenated partial ITS, matK, and rbcL barcodes. The numbers on the branches represent bootstrap values from 1000 replicates. The bolded words on the far right indicate subfamilies and an incertae sedis taxon. DOI: 10. 5614/3bio. 3Bio Journal of Biological Science. Technology and Management 8. : 28-43 Endlessa et al. Table 3. Consistency and retention indices of Lamiaceae maximum parsimony phylogenetic trees. DNA Barcodes Consistency Index (CI) Retention Index (RI) ITS ITS matK rbcL reliable identification of authentic species and their likely substitutes in raw materials and finished products. Such authentication is crucial, as adulteration may compromise the purity and efficacy of herbal products . Furthermore, considering that different Lamiaceae taxa possess distinct secondary metabolite profiles, adulteration may also introduce toxic compounds that raise safety concerns . The issue is exemplified by Scutellaria baicalensis, widely used in traditional cancer treatments for its wogonin content . High market demand has resulted in the adulteration of its commercial products with S. Similarly, peppermint oil from Mentha piperita, which contains menthone and menthofuran that are beneficial against oxidative stress and inflammation . , . , is frequently substituted with menthol oil from M. These substitutions may involve differences in chemical compositions and affect therapeutic Therefore. DNA-based phylogenetic analysis enables reliable species authentication and adulterant detection. Conclusion In the present study, individual partial DNA barcodes (ITS, matK, and rbcL) consistently grouped closely related genera, confirming their ability to reconstruct relatively accurate Lamiaceae phylogenies. Nevertheless, several species were separated from their respective genera. Concatenation of the three barcodes resolved these discrepancies, yielding improved phylogenetic resolution further supported by the ML tree model. The combined analysis clustered Lamiaceae into six well-supported monophyletic clades, largely consistent with previous studies, although some differences in phylogenetic placement were observed. These findings highlight the utility of multilocus DNA barcodes for clarifying evolutionary relationships in Lamiaceae and supporting species authentication and adulteration prevention in herbal Future research should expand taxon sampling, including but not limited to more incertae sedis genera and geographically diverse species. Incorporating longer DNA sequences, additional markers such as mitochondrial apocytochrome b (COB), and diverse phylogenetic tree models may further refine phylogenetic resolution. References