跳到主要内容

瑞兴生物RXBio

多组学技术服务与深度数据挖掘 · 从研究设计到实验验证的全流程支持
表观和遗传组 · 染色质可及性

ATAC-seq 染色质可及性测序

DNA 调控元件有没有活性,一个关键特征在于染色质可及性——这些元件在染色质结构里能不能被接近。ATAC-seq 利用 Tn5 转座酶"只往开放区域里插接头"的特性,把全基因组的开放染色质区域直接测出来。

  • 全基因组开放染色质
  • 起始量低至数万细胞
  • motif 识别转录因子
  • 核小体定位
  • 多组学联合
≥2×10⁵ 活细胞冻存前活性 >90%、复苏后 >85%
2013年建立斯坦福 Howard Y. Chang 团队开发
TF足迹转录因子结合位点与组蛋白修饰位点
低线粒体污染有效降低线粒体 DNA 数据污染
01

技术简介

DNA 调控元件活性的一个重要特征在于其染色质可及性,也就是这些元件在染色质结构中的可接近性和可操作性。只有那些活跃的、具有活性的调节元件才有可能被细胞内部的机制所识别并进一步表达。因此,染色质的可及性在转录调控的过程中扮演着举足轻重的角色,为研究转录因子如何通过调节染色质可及性来影响基因表达提供了关键线索。

2013 年,美国斯坦福大学 Howard Y. Chang 研究团队开发了用于研究染色质可及性的方法 ATAC-seq,其原理是通过转座酶 Tn5 容易结合在开放染色质的特性,然后对 Tn5 酶捕获到的 DNA 序列进行测序。ATAC-seq 由于所需细胞量少、实验简单,可在全基因组范围内检测染色质可及性,已成为相关研究的首选技术方法。

02

技术优势

起始量低

所需细胞数 10,000–50,000 个,甚至可以做单细胞水平研究,对稀有细胞亚群与微量分选样本友好。

流程简单

步骤简单、成功率高,转座反应在建库的同时完成片段化与接头插入,无需额外超声与末端修复。

特异性强

有效降低线粒体 DNA 数据污染,让有限的测序数据用在真正有信息的核基因组区域上。

重复性好

实验重复结果一致性好,适合多样本平行比较与大队列研究。

03

技术应用

1. 绘制开放染色质图谱

对于研究基因组结构和功能、疾病诊断和治疗等方面都具有重要的意义。

2. motif 分析识别转录因子

识别参与基因调控的转录因子,对于理解基因表达的复杂机制以及与基因相关疾病的发生机制非常重要。

3. 识别转录因子靶基因

识别转录因子调控的靶基因和功能元件,研究转录因子在细胞生长、分化、凋亡等生命活动中的重要作用。

4. 核小体定位分析

结合核小体预测分析,精准地确定核小体定位的规律和影响因素,为研究染色质结构和功能提供更全面的认识。

5. 多组学联合分析

获得基因表达与上游调控元件、转录因子的关联信息,解析转录因子靶向调控变化机制。

04

技术流程

ATAC-seq 的核心技术酶是 Tn5 转座酶,它可以将测序接头高效地插入到染色质可接近区域。然而,染色质的其它区域则与组蛋白紧密结合,以核小体的形式存在,形成一种紧密的结构,限制了转座酶的插入。Tn5 转座酶能够精准地将测序接头直接插入到这些开放区的双链 DNA 中,对 DNA 进行片段化。随后,这些 DNA 片段经过提取和 PCR 扩增等后续处理,最终制备出适合二代测序的文库。

通过测序和数据分析,我们可以进一步获得有关开放染色质的重要信息,这包括了转录因子(TF)的结合位点、组蛋白修饰位点等关键生物学信息。

1

细胞制备

收集活细胞,控制活性与细胞数

2

转座反应

Tn5 转座酶切割开放染色质并同时插入接头

3

DNA 纯化

纯化转座产物,去除蛋白与杂质

4

PCR 扩增

扩增目标片段并构建测序文库

5

高通量测序

文库质检后上机测序

6

数据分析

Peak 调用、motif 与核小体分析

ATAC-seq 流程:Tn5 转座酶将测序接头插入开放染色质区域,PCR 扩增后测序
ATAC-seq:Tn5 转座酶插入开放染色质 → 纯化 → PCR → 测序 → Peak 与 motif 分析
05

样本要求

活细胞活性是关键指标:冻存前活性 >90%,复苏后活性 >85%。

样本类型送样量说明
活细胞≥2×10⁵ 个冻存前活性 >90%、复苏后活性 >85%;细胞悬液冰袋运输,参考血液类样本
新鲜动物组织解离 · 心、脑、肾、肺≥100 mg新鲜组织取样后尽快解离;当天干冰运输到实验室当天处理。易存在坏死和钙化的组织建议先解离再提核,可去除部分杂质、降低上机堵芯片风险
新鲜动物组织解离 · 肝、脾、胃≥100 mg
新鲜动物组织解离 · 皮肤、血管、软骨、骨等坚韧组织≥200 mg
冻存动物组织抽核 · 心、脑、肾、肺≥200 mg冻存组织需液氮速冻后 −80 ℃ 保存、干冰运输,一个月内处理;所有冻存样本均需提 RNA 质检是否降解
冻存动物组织抽核 · 肝、脾、胃≥200 mg
冻存动物组织抽核 · 皮肤、血管、软骨、骨等坚韧组织≥400 mg
植物组织≥0.5 g选健康幼嫩组织;70% 酒精冲洗表面后液氮速冻,干冰运输
三个容易踩的坑:① 冻存组织的送样量要求明显高于新鲜组织(抽核路线,心脑肾肺从 100 mg 提高到 200 mg),因为提核过程有损耗;② 细胞悬液浓度建议 700–1200 个/µL、活性 >90%、结团率 <10%,结团严重的样本 Tn5 反应效率会下降;③ 冻存样本转移至干冰运输时避免室温停留过久,样本溶解会导致活性下降与 RNA 降解。完整版要求以最新版《送样规范》为准。
06

分析内容

标准分析 STANDARD

  1. 原始数据处理fastqc 质控、trimmomatic 去接头与低质量碱基、MultiQC 汇总;Clean 数据统计含 Q20 / Q30 / Dup / GC
  2. 比对统计Bowtie2 比对;统计 total / organelle(线粒体·叶绿体)/ multiple / unique(去细胞器)/ rmdup
  3. 酶切片段长度分布统计 Tn5 插入片段长度,看无核小体、单核小体、双核小体三个峰的占比
  4. 信号分布reads 在 TSS 上下游 3kb、基因本体上下游 3kb 的密度分布热图
  5. 样本关系PCA 主成分分析 + Pearson 相关性热图,组内 R² 需大于 0.8
  6. Peak calling 与峰特征MACS2(qvalue < 0.05,narrow / broad 可选);峰数量(含启动子区峰数)、峰中心上下游 3kb 平均信号、峰长度分布
  7. Peak 注释与组内重叠ChIPseeker 注释 promoter / exon / UTR / intron / distal intergenic 与距 TSS 距离;DiffBind 组内重叠 Venn
  8. 富集与 Motif启动子区(TSS ±3kb)峰基因做 GO 与 KEGG;HOMER motif 富集,含 motif 对应峰与已知转录因子 motif
  9. 差异 PeakDiffBind 合并修剪 + DESeq2,FDR ≤ 0.05;输出火山图、维恩图、差异峰热图与注释表

高级分析 ADVANCED

  1. 与 RNA-seq 联合可及性 × 基因表达相关性,预测 peak 靶基因
  2. 与 ChIP-seq / CUT&Tag 联合开放区域 × 因子结合峰,定位活跃调控元件
  3. 核小体定位与染色质状态由片段长度与信号周期性推断核小体排布
  4. 转录因子足迹分析在开放区内找 TF 保护的碱基足迹,推断实际占位
07

实测数据

下面 5 张图来自我司 ATAC-seq 项目的真实测序结果,对应五项必看的质控与结果指标——插入片段大小、reads 在基因上的分布、Peak 基因组占比、Peak 重复性、peak 区域 IGV 可视化。拿到数据先按这个顺序过一遍,前两项判断实验是否成功,后三项判断结果能不能用。

ATAC-seq 实测数据:插入片段大小分布,显示无核小体、单核小体与双核小体区间的峰型
图 1:插入片段大小分布——无核小体 / 单核小体 / 双核小体峰型
ATAC-seq 实测数据:reads 在基因上的分布,TSS 与 TES 上下游 3kb 的信号密度曲线与热图
图 2:reads 在基因上的分布——TSS / TES 上下游 3 kb 信号
ATAC-seq 实测数据:Peak 在基因组上的占比,各样本 Peak 在 promoter、UTR、exon、intron 与远端基因间区的百分比堆叠柱状图
图 3:Peak 在基因组上的占比——各功能区域百分比分布
ATAC-seq 实测数据:Peak 在 2 个样本中的重复性,样本间 Peak 重叠的维恩图
图 4:Peak 在 2 个样本中的重复性——样本间 Peak 重叠
ATAC-seq 实测数据:Reads 在 peak 区域的 IGV 可视化结果,多样本信号轨道与基因结构注释
图 5:Reads 在 peak 区域的 IGV 可视化——信号峰与基因结构对齐
这五张图怎么看:① 插入片段大小分布——ATAC 的正常峰型应能看到无核小体(<100 bp)、单核小体(约 200 bp)、双核小体(约 400 bp)的周期性起伏,这是 Tn5 成功切割开放染色质的直接证据;峰型平坦或只剩大片段,提示细胞状态或裂解条件有问题。② reads 在基因上的分布——信号应在 TSS 附近明显富集,符合开放染色质集中在启动子的预期。③ Peak 基因组占比——启动子区占比高说明富集有效;若大量落在远端基因间区且信号弱,需结合重复样本判断是否为背景。④ Peak 重复性——重复样本间的 Peak 重叠率是判断数据可用性的直接依据,重叠过低的数据不建议用于差异分析。⑤ IGV 可视化——把信号峰落到具体基因座人工复核,确认峰形是否锐利、是否与已知开放区域吻合。
08

应用案例

案例一
IF 45.8 Science · 泛癌染色质可及性

The chromatin accessibility landscape of primary human cancers

人类原发癌症的染色质可及性图谱

期刊:Science(IF 45.8,2024 JCR) · 发表时间:2018 年 10 月
阅读原文:doi.org/10.1126/science.aav1898

样本与方法:2018 年 10 月 26 日,Howard Y. Chang 教授团队在国际顶级学术期刊 Science 发表了染色质可及性研究的最新成果。利用 ATAC-seq 方法,结合癌症基因组图谱研究的现有数据,研究人员在 410 例囊括 23 种癌症类型的冷冻样本中,定位了 562,709 个可重复、转座酶可接近的染色质可及性位点。

  • 与已知调控元件高度一致:泛癌分析发现的染色质可及性峰值中,约有三分之二与过去发现的调控元件信息一致,表明 ATAC-seq 技术能够很好地重现过去的研究发现,同时还能发现大量新的染色质可及性敏感位点。
  • 可及性数据可用于肿瘤分型:利用 Distal elements 和 Promoters 数据分型与 RNA-seq 数据分型的比较,说明 ATAC-seq 数据可以用于区分不同肿瘤类型。
  • 发现新的肾癌亚型:ATAC-seq 数据对 KIRP(一种肾癌类型)的分子分型研究发现,图中红色代表了一种新的亚型,该亚型 MECOM 基因的染色质结合活性异常升高,使 MECOM 基因高表达;MECOM 的过表达与较差的总生存率显著相关。
  • peak 与转录因子和甲基化密切相关:通过不同分型中特异性 peak 的寻找,发现这些特异性 peak 以转录因子结合位点为主,并和甲基化密切相关,同时发现了不同肿瘤类型特异的转录因子。
  • 可及性-表达相关性预测:利用 ATAC-seq 数据和 RNA-seq 数据计算染色质可访问性与基因表达的相关性,一共预测了 8,552 个蛋白编码基因表达调控中对应的 peak,不同肿瘤类型中特异的 peak 不同。
09

常见问题

ATAC-seq 和 DNase-seq、MNase-seq 怎么选?
三者都测染色质开放程度,但 ATAC-seq 优势明显:所需细胞量最少(数万个细胞即可)、操作最简单(转座反应一步完成片段化与接头插入)、背景低(有效降低线粒体 DNA 污染)。DNase-seq 需要大量细胞且酶切条件敏感;MNase-seq 更适合核小体定位而非开放区域图谱。常规开放染色质研究首选 ATAC-seq。
需要多少细胞?组织样本可以吗?
活细胞 ≥2×10⁵ 个(冻存前活性 >90%、复苏后 >85%)。组织样本分两条路线:新鲜组织当天干冰运输、当天处理,心脑肾肺与肝脾胃 ≥100 mg、坚韧组织 ≥200 mg;冻存组织走抽核路线,心脑肾肺与肝脾胃 ≥200 mg、坚韧组织 ≥400 mg,干冰运输、一个月内处理。冻存路线送样量更高是因为提核过程有损耗。
线粒体 DNA 比例高是什么原因?
常见原因是样本中死细胞比例高或细胞状态差——死细胞的线粒体膜完整性下降,线粒体 DNA 更容易被 Tn5 接触。对策是提高细胞活性(冻存前 >90%)、控制细胞悬液结团率(<10%),必要时在转座反应前做死活细胞富集。我们的流程已针对线粒体污染做了优化。
能直接看出转录因子结合位点吗?
ATAC-seq 的 peak 是开放染色质区域,不是直接的转录因子结合位点。要推断结合的转录因子,需要对 peak 区域做 motif 分析——如果某转录因子的识别基序在 peak 中显著富集,就提示它可能结合在这里。要确证还需要 ChIP-seq 或 CUT&Tag。
需要多少生物学重复?
差异可及性分析建议每组至少 2–3 个生物学重复。ATAC-seq 的重复性通常优于 ChIP 类实验,但组间差异峰仍依赖统计一致性;只做开放染色质图谱描述时,每组 2 个样本也可以启动。
能和 RNA-seq、CUT&Tag 一起做吗?
这是标准组合。ATAC-seq 给出"哪些区域开放",CUT&Tag / ChIP-seq 给出"哪个因子结合在哪里",RNA-seq 给出"基因表达变了没有"。三套数据联合可以直接回答"染色质开放 → 因子结合 → 基因表达"这条调控链,也是解析转录因子靶向调控机制最完整的方案。
Get A Solution

获得技术方案

您有任何关于组学测序的项目问题,请留言给我们。技术团队将在工作日 24 小时内给您答复。

项目咨询

技术顾问对接具体课题,先讨论科学问题再谈检测清单

027-87050299

sales@rxbio.cc

寄样收样

李欣 13237124769
武汉市东湖高新区高新二路 388 号光谷生物医药加速器 18 栋 1-2 层

提交您的项目需求

填写以下信息,我们会尽快与您联系并给出技术方案建议。

提交即表示同意我们通过电话或邮件与您联系,仅用于项目沟通。