Wild-Grown and Cultivated Glechoma hederacea L.: Chemical Composition and Potential for Cultivation in Organic Farming Conditions

Glechoma hederacea L. is a medicinal plant that is known in traditional medicine for its anti-inflammatory, antibacterial, antiviral, and anticancer properties. This study evaluated the potential for commercial production of G. hederacea and compared the chemical composition and activity of 70% ethanol extracts and steam-distilled essential oils from wild-grown and cultivated G. hederacea collected in different harvesting periods. The main compounds identified in the 70% ethanol extracts were phenolic acids (chlorogenic and rosmarinic acids) and flavonoid O-glycosides. The essential oil varied in the three accessions in the range of 0.32–2.98 mL/kg−1 of dry weight. The extracts possessed potent antioxidant and anti-inflammatory properties in LPS-treated bone-marrow-derived macrophages. The results of flow cytometry show that extracts from different vegetation periods reduced the conversion of macrophages to the proinflammatory phenotype M1. The chemical composition varied the most with the different harvesting periods, and the most suitable periods were the flowering and vegetative phases for the polyphenolic compounds and essential oils, respectively. G. hederacea can be successfully grown under organic farming conditions, and cultivation does not significantly affect the chemical composition and biological activity compared to wild-grown plants.


Introduction
Glechoma hederacea L., or ground-ivy, is a medicinal plant that belongs to the Lamiaceae family, is native to Eurasia, and is widely distributed throughout the territory of Latvia in forests, shrubs, and roadsides [1]. The aerial parts of G. hederacea have been used in traditional medicine for centuries to treat colds, asthma, bronchitis, gastric diseases, diabetes, and inflammation [2,3]. G. hederacea provides anti-inflammatory, analgesic, antibacterial, antiviral, anticancer, diuretic, and antioxidant properties [4][5][6]. Its potential use in pharmacy applications [7] and as a food additive [8] has been studied. An in vivo study of G. hederacea hot water extract in a rat model indicated protective effects against cholestatic liver injuries [9]. The antioxidant activity of the hot water extract is significantly higher than that of vitamin C and Trolox [10]. The G. hederacea methanol extract exhibited an antibacterial inhibitory effect against 14 bacterial species [3]. The aqueous extract also reduced the expression of nuclear transcription factors that are related to inflammatory pathways and inflammatory cytokine production [9]. Both the ethanol and water extracts possessed

Seasonal Variation in G. hederacea during Field Trials
Meteorological conditions differed in the growing seasons of 2019, 2020, and 2021 (Priekuli meteorological station data). The total precipitation from the 1st of April until the 31st of October was highest in 2019 (528 mm), but was lowest in the same period in 2021 (472 mm). The sum of the active temperatures (>5 • C) was 1787, 1692, and 1760 • C in 2019, 2020, and 2021, respectively. The growing seasons of 2019 and 2020 were characterized by an optimal temperature and moisture regime in May, June, and July. This scenario resulted in the better development of G. hederacea, and plant regrowth allowed for three harvests in 2019 and four in 2020. Due to a very hot and dry June and July in 2021, plant regrowth was reduced, and the aerial part of the ground-ivy could only be harvested twice. Cultivated local accessions significantly varied in plant height and dry herb yield. In spring, G. hederacea develops flowering shoots. According to Slade and Hutching [21], at the flowering stage, plants grow vertically (up to 60 cm). In this study, the average plant height was 17.6 cm ( Table 1). The average height of the flowering shoots from the tested accessions varied between 13.5 cm (GH01) and 20.2 cm (GH03). The dry herb yield was highly correlated with the height of the flowering shoots; thus, the highest yield of 565.8 kg ha −1 was observed in accession GH03. After flowering, G. hederacea develops secondary aboveground stolons [22], and its aerial shoots only reach 10.1 cm in height on average. Due to the minimum clearance of the harvester being at 8 cm, none of the accessions have potential for commercial harvesting after the regrowth of vegetative shoots.
Under organic conditions, the ability to suppress weeds is an important trait because herbicides are prohibited from being used. Crop ground cover is negatively correlated with weed development; if the soil is covered with more crops, weed suppression is higher. The  Figure 1). Accession GH01 had the lowest plant height and yield and developed more slowly, which also contributed to its lower weed-suppression ability.  Under organic conditions, the ability to suppress weeds is an important trait because herbicides are prohibited from being used. Crop ground cover is negatively correlated with weed development; if the soil is covered with more crops, weed suppression is higher. The development intensity (soil coverage) of the tested local accessions significantly differed ( Figure 1). Accession GH01 had the lowest plant height and yield and developed more slowly, which also contributed to its lower weed-suppression ability.

Qualitative Analysis of G. hederacea 70% Ethanol Extracts
Qualitative analysis of the 70% ethanol extracts of G. hederacea was performed by UHPLC-HRMS/MS with an IT-TOF mass analyzer and a DAD detector, and the chromatograms are shown in Figure S1. Identification was based on high-resolution mass spectra, fragmentation, analysis of the available reference substances, and the literature data. The tentative identification was also confirmed by scanning the corresponding aglycone (quercetin, luteolin, kaempferol, and apigenin) masses in MRM mode, followed by a parent scan using a tandem mass spectrometer. The results of the qualitative analysis are summarized in Table 2 The main compounds identified in the 70% ethanol extracts of G. hederacea were phenolic acids (chlorogenic acid and rosmarinic acid) and flavonoid (quercetin, kaempferol, luteolin, and apigenin) O-glycosides, and a minor quantity of flavones (luteolin and apigenin) was also found. Several quercetin, luteolin, and apigenin glycosides in polar solvent extracts were previously discovered [15,18,22]. Small phenolic acids such as gallic, vanillic and ferulic acids described as being in G. hederacea [14,18] were not observed, nor were the alkaloids described by Kumarasamy et al. [23].
Three peaks (peaks 1-3) showed a signal at m/z 355 [M+H] + (C 16 H 18 O 9 ), corresponding to chlorogenic acid and its isomers. Peak 2 was identified as chlorogenic acid by comparing its retention time with that of a reference compound. However, the retention time of the chlorogenic acid isomers (peaks 1 and 3) did not correspond to the retention time of the available cryptochlorogenic acid reference compound. Peak 3 was tentatively identified as neochlorogenic acid, and the identification was based on the elution order of the chlorogenic acids in reversed-phase chromatography [24]. Peak 1 was tentatively identified as cis-chlorogenic acid. Chlorogenic acid was previously obtained in G. hederacea extracts [10,20]; in turn, chlorogenic acid isomers were reported for the first time.

Quantitative Analysis of the Main Phytocomponents in 70% Ethanol Extracts of G. hederacea
To determine the effect of the cultivation and harvesting time on the chemical composition of the plant, 10 major polyphenolic compounds (rosmarinic acid, chlorogenic acid, caffeic acid, rutin, apigenin 7-O-glucoside, luteolin 7-O-glucoside, apigenin, luteolin, hyperoside, and kaempferol 3-O-rutinoside) were selected, and their quantitative analysis was performed. The results are summarized in Tables 3 and 4. The major polyphenolic compound content in 70% ethanol extracts prepared from three accessions of wild-grown and cultivated plants harvested in the same growth phase (flowering shoot) and that in cultivated plants collected in different vegetation periods (May-September) were compared. Previously, Chou et al. [14] stated that rosmarinic acid, chlorogenic acid, caffeic acid, rutin, genistin, and ferulic acid were the most abundant phytochemicals in the hot water extract of G. hederacea and possessed potent antioxidant and anti-inflammatory properties. Our study is the first to provide data about quantitative analyses of apigenin 7-O-glucoside, luteolin 7-O-glucoside, apigenin, luteolin, hyperoside, and kaempferol 3-O-rutinoside in G. hederacea.

Comparison of the Content of the Main Phytocomponents in 70% Ethanol Extracts of G. hederacea Prepared from Wild-Grown and Cultivated Accessions
The content of 10 main polyphenolic compounds in 70% ethanol extracts prepared from the flowering shoots of wild and cultivated plant accessions was analyzed. The results obtained show that the contents of both kaempferol 3-O-rutinoside and rutin were significantly lower in the extracts prepared from two accessions (GH02 and GH03) of the cultivated plant. In contrast, the GH01 accession yielded a considerably higher content of these compounds in its extracts after the first year of cultivation. After the second year of cultivation, the kaempferol 3-O-rutinoside and rutin contents were slightly lower than those in the wild-grown plants. The content of other flavonoids in extracts prepared from wild-grown plants is similar to that yielded by cultivated accessions. However, after the second year of cultivation, a slight lowering in the concentration of the compounds was observed.
The content of two flavone derivatives, apigenin and luteolin, in extracts prepared from accession GH01 was higher after the first year of cultivation, while the other two accessions, GH02 and GH03, yielded concentrations similar to those obtained from the wild-grown plants. In the second year of cultivation, the apigenin and luteolin contents in all extracts were similar to those obtained from wild-grown plants.
The chlorogenic acid content in extracts prepared from cultivated flowering shoots in the first year was significantly higher than that in extracts prepared from wild-grown flowering shoots. However, in the second year of cultivation, the chlorogenic acid concentration in the extracts was similar to that of wild-grown plants. Although a similar trend was observed for rosmarinic acid, it was mainly due to the high content of this phenolic in extracts prepared from the GH01 accession after the first year of cultivation. In the second year of cultivation, the content of rosmarinic acid was similar in extracts prepared from all the accessions and reached the level of that in wild-grown plants. Cultivation did not affect the caffeic acid content in the extracts. The average content of the main phenolic compounds in wild-grown and cultivated plants is summarized in Tables 3 and 4.

Comparison of the Main Phytocomponents' Contents in the 70% Ethanol Extracts of G. hederacea Prepared from Cultivated Specimens That Were Harvested in Different Vegetation Periods
The content of flavone derivatives, apigenin, and luteolin was higher in the extracts prepared from flowering shoots. In extracts prepared from vegetative shoots harvested in July and September 2020, only trace amounts of apigenin were detected, and the luteolin level was below the detection limit in all samples except one. The extracts prepared from plants harvested in June 2020 showed a certain level of both compounds, although the level was considerably lower than in those prepared from flowering shoots. Neither apigenin nor luteolin were detected in extracts that were prepared from the vegetative shoots harvested in July 2021 (GH01, GH02, GH03), compared to the 97.0 µg/g and 39.5 µg/g of dry material from apigenin and luteolin, respectively, found in extracts prepared from flowering shoots harvested in June 2021. The content of both apigenin 7-O-glucoside and luteolin 7-O-glucoside in the extracts closely follows the pattern of apigenin and luteolin. Similarly, the concentration of kaempferol 3-O-rutinoside tended to be higher in extracts prepared from flowering shoots. However, the very low overall content of this compound in extracts prepared from cultivated accessions GH02 and GH03 makes the pattern less explicit compared to extracts prepared from accession GH01. The hyperoside content in extracts prepared from flowering shoots harvested in 2021 was significantly higher (p = 0.003) than that in vegetative shoots. In contrast, the results for extracts prepared from plant material harvested in 2020 scatter very much, so a statistically justified conclusion cannot be made ( Table 4).
The caffeic acid content in extracts prepared from the plants harvested in July 2020 was significantly higher (p = 0.02) than that in the extracts of plant material harvested in June 2020, despite both extracts being prepared from vegetative shoots. Similar contents of caffeic acid (294 ± 78 vs. 319 ± 50) were found in extracts prepared from material harvested in September. In contrast, the extracts prepared from the flowering shoots contained significantly less (p = 0.02) caffeic acid. Surprisingly, other phenolic acids showed the opposite tendency. Significantly (p = 0.009) less rosmarinic acid was found in the extracts prepared from vegetative shoots harvested in 2021 (3218 ± 592 µg/g) compared to those prepared from flowering shoots (5105 ± 362), with chlorogenic acid levels of 1461 ± 363 µg/g vs. 3390 ± 461 µg/g (p = 0.005), respectively (Table 3).

Quantitative Analysis of the Main Components of the Essential Oil of G. hederacea
As part of the secondary metabolism of plants, essential oils play an important role in a wide range of biological activities. Aromatic and volatile essential oil products can be used as natural additives to reduce oxidation and prevent inflammation, although their

Quantitative Analysis of the Main Components of the Essential Oil of G. hederacea
As part of the secondary metabolism of plants, essential oils play an important role in a wide range of biological activities. Aromatic and volatile essential oil products can be used as natural additives to reduce oxidation and prevent inflammation, although their chemical composition can be very variable due to their susceptibility to ecological, ontogenetic, climatic, postharvest, and intraspecies genetic factors [15]. The amount of essential oils and the composition of wild G. hederacea have been presented in several studies [15,[28][29][30]; however, no information on its commercial propagation was found. Previous research on the chemical composition of the essential oil of G. hederacea from Lithuania and Serbia [28][29][30] showed that terpenes represented the main portion of the oil. The major constituents were mostly sesquiterpene hydrocarbons and constituents with a germacrene skeleton. Germacrene D dominated all oils from the Vilnius district [30], while palmitic and linoleic acids were the main oils of the Serbian G. hederacea oil, along with germacrene D [29].
Hydrodistillation of G. hederacea yielded a pale yellow-green, pleasant-smelling essential oil with a dry weight concentration varying in the range of 0.32 to 2.98 mL·kg −1 (0.03-0.3%) for the plant material obtained in different vegetation periods (Table 5). Previous studies of plant material collected during the flowering period reported a lower yield of essential oils [28][29][30]. We observed that the potential to produce essential oils is slightly higher just after the flowering period, when the plant is busy spreading widely throughout its aboveground runners. The year of cultivation, the phenological stage, and the accession had a significant effect on the essential oil content; therefore, these factors must be taken into account regarding the quality of oil production.  The chemical compositions of the essential oils were determined according to their retention time and the spectrometric electronic library (NIST). The identities of the oil constituents were established using GC retention indices (RIs). In total, 64 compounds were identified in the essential oil of G. hederacea (Table S1). The main compounds identified in the oil were β-myrcene, β-ocimene, germacrene-D and germacrene-B, eucalyptol, and 1-octen-3-ol. The evaluation of the ten most dominant volatile compounds of the essential oils obtained from all accessions is summarized in Figure 3. According to the hierarchical clustering, the G. hederacea samples were clustered into three main groups. Group one contained relatively high amounts of β-myrcene, cis-β-ocimene, and 1-octen-3-ol. The samples in group one consisted of vegetative shoots of cultivated plants collected in 2019. Group two was constituted by high levels of germacene-D, germacene-B, and δ-guaiene essential oil compounds. Group three was rich in cumene, β-elemene, and bicyclogermacene. Groups two and three included both vegetative and flowering shoots over several seasons. Germacrene D was found to be one of the dominant components in G. hederacea samples analyzed in Lithuania, where germacrene D was detected at a rate of 20.7%. Our results show a similar trend in samples analyzed from North America (germacrene D, 19.4%) and Serbia (germacrene D, 7.3%), where one of the predominant chemical compounds is germacrene D [28,29]. Significant amounts of separated eucalyptol, also known as 1,8-cineole, should be noted among the dominant compounds. Eucalyptol is well known for its antiseptic and expectorant activity, as well as its antiviral activity [31]. Although the G. hederacea samples studied had a lower essential oil content than that of widely cultivated essential oil plants, such as chamomile, whose essential oil content exceeds 6 mL·kg −1 [32], peppermint at 8 mL·kg −1 [33], and salvia at 24 mL·kg −1 [34], it can be used as a potential source of niche essential oil rich in germacrene B, germacrene D, β-ocimene, eucalyptol, and β-myrcene. eucalyptol, also known as 1,8-cineole, should be noted among the dominant compounds. Eucalyptol is well known for its antiseptic and expectorant activity, as well as its antiviral activity [31]. Although the G. hederacea samples studied had a lower essential oil content than that of widely cultivated essential oil plants, such as chamomile, whose essential oil content exceeds 6 mL·kg −1 [32], peppermint at 8 mL·kg −1 [33], and salvia at 24 mL·kg −1 [34], it can be used as a potential source of niche essential oil rich in germacrene B, germacrene D, β-ocimene, eucalyptol, and β-myrcene.

Correlation of the Essential Oil Composition and Ethanol Extract Composition
The correlation of the variables can be determined based on their loading plots (Figure 4a). A small, large, or 90° angle implies positive correlation, negative correlation, or no correlation between two components, respectively. As indicated in Figure 4a, the contents of germacene D, germacene B, chlorogenic acid, rutin, bicyclogarmacene, δ-

Correlation of the Essential Oil Composition and Ethanol Extract Composition
The correlation of the variables can be determined based on their loading plots (Figure 4a). A small, large, or 90 • angle implies positive correlation, negative correlation, or no correlation between two components, respectively. As indicated in Figure 4a, the contents of germacene D, germacene B, chlorogenic acid, rutin, bicyclogarmacene, δ-guaiene, and kaempferol 3-O-rutinoside were correlated. This cluster of components had a weak or no correlation with another group of mutually correlated compounds-luteolin 7-O-glucoside, luteolin, apigenin 7-O-glucoside, apigenin, and cumene. In contrast to the flowering samples grouped more on the right side of the X axis in Figure 4c, both of these groups corresponded to the grouping of the vegetative samples on the left side. The flowering samples had more beta-myrcene, cis-beta-ocimene, and caffeic acid than other samples. The samples that were from different locations (Figure 4d) and were collected in different years (Figure 4b) formed overlapping clusters, indicating no differentiation due to these factors. The chemical composition is stable across different growing seasons and different vegetative shoots within seasons. Different chemotypes of the flowering and vegetative shoots can be used for different end-user purposes.
samples. The samples that were from different locations (Figure 4d) and were collected in different years (Figure 4b) formed overlapping clusters, indicating no differentiation due to these factors. The chemical composition is stable across different growing seasons and different vegetative shoots within seasons. Different chemotypes of the flowering and vegetative shoots can be used for different end-user purposes.

Total Content of Phenolic Compounds, Flavonoids, and DPPH Free Radical Scavenging Activity of G. hederacea Extracts
The capacity of the G. hederacea extracts to scavenge stable DPPH radicals, and the total phenolic and flavonoid contents in the extracts of the cultivated accessions, were higher after harvest in the later stages of vegetative growth (July-September) ( Table 6). Although the flowering phase is mentioned in the literature as the time for aerial part collection, scientific evidence to support this has never been published. Our results show that G. hederacea contains valuable substances and can also be harvested later in the growing season; vegetative shoots are also valuable and can be harvested several times from May to September. Our results are similar to those of the study by Varga et al. [20], in which the TPC and antioxidant capacity of the water extracts of various cultivated populations of G. hederacea in Hungary at different harvest times showed that the values

Total Content of Phenolic Compounds, Flavonoids, and DPPH Free Radical Scavenging Activity of G. hederacea Extracts
The capacity of the G. hederacea extracts to scavenge stable DPPH radicals, and the total phenolic and flavonoid contents in the extracts of the cultivated accessions, were higher after harvest in the later stages of vegetative growth (July-September) ( Table 6). Although the flowering phase is mentioned in the literature as the time for aerial part collection, scientific evidence to support this has never been published. Our results show that G. hederacea contains valuable substances and can also be harvested later in the growing season; vegetative shoots are also valuable and can be harvested several times from May to September. Our results are similar to those of the study by Varga et al. [20], in which the TPC and antioxidant capacity of the water extracts of various cultivated populations of G. hederacea in Hungary at different harvest times showed that the values of those collected during summer (July) were higher than those of early spring (April) or late fall (October). All extracts prepared from the cultivated samples showed higher activity than that of the wild samples. In the wild, accession GH03 showed the best results, but in the cultivation process in 2020 and 2021, the differences between the three accessions related to chemical composition were revealed, and GH03 showed slightly lower activity than GH01 or GH02. However, the antiradical activity and total phenolic and flavonoid contents in all cultivated samples during different vegetation periods were still slightly higher than those in wildgrown G. hederacea extracts, suggesting that the three accessions, from the perspective of chemical composition, have the potential for commercial cultivation.  The aerial parts of wild accessions (GH01, GH02, GH03) of G. hederacea originating from different regions of Latvia were analyzed. Data are represented as the mean ± SEM of three independent experiments carried out in duplicate. a The total phenolic content is expressed as the gallic acid equivalents per gram (mg GAE/g) of the lyophilized extract. b The total flavonoid content is expressed as the quercetin equivalents per gram (mg QE/g) of the lyophilized extract. The b IC 50 (µg/mL) value corresponds to the concentration that can scavenge 50% of DPPH free radicals.

The Effect of G. hederacea Extracts on Bone-Marrow-Derived Macrophage Polarization toward the M1 Phenotype
In the present study, we tested how G. hederacea extracts affect LPS-and IFN-γstimulated macrophage polarization of the proinflammatory phenotype M1, characterized by the presence of the surface markers CD80 and CD86. The number of M1-polarized macrophages increased four-fold compared to that in the untreated control ( Figure 5). Accession GH01 was selected for further experiments because it showed the most promising results after chemical composition analysis. The flow cytometry results indicate that extracts from different vegetation periods and at a concentration of 500 µg/mL reduced macrophage conversion of the proinflammatory phenotype M1 by 13-56%. The results of the study show that the G. hederacea extract generated from the July 2020 accessions had the strongest effect on the polarization of M1 macrophages ( Figure 5C). That the highest activity was seen in the extract from July 2020 correlates with the phytochemical composition data; this extract contained the highest amount of biologically active components. The ability to reduce the M1 macrophage population decreased as follows, in order from highest to lowest: July 2020 > September 2020 > May 2020 > June 2021 > May 2019 (wild-grown) > July 2021. In general, the G. hederacea extracts cultivated in 2020 showed the highest activity. According to the results of the MTT assay, the extract of G. hederacea was not toxic to BMDMs when applied for 24 h at concentrations ranging from 50 to 750 µg/mL ( Figure S2).
Rosmarinic acid and caffeic acid are phenolic acids found in a variety of plants, especially those of the Lamiaceae family [35]. Rosmarinic, caffeic, and chlorogenic acids were the main compounds identified in the extracts of G. hederacea; therefore, the activity of these substances in M1 macrophages was also examined. Rosmarinic acid and caffeic acid significantly decreased the LPS-induced expression of double-positive CD80 and CD86 cells ( Figure 5D). According to quantitative analyses, the rosmarinic acid content in the G. hederacea extracts was the highest among all the components. Furthermore, rosmarinic acid showed the highest ability to reduce the population of M1 macrophages compared to caffeic and chlorogenic acids. Rosmarinic acid at a concentration of 100 mM reduced the level of CD80-and CD86-positive cells by 29% compared with the LPS/IFN-γ control. Chou et al. [14] demonstrated that the hot water extract of G. hederacea inhibited NF-кB expression in LPS-stimulated RAW264.7 macrophages through the reduction of ROS levels and the downregulation of the expression of proinflammatory genes, which may be related to the presence of rosmarinic acid, chlorogenic acid, caffeic acid, rutin, genistin, and ferulic acid.  figure. (B,C) Expression of the proinfla cell surface markers CD80 and CD86 was analyzed by flow cytometry 24 h after treating with (GH01) extracts (250 µg/mL and 500 µg/mL) or (D) rosmarinic acid, chlorogenic acid, acid (each 100 mM) and LPS/IFN-γ. Data are represented as the mean ± SEM of three ind measurements from two parallel experiments. Differences between the measurements we using one-way ANOVA followed by Dunnett's test. * Significantly different from the L control (p < 0.05).
Rosmarinic acid and caffeic acid are phenolic acids found in a variety o especially those of the Lamiaceae family [35]. Rosmarinic, caffeic, and chlorogen were the main compounds identified in the extracts of G. hederacea; therefore, the of these substances in M1 macrophages was also examined. Rosmarinic acid an acid significantly decreased the LPS-induced expression of double-positive CD CD86 cells ( Figure 5D). According to quantitative analyses, the rosmarinic acid co the G. hederacea extracts was the highest among all the components. Furth rosmarinic acid showed the highest ability to reduce the population of M1 macr compared to caffeic and chlorogenic acids. Rosmarinic acid at a concentration of reduced the level of CD80-and CD86-positive cells by 29% compared with the L γ control. Chou et al. [14] demonstrated that the hot water extract of G. hederacea i NF-кB expression in LPS-stimulated RAW264.7 macrophages through the redu  figure. (B,C) Expression of the proinflammatory cell surface markers CD80 and CD86 was analyzed by flow cytometry 24 h after treating BMDMs with (GH01) extracts (250 µg/mL and 500 µg/mL) or (D) rosmarinic acid, chlorogenic acid, or caffeic acid (each 100 mM) and LPS/IFN-γ. Data are represented as the mean ± SEM of three independent measurements from two parallel experiments. Differences between the measurements were tested using one-way ANOVA followed by Dunnett's test. * Significantly different from the LPS/IFN-γ control (p < 0.05).

Plant Material
The aerial parts of three Glechoma hederacea habitats were collected for chemical analysis at the flowering stage during May 2019. The locations the wild accessions were obtained from are summarized in Table S2. The vouchers were deposited at the Institute for Environmental Solutions (IES) in Latvia under codes GH01, GH02, and GH03.

Field Trials and Preparation of G. hederacea Extracts
The field trials with G. hederacea were set up in two stages. First, randomly selected vegetative shoots were collected from three accessions in the wild in April 2019 and planted in an organically certified experimental field at the IES (57 • 19 11.7 N 25 • 19 18.8 E, 115 m altitude). The plot size was 6 m 2 (one replication), and the plant spacing was 30 × 75 cm. The aerial parts (vegetative shoots) were collected for chemical analysis after the plants were regrown on 1 July 2019 and 23 August 2019 (hereafter referred to as "cultivated") by cutting shoots approximately 8 cm above the soil surface in all plots and taking an average sample of 150 g for chemical analysis. In 2020, aerial parts from the same plots were collected four times-on the 18 May (flowering shoots), 1 June, 21 July and 15 September (vegetative shoots).
Vegetative shoots were collected from the same plots on 23 July 2020, and a new trial was carried out with the aim of evaluating the suitability of accessions for commercial cultivation. The field assays were set up as a randomized complete block design in four replications. The plot size was 16.2 m 2 (6 rows with 75 cm between rows). The shoots were cut into 5-15 long pieces before being embedded in the soil in rows.
In 2021, the plants' regrowth was slower, and the aerial parts were collected on the 1st of June (flowering branches) and the 23 July (vegetative branches). The aerial parts of each plot were cut approximately 8 cm high, simulating mechanical harvesting, and were collected and weighed, and the yield of the fresh part was calculated. The aerial parts collected from G. hederacea were air dried at room temperature. After drying the weighed herbs, the yield of the dry herbs was calculated (with 12% moisture). Additionally, the percentage of soil covering the accessions was measured every 2-3 weeks during the vegetation seasons of 2020 and 2021 ( Figure 6). A quadrocopter platform with vertical take-off and landing (VTOL) was used to collect a set of aerial images at a fixed flight altitude of 12 m over the experimental crop fields. The RGB sensor was mounted during the flights to obtain high-resolution image data. An automated Python code was developed to process the collected data. Soil coverage (%) by green plant segments was calculated for each separate plot and used for statistical analysis and visualization.
analysis at the flowering stage during May 2019. The locations the wild accessions were obtained from are summarized in Table S2. The vouchers were deposited at the Institute for Environmental Solutions (IES) in Latvia under codes GH01, GH02, and GH03.

Field Trials and Preparation of G. hederacea Extracts
The field trials with G. hederacea were set up in two stages. First, randomly selected vegetative shoots were collected from three accessions in the wild in April 2019 and planted in an organically certified experimental field at the IES (57°19′11.7′′ N 25°19′18.8′′ E, 115 m altitude). The plot size was 6 m 2 (one replication), and the plant spacing was 30 × 75 cm. The aerial parts (vegetative shoots) were collected for chemical analysis after the plants were regrown on 1 July 2019 and 23 August 2019 (hereafter referred to as "cultivated") by cutting shoots approximately 8 cm above the soil surface in all plots and taking an average sample of 150 g for chemical analysis. In 2020, aerial parts from the same plots were collected four times-on the 18 May (flowering shoots), 1 June, 21 July and 15 September (vegetative shoots).
Vegetative shoots were collected from the same plots on 23 July 2020, and a new trial was carried out with the aim of evaluating the suitability of accessions for commercial cultivation. The field assays were set up as a randomized complete block design in four replications. The plot size was 16.2 m 2 (6 rows with 75 cm between rows). The shoots were cut into 5-15 long pieces before being embedded in the soil in rows.
In 2021, the plants' regrowth was slower, and the aerial parts were collected on the 1st of June (flowering branches) and the 23 July (vegetative branches). The aerial parts of each plot were cut approximately 8 cm high, simulating mechanical harvesting, and were collected and weighed, and the yield of the fresh part was calculated. The aerial parts collected from G. hederacea were air dried at room temperature. After drying the weighed herbs, the yield of the dry herbs was calculated (with 12% moisture). Additionally, the percentage of soil covering the accessions was measured every 2-3 weeks during the vegetation seasons of 2020 and 2021 ( Figure 6). A quadrocopter platform with vertical take-off and landing (VTOL) was used to collect a set of aerial images at a fixed flight altitude of 12 m over the experimental crop fields. The RGB sensor was mounted during the flights to obtain high-resolution image data. An automated Python code was developed to process the collected data. Soil coverage (%) by green plant segments was calculated for each separate plot and used for statistical analysis and visualization.  The dried herbs from all replications were mixed, and an average sample of 150 g was taken for chemical analysis. Powdered and dried samples of G. hederacea were macerated with 70% ethanol solution in water at 1:10 w/v. The prepared solutions were incubated for 7 days in a dark, cool place and were frequently shaken. Afterward, the material was pressed, and the remaining solid was squeezed to remove all the remaining solvent. The obtained solutions were clarified by decantation and centrifugation. Subsequently, the extracts were concentrated using a rotary evaporator. Finally, the solutions underwent lyophilization. The powder was labeled and stored in a refrigerator at −20 • C prior to further analysis.
In vitro and ex vivo experiments with G. hederacea extracts were carried out using lyophilized plant material, which was dissolved in distilled water.
The following mass spectrometer operating parameters were utilized: electrospray ionization in positive and negative ionization modes, mass scan range (m/z)-from 120 to 1000; detector voltage-1.5 kV; nebulizing gas (N 2 ) flow-1.5 mL/min; ion accumulation time-100 ms; CDL temperature-250 • C; collision gas-argon and collision energy was set at 50%. LC-MS data were processed by LabSolutions software. A diode array detector was used to record UV-Vis spectra over a range from 210 nm to 800 nm.
70% ethanol extracts of G. hederacea were injected into the chromatographic system without preliminary processing.

UHPLC-MS/MS Analysis of Flavonoids
Quantitative analysis of flavones and flavonoid glycosides was performed using a Xevo TQ-S micro (Waters) tandem mass spectrometer that operated in the positive electrospray ionization and multiple reaction monitoring (MRM) modes. The MRM parameters of each compound were optimized by infusion in a mass spectrometer and are detailed in Table S3. Chromatographic separation was achieved on an Acquity BEH C18 column (2.1 × 100 mm, 1.7 µm, Waters) at a constant temperature of 30 • C using a Waters Acquity UPLC system with a gradient elution with mobile phases A (0.1% formic acid in water) and B (acetonitrile) at a flow rate of 0.4 mL/min. The gradient parameters were as follows: 5% B, 0.5 min-5% B; 8 min-98% B; 10 min-98% B; 11 min-5% B; 12 min-5% B. The autosampler was set at 10 • C, and the sample injection volume was 1 µL.
Before analysis, samples were diluted 10 or 100 times with reserpine solution (internal standard) (10 ng/mL) in 70% ethanol; for each sample, two measurements were taken. The concentration of individual components in the extracts was determined by a calibration curve (ranging from 50 ng/mL to 10 µg/mL for all analytes).

UHPLC-MS/MS Analysis of Phenolic Acids
Quantitative analysis of phenolic acids was performed using a Quattro Micro (Waters) tandem mass spectrometer that operated in the negative electrospray ionization mode. Multiple reaction monitoring (MRM) parameters are detailed in Table S4. Chromatographic separation was achieved on an Acquity BEH C18 column (2.1 × 50 mm, 1.7 µm, Waters) at a constant temperature of 30 • C using a Waters Acquity UPLC system with a gradient elution with mobile phases A (0.1% formic acid in water) and B (acetonitrile) at a flow rate of 0.25 mL/min. The gradient parameters were as follows: 5% B, 2.5 min-98% B; 4.5 min-98% B; 4.7 min-5% B; 6 min-5% B. The autosampler was set at 10 • C, and the sample injection volume was 5 µL. Samples before analyses were diluted 10 or 200 times with 70% ethanol; for each sample, two measurements were taken. The calibration concentrations ranged from 50 ng/mL to 5 µg/mL for all analytes.

Determination of Essential Oil Content
Plant samples dried at room temperature (~25 • C) were used for analysis. Essential oils were prepared using a Clevenger-type hydrodistillation apparatus. Fifteen grams of powdered herbal drug were transferred to a 500 mL flask, distilled water was added as the distillation liquid, and 0.50 mL of xylene was added to a graduated tube. The distillation was carried out at a rate of 3-4 mL·min −1 for 3 h. The essential oil yield (mL·kg −1 ) was calculated based on the dried weight of the samples. The oil dissolved in the organic layer was separated and dried over anhydrous sodium sulfate to eliminate moisture. The samples were preserved in a sealed amber glass vial at 4 • C until GC-MS analysis. The samples were replicated 3 times.

GC-MS Analysis of Essential Oils
For the test solution, 100 µL of the essential oil sample was diluted with 900 µL of cyclohexane and was mixed. Analyses were performed on an Agilent Technologies 7820A gas chromatograph coupled to Agilent 5977B mass selective detector (MSD) equipment. A nonpolar HP-5 capillary column (60 m × 0.25 mm, 0.25 µm film thicknesses) was coated with 5% phenyl and 95% methyl polysiloxane. The carrier gas was helium (He) with a split ratio of 1:100, and a flow rate of 1.3 mL/min was applied. The volume of injection was 3 µL. The temperature program began at 70 • C for 3 min and then increased at a rate of 10 • C/min to 290 • C; finally, 290 • C was maintained for 10 min. The injector temperature was 270 • C. Mass spectra were recorded at 70 eV. The mass range was from m/z 30-550. The ion source temperature was maintained at 230 • C. The components were identified based on their retention indices (which were determined with reference to a homologous series of C5-C24 n-alkanes) by comparing their mass spectra with those stored in the National Institute of Standards and Technology (NIST) MS search 2.2 library. Agilent MassHunter Qualitative Analysis 10.0 data acquisition software was applied to analyze the GC-MS data. The amount of separated compounds was calculated in the peak areas using the normalization method without correction factors.

Determination of the Total Phenolic Content
The total phenolic content (TPC) in the G. hederacea extracts was determined using the Folin-Ciocalteu colorimetric method described by Kähkönen et al. [36], with slight modifications. In brief, 20 µL of extract was added to a 96-well plate and mixed with 100 µL of 10% Folin-Ciocalteu reagent, followed by the addition of 80 µL 7.5% Na 2 CO 3 solution. After incubation at room temperature for 30 min in the dark with slight shaking, the absorbance at 765 nm was measured on a Hidex Sense microplate reader. Gallic acid was used as a standard for the calibration curve. The total phenolic content was expressed as mg of gallic acid equivalent (GAE) per g of lyophilized extract. All measurements were made in triplicate.

Determination of the Total Flavonoid Content
The total flavonoid content (TFC) was determined using the colorimetric method described by Wang et al. [37]. Briefly, 100 µL of sample solution was mixed with the same volume of 2% aluminum trichloride in methanol. Similarly, a blank was prepared by adding 100 µL of sample solution to methanol without AlCl 3 . After incubation at room temperature for 10 min, the absorbance at 415 nm was measured on a Hidex Sense microplate reader. The calibration curve was prepared using various concentrations of quercetin (0-250 µg/mL) dissolved in methanol. TFC was expressed as mg of quercetin equivalent (QE) per g of lyophilized extract. All measurements were made in triplicate.

2,2-Diphenyl-1-Picrylhydrazyl (DPPH) Free Radical Scavenging Assay
DPPH was used to assess the free radical scavenging (antioxidant) properties of the G. hederacea extracts. DPPH radical scavenging activity was measured according to Brand-Williams et al. [38], with some modifications. For the assay, 20 µL of extract diluted in water was mixed with 180 µL of DPPH in methanol (40 µg/mL) in the wells of a 96-well plate. The plate was kept in the dark at room temperature for 15 min. Decreases in absorbance at 517 nm were measured using a Hidex Sense microplate reader. Ascorbic acid solutions in the concentration range of 0-800 µg/mL were used as a standard. The extract was tested in a range of concentrations to establish the IC50 (the concentration that reduced the absorbance of DPPH by 50%). The radical scavenging activity was calculated using the following formula: DPPH radical scavenging activity % = [(A 0 -A 1 )/A 0 ] × 100, where A 0 is the absorbance of the control and A 1 is the absorbance of the sample.

Isolation of Bone-Marrow-Derived Macrophages
Bone-marrow-derived macrophages (BMDMs) were isolated from male C57BL6/J inbred mice (18-20 weeks old, Envigo, Netherlands). BMDMs were cultured for 6-7 days in RPMI-1640 medium with GlutaMAX supplemented with 10% FBS, 1% antibiotics (100 U/mL penicillin and 100 µg/mL streptomycin), and 10 ng/mL M-CSF. The experimental procedures were carried out in accordance with the guidelines of the European Community (2010/63/EU) and local laws and policies, and were approved by the Latvian Animal Protection Ethical Committee, Food and Veterinary Service, Riga, Latvia.
The Petri dish containing BMDMs was washed twice with HBSS. The cells were detached with 0.5% trypsin and placed in RPMI-1640 medium with 10% FBS and 1% antibiotics, and the cell suspension was centrifuged at 300× g at room temperature for 5 min. The cells were resuspended in RPMI-1640 medium supplemented with 10% FBS and 1% antibiotic, and seeded in 12-or 96-well plates. The cells were then incubated in an incubator at 37 • C for at least 1 h prior to the experiment.

Evaluation of Cell Viability with the MTT Assay
The viability of BMDMs after 24 h of incubation with different concentrations of G. hederacea extract was determined using the MTT assay. BMDMs were seeded in 96-well plates at a final concentration of 50 × 10 4 cells/mL. After incubation with the extract, the cells were further incubated with MTT solution (1 mg/mL) at 37 • C for 1-2 h. Subsequently, the medium was aspirated, and isopropanol was added to each well to dissolve the formazan crystals that formed during the incubation period. The absorbance was determined spectrophotometrically at 570 nm using a reference wavelength of 650 nm on a Hidex Sense microplate reader.

Treatment of Bone-Marrow-Derived Macrophages with Extract, Polarization toward the M1 Phenotype, and Analysis by Flow Cytometry
BMDMs were seeded at a density of 30 × 10 4 cells/well in 1 mL of culture medium. The cells were stimulated with 5 ng/mL LPS and 10 U/mL murine IFN-γ for macrophage polarization toward the M1 (proinflammatory) phenotype, together with GH01 extracts (250 µg/mL and 500 µg/mL) or 100 µM of rosmarinic acid, chlorogenic acid, and caffeic acid for 24 h.
The cells were then washed twice with HBSS and harvested with 0.5% trypsin. Then, an RPMI medium with 10% FBS and 1% antibiotic was added, and the cell suspension was centrifuged at 300× g for 5 min. The cells were then incubated with specific conjugated antibody mixtures (at a dilution of 1:200) for 45 min on ice in the dark. For M1-polarized cells, the mixture contained the following monoclonal antibodies: FITC-conjugated antimouse F4/80, phycoerythrin (PE)-conjugated anti-mouse CD86, and biotin-conjugated anti-mouse CD80. After staining, the expression of the markers was analyzed by flow cytometry (BD FACSMelody TM , BD Biosciences, San Jose, CA, USA).

Statistical Analysis
Three biological replicates per sample were used for experiments with BMDMs, and two technical replicates were analyzed per biological replicate.
Quantitative results are expressed as the mean ± standard error of the mean (SEM) of three independent experiments and were analyzed using the computer software Graph-Pad Prism 8. Statistical analyses were conducted using one-way analysis of variance (ANOVA), followed by the Dunnett test. An unpaired t-test was used for the MTT assay. The results of the MTT test are expressed as the mean ± standard deviation (SD). Values of p < 0.05 were considered significant. A heatmap with scaled data was created in the R package pheatmap [39]. Principal component analysis (PCA) was performed using the FactoMineR [40] and factoextra [41] packages in R software version 4.1.4.

Conclusions
No detailed guidance on the cultivation of G. hederacea has been published thus far. Several factors, such as genotype, harvesting time, growing location, and other conditions, could affect the chemical composition of the plant and its biological activity [20]. Our study on G. hederacea demonstrated differences between the accessions tested. The wild populations of G. hederacea under field conditions differed significantly in their plant height, yield, and development. G. hederacea could be harvested several times during the season. The data obtained clearly revealed the presence of chemical variety within the accessions of G. hederacea and indicated the chemical compositions were different between harvesting periods. From the point of view of commercial cultivation and based on the different chemical compositions of the flowering and vegetative shoots of G. hederacea, the most suitable harvesting time for polyphenolic compounds is the flowering phase, while for essential oils, it is the vegetative phase. To harvest G. hederacea parts that are rich with caffeic acid, the vegetative phase is more appropriate. The changes in the content of flavonoids in G. hederacea plant extracts are more pronounced with changes in the time of plant collection rather than in the phenolic acid (chlorogenic acid, rosmarinic acid, and caffeic acid) content. G. hederacea can be successfully grown under organic farming conditions, and cultivation does not significantly affect its chemical composition or biological activity compared to wildgrown plants. Accession GH03 was the most suitable sample for commercial cultivation within the tested samples and cultivation seasons.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/plants11060819/s1. Figure S1: Representative chromatograms of G. hederacea 70% ethanol extract; Figure S2: Effects of G. hederacea extract on bone-marrow-derived macrophage viability measured by the MTT assay; Table S1: Volatile compound composition (%) of the essential oil of dry plants of G. hederacea extract harvested during 2019-2021; Table S2: Locations of the collected wild G. hederacea accessions; Table S3: MRM parameters applied for the analysis of flavonoids in G. hederacea extract; Table S4: MRM parameters applied for the analysis of phenolic acids in G. hederacea extract.