De novo sequencing of Hypericum perforatum transcriptome to identify potential genes involved in the biosynthesis of active metabolites

Single molecule real-time sequencing data sets of Hypericum perforatum L. plantlets and cell suspension cultures

Hypericum is a large genus that includes more than 500 species of pharmacological, ecological and conservation value. Although latest advances in sequencing technologies were extremely exploited for generating and assembling genomes of many living organisms, annotated whole genome sequence data is not publicly available for any of the Hypericum species so far. Bioavailability of secondary metabolites varies for different tissues and the data derived from different cultures will be a valuable tool for comparative studies. Here, we report the single molecule real-time sequencing (SMRT) data sets of Hypericum perforatum L. plantlets and cell suspension cultures for the first time. Sequencing data from cell suspension cultures yielded more than 33,000 high-quality transcripts from 20 Gb of raw data, while more than 55,000 high-quality transcripts were obtained from 35 Gb of raw data from plantlets. This dataset is a valuable tool for comparative transcriptomic analysis and will help to understand the unknown biosynthetic pathways of high medicinal value in the Hypericum genus.

Transcriptional responses of Hypericum perforatum cells to Agrobacterium tumefaciens and differential gene expression in dark glands

Hypericum perforatum L. (St. John's wort) is a well-known medicinal plant that possesses secondary metabolites with beneficial pharmacological properties. However, improvement in the production of secondary metabolites via genetic manipulation is a challenging task as H. perforatum remains recalcitrant to Agrobacterium tumefaciens-mediated transformation. Here, the transcripts of key genes involved in several plant defence responses (secondary metabolites, RNA silencing, reactive oxygen species (ROS) and specific defence genes) were investigated in H. perforatum suspension cells inoculated with A. tumefaciens by quantitative real-time PCR. Results indicated that key genes from the xanthone, hypericin and melatonin biosynthesis pathways, the ROS-detoxification enzyme HpAOX, as well as the defence genes Hyp-1 and HpPGIP, were all upregulated to rapidly respond to A. tumefaciens elicitation in H. perforatum. By contrast, expression levels of genes involved in hyperforin and flavonoid biosynthesis pathways were markedly downregulated upon A. tumefaciens elicitation. In addition, we compared the expression patterns of key genes in H. perforatum leaf tissues with and without dark glands, a major site of secondary metabolite production. Overall, we provide evidence for the upregulation of several phenylpropanoid pathway genes in response to elicitation by Agrobacterium, suggesting that production of secondary metabolites could modulate H. perforatum recalcitrance to A. tumefaciens-mediated transformation.

Whole-genome sequence data of Hypericum perforatum and functional characterization of melatonin biosynthesis by N-acetylserotonin O-methyltransferase

Hypericum perforatum is among the most commonly used herbal remedies and supplements. The aerial plant parts are often used to treat depression. Due to the lack of genomic information of H. perforatum, the gene networks regulating secondary metabolite synthesis remain unclear. Here, we present a high-quality genome for H. perforatum with a 2.3-Mb scaffold N50. The draft assembly covers 91.9% of the predicted genome and represents the fourth sequenced genus in the order Malpighiales. Comparing this sequence with model or related species revealed that Populus trichocarpa and Hevea brasiliensis could be grouped into one branch, while H. perforatum and Linum usitatissimum are grouped in another branch. Combined with transcriptome data, 40 key genes related to melatonin, hyperforin, and hypericin synthesis were screened and analyzed. Five N-acetylserotonin O-methyltransferases (HpASMT1-HpASMT5) were cloned and functionally characterized. Purified HpASMT3 protein converted N-acetylserotonin into melatonin with a V_max of about 1.35 pkat/mg protein. HpASMT1 and HpASMT3 overexpression in Arabidopsis mutants caused 1.5-2-fold higher melatonin content than in mutant and wild-type plants. The endogenous reactive oxygen species (ROS) in transgenic plants was significantly lower than ROS in mutant and wild-type plants, suggesting higher drought tolerance. The obtained genomic data offer new resources for further study on the evolution of Hypericaceae family, but also provide a basis for further study of melatonin biosynthetic pathways in other plants.

Comparative transcriptome among Euscaphis konishii Hayata tissues and analysis of genes involved in flavonoid biosynthesis and accumulation

BACHGROUND

Euscaphis konishii Hayata, a member of the Staphyleaceae Family, is a plant that has been widely used in Traditional Chinese Medicine and it has been the source for several types of flavonoids. To identify candidate genes involved in flavonoid biosynthesis and accumulation, we analyzed transcriptome data from three E. konishii tissues (leaf, branch and capsule) using Illumina Hiseq 2000 platform.

RESULTS

A total of 91.7, 100.3 and 100.1million clean reads were acquired for the leaf, branch and capsule, respectively; and 85,342 unigenes with a mean length of 893.60 bp and N50 length of 1307 nt were assembled using Trinity program. BLASTx analysis allowed to annotate 40,218 unigenes using public protein databases, including NR, KOG/COG/eggNOG, Swiss-Prot, KEGG and GO. A total of 14,291 (16.75%) unigenes were assigned to 128 KEGG pathways, and 900 unigenes were annotated into 22 KEGG secondary metabolites, including flavonoid biosynthesis. The structure enzymes involved in flavonoid biosynthesis, such as phenylalanine ammonia lyase, cinnamate 4-hydroxylase, 4-coumarate CoA ligase, shikimate O-hydroxycinnamoyltransferase, coumaroylquinate 3'-monooxygenase, caffeoyl-CoA O-methyltransferase, chalcone synthase, chalcone isomerase, flavanone 3-hydroxylase, flavonoid 3'-hydroxylase, flavonoid 3',5'-hydroxylase, flavonolsynthese, dihydroflavonol 4-reductase, anthocyanidinreductase, leucoanthocyanidin dioxygenase, leucoanthocyanidin reductase, were identified in the transcriptome data, 40 UDP-glycosyltransferase (UGT), 122 Cytochrome P450 (CYP) and 25 O-methyltransferase (OMT) unigenes were also found. A total of 295 unigenes involved in flavonoid transport and 220 transcription factors (97 MYB, 84 bHLH and 39 WD40) were identified. Furthermore, their expression patterns among different tissues were analyzed by DESeq, the differentially expressed genes may play important roles in tissues-specific synthesis, accumulation and modification of flavonoids.

CONCLUSION

We present here the de novo transcriptome analysis of E. konishii and the identification of candidate genes involved in biosynthesis and accumulation of flavonoid. In general, these results are an important resource for further research on gene expression, genomic and functional genomics in E. konishii and other related species.

Comparative Transcriptome Analysis Identifies Genes Involved in Diosgenin Biosynthesis in Trigonella foenum-graecum L

Trigonella foenum-graecum L. (fenugreek) is a valuable resource of producing diosgenin which serves as a substrate for synthesizing more than two hundred kinds of steroidal drugs. Phytochemical analysis indicated that methyl jasmonate (MeJA) efficiently induced diosgenin biosynthesis in fenugreek seedlings. Though early steps up to cholesterol have recently been elucidated in plants, cytochrome P450 (CYP)- and glycosyltransferase (GT)-encoding genes involved in the late steps from cholesterol to diosgenin remain unknown. This study established comparative fenugreek transcriptome datasets from the MeJA-treated seedlings and the corresponding control lines. Differential gene expression analysis identified a number of MeJA-induced CYP and GT candidate genes. Further gene expression pattern analysis across a different MeJA-treating time points, together with a phylogenetic analysis, suggested specific family members of CYPs and GTs that may participate in the late steps during diosgenin biosynthesis. MeJA-induced transcription factors (TFs) that may play regulatory roles in diosgenin biosynthesis were also discussed. This study provided a valuable genetic resource to functionally characterize the genes involved in diosgenin biosynthesis, which will push forward the production of diosgenin in microbial organisms using a promising synthetic biology strategy.

Characterization of the Lycium barbarum fruit transcriptome and development of EST-SSR markers

Lycium barbarum, commonly known as goji, is important in Chinese herbal medicine and its fruit is a very important agricultural and biological product. However, the molecular mechanism of formation of its fruit and associated medicinal and nutritional components is unexplored. Moreover, this species lacks SSR markers due to lack of genomic and transcriptomic information. In this study, a total of 139,333 unigenes with average length of 1049 bp and N50 of 1579 bp are obtained by trinity assembly from Illumina sequencing reads. A total of 92,498 (66.38%) unigenes showed similarities in at least one database including Nr (46.15%), Nt (56.56%), KO (15.56%), Swiss-prot (33.34%), Pfam (33.43%), GO (33.62%) and KOG/COG (17.55%). Genes in flavonoid and taurine biosynthesis pathways were found and validated by RT-qPCR. A total of 50,093 EST-SSRs were identified from 38,922 unigenes, and 22,537 EST-SSR primer pairs were designed. Four hundred pairs of SSR markers were randomly selected to validate assembly quality, of which 352 (88%) were successful in PCR amplification of genomic DNA from 11 Lycium accessions and 210 produced polymorphisms. The polymorphic loci showed that the genetic similarity of the 11 Lycium accessions ranged from 0.50 to 0.99 and the accessions could be divided into 4 groups. These results will facilitate investigations of the molecular mechanism of formation of L. barbarum fruit and associated medicinal and nutritional components, and will be of value to novel gene discovery and functional genomic studies. The EST-SSR markers will be useful for genetic diversity evaluation, genetic mapping and marker-assisted breeding.

Transcriptional Responses of Chilean Quinoa (Chenopodium quinoa Willd.) Under Water Deficit Conditions Uncovers ABA-Independent Expression Patterns

HIGHLIGHTS R49 genotype displayed best performance on selected physiological parameters and highest tolerance to drought.R49 drought over-represented transcripts has exhibited 19% of genes (306 contigs) that presented no homology to published databases.Expression pattern for canonical responses to drought such as ABA biosynthesis and other genes induced in response to drought were assessed by qPCR. Global freshwater shortage is one of the biggest challenges of our time, often associated to misuse, increased consumption demands and the effects of climate change, paralleled with the desertification of vast areas. Chenopodium quinoa (Willd.) represents a very promising species, due to both nutritional content and cultivation under water constraint. We characterized drought tolerance of three Chilean genotypes and selected Genotype R49 (Salares ecotype) based upon Relative Water Content (RWC), Electrolyte Leakage (EL) and maximum efficiency of photosystem II (Fv/Fm) after drought treatment, when compared to another two genotypes. Exploratory RNA-Seq of R49 was generated by Illumina paired-ends method comparing drought and control irrigation conditions. We obtained 104.8 million reads, with 54 million reads for control condition and 51 million reads for drought condition. Reads were assembled in 150,952 contigs, were 31,523 contigs have a reading frame of at least 300 nucleotides (100 aminoacids). BLAST2GO annotation showed a 15% of genes without homology to NCBI proteins, but increased to 19% (306 contigs) when focused into drought-induced genes. Expression pattern for canonical drought responses such as ABA biosynthesis and other genes induced were assessed by qPCR, suggesting novelty of R49 drought responses.

Quality control of Hypericum perforatum L. analytical challenges and recent progress

Objectives

The most widely applied qualitative and quantitative analytical methods in the quality control of Hypericum perforatum extracts will be reviewed, including routine analytical tools and most modern approaches.

Key findings

Biologically active components of H. perforatum are chemically diverse; therefore, different chromatographic and detection methods are required for the comprehensive analysis of St. John's wort extracts. Naphthodianthrones, phloroglucinols and flavonoids are the most widely analysed metabolites of this plant. For routine quality control, detection of major compounds belonging to these groups seems to be sufficient; however, closer characterization requires the detection of minor compounds as well.

Conclusions

TLC and HPTLC are basic methods in the routine analysis, whereas HPLC-DAD is the most widely applied method for quantitative analysis due to its versatility. LC-MS is gaining importance in pharmacokinetic studies due to its sensitivity. Modern approaches, such as DNA barcoding, NIRS and NMR metabolomics, may offer new possibilities for the more detailed characterization of secondary metabolite profile of H. perforatum extracts.

Pistil Transcriptome Analysis to Disclose Genes and Gene Products Related to Aposporous Apomixis in Hypericum perforatum L

Unlike sexual reproduction, apomixis encompasses a number of reproductive strategies, which permit maternal genome inheritance without genetic recombination and syngamy. The key biological features of apomixis are the circumvention of meiosis (i.e., apomeiosis), the differentiation of unreduced embryo sacs and egg cells, and their autonomous development in functional embryos through parthenogenesis, and the formation of viable endosperm either via fertilization-independent means or following fertilization with a sperm cell. Despite the importance of apomixis for breeding of crop plants and although much research has been conducted to study this process, the genetic control of apomixis is still not well understood. Hypericum perforatum is becoming an attractive model system for the study of aposporous apomixis. Here we report results from a global gene expression analysis of H. perforatum pistils collected from sexual and aposporous plant accessions for the purpose of identifying genes, biological processes and molecular functions associated with the aposporous apomixis pathway. Across two developmental stages corresponding to the expression of aposporous apomeiosis and parthenogenesis in ovules, a total of 224 and 973 unigenes were found to be significantly up- and down-regulated with a fold change ≥ 2 in at least one comparison, respectively. Differentially expressed genes were enriched for multiple gene ontology (GO) terms, including cell cycle, DNA metabolic process, and single-organism cellular process. For molecular functions, the highest scores were recorded for GO terms associated with DNA binding, DNA (cytosine-5-)-methyltransferase activity and heterocyclic compound binding. As deregulation of single components of the sexual developmental pathway is believed to be a trigger of the apomictic reproductive program, all genes involved in sporogenesis, gametogenesis and response to hormonal stimuli were analyzed in great detail. Overall, our data suggest that phenotypic expression of apospory is concomitant with the modulation of key genes involved in the sexual reproductive pathway. Furthermore, based on gene annotation and co-expression, we underline a putative role of hormones and key actors playing in the RNA-directed DNA methylation pathway in regulating the developmental changes occurring during aposporous apomixis in H. perforatum.

Transcriptome analysis revealed the dynamic oil accumulation in Symplocos paniculata fruit

BACKGROUND

Symplocos paniculata, asiatic sweetleaf or sapphire berry, is a widespread shrub or small tree from Symplocaceae with high oil content and excellent fatty acid composition in fruit. It has been used as feedstocks for biodiesel and cooking oil production in China. Little transcriptome information is available on the regulatory molecular mechanism of oil accumulation at different fruit development stages.

RESULTS

The transcriptome at four different stages of fruit development (10, 80,140, and 170 days after flowering) of S. paniculata were analyzed. Approximately 28 million high quality clean reads were generated. These reads were trimmed and assembled into 182,904 non-redundant putative transcripts with a mean length of 592.91 bp and N50 length of 785 bp, respectively. Based on the functional annotation through Basic Local Alignment Search Tool (BLAST) with public protein database, the key enzymes involved in lipid metabolism were identified, and a schematic diagram of the pathway and temporal expression patterns of lipid metabolism was established. About 13,939 differentially expressed unigenes (DEGs) were screened out using differentially expressed sequencing (DESeq) method. The transcriptional regulatory patterns of the identified enzymes were highly related to the dynamic oil accumulation along with the fruit development of S. paniculata. In addition, quantitative real-time PCR (qRT-PCR) of six vital genes was significantly correlated with DESeq data.

CONCLUSIONS

The transcriptome sequences obtained and deposited in NCBI would enrich the public database and provide an unprecedented resource for the discovery of the genes associated with lipid metabolism pathway in S. paniculata. Results in this study will lay the foundation for exploring transcriptional regulatory profiles, elucidating molecular regulatory mechanisms, and accelerating genetic engineering process to improve the yield and quality of seed oil of S. paniculata.

Comparative Transcriptome Reconstruction of Four Hypericum Species Focused on Hypericin Biosynthesis

Next generation sequencing technology rapidly developed research applications in the field of plant functional genomics. Several Hypericum spp. with an aim to generate and enhance gene annotations especially for genes coding the enzymes supposedly included in biosynthesis of valuable bioactive compounds were analyzed. The first de novo transcriptome profiling of Hypericum annulatum Moris, H. tomentosum L., H. kalmianum L., and H. androsaemum L. leaves cultivated in vitro was accomplished. All four species with only limited genomic information were selected on the basis of differences in ability to synthesize hypericins and presence of dark nodules accumulating these metabolites with purpose to enrich genomic background of Hypericum spp. H. annulatum was chosen because of high number of the dark nodules and high content of hypericin. H. tomentosum leaves are typical for the presence of only 1-2 dark nodules localized in the apical part. Both H. kalmianum and H. androsaemum lack hypericin and have no dark nodules. Four separated datasets of the pair-end reads were gathered and used for de novo assembly by Trinity program. Assembled transcriptomes were annotated to the public databases Swiss-Prot and non-redundant protein database (NCBI-nr). Gene ontology analysis was performed. Differences of expression levels in the marginal tissues with dark nodules and inner part of leaves lacking these nodules indicate a potential genetic background for hypericin formation as the presumed site of hypericin biosynthesis is in the cells adjacent to these structures. Altogether 165 contigs in H. annulatum and 100 contigs in H. tomentosum were detected as significantly differentially expressed (P < 0.05) and upregulated in the leaf rim tissues containing the dark nodules. The new sequences homologous to octaketide synthase and enzymes catalyzing phenolic oxidative coupling reactions indispensable for hypericin biosynthesis were discovered. The presented transcriptomic sequence data will improve current knowledge about the selected Hypericum spp. with proposed relation to hypericin biosynthesis and will provide a useful resource of genomic information for consequential studies in the field of functional genomics, proteomics and metabolomics.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

A Perspective on Hypericum perforatum Genetic Transformation

Hypericum perforatum (St John's wort) is a reservoir of diverse classes of biologically active and high value secondary metabolites, which captured the interest of both researchers and the pharmaceutical industry alike. Several studies and clinical trials have shown that H. perforatum extracts possess an astounding array of pharmacological properties. These properties include antidepressant, anti-inflammatory, antiviral, anti-cancer, and antibacterial activities; and are largely attributed to the naphtodianthrones and xanthones found in the genus. Hence, improving their production via genetic manipulation is an important strategy. In spite of the presence of contemporary genome editing tools, genetic improvement of this genus remains challenging without robust transformation methods in place. In the recent past, we found that H. perforatum remains recalcitrant to Agrobacterium tumefaciens mediated transformation partly due to the induction of plant defense responses coming into play. However, H. perforatum transformation is possible via a non-biological method, biolistic bombardment. Some research groups have observed the induction of hairy roots in H. perforatum after Agrobacterium rhizogenes co-cultivation. In this review, we aim at updating the available methods for regeneration and transformation of H. perforatum. In addition, we also propose a brief perspective on certain novel strategies to improve transformation efficiency in order to meet the demands of the pharmaceutical industry via metabolic engineering.

Melatonin in Plants - Diversity of Levels and Multiplicity of Functions

Melatonin has been detected in numerous plant species. A particularly surprising finding concerns the highly divergent levels of melatonin that vary between species, organs and environmental conditions, from a few pg/g to over 20 μg/g, reportedly up to 200 μg/g. Highest values have been determined in oily seeds and in plant organs exposed to high UV radiation. The divergency of melatonin concentrations is discussed under various functional aspects and focused on several open questions. This comprises differences in precursor availability, catabolism, the relative contribution of isoenzymes of the melatonin biosynthetic pathway, and differences in rate limitation by either serotonin N-acetyltransferase or N-acetylserotonin O-methyltransferase. Other differences are related to the remarkable pleiotropy of melatonin, which exhibits properties as a growth regulator and morphogenetic factor, actually debated in terms of auxin-like effects, and as a signaling molecule that modulates pathways of ethylene, abscisic, jasmonic and salicylic acids and is involved in stress tolerance, pathogen defense and delay of senescence. In the context of high light/UV intensities, elevated melatonin levels exceed those required for signaling via stress-related phytohormones and may comprise direct antioxidant and photoprotectant properties, perhaps with a contribution of its oxidatively formed metabolites, such as N (1)-acetyl-N (2)-formyl-5-methoxykynuramine and its secondary products. High melatonin levels in seeds may also serve antioxidative protection and have been shown to promote seed viability and germination capacity.

De Novo Assembly and Transcriptome Characterization of Canine Retina Using High-Throughput Sequencing

We performed transcriptome sequencing of canine retinal tissue by 454 GS-FLX and Ion Torrent PGM platforms. RNA-Seq analysis by CLC Genomics Workbench mapped expression of 10,360 genes. Gene ontology analysis of retinal transcriptome revealed abundance of transcripts known to be involved in vision associated processes. The de novo assembly of the sequences using CAP3 generated 29,683 contigs with mean length of 560.9 and N50 of 619 bases. Further analysis of contigs predicted 3,827 full-length cDNAs and 29,481 (99%) open reading frames (ORFs). In addition, 3,782 contigs were assigned to 316 KEGG pathways which included melanogenesis, phototransduction, and retinol metabolism with 33, 15, and 11 contigs, respectively. Among the identified microsatellites, dinucleotide repeats were 68.84%, followed by trinucleotides, tetranucleotides, pentanucleotides, and hexanucleotides in proportions of 25.76, 9.40, 2.52, and 0.96%, respectively. This study will serve as a valuable resource for understanding the biology and function of canine retina.

Removal of redundant contigs from de novo RNA-Seq assemblies via homology search improves accurate detection of differentially expressed genes

Background

For plant species with unsequenced genomes, cDNA contigs created by de novo assembly of RNA-Seq reads are used as reference sequences for comparative analysis of RNA-Seq datasets and the detection of differentially expressed genes (DEGs). Redundancies in such contigs are evident in previous RNA-Seq studies, and such redundancies can lead to difficulties in subsequent analysis. Nevertheless, the effects of removing redundancy from contig assemblies on comparative RNA-Seq analysis have not been evaluated.

Results

Here we describe a method for removing redundancy from raw contigs that were primarily created by de novo assembly of Arabidopsis thaliana RNA-Seq reads. Specifically, the contigs with the highest bit scores were selected from raw contigs by a homology search against the gene dataset in the TAIR10 database. The two existing methods for removal of redundancy based on contig length or clustering analysis used to eliminate redundancies from raw contigs. Contig number was reduced most effectively with the method based on homology search. In a comparative analysis of RNA-Seq datasets, DEGs detected in contigs that underwent redundancy removal via the homology search method showed the highest identity to the DEGs detected when the TAIR10 gene dataset was used as an exact reference. Redundancy in raw contigs could also be removed by a homology search against integrated protein datasets from several plant species other than A. thaliana. DEGs detected using contigs that underwent such redundancy-removed also showed high homology to DEGs detected using the TAIR10 gene dataset.

Conclusion

Here we describe a method for removing redundant contigs within raw contigs; this method involves a homology search against a gene or protein database. In principal, this method can be used with unsequenced plant genomes that lack a well-developed gene database. Redundant contigs were not removed adequately via either of two existing methods, but our method allowed for removal of all redundant contigs. To our knowledge, this is the first reported improvement in accurate detection of DEGs via comparative RNA-Seq analysis that involved preparation of a non-redundant reference sequence. This method could be used to rapidly and cost-effectively detect useful genes in unsequenced plants.

Andrographis paniculata transcriptome provides molecular insights into tissue-specific accumulation of medicinal diterpenes

Background

Kalmegh (Andrographis paniculata) has been widely exploited in traditional medicine for the treatment of infectious diseases and health disorders. Ent-labdane-related diterpene (ent-LRD) specialized (i.e., secondary) metabolites of kalmegh such as andrographolide, neoandrographolide and 14-deoxy-11,12-didehydroandrographolide, are known for variety of pharmacological activities. However, due to the lack of genomic and transcriptomic information, underlying molecular basis of ent-LRDs biosynthesis has remained largely unknown. To identify candidate genes of the ent-LRD biosynthetic pathway, we performed comparative transcriptome analysis using leaf and root tissues that differentially accumulate ent-LRDs.

Results

De novo assembly of Illumina HiSeq2000 platform-generated paired-end sequencing reads resulted into 69,011 leaf and 64,244 root transcripts which were assembled into a total of 84,628 unique transcripts. Annotation of these transcripts to the Uniprot, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Carbohydrate-Active Enzymes (CAZy) databases identified candidate transcripts of the ent-LRD biosynthetic pathway. These included transcripts that encode enzymes of the plastidial 2C-methyl-D-erythritol-4-phosphate pathway which provides C5 isoprenoid precursors for the ent-LRDs biosynthesis, geranylgeranyl diphosphate synthase, class II diterpene synthase (diTPS), cytochrome P450 monooxygenase and glycosyltransferase. Three class II diTPSs (ApCPS1, ApCPS2 and ApCPS3) that showed distinct tissue-specific expression profiles and are phylogenetically related to the dicotyledon ent-copalyl diphosphate synthases, are identified. ApCPS1, ApCPS2 and ApCPS3 encode for 832-, 817- and 797- amino acids proteins of 55–63 % identity, respectively. Spatio-temporal patterns of transcripts and ent-LRDs accumulation are consistent with the involvement of ApCPS1 in general (i.e., primary) metabolism for the biosynthesis of phytohormone gibberellin, ApCPS2 in leaf specialized ent-LRDs biosynthesis and ApCPS3 in root diterpene biosynthesis. Moreover, simple sequence repeats (SSRs) that might assist in genotyping and developing specific chemotypes were identified in transcripts of the specialized metabolic pathways, including ent-LRDs.

Conclusions

Comparative analysis of root and leaf transcriptomes disclosed novel genes of the ent-LRD biosynthetic pathway, including three class II diTPSs that showed discrete spatio-temporal expression patterns; thus, suggesting their participation into distinct diterpene metabolic pathways of kalmegh. Overall, these results will be useful in understanding molecular basis of the medicinal ent-LRDs biosynthesis and developing breeding strategies for improving their yields.

Electronic supplementary material

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

Next Generation Sequencing and Transcriptome Analysis Predicts Biosynthetic Pathway of Sennosides from Senna (Cassia angustifolia Vahl.), a Non-Model Plant with Potent Laxative Properties

Senna (Cassia angustifolia Vahl.) is a world's natural laxative medicinal plant. Laxative properties are due to sennosides (anthraquinone glycosides) natural products. However, little genetic information is available for this species, especially concerning the biosynthetic pathways of sennosides. We present here the transcriptome sequencing of young and mature leaf tissue of Cassia angustifolia using Illumina MiSeq platform that resulted in a total of 6.34 Gb of raw nucleotide sequence. The sequence assembly resulted in 42230 and 37174 transcripts with an average length of 1119 bp and 1467 bp for young and mature leaf, respectively. The transcripts were annotated using NCBI BLAST with 'green plant database (txid 33090)', Swiss Prot, Kyoto Encylcopedia of Genes & Genomes (KEGG), Cluster of Orthologous Gene (COG) and Gene Ontology (GO). Out of the total transcripts, 40138 (95.0%) and 36349 (97.7%) from young and mature leaf, respectively, were annotated by BLASTX against green plant database of NCBI. We used InterProscan to see protein similarity at domain level, a total of 34031 (young leaf) and 32077 (mature leaf) transcripts were annotated against the Pfam domains. All transcripts from young and mature leaf were assigned to 191 KEGG pathways. There were 166 and 159 CDS, respectively, from young and mature leaf involved in metabolism of terpenoids and polyketides. Many CDS encoding enzymes leading to biosynthesis of sennosides were identified. A total of 10,763 CDS differentially expressing in both young and mature leaf libraries of which 2,343 (21.7%) CDS were up-regulated in young compared to mature leaf. Several differentially expressed genes found functionally associated with sennoside biosynthesis. CDS encoding for many CYPs and TF families were identified having probable roles in metabolism of primary as well as secondary metabolites. We developed SSR markers for molecular breeding of senna. We have identified a set of putative genes involved in various secondary metabolite pathways, especially those related to the synthesis of sennosides which will serve as an important platform for public information about gene expression, genomics, and functional genomics in senna.

De novo sequencing of the Hypericum perforatum L. flower transcriptome to identify potential genes that are related to plant reproduction sensu lato

Background

St. John’s wort (Hypericum perforatum L.) is a medicinal plant that produces important metabolites with antidepressant and anticancer activities. Recently gained biological information has shown that this species is also an attractive model system for the study of a naturally occurring form of asexual reproduction called apomixis, which allows cloning plants through seeds. In aposporic gametogenesis, one or multiple somatic cells belonging to the ovule nucellus change their fate by dividing mitotically and developing functionally unreduced embryo sacs by mimicking sexual gametogenesis. Although the introduction of apomixis into agronomically important crops could have revolutionary implications for plant breeding, the genetic control of this mechanism of seed formation is still not well understood for most of the model species investigated so far. We used Roche 454 technology to sequence the entire H. perforatum flower transcriptome of whole flower buds and single flower verticils collected from obligately sexual and unrelated highly or facultatively apomictic genotypes, which enabled us to identify RNAs that are likely exclusive to flower organs (i.e., sepals, petals, stamens and carpels) or reproductive strategies (i.e., sexual vs. apomictic).

Results

Here we sequenced and annotated the flower transcriptome of H. perforatum with particular reference to reproductive organs and processes. In particular, in our study we characterized approximately 37,000 transcripts found expressed in male and/or female reproductive organs, including tissues or cells of sexual and apomictic flower buds. Ontological annotation was applied to identify major biological processes and molecular functions involved in flower development and plant reproduction. Starting from this dataset, we were able to recover and annotate a large number of transcripts related to meiosis, gametophyte/gamete formation, and embryogenesis, as well as genes that are exclusively or preferentially expressed in sexual or apomictic libraries. Real-Time RT-qPCR assays on pistils and anthers collected at different developmental stages from accessions showing alternative modes of reproduction were used to identify potential genes that are related to plant reproduction sensu lato in H. perforatum.

Conclusions

Our approach of sequencing flowers from two fully obligate sexual genotypes and two unrelated highly apomictic genotypes, in addition to different flower parts dissected from a facultatively apomictic accession, enabled us to analyze the complexity of the flower transcriptome according to its main reproductive organs as well as for alternative reproductive behaviors. Both annotation and expression data provided original results supporting the hypothesis that apomixis in H. perforatum relies upon spatial or temporal mis-expression of genes acting during female sexual reproduction. The present analyses aim to pave the way toward a better understanding of the molecular basis of flower development and plant reproduction, by identifying genes or RNAs that may differentiate or regulate the sexual and apomictic reproductive pathways in H. perforatum.

Electronic supplementary material

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

Transcriptome sequencing of Codonopsis pilosula and identification of candidate genes involved in polysaccharide biosynthesis

Background

Codonopsis pilosula (Franch.) Nannf. is one of the most widely used medicinal plants. Although chemical and pharmacological studies have shown that codonopsis polysaccharides (CPPs) are bioactive compounds and that their composition is variable, their biosynthetic pathways remain largely unknown. Next-generation sequencing is an efficient and high-throughput technique that allows the identification of candidate genes involved in secondary metabolism.

Principal Findings

To identify the components involved in CPP biosynthesis, a transcriptome library, prepared using root and other tissues, was assembled with the help of Illumina sequencing. A total of 9.2 Gb of clean nucleotides was obtained comprising 91,175,044 clean reads, 102,125 contigs, and 45,511 unigenes. After aligning the sequences to the public protein databases, 76.1% of the unigenes were annotated. Among these annotated unigenes, 26,189 were assigned to Gene Ontology categories, 11,415 to Clusters of Orthologous Groups, and 18,848 to Kyoto Encyclopedia of Genes and Genomes pathways. Analysis of abundance of transcripts in the library showed that genes, including those encoding metallothionein, aquaporin, and cysteine protease that are related to stress responses, were in the top list. Among genes involved in the biosynthesis of CPP, those responsible for the synthesis of UDP-L-arabinose and UDP-xylose were highly expressed.

Significance

To our knowledge, this is the first study to provide a public transcriptome dataset prepared from C. pilosula and an outline of the biosynthetic pathway of polysaccharides in a medicinal plant. Identified candidate genes involved in CPP biosynthesis provide understanding of the biosynthesis and regulation of CPP at the molecular level.

Reference genes selection and normalization of oxidative stress responsive genes upon different temperature stress conditions in Hypericum perforatum L

Reverse transcription-quantitative real-time PCR (RT-qPCR) is a widely used technique for gene expression analysis. The reliability of this method depends largely on the suitable selection of stable reference genes for accurate data normalization. Hypericum perforatum L. (St. John's wort) is a field growing plant that is frequently exposed to a variety of adverse environmental stresses that can negatively affect its productivity. This widely known medicinal plant with broad pharmacological properties (anti-depressant, anti-tumor, anti-inflammatory, antiviral, antioxidant, anti-cancer, and antibacterial) has been overlooked with respect to the identification of reference genes suitable for RT-qPCR data normalization. In this study, 11 candidate reference genes were analyzed in H. perforatum plants subjected to cold and heat stresses. The expression stability of these genes was assessed using GeNorm, NormFinder and BestKeeper algorithms. The results revealed that the ranking of stability among the three algorithms showed only minor differences within each treatment. The best-ranked reference genes differed between cold- and heat-treated samples; nevertheless, TUB was the most stable gene in both experimental conditions. GSA and GAPDH were found to be reliable reference genes in cold-treated samples, while GAPDH showed low expression stability in heat-treated samples. 26SrRNA and H2A had the highest stabilities in the heat assay, whereas H2A was less stable in the cold assay. Finally, AOX1, AOX2, CAT1 and CHS genes, associated with plant stress responses and oxidative stress, were used as target genes to validate the reliability of identified reference genes. These target genes showed differential expression profiles over time in treated samples. This study not only is the first systematic analysis for the selection of suitable reference genes for RT-qPCR studies in H. perforatum subjected to temperature stress conditions, but may also provide valuable information about the roles of genes associated with temperature stress responses.

Hypericin biosynthesis in Hypericum hookerianum Wight and Arn: investigation on biochemical pathways using metabolite inhibitors and suppression subtractive hybridization

The biochemical pathway to hypericin biosynthesis is presumed to be polyketide synthase (PKS) mediated, but it has not been experimentally validated, and no alternate route (chorismate/o-succinylbenzoate pathway) has been analyzed. We report here our earlier developed auxin inducible culture systems of Hypericum hookerianum as a model, to study the metabolic pathway to hypericin synthesis. Inhibitors of the alternate pathway at varying concentrations showed steady synthesis of total hypericins with means of 2.80±0.22, 18.75±0.01; 16.39±3.75, 29.60±1.90 (mevinolin) 2.53±0.10, 18.12±0.56; 0.14±0.01, 14.28±1.11 (fosmidomycin) and 2.7±0.35, 18.75±0.61; 0.14±0.01, 12.80±1.09 mg g(-1) DW (glyphosate) in the control and auxin-induced shoot and shoot-forming callus cultures, respectively. SSH analysis classified the differentially expressed sequences into protein synthesis (38%), modification (20%), electron transport (9%) and remaining as unclassified (11%) and unknown proteins (22%). Functional annotation of sequences indicates the presence of additional protein components besides PKS activity. Our results demonstrate direct biochemical and molecular evidence of PKS hypothesis of hypericin biosynthesis for the first time.

An insight into the genes involved in secoiridoid biosynthesis in Gentiana macrophylla by RNA-seq

The dried root of Gentiana macrophylla is a well-known traditional Chinese herbal medicine for treating jaundice, hepatitis, and stomachic and choleretic ailments. However, natural sources are now in short supply. A lack of information about its genetic background has been a great hindrance to producing its active constituents via genetic engineering. We performed RNA-seq to obtain 42,918 unigenes (average length = 667 bp) in its transcriptome. Of these, 32,141 (74.89 %) were annotated and 2,339 unigenes were assigned to secondary-metabolite pathways. In all, 114 putative unigenes involved in secoiridoid biosynthesis were identified in our transcriptome library. A Blast X search against the Arabidopsis gene regulatory information server showed that 4,413 unigenes are homologous to transcription factor genes from Arabidopsis. Organ-specific genes and candidate gene expression profiles were also investigated with digital gene expression technology. Quantitative PCR was used to verify the expression patterns of several novel transcripts involved in secoiridoid biosynthesis. Our results not only enrich the gene resource but will also benefit research into the molecular genetics and functional genomics of this species.

Characterization of Liaoning cashmere goat transcriptome: sequencing, de novo assembly, functional annotation and comparative analysis

Background

Liaoning cashmere goat is a famous goat breed for cashmere wool. In order to increase the transcriptome data and accelerate genetic improvement for this breed, we performed denovo transcriptome sequencing to generate the first expressed sequence tag dataset for the Liaoning cashmere goat, using next-generation sequencing technology.

Results

Transcriptome sequencing of Liaoning cashmere goat on a Roche 454 platform yielded 804,601 high-quality reads. Clustering and assembly of these reads produced a non-redundant set of 117,854 unigenes, comprising 13,194 isotigs and 104,660 singletons. Based on similarity searches with known proteins, 17,356 unigenes were assigned to 6,700 GO categories, and the terms were summarized into three main GO categories and 59 sub-categories. 3,548 and 46,778 unigenes had significant similarity to existing sequences in the KEGG and COG databases, respectively. Comparative analysis revealed that 42,254 unigenes were aligned to 17,532 different sequences in NCBI non-redundant nucleotide databases. 97,236 (82.51%) unigenes were mapped to the 30 goat chromosomes. 35,551 (30.17%) unigenes were matched to 11,438 reported goat protein-coding genes. The remaining non-matched unigenes were further compared with cattle and human reference genes, 67 putative new goat genes were discovered. Additionally, 2,781 potential simple sequence repeats were initially identified from all unigenes.

Conclusion

The transcriptome of Liaoning cashmere goat was deep sequenced, denovo assembled, and annotated, providing abundant data to better understand the Liaoning cashmere goat transcriptome. The potential simple sequence repeats provide a material basis for future genetic linkage and quantitative trait loci analyses.

Transcriptome analysis of yellow horn (Xanthoceras sorbifolia Bunge): a potential oil-rich seed tree for biodiesel in China

Background

Yellow horn (Xanthoceras sorbifolia Bunge) is an oil-rich seed shrub that grows well in cold, barren environments and has great potential for biodiesel production in China. However, the limited genetic data means that little information about the key genes involved in oil biosynthesis is available, which limits further improvement of this species. In this study, we describe sequencing and de novo transcriptome assembly to produce the first comprehensive and integrated genomic resource for yellow horn and identify the pathways and key genes related to oil accumulation. In addition, potential molecular markers were identified and compiled.

Methodology/Principal Findings

Total RNA was isolated from 30 plants from two regions, including buds, leaves, flowers and seeds. Equal quantities of RNA from these tissues were pooled to construct a cDNA library for 454 pyrosequencing. A total of 1,147,624 high-quality reads with total and average lengths of 530.6 Mb and 462 bp, respectively, were generated. These reads were assembled into 51,867 unigenes, corresponding to a total of 36.1 Mb with a mean length, N50 and median of 696, 928 and 570 bp, respectively. Of the unigenes, 17,541 (33.82%) were unmatched in any public protein databases. We identified 281 unigenes that may be involved in de novo fatty acid (FA) and triacylglycerol (TAG) biosynthesis and metabolism. Furthermore, 6,707 SSRs, 16,925 SNPs and 6,201 InDels with high-confidence were also identified in this study.

Conclusions

This transcriptome represents a new functional genomics resource and a foundation for further studies on the metabolic engineering of yellow horn to increase oil content and modify oil composition. The potential molecular markers identified in this study provide a basis for polymorphism analysis of Xanthoceras, and even Sapindaceae; they will also accelerate the process of breeding new varieties with better agronomic characteristics.

Bioinformatics challenges in de novo transcriptome assembly using short read sequences in the absence of a reference genome sequence

Plant natural product research can be facilitated through genome and transcriptome sequencing approaches that generate informative sequence and expression datasets that enable characterization of biochemical pathways of interest. As the overwhelming majority of plant-derived natural products are derived from species with little, if any, sequence and/or genomic resources, the ability to perform whole genome shotgun sequencing and assembly has been and will continue to be transformative as access to a genome sequence provides molecular resources and a context for discovery and characterization of biosynthetic pathways. Due to the reduced size and complexity of the transcriptome relative to the genome, transcriptome sequencing provides a rapid, inexpensive approach to access gene sequences, gene expression abundances, and gene expression patterns in any species, including those that lack a reference genome sequence. To date, successful applications of RNA sequencing in conjunction with de novo transcriptome assembly has enabled identification of new genes in an array of biochemical pathways in plants. While sequencing technologies are well developed, challenges remain in the handling and analysis of transcriptome sequences. In this Highlight article, we provide an overview of the bioinformatics challenges associated with transcriptome analyses using short read sequences and how to address these issues in plant species that lack a reference genome.

CathaCyc, a metabolic pathway database built from Catharanthus roseus RNA-Seq data

The medicinal plant Madagascar periwinkle (Catharanthus roseus) synthesizes numerous terpenoid indole alkaloids (TIAs), such as the anticancer drugs vinblastine and vincristine. The TIA pathway operates in a complex metabolic network that steers plant growth and survival. Pathway databases and metabolic networks reconstructed from 'omics' sequence data can help to discover missing enzymes, study metabolic pathway evolution and, ultimately, engineer metabolic pathways. To date, such databases have mainly been built for model plant species with sequenced genomes. Although genome sequence data are not available for most medicinal plant species, next-generation sequencing is now extensively employed to create comprehensive medicinal plant transcriptome sequence resources. Here we report on the construction of CathaCyc, a detailed metabolic pathway database, from C. roseus RNA-Seq data sets. CathaCyc (version 1.0) contains 390 pathways with 1,347 assigned enzymes and spans primary and secondary metabolism. Curation of the pathways linked with the synthesis of TIAs and triterpenoids, their primary metabolic precursors, and their elicitors, the jasmonate hormones, demonstrated that RNA-Seq resources are suitable for the construction of pathway databases. CathaCyc is accessible online (http://www.cathacyc.org) and offers a range of tools for the visualization and analysis of metabolic networks and 'omics' data. Overlay with expression data from publicly available RNA-Seq resources demonstrated that two well-characterized C. roseus terpenoid pathways, those of TIAs and triterpenoids, are subject to distinct regulation by both developmental and environmental cues. We anticipate that databases such as CathaCyc will become key to the study and exploitation of the metabolism of medicinal plants.

De novo assembly, gene annotation and marker development using Illumina paired-end transcriptome sequences in celery (Apium graveolens L.)

Background

Celery is an increasing popular vegetable species, but limited transcriptome and genomic data hinder the research to it. In addition, a lack of celery molecular markers limits the process of molecular genetic breeding. High-throughput transcriptome sequencing is an efficient method to generate a large transcriptome sequence dataset for gene discovery, molecular marker development and marker-assisted selection breeding.

Principal Findings

Celery transcriptomes from four tissues were sequenced using Illumina paired-end sequencing technology. De novo assembling was performed to generate a collection of 42,280 unigenes (average length of 502.6 bp) that represent the first transcriptome of the species. 78.43% and 48.93% of the unigenes had significant similarity with proteins in the National Center for Biotechnology Information (NCBI) non-redundant protein database (Nr) and Swiss-Prot database respectively, and 10,473 (24.77%) unigenes were assigned to Clusters of Orthologous Groups (COG). 21,126 (49.97%) unigenes harboring Interpro domains were annotated, in which 15,409 (36.45%) were assigned to Gene Ontology(GO) categories. Additionally, 7,478 unigenes were mapped onto 228 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Large numbers of simple sequence repeats (SSRs) were indentified, and then the rate of successful amplication and polymorphism were investigated among 31 celery accessions.

Conclusions

This study demonstrates the feasibility of generating a large scale of sequence information by Illumina paired-end sequencing and efficient assembling. Our results provide a valuable resource for celery research. The developed molecular markers are the foundation of further genetic linkage analysis and gene localization, and they will be essential to accelerate the process of breeding.

Tapping the promise of genomics in species with complex, nonmodel genomes

Genomics is enabling a renaissance in all disciplines of plant biology. However, many plant genomes are complex and remain recalcitrant to current genomic technologies. The complexities of these nonmodel plant genomes are attributable to gene and genome duplication, heterozygosity, ploidy, and/or repetitive sequences. Methods are available to simplify the genome and reduce these barriers, including inbreeding and genome reduction, making these species amenable to current sequencing and assembly methods. Some, but not all, of the complexities in nonmodel genomes can be bypassed by sequencing the transcriptome rather than the genome. Additionally, comparative genomics approaches, which leverage phylogenetic relatedness, can aid in the interpretation of complex genomes. Although there are limitations in accessing complex nonmodel plant genomes using current sequencing technologies, genome manipulation and resourceful analyses can allow access to even the most recalcitrant plant genomes.

Engineering plant metabolism into microbes: from systems biology to synthetic biology

Plant metabolism represents an enormous repository of compounds that are of pharmaceutical and biotechnological importance. Engineering plant metabolism into microbes will provide sustainable solutions to produce pharmaceutical and fuel molecules that could one day replace substantial portions of the current fossil-fuel based economy. Metabolic engineering entails targeted manipulation of biosynthetic pathways to maximize yields of desired products. Recent advances in Systems Biology and the emergence of Synthetic Biology have accelerated our ability to design, construct and optimize cell factories for metabolic engineering applications. Progress in predicting and modeling genome-scale metabolic networks, versatile gene assembly platforms and delicate synthetic pathway optimization strategies has provided us exciting opportunities to exploit the full potential of cell metabolism. In this review, we will discuss how systems and synthetic biology tools can be integrated to create tailor-made cell factories for efficient production of natural products and fuel molecules in microorganisms.