De novo assembly and characterization of Sophora japonica transcriptome using RNA-seq

Comparative Transcriptome Analysis of Sophora japonica (L.) Roots Reveals Key Pathways and Genes in Response to PEG-Induced Drought Stress under Different Nitrogen Conditions

Sophora japonica is a native leguminous tree species in China. The high stress tolerance contributes to its long lifespan of thousands of years. The lack of genomic resources greatly limits genetic studies on the stress responses of S. japonica. In this study, RNA-seq was conducted for S. japonica roots grown under short-term 20% polyethylene glycol (PEG) 6000-induced drought stress under normal N and N starvation conditions (1 and 0 mM NH4NO3, respectively). In each of the libraries, we generated more than 25 million clean reads, which were then de novo assembled to 46,852 unigenes with an average length of 1310.49 bp. In the differential expression analyses, more differentially expressed genes (DEGs) were found under drought with N starvation than under single stresses. The number of transcripts identified under N starvation and drought in S. japonica was nearly the same, but more upregulated genes were induced by drought, while more downregulated genes were induced by N starvation. Genes involved in “phenylpropanoid biosynthesis” and “biosynthesis of amino acids” pathways were upregulated according to KEGG enrichment analyses, irrespective of the stress treatments. Additionally, upregulated N metabolism genes were enriched upon drought, and downregulated photosynthesis genes were enriched under N starvation. We found 4,372 and 5,430 drought-responsive DEGs under normal N and N starvation conditions, respectively. N starvation may aggravate drought by downregulating transcripts in the “carbon metabolism”, “ribosome”, “arginine biosynthesis pathway”, “oxidative phosphorylation” and “aminoacyl-tRNA biosynthesis” pathways. We identified 78 genes related to N uptake and assimilation, 38 of which exhibited differential expression under stress. A total of 395 DEGs were categorized as transcription factors, of which AR2/ERF-ERF, WRKY, NAC, MYB, bHLH, C3H and C2C2-Dof families played key roles in drought and N starvation stresses. The transcriptome data obtained, and the genes identified facilitate our understanding of the mechanisms of S. japonica responses to drought and N starvation stresses and provide a molecular foundation for understanding the mechanisms of its long lifespan for breeding resistant varieties for greening.

RNA-Seq De Novo Assembly and Differential Transcriptome Analysis of Korean Medicinal Herb Cirsium japonicum var. spinossimum

Cirsium japonicum belongs to the Asteraceae or Compositae family and is a medicinal plant in Asia that has a variety of effects, including tumour inhibition, improved immunity with flavones, and antidiabetic and hepatoprotective effects. Silymarin is synthesized by 4-coumaroyl-CoA via both the flavonoid and phenylpropanoid pathways to produce the immediate precursors taxifolin and coniferyl alcohol. Then, the oxidative radicalization of taxifolin and coniferyl alcohol produces silymarin. We identified the expression of genes related to the synthesis of silymarin in C. japonicum in three different tissues, namely, flowers, leaves, and roots, through RNA sequencing. We obtained 51,133 unigenes from transcriptome sequencing by de novo assembly using Trinity v2.1.1, TransDecoder v2.0.1, and CD-HIT v4.6 software. The differentially expressed gene analysis revealed that the expression of genes related to the flavonoid pathway was higher in the flowers, whereas the phenylpropanoid pathway was more highly expressed in the roots. In this study, we established a global transcriptome dataset for C. japonicum. The data shall not only be useful to focus more deeply on the genes related to product medicinal metabolite including flavolignan but also to study the functional genomics for genetic engineering of C. japonicum.

Identification of cytochrome P450 monooxygenase genes and their expression in response to high temperature in the alligatorweed flea beetle Agasicles hygrophila (Coleoptera: Chrysomelidae)

Cytochrome P450 monooxygenases (P450s) are a large class of enzymes that play essential roles in metabolic processes such as hormone synthesis and the catabolism of toxins and other chemicals in insects. In the present study, we identified 82 P450 genes using comprehensive RNA sequencing in the flea beetle Agasicles hygrophila, and all of the sequences were validated by cloning and sequencing. Phylogenetic analysis showed that the P450 genes in A. hygrophila fell into the mitochondrial clan, CYP2 clan, CYP3 clan and CYP4 clan and were classified into 20 families and 48 subfamilies. Most A. hygrophila P450 genes had high sequence homology with those from other coleopteran insects. To understand the effects of high temperatures on the metabolic processes of female and male adults, we studied the effects of two temperature regimes (constant temperature of 28 °C for 20 h with a 4-h period of high temperatures of 30 °C and 39 °C) on the expression levels of P450 genes in A. hygrophila using RT-PCR and qRT-PCR. The results showed that there were no differences in expression in 30 P450 genes between the control and high-temperature-treated A. hygrophila adults, while 22 P450 genes showed up-regulated expression and 19 P450 genes were down-regulated in A. hygrophila female adults after high-temperature treatment. For A. hygrophila male adults exposed to high temperatures, we found that 8 P450 genes had higher expression levels and 12 P450 genes had lower expression levels under the same conditions. The P450 genes are candidates that showed significantly different expression levels after high-temperature treatments in A. hygrophila adults, and further studies are needed to determine their possible roles in metabolic processes during the response to elevated temperatures.

Comparative transcriptome analysis to identify putative genes involved in thymol biosynthesis pathway in medicinal plant Trachyspermum ammi L

Thymol, as a dietary monoterpene, is a phenol derivative of cymene, which is the major component of the essential oil of Trachyspermum ammi (L.). It shows multiple biological activities: antifungal, antibacterial, antivirus and anti-inflammatory. T. ammi, commonly known as ajowan, belongs to Apiaceae and is an important medicinal seed spice. To identify the putative genes involved in thymol and other monoterpene biosynthesis, we provided transcriptomes of four inflorescence tissues of two ajowan ecotypes, containing different thymol yield. This study has detected the genes encoding enzymes for the go-between stages of the terpenoid biosynthesis pathways. A large number of unigenes, differentially expressed between four inflorescence tissues of two ajowan ecotypes, was revealed by a transcriptome analysis. Furthermore, differentially expressed unigenes encoding dehydrogenases, transcription factors, and cytochrome P450s, which might be associated with terpenoid diversity in T. ammi, were identified. The sequencing data obtained in this study formed a valuable repository of genetic information for an understanding of the formation of the main constituents of ajowan essential oil and functional analysis of thymol-specific genes. Comparative transcriptome analysis led to the development of new resources for a functional breeding of ajowan.

Characterization and analysis of a de novo transcriptome from the pygmy grasshopper Tetrix japonica

The pygmy grasshopper Tetrix japonica is a common insect distributed throughout the world, and it has the potential for use in studies of body colour polymorphism, genomics and the biology of Tetrigoidea (Insecta: Orthoptera). However, limited biological information is available for this insect. Here, we conducted a de novo transcriptome study of adult and larval T. japonica to provide a better understanding of its gene expression and develop genomic resources for future work. We sequenced and explored the characteristics of the de novo transcriptome of T. japonica using Illumina HiSeq 2000 platform. A total of 107 608 206 paired-end clean reads were assembled into 61 141 unigenes using the trinity software; the mean unigene size was 771 bp, and the N50 length was 1238 bp. A total of 29 225 unigenes were functionally annotated to the NCBI nonredundant protein sequences (Nr), NCBI nonredundant nucleotide sequences (Nt), a manually annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of putative genes that are potentially involved in pigment pathways, juvenile hormone (JH) metabolism and signalling pathways were identified in the T. japonica transcriptome. Additionally, 165 769 and 156 796 putative single nucleotide polymorphisms occurred in the adult and larvae transcriptomes, respectively, and a total of 3162 simple sequence repeats were detected in this assembly. This comprehensive transcriptomic data for T. japonica will provide a usable resource for gene predictions, signalling pathway investigations and molecular marker development for this species and other pygmy grasshoppers.

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.

De Novo Transcriptome Analysis of Medicinally Important Plantago ovata Using RNA-Seq

Plantago ovata is an economically and medicinally important plant of the family Plantaginaceae. It is used extensively for the production of seed husk for its application in pharmaceutical, food and cosmetic industries. In the present study, the transcriptome of P. ovata ovary was sequenced using Illumina Genome Analyzer platform to characterize the mucilage biosynthesis pathway in the plant. De novo assembly was carried out using Oases followed by velvet. A total of 46,955 non-redundant transcripts (≥100 bp) using ~29 million high-quality paired end reads were generated. Functional categorization of these transcripts revealed the presence of several genes involved in various biological processes like metabolic pathways, mucilage biosynthesis, biosynthesis of secondary metabolites and antioxidants. In addition, simple sequence-repeat motifs, non-coding RNAs and transcription factors were also identified. Expression profiling of some genes involved in mucilage biosynthetic pathway was performed in different tissues of P. ovata using Real time PCR analysis. The study has resulted in a valuable resource for further studies on gene expression, genomics and functional genomics in P. ovata.

De novo assembly and characterisation of the field pea transcriptome using RNA-Seq

BACKGROUND

Field pea (Pisum sativum L.) is a cool-season grain legume that is cultivated world-wide for both human consumption and stock-feed purposes. Enhancement of genetic and genomic resources for field pea will permit improved understanding of the control of traits relevant to crop productivity and quality. Advances in second-generation sequencing and associated bioinformatics analysis now provide unprecedented opportunities for the development of such resources. The objective of this study was to perform transcriptome sequencing and characterisation from two genotypes of field pea that differ in terms of seed and plant morphological characteristics.

RESULTS

Transcriptome sequencing was performed with RNA templates from multiple tissues of the field pea genotypes Kaspa and Parafield. Tissue samples were collected at various growth stages, and a total of 23 cDNA libraries were sequenced using Illumina high-throughput sequencing platforms. A total of 407 and 352 million paired-end reads from the Kaspa and Parafield transcriptomes, respectively were assembled into 129,282 and 149,272 contigs, which were filtered on the basis of known gene annotations, presence of open reading frames (ORFs), reciprocal matches and degree of coverage. Totals of 126,335 contigs from Kaspa and 145,730 from Parafield were subsequently selected as the reference set. Reciprocal sequence analysis revealed that c. 87% of contigs were expressed in both cultivars, while a small proportion were unique to each genotype. Reads from different libraries were aligned to the genotype-specific assemblies in order to identify and characterise expression of contigs on a tissue-specific basis, of which 87% were expressed in more than one tissue, while others showed distinct expression patterns in specific tissues, providing unique transcriptome signatures.

CONCLUSION

This study provided a comprehensive assembled and annotated transcriptome set for field pea that can be used for development of genetic markers, in order to assess genetic diversity, construct linkage maps, perform trait-dissection and implement whole-genome selection strategies in varietal improvement programs, as well to identify target genes for genetic modification approaches on the basis of annotation and expression analysis. In addition, the reference field pea transcriptome will prove highly valuable for comparative genomics studies and construction of a finalised genome sequence.

De Novo Assembly and Annotation of the Chinese Chive (Allium tuberosum Rottler ex Spr.) Transcriptome Using the Illumina Platform

Chinese chive (A. tuberosum Rottler ex Spr.) is one of the most widely cultivated Allium species in China. However, minimal transcriptomic and genomic data are available to reveal its evolution and genetic diversity. In this study, de novo transcriptome sequencing was performed to produce large transcript sequences using an Illumina HiSeq 2000 instrument. We produced 51,968,882 high-quality clean reads and assembled them into 150,154 contigs. A total of 60,031 unigenes with an average length of 631 bp were identified. Of these, 36,523 unigenes were homologous to existing database sequences, 35,648 unigenes were annotated in the NCBI non-redundant (Nr) sequence database, and 23,509 unigenes were annotated in the Swiss-Prot database. A total of 26,798 unigenes were assigned to 57 Gene Ontology (GO) terms, and 13,378 unigenes were assigned to Cluster of Orthologous Group categories. Using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, we mapped 21,361 unigenes onto 128 pathways. Furthermore, 2,125 sequences containing simple sequence repeats (SSRs) were identified. This new dataset provides the most comprehensive resource currently available for gene expression, gene discovery, and future genomic research on Chinese chive. The sequence resources developed in this study can be used to develop molecular markers that will facilitate further genetic research on Chinese chive and related species.

RNA-Seq Based De Novo Transcriptome Assembly and Gene Discovery of Cistanche deserticola Fleshy Stem

Backgrounds

Cistanche deserticola is a completely non-photosynthetic parasitic plant with great medicinal value and mainly distributed in desert of Northwest China. Its dried fleshy stem is a crucial tonic in traditional Chinese medicine with roles of mainly improving male sexual function and strengthening immunity, but few mechanistic studies have been conducted partly due to the lack of genomic and transcriptomic resources.

Results

In this study, we performed deep transcriptome sequencing in fleshy stem of C. deserticola, and about 80 million reads were generated using Illumina pair-end sequencing on HiSeq2000 platform. Using trinity assembler, we obtained 95,787 transcript sequences with transcript lengths ranging from 200bp to 15,698bp, having an average length of 950 bases and the N50 length of 1,519 bases. 63,957 transcripts were identified actively expressed with FPKM ≥ 0.5, in which 30,098 transcripts were annotated with gene descriptions or gene ontology terms by sequence similarity analyses against several public databases (Uniprot, NR and Nt at NCBI, and KEGG). Furthermore, we identified key enzyme genes involved in biosynthesis of lignin and phenylethanoid glycosides (PhGs) which are known to be the primary active ingredients. Four phenylalanine ammonia-lyase (PAL) genes, the first key enzyme in lignin and PhG biosynthesis, were identified based on sequences comparison and phylogenetic analysis. Two biosynthesis pathways of PhGs were also proposed for the first time.

Conclusions

In all, we completed a global analysis of the C. deserticola fleshy stem transcriptome using RNA-seq technology. A collection of enzyme genes related to biosynthesis of lignin and phenylethanoid glysides were identified from the assembled and annotated transcripts, and the gene family of PAL was also predicted. The sequence data from this study will provide a valuable resource for conducting future phenylethanoid glysides biosynthesis researches and functional genomic studies in this important medicinal plant.

Rapid transcriptome sequencing of an invasive pest, the brown marmorated stink bug Halyomorpha halys

Background

Halyomorpha halys (Stål) (Insecta:Hemiptera;Pentatomidae), commonly known as the Brown Marmorated Stink Bug (BMSB), is an invasive pest of the mid-Atlantic region of the United States, causing economically important damage to a wide range of crops. Native to Asia, BMSB was first observed in Allentown, PA, USA, in 1996, and this pest is now well-established throughout the US mid-Atlantic region and beyond. In addition to the serious threat BMSB poses to agriculture, BMSB has become a nuisance to homeowners, invading home gardens and congregating in large numbers in human-made structures, including homes, to overwinter. Despite its significance as an agricultural pest with limited control options, only 100 bp of BMSB sequence data was available in public databases when this project began.

Results

Transcriptome sequencing was undertaken to provide a molecular resource to the research community to inform the development of pest control strategies and to provide molecular data for population genetics studies of BMSB. Using normalized, strand-specific libraries, we sequenced pools of all BMSB life stages on the Illumina HiSeq. Trinity was used to assemble 200,000 putative transcripts in >100,000 components. A novel bioinformatic method that analyzed the strand-specificity of the data reduced this to 53,071 putative transcripts from 18,573 components. By integrating multiple other data types, we narrowed this further to 13,211 representative transcripts.

Conclusions

Bacterial endosymbiont genes were identified in this dataset, some of which have a copy number consistent with being lateral gene transfers between endosymbiont genomes and Hemiptera, including ankyrin-repeat related proteins, lysozyme, and mannanase. Such genes and endosymbionts may provide novel targets for BMSB-specific biocontrol. This study demonstrates the utility of strand-specific sequencing in generating shotgun transcriptomes and that rapid sequencing shotgun transcriptomes is possible without the need for extensive inbreeding to generate homozygous lines. Such sequencing can provide a rapid response to pest invasions similar to that already described for disease epidemiology.

Electronic supplementary material

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