生物信息学高通量测序技术和数据分析陈润生院士省公共课一等奖全国赛课获奖课件.pptx

《生物信息学高通量测序技术和数据分析陈润生院士省公共课一等奖全国赛课获奖课件.pptx》由会员分享，可在线阅读，更多相关《生物信息学高通量测序技术和数据分析陈润生院士省公共课一等奖全国赛课获奖课件.pptx（45页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、生物信息学陈小伟中国科学院生物物理研究所.09.30高通量测序技术及数据分析介绍第1页高通量测序技术及数据分析介绍背景介绍第一代测序技术第二代（高通量）测序技术基因芯片与高通量测序比较高通量测序技术应用高通量测序数据分析概览高通量测序数据质量评定与过滤基因组测序RNA-seqChIP-seqUCSC Genome Bioinformatics第2页背景介绍第3页背景介绍第一代测序技术Sanger测序法链终止法双脱氧终止法1975年Transcriptionhttp:/ or Illumina 1.8 and later.UsingaPhredscaleencodedusingASCII33to

2、93.ThisisthestandardforfastqformatsexceptfortheearlyIlluminadataformats(thischangedwithversion1.8oftheIlluminaPipeline).Illumina Pipeline 1.2 and earlier.UsingaSolexa/Illuminascale(-5to40)usingASCII59to104.TheWorkbenchautomaticallyconvertsthesequalityscorestothePhredscaleonimportinordertoensureacomm

3、onscaleforanalysesacrossdatasetsfromdifferentplatforms(seedetailsontheconversionnexttothesamplebelow).Illumina Pipeline 1.3 and 1.4.UsingaPhredscaleusingASCII64to104.Illumina Pipeline 1.5 to 1.7.UsingaPhredscaleusingASCII64to104.Values0()and1(A)arenotusedanymore.Value2(B)hasspecialmeaningandisusedas

4、atrimclipping.ThismeansthatwhenselectingIlluminaPipeline1.5andlater,thereadsaretrimmedwhenaBisencounteredintheinputfileiftheTrimreadsoptionischecked.363939393939393939393839393636343429312202019191938383836363636363630323535第9页基因芯片与高通量测序比较第10页芯片与测序比较基因芯片约历史，技术比较成熟，成本相对较低原理探针，互补配正确标准靶序列用荧光标识经过荧光强度间接反

5、应靶序列数量应用检测已知基因表示水平检测SNP位点基因型检测CNV第11页芯片与测序比较高通量测序约历史，发展快速，成本逐步降低原理边合成边测序碱基用荧光基团标识直接测定碱基序列应用全基因组测序转录组测序(smallRNAseq,RNA-seq)，能够检测已知基因表示水平，能够发觉全新转录本ChIP-seqCLIP-seq,第12页芯片与测序比较用高通量测序技术和基因芯片技术检测基因表示Malone,J.H.,andOliver,B.().Microarrays,deepsequencingandthetruemeasureofthetranscriptome.BMCBiol9,34.第13页

6、高通量测序技术应用第14页测序应用高通量测序数据分析概览第15页测序应用QualityAssessmentRawDataFastQC;fastx_quality_statsRemoveadaptor/linkerfastx_trimmerfastx_clipperSplitaccordingtobarcodefastx_barcode_splitter.plfastx_trimmerQualityControlfastq_quality_trimmerfastq_quality_filterFurtherAnalysis高通量测序数据质量评定与过滤FastQCFASTX-Toolkit第16页

7、测序应用全基因组de novo测序第一期：基因组调研图整体测序深度不低于20倍覆盖度。进行初步数据分析，对基因组大小，GC含量等做出初步评定，确定框架图梯度文库构建详细策略第二期：基因组框架图基因组覆盖度到达90%以上，基因区覆盖度到达95%以上，单碱基错误率到达1万分之一以内，整体测序覆盖深度不低于60倍覆盖度。同时对框架图进行基本基因注释和功效注释，和简单比较基因组学分析。第三期：基因组精细图基因组覆盖度到达95%以上，基因区覆盖度到达98%以上，单碱基错误率到达10万分之一以内，整体基因组覆盖度不低于100倍，ScaffoldN50大小不低于300Kb，对基因组精细图进行详细基因注释，基

8、因功效注释，基因代谢路径注释和比较基因组学分析。第17页全基因组denovo测序数据拼接组装算法流程DeBruijnGraph(德布鲁因图)Read：AGATACTk-merAGAGATATATACACTAGAGATATATACACT测序应用第18页全基因组重测序(外显子组测序)算法流程发觉遗传变异（SNP，indel等）测序数据与参考基因组做比对重新校对测序质量打分每一个基因型先验概率对基因型做推断计算每一个基因型概率测序应用第19页测序应用转录组测序Small RNA seq检测small RNA（主要是miRNA）表示水平发觉新small RNARNA-seqPoly(A)检测蛋白质编码

9、基因可变剪切体及表示水平Total RNA(except rRNA)检测mRNA及long noncoding RNA表示水平发觉新long noncoding RNA数据分析工具Bowtie(http:/bowtie- StevenSalzberg,UniversityofMarylandLiorPachter,UniversityofCalifornia,BerkeleyPostdoc JoinRinnslab,TheBroadInstitute第21页测序应用Overview of TopHat第22页测序应用SplicingJunctionsExon skippingorcassett

10、e exonMutually exclusive exonsAlternative donor siteAlternative acceptor siteIntron retention第23页测序应用TopHat:DiscoveringsplicejunctionsTopHatv1.0.7earlierseed-and-extendalignmentTopHatv1.0.7andlaterSupposeSisareadoflengthlthatcrossesasplicejunctionsplitsSintonsegments,n=floor(l/k),(k=25bp)mapsthesegm

11、entss1,snwithBowtietothegenomesegmentssi,si+1thatbothaligntothegenome,butnotadjacentlyasegmentsifailstoalignbecauseitcrossesasplicejunction,butsi-1andsi+1arealigned.第24页测序应用TopHat:DiscoveringsplicejunctionsTopHatv1.0.7andlaterasegmentsifailstoalignbecauseitcrossesasplicejunction,butsi-1andsi+1areali

12、gned.si-1Si+1si12345678910111213141516171819202122232425simbpk-mbpm=1,24m=12第25页测序应用OverviewofCufflinks第26页测序应用转录本拼接算法中包括到概念偏序关系与偏序集合 PartialorderandPartiallyorderedset偏序关系偏序（亦称半序）关系是定义在集合上一个序结构，是集合上满足一定条件二元关系。直观说，偏序指集合中仅有部分部分组员之间能够排序。全序关系在集合A中，存在偏序关系“”，假如对于任意aA,bA,有a b或b a，即A中每对元素都满足关系“”，则集合A上偏序“”是

13、全序或线性次序。直观来说，全序指集合中全体组员之间都能够进行比较，能够排出全部元素次序。偏序集合指配置了偏序关系集合第27页测序应用转录本拼接算法中包括到概念偏序关系非严格偏序，自反偏序给定集合S，“”是S上二元关系，若“”满足：自反性自反性：aS，有aa；反对称性反对称性：a，bS，ab且ba，则a=b；传递性传递性：a，b，cS，ab且bc，则ac；则称“”是S上非严格偏序或自反偏序严格偏序，反自反偏序给定集合S，“”是S上二元关系，若“”满足：反自反性反自反性：aS，有aa；非对称性非对称性：a，bS，abba；传递性传递性：a，b，cS，ab且bc，则ac；则称“Blat第41页测序应用UCSC Genome Bioinformaticshttp:/genome.ucsc.edu/查看特定序列在基因组上位置第42页测序应用UCSC Genome Bioinformaticshttp:/genome.ucsc.edu/查看特定序列在基因组上位置Zoom out 3x第43页生物信息学入门几点提议学习在线生物信息学惯用工具了解惯用数据库学习一门编程语言学习linux操作系统第44页谢谢！第45页