Comparative omics-based characterization, phylogeny and melatonin-mediated expression analyses of GDSL genes in pitaya (Selenicereus undatus L.) against multifactorial abiotic stresses

The GDSL gene family plays diverse roles in plant growth and development. Despite its significance, the functions of the GDSL in the pitaya plant are still unknown. Pitaya (Selenicereus undatus L.) also called Hylocereus undatus (Hu), belongs to the family Cactaceae and is an important tropical plant that contains high dietary fibers and antioxidants. In the present investigation, we screened 91 HuGDSL genes in the pitaya genome by conducting a comprehensive computational analysis. The phylogenetic tree categorized HuGDSL genes into 9 distinct clades in combination with four other species. Further, 29 duplicate events were identified of which 12 were tandem, and 17 were segmental. The synteny analysis revealed that segmental duplication was more prominent than tandem duplication among these genes. The majority of duplicated gene pairs (95%) indicate their Ka/Ks ratios ranging from 0.1 to 0.3, which shows that maximum HuGDSL genes were under purifying selection pressure. The cis-acting element in the promotor region contains phytohormones such as auxin, gibberellin, jasmonic acid, and abscisic acid abundantly. Finally, the HuGDSL gene expression pattern under single and multiple stresses was analyzed via; RNA-seq. We select ten stress-responsive HuGDSL genes for RT-qPCR validation. After careful investigation, we identified five HuGDSL candidate genes (HuGDSL-1/3/55/59, and HuGDSL-78) based on RNA-seq, and RT-qPCR data that showed enhanced expression in stress and melatonin-applied seedlings. This study represents valuable insights into maintaining pitaya growth and development by preparing stress-resilient pitaya genotypes through modern biotechnological techniques.