Insights into Taxol® biosynthesis by endophytic fungi

There have been two hundred reports that endophytic fungi produce Taxol®, but its production yield is often rather low. Although considerable efforts have been made to increase Taxol/taxanes production in fungi by manipulating cocultures, mutagenesis, genome shuffles, and gene overexpression, little is known about the molecular signatures of Taxol biosynthesis and its regulation. It is known that some fungi have orthologs of the Taxol biosynthetic pathway, but the overall architecture of this pathway is unknown. A biosynthetic putative gene homology approach, combined with genomics and transcriptomics analysis, revealed that a few genes for metabolite residues may be located on dispensable chromosomes. This review explores a number of crucial topics (i) finding biosynthetic pathway genes using precursors, elicitors, and inhibitors; (ii) orthologs of the Taxol biosynthetic pathway for rate-limiting genes/enzymes; and (iii) genomics and transcriptomics can be used to accurately predict biosynthetic putative genes and regulators. This provides promising targets for future genetic engineering approaches to produce fungal Taxol and precursors. KEY POINTS: • A recent trend in predicting Taxol biosynthetic pathway from endophytic fungi. • Understanding the Taxol biosynthetic pathway and related enzymes in fungi. • The genetic evidence and formation of taxane from endophytic fungi.

The first transcriptome sequencing and data analysis of the Javan mahseer (Tor tambra)

The Javan mahseer (Tor tambra) is one of the most valuable freshwater fish found in Tor species. To date, other than mitogenomic data (BioProject: PRJNA422829), genomic and transcriptomic resources for this species are still lacking which is crucial to understand the molecular mechanisms associated with important traits such as growth, immune response, reproduction and sex determination. For the first time, we sequenced the transcriptome from a whole juvenile fish using Illumina NovaSEQ6000 generating raw paired-end reads. De novo transcriptome assembly generated a draft transcriptome (BUSCO5 completeness of 91.2% [Actinopterygii_odb10 database]) consisting of 259,403 putative transcripts with a total and N50 length of 333,881,215 bp and 2283 bp, respectively. A total count of 77,503 non-redundant protein coding sequences were predicted from the transcripts and used for functional annotation. We mapped the predicted proteins to 304 known KEGG pathways with signal transduction cluster having the highest representation followed by immune system and endocrine system. In addition, transcripts exhibiting significant similarity to previously published growth-and immune-related genes were identified which will facilitate future molecular breeding of Tor tambra.

Prospects of next generation sequencing in lentil breeding

Lentil is an important food legume crop that has large and complex genome. During past years, considerable attention has been given on the use of next generation sequencing for enriching the genomic resources including identification of SSR and SNP markers, development of unigenes, transcripts, and identification of candidate genes for biotic and abiotic stresses, analysis of genetic diversity and identification of genes/ QTLs for agronomically important traits. However, in other crops including pulses, next generation sequencing has revolutionized the genomic research and helped in genomic assisted breeding rapidly and cost effectively. The present review discuss current status and future prospects of the use NGS based breeding in lentil.

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.

Global gene expression and pigment analysis of two contrasting flower color cultivars of Canna

Development of flower color in plants is a complex process. Among others, it is an important trait for ornamental flowering plants. Canna is a flowering ornamental plant of family Cannaceae. To understand the molecular mechanism of flower color development in Canna, RNA sequencing from flower tissues of two contrasting flower color cultivars, Red President (RP) and Tropical Sunrise (TS) was performed. More than 27.0 million and 19.0 million clean reads were obtained from RP and TS, respectively. The combined clean reads were assembled into 147,295 unigenes. The Canna unigenes showed maximum homology with Populus trichocarpa (26.79%). A total of 2702 unigenes expressed differentially between the two cultivars of which 1972 were up-regulated and 730 were down-regulated in RP. Phenylpropanoid and flavonoid biosynthetic processes were the significant processes in RP. Expression of a vast number of transcription factors including MYB, bHLH, ARF, and WRKY were higher in RP than TS. The expression analysis of RNA sequencing data was validated by qRT-PCR analysis. Further, concentration of measured anthocyanidins and flavonols were very low or absent in TS, corroborating largely with our transcriptome data. These findings may help in understanding flower color development in Canna and in future crop breeding program.

Transcript profiling and gene expression analysis under drought stress in Ziziphus nummularia (Burm.f.) Wright & Arn

Drought is one of the prime abiotic stresses responsible for limiting agricultural productivity. A number of drought responsive genes have been isolated and functionally characterized but these studies have been restricted to a few model plant systems. Very few drought responsive genes have been reported till date from non model drought tolerant plants. The present study aimed at identifying differentially expressed genes from a drought tolerant, non-model plant, Ziziphus nummularia (Burm.f.) Wight & Arn. One month old seedlings of Z. nummularia were subjected to drought stress by 30% Polyethylene glycol (PEG 6000) treatment for 6, 12, 24, 48 and 72 h. A significant reduction in RWC and increase in proline was observed at 24 h and 48 h of treatment. Suppression subtractive hybridization (SSH) library was constructed with drought stressed seedlings after 24 h and 48 h of PEG 6000 treatment. A total of 142 and 530 unigenes from 24 h and 48 h library were identified respectively. Gene ontology studies revealed that about 9.78% and 15.07% unigenes from 24 h and 48 h SSH libraries were expressed in "response to stress". Fifteen putative drought responsive genes identified in SSH library were validated for drought responsive differential expression by RT-qPCR. Significant changes in fold expressions were observed with time in the treated samples compared to the control. A heat map revealing the expression profile of genes was constructed by hierarchical clustering. Various genes identified in SSH libraries can serve as a resource for marker discovery and selection of candidate genes to improve drought tolerance in other susceptible crops.

De novo transcriptome analysis of Rhododendron molle G. Don flowers by Illumina sequencing

Rhododendron molle G. Don occupies an important phylogenetic node in the genus rhododendron with unique yellow flower and medicinal functions. However, only limited genetic resources and their genome information are available for the generation of rhododendron flowers. The next generation sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provide a turning point for rhododendron research. Our goal is to use the genetic information to facilitate the relevant research on flowering and flower color formation in R. molle. In total, 66,026 unigenes were identified, among which 31,298 were annotated in the NCBI non-redundant protein database and 22,410 were annotated in the Swiss-Prot database. Of these annotated unigenes, 9490 and 18,680 unigenes were assigned to clusters of orthologous groups and gene ontology categories, respectively. A total of 7177 genes were mapped to 118 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. In addition, 8266 simple sequence repeats (SSRs) were detected, and these SSRs will undoubtedly benefit rhododendron breeding work. Metabolic pathway analysis revealed that 32 unigenes were predicted to be involved in carotenoid biosynthesis. Our transcriptome revealed 32 engines that encode key enzymes in the carotenoid biosynthesis pathway, including PSY, PDS, LCYB, LCYE, etc. The content of β-carotene was much higher than the other carotenoids throughout the flower development. It was consistent with the key genes expression level in the carotenoid biosynthesis pathway by the Illumina expression profile analysis and the qRT-PCR analysis. Our study identified genes associated with carotenoid biosynthesis in R. molle and provides a valuable resource for understanding the flowering and flower color formation mechanisms in R. molle.

Comparative Transcriptome Analysis Provides Insights into Differentially Expressed Genes and Long Non-Coding RNAs between Ovary and Testis of the Mud Crab (Scylla paramamosain)

The molecular mechanism underlying sex determination and gonadal differentiation of the mud crab (Scylla paramamosain) has received considerable attention, due to the remarkably biological and economic differences between sexes. However, sex-biased genes, especially non-coding genes, which account for these differences, remain elusive in this crustacean species. In this study, the first de novo gonad transcriptome sequencing was performed to identify both differentially expressed genes and long non-coding RNAs (lncRNAs) between male and female S. paramamosain by using Illumina Hiseq2500. A total of 79,282,758 and 79,854,234 reads were generated from ovarian and testicular cDNA libraries, respectively. After filtrating and de novo assembly, 262,688 unigenes were produced from both libraries. Of these unigenes, 41,125 were annotated with known protein sequences in public databases. Homologous genes involved in sex determination and gonadal development pathways (Sxl-Tra/Tra-2-Dsx/Fru, Wnt4, thyroid hormone synthesis pathway, etc.) were identified. Three hundred and sixteen differentially expressed unigenes were further identified between both transcriptomes. Meanwhile, a total of 233,078 putative lncRNAs were predicted. Of these lncRNAs, 147 were differentially expressed between sexes. qRT-PCR results showed that nine lncRNAs negatively regulated the expression of eight genes, suggesting a potential role in sex differentiation. These findings will provide fundamental resources for further investigation on sex differentiation and regulatory mechanism in crustaceans.

Transcriptome analysis of the Chinese giant salamander (Andrias davidianus) using RNA-sequencing

The Chinese giant salamander (Andrias davidianus) is an economically important animal on academic value. However, the genomic information of this species has been less studied. In our study, the transcripts of A. davidianus were obtained by RNA-seq to conduct a transcriptomic analysis. In total 132,912 unigenes were generated with an average length of 690 bp and N50 of 1263 bp by de novo assembly using Trinity software. Using a sequence similarity search against the nine public databases (CDD, KOG, NR, NT, PFAM, Swiss-prot, TrEMBL, GO and KEGG databases), a total of 24,049, 18,406, 36,711, 15,858, 20,500, 27,515, 36,705, 28,879 and 10,958 unigenes were annotated in databases, respectively. Of these, 6323 unigenes were annotated in all database and 39,672 unigenes were annotated in at least one database. Blasted with KEGG pathway, 10,958 unigenes were annotated, and it was divided into 343 categories according to different pathways. In addition, we also identified 29,790 SSRs. This study provided a valuable resource for understanding transcriptomic information of A. davidianus and laid a foundation for further research on functional gene cloning, genomics, genetic diversity analysis and molecular marker exploitation in A. davidianus.

De novo transcriptome assembly for the lung of the ornamented pygmy frog (Microhyla fissipes)

Microhyla fissipes, belonging to Neobatrachia, is a new model organism to study developmental biology, adaptive mechanisms from aquatic to terrestrial life, environmental toxicology, and human disease. M. fissipes use of lungs soon after hatching makes it extremely valuable for the study of lung function and development mechanisms. However, our knowledge of genes and pathways associated with lung development in M. fissipes is very limited. In this study, we conducted de novo transcriptome assembly for the lung of M. fissipes using the Illumina HiSeq4000 platform. We obtained approximately 9.0 GB clean data from the lung of the stage 28 tadpole with lung inflation (NCBI accession numbers: SRP107055). De novo transcriptome assembly identified 209,358 transcripts and 93,813 unigenes. In addition, BLASTX against NR, NT, KO, SwissProt, PFAM, GO and KOG databases were used to annotate all the 93,813 unigenes. This study provides the transcriptome and functional annotation of genes in M. fissipes lung development, which will be useful for comparative transcriptome analyses and promote research into mechanism of lung development in anuran.

Comparative transcriptome analysis reveals candidate genes for the biosynthesis of natural insecticide in Tanacetum cinerariifolium

Background

Pyrethrins are monoterpenoids and consist of either a chrysanthemic acid or pyrethric acid with a rethrolone moiety. Natural pyrethrins are safe and eco-friendly while possessing strong insecticidal properties. Despite such advantages of commercial value coming with the eco-friendly tag, most enzymes/genes involved in the pyrethrin biosynthesis pathway remain unidentified and uncharacterized. Since the flowers of Tanacetum cinerariifolium are rich in major pyrethrins, next generation transcriptome sequencing was undertaken to compare the flowers and the leaves of the plant de novo to identify differentially expressed transcripts and ascertain which among them might be involved in and responsible for the differential accumulation of pyrethrins in T. cinerariifolium flowers.

Results

In this first tissue specific transcriptome analysis of the non-model plant T. cinerariifolium, a total of 23,200,000 and 28,500,110 high quality Illumina next generation sequence reads, with a length of 101 bp, were generated for the flower and leaf tissue respectively. After functional enrichment analysis and GO based annotation using public protein databases such as UniRef, PFAM, SMART, KEGG and NR, 4443 and 8901 unigenes were identified in the flower and leaf tissue respectively. These could be assigned to 13344 KEGG pathways and the pyrethrin biosynthesis contextualized. The 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was involved in the biosynthesis of acid moiety of pyrethrin and this pathway predominated in the flowers as compared to the leaves. However, enzymes related to oxylipin biosynthesis were found predominantly in the leaf tissue, which suggested that major steps of pyrethrin biosynthesis occurred in the flowers.

Conclusions

Transcriptome comparison between the flower and leaf tissue of T. cinerariifolium provided an elaborate list of tissue specific transcripts that was useful in elucidating the differences in the expression of the biosynthetic pathways leading to differential presence of pyrethrin in the flowers. The information generated on genes, pathways and markers related to pyrethrin biosynthesis in this study will be helpful in enhancing the production of these useful compounds for value added breeding programs. Related proteome comparison to overlay our transcriptome comparison can generate more relevant information to better understand flower specific accumulation of secondary metabolites in general and pyrethrin accumulation in particular.

Electronic supplementary material

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

Identification of candidate genes involved in wax deposition in Poa pratensis by RNA-seq

Background

The cuticular wax plays important roles in plant resistance to various biotic and abiotic stresses. Understanding the synthesis and secretion of cuticular waxes is necessary in utilizing cuticular waxes to improve crop productivity and plant ecological adaptation. Due to the lack of genomic resources, little genetic research on cuticular wax deposition has been focused on Poa pratensis, a perennial forage and turf grass species that is widely distributed under various habitats. In this study, we performed de novo transcriptome sequencing to explore differentially expressed genes between the leaf non-elongation zone (NEZm) and the emerged blade zone (EBZ) and to identify genes related to cuticular wax deposition.

Results

A total of 77,707,414 high quality reads were obtained from llumina HiSeq 2500 platform, which were then assembled into 106,766 unigenes. Among them, 6019 unigenes showed significant differences in expression between NEZm and EBZ. In our assembled sequences, 3087 SSRs molecular markers were discovered. All the unigenes were searched against the NR, Swissprot, GO, COG, and KEGG databases using BLAST program for functional annotation. From 3156 unigenes with more expression in NEZm compared to EBZ, a number of unigenes involved in very long chain fatty acids (VLCFAs) and cuticular wax biosynthesis, transportation and regulation were identified. Several unigenes related to defense response and epidermal patterning were also found. Twelve putative genes involved in VLCFAs and cuticular wax biosynthesis were further analyzed for their expressions using qRT-PCR.

Conclusions

The transcriptome of P. pratensis leaf was deep sequenced, de novo assembled and annotated, and the candidate genes potentially involved in VLCFAs and cuticular wax biosynthesis, secretion and regulation in P. pratensis were identified. This provides fundamental genetic resources in improving plant adaptation to abiotic and biotic stresses.

Electronic supplementary material

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

RNA-Seq analysis for indigo biosynthesis pathway genes in Indigofera tinctoria and Polygonum tinctorium

Natural indigo is the most important blue dye for textile dyeing and valuable secondary metabolite biosynthesized in Indigofera tinctoria and Polygonum tinctorium plants. Present investigation is made to generation of gene resource for pathway enrichment and to understand possible gene expression involved in indigo biosynthesis. The data about raw reads and the transcriptome assembly project has been deposited at GenBank under the accessions SRA180766 and SRX692542 for I. tinctoria and P. tinctorium, respectively.

Analysis of Jatropha curcas transcriptome for oil enhancement and genic markers

Oil-rich seeds of Jatropha curcas are being focussed as a source of bio-diesel. However, prior to its industrial use, a lot of crop improvement efforts are required in Jatropha. Availability of a large number of EST sequences of Jatropha in public domain allow identification of candidate genes for several agronomic characters including oil content in seeds. Here, we have analysed 42,477 ESTs of Jatropha spanning 22.9 Mbp for microsatellites and fatty acid metabolism related sequences. Unigene sequences were built using CAP 3 programme resulted in 12,358 contigs and 5,730 singlets. Nearly, 8 % unigenes showed presence of microsatellites, slightly over-represented compared to their occurrence in ESTs. Most of the microsatellites were either di- or tri-nucleotide repeats, while other categories of tetra-, penta- and hexa-nucleotide repeats together constituted ~4 % of total microsatellites. Assessment of functional relevance of unigenes was carried out using Blast2GO using its default settings. The overall sequence similarity level against sequences in 'nr' database was >80 %. A total of 931 sequences that participated in any of the pathways related to fatty acid or lipid metabolism were found at GO level 6. Among these, GO terms "Fatty acid metabolic process" and "Fatty acid biosynthetic process" were most over-represented. Overall, our work has due relevance in identifying molecular markers for the candidate genes for oil content in Jatropha seeds, and will prove to be an important reference for further studies for identification of trait specific markers in Jatropha.

De novo transcriptome analysis of Liriodendron chinense petals and leaves by Illumina sequencing

Liriodendron chinense (Hemsl.) Sarg is an endangered species and occupies a pivotal position in phylogenetic studies of flowering plants, while its genomic resources are limited. In this study, we performed transcriptome sequencing for L. chinense petals and leaves using the Illumina paired-end sequencing technique. Approximately 17.02-Gb clean reads were obtained, and de novo assembly generated 87,841 unigenes, with an average length of 778 bp. Of these, there were 65,535 (74.61%) unigenes with significant similarity to publically available plant protein sequences. There were 3386 genes identified as significant differentially expressed between petals and leaves, among them 2969 (87.68%) were up-regulated and 417 (12.31%) down-regulated in petals. Metabolic pathway analysis revealed that 25 unigenes were predicted to be responsible for the biosynthesis of carotenoids, with 7 genes differentially expressed between these two tissues. This report is the first to identify genes associated with carotenoid biosynthesis in Liriodendron and represents a valuable resource for future genomic studies on the endangered species L. chinense.

Mining functional microsatellites in legume unigenes

Highly polymorphic and transferable microsatellites (SSRs) are important for comparative genomics, genome analysis and phylogenetic studies. Development of novel species-specific microsatellite markers remains a costly and labor-intensive project. Therefore, interest has been shifted from genomic to genic markers owing to their high inter-species transferability as they are developed from conserved coding regions of the genome. This study concentrates on comparative analysis of genic microsatellites in nine important legume (Arachis hypogaea, Cajanus cajan, Cicer arietinum, Glycine max, Lotus japonicus, Medicago truncatula, Phaseolus vulgaris, Pisum sativum and Vigna unguiculata) and two model plant species (Oryza sativa and Arabidopsis thaliana). Screening of a total of 228090 putative unique sequences spanning 219610522 bp using a microsatellite search tool, MISA, identified 12.18% of the unigenes containing 36248 microsatellite motifs excluding mononucleotide repeats. Frequency of legume unigene-derived SSRs was one SSR in every 6.0 kb of analyzed sequences. The trinucleotide repeats were predominant in all the unigenes with the exception of C. cajan, which showed prevalence of dinucleotide repeats over trinucleotide repeats. Dinucleotide repeats along with trinucleotides counted for more than 90% of the total microsatellites. Among dinucleotide and trinucleotide repeats, AG and AAG motifs, respectively, were the most frequent. Microsatellite positive chickpea unigenes were assigned Gene Ontology (GO) terms to identify the possible role of unigenes in various molecular and biological functions. These unigene based microsatellite markers will prove valuable for recording allelic variance across germplasm collections, gene tagging and searching for putative candidate genes.