从生物实验角度理解芯片生物信息分析

  • 格式:pdf
  • 大小:386.74 KB
  • 文档页数:13

生物芯片技术中的生物信息浅谈

博奥生物公司 张亮 2005年7月 生物芯片这门强大的生物分析工具,以它低消耗、高通量的巨大优势正在生

物医学的各个领域得到越来越多的应用。生物芯片从在它的反应原理来说,可以

分为两大类,具体分类可参见图1-1。

图1-1. 生物芯片的主要类别

结合现有的生物芯片主要类别,这里我们给生物芯片下一个定义。生物芯片:

指能储藏大量生物信息或快速并行处理多个生物样品的微器件,它的加工运用了

微电子工业中十分成熟的光学光刻技术和微机电系统加工中所采用的各种方法,

所处理的对象是生物样品,故称之为生物芯片。

正是由于生物芯片的高通量优势,一次实验能产生海量的数据,因此如何充

挖掘生物芯片产生的数据正在成为生物芯片应用的瓶颈。鉴于目前对生物芯片的

生物信息分析多集中于分析DNA芯片产生的数据,这里我们以DNA芯片为例

来谈谈生物芯片技术中的生物信息分析。我们将顺着DNA芯片的技术过程,来

谈谈生物信息在DNA芯片技术中的应用 一、在DNA芯片的制备过程中生物信息的应用

DNA芯片分为cDNA芯片和寡聚核苷酸芯片(Oligo芯片)。cDNA芯片是

通过PCR扩增的方法把靶序列扩增出来,然后把PCR产物用芯片点样仪点在经

过化学修饰的基片上(通常是75×25mm标准规格的玻璃片);cDNA芯片的一

个重要特点是:并不一定需要知道靶序列的序列信息,例如未经测序的cDNA

克隆都可以制备成cDNA芯片。Oligo芯片是根据靶基因的核酸序列而人工设计

并化学合成探针,目前Oligo芯片包括美国Affymetrix公司的基因芯片平台,它

的技术是在硅基质载体上原位合成25 mer长度的寡聚核苷酸,并且每个基因有

多条寡聚核苷酸对应,这些寡聚核苷酸和检测的基因序列之间是完全匹配和非完

全匹配的关系;另外一种Oligo芯片是针对每个基因先人工化学合成比较长的寡

聚核苷酸,一般是50-90个碱基,然后把这些合成好的寡聚核苷酸配制成一定

浓度的溶液,用芯片点样仪点到经过化学修饰的玻璃片上,我们把这类芯片称之

为长Oligo芯片。在一般实验室,拥有了芯片点样仪和芯片激光扫描仪后都可以

自行制备cDNA芯片和长Oligo芯片。

制备cDNA芯片从价格上要相对比Oligo芯片便宜,目前在制备一些微生物

全基因组的基因芯片时,常采取cDNA芯片的制备方法。比如有一种微生物的

全基因组序列,每个基因编码蛋白的开放阅读框(open reading frame, ORF)序

列都已经给出,格式如下:

>XC4319

ATGACCACCCACGACGACGAACGCAAGGCGCTGCGCCGGTTGTTGCGGCAGATCGAG

CGCCC……

>XC4320

ATGAATACTCCTGTCCCCGTCGTCCGACTCAAGAATGCGTGGCGCTCCAGCCACCCGT

GGAT…… 在生物实验工作者的日常工作中,要单独设计每个基因的PCR扩增引物是

很简单的,可以利用Primer3.0的引物设计程序 (http://frodo.wi.mit.edu/primer3 /primer3_code.html) 或者商用的Primer Premier 5.0等引物设计软件;但对于DNA

芯片技术,成千上万个PCR产物会在同一个反应条件下进行杂交反应,因此设计某个微生物全部基因对应的PCR引物,就需要生物信息工作者参与。引物设

计的核心仍可使用Primer3.0的引物设计程序,但需要使用Perl或Java编写特异

性分析程序,对Primer3.0设计的引物进行特异性分析,把生物实验工作者要求

的参数,例如PCR产物的长度、引物的退火温度、对基因序列同源性高的基因

如何扩增差异最大的区段等条件,输入到生物信息者编写的批处理程序中。生物

芯片北京国家工程研究中心的生物信息工作者可以在48小时内设计出5000个基

因对应的PCR扩增引物,给出的结果格式如下:

Name Sequence (5’-3’) Tm Product length

XC4319-F ACGACGACGAACGCAAGGCGC 67.2 244

XC4319-R GGCCACTGCCGTGCCCACAC 68.7

XC4320-F GACTCAAGAATGCGTGGCGC 62.5 585

XC4320-R GCTGGTCGGCAAAGAAACCG 62.5

……

若是制备70mer 的长Oligo芯片,那么在设计70mer的寡聚核苷酸时需要遵

循的原理包括

1. 所有的Oligo探针的Tm值在78℃±5℃;

2. 所有的Oligo探针尽量在靠近基因序列3'端的1000个碱基内,以便以PolyA

为引物做反转录时,探针的位置尽可能位于反转录产生的cDNA序列内;

3. 在Oligo探针中,不能有连续7个相同的碱基;

4. 所有的探针不能有长于9个碱基的自身发卡结构;

5. 每个设计出来的Oligo探针都要和所有的人的基因进行序列比对,不能超过

70%的同源性; 并且和人的其他的基因不能存在连续20个连续完全匹配的

碱基。

等等。当然,设计长Oligo探针过程也是通过生物信息工作者编写的程序由

计算机自动完成。 二、在杂交图片转化为数据信号过程中生物信息的应用

DNA芯片制备好后,可以把两个RNA样品经过系列的酶学反应分别标记上

具有不同吸收和激发波长的荧光素,例如cy5和cy3荧光素,然后用双通道激光

扫描仪分别获得cy5和cy3通道的杂交图片。在扫描的时候,要注意调节cy5和

cy3通道的扫描参数,主要是激光光源的强度和光电倍增管放大倍数,使芯片上

绝大多数的点cy5和cy3的扫描强度大致相同。杂交图片上基因的杂交原始信号

是以216灰阶Tif图片格式进行保存的。可以用基因芯片杂交图片分析软件,例

如美国Axon公司GenPix Pro软件,美国PE公司的QuantArray软件,美国

BioDiscovery公司的ImaGene软件,北京博奥生物芯片公司的SpotData软件都

可把Tif的灰阶图像信息转化为数据信号(从0-65536)。例如某个X基因在cy5

标记的处理样品通道中得到的杂交信号是55000,在cy3标记的对照样品通道得

到的杂交信号是2500。这样可以初步知道此X基因在处理样品中的表达丰度是

对照样品的22倍。但实际上,这个比值是原始的比值,还需要经过归一化处理。

双通道基因芯片技术的归一化往往是针对信号的比值来进行的。芯片的原始数据

之所以要经过归一化处理是由DNA芯片实验中存在着系统误差,系统误差的来

源包括荧光染料的差异性,例如Cy5和Cy3有不同的量子产率,对光猝灭的敏

感性不同,在反转录过程中标记cDNA的效率也不一样,还包括试验者操作不

完全一致性,以及在扫描时不同荧光通道的硬件扫描设置不同等因素。目前在双

通道基因芯片技术中常用的归一化方法有:1)看家基因归一化;2)全部基因信

号值归一化;3)基于统计学方法LOWESS程序的归一化;4)利用外标基因进

行归一化。

1) 看家基因归一化是在基因芯片技术初期采用得最广泛的归一化方法,它是从

传统的Northern印迹方法继承而来,基本的思想是先确定在所比较的生物样

品中表达不发生变化的基因,例如3-磷酸甘油脱氢酶编码基因(GAPD),在

原始的数据中,GAPD的荧光信号值在处理样品/对照样品=1.4,那么得到的

X基因的处理样品/对照样品比值最初为22,归一化后为22/1.4=15.7倍。用

看家基因进行归一化的缺陷是:在一些实验条件下,看家基因表达也常会发

生变化,因此若选择的看家基因自身表达发生了变化,那么得到的芯片数据就是错误的;另外,由于看家基因的表达丰度相对比较高,因此杂交信号比

较强,那么由看家基因得到的校正系数对于那些丰度中等或更弱的基因来说

并不是等同的,导致归一化精确性差。

2) 全部基因信号值归一化方法;它是基于从整体上说,所比较的两个生物样品

基因表达总体上水平是一样的,因此在cy5通道中得到的基因所有杂交信号

之和应该和cy3杂交信号之和相等;例如在一次实验中cy5通道中所有的基

因信号值之和为90,000,000,在cy3通道中所有基因信号值之和为

60,000,000,那么cy3通道的信号都要乘以1.5的系数。此归一化的思想和看

家基因归一化方法一样,实际也是一种线性归一化方法;而在一个通道上的

基因杂交信号是从高到低分步的,这些杂交信号不同的基因归一化的系数应

该是不一样的,因此任何基因都乘以同一个系数会影响数据的精确性。

3) 基于统计学方法LOWESS程序的归一化是在2)的基础上发展而来,它的基

本数学模型是:设R-红色荧光信号强度,即Cy5荧光信号的强度;G-绿色荧

光信号强度,即Cy3荧光信号的强度;此方法认为:log2R/G-c(A) =

log2R/(k(A)G);引进参数M= log2R/G,A=1/2× log2RG,c(A) 不是一个常数,

而是一个与R,G相关的变量[1]。可以通过统计学软件Splus语言中的lowess

fit命令来作出 M vs. A的散点图。此方法归一化的基本思想是若没有系统误

差的影响,cy5通道中得到的基因所有杂交信号之和应该和cy3杂交信号之

和一致,即R=G,则M vs. A的散点图的分布应该趋向于M=0这根轴,参见

图1-2。对于不同信号强度的基因用lowess fit进行校正以后,可以设置不同

置信区间,找出此置信区间内表达真正发生变化的基因。一般当芯片上有少

于10%的基因表达发生变化时,也可以认为M vs. A的散点图的分布趋向于

M=0,用这种归一化方法是可行的。在用LOWESS归一化方法时,还可以

考虑玻片不同杂交区域的影响,即空间位置的影响。目前芯片杂交一般在盖

玻片下进行,液体的流动性不好,常造成不同杂交区域信号强度不一,另外,

一个芯片上的DNA样品在点样时往往是由多根针点样,每一根针点的样就

组成空间位置上相对靠近的一个小亚阵。不同针的点样量和点的大小会存在

差异,因此需要先进行每个小亚阵中基因的归一化,即用上述方法认为每个

小亚阵上的基因(往往超过400个)的Cy5和Cy3的荧光信号存在相同的系