1、中科院生物信息学期末考试复习题陈润生老师部分:1.什么是生物信息学,如何理解其含义?为什么在大规模测序研究中,生物信息学至关重要?答:生物信息学有三个方面的含义:1) 生物信息学是一个学科领域,包含着基因组信息的获取、处理、存储、分配 、分析和解释的所有方面,是基因组研究不可分割的部分。2) 生物信息学是把基因组 DNA 序列信息分析作为源头,破译隐藏在 DNA 序列中的遗传语言,特别是非编码区的实质;同时在发现了新基因信息之后进行蛋白质空间结构模拟和预测;其本质是识别基因信号。3) 生物信息学的研究目标是揭示“基因组信息结构的复杂性及遗传语言的根本规律”。它是当今自然科学和技术科学领域中“基
2、因组、“信息结构”和“复杂性”这三个重大科学问题的有机结合。生物信息学是把基因组 DNA 序列信息分析作为源头,找到基因组序列中代表蛋白质和 RNA基因的编码区;同时阐明基因组中大量存在的非编码区的信息实质,破译隐藏在 DNA 序列中的遗传语言规律:在此基础上,归纳、整理与基因组遗传信息释放及其调控相关的转录谱和蛋白谱数据,从而认识代谢、发育、分化、进化的规律。同时在发现了新基因信息之后,其还利用基因组中编码区信息进行蛋白空间结构模拟和蛋白功能预测,并将此类信息与生物体和生命过程中的生理生化信息结合,阐明其分子机制,最终进行蛋白、核酸分子设计、药物设计、个体化医疗保健设计。2.如何利用数据库信
3、息发现新基因,基本原理?答:利用数据库资源发现新基因,根据数据源不同,可分 2 种不同的查找方式:1) 从大规模基因组测序得到的数据出发,经过基因识别发现新基因:(利用统计,神经网络,分维,复杂度,密码学,HMM,多序列比对等方法识别特殊序列,预测新 ORF。但因为基因组中编码区少,所以关键是“数据识别”问题。)利用大规模拼接好的基因组,使用不同数据方法,进行标识查找,并将找到的可能的新基因同数据库中已有的基因对比,从而确定是否为新基因。可分为:基于信号,如剪切位点、序列中的启动子与终止子等。基于组分,即基因家族、特殊序列间比较,Complexity analysis,Neural Netwo
4、rk2) 利用 EST 数据库发现新基因和新 SNPs:(归属于同一基因的 EST 片断一定有 overlapping,通过 alignment 可组装成一完整的基因,但 EST 片断太小,不存在数据来源,主要是拼接问题)数据来源于大量的序列小片段,EST 较短,故关键在正确拼接。方法有基因组序列比对、拼接、组装法等。经常采用 SiClone 策略。其主要步骤有:构建数据库;将序列纯化格式标准化;从种子库中取序列和大库序列比对;延长种子序列,至不能再延长;放入contig 库构建若干数据库:总的纯化的 EST 数据库,种子数据库,载体数据库,杂质、引物数据库,蛋白数据库,cDNA 数据库;用所
5、用种子数据库和杂质、引物数据库及载体数据库比对,去除杂质;用种子和纯化的 EST 数据库比对用经过一次比对得到的长的片段和蛋白数据库、cDNA 数据库比较,判断是否为已有序列,再利用该大片段与纯化的 EST 数据库比对,重复以上步骤,直到序列不能再延伸;判断是否为全长 cDNA 序列。(利用 EST 数据库:原理:当测序获得一条 EST 序列时,它来自哪一个基因的哪个区域是未知的(随机的),所以属于同一个基因的不同 EST 序列之间常有交叠的区域。根据这种 “交叠”现象,就能找出属于同一个基因的所有 EST 序列,进而将它们拼接成和完整基因相对应的全长 cDNA 序列。而到目前为止,公共EST
6、 数据库(dbEST)中已经收集到约 800 万条的人的 EST 序列。估计这些序列已覆盖了人类全部基因的 95%以上,平均起来每个基因有 10 倍以上的覆盖率。)3.用蛋白或核酸序列数据库研究生物演化的主要步骤是什么?当前的困难是什么,如何克服?(核酸或氨基酸序列进行进化研究要进行哪些计算步骤?当前遇到什么问题?怎样解决?)答:计算步骤,构建系统进化树,其主要步骤如下:1) 序列相似性比较。就是将待研究序列与 DNA 或蛋白质序列库进行比较,用于确定该序列的生物属性,也就是找出与此序列相似的已知序列是什么。完成这一工作只需要使用两两序列比较算法。常用的程序包有 BLAST、FASTA 等;2
7、) 序列同源性分析。是将待研究序列加入到一组与之同源,但来自不同物种的序列中进行多序列同时比较,以确定该序列与其它序列间的同源性大小。这是理论分析方法中最关键的一步。完成这一工作必须使用多序列比较算法。常用的程序包有 CLUSTAL 等;3) 构建系统进化树。根据序列同源性分析的结果,重建反映物种间进化关系的进化树。为完成这一工作已发展了多种软件包,如 PYLIP、MEGA 等;4) 稳定性检验。为了检验构建好的进化树的可靠性,需要进行统计可靠性检验,通常构建过程要随机地进行成百上千次,只有以大概率(70 以上)出现的分支点才是可靠的。通用的方法使用 Bootstrap 算法。【(1.序列相似
8、性比较: 就是将待研究序列与 DNA 或蛋白质序列库进行比较,用于确定该序列的生物属性,也就是找出与此序列相似的已知序列是什么,完成这一工作只需要使用两两序列比较算法。常用的序列包有 BBLAST、FASTA 等;(2. 序列同源性分析:将待研究序列加入到一组与之同源,但来自不同物种的序列中进行多序列同时比较,以确定该序列与其他序列间的同源性大小,这是理论分析方法中最关键的一步,完成这一工作必须使用多序列比较算法,常用的程序包有 CLUSTAL 等;(3.构建系统进化树:根据序列同源性分析的结果,重建反应物种间进化关系的进化树,为完成这一工作,已发展了多种软件包,如 PYLIP、MEGA 等(
9、4.稳定性检验:为了检验构建好的进化树的可靠性,需要进行统计可靠性检验,通常构建过程要随机地进行成百上千次,只有以大概率(70% 以上)出现的分支点才是可靠的。通用的方法使用 Bootstrap 算法,相应的软件已包括在构建系统进化树所用的软件包当中。 】当前的主要困难:是发现了基因的横向迁移(LGT)现象,即进化程度不同的物种间存在着遗传信息基因的传递,如果拿迁移的基因做进化分析就会出错。克服 LGT 的方法(可能的解决途径):1) 纵向思路:选择垂直进化而来的序列进行研究,即去除横向迁移的数据库,如COG 数据库;2) 横向思路:发展基于完整基因组构建进化树,即使用全基因组数据库进行基因组
10、水平上的对比;利用生物体的蛋白质组构建进化树。选取特征对比,不同长度的序列字符串进行对比后,对照其 genome 进行归一化;ORF 对比,将 all predicted ORF 采用 COG 的分类规则进行分类,再构建进化树4.什么是 SNP?为什么 SNP 的研究是重要的?SNP 研究有哪些优点?举出 23个 SNP 相关的网站。答:SNP 是指单核苷酸多态性,主要是指在基因组水平上由单个核苷酸的变异所引起的DNA 序列多态性,代表了基因组水平上遗传密码的变异,由于这种变异很多以单碱基突变的形式出现,因此称为单核苷酸多态性;它反映了不同个体间、正常与异常个体之间基因组上的差别,现在这个概念
11、有所扩大,不限于一个核苷酸的差异。重要性:因为 SNP 研究是基因组领域理论成果走向应用的关键步骤,是联系基因型和表现型之间关系的桥梁,是研究人类基因组计划走向应用的重要步骤。优点:(1)SNP 在基因组中分布相当广泛,使人们有机会发现与各种疾病相关的基因组突变;(2 )不直接导致疾病基因表达的 SNP,与某些疾病基因相邻,成为重要标记,有助于发现疾病基因(3 )从实验操作来看,通过 SNP 发现疾病相关基因突变,比通过家系发现更加容易。(4 )基础研究中非常重要,如对 Y 染色体 SNP 分析有重要成果。SNP 的特点:1.位点丰富2.具有代表性3.遗传稳定性4.易于进行自动化,规模化分析,
12、缩短了研究时间SNP 研究的意义:通过大批量、高通量的 SNP 的发现与鉴定,人类 SNPHaplotype 遗传图谱的构建,在连锁不平衡基础上的关联分析等,有望为人类致命基因的寻找和疾病的防治提供快速和有效的途径,一系列发现和检测 SNP 的方法,构建图谱的策略,及连锁不平衡和关联分析等技术,正在动植物研究领域中受到广泛的关注,毫无疑问将在分子和群体遗传、动植物育种和生物进化等研究领域中发挥越来越大的作用。SNP 相关的一些网站:1) SNP Consortiums database(http:/snp.cshl.org/index.html)2) NCBI SNP database 将这些
13、数据进行整理,去掉冗余,使每个 SNP 都是唯一的。此时的SNP 被称为 reference SNP 或refSNP。(( http:/www.ncbi.nlm.nih.gov/SNP/overview.html) 3) The Human Genic Bi-Allelic Sequences Database(HGBASE) 这一数据库收录了人基因组中所有已知的序列变化,包括:SNPs、序列的插入和缺失(Indels) 、简单重复序列等。(http:/hgbase.cgr.ki.se/)4) The Human Gene Mutation Database(HGMD)(http:/www.h
14、gmd.org/) 5) The Protein Mutant Database(PMD),蛋白突变数据库。收录了蛋白质特定位点的氨基酸突变信息,以及这些突变对蛋白质结构功能的影响。( http:/pmd.ddbj.nig.ac.jp/)6) The Allele Frequency Database(ALFRED):人类群体等位基因频率数据库,http:/alfred.med.yale.edu/alfred/index.asp5. 什么是系统生物学?系统生物学对生命科学概念上的发展?系统生物学对生物功能实现的理解有何本质变化?系统生物学的研究思路是什么?答:系统生物学是指在系统的层面上研究生
15、命活动。(研究一个生物系统中所有组成成分的构成,以及特定条件下组分间互作关系。)【 系统生物学就是自基因组研究以来,各个层次的所有资料和数据(包括基因组测序数据,功能基因组数据,蛋白质三维结构信息以及相互作用的数据等)的整合,以及这些整合数据为基础建立数学模型,再以这些模型模拟仿真研究生命活动的影响之后生命活动的反应以及变化 】包含三个相互衔接的组成(三部曲):整合数据,即整合所有各个层次(DNA 水平,RNA 水平,蛋白质水平,蛋白质相互作用水平)的信息数据;系统建模,即用这些信息构建描绘生命活动的数学模型;预测未知,即用这个模型预测生命未来的发展及外界干扰后系统的变异(生命活动及外界因素变
16、化对其产生的影响)。学术概念上的发展主要有:传统生物学是从基因组序列到结构,再到功能,而它从各个层次的相互作用到网络,再到功能。与以往不同的是,系统生物学一开始就考虑元件之间的相互作用,把整个生命活动作为网络,考虑其相互作用。1) 研究思路的变化:传统的分子生物学研究步骤一般为:DNA 序列蛋白结构蛋白功能(一维),而系统生物学是在二维的角度研究生命科学,即:相互作用网络功能,是由一组基因产生并相互作用共同实现的。2) 看待生命活动本质的变化:因为没有一个生命活动是靠一个基因完成的,生命活动是一组基因相互作用实现的,这种相互作用形成一个网络,既包括每个单元的结构,又包括单元与单元之间的相互作用
17、。因此,系统生物学不仅考虑每个基因的活动,还描述了基因间的相互作用并导致了网络的产生。(系统生物学与传统生物学看待生命活动有着本质的不同: 系统生物学认为生命活动是由一组基因及其相互作用来实现其过程的,这种相互作用形成了一个网络,既包括每个单元的结构,又包括单元与单元之间的相互作用,因此在考虑结构的过程中考虑其结构间的相互作用,一组一组地研究。而传统的分子生物学考虑的只有结构,是一个一个地去研究。)其对生物功能实现的理解发生了本质性变化:它不仅考虑单个分子而且考虑其间相互作用,把整个生命活动作为一个相互作用的网络来研究其功能,基因组只是网络中的一部分,只有通过相互作用的网络才能体现功能;通过系
18、统地整合生物过程不同阶段的分散数据,如基因组,转录组,蛋白组,代谢组,可以对复杂的生物过程,如折叠、信号传导途径、代谢途径更好地模拟,研究生物过程的动态变化;它不仅全息的了解复杂的生命系统中的所有成分以及他们之间的动态联系,还可以预测如果这个系统一旦受到了刺激和外界干扰,系统未来的行为是什么。系统生物学与传统生物学有什么不同:区别:传统生物学:序列结构功能,只考虑单个个体,单个 gene,单个蛋白质系统生物学:相互作用网络功能,除考虑单个个体,单个 gene,还考虑个体与个体之间的相互作用,把整个生命活动作为一个网络来考查它们的相互作用。(传统分子生物学是从基因组中发现特殊序列,即基因,然后找
19、到基因编码的蛋白,再通过测知其结构,而知其功能。而系统生物学研究是从各个层次的相互作用到网络,再到功能。系统生物学不仅考虑单个分子,而且考虑其间相互作用,认为生命活动由大量相互作用的结构单元组成,这些结构单元形成网络。基因组只是网络中的一部分,只有通过相互作用的网络才能体现功能。它不仅全息的了解复杂的生命系统中的所有成分以及他们之间的动态联系,还可以预测如果这个系统一旦受到了刺激和外界干扰,系统未来的行为是什么。)系统生物学与分子生物学有什么不同:区别:分子生物学:序列结构功能,只考虑单个 gene,单个蛋白质系统生物学:是研究生物系统组成成分的构成与相互关系的结构、动态与发生,以系统论和实验
20、、计算方法整合研究为特征的生物学。系统生物学不同于以往仅仅关心个别的基因和蛋白质的分子生物学,在于研究细胞信号传导和基因调控网路、生物系统组成之间相互关系的结构和系统功能的涌现。系统生物学的研究思路(研究流程):1.针对选定生物系统进行实验设计,了解系统所有组成成分:基因,RNA,蛋白,膜脂等2.通过系统行为动力学分析,总结系统设计和控制规律3.通过总结规律来提出新的实验设计,验证系统模拟的正确性【分子生物学与系统生物学的区别与联系? 答:二者的区别和联系主要从宏观和微观上讲。分子生物学的研究采用典型的还原论方法,研究对象主要是分子水平上的,即生物系统中的大分子、信号分子的结构、生化性质以及功
21、能,基因表达过程中的调控,以及 DNA 重组。分子生物学只研究系统的组成元素,最后给出系统的组成元素清单,它是系统生物学的基础,但它的研究结果只能解释生物系统的微观或局部现象,无法说明系统整体所具有的功能从何而来。而系统生物学作为一个整体,表现出完善的整体行为,而组成系统的细胞、基因、蛋白质等只能作为系统的一个构件、一个元素、通常情况下它无法表现出“系统”行为。系统生物学与分子生物学研究对象不同,系统生物学研究的是系统整体,研究由系统元素形成有功能的整体所依赖的组织方式和潜藏规则,它同时研究系统的不同层次,以及他们之间的相互作用关系,并将这些整合起来深刻挖掘系统整体的功能形成机制。系统生物学虽
22、然在研究对象上与分子生物学不同,但他们之间并不是完全不相关的,系统生物学的研究离不开分子生物学研究所给出的大量资料和数据,正是依赖这些,系统生物学才有了建模的基础。同时分子生物学的研究结果只有通过系统生物学进行整合才能从理论上对系统的宏观性质达到定性定量的理解,反过来,系统生物学的研究成果也可以用来指导分子生物学的实验设计。因此二者之间其实是相互补充的,只有结合起来,才能充分认识生命现象。】6.( 1)什么是非编码序列,非编码 RNA,非编码基因?(2)以人的基因组为例回答:在基因组中有多少非编码序列,有多少存在转录本,举 23 个非编码核酸的生物学功能?答:(1)非编码序列是基因组中不编码蛋
23、白质和多肽的序列;(基因组中不归属于基因调控元件,稳定元件之外的,也无明确生物学功能意义的基因序列统称为非编码序列,即不编码蛋白质同时也无明确生物学功能的序列)非编码 RNA 是指来自基因组的非编码的转录元件,即基因组中非编码序列的转录产物/ 转录本;非编码基因指那些具有明确生物学功能的非编码 RNA 在基因组上非编码序列上的位置,即功能性的非编码 RNA 对应基因组上的位置称为非编码基因;(2)人类基因组中 9798%的序列是非编码序列,有 70%80%存在转录本,非编码核酸的生物学功能:1) Xist:X-inactivation(X 染色体失活)是哺乳动物的一种剂量补偿机制,其中一半拷贝
24、转录被抑制从而失活,抑制转录是通过一个 2kb 的非编码 RNA(Xist RNA)实现的,xist RNA 装配在失活 X 染色体的外侧,引起结构改变导致失活;2) Small RNA and RNAi: RNAi 是由 RNA(siRNA 、microRNA)导致的转录后基因沉默现象,如由双链小 RNA 引起的干扰和转录后基因沉默现象,在植物病毒抗性和线虫中的转座子沉默;一些小核 RNA 调控基因转录。(单链易降解,但发现细胞中存在另一种 pathway,双链小 RNA 进入细胞后结合组蛋白形成复合体,该复合体和识别并降解 target)3) piRNA(具有大量转录本,功能不详)和 Pr
25、ions(生物复杂度到一定程度后会出现发病情况,可能和非编码 RNA 有关)等。7.什么是基因组中的非编码区?请以人类基因组为例,说明:(1)非编码区所占的比例?(2)按在基因组中的位置(组成)(功能)区分,非编码序列有哪些组分?它们所占比例如何?(3)按序列编码特征区分,非编码序列有哪些组分?它们所占比例如何?(4)请说明非编码区研究的重要性(可以举出一、两个典型非编码序列作为例子)答:基因组中不能编码蛋白质的区段叫做非编码区。非编码区位于编码区前后,同属于一个基因,控制基因的表达和强弱 。(1)人类非编码区占 9798%(2)按照在基因组中的位置(组成)来分,各个组分占基因组的份额:编码基
26、因(编码蛋白质和 tRNA、rRNA):1.52% ;Intron(广义):25% ;端粒、中心粒等特定位置:12% ;基因间序列:6070% ;按照在基因组中的功能区分,各个组分占基因组的份额:功能蛋白质基因 1.7%,功能 RNA 基因 0.5%,总共大约 13% ;内含子:24% ;Satellite DNA(主要分布在中心粒和端粒 ): 12% ;基因间序列(Intergene DNA):6070% ;(3)按照序列特征区分,各个组分占基因组的份额:编码区(包括编码蛋白质和 tRNA 和 rRNA 的基因)占总基因组的 2% ;非编码区占到 98%:其中:简单重复序列:12% ;散在重
27、复序列:45% ;假基因:1% ;非编码非重复序列:3540% ;(4)举例:非编码基因:1.SINE 作为调节源,调节基因重组、交换,丰富多样性,获得新功能;2. 鸡溶菌酶基因中,位于编码区上游的 CRI 元件起着转录沉默子的作用;3.ncDNA 产物有重要生物学功能,如 tmRNA 介导错误翻译蛋白的降解 RNAi导致基因沉默非编码基因产物的功能:smallRNA 是 ncDNA 产物,是机体固有的,例如:microRNA,SiRNA 小 RNA 对染色质的形状有关,也可直接关闭或删除部分DNA。 NcRNA 起着非常重要的生物学功能,如影响发育过程,调节转录、影响染色体复制、对 RNA
28、加工修饰、影响 mRNA 稳定性进而影响翻译、甚至影响蛋白降解转运;Xist 介导 X 染色体失活是通其编码的一个大的剪接过的多聚 A 非编码产物进行的。(长链非编码 RNA(lncRNA)是一类转录本长度超过 200nt 的 RNA 分子,它们并不编码蛋白,而是以 RNA 的形式在多种层面上(表观遗传调控、转录调控以及转录后调控等)调控基因的表达水平。 lncRNA 起初被认为是基因组转录的“噪音”,是 RNA 聚合酶 II 转录的副产物,不具有生物学功能。然而,近年来的研究表明,lncRNA 参与了 X 染色体沉默,基因组印记以及染色质修饰,转录激活,转录干扰,核内运输等多种重要的调控过程
29、,lncRNA 的这些调控作用也开始引起人们广泛的关注。哺乳动物基因组序列中约 4%9%的序列产生的转录本是lncRNA(相应的蛋白编码 RNA 的比例是 1%),虽然近年来关于 lncRNA 的研究进展迅猛,但是绝大部分的 lncRNA 的功能仍然是不清楚的。)( 已有的研究结果表明,在高等生物中,小分子非编码 RNA 在干细胞干性维持、胚胎发育、细胞分化、凋亡、代谢、信号传导、感染以及免疫应答等几乎所有重要生命活动中发挥关键的调控作用,提示生物体内可能存在着由 RNA 介导的遗传信息表达调控网络。)8.精准医学的重大意义是什么?实现精准医学的重要基础是什么?精准医学的重大意义;精准医学有可
30、能导致医疗体系本质上的转变,把目前的医疗体系由诊断治疗过渡到健康保障,使得健康体系的关口前移,有可能产生新兴产业。健康人可以通过组学等一系列研究,对现在的健康作以评估。在健康检查的基础上,对未来可能导致疾病的部分进行干预,使得能够延缓疾病的发生,或者排除某些疾病的发生,使得健康得以保障。实现精准医学的重要基础:1.必须获取分子水平上的数据信息,并挖掘其内涵,在挖掘组学数据时,一定要使用大数据分析技术,因此是大数据与组学的交汇。组学包括基因组,转录组,蛋白质组,代谢组;大数据包括人群和队列2.建立分子水平上的知识与宏观疾病表型的联系,即基因型和表型的关联,搭建分子水平信息和疾病间的桥梁,在搭建桥
31、梁时,生物信息学,生物网络,系统生物学的知识是其核心知识。3.在此基础上,融合临床检验,影像学等指标,使得医学做得更加精准。【定义:精准医学是以个体化医疗为基础、随着基因组测序技术快速进步以及生物信息与大数据科学的交叉应用而发展起来的新型医学概念与医疗模式。本质上:是通过基因组、蛋白质组等组学技术和医学前沿技术,对于大样本人群与特定疾病类型进行生物标志物的分析与鉴定、验证与应用,从而精确寻找到疾病的原因和治疗的靶点,并对一种疾病不同状态和过程进行精确亚分类,最终实现对于疾病和特定患者进行个性化精准治疗的目的,提高疾病诊治与预防的效益。精准医学是因人因病而异的、更加精确的个体化医疗,其进步之处是
32、将人们对疾病机制的认识与生物大数据和信息科学相交叉,精确进行疾病分类及诊断,为疾病患者提供更具针对性和有效性的防疗措施,最终目的是更好地为患者服务。与个体化医疗相比,精准医疗更重视“病”的深度特征和“药”的高度精准性;是在对人、病、药深度认识基础上,形成的高水平医疗技术。精准医学实现了从诊断治疗到健康保障这一本质性转变。精准医学包括精准诊断和精准治疗,而“迈向精准医学”需要构造的生物医学知识网络是建立在系统生物学的基础之上。实施精准医学计划的战略意义总共有 4 点:提高疾病诊治水平,惠及民生与国民健康;推动医学科技前沿发展,增强国际竞争力;发展医药生物技术,促进医疗体制改革; 形成经济新增长点
33、,带动大健康产业发展。】【有可能将基因组变异作为疾病诊断,精准医学导致医疗体系本质的转变,把目前阶段治疗过渡到健康保障,使健康体系关口前移,在健康筛查基础上,排除疾病发生。就是评估-检查-干预的过程。基础:1 获取分子水平上数据信息,挖掘信息内容,发展大数据新算法,理论技术如组学的信息。2 建立分子水平知识宏观疾病表型关联,搭建分子水平信息与疾病的桥梁。问题:样本量少,有效治疗事件频率低,疾病相关复杂网络构建分析的困难。】陈小伟老师部分:1.芯片间标准化的方法:排序:每列由小到大排序,而可以得到每一组基因表达量的真实值求排的平均值作为标准值基本方法:芯片间标准化的目的是基于 Gene1Gene
34、5 五个基因表达量理论的和应该保持恒定,即 S1S3 三列每一列的和是相等的。但实际测定过程中不可能完全相等,因此将这种不等归结于每一组芯片自身的差异而进行芯片间标准化,基本步骤为通过排序取平均重新排序的方法消除芯片间误差,从而可以得到每一组基因表达量的真实值。(老师给的这组芯片基因完全相同的情况下 S3 一列数据明显偏高,通过这种标准化实现了芯片间差异的消除)。【Quantile 归一化过程:首先假设不同芯片整体分布一致,归一化后芯片的分布一样。下图四个部分代表四步,行代表基因,列代表样本,图一对每个列的表达值排序,图二计算每行的平均值,图三用每行计算的平均值代替该行的原值,图四将排序后的行
35、恢复到未排序前的位置。】 2. FDR 控制假阳性的方法BenjaminiHochberg procedure基本方法:对于 m 个独立的样本,其 p-value 记为 pi,i=1,2,3m;(1)对所有的 p-value 进行从小到大排序 p(1)p(2) p(m);(2)对于一个给定的 (此时的 即为统计里的显著水平,范围 01,通常取 0.05),找到最大的 k 值,满足 ;(3)拒绝从 p(1)p(k)的无效假设 H0(即表示 p(1)p(k)表达量存在显著差异)。重排:按颜色重排计算方法 1(=0.05):P(4)=0.030.05*5/6=0.041;k=4. 即 G2, G6,
36、 G5, G4 差异表达,FDR0.05计算方法 2(q-value 法):根据 可以推出 因此直接计算并与 进行对比即可:由于 G3 的 q-value 大于 0.05,因此 G2, G6, G5, G4 差异表达。【FDR 过程,如何控制 FDR:首先,FDR 过程是为了控制假阳性率的过程。假阳性指样本本质为假但判定为真。比如在找到一组差异表达的基因之后,我们要考虑这个差异是否够显著,即假阳性率是否足够低。BH 过程是 FDR 矫正的一种,首先对检验得到的 P 值进行排序,然后从 1 开始增加找K 的值,使 K 满足 ,其中 m 为个数,a 一般取 0.05 或 0.1。所有满足 P 值的基因认为表达有显著差异且假阳性不超过 a。】3. 转录本表达量的表示方法( RPKM:Reads Per Kilobase of transcript per Million mapped reads):
Copyright © 2018-2021 Wenke99.com All rights reserved
工信部备案号:浙ICP备20026746号-2
公安局备案号:浙公网安备33038302330469号
本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。