Comparative transcriptomics analysis of tolerant and sensitive genotypes reveals genes involved in the response to cold stress in bitter gourd (Momordica charantia L.).

Bitter gourd is an economically important horticultural crop for its edible and medicinal value. However, the regulatory mechanisms of bitter gourd in response to cold stress are still poorly elucidated. In this study, phytohormone determination and comparative transcriptome analyses in XY (cold-tolerant) and QF (cold-sensitive) after low temperature treatment were conducted. Under cold stress, the endogenous contents of abscisic acid (ABA), jasmonic acid (JA) and salicylic acid (SA) in XY were significantly increased at 24 h after treatment (HAT), indicating that ABA, JA and SA might function in regulating cold resistance. RNA-seq results revealed that more differentially expressed genes were identified at 6 HAT in QF and 24 HAT in XY, respectively. KEGG analysis suggested that the plant hormone signal transduction pathway was significantly enriched in both genotypes at all the time points. In addition, transcription factors showing different expression patterns between XY and QF were identified, including CBF3, ERF2, NAC90, WRKY51 and WRKY70. Weighted gene co-expression network analysis suggested MARK1, ERF17, UGT74E2, GH3.1 and PPR as hub genes. These results will deepen the understanding of molecular mechanism of bitter gourd in response to cold stress and the identified genes may help to facilitate the genetic improvement of cold-resistant cultivars.