Abstract
Poplar is a major afforestation tree species in flood-prone areas. Here, we compared molecular and physiological responses in the roots of two full-sib poplar clones, LS1 (flood-tolerant) and LS2 (flood-susceptive), subjected to stagnant flooding using transcript and metabolite profiling. LS1 displayed less phenotypic damage and superior leaf gas exchange and plant growth compared with those of LS2. We concluded that three characteristics might contribute to the differences in flood tolerance between LS1 and LS2. First, fermentation was initiated through lactic dehydrogenation in LS1 roots under flooding and subsequently dominated by alcohol fermentation. However, lactic dehydrogenase was persistently active in flooded LS2. Second, 13 differentially expressed genes associated with energy and O2 consumption processes under soil flooding had lower transcript levels in LS1 than those in LS2, which might contribute to better energy-/O2-saving abilities and behaviours in flood-tolerant LS1 than those in flood-susceptible LS2 under hypoxic stress. Third, LS1 possessed increased reactive oxygen species scavenging abilities compared with those of LS2 under edaphic flooding. Our data are a valuable contribution to understanding the mechanisms involved in the flood tolerance of poplar.
Similar content being viewed by others
Introduction
The growth of plants is usually negatively affected by flooding stress. Because of the characteristics of fast growth and flood tolerance, poplar has served as a major afforestation tree species in flood-prone areas to utilize these flood-affected soils1. However, flooding stress still negatively affects poplar growth and survival in these regions2,3. Thus, understanding the mechanism involved in the flood tolerance of poplar is important, which may facilitate developing poplar clones with improved flood tolerance.
Plants evolved various mechanisms to withstand flooding stress, including features of morphology and anatomy, metabolism and molecular transcriptional regulation4,5,6,7,8,9. The primary metabolic responses to flooding in plants include the stimulation of fermentative pathways and an increase in glycolytic flux, indicated by elevated transcript abundance, increased activities of pyruvate decarboxylase and alcohol dehydrogenase (ADH), and highly accumulated product (acetaldehyde and ethanol) contents in flooded roots10,11. Because of the energy limitation that results from the switch from respiration to fermentation, ammonium and nitrate uptake in roots is also often strongly impaired by flooding12. Transcript abundance profiles demonstrate that many genes involved in amino acid biosynthesis and degradation are differentially expressed in flooded trees13,14,15. Many of the amino acids that increase under hypoxia are formed from pyruvate (e.g., leucine, valine and alanine) or other intermediates of glycolysis (e.g., tyrosine, serine and glycine). By contrast, contents of amino acids derived from tricarboxylic (TCA) cycle components (e.g., glutamine, asparagine, glutamate and aspartate) are often lower in hypoxic than in normoxic roots10,12,13. Reactive oxygen species (ROS) generating and scavenging systems also play important roles in the ability of plants to respond to various abiotic and biotic stresses16. In addition to the role in oxidative damage, ROS also act as signal molecules for adaptive stress responses17. In the adaptation of plants to flooding, the group VII Ethylene Response Factor genes (ERF-VIIs) play pivotal roles through regulation of anaerobic gene expression and antithetical survival strategies18.
Characteristics of morphology, anatomy, eco-physiology and growth of poplar under hypoxia stress are well characterized in the previous literature, accompanied by a few studies on the transcriptome and metabolome19,20,21. Kreuzwieser (2009)11 revealed metabolite changes occur in both leaves and roots of flooded grey poplar (Populus × canescens), but changes in transcript abundance were restricted to the roots. Additionally, increased glycolysis, induced fermentation and inhibited energy-intensive processes are crucial for grey poplar (relatively flood-tolerant) to grow better than Arabidopsis (relatively flood-sensitive) under low oxygen conditions. Compared with cotton and Arabidopsis, a waterlogging-tolerant poplar clone displayed a different transcriptome under early stages of hypoxia stress14. The tolerance to flooding differs not only between species but also between different poplar clones22. To date, the mechanism that leads to those differences between poplar clones has been ambiguous. Peng (2013)22 clarified that roots play a vital role in flood-tolerance of poplar demonstrated by reciprocal grafting. According to previous reports, the root is the most important organ in response to flooding stress in poplar13,20,23. Therefore, molecular and physiological responses in roots of two poplar clones following root hypoxia were compared using transcript and metabolite profiling in this study. The two plant materials (flood-tolerant LS1 and flood-susceptible LS2) were full-sib clones and differed in flood tolerance. Under 15-day partial submergence treatment followed by 3-day recovery, the total biomass of LS1 and LS2 decreased by 6.38% and 32.01%, respectively, compared with their watered control20. After partial submergence for 21 days and recovery for 6 days, the reduction in the rate of total biomass for LS1 and LS2 was 24.4% and 59.9%, respectively23. We aimed to detect the different mechanisms of the two clones in response to flooding stress and to identify candidate genes involved in flood tolerance that may be used in molecular breeding programmes to improve flood tolerance of poplar.
Results
Phenotype, leaf gas exchange, chlorophyll fluorescence and plant growth
Soil flooding negatively influenced the phenotype and most parameters of photosynthesis, chlorophyll fluorescence and growth in both LS1 and LS2 (Fig. 1, Table 1). By the end of the experiment, all plants survived. All watered seedlings grew vigorously without obvious phenotype damage, but for the flooded plants, different levels of phenotype damage occurred between LS1 and LS2 (Fig. 1a). The flooded LS1 grew more new leaves and showed less leaf defoliation than the LS2 during the 18-day treatment (Fig. 1a,b). Height and root-collar diameter growth in the two clones were both dramatically restrained by soil flooding (p < 0.05, Fig. 1b). In LS2_FL, net photosynthesis rate (Pn, p < 0.05), potential efficiency of primary conversion of light energy of PSII (Fv/Fm, p < 0.01) and the ratio of variable fluorescence to initial fluorescence (Fv/Fo, p < 0.01) values were significantly lower than those of LS2_CK (Fig. 1c). However, in LS1_FL, non-significant changes were observed with the exception of a decrease in Pn on day 7 (p < 0.01, Fig. 1c). Stomatal conductance (Gs) was not noticeably affected by flooding stress in either LS1 or LS2 (Fig. 1c). A reduction in rates due to flooding stress in parameters of leaf gas exchange, chlorophyll fluorescence and plant growth was markedly different (p < 0.05) between LS1 and LS2 at every time point, with exceptions of Gs and 7-day Pn. Therefore, LS1 performed better and exhibited more growth than LS2 under soil flooding stress.
Principal component analysis and gene expression
The major characteristics of all libraries are summarized in Supplementary Table 2. Waterlogging stress caused changes in global gene expression levels in both LS1 and LS2 roots (Supplementary Fig. 1). Compared with their controls, a total of 1,161 and 959 genes were only expressed in flooded LS1 and LS2 roots, respectively (Fig. 2a). LS1_FL1 had 74 more differentially expressed genes (DEGs) than LS2_FL1, and LS1_FL7 had 86 fewer DEGs than LS2_FL7 (Fig. 2b). Several clone-specific expression genes (Fig. 2a) and DEGs between LS1 and LS2 were also investigated (Fig. 2b). Compared with LS2, DEG numbers in LS1 on days zero, 1 and 7 were 904, 773 and 1,375, respectively.
Quantitative real-time polymerase chain reaction validation of the DEGs from the RNA-seq analysis
The correlations between results of RNA-seq and quantitative real-time polymerase chain reaction (qRT-PCR) were evaluated by the expression of the selected 18 DEGs in Supplementary Table 1, including ten transcription factor genes and four genes encoding some important enzymes involved in plant response to abiotic stress. Scatterplots were generated by comparing the log2-fold changes determined by RNA-seq and qRT-PCR in 1-day flooding versus CK and 7-day flooding versus CK in each clone. These transcripts displayed similar expression patterns between qRT-PCR and RNA-seq experiments illustrated by 0.895 and 0.925 correlation coefficients between the two sets of data for LS1 and LS2, respectively (Supplementary Fig. 2). These results confirmed the accuracy and reproducibility of the RNA-seq profiles in this study.
Functional annotation and classification of DEGs
The gene ontology (GO) terms of the unigenes were categorized into molecular functions, cellular components and biological processes. A total of 12,947 unigenes were assigned in the secondary category in accordance with the GO assignments (Supplementary Table 4). Compared with the control, a total of 4,695 (accounting for 36.26% of the 12,947 unigenes), 5,740 (44.33%), 4,190 (32.36%) and 5,515 (42.60%) DEGs in LS1_FL1, LS1_FL7, LS2_FL1 and LS2_FL7, respectively, were assigned to GO terms (Supplementary Table 5). A total number of 49 GO terms were significantly enriched in either LS1 or LS2 under flooding treatment (corrected p < 0.05), including 30 GO terms in both (Fig. 3a). Most significantly enriched GO terms had more down-regulated genes than up-regulated genes (Supplementary Table 5). Fourteen GO terms were specifically enriched in LS1_FL7. One and four GO terms were specifically significantly enriched in LS2_FL1 and LS2_FL7, respectively. Cell wall biogenesis, an energy-consuming biological process, was only significantly enriched in LS1_FL7 with 6 up-regulated DEGs and 22 down-regulated DEGs (Supplementary Table 5).
DEGs were also classified into pathways of poplar metabolism based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Fig. 3b; Supplementary Table 6). In total, 9,906 genes were assigned into KEGG pathways; 4,551 (accounting for 45.94% of the 9,906 genes), 5,575 (56.28%), 4,084 (41.23%) and 5,366 (54.17%) DEGs were, respectively, classified into pathways in LS1_FL1, LS1_FL7, LS2_FL1 and LS2_FL7, and the biosynthesis of secondary metabolites was the most enriched pathway. Most enriched pathways in the two clones were classified into carbon metabolism, nitrogen metabolism, lipid metabolism and energy metabolism, including glycolysis, fermentation, fatty acid metabolism and amino acid metabolisms. The top 10 DEGs and corresponding KEGG pathways are listed in Supplementary Table 7.
Effects on carbon metabolism
Soil flooding caused increased glycolysis, inhibited O2-demanding processes (e.g., TCA cycle), and activated alcohol and lactic fermentation in the present study (Supplementary Table 6). Most genes associated with glycolytic processes were up-regulated by hypoxia stress. For example, phosphohexokinase (EC 2.7.1.11) gene POPTR_0006s25170 and pyruvate kinase (EC 2.7.1.40) gene POPTR_0008s00310 were up-regulated more than 48-fold and 10-fold, respectively, in LS1_FL and LS2_FL, accompanied with dramatically increased soluble sugar (sucrose, fructose and glucose) contents on day 7 (Fig. 4a). Expression levels of almost half the genes of the TCA cycle were decreased by 1-day flooding and further decreased on day 7. Pyruvate decarboxylase (EC 4.1.1.1) gene POPTR_0011s07630 and ADH (EC 1.1.1.1) gene POPTR_0002s07290, two key genes associated with ethanol fermentation, were activated (>9.5-fold and 19.6-fold, respectively) under flooding. In parallel with the transcript levels, ADH activities increased significantly on the 1st day in both clones (p < 0.05) and were highest on the 7th day (p < 0.05, Fig. 5a). For lactate dehydrogenase (LDH), activities increased significantly in both LS1_FL and LS2_FL (p < 0.05), with the exception of LS1_FL7 (Fig. 5a).
In this study, the sucrose degradation patterns were strongly affected in the roots of LS1_FL and LS2_FL induced by sucrose synthase (EC 2.4.1.13) and invertase (EC 3.2.1.26). For sucrose synthase, POPTR_0002s20340 (SUSY3) was up-regulated more than 21.6-fold on the 1st day and 6-fold on the 7th day in the two clones (Supplementary Table 6). For invertase, one gene POPTR_0003s11210 was downregulated in both LS1_FL (>27-fold) and LS2_FL (>6.4-fold). Another invertase-encoding gene POPTR_0006s22710 was up-regulated (>6.6-fold) in LS2_FL but was unchanged in LS1_FL (Supplementary Table 6). Some sucrose transporter genes (SUT, SWEET and TMT) were also differentially expressed between different treatments or clones (Supplementary Table 8).
Effects on nitrogen metabolism
In comparison with the control, NO3− uptake in both LS1_FL and LS2_FL roots was most likely negatively affected as indicated by the 5–6 significantly down-regulated NO3− transporter genes and decreased soluble protein contents (Supplementary Table 8; Fig. 4b). Two NO3− transporter genes, POPTR_0015s09290 (PTNRT2–6, >5.5-fold) and POPTR_0015s09310 (PTRNTR2–3, >23.6-fold), were upregulated on day 7. Different response levels also occurred between LS1 and LS2. Compared with LS2, LS1 showed four significantly down-regulated transporter genes on day zero and three significantly upregulated transporter genes on day 7 (Supplementary Table 8). For NH4+ uptake, numbers of transporter genes that showed increased transcript abundance were markedly greater than those of down-regulated transporter genes in both flooded LS1 and LS2 (Supplementary Table 8). Compared with LS2, LS1 showed six significantly up-regulated transporter genes on day 7.
A wide variety of amino acid contents in LS1 and LS2 roots changed markedly during flooding treatment, together with the transcript levels of corresponding genes associated with amino acid metabolism (Supplementary Table 6; Fig. 4c). Most amino acids closely derived from pyruvate (e.g., Ala, Val) and intermediates of glycolysis (e.g., Gly, Tyr) were highly accumulated, whereas contents of most amino acids derived from intermediates of the TCA cycle (e.g., Glu, Asp and Asn) decreased, with the exception of increased Tyr and Val in LS2_FL7 and Asn in LS2_FL1. These results supported the responses of carbon metabolism (i.e., glycolysis and TCA cycle) to flooding in LS1 and LS2 roots at the transcript level.
Different responses to flooding between LS1 and LS2 roots were also observed in amino acid contents (Fig. 4c). For example, in comparison with the controls, Arg content increased non-significantly in both LS1_FL1 and LS1_FL7 but increased significantly in LS2_FL1 (p < 0.05) and was unchanged in LS2_FL7. Glu content decreased significantly in both LS1_FL1 and LS1_FL7 (p < 0.01), whereas the content significantly decreased in LS2_FL1 (p < 0.01) and increased in LS2_FL7 (p < 0.01). For Phe content, a significant increase was observed in LS1_FL (p < 0.05) but not in LS2_FL.
Different transcript levels of genes involved in amino acid metabolism were also detected between LS1_FL and LS2_FL (Supplementary Table 6). In Arg metabolism, expression of the nitric-oxide synthase (NOS, EC 1.14.13.39) gene POPTR_0001s00700 was up-regulated in LS2_FL1 (2.4-fold) but uninfluenced in LS1_FL1. In Phe catabolism, transcription of the 4-hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27) gene POPTR_0002s05840 was up-regulated in LS2_FL roots (5.2-fold and 2.1-fold on the 1st and 7th day, respectively) but was uninfluenced in LS1_FL roots, which explained the response of the Phe contents in the two clones (Fig. 4c).
Effects on reactive oxygen species generation and scavenging
Soil flooding affected the systems of ROS generation and scavenging in both LS1 and LS2, and the two clones responded differently (Supplementary Table 6; Fig. 5b). For antioxidant enzymes, superoxide dismutase (SOD, EC 1.15.1.1) activities increased in both LS1_FL and LS2_FL, but the increase was significantly greater in LS1_FL (p < 0.01). Correspondingly, the SOD gene POPTR_0013s05350 was also up-regulated in both LS1_FL (>5.9-fold) and LS2_FL (4.2-fold). Peroxidase (POD) activities were markedly elevated in LS1_FL (p < 0.05) and LS2_FL1 (p < 0.05) but not in LS2_FL7. At the gene transcription level, the POD (EC 1.11.1.7) gene POPTR_0012s00700 was up-regulated more than 58.4-fold and 36.4-fold in LS1_FL and LS2_FL, respectively. However, another POD gene POPTR_0004s01510 was downregulated more than 21.6-fold and 30.0-fold in LS1_FL and LS2_FL, respectively. LS2_FL7 displayed higher CAT (EC 1.11.1.6) activity than that of LS1_FL7 (p < 0.01). The catalase (CAT) gene POPTR_0002s01080 was also up-regulated more in LS2 than that in LS1 with a fold change of 2.6, 1.5, 1.8 and 1.1 in LS2_FL1, LS2_FL7, LS1_FL1 and LS1_FL7, respectively. Glutathione reductase (GR, EC 1.8.1.7) activities were, respectively, 109.8% and 144.0% of the control in LS1_FL1 and LS1_FL7 and 130.7% and 67.6% of the control in LS2_FL1 and LS2_FL7. Correspondingly, the GR gene POPTR_0015s04650 transcript level was up-regulated 6.9-fold, 4.4-fold, 6.5-fold and 3.7-fold. The glutathione peroxidase (GPX, EC 1.11.1.9) gene POPTR_0003s12620 was up-regulated in both LS1_FL (>1.6-fold) and LS2_FL (>2.0-fold). For non-enzymatic antioxidants, glutathione (GSH) contents were, respectively, 132.1% and 137.2% of the control in LS1_FL1 and LS1_FL7 and 163.1% and 110.1% of the control in LS2_FL1 and LS2_FL7. Ascorbate (Asc) dramatically accumulated in LS1_FL7 (p < 0.01, 258.9% of the control) and LS2_FL7 (p < 0.01, 302.4% of the control) but not in LS1_FL1 and LS2_FL1. For H2O2, the content was only significantly elevated in LS2_FL1 (p < 0.01, 145.2% of the control) but not in LS1_FL and LS2_FL7.
DEGs associated with ATP/O2 metabolism and transcription factors
Thirteen DEGs associated with ATP/O2 metabolism differentially responded to flooding stress between LS1 and LS2 (Fig. 6a), and ATP/O2 consumption genes were less expressed in LS1_FL, which might contribute to the difference in flood tolerance. The ATP phosphoribosyltransferase (EC 2.4.2.17) gene POPTR_0019s08520, involved in ATP generation, was only up-regulated (2-fold) in LS1_FL7. The other 12 genes involved in O2/ATP consumption all had lower transcript levels in LS1 roots than those in LS2 roots under flooding stress, including the invertase gene POPTR_0006s22710, 3′-phosphoadenosine 5′-phosphosulfate synthase (EC 2.7.7.4/2.7.1.25) gene POPTR_0008s15880, acetyl-CoA acyltransferase 1 gene (EC 2.3.1.16) POPTR_0001s14280, acyl-CoA oxidase (ACX) gene POPTR_0016s12540, NOS gene POPTR_0001s00700, carbamoyl-phosphate synthase gene POPTR_0010s02380 and 4-hydroxyphenylpyruvate dioxygenase gene POPTR_0002s05840, among others.
A total number of 1,038 and 1,006 transcription factors (belonging to 52 families) were differentially expressed in flooded LS1 and LS2 compared with their controls, respectively. AP2/ERF, bHLH, MYB, NAC and bZIP were the top 5 families that contained the most DEGs in both LS1 and LS2. Expression patterns of several transcription factors differentially responded to flooding stress between LS1 and LS2, and 20 (from 17 families) are listed in Fig. 6b, including AP2/ERF, MYB, WRKY, NAC, and HSF, among others. The 20 genes were differentially expressed between LS1 and LS2 on both day 1 and day 7 but not on day zero.
Discussion
In comparison with LS1, LS2 showed more phenotypic damage and greater reductions in values of Pn, Fv/Fm, Fv/Fo, height growth and root-collar diameter growth under flooding stress (Fig. 1, Table 1), which indicated a more severe flood-injury and an inferior flood tolerance23,24. Generally, restrictions on both stomatal (stomatal closure, indicated by Gs) and non-stomatal factors (factors except for stomatal closure, i.e., biochemical reactions of photosynthesis) can result in decreased photosynthesis efficiency under abiotic stress25. In our study, Pn values decreased significantly in both flooded LS1 and LS2, whereas their Gs values were reduced non-significantly. The decrease in Pn values might be primarily caused by non-stomatal factors, such as the reduced activity of ribulose-1,5-bisphosphate carboxylase (EC 4.1.1.39), glycolate oxidase (EC 1.1.3.15) and phosphoglycolate phosphatase (EC 3.1.3.18), destruction of chloroplast membrane structure, and inhibition of photosynthetic electron transport and PSII activity26,27. A similar result was also reported in long-term flooded (up to 5 d) maize plants (Zea mays L.), and soil flooding led to a decrease in photosynthesis and ribulose-1,5-bisphosphate carboxylase activity without significant changes in the values of stomatal conductance26.
A global analysis of the transcriptome facilitates the identification of systemic gene expression and regulatory mechanisms for tolerance to a certain stress in a plant28,29,30. The DEGs identified in our study between watered and flooded plants were mostly enriched in glycolysis, fermentation, fatty acid metabolism and amino acid metabolism, consistent with the performance of grey poplar under hypoxia13.
The increased glycolysis, stimulated alcohol and lactic fermentation, and inhibited TCA cycle pathways in this study are consistent with the transcript profile responses of rice, Arabidopsis, and poplar under hypoxia31. Soil flooding stress led to limited carbon entry into the TCA cycle and directed carbon entry into fermentation13, which was also indicated by the transcripts of genes involved in these biological processes in this study (Supplementary Table 6). Because roots usually suffer an energy crisis under hypoxia due to insufficient ATP produced by the TCA cycle, to compensate, glycolysis is heightened, and aerobic respiration transforms into fermentation to maintain ATP production13. Therefore, higher transcript levels of pyruvate decarboxylase, pyruvate kinase and ADH (Supplementary Table 6), coupled with stronger ADH activity (Fig. 5a), suggested a superior ability of ATP production in LS1_FL than that in LS2_FL, in parallel with their difference in flood tolerance. Different fermentation patterns also occurred in hypoxic LS1 and LS2 roots (Fig. 5a). In LS1_FL, fermentation was initiated by activating lactate fermentation and then quickly dominated by alcohol fermentation, elucidated by the activities of LDH and ADH (Fig. 5a); this process is similar to that in reports with grey poplar13, Arabidopsis32 and maize33 and consistent with the pH-stat hypothesis34. In LS2_FL, however, lactate fermentation played a sustained role during flooding treatment (Fig. 5a). A more efficient fermentation pathway was previously reported to partly explain the improved waterlogging tolerance of flood-tolerant oak15, and excessive lactic dehydrogenation could cause continual cytoplasmic acidification, inhibited glycolysis and irreversible cell injury34.
Expression of genes associated with nitrogen metabolism in roots was negatively influenced in both LS1_FL and LS2_FL. The decreased expression levels of NO3− transporter genes and significantly reduced soluble protein contents (Fig. 4b) are consistent with previous observations of reduced nitrogen uptake in flooded poplar12,20. The internal transport of NO3− might be more active, whereas NO3− uptake was reduced in roots, as suggested by the upregulated transcript levels of PTNRT2-6. PTNRT2-6 performs a role in the transport of NO3− from stored pools, such as vacuoles, to the cytoplasm35. PTRNTR2-3 might also have a similar function to that of PTNRT2-6, as suggested by the similar response in expression to root hypoxia. The upregulated NH4+ transporter genes could be involved in root internal transport, such as NH4+ loading into xylem14.
In our study, the contents of a wide variety of amino acids responded rapidly to flooding stress in both LS1 and LS2. Consistent with the results in hypoxic grey poplar13 and Arabidopsis36, the most intensively accumulated amino acids were closely derived from pyruvate (e.g., Ala, Val) and glycolysis intermediates (e.g., Gly, Tyr), whereas the amino acids with the greatest decrease in contents were derived from intermediates of the TCA cycle (e.g., Glu, Asp and Asn). Different responses to flooding also occurred between LS1 and LS2 in amino acid contents and the corresponding transcripts (Supplementary Table 6; Fig. 4c). For example, Arg content and the transcript level of the NOS gene POPTR_0001s00700 both increased significantly (p < 0.05) in LS2_FL1 but not in LS1_FL1 (Fig. 4c, Supplementary Table 6). NOS participates in Arg consumption and nitric oxide (NO) formation37. Thus, the increase in POPTR_0001s00700 expression might be a feedback regulation for Arg accumulation in LS2_FL1. Additionally, increased NOS activity might also lead to increased NO emission in LS2_FL, as reported in several flood-sensitive tree species38,39. Compared with that of LS2_FL, Phe was significantly accumulated in LS1_FL (p < 0.01, Fig. 4c), and the lower Phe content in LS2_FL might result from the elevated transcript level of POPTR_0002s05840 in LS2_FL. This gene encodes 4-hydroxyphenylpyruvate dioxygenase, an important enzyme in catalysing Phe catabolism in most plant organs40. In comparison with LS1_FL, upregulated POPTR_0002s05840 expression in LS2_FL would expedite Phe consumption, consistent with the Phe contents observed.
LS1 had stronger ROS-scavenging abilities than those of LS2 under soil flooding in this study, which could contribute to its superior flood tolerance. A wide variety of enzymes and non-enzymatic antioxidants participate in the scavenging of ROS. In enzymatic systems, SOD constitutes the first line of defence against ROS by dismutating O2− to H2O2, and then H2O2 is decomposed by POD and CAT41. Compared with those in LS2_FL, the higher activities of SOD and POD in LS1_FL could contribute to its higher oxidation resistance under hypoxia. As one of the ROS, H2O2 significantly accumulated in LS2_FL1 (p < 0.01, 145.2% of the control) but not in LS1_FL1. The excess H2O2 could cause severe damage to cell membrane systems and flood-injury in LS242. SOD, L-aspartate oxidase and ACX were pivotal enzymes involved in H2O2 production in this study43,44,45,46. The L-aspartate oxidase gene POPTR_0001s34980 and ACX gene POPTR_0016s12540 were more upregulated in LS2_FL1 than in LS1_FL1, which might increase the activities of the two enzymes and contribute to the high H2O2 content in LS2_FL1. POD, CAT, Asc peroxidase, Asc and GSH were key enzymatic/non-enzymatic antioxidants associated with H2O2 degradation in the two clones47,48,49. In plants, the primary substrate for reductive detoxification of H2O2 is Asc, and GSH primarily participates in the re-reduction of Asc rather than H2O2 degradation50. In the Asc-GSH cycle, GSH can be oxidized to GSH disulfide and be reduced from GSH disulfide by GR51. In our study, GR activity and GSH content both dramatically increased in LS2_FL1 (Fig. 5b). Previous research on the SOD-Asc-GSH cycle highlighted that when SOD activity is low, non-enzymatic antioxidants (e.g., GSH) always play an important role in O2− scavenging and induce additional H2O2 production51. Thus, considering the low SOD activity in LS2_FL1, high GSH content most likely contributed to the over-accumulation of H2O251. On day 7, H2O2 in flooded LS2 recovered to the control level under the effect of reduced GR activity and GSH content, increased SOD and CAT activities, and increased Asc content (Fig. 5b).
Plants commonly suffer from O2 deprivation under waterlogging, and the hypoxic condition can cause an energy crisis and induce down-regulated energy consumption10,52,53. We concluded that 13 genes differentially responded to partial submergence between LS1 and LS2, whose products primarily participated in energy/O2-related metabolism in amino acid metabolism, sucrose degradation, sulfur reduction, and beta-oxidation of fatty acids, among others. For example, POPTR_0006s22710 encodes invertase and participates in sucrose degradation. Soil flooding stress can suppress invertase expression and promote sucrose synthase expression, switching to the more energy saving pathway13. In our results, the transcript level of one invertase gene POPTR_0003s11210 showed greater decrease in LS1_FL than that in LS2_FL (p < 0.05), and another invertase gene POPTR_0006s22710 only increased in LS2_FL (p < 0.05) but not in LS1_FL (Supplementary Table 6), which would cause higher enzyme activity and more energy consumption in LS2_FL than in LS1_FL. In comparison with LS2_FL, lower transcript levels of two genes (POPTR_0016s12540 and POPTR_0001s14280) associated with beta-oxidation of fatty acids might cause a lower activity of beta-oxidation in LS1_FL roots. Beta-oxidation of fatty acids consists of a series of O2-requiring processes located in the peroxisomes of plants54. Thus, lower beta-oxidation activity in LS1 would also attenuate O2 consumption under hypoxia. Other genes listed in Fig. 6a also increased or were not inhibited in flooded LS2 roots, which might lead to over-consumption of energy and O2, resulting in severe energy starvation and O2 deficiency. By contrast, those genes involved in energy/O2 consumption were inhibited in flooded LS1 and contributed to maintaining available energy and O2 under hypoxic conditions.
Years of studies on biotic and abiotic stresses show that transcription factors play very important roles in regulating plant stress resistance. Members of families AP2/ERF, WRKY, NAC and HSF participate in regulating flood tolerance of plants27,55,56,57,58. In our study, 20 transcription factors were identified that were differentially expressed between flooded LS1 and LS2, including AP2/ERF, MYB, WRKY, NAC, HSF, GATA, GRF, bZIP, B3, M-type, Dof, C2H2, GRAS, MIKC, HD-ZIP, MYB_related and LBD family genes (Fig. 6b). Family AP2/ERF was well characterized, and abundant AP2/ERF family genes participate in regulating flood tolerance in rice, Arabidopsis and soybean, such as Sub1A, AtERF73/HRE1, SNORKEL, and OsEATB, among others4,5,59,60,61,62. Hypoxia induction of ADH expression levels was reduced in three AtERF73/HRE1 knockout Arabidopsis lines, whereas the expression of ADH increased in the AtERF73/HRE1 overexpression line, which suggested that the Arabidopsis gene AtERF73/HRE1 plays a central role in anaerobic respiration regulation under hypoxia60. Expression of rice genes SNORKEL1 and SNORKEL2 can be induced by ethylene accumulation under a submerged condition, and products of the two genes then triggered remarkable internode elongation via gibberellin, allowing rice to grow in flooding terrain61. In this study, three DEGs, POPTR_0001s10300, POPTR_0017s11470 and POPTR_0019s13330, were also classified into family AP2/ERF and might also play essential roles in poplar responding to root hypoxia stress. Further work is required to determine whether these DEGs regulate flood tolerance of poplar.
Materials and Methods
Plant materials and growth conditions
Two full-sib clones, LS1 (flood-tolerant) and LS2 (flood-susceptible), originating from P. deltoides cv. Lux ex. I-69/55 (flood-tolerant) × P. simonii (flood-susceptible) were used as the materials in the present study1. Four-week-old stem saplings of LS1 and LS2 were transferred from sterile agar cultures into plastic pots (14 cm × 10 cm × 13 cm) containing a mixture of vermiculite and culture soil (1:3, v/v). The soil (pH 6.0) consisted of 2–5% N, P2O5 and K2O, with more than 20% organic matter (dry weight). After establishing growth after transplantation, all plants were watered with Hoagland solution weekly and tap water twice a week. All plants were grown in a room with a temperature of 25 ± 2 °C under a photoperiod of 16 h light and 8 h dark with a light intensity of 40 μmol m−2 s−1 provided by cool white fluorescent lamps and with 70–80% humidity. All plants with a height of 35–40 cm were randomly assigned to one of two water regime treatments: (1) watered (control check, CK) and (2) flooded (FL). The plants were divided into the two groups, and the experimental design is illustrated in Fig. 7. Pots of the control group had four drainage holes in the bottom. The plants were watered using tap water daily to maintain soil moisture at field capacity. The flooded plants were partially submerged by stagnant flooding to a height of 10 cm above the soil surface in tanks63.
After 0, 1, and 7 d of stagnant flooding treatment, 10 trees per treatment per clone were harvested. Fine roots were immediately frozen in liquid nitrogen and kept at −80 °C until utilization for RNA-sequencing (RNA-seq) and metabolism analysis.
For the convenience of description, we defined watered LS1, 1-day flooded LS1, 7-day flooded LS1, and 15-day flooded LS1 followed by 3-day recovery as LS1_CK, LS1_FL1, LS1_FL7 and LS1_FL18, respectively. In LS2, the corresponding definitions were LS2_CK, LS2_FL1, LS2_FL7 and LS2_FL18, respectively. LS1_FL1 and LS1_FL7 composed LS1_FL, and LS2_FL1 and LS2_FL7 composed LS2_FL.
Measurement of leaf gas exchange and chlorophyll fluorescence
The 5th fully expanded and mature leaf from the top of the stem was chosen to measure leaf gas exchange and chlorophyll fluorescence23. Five plants per treatment per clone were measured between 9:00 am. and 11:00 am. using an LI-6400 photosynthesis system (LI-COR Inc., Lincoln, NE, USA) with a standard LI-COR gas exchange chamber. A 1500 μmol m−2 s−1 light intensity of illumination was provided by red diodes (6400–02 LED Source), and the gas flow rate was set as 500 μmol s−1. The gas exchange measurements included Pn and Gs with transpiration rate, intercellular CO2 concentration, atmospheric CO2 concentration and the corresponding ambient environmental conditions, such as temperature, relative humidity and photosynthetically active radiation. Chlorophyll fluorescence (Fv/Fm and Fv/Fo) of leaves of five plants per treatment per clone was measured using an LI-6400 fluorescence system (LI-COR Inc., USA) after a 20-min dark adaptation. To compare the effects of flooding treatment on plants, changes in all growth and ecophysiological parameters under flooding were calculated as follow: change rate = (CK − FL)/CK, where CK represents the value of the control, and FL represents the value in the flooding treatment.
RNA-sequencing
Total RNA was extracted and purified using Trizol (Invitrogen, Canada), and 3 μg of purified RNA per sample was used for library construction. Briefly, the first-strand cDNA was synthesized using a random hexamer primer and M-MuLV reverse transcriptase (RNase H); the second-strand cDNA synthesis was subsequently performed using DNA polymerase I and RNase H; and the cDNA fragments were ligated with adapters. These products were purified and amplified by PCR to create the final cDNA library. The clustering of the index-coded samples was performed on a cBot cluster generation system using TruSeq SR cluster kit v3-cBot-HS (Illumina) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq 2000 platform, and 100 bp single-end reads were generated (Novogene Bioinformatics Institute, China). Two biological replicates per clone per treatment were used in this study, and each biological replicate included five well-mixed individual plants.
RNA-seq data analysis
Raw data were utilized to obtain clean data by removing reads containing adapter sequences, reads containing poly-N and low quality reads. All the downstream analyses were performed based on the clean data. The genome of P. trichocarpa64 was used as the reference genome for sequence read alignment and identification. For unigenes, Bowtie v0.12.9 was used to align single-end clean reads to the unigene sequences. HTSeq v0.5.4p3 was used to count the read numbers mapped to each gene. Then, the reads per kilobase of exon model per million mapped reads (RPKM) of each gene was calculated based on the length of the gene and read counts mapped to the gene65. In this study, we defined genes with RPKM ≥1 as expressed genes. Principal component analysis was performed using RPKM of 12 RNA-seq libraries in the SPSS 19 statistical software package (IBM Co., USA).
Differential expression analysis of the two conditions was performed using the DESeq R package (1.10.1; TNLIST, China)66. The resulting p-values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate67. Genes with an adjusted p-value (padj) <0.05 and |log2-fold change| ≥1 were assigned as differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis of the DEGs was implemented by the GOseq R package in which gene length bias was corrected. GO terms with a corrected p-value < 0.01 were considered significantly enriched by DEGs. KOBAS software (KOBAS, UK) was used to test the statistical enrichment of DEGs in KEGG pathways.
Validation by qRT-PCR
Eighteen DEGs were selected to validate all gene expression results from RNA-seq by qRT-PCR (Supplementary Table 1). Five independent biological repeats were conducted, and RNA was reverse transcribed using PrimescriptTM RT master mix reagent (Takara, China). The qRT-PCR was performed using SYBR® Premix Ex Taq™ II (Takara, China) on a 7500 fast real-time PCR system (Applied Biosystems, NY). Three technical replicates for each gene in all samples were used in qRT-PCR. Some primers were obtained from the literature (Supplementary Table 1), and the others were designed using Premier 5.0 (Premier Biosoft, CA). 18 S rRNA was used as an internal control to normalize all data68, and the relative expression of a gene pair was calculated using the 2−ΔΔCT method and presented by log2-fold change values69.
Gas chromatography-mass spectrometer analysis
The same samples used for RNA-seq were also used for gas chromatography-mass spectrometer (GC-MS) analysis. Metabolites from poplar roots (100 ± 10 mg) were extracted, derived and profiled according to the method used in barley70 with some modifications. Data were analysed using the R software platform (http://cran.r-project.org/) and Xcalibur software. Identification of compounds was based on a comparison of their mass spectra, retention indices (RIs) and retention times with the authentic standards and published data, in addition to the standard mass spectra published in the NIST201471. Ribitol was used as an internal standard for semi-quantification analysis. By comparing the GC-peak area of each compound with that of the internal standard, relative units were used to express the contents of the metabolites72.
Measurement of biochemical parameters
Measurements of physiological-biochemical parameters of fine roots were conducted (four replicates per treatment per clone) to clarify the metabolisms of energy and ROS, two important pathways that respond to abiotic stress in plants. Tissue homogenate was prepared by grinding 0.5 g of fine roots in 4.5 ml of normal saline. The parameters were measured according to the instruction manuals of the reagent kits (Nanjing Jiancheng Bioengineering Institute, China), including contents of soluble protein (Bradford method), GSH (DTNB method) and H2O2 and activities of CAT, GR, POD, SOD (Xanthine oxidase method), ADH and LDH. H2O2 contents were determined spectrophotometrically by monitoring the amount of complex compound (H2O2 and molybdic acid) at a wavelength of 405 nm. CAT activity was measured by monitoring the decrease in the amount of complex compound due to H2O2 decomposition at 405 nm. GR activity was determined by detecting the oxidation of nicotinamide adenine dinucleotide phosphate (NADPH) at 340 nm in the presence of glutathione disulfide and expressed as nmol NADPH oxidized per mg of protein per minute. POD activity was measured following the principle described by Montavon73. ADH activity was detected spectrophotometrically at 340 nm based on nicotinamide adenine dinucleotide (NAD+) reduction catalysed by ADH. In the LDH activity measurement, LDH in tissue homogenate was used to catalyse pyruvic acid generation. Then, pyruvic acid content was measured by the 2,4-dinitrophenylhydrazine method using an automatic microplate reader (Infinite M200, Austria). LDH activity was finally calculated from pyruvic acid production.
Statistical analyses
The data on gas exchange, chlorophyll fluorescence and growth were subjected to analysis of variance (ANOVA). The data of enzyme activities, metabolites of ROS and soluble protein were subjected to Student’s t-test between flooding treatment and control and between LS1 and LS2 at every time point. The data of carbohydrates and amino acids were subjected to Student’s t-test between flooding treatment and control at every time point. All analyses were conducted using the SAS statistical software package version 9.0 (SAS Institute Inc., USA).
Accession code
Sequence reads of transcriptome sequencing have been deposited in the NCBI sequence read archive under accession number SRP127249.
References
Du, K. B., Xu, L., Wu, H., Tu, B. K. & Zheng, B. Ecophysiological and morphological adaption to soil flooding of two poplar clones differing in flood-tolerance. Flora 207, 96–106, https://doi.org/10.1016/j.flora.2011.11.002 (2012).
Gong, J. R., Zhang, X. S., Huang, Y. M. & Zhang, C. L. The effects of flooding on several hybrid poplar clones in Northern China. Agroforest Syst 69, 77–88 (2007).
Štícha, V., Macků, J. & Nuhlíček, O. Effect of permanent waterlogging on the growth of poplar clones MAX 4, MAX 5 (J-104, J-105) (Populus maximowiczii A. Henry × P. nigra Linnaeus) and evaluation of wood moisture content in different stem parts - short communication. Journal of Forest Science (2016).
Kenong, X. et al. SubIA Is an Ethylene Response Factor-Like Gene that Confers Submergence Tolerance to Rice. Nature 442, 705 (2006).
Syed, N. H. et al. Core clock, SUB1, and ABAR genes mediate flooding and drought responses via alternative splicing in soybean. J Exp Bot 66, 7129 (2015).
Lone, A. A., Khan, M. H., Dar, Z. A. & Wani, S. H. Breeding strategies for improving growth and yield under waterlogging conditions in maize: a review. Maydica (2016).
Yin, X., Hiraga, S., Hajika, M., Nishimura, M. & Komatsu, S. Transcriptomic analysis reveals the flooding tolerant mechanism in flooding tolerant line and abscisic acid treated soybean. Plant Mol Biol 93, 479 (2017).
Yin, X. & Komatsu, S. Comprehensive analysis of response and tolerant mechanisms in early-stage soybean at initial-flooding stress. Journal of Proteomics (2017).
Loreti, E., van Veen, H. & Perata, P. Plant responses to flooding stress. Curr Opin Plant Biol 33, 64–71 (2016).
Kreuzwieser, J. & Rennenberg, H. Molecular and physiological responses of trees to waterlogging stress. Plant, cell & environment 37, 2245–2259, https://doi.org/10.1111/pce.12310 (2014).
Reeksting, B. J., Olivier, N. A. & van den Berg, N. Transcriptome responses of an ungrafted Phytophthora root rot tolerant avocado (Persea americana) rootstock to flooding and Phytophthora cinnamomi. Bmc Plant Biol 16 (2016).
Kreuzwieser, J., Furniss, S. & Rennenberg, H. Impact of waterlogging on the N-metabolism of flood tolerant and non-tolerant tree species. Plant Cell and Environment 25, 1039–1049 (2002).
Kreuzwieser, J. et al. Differential Response of Gray Poplar Leaves and Roots Underpins Stress Adaptation during Hypoxia. Plant Physiol 149, 461–473 (2009).
Christianson, J. A., Llewellyn, D. J., Dennis, E. S. & Wilson, I. W. Comparisons of early transcriptome responses to low-oxygen environments in three dicotyledonous plant species. Plant signaling & behavior 5, 1006–1009, https://doi.org/10.1093/pcp/pcp163 (2010).
Le Provost, G. et al. Role of waterlogging-responsive genes in shaping interspecific differentiation between two sympatric oak species. Tree Physiol 32, 119–134 (2012).
Vuleta, A., Jovanovic, S. M. & Tucic, B. Adaptive flexibility of enzymatic antioxidants SOD, APX and CAT to high light stress: The clonal perennial monocot Iris pumila as a study case. Plant Physiol Bioch 100, 166–173 (2016).
Zhou, L. N. et al. Antioxidative Responses and Morpho-anatomical Alterations for Coping with Flood-Induced Hypoxic Stress in Grass Pea (Lathyrus sativus L.) in Comparison with Pea (Pisum sativum). J Plant Growth Regul 35, 690–700 (2016).
Van, V. H. et al. Group VII ethylene response factor diversification and regulation in four species from flood-prone environments. Plant Cell & Environment 37, 2421–2432 (2014).
Yu, B., Zhao, C. Y., Li, J., Li, J. Y. & Peng, G. Morphological, physiological, and biochemical responses of Populus euphratica to soil flooding. Photosynthetica 53, 110–117 (2015).
Peng, Y., Zhou, Z., Tong, R., Hu, X. & Du, K. Anatomy and ultrastructure adaptations to soil flooding of two full-sib poplar clones differing in flood-toleranc. Flora 233, 9 (2017).
Bejaoui, Z., Mguis, K., Abassi, M., Albouchi, A. & Lamhamedi, M. S. Involvement of Carbohydrates in Response to Preconditioning Flooding in Two Clones of Populus deltoides Marsh. × P. nigra L. J Plant Growth Regul 35, 492–503 (2016).
Guo, X. Y., Huang, Z. Y., Xu, A. C. & Zhang, X. S. A comparison of physiological, morphological and growth responses of 13 hybrid poplar clones to flooding. Forestry 84, 1–12 (2011).
Peng, Y. J. et al. Roots play a vital role in flood-tolerance of poplar demonstrated by reciprocal grafting. Flora 208, 479–487 (2013).
Du, K. B., Shen, B. X., Xu, L. & Tu, B. K. Estimation of genetic variances in flood tolerance of poplar and selection of resistant F1 generations. Agroforest Syst 74, 243–257 (2008).
Mielke, M. S., De Almeida, A. A. F., Gomes, F. P., Mangabeira, P. A. O. & Silva, D. D. C. Effects of soil flooding on leaf gas exchange and growth of two neotropical pioneer tree species. New Forest 29, 161–168 (2005).
Yordanova, R. Y. & Popova, L. P. Flooding-induced changes in photosynthesis and oxidative status in maize plants. Acta Physiol Plant 29, 535–541 (2007).
Chen, W. et al. Identification and Comparative Analysis of Differential Gene Expression in Soybean Leaf Tissue under Drought and Flooding Stress Revealed by RNA-Seq. Front Plant Sci 7 (2016).
Qi, X. H., Xu, X. W., Lin, X. J., Zhang, W. J. & Chen, X. H. Identification of differentially expressed genes in cucumber (Cucumis sativus L.) root under waterlogging stress by digital gene expression profile. Genomics 99, 160–168 (2012).
Ren, C. G. et al. Elucidation of the molecular responses to waterlogging in Sesbania cannabina roots by transcriptome profiling. Scientific Reports 7, 9256 (2017).
Yang, Z. et al. Transcriptome Analysis of Two Species of Jute in Response to Polyethylene Glycol (PEG)- induced Drought Stress. Scientific Reports 7 (2017).
Narsai, R., Rocha, M., Geigenberger, P., Whelan, J. & van Dongen, J. T. Comparative analysis between plant species of transcriptional and metabolic responses to hypoxia. New Phytologist 190, 472–487 (2011).
Klok, E. J. et al. Expression Profile Analysis of the Low-Oxygen Response in Arabidopsis Root Cultures. The Plant Cell 14, 14 (2002).
Christopher, M. E. & Good, A. G. Characterization of Hypoxically Inducible Lactate Dehydrogenase in Maize. Plant Physiol 112, 1015–1022 (1996).
Rroberts, J. K. M., Calus, J., Wemmer, D., Walbot, V. & Jardetzky, O. Mechanism of cytoplasmic p H regulation in hypoxic maize root tips and its role in survival under hypoxia. Proceedings of the National Academy of Sciences of the United States of America 81, 5 (1984).
Von Wittgenstein, N. J. J. B. Nitrogen transporters: comparative genomics, transport activity, and gene expression of NRTs and AMTs in Black Cottonwood (Populus trichocarpa). Poplar (2013).
Branco-Price, C., Kaiser, K. A., Jang, C. J. H., Larive, C. K. & Bailey-Serres, J. Selective mRNA translation coordinates energetic and metabolic adjustments to cellular oxygen deprivation and reoxygenation in Arabidopsis thaliana. Plant J 56, 743–755 (2008).
Modolo, L. V., Cunha, F. Q., Braga, M. R. & Salgado, I. Nitric oxide synthase-mediated phytoalexin accumulation in soybean cotyledons in response to the Diaporthe phaseolorum f. sp meridionalis elicitor. Plant Physiol 130, 1288–1297 (2002).
Del Rio, L. A., Corpas, F. J. & Barroso, J. B. Nitric oxide and nitric oxide synthase activity in plants. Phytochemistry 65, 783–792 (2004).
Liu, B., Rennenberg, H. & Kreuzwieser, J. Hypoxia induces stem and leaf nitric oxide (NO) emission from poplar seedlings. Planta 241, 579–589 (2015).
Fritze, I. M. et al. The crystal structures of Zea mays and Arabidopsis 4-hydroxyphenylpyruvate dioxygenase. Plant Physiol 134, 1388–1400 (2004).
Wang, K. H. & Jiang, Y. W. Antioxidant responses of creeping bentgrass roots to waterlogging. Crop Sci 47, 232–238 (2007).
Yordanova, R. Y., Christov, K. N. & Popova, L. P. Antioxidative enzymes in barley plants subjected to soil flooding. Environ Exp Bot 51, 93–101 (2004).
Rajhi, I. et al. Identification of genes expressed in maize root cortical cells during lysigenous aerenchyma formation using laser microdissection and microarray analyses. New Phytologist 190, 351–368 (2011).
Arora, A., Sairam, R. K. & Srivastava, G. C. Oxidative stress and antioxidative system in plants. Curr Sci India 82, 1227–1238 (2002).
Bossi, R. T., Negri, A., Tedeschi, G. & Mattevi, A. Structure of FAD-bound L-aspartate oxidase: Insight into substrate specificity and catalysis. Biochemistry 41, 3018–3024 (2002).
Li, C. Y. et al. Role of beta-oxidation in jasmonate biosynthesis and systemic wound signaling in tomato. Plant Cell 17, 971–986 (2005).
Yoshida, K., Kaothien, P., Matsui, T., Kawaoka, A. & Shinmyo, A. Molecular biology and application of plant peroxidase genes. Appl Microbiol Biot 60, 665–670 (2003).
Yang, T. & Poovaiah, B. W. Hydrogen peroxide homeostasis: Activation of plant catalase by calcium/calmodulin. Proceedings of the National Academy of Sciences of the United States of America 99, 4097–4102 (2002).
Caverzan, A. et al. Plant responses to stresses: Role of ascorbate peroxidase in the antioxidant protection. Genet Mol Biol 35, 1011–1019 (2012).
Noctor, G. et al. Glutathione: biosynthesis, metabolism and relationship to stress tolerance explored in transformed plants. J Exp Bot 49, 623–647 (1998).
Polle, A. Dissecting the superoxide dismutase-ascorbate-glutathione-pathway in chloroplasts by metabolic modeling. Computer simulations as a step towards flux analysis. Plant Physiol 126, 445–462 (2001).
Bailey-Serres, J. & Voesenek, L. A. C. J. Flooding stress: Acclimations and genetic diversity. Annu Rev Plant Biol 59, 313–339 (2008).
Carnus, J. M. et al. Planted forests and biodiversity. J Forest 104, 65–77 (2006).
Poirier, Y., Antonenkov, V. D., Glumoff, T. & Hiltunen, J. K. Peroxisomal beta-oxidation - A metabolic pathway with multiple functions. Bba-Mol Cell Res 1763, 1413–1426 (2006).
Bailey-Serres, J. et al. Making sense of low oxygen sensing. Trends Plant Sci 17, 129–138 (2012).
Raineri, J., Ribichich, K. F. & Chan, R. L. The sunflower transcription factor HaWRKY76 confers drought and flood tolerance to Arabidopsis thaliana plants without yield penalty. Plant Cell Rep 34, 2065–2080 (2015).
Rauf, M. et al. NAC Transcription Factor SPEEDY HYPONASTIC GROWTH Regulates Flooding-Induced Leaf Movement in Arabidopsis. Plant Cell 25, 4941–4955 (2013).
Banti, V., Mafessoni, F., Loreti, E., Alpi, A. & Perata, P. The Heat-Inducible Transcription Factor HsfA2 Enhances Anoxia Tolerance in Arabidopsis. Plant Physiol 152, 1471–1483 (2010).
Jung, K. H. et al. The Submergence Tolerance Regulator Sub1A Mediates Stress-Responsive Expression of AP2/ERF Transcription Factors. Plant Physiol 152, 1674–1692 (2010).
Yang, C. Y., Hsu, F. C., Li, J. P., Wang, N. N. & Shih, M. C. The AP2/ERF Transcription Factor AtERF73/HRE1 Modulates Ethylene Responses during Hypoxia in Arabidopsis. Plant Physiol 156, 202–212 (2011).
Hattori, Y. et al. The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature 460, 1026–U1116 (2009).
Qi, W. W. et al. Rice Ethylene-Response AP2/ERF Factor OsEATB Restricts Internode Elongation by Down-Regulating a Gibberellin Biosynthetic Gene. Plant Physiol 157, 216–228 (2011).
Sasidharan, R. et al. Community recommendations on terminology and procedures used in flooding and low oxygen stress research. New Phytologist 214, 1403–1407 (2017).
Tuskan, G. A. et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313, 1596–1604 (2006).
Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L. & Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods 5, 621–628, https://doi.org/10.1038/nmeth.1226 (2008).
Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol 11 (2010).
Benjamini, Y. & Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society 57, 12 (1995).
Xu, M., Zhang, B., Su, X. H., Zhang, S. G. & Huang, M. R. Reference gene selection for quantitative real-time polymerase chain reaction in Populus. Anal Biochem 408, 337–339 (2011).
Livak, K. J. & Schmittgen, T. D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−delta delta CT Method. Methods 25, 7 (2001).
Wu, D. Z. et al. Tissue Metabolic Responses to Salt Stress in Wild and Cultivated Barley. Plos One 8 (2013).
Cai, X. et al. Analysis of aroma-active compounds in three sweet osmanthus (Osmanthus fragrans) cultivars by GC-olfactometry and GC-MS. J Zhejiang Univ-Sc B 15, 638–648 (2014).
Kaseleht, K., Leitner, E. & Paalme, T. Determining aroma-active compounds in Kama flour using SPME-GC/MS and GC-olfactometry. Flavour Frag J 26, 122–128 (2011).
Montavon, P. & Bortlik, K. Evolution of robusta green coffee redox enzymatic activities with maturation. Journal of agricultural and food chemistry 52, 3590–3594, https://doi.org/10.1021/jf0308004 (2004).
Acknowledgements
This study was supported by the National Natural Science Foundation of China (Program Nos 31570665 and 31100458), Major technical innovation project of Hubei Province (Program No. 2016ABA111), and Fundamental Research Funds for the Central Universities (Program Nos 201710504036 and 201310504023). The authors are grateful to AJE (American Journal Experts) for improving the English editing in this paper.
Author information
Authors and Affiliations
Contributions
K.B.D. and Y.J.P. conceived and designed the study. Y.J.P., Z.X.Z. and K.B.D. analyzed the data and drafted the manuscript. Y.J.P., X.L.Y., Z.Z. and X.Y.Z. carried out the experiments. All authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Peng, Y., Zhou, Z., Zhang, Z. et al. Molecular and physiological responses in roots of two full-sib poplars uncover mechanisms that contribute to differences in partial submergence tolerance. Sci Rep 8, 12829 (2018). https://doi.org/10.1038/s41598-018-30821-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-30821-y
This article is cited by
-
Yield, cell structure and physiological and biochemical characteristics of rapeseed under waterlogging stress
BMC Plant Biology (2024)
-
Physiological and transcriptional regulation in Taxodium hybrid ‘Zhongshanshan’ leaves in acclimation to prolonged partial submergence
Planta (2023)
-
Distinctive Characteristics of Nitrogen Metabolism in Teosinte Versus Modern Maize for Water-Logging Tolerance
Journal of Plant Growth Regulation (2023)
-
Comparative transcriptomic analysis uncovers conserved pathways involved in adventitious root formation in poplar
Physiology and Molecular Biology of Plants (2021)