Integrated Analysis of Transcriptome and Metabolome Reveals Molecular Mechanisms of Rice with Different Salinity Tolerances

Rice is a crucial global food crop, but it lacks a natural tolerance to high salt levels, resulting in significant yield reductions. To gain a comprehensive understanding of the molecular mechanisms underlying rice’s salt tolerance, further research is required. In this study, the transcriptomic and metabolomic differences between the salt-tolerant rice variety Lianjian5 (TLJIAN) and the salt-sensitive rice variety Huajing5 (HJING) were examined. Transcriptome analysis revealed 1518 differentially expressed genes (DEGs), including 46 previously reported salt-tolerance-related genes. Notably, most of the differentially expressed transcription factors, such as NAC, WRKY, MYB, and EREBP, were upregulated in the salt-tolerant rice. Metabolome analysis identified 42 differentially accumulated metabolites (DAMs) that were upregulated in TLJIAN, including flavonoids, pyrocatechol, lignans, lipids, and trehalose-6-phosphate, whereas the majority of organic acids were downregulated in TLJIAN. The interaction network of 29 differentially expressed transporter genes and 19 upregulated metabolites showed a positive correlation between the upregulated calcium/cation exchange protein genes (OsCCX2 and CCX5_Ath) and ABC transporter gene AB2E_Ath with multiple upregulated DAMs in the salt-tolerant rice variety. Similarly, in the interaction network of differentially expressed transcription factors and 19 upregulated metabolites in TLJIAN, 6 NACs, 13 AP2/ERFs, and the upregulated WRKY transcription factors were positively correlated with 3 flavonoids, 3 lignans, and the lipid oleamide. These results suggested that the combined effects of differentially expressed transcription factors, transporter genes, and DAMs contribute to the enhancement of salt tolerance in TLJIAN. Moreover, this study provides a valuable gene–metabolite network reference for understanding the salt tolerance mechanism in rice.


Introduction
Soil salinization is one of the most significant abiotic stresses that affects various crop growth, development, and yield worldwide [1]. Except for in coastal areas, improper irrigation, the use of poor-quality water, and climate change particularly aggravate the salinization of arid and semi-arid soil in inland areas, which generally impairs the ability of the root to absorb water, leading to the disruption of metabolic processes and reduced photosynthetic efficiency [2,3]. Under salinity conditions, osmotic, ionic, and oxidative stresses affect most plants, inducing morphological and biochemical changes [4]. Plants respond to salt stress through many strategies, including the selective preferential uptake and exclusion of K + ions over Na + ions [5] and the elimination of reactive oxygen species (ROS) CCX2 and CCX5_Ath, AAP3_Ath, AB2E_Ath and OCT7_Ath) in salt-tolerant rice TLJIAN were positively correlated with multiple upregulated differential metabolites in this salttolerant rice. Furthermore, the upregulated transcription factors, such as NAC family, AP2/ERF family, and WRKY, in the TLJIAN were also positively correlated with the coordinately upregulated three flavonoid metabolites (dihydrokaempferide, isovitexin-8-O-xyloside, kaempferol -3-O-arabinoside-7-O-rhamnoside), three lignans and coumarins metabolites (matairesinol, matairesinoside, and lariciresinol-4 -O-glucoside), as well as one lipid metabolite (Oleamide) in this salt-tolerant variety. These findings will provide technical support for the exploration of the potential mechanism of salt tolerance and the selection of salt-tolerant rice varieties.

Transcriptome Differences of Rice Varieties with Different Salinity Tolerance
To investigate the disparity in salt tolerance between two types of rice, namely salttolerant Lianjian 5 (TLJIAN) and salt-sensitive Huajing 5 (HJING), a comparative analysis of the transcriptomes was conducted. A total of 1518 differentially expressed genes (DEGs) were identified in the two types of rice cultivated in saline-alkali land, of which 882 genes were upregulated and 636 genes were downregulated in salt-tolerant rice (Table S1). Among the DEGs, 935 of them had clear annotation information, including 66 transcription factor TFs, 33 transporters, and 61 receptor-like-kinase. Interestingly, 46 of these genes have been previously reported to be associated with salt-tolerance (Table S2).
To gain a deeper understanding of the resulting DEGs functionality, MapMan software was used to locate these DEGs to the modules and biochemical or molecular biology pathways of the Arabidopsis Seed-Molecular Networks blueprint. In general, almost all differentially expressed transcription factors, such as NAC, WRKY, MYB, and AP2/ERE family, were upregulated in salt-tolerant rice, and OsbHLH6, OsWRKY22, ERF17_Ath, OsDREB1B, OsDREB1C, OsDREB1F, and OsDREB1G were significantly upregulated among the differentially expressed transcription factors in the major classes. The genes of stressrelated DEGs, such as cold, drought/salt, and heat, were also upregulated in salt-tolerant rice, and the precursor synthesis genes of cell wall modules were all upregulated in salttolerant rice ( Figure 1, Table 1). Both receptor-like protein kinase FERONIA (Os04g0586500, Os01g0769700) were upregulated in salt-tolerant rice ( Figure 1). In addition, differentially expressed transporter genes AB2E_Ath, OCT7_Ath, HIP7_Ath, NRAM6, and AB3C_Ath were significantly upregulated in salt-tolerant rice ( Table 2).
Compared to salt-sensitive differentially expressed genes, all three small molecules of heat shock protein were upregulated in salt-tolerant rice; all Early Responsive to Dehydration, Scarecrow-like Protein, and Accelerated Cell Death were also upregulated in salt-tolerant rice, but photosynthesis, fermentation genes and five of the seven genes involved in cell wall modification were downregulated in salt-tolerant rice. Peroxidase, UDP-glycosyltransferase, and Xylanase inhibitor were all raised, whereas 8 of the 9 genes involved in auxin synthesis and signaling were downregulated, as was the case for gibberellin, brassinosteroid, and jasmonate ( Figure 1).
To gain insight into the biological significance of DEGs, GO and KEGG enrichment analyses were performed. GO enrichment analysis showed that DEGs were enriched in response to stress, maintaining the integrity of cell anatomy, catalytic activity, and transcriptional regulatory activity (Figure 2A), as well as enriched in response to environmental signal response processes, especially those related to immune responses ( Figure 2B). KEGG pathway enrichment analysis further confirmed that DEGs were mainly enriched in anabolic pathways such as environmental signaling response plant-pathogen interaction, plant hormone signaling, MAPK signaling pathway, and amino sugar and nucleotide sugar metabolism ( Figure 2C).
Gene locus in bold means gene reported in previous studies, and the number in brackets represents the related reference. Log 2 FC means log 2 foldchange (TLJIAN/HJING).

Metabolome Differences of Rice Varieties with Different Salt Tolerance
Plants required the production of numerous metabolites during growth and development to adapt to environmental changes and withstand various stresses [46]. To examine the metabolism differences between TLJIAN and HJING rice varieties under salt stress, UPLC-MS/MS was used to perform metabolomic analysis. A total of 1336 metabolites were detected in both types of rice, with 42 metabolites showing upregulation (|log2FC| ≥ 1, VIP > 1) in the salt-tolerant TLJIAN 5 compared with salt-sensitive HJING. Notably, the flavonoids dihydrokaempferide and isovitexin-8-O-xyloside, the lignans matairesinol and matairesinoside, the amino acid derivative N-Acetyl-L-Tryptophan, and the alkaloids octadecadienoic acid amide were significantly upregulated in the salt-tolerant rice. In addition, pyrocatechol, catechin, trehalose 6-phosphate, naringin, quercetin-5-O-β-D-glucoside, and lysophosphatidylcholine (20:5 and 16:4) were also found to be upregulated in the salt-tolerant rice, whereas the majority of organic acids were downregulated in salt-tolerant rice (Table 3).

Conjoint Analysis of Transcriptomic Difference and Metabolomic Fluctuation Accompanied with Different Salt Tolerance in Rice
In order to understand the difference mechanism of salt stress response between the two types of rice, a conjoint analysis of its transcriptome and metabolome data was performed in this study. The differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were mapped onto the KEGG pathway map to identify the pathways enriched with DEGs and DAMs. The major pathway with the highest enrichment was found to be amino acid biosynthesis, starch and sucrose metabolism, and phenylpropane biosynthesis pathways (Figure 3). Among them, amino acid biosynthesis was the pathway with the highest degree of differential enrichment of the KEGG pathway, including 11 DAMs and DEGs, respectively. It could be seen that, during the process of salt stress, the process of amino acid biosynthesis had the greatest influence on the root system of rice, and played an important role in regulating the process of salt tolerance stress in rice. To further study the interaction between transcriptome and metabolome more systematically, gene-metabolite combinations with coefficient correlation (PCC > 0.80 and p-value < 0.05) were screened, and then stress-tolerance-related genes and metabolites were selected to obtain 29 differentially expressed transporter genes (Table S3), 63 differentially expressed transcription factors, and 19 differentially expressed upregulated metabolites (Table S4). Thus, transporter gene-metabolite and transcription-factor-metabolite-related network maps were constructed to show the relationship between salt-tolerance-related transcription factors and transporter genes and differential metabolite accumulation in the transcriptome under salt stress.

Discussion
Under environmental stress, plants have evolved a variety of protective mechanisms, such as changing the content of related metabolites and regulating the expression of DEGs to resist stress damage, but a comprehensive understanding of their complex protection mechanisms is still an increasingly urgent need [47,48]. Currently, a range of omics techniques have emerged as effective strategies for systematically revealing the stress response mechanism [46]. Genes responsive to salinity have been categorized based on their functions, including transcription factors (TF), transporters, carbohydrates, protein metabolism, energy metabolism, hormones, reactive oxygen network components, osmoprotectants, cell walls, and signal transduction components [49]. In this study, a comparative transcriptome comparative analysis of different salt-tolerant rice varieties showed that 1518 DEGs of salt-tolerant TLJIAN responded to salt stress compared to salt-sensitive HJING, and these DEGs were functionally categorized into the aforementioned salt-responsive clusters (Figure 1, Tables 2 and S2).
The intensive studies indicated that various TFs are involved in improving salinity tolerance in rice [24]. Major Transcription factors (TFs) families regulated salt tolerance in rice such as dehydration-responsive element binding protein (DREB), ABA-responsive element binding protein/factor, and NAC [9]. The expression of the MYB, bHLH, NAC, ERF, WRKY, and bZIP transcription factors was altered under salt stress [50]. Intriguingly, almost all differentially expressed transcription factors, such as the NAC, WRKY, MYB, and AP2/ERE family, were upregulated in salt-tolerant rice, and several of them were demonstrated to be positive controllers of salinity tolerance in previous studies (Figure 1, Tables 2 and S2). ONAC022, a stress-responsive transcriptional activator, improved drought and salt stress tolerance in rice through modulating an ABA-mediated pathway [32]. OsNAC14 mediated drought tolerance by recruiting factors involved in DNA damage repair and defense response resulting in improved tolerance to drought [51]. The important regulator of ethylene signal transduction, a transcription factor of the AP2/ERF family, also played an important role in plant growth and development and responded to biotic and abiotic stresses, such as the salt-treatment-induced expression of five ERF genes (SpERF) in tomatoes [13]. In this study, the AP2/ERF family showed an enhanced expression of members (including ERF17_Ath, OsDREB1B, OsDREB1C, OsDREB1F, and OsDREB1G) in salt-tolerant rice (Table 1). OsDREB1C-overexpressing rice boosted grain yields, shortened the growth duration, and conferred higher nitrogen use efficiency and early flowering [52]. The overexpression of a rice OsDREB1F gene increased salt, drought, and low temperature tolerance in both Arabidopsis and rice through the ABA-dependent pathway [34].
Besides these positive regulators, certain genes act as the negative controllers of salinity tolerance in rice [24]. Notably, salt-tolerant rice TLJIAN showed the downregulation of several known salt-tolerance-negative regulator genes, including Drought and Salt Tolerance Coactivator 1 (DCA1) (Table S2), a negative regulator of salt stress, which is consistent with the report where the downregulation of DCA1 enhanced drought and salt tolerance [42]. Incremental studies suggest that numerous TFs and their interacting proteins have been implicated in rice salt stress response via regulating a series of signaling pathways [11]. Therefore, further exploration of these regulons, their interacting proteins, and target genes will contribute to understanding the network of signaling pathways involved in salt stress response and accelerate the elucidation of salt-tolerance mechanisms [11].
In addition, compared with salt-sensitive rice, many upregulated transporter genes in response to salt stress were found in salt-tolerant rice, among which two ABC transporter (AB3C_Ath and AB2E_Ath) genes were significantly upregulated. In plants, ABC transporters play an important role in stress tolerance [53]. Transgenic Arabidopsis plants that overexpress an ATP-binding cassette (ABC) transporter, AtABCG36/AtPDR8, were more resistant to drought and salt stress and grew to higher shoot fresh weight (FW) than the wild-type [54]. The expression of the ABCG gene in xylem, leaf, and apical buds was upregulated under NaCl stress in Betula halophila [53]. Therefore, the investigation of whether these transporters is regulated by salt-tolerant rice enhanced TFs (e.g., NACs, OsDREB1 subfamily) under salt stress is of particular interest and warrants further research. Regarding the receptor-like protein kinase, two FERONIA (FER) members (Os04g0586500, Os01g0769700) were coordinately upregulated in salt-tolerant rice (Figure 1). FERONIA coordinated plant growth and salt tolerance via the phosphorylation of phytochrome B [55]. Under low temperature stress conditions, the FERONIA receptor kinase MdMRLK2 could enhance the cold tolerance of apple plants by increasing the accumulation of osmoregulatory substances in cells, the synthesis of anthocyanins, and maintaining the stability of cell wall components [56].
Some recent studies in metabolomics have reported that salt stress induced the accumulation of some compounds in D. officinale leaves, especially flavonoids, sugars, and alkaloids, which may play an important role in the salt-stress responses of leaf tissues from D. officinale [22]. Soil salt-alkaline stress enhanced the metabolites accumulation of amino acids, alkaloids, organic acids, and polyphenols contents in wolfberry fruits [21]. Furthermore, Wei et al. demonstrated that the overexpression of the phenylcoumaran benzylic ether reductase gene (PtPCBER) in the phenylpropane metabolic pathway improved the salt tolerance of poplar through lignan-mediated reactive oxygen species clearance [57]. Studies of mungbean reported that the contents of phenylpropanoid-derived metabolites, including catechin, chlorogenic acid, isovitexin, p-coumaric acid, syringic acid, ferulic acid, and vitexin, significantly increased under salinity conditions [58]. The salt stress study of Arabidopsis found that MYB111 was a positive regulator in salt stress response, and that the additions of bioflavonoids and chalcone/dihydrokaempferole/quercetin were able to rescue the loss of salt tolerance in myb111 mutants, which proved that flavonoids were crucial against salt stress [59].
Sugars are energy transporters and essential hydrophilic solutes that protect cellular proteins and membranes from extreme environmental challenges. Trehalose, an important osmolyte and osmo-protectant, and its intermediate compounds effectively modulate salt response and salt tolerance in different plants [60]. The variations in trehalose and trehalose-6-phosphate concentrations impact many biological functions, including adaptive stress responses [61]. In this study, differential metabolites in the flavonoids dihydrokaempferide, isovitexin-8-O-xyloside and catechin, and lignans matairesinol and matairesinoside, as well as trehalose-6-phosphate, were up-regulated in salt-tolerant rice (TLJIAN); meanwhile, the majority of organic acids were downregulated in salt-tolerant rice (Table 3). Under stress circumstances, the reduction in organic acid intermediates may be interpreted as a general consequence of decreased de novo CO 2 assimilation of carbon dioxide under stress because of stomatal limitation [25]. When compared to the salt-sensitive rice IR64, the salt-tolerant variety FL478 displayed more significant increases in amino acids and sugars, but more noticeable decreases in organic acids in both leaves and roots [62]. In Lotus japonicus, the simultaneous decrease in the concentration of most organic acids may serve multiple purposes in maintaining vitality under salt stress [25]. Therefore, the function of fluctuated metabolites between TLJIAN and HJING identified here in rice salt tolerance deserves in-depth study.
To better understand the relationship between differential genes and differential metabolites and their role in salt tolerance, combined transcriptome and metabolome analysis may be the best approach [63]. The combined transcriptome and LC-ESI-MS/MS metabolomics analysis revealed that melatonin enhanced plant tolerance to multiple stresses by mediating flavonoid biosynthesis, providing new ideas for studying the crosstalk between melatonin and flavonoids [64]. Previous studies have reported that trehalose-6phosphate synthase (TPS) plays a crucial role in the biosynthesis of trehalose. Furthermore, rice OsTPS1 and OsTPS8 have been found to positively regulate salt tolerance by enhancing the accumulation of trehalose and proline in rice overexpression plants [65,66]. In this study, an up-regulation of trehalose-6-phosphate was observed in salt-tolerant TLJIAN (Table 3), concomitant with the increased expression of trehalose-6-phosphate synthase OsTPS1 (Table S4).
Transcriptome and metabolomic analysis found that genes involved in the synthesis pathways of phenylpropanin, ferulic acid (FA), and spermidine (Spd) in tomatoes, especially transcription factors such as MYB, Dof, BPC, and AP2/ERF, were expressed under salt stress [67]. Combined transcriptome and metabolome analysis under watermelon salt stress showed that genes involved in phenylpropanin synthesis, hormone signaling, and carbohydrate synthesis pathways, as well as bHLH family transcription factors, played an important role in improving the salt tolerance of watermelon [68]. In this study, the metabolites and genes with significant differences were analyzed for correlation, and by constructing a network map of transporters/transcription factors-metabolites. It was found that transporters, including cation/calcium exchange proteins OsCCX2 and CCX5_Ath, amino acid transporter AAP3_Ath, ABC transporter AB2E_Ath, and organic cation/carnitine transporter OCT7_Ath, were significantly upregulated under salt stress ( Figure 4, Table 2).  Table 1). These results suggested that the transporters and transcription factors helped to improve the salt tolerance of TLJIAN rice.

Plant Materials and Salt Treatment
Two japonica rice (O. sativassp. japonica), including the salt-tolerant rice variety Lianjian 5 (TLJIAN) and salt-sensitive variety Huajing 5 (HJING), were cultivated in the Qingkou salt field experimental base (salt content in soil > 0.5%, electrical conductivity EC~2.32 S m −1 ) at Jiangsu Lianyungang Academy of Agricultural Sciences in 2022. Rice seeds were sown on May 15 and transplanted on June 15. During the growth period, except the tillering-hearty period, the water layer was maintained in the field using the irrigation water (salt content~0.15%). The roots were sampled on August 28 (7 days after the heading stage), placed in sterilized 5 mL centrifuge tubes, and stored in a −80 • C refrigerator before extracting RNA, and three biological replicates were set up per sample.

RNA Extraction and Transcriptome Sequencing
First, TRIzol reagent was used to extract high-quality RNA (no degradation and no DNA contamination) from each sample, then a NanoPhotometer spectrophotometer was used to measure RNA purity (OD 260/280 and OD 260/230 ratios), then the RNA concentration was accurately quantified using a Qubit2.0 Fluorometer, and finally Agilent 2100 was used as a Bioanalyzer to accurately detect RNA integrity. Illumina's NEBNext ® kit was used to build the library with total RNA > 1 ug as the starting RNA. We enriched mRNA with polyA tails using Oligo(dT) magnetic beads. A fragmentation buffer was then added to break the RNA into short fragments, and double-stranded cDNA was synthesized using the short RNA as a template, with three biological replicates. After pooling different libraries according to the effective concentration and 5 Gb data volume, high-throughput sequencing was performed using the Illumina NovaSeq 6000 platform to generate 150 bp paired end sequencing data.
The sequence of the new gene was extracted from the genome, and the new gene was associated with KEGG, GO, NR, using diamond software. Swiss-Prot, trEMBL, KOG database sequence alignment results were annotated with Evalue le-5 [72]. Plant transcription factor prediction used iTAK software, which integrated two databases, PnTFDB and PlantTFDB [73]. The gene alignment was calculated using featureCounts (v1.6.2), and then the FPKM of each gene was calculated based on the gene length. Gene differential expression analysis between the two groups was performed using DESeq2 [74]. The adjusted p values < 0.05 and |log2foldchange| ≥ 1 as a threshold for identifying differentially expressed genes [74].
To analyze the pathways and functions of DEGs related to salt stress, MapMan (v3.6.0 RC1) was utilized to annotate and visualize DEGs on metabolic pathways [75,76].

Enrichment Analysis of Differentially Expressed Genes
Enrichment analysis was performed based on hypergeometric testing, and a hypergeometric distribution test of KEGG was performed on pathway units and genetic annotation was performed [77]. Using GO-Term in the GO database, a hypergeometric test of GO was performed to find GO-Term, which was significantly enriched in differentially expressed genes compared to the entire genome background [78]. Diamond software was used to align protein sequences or cDNA sequences to the KOG database, and then extract the annotations of the KOG database [79].

Sample Preparation and Extraction for Widely Targeted Metabolomics
The three biological samples from TLJIAN and HJING, respectively, were freeze-dried with a vacuum freeze dryer (Scientz-100F), then the freeze-dried samples were crushed using a mixer with zirconia beads (MM 400, Retsch) at 30 Hz for 1.5 min. Next, 50 mg of lyophilized powder was dissolved in 1.2 mL of 70% methanol solution and vortexed for 30 s every 30 min for a total of 6 times. After centrifugation at 12,000 rpm for 3 min, the extract was filtered prior to UPLC-MS/MS analysis.

Metabolomic Analysis
For the quantitative analysis of metabolites, we employed a multiple reaction monitoring (MRM) method using the triple-quadrupole mass spectrometry. QQQ scanning used MRM mode with collision gas (nitrogen) set to medium. Through further DP (decomposition potential) and CE (collision energy) optimization, the DP and CE of each MRM ion pair are obtained. A specific set of MRM ion pairs was monitored during each period based on the metabolites eluted in each epoch [79].
Based on the NIST compound mass spectrometry library and the metabolome database of the Wuhan Metware Company (www.metware.cn, 1 November 2022), the mass spectrometry data were processed using software Analyst 1.6.3 to qualitatively quantify the metabolites of the samples. After the qualitative and quantitative analysis of the detected metabolites, combined with the grouping of specific samples, on the one hand, based on the OPLS-DA results, from the variable importance projection (VIP) of the obtained multivariate analysis OPLS-DA model, the metabolites with differences between different materials were preliminarily screened. On the other hand, the fold change (FC) of metabolites in each group was compared to further screen for differential metabolites. Differentially accumulated metabolites (DAMs) between two rice varieties were determined using VIP ≥ 1 and absolute Log2FC (fold change) ≥ 1 [10,63].

Combined Analysis of Transcriptome and Metabolome
To better understand the interactions between transcriptome and metabolome in rice roots under salt stress, for each differentially accumulated metabolite (DAM) (|log2FC| ≥ 1, VIP > 1), all DEGs were submitted to analyze metabolite-gene correlation according to the Pearson correlation coefficients (PCCs) [63].
DEGs and DAMs were mapped simultaneously to the KEGG pathway database, further according to the KEGG enrichment analysis results of DEGs and DAMs, the KEGG metabolic pathway that enriched the DEGs and DAMs at the same time was obtained. Then, the functional analysis and the correlation analysis of DEGs and DAMs were combined to screen key metabolic pathways, genes, and metabolites [10].

Conclusions
The combined analysis of transcriptome and metabolome in this study provides insights into the molecular mechanisms of salt tolerance in rice landrace TLJIAN. The results suggest that the molecular mechanisms of salt tolerance in rice plants might involve amino acid biosynthesis, starch and sucrose metabolism, and phenylpropane biosynthesis pathways. Additionally, the interaction network analysis reveals that the significantly upregulated transporter genes OsCCX2, CCX5_Ath, AAP3_Ath, AB2E_Ath, and OCT7_Ath in TLJIAN were positively correlated with multiple upregulated differential metabolites. The upregulated transcription factors, such as NAC family, AP2/ERF family, and WRKY, in salt-tolerant rice were positively correlated with specific flavonoid metabolites, lignans and coumarins metabolites, and lipid metabolites, suggesting that the differentially expressed transcription factors, transporter genes, and metabolites played a crucial role in the resistance of salt-tolerant varieties to salt stress. Overall, this study highlights the differences in transcriptome and metabolome between rice with different salinity tolerances, and provides a gene-metabolite network reference for revealing the mechanism of rice salt tolerance.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/plants12193359/s1, Table S1: Differential expression of salt-stressresponsive genes between salt-tolerant rice TLJIAN and salt-sensitive rice HJING. Table S2: Statistics of the reported genes responding to salt stress. Table S3: The salt-stress-responsive differential transporter genes and differential accumulation metabolites between TLJIAN and HJING. Table S4: The salt-stress-responsive differential transcription factors and differential accumulation metabolites between TLJIAN and HJING. Funding: This research was funded by the Lianyungang "521Project" scientific research project and the Financial Grant Support Program of Lianyungang City (QNJJ2201).

Data Availability Statement:
The data generated or analyzed during this study are available in the NCBI database repository (BioProject accession: PRJNA994662).