Abstract
Abiotic stress negatively impacts cassava (Manihot esculenta) growth and yield. Several molecular mechanisms of plant response to cold and drought have been identified and described in the literature, however, little is known about the crosstalk of the responses of cassava to these two stresses. To elucidate this question, transcriptome analysis of cassava seedlings under cold or PEG-simulated drought stress treatment was performed. Our results showed that 6103 and 7462 transcripts were significantly regulated by cold and drought stress, respectively. Gene Ontology annotation revealed that the abscisic and jasmonic acid signaling pathways shared between the two stresses responses. We further identified 2434 common differentially expressed genes (DEGs), including 1130 up-regulated and 841 down-regulated DEGs by the two stresses. These co-induced or co-suppressed genes are grouped as stress signal perception and transduction, transcription factors (TFs), metabolism as well as transport facilitation according to the function annotation. Furthermore, a large proportion of well characterized protein kinases, TF families and ubiquitin proteasome system related genes, such as RLKs, MAPKs, AP2/ERFBPs, WRKYs, MYBs, E2 enzymes and E3 ligases, including three complexes of interacting proteins were shown as key points of crosstalk between cold and drought stress signaling transduction pathways in a hierarchical manner. Our research provides valuable information and new insights for genetically improving the tolerance of crops to multiple abiotic stresses.