De novo transcriptome sequencing of Momordica cochinchinensis to identify genes involved in the carotenoid biosynthesis

Comparative transcriptome analysis reveals an insight into the candidate genes involved in anthocyanin and scent volatiles biosynthesis in colour changing flowers of Combretum indicum

Combretum indicum is a widely cultivated ornamental species displaying the distinct phenomenon of floral colour change. Flowers display a gradual colour change from white to red, attributed to increased cyanidin 3-O glucoside in petal tissues. The differently coloured flowers also emanate a complex blend of VOCs with trans-linalool oxide (furanoid) as the major compound in the emission profile. To understand molecular mechanisms regulating floral colour shifts and scent biosynthesis, we performed Illumina transcriptome sequencing, including de novo assembly and functional annotation, for the two stages of floral maturation (white and red). Homology analysis with functional classification identified 84 and 42 candidate genes associated with pigment and scent biosynthesis, respectively. Genes encoding transcription factors, such as MYB, ERF, WD40, WRKY, NAC, bHLH and bZIP, that play critical roles in regulating specialized metabolism were also identified in the transcriptome data. Differences in expression of genes were consistent with accumulation patterns of anthocyanins in the two different flower colours. A clear upregulation of flavonoid biosynthesis genes in red flower tissue is associated with increased pigment content. RT-qPCR-based expression analyses gave results consistent with the RNA-Seq data, suggesting the sequencing data are consistent and reliable. This study presents the first report of genetic information for C. indicum. Gene sequences generated from RNA-Seq, along with candidate genes identified by pathway mapping and their expression profiles, provide a valuable resource for subsequent studies towards molecular understanding of specialized metabolism in C. indicum flowers.

Transcriptome sequencing and DEG analysis in different developmental stages of floral buds induced by potassium chlorate in Dimocarpus longan

Potassium chlorate can promote off-season flowering in longan, but the molecular mechanisms are poorly understood. In this study, four-year-old 'Shixia' longan trees were injected in the trunk with potassium chlorate, and terminal buds were sampled and analyzed using transcriptomics and bioinformatics tools. To generate a reference longan transcriptome, we obtained 207,734 paired-end reads covering a total of 58,514,149 bp, which we assembled into 114,445 unigenes. Using this resource, we identified 3,265 differentially expressed genes (DEGs) that were regulated in longan terminal buds in response to potassium chlorate treatment for 2, 6 or 30 days, including 179 transcription factor genes. By reference to the Arabidopsis literature, we then defined 38 longan genes involved in flowering, from which we constructed the longan flowering pathway. According to RNA-seq data, at least 24 of these genes, which participate in multiple signaling pathways, are involved in potassium chlorate-stimulated floral induction, and the differential regulation in terminal buds of ten floral pathway genes (GI, CO, GID1, GA4, GA5, FLC, AP1, LFY, FT and SOC1) was confirmed by qRT-PCR. These data will contribute to an improved understanding of the functions of key genes involved in longan floral induction by potassium chlorate.

Comparative transcriptome analysis of Lupinus polyphyllus Lindl. provides a rich molecular resource for research on coloration mechanism

Lupinus polyphyllus is rich in color, making it a well-known horticultural ornamental plant. However, little is known about the genes related to anthocyanin and carotenoid biosynthesis in L. polyphyllus. In this study, transcriptome sequencing was performed on eight different colors of L. polyphyllus. A total of 1.13 billion clean reads were obtained and assembled into 89,124 unigenes, which were then aligned with six databases, resulting in the identification of 54,823 annotated unigenes. Among these unigenes, 76 and 101 were involved in the biosynthetic pathway of carotenoids and anthocyanins, respectively. In addition, 505 transcription factors were revealed, which belonged to the MYB, R2R3-MYB, NAC, bHLH, and WD40 families. A total of 6,700 differentially expressed genes (DEGs) were obtained by comparative transcriptome analysis. Among them, 17 candidate unigenes (four carotenoid genes, seven anthocyanin genes, and six TFs) were specifically up-regulated for one or more colors of L. polyphyllus. Eight representative candidate unigenes were analyzed by qRT-PCR. The findings enrich the transcriptome database of lupine, and provide a rich molecular resource for research on the coloration mechanism of L. polyphyllus.

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.

Identification of Key Genes Controlling Carotenoid Metabolism during Apricot Fruit Development by Integrating Metabolic Phenotypes and Gene Expression Profiles

To explore the metabolic basis of carotenoid accumulation in different developmental periods of apricot fruits, targeted metabonomic and transcriptomic analyses were conducted in four developmental periods (S1-S4) in two cultivars (Prunus armeniaca cv. "Kuchebaixing" with white flesh and P. armeniaca cv. "Shushangganxing" with orange flesh) with different carotenoid contents. 14 types of carotenes and 27 types of carotene lipids were identified in apricot flesh in different developmental periods. In S3 and S4, the carotenoid contents of the two cultivars were significantly different, and β-carotene and (E/Z)-phytoene were the key metabolites that caused the difference in the total carotenoid content between the examined cultivars. Twenty-five structural genes (including genes in the methylerythritol 4-phosphate and carotenoid biosynthesis pathways) related to carotenoid biosynthesis were identified among the differentially expressed genes in different developmental periods of the two cultivars, and a carotenoid metabolic pathway map of apricot fruits was drawn according to the KEGG pathway map. The combined analysis of carotenoid metabolism data and transcriptome data showed that PSY, NCED1, and CCD4 were the key genes leading to the great differences in the total carotenoid content. The results provide a new approach to study the synthesis and accumulation of carotenoids in apricot fruits.

Identification and characterization of apocarotenoid modifiers and carotenogenic enzymes for biosynthesis of crocins in Buddleja davidii flowers

Crocetin biosynthesis in Buddleja davidii flowers proceeds through a zeaxanthin cleavage pathway catalyzed by two carotenoid cleavage dioxygenases (BdCCD4.1 and BdCCD4.3), followed by oxidation and glucosylation reactions that lead to the production of crocins. We isolated and analyzed the expression of 12 genes from the carotenoid pathway in B. davidii flowers and identified four candidate genes involved in the biosynthesis of crocins (BdALDH, BdUGT74BC1, BdUGT74BC2, and BdUGT94AA3). In addition, we characterized the profile of crocins and their carotenoid precursors, following their accumulation during flower development. Overall, seven different crocins, crocetin, and picrocrocin were identified in this study. The accumulation of these apocarotenoids parallels tissue development, reaching the highest concentration when the flower is fully open. Notably, the pathway was regulated mainly at the transcript level, with expression patterns of a large group of carotenoid precursor and apocarotenoid genes (BdPSY2, BdPDS2, BdZDS, BdLCY2, BdBCH, BdALDH, and BdUGT Genes) mimicking the accumulation of crocins. Finally, we used comparative correlation network analysis to study how the synthesis of these valuable apocarotenoids diverges among B. davidii, Gardenia jasminoides, and Crocus sativus, highlighting distinctive differences which could be the basis of the differential accumulation of crocins in the three species.

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.

Transcriptomic analysis of selenium accumulation in Puccinellia distans (Jacq.) Parl., a boron hyperaccumulator

Selenium (Se) is present in a wide variety of natural and man-made materials on Earth. Plants are able to tolerate and (hyper)accumulate Se to different extents. In fact, some species can tolerate and accumulate multiple elements. Puccinellia distans (P. distans), weeping alkali grass, is known to hyperaccumulate extreme concentrations of boron and tolerate high levels of salinity, therefore, we investigated the Se accumulation and tolerance capacities of this species. In addition, P. distans' Se tolerance mechanism was studied using a transcriptomic approach. The results of this study indicated that, when grown in a hydroponic system containing 80 or 120 μM Se, P. distans shoots accumulated from 1500 to 2500-fold more Se than plants grown without the element. Thus, P. distans was discovered to be a novel Se accumulator plant. RNA sequencing results and biochemical analyses helped to shed light on the Se tolerance and accumulation mechanism of P. distans. Here, we suggest that upregulation of Se assimilation and stress response genes may be due to induction of jasmonic acid signaling. In addition, we propose that the cell wall may play an important role in restriction of Se movement to the cytoplasm. Also, we hypothesize that Se accumulates in cells by sequestration of selenate in the vacuole.

Berry transcriptome: insights into a novel resource to understand development dependent secondary metabolism in Withania somnifera (Ashwagandha)

Withania somnifera (Ashwagandha) is considered as Rasayana in Indian systems of medicine. This study reports a novel transcriptome of W. somnifera berries, with high depth, quality and coverage. Assembled and annotated transcripts for nearly all genes related with the withanolide biosynthetic pathway were obtained. Tissue-wide gene expression analysis reflected almost similar definitions for the terpenoid pathway in leaf, root and berry tissues with relatively higher abundance of transcripts linked to steroid, phenylpropanoid metabolism as well as flavonoid metabolism in berries. The metabolome map generated from the data embodied transcripts from 143 metabolic pathways connected together and mediated collectively by about 1792 unique enzyme functions specific to berry, leaf and root tissues, respectively. Transcripts specific to cytochrome p450 (CYP450), methyltransferases and glycosyltransferases were distinctively located in a tissue specific manner and exhibited a complex network. Significant distribution of transcription factor genes such as MYB, early light inducible protein (ELI), minichromosome maintenance1, agamous, deficiens and serum response factor (MADS) and WRKY etc. was observed, as the major transcriptional regulators of secondary metabolism. Validation of the assembly was ascertained by cloning WsELI, which has a light dependent regulatory role in development. Quantitative expression of WsELI was observed to be considerably modulated upon exposure to abiotic stress and elicitors. Co-relation of over-expression of WsELI, may provide new aspects for the functional role of ELI proteins in plants linked to secondary metabolism. The study offers the first comprehensive and comparative evaluation of W. somnifera transcriptome data between the three tissues and across other members of Solanaceae (Nicotiana, Solanum and Capsicum) with respect to major pathways and their metabolome regulation.

Analysis of Centranthera grandiflora Benth Transcriptome Explores Genes of Catalpol, Acteoside and Azafrin Biosynthesis

Cardiovascular diseases (CVDs) are a major cause of health loss in the world. Prevention and treatment of this disease by traditional Chinese medicine is a promising method. Centranthera grandiflora Benth is a high-value medicinal herb in the prevention and treatment of CVDs; its main medicinal components include iridoid glycosides, phenylethanoid glycosides, and azafrin in roots. However, biosynthetic pathways of these components and their regulatory mechanisms are unknown. Furthermore, there are no genomic resources of this herb. In this article, we provide sequence and transcript abundance data for the root, stem, and leaf transcriptome of C. grandiflora Benth obtained by the Illumina Hiseq2000. More than 438 million clean reads were obtained from root, stem, and leaf libraries, which produced 153,198 unigenes. Based on databases annotation, a total of 557, 213, and 161 unigenes were annotated to catalpol, acteoside, and azafrin biosynthetic pathways, respectively. Differentially expressed gene analysis identified 14,875 unigenes differentially enriched between leaf and root with 8,054 upregulated genes and 6,821 downregulated genes. Candidate MYB transcription factors involved in catalpol, acteoside, and azafrin biosynthesis were also predicated. This work is the first transcriptome analysis in C. grandiflora Benth which will aid the deciphering of biosynthesis pathways and regulatory mechanisms of active components.

Transcriptome assembly and abiotic related gene expression analysis of date palm reveal candidate genes involved in response to cadmium stress

Date palm in Tunisia is of major economic importance but are also factors of social, environmental and economic stability. An annotated assembly of the transcriptome of cultivar Deglet Nour was reported. RNA was isolated from plant Cd-contaminated leaves, and 37,049 unique Illumina RNA-seq reads were used in a transcriptome assembly. The draft transcriptome assembly consists of 6789 contigs and 17.285 singletons with a means length of 858 bp and 1.042 bp, respectively. The final assembly was functionally annotated using Blast2GO software, allowing the identification of putative genes controlling important agronomic traits. The annotated transcriptome data sets were used to query all known Kyoto Encyclopedia of Genes and Genomes pathways. The most represented molecular functions and biological processes were nucleotide binding and transcription, transport and response to stress and abiotic and biotic stimuli. A prediction of the genes interaction network was proposed by selecting corresponding functionally similar genes from Arabidopsis datasets, downloaded by GeneMANIA version 2.1. Several Cd-responsive genes expression was monitored in in vitro isolated explant of Cd stressed Deglet Nour. Some chelators encoding genes were upregulated confirming in silico findings. Genes encoding HMs transporters in date palm showed expression enhancement more pronounced after 20 days of exposure. P. dactylifera transcriptome provides a valuable resource for future functional analysis of candidate genes involved in metal stress response.

Transcriptome Profiling Provides Insight into the Genes in Carotenoid Biosynthesis during the Mesocarp and Seed Developmental Stages of Avocado (Persea americana)

Avocado (Persea americana Mill.) is an economically important crop because of its high nutritional value. However, the absence of a sequenced avocado reference genome has hindered investigations of secondary metabolism. For next-generation high-throughput transcriptome sequencing, we obtained 365,615,152 and 348,623,402 clean reads as well as 109.13 and 104.10 Gb of sequencing data for avocado mesocarp and seed, respectively, during five developmental stages. High-quality reads were assembled into 100,837 unigenes with an average length of 847.40 bp (N50 = 1725 bp). Additionally, 16,903 differentially expressed genes (DEGs) were detected, 17 of which were related to carotenoid biosynthesis. The expression levels of most of these 17 DEGs were higher in the mesocarp than in the seed during five developmental stages. In this study, the avocado mesocarp and seed transcriptome were also sequenced using single-molecule long-read sequencing to acquired 25.79 and 17.67 Gb clean data, respectively. We identified 233,014 and 238,219 consensus isoforms in avocado mesocarp and seed, respectively. Furthermore, 104 and 59 isoforms were found to correspond to the putative 11 carotenoid biosynthetic-related genes in the avocado mesocarp and seed, respectively. The isoform numbers of 10 out of the putative 11 genes involved in the carotenoid biosynthetic pathway were higher in the mesocarp than those in the seed. Besides, alpha- and beta-carotene contents in the avocado mesocarp and seed during five developmental stages were also measured, and they were higher in the mesocarp than in the seed, which validated the results of transcriptome profiling. Gene expression changes and the associated variations in gene dosage could influence carotenoid biosynthesis. These results will help to further elucidate carotenoid biosynthesis in avocado.

Behind the Myth of the Fruit of Heaven, a Critical Review on Gac (Momordica cochinchinensis Spreng.) Contribution to Nutrition

Gac, Momordica cochinchinensis (Lour.) Spreng. belongs to the Cucurbitaceae family. It is more considered as a super fruit. The demand for this plant is growing in countries where its reputation is high, including traditional countries of gac culture and countries fond of super fruits and food supplements. In these latter countries, the industrial strategy aims at producing high added value in food supplements or nutritional rich preparations. However, when marketing is not the driving force and claims have to be related to scientific data, the situation of gac is less "heavenly", mainly because its most remarkable properties are in the field of micronutrients. These latter components are indeed very important for health but their supplementation on healthy populations brings no significant advantage. This paper proposes to review aspects important for the nutritional reputation of this plant: where it comes from, how it is cultured to have an optimal nutritional composition, what is its composition and how it can impact health of consumers, in which products it is used and what are the regulations to use it in different markets. One important goal of this review is to give a critical and scientific approach to confirm data on this fruit, which has been promoted by marketing departments injecting so many wrong and unverified information. Missing data will be highlighted and potential positive applications are proposed all along the text.

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of β-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

De Novo Assembled Transcriptome Analysis and Identification of Genic SSR Markers in Red-Flowered Strawberry

Red-flowered strawberry is a new ornamental flower derived from intergeneric hybridization (Fragaria × Potentilla). To date, few molecular markers have been reported for this plant. RNA sequencing provides a relatively fast and low-cost approach for large-scale detection of simple sequence repeats (SSRs). In the present study, we profiled the transcriptome of red-flowered strawberry by Illumina HiSeq 2500 to identify SSRs related to petal color. Based on 2 million clean reads of red and white flowers from red-flowered strawberry hybrids, we assembled 91,835 unigenes with an average length of 717 bp. After functional annotation and prediction, there were 47,058 unigenes; of these, 26,861 had a gene ontology annotation, with 14,264 SSR loci. Mononucleotide SSRs were the predominant repeat type (47.20%, n = 6724), followed by di- (32.50%, n = 4641), tri- (19.10%, n = 2729), tetra- (0.90%, n = 132), hexa- (0.2%, n = 21), and penta- (0.10%, n = 16) nucleotide repeats. The most frequent di-, tri-, and tetra-nucleotide repeats were AG/CT, AAG/CTT, and AAAG/CTTT, respectively. PCR amplification with 105 SSR primer pairs yielded four bands specific to red flowers, namely UgRFsr57622, UgRFsr94149, UgRFsr40142, and UgRFsr54608; corresponding 4 trait-specific markers were found to co-segregate with white and red flower color in hybrid population, demonstrating that the genic SSR marker is useful to discriminate between white and red flowers in strawberry. Markers to discriminate flower color in red-flowered strawberry will be useful for early selection of progeny and for breeding management.

De novo transcriptome sequencing and metabolite profiling analyses reveal the complex metabolic genes involved in the terpenoid biosynthesis in Blue Anise Sage (Salvia guaranitica L.)

Many terpenoid compounds have been extracted from different tissues of Salvia guaranitica. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 4 Gb of raw data were generated from the transcriptome of S. guaranitica leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 32 million, comprising 186 million of high-quality nucleotide bases. A total of 61,400 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 267 unigenes that are putatively involved in terpenoid metabolism (including, 198 mevalonate and methyl-erythritol phosphate (MEP) pathways, terpenoid backbone biosynthesis genes and 69 terpene synthases genes). Moreover, three terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic Arabidopsis; most transgenic Arabidopsis plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich blue anise sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. guaranitica.

A transcriptome screen for positive selection in domesticated breadfruit and its wild relatives (Artocarpus spp.)

PREMISE OF THE STUDY

Underutilized crops, such as breadfruit (Artocarpus altilis, Moraceae) have the potential to improve global food security. Humans have artificially selected many cultivars of breadfruit since its domestication began approximately 3500 years ago. The goal of this research was to identify transcriptomic signals of positive selection and to develop genomic resources that may facilitate the development of improved breadfruit cultivars in the future.

METHODS

A reference transcriptome of breadfruit was assembled de novo and annotated. Twenty-four transcriptomes of breadfruit and its wild relatives were generated and analyzed to reveal signals of positive selection that may have resulted from local adaptation or natural selection. Emphasis was placed on MADS-box genes, which are important because they often regulate fruiting timing and structures, and on carotenoid biosynthesis genes, which can impact the nutritional quality of the fruit.

KEY RESULTS

Over 1000 genes showed signals of positive selection, and these genes were enriched for localization to plastids. Nucleotide sites and individuals under positive selection were discovered in MADS-box genes and carotenoid biosynthesis genes, with several sites located in cofactor or DNA-binding domains. A McDonald-Kreitman test comparing wild to cultivated samples revealed selection in one of the carotenoid biosynthesis genes, abscisic acid 8'-hydroxylase 3.

CONCLUSIONS

This research highlights some of the many genes that may have been intentionally or unintentionally selected for during the human-mediated dispersal of breadfruit and stresses the importance of conserving a varied germplasm collection. It has revealed candidate genes for further study and produced new genomic resources for breadfruit.

Identification of 'Xinlimei' radish candidate genes associated with anthocyanin biosynthesis based on a transcriptome analysis

Radish is an economically important vegetable crop belonging to the family Brassicaceae. The high anthocyanin content of the 'Xinlimei' radish roots has been associated with diverse health benefits. However, there is a lack of transcript-level information regarding anthocyanin biosynthesis. In the present study, the 'Xinlimei' radish root transcriptome was analyzed by RNA sequencing at five developmental stages. A total of 222,384,034 clean reads were obtained and 32,253 unigenes were annotated. Expression profiles revealed 10,890 differentially expressed genes (DEGs) among the five analyzed libraries. The DEGs were predominantly involved in KEGG pathways related to the biosynthesis of phenylpropanoids, flavonoids, flavone, and flavonol. The transcriptome data revealed 44 structural and 182 transcription factor genes (TFs) associated with the anthocyanin biosynthetic pathway. Ten structural genes (i.e., 4CL3, CHSB2, CHS1, CHS3, F3H1, F3'H, DFR, DFR1, ANS, and UFGT) and two MYB genes, which were highly and differentially expressed during root development, may be critical for anthocyanin biosynthesis. Additionally, the co-expression of TFs and structural genes was analyzed. Three structural genes (i.e., DFR, ANS, and UFGT) were validated by molecular cloning. The qRT-PCR results indicated that the expression profiles of DEGs were generally consistent with the high-throughput sequencing results. These findings helped identify candidate genes involved in anthocyanin biosynthesis and may be useful for clarifying the molecular mechanism underlying the accumulation of anthocyanins in radish roots.

Expression of Carotenoid Biosynthetic Genes and Carotenoid Biosynthesis during Seedling Development of Momordica charantia

Carotenoids belong to a large group of secondary metabolites, and have pivotal roles in plants, including photosynthesis and phytohormone synthesis, pigmentation, and membrane stabilization. Additionally, carotenoids are potent antioxidants, and their health benefits are becoming increasingly prominent. In recent years, carotenoids have been studied in many plants. Furthermore, gene expression, as well as carotenoid accumulation in different parts of the bitter melon, has been investigated; however, it has not been studied in bitter melon seedlings. In this study, carotenoid accumulation and transcript levels of McGGPPS1, McGGPPS2, McPSY, McPDS, McZDS, McLCYB, McLCYE1, McLCYE2, McCXHB, and McZEP, involved in carotenoid biosynthesis, were analyzed during seedling development using HPLC and qRT-PCR. The major carotenoids that accumulated in the bitter melon seedlings were lutein and E-β-carotene. The expression of most carotenoid biosynthetic genes increased during seedling development, consistent with the accumulation of violaxanthin, lutein, zeaxanthin, β-cryptoxanthin, 13Z-β-carotene, E-β-carotene, and 9Z-β-carotene in bitter melon seedlings. The results of this study provide a firm basis for comprehending the link between gene expression and carotenoid concentration in bitter melon seedlings.

Transcriptome profiling of pumpkin (Cucurbita moschata Duch.) leaves infected with powdery mildew

Cucurbit powdery mildew (PM) is one of the most severe fungal diseases, but the molecular mechanisms underlying PM resistance remain largely unknown, especially in pumpkin (Cucurbita moschata Duch.). The goal of this study was to identify gene expression differences in PM-treated plants (harvested at 24 h and 48 h after inoculation) and untreated (control) plants of inbred line “112–2” using RNA sequencing (RNA-Seq). The inbred line “112–2” has been purified over 8 consecutive generations of self-pollination and shows high resistance to PM. More than 7600 transcripts were examined in pumpkin leaves, and 3129 and 3080 differentially expressed genes (DEGs) were identified in inbred line “112–2” at 24 and 48 hours post inoculation (hpi), respectively. Based on the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database and GO (Gene Ontology) database, a complex regulatory network for PM resistance that may involve hormone signal transduction pathways, transcription factors and defense responses was revealed at the transcription level. In addition, the expression profiles of 16 selected genes were analyzed using quantitative RT-PCR. Among these genes, the transcript levels of 6 DEGs, including bHLH87 (Basic Helix-loop-helix transcription factor), ERF014 (Ethylene response factor), WRKY21 (WRKY domain), HSF (heat stress transcription factor A), MLO3 (Mildew Locus O), and SGT1 (Suppressor of G-Two Allele of Skp1), in PM-resistant “112–2” were found to be significantly up- or down-regulated both before 9 hpi and at 24 hpi or 48 hpi; this behavior differed from that observed in the PM-susceptible material (cultivar “Jiujiangjiaoding”). The transcriptome data provide novel insights into the response of Cucurbita moschata to PM stress and are expected to be highly useful for dissecting PM defense mechanisms in this major vegetable and for improving pumpkin breeding with enhanced resistance to PM.

Transcriptome and metabolome profiling of Narcissus pseudonarcissus 'King Alfred' reveal components of Amaryllidaceae alkaloid metabolism

Amaryllidaceae alkaloids (AAs) represent a diverse class of plant specialized metabolites and many display potent pharmacological activities. The AA metabolic pathway is poorly understood and resources are minimal. To enable AA pathway elucidation and novel biosynthetic enzymes discovery, we generated comprehensive metabolomic and corresponding transcriptomic datasets from different tissues of Narcissus pseudonarcissus ‘King Alfred’. In this study, we performed untargeted UPLC-QTOF-MS metabolite analysis from different tissues, which generated exhaustive list of compounds, including several AAs, most predominant and diverse in bulbs. RNA sequencing of N. pseudonarcissus ‘King Alfred’ bulbs yielded 195,347 transcripts, after assembly. Top expressed genes belong to process like metabolism, survival, and defense including alkaloid biosynthetic genes. The transcriptome contained complete sequences for all proposed genes encoding AA-biosynthetic enzymes such as tyrosine decarboxylase (TYDC1 and TYDC2), phenylalanine ammonia-lyase (PAL1 and PAL2) and phenolic acids hydroxylases (C4H and C3H) to name a few. Furthermore, transcriptome data were validated using RT-qPCR analysis and expression study in different tissues of N. pseudonarcissus ‘King Alfred’ was performed. Here, we present the first comprehensive metabolome and transcriptome study from N. pseudonarcissus ‘King Alfred’ providing invaluable resources for metabolic engineering and biotechnological applications.

Transcriptome and metabolite analyses reveal the complex metabolic genes involved in volatile terpenoid biosynthesis in garden sage (Salvia officinalis)

A large number of terpenoid compounds have been extracted from different tissues of S. officinalis. However, the molecular genetic basis of terpene biosynthesis pathways is virtually unknown. In this study, approximately 6.6 Gb of raw data were generated from the transcriptome of S. officinalis leaves using Illumina HiSeq 2000 sequencing. After filtering and removing the adapter sequences from the raw data, the number of reads reached 21 million, comprising 98 million of high-quality nucleotide bases. 48,671 unigenes were assembled de novo and annotated for establishing a valid database for studying terpenoid biosynthesis. We identified 135 unigenes that are putatively involved in terpenoid metabolism, including 70 mevalonate and methyl-erythritol phosphate pathways, terpenoid backbone biosynthesis genes, and 65 terpene synthase genes. Moreover, five terpene synthase genes were studied for their functions in terpenoid biosynthesis by using transgenic tobacco; most transgenic tobacco plants expressing these terpene synthetic genes produced increased amounts of terpenoids compared with wild-type control. The combined data analyses from the transcriptome and metabolome provide new insights into our understanding of the complex metabolic genes in terpenoid-rich sage, and our study paves the way for the future metabolic engineering of the biosynthesis of useful terpene compounds in S. officinalis.

Medically important carotenoids from Momordica charantia and their gene expressions in different organs

Carotenoids, found in the fruit and different organs of bitter melon (Momordica charantia), have attracted great attention for their potential health benefits in treating several major chronic diseases. Therefore, study related to the identification and quantification of the medically important carotenoid metabolites is highly important for the treatment of various disorderes. The present study involved in the identification and quantification of the various carotenoids present in the different organs of M. charantia and the identification of the genes responsible for the accumulation of the carotenoids with respect to the transcriptome levels were investigated. In this study, using the transcriptome database of bitter melon, a partial-length cDNA clone encoding geranylgeranyl pyrophosphate synthase (McGGPPS2), and several full-length cDNA clones encoding geranylgeranyl pyrophosphate synthase (McGGPPS1), zeta-carotene desaturase (McZDS), lycopene beta-cyclase (McLCYB), lycopene epsilon cyclases (McLCYE1 and McLCYE2), beta-carotene hydroxylase (McCHXB), and zeaxanthin epoxidase (McZEP) were identified in bitter melon. The expression levels of the mRNAs encoding these eight putative biosynthetic enzymes, as well as the accumulation of lycopene, α-carotene, lutein, 13Z-β-carotene, E-β-carotene, 9Z-β-carotene, β-cryptoxanthin, zeaxanthin, antheraxanthin, and violaxanthin were investigated in different organs from M. charantia as well as in the four different stages of its fruit maturation. Transcripts were found to be constitutively expressed at high levels in the leaves where carotenoids were also found at the highest levels. Collectively, these results indicate that the putative McGGPPS2, McZDS, McLCYB, McLCYE1, McLCYE2, and McCHXB enzymes might be key factors in controlling carotenoid content in bitter melon. In conclusion, the over expression of the carotenoid biosynthetic genes from M. charantia crops to increase the yield of these medically important carotenoids.

Accumulation of Charantin and Expression of Triterpenoid Biosynthesis Genes in Bitter Melon (Momordica charantia)

Charantin, a natural cucurbitane type triterpenoid, has been reported to have beneficial pharmacological functions such as anticancer, antidiabetic, and antibacterial activities. However, accumulation of charantin in bitter melon has been little studied. Here, we performed a transcriptome analysis to identify genes involved in the triterpenoid biosynthesis pathway in bitter melon seedlings. A total of 88,703 transcripts with an average length of 898 bp were identified in bitter melon seedlings. On the basis of a functional annotation, we identified 15 candidate genes encoding enzymes related to triterpenoid biosynthesis and analyzed their expression in different organs of mature plants. Most genes were highly expressed in flowers and/or fruit from the ripening stages. An HPLC analysis confirmed that the accumulation of charantin was highest in fruits from the ripening stage, followed by male flowers. The accumulation patterns of charantin coincide with the expression pattern of McSE and McCAS1, indicating that these genes play important roles in charantin biosynthesis in bitter melon. We also investigated optimum light conditions for enhancing charantin biosynthesis in bitter melon and found that red light was the most effective wavelength.

Identification of critical genes associated with lignin biosynthesis in radish (Raphanus sativus L.) by de novo transcriptome sequencing

Radish is an important root vegetable crop with high nutritional, economic, and medicinal value. Lignin is an important secondary metabolite possessing a great effect on plant growth and product quality. To date, lignin biosynthesis-related genes have been identified in some important plant species. However, little information on characterization of critical genes involved in plant lignin biosynthesis is available in radish. In this study, a total of 71,148 transcripts sequences were obtained from radish root, of which 66 assembled unigenes and ten candidate genes were identified to be involved in lignin monolignol biosynthesis. Full-length cDNA sequences of seven randomly selected genes were isolated and sequenced from radish root, and the assembled unigenes covered more than 80% of their corresponding cDNA sequences. Moreover, the lignin content gradually accumulated in leaf during the developmental stages, and it increased from pre-cortex to cortex splitting stage, followed by a decrease at thickening stage and then increased at mature stage in root. RT-qPCR analysis revealed that all these genes except RsF5H exhibited relatively low expression level in root at thickening stage. The expression profiles of Rs4CL5, RsCCoAOMT1, and RsCOMT genes were consistent with the changes of root lignin content, implying that these candidate genes may play important roles in lignin formation in radish root. These findings would provide valuable information for identification of lignin biosynthesis-related genes and facilitate dissection of molecular mechanism underlying lignin biosynthesis in radish and other root vegetable crops.

Comparison of transcriptional expression patterns of carotenoid metabolism in 'Cabernet Sauvignon' grapes from two regions with distinct climate

The downstream flux of carotenoid metabolism in grape berries includes the biosynthesis of norisoprenoids, a group of important aroma compounds, and the production of ABA, a well-known plant hormone. This study focused on the transcriptional profiling comparison of genes participating in the biosynthesis of carotenoids, norisoprenoids, and ABA in Vitis vinifera 'Cabernet Sauvignon' grapes at pea size, veraison, and ripening stages. The grapes were obtained from Changli (CL, eastern China) and Gaotai (GT, western China) regions and analyzed using RNA-sequencing technology. The transcripts required for the carotenoid biosynthesis pathway showed a coordinated expression pattern, mainly expressing at green stage for CL and at veraison for GT, respectively. However, the carotenoid content evolution was not coincident with the timing and pattern of related gene expressions, since more carotenoids were accumulated at veraison in CL relative to two weeks before veraison in GT. Interestingly, norisoprenoid content was higher in GT than in CL, particularly at veraison and ripening, while the key gene encoding carotenoid cleavage dioxygenases, VvCCD1, showed an inverse relationship within the two regions. Higher flux was expected through the carotenoid pathway into ABA production in GT, based on the higher expression level of 9-cis-epoxycarotenoid dioxygenase and drought growing conditions. Most components involved in ABA and ethylene signaling showed distinct expression profiles in the two regions. These results revealed that downstream flux of carotenoid metabolism in grape berries showed regional differences. This study lays a foundation for future research to explore the molecular basis of climatic influences on carotenoid, norisoprenoid, and ABA biosynthesis.

De novo Transcriptome Assembly of Chinese Kale and Global Expression Analysis of Genes Involved in Glucosinolate Metabolism in Multiple Tissues

Chinese kale, a vegetable of the cruciferous family, is a popular crop in southern China and Southeast Asia due to its high glucosinolate content and nutritional qualities. However, there is little research on the molecular genetics and genes involved in glucosinolate metabolism and its regulation in Chinese kale. In this study, we sequenced and characterized the transcriptomes and expression profiles of genes expressed in 11 tissues of Chinese kale. A total of 216 million 150-bp clean reads were generated using RNA-sequencing technology. From the sequences, 98,180 unigenes were assembled for the whole plant, and 49,582~98,423 unigenes were assembled for each tissue. Blast analysis indicated that a total of 80,688 (82.18%) unigenes exhibited similarity to known proteins. The functional annotation and classification tools used in this study suggested that genes principally expressed in Chinese kale, were mostly involved in fundamental processes, such as cellular and molecular functions, the signal transduction, and biosynthesis of secondary metabolites. The expression levels of all unigenes were analyzed in various tissues of Chinese kale. A large number of candidate genes involved in glucosinolate metabolism and its regulation were identified, and the expression patterns of these genes were analyzed. We found that most of the genes involved in glucosinolate biosynthesis were highly expressed in the root, petiole, and in senescent leaves. The expression patterns of ten glucosinolate biosynthetic genes from RNA-seq were validated by quantitative RT-PCR in different tissues. These results provided an initial and global overview of Chinese kale gene functions and expression activities in different tissues.

Transcriptome Analysis of Salicylic Acid Treatment in Rehmannia glutinosa Hairy Roots Using RNA-seq Technique for Identification of Genes Involved in Acteoside Biosynthesis

Rehmannia glutinosa is a common bulk medicinal material that has been widely used in China due to its active ingredients. Acteoside, one of the ingredients, has antioxidant, antinephritic, anti-inflammatory, hepatoprotective, immunomodulatory, and neuroprotective effects, is usually selected as a quality-control component for R. glutinosa herb in the Chinese Pharmacopeia. The acteoside biosynthesis pathway in R. glutinosa has not yet been clearly established. Herein, we describe the establishment of a genetic transformation system for R. glutinosa mediated by Agrobacterium rhizogenes. We screened the optimal elicitors that markedly increased acteoside accumulation in R. glutinosa hairy roots. We found that acteoside accumulation dramatically increased with the addition of salicylic acid (SA); the optimal SA dose was 25 μmol/L for hairy roots. RNA-seq was applied to analyze the transcriptomic changes in hairy roots treated with SA for 24 h in comparison with an untreated control. A total of 3,716, 4,018, and 2,715 differentially expressed transcripts (DETs) were identified in 0 h-vs.-12 h, 0 h-vs.-24 h, and 12 h-vs.-24 h libraries, respectively. KEGG pathway-based analysis revealed that 127 DETs were enriched in "phenylpropanoid biosynthesis." Of 219 putative unigenes involved in acteoside biosynthesis, 54 were found to be up-regulated at at least one of the time points after SA treatment. Selected candidate genes were analyzed by quantitative real-time PCR (qRT-PCR) in hairy roots with SA, methyl jasmonate (MeJA), AgNO₃ (Ag⁺), and putrescine (Put) treatment. All genes investigated were up-regulated by SA treatment, and most candidate genes were weakly increased by MeJA to some degree. Furthermore, transcription abundance of eight candidate genes in tuberous roots of the high-acteoside-content (HA) cultivar QH were higher than those of the low-acteoside-content (LA) cultivar Wen 85-5. These results will pave the way for understanding the molecular basis of acteoside biosynthesis in R. glutinosa, and can serve as a basis for future validation studies.

Transcriptional Profiling and Identification of Heat-Responsive Genes in Perennial Ryegrass by RNA-Sequencing

Perennial ryegrass (Lolium perenne) is one of the most widely used forage and turf grasses in the world due to its desirable agronomic qualities. However, as a cool-season perennial grass species, high temperature is a major factor limiting its performance in warmer and transition regions. In this study, a de novo transcriptome was generated using a cDNA library constructed from perennial ryegrass leaves subjected to short-term heat stress treatment. Then the expression profiling and identification of perennial ryegrass heat response genes by digital gene expression analyses was performed. The goal of this work was to produce expression profiles of high temperature stress responsive genes in perennial ryegrass leaves and further identify the potentially important candidate genes with altered levels of transcript, such as those genes involved in transcriptional regulation, antioxidant responses, plant hormones and signal transduction, and cellular metabolism. The de novo assembly of perennial ryegrass transcriptome in this study obtained more total and annotated unigenes compared to previously published ones. Many DEGs identified were genes that are known to respond to heat stress in plants, including HSFs, HSPs, and antioxidant related genes. In the meanwhile, we also identified four gene candidates mainly involved in C₄ carbon fixation, and one TOR gene. Their exact roles in plant heat stress response need to dissect further. This study would be important by providing the gene resources for improving heat stress tolerance in both perennial ryegrass and other cool-season perennial grass plants.

Transcriptome Sequencing of Gynostemma pentaphyllum to Identify Genes and Enzymes Involved in Triterpenoid Biosynthesis

G. pentaphyllum (Gynostemma pentaphyllum), a creeping herbaceous perennial with many important medicinal properties, is widely distributed in Asia. Gypenosides (triterpenoid saponins), the main effective components of G. pentaphyllum, are well studied. FPS (farnesyl pyrophosphate synthase), SS (squalene synthase), and SE (squalene epoxidase) are the main enzymes involved in the synthesis of triterpenoid saponins. Considering the important medicinal functions of G. pentaphyllum, it is necessary to investigate the transcriptomic information of G. pentaphyllum to facilitate future studies of transcriptional regulation. After sequencing G. pentaphyllum, we obtained 50,654,708 unigenes. Next, we used RPKM (reads per kilobases per million reads) to calculate expression of the unigenes and we performed comparison of our data to that contained in five common databases to annotate different aspects of the unigenes. Finally, we noticed that FPS, SS, and SE showed differential expression of enzymes in DESeq. Leaves showed the highest expression of FPS, SS, and SE relative to the other two tissues. Our research provides transcriptomic information of G. pentaphyllum in its natural environment and we found consistency in unigene expression, enzymes expression (FPS, SS, and SE), and the distribution of gypenosides content in G. pentaphyllum. Our results will enable future related studies of G. pentaphyllum.

Comparative Transcriptome Analyses Reveal a Special Glucosinolate Metabolism Mechanism in Brassica alboglabra Sprouts

Brassica sprouts contain abundant phytochemicals, especially glucosinolates (GSs). Various methods have been used to enhance GS content in sprouts. However, the molecular basis of GS metabolism in sprouts remains an open question. Here we employed RNA-seq analysis to compare the transcriptomes of high-GS (JL-08) and low-GS (JL-09) Brassica alboglabra sprouts. Paired-end Illumina RNA-seq reads were generated and mapped to the Brassica oleracea reference genome. The differentially expressed genes were analyzed between JL-08 and JL-09. Among these, 1477 genes were up-regulated and 1239 down-regulated in JL-09 compared with JL-08. Enrichment analysis of these differentially expressed genes showed that the GS biosynthesis had the smallest enrichment factor and the highest Q-value of all metabolic pathways in Kyoto Encyclopedia of Genes and Genomes database, indicating the main metabolic difference between JL-08 and JL-09 is the GS biosynthetic pathway. Thirty-seven genes of the sequenced data were annotated as putatively involved in GS biosynthesis, degradation, and regulation, of which 11 were differentially expressed in JL-08 and JL-09. The expression level of GS degradation enzyme myrosinase in high-GS JL-08 was lower compared with low-GS JL-09. Surprisingly, in high-GS JL-08, the expression levels of GS biosynthesis genes were also lower than those in low-GS JL-09. As the GS contents in sprouts are determined by dynamic equilibrium of seed stored GS mobilization, de novo synthesis, degradation, and extra transport, the result of this study leads us to suggest that efforts to increase GS content should focus on either raising GS content in seeds or decreasing myrosinase activity, rather than improving the expression level of GS biosynthesis genes in sprouts.

Regulation of Carotenoid Biosynthesis During Fruit Development

Carotenoids are recognized as the main pigments in most fruit crops, providing colours that range from yellow and pink to deep orange and red. Moreover, the edible portion of widely consumed fruits or their derived products represent a major dietary source of carotenoids for animals and humans. Therefore, these pigments are crucial compounds contributing to fruit aesthetic and nutritional quality but may also have protecting and ecophysiological functions in coloured fruits. Among plant organs, fruits display one of the most heterogeneous carotenoids patterns in terms of diversity and abundance. In this chapter a comprehensive list of the carotenoid content and profile in the most commonly cultivated fleshy fruits is reported. The proposed fruit classification systems attending to carotenoid composition are revised and discussed. The regulation of carotenoids in fruits can be rather complex due to the dramatic changes in content and composition during ripening, which are also dependent on the fruit tissue and the developmental stage. In addition, carotenoid accumulation is a dynamic process, associated with the development of chromoplasts during ripening. As a general rule, carotenoid accumulation is highly controlled at the transcriptional level of the structural and accessory proteins of the biosynthetic and degradation pathways, but other mechanisms such as post-transcriptional modifications or the development of sink structures have been recently revealed as crucial factors in determining the levels and stability of these pigments. In this chapter common key metabolic reactions regulating carotenoid composition in fruit tissues are described in addition to others that are restricted to certain species and generate unique carotenoids patterns. The existence of fruit-specific isoforms for key steps such as the phytoene synthase, lycopene β-cyclases or catabolic carotenoid cleavage dioxygenases has allowed an independent regulation of the pathway in fruit tissues and a source of variability to create novel activities or different catalytic properties. Besides key genes of the carotenoid pathway, changes in carotenoid accumulation could be also directly influenced by differences in gene expression or protein activity in the pathway of carotenoid precursors and some relevant examples are discussed. The objective of this chapter is to provide an updated review of the main carotenoid profiles in fleshy fruits, their pattern of changes during ripening and our current understanding of the different regulatory levels responsible for the diversity of carotenoid accumulation in fruit tissues.

De novo Taproot Transcriptome Sequencing and Analysis of Major Genes Involved in Sucrose Metabolism in Radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is an important annual or biennial root vegetable crop. The fleshy taproot comprises the main edible portion of the plant with high nutrition and medical value. Molecular biology study of radish begun rather later, and lacks sufficient transcriptomic and genomic data in pubic databases for understanding of the molecular mechanism during the radish taproot formation. To develop a comprehensive overview of the 'NAU-YH' root transcriptome, a cDNA library, prepared from three equally mixed RNA of taproots at different developmental stages including pre-cortex splitting stage, cortex splitting stage, and expanding stage was sequenced using high-throughput Illumina RNA sequencing. From approximately 51 million clean reads, a total of 70,168 unigenes with a total length of 50.28 Mb, an average length of 717 bp and a N50 of 994 bp were obtained. In total, 63,991 (about 91.20% of the assembled unigenes) unigenes were successfully annotated in five public databases including NR, GO, COG, KEGG, and Nt. GO analysis revealed that the majority of these unigenes were predominately involved in basic physiological and metabolic processes, catalytic, binding, and cellular process. In addition, a total of 103 unigenes encoding eight enzymes involved in the sucrose metabolism related pathways were also identified by KEGG pathway analysis. Sucrose synthase (29 unigenes), invertase (17 unigenes), sucrose-phosphate synthase (16 unigenes), fructokinase (17 unigenes), and hexokinase (11 unigenes) ranked top five in these eight key enzymes. From which, two genes (RsSuSy1, RsSPS1) were validated by T-A cloning and sequenced, while the expression of six unigenes were profiled with RT-qPCR analysis. These results would be served as an important public reference platform to identify the related key genes during taproot thickening and facilitate the dissection of molecular mechanisms underlying taproot formation in radish.

De novo assembly and functional annotation of Myrciaria dubia fruit transcriptome reveals multiple metabolic pathways for L-ascorbic acid biosynthesis

Background

Myrciaria dubia is an Amazonian fruit shrub that produces numerous bioactive phytochemicals, but is best known by its high L-ascorbic acid (AsA) content in fruits. Pronounced variation in AsA content has been observed both within and among individuals, but the genetic factors responsible for this variation are largely unknown. The goals of this research, therefore, were to assemble, characterize, and annotate the fruit transcriptome of M. dubia in order to reconstruct metabolic pathways and determine if multiple pathways contribute to AsA biosynthesis.

Results

In total 24,551,882 high-quality sequence reads were de novo assembled into 70,048 unigenes (mean length = 1150 bp, N50 = 1775 bp). Assembled sequences were annotated using BLASTX against public databases such as TAIR, GR-protein, FB, MGI, RGD, ZFIN, SGN, WB, TIGR_CMR, and JCVI-CMR with 75.2 % of unigenes having annotations. Of the three core GO annotation categories, biological processes comprised 53.6 % of the total assigned annotations, whereas cellular components and molecular functions comprised 23.3 and 23.1 %, respectively. Based on the KEGG pathway assignment of the functionally annotated transcripts, five metabolic pathways for AsA biosynthesis were identified: animal-like pathway, myo-inositol pathway, L-gulose pathway, D-mannose/L-galactose pathway, and uronic acid pathway. All transcripts coding enzymes involved in the ascorbate-glutathione cycle were also identified. Finally, we used the assembly to identified 6314 genic microsatellites and 23,481 high quality SNPs.

Conclusions

This study describes the first next-generation sequencing effort and transcriptome annotation of a non-model Amazonian plant that is relevant for AsA production and other bioactive phytochemicals. Genes encoding key enzymes were successfully identified and metabolic pathways involved in biosynthesis of AsA, anthocyanins, and other metabolic pathways have been reconstructed. The identification of these genes and pathways is in agreement with the empirically observed capability of M. dubia to synthesize and accumulate AsA and other important molecules, and adds to our current knowledge of the molecular biology and biochemistry of their production in plants. By providing insights into the mechanisms underpinning these metabolic processes, these results can be used to direct efforts to genetically manipulate this organism in order to enhance the production of these bioactive phytochemicals.

The accumulation of AsA precursor and discovery of genes associated with their biosynthesis and metabolism in M. dubia is intriguing and worthy of further investigation. The sequences and pathways produced here present the genetic framework required for further studies. Quantitative transcriptomics in concert with studies of the genome, proteome, and metabolome under conditions that stimulate production and accumulation of AsA and their precursors are needed to provide a more comprehensive view of how these pathways for AsA metabolism are regulated and linked in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2225-6) contains supplementary material, which is available to authorized users.

De novo transcriptome sequencing in Bixa orellana to identify genes involved in methylerythritol phosphate, carotenoid and bixin biosynthesis

Background

Bixin or annatto is a commercially important natural orange-red pigment derived from lycopene that is produced and stored in seeds of Bixa orellana L. An enzymatic pathway for bixin biosynthesis was inferred from homology of putative proteins encoded by differentially expressed seed cDNAs. Some activities were later validated in a heterologous system. Nevertheless, much of the pathway remains to be clarified. For example, it is essential to identify the methylerythritol phosphate (MEP) and carotenoid pathways genes.

Results

In order to investigate the MEP, carotenoid, and bixin pathways genes, total RNA from young leaves and two different developmental stages of seeds from B. orellana were used for the construction of indexed mRNA libraries, sequenced on the Illumina HiSeq 2500 platform and assembled de novo using Velvet, CLC Genomics Workbench and CAP3 software. A total of 52,549 contigs were obtained with average length of 1,924 bp. Two phylogenetic analyses of inferred proteins, in one case encoded by thirteen general, single-copy cDNAs, in the other from carotenoid and MEP cDNAs, indicated that B. orellana is closely related to sister Malvales species cacao and cotton. Using homology, we identified 7 and 14 core gene products from the MEP and carotenoid pathways, respectively. Surprisingly, previously defined bixin pathway cDNAs were not present in our transcriptome. Here we propose a new set of gene products involved in bixin pathway.

Conclusion

The identification and qRT-PCR quantification of cDNAs involved in annatto production suggest a hypothetical model for bixin biosynthesis that involve coordinated activation of some MEP, carotenoid and bixin pathway genes. These findings provide a better understanding of the mechanisms regulating these pathways and will facilitate the genetic improvement of B. orellana.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2065-4) contains supplementary material, which is available to authorized users.

De novo sequencing of root transcriptome reveals complex cadmium-responsive regulatory networks in radish (Raphanus sativus L.)

Cadmium (Cd) is a nonessential metallic trace element that poses potential chronic toxicity to living organisms. To date, little is known about the Cd-responsive regulatory network in root vegetable crops including radish. In this study, 31,015 unigenes representing 66,552 assembled unique transcripts were isolated from radish root under Cd stress based on de novo transcriptome assembly. In all, 1496 differentially expressed genes (DEGs) consisted of 3579 transcripts were identified from Cd-free (CK) and Cd-treated (Cd200) libraries. Gene Ontology and pathway enrichment analysis indicated that the up- and down-regulated DEGs were predominately involved in glucosinolate biosynthesis as well as cysteine and methionine-related pathways, respectively. RT-qPCR showed that the expression profiles of DEGs were in consistent with results from RNA-Seq analysis. Several candidate genes encoding phytochelatin synthase (PCS), metallothioneins (MTs), glutathione (GSH), zinc iron permease (ZIPs) and ABC transporter were responsible for Cd uptake, accumulation, translocation and detoxification in radish. The schematic model of DEGs and microRNAs-involved in Cd-responsive regulatory network was proposed. This study represents a first comprehensive transcriptome-based characterization of Cd-responsive DEGs in radish. These results could provide fundamental insight into complex Cd-responsive regulatory networks and facilitate further genetic manipulation of Cd accumulation in root vegetable crops.

Transcriptome profiling provides new insights into the formation of floral scent in Hedychium coronarium

BACKGROUND

Hedychium coronarium is a popular ornamental plant in tropical and subtropical regions because its flowers not only possess intense and inviting fragrance but also enjoy elegant shape. The fragrance results from volatile terpenes and benzenoids presented in the floral scent profile. However, in this species, even in monocots, little is known about the underlying molecular mechanism of floral scent production.

RESULTS

Using Illumina platform, approximately 81 million high-quality reads were obtained from a pooled cDNA library. The de novo assembly resulted in a transcriptome with 65,591 unigenes, 50.90% of which were annotated using public databases. Digital gene expression (DGE) profiling analysis revealed 7,796 differential expression genes (DEGs) during petal development. GO term classification and KEGG pathway analysis indicated that the levels of transcripts changed significantly in "metabolic process", including "terpenoid biosynthetic process". Through a systematic analysis, 35 and 33 candidate genes might be involved in the biosynthesis of floral volatile terpenes and benzenoids, respectively. Among them, flower-specific HcDXS2A, HcGPPS, HcTPSs, HcCNL and HcBCMT1 might play critical roles in regulating the formation of floral fragrance through DGE profiling coupled with floral volatile profiling analyses. In vitro characterization showed that HcTPS6 was capable of generating β-farnesene as its main product. In the transcriptome, 1,741 transcription factors (TFs) were identified and 474 TFs showed differential expression during petal development. It is supposed that two R2R3-MYBs with flower-specific and developmental expression might be involved in the scent production.

CONCLUSIONS

The novel transcriptome and DGE profiling provide an important resource for functional genomics studies and give us a dynamic view of biological process during petal development in H. coronarium. These data lay the basis for elucidating the molecular mechanism of floral scent formation and regulation in monocot. The results also provide the opportunities for genetic modification of floral scent profile in Hedychium.

De novo assembly of transcriptome sequencing in Caragana korshinskii Kom. and characterization of EST-SSR markers

Caragana korshinskii Kom. is widely distributed in various habitats, including gravel desert, clay desert, fixed and semi-fixed sand, and saline land in the Asian and African deserts. To date, no previous genomic information or EST-SSR marker has been reported in Caragana Fabr. genus. In this study, more than two billion bases of high-quality sequence of C. korshinskii were generated by using illumina sequencing technology and demonstrated the de novo assembly and annotation of genes without prior genome information. These reads were assembled into 86,265 unigenes (mean length = 709 bp). The similarity search indicated that 33,955 and 21,978 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Among these annotated unigenes, 26,232 a unigenes were separately assigned to Gene Ontology (GO) database. When 22,756 unigenes searched against the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database, 5,598 unigenes were assigned to 5 main categories including 32 KEGG pathways. Among the main KEGG categories, metabolism was the biggest category (2,862, 43.7%), suggesting the active metabolic processes in the desert tree. In addition, a total of 19,150 EST-SSRs were identified from 15,484 unigenes, and the characterizations of EST-SSRs were further compared with other four species in Fabraceae. 126 potential marker sites were randomly selected to validate the assembly quality and develop EST-SSR markers. Among the 9 germplasms in Caranaga Fabr. genus, PCR success rate were 93.7% and the phylogenic tree was constructed based on the genotypic data. This research generated a substantial fraction of transcriptome sequences, which were very useful resources for gene annotation and discovery, molecular markers development, genome assembly and annotation. The EST-SSR markers identified and developed in this study will facilitate marker-assisted selection breeding.

Transcriptome analysis reveals in vitro cultured Withania somnifera leaf and root tissues as a promising source for targeted withanolide biosynthesis

Background

The production of metabolites via in vitro culture is promoted by the availability of fully defined metabolic pathways. Withanolides, the major bioactive phytochemicals of Withania somnifera, have been well studied for their pharmacological activities. However, only a few attempts have been made to identify key candidate genes involved in withanolide biosynthesis. Understanding the steps involved in withanolide biosynthesis is essential for metabolic engineering of this plant to increase withanolide production.

Results

Transcriptome sequencing was performed on in vitro adventitious root and leaf tissues using the Illumina platform. We obtained a total of 177,156 assembled transcripts with an average unigene length of 1,033 bp. About 13% of the transcripts were unique to in vitro adventitious roots but no unique transcripts were observed in in vitro-grown leaves. A putative withanolide biosynthetic pathway was deduced by mapping the assembled transcripts to the KEGG database, and the expression of candidate withanolide biosynthesis genes -were validated by qRT PCR. The accumulation pattern of withaferin A and withanolide A varied according to the type of tissue and the culture period. Further, we demonstrated that in vitro leaf extracts exhibit anticancer activity against human gastric adenocarcinoma cell lines at sub G1 phase.

Conclusions

We report here a validated large-scale transcriptome data set and the potential biological activity of in vitro cultures of W. somnifera. This study provides important information to enhance tissue-specific expression and accumulation of secondary metabolites, paving the way for industrialization of in vitro cultures of W. somnifera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1214-0) contains supplementary material, which is available to authorized users.

Structural characterization of highly glucosylated crocins and regulation of their biosynthesis during flower development in Crocus

Crocin biosynthesis in Crocus has been proposed to proceed through a zeaxanthin cleavage pathway catalyzed by carotenoid cleavage dioxygenase 2 (CCD2), and followed by glucosylation reactions catalyzed by CsGT2 (UGT74AD1). In Crocus ancyrensis flowers, crocins with eight (crocin-1), seven (crocin-2), and six glucose (crocin-3) moieties accumulated both in stigma and tepals. We have characterized the structure of these highly glucosylated crocins and follow up their accumulation by high-resolution liquid chromatography coupled with diode array detector along the development of both tissues, and coupled to the isolation and analysis of the expression of eighteen genes (PSY-I, PSY-II, PDS-(I-V), ISO-ZDS, ZDS, CtrISO, LYC-I and II, BCH, CaCCD2, UGT74AD2-5) related with the apocarotenoid metabolism in C. ancyrensis tepals and stigmas. Structure elucidation of crocin-1 and crocin-2 was done by the combined use of 1D and 2D [(1)H, (1)H] (gCOSY and TOCSY and ROESY) and [(1)H-(13)C] NMR experiments, revealing that for crocin-1 was all-trans-crocetin O-[β-D- Glucopyranosyl)-(1→4)-(β-D-glucopyranosyl)-(1→2)]-O-[β-D-glucopyranosyl-(1→6)]-β-D-glucopyranosyl diester, while crocin-2 showed an identical structure except for the absence of one glucose residue in one end of the molecule. Crocins accumulation was not synchronically regulated in stigma and tepals, although in both cases crocins accumulation parallels tissue development, decreasing at anthesis. The expression of the carotenogenic genes PSY, ZDS-V, BCH, and LCY-II was correlated with crocins accumulation. In addition, CaCCD2 and only one of the four glucosyltransferase encoding genes, UGT74AD2, were highly expressed, and the expression was correlated with high levels of crocins accumulation in stigma and tepals.

Comparison of distinct transcriptional expression patterns of flavonoid biosynthesis in Cabernet Sauvignon grapes from east and west China

Flavonoids make a very important contribution to the organoleptic qualities of grapes and wines. In this work these were analyzed in Cabernet Sauvignon grown in Changli, Hebei Province in east China and Gaotai, Gansu Province in west China. These regions have distinctly different climates contributing to their different 'terroir'. RNA sequencing was performed to trace transcriptome changes in Cabernet Sauvignon berries at pea size, veraison and ripening, corresponding to E-L 31, 35 and 38. The accumulation of flavonols, flavan-3-ols and anthocyanins together with the expression of relevant genes were analyzed and compared between the two regions. The biosynthesis patterns were similar between two regions, but more flavonols, anthocyanins, and tri-hydroxylated flavonoids accumulated in grapes from Gaotai before berry harvest, possibly due to the higher transcript levels of the genes that encode biosynthetic enzymes and their potential candidate transcription factors. The lower levels of flavan-3-ols, mainly (-)-epigallocatechin, in the pre-veraison grapes from Changli, might be due to limited flow of carbon to the F3'5'H branch pathway, as the ratio of F3'5'H to F3'H was lower in these berries from Changli. It is suggested that the combination of climatic factors profoundly affect the flavonoid pathway in grapes from China, providing regionally specific metabolism patterns.

RNA-seq analysis of Rubus idaeus cv. Nova: transcriptome sequencing and de novo assembly for subsequent functional genomics approaches

Using Illumina sequencing technology, we have generated the large-scale transcriptome sequencing data containing abundant information on genes involved in the metabolic pathways in R. idaeus cv. Nova fruits. Rubus idaeus (Red raspberry) is one of the important economical crops that possess numerous nutrients, micronutrients and phytochemicals with essential health benefits to human. The molecular mechanism underlying the ripening process and phytochemical biosynthesis in red raspberry is attributed to the changes in gene expression, but very limited transcriptomic and genomic information in public databases is available. To address this issue, we generated more than 51 million sequencing reads from R. idaeus cv. Nova fruit using Illumina RNA-Seq technology. After de novo assembly, we obtained 42,604 unigenes with an average length of 812 bp. At the protein level, Nova fruit transcriptome showed 77 and 68 % sequence similarities with Rubus coreanus and Fragaria versa, respectively, indicating the evolutionary relationship between them. In addition, 69 % of assembled unigenes were annotated using public databases including NCBI non-redundant, Cluster of Orthologous Groups and Gene ontology database, suggesting that our transcriptome dataset provides a valuable resource for investigating metabolic processes in red raspberry. To analyze the relationship between several novel transcripts and the amounts of metabolites such as γ-aminobutyric acid and anthocyanins, real-time PCR and target metabolite analysis were performed on two different ripening stages of Nova. This is the first attempt using Illumina sequencing platform for RNA sequencing and de novo assembly of Nova fruit without reference genome. Our data provide the most comprehensive transcriptome resource available for Rubus fruits, and will be useful for understanding the ripening process and for breeding R. idaeus cultivars with improved fruit quality.

De novo transcriptome analysis of the Siberian apricot (Prunus sibirica L.) and search for potential SSR markers by 454 pyrosequencing

The Siberian apricot, an economically and ecologically important plant in China, contains seeds high in oil and can grow on marginal land. Although this species has multiple purposes and may be a feedstock of biofuel in China, transcriptome information and molecular research on this species remain limited. RNA-Seq technology has been widely applied to transcriptomics, genomics and the development of molecular markers, and functional gene studies. In this study, we obtained 1,243,067 high-quality reads with a mean size of 425 bp in a single run, totaling 528.4 Mb of sequence data using 454 GS FLX Titanium sequencing. All reads were assembled de novo into 46,940 unigenes with a mean size of 651 bp (range: 45-5566 bp). Assembled unigenes were annotated in multiple public databases based on similarity alignments to genes and proteins. 191 unigenes involving in lipid biosynthesis and metabolism were found, among them, expression patterns of two desaturase enzymes were analyzed by quantitative real-time polymerase chain reaction (qRT-PCR), based on six tissues from Siberian apricot, the seeds had the highest expression. 7304 simple sequence repeats (SSR) were identified from 6509 unigenes, a total of 9930 primer pairs were designed, 50 primer pairs were randomly selected to validate of the usefulness, and 24 (48%) primer pairs produced bands of the expected size. These data provide a base of sequence information to improve agronomic characters and molecular marker-assisted breeding to alter the composition of fatty acids in seeds from this plant, and hence, facilitate its utilization as a future biodiesel feedstock.

De-novo RNA sequencing and metabolite profiling to identify genes involved in anthocyanin biosynthesis in Korean black raspberry (Rubus coreanus Miquel)

The Korean black raspberry (Rubus coreanus Miquel, KB) on ripening is usually consumed as fresh fruit, whereas the unripe KB has been widely used as a source of traditional herbal medicine. Such a stage specific utilization of KB has been assumed due to the changing metabolite profile during fruit ripening process, but so far molecular and biochemical changes during its fruit maturation are poorly understood. To analyze biochemical changes during fruit ripening process at molecular level, firstly, we have sequenced, assembled, and annotated the transcriptome of KB fruits. Over 4.86 Gb of normalized cDNA prepared from fruits was sequenced using Illumina HiSeq™ 2000, and assembled into 43,723 unigenes. Secondly, we have reported that alterations in anthocyanins and proanthocyanidins are the major factors facilitating variations in these stages of fruits. In addition, up-regulation of F3'H1, DFR4 and LDOX1 resulted in the accumulation of cyanidin derivatives during the ripening process of KB, indicating the positive relationship between the expression of anthocyanin biosynthetic genes and the anthocyanin accumulation. Furthermore, the ability of RcMCHI2 (R. coreanus Miquel chalcone flavanone isomerase 2) gene to complement Arabidopsis transparent testa 5 mutant supported the feasibility of our transcriptome library to provide the gene resources for improving plant nutrition and pigmentation. Taken together, these datasets obtained from transcriptome library and metabolic profiling would be helpful to define the gene-metabolite relationships in this non-model plant.

RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis

BACKGROUND

One of the concerns of assembling de novo transcriptomes is determining the amount of read sequences required to ensure a comprehensive coverage of genes expressed in a particular sample. In this report, we describe the use of Illumina paired-end RNA-Seq (PE RNA-Seq) reads from Hevea brasiliensis (rubber tree) bark to devise a transcript mapping approach for the estimation of the read amount needed for deep transcriptome coverage.

FINDINGS

We optimized the assembly of a Hevea bark transcriptome based on 16 Gb Illumina PE RNA-Seq reads using the Oases assembler across a range of k-mer sizes. We then assessed assembly quality based on transcript N50 length and transcript mapping statistics in relation to (a) known Hevea cDNAs with complete open reading frames, (b) a set of core eukaryotic genes and (c) Hevea genome scaffolds. This was followed by a systematic transcript mapping process where sub-assemblies from a series of incremental amounts of bark transcripts were aligned to transcripts from the entire bark transcriptome assembly. The exercise served to relate read amounts to the degree of transcript mapping level, the latter being an indicator of the coverage of gene transcripts expressed in the sample. As read amounts or datasize increased toward 16 Gb, the number of transcripts mapped to the entire bark assembly approached saturation. A colour matrix was subsequently generated to illustrate sequencing depth requirement in relation to the degree of coverage of total sample transcripts.

CONCLUSIONS

We devised a procedure, the "transcript mapping saturation test", to estimate the amount of RNA-Seq reads needed for deep coverage of transcriptomes. For Hevea de novo assembly, we propose generating between 5-8 Gb reads, whereby around 90% transcript coverage could be achieved with optimized k-mers and transcript N50 length. The principle behind this methodology may also be applied to other non-model plants, or with reads from other second generation sequencing platforms.

De novo transcriptome sequencing of radish (Raphanus sativus L.) and analysis of major genes involved in glucosinolate metabolism

Background

Radish (Raphanus sativus L.), is an important root vegetable crop worldwide. Glucosinolates in the fleshy taproot significantly affect the flavor and nutritional quality of radish. However, little is known about the molecular mechanisms underlying glucosinolate metabolism in radish taproots. The limited availability of radish genomic information has greatly hindered functional genomic analysis and molecular breeding in radish.

Results

In this study, a high-throughput, large-scale RNA sequencing technology was employed to characterize the de novo transcriptome of radish roots at different stages of development. Approximately 66.11 million paired-end reads representing 73,084 unigenes with a N50 length of 1,095 bp, and a total length of 55.73 Mb were obtained. Comparison with the publicly available protein database indicates that a total of 67,305 (about 92.09% of the assembled unigenes) unigenes exhibit similarity (e –value ≤ 1.0e^-5) to known proteins. The functional annotation and classification including Gene Ontology (GO), Clusters of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the main activated genes in radish taproots are predominately involved in basic physiological and metabolic processes, biosynthesis of secondary metabolite pathways, signal transduction mechanisms and other cellular components and molecular function related terms. The majority of the genes encoding enzymes involved in glucosinolate (GS) metabolism and regulation pathways were identified in the unigene dataset by targeted searches of their annotations. A number of candidate radish genes in the glucosinolate metabolism related pathways were also discovered, from which, eight genes were validated by T-A cloning and sequencing while four were validated by quantitative RT-PCR expression profiling.

Conclusions

The ensuing transcriptome dataset provides a comprehensive sequence resource for molecular genetics research in radish. It will serve as an important public information platform to further understanding of the molecular mechanisms involved in biosynthesis and metabolism of the related nutritional and flavor components during taproot formation in radish.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-836) contains supplementary material, which is available to authorized users.

Transcriptome comparison reveals key candidate genes responsible for the unusual reblooming trait in tree peonies

Tree peonies are important ornamental plants worldwide, but growing them can be frustrating due to their short and concentrated flowering period. Certain cultivars exhibit a reblooming trait that provides a valuable alternative for extending the flowering period. However, the genetic control of reblooming in tree peonies is not well understood. In this study, we compared the molecular properties and morphology of reblooming and non-reblooming tree peonies during the floral initiation and developmental processes. Using transcriptome sequencing technology, we generated 59,275 and 63,962 unigenes with a mean size of 698 bp and 699 bp from the two types of tree peonies, respectively, and identified eight differentially expressed genes that are involved in the floral pathways of Arabidopsis thaliana. These differentially regulated genes were verified through a detailed analysis of their expression pattern during the floral process by real time RT-PCR. From this combined analysis, we identified four genes, PsFT, PsVIN3, PsCO and PsGA20OX, which likely play important roles in the regulation of the reblooming process in tree peonies. These data constitute a valuable resource for the discovery of genes involved in flowering time and insights into the molecular mechanism of flowering to further accelerate the breeding of tree peonies and other perennial woody plants.

De novo assembly and characterization of fruit transcriptome in Litchi chinensis Sonn and analysis of differentially regulated genes in fruit in response to shading

BACKGROUND

Litchi (Litchi chinensis Sonn.) is one of the most important fruit trees cultivated in tropical and subtropical areas. However, a lack of transcriptomic and genomic information hinders our understanding of the molecular mechanisms underlying fruit set and fruit development in litchi. Shading during early fruit development decreases fruit growth and induces fruit abscission. Here, high-throughput RNA sequencing (RNA-Seq) was employed for the de novo assembly and characterization of the fruit transcriptome in litchi, and differentially regulated genes, which are responsive to shading, were also investigated using digital transcript abundance(DTA)profiling.

RESULTS

More than 53 million paired-end reads were generated and assembled into 57,050 unigenes with an average length of 601 bp. These unigenes were annotated by querying against various public databases, with 34,029 unigenes found to be homologous to genes in the NCBI GenBank database and 22,945 unigenes annotated based on known proteins in the Swiss-Prot database. In further orthologous analyses, 5,885 unigenes were assigned with one or more Gene Ontology terms, 10,234 hits were aligned to the 24 Clusters of Orthologous Groups classifications and 15,330 unigenes were classified into 266 Kyoto Encyclopedia of Genes and Genomes pathways. Based on the newly assembled transcriptome, the DTA profiling approach was applied to investigate the differentially expressed genes related to shading stress. A total of 3.6 million and 3.5 million high-quality tags were generated from shaded and non-shaded libraries, respectively. As many as 1,039 unigenes were shown to be significantly differentially regulated. Eleven of the 14 differentially regulated unigenes, which were randomly selected for more detailed expression comparison during the course of shading treatment, were identified as being likely to be involved in the process of fruitlet abscission in litchi.

CONCLUSIONS

The assembled transcriptome of litchi fruit provides a global description of expressed genes in litchi fruit development, and could serve as an ideal repository for future functional characterization of specific genes. The DTA analysis revealed that more than 1000 differentially regulated unigenes respond to the shading signal, some of which might be involved in the fruitlet abscission process in litchi, shedding new light on the molecular mechanisms underlying organ abscission.

De novo assembly and characterization of the transcriptome, and development of SSR markers in wax gourd (Benicasa hispida)

BACKGROUND

Wax gourd is a widely used vegetable of Cucuribtaceae, and also has important medicinal and health values. However, the genomic resources of wax gourd were scarcity, and only a few nucleotide sequences could be obtained in public databases.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we examined transcriptome in wax gourd. More than 44 million of high quality reads were generated from five different tissues of wax gourd using Illumina paired-end sequencing technology. Approximately 4 Gbp data were generated, and de novo assembled into 65,059 unigenes, with an N50 of 1,132 bp. Based on sequence similarity search with known protein database, 36,070 (55.4%) showed significant similarity to known proteins in Nr database, and 24,969 (38.4%) had BLAST hits in Swiss-Prot database. Among the annotated unigenes, 14,994 of wax gourd unigenes were assigned to GO term annotation, and 23,977 were found to have COG classifications. In addition, a total of 18,713 unigenes were assigned to 281 KEGG pathways. Furthermore, 6,242 microsatellites (simple sequence repeats) were detected as potential molecular markers in wax gourd. Two hundred primer pairs for SSRs were designed for validation of the amplification and polymorphism. The result showed that 170 of the 200 primer pairs were successfully amplified and 49 (28.8%) of them exhibited polymorphisms.

CONCLUSION/SIGNIFICANCE

Our study enriches the genomic resources of wax gourd and provides powerful information for future studies. The availability of this ample amount of information about the transcriptome and SSRs in wax gourd could serve as valuable basis for studies on the physiology, biochemistry, molecular genetics and molecular breeding of this important vegetable crop.