Global Responses of Autopolyploid Sugarcane Badila (Saccharum officinarum L.) to Drought Stress Based on Comparative Transcriptome and Metabolome Profiling

Sugarcane (Saccharum spp. hybrid) is frequently affected by seasonal drought, which causes substantial declines in quality and yield. To understand the drought resistance mechanisms of S. officinarum, the main species of modern sugarcane, at a molecular level, we carried out a comparative analysis of transcriptome and metabolome profiling of the sugarcane variety Badila under drought stress (DS). Compared with control group (CG) plants, plants exposed to DS had 13,744 (6663 up-regulated and 7081 down-regulated) differentially expressed genes (DEGs). GO and KEGG analysis showed that the DEGs were enriched in photosynthesis-related pathways and most DEGs had down-regulated expression. Moreover, the chlorophyll content, photosynthesis (Photo), stomatal conductance (Cond), intercellular carbon dioxide concentration (Ci) and transpiration rate (Trmmol) were sharply decreased under DS. These results indicate that DS has a significant negative influence on photosynthesis in sugarcane. Metabolome analysis identified 166 (37 down-regulated and 129 up-regulated) significantly regulated metabolites (SRMs). Over 50% of SRMs were alkaloids, amino acids and their derivatives, and lipids. The five most significantly enriched KEGG pathways among SRMs were Aminoacyl-tRNA biosynthesis, 2-Oxocarboxylic acid metabolism, Biosynthesis of amino acids, Phenylalanine metabolism, and Arginine and proline metabolism (p < 0.05). Comparing CG with DS for transcriptome and metabolome profiling (T_CG/DS and M_CG/DS, respectively), we found three of the same KEGG-enriched pathways, namely Biosynthesis of amino acids, Phenylalanine metabolism and Arginine and proline metabolism. The potential importance of Phenylalanine metabolism and Arginine and proline metabolism was further analyzed for response to DS in sugarcane. Seven SRMs (five up-regulated and two down-regulated) and 60 DEGs (17 up-regulated and 43 down-regulated) were enriched in Phenylalanine metabolism under DS, of which novel.31261, Sspon.04G0008060-1A, Sspon.04G0008060-2B and Sspon.04G0008060-3C were significantly correlated with 7 SRMs. In Arginine and proline metabolism, eight SRMs (seven up-regulated and one down-regulated) and 63 DEGs (32 up-regulated and 31 down-regulated) were enriched, of which Sspon.01G0026110-1A (OAT) and Sspon.03G0002750-3D (P5CS) were strongly associated with proline (r > 0.99). These findings present the dynamic changes and possible molecular mechanisms of Phenylalanine metabolism as well as Arginine and proline metabolism under DS and provide a foundation for future research and sugarcane improvement.