Skip to main content

Salinity stress induces the production of 2-(2-phenylethyl)chromones and regulates novel classes of responsive genes involved in signal transduction in Aquilaria sinensis calli



Agarwood, is a resinous portion derived from Aquilaria sinensis, has been widely used in traditional medicine and incense. 2-(2-phenylethyl)chromones are principal components responsible for the quality of agarwood. However, the molecular basis of 2-(2-phenylethyl)chromones biosynthesis and regulation remains almost unknown. Our research indicated that salt stress induced production of several of 2-(2-phenylethyl)chromones in A. sinensis calli. Transcriptome analysis of A. sinensis calli treated with NaCl is required to further facilitate the multiple signal pathways in response to salt stress and to understand the mechanism of 2-(2-phenylethyl)chromones biosynthesis.


Forty one 2-(2-phenylethyl)chromones were identified from NaCl-treated A. sinensis calli. 93 041 unigenes with an average length of 1562 nt were generated from the control and salt-treated calli by Illmunina sequencing after assembly, and the unigenes were annotated by comparing with the public databases including NR, Swiss-Prot, KEGG, COG, and GO database. In total, 18 069 differentially expressed transcripts were identified by the transcriptome comparisons on the control calli and calli induced by 24 h or 120 h salinity stress. Numerous genes involved in signal transduction pathways including the genes responsible for hormone signal transduction, receptor-like kinases, MAPK cascades, Ca2+ signal transduction, and transcription factors showed clear differences between the control calli and NaCl-treated calli. Furthermore, our data suggested that the genes annotated as chalcone synthases and O-methyltransferases may contribute to the biosynthesis of 2-(2-phenylethyl)chromones.


Salinity stress could induce the production of 41 2-(2-phenylethyl)chromones in A. sinensis calli. We conducted the first deep-sequencing transcriptome profiling of A. sinensis under salt stress and observed a large number of differentially expressed genes in response to salinity stress. Moreover, salt stress induced dynamic changes in transcript abundance for novel classes of responsive genes involved in signal transduction, including the genes responsible for hormone signal transduction, receptor-like kinases, MAPK cascades, Ca2+ signal transduction, and transcription factors. This study will aid in selecting the target genes to genetically regulate A. sinensis salt-stress signal transduction and elucidating the biosynthesis of 2-(2-phenylethyl)chromones under salinity stress.


Aquilaria sinensis is a tropical evergreen tree widely distributed in Fujian, Guangdong, Guangxi and Hainan provinces in China and the other countries such as Vietnam, India, Indonesia, Malaysia, and Thailand [1]. Under stress conditions such as infected by fungi or wounded by wind, lighting, and bited by insects, resin-impregnated heartwoods are slowly forming in the trunk and branches of A. sinensis [24]. Those resinous heartwoods are commercially called agarwood which has been long-term used as an anti-emetic, digestive, and sedative agent in traditional medicines, and also as incense and peculiar perfume [1]. However, the production of agarwood always takes decades in natural processes, and the natural Aquilaria forests have been seriously destroyed in all countries because of the high value and great demand of agarwood. Therefore, A. sinensis has been listed in Appendix II of the Convention on Internal Trade in Endangered Species of Wild Fauna and Flora [5]. Under these circumstances, Aquilaria trees were cultivated for production of pharmaceutically important and commercially valuble agarwood using artificial methods such as burn-chisel-drill, trunk pruning, and fungi inoculation [3]. However, production of agarwood using artificial methods still takes long time, and the products are always with low quality.

Previous investigations revealed that 2-(2-phenylethyl)chromones are the principal components of argarwood [69]. There are more than 100 congeners of 2-(2-phenylethyl)chromones have been reported [10], and many of 2-(2-phenylethyl)chromones have potentially pharmacological activities including neuroprotective activity, cytotoxic activity, antibacterial activity, AchE inhibitory, anti-inflammatory activity and antioxidatic activity [7, 1114]. However, the biosynthesis and regulation of 2-(2-phenylethyl)chromones remains completely unknown.

Agarwood-producing plants are timber species which take a considerably long time to grow and the resinous portion is formed inside of the wood. It makes studies using fresh plants difficult and inconvenient. Thus, establishing calli and cell suspension cultures of A. sinensis with high production of 2-(2-phenylethyl)chromones, the principal components of agarwood, would be undoubtedly useful for the studies on the mechanism of agarwood formation [15, 16]. It has been reported that salicylic acid and the crude extracts of fungi could elicit the production of 2-(2-phenylethyl)chromones in the calli and cell suspension cultures of A. sinensis [15, 16]. We are, recently, focusing on exploring the mechanism of agarwood formation, establishment of effective method which can be used to induce the production of 2-(2-phenylethyl)chromones in calli and cell suspension is therefore critically important. Surprisingly, we firstly found that salinity stress induced the production of structurally diverse 2-(2-phenylethyl)chromones in A. sinensis calli and suspension cells, suggesting that 2-(2phenylethyl)chromones might be responsible to salt stress responses. Identification of these 2-(2-phenylethyl)chromones produced in salt-treated calli and suspension cells would be useful for further research on the biological functions of 2-(2-phenylethyl)chromones in stress responses and the mechanism of agarwood formation.

On the other hand, plants integrate complex signal pathways that may cross-talk and diverge at various steps in response to salt stress [17]. High salinity stress induces the biosynthesis of hormones to regulate the expression of specific genes and metabolites including the most important stress-responsive hormone abscisic acid (ABA) [18]. Salinity stress causes water deficit and osmotic stress, enriching the production of ABA in shoots and roots [19, 20]. The accumulation of ABA can alleviate the inhibitory influence of salinity stress on photosynthesis and growth [21]. Some other phytohormones such as salicylic acid (SA) and brassinosteroids (BR), also participate in plant responses to abiotic stress [22, 23]. Cross-talk among Ca2+ signaling pathways and mitogen-activated protein kinase(MAPK) cascades in salt stress responses have been recently been reported [2426]. Moreover, novel classes of transcription factor family members viral for signal transduction are induced by salt tress, including bZIP, WRKY, AP2/ERF and NAC families which facilitate the expression levels of various genes that eventually influence plant tolerance of salt stresses [2731]. Previous research indicated that the transcriptional expression of bZIP genes were enriched in salt-sensitive wheat variety under salt stress, but decreased in salt-tolerant cultivar [28]. In Arabidopsis, salt stress induced the expression of At WRKY8 which directly binds with the promoter of RD29A [29]. Ap2/ERF family members of rice play a significant role in salinity stress response [30]. Over expression of a NAC transcription factor family member in rice and wheat confers salt tolerance [31]. Although traditional forward and genetic approaches can provide valuable insights to salt stress responses, technical limitations may prevent further research. Genome-wide transcriptome analyses have dramatically improved the efficiency of salt stress-related gene discovery [26, 32]. In Arabidopsis, more than 20 % of the transcriptome was observed regulating under salinity stress using transcriptome analysis [32]. However, no systematic consensus on the specific classes of genes corresponding to particular signaling events in response to salt stress has been established so far. Identification and characterization of the key factors for salt stress-response signaling pathways will be meaningful for further understanding the mechanism of stress responses and the biosynthesis of specific secondary metabolites.

Herein, NaCl was demonstrated to be an ideal elicitor to induce the production of 2-(2-phenylethyl)chromones in A. sinensis calli. Using LC-MS-IT-TOF, 41 phenylethylchromones were identified from NaCl-treated A. sinensis calli. In order to elucidate the possible mechanisms of salt stress responses of A. sinensis, transcriptome sequencing was performed using Illumina sequencing technology, and the data was analyzed to identify the differentially and specifically expressed transcripts of salt-regulated genes. Concurrently, the novel classes of NaCl-responsive genes relevant to signal transduction in response to salt stress were characterized. The results provided valuable insights for further studies on the mechanism of salt stress signaling transduction and agarwood formation.

Results and discussion

Salt stress induced the production of 2-(2-phenylethyl)chromones in A.sinensis calli

To study the effects of different NaCl concentration on the biosynthesis of 2-(2-phenylethyl)chromones, 75 mM, 150 mM and 300 mM NaCl was applied to the media and 2-(2-phenylethyl)chromones were measured by LC-MS-IT-TOF system at 10 days (Fig. 1a). The peaks of tentative 2-(2-phenylethyl)chromones in the BPC profiles of calli extracts indicated that the most species and contents of 2-(2-phenylethyl)chromones were induced by 150 mM NaCl (Fig. 1b). No 2-(2-phenylethyl)chromones were produced in the control calli (no NaCl supply) (Fig. 1a and b). Our experiments indicated that the accumulation of 2-(2-phenylethyl)chromones kept constantly increasing until four weeks in NaCl-treated calli. Therefore, extracts of A.sinensis calli treated with 150 mM NaCl for 4 weeks were analyzed by LCMS-IT-TOF. The BPC profiles of A.sinensis calli extracts and mixed standards comprising 33 known 2-(2-phenylethyl)chromones isolated from agarwood are shown in Fig. 1c and d. Forty one 2-(2-phenylethyl)chromones were putatively identified on the basis of their UV and MS data, and 13 of them were unambiguously identified by comparing their retention time (R t ) on HPLC chromatogram, UV and MS data with those of authentic compounds. The other 28 compounds were tentatively identified by their predicted molecular formulas deduced from their HRESIMS data, and further confirmed by comparison of their MS/MS data with those of in literature [10]. All the data of 2-(2-phenylethyl)chromones identified from NaCl-treated A. sinensis calli are summarized in Table 1, including R t , molecular formula, calculated and experimental molecular weight (m/z), error (a relative error between calculated value and measured value) in generated molecular formula, and MS/MS data. The structures of 13 unambiguously identified 2-(2-phenylethyl)chromones from NaCl-treated A.sinensis calli are represented in Fig. 2. Previous studies showed that crude extracts of Melanotus flavolives (B.etc.) Sing. only induced four 2-(2-phenylethyl)chromones in A.sinensis cell suspension cultures [15]. In this study, we firstly used the salt treatment which is the most important abiotic stress to produce a lot of 2-(2-phenylethyl)chromones. These results indicated that salt stress was the effective method to induce the production of 2-(2-phenylethyl)chromones in calli.

Fig 1
figure 1

Analysis of 2-(2-phenylethyl)chromones from NaCl-treated A. sinensis calli by LC-DAD-IT-TOF-MS system. a Calli treated with different consistence of NaCl at 10dpi. b Effects of treatment with different consistence of NaCl on the production of 2-(2-phenylethyl)chromones. c BPCs of salt-treated Aquilaria calli extracts. d BPCs of mixed 2-(2-phenylethyl) chromone standards isolated from agarwood. BPCs: 50–1000 m/z

Table 1 Identified and tentative 2-(2-phenylethyl)chromones compounds from the salt-treated Aquilaria calli
Fig. 2
figure 2

The structures of 2-(2-phenylethyl)chromones identified with the standards

Optimization of Illumina sequencing timing

Previous experiments indicated that the cell were almost died at 10 days, in order to determine the best time for tanscriptome analysis, cell activity and 2-(2-phenylethyl)chromone accumulation in calli treated with 150 mM NaCl was investigated at 0 h, 24 h,72 h, 120 h, 168 h and 216 h. The results revealed that the activities of the cells treated with NaCl decreased significantly at 24 h and 120 h, and the activities reached to 30 % and 5 % of the control, respectively (Fig. 3a). We also discovered that the main 2-(2-phenylethyl)chromones detected from the calli were 6,7-dimethoxy-2-[2-(4'-methoxyphenyl)ethyl] chromone (compound 35) and 6,7-dimethoxy-2-(2-phenylethyl) chromone (compound 36) during the early inducing period (Fig. 3b). These two compounds are also the main 2-(2-phenylethyl)chromones in agarwood and kept increasing in the wood tissues of A. sinensis with the time of fungal infection [27]. Therefore, the occurrence of these two 2-(2-phenylethyl)chromones could be used as an important indicator for studies on the formation of agarwood. Compound 36 was firstly detected in the NaCl-treated calli for 24 h (Fig. 3b), and then the contents of compounds 35 and 36 increased constantly. The production of 2-(2-phenylethyl)chromones 35 and 36 increased remarkably at 120 h (Fig. 3b), continuously, high-quality and sufficient RNA could be isolated from the calli which were treated with 150 mM NaCl until 120 h. Therefore, three cDNA libraries from the control and induced A.sinensis calli treated with salt at 24 h and 120 h were constructed using Illumina sequencing.

Fig. 3
figure 3

Effects of NaCl treatment on the cell activity of A. sinensis calli and production of 2-(2-phenylethyl)chromones at different time points. a Relative cell activity of calli exposed to 150 mM NaCl treatment. Values was means standard error (n = 3). Means denoted by the same letter did not significantly differ at P <0.05 according to Duncan’s multiple range tests. b BPCs of 150 mM NaCl treated A. sinensis calli extracts at 0 h, 24 h and 120 h

Transcriptome sequence assembly and annotation of the unigenes

Three cDNA libraries which were generated with mRNA from control calli, and calli induced by 150 mM NaCl at 24 h and 120 h were assembled and annotated. A total of 68 962 124, 70 631 522 and 70 951 038 clean reads for control (designated control) and induced calli which treated with 150 mM NaCl for 24 h (designated induced-24 h) and 120 h (designated induced-120 h) were generated after removal of the adaptors and unknown or low-quality reads, giving a total of 6 206 591 160 nt, 6 356 836 980 nt and 6 385 593 420 nt for the control, induced-24 h and induced-120 h library (Table 2). After complete assembly of the reads, 104 316, 99 429 and 98 697 contigs with median contig size of 476 nt, 466 nt and 474 nt, were yielded from the control, induced-24 h and induced-120 h library, respectively. Further assembly analysis showed that the control, induced-24 h and induced-120 h library consisted of 91 835, 83 674 and 83 674 unigenes (Table 2), respectively. However, there were 93 041 unigenes with a mean length of 1562 nt were generated from three libraries. The length distribution of unigenes was shown in the Additional file 1: Table S1.

Table 2 Summary statistics for sequencing and sequence assembly for three libraries prepared from the control and salt-treated calli

Functions of the unigenes were annotated by BLASTX based on sequence similarity to sequences in the public databases, including NR, Swiss-Prot, KEGG, COG and GO database, and then aligned to the nucleotide database NT (E-value ≤ 1.0e−5) by BLASTN. There were 29 387 unigenes matched to one or more database and a total of 65585 unigenes were annotated. NR classification results revealed that there were 64092 unigenes matched to this database and the unigenes of A.sinensis shared 23.2 % and 20.9 % similarity to the homologs with Vitis vinifera and Ricinus communis by BLASTX annotation, followed 15.5 % and 14.5 % for Populus balsamifera and Amygdalus persica, respectively (Fig. 4). To further classify the function of the total unigenes, the unigenes were annotated by the GO, COG and KEGG database. GO analysis indicated that the unigenes were grouped into three main categories (biological process, cellular component and molecular function), which together include 55 function classes (Additional file 2: Figure S1; See supporting information). 29395 unigenes had COG annotations and were distributed in 25 clusters including the largest group “General function prediction”, followed by the group of “replication, recombination and repair” and “transcription” (Additional file 3: Figure S2; See supporting information). Gene annotation and pathway mapping in KEGG database indicated that 40552 unigenes were distributed in 128 KEGG pathways. The top three KEGG pathways which include the largest number of unigenes were metabolic pathways, biosynthesis of secondary metabolites and plant-pathogen interaction (Additional file 4: Table S2; See supporting information).

Fig. 4
figure 4

Species distribution of unigenes by BLASTX annotation. The figure indicated species distribution of unigenes BLASTX annotation with a cut-off E-value of 1.0E5. Different color showed different species

With the advent of sequencing, genome-wide analyses become available in many plant species and have significantly improved the efficiency of gene discovery. However, no genomic data is available for A.sinensis. The Illumina technology have been the first choice in field of transcriptome sequencing studies since 2012, owing to the increasing in sequencing length of reads to 150 bp or more [33]. In this study, we carried out de novo transcriptome assembly of A.sinensis calli. To date, the transcriptome information of A.sinensis was acquired by Roche 454 GS platform [2]. However, our de novo transcriptome analysis generated greater depth of sequencing, obtaining more complete coverage of the transcriptome comparing with 454 pyrosequencing in previous study. In this study, our transcriptome assembly was compared with the previous published transcriptomes using 454 pyrosequencing [2]. As shown in Additional file 1: Table S1, more than 57 % unigenes were greater than 1 kb, and more than 73 % unigenes were greater than 500 bp. However, over 70 % of the unigenes assembled by 454 GS platform were 200 bp and 600 bp long. These results indicated that our seq-RNA assembly captured larger portion of the transcriptome of A.sinensis. Owing to lack of genomic resources for A.sinensis, the proportions of unigenes which significantly corresponding to the known proteins in GenBank were considered as another useful metric. Nearly 68.89 % of our unigenes had matched in NR protein database, and this value was higher than 42.8 % reported in 454 GS platform assembly [2]. Thus, the de novo assembly of mRNA-seq will significantly improve the genome annotation of A.sinensis and be used for further study on the functional members of gene families.

Functional analysis of differentially expressed genes

To investigate the gene expression changes in response to salinity stress, the gene expression level was measured in fragment per kilobase of exon per million fragments mapped (FRKM) and the false discovery rate [FDR] < 0.001 and the absolute value of log2 Ratio ≥ 1 was used as a threshold to estimate the statistical significance of transcript expression. A total of 18069 differentially expressed genes were identified, including 454, 940 and 220 expressed uniquely in control, induced-24 h and induced-120 h libraries, respectively; 10 881 unigenes were expressed in three libraries but at different levels (Fig. 5a). After induced by salt treatment for 24 h, 5313 genes were induced while 10266 genes were down-regulated. However, 2898 genes were up-regulated and 6268 genes were down-regulated after induced by salt stress for 120 h (Fig. 5b).

Fig. 5
figure 5

Comparison of the unigenes expressed in the salt-treated and control labraries. a Venn diagram showing the unigenes expressed in the control and salt-treated calli. b Changes in gene expression profiles among the control, induced-24 h and induced-120 h calli libraries. The number of up-regulated and down-regulated DEGs between the control and induced-24 h, control and induced-120 h libraries were shown. We judged the significance of gene expression difference with the threshold FDR(False Discovery Rate) ≤ 0.001 and the absolute value of log2Ratio ≥ 1

The differentially expressed genes from the control and induced libraries were organized in functional categories by GO enrichment analysis (Fig. 6). According to the sequence homologies, the differentially expressed genes were assigned to three principal categories: biological process, cellular component and molecular function, which contained 22, 17 and 14 functional groups. Among these groups, the induced-24 h VS control and induced-120 h VS control comparisons had similar distribution of gene functions in biological process, cellular component and molecular function. However, the GO terms “cell periphery”, “extracellular region” and “external encapsulating structure” were significantly enriched after inducing by salt stress for 24 h, whereas “cell periphery”, “anchored to membrane” and “regulation of meristem growth” were primarily enriched after inducing by salt stress for 120 h (Fig. 6). We also found that a high percentage of transcripts after inducing by salt for 24 h and 120 h fell into the functional groups: “cell”, “cell part”, “cellular process”, “metabolic process” and “organelle part”(Fig. 6). To further analysis the function of differentially expressed transcripts, the differentially expressed genes were mapped in the KEGG database. The KEGG analysis indicated that 6621 DEGs with pathway annotation were distributed in 125 KEGG pathways after salt treatment for 24 h, and among these pathways, 31 pathways with P-value ≤ 0.01 were significantly influenced (Additional file 5: Table S3; See supporting information). After inducing 120 h by 150 mM NaCl, 4168 annotated DEGs transport involved in 123 KEGG pathways, and 40 pathways with P-value ≤ 0.01 were significantly enriched (Additional file 6: Table S4; See supporting information). Notably, remarkable enrichment was observed in plant-pathogen interaction pathway, stilbenoid, diarylheptanoid and gingerol biosynthesis pathway, plant hormone signal transduction pathway and pheylpropanoid biosynthesis after inducing by salt for 24 h and 120 h (Table 3). These results suggested that the transcriptome of A.sinensis calli was remarkably affected in response to salinity stress and provided resources for screening genes required for salinity stress and agarwood formation.

Fig. 6
figure 6

Gene ontology (GO) functional enrichment analysis of differentially expressed genes in salt-treated calli relative to control calli. The functions of DEGs covered three principal categories : biological process, cellular component and molecular function

Table 3 Significantly enriched KEGG pathways induced in A.sinensis in response to salt stress

Confirmation of differentially expressed candidate transcripts by quantitative real-time PCR (qRT-PCR)

To confirm the reliability of the RNA-seq results, a total of 26 candidate genes were selected for qRT-PCR analysis with specific primers (Additional file 7: Table S5; See supporting information). The analysis results suggested that all 26 DEGs selected had the same expression trends as Illumia-Solexa sequencing (Fig. 7a). For example, both qRT-PCR and RNA-seq analysis indicated that the Mitogen-activated protein kinase kinase kinase (MAPKKKA, MAPKKK2 and MAPKKK3), calmodulin, WRKY transcription factors (WRKY39, WRKY40 and WRKY75), caffeoyl-CoA-O-methyltransferase, and chalcone synthase 1(CHS1) were significantly more highly expressed in salt-treated calli compared with the control calli. Otherwise, suppression of auxin influx carrier and auxin response factor 4 by salt treatment indicated by RNA-seq analysis was verified by qRT-PCR analysis. Furthermore, a high correlation (R2 = 0.8443) was detected between RNA-seq and qRT-PCR (Fig. 7b). These results demonstrated the changes in the gene expression analyzed by RNA-seq reflecting the practical transcriptome difference between the control and salt-treatment calli.

Fig. 7
figure 7

Validation of the relative expression levels of differential expression transcripts by quantitative RT-PCR (qRT-PCR) analysis. a Expression profiles of the selected DEGs in the salt-treated calli relative to the control calli, as determined by qRT-PCR(24 h: blue; 120 h:green) and RNA-seq(24 h:red; 120 h: purple). The x-axis indicated the annotation of the selected DEGs. The y-axis indicated the normalized expression level of the genes. A: Calmodulin 1; B:Calcium-binding protein CML37; C: Calcium-dependent protein kinase 10; D: Calcium-dependent protein kinase 13-like; E: Mitogen-activated protein kinase kinase kinase A; F: Mitogen-activated protein kinase kinase kinase 2; G: Mitogen-activated protein kinase kinase kinase 2; H: Auxin influx carrier; I: Auxin response factor 4; J: G-type lectin S-receptor-like serine/threonine-protein kinase; K: Cysteine-rich receptor-like protein kinase 25; L: LRR receptor-like serine/threonine-protein kinase FLS2; M: WRKY transcription factor 75; N: WRKY transcription factor 40; O:WRKY transcription factor 29; P: Ethylene-responsive transcription factor ERF(AP2/ERF); Q: MYB-related protein MYB4; R: MYB superfamily protein 1; S: Methyltransferase PMT15; T: Caffoyl-CoA-O-methyltransferase; U: Caffeic acid 3-O-methyltransferase; V: Chalcone synthase; W: Respiratory burst oxidase homolog protein A; X:Respiratory burst oxidase homolog protein B; Y: Respiratory burst oxidase homolog protein D; Z: Pathogenesis-related protein STH-2. The transcriptional level of the selected genes was performed by qRT-PCR with three biological replications and action was used as an internal reference. Error barsrepresent the standard deviations of qRT-PCR signals (n ≥ 3). b Correlation of the expression ratio of selected DEGs analyzed by qRT-PCR and RNA-seq

Salinity stress induced a complex hormone signal pathway

Salinity stress in plant induced hormone-independent signaling pathway and hormone biosynthesis [34]. KEGG enrichment analysis showed that 2 042 out of 40 402 A.sinensis calli genes annotated as being related to hormone signal transduction pathways were detected in three libraries and 396 of these genes were differentially expressed in calli under salinity stress. Among the hormone signal transduction related DEGs, 130 were up-regulated and 202 down-regulated at 24 h induced by salt treatment, while 74 of genes were up-regulated and 163 were down-regulated at 120 h induced by salt treatment, a total of 175 were co-regulated in 24 h and 120 h for salt stress. Hierarchical clustering of the differentially expressed hormone-related genes indicated overall difference at control, 24 h and 120 h in response to salt stress (Fig. 8a). Pathway enrichment analysis indicated that 8 of hormone-independent signaling pathways including ABA, cytokinine, auxin, brassinosteroid, jasmonic acid, salicylic acid, ethylene and gibberellin pathways were induced in calli after salt treatment for 24 h and 120 h (Table 4). However, more than 70 % of DEGs related to auxin, cytokinine and ethylene pathways were down-regulated, whereas more than 50 % of DEGs required for gibberellin, salicylic acid were up-regulated. These results indicated that salinity stress induced a complex hormone signal transduction pathway.

Fig. 8
figure 8

Hierarchical clustering of differentially expressed genes involved in plant hormone signal transduction pathways, receptor-like kinases and transcription factors. Hierarchical clustering of differential expression transcripts among control, induced-24 h and induced-120 h libraries was conducted by the software PermutMatrix v1.93 and based on log2(FPKM) data for the intensity of expression of each differentially expressed genes. a Herarchial clustering of differential expression profiles related to plant hormone signal transduction pathways. b Herarchial clustering of differentially expressed genes annotated as receptor-like kinases. c Herarchial clustering of differentially expressed genes annotated as transcription factors. FPKM: numbers of fragments per kilobase of exon per million fragments mapped

Table 4 The number of differentially expressed genes involved in in the hormone signal transduction pathways

High salinity elicits rapid and constantly changes in gene expression overlapping with responses to the hormone [34]. Abscisic acid is an essential hormone that mediates plant growth in response to salinity stress through regulating the expression of many genes which encode various proteins vital for biochemical and physiological processes. ABA mediates the salt-stress signaling transduction pathway through the important signaling components, including the ABA receptor PYR/PYL family proteins, the group of protein phosphate 2C, the members of serine/threonine –protein kinase SnRK2, and ABA responsive element binding proteins [35, 36]. In A.sinensis calli during salinity stress, the proteins of ABA receptor PYR/PYL family, phosphate 2C family and ABA responsive element binding proteins family showed significantly differential expression patterns, however, the transcripts which annotated as serine/threonine-protein kinase SnRK2 showed no obvious changes(Table 4). The comparative analysis of expression profiles indicated a crosstalk between salt stress and ABA signal pathway. As well as the well-known stress-responsive ABA, other phytohormones were also responsible for salt stress. Recent studies indicated that alterations in cytokinine signaling influenced plant response to abiotic stress, seed germination, cell division, and floral initiation [37, 38]. In Arabidopsis, endogenous cytokinin overproduction decreased salt stress resistance, however, low cytokinin levels promoted salt-induced senescence in tomato [39]. Our transcriptome analysis indicated that the majority of receptor kinases and response regulators responsible for cytokinine signal transduction were significantly decreased under salt treatment (Table 4). For example, the majority of changed A.sinensis CRE1 genes which are the important cytokinin receptors were down-regulated under salt stress. Down-stream genes such as ARR-B genes indicated both up- or down-regulated under salt stress at different time point. Plants utilized auxin signaling to modify root development when responding to diverse biotic and abiotic signals including salinity stress [37, 38]. In Arabidopsis, salt stress inhibits the PIN2 (an auxin efflux carriers) expression which is known to dominate gravitropic root response by monitoring bicipital auxin transport [40]. A total of 72 genes required for the auxin signal pathway were differentially expressed under salt treatment in calli. The majority of changed Aquilaria auxin influx carrier and auxin response factor genes were down-regulated under salt stress conditions in calli, however, the genes annotated as auxin-responsive protein (IAA), auxin responsive GH3 gene family and SAUR family protein were both up- and down- regulated during salt treatment (Table 4). In addition, many Aquilaria brassinoteroid signal transduction genes, such as brassinoteroid resistant1/2 and brassinoteroid insensitive 1-associated receptor kinases (BAK1) family of genes, were also differently expressed under salt stress. In particular, the brassinoteroid resistant1/2 family genes were down-regulated during salt treatment, however, the BAK1 family members which acting on serine/threonine containing substrate and required for brassinosteroid signal transduction [41], were both up- and down-regulated during salt stress in calli at different time point (Table 4). Furthermore, families of jasmonate ZIM domain-containing protein (JAZ), transcription factor MYC2 and TGA which involved in jasmonic acid and salicylic acid signal transduction [42], respectively, were also different expressed under salt stress in calli, with majority up-regulated. Besides ABA, cytokinine, auxin, brassinoteroid, jasmonic acid and salicylic acid signal pathways in response to salinity stress, other hormones including ethylene and gibberellin also play an important role for salt stress signaling transduction. The comparative analysis of expression profiles indicated that salinity stress induced complex hormone signal pathways.

Salt stress regulated the genes associated with signal transduction: receptor-like kinases, Ca2+ -binding proteins and MAPK proteins

An appropriate defense response to salt stress initially requires recognition of the stress. One of the largest groups of receptor proteins related to the responses to various biotic and abiotic stresses is the receptor-like kinase superfamily which possesses an amino-terminal signal sequence, a transmembrane segment and intracellular serine/threonine kinase domain [4345]. We identified 688 DEGs annotated as receptor-like kinases induced in the calli under salinity stress (Additional file 8: Table S6, See supporting information). A total of 525 DEGs identified at 24 h in calli after induced by salt stress were annotated as receptor-like kinases, while a total of 409 DEGs categorized as receptor-like kinases were induced by salt treatment at 120 h. Among these DEGs, 265 were up-regulated and 263 were down-regulated induced by salt treatment after 24 h, whereas, 181 were up-regulated and 230 were down-regulated at 120 h after induced by salt, however, with 279 DEGs co-regulated in 24 h and 120 h after induced by salt. Hierarchical clustering of the differentially expressed genes annotated as receptor-like kinases suggested overall difference at control, 24 h and 120 h in the calli induced by salinity stress (Fig. 8b).

At least four receptor-like kinase families have been reported to be induced in stress-responsive genes: LRR receptor-like serine/threonine-protein kinase (LRR-RLK), proline-rich receptor-like kinase (PERK), cysteine-rich receptor-like protein kinase (CRR-RLK), and G-type lectin S-receptor-like serine/threonine-protein kinase(SRK) (Additional file 8: Table S6, See supporting information). LRR-RLK family constitutes the largest group among the receptor-like kinases, and regulates various plant processes of plant growth and development underlying the abiotic stress and biotic stress response [4446]. For instance, FLS2 which is the member of LRR-RLK, recognizes the bacterial flagellin fragment and mediates defense responses in plants [47]. Brassinosteroid insensitive 1-associated receptor kinase (BRI) which is a member of LRR-RLK and is responsible for BR hormones, regulate plant tolerance to abiotic and biotic stress [46]. Consistent with these, we observed that a total of 262 differential expressed transcripts were annotated as LRR-RLK family members, including 88 of FLS2 subgroup members and 16 of BRI subgroup members (Additional file 8: Table S6, See supporting information). For instance, a FLS2 subfamily member, were induced more than 4-fold by NaCl at 24 h and 120 h in calli, and we confirmed this by qRT-PCR. PERK family members act as sensors/receptors in cell wall and monitor changes to the cell wall during plant exposure to abiotic and biotic stresses [48]. The first characterized PERK member was Brassica napus PERK1, which expressed in different tissues and rapidly induced by wounding [49]. Another member of PERK family, PERK4, was a positive regulator in ABA and salt stress reponses. We observed that 92 of transcripts for PERK family members detected were regulated in response to NaCl treatment and about half of transcripts were up-regulated. However, we observed the majority of the CRR-RLK and SRK family members were up-regulated under salt stress in different time point, as expected for these two families of stress-inducible receptor-like kinases(Additional file 8: Table S6, See supporting information). Many of other NaCl-responsive receptor-like kinases shown in the table are potentially regulators of the NaCl-stress response in A.thaliana. These results suggested that the transcription of the majority of receptor-like genes is controlled and modulated by salt stress, implying their important roles in salt stresses signaling transduction. The downstream events of receptor-like kinase signaling relate to protein phosphorylation via a mitogen-activated protein kinase (MAPK, MAPKK, MAPKKK) cascade [50, 51]. In our analysis, salt stress influenced the expression pattern of several mitogen-activated protein kinase genes. A total of 26 differentially expressed transcripts involved in MAPK pathways and up-regulated transcripts were especially predominant in the mitogen-activated protein kinase kinase kinase family (Additional file 8: Table S6, See supporting information). For instance, three mitogen activated protein kinase kinase kinases that encoded MAPKKKA, MAPKKK2 and MAPKKK3 were transcriptionally enriched (Unigene19214_All, Unigene3063_All, CL5738.contig1_All) at 24 h and 120 h under salt treatment (Table 5), and we conformed this by qRT-PCR. However, most members of mitogen activated protein kinase kinase and mitogen activated protein kinase genes including the protein encode MPK3 and MEKK2 in A.thaliana that play an important role for stresses response were down-regulated at two time point[51, 52] (Table 5). Calcium ions (Ca2+) have been demonstrated as a second messenger in plant signaling pathways including response to salt stress. It has been reported that salt treatment trigger a transient increase in Ca2+ level [53]. The activation of MAPK through stresses is always dependent on binding of Ca2+-binding proteins including calmodulins (CaM), calcium-binding protein(CML), calcium-dependent protein kinases (CDPKs) and calcineurin B-like protein(CBL) (Additional file 8: Table S6, See supporting information) [5356]. In our analysis, 81 differentially expressed unigenes which annotated as Ca2+ dependent proteins were responsive to salt stress signal transduction and up-regulated transcripts were especially predominant in CML proteins at 24 h and 120 h under salt stress. Aquilaria calli genes encoding Ca2+ dependent proteins including CaM1, CML37, CML27, CML29 and CDPK1 which are responsible for plant innate immunity [5456] (Table 5), showed a high expression, suggesting these proteins might be a component of Ca2+ signaling to modulate plant defense response against salt stress. Thus, receptor-like kinases, MAPK cascades and Ca2+ dependent proteins involved in salt stress signal transduction.

Table 5 Expression patterns of unigenes annotated as members of MAPK cascades and Ca2+ signal pathways

Salt stress regulated the expression of transcription factors genes in calli

Transcription factors are essential for regulation of gene expression, biotic and abiotic stress responses, and signal transduction. The Arabidopsis genome possesses at least 1819 predicted transcription factors that have been classified into 56 families. However, there were 428 predicted transcription factors were regulated by salt treatment. In A.sinensis, a total of 598 DEGs were annotated as transcription factors. 512 DEGs were induced by 150 mM NaCl at 24 h, with 264 DEGs were up-regulated and 248 DEGs were down-regulated, whereas 330 DEGs were induced by NaCl at 120 h, 128 were up-regulated and 202 were down-regulated, however, only 220 DEGs were co-regulated at 24 h and 120 h (Additional file 9: Table S7, See supporting information). Hierarchical clustering of DEGs annotated as transcription factors suggested overall difference in control, induced-24 h and induced-120 h libraries (Fig. 8c). At least seven transcription factor families have been reported to be enriched in stress-responsive genes: AP2/ERF(Apetala-2/EREBP), MYB(Myeloblastosis), WRKY(named after the wrky amino acid motif), bHLH(basic helix-loop-helix), HOX(Homeodomain-containing transcription factor), NAC(NAC domain protein), and HSF(Heat Shock Factor) (Table 6). AP2/ERF is one of the largest transcription factor family involved in the response to salt stress, which can be further classified into four subfamilies: ERF, DREB, AP2, RAV [57]. Of these, ERFs and AP2s expression levels were found to be correlated with salt tolerance and activate the downstream salt-responsive genes. In A.sinensis, we observed 73 of AP2/ERF family transcripts were regulated by salt stress (Additional file 9: Table S7, See supporting information), and the up-regulated transcripts were detected in three AP2/ERF family’s subgroups(AP2, ERF, DREB) at two time point (Table 6). However, the ERF subfamily members were found to be significantly induced by salt stress at 24 h and 120 h in A.sinensis calli (Table 6). Recent analysis of the MYB family members expression level concluded that most of MYB-transcription factor are responsible for stresses and hormones [58]. Consistent with this analysis, in Arabidopsis, at least one third of the 84 predicted MYB transcription factors were enriched by NaCl at the transcript abundance level, and in wheat 16 of 60 detectable MYB family members were induced by high salt at transcript level [32, 59]. Similarly, we observed that a total of 71 DEGs were annotated as MYB family members and 14 of 24 co-regulated transcripts were enriched by NaCl (Additional file 9: Table S7, See supporting information). MYB1(Unigene9921_All) and MYB4(CL3573.Contig1_All), were induced more than 4-fold by NaCl at each of two time point (Table 6), and we confirmed this by qRT-PCR. In addition, the other MYB family members such as MYB75, MYB39, MYB78 and MYB21 were up-regulated by NaCl. The majority of WRKY family members which contain zinc-finger motifs are known to involve in biotic and abiotic response [60]. In Arabidopsis, more than 1/2 WRKY transcription factors were induced by NaCl, and the members of WRKY17, WRKY25 and WRKY33 were significantly enriched [32]. We observed that 58 DEGs were annotated as WRKY family transcription factors and the majority of the WRKY family members were up-regulated by salt treatment at 24 h and 120 h (Additional file 9: Table S7, See supporting information). We also used qRT-PCR to confirm these observations for WRKY29, WRKY40 and WRKY75 (Table 6), which were significantly enriched by salt stress. Three other large families of transcription factors with distinct DNA binding motifs are bHLH, HOX and NAC genes [32]. All of these families contained NaCl-responsive genes, however, the members of HOX and NAC families were remarkably enriched in NaCl-induced transcripts. Furthermore, most of the HSF family members were enriched by salt stress [61] (Table 6), as expected for this family of stress-responsive transcription factors. Many other transcription families including MYC, GRAS, TGA, TCP and GATA (Additional file 9: Table S7, See supporting information), contained NaCl-responsive genes and played important role for the physiological response to NaCl treatment. The large number of up- or down-regulated transcription factor genes we observed, is consistent with the existence of a complex signal transduction network underlying the response to salt stress.

Table 6 Transcription factors co-up-regulated by salt stress at different time in A.sinensis calli

Putative genes involved in biosynthesis of 2-(2-phenylethyl)chromones in calli under salinity stress

Our study showed that salinity stress induced the production of 2-(2-phenylethyl)chromones in A.sinensis calli. These results suggested that the activities of enzymes which are responsible for biosynthesis of 2-(2-phenylethyl)chromones were enriched by salinity stress. The RNA-seq analysis indicated that the diarylheptanoid and ginerol biosynthesis pathways were significantly enriched, however, curcumin synthase which is a plant polyketide synthase and is the most important synthase for diarylheptanoid and ginerol biosynthesis [62], could not be annotated the three libraries including control, induced-24 h and induced-120 h. According to the structure of 2-(2-phenylethyl)chromones, we proposed that the biosynthesis of 2-(2-phenylethyl)chromones might be similar to diarylheptanoids, and polyketide synthases (PKSs) except curcumin synthase may play a vital role in the biosynthesis of 2-(2-phenylethyl)chromones [63]. The transcriptome data indicated that the expression of chalcone synthase which is a member of type III polyketide synthase was affected after induced by salt for 24 h and 120 h (Additional file 10: Table S8, See supporting information). Among these genes, the chalcone synthase 1 (CHS1) was up-regulated in calli after induced for 24 h and 120 h; however, the expression of CHS2 was reduced at 120 h (Table 7). The phylogenetic tree was constructed according to the PKSs amino acid sequences by Neighbor-Joining method using MEGA 6.0 software to investigate the evolutionary relationship and distribution among PKSs from different species. The phylogenetic tree analysis allow the classification of PKSs in higher plants in two different clusters called CHS(chalcone-prodcuing)and non-CHS (nonchalcone-producing) cluster, and AsCHS1 was grouped in non-CHS cluster, including curcuminoid synthase (CUS) from monocots Oriza sativa, diketide-CoA synthase (DCS) and curcumin synthase (CURS) from monocots Curcuma longa, and AsCHS2 was grouped in CHS cluster [63] (Additional file 11: Figure S3). These results suggested that AsCHS1 might possess an unusual catalytic potential and could be the important enzyme for biosynthesis of 2-(2-phenylethyl)chromones. 2-(2-phenylethyl)chromones always display O-methyl groups on the phenylethylchromone scaffold. It suggested that O-methyltransferase plays a critical role in transferring O-methyl groups to 2-(2-phenylethyl)chromone scaffolds. Recently, transcriptome analysis result from the medical plant Glaucium flavum indicated that the O-methyltransferase (OMTs) were required for benzylisoquinoline alkaloid (which displays four O-methyl groups at C6, C7, C3’and C4’ on the benzylisoquinoline scaffold) biosynthesis [64]. Transcriptome resource analysis of Aquilaria calli revealed the biosynthesis of 2-(2-phenylethyl)chromones required the members of flavonol 3-O-methyltransferase and caffeoyl-CoA O-methyltransferase families transfer O-methyl groups to the 2-(2-phenylethyl)chromones scaffolds (Table 7). RNA-seq analysis showed that a total of 21 DEGs were annotated as flavonol 3-O-methyltransferase and a total of 8 DEGs were annotated as caffeoyl-CoA O-methyltransferase (Additional file 10: Table S8, See supporting information). More than half of transcripts annotated as O-methyltransferase were enriched by salinity stress at 24 h and 120 h. These results demonstrated that chalcone synthase and O-methyltransferase play a critical role in biosynthesis of 2-(2-phenylethyl)chromones.

Table 7 Expression patterns of genes involved in biosynthesis of 2-(2-phenylethyl)chromones


We firstly identified 41 2-(2-phenylethyl)chromones from the salt-treated A.sinensis calli by LC-MS system. To find the transcripts responsible salinity stress, a comprehensive transcriptome analysis was conducted from the control calli and induced calli with 150 mM NaCl treatment. A total of 93 041 unigenes with an average length of 1562 nt were obtained, and they were annotated by comparing them with the public database including NR, Swiss-Prot, KEGG, COG and GO database. Our transcriptional analysis revealed numerous genes that were differentially expressed at 24 h and 120 h under salt stress in calli. The differentially expressed candidate genes from RNA-seq were confirmed by qRT-PCR analysis. We found that numerous genes involved in hormone signal transduction, MAPK cascades signal transduction, Ca2+ signal transduction, or encoding receptor-like kinase and transcription factors showed different expressions between the control and salt-treated calli. These differences suggested that a complex signal pathway were induced in response to salinity stress. In addition, our RNA-seq analysis observed that chalcone synthase and O-methyltransferase may regulate the biosynthesis of 2-(2-phenylethyl)chromones. These results provide a better understanding of salt stresses signal transduction in A.sinensis calli, and clues to facilitate the mechanism of agarwood formation.


Plant material and chemical treatment

The A.sinensis leaves(Zhongshan, Guangdong province, China) were cut into 1 cm diameter, surface-sterilized by Clorox (2.5 %) for 10 min, disinfected using 70 % ethanol for 30 s, and washed with sterile distilled water for 4 times. The treated leaves pieces were plated on the MS medium with 2 μg/mL naphthalene-1-acetic acid (NAA) and 1 μg/mL 6-benzylaminopurine (6-BA) for inducing callus. After incubation at 25 °C for one month in dark, callus were sub-cultured onto the fresh MS medium containing 2 μg/mL NAA, 1 μg/mL 6-BA, 1 μg/mL dichlorophenoxyacetic acid (2,4-D) and 1 μg/mg kinetin(KT) every month. To identify and speculate 2-(2-phenylethyl)chromones, the calli were transferred to the medium containing 150 mM NaCl and harvested every 4 weeks for analysis by LCMS-IT-TOF. To optimize the sequencing timing, 150 mM NaCl treated calli were harvested at 0 h, 24 h, 72 h, 120 h, 168 h, and 216 h. For salt treatment, NaCl was added to the medium with a final concentration of 75 mM, 150 mM, and 300 mM. The calli subcultured for 4 weeks after inoculation were transferred onto these medium with or without NaCl, and were harvested at 10 days. Following salt treatment, the calli were collected and washed for 5 min by distilled water.

Analysis of 2-(2-phenylethyl)chromones of treated Aquilaria calli

The calli were dried at 65 °C and extracted with 95 % methanol (0.1 g dried calli with 1 mL 95 % methanol) at room temperature in an ultrasonic bath for half an hour. The extracts were then centrifuged at 12 000 rpm for one hour at 4 °C and the supernatant were analyzed with LC-DAD-IT-TOF-MS system (Shimadzu, Kyoto, Japan). For liquid Chromatography analysis, Chromatographic separations were carried out on an Agilent SB-C18 column (250 × 4.6 mm i.d., particle size 5 μm, Agilent Technologies. Palo Alto, CA, USA). The injection volume was set at 10 μL. The mobile phase was composed of acetonitrile (A) and 0.1 % aqueous formic acid (B), and delivered at 1.0 mL/min following the gradient program: 0 – 20 min, 10 % – 20 % A; 20 – 35 min, 20 % – 25 % A; 35 – 55 min, 25 % – 35 % A; 55 – 70 min, 35 % – 38 % A; 70 – 90 min, 38 % – 50 % A; 90 – 105 min, 50 – 70 % A; 105 – 120 min, isocratic 70 % A; 120 – 130 min, 70 % – 90 % A. UV absorption over 190–400 nm was recorded by DAD module. In mass spectrometer domain, an hybrid ion trap-time-of-flight mass spectrometer (IT-TOF-MS, Shimadzu) equipped with an electrospray ionization (ESI) interface was connected to LC system via a PEEK tube (0.13 mm i.d.) to perform high-resolution tandem mass spectrometry. The accurate ion axis was calibrated using the sodium trifluoroacetate (TFA) clusters as reference. Positive mass spectra were recorded in the full scan and automatic multiple stage fragmentation scan modes over a range of m/z 100 – 1000 for all MS1, MS2 spectra acquisition. The accuracy of the assigned chemical formula was determined using mass difference tolerance of ± 5 ppm, which was calculated by the deviation between the experimental mass and calculated mass.

RNA extraction and cDNA library preparation

Total RNA was extracted from the control and NaCl-induced calli at 24 h, 120 h, using the total RNA extraction kit (Norgen, Cat 72200) according to the manufacturer’s recommendations. The RNA concentration and integrity were analyzed by Nanodrop2000 (ND-100 Spectrophotometer, Peqlab). The RNA from the induced calli (three replicates) was mixed in equal proportion to generate the pool for the induced library, and the RNA from the control callus (three replicates from control calli at 24 h and three replicates from control calli at 120 h) was mixed to a pool for the control library.

Library construction and De novo transcriptome sequencing

Total RNA from control and induced calli was purified using oligo(dT) magnetic beads and fragmented into small pieces according to the manufacturer’s instruction (Illumina, San Diego, CA USA). The cleaved fragments were then used to synthesize first-strand cDNA with random hexamer primiers and reverse transcriptase (Invitrogen, Carlsbad, CA USA). Second-strand cDNA was synthesized using DNA polymerase I (Invitrogen, Carlsbad, CA USA), dNTPs and RNase H. Following, the cDNA fragment were amplified by PCR and purified after removal of the ligation adaptors and end-repair to construct the final library. Three libraries with the insert size ranging from 150 to 250 bp were sequenced on the Illumina HiSeq-2000 platform using the 100 bp paired-end approach. In total, 74,881,134, 80,382,270 and 77,355,892 raw reads were generated by Solexa/Illumina sequencing in three libraries respectively, which had been deposited in Sequence Read Archive database in NCBI, Accession No. SRA319923.

De novo Transcriptome Assembly and Annotation of Unigene Functions

After removal of the adaptor sequences, reads in low quality (more than 20%of reads having quality value ≤ 10) and reads with unknown nucleotides larger than 5 %, the clean reads assembly was carried out with short reads assembling program- Trinity program [65]. Firstly, Trinity assembled the reads with a certain length of overlap into contigs, clustered contigs into clusters and constructed de Brujin graphs for each cluster which represented the full transcriptional complexity for a given gene, and then partitioned the full read set among these graphs. Finally, Trinity processed the individual graphs, traced the paths that the reads took within the graph, ultimately programed full-length transcripts, and tested apart transcripts that corresponded to paralogous genes. The final sequences assembled by trinity were defined as ‘unigenes’. The Trinity unigenes of three libraries were further assembled with TGICL software to get non-redundant unigenes for further analysis[66]. The functions of the unigenes were aligned using BLASTX searches (with an E-value threshold of 10−5) to the public protein database including the NCBI non-redundant database (NR) (, the Swiss-Pot database (, Cluster of Orthologous Groups(COG) database (, and Kyoto Encyclopedia of Genes and Genomes (KEGG) database ( The proteins with highest sequence similarity were retrieved to analysis. COG matched each annotated sequence to an ancient conserved domain to predict and classify possible functions, while KEGG produced the inner-cell metabolic pathways. Based on NR annotation, GO functional annotation and further functional classification were carried out with Blast2GO program ( [67].

Functional Analysis of differentially expressed genes

The unigenes expressed value and transcript levels were calculated by FPKM method [68]. The rigorous algorithm which described by Stephane Audic was used to identify differentially expressed genes between different samples [69], and False discovery rare(FDR) method was applied to analyze the threshold of the P-value in multiple test. We choose “FDR ≤ 0.001 and the absolute value of log2 Ratio ≥ 1” as the threshold to identify the significance of gene expression difference. All differentially expressed genes were performed to GO enrichment analysis, and were mapped to KEGG database to identify significantly enrich metabolic pathways or signal transduction pathways.

Quantitative real-time PCR (qRT-PCR) analysis

To investigate the mRNA expression of putative genes involved in phenylethyl chromones biosynthesis and salt stress, qRT-PCR was performed using the Bio-Rad Real-time System and SYBR Green PCR Master kit according to the manufacturer’s instruction. Total RNA from different samples was reverse transcribed using the MLV kit (Sigma) according the manufacturer’s recommendations. 10 ng template cDNA was used in 20 μL PCR reaction mixture containing 10 μL fast SYBR Green master mix (Invitrogen) and 0.2 μM gene-specific primer (Additional file 7: Table S5; See supporting information). A DNA Engine Option 2 thermal cycler was carried out as the following program: one cycle of 95 °C for 3 min and 40 cycles of 95 °C for 30 s and 60 °C for 30 s. The expression of candidate genes was normalized to the Actin gene which was normally used as an internal control [4]. The relative gene expression level was normalized by rasing 2 to the power of the negative value of ∆∆Ct for each sample [70]. For all the qRT-PCR analysis, three biological experiments and three experimental replicates were calculated to analyze the candidate gene expression.



Basic Local Aligment Search Tool


CRT binding transcription factor


chalcone synthase


clusters of orthologous groups of proteins


differentially expression genes


false discovery rate


gene ontology


The KEGG resource for deciphering the genome


ncbi nr database


quantitative real-time PCR


RNA sequencing


transcription factor


  1. Naef R. The volatile and semi-volatile constituents of agarwood, the infected heartwood of Aquilaria species: A review. Flavour Frag J. 2011;26(2):73–89.

    Article  CAS  Google Scholar 

  2. Xu Y, Zhang Z, Wang M, Wei J, Chen H, Gao Z, Sui C, Luo H, Zhang X, Yang Y, et al. Identification of genes related to agarwood formation: transcriptome analysis of healthy and wounded tissues of Aquilaria sinensis. BMC Genomics. 2013;14:227–43.

  3. Zhang XL, Liu YY, Wei JH, Yang Y, Zhang Z, Huang JQ, Chen HQ, Liu YJ. Production of high-quality agarwood in Aquilaria sinensis trees via whole-tree agarwood-induction technology. Chinese Chem Lett. 2012;23(6):727–30.

  4. Gao ZH, Wei JH, Yang Y, Zhang Z, Zhao WT. Selection and validation of reference genes for studying stress-related agarwood formation of Aquilaria sinensis. Plant Cell Rep. 2012;31(9):1759–68.

    Article  CAS  PubMed  Google Scholar 

  5. China pharmacopoeia Committee. The Pharmacopoeia of People’s Republic of China(I). Beijing: Chemical Industry Press; 2010.

    Google Scholar 

  6. Chen HQ, Wei JH, Yang JS, Zhang Z, Yang Y, Gao ZH, Sui C, Gong B. Chemical constituents of agarwood originating from the endemic genus Aquilaria plants. Chem Biodivers. 2012;9(2):236–50.

  7. Chen D, Xu ZR, Chai XY, Zeng KW, Jia YX, Bi D, Ma ZZ, Tu PF. Nine 2-(2-phenylethyl)chromone derivatives from the resinous wood of Aquilaria sinensis and their inhibition of LPS-induced NO production in RAW 264.7 cells. Eur J Org Chem. 2012;27:5389–97.

  8. Ishihara M, Tsuneya T, Uneyama K. Fragrant sesquiterpenes from agarwood. Phytochemistry. 1993;33(5):1147–55.

    Article  CAS  Google Scholar 

  9. Kumeta Y, Ito M. Characterization of δ-guaiene synthases from cultured cells of Aquilaria, responsible for the formation of the sesquiterpenes in agarwood. Plant Physiol. 2010;154(4):1998–2007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Ibrahim SRM, GAM. Natural occurring 2-(2-phenylethyl) chromones, structure elucidation and biological activities. Nat Prod Res. 2015;29(16):1–32.

    Article  Google Scholar 

  11. Yoon JS, Lee MK, Sung SH, Kim YC. Neuroprotective 2-(2-phenylethyl)chromones of Imperata cylindrica. J Nat Prod. 2006;69(2):290–1.

    Article  CAS  PubMed  Google Scholar 

  12. Liu J, Wu J, Zhao YX, Deng YY, Mei WL, Dai HF. A new cytotoxic 2-(2-phenylethyl)chromone from Chinese eaglewood. Chinese Chem Lett. 2008;19(8):934–6.

    Article  CAS  Google Scholar 

  13. Yang DL, Wang H, Guo ZK, Dong WH, Mei WL, Dai HF. A new 2-(2-phenylethyl)chromone derivative in Chinese agarwood 'Qi-Nan' from Aquilaria sinensis. J Asian Nat Prod Res. 2014;16(7):770–6.

    Article  CAS  PubMed  Google Scholar 

  14. Li W, Cai CH, Dong WH, Guo ZK, Wang H, Mei WL, Dai HF. 2-(2-Phenylethyl)chromone derivatives from Chinese agarwood induced by artificial holing. Fitoterapia. 2014;98:117–23.

  15. Qi SY, He ML, Lin LD, Zhang CH, Hu LJ, Zhang HZ. Production of 2-(2-phenylethyl)chromones in cell suspension cultures of Aquilaria sinensis. Plant Cell Tiss Org. 2005;83(2):217–21.

    Article  CAS  Google Scholar 

  16. Okudera Y, Ito M. Production of agarwood fragrant constituents in Aquilaria calli and cell suspension cultures. Plant Biotechnol. 2009;26(3):307–15.

    Article  CAS  Google Scholar 

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

  18. Kempa S, Krasensky J, Dal Santo S, Kopka J, Jonak C. A central role of abscisic acid in stress-regulated carbohydrate metabolism. PLoS One. 2008;3(12):e3935.

    Article  PubMed  PubMed Central  Google Scholar 

  19. He T, Cramer GR. Abscisic acid concentrations are correlated with leaf area reductions in two salt-stressed rapid-cycling Brassica species. Plant Soil. 1996;179(1):25–33.

    Article  CAS  Google Scholar 

  20. Cramer GR, Quarrie SA. Abscisic acid is correlated with the leaf growth inhibition of four genotypes of maize differing in their response to salinity. Func Plant Biol. 2002;29(1):111–5.

    Article  CAS  Google Scholar 

  21. Popova LP, Stoinova ZG, Maslenkova LT. Involvement of abscisic acid in photosynthetic process in Hordeum vulgare L. during salinity stress. J Plant Growth Regul. 1995;14(4):211–8.

    Article  CAS  Google Scholar 

  22. Fragnire C, Serrano M, Abou-Mansour E, Métraux JP, Haridon FL. Salicylic acid and its location in response to biotic and abiotic stress. FEBS Lett. 2011;585(12):1847–52.

    Article  Google Scholar 

  23. Clause SD, Sasse JM. Brassinosteroids: essential regulators of plant growth and development. Annu Rev Plant Biol. 1998;49:427–51.

    Article  Google Scholar 

  24. Chinnusamy V, Schumaker K, Zhu JK. Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. J Exp Bot. 2004;55(395):225–36.

    Article  CAS  PubMed  Google Scholar 

  25. Huang GT, Ma SL, Bai LP, Zhang L, Ma H, Jia P, et al. Signal transduction during cold, salt, and drought stresses in plants. Mol Biol Rep. 2012;39(2):969–87.

    Article  PubMed  Google Scholar 

  26. Luo JL, Tang SH, Peng XJ, Yan XH, Zeng XH, Li J, Li XF, Wu G. Elucidation of cross-talk and specificity of early response mechanisms to salt and PEG-simulated drought stresses in Brassica napus using comparative proteomic analysis. PLoS One. 2015;10(10):e0138974.

  27. Gupta B, Huang B. Mechanism of salinity tolerance in plants: physiological, biochemical, and molecular characterization. Int J Genomics. 2014;2014:701596.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Johnson RR, Wagner RL, Verhey SD, Walker-Simmons MK. The abscisic acid-responsive kinase PKABA1 interacts with a seed-specific abscisic acid response element-binding factor, TaABF, and phosphorylates TaABF peptide sequences. Plant Physiol. 2002;130(2):837–46.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Hu Y, Chen L, Wang H, Zhang L, Wang F, Yu D. Arabidopsis transcription factor WRKY8 functions antagonistically with its interacting partner VQ9 to modulate salinity stress tolerance. Plant J. 2013;74:730–45.

    Article  CAS  PubMed  Google Scholar 

  30. Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012;1819(2):86–96.

    Article  CAS  PubMed  Google Scholar 

  31. 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(2):331–45.

    Article  CAS  PubMed  Google Scholar 

  32. Jiang Y, Deyholos MK. Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006;6:25.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Wang Y, Tao X, Tang XM, Xiao L, Sun JL, Yan XF, Li D, Deng HY, Ma XR. Comparative transcriptome analysis of tomato (Solanum lycopersicum) in response to exogenous abscisic acid. BMC Genomics. 2013;14:841.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Peleg Z, Blumwald E. Hormone balance and abiotic stress tolerance in crop plants. Curr Opin Plant Biol. 2011;14(3):1290–5.

    Article  Google Scholar 

  35. Ma Y, Szostkiewicz I, Korte A, Moes D, Yang Y, Christmann A, Grill E. Regulators of PP2C phosphatase activity function as abscisic acid sensors. Science. 2009;324(5930):1064–8.

    CAS  PubMed  Google Scholar 

  36. Sang-Youl Park PF, Nishimura N, Jensen DR, Fujii H, Zhao Y, Lumba S, Rodrigues A, Tsz-fung F. Chow, Abscisic acid inhibits PP2Cs via the PYR/PYL family of ABA-binding START proteins. Science. 2009;324(5930):1068–71.

  37. Ghanem ME, Albacete A, Martinez-Andujar C, Acosta M, Romero-Aranda R, Dodd IC, Lutts S, Perez-Alfocea F. Hormonal changes during salinity-induced leaf senescence in tomato (Solanum lycopersicum L.). J Exp Bot. 2008;59(11):3039–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Galvan-Ampudia CS, Testerink C. Salt stress signals shape the plant root. Curr Opin Plant Biol. 2011;14(3):296–302.

    Article  CAS  PubMed  Google Scholar 

  39. Rowe JH, Topping JF, Liu J, Lindsey K. Abscisic acid regulates root growth under osmotic stress conditions via an interacting hormonal network with cytokinin, ethylene and auxin. New Phytol. 2016. doi:10.1111.

    PubMed  Google Scholar 

  40. Xu WF, Jia LG, Baluska F, Ding GC, Shi WM, Ye NH, Zhang JH. PIN2 is required for the adaptation of Arabidopsis roots to alkaline stress by modulating proton secretion. J Exp Bot. 2012;63(17):6105–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Postma J, Liebrand TW, Bi G, Evrard A, Bye RR, Mbengue M, Kuhn H, Joosten MH2, Robatzek S. Avr4 promotes Cf-4 receptor-like protein association with the BAK1/SERK3 receptor-like kinase to initiate receptor endocytosis and plant immunity. New Phytol. 2016;210(2):627–42.

  42. Hossain MA, Munemasa S, Uraji M, Nakamura Y, Mori IC, Murata Y. Involvement of endogenous abscisic acid in methyl jasmonate-induced stomatal closure in Arabidopsis. Plant Physiol. 2011;156(1):430–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Shiu SH, Bleecker AB. Plant receptor-like kinase gene family: diversity, function, and signaling. Sci STKE. 2001;2001:re22.

    CAS  PubMed  Google Scholar 

  44. Shiu SH, Bleecker AB. Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis. Plant Physiol. 2003;132:530–43.

    Article  CAS  PubMed  Google Scholar 

  45. Marshall A, Aalen RB, Audenaert D, Beeckman T, Broadley MR, Butenko MA, Cano-Delgado AI, de Vries S, Dresselhaus T, Felix G, et al. Tackling drought stress: receptor-like kinases present new approaches. Plant Cell. 2012;24(6):2262–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Afzal AJ, Wood AJ, Lightfoot DA. Plant receptor-like serine threonine kinases: roles in signaling and plant defense. Mol Plant Microbe In. 2008;21(5):507–17.

    Article  CAS  Google Scholar 

  47. Giranton JL, Dumas C, Cock JM, Gaude T. The integral membrane S-locus receptor kinase of Brassica has serine/threonine kinase activity in a membranous environment and spontaneously forms oligomers in planta. Proc Natl Acad Sci U S A. 2000;97:3759–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Humphrey TV, Haasen KE, Aldea-Brydges MG, Sun H, Zayed Y, Indriolo E, et al. PERK-KIPK-KCBP signalling negatively regulates root growth in Arabidopsis thalian a. J Exp Bot. 2015;66(1):71–83.

  49. Silva NF, Goring DR. The proline-rich, extensin-like receptor kinase-1 (PERK1) gene is rapidly induced by wounding. Plant Mol Biol. 2002;50(4–5):667–85.

    Article  CAS  PubMed  Google Scholar 

  50. Li CH, Wang G, Zhao JL, Zhang LQ, Ai LF, Han YF, Sun Y. The receptor-like kinase SIT1 mediates salt sensitivity by activating MAPK3/6 and regulating ethylene homeostasis in rice. Plant Cell. 2014;26(6):2538–53.

  51. Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, Boller T, Ausubel FM, Sheen J. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002;415(6875):977–83.

    Article  CAS  PubMed  Google Scholar 

  52. Opdenakker K, Remans T, Vangronsveld J, Cuypers A. Mitogen-Activated Protein (MAP) kinases in plant metal stress: regulation and responses in comparison to other biotic and abiotic stresses. Intl J Mol Sci. 2012;13(6):7828–53.

    Article  CAS  Google Scholar 

  53. Wurzinger B, Mair A, Pfister B, Teige M. Cross-talk of calcium-dependent protein kinase and MAP kinase signaling. Plant Signal Behav. 2011;6(1):8–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Luan S, Kudla J, Rodriguez-Concepcion M, Yalovsky S, Gruissem W. Calmodulins and calcineurin B-like proteins: calcium sensors for specific signal response coupling in plants. Plant Cell. 2002;14(Suppl):S389–400.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Isope P, Wilms CD, Schmidt H. Editorial: determinants of synaptic information transfer: from Ca2+ binding proteins to Ca2+ signaling domains. Front Cell Neurosci. 2016;10:69. doi:10.3389.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Campo S, Baldrich P, Messeguer J, Lalanne E, Coca M, San SB. Overexpression of a calcium-dependent protein kinase confers salt and drought tolerance in rice by preventing membrane lipid peroxidation. Plant Physiol. 2014;165(2):688–704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006;140:411–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Chen YH, Yang XY, He K, Liu MH, Li JG, Gao ZF, Lin ZQ, Zhang YF, Wang XX, et al. The MYB transcription factor superfamily of Arabidopsis: Expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2006;60:107–24.

    Article  CAS  Google Scholar 

  59. Zhang LC, Zhao GY, Jia JZ, Liu X, Kong XY. Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress. J Exp Bot. 2012;63(1):203–14.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Dong JX, Chen CH, Chen ZX. Expression profiles of the Arabidopsis WRKY gene superfamily during plant defense response. Plant Mol Biol. 2003;51(1):21–37.

    Article  CAS  PubMed  Google Scholar 

  61. Busch W, Wunderlich M, Schoffl F. Identification of novel heat shock factor-dependent genes and biochemical pathways in Arabidopsis thaliana. Plant J. 2005;41(1):1–14.

    Article  CAS  PubMed  Google Scholar 

  62. Morita H, Wanibuchi K, Nii H, Kato R, Sugio S, Abe I. Structural basis for the one-pot formation of the diarylheptanoid scaffold by curcuminoid synthase from Oryza sativa. Proc Natl Acad Sci U S A. 2010;107(46):19778–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Abe I, Morita H. Structure and function of the chalcone synthase superfamily of plant type III polyketide synthases. Nat Prod Rep. 2010;27(6):809–38.

    Article  CAS  PubMed  Google Scholar 

  64. Farrow SC, Facchini PJ. Dioxygenases catalyze O-demethylation and O, O-demethylenation with widespread roles in benzylisoquinoline alkaloid metabolism in opium poppy. J Biol Chem. 2013;288(40):28997–9012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19:651–2.

    Article  CAS  PubMed  Google Scholar 

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

  68. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  PubMed  Google Scholar 

  69. Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95.

    CAS  PubMed  Google Scholar 

  70. Zhang X, Liao M, Chang D, Zhang F. Comparative transcriptome analysis of the Asteraceae halophyte Karelinia caspica under salt stress. BMC Res Notes. 2014;7:927.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank the Beijing Genomic Institute for the Illumina sequencing of Aquilaria sinensis calli.


This work was supported by the Beijing Natural Science Foundation (No.7154216), the National Natural Science Foundation of China (No. 81503262) and the Doctoral Foundation of Baotou Medical College(No. BSJJ201607).

Availability of data and materials

Not applicable.

Authors’ contributions

XHW carried out differential gene expression analysis, qRT-PCR and drafted the manuscript. BWG identified and speculated 2-(2-phenylethyl)chromones from salt-treated A.sinensis calli, and performed statistical analysis and visualization of results. XL, ZXZ and HYF provided the suggestion for editing and revising the manuscript. XJD, ZL and WJ participated in RNA-Seq analysis. Shepo Shi and Pengfei Tu initiated the project and supervised the work throughout. All authors read and approved the final manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Pengfei Tu.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Consent to publish

Not applicable.


Not applicable.

Additional information

Xiaohui Wang and Bowen Gao are first author. Shepo Shi and Pengfei Tu are corresponding author

Additional files

Additional file 1: Table S1.

The length distribution of unigenes. (DOCX 18 kb)

Additional file 2: Figure S1.

Gene ontology (GO) functional enrichment analysis of assembled unigenes. A total of 53514 matched unigenes were classfied into three principal categories: biological process, cellular component and molecular function. (PPTX 165 kb)

Additional file 3: Figure S2.

COG functional classification of all unigenes sequence. 29 387(31.59 %) transcripts showed significant similarity to sequences in the COG databases and were clustered into 25 categories. (PPTX 213 kb)

Additional file 4: Table S2.

The KEGG analysis of all unigenes. (XLSX 17 kb)

Additional file 5: Table S3.

The KEGG analysis of differentially expressed genes at 24 h in salt-treated calli. (XLSX 20 kb)

Additional file 6: Table S4.

The KEGG analysis of differentially expressed genes at 120 h in salt-treated calli. (XLSX 19 kb)

Additional file 7: Table S5.

The specific primers of 26 candidate genes. (XLSX 15 kb)

Additional file 8: Table S6.

The number of differentially expression genes related to receptor kinase, MAPK pathway and Ca2+ signal pathway. (DOCX 22 kb)

Additional file 9: Table S7.

The number of differentially expressed genes related to transcription factors. (DOCX 20 kb)

Additional file 10: Table S8.

The number of differentially expressed ungenes involved in biosynthesis of 2-(2-phenylethyl)chromones. (DOCX 19 kb)

Additional file 11: Figure S3

Phylogenetic tree analysis of plant type III PKSs. The tree was constructed by neighbor-joining algorithm and the reliability of the tree was measured by bootstrap analysis with 1000 replicates. The indicated scale represents 0.05 amino acid substitutions per site. (PPTX 150 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( 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

Wang, X., Gao, B., Liu, X. et al. Salinity stress induces the production of 2-(2-phenylethyl)chromones and regulates novel classes of responsive genes involved in signal transduction in Aquilaria sinensis calli. BMC Plant Biol 16, 119 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: