Uncovering the multi-level response of Glycine max L. to the application of allelopathic biostimulant from Levisticum officinale Koch

The interest expressed by the agriculture in the category of innovative biostimulants is due to the intensive search for natural preparations. Our study is the first ever to report a complex approach to the use of allelopathic extracts from Levisticum officinale Koch. roots in soybean cultivation, includes analyses of morphological observations, and analyses of biochemical indicators. Hot method of aqueous extraction was applied. The extracts were administered via foliar application and soil treatment. Lovage extracts had high contents of polyphenolic compounds and rich micro- and macroelemental composition. The infusions did not contain gibberellic acid and indole-3-acetic acid but the abscisic acid and saccharose, glucose, and fructose were found. The extracts modified soybean plant physiology, as manifested by changes in biometric traits. Plants responded positively by increased yield. Seeds from the treated plants had higher contents of micro- and macroelements, as well as total concentrations of lipids (with a slight decrease in protein content). In addition, they featured changes in their amino acid profile and fatty acid composition. The application of allelopathic biostimulant caused increased concentrations of isoflavones and saponins. The natural biostimulants from Levisticum officinale may become a valuable tool in the sustainable agriculture.

Because the allelochemicals cause no residual nor toxic effects, they are the perfect base components of innovative biostimulants 7 . Thus, the major aim of the allelopathic research is firstly to observe the effects of such preparations and, afterward, to exploit them in agricultural production. Only this approach will allow reducing or eliminating agri-chemicals as well as putting effective methods for the sustainable development of agricultural production and ecological systems into practice 4,[8][9][10] . Thus, research efforts in this area should focus on the search for and study of potentially allelopathic plant species producing more allelochemicals. In addition, initiatives should be undertaken to obtain such forms of preparations and methods of their application that would promote plant growth and crop yield, taking into account mechanisms of their action and their biochemical composition. This may be the right direction in the long run, and following this direction will allow achieving sustainable development of agriculture, environmental and food safety, as well as the important economic stability for farmers [11][12][13][14][15][16] .
Our earlier results indicated that the extracts produced from Artemisia absinthium can be successfully used as effective allelopathic biostimulant in soybean cultivation. The application of such bioproduct increased soybean yield, while their positive effect can be attributed to their phenolic compounds as well as micro-and macroelements. However further research is needed to thoroughly characterize the chemical composition and the mechanisms of action of the allelopathic biostimulants. One of the solutions to these challenges may be seen also in prototyping medicinal plants, which is the first step to create biostimulants of new generation. This approach will not only bring benefits for the natural environment, but will also represent some kind of a designed agronomic tool for ensuring optimal yield and economic profits 17 .
Levisticum officinale Koch., growing wild in most parts of the world, including Asia and Europe, can be indicated as a plant exhibiting the allelopathic potential 18 . It belongs to the family Apiaceae, which includes many culinary and medicinal plants that contain multiple compounds exerting broad and varied biological effects 19 . For this reason, lovage can be a candidate for producing an allelochemical biostimulant. However, such a biostimulant cannot be tested only under laboratory conditions. Its field trials should be performed with crops especially sensitive to abiotic stresses, like e.g. soybean and bean. Considering the wide applicability of soybean seeds in the food industry, the greatest challenges seem to be yield quality control and modelling. Hence, the use of an allelochemical preparation in the form of extracts from Levisticum officinale can be perceived as an agronomic practice useful in modelling soybean yield by the plant's secondary metabolism. By this means, it represents a novel, unconventional, and environmentally-friendly approach coupling allelopathy principles and biostimulant use that not only increases crop yield but also enhances phytochemical biosynthesis in plant response to the biotic stress. Obviously, plant response to an allelopathic biostimulant will depend on many factors, including plant physiological condition as well as the dose and application method of the preparation 20 . However, the innovative character of our research lies not only in the development of a novel biostimulating formulation but also in its dualistic nature. This formulation will also have features of an allelopathic preparation, i.e., it will be a 'multi-component balanced system of biologically-active substances of metabolic origin based on plant raw materials exhibiting a broad spectrum of bio-activities' 21 .
Our study includes analyses of the allelopathic potential of extracts from Levisticum officinale, meaningful morphological observations, and analyses of basic metabolic activities and biochemical indicators of soybean plants' defense systems. We have advanced a hypothesis that the allelopathic biostimulant produced from roots of Levisticum officinale would promote the growth and seed yield of soybean (Glycine max L.) due to the mechanisms based on the improved key physiological and biochemical process. This assumption was based on the fact that the biotic allelopathic stress is multi-dimensional and can induce various responses at different levels, including biochemical, physiological, morphological, and-ultimately-ecological one in plants treated with allelochemicals.

Results and discussion
Chemical composition of Levisticum officinale water extracts. Phenolic compounds identified in Levisticum officinale water extracts. The quantitative identification of infusions from Levisticum officinale roots revealed fourteen polyphenolic compounds (having various retention times) ( Table 1). The major phenolics identified included quercetin 3-O-rutinoside (118.54 μg/mL) and equisetumpyrone (58.87 μg/mL of chlorogenic acid, ascribed to its deprotonated pseudomolecular ion [M-H] at m/z 423 and its fragments at m/z 173 and 161, respectively, at the absorption maximum of 255-352 nm). The water extracts were also rich in chlorogenic acid (39.28 μg/mL) and neochlorogenic acid (16.17 μg/mL).
The extracts contained also four quercetin derivatives, identified based on deprotonated pseudomolecular ions [M-H] with m/z 301 and showing UV absorption maxima at 255-352 nm. The UPLC-PDA-MS/MS analysis of Levisticum officinale root extract revealed also the presence of two quercetin derivates (quercetin 3-O-glucuronide-pentoside and quercetin 3-O-glucosyl-rhamnosyl-(p-coumaroyl)hexoside)), that were ascribed to the deprotonated pseudomolecular ion [M-H] at m/z 433 and 771, respectively, and to their fragments at m/z 301 and 609. Their UV absorption maximum fitted in the following ranges 259-350 nm and 255-314 nm, respectively. The listed in Table 1 compounds with identical spectral characteristic, measured as a mass to charge ratio (m/z) at 609 indicates the presence of constitutional isomers. As a result of fragmentation patterns, these compounds break down and release a radical ion at m/z 301. This ion indicates the presence in the aglycone part the quercetin-flavonol, very common in nature polyphenolic compound.
The study results demonstrate high concentrations of phenolic compounds in the extracts examined, reaching 259.68 μg/mL on average, that could support natural strategies of soybean plants for coping with adverse stress factors. Besides, such a composition of the extract (i.e., a high concentration of allelochemicals) can increase the availability of nutrients for plants, because-by being the so-called chelating agents-phenolic compounds indirectly contribute to the increased solubility of nutrients in the soil 22 . Even though the higher concentrations of secondary metabolites (also in the form of phenolic compounds) may increase plant resistance, the production www.nature.com/scientificreports/ of secondary metabolites itself assumes some costs borne by the plants [23][24][25][26][27] . The precise identification of the costs of allelopathy borne by plants remains controversial and disputable 26,28,29 . Scientists have been inconclusive in this regard mainly due to a lack of experimental research addressing this issue. Results of sparse studies depend on the bio-system examined and cannot be extrapolated 30,31 . Nevertheless, phenolic compounds can be found in all plant organs and, therefore, protect it against sun radiation and attacks of pathogens and herbivores, which justifies the use of polyphenols-rich extracts in allelopathic preparations. Their presence in extracts from Levisticum officinale induces some type of biotic stress in soybean plants, thereby contributing to their protection against abiotic stress conditions, including particularly oxidative stress 32 . It needs to be emphasized that the production and removal of reactive oxygen species in a cell play an essential role in allelopathic effects. After the administration of allelochemicals, the acceptor plants can rapidly synthesize reactive oxygen species at the contact site with the allelopathic preparation 33,34 , and modify activities of such antioxidative enzymes as superoxide dismutase and peroxidase 35,36 , or ascorbic acid peroxidase 37 . Therefore, this composition of the Levisticum officinale extract contributes directly to counteracting oxidative stress. As explained by Shearer et al. 38 , the allelopathic interactions can lead to changes in signal transduction and to the imbalance between the production of reactive oxygen species and antioxidative capabilities of acceptor plants. According to Van Oosten et al. 39 , soybean plants have stress-sensitive phenotypes. The phenotype of these plants can be changed due to the exogenous use of various biostimulants, especially those showing antioxidant and free radical scavenging properties, as well as the ability to chelate metals 40 . Thus, foliar application of Levisticum officinale extracts, rich in polyphenolic compounds and micro-and macroelements, could cause their rapid uptake by the leaves and increased accumulation in meristematic tissues. According to Park et al. 41 , the uptake of bioactive compounds in the most sensitive plant tissues has a great importance in cases of thermal and oxidative stress, to which crops are exposed. Exogenous use of allelopathic extracts with compatible solutes such as polyphenols or micro-and macroelements can affect the flow of ions across plant membranes 42 . However, it should be emphasized that the cause-and-effect relationship between the application of biostimulants, phenotypic evaluation of the treated plants and the accumulation of solutes and protection against stress is still not fully understood 43 . According to Nephali et al. 44 even more careful research doesn't always show a clear answer. It turns out that changes on the transcriptome or proteome levels do not always correspond to changes in biochemical phenotypes. Additionally, they may not fully reflect the exact biochemical condition of plants in response to stress and the application of biostimulants. It should be mentioned that the suggested resistance of plants to oxidative stress, resulting from the composition of the applied allelopathic biostimulants, could also be the result of gene activation, e.g. the ALDH7 gene, associated with oxidative stress (protective function for plants) 45 . Research by Wang et al. 46 showed the existence of 51 co-functional networks of GmALDH genes along with 1146 co-functional connections that are involved in the processes of lipid metabolism, photosynthesis, proline catabolism and catabolism of small molecules. The application of allelopathic biostimulants led to changes in the concentration of these components in soybeans, which may probably be the result of shaping the expression of these genes. In addition, the activation of these genes, which support soybean against stress conditions, may affect the resistance of soybean plants to oxidative stress, but without changes in plant phenotypes 45 . Although, in own research it was shown that the number of pods and plant height changed, however, further studies are required to indicate the mechanism responsible for the observed phenotypic changes in soybean plants.  Table 2 presents the multielemental composition of water extracts prepared from powdered roots of Levisticum officinale. Results obtained provide an important clue, indicating that the extracts were highly diversified in their both macro-and microelemental composition. The infusions were especially rich in macroelements. The allelopathic extract had high concentrations of calcium, potassium, magnesium, and sulfur. The extracts obtained from a potentially allelopathic plant contained iron, copper, manganese, zinc, and molybdenum. Interesting outcomes of analyses were achieved in the case of molybdenum and zinc. Their levels are especially controlled in soybean biostimulation not only because of their yield-promoting effect but also because of their influence on the appropriate course process of nodulation and, thus, on nitrogen assimilation. While molybdenum is an essential element of an enzymatic complex-nitrogenase, the yield-promoting effect of zinc is primarily due to the activities of plant hormones (auxins in particular), whose deficiency results in growth inhibition 47 . Our previous research on the use of Artemisia absinthium extracts also proved that this plant has allelopathic potential. It was shown that the allelopathic biostimulant produced from this plant influenced the development and yielding of soybean. However, compared to Levisticum officinale extracts, it was less rich in macro-and microelements. Particularly significant differences were observed in the levels of calcium, magnesium, phosphorus and sulfur. Nevertheless, both preparations showed a biostimulating effect 17 . Summing up this stage of the present study, it can be concluded that the extracts from Levisticum officinale contained a pool of elements essential for the successful growth and development of plants. Plant response to the application of the allelopathic extracts was also associated with the effects of micro-and macroelements on certain vital processes, considering the proven capability of selected elements for stimulating the mechanisms of plant resistance to biotic and abiotic stresses or aiding the uptake of other nutrients. In addition, such a multielemental composition of the extracts from Levisticum officinale was proved to even compensate for or alleviate the toxic effects of other elements or induce adaptive responses of plants to adverse abiotic stresses 48 .
The research showed that the Levisticum officinale extract was characterized by a low acid pH. The water extracts produced contained three carbohydrates, i.e., saccharose, glucose, and fructose ( Table 2). Infusions from Levisticum officinale roots were especially rich in saccharose. Nearly three-fold lower concentration was noted for fructose and 17-fold lower for glucose. This carbohydrate composition of the allelopathic extracts affects their bioactivity. According to the latest theory presented by Kumar et al. 49 , the positive and initiating effects of extracts from Levisticum officinale can be due to the induction of the so-called 'sweet resistance' by sugars, namely plant resistance controlled by saccharide compounds that enhances the defense system 49 . By this means, this theory also assumes that the alternative carbohydrate-rich extracts influence the primary antioxidative system of plants by counteracting stress through ensuring cellular equilibrium between redox reactions. This process triggers beneficial changes in the defense system of plants by initiating and inducing their tolerance of stress 50 . The 'sweet resistance' is one of the newest concepts in plant stress physiology, while only a few carbohydrates have Table 2. Multielemental composition of water extracts from Levisticum officinale obtained by hot extraction (infusion) (average ± SD). < LLD-below limit of detection. www.nature.com/scientificreports/ so far been classified as growth and development promoters. The results of our previous research on prototyping extracts from Artemisia absinthium L. for their biostimulating properties showed that these preparations also contained sugars. However, the concentration of individual of them was lower compared to the extracts from Levisticum officinale. However, both allelopathic extracts stimulated the yield potential of soybean in three-year field cultivation 17 . The composition of the tested extracts from Levisticum officinale shows they are a complicated system of active compounds, with carbohydrate content being one of their strong points. This is due to the multi-faceted role of carbohydrates in plant growth and development. Following the presented theories, it may be concluded that an appropriate content of sugars in allelopathic extracts would determine their biostimulating activity and plant response at the physiological and biochemical levels.
Phytohormones in water extracts from Levisticum officinale. The extract from Levisticum officinale contains a wide range of phytochemicals. Notwithstanding to the mineral nutrient elements, polyphenols and sugars, growth enhancing hormones have not been found in tested extract. Analysis showed that extract did not contain gibberellic acid and indole-3-acetic acid (Table 3). On the other hand, the Levisticum officinale water extract contained the abscisic acid in concentration 0.58 ± 0.05 µg/ml. Plant hormones, like abscisic acid, are signaling phytohormones with different regulatory roles in plant metabolism and adaptation to abiotic stresses 51 . ABA stimulates stomatal closure and minimizes water loss by transpiration 52 . Therefore, exogenous ABA application can have a great interest in plant growth biostimulation 53 . The results of our research indicates that allelopathic extract from Levisticum officinale elicit the response of soybean plants not only due to the presence of phenolic compounds and multielemental composition but also to ABA concentration. According to Mahdavi et al. 54 and Zhang et al. 55 the application of abscisic acid both under optimal and sub-optimal growth conditions can increase water potential and improve chlorophyll concentration in soy bean. Nowadays the role of ABA on the reduction of stress in protein-rich crops has received considerable attention.

Response of soybean plants to biostimulating allelopathic extracts. Soybean yield and biomet-
ric traits. The spraying of soybean plants with an extract from Levisticum officinale promoted their growth, leading to plant height increase by 13% compared to the control combination ( Table 4). The application of the allelopathic extracts into the soybean rhizosphere also improved plant growth, but only by 7% compared to the control combination. The study results showed also that the response of soybean plants to the application of extracts did not affect the first pod location. The foliar application of the allelopathic extract from Levisticum officinale increased pod number by 23%, whereas its application into the soil-by 17%, compared to the control combination.
The response of plants to the biotic stress induced by the allelopathic biostimulant was negative, as evidenced by the 1000 seed weight which decreased significantly in the combinations involving plant treatment with the extracts tested. However, this negative plants' response had no unbeneficial effect on seed yield, which was higher upon the foliar application of the extract from Levisticum officinale (by 23% compared to the control combination) ( Table 4). The results of our study can be explained by the composition of the allelopathic biostimulants which contained significant amounts of not only polyphenolic compounds but also micro-and macroelements. As reported by Salwa et al. 56 , the microelements of the preparations could have a key influence on crop growth and yield because their availability for plants at particular stages of their phenological development determines the final yield. Table 3. Phytohormones in water extracts from Levisticum officinale obtained by boiling (infusion) (average ± SD). < LLD-below limit of detection. www.nature.com/scientificreports/ However, this response resulted in their increased yield, probably due to the induced activity of their immune system 57,58 , which-in response to the use of allelochemicals-activated many detoxifying enzymes and transporters in order to facilitate toxin inactivation and elimination from plants and to adjust respective metabolic processes 3 .
Multielemental composition of seeds from soybean plants treated with extracts from Levisticum officinale. Table 5 presents the multielemental composition of seeds of soybean treated with infusions from Levisticum officinale roots. The results of analyses demonstrate that the seed yield was highly diversified in terms of macro-and microelemental composition. Inducing the biotic stress in plants through the application of the allelochemical extracts increased concentrations of macroelements in seeds compared to the control treatment. A small decrease was only observed for sulfur concentration.
Higher, on average, contents of microelements were determined in the seeds from plants treated with the allelopathic extracts from Levisticum officinale in the form of infusions.
The study results presented above indicate that, at this stage of evaluation, the activity of extracts from lovage roots allows classifying them to the group of allelopathic biostimulants. Undoubtedly, they increased the effectiveness of nutrient uptake and assimilation, and improved quality traits of soybean yield. However, such properties of the allelopathic extracts could be achieved only when they were administered by spraying.
Protein contents, its amino acid and fatty acids composition in seeds from soybean plants treated with extracts from Levisticum officinale. Protein contents and its amino acid composition are presented in Table 6. Globulins were the dominant fraction of soybean seeds. The highest content of this fraction was determined for the control seeds, while the foliar application of infusions caused a slight decreased in this fraction by c.a. 4%. Contrary the seed of plants sprayed with infusion were characterized by higher amount of albumins. In fact, changes in globulins fraction were reflected in slightly lower total content of proteins. It should be noted that content and quality of proteins is affected by both genetic (e.g. variety) and breeding factors (e.g. fertilization, stress). An increase in albumin fraction may be due to over-synthesis of anti-antioxidant enzymes and pathogen-related proteins usually observed after application of natural inductors 59 . On the other hand, a slight decrease in globulins, being a main storage fraction in soy-bean, may arise from utilization of energy in development of systemic acquired resistance rather than in production of storage materials.
The soybean seeds differed also in their amino acid profile. In spite of lower contents of total proteins the seeds obtained after treatments were characterized by increased amounts of aminoacids. It may be explained by the fact that the control seeds contain higher amount of globulins that in soy-beans are built with up to 40% of sugar moieties 60 . Significant changes were noted in contents of all amino acids, presumably due to the effect of the plant extracts on this metabolic pathway in plants. In turn, proline content was observed to increase. This Table 5. Effect of soybean treatment with extracts from Levisticum officinale on the content of elements in its seeds (means of 2018-2020 ± SD). < LLD-below limit of detection; Means in the rows, concerning the selected traits, followed by different letters are significantly different at p < 0.05. www.nature.com/scientificreports/ amino acid can be found in various plants, including especially those exposed to stress conditions. It has been reported to serve many physiological functions, like osmoregulation or energy and nitrogen absorption. It has also been shown to serve as an ageing signaler and the so-called stress sensor 61 . The application of Levisticum officinale extracts also induced arginine accumulation. It needs to be emphasized that poor transamination of arginine can lead to amino acid degradation and promote ammonia production 62-64 . Table 6. Effect of soybean treatment with extracts from Levisticum officinale on the protein content, its amino acid profile and fatty acid composition (means of 2018-2020 ± SD). Means in the rows, concerning the selected traits, followed by different letters are significantly different at p < 0.05. www.nature.com/scientificreports/ The branched-chain amino acids (isoleucine, valine, and leucine) occurred in higher concentrations in the seeds from crops in which the allelopathic extracts were administered through foliar application. Furthermore, the study results showed varied responses of plants to the treatment with allelopathic extracts manifested in the levels of alanine, aspartic acid, and glutamic acids in their seeds. Higher contents of these amino acids were determined in the seeds from plants treated with the extracts than in the control seeds. In response to the stress factor, these amino acids act as active precursors of other macromolecules because glutamic acid is involved in nitrogen metabolism and TCA cycle, whereas aspartic acid and alanine in the transamination. The differences in the amino acid composition of soybean seeds could be due to the coupled effect of various factors; however contents of amino acids in soybean seeds were not necessarily correlated with the total protein content of the seeds.
The results of the study concerning the fatty acid composition of soybean seeds point to the ambiguous response of the plants to the biotic stress induced by the application of extracts from Levisticum officinale ( Table 6). The relative percentage of fatty acids differed among soybean seeds depending on the combination tested. Concentrations of some of the fatty acids increased, while those of the others decreased in soybean seed.
However, it needs to be emphasized that the concentration of unsaturated fatty acids in soybean seeds is one of key indicators of their quality, considering their final application and market price. Oleic, linolenic, linoleic, and elaidic acids are regarded as extremely important and valuable unsaturated fatty acids of soybean seeds. The use of extracts from Levisticum officinale in the form of spraying resulted in slightly increased concentrations of linoleic and linoleaidic acids, compared to the control seeds. Somehow different observations were made for the levels of oleic and elaidic acids. Their concentrations increased in the seeds from one combination, namely that where the allelopathic extracts were applied by spraying. In the case of α-linolenic acid, soybean treatment with infusions from Levisticum officinale decreased its content in the seeds. Soybean plants responded to the treatment with extracts via the enhanced accumulation of unsaturated palmitoleic acid in seeds, compared to the control samples. Interesting dependencies were observed for cis-13,16-docosadienoic acid, which was not detected in the control seeds, but appeared in the seeds from plants treated with Levisticum officinale extracts.
The study results show also the effect of various methods of extract application on the levels of saturated fatty acids in soybean seeds. The content of palmitic acid was observed to decrease, whereas contents of stearic and arachidic acids-to increase upon soybean plant treatment with Levisticum officinale extracts. In turn, lignoceric acid present in control seeds was not detected in the seeds from plants treated with the allelopathic preparations.
In the case of soybean seeds, important are not only the contents of individual fatty acids but also the ratio of unsaturated to saturated fatty acids. The plant response to the application of infusions from Levisticum officinale was manifested by a higher content of saturated fatty acids. In turn, plant watering with the extracts stimulated an increase in the content of monounsaturated fatty acids in the seeds, compared to the control samples. While, contents of polyunsaturated fatty acids increased as a result of the foliar application of the extracts in soybean cultivation.
Soybean plants' response to the biotic stress was also manifested in changes in the unsaturated to saturated fatty acids ration in seeds. They reached 4.98 and 4.94, compared to 4.99 noted in the control samples. The lower ratios are indicative of the lower percentages of oleic, linoleic, and linoleic, and higher percentages of palmitic and stearic acids synthesized in soybean seeds. An increased percentage of unsaturated fatty acids compared to the saturated fatty acids shows that plant treatment with the extracts had a small negative impact on the nutritional value of soybean seeds.
The increased contents of unsaturated fatty acids in soybean seeds can be ascribed to the effect of extracts from Levisticum officinale on the enhancement of biosynthesis and desaturation enzymes of the fatty acid metabolic pathway. This beneficial effect could be due to the contents of macro-and microelements because certain microelements, like iron, are essential for the synthesis of ferredoxin, which plays an important role in NADPH as an electron donor to stearyl desaturase in the higher plants 65 .
The application of extracts modified also soybean response expressed in the levels of omega 3, 6, and 9 acids. In the case of omega 3 acids, plant treatment with infusions contributed to a slight decrease in their content in the seeds, compared to the control treatment. The level of omega 6 acids increased in the seeds from plants treated with the extracts by spraying. In turn, the content of omega 9 acids tended to increase, compared to the control, only upon extracts application into the soil.
Isoflavonoid and saponins compounds identified in seeds from soybean plants treated with extracts from Levisticum officinale. Plant response to the biotic stress, manifested in the accumulation of daidzin, genistin, and glycitin isoflavones, and their corresponding conjugated forms is presented in Table 7. The results of our study demonstrate that the application of allelopathic biostimulants in soybean cultivation caused changes in the contents of individual isoflavones. The differences observed were due to the various methods of biostimulant application. The total concentration of these bioactive compounds in seeds increased after plant treatment with infusions in the form of spraying and watering (an increase by over 55% and 36%, respectively, compared to the control samples). The assessment of soybean plant response to the use of allelopathic biostimulants revealed that they responded to this type of biotic stress with an increased concentration of individual isoflavones in seeds.
The analyses demonstrated malonylglucosides (malonylgenistin and malonyldaidzin), followed by genistin to be the major forms of isoflavones in soybean seeds. In turn, the lowest content was determined for acetyl genistin.
According to Kim et al. 66 , the acetyl forms are only products of degradation formed from malonate. In these authors' opinion, it is also likely that all isoflavonoids accumulated in the seeds were earlier in the form of malonyl. The results of our study demonstrate that the application of extracts from Levisticum officinale influenced the synthesis of isoflavonoids, probably because soybean seeds were the main site of their synthesis. However, the synthesis process itself was initiated by active compounds of the allelopathic biostimulant. www.nature.com/scientificreports/ indicates that some pool of isoflavonoids in seeds can result from their transport from other morphological parts of the plants 66,67 . Nevertheless, as emphasized by Kim et al. 66 , the accumulation of flavone compounds in soybean seeds is also largely dependent on the cultivar, meteorological conditions, and plant exposure to stress factors at the stage of seed growth and development 66,68,69 . Flavonoids represent a large class of plant secondary metabolites, exhibiting antioxidative effects in the higher plants, which are exposed to or challenged with various stress-inducing factors [70][71][72] . This hypothesis may explain Table 7. Individual isoflavonoid and saponins compounds identified by UPLC-PDA-MS/MS in soybean seeds (means of 2018-2020 ± SD). Means in the rows, concerning the selected traits, followed by different letters are significantly different at p < 0.05. www.nature.com/scientificreports/ changes in their concentration in soybean seeds caused by the extracts from Levisticum officinale. Biotic stress induction in plants by the allelochemicals could affect the activities of enzymes involved in the oxidation process of polyphenols 73 . At this stage of the study, we may also conclude that the preparation tested revealed the traits of a biostimulant because it influenced the metabolic pathways in soybean plants. Plant response to this type of stress was manifested in the beneficial change in the concentrations of biologically-active compounds in the seeds. This positive response could be due to the fact that the plants exposed to severe stress conditions accumulate flavonoids substituted by a dihydroxy-B-ring that are effective scavengers of reactive oxygen species 74 .
The analyzed soybean seeds contained many types of saponins, classified to group A, B, and E ( Table 7). The application of the allelopathic water extracts in the form of plant spraying and watering contributed to an increased content of saponins in the seeds. In soybean seeds tested, we detected the presence of soyasaponin I, II, III, and IV. The highest concentration was reported for soyasaponin I. Soybean seeds contained also saponin BeA, which is a saccharide, 2,3-dihydro-2,5-dihydroxy-6-methyl-4H-pyran-4-one (DDMP), attached through an acetal linkage to the C-22 hydroxyl of algycones of soyasaponin I. Results of previous investigations 75,76 also show that soyasaponins I, II, III, and IV are saponin forms conjugated with DDMP.
The analysis of soybean seed composition revealed additionally 4 types of group A saponins, called A1, A2, A3, and A4 (having retention times of 3.53, 3.66, 3.70, and 3.94 min, respectively). The highest concentration was determined for soyasaponin A3. This group of saponins was also affected by plant treatment with the allelopathic extracts, which resulted in their increased levels in the seeds.
Among the acetyl-soyasaponins A, the largest group in soybean seeds was represented by acetyl-soyasaponins A1. Results of our analyses demonstrate also the presence of B group saponins (7 groups) containing DDMP, called soyasaponins Bd, Be, BdA, BeA as well as βa, γg, and γa. The most abundant group in soybean seeds was represented by soyasaponins Bd, whereas the smallest one by soyasaponins βa. Results of our study demonstrate that the use of extracts from Levisticum officinale, in the form of both foliar application and soil treatment, increased the content of this group of saponins in soybean seeds.
Our study also showed that the application of extracts from lovage affected the isoflavone to soyasaponin ratio. Plant response to this type of biotic stress was manifested in its value increase from 0.046 in the control samples to 0.059 in the seeds from plants sprayed with the extract; with the maximal reported value at 0.061. However, according to Kim et al. 66 , changes in the contents of flavonoids and saponins are not correlated. Nevertheless, the most important finding from our study is that the application of the completely natural, allelopathic biostimulant offers an agronomic tool enabling to increase contents of compounds, whose bio-activity determines the healthpromoting values of soybean products 66 .
It is noteworthy that the increased levels of saponins determined in soybean seeds are due to complicated reactions proceeding in plants also under the influence of stress factors, mainly as a consequence of modified gene expression. As reported by Shimoyamada et al. 77 , the Sg-3 gene encodes glucosyltransferase, which afterward catalyzes the glucosylation of the galactosyl group of saponin A. In turn, Tsukamoto et al. 78 have proved that the Sg-4 gene is an arabinosyltransferase which catalyzes the arabinosylation of glucuronic acid residue attached to C-3 of soyapapolyols. However, the above authors have emphasized that these glycosyltransferases have not been identified so far. A study conducted by Kurosawa et al. 79 supported these observations, as these authors have demonstrated that soybean glucuronosyltransferase is a specific enzyme of the UDP-glucuronic acid as a donor and of soyapogenols as an acceptor. In addition, Kurosawa et al. 79 confirmed that this enzyme played an important role in the biosynthesis of saccharide chains in soyasaponins. In turn, Shibuya et al. 80 characterized two soybean glycosyltransferases, i.e. GmSGT2 and GmSGT3. An essential conclusion from their study is that GmSGT3 transfers the rhamnosyl group from UDP-rhamnose not only to soyasaponin III but also to soyasaponin γa. In addition, studies have shown that the application of the allelopathic extract affected the biosynthesis of the main soyasaponin in seed, i.e. Soyasaponin. According to the Shibuya et al. 80 this compound is biosynthesized by the sequential addition of glucuronic acid, galactose, and rhamnose to soyasapogenol B. The results of research conducted by these authors, indicate that UDP-galactose (GmSGT2) and UDP-rhamnose (GmSGT3) are involved in the biosynthesis of soyasaponin I in soybean. Therefore, further research on allelopathic biostimulators is necessary, showing their influence on the activity of genes involved in the biosynthesis of these bioactive compounds in Glycine max. Although still scarce information is available on the molecular mechanisms underlying the glycosylation process, which play the key role in soyasaponin biosynthesis 81,82 . The results obtained in our study allow us to conclude that the allelopathic formulations induced the metabolic response of plants, manifested by changes in the quality and quantity of saponins in the seeds.
The exogenous application of extracts from Levisticum officinale triggered varied responses of soybean plants, depending on the method of their administration. Its extracts had high contents of polyphenolic compounds and rich micro-and macroelemental composition. The application of the allelopathic extracts modified soybean plant physiology, as manifested by changes in their biometric traits. Besides, soybean plants responded positively to this type of biostimulant by increased yield. The allelochemicals contained in the extracts modified the basic metabolic activities and biochemical indicators of soybean. Seeds obtained from the treated plants had higher contents of individual micro-and macroelements, as well as total concentrations of lipids (with a slight decrease in protein content). In addition, they featured changes in their amino acid profile and fatty acid composition. Analysis showed that extract did not contain growth regulators such as gibberellic acid and indole-3-acetic acid but the abscisic acid was found. The most important outcome of this study includes the increased contents of bioactive compounds, as the seeds produced by plants treated with the allelopathic biostimulant showed significantly increased concentrations of isoflavones and saponins. Considering the multi-level response of soybean plants, the foliar application of the analyzed allelopathic extracts is strongly recommended.
All observations made over the three-year study performed in real field conditions enable concluding that soybean plants responded positively to the biotic stress induced by the allelopathic biostimulant. www.nature.com/scientificreports/ To recapitulate, the natural biostimulants based on allelopathic extracts from Levisticum officinale may become a valuable tool in both the sustainable and organic agriculture. However, our research needs to be continued to fully exploit and elucidate their effects. The genetic approach would seem indispensable to identify the key enzymes and genes involved in plant response, which will, in turn, pave the way to the agricultural applications of these allelopathic biostimulants.

Methods
All methods were carried out in accordance with relevant guidelines and regulations. pdf). Infusions from lovage were prepared by the hot extraction method, i.e. 5 g of Levisticum officinale were added to 250 mL of distilled water. The solution was boiled in a water bath for 30 min. Then the solution was left in a dark place for 48 h, at the temperature 4 °C. Afterwards, extracts were centrifuged at 4250 rpm for 5 min and filtered through a blotted filter paper Whatman no. 1 84 . An active acidity (pH) was measured as well using a VOLTCRAFT KBM-110 m (Conrad Electronic SE, Hirschau, Germany) with a pH electrode.

Plant
Chemical composition of infusions from Levisticum officinale. Extracts from roots of Levisticum officinale, were determined for their chemical composition. Their contents of elements and sugars as well as composition and contents of phenolic compounds and phytohormones concentration were analyzed. Determination of mineral content in extract was assessed using the method of Zaguła et al. 85 by ICP-OES analysis. Identification and quantification of polyphenols in extract was studied exactly as previously described by Oszmiański et al. 86 . by the UPLC-PDA-MS method. Sugars in extracts from Levisticum officinale were evaluated based on the EN 12,630, 1999 standard and procedure of Pereira da Costa and Conte-Junior 87 , using the HPLC system. The phytohormones analysis in extracts from Levisticum officinale was assessed according to the procedure of Šimura et al. 88 , by the UPLC-PDA-TQD with an ESI method.
Plant yielding and seeds chemical composition determination. The seed yield and fat concentration in soybean seeds were determined. Total fat content was analyzed based on the acid hydrolysis method 89 . Determination of fatty acids in soybean seeds following the methods of Zhang et al. 90 , using gas chromatograph. Protein sequential fractionation was done using the procedure described by 91 . The protein fractions content was determined according to Bradford 92 using a calibration curve of albumins, globulins, prolamins, and glutenins of beans as standard protein. Hydrolysis of protein into amino acids has been carried out according to Davies and Thomas 93 , using HPLC analysis. Extraction and purification procedure for determination of isoflavones and saponins in soybean seeds was determined following the procedures of Kapusta et al. 94 . The isoflavones were identified by an UPLC Aquity system consisting of binary solvent manager, sample manager, column manager, PDA detector and tandem quadrupole mass spectrometer (TQD) with electrospray ionization mode (ESI) 95 . Saponins in soybean seeds were qualified and quantitated using the same system mentioned above, based on the method described by Jervis et al. 96 . The measurements of microelements and macroelements in soybean seeds were performed by ICP-OES (Inductively Coupled Plasma Optical Emission Spectrometers, Thermo iCAP Dual 6500, USA) 85 .