Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 23 September 2022
Sec. Plant Abiotic Stress
This article is part of the Research Topic Halophytes: Salt Stress Tolerance Mechanisms and Potential Use View all 14 articles

Transcriptome analysis reveals molecular mechanisms underlying salt tolerance in halophyte Sesuvium portulacastrum

Dan Wang,&#x;Dan Wang1,2Nan Yang,&#x;Nan Yang1,2Chaoyue Zhang,Chaoyue Zhang1,2Weihong He,Weihong He1,2Guiping Ye,Guiping Ye1,2Jianjun Chen
Jianjun Chen3*Xiangying Wei,
Xiangying Wei1,2*
  • 1Institute of Oceanography, College of Geography and Oceanography, Minjiang University, Fuzhou, China
  • 2Fuzhou Institute of Oceanography, Fuzhou, China
  • 3Department of Environmental Horticulture, Mid-Florida Research and Education Center, Institute of Food and Agricultural Sciences, University of Florida, Apopka, FL, United States

Soil salinity is an important environmental problem that seriously affects plant growth and crop productivity. Phytoremediation is a cost-effective solution for reducing soil salinity and potentially converting the soils for crop production. Sesuvium portulacastrum is a typical halophyte which can grow at high salt concentrations. In order to explore the salt tolerance mechanism of S. portulacastrum, rooted cuttings were grown in a hydroponic culture containing ½ Hoagland solution with or without addition of 400 mM Na for 21 days. Root and leaf samples were taken 1 h and 21 days after Na treatment, and RNA-Seq was used to analyze transcript differences in roots and leaves of the Na-treated and control plants. A large number of differentially expressed genes (DEGs) were identified in the roots and leaves of plants grown under salt stress. Several key pathways related to salt tolerance were identified through KEGG analysis. Combined with physiological data and expression analysis, it appeared that cyclic nucleotide gated channels (CNGCs) were implicated in Na uptake and Na+/H+ exchangers (NHXs) were responsible for the extrusion and sequestration of Na, which facilitated a balance between Na+ and K+ in S. portulacastrum under salt stress. Soluble sugar and proline were identified as important osmoprotectant in salt-stressed S. portulacastrum plants. Glutathione metabolism played an important role in scavenging reactive oxygen species. Results from this study show that S. portulacastrum as a halophytic species possesses a suite of mechanisms for accumulating and tolerating a high level of Na; thus, it could be a valuable plant species used for phytoremediation of saline soils.

Introduction

Soil salinization is a soil degradation process and considered as one of the most important global issues (Zelm et al., 2020). Excessive salt, mainly sodium (Na) can affect the decomposition of soil aggregates, deplete soil nutrients, and lead to deteriorate soil structure and soil quality, thereby threatening crop production and the environment (Qadir et al., 2001). Globally, about 831 million hectares of land were affected by salinization (Butcher et al., 2016), of which 23.35% is located in Asia (FAO, 2015). Currently, soil salinization takes up to 1.5 million hectares of farmland per year out of crop production (FAO, 2015). Land restoration under the influence of severe salinization is the key to land resource management and sustainable development. Among the strategies, phytoremediation is regarded as an important strategy for the restoration of saline soils because of its cost-effective and environmentally friendly characteristics (Ashraf et al., 2010; Imadi et al., 2016). Halophytes are naturally salt-tolerant plants and have been used in the reclamation of saline-alkali land based on the characteristics of salt accumulation in their leaves and stems (Ayyappan and Ravindran, 2014; Hasanuzzaman et al., 2014; Hayat et al., 2020). Exploring salt tolerance mechanisms of halophytes could assist in our effort on the use of this group of plants for remediation of saline soils.

Halophytes usually live in intertidal zone or inland saline soil. Unlike glycophytes, halophytes can complete their life cycle when grown in seawater or high-salt soils (Chen et al., 2016). Salt stress generally causes excessive accumulation of Na+ in plants, resulting in osmotic stress and ion toxicity and inhibiting plant growth (Hasegawa et al., 2000; Rana and Mark, 2008). To adapt to the high-salt environments, while sharing the salt tolerance mechanism with glycophytes, halophytes have evolved special features. For example, osmotic protection is a general response to salt that maintains cell osmotic pressure and turgor pressure through osmotic regulation or ion membrane transport (Singh et al., 2015). Primary and secondary metabolites, including proline and soluble sugars, play an osmotic adjustment function under salt stress (Chen and Jiang, 2010). On the other hand, to deal with ionic toxicity and maintain Na+/K+ homeostasis in the cell, plants develop mechanisms by limiting the excessive accumulation of Na+ in the cytoplasm, including the restriction of the entry of Na+ into cells, removal of Na+ from the cells, and storing excess Na+ into the vacuole (Zhao et al., 2020). Some important ion transporters include SOS1, HKT1, and NHXs, are involved in the removal of Na+ from the cytoplasm, the transport of Na+ from root cells to xylem, and ion isolation in vacuoles (Shi et al., 2000; Davenport et al., 2010; Apse et al., 2019). Compared with glycophytes, most halophytes can accumulate more Na+ in the aerial portion under salt stress and can maintain a better Na+/K+ ratio (Orsini et al., 2010; Flowers et al., 2015).

Sesuvium portulacastrum, as a mangrove companion plant, is a typical halophyte in the family Aizoaceae and usually grows on wet sand, such as beaches (Lonard and Judd, 1997; He et al., 2022; Zhang et al., 2022). Sesuvium portulacastrum plants are tolerant to salt, drought, and heavy metal stresses and widely used in phytoremediation projects including salt-alkali soil restoration in coastal areas and coastal sand fixation (Lokhande et al., 2013). Early studies have shown that S. portulacastrum is a facultative halophyte and salt accumulator, suggesting that S. portulacastrum plants can grow well in not only non-saline but also saline soils. Appropriate saline levels could even increase the photosynthetic rate of S. portulacastrum plants and promote their growth (Peng et al., 2018). A high salt concentration (200 mM NaCl) significantly alleviated Cd toxicity symptoms in S. portulacastrum by both limiting Cd uptake and compartmenting Cd inside plant tissues (Mariem et al., 2014). Sesuvium portulacastrum plants are able to sustain their growth by isolating salt ions into the vacuole to maintain the osmotic balance between the vacuole and the cytoplasm (Ramani et al., 2006). A quantitative proteomic analysis identified 96 salt-responsive proteins implicated in salt stress in S. portulacastrum, and these proteins are highly involved in ion binding, proton transport, photosynthesis, and ATP synthesis. Under high salinity conditions, the expression of Na+/H+ antiporter and ATP synthase subunits was activated (Yi et al., 2014). Enzymatic analyses showed that NaCl could induce S. portulacastrum P-ATPase and V-ATPase activities, but V-PPase activity was inhibited (Ding et al., 2022). Several genes in S. portulacastrum have been demonstrated to promote salt tolerance in plants. For example, an aquaporin gene SpAQP1 from S. portulacastrum, increased salt tolerance in transgenic tobacco (Chang et al., 2016). Heterologous overexpression of SpSOS1 and SpAHA1 genes improved salt tolerance in yeast and Arabidopsis (Zhou et al., 2015, 2018). With the increasing recognition of S. portulacastrum as an important species for phytoremediation of saline soils and water body, there is a need to further study of this species in tolerance to saline environments. Understanding its tolerance to salinity at the molecular level would enable us to better utilize this and other halophytic plants for remediation of saline soils.

The objectives of this study were to (1) evaluate morphological responses of S. portulacastrum plants to different concentrations of Na in a hydroponic culture, (2) analyze transcript changes of the plants after exposure to a high level of Na in contrast with those without Na treatment through RNA-Seq, (3) identify important genes involved in the tolerance of high Na concentration, (4) outline the major strategies of the plants in adaptation to saline conditions and maintenance of plant growth. It was anticipated that this effort could provide a global view of the transcript changes in and major routes for S. portulacastrum to adapt to saline environments, and such information could assist in our effort on improving crop tolerance of salt stress.

Materials and methods

Plant growth and salt treatments

The shoots of S. portulacastrum were initially collected from Putian, a coastal city in Fujian Province, China. They were potted in containers filled with a substrate composed of 40% peat, 20% pine bark, 20% perlite, and 20% sand based on volume (Zhonghe Agriculture, Huaian, China) and grown in a greenhouse at the Minjiang University, Fuzhou, Fujian, China. Tip cuttings with four nodes and about 10 leaves were made. After removing leaves from the lowest node, they were rooted in ½ Hoagland solution (Hoagland and Synder, 1933). Uniform rooted cuttings were selected, and four cuttings were planted in each 3-L container filled with 2 L ½ Hoagland solution supplemented with NaCl resulting in Na concentrations at 0, 100, 200, 300, 400, 500, and 600 mM, respectively. The plants were grown in a growth room with a constant temperature of 24°C under a light intensity about 400 μmol/m2/s and a photoperiod of 16 h. The experiment was arranged as a completely randomized design with three replications. Plant growth was monitored every 3 days for 21 days by taking each plant out of the solution culture and carefully blotting with paper towel, and fresh weight and root numbers were recorded (Supplementary Table 1). Plant growth photos were taken on days 0, 12, and 21, and growth curves based on root numbers and fresh weight produced under different concentrations of Na were drawn at the end of experiment.

RNA extraction and sequencing

Based on the results of above experiment, 400 mM Na was selected for the second solution culture experiment to determine transcript changes in S. portulacastrum plants. Rooted cuttings were grown in the same containers filled with ½ Hoagland solution devoid of Na or containing 400 mM Na. The experiment was set as a completely randomized experiment with three replications. Root and leaf samples were taken 1 h and 21 days after the initiation of the experiment. The samples were taken from three containers per treatment, respectively, and frozen immediately in liquid nitrogen. Thus, there were three biological replications for each treatment. After 21 days, entire plants were removed from each container, washed with deionized water three times, blotted with paper towel. Leaves and roots were separated, placed in paper bags, and dried at 105°C for 2 h, and then dried at 80°C for 48 h. Dry weight of each sample was weighed, which were used for analysis of Na+ and K+ contents.

A portion of the frozen leaf and root samples were sent to Novogene Co. (Beijing, China) for RNA extraction and sequencing. RNA was isolated from roots and leaves, respectively, and they were sequenced separately. There were three biological replicates for each sample. The RNA integrity and total concentrations were analyzed using Agilent 2100 BioAnalyzer (Agilent Technologies, Palo Alto, CA, United States). 1 μg of RNA from each sample was taken for library construction. Agilent 2100 BioAnalyzer was used to determine the insert size of the library. The qRT-PCR was used to accurately quantify the library effective concentration to ensure the quality of the library. The prepared libraries were sequenced by the Illumina HiSeq 2000 sequencing system (Novogene Co., Beijing, China) to generate 150 bp paired end readings. To obtain the clean data, reads with adapter, reads containing N and low-quality reads (reads with a Qphred <=20 base number) were removed from the raw data. The sequence data was submitted to the National Center for Biotechnology Information (NCBI) with an accession number of PRJNA848266.

Differential expression analysis

Trinity (v2.5.1; Grabherr et al., 2011) was used for transcriptome assembly. Clean reads were mapped on the transcriptome to obtain the read count of each gene using the RSEM software (Dewey and Li, 2011). The read count of genes was converted to FPKM (fragments per kilobase of transcript per million fragments mapped). Differential expression analysis of the two treatments (0.0 mM vs. 400.0 mM at each sampling time) was performed using DESeq2 R package (v1.20.0) based on padj <0.05 and | log2FoldChange | > 1 (Love et al., 2014). KEGG pathway enrichment analysis of the differential expression genes was conducted using KOBAS software (version 3.0; http://kobas.cbi.pku.edu.cn/). Gene functional annotation was performed using the following five databases: Non-Redundant Protein Sequence Database (Nr), Nucleotide Sequence Database (Nt), Protein family (Pfam), Clusters of Orthologous Groups of proteins (KOG/COG), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO). Furthermore, the top 5,000 differentially expressed genes of FPKM were selected for WGCNA analysis. The WGCNA R package was used for co-expression network analysis (Langfelder and Horvath, 2009).

qRT-PCR analysis

To verify the expression of identified genes, RNA was extracted from roots and leaves of S. portulacastrum grown in 0 or 400 mM Na (A portion of the remaining frozen samples mentioned above) using an Omega plant RNA kit (Omega, Bio-Tek Inc. Norcross, GA, United States). After analysis of the quality of RNA, the first strand of cDNA was synthesized with a FastKing RT Kit (Tiangen, Beijing, China) using 1 μg RNA. For qRT-PCR analysis, the cDNA was diluted fivefold with water. The reaction solution for qRT-PCR included 0.8 μl forward primer and 0.8 μl reverse primer, 1 μl cDNA, 10 μl of 2 × SYBR Premix Ex Taq (Takara, Japan) and 7.4 μl deionized water. The reaction program was composed of 95°C for 3 min, then 95°C for 10 s, 60°C for 15 s, and 72°C for 20 s with a total of 45 cycles. qRT-PCR was carried out on an iCycler iQ5 thermal cycler (Bio-Rad, Hercules, CA, United States). The S. portulacastrum SpGAPDH was used as internal control (Zhou et al., 2015). Three biological replicates were performed for each sample. Gene relative expression levels were calculated using the 2−∆∆Ct method (Livak and Schmittgen, 2002). Primers used for the study are listed in Supplementary Table 2.

Analysis of reduced glutathione, soluble sugar, proline, K, and Na

The contents of reduced glutathione (GSH), soluble sugar, and proline in leaves and roots of S. portulacastrum were measured according to the instructions of the commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Fresh plant samples were ground with liquid nitrogen. Precisely weighted each sample was added to PBS buffer with a ratio of 1:9 (g:ml) and vortexed. The solution was centrifuged for 10 min (2,500 rpm), and resultant supernatant was used for testing. Reagents were added according to the manufacturer instructions. The contents of GSH, soluble sugar, and proline were analyzed by a microplate reader with the absorbance at 405, 620, and 520 nm, respectively.

For analysis of Na+ and K+, the aforementioned dried leaf and root samples harvested at the end of the experiment were digested using the method described by Haynes (1980). The inductively coupled plasma (ICP)-atomic emission spectrometry (Perkin Elmer Instruments, Shelton, CT, United States) was used to measure Na+ and K+ content (Plank, 1992). Three biological replicates were analyzed per treatment.

Results

Plant growth responses to different concentrations of Na

Rooted cuttings of S. portulacastrum were able to sustain their growth over a 21-day period regardless of Na concentrations (Figure 1A). Fresh weights of plants exposed to Na at a concentration ranging from 100 to 500 mM were either comparable to or higher than those of the control plants (Figure 1B). The highest fresh weights were produced when plants were exposed to 100 mM Na. However, when exposed to 600 mM Na, the magnitude of fresh weight increase over the time was significantly lower than the other treatments, indicating growth suppression. Root numbers of plants grown at 0 and 100 mM Na were the highest, followed by plants grown with Na at 200–400 mM (Figure 1C). Root numbers of plants grown at 500 mM Na were significantly lower than those grown in the other treatments except for plants grown at 600 mM that had the lowest number of roots.

FIGURE 1
www.frontiersin.org

Figure 1. The growth of Sesuvium portulacastrum in a hydroponic culture. The phenotype of S. portulacastrum plants grown in different concentrations of NaCl over 0, 12, and 21 days (A), fresh weight (B), and root numbers (C) of plants over 21 days of growth. Data were means of three replications. Different letters at the end of lines indicate significant differences among treatments based on Tukey HSD test at p < 0.05 level.

Global transcriptome analysis of S. portulacastrum under salt stress

Since the fresh weights and root numbers were similar between plants grown at 0 mM and 400 mM Na (Figure 1), leaf and root samples of plants grown at the two levels of Na were analyzed using RNA-Seq. We obtained 152.28 Gb of filtered sequence data, of which 90.76% of the data quality reached Q30 or higher (Supplementary Table 3). Because the genome of S. portulacastrum has not been released yet, the reads were mapped to the transcriptome assembled with Trinity (v2.5.1). The BUSCO software was used to evaluate the quality of “cluster.fasta” file obtained by assembling. The comparison rate of complete transcripts (single copy and duplicated copies) reached 55.2% (Supplementary Table 4). All the genes were clustered according to the correlation coefficient (r2). Genes from roots and leaves were separately clustered, and three biological replicates of the same treatment were clustered together (Figure 2A). The results showed that replicate samples were highly correlated. In order to identify differentially expressed salt-responsive genes in roots and leaves, respectively, we compared the expression levels between the following pairs: R_S_1h vs. R_M_1h (R_S_1h represented roots treated with Na for 1 h, and R_M_1h was roots without exposure to Na for 1 h); R_S_21d vs. R_M_21d (R_S_21d, roots treated with Na for 21 days, and R_M_21d, roots without Na treatment for 21 days); L_S_1h vs. L_M_1h (L_S_1h, leaves of plants treated with Na for 1 h, and L_M_1h, leaves of plants without Na treatment); and L_S_21d vs. L_M_21d (L_S_21d, leaves of plants treated with Na for 21 days, and L_M_21d, leaves of plants without Na treatment for 21 days). Compared with the control, a total of 13,013 differentially expressed genes (DEGs) were identified in the roots, of which 8,135 were up-regulated and 4,878 were down-regulated. Whereas 1652 DEGs were identified in leaves, of which 929 were up-regulated, and 723 were down-regulated after 1 h of salt treatment. For plants exposed to 400 mM for 21 days, 18,282 DEGs were identified in the roots (7,584 up-regulated and 10,698 down-regulated) and 3,318 DEGs were identified in leaves (624 were up-regulated, and 2,694 were down-regulated; Figure 2B). These results indicated that the expression of salt-responsive genes varied significantly between roots and leaves, and the expression levels were affected by the duration of the treatments.

FIGURE 2
www.frontiersin.org

Figure 2. RNA-Seq data and DEGs resulted from S. portulacastrum plants treated and untreated with 400 mM NaCl. (A) Hierarchical clustering of 24 samples based on the correlation coefficient (r2), where R_M_1h and R_M_21d were samples of roots (1 h and 21 days) without NaCl treatment; L_M_1h and L_M_21d were samples of leaves (1 h and 21 days) without NaCl treatment; R_S_1h and R_S_21d were samples of roots (1 h and 21 days) with NaCl treatment; and L_S_1h and L_S_21d were samples of leaves (1 h and 21 days) with NaCl treatment. (B) Statistics of DEGs numbers in roots and leaves.

Stress-responsive genes of S. portulacastrum enriched under salt stress

To further explore the changes in gene expression in S. portulacastrum under salt stress, GO enrichment analysis was performed with DEGs by pair comparisons: R_S_1h vs. R_M_1h; R_S_21d vs. R_M_21d; L_S_1h vs. L_M_1h; and L_S_21d vs. L_M_21d (Figures 3AD). Several GO terms related to stress responses were identified in both roots and leaves, such as response to stimulus, response to stress, and response to oxidative stress. But some GO terms were not identified in leaves after salt treatment for 1-h and 21-day, such as response to reactive oxygen species and cellular response to oxidative stress (Figure 3). More DEGs were identified in roots (Figures 3A,B) than leaves (Figures 3C,D), indicating the numbers of stress-responsive genes in roots were significantly higher than those in leaves.

FIGURE 3
www.frontiersin.org

Figure 3. GO analysis of DEGs in roots and leaves of plants treated and untreated with 400 mM Na. The GO terms were identified based on DEGs of different comparisons: R_S_1h vs. R_M_1h (A); R_S_21d vs. R_M_21d (B); L_S_1h vs. L_M_1h (C); L_S_21d vs. L_M_21d (D). The bars represent the number of DEGs in each GO term.

Genes related to ion transport in S. portulacastrum

A large number of DEGs related to ion transporters were identified (Table 1), which included four Na+/K+ transporters (HKTs), 10 shaker-type K+ channels (AKT1s), two high affinity potassium (HAKs), 12 cyclic nucleotide gated channels (CNGCs), 47 aquaporins, and nine glutamate receptors (GLRs). At the same time, multiple DEGs involved in Na+ extrusion and sequestration were identified, including four Na+/H+ exchangers (NHXs) and 71 H+-ATPases. In addition, multiple cation channel genes were detected, including five K+ efflux antiporters (KEAs), three K+ uptake permeases (KUPs), two pore potassium channels (TPKs), two Ca2+/H+ antiporters, and nine cation/H+ antiporters (Supplementary Table 5). In general, the expression of DEGs involved in Na+ extrusion and sequestration increased after being exposed to Na stress.

TABLE 1
www.frontiersin.org

Table 1. Differentially expressed ion transporters in roots and leaves under salt stress.

The analysis of Na+ and K+ contents in roots and leaves showed that the Na+ content increased significantly in roots and leaves after salt treatment. Na+ contents in roots and leaves reached 58.11 and 152.05 mg/g, respectively, after 21 days of salt treatment. Compared with the control samples at the same period, the Na content in the roots and leaves increased by 42.9 folds and 2.1 folds, respectively (Figure 4A). The K+ contents in S. portulacastrum decreased after 21 days of salt treatment, and the contents in roots and leaves were 17.60 and 27.37 mg/g, respectively. Compared with the control samples at the same period, the K content in the roots and leaves was reduced by 0.5 and 0.43 folds, respectively (Figure 4B). We further used qRT-PCR to analyze the expression levels of ion transport related genes in roots. The expression level of CNGC increased significantly after 1 h of salt stress, and decreased after 21-day. The expression level of NHX was significantly increased relative to the control after salt treatment for 1-h and 21-day (Figure 4C). In addition, the expression of K+ transport related genes AKT and HTK1, was inhibited under salt stress (Figure 4C). The quantitative expression trends of these genes were consistent with the trends in RNA-Seq (Figure 4D).

FIGURE 4
www.frontiersin.org

Figure 4. Tissue concentrations of Na and K and qRT-PCR analysis of relevant gene expression. (A,B) The relative content of Na and K in roots and leaves of plants after 1 h and 21 days of growth with 400 mM Na compared to the R_M sample. R_M and L_M mean roots and leaves of plants unexposed to Na. R_S and L_S mean roots and leaves of plants exposed to 400 mM Na. (C) qRT-PCR analysis of the expression of Na+ and K+ transport related genes. The error bars were calculated from three biological replicates. Mean differences were analyzed based on Tukey’s HSD test at p < 0.01 (**) levels. (D) The expression of genes (CNGC, NHX, AKT, and HKT1) in roots of plants treated with or without 400 mM Na. The data represent log2 FPKM values based on RNA-Seq analysis.

Metabolic pathways related to salt stress in S. portulacastrum

In order to explore metabolic pathways involved in salt stress, a pathway enrich analysis was conducted using KEGG. We analyzed DEGs in roots of S. portulacastrum under salt stress and found that 20 pathways were significantly enriched, including photosynthesis – antenna proteins, alpha-linolenic acid metabolism, lysine biosynthesis, and GSH metabolism (Figures 5A,B). Among them, GSH metabolism was among the most enriched pathways. As a result, the expression levels of genes involved in GSH metabolism were analyzed. Results showed that after 1-h salt treatment, a large number of glutathione S-transferase (GST) genes were highly up-regulated in roots. In addition, genes involved in catalyzing the production of oxidized glutathione (GSSG) from GSH, such as isocitrate dehydrogenase (IDH), 6-phosphogluconate dehydrogenase (PGD), glucose-6-phosphate 1-dehydrogenase (G6PD), and glutathione peroxidase (GPX) were mostly up-regulated (Figure 5C). However, after 21-day salt treatment, the expression of these DEGs in roots were largely down-regulated compared to those after 1-h salt treatment (Figure 5C). We further analyzed the GSH content and found that salt treatment led to a reduction of the GSH content in the S. portulacastrum roots (Figure 5D).

FIGURE 5
www.frontiersin.org

Figure 5. KEGG analysis of DEGs in roots and expression analysis of DEGs involved in “glutathione (GSH) metabolism.” (A) The top 20 pathways that were significantly enriched in roots after 1 h of salt treatment (400 mM Na). (B) The top 20 pathways that were significantly enriched in roots after 21 days of salt treatment. (C) Expression analysis of DEGs in roots involved in “Glutathione metabolism” after 1 h and 21 days of salt treatment. (D) The relative content of GSH in roots and leaves compared to the R_M sample. R_M and L_M mean roots and leaves of plants untreated with Na. R_S and L_S mean roots and leaves of plants treated with 400 mM Na. The error bars were calculated from three biological replicates. Mean separation was based on Tukey HSD test at p < 0.01 (**) levels.

The KEGG pathway enrichment was also used for analysis of DEGs in leaves, and 20 significantly enriched pathways were identified, including plant hormone signal transduction, phenylpropanoid biosynthesis, starch and sucrose metabolism, and GSH metabolism (Figures 6A,B). Considering the importance of plant hormones in abiotic stress, we analyzed the expression changes of DEGs involved in the plant hormone signal transduction. After 1-h salt treatment, DEGs related to signaling network of auxin: auxin influx carrier (AUX), auxin-responsive protein (IAA), gibberellin (gibberellin receptor, GID1), abscisic acid (protein phosphatase 2C, PP2C), and ABA responsive element binding factor (ABF) were highly up-regulated. While DEGs involved in ethylene (ethylene-responsive transcription factor 1, ERF1), brassinosteroid (BR-signaling kinase, BSK), and xyloglucosyl transferase (TCH) were down-regulated. On the other hand, the numbers of DEGs involved in plant hormone signal transduction were lower after 21-day salt treatment than those after 1-h salt treatment (Figure 6C). Soluble sugars play an important role in the osmotic adjustment of plant salt tolerance. Soluble sugar in leaves of S. portulacastrum increased 1.5 fold after 1 h of salt treatment but decreased 0.72 fold after 21 days of salt treatment (Figure 6E). These results were consistent with the expression trend of genes related to starch and sucrose metabolism in the transcriptome (Figure 6D).

FIGURE 6
www.frontiersin.org

Figure 6. KEGG analysis of DEGs in leaves and expression analysis of DEGs involved in “plant hormone signal transduction” and “starch and sucrose metabolism.” (A) The top 20 pathways that were significantly enriched in leaves of plants after 1 h of salt treatment (400 mM Na). (B) The top 20 pathways that were significantly enriched in leaves after 21 days of salt treatment. (C,D) Expression analysis of DEGs in leaves involved in “plant hormone signal transduction” and “starch and sucrose metabolism” after 1 h and 21 days of salt treatment. (E) The relative content of soluble sugar content in roots and leaves compared to the R_M sample. R_M and L_M mean roots and leaves of plants untreated with Na. R_S and L_S mean roots and leaves of plants treated with 400 mM Na. The error bars were calculated from three biological replicates. Mean separation was based on Tukey HSD test at p < 0.01 (**) levels.

Our results showed that genes involved in proline metabolism were differentially expressed (Figures 7A,B) and proline content in S. portulacastrum significantly increased after 21-day salt stress (Figure 7C). Proline is an important osmolyte in higher plants, which is synthesized from glutamic acid and ornithine. ∆-1pyrroline-5-carboxylate synthetase (P5CS) and 1-pyrroline-5-carboxylate reductase (P5CR) catalyze the biosynthesis of proline from glutamic acid (Székely et al., 2008). In catabolism, proline dehydrogenase (PDH) degrades proline as pyrroline-5-carboxylate (P5C; Peng et al., 1996). In roots, the increased content of proline was mainly due to the up-regulated expression of P5CS, P5CR and ornithine-oxo-acid transaminase (rocD), which catalyzed the production of proline from ornithine. At the same time, the down-regulated expression of PDH might inhibit the degradation of proline (Figure 7A). The increase in proline content in leaves could be due to the inhibition of PDH expression (Figure 7B). Additionally, genes in regulation of polyamine metabolism were also differentially expressed (Figures 7D,E).

FIGURE 7
www.frontiersin.org

Figure 7. Expression analysis of DEGs involved in proline metabolism. (A) Expression analysis of DEGs in S. portulacastrum roots involved in proline metabolism after 1 h (R_S_1h vs. R_M_1h) and 21 days (R_S_21d vs. R_M_21d) of salt treatment (400 mM Na). (B) Expression analysis of DEGs in S. portulacastrum leaves involved in proline metabolism after 1 h (L_S_1h vs. L_M_1h) and 21 days (L_S_21d vs. L_M_21d) of salt treatment. (C) The relative content of proline in roots and leaves compared to the R_M sample. R_M and L_M mean roots and leaves under control condition. R_S and L_S mean roots and leaves under salt stress. The error bars were calculated from three biological replicates. Mean separation was based on Tukey HSD test at p < 0.01 (**) levels. (D,E) Expression analysis of DEGs in S. portulacastrum roots and leaves respectively involved in polyamines metabolism after 1 h and 21 days of salt treatment.

Differential expression of transcription factors under salt stress

A large number of transcription factors (TFs) were differentially expressed after the exposure to Na. In roots of S. portulacastrum, 367 TFs and 452 TFs were differentially expressed 1-h and 21-day after salt treatment, respectively (Supplementary Table 6). The top five most differentially expressed transcription factors were genes in MYB, C2H2, AP2/ERF, WRKY, and bZIP families (Figure 8A). Compared with roots, the number of differentially expressed TFs in S. portulacastrum leaves was significantly lower than in roots. There were 367 and 452 DEGs in roots compared to 77 and 111 in leaves after 1 h and 21 days, respectively. The top five abundant TF families in leaves were AP2/ERF, MYB, WRKY, NAC, and bHLH (Figure 8A). The analysis of differentially expressed TFs showed that more TFs in the roots than leaves of S. portulacastrum after 1-h salt treatment, and their expressions were mostly up-regulated (Figure 8B). However, after 21 days of salt treatment, the number of down-regulated TFs increased significantly (Figure 8C).

FIGURE 8
www.frontiersin.org

Figure 8. Analysis of differentially expressed TFs in roots and leaves of S. portulacastrum. (A) Venn diagrams for salt-induced TFs in roots and leaves of S. portulacastrum. R_1h and R_21d means differentially expressed TFs in roots after salt treatment (400 mM Na) for 1 h and 21 days. L_1h and L_21d means differentially expressed TFs in leaves after salt treatment after 1 h and 21 days. The TF families and gene numbers (n > 4) were shown in the right. (B) Expression analysis of the TFs in each comparison.

Identification of coexpression modules related to salt tolerance

In order to gain a better understanding of the relationships among different genes implicated in S. portulacastrum responses to salt stress, the WGCNA analysis was carried out. The top 5000 DEGs with the average FPKM were divided into 15 modules, which were presented in the cluster dendrogram (Figure 9A). Each module was marked with a different color. Further correlation analysis between the modules and the samples was carried out (Figure 9B). The results showed that the yellow module, with 431 identified genes, was highly associated with the R_S_1h sample (Supplementary Table 7). The brown module (438 genes) was highly associated with R_S_21d sample. The purple module, representing 144 genes, was highly correlated with L_S_1h sample. The cyan module, containing 45 genes, was highly related to the L_S_21d sample.

FIGURE 9
www.frontiersin.org

Figure 9. WGCNA analysis of genes in roots and leaves of S. portulacastrum. (A) Hierarchical cluster tree showing co-expression modules identified by WGCNA. Each leaf in the tree represents one gene. The major tree branches constitute 15 modules, labeled with different colors. (B) Module–sample association. Each row corresponds to a module, labeled with a color as in (A). The numbers in each cell at the row-column intersection indicates the correlation coefficient between the module and the sample. (C) Heatmaps showing the selected hub genes in each module.

The WGCNA was also be used to screen the hub genes in the module based on the KME (eigengene connectivity) value. The hub gene means a high degree of connection with the genes in the module, which may play an important role in the salt tolerance of S. portulacastrum. Heatmap showed that the yellow module-specific genes were over-expressed in R_S_1h (Figure 9C). Cytochrome P450 71A1, UDP-glycosyltransferase 88B1, cellulose synthase-like protein G2, chitin elicitor receptor kinase 1, and chitin synthase were identified as candidate hub genes in this module. The brown module genes were over-expressed in R_S_21d (Figure 9C). Genes encoding peroxidase 12, cysteine-rich receptor-like protein kinase, clathrin interactor 1, CO (2)-response secreted protease, and CNGC15b were identified as candidate hub genes for this module. The purple module genes were over-expressed in L_S_1h (Figure 9C). Glutaredoxin-C6, fructose-bisphosphate aldolase 1, ultraviolet-B receptor UVR8, heat shock factor protein HSF30, and cytochrome P450 4d2 were identified as candidate hub genes in this module. The cyan module genes were over-expressed in L_S_21d (Figure 9C). Genes encoding bifunctional purple acid phosphatase 26, vesicle transport v-SNARE 12, telomere capping C-terminal wHTH, translation factor SUI1 homolog, and beta-3 tubulin were identified as candidate hub genes for this module.

Discussion

Land salinization is one of the most important environmental problems in the world. It was estimated that about 10% of the land surface and 50% of irrigated land have been affected by salinization (Ruan et al., 2010). Annual losses due to salt-affected land in agricultural production are more than US$12 billion and rising (Flowers et al., 2010). Salinity is now a key restraint to crop productivity (Shahzad et al., 2019). On the other hand, world food production is projected to increase between 50% and 70% by 2050 to match the estimated growth of population (Shabala et al., 2013). It is certain that the current arable land is far less sufficient for crop production. A solution to expand crop production is to use halophytic plants to remediate soil salinity and convert the soils to cultivable land (Karakas et al., 2020; Munir et al., 2022). Thus, an in-depth understanding of halophytic plant tolerance to salt stress would enable us to better use this group of plants for remediation of saline soils. Sesuvium portulacastrum is a typical halophyte, but its salt tolerance mechanisms remain largely unclear. Our study showed that S. portulacastrum plants when exposed to 400 mM Na were able to produce shoot fresh weight and root numbers comparable to those of control plants (Figure 1). The ability of S. portulacastrum to thrive in such a higher concentration of Na was elucidated through RNA-Seq and qRT-PCR analyses. Our results showed that S. portulacastrum has developed multiple strategies to cope with saline growth conditions.

Regulation of cellular Na+/K+ homeostasis

Maintaining a proper intracellular K+/Na+ ratio is critical for plants to tolerate salt stress. Sesuvium portulacastrum is able to deliberately control K+ and Na+ levels when grown under high salt stress. Compared with the control samples, higher content of Na+ was detected in the leaves and roots of S. portulacastrum after 21 days of salt treatment. At the same time, the K+ content in the samples of the treatment group was significantly reduced (Figures 4A,B). Like other halophytes, S. portulacastrum can maintain a higher Na+/K+ ratio in roots and leaves. This is likely due to the expression of ion transport related genes. We have identified a large number of ion transport related genes that were differentially expressed in S. portulacastrum under salt stress. The expression of CNGC genes in roots after 1 h of Na treatment was up-regulated (Table 1). This suggests that CNGC gene was involved in Na uptake in the early stage of salt stress. After 21 days of salt treatment, the expression of CNGC was down-regulated, which could inhibit Na uptake at the later stage of salt stress (Figure 4C). Recent studies have shown that increased transcript levels of a few AtCNGCs in Arabidopsis roots or shoots after exposed to high levels of NaCl (Duszyn et al., 2019). For example, AtCNGC20 was up-regulated in roots after exposure to Na; AtCNGC5 and AtCNGC17 were reported to be implicated in salt tolerance (Guo et al., 2008; Mian et al., 2011; Ladwig et al., 2015; Massange-Sánchez et al., 2016; Vadim, 2018). Our speculation is that ion transporter, such as CNGCs regulate Na absorption, while NHXs may play an important role in control of Na accumulation. From our result, NHX and H+-ATPase genes were up-regulated in S. portulacastrum root under salt stress (1 h and 21 days). NHX is one of the important families participating in Na+ extrusion and sequestration. NHXs located on the plasma membrane of the cell, such as NHX7/SOS1 can exclude Na+ from the root system (Ji et al., 2013). In addition, NHX1 and NHX2 located on the vacuole membrane may participate in the compartmentalization of Na+ in the vacuole and the K+ balance in the vacuole (Jiang et al., 2010). At the same time, proton pump (H+-ATPase) located on the plasma membrane of the cell provided power for Na+ extrusion and sequestration (Gaxiola et al., 2007). In addition, several ion transport related genes in S. portulacastrum had been reported. Previous studies found that heterologous overexpression of SpSOS1 and SpAHA1 genes improved salt tolerance in yeast and Arabidopsis (Zhou et al., 2015, 2018). These genes were also identified in our study (Supplementary Table 5). In addition, more ion transport related genes were identified in our study, such as AKT and HKT. The expression of AKT and HKT in roots of S. portulacastrum was down-regulated after salt stress. But in halophytic turf grass, SvHKT1;1 gene was up-regulated in response to high concentrations of NaCl (500 mM), (Kawakami et al., 2020). This is in contrast to the results in the S. portulacastrum. We think that CNGCs and NHXs are important gene families implicated in the regulation of Na+ uptake, extrusion, and sequestration, and their coordinate actions could maintain Na+/K+ at appropriate ratios and lead to Na tolerance in S. portulacastrum. Further research is warranted to confirm these propositions.

Changes in osmolytes

The loss of water due to the decrease in osmotic pressure is one of the great constraints faced by plants growing in saline soils. Accumulation of compatible osmolytes is an important strategy for plants to cope with osmotic stress (Apse and Blumwald, 2002). Different species have different osmolyte profiles under salt stress. Halophytes usually accumulate one or more than one compatible osmolytes, such as proline, glycine betaine, and sorbitol. Previous studies showed that under salt and drought conditions, the proline content in the callus and axillary buds of S. portulacastrum significantly increased (Lokhande et al., 2010). Our study also found that, compared with the control samples, the proline content in S. portulacastrum roots and leaves significantly increased under short-term (1-h) and long-term (21-day) salt treatments (Figure 7C). According to transcriptome data, we found that the accumulation of proline in roots depends on the up-regulated expression of genes related to proline synthesis (P5CS, P5CR, rocD) on the one hand, and on the other hand from the down-regulated expression of genes related to degradation (PDH; Figure 7A). However, the accumulation of proline in leaves is more dependent on the inhibition of degradation-related genes (Figure 7B). Soluble sugar is another important osmolyte under salt stress, which can protect specific macromolecules or sustain membrane stability (Bartels and Sunkar, 2005). Other studies have demonstrated that halophyte Thellungiella under salt treatment accumulated high content of soluble sugars (Wang et al., 2013). From the result of KEGG analysis, the pathway of starch and sucrose metabolism was enriched in leaves of S. portulacastrum under salt stress (Figures 6A,B). Further study showed that the soluble sugar content in leaves of S. portulacastrum increased after 1 h of salt stress but decreased after 21d of salt stress (Figure 6E). Thus, proline and soluble sugar may play an important role in maintaining the osmotic balance of S. portulacastrum in a high-salt environment.

Polyamines play significant roles in regulating plant defense responses to various environmental stresses, including salt stress, metal toxicity, and oxidative stress (Zhong et al., 2020). Exogenous application of putrescine (Put) improved plant salt tolerance (Chattopadhayay et al., 2002; Gill and Tuteja, 2010). In plants, the first rate-limiting reaction of Put synthesis is either arginine catalyzed by arginine decarboxylase (ADC) or ornithine catalyzed by ornithine decarboxylase (ODC). In this study, after 1 h of salt stress, the expression of several ODC genes in the root of S. portulacastrum increased, which catalyzed the biosynthesis of putrescine from ornithine. The expression of polyamine oxidase 4 (PAO4) was up-regulated (Figure 7D). After 1 h of salt stress, the polyamine oxidase gene MAPO was up-regulated in leaves of S. portulacastrum (Figure 7E). PAO4 is able to catalyze the biosynthesis of putrescine from spermine and spermidine (Alcázar et al., 2011). After 21 days of salt stress, the expression of spermidine synthase genes (speEs) in S. portulacastrum was up-regulated, which catalyzes putrescine to spermine and spermidine (Figure 7D). This suggests that the balance of polyamines plays a role in the salt tolerance of the S. portulacastrum.

ROS homeostasis

Salt stress can induce excessive ROS accumulation in plants which leads to cell damage and cell death. The tolerance of halophytes to high salinity depends to a certain extent on their ability to maintain ROS homeostasis (Jayakumar et al., 2014). Synthesis of antioxidant metabolites is an important strategy for plants to eliminate excess ROS under high salt stress. Ascorbate and GSH are important antioxidant metabolites in plants. Halophytes are found to have a more efficient ascorbate-glutathione cycle. For example, the ascorbate and GSH content of halophyte Lycopersicon pennellii was significantly higher than that of glycophyte relative L. esculentum (Shalata et al., 2001). However, the content of GSH in S. portulacastrum decreased after 1-h of salt treatment (Figure 5D). This may be due to the oxidation of GSH to GSSG. RNA-Seq data also showed that the expression of genes involved in catalyzing the oxidation of GSH into GSSG was significantly up-regulated (Figure 5C). Another important way to scavenge ROS in halophytes is through enzymatic agents, such as SOD (Grene et al., 2002), CAT (Willekens et al., 2014), APX (Shigeru et al., 2002), and POX (Fagerstedt et al., 2010). From the WGCNA analysis of this study, POX was identified as one of the hub genes in the salt-stressed modules in roots (Figure 9C). POXs are isoenzyme that removes H2O2 in the outer plastid space (Fagerstedt et al., 2010). Previous studies showed that increased activity of POXs was found in halophytes, such as Cakile maritima and Hordeum marinum (Seckin et al., 2010; Ellouzi et al., 2011). Our results also indicated that the expression of POX genes was up-regulated in roots of S. portulacastrum under salt stress (Figure 9C). Therefore, we believe that the GSH/GSSH cycle and POX play an important role in the scavenging ROS induced by salt in S. portulacastrum.

Changes in tissue-special gene modules

Another finding in this study is tissue-special gene modules identified by WGCNA analysis (Figure 9). WGCNA analysis aimed to identify co-expressed gene modules and to explore the core genes in the modules, through which a hierarchical clustering tree was constructed by calculating the correlation coefficient between any two genes. Different branches of the clustering tree represent different gene modules, and different colors represent different module. Based on the weighted correlation coefficients of genes, genes were classified according to their expression patterns, and genes with similar patterns were grouped into one module. The identification of hub genes in relation to modules provided a new and important perspective for understanding the salt tolerance mechanism in S. portulacastrum. In roots of S. portulacastrum, the expression of genes related to lignin synthesis in the phenylpropane metabolic pathway (cytochrome P450 71A1-like), cellulose synthesis (cellulose synthase-like protein G2), endocytosis (clathrin interactor 1-like), and perceive CO2 concentration to regulate the pathway of stomatal development CO2-response secreted protease-like) were up-regulated. An important adaptation of plants to salt stress is differential regulation of growth, accompanied by dynamic changes and rearrangement of plant cell walls, and cellulose and lignin were the structural composition of cell walls (Cosgrove, 2015; Tenhaken, 2015; Yan et al., 2021). This may indicate that the strong cell wall regulation ability of S. portulacastrum enhances its salt tolerance. In leaves of S. portulacastrum, the expression level of vesicle transport related genes (vesicle transport v-SNARE 12) and tubulin (beta-3 tubulin) were up-regulated under salt stress. Vesicle transport and tubulin play important roles in plant salt tolerance. Under high salinity conditions, regulating the reorganization of cortical microtubules is the key to the survival of plant cells (Endler et al., 2015). In Arabidopsis, ethylene signal promotes the reorganization of cortical microtubules induced by salt stress (Dou et al., 2018). This inspired us to think what role does tubulin play in promoting salt tolerance in halophytes? These questions will be part of our follow-up research.

The mechanisms of salt tolerance in halophytes have been continuously studied over the years. Most studies focus on osmotic protection and ion homeostasis. Some halophytes maintain Na+/K+ balance through salt secretion. For example, in halophyte quinoa, salts are deposited in epidermal bladder cells (EBCs; Kiani-Pouya et al., 2017). In S. portulacastrum, the maintenance of Na+/K+ balance is mainly dependent on the storage of the excess Na in the vacuole. Halophytes maintain a better Na/K ratio by removing Na from the cytoplasm, transporting Na from root cells to the xylem, and sequestrating Na in the vacuole in which SOS1, HKT1, and NHXs ion transporters play important roles (Flowers et al., 2015). In this study, more ion transport related genes (GLRs, CNGCs, AKTs, HAKs, KEAs, KUPs and TPKs) were identified in response to salt stress. In addition, the results of WGCNA analysis suggest that cell wall formation, vesicle transport and tubulin may play significant roles in salt tolerance in S. portulacastrum.

Conclusion

Sesuvium portulacastrum is a halophytic species and has been considered a valuable plant to be used for phytoremediation of saline soils. After 400 mM Na treatment, Na contents in roots and leaves were much higher than K contents in S. portulacastrum. A large number of DEGs were detected from RNA-Seq data. A variety of ion transport related genes (CNGCs, NHXs, AKTs, HKTs) are implicated in Na uptake and accumulation, of which CNGCs could be involved in Na uptake, and NHXs could be responsible for controlling Na extrusion and sequestration. Their coordinated action could maintain Na+/ K+ homeostasis in the S. portulacastrum. Soluble sugar and proline play an important role in the osmotic regulation of S. portulacastrum under salt stress. Glutathione metabolism and POX participated in scavenging reactive oxygen species under salt stress. Additionally, transcription factors and plant hormones including auxin, gibberellin, abscisic acid, and ethylene are involved in regulating the salt tolerance of S. portulacastrum. Our study shows that S. portulacastrum has developed a suite of mechanisms for accumulation and tolerance of Na, suggesting that it could be an important species to be used for remediation of saline soils (Figure 10).

FIGURE 10
www.frontiersin.org

Figure 10. Schematic presentation of the potential mechanism of salt stress tolerance in S. portulacastrum. The yellow arrows in the figure represent the increased expression of related genes. Blue arrows represent decreased gene expression. The mentioned genes were selected from plants after 1 h of Na treatment.

Data availability statement

The datasets presented in this study were deposited in the NCBI repository with an accession number of PRJNA848266, available at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA848266.

Author contributions

DW: conceptualization, data curation, formal analysis, investigation, methodology, and writing–original draft. NY: data curation, methodology, software, and visualization. CZ: methodology. WH: methodology. GY: methodology and software. JC: conceptualization and writing–review and editing. XW: conceptualization, funding acquisition, project administration, resources, supervision, and writing–review and editing. All authors contributed to the article and approved the submitted version.

Funding

The authors are thankful to the National Natural Science Foundation of China (32102331), the Natural Science Foundation of Fujian Province (2022J011141 and 2020J01867) and Department of Education, Fujian Province (JAT200426) for providing the financial support.

Acknowledgments

We are grateful to all the members of our team for helpful discussion on this project.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.973419/full#supplementary-material

References

Alcázar, R., Bitrián, M., Bartels, D., Koncz, C., Altabella, T., and Tiburcio, A. F. (2011). Polyamine metabolic canalization in response to drought stress in Arabidopsis and the resurrection plant Craterostigma plantagineum. Plant Signal. Behav. 6, 243–250. doi: 10.4161/psb.6.2.14317

PubMed Abstract | CrossRef Full Text | Google Scholar

Apse, M. P., and Blumwald, E. (2002). Engineering salt tolerance. Cur. Opin. Biotech. 13, 146–150. doi: 10.1016/S0958-1669(02)00298-7

CrossRef Full Text | Google Scholar

Apse, M. P., Aharon, G. S., Snedden, W. A., and Blumwald, E. (2019). Salt tolerance conferred by overexpression of a vacuolar Na+/H+ antiport in Arabidopsis. Science 285, 1256–1258. doi: 10.1126/science.285.5431.1256

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashraf, M. Y., Ashraf, M., Mahmood, K., Akhter, J., Hussain, F., and Arshad, M. (2010). “Phytoremediation of saline soils for sustainable agricultural productivity,” in Plant Adaptation and Phytoremediation. eds. M. Ashraf, M. Ozturk, and M. Ahmad (Dordrecht: Springer Netherlands), 335–355.

Google Scholar

Ayyappan, D., and Ravindran, K. C. (2014). Potentiality of suaeda monoica forsk. A salt marsh halophyte on bioaccumulation of heavy metals from tannery effluent. Int. J. Modn. Res. Revs. 2, 267–274. doi: 10.5829/idosi.wasj.2013.28.12.2136

CrossRef Full Text | Google Scholar

Bartels, D., and Sunkar, R. (2005). Drought and Salt Tolerance in Plants. Crit. Rev. Plant Sci. 24, 23–58. doi: 10.1080/07352680590910410

CrossRef Full Text | Google Scholar

Butcher, K., Wick, A. F., DeSutter, T., Chatterjee, A., and Harmon, J. (2016). Soil salinity: a threat to global food security. Agron. J. 108, 2189–2200.

Google Scholar

Chang, W., Liu, X., Zhu, J., Fan, W., and Zhang, Z. (2016). An aquaporin gene from halophyte Sesuvium portulacastrum, SpAQP1, increases salt tolerance in transgenic tobacco. Plant Cell Rep. 35, 385–395. doi: 10.1007/s00299-015-1891-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chattopadhayay, M. K., Tiwari, B. S., Chattopadhyay, G., Bose, A., Sengupta, D. N., and Ghosh, B. (2002). Protective role of exogenous polyamines on salinity-stressed rice (Oryza sativa) plants. Physiol. Plant. 116, 192–199. doi: 10.1034/j.1399-3054.2002.1160208.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., and Jiang, J. G. (2010). Osmotic adjustment and plant adaptation to environmental changes related to drought and salinity. Environ. Rev. 18, 309–319. doi: 10.1139/A10-014

CrossRef Full Text | Google Scholar

Chen, T., Yuan, F., Song, J., and Wang, B. (2016). Nitric oxide participates in waterlogging tolerance through enhanced adventitious root formation in the euhalophyte Suaeda salsa. Funct. Plant Biol. 43, 244–253. doi: 10.1071/FP15120

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosgrove, D. J. (2015). Plant expansins: diversity and interactions with plant cell walls. Curr. Opin. Plant Biol. 25, 162–172. doi: 10.1016/j.pbi.2015.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Davenport, R. J., Muoz-Mayor, A., Jha, D., Essah, P. A., Rus, A., and Tester, M. (2010). The Na+ transporter AtHKT1;1 controls retrieval of Na+ from the xylem in Arabidopsis. Plant Cell Environ. 30, 497–507. doi: 10.1111/j.1365-3040.2007.01637.x

CrossRef Full Text | Google Scholar

Dewey, C. N., and Li, B. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf. 12:323. doi: 10.1186/1471-2105-12-323

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, G., Yang, Q., Ruan, X., Si, T., Yuan, B., Zheng, W., et al. (2022). Proteomics analysis of the effects for different salt ions in leaves of true halophyte Sesuvium portulacastrum. Plant Physiol. Biochem. 170, 234–248. doi: 10.1016/j.plaphy.2021.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dou, L., He, K. K., Higaki, T., Wang, X., and Mao, T. (2018). Ethylene signaling modulates cortical microtubule reassembly in response to salt stress. Plant Physiol. 176, 2071–2081. doi: 10.1104/pp.17.01124

PubMed Abstract | CrossRef Full Text | Google Scholar

Duszyn, M., Świeżawska, B., Szmidt-Jaworska, A., and Jaworski, K. (2019). Cyclic nucleotide gated channels (CNGCs) in plant signaling: current knowledge and perspectives. J. Plant Physiol. 241:153035. doi: 10.1016/j.jplph.2019.153035

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellouzi, H., Hamed, K. B., Cela, J., Munné-Bosch, S., and Abdelly, C. (2011). Early effects of salt stress on the physiological and oxidative status of Cakile maritima (halophyte) and Arabidopsis thaliana (glycophyte). Physiol. Plant. 142, 128–143. doi: 10.1111/j.1399-3054.2011.01450.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Endler, A., Kesten, C., Schneider, R., Zhang, Y., Ivakov, A., Froehlich, A., et al. (2015). A mechanism for sustained cellulose synthesis during salt stress. Cells 162, 1353–1364. doi: 10.1016/j.cell.2015.08.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Fagerstedt, K. V., Kukkola, E. M., Koistinen, V., Takahashi, J., and Marjamaa, K. (2010). Cell wall lignin is polymerised by class III secretable plant peroxidases in Norway spruce. J. Integr. Plant Biol. 52, 186–194. doi: 10.1111/j.1744-7909.2010.00928.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Flowers, T. J., Galal, H. K., and Bromham, L. (2010). Evolution of halophytes: multiple origins of salt tolerance in land plants. Funct. Plant Biol. 37, 604–612. doi: 10.1038/79602

CrossRef Full Text | Google Scholar

Flowers, T. J., Rana, M., and Colmer, T. D. (2015). Sodium chloride toxicity and the cellular basis of salt tolerance in halophytes. Ann Bot. 3, 419–431. doi: 10.1093/aob/mcu217

CrossRef Full Text | Google Scholar

Food and Agricultural Organization (FAO) (2015). Status of the World’s Soil Resources. FAO, Rome.

Google Scholar

Gaxiola, R. A., Palmgren, M. G., and Schumacher, K. (2007). Plant proton pumps. FEBS Lett. 581, 2204–2214. doi: 10.1016/j.febslet.2007.03.050

CrossRef Full Text | Google Scholar

Gill, S. S., and Tuteja, N. (2010). Polyamines and abiotic stress tolerance in plants. Plant Signal. Behav. 5, 26–33. doi: 10.4161/psb.5.1.10291

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yasssour, M., Levin, J. Z., Thompson, D. A., Ido, A., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Grene, A. R., Neval, E., and Heath, L. S. (2002). Role of superoxide dismutases (SODs) in controlling oxidative stress in plants. J. Exp. Bot. 53, 1331–1341. doi: 10.1093/jexbot/53.372.1331

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, K. M., Babourina, O., Christopher, D. A., Borsics, T., and Rengel, Z. (2008). The cyclic nucleotide-gated channel, AtCNGC10, influences salt tolerance in Arabidopsis. Physiol. Plant. 134, 499–507. doi: 10.1111/j.1399-3054.2008.01157.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasanuzzaman, M., Nahar, K., Alam, M. M., Bhowmik, P. C., and Fujita, M. (2014). Potential use of halophytes to remediate saline soils. J. Biomed. Biotechnol. 2014:589341. doi: 10.1155/2014/589341

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasegawa, P. M., Bressan, R. A., Zhu, J. K., and Bohnert, H. J. (2000). Plant cellular and molecular responses to high salinity. Annu. Rev. Plant Physiol. Plant Mol. Biol. 51, 463–499. doi: 10.1146/annurev.arplant.51.1.463

CrossRef Full Text | Google Scholar

Hayat, K., Bundschuh, J., Jan, F., Menhas, S., Hayat, S., Haq, F., et al. (2020). Combating soil salinity with combining saline agriculture and phytomanagement with salt-accumulating plants. Crit. Rev. Env. Sci. Tec. 50, 1085–1115. doi: 10.1080/10643389.2019.1646087

CrossRef Full Text | Google Scholar

Haynes, R. (1980). A comparison of two modified Kjeldahl digestion techniques for multi-element plant analysis with conventional wet and dry ashing methods. Commun. Soil Sci. Plant Anal. 11, 459–467. doi: 10.1080/00103628009367053

CrossRef Full Text | Google Scholar

He, W., Wang, D., Yang, N., Cao, D., Chen, X., Chen, J., et al. (2022). In vitro shoot culture of Sesuvium portulacastrum: an important plant for phytoremediation. Agriculture 12:47. doi: 10.3390/agriculture12010047

CrossRef Full Text | Google Scholar

Hoagland, D. R., and Synder, W. C. (1933). Nutrition of strawberry plants under controlled conditions. Proc. Am. Soc. Hortic. Sci. 30, 288–294.

Google Scholar

Imadi, S. R., Shah, S. W., Kazi, A. G., Azooz, M., and Ahmad, P. (2016). “Phytoremediation of saline soils for sustainable agricultural productivity,” in Plant Metal Interaction. ed. P. Ahmad (New York, NY: Elsevier), 455–468.

Google Scholar

Jayakumar, B., Ana, R. M., and Sergey, S. (2014). ROS homeostasis in halophytes in the context of salinity stress tolerance. J. Exp. Bot. 65, 1241–1257. doi: 10.1093/jxb/ert430

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, X., Leidi, E. O., and Pardo, J. M. (2010). How do vacuolar NHX exchangers function in plant salt tolerance? Plant Signal. Behav. 5, 792–795. doi: 10.4161/psb.5.7.11767

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, H., Pardo, J. M., Batelli, G., Van, O. M. J., Bressan, R. A., and Li, X. (2013). The salt overly sensitive (SOS) pathway: established and emerging roles. Mol. Plant 6, 275–86. doi: 10.1093/mp/sst017

CrossRef Full Text | Google Scholar

Karakas, S., Dikilitas, M., and Tıpırdamaz, R. (2020). “Phytoremediation of salt-affected soils using halophytes,” in Handbook of Halophytes. ed. M. N. Grigore (Cham: Springer).

Google Scholar

Kawakami, Y., Imran, S., Katsuhara, M., and Tada, Y. (2020). Na+ transporter svhkt 1;1 from a halophytic turf grass is specifically upregulated by high Na+ concentration and regulates shoot Na+ concentration. Int. J. Mol. Sci. 21:6100. doi: 10.3390/ijms21176100

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiani-Pouya, A., Roessner, U., Jayasinghe, N. S., Lutz, A., Rupasinghe, T., Bazihizina, N., et al. (2017). Epidermal bladder cells confer salinity stress tolerance in the halophyte quinoa and Atriplex species. Plant Cell Environ. 40, 1900–1915. doi: 10.1111/pce.12995

PubMed Abstract | CrossRef Full Text | Google Scholar

Ladwig, F., Dahlke, R. I., Stuhrwohldt, N., Hartmann, J., Harter, K., and Sauter, M. (2015). Phytosulfokine regulates growth in Arabidopsis through a response module at the plasma membrane that includes CYCLIC NUCLEOTIDE-GATED CHANNEL17, H+-ATPase, and BAK1. Plant Cell 27, 1718–1729. doi: 10.1105/tpc.15.00306

PubMed Abstract | CrossRef Full Text | Google Scholar

Langfelder, P., and Horvath, S. (2009). WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. 9:559. doi: 10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2002). Analysis of relative gene expression data using real-time quantitative PCR. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

CrossRef Full Text | Google Scholar

Lokhande, V. H., Gor, B. K., Desai, N. S., Nikam, T. D., and Suprasanna, P. (2013). Sesuvium portulacastrum, a plant for drought, salt stress, sand fixation, food and phytoremediation. A review. Agron. Sustain. Dev. 33, 329–348. doi: 10.1007/s13593-012-0113-x

CrossRef Full Text | Google Scholar

Lokhande, V. H., Nikam, T. D., and Penna, S. (2010). Differential osmotic adjustment to iso-osmotic NaCl and PEG stress in the in vitro cultures of Sesuvium portulacastrum (L.). J. Crop Sci. Biotech. 13, 251–256. doi: 10.1007/s12892-010-0008-9

CrossRef Full Text | Google Scholar

Lonard, R. I., and Judd, F. W. (1997). The biological flora of coastal dunes and wetlands. Sesuvium portulacastrum (L.) L. J. Coastal Res. 13, 96–104. doi: 10.2307/4298595

CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Bio. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Mariem, W., Kilani, B. R., Benet, G., Abdelbasset, L., Stanley, L., Charlotte, P., et al. (2014). How does NaCl improve tolerance to cadmium in the halophyte Sesuvium portulacastrum? Chemosphere 117, 243–250. doi: 10.1016/j.chemosphere.2014.07.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Massange-Sánchez, J. A., Palmeros-Suárez, P. A., Espitia-Rangel, E., Rodríguez-Arévalo, I., Sánchez-Segura, L., Martínez-Gallardo, N. A., et al. (2016). Overexpression of grain amaranth (Amaranthus hypochondriacus) AhERF or AhDOF transcription factors in Arabidopsis thaliana increases water deficit- and salt-stress tolerance, respectively, via contrasting stress-amelioration mechanisms. PLoS One 11:e0164280. doi: 10.1371/journal.pone.0164280

PubMed Abstract | CrossRef Full Text | Google Scholar

Mian, A., Oomen, R., Isayenkov, S., Sentenac, H., Maathuis, F., and Véry, A. (2011). Over-expression of an Na+-and K+-permeable HKT transporter in barley improves salt tolerance. Plant J. 68, 468–479. doi: 10.1111/j.1365-313X.2011.04701.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Munir, N., Hasnaina, M., Roessnerb, U., and Abideen, Z. (2022). Strategies in improving plant salinity resistance and use of salinity resistant plants for economic sustainability. Crit. Rev. Environ. Sci. Technol. 52, 2150–2196. doi: 10.1080/10643389.2021.1877033

CrossRef Full Text | Google Scholar

Orsini, F., DUrzo, M. P., Inan, G., Serra, S., Oh, D. H., Mickelbart, M. V., et al. (2010). A comparative study of salt tolerance parameters in 11 wild relatives of Arabidopsis thaliana. J. Exp. Bot. 61, 3787–3798. doi: 10.1093/jxb/erq188

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, C., Chang, L., Yang, Q., Tong, Z., Wang, D., Tan, Y., et al. (2018). Comparative physiological and proteomic analyses of the chloroplasts in halophyte Sesuvium portulacastrum under differential salt conditions. J. Plant Physiol. 232, 141–150. doi: 10.1016/j.jplph.2018.10.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Z., Lu, Q., and Verma, D. (1996). Reciprocal regulation of delta 1-pyrroline-5-carboxylate synthetase and proline dehydrogenase genes controls proline levels during and after osmotic stress in plants. Mol. Gen. Genet. 253, 334–341. doi: 10.1007/pl00008600

PubMed Abstract | CrossRef Full Text | Google Scholar

Plank, C. O. (1992). Plant analysis reference procedures for the southern region of the United States. South. Coop. Ser. Bull. 368:9. Athens, GA: The Georgia Agricultural of Experiment Stations.

Google Scholar

Qadir, M., Ghafoor, A., and Murtaza, G. (2001). Amelioration strategies for saline soils: a review. Land Degrad. Dev. 12, 357–386. doi: 10.1002/1099-145X(200011/12)11:63.0.CO;2-S

CrossRef Full Text | Google Scholar

Ramani, B., Reeck, T., Debez, A., Stelzer, R., Huchzermeyer, B., Schmidt, A., et al. (2006). Aster tripolium L. and Sesuvium portulacastrum L.: two halophytes, two strategies to survive in saline habitats. Plant Physiol. Biochem. 44, 395–408. doi: 10.1016/j.plaphy.2006.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Rana, M., and Mark, T. (2008). Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 59, 651–681. doi: 10.1146/annurev.arplant.59.032607.092911

CrossRef Full Text | Google Scholar

Ruan, C. J., da Silva, J. A. T., Mopper, S., Qin, P., and Lutts, S. (2010). Halophyte improvement for a salinized world. Crit. Res. Plant Sci. 29, 329–359. doi: 10.1080/07352689.2010.524517

CrossRef Full Text | Google Scholar

Seckin, B., Turkan, I., Sekmen, A. H., and Ozfidan, C. (2010). The role of antioxidant defense systems at differential salt tolerance of Hordeum marinum Huds. (sea barleygrass) and Hordeum vulgare L. (cultivated barley). Environ. Exp. Bot. 69, 76–85. doi: 10.1016/j.envexpbot.2010.02.013

CrossRef Full Text | Google Scholar

Shabala, S., Hariadi, Y., and Jacobsen, S. E. (2013). Genotypic difference in salinity tolerance in quinoa is determined by differential control of xylem Na+ loading and stomatal density. J. Plant Physiol. 170, 906–914. doi: 10.1016/j.jplph.2013.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahzad, B., Fahad, S., Tanveer, M., Saud, S., and Khan, I. A. (2019). “Plant responses and tolerance to salt stress,” in Approaches for Enhancing Abiotic Stress Tolerance in Plants. eds. M. Hasanuzzaman, K. Nahar, M. Fujita, H. Oku, and M. T. Islam (Boca Raton, Florida, USA: CRC press,Taylor & Francis), 61–77.

Google Scholar

Shalata, A., Mittova, V., Volokita, M., Guy, M., and Tal, M. (2001). Response of the cultivated tomato and its wild salt-tolerant relative Lycopersicon pennellii to salt-dependent oxidative stress: the root antioxidative system. Physiol. Plant. 112, 487–494. doi: 10.1080/10715760290006402

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, H., Ishitani, M., Kim, C., and Zhu, J. K. (2000). The Arabidopsis thaliana salt tolerance gene SOS1 encodes a putative Na+/H+ antiporter. Proc. Natl. Acad. Sci. U. S. A. 97, 6896–6901. doi: 10.1073/pnas.120170197

PubMed Abstract | CrossRef Full Text | Google Scholar

Shigeru, S., Takahiro, I., Masahiro, T., Yoshiko, M., Toru, T., Yukinori, Y., et al. (2002). Regulation and function of ascorbate peroxidase isoenzymes. J. Exp. Bot. 53, 1305–1319. doi: 10.1093/jexbot/53.372.1305

CrossRef Full Text | Google Scholar

Singh, M., Kumar, J., Singh, S., Singh, V. P., and Prasad, S. M. (2015). Roles of osmoprotectants in improving salinity and drought tolerance in plants: a review. Rev. Environ. Sci. Biotechnol. 14, 407–426. doi: 10.1007/s11157-015-9372-8

CrossRef Full Text | Google Scholar

Székely, G., Brahám, E., Csépl, G., Rigó, G., Zsigmond, L., Csiszár, J., et al. (2008). Duplicated P5CS genes of Arabidopsis play distinct roles in stress regulation and developmental control of proline biosynthesis. Plant J. 53, 11–28. doi: 10.1111/j.1365-313X.2007.03318.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tenhaken, R. (2015). Cell wall remodeling under abiotic stress. Front. Plant Sci. 5:771. doi: 10.3389/fpls.2014.00771

PubMed Abstract | CrossRef Full Text | Google Scholar

Vadim, D. (2018). ROS-activated ion channels in plants: biophysical characteristics, physiological functions and molecular nature. Int. J. Mol. Sci. 19:1263. doi: 10.3390/ijms19041263

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Chang, L., Wang, B., Wang, D., Li, P., Wang, L., et al. (2013). Comparative proteomics of Thellungiella halophila leaves under different salinity revealed chloroplast starch and soluble sugar accumulation played important roles in halophyte salt tolerance. Mol. Cell. Proteomics 12, 2174–2195. doi: 10.1074/mcp.M112.022475

PubMed Abstract | CrossRef Full Text | Google Scholar

Willekens, H., Chamnongpol, S., Davey, M., Schraudner, M., Langebartels, C., Montagu, M. V., et al. (2014). Catalase is a sink for H2O2 and is indispensable for stress defence in C3 plants. EMBO J. 16, 4806–4816. doi: 10.1093/emboj/16.16.4806

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, J., Liu, Y., Yang, L., He, H., Huang, Y., Fang, L., et al. (2021). Cell wall β-1, 4-galactan regulated by BPC1/BPC2-GALS1 module aggravates salt sensitivity in Arabidopsis thaliana. Mol. Plant 14, 411–425. doi: 10.1016/j.molp.2020.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, X., Sun, Y., Yang, Q., Guo, A., Chang, L., Wang, D., et al. (2014). Quantitative proteomics of Sesuvium portulacastrum leaves revealed that ion transportation by V-ATPase and sugar accumulation in chloroplast played crucial roles in halophyte salt tolerance. J. Proteome 99, 84–100. doi: 10.1016/j.jprot.2014.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zelm, E., Zhang, Y., and Testerink, C. (2020). Salt tolerance mechanisms of plants. Annu. Rev. Plant Biol. 71, 403–433. doi: 10.1146/annurev-arplant-050718-100005

CrossRef Full Text | Google Scholar

Zhang, C., Wang, D., He, W., Liu, H., Chen, J., Wei, X., et al. (2022). Sesuvium portulacastrum-mediated removal of nitrogen and phosphorus affected by sulfadiazine in aquaculture wastewater. Antibiotics 11:68. doi: 10.3390/antibiotics11010068

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Zhang, H., Song, C., Zhu, J. K., and Shabala, S. (2020). Mechanisms of plant responses and adaptation to soil salinity. Innovation (N Y) 1:100017. doi: 10.1016/j.xinn.2020.100017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, M., Song, R., Wang, Y., Shu, S., and Guo, S. (2020). TGase regulates salt stress tolerance through enhancing bound polyamines-mediated antioxidant enzymes activity in tomato. Environ. Exp. Bot. 179:104191. doi: 10.1016/j.envexpbot.2020.104191

CrossRef Full Text | Google Scholar

Zhou, Y., Yin, X., Duan, R., Hao, G., Guo, J., Jiang, X., et al. (2015). SpAHA1 and SpSOS1 coordinate in transgenic yeast to improve salt tolerance. PLoS One 10:e0137447. doi: 10.1371/journal.pone.0137447

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Yin, X., Wan, S., Hu, Y., Xie, Q., Li, R., et al. (2018). The Sesuvium portulacastrum plasma membrane Na+/H+ antiporter SpSOS1 complemented the salt sensitivity of transgenic Arabidopsis sos 1 mutant plants. Plant Mol. Bio. Rep. 36, 553–563. doi: 10.1007/s11105-018-1099-6

CrossRef Full Text | Google Scholar

Keywords: salt tolerance, Sesuvium portulacastrum, transcriptomic analysis, differentially expressed genes, ion transport

Citation: Wang D, Yang N, Zhang C, He W, Ye G, Chen J and Wei X (2022) Transcriptome analysis reveals molecular mechanisms underlying salt tolerance in halophyte Sesuvium portulacastrum. Front. Plant Sci. 13:973419. doi: 10.3389/fpls.2022.973419

Received: 20 June 2022; Accepted: 23 August 2022;
Published: 23 September 2022.

Edited by:

Amr Adel Elkelish, Suez Canal University, Egypt

Reviewed by:

Ganesh Chandrakant Nikalje, R. K. Talreja College of Arts, Science and Commerce, India
Bhaskar Gupta, Government General Degree College, Singur, India

Copyright © 2022 Wang, Yang, Zhang, He, Ye, Chen and Wei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jianjun Chen, ampjaGVuQHVmbC5lZHU=; Xiangying Wei, eGlhbmd5aW5nd2VpQG1qdS5lZHUuY24=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.