Comparative transcriptome analysis provides insights into molecular mechanisms for parthenocarpic fruit development in eggplant (Solanum melongena L.)

Abstract

Genetic control of parthenocarpy, a desirable trait in edible fruit with hard seeds, has been extensively studied. However, the molecular mechanism of parthenocarpic fruit development in eggplant (Solanum melongena L.) is still unclear. To provide insights into eggplant parthenocarpy, the transcriptomic profiles of a natural parthenocarpic (PP05) and two non-parthenocarpic (PnP05 and GnP05) eggplant lines were analyzed using RNA-sequencing (RNA-seq) technology. These sequences were assembled into 38925 unigenes, of which 22683 had an annotated function and 3419 were predicted as novel genes or from alternative splicing. 4864 and 1592 unigenes that were identified as DEGs between comparison groups PP05 vs PnP05 and PP05 vs GnP05, respectively. 506 common DEGs were found contained in both comparison groups, including 258 up-regulated and 248 down-regulated genes. Functional enrichment analyses identified many common or specific biological processes and gene set potentially associated with plant development. The most pronounced findings are that differentially regulated genes potentially-related with auxin signaling between parthenocarpic and non-parthenocarpic eggplants, e.g. calcium-binding protein PBP1 and transcription factor E2FB, which mediate the auxin distribution and auxin-dependent cell division, respectively, are up-regulated in the PP05; whereas homologs of GH3.1 and AUX/IAA, which are involved in inactivation of IAA and interference of auxin signaling, respectively, are down-regulated in PP05. Furthermore, gibberellin and cytokinin signaling genes and genes related to flower development were found differentially regulated between these eggplant lines. The present study provides comprehensive transcriptomic profiles of eggplants with or without parthenocarpic capacity. The information will deepen our understanding of the molecular mechanisms of eggplant parthenocarpy. The DEGs, especially these filtered from PP05 vs PnP05 + GnP05, will be valuable for further investigation of key genes involved in the parthenocarpic fruit development and genomics-assisted breeding.

Fig 1. Principal component analysis of four replicates of three eggplant PP05, PnP05 and GnP05.

Fig 2. Venn diagram depicting the numbers of DEGs found in these three eggplant lines.

Fig 3. Histograms showing GO functional analysis of the up-regulated (A) and down-regulated (B) common DEGs.

The x-axis indicates the names of the eight most abundant classes in each of the three main GO categories. The y-axis indicates the numbers of common DEGs.

Fig 4. Statistics of KEGG enrichment of the up-regulated (A) and down-regulated (B) common DEGs.

Histograms indicate 10 most abundant KEGG pathways.