基因组学

chipseeker包的一些参数上的理解

2021-02-14  本文已影响0人  晓颖_9b6f

关于tssRegion参数:

理解:
1、
TSS是: 转录时,mRNA链第一个核苷酸相对应DNA链上的碱基
它的上游是promoter
它的下游就是相应的基因

2、做注释的时候,需要找到离TSS最近的基因,离结合位点最近的基因更有可能被调控
所以annotatePeak就是可以查看peak上下游某个范围内(比如说-3k到3k的距离)都有什么基因。
只要peaks和TSS上下游的定义的区域,比如说-3k到3k的距离,有交集,那么这个peaks就可以被定义为promoter。

3、chipseeker的使用报告里面是写着,默认的tssRegion是+-3000的范围。


image.png

4、但是我在网上找到的很多流程,对于tssRegion的范围,
tssRegion = c(-2500, 2500)
有的是1000
有的是3000
也是有2000的?
不过也没说,为什么要取这个范围值???

5、我的上一次流程里面,x用的是1000。
是想看到,TSS附近1000以内的基因。

b1b012a32a905abe9ade662ef09025e.png e690edde5530a3dba25bb5d448ece98.png

6、chipseeker的使用报告:seq2gene函数,用于将基因组区域与基因进行多对多映射。它考虑了宿主基因(外显子/内含子),启动子区域和基因间区域的侧翼基因,这些基因可能受顺式调控。此功能旨在将编码和非编码基因组区域链接到编码基因,并促进功能分析。
flankDistance = 3000,是想看左右3000的距离上有什么基因?

gene <- seq2gene(peak, tssRegion = c(-1000, 1000), flankDistance = 3000, TxDb=txdb)

7、总的来说,网上找了这么些资料,感觉还是不太理解这些范围的调节选择?
①、annotatePeak的时候,tssRegion究竟取多少才好?一般是用默认的3000。
如果调大一点,调到5000,找到的基因多一点吗?
如果调小一点,调到1000,是基因少一点,准确性高一点?

②、后面在做KEGG的时候,tssRegion ,选择,和前面annotatePeak,不需要一致。
但是选择多少比较合适?

上一篇下一篇

猜你喜欢

热点阅读