做完轉(zhuǎn)錄組分析完之后萤晴,一般會找到一些差異表達的基因吐句。差異表達基因中如果有轉(zhuǎn)錄因子的話,可以研究一下轉(zhuǎn)錄因子調(diào)控哪些基因的表達店读。要預測轉(zhuǎn)錄因子調(diào)控的靶基因嗦枢,需要兩個步驟:
1. 首先需要知道轉(zhuǎn)錄因子的結(jié)合序列特征
2. 基于轉(zhuǎn)錄因子的結(jié)合序列特征去基因上游的promoter區(qū)搜索,如果能搜索到結(jié)合特征的序列屯断,那么該轉(zhuǎn)錄因子就有可能調(diào)控該靶基因文虏。
今天我們就介紹采用JASPAR網(wǎng)站來完成轉(zhuǎn)錄因子調(diào)控靶基因的預測
轉(zhuǎn)錄因子結(jié)合序列特征的查找
1. 轉(zhuǎn)錄因子結(jié)合序列特征,這個可以采用JASPAR網(wǎng)站(http://jaspar.genereg.net/)中的"Profile Inference" 功能實現(xiàn)裹纳。輸入轉(zhuǎn)錄因子的氨基酸序列择葡,如下圖所示:
2. 預測的結(jié)果可能會很多紧武,一般選擇E-value最小的結(jié)合Maxtrix作為轉(zhuǎn)錄因子的結(jié)合位點特征剃氧,如下圖所示:
3. 查看結(jié)合序列的詳細信息:可以下載結(jié)合位點的特征信息, 在本地進行轉(zhuǎn)錄因子調(diào)控基因的批量預測阻星,也可以采用JASPAR 網(wǎng)站提供的scan功能朋鞍,進行少量基因的調(diào)控預測(下面會詳述)。
JASPAR預測轉(zhuǎn)錄因子靶基因
少量基因啟動子區(qū)的結(jié)合位點預測可以采用網(wǎng)站提供是“Scan” 實現(xiàn)妥箕。具體步驟如下:
1. 首先搜索轉(zhuǎn)錄因子結(jié)合位點: 如上一步"Profile Inference"預測到的轉(zhuǎn)錄因子結(jié)合位點特征
2. 選擇要搜索的轉(zhuǎn)錄因子結(jié)合位點特征
3. 輸入啟動子區(qū)的序列滥酥,點擊“Scan” 運行
4. 預測結(jié)果如下圖所示,會顯示轉(zhuǎn)錄因子結(jié)合位點在promoter區(qū)域的位置信息和打分信息:
當然這個網(wǎng)站畦幢,只適合對轉(zhuǎn)錄因子進行少量的靶基因預測坎吻,如果需要對大量的靶基因進行預測,需要將轉(zhuǎn)錄因子結(jié)合位點的PFM文件下載到本地宇葱,采用一些本地的軟件進行大規(guī)模的預測瘦真。
更多生物信息課程:
1. 文章越來越難發(fā)?是你沒發(fā)現(xiàn)新思路黍瞧,基因家族分析發(fā)2-4分文章簡單快速诸尽,學習鏈接:基因家族分析實操課程、基因家族文獻思路解讀
2. 轉(zhuǎn)錄組數(shù)據(jù)理解不深入印颤?圖表看不懂您机?點擊鏈接學習深入解讀數(shù)據(jù)結(jié)果文件,學習鏈接:轉(zhuǎn)錄組(有參)結(jié)果解讀年局;轉(zhuǎn)錄組(無參)結(jié)果解讀
3. 轉(zhuǎn)錄組數(shù)據(jù)深入挖掘技能-WGCNA际看,提升你的文章檔次,學習鏈接:WGCNA-加權(quán)基因共表達網(wǎng)絡(luò)分析
4. 轉(zhuǎn)錄組數(shù)據(jù)怎么挖掘矢否?學習鏈接:轉(zhuǎn)錄組標準分析后的數(shù)據(jù)挖掘仲闽、轉(zhuǎn)錄組文獻解讀
5.?微生物16S/ITS/18S分析原理及結(jié)果解讀、OTU網(wǎng)絡(luò)圖繪制兴喂、cytoscape與網(wǎng)絡(luò)圖繪制課程
6. 生物信息入門到精通必修基礎(chǔ)課蔼囊,學習鏈接:linux系統(tǒng)使用焚志、perl入門到精通、perl語言高級畏鼓、R語言畫圖
7. 醫(yī)學相關(guān)數(shù)據(jù)挖掘課程酱酬,不用做實驗也能發(fā)文章,學習鏈接:TCGA-差異基因分析云矫、GEO芯片數(shù)據(jù)挖掘膳沽、GSEA富集分析課程、TCGA臨床數(shù)據(jù)生存分析让禀、TCGA-轉(zhuǎn)錄因子分析挑社、TCGA-ceRNA調(diào)控網(wǎng)絡(luò)分析
8.其他課程鏈接:二代測序轉(zhuǎn)錄組數(shù)據(jù)自主分析、NCBI數(shù)據(jù)上傳巡揍、二代測序數(shù)據(jù)解讀痛阻。
i