De novo sequencing and comparative transcriptome analysis of adventitious root development induced by exogenous indole-3-butyric acid in cuttings of tetraploid black locust
- Jine Quan†1, 2,
- Seng Meng†1,
- Erhui Guo2,
- Sheng Zhang1,
- Zhong Zhao1Email author and
- Xitian Yang1, 2Email author
© The Author(s). 2017
Received: 28 July 2016
Accepted: 3 February 2017
Published: 16 February 2017
Indole-3-butyric acid (IBA) is applied to the cuttings of various plant species to induce formation of adventitious roots (ARs) in commercial settings. Tetraploid black locust is an attractive ornamental tree that is drought resistant, sand tolerant, can prevent sand erosion and has various commercial uses. To further elucidate the mechanisms of AR formation, we used Illumina sequencing to analyze transcriptome dynamics and differential gene expression at four developmental stages in control (CK) and IBA-treated groups.
The short reads were assembled into 127,038 unitranscripts and 101,209 unigenes, with average lengths of 986 and 852 bp. In total, 10,181 and 14,924 differentially expressed genes (DEGs) were detected in the CK and IBA-treated groups, respectively. Comparison of the four consecutive developmental stages showed that 282 and 260 DEGs were shared between IBA-treated and CK, suggesting that IBA treatment increased the number of DEGs. We observed 1,721 up-regulated and 849 down-regulated genes in CI vs. II, 849 up-regulated and 836 down-regulated genes in CC vs. IC, 881 up-regulated and 631 down-regulated genes in CRP vs. IRP, and 5,626 up-regulated and 4,932 down-regulated genes in CAR vs. IAR, of which 25 up-regulated DEGs were common to four pairs, and these DEGs were significantly up-regulated at AR. These results suggest that substantial changes in gene expression are associated with adventitious rooting. GO functional category analysis indicated that IBA significantly up- or down-regulated processes associated with regulation of transcription, transcription of DNA dependent, integral to membrane and ATP binding during the development process. KEGG pathway enrichment indicated that glycolysis/gluconeogenesis, cysteine and methionine metabolism, photosynthesis, nucleotide sugar metabolism, and lysosome were the pathways most highly regulated by IBA. We identified a number of differentially regulated unigenes, including 12 methionine-related genes and 12 ethylene-related genes, associated with the KEGG pathway cysteine and methionine metabolism. The GO enrichment, pathway mapping, and gene expression profile analyses revealed molecular traits for root induction and initiation.
Our study presents a global view of the transcriptomic profiles of tetraploid black locust cuttings in response to IBA treatment and provides new insights into the fundamental mechanisms associated with auxin-induced adventitious rooting.
Tetraploid black locust (Robinia pseudoacacia L.) is an attractive ornamental tree that has various commercial uses. It exhibits fast growth, drought resistance, saline-alkaline tolerance and low soil fertility requirements, and it is the primary tree species used as a wind break and for sand fixation and for soil and water conservation in the Loess Plateau Region . Tetraploid black locust does not root well and thus is difficult to plant in a variety of environments. However, the phytohormone auxin can promote the formation of adventitious roots (ARs) in cuttings [1, 2]. The formation of ARs relies on the method used. In agricultural practice, plant loss is usually caused by AR formation and slow-rooting cuttings. Therefore, it is believed that the formation of ARs is necessary for the smooth spreading of cuttings of tetraploid black locust.
AR formation is a highly complex regenerative process that is influenced by numerous internal and external factors, including environmental conditions, phytohormones and nutritional status [3–6]. Auxin is a crucial phytohormone that promotes AR formation in cuttings . However, the mechanisms underlying the role of auxin in AR formation are only superficially understood, and the lack of details at the molecular level limits improvements to cutting propagation.
The auxin IBA is widely used in woody plant propagation to induce rooting. Although indole-3-acetic acid (IAA) is the primary native auxin in plants, IBA is more effective in promoting ARs [8–11]. For example, treatment of cuttings with IBA significantly improves rooting rates in Terminalia tomentosa , Pinus contorta [12, 13], Malus pumila [14, 15], and Pinus radiate . Brinker et al.  reported that IBA induces the expression of genes involved in cell replication and cell-wall weakening but inhibits genes related to auxin transport, photosynthesis and cell-wall synthesis during P. contorta root initiation. Thus, the processes that occur in cuttings after IBA treatment, and particularly the functions of IBA-regulated genes, should be elucidated. These results also indicate that IBA may directly or indirectly induce the formation and differentiation of root primordia. We previously demonstrated that 5.4 mmol/L IBA significantly increases the rooting rate of tetraploid black locust hardwood cuttings to approximately 80% [1, 2], in contrast to approximately 2% AR formation in CK cuttings. To explore the significant impact of IBA on rooting, our recent studies have primarily focused on tetraploid black locust at the anatomical, physiological and biochemical levels [1, 2]. Developments in molecular biology and proteomics techniques have allowed us to investigate IBA-induced AR development in tetraploid black locust via homology cloning, quantitative real-time PCR (qPCR), two-dimensional electrophoresis and protein bio-mass spectrometry (MALDI-TOF, Q-TOF), resulting in the isolation and identification of hundreds of IBA-response-related genes and proteins [18–20]. The studies cited above represent the first explorations of genes involved in AR formation in tetraploid black locust, but transcriptomic information and identification of many genes related to IBA-induced AR development is scarce. The molecular mechanism of rooting in tetraploid black locust is complex, and the mechanism by which IBA promotes the formation of ARs in cuttings remains unclear due to the lack of transcriptomic and genomic information. Therefore, our study represents a necessary acceleration of the acquisition of transcriptomes related to IBA-induced AR development in tetraploid black locust cuttings. Transcriptomic studies of IBA-induced AR development have been conducted in Camellia sinensis , Petunia hybrid , Pinus contorta , and Morus alba L. . These studies primarily focused on stem cuttings and involved techniques such as Illumina sequencing , massively parallel signature sequencing , EST analysis , microarray analysis, and suppression subtraction hybridization . However, the transcriptome of tetraploid black locust has not been studied yet.
The rapid development of next-generation sequencing (NGS) technology has improved the efficiency and reduced the cost of Illumina/Solexa sequencing technology. The results of Illumina/Solexa sequencing are highly reproducible, both technically and biologically [23, 25]. Illumina/Solexa sequencing technology is the most widely used NGS technology for the de novo sequencing and analysis of the transcriptomes of non-model organisms.
In this study, we used the Illumina sequencing platform to analyze CK and IBA-treated cuttings of tetraploid black locust in order to identify new genes involved in the IBA-induced formation of ARs in cuttings and to obtain deeper insight into the mechanism of propagation in tetraploid black locust. The application of Illumina next-generation sequencing provides more transcripts to facilitate further genomic studies of tetraploid black locust. This study presents the transcriptome for IBA-treated cuttings and provides a genetic resource for improving woody plant propagation.
IBA induced adventitious root formation in tetraploid locust
Illumina sequencing, de novo assembly and sequence analysis
Summary for RNA-Seq data of tetraploid black locust
Transcriptome changes during AR development in tetraploid black locust cuttings
Moreover, the expression levels of 282 DEGs (group 23) in the IBA-treated cuttings and of 260 DEGs (group 5) in CK increased continuously. The expression of 61 DEGs increased continuously in both IBA-treated and CK cuttings (Additional file 4). These results suggest that a substantial alteration of gene expression is associated with adventitious rooting. For example, the expression of the beta HLH protein 93 (Trb7468921), early-responsive to dehydration stress protein (ERD4) (Trb4932701), ethylene responsive element binding factor 1 (Trb5562722), peroxidase 2 (Trb4932701), and zinc-finger protein 1 (Trb8690271) DEGs increased continuously in the IBA-treated cuttings. In addition, the expression of the C2H2-like zinc finger protein (Trb5430271), cold, circadian rhythm, and RNA binding 1 (Trb6574322), embryo-specific protein 3 (Trb4689243), and heat shock protein 60 (Trb7645302) DEGs increased continuously in the CK cuttings. Moreover, the expression of the ACC oxidase 1 (Trb7020711), basic helix-loop-helix (bHLH) DNA-binding protein 5 (Trb8657204), cell wall/vacuolar inhibitor of fructosidase 2 (Trb5894705), LOB domain-containing protein 41 (Trb6654701), and CAP160 protein (Trb3457711) DEGs increased in both IBA-treated and CK cuttings.
DEGs in IBA effects on adventitious root development
DEGs in response to IBA during the four developmental stages of adventitious rooting
Analysis using the Z-score method suggested that in each of the four developmental stages, the expression levels of 16,325 genes differed significantly between the transcriptomes of IBA and CK at P < 0.05 and fold-change > 2. Of these, 2570 unigenes, 1721 up-regulated and 849 down-regulated, were differentially expressed in CI vs. II; 1,685 unigenes, 849 up-regulated and 836 down-regulated, were differentially expressed in CC vs. IC; 1,512 unigenes, 881 up-regulated and 631 down-regulated, were differentially expressed in CRP vs. IRP; and 10,558 unigenes, 5,626 up-regulated and 4,932 down-regulated, were differentially expressed in CAR vs. IAR (Fig. 4). In addition, 10,588 genes (64.86% of the total) showed differential expression between CK and IBA treatment only during the last stage. The identification of more differentially expressed genes in the final stage might be related to the greater distinction of this developmental stage. Moreover, in all four developmental stages, there were more up-regulated genes than down-regulated genes. These results suggest that IBA treatment increased the number of genes that were up-regulated to promote adventitious rooting.
Only 25 up-regulated DEGs were detected in all four stages (Additional file 7). These included seven individual genes, namely, ACC oxidase 1 (ACO1, Trb70207), arabinogalactan protein 22 (AGP22, Trb44105), flavodoxin-like quinone reductase 1 (FQR1, Trb50725), multidrug resistance-associated protein 5 (AtMRP5, Trb74407), photosystem II light harvesting complex gene 2 (LHCB2, Trb78279), response regulator 9 (ARR9, Trb73071) and uclacyanin 3 (UCC3, Trb51917), and ten gene families (13 genes), namely, an ARM repeat superfamily protein (Trb61444), an auxin-responsive family protein (Trb63108), a bifunctional inhibitor/lipid-transfer protein/seed storage 2S albumin super family protein (Trb55916), a C2H2-type zinc finger family protein (Trb25048), an NmrA-like negative transcriptional regulator family protein (Trb76075), a P-loop-containing nucleoside triphosphate hydrolase super family protein (Trb68969), a rhodanese/cell cycle control phosphatase super family proteins (Tr73938), two thiamine pyrophosphate-dependent pyruvate decarboxylase family protein (Trb64980 and Trb66324), three wound-responsive family proteins (Trb51890, Trb68418 and Trb51890) and a WRKY family transcription factor (Trb62764), as well as five proteins of unknown function (Trb69581, Trb64606, Trb43851, Trb73653 and Trb55166). These genes likely play an important role in the development of adventitious rooting in response to IBA treatment .
Additionally, several DEGs encode auxin-related products, such as indole-3-acetic acid synthetase (Trb7222702) and serine/threonine-protein kinase protein (Trb7444204). Six of these DEGs encode a1-aminocyclopropane-1-carboxylate oxidase (Trb7020714), which was previously reported to be involved in ethylene biosynthesis .
GO enrichment analysis
The GO analysis of the DEGs in the four developmental stages between IBA-treated and CK revealed that most of the encoded products were associated with the following GO categories: regulation of transcription, transcription of DNA dependent, integral to membrane and ATP binding. The most common categories associated with the AR stage were integral to membrane and ATP binding activity (Additional file 8). Several genes showed highly significant differences, including cytokinin dehydrogenase (Trb5994202), 1-aminocyclopropane-1-carboxylate oxidase (Trb7341803), 4-coumarate: CoA ligase-like 9 (Trb6655606), ubiquitin protein ligase RGLG1-like (Trb5586102), S-adenosylmethionine synthase (Trb5508801 and Trb9018001), spermidine synthase, S-adenosylmethionine decarboxylase (Trb6300002 and Trb4357121) and zeaxanthin epoxidase (Trb7243506). These auxin-responsive genes have important functions in AR formation (Additional file 9). In addition, highly significant differences in gene expression were detected, most of which were distributed among the following GO terms: cytokines in metabolic process, ethylene biosynthetic process, jasmonic acid biosynthetic process, auxin metabolic process, S-adenosylmethionine biosynthetic process, spermidine biosynthetic process and abscisic acid biosynthetic process.
The top 50 most significantly up- and down-regulated GO categories during the IBA-treated development process are listed in Additional file 10. The results revealed that the GO classifications associated with organism development, such as AR development, xylem development, phloem development, post-embryonic root development, and organ development, were significantly up-regulated in CI compared with II but significantly down-regulated at stage IRP compared with IAR. Hormone-related pathways, such as ethylene-mediated signaling pathway and auxin metabolic process, were significantly up-regulated in CI compared with II, whereas genes with the GO classifications abscisic acid biosynthetic process, jasmonic acid biosynthetic process and cytokinin metabolic process were significantly up-regulated at stage IRP compared with IAR.
KEGG pathway enrichment analysis
To further determine which biological pathways were significantly (P ≤ 0.05) modulated during AR formation, KEGG pathway enrichment was performed using the KEGG Automatic Annotation Server (KAAS)  to reveal KEGG pathway enrichment in the transcriptomes of IBA-treated and CK cuttings at four developmental stages. Between stages CI and II, thirteen KOs were significantly down-regulated or up-regulated, including glycolysis/gluconeogenesis, cysteine and methionine metabolism, photosynthesis, amino sugar and nucleotide sugar metabolism, lysosome, and starch and sucrose metabolism. Between stages CC and IC, five KOs were significantly down-regulated or up-regulated, including alanine, aspartate and glutamate metabolism, bacterial secretion system, butanoate metabolism, and starch and sucrose metabolism. Further, when comparing stages CRP and IRP, seven KOs were significantly down-regulated or up-regulated, including Alanine, aspartate and glutamate metabolism, beta-Alanine metabolism, and Type I diabetes mellitus. Meanwhile, comparing stages CAR and IAR, twenty-five KOs were significantly down-regulated or up-regulated, including glycolysis/gluconeogenesis, carbon fixation in photosynthetic organisms, alzheimer’s disease, starch and sucrose metabolism, cysteine and methionine metabolism, arginine and proline metabolism, MAPK signaling pathway, ABC transporters, valine, leucine and isoleucine degradation, antigen processing and presentation, nitrogen metabolism, tryptophan metabolism and selenoamino acid metabolism (Additional file 11). These results suggest that significant metabolic changes occur during the period of AR formation.
To gain insight into the differential KOs specifically induced in cuttings by IBA at the four developmental stages, the pairs of consecutive stages in the IBA-treated group were compared: II vs. IC, IC vs. IRP, and IRP vs. IAR. Forty KOs were significantly down-regulated or up-regulated between stage II and stage IC, including glycolysis/gluconeogenesis, carbon fixation in photosynthetic organisms, pentose phosphate pathway, pyruvate metabolism, amino sugar and nucleotide sugar metabolism, cysteine and methionine metabolism, and alanine aspartate and glutamate metabolism. These results indicate that glycolysis/gluconeogenesis and carbon fixation in photosynthetic organisms were significantly up-regulated by IBA treatment in stage IC relative to stage II. In stage IRP compared with stage IC, the KO Phenylalanine metabolism was significantly down-regulated, and the group Starch and sucrose metabolism was significantly up-regulated. Phenylpropanoids contribute to plant defenses as inducible chemical barriers or as signaling molecules [35–37]. From stage IRP to stage IAR, fifteen KOs were significantly down-regulated or up-regulated, including alanine, aspartate and glutamate metabolism, carbon fixation in photosynthetic organisms, glycolysis/gluconeogenesis, arginine and proline metabolism, cysteine and methionine metabolism, pentose phosphate pathway, oxidative phosphorylation, and photosynthesis. cysteine and methionine metabolism, as well as ethylene pathway and associated polyamines, might play important roles in IBA-induced adventitious rooting . Furthermore, the results further suggest that IBA increased the number of KOs and genes exhibiting changes in expression during the early stages of rooting.
Verification of DEGs during the four developmental stages of adventitious rooting
The genes threonine kinase (Trb6720901), auxin-repressed protein (Trb2435601), ethylene responsive transcription factor (Trb6253801), auxin-induced in root cultures protein (Trb6310802), auxin-responsive protein IAA (Trb8029901), and SAMDC (Trb7651101) beta-D-xylosidase (Trb6237301) remained highly expressed during AR development for all cutting stages, and the expression of these genes was higher in IBA-treated cuttings than in CK cuttings in all four stages. We observed the highest expression of flowering promoting factor protein (Trb5566601) in the callus induction phase of IBA-treated cuttings, as well as high expression in the AR formation phase, and the expression of this gene in the AR formation phase was higher in IBA-treated cuttings than in CK cuttings.
We observed the highest expression levels of heat shock cognate protein (Trb5873209), zinc finger protein (Trb6191801), ACO (Trb7020711), heat shock cognate protein (Trb5873209), SAMS (Trb6301201) wound-induced protein (Trb7602201), and aspartic protease in guard cell (Trb7014201) in the AR formation phase of IBA-treated cuttings, and the expression of these genes in the AR formation phase was higher in IBA-treated cuttings than in CK cuttings. The most highly expressed protein detected was phloem protein (Trb6463501) in the root primordia formation phase of IBA-treated cuttings, as well as high expression in the AR formation phase, and expression of this gene in the AR formation phase was higher in IBA-treated cuttings than in CK cuttings. We observed the highest expression levels of Phloem WRKY transcription factor (Trb6276401), 2,4-D inducible glutathione S-transferase (Trb5981501), and peroxidase (Trb6791601) during the initiation formation phase of IBA-treated cuttings, as well as high expression in the AR formation phase, and expression of these genes in the AR formation phase was higher in IBA-treated cuttings than in CK cuttings. These results demonstrate that IBA might directly or indirectly regulate the expression of the above genes during AR development in tetraploid black locust. Several members of the threonine kinase, auxin-repressed protein, ethylene responsive transcription factor, auxin induced in root cultures protein, auxin-responsive protein IAA, ACO and SAMDC beta-D-xylosidase families have been identified and shown to mediate adventitious rooting [22, 41]. However, during the AR formation phase, the expression levels of the genes cytochrome p450 (Trb5844701), Delta-l-pyrroline-5-carboxylate synthetase (Trb6734101), and transcription factor bHLH (Trb6501202) were higher in CK cuttings than in IBA-treated cuttings. These results show that these genes were less affected by IBA treatment during the stages of AR development. Instead, cytochrome p450 and transcription factor bHLH genes are involved in responses to stress, such as drought and high salinity, thus leading to efficient adventitious rooting .
Unigene determination of transcriptome sequencing during AR development in tetraploid black locust cuttings
Illumina RNA-seq technology has been extensively used for model plant transcriptome sequencing  with reference genome data or for non-model plants  without reference genomic information. In this study, the Illumina HiSeq 2000 platform was used to perform a de novo transcriptome sequencing analysis of the tetraploid black locust cuttings to better understand the gene expression changes during adventitious rooting. Pooled RNA samples from IBA-treated and CK cuttings sampled at four time points after AR excision were used to construct cDNA libraries for deep sequencing. In this sequencing, approximately 33.74 million to 46.81 million paired-end clean reads were obtained from the IBA-treated and CK cuttings at the four time points. After de novo assembly, we obtained 101,209 unigenes with a mean length of 986 bp, which is longer than has been reported previously in studies using the same technology [26, 28, 43, 44].
We identified a total of 10,181 and 14,924 DEGs (fold change > 2) were detected in at least one of the four stages in the CK and IBA treated cuttings, respectively. The genes that exhibited differences in expression between two consecutive rooting stages are shown in Fig. 4. A comparison of tetraploid black locust cuttings in stage II and stage IC revealed that 8,976 genes were differentially expressed, of which 3,948 were down-regulated and 5,028 were up-regulated in stage IC. Using a DNA microarray method, Rigal et al.  studied gene expression changes during adventitious rooting in the model tree Populus trichocarpa. Their results indicated that 5,781 genes were differentially expressed in the organization of the AR primordium; 6,538 genes were differentially expressed during primordium differentiation; and 1,146 genes were differentially expressed between these two stages . In another similar study using cDNA microarrays, Brinker et al.  identified 220 genes whose expression changed significantly during root development in hypocotyl cuttings of Pinus contorta . The results obtained suggest that RNA-Seq is a sensitive, low-cost, and accurate method for deep-sequencing the transcriptomes of plant without available genomic information, and this method was able to identify more DEGs during the early stages of adventitious rooting relative to the results of DNA microarrays. This technology also enables the precise elucidation of transcripts in the samples.
The associations between the GO terms and the lists of DEGs were investigated using GO functional enrichment analysis. Significantly enriched GO terms for the DEGs in the II vs. CI, IC vs. CC, IRP vs. CRP and IAR vs. CAR comparisons included cellular component, cell part, membrane and membrane-bound organelle (Fig. 5). Among biological processes, the highest number of unigenes belonged to cellular metabolic process, and the top three classes of genes were primary metabolic process, macromolecular metabolic process and response to stimulus. Additionally, the molecular function and cellular metabolic process terms contained the highest numbers of unigenes, followed by metabolic process, macromolecular metabolic process and response to stimulus. The GO functional enrichment analysis and unigene expression abundance data will provide useful information for the identification of genes involved in AR development in tetraploid black locust.
Genes involved in cysteine and methionine metabolism were significantly regulated by IBA during AR formation
We further examined the genes encoding proteins involved in Cysteine and methionine metabolism during the process of adventitious rooting. We identified a number of unigenes (fold change > 2), including 12 methionine-related genes, and 12 ethylene-related genes, associated with the KEGG pathway Cysteine and methionine metabolism (Additional file 11). Among those 62 genes, 60 were identified as auxin-related. The genes down-regulated at stage II included a threonine synthase gene and four S-adenosylmethionine decarboxylase genes; between stages II and IC, the down-regulated genes included six S-adenosylmethionine synthetase-like genes and an aspartate kinase-like gene. Compared with stages II and IC, a total of 21 genes, including 11 up-regulated and 10 down-regulated, were identified as differentially regulated at stages IRP and IAR. The up-regulated genes were primarily members of the S-adenosylmethionine synthetase family, while the down-regulated genes were mostly members of the adenosylhomocysteinase family. S-adenosylmethionine decarboxylase has been known to function as an auxin carrier complex in cellular auxin efflux and influx . These results indicate that most of the genes related to methionine metabolism were significantly up-regulated by IBA treatment. The up-regulation of aspartate aminotransferase and the DNA methyltransferase during stages II and IC suggests that auxin transport occurs in these stages. In other studies, the expression levels of spermidine synthase and S-adenosylmethionine decarboxylase, which are essential for AR formation , were up-regulated by IBA treatment. Among the ethylene-related genes, 17 DEGs were identified at stage IRP, with 15 up-regulated and two down-regulated, and 11 DEGs were identified at stage IAR, with eight up-regulated and three down-regulated. Compared with stages II and IC, stages IRP and IAR showed a total of 20 DEGs, with 9 up-regulated and 11 down-regulated.
The DEGs in each developmental stage in the CK and IBA treatments
To further understand the roles of the DEGs in each developmental stage in the CK and IBA-treated cuttings, we further analyzed the expression levels of a total of 25 genes that showed differential expression in the samples. Among these, several genes were differentially expressed in the AR formation phase, with higher expression in IBA-treated cuttings than in CK, including ACC oxidase 1 (ACO1). Interestingly, the expression of these genes is induced by exogenous auxin in the rooting-competent cuttings of two distantly related forest species, and there are some similar reports that ACO genes lead to AR formation in mung bean and tomato [46, 47]. It appears that IBA-induced ethylene production may contribute to the stimulation of adventitious rooting . Arabinogalactan protein 22 (AGP22) and other AGPs are extracellular proteoglycans that are implicated in many plant growth and developmental processes; for example, AtAGP30 is a non-classical AGP core protein from Arabidopsis that is expressed only in roots . Laskowski et al. report that FQR1 is a novel primary auxin-response gene that encodes a flavin nucleotide-binding flavodoxin-like quinone reductase, and accumulation of FQR1 mRNA can be detected in roots and root cultures . AtMRP5, as a newly identified member of the ABC transporter superfamily, controls root development in Arabidopsis thaliana, according to Gaedeke et al. . Further, Zhang et al. and To et al. have also reported that the genes ARR8/ARR9 are expressed in the root and that they act partially redundantly to negatively regulate the response of roots to exogenous cytokinin during root meristem initiation [51, 52]. Plants possess a superfamily of arm-repeat proteins, which were reported to be involved in abscisic acid signaling by Kim et al.  and in lateral root development by Coates .
The DEGs detected encompassed ten gene families. Among the auxin-responsive factor (ARF) family proteins, ARF2 and ARF5 were up-regulated both in CI compared with II and in CAR compared with IAR. Two auxin-responsive factor genes have been found to be involved in auxin signaling and to regulate adventitious rooting . Zinc finger family proteins, specifically of the basic leucine zipper and C2H2 types, were up-regulated in CI compared with II and in CRP compared with IRP. Similar reports have indicated that basic leucine zipper genes promote the modulation of meristems and primordial development. The two wound-responsive family genes and two WRKY transcription factor family genes were significantly down-regulated in CAR compared with IAR and were significantly regulated by IBA. However, a conflicting report by Li et al. indicated that zinc finger proteins and WRKY were down-regulated by IBA treatment during the early stages of adventitious rooting . These results indicate that these differentially expressed genes were directly or indirectly regulated by IBA.
Transcriptome sequence data for CK and IBA-treated tetraploid black locust cuttings at four developmental stages were obtained using the Illumina sequencing method, with subsequent de novo assembly. Despite the economic importance of this tree, the tetraploid black locust genome is not publicly available, and sequence data for tetraploid black locust are limited. Our study generated the first large-scale transcriptome dataset of tetraploid black locust for CK and IBA-treated cuttings. Additionally, the types and quantities of the genes expressed in CK and IBA-treated tetraploid black locust cuttings, along with their functions, classifications, and metabolic pathways, were revealed for the first time. In total, 10,181 and 14,924 DEGs were detected in CK and IBA-treated cuttings, respectively. Our study presents a global view of transcriptome dynamics and differential gene expression analysis for IBA-treated tetraploid black locust cuttings over four developmental stages and provides new insights into the fundamental mechanisms associated with auxin-induced adventitious rooting. Our data constitute a valuable resource for genomic investigations of AR formation in tetraploid black locust cuttings for improving rooting of difficult-to-root varieties.
Plant materials and IBA-treated adventitious root development
Cuttings of tetraploid black locust (Robinia pseudoacacia L.) were collected from a 1-year-old seed bearer at the nursery of Northwest Agriculture and Forestry University, Yangling, China. The preparation, experimental plan and IBA-treated cuttings were performed according to the methods described by Quan et al. [18, 19]. Sub-terminal parts of stems 15 cm in length and 10–12 mm in diameter were collected from the same cultivar. The basal 2.5 cm of each cutting was then placed for 4 h in water for the CK or in 5.4 mM IBA as the auxin treatment. The cuttings of tetraploid black locust were subsequently placed on a bench in a glasshouse, and 5 cm portions of the basal parts of cuttings were buried in sand. Cuttings were managed and randomly selected according to the method described by Quan et al. . The 2-cm-long basal stem regions of each cutting, where ARs develop, were cut and collected separately at 0, 15, 20, and 25 days after planting, and the harvested samples were designated control cuttings sampled at stages CI, CC, CRP and CAR for the CK and at stages II, IC, IRP and IAR for the IBA-treated cuttings. Ten cuttings randomly selected from each group were pooled. For each treatment and time point, three biological replicates, each consisting of cuttings from ten base stems, were collected, frozen immediately in liquid nitrogen, and stored at −80 °C prior to RNA extraction.
RNA separation, cDNA library preparation and Illumina sequencing
Using the manufacturer’s protocol (TIANGEN Biotech, Beijing), the RNA prep Pure Plant Kit was applied to isolate total RNA from each sample. To prevent contamination with genomic DNA, RNA samples were processed with RNase-free DNase I (Takara, Japan). We used an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) to obtain and quantify RNAs. Thereafter, we used ethidium bromide staining and denatured agarose gel electrophoresis, to inspect the completeness of RNAs. Subsequently, the analysis only included RNA integrity numbers ≥7.5, RNA 28S:18S ratios higher than one and RNA samples with A260/A280 ratios between 1.9 and 2.10.
Based on the manufacturer’s guidelines, Dynal oligo (dT)25 beads were used to separate poly(A)+ RNA from total RNA of tetraploid black locust in order to conduct Illumina sequencing. We purified the poly (A)+ RNA and then cleaved the mRNA into short pieces in a cleavage buffer. Using the random hexamer primer N6, SuperScript III reverse transcriptase and short pieces as templates, we synthesized the first-strand cDNAs. Next, RNase H, buffer, DNA polymerase I and dNTPs were used for second-strand cDNA synthesis. End repair of the double-stranded cDNAs was performed using T4 polynucleotide kinase, T4 DNA polymerase and DNA polymerase I Klenow fragment. Then, T4 DNA ligase was used to ligate the fragments to adapters. The QiaQuick PCR extraction kit was used to isolate the ligated fragments (200 ± 25 bp), which were then eluted in EB buffer and analyzed using agarose gel electrophoresis. Then, we chose appropriate segments as samples for PCR amplification. Eventually, we used an Illumina HiSeq™ 2000 to sequence the library at LC Sciences in Hangzhou, China. All technical steps were performed twice.
De novo assembly and functional annotation
We discarded empty reads, adapter sequences, and poor sequences from the clean raw reads (with unknown sequences (N) or sequences less than 25 bp). Next, short reads were processed in triplicate for de novo assembly of transcriptomes. Then, the triplicate samples were used to collect clean reads into non-repeating transcripts. To enhance the quality of the reads, we also deleted short sequences (<200 bp). Using the Kyoto Encyclopedia of Genes and Genomes (KEGGs) database with an E-value cut-off of 10−5 and the Clusters of Orthologous Groups of proteins (COGs), Swiss-Prot, Pfam and NCBI non-redundant protein (Nr) databases, we identified the most important sequences through BLAST searches and annotation. We used Blast2go software to perform functional annotation of Gene Ontology (GO) (www.geneontology.org).
Identifying differentially expressed genes
We applied RPKM (reads per kilobase per million reads), as the standard value of gene expression , to analyze differential gene expression. Using the IDEG6 web tool (http://telethon.bio.unipd.it/bioinfo/IDEG6_form/) , we performed statistical comparisons between the CK samples and the IBA-treated samples in terms of RPKM values. In many experiments, the importance of changes in gene expression was assessed based on the P-value, limiting the false discovery rate (FDR) at 0.01. In comparisons between libraries, genes showing fold-change >2 with P ≤0.05 were defined as differentially expressed. For up- and down-regulated unigenes, we used WEGO to compare their GO classifications .
Real-time quantitative PCR (RT-qPCR)
For RT-qPCR, we chose 21 candidate genes (Additional file 12) that showed differential expression to verify the Illumina sequencing data. Reverse transcription was performed using the PrimeScript RT reagent kit (Takara, Dalian, China) and 1 μg of RNA. The products of reverse transcription were diluted to 150 ng/μl in double-distilled water. Primer Premier 6.0 was used to generate antisense and sense primers. Using TaKaRa SYBR Premix Ex Taq II (Perfect Real Time), we performed RT-qPCR in a 20 μl reaction volume with a RC96 Real-Time PCR detection system (Rocipe Laboratories, Hercules, CA, USA). The 18S rRNA gene was used as the endogenous reference gene. After 2 min of predenaturation at 94 °C to measure fluorescence, RT-qPCR was performed with forty-five cycles of denaturation at 94 °C for 10, annealing at 60 °C for 15 s and elongation at 72 °C for 30 s. To ensure the absence of primer dimers, the melting curve was analyzed in the first cycle. The 2-△△Ct method  was used to confirm the relative expression of target genes. Each sample was analyzed in triplicate. We used Origin (Version 8.0) to analyze the data and then log-transformed all data obtained from the q-PCR analysis.
This study was supported by the National Natural Science Foundation of China (30972352) and the Innovation Projects in Institute of Soil and Water Conservation, Northwest A&F University (A315021398). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The raw sequencing data, mapped data, and data for visualization of the transcriptome sequence data for the four developmental stages of CK and IBA-treated tetraploid black locust cuttings have been deposited in the Gene Expression Omnibus (GEO) at the National Center for Biotechnology Information (NCBI) under accession number GSE80126.
Conceived and designed the experiments: ZZ, XY. Performed the experiments: JQ, SM. Analyzed the data: JQ. Contributed reagents/materials/analysis tools: ZZ. Wrote the paper: JQ. Interpretation of data and revised the paper: EG, SZ. All authors contributed to the preparation of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Wang XL, Zhao Z. Seasonal variation in rooting of the cuttings from tetraploid locust in relation to nutrients and endogenous plant hormones of the shoot. Turkish J Agr For. 2012;36:257–66.Google Scholar
- Wang XL, Zhao Z, Quan JE. Indole-3-butyric acid on rooting and endogenous plant hormones in tetraploid and diploid Robinia pseudoacacia hardwood cuttings. Phyton. 2011;80:93–100.Google Scholar
- Ahkami AH, Lischewski S, Haensch KT, Porfirova S, Hofmann J, Rolletschek H, et al. Molecular physiology of adventitious root formation in Petunia hybrida cuttings: involvement of wound response and primary metabolism. New Phytol. 2009;181:613–25.View ArticlePubMedGoogle Scholar
- da Costa CT, de Almeida MR, Ruedell CM, Schwambach J, Maraschin FS, Fett-Neto AG. When stress and development go hand in hand: main hormonal controls of adventitious rooting in cuttings. Front Plant Sci. 2013;4:133.View ArticlePubMedPubMed CentralGoogle Scholar
- Geiss G, Gutierrez L, Bellini C. Adventitious root formation: new insights and perspective. In: Beeckman T, editor. Root development-annual plant reviews. London: Wiley; 2009. p. 127–56.Google Scholar
- Tsanakas GF, Manioudaki ME, Economou AS, Kalaitzis P. De novo transcriptome analysis of petal senescence in Gardenia jasminoides Ellis. BMC Genomics. 2014;15:554.View ArticlePubMedPubMed CentralGoogle Scholar
- Cooper WC. Hormones in relation to root formation on stem cuttings. Plant Physiol. 1935;10:789–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Heloir MC, Kevers C, Hausman JF, Gaspar T. Changes in the concentrations of auxins and polyamines during rooting of in-vitro-propagated walnut shoots. Tree Physiol. 1996;16:515–9.View ArticlePubMedGoogle Scholar
- Ludwig-Müller J, Vertocnik A, Town CD. Analysis of indole-3-butyric acid-induced adventitious root formation on Arabidopsis stem segments. J Exp Bot. 2005;56:2095–105.View ArticlePubMedGoogle Scholar
- Rout GR. Effect of auxins on adventitious root development from single node cuttings of Camellia sinensis (L.) kuntze and associated biochemical changes. Plant Growth Regul. 2006;48:111–7.View ArticleGoogle Scholar
- Wei K, Wang L, Cheng H, Zhang C, Ma C, Zhang L, et al. Identification of genes involved in indole-3-butyric acid-induced adventitious root formation in nodal cuttings of Camellia sinensis (L.) by suppression subtractive hybridization. Gene. 2013;514:91–8.View ArticlePubMedGoogle Scholar
- Lindroth AM, Kvarnheden A, von Arnold S. Isolation of a PSTAIRE CDC2 cDNA from Pinus contorta and its expression during adventitious root development. Plant Physiol Biochem. 2001;39:107–14.View ArticleGoogle Scholar
- Lindroth AM, Saarikoski P, Flygh G, Clapham D, Gronroos R, Thelander M, et al. Two S-adenosylmethionine synthetase-encoding genes differentially expressed during adventitious root development in Pinus contorta. Plant Mol Biol. 2001;46:335–46.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Sedira M, Butler E, Gallagher T, Welander M. Verification of auxin-induced gene expression during adventitious rooting in rol B-transformed and untransformed apple Jork 9. Plant Sci. 2005;168:1193–8.View ArticleGoogle Scholar
- Ricci A, Rolli E, Dramis L, Diaz-Sala C. N, N’-bis-(2,3-Methylenedioxyphenyl)urea and N, N’-bis-(3,4-methylenedioxyphenyl)urea enhance adventitious rooting in Pinus radiata and affect expression of genes induced during adventitious rooting in the presence of exogenous auxin. Plant Sci. 2008;175:356–63.View ArticleGoogle Scholar
- Brinker M, van Zyl L, Liu W, Craig D, Sederoff RR, Clapham DH, et al. Microarray analyses of gene expression during adventitious root development in Pinus contorta. Plant Physiol. 2004;135:1526–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang S, Zhao Z, Zhang L, Zhou Q. Comparative proteomic analysis of tetraploid black locust (Robinia pseudoacacia L.) cuttings in different stages of adventitious root development. Trees. 2014;12:223–34.Google Scholar
- Quan J, Zhang S, Zhang C, Meng S, Zhao Z, Xu X. Molecular cloning, characterization and expression analysis of the SAMS gene during adventitious root development in IBA-induced tetraploid black locust. PLoS One. 2014;9:699–709.Google Scholar
- Quan J, Zhang C, Zhang S, Meng S, Zhao Z, Xu X. Molecular cloning and expression analysis of the MTN gene during adventitious root development in IBA-induced tetraploid black locust. Gene. 2014;155:140–50.View ArticleGoogle Scholar
- Wei K, Wang LY, Wu LY, Zhang CC, Li HL, Tan LQ, et al. Transcriptome analysis of indole-3-butyric acid-induced adventitious root formation in nodal cuttings of Camellia sinensis (L.). PLoS One. 2014;9:e107201.View ArticlePubMedPubMed CentralGoogle Scholar
- Ahkami A, Scholz U, Steuernagel B, Strickert M, Haensch KT, Druege U, et al. Comprehensive transcriptome analysis unravels the existence of crucial genes regulating primary metabolism during adventitious root formation in Petunia hybrida. PLoS One. 2014;9, e100997.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang H, Tong W, Feng L, Jiao Q, Long L, Fang R, et al. De novo transcriptome analysis of mulberry (Morus L.) under drought stress using RNA-seq technology. Bioorg Khim. 2014;40:458–67.PubMedGoogle Scholar
- Gupta S, Shi X, Lindquist IE, Devitt N, Mudge J, Rashotte AM. Transcriptome profiling of cytokinin and auxin regulation in tomato root. J Exp Bot. 2013;64:695–704.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng Z, Zhang C, Zhang Y, Hu T, Mu S, Li X, et al. Transcriptome sequencing and analysis of the fast growing shoots of moso bamboo (Phyllostachys edulis). PLoS One. 2013;8, e78944.View ArticlePubMedPubMed CentralGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18:1509–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Romualdi C, Bortoluzzi S, D’Alessi F, Danieli GA. IDEG6: a web tool for detection of DEGs in multiple tag sampling experiments. Physiol Genomics. 2003;12:159–62.View ArticlePubMedGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34:W293–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Dang ZH, Zheng LL, Wang J, Gao Z, Wu SB, Qi Z, et al. Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna. BMC Genomics. 2013;14:29.View ArticlePubMedPubMed CentralGoogle Scholar
- Rigal A, Yordanov YS, Perrone I, Karlberg A, Tisserant E, Bellini C, et al. The AINTEGUMENTA LIKE1 homeotic transcription factor PtAIL1 controls the formation of adventitious root primordia in poplar. Plant Physiol. 2012;160:1996–2006.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams D, Yang S. Ethylene biosynthesis: identification of 1-aminocyclopropane-1-carboxylic acid as an intermediate in the conversion of methionine to ethylene. Proc Natl Acad Sci U S A. 1979;76:170–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Lister R, Gregory BD, Ecker JR. Next is now: new technologies for sequencing of genomes, transcriptomes, and beyond. Curr Opin Plant Biol. 2009;12:107–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35 Suppl 2:W182–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Mochida K, Shinozaki K. Advances in omics and bioinformatics tools for systems analyses of plant functions. Plant Cell Physiol. 2011;52:2017–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Abu-Abied M, Szwerdszarf D, Mordehaev I, Levy A, Rogovoy O, Belausov E, et al. Microarray analysis revealed upregulation of nitrate reductase in juvenile cuttings of Eucalyptus grandis, which correlated with increased nitric oxide production and adventitious root formation. Plant J. 2012;71:787–99.View ArticlePubMedGoogle Scholar
- Dixon RA, Achnine L, Kota P, Liu CJ, Reddy MSS, Wang L. The phenylpropanoid pathway and plant defence-a genomics perspective. Mol Plant Pathol. 2002;3:371–90.View ArticlePubMedGoogle Scholar
- Davin LB, Jourdes M, Patten AM, Kim KW, Vassao DG, Lewis NG. Dissection of lignin macromolecular configuration and assembly: comparison to related biochemical processes in allyl/propenyl phenol and lignin biosynthesis. Nat Prod Rep. 2008;25:1015–90.View ArticlePubMedGoogle Scholar
- Lim CC, Liu JZ, Pua EC. Characterization of S-adenosylmethionine synthetase genes and its expression is associated with ethylene synthesis in mustard (Brassica juncea). Physiol Plantarum. 2002;116:522–30.View ArticleGoogle Scholar
- Villacorta-Martín C, Sánchez-García AB, Villanova J, Cano A, van de Rhee M, de Haan J, et al. Gene expression profiling during adventitious root formation in carnation stem cuttings. BMC Genomics. 2015;16(1):1.View ArticleGoogle Scholar
- Yu K, Xu Q, Da X, Guo F, Ding Y, Deng X. Transcriptome changes during fruit development and ripening of sweet orange (Citrus sinensis). BMC Genomics. 2012;13(1):1.View ArticleGoogle Scholar
- Albers E. Metabolic characteristics and importance of the universal methionine salvage pathway recycling methionine from 5’-methylthioadenosine. IUBMB Life. 2009;61:1132–42.View ArticlePubMedGoogle Scholar
- Wang Y, Tao X, Tang XM, Xiao L, Sun JL, Yan XF, et al. Comparative transcriptome analysis of tomato (Solanum lycopersicum) in response to exogenous abscisic acid. BMC Genomics. 2013;14:841.View ArticlePubMedPubMed CentralGoogle Scholar
- Perrot-Rechenmann C. Cellular responses to auxin: division versus expansion. Cold Spring Harb Perspect Biol. 2010;2:a001446.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim SB, Yu JG, Lee GH, Park YD. Characterization of Brassica rapa sadenosyl-L-methionine synthetase gene including its roles in biosynthesis pathway. Hortic Environ Biotechnol. 2012;53:57–65.View ArticleGoogle Scholar
- Li YH, Zou MH, Feng BH, Huang X, Zhang Z, Sun GM. Molecular cloning and characterization of the genes encoding an auxin efflux carrier and the auxin influx carriers associated with the adventitious root formation in mango (Mangifera indica L.) cotyledon segments. Plant Physiol Biochem. 2012;55:33–42.View ArticlePubMedGoogle Scholar
- Li SW, Shi RF, Leng Y, Zhou Y. Transcriptomic analysis reveals the gene expression profile that specifically responds to IBA during adventitious rooting in mung bean seedlings. BMC Genomics. 2016;17:43.View ArticlePubMedPubMed CentralGoogle Scholar
- Negi S, Sukumar P, Liu X, Cohen JD, Muday GK. Genetic dissection of the role of ethylene in regulating auxin-dependent lateral and adventitious root formation in tomato. Plant J. 2010;61:3–15.View ArticlePubMedGoogle Scholar
- Arjon J, Van H, Keith R. AtAGP30, an arabinogalactan-protein in the cell walls of the primary root, plays a role in root regeneration and seed germination. Plant J. 2003;36:256–70.View ArticleGoogle Scholar
- Laskowski MJ, Dreher KA, Gehring MA, Abel S, Gensler AL, Sussex IM. FQR1, a novel primary auxin-response gene, encodes a flavin mononucleotide-binding quinone reductase. Plant Physiol. 2002;128:578–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Gaedeke N, Klein M, Kolukisaoglu U, Forestier C, Müller A, Ansorge M, et al. The Arabidopsis thaliana ABC transporter AtMRP5 controls root development and stomata movement. EMBO J. 2001;20:1875–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang W, To JP, Cheng CY, Schaller GE, Kieber JJ. Type-A response regulators are required for proper root apical meristem function through post-transcriptional regulation of PIN auxin efflux carriers. Plant J. 2011;68:1–10.View ArticlePubMedGoogle Scholar
- To JPC, Haberer G, Ferreira FJ, Deruère J, Mason MG, Schaller GE. Type-A Arabidopsis response regulators are partially redundant negative regulators of cytokinin signaling. Plant Cell. 2004;16:658–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim S, Choi H, Ryu HJ, Ji HP, Kim MD, Kim SY. ARIA, an Arabidopsis arm repeat protein interacting with a transcriptional regulator of abscisic acid-responsive gene expression, is a novel abscisic acid signaling component. Plant Physiol. 2004;136:3639–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Coates JC, Laplaze L, Haseloff J. Armadillo-related proteins promote lateral root development in Arabidopsis. Proc Natl Acad Sci U S A. 2006;103:1621–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Negi S, Ivanchenko MG, Muday GK. Ethylene regulates lateral root formation and auxin transport in Arabidopsis thaliana. Plant J. 2008;55:175–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticlePubMedGoogle Scholar
- Conesa A, Gotz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.View ArticlePubMedGoogle Scholar