Rigidoporus microporus is the fungus accountable for the white root rot disease that is detrimental to the rubber tree, Hevea brasiliensis. The pathogenicity mechanism of R. microporus and the identity of the fungal proteins and metabolites involved during the infection process remain unclear. In this study, the protein and metabolite profiles of two R. microporus isolates, Segamat (SEG) and Ayer Molek (AM), were investigated during an in vitro interaction with H. brasiliensis. The isolates were used to inoculate H. brasiliensis clone RRIM 2025, and mycelia adhering to the roots of the plant were collected for analysis. Transmission electron microscope (TEM) images acquired confirms the hyphae attachment and colonization of the mycelia on the root of the H. brasiliensis clones after 4 days of inoculation. The protein samples were subjected to 2-DE analysis and analyzed using MALDI-ToF MS/MS, while the metabolites were extracted using methanol and analyzed using LC/MS-QTOF. Based on the differential analyses, upregulation of proteins that are essential for fungal evolution such as malate dehydrogenase, fructose 1,6-biphosphate aldolase, and glyceraldehyde-3-phosphate dehydrogenase hints an indirect role in fungal pathogenicity, while metabolomic analysis suggests an increase in acidic compounds which may lead to increased cell wall degrading enzyme activity. Bioinformatics analyses revealed that the carbohydrate and amino acid metabolisms were prominently affected in response to the fungal pathogenicity. In addition to that, other pathways that were significantly affected include "Protein Ubiquitination Pathway," Unfolded Protein Response," "HIFα Signaling," and "Sirtuin Signaling Pathway." The identification of responsive proteins and metabolites from this study promotes a better understanding of mechanisms underlying R. microporus pathogenesis and provides a list of potential biological markers for early recognition of the white root rot disease.
Natural rubber or latex from the Hevea brasiliensis is an important commodity in various economic sectors in today's modern society. Proteins have been detected in latex since the early twentieth century, and they are known to regulate various biological pathways within the H. brasiliensis trees such as the natural rubber biosynthesis, defence against pathogens, wound healing, and stress tolerance. However, the exact mechanisms of the pathways are still not clear. Proteomic analyses on latex have found various proteins and revealed how they fit into the mechanisms of the biological pathways. In the past three decades, there has been rapid latex protein identification due to the improvement of latex protein extraction methods, as well as the emergence of two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). In this manuscript, we reviewed the methods of latex protein extraction that keeps on improving over the past three decades as well as the results of numerous latex protein identification and quantitation.
Hevea brasiliensis is exploited for its latex production, and it is the only viable source of natural rubber worldwide. The demand for natural rubber remains high due its high-quality properties, which synthetic rubber cannot compete with. In this paper, we present transcriptomic data and analysis of three H. brasiliensis clones using tissue from latex and bark tissues collected from 10-year-old plant. The combined, assembled transcripts were mapped onto an H. brasiliensis draft genome. Gene ontology analysis showed that the most abundant transcripts related to molecular functions, followed by biological processes and cellular components. Simple sequence repeats (SSR) and single nucleotide polymorphisms (SNP) were also identified, and these can be useful for selection of parental and new clones in a breeding program. Data generated by RNA sequencing were deposited in the NCBI public repository under accession number PRJNA629890.
In this study, RNA sequencing of several Hevea brasiliensis clones grown in Malaysia with different annual rubber production yields and disease resistance was performed on the Illumina platform. A total of 29,862,548 reads were generated, resulting in 101,269 assembled transcripts that were used as the reference transcripts. A similarity search against the non-redundant (nr) protein databases presented 83,771 (83%) positive BLASTx hits. The transcriptome was annotated using gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Pfam database. A search for putative molecular markers was performed to identify single-nucleotide polymorphisms (SNPs). Overall, 3,210,629 SNPs were detected and a total of 1314 SNPs associated with the genes involved in MVA and MEP pathways were identified. A total of 176 SNP primer pairs were designed from sequences that were related to the MVA and MEP pathways. The transcriptome of RRIM 3001 and RRIM 712 were subjected to pairwise comparison and the results revealed that there were 1262 significantly differentially expressed genes unique to RRIM 3001, 1499 significantly differentially expressed genes unique to RRIM 712 and several genes related to the MVA and MEP pathways such as AACT, HMGS, PMK, MVD, DXS and HDS were included. The results will facilitate the characterization of H. brasiliensis transcriptomes and the development of a new set of molecular markers in the form of SNPs from transcriptome assembly for the genotype identification of various rubber varieties with superior traits in Malaysia.
Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.
The natural rubber latex extracted from the bark of Hevea brasiliensis plays various important roles in modern society. Post-translational modifications (PTMs) of the latex proteins are important for the stability and functionality of the proteins. In this study, latex proteins were acquired from the C-serum, lutoids, and rubber particle layers of latex without using prior enrichment steps; they were fragmented using collision-induced dissociation (CID), higher-energy collisional dissociation (HCD), and electron-transfer dissociation (ETD) activation methods. PEAKS 7 were used to search for unspecified PTMs, followed by analysis through PTM prediction tools to crosscheck both results. There were 73 peptides in 47 proteins from H. brasiliensis protein sequences derived from UniProtKB were identified and predicted to be post-translationally modified. The peptides with PTMs identified include phosphorylation, lysine acetylation, N-terminal acetylation, hydroxylation, and ubiquitination. Most of the PTMs discovered have yet to be reported in UniProt, which would provide great assistance in the research of the functional properties of H. brasiliensis latex proteins, as well as being useful biomarkers. The data are available via the MassIVE repository with identifier MSV000082419.
The stabilization mechanism of natural rubber (NR) latex from Hevea brasiliensis was studied to investigate the components involved in base-catalyzed ester hydrolysis, namely, hydrolyzable lipids, ammonia, and the products responsible for the desired phenomenon observed in ammonia-preserved NR latex. Latex stability is generally thought to come from a rubber particle (RP) dispersion in the serum, which is encouraged by negatively charged species distributed on the RP surface. The mechanical stability time (MST) and zeta potential were measured to monitor field latices preserved in high (FNR-HA) and low ammonia (FNR-LA) contents as well as that with the ester-containing components removed (saponified NR) at different storage times. Amounts of carboxylates of free fatty acids (FFAs), which were released by the transformation and also hypothesized to be responsible for the like-charge repulsion of RPs, were measured as the higher fatty acid (HFA) number and corroborated by confocal laser scanning microscopy (CLSM) both qualitatively and quantitatively. The lipids and their FFA products interact differently with Nile red, which is a lipid-selective and polarity-sensitive fluorophore, and consequently re-emit characteristically. The results were confirmed by conventional ester content determination utilizing different solvent extraction systems to reveal that the lipids hydrolyzed to provide negatively charged fatty acid species were mainly the polar lipids (glycolipids and phospholipids) at the RP membrane but not those directly linked to the rubber molecule and, to a certain extent, those suspended in the serum. From new findings disclosed herein together with those already reported, a new model for the Hevea rubber particle in the latex form is proposed.
Hevea brasiliensis is one of the popular industrial crops in Malaysia better known as rubber tree belongs to the family Euphorbiaceace. From more than 12 species of Hevea, only Hevea brasiliensis is economically exploited because the milky latex extracted from the tree is the primary source of natural rubber. As in other crops, various plant physiological conditions and pathogenic diseases influence rubber production. Brown bast is one of the most serious threats to natural rubber production. In general, high-yielding clones of rubber tree are often considered to be more susceptible to this physiological disorder also commonly termed Tapping Panel Dryness (TPD). It is estimated that brown bast leads to approximately 15-20% decrease in yield. There is no known cure for brown bast yet. However, many plantation practices manage brown bast in rubber by giving tapping rest and changing tapping panel. Hence, this review condenses the causal of brown bast, symptoms of diseases and also control of brown bast affected Hevea tree.
Natural rubber (polyisoprene) from the rubber tree Hevea brasiliensis is synthesized by specialized cells called laticifers. It is not clear how rubber particles arise, although one hypothesis is that they derive from the endoplasmic reticulum (ER) membrane. Here we cloned the genes encoding four key proteins found in association with rubber particles and studied their intracellular localization by transient expression in Nicotiana benthamiana leaves. We show that, while the cis-prenyltransferase (CPT), responsible for the synthesis of long polyisoprene chains, is a soluble, cytosolic protein, other rubber particle proteins such as rubber elongation factor (REF), small rubber particle protein (SRPP) and Hevea rubber transferase 1-REF bridging protein (HRBP) are associated with the endoplasmic reticulum (ER). We also show that SRPP can recruit CPT to the ER and that interaction of CPT with HRBP leads to both proteins relocating to the plasma membrane. We discuss these results in the context of the biogenesis of rubber particles.
Chemical interesterification of rubber seed oil has been investigated for four different designed orifice devices in a pilot scale hydrodynamic cavitation (HC) system. Upstream pressure within 1-3.5bar induced cavities to intensify the process. An optimal orifice plate geometry was considered as plate with 1mm dia hole having 21 holes at 3bar inlet pressure. The optimisation results of interesterification were revealed by response surface methodology; methyl acetate to oil molar ratio of 14:1, catalyst amount of 0.75wt.% and reaction time of 20min at 50°C. HC is compared to mechanical stirring (MS) at optimised values. The reaction rate constant and the frequency factor of HC were 3.4-fold shorter and 3.2-fold higher than MS. The interesterified product was characterised by following EN 14214 and ASTM D 6751 international standards.
Natural rubber has unique physical properties that cannot be replaced by products from other latex-producing plants or petrochemically produced synthetic rubbers. Rubber from Hevea brasiliensis is the main commercial source for this natural rubber that has a cis-polyisoprene configuration. For sustainable production of enough rubber to meet demand elucidation of the molecular mechanisms involved in the production of latex is vital. To this end, we firstly constructed rubber full-length cDNA libraries of RRIM 600 cultivar and sequenced around 20,000 clones by the Sanger method and over 15,000 contigs by Illumina sequencer. With these data, we updated around 5,500 gene structures and newly annotated around 9,500 transcription start sites. Second, to elucidate the rubber biosynthetic pathways and their transcriptional regulation, we carried out tissue- and cultivar-specific RNA-Seq analysis. By using our recently published genome sequence, we confirmed the expression patterns of the rubber biosynthetic genes. Our data suggest that the cytoplasmic mevalonate (MVA) pathway is the main route for isoprenoid biosynthesis in latex production. In addition to the well-studied polymerization factors, we suggest that rubber elongation factor 8 (REF8) is a candidate factor in cis-polyisoprene biosynthesis. We have also identified 39 transcription factors that may be key regulators in latex production. Expression profile analysis using two additional cultivars, RRIM 901 and PB 350, via an RNA-Seq approach revealed possible expression differences between a high latex-yielding cultivar and a disease-resistant cultivar.
The natural rubber latex extracted from the bark of Hevea brasiliensis plays various important roles in today's modern society. Following ultracentrifugation, the latex can be separated into 3 layers: C-serum, lutoids, and rubber particles. Previous studies have shown that a large number of proteins are present in these 3 layers. However, a complete proteome for this important plant is still unavailable. Protein sequences have been recently translated from the completed draft genome database of H. brasiliensis, leading to the creation of annotated protein databases of the following H. brasiliensis biosynthetic pathways: photosynthesis, latex allergens, rubberwood formation, latex biosynthesis, and disease resistance. This research was conducted to identify the proteins contained within the latex by way of de novo sequencing from mass spectral data obtained from the 3 layers of the latex. Peptides from these proteins were fragmented using collision-induced dissociation, higher-energy collisional dissociation, and electron-transfer dissociation activation methods. A large percentage of proteins from the biosynthetic pathways (63% to 100%) were successfully identified. In addition, a total of 1839 unique proteins were identified from the whole translated draft genome database (AnnHBM).
Ammonium-enriched skim latex serum - used for absorption of contaminating ammonia gas - when composted with other rubber tree wastes, is promising as a good compost. The objective of this research was to utilize ammonium-enriched skim latex serum (S) as a raw composting ingredient after being combined with para sawdust (W1) and para rubber leaves (W2). Several ratios of S, W1 and W2 were experimented in a 15L composting vessel to determine the most effective compost. The best ratio was found to be 3:1:3 by weight at 12-day retention. The modified 30 L composting reactor employed with the derived optimum mixing conditions yielded N, P and K of 2.40, 1.51 and 14.84 %w/w. The growth of Brassica alboglabra after application of this compost combined with a chemical fertilizer generated the highest fresh weight (4.48 g/plant). Thus, compost from these wastes could be used as a fertilizer and logically should contribute to cost saving of waste disposal.
Corynespora cassiicola is a common plant pathogen that causes leaf spot disease in a broad range of crop, and it heavily affect rubber trees in Malaysia (Hsueh, 2011; Nghia et al., 2008). The isolation of UM 591 from a patient's contact lens indicates the pathogenic potential of this dematiaceous fungus in human. However, the underlying factors that contribute to the opportunistic cross-infection have not been fully studied. We employed genome sequencing and gene homology annotations in attempt to identify these factors in UM 591 using data obtained from publicly available bioinformatics databases. The assembly size of UM 591 genome is 41.8 Mbp, and a total of 13,531 (≥99 bp) genes have been predicted. UM 591 is enriched with genes that encode for glycoside hydrolases, carbohydrate esterases, auxiliary activity enzymes and cell wall degrading enzymes. Virulent genes comprising of CAZymes, peptidases, and hypervirulence-associated cutinases were found to be present in the fungal genome. Comparative analysis result shows that UM 591 possesses higher number of carbohydrate esterases family 10 (CE10) CAZymes compared to other species of fungi in this study, and these enzymes hydrolyses wide range of carbohydrate and non-carbohydrate substrates. Putative melanin, siderophore, ent-kaurene, and lycopene biosynthesis gene clusters are predicted, and these gene clusters denote that UM 591 are capable of protecting itself from the UV and chemical stresses, allowing it to adapt to different environment. Putative sterigmatocystin, HC-toxin, cercosporin, and gliotoxin biosynthesis gene cluster are predicted. This finding have highlighted the necrotrophic and invasive nature of UM 591.
The natural rubber of Para rubber tree, Hevea brasiliensis, is the main crop involved in industrial rubber production due to its superior quality. The Hevea bark is commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. The laticifer is well defined in the aspect of morphology; however, only some genes associated with its development have been reported. We successfully induced secondary laticifer in the jasmonic acid (JA)-treated and linolenic acid (LA)-treated Hevea bark but secondary laticifer is not observed in the ethephon (ET)-treated and untreated Hevea bark. In this study, we analysed 27,195 gene models using NimbleGen microarrays based on the Hevea draft genome. 491 filtered differentially expressed (FDE) transcripts that are common to both JA- and LA-treated bark samples but not ET-treated bark samples were identified. In the Eukaryotic Orthologous Group (KOG) analysis, 491 FDE transcripts belong to different functional categories that reflect the diverse processes and pathways involved in laticifer differentiation. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) and KOG analysis, the profile of the FDE transcripts suggest that JA- and LA-treated bark samples have a sufficient molecular basis for secondary laticifer differentiation, especially regarding secondary metabolites metabolism. FDE genes in this category are from the cytochrome (CYP) P450 family, ATP-binding cassette (ABC) transporter family, short-chain dehydrogenase/reductase (SDR) family, or cinnamyl alcohol dehydrogenase (CAD) family. The data includes many genes involved in cell division, cell wall synthesis, and cell differentiation. The most abundant transcript in FDE list was SDR65C, reflecting its importance in laticifer differentiation. Using the Basic Local Alignment Search Tool (BLAST) as part of annotation and functional prediction, several characterised as well as uncharacterized transcription factors and genes were found in the dataset. Hence, the further characterization of these genes is necessary to unveil their role in laticifer differentiation. This study provides a platform for the further characterization and identification of the key genes involved in secondary laticifer differentiation.
Hevea brasiliensis Muell. Arg, a member of the family Euphorbiaceae, is the sole natural resource exploited for commercial production of high-quality natural rubber. The properties of natural rubber latex are almost irreplaceable by synthetic counterparts for many industrial applications. A paucity of knowledge on the molecular mechanisms of rubber biosynthesis in high yield traits still persists. Here we report the comprehensive genome-wide analysis of the widely planted H. brasiliensis clone, RRIM 600. The genome was assembled based on ~155-fold combined coverage with Illumina and PacBio sequence data and has a total length of 1.55 Gb with 72.5% comprising repetitive DNA sequences. A total of 84,440 high-confidence protein-coding genes were predicted. Comparative genomic analysis revealed strong synteny between H. brasiliensis and other Euphorbiaceae genomes. Our data suggest that H. brasiliensis's capacity to produce high levels of latex can be attributed to the expansion of rubber biosynthesis-related genes in its genome and the high expression of these genes in latex. Using cap analysis gene expression data, we illustrate the tissue-specific transcription profiles of rubber biosynthesis-related genes, revealing alternative means of transcriptional regulation. Our study adds to the understanding of H. brasiliensis biology and provides valuable genomic resources for future agronomic-related improvement of the rubber tree.
Pretreatment of the high free fatty acid rubber seed oil (RSO) via esterification reaction has been investigated by using a pilot scale hydrodynamic cavitation (HC) reactor. Four newly designed orifice plate geometries are studied. Cavities are induced by assisted double diaphragm pump in the range of 1-3.5 bar inlet pressure. An optimised plate with 21 holes of 1mm diameter and inlet pressure of 3 bar resulted in RSO acid value reduction from 72.36 to 2.64 mg KOH/g within 30 min of reaction time. Reaction parameters have been optimised by using response surface methodology and found as methanol to oil ratio of 6:1, catalyst concentration of 8 wt%, reaction time of 30 min and reaction temperature of 55°C. The reaction time and esterified efficiency of HC was three fold shorter and four fold higher than mechanical stirring. This makes the HC process more environmental friendly.
A study has been carried out to characterize hydrocarbons emitted from the burning of three tropical wood species. The woods were burned to ember and smoke aerosols emitted were sampled using high volume sampler fitted with a pre-cleaned glass fibre filters. Hydrocarbons were extracted using ultrasonic agitation with dichloromethane-methanol (3:1 v/v) as solvent and the extracts obtained were then fractionated on silica-alumina column. Detection and quantification of aliphatic and polycyclic aromatic hydrocarbons (PAHs) compounds were carried out using GC-MS. The results indicated that the major aliphatic hydrocarbons characterized were straight chain n-alkanes in the range of C12-C35 with Cmax in the range of C27-C33. Rhizophora apiculata and Hevea brasiliensis wood smoke exhibited a weak odd to even carbon number predominance with carbon preference index (CPI) values greater than one whereas Melaleuca cajuputi wood smoke aerosols did not exhibit similar pattern with CPI obtained close to one. The results obtained also indicated that burning of these wood resulted in formation of PAHs compounds in their smoke aerosols with predominance of three to four rings PAHs over the two, five and lesser of six rings PAHs. PAH diagnostic ratios calculated except for Flan/(Flan+Py) and Indeno/(Indeno+BgP) were consistent with the ratios suggested for wood combustion source as reported in literatures. In the case of the latter, two diagnostic ratios, the values were generally lower than the range normally reported for wood combustion.
Rubber tree (Hevea brasiliensis) is an important crop in tropical regions of China. In October 2013, a new stem rot disease was found on cv. Yunyan77-4 at a rubber tree plantation in Hekou, Yunnan Province. There were about 100 plants, and diseased rubber trees accounted for 30% or less. Initially, brown-punctuate secretion appeared on the stem, which was 5 to 6 cm above the ground. Eventually, the secretion became black and no latex produced from the rubber tree bark. After removing the secretion, the diseased bark was brown putrescence, but the circumambient bark was normal. Upon peeling the surface bark, the inner bark and xylem had brown rot and was musty. The junction between health and disease was undulate. On the two most serious plants, parts of leaves on the crown were yellow, and the root near the diseased stem was dry and puce. The pathogen was isolated and designated HbFO01; the pathogenicity was established by following Koch's postulates. The pathogen was cultivated on a potato dextrose agar (PDA) plate at 28°C for 4 days. Ten plants of rubber tree cv. Yunyan77-4 were selected from a disease-free plantation in Haikou, Hainan Province, and the stem diameter was about 7 cm. The bark of five plants was peeled, and one mycelium disk with a diameter of 1 cm was inserted into the cut and covered again with the bark. The other five plants were treated with agar disks as controls. The inoculation site was kept moist for 2 days, and then the mycelium and agar disk were removed. On eighth day, symptoms similar to the original stem lesions were observed on stems of inoculated plants, while only scars formed on stems of control plants. The pathogen was re-isolated from the lesions of inoculated plants. On PDA plates, the pathogen colony was circular and white with tidy edges and rich aerial hyphae. Microscopic examination showed microconidia and chlamydospores were produced abundantly on PDA medium. The falciform macroconidia were only produced on lesions and were slightly curved, with a curved apical cell and foot shaped to pointed basal cell, usually 3-septate, 16.2 to 24.2 × 3.2 to 4.0 μm. Microconidia were produced in false heads, oval, 0-septate, 6.2 to 8.2 × 3.3 to 3.8 μm, and the phialide was cylindrical. Chlamydospores were oval, 6.4 to 7.2 × 3.1 to 3.8 μm, alone produced in hypha. Morphological characteristics of the specimen were similar to the descriptions for Fusarium oxysporum (2). Genomic DNA of this isolate was extracted with a CTAB protocol (4) from mycelium and used as a template for amplification of the internal transcribed spacer (ITS) region of rDNA with primer pair ITS1/ITS4 (1). The full length of this sequence is 503 nt (GenBank Accession No. KJ009335), which exactly matched several sequences (e.g., JF807394.1, JX897002.1, and HQ451888.1) of F. oxysporum. Williams and Liu had listed F. oxysporum as the economically important pathogen of Hevea in Asia (3), while this is, to our knowledge, the first report of stem rot caused by F. oxysporum on rubber tree in China. References: (1) D. E. L. Cooke et al. Fungal Genet. Biol. 30:17, 2000. (2) J. F. Leslie and B. A. Summerell. The Fusarium Laboratory Manual, 2006. (3) T. H. Williams and P. S. W. Liu. A host list of plant diseases in Sabah, Malaysia, 1976. (4) J. R. Xu et al. Genetics 143:175, 1996.
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.