NGS测序技术与分析流程图PPT课件.ppt

第一代测序技术Sanger测序1977年，桑格测定了第一个基因组序列，是噬菌体X174的，全长5375个碱基1.Sanger法核心原理是：由于ddNTP的2和3都不含羟基，其在DNA的合成过程中不能形成磷酸二酯键，因此可以用来中断DNA合成反应，在4个DNA合成反应体系中分别加入一定比例带有放射性同位素标记的ddNTP（分为：ddATP,ddCTP,ddGTP和ddTTP），通过凝胶电泳和放射自显影后可以根据电泳带的位置确定待测分子的DNA序列2.新一代测序介绍三大测序平台的前世今生Lynx MPSSSolexaABI SOLiD454Ion TorrentHelicosSMRTIllumina SolexaRoche 454Polony SeqABI Ion TorrentPacific Biosciences 3.Roche-454454公司可谓新一代测序技术的奠基人。2005年底，454公司推出了革命性的基于焦磷酸测序法的超高通量基因组测序系统GenomeSequencer20System，被Nature杂志以里程碑事件报道，开创了边合成边测序（sequencing-by-synthesis）的先河。之后，454公司被罗氏诊断公司以1.55亿美元收购。454测序原理4.测序实验流程：1、文库制备：根据样品的种类和实验目的，将基因组DNA/cDNA片段化处理至400-800bp间，经末端修复与特异性接头连接等修饰后变性处理回收单链的DNA（sstDNA）；然后和两个44个碱基长的衔接子（adaptor）A、B进行平端连接。A、B衔接子各自含有20个碱基的PCR引物序列、20个碱基的测序引物序列和4个碱基的对照序列（TACG），除此之外，B衔接子的5端还标记有一个生物素基团，供后续的分离合适的测序模板使用。5.测序实验流程：2、EmulsionPCR：特定比例的单链DNA文库被固定在特别设计的DNA捕获磁珠上，使大部分磁珠磁珠携带了一个独特的单链DNA片断。磁珠结合的文库被扩增试剂乳化，形成油包水的混合物，每个独特的片断在自己的微反应器里进行独立的扩增，而不受其他的竞争性或者污染性序列的影响。整个片段文库的扩增平行进行。扩增后产生了几百万个相同的拷贝。随后，乳液混合物被打破，扩增后仍结合在磁珠上的片段既可被回收纯化用于后续的测序实验；6.测序实验流程：3、测序反应：携带DNA的珠子与其他反应物混合物，随后放入PTP板中进行后继的测序。PTP孔的直径（29um）只能容纳一个珠子（20um）。然后将PTP板放置在GSFLX中，测序开始。每一个与模板链互补的核苷酸的添加都会产生化学发光的信号，并被CCD照相机所捕获；7.测序实验流程：4、数据分析：GSFLX系统在10小时的运行当中可获得100多万个读长，读取超过4-6亿个碱基信息，通过GSFLX系统提供两种不同的生物信息学工具对测序数据进行分析。8.454Pyrosequencing9.10.454的特点与主要应用读长较长，400600bp通量较低，1Run1M序列，400600Mb相对成本较高主要应用：denovo测序11.IlluminasolexaSolexa测序原理12.Illuminasolexa13.IlluminaSolexa桥式PCRdioldiol1st cycle denaturation1st cycle annealingdioldioln=35total1st cycle extensiondioldioldioldiol2nd cycle denaturation2nd cycle annealingdioldioldioldioldioldiol2nd cycle extension14.IlluminaSolexaBaseCalling123789456T T T T T T T G T T G C T A C G A T 15.Solexa的特点与主要应用读长较短，100150bp通量高，25G每天，120-150G每Run主要应用：RNA测序、表观遗传学研究16.ABISOLiDSOLiDSequencingbyOligoLigation/DetectionOligo连接测序：通过连接酶连接，再对oligo上荧光基团进行检测SOLiD 5500 xl17.ABISOLiD测序前期制备A 样品片段化磁珠连接B 乳化PCR3末端修饰C 磁珠富集转到测序玻片18.ABISOLiD测序原理19.ABISOLiD荧光结合和结果示例SRR029969.1 VAB_5551_12_381_F3 length=35T11.0203.3.1113211010332111302330201+SRR029969.1 VAB_5551_12_381_F3 length=35!36!8/8:!:!4626=(8.)43),(95,A.SOLiD Oligo荧光基团模式图B.SOLiD 测序结果示例（Color Space)20.21.SOLiD的特点与主要应用读长较短，50-75bp精度高，可达Q40通量高，20-30G每天，1Run可达120G主要应用：基因组重测序、SNP检测等22.IontorrentIonTorrent测序原理23.精度Examples:90%confidence(10%errorrate)=Q1099%confidence(1%errorrate)=Q2099.9%confidence(.1%errorrate)=Q3024.三种平台的技术差异平台454454SolexaSolexaSOLiDSOLiDPCR磁珠乳化PCR桥式PCR磁珠乳化PCR测序载体磁珠玻片玻片测序方式焦磷酸、荧光可逆终止物、荧光连接酶、荧光结果序列FastQFastQCSFastQ25.三种平台的效能参数差异平 台读长通量周期精度Solexa HiSeq2000Single-end:1 x 35 bpPaired-end:2 x 50 bpPaired-end:2 x 100 bp25 Gb/d1.5d4d8d50 bp 85%以上Q30100 bp 80%以上Q30SOLiD 5500 xlSingle-end:75 bpPaired-end:75 x 35 bpMate-pair:60 x 60 bp20 30 Gb/d1d/1lane7d/12 lane7d/12 laneQ40454 GS FLX400-600 bp400 600 Mb/Run10hQ2026.转录组测序测序流程全转录组总RNApolyT富集mRNA去除rRNANon-coding RNA转录组mRNARNA片段化、纯化、检测产量连接两端接头序列逆转录生成cDNA选择适当长度cDNA进行扩增纯化扩增产物，评估产量上机进行高通量测序27.转录组主要分析内容无参考序列转录组分析内容有参考序列转录组分析内容1 测序数据产量统计，数据成分和质量评估；2 Contig及Scaffold长度分布3 Unigene的长度分布和功能注释，GO分类，Pathway分析，差异表达分析4 蛋白功能预测与分类，差异表达基因GO富集和 Pathway富集分析。1 基本数据统计，比对参考序列2 序列在基因组上在分布3 测序深度分析、随机性评估和基因差异表达分析4 新基因预测，基因可变剪接鉴定和基因融合鉴定等。28.小RNA测序29.果蝇三种组织中MiRNA表达情况30.MiRNA表达模式分析31.新miRNA预测32.miRNA编辑情况分析与统计33.分析软件介绍质量控制软件(QualityControl)定位软件(Mappingthereads)已知基因（差异）表达评估软件(DifferentialExpression)新基因鉴定软件(Newgeneidentification)可视化展示软件(Virsulization)基因功能注释软件(GOtermorKEGGpathwayanalysis)34.数据格式示例Fastq格式ColorSpace格式35.Phredscore36.质量控制37.利用galaxy进行原始数据处理https:/main.g2.bx.psu.edu/38.GalaxyHistoryVGNFASTQhttps:/main.g2.bx.psu.edu/u/jjv5/h/unnamed-history-139.http:/222.200.187.83/nongpost/linux_lecture/msg.php40.知识回顾Knowledge Knowledge ReviewReview41.谢谢！放映结束感谢各位的批评指导！让我们共同进步42.