Genome-wide expression analysis of salt-stressed diploid and autotetraploid Paulownia tomentosa

Paulownia tomentosa is a fast-growing tree species with multiple uses. It is grown worldwide, but is native to China, where it is widely cultivated in saline regions. We previously confirmed that autotetraploid P. tomentosa plants are more stress-tolerant than the diploid plants. However, the molecular mechanism underlying P. tomentosa salinity tolerance has not been fully characterized. Using the complete Paulownia fortunei genome as a reference, we applied next-generation RNA-sequencing technology to analyze the effects of salt stress on diploid and autotetraploid P. tomentosa plants. We generated 175 million clean reads and identified 15,873 differentially expressed genes (DEGs) from four P. tomentosa libraries (two diploid and two autotetraploid). Functional annotations of the differentially expressed genes using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases revealed that plant hormone signal transduction and photosynthetic activities are vital for plant responses to high-salt conditions. We also identified several transcription factors, including members of the AP2/EREBP, bHLH, MYB, and NAC families. Quantitative real-time PCR analysis validated the expression patterns of eight differentially expressed genes. Our findings and the generated transcriptome data may help to accelerate the genetic improvement of cultivated P. tomentosa and other plant species for enhanced growth in saline soils.


Introduction
Plants, as sessile organisms, must cope with a multitude of biotic and abiotic stresses throughout their life cycle. Among the abiotic stresses (e.g., high salinity, extreme temperatures, drought, heavy metal contamination, and nutrient deficiency), salt stress is the most serious threat to sustainable agricultural productivity worldwide, but especially in developing countries [1]. High salinity affects more than 800 million ha of land, which represents more than 6% of the total global land area. This area includes 45 million ha (i.e., approximately 20%) of the 230 million ha of irrigated land [2]. Soil salinity affects plants in two ways: One, high salt

Plant materials and salt treatments
All plant materials used in this study were obtained from the Institute of Paulownia, Henan Agricultural University, Zhengzhou, Henan Province, China. The tissue-cultured diploid and autotetraploid P. tomentosa seedlings were grown for 30 days at 25 ± 2˚C on half-strength Murashige and Skoogmedium [28] with a 16-h photoperiod (light intensity: 130 μmol m −2 s −1 ). The test-tube plantlets were then transferred outdoors into nutritive bowls containing normal garden soil supplemented with 0.11% (quality score) NaCl solution, and grown for another 30 days. Consistently growing healthy seedlings were transferred to larger nutritive bowls with trays underneath. After 50 days, the uniformly growing seedlings were irrigated with 0% (control) or 0.4% (quality score) NaCl solutions. For the salt treatment, NaCl was weighed and divided into three equal parts, which were dissolved in water and poured into the bowls. The water that accumulated in the trays below the bowls was poured back into the bowls. The seedlings were treated with the control or salt solution every 3 days. After all of the salt water was added back into the bowls, the plants were watered every 2 days to maintain the soil moisture content at 75%. After 15 days, the second pair of leaves (fully expanded leaves) from the apex shoot of the control and salt-treated diploid and autotetraploid P. tomentosa plants were collected and immediately frozen in liquid nitrogen. The samples were stored at −80 ± 2˚C. The control diploid and autotetraploid P. tomentosa samples were named PT2 and PT4, respectively, while the salt-treated diploid and autotetraploid samples were named PT2S and PT4S, respectively.

RNA extraction, cDNA library preparation, and sequencing
Total RNA was extracted from PT2, PT4, PT2S, and PT4S leaf samples using the Plant RNA Isolation Kit (AutoLab, Beijing, China) and concentrated using an RNeasy MinElute Cleanup Kit (Qiagen, Valencia, CA, USA). The RNA was then treated with DNase I to eliminate any contaminating genomic DNA. Oligo-d(T) magnetic beads were used to isolate mRNA from the purified total RNA. The mRNA was fragmented in fragmentation buffer (Life Technologies, Beijing, China), and then used as the template for cDNA synthesis. Short fragments were purified and resolved with ethidium bromide buffer for the subsequent end-repair and A (adenine)-tailing. The short fragments were linked to adapters, and analyzed by agarose gel electrophoresis. The suitable fragments were selected as templates for qRT-PCR. The quality and quantity of the sample libraries were analyzed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). The qRT-PCRs were conducted using an ABI StepOne-Plus Real-Time PCR System (ABI, New York, NY, USA). Finally, the libraries were sequenced on an IlluminaHiSeq™ 2000 platform.

Quality control and alignment/mapping of clean reads
The original image data were transformed into sequence data via base calling. These raw reads were filtered using the SOAPnuke program, as previously described [20], and the remaining "clean" reads were used in the downstream bioinformatics analyses. The quality of the clean reads data was assessed using base composition and quality distribution charts. The remaining reads were aligned to the P. fortunei genome and gene sequences using BWA and Bowtie, respectively. The statistics of alignment results will be presented for each reference. During the RNA-sequencing experiment, transcripts were chemically fragmented and then sequenced. The distribution of reads on the genes was used to evaluate the randomness of the RNA fragmentation associated with the constructed libraries [29].

Gene expression analysis
Gene expression levels were determined with the RSEM (RNA-Seq by Expectation Maximization) software package. The FPKM method (fragments per kb per million fragments) was used to calculate expression levels and to identify DEGs [30]. In multiple tests, the false discovery rate (FDR) was applied to determine the threshold p-value [31]. FDR 0.001 and absolute log 2 ratio ! 1 were used as the threshold to evaluate the significance of gene expression differences.
Correlations between two samples were assessed based on the FPKM result. Ideally, the squared correlation value should be ! 0.92 according to the standard recommended by the ENCODE. The distances between expressed genes were calculated according to the Euclidean method. The sum of squared deviations algorithm was used to calculate the distance between samples so that a cluster tree could be built. The DEGs between two libraries were identified using the Audic-Claverie statistic [32]. The DEGs were functionally annotated based on GO terms and KEGG pathways.

Analysis of differentially expressed genes
The DEGs were first mapped to GO terms in the GO database, and the number of genes for each term was calculated. Then, the hypergeometric test was used to identify significantly enriched GO terms (http://www.geneontology.org/). The calculated p-values were adjusted using the Bonferroni correction [33], with a corrected p-value 0.05 used as the threshold. GO terms that fulfilled, this condition were defined as significantly enriched, and were considered to describe the main biological functions of the DEGs. The WEGO program [34] was used to functionally classify the DEGs and to assess the distribution of gene functions in P. tomentosa at a macro level. The KEGG database was used to identify enriched pathways associated with the DEGs [33].
TFs generally contain a DNA-binding domain and a trans-acting functional domain [35]. We used the hmmsearch program to search the HMMS database for domain characteristics, and to predict if a gene encoded a TF. Additionally, the TopHat program was used [36] to analyze the various types of alternative splicing (AS)-related clean reads mapped to the P. fortunei reference genome in the four samples. The AS events that were identified in both replicates were considered as stable events.

Comparison of gene expression profiles among different samples
The following pairwise comparisons of DEGs were evaluated ( PT2 and PT4 represent the diploid and autotetraploid P. tomentosa with salt-untreated, PT2S and PT4S represents the diploid and autotetraploid P. tomentosa with salt-treated (0.4% Nacl).(1) Differentially expressed genes (DEGs) were co-up-regulation between the PT2S vs. PT2 and PT4S vs. PT4 comparison to screen for gene from both diploid and autotetraploid P. tomentosa after salt treatment were strengthened. Quantitative real-time polymerase chain reaction analysis of differentially expressed genes Total RNA that was extracted from the PT2 and PT4 leaf samples as described above. Three independent biological samples of each were employed in this analysis. The SsoAdvanced™ SYBR 1 Green Supermix (Bio-Rad, Hercules, CA, US) was used for the qRT-PCRs. The amplification reactions were performed as follows: 95˚C for 3 min, followed by 40 cycles of 95˚C for 10 s and 55˚C for 30 s, and derived melting curves were obtained. The reactions were performed on a CFX96™ Real-Time System (Boi-Rad, Hercules, CA, USA) and all reactions were run in triplicate. 18S rRNA was used as an internal reference gene. The relative expression levels were calculated using the 2 -ΔΔCt method. The primers used for the qPCR reactions were listed in Table 1.

Sequencing of mRNA and alignments with the reference genes and genome
Approximately 187 million raw reads were generated for the four sequenced libraries (PT2, PT2S, PT4, and PT4S). Low-quality reads as well as reads with adapter sequences or several unknown bases were eliminated. The quality of the remaining clean reads was assessed using base composition and quality distribution charts. We observed a balanced base composition, with a T curve that was in accordance with the A curve, as well as a satisfactory base composition, high quality sequences, and base ratios that were mostly > 20 (Fig 2). The clean reads were evenly distributed on the genes as indicated by the assessment of the randomness. A total of 175 million clean reads (43, . These results suggested that high-quality sequencing data were obtained. The clean reads were mapped to the P. fortunei genes and genome with mapping means of 48.92% and 71.63% respectively ( Table 2). We observed that 8.95-17.62% of the total mapped clean reads were aligned to two or more positions, and were considered multi-position matches. Multiposition matches can be problematic during transcriptome-level studies, because they are affected by read complexity and length, Therefore, we excluded the multi-position reads from subsequent analyses to decrease the error rate. The unmapped reads may represent novel genes.

Overall analysis of gene expression
The correlation value for PT4 and PT4S was > 0.86, indicating the molecular factors responsive to salt stress partially overlapped. The correlations among the other samples are indicated in Fig 3A and S1 Table. The correlation values were close to 1 for the same samples in almost every experiment, indicating that our sequencing data were reliable and the selected samples were reasonable. The clustering of the expression profiles for the salt stress and control treatments of the four samples is presented in Fig 3B. The dendrogram revealed that PT2S and PT4S sequences tended to cluster together. We detected 4,423 DEGs between PT2S and PT2, and 4,104 DEGs between PT4S and PT4. Additionally, 2,946 and 4,400 DEGs were observed for the comparisons between PT4S and PT2S and between PT4 and PT2, respectively (Fig 4A,  S2 Table).
In the pairwise comparisons of specific sets of DEGs, we detected 1,717 DEGs for Comparison A, which may be related to P. tomentosa responses to saline conditions ( Fig 4B, S3 Table); and 160 DEGs for Comparison B, which may be relevant to the differences in the salt-induced responses of PT2 and PT4. Furthermore, Comparisons C and D revealed 589 (17.50%) and 900 (26.74%) DEGs, respectively ( Fig 4C, S3 Table), which may help to explain why PT4 is more salt-tolerant than PT2.

Analysis of differentially expressed genes using the GO and KEGG databases
The DEGs detected in Comparisons A, B, C, and D (Fig 1)    The DEGs detected in Comparisons A, B, C, and D were analyzed using the KEGG database to identify associated pathways. A KEGG enrichment analysis of DEGs from Comparisons A, B, C, and D involved 1,062, 103, 364, and 581 DEGs, accounting for 61.82%, 63.98%, 61.69%, and 64.48% of the total DEGs, respectively. The enriched pathways for the top 20 DEGs from Comparisons A, B, C, and D are provided in Fig 5. The top 20 DEGs were associated mainly with metabolic pathways (ko01100) and the biosynthesis of secondary metabolites (ko01110), as well as plant hormone signal transduction (ko04075), photosynthesis (ko00195), isoflavonoid biosynthesis (ko00943), and nitrogen metabolism (ko00910).

Differentially expressed genes involved in alternative splicing
AS affects many plant physiological processes, including responses to abiotic and biotic stresses [39]. There are seven main types of AS: Exon skipping (ES), Intron retention (IR), Alternative 5 0 splice site (A5SS), and Alternative 3 0 splice site (A3SS), Alternative first exon (AFE), Alternative last exon (ALE), and Mutually exclusive exon (MEX) (Fig 6). The AS of eukaryotic genes mainly involves the first four types. The last three types were excluded in this study because of the potential for considerable false-positive results. Among the four included AS categories, we observed that IR was the most common (73,460), followed by A3SS (33,961), A5SS (20032), and ES (13500) (Fig 7). That IR was the primary AS category was consistent with the findings of previous studies [40,41]. Of the four main AS types, A5SS, A3SS, and IR were significantly enhanced in response to saline conditions. Furthermore, IR was especially enhanced in the PT4S vs. PT2S and PT4 vs. PT2 comparisons. These results implied that IR is influenced by the ploidy level.

Differentially expressed genes involved in plant hormone signal transduction
Plant hormone signal transduction is crucial for plant responses to adverse environmental conditions. We detected 157 DEGs related to plant hormone signal transduction. In Comparison A, 90 DEGs were related to signal transductions involving several hormones, including ABA, GA, IAA, CK, ET, JA, SA, and BR, and in Comparisons B, C, and D 7, 20, and 40 DEGs were related to plant hormone signal transduction, respectively (S7 Table).

Differentially expressed genes involved in photosynthesis
Exposure to high-salt conditions impairs osmotic homeostasis as well as the cellular ionic and redox balance, and also inhibits photosynthetic activities. In this study, the "photosynthesis, light reaction", "photosynthetic membrane", "chloroplast", and "chloroplast part" GO terms were enriched in Comparisons A, B, C, and D. The GO analyses indicated that 228 DEGs were associated with "photosynthesis". The KEGG pathway analyses revealed that photosystem II P680 reaction center D1 protein (psbA), photosystem II oxygen-evolving enhancer protein 3 (psbQ), photosystem I P700 chlorophyll a apoprotein A1 (psaA), photosystem I subunit PsaN (psaN), ferredoxin (petF), and ferredoxin-NADP + reductase (petH) were significantly enriched in the four comparisons (S8 Table).

Validation of differentially expressed genes by quantitative real-time polymerase chain analysis
To further validate the reliability and reproducibility of the gene expression results from the RNA-Seq data, eight DEGs from the four libraries were selected for qRT-PCR analysis. The qRT-PCR results for the eight DEGs reflected the same trends (up-and down-regulation trends) as those acquired from the RNA-Seq data using the FPKM method (Fig 8A). The disaccord between the transcriptome and qRT-PCR results is likely because transcriptome data better reflect small changes in gene expression and better detection of low-abundant than qRT-PCR, as reported previously [42]. Overall, the qRT-PCR results support the reliability of the relative values provided by the RNA-Seq analysis. The qRT-PCR results are related to the DEG data. The expression trends of all 8 genes from RNA-Seq and qRT-PCR analyses were largely consistent, demonstrating the correlation analysis of DEGs under salt stress based on RNA-Seq and qRT-PCR data (Fig 8B). Genome-wide expression analysis of salt-stressed Paulownia tomentosa

Discussion
Salinity is a major environmental stressor that constrains tree growth in arid and semi-arid regions [43]. High-salt conditions can induce morphological and physiological changes in P. tomentosa leaves [20]. Currently, only limited information is available regarding the mechanism underlying the salt tolerance of P. tomentosa. In this study, the second pairs of leaves from the apex shoot of PT4S, PT4, PT2S, and PT2 seedlings were analyzed by RNA-Seq to elucidate the complex mechanisms that regulate the responses of diploid and autotetraploid P. tomentosa to salt-stress conditions. The transcriptome data yielded 175 million clean reads, many of which were differentially expressed. Overall, the control plants had fewer total mapped reads than the salt-treated plants (i.e., PT2 < PT2S and PT4 < PT4S) ( Table 2). Additionally, some overlapping DEGs were detected. A comparison of the control and salt-treated P. tomentosa plants of different ploidy levels decreased the number of DEGs. This increased the likelihood that the identified genes were related to plant responses and adaptations to saline conditions. "Plant hormone signal transduction" (ko07075) and "Photosynthesis" (ko00195) were considerably enriched metabolic pathways. The enriched TFs may be useful for clarifying the molecular mechanism underlying the salt tolerance of the diploid and autotetraploid P. tomentosa lines.

Transcription factors involved in salt stress responses
Transcription factors are important for the acclimation of plants to extreme environmental stresses, including high-salt conditions. In this study, members of the bHLH, MYB, NAC, GRAS, WRKY, and AP2/EREBP families were the most abundant stress-or ploidy-related TFs. Some calcium-dependent TFs are regulated by the direct binding of a Ca 2+ ion or the Ca 2+ /calmodulin (CaM) complex or by post-translation modifications mediated by Ca 2+ or the Ca 2+ /CaM complex [44]. AtNIG1 (Arabidopsis thaliana NaCl-inducible gene 1), which is a basic helix-loop-helix (bHLH) TF, was the first identified salt stress-responsive Ca 2+ -binding TF involved in salt stress signaling in plants according to a suppression subtractive hybridization analysis [45]. In the current study, we identified 62 up-regulated and 33 down-regulated bHLH genes that were induced by salt stress, which is consistent with the results of a previous study in A.thaliana [46]. The MYB super family, which includes MYB and MYB-related TFs, was the largest TF family in our study. MYB TFs have regulatory roles in plant developmental processes and defense responses [47]. A two-repeat MYB (i.e., MYB2), which contains a Ca 2 + -dependent CaM-binding domain, regulates the expression of dehydration-and salt-responsive genes in A. thaliana [48]. Additionally, the over-expression of OsMYB3R-2 in A. thaliana [49] and OsMYB48-1 in rice [50] reportedly enhanced tolerance to high-salt conditions. We detected 203 up-regulated and 165 down-regulated MYB genes in our transcriptome data. The NAC TFs, which form one of the largest plant-specific TF families, are essential for responses to various abiotic stresses. A NAC domain-containing TF was identified as a CaM-binding protein. Analyses of transcriptional regulation mediated by a CaM-binding NAC protein helped characterize the function of the Ca 2+ /CaM complex [51]. The over-expression of TaNAC2 in A. thaliana plants led to increased tolerance to salt, drought, and freezing stresses [52]. We observed that NAC expression levels were mostly up-regulated in salt-stressed P. tomentosa plants. This is consistent with the results of a previous study in cotton [53]. These DEGs may be relevant for future studies of P. tomentosa salt tolerance.

Plant hormone signal transduction pathways involved in salt stress responses
Phytohormones influence plant responses to abiotic stresses. For example, to optimize the chance of survival, phytohormone production may be altered to decrease plant growth, enabling the plant to divert additional resources to maintain responses to environmental stresses [54]. Indeed, the perception of a stress signal causes the phytohormone-related signal transduction cascades in plants to decrease to baseline levels [55,56].
ABA is a lipophilic plant hormone with key roles in signaling and adaptation to abiotic stresses, such as drought and high salinity [57]. The ABA signal transduction pathway involves the activities of PP2C (Proteinphosphatases 2C), PYR/PYL (abscisic acid receptor PYR/PYL family), SnRK2 (serine/threonine-protein kinase SRK2), and ABF (ABA responsive element binding factor). This pathway is initiated when ABA binds to the pyrabactin resistance-like/ regulatory component of ABA receptors (RCAR). Which results in the inactivation of protein phosphatase PP2C and activation of SnRK2-type kinases, ultimately inducing stomatal closure [58]. As a positive regulator of ABA signaling, SnRK2 is crucial for abiotic stress responses in plants [59]. In this study, SnRK2 expression was up-regulated in the PT4 vs. PT2 and PT4S vs. PT4 comparisons (S7 Table), similar to the results of a previous study of maize [60]. There is genetic evidence that A-type PP2Cs are negative regulators of ABA signaling in A.thaliana [61]. Under saline conditions, ABF2 expression is reportedly up-regulated in grape leaves and roots in response to increasing salt concentrations [62]. Thus, the expression of salt-responsive genes may be related to the down-regulation of PP2C and up-regulation of ABF associated with ABA signal transduction. In this study, PP2C expression (PAU005151.1, PAU008363.1, and PAU024054.1) was down-regulated in the PT4 vs. PT2 and PT4S vs. PT2S comparisons. In contrast, ABF (PAU013700.1, PAU023164.1, and PAU013700.1) expression levels increased in salt-stressed PT2S and PT4S plants.
Cytokinins (CKs) are a class of phytohormones involved in various physiological events, including responses to environmental stresses. Recently, CKs were revealed to regulate plant adaptations to salt stress during growth and development [63], and increased CK levels in seeds were observed to increase plant tolerance to salt stress [64]. The CRE1 histidine kinase (the cytokinin response 1), which contains an extracellular domain called CHASE, is a CK receptor. The CK signaling pathway involves several steps [65]. First, CRE1 is activated when CK binds to the CHASE domain and isauto-phosphorylated. Then, histidine-containing phosphotransfer proteins (AHPs) are phosphorylated by the activated CRE1,and migrate from the cytoplasm to the nucleus, where they transfer the phosphate group to (the response regulator proteins) ARRs. The expression of CRE1 (one of the receptor histidine kinases) is induced by various stresses, suggesting that the histidine kinase CK receptors are crucial for CK and stress responses in A. thaliana [66]. In this study, CRE1 expression (PAU010104.1 and PAU010665.1) was up-regulated in the PT2S vs. PT2 and PT4S vs. PT4 comparisons, which is in agreement with the aforementioned results. Furthermore, CRE1 expression (PAU021126.1) was also up-regulated in the PT4 vs. PT2 and PT4S vs. PT2S comparisons, which may be related to differences in ploidy levels. Nishiyama et al. [67] reported that among genes related to CK signaling, the expression of a type-B ARR gene was down-regulated, whereas AHP4 was up-regulated in response to salt treatments. This implies that CK signaling might be inhibited by salt stress. We observed that type-B ARR expression (PAU012990.1) levels were down-regulated in the PT2S vs. PT2 and PT4S vs. PT4 comparisons, which is consistent with the aforementioned findings. Type-B ARR expression (PAU003823.1, PAU001788.1, PAU022082.1, and PAU022609.1) was down-regulated in the PT4 vs. PT2 and PT4S vs. PT2S comparisons, which may be related to the variability in ploidy levels. The expression of AHP (PAU029211.1) was up-regulated in the PT2S vs. PT2 and PT4S vs. PT4 comparisons, which is consistent with the previously described results.
In this study, we identified 157 DEGs related to plant hormone signal transduction. The results of this study along with previously reported findings revealed that autotetraploid P. tomentosa was more salt-tolerant than diploid P. tomentosa. We observed that all the analyzed plant hormones were directly or indirectly involved in regulating plant responses to salt stress. The activities of plant genes and hormones are strongly linked, as indicated by the fact that some plant genes, which are essential for activating plant hormones and other genes, are activated by phytohormones.

Photosynthetic activities involved in salt stress responses
Salt stress decreases plant growth and productivity by disrupting physiological processes [68], especially photosynthetic activities [69]. In this study, several DEGs were involved in photosynthesis, including genes encoding the photosystem II P680 reaction center D1 protein (psbA), photosystem II oxygen-evolving enhancer protein 3 (psbQ), photosystem I P700 chlorophyll a apoprotein A1 (psaA), ferredoxin (petF), and ferredoxin-NADP + reductase (petH). Thylakoid membrane proteins were influenced by salt stress in Synechococcus sp.PCC 7942 [70]. The photosystem II D1 protein in thylakoid membranes is encoded by a gene in the plastid genome (cpDNA). This protein transforms radiant energy via the oxidation of water and reduction of plastoquinone. In Synechocystis sp., salt stress suppresses the repair of photosystem II by inhibiting the activities of the transcriptional and translational machinery [71]. Additionally, the D1 protein is sensitive to environmental stresses [72]. We observed that the expression of the gene (PAU015578.1), which encoded a D1 protein, was down-regulated in the PT4S vs. PT4 and PT2S vs. PT2 comparisons (S8 Table). This is consistent with the previous observation that D1 levels decreased because of salt stress-induced photoinhibition in Brassica juncea [73]. The psbQ peripheral protein is an important part of the photosystem II complex. It is also the most diverse extrinsic photosystem II protein in higher plants, and can maintain the integrity of photosystem II under high-salt conditions [74]. The psaQ expression level was up-regulated in the PT4 vs. PT2 and PT4S vs. PT2S comparisons, which may be closely related to the differences in ploidy. The psaA gene encodes photosystem I P700 chlorophyll a apoprotein A1, which is one of two large core subunits of photosystem I. It carries a set of cofactors required for a functional electron transport chain through photosystem I. In this study, the abundance of the psaA transcript (PAU029948.1) was up-regulated in the PT2S vs.
PT2 comparison. This is in agreement with the results of a previous study of Salicorniabigelovii Torr [75]. Electrons derived from water by the oxygen-evolving complex of photosystem II are transferred to NADP + , leading to the production of NADPH. The electrons are passed along the photosynthetic electron transport chain via the plastoquinone, cytochrome b 6 f complex, plastocyanin, photosystem I, ferredoxin, and ferredoxin-NADP + oxidoreductase [76].

Conclusions
Here we report for the first time, comparisons of the effects of salt stress on diploid and autotetraploid P. tomentosa plants using Illumina sequencing technology based on the P. fortunei genome as a reference. Overall, 85.7 and 125 million clean reads were mapped to P. fortunei gene and genome sequences, respectively. We identified 15,873 million DEGs among the control and salt-stressed P. tomentosa plants. Some DEGs were functionally associated with plant hormone signal transductions, photosynthesis, and other metabolic pathways. Notably, some TF genes, including NAC, MYB, bHLH, GRAS, WRKY, and AP2/EREBP, were responsive to salt stress and/or ploidy levels. We also identified several DEGs related to AS which were affected by salt stress or ploidy levels. Our results may be useful for future efforts to elucidate the molecular mechanisms underlying salt tolerance as well as the gene regulatory networks of P. tomentosa plants.