Mapping Accuracy of Short Reads from Massively Parallel Sequencing and the Implications for Quantitative Expression Profiling

Background

Massively parallel sequencing offers an enormous potential for expression profiling, in particular for interspecific comparisons. Currently, different platforms for massively parallel sequencing are available, which differ in read length and sequencing costs. The 454-technology offers the highest read length. The other sequencing technologies are more cost effective, on the expense of shorter reads. Reliable expression profiling by massively parallel sequencing depends crucially on the accuracy to which the reads could be mapped to the corresponding genes.

Methodology/Principal Findings

We performed an in silico analysis to evaluate whether incorrect mapping of the sequence reads results in a biased expression pattern. A comparison of six available mapping software tools indicated a considerable heterogeneity in mapping speed and accuracy. Independently of the software used to map the reads, we found that for compact genomes both short (35 bp, 50 bp) and long sequence reads (100 bp) result in an almost unbiased expression pattern. In contrast, for species with a larger genome containing more gene families and repetitive DNA, shorter reads (35–50 bp) produced a considerable bias in gene expression. In humans, about 10% of the genes had fewer than 50% of the sequence reads correctly mapped. Sequence polymorphism up to 9% had almost no effect on the mapping accuracy of 100 bp reads. For 35 bp reads up to 3% sequence divergence did not affect the mapping accuracy strongly. The effect of indels on the mapping efficiency strongly depends on the mapping software.

Conclusions/Significance

In complex genomes, expression profiling by massively parallel sequencing could introduce a considerable bias due to incorrectly mapped sequence reads if the read length is short. Nevertheless, this bias could be accounted for if the genomic sequence is known. Furthermore, sequence polymorphisms and indels also affect the mapping accuracy and may cause a biased gene expression measurement. The choice of the mapping software is highly critical and the reliability depends on the presence/absence of indels and the divergence between reads and the reference genome. Overall, we found SSAHA2 and CLC to produce the most reliable mapping results.     

兰州熊蜂蜂王头部产卵相关microRNA生物信息学分析

【目的】头部是熊蜂感觉和取食中心,也是生理行为的指挥中心。microRNA（miRNA）参与重要生理行为的调控。本研究初步探索了产卵和未产卵兰州熊蜂Bombus lantschouensis蜂王头部差异表达miRNA及其靶基因,以期探明蜂王头部miRNA在产卵过程中作用。【方法】以本土优良蜂种兰州熊蜂 B. lantschouensis 为材料,利用Solexa高通量测序技术对产卵和未产卵蜂王头部miRNA进行测序,运用生物信息学软件对miRNA及其靶基因进行预测;并应用qPCR技术验证极显著差异表达的miRNA。【结果】共获得兰州熊蜂产卵和未产卵蜂王头部miRNA的clean data数,分别为14 228 864和21 431 031 条reads;长度分析表明,在19~24 nt序列中22 nt序列数量最多,分别占产卵和未产卵蜂王序列的45.8%和45.1%。miRNA预测结果显示,共获得297个miRNA,其中270个为已知miRNA,有92个存在于蜜蜂中,其余178个分别存在于蜜蜂以外的其他物种中;另外27个为新的miRNA。在蜜蜂92个已知的miRNA中,在产卵蜂王头部中表达的共有91个,特异表达的有2个;在未产卵蜂王头部中表达的共有90个,特异表达的有1个。27个新的miRNA中,在产卵蜂王头部中表达的共有22个,特异表达的有2个;未产卵蜂王头部中表达的共有25个,特异表达的有5个。miRNA差异表达分析表明,在产卵与未产卵蜂王头部,共有8个已知miRNA表达量达到差异极显著水平（P<0.01）。qPCR结果表明,这8个表达差异极显著的miRNA均存在于熊蜂体内。在8个差异极显著的miRNA中,3个miRNA进行靶基因预测,发现它们共同调控9个靶基因,这些基因主要参与GTP结合和转录调控。【结论】本研究首次利用高通量测序技术鉴定了熊蜂头部组织的miRNA,并发现蜂王产卵与否受到miRNA差异表达的调控。结果为今后深入开展熊蜂miRNA功能验证及产卵调控机制研究奠定了基础。     

拟南芥AtLH基因过量表达引起叶向下性卷曲和晚开花

AtLH基因是BcpLH基因在拟南芥(Arapsis thaliana L.)中的同源基因,含有两个编码双链RNA结合蛋白的结构域.在大白菜叶球发育过程中,BcpLH基因与包叶的卷曲有关.为研究AtLH的基因对叶卷曲这一重要生物学现象的调控作用,构建了35S:AtLH基因的正义表达载体并转化拟南芥.与野生型比较页言,转基因植株的花和叶中AtLH的表达量有显著增加,成为AtLH基因过量的植株.这些植株的莲座叶向外或向下卷曲,呈现明显的偏上性生长;而且抽苔和开花时间延迟;在营养生长期其短缩茎的叶腑处着生数个侧茎,表现为顶端优势减弱;在生殖生长期二级花序减少使得主花序更加发达,表现为顶端优势增强,转基因植株对激素的敏感性改变,IAA剌激根生长的作用增强,ABA抑制根生长的作用减弱.由此可见,AtLH基因的过量表达可引起转基因植株的叶片向下卷曲.