RNA-seq-based transcriptome profiling of early fruit development in Chieh-qua and analysis of related transcription factors

Chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How.) fruit development starts post pollination. With the continuous expansion of the fruit, the soluble solid content of the fruit decreases. Because there are no reports on the early development of Chieh-qua fruit, this study compared fruit transcriptomes at 0-, 3-, and 7 day post pollination (dpp). 104,747 unigenes were assembled from clean reads and compared using six public databases for similarity searching. Compared with those of 0 dpp (C), there were differences in the expression of 12,982 and 6541 genes in the fruit tissue at 3 dpp and 7 dpp, respectively. Compared with 3 dpp (B), there were 14,314 differentially expressed genes in the fruit at 7 dpp (A). Based on the analysis of transcription factors, 213 nucleotides in the MYB superfamily were identified; among them, 94 unigenes of the MYB superfamily were differentially expressed at the three stages. In the pairwise comparison of differential expression, eight unigenes (Gene_id: TRINITY_DN32880_c1_g2, TRINITY_DN35142_c2_g2, TRINITY_DN32454_c11_g6, TRINITY_DN34105_c2_g7, TRINITY_DN32758_c3_g3, TRINITY_DN33604_c4_g10, TRINITY_DN34466_c3_g1, TRINITY_DN35924_c3_g2) were homologous to those of MYB59, MYB-GT3b, MYB18, MYB4, MYB108, MYB306, MYB340, and MYB-bHLH13. These unigenes differed significantly among the three stages. Furthermore, MYB59 and MYB18 exhibited higher expression at 7 dpp. MYB4, MYB-GT3b, MYB108, and MYB306 showed the highest expression levels in fruits at 3 dpp. In addition, MYB340 and MYB-bHLH13 showed higher expression levels during the unpollinated stage. MYB59, MYB-GT3b, MYB18, MYB4, MYB108, MYB306, MYB340, and MYB-bHLH13 may play crucial roles in Chieh-qua fruit development, defense, and blossoming. This study provides a basis for further investigation of MYB superfamily genes involved in early fruit expansion in chieh-qua.