Ionomic and Metabolomic Analyses Reveal Different Response Mechanisms to Saline–Alkali Stress Between Suaeda salsa Community and Puccinellia tenuiflora Community

Soil salinization imposes severe stress to plants, inhibits plant growth, and severely limits agricultural productivity and land utilization. The response of a single plant to saline-alkali stress has been well investigated. However, the plant community that usually works as a group to defend against saline–alkali stress was neglected. To determine the functions of plant community, in our current work, Suaeda salsa (S. salsa) community and Puccinellia tenuiflora (P. tenuiflora) community, two communities that are widely distributed in Hulun Buir Grassland in Northeastern China, were selected as research objects. Ionomic and metabolomic were applied to compare the differences between S. salsa community and P. tenuiflora community from the aspects of ion transport and phenolic compound accumulation, respectively. Ionomic studies demonstrated that many macroelements, including potassium (K) and calcium (Ca), were highly accumulated in S. salsa community whereas microelement manganese (Mn) was highly accumulated in P. tenuiflora community. In S. salsa community, transportation of K to aboveground parts of plants helps to maintain high K+ and low Na+ concentrations whereas the accumulation of Ca triggers the salt overly sensitive (SOS)-Na+ system to efflux Na+. In P. tenuiflora community, enrichment of Mn in roots elevates the level of Mn-superoxide dismutase (SOD) and increases the resistance to saline–alkali stress. Metabolomic studies revealed the high levels of C6C1-compounds and C6C3C6-compounds in S. salsa community and also the high levels of C6C3-compounds in P. tenuiflora community. C6C1-compounds function as signaling molecules to defend against stress and may stimulate the accumulation of C6C3C6-compounds. C6C3-compounds contribute to the elimination of free radicals and the maintenance of cell morphology. Collectively, our findings determine the abundance of phenolic compounds and various elements in S. salsa community and P. tenuiflora community in Hulun Buir Grassland and we explored different responses of S. salsa community and P. tenuiflora community to cope with saline–alkali stress. Understanding of plant response strategies from the perspective of community teamwork may provide a feasible and novel way to transform salinization land.


INTRODUCTION
According to the statistics of the United Nations Educational, Scientific and Cultural Organization (UNESCO) and the Food and Agriculture Organization of the United Nations (FAO), there are over 900 million ha of land with varying degrees of salinization worldwide. These salinized lands span more than 100 countries and account for 6.5% of the world's total land area as well as 60% of the world's total arable lands (Xia et al., 2013;Xu et al., 2020). More seriously, the area of salinized land is estimated to expand at a rate of 1.5 million ha per year. It is predicted that up to 50% of the land will be lost by 2050 (Xia et al., 2013;Xu et al., 2020).
Hulun Buir Grassland located in Northeast China is severely threatened by salinization due to man-made destruction and climate impact. In Hulun Buir Grassland, soil pH is higher than 9.8 and is classified as "saline-alkali" soil (Wang et al., 2007;Kobayashi et al., 2012). Saline-alkali soil normally means NaCl, Na 2 CO 3 , NaHCO 3 , Na 2 SO 4 , and NaOH (Yin et al., 2019). Saline-alkali stress leads to more serious damage than salt or alkali alone (Yu et al., 2013), impairs the physiological functions of plants, disturbs metabolic pathways, and adversely affects the growth of plants Ye et al., 2019). The tolerance or defense responses against saline-alkali stress are composed of many cascades, including signal transduction, ion and osmotic homeostasis, antioxidant system, primary and secondary metabolisms, and transcription regulation. A large number of saline-alkali tolerant plants, such as S. salsa, P. tenuiflora, Mesembryanthemum crystallinum, and halophytes such as Tamarix hispida, Leymus chinensis, Populus nigra, Limonium bicolor, and Puccinellia distans have been investigated to explore the underlying tolerance or defense mechanisms (Kore-Eda et al., 2004;Taji et al., 2004;Wang et al., 2004;Ban, 2007;Wang et al., 2014;Zhao et al., 2016;Yin et al., 2019).
The majority of these researches focus on a single plant. However, in the natural ecological environment, plants resist saline-alkali stress via forming plant communities. S. salsa community and P. tenuiflora community were the dominant community types in Hulun Buir Grassland. S. salsa community is composed of its main constituent S. salsa, and also Atriplex patens and Polygonum sibiricum. P. tenuiflora community is composed of its main constituent P. tenuiflor and also Plantago maritima and Carex reptabunda. Investigating a single plant may result in loss of critical information, and that plant community functions as an entire group was particularly important. However, the systematic roles of plant community under saline-alkali stress have not been clearly reported.
Currently, ionomic and metabolomic-based technologies have been broadly employed to clarify the biological responses of various plants to saline-alkali stress, such as Hordeum vulgare, Oryza sativa, Suaeda corniculata, Glycine max, and Zea mays (Heleno et al., 2015;Pang et al., 2016;Guo et al., 2017;Yang et al., 2017). The application of ionomic and metabolomic offers a rational way to reveal changes under saline-alkali stress. It accelerates the understanding of the response mechanisms of plant salinization tolerance.
To overcome these limitations, instead of investigating a single plant, S. salsa community and P. tenuiflora community, two essential communities in Hulun Buir Grassland, were investigated here. Ionomic and metabolomic were used to compare the differences between S. salsa community and P. tenuiflora community. By studying these plant communities in groups, we attempt to build a response model to salinealkali stress and decipher the resistance strategies of plants to saline-alkali stress from the perspective of community teamwork.
Soils were collected along the vertical length of 20 cm in depth and sieved through a 2-mm nylon sieve. Then 0.2 g dry sample of soil was treated with HF-HClO 4 -HNO 3 method (Hseu et al., 2002;Bielicka-Daszkiewicz et al., 2012). The pH of the soil was measured with a glass electrode pH meter (pHM-2000, Rikakikai Co., Japan) in the saturation phase. The contents of carbonates and bicarbonates were measured by titration with HCl. The contents of chlorides were measured by titration with silver nitrate. The amounts of sulfates were measured by colorimetry (UV-140-02, Shimadzu, Japan) .

Determination of Element Accumulation in Plants
Harvested plants were separated into leaves, stems, and roots and dried for 48 h at 60 • C. Plant tissues were digested in concentrated HNO 3 (95%) using a graphite plate (EH45A plus) at 130 • C. The contents of potassium (K), calcium (Ca), sodium (Na), magnesium (Mg), manganese (Mn), copper (Cu), zinc (Zn), boron (B), nickel (Ni), and molybdenum (Mo) were determined by atomic absorption. Bioaccumulation factors (BFs) and transfer factors (TF) were calculated based on following formulas : Mplant represents mass of detected elements in plant, Mmedium represents mass of detected elements in medium, Mleaf represents mass of detected elements in leaves, and Mroot represents mass of detected elements in roots.

Determination of Primary Metabolisms
Samples weighting 60 mg were extracted as previously described (Chen et al., 2017). After extraction, 1 µL prepared sample solution was injected into the Agilent 7890A-5975C gas chromatography-mass spectrometry (GC-MS) system (Agilent Corporation, United States). GC-MS analysis was carried out on a non-polar DB-5 capillary column (30 m × 250 µm I.D., J&W Scientific, Folsom, CA). High purity of helium was used as the carrier gas at a constant flow rate of 1.0 mL/min. The temperatures of injection and ion source were set to 305 and 230 • C, respectively. Electron impact ionization (-70 eV) at full scan mode (m/z 30-600) was used, with an acquisition rate of 20 spectrum/second in the MS setting. QC sample was prepared by mixing aliquots of tissues samples to be a pooled sample. Acquired data were analyzed by ChromaTOF software (v 4.34, LECO, St Joseph, MI). Internal standards and any known pseudo positive peaks were removed from the obtained data set. Data set was normalized using sum intensity of peaks in each sample. Obtained three-dimensional data sets included sample information, retention time, and peak intensities.

Statistical Analysis
GC-MS and LC-qTOF-MS results were imported into the SIMCA-P14 software package (Umetrics, Umeå, Sweden). Orthogonal partial least squares discriminant analysis (OPLS-DA) was performed to visualize differences between groups after mean centering and unit variance scaling. Significantly different metabolites were screened by the multivariate statistical method and Student's t-test (VIP > 1.0 and P-value < 0.05). A 7-round crossvalidation was applied, with 1/7 of the samples being excluded from the mathematical model in each round to avoid overfitting. Significantly different metabolites were subjected to Kyoto Enrichment of Genes and Genomes (KEGG) enrichment. SPSS version 24.0 software (Chicago, IL, United States) was applied to calculate the score of principal component "Q." A log2 transformation was conducted to improve data normality. Min-max normalization was additionally tested. Hierarchical clustering analysis was performed with R-3.2. Heatmaps, histograms, and pathway maps were drawn with R-3.2 language software, GraphPad Prism8, and Visor, respectively.

RESULTS
The Saline-Alkali Degree of the Rhizosphere Community Suaeda salsa community ( Figure 1A) and P. tenuiflora community ( Figure 1B) are pioneer communities for salinealkali land improvement in Hulun Buir Grassland. Measured indicators of the rhizosphere soils around S. salsa and P. tenuiflora community demonstrated that both communities showed high pH values and could be defined as salinization soils according to soil alkalization classification standard (SSC) ( Table 1). The amounts of anions indicated that S. salsa community rhizosphere soil could be identified as the salinealkali soil of chloride whereas P. tenuiflora community could be identified as saline-alkali soil of sulfate.
Moreover, these recorded parameters of rhizosphere soils implied that S. salsa community is facing more severe alkalinity stress than P. tenuiflora community. It was consistent with the observations that the rhizosphere soil around P. tenuiflora seemed to be more fertile and less exposed than the rhizosphere soil around S. salsa (Figure 1). S. salsa community even presented "leucophylline" and harden into a lump.

The Accumulation of Main Elements in S. salsa Community and P. tenuiflora Community
To explore the influences of elements on plant communities on saline-alkali stress, we tested 11 kinds of main elements in plants and rhizosphere soils, including K, Ca, Na, Mg, Fe, Zn, Cu, Ni, Mn, Mo, and B. Among them K, Ca, Na, Mg, and Mn, were discovered to be significantly different between S. salsa community and P. tenuiflora community by the analysis of OPLS-DA (Supplementary Figure 1A and Supplementary Table 1). Calculation of Q-values of these significantly different elements showed that Q-values of K, Na, Mg, and Ca were obviously higher in S. salsa community as compared with in P. tenuiflora community (Figure 2). The only microelement that was significantly accumulated in P. tenuiflora compared with S. salsa community was Mn (Figure 2). These observations implied that S. salsa community has stronger capabilities of element adsorption, mobilization, and utilization under salinealkali stress.
In addition, the contents of K, Ca, Na, Mg, and Mn in the organs of different kinds of plants in S. salsa community and P. tenuiflora community were measured. In S. salsa community,   Na was mainly accumulated in the leaves of S. salsa and P. sibiricum ( Figure 3A). In P. tenuiflora community, Na was enriched in the stems of P. maritima and the roots of C. reptabunda ( Figure 3B). Calculated BFs showed that community plants obtained abundant Na from soils (BFs > 1) (Figures 3C,D). Plants in S. salsa community had relatively higher TF values as compared with plants in P. tenuiflora community, suggesting that plants in the S. salsa community had stronger Na transport capacities (Figures 3C,D). Measurement of the contents of K demonstrated that K was mainly enriched in aboveground parts of plants of S. salsa community ( Figure 3E). The amount of K was less in P. tenuiflora community and significantly accumulated in the stems of plants in P. tenuiflora community ( Figure 3F). Therefore, the transport capacity of K in P. tenuiflora community might be lower than FIGURE 2 | The Q-values of significantly different elements between S. salsa community and P. tenuiflora community. that in S. salsa community. BFs in all tested S. salsa community plants and C. reptabunda in P. tenuiflora community were higher than 1, indicating that these plants are K hyperaccumulator plants (Figures 3G,H). The high TF values of these plants suggested that S. salsa, A. patens, C. reptabunda, and P. tenuiflora have high K transport capacity from roots to aboveground parts of plants (Figures 3G,H). K/Na ratio, an important index of plants to saline-alkali stress, was further calculated. High K/Na ratios were observed in all plants expect for P. maritime in P. tenuiflora community, demonstrating that these plants ensure the absorption of K under saline-alkali stress ( Figure 3I). Moreover, higher K/Na values in S. salsa community as compared with in P. tenuiflora community implied that the mobilization and transportation of K in S. salsa community are more active to saline-alkali stress whereas plants in P. tenuiflora community are more tolerant to saline-alkali stress.
Magnesium was found to be mainly concentrated in all organs of plants in S. salsa community ( Figure 3J) and the roots of plants in P. tenuiflora community (Figure 3K). High BFs and low TF values of these plants revealed that S. salsa community and P. tenuiflora community plants can absorb a large number of Mg from soil and store Mg in roots (Figures 3L,M). Ca seemed to be equally distributed in leaves, stems, and roots of plants in S. salsa community ( Figure 3N) as well as the leaves and roots in P. maritima in P. tenuiflora community (Figure 3O). BFs and TF values revealed that Ca mainly plays biological roles in leaves and roots in community plants (Figures 3P,Q).
Different from the low abundance of Mn in S. salsa community (Figure 3R), Mn was discovered to be accumulated in the plants in P. tenuiflora community, especially roots ( Figure 3S). The patterns of BFs and TF values were also different in S. salsa community and P. tenuiflora community Frontiers in Plant Science | www.frontiersin.org (Figures 3T,U). High levels of TF values in S. salsa community indicated that Mn is mainly transported to aboveground parts of plants in S. salsa community (Figure 3T). On the other hand, low levels of TF values in P. tenuiflora community suggested that the absorbance and transportation of Mn to aboveground organs is less robust in P. tenuiflora community (Figure 3U).

The Response of Primary Metabolisms to Saline-Alkali Stress
GC-MS detection of primary metabolisms also revealed the differences between S. salsa community and P. tenuiflora community. Among 226 measured metabolisms, a total of 68 primary metabolites showed significantly different expression patterns between S. salsa community and P. tenuiflora community by analysis of OPLS-DA (Supplementary Figure 1B). These primary metabolisms were classified to seven sugars, seven amino acids, 12 alcohols, six esters, four amines, and 32 acids (including 10 phenolic acids and 22 non-phenolic compounds) (Supplementary Table 2). Q-values of primary metabolisms were calculated based on each class. Q-values of amino acids, alcohols, esters, and acids were significantly higher in S. salsa community ( Figure 4A). Q-values of sugars were similar in S. salsa community and P. tenuiflora community ( Figure 4A). Notably, soluble sugars such as tagatose, D-talose, and trehalose were more accumulated in P. tenuiflora community (Supplementary Table 2). It suggested that soluble sugars play an important role in P. tenutiflora community to respond to saline-alkali stress.

The Response of Phenolic Compounds to Saline-Alkali Stress
Considering the importance of acids, especially phenolic compounds, under saline-alkali stress, phenolic compounds were investigated in detail. A total of 34 important phenolic compounds were detected. A total of 24 phenolic compounds were found to be accumulated in our community plants and a total of 11 significantly different metabolisms were screened by OPLS-DA (Supplementary Figure 1C). Differences of acids between S. salsa community and P. tenuiflora community were mainly due to differences in phenolic compounds ( Figure 4B). These differentially expressed phenolic compounds could be classified to C6C1-compounds, C6C3compounds, and C6C3C6-compounds. C6C1-compounds and C6C3C6-compounds were more accumulated in S. salsa community whereas C6C3-compounds were more accumulated in P. tenuiflora community (Figure 4B).
The distributions of C6C1-compounds in plant organs showed that C6C1-compound protocatechuic acid was mainly gathered in the stems and roots of A. patens and P. sibiricum in S. salsa community ( Figure 5A). In P. tenuiflora community, protocatechuic acid was more accumulated in the roots of P. maritima (Figure 5B). Another C6C1-compound, p-hydroxybenzoic acid, was mainly accumulated in the aboveground of P. sibiricum in S. salsa community and also in the roots of C. reptabunda in P. tenuiflora community (Figures 5C,D). C6C3-compound ferulic acid was mainly concentrated in the aboveground part of our community plants (Figures 5E-G). It was significantly accumulated in P. tenuiflora community, especially the leaves of C. reptabunda and P. tenuiflora (Figure 5F). C6C3C6-compounds were evidently enriched in the aboveground part of A. patens and P. sibiricum in S. salsa community (Figures 5H-M). Catechin and kaempferol were only detected in P. sibiricum in S. salsa community. The amounts of naringenin and pelunidin were also higher in P. sibiricum in the S. salsa community.

DISCUSSION
Soil salinization is a major abiotic environmental factor that limits land utilization efficiency (Abdel and Tran, 2016). Given that Hulun Buir Grassland is influenced by obviously gradient salinization, here, we collected the rhizosphere soil of S. salsa community and P. tenuiflora community in Hulun Buir Grassland. We found that the rhizosphere soil of S. salsa community was chloride saline-alkali soil (Cl − and SO 4 2− ) and the rhizosphere soil of P. tenuiflora community was salinealkali sulfate soil (SO 4 2− and CO 3 2− ). Moreover, the unique features of S. salsa community and P. tenuiflora community were characterized by determining the amounts of elements and primary metabolisms, especially phenolic compounds (Figures 6, 7).

Element Response
Sodium is one of the most dominating toxic ions that induce soil salinization. Extremely accumulated Na can induce ionic imbalance and disrupt normal physiological metabolisms (Lin et al., 2018;He et al., 2020). Our observations showed that Na was significantly accumulated in plants in both S. salsa community and P. tenuiflora community, indicating that these plants can survive from excess Na from soils. Na was found to be transported to the aboveground parts of plants and accumulated in leaves and stems. High K + and low Na + concentrations are critical for maintaining osmotic pressure. The accumulation of Na in S. salsa community and P. tenuiflora community implied that the osmotic balance might be disturbed. Since there exists no special Na transportation channel, a competitive relationship between K (or other ions) and Na may be present during K/Na uptake under saline-alkali stress. Excessive Na + can damage metabolic pathways by replacing K + in key enzymatic reactions, causing the dysregulation of various biological processes . However, our results showed that K was also significantly accumulated in the aboveground parts of community plants, indicating that the accumulation of Na may not affect the absorption of K. The value of K/Na ratio suggested that community plants can grow normally under saline-alkali stress. Different from reducing Na uptake to alleviate injuries to plants, K uptake of community plants are seemingly not affected. Therefore, Na may be necessary for the growth of some halophytes.
Salt overly sensitive pathway is an important ionic stress signaling pathway to maintain Na/K balance by extruding Na + into the apoplast. Ca 2+ is a key signaling component in this pathway. Ca accumulation in tissues can instantly trigger the SOS-Na + system to exclude and diminish the damage caused by Na toxicity . The significant accumulation of Ca in plant leaves and roots may activate the SOS-Na + system. Moreover, accumulated Ca can enhance Ca 2+ binding to Ca 2+ sensor proteins and further support Na + detoxification . Ca not only contributes to the SOS-Na + signaling pathway, but also plays a pivotal role in cellular structure maintenance under saline-alkali stress. The accumulation of Ca helps to increase membrane stability under adverse conditions such as osmotic stress (Zeng et al., 2015). Ca accumulated in roots can function as a Ca storage and provide a Ca gradient for multiple biological processes . Our results implied that Ca may play an important role in plant saline-alkali tolerance by increasing the levels of Ca in roots, rather than by reducing the levels of Na.
It has been reported that saline-alkali stress can lead to decreased plant biomass and reduced photosynthesis Nie et al., 2018). Mg participates in chlorophyll biosynthesis and affects photosynthetic rate. Superfluous Na and high pH value can result in deficiency of Mg, affecting the structure of chloroplast and the activation of chlorophyllase enzyme. Our results showed that Mg was strongly absorbed from soil by community plants. The entry of Mg 2+ into root cells is passively mediated via Mg 2+ permeable channels. Salinealkali stress may activate Mg 2+ -uptaking channels in roots. Chloroplast membrane-localized Mg-chelatase H subunit has been demonstrated as a key putative ABA receptor with essential roles in ABA signaling. Plant roots can perceive and defend against saline-alkali stress via ABA signaling . The accumulation of Mg 2+ in roots of S. salsa community plants may regulate ABA signaling transduction and ABA receptor synthesis.
Manganese is also a key element involved in photosynthesis. Here, we identified the enrichment of Mn in P. tenuiflora community. Mn in leaves protects PSII cluster on photosynthetic electron transport whereas Mn in roots improves the resistance to saline-alkali stress by functioning as a component of Mn-SOD Zhu et al., 2020).
It is worth noting that our comparison of elements in S. salsa community and P. tenuiflora community demonstrated different element strategies between these community plants. In S. salsa community plants, macroelements play significant roles under stress. Specifically, K is transported to the aboveground part to maintain high K + and low Na + concentrations and support normal physiological activity. Ca is used for triggering the SOS-Na + system to efflux Na and improve membrane stability. Mg, as a Mg-chelatase H subunit of ABA receptor, plays a role in perceiving and defending against the saline-alkali stress. In P. tenuiflora community, the high abundance of microelement Mn in roots revealed the potential involvement of Mn-SOD in improving the resistance to saline-alkali stress.

Phenolic Compounds Response
Phenolic compounds play primary roles in the regulation of plant's life cycle and growth as well as the acceptance of the environmental challenge. They are the most widely distributed metabolites that are involved in environment responses (Liu et al., 2016;. Phenolic compounds are referred to as cyclic compounds containing at least one hydroxyl group attached to an aromatic ring and can be divided into three major groups, i.e., C6C1-compounds, C6C3-compounds, and C6C3C6compounds, according to the number and binding position of convertible hydroxyl groups on the aromatic chain. C6C1compounds, such as protocatechuic acid and hydroxybenzoic acid, have C6C1 carbon skeletons derived from cinnamic acid. C6C3-compounds, such as ferulic acid and p-hydroxycinnamic acid, have C6C3 carbon skeleton derived from phenylalanine. C6C3C6-compounds (also named flavonoids), such as naringenin and kaempferol, have C6C3C6 carbon skeletons (Liu et al., 2016;Wu et al., 2018;. Phenolic compounds are able to scavenge reactive oxygen species (ROS) and respond to stressful environmental conditions, such as FIGURE 5 | Relative contents of significantly different phenolic compounds in S. salsa community and P. tenuiflora community: (A,C) relative contents of C6C1-compounds in S. salsa community; (B,D) relative contents of C6C1-compounds in P. tenuiflora community; (E) relative contents of C6C3-compounds in S. salsa community; (F,G) relative contents of C6C3-compounds in P. tenuiflora community (H,I,K,L) relative contents of C6C3C6-compounds in S. salsa community (J,M) relative contents of C6C3C6-compounds in P. tenuiflora community; The accumulation of significantly different elements is summarized from 3 biological replicates and presented as the mean ± standard error. Different letters indicate significant differences among treatments (p < 0.05).
FIGURE 6 | Visualization of pathway maps of significantly different elements and phenolic metabolites in S. salsa community. Blue circle, highly expressed in S. salsa, purple circle, highly expressed in A. patents, red circle, highly expressed in P. sibiricum. The green circle, yellow circle, and black circle represent leaf, stem and root, respectively. drought, temperature extremes, and high ultraviolet radiation (Heleno et al., 2015). Some certain species of plants even develop phenolic compounds to inhibit the growth of other competitors (Heleno et al., 2015;.
Some previous studies have demonstrated that organic acid metabolisms can help many plants, such as Chloris virgata, Suaeda salsa, and Leymus chinensis, to adapt to alkali stress (Yang et al., 2008Lin et al., 2014). Similar to their observations, our results showed that 32 acids were enriched in differentially expressed primary metabolisms in S. salsa community and P. tenuiflora community, accounting for half of the significantly expressed metabolisms. The high levels of organic acids may be key adaptive mechanisms that help to stabilize the intracellular pH, neutralize excess toxic cations, and maintain ionic balance . The secretion of organic acids may make the rhizosphere to be more acidic, helping the rhizosphere to neutralize surrounding alkalinity . Our results demonstrated that phenolic compounds occupied a large proportion in differentially expressed acids between S. salsa community and P. tenuiflora community. C6C1-compounds and C6C3C6-compounds were more accumulated in S. salsa community. C6C1-compounds, usually induced by biotic elicitors, are commonly used as signaling molecules to defend stress. We speculate that in P. sibiricum, activated p-hydroxybenzoic acid functions as a signaling cascade and communicates with C6C3C6compounds. There may also be communications between p-hydroxybenzoic acid and protocatechuic acid. Communication signaling was magnified in the aboveground part of plants via p-hydroxybenzoic acid, and affected C6C3C6-compounds accumulated in the aboveground in P. sibiricum (Figure 7). Therefore, there may exist a continuity from C6C1-compounds to C6C3C6-compounds. C6C3C6-compounds accumulated in the epidermal layers of leaves and stems act as filters and absorb radiation in the UV-B portion of the spectrum (Agati et al., 2011;. We speculate that C6C3C6compounds block the toxicity of saline-alkali via a similar way. Excess Na + or ROS can be chelated in the epidermal layers by C6C3C6-compounds to protect cells. Increased C6C3C6compounds may be due to their strong defense and antioxidant FIGURE 7 | Visualization of pathway maps of significantly different elements and phenolic metabolites in P. tenuiflora community. Blue circle, highly expressed in C. reptabunda, pink circle, highly expressed in P. tenuiflora, orange circle, highly expressed in P. maritime. The green circle, yellow circle, and black circle represent leaf, stem and root, respectively. capacity . We speculate that the distributions of the carbon skeleton of phenolic compounds are "primed" by the mechanisms that involve metabolism reprogramming. It will help S. salsa community plants to cope with saline-alkali stress.
The relative contents of C6C1-compounds in P. tenuiflora community were far below their contents in S. salsa community, indicating that the responses of phenolic compounds in P. tenuiflora community are weaker than that in S. salsa community. Interestingly, the resources of the carbon skeleton flow more to the C6C3 pathway in P. tenuiflora community. The roles of C6C3-cmpounds are different from C6C3C6-compounds. C6C3-cmpounds may only have basic reducibility to eliminate free radicals and maintain cell morphology. Therefore, different response modes of phenolic compounds may be applied in S. salsa community and P. tenuiflora community. In S. salsa community, C6C1-compounds stimulate C6C3C6-compounds to block toxicity whereas in P. tenuiflora community, C6C1-compounds induce C6C3-compounds to tolerate stress.
Collectively, in our current study, by applying the community as a unit instead of focusing on a single plant, We achieved characteristics of plant communities, built a conceptual model of community plant cooperation under saline-alkali stress, and provided a novel way for the investigation of response strategies against saline-alkali stress.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.