基因芯片分析的理论与方法

基因芯片分析的理论与方法第1页，课件共72页，创作于2023年2月基因芯片概论第2页，课件共72页，创作于2023年2月基因芯片分析的概念基因芯片(Genechip)技术是指通过微阵列(Microarray)技术将高密度DNA片段阵列通过高速机器人或原位合成方式以一定的顺序或排列方式使其附着在如玻璃片等固相表面，以荧光标记的DNA探针，借助碱基互补杂交原理，进行大量的基因表达及监测等方面研究的技术。第3页，课件共72页，创作于2023年2月基因芯片分析的优点快速高通量(104-106)自动化使用的试剂少低成本第4页，课件共72页，创作于2023年2月基因芯片分析的应用范围AcademicresearchofgeneticdiseasesCancerPrenatalgenetics（产前遗传学研究）GeneralgeneticdiseasesInfectiousdiseasesDrugdiscoveryAnimalfarming/veterinaryIndustrial(fermentation)Environmental第5页，课件共72页，创作于2023年2月发展的方向:从多角度研究揭示生命现象Stages:developmental,transformation,timeftertreatment,etc.Treatments(stimuli,drugs,nutrients,etc.)Physiologicalstates(stressed,fasting,etc.)Tissuedistribution(position,3D)第6页，课件共72页，创作于2023年2月基因芯片的发展是推动系统生物学发展的动力第7页，课件共72页，创作于2023年2月利用基因芯片研究生命现象的测略Factorsinvolved=>ComponentsOrderofevents=>PathwaysInteractions=>Circuit

KEGGGO第8页，课件共72页，创作于2023年2月基因芯片的分类根据用途分类

geneexpressionpatternBiologicalSampleFunctional

Information第9页，课件共72页，创作于2023年2月基因芯片的分类Oligonucleotidearray–Synthesizedonachip（

Affymetrix）–Spotonasolidmatrix（

Compugen）cDNAarray（

Incyte）根据探针类型分类

expressiongenomicanalysiscDNA-Chip GenomicChip<2,000n>50,000n第10页，课件共72页，创作于2023年2月一些发展中的基因芯片技术平台利用生物分子的电物理特性进行基因表达监测：监测速度很快，适用于基因表大，蛋白质组及基因型的研究利用电场原理进行高密度芯片生产：基于适合用于生物学的集成电路，集成电路包含可以独立寻址的微电极阵列，结合特殊的液体流动系统，可以使大部分生物分子按照来自于计算机的数字指令运动。喷墨点样技术：以高度定位的形式把合成好的寡核苷酸分子喷点倒玻璃表面。寡核苷酸包被的微珠芯片平行信号测序技术：对基因表达进行定量分析第11页，课件共72页，创作于2023年2月基因芯片分析试验方法第12页，课件共72页，创作于2023年2月基因芯片分析的主要步骤第13页，课件共72页，创作于2023年2月cDNA基因芯片分析的主要步骤第14页，课件共72页，创作于2023年2月cDNA芯片分析的主要步骤SpotbyArrayspotter第15页，课件共72页，创作于2023年2月cDNA芯片分析的主要步骤HybridizingbyAutomatichybridizationprocessor第16页，课件共72页，创作于2023年2月cDNA芯片分析的主要步骤Laserscanner第17页，课件共72页，创作于2023年2月Oligonucleotidearray（GeneChip）LLLLLLLLcDNAAAAA总RNA的制备反转录体外转录生物素标记的cRNA片段化处理带标记的cRNA片断35-200bases0.5-2ug/ul起始用量5-10ug（IVT）操作流程（以真核生物为例）第18页，课件共72页，创作于2023年2月LLLLL标记的cRNA片断杂交混合液的制备EukaryoticHyb.ControlControlOligoB2杂交（16hour）数据分析扫描洗脱染色第19页，课件共72页，创作于2023年2月Oligonucleotidearray的特点1个平方厘米的面积至少可排列四十多万个探针合成区（“点”）基因2基因1cDNA基因2cDNA用于cDNA芯片的探针Oligoprobe基因1多个检测结果可以参考第20页，课件共72页，创作于2023年2月Oligonucleotidearray的优越性序列准确性高原位合成PCR扩增,点样起始TotalRNA1－10ug≥50ug均一的退火温度25mer300bp－3Kb特异性更高多段探针单个探针非特异性杂交≤2％30％cDNA芯片Affy芯片优势第21页，课件共72页，创作于2023年2月芯片分析数据的标准化Quantitation

data-qualityassessments背景处理：图像上各点的吸光度值包含了样品和背景信号，在提取数据前必须将背景扣除杂交点质量：由于点样或膜变形等原因目前较多的软件对杂交点的识别定位仍需要人为的调整数据的标准化：其目的是避免基因芯片实验中因系统差异造成芯片间数据比较的困难。最常用的是“看家基因”法，它预先选择一组表达水平不变的看家基因，计算出这组基因平均ratio值为1时的标准化系数，然后将其应用于全部的数据以达到标准化的目的第22页，课件共72页，创作于2023年2月ScatterplotsforaCy5-liver/Cy3-livercontrolhubridizationACy5-ovariantumor/Cy3-normalovaryhubridizationThevalueofCy3andCy5hydridizationsignalsfromeachclonewereplotteddirectlyontotheplot（卵巢癌）DataAnalysis-RevealtheDifference第23页，课件共72页，创作于2023年2月基因芯片数据分析和信息挖掘方法第24页，课件共72页，创作于2023年2月一.聚类分析第25页，课件共72页，创作于2023年2月聚类分析聚类分析是模式识别中一种非常有吸引力的方法，特别适用于模式分类数不知道的情况。从机器学习的角度来看，有两种基本的聚类分析：有教师聚类无教师聚类第26页，课件共72页，创作于2023年2月基因表达数据聚类分析一般包括以下几个步骤：（１）确定基因表达的数据（２）计算相似性矩阵，各个矩阵元素代表两个基因的表达是否相似（３）选择算法进行聚类分析（４）显示分析结果。第27页，课件共72页，创作于2023年2月对数据进行聚类分析之前，必须将包含在基因表达矩阵中的数据进行相似程度分析，并且对分析结果进行量化。通常情况下，相似往往被赋于一个较大的量化的值，而不相似则由一个较小的量化的值来表示。在实际计算中，往往以距离代替相似的概念，相似性度量被转化为两个基因表达模式之间的距离。距离越小，表达模式越相近，反之，则表达模式差异大。第28页，课件共72页，创作于2023年2月两个表达模式之间的关系(a)相似(b)变化趋势一致(c)两个基因的调控结果不一样或甚至相反第29页，课件共72页，创作于2023年2月聚类分析的目的可诱导基因是共表达的许多构成性表达的基因不受调节根据表达相似形排列基因功能相似的基因被聚为不同的类可以揭示细胞的生理状态可以帮助研究未知基因的功能第30页，课件共72页，创作于2023年2月聚类结果显示:

Cluster,Clusterviewer第31页，课件共72页，创作于2023年2月二.基因调控网络分析第32页，课件共72页，创作于2023年2月2.1基本理论第33页，课件共72页，创作于2023年2月定义一个基因网络(pathway)由一组生物分子（如基因、蛋白质）以及它们之间的相互作用构成，这些生物分子共同完成一些特定的细胞功能任务在实际分析过程中，往往以图这种数据结构表示基因网络，图中的节点代表基因或者蛋白质，而节点之间的连线代表基因、蛋白质之间的相互作用。基因网络描述了特定细胞或组织中的功能路径，如代谢、基因调控，信号传导等。第34页，课件共72页，创作于2023年2月原理基因表达实际上是细胞、组织、器官受遗传和环境影响的结果。一个基因的转录和表达由细胞的生化状态所决定，在一个基因的转录过程中，一组转录因子作用于该基因的启动子区域，控制该基因转录，而这些转录因子本身又是其它基因的产物。当一个基因通过转录、翻译形成功能基因产物后，它将改变细胞的生化状态，从而直接或间接地影响其它基因的表达，甚至影响自身的表达。多个基因的表达不断变化，使得细胞的生化状态不断地变化。第35页，课件共72页，创作于2023年2月一个基因的表达受其它基因的影响，而这个基因又会影响其它基因的表达，这种相互影响、相互制约关系构成了复杂的基因表达调控网络。基因表达数据之中隐含基因之间的相互作用关系，因而可以通过分析基因表达数据，构建基因调控网络第36页，课件共72页，创作于2023年2月应用BetterunderstandthephysiologyoftheorganismPredicttheeffectofadrug–PredicttheweaknessesorsideeffectsofadrugPredictgoodorbaddrugcombinationsUnderstanddiseaseprocessesGenefunctionassignmentMakea“perfect”simulationofcellularfunctiontouseasanexperimentalmodelManytools“feed”intothisunderstandingincludingmicroarrays第37页，课件共72页，创作于2023年2月2.2pathway数据库介绍第38页，课件共72页，创作于2023年2月Pathway数据库GOKEGGGENMAPPBIACARTAGeneNet细胞因子网络数据库自由基信号数据库第39页，课件共72页，创作于2023年2月关于GOGO(geneontology)是基因本体论联合会(GeneOntologyConsortium)所建立的数据库，旨在建立一个适用于各种物种的，对基因和蛋白功能进行限定和描述的，并随着研究的不断深入而更新的语言词汇标准。GO是多种生物学本体论语言中的一种，提供了三层结构的系统定义方式，用于描述基因产物的功能。GO可以被用来在小鼠基因组中查询和信号转导相关的基因产物第40页，课件共72页，创作于2023年2月GO发展了具有三级结构的标准语言（ontologies）。根据基因产物的相关分子功能，生物学途径，细胞学组件而给予定义，无物种相关性。第41页，课件共72页，创作于2023年2月关于KEGGKEGG(KyotoEncyclopediaofGenesandGenomes)（京都大学基因和基因组百科全书数据库）KEGG是系统分析基因功能的数据库，将基因组的信息与基因功能联系起来，旨在揭示生命现象的遗传与化学蓝图第42页，课件共72页，创作于2023年2月数据库KEGGPATHWAYDatabaseKEGGGENESDatabaseKEGGLIGANDDatabase用途搜寻pathway,产生可能的反应路径搜寻类似的基因序列搜寻类似的基因组序列搜寻类似的复合物结构,类似的多糖结构及类似的反应类别通路信息基因组信息化学信息第43页，课件共72页，创作于2023年2月KEGG中的通路分为五大类：

1．新陈代谢通路

2．遗传信息处理通路

3．环境信息处理通路

4．细胞内通路

5．人类病症相关通路PATHWAY提供所选定通路的参考资料、涉及的反应和在其他生物中的通路情况。由GIF等平板文件类型图示，方框显示为涉及的酶EC名称，圆滑框为反应类型，以实线和箭头连接反应物和方向，虚线指向预测的反应类型。第44页，课件共72页，创作于2023年2月第45页，课件共72页，创作于2023年2月关于GenMappGemapp(http://www.GenMAPP.org)

GenMAPP是一个以学术研究为基础的科学团体，提供100个以上的通路图.SiliconGenetics公司已经和GenMAPP组织达成了协议，把这些通路整合到了GeneSpring的通路库.利用GeneSpring分析的基因芯片数据可以用GenMAPP的通路可视化，使研究者研究某些基因在通路中的作用.第46页，课件共72页，创作于2023年2月第47页，课件共72页，创作于2023年2月关于BioCartaBioCarta:“Carta”中文意思是”地图”，BioCarta”即用图形来描述生物相关知识，形成生物学家所熟悉的生化调控通路(biochemicalpathways)。BioCarta公司为客户免费提供通路绘制软件,让他们根据自己的研究领域发表通路,建成数据库,通路中的蛋白都是按照HUGO的标准建立的.BioCarta已收集各物种共大約120,000基因及136個调控通路，数据库载不断增加，其信号传导通路(signaltransductionpathways)是当前最全的.第48页，课件共72页，创作于2023年2月第49页，课件共72页，创作于2023年2月GeneNet真核生物生理过程的调控网络(GeneNet)GeneNet展示了协调方式下基因的整体功能，良好的调控功能和对外界刺激的反应。该基因网具有以下基本功能：(a)描述执行特定生物功能时所涉及到的整体基因交互作用(b)描述基因的蛋白质编码(c)描述基因对外界刺激的传感途径(d)通过一组反馈自动地稳定网络参数或者使网络转换到新的功能状态(e)利用外部信号、激素和代谢产物作为生理学参数激发基因网纠正其作用第50页，课件共72页，创作于2023年2月该基因网描述了四类实体：

(a)细胞(组织，器官)(b)蛋白质

(c)基因

(d)物质该基因网还描述了实体间两种关系：(a)反应（reaction）,即通过交互作用产生新的实体或过程(b)调控（regulatory）事件,特定反应对实体的作用。第51页，课件共72页，创作于2023年2月第52页，课件共72页，创作于2023年2月2.3功能富集分析第53页，课件共72页，创作于2023年2月功能富集分析定义功能富集分析(FunctionalEnrichmentAnalysis)

又称功能聚类分析，借助于各种生物学信息数据库（如GO,KEGG,GENMAPP,BIACARTA,TRANSFAC,OMIM）和分析工具（如MAPPFinder,ArrayXPath）进行统

计分析，挖掘同差异表达或共表达这组感兴趣基因具有显著差异的功能类别。功能富集分析的统计原理是用超几何分布型来检验一组基因（共表达或差异表达）中某个功能类的显著性，并结合多重假设检验判断思想选择同这组基因显著相关功能类别

第54页，课件共72页，创作于2023年2月MAPPFinderMAPPFinderisatoolthatcreatesaglobalgene-expressionprofileacrossallareasofbiologybyintegratingtheannotationsoftheGeneOntology(GO)ProjectwiththefreesoftwarepackageGenMAPPhttp://www.GenMAPP.org.可以整合GO定义，产生整体的基因表达谱.输出可以搜索的浏览器,帮助用户快速识别代表性的差异表达基因在GO的定义。第55页，课件共72页，创作于2023年2月工作方式第56页，课件共72页，创作于2023年2月Z值得计算第57页，课件共72页，创作于2023年2月第58页，课件共72页，创作于2023年2月GOminerGominer:最初Version在算法上虽亚于GenMAPP,在再建立可视关系上(tree-likestructure和‘directedacyclicgraph)有其独特之处.今年对OriginalVersion进行较大改进,不仅可以富集significantGOcategories,还可以同时对多套芯片实验数据批量分析,控制假发现率FDR,还整合转录因子结合位点信息.综合而言,尤其适合TIME-COURSE功能富集分析需要建立当地Mysql数据库,建立JDB数据源

第59页，课件共72页，创作于2023年2月第60页，课件共72页，创作于2023年2月PathwayExplorerPathwayExplorer:providescomprehensiveandeasilyaccessiblerepresentationsofexpressionprofilesontomajorregulatory,metabolicandcellularpathways.TheintegratedpathwayresourcesincludeKEGG,BioCartaandGenMAPP.第61页，课件共72页，创作于2023年2月LocusLinkwasagainusedasrootidentifier.TheLocusLinksarelinkedwiththeuser-definedgeneidentifiergroups(UniGene,GeneOntology,GenBankand/orRefSeq),whichareusedthentoalignthemappedgeneIDs.https://pathwayexplorer.genome.tugraz.at第62页，课件共72页，创作于2023年2月第63页，课件共72页，创作于2023年2月ArrayXPathaweb-basedserviceformappingandvisualizingmicroarray

gene-expressiondataforintegratedbiologicalpathwayresources

Whenoneinputsgene-expressionclusters,ArrayXPath

producesalistofthebestmatchingpathwaysforeachcluster

appliedFisher'sexacttestandthefalsediscoveryrate

(FDR)toevaluatethestatisticalsignific