Selecting suitable internal reference genes is an important prerequisite for the application of quantitative real-time PCR (qRT-PCR). However, no systematic studies have been conducted on reference genes in luffa. In this study, seven reference genes were selected, and their expression levels in luffa plants exposed to various simulated abiotic stresses [i.e., cold, drought, heat, salt, H₂O₂, and abscisic acid (ABA) treatments] were analyzed by qRT-PCR. The stability of the reference gene expression levels was validated using the geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The results indicated thatEF-1αwas the most stably expressed and suitable reference gene overall and for the heat, cold, and ABA treatments. Additionally,UBQexpression was stable following the salt treatment, whereasTUBwas identified as a suitable reference gene for H₂O₂and drought treatments. The reliability of the selected reference genes was verified by analyzing the expression of copper/zinc superoxide dismutase (Cu/Zn-SOD) gene in luffa. When the most unstable reference genes were used for data normalizations, the resulting expression patterns had obvious biases when compared with the expression patterns for the most ideal reference genes used alone or combined. These results will be conducive to more accurate quantification of gene expression levels in luffa.

Fresh-cut luffa (Luffa cylindrica) fruits commonly undergo browning. However, little is known about the molecular mechanisms regulating this process. We used the RNA-seq technique to analyze the transcriptomic changes occurring during the browning of fresh-cut fruits from luffa cultivar ‘Fusi-3’. Over 90 million high-quality reads were assembled into 58,073 Unigenes, and 60.86% of these were annotated based on sequences in four public databases. We detected 35,282 Unigenes with significant hits to sequences in the NCBInr database, and 24,427 Unigenes encoded proteins with sequences that were similar to those of known proteins in the Swiss-Prot database. Additionally, 20,546 and 13,021 Unigenes were similar to existing sequences in the Eukaryotic Orthologous Groups of proteins and Kyoto Encyclopedia of Genes and Genomes databases, respectively. Furthermore, 27,301 Unigenes were differentially expressed during the browning of fresh-cut luffa fruits (i.e., after 1–6 h). Moreover, 11 genes from five gene families (i.e., PPO, PAL, POD, CAT, and SOD) identified as potentially associated with enzymatic browning as well as four WRKY transcription factors were observed to be differentially regulated in fresh-cut luffa fruits. With the assistance of rapid amplification of cDNA ends technology, we obtained the full-length sequences of the 15 Unigenes. We also confirmed these Unigenes were expressed by quantitative real-time polymerase chain reaction analysis. This study provides a comprehensive transcriptome sequence resource, and may facilitate further studies aimed at identifying genes affecting luffa fruit browning for the exploitation of the underlying mechanism.

The plant-specific WRKY transcription factor family members have diverse regulatory effects on the genes associated with many plant processes. Although the WRKY proteins in Arabidopsis thaliana and other species have been thoroughly investigated, there has been relatively little research on the WRKY family in Luffa cylindrica, which is one of the most widely grown vegetables in China. In this study, we performed a genome-wide analysis to identify L. cylindrica WRKY genes, which were subsequently classified and examined in terms of their gene structures, chromosomal locations, promoter cis-acting elements, and responses to abiotic stress. A total of 62 LcWRKY genes (471–2238 bp) were identified and divided into three phylogenetic groups (I, II, and III), with group II further divided into five subgroups (IIa, IIb, IIc, IId, and IIe) in accordance with the classification in other plants. The LcWRKY genes were unevenly distributed across 13 chromosomes. The gene structure analysis indicated that the LcWRKY genes contained 0–11 introns (average of 4.4). Moreover, 20 motifs were detected in the LcWRKY proteins with conserved motifs among the different phylogenetic groups. Two subgroup IIc members (LcWRKY16 and LcWRKY31) contained the WRKY sequence variant WRKYGKK. Additionally, nine cis-acting elements related to diverse responses to environmental stimuli were identified in the LcWRKY promoters. The subcellular localization analysis indicated that three LcWRKY proteins (LcWRKY43, LcWRKY7, and LcWRKY23) are localized in the nucleus. The tissue-specific LcWRKY expression profiles reflected the diversity in LcWRKY expression. The RNA-seq data revealed the effects of low-temperature stress on LcWRKY expression. The cold-induced changes in expression were verified via a qRT-PCR analysis of 24 differentially expressed WRKY genes. Both LcWRKY7 and LcWRKY12 were highly responsive to the low-temperature treatment (approximately 110-fold increase in expression). Furthermore, the LcWRKY8, LcWRKY12, and LcWRKY59 expression levels increased by more than 25-fold under cold conditions. Our findings will help clarify the evolution of the luffa WRKY family while also providing valuable insights for future studies on WRKY functions.

Plant-specific APETALA2/Ethylene-Responsive Factor (AP2/ERF) transcription factors are involved in the regulation of genes associated with the growth and developmental processes of numerous plants. Although AP2/ERF proteins from other species have been intensively studied, no studies have been reported on the AP2/ERF family of Luffa cylindrica, an important vegetable of the cucurbit family, and one of the most popular vegetables in the world. In this study, 133 genes (315–6696 bp) encoding LcAP2/ERF proteins with complete AP2/ERF domains were identified according to the luffa P93075 genome. These LcAP2/ERF genes were subsequently classified and analyzed for their gene structures, chromosomal distribution locations, promoter cis-acting elements, conserved structural domains of encoded proteins, and responses to abiotic stresses. The LcAP2/ERF genes were identified and divided into five phylogenetic groups (AP2, DREBs, ERFs, RAV, and soloists). These genes were unevenly distributed across 13 chromosomes. An analysis of gene structures indicated the LcAP2/ERF genes contained 0–11 introns (average of 4.4). Additionally, 16 motifs were identified in the LcAP2/ERF proteins that were conserved across different phylogenetic groups. Moreover, 11 cis-acting elements associated with response to the environment were analyzed in a 2000 bp region upstream of the LcAP2/ERF gene promoters. A transcriptome analysis involving RNA-seq data revealed tissue-specific LcAP2/ERF expression profiles and the diversity in LcAP2/ERF expression. The effects of low-temperature stress on LcAP2/ERF expression were determined. Furthermore, fruit-development-related and low-temperature-induced expressional changes were verified by RT-qPCR analyses of 14 differentially expressed LcAP2/ERF genes in luffa. Our findings will help clarify the evolution of the luffa AP2/ERF family, while also providing valuable insights for future studies on AP2/ERF functions.

Cyclic nucleotide-gated ion channels (CNGCs) are cell membrane channel proteins for calcium ions. They have been reported to play important roles in survival and in the responses to environmental factors in various plants. However, little is known about the CNGC family and its functions in luffa (Luffa cylindrica L.). In this study, a bioinformatics-based method was used to identify members of the CNGC gene family in L. cylindrica. In total, 20 LcCNGCs were detected, and they were grouped into five subfamilies (I, II, Ⅲ, IV-a, and IV-b) in a phylogenetic analysis with CNGCs from Arabidopsis thaliana (20 AtCNGCs) and Momordica charantia (17 McCNGCs). The 20 LcCNGC genes were unevenly distributed on 11 of the 13 chromosomes in luffa, with none on Chromosomes 1 and 5. The members of each subfamily encoded proteins with highly conserved functional domains. An evolutionary analysis of CNGCs in luffa revealed three gene losses and a motif deletion. An examination of gene replication events during evolution indicated that two tandemly duplicated gene pairs were the primary driving force behind the evolution of the LcCNGC gene family. PlantCARE analyses of the LcCNGC promoter regions revealed various cis-regulatory elements, including those responsive to plant hormones (abscisic acid, methyl jasmonate, and salicylic acid) and abiotic stresses (light, drought, and low temperature). The presence of these cis-acting elements suggested that the encoded CNGC proteins may be involved in stress responses, as well as growth and development. Transcriptome sequencing (RNA-seq) analyses revealed tissue-specific expression patterns of LcCNGCs in various plant parts (roots, stems, leaves, flowers, and fruit) and the upregulation of some LcCNGCs under low-temperature stress. To confirm the accuracy of the RNA-seq data, 10 cold-responsive LcCNGC genes were selected for verification by quantitative real-time polymerase chain reaction (RT-qPCR) analysis. Under cold conditions, LcCNGC4 was highly upregulated (>50-fold increase in its transcript levels), and LcCNGC3, LcCNGC6, and LcCNGC13 were upregulated approximately 10-fold. Our findings provide new information about the evolution of the CNGC family in L. cylindrica and provide insights into the functions of the encoded CNGC proteins.

Cauliflower is one of the most important popular vegetables in China. Cauliflower (Brassica oleracea L. var. botrytis) produces glucosinolates, which are secondary metabolites that have anti-cancer properties. The choice of suitable reference ...