Comparative Transcriptome Analysis of Two Sugarcane Cultivars in Response to Paclobutrazol Treatment

Sugarcane is an important crop across the globe, and the rapid multiplication of excellent cultivars is an important object of the sugarcane industry. As one of the plant growth regulators, paclobutrazol (PBZ) has been frequently used in the tissue culture of sugarcane seedlings. However, little is known about the molecular mechanisms of response to PBZ in this crop. Here, we performed a comparative transcriptome analysis between sensitive (LC05-136) and non-sensitive (GGZ001) sugarcane cultivars treated by PBZ at three time points (0 d, 10 d, and 30 d) using RNA sequencing (RNA-Seq). The results showed that approximately 70.36 Mb of clean data for each sample were generated and assembled into 239,212 unigenes. A total of 6108 and 4404 differentially expressed genes (DEGs) were identified within the sensitive and non-sensitive sugarcane cultivars, respectively. Among them, DEGs in LC05-136 were most significantly enriched in the photosynthesis and valine, leucine and isoleucine degradation pathways, while in GGZ001, DEGs associated with ion channels and plant-pathogen interaction were mainly observed. Notably, many interesting genes, including those encoding putative regulators, key components of photosynthesis, amino acids degradation and glutamatergic synapse, were identified, revealing their importance in the response of sugarcane to PBZ. Furthermore, the expressions of sixteen selected DEGs were tested by quantitative reverse transcription PCR (RT-qPCR), confirming the reliability of the RNA-seq data used in this study. These results provide valuable information regarding the transcriptome changes in sugarcane treated by PBZ and provide an insight into understanding the molecular mechanisms underlying the resistance to PBZ in sugarcane.