Skip to main content

Comparative transcriptome analysis reveals gene expression differences between two peach cultivars under saline-alkaline stress



Saline-alkaline stress is a major abiotic stress that is harmful to plant growth worldwide. Two peach cultivars (GF677 and Maotao) display distinct phenotypes under saline-alkaline stress. The molecular mechanism explaining the differences between the two cultivars is still unclear.


In the present study, we systematically analysed the changes in GF677 and Maotao leaves upon saline-alkaline stress by using cytological and biochemical technologies as well as comparative transcriptome analysis. Transmission electron microscopy (TEM) observations showed that the structure of granum was dispersive in Maotao chloroplasts. The biochemical analysis revealed that POD activity and the contents of chlorophyll a and chlorophyll b, as well as iron, were notably decreased in Maotao. Comparative transcriptome analysis detected 881 genes with differential expression (including 294 upregulated and 587 downregulated) under the criteria of |log2 Ratio| ≥ 1 and FDR ≤0.01. Gene ontology (GO) analysis showed that all differentially expressed genes (DEGs) were grouped into 30 groups. MapMan annotation of DEGs showed that photosynthesis, antioxidation, ion metabolism, and WRKY TF were activated in GF677, while cell wall degradation, secondary metabolism, starch degradation, MYB TF, and bHLH TF were activated in Maotao. Several iron and stress-related TFs (ppa024966m, ppa010295m, ppa0271826m, ppa002645m, ppa010846m, ppa009439m, ppa008846m, and ppa007708m) were further discussed from a functional perspective based on the phylogenetic tree integration of other species homologues.


According to the cytological and molecular differences between the two cultivars, we suggest that the integrity of chloroplast structure and the activation of photosynthesis as well as stress-related genes are crucial for saline-alkaline resistance in GF677. The results presented in this report provide a theoretical basis for cloning saline-alkaline tolerance genes and molecular breeding to improve saline-alkaline tolerance in peach.


Saline-alkaline stress adversely affects plant growth and development. Based on the Food and Agriculture Organization of the United Nations (FAO) research, at least 242 million hectares of saline-alkaline soils are distributed within the Eurasian Region [1].

High salt concentrations and high pH levels of saline-alkaline soil generate reactive oxygen species (ROS), repress photosynthesis, energy production, and lipid metabolism, and damage the plant cell membrane and intracellular components [2]. Upon saline-alkaline stress, a large cluster of gene expression is reprogrammed. Revealing the molecular mechanisms is helpful for saline-alkaline resistant plant breeding. In recent decades, several key genes such as OsLOL5, Gshdz4, and SsMT2, have been cloned and showed notable tolerance to saline-alkaline stress in transgenic plants [3,4,5]. Recently, with the rapid development of bioinformatics, transcriptome technology has become a powerful tool for elucidating the gene regulation of networks in many species [6, 7]. Using RNA-Seq technology, the effects of gene expression on saline-alkaline stress have been widely studied in alfalfa (Medicago sativa L.), flax (Linum usitatissinum L.), jujube (Ziziphus jujuba Mill.), black locust (Robinia pseudoacacia L.), and Chinese plum (Prunus salicina Lindl.), and thousands of differentially expressed genes have been detected [8,9,10,11].

Peach (Prunus persica L.) is an important deciduous fruit tree around the world. In China, the peach planting area is more than 700 thousand hectares, and the yield reaches more than one million tons per year [12]. The North China Plain and Sichuan Basin, which partly suffer from saline-alkaline stress, are the main productive areas of peach planting [13, 14]. Saline-alkaline stress inhabits peach plant growth and affects fruit development, finally causing a reduction of in yield. GF677 (P. amygdalus x P. persica), bred by the French Institut National de la Recherche Agronomique (INRA) in the 1960s, has higher salt and alkali tolerance than Maotao (P. persica) [15]. In 2017, Chen et al. identify more than one thousand genes differentially expressed between roots of GF677 and Maotao in response to saline-alkaline stresses [16]. Considering leaf chlorosis usually acts as an early symptom of this stress [17], we carried out cytological observations and comparative transcriptome analysis of GF677 and Maotao leaves under salt-alkali stress in the present study, to elucidate the molecular mechanism governing the salt-alkali tolerance of GF677. These findings will provide a theoretical basis for genetic improvement and breeding of peach under salt-alkali stress.

Materials and methods

Plant materials

The field experiment was performed in the Late White Peach Planting Base, Jianyang City, Sichuan Province (N30°30′50.97″, E104°26′35.21″). The pH of the test soil was 8.44 and the total salt content was 0.46% [15]. Two-years-old GF677 and Maotao plantlets were planted in the field with an inter-row spacing of 1 m × 4 m. As a control, two varieties were planted in a greenhouse with normal soil (pH = 7.05). Nine GF677 and nine Maotao plantlets were selected for further experiments. The leaves from each plantlet were grouped into two parts: one for RNA isolation and another for cytological and biochemical analysis.

Ultrastructural observation

Each cultivar leaf sample was sliced into several sections, and soaked in 2% glutaraldehyde overnight at 4 °C. After rinsing with ddH2O, the samples were stained with 1% OsO4 for 2 h. The stained samples were further washed with ddH2O and dehydrated with acetone. The dehydrated samples were embedded in an epoxy resin and then sliced into ultrathin sections for ultrastructural observation as described by Jiang et al. [7].

Physiological parameters measurement

Chlorophyll a and b

Leaves (0.2 g) of GF677 and Maotao were homogenized with 80% acetone. After centrifugation at 12,000 rpm for 10 min, the supernatant was transferred separately into a clean microcentrifuge tube. The chlorophyll a and chlorophyll b contents were measured with a spectrophotometer (TU-1810, Beijing, China) at 663 nm and 645 nm, respectively [7].

Iron content

Leaves (0.5 g) of GF677 and Maotao were washed with ddH2O. The total Fe concentrations were determined by inductively coupled plasma atomic emission spectroscopy (ICP-AES; Fisons ARL Accuris, Ecublens, Switzerland) [18].

Peroxidase activity

Leaves (0.5 g) of GF677 and Maotao were homogenized with liquid nitrogen, and peroxidase (POD) activity was determined according to the protocol of the Peroxidase Assay Kit (A084–3-1, Nanjing Jiancheng Bioengineering Institute, China).

RNA extraction and transcriptome sequencing

For each cultivar, nine plantlets were randomly selected for RNA extraction. Total RNA of each plantlet was extracted according to the protocol of the RNAprep Pure Plant Plus Kit (Tiangen Biotech Co., Ltd., Beijing, China). RNA from three plantlets was equally pooled for cDNA library construction and transcriptome sequencing. For library preparation, mRNA was first extracted using oligo (dT) and then broken into 200 bp segments by adding fragmentation buffer. The first-stranded cDNA was synthesized by random hexamers, and the second-stranded cDNA was obtained by adding dNTPs and DNA polymerase I. The purified second strand cDNA was finally amplified using PCR. The libraries were sequenced on an Illumina HiSeq™ 4000 platform (Majorbio, Shanghai, China), and all raw reads were submitted to NCBI with the accession number GSE100180.

Sequence data assembly

Clean reads were obtained by filtering adapters, poly-N and low-quality reads from the raw data. The reference genome of peach (version 2.1) was downloaded from the website!info?alias=Org_Ppersica. An index of the reference genome was built using Bowtie v2.0.6, and all clean reads were then mapped to the reference genome using TopHat2 (v. 2.0.9). The transcript abundance of mapped genes was normalized by the fragments per kilobase of exon per million fragments mapped using Cufflinks (v. 2.1.1).

Differentially expressed gene (DEG) identification and functional annotation

Bioconductor package edgeR (v. 3.0.8) was used to identify DEGs. FDR ≤ 0.01 and |log2 ratio| ≥ 1 were the criteria for recognizing the significance of the gene expression difference. For visual display the DEGs by MapMan (version 3.6.0RC1) [19], all gene IDs were transformed to version 1.0 format. Gene ontology (GO) analysis of the DEGs was implemented using the GOseq package.

qRT-PCR validation

Total RNA from the leaves of the two cultivars was extracted as previously described. Total RNA (1 μg) was reverse transcribed to cDNA according to the RT reagent kit (Takara Bio Inc., Japan). The reaction mixture was used as follows: 10 μl SYBR Green 127 Premix Ex Taq, 10 μl cDNA, 0.5 μM forward primer, 0.5 μM reverse primer. The amplification was performed on a Bio-Rad CFX96 real-time system. The UBQ10 gene was used for the normalization of the tested gene [20]. The primers used in the present research are shown in Table S1.

Phylogenetic analysis

The AtR2R3-MYB and AtbHLH protein sequences were downloaded from previous research [7]. The AtWRKY protein sequences were obtained from online databases ( Other species related sequences were obtained from NCBI (Table S2). The construction of the phylogenetic tree was performed by the online software Multiple Sequence Alignment and iTOL [7].

Statistical analysis

All analyses in the present research were carried out at least in triplicate. Based on the t-test, GraphPad Prism 5 was used for statistical analysis. Double stars showed a 0.01 significant difference.


Cytological and physiological features of two peach cultivars in tolerating saline-alkaline stress

There was no difference between GF677 and Maotao in normal soil (Fig. S1). However, the phenotypes of the two varieties differed upon saline alkaline field planting (Fig. 1a). In contrast to GF677, the leaves of Maotao were clear yellowish. The TEM analysis of leaves revealed a large structural change between the two cultivars. The chloroplast of GF677 had a compact granum structure, but the structure of the granum was dispersive in the Maotao chloroplast (Fig. 1b-c). Further chlorophyll determination showed that the contents of chlorophyll a or b in Maotao were only 22.3% or 15.9% in comparison with GF677. The POD enzyme activity and Fe content decreased by 50% in Maotao (Table 1).

Fig. 1
figure 1

Phenotypes (a) and chloroplast ultrastructures (b-c) of GF677 and Maotao under salt-alkali stress. b Chloroplast ultrastructures of GF677. c Chloroplast ultrastructures of Maotao. The scale bar was shown 1 μm. GR: granum; S: starch grain; O: osmiophile globule

Table 1 Physiological parameters of GF677 and Maotao under saline-alkaline stress

Illumina sequence and DEG identification

A total of 29,006,496 and 31,874,666 clean reads were obtained in Maotao and GF677, respectively. The Q20 was more than 96.00% and the Q30 was more than 90.00%, suggesting the high quality of the clean reads. Of all clean reads, 83.68% ~ 90.43% were mapped to the reference genome (Table 2). To identify the gene response to saline-alkaline stress, the criteria of FDR ≤ 0.01 and |log2 ratio| ≥ 1 were used as thresholds. Among the DEGs, 294 and 587 were significantly up- and downregulated GF677 versus Maotao (Fig. 2a). In this study, the upregulation represented genes highly expressed in GF677 and downregulated showed genes highly expressed in Maotao. To validate the expression of DEGs, 11 genes involved in photosynthesis (ppa007547, ppa012123), carbohydrate metabolism (ppa007458), antioxidation (ppa027053, ppa011202), metal transporter (ppa003097), stress (ppa008441), polyamine biosynthesis (ppa007732), and transcription factor (ppa007708, ppa010846, ppa016095) were selected for the validation of DEGs, and the results showed that their expression pattern was highly related to the DEG data (Fig. 2b). To functionally characterize the DEGs, GO analysis was performed. A total of 30 GO terms were grouped, including several stress-related terms, such as “transporter activity”, “antioxidant activity”, “response to stimulus” (Fig. 2c). In the “molecular function” part, nine groups were obtained, and the “catalytic activity” was the most abundant group. In the “cell component” part, nine groups were obtained, and the “cell parts” was the most abundant group. In the “biological process” part, 12 groups were obtained, and the “metabolic process” was the most abundant group (Fig. 2c).

Table 2 Summary of transcriptome sequencing data from leaves of GF677 (GY) and Maotao (MY) under saline-alkaline stress
Fig. 2
figure 2

Identification (a), verification (b) and classification (c) of DEGs. a DEGs identification. The red or blue column represented up-regulation or down-regulation, respectively. b DEGs validation. 11 DEGs were randomly selected for validation by qRT-PCR. c DEGs classification. All DESs were classified by Gene Ontology (GO). 1: transporter activity; 2: structural molecule activity; 3: signal transducer activity; 4: nucleic acid binding transcription factor activity; 5: molecular function regulator; 6: electron carrier activity; 7: catalytic activity; 8: binding; 9:antioxidant activity; 10: supra molecular fiber; 11: organelle part; 12: organelle; 13:membrane part; 14: membrane; 15: macromolecular complex; 16: extracellular region;17: cell part; 18: cell; 19: single-organism process; 20: signaling; 21: response to stimulus; 22: reproductive process; 23: reproduction; 24: regulation of biological process; 25: multi-organism process; 26: metabolic process; 27: localization; 28:cellular process; 29: cellular component organization or biogenesis; 30: biological regulation

Overview of metabolism processes by MapMan

To provide an overview of the DEGs at the metabolic process level, the MapMan tool was used. It was clearly shown that photosynthesis was upregulated in GF677, while cell wall degradation, secondary metabolism, and starch degradation were activated in Maotao (Fig. 3; Table S3).

Fig. 3
figure 3

MAPMAN visualization of DEGs related to metabolic pathways. a The log2 value of significant DEGs was imported and visualized in MapMan software. Red or blue color indicated up-regulation or down-regulation, respectively. The scale bar was shown from − 3 to 3

Category and verification of genes related to antioxidation, iron metabolism, transcription factor


Seven genes encoding glutathione S-transferase (GST) and six genes encoding POD were detected in DEGs (Fig. 4a). Of these genes, all seven GST genes and four POD genes were highly upregulated in GF677. ppa011202m encoding a GST and ppa027053m encoding a POD were further confirmed by qRT-PCR (Fig. 4b).

Fig. 4
figure 4

DEGs related to antioxidation, ion, transcription factor. The left showed heat map of DEGs (a) and the right showed qRT-PCR verification (b). DEG values displayed as heat map. Colors bar represented expression levels of each gene which were either up-regulated (red) or down-regulated (green). Error bars for qRT-PCR showed the standard deviation of three replicates. ** represented 0.01significant differences. Correlation analysis of Log 2 change values obtained from RNA-Seq and qRT-PCR

Iron metabolism

Eight iron-related genes were detected in the present research. Seven of them were upregulated, including two genes (ppa013552m, ppa019605m) encoding a metal ion transporter, one gene (ppa026143m) encoding calcium-regulated channel protein, one gene (ppa001921m) encoding an iron-sulfur cluster, one gene (ppa003816m) encoding an inorganic ion transport, one gene (ppa003097m) encoding an iron ion transmembrane transporter, and one gene (ppa009703m) encoding ferritin (Fig. 4a). Of these genes, ppa003097m was further confirmed by qRT-PCR (Fig. 4b).

Transcription factor

Five genes encoding bHLH, three genes encoding MYB, and three genes encoding WRKY were detected in this study. Interestingly, the expression of bHLH and MYB downregulated, in contrast to the upregulation of WRKY (Fig. 4a). Furthermore, ppa016095m encoding a bHLH, ppa010846m encoding an MYB, and ppa007708m encoding a WRKY were confirmed by qRT-PCR (Fig. 4b).

We further constructed phylogenetic trees by integrating Arabidopsis homologues and other species functional verification genes (Fig. 5). The results showed that for each phylogenetic tree, the peach genes were scattered into different clades (Fig. 5a-c).

Fig. 5
figure 5

Phylogenetic relationships of bHLH (a), R2R3-MYB (b), WRKY (c) related DEGs with Arabidopsis and other species homolog. The trees were constructed by neighboring-joining phylogeny test, and 1000 bootstrap replicates. The accession numbers for the genes are provided in Table S2. The blue circle represented ion or salt stress-related TFs in other species. The red circle represented DEGs identified in this study.


As a serious abiotic stress, saline-alkaline stress affects plant growth and reduces yield. In peach, comparative transcriptome analysis of roots of GF677 and Maotao suggests osmotic pressure increase and redox balance are crucial for GF677 response to saline-alkaline stresses [16]. Leaf chlorosis usually acts as an early symptom of this stress [10, 11]. Our previous research showed that the photosynthetic rate was significantly higher in the leaves of GF677 than in Maotao [21]. In the present research, TEM analysis revealed that the chloroplast granum structure was more intact and the content of chlorophyll a/b was higher, in GF677, suggesting that GF677 maintains higher photosynthesis efficiency under saline-alkaline stress (Fig. 1b-c, Table 1). Iron is an important ion that plays a key role in sustaining the activity of photosynthesis-related enzymes [22, 23]. In apple, Wang et al. [24] reported that Fe deficiency induces significant downregulation of genes involved in photosynthesis. In the present study, we show that iron content is higher in GF677, which is consistent with the upregulation of several iron and photosynthesis-related genes.

Saline-alkaline stress causes the accumulation of ROS, which severely damages such cell components as DNA, lipids, proteins, and sugars [25]. Detoxification of ROS is the most important way in plants to weaken oxidative harmful stress [26]. Of the ROS scavenging enzyme family, POD and GST are the key members [27, 28]. Previous studies have shown that these proteins usually play a synergetic role in scavenging ROS. In cotton, Li et al. identified that the POD and GST proteins are upregulated upon salt treatment, and both these upregulation effects are further verified by qRT-PCR [29]. In the present research, notably, we found that POD activity was higher in GF677 (Table 1). This result is in accordance with the upregulation of most POD-related genes. Additionally, we also found that all seven GST genes were upregulated in GF677 (Fig. 4). Considering that the transcript level of a gene can largely reflect its translation state generally, we speculate that GST and POD have a synergistic role in scavenging ROS in GF677, which is consistent with the previous study [16].

TFs are key regulators of gene expression and have a variety of important functions in the plant response to abiotic stress [30, 31]. The identification of TFs involved in saline-alkaline stress is crucial to reveal the innate molecular mechanisms. In the present research, 11 genes related to bHLH, MYB, and WRKY TFs were detected (Fig. 4). As the most extensive TF class in eukaryotes, bHLH is not only universally involved in plant growth and metabolism but also plays an important role in plant response to stress [32, 33]. Fan et al. [34] revealed that Ib subgroup bHLH genes (AtbHLH38, AtbHLH39, AtbHLH100, and AtbHLH101) can bind the promoter of SKB1 and negatively regulate the accumulation of iron. In the present study, two-bHLH encoding genes (ppa024966m and ppa010295m), together with AtbHLH38, AtbHLH39, AtbHLH100, and AtbHLH101, were grouped into a subgroup by phylogenetic analysis (Fig. 5a). This result suggests ppa024966m and ppa010295m may have a similar role as Ib subgroup bHLHs, which is further confirmed by the higher content of Fe in GF677. Babitha et al. [35] showed that overexpression of AtbHLH17 enhances tolerance to NaCl, mannitol and oxidative stress in transgenic lines. Notably, ppa0271826m and ppa002645m show high homology to AtbHLH17, suggesting their putative function in the regulation of saline alkaline stress. MYB is one of the largest TF families in plants [36]. The R2R3-MYB subfamily plays a key role in diverse biological processes, especially in response to various stresses [37, 38]. In soybean, overexpression of GmMYB3a negatively regulates saline-alkaline stress-related genes [39]. In Arabidopsis, Cui et al. [40] showed that overexpression of AtMYB20 can enhance salt tolerance by negatively regulating type 2C serine/threonine protein phosphatases. In the present study, the phylogenetic tree shows that ppa010846m and ppa009439m are homologous to AtMYB20 and GmMYB3a, indicating a putative role in the regulation of saline-alkaline stress (Fig. 5b). WRKY is a class of TFs unique to plants that are mainly involved in development and stress responses [41, 42]. In Arabidopsis, overexpression of AtWRKY46, GmWRKY13 or VvWRKY11 can both positively regulate salt and drought stress tolerance [43,44,45]. In Nicotiana benthamiana, overexpression of GhWRKY41 confers transgenic plant salt and drought stress tolerance [46]. In the present study, phylogenetic analysis shows that ppa008846m is highly similar to GmWRKY13 and VvWRKY11, while ppa007708m together with GhWRKY41 and AtWRKY46 are grouped into same subgroup (Fig. 5c), suggesting that it exhibits a similar function in the regulation of saline-alkaline stress.


In the present study, the biochemical, cytological and transcriptome differences between GF677 and Maotao were systematically analysed under saline-alkaline stress. The structure of chloroplast granum was intact in GF677, but dispersed in Maotao. Functional analysis of 881 DEGs showed that photosynthesis was activated, whereas cell wall degradation, secondary metabolism, and starch degradation were repressed in GF677. Based on the phylogenetic tree integration of other species homologues, several stress-related TFs were further functionally discussed. We speculate that the integrity of chloroplast structure and the activation of photosynthesis as well as stress-related genes, are crucial for saline-alkaline resistance in GF677. The results described in this report provide a theoretical basis for cloning saline-alkaline tolerance genes and molecular breeding for improving saline-alkaline tolerance in peach.

Availability of data and materials

We have provided detailed information about the materials and methods in our manuscript.



Differentially expressed gene


Food and Agriculture Organization of the United Nations


Gene ontology


Glutathione S-transferase


Institut Nationaldela Recherche Agronomique




Quantitative real-time PCR


Reactive oxygen species


Transmission electron microscopy


  1. Pankova EI, Vorobieva LA, Balyuk SA, Khasankhanova GM, Konyushkova MV, Yamnova IA. Salt-affected soils of the Eurasian region: diagnostics, criteria and distribution. In: Vargas R, Pankova EI, Balyuk SA, Krasilnikov PV, Khasankhanova GM, editors. Handbook for saline soil management. Rome: FAO Press; 2018. p. 1–15.

  2. Fu J, Liu Z, Li Z, Wang Y, Yang K. Alleviation of the effects of saline-alkaline stress on maize seedlings by regulation of active oxygen metabolism by Trichoderma asperellum. PLoS One. 2017;2:e0179617.

    Article  CAS  Google Scholar 

  3. Guan QJ, Ma HY, Wang ZJ, Wang ZY, Bu QY, Liu SK. A rice LSD1-like-type ZFP gene OsLOL5 enhances saline-alkaline tolerance in transgenic Arabidopsis thaliana, yeast and rice. BMC Genomics. 2016;17:142.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Cao L, Yu Y, DuanMu H, Chen C, Duan X, Zhu P, et al. A novel Glycine soja homeodomain-leucine zipper (HD-zip) I gene, Gshdz4, positively regulates bicarbonate tolerance and responds to osmotic stress in Arabidopsis. BMC Plant Biol. 2016;16:184.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Jin S, Xu C, Li G, Sun D, Li Y, Wang X, et al. Functional characterization of a type 2 metallothionein gene, SsMT2, from alkaline-tolerant Suaeda salsa. Sci Rep. 2017;7:17914.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Chen L, Qu H, Xia L, Liu Y, Jiang H, Sun Y, et al. Transcriptome profiling of the fertile parentand sterile hybrid in tea plant flower buds. Hereditas. 2019;156:12.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Jiang Y, Song HY, He JR, Wang Q, Liu J. Comparative transcriptome analysis provides global insight into gene expression differences between two orchid cultivars. PLoS One. 2018;13:e0200155.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. An YM, Song LL, Liu YR, Shu YJ, Guo CH. De Novo transcriptional analysis of alfalfa in response to saline-alkaline stress. Front Plant Sci. 2016;7:931.

    PubMed  PubMed Central  Google Scholar 

  9. Guo M, Li S, Tian S, Wang B, Zhao X. Transcriptome analysis of genes involved in defense against alkaline stress in roots of wild jujube (Ziziphus acidojujuba). PLoS One. 2017;12:e0185732.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Liu J, Wang Y, Li Q. Analysis of differentially expressed genes and adaptive mechanisms of Prunus triloba Lindl. Under alkaline stress. Hereditas. 2017;154:10.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Meng F, Luo Q, Wang Q, Zhang X, Qi Z, Xu F, et al. Physiological and proteomic responses to salt stress in chloroplasts of diploid and tetraploid black locust (Robinia pseudoacacia L.). Sci Rep. 2016;6:23098.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ji P. Research on international competitiveness and export influence factors of china’s peach industryMaster Thesis, North West Agriculture and Forestry University, College of Economics and Management (In English Abstract); 2011.

    Google Scholar 

  13. Yang Z, Wang BS. Present status of saline soil resources and countermeasures for improvement and utilization in China. Shandong Agric Sci. 2015;47:125–30 (In English Abstract).

    Google Scholar 

  14. Li Y, Cao K, Zhu G, Fang W, Chen C, Wang X, et al. Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history. Genome Biol. 2019;20:36.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Gullo G, Motisi A, Zappia R, Dattola A, Diamanti J, Mezzetti B. Rootstock and fruit canopy position affect peach [Prunus persica (L.) Batsch] (cv. Rich may) plant productivity and fruit sensorial and nutritional quality. Food Chem. 2014;153:234–42.

    Article  CAS  PubMed  Google Scholar 

  16. Chen D, Song H, Sun S, Li J, Tu M, Jiang G. Scanning electron microscopy observation and comparative transcriptome analysis of roots of GF677 and Maotao in response to saline-alkali stresses. Mol Plant Breed. 2017;15:1580–6 (In English Abstract).

    Google Scholar 

  17. Song HY. Effects of chlorophyll synthesis and photosynthetic characteristics between GF677 and wide peach on alkaline soilMaster Thesis, Sichuan Agricultural University, College of Horticulture (In English Abstract); 2017.

    Google Scholar 

  18. Lei GJ, Zhu XF, Wang ZW, Dong F, Dong NY, Zheng SJ. Abscisic acid alleviates iron deficiency by promoting root iron reutilization and transport from root to shoot in Arabidopsis. Plant Cell Environ. 2014;37:852–63.

    Article  CAS  PubMed  Google Scholar 

  19. Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37:914–39.

    Article  CAS  PubMed  Google Scholar 

  20. Tong Z, Gao Z, Wang F, Zhou J, Zhang Z. Selection of reliable reference genes for gene expression studies in peach using real-time PCR. BMC Mol Biol. 2009;10:71.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Tu M, Song H, Chen D, Li J, Sun S, Jiang G. Comparison of photosynthetic characteristics and chlorophyll fluorescence parameters between GF677 and wild peach on alkaline soil in Sichuan Hilly Basin. Mountain Res. 2018;1:153–62 (In English Abstract).

    Google Scholar 

  22. Varotto C, Maiwald D, Pesaresi P, Jahns P, Salamini F, Leister D. The metal ion transporter IRT1 is necessary for iron homeostasis and efficient photosynthesis in Arabidopsis thaliana. Plant J. 2002;31:589–99.

    Article  CAS  PubMed  Google Scholar 

  23. Pereira EG, Oliva MA, Rosado-Souza L, Mendes GC, Colares DS, Stopato CH, et al. Iron excess affects rice photosynthesis through stomatal and non-stomatal limitations. Plant Sci. 2013;201-202:81–92.

    Article  CAS  PubMed  Google Scholar 

  24. Wang YX, Hu Y, Zhu YF, Baloch AW, Jia XM, Guo AX. Transcriptional and physiological analyses of short-term iron deficiency response in apple seedlings provide insight into the regulation involved in photosynthesis. BMC Genomics. 2018;19:461.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Pandhair V, Sekhon BS. Reactive oxygen species and antioxidants in plants: an overview. J Plant Biochem Biotechnol. 2006;15:71–8.

    Article  CAS  Google Scholar 

  26. Choudhury FK, Rivero RM, Blumwald E, Mittler R. Reactive oxygen species, abiotic stress and stress combination. Plant J. 2017;90:856–67.

    Article  CAS  PubMed  Google Scholar 

  27. Guan T, Shen J, Fa Y, Su Y, Wang X, Li H. Resistance-breaking population of Meloidogyne incognita utilizes plant peroxidase to scavenge reactive oxygen species, thereby promoting parasitism on tomato carrying Mi-1 gene. Biochem Biophys Res Commun. 2017;482:1–7.

    Article  CAS  PubMed  Google Scholar 

  28. Kumar S, Trivedi PK. Glutathione S-transferases: role in combating abiotic stresses including arsenic detoxification in plants. Front Plant Sci. 2018;9:751.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Li W, Zhao F, Fang W, Xie D, Hou J, Yang X, et al. Identification of early salt stress responsive proteins in seedling roots of upland cotton (Gossypium hirsutum L.) employing iTRAQ-based proteomic technique. Front Plant Sci. 2015;6:732.

    PubMed  PubMed Central  Google Scholar 

  30. Ohama N, Sato H, Shinozaki K, Yamaguchi-Shinozaki K. Transcriptional regulatory network of plant heat stress response. Trends Plant Sci. 2017;22:53–65.

    Article  CAS  PubMed  Google Scholar 

  31. Dong X, Yang Y, Zhang Z, Xiao Z, Bai X, Gao J, et al. Genome-wide identification of WRKY genes and their response to cold stress in Coffea canephora. Forests. 2019;10:335.

    Article  Google Scholar 

  32. Ohashi-Ito K, Fukuda H. Functional mechanism of bHLH complexes during early vascular development. Curr Opin Plant Biol. 2016;33:42–7.

    Article  CAS  PubMed  Google Scholar 

  33. Sun X, Wang Y, Sui N. Transcriptional regulation of bHLH during plant response to stress. Biochem Biophys Res Commun. 2018;503:397–401.

    Article  CAS  PubMed  Google Scholar 

  34. Fan H, Zhang Z, Wang N, Cui Y, Sun H, Liu Y, et al. SKB1/PRMT5-mediated histone H4R3 dimethylation of Ib subgroup bHLH genes negatively regulates iron homeostasis in Arabidopsis thaliana. Plant J. 2014;77:209–21.

    Article  CAS  PubMed  Google Scholar 

  35. Babitha KC, Ramu SV, Pruthvi V, Mahesh P, Nataraja KN, Udayakumar M. Co-expression of AtbHLH17 and AtWRKY28 confers resistance to abiotic stress in Arabidopsis. Transgenic Res. 2013;22:327–41.

    Article  CAS  PubMed  Google Scholar 

  36. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15:573–81.

    Article  CAS  PubMed  Google Scholar 

  37. Roy S. Function of MYB domain transcription factors in abiotic stress and epigenetic control of stress response in plant genome. Plant Signal Behav. 2016;11:e1117723.

    Article  PubMed  CAS  Google Scholar 

  38. Hajiebrahimi A, Owji H, Hemmati S. Genome-wide identification, functional prediction, and evolutionary analysis of the R2R3-MYB superfamily in Brassica napus. Genome. 2017;60:797–814.

    Article  CAS  PubMed  Google Scholar 

  39. He Y, Yang X, Xu C, Guo D, Niu L, Wang Y, et al. Overexpression of a novel transcriptional repressor GmMYB3a negatively regulates salt-alkali tolerance and stress-related genes in soybean. Biochem Biophys Res Commun. 2018;498:586–91.

    Article  CAS  PubMed  Google Scholar 

  40. Cui MH, Yoo KS, Hyoung S, Nguyen HT, Kim YY, Kim HJ, et al. An Arabidopsis R2R3-MYB transcription factor, AtMYB20, negatively regulates type 2C serine/threonine protein phosphatases to enhance salt tolerance. FEBS Lett. 2013;587:1773–8.

    Article  CAS  PubMed  Google Scholar 

  41. 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 

  42. Jing Z, Liu Z. Genome-wide identification of WRKY transcription factors in kiwifruit (Actinidia spp.) and analysis of WRKY expression in responses to biotic and abiotic stresses. Genes Genomics. 2018;40:429–46.

    Article  CAS  PubMed  Google Scholar 

  43. Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, et al. Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008;6:486–503.

    Article  CAS  PubMed  Google Scholar 

  44. Liu H, Yang W, Liu D, Han Y, Zhang A, Li S. Ectopic expression of a grapevine transcription factor VvWRKY11 contributes to osmotic stress tolerance in Arabidopsis. Mol Biol Rep. 2011;38:417–27.

    Article  CAS  PubMed  Google Scholar 

  45. Ding ZJ, Yan JY, Li CX, Li GX, Wu YR, Zheng SJ. Transcription factor WRKY46 modulates the development of Arabidopsis lateral roots in osmotic/salt stress conditions via regulation of ABA signaling and auxin homeostasis. Plant J. 2015;84:56–69.

    Article  CAS  PubMed  Google Scholar 

  46. Chu X, Wang C, Chen X, Lu W, Li H, Wang X, et al. The cotton WRKY gene GhWRKY41 positively regulates salt and drought stress tolerance in transgenic Nicotiana benthamiana. PLoS One. 2015;10:e0143022.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references


The authors acknowledge the contributions of field staff. The authors sincerely appreciate the helpful and constructive comments made by reviewers of the draft manuscript.


This work was funded by National Peach Industry Technical System (CARS-31-Z-12), Sichuan Sci-Tech Support Plan (16ZC2865), Innovation Capability Upgrading Project of Sichuan’s Financial Department (2016ZYPZ-019), Sichuan Science and Technology Project (2019YJ0607), National Key Research and Development Program (2018YFD020141).

Author information

Authors and Affiliations



SS carried out the experiments and wrote the manuscript. HS, JL and CD participated and analyzed data. CD and MT participated in the collection of plant materials. GJ, GY, and ZZ designed the experiment of this study. All authors have read and approved the manuscript for publication.

Corresponding author

Correspondence to Zhiqin Zhou.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Table S1.

Primers used in this paper.

Additional file 2: Table S2.

Gene accession numbers and amino acid sequences used in this paper.

Additional file 3: Table S3.

Gene expression pattern for primary and secondary metabolisms.

Additional file 4: Figure S1.

Phenotypes (A) and chloroplast ultrastructures (B-C) of GF677 and Maotao in normal soil. (B) Chloroplast ultrastructures of GF677. (C) Chloroplast ultrastructures of Maotao. The scale bar was shown 1 μm. GR: granum; S: starch grain; O: osmiophile globule

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sun, S., Song, H., Li, J. et al. Comparative transcriptome analysis reveals gene expression differences between two peach cultivars under saline-alkaline stress. Hereditas 157, 9 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: