1. 程式人生 > 實用技巧 >GWAMA:gwas meta-analysis工具

GWAMA:gwas meta-analysis工具

下載和安裝

從官網(https://genomics.ut.ee/en/tools/gwama-download)上下載或:

wget http:/www.geenivaramu.ee/tools/sample.zip

unzip gwama.zip (解壓)

make (安裝)

執行命令:

../GWAMA --filelist gwama.in --output gwama_xp_wz_jp --quantitative

其中,

gwama.in為輸入檔案,內容如下

每一行對應一個人群的gwas分析結果,需要的檔案格式如下:

MARKERNAME代表SNP位點的名稱,EA代表effect allele, NEA

代表other allele,EAF代表EA的頻率, 對於連續型的性狀,要求輸入BETASE兩列,對於二分類的性狀,要求輸入OR, OR_95L, OR_95U3列,代表odd ratio以及95%置信區間。

以上是最基礎的資料,如果有額外的資料也可以加入,比如STRAND代表鏈的方向,IMPUTED代表該位點是否是填充得到的等

gwama_xp_wz_jp為輸出檔案的字首,輸出結果儲存在gwama_xp_wz_jp.out檔案中,內容如下

--quantitative 表示使用的表型資料是連續變數。

其它需要請看官方文件:https://genomics.ut.ee/en/tools/gwama-tutorial

為什麼需要做meta分析?

群體分層是GWAS研究中一個比較常見的假陽性來源。

具體來說,GWAS的核心在於SNP與表型的關聯,具體來說,我們觀測到每一個genome位點會有許多不同的SNP,如果某個SNP總是與某種表型/疾病同時出現,那麼我們可以推測這個SNP極有可能會導致這個表型/疾病。如果存在群體分層,那麼發現的SNP可能是群體的一個特徵而不是與表型相關,也就是假陽性位點。

當群體出現分層時,常規手段時將分層的群體獨立分析,最後做meta分析。

如何判斷群體分層?

用plink計算PCA,畫PC1和PC2在不同群體的散點圖,看是否將人群分開,如果能觀察到幾個group,就是出現群體分層。

什麼是meta分析?

GWAS meta analysis應該翻譯成整合分析,略區別於醫學的薈萃分析。GWAS中的meta analysis,是將不同的GWAS整合在一起,減少群體分層的影響,提升了power。(還是不知道原理)

一般用的軟體有:plink、GWAMA、METAL

關於Meta的在不同場景下的意思

Meta的字面意思:Meta means about the thing itself. It's seeing the thing from a higher perspective instead of from within the thing, like being self-aware

Metadata is "data that provides information about other data"

在計算機裡,meta就是描述資料的資料。

META是一種思想概念,一種抽象思維,用來描述資料的資料,比如有一張學生表,記錄著學生的基本資訊,我們通過表可以獲取學生資訊(資料),但是有時候也要得到表本身的資訊資料(比如表結構資訊:欄位名稱,欄位資料型別,長度等資訊),對於這種基礎資訊的描述,就會使用META的概念,使用META元資料來描述表本身。放到HTML中也是一樣的,HTML用來描述網頁資訊,但是HTML自己也有一些資訊(比如網頁標題,網頁描述,搜尋關鍵字),這些資訊也就稱之為HTML META資訊,並且HTML也定義了專門的META標籤。

在醫學裡,meta就是薈萃分析。

Meta分析是對具有相同目的且相互獨立的多個研究結果進行系統的綜合評價和定量分析的一種研究方法。它不僅僅是蒐集很多以往的研究成果,同時進行全面系統的質量評價,另外還要對符合條件的研究進行定量分析。

在微生物裡,巨集基因組學(Metagenomics),又稱元基因組學。是以特定生境中的整個微生物群落作為研究物件,無需分離培養,直接提取環境樣本DNA進行測序,研究環境微生物的群落結構、物種分類、系統進化、基因功能及代謝網路等,已廣泛應用於微生物領域。

Reference:

https://genomics.ut.ee/en/tools/gwama

https://www.cnblogs.com/chenwenyan/p/10912521.html

https://www.cnblogs.com/leezx/p/13369846.html