单细胞RNA测序(scRNA-seq)细胞分离与扩增

单细胞RNA测序入门可以查看以下文章:

单细胞RNA测序(scRNA-seq)工作流程入门

1. 单细胞的分离

如何获得单细胞,从而进行下一步的测序过程?具体有以下几种方法:

一、流式细胞仪(FACS)方法

常用的方法之一是使用流式细胞仪(FACS)方法进行分离。但是 FACS 也有缺点,首先需要较大的起始细胞量,同时可能导致有些孔中有两个,有些孔中没有细胞。

二、显微操作分离

显微操作提供了获得少量细胞的替代方法,我们可以通过直接特定位置摘取单个细胞进行测序,但是这种方法依赖于人工操作,无法实现高通量进行。

三、微流控技术

近些年微流控技术在单细胞分选上展现良好的应用,微流控技术可以对整个过程进行监控和处理,其缺点主要在于固定的微流控芯片等。通过多通道微流控自动处理技术进行处理,可以允许自动分离和处理细胞,同时控制细胞的数量和实验的误差

2. 全基因组扩增

对于一个细胞(哺乳动物细胞)而言,通常细胞DNA不到10pg,需要先对全基因组进行扩增才能开展后续的实验。全基因组扩增技术(WGA)主要是通过 PCR多重替换扩增技术(MDA)两者的组合来实现。

全基因组扩增策略

一、PCR扩增策略

PCR 策略主要是通过使用随机引物进行扩增,然后引入通用 接头(adaptor)。随机引物 PCR 方法相对于MDA,有更高水平的一致扩增, 但是覆盖度较低。

二、多重替换扩增技术(MDA)策略

MDA 方法则利用具有链置换活性的 DNA 聚合酶(一种滚环复制酶),进行恒温扩增, 可以处理长达 10 Kb 的序列,一旦接触到前面结合的序列后,聚合酶可以释放前者 DNA,从而释放用于进一步的扩增。MDA 方法由于 DNA 聚合酶的高保真性和高处理特性,可以实现较好的覆盖度也有利于 SNP 等精细实验,然而由于 GC 含量等影响,MDA 可能导致不一致扩增,Gcbias 比较大,因此 MDA 用于 CNV 检测等的假阳性概率比较大

三、MALBAC 策略

华人科学家谢晓亮等人基于前两种方法提出综合两种方法的 MALBAC 方法,实现了高覆盖度的一致性扩增方法。

通过增加互补末端引物使用滚环复制酶产生可以末端配对的 loop,被扩增出来的 DNA,由于末端配对,因此无法作为接下来的模板,从而可以实现近乎线性的扩增,经过五轮线性扩增后,然后通过 PCR 进行指数级扩增,相对于 MDA 方法 31%-61% 的丢失率,MALBAC 方法只有 1% 的丢失率,同时对于 SNP 和 CNV 的检测的假阳性更低

MALBAC 策略

3. 单细胞DNA-seq基本工作流

单细胞基因组测序后,需要剔除低质量的reads,比对到参考基因组,校正由GC偏好引起的reads数的误差,后续可以使用GATK工具进行Variant calling。

不同单细胞 DNA 测序的方法特征对比

不同单细胞 DNA 测序的方法

4. 单细胞RNA-seq基本工作流

单细胞 RNA-seq 的主要障碍是如何从单个细胞中获取足够的 RNA 用于下一步的测序,现有多种的方法扩增得到低于 pg 级的 mRNA。扩增mRNA的主要障碍是扩增的偏好(bias)问题,其会对来自不同基因的 mRNA 之间的差异比较产生干扰

一、单细胞RNA-seq扩增技术的发展

利用单细胞进行 RNA-seq 是 Surani 实验室在 2009 年开发出来的,扩增的方法主要是使用将 RNA pull down 下来,然后使用 polyT 引物反转录带 ployA 的 RNA,同时引物上含有特定的 anchor 序列

研究人员借鉴了基于 microarray 用于 single cell sequence 的技术后,修改了一些实验条件(比如延伸时间等),从而获得全长的 cDNA,进一步通过 SOLID 测序,我们可以看到具有非常好的稳定性。

2011 年科学家们使用了新的 cell barcode 方法,从而可以探究多个复杂细胞的特性,称为 STRT 技术,但是 STRT 技术对于 5’end 具有偏好性基于转座子酶,人们开发了升级版的 STRT2 技术

为了减少由于扩增引起的误差,人们在一些单细胞测序的步骤中增加了 UMI(unique molecular identifiers),UMIs 是由 4-10 个随机核苷酸组成的序列组成,在 mRNA 反转录后,加入到文库中,每一个 mRNA随机连上一个 UMI,因此可以计数不同的 UMI达到计数 mRNA 数量的目的

二、单细胞RNA-seq扩增技术的选择

如果我们想了解不同细胞间的异质性,那么需要减少实验的 bias,推荐使用 UMI 相应的技术,如果我们想了解 RNA 全长,包括 剪切位点 等信息,那么我们就需要选择包含全部转录长度的技术

三、UMI计数

如果我们使用 UMI 计数方法( UMI 计数,是指在建库过程中,通过在引物上,增加随机序列,则对于同一种 mRNA 连上同样的 UMI 概率几乎为 0,则可以忽略由于 PCR 造成的误差,对于一种 mRNA,测到的 UMI 数量可以近似看成 mRNA 的表达个数),推荐至少每一个转录本至少覆盖 3-4 倍确保一些低表达的 mRNA 不会因为 噪声而被忽略,一般而言,illumina 的一条 lane 中最终可以被使用的数据大概占 50%,另一些 reads 因为实验技术所限,为一些 接头adaptor 等无效序列。

假设想获得1000个转录本信息,以nextseq一个lane可获得200M条reads,需要进行检测细胞数量则= 200,000,000 / 1000 / 4 / 2 = 250个细胞。

5. 单细胞RNA-seq的数据挖掘

单细胞 RNA-seq 为我们了解组织的组成等提供非常好的工具,而单细胞测序的效果也取决于我们到底测多少细胞,一般而言几百个细胞的单细胞测序不仅仅捕获常见的细胞成分,也会捕获稀少类型的细胞。由于 RNA 降解,细胞的凋亡和实验误差等原因,一些数据将会被排除,不纳入数据分析环节,因此初始选择的细胞个数往往会大于我们目标测序的细胞个数

对于不同转录本,没有要求的话,推荐大家使用多转录本合并的序列。因为大部分单细胞测序主要是富集 3‘UTR 部分,那么如果使用 cufflink,就会造成很大的问题,比对完后,需要确定每个转录本的表达量,前面我们说到UMI用于 RNA 的计数,在计数前,还需要对不同的细胞进行滤过。

由于存在 RNA 降解等多种问题,需要考虑到去除一些低转录本等情况,同时还需要考虑去除低扩增效率的细胞。需要做的便是归一化,进而可以在不同细胞之间的进行横向比对。将所有的 cell 归一化到 reads 总数是中位数的细胞的 reads 总数中,或者将所有细胞归结成同样的 reads,被称为 down-sampling。

在获得不同细胞的不同基因的表达矩阵后,便可进行细胞分类注释及下游分析。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/503710.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

深度学习基础模型之Mamba

Mamba模型简介 问题:许多亚二次时间架构(运行时间复杂度低于O(n^2),但高于O(n)的情况)(例如线性注意力、门控卷积和循环模型以及结构化状态空间模型(SSM))已被开发出来,以解决 Transformer 在长…

LabVIEW转动设备故障诊断系统

LabVIEW转动设备故障诊断系统 随着工业自动化技术的不断进步,转动设备在电力、化工、船舶等多个行业中扮演着越来越重要的角色。然而,这些设备在长期运行过程中难免会出现故障,如果不能及时诊断和处理,将会导致生产效率下降&…

05. 【Android教程】Android 程序签名打包

在上一章,我们创建了自己的 Android 工程,并成功的在模拟器中运行起来。同时提到,工程目录中有一个 bin 目录,运行之后我们可以在此目录下找到我们的 apk。那么不难想到,我们在点“Run”之后,系统会编译我们…

[技术笔记] Flash选型之基础知识芯片分类

1、按照接口分类 分为 Serial串口Flash 和 Parallel并口Flash; 市场大量使用Serial Flash;价格便宜;已满足系统对数据读写速度的要求; Serial Flash已经可以代表 NOR Flash; 小知识: 1)在…

深度学习算法概念介绍

前言 深度学习算法是一类基于人工神经网络的机器学习方法,其核心思想是通过多层次的非线性变换,从数据中学习表示层次特征,从而实现对复杂模式的建模和学习。深度学习算法在图像识别、语音识别、自然语言处理等领域取得了巨大的成功&#xf…

Linux基础篇:VMware虚拟机3种常用的网络模式介绍

VMware虚拟机3种常用的网络模式介绍 VMware虚拟机提供了几种不同的网络连接模式,以满足不同场景下的网络需求。以下是VMware虚拟机的三种主要网络模式: 1.桥接模式(Bridged Mode)网卡名称VMnet0 桥接模式允许虚拟机直接连接到物…

鸿蒙OS开发实例:【瀑布流式图片浏览】

介绍 瀑布流式展示图片文字,在当前产品设计中已非常常见,本篇将介绍关于WaterFlow的图片浏览场景,顺便集成Video控件,以提高实践的趣味性 准备 请参照[官方指导],创建一个Demo工程,选择Stage模型熟读Har…

解决前后端通信跨域问题

因为浏览器具有同源策略的效应。 同源策略是一个重要的网络安全机制,用于Web浏览器中,以防止一个网页文档或脚本来自一个源(域、协议和端口),获取另一个源的数据。同源策略的目的是保护用户的隐私和安全,防…

PostgreSQL到Doris的迁移技巧:实时数据同步新选择!

PostgreSQL可以说是目前比较抢手的关系型数据库了,除了兼具多样功能和强大性能之外,还具备非常优秀的可扩展性,更重要的是它还开源,能火不是没有理由的。 虽然PostgreSQL很强大,但是它也有短板,相对于专业…

【Java数据结构】关于栈的操作出栈,压栈,中缀表达式,后缀表达式,逆波兰表达式详解

🔥个人主页:努力学编程’ 🔥内容管理:java数据结构 上一篇文章我们讲过了java数据结构的链表,对于链表我们使用了它的一些基本操作,完成了扑克牌小游戏的操作,如果你感兴趣的话,点…

数组类模板(进阶版)

目录 介绍: 分析: 实现: .hpp框架创建 .hpp函数内容 有参构造 拷贝构造: 重载 插入数据 删除数据 通过下标访问 获取数组大小 获取数组容量 析构函数 .cpp框架 int类型数据测试 char类型测试 总代码 .hpp代码 …

是德科技keysight N9000B 信号分析仪

181/2461/8938产品概述: 工程的内涵就是将各种创意有机地联系起来,并解决遇到的问题。 CXA 信号分析仪具有出色的实际性能,它是一款出类拔萃、经济高效的基本信号表征工具。 它的功能十分强大,为一般用途和教育行业的用户执行测试…

wireshark 使用

wireshark介绍 wireshak可以抓取经过主机网卡的所有数据包(包括虚拟机使用的虚拟网卡的数据包)。 环境安装 安装wireshark: https://blog.csdn.net/Eoning/article/details/132141665 安装网络助手工具:https://soft.3dmgame.com/down/213…

【LIMS】CMA与CNAS:中国认证体系中的两大支柱

目录 一、CMA:[中国计量认证](http://cma-cma.org.cn/)什么是CMA?CMA的作用 二、CNAS:[中国合格评定国家认可委员会](https://www.cnas.org.cn/)什么是CNAS?CNAS的作用 三、CMA与CNAS的关系相互促进共同目标 结语系列文章版本记录…

TCP网络协议栈和Posix网络部分API总结

文章目录 Posix网络部分API综述TCP协议栈通信过程TCP三次握手和四次挥手(看下图)三次握手常见问题?为什么是三次握手而不是两次?三次握手和哪些函数有关?TCP的生命周期是从什么时候开始的? 四次挥手通信状态…

git基本操作二(小白快速上手)

1、前言 接上篇我们接着来继续讲 2、.gitignore忽略文件 创建一个.gitignore文件,并将其置于项目的根目录下,Git将自动识别并根据该规则忽略相应的文件和目录。 # 忽略所有的 .log 文件 *.log# 但跟踪所有的 build.log 文件 !build.log# 忽略所有的 /lo…

lookup函数

lookup函数 单条件查询 示例 扩展多条件 扩展

文件的顺序读写——顺序读写函数——fgets、fgetc、fputs、 fputc

✨✨ 欢迎大家来到莉莉的博文✨✨ 🎈🎈养成好习惯,先赞后看哦~🎈🎈 目录 一、fgetc和fputc函数 1.1 fputc 1.2 fgetc 二、fputs和fgets函数 2.1 fputs函数 2.2 fgets函数 一、fgetc和fputc函数 1.1 fputc 返回类…

结构体类型,结构体变量的创建和初始化 以及结构中存在的内存对齐

一般结构体类型的声明 struct 结构体类型名 { member-list; //成员表列 }variable-list; //变量表列 例如描述⼀个学⽣: struct Stu { char name[20]; //名字 int age; //年龄 char sex[5]; //性别 }; //结构体变量的初始化 int main() { S…

鸿蒙OS开发实例:【工具类封装-页面路由】

import common from ohos.app.ability.common; import router from ohos.router 封装app内的页面之间跳转、app与app之间的跳转工具类 【使用要求】 DevEco Studio 3.1.1 Release api 9 【使用示例】 import MyRouterUtil from ../common/utils/MyRouterUtil MyRouterUtil…