Distribution Patterns of Essential Oil Terpenes in Native and Invasive Solidago Species and Their Comparative Assessment

The importance of invasive Solidago L. species to the environment creates a new approach to controlling their spread through the use of potentially high value raw materials. The aim of this study was to assess the distribution patterns of volatile compounds in the four Solidago spp., by identifying common and species-specific compounds with their potentials, and to confirm the origin of the spontaneous hybrid Solidago × niederederi on the basis of comparative assessment of essential oil (EO) profiles. Plant material in the flowering phase was collected in mixed populations from six different sites. The EOs were isolated separately from the leaf and the inflorescence samples by hydrodistillation for 3 h. The chemical analysis was performed by gas chromatography—mass spectrometry. Multivariate data analysis was employed to explain the interspecies relationships among Solidago spp. The results revealed the similarity among Solidago spp. EO profiles, which were dominated by monoterpenes and oxygenated compound fractions. Solidago spp. differed in species distinctive terpenes and their distribution between accessions and plant parts. Volatile compound patterns confirmed the origin of Solidago × niederederi between Solidago canadensis and Solidago virgaurea, with the higher contribution of alien species than native ones. Correct taxonomic identification of species is highly essential for the targeted collection of raw material from the wild for different applications. Solidago spp. can be considered to be underutilized sources of bioactive secondary metabolites.


Introduction
The use of plant products has grown remarkably in recent years, and research into natural products such as volatile terpenoids has become an important task for future human and animal well-being [1]. Natural products are generally easy to prepare and, due to their natural origin, are environmentally friendly and not financially challenging. In this respect, invasive species are of growing interest as a potential resource for obtaining high value-added products. Solidago canadensis L. (Canadian goldenrod) and Solidago gigantea Aiton (Giant goldenrod), native to North America, are considered to be one of the most aggressive plant invaders, which were introduced to Europe as garden ornamentals in the middle of the 18th century and began to spread in the 19th century [2]. Abandoned previously cultivated and disturbed areas contribute to a rapid and successful invasion and confirm the origin of the S. × niederederi taxon. The objectives of the study were: (1) to assess the distribution of volatile constituents in the populations of native and alien Solidago spp.; (2) to identify intraspecific and interspecific diversity in Solidago spp.; (3) to identify common and species distinctive volatiles; and (4) to confirm the origin of the spontaneous hybrid S. × niederederi on the basis of a multivariate comparative analysis of volatile profiles.

Chemical Profiles of Essential Oils
The EOs of four Solidago spp. were dominated by monoterpenes, with an average of 43.9-74.6% of the total composition of EOs in inflorescences and 39.7-69.1% in leaves. The mean percentage of sesquiterpenes ranged from 15.0 to 35.4% in inflorescences and from 25.5 to 38.2% in leaf EOs. Chemical profiles for EOs of Solidago spp. varied according to the contents of individual compounds and their distribution among accessions and plant parts. Examples of fingerprint profiles for inflorescence EOs are shown in Figure 1.

Chemical Profiles of Essential Oils
The EOs of four Solidago spp. were dominated by monoterpenes, with an average of 43.9-74.6% of the total composition of EOs in inflorescences and 39.7-69.1% in leaves. The mean percentage of sesquiterpenes ranged from 15.0 to 35.4% in inflorescences and from 25.5 to 38.2% in leaf EOs. Chemical profiles for EOs of Solidago spp. varied according to the contents of individual compounds and their distribution among accessions and plant parts. Examples of fingerprint profiles for inflorescence EOs are shown in Figure 1.
Solidago × niederederi leaf EOs contained seven compounds in concentrations above 10% in at least one sample and were considered as principal compounds. The major compounds such as α-pinene (1.1-25.8%), trans-verbenol (2.2-20.2%), bornyl acetate (4.2-20.2%), and caryophyllene oxide (2.1-38.1%) were found in all leaves in highly different concentrations. Verbenone and germacrene D were present in most EOs, averaging 5.2 and 8.8%, respectively, with the exception of the two EOs in which these compounds were predominant, accounting for 20.2 and 18.0%, respectively. Meanwhile, β-cubebene was detected only in two oils (16.4 and 23.9%) in which it was the first or second major compound. One leaf EO contained noteworthy concentrations of thymol (6.4%) and carvacrol (9.5%). In addition, 21 compounds were detected with a mean content of 1 to 5.0%, the most prominent of which were sesquiterpenes such as β-copaene, β-bisabolene, epoxyazulene, humulene epoxide II, and isospathulenol. Similar results for the composition of the major compounds in S. × niederederi EO have been recently published [26]. However, the results of only one plant accession were reported, making the comparison insufficient as the prevalence and number of compounds varied among samples.

Interspecific Differences
The results presented revealed similarities and differences in the frequency of distribution and contents of EO constituents among the four Solidago spp. Significant differences in monoterpene and sesquiterpene fractions were found between the inflorescences of the species studied (Table 1). Meanwhile, the chemical groups of compounds in leaves did not differ significantly between the four Solidago spp. (Table 2).
Oxygenated monoterpenes predominated in the EOs of S. canadensis and S. × niederederi inflorescences and S. virgaurea and S. × niederederi leaves. The highest proportion of sesquiterpenoids among the species was found in S. gigantea inflorescences and leaves. Mono-and sesquiterpenes prevailed in the S. virgaurea inflorescence EOs; however, no significant differences were found between the species for monoterpenes. The common principal constituents of the inflorescence and leaf EOs of all four Solidago species were α-pinene, bornyl acetate, and caryophyllene oxide. Species showed significant differences (p ≤ 0.05) in the accumulation of bornyl acetate and caryophyllene oxide in leaf EOs, but no differences in α-pinene were observed between species (Tables 1 and 2). The highest mean concentration of bornyl acetate was found in inflorescence and leaf EOs of both S. gigantea and S. canadensis, while caryophyllene oxide prevailed in S. × niederederi and S. virgaurea. The other major compound trans-verbenol was prevalent in all inflorescence EOs with the highest level (p ≤ 0.001) in S. canadensis and S. × niederederi. Meanwhile, the leaves of S. virgaurea and S. × niederederi had priority over the other two species in accumulation of trans-verbenol. Inflorescences of S. gigantea accumulated the highest (p ≤ 0.001) level of germacrene D compared to other species, but this compound did not differ significantly (p > 0.05) in the leaves between species. The inflorescences and leaves of S. gigantea were in priority to other species in accumulation of oxygenated sesquiterpenes as epoxyazulene, spathulenol, and isospathulenol. In addition, γ-muurolene, differently other species, was a common compound in all inflorescence EOs of S. gigantea, while other monoterpene hydrocarbons such as camphene and o-cymene were common in all leaf EOs of this species.
The inflorescence EOs of all species differed significantly in the mean concentrations of oxygenated monoterpenes, such as α-campholenal, trans-pinocarveol, cis-verbenol, pinocarvone, verbenone, and trans-carveol, with the highest levels and frequency found in S. canadensis and S. × niederederi, followed by S. virgaurea. The same compounds, with the exception of trans-carveol, differed significantly in leaf EOs, with the highest concentrations in S. virgaurea and S. × niederederi leaves. Meanwhile, trans-carveol was detected in small amounts only in S. canadensis and S. × niederederi leaves. In addition, S. virgaurea inflorescence EOs differed from other species in the highest levels and frequency of distribution in α-and β-copaene, cubebol, α-muurolene, and δ-cadinene. As a consequence, quantitative rather than qualitative differences were observed between the species EOs. However, 10 compounds common in more than 30% of all studied inflorescence and/or leaf EOs did not differ significantly between Solidago spp. (Tables 1 and 2). Among them, the most abundant were α-and β-pinene, limonene, bornyl acetate, germacrene D, and caryophyllene oxide.

Principal Component Analysis (PCA)
PCA was employed to explain the phytochemical relationships arising due to interand intraspecific differences between the four Solidago spp., using the selected EOs compounds. A scree plot criterion was applied to reduce the number of PCs for explaining the variance in the selected variables. A two-dimensional PCA square matrix model explained more than 53.3% of the total variance and was used to visualize the available patterns of Solidago spp. EOs profiles (Figure 2). PC3 explained only 6.8% of the total variance and had no significant effect on scores differentiation, so results were not presented. PC1 accounted for 29.4% of the total data set variance and showed high negative correlation with α-campholenal, trans-pinocarveol, cisand trans-verbenol, pinocarvone, myrtenal, myrtenol, verbenone, and trans-carveol and positive with carvone, germacrene D, γ-muurolene, γ-cadinene, spathulenol, and isospathulenol in inflorescences, and camphene, o-cymene, epoxyazulene, and (E)-nerolidol in leaves (Figure 2a). PC2 explained 23.6% of the total data set variance and was highly associated with positive loadings of αand β-copaene, α-cubebene, cubebol, α-muurolene, and δ-cadinene in inflorescences and with α-campholenal, trans-pinocarveol, cis-verbenol, and pinocarvone in leaves, as well as with negative loadings of bornyl acetate and isospathulenol in leaves.  Tables 1 and 2. Inflorescence variable numbers were marked in brown; leaf variable numbers were marked in green.
The PCA score plot model showed the arrangement of 40 EOs into two separate and two overlapping ellipses, each with a 95% confidence interval limit (Figure 2b). The group on the right-hand plot combined all S. gigantea EOs along the positive PC1. Variables with high PC1 loadings contributed the highest impact on the grouping of S. gigantea samples were germacrene D, γ-cadinene, γ-muurolene, and spathulenol in inflorescences, and camphene, o-cymene, epoxyazulene, (E)-nerolidol, and spathulenol in leaves. These compounds were shared among all S. gigantea EOs in the highest amounts compared to the other species studied. Conversely, variables with high PC2 loadings had a weak contribution on the grouping of EOs and were found in minor quantities in S. gigantea.
Solidago virgaurea EOs were clustered into a separate group on the upper positive side of the score plot, in distance from all other samples, indicating differences in their composition. The location of the samples can be explained by the same position of variables, which have a significant positive contribution to PC2. Variables with unit vectors close to each other were positively correlated, and their impact on the position of the samples was similar. Thus, α-and β-copaene, α-cubebene, cubebol, δ-cadinene, and α-muurolene in inflorescences, and α-campholenal, cis-verbenol, trans-verbenol, pinocarvone, and transpinocarveol in leaves, were common in S. virgaurea EOs and found in significantly higher amounts than in other species analyzed. Meanwhile, S. canadensis and S. × niederederi EOs were clustered into two partially overlapping ellipses, mainly in the left-hand score plot, showing the similarity of the volatile compound patterns. The arrangement of inflo- The PCA score plot model showed the arrangement of 40 EOs into two separate and two overlapping ellipses, each with a 95% confidence interval limit (Figure 2b). The group on the right-hand plot combined all S. gigantea EOs along the positive PC1. Variables with high PC1 loadings contributed the highest impact on the grouping of S. gigantea samples were germacrene D, γ-cadinene, γ-muurolene, and spathulenol in inflorescences, and camphene, o-cymene, epoxyazulene, (E)-nerolidol, and spathulenol in leaves. These compounds were shared among all S. gigantea EOs in the highest amounts compared to the other species studied. Conversely, variables with high PC2 loadings had a weak contribution on the grouping of EOs and were found in minor quantities in S. gigantea.
Solidago virgaurea EOs were clustered into a separate group on the upper positive side of the score plot, in distance from all other samples, indicating differences in their composition. The location of the samples can be explained by the same position of variables, which have a significant positive contribution to PC2. Variables with unit vectors close to each other were positively correlated, and their impact on the position of the samples was similar. Thus, α-and β-copaene, α-cubebene, cubebol, δ-cadinene, and α-muurolene in inflorescences, and α-campholenal, cis-verbenol, trans-verbenol, pinocarvone, and transpinocarveol in leaves, were common in S. virgaurea EOs and found in significantly higher amounts than in other species analyzed. Meanwhile, S. canadensis and S. × niederederi EOs were clustered into two partially overlapping ellipses, mainly in the left-hand score plot, showing the similarity of the volatile compound patterns. The arrangement of inflorescence EOs for both species coincided with a significant correlation of α-campholenal, transpinocarveol, pinocarvone, verbenone, cis-verbenol, trans-verbenol, myrtenal with PC1. In addition, the clustering of S. canadensis was influenced by myrtenol, trans-carveol, and carvone. Variables with high PC2 loading had no significant impact on EOs' arrangement, except for trans-verbenol in S. × niederederi leaves, showing similarity of this taxon to S. virgaurea. Meanwhile, a previous study of phenolic compounds showed greater chemical similarity of S. × niederederi to S. virgaurea than to S. canadensis [20]. Consequently, the phytochemical patterns complemented the evidence of S. × niederederi origin between native S. virgaurea and invasive S. canadensis, with the higher contribution of alien species than that of native ones. In addition, S. × niederederi EOs were much more scattered on the PCs space, indicating higher diversity than other species, suggesting that S. × niederederi is a continuously evolving taxon.

Discussion
Alien goldenrods are morphologically and phylogenetically close to each other, but differ in their ploidy level; S. gigantea is tetraploid (2n = 36), while S. canadensis, and S. × niederederi together with native S. virgaurea are diploids (2n = 18) [3,38]. The close relationships between Solidago spp. were reflected in the similarity of their phytochemical profiles. A comprehensive metabolomics approach to the different species indicated that the successful alien species had higher total number and more unique composition of secondary metabolites than their native congeners [39]. A comparison of the current and previous reports showed that our results are in agreement with previous reports for higher proportions of sesquiterpenoids in S. gigantea and hydrodrocarbons in S. virgaurea EOs [26,29,32,33]. Similar to our identification, the most abundant common compounds detected in the present study were also observed in previous studies on different species. Thus, the volatiles commonly found in various plant species have a high potential to accumulate in Solidago spp. EOs as well. On the other hand, there were compounds such as thymol and carvacrol, sporadically high levels of which were found in only a few S. × niederederi and S. canadensis EOs. Populations rich in thymol and carvacrol, compounds with a broad spectrum of biological activity [40][41][42], can be considered as a source of high potential raw material. In addition, Solidago spp. EOs differed significantly in some of the terpenes that can be considered as volatiles with great potential in chemophenetic studies of species. The first comparative study on terpenes as species distinctive compounds confirmed the origin of S. × niederederi as an interspecific taxon between S. canadensis and S. virgaurea. According to Orians [43], parental phytochemicals in hybrids tend to mostly express as either interme-diate or similar to one of the parent's compositions. The composition of S. × niederederi EOs was close to S. canadensis, one of parental species. The chromosome number may provide information about the hybrid origin of the species when it display allopolyploidy, but S. × niederederi exhibit a homoploid condition compared to its parental species and may backcross toward parental species. In this way, hybridization can increase the invasive capacity of goldenrods through gene introgression and significantly alter the ecosystems in which they grow [44]. However, species-specific compounds and chromosome number are not the main tools for hybrid identification. DNA fingerprinting techniques are the most reliable tools, but the use of additional phytochemical markers can provide insight into the ecological performance of hybrids and their further applications [45].
More often, the lower proportion of monoterpenes compared to sesquiterpenes [46] accounted for a higher proportion in the Solidago spp. EOs tested. Oxygenated monoterpenes have been proven to be the main phytotoxic active compounds in different plant EOs and have been highlighted as predictors of potential bioherbicides [47]. In addition, a tendency has been suggested that the monoterpene-rich EOs to promote higher phytotoxicity than sesquiterpene-rich EOs [48]. Comparative studies on the toxic activity of oxygenated monoterpenes revealed that the most active were alcohols, myrtenol and trans-pinocarveol, and ketones, verbenone and pinocarvone, which can be classified as predictors of the herbicidal activity of EOs [49][50][51]. Thus, S. canadensis and S. × niederederi, whose EOs differed from other species in some of oxygenated monoterpenes, suggested the potential of their raw materials for herbicidal activity. Meanwhile, terpene hydrocarbons have been found to be low phytotoxic EOs compounds [40,52]. On the other hand, Lawson et al. [36] reported that monoterpenes such as αand β-pinene, limonene, or myrcene showed weak antifungal activity. Recent findings reported that (E)-nerolidol and spathulenol, which were presented as species-specific sesquiterpenoids in S. gigantea EOs, revealed effective allelopathic and insecticidal effects with potential for the developing a new natural pesticide [53][54][55]. Many findings have demonstrated the biological activity of caryophyllene oxide [53,[55][56][57][58] that was a common compound in the presented Solidago spp. EOs. The potential biological activity of EOs is associated with the presence of high oxygenated compounds, as confirmed by a systematic review of phytotoxicity studies [59]. In this context, S. canadensis and S. × niederederi EOs, are of greatest interest in the development of new and safe bioproducts.
Considering the previous and presented results for Solidago spp. EOs, the prevalence of predominant and other compounds varied across different studies, and their comparison is not informative enough. Reports often provide single sample data that are difficult to summarize as a species-specific composition of volatiles. According to Zidorn [60], correct taxonomic identification, geographical location, plant harvesting season, plant parts, and other indirect factors are crucial in phytochemical studies, which often receive little attention. Exogenous factors or environmental regulated factors such as light, precipitation, growing site and soil are often considered to be the most important factors modifying the qualitative/quantitative composition of EOs [61]. Experiments have shown that plants exposed to drought stress increased the concentration of monoterpenes to protect plant cells from ROS damage [62]. Meanwhile, Caser et al. [63] found that drought increased the production of sesquiterpenoids and decreased monoterpenoids. According to Paulsen and Selmar [64], increased terpene synthesis is not supported by carbon allocation theory, but is attributed to changes in biomass production. In addition, differences in the chemical profiles of EOs are often explained in the context of the interaction of the metabolic inversion of the ratio of oxygen-free to oxygen-containing terpenes with the surrounding environment. Sesquiterpenes have been observed to be predominate during the dry season, while higher concentrations of sesquiterpenoids were found during the wet season [65]. The water scarcity increased the production of monoterpenoids and monoterpenes, while the opposite trend was observed for sesquiterpenes [66]. On the other hand, Tsusaka et al. [67], investigating the influence of genetic and environmental factors on sesquiterpenoids in Atractylodes lancea, (Thunb.) DC. found that the genotype had a greater effect on EO com-pounds than the conditions of the plant cultivation year. The volatiles were stable despite the changing growing conditions, but the absolute values of terpenoids were induced by the site of cultivation. Tardugno et al. [68] determined that the composition of Thymus vulgaris L. essential oils were highly influenced by the cultivating techniques. Over time, local environment leads to differences in metabolomics and the formation of ecotypes and chemotypes within a species [69].
The intraspecific differences in Solidago spp. volatile compounds observed in the present study can be explained by genotypic differences, as the plants grew under close conditions. Sexual reproduction helps maintain a high level of genetic and phenotypic diversity in goldenrods. Our previous study showed a high morphological diversity in Solidago spp. both between populations and between individual genets [3]. Similarly, high variability in volatiles and morphological characters was observed in wild Mentha longifolia L. accessions growing in the same field [70]. According to Zhao et al. [71], high genetic variations are characteristic within invasive and native areas of S. canadensis populations. Considerations suggested that the study of local populations makes it possible to identify intraspecific diversity that potentially reflects local genetic changes rather than the controversial dependence of terpenes synthesis under changing environmental conditions. The outstanding diversity in the goldenrods studied allows the selection of accessions in terms of the desired composition of EO volatiles. Phytochemical profiling of plant raw materials is an informative tool to learn about their potential for further development of new natural products.

Plant Material
Plant material of four Solidago spp. in the flowering phase was collected from six different sites in Vilnius district, Lithuania, in August 2018. Eighteen accessions of S. canadensis, seven of S. gigantea, nine of S. × niederederi, and three accessions of S. virgaurea were collected at least one kilometre apart from each other in abandoned dry grasslands and disturbed farmlands ( Table 3). The vegetation of the collection sites was characterized as semi-ruderal dry grassland dominated by plant communities of Agropyretea intermediirepentis and Artemisietea vulgaris. The EOs in habitats were sand or sandy loam, with low to moderate humus content (1.8-2.7%), rich in phosphorus (126-260 P 2 O 5 mg kg −1 ) and potassium (146-205 K 2 O mg kg −1 ), pH KCl varied from 5.8 to 7.1. The harvested plant material consisted of shoots of a single genet derived from a single seed. Individual genets were identified by phenological and morphological characteristics and by rhizome connections. The accessions of the same species were collected at least five meters apart from each other if more than one accession was collected from the same site. The plant material was dissected into inflorescences and leaves and dried separately at 25 • C. The botanical identification of species was based on morphological diagnostic characters such as the shape and size of inflorescences and ray flowers, stem color, and stem hairiness by Birutė Karpavičienė and Jolita Radušienė [3]. The specimens of evaluated Solidago spp. were deposited in the Herbarium of the Institute of Botany of Nature Research Centre (BILAS), Vilnius, Lithuania.

Isolation of Essential Oils
The plant material from 30 g of air-dried leaves and inflorescences was hydrodistilled separately for three hours using a Clevenger type apparatus. Each sample of yellowish EO was dried over anhydrous sodium sulphate and stored in a sealed vial at 4 • C until analysis. A sample preparation for chemical analysis included 1.0 µL of EO added into 1.0 mL of n-hexane following previous studies [23,33]. The essential oil content was calculated as relative percentages per 100 g of dry plant material.

Analysis of Essential Oils
The EOs analysis was performed using the GCMS-QP 2010 Ultra system equipped with a Shimazu autoinjector AOC-5000 (Shimadzu, Europa GmbH). A capillary column RXi-5MS (30 m × 0.25 mm i.d. × 0.25 film thickness µm) (Restek, Bellefonte, PA, USA) was used. The sample injection volume was 1 µL, a split ratio was 1:60 (v:v) and the split injector temperature was 260 • C. Helium was used as carrier gas with flow rate of 1.22 mL min −1 . The initial column temperature was 50 • C, held for 5 min and raised to 200 • C at the rate of 2 • C min −1 , then raised from 200 to 315 • C at the rate of 15 • C min −1 and held for 5 min. The detector ion source and interface temperatures were 200 • C and 280 • C, respectively. Mass spectra were acquired at an ionization voltage of 70 eV, a scan rate of 2500 m/z within the range of 29-500 m/z and a scan time of 0.2 s. The chromatographic analysis was run in triplicate.

Identification and Quantification of Components
The linear retention indices (LRI) of compounds were calculated using a homologous series of n-alkanes C8-C26 (Sigma-Aldrich, UK, purity >99.2%) injected at the beginning of the analysis and comparing the retention times of the eluted peaks with those of the alkanes [72]. Chromatographic data were analyzed using GC-MS solution software (Shimadzu, Europa GmbH). The EO constituents were identified by comparing the unique mass spectral fragmentation patterns of each peak with the mass spectral computer library database and those presented as standards in the NIST 14, FFNSC, WR10, and WR10R libraries, as well as comparing the obtained LRI with presented in NIST 14 datasets and reference [73] corresponding to the conditions for dimethylsilicone stationary phase with 5% phenyl groups. The relative percentages of analytes as the mean of the three runs were calculated from their peak areas in the chromatographic profiles without the use of correction factors corresponding to the conditions for the stationary phase of dimethylsilicone with 5% phenyl groups.
The repeatability and intermediate precision of analysis, expressed as relative standard deviation (RSD), was evaluated by performing the retention time and peak area values of five analytes in the same S. canadensis EO extract for intra-and inter-daily tests (Table 4). Table 4. Repeatability and intermediate precision on the relative retention time (t R ,) and peak area (A) of the five analytes in S. canadensis EOs expressed as relative standard deviation (RSD, %).

Analytes
Repeatability ( Repeatability was determined in five consecutive injections of EO in the same day. The RSD for the retention time ranged from 0.22 to 0.32% and for relative peak area from 0.33 to 0.45%. Intermediate precision was assessed by five injections over two different days, with RSD values ranging from 0.54 to 0.92% for retention time and 0.98 to 2.21% for peak areas. The accuracy of the quantification was satisfactory, as RSD value within the 3% range is generally considered as acceptable.

Data Analysis
Multivariate statistical analysis was performed using software package Statistica 10.0 (StatSoft Inc.). The Kruskal-Wallis ANOVA was used to determine the differences between species. Significant differences were specified by two-tailed test at p ≤ 0.05. Principal component analysis (PCA) was used to identify the similarities and differences between the EOs analyzed using statistically independent variables. The PCA was based on 24 inflorescence and 16-leaf standardized variables that different significantly between species and that represented constituents detected in 30% or more of EO samples in at least one species. Leaf and inflorescence data sets were pooled and used in PCA, resulting in more convincing results than separate leaf and inflorescence PCAs.

Conclusions
The frequency and percentage of distribution in the volatile constituents of Solidago spp. varied depending on the species, accessions, and plant parts. The principal compounds common to all Solidago spp. inflorescence and leaf EOs were α-and β-pinene, limonene, bornyl acetate, germacrene D, spathulenol, and caryophyllene oxide. Solidago spp. differed significantly in some of the distinctive terpenes that can be considered as compounds with high potential for chemophenetic and taxonomic studies of the genus. A comparison of volatile profiles for Solidago spp. confirmed the interspecific origin of S. × niederederi between S. canadensis and S. virgaurea with a higher metabolic contribution of alien species than native ones. The findings provide the bioprospecting of Solidago spp. as a source for specified composition of volatiles. The vast resources of invasive goldenrods are of great interest as a convenient and readily acceptable, underutilized source of natural bioactive compounds that can be used for different applications.
The combination of fingerprint and multivariate data analysis demonstrated a simplified assessment of the quality of wild plant materials. Correct species identification is essential for the development of raw material quality control protocols for the targeted collection and assessment of raw materials from wild populations. The screening of a relatively large number of plant accessions from mixed populations of different species allows for a more reasonable comparison of their volatile profiles and enables prediction of the most likely quality of raw materials harvested from the wild.