ribosome profiling是一种新出现的实验技术,它结合了高通量测序技术的应用,可以在全基因组水平上监测蛋白质的翻译水平,给研究细胞内蛋白质的翻译调控机制带来了新的可能。
目前ribo-seq建库测序的原理主要是通过对细胞使用翻译抑制剂,使正在翻译的核糖体固定在mrna序列或者起始位点上。裂解细胞后在细胞裂解液中加入rnase,消化不受核糖体保护的mrna,分离单一核糖体并提取纯化核糖体上未被消化的短片段mrna,进行建库测序和相应的数据分析。
康成生物丨数谱生物采用优化的建库方法,利用rrna remove kit去除非目标rna, 随后进行建库测序,该方法提高了基因的检出率和数据匹配率。同时,优化的核酸回收步骤也可以最大限度减少实验过程中rna的损失。
1 ribo-seq数据长度统计
高质量的ribo-seq文库的读长分布peaks主要集中在26-34nt左右范围(峰值最高点对应的rpf长度值随样本不同有差异),尤其是真核生物的细胞质核糖体。这里,我们根据测序的mrna的片段分布来评估实验富集的效率。
2 reads gene feature分布图
将比对到基因组上的reads在gene feature中的分布情况进行统计,定位区域分为5’utr,cds,3’utr。
3 principal component analysis 主成分分析
主成分分析(pca)是一种无监督特征学习的数据降维算法,通过样本的表达数据展示数据的分类情况。 利用 pca 分析可获得样本在实验组和对照组之间的直观分布情况,从而便于对离群样本进行检测和去除,并能找出相似度高的样本集合。我们采用所有样本中均数差别(anova)显著的基因(p value ≤ 0.05)做 pca 分析,绘制得到 pca 图(无重复样本使用全部 gene 做 pca 图)
pca 结果展示了样本间数据的分类情况,每种颜色代表各个分组,每个坐标轴代表一种主成分。
4 ribo-seq样品间相关性分析
样品间基因表达水平相关性是检验实验可靠性和样本选择是否合理的重要指标。相关系数越接近1,表明样品之间表达模式的相似程度越高。若样品中有生物学重复,通常生物重复相关系数要求较高。
样品相关系数heatmap图,蓝色深浅代表相关性高低。
5 差异te(translation efficiency翻译效率)基因聚类分析
ribo-c vs ribo-t 显著差异te基因聚类示例。差异表达基因聚类图。每行代表一个基因,每列代表一个样本。红色代表te上调显著差异基因,绿色代表te下调显著差异基因。
注:te的计算需要同时做原样本的rna seq。
显著差异te基因聚类示例
6 差异显著te基因go富集分析
差异表达基因的 go 富集分析结果(top 3)
go 分析(bp、cc、mf)结果中差异最显著的前 10 个 go 条目 bar 图。按照 p-value 从低到高排列,纵坐标表示 p-value(-log10 转化)。
7 差异显著te基因kegg富集分析
前 5 个差异表达显著基因的 pathway 结果
benjamini-hochberg false discovery rate (fdr) 校正 p_value