1.什么是生物信息学,如何理解其含义?答:生物信息学有三个方面的含义:1)生物信息学是一个学科领域,包含着基因组信息的获取、处理、存储、分配、分析和解释的所有方面。2)生物信息学是把基因组DNA序列信息分析作为源头,破译隐藏在DNA序列中的遗传语言,特别是非编码区的实质;同时在发现了新基因信息之后进行蛋白质空间结构模拟和预测;其本质是识别基因信号。3)生物信息学的研究目标是揭示“基因组信息结构的复杂性及遗传语言的根本规律”。它是当今自然科学和技术科学领域中“基因组、“信息结构”和“复杂性”这三个重大科学问题的有机结合。怎样理解生物信息学:生物信息学是把基因组DNA序列信息分析作为源头,找到基因组序列中代表蛋白质和RNA基因的编码区;同时阐明基因组中大量存在的非编码区的信息实质,破译隐藏在DNA序列中的遗传语言规律:在此基础上,归纳、整理与基因组遗传信息释放及其调控相关的转录谱和蛋白谱数据,从而认识代谢、发育、分化、进化的规律。其还利用基因组中编码区信息进行蛋白空间结构模拟和蛋白功能预测,并将此类信息与生物体和生命过程中的生理生化信息结合,阐明其分子机制,最终进行蛋白、核酸分子设计、药物设计、个体化医疗保健设计。2.如何利用数据库信息发现新基因,基本原理?答:利用数据库资源发现新基因,根据数据源不同,可分2种不同的查找方式:1)从大规模基因组测序得到的数据出发,经过基因识别发现新基因:利用大规模拼接好的基因组,使用不同数据方法,进行标识查找,并将找到的可能的新基因同数据库中已有的基因对比,从而确定是否为新基因。可分为:①基于信号,如剪切位点、序列中的启动子与终止子等。②基于组分,即基因家族、特殊序列间比较,Complexityanalysis,NeuralNetwork2)利用EST数据库发现新基因和新SNPs:数据来源于大量的序列小片段,EST较短,故关键在正确拼接。方法有基因组序列比对、拼接、组装法等。经常采用SiClone策略。其主要步骤有:构建数据库;将序列纯化格式标准化;从种子库中取序列和大库序列比对;延长种子序列,至不能再延长;放入contig库①构建若干数据库:总的纯化的EST数据库,种子数据库,载体数据库,杂质、引物数据库,蛋白数据库,cDNA数据库;②用所用种子数据库和杂质、引物数据库及载体数据库比对,去除杂质;③用种子和纯化的EST数据库比对④用经过一次比对得到的长的片段和蛋白数据库、cDNA数据库比较,判断是否为已有序列,再利用该大片段与纯化的EST数据库比对,重复以上步骤,直到序列不能再延伸;⑤判断是否为全长cDNA序列。(利用EST数据库:原理:当测序获得一条EST序列时,它来自哪一个基因的哪个区域是未知的(随机的),所以属于同一个基因的不同EST序列之间常有交叠的区域。根据这种“交叠”现象,就能找出属于同一个基因的所有EST序列,进而将它们拼接成和完整基因相对应的全长cDNA序列。而到目前为止,公共EST数据库(dbEST)中已经收集到约800万条的人的EST序列。估计这些序列已覆盖了人类全部基因的95%以上,平均起来每个基因有10倍以上的覆盖率。)3.用蛋白或核酸序列数据库研究生物演化的主要步骤是什么?当前的困难是什么,如何克服?答:构建系统进化树,其主要步骤如下:1)序列相似性比较。就是将待研究序列与DNA或蛋白质序列库进行比较,用于确定该序列的生物属性,也就是找出与此序列相似的已知序列是什么。完成这一工作只需要使用两两序列比较算法。常用的程序包有BLAST、FASTA等;2)序列同源性分析。是将待研究序列加入到一组与之同源,但来自不同物种的序列中进行多序列同时比较,以确定该序列与其它序列间的同源性大小。这是理论分析方法中最关键的一步。完成这一工作必须使用多序列比较算法。常用的程序包有CLUSTAL等;3)构建系统进化树。根据序列同源性分析的结果,重建反映物种间进化关系的进化树。为完成这一工作已发展了多种软件包,如PYLIP、MEGA等;4)稳定性检验。为了检验构建好的进化树的可靠性,需要进行统计可靠性检验,通常构建过程要随机地进行成百上千次,只有以大概率(70%以上)出现的分支点才是可靠的。通用的方法使用Bootstrap算法。当前的主要困难是发现了基因的横向迁移(LGT)现象,即进化程度不同的物种间存在着遗传信息基因的传递,如果拿迁移的基因做进化分析就会出错。克服LGT的方法:1)选择垂直进化而来的序列进行研究,即去除横向迁移的数据库,如COG数据库;2)使用全基因组数据库进行基因组水平上的对比;利用生物体的蛋白质组构建进化树。选取特征对比,不同长度的序列字符串进行对比后,对照其genome进行归一化;ORF对比,将allpredictedORF采用COG的分类规则进行分类,再构建进化树4.什么是SNP,为什么SNP的研究是重要的,举出2~3个SNP相关的网站。答:SNP是指单核苷酸多态性,代表了基因组水平上遗传密码的变异,由于这种变异很多以单碱基突变的形式出现,因此称为单核苷酸多态性;因为SNP研究是基因组领域理论成果走向应用的关键步骤,是联系基因型和表现型之间关系的桥梁,是研究人类基因组计划走向应用的重要步骤。SNP相关的一些网站:1)SNPConsortium'sdatabase()2)NCBISNPdatabase将这些数据进行整理,去掉冗余,使每个SNP都是唯一的。此时的SNP被称为referenceSNP或refSNP。(()3)TheHumanGenicBi-AllelicSequencesDatabase(HGBASE)这一数据库收录了人基因组中所有已知的序列变化,包括:SNPs、序列的插入和缺失(Indels)、简单重复序列等。()4)TheHumanGeneMutationDatabase(HGMD)()5)TheProteinMutantDatabase(PMD),蛋白突变数据库。收录了蛋白质特定位点的氨基酸突变信息,以及这些突变对蛋白质结构功能的影响。()6)TheAlleleFrequencyDatabase(ALFRED):人类群体等位基因频率数据库,什么是系统生物学?系统生物学对生命科学概念上的发展?答:系统生物学是指在系统的层面上研究生命活动。(研究一个生物系统中所有组成成分的构成,以及特定条件下组分间互作关系。)包含三个相互衔接的组成:整合数据,即整合所有各个层次(DNA水平,RNA水平,蛋白质水平,蛋白质相互作用水平)的信息数据;系统建模,即用这些信息构建描绘生命活动的数学模型;预测未知,即用这个模型预测生命未来的发展及外界干扰后系统的变异。概念上的发展主要有:1)研究思路的变化:传统的分子生物学研究步骤一般为:DNA序列→蛋白结构→蛋白功能(一维),而系统生物学是在二维的角度研究生命科学,即:相互作用→网络→功能,是由一组基因产生并相互作用共同实现的。2)看待生命活动本质的变化:因为没有一个生命活动是靠一个基因完成的,生命活动是一组基因相互作用实现的,这种相互作用形成一个网络,既包括每个单元的结构,又包括单元与单元之间的相互作用。因此,系统生物学不仅考虑每个基因的活动,还描述了基因间的相互作用并导致了网络的产生。6.(1)什么是非编码序列,非编码RNA,非编码基因?(2)以人的基因组为例回答:在基因组中有多少非编码序列,有多少存在转录本,举2~3个非编码核酸的生物学功能?答:(1)非编码序列是基因组中不编码蛋白质和多肽的序列;非编码RNA是基因组中非编码序列的转录产物/转录本;功能性的非编码RNA对应基因组上的位置称为非编码基因;(2)人类基因组中97~98%的序列是非编码序列,有70%~80%存在转录本,非编码核酸的生物学功能:1)Xist:X-inactivation(X染色体失活)是哺乳动物的一种剂量补偿机制,其中一半拷贝转录被抑制从而失活,抑制转录是通过一个2kb的非编码RNA(XistRNA)实现的,xistRNA装配在失活X染色体的外侧,引起结构改变导致失活;2)SmallRNAandRNAi:RNAi是由RNA(siRNA、microRNA)导致的转录后基因沉默现象,如由双链小RNA引起的干扰和转录后基因沉默现象,在植物病毒抗性和线虫中的转座子沉默;一些小核RNA调控基因转录。(单链易降解,但发现细胞中存在另一种pathway,双链小RNA进入细胞后结合组蛋白形成复合体,该复合体和识别并降解target)3)piRNA(具有大量转录本,功能不详)和Prions(生物复杂度到一定程度后会出现发病情况,可能和非编码RNA有关)等。1.芯片间标准化的方法:基本方法:芯片间标准化的目的是基于Gene1~Gene5五个基因表达量理论的和应该保持恒定,即S1~S3三列每一列的和是相等的。但实际测定过程中不可能完全相等,因此将这种不等归结于每一组芯片自身的差异而进行芯片间标准化,基本步骤为通过排序取平均重新排序的方法消除芯片间误差,从而可以得到每一组基因表达量的真实值。(老师给的这组芯片基因完全相同的情况下S3一列数据明显偏高,通过这种标准化实现了芯片间差异的消除)。2.FDR控制假阳性的方法——Benjamini–Hochbergprocedure基本方法:对于m个独立的样本,其p-value记为pi,i=1,2,3…m;(1)对所有的p-value进行从小到大排序p(1)≤p(2)≤p(m);(2)对于一个给定的α(此时的α即为统计里的显著水平,范围0~1,通常取0.05),找到最大的k值,满足p(k)≤α∙km;(3)拒绝从p(1)~p(k)的无效假设H0(即表示p(1)~p(k)表达量存在显著差异)。计算方法1(α=0.05):P(4)=0.030.05*4/6=0.033;P(5)=0.0450.05*5/6=0.041;k=4.即G2,G6,G5,G4差异表达,FDR0.05计算方法2(q-value法):根据p(k)≤α∙km可以推出α≤p(k)∙mk因此直接计算并与α进行对比即可:由于G3的q-value大于0.05,因此G2,G6,G5,G4差异表达。3.转录本表达量的表示方法(RPKM:ReadsPerKilobaseoftranscriptperMillionmappedreads):(1)RPKM的作用:RNA-seq是透过次世代定序的技术来侦测基因表现量的方法,在衡量基因表现量时,若是单纯以map到的read数来计算基因的表现量,在统计上是一件相当不合理的事,因为在随机抽样的情况下,序列较长的基因被抽到的机率本来就会比序列短的基因较高,如此一来,序列长的基因永远会被认为表现量较高,而错估基因真正的表现量,所以AliMortazavi等人在2008年提出以RPKM在估计基因的表现量假设一个物种的基因组上只有两个基因,基因G1的外显子长8Kb,基因G2的外显子长2Kb。对该物种的一个样本做RNA-seq,共得到23millions的read,其中能够比对到G1的read有16million个,能够比对到G2的有4million个.计算G1和G2的RPKM。Totalmappedreads=16million+4million=20millionG1:totalexonreads=16,000,000exonlength=8kbRPKM=16,000,000/(20*8)=100,000G2:totalexonreads=4,000,000exonlength=2kbRPKM=4,000,000/(20*2)=100,000(2)FPKM与