Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Comparative transcriptome analysis reveals K+ transporter gene contributing to salt tolerance in eggplant

BMC Plant Biology volume 19, Article number: 67 (2019) Cite this article

Abstract

Background

Soil salinization is one of the most crucial abiotic stresses that limit the growth and production of eggplant. The existing researches in eggplant were mostly focused on salt-induced morphological, biochemical and physiological changes, with only limited works centered on salt-response genes in eggplant at the transcriptomic level.

Results

Our preliminary work found that Zhusiqie (No.118) is salt-tolerant and Hongqie (No.30) is salt-sensitive. Consequently, they were re-named as ST118 and SS30, respectively. ST118 showed less damaged on growth and higher K+/Na+ ratios in leaves than SS30. Comparative-transcriptome analysis was used as a powerful approach to understand the salt-response mechanisms in the leaves and roots of SS30 and ST118. And it revealed that genotype-specific and organ-specific manners exist in eggplant in response to salt stress. Strikingly, the genotype-specific differentially expressed genes (DEGs) in ST118 were considered crucial to its higher salt-tolerance, because the expression patterns of common DEGs in the leaves/roots of the two eggplant genotypes were almost the same. Among them, some transcription factors have been reported to be in response to elevated external salinity, including the members of C2C2-CO-like, WRKY, MYB and NAC family. In addition, the AKT1, KAT1 and SOS1 were up-regulated only in the leaves of ST118. Furthermore, the complementation assays demonstrated that the salt-tolerances of both yeast and Arabidopsis akt1 mutants were enhanced by heterologous expression of SmAKT1.

Conclusion

The comparative-transcriptome analysis indicated that the salt-tolerance can be increased by higher transcript level of some genotype-specific genes. This work revealed that eggplants seem to be more inclined to absorb K+ rather than to exclude Na+ under salt stress conditions because seven K+ transporters were significantly up-regulated, while only one Na+ transporter was similarly regulated. Finally, the complementation assays of SmAKT1, which is genotype-specific up-regulated in ST118, suggest that the other TFs and K+ transport genes were worthy of future investigation for their functions in salinity tolerance.

Background

Soil salinity is one of the most important factors that limit plant growth, development, and productivity. According to the data from the FAO (Food and Agriculture Organization of the United Nations), food production should be increased by 70% in the world to meet the requirement of increasing population (http://www.fao.org/documents/card/en/c/a2128b09-361c-5468-9d93-2189cc430234/). In order to develop and utilize the salinized soil as much as possible, it is necessary to understand the salt-response mechanisms of crop plants.

To cope with salt stress, plants developed various protective mechanisms from the physiological and biochemical to the cellular and molecular level. On the molecular level, genes functioning in stress signaling, transcription regulation, ion transport and biosynthesis of specific metabolites are involved in responding to salt stress [1,2,3,4]. Transcription factors (TFs) involved in the regulation of salt-response could be activated by multiple signal transduction pathways in plants, such as the ABA-mediated signal [1, 5, 6]. Previous studies reported that the members of TF family genes were differentially expressed in response to elevated external salinity [7], including the APETALA2/ETHYLENE RESPONSE FACTOR (AP2/EREBP) [8], basic leucine zipper (bZIP) [6, 9], NAC [10, 11], basic helix–loop–helix (bHLH) [12], MYB [13,14,15] and WRKY [16, 17] gene families. In turn, these TFs could amplify the signals for gene regulation and promote the protective mechanisms in plants.

The major damage caused by excess salt was ion toxicity (mainly Na+) except water deficiency that is different from drought stress. The salt-overly-sensitive (SOS) signal transduction pathway has been described as crucial for cellular Na+ detoxification and maintaining intracellular ion homeostasis in plants [18,19,20]. However, excessive accumulation of Na+ under salt stress would be accompanied by K+ deficiency. Because of the similarity in physicochemical properties between Na+ and K+ (i.e. ionic radius and ion hydration energy), the root cells absorbed excessive Na+ instead of K+ under the saline soil [21], and ion homeostasis in plant cells could be destructed. K+ is one of the most important elements that is required by the key metabolic processes in the cytoplasm, including enzymatic reactions, protein synthesis, and ribosome functions [22]. Thus, K+/Na+ ratio is likely to be one of the key determinants of plant salt tolerance [22] and maintaining a high cytosolic K+/Na+ ratio is very important [23].

Eggplant (Solanum melongena L.) is an important greenhouse crop for out of season production and cultivated on more than 1.5Mha in the world [24]. Eggplant is considered as moderately sensitive to salinity with a very low threshold value [25, 26]. The existing researches in eggplant were focused on salt-induced morphological, biochemical and physiological changes [27,28,29,30]. However, there was limited work on salt-response genes in eggplant at the transcription level. Comparative genome and transcriptome have been extensively used as a powerful approach for discovering the genetic information involved in stress tolerance [31,32,33]. A number of transcriptomic comparisons have been done between salt-sensitive and salt-tolerant genotypes of plant species, such as Arabidopsis [34], Oryza sativa [3] and tomato [35]. Here, comparative transcriptome was used for the first time to explore the molecular mechanisms of salinity tolerance in eggplant.

In this work, the leaves and roots of two eggplant genotypes were exposed to salt-tolerant and comparative-transcriptome analysis under salt conditions. We successfully identified several TFs and ion transporters which might be crucial for the salt-tolerant eggplant genotype ST118 under salt conditions. In particular, a differentially expressed ion transporter was identified and functional verified which is potentially associated with eggplant responses and adaptability to salt conditions.

Methods

Plant material and growth conditions

Uniformly germinated eggplant seeds were selected and transplanted into growing trays with vermiculite and kept in growth chamber with 16/8 h light/dark photoperiod at 25/16 °C, respectively. About 1-months-old eggplant plants with four-true-leaves were treated with 200 mM NaCl. Roots and leaves for RNA extraction and ion content measurement were harvested at 0, 6, 12, 24, 48, 72, 168 h (7d) and 23d after stress treatments, respectively. After salt stress treatment for 23d, the phenotypic and physiological characteristics were inspected and measured. The samples at all the time points were used for ion content analysis. Based on the ion content difference between the two eggplant genotypes, samples collected at 0 h and 12 h were chosen for transcriptomic analysis.

The Arabidopsis (Arabidopsis thaliana) wild-type, mutant and transgenic plants used in this study were Columbia-0 ecotype (Col). The Arabidopsis seeds were germinated on Murashige and Skoog (MS) medium containing 0.8% (w/v) agar and 3% sucrose at 4 °C for 3 days. Then plates were incubated in a controlled-environment growth chamber. 3 days later, uniformly germinated seeds were chosen for low K+ or salt stress tests.

Ion content measurement

All the samples were dried at 105 °C for 30 min, and then kept at 75 °C for 4 days. The grinded samples were digested in 20 ml HNO3, then added 5 mL HClO4 at room temperature. After overnight digestion, HNO3 and HClO4 were removed by heating. The digested samples were diluted with ddH2O. The Na+ and K+ contents were then measured by inductively coupled plasma optical (ICP-AES, iCAP7600).

RNA extraction, library construction and illumina sequencing

Total RNA was extracted by the MiniBEST Universal RNA Extraction Kit (TaKaRa) according to the manufacturer’s instructions. The total RNA sample quality control (QC), library construction and sequencing on BGISEQ-500 was performed at Beijing Genomics Institute (BGI). The Agilent 2100 Bio analyzer (Agilent RNA 6000 Nano Kit) was used to do the total RNA sample QC, including RNA concentration, RIN value, 28S/18S and the fragment length distribution. The mRNA was enriched by magnetic beads with Oligo (dT) and then fragment the RNA and reverse transcription to double-strand cDNA (dscDNA) by N6 random primer. The synthesized cDNA was subjected to end-repair and then was 3′ adenylated. Adaptors were ligated to the ends of these 3′ adenylated cDNA fragments. The ligation products were purified and PCR amplification was performed to enrich the purified cDNA template using PCR primer. Lastly, the PCR products were denatured by heat and the single strand DNA was cyclized by splint oligo and DNA ligase. Then, the libraries were used for sequencing with the sequencing platform BGISEQ-500 (BGI), and the products were called as ‘raw reads’. All the generated raw sequencing reads were filtered to remove the low quality reads by the software SOAPnuke (BGI). After filteration, the remaining reads are called ‘Clean Reads’ and stored in FASTQ format.

Bioinformatics analysis

After QC analysis, the clean reads were assembled into Unigenes and mapped to the eggplant genome sequences (http://eggplant.kazusa.or.jp/) [36] by HISAT (Hierarchical Indexing for Spliced Alignment of Transcripts) [37]. The gene expression level was calculated with RSEM [38]. Pearson’s correlation was exploited to calculate the relevance between all samples [39]. The differentially expressed genes (DEGs) were detected with DEGseq [40], which is based on the Poisson distribution. Combining the strategies described by Y Benjamini and Y Hochberg [41] and JD Storey and R Tibshirani [42], the P-values was adjusted as Q-values. And the threshold of Q-values ≤0.001 and an absolute Log2Ratio value ≥1 among the three biological replicates were used to determine whether a gene was DEG. The sequences of DEGs were compared with the NCBI non-redundant (Nr) database to identify and annotate the obtained DEGs using Blast software [43, 44].

Gene ontology (GO) functional classification of the identified DEGs was performed using Blast2GO [45]. GO enrichment analysis of the DEGs was conducted according to the information from GO databases (http://wego.genomics.org.cn/). Then we calculate the false discovery rate (FDR) for each p-value, in general, the terms which FDR ≤ 0.01 are defined as significant enriched. As for transcription factor prediction, getorf was used to find ORF of each DEG and then ORF was aligned to TF domains (from PlntfDB) using hmmsearch [46].

Quantitative real-time PCR analysis

Total RNA was extracted by the MiniBEST Universal RNA Extraction Kit (TaKaRa). 500 ng RNA was transcribed into cDNA with the PrimeScript™ RT Master Mix (Perfect Real Time) (Takara). According to the manufacturer’s instructions of SYBR® Premix Ex Taq™ II (Tli RNaseH Plus) (Takara), qRT-PCR was performed on CFX Real Time PCR Detection System (BioRAD) using the following procedure: 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. The Smactin (Sme2.5_00072.1_g00003.1) from eggplant and AtACT2 (AT3G18780.1) were amplified in parallel as internal reference genes, respectively. The relative expression levels of the amplified products were analyzed using the comparative CT method based on CT values [47]. All primers used in this study are listed in Additional file 1.

Analysis of the protein structure

The protein sequences of AKT1 from 9 plant species were searched in the Pfam database (http://pfam.xfam.org/) and NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Then the sequences of each putative conserved domains were obtained using ClustalX (version 1.83) [48] and WebLogo 3 (http://weblogo.threeplusone.com/create.cgi). For phylogenetic analysis, ClustalX (version 1.83) and MEGA 6.0 [49] programs were used to construct neighbor-joining (NJ) tree with the following parameters: poission model, complete deletion and bootstrap (1000 replicates; random seed).

Yeast complementation

The coding sequences of SmAKT1 and AtAKT1 were constructed into pYES2.0 vector and transformed into the yeast strain R5421 (trk1△ trk2△), in which the two endogenous K+ transporter genes (TRK1, 2) were deleted. The yeast complementation assay was done as described by J Li, et al. [50]. After 5 days, all the plates were examined and photographed. Three independent experiments were performed.

Generation of Arabidopsis transgenic plants

Full-length coding sequence of SmAKT1 was constructed into the overexpression pHB vector. The construct was transformed into Arabidopsis akt1 mutant. The Arabidopsis was transformed by the floral dip method with Agrobacterium [51]. The T4 homozygous transgenic plants were used to examine the phenotype under low-K+ [50] or salt stress conditions. The expression of targeted genes in complementary plants was detected using qRT-PCR.

Salt tolerance assays of transgenic Arabidopsis

The K+ deficiency assay was done as described previously [52]. The phenotype was observed after low-K+ treatment for 7 d. For salt tolerance assays of transgenic Arabidopsis, 3-week-old wild-type, mutant and complementary plants were subjected to 200 mM NaCl treatment three times a week. The rosette leaves and roots of each Arabidopsis lines were collected for genes expression analysis after NaCl treatment for 0, 12 h and 7 days.

Chlorophyll a fluorescence

Chlorophyll a fluorescence of Arabidopsis leaves was determined with the pulse-amplitude-modulated chlorophyll fluorescence system (PAM; Heinz Walz, GmbH, Effeltrich, Germany). Plants were kept in darkness for 30 min to quantify photosystem II (PSII; Fv/Fm)) maximum efficiency using the saturation pulse method: Fv/Fm = (Fm – F0)/Fm [53, 54]. Data are the means of 6 replicates.

Data analysis

All data were presented as means with standard errors. The data were analyzed using SPSS 17.0 by one-way analysis of variance (ANOVA). Significance statistical analysis was calculated by Duncan’s Multiple Range test at significance levels of P < 0.05 and P < 0.01.

Results

Effect of salt stress on two eggplant genotypes

We investigated the salt tolerances of two eggplant genotypes, Hongqie (No.30) and Zhusiqie (No.118), four leaf-stage seedlings were irrigated with 200 mM NaCl. After 23 days, the phenotypic and physiological characteristics of No.30 and No.118 under salt conditions were compared with control. The salt-tolerance related trait values were dramatically reduced, and significant differences were detected between No.30 and No.118. As shown in Fig. 1, the plant height (PH), cross-cut length of stem, shoot dry weight (SDW) and root dry weight (RDW) were much reduced by salt stress in No.30 than No.118. However, no differences in the cross-cut width of stem and relative water content (RWC) between No.30 and No.118 were found.

Fig. 1
figure 1

Morphological trait performance of two eggplant genotypes measured under control and salt conditions. The cross-cut length and width represent stem development. Bars represent means ± SD of three biological replicates. Columns with different letters indicate significant differences at P < 0.05 (Duncan’s test)

Full size image

Furthermore, the concentration and distribution of Na+ and K+ were affected along with NaCl stress treatment for 0, 6, 12, 24, 48, 72, 168 h (7 days) and 23 days (Fig. 2 and Additional file 2). The Na+ concentrations increased significantly while K+ concentrations reduced in leaves and roots of both two eggplant genotypes (Additional file 2). As shown in Additional file 2, lower total K+ concentrations (total content covers leaves and roots) but higher total Na+ concentrations were found in No.118 than No.30, and this difference peaked after salt treatment for 12 h. However, the Na+[leaves]/Na+[roots] ratio increased less and the K+[leaves]/K+[roots] ratio increased more in No.118 than in No.30 (Fig. 2a). In addition, the K+/Na+ ratios were gradually reduced in both leaves and roots along with the salt treatment, while a higher decrease found in No.30 leaves and No.118 roots (Fig. 2b).

Fig. 2
figure 2

The distribution of Na+ and K+ and the change of K+/Na+ ratio in leaves/roots of two eggplant genotypes along with salt treatment. a The distribution of Na+ and K+ in leaves and roots. b The change of K+/Na+ ratio in leaves and roots. Three replicates were used in each time point, with three seedlings per replicate. Bars represent means ± SD of three biological replicates. Duncan’s Multiple Range test (*P < 0.05 and **P < 0.01) was used to analyze statistical significance

Full size image

These results suggested that K+ in No.118 was preferentially translocated into leaves, resulting in a higher K+/Na+ ratios. They also indicated that No.30 is more salt-sensitive than No.118. Therefore, we named the two eggplant genotypes as SS30 and ST118, respectively. Taken all together, the 0 h and 12 h time points were chosen for exploring the distribution mechanism of K+ and Na+, which might be a crucial point to explain the salt-tolerant difference between the two eggplant genotypes.

Identification of differentially expressed genes in SS30 and ST118 by RNA-seq

The leaves and roots were harvested from eggplants that have been treated with 200 mM NaCl for 0 h and 12 h, respectively. Using the BGISEQ-500 platform, an average about 24.11 M from 12 leaves samples and 23.77 M reads from 12 roots samples were generated, respectively (Additional file 3). In this project, the average mapping ratios with reference genome were 94.31% (leaves) and 89.85% (roots), the average mapping ratios with genes were 65.27% (leaves) and 58.47% (roots), and total of 50,956 (leaves) and 49,354 (roots) genes were detected. The differentially expressed genes (DEGs) were selected on the basis of DEGseq method with the following parameters: fold change ≥2 and adjusted P-value ≤0.001 [40]. A total of 5649 and 5927 DEGs were obtained in the leaves of SS30 and ST118 (Fig. 3a), while 1468 and 1202 DEGs were obtained in the roots of SS30 and ST118 (Fig. 3a). Subsequently, nine DEGs with different expression pattern were selected to validate the RNA-seq results by qRT-PCR (Additional file 4). Despite some differences, the general expression profiles were conserved between the RNA-seq and qRT-PCR data, which validates the former.

Fig. 3
figure 3

Overview and Venn diagrams of up- or down-regulated genes by salt stress in leaves and roots of both two eggplant materials at a level of ≥2-fold and adjusted P-value ≤0.001. a The total number of differentially expressed genes found in the leaves (L) and roots (R) of SS30 and ST118 in the comparison between salt-stressed (12 h) and non-stress treatments (0 h). b Four-way Venn diagram indicating that the DEGs were genotype-specific. The number of salt-up-regulated (c) and -down-regulated (d) genes found exclusively in the roots and leaves of two eggplant genotypes were analyzed

Full size image

Venn analysis showed that the DEGs identified in both SS30 and ST118 have relatively same expression patterns except 62/3723 (62 out of 3723) in leaves and 4/ 666 in roots (Additional file 5). Among these DEGs, 2664/997 and 243/419 genes were commonly up−/down-regulated in the leaves and roots of both genotypes (Additional file 5). On the other side, only 67/116 genes were commonly up−/down-regulated in both the leaves and roots of the two eggplant genotypes under salt conditions (Fig. 3c, d), indicating that the organ-specific manner adapt to salt stress, observed in rice [2] and Arabidopsis [55], also exists in eggplant.

As shown in Fig. 3b, many DEGs identified under salt conditions were genotype-specific, suggesting that the genotype-specific DEGs might be contributed to the phenotypic differences in salt-tolerance between SS30 and ST118. The Venn diagram in Fig. 3c, d showed that 1201/663 and 205/400 genes were exclusively up−/down-regulated in the leaves and roots of SS30, while 1780/376 and 203/176 genes were exclusively up−/down-regulated in the leaves and roots of ST118. Since ST118 is more salt tolerant than SS30, more attention has been paid to the genotype-specific DEGs in ST118 in the following sections.

Gene ontology analysis of DEGs

The functions of all the DEGs identified in this project were classified by the Gene Ontology (GO) assignments [56]. There were three GO categories including molecular function, biological process and cellular component in leaves and roots of both SS30 and ST118 (Additional file 6). The two largest subcategories found in the three GO categories were consistent, which were ‘metabolic process’ and ‘cellular process’ in the ‘biological process’ category, ‘cell’ and ‘cell part’ in the ‘cellular component’ category, and catalytic activity’ and ‘binding activity’ in the ‘molecular function’ category. Strikingly, the ‘response to stimulus’, ‘transporter activity’ and ‘transcription factor activity, protein binding’ were abundant GO terms.

GO terms enriched in the genotype-specific DEGs of SS30 or ST118 were identified using a threshold of P-value < 0.05. In the leaves of ST118, three GO terms were significantly enriched in the ‘molecular function’ category (Fig. 4a). Four GO terms were most abundant in the ‘molecular function’ and ‘biological process’ categories in the roots of ST118 (Fig. 4b). As for SS30, 17 GO terms were distributed in three GO categories including cellular component, molecular function and biological process in the leaves (Fig. 4c), but none GO terms were significantly enriched in the roots. Generally, the up-regulated genes enriched in both leaves and roots ST118 were much more than down-regulated genes, while opposite results showed in the molecular function and biological process categories of SS30. Compared with SS30 leaves, the genes related to ‘organic cyclic compound binding’ and ‘carbohydrate derivative binding’ activities were significantly enriched in ST118 leaves (Fig. 4a). In addition, genes related to ‘ion binding’ activity were exclusively enriched in the roots of ST118. These results suggested that genes with the binding sites for ion, inorganic or organic molecules might play important roles in response to salt stress.

Fig. 4
figure 4

GO classification of the genotype-specific DEGs in the leaves/roots of SS30 or ST118. The left y-axis shows the percentages of genes identified, and the right y-axis shows the gene number. The genes were categorized according to the annotation of GO, and the number of every category is displayed based on biological process, cellular components, and molecular functions. The enriched GO terms were identified using a threshold of P-value < 0.05

Full size image

Differentially expressed transcription factors in SS30 and ST118 caused by salt stress

10 transcription factors (TFs) were found through the analysis of the genes related to ‘heterocyclic compound binding’, ‘organic cyclic compound binding’ and ‘carbohydrate derivative binding’ activities. Considering the crucial role of TFs in response to salt stress, we highlighted the analysis on the TFs that were identified as DEGs in leaves and roots of both SS30 and ST118.

In leaves, a total of 413 TFs were identified as DEGs. Among the 413 TFs, 201 TFs were identified as DEGs in both SS30 and ST118 (named as ST-SS-L-inter), while 120 and 92 TFs were specifically identified as DEGs in SS30 (named as SS30-L-Spe) and ST118 (named ST118-L-Spe), respectively (Additional file 7a). As shown in Additional file 7a, the highest rates of induction by salt stress were observed for genes belonging to the AP2/EREBP, MYB, bHLH, WRKY, NAC, ABI3/VP1, C3H, GRAS and C2C2-Dof families. Strikingly, AP2-EREBP and MYB super-families were the largest in ST118-L-Spe and ST-SS-L-inter, while the WRKY super-family was the largest in the SS30-L-Spe. Members of these identified TFs have been reported to be associated with salt stress responses [8, 17, 57, 58]. Subsequently, the 413 TFs were searched against the Stress Responsive Transcription Factor Database (STIDB) in Arabidopsis [58] for salt-responsive genes. 43 TFs were identified as salt-response genes, including 10 MYBs, 10 NACs, 6 WRKYs, 5 AP2-EREBPs, 3 C2C2-CO-likes (COL), 3 TAZs, 2 bHLHs, 2 Tifys, C3H and G2-like (Table 1). Generally, the expression patterns in majority of TFs were same and most of them were down-regulated by salt stress (Table 1). Strikingly, in ST-SS-L-inter, Sme2.5_00556.1_g00019.1 annotated as WRKY was significantly up-regulated in ST118, but was significantly down-regulated in SS30. In addition, three TFs were slightly up-regulated in ST118 with 0.60~0.99 folds change, but significantly down-regulated in SS30 with − 3.10~ − 1.30 folds change (Table 1), including Sme2.5_03951.1_g00007.1 (annotated as MYB), Sme2.5_03886.1_g00002.1 (annotated as NAC) and Sme2.5_04464.1_g00002.1 (annotated as NAC).

Table 1 The salt-responsive TFs identified by searching against the Stress Responsive Transcription Factor Database (STIDB) in Arabidopsis in the leaves and roots of the two eggplant materials. SS_0h_L/R-Expression: the mean expression level of each TF after NaCl treatment for 0 h in the leaves/roots of SS30; log2Ratio(s2/s1): Log2(folds of mean expression in two time points); q-value: corrected p-value [41, 42]; Tair10: gene ID in Arabidosis corresponding to the TFs in eggplant
Full size table

In roots, 147 TFs were obtained including 58 ST-SS-R-inter TFs, 56 SS30-R-Spe and 33 ST118-R-Spe TFs (Additional file 7b), and the highest rates of TFs belong to AP2/EREBP, MYB and bHLH families. After searching against STIDB in Arabidopsis [58], 21 TFs were found to be salt-response genes (Table 1), including 6 MYBs, 4 NACs, 3 AP2-EREBPs, 3 COLs, 3 TAZs, bHLH and HSF. Among the 21 TFs, 11 TFs could also be identified as DEGs in leaves with the same expression pattern. Of the 11 TFs, one C2C2-CO-like family member was highly up-regulated by salt-stress with 9.1/3.5 and 7.4/4.9 folds in the leaves/roots of SS30 and ST118, respectively.

These results indicated that the basal salt-resistance mechanism was the same in eggplant varieties, but the specifically up-regulated TFs in SS118 might make a positive contribution to its salt-tolerance.

Identification of the DEGs related to ion transport in SS30 and ST118 under salt condition

K+/Na+ ratio is one of the key determinants of plant salt tolerance, and significant difference was found between SS30 and ST118. Although the ‘transporter activity’ category was enriched in both SS30 and ST118, the number and the members of genes were different. Analysis of these genes involved in the ‘transporter activity’ category showed that 43 DEGs belonged to ST-SS-L-inter, while 24 and 16 DEGs belonged to ST118-L-Spe and SS30-L-Spe, respectively.

In the ST-SS-L-inter category, five DEGs were identified as K+ transporter or K+ channel protein compared with the NCBI non-redundant (Nr) database [43] and all of them were upregulated by salt stress (Table 2). Except the five genes encoding K+ transporters or K+ channel proteins, another K+ transporter gene (AKT1) and K+ channel gene (KAT1) were specifically up-regulated by salt stress in ST118. Strikingly, the ‘salt overly sensitive’ (SOS1) gene was exclusively up-regulated in ST118 but was slightly down-regulated in SS30 (Table 2). However, no more genes related to K+ and Na+ homeostasis were found in SS30-L-spe. The specifically up-regulated expression of AKT1, KAT1 and SOS1 in ST-118 during salinity stress would be expected to stabilize the K+/Na+ ratio in leaves (Fig. 2b).

Table 2 The different expressed genes related to K+ and Na+ transport in the leaves of the two eggplant materials
Full size table

In the roots of both salt-tolerant and salt-sensitive eggplant varieties, none genes related to K+ and Na+ homeostasis was identified as DEG. Further analysis showed that the 8 ion transporter genes identified in leaves remained higher expression level in roots under both control and salt condition comparing with leaves, except for Sme2.5_00191.1_g00006.1 (KAT1) and Sme2.5_02726.1_g00002.1 (a K+ transporter gene) (Table 2). This might be the reason that salinity tolerance is more related to the fine tuning of the ion transporter genes rather than to significant up−/down-regulate these genes by salt stress in roots [3].

It is well known that SOS signaling pathway was the first demonstrated regulator in mediating Na+ extrusion in Arabidopsis and later in other plant species [23, 59,60,61]. Here, more genes closely related to K+ absorption than those related to Na+ extrusion were found to be up-regulated, indicating that K+ absorption is equally important with Na+ extrusion for maintaining K+ and Na+ homeostasis in plants under salt conditions.

Functional characterization of SmAKT1 in yeast and Arabidopsis under salt conditions

A series of studies showed that AKT1 plays an important role on resisting low-K+ stress in plants [62,63,64]. However, the function of AKT1 in eggplant under low K+-starvation and salt stress has not been report so far. The full-length amino acid sequences of the two identified AKT1s (Sme2.5_09079.1_g00001.1 and Sme2.5_00439.1_g00001.1) in eggplant together with AKT1 from the other nine plant species were aligned separately and a bootstrapped consensus neighbor-joining (NJ) tree was inferred for SmAKT1 (Additional file 8). As shown in Additional file 8, Sme2.5_00439.1_g00001.1 had the highest degree of similarity with AKT1s from the other plant species. In addition, the typically conserved domains of AKT1 were also found in an 884 amino acid polypeptide of Sme2.5_00439.1_g00001.1 (Additional file 9). Taken together, the Sme2.5_00439.1_g00001.1 could be named as SmAKT1.

Subsequently, the K+ transport activity of SmAKT1 was tested in the auxotrophic yeast mutant strain R5421 (trk1â–³, trk2â–³) [50, 65, 66] and Arabidopsis akt1 mutant [52, 62, 64], respectively. The complementation assays in yeast showed that along with the decline of K+ concentration, the growth of R5421 with empty vector was significantly depressed while both SmAKT1 and AtAKT1 could rescue the growth defect of R5421 mutant (Fig. 5a). In addition, the K+ deficiency symptoms phenotype of akt1 mutant was rescued in the two complementary Arabidopsis lines (akt1/SmAKT1–1 and akt1/SmAKT1–2), which displayed a similar phenotype with wild-type (Col) plants (Fig. 5b-d). These results suggested that SmAKT1 conferred significant K+ uptake in yeasts and Arabidopsis under low K+ concentrations condition.

Fig. 5
figure 5

Functional characterization of SmAKT1 in yeast and Arabidopsis under low K+ condition. a SmAKT1 and AKT1 complement the K+ uptake-deficient yeast mutant R5421 on AP medium containing different K+ concentrations. b Real-time PCR verification of SmAKT1 and AKT1 expression in different plant materials. c Phenotype comparison of wild-type Arabidopsis (Col), akt1 mutant and two complementary lines (akt1/SmAKT1–1 and akt1/SmAKT1–2) grown on MS and LK (100 mM K+) medium for 7 d. d Average Fv/Fm values of the whole plant. Bars represent means ± SD of three biological replicates. Columns with different letters indicate significant differences at P < 0.05 (Duncan’s test)

Full size image

In addition, the transformed yeasts were plated on AP medium containing 1, 5 or 10 mM KCl in combination with 100, 200 or 300 mM NaCl, and the yeasts expressing SmAKT1 and AtAKT1 were able to tolerate higher salt stress than the yeast with empty vector (Fig. 6a). In Arabidopsis, comparing with the WT, the growth of akt1 mutant was inhibited throughout development but was partly recovered in the two complementary lines under control condition (Fig. 6b). These results indicated that SmAKT1 was involved in responding to salt stress.

Fig. 6
figure 6

SmAKT1 is involved in response to salt stress in yeast and Arabidopsis under salt condition. a Expression of SmAKT1 and AKT1in yeast mutant strain R5421. Yeast cells were plated on AP medium containing various concentrations of Na+ (100, 200 and 300 mM) with different K+ concentrations (1, 5 and 10 mM). b Phenotype comparison of the four Arabidopsis lines after 200 mM NaCl treatment for 0 and 7 days. c Real-time quantitative PCR analysis of the expression pattern of SOS1, HKT1 and NHX in the four Arabidopsis lines treated without (control) or with 200 mM NaCl for 12 h and 7 days. Bars represent means ± SD of three biological replicates. Columns with different letters indicate significant differences at P < 0.05 (Duncan’s test)

Full size image

In order to further explore the potential molecular mechanisms underlying the above observations in Arabidopsis, the expression patterns of genes related to Na+ extrusion and transport were analyzed in the four plants under salt conditions with 200 mM NaCl (Fig. 6c). The expression patterns of SOS1, HKT1 and NHX1 in the two complementary lines were all the same with WT in both leaves and roots, while it was completely different with the akt1 mutant, except for NHX1 in leaves after 200 mM NaCl treatment for 0 h, 12 h and 7 days.

Taken together, we speculated that SmAKT1 could enhance the salt tolerance of plants not only through modulating K+ uptake, but also altering Na+ exclusion, transport and homeostasis under salt conditions.

Discussion

Control of the K+ and Na+ distribution is critical for salt-tolerance

In this study, the salt-tolerances of two eggplant genotypes were characterized. By comparison, the SS30 was more significantly affected than ST118 in the phenotypic and physiological attributes by salt stress, including PH, cross-cut length of stem, SDW, RDW and the concentration and distribution of Na+ and K+ (Fig. 1 and Additional file 2). These results were in analogy with previous studies in eggplant [30, 67, 68]. It was well known that K+/Na+ ratio in leaves is an important indicator to measure the salt-tolerance of plants [22, 23]. Here, the K+/Na+ ratio in ST118 was significantly higher in leaves but lower in roots compared with the SS30 (Fig. 2b). Although total K+ concentration was a bit lower in ST118 than in SS30, higher K+[leaves]/K+[roots] ratio were observed in ST118 than SS30. Conversely, higher total Na+ content but lower Na+[leaves]/Na+[roots] were observed in ST118 than SS30 (Fig. 2a and Additional file 2). These results suggested that ST118 preferentially translocated K+ from roots to leaves, but restricted Na+ accumulation in leaves in order to maintain a higher K+/Na+ ratios (Fig. 2b). Taken together, we speculated that the distribution mechanism of K+ and Na+ might be another key factor that determined the different salt-resistance of two eggplant genotypes.

Effect of salt stress on transcriptome changes in SS30 and ST118

Here, the comparative-transcriptome analysis between SS30 and ST118 was carried out in a way similar to previous studies in Arabidopsis [34], rice [3] and tomato [35]. Consistent with earlier studies in rice [2] and Arabidopsis [55], genotype-specific and organ-specific manners also existed in eggplant in response to salt stress (Fig. 3). Since the expression patterns of common DEGs in the leaves/roots of the two eggplant genotypes were almost same (Additional file 5), the genotype-specific DEGs in ST118 were likely responsible for the higher salt-tolerance.

The expressions of genes encoding 2 NACs, WRKY, MYB and COL transcription factors were found different between SS30 and ST118 (Table 1), which were valuable for further investigation in eggplants. Some studies have reported that the members of NAC [10, 11], MYB [13,14,15] and WRKY [16, 17] were involved in response to elevated external salinity. However, few studies on the function of COL family members in salinity tolerance have been reported so far. Although JH Min, et al. [69] reported that the AtCOL4-overexpressing plants were more tolerant to salt stress than the wild-type, most researches of the COL genes family focused on exploring its function on the flowering time of plants, such as OsCOL10 [70], OsCOL9 [71] and GhCOL1 [72]. In addition, previous studies reported that BTB/TAZ played an essential role during gametogenesis, and probably throughout plant development [73]. Recently, Q Zhao, et al. [74] reported that MdBT1/2 (a BTB/TAZ protein) interact with MdCUL3 to bridge the formation of the MdBTsMdCUL3 complex, which negatively modulates the degradation of the MdbHLH104 protein in response to changes in Fe status to maintain iron homeostasis in plants. And V Araus, et al. [75] reported that BT2 was the most central and connected gene in the nitrogen use efficiency (NUE) network in Arabidopsis and rice. Taken together, we thought that the 6 TFs are good candidates for further investigation of their role in salinity tolerance.

Candidate genes associated with K+ and Na+ homeostasis

Maintaining ion homeostasis is one of the key determinants for the plants survival under salt stress. The finding in this work that the Na+[leaves] /Na+[roots] increased less in ST118 than in SS30 along with salt treatment, indicating that ST118 may possess a mechanism to restrict the accumulation of Na+ in the leaves. The ‘salt overly sensitive’ (SOS) signaling pathway, including SOS1, SOS2 and SOS3 genes, has been proven to be important for plant salt tolerance [76, 77]. Among them, SOS1 was well known to be expressed in root epidermal cells and xylem parenchyma cells and was involved in extruding Na+ into the external medium and loading Na+ into the xylem for long-distance transport to leaves [61, 78, 79]. However, the SOS1 were expressed constitutively at higher levels in the roots of both eggplant genotypes. Strikingly, SOS1 was significantly up-regulated in the leaves of ST118 but was slightly down-regulated in SS30. Previous studies have been reported that SOS1 is also expressed in the xylem parenchyma in leaves but where its function is unclear so far. JK Zhu [18] speculated that the function of SOS1 in leaves may function to extrude Na+ from the xylem parenchyma cells into the apoplastic space of mesophyll cells.

Except SOS1, seven genes encoding K+ transporters or K+ channel proteins were identified in leaves to be up-regulated in response to salt stress (Table 2). It is worth noting that the genes encoding KAT1 and AKT1 were significantly up-regulated only in ST118. These results could partially explain the higher K+ [leaves] /K+ [roots] ratio in ST118 than in SS30 under salt stress. Similar with SOS1, AKT1 were expressed constitutively at higher levels in the roots of both two eggplant genotypes (Table 2). AKT1 was the first inward-rectifying K+ channel identified in Arabidopsis by functional complementation of yeast mutant strains defective in K+ transport system [80]. Moreover, a model of K+ uptake regulated by AKT1 in Arabidopsis and Oryza sativa under low-K+ conditions was proposed [50, 52]. In addition, previous studies showed that the osmotic- and drought-tolerance of rice could be enhanced by overexpression of OsAKT1 [81]. In fact, K+ deficiency would be accompanied by excessive accumulation of Na+ under salt stress. However, extensive researches were directed to the genes related to the influx, extrusion and accumulation of Na+ to improve K+/Na+ ratio in plants. Relatively limited studies focused on investigation of the AKT1 roles in maintaining K+ and Na+ homeostasis in plants under salt stress, especially in eggplants [82,83,84].

SmAKT1 is not only involved in modulating K+ uptake, but also in altering Na+ exclusion, transportation and homeostasis in Arabidopsis under salt conditions

In this study, more genes related to K+ uptake were identified as DEGs than those related to Na+ regulation (Table 2). Subsequently, the complementation assays in both yeast and Arabidopsis akt1 mutants demonstrated that SmAKT1 was involved in response to both low-K+ condition (Fig. 5) and salt conditions (Fig. 6a, b). Given the phenotype of K+ concentration and distribution in ST118 under salt stress, we speculated that SmAKT1 not only mediates K+ uptake in roots, but is also essential for maintaining long-distance transport and homeostasis of K+ in eggplants, which is similar with Arabidopsis [52] and Z. xanthoxylum [83]. In addition, the expression patterns of SOS1, HKT1 and NHX1, known as Na+ uptake and transport genes, were significantly changed in the Arabidopsis akt1 mutants and recovered in the complementary lines, when compared with the wild type under salt stress.

As described above, the functions of SOS1 in roots were Na+ extrusion and Na+ upload into the xylem [18, 78, 79]. However, the coordination mechanism of these two roles is not well understood. The other important Na+ transporter is HKT1, which acts in the retrieval of Na+ from the xylem to restrict the Na+ amount in the transpiration stream in roots [85] and uploading Na+ into the phloem for recirculation back to roots from the leaves [86]. In this work, the transcription level of SOS1 was significantly increased in roots but no change was observed in leaves of akt1 mutants after short-term salt treatment (12 h) (Fig. 6c). However, the transcription level of HKT1 was significantly decreased in roots but was increased in leaves at the same time (Fig. 6c). Taken together, we speculated that the akt1 mutant transported Na+ into leaves by SOS1, while the leaves restrict Na+ accumulation in leaves by the function of HKT1. And these opposing works by two different genes might be an important reason for its intolerance to salt stress. In contrast to the akt1 mutant, the wild type and the two complementary lines might have developed a mechanism to avoid Na+ from entering into leaves and to transfer Na+ into the apoplastic space of mesophyll cells as soon as possible. In addition, after being exposed to prolonged salt stress (7d), the wild type and the two complementary lines unload the Na+ from xylem by upregulating the expression of HKT1 and extrude it to soil solution by upregulating the expression of SOS1.

As for NHX1, it seems to be not closely to AKT1, and the function of NHX1 in plant leaves has been well studied while it was only partly understood in roots, which could transport the excessive Na+ to vacuole. Here, the NHX1 was down-regulated in roots of the three tolerant Arabidopsis lines at 12 h but up-regulated at 7 days (Fig. 6c). It could be explained as that the vacuole was the ultimate storage space for additional Na+.

Taken together, our results suggested that SmAKT1 is an important determinant for maintaining K+ and Na+ homeostasis in eggplant under salt stress, since it not only mediates K+ uptake, but also modulates Na+ uptake and transport systems.

Conclusion

In order to grow on saline soils, plants developed coordinated physiological traits throughout the lifecycle, among which the K+ and Na+ homeostasis is a key determinant to evaluate salt-tolerance. Here, comparative analysis of transcript levels in response to salt stress between salt-sensitive and salt-tolerant eggplant genotypes provided insights into key candidate genes related to salinity tolerance. The transcriptomic differences between SS30 and ST118 indicated the diversity of approaches to resist the challenge of salt stress. Further, the differently expressed TFs and ion transport genes were selectively analyzed, and the complementation assays demonstrated that SmAKT1 is an important regulator under salt conditions. Objectively, it also suggested that the other TFs and K+ transport genes were also worth further investigation for their functions in salinity tolerance. These data provides a foundation for elucidating the molecular networks underlying salt tolerance in eggplants.

Abbreviations

FAO:

Food and Agriculture Organization of the United Nations

SS30-L/R-Spe:

The DEGs specifically identified in the leaves/roots of SS30

ST118-L/R-Spe:

The DEGs specifically identified in the leaves/roots of ST118

STIDB:

The Stress Responsive Transcription Factor Database

ST-SS-L/R-inter:

The intersection of DEGs identified in the leaves/roots of both SS30 and ST118

References

  1. Kumar K, Kumar M, Kim S-R, Ryu H, Cho Y-G. Insights into genomics of salt stress response in rice. Rice. 2013;6:27.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Wang WS, Zhao XQ, Li M, Huang LY, Xu JL, Zhang F, Cui YR, Fu BY, Li ZK. Complex molecular mechanisms underlying seedling salt tolerance in rice revealed by comparative transcriptome and metabolomic profiling. J Exp Bot. 2016;67:405–19.

    Article  CAS  PubMed  Google Scholar 

  3. Cotsaftis O, Plett D, Johnson AAT, Walia H, Wilson C, Ismail AM, Close TJ, Tester M, Baumann U. Root-specific transcript profiling of contrasting Rice genotypes in response to salinity stress. Mol Plant. 2011;4:25–41.

    Article  CAS  PubMed  Google Scholar 

  4. Walia H, Wilson C, Zeng L, Ismail AM, Condamine P, Close TJ. Genome-wide transcriptional analysis of salinity stressed japonica and indica rice genotypes during panicle initiation stage. Plant Mol Biol. 2007;63:609–23.

    Article  CAS  PubMed  Google Scholar 

  5. Singh K, Foley RC, Onate-Sanchez L. Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002;5:430–6.

    Article  CAS  PubMed  Google Scholar 

  6. Todaka D, Nakashima K, Shinozaki K, Yamaguchi-Shinozaki K. Toward understanding transcriptional regulatory networks in abiotic stress responses and tolerance in rice. Rice (N Y). 2012;5:6.

    Article  Google Scholar 

  7. Golldack D, Luking I, Yang O. Plant tolerance to drought and salinity: stress regulating transcription factors and their functional significance in the cellular transcriptional network. Plant Cell Rep. 2011;30:1383–91.

    Article  CAS  PubMed  Google Scholar 

  8. Yang R, Liu J, Lin Z, Sun W, Wu Z, Hu H, Zhang Y. ERF transcription factors involved in salt response in tomato. Plant Growth Regul. 2018;84:1–10.

    Article  CAS  Google Scholar 

  9. Banerjee A, Roychoudhury A. Abscisic-acid-dependent basic leucine zipper (bZIP) transcription factors in plant abiotic stress. Protoplasma. 2017;254:3–16.

    Article  CAS  PubMed  Google Scholar 

  10. Takasaki H, Maruyama K, Kidokoro S, Ito Y, Fujita Y, Shinozaki K, Yamaguchi-Shinozaki K, Nakashima K. The abiotic stress-responsive NAC-type transcription factor OsNAC5 regulates stress-inducible genes and stress tolerance in rice. Mol Gen Genomics. 2010;284:173–83.

    Article  CAS  Google Scholar 

  11. Song SY, Chen Y, Chen J, Dai XY, Zhang WH. Physiological mechanisms underlying OsNAC5-dependent tolerance of rice plants to abiotic stress. Planta. 2011;234:331–45.

    Article  CAS  PubMed  Google Scholar 

  12. Mao K, Dong Q, Li C, Liu C, Ma F. Genome wide identification and characterization of apple bHLH transcription factors and expression analysis in response to drought and salt stress. Front Plant Sci. 2017;8:480.

    PubMed  PubMed Central  Google Scholar 

  13. Li C, Ng CK-Y, Fan L-M. MYB transcription factors, active players in abiotic stress signaling. Environ Exp Bot. 2015;114:80–91.

    Article  CAS  Google Scholar 

  14. Dai X, Xu Y, Ma Q, Xu W, Wang T, Xue Y, Chong K. Overexpression of an R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic Arabidopsis. Plant Physiol. 2007;143:1739–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Yang A, Dai X, Zhang W-H. A R2R3-type MYB gene, OsMYB2, is involved in salt, cold, and dehydration tolerance in rice. J Exp Bot. 2012;63:2541–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wu J, Chen J, Wang L, Wang S. Genome-wide investigation of WRKY transcription factors involved in terminal drought stress response in common bean. Front Plant Sci. 2017;8:380.

    PubMed  PubMed Central  Google Scholar 

  17. Jiang J, Ma S, Ye N, Jiang M, Cao J, Zhang J. WRKY transcription factors in plant responses to stresses. J Integr Plant Biol. 2017;59:86–101.

    Article  CAS  PubMed  Google Scholar 

  18. Zhu JK. Abiotic stress signaling and responses in plants. Cell. 2016;167:313–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Liu JP, Ishitani M, Halfter U, Kim CS, Zhu JK. The Arabidopsis thaliana SOS2 gene encodes a protein kinase that is required for salt tolerance. Proc Natl Acad Sci U S A. 2000;97:3730–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Quintero FJ, Ohta M, Shi HZ, Zhu JK, Pardo JM. Reconstitution in yeast of the Arabidopsis SOS signaling pathway for Na+ homeostasis. Proc Natl Acad Sci U S A. 2002;99:9061–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Rozhdestvenskii VI, Vil’iams MV, Tsvetkova IV, Lebedeva EV, Alekhina TP. Control of mineral nutrition of higher plants in biological life support systems. Kosm Biol Aviakosm Med. 1975;9:30–5.

    CAS  PubMed  Google Scholar 

  22. Shabala S, Cuin TA. Potassium transport and plant salt tolerance. Physiol Plant. 2008;133:651–69.

    Article  CAS  PubMed  Google Scholar 

  23. Zhu JK. Regulation of ion homeostasis under salt stress. Curr Opin Plant Biol. 2003;6:441–5.

    Article  CAS  PubMed  Google Scholar 

  24. Kantharajah A, Golegaonkar P. Somatic embryogenesis in eggplant. Sci Hortic. 2004;99:107–17.

    Article  Google Scholar 

  25. Jain RK, Dhawan RS, Sharma DR, Chowdhury JB. Salt-tolerance and proline accumulation: a comparative study in salt-tolerant and wild type cultured cells of eggplant. Plant Cell Rep. 1987;6:382–4.

    Article  CAS  PubMed  Google Scholar 

  26. Heuer B, Meiri A, Shalevet J. Salt tolerance of eggplant. Plant Soil. 1986;95:9–13.

    Article  CAS  Google Scholar 

  27. Akinci IE, Akinci S, Yilmaz K, Dikici H. Response of eggplant varieties (Solanum melongena) to salinity in germination and seedling stages. N Z J Crop Hortic Sci. 2004;32:193–200.

    Article  Google Scholar 

  28. Unlukara A, Kurunc A, Kesmez GD, Yurtseven E, Suarez DL. Effects of salinity on eggplant (Solanum Melongena L.) growth and evapotranspiration. Irrig Drain. 2010;59:203–14.

    Google Scholar 

  29. Yasar F, Ellialtioglu S, Kusvuran S. Ion and lipid peroxide content in sensitive and tolerant eggplant callus cultured under salt stress. Eur J Hortic Sci. 2006;71:169–72.

    CAS  Google Scholar 

  30. Abbas W, Ashraf M, Akram NA. Alleviation of salt-induced adverse effects in eggplant (Solanum melongena L.) by glycinebetaine and sugarbeet extracts. Sci Hortic. 2010;125:188–95.

    Article  CAS  Google Scholar 

  31. Roy NC, Altermann E, Park ZA, McNabb WC. A comparison of analog and next-generation transcriptomic tools for mammalian studies. Brief Funct Genomics. 2011;10:135–50.

    Article  CAS  PubMed  Google Scholar 

  32. Ward JA, Ponnala L, Weber CA. Strategies for transcriptome analysis in nonmodel plants. Am J Bot. 2012;99:267–76.

    Article  CAS  PubMed  Google Scholar 

  33. Xiao N, Gao Y, Qian H, Gao Q, Wu Y, Zhang D, Wang Z, Zhang X, Yu L, Li Y. Identification of genes related to cold tolerance and a functional allele that confers cold tolerance. Plant Physiol. 2018;177:1108–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wang YP, Yang L, Zheng ZM, Grumet R, Loescher W, Zhu JK, Yang PF, Hu YL, Chan ZL. Transcriptomic and physiological variations of three Arabidopsis ecotypes in response to salt stress. PLoS One. 2013;8:e69036.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Sun W, Xu XN, Zhu HS, Liu AH, Liu L, Li JM, Hua XJ. Comparative transcriptomic profiling of a salt-tolerant wild tomato species and a salt-sensitive tomato cultivar. Plant Cell Physiol. 2010;51:997–1006.

    Article  CAS  PubMed  Google Scholar 

  36. Hirakawa H, Shirasawa K, Miyatake K, Nunome T, Negoro S, Ohyama A, Yamaguchi H, Sato S, Isobe S, Tabata S, et al. Draft genome sequence of eggplant (Solanum melongena L.): the representative solanum species indigenous to the old world. DNA Res. 2014;21:649–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. Bmc Bioinformatics. 2011;12:323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Sedgwick P. Pearson’s correlation coefficient. Br Med J. 2012;344:e4483.

    Article  Google Scholar 

  40. Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26:136–8.

    Article  PubMed  CAS  Google Scholar 

  41. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57:289–300.

  42. Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100:9440–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5.

    Article  CAS  PubMed  Google Scholar 

  44. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ, Gapped BLAST. PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.

    Article  CAS  PubMed  Google Scholar 

  46. Mistry J, Finn RD, Eddy SR, Bateman A, Punta M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013;41:e121.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li J, Ren L, Gao Z, Jiang M, Liu Y, Zhou L, He Y, Chen H. Combined transcriptomic and proteomic analysis constructs a new model for light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.). Plant Cell Environ. 2017;40:3069–87.

    Article  CAS  PubMed  Google Scholar 

  48. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Li J, Long Y, Qi GN, Li J, Xu ZJ, Wu WH, Wang Y. The Os-AKT1 channel is critical for K+ uptake in Rice roots and is modulated by the Rice CBL1-CIPK23 complex. Plant Cell. 2014;26:3387–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Clough SJ, Bent AF. Floral dip: a simplified method for agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998;16:735–43.

    Article  CAS  PubMed  Google Scholar 

  52. Xu J, Li HD, Chen LQ, Wang Y, Liu LL, He L, Wu WH. A protein kinase, interacting with two calcineurin B-like proteins, regulates K+ transporter AKT1 in Arabidopsis. Cell. 2006;125:1347–60.

    Article  CAS  PubMed  Google Scholar 

  53. Genty B, Briantais J-M, Baker NR. The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence. Biochim Biophys Acta (BBA)-Gen Subj. 1989;990:87–92.

    Article  CAS  Google Scholar 

  54. Kitajima M, Butler WL. Quenching of chlorophyll fluorescence and primary photochemistry in chloroplasts by dibromothymoquinone. Biochim Biophys Acta. 1975;376:105–15.

    Article  CAS  PubMed  Google Scholar 

  55. Begara-Morales JC, Sanchez-Calvo B, Luque F, Leyva-Perez MO, Leterrier M, Corpas FJ, Barroso JB. Differential transcriptomic analysis by RNA-Seq of GSNO-responsive genes between Arabidopsis roots and leaves. Plant Cell Physiol. 2014;55:1080–95.

    Article  CAS  PubMed  Google Scholar 

  56. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Kim HJ, Nam HG, Lim PO. Regulatory network of NAC transcription factors in leaf senescence. Curr Opin Plant Biol. 2016;33:48–56.

    Article  PubMed  CAS  Google Scholar 

  58. Naika M, Shameer K, Mathew OK, Gowda R, Sowdhamini R. STIFDB2: an updated version of plant stress-responsive transcription factor database with additional stress signals, stress-responsive transcription factor binding sites and stress-responsive genes in Arabidopsis and rice. Plant Cell Physiol. 2013;54:e8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Shi H, Ishitani M, Kim C, Zhu JK. The Arabidopsis thaliana salt tolerance gene SOS1 encodes a putative Na+/H+ antiporter. Proc Natl Acad Sci U S A. 2000;97:6896–901.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Ji H, Pardo JM, Batelli G, Van Oosten MJ, Bressan RA, Li X. The salt overly sensitive (SOS) pathway: established and emerging roles. Mol Plant. 2013;6:275–86.

    Article  CAS  PubMed  Google Scholar 

  61. Olias R, Eljakaoui Z, Li J, De Morales PA, Marin-Manzano MC, Pardo JM, Belver A. The plasma membrane Na+/H+ antiporter SOS1 is essential for salt tolerance in tomato and affects the partitioning of Na+ between plant organs. Plant Cell Environ. 2009;32:904–16.

    Article  CAS  PubMed  Google Scholar 

  62. Hirsch RE, Lewis BD, Spalding EP, Sussman MR. A role for the AKT1 potassium channel in plant nutrition. Science. 1998;280:918–21.

    Article  CAS  PubMed  Google Scholar 

  63. Gierth M, Maser P, Schroeder JI. The potassium transporter AtHAK5 functions in K(+) deprivation-induced high-affinity K(+) uptake and AKT1 K(+) channel contribution to K(+) uptake kinetics in Arabidopsis roots. Plant Physiol. 2005;137:1105–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Spalding EP, Hirsch RE, Lewis DR, Qi Z, Sussman MR, Lewis BD. Potassium uptake supporting plant growth in the absence of AKT1 channel activity - inhibition by ammonium and stimulation by sodium. J Gen Physiol. 1999;113:909–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Gaber RF, Styles CA, Fink GR. TRK1 encodes a plasma membrane protein required for high-affinity potassium transport in Saccharomyces cerevisiae. Mol Cell Biol. 1988;8:2848–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Nakamura RL, Anderson JA, Gaber RF. Determination of key structural requirements of a K+ channel pore. J Biol Chem. 1997;272:1011–8.

    Article  CAS  PubMed  Google Scholar 

  67. Shahbaz M, Mushtaq Z, Andaz F, Masood A. Does proline application ameliorate adverse effects of salt stress on growth, ions and photosynthetic ability of eggplant (Solanum melongena L.)? Sci Hortic. 2013;164:507–11.

    Article  CAS  Google Scholar 

  68. Chartzoulakis KS, Loupassaki MH. Effects of NaCl salinity on germination, growth, gas exchange and yield of greenhouse eggplant. Agric Water Manag. 1997;32:215–25.

    Article  Google Scholar 

  69. Min JH, Chung JS, Lee KH, Kim CS. The CONSTANS-like 4 transcription factor, AtCOL4, positively regulates abiotic stress tolerance through an abscisic acid-dependent manner in Arabidopsis. J Integr Plant Biol. 2015;57:313–24.

    Article  CAS  PubMed  Google Scholar 

  70. Tan JJ, Jin MN, Wang JC, Wu FQ, Sheng PK, Cheng ZJ, Wang JL, Zheng XM, Chen LP, Wang M, et al. OsCOL10, a CONSTANS-like gene, functions as a flowering time repressor downstream of Ghd7 in Rice. Plant Cell Physiol. 2016;57:798–812.

    Article  CAS  PubMed  Google Scholar 

  71. Liu H, Gu FW, Dong SY, Liu W, Wang H, Chen ZQ, Wang JF. CONSTANS-like 9 (COL9) delays the flowering time in Oryza sativa by repressing the Ehd1 pathway. Biochem Biophys Res Commun. 2016;479:173–8.

    Article  CAS  PubMed  Google Scholar 

  72. Cai DR, Liu H, Sang N, Huang XZ. Identification and characterization of CONSTANS-like (COL) gene family in upland cotton (Gossypium hirsutum L.). PLoS One. 2017;12:e0179038.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  73. Robert HS, Quint A, Brand D, Vivian-Smith A, Offringa R. BTB and TAZ domain scaffold proteins perform a crucial function in Arabidopsis development. Plant J. 2009;58:109–21.

    Article  CAS  PubMed  Google Scholar 

  74. Zhao Q, Ren YR, Wang QJ, Wang XF, You CX, Hao YJ. Ubiquitination-related MdBT scaffold proteins target a bHLH transcription factor for Iron homeostasis. Plant Physiol. 2016;172:1973–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Araus V, Vidal EA, Puelma T, Alamos S, Mieulet D, Guiderdoni E, Gutiérrez RA. Members of BTB gene family regulate negatively nitrate uptake and nitrogen use efficiency in Arabidopsis thaliana and Oryza sativa. Plant Physiol. 2016;171:1523–1532.

  76. Quan R, Lin H, Mendoza I, Zhang Y, Cao W, Yang Y, Shang M, Chen S, Pardo JM, Guo Y. SCABP8/CBL10, a putative calcium sensor, interacts with the protein kinase SOS2 to protect Arabidopsis shoots from salt stress. Plant Cell. 2007;19:1415–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53:247–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Zhu M, Shabala L, Cuin TA, Huang X, Zhou M, Munns R, Shabala S. Nax loci affect SOS1-like Na+/H+ exchanger expression and activity in wheat. J Exp Bot. 2016;67:835–44.

    Article  CAS  PubMed  Google Scholar 

  79. Shi H, Quintero FJ, Pardo JM, Zhu JK. The putative plasma membrane Na(+)/H(+) antiporter SOS1 controls long-distance Na(+) transport in plants. Plant Cell. 2002;14:465–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Sentenac H, Bonneaud N, Minet M, Lacroute F, Salmon JM, Gaymard F, Grignon C. Cloning and expression in yeast of a plant potassium ion transport system. Science. 1992;256:663–5.

    Article  CAS  PubMed  Google Scholar 

  81. Ahmad I, Mian A, Maathuis FJM. Overexpression of the rice AKT1 potassium channel affects potassium nutrition and rice drought tolerance. J Exp Bot. 2016;67:2689–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Duan HR, Ma Q, Zhang JL, Hu J, Bao AK, Wei L, Wang Q, Luan S, Wang SM. The inward-rectifying K+ channel SsAKT1 is a candidate involved in K+ uptake in the halophyte Suaeda salsa under saline condition. Plant Soil. 2015;395:173–87.

    Article  CAS  Google Scholar 

  83. Ma Q, Hu J, Zhou XR, Yuan HJ, Kumar T, Luan S, Wang SM. ZxAKT1 is essential for K+ uptake and K+/Na+ homeostasis in the succulent xerophyte Zygophyllum xanthoxylum. Plant J. 2017;90:48–60.

    Article  CAS  PubMed  Google Scholar 

  84. Ardie SW, Liu SK, Takano T. Expression of the AKT1-type K+ channel gene from Puccinellia tenuiflora, PutAKT1, enhances salt tolerance in Arabidopsis. Plant Cell Rep. 2010;29:865–74.

    Article  CAS  PubMed  Google Scholar 

  85. Sunarpi HT, Motoda J, Kubo M, Yang H, Yoda K, Horie R, Chan WY, Leung HY, Hattori K, et al. Enhanced salt tolerance mediated by AtHKT1 transporter-induced Na unloading from xylem vessels to xylem parenchyma cells. Plant J. 2005;44:928–38.

    Article  CAS  PubMed  Google Scholar 

  86. Berthomieu P, Conejero G, Nublat A, Brackenbury WJ, Lambert C, Savio C, Uozumi N, Oiki S, Yamada K, Cellier F, et al. Functional analysis of AtHKT1 in Arabidopsis shows that Na+ recirculation by the phloem is crucial for salt tolerance. EMBO J. 2003;22:2004–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We thank Richard Gaber (Northwestern University) and Weihua Wu (China Agricultural University) for providing yeast strains R5421 that were used in the yeast complementation experiments. We thank Weihua Wu (China Agricultural University) for providing the Arabidopsis akt1 mutant.

Funding

The financial of this research was supported by the National Natural Science Foundation of China (31471870) and the Science Project Supported by Science and Technology Commission of Shanghai Municipality (18391900500). Funding body had no role in the design of the study and collection, analysis, interpretation of data and in writing the manuscript.

Availability of data and materials

The raw sequencing data were deposited in the NCBI Short Read Archive (SRA) database (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/) under the accession number SRP151507. The other supporting data were included as additional files.

Author information

Authors and Affiliations

  1. School of Agriculture and Biology, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai, 200240, China

    Jing Li, Zhen Gao, Lu Zhou, Linzhi Li, Junhao Zhang, Yang Liu & Huoying Chen

Authors
  1. Jing Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zhen Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Lu Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Linzhi Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Junhao Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yang Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Huoying Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

H-YC, YL, ZG and JL conceived and designed the experiments; JL and ZG performed experiments and analyzed data; LZ, L-ZL and J-HZ assisted with the experiments; JL wrote the manuscript, and H-YC and YL revised. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Yang Liu or Huoying Chen.

Ethics declarations

Ethics approval and consent to participate

Not applicable, as this study did not involve human or animal subjects, and the seeds of two eggplant cultivars were stored in School of Agriculture and Biology, Shanghai Jiao Tong University. The seeds of Arabidopsis akt1 mutant for were provided by Prof. Weihua Wu (China Agricultural University).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. List of primers sequences used in this study. (DOCX 32 kb)

Additional file 2:

Figure S1. The K+ (a) and Na+ (b) content in leaves and roots of two eggplant genotypes along with 200 mM NaCl treatment. DW represents dry weight. Three replicates were used in each time point, with three seedlings per replicate. Bars represent means ± SD of three biological replicates. Duncan’s Multiple Range test (*P < 0.05 and **P < 0.01) was used to analyze statistical significance. (DOCX 391 kb)

Additional file 3:

Table S2. Summary statistics of sequencing and assembly. Tissue: The tissue of eggplant seedling; Samples: Sample names; Total Clean Reads(Mb): The reads amount after filtering, Unit: Mb; Clean Reads Ratio(%): The ratio of the amount of filtered clean reads; Total Mapping Ratio: The percentage of mapped reads; Uniquely Mapping Ratio: The percentage of uniquely mapped reads (%); Expressed Gene No.: The amount of expressed genes; SS represents salt sensitive eggplant SS30; ST represents salt tolerant eggplant ST118; 0 h and 12 h represent the time after NaCl treatment; L: leaves; R: Roots. (DOCX 40 kb)

Additional file 4:

Figure S2. Validation of RNA-seq data in leaves (a) and roots (b) using qRT-PCR. (DOCX 1923 kb)

Additional file 5:

Figure S3. Four-way Venn diagram indicating the number of salt-up-regulated and -down-regulated genes found exclusively in the leaves (a) and roots (b) of two eggplant genotypes in the comparison between salt-stressed and non-stress treatments. (DOCX 468 kb)

Additional file 6:

Figure S4. GO classification of up- and down-regulated genes in leaves or roots of SS30 or ST118. (DOCX 1470 kb)

Additional file 7:

Figure S5. Overview the salt-up- or down-regulated TFs in the leaves and roots of both two eggplant genotypes at a level of ≥2-fold and adjusted P-value ≤0.001. (DOCX 3480 kb)

Additional file 8:

Figure S6. Phylogenetic relationships of the two SmAKT1s with AKT1 from other species. Protein sequences of AKT1 were analyzed using MEGA7.0 and the Neighbor-Joining method with 1000 bootstrap replicates. (DOCX 258 kb)

Additional file 9:

Figure S7. The conserved domains in across AKT1 proteins. The overall height of each stack indicates the conservation of the sequence at that position, whereas the height of letters within each stack represents the relative frequency of the corresponding amino acid. (DOCX 1146 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, J., Gao, Z., Zhou, L. et al. Comparative transcriptome analysis reveals K+ transporter gene contributing to salt tolerance in eggplant. BMC Plant Biol 19, 67 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-1663-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-1663-8

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us