De Novo transcriptome assembly (NGS) of Curcuma longa L. rhizome reveals novel transcripts related to anticancer and antimalarial terpenoids

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

De novo transcriptome assembly, gene annotation, and EST-SSR marker development of an important medicinal and edible crop, Amomum tsaoko (Zingiberaceae)

BACKGROUND

Amomum tsaoko is a medicinal and food dual-use crop that belongs to the Zingiberaceae family. However, the lack of transcriptomic and genomic information has limited the understanding of the genetic basis of this species. Here, we performed transcriptome sequencing of samples from different A. tsaoko tissues, and identified and characterized the expressed sequence tag-simple sequence repeat (EST-SSR) markers.

RESULTS

A total of 58,278,226 high-quality clean reads were obtained and de novo assembled to generate 146,911 unigenes with an N50 length of 2002 bp. A total of 128,174 unigenes were successfully annotated by searching seven protein databases, and 496 unigenes were identified as annotated as putative terpenoid biosynthesis-related genes. Furthermore, a total of 55,590 EST-SSR loci were detected, and 42,333 primer pairs were successfully designed. We randomly selected 80 primer pairs to validate their polymorphism in A. tsaoko; 18 of these primer pairs produced distinct, clear, and reproducible polymorphisms. A total of 98 bands and 96 polymorphic bands were amplified by 18 pairs of EST-SSR primers for the 72 A. tsaoko accessions. The Shannon's information index (I) ranged from 0.477 (AM208) to 1.701 (AM242) with an average of 1.183, and the polymorphism information content (PIC) ranged from 0.223 (AM208) to 0.779 (AM247) with an average of 0.580, indicating that these markers had a high level of polymorphism. Analysis of molecular variance (AMOVA) indicated relatively low genetic differentiation among the six A. tsaoko populations. Cross-species amplification showed that 14 of the 18 EST-SSR primer pairs have transferability between 11 Zingiberaceae species.

CONCLUSIONS

Our study is the first to provide transcriptome data of this important medicinal and edible crop, and these newly developed EST-SSR markers are a very efficient tool for germplasm evaluation, genetic diversity, and molecular marker-assisted selection in A. tsaoko.

Designing a Novel Functional Peptide With Dual Antimicrobial and Anti-inflammatory Activities via in Silico Methods

As spider venom is composed of various bioactive substances, it can be utilized as a platform for discovering future therapeutics. Host defense peptides are great candidates for developing novel antimicrobial agents due to their multifunctional properties. In this study, novel functional peptides were rationally designed to have dual antibacterial and anti-inflammatory activities with high cytocompatibility. Based on a template sequence from the transcriptome of spider Agelena koreana, a series of via in silico analysis were conducted, incorporating web-based machine learning tools along with the alteration of amino acid residues. Two peptides, Ak-N’ and Ak-N’m, were designed and were subjected to functional validation. The peptides inhibited gram-negative and gram-positive bacteria by disrupting the outer and bacterial cytoplasmic membrane. Moreover, the peptides down-regulated the expression of pro-inflammatory mediators, tumor necrosis factor-α, interleukin (IL)-1β, and IL6. Along with low cytotoxicity, Ak-N’m was shown to interact with macrophage surface receptors, inhibiting both Myeloid differentiation primary response 88-dependent and TIR-domain-containing adapter-inducing interferon-β-dependent pathways of Toll-like receptor 4 signaling on lipopolysaccharide-stimulated THP-1-derived macrophages. Here, we rationally designed functional peptides based on the suggested in silico strategy, demonstrating new insights for utilizing biological resources as well as developing therapeutic agents with enhanced properties.

Development of EST-SSR markers based on transcriptome and its validation in ginger (Zingiber officinale Rosc.)

Ginger (Zingiber officinale Rosc.) is an economically important and valuable spice crop around the world. It is used as food, spice, condiment, and medicine. A considerable extent of genetic diversity in ginger occurs in the Western Ghats and North-Eastern India. However, genetic diversity studies at the molecular level in ginger is limited due to limited availability of genetic and genomic information. In the present study, for the first time, we have identified and validated expressed sequence tag (EST)-simple sequence repeat (SSR) markers from ginger. We obtained 16,790 EST-SSR loci from 78987 unigenes, and 4597 SSR loci in the predicted 76929 coding sequences from RNA-Seq assembled contigs of ginger through Illumina paired-end sequencing. Gene ontology results indicate that the unigenes with SSR loci participate in various biological processes such as metabolism, growth, and development in ginger. One hundred and twenty-five primer pairs were designed from unigenes and coding sequences. These primers were tested for PCR optimization, characterization, and amplification and identified 12 novel EST-SSR markers. Twelve flanking polymorphic EST-SSR primers were validated using 48 ginger genotypes representing North-Eastern India and different eco-geographical adaptations by PCR amplification and allele sizing through capillary electrophoresis. Twelve EST-SSR primers generated a total of 111 alleles with an average of 9.25 alleles per locus and allele sizes ranging between 115-189bp. This study implies that the SSR markers designed from transcriptome sequences provides ample EST-SSR resources, which are helpful for genetic diversity analysis of Zingiberaceae species and molecular verification of ginger genotypes.

Genome sequencing of turmeric provides evolutionary insights into its medicinal properties

Curcuma longa, or turmeric, is traditionally known for its immense medicinal properties and has diverse therapeutic applications. However, the absence of a reference genome sequence is a limiting factor in understanding the genomic basis of the origin of its medicinal properties. In this study, we present the draft genome sequence of C. longa, belonging to Zingiberaceae plant family, constructed using 10x Genomics linked reads and Oxford Nanopore long reads. For comprehensive gene set prediction and for insights into its gene expression, transcriptome sequencing of leaf tissue was also performed. The draft genome assembly had a size of 1.02 Gbp with ~70% repetitive sequences, and contained 50,401 coding gene sequences. The phylogenetic position of C. longa was resolved through a comprehensive genome-wide analysis including 16 other plant species. Using 5,388 orthogroups, the comparative evolutionary analysis performed across 17 species including C. longa revealed evolution in genes associated with secondary metabolism, plant phytohormones signaling, and various biotic and abiotic stress tolerance responses. These mechanisms are crucial for perennial and rhizomatous plants such as C. longa for defense and environmental stress tolerance via production of secondary metabolites, which are associated with the wide range of medicinal properties in C. longa.

The Role of Chromatographic and Electromigration Techniques in Foodomics

Foodomics is the discipline aimed at studying the prevention of diseases by food, identifying chemical, biological and biochemical food contaminants, determining changes in genetically modified foods, identifying biomarkers able to confirm the authenticity and quality of foods or studying the safety, quality and traceability of foods, among other issues. It is mainly based on the use of genomic, transcriptomic, proteomic and metabolomic tools, among others, in order to understand the effect of food on animals and humans at the level of genes, messenger ribonucleic acid, proteins and metabolites. Since the first definition of Foodomics, a reasonable number of works have shown the extremely high possibilities of this discipline, which is highly based on the use of advanced analytical hyphenated techniques - especially for proteomics and metabolomics. This book chapter aims at providing a general description of the role of chromatographic and electromigration techniques that are currently being applied to achieve the main objectives of Foodomics, particularly in the proteomic and metabolomic fields, since most published works have been focused on these approaches, and to highlight relevant applications.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

From time immemorial, humans have exploited plants as a source of food and medicines. The World Health Organization (WHO) has recorded 21,000 plants with medicinal value out of 300,000 species available worldwide. The promising modern "multi-omics" platforms and tools have been proven as functional platforms able to endow us with comprehensive knowledge of the proteome, genome, transcriptome, and metabolome of medicinal plant systems so as to reveal the novel connected genetic (gene) pathways, proteins, regulator sequences and secondary metabolite (molecule) biosynthetic pathways of various drug and protein molecules from a variety of plants with therapeutic significance. This review paper endeavors to abridge the contemporary advancements in research areas of multi-omics and the information involved in decoding its prospective relevance to the utilization of plants with medicinal value in the present global scenario. The crosstalk of medicinal plants with genomics, transcriptomics, proteomics, and metabolomics approaches will be discussed.

The molecular mechanisms of Chlorella sp. responding to high CO2: A study based on comparative transcriptome analysis between strains with high- and low-CO2 tolerance

High CO₂ acclimation for microalgae has attracted large research attention owing to the usefulness of microalgae in bio-sequestration of CO₂ from the emission source. In this study, one high CO₂ tolerant (LAMB 31) and non-tolerant (LAMB 122) Chlorella sp. strains were transferred from air to 40% CO₂, during which four time points were chosen for comparative transcriptome analysis. Gene changes started in the lag phase (T1) of population growth with more genes (7889) upregulated in LAMB 31 than in LAMB 122 (1092). Further function enrichments indicated: In LAMB 31, up-regulation of genes in cyclic electron transportation, F-type ATPase and Calvin cycle were associated with the enhancement of carbon fixation abilities; upregulation of genes in phosphorylation together with V-ATPase, which contributed to cytoplasmatic pH stability; Lastly, enhancement of carbon metabolisms including TCA cycle and glycolysis accelerated the consumption of cellular organic carbon. Most of the genes in these pathways and processes showed downregulation in LAMB 122. This study disclosed the most complete transcriptional molecular mechanisms of Chlorella sp. responding to high CO₂ by combining CO₂ fixation, transportation, and metabolic processes. The results provided valuable genetic information for future screening and breeding of microalgae with high-CO₂ tolerance for more efficient CO₂ bio-sequestration.

Associating gene expressions with curcuminoid biosynthesis in turmeric

BACKGROUND

Biologically important curcuminoids viz curcumin, demethoxycurcumin, and bisdemethoxycurcumin in turmeric rhizome have several health benefits. Pharmaceutical industries synthesize curcuminoids manipulating gene expressions in vitro or in vivo because of their medicinal importance. In this experiment, we studied the gene expressions involved in the curcuminoid biosynthesis pathway in association with curcuminoid yield in turmeric rhizome to study the impact of individual gene expression on curcuminoid biosynthesis.

RESULTS

Gene expressions at the different growth stages of turmeric rhizome were association tested with respective curcuminoid contents. Gene expression patterns of diketide-CoA synthase (DCS) and multiple curcumin synthases (CURSs) viz curcumin synthase 1 (CURS1), curcumin synthase 2 (CURS2), and curcumin synthase 3 (CURS3) were differentially associated with different curcuminoid contents. Genotype and growth stages showed a significant effect on the gene expressions resulting in a significant impact on curcuminoid balance in turmeric rhizome. DCS and CURS3 expression patterns were similar but distinct from CURS1 and CURS2 expression patterns in turmeric rhizome. DCS expression had a significant positive and CURS3 expression had a significant negative association with curcumin, demethoxycurcumin, bisdemethoxycurcumin, and total curcuminoid in turmeric rhizome. CURS1 expression had a negative association with curcumin whereas CURS2 expression showed a positive association with demethoxycurcumin.

CONCLUSIONS

The effects of DCS and CURS expressions are not always positive with different curcuminoid contents in turmeric rhizome. DCS expression has a positive and CURS3 expression has a negative association with curcuminoids. CURS1 and CURS2 are also associated with curcumin and demethoxycurcumin synthesis. This mechanism of co-expression of diketide-CoA synthase and multiple curcumin synthases in turmeric rhizome has a significant effect on curcuminoid balance in different turmeric cultivars. Further experiment will explore more insights; however, association-based test results from this experiment can be useful in improving curcuminoid yield in turmeric rhizome or in vitro through the application of genetic engineering and biotechnology. Associating gene expressions with curcuminoid biosynthesis in turmeric.

Transcriptome analysis of Curcuma wenyujin from Haikou and Wenzhou, and a comparison of the main constituents and related genes of Rhizoma Curcumae

For more than a thousand years, Rhizoma Curcumae (known as E zhu), a Chinese herbal medicine, has been used to eradicate blood stasis and relieve aches. The plant Curcuma wenyujin, which is grown primarily in Wenzhou, China, is considered the best source of Rhizoma Curcumae. In this study, we sought to ascertain differences in transcript profiles of C. wenyujin grown in traditional (Wenzhou) and recently established (Haikou) production areas based on Illumina and RNA (RNA-seq) sequencing. We also examined differences in the main components of the volatile oil terpene; curcumin, polysaccharide, and starch constituents and related genes in the corresponding pathways, in C. wenyujin cultivated in the two production areas. We accordingly found that the essential oil (2.05%), curcumin (1.46%), and polysaccharide (8.90%) content in Wenzhou rhizomes was higher than that in the rhizomes of plants from Haikou (1.60%, 0.91%, and 6.15%, respectively). In contrast, the starch content of Wenzhou rhizomes (17.0%) was lower than that of Haikou rhizomes (23.8%). Furthermore, we detected significant differences in the oil components of Haikou and Wenzhou rhizomes, with curzerene (32.34%), curdione (21.35%), and germacrene B (9.39%) being the primary components of the essential oil derived from Wenzhou rhizomes, and curzerene (20.13%), curdione (14.73%), and cineole (9.76%) being the main constituents in Haikou rhizomes. Transcriptome and qPCR analyses revealed considerable differences in gene expression between Wenzhou and Haikou rhizomes. The expression of terpene, curcumin, and polysaccharide pathway-related genes in Wenzhou rhizomes was significantly up-regulated, whereas the expression of starch-associated genes was significantly down-regulated, compared with those in Haikou rhizomes. Difference in the content of terpene, curcumin, polysaccharides, and starch in rhizomes from the two production areas could be explained in terms of differences in expression of the related genes.

Transcriptome Sequencing and Development of Novel Genic SSR Markers From Pistacia vera L

In this study, we aimed to develop novel genic simple sequence repeat (eSSR) markers and to study phylogenetic relationship among Pistacia species. Transcriptome sequencing was performed in different tissues of Siirt and Atl cultivars of pistachio (Pistacia vera). A total of 37.5-Gb data were used in the assembly. The number of total contigs and unigenes was calculated as 98,831, and the length of N50 was 1,333 bp after assembly. A total of 14,308 dinucleotide, trinucleotide, tetranucleotide, pentanucleotide, and hexanucleotide SSR motifs (4–17) were detected, and the most abundant SSR repeat types were trinucleotide (29.54%), dinucleotide (24.06%), hexanucleotide (20.67%), pentanucleotide (18.88%), and tetranucleotide (6.85%), respectively. Overall 250 primer pairs were designed randomly and tested in eight Pistacia species for amplification. Of them, 233 were generated polymerase chain reaction products in at least one of the Pistacia species. A total of 55 primer pairs that had amplifications in all tested Pistacia species were used to characterize 11 P. vera cultivars and 78 wild Pistacia genotypes belonging to nine Pistacia species (P. khinjuk, P. eurycarpa, P. atlantica, P. mutica, P. integerrima, P. chinensis, P. terebinthus, P. palaestina, and P. lentiscus). A total of 434 alleles were generated from 55 polymorphic eSSR loci with an average of 7.89 alleles per locus. The mean number of effective allele was 3.40 per locus. Polymorphism information content was 0.61, whereas observed (Ho) and expected heterozygosity (He) values were 0.39 and 0.65, respectively. UPGMA (unweighted pair-group method with arithmetic averages) and STRUCTURE analysis divided 89 Pistacia genotypes into seven populations. The closest species to P. vera was P. khinjuk. P. eurycarpa was closer P. atlantica than P. khinjuk. P. atlantica–P. mutica and P. terebinthus–P. palaestina pairs of species were not clearly separated from each other, and they were suggested as the same species. The present study demonstrated that eSSR markers can be used in the characterization and phylogenetic analysis of Pistacia species and cultivars, as well as genetic linkage mapping and QTL (quantitative trait locus) analysis.

The GAMYB gene in rye: sequence, polymorphisms, map location, allele-specific markers, and relationship with α-amylase activity

BACKGROUND

Transcription factor (TF) GAMYB, belonging to MYB family (named after the gene of the avian myeloblastosis virus) is a master gibberellin (GA)-induced regulatory protein that is crucial for development and germination of cereal grain and involved in anther formation. It activates many genes including high-molecular-weight glutenin and α-amylase gene families. This study presents the first attempt to characterize the rye gene encoding GAMYB in relation to its sequence, polymorphisms, and phenotypic effects.

RESULTS

ScGAMYB was mapped on rye chromosome 3R using high-density Diversity Arrays Technology (DArT)/DArTseq-based maps developed in three mapping populations. The ScGAMYB sequences were identified in RNA-seq libraries of four rye inbred lines. The transcriptome used for the search contained almost 151,000 transcripts with a median contig length of 500 nt. The average amount of total base raw data was approximately 9 GB. Comparative analysis of the ScGAMYB sequence revealed its high level of homology to wheat and barley orthologues. Single nucleotide polymorphisms (SNPs) detected among rye inbred lines allowed the development of allele specific-PCR (AS-PCR) markers for ScGAMYB that might be used to detect this gene in wide genetic stocks of rye and triticale. Segregation of the ScGAMYB alleles showed significant relationship with α-amylase activity (AMY).

CONCLUSIONS

The research showed the strong similarity of rye GAMYB sequence to its orthologues in other Graminae and confirmed the position in the genome consistent with the collinearity rule of cereal genomes. Concurrently, the ScGAMYB coding sequence (cds) showed stronger variability (24 SNPs) compared to the analogous region of wheat (5 SNPs) and barley (7 SNPs). The moderate regulatory effect of ScGAMYB on AMY was confirmed, therefore, ScGAMYB was identified as a candidate gene for partial control of α-amylase production in rye grain. The predicted structural protein change in the aa region 362-372, caused by a single SNP (C/G) at the 1100 position in ScGAMYB cds and single aa sequence change (S/C) at the 367 position, is the likely cause of the differences in the effectiveness of ScGAMYB regulatory function associated with AMY. The development of sequence-based, allele-specific (AS) PCR markers could be useful in research and application.

De novo assembly and characterization of transcriptome in the medicinal plant Euphorbia jolkini

BACKGROUND

Euphorbia jolkini, a medicinal herb that grows on the warm beaches in Japan and South Korea, is known to be used for traditional medicines to treat a variety of ailments, including bruises, stiffness, indigestion, toothache, and diabetes.

OBJECTIVE

It is to analyze the whole transcriptome and identify the genes related to the phenylpropanoid biosynthesis in the medicinally important herb E jolkini.

METHODS

Paired-end Illumina HiSeq™ 2500 sequencing technology was employed for cDNA library construction and Illumina sequencing. Public databases like TAIR (The Arabidopsis Information Resource), Swissprot and KEGG (Kyoto Encyclopedia of Genes and Genomes) were used for annotations of unigenes obtained.

RESULTS

The transcriptome of E. jolkini generated 139,215 assembled transcripts with an average length of 868 bp and an N50 value of 1460 bp that were further clustered using CD-HIT into 93,801 unigenes with an average length of 847 bp (N50-1410 bp). Sixty-three percent of the coding sequences (CDS) were annotated from the longest open reading frame (ORF). A remarkable percentage of unigenes were annotated against various databases. The differentially expressed gene analysis revealed that the expression of genes related to the terpenoid backbone biosynthesis pathway was higher in the flowers, whereas that of genes related to the phenylpropanoid biosynthesis pathway was both up- and downregulated in flowers and leaves. A search of against the transcription factor domain found 1023 transcription factors (TFs) that were from 54 TF families.

CONCLUSION

Assembled sequences of the E. jolkini transcriptome are made available for the first time in this study E. jolkini and lay a foundation for the investigation of secondary metabolite biosynthesis.

Antibacterial and Anti-Inflammatory Effects of Novel Peptide Toxin from the Spider Pardosa astrigera

The prevalence of antibiotic-resistant bacteria has become an immediate threat to public health. Antimicrobial peptides are attracting attention as a new source of antibiotics due to their ability to prevent drug-resistances with fewer side effects. Spider venom is composed of various bioactive substances with multiple functionalities such as antimicrobial and anti-inflammatory effects. Here, RNA sequencing was conducted on the venom gland of the spider Pardosa astrigera, and a potential toxin peptide with antibacterial properties was selected via homology and in silico analysis. A novel toxin, Lycotoxin-Pa4a, inhibited both gram-negative and gram-positive bacteria by disrupting the outer and bacterial cytoplasmic membrane. Moreover, the peptide downregulated the expression of pro-inflammatory mediators while upregulating the level of anti-inflammatory cytokine by inactivating mitogen-activated protein kinase signaling in a lipopolysaccharide-stimulated murine macrophage cell line. In this research, we identified a novel peptide toxin, Lycotoxin-pa4a, with antibacterial and anti-inflammatory properties, suggesting its potential for the development of a new antibiotics, as well as offering insights into the utilization of biological resources.

De novo transcriptome analysis unravels tissue-specific expression of candidate genes involved in major secondary metabolite biosynthetic pathways of Plumbago zeylanica: implication for pharmacological potential

KEY MESSAGE

The present study provides comparative transcriptome analysis, besides identifying functional secondary metabolite genes of Plumbago zeylanica with pharmacological potential for future functional genomics, and metabolomic engineering of secondary metabolites from this plant towards diversified biomedical applications.

ABSTRACT

Plumbago zeylanica is a widely used medicinal plant of the traditional Indian system of medicine with wide pharmacological potential to treat several disorders. The present study aimed to carry out comparative transcriptome analysis in leaf and root tissue of P. zeylanica using Illumina paired end sequencing to identify tissue-specific functional genes involved in the biosynthesis of secondary metabolites, contributing to its therapeutic efficacy. De novo sequencing assembly resulted in the identification of 62,321 "Unigenes" transcripts with an average size of 1325 bp. Functional annotation using BLAST2GO resulted in the identification of 50,301 annotated transcripts (80.71%) and GO assigned to 18,814 transcripts. KEGG pathway annotation of the "Unigenes" revealed that 2465 transcripts could be assigned to 242 KEGG pathway maps wherein the number of transcripts involved in secondary metabolism was distinct in root and leaf transcriptome. Among the secondary metabolite biosynthesis pathways, the cluster of "Unigenes" encoding enzymes of 'Phenylpropanoid biosynthesis pathway' represents the largest group (84 transcripts) followed by 'Terpenoid Backbone biosynthesis' (48 transcripts). The transcript levels of the candidate unigenes encoding key enzymes of phenylpropanoid (PAL, TAL) and flavanoid biosynthesis (CHS, ANS, FLS) pathways were up-regulated in root, while the expression levels of candidate "Unigenes" transcript for monoterpenoid (DXS, ISPF), diterpenoid biosynthesis (SPS, SDS) and indole alkaloid pathways (STR) were significantly higher in leaf of P. zeylanica. Interestingly, validation of differential gene expression profile by qRT-PCR also confirmed that candidate "Unigenes" enzymes of phenylpropanoid and flavonoid biosynthesis were highly expressed in the root, while the key regulatory enzymes of terpenoid and indole alkaloid compounds were up-regulated in the leaf, suggesting that (differences in) the levels of these functional genes could be attributed to the (differential) pharmacological activity (between root and leaf) in tissues of P. zeylanica.

De novo Transcriptome Assembly, Gene Annotation and SSR Marker Development in the Moon Seed Genus Menispermum (Menispermaceae)

The moonseed genus Menispermum L. (Menispermaceae) is disjunctly distributed in East Asia and eastern North America. Although Menispermum has important medicinal value, genetic and genomic information is scarce, with very few available molecular markers. In the current study, we used Illumina transcriptome sequencing and de novo assembly of the two Menispermum species to obtain in-depth genetic knowledge. From de novo assembly, 53,712 and 78,921 unigenes were generated for M. canadense and M. dauricum, with 37,527 (69.87%) and 55,211 (69.96%) showing significant similarities against the six functional databases, respectively. Moreover, 521 polymorphic EST-SSRs were identified. Of them, 23 polymorphic EST-SSR markers were selected to investigate the population genetic diversity within the genus. The newly developed EST-SSR markers also revealed high transferability among the three examined Menispermaceae species. Overall, we provide the very first transcriptomic analyses of this important medicinal genus. In addition, the novel microsatellite markers developed here will aid future studies on the population genetics and phylogeographic patterns of Menispermum at the intercontinental geographical scale.

Genetic Divergence and Relationship Among Opisthopappus Species Identified by Development of EST-SSR Markers

Opisthopappus Shih is an endemic and endangered genus restricted to the Taihang Mountains that has important ornamental and economic value. According to the Flora Reipublicae Popularis Sinicae (FRPS, Chinese version), this genus contains two species (Opisthopappus longilobus and Opisthopappus taihangensis), whereas in the Flora of China (English version) only one species O. taihangensis is present. The interspecific phylogenetic relationship remains unclear and undefined, which might primarily be due to the lack of specific molecular markers for phylogenetic analysis. For this study, 2644 expressed sequence tag-simple sequence repeats (EST-SSRs) from 33,974 unigenes using a de novo transcript assembly of Opisthopappus were identified with a distribution frequency of 7.78% total unigenes. Thereinto, mononucleotides (1200, 45.39%) were the dominant repeat motif, followed by trinucleotides (992, 37.52%), and dinucleotides (410, 15.51%). The most dominant trinucleotide repeat motif was ACC/GGT (207, 20.87%). Based on the identified EST-SSRs, 245 among 1444 designed EST-SSR primers were selected for the development of potential molecular markers. Among these markers, 63 pairs of primers (25.71%) generated clear and reproducible bands with expected sizes. Eventually, 11 primer pairs successfully amplified all individuals from the studied populations. Through the EST-SSR markers, a high level of genetic diversity was detected between Opisthopappus populations. A significant genetic differentiation between the O. longilobus and O. taihangensis populations was found. All studied populations were divided into two clusters by UPGMA, NJ, STRUCTURE, and PCoA. These results fully supported the view of the FRPS, namely, that O. longilobus and O. taihangensis should be regarded as two distinct species. Our study demonstrated that transcriptome sequences, as a valuable tool for the quick and cost-effective development of molecular markers, was helpful toward obtaining comprehensive EST-SSR markers that could contribute to an in-depth assessment of the genetic and phylogenetic relationships between Opisthopappus species.

De novo transcriptome of Gymnema sylvestre identified putative lncRNA and genes regulating terpenoid biosynthesis pathway

Gymnema sylvestre is a highly valuable medicinal plant in traditional Indian system of medicine and used in many polyherbal formulations especially in treating diabetes. However, the lack of genomic resources has impeded its research at molecular level. The present study investigated functional gene profile of G. sylvestre via RNA sequencing technology. The de novo assembly of 88.9 million high quality reads yielded 23,126 unigenes, of which 18116 were annotated against databases such as NCBI nr database, gene ontology (GO), KEGG, Pfam, CDD, PlantTFcat, UniProt & GreeNC. Total 808 unigenes mapped to 78 different Transcription Factor families, whereas 39 unigenes assigned to CYP450 and 111 unigenes coding for enzymes involved in the biosynthesis of terpenoids including transcripts for synthesis of important compounds like Vitamin E, beta-amyrin and squalene. Among them, presence of six important enzyme coding transcripts were validated using qRT-PCR, which showed high expression of enzymes involved in methyl-erythritol phosphate (MEP) pathway. This study also revealed 1428 simple sequence repeats (SSRs), which may aid in molecular breeding studies. Besides this, 8 putative long non-coding RNAs (lncRNAs) were predicted from un-annotated sequences, which may hold key role in regulation of essential biological processes in G. sylvestre. The study provides an opportunity for future functional genomic studies and to uncover functions of the lncRNAs in G. sylvestre.

The transcriptome enables the identification of candidate genes behind medicinal value of Drumstick tree (Moringa oleifera)

Moringa oleifera is a plant well-known for its nutrition value, drought resistance and medicinal properties. cDNA libraries from five different tissues (leaf, root, stem, seed and flower) of M. oleifera cultivar Bhagya were generated and sequenced. We developed a bioinformatics pipeline to assemble transcriptome, along with the previously published M. oleifera genome, to predict 17,148 gene models. Few candidate genes related to biosynthesis of secondary metabolites, vitamins and ion transporters were identified. Expressions were further confirmed by real-time quantitative PCR experiments for few promising leads. Quantitative estimation of metabolites, as well as elemental analysis, was also carried out to support our observations. Enzymes in the biosynthesis of vitamins and metabolites like quercetin and kaempferol are highly expressed in leaves, flowers and seeds. The expression of iron transporters and calcium storage proteins were observed in root and leaves. In general, leaves retain the highest amount of small molecules of interest.

De novo assembly of transcriptomes, mining, and development of novel EST-SSR markers in Curcuma alismatifolia (Zingiberaceae family) through Illumina sequencing

Curcuma alismatifolia widely used as an ornamental plant in Thailand and Cambodia. This species of herbaceous perennial from the Zingiberaceae family, includes cultivars with a wide range of colours and long postharvest life, and is used as an ornamental cut flower, as a potted plant, and in exterior landscapes. For further genetic improvement, however, little genomic information and no specific molecular markers are available. The present study used Illumina sequencing and de novo transcriptome assembly of two C. alismatifolia cvs, 'Chiang Mai Pink' and 'UB Snow 701', to develop simple sequence repeat markers for genetic diversity studies. After de novo assembly, 62,105 unigenes were generated and 48,813 (78.60%) showed significant similarities versus six functional protein databases. In addition, 9,351 expressed sequence tag-simple sequence repeats (EST-SSRs) were identified with a distribution frequency of 12.5% total unigenes. Out of 8,955 designed EST-SSR primers, 150 primers were selected for the development of potential molecular markers. Among these markers, 17 EST-SSR markers presented a moderate level of genetic diversity among three C. alismatifolia cultivars, one hybrid, three Curcuma, and two Zingiber species. Three different genetic groups within these species were revealed using EST-SSR markers, indicating that the markers developed in this study can be effectively applied to the population genetic analysis of Curcuma and Zingiber species. This report describes the first analysis of transcriptome data of an important ornamental ginger cultivars, also provides a valuable resource for gene discovery and marker development in the genus Curcuma.

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.

Molecular authentication of Anthemis deserti Boiss. (Asteraceae) based on ITS2 region of nrDNA gene sequence

The dried plant material of medicinally important Anthemis deserti Boiss. (family: Asteraceae) especially when it remains in the powdered form often look similar to Anthemis melampodina Del.; and therefore, difficult to distinguish, finally lead to chances of adulteration. The adulteration in medicinal plants effects on the efficacy of the drugs. The molecular authentication of herbal plant materials such as based on the internal transcribed spacer 2 (ITS2) sequences of nuclear ribosomal DNA (nrDNA) is considered as more reliable method compared to other the biochemical or histological methods. The present study aims to molecular authentication ofA. deserti based on molecular phylogenetic analyses of ITS2 gene sequence of nrDNA region. The ITS2 region of nrDNA of A. deserti were sequenced, and the molecular phylogenetic analyses were performed together with the GenBank sequences. The Maximum Parsimony tree revealed the close relationships of A. deserti with A. melampodina; however, the Neighbor-Joining and Maximum Likelihood tree clearly revealed that A. deserti is distinct from A. melampodina, which is also supported by the differences in nucleotides at five diffident positions (i.e. 22, 28, 87, 175 and 198) in the DNA sequence alignment.

Transcriptome analysis, using RNA-Seq of Lomandra longifolia roots infected with Phytophthora cinnamomi reveals the complexity of the resistance response

The plant pathogen Phytophthora cinnamon the causal agent of disease in numerous species, is a major threat to natural vegetation and has economic impacts in agriculture. The pathogen principally invades the root system, which, in susceptible species, is rapidly colonised and functionally destroyed. Few species are resistant, however, where resistance is expressed the pathogen is restricted to small, localised lesions. The molecular mechanisms that underpin this response in resistant species are not well understood. Lomandra longifolia, an Australian native species, is highly resistant to P. cinnamomi. In an earlier study, we showed induction of resistance-related components such as callose, lignin and hydrogen peroxide (H₂ O₂ ) in L. longifolia roots that had been inoculated with P. cinnamomi. Here, in order to further identify, during the very early stages of infection, the molecular components and regulatory networks that may trigger resistance, a comprehensive root transcriptome analysis was performed using next generation sequencing. Overall, 18 cDNA libraries were produced generating 52.8 GB 126 base pair reads, which were de novo assembled into contigs. Differentially expressed genes (DEGs) were identified allowing the identification of infection-responsive candidate genes that were putatively related to resistance, and from this set ten were selected for qRT-PCR to validate the RNA-Seq expression value. Further analysis of individual candidates revealed that many were involved in PAMP-triggered immunity (PTI; pattern recognition receptors, glutathione S-transferase, callose synthases, pathogenesis-related protein-1, mitogen activated protein kinases) and effector-triggered immunity (ETI) (NBS-LRR, signalling genes, transcription factors and anti-pathogenic compound synthase genes). As these candidate genes or mediated components activate different defence signalling systems, they may have potential for investigation of novel approaches to disease control and in transgenic approaches for improvement, in susceptible species, of resistance to P. cinnamomi.

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.

Molecular insights into photosynthesis and carbohydrate metabolism in Jatropha curcas grown under elevated CO2 using transcriptome sequencing and assembly

Jatropha curcas L. (Family - Euphorbiaceae) is a perennial tree of special interest due to its potential as a biofuel plant with high carbon sequestration. In this study, physiological investigations coupled with transcriptomics in relation to photosynthesis were evaluated in Jatropha grown under ambient (395 ppm) and elevated (550 ppm) CO2 atmosphere. Morphophysiological analysis revealed that Jatropha sustained enhanced photosynthesis during its growth under elevated CO2 for one year which might be linked to improved CO2 assimilation physiology and enhanced sink activity. We sequenced and analyzed the leaf transcriptome of Jatropha after one year of growth in both conditions using Illumina HiSeq platform. After optimized assembly, a total of 69,581 unigenes were generated. The differential gene expression (DGE) analysis revealed 3013 transcripts differentially regulated in elevated CO2 conditions. The photosynthesis regulatory genes were analysed for temporal expression patterns at four different growth phases which highlighted probable events contributing to enhanced growth and photosynthetic capacity including increased reducing power, starch synthesis and sucrose mobilization under elevated CO2. Overall, our data on physiological and transcriptomic analyses suggest an optimal resource allocation to the available and developing sink organs thereby sustaining improved photosynthetic rates during long-term growth of Jatropha under CO2 enriched environment.

High-quality draft genome sequence of a biofilm forming lignocellulolytic Aspergillus niger strain ATCC 10864

Filamentous fungus Aspergillus niger has high industrial value due to their lignocellulolytic enzyme activities and ATCC 10864 is one of the few type strains of A. niger which has a unique biofilm forming capability. Here we report the first draft genome sequence of A. niger ATCC 10864 strain. The genome of A. niger ATCC 10864 is 36,172,237 bp long and comprise of 310 scaffolds with 49.5% average GC content. A total of 10,804 protein-coding genes were predicted among which 10,761 genes were with putative functions. A. niger ATCC 10864 genome coded for 709 putative carbohydrate active enzyme families distributed in six functional categories and among them glycoside hydrolases (GHs) represent the most number of families (279). Genes that include pepA, brlA, exgA, LaeA, rodA, GCN have also been identified in this study, which may play a role in biofilm formation. This high-quality draft genome sequence will facilitate our understanding of the mechanisms behind fungal biofilm formation and higher lignocellulolytic enzyme production.

Genome and Transcriptome sequence of Finger millet (Eleusine coracana (L.) Gaertn.) provides insights into drought tolerance and nutraceutical properties

Background

Finger millet (Eleusine coracana (L.) Gaertn.) is an important staple food crop widely grown in Africa and South Asia. Among the millets, finger millet has high amount of calcium, methionine, tryptophan, fiber, and sulphur containing amino acids. In addition, it has C4 photosynthetic carbon assimilation mechanism, which helps to utilize water and nitrogen efficiently under hot and arid conditions without severely affecting yield. Therefore, development and utilization of genomic resources for genetic improvement of this crop is immensely useful.

Results

Experimental results from whole genome sequencing and assembling process of ML-365 finger millet cultivar yielded 1196 Mb covering approximately 82% of total estimated genome size. Genome analysis showed the presence of 85,243 genes and one half of the genome is repetitive in nature. The finger millet genome was found to have higher colinearity with foxtail millet and rice as compared to other Poaceae species. Mining of simple sequence repeats (SSRs) yielded abundance of SSRs within the finger millet genome. Functional annotation and mining of transcription factors revealed finger millet genome harbors large number of drought tolerance related genes. Transcriptome analysis of low moisture stress and non-stress samples revealed the identification of several drought-induced candidate genes, which could be used in drought tolerance breeding.

Conclusions

This genome sequencing effort will strengthen plant breeders for allele discovery, genetic mapping, and identification of candidate genes for agronomically important traits. Availability of genomic resources of finger millet will enhance the novel breeding possibilities to address potential challenges of finger millet improvement.

Electronic supplementary material

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

Spatial transcriptome analysis provides insights of key gene(s) involved in steroidal saponin biosynthesis in medicinally important herb Trillium govanianum

Trillium govanianum, an endangered medicinal herb native to the Himalaya, is less studied at the molecular level due to the non-availability of genomic resources. To facilitate the basic understanding of the key genes and regulatory mechanism of pharmaceutically important biosynthesis pathways, first spatial transcriptome sequencing of T. govanianum was performed. 151,622,376 (~11.5 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 69,174 transcripts. Functional annotation with multiple public databases identified array of genes involved in steroidal saponin biosynthesis and other secondary metabolite pathways including brassinosteroid, carotenoid, diterpenoid, flavonoid, phenylpropanoid, steroid and terpenoid backbone biosynthesis, and important TF families (bHLH, MYB related, NAC, FAR1, bZIP, B3 and WRKY). Differentially expressed large number of transcripts, together with CYPs and UGTs suggests involvement of these candidates in tissue specific expression. Combined transcriptome and expression analysis revealed that leaf and fruit tissues are the main site of steroidal saponin biosynthesis. In conclusion, comprehensive genomic dataset created in the current study will serve as a resource for identification of potential candidates for genetic manipulation of targeted bioactive metabolites and also contribute for development of functionally relevant molecular marker resource to expedite molecular breeding and conservation efforts in T. govanianum.

De novo sequencing and comparative analysis of leaf transcriptomes of diverse condensed tannin-containing lines of underutilized Psophocarpus tetragonolobus (L.) DC

Condensed tannin (CT) or proanthocyanidin (PA) is a unique group of phenolic metabolite with high molecular weight with specific structure. It is reported that, the presence of high-CT in the legumes adversely affect the nutrients in the plant and impairs the digestibility upon consumption by animals. Winged bean (Psophocarpus tetragonolobus (L.) DC.) is one of the promising underutilized legume with high protein and oil-content. One of the reasons for its underutilization is due to the presence of CT. Transcriptome sequencing of leaves of two diverse CT-containing lines of P. tetragonolobus was carried out on Illumina Nextseq 500 sequencer to identify the underlying genes and contigs responsible for CT-biosynthesis. RNA-Seq data generated 102586 and 88433 contigs for high (HCTW) and low CT (LCTW) lines of P. tetragonolobus, respectively. Based on the similarity searches against gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) database revealed 5210 contigs involved in 229 different pathways. A total of 1235 contigs were detected to differentially express between HCTW and LCTW lines. This study along with its findings will be helpful in providing information for functional and comparative genomic analysis of condensed tannin biosynthesis in this plant in specific and legumes in general.

Natural products against Alzheimer's disease: Pharmaco-therapeutics and biotechnological interventions

Alzheimer's disease (AD) is a severe, chronic and progressive neurodegenerative disease associated with memory and cognition impairment ultimately leading to death. It is the commonest reason of dementia in elderly populations mostly affecting beyond the age of 65. The pathogenesis is indicated by accumulation of the amyloid-beta (Aβ) plaques and neurofibrillary tangles (NFT) in brain tissues and hyperphosphorylation of tau protein in neurons. The main cause is considered to be the formation of reactive oxygen species (ROS) due to oxidative stress. The current treatment provides only symptomatic relief by offering temporary palliative therapy which declines the rate of cognitive impairment associated with AD. Inhibition of the enzyme acetylcholinesterase (AChE) is considered as one of the major therapeutic strategies offering only symptomatic relief and moderate disease-modifying effect. Other non-cholinergic therapeutic approaches include antioxidant and vitamin therapy, stem cell therapy, hormonal therapy, use of antihypertensive or lipid-lowering medications and selective phosphodiesterase (PDE) inhibitors, inhibition of β-secretase and γ-secretase and Aβ aggregation, inhibition of tau hyperphosphorylation and intracellular NFT, use of nonsteroidal anti-inflammatory drugs (NSAIDs), transition metal chelators, insulin resistance drugs, etanercept, brain-derived neurotrophic factor (BDNF) etc. Medicinal plants have been reported for possible anti-AD activity in a number of preclinical and clinical trials. Ethnobotany, being popular in China and in the Far East and possibly less emphasized in Europe, plays a substantial role in the discovery of anti-AD agents from botanicals. Chinese Material Medica (CMM) involving Chinese medicinal plants has been used traditionally in China in the treatment of AD. Ayurveda has already provided numerous lead compounds in drug discovery and many of these are also undergoing clinical investigations. A number of medicinal plants either in their crude forms or as isolated compounds have exhibited to reduce the pathological features associated with AD. In this present review, an attempt has been made to elucidate the molecular mode of action of various plant extracts, phytochemicals and traditional herbal formulations investigated against AD as reported in various preclinical and clinical tests. Herbal synergism often found in polyherbal formulations were found effective to combat disease heterogeneity as found in complex pathogenesis of AD. Finally a note has been added to describe biotechnological improvement, genetic and genomic resources and mathematical and statistical techniques for empirical model building associated with anti-AD plant secondary metabolites and their source botanicals.

Transcriptome Analysis of the Tadpole Shrimp (Triops longicaudatus) by Illumina Paired-End Sequencing: Assembly, Annotation, and Marker Discovery

The tadpole shrimp (Triops longicaudatus) is an aquatic crustacean that helps control pest populations. It inhabits freshwater ponds and pools and has been described as a living fossil. T. longicaudatus was officially declared an endangered species South Korea in 2005; however, through subsequent protection and conservation management, it was removed from the endangered species list in 2012. The limited number of available genetic resources on T. longicaudatus makes it difficult to obtain valuable genetic information for marker-aided selection programs. In this study, whole-transcriptome sequencing of T. longicaudatus generated 39.74 GB of clean data and a total of 269,822 contigs using the Illumina HiSeq 2500 platform. After clustering, a total of 208,813 unigenes with an N₅₀ length of 1089 bp were generated. A total of 95,105 unigenes were successfully annotated against Protostome (PANM), Unigene, Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases using BLASTX with a cut-off of 1E−5. A total of 57,731 unigenes were assigned to GO terms, and 7247 unigenes were mapped to 129 KEGG pathways. Furthermore, 1595 simple sequence repeats (SSRs) were detected from the unigenes with 1387 potential SSR markers. This is the first report of high-throughput transcriptome analysis of T. longicaudatus, and it provides valuable insights for genetic research and molecular-assisted breeding of this important species.

Transcriptomic Analysis of the Endangered Neritid Species Clithon retropictus: De Novo Assembly, Functional Annotation, and Marker Discovery

An aquatic gastropod belonging to the family Neritidae, Clithon retropictus is listed as an endangered class II species in South Korea. The lack of information on its genomic background limits the ability to obtain functional data resources and inhibits informed conservation planning for this species. In the present study, the transcriptomic sequencing and de novo assembly of C. retropictus generated a total of 241,696,750 high-quality reads. These assembled to 282,838 unigenes with mean and N50 lengths of 736.9 and 1201 base pairs, respectively. Of these, 125,616 unigenes were subjected to annotation analysis with known proteins in Protostome DB, COG, GO, and KEGG protein databases (BLASTX; E ≤ 0.00001) and with known nucleotides in the Unigene database (BLASTN; E ≤ 0.00001). The GO analysis indicated that cellular process, cell, and catalytic activity are the predominant GO terms in the biological process, cellular component, and molecular function categories, respectively. In addition, 2093 unigenes were distributed in 107 different KEGG pathways. Furthermore, 49,280 simple sequence repeats were identified in the unigenes (>1 kilobase sequences). This is the first report on the identification of transcriptomic and microsatellite resources for C. retropictus, which opens up the possibility of exploring traits related to the adaptation and acclimatization of this species.

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.

High-throughput sequencing reveals differential regulation of miRNAs in fenoxaprop-P-ethyl-resistant Beckmannia syzigachne

Non-target site resistance (NTSR) to herbicides is an increasing concern for weed control. The majority of previous studies have focused on metabolic resistance mechanisms of NTSR, but no research exists on gene regulation mechanisms behind herbicide resistance, such as microRNA (miRNA). Here, we identified 3 American sloughgrass (Beckmannia syzigachne Steud.) populations containing fenoxaprop-P-ethyl-resistant plants. We then constructed small RNA libraries and subjected them to deep sequencing and bioinformatics analyses. Forty known and 36 potentially novel, predicted miRNAs were successfully identified. Of these, we identified 3 conserved, predicted candidate NTSR-determinant miRNAs and their potential corresponding target genes, as well as 4 novel potential miRNAs with high count. Target gene prediction and annotation indicated that these 7 differentially expressed miRNAs potentially play a role in regulating specific stress-responsive genes, very likely related to herbicide resistance. Expression profiles were determined with quantitative real-time PCR. The present study is a novel, large-scale characterization of weed miRNAs. The results should further our understanding of miRNA expression profiles associated with herbicide resistance, allowing for the development of more effective weed management strategies.

Accessing Nature's diversity through metabolic engineering and synthetic biology

In this perspective, we highlight recent examples and trends in metabolic engineering and synthetic biology that demonstrate the synthetic potential of enzyme and pathway engineering for natural product discovery. In doing so, we introduce natural paradigms of secondary metabolism whereby simple carbon substrates are combined into complex molecules through “scaffold diversification”, and subsequent “derivatization” of these scaffolds is used to synthesize distinct complex natural products. We provide examples in which modern pathway engineering efforts including combinatorial biosynthesis and biological retrosynthesis can be coupled to directed enzyme evolution and rational enzyme engineering to allow access to the “privileged” chemical space of natural products in industry-proven microbes. Finally, we forecast the potential to produce natural product-like discovery platforms in biological systems that are amenable to single-step discovery, validation, and synthesis for streamlined discovery and production of biologically active agents.

Comparative transcriptome analysis of Gastrodia elata (Orchidaceae) in response to fungus symbiosis to identify gastrodin biosynthesis-related genes

BACKGROUND

Gastrodia elata Blume (Orchidaceae) is an important Chinese medicine with several functional components. In the life cycle of G. elata, the orchid develops a symbiotic relationship with two compatible mycorrhizal fungi Mycena spp. and Armillaria mellea during seed germination to form vegetative propagation corm and vegetative growth to develop tubers, respectively. Gastrodin (p-hydroxymethylphenol-beta-D-glucoside) is the most important functional component in G. elata, and gastrodin significantly increases from vegetative propagation corms to tubers. To address the gene regulation mechanism in gastrodin biosynthesis in G. elata, a comparative analysis of de novo transcriptome sequencing among the vegetative propagation corms and tubers of G. elata and A. mellea was conducted using deep sequencing.

RESULTS

Transcriptome comparison between the vegetative propagation corms and juvenile tubers of G. elata revealed 703 differentially expressed unigenes, of which 298 and 405 unigenes were, respectively up-regulated (fold-change ≥ 2, q-value < 0.05, the trimmed mean of M-values (TMM)-normalized fragments per kilobase of transcript per Million mapped reads (FPKM) > 10) and down-regulated (fold-change ≤ 0.5, q-value <0.05, TMM-normalized FPKM > 10) in juvenile tubers. After Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 112 up-regulated unigenes with KEGG Ortholog identifiers (KOids) or enzyme commission (EC) numbers were assigned to 159 isogroups involved in seventy-eight different pathways, and 132 down-regulated unigenes with KOids or EC numbers were assigned to 168 isogroups, involved in eighty different pathways. The analysis of the isogroup genes from all pathways revealed that the two unigenes TRINITY_DN54282_c0_g1 (putative monooxygenases) and TRINITY_DN50323_c0_g1 (putative glycosyltransferases) might participate in hydroxylation and glucosylation in the gastrodin biosynthetic pathway.

CONCLUSIONS

The gene expression of the two unique unigenes encoding monooxygenase and glycosyltransferase significantly increases from vegetative propagation corms to tubers, and the molecular basis of gastrodin biosynthesis in the tubers of G. elata is proposed.

Establishing a herbicide-metabolizing enzyme library in Beckmannia syzigachne to identify genes associated with metabolic resistance

Non-target site resistance (NTSR) to herbicides is an increasing concern for weed control. Metabolic herbicide resistance is an important mechanism for NTSR. However, little is known about metabolic resistance at the genetic level. In this study, we have identified three fenoxaprop-P-ethyl-resistant American sloughgrass (Beckmannia syzigachne Steud.) populations, in which the molecular basis for NTSR remains unclear. To reveal the mechanisms of metabolic resistance, the genes likely to be involved in herbicide metabolism (e.g. for cytochrome P450s, esterases, hydrolases, oxidases, peroxidases, glutathione S-transferases, glycosyltransferases, and transporter proteins) were isolated using transcriptome sequencing, in combination with RT-PCR (reverse transcription-PCR) and RACE (rapid amplification of cDNA ends). Consequently, we established a herbicide-metabolizing enzyme library containing at least 332 genes, and each of these genes was cloned and the sequence and the expression level compared between the fenoxaprop-P-ethyl-resistant and susceptible populations. Fifteen metabolic enzyme genes were found to be possibly involved in fenoxaprop-P-ethyl resistance. In addition, we found five metabolizing enzyme genes that have a different gene sequence in plants of susceptible versus resistant B. syzigachne populations. These genes may be major candidates for herbicide metabolic resistance. This established metabolic enzyme library represents an important step forward towards a better understanding of herbicide metabolism and metabolic resistance in this and possibly other closely related weed species. This new information may help to understand weed metabolic resistance and to develop novel strategies of weed management.

Sequencing, De novo Assembly, Functional Annotation and Analysis of Phyllanthus amarus Leaf Transcriptome Using the Illumina Platform

Phyllanthus amarus Schum. and Thonn., a widely distributed annual medicinal herb has a long history of use in the traditional system of medicine for over 2000 years. However, the lack of genomic data for P. amarus, a non-model organism hinders research at the molecular level. In the present study, high-throughput sequencing technology has been employed to enhance better understanding of this herb and provide comprehensive genomic information for future work. Here P. amarus leaf transcriptome was sequenced using the Illumina Miseq platform. We assembled 85,927 non-redundant (nr) "unitranscript" sequences with an average length of 1548 bp, from 18,060,997 raw reads. Sequence similarity analyses and annotation of these unitranscripts were performed against databases like green plants nr protein database, Gene Ontology (GO), Clusters of Orthologous Groups (COG), PlnTFDB, KEGG databases. As a result, 69,394 GO terms, 583 enzyme codes (EC), 134 KEGG maps, and 59 Transcription Factor (TF) families were generated. Functional and comparative analyses of assembled unitranscripts were also performed with the most closely related species like Populus trichocarpa and Ricinus communis using TRAPID. KEGG analysis showed that a number of assembled unitranscripts were involved in secondary metabolites, mainly phenylpropanoid, flavonoid, terpenoids, alkaloids, and lignan biosynthetic pathways that have significant medicinal attributes. Further, Fragments Per Kilobase of transcript per Million mapped reads (FPKM) values of the identified secondary metabolite pathway genes were determined and Reverse Transcription PCR (RT-PCR) of a few of these genes were performed to validate the de novo assembled leaf transcriptome dataset. In addition 65,273 simple sequence repeats (SSRs) were also identified. To the best of our knowledge, this is the first transcriptomic dataset of P. amarus till date. Our study provides the largest genetic resource that will lead to drug development and pave the way in deciphering various secondary metabolite biosynthetic pathways in P. amarus, especially those conferring the medicinal attributes of this potent herb.

De novo Assembly of Leaf Transcriptome in the Medicinal Plant Andrographis paniculata

Andrographis paniculata is an important medicinal plant containing various bioactive terpenoids and flavonoids. Despite its importance in herbal medicine, no ready-to-use transcript sequence information of this plant is made available in the public data base, this study mainly deals with the sequencing of RNA from A. paniculata leaf using Illumina HiSeq™ 2000 platform followed by the de novo transcriptome assembly. A total of 189.22 million high quality paired reads were generated and 1,70,724 transcripts were predicted in the primary assembly. Secondary assembly generated a transcriptome size of ~88 Mb with 83,800 clustered transcripts. Based on the similarity searches against plant non-redundant protein database, gene ontology, and eukaryotic orthologous groups, 49,363 transcripts were annotated constituting upto 58.91% of the identified unigenes. Annotation of transcripts-using kyoto encyclopedia of genes and genomes database-revealed 5606 transcripts plausibly involved in 140 pathways including biosynthesis of terpenoids and other secondary metabolites. Transcription factor analysis showed 6767 unique transcripts belonging to 97 different transcription factor families. A total number of 124 CYP450 transcripts belonging to seven divergent clans have been identified. Transcriptome revealed 146 different transcripts coding for enzymes involved in the biosynthesis of terpenoids of which 35 contained terpene synthase motifs. This study also revealed 32,341 simple sequence repeats (SSRs) in 23,168 transcripts. Assembled sequences of transcriptome of A. paniculata generated in this study are made available, for the first time, in the TSA database, which provides useful information for functional and comparative genomic analysis besides identification of key enzymes involved in the various pathways of secondary metabolism.

Authentication of medicinal plants by DNA markers

Medicinal plants have been used worldwide for centuries to maintain health and to treat diseases, more so chronic diseases. However, adulteration and use of spurious materials as substitutes have become a major concern for users and industry for reasons of safety and efficacy. Therefore, authentication of medicinal plants is of utmost importance. Morphological, anatomical, chemical and DNA markers solve the problem by differentiating the genuine material from the adulterants, substitutes and spurious drugs. DNA markers use nucleotide sequences to identify species; it takes preference over the other two markers being not age dependent, tissue specific and having a higher discriminating power. Therefore, characterization of plants with such markers is an ideal approach for identification of medicinal plant species and populations/varieties of the same species. Availability of certified taxonomic specimens in herbaria is certainly required for unambiguous confirmation through final visual comparison and analysis.

De Novo Transcriptome Assembly in Firmiana danxiaensis, a Tree Species Endemic to the Danxia Landform

Many Firmiana species are locally endemic, providing an interesting system for studying adaptation and speciation. Among these species, F. danxiaensis is a tree species endemic to Mount Danxia in Guangdong, China, which is an area known for presenting the Danxia landform. How F. danxiaensis could have adapted to the stressful environment of rocky cliffs covered with barren soils in the Danxia landform is still unknown. In this study, we performed de novo assembly of the transcriptome of F. danxiaensis, obtaining 47,221 unigenes with an N50 value of 987 bp. Homology analysis showed that 32,318 of the unigenes presented hits in the NCBI non-redundant database, and 31,857 exhibited significant matches with the protein database of Theobroma cacao. Gene Ontology (GO) annotation showed that hundreds of unigenes participated in responses to various stresses or nutritional starvation, which may help us to understand the adaptation of F. danxiaensis to Danxia landform. Additionally, we found 263 genes related to responses to Cd, partially explaining the high accumulation of Cd observed in Firmiana species. The EuKaryotic Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations revealed many genes playing roles in the biosynthesis of secondary metabolites and environmental adaptation, which may also contribute to the survivor and success of Firmiana species in extreme environments. Based on the obtained transcriptome, we further identified a Firmiana-specific whole-genome duplication event that occurred approximately 20 Mya, which may have provided raw materials for the diversification of Firmiana species.

De Novo Assembly and Characterization of Tall Fescue Transcriptome under Water Stress

Water stress is a fundamental problem for tall fescue [Lolium arundinaceum (Schreb.) Darbysh.] cultivation in the south-central United States. Genetic improvement of tall fescue for water-stress tolerance is the key strategy for improving its persistence in the region. Genotypes with contrasting characteristics for relative water content and osmotic potential were identified from a tall fescue population. Transcriptome profiling between water-stress-tolerant (B400) and water-stress-susceptible (W279) genotypes was performed to unravel the genetic regulatory mechanism of water-stress responses in tall fescue. RNA samples from leaf, shoot, root, and inflorescence were pooled and sequenced through Illumina paired-end sequencing. A total of 199,399 contigs were assembled with an average length of 585 bp. Between the two genotypes, 2986 reference transcripts (RTs) were significantly differentially expressed and 1048 of them could be annotated and found to associate with metabolic pathways and enzyme coding genes. In total, 175 differentially expressed RTs were reported for various stress-related functions. Among those, 65 encoded kinase proteins, 40 each encoded transposons, and transporter proteins were previously reported to be involved with abiotic stress responses. A total of 6348 simple sequence repeats and 6658 single-nucleotide polymorphisms were identified in the contig sequences. Primers were developed from the corresponding sequences, which might be used as candidate gene markers in tall fescue. This study will lead to identification of genes or transcription factors related to water-stress tolerance and development of a comprehensive molecular marker system to facilitate marker-assisted breeding in tall fescue.

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.

Transcriptome profiling shows gene regulation patterns in a flavonoid pathway in response to exogenous phenylalanine in Boesenbergia rotunda cell culture

Background

Panduratin A extracted from Boesenbergia rotunda is a flavonoid reported to possess a range of medicinal indications which include anti-dengue, anti-HIV, anti-cancer, antioxidant and anti-inflammatory properties. Boesenbergia rotunda is a plant from the Zingiberaceae family commonly used as a food ingredient and traditional medicine in Southeast Asia and China. Reports on the health benefits of secondary metabolites extracted from Boesenbergia rotunda over the last few years has resulted in rising demands for panduratin A. However large scale extraction has been hindered by the naturally low abundance of the compound and limited knowledge of its biosynthetic pathway.

Results

Transcriptome sequencing and digital gene expression (DGE) analysis of native and phenylalanine treated Boesenbergia rotunda cell suspension cultures were carried out to elucidate the key genes differentially expressed in the panduratin A biosynthetic pathway. Based on experiments that show increase in panduratin A production after 14 days post treatment with exogenous phenylalanine, an aromatic amino acid derived from the shikimic acid pathway, total RNA of untreated and 14 days post-phenylalanine treated cell suspension cultures were extracted and sequenced using next generation sequencing technology employing an Illumina-Solexa platform. The transcriptome data generated 101, 043 unigenes with 50, 932 (50.41%) successfully annotated in the public protein databases; including 49.93% (50, 447) in the non-redundant (NR) database, 34.63% (34, 989) in Swiss-Prot, 24,07% (24, 316) in Kyoto Encyclopedia of Genes and Genomes (KEGG) and 16.26% (16, 426) in Clusters of Orthologous Groups (COG). Through DGE analysis, we found that 14, 644 unigenes were up-regulated and 14, 379 unigenes down-regulated in response to exogenous phenylalanine treatment. In the phenylpropanoid pathway leading to the proposed panduratin A production, 2 up-regulated phenylalanine ammonia-lyase (PAL), 3 up-regulated 4-coumaroyl:coenzyme A ligase (4CL) and 1 up-regulated chalcone synthase (CHS) were found.

Conclusions

This is the first report of Boesenbergia rotunda de novo transcriptome data that could serve as a reference for gene or enzyme functional studies in the Zingiberaceae family. Although enzymes that are directly involved in the panduratin A biosynthetic pathway were not completely elucidated, the data provides an overall picture of gene regulation patterns leading to panduratin A production.

Electronic supplementary material

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

Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding

The transcript pool of a plant part, under any given condition, is a collection of mRNAs that will pave the way for a biochemical reaction of the plant to stimuli. Over the past decades, transcriptome study has advanced from Northern blotting to RNA sequencing (RNA-seq), through other techniques, of which real-time quantitative polymerase chain reaction (PCR) and microarray are the most significant ones. The questions being addressed by such studies have also matured from a solitary process to expression atlas and marker-assisted genetic enhancement. Not only genes and their networks involved in various developmental processes of plant parts have been elucidated, but also stress tolerant genes have been highlighted. The transcriptome of a plant with altered expression of a target gene has given information about the downstream genes. Marker information has been used for breeding improved varieties. Fortunately, the data generated by transcriptome analysis has been made freely available for ample utilization and comparison. The review discusses this wide variety of transcriptome data being generated in plants, which includes developmental stages, abiotic and biotic stress, effect of altered gene expression, as well as comparative transcriptomics, with a special emphasis on microarray and RNA-seq. Such data can be used to determine the regulatory gene networks, which can subsequently be utilized for generating improved plant varieties.

Exploring the genes of yerba mate (Ilex paraguariensis A. St.-Hil.) by NGS and de novo transcriptome assembly

Yerba mate (Ilex paraguariensis A. St.-Hil.) is an important subtropical tree crop cultivated on 326,000 ha in Argentina, Brazil and Paraguay, with a total yield production of more than 1,000,000 t. Yerba mate presents a strong limitation regarding sequence information. The NCBI GenBank lacks an EST database of yerba mate and depicts only 80 DNA sequences, mostly uncharacterized. In this scenario, in order to elucidate the yerba mate gene landscape by means of NGS, we explored and discovered a vast collection of I. paraguariensis transcripts. Total RNA from I. paraguariensis was sequenced by Illumina HiSeq-2000 obtaining 72,031,388 pair-end 100 bp sequences. High quality reads were de novo assembled into 44,907 transcripts encompassing 40 million bases with an estimated coverage of 180X. Multiple sequence analysis allowed us to predict that yerba mate contains ∼ 32,355 genes and 12,551 gene variants or isoforms. We identified and categorized members of more than 100 metabolic pathways. Overall, we have identified ∼ 1,000 putative transcription factors, genes involved in heat and oxidative stress, pathogen response, as well as disease resistance and hormone response. We have also identified, based in sequence homology searches, novel transcripts related to osmotic, drought, salinity and cold stress, senescence and early flowering. We have also pinpointed several members of the gene silencing pathway, and characterized the silencing effector Argonaute1. We predicted a diverse supply of putative microRNA precursors involved in developmental processes. We present here the first draft of the transcribed genomes of the yerba mate chloroplast and mitochondrion. The putative sequence and predicted structure of the caffeine synthase of yerba mate is presented. Moreover, we provide a collection of over 10,800 SSR accessible to the scientific community interested in yerba mate genetic improvement. This contribution broadly expands the limited knowledge of yerba mate genes, and is presented as the first genomic resource of this important crop.

De novo assembly of the perennial ryegrass transcriptome using an RNA-Seq strategy

BACKGROUND

Perennial ryegrass is a highly heterozygous outbreeding grass species used for turf and forage production. Heterozygosity can affect de-Bruijn graph assembly making de novo transcriptome assembly of species such as perennial ryegrass challenging. Creating a reference transcriptome from a homozygous perennial ryegrass genotype can circumvent the challenge of heterozygosity. The goals of this study were to perform RNA-sequencing on multiple tissues from a highly inbred genotype to develop a reference transcriptome. This was complemented with RNA-sequencing of a highly heterozygous genotype for SNP calling.

RESULT

De novo transcriptome assembly of the inbred genotype created 185,833 transcripts with an average length of 830 base pairs. Within the inbred reference transcriptome 78,560 predicted open reading frames were found of which 24,434 were predicted as complete. Functional annotation found 50,890 transcripts with a BLASTp hit from the Swiss-Prot non-redundant database, 58,941 transcripts with a Pfam protein domain and 1,151 transcripts encoding putative secreted peptides. To evaluate the reference transcriptome we targeted the high-affinity K+ transporter gene family and found multiple orthologs. Using the longest unique open reading frames as the reference sequence, 64,242 single nucleotide polymorphisms were found. One thousand sixty one open reading frames from the inbred genotype contained heterozygous sites, confirming the high degree of homozygosity.

CONCLUSION

Our study has developed an annotated, comprehensive transcriptome reference for perennial ryegrass that can aid in determining genetic variation, expression analysis, genome annotation, and gene mapping.

De novo sequencing and comparative analysis of holy and sweet basil transcriptomes

Background

Ocimum L. of family Lamiaceae is a well known genus for its ethnobotanical, medicinal and aromatic properties, which are attributed to innumerable phenylpropanoid and terpenoid compounds produced by the plant. To enrich genomic resources for understanding various pathways, de novo transcriptome sequencing of two important species, O. sanctum and O. basilicum, was carried out by Illumina paired-end sequencing.

Results

The sequence assembly resulted in 69117 and 130043 transcripts with an average length of 1646 ± 1210.1 bp and 1363 ± 1139.3 bp for O. sanctum and O. basilicum, respectively. Out of the total transcripts, 59648 (86.30%) and 105470 (81.10%) from O. sanctum and O. basilicum, and respectively were annotated by uniprot blastx against Arabidopsis, rice and lamiaceae. KEGG analysis identified 501 and 952 transcripts from O. sanctum and O. basilicum, respectively, related to secondary metabolism with higher percentage of transcripts for biosynthesis of terpenoids in O. sanctum and phenylpropanoids in O. basilicum. Higher digital gene expression in O. basilicum was validated through qPCR and correlated to higher essential oil content and chromosome number (O. sanctum, 2n = 16; and O. basilicum, 2n = 48). Several CYP450 (26) and TF (40) families were identified having probable roles in primary and secondary metabolism. Also SSR and SNP markers were identified in the transcriptomes of both species with many SSRs linked to phenylpropanoid and terpenoid pathway genes.

Conclusion

This is the first report of a comparative transcriptome analysis of Ocimum species and can be utilized to characterize genes related to secondary metabolism, their regulation, and breeding special chemotypes with unique essential oil composition in Ocimum.

Electronic supplementary material

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