[转载]第二代测序技术(NGS)漫谈
NGS之基础篇
2001年,美、英、法、德、日、中六国合作,历时十年,耗资数十亿美元的人类基因组计划(Human Genome Project,HGP)宣告完成。转眼又是十年过去,在此期间,各国科学家仍在为解读基因的密码而不懈努力,这其中最大的突破,就是第二代测序技术的推 出。HGP的顺利完成证明了我们有能力对自身的遗传信息进行研究,然而,高昂的成本、漫长的时间、巨大的人力需求,无不限制着对遗传密码的进一步认识。从 HGP开始的第一天期,科学家们就在寻求更好的方法来对基因组进行研究,“鸟枪法”就是其中之一。2006年,美国X大奖基金会(www.xprize.org) 设立了奖金高达1000万美元的基因组Archon X大奖,旨在奖励第一个在10天内以低于100万美元的成本完成100个人类基因组测序的民间团队。而罗氏(Roche)、应用生物系统(Applied Biosystems,ABI)、Illumina三家公司先后推出了各自的第二代高通量测序平台,成为NGS领域的领头羊。
2005年底,454公司推出第一个基于焦磷酸测序原理的高通量基因组测序系统——Genome Sequencer 20 System,这是核酸测序技术发展史上里程碑式的事件。随后,罗氏公司以1.55亿美元收购了454公司,并在2006年推出了更新的GS FLX测序系统,该系统可在10小时的运行中获得100万条读长(reads),4~6亿个碱基信息(base pair),且准确率达到99%以上。2008年,GS FLX系统再次升级,通量提高了5倍,读长和准确率也有所增加。虽然454 GS测序平台也许不是市场占有率最高的测序仪,但截至2011年3月,利用该系统进行研究的论文已发表超过1000余篇,而它在读长上的优势明显胜于另两 套系统,因此在从头测序(de novo)和宏基因组测序(meta genome)方面有着不可替代的地位。
2006年,Solexa公司也推出了自己的NGS系统——Genome Analyzer,简称GA。这套基于DNA簇(DNA cluster)、桥式PCR(Bridge PCR)和可逆阻断(Reversible terminator)等核心技术的系统具有高通量、低错误率、低成本、应用范围广等优点。2007年,Illumina公司以6亿美元的高价收购了 Solexa,使GA得以商品化。GA最早期的版本一次运行可获得1Gb的数据,因此也有1Gb Analyzer的含义,而最新的Hiseq2000平台则能够在10天的运行中获得300Gb以上的数据,读取的碱基长度达到150bp左右。更有消息 称,Illumina已完成了600Gb的运行测试并在部分客户中开展了前期体验,Tb(1000Gb)级的测试Run也将于年内进行。据不完全统 计,Illumina公司已售出超过600台/套GA IIx和Hiseq2000平台,2010年仅深圳华大基因研究院一家就购买了128台Hiseq2000,一举成为全球最大的基因组测序与分析中 心,Illumina公司在测序领域的影响力由此可见一斑。
在Sanger测序时代,美国应用生物系统公司(ABI)一直是该行业的龙头老大,其垄断地位无人能撼,从早期的377到全自动化的3730xl,ABI 的测序仪被广泛应用在基因组学研究的各个方面。然而在第二代测序技术迅猛发展之初,ABI起步较晚,显得有些漫不经心。直到2005年454公司推出GS 平台,ABI的领先地位受到威胁,这才开始发力,迅速收购了研发NGS的一家小公司Agencourt,并于2007年推出了它的SOLiD测序平台。此 后SOLiD不断升级,目前已到SOLiD 5版本(SOLiD 5500xl)。SOLiD的全称是Sequencing by Oligo Ligation Detection,即寡聚物连接检测测序,其基本原理是通过荧光标记的8碱基单链DNA探针与模板配对连接,发出不同的荧光信号,从而读取目标序列的碱 基排列顺序。在该方法下,目标序列的所有碱基都被读取了两遍,因此SOLiD最大的优势就是它的高准确率。据悉,SOLiD 5平台的测序通量已达到30Gb/天,成本低于60美元/Gb,准确率高达99.99%。并且由于SOLiD系统采用的不是PCR反应进行DNA合成与测 序,因此对于高GC含量的样本,SOLiD系统具有非常大的优势。
宣传视频:
Solexa
http://www.illumina.com/Media/flash_player.ilmn?dirname=systems&swfname=GA_workflow_vid&width=780&height=485&iframe
SOLiD
http://media.invitrogen.com.edgesuite.net/ab/applications-technologies/solid/solid-5500.html
NGS之进阶篇
454的焦磷酸测序原理,简单来说就是利用DNA聚合酶、ATP硫酸化酶、荧光素酶和双磷酸酶的协同作用,将PCR反应每一个碱基(dNTP)的延伸与一 次荧光信号的释放偶联起来,通过记录荧光信号的有无和强度,达到实时测定DNA序列的目的。在454测序仪中,A、T、G、C四种碱基是分别存储在单独的 试剂瓶中的,每步反应四种碱基依次加入反应池,当碱基配对结合,就会释放出一个焦磷酸(PPi),而这个焦磷酸在酶的作用下,将荧光素氧化成氧化荧光素, 并发出光信号,从而读取出这一位置的碱基信息。
454测序仪的整个实验步骤可大致概括为:样品处理、文库制备、emPCR、反应板准备、上机测序。样品处理主要是针对大片段的DNA分子,如基因组 DNA、Fosmid或BAC质粒等,利用超声或氮气打断将这些DNA分子片段化,然后采用琼脂糖凝胶电泳回收或磁珠纯化,选择500-800bp的 DNA片段。对于非编码RNA或PCR产物,则不需要这一步骤。文库制备包括接头连接和磁珠纯化两步,454的文库接头分A、B两种,各44bp,由 20bp的PCR引物、20bp的测序引物及4bp(TCAG)的“key”碱基构成,其中B接头的5’端带有生物素(Biotin)标记,用于磁珠纯化 步骤。经过磁珠结合与DNA变性之后,只有A+目的片段+B形式的连接产物得以富集,另两种形式(AA、BB)的产物都被去除。emPCR是454测序的 一个关键步骤,将富集到的文库与测序磁珠、各反应物混合,加入特定的矿物油和表面活性剂,再利用振荡器剧烈振荡,使反应体系形成油包水(water- in-oil)的稳定乳浊液。在理想条件下,每一个液滴,或称微反应器(microreactor)中将只包含一个磁珠和一条单链DNA,通过控制该步骤 的条件,1mL乳液中可以形成至少10的6次方个理想的微反应器。经过PCR扩增后,每一个磁珠上将形成密集的DNA簇,这些DNA序列完全相同,即可用 于后续的步骤。454测序的反应板称为PTP(Pico Titer Plate),含有350万个由光纤组成的小孔,每个孔的直径为29μm,而测序磁珠的直径为20μm,因此每个孔中仅能容纳一个磁珠。将磁珠与测序试剂 加入PTP中,使之可用于上机测序。测序步骤如前所述,四种碱基在泵的控制下依次加入反应板,反应完成后再洗去,每延伸一个或若干个碱基,就会发出一次光 信号,通过记录信号的有无和强度,即可测定DNA序列。
454测序准确度较高,当读长超过400bp时,其准确性仍能达到99%以上,主要的错误来自于同聚物,即相同碱基的连续延伸,如ATTTG这样一段序 列,A和G的读取没有问题,但T只记录了一次光信号,仅信号强度与ATG序列的T有所不同,因此同聚物越长,可能产生的误差就越大。目前,由于454测序 仪在读长上的明显优势,它在大基因组从头测序(de novo)、转录组分析、基因组结构分析等领域有着广泛的应用。
以下图片:1. 454 GS FLX测序仪外观;2. 3(4)种连接产物的磁珠纯化;3. PTP板与测序磁珠;4. 454 GS测序峰图;5. 上机准备实拍图。
NGS之进阶篇二:Solexa
Illumina Solexa高通量测序平台可以说是目前“测序界”应用最广泛的NGS平台,它在兼容性、操作性和成本方面有着较大的优势,第一个亚洲人基因组“炎黄一 号”、第一个非洲人基因组、熊猫基因组、家蚕基因组甲基化图谱等等,都是在该平台的支持下完成的。从最初的GA到GA IIx,又更新换代为现在的Hiseq2000,Solexa测序平台的通量由1Gb/run一路提升至300Gb/run,预计在年内还将实现1Tb /run的升级,测序读长也从几十bp提高到150bp。当年的人类基因组计划用了10年时间完成了一个基因组的精细图,而现在使用Hiseq2000测 序仪只需10天左右的时间,即可完成至少3个人类全基因组的测序工作,NGS测序能力的飞速发展甚至超过了IT届的摩尔定律。
与454一样,Solexa测序平台所采用的也是SBS(Sequencing-by-Synthesis,边合成边测序)的方式,并且也利用光信号收集 信息。有所不同的是,Solexa并没有采用间接反应(焦磷酸氧化荧光素)的形式激发光信号,而是直接在dNTP上连接荧光基团和阻断基团,通过“去阻断 —延伸—激发荧光—切割荧光基团—去阻断”这样一个循环的方法来依次读取目的DNA上的碱基排列顺序。如下图所示,该原理在基础篇的Solexa宣传视频 中亦有提及。由于采用了可逆阻断技术(即在dNTP上连接可剪切的阻断基团),Solexa测序的每一步只延伸一个碱基,不会出现类似于454测序的同聚 物影响准确性的问题,因此其单碱基准确性较高,但随着读长的增加,荧光信号会有所衰弱,所以越“靠后”的碱基准确性会逐渐降低,这也是Solexa测序读 长受限的一个主要因素。
Solexa平台的应用范围极广,几乎囊括了目前基因组学研究的所有方面,例如基因组从头测序(de novo)、重测序(re-sequencing)、基因组结构分析、转录组测序、表达谱分析、小RNA及非编码RNA测序、表观遗传学研究等等。然而, 应用该平台的核心过程是大致相同的,这也为它的兼容性提供了很大的便利。Solexa测序的实验流程主要包括:样品处理、文库制备、芯片准备及上级测序。 根据实验目的和样品来源的不同,Solexa测序的样品处理也有所不同,基因组DNA需进行打断及片段选择,total RNA需富集mRNA或小RNA,mRNA也要进行片段化处理,ChIP(染色质免疫共沉淀)、甲基化及PCR产物等都有各自的处理方式,需要实验人员根 据情况进行选择。以基因组DNA测序为例,在获得样本后,首先需要对DNA进行检测,保证样本的浓度和完整性符合实验要求,然后使用超声或氮气打断,将这 些DNA分子片段化,这步与454的样品处理类似,但不需要收集特定范围的DNA片段即可用于下一步实验。文库制备过程可分为末端修复、3’端腺苷化、接 头连接、片段选择、PCR扩增以及文库纯化六个步骤。末端修复是将片段化的DNA分子在几种酶的作用下,补齐随机打断造成的黏性末端而成为平末端。3’腺 苷化目的是在DNA双链的3’端加上dATP,以便于下一步的接头连接。Solexa文库制备所用的接头(adapter)是一个一端互补,另一端开叉Y 字形DNA片段,互补的部分5’端有一个T,可与3’腺苷化的DNA进行T-A连接,开叉的部分则是为了通过PCR扩增引入测序引物P5和P7的互补序 列。接头连接完成后,再采用采用琼脂糖凝胶电泳回收或磁珠纯化获得特定大小的片段(通常为目的片段大小+120bp),如构建200bp文库,则回收 320±20bp的DNA。之后再进行PCR扩增和纯化,即得到可用与上机测序的文库。此时的文库DNA由待测序列、测序接头、引物互补序列及index 标签序列(如非index文库则无),如下图所示。
芯片准备与上机测序主要由机器完成,在此不做赘述,若想做进一步了解,可访问Illumina官方网站的技术支持:http://www.illumina.co ... chnology.ilmn。
同样,最后献上几张图片:1. cBot芯片制备仪与Hiseq2000测序仪;2. Solexa GAII测序芯片;3. Solexa测序光信号采集直观图;4. Solexa测序流程与原理示意图。