目前主流三代测序平台除了Oxford 家的 Nanopore,还有 Pacific Biosciences(简称 PacBio)公司的 Single Molecule Real-Time(SMRT)Sequencing。该平台的优势在于:

  • 在不会影响吞吐量和准确性的前提下,提供目前最长的 25 kb 的 Reads 长度

  • 如果不含系统误差,准确度可达 99.999%,这样高质量的 Reads 可以解析几乎所有类型变体,从头组装高质量基因组

  • 可测取富含 AT 或 GC 区域,高度重复序列,回文序列等,不会产生 GC 的较大偏差

  • 可直接测取化学修饰,在表观遗传学中有重要应用

吃个瓜,2018年11月1日,Illumina 同意以 12 亿美元现金收购 PacBio 和其三代测序技术。

但是,去年 Illumina 放弃了收购计划,摊手。

接下来,我们看看它如何巧妙地完成这样的高质量三代测序。

1 基本原理

边合成边测序,与前文我们说的 Illumina 的基本测序原理一样。

2 构建文库

将样本中的 DNA 或 RNA 分子提取后,构建如下的哑铃状分子结构:

  • 黄色,紫色:双链 DNA 分子

  • 蓝色:接头(Adapter)

将文库分子展开,一个完整的圆环出现在我们眼前:

这种结构有利于进行周而复始的滚环复制,我们后文会讲这种复制方式的好处。

将样本中所有的DNA片段都构建哑铃状分子结构,组成的集合就叫文库(SMRTbell Library),随后,它们会被放到测序芯片中。

3 测序芯片

以 RSII 测序平台为例,测序仪芯片(SMRT Cell)长这样:

放大后:

上面整齐排列着15万个直径为70纳米的测序微孔(Zero-Model Waveguides,ZMWs)。

4 上机测序

1、构建测序复合物

测序复合物:聚合酶,测序模板,测序引物

2、复合物撒入测序小孔

3、固定测序复合物

由于聚合酶加了生物素,在芯片玻璃底板有链酶亲和素。利用生物素和链酶亲和素的亲和力,包含聚合酶的测序复合物会被固定在玻璃底板。

4、构建带有荧光基团的 dNTP

在芯片溶液中含有许多游离 dNTP,所谓游离 dNTP 就是随机飘在溶液中的 dNTP。

ATGC 四种碱基的 dNTP,在磷酸基团上分别带有四种颜色的荧光基团。

5、边合成边测序

在合成时,游离的 dNTP 被固定在底板上的酶捕获,激发光会从玻璃板底部发出。

怎么保证每次测取一个碱基?

由于测序小孔直径很小,激发光的穿透能力会逐渐衰减,只能在小孔中传输很短的距离,所以只有当 dNTP 足够靠近底部,荧光基团才会被激发光照到,发出荧光。当然,其他的游离 dNTP,虽然也有可能飘到小孔底部被激发光照到,但这种情况极少。

在一个碱基合成结束后,带有荧光基团的磷酸基团会从 dNTP 上掉落,发生猝灭,不影响其他碱基的信号检测。

在发生测序的小孔有各自的 DNA 片段和测序复合物,同一时间发出不同颜色的激发光,机器会检测到如下的光信号,实际同时会得到多达几万个光点。

重复上述步骤,经过计算机分析光谱,最终我们拿到样本的测序文件。SMRT Sequencing 测序过程中,每秒读取三个碱基,一个小时可检测大约一万多碱基。

6、检测碱基甲基化

有意思的是,在 SMRT Sequencing 测序过程中,可以直接测到碱基被修饰的状态,聚合酶遇到碱基上带有甲基化的碱基,合成速度会明显变慢,而且光谱也会发生改变。

因此,SMRT Sequencing 可以检测到碱基的甲基化修饰情况。

5 测序模型

SMRT 测序有如下两种测序模式:

1、Circular Consensus Sequencing (CCS)

说这种测序模型前,就不得不提三代测序最大的缺点:碱基读取不准,错误率在12.5%,也就是说,每读取八个碱基,就会读错一个。

好在碱基读取错误是随机的,如果重新读一遍同样位置的碱基,不一定会发生同样的错误。

如果对同一个序列,多测几遍,那么这些读错的碱基就能矫正过来。

前边提到的滚环复制的优势就来了,我们可以利用测序复合物在环状文库分子循环测序同一个片段来消除错误率。

这种测序模型,复制出的 Reads叫 HiFi Reads,测序准确率 > 99%。

2、Continuous Long Read (CLR) Sequencing

这种测序的优势在于可以读取更长的 Reads。

6 其他影响因素

1、GC bias 影响

什么是 GC bias?

PCR 时,如果模板里的 G、C 碱基含量高,PCR 效率低,A、T碱 基含量高,PCR 效率高。一般测序过程,如二代测序,都会有大量的 PCR 过程。这样就会有一个问题,G、C 含量高的片段,读到的 Reads 数少。

SMRT 在测序过程中,没有 PCR 过程,因此富含 GC 含量高,含量低的 Reads 片段都会有相似的概率被测序,所以三代测序中的 GC Bias 影响小。

2、读长的限制因素
  • DNA 模板断裂,用激发光长时间照射 DNA 链时,会发生断裂,DNA 链会从酶上掉下来,测序终止。

  • 酶变性,酶被长时间照射时,酶会变性,失去聚合酶活性,测序终止。

  • 文库序列短,如果做文库序列片段大于 20~30 K ,且保证质量的文库是有技术难度的

3、测序通量

目前,主流的测序平台由三种,各有利弊,可以根据自己的课题来选择。

以 RSII 为例,将测序复合物,随机撒到 15 万个小孔中,正好有一个复合物进入到单个小孔的概率符合泊松分布。理论情况是

  • 1/3 的小孔中有一个测序复合物,正常信号

  • 1/3 的小孔什么都没有,无信号

  • 1/3 的小孔中有两个以上的测序复合物,杂乱信号

五万个小孔 * 10kb,所以一张芯片大约会产出 500M 的数据。


参考:

  1. https://www.cnbc.com/2020/01/02/illumina-abandons-1point2-billion-deal-to-buy-rival-pacific-biosciences.html

  2. https://www.youtube.com/watch?v=rUKhfITd2CA

  3. https://www.youtube.com/watch?v=NHCJ8PtYCFc&t=1s

文末友情推荐

•《2021生信学习班起航,先送福利》•《96核心384G内存的超级服务器(共享)使用权一年

本文分享自微信公众号 - 生信科技爱好者(bioitee)。
如有侵权,请联系 support@oschina.cn 删除。
本文参与“OSC源创计划”,欢迎正在阅读的你也加入,一起分享。

图解三代测序(SMRT Sequencing)的更多相关文章

  1. 4、Brief primer and lexicon for PacBio SMRT sequencing

    转载:http://pacbiofileformats.readthedocs.io/en/5.1/Primer.html 转载:http://pacbiofileformats.readthedoc ...

  2. 用单分子测序(single-molecule sequencing)和局部敏感哈希(locality-sensitive hashing)来组装大型基因组

    Assembling large genomes with single-molecule sequencing and locality-sensitive hashing 好好读读,算法系列的好文 ...

  3. 三代PacBio reads纠错 - 专题

    三代纠错的重要性不言而喻,三代的核心优势就是长,唯一的缺点就是错误率高,但好就好在错误是随机分布的,可以通过算法解决,这也就是为什么现在有这么多针对三代开发的纠错工具. 纠错和组装是分不开的,纠错就是 ...

  4. DNA甲基化测序方法介绍

    DNA甲基化测序方法介绍 甲基化 表观遗传学 DNA 甲基化是表观遗传学(Epigenetics)的重要组成部分,在维持正常细胞功能.遗传印记.胚胎发育以及人类肿瘤发生中起着重要作用,是目前新的研究热 ...

  5. 测序原理 - PacBio技术资料

    手头有一套完整的PacBio技术资料,会慢慢的总结到博客上. 写在前面:PacBio公司主要有两个测序平台一个是RS,一个是最新的Sequel,下面如果没有指明则是在讲RS平台. SMRT测序技术总览 ...

  6. smrt analysis 软件安装

    pacbio 公司把三代测序常用软件整合起来,做成了一个在线的分析软件,类似于galaxy; 安装过程如下,参考官网给的指导手册 http://www.pacb.com/wp-content/uplo ...

  7. Pacbio 纯三代组装复活草基因组

    对于植物等真核生物基因组来说,重复序列, 多倍体,高杂合度等特征在利用二代数据进行组装的时候都会有很大的问题: 利用二代数据组装出来的基因组,大多达不到完成图的水准,通常只是覆盖到编码蛋白的基因区域, ...

  8. 第三代PacBio测序技术的测序原理和读长

    针对PacBio单分子测序——第三代测序技术的测序原理和读长     DNA基因测序技术从上世纪70年代起,历经三代技术后,目前已发展成为一项相对成熟的生物产业.测序技术的应用也扩展到了生物.医学.制 ...

  9. 植物基因组|注释版本问题|重测序vs泛基因组

    生命组学: 细菌和其他物种比,容易发生基因漂移,duplication和重排. 泛基因组学研究的一般思路是通过comparison找到特殊基因区域orspecific gene,研究其调控机制(即通过 ...

  10. Iso-Seq学习

    SMRT portal安装教程: http://www.pacb.com/wp-content/uploads/2015/09/SMRT-Analysis-Software-Installation- ...

随机推荐

  1. cost function 成本函数

    cost function 成本函数 cost function-成本函数 1.目标 :实现和探索具有一个变量的线性回归的成本函数. import numpy as np %matplotlib wi ...

  2. 一连串div跟随鼠标移动

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  3. 重新编译Unity Mono遇到的坑,大坑,巨坑!!!

    因为项目需要,要整一个DLL热更功能,本身也想对DLL进行加密,所以需要修改Mono的代码,并重新编译.参考了网上的教程后,决定在linux下做重编的工作. 然而在实际执行的时候,一开始就碰到了巨大的 ...

  4. 重磅!AWS升级对Apache Hudi的集成

    全球最大云厂商AWS的 Athena 团队又更新了 Athena 与 Apache Hudi 的集成,以支持新功能及最新的 0.8.0 社区版本.早在Apache Hudi还处于孵化阶段时,AWS A ...

  5. ggplot2图形可视化应用集锦

    数据可视化就是将我们从数据中探索的信息与图形要素对应起来的过程.数据可视化,先要理解数据,再去掌握可视化的方法,这样才能实现高效的数据可视化.数据可视化技术的基本思想,是将数据库中每一个数据项作为单个 ...

  6. [Java]变量及其初始化 与 类对象的初始化

    1 变量 1.1 变量的[定义] 1.2 变量的[作用域] 1.3 变量的[初始值] 1.4 补充:缓存变量 1.5 变量的[分类]与[未初始化情况] 2 类对象 2.1 类对象的初始化/构造过程 1 ...

  7. LeeCode 栈与队列问题(一)

    LeeCode 20: 有效的括号 题目描述 给定一个只包括 '(', ')', '{', '}', '[', ']' 的字符串s,判断字符串是否有效. 有效字符串满足: 左括号必须用相同类型的右括号 ...

  8. 访问不到tomcat下的静态资源文件

    问题:idea+tomcat 测试Servlet,居然访问不到与WEB-INF同级的静态资源1.html.又不是访问WEB-INF里的文件.... 分析:可能是开发目录中的文件没有自动复制到部署目录. ...

  9. TEC-6计算机组成原理实验(简图)

    TEC-6计算机组成原理实验

  10. 面试最常问的数组转树,树转数组 c++ web框架paozhu实现

    刚毕业同学,找工作常被问 二维数组转树,树转二维数组 需要支持无限层级实现,如果你了解这个语言那么实现起来还要一番思考 c++ web框架 paozhu使用 需要实现数据库表数据到前台菜单实现,就是这 ...