基因組組裝項目的十二步建議


1. 建立項目團體

多機構合作,數據和利益共享。

2. 收集目標基因組信息

考慮的因素:
基因組大小、倍性、雜合性、GC含量和重復。

數據庫查詢:
fungi (http://www.zbi.ee/fungalgenomesize)
animals (http://www.genomesize.com)
plants (http://data.kew.org/cvalues)

估計:
流式細胞儀和kmer頻率分布(建議兩種都用)。

3. 設計最佳實驗流程

高質量染色體水平的參考基因組是關鍵。
質控:reads長度、錯誤率、深度、覆蓋度、文庫等。

有錢:PacBio/ONT + Hi-C
沒錢:Illumina/10X GC(genomics chrominum) + Hi-C

從頭組裝:一般是完全denovo。
參考基因組輔助:利用近緣物種作為參考和指導進行組裝,該方法對數據和計算量較小,但是現有參考基因組可能有錯誤和重排。

目的:構建一致的單倍型或定相單倍型的染色體水平組裝。一般的組裝是將2條序列整合為1個單倍型,因此不能得到二倍體信息。

選擇合適的工具和流程:考慮組裝的質量和連續性,包括速度和敏感性。

三代組裝工具網站:
LRS-DB https://long-read-tools.org/

常用的組裝工具軟件:
image.png

4. 選擇最佳測序平台和准備文庫

文庫制備的兩個考慮:目標基因組大小、測序樣本數。

reads: 短(Illumina, 454, SOLiD, MGI, Ion Torrent),長(ONT and PacBio)或混合(hybrid) read

5. 選擇最佳DNA來源和提取方法

不含雜質。
最低量要求:
Illumina 和 10xGC > 3 ng, PacBio > 20 μg, ONT > 1 μg, BioNano > 200 ng, Dovetail > 5 μg 。
三代平均DNA長度>25 kb。
使用核與細胞器DNA比率更高的組織。
純化DNA的測量/定量可使用分光光度法和基於熒光的方法。

6. 檢查計算資源與要求

數據量、基因組大小、雜合率和倍性等對內存
需求、CPU數量和計算成本成幾何增加。
可選擇雲計算合理分配。

7. 選擇最佳計算設計和流程

三種選擇:
(1)最大化內部員工或協作
(2)從服務外包提供者
(3)模擬具有不同設置的數據

8. 基因組組裝

推薦的基因組組裝和注釋流程圖:
image.png

強烈建議使用BioNano和Hi-C數據來達到染色體級組裝,因為這兩種方法可通過驗證初始組裝的完整性,糾正方向錯誤,排序scaffolds來完善結果。

9. 在注釋前檢查組裝質量

在鳥槍法時代,denovo依賴於於算法和試驗設計。reads長度、文庫大小、reads准確性和基因組復雜性等決定了組裝的准確性和連續性。

質量評估:

  • 組裝大小
  • 組裝連續性(N50,NG50,NA50,NGA50)
  • 重疊群contig數目和(平均)長度
  • 組裝可能性得分(通過reads比對每一個候選組裝來計算)
  • 組裝完整度(BUSCO得分或RNAseq mapping)
  • 其他:QTL、ESTs、熒光原位雜交、BAC克隆、染色體水平遺傳圖譜。

三個最重要的指標:連續性、准確性、完整性

方法:三代/10XGC,BioNano,Hi-C數據;軟件LR_Gapcloser。

10. 基因組注釋

注釋內容:

  • 識別非編碼區:重復序列、轉座子。
  • 識別編碼區(稱為基因預測):內含子、外顯子、CDS、5/3 UTR。
  • 附加這些元素的生物學信息。

注釋的方法:

  • 手動注釋:耗時昂貴,需要獲得准確的基因模型和基因集。
  • 自動注釋:置信度和可靠性低(通常基於直系同源物種,不同數據庫數據不同)。
  • 半自動注釋:集成不同的結果獲得一致的注釋,平衡了手動和自動方法。

結合比對EST、RNAseq、蛋白序列作為外部基因組組裝證據。

結合方法和結果(尤其是MAKER,BRAKER和String-Tie)可以有效地提高注釋預測的數量和准確性(尤其是對孤兒基因和其他年輕基因)。

功能注釋GO等。

在線基因組注釋工具:
image.png

命令行注釋工具:
image.png
image.png

非編碼RNA注釋:
image.png

重復序列注釋:
image.png

11. 建立一種可查詢和可共享的輸出格式

公共數據庫 or 自建數據庫?

12. 分發社區來優化組裝和注釋

不同版本軟件結果不同,為確保穩定,數據可重復,需持續維護和更新。

植物社區示例:
https://nbenth.com/annotator/index,
https://solgenomics.net
https://www.helmholtz-muenchen.de/pgsb

動物社區示例:
http://www.slimsuite.unsw.edu.au/servers/apollo.php
https://bovinegenome.elsiklab.missouri.edu
http://www.gmgi.org/genomics-fish-shellfish
https://www.sanger.ac.uk/science/data/vertebrate-genomes-sequencing

對於初學者的基因組組裝和注釋流程的建議

不建議純二代組裝。
純三代或混合組裝方法:
image.png

此文太多廢話,慎讀~~~

文獻來源: Hyungtaek JungID et al. Twelve quick steps for genome assembly and annotation in the classroom. PLoS Comput Biol. 2020 Nov 12;16(11):e1008325.


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM