Phytochemical Composition and Antioxidant, Anti-Acetylcholinesterase, and Anti-α-Glucosidase Activity of Thymus carnosus Extracts: A Three-Year Study on the Impact of Annual Variation and Geographic Location

Thymus carnosus Boiss. is a near-threatened species, and, as for many species, its potential for medicinal purposes may be lost if measures towards plant protection are not taken. A way of preserving these species is to increase knowledge about their medicinal properties and economic potential. Thus, with the objective of studying the potentiality of introducing T. carnosus as a crop, the stability of the phytochemical profile of T. carnosus was studied during a period of three years by comparing the phytochemical profile of extracts obtained from plants harvested in two different edaphoclimatic locations, as well as by comparing the respective bioactivities, namely, antioxidant, antidiabetic, antiaging, and neuroprotective activities. It was reported, for the first time, the effect of annual variation and geographic location in the phytochemical composition of aqueous decoction and hydroethanolic extracts of T. carnosus. In addition, the presence of two salvianolic acid B/E isomers in T. carnosus extracts is here described for the first time. Despite the variations in phytochemical composition, according to harvesting location or year, T. carnosus extracts maintain high antioxidant activity, assessed by their capacity to scavenge ABTS•+, •OH , NO•, O2•− radicals, as well as to prevent β-carotene bleaching. All extracts presented significant potential to inhibit acetylcholinesterase (AChE), tyrosinase, and α-glucosidase, denoting neuroprotective, anti-aging, and anti-diabetic potential. In conclusion, the vegetative stage and location of harvest are key factors to obtain the maximum potential of this species, namely, a phytochemical profile with health benefit bioactivities.


Introduction
Nowadays, the food, pharmaceutical, and cosmetic industries have a growing and constant demand for new products and ingredients with health promoting effects, with an emphasis on natural products. In addition to the general low cost and safe use, the consumption of natural products attracts consumers [1]. Medicinal and aromatic plants are dependent on other factors, such as the geographic location and edaphoclimatic conditions, which may induce changes in the phytochemical composition.
In the Thymus genus, most studies analysing the effects of edaphoclimatic conditions, inter-year climate changes, and vegetative stage are mainly towards the composition of essential oils, being reported that these factors modulate the essential oils' phytochemical composition, and therefore their bioactivities, as described for Thymus vulgaris L. [16], Thymus pulegioides L. [17], Thymus pallescens Noë. [18], or Thymus hyemalis Lange [19]. For Thymus vulgaris [16], in addition to the phytochemical composition, a change in antioxidant and antibacterial activities was also observed. Regarding thyme extracts, using Thymus longicaulis C. Presl collected in different seasons, it was observed that its hydro-methanolic (1:1; % v/v) extract's polyphenolic composition presented significant variances [20]. The main compound, RA, varied from 12.97 to 3029.56 µg/mL (in quercetin equivalents) in plants harvested from July to October of the same year. However, other phenolic acids also showed variations with the time of harvest and, for example, salvianolic acid K showed a concentration variation inverse to that of RA, and this was also observed for other phenolic acids. The seasonal variance effect was also clearly observed in the tested bioactivities, with significant differences for anti-inflammatory, anti-proliferative, and antioxidant activities [20]. The variation of pentacyclic triterpenoids concentration, such as OA and UA (compounds present in high amounts in T. carnosus HE extract [9]), induced by the vegetative stage, has been described in various thyme species. Analysed in methanolic extracts, small variations in OA and UA content through the various vegetative phases was reported in Thymus praecox ssp. arcticus Opiz. extracts, while, in other species, such as Thymus pulegioides at the end of vegetative stage OA and UA content, was 2.2 and 2.98 times higher than the value at fruit maturation stage [21]. Interestingly, the stage with higher content in each pentacyclic triterpenoids is dependent on the species under study [21], highlighting the need to understand the best conditions for each species.
Therefore, considering T. carnosus as a potential crop, the aim of this work was to study the stability of T. carnosus phytochemical profile over a period of three years to compare the phytochemical profile of extracts obtained from plants harvested in two different edaphoclimatic conditions, as well as to compare the respective bioactivities, namely, antioxidant, anti-diabetic, anti-aging, and neuroprotective activities. For this reason, in the present study, aerial parts of T. carnosus were collected in November, a post-flowering stage, in which the phytochemical profile of the plant may reflect the environmental stresses experienced in the previous months. Aerial parts were collected both at its natural habitat, where it grows as an endemic wild plant, as well as at UTAD's botanical garden, where the plant has adapted to a different climatic condition. The first location generally presents higher average temperatures in the last trimester of the year when compared to UTAD's botanical garden. In addition, precipitation is usually higher in northern Portugal, which may also induce variations in the phytochemical profile. The composition of aqueous and hydroethanolic extracts in phenolic and terpenoid compounds was assessed by chromatographic methodologies and correlated to its potential as an antioxidant, anti-diabetic, anti-aging, and/or neuroprotective agent.

Preparation of Extracts
Aqueous and hydroethanolic extracts were obtained as described by Martins-Gomes et al. (2018) [9], using aqueous decoction extraction (AD) and exhaustive hydroethanolic (HE) extraction procedures, respectively. Briefly, 0.5 g of lyophilized, ground plant material were used for both methods. To obtained AD extracts, 150 mL of distilled water were added to the plant material, followed by heating to 100 • C, where it was maintained for 20 min, under agitation. After this period, the mixture was allowed to cool down, to room temperature, and then was filtered. HE exhaustive extraction comprised a three-step sequential extraction method of the plant material with 50 mL of an ethanol:water solution (80:20, % v/v), each of the steps being under agitation (orbital shaker; 150 rpm) for one hour and then centrifuged (7000 rpm, Sigma Centrifuges 3-30 K, St. Louis, MO, USA). The three supernatants were collected, combined, and filtered. Both extracts were filtered twice (Whatman nº 4 filter and fiberglass filter (1.2 µm; acquired from VWR International Ltd., Alfragide, Portugal)) and concentrated to 100 mL in a rotary evaporator (35 • C), the step in which the ethanol was removed from the HE extract [9]. These methodologies were repeated three times, and all extracts were frozen and lyophilized, followed by weighing for yield calculation and proper storing until further analysis.

Total Phenolic Compounds, Total Flavonoids and Ortho-Diphenols Content
Total phenolic compounds content (TPC), total flavonoid content (TFC), and orthodiphenol content (ODC) were quantified using colorimetric reactions based on Folin-Ciocalteau reagent, molybdenum complexation, and aluminium complexation, respectively. All methodologies were performed as described by Taghouti, et al. 2020 [23]. TPC and ODC were expressed as caffeic acid equivalents (mg CA eq./g lyophilized plant or mg Figure 1. Identification of harvesting locations and climatic conditions in which the Thymus carnosus used in this work grew. Geographic location of T. carnosus harvest points, elevation from sea level and distance to coastline, nearest location in straight-line distance (A), and annual (B) and mensal (C) variation of average temperature (in °C; yellow dots) and accumulated precipitation for 2018-2020 period in L1 (blue bars) and L2 (green bars).

Preparation of Extracts
Aqueous and hydroethanolic extracts were obtained as described by Martins-Gomes et al. (2018) [9], using aqueous decoction extraction (AD) and exhaustive hydroethanolic (HE) extraction procedures, respectively. Briefly, 0.5 g of lyophilized, ground plant material were used for both methods. To obtained AD extracts, 150 mL of distilled water were added to the plant material, followed by heating to 100 °C, where it was maintained for 20 min, under agitation. After this period, the mixture was allowed to cool down, to room temperature, and then was filtered. HE exhaustive extraction comprised a three-step sequential extraction method of the plant material with 50 mL of an ethanol:water solution (80:20, % v/v), each of the steps being under agitation (orbital shaker; 150 rpm) for one hour and then centrifuged (7000 rpm, Sigma Centrifuges 3-30 K, St. Louis, MO, USA). The three supernatants were collected, combined, and filtered. Both extracts were filtered twice (Whatman nº 4 filter and fiberglass filter (1.2 µ m; acquired from VWR International Ltd., Alfragide, Portugal)) and concentrated to 100 mL in a rotary evaporator (35 °C), the step in which the ethanol was removed from the HE extract [9]. These methodologies were repeated three times, and all extracts were frozen and lyophilized, followed by weighing for yield calculation and proper storing until further analysis.

Total Phenolic Compounds, Total Flavonoids and Ortho-Diphenols Content
Total phenolic compounds content (TPC), total flavonoid content (TFC), and orthodiphenol content (ODC) were quantified using colorimetric reactions based on Folin-Ciocalteau reagent, molybdenum complexation, and aluminium complexation, respectively. All methodologies were performed as described by Taghouti, et al. 2020 [23]. TPC and ODC were expressed as caffeic acid equivalents (mg CA eq./g lyophilized plant or mg CA eq./g extract), and TFC was expressed as catechin equivalents (mg C eq/g lyophilized plant or mg C eq/g extract).

Phytochemical Composition Profiling and Quantification by HPLC-DAD and HPLC-ESI-MS n
Individual phenolic compounds, oleanolic acid, and ursolic acid identification and quantification were performed by RP-HPLC-DAD analysis using a Vanquish Core HPLC system (Thermo Fisher Scientific, Waltham, MA, USA) equipped with auto-sampler, pump, column compartment, and diode array detector. Chromatographic separation was performed using a C18 column (Merck Purospher ® STAR, Hibar ® C18; 250 mm × 4.6 mm; particle size 5 µm), with an injection volume of 100 µL, and the temperature kept at 40 • C, and the flow rate was 0.5 mL/min.
UV/Vis detection was performed at 200-600 nm, being 280 nm and 325 nm, which were used for phenolic compound quantification, and 210 nm was used for terpenoid quantification. Chromeleon software (Version 7.3; Dionex, USA) was used for data acquisition, peak integration, and analysis.
RP-HPLC-ESI-MS n analysis was performed for accurate phenolic compounds identification using a Thermo Scientific system equipped with a Finnigan Surveyor Plus auto-sampler, pump, LXQ Linear ion trap detector, and a photodiode array detector. The elution system, column, solvents, temperature, flow rate, injection volume, and detection parameters were performed as described by Martins-Gomes et al. (2018) [9].
The identification of individual phenolic compounds present in T. carnosus extracts was based on the data acquired from UV-VIS and mass spectrometry analysis, as well as retention time comparison with commercial standards and/or literature data. Oleanolic and ursolic components were identified only by HPLC-DAD by comparison to their respective commercial standards. Phytochemicals' quantification was performed based on calibration curves of commercial standards, if available, or using the aglycones or standard compounds with structural similarity to commercial standards. Caffeic acid (CA; PubChem CID: 689043) was quantified as its respective standard. Luteolin and apigenin derivatives were quantified as luteolin-7-O-glucoside (L-7-G; PubChem CID: 5280637); quercetin derivatives were quantified as quercetin-3-O-glucoside (Pubchem CID 25203368); eriodyctiol derivatives were quantified as eriodyctiol-7-O-glucoside (Pubchem CID 13254473); RA, SAA iso, SAK, and salvianolic acid K isomer (SAK iso) were quantified as RA (PubChem CID: 5281792).

In Vitro Antioxidant Activity Assessment
In vitro radical scavenging capacity of T. carnosus aqueous and hydroethanolic extracts was evaluated using ABTS (ABTS •+ ) and superoxide (O 2 •− ) radicals scavenging and βcarotene bleaching assays. ABTS •+ scavenging assay was performed, as described by Taghouti et al. (2018) [24] and expressed as mmol Trolox equivalent/g dry plant. Trolox was also used as a positive control (IC 50 = 0.24±0.01 mg/mL) Regarding O 2 •− scavenging by T. carnosus extracts, 6.7 µL of extracts were added to 193.3 µL of the reaction mixture (174 µL of phosphate buffer (50 mM; pH 8), 12.86 µL of nitro blue tetrazolium solution (NBT; 4 mM) and 6.43 µL of hypoxanthine solution (4 mM)), and the mixture was incubated 2 min at 37 • C. The reaction was initiated with the addition of xanthine oxidase solution (20 µL at 0.04 U/mL; in 50 mM phosphate buffer (pH 8) supplemented with 500 µM EDTA). The absorbance was first measured immediately after enzyme addition (blank) at 570 nm (Multiskan EX microplate reader (MTX Labsystems; Bradenton, Florida, USA)). After a 20 min incubation at 37 • C, 20 µL of HCl (0.6 M) were added to stop the reaction, followed by a second absorbance measurement at 570 nm. Rosmarinic acid was used as a positive control (95.71 ± 8.55% inhibition at 120 µg/mL).
For β-carotene bleaching assay, the emulsion was prepared by adding 500 mg of between 20 to 250 µL of β-carotene solution (2 mg/mL solution; in chloroform) and was followed by mixing in a round-bottom evaporation flask [25]. After evaporating the solvent in a rotary evaporator (35 • C), 25 mg of linoleic acid and 50 mL of distilled water were added, in this order. To produce the emulsion, the mixture was then gently homogenized using the rotary evaporator (rotary motion with no vacuum) at room temperature. The assay was carried out in a 96-well microplate, in which 50 µL of the extracts were added to 250 µL of the emulsion, followed by blank measurement at 450 nm. After 2 h incubation at 50 • C, the microplates were placed over ice to stop the reaction, in the dark, for two minutes, followed by a second absorbance measurement. Hydroethanolic extracts were dissolved in 10% (v/v) DMSO solution, and then they were tested to assure no interference with the assay. Rosmarinic acid was used as a positive control (IC 50 = 22.05 ± 1.02 µg/mL).
Hydroxyl ( • OH) and nitric oxide (NO • ) radicals scavenging assays were only performed for aqueous extract, due to ethanol interference, as the HE extracts are not fully water-soluble. Both assays were performed as described by Taghouti et al. (2020) [23]. With the exception of ABTS •+ scavenging assay (tested at 1 mg/mL), a range of concentrations of the extracts (0.1-1 mg/mL) was analyzed, and results are expressed as inhibition percentage and IC 50 , calculated according to equation 1. Distilled water was used as the negative control, and rosmarinic acid was used as a positive control for • OH scavenging assay without EDTA (43.27 ± 3.50% inhibition at 45 µg/mL) and for the NO • scavenging assay (44.43 ± 2.62% inhibition at 15 µg/mL).

Enzymatic Inhibition Assays
T. carnosus aqueous decoction and hydroethanolic extracts were studied for their potential neuroprotective, anti-aging, and anti-diabetic activities. These bioactivities were evaluated based on the capacity to inhibit key enzymes of target metabolic pathways. Acetylcholinesterase (AChE) and tyrosinase inhibition were evaluated for neuroprotection, tyrosinase and elastase for anti-aging activity, and α-amylase and α-glucosidase for antidiabetic activity.

Statistical Analysis
The experimental assays were performed for all the extracts obtained in each extraction method, with three experimental repetitions for each extract. Analysis of variance (ANOVA), followed by Tukey's multiple tests, were performed to analyze statistically significant differences. Correlations were evaluated using Pearson's coefficient (significant if p < 0.05). The IC 50 values were obtained from the dose-response assays described above and calculated, as described by Silva, et al. 2019 [26]. Principal component analysis (PCA) was used to evaluate inter-year variance of the individual phenolic components and performed as described by Ferreira, et al. 2020 [27]. The correlation between individual phytochemicals and antioxidant or enzymatic inhibition activities was performed through orthogonal partial least squares-discriminant analysis (OPLS-DA), as described by Martins, et al. 2022 [28]. Statistical analyses and graphic design were performed using Statistica (Version 14; TIBCO Software Inc., California, USA), SIMCA software (Version 14.1. Umetrics, Umea, Sweden), GraphPad Prism (Version 8; GraphPad Software Inc, California, USA), and Microsoft Office Excel (Microsoft Corporation, Washington, DC, USA).

Results and Discussion
Medicinal and aromatic plants' phytochemical composition is known to present high heterogeneity, even in species belonging to same genus. An example is the Thymus genus, since several species, such as T. pulegioides [24], Thymus zygis Loefl. ex L. [29], Thymus fragrantissimus [30], Thymus mastichina L. [23], Thymus × citriodorus (Pers.) Schreb. [31], and T. vulgaris [31], were harvested in the same location, grown in the same conditions, and whose phytochemicals were extracted using the same methodologies, and the extracts presented different yields, phytochemical profiles, and bioactivities. In addition to interspecies genetic variations, factors, such the vegetative phase and edaphoclimatic conditions, play a critical role in phytochemical composition variation [27,[32][33][34]. The latter is being widely discussed in light of climate changes induced by global warming. The effect of edaphoclimatic factors, such as temperature, precipitation or soil chemistry, and moisture, is well established as a determinant of secondary metabolite production [27,[32][33][34].
Within the Thymus genus, several studies have been performed to evaluate these variations, T. vulgaris and its essential oils being the most frequently addressed, given its significant economic impact. Lemos et al. (2017) studied the seasonal variance of T. vulgaris' essential oil from plants harvested in Brazil between July 2012 and July 2013, and it was observed that the October harvest presented higher antioxidant and antimicrobial activity, as well as an increase of 1.36 times in thymol and 1.85 times in p-cymene content, the major phytochemicals [16]. Additionally, using essential oils in T. vulgaris' and T. hyemalis, Jordán et al. (2006) have addressed the effect of the vegetative cycle on the phytochemical profile [35], while Pirbalouti et al. (2013) evaluated wild and cultivated samples of T. daenensis essential oil to ascertain the adaptability to crops, where this species produced higher contents of carvacrol or thymol under wild or cultivated growth, respectively [36]. As stated above, the effect of edaphoclimatic parameters on thyme extracts' phytochemical composition is poorly described. In the present research, we provided new data on the composition variation of aqueous and hydroethanolic extracts of Thymus carnosus over a three-year period, comparing both wild plants (harvested at location 1: L1) and plants cultivated in a botanical garden (harvested at location 2: L2), being the geographical locations shown in Figure 1A. L1 corresponds to plants grown in natural conditions, in sand dunes near to the coastline (229 m), and with low elevation (19 m), while L2 corresponds to plants originating in Arrábida National Park (L1) that were cultivated at the botanical garden of the University of Trás-os-Montes and Alto Douro, at an altitude of 451 m and 79 km from the coastline, which adapted to the northern inland climate and soil over a 12-year period. Regarding climate parameters, L2 presents the highest temperature variation, registering lower average temperatures in the winter months and higher average temperatures in summer months when compared to L1, but it overall presents a lower annual average temperature for the 2018-2020 period, as seen in Figure 1B,C. When considering the accumulated precipitation, L2 registered a significantly higher value than L1 ( Figure 1B,C).

Extraction Yield, Total Phenolic, Total Flavonoid, and Ortho-Diphenols Content
The extraction yields and results concerning total phenolic (TPC), total flavonoid (TFC), and ortho-diphenols (ODC) content, assessed for all extracts of T. carnosus harvested in L1 and L2, in the 2018-2020 period, are presented in Table 1. The extraction yield values ranged from 17.82% to 25.43%, being both the lowest and highest yields obtained for hydroethanolic extracts. With the exception of 2020's harvest at L1, all HE extracts present higher yield compared to the respective AD extracts, as reported for other thyme extracts using the same extraction methods [23,31]. The geographical location also affected the variation of extraction yields, with extracts from L2 presenting higher mean yield values than L1, this difference being more notorious in 2018's harvest. Concerning inter-year variance, the major differences were observed in 2020, where both HE extracts and AD extracts from L2 presented a decrease in the extraction yield.
T. carnosus AD extracts present higher TPC than the respective HE extract (Table 1), unlike the described for other Thymus extracts obtained using the same extraction methodologies (e.g., T. pulegioides [24] or T. mastichina [23]), in which TPC, assessed by the Folin-Ciocalteau method, was higher in HE extracts. The exception is L2 harvest of 2018, in which no significant differences were found between AD and HE extracts, identical to that reported for T. fragrantissimus [30]. However, Folin-Ciocalteau method limitations for TPC quantification are well described, since the reagent can be reduced by other chemical components in the extracts, such as, for example, proteins, thiols, carbohydrates, and amino-acids, where it is likely that AD extracts contain other reducing compounds, which might contribute to TPC overestimation [37]. This observation is supported by ODC and TFC quantification (Table 1), where most HE extracts have significantly higher ODC and TFC contents compared to the respective AD extract, L1 2019 s harvest being an exception. In fact, AD-L1-2019 presented higher TPC, ODC, and TFC than the remaining AD extracts, the ODP and TFC contents being in line with the respective HE extracts, and this extract was highlighted, even having the second lowest extraction yield within AD extracts (Table 1). Regarding the harvest location, overall L1 s extracts present the highest TPC, ODC, and TFC when compared to L2. Inter-year effect proved to induce variations in the extracts' phytochemical composition. As an example, ODC and TFC inter-year variation can be considered for L1-AD extract, as it is seen that from the first (2018) to the second (2019) year, where the content increased (in mg/g dry plant) 1.40 and 1.34 times, and then it reduced 0.94 and 0.86 times from 2019 to 2020, for ODC and TFC, respectively. In an opposing trend, L2 s HE extract presented the highest ODC and TFC in 2018, being followed by a reduction in the following years, which then present similar values. Table 1. Extraction yields and total phenolics, ortho-diphenols, and flavonoid content of Thymus carnosus extracts.

Year
Loc. E. M.

Extraction Yield
Total Phenols (mg Caffeic Acid eq./g) Ortho-Diphenols (mg Caffeic Acid eq./g) Total Flavonoids (mg Catechin eq./g) Thus, considering these results, geographical location induces a clearer pattern in the presented data, while climate variations present less predictable results. Regarding HE extracts, a previous report presented TPC values of 41.89 and 45.47 mg/g dry plant for flowering (July) and post-flowering (October) phases, respectively, both higher than the values here presented for T. carnosus HE extracts (Table 1).

Profiling and Quantification of Individual Compounds by HPLC-DAD and HPLC-ESI-MS n
To better understand the variations in the phytochemical composition, HPLC-DAD-MS n analysis was performed to analyze the variation of the main components of extracts, and that will be relevant to correlate with bioactivities. In Figure 2, HPLC-DAD quantification of total phenolic acids, total flavonoids, total phenols, and total terpenoids, for HE ( Figure 2B) and AD ( Figure 2C) extracts, was obtained from chromatograms, such as those presented in Figure 2A. The identification of individual compounds was performed with HPLC-ESI-MS n and by comparison of the literature. When considering the sum of all identified and quantified compounds by HPLC-DAD, L1-2020-HE arises as the harvest with the higher content in both phenolics and terpenoids ( Figure 2B), while, in AD extracts ( Figure 2C), L1-2019-AD is the harvest, presenting higher content in phenolic compounds, as it was also observed in Table 1 for TPC, ODP, and TFC contents in AD extracts. Martins-Gomes, et al. 2018 [9] reported the total phenolics contents assessed by HPLC-DAD, as well as for extracts of T. carnosus harvested in July and October 2015, at flowering and postflowering stages, respectively, with values ranging from 42.92-60.18 mg/g dry AD extract and 146.09-166.45 mg/g dry HE extract. In Table 2, values range from 37.96-99.86 mg/g dry AD extracts, thus most extracts are within the range of those previously reported. However, L1-2019-AD harvest presents a phenolic content above the average, as observed through HPLC-DAD quantification ( Figure 2C and Supplementary Table S2). Regarding HE extracts, in the present research, we report phenolic contents within 87.17-122.36 mg/g dry HE extract (Figure 2B and Supplementary Table S1), and all values are lower than the ones previously reported by Martins-Gomes, et al. 2018 [9] for the 2015 harvest. Given the extractability limit of aqueous extraction method, we observe that the contents here presented for AD extracts are in-line with the previous reported [9], while the increased extractability of exhaustive hydroethanolic extraction allows the extraction of all phenolic compounds, revealing that, when compared to July and October 2015 harvests, a harvest in November presents a decrease in total extractable phenolic compounds, most likely arising from changes in the plants' vegetative phases.   Abbreviations: AD: aqueous decoction and HE: hydroethanolic extractions; Loc.: location; E.M.: extraction method; Tukey's post hoc test was performed to analyze significant statistical differences between extraction methods (*), between harvest years for the same location (different capital letters), and between locations for the same year (different lowercase letters) if p < 0.05. Results are presented as mean ± standard deviation (n = 3). Trolox was used as the positive control for ABTS •+ scavenging, and RA was used for all the other assays, please see Section 2.4.
In Figure 3 (phenolic acids), Figure 4 (flavonoids), Figure 5 (terpenoids), and in Supplementary Tables S1 and S2, we present the quantification of individual compounds identified in T. carnosus extracts from L1 and L2 harvests. Overall, T. carnosus extracts present a similar profile to extracts reported previously [9] regarding phenolic acids with the presence of RA and salvianolic acids as major components, the identification of two salvianolic acid B/E isomers being identified, and T. mastichina [23] and T. zygis [29] extracts were also identified, but they were not previously described in T. carnosus extracts.
Regarding the effect of the location, edaphoclimatic factors may play a significant role. As described above, L1 is the natural habitat of T. carnosus. In fact, extracts obtained from harvests in L1 generally present higher content in phytochemicals. In L2, the average temperature in October and November is lower when compared to L1, which could induce an earlier end of the vegetative phase, thus reducing phytochemicals' production. A second hypothesis is the effect of drought stress, which is linked to increased secondary metabolite production [38]. Within Portugal's various edaphoclimatic zones, L1 is within an upper thermomediterranean and dry sub-humid zone, with higher aridity and less precipitation, while L2 is within lower supramediterranean and humid/upper mesomediterranean and humid zones [39]. Thus, the higher availability of water in plants adapted to L2 s climate may justify, in part, the decrease in phytochemicals' production.  [9] for plants harvested in July and October 2015 (0.92-5.76 mg/g AD extracts and 8.57-9.02 mg/g HE extracts). Contributing to this higher content in flavonoids, the presence of two apigenin derivatives, two eriodictyol derivatives, two quercetin derivatives, and six luteolin derivatives, from which quercetin-O-hexoside and luteolin-O-hexoside are highlighted as major compounds within the extracts (Figure 4), can be highlighted. Concerning HE extracts, quercetin-O-hexoside presented a larger variation in its content in L1. While 2018 harvest had a similar content to L2 harvests in the three years, for the 2020 harvest at L1, we report the highest value for T. carnosus extracts with 25.92 mg/g extract. In accordance with the compounds described above, AD extracts obtained from 2019 harvest at L1 presented not only the highest concentration of quercetin-O-hexoside, but also of luteolin-O-hexoside, this being last the major phenolic compound of the 2020 harvest at both L1 and L2. As expected, the highest concentrations were obtained for HE extracts, L1-2020-HE once again being highlighted as the extract with higher content in luteolin-O-hexoside (19.63 mg/g), and overall, luteolin-O-hexoside content in HE extracts ranged between 13.04 and 19.63 mg/g. The concentration of this compound in other Thymus species HE extracts was reported for T. carnosus [9] (4.61 mg/g extract), T. pulegioides [24] (6.31 mg/g extract), and T. zygis [29] (3.64 mg/g extract), which were lower than the   Results presented in Figures 3 and 4 and Supplementary Tables S1 and S2 were used to perform a statistical analysis of inter-year variance based on the quantification of its phenolic individual components through principal component analysis (PCA), which is presented in Figure 5. In Figure 5A,B we present the scatter plots for each extract (and its replicates) as a function of two PCs, explaining 42.45% (PC1) and 19.23% (PC2) of the variation for HE extracts, as well as 64.31% (PC1) and 20.60% (PC2) for AD extracts. PC2 correlates positively with quercetin derivatives in both HE and AD extracts. Figure 5A represents the variance of HE extracts, in which it is possible to observe a major cluster, In a previous report, T. carnosus extracts' phytochemical composition, during flowering and post-flowering stages (July and October 2015, respectively), presented salvianolic acids' A isomer (SAA isomer) and K (SAK), determined as the major phenolic compounds of both AD and HE extract, ranging from 14.87-27.50 mg/g and 61.92-67.34 mg/g of SAA isomer, as well as 12.53-19.66 mg/g and 38.51-65.33 mg/g of SAK, in AD and HE extracts, respectively [9]. The values for inter-year variance of phenolic acids (Figure 3), for harvests in November, report SAA isomer contents ranging from 6.57-12.24 mg/g in AD extract and 12.83-17.90 mg/g in HE extracts, both lower than the ones from October 2015 harvest [9]. of its higher RA, salvianolic acids, caffeic acid, and quercetin derivatives content. L1-2020-AD correlates poorly with PC2, being that these extracts are differentiated by a higher content in luteolin derivatives. L1-2018-AD, presents higher content of apigenin-(6,8)-Cdiglucoside, the second highest caffeic acid concentration, but lower contents of quercetin derivatives, which assume a distinct position from the other L1 harvests. Of all the analyzed extracts, the AD extracts from L2 present the least variation among them. Concerning the pentacyclic diterpenoids, (Figure 6 and Supplementary Table S1), we observed that L1 overall presented a higher content of both OA and UA when compared to L2 (with exception of UA content in 2019), and within the 2018-2020 period, 2020′s harvest presented the highest content of terpenoids (OA: 55.80 mg/g and UA: 51.84 mg/g). Oleanolic (OA) and ursolic (UA) acid have been described in T. carnosus HE extracts for Regarding RA content, T. carnosus AD extracts content in RA ranges between 4.48 and 20.12 mg/g extract, while HE extracts range between 19.99 and 29.94 mg/g, being the major phenolic acid in these extracts ( Figure 3). Considering T. carnosus extracts from flowering and post-flowering stage, RA was less predominant when compared to salvianolic acids, with a content of 0.16-4.40 mg/g AD extract and 27.84-29.07 mg/g HE extract [9]. We hypothesize that the much higher amount of salvianolic acids in the plants harvested in July and October 2015 limited RA extraction in AD extracts, whilst, in the present research, a lower content in other phenolic acids allowed a higher RA content in AD extracts, since the overall maximum extractable content (evaluated through HE extracts) is lower or similar to the previous report [9]. Regarding inter-year variance, RA content in HE extracts is overall stable in the three-year period in both locations. On the other hand, AD extracts present an increasing pattern for L2, where RA content increases 1.8 times between 2018 and 2020, whilst, in L1, there were no significant differences between these two years, but L1-2019 presented an increase of 1.83 times in RA content, being much richer in this phenolic acid, as was also observed for salvianolic acids. Raudone, et al. 2017 [21] studied the variation of rosmarinic acid through the various vegetative stages of 8 thyme species (70% hydroethanolic extracts). In all species, the content in RA was decreased in the end of vegetative phase, being, for example, Thymus praecox ssp. arcticus RA's content 6.95 times higher in May/June, when compared to August/September [21]. Thus, T. carnosus RA content at L2 (19.99-23.24 mg/g extract) increased in plants harvested in November, although it decreased when compared to July (29.07 mg/g extract [9]; between 1.25 and 1.45 times higher), and it was is less expressive when compared to the species described by Raudone, et al. 2017 [21], harvested in September, revealing the stability of RA content in T. carnosus extracts even at later harvests.
Nevertheless, the major difference between the data here presented and the previous report is concerning the flavonoids content, presented in Figure 4 and Supplementary Tables S1 and S2. The extracts from L1 and L2 (2018-2020 harvests) exhibited a flavonoid content ranging from 17.24 to 45.69 mg/g for AD extracts and 34.94 to 60.63 mg/g for HE extracts, higher values than the ones reported by Martins-Gomes, et al. 2018 [9] for plants harvested in July and October 2015 (0.92-5.76 mg/g AD extracts and 8.57-9.02 mg/g HE extracts). Contributing to this higher content in flavonoids, the presence of two apigenin derivatives, two eriodictyol derivatives, two quercetin derivatives, and six luteolin derivatives, from which quercetin-O-hexoside and luteolin-O-hexoside are highlighted as major compounds within the extracts (Figure 4), can be highlighted. Concerning HE extracts, quercetin-O-hexoside presented a larger variation in its content in L1. While 2018 harvest had a similar content to L2 harvests in the three years, for the 2020 harvest at L1, we report the highest value for T. carnosus extracts with 25.92 mg/g extract.
In accordance with the compounds described above, AD extracts obtained from 2019 harvest at L1 presented not only the highest concentration of quercetin-O-hexoside, but also of luteolin-O-hexoside, this being last the major phenolic compound of the 2020 harvest at both L1 and L2. As expected, the highest concentrations were obtained for HE extracts, L1-2020-HE once again being highlighted as the extract with higher content in luteolin-O-hexoside (19.63 mg/g), and overall, luteolin-O-hexoside content in HE extracts ranged between 13.04 and 19.63 mg/g. The concentration of this compound in other Thymus species HE extracts was reported for T. carnosus [9] (4.61 mg/g extract), T. pulegioides [24] (6.31 mg/g extract), and T. zygis [29] (3.64 mg/g extract), which were lower than the concentration presented in this research, while being in line with that reported for T. mastichina [23] (20.85 mg/g extract), which also presented similar values for quercetin-Ohexoside (20.34 mg/g extract) to those here reported for T. carnosus HE extracts.
Results presented in Figures 3 and 4 and Supplementary Tables S1 and S2 were used to perform a statistical analysis of inter-year variance based on the quantification of its phenolic individual components through principal component analysis (PCA), which is presented in Figure 5. In Figure 5A,B we present the scatter plots for each extract (and its replicates) as a function of two PCs, explaining 42.45% (PC1) and 19.23% (PC2) of the variation for HE extracts, as well as 64.31% (PC1) and 20.60% (PC2) for AD extracts. PC2 correlates positively with quercetin derivatives in both HE and AD extracts. Figure 5A represents the variance of HE extracts, in which it is possible to observe a major cluster, with clear separation of L1-2019 and L1-2020-HE. This arises from the fact that both extracts present higher content in quercetin-O-hexoside, but L1-2019-HE does not present the same pattern for the RA and SAA isomers, which have less correlation with PC1. The remaining HE extracts are mainly grouped within a major cluster, where L2-2018 presents the major correlation to PC1 and lowest correlation to PC2, thus representing a strong correlation with apigenin-(6,8)-C-diglucoside and with one eriodictyol derivative and caffeic acid.
Regarding AD extracts ( Figure 5B,D), a larger dispersion is observed between harvests, with a cluster being formed by L2 harvests in the 2018-2020 period, and then a strong dispersion is observed within L1 in the three-year harvest. L1-2019-AD clearly stands out as varying according to PC2, and it is correlated negatively with PC1 as a result of its higher RA, salvianolic acids, caffeic acid, and quercetin derivatives content. L1-2020-AD correlates poorly with PC2, being that these extracts are differentiated by a higher content in luteolin derivatives. L1-2018-AD, presents higher content of apigenin-(6,8)-C-diglucoside, the second highest caffeic acid concentration, but lower contents of quercetin derivatives, which assume a distinct position from the other L1 harvests. Of all the analyzed extracts, the AD extracts from L2 present the least variation among them.
Concerning the pentacyclic diterpenoids, (Figure 6 and Supplementary Table S1), we observed that L1 overall presented a higher content of both OA and UA when compared to L2 (with exception of UA content in 2019), and within the 2018-2020 period, 2020 s harvest presented the highest content of terpenoids (OA: 55.80 mg/g and UA: 51.84 mg/g). Oleanolic (OA) and ursolic (UA) acid have been described in T. carnosus HE extracts for the first time by Martins-Gomes, et al. 2018 [9], who reported a high content of both pentacyclic triterpenoids for both flowering stage (OA: 39.43 mg/g and UA: 75.17 mg/g) and post-flowering stage (OA: 41.74 mg/g and UA: 75.76 mg/g). A second conclusion can be retrieved when comparing vegetative phases, where, while OA content is similar to those presented by Martins-Gomes, et al. 2018 [9], the harvests of July and October favour a higher content in UA, whose highest content was 1.45 times higher than the observed for L1-2020-HE. Therefore, changes in the vegetative phase may induce differences in the production of such phytochemicals, that, similar to the ones observed for the phenolic acids described above, have a reduced concentration in November. Nevertheless, this comparison arises from an initial study from Martins-Gomes, et al. 2018 [9] with a single year analyzed. Additional studies for inter-year variation are essential to ascertain T. carnosus extracts' phytochemical variation in response to vegetative phase changes and climate adaptation. It is possible to observe that November also benefits the production of flavonoids derivatives, and L1 (the location from where the species is native) presents overall higher content in the phytochemicals here analyzed, most likely a result of a better adaptation to the edaphoclimatic conditions.  [9], the harvests of July and October favour a higher content in UA, whose highest content was 1.45 times higher than the observed for L1-2020-HE. Therefore, changes in the vegetative phase may induce differences in the production of such phytochemicals, that, similar to the ones observed for the phenolic acids described above, have a reduced concentration in November. Nevertheless, this comparison arises from an initial study from Martins-Gomes, et al. 2018 [9] with a single year analyzed. Additional studies for inter-year variation are essential to ascertain T. carnosus extracts' phytochemical variation in response to vegetative phase changes and climate adaptation. It is possible to observe that November also benefits the production of flavonoids derivatives, and L1 (the location from where the species is native) presents overall higher content in the phytochemicals here analyzed, most likely a result of a better adaptation to the edaphoclimatic conditions. The variation of oleanolic and ursolic acid was previously evaluated in the Thymus genus, in a single year harvest of plants grown in the same location, to study the variation induced by the vegetative stage in various species. In the various Thymus species analyzed by Raudone, et al. 2017 [21], Thymus sibtorpii Benth., Thymus austriacus Bernh. ex Rchb., Thymus × oblongifolius Opiz, and Thymus × citriodorus presented the highest content of OA and UA during the flowering stage, while Thymus serpyllum L., Thymus pulegioides, and Thymus longicaulis presented the highest content at the end of vegetative phase. Thus, we observe that pentacyclic terpenoids variation is also dependent on inter-species variation motivated by genetic variations, in addition to the edaphoclimatic parameters.

In Vitro Antioxidant Activity Assessment
The variation of herbal extracts' antioxidant capacity induced by biotic and abiotic factors has been previously addressed for several species, as the antioxidant activity is one of the main sought-after bioactivities for these products. We report here, for the first time, the study of T. carnosus extracts antioxidant activity variation as a result of different edaphoclimatic parameters and inter-year variance. The results for ABTS •+ , • OH , NO • , O2 •− , and β-carotene bleaching assays are presented in Table 2.
ABTS •+ scavenging was used as a reference antioxidant assay to compare both the extracts here reported, but also to allow the comparison with other species. In Table 2, it The variation of oleanolic and ursolic acid was previously evaluated in the Thymus genus, in a single year harvest of plants grown in the same location, to study the variation induced by the vegetative stage in various species. In the various Thymus species analyzed by Raudone, et al. 2017 [21], Thymus sibtorpii Benth., Thymus austriacus Bernh. ex Rchb., Thymus × oblongifolius Opiz, and Thymus × citriodorus presented the highest content of OA and UA during the flowering stage, while Thymus serpyllum L., Thymus pulegioides, and Thymus longicaulis presented the highest content at the end of vegetative phase. Thus, we observe that pentacyclic terpenoids variation is also dependent on inter-species variation motivated by genetic variations, in addition to the edaphoclimatic parameters.

In Vitro Antioxidant Activity Assessment
The variation of herbal extracts' antioxidant capacity induced by biotic and abiotic factors has been previously addressed for several species, as the antioxidant activity is one of the main sought-after bioactivities for these products. We report here, for the first time, the study of T. carnosus extracts antioxidant activity variation as a result of different edaphoclimatic parameters and inter-year variance. The results for ABTS •+ , • OH, NO • , O 2 •− , and β-carotene bleaching assays are presented in Table 2.
ABTS •+ scavenging was used as a reference antioxidant assay to compare both the extracts here reported, but also to allow the comparison with other species. In Table 2, it is possible to observe that ABTS •+ scavenging results ranged between 0.15 and 0.25 mmol Trolox equivalent/DP. The highest scavenging value was obtained for L1-AD-2019, which also presented the highest TPC content (Table 1), and thus it correlated with the ABTS •+ scavenging for a higher content in reducing compounds, including phenolics and polysaccharides, for example. The lowest value was obtained for L2-HE-2020. Martins-Gomes, et al. 2018 [9], for T. carnosus, harvested in July October 2015, reported inhibitions of 0.14 and 0.21 mmol Trolox eq./DP, for AD and HE extracts, respectively, thus being comparable to the results here obtained. When compared to other Thymus species, our results align with the values reported for HE extracts T. vulgaris (0.22 mmol Trolox eq./DP) [31] and T. mastichina (0.20 mmol Trolox eq./DP) [23], whilst T. zygis (0.25 mmol Trolox eq./DP) presented higher inhibition for HE extracts. Concerning AD extracts, the T. carnosus extracts here presented performed better than T. x citriodorus (0.11 mmol Trolox eq./DP) [31] and T. mastichina (0.081 mmol Trolox eq./DP) [23], species commonly used for human consumption in herbal teas and as condiments. Concerning the inter-year variation, in AD extracts, L1-2018 presented an ABTS •+ scavenging activity significantly lower than the respective harvests in 2019 and 2020, whilst no significant variations were observed for L2. For the HE extracts, L1 s harvests did not present significant variations between years, however, in L2 s harvests, the 2018-HE-L2 extract presented significantly higher ABTS •+ scavenging. When comparing the variation induced by the different locations, AD extracts scavenging was only statistically different in 2020, where L1-2020-AD presented higher antioxidant activity. In the HE extracts, significant differences were observed in 2018, as L2-2018-HE presented higher scavenging.
Aiming to study the correlation between the various individual phytochemicals and each bioactivity, score-plots similar to those presented in Figure 5 for PCA were obtained using the OPLS-DA (orthogonal partial least squares discriminant analysis) model. This analysis allows an easier interpretation and presents advantages to data sets comprising a higher number of observations than the number of variables, being proposed as a promising tool to highlight statistically significant components in comparable data sets [28,40]. The OPLS-DA model uses an orthogonal signal correction filter to identify the variations related to the prediction of a quantitative response from the variations not related to the prediction [40]. We present the correlations obtained by the OPLS-DA model when significant correlations were observed, and the models present an adequate degree of validation (models' validation is presented in the Supplementary material). Using data obtained in ABTS •+ scavenging assays, score-plots were obtained using the OPLS-DA model, which are presented in Figure 7. Most of phytochemicals identified and quantified in T. carnosus extracts positively correlate with the ABTS •+ scavenging ( Figure 7A). Quercetin-O-hexoside-hexuronide, luteolin-O-hexoside-pentoside, and caffeic acid present the lowest correlation with this antioxidant activity in HE extracts. CA content does not present a large variation influenced by either year or location. Quercetin-O-hexoside-hexuronide is present at the highest concentration in L1-2019-HE, L1-2020-HE, and L2-2020-HE, which are within the extracts that present lowest ABTS •+ scavenging, thus being the components that least explain the increase in scavenging as provided by OPLS-DA model. Eriodictyol-O-hexoside isomer 2 and luteolin-O-hexoside-hexoside isomer 1 are the components which better explain the variation observed. Regarding AD extracts ( Figure 7B), acetyl-luteolin-Ohexoside-pentoside, SAK, and SAA iso were the phytochemicals with highest contribution to ABTS •+ scavenging variance. L1-2019-AD and L1-2020-AD present the highest content of these three phytochemicals among the AD extracts studied, thus supporting the correlation.
In addition to ABTS •+ , the scavenging of biologically relevant radicals was also performed, namely, for radicals such as • OH, NO • , and O 2 •− , as well as β-carotene bleaching as a model for lipid peroxidation. The capacity of extracts in • OH scavenging in the absence of EDTA was dependent on their geographical origin, as L1 presented higher activity, obtaining inhibitions greater than 50%, and L1-AD-2019 produced the lowest IC 50 (0.8 mg/mL; Table 2), once again correlating to the higher content in phenolic compounds.
When analyzing inter-year variance, L1-2020-AD presented a significantly lower inhibition when compared to 2018 and 2019, whilst, in L2, 2019 s harvest presented a higher inhibition, although not significantly higher than L2-2018-AD and L2-2020-AD. In the presence of EDTA, all extracts performed in a similar manner, with approximately 30% of inhibition as the maximum value obtained in these conditions. In similar manner to that described for ABTS, the OPLS-DA model was applied to data obtained from • OH scavenging assay in the absence of EDTA and is presented in Figure 7C. All phytochemicals correlate positively with • OH scavenging, with phenolic acids (by this order: CA, RA, SAA iso and SAK) being the components with the highest influence in • OH scavenging, in line with the phytochemical composition obtained by HPLC-DAD, where overall extracts from L1 present higher content in these phenolic acids.
acid, ursolic acid, SAK, and acetyl-luteolin-O-hexoside-pentoside, positively influence the extracts' ability to inhibit the bleaching. In AD extracts, with a lesser content in these phytochemicals, luteolin-O-hexoside-pentoside and SAA isomer present the highest contributions to explain the variance observed.
In sum, both extraction methods, year of harvest and location, induce significant variations in the overall antioxidant capacity of T. carnosus extracts. In addition, this activity, although reported for the various tested radicals, reveals that the differences in each extract's composition (within each extraction method) induce variations in the scavenging ability, highlighting a capacity that is dependent of the radical selected. Considering its health-promoting effects, it is here demonstrated that T. carnosus bioactive components might be effective tolls to counter the oxidative stress activity in biological systems, which still requires further studies.   (Figures S2-S4).
The capacity of extracts in NO • scavenging followed a pattern identical to • OH scavenging, as presented in Table 2. L1-AD-2019 produced both the highest inhibition (73.31%) and the lowest IC 50 (0.57 mg/mL) values. Similarly, L1 extracts presented higher ability to scavenge this radical, most likely due to higher content of phenolic compounds. When compared to a previous report of NO • scavenging by T. carnosus extracts [9] (41.79%; 1 mg/mL), our extracts presented higher potential, all achieving IC 50 values bellow 1 mg/mL. Considering the different harvests in the same location, extracts from plants harvested in 2019 presented higher inhibition, at 1 mg/mL, as well as lower IC 50 values (Table 2). OPLS-DA analysis for NO • scavenging ( Figure 7D) also revealed a similar pattern to ABTS •+ scavenging, with all compounds correlating with the activity, supported by the higher content in phytochemicals observed in L1-AD-2019, producing the highest inhibition and lowest IC 50 . The compound with higher correlation is also identical to ABTS assay, acetyl-luteolin-O-hexoside-pentoside, whose highest content is observed in L1-2019-AD.
In this research, we report, for the first time, T. carnosus extracts' ability to scavenge O 2 •− , as well as the potential to reduce lipid peroxidation ( Table 2). All extracts displayed significant capacity to scavenge O 2 •− using the xanthine oxidase/NBT assay. Overall, AD extracts present higher inhibition mean values, but these are only significant in the L1-2019 harvest. The inhibition values ranged between 31.41% (L1-HE-2018) and 49.72% (L1-AD-2019) for the most efficient extract ( Table 2). Unlike for • OH or NO • radicals' scavenging, there is no pattern for the location effect concerning the extracts' scavenging of O 2 •− , since, when comparing locations, L1-2018 performed more poorly than L2-2018, and, when comparing inter-year effect, it presented less activity than L1-2019 and L1-2020. L2-2018 presented less O 2 •− scavenging than the respective harvest at L1 and the other L2 harvests. T. fragrantissimus extracts capacity to scavenge this radical indicated an inhibition of 48.81% for AD extracts, as well as 49.63% for HE extracts [30], values similar to those obtained in the 2019 harvest (Table 2), which are also in line with results obtained for T. zygis methanolic extracts (40.3%) [41] and for T. vulgaris aqueous extracts (~45%) [42].
Concerning lipid peroxidation, as seen in Table 2, all extracts were able to greatly inhibit β-carotene bleaching, with the most effective extract being L1-HE-2020, as well as being the extract with higher phenolic and terpenoid compounds quantified by HPLC-DAD, with the maximum inhibition of 96.25% at 1 mg/mL. Overall, HE extracts presented higher maximum inhibitions, with the exception of L2-2018, but the IC 50 revealed a different pattern. Only in L1-2019 and L1-2020 were significant differences observed concerning the IC 50 values, which were found between extraction methods, with AD extracts presenting lower IC 50 values. Thus, HE extracts, most likely due to a higher phenolic content, generally achieve higher inhibition, whilst AD extracts achieve about 50% of inhibition with equal or slightly lower concentrations ( Table 2). Concerning inter-year variation, in both AD and HE extracts, an increase in β-carotene bleaching inhibition is observed with the year of harvest, as 2018 < 2019 < 2020, whilst, in L2, the harvests of 2019 present higher inhibition. OPLS-DA model analysis ( Figure 7E) revealed that, in HE extracts, compounds with increased hydrophobicity among the components of the extracts, namely, oleanolic acid, ursolic acid, SAK, and acetyl-luteolin-O-hexoside-pentoside, positively influence the extracts' ability to inhibit the bleaching. In AD extracts, with a lesser content in these phytochemicals, luteolin-O-hexoside-pentoside and SAA isomer present the highest contributions to explain the variance observed.
In sum, both extraction methods, year of harvest and location, induce significant variations in the overall antioxidant capacity of T. carnosus extracts. In addition, this activity, although reported for the various tested radicals, reveals that the differences in each extract's composition (within each extraction method) induce variations in the scavenging ability, highlighting a capacity that is dependent of the radical selected. Considering its health-promoting effects, it is here demonstrated that T. carnosus bioactive components might be effective tolls to counter the oxidative stress activity in biological systems, which still requires further studies.
Various extracts of Gingko biloba obtained from material harvested in various locations in India over three different seasons showed that, in addition to a direct correlation between the phenolic content and antioxidant activity, the highest values of both parameters were reported in the autumn, where the antioxidant activity was expressed as ABTS and DPPH radicals scavenging [45]. In Thymus spp., variation of antioxidant activity induced by seasonal changes and developmental stages was also reported. Using essential oils from T. vulgaris harvested in Brazil over a twelve months period, Lemos, et al. 2017 [16] observed that ABTS and DPPH radicals scavenging was higher in Spring, correlating with a high content of thymol and carvacrol. Regarding the developmental stages, aqueous extracts from Thymus hirtus Boiss. and Reut., harvested in Tunisia, in a single year, in various vegetative phases and rich in catechin and epicatechin, were reported to present the highest antioxidant activity at flowering stage [46]. In T. longicaulis methanolic extracts, rich in RA, salvianolic acid K and luteolin derivatives, having a more similar phytochemical profile to the extracts of T. carnosus here reported, it was observed that, in a nine months period, the scavenging activity of ABTS and DPPH radical varied greatly, with IC 50 values ranging 8.91-72.01 µg/mL and 9.50-64.61 µg/mL for ABTS and DPPH radicals, respectively, which provide a new insight into a significant variation of the antioxidant activity dependent of the date of harvest [20]. Nevertheless, even extracts from commonly consumed species, such as T. vulgaris, present reduced information regarding the variation of bioactivities induced by changes in the harvests' location and date.

Enzymatic Inhibition Assays
In Table 3 are presented the results for T. carnosus extracts' inhibition of key enzymes often described as druggable targets aiming neuroprotection (AChE and tyrosinase), antiaging (tyrosinase and elastase), and diabetes management through lowering the intestinal absorption of glucose (α-amylase and α-glucosidase). With the exception of reports on T. carnosus essential oil activity in AChE inhibition [47], we present here, for the first time, the evaluation of T. carnosus AD and HE extracts' anti-enzymatic activity described above. Abbreviations: AD: aqueous decoction and HE: hydroethanolic extractions; AChE: acetylcholinesterase; Loc.: location; E.M.: extraction method; Tukey's post hoc test was performed to analyze significant statistical differences between extraction methods (*), between harvest years for the same location (different capital letters), and between locations for the same year (different lowercase letters) if p < 0.05. Results are presented as mean ± standard deviation (n = 3). Quercetin was used as a positive control for AChE and elastase assay, kojic acid for tyrosinase assay, and acarbose for α-amylase and α-glucosidase assays, please see Section 2.5.
Concerning potential neuroprotection effect, all extracts were effective in inhibiting AChE and tyrosinase activity, the most efficient extracts being L1-2019-AD for AChE and L2-HE-2020 for tyrosinase. All AD extracts presented significantly higher capacity to inhibit AChE when compared to the respective HE extract, despite the lower content in phenolic compounds (quantified by HPLC-DAD; Figures 2-4). This may be explained by the presence of other water-soluble non-phenolic components that are extracted by the AD method that present the capacity to inhibit this enzyme. T. carnosus essential oil has been reported to inhibit AChE activity, with an IC 50 of 0.72 mg/mL [47], and thus it is in-line with the extracts here present, whose best inhibition value at 1 mg/mL (L1-2019-AD) was 61.47% inhibition. Rich in RA, T. pulegioides extracts presented a higher efficiency in inhibiting AChE, being able to reduce the enzyme activity in almost 90% at 0.5 mg/mL [24]. T. vulgaris ethanolic extracts produced an inhibition above 75% at 1 mg/mL [48], which is also higher than the values here reported (Table 3). A similar comparison can be performed for T. serpyllum ethanolic and aqueous extracts, which inhibited 50% of AChE activity at 0.25 and 0.35 mg/mL, respectively [49]. On the other hand, T. fragrantissimus extracts produced lower inhibitions, being able to inhibit only 27.30% at 0.5 mg/mL [30].
The OPLS-DA model was used to evaluate the influence of individual phytochemicals present in T. carnosus extract in AChE inhibition variation, as is presented in Figure 8A for HE extracts and Figure 8B for AD extracts. For HE extracts, eriodictyol-O-hexoside isomer 2, luteolin-O-hexoside-hexoside isomer 1, acetyl-luteolin-O-hexoside-pentoside, and the two quercetin derivatives are the individual components which most contribute to explain the increase in AChE inhibition. Quercetin was used as the positive control for this assay, as it inhibits AChE activity significantly, and thus quercetin's glycoside derivatives identified in T. carnosus HE extracts likely also present inhibitory activity. In AD extracts, luteolin-O-hexoside-hexoside isomer 1 is also a major contributor to the variance observed in AChE inhibition, namely, in its increase.
Considering the location effect, L1, AD, and HE extracts, from the 2018 and 2019 harvests, produce greater AChE inhibitions. For the 2020 harvest, identical AChE inhibition was obtained with L1 and L2 HE extracts (p > 0.05). However, L2-AD-2020 has a greater capacity to inhibit AChE than L1-AD-2020. Regarding the inter-year variation, L1-HE extracts produced no variation in AChE inhibition, while L1-AD extracts harvested in 2020 produced a significant decrease in AChE inhibition. For L2 s extracts, AD extracts' anti-AChE inhibition increases in 2019 and 2020, when compared to 2018, and HE extracts also present significantly higher inhibition in 2020.
The inhibition of tyrosinase activity by T. carnosus extracts is dependent on their phytochemical content, since the HE extracts always present higher inhibition than the AD extracts. However, L1-2020-HE, the extract with higher phenolic acids, flavonoids, and terpenoids quantified by HPLC-DAD did not produce the highest inhibitory effect. Overall, the extracts present potential for neuroprotection, with efficacy varying based on location and year of harvest, but without a defined pattern (Table 3). T. fragrantissimus AD extracts were also able to inhibit tyrosinase (56.30% at 0.5 mg/mL) [30], while T. pulegioides presented a much higher inhibition capacity (93.72% at 0.5 mg/mL) [24]. Both species showed higher potential than T. carnosus, whose extracts were able to inhibit a maximum of 42.54% for L2-2020-HE (Table 3).
T. carnosus extracts presented a poor efficacy regarding elastase inhibition, thus not being the primary potential choice for anti-aging applications, unlike T. fragrantissimus HE extracts that inhibited 48.69% of elastase activity [30]. T. carnosus extracts here presented (Table 3) showed a pattern similar to T. pulegioides [24], where no elastase inhibition capacity was observed for HE extracts, but AD extracts were able to reduce elastase activity. The highest inhibition obtained in the present research was 7.31% using L1-2019-AD (Table 3), which is also the best extract for AChE inhibition. As determined by OPLS-DA model (Figure 8), the higher contribute to α-glucosidase inhibition variation produced by HE extracts is due to the presence of luteolin-O-hexoside-pentoside, eriodictyol-O-hexoside isomer 1, and quercetin-O-hexoside-hexuronide, which may be connect to the presence of sugar residues in these glycoside derivatives. Future analysis of the role of individual compounds in this bioactivity should address the effect of the presence of sugar residues in comparison to the aglycone to ascertain if these moieties have higher affinity to the enzyme's active site. In fact, L2-2018-HE, the HE extract with higher α-glucosidase inhibition, is also the extract with higher content of luteolin-O-hexoside-pentoside and of the first eriodictyol-O-hexoside isomer identified, which correlates with the OPLS-DA model. For AD extracts, the compounds whose variation better explain the inhibition pattern are acetyl-luteolin-O-hexoside-pentoside and the salvianolic acids, factors correlated with the increase in phytochemicals observed in L1-2019-AD, which produced a significantly higher α-glucosidase inhibition. As with the results of tyrosinase inhibition, the inhibition of α-amylase activity is dependent on the phenolic and terpenoids content, since the inhibition induced by HE extracts is significantly greater than the respective AD extracts. Nevertheless, once again L1-2020-HE was not the best performing extract, suggesting that the overall composition and ratio between all components may play a role, which are opposed to the simple higher content in each individual component. The extracts were not able to induce a reduction in α-amylase activity greater than 6%, which are results in line with those obtained by Taghouti et al. (2018) using T. pulegioides extracts [24]. On the other hand, T. carnosus' potential for anti-diabetic applications may be achievable through α-glucosidase inhibition, where a maximum inhibition of 27.37% was achieved for L1-2019-AD.
Regarding the location, AD extracts from L1 (in all harvests) produced, on average, higher inhibition of tyrosinase activity than AD extracts from L2. While, for HE extracts, significantly higher inhibition was produced by extracts from L2 harvests in 2018 and 2019, but not in 2020. Regarding the harvest year variation, it was observed that AD-2019-L1 extract produced higher tyrosinase inhibition than the other AD extracts, while the best HE extract was from L2 during the harvest of 2020 when compared to the other harvest years.
The extraction method influenced anti-glucosidase activity, as well as L2-HE extracts from 2018 and 2020 harvests, which show higher content in phytochemicals and also show higher inhibition capacity. In an opposed pattern, L1-2019-AD produced significantly higher α-glucosidade inhibition when compared to the respective HE extract ( Table 3).
The contribute of individual phytochemicals to α-glucosidase inhibition variation was also evaluated by OPLS-DA model, and it presented in Figure 8C for HE extracts and Figure 8D for AD extracts.
As determined by OPLS-DA model (Figure 8), the higher contribute to α-glucosidase inhibition variation produced by HE extracts is due to the presence of luteolin-O-hexosidepentoside, eriodictyol-O-hexoside isomer 1, and quercetin-O-hexoside-hexuronide, which may be connect to the presence of sugar residues in these glycoside derivatives. Future analysis of the role of individual compounds in this bioactivity should address the effect of the presence of sugar residues in comparison to the aglycone to ascertain if these moieties have higher affinity to the enzyme's active site. In fact, L2-2018-HE, the HE extract with higher α-glucosidase inhibition, is also the extract with higher content of luteolin-O-hexoside-pentoside and of the first eriodictyol-O-hexoside isomer identified, which correlates with the OPLS-DA model. For AD extracts, the compounds whose variation better explain the inhibition pattern are acetyl-luteolin-O-hexoside-pentoside and the salvianolic acids, factors correlated with the increase in phytochemicals observed in L1-2019-AD, which produced a significantly higher α-glucosidase inhibition.
At 1 mg/mL, T. carnosus extracts presented an α-glucosidase inhibition~2.16 times and 2.71 times higher than T. fragrantissimus [30] and T. pulegioides [24] extracts at 0.5 mg/mL (the only concentration tested). Aqueous extracts of T. vulgaris (concentration not mentioned) were reported to inhibit 4% and 20% of α-amylase and α-glucosidase, respectively [50], which are also in-line with the results were reported. Therefore, T. carnosus extracts present a moderate potential for the reduction of post-prandial blood sugar uptake by inhibiting α-glucosidase. Moreover, despite the variations induced by the unique phytochemical composition of each extract, influenced by harvest year, location, and extraction method, all extracts present neuroprotective and anti-diabetic activity, and thus T. carnosus can be a reliable functional food for these applications. Although there are few reports for salvianolic acids' ability to inhibit these key enzymes, salvianolic acid B presents potential to inhibit AChE, and salvianolic acid C has been shown to inhibit α-glucosidase [51,52]. RA and pentacyclic triterpenoids have been shown to inhibit AChE, α-glucosidase, and α-amylase activity, and thus the harvest at earlier vegetative stages may increase this potential effect [48,[53][54][55]. Thus, it is also relevant for future studies to further unveil the potential effects of T. carnosus extracts bioactivities, especially at the gastrointestinal tract, since AD extracts present a wide array of potential bioactivities, and this is reported in the present research.

Conclusions
Considering the high number of species used for its health-promoting bioactivities and the large number of pharmaceutical products obtained from plants, it is essential to increase the screening for new medicinal plants and phytochemicals with pharmaceutical potential. In addition, due to the increasing number of species threatened by climate and anthropogenic factors, there is a risk of losing the potential applications of a high number of species that were never considered for screening. Aiming to increase the interest in T. carnosus and encourage a sustainable crop, in the present research we complement the previous work related to T. carnosus extracts phytochemical composition and antiproliferative activity by further analyzing its neuroprotective, anti-aging, and antidiabetic activity, complementing the antioxidant activity and evaluating the effect of edaphoclimatic parameters in the phytochemical composition and bioactivities. In addition to a unique phytochemical composition, T. carnosus extracts are a promising source of rosmarinic and salvianolic acids, as well as glycoside derivatives of various flavonoids, where the vegetative stage plays a significant role in the concentration of each phenolics subclass.
Despite the variations in the phytochemical profile, the extracts maintained significant antioxidant and neuroprotective activities and moderate anti-diabetic activity in the three-year period, thus revealing the potential to induce health-promoting activities when included in the diet, as demonstrated by the AD extracts, but also as a source of phytochemicals. Further studies on T. carnosus extracts should be performed to complement the information available, such as studies on its antioxidant and anti-diabetic potential in in vitro and in vivo experimental models.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/antiox12030668/s1, Table S1: Phytochemical composition of T. carnosus hydroethanolic (HE) extracts obtained by HPLC/DAD-ESI/MS n analysis.; Table S2: Phytochemical composition of T. carnosus aqueous decoction (AD) extracts obtained by HPLC/DAD-ESI/MS n analysis. Figure S1. Chromatograms of hydroethanolic (A) and aqueous decoction (B) extracts of T. carnosus obtained by HPLC-DAD. This are examples of the chromatograms used for phytochemicals' quantification, obtained at 280 nm for polyphenolic compounds, and the chromatogram used for terpenoids quantification at 210 nm (C) in the hydroethanolic extract. For peak identification please refer to Supplementary Tables S1 and S2. Compounds not identified in these chromatograms were identified in HPLC-ESI-MS n analysis but not in the HPLC-DAD chromatogram used for quantification.; Figure S2. Performance statistics for overall model, predictive components and orthogonal components for the correlation between ABTS •+ scavenging and T. carnosus extracts' phytochemicals (A: HE extract; B: AD extract) and respective plots of model's validation. Sig.: significance. Figure S3: Performance statistics for overall model, predictive components and orthogonal components for the correlation between • OH (A) and NO • (B) scavenging and T. carnosus AD extracts' phytochemicals and respective plots of model's validation; Figure S4: Performance statistics for overall model, predictive components and orthogonal components for the correlation between β-carotene bleaching inhibition and T. carnosus extracts' phytochemicals (A: HE extract; B: AD extract) and respective plots of model's validation; Figure S5: Performance statistics for overall model, predictive components and orthogonal components for the correlation between acetylcholinesterase inhibition and T. carnosus extracts' phytochemicals (A: HE extract; B: AD extract) and respective plots of model's validation; Figure S6: Performance statistics for overall model, predictive components and orthogonal components for the correlation between α-glucosidase inhibition and T. carnosus extracts' phytochemicals (A: HE extract; B: AD extract) and respective plots of model's validation.