Biological Activity and NMR-Fingerprinting of Balkan Endemic Species Stachys thracica Davidov

Stachys thracica Davidov is a Balkan endemic species distributed in Bulgaria, Greece, and Turkey. In Bulgaria, it is classified as “rare” and is under the protection of the Bulgarian biodiversity law. The aim of our study was to develop an efficient protocol for ex situ conservation of S. thracica and to perform comparative NMR-based metabolite profiling and bioactivity assays of extracts from in situ grown, in vitro cultivated, and ex vitro acclimated plants. Micropropagation of S. thracica was achieved by in vitro cultivation of mono-nodal segments on basal MS medium. Ex vitro adaptation was accomplished in the experimental field with 83% survival while conserved genetic identity between in vitro and ex vitro plants as shown by the overall sequence-related amplified polymorphism marker patterns was established. Verbascoside, chlorogenic acid, and trigonelline appeared the main secondary metabolites in in situ, in vitro cultivated, and ex vitro acclimated S. thracica. High total phenolic and flavonoid content as well as antioxidant and radical scavenging activity were observed in in situ and ex vitro plants. Further, the anti-inflammatory activity of S. thracica was tested by hemolytic assay and a high inhibition of the complement system was observed. Initiated in vitro and ex vitro cultures offer an effective tool for the management and better exploitation of the Stachys secondary metabolism and the selection of lines with high content of bioactive molecules and nutraceuticals.


Introduction
Stachys L., or woundwort, is a subcosmopolitan genus of herbs and shrubs that comprises more than 300 species [1] and is considered as one of the largest genera among the Liliaceae family. The plants from the genus are distributed in temperate and tropical regions of the Mediterranean, Asia, America, and Southern Africa. In Bulgaria, 22 species of the genus Stachys are naturally distributed, as five of them are under the protection of the Bulgarian biodiversity law [2]. Since ancient times, species of this genus have been used in traditional medicine under the form of extracts, decoctions, and ointments for the treatment of stomach disorders, genital tumors, skin inflammations, sclerosis of the spleen, cough, and ulcers [3,4]. In recent years, several studies on different taxa of genus Stachys have demonstrated that woundworts exert various biological effects, such as antioxidant,

In Vitro Cultivation and Ex Vitro Acclimation of S. thracica
Bulgaria is characterized by a high level of endemism-approximately 270 species are considered Balkan endemics [16]. Most of those species have medicinal qualities and are a potential source of nutraceuticals and functional ingredients but as yet they have not been thoroughly studied. Micropropagation (in vitro cultivation) represents a biotechnological approach for ex situ conservation of rare and endangered plants, which can be carried out in a controlled environment regardless of seasons and climate changes [17]. Low seed germination of S. thracica seeds was found, as after 14 days only 20% of the seeds cultivated on 0.7% water agar germinated, while these inoculated on half-strength MS medium did not germinate. The sprouting seedlings were then transferred to a basal MS medium and cultivated under controlled environmental conditions. After one month, approximately, 78% of the explants showed shoot proliferation and produced 11.14 ± 0.73 shoots per explant. The regenerated micro-plants managed to form thick and healthy roots and plentiful leaf biomass ( Figure 1). Metabolites 2022, 12, x FOR PEER REVIEW 3 of 14 The regenerated in vitro plants with 1 cm long shoots and a well-developed root system were transferred into pots with a sterile soil mixture. The first stage of acclimation was performed in a phytotron chamber under controlled environmental conditions and a gradual decrease in relative humidity. The survival rate of S. thracica acclimated plants appeared to be relatively high-83%. In the next step, the ex vitro adapted plants were transferred to a greenhouse and adapted for a period of one month followed by transfer to normal garden soil in the experimental field ( Figure 1). The percentage of surviving plants remained unchanged.
Successful in vitro regeneration on hormone free MS medium and subsequent ex vitro adaptation has been previously reported for S. maritima [18]. However, in S. leucoglossa and S. annua, MS medium supplemented with 0.5 mg L −1 6-benzyladenine was most effective for shoot proliferation and produced 5.61 ± 1.15 and 4.5 ± 0.54 shoots per explant respectively [19,20]. Root development in in vitro cultivated S. leucoglossa was established in 60% of the explants cultivated on WPM media supplemented with 0.5 mg L −1 indole-3butyric acid [20].

Genetic Stability of In Vitro Cultivated and Ex Vitro Adapted S. thracica
The SRAP approach, developed by Li and Quiros [21], is a simple, effective, and reliable marker system that can be adapted for a variety of purposes. The SRAP markers have been successfully applied for analysis of genetic diversity in plant genetic resources and natural populations from the Lamiaceae family [22,23]. Alekseeva et al. [23] studied genetic diversity in Origanum vulgare subsp. hirtum and reported that the SRAP assay is sensitive in the detection of heterogeneity between individuals within populations of the same species. As well, comparison of SRAP profiles between different lavender varieties indicated high genetic diversity [22]. However, it has been pointed out that in the case of vegetatively propagated plants from one lavender variety, SRAP analysis with three primer pairs displayed the same pattern of SRAP fragments. In our work, a total of 16 combinations of SRAP primer pairs was selected for successful PCR amplification and generation of fragments for the polymorphism study of S. thracica plant variants. These primers showed sufficient efficiency to distinguish between S. thracica and other Stachys species (unpublished data). On the other hand, the comparison between in vitro and respective ex vitro S. thracica individuals demonstrated the presence of overall 496 different alleles The regenerated in vitro plants with 1 cm long shoots and a well-developed root system were transferred into pots with a sterile soil mixture. The first stage of acclimation was performed in a phytotron chamber under controlled environmental conditions and a gradual decrease in relative humidity. The survival rate of S. thracica acclimated plants appeared to be relatively high-83%. In the next step, the ex vitro adapted plants were transferred to a greenhouse and adapted for a period of one month followed by transfer to normal garden soil in the experimental field ( Figure 1). The percentage of surviving plants remained unchanged.
Successful in vitro regeneration on hormone free MS medium and subsequent ex vitro adaptation has been previously reported for S. maritima [18]. However, in S. leucoglossa and S. annua, MS medium supplemented with 0.5 mg L −1 6-benzyladenine was most effective for shoot proliferation and produced 5.61 ± 1.15 and 4.5 ± 0.54 shoots per explant respectively [19,20]. Root development in in vitro cultivated S. leucoglossa was established in 60% of the explants cultivated on WPM media supplemented with 0.5 mg L −1 indole-3butyric acid [20].

Genetic Stability of In Vitro Cultivated and Ex Vitro Adapted S. thracica
The SRAP approach, developed by Li and Quiros [21], is a simple, effective, and reliable marker system that can be adapted for a variety of purposes. The SRAP markers have been successfully applied for analysis of genetic diversity in plant genetic resources and natural populations from the Lamiaceae family [22,23]. Alekseeva et al. [23] studied genetic diversity in Origanum vulgare subsp. hirtum and reported that the SRAP assay is sensitive in the detection of heterogeneity between individuals within populations of the same species. As well, comparison of SRAP profiles between different lavender varieties indicated high genetic diversity [22]. However, it has been pointed out that in the case of vegetatively propagated plants from one lavender variety, SRAP analysis with three primer pairs displayed the same pattern of SRAP fragments. In our work, a total of 16 combinations of SRAP primer pairs was selected for successful PCR amplification and generation of fragments for the polymorphism study of S. thracica plant variants. These primers showed sufficient efficiency to distinguish between S. thracica and other Stachys species (unpublished data). On the other hand, the comparison between in vitro and respective ex vitro S. thracica individuals demonstrated the presence of overall 496 different alleles without any variation in the SRAP profiles, which suggests preserved genetic stability during the process of micropropagation (Supplementary Table S1).

NMR-Based Metabolite Profiling during S. thracica Ex Situ Conservation
In order to reveal the metabolic alterations in S. thracica plants during the process of ex situ conservation an 1 H NMR in combination with 2D NMR techniques (J-resolved, COSY, TOCSY, and HSQC) were applied. In total 15 individual compounds were unambiguously assigned to the abundant signals including carbohydrates, amino acids, organic acids, phenolic compounds, and alkaloids. No major qualitative differences were observed in the metabolic profiles of in situ, in vitro cultivated, and ex vitro adapted plants except for the two unidentified phenolic compounds, whose signals were detected only in the NMR spectra of ex vitro acclimated plants ( Figure 2, Table 1). Verbascoside and leucoseptoside A along with chlorogenic acid and trigonelline appeared to be the main secondary metabolites in the extracts of all samples. Phenylethanoid glycosides are of common occurrence in the members of Stachys genus [4]. Verbascoside and leucoseptoside A have also been reported in S. officinalis, S. recta, S. affinis, S. alpina subsp. dinarica, S. anisochila, S. beckeana, S. byzantine, S. plumose, S. iva, S. candida, S. schtschegleevii, and others [1,4]. It has been found that verbascoside possesses a wide range of biological properties including antioxidant, anti-inflammatory, cytotoxic, antimicrobial, anti-thrombotic, and wound healing [24]. The identification of phenylethanoids was confirmed by comparison with authentic samples. without any variation in the SRAP profiles, which suggests preserved genetic stability during the process of micropropagation (Supplementary Table S1).

NMR-Based Metabolite Profiling during S. thracica Ex Situ Conservation
In order to reveal the metabolic alterations in S. thracica plants during the process of ex situ conservation an 1 H NMR in combination with 2D NMR techniques (J-resolved, COSY, TOCSY, and HSQC) were applied. In total 15 individual compounds were unambiguously assigned to the abundant signals including carbohydrates, amino acids, organic acids, phenolic compounds, and alkaloids. No major qualitative differences were observed in the metabolic profiles of in situ, in vitro cultivated, and ex vitro adapted plants except for the two unidentified phenolic compounds, whose signals were detected only in the NMR spectra of ex vitro acclimated plants ( Figure 2, Table 1). Verbascoside and leucoseptoside A along with chlorogenic acid and trigonelline appeared to be the main secondary metabolites in the extracts of all samples. Phenylethanoid glycosides are of common occurrence in the members of Stachys genus [4]. Verbascoside and leucoseptoside A have also been reported in S. officinalis, S. recta, S. affinis, S. alpina subsp. dinarica, S. anisochila, S. beckeana, S. byzantine, S. plumose, S. iva, S. candida, S. schtschegleevii, and others [1,4]. It has been found that verbascoside possesses a wide range of biological properties including antioxidant, anti-inflammatory, cytotoxic, antimicrobial, anti-thrombotic, and wound healing [24]. The identification of phenylethanoids was confirmed by comparison with authentic samples. Chlorogenic acid is another major metabolite present with high intensity in the NMR spectra of all S. thracica samples. It is an important and biologically active dietary polyphenol found in the Asteraceae and Lamiaceae families [25]. Chlorogenic acid offers valuable therapeutic properties including antioxidant activity, antibacterial, hepatoprotective, cardioprotective, anti-inflammatory, antipyretic, neuroprotective, anti-obesity, antiviral, anti-microbial, anti-hypertension, and a central nervous system stimulation [25]. It has been also found that chlorogenic acid could modulate lipid metabolism and glucose in both genetic and healthy metabolic disorders [25].
The biological activities of trigonelline, a vitamin B6 derivative, have been thoroughly evaluated especially regarding diabetes and central nervous system disease. It possesses also hypoglycaemic, hypolipidemic, neuroprotective, antimigraine, sedative, Chlorogenic acid is another major metabolite present with high intensity in the NMR spectra of all S. thracica samples. It is an important and biologically active dietary polyphenol found in the Asteraceae and Lamiaceae families [25]. Chlorogenic acid offers valuable therapeutic properties including antioxidant activity, antibacterial, hepatoprotective, cardioprotective, anti-inflammatory, antipyretic, neuroprotective, anti-obesity, antiviral, anti-microbial, anti-hypertension, and a central nervous system stimulation [25]. It has been also found that chlorogenic acid could modulate lipid metabolism and glucose in both genetic and healthy metabolic disorders [25]. Table 1. Chemical shifts (δ) and coupling constants (J) of the metabolites, identified in S. thracica plants by the relevant 1D and 2D NMR spectra.

Amino acids
Alanine

Organic acids
Acetic acid

Phenolic acids
Chlorogenic acid ++ ++ +++ The biological activities of trigonelline, a vitamin B6 derivative, have been thoroughly evaluated especially regarding diabetes and central nervous system disease. It possesses also hypoglycaemic, hypolipidemic, neuroprotective, antimigraine, sedative, and memoryimproving activities, and may also reduce diabetic auditory neuropathy and platelet aggregation. It has been established that trigonelline affects β-cell regeneration, insulin secretion, activities of enzymes related to glucose metabolism, reactive oxygen species, axonal extension, and neuron excitability [26].

Comparative Determination of Total Phenols and Flavonoids in In Situ, In Vitro Cultivated, and Ex Vitro Adapted Plants
Phenolic compounds are the most diverse and widely distributed secondary metabolites in plant species with pronounced antioxidant activity. They are synthesized mainly in the presence of adverse environmental factors (biotic and abiotic) and their amount varies significantly depending on the growth conditions. The determination of the total phenolic and flavonoid content in in situ, in vitro cultivated and ex vitro adapted S. thracica plants provides information on the quantitative changes in metabolic content in plants of the same genotype which are cultivated under different conditions. The highest phenolic content was established in in situ wild plants, followed by ex vitro adapted S. thracica (171 ± 4.5 and 161.6 ± 1.24 µg GAE mg −1 extract, respectively; Figure 3a). Reverse dependence was observed in the flavonoid content, which appeared highest in ex vitro acclimated plants followed by in situ grown (41.7 ± 0.68 and 29.7 ± 0.65 µg QE mg −1 extract, respectively; Figure 3b). A nearly fourfold decrease in the content of total phenols and a twofold decrease of flavonoids in in vitro cultivated S. thracica compared to wild growing genotypes was observed ( Figure 3). It is likely that the aseptic culture conditions in controlled environment and mixotrophic nutrition affect the quantity of phenolic compounds in micropropagated S. thracica plants.

Comparative Determination of Total Phenols and Flavonoids in In Situ, In Vitro Cultivated, and Ex Vitro Adapted Plants
Phenolic compounds are the most diverse and widely distributed secondary metabolites in plant species with pronounced antioxidant activity. They are synthesized mainly in the presence of adverse environmental factors (biotic and abiotic) and their amount varies significantly depending on the growth conditions. The determination of the total phenolic and flavonoid content in in situ, in vitro cultivated and ex vitro adapted S. thracica plants provides information on the quantitative changes in metabolic content in plants of the same genotype which are cultivated under different conditions. The highest phenolic content was established in in situ wild plants, followed by ex vitro adapted S. thracica (171 ± 4.5 and 161.6 ± 1.24 μg GAE mg -1 extract, respectively; Figure 3a). Reverse dependence was observed in the flavonoid content, which appeared highest in ex vitro acclimated plants followed by in situ grown (41.7 ± 0.68 and 29.7 ± 0.65 μg QE mg -1 extract, respectively; Figure 3b). A nearly fourfold decrease in the content of total phenols and a twofold decrease of flavonoids in in vitro cultivated S. thracica compared to wild growing genotypes was observed ( Figure 3). It is likely that the aseptic culture conditions in controlled environment and mixotrophic nutrition affect the quantity of phenolic compounds in micropropagated S. thracica plants. Strong inhibition of metabolic activity was also observed during in vitro cultivation of Lamium album and Achillea thracica, however, it recovered upon ex vitro adaptation [17,27]. In several plant model systems, a significant decrease in the concentration of total phenolic and flavonoid content was also found in in vitro cultivated plants compared to those grown in situ or adapted ex vitro [28,29].

Antioxidant and Radical Scavenging Activity of S. thracica during the Process of Ex Situ Conservation
Phenolics have been recognized as essential antioxidant agents due to their structural characteristics and chemical behaviour. Based on the hydrogen-donating ability, they may act as free-radical scavengers and, consequently, exert a protective effect against reactive oxygen species (ROS) [30]. A positive correlation between the phenolic content and total antioxidant activity (TAA) was observed as the highest TAA was established in in situ and ex vitro acclimated plants (0.190 ± 0.005 and 0.197 ± 0.006 mM α-tocopherol g -1 extract, respectively). Approximately, two times lower TAA was found in in vitro cultivated S. thracica plants (0.104 ± 0.004 mM α-tocopherol g -1 extract; Figure 4a). Similar dependence was observed in the FRAP assay. The ferric reducing activity was highest in in situ wild and ex vitro adapted plants (2.23 ± 0.05 and 2.07 ± 0.07 mM Fe 2+ , respectively), and Strong inhibition of metabolic activity was also observed during in vitro cultivation of Lamium album and Achillea thracica, however, it recovered upon ex vitro adaptation [17,27]. In several plant model systems, a significant decrease in the concentration of total phenolic and flavonoid content was also found in in vitro cultivated plants compared to those grown in situ or adapted ex vitro [28,29].

Antioxidant and Radical Scavenging Activity of S. thracica during the Process of Ex Situ Conservation
Phenolics have been recognized as essential antioxidant agents due to their structural characteristics and chemical behaviour. Based on the hydrogen-donating ability, they may act as free-radical scavengers and, consequently, exert a protective effect against reactive oxygen species (ROS) [30]. A positive correlation between the phenolic content and total antioxidant activity (TAA) was observed as the highest TAA was established in in situ and ex vitro acclimated plants (0.190 ± 0.005 and 0.197 ± 0.006 mM α-tocopherol g −1 extract, respectively). Approximately, two times lower TAA was found in in vitro cultivated S. thracica plants (0.104 ± 0.004 mM α-tocopherol g −1 extract; Figure 4a). Similar dependence was observed in the FRAP assay. The ferric reducing activity was highest in in situ wild and ex vitro adapted plants (2.23 ± 0.05 and 2.07 ± 0.07 mM Fe 2+ , respectively), and threefold lower in in vitro cultivated plants (0.598 ± 0.020 mM Fe 2+ ; Figure 4b). In addition, the extracts from in situ and ex vitro acclimated plants have nearly two times higher antioxidant potential of the used standard α-tocopherol (Figure 4b). In DPPH assay all the three methanolic extracts exhibited concentration dependant radical scavenging activity. The maximum inhibition of the DPPH free radical was 76% and 74% at a concentration of 80 µg mL −1 for in situ and ex vitro adapted plants and 64% for in vitro cultivated S. thracica at the highest tested concentration-150 µg mL −1 (Figure 4c). The lowest concentrations at which 50% inhibition (IC 50 ) of the DPPH radical was observed were 8.9 µg mL −1 for in situ, 24.4 µg mL −1 for ex vitro adapted, and 96.3 µg mL −1 for in vitro cultivated plants. The highest ABTS radical scavenging activity was detected in in situ wild and ex vitro adapted plants (1.02 ± 0.02 and 0.91 ± 0.03 mg TE mg −1 extract, respectively) and fivetimes lower activity was found in in vitro cultivated S. thracica (0.18 ± 0.02 mg TE mg −1 extract; Figure 4d).  Figure 4b). In addition, the extracts from in situ and ex vitro acclimated plants have nearly two times higher antioxidant potential of the used standard α-tocopherol (Figure 4b). In DPPH assay all the three methanolic extracts exhibited concentration dependant radical scavenging activity. The maximum inhibition of the DPPH free radical was 76% and 74% at a concentration of 80 μg mL -1 for in situ and ex vitro adapted plants and 64% for in vitro cultivated S. thracica at the highest tested concentration-150 μg mL -1 (Figure 4c). The lowest concentrations at which 50% inhibition (IC50) of the DPPH radical was observed were 8.9 μg mL -1 for in situ, 24.4 μg mL -1 for ex vitro adapted, and 96.3 μg mL -1 for in vitro cultivated plants. The highest ABTS radical scavenging activity was detected in in situ wild and ex vitro adapted plants (1.02 ± 0.02 and 0.91 ± 0.03 mg TE mg -1 extract, respectively) and five-times lower activity was found in in vitro cultivated S. thracica (0.18 ± 0.02 mg TE mg -1 extract; Figure  4d). Most likely the high antioxidant potential of in situ and ex vitro adapted S. thracica plants is due to the presence of verbascoside and chlorogenic acid, which are the main secondary metabolites found in the extracts through NMR profiling (Table 1). There is considerable evidence showing that phenylethanoids and phenolic acid are powerful antioxidants by scavenging ROS directly or acting as chain-breaking peroxyl radical scavengers [1,31]. It has been established that the major components, responsible for the high radical scavenging potential, in the methanolic extract of S. officinalis are verbascoside and chlorogenic acid, comprising 69% of the total antioxidant activity [31]. Methanolic extracts of aerial flowering parts of four Stachys taxa (S. alpina subsp. dinarica, S. anisochila, S. beckeana and S. plumose) were studied for their antioxidant activity and high correlations Most likely the high antioxidant potential of in situ and ex vitro adapted S. thracica plants is due to the presence of verbascoside and chlorogenic acid, which are the main secondary metabolites found in the extracts through NMR profiling (Table 1). There is considerable evidence showing that phenylethanoids and phenolic acid are powerful antioxidants by scavenging ROS directly or acting as chain-breaking peroxyl radical scavengers [1,31]. It has been established that the major components, responsible for the high radical scavenging potential, in the methanolic extract of S. officinalis are verbascoside and chlorogenic acid, comprising 69% of the total antioxidant activity [31]. Methanolic extracts of aerial flowering parts of four Stachys taxa (S. alpina subsp. dinarica, S. anisochila, S. beckeana and S. plumose) were studied for their antioxidant activity and high correlations between total phenolic content, TAA and DPPH scavenging activity was established, indicating that polyphenols are the main antioxidants [12].

Anti-Inflammatory Activity of S. thracica
The anti-inflammatory potential was examined by monitoring the effect of methanolic extracts of S. thracica on the complement hemolytic activity via the classical pathway (CP). The complement system is an innate effector mechanism with the intrinsic ability to initiate local inflammation at the sites of accumulation of ligands for its classical recognition molecule C1q, typically immune complexes, as well as to augment significantly an ongoing inflammatory response by the generation of plethora of pro-inflammatory mediators in a cascade mode. Therefore, inhibiting the complement activation (e.g., CP) would produce a profound anti-inflammatory effect. After treatment with the S. thracica extracts a dosedependent inhibition on the CP was observed. The extracts from in situ and ex vitro adapted plants exhibited a similar degree of inhibition reaching up to 94% and 97%, respectively, at a concentration of 2000 µg mL −1 , while those from in vitro cultivated plants reached a maximum of 69% ( Figure 5). The lowest concentrations at which 50% inhibition (IC50) of the hemolysis was observed were 351.7 µg mL −1 for in situ, 358.5 µg mL −1 for ex vitro adapted, and 872 µg mL −1 for in vitro cultivated plants. The extracts from in situ and ex vitro acclimated plants had not only similar trends of inhibition on the CP but also a comparable anti-inflammatory effectiveness. The extract from in vitro cultivated S. thracica was 2.5-fold less active than the rest. It is very apparent that environmental factors significantly reduce the concentration of phenolic compounds and hence the biological activity of in vitro cultivated Thracian woundwort, while in ex vitro conditions, the biosynthetic capacity is restored, and this leads to an increase in bioactivity of the plants.
between total phenolic content, TAA and DPPH scavenging activity was established, indicating that polyphenols are the main antioxidants [12].

Anti-Inflammatory Activity of S. thracica
The anti-inflammatory potential was examined by monitoring the effect of methanolic extracts of S. thracica on the complement hemolytic activity via the classical pathway (CP). The complement system is an innate effector mechanism with the intrinsic ability to initiate local inflammation at the sites of accumulation of ligands for its classical recognition molecule C1q, typically immune complexes, as well as to augment significantly an ongoing inflammatory response by the generation of plethora of pro-inflammatory mediators in a cascade mode. Therefore, inhibiting the complement activation (e.g., CP) would produce a profound anti-inflammatory effect. After treatment with the S. thracica extracts a dose-dependent inhibition on the CP was observed. The extracts from in situ and ex vitro adapted plants exhibited a similar degree of inhibition reaching up to 94% and 97%, respectively, at a concentration of 2000 μg mL -1 , while those from in vitro cultivated plants reached a maximum of 69% ( Figure 5). The lowest concentrations at which 50% inhibition (IC50) of the hemolysis was observed were 351.7 μg mL -1 for in situ, 358.5 μg mL -1 for ex vitro adapted, and 872 μg mL -1 for in vitro cultivated plants. The extracts from in situ and ex vitro acclimated plants had not only similar trends of inhibition on the CP but also a comparable anti-inflammatory effectiveness. The extract from in vitro cultivated S. thracica was 2.5-fold less active than the rest. It is very apparent that environmental factors significantly reduce the concentration of phenolic compounds and hence the biological activity of in vitro cultivated Thracian woundwort, while in ex vitro conditions, the biosynthetic capacity is restored, and this leads to an increase in bioactivity of the plants. Maleki-Dizaji et al. [32] reported on the anti-inflammatory properties of a hydro alcoholic extract from S. schtschegleevii and associated this activity with the presence of phenylethanoid glycosides (verbascoside and betonyoside F) and flavonoids. It has been found that the polyphenol-rich extract from S. officinalis inhibited the enzymes lipoxygenase and cyclooxygenase-2 with IC50-1.22 μg mL -1 and 10.1 μg mL -1 , respectively [33]. Haznagy-Radnai et al. [6] also reported high anti-inflammatory activity of aqueous extracts from S. alpina, S. germanica, S. officinalis, and S. recta in the carrageenan-induced paw Maleki-Dizaji et al. [32] reported on the anti-inflammatory properties of a hydro alcoholic extract from S. schtschegleevii and associated this activity with the presence of phenylethanoid glycosides (verbascoside and betonyoside F) and flavonoids. It has been found that the polyphenol-rich extract from S. officinalis inhibited the enzymes lipoxygenase and cyclooxygenase-2 with IC 50 -1.22 µg mL −1 and 10.1 µg mL −1 , respectively [33]. Haznagy-Radnai et al. [6] also reported high anti-inflammatory activity of aqueous extracts from S. alpina, S. germanica, S. officinalis, and S. recta in the carrageenan-induced paw oedema in rats. The extracts showed greater potency compared to the positive control diclofenac-Na.

Chemicals
Murashige & Skoog medium, plant agar, and sucrose used for in vitro cultivation were purchased from Duchefa Biochemie (Haarlem, The Netherlands). The organic solvents used for extraction, the reagents, and the standards used for the determination of total phenols and flavonoids, as well as the DPPH, ABTS, and FRAP free radicals were purchased from Sigma-Aldrich (Madrid, Spain). CD 3 OD and D 2 O came from Deutero GmbH (Kastellaun, Germany). The sensitized red sheep erythrocytes, hemolysin and the guinea pig serum were purchased from BulBio (Sofia, Bulgaria).

Plant Material and Culture Conditions
Wild grown in situ Stachys thracica Davidov plants were collected from their natural habitat (village of Sinemorets, Tsarevo municipality, Bulgaria) in the period of blooming in June and seeds in September, with the permission of the Ministry of Environment and Water of Bulgaria. A voucher specimen SO107847 was deposited in the Herbarium of Sofia University "St. Kliment Ohridski". In Vitro shoot culture was induced by sterilization of 100 ripe dry seeds with 70% ethanol as described by Yordanova et al. [27]. After that, sterilized seeds were inoculated on half-strength Murashige and Skoog (MS) medium [34] and on 0.7% water agar (w/v). After 14 days the seedlings were transferred on MS medium with 3% (w/v) sucrose and 0.7% (w/v) agar and cultivated under controlled environmental conditions (80 µmol m −2 s −1 photosynthetic active radiation, cool white fluorescent TL-D 36W/54-765 1SL/25 Philips, photoperiod 16 h light/8 h dark, 25 ± 1 • C, 50-60% moderate humidity). The plants were micropropagated twice over a period of 30 days before performing further experiments.

Ex Vitro Acclimation
In Vitro regenerated S. thracica plants with well-developed root systems were transferred to plastic pots containing a sterile soil mixture (peat:coconut-fibers:sand = 2:1:1). The first stage of ex vitro adaptation was performed for one month in a phytotron chamber (POL-EKO APARATURA SP.J.A. Polok-Kowalska KK 350 STD 1400 W) with 16/8 h light/dark, 100 µmol m −2 s −1 PPFD, 22 ± 2 • C, the relative humidity was decreased from 90% to 60% every week. In the next step, the ex vitro adapted plants were transferred to a greenhouse and adapted for a period of one month followed by transfer to normal garden soil at the experimental field of Sofia University "St. Kliment Ohridski". After one year of acclimation to the field conditions, newly formed fully expanded leaves from the 2nd or 3rd nodes of the stem of ex vitro plants in the period of blooming were harvested and used for further analyses and NMR-metabolic profiling.

Genetic Stability Assay by SRAP Markers
To check genetic stability between in vitro and respective ex vitro adapted plants, the sequence-related amplified polymorphism (SRAP) approach was applied according to [21]. Frozen leaf material was used for genomic DNA purification and for subsequent PCR amplification, the primer pair sequences (ME1 to10 for forward primers targeting exons, and EM1 to 10 reverse primers targeting DNA non-coding region) used by Zagorcheva et al. [22] were utilized. In total, 16 SRAP primer pairs, which showed the highest number of peaks and clarity of the electropherogram, were selected for further use in the current study. The following primer pair combinations were used: ME1 + EM3; ME1 + EM5; ME1 + EM6; ME1 + EM7; ME1 + EM10; ME3 + EM2; ME3 + EM5; ME4 + EM5; ME6 + EM3; ME6 + EM5; ME7 + EM2; ME8 + EM1; ME8 + EM8; ME10 + EM4; ME10 + EM7; ME10 + EM9. The forward ME primers were 5 end labelled with FAM (carboxyfluorescein) dye. PCR reactions were carried out with 50 ng genomic DNA, while the PCR amplification steps and fragment analysis were performed as described by [22]. The analysis of the length of the generated SRAP fragments was estimated by the GeneMapper Analysis Software v4.0 (Thermo Fisher Scientific Inc., Waltham, MA, USA). The threshold for peak detection was set to 400 relative fluorescence units. Fragment analysis was carried out for fragment sizes in the range of 50 to 950 bp.

Extraction Procedure and NMR Analyses
Extraction of the plant material and NMR analysis were described in detail by Zahmanov et al. [35]. Proton ( 1 H) as well as 2D NMR spectra (J-resolved, COSY, HSQC, and TOCSY) were recorded at 25 • C on an AVII+ 600 spectrometer (Bruker, Karlsruhe, Germany), operating at a proton NMR frequency of 600.01 MHz [35]. Deuterated methanol was used for internal lock. The resulting 1 H NMR spectra for each sample were phased, baseline corrected, and referenced to the residual signal of methanol-d 4 at 3.30 ppm, by running TopSpin software (3.6.5, Bruker BioSpin Group).

Methanolic Extract Preparation
Three grams (3 g) of finely powdered dry plant material from aboveground parts of in situ grown, in vitro cultivated, and ex vitro adapted S. thracica was subjected to triple extraction with 30 mL chloroform (Sigma-Aldrich, Madrid, Spain) in an ultrasonic bath for 10 min. At the next step the biomass was extracted three times by distillation with methanol for 30 min. The final plant extract from each variant was concentrated through vacuum evaporator (IKA, Germany) and dried to constant dry weight. The yields of extracts from in situ, in vitro cultivated, and ex vitro adapted plants were 13.8%, 28.46%, and 13.6% respectively.

Determination of Total Phenolic and Flavonoid Content
The total phenolic content of S. thracica methanol extracts was determined using Folin-Ciocalteu reagent, according to the methodology proposed by Singleton et al. [36]. The determination of flavonoids was performed following the methodology of Chang et al. [37]. The content of polyphenols was quantified by a standard curve using gallic acid as a standard, and expressed as µg gallic acid equivalents per mg (µg GAE mg −1 extract). The flavonoid concentration was quantified by a standard curve using quercetin as a standard, and expressed as µg quercetin equivalents per mg extract (µg QE mg −1 extract).

Total Antioxidant Activity
Total antioxidant activity (TAA) was determined by the method of [38] with modifications. The incubation medium contains 0.25 mL extract (1 mg mL −1 ), 2.5 mL reagent (0.6 M H 2 SO 4 , 28 mM KH 2 PO 4 , and 4 mM (NH 4 ) 2 MoO 4 ). A control containing 0.25 mL of methanol and 2.5 mL of reagent was set in parallel and incubated under the same conditions. The samples were incubated in a water bath for 90 min at 95 • C. The reaction was quenched after placing the samples on ice and subsequent tempering to room temperature. Absorbance was measured at 695 nm on a spectrophotometer, Shimadzu 1800 UV. The TAA was expressed as mM α-tocopherol per gram extract (mM α-tocopherol g −1 extract).

DPPH Radical Scavenging Activity
The radical scavenging activity of methanolic extracts was determined using stable DPPH (1,1-diphenyl-2-picrylhydrazyl) free radical according to [39]. Eight different concentrations of the selected extracts were tested. Briefly, 0.5 mL of 0.1 mM DPPH solution, dissolved in methanol was added to 1.5 mL of sample. Methanol was used as a blank and a mixture of methanol (1.5 mL) with DPPH (0.5 mL) as a positive control. The samples were incubated for 30 min in the dark and the absorbance was measured at 517 nm. The results are given as maximum % inhibition and 50% inhibition (IC 50 ) concentration in µg mL −1 .

ABTS Radical Scavenging Activity
The ability of the different extracts to reduce the ABTS (2,20-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) free radical was determined by the method described by Re et al. [40]. The active ABTS free radical was produced by the reaction of 7 mM ABTS and 2.5 mM K 2 S 2 O 8 in water. The solution is stored in the dark at room temperature for 12 to 16 h. The working ABTS solution was prepared by dilution of the active radical in methanol until the absorbance reached 0.7 ± 0.03 at 734 nm. An aliquot of 20 µL of each extract was mixed with 2.0 mL of the ABTS solution and after 1 min the absorbance was measured on a spectrophotometer, Shimadzu UV1800. Trolox was used as a standard. The results were expressed as mg Trolox equivalents per mg extract (mg TE mg −1 extract).

FRAP Assay
The ability of the extracts to reduce ferrous-containing radical was determined according to the method of [41]. The stock solution of the FRAP (Ferric Reducing Antioxidant Power) reagent was prepared by mixing 10 parts of 300 mM acetate buffer (pH = 3.6), 1 part of 10 mM TPTZ dissolved in 40 mM HCl, and 1 part of 20 mM FeCl 3 . The reaction mixture contained 1.8 mL of FRAP reagent, 0.2 mL of dH 2 O, and 0.040 mL of extract. Ten different concentrations of methanolic extracts were tested. After 30 min incubation at 37 • C the absorbance was measured at 593 nm on a Shimadzu UV1800 spectrophotometer. Ferric-reducing antioxidant capacity was represented as mM Fe 2+ using an FeSO 4 standard curve.

Microtitre Hemolytic Complement Assay
The hemolytic complement assay was performed in 96-well flat-bottom microtiter plates. The dried extracts were dissolved in 3% dimethyl sulfoxide (DMSO) and further diluted with barbitone buffered saline, pH 7.2, containing 0.15 mM Ca 2+ (BBS). The assay was performed on 6% sheep erythrocyte (SE) suspension sensitized with rabbit polyclonal anti-SE serum (hemolysin, BulBio, Sofia, Bulgaria) and guinea pig complement (BulBio, Sofia, Bulgaria). Preliminary titration of sera was performed to determine the dilution producing 50% haemolysis of target erythrocytes. The SE (25 µL/well) were sensitized by 30 min incubation with hemolysin (dilution 1:1600, 25 µL/well) at 37 • C. After that complement (125 µL/well of appropriate dilution) and increasing amounts of the analyzed plant extracts (100 µL/well) were added and incubated for 1 h at 37 • C. Next, the microtiter plates were centrifuged at 1000× g for 5 min, and 200 µL of the supernatant from each well was transferred to new 96-well flat-bottom microtiter plates and the absorbance at 540 nm was measured by an ELISA reader (DR-200B, Hiwell Diatek Instruments, Jiangsu, China). Each assay was carried out in triplicate.

Data Analysis
The data are presented as mean ± SE of at least 12 scores (3 repetitions per variant in each of 4 independent sets of experiments). One-way ANOVA followed by Holm-Sidak statistical test with significance level 0.001 was performed with Sigma Plot 11.0 software to estimate the difference between all the variants.

Conclusions
The protocol for micropropagation of S. thracica comprises germination of sterilized ripe dry seeds on water agar and cultivation of sprouting seedlings on hormone free MS medium with 30 g L −1 sucrose for 30 days. Optimal induction of shoot proliferation from mono-nodal segments was achieved on MS medium without supplement of growth regulators. The ex vitro adaptation was accomplished on the experimental field and conserved genetic identity as shown by the overall sequence-related amplified polymorphism DNA marker patterns was established during the process of ex situ conservation. No significant changes in the NMR metabolic profiles, as measured by NMR, were found and verbascoside, leucoseptoside A, as well as chlorogenic acid and trigonelline were the most abundant metabolites in the extracts of S. thracica. As usual the aseptic conditions significantly reduced the quantity of total phenols and flavonoids and hence the antioxidant, radical scavenging, and anti-inflammatory potential of the extracts from in vitro cultivated plants. However, the biosynthetic potential and the associated biological activity were restored after adaptation of the plants to ex vitro conditions. The successful initiation of in vitro and ex vitro cultures is an alternative biotechnological approach for preservation of S. thracica which will also offer an effective tool for the management and better exploitation of the Stachys secondary metabolism and the selection of lines with high content of nutraceuticals and pharmaceutically valuable molecules.