Fruit expansion is an essential and very complex biological process. Regulatory roles of microRNAs (miRNAs) and miRNA–mRNA modules in the cucumber fruit expansion are not yet to be investigated. In this work, 1253 known and 1269 novel miRNAs were identified from nine cucumber fruit small RNA (sRNA) libraries through high-throughput sequencing. A total of 105 highly differentially expressed miRNAs were recognized in the fruit on five days post anthesis with pollination (EXP_5d) sRNA library. Further, expression patterns of 11 differentially expressed miRNAs were validated by quantitative real-time PCR (qRT-PCR). The expression patterns were similar to sRNAs sequencing data. Transcripts of 1155 sequences were predicted as target genes of differentially expressed miRNAs by degradome sequencing. Gene Ontology (GO) enrichment showed that these target genes were involved in 24 biological processes, 15 cell components and nine molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that these target genes were significantly enriched in 19 pathways and the enriched KEGG pathways were associated with environmental adaptation, signal transduction and translation. Based on the functional prediction of miRNAs and target genes, our findings suggest that miRNAs have a potential regulatory role in cucumber fruit expansion by targeting their target genes, which provide important data for understanding the miRNA-mediated regulatory networks controlling fruit expansion in cucumber. Specific miRNAs could be selected for further functional research and molecular breeding in cucumber.

DELLA gene family plays a key role in regulating plant development and responding to stress. Currently, many DELLA family members have been identified in plants, however, information on DELLA genes in pumpkin (Cucurbita moschata) is scarce. In this study, physical and chemical properties, gene structure cis-regulatory elements and expression of CmoDELLA genes were examined in pumpkin. We found that seven CmoDELLA genes were identified in pumpkin, and they were unevenly classified into five chromosomes. CmoDELLA proteins were relatively unstable and their secondary structures were mainly made up α-helix and random coil. All seven CmoDELLA proteins contained typical DELLA domain and GRAS domain, however, motif numbers between CmoDELLA proteins were unevenly distributed, implying the complex evolution and functional diversification of CmoDELLA proteins. Cis-regulatory elements analysis revealed that CmoDELLA genes might play an essential role in regulating plant growth and development, and response to stress in pumpkin. Transcriptome data in the roots, stems, leaves and fruits demonstrated that CmoDELLA2, CmoDELLA3 and CmoDELLA7 were related to the stems development, CmoDELLA1, CmoDELLA4, CmoDELLA5 and CmoDELLA6 were associated with the fruits development. Furthermore, we found that CmoDELLA1 and CmoDELLA5 were up-regulated under NaCl stress. CmoDELLA1, CmoDELLA2, CmoDELLA3, CmoDELLA5, CmoDELLA6 and CmoDELLA7 were remarkably induced under waterlogging stress. While, all of the 7 CmoDELLA genes showed significantly induced expression under cold stress. The expression patterns under abiotic stress suggested that CmoDELLA genes might mediate the stress response of pumpkin to NaCl, waterlogging and cold, however, the functions of different CmoDELLA genes varied under different stress. Overall, our study provides valuable information for further research about the potential functions and regulatory networks of CmoDELLA genes in pumpkin.

The miR395 plays an indispensable role in biochemical processes by regulating their target genes. However, little is known about the roles of miR395 in cucumber fruit expansion and response to abiotic stresses. Here, 4 Csa-miR395s and 8 corresponding target genes were identified in the cucumber genome. Csa-miR395s were all located on the same chromosome (Chr 5). Csa-miR395a/b/c and Csa-miR395d were distributed in different branches without a closer genetic relationship. Massive cis-acting elements, including light, phytohormone, and stress response elements, were detected in the promoter regions of Csa-MIR395s, indicating that Csa-miR395s might be involved in complex regulatory networks to control cucumber growth and development and stress response. In addition, Csa-miR395a/b/c shared the same target genes, and Csa-miR395d had its specific target genes. Tissue-specific expression analysis showed that Csa-miR395a/b/c were all expressed in the leaf, root, ovary, and expanded fruit of cucumber and highly expressed in the expanded fruits compared to the ovary, while Csa2G215520 and Csa1G502860 (target genes of Csa-miR395a/b/c) presented a downregulated trend in the expanded fruit compared to the ovary. Meanwhile, the protein co-expression network revealed that these target genes had interactions in sulfur metabolism. These results suggested that Csa-miR395a/b/c targeting Csa2G215520 and Csa1G502860 might promote cucumber fruit expansion by affecting sulfur metabolism. Additionally, Quantitative Real-time PCR analysis validated that Csa-miR395s could be regulated by NaCl stress, and Csa-miR395a/b/c could respond to PEG stress, which further confirmed the reliability of cis-acting elements data. Taken together, our results could be helpful for further exploration of the functions of miR395s in cucumber fruit expansion and response to abiotic stresses.

The transcription factor bHLH gene family plays fundamental roles in plant development and mitigating diverse biotic and abiotic stresses. However, the information of bHLH genes in pumpkin (Cucurbita moschata) is still unknown. In this current study, 222 CmobHLH genes were identified and mapped onto different chromosomes through bioinformatics analysis in pumpkin. CmobHLH and AtbHLH proteins could be classified into 19 subfamilies according to the phylogenetic tree. CmobHLH proteins within the same subfamily had similar motif composition and gene structures. Gene ontology (GO), cis-regulatory elements (CREs) and protein–protein interaction analyses suggested the potential regulatory roles of CmobHLH genes during the plant development process and abiotic stresses response in pumpkin. Tissue expression patterns based on transcriptome data demonstrated that CmobHLH genes were involved in pumpkin development process, and they had unique functions in different tissues. The expression patterns of five selected CmobHLH genes after exposure to abiotic stresses showed that the CmobHLH genes played varied roles in the stress responses of pumpkin to NaCl, waterlogging, cold, ABA and drought. In brief, these findings offer important information for further functional research of CmobHLH genes and resistance breeding in pumpkin.