2.1.1. Phenological Profiling by HPLC-PDA
A total of 43 compounds, involving major groups of secondary metabolites—phenolic and triterpenic compounds—were found in extracts of lingonberry leaves. Being an evergreen species, lingonberry keeps its leaves thick and waxy even during the coldest months
[4][15]. Our results showed that most compounds were detected at a quantifiable level throughout the year thus proving continuous biosynthesis and storage of secondary metabolites in lingonberry leaves. Contents of all compounds varied along the year and different groups of phenolic and triterpenic compounds showed different patterns of seasonal variation (accumulation patterns of particular constituents in young and old lingonberry leaves, collected throughout one year, are presented in
Supplementary Tables S1–S7). The content of secondary metabolites varies due to modifications in gene expression or their encoded protein activity involved in secondary metabolism in the presence of different developmental stages, stresses, and adaptation to them
[18]. The production of different compounds is not evenly affected by these factors, depending on the main biosynthesis pathway. The current understanding of metabolic pathways of all compounds has not been comprehensively defined yet and is still the topic of genetic studies. However, it has now been proposed that phenolic compounds synthesis is a combinatorial biosynthesis, including shikimate and acetate pathways, whereas triterpenic compounds are basically formed by the mevalonate pathway with some alternative routes and many different enzymes involved
[19][20][21]. Since phenological and physiological factors influence the biosynthesis of plant secondary metabolites, seasonality and harvest time stand as critical factors determining the chemical composition
[22].
The variance of sum contents of different subgroups of phenolics in young and old lingonberry leaves, collected during one year period, is shown in
Figure 1. The most important contributor to the leaves extracts was arbutin, accounting for 42–60% of total identified phenolics at different collection dates. Strong correlation between arbutin level and total identified phenolics (
r = 0.987,
p < 0.05) resulted in similar fluctuation curves (
Figure 1A,B). The content of arbutin in old lingonberry leaves was quite stable with the lowest coefficient of variation (CV)—7%, compared to other identified phenolic compounds. Nevertheless, the highest (
p < 0.05) levels were detected at the end of the summer and autumn. Levels of arbutin in young current-year lingonberry leaves were similar to those found in old previous-year leaves, except for a notable decrease in amount during the May–August months. These basic findings are partly consistent with a previous study showing that arbutin level is highest in autumn, suddenly decreasing during the burst of buds and slowly increasing during flowering
[23].
Figure 1. Variation in sums of (
A) all identified phenolics, (
B) arbutin, (
C) arbutin derivatives, (
D) catechins, (
E) proanthocyanidins, (
F) flavonol aglycones, (
G) kaempferol glycosides, (
H) quercetin glycosides, and (
I) phenolic acids (μg/g DW) in young (light line) and old (dark line) lingonberry leaves during one year testing period. Abbreviations of months have been defined in
Appendix A.
Similar accumulation and content decreasing patterns were obtained for other important contributors—catechins (the sum of contents of (+)-catechin and (–)-epicatechin) and proanthocyanidins (procyanidins A1, A2, A4, B1, B2, B3, C1), which made up to 15% and 24% of total identified phenolics, respectively (
Figure 1D,E). These groups were characterized by high CV, with the greatest variation (50%) of catechins in current-year lingonberry leaves. The highest (
p < 0.05) levels of most proanthocyanidins and catechins were observed in 16 November–11 January samples, while the content of these compounds was about 1.5-fold and 3-fold lower in old and young lingonberry leaves taken on 1 May–24 August, respectively. These results go beyond the previous report about increasing biosynthesis of flavanol monomers and oligomers in lingonberry leaves from May to September, and which found the greatest amounts in autumn
[24].
In contrast to these accumulation patterns, arbutin derivatives (hydroquinone and 2-
O-caffeoylarbutin), minor contributors—flavonol aglycones (kaempferol, quercetin) and kaempferol glycosides (nicotiflorin, astragalin, afzelin)—were present in an intermittently increasing trend from 1 May till 24 August in samples of young leaves, whereas in previous-year lingonberry leaves they appeared in more constant mode with the highest (
p < 0.05) concentrations on 7 February (3266.0 μg/g dry weight (DW)), 15 July (99.1 μg/g DW), and 20 September (262.7 μg/g DW) in case of arbutin derivatives, flavonol aglycones, and kaempferol glycosides, accordingly (
Figure 1C,F,G). Similarly, Bujor et al.
[24] reported that arbutin derivative—2-
O-caffeoylarbutin in lingonberry leaves exhibited a seasonal increase, that started in May.
The metabolic pattern of quercetin glycosides (quercitrin, quercetin-HMG-rhamnoside, 6′′-
O-acetylisoquercitrin, rutin, hyperoside, isoquercitrin, reynoutrin, guaiaverin, and avicularin) and phenolic acids (chlorogenic, cryptochlorogenic, neochlorogenic, and
p-coumaric acids) was markedly distinguished by the lowest (
p < 0.05) amounts in the 1 May samples of current-year leaves and sharply increased up to 13,443.2 μg/g DW on 16 November in case of quercetin glycosides and up to 2851.3 μg/g DW on 24 August in case of phenolic acids (
Figure 1H,I). In old lingonberry leaves, levels of quercetin glycosides (CV = 10%) ranged between 7725.2 μg/g DW on 4 April and 11,823.3 μg/g DW on 29 November, and the content of phenolic acids (CV = 15%) was in a range of 1245.8–2135.4 μg/g DW (1 May–29 November). The same increasing tendency of flavonol glycosides was obtained previously when assessing contents of flavonol glycosides in leaves of another
Vaccinium member—bilberry— in May, July, and September
[25].
Overall, in many cases levels of phenolics in old lingonberry leaves fluctuated but showed no discernible trend, whereas phenolics in current-year leaves showed sharply increases or decreases, outlining a major annual effect during the growth season. The observed changes in phenolic compounds during spring–summer seasons, preeminent in the beginning of May could be attributed to the production of new current-year leaves and intensified biosynthesis of secondary metabolites. A similar annual effect was reached by Sommavilla et al.
[26], which analyzed seasonal variation of phenolics in leaves of
Celtis australis L. It seems likely that the initial level of phenolics may be even higher than in previous-year leaves or be significantly lower with a slow increasing tendency, highly depending on the group of phenolic compounds. Seasonal variance is expected to be more stable in previous-years leaves due to slower metabolism processes, decay of some compounds, and also because young leaves, which are on the upper position of the stem, are more likely to be affected by UV-screening, acclimatization to light conditions, and other biotic and abiotic elicitors
[27]. A recent study calculated about 1.5–8-fold lower CVs of most phenolics in old leaves, compared to young leaves. CV in young leaves varied between 14 and 111%, whereas in old leaves it ranged between 7 and 104% with the highest variation of hydroquinone in both cases, followed by
p-coumaric acid, kaempferol, nicotiflorin, procyanidin A4, rutin, and (+)-catechin. It can thus be suggested that the metabolism of these compounds is strongly affected by environmental and phenological factors. It was reported in the literature that with the growth of evergreen plants, the content and yield of secondary metabolites tend to be higher and higher
[18]. Our results showed that average annual amounts of particular compounds in old leaves were slightly higher (1.1–1.2-fold) than those in young leaves in cases of arbutin, catechins, and proanthocyanidins, and lower (1.2–1.5-fold) when considering amounts of arbutin derivatives, flavonols, and phenolic acids. The reason for these partly contradictory results is still not entirely clear, but it can be supported by variable decay processes of compounds and different adaptations to weather conditions.
Contents of neutral triterpenoids, triterpenoid acids, sterols, and total identified triterpenoids in lingonberry leaves were also variable according to collection time (
Figure 2 represents the variance of sum contents of different subgroups triterpenoids in young and old lingonberry leaves, collected during a one-year period). The accumulation pattern of total identified triterpenoids was mainly influenced by contents of neutral triterpenoids (the sum of contents of betulin, erythrodiol, uvaol, lupeol, α-amyrin, β-amyrin, and friedelin), which made up 40–69% of total identified triterpenoids at different collection dates (
r = 0.972,
p < 0.05) and can be considered as a prevailing group in lingonberry leaves with the principal component—α-amyrin (
Figure 2A,B). The highest (
p < 0.05) contents of most neutral triterpenoids in young leaves emerged in colder months (16 November–7 February) with the greatest total content of 1472.0 μg/g DW on 28 December, whereas the accumulation peak in old leaves was reached earlier—on 6 September (total content of 2714.7 μg/g DW). The lowest (
p < 0.05) levels of neutral triterpenoids from 29 May–13 June and 28 December–11 January were observed in young and old leaves, respectively. These findings correlate favorably with Szakiel et al.
[12] and further support the concept that the content of neutral triterpenoids is considerably higher in old leaves during warm months, while in winter, content in young leaves reaches or even surpasses neutral triterpenoids level in old leaves.
Figure 2. Variation in sums of (A) all identified triterpenoids, (B) neutral triterpenoids, (C) triterpenoid acids, and (D) β-sitosterol (μg/g DW) in young (light line) and old (dark line) lingonberry leaves during one year testing period.
Very similar fluctuation patterns of triterpenoid acids (maslinic, corosolic, betulinic, oleanolic, and ursolic acids) harvested throughout one year were found in young and old lingonberry leaves, but contents were more stable—approx. 1.5 and 1.3-fold lower values of CV, respectively, compared to neutral triterpenoids, were calculated (
Figure 2C). The maximum sum values of triterpenoid acids with predominant ursolic acid out of this group occurred on 13 December (927.0 μg/g DW) and the lowest one was registered on 29 May (255.1 μg/g DW) in the case of current-year lingonberry leaves. Meanwhile, concentrations in old lingonberry leaves ranged between 724.8 and 1285.5 μg/g DW in samples taken on 21 March and 19 October, respectively. Rahajanirina et al.
[28] indicated that the biosynthesis of triterpenoid acids, as well as other triterpenoids, may be enhanced during the cold and rainy months, while summer may be less favorable for biosynthetic processes.
Further data would be needed to determine exactly how the phenological stage affects the accumulation of all sterols, found in lingonberry leaves, but our preliminary study of variance of β-sitosterol, which has been reported to be the main compound out of this group
[29] suggest that biosynthesis of sterols intensifies in late autumn. The highest (
p < 0.05) concentration of sterols (β-sitosterol) was observed from 19 October–1 November (up to 390.7 μg/g DW) and 16–29 November (up to 431.2 μg/g DW) in samples of old and young lingonberry leaves, accordingly (
Figure 2D). Different periods of year were outlined because of low contents; about a 4-fold lower amount of β-sitosterol, compared to its highest value, was found in samples harvested on 11 January in case of old lingonberry leaves and 8-fold lower on 29 May in the case of current-year leaves.
Our findings on seasonal variance clearly propose that biosynthesis of most triterpenoids appears to be most intensive in current-year leaves during late autumn–winter months, before production of new vegetative buds, whereas old leaves display an increasing tendency until the winter. The fact that old leaves are located in the bottom of the stem—usually beneath the snowpack during winter should be taken into consideration
[27]. Some researchers claim that triterpenoid accumulation in leaves during snow covering acquires an important protective role against herbivores and fungal pathogens
[30]. It is also noteworthy that the average annual amounts of all tested groups of triterpenoids in previous-year leaves were higher (1.2–1.8-fold) than those found in current-year leaves, reaching a similar level only in December–January. In contrast to phenolics, biosynthesis of triterpenoids seems to have an increasing trend with the increase of growth years, with likely slower decay processes and sustained metabolic activity. In addition, a greater variance of triterpenoids was found. Mean values of CV were twice and nearly 1.5-fold higher in old and young lingonberry leaves, accordingly. The highest variance due to phenological stages of the following individual triterpenoids was established: lupeol (CV = 70%), erythrodiol (68%), and maslinic acid (65%), while the most constant values were observed for ursolic (31%) and oleanolic acids (36%), suggesting different accumulation and adaptation processes to the growing season. Considering similarities to phenolics, the same annual increasing pattern of triterpenoids could be seen from May.
2.1.3. Environmental Factors Analysis
Variation in phenolic and triterpenoids content may be affected not only by phenological stage and age of leaves but additionally by seasonal factors—temperature, sunshine duration, precipitation, humidity, and plant defense mechanisms to climatic conditions
[32][33]. In response to weather fluctuations or negative abiotic factors, plants are able to adapt to the changes and trigger variation in secondary metabolites by their displacement between tissues and intensified accumulation
[34]. Since environmental factors are crucial determinants for the changes in plant secondary metabolites, the contents of total identified (
Figure 4) and individual compounds were examined for any correlations to known climatic conditions.
Figure 4. Correlations between sunshine duration (h), temperature (°C), precipitation (mm), and humidity (%) throughout the year with total identified phenolic and triterpenoid contents (µg/g DW) in lingonberry leaves.
Among all abiotic factors, light plays an irreplaceable role in regulating plant phenology due to its direct relation to photosynthetic activity
[35]. Light provides the energy, which is required for photosynthesis, promoting plant growth, and inducing or regulating plant metabolism through photoreceptors
[32]. Most commonly moderate light intensity has a stimulatory effect on the formation of secondary metabolites, however, the effect can be different, depending on the plant, organ, and group of compounds
[18]. The data analyses showed that levels of total phenolics and triterpenoids were negatively correlated with the temperature throughout the year (
r = −0.898 and −0.671, respectively,
p < 0.05). The strongest negative correlation to sunshine duration was found with proanthocyanidins and catechins from the phenolics group and with triterpenic acids from the triterpenoids group (
r = −0.926 and −0.883, −0.790 respectively,
p < 0.05). It is still controversial why plants growing under shading conditions may experience an increase of secondary metabolites. Excessive levels of solar radiation may reduce carbon metabolism and cause photoinhibition by impairing the photosynthetic reaction centers of the chloroplasts
[35]. Among all identified compounds, only kaempferol glycosides and flavonol aglycones correlated positively with sunshine duration (
r = 0.683 and 0.611, respectively,
p < 0.05), suggesting greater demand for light energy for their biosynthesis, compared to other constituents. In line with our study, a shorter sunshine duration caused a double reduction in the content of flavonol aglycones in aerial parts of
Xantium spp. earlier
[36]. Correlation results for sunshine duration indicated that it may be possible to adjust the desirable accumulation of secondary metabolites by adjusting light and dark regimes.
Plants’ growth and development are also directly linked with air temperature, which regulates the responding genes. A temperature range of 17–25 °C is generally considered optimal for the maintenance of plant cells and productivity; however, each plant species may exhibit optimum growth and metabolism under different temperatures
[18][32]. Some plants, like lingonberry, can survive low-temperature stress, tolerate –40 °C or even lower temperatures, by inducing different physiological, biochemical, and molecular changes and thus increasing levels of particular compounds
[4][14][21]. Cold acclimation phenomenon in our study can be most clearly evidenced by the negative correlation of total identified phenolics with air temperature (
r = −0.736,
p < 0.05) and the strongest correlations with catechins, proanthocyanidins, and arbutin out of the phenolic group (
r > −0.700,
p < 0.05). In line with our results, an adverse secondary metabolism response to higher temperature was reported to be expected in grapes, which are widely known as a rich source of proanthocyanidins and catechins
[37]. Concerning the content of total triterpenoids, no significant correlation was found to temperature, only betulin as an individual triterpenoid was distinguished by a strong positive correlation (
r = 0.781,
p < 0.05).
Drought stress is known as a further key environmental factor profoundly affecting plant metabolism by regulating carbon allocation from the roots and nutrient circulation
[38]. Severe water deficit has been considered to reduce plant development and decrease biomass production, but as reported recently, plants may indeed accumulate higher contents of secondary metabolites as a response to drought
[14][39]. A moderate negative correlation was found between total phenolics content and precipitation level throughout the year (
r = −0.602,
p < 0.05) with the strongest correlation with catechins out of this group (
r = −0.766,
p < 0.05). It was reported that deficient water uptake induced high increases in contents of (+)-catechin and (–)-epicatechin in
Crataegus spp. leaves
[40], thus supporting the idea that drought may promote the production of tannins in plants. Some positive significant correlations between arbutin derivatives, kaempferol glycosides, phenolic acids, and flavonol aglycones to precipitation level (
r up to 0.866,
p < 0.05) were have also been found, suggesting that high water uptake is favorable for circulation and metabolism of these compound. Meanwhile, rainfall level was not closely linked to the contents of triterpenoids, and this statement is consistent with what has been found in a previous study by Alqahtani et al.
[22].
Even though air humidity has more impact on the composition of volatile compounds, it also can shift the metabolism of other secondary metabolites even in perennial plants with wax layer leaves, as in lingonberries
[38]. High humidity can reduce growth rate, the biomass of leaves, and bud size, and result in an earlier burst of buds, but simultaneously may lower hydrophobicity of leaf surface, exposing the leaves for fungal pathogen attacks
[41]. Hence, positive changes in contents of phenolics and triterpenoids under higher air humidity can be explained as a defense mechanism. This concept ties well with our study, wherein strong positive correlations were found between total phenolic or triterpenoid contents and air humidity (
r = 0.742, and 0.814, respectively,
p < 0.05) with strongest correlations to contents of proanthocyanidins, arbutin, and triterpenoid acids (
r > 0.700,
p < 0.05).