OrthoFinder如何尋找同源基因

OrthoFinder如何尋找同源基因，針對這個問題，這篇文章詳細(xì)介紹了相對應(yīng)的分析和解答，希望可以幫助更多想解決這個問題的小伙伴找到更簡單易行的方法。

OrthoFinder與其它工具比較的結(jié)果見下圖

1. 下載安裝

        https://github.com/davidemms/OrthoFinder

        可以下載之后源碼安裝，也可以用conda安裝：

conda install -y orthofinder

        ps：強(qiáng)烈建議安裝DIAMOND

2. 準(zhǔn)備輸入數(shù)據(jù)

        OrthoFinder所需的輸入數(shù)據(jù)很簡單，把每個物種的蛋白序列放進(jìn)單獨的fasta文件中，然后把這些fasta文件放到一個目錄下。fasta文件命名為對應(yīng)的物種名。

3. 運行

orthofinder -f Dataset_ directory

       如果你想更改線程數(shù)，使用-t參數(shù)即可修改。默認(rèn)的比對工具是DIAMOND，你也可以通過-S指定blast等其他工具。其他參數(shù)詳情可以運行 orthofinder -h 看到。

 4. 結(jié)果文件

        在結(jié)果文件夾中，Orthogroups文件夾里面有所有的同源基因信息，還貼心的單獨給出了單拷貝同源基因信息。Gene_Trees 和 Species_Tree 文件夾分別是單獨的同源基因構(gòu)建的tree已經(jīng)整合所有同源基因構(gòu)建的物種樹。

Note:

1. 可以指定-M參數(shù)來指定物種樹的構(gòu)建算法。OrthoFinder默認(rèn)的方法是STAG算法。STAG整合了所有的同源基因（包括多拷貝基因），這種方法特別適合物種遺傳距離較遠(yuǎn)，單拷貝同源基因很少甚至沒有的情況。

同樣，可以指定只使用單拷貝基因來構(gòu)建物種樹，-M raxml 就會調(diào)用Raxml進(jìn)行構(gòu)建。

2. Orthogroups文件夾中的 Orthogroups.GeneCount.tsv 統(tǒng)計了同源基因在每個物種中的數(shù)目，可以利用這個文件很方便的挑選我們需要的基因。

關(guān)于OrthoFinder如何尋找同源基因問題的解答就分享到這里了，希望以上內(nèi)容可以對大家有一定的幫助，如果你還有很多疑惑沒有解開，可以關(guān)注億速云行業(yè)資訊頻道了解更多相關(guān)知識。