Drought is the main abiotic stress that constrains sugarcane growth and production. To understand the molecular mechanisms that govern drought stress, we performed a comprehensive comparative analysis of physiological changes and transcriptome dynamics related to drought stress of highly drought-resistant (ROC22, cultivated genotype) and weakly drought-resistant (Badila, wild genotype) sugarcane, in a time-course experiment (0 h, 4 h, 8 h, 16 h and 32 h). Physiological examination reviewed that ROC22, which shows superior drought tolerance relative to Badila, has high performance photosynthesis and better anti-oxidation defenses under drought conditions. The time series dataset enabled the identification of important hubs and connections of gene expression networks. We identified 36,956 differentially expressed genes (DEGs) in response to drought stress. Of these, 15,871 DEGs were shared by the two genotypes, and 16,662 and 4423 DEGs were unique to ROC22 and Badila, respectively. Abscisic acid (ABA)-activated signaling pathway, response to water deprivation, response to salt stress and photosynthesis-related processes showed significant enrichment in the two genotypes under drought stress. At 4 h of drought stress, ROC22 had earlier stress signal transduction and specific up-regulation of the processes response to ABA, L-proline biosynthesis and MAPK signaling pathway–plant than Badila. WGCNA analysis used to compile a gene regulatory network for ROC22 and Badila leaves exposed to drought stress revealed important candidate genes, including several classical transcription factors: NAC87, JAMYB, bHLH84, NAC21/22, HOX24 and MYB102, which are related to some antioxidants and trehalose, and other genes. These results provide new insights and resources for future research and cultivation of drought-tolerant sugarcane varieties.

Erianthus arundinaceus is a valuable gene reservoir for sugarcane improvement. However, insufficient molecular markers for high-accuracy identification and tracking of the introgression status of E. arundinaceus chromatin impede sugarcane breeding. Fortunately, suppression subtractive hybridization (SSH) technology provides an excellent opportunity for the development of high-throughput E. arundinaceus-specific molecular markers at a reasonable cost. In this study, we constructed a SSH library of E. arundinaceus. In total, 288 clones of E. arundinaceus-specific repetitive sequences were screened out and their distribution patterns on chromosomes were characterized by fluorescence in situ hybridization (FISH). A subtelomeric repetitive sequence Ea086 and a diffusive repetitive sequence Ea009, plus 45S rDNA-bearing E. arundinaceus chromosome repetitive sequence EaITS were developed as E. arundinaceus-specific molecular markers, namely, Ea086-128, Ea009-257, and EaITS-278, covering all the E. arundinaceus chromosomes for high-accuracy identification of putative progeny. Both Ea086-128 and Ea009-257 were successfully applied to identify the authenticity of F₁, BC₁, BC₂, BC₃, and BC₄ progeny between sugarcane and E. arundinaceus. In addition, EaITS-278 was a 45S rDNA-bearing E. arundinaceus chromosome-specific molecular marker for rapid tracking of the inherited status of this chromosome in a sugarcane background. Three BC₃ progeny had apparently lost the 45S rDNA-bearing E. arundinaceus chromosome. We reported herein a highly effective and reliable SSH-based technology for discovery of high-throughput E. arundinaceus-specific sequences bearing high potential as molecular markers. Given its reliability and savings in time and efforts, the method is also suitable for development of species-specific molecular markers for other important wild relatives to accelerate introgression of wild relatives into sugarcane.