Transcriptome Analysis of Pecan (Carya illinoinensis) Differentially Expressed Genes in Response to Drought Stress

Abstract

Pecan (Carya illinoinensis) is an economically important nut tree that is greatly affected by drought, limiting its production and distribution. Although the complete genome of the pecan was published several years ago, the molecular basis of the pecan’s response to drought remains unclear. In this study, we analyzed the high-throughput transcriptome data for pecans under 3, 6, 9, 12, and 15 days of drought stress compared with the controls. A total of 12,893 differentially expressed genes (DEGs) were identified under drought stress, with 11,684 of them showing significant changes after 15 d of drought treatment. Among these, 4448 genes were up-regulated while 7226 were down-regulated. The trend analysis revealed that DEGs could be classified into 20 clusters. Surprisingly, the majority of genes (6148) showed a gradual down-regulation, and 3683 genes showed a gradual up-regulation in response to drought. Gene ontology enrichment analysis showed that the DEGs were mainly enriched in biological processes. The KEGG pathway enrichment results indicated that the DEGs were mainly enriched in several pathways, including metabolic pathways, the biosynthesis of secondary metabolites, and plant hormone signal transduction processes. Among the DEGs, 457 protein kinase and 734 transcription factor genes were shown to be drought-responsive and may play key roles in the response to drought, and the expression patterns of selected candidate genes were further validated using quantitative real-time PCR. Collectively, these findings highlighted the multiple processes in pecans under drought stress and provided valuable insights into the further investigation of the functions of drought stress-responsive genes and the molecular basis of the pecan drought stress response.