Effects of drought and rehydration on root gene expression in seedlings of Pinus massoniana Lamb

The mechanisms underlying plant response to drought involve the expression of numerous functional and regulatory genes. Transcriptome sequencing based on the second- and/or third-generation high-throughput sequencing platforms has proven to be powerful for investigating the transcriptional landscape under drought stress. However, the full-length transcriptomes related to drought responses in the important conifer genus Pinus L. remained to be delineated using the third-generation sequencing technology. With the objectives of identifying the candidate genes responsible for drought and/or rehydration and clarifying the expression profile of key genes involved in drought regulation, we combined the third- and second-generation sequencing techniques to perform transcriptome analysis on seedling roots under drought stress and rewatering in the drought-tolerant conifer Pinus massoniana Lamb. A sum of 294,114 unique full-length transcripts were produced with a mean length of 3217 bp and N50 estimate of 5075 bp, including 279,560 and 124,438 unique full-length transcripts being functionally annotated and Gene Ontology enriched, respectively. A total of 4076, 6295 and 18,093 differentially expressed genes (DEGs) were identified in three pair-wise comparisons of drought-treatment versus control transcriptomes, including 2703, 3576 and 8273 upregulated and 1373, 2719 and 9820 downregulated DEGs, respectively. Moreover, 157, 196 and 691 DEGs were identified as transcription factors in the three transcriptome comparisons and grouped into 26, 34 and 44 transcription factor families, respectively. Gene Ontology enrichment analysis revealed that a remarkable number of DEGs were enriched in soluble sugar-related and cell wall-related processes. A subset of 75, 68 and 97 DEGs were annotated to be associated with starch, sucrose and raffinose metabolism, respectively, while 32 and 70 DEGs were associated with suberin and lignin biosynthesis, respectively. Weighted gene co-expression network analysis revealed modules and hub genes closely related to drought and rehydration. This study provides novel insights into root transcriptomic changes in response to drought dynamics in Masson pine and serves as a fundamental work for further molecular investigation on drought tolerance in conifers.