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

Comprehensive analysis of dynamic gene expression and investigation of the roles of hydrogen peroxide during adventitious rooting in poplar

Abstract

Background

Adventitious roots (ARs) are roots that are generated from nonrooting tissues. ARs are usually produced both during normal development and in response to stress conditions, such as flooding, nutrient deprivation, heavy metal stress and wounding. The ability of plants to form ARs is a key trait that enables plant propagation, especially for most tree species.

Results

Here, the kinetics of AR formation in a tissue culture of a hybrid variety of poplar were investigated. AR formation mainly occurred during the first 8 days and both pre- and newly- formed primordia contributed to AR formation in poplar by histological study. RNA-Seq-based transcriptome analysis was performed for stem bases collected at 0, 2, 4, 6 and 8 days after excision (DAE). Based on the data, the expression patterns of 8 phytohormone-related genes were investigated, and their influences on AR formation were considered. Subsequent gene expression cluster analysis showed a number of biological processes involved in AR formation. Among these biological pathways, genes involved in H2O2 homeostasis showed enrichment in one cluster that was highly upregulated from DAE0 to DAE8. Pharmacological assay confirmed that an appropriate content of H2O2 in stem bases could accelerate the formation of ARs in poplar.

Conclusions

Based on the results of this study, we were able to predict a regulatory network for 7 phytohormones that are involved in poplar AR formation. The influence of H2O2 on AR formation was also confirmed. These results enhance our understanding of the regulation of AR formation in tree species.

Background

Adventitious roots (ARs) are roots that are generated both in vitro and in field cultivation from nonrooting tissues such as stems, leaves and hypocotyls. ARs are usually produced both during normal development and in response to stress conditions, such as flooding, nutrient deprivation, heavy metal stress and wounding [1]. The formation of ARs is one of the most important biological processes in most trees and some herbaceous plants, which are vegetative propagated by semi-hardwood cuttings or by microcuttings if under an in vitro system. Therefore, in some plants, the ability to form ARs is the key trait enabling their propagation. During the past several decades, many studies have focused on determining the mechanisms involved in AR formation in herbaceous plants such as Arabidopsis thaliana, Oryza sativa and Petunia hybrids [2, 3]. Generally, three factors are attributed to controlling AR formation, including genetic, environmental and physiological factors (e.g., phytohormones, phenolic compounds and nitrogen oxide) [2, 4]. All 3 of these factors and their components interact with each other and form a complex network that regulates AR formation. However, the understanding of AR formation in tree species is relatively limited. Though the details of this biological process may be shared between herbaceous and woody plants to some extent, there are considerable differences between these two types of plants; for example, the ability of trees to form ARs usually declines with age; while there is absent of juvenility for most herbaceous plants [5,6,7]. Thus, an investigation of the mechanisms involved in AR formation in trees is necessary.

Phytohormones, also known as plant hormones, play important roles in plant growth and environmental responses. Phytohormones act as the most important endogenous factors controlling the process of AR formation in plants. Among the diverse phytohormones, auxin has been confirmed to play a central role in AR formation in many studies [8]. Generally, the tissues that produce ARs require high auxin levels, and this high content of auxin is usually dependent on polar auxin transport (PAT) [9, 10]. Interestingly, excessively high auxin levels have been observed to have an inhibitory role at later stages of AR formation [11]. Cytokinins, which compose another key phytohormone group, regulate the division and differentiation of plant cells. In the early stages of AR induction in carnation, accumulation of cytokinins increases, while the accumulation decreases during the middle stages [12]. In addition to these two phytohormones, abscisic acid (ABA), brassinosteroids (BRs), ethylene (ET), gibberellins (GAs), jasmonic acid (JA) and salicylic acid (SA) have also been found to play roles during AR formation. Generally, a large and complex network is formed among these phytohormones, and this network complicates the phytohormones’ influence on AR formation. For example, interactions between ET, GAs and ABA regulate the emergence and growth rate of ARs in deepwater rice [13]. GAs inhibit adventitious rooting in poplar and Arabidopsis by affecting auxin transport [14]. In petunia, JA acts as a positive regulator of AR formation in cuttings [15]. Since there are interactions among all these phytohormones, BRs and SA may also be involved in the regulation of AR formation. However, the understanding of how plants adjust and balance the influences of different phytohormones is limited. Thus, in summary, the roles of phytohormones in AR formation are complex, and knowledge regarding the regulatory network is still scarce.

Among the endogenous factors that affect AR formation, reactive oxygen species (ROS) also play important roles in regulating AR formation. Generally, ROS are considered signaling molecules and have been reported to regulate stomatal movements and biotic and abiotic stress responses as well as many aspects of plant development, including lateral root formation [16,17,18,19,20]. Recently, some studies have focused on the roles of ROS in AR formation. In mung bean, hydrogen peroxide (H2O2), one of the main types of ROS, has been predicted by RNA-Seq analysis to positively regulate AR formation by modulating a number of genes [21]. On olive adventitious rooting, an AOX protein (key enzyme on the alternative respiratory pathway) was found to putatively regulate the ROS homeostasis during the first stages of the adventitious rooting process [22]. Exogenous application of H2O2 has also been shown to increase the ability of marigold to form adventitious roots [23, 24]. Compared with that of other endogenous factors, the understanding of the role of ROS in AR formation is quite limited.

Poplar is a model tree species due to the small size of its genome, the early release of a high-quality draft genome and the availability of a highly efficient transgenic platform. Usually, poplar is considered as easy-rooting plant, however, different types of species still varied for this trait [25]. Recently, several studies have been performed on poplar AR formation. For example, transcripts involved in root and AR development were uncovered by expressed sequence tags sequencing [26]. The cytokinin type-B response regulator PtRR13 was also found to act as a negative regulator of adventitious root development in Populus [27]. Subsequently, a transcription factor in the AP2 family was also found to play a role during AR formation [28]. Several WUSCHEL-related HOMEOBOX (WOX) transcription factors have also been reported to regulate key developmental processes in AR formation in poplar [29, 30]. However, there is still insufficient understanding of poplar AR formation compared with that of AR formation in herbaceous model plants. Here, in order to increase our understanding to the regulation of poplar ARs formations and screen of genes or endogenous factors that are responsible for this biological processes, we investigated the kinetics of poplar AR development and the dynamic gene expression during AR formation by RNA-Seq based transcriptomic analysis in poplar. Based on these data, the expression patterns of phytohormone-related genes were investigated, and their influences on AR formation were considered. Subsequent gene expression cluster analysis showed that genes involved in H2O2 pathways may contribute to poplar AR formation. Finally, pharmacological assay confirmed the influence of H2O2 on AR formation in poplar. Thus, in this study, we provide novel information on the regulation of AR formation in tree species.

Results

Poplar AR formation in WPM

To investigate the kinetics of poplar AR formation, 30 microcuttings with 2–3 upper leaves and an apical bud (ca. 2–3 cm) of “NL895” were cultured in WPM in the absence of plant hormones and other growth regulators. The first visible AR was observed at 5 days after excision (DAE5), and > 80% of the microcuttings showed roots at DAE8. The number of ARs from DAE0 to DAE8 is shown in Fig. 1a. There were no visible ARs from DAE1 to DAE4, while 0.8 ± 0.3, 1.6 ± 0.4, 2.8 ± 0.7 and 3.2 ± 0.5 ARs could be observed at DAE5, DAE6, DAE7, and DAE8, respectively. After DAE8, the total number of ARs increased slightly, and the length of the ARs increased sharply (data not shown). The typical AR formation processes from DAE0 to DAE8 are illustrated in Fig. 2. According to the data showed in Figs. 1a and 2, it is easy to conclude that AR development mainly occurs inside stem cuttings at the first 4 days; while visible AR elongation mainly occurs from DAE4 to DAE8 for poplar cuttings growing in WPM.

Fig. 1
figure 1

The kinetics of poplar AR development from DAE0 to DAE8 (a), mean AR from DAE0 to DAE8; (b), mean AR primordium percent investigated by anatomy. At least the basal region of 15 plantlets was collected for each time point and at least 10 cross sections (> 50 cross sections were prepared and 10 were randomly selected) were used per cutting. The AR percentage was calculated considering the number of sections with AR primordia in each cutting in comparison with the total number of sections analyzed per cutting. Significant difference among multiple comparisons were labeled as different letters above the dot in a and bar in b (P = 0.05)

Fig. 2
figure 2

AR development from DAE0 to DAE8 in WPM

When using anatomical method to investigate poplar AR formation in WPM, we found there were different AR primordia rate among the 5 tested stages. The percent are DAE0 (30.4% ± 7.0%), DAE2 (35.2% ± 13.7%), DAE4 (45.7% ± 17.9%), DAE6 (60.0% ± 12.3%) and DAE8 (55.3% ± 11.0%) (Fig. 1b). With multiple comparison (P = 0.05), significant difference among DAE0, DAE2 and DAE4; while difference between DAE6 and DAE8 didn’t reach significant level. The typical sections from DAE0 to DAE8 are illustrated in Fig. 3. Clearly, a relative high rate of section could be observed AR primordia in stem cortex at DAE0 (30.4 ± 7.0%), it suggested that there were AR primordia formed before the preparation of cuttings. Additionally, AR primordia rate also increased after DAE0 (except DAE8, the slightly decreased percent might be attributed to difficulty in counting primordia due to large number of mature AR), it suggested that poplar cuttings also could form new AR primordia during its development. Therefore, both pre- and newly- formed AR primordia exist in poplar cuttings.

Fig. 3
figure 3

Anatomy of AR formation from in WPM from DAE0 to DAE8 Bar = 20 μm

Time course of gene expression during AR formation

To determine the gene expression profiles during AR formation, the 0.5-cm cutting bases inserted in the WPM (DAE0, DAE2, DAE4, DAE6 and DAE8) were collected for each time point. In total, 15 RNA-Seq libraries were created (three biological samples for each time point) using the collected samples. A total of approximately 284 million 150-bp paired-end (PE) reads were obtained (Table 1). After quality control and filtering, 84.0 G of clean data remained, and the average amount of data in the samples was 5.60 ± 0.61 G (Table 1). The smallest sample was DAE4–1 with 4.83 G, while the largest sample was DAE2–2 with 6.29 G of clean data (Table 1). The clean reads constituted 98.5% of all raw reads. The high quality and sufficient quantity of the RNA-Seq data supported further gene expression analysis.

Table 1 Statistical summary of RNA-Seq libraries

All the clean reads were mapped onto the reference genome, and the gene expression numbers were calculated by TopHat and Cufflinks [31, 32]. To identify differentially expressed genes, gene expression at the four later time points, DAE2, DAE4, DAE6 and DAE8, was compared with that at DAE0. Moreover, consecutive time points were also compared, including DAE2 and DAE4, DAE4 and DAE6, and DAE6 and DAE8. Therefore, a total of 7 paired comparisons were conducted. According to the criteria for differentially expressed genes mentioned in the Materials and Methods section, there were 7483 nonredundant differentially expressed genes (Additional file 1: Table S2). When these differentially expressed genes were subjected to cluster analysis, samples from the same time point clearly tended to cluster together (Additional file 2: Figure S1). These results suggested that the samples collected were of high quality. Moreover, DAE8 and DAE0 showed less distance than the other time points and could be grouped in one subclade, while DAE2, DAE4, and DAE6 were in another subclade. The existence of 2 expression subclades is consistent with the above result regarding the kinetics of poplar AR formation. AR formation mainly occurred from DAE2 to DAE6, while the stem base returned to normal growth by DAE8.

When the differentially expressed genes for all the selected comparisons were examined, there were 4882, 4774, 4354, 3129, 940, 1173 and 2155 genes showing differential expression between DAE0 and DAE2, DAE0 and DAE4, DAE0 and DAE6, DAE0 and DAE8, DAE2 and DAE4, DAE4 and DAE6, and DAE6 and DAE8, respectively (Additional file 3: Figure S2). The highest number of differentially expressed genes was found between DAE0 and DAE2, and the lowest number of differentially expressed genes was found between DAE2 and DAE4. Of the 4 consecutive comparisons (DAE0 and DAE2, DAE2 and DAE4, DAE4 and DAE6, and DAE0 and DAE8), DAE0 and DAE2 had the highest number of differentially expressed genes, while DAE2 and DAE4 had the lowest number. Since consecutive comparisons are likely reliable reflections of the dynamic changes in gene expression during AR development, this result might suggest that the period of time from 0 to 2 DAE is the most important stage for AR formation.

To validate the gene expression data from RNA-Seq analysis, 10 genes showing patterns of differential expression among the 5 time points were randomly selected for RT-qPCR assay. The relative fold changes between each time point and the control samples (DAE0) were calculated by setting expression at DAE0 equal to 1 (Fig. 4). To compare the gene expression profiles of the RNA-Seq and RT-qPCR data, a correlation analysis was conducted. The two types of expression data were highly consistent, with an R2 value equal to 0.86 and a p-value less than 0.0001. This high R2 value clearly suggested that the expression results were very reliable for further analyses.

Fig. 4
figure 4

Gene expression analysis of the 10 selected genes by RT-qPCR assay. The FPKM values of 9 genes, excluding Potri.001G404100, can be found in Additional file 1: Table S2. The FPKM values of Potri.001G404100 for DAE0, DAE2, DAE4, DAE6, and DAE8 are 23.1, 27.9, 23.6, 23.2 and 13.2, respectively. Because these genes were not classified as differentially expressed genes, their FPKM values are not in Additional file 1: Table S2. Two stars represent the relative expression fold is different between the tested time point and DAE0 under the significant level of P value equal to 0.01; while one star represents the relative expression fold is different between the tested time point and DAE0 under the significant level of P value equal to 0.05

Investigation of the expression patterns of phytohormone-related genes

Phytohormones are regarded as key regulators of AR formation. During AR formation, different hormones interact with each other and form a large and complex network that regulates this dynamic biological process. To gain insight into the roles of phytohormones in AR formation, we investigated phytohormone-related genes in poplar with homology searching (see Materials and Methods). In total, 807 genes were predicted to be involved in biosynthesis, transport, metabolism and reception of the 8 phytohormones (Additional file 4: Table S3). These 8 phytohormones included ABA, auxin, BRs, cytokinins, ET, GAs, JA and SA, and the corresponding numbers of related genes were 75, 201, 291, 29, 42, 58, 82, and 29, respectively (Table 2). No genes were predicted to be involved in BR transport, ET transport and metabolism, GA transport and metabolism, JA transport and SA reception (Table 2). When the expression patterns for these phytohormone-related genes were examined, 180 genes showed differential expression during AR formation; their normalized expression values are shown in Additional file 5: Table S4.

Table 2 The number of phytohormone-related genes in poplar

We therefore summarized the dynamic gene expression that occurred over the 5 AR time points. Among the phytohormone biosynthesis genes, clear patterns could be observed for 6 hormones (Fig. 5). No auxin and SA biosynthesis genes showed differential expression during AR formation. The expression of ABA biosynthesis genes was reduced at DAE2, DAE4 and DAE6 and was then slightly elevated to normal levels by DAE8. The expression of BR biosynthesis genes increased sharply to a very high level by DAE2 and then slightly decreased at subsequent time points; however, the expression of these genes from DAE2 to DAE8 remained higher than at DAE0. In contrast, the expression of ET and GA biosynthesis genes showed an opposite trend. Gene expression in this category decreased sharply to a very low level by DAE2 and then slightly increased at subsequent time points. Cytokinin biosynthesis genes showed the highest and lowest expression at DAE4 and DAE8, respectively, while their expression remained at normal levels at the other 3 time points. The expression of JA biosynthesis genes was reduced at all stages compared to that at DAE0.

Fig. 5
figure 5

Boxplot of FPKM values for differentially expressed genes related to phytohormone biosynthesis. FPKM represents fragments per kilobase of transcript per million fragments mapped. Each boxplot consists of the Q1 (1/4 quartiles, the bottom of the box), Q3 (3/4 quartiles, the top bottom of the box), the median Q2 (the line inside into the box), Q3 + 1.5 IQR (interquartile range, Q3-Q1, the highest line outside the box) and Q1–1.5 IQR (the lowest line outside the box). Values larger than Q3 + 1.5 IQR or less than Q1–1.5 IQR were considered as outliers and they were plotted as dots outside of the Q3 + 1.5 IQR or Q1–1.5 IQR. If Q3 + 1.5 IQR is less than Q3 or Q1–1.5 IQR is less than Q1, their corresponding lines will not be shown. ABA represents abscisic acid, BR represents brassinosteroids, ET represents ethylene, GA represents gibberellins and JA represents jasmonic acid

Among the phytohormone transport genes, only 3 hormones, ABA, auxin and cytokinins, had data that enabled an investigation of their expression patterns (Additional file 6: Figure S3). The average expression of all 3 phytohormone transport genes was higher at DAE0 than other stages. The plant material used for the stage of DAE0 also represents a tissue of a plant at normal growth status. Thus, the highest average expression for those phytohormone transport genes at DAE0 suggested that transport of ABA, auxin and cytokinins is required for poplar normal growth. The transport of the 3 phytohormones might also support the requirement of the initiation of AR development in poplar. At the stage of DAE2, the average expression of these phytohormone transport genes decreased sharply might be attributed to the excision of shootings from mother plants. The steadily increase of the average expression of these phytohormone transport genes from DAE2 to DAE8 suggested that ABA, auxin and cytokinins were involved more and more for poplar AR development. Since the processes of biosynthesis and transport are the source of endogenous phytohormones in plants, the dynamic expression patterns of these 7 phytohormone biosynthesis and transport genes might reflect their roles during AR formation.

Among the phytohormone metabolism genes, those related to BRs, cytokinins and SA had elevated expression from DAE2 to DAE8 compared to their expression at DAE0 (Additional file 7: Figure S4). Among the phytohormone receptor genes, those related to 6 hormones had data that enabled an investigation of their expression patterns (Additional file 8: Figure S5). BR and cytokinin receptor genes showed similar expression patterns, decreasing sharply by DAE2 and then increasing at subsequent stages; however, their expression at DAE8 remained lower than at DAE0. ABA and GA receptor genes also shared similar expression patterns, which demonstrated trends opposite those of the BR and cytokinin receptor genes. Auxin and ET receptor genes also showed opposite patterns to those of the BR and cytokinin receptor genes; the expression of auxin receptor genes first decreased at DAE2 and DAE4 and then increased by DAE6 and DAE8 to a higher level than at DAE0. Usually, both negative and positive regulation exists for metabolism and receptor genes when they respond to the corresponding phytohormones, therefore, we can’t conclude roles for these genes during poplar AR formation.

Functional annotation of genes differentially expressed during AR formation by gene clustering

To comprehensively investigate other factors affecting AR formation in poplar, we classified the differentially expressed genes by clustering them based on their dynamic expression patterns. With the K-means algorithm, 12 clusters were formed for these 7483 differentially expressed genes. There were 551, 743, 304, 479, 457, 454, 854, 1567, 305, 896, 638 and 235 genes in clusters 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 and 12, respectively (Additional file 1: Table S2). The expression profiles for each cluster are shown in Fig. 6. The genes in clusters 1, 4, 5, 6 and 7 showed upregulation from DAE2 to DAE8 compared to their expression at DAE0, while the genes in clusters 8, 9, 10 and 12 showed downregulation at all time points compared to their expression at DAE0 (Fig. 6). However, all 12 clusters had unique expression patterns, with the highest and lowest mean expression values at different AR formation stages. The large number of clusters and the diverse expression patterns suggest that AR formation in poplar is a complex biological process.

Fig. 6
figure 6

Clustering of genes differentially expressed during AR formation. FPKM represents fragments per kilobase of transcript per million fragments mapped

The grouping of genes in the same cluster usually indicates that those genes are involved in similar biological processes. Therefore, functional enrichment analyses were performed for each gene cluster. Ten of the 12 clusters showed enrichment of at least one biological function (Additional file 9: Table S5); only clusters 3 and 9 failed to show any enrichment. The most enriched biological process for cluster 1 was “photosynthesis”, and the most enriched process for cluster 5 was “photosynthesis, light harvesting”. The gene expression patterns in these 2 clusters showed steady upregulation (Fig. 6). These results might suggest that photosynthesis took place during AR formation and that photosynthesis would be beneficial for poplar AR formation. Additionally, “chitin catabolic process”, “response to biotic stimulus”, and “defense response” were enriched in clusters 7 and 11. Previously, genes with these functions have been reported to be involved in plant responses to disease [33, 34]. Both clusters 7 and 11 showed peak expression at DAE2 before the expression gradually decreased to normal levels; this result might suggest that some disease response genes also play roles during AR formation. Since KEGG analyses more directly indicate gene function and their involved pathways, we carefully examined results from KEGG analysis for each cluster. In total, five KEGG pathways were enriched in 3 gene clusters (Table 3). “Tubulin beta” was enriched in cluster 2; this result is similar to the result described above for cluster 2 [35, 36]. Three KEGG pathways, “peroxidase”, “aquaporin PIP” and “aconitate hydratase 1/homoaconitase”, were enriched in cluster 6. “Large subunit ribosomal protein L39e” was enriched in cluster 12. Among these 4 KEGG pathways, “peroxidase” is usually responsible for ROS homeostasis, and ROS usually act as signaling molecules during plant development. Therefore, genes encoding peroxidase could possibly be involved in AR formation in poplar.

Table 3 Functional enrichment for gene clusters by KEGG analysis

Investigation of the roles of H2O2 in poplar AR formation

To test the above hypothesis, we mainly focused on genes in the peroxidase pathway. A total of 7 differentially expressed genes encoding peroxidase were grouped in cluster 6, and most of these genes showed increased expression after DAE0 and peak expression at DAE6 (Table 3). The gene Potri.005G195600 was selected for RT-qPCR assay. The results were very consistent with the RNA-Seq results (Fig. 4). After RT-qPCR analysis, the H2O2 content in the basal stems of tissue culture plantlets was tested. H2O2 accumulation clearly increased during the entire process of AR formation. The H2O2 content was 53.5 ± 2.0, 86.8 ± 5.2, 97.3 ± 5.3, 110.2 ± 6.7 and 374.8 ± 41.3 mmol/g.prot for DAE0, DAE2, DAE4, DAE6, and DAE8, respectively (Fig. 7). Through variance analysis and multiple comparisons, the content of H2O2 at DAE2, DAE4 and DAE6 was revealed to be significantly higher than at DAE0 (P < 0.05), while the H2O2 content at these 3 time points was significantly lower than at DAE8 (P < 0.05). Therefore, we concluded that the H2O2 content showed three levels during AR formation: the H2O2 content was at a baseline level at the beginning time point; it increased and remained at a relatively high level from DAE2 to DAE6, and it then increased sharply at DAE8.

Fig. 7
figure 7

H2O2 accumulation during AR formation. Five to 15 stem bases were harvested to form one biological samples (at least 0.1 g) and 6 biological samples were prepared for one treatment. The multiple comparisons was performed under the significant level of P value equal to 0.05. The error bar represents standard deviation for each treatment

To further understand the roles played by H2O2 during poplar AR formation, we performed hydroponic culture by using semi-hardwood cuttings of poplar “NL895” as plant material. In the preliminary experiment, five cuttings with uniform growth status were grown in hydroponic solution supplemented with 0, 500 or 1000 μM H2O2. Cuttings grown in 500 μM H2O2 produced markedly more ARs at DAE11 than cuttings grown in the other H2O2 concentrations (Additional file 10: Figure S6). Then, an average of 30 cuttings with uniform growth status were grown in each of 5 hydroponic solutions (Note, only 15 data were finally obtained for 1000 μM H2O2 treatment), and the number of ARs for each cutting was recorded every day. The number of ARs was recorded after 11, 15, 19 and 23 days of growth in hydroponic solution supplemented with 0, 250, 500, 750 or 1000 μM H2O2. The results showed that the average number of ARs was significantly affected by the concentration of H2O2 (Fig. 8). From days 11 to 19, the average AR number for the cuttings grown in the solution with 500 μM H2O2 was higher than that for the cuttings grown in the other solutions (P < 0.05), while the differences between the average AR numbers seemed to decline after day 19. Additionally, the average AR numbers of the cuttings grown in the solutions with 750 and 1000 μM H2O2 were lower than those of the cuttings grown in the untreated control solution (0 μM). Therefore, three findings were evident from these data. First, treatment with exogenous H2O2 could affect AR formation; second, AR developed faster in the 500 μM solution than in the other solutions; and third, higher concentrations of exogenous H2O2 slowed AR formation in poplar. All these results also suggest that appropriate levels of H2O2 production can accelerate AR formation in poplar. To validate the above results, poplar “NL895” cuttings were also grown in hydroponic solution with 0.5 mM DMTU, which is a chemical H2O2 scavenger. By day 8, the H2O2 content in the stem base of cuttings that were grown in DMTU solution was 142.7 ± 7.3 (n = 30), while that in the stem base of untreated controls was 367.7 ± 3.2 (n = 30). This revealed that DMTU reduced the accumulation of H2O2 (T test, P < 0.05). The average number of ARs in the cuttings grown in hydroponic solution with 0.5 mM DMTU was significantly lower than that in the cuttings grown in control solution at DAE12 (Fig. 9). The average number of ARs in DMTU was 0.8 ± 1.1, while that in the untreated controls was 6.6 ± 1.3. This result supports the finding that appropriate levels of H2O2 production can accelerate AR formation in poplar. Though AR developed faster in solutions with appropriate concentrations of H2O2, the influence of H2O2 on AR length seemed to follow a different trend. At day 23, the average AR lengths for the cuttings grown in the solutions supplemented with 0, 250, 500, 750 and 1000 μM H2O2 were 2.2 ± 1.5, 2.1 ± 0.9, 0.5 ± 0.6, 0.6 ± 0.5 and 0.7 ± 0.8 cm, respectively (Fig. 10). The solutions supplemented with 500, 750 and 1000 μM H2O2 had cuttings with significantly lower average AR lengths than the 2 solutions supplemented lower concentrations (P < 0.05). Although the average AR length was not significantly different between cuttings grown in solutions supplemented with 0 and 250 μM H2O2, the absolute number of ARs was lower in cuttings from the solution with 250 μM H2O2. This result might suggest that H2O2 negatively regulates AR elongation in poplar.

Fig. 8
figure 8

AR numbers during hydroponic culture of poplar cuttings with different concentrations of supplemented H2O2. There were 15 plants for 1000 μM H2O2 treatment; while 30 plants for all other treatments. AR number was presented by mean value of each treatment. The multiple comparisons was performed under the significant level of P value equal to 0.05. The error bar represents standard deviation for each treatment

Fig. 9
figure 9

Aspect of poplar cuttings 12 days after establishment in hydroponic culture with or without DMTU

Fig. 10
figure 10

AR length of poplar cuttings 12 days after establishment under hydroponic culture with different concentrations of H2O2.. The plant number for each treatment were the same as Fig. 8. AR length was presented by mean value of each treatment. The multiple comparisons was performed under the significant level of P value equal to 0.05. The error bar represents standard deviation for each treatment

Discussion

Complex biological processes occur in poplar AR formation

Poplars is considered as easy-rooting species and there is still complex biological processes during AR formation [1]. In previous studies, the development of poplar AR was classified into 3 stages: primordium organization, primordium differentiation and root elongation [4, 7]. As shown in Figs. 1b and 3, our results were partially consistent with this conclusion. Our data showed that primordium organization and differentiation also happened before the excision of cuttings. Both pre- and newly- formed primordia contribute to AR formation for poplar “NL895”. This result is consistent with previous reports [37,38,39]. Some cuttings produced ARs from DAE5 to DAE8, while others did not produce ARs over the same time points. Newly-formed AR primordium could also be observed at all 5 stages. This result suggested that there were no clear boundaries between stages and that the 3 processes (primordium organization, primordium differentiation and root elongation) might exist simultaneously at the same time point within one stem cutting after DAE4. However, the time points of the major developmental events were clear; primordium organization and differentiation mainly occurred from DAE0 to DAE4, and root elongation mainly occurred from DAE4 to DAE8. Therefore, we classified AR formation in poplar “NL895” into 2 stages, early and late. The early stage is from DAE0 to DAE4 and AR development occurs at the inside of stems; while the late stage is from DAE4 to DAE8 and AR development occurs at the outside of stems. According to the above analysis, we can conclude that 1) the key stage for AR formation in WPM is in the first 8 days after excision; 2) both pre- and newly- formed AR primordia contribute to AR formation for poplar; 3) the production of newly formed AR primordia occurs at most of the whole stage from DAE0 to DAE8. Therefore, our data demonstrated that the development of poplar ARs is a dynamic and complex biological process.

Multiple factors show influence on poplar AR formation

Phytohormones are the most important factors affecting AR formation. Among diverse phytohormones, auxin has been confirmed to play a central role in AR formation in many studies [8,9,10]. Since transcriptomic analysis has been proved to be a powerful strategy for uncover underlying mechanisms for a given biological process, many studies have been focused on root or AR formation by using this method [22, 25, 26]. In this study, we also uncovered some underlying influencing factors for poplar AR regulation. According to data illustrated in Fig. 5 and Additional file 6: Figure S3, Additional file 7: Figure S4 and Additional file 8: Figure S5, there were no genes related to auxin biosynthesis showing differential expression during AR formation in this study; while 50 genes related to auxin transport showing differential expression. Therefore, this data suggested that there was no local biosynthesis occurring during poplar AR formation and the accumulation of auxin was mainly attributed to transport. This result is similar to previous reports [10]. According to the expression profiles of the phytohormone biosynthesis- and transport-related genes and the accumulation characteristics (local biosynthesis or transport from other sites), we were able to deduce the roles of 7 investigated phytohormones, including ABA, auxin, BRs, cytokinins, ET, GAs and JA. These 7 phytohormones could clearly be classified into 5 groups. The first group included only one phytohormone, ABA. Local ABA biosynthesis decreased from DAE2 to DAE6 and then increased by DAE8. This pattern suggested that ABA negatively regulates primordium organization and primordium differentiation but positively regulates root elongation during AR formation. The second group included one phytohormones, BRs. This phytohormone positively regulated all 3 AR developmental stages, but their effects decreased late in AR development. The third group included one phytohormone, cytokinins. The roles of cytokinins in AR regulation are complex compared to those of other phytohormones. The expression of cytokinin biosynthesis genes had very large variation from DAE0 to DAE8; the expression levels were slightly decreased by DAE2, sharply elevated by DAE4, and then decreased by DAE6 and DAE8. This complex expression profile suggested a complex role in which cytokinins act as positive regulators at some stages and as negative regulators at the other stages. The fourth group included 3 phytohormones, ET, GAs and JA. All of these phytohormones exerted negative effects on AR formation at all 3 stages. The fifth included one phytohormones, auxin. Its accumulation mainly depended on transport and it positively regulates AR formation at both early and late stages.

Besides phytohormones, other factors could be deduced to show influence on poplar AR formation according to the RNA-Seq results. According to the data listed in Additional file 9: Table S5, some of the results agreed strongly with those of previous studies. For example, the “microtubule” function was enriched in cluster 2 [35, 36], and “response to hormone” and “response to wounding” were enriched in clusters 10 and 11 [33, 34], respectively. Moreover, some interesting and novel information on poplar AR formation was also revealed. The photosynthesis and carbon fixation pathways enriched in gene expression cluster 1 and the expression of this cluster was steadily increased from DAE0 to DAE8 (Fig. 6 and Additional file 9: Table S5). This result might suggest that sugar-signaling pathway is involved in poplar AR formation. Similarly, wounding responsive, biotic stimulus and defense regulation pathways showed enrichment in cluster 11 and the expression of this cluster was sharply increased at DAE0 (Fig. 6 and Additional file 9: Table S5). The related gene expression was abundant at early stage of poplar AR formation might suggest that the excision of poplar cuttings activated wounding signal and subsequently inducted AR initiation. This prediction is very similar to previous reports [26]. According to result illustrated in Fig. 6 and Additional file 9: Table S5, more other factors that show potential involvement in poplar AR formation also could be found, such as secondary metabolites, energy, nitric oxide and redox. All these results are similar to a number of previous studies [23, 24, 40]. Therefore, we could conclude that multiple factors show influence on poplar AR formation.

H2O2 could affect poplar AR formation

In this study, we investigated the influence of H2O2 on poplar AR formation by using pharmacological tests. The accumulation of H2O2 was increased after DAE0 (Fig. 7). However, we also could observe that H2O2 content at DAE8 was also significantly higher than DAE2 to DAE6 (P < 0.05). It suggested that there were 3 levels of H2O2 content from DAE0 to DAE8, with the lowest, moderate and highest at DAE0, DAE2 to DAE6 and DAE8, respectively. Since there were visible ARs at DAE8 and we couldn’t completely exclude the influence of newly formed ARs on H2O2 content assay, thus, the much higher H2O2 accumulation at DAE8 might be partially attributed to the newly formed ARs (poplar roots generated higher H2O2 content than stems, data not shown). However, the accumulation of H2O2 was still higher at DAE2 to DAE6 than at DAE0 (from 53.5 ± 2.0 to 110.2 ± 6.7 mmol/g.prot), it revealed that accumulation of H2O2 was indeed increased during the early stage of poplar AR formation. In previous studies, H2O2 was also found to be involved in AR or lateral root formation, such as in Mung bean [21], Marigold [23, 24] and Alfalfa [17]. Comparing with these studies, our result is highly consistent with previous reports.

The data clearly revealed that an appropriate content of H2O2 could accelerate the formation of AR, while higher H2O2 content could slow AR formation and reduce root elongation. H2O2 is a type of reactive oxygen species that has various roles in many biological processes, including responses to biotic and abiotic stresses and plant growth [40]. However, the underlying mechanisms of its influence on poplar AR formation are unclear. Some previous studies have reported a relationship between H2O2 and a number of phytohormones [41]. Here, we also predicted the roles of 7 phytohormones and H2O2 during poplar AR formation. Since the accumulation of H2O2 increased steadily during the whole AR formation process and since two phytohormones, auxin and BRs, positively regulated AR formation, these results might suggest that there was direct or indirect regulation among H2O2, auxin and BRs. In A. thaliana, ROS act as downstream regulators of auxin in the promotion of root hair development [42]. Based on our data, we could also speculate that H2O2 was regulated by auxin and BRs during AR formation. An appropriate H2O2 content promoted AR formation. However, there are 3 stages for poplar AR formation, primordium organization, primordium differentiation and elongating roots. The H2O2 concentrations at these 3 stages were significantly different (P < 0.05). Thus, we cannot conclude whether H2O2 plays a role at all 3 stages or only at a special stage. Therefore, more research should be performed to comprehensively determine the influence of H2O2 on AR formation in poplar.

Conclusions

Here, we were able to provide a global view of the roles of phytohormones and H2O2 during poplar AR formation by using RNA-Seq and pharmacological testing. The results showed a large and complex regulatory network among diverse phytohormones and H2O2. The knowledge gained from this study will help us to understand the underlying mechanisms of AR formation.

Methods

Plant materials, growth conditions and treatments

A hybrid variety of poplar “NL895” (Populus. euramericana), which forms ARs easily both in tissue culture and field cultivation, was used as plant material in this study. This variety was obtained from the tree germplasm repository of Huazhong Agricultural University (HZAU) under legal permission. For the investigation of the kinetics of poplar AR formation and for the preparation of RNA for high-throughput sequencing, tissue culture plantlets of “NL895” were used as plant material. Woody plants with different juvenile periods exhibit very different AR-forming abilities; therefore, two strategies were employed to avoid this influence. First, all microcuttings with undifferentiated ARs were obtained from apical buds that were generated from the first subculture of the differentiated poplar explants. Second, only microcuttings with 2–3 upper leaves and an apical bud (ca. 2–3 cm) were propagated in woody plant medium (WPM) for further analyses. The WPM was prepared according to previous report [43]. Additionally, all microcuttings that were inserted similar deep into the WPM (ca. 0.5 cm). All the tissue culture experiments were conducted at 25 °C with a 16/8 h photoperiod. To prepare samples for RNA extraction, the 0.5 cm cutting base was collected from each plantlet. Ten to 15 cutting bases were used for one biological sample. Three biological samples were prepared for each time point. Samples were collected at 5 time points: 0, 2, 4, 6 and 8 days after the cuttings were inserted into WPM. These time points were termed DAE0, DAE2, DAE4, DAE6 and DAE8, respectively. At DAE0, the cuttings had not experienced any treatment but were simply excised from the explants.

Test the influence of H2O2 and N,N′-dimethylthiourea (DMTU) on AR formation

An hydroponic system was conducted to test the influence of H2O2 and DMTU on poplar AR formation. For the preparation of poplar semi-hardwood cuttings used in hydroponic propagation, young plantlets growing under in vitro conditions, with 6 to 10 cm in height, were acclimatized and then transplanted into 1-L soil pots. These plantlets were then grown for approximately 1 month until the plants had 10–15 leaves. Cuttings with 2–4 upper leaves and an apical bud were taken from the 1-month-old poplar plants. The hydroponic solution was prepared with double-sterilized water.

For the test of the influence of H2O2 on poplar AR formation, cuttings were inserted into the hydroponic solution with 0, 250, 500, 750 or 1000 μM H2O2. For the test of the influence of DMTU on poplar AR formation, cuttings were inserted into the hydroponic solution with 0.5 mM or without DMTU.

The number of visible ARs was recorded for one month for both experiments. For data collection and analysis, the AR length and number were recorded for each individual cutting and then average values were generated for each treatment. All hydroponic was conducted at 25 °C with a 16/8 h photoperiod. The hydroponic solution was changed every 2 days. An average of 30 cuttings were grown for each treatment during hydroponic propagation supplemented with H2O2 and DMTU, except the 1000 μM H2O2 treatment had 15 final collected data.

Anatomical and microscopic analysis

For histological analysis, the same method and growth condition as RNA-Seq preparation were employed for collecting the 0.5 cm cutting base. At least 15 cutting bases were collected for each time points. The collected tissues were fixed in FAA (50% ethanol, 5% glacial acetic and 4% formaldehyde) and then treated with protocol described in a previous study [34]. Sections (15 μm) were made with rotary microtome (Leica RM2255) and stained with Safranin (1%) and fast-green (1%). At least 10 cross sections for one 0.5 cm cutting base were randomly selected for AR primordium counting. All the sections were observed and photographed with an Olympus BX51 microscope. The software FV10-ASW 3.1 Viewer (Olympus Support) was used to export photographs.

H2O2 content assay

For H2O2 content assay, we followed protocol described in Bio-101 with slightly modification [44]. Briefly, the stem bases that were inserted into WPM or hydroponic solution were collected. Five to 15 stem bases were harvested to form one biological samples (at least 0.1 g) and 6 biological samples were prepared for one treatment. The sample (0.1 g) was ground to powder in liquid nitrogen and suspended with 4 °C phosphate buffer with thoroughly vortex. The mixture was centrifuged at 12000 rpm at 4 °C. The supernatant was used to test concentration of H2O2 with Amplex Red hydrogen peroxide/peroxidase assay kit (Molecular Probes, USA, CAS: A22188) according to its introduction. Six technique replications were set for each biological sample.

Isolation of total RNA and whole transcriptome analysis

Total RNA was isolated using an RNeasy Plant Mini Kit according to the manufacturer’s instructions (DP432, TIANGEN Biotech (Beijing) Co., Ltd., Beijing, China). The quality of all RNA samples was examined by performing agarose gel electrophoresis and examination of the A260/280 wavelength ratio. The RNA integrity number (RIN) of all samples were also tested by an Agilent Bioanalyzer 2100. All samples had sufficient quality were used for RNA-Seq library construction according to a previously described protocol [45]. The constructed libraries were sent for sequencing using a HiSeq 2500 platform.

Bioinformatic analysis

The raw RNA-Seq reads were quality trimmed and quality filtered using Trimmomatic software with the default parameter settings [46]. The P. trichocarpa genome version 3.0 was used as the reference genome in this study [47]. Clean reads from each sample were mapped onto the reference genome using TopHat software with the default parameters [31, 32]. Values of fragments per kilobase million (FPKM) were obtained by the software Cufflinks, and the differential expression was analyzed by Cuffdiff [31, 32]. A false discovery rate (FDR) less than 0.05 and a fold change (FC) greater than 2.0 were set as the thresholds for differentially expressed genes. Differentially expressed genes with FPKM values less than 10 at all 5 time points were excluded from further analysis.

For the annotation of phytohormone-related genes in poplar, a homology searching strategy was employed. First, phytohormone-related genes in A. thaliana were obtained from the Arabidopsis Hormone Database (http://ahd.cbi.pku.edu.cn/) [48], and genes with functions related to biosynthesis, transport, metabolism and receptors were kept. Second, a BLAST search of poplar genes that were homologous with the retained Arabidopsis hormone genes was conducted by using the command “blastp” in a local BLAST database (ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/). The predicted proteins from the P. trichocarpa genome version 3.0 were used as queries. BLAST hits with scores greater than 250 (E-values of approximately e− 60) were kept. The biological processes of the retained poplar query proteins were annotated based on the functions of their top hit counterparts in the Arabidopsis Hormone Database.

To cluster the differentially expressed genes, the z-score (zero-mean normalization) of each FPKM value was determined with the “scale” command in R software (https://www.r-project.org/). The z-scores of all the differentially expressed genes were then used to form clusters with the “k-means” command in R software. The number of clusters was determined as previously described [49]. Functional annotation of the differentially expressed genes in each cluster was performed online with PlantGSEA (Plant GeneSet Enrichment Analysis Toolkit, http://structuralbiology.cau.edu.cn/PlantGSEA/) [50].

Gene expression analysis

Ten differentially expressed genes were selected for RT-qPCR assay. The primers for these 10 randomly selected differential expressed genes are listed in Additional file 11: Table S1. RNA samples used for RNA-Seq analysis were used for gene expression studies. The concentration of each RNA sample was quantified with NanoDrop 8000 (Thermo Fisher Scientific CN, Shanghai, China) and equal amount of RNA was used for cDNA synthesis. A One-Step gDNA Removal and cDNA Synthesis SuperMix kit (AU311–02; Trans Biotech, Beijing, China) was used for cDNA synthesis. The PCR products of the 10 primer pairs were sent for performing Sanger sequencing to confirm amplification of the expected amplicons before the conduct of RT-qPCR assay. A LightCycler 96 (Roche) platform and FastStart Essential DNA Green Master mix (Roche Molecular Systems, Inc., China) were used for all RT-qPCR assays. Three technical replicates were performed for each sample. The standard curves for each primer pairs were examined to determine the amplification efficiency of a certain gene. Only primer pairs that showed the expected amplicons, PCR efficiency estimate value from 1.9 to 2.1 and unique melting curves (see below) were used for gene expression analysis. The cycling conditions were set as followed: 95 °C for 10 min and 40 cycles of 95 °C for 10 s and 65 °C for 60 s. Fluorescence data were collected during the 65 °C step and analyzed with the LightCycler 96 Manager. After the running of cycling, melting curves were also examined by an addtiontional program of 95 °C for 10 s, 65 °C for 60 s and 97 °C for 1 s. Sample quantification cycle (Cq) values were determined and were standardized relative to the expression of the two reference genes (ACTIN and UBQL, Additional file 11: Table S1). The 2–ΔΔCq method was used to calculate the relative gene expression based on the RT-qPCR data [51]. The expression level at DAE0 for each tested gene was set as 1. All data analysis was performed by LightCycler® 96 Software (Version 1.1.0.1320, Roche, Shanghai, China).

Statistical analysis

All the Statistical analysis was conducted by using the R software (https://www.r-project.org/). The function ‘aov’ was used for variance analysis and all multiple comparisons was performed under the significant level of P value equal to 0.05.

Abbreviations

ABA:

Abscisic acid

AR:

Adventitious root

BR:

Brassinosteroid

DAE:

Days after excision

DMTU:

Dimethylthiourea

ET:

Ethylene

FPKM:

Fragments per kilobase million

GA:

Gibberellin

JA:

Jasmonic acid

PAT:

Polar auxin transport

ROS:

Reactive oxygen species

SA:

Salicylic acid

References

  1. Steffens B, Rasmussen A. The physiology of adventitious roots. Plant Physiol. 2016;170(2):603–17.

    CAS  PubMed  Google Scholar 

  2. Bellini C, Pacurar DI, Perrone I. Adventitious roots and lateral roots: similarities and differences. Annu Rev Plant Biol. 2014;65:639–66.

    CAS  PubMed  Google Scholar 

  3. Druege U, Franken P, Lischewski S, Ahkami AH, Zerche S, Hause B, Hajirezaei MR. Transcriptomic analysis reveals ethylene as stimulator and auxin as regulator of adventitious root formation in petunia cuttings. Front Plant Sci. 2014;5:494.

    PubMed  PubMed Central  Google Scholar 

  4. Zhao X, Zheng H, Li S, Yang C, Jiang J, Liu G. The rooting of poplar cuttings: a review. New For. 2014;45(1):21–34.

    Google Scholar 

  5. Abu-Abied M, Szwerdszarf D, Mordehaev I, Levy A, Stelmakh OR, Belausov E, Yaniv Y, Uliel S, Katzenellenbogen M, Riov J, 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(5):787–99.

    CAS  PubMed  Google Scholar 

  6. Pizarro A, Diaz-Sala C. Cellular dynamics during maturation-related decline of adventitious root formation in forest tree species. Physiol Plant. 2018;165(1):73–80.

    PubMed  Google Scholar 

  7. Legue V, Rigal A, Bhalerao RP. Adventitious root formation in tree species: involvement of transcription factors. Physiol Plant. 2014;151(2):192–8.

    CAS  PubMed  Google Scholar 

  8. Pacurar DI, Perrone I, Bellini C. Auxin is a central player in the hormone cross-talks that control adventitious rooting. Physiol Plant. 2014;151(1):83–96.

    CAS  PubMed  Google Scholar 

  9. Garrido G, Ramon Guerrero J, Angel Cano E, Acosta M, Sanchez-Bravo J. Origin and basipetal transport of the IAA responsible for rooting of carnation cuttings. Physiol Plant. 2002;114(2):303–12.

    CAS  PubMed  Google Scholar 

  10. Ahkami AH, Melzer M, Ghaffari MR, Pollmann S, Ghorbani Javid M, Shahinnia F, Hajirezaei MR, Druege U. Distribution of indole-3-acetic acid in Petunia hybrida shoot tip cuttings and relationship between auxin transport, carbohydrate metabolism and adventitious root formation. Planta. 2013;238(3):499–517.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. De Klerk GJ, Van der Krieken W, De Jong JC. Review - the formation of adventitious roots: new concepts, new possibilities. In Vitro Cell Dev-Pl. 1999;35(3):189–99.

    Google Scholar 

  12. Agullo-Anton MA, Ferrandez-Ayela A, Fernandez-Garcia N, Nicolas C, Albacete A, Perez-Alfocea F, Sanchez-Bravo J, Perez-Perez JM, Acosta M. Early steps of adventitious rooting: morphology, hormonal profiling and carbohydrate turnover in carnation stem cuttings. Physiol Plant. 2014;150(3):446–62.

    CAS  PubMed  Google Scholar 

  13. Steffens B, Wang JX, Sauter M. Interactions between ethylene, gibberellin and abscisic acid regulate emergence and growth rate of adventitious roots in Deepwater rice. Planta. 2006;223(3):604–12.

    CAS  PubMed  Google Scholar 

  14. Mauriat M, Petterle A, Bellini C, Moritz T. Gibberellins inhibit adventitious rooting in hybrid aspen and Arabidopsis by affecting auxin transport. Plant J. 2014;78(3):372–84.

    CAS  PubMed  Google Scholar 

  15. Lischweski S, Muchow A, Guthorl D, Hause B. Jasmonates act positively in adventitious root formation in petunia cuttings. BMC Plant Biol. 2015;15:229.

    PubMed  PubMed Central  Google Scholar 

  16. van Buer J, Cvetkovic J, Baier M. Cold regulation of plastid ascorbate peroxidases serves as a priming hub controlling ROS signaling in Arabidopsis thaliana. BMC Plant Biol. 2016;16(1):163.

    PubMed  PubMed Central  Google Scholar 

  17. Chen Z, Gu Q, Yu X, Huang L, Xu S, Wang R, Shen W, Shen W. Hydrogen peroxide acts downstream of melatonin to induce lateral root formation. Ann Bot. 2018;121(6):1127–36.

    PubMed  PubMed Central  Google Scholar 

  18. O'Brien JA, Daudi A, Finch P, Butt VS, Whitelegge JP, Souda P, Ausubel FM, Bolwell GP. A peroxidase-dependent apoplastic oxidative burst in cultured Arabidopsis cells functions in MAMP-elicited defense. Plant Physiol. 2012;158(4):2013–27.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Bindschedler LV, Dewdney J, Blee KA, Stone JM, Asai T, Plotnikov J, Denoux C, Hayes T, Gerrish C, Davies DR, et al. Peroxidase-dependent apoplastic oxidative burst in Arabidopsis required for pathogen resistance. Plant J. 2006;47(6):851–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Gao XQ, Liu CZ, Li DD, Zhao TT, Li F, Jia XN, Zhao XY, Zhang XS. The Arabidopsis KINbetagamma subunit of the SnRK1 complex regulates pollen hydration on the stigma by mediating the level of reactive oxygen species in pollen. PLoS Genet. 2016;12(7):e1006228.

    PubMed  PubMed Central  Google Scholar 

  21. Li SW, Leng Y, Shi RF. Transcriptomic profiling provides molecular insights into hydrogen peroxide-induced adventitious rooting in mung bean seedlings. BMC Genomics. 2017;18(1):188.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Velada I, Grzebelus D, Lousa D, Soares CM, Santos Macedo E, Peixe A, Arnholdt-Schmitt B, Cardoso HG. AOX1-Subfamily Gene Members in Olea europaea cv. "Galega Vulgar"-Gene Characterization and Expression of Transcripts during IBA-Induced in Vitro Adventitious Rooting. Int J Mol Sci. 2018;19(2):597.

    PubMed Central  Google Scholar 

  23. Liao WB, Huang GB, Yu JH, Zhang ML. Nitric oxide and hydrogen peroxide alleviate drought stress in marigold explants and promote its adventitious root development. Plant Physiol Biochem. 2012;58:6–15.

    CAS  PubMed  Google Scholar 

  24. Liao WB, Huang GB, Yu JH, Zhang ML, Shi XL. Nitric oxide and hydrogen peroxide are involved in indole-3-butyric acid-induced adventitious root development in marigold. J Hortic Sci Biotechnol. 2011;86(2):159–65.

    CAS  Google Scholar 

  25. Zalesny RS, Zalesny JA: Selecting Populus with different adventitious root types for environmental benefits, fiber, and energy. In: Niemi K, Scagel C (eds) Adventitious root formation in forest trees and woody horticultural crops—from genes to applications 2009, Chap. 18(Res Signpost):359–384.

  26. Li K, Liang Y, Xing L, Mao J, Liu Z, Dong F, Meng Y, Han M, Zhao C, Bao L, et al. Transcriptome Analysis Reveals Multiple Hormones, Wounding and Sugar Signaling Pathways Mediate Adventitious Root Formation in Apple Rootstock. Int J Mol Sci. 2018:19(8).

  27. Ramirez-Carvajal GA, Morse AM, Dervinis C, Davis JM. The cytokinin type-B response regulator PtRR13 is a negative regulator of adventitious root development in Populus. Plant Physiol. 2009;150(2):759–71.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Rigal A, Yordanov YS, Perrone I, Karlberg A, Tisserant E, Bellini C, Busov VB, Martin F, Kohler A, Bhalerao R, et al. The AINTEGUMENTA LIKE1 homeotic transcription factor PtAIL1 controls the formation of adventitious root primordia in poplar. Plant Physiol. 2012;160(4):1996–2006.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Xu M, Xie W, Huang M. Two WUSCHEL-related HOMEOBOX genes, PeWOX11a and PeWOX11b, are involved in adventitious root formation of poplar. Physiol Plant. 2015;155(4):446–56.

    CAS  PubMed  Google Scholar 

  30. Li J, Zhang J, Jia H, Liu B, Sun P, Hu J, Wang L, Lu M. The WUSCHEL-related homeobox 5a (PtoWOX5a) is involved in adventitious root development in poplar. Tree Physiol. 2018;38(1):139–53.

    PubMed  Google Scholar 

  31. Ghosh S, Chan CK. Analysis of RNA-Seq data using TopHat and cufflinks. Methods Mol Biol. 2016;1374:339–61.

    CAS  PubMed  Google Scholar 

  32. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7(3):562–78.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Ahkami AH, Lischewski S, Haensch KT, Porfirova S, Hofmann J, Rolletschek H, Melzer M, Franken P, Hause B, Druege U, et al. Molecular physiology of adventitious root formation in Petunia hybrida cuttings: involvement of wound response and primary metabolism. New Phytol. 2009;181(3):613–25.

    CAS  PubMed  Google Scholar 

  34. Zhou Z, Tian Y, Cong P, Zhu Y. Functional characterization of an apple (Malus x domestica) LysM domain receptor encoding gene for its role in defense response. Plant science : An Int J Exp Plant Biol. 2018;269:56–65.

    CAS  Google Scholar 

  35. Abu-Abied M, Rogovoy Stelmakh O, Mordehaev I, Grumberg M, Elbaum R, Wasteneys GO, Sadot E. Dissecting the contribution of microtubule behaviour in adventitious root induction. J Exp Bot. 2015;66(9):2813–24.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Abu-Abied M, Szwerdszarf D, Mordehaev I, Yaniv Y, Levinkron S, Rubinstein M, Riov J, Ophir R, Sadot E. Gene expression profiling in juvenile and mature cuttings of Eucalyptus grandis reveals the importance of microtubule remodeling during adventitious root formation. BMC Genomics. 2014;15:826.

    PubMed  PubMed Central  Google Scholar 

  37. Smith NG, Wareing PF. The distribution of latent root primordia in stems of Populus×robusta, and factors affecting the emergence of preformed roots from cuttings. Forestry. 1972;45(2):197–209.

    Google Scholar 

  38. Carlson MC. The formation of nodal adventitious roots in Salix Cordata. Am J Bot. 1938;25(9):721–5.

    Google Scholar 

  39. Carlson MC. Nodal adventitious roots in willow stems of different ages. Am J Bot. 1950;37(7):555–61.

    Google Scholar 

  40. Finkel T, Holbrook NJ. Oxidants, oxidative stress and the biology of ageing. Nature. 2000;408(6809):239–47.

    CAS  PubMed  Google Scholar 

  41. Singh R, Singh S, Parihar P, Mishra RK, Tripathi DK, Singh VP, Chauhan DK, Prasad SM. Reactive oxygen species (ROS): beneficial companions of Plants' developmental processes. Front Plant Sci. 2016;7:1299.

    PubMed  PubMed Central  Google Scholar 

  42. Mangano S, Denita-Juarez SP, Choi HS, Marzol E, Hwang Y, Ranocha P, Velasquez SM, Borassi C, Barberini ML, Aptekmann AA, et al. Molecular link between auxin and ROS-mediated polar growth. P Natl Acad Sci USA. 2017;114(20):5289–94.

    CAS  Google Scholar 

  43. Mccown BH, Lloyd G: Woody Plant medium (wpm) - a mineral nutrient formulation for microculture of Woody Plant-species. Hortscience 1981, 16(3):453–453.

  44. You J, Fang Y, Xiong L: Reactive oxygen detection. Bio-101 2018, e1010170.

  45. Masters TL, Hilker CA, Jeraldo PR, Bhagwate AV, Greenwood-Quaintance KE, Eckloff BW, Chia N, Hanssen AD, Abdel MP, Yao JZ, et al. Comparative evaluation of cDNA library construction approaches for RNA-Seq analysis from low RNA-content human specimens. J Microbiol Methods. 2018;154:55–62.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & gray). Science. 2006;313(5793):1596–604.

    CAS  PubMed  Google Scholar 

  48. Jiang Z, Liu X, Peng Z, Wan Y, Ji Y, He W, Wan W, Luo J, Guo H. AHD2.0: an update version of Arabidopsis hormone database for plant systematic studies. Nucleic Acids Res. 2011;39(Database issue):D1123–9.

    CAS  PubMed  Google Scholar 

  49. Steinley D, Brusco MJ. Choosing the number of clusters in kappa-means clustering. Psychol Methods. 2011;16(3):285–97.

    PubMed  Google Scholar 

  50. Yi X, Du Z, Su Z. PlantGSEA: a gene set enrichment analysis toolkit for plant community. Nucleic Acids Res. 2013;41(Web Server issue):W98–103.

    PubMed  PubMed Central  Google Scholar 

  51. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(−Delta Delta C) method. Methods. 2001;25(4):402–8.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors are grateful to American Journal Experts (AJE) for improving the English in this paper.

Funding

This work was supported by the National Natural Science Foundation of China (NSFC accession No. 31670651). The funding agency was not involved in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

All of the raw data have been deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/sra/) under project number PRJNA516245.

Author information

Authors and Affiliations

Authors

Contributions

NW organized the entire project, and NW, YZ, ZX, CZ, ML and WX performed the experiments and data analysis. NW wrote and edited this manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nian Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S2. Genes differentially expressed during poplar AR formation. The numbers indicate FPKM values. (XLSX 507 kb)

Additional file 2:

Figure S1. Dendrogram of RNA-Seq samples used for time course transcriptome analysis. (PPTX 125 kb)

Additional file 3:

Figure S2. Number of genes that are differentially expressed during poplar AR formation. (PPTX 9022 kb)

Additional file 4:

Table S3. Phytohormone-related genes in poplar (XLSX 30 kb)

Additional file 5:

Table S4. Normalized expression values of differentially expressed phytohormone-related genes (XLSX 28 kb)

Additional file 6:

Figure S3. Boxplot of FPKM values for differentially expressed genes related to phytohormone transport. FPKM represents fragments per kilobase of transcript per million fragments mapped. Each boxplot consists of the Q1 (1/4 quartiles, the bottom of the box), Q3 (3/4 quartiles, the top bottom of the box), the median Q2 (the line inside into the box), Q3 + 1.5 IQR (interquartile range, Q3-Q1, the highest line outside the box) and Q1–1.5 IQR (the lowest line outside the box). Values larger than Q3 + 1.5 IQR or less than Q1–1.5 IQR were considered as outliers and they were plotted as dots outside of the Q3 + 1.5 IQR or Q1–1.5 IQR. If Q3 + 1.5 IQR is less than Q3 or Q1–1.5 IQR is less than Q1, their corresponding lines will not be shown. (PPTX 38 kb)

Additional file 7:

Figure S4. Boxplot of FPKM values for differentially expressed genes related to phytohormone metabolism. The figure legends are similar to Additional file 3: Figure S3. (PPTX 44 kb)

Additional file 8:

Figure S5. Boxplot of FPKM values for differentially expressed genes related to phytohormone receptors. The figure legends are similar to Additional file 7: Figure S3. (PPTX 45 kb)

Additional file 9:

Table S5. Functional enrichment of genes in the 12 clusters. (XLSX 15 kb)

Additional file 10:

Figure S6. Hydroponic culture of poplar cuttings with different concentrations of supplemented H2O2 at DAE11. Panels a, b, and c represent poplar cuttings grown in 0, 500, and 1000 μM H2O2, respectively. Bar = 1 cm. (PPTX 579 kb)

Additional file 11:

Table S1. Primers used for RT-qPCR assay. (XLSX 9 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Xiao, Z., Zhan, C. et al. Comprehensive analysis of dynamic gene expression and investigation of the roles of hydrogen peroxide during adventitious rooting in poplar. BMC Plant Biol 19, 99 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-1700-7

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords