Nanopore測序筆記

原理:基于不同建立經(jīng)過納米孔時孕讳,電流信號的差異
典型應(yīng)用:

  • 讀長比較長匠楚,利于拼接
  • 解決動植物多倍體高重復(fù)、高雜合的特性
  • 也適用于細菌拼接
  • 對RNA直接測序厂财,無需反轉(zhuǎn)錄芋簿,直接識別RNA的堿基修飾
  • 適合檢測大片段結(jié)構(gòu)變異
  • 測序長度越長對定相分析越有利

測序設(shè)備

設(shè)備名稱 特點
minION 入門選擇,手持式璃饱,數(shù)據(jù)量30Gb
GridION 相當(dāng)于5張芯片外加一臺服務(wù)器
PromethION 24或48個測序芯片与斤,數(shù)據(jù)量達到Tb級別
MinION Mk1C 可連接藍牙或者Sim卡
Flongle MinION和GridION的轉(zhuǎn)換器
MinIT 預(yù)裝Linux系統(tǒng)
VolTRAX 自動化文庫制備

建庫測序
對待測DNA的一些處理,是一個模板過程。

  • DNA提取
    評估核酸質(zhì)量和完整性
    可以使用Nanodrop儀器進行檢驗撩穿。
    DNA :A260/A280 = ~1.8 磷支,A260/A230=~2.0-2.2
    RNA : A260/A280 = ~2.0 ,A260/A230=~2.0-2.2
    評估片段大小
    使用Agilent Bioanalyzer/Tapestation食寡,PFGE雾狈,fragment analyzer,F(xiàn)EMTO pulse
    nanopore建庫
    建庫圖

    給待測DNA兩側(cè)先加上A堿基抵皱,使平末端變成粘性末端善榛,然后再加上Y接頭以及馬達蛋白。
    測序錯誤
    堿基識別的非準確性呻畸,每秒讀數(shù)也不恒定移盆。信號讀取的頻率無法做到一個堿基讀一個,而是多個堿基一起讀一個信號伤为,相當(dāng)于一個分類器咒循。

nanopore測序文件可以生成fast5文件和fastq文件
fast5文件: 大而全,耗存儲大钮呀,適合做堿基校正和表觀遺傳學(xué)分析剑鞍。記錄著設(shè)備運行時全部的信息,包括捕獲的電信號值爽醋,設(shè)備運行時間蚁署,電壓,溫度等等信息蚂四。

fastq文件:常用,廠商一般給的遂赠,包含堿基和堿基質(zhì)量
流程
Basecalling:將原始測序文件轉(zhuǎn)換為堿基序列的過程
久妆,如guppy工具主要包括三大功能:
1、堿基識別
這是它的主要功能跷睦,Guppy使用神經(jīng)網(wǎng)絡(luò)算法(RNN)將電信號數(shù)據(jù)解讀為DNA或RNA的五個標準堿基筷弦,即腺嘌呤,鳥嘌呤抑诸,胞咳淀烂琴,胸腺咳淀和尿密淀。
2蜕乡、拆分條碼數(shù)據(jù)
如果建庫時添加了barcode條形碼奸绷,那么在序列開端和尾端均有Oxford Nanopore Technologies提供的條形碼序列,Guppy數(shù)據(jù)拆分功能基于識別出的堿基序列层玲。這個原理類似于illumina basecalling中的index拆分号醉。
3反症、比對
將測序數(shù)據(jù)與參考序列進行比對,其實這個是調(diào)用minimap2做的畔派,而且選項參數(shù)都是內(nèi)置的铅碍,無法調(diào)整,所以父虑,這個功能還是不要使用了该酗,一個軟件還是把自己該做的事情做好。
4士嚎、甲基化修飾位點檢測
如果某個堿基發(fā)生了甲基化修飾呜魄,那么在堿基識別過程中就會表現(xiàn)出一些異常信息,這些離群信息有可能就是潛在的甲基化修飾位點莱衩,guppy目前也可以進行堿基修飾的檢測爵嗅,不過這個需要學(xué)習(xí)集,不同物種需要單獨的學(xué)習(xí)集笨蚁,目前只試驗了大腸桿菌睹晒,人等模式生物,如果有需要可以嘗試進行檢測括细。
下載和安裝參考官網(wǎng)

使用minion_qc來進行統(tǒng)計繪圖

常用軟件

軟件名 功能 官方網(wǎng)址
bioconda 生物軟件商店 http://bioconda.github.io/
sratoolkit SRA數(shù)據(jù)管理 https://www.ncbi.nlm.nih.gov/sra/docs/toolkitsoft
flappie 堿基識別 https://github.com/nanoporetech/flappie
albacore 堿基識別 https://github.com/Albacore/albacore
guppy 堿基識別 nanopore社區(qū)
MINKNOW 測序操控軟件 nanopore社區(qū)
hdfview 查看fast5 https://www.hdfgroup.org/downloads/hdfview/
ont_fast5_api 拆分合并fast5 https://github.com/nanoporetech/ont_fast5_api
pomoxis 綜合分析 https://github.com/nanoporetech/pomoxis
minion_qc 質(zhì)控繪圖 https://github.com/roblanf/minion_qc
Filtlong 過濾 https://github.com/rrwick/Filtlong
seqkit 序列處理工具包 https://bioinf.shenwei.me/seqkit/
nanopack 數(shù)據(jù)質(zhì)控顧慮包 https://github.com/wdecoster/nanopack
NanoPlot 質(zhì)控繪圖 https://github.com/wdecoster/NanoPlot
Porechop 過濾 https://github.com/rrwick/Porechop
MaSuRCA 拼接 https://github.com/alekseyzimin/masurca
Tulip 拼接 https://github.com/Generade-nl/TULIP
spades 拼接 https://github.com/ablab/spades
canu 拼接 https://github.com/marbl/canu
NextDenovo 拼接 https://github.com/Nextomics/NextDenovo
RaGOO 參考引導(dǎo) https://github.com/malonge/RaGOO
Unicycler 拼接 https://github.com/rrwick/Unicycler
MECAT 拼接 https://github.com/xiaochuanle/MECAT
NECAT 拼接 https://github.com/xiaochuanle/NECAT
miniasm 拼接 https://github.com/lh3/miniasm
WTDBG2 拼接 https://github.com/ruanjue/wtdbg2
SMARTdenovo 拼接 https://github.com/ruanjue/smartdenovo
RA 拼接 https://github.com/rvaser/ra
gfatools gfa文件處理 https://github.com/lh3/gfatools
medaka 糾錯 https://github.com/nanoporetech/medaka
pilon 糾錯 https://github.com/broadinstitute/pilon/wiki
racon 糾錯 https://github.com/lbcb-sci/racon
NextPolish 糾錯 https://github.com/Nextomics/NextPolish
quast 拼接評估 http://bioinf.spbau.ru/quast
busco 拼接評估 https://gitlab.com/ezlab/busco
bwa 短序列比對 https://github.com/lh3/bwa
samtools sam格式處理 https://github.com/lh3/samtools
minimap2 比對 https://github.com/lh3/minimap2
NGMLR 比對 https://github.com/philres/ngmlr
last 比對 http://last.cbrc.jp/
GraphMap 比對 https://github.com/isovic/graphmap
MarginAlign 比對 https://github.com/benedictpaten/marginAlign
MUMMER 全局比對 http://mummer.sourceforge.net/
megan 物種鑒定 http://ab.inf.uni-tuebingen.de/software/megan6/download/
diamond blast比對 https://github.com/bbuchfink/diamond
seqtk 序列處理工具 https://github.com/lh3/seqtk
taxonkit 物種分類鑒定 https://bioinf.shenwei.me/taxonkit/
Centrifuge 物種分類 https://github.com/DaehwanKimLab/centrifuge
Nanopolish 綜合分析 https://github.com/nanoporetech/nanopolish
NanoSV 變異檢測 https://github.com/mroosmalen/nanosv
Sniffles 變異檢測 https://github.com/fritzsedlazeck/Sniffles
NanoVar 變異檢測 https://github.com/benoukraflab/NanoVar
tombo 甲基化 https://github.com/nanoporetech/tombo
bcftools vcf格式處理 https://github.com/samtools/bcftools
Truvari 變異檢測 https://github.com/spiralgenetics/truvari/
Pinfish 異構(gòu)體鑒別 https://github.com/nanoporetech/pinfish
pychopper 異構(gòu)體鑒別 https://github.com/nanoporetech/pychopper
tablet 可視化 https://ics.hutton.ac.uk/tablet/
IGV 可視化 http://software.broadinstitute.org/software/igv/

以上內(nèi)容整理至知乎https://zhuanlan.zhihu.com/p/102392867
原作者 公眾號 基因?qū)W苑

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末伪很,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子奋单,更是在濱河造成了極大的恐慌锉试,老刑警劉巖,帶你破解...
    沈念sama閱讀 219,589評論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件览濒,死亡現(xiàn)場離奇詭異呆盖,居然都是意外死亡,警方通過查閱死者的電腦和手機贷笛,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,615評論 3 396
  • 文/潘曉璐 我一進店門应又,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人乏苦,你說我怎么就攤上這事株扛。” “怎么了汇荐?”我有些...
    開封第一講書人閱讀 165,933評論 0 356
  • 文/不壞的土叔 我叫張陵洞就,是天一觀的道長。 經(jīng)常有香客問我拢驾,道長奖磁,這世上最難降的妖魔是什么改基? 我笑而不...
    開封第一講書人閱讀 58,976評論 1 295
  • 正文 為了忘掉前任繁疤,我火速辦了婚禮咖为,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘稠腊。我一直安慰自己躁染,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 67,999評論 6 393
  • 文/花漫 我一把揭開白布架忌。 她就那樣靜靜地躺著吞彤,像睡著了一般。 火紅的嫁衣襯著肌膚如雪叹放。 梳的紋絲不亂的頭發(fā)上饰恕,一...
    開封第一講書人閱讀 51,775評論 1 307
  • 那天,我揣著相機與錄音井仰,去河邊找鬼埋嵌。 笑死,一個胖子當(dāng)著我的面吹牛俱恶,可吹牛的內(nèi)容都是我干的雹嗦。 我是一名探鬼主播,決...
    沈念sama閱讀 40,474評論 3 420
  • 文/蒼蘭香墨 我猛地睜開眼合是,長吁一口氣:“原來是場噩夢啊……” “哼了罪!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起聪全,我...
    開封第一講書人閱讀 39,359評論 0 276
  • 序言:老撾萬榮一對情侶失蹤泊藕,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后荔烧,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體吱七,經(jīng)...
    沈念sama閱讀 45,854評論 1 317
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 38,007評論 3 338
  • 正文 我和宋清朗相戀三年鹤竭,在試婚紗的時候發(fā)現(xiàn)自己被綠了踊餐。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 40,146評論 1 351
  • 序言:一個原本活蹦亂跳的男人離奇死亡臀稚,死狀恐怖吝岭,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情吧寺,我是刑警寧澤窜管,帶...
    沈念sama閱讀 35,826評論 5 346
  • 正文 年R本政府宣布,位于F島的核電站稚机,受9級特大地震影響幕帆,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜赖条,卻給世界環(huán)境...
    茶點故事閱讀 41,484評論 3 331
  • 文/蒙蒙 一失乾、第九天 我趴在偏房一處隱蔽的房頂上張望常熙。 院中可真熱鬧,春花似錦碱茁、人聲如沸裸卫。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,029評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽墓贿。三九已至,卻和暖如春蜓氨,著一層夾襖步出監(jiān)牢的瞬間聋袋,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,153評論 1 272
  • 我被黑心中介騙來泰國打工穴吹, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留舱馅,地道東北人。 一個月前我還...
    沈念sama閱讀 48,420評論 3 373
  • 正文 我出身青樓刀荒,卻偏偏與公主長得像代嗤,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 45,107評論 2 356

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

  • 最近得知王通老師Nanopore測序技術(shù)的培訓(xùn)課在B站免費了,那一定是要跟進一波了废封。王老師在2019年下半年在公眾...
    周運來就是我閱讀 2,466評論 0 7
  • 軟件名功能官方網(wǎng)址biocondabioconda商店http://bioconda.github.io/[htt...
    wo_monic閱讀 2,138評論 0 7
  • 基因組組裝 基因組組裝一般分為三個層次卤档,contig, scaffold和chromosomes. contig表...
    xuzhougeng閱讀 30,214評論 3 122
  • 移步github 測序原理 共有的特點: 優(yōu)點:無需前期擴增,不引入偏向性;讀長長缺點:錯誤率高; 10X Gen...
    UnderStorm閱讀 23,168評論 2 54
  • 推薦指數(shù): 6.0 書籍主旨關(guān)鍵詞:特權(quán)、焦點渠鸽、注意力、語言聯(lián)想柴罐、情景聯(lián)想 觀點: 1.統(tǒng)計學(xué)現(xiàn)在叫數(shù)據(jù)分析徽缚,社會...
    Jenaral閱讀 5,721評論 0 5