ABSTRACT
BACKGROUND: One of the concerns of assembling de novo transcriptomes is determining the amount of read sequences required to ensure a comprehensive coverage of genes expressed in a particular sample. In this report, we describe the use of Illumina paired-end RNA-Seq (PE RNA-Seq) reads from Hevea brasiliensis (rubber tree) bark to devise a transcript mapping approach for the estimation of the read amount needed for deep transcriptome coverage. FINDINGS: We optimized the assembly of a Hevea bark transcriptome based on 16 Gb Illumina PE RNA-Seq reads using the Oases assembler across a range of k-mer sizes. We then assessed assembly quality based on transcript N50 length and transcript mapping statistics in relation to (a) known Hevea cDNAs with complete open reading frames, (b) a set of core eukaryotic genes and (c) Hevea genome scaffolds. This was followed by a systematic transcript mapping process where sub-assemblies from a series of incremental amounts of bark transcripts were aligned to transcripts from the entire bark transcriptome assembly. The exercise served to relate read amounts to the degree of transcript mapping level, the latter being an indicator of the coverage of gene transcripts expressed in the sample. As read amounts or datasize increased toward 16 Gb, the number of transcripts mapped to the entire bark assembly approached saturation. A colour matrix was subsequently generated to illustrate sequencing depth requirement in relation to the degree of coverage of total sample transcripts. CONCLUSIONS: We devised a procedure, the "transcript mapping saturation test", to estimate the amount of RNA-Seq reads needed for deep coverage of transcriptomes. For Hevea de novo assembly, we propose generating between 5-8 Gb reads, whereby around 90% transcript coverage could be achieved with optimized k-mers and transcript N50 length. The principle behind this methodology may also be applied to other non-model plants, or with reads from other second generation sequencing platforms.
Subject(s)
Gene Expression Profiling/methods , Gene Expression Regulation, Plant , Hevea/genetics , RNA, Messenger/genetics , RNA, Plant/genetics , Transcriptome , Databases, Genetic , Gene Library , Genes, Plant , Hevea/chemistry , High-Throughput Nucleotide Sequencing , Open Reading Frames , Plant Bark/metabolism , Plant Leaves/metabolism , Plant Proteins/genetics , RNA, Messenger/biosynthesis , RNA, Messenger/chemistry , RNA, Messenger/isolation & purification , RNA, Plant/biosynthesis , RNA, Plant/chemistry , RNA, Plant/isolation & purification , Reproducibility of ResultsABSTRACT
The cytosolic mevalonate (MVA) pathway in Hevea brasiliensis latex is the conventionally accepted pathway which provides isopentenyl diphosphate (IPP) for cis-polyisoprene (rubber) biosynthesis. However, the plastidic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway may be an alternative source of IPP since its more recent discovery in plants. Quantitative RT-PCR (qRT-PCR) expression profiles of genes from both pathways in latex showed that subcellular compartmentalization of IPP for cis-polyisoprene synthesis is related to the degree of plastidic carotenoid synthesis. From this, the occurrence of two schemes of IPP partitioning and utilization within one species is proposed whereby the supply of IPP for cis-polyisoprene from the MEP pathway is related to carotenoid production in latex. Subsequently, a set of latex unique gene transcripts was sequenced and assembled and they were then mapped to IPP-requiring pathways. Up to eight such pathways, including cis-polyisoprene biosynthesis, were identified. Our findings on pre- and post-IPP metabolic routes form an important aspect of a pathway knowledge-driven approach to enhancing cis-polyisoprene biosynthesis in transgenic rubber trees.