RNA sequencing read depth requirement for optimal transcriptome coverage in Hevea brasiliensis

Area of Science:

• Genomics
• Bioinformatics
• Molecular Biology

Background:

• De novo transcriptome assembly requires sufficient read sequences for comprehensive gene coverage.
• Estimating the necessary read amount is a key challenge in transcriptome sequencing projects.
• Illumina paired-end RNA-Seq (PE RNA-Seq) data from Hevea brasiliensis bark was utilized.

Purpose of the Study:

• To devise a method for estimating the optimal number of RNA-Seq reads for deep transcriptome coverage.
• To establish a transcript mapping approach for assessing sequencing depth requirements.
• To optimize de novo transcriptome assembly for Hevea brasiliensis.

Main Methods:

• Optimized transcriptome assembly using the Oases assembler with varying k-mer sizes on 16 Gb of Illumina PE RNA-Seq reads.
• Assessed assembly quality using transcript N50 length and mapping statistics against known cDNAs, core eukaryotic genes, and genome scaffolds.
• Performed systematic transcript mapping by aligning sub-assemblies from incremental read amounts to the full assembly to determine saturation levels.

Main Results:

• Transcriptome assembly was optimized using 16 Gb of reads, with quality assessed via N50 and mapping statistics.
• Transcript mapping levels approached saturation as read amounts increased towards 16 Gb.
• A color matrix was generated to visualize sequencing depth requirements versus transcript coverage.

Conclusions:

• A "transcript mapping saturation test" was developed to estimate RNA-Seq read quantities for deep transcriptome coverage.
• For Hevea de novo assembly, 5-8 Gb of reads are proposed to achieve approximately 90% transcript coverage.
• The methodology is applicable to other non-model plants and sequencing platforms.