Transcriptome Profiles Reveal Key Regulatory Networks during Single and Multifactorial Stresses Coupled with Melatonin Treatment in Pitaya (Selenicereus undatus L.)

In response to evolving climatic conditions, plants frequently confront multiple abiotic stresses, necessitating robust adaptive mechanisms. This study focuses on the responses of <i>Selenicereus undatus</i> L. to both individual stresses (cadmium; Cd, salt; S, and drought; D) and their combined applications, with an emphasis on evaluating the mitigating effects of (M) melatonin. Through transcriptome analysis, this study identifies significant gene expression changes and regulatory network activations. The results show that stress decreases pitaya growth rates by 30%, reduces stem and cladode development by 40%, and increases Cd uptake under single and combined stresses by 50% and 70%, respectively. Under stress conditions, enhanced activities of H₂O₂, POD, CAT, APX, and SOD and elevated proline content indicate strong antioxidant defenses. We identified 141 common DEGs related to stress tolerance, most of which were related to AtCBP, ALA, and CBP pathways. Interestingly, the production of genes related to signal transduction and hormones, including abscisic acid and auxin, was also significantly induced. Several calcium-dependent protein kinase genes were regulated during M and stress treatments. Functional enrichment analysis showed that most of the DEGs were enriched during metabolism, MAPK signaling, and photosynthesis. In addition, weighted gene co-expression network analysis (WGCNA) identified critical transcription factors (WRKYs, MYBs, bZIPs, bHLHs, and NACs) associated with antioxidant activities, particularly within the salmon module. This study provides morpho-physiological and transcriptome insights into pitaya's stress responses and suggests molecular breeding techniques with which to enhance plant resistance.