Weighted Gene Co-Expression Network Analysis Uncovers Core Drought Responsive Genes in Pecan (Carya illinoinensis)

Drought severely affects the growth and production of pecan (Carya illinoinensis), while genes conferred drought adaptation are yet to be fully elucidated. Here, an in-depth exploration of the two different RNA-seq projects regarding drought stress (designated as P1 and P2) was performed via weighted gene co-expression network analysis. For the two projects, there existed one pair of modules (P1 turquoise module and P2 blue module) that was probably associated with drought resistance, as the paired modules both exhibited an increased expression profile with increasing water shortage stress and were annotated to be involved in oxidative stress response and the signaling pathways of abscisic acid and jasmonic acid. There were 441 and 1258 hub genes in the P1 turquoise module and P2 blue module, respectively, among which, 140 were overlapped and thus were recognized as core drought responsive genes. An additional drought stress experiment was conducted for RT-qPCR validation, and the results showed that the 20 core genes selected for detection were highly responsive to water deficit. Together, our results will be helpful for understanding the molecular mechanism of drought response and improving drought resistance in pecan.