Next Article in Journal
Death Processes in Bovine Theca and Granulosa Cells Modelled and Analysed Using a Systems Biology Approach
Next Article in Special Issue
Overexpression of OsCASP1 Improves Calcium Tolerance in Rice
Previous Article in Journal
Botulinum Toxin: From Poison to Possible Treatment for Spasticity in Spinal Cord Injury
Previous Article in Special Issue
Peptide-Based Identification of Phytophthora Isolates and Phytophthora Detection in Planta
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 9
10.3390/ijms22094882
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea)

by
Efi Sarri
1,
Aikaterini Termentzi
2,
Eleni M. Abraham
3,
George K. Papadopoulos
1,
Eirini Baira
2,
Kyriaki Machera
2,
Vassilis Loukas
1,
Fotios Komaitis
4 and
Eleni Tani
1,*
1
Department of Crop Science, Laboratory of Plant Breeding and Biometry, Agricultural University of Athens, Iera Odos 75, 11855 Athens, Greece
2
Laboratory of Pesticides’ Toxicology, Department of Pesticides Control and Phytopharmacy, Benaki Phytopathological Institute, 8 St. Delta Street, Kifissia, 14561 Athens, Greece
3
Faculty of Agriculture, Forestry and Natural Environment, School of Forestry and Natural Environment, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
4
Department of Biotechnology, Laboratory of Molecular Biology, Agricultural University of Athens, Iera Odos 75, 11855 Athens, Greece
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(9), 4882; https://doi.org/10.3390/ijms22094882
Submission received: 5 April 2021 / Revised: 26 April 2021 / Accepted: 2 May 2021 / Published: 5 May 2021
(This article belongs to the Special Issue Plant Health and Food Security)
Browse Figures
Review Reports Versions Notes

Abstract

:
Increased soil salinity, and therefore accumulation of ions, is one of the major abiotic stresses of cultivated plants that negatively affect their growth and yield. Among Medicago species, only Medicago truncatula, which is a model plant, has been extensively studied, while research regarding salinity responses of two important forage legumes of Medicago sativa (M. sativa) and Medicago arborea (M. arborea) has been limited. In the present work, differences between M. arborea, M. sativa and their hybrid Alborea were studied regarding growth parameters and metabolomic responses. The entries were subjected to three different treatments: (1) no NaCl application (control plants), (2) continuous application of 100 mM NaCl (acute stress) and (3) gradual application of NaCl at concentrations of 50-75-150 mM by increasing NaCl concentration every 10 days. According to the results, M. arborea maintained steady growth in all three treatments and appeared to be more resistant to salinity. Furthermore, results clearly demonstrated that M. arborea presented a different metabolic profile from that of M. sativa and their hybrid. In general, it was found that under acute and gradual stress, M. sativa overexpressed saponins in the shoots while M. arborea overexpressed saponins in the roots, which is the part of the plant where most of the saponins are produced and overexpressed. Alborea did not perform well, as more metabolites were downregulated than upregulated when subjected to salinity stress. Finally, saponins and hydroxycinnamic acids were key players of increased salinity tolerance.

1. Introduction

Soil salinization is a major problem facing agriculture and contributes to reducing productivity of arable lands mainly in arid and semi-arid areas [1]. According to FAO 1997 [2], soil with an electrical conductivity of the saturation extract (ECe) more than 4 dS m−1 is defined as saline. The problem is steadily increasing and it is predicted that more than 50% of the arable land is going to be salinized by 2050 [3] as a consequence of the expected increase of irrigated land as well as due to the expected global climate change [4]. The effect of salinity on crops depends on species, cultivars, environmental factors and plant growth stage [5]. The loss of crop productivity in saline soil is a consequence of the negative impact of salinity on plant growth and development [6].
The depression of plant growth under high levels of NaCl concentrations is a result of the reduction of water potential (hyperosmotic stress) and the accumulation of ions to a toxic level for their growth (hyperionic stress) [7]. Plants respond to these stresses by modifying their physiological and biochemical processes. Thus, researchers have focused on plants’ signaling perception and adaptation to an adverse environment by monitoring gene regulation, expression of functional proteins and accumulation of small molecules (i.e., metabolites) [5,8,9]. Several studies have implicated the role of primary metabolites such as sugars, amino acids and organic compounds in osmotic regulation. On the other hand, secondary metabolites being the end products of primary metabolites are more species specific and are associated with plant protection by having multiple functions (e.g., acting as antioxidants, reactive oxygen species (ROS) scavengers, regulatory molecules) [10]. Although the synthesis of secondary metabolites such as phenols, saponins, flavonoids, carotenoids, lignins, etc. generally increases when plants are cultivated under salt stress conditions [11], special attention has been given only to afew target compounds in previous studies [12,13]. Thus, the association of antioxidant activity and accumulation of phenolic compounds under salt stress has been confirmed by several studies [14,15]. Other secondary metabolites which were found in high abundance and are possibly involved in salt tolerance by rendering a more efficient antioxidant capacity are flavonoids [16], polyphenols, tannins, and anthocyanins [17]. In a recent work, Yang et al. (2017) [18] used saponin as a priming agent to improve the seed germination of quinoa plants in saline conditions due to the fact that saponins can act as ROS scavengers. In this regard, by analyzing the changes of metabolites at the whole-metabolome scale [19] and by using these metabolic profile changes in combination with other ‘omic’ analyses, such as genome, transcriptome and proteome analysis, one can elucidate the regulatory networks and identify biomarkers that regulate stress responses and are useful in plant improvement [20,21].
Medicago sativa L. (lucerne, alfalfa) (M. sativa) is a perennial legume of high importance in animal feeding. It is the main fodder crop in many countries worldwide for its high nutritive value as well as its contribution to soil fertility due to nitrogen fixation. Beyond its use as animal feed, M. sativa is also important as cover crop [22] and as a source of biologically active compounds (secondary metabolites) for the agriculture and pharmaceutical industries [23]. On the other hand, Medicago arborea L. (tree medick) (M. arborea) is an evergreen shrub that mainly spreads in the Mediterranean basin and it is used in animal feeding mostly during summer and winter. It is considered the oldest species of the genus and is more tolerant to abiotic stresses compared to alfalfa [24,25]. The Alborea hybrid was created by crossing the species M. sativa and M. arborea in order to meet the great needs for animal feed and high biomass production as well as resilience to environmental changes [26].
There are only a few studies on the metabolite profile of the Medicagos species including the model legume M. truncatula. In particular, Barsch et al. (2006) [27] determined the characteristic metabolites for M. truncatula root nodules. Recently, pathway-specific metabolome analysis using 18O2-labeled M. truncatula was reported by Kera et al. (2018) [28], while Dickinson et al. (2018) [29] combined transcriptomic techniques and k-means clustering in metabolomics to detect biotic and abiotic stress markers for M. truncatula. There is only one work regarding the changes of primary metabolites in M. sativa when subjected to salt stress treatment [30]. To our knowledge, there is no research regarding the modification of secondary metabolite profiles under salt stress of Medicagos. In this regard, the aim of the present study was the identification of the secondary metabolite profile in shoots and roots of M. sativa, M. arborea and their hybrid Alborea under salinity conditions, after application of two salt treatments: (a) continuous application (acute stress) of 100 mM NaCl, and (b) a gradual increase in salt concentration up to 100 mM.

2. Results

2.1. Growth Parameters and Salinity Sensitivity Index

The studied growth parameters of seedlings (stem height, stem elongation rate, salinity sensitivity index) were significantly affected by salt treatments, entries and dates of treatments (Table 1). The only exception was the stem height and the stem elongation rate (SER) that were not affected by the entries. In all cases, the interaction between salt treatments and entries was significant, indicating different responses of the entries to salt treatments (Table 1). The stem height and SER of the hybrid (Al) were significantly reduced under both salt stress of 50-75-100 mM and salt shock of 100 mM (Figure 1). Regarding the parental species, the stem height and SER of M. sativa were significantly reduced only under salt shock of 100 mM (Figure 1). On the other hand, the stem height and SER of M. arborea tended to decrease under both salt stress and shock, but this did not show statistical differences (Figure 1). It has to be mentioned that the stem height and SER of the Alborea population were the highest under control while both parameters did not significantly differ among the entries under salt stress. However, the population of M. arborea had the highest stem height and SER under salt shock (Figure 1). As a result, the highest salinity sensitivity index (absolute number) was recorded for the population of Alborea, while the lowest was recorded under salt stress for the population of M. sativa and under salt shock for the population of M. arborea (Figure 2).

2.2. Metabolomics Analysis

2.2.1. Comparison among the Entries

Regarding the secondary metabolic profiling, PCA analysis clearly demonstrated that M. arborea differed from the other two entries (M. sativa and Alborea) by comparing either all treatments or each individual treatment between entries in both shoots and roots (Figure 3, R1-S3). On the other hand, the separation of the Alborea and M. sativa samples in all categories compared was much more blurred and in some cases the permutation tests showed no differentiation of the two groups. We drew the same conclusions by conducting PLS-DA. Permutation tests, allowing 100 permutations, indicated good predictability and goodness offit for PLS-DA models in all cases (details are omitted).
Specifically, for all treatments (Figure 3, R1 & S1), it was clearly shown that the PC1 axis, which is a linear combination of secondary metabolites for all treatments and interprets 38.3% and 36.5%, respectively, of the total variability of the phenomenon, differentiated the species as follows: in one group M. arborea and in another group M. sativa and Alborea. The metabolites that differentiated M. arborea from the other entries were mainly saponins (e.g., medicagenic acid, caryophyllogenic type), (Supplementary Table S1).
However, at 100 mM NaCl (Figure 3, R2, S2) and at 50-75-100 mM NaCl a tendency to distinguish the metabolic profile of M. sativa from Alborea was evident. More specifically, for the gradual stress treatment of 50-75-100 mM NaCl for the shoots, phenolic compounds and some triterpenic saponins were highly induced for all species. Nevertheless, for M. arborea, hydroxycinnamic acids, a class of phenolic compounds, were overexpressed. Moreover, some flavonoids and saponins were upregulated as well. On the other hand, M. sativa and Alborea upregulated mainly different triterpenic saponins and flavonoids, compared to M. arborea, (Supplementary Table S1). For the roots, under gradual salinity stress, the VIPs for M. arborea were exclusively saponins (mainly triterpenic) while for M. sativa and Alborea, apart from saponins, some phenols and flavonoids were also detected (Table 2). It is noteworthy that the secondary metabolite that mostly characterized the Alborea roots extracts under gradual salinity was tentatively identified as a diphenyl, most probably a lignan of eudesmin type, and was absent from M. arborea.
In general, under acute stress, the VIPs at the shoots for the M. arborea metabolomic profile were like those of gradual salinity stress, which implies a similar way of reaction. In particular, M. arborea plants were still the only ones that produced hydroxycinnamic acids (caffeic, ferulic, cinnamic acids). Those included mainly flavonoid glycosides. However, under acute stress, M. sativa significantly differed from Alborea by overexpressing some unique compounds such as a flavonoid compound, probably an apigenin derivative, and two saponins. On the other hand, Alborea overexpressed a special category of phenols, a lignan, that was also abundant in the roots (Table 3).
Concerning the metabolites accumulated in roots under acute stress, in M. arborea many triterpenic saponins, that are species specific, mainly glycosides of the aglycons zahnic and medicagenic acids, were present. On the other hand, in M. sativa roots under acute stress, no significant accumulation of secondary metabolites was detected, whereas Alborea reacted under acute stress by increasing its characteristic lignans that have been previously mentioned (Table 3). In general, it seems that under acute stress, M. sativa overexpressed saponins in the shoots while M. arborea overexpressed saponins in the roots. Alborea, on the other hand, displayed an increased production of lignans and phenyl tetrahydrofurans mainly in the roots.
As mentioned above, the Alborea and M. sativa species roughly differed when it came to their secondary metabolite content in both shoots and roots. The permutation tests showed only small differentiations presented in Supplementary Table S2. The fact that M. arborea expressed different metabolites compared to the other two entries was also shown by the heat maps (Figure 4 and Figure 5). Thus, the metabolites that were upregulated in M. arborea (in blue), were downregulated in the other two entries (in red) and vice versa.
Specifically, for Figure 4 (roots), it was evident that M. arborea in both stressed and control plants up- or downregulated roughly the same group of metabolites. For Alborea and M. sativa no clear trend was detected. However, for M. sativa, more metabolites were accumulated in gradual stress than in acute stress compared to control. The same results are captured in Figure 5.

2.2.2. Comparison between Roots and Shoots

A distinct set of secondary metabolites was expressed in the roots compared to shoots (Figure 6). More specifically, for the M. arborea and Alborea, a clear difference between the metabolites produced in the shoots compared to roots was demonstrated, regardless of the treatment (Figure 6a,c). For M. sativa, however, in addition to a clear separation of metabolites between roots and shoots, a differentiation of secondary metabolites was apparent, even for each treatment separately (Figure 6b).
From the comparison of the Table 2a,b and Table 3a,b, it is obvious that roots and shoots of the same treatments produce different secondary metabolites that characterize salinity response. In general, roots are the parts of the plants where most of the saponins are produced and overexpressed, especially under salinity stress. Specifically, for M. arborea species, saponins were almost exclusively found in roots and this profile is even more pronounced under the salinity treatment. On the other hand, in the shoots other types of compounds, mainly phenols, were produced. Hydroxycinnaminc acids were exclusively found in the shoots. Flavonoid content was also more intense in the shoot extracts. Lignans on the other hand were equally characterized in the Alborea extracts in shoots and roots. Saponins in the aerial parts were usually less intensely expressed as peaks in the extract and sometimes they differed slightly, mainly in the glycosylation level.

2.3. Secondary Metabolomic Changes in Response to Salt Treatment

As shown in Figure 7 and Supplementary Tables S3 and S4, there were differences in the expression of secondary metabolites, in some cases, between control plants and plants subjected to salinity stress for each entry. In roots of M. arborea, a reduction in some saponins, flavonoids and triterpenic acids in plants subjected to 100 mM NaCl and 50-75-100 mM NaCl, compared to control plants, was detected. On the other hand, the expressions of benzyl tetrahydrofurans, lignans and phenols were significantly higher in the plants under gradual stress compared to the controls. For the other treatments, there were no significant differences in roots or shoots. According to our results, M. arborea altered the accumulation of secondary metabolites when subjected to salinity stress, only in roots and not in shoots.
For M. sativa, flavonoids and phenols were decreased compared to control plants, while lignans, benzyl tetrahydrofurans and phenols were increased in roots. In the shoots, biphenyl and flavonoids were decreased, while benzyl tetrahydrofuran, phenol and lignans were increased in plants receiving 50-75-100 mM NaCl compared to control plants.
Finally, Alborea plants, which were subjected to gradual stress, showed higher levels of a specific flavonoid in roots compared to control plants. On the other hand, the comparison between control plants and plants subjected to acute stress showed reduced levels for saponins and triterpenic acids while an increase was observed for benzyl tetrahydrofuran and lignans for plants receiving 100 mM NaCl. Regarding the shoots, there was an overexpression of benzyl tetrahydrofuran, triterpenic acid, flavonoids and lignans in the plants that were subjected to acute stress compared to the control plants, while on the contrary there was a decrease in phenol accumulation in the plants that received 100 mM NaCl compared to control. For the shoots of Alborea plants that were subjected to gradual stress of 50-75-100 mM NaCl, overexpression of benzyl tetrahydrofurans and lignans were observed in relation to the control plants.

3. Discussion

The accumulation of ions, due to increased salinity in the soil, is one of the major abiotic stresses of cultivated plants that negatively affect their productivity. Among Medicago species, the responses to salt stress of Medicago truncatula, which is a model plant, and Medicago sativa L., have been extensively studied [31,32] and to a lesser extent so have the responses of Medicago arborea [25,33,34]. Nevertheless, research of the metabolomic profile of the Medicago species upon salt stress has been very limited. In the present study the comparative alterations of the secondary metabolomic profile in roots and shoots of M. sativa, M. arborea and their hybrid Alborea were investigated under acute salinity stress of 100 mM NaCl and gradual increase in salt concentration up to 100 mM NaCl.
According to the results of the present study, the hybrid in terms of seedling growth performed better compared to parental species under well-growing conditions but it was more sensitive under salinity stress. On the other hand, regarding the seedling growth of the parental species, M. sativa and M. arborea were not affected by the gradual application of NaCl but reduced about 20% and 10%, respectively, in their growth under acute salinity stress. These results are consistent with the results of our previous experiment [25] as well as the results of a preliminary experiment in order to confirm their performance under salt conditions.
The secondary metabolic profiles of M. arborea in shoots and roots were clearly differentiated from that of M. sativa and their hybrid under both control and salinity stress conditions. Generally, saponins were the main metabolites that discriminated M. arborea from M. sativa and their hybrid. In particular, under salinity stress M. arborea exclusively accumulated in the roots saponins such as medicagenic acid and zahnic acid while in the shoots had few specific flavonoids and phenolic compounds (e.g., hydroxycinnamic acid). On the other hand, the metabolic profile of the hybrid was more similar to that of M. sativa. There was a clear differentiation between them mainly under acute salinity stress conditions where M. sativa upregulated few unique flavonoid and saponin compounds in the shoots while Alborea produced lignans in the roots. It was noteworthy that Alborea overexpressed a characteristic lignan in both shoots and roots under stress conditions while M. sativa did not accumulate any secondary metabolite in the roots under acute salinity stress.
Accumulation of secondary metabolites is of particular interest, as they usually regulate responses to environmental stressors, such as high and low temperatures, drought, salinity, pathogen attack etc. [35,36,37]. Very recent findings unravel the multifunctional role of secondary metabolites, either as plant growth regulators or as primary metabolites in a broad sense [36].
Among the Medicago species, the research about the secondary metabolites has mainly focused on M. sativa and Μ. truncatula. On the other hand, information on the profile of Μ. arborea secondary metabolites is limited. The majority of the studies emphasized saponins and flavonoids which are of highest interest because of their pharmaceutical use [38]. Saponins of Medicago species are triterpenoid compounds and more than thirty-three different saponins containing one or more sugar chain units have been identified in M. sativa. Among them, medicagenic acid, hederagenin, zahnic acid and soyasapogenols A and B are the main saponins of M. sativa aerial parts [39,40,41]. Nevertheless, the content of medicagenic acid in several cases is higher in roots than in aerial parts [42]. This was the case in the present study. The accumulation of the aforementioned metabolites was characteristic of M. arborea roots rather than M. sativa aerial parts.
Accumulation of such compounds has been associated with a plethora of biological activities, including seed germination, allelopathy, poor digestibility and beneficial anti-pathogen (including fungi and insects) properties [43]. Moreover, M. arborea has high allelopathic potential that affects the growth and germination of other species [44]. It is possible that this characteristic of M. arborea is associated with the accumulation of these compounds in its roots in the present study. Furthermore, the presence of the aforementioned metabolites and especially zahnic acid in shoots of Medicago is related tolow palatability of forages [23]. The M. sativa population used in the present study comes from breeding and as a consequence it contains a low amount of anti-nutritional factors in its aerial parts.
Although reports in Medicago species relate saponins to biotic stress tolerance, several reports in other species highlight their importance in abiotic stress responses [45,46]. Special attention should be given to triterpenoid compounds and their accumulation upon abiotic stresses, as they may contribute to fortification of plant membranes. Moses et al., 2014 and Pulvento et al., 2012 [43,47] studied the response of quinoa plants (a salt-tolerant species) under conditions of increased salinity and found that the plants had accumulated high levels of saponins compared to control plants. Another study that correlated saponins with plant salinity tolerance was conducted by Oku et al., 2003 [48] in which the concentration of triterpenoids increased significantly under high salinity in both shoots and roots of Kandeliacandel and Bruguieragymnorrhiza.
As previously mentioned, Μ. arborea changed the profile of secondary metabolites in shoots mainly by over producing hydroxycinnamic acids. The three polyphenols (caffeic acid, ferulic acid and cinnamic acid) belong to hydroxycinnamic acids (HCAs) and are all members of the biochemical pathway leading to lignin formation. Additionally, ferulic acid, apart from acting as an ROS scavenger, can positively regulate the activity of other scavenging enzymes and can also inhibit the expression of enzymes involved in ROS generation [39,40].
Regarding secondary metabolites, changes in each entry in relation to its control, M. sativa, deeply altered the accumulation of metabolites in both shoots and roots (mainly phenolic compounds), to face the harmful effects caused by salt stress. On the other hand, M. arborea showed modifications in the levels of few metabolites, mainly phenolic compounds (phenol tetrahydrofurans and lignans), only in the roots while no modification of secondary metabolites took place in shoots. Interestingly, both species upregulated these compounds under gradual NaCl application. Alborea upregulated only two groups of secondary metabolites and downregulated several other groups under acute stress and only one metabolite under gradual stress in roots (Supplementary Figures S1 and S2). It seems that the root system changes in metabolomic profile are key players in salt tolerance. Roots contact the soil directly and should react immediately to avoid the detrimental effects of soil salinity [49]. When plants are subjected to salt stress, they can resist or reduce the damage caused by increasing osmotic adjustment, maintaining cell wall integrity and improving ROS scavenging [50,51,52]. Several studies have highlighted that phenolic compounds maintain high antioxidant activity and are involved in ROS detoxification and prevention from oxidative damage under salt stress [53,54,55]. More specifically, lignans are involved in lignification and cell wall synthesis and they exhibit antioxidant capacity, thus they have a positive role in stress tolerance [39,56,57]. It is known that lignans have remarkable ecological functions in plants, providing protection against herbivores and microorganisms [58]. However the accumulation of lignans in roots of Alborea did not provide the expected tolerance to salinity stress.
Finally, secondary metabolites such as phenolics and flavonoids might act as hormones that ultimately elicit an improved salt tolerance [59]. In several studies, salt tolerance was tightly linked to higher accumulation of phenolic compounds [60,61,62]. A very popular strategy to induce the production of several polyphenols is by application of stress hormones [39,43]. Given that M. arborea’s seedling growth was superior compared to the other two entries under acute stress it seems that its production of secondary metabolites, and in particular the production of saponins and polyphenols, had a positive effect under salt stress conditions without metabolic costs. Plants can minimize the yield costs of secondary metabolite production either by controlled recycling of the resulting compounds or by evolving the use of same biosynthetic machinery for multifunction purposes [36]. Gaining information considering the metabolomic pathways under salinity stress of both parents (and especially M. arborea) and the ways that they recycle or use the produced secondary metabolites not only to combat salt stress but for multiple purposes may be very valuable for genetic engineering Alborea to breed a more tolerant and high yielding offspring.

4. Materials and Methods

4.1. Plant Material

The M. arborea parents used to develop Alborea were originally collected from Greece. The M. sativa parent population used in this study was a hybrid of the two M. sativa parents. Two M. sativa parents were used in crosses of M. sativa × M. arborea to produce 27 initial Alborea hybrids [26]. The Alborea population used in this study was developed from intercrosses of 20 initial hybrids using methods described by Irwin et al. (2015) [63]. Seeds were obtained from E. Bingham, Agronomy Department, Univ. Wisconsin-Madison, USA.

4.2. Seed Pretreatment

Seeds of M. sativa, M. arborea, and Alborea were scarified with absolute sulfuric acid for 10 min and then rinsed thoroughly with sterile distilled water. Bleach solution (3%) was added for 1.5–2 min and then rinsed thoroughly. The scarified seeds were placed on petri dishes containing 0.5% agar at 4 °C overnight, and then transferred to 20 °C (dark) for 3–4 days.

4.3. Growth Conditions

From each species and their hybrid (hereafter entries), 60 seeds were planted, after pre-germination, in pots of 8.5 cm height and 10 cm in diameter. Each pot contained peat, Kronos N 50–300 mg/L, P2O5 80–300 mg/L, K2O 80–300 mg/L, pH 5–6.5, salinity <1.75 g/L and one seed. Then all the pots were placed in a growth chamber at a constant temperature of 25 °C, and a photoperiod of 16 h of light and 8 h of darkness, where they were regularly watered and fertilized each week for a month during the salinity treatments. For each treatment, 20 plants were used for each species. The completely randomized block design was used (block to treatments).

4.4. Salt Stress Treatments

The plant material was subjected to three different treatments: (1) no NaCl application (control plants), (2) continuous application of 100 mM NaCl (acute stress) and (3) gradual application of NaCl at concentrations of 50-75-100 mM with increased concentration every 10 days.

4.5. Growth Characteristics Measurements

The length of the stems from the base to the tip was measured in each seedling at an interval of three days through the duration of the salt treatments. Stem elongation rate (SER) was estimated as SER = (T2 − T1)/t, where T1 and T2 represent the stem length at the beginning and at the end of a time t, respectively. Additionally, the salinity sensitivity index (IS) based on the stem length was estimated according to the formula proposed by [64], IS = (Hs − Ht)/Ht × 100, in which Hs and Ht represent the values of stem length of the salt-stressed and control plants, respectively.

4.6. Metabolites Extraction

The plant tissue was ground with nitrogen for the extraction of metabolites. Then, in 20 mg of plant tissue, 395 μL of methanol and 5 μL of ribitol were added. The samples were incubated for 15 min at 70 °C with continuous shaking and sealed with parafilm. Next, 200 μL of chloroform was added, followed by vortex. The samples were then maintained at 37 °C for 5 min and continued shaking. Then, 400 μL of water was added and then vortexed. Samples were centrifuged for 5 min at 13,000 rpm and 250 μL of the supernatant phase was transferred to a new Eppendorf. Finally, the samples were evaporated in nitrogen gas.

4.7. Metabolomic Analysis

The high-resolution mass spectrometry analyses of the extracts for the dereplication and the metabolomics experiments were performed on an UHPLC-HRMS/MS (Ultra High-Performance Liquid Chromatography–High-Resolution Mass Spectrometry) Orbitrap Q-Exactive platform (Thermo Scientific, Erlangen, Germany). The ultra-high performance liquid chromatography system was a Dionex Ultimate 3000 UHPLC system (Thermo Scientific, Germany). A Hypersil Gold UPLC C18 (2.1 × 100 mm, 1.9 μm) reversed phased column (Thermo Scientific, Germany) was used for the separations.
The mobile phase consisted of solvents A: aqueous 0.1% (v/v) formic acid and B: acetonitrile. Both formic acid and acetonitrile were purchased from Merck (Darmstadt, Germany).
Different gradient elutions were performed for positive and negative ion mode detection and after optimization of the chromatography the gradient applied was: T = 0 min, 5% B; T = 3 min, 5% B, T = 21 min, 95% B, T = 23 min, 95% B, T = 24 min, 5% B; T = 30 min, 5% B. The flow rate was 0.220 mL/min and the injection volume 5 μL. The column temperature was kept at 40 °C while the sample tray temperature was set at 10 °C.
The ionization was performed at HESI (Heated ElectroSpray Ionization), at both positive and negative modes. The conditions for the HRMS for both negative and positive ionization modes were set as follows: capillary temperature, 350 °C; spray voltage, 2.7 kV; S-lense Rf level, 50 V; sheath gas flow, 40 arb. units; aux gas flow, 5 arb. units; aux. gas heater temperature, 50°C. Analysis was performed using the Fourier transform mass spectrometry mode of the LTQ orbitrap (FTMS) in the full scan ion mode, applying a resolution of 70,000, while acquisition of the mass spectra was performed in every case using the centroid mode. The data-dependent acquisition capability was also used at 35,000 resolution, allowing for MS/MS fragmentation of the three most intense ions of every peak exceeding the predefined threshold applying a 10 s dynamic exclusion. Stepped normalized collision energy was set at 40, 60 and 100. Data acquisition and analysis was completed employing Xcalibur 2.1. The dereplication was performed by Compound Discoverer software (Thermo) where spectral libraries of natural products were used (m/z cloud, in house libraries).
The analyses of the samples followed a classical metabolomic sequence where all samples were triplicated. Quality control (QC) samples were created as pooled samples and were injected periodically, every 6 samples, throughout the run to assess instrument stability.
Following the overall chemical profiling of samples, compounds responsible for group clustering were further investigated. Variable Importance in Projection (VIPs) and coefficient values extracted from the Partial Least Squares-Discriminant Analysis (PLS-DA) revealed the most important features in positive ESI (ElectroSpray Ionization) mode, responsible for the differentiation of several comparisons that were performed.
The structural elucidation of features responsible for sample clustering was performed by comparison of the chromatographic and spectrometric features of each peak in comparison with data from literature and the standard compounds. The Orbitrap HR-UPLC-MS2 profiling was performed in both positive and negative modes. The high resolving power (70,000 at the full scan experiments and 35,000 for the MS/MS fragments) of the Q-Exactive Orbitrap analyser in correlation to the accurate mass measurements (m < 1 ppm for both full scan and MS/MS ions) assured the identification of the VIPs compounds responsible for the differentiation of the two groups with high confidence. The suggested EC (Elemental Composition) for molecular ions and MS/MS fragments, as well as the respective RDBeq (Ring Double Bond Equivalents) further assisted the safe identification process. However, due to the great similarity of the compounds found in the extracts, in several cases the structural elucidation was not feasible. Instead, a suggested annotation is given, while the general category in which the compound belongs is assigned, beyond any question.

4.8. Statistical Analysis

Differences between the entries for growth parameters and salinity sensitivity index were conducted using analysis of variance (ANOVA) followed by Tukey’s means post-comparison tests at 0.05 significance level. Statistical analysis was performed using IBM SPSS Statistics 23 for Windows (SPSS Inc. Chicago, IL, USA).
Spectrometric and chromatographic data were imported to Compound Discoverer 2.1 (Thermo Scientific) prior to statistical analysis. The appropriate set up for peak detection, deconvolution, deisotoping, alignment and gap filling procedures was applied. For the gap filling the signal to noise threshold was set at 1.5. The normalization was QC-based and the regression model used was linear. Blanks were excluded and the type was constant mean. The appropriate settings were also used for the composition prediction, setting the mass tolerance at 5 ppm. Principal component analysis (PCA) was used to group both species and their hybrids based on secondary metabolites and heatmaps were used to investigate expression differences of secondary metabolites.
The MS (mass spectrometry) data were also subjected to multivariate statistical analysis (PLS-DA) using SIMCA P + 11.5 software (Umetrics, Umea, Sweden). Permutation testing employing 100 random permutations was applied to validate the statistical models. The optimal number of principal components for all models resulted from the R2 and Q2 values. The variables that exhibited a VIP scoring greater than 1 were verified by t-test using a p-value at ≤0.05. These variables were annotated using the Compound Discoverer software in correlation with in-house libraries of natural products, applying m/z to lerance of 5 ppm and taking into consideration the isotopic and MS/MS fragmentation pattern. Heatmaps were constructed by using XLSTAT software.

5. Conclusions

According to our results the parental lines responded better compared to their hybrid under salinity stress conditions. Regarding metabolomic analysis, M. arborea upregulated saponins in the roots under both gradual and acute stress while in the shoots mostly phenolic compounds (hydroxycinnamic acids) were upregulated. Thus, these secondary metabolites may play a role in M. arborea’s stable performance under salinity stress (Supplementary Table S5). These findings are useful for future studies to decipher the complex plant metabolome in order to improve abiotic stress tolerance especially in such high-importance legumes. Lastly, future breeding efforts should focus on engineering biochemical pathways leading to the specific secondary metabolite synthesis (either phenolic compounds or saponins) specifically in the roots of Alborea plants.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22094882/s1.

Author Contributions

Conceptualization, E.T. and E.M.A.; methodology, E.S., V.L., F.K.; software, A.T.; validation, E.T., E.M.A., E.S., G.K.P., A.T., E.B.; formal analysis, E.T., E.M.A., A.T., G.K.P.; investigation, E.T.; resources, E.T.; data curation, A.T., G.K.P.; writing—original draft preparation, E.T., E.M.A., E.S., A.T.; writing—review and editing, E.T., E.M.A., E.S., A.T., G.K.P., E.B., K.M.; visualization, E.T., E.M.A., E.S., A.T.; supervision, E.T.; project administration, E.T.; funding acquisition, E.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the MSc program of Agricultural University of Athens.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting reported results can be found at: https://drive.google.com/drive/folders/1j4lARpnnTZjdtA-h2ihUHqQRa2Y4VOql?usp=sharing (accessed on 4 May 2021).

Acknowledgments

This work has been partially financed by the European Project Ethnoherbs, (Η2020-MSCA-RISE-2018-823973). We would like to thank E. Bingham for providing us the genetic material used for the study and A. Tornesaki for her assistance during metabolites extraction. Finally we would like to thank Ioanna Kendrick for English editing.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Pitman, M.G.; Läuchli, A. Global impact of salinity and agricultural ecosystems. In Salinity: Environment-Plants-Molecules; Springer: Dordrecht, The Netherlands, 2004; pp. 3–20. [Google Scholar]
  2. FAO. State of The Worlds Forests; FAO: Rome, Italy, 1997. [Google Scholar]
  3. Jamil, A.; Riaz, S.; Ashraf, M.; Foolad, M.R. Gene expression profiling of plants under salt stress. Crit. Rev. Plant Sci. 2011, 30, 435–458. [Google Scholar] [CrossRef]
  4. Shrivastava, P.; Kumar, R. Soil salinity: A serious environmental issue and plant growth promoting bacteria as one of the tools for its alleviation. Saudi J. Biol. Sci. 2015, 22, 123–131. [Google Scholar] [CrossRef] [Green Version]
  5. Chandna, R.; Azooz, M.; Ahmad, P. Recent advances of metabolomics to reveal plant response during salt stress. Salt Stress Plants Signal. Omics Adapt. 2014, 1–14. [Google Scholar] [CrossRef]
  6. Guo, R.; Shi, L.; Yang, Y. Germination, growth, osmotic adjustment and ionic balance of wheat in response to saline and alkaline stresses. Soil Sci. Plant. Nutr. 2009, 55, 667–679. [Google Scholar] [CrossRef]
  7. Apse, M.P.; Aharon, G.S.; Snedden, W.A.; Blumwald, E. Salt tolerance conferred by overexpression of a vacuolar Na+/H+ antiport in Arabidopsis. Science 1999, 285, 1256–1258. [Google Scholar] [CrossRef] [PubMed]
  8. Isah, T. Stress and defense responses in plant secondary metabolites production. Biol. Res. 2019, 52, 39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Zhu, J.-K. Abiotic stress signaling and responses in plants. Cell 2016, 167, 313–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Lucini, L.; Borgognone, D.; Rouphael, Y.; Cardarelli, M.; Bernardi, J.; Colla, G. Mild potassium chloride stress alters the mineral composition, hormone network, and phenolic profile in artichoke leaves. Front. Plant Sci. 2016, 7, 948. [Google Scholar] [CrossRef] [Green Version]
  11. Ramalingam, A.; Kudapa, H.; Pazhamala, L.T.; Weckwerth, W.; Varshney, R.K. Proteomics and Metabolomics: Two emerging areas for legume improvement. Front. Plant Sci 2015, 6, 1116. [Google Scholar] [CrossRef] [Green Version]
  12. Borgognone, D.; Cardarelli, M.; Rea, E.; Lucini, L.; Colla, G. Salinity source-induced changes in yield, mineral composition, phenolic acids and flavonoids in leaves of artichoke and cardoon grown in floating system. J. Sci. Food Agric. 2014, 94, 1231–1237. [Google Scholar] [CrossRef]
  13. Colla, G.; Rouphael, Y.; Cardarelli, M.; Svecova, E.; Rea, E.; Lucini, L. Effects of saline stress on mineral composition, phenolic acids and flavonoids in leaves of artichoke and cardoon genotypes grown in floating system. J. Sci. Food Agric. 2013, 93, 1119–1127. [Google Scholar] [CrossRef]
  14. Ksouri, R.; Megdiche, W.; Debez, A.; Falleh, H.; Grignon, C.; Abdelly, C. Salinity effects on polyphenol content and antioxidant activities in leaves of the halophyte Cakile maritima. Plant. Physiol. Biochem. 2007, 45, 244–249. [Google Scholar] [CrossRef] [PubMed]
  15. Waskiewicz, A.; Muzolf-Panek, M.; Golinski, P. Phenolic Content Changes in Plants under Salt Stress. In Ecophysiology and Responses of Plants under Salt Stress; Springer: New York, NY, USA, 2013; pp. 283–314. [Google Scholar]
  16. Edreva, A.; Velikova, V.; Tsonev, T.; Dagnon, S.; Gürel, A.L.; Aktas, L. Stress-protective role of secondary metabolites: Diversity of functions and mechanisms. Gen. Appl. Plant. Physiol. 2007, 34, 67–78. [Google Scholar]
  17. Ramakrishna, A.; Ravishankar, G.A. Influence of abiotic stress signals on secondary metabolites in plants. Plant. Signal. Behav 2011, 6, 1720–1731. [Google Scholar] [CrossRef]
  18. Yang, A.; Akhtar, S.; Iqbal, S.; Qi, Z.; Saddiq, M.; Jacobsen, S.-E. Saponin seed priming improves salt tolerance in quinoa. J. Agron. Crop. Sci. 2017, 204. [Google Scholar] [CrossRef]
  19. Fukushima, A.; Kusano, M. Recent progress in the development of metabolome databases for plant systems biology. Front. Plant Sci 2013, 4, 73. [Google Scholar] [CrossRef] [Green Version]
  20. Cai, Z.; Chen, H.; Chen, J.; Yang, R.; Zou, L.; Wang, C.; Chen, J.; Tan, M.; Mei, Y.; Wei, L.; et al. Metabolomics characterizes the metabolic changes of Lonicerae Japonicae Flos under different salt stresses. PLoS ONE 2020, 15, e0243111. [Google Scholar] [CrossRef]
  21. Jha, U.C.; Bohra, A.; Jha, R.; Parida, S.K. Salinity stress response and ‘omics’ approaches for improving salinity stress tolerance in major grain legumes. Plant. Cell Rep. 2019, 38, 255–277. [Google Scholar] [CrossRef] [PubMed]
  22. Fageria, N.K.; Baligar, V.C.; Bailey, B.A. Role of cover crops in improving soil and row crop productivity. Commun. Soil Sci. Plant. Anal. 2005, 36, 2733–2757. [Google Scholar] [CrossRef]
  23. Rafińska, K.; Pomastowski, P.; Wrona, O.; Górecki, R.; Buszewski, B. Medicago sativa as a source of secondary metabolites for agriculture and pharmaceutical industry. Phytochem. Lett. 2017, 20, 520–539. [Google Scholar] [CrossRef]
  24. Lefi, E.; Conesa, M.À.; Cifre, J.; Gulías, J.; Medrano, H. Dry matter allocation in Medicago arborea and Medicago citrina in response to drought and defoliation. Crop. Pasture Sci. 2012, 63, 179–189. [Google Scholar] [CrossRef]
  25. Tani, E.; Sarri, E.; Goufa, M.; Asimakopoulou, G.; Psychogiou, M.; Bingham, E.; Skaracis, G.N.; Abraham, E.M. Seedling growth and transcriptional responses to salt shock and stress in Medicago sativa L., Medicago arborea L., and their hybrid (Alborea). Agronomy 2018, 8, 231. [Google Scholar] [CrossRef] [Green Version]
  26. Bingham, E.; Armour, D.; Irwin, J. The hybridization barrier between herbaceous Medicago sativa and woody M. arborea is weakened by selection of seed parents. Plants 2013, 2, 343–353. [Google Scholar] [CrossRef] [Green Version]
  27. Barsch, A.; Carvalho, H.G.; Cullimore, J.V.; Niehaus, K. GC-MS based metabolite profiling implies three interdependent ways of ammonium assimilation in Medicago truncatula root nodules. J. Biotechnol. 2006, 127, 79–83. [Google Scholar] [CrossRef] [PubMed]
  28. Kera, K.; Fine, D.D.; Wherritt, D.J.; Nagashima, Y.; Shimada, N.; Ara, T.; Ogata, Y.; Sumner, L.W.; Suzuki, H. Pathway-specific metabolome analysis with (18)O(2)-labeled Medicago truncatula via a mass spectrometry-based approach. Metabolomics 2018, 14, 71. [Google Scholar] [CrossRef] [Green Version]
  29. Dickinson, E.; Rusilowicz, M.J.; Dickinson, M.; Charlton, A.J.; Bechtold, U.; Mullineaux, P.M.; Wilson, J. Integrating transcriptomic techniques and k-means clustering in metabolomics to identify markers of abiotic and biotic stress in Medicago truncatula. Metabolomics 2018, 14, 126. [Google Scholar] [CrossRef] [Green Version]
  30. Li, J.; Essemine, J.; Shang, C.; Zhang, H.; Zhu, X.; Yu, J.; Chen, G.; Qu, M.; Sun, D. Combined proteomics and metabolism analysis unravels prominent roles of antioxidant system in the prevention of Alfalfa (Medicago sativa L.) against salt stress. Int. J. Mol. Sci. 2020, 21, 909. [Google Scholar] [CrossRef] [Green Version]
  31. Gao, Y.; Long, R.; Kang, J.; Wang, Z.; Zhang, T.; Sun, H.; Li, X.; Yang, Q. Comparative proteomic analysis reveals that antioxidant system and soluble sugar metabolism contribute to salt tolerance in Alfalfa (Medicago sativa L.) leaves. J. Proteome Res. 2019, 18, 191–203. [Google Scholar] [CrossRef]
  32. Kang, Y.; Torres-Jerez, I.; An, Z.; Greve, V.; Huhman, D.; Krom, N.; Cui, Y.; Udvardi, M. Genome-wide association analysis of salinity responsive traits in Medicago truncatula. Plant Cell Environ. 2019, 42, 1513–1531. [Google Scholar] [CrossRef] [Green Version]
  33. Boughalleb, F.; Denden, M.; Tiba, B.B. Photosystem II photochemistry and physiological parameters of three fodder shrubs, Nitraria retusa, Atriplex halimus and Medicago arborea under salt stress. Acta Physiol. Plant. 2009, 31, 463–476. [Google Scholar] [CrossRef]
  34. Boughalleb, F.; Denden, M.; Tiba, B.B. Anatomical changes induced by increasing NaCl salinity in three fodder shrubs, Nitraria retusa, Atriplex halimus and Medicago arborea. Acta Physiol. Plant. 2009, 31, 947–960. [Google Scholar] [CrossRef]
  35. Aissat, A.; Mehdadi, Z.; Leogrande, R.; Stellacci, A.M. Characterization of Medicago arborea L. response to water and salt stress. J. Stress Physiol. Biochem. 2019, 15, 103–113. [Google Scholar]
  36. Bennett, R.; Wallsgrove, R. Secondary metabolites in plant defence mechanisms. New Phytol. 1994, 127, 617–633. [Google Scholar] [CrossRef]
  37. Erb, M.; Kliebenstein, D.J. Plant secondary metabolites as defenses, regulators, and primary metabolites: The blurred functional trichotomy. Plant Physiol. 2020, 184, 39–52. [Google Scholar] [CrossRef]
  38. Zaynab, M.; Fatima, M.; Abbas, S.; Sharif, Y.; Umair, M.; Zafar, M.H.; Bahadar, K. Role of secondary metabolites in plant defense against pathogens. Microb. Pathog. 2018, 124, 198–202. [Google Scholar] [CrossRef] [PubMed]
  39. Šamec, D.; Karalija, E.; Šola, I.; Vujčić Bok, V.; Salopek-Sondi, B. The role of polyphenols in abiotic stress response: The influence of molecular structure. Plants 2021, 10, 118. [Google Scholar] [CrossRef] [PubMed]
  40. Zduńska, K.; Dana, A.; Kolodziejczak, A.; Rotsztejn, H. Antioxidant properties of ferulic acid and its possible application. Ski. Pharmacol. Physiol. 2018, 31, 332–336. [Google Scholar] [CrossRef] [PubMed]
  41. Carelli, M.; Biazzi, E.; Tava, A.; Losini, I.; Abbruscato, P.; Depedro, C.; Scotti, C. Sapogenin content variation in Medicago inter-specific hybrid derivatives highlights some aspects of saponin synthesis and control. New Phytol. 2015, 206, 303–314. [Google Scholar] [CrossRef]
  42. Faizal, A.; Geelen, D. Saponins and their role in biological processes in plants. Phytochem. Rev. 2013, 12, 877–893. [Google Scholar] [CrossRef]
  43. Mazahery-Laghab, H.; Yazdi-Samadi, B.; Bagheri, M.; Bagheri, A.R. Alfalfa (Medicago sativa L.) shoot saponins: Identification and bio-activity by the assessment of aphid feeding. Br. J. Nutr. 2011, 105, 62–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Moses, T.; Papadopoulou, K.K.; Osbourn, A. Metabolic and functional diversity of saponins, biosynthetic intermediates and semi-synthetic derivatives. Crit. Rev. Biochem. Mol. Biol. 2014, 49, 439–462. [Google Scholar] [CrossRef] [Green Version]
  45. Basyuni, M.; Baba, S.; Kinjo, Y.; Putri, L.A.; Hakim, L.; Oku, H. Salt-dependent increase in triterpenoids is reversible upon transfer to fresh water in mangrove plants Kandelia candel and Bruguiera gymnorrhiza. J. Plant. Physiol. 2012, 169, 1903–1908. [Google Scholar] [CrossRef] [PubMed]
  46. Gomez-Pando, L.; Álvarez-Castro, R.; Eguiluz de la Barra, A. Short Communication: Effect of salt stress on Peruvian germplasm of Chenopodium quinoa Willd.: A promising crop. J. Agron. Crop. Sci. 2010, 196, 391–396. [Google Scholar] [CrossRef]
  47. Pulvento, C.; Riccardi, M.; Lavini, A.; Iafelice, G.; Marconi, E.; d’Andria, R. Yield and quality characteristics of quinoa grown in open field under different saline and non-saline irrigation regimes. J. Agron. Crop Sci. 2012, 198, 254–263. [Google Scholar] [CrossRef]
  48. Oku, H.; Baba, S.; Koga, H.; Takara, K.; Iwasaki, H. Lipid composition of mangrove and its relevance to salt tolerance. J. Plant Res. 2003, 116, 37–45. [Google Scholar] [CrossRef] [PubMed]
  49. Li, M.; Guo, R.; Jiao, Y.; Jin, X.; Zhang, H.; Shi, L. Comparison of salt tolerance in soja based on metabolomics of seedling roots. Front. Plant Sci. 2017, 8, 1101. [Google Scholar] [CrossRef] [PubMed]
  50. Widodo; Patterson, J.H.; Newbigin, E.; Tester, M.; Bacic, A.; Roessner, U. Metabolic responses to salt stress of barley (Hordeum vulgare L.) cultivars, Sahara and Clipper, which differ in salinity tolerance. J. Exp. Bot. 2009, 60, 4089–4103. [Google Scholar] [CrossRef] [Green Version]
  51. Wu, D.; Cai, S.; Chen, M.; Ye, L.; Chen, Z.; Zhang, H.; Dai, F.; Wu, F.; Zhang, G. Tissue metabolic responses to salt stress in wild and cultivated barley. PLoS ONE 2013, 8, e55431. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, J.; Yang, D.; Li, M.; Shi, L. Metabolic profiles reveal changes in wild and cultivated soybean seedling leaves under salt stress. PLoS ONE 2016, 11, e0159622. [Google Scholar] [CrossRef] [Green Version]
  53. Ouhibi, C.; Attia, H.; Rebah, F.; Msilini, N.; Chebbi, M.; Aarrouf, J.; Urban, L.; Lachaal, M. Salt stress mitigation by seed priming with UV-C in lettuce plants: Growth, antioxidant activity and phenolic compounds. Plant Physiol. Biochem. 2014, 83, 126–133. [Google Scholar] [CrossRef]
  54. Placines, C.; Castañeda-Loaiza, V.; João Rodrigues, M.; Pereira, C.G.; Stefanucci, A.; Mollica, A.; Zengin, G.; Llorent-Martínez, E.J.; Castilho, P.C.; Custódio, A.L. Phenolic profile, toxicity, enzyme inhibition, in silico studies, and antioxidant properties of Cakile maritima Scop. (Brassicaceae) from Southern Portugal. Plants 2020, 9, 142. [Google Scholar] [CrossRef] [Green Version]
  55. Rahman, A.; Hossain, M.S.; Mahmud, J.A.; Nahar, K.; Hasanuzzaman, M.; Fujita, M. Manganese-induced salt stress tolerance in rice seedlings: Regulation of ion homeostasis, antioxidant defense and glyoxalase systems. Physiol. Mol. Biol. Plants 2016, 22, 291–306. [Google Scholar] [CrossRef] [PubMed]
  56. Gupta, P.; De, B. Differential responses of cell wall bound phenolic compounds in sensitive and tolerant varieties of rice in response to salinity. Plant Signal. Behav. 2017, 12, e1379643. [Google Scholar] [CrossRef]
  57. Sharma, A.; Shahzad, B.; Rehman, A.; Bhardwaj, R.; Landi, M.; Zheng, B. Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress. Molecules 2019, 24, 2452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Yamauchi, S.; Taniguchi, E. Synthesis and insecticidal activity of lignan analogs (II). Biosci. Biotechnol. Biochem. 1992, 56, 412–417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Abdel-Farid, I.B.; Marghany, M.R.; Rowezek, M.M.; Sheded, M.G. Effect of salinity stress on growth and metabolomic profiling of Cucumis sativus and Solanum lycopersicum. Plants 2020, 9, 1626. [Google Scholar] [CrossRef] [PubMed]
  60. Ashraf, M.; Ali, Q. Response of two genetically diverse wheat cultivars to salt stress at different growth stages: Leaf lipid peroxidation and phenolic contents. Pak. J. Bot. 2010, 42, 559–565. [Google Scholar]
  61. Haghighi, Z.; Karimi, N.; Modarresi, M.; Mollaei, S. Enhancement of compatible solute and secondary metabolites production in Plantago ovata Forsk. by salinity stress. J. Med. Plants Res. 2012, 6, 3495–3500. [Google Scholar]
  62. Mohamad, A.H.; Marina, M.F.; Francisco José, R.S.; Oscar, V.; Monica, B. Effects of salt and water stress on plant growth and on accumulation of osmolytes and antioxidant compounds in cherry tomato. Not. Bot. Horti Agrobot. Cluj-Napoca 2015, 43. [Google Scholar] [CrossRef]
  63. Irwin, J.; Sewell, J.; Woodfield, D.; Bingham, E. Restructuring Lucerne (Medicago sativa) through introgression of the Medicago arborea genome. Agric. Sci. 2015, 28, 40–46. [Google Scholar]
  64. Hamrouni, L.; Abdallah, F.B.; Abdelly, C.; Ghorbel, A. In vitro culture: A simple and efficient way for salt-tolerant grapevine genotype selection. C. R. Biol. 2008, 331, 152–163. [Google Scholar] [CrossRef] [PubMed]
Ijms 22 04882 g001 550
Figure 1. The stem height and the stem elongation rate (SER) of M. arborea (M. ar), M. sativa (M. sa), and Alborea (Al) under control, 50-75-100 mM and 100 mM. The vertical bars indicate the mean ± standard error (SE). * The different small letters indicate significant differences for the same entry among treatments; different capital letters indicate significant differences for the same treatment among the entries at p < 0.05 (Tukey’s test).
Figure 1. The stem height and the stem elongation rate (SER) of M. arborea (M. ar), M. sativa (M. sa), and Alborea (Al) under control, 50-75-100 mM and 100 mM. The vertical bars indicate the mean ± standard error (SE). * The different small letters indicate significant differences for the same entry among treatments; different capital letters indicate significant differences for the same treatment among the entries at p < 0.05 (Tukey’s test).
Ijms 22 04882 g001
Ijms 22 04882 g002 550
Figure 2. The salinity sensitivity index of M. arborea (M. ar), M. sativa (M. sa), and Alborea (Al) under 50-75-100 mM and 100 mM salt treatments. The vertical bars indicate the mean ± SE. The different letters refer to the significant differences at p < 0.05 (Tukey’s test). * The different small letters refer to the significant differences for the same entry among treatments; different capital letters indicate significant differences for the same treatment among the entries at p < 0.05 (Tukey’s test).
Figure 2. The salinity sensitivity index of M. arborea (M. ar), M. sativa (M. sa), and Alborea (Al) under 50-75-100 mM and 100 mM salt treatments. The vertical bars indicate the mean ± SE. The different letters refer to the significant differences at p < 0.05 (Tukey’s test). * The different small letters refer to the significant differences for the same entry among treatments; different capital letters indicate significant differences for the same treatment among the entries at p < 0.05 (Tukey’s test).
Ijms 22 04882 g002
Ijms 22 04882 g003 550
Figure 3. Principal component analysis (PCA) score plots of secondary metabolic profiles in Medicago entries: (R1 & S1) PCA score plots of the average for all the treatments, i.e., control, 100 mM, 50-75-100 mM, in roots and shoots, (R2 & S2) PCA score plot for, 100 mM in roots and shoots, (R3 & S3) PCA score plot for 50-75-100 mM in roots and shoots. M. sativa (light blue), M. arborea (yellow) and Alborea (dark blue).
Figure 3. Principal component analysis (PCA) score plots of secondary metabolic profiles in Medicago entries: (R1 & S1) PCA score plots of the average for all the treatments, i.e., control, 100 mM, 50-75-100 mM, in roots and shoots, (R2 & S2) PCA score plot for, 100 mM in roots and shoots, (R3 & S3) PCA score plot for 50-75-100 mM in roots and shoots. M. sativa (light blue), M. arborea (yellow) and Alborea (dark blue).
Ijms 22 04882 g003
Ijms 22 04882 g004 550
Figure 4. Heat map analysis of UHPLC-HRMS/MS metabolomics in roots. The row represents each metabolite and the column represents each treatment. Metabolites significantly decreased are displayed in red, while metabolites significantly increased are displayed in blue.
Figure 4. Heat map analysis of UHPLC-HRMS/MS metabolomics in roots. The row represents each metabolite and the column represents each treatment. Metabolites significantly decreased are displayed in red, while metabolites significantly increased are displayed in blue.
Ijms 22 04882 g004
Ijms 22 04882 g005 550
Figure 5. Heat map analysis of UHPLC-HRMS/MS metabolomics in shoots. The row represents each metabolite, and the column represent each treatment. Metabolites significantly decreased are displayed in red, while metabolites significantly increased are displayed in blue.
Figure 5. Heat map analysis of UHPLC-HRMS/MS metabolomics in shoots. The row represents each metabolite, and the column represent each treatment. Metabolites significantly decreased are displayed in red, while metabolites significantly increased are displayed in blue.
Ijms 22 04882 g005
Ijms 22 04882 g006 550
Figure 6. Principal component analysis (PCA) score plots of secondary metabolic profiles in roots and shoots for (a) M. arborea (M. ar), (b) M. sativa (M. sa) and (c) Alborea (Al): (a–c). Shoots (yellow) and roots (dark blue).
Figure 6. Principal component analysis (PCA) score plots of secondary metabolic profiles in roots and shoots for (a) M. arborea (M. ar), (b) M. sativa (M. sa) and (c) Alborea (Al): (a–c). Shoots (yellow) and roots (dark blue).
Ijms 22 04882 g006
Ijms 22 04882 g007 550
Figure 7. Relative fold changes of differential secondary metabolites for (A) roots and (B) shoots of M. sativa (Sa), M. arborea (Ar) and Alborea (Al) under salt stress and salt shock treatments. Btet = benzyl tetrahydrofuran, Bphe = Biphenyl, Phe = Phenol, Fla= Flavonoid, Lig = Lignan, Sa = Saponins, Trac = Triterpenic acid, No = No name.
Figure 7. Relative fold changes of differential secondary metabolites for (A) roots and (B) shoots of M. sativa (Sa), M. arborea (Ar) and Alborea (Al) under salt stress and salt shock treatments. Btet = benzyl tetrahydrofuran, Bphe = Biphenyl, Phe = Phenol, Fla= Flavonoid, Lig = Lignan, Sa = Saponins, Trac = Triterpenic acid, No = No name.
Ijms 22 04882 g007
Table
Table 1. Statistical significance of F ratios from the analysis of variance for stem height of seedlings, stem elongation rate (SER) and sensitivity index.
Table 1. Statistical significance of F ratios from the analysis of variance for stem height of seedlings, stem elongation rate (SER) and sensitivity index.
Source of VariationHeightSERSalinity Sensitivity Index
Salt (A)p ˂ 0.05p ˂ 0.05p ˂ 0.05
Species (B)nsnsp ˂ 0.05
Dates (C)p ˂ 0.05p ˂ 0.05p ˂ 0.05
AXB (Interaction)p ˂ 0.05p ˂ 0.05p ˂ 0.05
AXC (Interaction)p ˂ 0.05nsns
BXC (Interaction)nsp ˂ 0.05p ˂ 0.05
AXBXC (Interaction)ns *p ˂ 0.05ns
* ns: not significant at 0.05 level.
Table
Table 2. VIP features (secondary metabolites) with their retention times (RT), molecular weight, molecular formulas, their tentative identification and their statistical measurements that characterize each species subjected to gradual salinity stress.
Table 2. VIP features (secondary metabolites) with their retention times (RT), molecular weight, molecular formulas, their tentative identification and their statistical measurements that characterize each species subjected to gradual salinity stress.
a. For Aerial Parts (Shoots)
RT [min]Molecular WeightFormulaCategoryLog2 Fold ChangeUpregulatedVIP Score
7.78356.0742C15H16O10Caffeic acid derivative8.01Arb1.08
8.28194.0579C10H10O4Ferulic acid6.93Arb1.52
8.31370.0899C16H18O10Ferulic acid derivative7.8Arb1.07
9.39798.1487C33H34O23Flavonoid (apigenin derivative)10.08Arb7.59
9.55828.1593C34H36O24Flavonoid7.61Arb1.08
10.32828.1746C38H36O21Flavonoid9.36Arb2.39
12.121090.5189C52H82O24Saponin4.82Arb3.80
12.891206.5670C57H90O27Saponin6.34Arb2.07
13.441190.5720C57H90O26Saponin4.77Arb2.07
15.051014.5392C51H82O20Triterpenic saponin7.22Arb1.52
15.33912.5081C47H76O17Soyasaponin4.84Arb1.98
16.21910.4925C47H74O17Triterpenic saponin6.15Arb
8.50432.1267C18H24O12Iridoid−6.49Sat, Alb8.68
9.04622.1170C27H26O17Flavonoid−6.98Sat, Alb3.58
9.59974.1958C43H42O26Flavonoid−9.21Alb5.02
11.07798.1640C37H34O20Flavonoid−6.71Sat, Alb5.69
11.191106.5505C53H86O24Triterpenic saponin−7.84Sat1.01
11.43768.1534C36H32O19Flavonoid−6.89Sat1.34
11.631090.5555C53H86O23Triterpenic saponin−8.71Sat, Alb2.47
11.921414.6249C64H102O34Triterpenic saponin−6.62Sat, Alb1.20
11.531118.5510C54H86O24Triterpenic saponin−7.21Alb1.41
11.741546.6690C62H114O43Triterpenic saponin−7.06Sat, Alb1.78
11.751414.6250C57H106O39Triterpenic saponin−5.72Alb1.04
14.98912.5081C47H76O17Triterpenic saponin−7.9Sat, Alb1.98
15.391028.5190C51H80O21Triterpenic saponin−7.32Sat, Alb1.29
15.79488.2197C30H32O6Triterpenic acid−10.57Sat4.04
b. Roots
RT [min]Molecular WeightFormulaCategoryLog2 Fold ChangeUpregulated inVIP Score
12.38664.3829C36H56O11Phytolaccoside type triterpenic saponin8.04Arb2.21
12.70644.3829C36H56O11Phytolaccoside type triterpenic saponin5.61Arb2.22
12.82648.3879C36H56O10Triterpenic saponin6.06Arb1.14
12.831220.5478C57H88O28Triterpenic saponin7.89Arb2.18
12.831088.5044C52H80O24Triterpenic saponin8.4Arb2.63
12.83956.4618C47H72O20Triterpenic saponin7.57Arb1.01
12.97678.3620C36H54O12Saponin (monoglycosilated)6.23Arb2.15
12.97956.4619C47H72O20Triterpenic saponin6.86Arb1.45
13.40620.3930C35H56O9Steroidal Saponin7.95Arb1.51
13.41782.4457C41H66O14Steroidal Saponin7.11Arb1.31
13.78706.3936C38H58O12Saponin6.97Arb1.14
14.321176.5935C57H92O25Saponin4.03Arb1.66
15.331072.5455C53H84O22Saponin4.55Arb2.67
15.84926.4881C47H74O18Triterpenic saponin5.38Arb1.65
10.72504.1271C24H24O12Flavonoid−6.05Alb1.15
11.18516.0907C24H20O13Dibenzofuran−4.7Sat2.10
11.541118.5513C54H86O24Saponin−7.84Sat1.12
14.71942.5189C48H78O18Soyasaponin I−7.38Sat, Alb3.54
15.28502.3295C30H46O6Medicagenic acid (aglycon)−7.59Alb3.29
15.49438.3132C29H42O3Phenol−7.84Sat, Alb1.64
15.49750.3823C39H58O14Saponin−7.59Sat, Alb7.19
15.49456.3240C29H44O4Triterpenic acid−8.15Sat, Alb1.66
15.49686.3662C38H54O11Triterpenic acid derivative−7.67Sat, Alb1.16
15.79386.1727C22H26O6Lignan−11.99Alb (vs. Sat)13.72
16.271026.50357C51H78O21Steroidal saponin−7.81Sat, Alb1.37
Table
Table 3. VIP features (secondary metabolites) with their tetention times (RT), molecular weight, molecular formulas, their tentative identification, and their statistical measurements that characterize each species subjected to acute salinity stress.
Table 3. VIP features (secondary metabolites) with their tetention times (RT), molecular weight, molecular formulas, their tentative identification, and their statistical measurements that characterize each species subjected to acute salinity stress.
a. For Aerial Parts (Shoots)
RT [min]FormulaMolecular WeightCategoryLog2 Fold ChangeUpregulatedVIP Score
7.78C15H16O10356.0741Caffeic acid derivative8.25Arb1.21
8.28C10H10O4194.0579Ferulic acid7.84Arb1.65
8.31C16H18O10370.0898Cinnamic acid derivative7.95Arb1.11
9.40C33H34O23798.1483Flavonoid (e.g., apigenin glycoside)10.15Arb6.63
9.55C34H36O24828.1591Flavonoid (e.g., diosmetin glycoside)8.23Arb1.67
10.46C40H40O23888.1954Flavonoid (flavonol glycoside)10.21Arb1.73
10.68C39H38O22858.1851Flavonoid (flavonol glycoside)8.38Arb1.92
12.90C57H90O271206.5664Medicagenic acid saponin7.52Arb1.74
8.40C33H34O23798.1492Flavonoid (e.g., apigenin glycoside)−7.03Alb, Sat5.15
8.51C18H24O12432.1268Iridoid−7.36Alb, Sat7.11
9.04C27H26O17622.1173Flavonoid (e.g., apigenin glycoside)−8.65Alb, Sat4.92
9.59C43H42O26974.1961Flavonoid (e.g., apigenin glycoside)−7.25Alb, Sat7.31
10.83C21H18O11455.0850Flavonoid (e.g., flavone glycoside)−5.03Sat1.14
11.07C37H34O20798.1641Flavonoid (e.g., apigenin glycoside)−5.53Sat, Alb5.73
11.28C17H14O6314.0790Flavonoid−7.69Alb (vs. Sat)2.21
11.63C53H86O231090.5560Triterpenic saponin−6.85Alb, Sat2.67
11.91C64H102O341414.6255Saponin−7.43Sat2.40
12.13C48H76O21988.4877Triterpenic saponin−6.89Alb, Sat1.24
12.48C69H110O371530.6740Saponin−6.77Sat1.01
15.796C22H26O6386.1727Lignan−11.57Alb (vs. Sat)16.01
b. Roots
RT [min]Molecular WeightFormulaCategoryLog2 Fold ChangeUpregulatedVIP Score
12.031236.5425C57H88O29Zahnic acid saponin8.03Arb1.55
12.081104.4988C52H80O25Zahnic acid saponin7.33Arb1.03
12.151222.5630C57H90O28Zahnic acid saponin7.21Arb1.14
12.831220.5476C57H88O28Medicagenic acid saponin6.81Arb2.13
12.831088.5040C52H80O24Medicagenic acid saponin8.16Arb3.04
12.83956.4616C47H72O20Medicagenic acid saponin5.77Arb2.46
12.97648.3876C36H56O10Medicagenic acid saponin6.11Arb1.58
12.971106.5337C49H86O27Medicagenic acid saponin5.34Arb1.087
12.97648.3876C36H56O10Medicagenic acid saponin6.48Arb1.36
12.97956.4616C47H72O20Medicagenic acid saponin5.77Arb2.15
12.97632.3563C35H52O10Medicagenic acid saponin6.73Arb1.04
12.97678.3618C36H54O12Medicagenic acid saponin5.82Arb2.64
12.97824.4196C42H64O16Medicagenic acid saponin6.74Arb2.22
13.41470.3398C30H46O4Triterpenic acid7.57Arb2.27
13.48 C48H70O22Steroidic saponin−8.27Sat−2.87
15.50750.3825C39H58O14Steroidic saponin−6.48Alb5.46
15.50502.3296C30H46O6Triterpenic acid (e.g., medicagenic)−7.56Sat, Alb−3.06
15.82284.1258C14H20O6Lignan−11.08Alb1.97
15.82386.1727C22H26O6Lignan−11.99Alb11.93
15.82266.1152C14H18O6Simple phenolic acid−10.19Alb1.85
16.50486.3350C30H46O5Triterpenic acid (e.g., quiallic)−5.17Sat1.51
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sarri, E.; Termentzi, A.; Abraham, E.M.; Papadopoulos, G.K.; Baira, E.; Machera, K.; Loukas, V.; Komaitis, F.; Tani, E. Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea). Int. J. Mol. Sci. 2021, 22, 4882. https://doi.org/10.3390/ijms22094882

AMA Style

Sarri E, Termentzi A, Abraham EM, Papadopoulos GK, Baira E, Machera K, Loukas V, Komaitis F, Tani E. Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea). International Journal of Molecular Sciences. 2021; 22(9):4882. https://doi.org/10.3390/ijms22094882

Chicago/Turabian Style

Sarri, Efi, Aikaterini Termentzi, Eleni M. Abraham, George K. Papadopoulos, Eirini Baira, Kyriaki Machera, Vassilis Loukas, Fotios Komaitis, and Eleni Tani. 2021. "Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea)" International Journal of Molecular Sciences 22, no. 9: 4882. https://doi.org/10.3390/ijms22094882

APA Style

Sarri, E., Termentzi, A., Abraham, E. M., Papadopoulos, G. K., Baira, E., Machera, K., Loukas, V., Komaitis, F., & Tani, E. (2021). Salinity Stress Alters the Secondary Metabolic Profile of M. sativa, M. arborea and Their Hybrid (Alborea). International Journal of Molecular Sciences, 22(9), 4882. https://doi.org/10.3390/ijms22094882

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop