Next Article in Journal
microRNAs as Early Biomarkers of Alzheimer’s Disease: A Synaptic Perspective
Next Article in Special Issue
Metabolomics Intervention Towards Better Understanding of Plant Traits
Previous Article in Journal
Cytokine-Mediated Crosstalk between Immune Cells and Epithelial Cells in the Gut
Previous Article in Special Issue
Metabolite Profiling of Alangium salviifolium Bark Using Advanced LC/MS and GC/Q-TOFTechnology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomics and DNA-Based Authentication of Two Traditional Asian Medicinal and Aromatic Species of Salvia subg. Perovskia

1
Department of Pharmaceutical Biotechnology, Wroclaw Medical University, Borowska 211A, 50-556 Wroclaw, Poland
2
Botanical Garden of Medicinal Plants, Wroclaw Medical University, Jana Kochanowskiego 14, 51-601 Wroclaw, Poland
3
Department of Horticulture, College of Agriculture, Isfahan University of Technology, Isfahan 841583111, Iran
4
Core Research Facility (CRF), Isfahan University of Medical Sciences, Hezar Jerib Street, Isfahan 8174673461, Iran
5
Department of Pharmaceutical Biology and Botany, Wroclaw Medical University, Borowska 211, 50-556 Wroclaw, Poland
*
Author to whom correspondence should be addressed.
Cells 2021, 10(1), 112; https://doi.org/10.3390/cells10010112
Submission received: 25 November 2020 / Revised: 2 January 2021 / Accepted: 7 January 2021 / Published: 9 January 2021
(This article belongs to the Special Issue Metabolomics in Plant Research)

Abstract

:
Subgenus Perovskia of the extended genus of Salvia comprises several Central Asian medicinal and aromatic species, of which S. yangii and S. abrotanoides are the most widespread. These plants are cultivated in Europe as robust ornamentals, and several cultivars are available. However, their medicinal potential remains underutilized because of limited information about their phytochemical and genetic diversity. Thus, we combined an ultra-high performance liquid chromatography quadrupole time of flight mass spectrometry (UHPLC-QTOF-MS) based metabolomics with DNA barcoding approach based on trnH-psbA and ITS2 barcodes to clarify the relationships between these two taxa. Metabolomic analysis demonstrated that aerial parts are more similar than roots and none of the major compounds stand out as distinct. Sugiol in S. yangii leaves and carnosic acid quinone in S. abrotanoides were mostly responsible for their chemical differentiation, whereas in roots the distinction was supported by the presence of five norditerpenoids in S. yangii and two flavonoids and one norditerpenoid in S. abrotanoides. To verify the metabolomics-based differentiation, we performed DNA authentication that revealed S. yangii and S. abrotanoides to be very closely related but separate species. We demonstrated that DNA barcoding coupled with parallel LC-MS profiling constitutes a powerful tool in identification of taxonomically close Salvia species.

Graphical Abstract

1. Introduction

The subgenus Perovskia Kar. comprises several species with an ethnobotanical value and medicinal potential. These shrubby, aromatic plants are widespread throughout Central Asia, in an area that includes Iran, Pakistan, Afghanistan, Turkmenistan, Uzbekistan, and Tibet. The ranges of most of the accepted species are limited to the arid montane regions in Central Asia, but two of them, Salvia abrotanoides Karel. (previously known as Perovskia abrotanoides Kar.) and Salvia yangii B.T. Drew (formerly Perovskia atriplicifolia Benth.), are distributed over vast areas from Western Iran and Pakistan and as far as Tibet and Xinjiang, in China [1,2,3,4,5,6,7,8,9].
There are a lot of discrepancies regarding the taxonomy of the subgenus Perovskia. According to primary classification based on flower morphology, species belonging to Perovskia together with four other genera, Dorystaechas, Meriandra, Rosmarinus, and Zhumeria, have previously been considered closely related to Salvia and were treated together as part of the subtribe Salviinae in the most recent and comprehensive conspectus of Lamiaceae [10].
Later, the field of molecular phylogenetics fundamentally altered the basis of classification of organisms. In 2004, Walker et al., using the rbcL and trnL-F gene regions and broad sampling within Salvia, were the first to reveal that Salvia was non-monophyletic [11]. They found Perovskia, Dorystaechas, and Rosmarinus, as well as Mentha L., Origanum L., and Thymus L., and later also Meriandra and Zhumeria embedded within Salvia [11,12].
In 2017, a broad phylogenetic study was undertaken using nuclear markers such as the internal transcribed spacer from nuclear ribosomal DNA (nrITS) and plastid marker rpl32-trnL [13]. This study confirmed the polyphyletic character of Salvia but also showed the existence of four distinct, well supported evolutionary lineages (Clades I-–IV) within Salvia and the inclusion of five additional non-Salvia genera, namely, Dorystaechas, Meriandra, Perovskia, Rosmarinus, and Zhumeria. Based on the plastid data, Perovskia and Rosmarinus were found in a trichotomy with Clade I (Salvia s.s.). The authors reconstructed the evolution and biogeography of Salvia and found that the most recent common ancestor of Perovskia and Rosmarinus, which are sister to Salvia Clade I, diverged most likely during the early to middle Miocene. Based on current knowledge, Will and Claßen-Bockhoff propose to split the former large genus Salvia into six genera, each supported by geographical distribution, morphology, and karyology, and suggest that they could be elevated to generic rank after revision.
In a related study, low-copy nuclear markers (PPR-AT3G09060, GBSSI) were used to further assess relationships of Salvia and related genera within Salviinae [14]. Based upon these results, previous phylogenetic findings, taxonomic, morphological, and practical considerations, the authors concluded that the botanical community would be better served by keeping the genus Salvia as traditionally circumscribed, but with the inclusion of 15 species in the five small embedded genera: Perovskia, Dorystaechas, Meriandra, Rosmarinus, and Zhumeria. Subsequently, to implement this approach and unify all genera within Salvia, the authors provided an updated nomenclatural revision, according to which former Perovskia abrotanoides is now circumscribed as Salvia abrotanoides (Kar.) Sytsma and former Perovskia atriplicifolia as Salvia yangii B.T. Drew.
Salvia yangii B.T. Drew, commonly called Russian sage, is a silver-grey herbaceous perennial half-shrubby plant. All parts of S. yangii have silvery hairs, including the flower heads. This plant is native to Central Asia in an area that includes Afghanistan, Eastern Iran, Pakistan, and Tibet [1,15]. Successful over a wide range of climate and soil conditions, S. yangii has become popular and is widely planted. Several cultivars/varieties have been developed and cultivated worldwide for ornamental purposes, mainly due to the bountiful blue to violet blossoms arranged in branched panicles.
Salvia abrotanoides Kar., called Caspian Russian sage, is also an aromatic erect shrub that occurs in mountains at an elevation of 2200 to 4200 m above sea level from Central and Northeastern Iran across Northern Pakistan to Northwestern India [7].
Both species are closely related and may share very similar morphology, except the leaf shape, which is the main determinative feature, according to many identification keys [1,8,16]. Numerous publications also inform that the characters of the inflorescence, flower, and seeds are similar in both species and they only differ in leaf shape [17,18,19]. The leaf shape is not a reliable determination feature though, as it is very variable within the species, especially in S. yangii.
Species of subgenus Perovskia range from southwestern to central Asia, with S. abrotanoides being more distributed to the west and S. yangii to the east. The ranges of the two species overlap in Northeastern Iran, Afghanistan, Pakistan, and Nepal [20]. It is suggested that both species can hybridise in natural habitats [1]. Plants grown from seeds originating from botanical gardens or other plant collections where both species grow in short distance from each other also might be crossed. As both species are able to hybridise, it is not clear which cultivars are variants of S. yangii or hybrids between S. abrotanoides and S. yangii [18,19]. This makes the seemingly simple determination of Perovskia species unexpectedly complicated.
Both S. yangii and S. abrotanoides are used as a folk medicinal herb in the areas of natural occurrence. S. yangii is a folk medicinal plant used in Pakistan, Turkmenistan, and Afghanistan [21,22]. The plant has antibacterial activity and is also used as a cooling medicine in the treatment of fever [3,23]. It is known for its flavouring qualities and it is also smoked as a euphoriant and used as a decoction for chronic dysentery [21]. In Uzbekistan and Kyrgyzstan, the aboveground parts are used to heal wounds, and a decoction is used to treat scabies [17]. In traditional Tibetan and Chinese medicine, it has been claimed that S. yangii is a powerful analgesic and parasiticide agent [9,24].
S. abrotanoides, locally known as brazambol, is used to treat leishmaniasis in Iranian folk medicine. Villagers in the Isfahan Province apply a poultice, made of crushed roots of the plant, water, sesame oil, and wax, on lesions caused by cutaneous leishmaniasis [25,26]. The plant is also used locally for treatment of typhoid, fever, headache, gonorrhoea, vomiting, motion, toothache, cardiovascular diseases, liver fibrosis, painful urination, cough, and as a sedative, analgesic and antiseptic drug [21,27,28,29,30]. Also, the essential oil of S. abrotanoides is effective as an antimicrobial agent [31].
Species of subgenus Perovskia contain a wide array of natural compounds belonging to different classes [22]. Most of the published data refer to the composition and activities of essential oils from aerial parts [32,33,34,35,36,37,38,39,40,41], indicating a significant variability associated with geographic origin and environmental factors Aerial parts contain also lignans, triterpenes, steroids, and their glycosides, as well as phenolic compounds, including considerable amounts of rosmarinic acid, which is a typical product of the Nepetoideae subfamily of Lamiaceae, known for its numerous health-promoting properties [2,3,7,15,42,43].
Roots of species belonging to the subgenus Perovskia are abundant in red-colored quinoid norabietanoids called tanshinones and are considered an alternative source of these valuable pharmacologically active compounds [25,44]. Tanshinones are the main constituents of the well-known Traditional Chinese Medicine herbal drug ‘Danshen’ (Salvia miltiorrhiza rhizome et radix a pharmacopoeial drug of China and Europe) [45,46,47,48].
The research performed to date on plants of subgenus Perovskia concentrated on either roots [25,44] or aerial parts [7,49,50,51]. There are also reports on phytochemical analysis of the air-dried whole plants [2,3,4,9,24,52]. Plants were collected from their natural habitat or were cultivated. Various methods and solvents were used by researchers to extract natural products of interest which often makes the results hardly comparable. Therefore, we aimed at comprehensive analysis of phytochemical profiles of non-volatile compounds in both roots and leaves of S. yangii and S. abrotanoides cultivated in Europe under identical climatic and edaphic conditions (Wroclaw, Poland). The ethnopharmacological importance of S. yangii and S. abrotanoides has prompted us to reanalyse their chemical constituents with potential or proven bioactivity. Here, for the first time, we offer a detailed comparison of phytochemical non-volatile profile in roots and leaves of S. yangii and S. abrotanoides together with the analysis of a mutual genetic relationship between the two species. The proposed approach of a comparative analysis of closely related plants of a similar but not identical chemical profile may bring tangible results which might potentially pinpoint the outlook for follow-up work through genetic approaches, transcriptomics, proteomics, and other applications. Choosing non-model but potentially important medicinal plants allowed us to test this methodology and prove that the slight differences in morphological and metabolic phenotypes can be verified by molecular analyses, thus facilitating the proper material identification and accelerating the selection of desired traits, breeding, and herb quality control and improvement.

2. Materials and Methods

2.1. Plant Material

Salvia abrotanoides and S. yangii plants were grown at the experimental field in the Botanical Garden of Medicinal Plants at the Wroclaw Medical University (Poland). The geographical location of the Garden is as follows: latitude 51°07′03″ N, longitude 17°04′27″ E, altitude 117 m a.s.l. S. abrotanoides was introduced from its original habitat in Iranian Karkas Mountains (coordinates: 33°24′57.3″ N 51°46′45.4″ E). Seeds of S. yangii were obtained from the collection of the Botanical Garden in Essen (Germany). Cultivars “Lacey Blue” and “Blue Spire” are available commercially as ornamental plants. Plants analyzed in this study were cultivated clonally via vegetative propagation since introduced into the Botanical Garden of Medicinal Plants at the Wroclaw Medical University. Plant voucher specimens were deposited in the herbarium of the Botanical Garden of Medicinal Plants at the Wroclaw Medical University (reference numbers: S. abrotanoides: P-122; S. yangii: P-123; cultivar “Blue Spire”: P-124; cultivar “Lacey Blue”: P-125).
Leaves and roots of S. abrotanoides and wild type S. yangii (Appendix A, Figure A1) were collected for phytochemical analysis and were immediately frozen in liquid nitrogen. Harvested leaf samples were finely ground using a mortar and pestle with liquid nitrogen. Root samples were finely ground using an automatic/electric grinder Analysette 3 Spartan and Pulverisette 0 Cryo-Box (Fritsch GmbH, Idar-Oberstein, Germany) filled with liquid nitrogen.
Three accessions of S. yangii (including wild type and cultivars “Lacey Blue” and “Blue Spire”) and a single accession of vegetatively propagated S. abrotanoides were subject to the molecular study. For phylogenetic analysis, a single, fresh, mature leaf of each plant was harvested and ground to a fine powder in liquid nitrogen.
All samples were stored at −80 °C. The harvest of leaves and roots was done during the flowering season (August 2016).

2.2. Sample Preparation for Qualitative Analysis

For direct qualitative analysis, 100 mg of lyophilized plant material was extracted with 5 mL of MeOH: H2O (8:2 v/v) in an ultrasonic bath for 20 min at 25 °C. The supernatant was collected and evaporated to dryness. Extracts were suspended in 2 mL of 0.2% formic acid (FA) in H2O prior to preliminary purification and concentration using solid phase extraction (SPE). SPE was performed using an Oasis Max 1cc Vac Cartridge, 500 mg Sorbent per Cartridge, 30 µm Particle Size (Waters Corporation, Milford, MA, USA) conditioned with MeOH (4 mL) and H2O (4 mL). The sample was loaded and cartridges washed with 5 mL water in 0.1% FA. The phenolic compounds were eluted with 5 mL of 80% MeOH in 0.1% FA. Samples were re-evaporated to dryness and dissolved in 25% of MeOH to obtain 1 mg dry weight/mL baseline plant material concentration.

2.3. Liquid Chromatography–Mass Spectrometry Analysis

Ultra-high performance liquid chromatography quadrupole time of flight mass spectrometry (UHPLC-QTOF-MS) estimation of metabolite contents was carried out on a Thermo Ultimate 3000 RS (Thermo Fischer Scientific, Waltham, MA, USA) chromatographic system coupled to a Bruker Compact (Bruker, Billerica, MA, USA) quadrupole time of flight (QTOF) mass spectrometer, consisting of a binary pump, sample manager, column manager, and the PDA/diode array detector. Separations were performed on a Kinetex C18 column (2.1 × 100 mm, 2.6 μm, Phenomenex, Torrance, CA, USA), with mobile phase A consisting of 0.1% (v/v) FA in water and mobile phase B consisting of 0.1% (v/v) FA in acetonitrile. The samples (2 μL) were injected, and a linear gradient from 15% to 70% phase B in phase A was used for 27 min to separate all main compounds. The flow rate was 0.3 mL/min, and the column was held at 30 °C. Mass spectra were acquired in positive-ion mode over a mass range from m/z 100 to 1500 with 5 Hz frequency. Ultrapure nitrogen was used as drying and nebulizer gas and argon was used as collision gas. The collision energy was set automatically from 15 to 75 eV, depending on the m/z of the fragmented ion. Operating parameters of the ESI ion source were as follows: capillary voltage 3 kV, dry gas flow 6 L/min, dry gas temperature 200 °C, nebulizer pressure 0.7 bar, collision radio frequency 700.0 V, transfer time 100.0 μs, and pre pulse storage 7.0 μs. Acquired data were calibrated internally with sodium formate introduced to the ion source at the beginning and end of each separation via a 20 μL loop. Processing of spectra was performed with Bruker DataAnalysis 4.3 software. UV DAD data were collected at 190–450 nm, and chromatograms acquired at UV 280 nm were smoothed by use of the Savitzky−Golay algorithm (window width 5 points, one iteration). Then, the detected peaks at retention time and the corresponding MS/MS were integrated. All determinations were done in triplicate (biological replication meaning a different extract n = 3); each biological replication was measured at least twice as a technical replication.
ProfileAnalysis software (version 2.1, Bruker Daltonik GmbH, Bremen, Germany) was used to preprocess the raw UHPLC-QTOF-MS data. ProfileAnalysis parameters were set as follows: advanced bucket generation with retention time range of 1.0–24.0 min, mass range of 100–800 m/z, normalization to sum of peaks, background subtraction, and time alignment; then, analyses were processed with the Find Molecular Futures (FMF) function to create compounds with a signal-to-noise threshold of 3 for peak detection. The generated bucket table consisting of m/z–retention time pairs and respective compound intensity was exported and uploaded to Metabolists 5.0 (http://www.MetaboAnalyst.ca/) for multivariate statistical analysis. The initial and final retention times were set for data collection. Data were imported into Metabolists 5.0 (http://www.MetaboAnalyst.ca/) online software to estimate missing values and to filter and normalize data (normalization by the median). No transformation was generalized, and the data matrix was mean-centered and divided by the square root of the standard deviation of each variable (Pareto scaling). The obtained data matrix was introduced into SIMCA-P+ 16.01 (Sartorius AG, Göttingen, Germany) software for in-review multivariate statistical analysis of principal component analysis (PCA). The PCA score plot was used to present a natural correlation between the observations. To identify differential compounds, the PLS-DA (Partial Least Squares Discriminant Analysis) model was used to explore differences in depth between the profile metabolome of Salvia yangii and S. abrotanoides The PLS-DA model with VIP values (VIP ≥ 1.0) and |p(corr)| ≥ 0.5 was used to select differential compounds (PLS-DA data not shown).

2.4. Genomic DNA Extraction

DNA was extracted with the modified Doyle & Doyle method [53]. Spectrometric DNA concentration measurements were conducted with a NanoDrop ND 2000C (Thermo Fisher Scientific, Waltham, MA, USA). If needed, DNA isolates were purified with the Syngen DNA clean-up Kit (Syngen Biotech, Wroclaw, Poland) and the DNA concentration was measured again.

2.5. DNA Barcoding

ITS rDNA regions were amplified using primers ITS1F [54] and ITS4 [55] while the trnH-psbA chloroplast intergenic spacer was amplified with trnHf_05 [56] and psbA3_f [57] primers (Table 1). PCR was performed in a T-100 thermal cycler (BIO-RAD Laboratories, Hercules, CA, USA) in 20 μL volumes containing template DNA and was set up according to the Q5™ High-Fidelity DNA Polymerase manufacturer’s protocol (Q5™ High-Fidelity DNA Polymerase, New England BioLabs Ltd., Hitchin, UK). Thermal cycling conditions were as follows: 95 °C for 30 s; 35 cycles of 95 °C for 10 s; 60 °C (for ITS1F&ITS4) or 62 °C (for trnHf_05&psbA3_f) for 30 s; 72 °C for 45 s and 72 °C for 2 min for final extension. PCR products were visualized by gel electrophoresis in 2% agarose gel in TEA buffer with 50 bp DNA Ladder (New England Biolabs Ltd., Hitchin, UK) and stained with SimplySafe (EURx, Gdansk, Poland). Amplified DNA products were purified with the Syngen Gel/PCR Mini Kit (Syngen Biotech, Wroclaw, Poland).
Sanger sequencing was carried out with the BrilliantDye™ Terminator v3.1 Kit (Nimagen B.V., Nijmegen, The Netherlands). Reactions were set with a fourfold dilution of the reaction premix and the addition BrilliantDye® Terminator 5X Sequencing Buffer (Nimagen B.V., Nijmegen, The Netherlands) according to the producer’s instructions. Sequencing products were precipitated with ethanol, dissolved in TSR (Hi-Di Formamide) (Thermo Fisher Scientific, Waltham, MA, USA) and then separated by capillary electrophoresis on an Applied Biosystems™ 310 Genetic Analyzer (Thermo Fisher Scientific, Waltham, MA, USA). Two reads were collected for each sample. Eight sequences were generated for molecular and phylogenetic analyses performed along with additional sequences obtained from NCBI (https://www.ncbi.nlm.nih.gov).
Sequence quality was checked with the Sanger Quality Check App (Thermo Fisher Scientific, Waltham, MA, USA). Forward and reverse sequencing reads for each marker region were assembled into a contig using the BioEdit Sequence Alignment Editor [58]. The identity of all obtained sequences was confirmed as ITS or trnH-psbA through similarity to published sequences using the BLAST algorithm. The ITS2 fragment was extracted with ITSx [59].
Alignments were assembled separately for each marker with the ClustalW algorithm under default conditions in MEGA 7 [60] and were then entered into the program Gblocks v.091b to objectively eliminate poorly aligned positions and divergent regions [61]. Then, a final alignment was assembled by a combination of both alignments calculated for individual markers. The relationship among the taxa was estimated using the Neighbor-Joining (NJ) method [62] and presented as an -rooted phylogram. Due to the large number of gaps in the dataset, a Pairwise deletion Gaps/Missing data treatment was applied. The evolutionary distances were computed using the Maximum Composite Likelihood method [63] according to the general recommendation. To estimate reliability of the tree bootstrap, a 1000-replication test was used [64]. The analysis involved 29 nucleotide sequence combinations including L. maculatum and L. amplexicaule as an outgroup.

3. Results

3.1. Metabolic Analysis

Using UHPLC-QTOF-MS analysis, thirty-eight different compounds were found in the methanolic extract of S. yangii leaves, while leaves of S. abrotanoides comprised thirty-six metabolites extracted with methanol (Table 2; Figure 1). The analysis showed the presence of 14 different phenylpropanoids (3–6, 8–11,14,15,17,18, 21, 27) in leaves of both S. yangii and S. abrotanoides, including rosmarinic acid (11), considered as one of the chemotaxonomic markers in the Nepetoideae subfamily [22]. In total, 13 diterpenoids were detected in S. yangii leaves (20, 22, 23, 28–30, 32, 34, 35, 37, 38, 40, 41) while leaves of S. abrotanoides contained 15 different diterpenoids, out of which carnosic acid quinone (31), isorosmanol (33), and trilobinol (39) were found in S. abrotanoides and not in S. yangii. On the other hand, rosmaridiphenol (35) and sugiol (38) were detected exclusively in S. yangii leaves. Among all detected diterpenoids, carnosol (30) turned out to be the most abundant in the leaves of both species as it yielded the highest peak (Figure 1). Both species were found to contain a phenantrene, methoxy-8,11,13-abietatrien-20,11-olide (36), in their leaves. Interestingly, three metabolites identified as triterpenoids (24–26) were detected exclusively in leaves of S. yangii. Seven compounds found in leaves of both species were not determined (1, 2, 7, 12, 13, 16, 19).
UHPLC-QTOF-MS analysis of methanolic extracts of S. yangii roots revealed the presence of thirty-five different compounds, while in roots of S. abrotanoides, only twenty-six metabolites were detected (Table 3; Figure 2). In roots of S. yangii, six phenylpropanoids were found (1, 3, 5, 9–11), out of which three, namely, methoxytaxifolin (1), 2-O-p-coumaroyltartronic acid (3), and melitric acid B (10), were found only in S. yangii and not in S. abrotanoides. Roots of S. abrotanoides contained four phenylpropanoids, one of which was found exclusively in this species: a flavonoid glucoside–hesperidin (7). Rosmarinic acid (9) and salviaflaside (5), a rosmarinic acid 3′ glucoside, were detected in both species as well as the 8-methylchromen-4-one (12), a phenolic compound belonging to benzopyrans, and salvianolic acid L (6), belonging to phenylnaphthalenes. Chrysoeriol 7-rutinoside (8), which belongs to fatty alcohols/aliphatic alcohols, was detected exclusively in the roots of S. abrotanoides. The most abundant group of non-volatile metabolites in roots was diterpenoids, represented by 18 compounds in S. yangii and by 13 metabolites in S. abrotanoides. Twelve diterpenoids (14, 15, 17–21, 31, 34, 35, 37, 38) were common to both species, including the most abundant, cryptotanshinone (31). Six diterpenoids were unique in the metabolite profile of S. yangii roots, namely: didehydrotanshinone IIa (25), acetyloxycryptotanshinone (26), isograndifoliol (27), grandifoliol (28), didehydroacetyloxycryptotanshinone (30), and ketoisogrndifoliol (32). An undetermined cryptotanshinone derivative (29) was found exclusively in S. abrotanoides roots. Nine compounds were not determined in methanolic extracts of roots (2, 4, 13, 16, 22, 23, 24, 33, 36) out of which three were detected only in S. yangii (2, 24, 36).

3.2. Statistical Analysis of UHPLC-QTOF-MS Data

In order to reveal group patterns in non-volatile metabolic profiles of analyzed leaves and roots, an unsupervised PCA was applied. The PC loadings revealed a clear trend to form two clusters among all tested samples, either in leaves or roots, which represent the data in the space formed by the principal components (Figure 3 and Figure 4).
In leaves, the first two principal components (PCs) explained 47.3% and 17.1% of the variation in the spectral data, respectively. On the right side of the score plots (green dots) a strong cluster of S. abrotanoides samples formed, while on the left side (blue dots) S. yangii samples were well separated (Figure 3a). According to the loadings plot (Figure 3b) these two groups were significantly positively correlated with particular compounds detected in leaves. The PC loadings analyses confirmed that the content of sugiol (38) in S. yangii leaves and of carnosic acid quinone (31) in leaves of S. abrotanoides were mostly responsible for their chemical differentiation and revealed a clear separation of samples derived from these two species.
In roots, the first principal component explained 65.3% of the variation in the spectral data and the second, 10.2%. Samples of S. yangii, represented by blue dots, were closely grouped on the right side of the plot, while on the left side (green dots), S. abrotanoides samples occurred (Figure 4a). The loadings plot revealed that these two groups were significantly positively correlated with particular compounds found in roots (Figure 4b). According to the PC loadings, the contents of isograndifoliol (27), acetyloxycryptotanshinone (26), ketoisograndofoliol (32), oxocryptotanshinone (20) and grandifoliol (28) in S. yangii roots and the presence of cryptotanshinone derivative (29), chrysoeriol 7-rutinoside (8), and hesperidin (7) in the roots of S. abrotanoides were mostly responsible for their chemical differentiation and allowed distinguishing the root samples from the two analyzed Salvia species.

3.3. Genetic Relationship Analysis

Sanger sequencing resulted in four rDNA ITS2 complete sequences (MT599312, MT599313, MT599314, MT599315) and four trnH-psbA intergenic spacer sequences (MT815871, MT815872, MT815873, MT815874) obtained for the cultivars S.OBRL.01, S.OBRL.03, S.OBRL.11 and S.OBRL.12, respectively (Table 4). All sequences were subject to BLAST analysis. Due to the lack of publicly available molecular data of the trnH-psbA barcode in subgenus Perovskia, BLAST analysis of the sequenced trnH-psbA region was not able to indicate similarity to any species from the Perovskia subgenus. Instead, the closest homolog of S.OBRL.01 and S.OBRL.11 trnH-psbA regions indicated by BLAST analysis was the trnH-psbA sequence of Salvia miltiorrhiza, while for S.OBRL.03 and S.OBRL.12, Salvia californica and Salvia chionopeplica were indicated as the closest homologs, respectively. BLAST analysis of the abundantly available data in the ITS2 region revealed high similarity between Salvia abrotanoides and Salvia yangii (Table 4).
NJ analysis was performed to estimate the relationship among the studied taxa (S. yangii, including wild type and cultivars “Lacey Blue” and “Blue Spire”, and S. abrotanoides) in comparison with 23 related Salvia species (Table 5). Lamium maculatum and Lamium amplexicaule were applied as an outgroup. Finally, 50 additional sequences used for the NJ analysis were downloaded from GenBank (25 sequences of ITS2 and 25 sequences of trnH-psbA). The manually adjusted alignment of 29 combined ITS2 and 29 trnH-psbA sequences was 540 characters long with 361/540 variable sites, 301/540 Parsimony-informative sites.
According to the phylogram reconstructed based on ITS2 and trnH-psbA regions, the subgenus Perovskia was separated from other sages, yet clearly embedded in Salvia, (Figure 5). The topology of the phylogram exhibiteds low interspecies genetic variation between subgenus Perovskia representatives indicated by short branches of the tree.
Despite the close genetic relationship, S. abrotanoides and S. yangii are suggested to be separate species. Moreover, intraspecies genetic variation within the S. yangii taxon, between “wild type” S. yangii and its cultivars (‘Lacey Blue’ and ‘Blue Spire’), was observed. Furthermore, the intermediate location of S. yangii ‘Blue Spire’ between the two Salvia species suggests that the ‘Blue Spire’ cultivar might have been created through hybridization between S. abrotanoides and S. yangii during the cultivation process.

4. Discussion

In this work, we have provided a comprehensive analysis of the phytochemical profile of non-volatile compounds separately for roots and leaves of S. yangii and S. abrotanoides cultivated in Europe (Wroclaw, Poland). A reliable comparison was possible due to ensuring the same cultivation conditions such as: latitude, weather conditions etc., the same sampling method (immediate freezing in liquid nitrogen), and time (middle of the growing season) as well as the same extraction procedure (MeOH: H2O; 8:2 v/v) and method for chemical analysis (UHPLC-QTOF-MS). Despite ensuring identical cultivation conditions and methodology, a detailed comparison of the phytochemical non-volatile profile in roots and leaves of S. yangii and S. abrotanoides resulted in detecting several metabolites present in one species and absent in the other. The PC loadings analyses showed that sugiol in S. yangii leaves and carnosic acid quinone in leaves of S. abrotanoides were mostly responsible for their chemical differentiation. In roots, isograndifoliol, acetyloxycryptotanshinone, ketoisograndofoliol, oxocryptotanshinone, and grandifoliol in S. yangii as well as chrysoeriol 7-rutinoside, hesperidin, and one unidentified cryptotanshinone derivative in S. abrotanoides allowed for chemical distinguishing of root samples from analyzed species.
As noticed before, most of the published data from aerial parts of species belonging to subgenus Perovskia refer to the composition and activities of essential oils, indicating a significant variation [67]. A thorough review of the chemical profile of essential oils from subgenus Perovskia is provided by Mohammadhosseini et al. [22]. There are very few reports, however, of non-volatile constituents from aerial parts of species belonging to subgenus Perovskia, as a vast majority of them present an analysis of a whole-plant material.
Moreover, the influence of geographic origin and environmental factors influenced the phytochemical profile in both the volatile and non-volatile subset. The extent of such influences in individual cases is hard to assess without extensive experimental research. For example, a light quality modified metabolic profile of cultivated S. abrotanoides [68].
Tarawneh et al. [51] isolated four flavoinoids from dried leaves of S. yangii, which were not present in our isolates, except for the hydroxy-trimethoxyflavone, detected in leaves of both S. yangii and S. abrotanoides. Additionally, we detected trihydroxy-dimethoxyflavone, not found by Tarawneh et al. [51].
Amongst ten compounds found in methanolic extracts of air-dried aerial parts of S. abrotanoides [7], we detected rosmarinic acid and hesperidin, both present in the leaves of S. yangii and S. abrotanoides. Interestingly, in our study, hesperidin was also found to be present in the roots of S. abrotanoides. Two new tertracyclic diterpens, abrotandiol and abrotanone, found by Khaliq et al. [7] were absent in our extracts. The dihydroxy-methoxymethylflavone, however, found in our leaf isolates, might be the same compound as cismaritin found by Khaliq et al. [7]. Confirming that would require additional NMR analysis.
Amongst diterpenoids found in n-hexane and ethyl acetate extracts of aerial parts of S. abrotanoides [69], carnosol, miltirone, rosmanol, and epi-rosmanol were found in methanol leaf extracts of both S. yangii and S. abrotanoides analyzed here. Miltiodiol, however, has been previously isolated from roots [25,44], which was confirmed by our study. Tabefam et al. [69] also detected ferruginol and dehydroferruginol, both belonging to abietane diterpenoids, which were absent in our isolates. Instead, we detected oxoferruginol (sugiol), present in leaves of S. yangii, while leaves of S. abrotanoides contained hydroxyferruginol (trilobinol).
Interestingly, our work has shown that the accumulation of triterpenoids, namely, 2-hydroxy-3-oxo-12-oleanen-28-oic acid, 3-oxo-12,18-ursadien-28-oic acid, and trihydroxy-12-ursen-28-oic acid, takes place exclusively in leaves of S. yangii and not in S. abrotanoides.
In a previous study of S. yangii roots (air-dried crude material, extraction with n-hexane), eight diterpenoids were detected [44]. Extraction with 80% methanol performed in this study allowed for the detection of 18 diterpenoids in S. yangii roots and 15 in roots of S. abrotanoides. 1βOH-cryptotanshinone, cryptotanshinone, tanshinone IIat, and arucadiol (a synonym of miltiodiol) found by Slusarczyk at al. [44] were also detected in our methanolic root extracts in both S. yangii and S.abrotanoides. Miltirone, however, was present only in the methanolic root extract of S. yangii and not in S. abrotanoides, which corroborates previous results [25,44]. Two dehydrogenation products, didehydrotanshinone and didehydromiltirone, which were previously detected in n-hexane root extract of S. yangii, were not present in our methanolic root extract, neither in S. yangii nor in S. abrotanoides. Interestingly, isograndifoliol, reported recently in n-hexane root extract of S. yangii [66] was also detected in the methanolic root extract of S. yangii and was proved to contribute most to chemical differentiation between S. yangii and S. abrotanoides.
Amongst diterpenoids detected in the EtOAc extracts of native dried roots of S. abrotanoides [25], cryptotanshinone, 1βOH-cryptotanshinone, and 1-oxocryptotanshinone were found in methanolic root extracts of both S. yangii and S. abrotanoides. 1-oxomiltirone, found in both previous studies [25,44], was not detected in this work.
The phytochemistry and chemotaxonomic significance of metabolites found in subgenus Perovskia representatives was reviewed by Mohammadhosseini et al. [22].
In this study, we also aimed at analysis of the mutual genetic relationship between S. abrotanoides and S. yangii rather than their taxonomical position within the Salvia subgenera. To our knowledge, no previous work has done that with using DNA barcoding. Usually, a single-entry representative of Perovskia subgenus was used for phylogenetic tree reconstructions of Salvia [11,12,13,14,70]. There are two reports of employing molecular markers (ISSR markers) to investigate genetic variability among populations from subgenus Perovskia growing wild in Iran [71,72]. Heshemifar and Rahimmalek [72] analysed 63 accessions of S. abrotanoides from 16 populations and four accessions of S. yangii. Thirteen ISSR markers were shown to have a relatively high efficiency in the assessment of population diversity, which was in most cases confirmed by morphological classifications. In the study of Pourhosseini at al. [71], twelve S. abrotanoides populations were analyzed and a cultivated population of S. yangii was used as an outgroup. The authors have shown that nine ISSR markers are able to genetically distinguish between all studied populations and that the S. yangii population was clearly separated from S. abrotanoides populations, proving the effectiveness of this DNA fingerprinting method in studying the genetic relationship within the Perovskia subgenus. Despite the suitability of the ISSR markers for studying genetic variation among species within subgenus Perovskia, we have chosen DNA barcoding as it is regarded as a more sophisticated, standardized, and worldwide acknowledged method, recommended for genetic relationship analysis of land plants [73,74,75].
Although the DNA barcoding has been endorsed as the method of choice for land plant species identification, no consensus has been reached on which plant DNA sequences should be regarded as barcodes. The generally low rate of nucleotide substitution in plant mitochondrial genomes precludes the use of cytochrome oxidase 1 (CO1), a potential equivalent for the standard animal barcode, as a universal plant barcode [76,77]. Instead, the search for a plant barcode has involved looking outside the mitochondrial genome and from the outset, many researchers have accepted that multiple markers will be required to obtain adequate species discrimination [78]. Consequently, the Consortium for the Barcode of Life (CBOL) Plant Working Group has recommended a core-barcode consisting of portions of two plastid-coding regions, rbcL+matK, to be supplemented with additional markers as required [75]. Limitations of the rbcL+matK barcodes and idiosyncratic performance of different markers in different taxonomic groups were well discussed by Fazekas et al. [77]. Using additional markers was also strongly suggested by these authors. Beyond the core rbcL+matK barcode, the most widely used plastid barcoding marker is the intergenic spacer trnH-psbA. This region is straightforward to amplify across land plants, and is one of the more variable intergenic spacers in plants [79]. Another suggestion is supplementing the phylogenetic analysis of plants with the internal transcribed spacer from nuclear ribosomal DNA (nrITS) [80,81]. A 20% gain in discriminatory power for full ITS and 10–15% gain even when just the ITS2 region is used is an observed and appreciable benefit [82].
Following these recommendations, we have chosen four most widely used and routinely included plant DNA barcodes out of seven leading and 15 known plant sequences considered as barcodes [78]. These were as follows: two plastid coding regions rbcL and matK, plastid intergenic spacer trnH-psbA, and nuclear internal transcribed spacer ITS2.
We have amplified, sequenced, and uploaded to GeneBank four matK (MT663156, MT663157, MT663158, MT663159) and four rbcL (MT663160, MT663161, MT663162, MT663163) fragments obtained for cultivars: S.OBRL.01, S.OBRL.03, S.OBRL.11, and S.OBRL.12, respectively. However, they were not used for phylogenetic analysis as they possessed little (about 1%) genetic diversity within the Salvia genus. We have also obtained novel DNA sequences for trnH-psbA and ITS2 barcodes for our cultivars: S.OBRL.01, S.OBRL.03, S.OBRL.11 and S.OBRL.12, which ensured significantly higher species discrimination success in our study.
In our approach of an extended marker framework for DNA barcodes, a combination of trnH-psbA and ITS2 barcodes showed sufficient discriminatory power to identify analysed plant organisms to the species level. The phylogram tree constructed with trnH-psbA and ITS2 barcodes has shown the existence of a separate Perovskia clade, which is clearly embedded within Salvia, which corroborates the findings of Drew et al. [14], who formally transferred former non-Salvia genera, such as Perovskia, to Salvia and created a subgeneric designation within Salvia for each of these former genera.
In the present form, our molecular analysis may suggest that S. abrotanoides and S. yangii are very closely related but regarded as two separate species. To make this conclusion fully justified, multiple accessions of S. abrotanoides would have to be included into the analysis. Unfortunately, at present, we do not have access to other taxonomically undoubtful accessions of S. abrotanoides; thus, we decided to rely on what is available. Despite a suggested close genetic relationship, S. abrotanoides and S. yangii possess distinct phytochemical profiles including qualitative differences within the non-volatile chemicals, which make them a potentially useful model for further analysis. Interestingly, the LC-MS profiles of aerial parts are apparently more similar, and none of the major compounds stand out as distinct, but the PCA processing revealed an obvious separation between the species. Hence, presence or absence of the minor compounds or their varied proportions is also important with such examples as the abovementioned sugiol/trilobinol pair. In contrast, the roots of both species are phytochemically distinct at first sight, and chemometric analysis has verified this fact. Isograndifoliol (Figure 2, peak 27), a hydroxylated norditerpenoid present only in S. yangii, is a rarely detected compound that has a high pharmacological potential [66], and its presence in a relatively high amount points out S. yangii as an especially valuable medicinal plant. Sugiol, even if a minor compound only, also has significant pharmacological potential, reviewed recently by Bajpai et al. 2020 [83], that can, for example, attenuate STAT3-dependent proliferation and unlocks apoptosis pathways in cancer cells and can block virus entry to cells. Whether or not this compound can contribute to the overall activity of S. yangii leaf as a drug has not been explored so far. On the other hand, carnosic acid quinone, a phytochemical marker of S. abrotanoides leaves, was reported as a highly potent factor in antioxidant reactions in S. rosmarinus and S. officinalis [84].
Yet, the phytochemical profiling alone is not sufficient to prove the taxonomic differentiation, especially if easily accessible organs (leaves) would be used for sampling. On the other hand, the results from molecular studies are independent from the sample type and coupled with parallel LC-MS profiling constitute a powerful tool in monitoring of taxonomically close and probably easily crossing species.
Qualitative chemical differences between closely related plant species might be of a different etiology. They might result directly from genetic (and/or epigenetic) differences between species or be an outcome of various processes occurring at the transcriptional and post-transcriptional levels. Comparative analysis of S. abrotanoides and S. yangii transcriptomes might help discover particular transcripts responsible for detected chemical differences. Performing genetic crosses between S. abrotanoides and S. yangii would create a segregating population which, through parallel phenotyping and genotyping, might serve for tracing genetic loci potentially responsible for the observed chemical profiles. A similar approach in Papaver somniferum allowed for discovering that ten genes from the noscapine biosynthesis pathway are clustered in the poppy genome [85]. In fact, there are numerous reports published to date proving that many secondary metabolites are a product of gene-clusters [86]. Once a certain target genetic locus is found, it could be subject to different biotechnological approaches, such as gene cloning, gene silencing or heterologous expression of proteins, in order to perform a functional analysis of target genes. In the long run, this may lead to developing applied research aiming at optimization of the chemical profile of Salvia, which may in turn become an alternative source of bioactive chemicals such as tanshinones.

Author Contributions

Conceptualization, M.B.; methodology, M.B.; validation, M.B.; formal analysis, B.P. and S.Ś.; investigation, M.B., B.P., M.S. and S.Ś.; resources, K.J., M.R., S.G. and A.M.; data curation, M.B.; writing—original draft preparation, M.B.; writing—review and editing, A.M.; visualization, M.B.; supervision, A.M.; project administration, M.B.; funding acquisition, M.R., A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Wroclaw Medical University, grant number, SUB.D033.19.009, and by the Ministry of Science and Higher Education grant-in-aid for Special Research Facility, grant number SPUBD.033.19.002.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

DNA sequence data generated in this study has been deposited in GeneBank at accession numbers: MT599312, MT599313, MT599314, MT599315, MT815871, MT815872, MT815873, MT815874. The UHPLC-QTOF-MS raw data presented in this study are available on request from the corresponding author. The data are not publicly available due to the institutional policy of the authors’ employer.

Acknowledgments

The authors thank Agata Biernat for cultivating the plants in the Botanical Garden of Medicinal Plants. We would like to thank Hanna Czapor-Irzabek for technical assistance in UHPLC-QTOF-MS data analysis. Measurements were carried out in the Laboratory of Elemental Analysis and Structural Research at Wroclaw Medical University.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Plants growing at the site of the experiment: Salvia yangii (a), Salvia abrotanoides (b) and mature leaves with fragments of inflorescence of Salvia yangii (c), Salvia abrotanoides (d); all photographs are from the collection of the Botanical Garden of Medicinal Plants, Wroclaw.
Figure A1. Plants growing at the site of the experiment: Salvia yangii (a), Salvia abrotanoides (b) and mature leaves with fragments of inflorescence of Salvia yangii (c), Salvia abrotanoides (d); all photographs are from the collection of the Botanical Garden of Medicinal Plants, Wroclaw.
Cells 10 00112 g0a1

References

  1. Flora of Pakistan. Available online: http://www.efloras.org/flora_page.aspx?flora_id=5 (accessed on 17 January 2020).
  2. Perveen, S.; Khan, S.B.; Malik, A.; Tareen, R.B.; Nawaz, S.A.; Choudhary, M.I. Phenolic constituents from Perovskia atriplicifolia. Nat. Prod. Res. 2006, 20, 347–353. [Google Scholar] [CrossRef]
  3. Perveen, S.; Malik, A.; Noor, A.T.; Tareen, R.B. Pervosides A and B, new isoferulyl glucosides from Perovskia atriplicifolia. J. Asian Nat. Prod. Res. 2008, 10, 1105–1108. [Google Scholar] [CrossRef]
  4. Ahmad, I.; Fozia; Waheed, A.; Tahir, N.B.; Rais, A.K. Anti-inflammatory constituents from Perovskia atriplicifolia. Pharm. Biol. 2015, 53, 1628–1631. [Google Scholar] [CrossRef]
  5. Rechinger, K.H. Flora Iranica; Akademische Druck-U. Verlagsanstalt: Graz, Austria, 1982. [Google Scholar]
  6. Mozaffarian, V. A Dictionary of Iranian Plant Names; Farhang Moaser: Tehran, Iran, 1996. [Google Scholar]
  7. Khaliq, S.; Volk, F.J.; Frahm, A.W. Phytochemical investigation of Perovskia abrotanoides. Planta Med. 2007, 73, 77–83. [Google Scholar] [CrossRef]
  8. Flora of China. Available online: http://www.efloras.org/flora_page.aspx?flora_id=2 (accessed on 17 January 2020).
  9. Gao, L.; Zhou, J.; Zhu, L.Y.; Zhang, J.R.; Jing, Y.X.; Zhao, J.W.; Huang, X.Z.; Li, G.P.; Jiang, Z.Y.; Xue, D.Y. Four New Diterpene Glucosides from Perovskia atriplicifolia. Chem. Biodivers. 2017, 14, 1–7. [Google Scholar] [CrossRef]
  10. Harley, R.M.; Atkins, S.; Budantsev, A.L.; Cantino, P.D.; Conn, B.J.; Grayer, R.; Harley, M.M.; De Kok, R.; Krestovskaja, T.; Morales, R.; et al. Labiatae. In Flowering Plants · Dicotyledons. The Families and Genera of Vascular Plants; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  11. Walker, J.B.; Sytsma, K.J.; Treutlein, J.; Wink, M. Salvia (Lamiaceae) is not monophyletic: Implications for the systematics, radiation, and ecological specializations of Salvia and tribe Mentheae. Am. J. Bot. 2004, 91, 1115–1125. [Google Scholar] [CrossRef]
  12. Walker, J.B.; Sytsma, K.J. Staminal evolution in the genus Salvia (Lamiaceae): Molecular phylogenetic evidence for multiple origins of the staminal lever. Ann. Bot. 2007, 100, 375–391. [Google Scholar] [CrossRef]
  13. Will, M.; Claßen-Bockhoff, R. Time to split Salvia s.l. (Lamiaceae)—New insights from Old World Salvia phylogeny. Mol. Phylogenet. Evol. 2017, 109, 33–58. [Google Scholar] [CrossRef]
  14. Drew, B.T.; González-Gallegos, J.G.; Xiang, C.L.; Kriebel, R.; Drummond, C.P.; Walker, J.B.; Sytsma, K.J. Salvia united: The greatest good for the greatest number. Taxon 2017, 66, 133–145. [Google Scholar] [CrossRef] [Green Version]
  15. Perveen, S.; Malik, A.; Tareen, R.B. Structural determination of atricins A and B, new triterpenes from Perovskia atriplicifolia, by 1D and 2D NMR spectroscopy. Magn. Reson. Chem. 2009, 47, 266–269. [Google Scholar] [CrossRef]
  16. Grubov, V.I. Plants of Central Asia: Plant Collection from China and Mongolia; Science Publishers, Inc.: Enfield, NH, USA; Plymouth, UK, 2002. [Google Scholar]
  17. Eisenman, S.W.; Zaurov, D.E.; Struwe, L. Medicinal Plants of Central Asia: Uzbekistan and Kyrgyzstan; Springer: Berlin/Heidelberg, Germany, 2013; ISBN 9781461439127. [Google Scholar]
  18. Grant, M. RHS Plant Trials and Assessments, Perovskia; Supplementary to RHS Trials Award; Royal Horticultural Society: London, UK, 2007; pp. 1–6. [Google Scholar]
  19. Huxley, A. The New Royal Horticultural Society Dictionary of Gardening, 1st ed.; Nature Pub Group: London, UK, 1999; ISBN 1561592404. [Google Scholar]
  20. Plants of the World Online. Available online: http://www.plantsoftheworldonline.org/ (accessed on 26 January 2020).
  21. Tareen, R.B.; Bibi, T.; Khan, M.A.; Ahmad, M.; Zafar, M. Indigenous knowledge of folk medicine by the women of Kalat and Khuzdar regions of Balochistan, Pakistan. Pak. J. Bot. 2010, 42, 1465–1485. [Google Scholar]
  22. Mohammadhosseini, M.; Venditti, A.; Akbarzadeh, A. The genus Perovskia Kar.: Ethnobotany, chemotaxonomy and phytochemistry: A review. Toxin Rev. 2019, 1–22. [Google Scholar] [CrossRef]
  23. Baquar, S.R. Medicinal and Poisonous Plants of Pakistan, 1st ed.; Printas Press: Karachi, Pakistan, 1989. [Google Scholar]
  24. Jiang, Z.Y.; Yu, Y.J.; Huang, C.G.; Huang, X.Z.; Hu, Q.F.; Yang, G.Y.; Wang, H.B.; Zhang, X.Y.; Li, G.P. Icetexane diterpenoids from Perovskia atriplicifolia. Planta Med. 2015, 81, 241–246. [Google Scholar] [CrossRef] [Green Version]
  25. Sairafianpour, M.; Christensen, J.; Stærk, D.; Budnik, B.A.; Kharazmi, A.; Bagherzadeh, K.; Jaroszewski, J.W. Leishmanicidal, antiplasmodial, and cytotoxic activity of novel diterpenoid 1,2-quinones from Perovskia abrotanoides: New source of tanshinones. J. Nat. Prod. 2001, 64, 1398–1403. [Google Scholar] [CrossRef]
  26. Jaafari, M.R.; Hooshmand, S.; Samiei, A.; Hossainzadeh, H. Evaluation of leishmanicidal effect of Perovskia abrotanoides Karel. Root extract by in vitro leishmanicidal assay using promastigotes of Leishmania major. Pharmacologyonline 2007, 1, 299–303. [Google Scholar]
  27. Kumar, G.P.; Gupta, S.; Murugan, M.P.; Bala Singh, S. Ethnobotanical studies of Nubra valley—A cold arid zone of Himalaya. Ethnobot. Leafl. 2009, 13, 752–765. [Google Scholar]
  28. Moallem, S.A.; Niapour, M. Study of embryotoxicity of Perovskia abrotanoides, an adulterant in folk-medicine, during organogenesis in mice. J. Ethnopharmacol. 2008, 117, 108–114. [Google Scholar] [CrossRef]
  29. Ballabh, B.; Chaurasia, O.P.; Ahmed, Z.; Singh, S.B. Traditional medicinal plants of cold desert Ladakh-Used against kidney and urinary disorders. J. Ethnopharmacol. 2008, 118, 331–339. [Google Scholar] [CrossRef]
  30. Hosseinzadeh, H.; Amel, S. Antinociceptive Effect of the Aerial Parts of Perovskia abrotanoides Extracts in Mice. Med. J. Iran. Hosp. 2001, 4, 15–17. [Google Scholar]
  31. Mahboubi, M.; Kazempour, N. The antimicrobial activity of essential oil from Perovskia abrotanoides karel and its main components. Indian J. Pharm. Sci. 2009, 71, 343–347. [Google Scholar] [CrossRef] [Green Version]
  32. Sardashti, A.R.; Valizadeh, J.; Adhami, Y. Variation in the essential oil composition of Perovskia abrotonoides of different growth stage in Baluchestan. Middle East J. Sci. Res. 2013, 13, 781–784. [Google Scholar] [CrossRef]
  33. Sajjadi, S.E.; Mehregan, I.; Khatamsaz, M.; Asgari, G. Chemical composition of the essential oil of Perovskia abrotanoides Karel. growing wild in Iran. Flavour Fragr. J. 2005, 20, 445–446. [Google Scholar] [CrossRef]
  34. Oreizi, E.; Rahiminejad, M.R.; Asghari, G. Influence of environment on glandular trichomes and composition of essential oil of Perovskia abrotanoides Karel. Jundishapur J. Nat. Pharm. Prod. 2014, 9, e16432. [Google Scholar] [CrossRef] [Green Version]
  35. Pourmortazavi, S.M.; Sefidkon, F.; Hosseini, S.G. Supercritical carbon dioxide extraction of essential oils from Perovskia atriplicifolia Benth. J. Agric. Food Chem. 2003, 51, 5414–5419. [Google Scholar] [CrossRef]
  36. Morteza-Semnani, K. The essential oil composition of Perovskia abrotanoides from Iran. Pharm. Biol. 2004, 42, 214–216. [Google Scholar] [CrossRef]
  37. Jassbi, A.R.; Ahmad, V.U.; Tareen, R.B. Constituents of the essential oil of Perovskia atriplicifolia Benth. Flavour Fragr. J. 1999, 14, 38–40. [Google Scholar] [CrossRef]
  38. Inouye, S.; Uchida, K.; Yamaguchi, H.; Miyara, T.; Gomi, S.; Amano, M. Volatile aroma constituents of three labiatae herbs growing wild in the karakoram-himalaya district and their antifungal activity by vapor contact. J. Essent. Oil Res. 2001, 13, 68–72. [Google Scholar] [CrossRef]
  39. Erdemgil, F.Z.; Ilhan, S.; Korkmaz, F.; Kaplan, C.; Mercangöz, A.; Arfan, M.; Ahmad, S. Chemical composition and biological activity of the essential oil of Perovskia atriplicifolia from Pakistan. Pharm. Biol. 2007, 45, 324–331. [Google Scholar] [CrossRef]
  40. Dabiri, M.; Sefidkon, F. Analysis of the essential oil from aerial parts of Perovskia atriplicifolia Benth. at different stages of plant growth. Flavour Fragr. J. 2001, 16, 435–438. [Google Scholar] [CrossRef]
  41. Ashraf, S.N.; Zubair, M.; Rizwan, K.; Tareen, R.B.; Rasool, N.; Zia-Ul-Haq, M.; Ercisli, S. Compositional studies and biological activities of Perovskia abrotanoides kar. oils. Biol. Res. 2014, 47, 1–9. [Google Scholar] [CrossRef] [Green Version]
  42. Senol, F.S.; Ślusarczyk, S.; Matkowski, A.; Pérez-Garrido, A.; Girón-Rodríguez, F.; Cerón-Carrasco, J.P.; Den-Haan, H.; Peña-García, J.; Pérez-Sánchez, H.; Domaradzki, K.; et al. Selective in vitro and in silico butyrylcholinesterase inhibitory activity of diterpenes and rosmarinic acid isolated from Perovskia atriplicifolia Benth. and Salvia glutinosa L. Phytochemistry 2017, 133, 33–44. [Google Scholar] [CrossRef] [PubMed]
  43. Petersen, M.; Simmonds, M.S.J. Rosmarinic acid. Phytochemistry 2003, 62, 121–125. [Google Scholar] [CrossRef]
  44. Slusarczyk, S.; Topolski, J.; Domaradzki, K.; Adams, M.; Hamburger, M.; Matkowski, A. Isolation and fast selective determination of nor-abietanoid diterpenoids from Perovskia atriplicifolia roots using LC-ESI-MS/MS with multiple reaction monitoring. Nat. Prod. Commun. 2015, 10, 1149–1152. [Google Scholar] [CrossRef] [PubMed]
  45. Ge, X.; Wu, J. Tanshinone production and isoprenoid pathways in Salvia miltiorrhiza hairy roots induced by Ag + and yeast elicitor. Plant Sci. 2005, 168, 487–491. [Google Scholar] [CrossRef]
  46. Bulgakov, V.P.; Inyushkina, Y.V.; Fedoreyev, S.A. Rosmarinic acid and its derivatives: Biotechnology and applications. Crit. Rev. Biotechnol. 2012, 32, 203–217. [Google Scholar] [CrossRef]
  47. Lu, S. Compendium of Plant Genomes; Springer Nature Switzerland AG: Cham, Switzerland, 2019; ISBN 978-3-030-24715-7. [Google Scholar]
  48. Wang, J.W.; Wu, J.Y. Tanshinone biosynthesis in Salvia miltiorrhiza and production in plant tissue cultures. Appl. Microbiol. Biotechnol. 2010, 88, 437–449. [Google Scholar] [CrossRef]
  49. Aoyagi, Y.; Takahashi, Y.; Satake, Y.; Takeya, K.; Aiyama, R.; Matsuzaki, T.; Hashimoto, S.; Kurihara, T. Cytotoxicity of abietane diterpenoids from Perovskia abrotanoides and of their semisynthetic analogues. Bioorganic Med. Chem. 2006, 14, 5285–5291. [Google Scholar] [CrossRef]
  50. Karimzadeh, S.M.; Moridi Farimani, M.; Amiri, M.S.; Tabefam, M.; Alilou, M.; Stuppner, H. Perovskanol, a new sesquiterpenoid with an unprecedented skeleton from Perovskia abrotanoides. Nat. Prod. Res. 2019, 1–5. [Google Scholar] [CrossRef]
  51. Tarawneh, A.; León, F.; Pettaway, S.; Elokely, K.M.; Klein, M.L.; Lambert, J.; Mansoor, A.; Cutler, S.J. Flavonoids from Perovskia atriplicifolia and Their in Vitro Displacement of the Respective Radioligands for Human Opioid and Cannabinoid Receptors. J. Nat. Prod. 2015, 78, 1461–1465. [Google Scholar] [CrossRef]
  52. Liu, W.X.; Zhao, J.W.; Zuo, A.X.; Yang, Z.; Gao, L.; Zhou, M.; Jiang, Z.Y. Two novel terpenoids from the cultured Perovskia atriplicifolia. Fitoterapia 2018, 130, 152–155. [Google Scholar] [CrossRef]
  53. Borges, A.; Rosa, M.; Recchia, G.; Queiroz-Silva, J.; Bressan, E.; Veasey, E. CTAB methods for DNA extraction of sweetpotato for microsatellite analysis. Sci. Agric. 2009, 66, 529–534. [Google Scholar] [CrossRef]
  54. Gardes, M.; Bruns, T.D. ITS primers with enhanced specificity for basidiomycetes—Application to the identification of mycorrhizae and rusts. Mol. Ecol. 1993, 2, 113–118. [Google Scholar] [CrossRef] [PubMed]
  55. White, T.J.; Bruns, S.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protoc. Guide Methods Appl. 1990, 18, 315–322. [Google Scholar]
  56. Tate, J.A.; Simpson, B.B. Paraphyly of Tarasa (Malvaceae) and Diverse Origins of the Polyploid Species. Syst. Bot. 2003, 28, 723–737. [Google Scholar] [CrossRef]
  57. Sang, T.; Crawford, D.J.; Stuessy, T.F. Chloroplast DNA phylogeny, reticulate evolution, and biogeography of Paeonia (Paeoniaceae). Am. J. Bot. 1997, 84, 1120–1136. [Google Scholar] [CrossRef] [Green Version]
  58. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  59. Bengtsson-Palme, J.; Ryberg, M.; Hartmann, M.; Branco, S.; Wang, Z.; Godhe, A.; De Wit, P.; Sánchez-García, M.; Ebersberger, I.; De Sousa, F.; et al. Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol. Evol. 2013, 4, 914–919. [Google Scholar] [CrossRef]
  60. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  61. Dereeper, A.; Guignon, V.; Blanc, G.; Audic, S.; Buffet, S.; Chevenet, F.; Dufayard, J.F.; Guindon, S.; Lefort, V.; Lescot, M.; et al. Phylogeny.fr: Robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36, 465–469. [Google Scholar] [CrossRef]
  62. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  63. Tamura, K.; Nei, M.; Kumar, S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. USA 2004, 101, 11030–11035. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Logacheva, M.D.; Valiejo-Roman, C.M.; Pimenov, M.G. Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Plant Syst. Evol. 2008, 270, 783–791. [Google Scholar] [CrossRef]
  65. Yan, X. Dan Shen (Salvia miltiorrhiza) in Medicine; Volume 2. Pharmacology and Quality Control; Springer: Dordrecht, The Netherlands, 2015; ISBN 978-94-017-9462-6. [Google Scholar]
  66. Slusarczyk, S.; Sezer Senol Deniz, F.; Abel, R.; Pecio, Ł.; Pérez-Sánchez, H.; Cerón-Carrasco, J.P.; Den-Haan, H.; Banerjee, P.; Preissner, R.; Krzyzak, E.; et al. Norditerpenoids with selective anti-cholinesterase activity from the roots of Perovskia atriplicifolia benth. Int. J. Mol. Sci. 2020, 21, 4475. [Google Scholar] [CrossRef] [PubMed]
  67. Ghaffari, Z.; Rahimmalek, M.; Sabzalian, M.R. Variations in Essential Oil Composition and Antioxidant Activity in Perovskia abrotanoides Kar. Collected from Different Regions in Iran. Chem. Biodivers. 2018, 15. [Google Scholar] [CrossRef] [PubMed]
  68. Ghaffari, Z.; Rahimmalek, M.; Sabzalian, M.R. Variation in the primary and secondary metabolites derived from the isoprenoid pathway in the Perovskia species in response to different wavelengths generated by light emitting diodes (LEDs). Ind. Crops Prod. 2019, 140, 111592. [Google Scholar] [CrossRef]
  69. Tabefam, M.; Farimani, M.M.; Danton, O.; Ramseyer, J.; Kaiser, M.; Ebrahimi, S.N.; Salehi, P.; Batooli, H.; Potterat, O.; Hamburger, M. Antiprotozoal Diterpenes from Perovskia abrotanoides. Planta Med. 2018, 84, 913–919. [Google Scholar] [CrossRef]
  70. Hu, G.X.; Takano, A.; Drew, B.T.; Liu, E.-D.; Soltis, D.E.; Soltis, P.S.; Peng, H.; Xiang, C.L. Phylogeny and staminal evolution of Salvia (Lamiaceae, Nepetoideae) in East Asia. Ann. Bot. 2018, 122, 649–668. [Google Scholar] [CrossRef]
  71. Pourhosseini, S.H.; Hadian, J.; Sonboli, A.; Nejad Ebrahimi, S.; Mirjalili, M.H. Genetic and Chemical Diversity in Perovskia abrotanoides Kar. (Lamiaceae) Populations Based on ISSRs Markers and Essential Oils Profile. Chem. Biodivers. 2018, 15. [Google Scholar] [CrossRef]
  72. Hashemifar, Z.; Rahimmalek, M. Genetic structure and variation in Perovskia abrotanoides Karel and P. atriplicifolia as revealed by molecular and morphological markers. Sci. Hortic. 2018, 230, 169–177. [Google Scholar] [CrossRef]
  73. Sucher, N.J.; Hennell, J.R.; Carles, M.C. Plant DNA Fingerprinting and Barcoding; Springer: Berlin/Heidelberg, Germany, 2013; ISBN 9781627033329. [Google Scholar]
  74. Chase, M.W.; Cowan, R.S.; Hollingsworth, P.M.; Van Den Berg, C.; Madriñán, S.; Petersen, G.; Seberg, O.; Jørgsensen, T.; Cameron, K.M.; Pedersen, N.; et al. A Proposal for a Standardised Protocol to Barcode All Land Plants. Taxon 2007, 56, 295–299. [Google Scholar] [CrossRef]
  75. CBOL Plant Working Group. A DNA barcode for land plants. Proc. Natl. Acad. Sci. USA 2009, 106, 12794–12797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Hebert, P.D.N.; Cywinska, A.; Ball, S.L.; DeWaard, J.R. Biological identifications through DNA barcodes. Proc. R. Soc. B Biol. Sci. 2003, 270, 313–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Fazekas, A.J.; Burgess, K.S.; Kesanakurti, P.R.; Graham, S.W.; Newmaster, S.G.; Husband, B.C.; Percy, D.M.; Hajibabaei, M.; Barrett, S.C.H. Multiple multilocus DNA barcodes from the plastid genome discriminate plant species equally well. PLoS ONE 2008, 3, e2802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Hollingsworth, P.M.; Graham, S.W.; Little, D.P. Choosing and using a plant DNA barcode. PLoS ONE 2011, 6, e19254. [Google Scholar] [CrossRef] [PubMed]
  79. Shaw, J.; Lickey, E.B.; Schilling, E.E.; Small, R.L. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The Tortoise and the hare III. Am. J. Bot. 2007, 94, 275–288. [Google Scholar] [CrossRef] [Green Version]
  80. Li, D.Z.; Gao, L.M.; Li, H.T.; Wang, H.; Ge, X.J.; Liu, J.Q.; Chen, Z.D.; Zhou, S.L.; Chen, S.L.; Yang, J.B.; et al. Comparative analysis of a large dataset indicates that internal transcribed spacer (ITS) should be incorporated into the core barcode for seed plants. Proc. Natl. Acad. Sci. USA 2011, 108, 19641–19646. [Google Scholar] [CrossRef] [Green Version]
  81. Chen, S.; Yao, H.; Han, J.; Liu, C.; Song, J.; Shi, L.; Zhu, Y.; Ma, X.; Gao, T.; Pang, X.; et al. Validation of the ITS2 region as a novel DNA barcode for identifying medicinal plant species. PLoS ONE 2010, 5, e8613. [Google Scholar] [CrossRef]
  82. Hollingsworth, P.M. Refining the DNA barcode for land plants. Proc. Natl. Acad. Sci. USA 2011, 108, 19451–19452. [Google Scholar] [CrossRef] [Green Version]
  83. Bajpai, V.K.; Sonwal, S.; Hwang, S.-K.; Shukla, S.; Khan, I.; Dey, D.K.; Chen, L.; Simal-Gandara, J.; Xiao, J.; Huh, Y.S.; et al. Sugiol, a diterpenoid: Therapeutic actions and molecular pathways involved. Pharmacol. Res. 2020, 105313. [Google Scholar] [CrossRef]
  84. Masuda, T.; Inaba, Y.; Maekawa, T.; Takeda, Y.; Tamura, H.; Yamaguchi, H. Recovery mechanism of the antioxidant activity from carnosic acid quinone, an oxidized sage and rosemary antioxidant. J. Agric. Food Chem. 2002, 50, 5863–5869. [Google Scholar] [CrossRef]
  85. Winzer, T.; Gazda, V.; He, Z.; Kaminski, F.; Kern, M.; Larson, T.R.; Li, Y.; Meade, F.; Teodor, R.; Vaistij, F.E.; et al. A Papaver somniferum 10-gene cluster for synthesis of the anticancer alkaloid noscapine. SOM. Science 2012, 336, 1704–1708. [Google Scholar] [CrossRef] [Green Version]
  86. Nützmann, H.W.; Osbourn, A. Gene clustering in plant specialized metabolism. Curr. Opin. Biotechnol. 2014, 26, 91–99. [Google Scholar] [CrossRef] [Green Version]
Figure 1. BPC chromatogram of the main metabolites identified in the MeOH extract of leaves of S. yangii and S. abrotanoides, harvested in the middle of the season.
Figure 1. BPC chromatogram of the main metabolites identified in the MeOH extract of leaves of S. yangii and S. abrotanoides, harvested in the middle of the season.
Cells 10 00112 g001
Figure 2. BPC chromatogram of the main metabolites identified in the MeOH extract of roots of S. yangii and S. abrotanoides, harvested in the middle of the season.
Figure 2. BPC chromatogram of the main metabolites identified in the MeOH extract of roots of S. yangii and S. abrotanoides, harvested in the middle of the season.
Cells 10 00112 g002
Figure 3. Principal component analysis presenting the relationship between the chemical profile of the main compounds in the leaf and species annotation of the analyzed plants; S. yangii (blue) and S. abrotanoides (orange) (a), with corresponding loadings (PCA model: R2XCum = 0.767, Q2Cum = 0.45) (b). Each loading number corresponds to the peak number identified in the BPC chromatogram.
Figure 3. Principal component analysis presenting the relationship between the chemical profile of the main compounds in the leaf and species annotation of the analyzed plants; S. yangii (blue) and S. abrotanoides (orange) (a), with corresponding loadings (PCA model: R2XCum = 0.767, Q2Cum = 0.45) (b). Each loading number corresponds to the peak number identified in the BPC chromatogram.
Cells 10 00112 g003
Figure 4. Principal component analysis presenting the relationship between the chemical profile of the main compounds in the root and species annotation of the analyzed plants; S. yangii (blue) and S. abrotanoides (orange) (a), with corresponding loadings (PCA model: R2XCum = 0.755, Q2Cum = 0.566) (b). Each loading number corresponds to the peak number identified in the BPC chromatogram.
Figure 4. Principal component analysis presenting the relationship between the chemical profile of the main compounds in the root and species annotation of the analyzed plants; S. yangii (blue) and S. abrotanoides (orange) (a), with corresponding loadings (PCA model: R2XCum = 0.755, Q2Cum = 0.566) (b). Each loading number corresponds to the peak number identified in the BPC chromatogram.
Cells 10 00112 g004
Figure 5. Genetic relationship analysis by the Neighbor-Joining method based on the Maximum Composite Likelihood model and Bootstrap phylogeny test (1000 replicates). Support values are shown above nodes (values lower than 50 are not presented). The branch length is shown below branches. The analysis involved 29 nucleotide sequence combinations. Evolutionary analyses were conducted in MEGA7.
Figure 5. Genetic relationship analysis by the Neighbor-Joining method based on the Maximum Composite Likelihood model and Bootstrap phylogeny test (1000 replicates). Support values are shown above nodes (values lower than 50 are not presented). The branch length is shown below branches. The analysis involved 29 nucleotide sequence combinations. Evolutionary analyses were conducted in MEGA7.
Cells 10 00112 g005
Table 1. Primers used for amplification and Sanger sequencing of the ITS rDNA region and trnH-psbA intergenic spacer.
Table 1. Primers used for amplification and Sanger sequencing of the ITS rDNA region and trnH-psbA intergenic spacer.
PrimerSequence (5′ to 3′)Reference
ITS1FCTTGGTCATTTAGAGGAAGTAA[54]
ITS4TCCTCCGCTTATTGATATGC[55]
psbA3_fGTTATGCATGAACGTAATGCTC[57]
trnHf_05CGCGCATGGTGGATTCACAATCC[56]
Table 2. Characterization of identified compounds in leaves of S. yangii and S. abrotanoides, UHPLC–QTOF-MS in positive ion mode.
Table 2. Characterization of identified compounds in leaves of S. yangii and S. abrotanoides, UHPLC–QTOF-MS in positive ion mode.
No.CompoundRTUVm/z[M + H]+FormulaMS2 Main-Ion (Relative Intensity %)MS2 Fragments (Relative Intensity %)Salvia yangiiSalvia abrotanoidesReferences
1ND0.9-365.1123C13H16O10169(100)271(45), 211(23), 151(18)******
2ND1.7237, 285, 323369.1184C17H20O9207(100) ******
3Quercetin 3-rutinoside2.7282, 342611.1609C27H30O16303(100)465(4)******HMDB0003249
4Luteolin-O-Hex-Dhex4.6260, 343595.1660C27H30O15287(100)449(3)*****MetFrag
5Luteolin-O-Pen-Glucuronide5.1260, 343595.1660C27H30O15287(100)395(5), 419(3)******MetFrag
6Quercetin 4′-glucoside5.3260, 348465.1025C21H20O12303(100) ****HMDB0037932
7ND6.1279, 324603.2056C31H38O12441(100) ******
8Apigenin-O-Hex-Dhex6.7282, 323579.1706C27H31O14271(100)433(3)*****MetFrag
9Hesperidin6.9284611.1970C27H30O16303(100)263(5), 449(2), 369(2)****HMDB0003265
10Diosmetin-O-Hex-Dhex7.2285, 324609.1815C28H32O15301(100)463(3)******HMDB0029548
11Rosmarinic acid7.4289, 328361.0927C18H16O8163(100)181(42), 325(15), 283(10)*****Standard
12ND8.6219, 279503.1915C26H30O10279(100)341(32), 311(14), 291(25)*****
13ND8.8219, 283, 323435.1267C21H22O10391(100)349(15)*****
14Isomargaritene9.2-593.1870C28H32O14285(100)447(12), 429(10)****HMDB0037415
15Poncirin (Isosakuranetin-7-O-beta-D-neohesperidoside)9.2-595.2017C28H34O14287(100)269(2)****HMDB0037487
16ND9.3-503.1920C26H30O10279(100)311(90), 291(80), 251(70), 259(20)****
17Dihydrokaempferide 7-rhamnoside10.5282, 320449.1428C22H24O10285(100) *****HMDB0040561
18Trihydroxy-dimethoxyflavone11.4334331.0813C17H14O7270(100)298(63), 286(4)******HMDB0128752
19ND11.9-631.2018C31H34O14291(100)323(65), 259(8)******
20Isorosmanol derivative12.6-347.1849C20H26O5301(100)273(15), 241(8), 231(7)****MetFrag
21Dihydroxy-methoxymethylflavone13.1275, 333315.0861C17H14O6254(100)282(61), 300(13), 271(7), 285(2)****MetFrag
22Rosmanol13.4288347.1858C20H26O5301(100)273(48), 241(27), 283(20)****MetFrag
23Epi-rosmanol14-347.1856C20H26O6229(100)273(42), 283(25), 231(18),****HMDB0035812
242-Hydroxy-3-oxo-12-oleanen-28-oic acid15.4-471.3471C30H46O4453(100)425(42), 261(18), 201(12)*** HMDB0040494
253-Oxo-12,18-ursadien-28-oic acid16-453.3377C30H44O3407(100)201(42), 435(10), 205(10)*** HMDB0037065
26Trihydroxy-12-ursen-28-oic acid16.6-489.3579C30H48O5407(100)471(25), 435(14)*** HMDB0036961
27Hydroxy-trimethoxyflavone16.8331329.1016C18H16O6296(100)329(15), 284(4)*****HMDB0040719
2811-Methylrosmanol17-361.2020C21H28O4329(100)301(48), 283(15)****MetFrag
297-Methylrosmanol17.5-361.2045C21H28O5273(100)245(90) 301(67), 237(50), 283(41)*****HMDB0035813
30Carnosol17.8284331.1903C20H26O4285(100)215(45), 267(42), 243(21), 331(20)******HMDB0002121
31Carnosic acid quinone18.2284331.1949C20H26O5285(100)303(42), 267(25), 243(18), 225(11), 215(10), 192(8) ***MetFrag
3212-Hydroxy-7-oxo-8,11,13-abietatrien-18-al18.9 315.1938C20H26O3287(100) ****HMDB0040746
33Isorosmanol19.5-347.1858C20H26O5301(100)231(91), 259(85), 283(45), 255(21) ***HMDB0036661
3411,12-Dimethylrosmanol19.9-375.2166C22H30O5287(100)329(45), 269(22), 217(15), 191(8)******HMDB0040525
35Rosmaridiphenol20.1-317.2104C20H28O3181(100)299(15), 273(14), 261(8)*** HMDB0037233
36Methoxy-8,11,13-abietatrien-20,11-olide20.3-329.2112C21H28O3217(100)191(42), 205(33), 231(17), 245(15)******HMDB0038391
37Sageone20.4-301.1789C19H24O3259(100)241(27), 272(20), 283(18), 303(13), (216(4)****HMDB0038684
38Sugiol21.1205, 223, 282301.2153C20H28O2272(100)191(60), 217(51), 203(42),*** HMDB0036564
39Trilobinol21.2205, 223, 282301.2153C20H28O2272(100)191(85), 217(31), 203(22), ***HMDB0038702
40Carnosic acid21.6-333.2056C20H28O4315(100)287(65), 273(41), 259(12), 245(5), 233(2)****Standard
41Miltirone22.2-283.1678C19H22O2225(100)268(40), 240(28), 253(20)*****[44]
ND: not determined, * ranges from 0.01–0.1 mg/g DM (Dry Mass), ** 0.1–1 mg/g DM and *** above 1 mg/g DM HMDB ID: The Human Metabolome Database MetFrag-MetFrag Online DataBase.
Table 3. Characterization of identified compounds in roots of S. yangii and S. abrotanoides, UHPLC–QTOF-MS in positive ion mode.
Table 3. Characterization of identified compounds in roots of S. yangii and S. abrotanoides, UHPLC–QTOF-MS in positive ion mode.
No.CompoundRTUVm/z[M + H]+FormulaMS2 Main-Ion (Relative Intensity %)MS2 Fragments (Relative Intensity %)Salvia yangiiSalvia abrotanoidesReferences
1Methoxytaxifolin1.3283, 323335.0763C16H14O8271(100)253(23), 225(14), 197(8)* MetFrag
2ND1.8217, 241, 292, 323549.1594C26H28O13531(100)313(8), 251(4), 295(3)*
32-O-p-Coumaroyltartronic acid2.3-267.0483C12H10O7251(100)223(23), 206(15)* MetFrag
4ND5.4287, 319511.2535 349(100)331(47), 258(41), 245(25), 227(17), 211(10)******
5Salviaflaside5.6288, 319523.1439C24H26O13163(100)325(45), 287(22), 361(10)******[65]
6Salvianolic acid L6.5256, 285, 315, 350719.1599C36H30O16521(100)295(25), 493(18), 221(12), 249(10)******[65]
7Hesperidin6.9285611.1982C28H34O15303(100)449(10), 263(4) ***HMDB03265
8Chrysoeriol 7-rutinoside7.3333609.1812C28H32O16301(100)463(10) ***HMDB37453
9Rosmarinic acid7.5288, 319361.0924C18H16O8163(100)181(42), 325(15), 283(10)******[44]
10Melitric acid B7.9285, 329521.1062C27H20O11295(100)493(22), 457(14), 329(6)*** HMDB0040680
113-(3,4-dimethoxyphenyl)-5-hydroxy-7-methoxy-8-methyl-3,4-dihydro-2H-1-benzopyran-4-one9.8-345.1334C19H20O6327(100)283(45), 268(15), 255(8), 201(4)****HMDB0129573
123-(3,4-dimethoxyphenyl)-5-hydroxy-7-methoxy-8-methylchromen-4-one11-343.1170C19H18O6325(100)281(51), 266(18), 251(11), 238(10), 211(4)****HMDB0129572
13ND11.6218, 339315.0859C17H14O6267(100)225(15)***
14Isograndifoliol derivative11.8-303.1940C19H26O3272(100)203(18), 256(14), 215(7)*****[66]
15Tanshinone V12.1-315.1210C18H18O5297(100)264(84), 277(21), 236(18)*****[42]
16ND12.1 295.1317 267(100)252(41), 237(28), 225(12)****
17Dihydroxycryptotanshinon12.3-329.1381C19H20O5265(100)283(65), 250(18), 311(14), 237(8)****[44]
18OH-Tanshindiol A12.4-331.1540C19H22O5295(100)277(42), 267(41), 249(32), 235(21), 225(10)****[65]
19OH-cryptotanshinone13.5-313.1435C19H20O4267(100)295(81), 280(74), 249(41), 225(25), 197(20)****[44]
20Oxocryptotanshinone13.7-311.1279C19H18O4265(100)283(74), 250(51), 295(14), 237(10), 197(8)****[44]
211-βOH-cryptotanshinone13.8220, 267313.1427C19H20O4267(100)295(81), 280(74), 249(41), 262(25), 197(20), 225(8)****[44]
22ND14.2 329.1388C19H20O5311(100)285(24), 267(15), 252(8), 185(2)****
23ND14.4-331.1898C19H20O4311(100)267(14), 283(10), 225(8), 185(5)****
24ND15.5-327.1228C19H18O5309(100)265(15), 250(14), 223(11), 281(8)***
25Didehydrotanshinone IIa16.1222, 273293.1163C19H17O3278(100)263(45), 247(41), 232(38), 219(18), 204(5)*** [44]
26Acetyloxy cryptotanshinone16.3-355.1548C21H22O5295(100)280(61), 277(48), 267(38), 262(32), 249(15), 247(8)*** [66]
27Isograndifoliol16.7-303.1956C19H26O3272(100)203(41), 215(25), 257(21), 187(18), 229(15)*** [66]
28Grandifoliol16.8-289.1793C18H24O3271(100)256(45), 241(18)***
29Cryptotanshinone derivative16.9-297.1484C19H20O3269(100)254(21), 239(25), 227(18), 209(3) ***[42]
30Didehydro acetyloxy cryptotanshinone18.6-353.1388C21H20O5315(100)293(85), 275(45), 247(28)*** [65]
31Cryptotanshinone18.9222, 264297.1469C19H20O3269(100)254(41), 251(38), 279(21), 282(18)****[44]
32Ketoisograndifoliol19.4-301.1796C19H24O3272(100)203(45), 215(32), 229(18), 257(12)*** MetFrag
33ND20.4-287.1640C18H23O3269(100)254(45), 213(22), 199(18)****
34Tanshinone IIa21.7-295.1326C19H19O3277(100)249(71), 262(56), 252(41), 266(20)*****[44]
35Miltiodiol (Arucadiol)22-299.1636C19H22O3271(100)203(41), 253(24), 183(7)*****[44]
36ND22.1-269.1520C18H20O2254(100)239(14), 225(12), 215(6)***
37Miltirone22.2-283.1689C19H22O2225(100)268(41), 240(32),253(28)*****[44]
386,7-Didehydroferruginol22.9-285.2197C20H28O271(100)201(12), 272(8)*****
ND: not determined, * ranges from 0.01–0.1 mg/g DM (Dry Mass), ** 0.1–1 mg/g DM and *** above 1 mg/g DM HMDB ID: The Human Metabolome Database MetFrag-MetFrag Online DataBase.
Table 4. Results of BLAST analysis.
Table 4. Results of BLAST analysis.
Species (Voucher)DNA RegionAccession Number (Query)BLAST ResultAccession Number (Result)Identity
Salvia yangii “Lacey blue” (S.OBRL.01)ITS2MT599312Perovskia atriplicifolia
(Salvia yangii)
KJ584242.199.76%
Salvia yangii “Lacey blue” (S.OBRL.01)trnH-psbAMT815871Salvia miltiorrhizaKC473184.195.51%
Salvia yangii “Blue spire” (S.OBRL.03)ITS2MT599313Perovskia atriplicifolia
(Salvia yangii)
KJ584242.199.63%
Salvia yangii “Blue spire” (S.OBRL.03)trnH-psbAMT815872Salvia californicaKP852621.196.61%
Salvia yangii (S.OBRL.11)ITS2MT599314Salvia yangiiDQ667223.199.62%
Salvia yangii (S.OBRL.11)trnH-psbAMT815873Salvia miltiorrhizaKC473184.196.83%
Salvia abrotanoides (S.OBRL.12)ITS2MT599315Salvia yangiiDQ667223.199.64
Salvia abrotanoides (S.OBRL.12)trnH-psbAMT815874Salvia chionopeplicaKP852626.195.93%
Table 5. DNA sequences of the ITS2 and trnH-psbA regions of 23 Salvia species used for the NJ analysis. Sequences of Lamium maculatum and Lamium amplexicaule were used as an outgroup.
Table 5. DNA sequences of the ITS2 and trnH-psbA regions of 23 Salvia species used for the NJ analysis. Sequences of Lamium maculatum and Lamium amplexicaule were used as an outgroup.
No.SpeciesSpecimen VoucherITS2trnH-psbA
1Salvia aethiopisx142 *DQ667272.1DQ667370.1
2Salvia amplexicaulisC756MG824168.1MG823900.1
3Salvia dorystaechasx108 *DQ667252.1DQ667360.1
4Salvia miltiorrhizaSMI4JQ934135.1JQ934199.1
5Salvia nipponicaY. Ibaragi s. nAB295103.1MG823977.1
6Salvia nubicolaFLPH 12-124MG824236.1MG823979.1
7Salvia officinalisPS1700MT02FJ883522.1FJ513122.1
8Salvia paramiltiorrhizaSPAR 3JQ934142.1JQ934206.1
9Salvia przewalskiiSPRZ 1JQ934153.1JQ934217.1
10Salvia rosmarinusx074 *DQ667241.1FJ493283.1
11Salvia sclareaPS1701MT01FJ883529.1FJ513083.1
12Salvia koyamaen/dMK425922.1MG823960.1
13Salvia campanulataPS1727MT01FJ883500.1FJ513135.1
14Salvia bulleyanaYin et al. 1314MG824179.1FJ513128.1
15Salvia taraxacifoliax001 *DQ667209.1DQ667337.1
16Salvia judaicaC769MG824218.1MG823957.1
17Salvia lyrataG.X. Hu QT011MG824225.1MG823966.1
18Salvia disjunctaJ. Walker 3018 (WIS)MF622126.1MF623945.1
19Salvia moniliformisCrone 15/9/00 (MJG)MF622160.1MF623979.1
20Salvia polystachyaMexico. Morelos F. Sazatornil 4 (MEXU)DQ667292.1MF663893.1
21Salvia patensx109 *DQ667253.1DQ667361.1
22Salvia funereaJBW 3131KP852812.1KP852638.1
23Salvia aristatax170 *DQ667280.1DQ667375.1
24Lamium maculatumn/dKF055055.1HE966679.1
25Lamium amplexicaulen/dJX073976.1HQ966679.1
n/d—not determined * isolate numbers retrieved from Walker and Sytsma, 2007 [12].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bielecka, M.; Pencakowski, B.; Stafiniak, M.; Jakubowski, K.; Rahimmalek, M.; Gharibi, S.; Matkowski, A.; Ślusarczyk, S. Metabolomics and DNA-Based Authentication of Two Traditional Asian Medicinal and Aromatic Species of Salvia subg. Perovskia. Cells 2021, 10, 112. https://doi.org/10.3390/cells10010112

AMA Style

Bielecka M, Pencakowski B, Stafiniak M, Jakubowski K, Rahimmalek M, Gharibi S, Matkowski A, Ślusarczyk S. Metabolomics and DNA-Based Authentication of Two Traditional Asian Medicinal and Aromatic Species of Salvia subg. Perovskia. Cells. 2021; 10(1):112. https://doi.org/10.3390/cells10010112

Chicago/Turabian Style

Bielecka, Monika, Bartosz Pencakowski, Marta Stafiniak, Klemens Jakubowski, Mehdi Rahimmalek, Shima Gharibi, Adam Matkowski, and Sylwester Ślusarczyk. 2021. "Metabolomics and DNA-Based Authentication of Two Traditional Asian Medicinal and Aromatic Species of Salvia subg. Perovskia" Cells 10, no. 1: 112. https://doi.org/10.3390/cells10010112

APA Style

Bielecka, M., Pencakowski, B., Stafiniak, M., Jakubowski, K., Rahimmalek, M., Gharibi, S., Matkowski, A., & Ślusarczyk, S. (2021). Metabolomics and DNA-Based Authentication of Two Traditional Asian Medicinal and Aromatic Species of Salvia subg. Perovskia. Cells, 10(1), 112. https://doi.org/10.3390/cells10010112

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop