Cytoscape構(gòu)建網(wǎng)絡圖實操

基礎理論

網(wǎng)絡圖

網(wǎng)絡圖非常常見,不僅被用在生物信息學分析积仗,在生活中也很常見疆拘,如民航航線圖、食物鏈寂曹、基因調(diào)控網(wǎng)絡都是網(wǎng)絡圖經(jīng)典的例子哎迄。


南航航線圖

食物鏈

基因調(diào)控網(wǎng)絡

網(wǎng)絡圖應用價值

  • 直觀形象的呈現(xiàn)數(shù)據(jù)關(guān)系
  • 高效的挖掘和解析數(shù)據(jù)(WGCNA后模塊內(nèi)大量基因的篩選)

網(wǎng)絡圖的構(gòu)成基礎

  • 節(jié)點(node):表示元素
    節(jié)點屬性:大小、形狀稀颁、顏色芬失、標簽楣黍、邊框可以表示表達量匾灶、差異倍數(shù)、RNA類型租漂、pathway分類阶女、基因名稱等。
  • 線(edge):表示元素之間的關(guān)系
    線屬性:粗細哩治、類型秃踩、顏色、方向可以表示:相關(guān)性強弱业筏、相關(guān)性正負憔杨、靶向關(guān)系、相關(guān)性顯著性等蒜胖。

重要概念

  • 連通性(Degree):一個節(jié)點擁有的線的數(shù)量消别。
  • 核心元素(Hub gene):位于調(diào)控網(wǎng)絡中心,也就是連通性較高的元素台谢。

數(shù)據(jù)準備

建立關(guān)聯(lián)

  • 基于已有成果獲得寻狂。(string數(shù)據(jù)庫)
  • 基于表達量/豐度的相關(guān)性。(WGCNA)
  • 基于序列的堿基互補關(guān)系朋沮。(miRNA-mRNA)
  • 基于功能分類關(guān)系蛇券。(富集分析)

數(shù)據(jù)要求

  • 格式:
    tsv文件(Tab Separated Values)
    csv文件(Comma )
    xls,xlsx文件

  • 內(nèi)容:
    edge文件:必須文件,包含節(jié)點關(guān)系關(guān)鍵數(shù)據(jù)纠亚。如下圖所示塘慕,第一列為起始點的基因名,第二列為終止點的基因名菜枷,第三列為TOM值(這個數(shù)據(jù)是WGCNA的數(shù)據(jù)結(jié)果苍糠,因此我用了TOM值,也可以是其他的數(shù)據(jù)啤誊,如兩點間相關(guān)性岳瞭。)WGCNA分析有點復雜,有朝一日我整理好后再發(fā)給大家參考蚊锹。


    edge文件

node文件:非必須瞳筏,按需準備。如下圖所示牡昆,第一列是點的基因名姚炕,注意對應edge文件中的點的信息,我下圖展示的是我給每個基因Ensembl ID添加的Gene Symbol信息和注釋信息丢烘,也可以是其他的數(shù)據(jù)柱宦,如基因的表達量等。


node文件

軟件基本操作

數(shù)據(jù)導入

  • 導入edge數(shù)據(jù)
    點擊下圖所示的按鍵播瞳,導入edge數(shù)據(jù)掸刊。


    導入數(shù)據(jù)

    導入后需要選擇每一列的數(shù)據(jù)類型,主要設置邊的起始位置和結(jié)束位置赢乓。


    數(shù)據(jù)類型選擇
  • 導入node數(shù)據(jù)
    點擊對應位置忧侧,導入node數(shù)據(jù),要確保數(shù)據(jù)和edge數(shù)據(jù)的命名一致牌芋。


    導入數(shù)據(jù)

生成網(wǎng)絡圖

導入edge后蚓炬,在圖片編輯框中就可以看到生成的最基礎的網(wǎng)絡圖,我們下一步需要做的就是對它進行美化躺屁,也可以進一步通過Cytoscape進行數(shù)據(jù)挖掘肯夏。


最初版網(wǎng)絡圖

圖片美化

對于最初版的網(wǎng)絡圖,我相信沒有人會覺得好看犀暑,也絕對滿足不了雜志的要求驯击,那么,我們就動手美化它母怜。
其實Cytoscape的可視化界面非常人性化余耽,按照我下面的批注,大家都點一點苹熏,應該很快就能掌握碟贾。


參數(shù)說明

需要特別解釋一下“參數(shù)賦值”和“點/線排列形式”的用法币喧。

  • 參數(shù)賦值:比如,想要用圈的大小表示基因的連通性大小袱耽、圈的顏色表示基因的表達量高低杀餐、線的粗細表示權(quán)重值的大小等,都可以通過給各個區(qū)域賦值達成朱巨。在Cytoscape中史翘,點和線的顏色、粗細都可以賦值冀续,點擊相應的參數(shù)琼讽,賦予你想要展示的數(shù)據(jù)即可。
  • 點/線排列形式:在Cytoscape中洪唐,所有的點都是可以自己手動拖動的钻蹬,當數(shù)據(jù)量很大時,非常辛苦凭需。因此问欠,軟件有自動排布功能,點擊“l(fā)ayout”就可以看到粒蜈,自動排布形式主要有:矩形排布顺献、圓形排布、層級排布枯怖,每一種都有適用范圍注整,如果展示基因相互關(guān)系建議用圓形排布,點一下試試就知道嫁怀。

Tips:選擇數(shù)據(jù)點后(手動選擇結(jié)合下面要說到的篩選選擇)设捐,再點擊排布借浊,就可以把雜亂無章的圖塘淑,變得整齊且能說明問題,如下圖所示蚂斤。

網(wǎng)上不知名數(shù)據(jù)

數(shù)據(jù)挖掘

Cytoscape的數(shù)據(jù)挖掘主要是用到“篩選”工具存捺。網(wǎng)絡圖常見的篩選方式是利用連通性篩選關(guān)鍵基因,Cytoscape可以直接計算點的連通性(K)曙蒸,點擊菜單欄Tools--NetworkAnalyzer--Network analysis--Analyze Network捌治,然后選擇自己數(shù)據(jù)對應的類型(有/無方向),點擊確認纽窟,連通性數(shù)據(jù)就會出現(xiàn)在“點信息頁”肖油,列名為“degree”。


選擇數(shù)據(jù)類型

點擊篩選模塊臂港,新建篩選條件森枪,在下拉菜單中選擇“degree”视搏,即可用連通性作為篩選條件,選擇滿足連通性要求的點县袱。如果圖片中點和線特別多浑娜,可以用該方法選擇連通性較高的點新建畫布,只對這些連通性高的基因做圖式散。當然筋遭,篩選條件不僅僅是K值,也可以用WGCNA分析里的TOM值暴拄,或者基因的表達量等漓滔,只要是數(shù)值型變量就行,大家嘗試一下就明白了乖篷。


篩選

結(jié)果展示
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末次和,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子那伐,更是在濱河造成了極大的恐慌踏施,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,651評論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件罕邀,死亡現(xiàn)場離奇詭異畅形,居然都是意外死亡,警方通過查閱死者的電腦和手機诉探,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,468評論 3 392
  • 文/潘曉璐 我一進店門日熬,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人肾胯,你說我怎么就攤上這事竖席。” “怎么了敬肚?”我有些...
    開封第一講書人閱讀 162,931評論 0 353
  • 文/不壞的土叔 我叫張陵毕荐,是天一觀的道長。 經(jīng)常有香客問我艳馒,道長憎亚,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,218評論 1 292
  • 正文 為了忘掉前任弄慰,我火速辦了婚禮第美,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘陆爽。我一直安慰自己什往,他們只是感情好,可當我...
    茶點故事閱讀 67,234評論 6 388
  • 文/花漫 我一把揭開白布慌闭。 她就那樣靜靜地躺著别威,像睡著了一般第献。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上兔港,一...
    開封第一講書人閱讀 51,198評論 1 299
  • 那天庸毫,我揣著相機與錄音,去河邊找鬼衫樊。 笑死飒赃,一個胖子當著我的面吹牛,可吹牛的內(nèi)容都是我干的科侈。 我是一名探鬼主播载佳,決...
    沈念sama閱讀 40,084評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼臀栈!你這毒婦竟也來了蔫慧?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 38,926評論 0 274
  • 序言:老撾萬榮一對情侶失蹤权薯,失蹤者是張志新(化名)和其女友劉穎姑躲,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體盟蚣,經(jīng)...
    沈念sama閱讀 45,341評論 1 311
  • 正文 獨居荒郊野嶺守林人離奇死亡黍析,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,563評論 2 333
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了屎开。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片阐枣。...
    茶點故事閱讀 39,731評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖奄抽,靈堂內(nèi)的尸體忽然破棺而出蔼两,到底是詐尸還是另有隱情,我是刑警寧澤逞度,帶...
    沈念sama閱讀 35,430評論 5 343
  • 正文 年R本政府宣布额划,位于F島的核電站,受9級特大地震影響第晰,放射性物質(zhì)發(fā)生泄漏锁孟。R本人自食惡果不足惜彬祖,卻給世界環(huán)境...
    茶點故事閱讀 41,036評論 3 326
  • 文/蒙蒙 一茁瘦、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧储笑,春花似錦甜熔、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,676評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽啃洋。三九已至芦瘾,卻和暖如春烦磁,著一層夾襖步出監(jiān)牢的瞬間措左,已是汗流浹背陨帆。 一陣腳步聲響...
    開封第一講書人閱讀 32,829評論 1 269
  • 我被黑心中介騙來泰國打工锈颗, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留诵闭,地道東北人炼团。 一個月前我還...
    沈念sama閱讀 47,743評論 2 368
  • 正文 我出身青樓,卻偏偏與公主長得像疏尿,于是被迫代替她去往敵國和親瘟芝。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 44,629評論 2 354