Transcriptome profiling of Cucumis melo fruit development and ripening

Transcriptomics reveal useful resources for examining fruit development and variation in fruit size in Coccinia grandis

Introduction

The Cucurbitaceae family comprises many agronomically important members, that bear nutritious fruits and vegetables of great economic importance. Coccinia grandis, commonly known as Ivy gourd, belongs to this family and is widely consumed as a vegetable. Members of this family are known to display an impressive range of variation in fruit morphology. Although there have been studies on flower development in Ivy gourd, fruit development remains unexplored in this crop.

Methods

In this study, comparative transcriptomics of two Ivy gourd cultivars namely "Arka Neelachal Kunkhi" (larger fruit size) and "Arka Neelachal Sabuja" (smaller fruit size) differing in their average fruit size was performed. A de novo transcriptome assembly for Ivy gourd was developed by collecting fruits at different stages of development (5, 10, 15, and 20 days after anthesis i.e. DAA) from these two varieties. The transcriptome was analyzed to identify differentially expressed genes, transcription factors, and molecular markers.

Results

The transcriptome of Ivy gourd consisted of 155205 unigenes having an average contig size of 1472bp. Unigenes were annotated on publicly available databases to categorize them into different biological functions. Out of these, 7635 unigenes were classified into 38 transcription factor (TF) families, of which Trihelix TFs were most abundant. A total of 11,165 unigenes were found to be differentially expressed in both the varieties and the in silico expression results were validated through real-time PCR. Also, 98768 simple sequence repeats (SSRs) were identified in the transcriptome of Ivy gourd.

Discussion

This study has identified a number of genes, including transcription factors, that could play a crucial role in the determination of fruit shape and size in Ivy gourd. The presence of polymorphic SSRs indicated a possibility for marker-assisted selection for crop breeding in Ivy gourd. The information obtained can help select candidate genes that may be implicated in regulating fruit development and size in other fruit crops.

The maturation profile triggers differential expression of sugar metabolism genes in melon fruits

Melons are commercially important crops that requires specific quality attributes for successful commercialization, including accumulation of sugars, particularly sucrose. This trait can be influenced by various factors, such as the type of ripening. Cucumis melo L. is an ideal species for studying sugar metabolism because it has both climacteric and non-climacteric cultivars. Thus, this study aimed to examine the gene expression of sucrose metabolism candidates using RT-qPCR, in conjunction with postharvest physiological analyzes and high-performance liquid chromatography-based sugar quantification, in the melon cultivars 'Gaúcho' (climacteric) and 'Eldorado' (non-climacteric). The results showed that sucrose synthase 1 played a role in the synthesis and accumulation of sucrose in both cultivars, whereas sucrose synthase 2 was more highly expressed in 'Gaúcho', contributing to lower hexose content. Invertase inhibitor 1 was more highly expressed in 'Eldorado' and may be involved in sugar-induced maturation. Neutral α-galactosidase had distinct functions, playing a role in substrate synthesis for the growth of young 'Eldorado' fruits, whereas in mature 'Gaúcho' fruits it participated in the metabolism of raffinose family oligosaccharides for sucrose accumulation. The expression of trehalose-6-phosphate synthase genes indicated a greater involvement of these enzymes in the sugar regulation in 'Gaúcho' melons. These findings shed light on the intraspecific differences related to fruit quality attributes in different types of maturation and contribute to a deeper understanding of the underlying molecular mechanisms involved in the metabolism of sugars in melons, which can inform breeding programs aimed at improving fruit quality attributes in this crop.

Time-Series Transcriptome of Cucumis melo Reveals Extensive Transcriptomic Differences with Different Maturity

As the most important melon cultivar grown in the north-western provinces of China, Hami melon (Cucumis melo) produces large edible fruits that serve as an important dietary component in the world. In general, as a climacteric plant, melon harvested at 60% maturity results in a product with bad quality, while the highest-quality product can be guaranteed when harvesting at 90% maturity. In order to clarify the genetic basis of their distinct profiles of metabolite accumulation, we performed systematic transcriptome analyses between 60% and 90% maturity melons. A total of 36 samples were sequenced and over 1.7 billion reads were generated. Differentially expressed genes in 60% and 90% maturity melons were detected. Hundreds of these genes were functionally enriched in the sucrose and citric acid accumulation process of C. melo. We also detected a number of distinct splicing events between 60% and 90% maturity melons. Many genes associated with sucrose and citric acid accumulation displayed as differentially expressed or differentially spliced between different degrees of maturity of Hami melons, including CmCIN2, CmSPS2, CmBGAL3, and CmSPS2. These results demonstrate that the phenotype pattern differences between 60% and 90% maturity melons may be largely resulted from the significant transcriptome regulation.

Basic leucine zipper gene VvbZIP61 is expressed at a quantitative trait locus for high monoterpene content in grape berries

The widely appreciated muscat flavor of grapes and wine is mainly attributable to the monoterpenes that accumulate in ripe grape berries. To identify quantitative trait loci (QTL) for grape berry monoterpene content, an F1 mapping population was constructed by a cross between two grapevine genotypes, one with neutral aroma berries (cv. 'Beifeng') and the other with a pronounced muscat aroma (elite Vitis vinifera line '3-34'). A high-density genetic linkage map spanning 1563.7 cM was constructed using 3332 SNP markers that were assigned to 19 linkage groups. Monoterpenes were extracted from the berry of the F1 progeny, then identified and quantified by gas chromatography-mass spectrometry. Twelve stable QTLs associated with the amounts of 11 monoterpenes in berries were thus identified. In parallel, the levels of RNA in berries from 34 diverse cultivars were estimated by RNA sequencing and compared to the monoterpene content of the berries. The expression of five genes mapping to stable QTLs correlated well with the monoterpene content of berries. These genes, including the basic leucine zipper VvbZIP61 gene on chromosome 12, are therefore considered as potentially being involved in monoterpene metabolism. Overexpression of VvbZIP61 in Vitis amurensis callus through Agrobacterium-mediated transformation significantly increased the accumulation of several monoterpenes in the callus, including nerol, linalool, geranial, geraniol, β-myrcene, and D-limonene. It is hypothesized that VvbZIP61 expression acts to increase muscat flavor in grapes. These results advance our understanding of the genetic control of monoterpene biosynthesis in grapes and provide important information for the marker-assisted selection of aroma compounds in grape breeding.

Genome-Wide Identification and Expression Patterns of Cucumber Invertases and Their Inhibitor Genes

Invertases and their inhibitors play important roles in sucrose metabolism, growth and development, signal transduction, and biotic and abiotic stress tolerance in many plant species. However, in cucumber, both the gene members and functions of invertase and its inhibitor families remain largely unclear. In this study, in comparison with the orthologues of Citrullus lanatus (watermelon), Cucumis melo (melon), and Arabidopsis thaliana (Arabidopsis), 12 invertase genes and 12 invertase inhibitor genes were identified from the genome of Cucumis sativus (cucumber). Among them, the 12 invertase genes were classified as 4 cell wall invertases, 6 cytoplasmic invertases, and 2 vacuolar invertases. Most invertase genes were conserved in cucumber, melon, and watermelon, with several duplicate genes in melon and watermelon. Transcriptome analysis distinguished these genes into various expression patterns, which included genes CsaV3_2G025540 and CsaV3_2G007220, which were significantly expressed in different tissues, organs, and development stages, and genes CsaV3_7G034730 and CsaV3_5G005910, which might be involved in biotic and abiotic stress. Six genes were further validated in cucumber based on quantitative real-time PCR (qRT-PCR), and three of them showed consistent expression patterns as revealed in the transcriptome. These results provide important information for further studies on the physiological functions of cucumber invertases (CSINVs) and their inhibitors (CSINHs).

Phloem unloading in cultivated melon fruits follows an apoplasmic pathway during enlargement and ripening

Melon (Cucumis melo L.) has a long history of cultivation worldwide. During cultivation, domestication, and selection breeding, the sugar content of mature melon fruits has been significantly increased. Compared with unsweet melon and wild melon, rapid sucrose accumulation can occur in the middle and late stages of sweet melon fruit development. The phloem unloading pathway during the evolution and development of melon fruit has not been identified and analyzed. In this study, the phloem unloading pathway and the function of related sugar transporters in cultivated and wild melon fruits were analyzed by CFDA [5(6)-carbofluorescein diacetate] and esculin tracing, cytological pathway observation, qRT-PCR, and gene function analysis, etc. Results show that the phloem unloading pathway of wild melon fruit is largely symplastic, whereas the phloem unloading pathway of cultivated melon fruit shifts from symplastic to apoplasmic during development. According to a fruit grafting experiment, the fruit sink accumulates sugars independently. Correlation analysis showed that the expression amounts of several sucrose transporter genes were positively correlated with the sucrose content of melon fruit. Furthermore, CmSWEET10 was proved to be a sucrose transporter located on the plasma membrane of the phloem and highly expressed in the premature stage of sweet melon fruits, which means it may be involved in phloem apoplast unloading and sucrose accumulation in sweet melon fruits. Finally, we summarize a functional model of related enzymes and sugar transporters involved in the apoplast unloading of sweet melon fruits during enlargement and sucrose accumulation.

Combined Metabolome and Transcriptome Analysis Elucidates Sugar Accumulation in Wucai (Brassica campestris L.)

Wucai (Brassica campestris L.) is a leafy vegetable that originated in China, its soluble sugars accumulate significantly to improve taste quality during maturation, and it is widely accepted by consumers. In this study, we investigated the soluble sugar content at different developmental stages. Two periods including 34 days after planting (DAP) and 46 DAP, which represent the period prior to and after sugar accumulation, respectively, were selected for metabolomic and transcriptomic profiling. Differentially accumulated metabolites (DAMs) were mainly enriched in the pentose phosphate pathway, galactose metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, and fructose and mannose metabolism. By orthogonal projection to latent structures-discriminant s-plot (OPLS-DA S-plot) and MetaboAnalyst analyses, D-galactose and β-D-glucose were identified as the major components of sugar accumulation in wucai. Combined with the transcriptome, the pathway of sugar accumulation and the interact network between 26 DEGs and the two sugars were mapped. CWINV4, CEL1, BGLU16, and BraA03g023380.3C had positive correlations with the accumulation of sugar accumulation in wucai. The lower expression of BraA06g003260.3C, BraA08g002960.3C, BraA05g019040.3C, and BraA05g027230.3C promoted sugar accumulation during the ripening of wucai. These findings provide insights into the mechanisms underlying sugar accumulation during commodity maturity, providing a basis for the breeding of sugar-rich wucai cultivars.

New Insights into MdSPS4-Mediated Sucrose Accumulation under Different Nitrogen Levels Revealed by Physiological and Transcriptomic Analysis

Nitrogen nutrition participates in many physiological processes and understanding the physiological and molecular mechanisms of apple responses to nitrogen is very significant for improving apple quality. This study excavated crucial genes that regulates sugar metabolism in response to nitrogen in apples through physiology and transcriptome analysis, so as to lay a theoretical foundation for improving fruit quality. In this paper, the content of sugar and organic acid in apple fruit at different developmental periods under different nitrogen levels (0, 150, 300, and 600 kg·hm^-2) were determined. Then, the transcriptomic analysis was performed in 120 days after bloom (DAB) and 150 DAB. The results showed that the fructose and glucose content were the highest at 120 DAB under 600 kg·hm^-2 nitrogen level. Meanwhile, different nitrogen treatments decreased malate content in 30 and 60 DAB. RNA-seq analysis revealed a total of 4537 UniGenes were identified as differentially expressed genes (DEGs) under nitrogen treatments. Among these DEGs, 2362 (52.06%) were up-regulated and 2175 (47.94%) were down-regulated. The gene co-expression clusters revealed that most DEGs were significantly annotated in the photosynthesis, glycolysis/gluconeogenesis, pyruvate metabolism, carbon metabolism, carbon fixation in photosynthetic organisms and plant hormone signal transduction pathways. The key transcription factor genes (ERF, NAC, WRKY, and C2H2 genes) were differentially expressed in apple fruit. Sugar and acid metabolism-related genes (e.g., HXK1, SPS4, SS2, PPC16-2, and MDH2 genes) exhibited significantly up-regulated expression at 120 DAB, whereas they were down-regulated at 150 DAB. Furthermore, the MdSPS4 gene overexpression positively promoted sucrose accumulation in apple callus and fruit. In conclusion, the combinational analysis of transcriptome and the functional validation of the MdSPS4 gene provides new insights into apple responses to different nitrogen levels.

Multispecies transcriptomes reveal core fruit development genes

During angiosperm evolution there have been repeated transitions from an ancestral dry fruit to a derived fleshy fruit, often with dramatic ecological and economic consequences. Following the transition to fleshy fruits, domestication may also dramatically alter the fruit phenotype via artificial selection. Although the morphologies of these fruits are well documented, relatively less is known about the molecular basis of these developmental and evolutionary shifts. We generated RNA-seq libraries from pericarp tissue of desert tobacco and both cultivated and wild tomato species at common developmental time points and combined this with corresponding, publicly available data from Arabidopsis and melon. With this broadly sampled dataset consisting of dry/fleshy fruits and wild/domesticated species, we applied novel bioinformatic methods to investigate conserved and divergent patterns of gene expression during fruit development and evolution. A small set of 121 orthologous "core" fruit development genes show a common pattern of expression across all five species. These include key players in developmental patterning such as orthologs of KNOLLE, PERIANTHIA, and ARGONAUTE7. GO term enrichment suggests that these genes function in basic cell division processes, cell wall biosynthesis, and developmental patterning. We furthermore uncovered a number of "accessory" genes with conserved expression patterns within but not among fruit types, and whose functional enrichment highlights the conspicuous differences between these phenotypic classes. We observe striking conservation of gene expression patterns despite large evolutionary distances, and dramatic phenotypic shifts, suggesting a conserved function for a small subset of core fruit development genes.

Identification of Key Gene Networks Controlling Soluble Sugar and Organic Acid Metabolism During Oriental Melon Fruit Development by Integrated Analysis of Metabolic and Transcriptomic Analyses

Oriental melon (Cucumis melo var. acidulus) is one of the most economically important fruit crops worldwide. To elucidate the molecular basis related to soluble sugar and organic acid metabolism in the fruits of two oriental melon cultivars with different sweetness, we performed integrated metabolomic and transcriptomic analyses of the fruits of 'Tianbao' (A) with high sweetness and 'Xiaocuigua' (B) with low sweetness at different ripening stages. The high accumulation of sucrose, D-glucose, D-(+)-raffinose, and the relatively lower citric acid and malic acid might contribute to the sweet taste of A. By screening the differentially expressed genes (DEGs) and correlation analysis of the DEGs and differentially accumulated metabolites, we deduced that the B cultivar might promote the conversion of glucose and fructose into intermediate compounds for downstream processes such as glycolysis. The tricarboxylic acid (TCA) cycle might also be enhanced compared to A, thus resulting in the differential accumulation of soluble sugars and organic acids, ultimately causing the taste difference between the two oriental melon cultivars. Our finding provides important information for further exploring the metabolic mechanisms of soluble sugars and organic acids in oriental melon.

Comparative transcriptomic analysis of two Cucumis melo var. saccharinus germplasms differing in fruit physical and chemical characteristics

BACKGROUND

Hami melon (Cucumis melo var. saccharinus) is a popular fruit in China because of its excellent taste, which is largely determined by its physicochemical characteristics, including flesh texture, sugar content, aroma, and nutrient composition. However, the mechanisms by which these characteristics are regulated have not yet been determined. In this study, we monitored changes in the fruits of two germplasms that differed in physicochemical characteristics throughout the fruit development period.

RESULTS

Ripe fruit of the bred variety 'Guimi' had significantly higher soluble sugar contents than the fruit of the common variety 'Yaolong.' Additionally, differences in fruit shape and color between these two germplasms were observed during development. Comparative transcriptome analysis, conducted to identify regulators and pathways underlying the observed differences at corresponding stages of development, revealed a higher number of differentially expressed genes (DEGs) in Guimi than in Yaolong. Moreover, most DEGs detected during early fruit development in Guimi were associated with cell wall biogenesis. Temporal analysis of the identified DEGs revealed similar trends in the enrichment of downregulated genes in both germplasms, although there were differences in the enrichment trends of upregulated genes. Further analyses revealed trends in differential changes in multiple genes involved in cell wall biogenesis and sugar metabolism during fruit ripening.

CONCLUSIONS

We identified several genes associated with the ripening of Hami melons, which will provide novel insights into the molecular mechanisms underlying the development of fruit characteristics in these melons.

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

Extreme sugar accumulation in late fig ripening is accompanied by global changes in sugar metabolism and transporter gene expression

Female fig (Ficus carica L.) fruit are characterized by a major increase in volume and sugar content during the final week of development. A detailed developmental analysis of water and dry matter accumulation during these final days indicated a temporal separation between the increase in volume due to increasing water content and a subsequent sharp increase in sugar content during a few days. The results present fig as an extreme example of sugar import and accumulation, with calculated import rates that are one order of magnitude higher than those of other sugar-accumulating sweet fruit species. To shed light on the metabolic changes occurring during this period, we followed the expression pattern of 80 genes encoding sugar metabolism enzymes and sugar transporter proteins identified in fig fruit. A parallel comparison with male fig fruits, which do not accumulate sugar during ripening, highlighted the genes specifically related to sugar accumulation. Tissue-specific analysis indicated that the expression of genes involved in sugar metabolism and transport undergoes a global transition.

Flavonoid biosynthesis in four Dendrobium species based on transcriptome sequencing and metabolite analysis

BACKGROUND

Dendrobium is a genus of plants used as traditional Chinese herbal medicines, with high economic and medicinal value.

METHODS AND RESULTS

To reveal the mechanism of flavonoid biosynthesis in Dendrobium, the metabolites and transcriptomes of four Dendrobium species (D. chrysotoxum, D. nobile, D. fimbriatum, and D. denneanum) were analyzed comprehensively. Ultra-high-performance liquid chromatography-tandem mass spectrometry analysis revealed ten flavonoid compounds in Dendrobium. In total, 100,096 unigenes were obtained from the transcript database of the four Dendrobium species. Among the identified differentially expressed genes, 51 were associated with flavonoid biosynthesis, and 670 differentially expressed transcription factors were predicted, including 194 MYB, 87 bHLH, and 100 WRKY family transcription factors, respectively. Transcriptome analysis showed that the expression levels of structural genes such as chalcone synthase (CHS), cinnamate-4-hydroxylase (C4H), and flavonoid 3'-hydroxylase (F3'H) were lower in D. chrysotoxum, D. nobile, and D. fimbriatum than those in D. denneanum, which may be the main reason for the low flavonoid contents in D. chrysotoxum, D. nobile, and D. fimbriatum.

CONCLUSIONS

The expression level of structural genes corresponded to the accumulation level of flavonols in the different Dendrobium species. The results deepen the understanding of the molecular mechanism of flavonoid biosynthesis in Dendrobium and provide novel insights into the synthesis and accumulation of flavonoids in Dendrobium.

Alternative splicing during fruit development among fleshy fruits

BACKGROUND

Alternative splicing (AS) is an important mechanism of posttranscriptional modification and dynamically regulates multiple physiological processes in plants, including fruit ripening. However, little is known about alternative splicing during fruit development in fleshy fruits.

RESULTS

We studied the alternative splicing at the immature and ripe stages during fruit development in cucumber, melon, papaya and peach. We found that 14.96-17.48% of multiexon genes exhibited alternative splicing. Intron retention was not always the most frequent event, indicating that the alternative splicing pattern during different developmental process differs. Alternative splicing was significantly more prevalent at the ripe stage than at the immature stage in cucumber and melon, while the opposite trend was shown in papaya and peach, implying that developmental stages adopt different alternative splicing strategies for their specific functions. Some genes involved in fruit ripening underwent stage-specific alternative splicing, indicating that alternative splicing regulates fruits ripening. Conserved alternative splicing events did not appear to be stage-specific. Clustering fruit developmental stages across the four species based on alternative splicing profiles resulted in species-specific clustering, suggesting that diversification of alternative splicing contributes to lineage-specific evolution in fleshy fruits.

CONCLUSIONS

We obtained high quality transcriptomes and alternative splicing events during fruit development across the four species. Dynamics and nonconserved alternative splicing were discovered. The candidate stage-specific AS genes involved in fruit ripening will provide valuable insight into the roles of alternative splicing during the developmental processes of fleshy fruits.

Transition of aromatic volatile and transcriptome profiles during melon fruit ripening

Melon (Cucumis melo L.) is an important diploid crop with a wide variety of flavors due to its distinct aromatic volatile organic compounds (VOC). To understand the development of VOC profiles during fruit development, we performed metabolomic and transcriptomic analysis of two cantaloupe varieties over the course of fruit development. A total of 130 metabolites were detected in fruit samples, and 449014207 reads were mapped to the melon genome. A total of 4469 differentially expressed genes in fruits were identified and used to visualize the transition of VOC and transcriptomic profiles during the fruit development. A shift of VOC profiles in both varieties was observed from early-fruit profiles enriched in C5-C8 lipid-derived VOCs to late-fruit profiles abundant in C9 lipid-derived VOCs, apocarotenoids, and esters. The shift coincided with the expression of specific isoforms of lipid and carotenoid metabolizing enzymes as well as transcription factors involved in fruit ripening, metabolite regulation, and hormone signaling.

Transcriptomic analysis of a near-isogenic line of melon with high fruit flesh firmness during ripening

BACKGROUND

A near-isogenic line (NIL) of melon (SC10-2) with introgression in linkage group X was studied from harvest (at firm-ripe stage of maturity) until day 18 of postharvest storage at 20.5 °C together with its parental control ('Piel de Sapo', PS).

RESULTS

SC10-2 showed higher flesh firmness and whole fruit hardness but lower juiciness than its parental. SC10-2 showed a decrease in respiration rate accompanied by a decrease in ethylene production during ripening, both of which fell to a greater extent than in PS. The introgression affected 11 volatile organic compounds (VOCs), the levels of which during ripening were generally higher in SC10-2 than in PS. Transcriptomic analysis from RNA-Seq revealed differentially expressed genes (DEGs) associated with the effects studied. For example, 909 DEGs were exclusive to the introgression, and only 23 DEGs were exclusive to postharvest ripening time. Major functions of the DEGs associated with introgression or ripening time were identified by cluster analysis. About 37 genes directly and/or indirectly affected the delay in ripening of SC10-2 compared with PS in general and, more particularly, the physiological and quality traits measured and, probably, the differential non-climacteric response. Of the former genes, we studied in more detail at least five that mapped in the introgression in linkage group (LG) X, and 32 outside it.

CONCLUSION

There is an apparent control of textural changes, VOCs and fruit ripening by an expression quantitative trait locus located in LG X together with a direct control on them due to genes presented in the introgression (CmTrpD, CmNADH1, CmTCP15, CmGDSL esterase/lipase, and CmHK4-like) and CmNAC18. © 2020 Society of Chemical Industry.

Comprehensive analysis of microRNAs in the mantle central and mantle edge provide insights into shell formation in pearl oyster Pinctada fucata martensii

MicroRNAs (miRNAs) are a class of non-coding RNA molecules with post-transcriptional regulatory activity in various biological processes. Pearl oyster Pinctada fucata martensii is one of the main species cultured for marine pearl production in China and Japan. In this study, we constructed two small RNA libraries of mantle central (MC) and mantle edge (ME) from P. f. martensii and obtained 24,175,537 and 21,593,898 clean reads, respectively. A total of 258 miRNAs of P. f. martensii (Pm-miRNA) were identified, and 93 differentially expressed miRNAs (DEMs) including 49 known Pm-miRNAs and 44 novel Pm-miRNAs were obtained from the MC and ME. The target transcripts of these DEMs were obviously enriched in neuroactive ligand-receptor interaction pathway, and others. After over-expression of Pm-miR-124 and Pm-miR-9a-5p in the MC by mimic injection into the muscle of P. f. martensii, nacre exhibited a disorderly growth as detected by scanning electron microscopy. Pm-nicotinic acetylcholine receptor alpha subunit, Pm-neuropeptide Y and Pm-chitin synthase were investigated as the targets of Pm-miR-124; and Pm-tumor necrosis factor receptor associated factor 2 and Pm-chitin synthase were investigated as the targets of Pm-miR-9a-5p. These predicted target transcripts were down-regulated after the over-expression of Pm-miR-124 and Pm-miR-9a-5p in MC. This study comprehensively analyzed the miRNAs in mantle tissues to enhance our understanding of the regulatory mechanism underlying shell formation.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

Metabolic and transcriptomic analysis of two Cucurbita moschata germplasms throughout fruit development

BACKGROUND

Pumpkins (Cucurbita moschata; Cucurbitaceae) are valued for their fruits and seeds and are rich in nutrients. Carotenoids and sugar contents, as main feature of pumpkin pulp, are used to determine the fruit quality.

RESULTS

Two pumpkin germplasms, CMO-X and CMO-E, were analyzed regarding the essential quality traits such as dry weight, soluble solids, organic acids, carotenoids and sugar contents. For the comparison of fruit development in these two germplasms, fruit transcriptome was analyzed at 5 different developmental stages from 0 d to 40 d in a time course manner. Putative pathways for carotenoids biosynthesis and sucrose metabolism were developed in C. moschata fruit and homologs were identified for each key gene involved in the pathways. Gene expression data was found consistent with the accumulation of metabolites across developmental stages and also between two germplasms. PSY, PDS, ZEP, CRTISO and SUS, SPS, HK, FK were found highly correlated with the accumulation of carotenoids and sucrose metabolites, respectively, at different growth stages of C. moschata as shown by whole transcriptomic analysis. The results of qRT-PCR analysis further confirmed the association of these genes.

CONCLUSION

Developmental regulation of the genes associated with the metabolite accumulation can be considered as an important factor for the determination of C. moschata fruit quality. This research will facilitate the investigation of metabolic profiles in other cultivars.

Transcriptome profiling of non-climacteric 'yellow' melon during ripening: insights on sugar metabolism

BACKGROUND

The non-climacteric 'Yellow' melon (Cucumis melo, inodorus group) is an economically important crop and its quality is mainly determined by the sugar content. Thus, knowledge of sugar metabolism and its related pathways can contribute to the development of new field management and post-harvest practices, making it possible to deliver better quality fruits to consumers.

RESULTS

The RNA-seq associated with RT-qPCR analyses of four maturation stages were performed to identify important enzymes and pathways that are involved in the ripening profile of non-climacteric 'Yellow' melon fruit focusing on sugar metabolism. We identified 895 genes 10 days after pollination (DAP)-biased and 909 genes 40 DAP-biased. The KEGG pathway enrichment analysis of these differentially expressed (DE) genes revealed that 'hormone signal transduction', 'carbon metabolism', 'sucrose metabolism', 'protein processing in endoplasmic reticulum' and 'spliceosome' were the most differentially regulated processes occurring during melon development. In the sucrose metabolism, five DE genes are up-regulated and 12 are down-regulated during fruit ripening.

CONCLUSIONS

The results demonstrated important enzymes in the sugar pathway that are responsible for the sucrose content and maturation profile in non-climacteric 'Yellow' melon. New DE genes were first detected for melon in this study such as invertase inhibitor LIKE 3 (CmINH3), trehalose phosphate phosphatase (CmTPP1) and trehalose phosphate synthases (CmTPS5, CmTPS7, CmTPS9). Furthermore, the results of the protein-protein network interaction demonstrated general characteristics of the transcriptome of young and full-ripe melon and provide new perspectives for the understanding of ripening.

Comparative Transcriptomic Analysis of Two Bottle Gourd Accessions Differing in Fruit Size

The bottle gourd (Lagenaria siceraria) is an important horticultural and medicinal crop with high nutritional value. This study aimed at examining the molecular regulation of fruit size in bottle gourd. We performed transcriptome sequencing of two bottle gourd cultivars differing in their fruit size. The average fruit length and weight of the cultivar Hang (39.48 cm/624.4 g) were higher than those of the cultivar USA (10.34 cm/152.8 g) at maturity. Transcriptome sequencing and assembly resulted in 89,347 unigenes. A total of 1250 differentially expressed genes (DEG) were found between the two cultivars, including 422 upregulated genes and 828 downregulated genes in Hang as compared to USA. Genes related to cell wall metabolism, phytohormones, cell cycle, and cell division showed significant differential expression between the two cultivars. DEGs encoding transcription factors (TF) from nine TF families were also identified. The ethylene response factor family was the most enriched among these families. Our study provides a basis for further investigations of the molecular regulation of fruit size in bottle gourd.

Transcriptome Analysis of Elm (Ulmus pumila) Fruit to Identify Phytonutrients Associated Genes and Pathways

Plant fruit is an important source of natural active phytonutrients that are profitable for human health. Elm (Ulmus pumila) fruit is considered as natural plant food in China that is rich in nutrients. In the present study, high-throughput RNA sequencing was performed in U. pumila edible fruits and leaves and 11,386 unigenes were filtered as dysregulated genes in fruit samples, including 5231 up- and 6155 downregulated genes. Hundreds of pathways were predicted to participate in seed development and phytonutrient biosynthesis in U. pumila by GO, MapMan, and KEGG enrichment analysis, including “seed maturation”, “glycine, serine, and threonine metabolism” and “phenylpropanoid biosynthesis”. ABA-mediated glucose response-related ethylene-activated signaling pathway (e.g., ABI4) were supposed to associate with elm fruit development; unsaturated fatty acids pathway (e.g., ACX2 and SAD) were predicted to participate in determination of fatty acid composition in elm fruit; flavonoid and coumarins biosynthesis (e.g., CYP98A3 and CCoAOMT1) were demonstrated to correlate with the bioactivity of elm fruits in human cancer and inflammation resistance. To provide more information about fruit developmental status, the qRT-PCR analysis for key genes of “phenylpropanoid biosynthesis” and “alpha-Linolenic acid metabolism” were conducted in samples of young fruits, ripe fruit, old fruit, and leaves. Two biosynthetic pathways for unsaturated fatty acid and Jasmonic acid (JA) were deduced to be involved in fruit development in U. pumila and the phenylpropanoid glycoside, syringin, was speculated to accumulate in the early development stages of elm fruit. Our transcriptome data supports molecular clues for seed development and biologically active substances in elm fruits.

Genome-wide identification and characterization of long non-coding RNAs involved in fruit ripening and the climacteric in Cucumis melo

BACKGROUND

Cucumis melo is a suitable study material for investigation of fruit ripening owing to its climacteric nature. Long non-coding RNAs have been linked to many important biological processes, such as fruit ripening, flowering time regulation, and abiotic stress responses in plants. However, knowledge of the regulatory roles of lncRNAs underlying the ripening process in C. melo are largely unknown. In this study the complete transcriptome of Cucumis melo L. cv. Hetao fruit at four developmental stages was sequenced and analyzed. The potential role of lncRNAs was predicted based on the function of differentially expressed target genes and correlated genes.

RESULTS

In total, 3857 lncRNAs were assembled and annotated, of which 1601 were differentially expressed between developmental stages. The target genes of these lncRNAs and the regulatory relationship (cis- or trans-acting) were predicted. The target genes were enriched with GO terms for biological process, such as response to auxin stimulus and hormone biosynthetic process. Enriched KEGG pathways included plant hormone signal transduction and carotenoid biosynthesis. Co-expression network construction showed that LNC_002345 and LNC_000154, which were highly expressed, might co-regulate with mutiple genes associated with auxin signal transduction and acted in the same pathways. We identified lncRNAs (LNC_000987, LNC_000693, LNC_001323, LNC_003610, LNC_001263 and LNC_003380) that were correlated with fruit ripening and the climacteric, and may participate in the regulation of ethylene biosynthesis and metabolism and the ABA signaling pathway. A number of crucial transcription factors, such as ERFs, WRKY70, NAC56, and NAC72, may also play important roles in the regulation of fruit ripening in C. melo.

CONCLUSIONS

Our results predict the regulatory functions of the lncRNAs during melon fruit development and ripening, and 142 highly expressed lncRNAs (average FPKM > 100) were identified. These lncRNAs participate in the regulation of auxin signal transduction, ethylene, sucrose biosynthesis and metabolism, the ABA signaling pathway, and transcription factors, thus regulating fruit development and ripening.

De novo assembly of Persea americana cv. 'Hass' transcriptome during fruit development

BACKGROUND

Avocado (Persea americana Mill.) is a basal angiosperm from the Lauraceae family. This species has a diploid genome with an approximated size of ~ 920 Mbp and produces a climacteric, fleshy and oily fruit. The flowering and fruit set are particularly prolonged processes, lasting between one to three months, generating important differences in physiological ages of the fruit within the same tree. So far there is no detailed genomic information regarding this species, being the cultivar 'Hass' especially important for avocado growers worldwide. With the aim to explore the fruit avocado transcriptome and to identify candidate biomarkers to monitore fruit development, we carried out an RNA-Seq approach during 4 stages of 'Hass' fruit development: 150 days after fruit set (DAFS), 240 DAFS, 300 DAFS (harvest) and 390 DAFS (late-harvest).

RESULTS

The 'Hass' de novo transcriptome contains 62,203 contigs (x̅=988 bp, N50 = 1050 bp). We found approximately an 85 and 99% of complete ultra-conserved genes in eukaryote and plantae database using BUSCO (Benchmarking Universal Single-Copy Orthologs) and CEGMA (Core Eukaryotic Gene Mapping Approach), respectively. Annotation was performed with BLASTx, resulting in a 58% of annotated contigs (90% of differentially expressed genes were annotated). Differentially expressed genes analysis (DEG; with False Discovery Rate ≤ 0.01) found 8672 genes considering all developmental stages. From this analysis, genes were clustered according to their expression pattern and 1209 genes show correlation with the four developmental stages.

CONCLUSIONS

Candidate genes are proposed as possible biomarkers for monitoring the development of the 'Hass' avocado fruit associated with lipid metabolism, ethylene signaling pathway, auxin signaling pathway, and components of the cell wall.

The kinome of pineapple: catalog and insights into functions in crassulacean acid metabolism plants

BACKGROUND

Crassulacean acid metabolism (CAM) plants use water 20-80% more efficiently by shifting stomata opening and primary CO₂ uptake and fixation to the nighttime. Protein kinases (PKs) play pivotal roles in this biological process. However, few PKs have been functionally analyzed precisely due to their abundance and potential functional redundancy (caused by numerous gene duplications).

RESULTS

In this study, we systematically identified a total of 758 predicted PK genes in the genome of a CAM plant, pineapple (Ananas comosus). The pineapple kinome was classified into 20 groups and 116 families based on the kinase domain sequences. The RLK was the largest group, containing 480 members, and over half of them were predicted to locate at the plasma membrane. Both segmental and tandem duplications make important contributions to the expansion of pineapple kinome based on the synteny analysis. Ka/Ks ratios showed all of the duplication events were under purifying selection. The global expression analysis revealed that pineapple PKs exhibit different tissue-specific and diurnal expression patterns. Forty PK genes in a cluster performed higher expression levels in green leaf tip than in white leaf base, and fourteen of them had strong differential expression patterns between the photosynthetic green leaf tip and the non-photosynthetic white leaf base tissues.

CONCLUSIONS

Our findings provide insights into the evolution and biological function of pineapple PKs and a foundation for further functional analysis of PKs in CAM plants. The gene duplication, expression, and coexpression analysis helped us to rapidly identify the key candidates in pineapple kinome, which may play roles in the carbon fixation process in pineapple and help engineering CAM pathway into C3 crops for improved drought tolerance.

Melonet-DB, a Grand RNA-Seq Gene Expression Atlas in Melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important Cucurbitaceae crop produced worldwide, exhibiting wide genetic variations and comprising both climacteric and non-climacteric fruit types. The muskmelon cultivar "'Earl's favorite Harukei-3 (Harukei-3)"' known for its sweetness and rich aroma is used for breeding of high-grade muskmelon in Japan. We conducted RNA sequencing (RNA-seq) transcriptome studies in 30 different tissues of the 'Harukei-3' melon. These included root, stems, leaves, flowers, regenerating callus and ovaries, in addition to the flesh and peel sampled at seven stages of fruit development. The expression patterns of 20,752 genes were determined with fragments per kilobase of transcript per million fragments sequenced (FPKM) >1 in at least one tissue. Principal component analysis distinguished 30 melon tissues based on the global gene expression profile and, further, the weighted gene correlation network analysis classified melon genes into 45 distinct coexpression groups. Some coexpression groups exhibited tissue-specific gene expression. Furthermore, we developed and published web application tools designated "'Gene expression map viewer"' and "'Coexpression viewer"' on our website Melonet-DB (http://melonet-db.agbi.tsukuba.ac.jp/) to promote functional genomics research in melon. By using both tools, we analyzed melon homologs of tomato fruit ripening regulators such as E8, RIPENING-INHIBITOR (RIN) and NON-RIPENING (NOR). The "'Coexpression viewer"' clearly distinguished fruit ripening-associated melon RIN/NOR/CNR homologs from those expressed in other tissues. In addition, several other MADS-box, NAM/ATAF/CUC (NAC) and homeobox transcription factor genes were identified as fruit ripening-associated genes. Our tools provide useful information for research not only on melon but also on other fleshy fruit plants.

Genome-wide identification of Hami melon miRNAs with putative roles during fruit development

MicroRNAs represent a family of small endogenous, non-coding RNAs that play critical regulatory roles in plant growth, development, and environmental stress responses. Hami melon is famous for its attractive flavor and excellent nutritional value, however, the mechanisms underlying the fruit development and ripening remains largely unknown. Here, we performed small RNA sequencing to investigate the roles of miRNAs during Hami melon fruit development. Two batches of flesh samples were collected at four fruit development stages. Small RNA sequencing yielded a total of 54,553,424 raw reads from eight libraries. 113 conserved miRNAs belonging to 30 miRNA families and nine novel miRNAs comprising nine miRNA families were identified. The expression of 42 conserved miRNAs and three Hami melon-specific miRNAs significantly changed during fruit development. Furthermore, 484 and 124 melon genes were predicted as putative targets of 29 conserved and nine Hami melon-specific miRNA families, respectively. GO enrichment analysis were performed on target genes, "transcription, DNA-dependent", "rRNA processing", "oxidation reduction", "signal transduction", "regulation of transcription, DNA-dependent", and "metabolic process" were the over-represented biological process terms. Cleavage sites of six target genes were validated using 5' RACE. Our results present a comprehensive set of identification and characterization of Hami melon fruit miRNAs and their potential targets, which provide valuable basis towards understanding the regulatory mechanisms in programmed process of normal Hami fruit development and ripening. Specific miRNAs could be selected for further research and applications in breeding practices.

Comparative transcriptome analysis of two contrasting watermelon genotypes during fruit development and ripening

BACKGROUND

Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai] is an economically important crop with an attractive ripe fruit that has colorful flesh. Fruit ripening is a complex, genetically programmed process.

RESULTS

In this study, a comparative transcriptome analysis was performed to identify the regulators and pathways that are involved in the fruit ripening of pale-yellow-flesh cultivated watermelon (COS) and red-flesh cultivated watermelon (LSW177). We first identified 797 novel genes to extend the available reference gene set. Second, 3958 genes in COS and 3503 genes in LSW177 showed at least two-fold variation in expression, and a large number of these differentially expressed genes (DEGs) during fruit ripening were related to carotenoid biosynthesis, plant hormone pathways, and sugar and cell wall metabolism. Third, we noted a correlation between ripening-associated transcripts and metabolites and the key function of these metabolic pathways during fruit ripening.

CONCLUSION

The results revealed several ripening-associated actions and provide novel insights into the molecular mechanisms underlying the regulation of watermelon fruit ripening.