2.1. Salt-Responsive Genes Identified by Microarray Experiment and Semi-Quantitative RT-PCR Validation
According to three screening criteria, a total of 2562 cDNAs were isolated, whose expression changed more than two-fold under stress (
Figure 1). Their corresponding genes were defined as salt-responsive. This suggested that approximately 17% of the tested cDNAs were salt-responsive, similar to results reported in wheat [
13]. A total of 1713 cDNAs were up-regulated and 849 cDNAs were down-regulated, among which the expression of 1431 cDNAs changed two- to three-fold, 493 cDNAs changed three- to four-fold, 245 cDNAs changed four- to five-fold, and 393 cDNAs changed greater than five-fold (
Figure 1). There are accumulated reports that genes which observe a fold-change of two may show no, or very little, differential expression due to the technical problems. The genes which showed changes in expression of over four-fold might be more possible to be involved in salt stress response in
Malus zumi. As Zhao
et al [
14] describe, The 235 cDNAs which changed four-to five-fold, we functionally classified, according to the Gene Ontology (GO) prediction, into 11 catagories (
Table 1): signal transduction (6%), ROS elimination (6%), osmoprotection (5%), cell maintenance and development (6%), photosynthesis (16%), transporter (5%), metabolism (17%), stress tolerance (7%), protein related (7%), others (4%), and unknown (21%) (
Table 1). A complete list of the functional classification was shown in
Supplementary file (Table S1).
The microarray results were validated by semi-quantitative RT-PCRs performed on 18 selected genes (
Figure 2).
MzTSL,
MzILR,
MzCIPK6,
MzIAA,
MzSCL,
MzSTO,
MzGTL, and
MzDREB1 were regulatory genes that function in the upstream.
MzRCA,
MzLhcb2,
MzRPE,
MzSDH,
MzHSP70,
MzNTR,
MzRD22,
MzPIP2, and
MzTIP2 were functional genes at the downstream level. As the maintaining photosynthesis efficiency under salt stress might be the one of the major mechanisms in
Malus zumi, we chose three genes (
MzRCA,
MzLhcb2,
MzRPE) involved in photosynthesis. In addition, we also select an unknown gene
MzUK1 for expression confirmation. In the microarray,
MzTSL,
MzILR,
MzIAA,
MzSTO,
MzRCA,
MzLhcb2,
MzHSP70,
MzNTR,
MzRD22,
MzPIP2,
MzTIP2, and
MzUK1 were induced by salt stress, when
MzDREB1,
MzGTL,
MzSCL,
MzCIPK6, and
MzSDH were suppressed by salt stress (
Supplementary file (Table S1)). With the exception of
MzNTR (sodium dependent phosphate transporter; BankIt1498107),
MzDREB1 (AP2 transcription factor; BankIt1494547),
MzCIPK6 (CIPK6; BankIt1494666), and
MzGTL (GT-like trihelix DNA-binding protein; BankIt1498102), the RT-PCR results of the other 14 tested genes (78%) were consistent with the microarray data. Similar levels of validation are reported in other microarray studies [
15]. This demonstrated the satisfactory quality of our experimental procedures.
Among the eight genes encoding regulatory proteins,
MzSCL (GRAS family transcription factor; BankIt1494667) was down-regulated and the other seven genes were up-regulated by salt-stress.
MzSCL belongs to GRAS transcription factor family. Its expression is reduced steadily under salt stress (
Figure 2). The seven salt-inducible genes could be divided into two groups according to the time of peak expression. The first group included
MzSTO (SALT TOLERANCE Homolog Protein; BankIt1493633),
MzDREB1,
MzILR (IAA-LEUCINE RESISTANT3; BankIt1495536), and
MzTSL (Tousled-like Serine/threonine Kinase; BankIt1494726), which responded to salt stress immediately after treatment and reached the peak expression as early as two hours (
Figure 2).
MzSTO,
MzIAA (BankIt1493693), and
MzGTL (BankIt1498102) of the second group were up-regulated after treatment and reached the expression peak at four hours (
Figure 3). The early response of the tested transcription factors and kinases suggested their possible roles in salt-stress signaling.
We also selected ten genes not involved in transcriptional regulation for RT-PCR analysis. A gene of unknown function (BankIt1498114) and
MzSDH (sorbitol dehydrogenase; BankIt1495577) were salt suppressed (
Figure 2). The gene of unknown function went down sharply and reached the lowest expression at 1.5 h after treatment. In contrast to the previous reports, three genes (
MzRCA; Ribulose-1, 5-bisphosphate carboxylase; BankIt1495565),
MzLhcb2 (Light-harvesting complex II protein; BankIt1495570), and
MzRPE (Ribulose-phosphate 3-epimerase; BankIt1498105)) involved in photosynthesis increased their expression steadily after salt treatment. The expression of
MzPIP1 (Plasma Membrane Intrinsic Protein 1; BankIt1494672),
MzTIP2 (Tonoplast Intrinsic Protein 2; BankIt1494673),
MzHSP (Heat Shock Protein 70; BankIt1498107), and
MzRD22 (RD22; BankIt1495591) were also induced by salt stress. The expression peak of
MzTIP2 appeared at eight hours, while the peak time of the other three genes appeared later than eight hours. A slight increase of expression was found for
MzNTR under stress.