Integrated analysis of transcriptome and small RNAome reveals regulatory network of rapid and long-term response to heat stress in Rhododendron moulmainense

The Rhododendron plays an important role in maintaining ecological balance. However, it is difficult to domesticate for use in urban ecosystems due to their strict optimum growth temperature condition, and its evolution and adaptation are little known. Here, we combined transcriptome and small RNAome to reveal the rapid response and long-term adaptability regulation strategies in Rhododendron moulmainense under high-temperature stress. The post-transcriptional gene regulatory pathway plays important roles in stress response, in which the protein folding pathway is rapidly induced at 4 h after heat stress, and alternative splicing plays an important role in regulating gene expression at 7 days after heat stress. The chloroplasts oxidative damage is the main factor inhibiting photosynthesis efficiency. Through WGCNA analysis, we identified gene association patterns and potential key regulatory genes responsible for maintaining the ROS steady-state under heat stress. Finally, we found that the sRNA synthesis pathway is induced under heat stress. Combined with small RNAome, we found that more miRNAs are significantly changed under long-term heat stress. Furthermore, MYBs might play a central role in target gene interaction network of differentially expressed miRNAs in R. moulmainense under heat stress. MYBs are closely related to ABA, consistently, ABA synthesis and signaling pathways are significantly inhibited, and the change in stomatal aperture is not obvious under heat stress. Taken together, we gained valuable insights into the transplantation and long-term conservation domestication of Rhododendron, and provide genetic resources for genetic modification and molecular breeding to improve heat resistance in Rhododendron.