Essential Oils of Artemisia frigida Plants (Asteraceae): Conservatism and Lability of the Composition

Plants of arid regions have adapted to harsh environments during the long span of their evolution and have developed a set of features necessary for their survival in water-limited conditions. Artemisia frigida Willd. (Asteraceae) is a widely distributed species possessing significant cenotic value in steppe ecosystems due to its high frequency and abundance. This study examines different patterns of formation of essential oil composition in A. frigida plants under the influence of heterogeneous factors, including climate and its integral characteristics (HTC, Cextr, SPEI and others). The work is based on the results of our research conducted in Russia (Republic of Buryatia, Irkutsk region), Mongolia, and China, from 1998 to 2021. A total of 32 constant compounds have been identified in the essential oil of A. frigida throughout its habitat range in Eurasia, from Kazakhstan to Qinghai Province, China. Among them, camphor, 1,8-cineol and bornyl acetate are the dominant components, contained in 93–95% of the samples. Among the sesquiterpenoids, germacrene D is the dominant component in 67% of the samples. The largest variability within the composition of the essential oils of A. frigida is associated with significant differences in the climatic parameters when plants grow in high-altitude and extrazonal conditions.


Introduction
Currently, the issue of biota response to ongoing climate warming, along with the back-reaction between "biota" and "climate", is at the cutting edge of global biology and ecology.There is no doubt that climate change is happening and will continue to happen in the foreseeable future.Vegetation is forced to respond to biotically significant changes in climate, whether they are significant once, are episodic oscillations, or sustainable long-term trends.Therefore, a key scientific problem is the study of mechanisms of plant adaptation to the conditions of short-term or long-term exposure to adverse environmental factors.
During the current climate warming, arid territories are occupying increasingly extensive areas.Plants in arid regions have adapted to harsh environmental conditions during the long span of their evolution, and have developed a set of specific morphological, physiological and phenological features necessary for their survival in water-limited conditions.The study of these plants growing in natural populations will allow to demonstrate the stability of the plant locus with respect to drought [1].
Artemisia frigida Willd.(Asteraceae) has significant cenotic value as dominant or codominant in the composition of the cryophyte, true and desert steppes of Siberia, the Volga region, Kazakhstan, Central Asia, Mongolia and North America.The high plasticity of the species, as well as its ability of rapid vegetative reproduction through the development of creeping and quickly rooting aboveground shoots, determine the competitiveness of A. frigida plants when pasture degradation is intensifying [2].Artemisia frigida plants are heat-resistant, being able to survive in environments comparable to habitats of typical desert plants with temperatures reaching 56-57 • C [3].
In conditions of high air temperatures and aridity of the territory, volatile terpenoid compounds act as the first line of defense for the photosynthetic apparatus.Thus, the fumigation of plants with monoterpenes increases their thermal resistance.Monoterpenes, as well as many volatile sesquiterpenes found in plants, quickly react with reactive oxygen species (ROS) [4], thereby protecting the membranes of plant cells from oxidative damage.Sesquiterpenes can mitigate damage even more effectively than isoprene and monoterpenes.For example, the sesquiterpene (E)-β-caryophyllene is 43 times more reactive with ozone than the monoterpene limonene.Two other sesquiterpenes, (E)-β-farnesene and (E)-α-bergamotene, have a high potential for protecting against abiotic stress [5].The protective mechanism is associated both with direct reactions of terpenoids with oxidants, either intracellularly or at the boundary between the leaf and the atmosphere/ boundary layer, stabilizing the membrane, and with indirect changes in ROS signal transduction [6].Plants emit an enormous number of terpenoids annually, reaching huge quantities (500 teragrams) [7].Climate changes affect terpenoid emissions and have significant consequences for the structure and functioning of the biosphere [8].In addition to high temperatures, terpenoid emissions can also be induced by herbivore feeding, pathogen attack, radiation, salt stress and other abiotic stresses [9,10].
Terpenoid compounds are obtained from plants in the form of complex mixtures known as essential oils (EOs).EOs possess antibacterial, antifungal, herbicidal and other properties, and therefore have found wide applications in the pharmaceutical, food, agricultural, cosmetic and healthcare industries [11,12].EOs are complex mixtures of both native components and artifact compounds that inevitably form because of such processes as storage, drying of plant material, extraction and analysis.Their composition is related to many variables, such as cultivar genotype [13], the geographical area of production [14], harvest year [15], and extraction system [16].When obtaining EOs under the same conditions of storage and drying of raw materials, and using the same method of extraction and analysis of EOs under strictly defined parameters, we can observe the influence of genetic and epigenetic factors on the composition.
The aim of this study is to determine the patterns of the formation of the composition of the essential oil of Artemisia frigida Willd.plants, and to give possible reasons for its conservatism and lability.The isolated EOs were light, mobile liquids of a pale yellow color with a yield ranging from 0.3 to 1.0% (v/w).A total of 179 compounds were identified in the samples of the essential oil, with the majority being monoterpenoids and sesquiterpenoids (Supplementary Materials, Table S1).

Results and Discussion
For data analysis convenience, the entire variety of identified compounds were divided into four groups, based on their frequency of occurrence in the 27 tested samples: (1)-components found in 75% or more of the samples; (2)-in 50-74% of the samples; (3)-in 26-49% of the samples; and (4)-in 25% or less of the samples.
The second group consisted of 19 compounds.Such components as para-trans-menth-2-en-1-ol, myrtenol, bicyclogermacrene, and silphiperfol-6-en-5-one were identified in 20 out of 27 samples.The remaining 15 compounds were found in 14-19 samples.The total mass fraction of compounds of the second group was 0.96-29.34%.The content of each of the components, with the exception of myrtenol, selinene, beta-and isocaryophyllene in individual samples, did not exceed 2-3%.
The fourth group was the most numerous (101 components).However, fifty compounds were found only in one of the samples, and only five (nonanal, nerol, transchrysanthenol, trans-sabinene hydrate, 7-epi-silphiperfol-5-ene) were identified in six samples.The total proportion of components in this group varied from 0.72% to 59.93%.The mass fraction of most compounds in this group did not exceed 1-2%, except for two cases: the proportion of α-terpineol acetate in the sample "BURi 30.08.2017" was 6.30%, and that of α-thujone in the sample "OLNs 05.08.2017" was 10.52%.
The differences in the composition of the dominant components in the samples "CHN, 2016 4333 m" and "BUR 1343 m 17.08.2016"can be explained by the growth of plants in areas of high altitudes.The composition of the sample "BURroad 19.08.2021" components may be related to the fact that plants grew close to the federal road R258 "Baikal", with intensive traffic.The reasons for the differences in the composition of EOs in other samples are currently unclear.In any case, it is obvious that the composition of components and the content of their EOs are changeable under the influence of different factors.In addition, let us consider some of the factors that influence the composition of EOs of Artemisia frigida in more detail.

Composition of Essential Oils of A. frigida during Ontogenesis and Phenological Development of Plants
It has been established that changes occur in the biosynthesis of plants during ontogenesis in A. frigida.Indeed, it has been shown with potted crops that during the first year of plant life, changes occur in nitrogen metabolism.Specifically, the nitrogen content in shoots and roots decreases throughout the vegetative period [17].However, the genetic basis of these changes in biosynthesis during ontogenesis is currently unknown.An interesting example of this is the encoding of the biosynthesis of a large number of sesquiterpenes in patchouli plants with only five genes, which may be indirect evidence of the importance of epigenetic mechanisms in the biosynthesis of terpenoid compounds in plants [18].The first step towards differentiating the genetic and epigenetic nature of the components of EOs may be the study of their dynamics during plant development, as it is known that the quantity and composition of EOs change not only in space but also during plant ontogeny [19,20].Unfortunately, there are no data on the dynamics of the qualitative and quantitative composition of the essential oil at all stages of A. frigida ontogeny.Only the composition of the essential oil in the latent and generative phases of age development has been studied.
The latent phase of plant ontogeny is represented by fruits.In A. frigida, the fruit seeds contain EOs that act as anti-attractants for herbivore feeding, inhibitors of seed germination, sprout growth and reproduction of other species [21,22].The same components identified in the EOs obtained from the aboveground parts are also found in the composition of the fruit EOs.The main components of the fruit EOs are compounds with allelopathic activity: camphor (22.58%), borneol (22.46%), 1,8-cineol (25.11%), terpinene-4 (4.79%) and bornyl acetate (3.57%).The content of the remaining components does not exceed 1.5% (Table 1).Interestingly, unlike the aboveground parts of adult A. frigida plants in the generative phase of ontogeny, most of the oil in the fruits consists of monoterpenoids, which are characterized by greater volatility (91.33%).In our opinion, monoterpenoids facilitate the survival of the fruits, including protection against pests such as invertebrates (insects and others) and fungi and bacteria, as more volatile compounds interact more quickly with the surrounding environment during the early stages of seed germination and require less energy for their biosynthesis.The composition of the essential oil changes during different phases of plant phenological development.For example, in the study [23], it was shown that as the Lavandula × intermedia flower develops, the camphor content in the essential oil remains unchanged, but the concentration of borneol (a precursor to camphor) in the tissues increases.However, nowadays there is no direct molecular link between plant development and terpene biosynthesis that has been revealed [18].
The compositions of the EOs of A. frigida plants collected in the same populations but in different phenological phases were studied.The EOs were isolated from plants collected in three populations of the Ivolginsky district of the Republic of Buryatia in 2021: steppe ("BURi 27.07.2021","BURi 20.08.2021", "BURi 17.09.2021"),roadside ("BURroad 19.08.2021", "BURroad 17.09.2021")and hilltop ("BURhill 27.07.2021","BURhill 19.08.2021").The plants collected in July were in the budding phase, in August the flowering phase, and in September the fruiting phase.
During the transition from one phase to another, the yield of essential oils did not change, but the composition slightly changed.Compounds in the first group (except β-pinene, β-myrcene) were found in all stages of plant development, and only their proportion in the essential oil changed.However, there was no specific trend of decrease or increase in the transition from one phenophase to another, meaning that the direction of change was different for each population.Monoterpenoids such as β-pinene and βmyrcene were identified in the flowering ("BURroad 19.08.2021", "BURhill 19.08.2021") and fruiting phases ("BURroad 17.09.2021","BURi 17.09.2021").Compounds of the second group were also found in all phenophases, except myrtenyl acetate.The latter was only identified in the sample "BURi 20.08.2021".This sample also, unexpectedly, had a high content of the monoterpenoid alcohol myrtenol (4.27%).In the sample "BURhill 27.07.2021",silphiperfol-6-en-5-one (4.43%) and isocaryophyllene (7.11%) differed significantly in content compared to other samples, and the latter was also found in significant amounts in the sample "BURroad 19.08.2021" (10.12%).Compounds in the third group can be found in all phenophases (β-elemene, pentadecanal, etc.) as well as in specific phenophases (for example, β-copaene in the budding phase, and cabreuva oxide B in budding and flowering).Some of them (silphiperfol-6-ene, aristolochene, 4,5-di-epi-, teointermedeol, β-elemene, etc.) may have a content higher than 1%, reaching, for example, 5.33% for β-elemene.The majority of components in the fourth group were found in only one sample.Mintoxide, trans-sabinene hydrate and eugenol were identified in all phenophases.In the budding and flowering phases, 1,8-, 2,3-dehydro-cineol, and chrysanthenone were detected, and in the flowering and fruiting phases salviadienol was found.Davanafuran isomer and silphiperfola-4,7(14)-diene were detected in the budding and fruiting phases (Supplementary Materials, Table S3).
The dynamics of the content of major components (with a content of 3% and more) of the essential oil from budding to fruiting was considered using a specific steppe population as an example.The main components included compounds from the first group: bornyl acetate, terpinen-4-ol, α-terpineol, camphor, spathulenol, 1,8-cineol and borneol, as well as one component each from the second (myrtenol) and third (neointermedeol) groups.The changes in the content of camphor, spathulenol and borneol were insignificant.However, the proportion of terpinen-4-ol decreased by half during the transition from budding to flowering, while α-terpineol increased 2.6 times during fruiting compared to other phenophases.The changes in the content of components from the labile part of the A. frigida essential oil were more significant.For example, the content of myrtenol increased by 8.2-13.3times during the flowering phase.Neointermedeol was not detected during the flowering phase, while its proportion was approximately the same during the budding and fruiting phases (3.30-3.61%)(Figure 1).detected during the flowering phase, while its proportion was approximately the same during the budding and fruiting phases (3.30-3.61%)(Figure 1).As for the trends in the content of individual compounds in different phenophases, four patterns can be identified, three monomodal (Nos 1-3) and one bimodal (No 4): (1) the content of the compound is highest during the budding phase and then gradually decreases or remains at the same level (camphor, spathulenol, and terpinen-4-ol); (2) the content reaches its peak during the flowering phase (bornyl acetate, 1,8-cineol, borneol, and myrtenol); (3) the content of the compound is highest during the fruiting phase (terpineol, α-); and (4) the compound is present during the budding and fruiting phases but absent during flowering (neointermedeol).The presence of different patterns in the essential oil content during phenological development demonstrates the complexity of adaptive mechanisms.
Thus, the composition of the essential oil changed during plant development, particularly during phenological development.Compounds from the first and second groups were found in almost all phenophases, but their content in the essential oil varied.They formed a constant part of the essential oil, but no trends in the dynamics of their content were observed.It can be assumed that the formation of compounds in the constant part of the essential oil is predetermined at the genetic level, and thus this set is an immanent property of the species.On the other hand, differences in the composition of compounds from the third and fourth groups in different phases of plant development can be explained by the formation of new functional organs in which additional substances are synthesized [24].Their formation might be associated with the species' biosynthetic reaction to certain conditions, which accordingly has an epigenetic nature, as well as with the transformations of native compounds during the processes of drying, extraction and analysis.

Composition of Essential Oils of A. frigida from the Same Population Depends on the Weather Conditions in Different Years
The differences in the compositions of EOs obtained from plants of the same population at different phenophases are balanced compared to the differences between the compositions of EOs of the same population in different years, i.e., under different weather conditions.Thus, on the biplot (PC1-PC2), the projections of samples of EOs were divided As for the trends in the content of individual compounds in different phenophases, four patterns can be identified, three monomodal (Nos 1-3) and one bimodal (No 4): (1) the content of the compound is highest during the budding phase and then gradually decreases or remains at the same level (camphor, spathulenol, and terpinen-4-ol); (2) the content reaches its peak during the flowering phase (bornyl acetate, 1,8-cineol, borneol, and myrtenol); (3) the content of the compound is highest during the fruiting phase (terpineol, α-); and (4) the compound is present during the budding and fruiting phases but absent during flowering (neointermedeol).The presence of different patterns in the essential oil content during phenological development demonstrates the complexity of adaptive mechanisms.
Thus, the composition of the essential oil changed during plant development, particularly during phenological development.Compounds from the first and second groups were found in almost all phenophases, but their content in the essential oil varied.They formed a constant part of the essential oil, but no trends in the dynamics of their content were observed.It can be assumed that the formation of compounds in the constant part of the essential oil is predetermined at the genetic level, and thus this set is an immanent property of the species.On the other hand, differences in the composition of compounds from the third and fourth groups in different phases of plant development can be explained by the formation of new functional organs in which additional substances are synthesized [24].Their formation might be associated with the species' biosynthetic reaction to certain conditions, which accordingly has an epigenetic nature, as well as with the transformations of native compounds during the processes of drying, extraction and analysis.

Composition of Essential Oils of A. frigida from the Same Population Depends on the Weather Conditions in Different Years
The differences in the compositions of EOs obtained from plants of the same population at different phenophases are balanced compared to the differences between the compositions of EOs of the same population in different years, i.e., under different weather conditions.Thus, on the biplot (PC1-PC2), the projections of samples of EOs were divided into two parts.One part was grouped, regardless of the different phenophases of the samples, of samples from 2021, while the other part consisted of samples from other years (Figure 2).Indeed, the year 2021 significantly differed from others in terms of temperature and moisture conditions (hydrothermal coefficient-HTC 1.6), being more humid.In contrast, the years 2014-2020 were drier (HTC 0.3-0.9).The distance between samples collected in June and July 2015 ("15-9", "15-18") and July and August 2017 ("BURi 01.07.2017", "BURi 30.08.2017") can be explained by differences in climatic parameters from the beginning of vegetation to the date of collection.For example, in 2015, the average daily temperature from May to the date of plant collection (average t veg. ) for the "15-9" sample was 13.8 • C, while for "15-18" it was 16.1 • C. The cumulative rainfall in May until the date of plant collection (∑R veg ) for the "15-9" sample was 29.1 mm, while for "15-18" it was 62.5 mm.Correspondingly, the extremeness coefficients (K extr.veg. ) differed almost twofold (0.5 and 0.3, respectively).In contrast, in 2021, temperature fluctuations and precipitation were more evenly distributed; climatic parameters (average t veg ., ∑Rveg) did not differ significantly at different collection times, and the extremeness coefficient (K extr.veg. ) from budding to fruiting had a value of 0.1 (Supplementary Materials, Table S4).
into two parts.One part was grouped, regardless of the different phenophases of the samples, of samples from 2021, while the other part consisted of samples from other years (Figure 2).Indeed, the year 2021 significantly differed from others in terms of temperature and moisture conditions (hydrothermal coefficient-HTC 1.6), being more humid.In contrast, the years 2014-2020 were drier (HTC 0.3-0.9).The distance between samples collected in June and July 2015 ("15-9", "15-18") and July and August 2017 ("BURi 01.07.2017", "BURi 30.08.2017") can be explained by differences in climatic parameters from the beginning of vegetation to the date of collection.For example, in 2015, the average daily temperature from May to the date of plant collection (average tveg.) for the "15-9" sample was 13.8 °C, while for "15-18" it was 16.1 °C.The cumulative rainfall in May until the date of plant collection (∑Rveg) for the "15-9" sample was 29.1 mm, while for "15-18" it was 62.5 mm.Correspondingly, the extremeness coefficients (Kextr.veg.)differed almost twofold (0.5 and 0.3, respectively).In contrast, in 2021, temperature fluctuations and precipitation were more evenly distributed; climatic parameters (average tveg., ∑Rveg) did not differ significantly at different collection times, and the extremeness coefficient (Kextr.veg.) from budding to fruiting had a value of 0.1 (Supplementary Materials, Table S4).Correlation analysis showed a high (correlation coefficient values from 0.7 to 0.9) and very high (from 0.9 to 1.0) degree of association between the values of climatic parameters and some components of the essential oil.For example, the correlation coefficient values (0.92-0.98) between HTC values and the content of spathulenol, caryophyllene oxide and silphiperfol-6-en-5-one indicated a very high positive correlation between them.The values of average t veg .showed a high positive correlation with the content of borneol and a negative correlation with silphiperfol-6-en-5-one.The values of ∑Rveg positively correlate with the content of spathulenol and caryophyllene oxide, while K extr.veg.correlated positively with αand γterpinene (Supplementary Materials, Table S5).Thus, the mass content of compounds from the first and second groups in the EOs of plants from the same population depends on the weather conditions in the year of raw material collection.

Composition of Essential Oils of A. frigida from Different Populations in Buryatia Depends Not Only on Climatic Conditions
PCA analysis was conducted on data obtained from other populations (in Selenginsky and Eravninsky districts) and compared with the EOs of plants growing in the Ivolginsky district.The differences in the composition of EOs of plants collected in different years were less pronounced than the differences between the compositions of EOs of plants from different locations (Figure 3).Correlation analysis showed a high (correlation coefficient values from 0.7 to 0.9) and very high (from 0.9 to 1.0) degree of association between the values of climatic parameters and some components of the essential oil.For example, the correlation coefficient values (0.92-0.98) between HTC values and the content of spathulenol, caryophyllene oxide and silphiperfol-6-en-5-one indicated a very high positive correlation between them.The values of average tveg.showed a high positive correlation with the content of borneol and a negative correlation with silphiperfol-6-en-5-one.The values of ∑Rveg positively correlate with the content of spathulenol and caryophyllene oxide, while Kextr.veg.correlated positively with α-and γ-terpinene (Supplementary Materials, Table S5).Thus, the mass content of compounds from the first and second groups in the EOs of plants from the same population depends on the weather conditions in the year of raw material collection.

Composition of Essential Oils of A. frigida from Different Populations in Buryatia Depends Not Only on Climatic Conditions
PCA analysis was conducted on data obtained from other populations (in Selenginsky and Eravninsky districts) and compared with the EOs of plants growing in the Ivolginsky district.The differences in the composition of EOs of plants collected in different years were less pronounced than the differences between the compositions of EOs of plants from different locations (Figure 3).The samples collected in the Eravninsky district differed the most.They were characterized by relatively low content of camphor (up to 5.5%) and noticeable amounts of certain sesquiterpenoids (β-bisabolene, caryophyllene, caryophyllene oxide, and spathulenol).In contrast, the composition of EOs from plants in the Ivolginsky and Selenginsky districts differed to a lesser extent.EOs of samples from the Selenginsky district contain higher amounts of cymol and para-and bornyl acetate compared to samples from the Ivolginsky district.
Linear regression models show that only the content of trans-carveol, para-, transmenth-2-en-1-ol, camphor and bornyl acetate in the EOs has been changed under the influence of climatic parameters.For example, the content of trans-carveol was influenced by the climatic extremeness of the year.Among the dominant components, the content of camphor in the essential oil linked to changes in climatic parameters (Table 2).The obtained results showed that the influence of climatic factors on the composition of EOs in plants was varying, which was also shown for plants of the genus Nepeta [26] and Lavandula [27].Significant differences in the composition of EOs of plants from the Eravninsky district compared to the Selenginsky and Ivolginsky districts may be related to the fact that negative values of the average annual temperature lead to the existence of permafrost in the soils of the Eravninsky district, where permafrost has discontinuous distribution.In contrast, the Selenginsky and Ivolginsky districts belong to the zone of the only isolated distribution of permafrost.Thus, plants in the Eravninsky district experience a greater influence of low soil temperatures, which in turn is reflected in the composition of EOs.There is information in the literature about the influence of soil type [28] and substrate composition [29], but there are no data on the effect of soil temperature or permafrost on the composition of the essential oil, although the soil temperature affects the content of secondary metabolites in the plant.For example, a higher maximum soil temperature resulted in lower phenol content, while a higher minimum soil temperature and shorter harvest age increased the total flavonoids in Java tea (Orthosiphon aristatus B.) [30].
Previously, it was shown for Artemisia vulgaris L. [41], A. annua L. [42], and A. jacutica Drob.[43] that under the influence of macroclimatic conditions, distinct "western" and "eastern" chemotypes of EOs were formed, which differed in their component composition.In the case of A. frigida, chemotypic differentiation was not observed, as there was no clear discrepancy between the essential oil samples obtained from "western" and "eastern" populations.In general, according to the PCA analysis in the case of A. frigida, most samples, regardless of their place of growth, did not show significant differences in composition.Nevertheless, there were 13 components that were most noticeably linked to climatic conditions, which was manifested in their greatest contribution to the distribution of samples on the biplot (Figure 4).
It can be noted that the majority of "western" samples were characterized by a relatively high content of monoterpenoids-camphor, 1,8-cineol and camphene.And for most "eastern" samples, in addition to monoterpenoids-bornyl acetate, terpineol-4, αterpineol, and myrtenol-sesquiterpenoids such as spathulenol and caryophyllene oxide were presented in noticeable concentrations.However, even these features are not enough to allow for the identification of specific chemotypes.The reason may lie in the biological characteristics of A. frigida, which often grows on sandy and stony soils with relatively sparse vegetation cover [44].Meanwhile, A. frigida plants from Western and Eastern Siberia are in similar conditions in terms of heat and moisture supply in their habitats.Thus, the values of hydrothermal coefficient (HTC) calculated from meteorological data from Western and Eastern Siberia ranged from 0.5 to 1.4, i.e., from "very arid" to "humid" conditions in both regions (Table 3).However, 10 samples of EOs of plants growing in areas where climatic parameters differed from typical ones formed a single locus ("atypical" group, Figure 4).The differences were associated with the influence of high altitude, as well as the proximity of roads, or increased insolation.This locus was formed by plant samples from China, growing at altitudes of 2735-4333 m, Mongolia, growing at 1327 m and 1854 m, Buryatia, growing at 1343 m above sea level (a.s.l.), and plants growing near roads ("ALT1road", "BURroad 19.08.2021"), as well as a sample of plants, "BURss 24.08.2016",located on the southern slope of the Soldatsky Ridge in Buryatia.It can be noted that the majority of "western" samples were characterized by a relatively high content of monoterpenoids-camphor, 1,8-cineol and camphene.And for most "eastern" samples, in addition to monoterpenoids-bornyl acetate, terpineol-4, α-terpineol, and myrtenol-sesquiterpenoids such as spathulenol and caryophyllene oxide were If the differences in the composition of EOs in plants due to the influence of altitude (which determines climatic parameters) and the proximity of roads do not require any special explanations, then the influence of the slope exposure, which is the habitat of the plants, should be considered in more detail.Essential oil samples "OLN4", "BUR10" [25] and "BURss 24.08.2016" were obtained from plants collected on slopes of different exposures.In the sharply continental climate of semiarid regions, including Transbaikalia and the Cis-Baikal area, the ecological conditions on different slopes differ significantly.Southern slopes receive more heat during the vegetative period.On southern insolated slopes the evaporation is much stronger, and plant communities experience a moisture deficit [45][46][47][48].Thus, it is confirmed again that the composition of essential oil is determined by climatic parameters, in particular the ratio of heat to moisture.
Samples from plants growing near the road ("ALT1road", "BURroad 19.08.2021") were most similar in composition and content of dominant components to the samples from "typical" habitats.Comparative analysis of the composition of EOs from A. frigida in "typical" and "atypical" habitats showed that 22 compounds were identified in half or more of the samples from both groups.These included camphor, 1,8-cineol, germacrene D and others.Some compounds found in the majority of EOs from "typical" habitats were found in fewer samples from "atypical" habitats.Borneol was found in all samples from "typical" habitats and in only 4 out of 10 samples from "atypical" habitats.Myrtenol was found in only one of the samples from "atypical" habitats.On the other hand, (E)-β-farnesene was found in 7 out of 10 samples from "atypical" habitats, while it occurred in only 20% of the EOs from plants in "typical" habitats.Among the other 261 components identified in less than half of the samples, 109 were found only in the EOs from "typical" habitats, 27 were found in "atypical" habitats, and 112 were common to plants from both groups.For example, α-bisabolol was specific to 18% of the samples from "typical" habitats, and yomogi alcohol was found in 3 out of 10 samples from "atypical" habitats, while trans-piperitol was identified in the EOs of both groups.

Conservatism s and Lability of the Composition of A. frigida Essential Oils
Based on the abovementioned, it was possible to distinguish the "core" of the A. frigida essential oil composition, which referred to those components that were formed as part of the genetic programs of terpenoid biosynthesis, regardless of the plant's habitat.A total of 24 such components formed the "core".However, the content of some of these components in the essential oil could vary during plant development and under the influence of climatic conditions.The total content of "core" components in the EOs of plants from "typical" habitats ranged from 48.31% to 92.60%, with an average value of 75.18%, while in the EOs of plants from "atypical" habitats it ranged from 4.82% to 77.17%, with an average value of 48.43%.The fluctuations in the content of camphor, 1,8-cineol, borneol and germacrene D were particularly noticeable.
Surrounding the "core" was a "cloud" of components whose formation was influenced by external factors.This "cloud" consisted of two concenters.The first concentric circle (11 compounds) included compounds that have been identified in samples from both groups but were found in the majority of essential oil samples from either "typical" or "atypical" habitats.The second concentric circle (248 compounds) comprised compounds that were found in a smaller number of samples from both groups (Figure 5).In addition to terpenoid compounds, the second concentric circle included phenylpropenes (eugenol, etc.) and acyclic compounds (capillene, etc.).The conservatism of the qualitative composition of the "core" was largely related to the characteristics of the compounds' biosynthesis.For example, one of the dominant components of A. frigida EOs, 1,8-cineol, was always accompanied by myrcene, limonene, sabinene, α-, β-pinene, α-terpineol [49], terpinolene, (E)-β-ocimene, and a-thujene [50,51] The conservatism of the qualitative composition of the "core" was largely related to the characteristics of the compounds' biosynthesis.For example, one of the dominant components of A. frigida EOs, 1,8-cineol, was always accompanied by myrcene, limonene, sabinene, α-, β-pinene, α-terpineol [49], terpinolene, (E)-β-ocimene, and a-thujene [50,51] in small quantities, because they are formed under the action of the multiproduct enzyme cineole synthase.These components are known as "cineole cassette".
Borneol and camphene are formed from a common precursor, bornyl diphosphate, and are catalyzed by the multiproduct enzyme BOR-synthases [52].Other compounds, such as borneol and camphor, are interconnected.The intermediate bornyl diphosphate is hydrolyzed to borneol through the catalytic activity of bornyl diphosphate synthase and borneol synthase, and then borneol is oxidized to form camphor through the catalytic activity of bornyl dehydrogenase.Previously it was considered that the enzymatic reaction catalyzed by γ-terpinene synthase lead to the formation of γ-terpinene, which further transformed into p-cymene and the phenolic terpene isomers thymol and carvacrol [53].Now it is established that γ-terpinene is oxidized by monooxygenases of the cytochrome P450 subfamily CYP71D to form unstable intermediate cyclohexadienols, which are then dehydrogenated by short-chain dehydrogenase/reductase (SDR) to the corresponding thymol and carvacrol.The formation of p-cymene occurs due to spontaneous rearrangement of the dienolic intermediate products because of their instability under aqueous conditions [54].
According to recent studies, both 2E,6E-farnesyl pyrophosphate (E,E-FPP) and its isomer Z,Z-FPP are substrates for sesquiterpene synthases [55].Acyclic sesquiterpenoids, such as (E)-β-farnesene and (E)-nerolidol, are formed through farnesyl and nerolidol cation, respectively.The biosynthesis of cyclic sesquiterpenoids is considered in terms of biogenetic pathways.In the essential oil of A. frigida, three groups (chemotypes) of sesquiterpenoids can be distinguished: germacrane, bisabolane and humulane sesquiterpenoids [56], with a predominance of germacrane compounds (Supplementary Materials, Table S7).The main compound of the germacrane chemotype in the EOs of A. frigida is germacrene D, which serves as a biogenetic precursor for many other sesquiterpene structures [57], including α-copaene, bicyclogermacrene, β-selinene and β-elemene.Spathulenol is formed from bicyclogermacrene during the extraction, storage and analysis of the essential oil [58].Humulene, caryophyllene, caryophyllene oxide and triquinanes are representatives of the humulane series.All of these compounds were found in most samples and were part of the "core" and "concentric circle 1".The exception was the group of triquinanes, which were found in specific samples.Among them, silphiperfol-6-en-5-one was the most frequently identified compound (48% of "typical" samples, and 3 out of 10 "atypical" samples).Bisabolane sesquiterpenoids were not represented in the "core", and β-bisabolene was a part of "concentric circle 1".The rest were included in "concentric circle 2", but their content in the EOs of high-altitude plants was quite high.The total content of bisabolane sesquiterpenoids in samples from Qinghai province (China) ranged from 8.70% to 45.91%."Concentric circle 2" included compounds formed as a result of further modification of terpenoids from the "core" and "concentric circle 1".The presence of irregular monoterpenoids was observed in "concentric circle 2", and they were formed through condensation of C5 units in a "non-head-to-tail orientation", such as yomogi alcohol, lavandulol, artemisia alcohol and artemisia ketone.
Mono-and sesquiterpenoid compounds are multifunctional compounds.Their role in biotic interactions, such as allelochemicals, allelopathic compounds, antifeedants, attractants, and repellents, and their importance in plant defense against bacterial pathogens and fungi, has been well studied.At the same time, it is known that their emission is closely related to changes in abiotic factors in their environment.The existence of special biosynthetic pathways and the preservation of their emission, despite their "expensive" cost for plants (1-2% of photosynthetic carbon fixation) has led to the proposal of a model of a "unified biochemical mechanism for multiple physiological stressors".It has been shown that volatile isoprenoids participate in direct reactions with oxidants and can also modify the transfer of ROS signals and/or contribute to membrane stabilization [4].For example, α-bisabolol demonstrates high antioxidant activity in vitro [59].Overall, the role of isoprene in plant thermotolerance has been extensively studied, while the role of volatile terpenoids has been less explored [60].It has been shown that terpenes can protect plants from heat damage, and there is a correlation between monoterpene emission and heat stress resistance [61].In our opinion, the emission of terpenoids by A. frigida explains its thermotolerance.The remarkable fact is that the "core" predominantly consists of compounds (pinene, 1,8-cineol, camphor, spatulenol, etc.) with pronounced allelopathic activity [20].At the same time, β-caryophyllene is involved in attracting natural enemies of herbivorous insects [62] and in the defense against microbial pathogens [63].
The components of "concentric circle 1" can act as allelochemicals (para-menth-2-en-1ol) or can have other activities.For example, (E)-β-farnesene on one hand is secreted by aphids as a signaling pheromone, and on the other hand it is produced by infected plants to attract natural enemies of aphids [64].
On the other hand, the composition of compounds in "concentric circle 2" is determined by the conditions of specific plant habitats, including abiotic factors and interactions with other plants, microorganisms, insects and herbivores, and even humans.Pathogen attack can shift biosynthesis towards the production of sesquiterpenoid phytoalexins [65].Triquinanes have antibiotic properties [66].The increase in the content of copaene in EOs, an intermediate compound in the pathway of terpenoid alkaloids' biosynthesis, may also be associated with protection against mammals [67].
Changes in the content of compounds in EOs are associated with their emission.In high-altitude conditions, for example, there seems to be a faster conversion of borneol into camphor, resulting in a tendency for a decrease in borneol content and an increase in camphor content in the essential oil.Since we extract a pool of accumulated terpenoids from the tissues in the form of essential oil, we cannot definitively say whether there is an enhancement of camphor and borneol emission or if their biosynthesis is blocked at very high altitudes (4333 m above sea level).It seems that the biosynthetic pathways of the "core" components are maintained regardless of plant habitat conditions, while the rate of their emission changes.Specific factors, often biotic, influencing plants, affect the components of "concentric circle 1 and 2".
Thus, mono-and sesquiterpenoids serve as the first line of defense for plants against biotic and abiotic stresses.They are essential for plant survival, and their immense diversity is shaped by factors that have both a genetic and epigenetic nature.Overall, we can observe a notable change in the compositional profile of A. frigida EOs under the influence of climatic parameters.

Plant Materials
The plant material was collected in the natural steppe habitats with chestnut soils within the territory of Russia (Republic of Buryatia) from 2016 to 2021 (Table 4).The collected raw material was dried to an air-dry state (2 weeks).Herbarium specimens are stored in the Laboratory of the Chemistry of Natural Systems at Baikal Institute of Nature Management, Siberian Branch, Russian Academy of Sciences.

Isolation of Essential Oils
A total of 30~40 g of the powdered plant material was subjected to hydrodistillation for 3 h using a Clevenger-type collector apparatus.The essential oils extracted in this research were obtained from the aerial parts of the plant.All experiments were conducted in triplicate, and the results were reported based on the air-dry weight.

GC-MS Analysis of Essential Oils
The composition of the EOs was determined using gas chromatography-mass spectrometry (GC-MS) analysis after obtaining EOs.An Agilent 6890 gas chromatograph (Agilent Technologies, Santa Clara, CA, USA) equipped with an HP 5973 quadrupole mass selective detector (Hewlett-Packard, Palo Alto, CA, USA) and an HP-5MS capillary column (30 m × 0.25 mm × 0.2 µm; Hewlett-Packard) was used.The analysis was conducted with an electron impact at 70 eV and with the electron multiplier set at 2200 V. Helium was used as the carrier gas at a flow rate of 1 mL/min.The temperature of the column was programmed to increase from 50 to 240 • C at a rate of 4 • C/min, with an initial hold of 5 min, and then from 240 to 280 • C at a rate of 20 • C/min, with a final hold of 5 min.The injector and detector temperatures were set at 280 and 250 • C, respectively.The purity of the helium gas used was 99.999%.The column pressure was set to 52.8 × 10 3 Pa, and the split ratio was 60:1.

Essential Oil Compound Identification
The identification of essential oil constituents was conducted by comparing their mass spectra with standards from the NIST and with calculated linear retention indices (RIs) found in the literature [58].The RIs were obtained by co-injecting a mixture of linear C8-C20 hydrocarbons.Identification was considered valid when the computer matching with the mass spectral libraries had a probability of over 90%.The contents of the oil components were determined by their relative amounts (%) and expressed as a percentage peak area relative to the total peak area, based on the GC-MS analyses of the oils.

Climatic Conditions
To characterize the main climatic parameters, data provided by the All-Russia Research Institute of Hydrometeorological Information (World Data Center) were used.The downloaded data included air temperatures, daily precipitation amounts, monthly precipitation sums, monthly air temperatures, number of days with precipitation equal to or greater than 1 mm, average monthly relative air humidity, monthly sunlight hours and more.Additionally, integrated indicators such as Selyaninov's hydrothermal coefficient (HTC), temperature-humidity extremeness coefficient (Cextr) and the Standardized Precipitation Evapotranspiration Index (SPEI) were calculated.
The HTC was estimated to assess the levels of heat and humidity in specific regions using the formula HTC = R × 10/Σt.In this equation, R represents the precipitation sum in millimeters for June-July at mean daily temperatures above +10 • C, and Σt denotes the sum of mean daily temperatures in degrees Celsius ( • C) for the same period.
The C extr coefficient was calculated to compare the plant habitat conditions across different years.This coefficient represents the ratio of the mean temperature value for a specific period to the total precipitation for the same period, using the formula C extr = t • C/mm.Several parameters were estimated, including the sum of daily temperatures (Σt veg ) and the sum of precipitation (ΣR veg ) from the start of the vegetative phase (May) until the plant harvesting.Additionally, the sum of daily temperatures (Σt year ) and the sum of precipitation (ΣR year ) from August of the previous year to August of the year of plant harvesting were calculated.The temperature-humidity extremeness coefficients were estimated for both the period from the start of the vegetative phase until plant harvesting (C extr.veg ), and for the entire year (C extr.year ) .To simplify the calculations, the resulting values were multiplied by 100 and 1000, respectively, and rounded to whole numbers.
The Standardized Precipitation Evapotranspiration Index (SPEI) was calculated using the SPEI package.This index was estimated for both the calendar year (January-December) (SPEI cal ) and the period from August of the previous year until August of the year of plant harvesting (SPEI year ).

Statistical Analysis
To determine the influence of the main climatic parameters (air temperature, atmospheric precipitation, etc.) and their integral characteristics on the composition of plant essential oil, two types of linear regression models were constructed [68].The first type of model was represented by the equation y = a*x + b.This type of model assumed that the content of a certain component of the essential oil, y, was linearly dependent on a certain parameter, x.The equation for the second type of model was z = a*x + b*y + c.Similar to the first case, a linear dependence was assumed between the content of a certain component of the essential oil, y, but this time on two parameters.The modeling was performed using the R programming language interpreter version 4.2.1 for the x86_64-pc-linux-gnu platform.The linear regression models were built using built-in R functions such as lm function and operations on data structures of the data.frametype.The Gauss-Markov hypothesis testing was conducted using the gvlma function from the "gvlma" package [69].Correlation analysis was performed using the Excel2016 software package.For visualization purposes the data of the component composition of the essential oil were processed

2 . 1 .
Composition of Essential Oils of Artemisia frigida (Based on Plants Growing within the Territory of Buryatia (Russia) and Qinghai (China) from 2016-2021

Figure 1 .
Figure 1.Content of major components of EOs of A. frigida plants during budding phase (green), flowering phase (blue) and fruiting phase (yellow).

Figure 1 .
Figure 1.Content of major components of EOs of A. frigida plants during budding phase (green), flowering phase (blue) and fruiting phase (yellow).

Figure 3 .
Figure 3. Biplot (PC1-PC2) of multi-year data on the content of constant compounds in the EOs of Artemisia frigida Willd.plants collected in the Ivolginsky, Selenginsky, and Eravninsky districts of Buryatia (Principal Component Analysis).Sample codes are displayed in red color, districtsblue circles.

Figure 4 .
Figure 4. Biplot (PC1-PC2) of data on the content of constant compounds in the EOs of Artemisia frigida Willd.plants from different regions of the world (Principal Component Analysis).Designation "abbreviation and serial number," "abbreviation and code," "abbreviation, year"-published data, including own data from Russia (Republic of Buryatia, Irkutsk region) and Mongolia.Only designated with a code and collection date-data of 2016-2021.

Figure 4 .
Figure 4. Biplot (PC1-PC2) of data on the content of constant compounds in the EOs of Artemisia frigida Willd.plants from different regions of the world (Principal Component Analysis).Designation "abbreviation and serial number", "abbreviation and code", "abbreviation, year"-published data, including own data from Russia (Republic of Buryatia, Irkutsk region) and Mongolia.Only designated with a code and collection date-data of 2016-2021.

Figure 5 .
Figure 5. Conservative and labile components of essential oils of A. frigida.

Table 2 .
Content of components in the essential oil of Artemisia frigida Willd.plants growing in the Ivolginsky, Selenginsky and Eravninsky districts of the Republic of Buryatia, according to linear regression models.

Table 4 .
Collection information for Artemisia frigida used in this study for the years 2016-2021.