如何预测转录因子靶基因!

做完转录组分析完之后,一般会找到一些差异表达的基因。差异表达基因中如果有转录因子的话,可以研究一下转录因子调控哪些基因的表达。要预测转录因子调控的靶基因,需要两个步骤:

1. 首先需要知道转录因子的结合序列特征

2. 基于转录因子的结合序列特征去基因上游的promoter区搜索,如果能搜索到结合特征的序列,那么该转录因子就有可能调控该靶基因。

今天我们就介绍采用JASPAR网站来完成转录因子调控靶基因的预测

转录因子结合序列特征的查找

1. 转录因子结合序列特征,这个可以采用JASPAR网站(http://jaspar.genereg.net/)中的"Profile Inference" 功能实现。输入转录因子的氨基酸序列,如下图所示:

2. 预测的结果可能会很多,一般选择E-value最小的结合Maxtrix作为转录因子的结合位点特征,如下图所示:

3. 查看结合序列的详细信息:可以下载结合位点的特征信息, 在本地进行转录因子调控基因的批量预测,也可以采用JASPAR 网站提供的scan功能,进行少量基因的调控预测(下面会详述)。

JASPAR预测转录因子靶基因

少量基因启动子区的结合位点预测可以采用网站提供是“Scan” 实现。具体步骤如下:

1. 首先搜索转录因子结合位点: 如上一步"Profile Inference"预测到的转录因子结合位点特征

2. 选择要搜索的转录因子结合位点特征

3. 输入启动子区的序列,点击“Scan” 运行

4. 预测结果如下图所示,会显示转录因子结合位点在promoter区域的位置信息和打分信息:

当然这个网站,只适合对转录因子进行少量的靶基因预测,如果需要对大量的靶基因进行预测,需要将转录因子结合位点的PFM文件下载到本地,采用一些本地的软件进行大规模的预测。

更多生物信息课程:

1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4分文章简单快速,学习链接:基因家族分析实操课程基因家族文献思路解读

2. 转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接:转录组(有参)结果解读转录组(无参)结果解读

3. 转录组数据深入挖掘技能-WGCNA,提升你的文章档次,学习链接:WGCNA-加权基因共表达网络分析

4. 转录组数据怎么挖掘?学习链接:转录组标准分析后的数据挖掘转录组文献解读

5. 微生物16S/ITS/18S分析原理及结果解读OTU网络图绘制cytoscape与网络图绘制课程

6. 生物信息入门到精通必修基础课,学习链接:linux系统使用perl入门到精通perl语言高级R语言画图

7. 医学相关数据挖掘课程,不用做实验也能发文章,学习链接:TCGA-差异基因分析GEO芯片数据挖掘GSEA富集分析课程TCGA临床数据生存分析TCGA-转录因子分析TCGA-ceRNA调控网络分析

8.其他课程链接:二代测序转录组数据自主分析NCBI数据上传二代测序数据解读

i

推荐阅读更多精彩内容