Joint transcriptomic and metabolomic analysis provides new insights into drought resistance in watermelon (Citrullus lanatus)

Introduction: Watermelon is an annual vine of the family Cucurbitaceae. Watermelon plants produce a fruit that people love and have important nutritional and economic value. With global warming and deterioration of the ecological environment, abiotic stresses, including drought, have become important factors that impact the yield and quality of watermelon plants. Previous research on watermelon drought resistance has included analyzing homologous genes based on known drought-responsive genes and pathways in other species. Methods: However, identifying key pathways and genes involved in watermelon drought resistance through high-throughput omics methods is particularly important. In this study, RNA-seq and metabolomic analysis were performed on watermelon plants at five time points (0 h, 1 h, 6 h, 12 h and 24 h) before and after drought stress. Results: Transcriptomic analysis revealed 7829 differentially expressed genes (DEGs) at the five time points. The DEGs were grouped into five clusters using the k-means clustering algorithm. The functional category for each cluster was annotated based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; different clusters were associated with different time points after stress. A total of 949 metabolites were divided into 10 categories, with lipids and lipid-like molecules accounting for the most metabolites. Differential expression analysis revealed 22 differentially regulated metabolites (DRMs) among the five time points. Through joint analysis of RNA-seq and metabolome data, the 6-h period was identified as the critical period for watermelon drought resistance, and the starch and sucrose metabolism, plant hormone signal transduction and photosynthesis pathways were identified as important regulatory pathways involved in watermelon drought resistance. In addition, 15 candidate genes associated with watermelon drought resistance were identified through joint RNA-seq and metabolome analysis combined with weighted correlation network analysis (WGCNA). Four of these genes encode transcription factors, including bHLH (Cla97C03G068160), MYB (Cla97C01G002440), HSP (Cla97C02G033390) and GRF (Cla97C02G042620), one key gene in the ABA pathway, SnRK2-4 (Cla97C10G186750), and the GP-2 gene (Cla97C05G105810), which is involved in the starch and sucrose metabolism pathway. Discussion: In summary, our study provides a theoretical basis for elucidating the molecular mechanisms underlying drought resistance in watermelon plants and provides new genetic resources for the study of drought resistance in this crop.