【轉(zhuǎn)】單倍型基因組組裝方法

1. 什么是單倍型?

image

同源染色體:同源染色體,一個來自母本腔寡,一個來自于父本。

單倍型:單倍體基因型的簡稱掌唾。遺傳學(xué)上指在單條染色體上一系列遺傳變異位點的組合放前。

2. 單倍型組裝的意義?

目前糯彬,大多數(shù)二倍體基因組組裝都忽略了同源染色體之間的差異凭语,將基因組組裝成一個假的單倍體序列,這是二倍體類型的組裝的人為共識撩扒。這種人為的共識可能導(dǎo)致基因注釋的不精確和生物學(xué)解釋的錯誤似扔。

為了深入研究的需要吨些,更多的物種需要將來自父母的遺傳信息都獲得,因此參考基因組就需要獲得兩個單倍體基因組炒辉,也就是單倍型基因組豪墅。

目前單倍型技術(shù)主要應(yīng)用領(lǐng)域包括:

  • 在醫(yī)學(xué)上探索致病機理,挖掘致病基因黔寇,尋找疾病治療新方法偶器;
  • 在群體遺傳學(xué)上分析等位基因間差異,追蹤個體親緣關(guān)系缝裤,了解生物遷徙模式和進化歷史屏轰;
  • 在農(nóng)業(yè)上發(fā)掘優(yōu)異等位基因變異,探索雜種優(yōu)勢理論等憋飞。

3. 如何進行單倍型組裝霎苗?

早期已經(jīng)提出了幾種算法來生成單倍型解析的程序集,也稱為分階段程序集榛做。FALCON-Unzip叨粘,Supernova等使用相對短距離的序列數(shù)據(jù)進行定相,但只能解析高達9Mb的單倍型人類樣品瘤睹。這些方法無法逐步完成著絲粒或長重復(fù)答倡。擴展FALCON-UnzipFALCON-Phase使用Hi-C連接相控序列模塊轰传,可以生成更長的單倍型,但無法實現(xiàn)染色體長的定相瘪撇。

近年出現(xiàn)了幾種有效的單倍型組裝方法获茬。

方法1:Trio-binning (Illumina+Pacbio)

由美國國家人類基因組研究所、Pacific Biosciences公司及阿德萊德大學(xué)等單位的研究人員開發(fā)倔既,發(fā)表在2018年10月22日的Nature Biotechnology雜志上恕曲。

Trio binning首先使用來自兩個親本基因組的高精度短讀長數(shù)據(jù)將子代的長讀長序列劃分為單倍型特異性的集合,然后每個單倍型獨立組裝渤涌,形成一個完整的二倍體重建佩谣。

組裝方法

image
  • 1)測序:兩個親本分別二代Illumina測序,對F1代進行三代PacBio測序实蓬。
  • 2)分割三代數(shù)據(jù):使用兩個親本的二代數(shù)據(jù)獲取單倍型特異性k-mers茸俭。利用特異性k-mer將三代數(shù)據(jù)分割
  • 3)利用分割的數(shù)據(jù)分別組裝

優(yōu)缺點

Trio binning是一種簡便、準確安皱、高效的二倍體參考基因組組裝方法调鬓。在擬南芥、人類及牛單倍型組裝中表現(xiàn)良好酌伊,但Trio binning對樣本具有很高的要求腾窝,必須能夠獲取雙親的二代數(shù)據(jù)。

在進行數(shù)據(jù)分割時一部分雜合子reads不能明確地劃分為親本單倍型:如果雙親在某個位點上都是雜合,那么這個位點無法給reads提供有效的kmer信息虹脯,并且不能被唯一地分配給一個親本單倍型驴娃;同樣如果父本在一個位點是雜合子,而母本是純合的归形,從母本單倍型來看也不能分割托慨。在標準的trio-binning中,不能被區(qū)分的雜合reads在兩個親本數(shù)據(jù)集中都會使用暇榴。因此厚棵,這兩個等位基因可能存在于一個單倍型組合中,并引入錯誤蔼紧。另外還可能存在將reads錯誤劃分到其中一個親本的情況婆硬。

方法2:DipAsm(HiFi+Hi-C)

由李恒、Evan E. Eichler奸例、George M. Church等人聯(lián)合開發(fā)的新的基因組組裝方法彬犯,發(fā)表在2020年12月7日的Nat Biotechnol 雜志上。

DipAsm使用HiFi數(shù)據(jù)和Hi-C數(shù)據(jù)查吊,可以在1天之內(nèi)生成染色體規(guī)模的分相組裝谐区,具有98-99%的準確性。

組裝方法

image
  • 1)Peregrine:Pacbio HiFi reads 組裝成不分相位的congtig逻卖;
  • 2)HiRise / 3D-DNA:結(jié)合Hi-C數(shù)據(jù)生成不分相位的scaffold宋列;
  • 3)DeepVariant :Hi-C reads mapping到scaffold上call 雜合snp;
  • 4)WhatsHap和HapCUT2:處理HiFi reads和Hi-C reads雜合snp评也;
  • 5)WhatsHap:根據(jù)相位分割reads炼杖;
  • 6)Peregrine:分割的reads分別組裝。

優(yōu)缺點

DipAsm將促進高質(zhì)量的精準醫(yī)學(xué)以及個體單倍型變異和種群多樣性的研究盗迟,但DipAsm使用SNP信息進行定相坤邪,這對于長度長數(shù)據(jù)準確性要求高,也就是需要使用 PacBio HiFi罚缕,否則將增加SNP的錯誤率艇纺,部分涉及長SV的高度雜合區(qū)域會出現(xiàn)錯誤。

方法3:strand seq + long reads

由德國杜塞爾多夫海因里嫌实·海涅大學(xué)Tobias Marschall和美國華盛頓大學(xué)Evan E. Eichler合作喂饥,使用單細胞鏈測序和長讀取實現(xiàn)了親本數(shù)據(jù)非依賴的全階段人基因組組裝,2020年12月7日發(fā)表在Nature Biotechnology上肠鲫。

組裝方法:

Strand-seq具有三個重要功能:

  • 1)它可以按染色體對reads或contig進行排序员帮;
  • 2)它可以定序和定向contig;
  • 3)它提供了一個染色體范圍內(nèi)的相位信號导饲,而與物理距離無關(guān)捞高。
image

步驟:

    1. 使用長reads組裝成不分相位的contig氯材;
    1. Strand-seq比對到congtig上,contig進行排序和連接硝岗,形成染色體氢哮;
    1. 長reads mapping到 染色體上call SNVs;
    1. WhatsHap:利用SNVs 分割長reads型檀;
    1. wtdbg2 , Flye, Canu or Peregrine2:分割的reads分別組裝冗尤。

優(yōu)缺點:

組裝準確(質(zhì)量值> 40)且高度連續(xù)(contig N50> 23 Mbp)、轉(zhuǎn)換錯誤率低(0.17%)胀溺、并可提供了全相單核苷酸變體裂七、插入缺失和結(jié)構(gòu)變體等。

Strand-seq是一種單細胞技術(shù)仓坞,它不需要親本或配子背零,這種技術(shù)利用基因圖譜技術(shù)對染色體、單倍型和scaffold的長序列進行聚類无埃; 然而徙瓶,生成Strand-seq數(shù)據(jù)的困難限制了它在少數(shù)模型物種中的應(yīng)用。

轉(zhuǎn)自《百邁克》推文

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末嫉称,一起剝皮案震驚了整個濱河市侦镇,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌织阅,老刑警劉巖虽缕,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異蒲稳,居然都是意外死亡,警方通過查閱死者的電腦和手機伍派,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進店門江耀,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人诉植,你說我怎么就攤上這事祥国。” “怎么了晾腔?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵舌稀,是天一觀的道長。 經(jīng)常有香客問我灼擂,道長壁查,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任剔应,我火速辦了婚禮睡腿,結(jié)果婚禮上语御,老公的妹妹穿的比我還像新娘。我一直安慰自己席怪,他們只是感情好应闯,可當我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著挂捻,像睡著了一般碉纺。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上刻撒,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天骨田,我揣著相機與錄音,去河邊找鬼疫赎。 笑死盛撑,一個胖子當著我的面吹牛,可吹牛的內(nèi)容都是我干的捧搞。 我是一名探鬼主播抵卫,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼胎撇!你這毒婦竟也來了介粘?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤晚树,失蹤者是張志新(化名)和其女友劉穎姻采,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體爵憎,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡慨亲,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了宝鼓。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片刑棵。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖愚铡,靈堂內(nèi)的尸體忽然破棺而出蛉签,到底是詐尸還是另有隱情,我是刑警寧澤沥寥,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布碍舍,位于F島的核電站,受9級特大地震影響邑雅,放射性物質(zhì)發(fā)生泄漏片橡。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一淮野、第九天 我趴在偏房一處隱蔽的房頂上張望锻全。 院中可真熱鬧狂塘,春花似錦、人聲如沸鳄厌。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽了嚎。三九已至泪漂,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間歪泳,已是汗流浹背萝勤。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留呐伞,地道東北人敌卓。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓,卻偏偏與公主長得像伶氢,于是被迫代替她去往敵國和親趟径。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內(nèi)容