The reads were pooled together and assembled into reference sequences, yielding 72,732 contigs representing 272,053,772 total assembled bases for the P. ginseng transcriptome.
Last, through simulations with varying read lengths, we show that the optimal length for RNA-Seq is around 20 25 bases for the mouse and maize transcriptomes when our method is used to handle multireads.
For each transcriptome, approximately 99% of the genes had at least a partial hit in the tef genome while 96-99% of the bases in the transcriptomes could be aligned to the genome.
Our results indicate that a read length between 20 and 25 bases gives the highest accuracy for the transcriptomes and throughputs considered.
The percentages of Q20 bases for the clean reads in all eight transcriptomes were all above 96% (Additional file 3).
The percentage of Q20 bases for the clean reads in the four sub-transcriptomes were all ∼93% (Supplementary Table S1).
Actually, the 204 genes showed appreciable changes in expression patterns between dogs and wolves than others for two different measurements: Absolute expression change and fold change (P = 0.022 and P = 0.005, respectively) (fig. 1 B), based on the transcriptome data for the frontal cortex (Albert et al. 2012).
Based on the transcriptome results obtained for the well-studied stringent response of E. coli[ 7, 53] and on earlier reports [ 54], a major event of the stringent response consists in the repression of the translation apparatus, including ribosomal proteins.
These studies will benefit from high quality assemblies with few singletons (less than 26% of bases for each assembled transcriptome are present in unassembled singleton ESTs) and effective transcript discovery (over 6,500 of our putative orthologs cover at least 50% of the corresponding model silkworm gene).
Based on the transcriptome changes, we hereby propose a putative pathway for carbon metabolism in conidia exposed to neutrophils (Fig. 8).
