The influence of environmental condition on the creation of organic compounds in Pinus sylvestris L. rhizosphere, roots and needles

Studied organic molecules in Pinus sylvestris L. seem to have acted as a safety net for metal transport, chelation and sequestration, allowing adaptation and growth under highly polluted conditions. Pinus sylvestris L. is known for its ability to survive in areas of highly elevated metal pollution, such as flotation tailings. The aim of the study was to estimate the content of selected organic molecules (including aliphatic low molecular weight organic acids (ALMWOAs), phenolic compounds and terpenes) and the physiological mechanisms underlying differences in metal/metalloid tolerance of P. sylvestris growing in unpolluted (soil) and polluted (flotation tailings) areas. The dominant ALMWOAs in rhizosphere soil extracts were citric acid followed by malic and oxalic acids, whereas in flotation tailings malic and oxalic acids. In roots and needles, the content of ALMOWAs was significantly higher in P. sylvestris L. tissue from flotation tailings in comparison to soil. Phenolic compounds were detected only in roots and needles, with a generally higher content of nearly all detected compounds from flotation tailings. The composition of roots did not contain all the compounds detected in needles. The profile of needles additionally contained four hydroxybenzoic, protocatechuic and salicylic acids. In pine needles, 24 volatile terpenes were identified in total. The content of these compounds in pine needles from the polluted area was markedly different from the unpolluted area. The dominant volatile monoterpenes in P. sylvestris L. needles from the unpolluted area was three carene, while in pine needles from the polluted area monoterpenes α-pinene was dominant.


Introduction
Soil contamination with toxic trace elements is widespread and has a significant impact on such specific ecosystem functions as found in soil (Huang et al. 2016). It is particularly Communicated by E. Magel .
* Zuzanna Magdziak zuzanna.magdziak@up.poznan.pl; zuza@positiv.pl associated with anthropogenic activities which currently pose a serious problem for ecological equilibrium, adversely affecting the maintenance of environments and proving destructive for flora and fauna (Rocha et al. 2016). Several species of woody plants, such as Salix and Populus (Tlustoš et al. 2007;Vamerali et al. 2009;, Quercus robur L. and Acer platanoides L. (Budzyńska et al. 2017a), Ulmus laevis Pall (Budzyńska et al. 2017b) and Betula pendula Roth. (Mleczek et al. 2016), have been shown to have a high potential for phytoextraction, uptake and accumulation of metals present in the soil matrix and can be a suitable choice for dendroremediation of extremely polluted areas. Pinus sylvestris L. (P. sylvestris) is a common woody species that can be used for restoring degraded soil ecosystems because of its adaptive capabilities (Placek et al. 2016). According to , P. sylvestris may survive in areas with high levels of metal pollution, such as flotation tailings. In studied specimens of P. sylvestris originating from polluted areas a significantly lower rate of phytoextration of selected elements was recorded in comparison with unpolluted areas ). Nevertheless, the P. sylvestris was able to survive, grow and develop. Differences in biomass and phytoextraction ability clearly demonstrate the influence of environmental conditions of growth, but at the same time indicate the importance of physiological changes in the studied plants. In the light of the diverse potential of the studied plants and limited literature data on the creation and exudation by woody plant roots of different organic compounds (aliphatic low molecular weight organic acids (ALMWOAs), enzymes, amino acids, phenolic compounds, simple and complex sugars, vitamins, purines, proteins and flavonoids) into the rhizosphere (Ryan et al. 2001;Adeleke et al. 2017;Magdziak et al. 2017), we investigated the effects of flotation tailings in organic biomolecule activity in the rhizosphere. The rhizosphere is the most important soil zone, where fundamental processes occur that are responsible for plant functioning. Any change in the soil chemistry, including metal/metalloid concentration, can impact the cycling of carbon and other nutrients (Huang et al. 2016). The above molecules, especially ALM-WOAs and phenolic compounds, are essential factors for nutrient acquisition (Dinh et al. 2017). They are present in the most abundance and are most reactive with metals (Koo et al. 2010); they play a role in alleviation of anaerobic stress in roots as well as mineral weathering (Adeleke et al. 2017). They additionally influence several soil processes, e.g. sorption and desorption (Wang et al. 2015), oxidation and reduction (Blaylock and James 1994) and precipitation and dissolution (Zhou et al. 2007). Moreover, in some of the previously mentioned studies, ALMWOAs and phenolic compounds have been found in plant tissues (mainly, roots and leaves) (Drzewiecka et al. 2017;Magdziak et al. 2017), where next to the rhizosphere, plants use ALMWOAs to transport, sequestrate and prevent cytoplasmic precipitation of toxic elements in cells, or in the case of phenolic compounds, in participation in adaptation and detoxification mechanisms mainly related to their structure and antioxidant properties (Ivanov et al. 2012;Jiang et al. 2017a, b;Benbettaieb et al. 2018). Plant cells also produce many other groups of chemical compounds. Extractives include, e.g. essential oils, which are multicomponent mixtures of mono-, sesqui-and diterpene compounds or phenylpropane derivatives (phenolics). They take various forms: hydrocarbons, alcohols, aldehydes, ketones, esters or ethers. Furthermore, essential oils may contain sulfur and nitrogen substances of coumarins. The percentage shares of individual components vary and depend on many factors, i.e., plant ontogenesis, geographical region, growing and harvest conditions, storage method and preparation for further processing (Gonçalves et al. 2003;Silva et al. 2003;Dob et al. 2005;Bakkali et al. 2008;Gilles et al. 2010;Stefanakis et al. 2013).
It is known that the chemical composition of oils depends on many factors (Grochowski 1990;Głowacki 1994). Unfavorable growth conditions, such as soil or air pollutions lead to a biochemical tree response to stress (including changes in the composition of essential oils) (Chojnacki and Cichy 1995). However, up to now, no studies have been able to show a clear relationship between the composition of the oils and toxic metal at high content level in the substrate on which the trees grow (Kainulainen et al. 1992;Supuka and Berta 1998;Fuksman 2002).
More studies are necessary to understand the mechanisms triggered by woody plants. Among toxic metals described in literature data, cadmium (Cd), lead (Pb) and zinc (Zn) are often presented as metal pollution stress which activates defense mechanisms related to changes in the creation and exudation of organic molecules. On the other hand, there is limited data on the response of woody species to their exposure to highly toxic elements, such as arsenic (As), mercury (Hg) or thallium (Tl).
For this reason, the present study extends the research carried out by , where P. sylvestris trees were grown in soil and flotation tailings characterised by specific chemical properties. The aim of the study was to estimate the content of selected organic molecules and physiological mechanisms underlying differences in metal/metalloid tolerance in P. sylvestris growing on soil (unpolluted area) and extremely contaminated flotation tailings (polluted area). The investigation focused on an evaluation of changes in aliphatic ALMWOAs, phenolic compound content in the rhizosphere and/or roots, as well as ALMWOAs, phenolic compounds and terpenes in needles. Such a determination is essential to verify how organic molecules might influence an increase in bioavailability and the accumulation of metals/metalloids in P. sylvestris organs. The above molecules were chosen as biochemical parameters of plant reaction to elevated concentrations of toxic elements in flotation tailings due to their accumulation, chelation, regulation and translocation as a probable mechanism of their detoxication as well as an antioxidant function (Viehweger 2014;Dinh et al. 2017).

Characteristics of P. sylvestris L. specimens and area description
The experimental materials were five specimens of 9-yearold P. sylvestris L. collected from two different experimental areas; an unpolluted area with soil characterised by an element concentration similar to the geochemical background of Polish soils and an area polluted by tailings from copper ore flotation processes. The trees had grown in Dystric Arenosols (unpolluted area); in a mixed forest in the Greater Poland province (mid-west region of Poland) and in Spolic Technosols (polluted area) in a flotation tailings disposal area in the Lower Silesian province (southwest region of Poland). They were not only characterised by similar height, but also by the size of their upper branches. From each experimental area five specimens of P. sylvestris L. were collected and analyzed. A detailed description of this experiment and characteristics of the mineral composition of the experimental areas was published in our previous study ).

Sample collection
In accordance with the method of Hammer and Keller (2002), the root zone from the unpolluted area (soil) and the polluted area (flotation tailings) was sampled from the surroundings of the P. sylvestris L. plant roots. Shaking and gentle cleansing by hand detached the soil attached to the roots. Samples of the root zone as well as soft roots (unheated) were separately preserved in polyethylene bags. Needles were collected from both the bottom and top of P. sylvestris L. trees and immediately placed in a portable refrigerator. Secured samples were then transported to the laboratory.

Preparation of rhizosphere, roots and needle samples of P. sylvestris L. for ALMWOAs and phenolic compound analysis
Environmental subsamples collected from the two significantly different environmental surfaces were carefully cleaned of the residue of roots and any other extraneous materials, homogenised, dried at room temperature, sieved by a nylon fibre sieve (< 1 mm) and stored for subsequent analysis.
Roots were immersed in 0.01 M HCl cold solution in order to eliminate trace elements adsorbed at the root surface (Adeniji et al. 2010), washed with cold deionised water and gently dried on a filter paper to remove excess water. Samples of P. sylvestris L. roots and needles (~ 1.0 g), ground to powder in a mortar chilled using liquid nitrogen were collected in 50 mL centrifuge tubes and stored frozen (− 80 °C) until analysis. The extraction method for ALMWOAs and phenolic compounds analysis in all the studied matrices was presented in detail by Magdziak et al. (2017) and Gąsecka et al. (2017). The obtained solutions were evaporated to dryness and stored frozen (− 80 °C) until analysis.
Samples prepared from the rhizosphere, roots and needles before HPLC analysis were dissolved in 1 mL of deionised water, centrifuged and filtered through filters of 0.22 μm immediately prior to chromatographic analysis. For the determination of ALMWOAs and phenolic compounds 10 μL of liquor was injected onto the HPLC column C 18 according to the method presented by Magdziak et al. (2017) and UPLC phenolic compounds analyses as described by Gąsecka et al. (2017).

Preparation of samples of P. sylvestris L. needles for terpenes analysis
Until analyses the needles were stored in a freezer at a temperature of − 18 °C after which they were cut manually to approx. 1-2 mm in length. Samples of 0.1 g were weighed and transferred to 15 mL glass vials equipped with a silicone-Teflon septum. In order to confirm the equilibrium of compounds in the headspace over the material each sample was subjected to a 10-min preincubation at a temperature of 40 °C in a water bath.
Microextraction on CAR/PDMS (carboxen/polydimethylsiloxane) fibre was run after the needles were placed into the vial through the septum. Absorption lasted for 10 min and, similarly to preincubation, it was run at 40 °C. Thermal desorption was run after the needles were introduced to the gas chromatograph injector heated to 230 °C. Desorption time was 5 min.
The GC oven temperature was set at 50 °C for 4 min and then programmed to 240 °C for 5 min at a rate of 10 °C/min, using He as a carrier gas (60 kPa). The injector and detector temperatures were maintained at 220 °C and 250 °C, respectively.
Volatile compounds from pine needles were analyzed by solid phase microextraction (SPME). The assay principle consists in the sorption of microscopic amounts of organic compounds in a thin, cylindrical layer of the stationary phase, which covers glass or quartz fibres. In analytical practice the headspace technique (HS-SPME) is most commonly applied, in which the fibre is placed in the headspace over the tested sample for a specific period of time, after which the fibre is placed within a needle. Sorption of volatile compounds depends on their affinity to the stationary phase. Silicone phases are the most frequently used group. The next stage consists in the desorption of the analyte to the gaseous phase which takes place after the fibre is placed in the injection port of a gas chromatograph (Lord and Pawliszyn 2000).
Chromatographic analysis was run in a GC/MS TRACE 1300 gas mass chromatograph (Thermo-Scientific) equipped with a DB 5 column (30 m × 0.25 mm i.d., film thickness 0.25 μm), a mass spectrometer with a single quadrupole and an ionization voltage of 70 eV, m/z scan range 35-350 Da. Qualitative analysis was based on a comparison of retention times and indices with NIST mass spectra libraries and other corresponding data (Adams 2007). Quantities of individual VOC components were calculated as relative concentrations (peak area percentages).

Statistical analysis
Statistical analysis was done using STATISTICA 10 and consisted of ANOVA followed by the post hoc Tukey's test. For the comparison of values characterising the studied material from two independent areas, the significances (*P < 0.05, **P < 0.01, ***P < 0.001) between control and treated plants were determined using a Student's Test. Principal Component Analysis, PCA, was performed to present the relationships between independent variables (the content of separately ALMWOAs, phenolic compounds and terpenes). To determine if the average content of individual ALMWOAs and phenolic compounds in two variants: root vs needle on the bottom and root vs needle on the top are equal, a t test with Welch correction for unequal variance (t test () in R) was used. These analyzes were performed independently for plants from unpolluted and polluted areas.

Results
ALMWOAs and concentration of phenolic compounds in rhizosphere zone of P. sylvestris L.
Concentration of ALMWOAs was lower in the rhizosphere zone than in roots and needles ( Table 1). The profile and concentration of acids were strictly dependent on the condition of P. sylvestris growth.
The total amount of identified organic acids in the rhizosphere zone of both soil and flotation tailings was almost the same. However, in the case of profile, the dominant acids in rhizosphere soil extracts from the unpolluted area were citric acid (~ 41.9 μg kg −1 dry weight (DW)) followed by malic and oxalic acids (~ 16.8 and 16.7 μg kg −1 DW, respectively), whereas from the flotation tailings root zone, malic and oxalic acids (~ 35.6 and 46.2 μg kg −1 DW, respectively). The rest of the acids were found at lower concentration and fumaric acid in both root zones was below the limit of detection. Phenolic compounds in the rhizosphere of P. sylvestris L. were also below the detection limit.

ALMWOAs content in roots and needles of P. sylvestris L.
Six ALMWOAs were identified in the root and needles extracts (acetic, citric, fumaric, malic, oxalic and succinic) of P. sylvestris L. plants. As shown in Table 1, the content of organic acids in the roots of pine had a tendency to increase. Among the detected acids, acetic, citric, malic and oxalic predominated in roots (Table 1). The amount of the aforementioned acids was significantly higher in P. sylvestris L. roots from flotation tailings in comparison to the unpolluted area. It should be noted that a significant increase was observed for citric and malic acids (~ 5.2 and ~ 8.5 times more, respectively, in comparison to control plants from soil), while, particularly in the case of acetic and oxalic acids, the situation was opposite (reduction was observed more than ~ 2 and 1.6 times, respectively). Fumaric acid was present in roots from unpolluted soil, while in roots of trees from flotation tailings it was below the limit of detection. The sum of ALMWOAs content in fine root samples increased ~ 2.4 times for plants harvested from the polluted area.
The dominant ALMWOAs present in P. sylvestris L. needles from the bottom and top of the control plants were acetic and oxalic (Table 1). The content of these acids decreased in whole needles of pine trees grown in the polluted area in comparison with the control plants. The exception were needles collected from the bottom in pine growing on flotation tailings, where the content of oxalic acid increased almost ~ 1.6 times. In the P. sylvestris L. needles collected from flotation tailings, the dominant acids were citric, malic and succinic. However, their content exhibited several differences related especially to the height of needles (bottom or top) and the area conditions of growth and accounted for approximately 14-37% of the total acid content. A significant effect on the citric acid creation in the bottom and top needles was noted, where their content increased in needles sampled from P. sylvestris L. growing on flotation tailings (more than ~ 3.5 times). In the case of malic and succinic acids, their content also increased significantly, but only in needles from the bottom (~ 5.9 and ~ 4 times higher, respectively). The content of these acids in needles from the top increased slightly (~ 1.18 and ~ 1.8, respectively). The lowest content of all the analyzed acids was observed for fumaric acid. In the case of needles collected from the bottom, its content increased in comparison to the control plants, where its content significantly decreased (~ twofold) in needles collected from the top.
More, for a graphical presentation of the obtained results and relationships between ALMWOAs determined in roots and needles of plants collected from unpolluted and polluted areas, a Principal Component Analysis (PCA) was performed. In the analysis of ALMWOAs, 64.75% (45.70 + 19.05) of variability was explained, what good reflects of observed relationships. PCA analysis allowed to show the high concentration of malonic acid in the rhizosphere of P. sylvestri from the unpolluted area ( Fig. 1).

Fig. 1
Principal component analysis for low molecular weight organic acids contents in rhizosphere and studied Pinus sylvestris L. organs (U Unpolluted, P Polluted). ALMWOAs: 1-acetic; 2-citric; 3-fumaric; 4-malic; 5-malonic; 6oxalic; 7-succinic; 8-sum Moreover, the higher content of malonic, succinic and sum of ALMWOAs in needles on the bottom and needles on the top, than in the same tissues of P. silvestris collected from the unpolluted area was observed. It is worth to underline a high content of citric acid in roots of plants from polluted than the unpolluted area was also recorded.
Besides, the content of ALMWOAs in roots was compared with their content in needles (both on the bottom and the top, separately) for plants collected from the unpolluted and polluted area. As it is shown in Table 2, significant differences were especially found for plants growing at the unpolluted area for almost all detectable ALMWOAs, but at different levels of significance.
However, no significant differences were found only for acetic acid determined in the root and the needle on the bottom and also citric acid and the sum of ALMWOAs in the root and the needle on the top, for plants collected from the polluted area.
The dominant was p-coumaric acid. Gallic and vanilic acids were also abundant. The content of the acids ranged from 10.5 to 79.9 μg g −1 DW. The content of phenolic compounds found in the contaminated area was higher; the Acetic * *** n.s *** Citric * *** ** n.s Fumaric *** *** * * Malic *** *** * ** Malonic nd nd nd nd Oxalic *** *** *** ** Succinic *** *** * ** Sum *** *** *** n.s increases were significant for caffeic, chlorogenic, ferulic, p-coumaric and sinapic acids and apigenin, katechin, kaempferol, quercetin, rutin and vitexin. The profile of needles at the bottom included phenolic compounds from roots enriched by 4-HBA, protocatechuic and salicylic acids. The content of phenolic compounds exceeded 1 μg g −1 DW for nearly all components (besides kaempferol and naringenin) up to ~ 31 μg g −1 DW. Among phenolic compounds protocatechuic acid > siringic acids > chlorogenic acid > sinapic > luteonin > quercetin > rutin were dominant and exceeded 10 μg g −1 DW. However, a significant increase in the content of phenolic compounds in the flotation tailings area was noted for 4-HBA, caffeic, ferulic, p-coumaric, protocatechiuc, salicylic (only in contaminated area), sinapic, syringic and t-cinnamic acids as well as apigenin, luteolin, naringenin and quercetin. The phenolic composition of top needles was similar to those at the bottom with the exception of 4-HBA and chlorogenic acid, which did not occur in the unpolluted area. The content of all the detected compounds was above 1 μg g −1 DW, but did not exceed 45 μg g −1 DW. A significantly higher content of phenolic compounds was detected in the flotation tailings, excluding apigenin. A reduction of the content of ferulic and galic acids was noted in the polluted area. The PCA analyses were also performed for phenolic compounds to show similarities and differences in their content in the studied trees of P. sylvestris, where 83.95% (66.59 + 17.36) of variability was explained. Dominant phenolic compounds in roots of P. silvestris from polluted and unpolluted areas were p-coumaric acid and kaempferol, while in case of needles on the bottom from both areas acids: protocatechuic, chlorogenic, synapic and also luteolin, quercetin and rutin (Fig. 2). Additionally, needles from the polluted area showed higher content of the majority of phenolic compounds than those from the unpolluted area.
Significant differences between the content of majority phenolic compounds determined in roots and needles (on the bottom and the top, separately) of plants collected from the unpolluted and polluted area were observed. According to data in Table 4, significant differences were found in almost all cases besides for sinapic acid in the roots and the needles on the top of plants from the unpolluted and p-cumaric acid for root vs. needle on the bottom from polluted area.  Table 5, the characteristics of terpenes emitted from the needles of pines growing in diverse forest area types are shown.
A total of 24 compounds were identified, present in pine needles collected from the polluted and unpolluted area. However, significant differences were found in terms of amounts of individual terpenes (Fig. 4).
In pine needles from the polluted area α-pinene was the dominant compound. In needles collected from tree tops its mean content was 27.9%, while in needles sampled from the middle sections of trees it was 23.3%. In terms of the amounts emitted from these needles 3-carene ranked second. In needles collected both from the tops and middle sections of trees growing in the contaminated area the levels of this compound were very similar, amounting to 16.3 and 16.7%, respectively. Other monoterpenes were dominant in needles of pines growing in the unpolluted area. In needles from both the central and apical parts of trees 3-cerene was emitted in the greatest amounts, at 27.2 and 24.4%, respectively. In turn, α-pinene ranked second in terms of the emitted amounts among monoterpenes emitted from pine needles from the unpolluted area. The levels were 17.2% from needles in the apical parts and 16.7% from needles in the middle sections of trees. Differences between emissions of α-pinene from needles of pines growing in the polluted and unpolluted areas amounted to approx. 10% (27.9 vs. 17.2% and 16.9 vs. 23.3%). There are significant differences between the content of dominant terpenes, as confirmed by statistical analysis (Table 5). Greater amounts of β-pinene were also emitted by needles of pines from the polluted area (2.93 and 3.31%) in comparison to those from the unpolluted area (2.53 and 2.64%). The emission of bornyl acetate from pine needles from the polluted area (tree top-1.45% and middle sections-1.95%) was almost threefold greater than the emission of this compound from pine needles from the unpolluted area (tree top-0.59% and middle sections-0.67%).
In case of relationships between terpenes determined in P. sylvestris needles (on the bottom and the top) collected from studied areas, the PCA analysis 94.86% (66.75 + 28.11) of variability was explained, what showed a good reflect of observed relationships. From the unpolluted area, the higher content of 3-carene, 2-hexenal, 4(10)thujene, ς-terpinene and cubedol than in needles from the polluted area was observed (Fig. 3). An opposite situation for β-pinene, boryl acetate, germacrene D, camphene and tricyclene was recorded (Fig. 4).

Discussion
Different conditions of growing trees affect not only their growth and development but also the diversity of the amount and composition of biologically active compounds (ALMWOAs, phenolic compounds, terpenes, proteins, sugars) secreted to the rhizosphere, as well as created and accumulated in the roots or leaves or/and needles of trees (Viehweger 2014). In turn, the tree species is a decisive factor in the dendroremediation process and one that determines the essence of the plant's response to metals. Therefore, our research focused on the analysis of organic acids, phenolic compounds and terpenes as a response to the high content of metals in post-flotation wastes. Such studies are necessary, but what is most important for P. sylvestris is really new.
The presented work showed that P. sylvestris responded differently, a fact which was strictly dependent on the physicochemical parameters of the soil where P. sylvestris grew. Metal/metaloid (As, Hg or Tl) present in flotation tailings were the main metals influencing the exudation of ALM-WOAs by P. sylvestris roots; metal stress having been frequently reported to influence the exudation of ALMWOAs (Qin et al. 2007;Mucha et al. 2010;Meier et al. 2012). These elements have been previously shown to affect the exudation of oxalic, malonic and formic acids and were found as dominant in Ulmus leavis Pall. ) and oxalic, malic and succinic acids in Pteris vittata (P. vittata) (Das et al. 2017). Acids exudation has been frequently considered to be a resistance and detoxifying mechanism against metal contamination. This stimulation might therefore have been triggered to reduce metal toxicity by, for instance, complexing their ions. However, in the case of P. sylvestris, the significantly higher concentration of metals in flotation tailings inhibited the exudation of oxalic and malic acids, but highly induced citric acid secretion. The total concentrations of the analyzed acids were lower in comparison to the rhizosphere from the controlled area. In the case of oxalic acid a possible explanation of its decrease in concentration is exudation accompanied by plant uptake of the remaining ALMWOAs present in solution (Rocha et al. 2016), and in the case of malic acid being released, P. sylvestris roots will have taken up the metal-malic complexes formed during their growth. Nevertheless, an increased amount of citric acid in the rhizosphere could plays a detoxifying role at the level of the root zone, as well as increasing its amount in the root at the cellular level. Such a situation is highly probable according to literature data, but also in accordance with the results obtained from the analysis of acids in the roots. Metal-ALMWOAs complexes have been found in soil plant tissue (Walker and Welch 1987;Rellán-Álvarez et al. 2010;Rocha et al. 2016) and ALMWOAs are well-known for their ability to form complexes with metal ions, being involved in several processes such as metal tolerance, metal transport through xylem and metal sequestration in vacuoles (Clemens 2001). Therefore, in the presented study, ALMWOAs were determined not only in the root zone, but also in roots, thereby supporting the possible explanations for the content of the studied acids.
Woody plants can tolerate high metal pollution, utilising different defence mechanisms ). The phenomenon of rhizodesposition has been shown to be a plant defence mechanism against metal toxicity (Miyasaka et al. 1991;Ryan et al. 2001;Pinto et al. 2008) and ALMOWAs obtained in roots could be the main molecules for decreasing metal toxicity in P. sylvestris roots ). Moreover, the higher content of ALM-WOAs obtained in pine roots may indicate possibilities of compartmentation of metals by the studied ALMWOAs in roots and their probable role as ligands, confirming previous research Potdukhe et al. 2018). However, it should also be noted that concentration of ALMWOAs in the rhizosphere as well as in roots is strictly correlated, hence it may be simultaneously assumed that their ability to tolerate metal is also correlated. What is also important in the dendroremediation process is that pines have been shown to uptake increasing metal content in their roots as exposure to metal concentration increases , with more than two times the content of ALMWOAs in roots in comparison to the control.
These data have also been extended with respect to the content of acids in the pine needles. Regardless of the height from which the needles were removed (the bottom or the top part of the tree) the same profile of analyzed acids was found, although with a significant difference in the content. In the case of needles collected from the bottom, significantly higher concentrations of ALMWOAs were determined as compared to needles growing at the top of the studied P. sylvestris. In addition, the needles from the bottom part of the pines growing on the flotation tailings were characterised by a significantly increased content of citric, fumaric, malic, oxalic and succinic acids, while the needles from the top part showed a significant increase of citric acid only. In studies previously carried out and described by  it was also found that needles from the bottom part accumulate significantly higher concentrations of metals present in the floatation tailings, and from the top part, significantly less. Therefore, the marked increase in the total acid content within the needles from the bottom part of P. sylvestris (more than 2.5 times), could be a specific pine response to stress and may be a potential mechanism of detoxification of accumulated metals (Ma et al. 1997;Meier et al. 2012). The obtained results could also confirm previous literature reports on woody plants where a significant increase in the content of metals in the above-ground parts of studied woody plants was correlated with an increased content of metals/metalloids in the substrate . The studied ALMWOAs could maintain a balance in metal homeostasis and keep the toxicity within physiological limits in whole plants. ALMWOAs could also be used as chelates to transport, sequestrate and protect the pine from the negative effects of their growth and development (Kutrowska and Szelag 2014), which may be indicated by significant differences in the number of marked acids in the roots and their amount in the needles.
Changes in the quantity and content of phenolics in different plant tissues are reported under metal stress because they play an important role in defense and adaptive mechanisms (Martinez et al. 2016;Pradas del Real et al. 2017;Ullah et al. 2018;Drzewiecka et al. 2018). The study on P. sylvestris shows the quantitative and qualitative changes in roots and needles. Moreover significant differences in nerealy all phenolic matabolites content were detected between needles and roots in polluted and unpolluted area. The phenolic compounds detected in P. sylvestris roots occurred predominantly in small amounts not exceeding 10 μg g −1 DW. Only three components (gallic acid, p-coumaric acid and vanillic acid) archived a content > 10 μg g −1 DW. The dominant acids were p-coumaric and gallic acid. However, only the changes in p-coumaric acid in roots from different localizations were significant. Data concerning phenolic compounds in roots under metal stress are rare. Most authors have focused on total phenolic contents, pointing out the changes that occur under metal stress. A hydroponic experiment with R. communis confirmed the increase of phenolic compounds in roots and leaves under Cd and iron (Fe)-stress; moreover, a correlation between phenolics and Cd was also confirmed (Ullah et al. 2018). Additionally, a lower content of phenolics in roots than in leaves has been recorded, although no profiling of the phenolics was conducted (Ullah et al. 2018). The induction of total phenolics in roots has also been noted in Phaseolus vulgaris under Cu treatements.
The significant increase of phenolic compounds in P. sylvestris from the contaminated area confirmed their contribution in defence mechanisms. Owing to their chemical structure, phenolic compounds are able to directly scavenge reactive oxygen species ROS, chelate ions of transition metals and suppress lipid peroxidation, making the diffusion of harmful radicals and peroxidative reactions difficult (Arora et al. 2000;Michalak 2006;Younis et al. 2018). ROS are generated via Fenton-Haber-Weiss reactions; however, some metals can induce the expression of enzymes such as lipoxygenases, which indirectly generate oxidation of polyunsaturated fatty acids (Kováčik et al. 2014). According to mechanisms proposed by Jiang et al. (2017a, b) for mangrove tissues and Cd-stress, the action of polyphenols may act in the apoplast and cytoplasm and include adsorption and lignification. Enhanced lignification is a natural barrier for metal ions which immobilise metals in the cell wall (Uraguchi et al. 2006). According to Jiang et al. (2017a, b), phenolic compounds are able to affect the bioactivity and translocation of toxic metals in plant tissues because the ability to chelate metals is associated with the presence of carboxyl and hydroxyl groups in the structure (Michalak 2006) and further translocation from cytoplasm to vacuole. Pradas del Real et al. (2017) pointed out that the increase in the contents of some phenolic compounds (quercetin-3-O-(2,6-di-O-rhamnosyl-galactoside and apiin) in the roots of Silene vulgaris under chrome (Cr)-stress was related to a defense mechanism based on free radical scavenging and chelation of Cr. The confirmed changes in the content of phenolic compounds in P. sylvestris roots in the contaminated area suggested that acids are important compounds in defence mechanisms (caffeic, chlorogenic, ferulic, p-coumaric and sinapic) and some flavonoids (apigenin, katechin, kaempferol rutin, quercetin and vitaxin). In P. sylvestris roots under Zn stress phenolic changes caused the opposite reaction, related to enhanced lignification (Ivanov et al. 2012). The results of our study pointed to a significant increase of phenolic compounds, suggesting they have a prevailing role in the defence mechanism connected with antioxidant mechanisms rather than a structural role associated with lignification, although the precursors of lignin are present in roots (caffeic, p-coumaric, chlorogenic, ferulic and sinapic acids). A very surprising fact is that gallic acid is one of the most dominant compounds in roots, but changes in its content are not significant, even though gallic acid is known to be a strong antioxidant with the ability to directly scavenge ROS.
The study of López-Orenes et al. (2018) on Pinus halepensis (P. halepensis) growing on mine tailings suggests that phenolics could participate in the defense mechanism as chelating ligands playing a role in the detoxification and accumulation of metals and metalloids.
The phenolic composition of needles depends on their position in the crown, contrary to observations of Norway spruce (Kopačková et al. 2015). This could be a result of the age of the needle as demonstrated by (Kopačková et al. 2015). Little is known about the phenolic composition of needles. Roitto et al. (2005) found flavonoids in needles of P. sylvestris and the content alternated in their responses to metal content (Ni and Cu) in soil. The results of our study confirmed a wide range of phenolic acids and flavonoids in needles on the bottom and the top with significant changes in the content of most of these components, especially in the top needles. The significant enhancement of phenolic acids and flavonoids pointed to their significant role in defense mechanisms, most likely associated with their ability to scavenge ROS (e.g. gallic acid, caffeic acid, p-coumaric acid (Benbettaieb et al. 2018) and functioning as ligands, protocatechiuc acid (Kováčik et al. 2010). In needles salicylic acid is present and the increase in its content in the contaminated area was observed as a result of its function in response mechanisms to environmental stressors, including metal stress (Arasimowicz-Jelonek et al. 2014;Drzewiecka et al. 2018). The reduction in ferulic and galic acids was probably connected with lignification, as also confirmed by (Ivanov et al. 2012).
Differing growth conditions for trees affect not only their growth and development but also result in varied amounts and composition of essential oils in needles and other plant parts. Both nutrients and pollutants found in the substrate influence the quantities and quality of the byproducts produced by plants (essential oils, tannins, resins, waxes) (Sembratowicz et al. 2008;Kupcinskiene et al. 2008;Sandre et al. 2014). Werner et al. (2004) recorded terpenoid contents in the roots of P. sylvestris seedlings growing in a substrate contaminated with heavy metals and stated that in the case of polluted substrate roots showed significantly lower terpenoid levels. Kupcinskiene et al. (2008) reported that the most heavily polluted environments trigger an increase in the concentrations of sabinene and β-pinene in needles of pines growing in the vicinity of industrial facilities. Along transects an increase was observed in the amount of some diterpenes and a decrease in the components of the shorter chain essential oils. In this study, lower amounts of certain diterpenes were recorded, i.e. ς-muurolene, cadina-1(10),4-diene and cubedol in needles of pines growing in the heavily contaminated environment. Judzentiene et al. (2006) reported the dominance of two basic monoterpenes (α-pinene and 3-carene) in oils from needles of pines growing in an environment with elevated ammonia levels in the air. However, it was stated in that study that the contents of a terpene, 3-carene (22.7-33.7%), exceeded those of α-pinene (19.8-35.0%). In contrast, needles from the polluted area showed an opposite relationship as the amounts of α-pinene (23.23-26.57%) exceeded those of 3-carene (16.57-16.89%). Similarly, to the authors of this study, Dziri and Hosni (2012) also observed an increase in the contents of β-pinene (29.5% and 22.0%) in needles of P. halepensis from a polluted area in comparison to control samples (11.6%). Supuka et al. (1997) recorded lesser amounts of terpenes in essential oils from needles collected in the area of Nitra than in oils from needles sampled in the arboretum. Moreover, like the findings of this study, analyses of the material collected in Nitra revealed an increase in β-pinene contents in needles from an urban area in comparison to needles collected from an unpolluted area. Urban conditions caused a changed proportion of α-pinene in relation to the following terpenes: α-phelandrene, 1.8-cineole, terpineol, citral and carvone (Supuka and Berta 1998). These changes discovered in the contents of terpenes in needles of Pinus strobus L. have shown that the group of secondary metabolites is an important biochemical marker of the environmental impact on woody plants. An increase in α-pinene contents in pine needles from an area polluted by dusts emitted by a cement plant was reported by Hosni et al. (2014). In the control samples those authors recorded 48.1% of this compound, while in samples from the polluted area it was 62.2% and 63.8%, respectively.
All the cited literature data are concerned with chemical changes in essential oils of plants growing in contaminated environments. Thus, the results of the performed experiment are consistent with the data published by other authors. According to Oviasogie et al. (2009) production of essential oil is an indicator of plant adaptation to habitat conditions and helps plants to easily adjust to environmental stress conditions, among others: drought, intense radiation, high temperature, high heavy metal contents. Natural selection favours the survival of plants with a higher concentration of essential oils which enhance their adaptive value (Stevovič et al. 2011). This plant defense mechanism caused an increase in the content of essential oils in the examined needles of P. sylvestris L. growing in a polluted area. However, the role of selected monoterpenes, especially α-pinene, in this mechanism requires further research.
Another interesting aspect related to monoterpens is their role in biotic interactions and soil processes. Ludley et al. (2009) analyzed the monoterpene content and distribution in the litters and roots of three species of conifer: Picea abies, Picea sitchensis and Pinus sylvestris. They showed an environmentally relevant content of α-or ß-pinene sufficient to increase mycorrhizal colonization of Picea abies root tips and to decrease the respiration rate of two species of saprotrophic fungi, in relatively natural substrata. Smolander et al. (2006) investigated the impact of monoterpens on soil microbes. According to authors there were differences in the response of both microbial biomass C and N to different monoterpenes, α-and ß-pinene being the most inhibitory whereas myrcene inhibited only microbial biomass N in the presence of arginine. Generally, they stated that monoterpens have negative effects on soil N transformations but may serve as a carbon and energy source for some soil microbes.

Conclusion
The obtained results reveal how the growth environment has a significant influence on the content of organic compounds in the rhizosphere, roots and needles. In the case of the concentration of ALMWOAs secreted into soil and into highly toxic flotation tailings, the root zone was similar. However, in roots and needles their value was significantly higher in tissue from the polluted area. Such an observation indicates that ALMWOAs may be one of the most important factors in P. sylvestris species, contributing to metal/ metalloid detoxification, involved in several processes such as metal tolerance, metal transport through xylem and metal sequestration in vacuoles. Moreover, the changes in the phenolic profile of roots and needles point to a major role in the defence mechanism associated with antioxidant mechanisms based on the scavenging ROS, functioning as ligands rather than having a structural role connected with lignification, although the precursors of lignin are present in roots (caffeic, p-coumaric, chlorogenic, ferulic and sinapic acids. The analysis of the results of volatile terpenes in pine needles growing in substrate contaminated with sludge and uncontaminated soil indicates a different dominance of the main monoterpenes. The contaminated soil had a significant effect on the reduction of 3-carene and the increase of α-pinene in pine needles.
Author contribution statement PG, MM and AM chose experimental areas and designed the study; WJC, BM, TK, PN and MM selected appropriate research material for testing; ZM, MG, BW, MZ analysed samples and took part in writing the manuscript, MG, AMP performed the statistical analysis; PN methodology. All authors contributed to writing the final version of the manuscript.