我们为生物学家提供了一种无需编码的工作流程,只需使用基于浏览器的工具即可识别组织特异性基因增强子。我们的方案利用公开的 H3K4me1/H3K27ac 组蛋白标记和 Hi-C 数据,使没有编程专业知识的研究人员能够访问、分析和识别与其目标基因相关的潜在调节元件。
方法文章
我们为生物学家提供了一种无需编码的工作流程,只需使用基于浏览器的工具即可识别组织特异性基因增强子。我们的方案利用公开的 H3K4me1/H3K27ac 组蛋白标记和 Hi-C 数据,使没有编程专业知识的研究人员能够访问、分析和识别与其目标基因相关的潜在调节元件。
增强子是调节基因表达的 DNA 区域。增强子内的突变可导致异常基因调控,从而导致疾病。因此,鉴定调节特定组织中基因活性的增强子对于理解疾病的遗传基础至关重要。然而,增强子很难识别,因为它们不编码蛋白质。虽然有许多增强子存储库和鉴定工具可用,但这些工具的复杂性可能会给生物学家带来挑战。为了帮助生物学家使用这些资源,我们提出了一种生物学家友好的方案 (https://github.com/Ramialison-Lab/EnhancerWorkflow),该方案利用现有的基于网络的基因组学数据,如 H3K4me1 和 H3K27ac 组蛋白标记和染色质构象分析 (Hi-C) 数据来发现与增强子活性的目标组织中感兴趣的基因 (GoI) 相关的增强子。该协议完全基于 Web,不需要最终用户的编程技能。我们通过表征调节 TBX5 的候选增强子来证明这种方法的实用性,TBX5 是心脏发育的关键基因。该方案有助于识别左心室中与该基因相关的增强子。
增强子是调节基因转录、发育和细胞分化的非编码 DNA 区域 1,2。增强子突变可导致各种疾病,包括发育障碍、癌症和其他遗传病 3,4,5,6。因此,了解增强子对于了解基因表达、突变和疾病至关重要。
要了解增强子如何与其靶基因相互作用,确定它们在基因组中的位置非常重要。然而,识别增强子的位置并不总是那么简单,因为增强子既可以位于转录起始位点 (TSS) 附近,也可以位于更远的地方,跨越数十到数百千碱基 2,7,8,9。
尽管增强子的基因组位置不可预测,但它们显示出独特的生化和结构特征,因此可以系统地追踪它们。通常,增强子往往在基因间和内含子区域富集,在外显子中发现少量增强子 2,8。它们通常以特异性组蛋白修饰和转录因子结合为标志,这定义了它们的调节作用并决定了它们在不同发育阶段和组织中的时空活动10,11。
ChIP-seq 用于鉴定转录因子结合位点 (TFBS) 和组蛋白修饰标记,例如增强子标志、H3K4me1、活性增强子标志、H3K27ac 和 H3K4me3 标记,这些标记在启动子区域1、12、13、14、15 处富集。染色质构象捕获 (3C) 技术及其衍生物,如 4C、5C、Hi-C 和 ChIA-PET,用于绘制遥远基因组区域之间的物理相互作用。3C 靶向特定组织中的特定相互作用,而 Hi-C 则提供跨细胞类型的全基因组结构16,17。
除了目前的方法外,还开发了专门的方法来表征增强子,包括合并增强子数据库,例如 EnhancerAtlas 或 EnhancerFinder18,19。然而,这些工具通常需要研究人员整合多个数据集来研究许多组织中的多个增强子,这对于没有生物信息学和数据挖掘经验的生物学家来说可能会让人不知所措。
在这里,我们描述了一种用户友好的协议来选择增强子,完全基于现有的 Web 工具。这使研究人员能够查询目标基因 (GoI) 并检索相应的增强子。此处的方案根据一组特定的标准选择增强子:组蛋白修饰、染色质相互作用和组织特异性 1,12,13,14,15,16,17,20,21。与位于基因22 之间的基因组区域的基因间增强子相比,在内含子中发现的增强子更有可能显示组织特异性活性。为了确保全面覆盖潜在的活性增强子,我们定义了两个相邻 GoI 之间的搜索范围,以增加捕获位于基因体外部的调节元件的可能性。我们使用增强子特异性表观遗传标志 H3K4me1 和活性增强子标志 H3K27ac 来列出候选增强子。然后根据 Hi-C 数据对这些候选子进行提炼,保留与相应启动子发生物理相互作用的增强子。该协议旨在指导生物学家仅使用公开可用的基于 Web 的工具完成增强子鉴定过程。通过整合表观遗传学和染色质相互作用数据,此处描述的方法提供了一个实用的框架,以生成有关潜在增强子的假设,以便进一步实验验证。
注意:https://github.com/Ramialison-Lab/EnhancerWorkflow 上提供了分步演练。 表 1 和 表 2 总结了协议中使用的数据。 补充文件 1 中提供了故障排除。
1. 定位 GoI(图 1)
2. 定义增强子检测区域(图 2)
3. 组蛋白标记分析(图 3)
4. 染色质构象捕获 (Hi-C) 分析(图 4)
为了说明所提出的方案的使用,我们研究了人类心脏中的 TBX5 基因,使用涉及 H3K4me1、H3K27ac 和 Hi-C 数据的综合工作流程探索 TBX5 相关增强子。 TBX5 是一种有助于肢体和心脏发育的基因,包括四个腔室的形成和隔膜分离24。该基因突变是 Holt-Oram 综合征 (HOS) 的主要原因,该综合征会导致肢体异常和先天性心脏病 (CHD),包括隔膜缺陷24。 TBX5 相关心脏增强子的突变可能对 CHD24 产生严重影响。之前的一项研究在人心脏特异性组织中发现了三种已知的 TBX5 增强子——即“增强子 2”、“增强子 9”和“增强子 16”(补充文件 2),它们被证明在转基因小鼠中具有比较表型23。
我们研究了 RBM19 和 TBX3 之间富含 H3K4me1 和 H3K27ac 的区域,这是人类 TBX5 下游和上游的两个侧翼基因,以检索 TBX5 基因座处的假定增强子(图 1 和 图 2)。为了识别心脏特异性增强子,选择了心肌细胞。将推定的 TBX5 心脏增强子区域检索为坐标 (chr12: start-end),并确定了 22 个 H3K4me1 和 H3K27ac 相关区域 (图 3 和 补充文件 3)。从 EnsEMBL 基因组数据库中检索推定的 TBX5 心脏增强子,以交叉引用 4DNucleome 数据库中保存的 Hi-C 数据(图 4)。这样做是为了衡量潜在增强子和 TBX5 心脏启动子之间可能的相互作用。按照此处描述的方案,证实 22 个基因组区域中的 21 个与心肌细胞中的 TBX5 启动子 (chr12: 114400143-114410103) 相互作用(补充文件 4)。有一个区域与启动子没有物理相互作用(图 4,步骤 4.8)。最后,我们将该方案与这些经过生物学验证的增强子和当前心脏增强子的金标准数据库 VISTA 心脏增强子浏览器进行了比较,并揭示了数据库25 当前未捕获的其他增强子。
我们将此处介绍的方案检索到的 21 个 TBX5 增强子与现有数据库进行了交叉比较。我们从 VISTA 心脏增强器浏览器 (补充文件 5)25 中检索到 4 个 TBX5 增强子。在 VISTA 鉴定的 4 个心脏增强子中,3 个增强子 hs2329 、 mm1282 和 m370 与该基于 Web 的增强子检测方案鉴定的区域重叠(图 5)。每个预测的增强子还与先前实验验证的增强子共享基因组区域,这些增强子来自 Smemo 等人23、增强子 2 (chr12:114025907-114026275, GRCh38) 和增强子 16 (chr12:114415466-114420433, GRCh38),而它们没有显示与增强子 9 (chr12:114263402-114266886, GRCh38) 的重叠。作为 VISTA 鉴定的增强子之一,hs498 与该方案或 Smemo 等人经实验验证的增强子23 的任何预测增强子不重叠(图 5),即使该区域显示与 H34Kme1 标记部分重叠(图 5)。同样,增强子 9 与该管道预测的增强子不重叠,但与 H3K4me1 标记相关(图 5)。

图 1:在 EnsEMBL 基因组浏览器中查找 GoI 的分步指南。 用户首先打开 EnsEMBL 主页 (1.1),选择物种 (Human),然后在搜索栏中输入基因 (1.2-1.3)。从结果列表中,选择适当的基因 ID (1.4),这将打开基因摘要页面。然后,用户单击 Region in Detail 超链接 (1.5) 以可视化 GoI 周围的基因组区域,包括相邻元件和调节特征。 请单击此处查看此图的较大版本。

图 2:使用 EnsEMBL 基因组浏览器定义 GoI 周围的增强子检测区域。 为了定义增强子检测区域,使用 来自 GENCODE 轨道的基本基因注释 识别 GoI 两侧的两个相邻基因,其中基因显示为用合并的 EnsEMBL/Havana 注释标记的深黄色块。每个基因的转录方向由基因名称 (2.1 ) 旁边的箭头 (< 或 >) 表示。要选择相邻基因之间的基因间区域,请单击并拖动感兴趣的区域,然后在弹出框中选择 Jump to region 以放大 (2.2)。要添加与监管或增强器相关的注释,请单击 Add/remove tracks (2.3)。 请单击此处查看此图的较大版本。

图 3:使用 EnsEMBL 基因组浏览器在增强子检测区域中配置组蛋白修饰轨迹。 在左侧工具栏中,单击 “配置此页面 (3.1)”以访问轨道配置面板,并导航到“调节”部分 (3.2) 下的“按细胞/组织划分的活动”。在打开的选项卡中,选择“实验”部分 (3.3),然后使用细胞/组织搜索栏查找并选择您感兴趣的组织(心肌细胞)(3.4)。在组蛋白标记面板 (3.5) 中,启用 H3K4me1 和 H3K27ac 作为活性增强子标记,启用 H3K4me3 作为启动子标记,然后单击“配置轨迹显示”(3.6)。确认轨迹选择后,单击“查看轨迹”(3.7) 返回基因组查看器。组蛋白标记峰现在在检测区域 (3.8) 中显示为相应组织标记下的彩色块(黄色:H3K4me1,蓝色:H3K27ac,橙色:H3K4me3)。“Hists & Pols”弹出窗口包含该区域的基因组坐标,以碱基对(chr:start-end)为单位,可以复制并保存以供下游分析。点击轨道中的彩色元素后,会出现一个 “Hists & Pols” 弹出窗口。弹出窗口包含碱基对中区域的基因组坐标(例如,启动子区域的 chr12:11443450-114451611),可以复制并保存以供下游分析 (3.8)。同样,要提取候选增强子,请优先考虑 H3K4me1 和 H3K27ac 峰重叠的区域,如跨轨迹的峰和框的垂直对齐 (3.9) 所示。可以通过单击重叠区域的框直接选择重叠区域,也可以通过手动单击并拖动对齐的峰来定义区域(例如,chr12:114400143-114410103 表示活动候选区域)。弹出窗口中显示的坐标应以 BED 格式保存,以便进行下游验证或可视化。 请单击此处查看此图的较大版本。

图 4:使用来自 4D 核组数据门户的 Hi-C 热图可视化启动子-增强子染色质相互作用。 4D 核组数据门户主页显示一个堆叠条形图,按生物体汇总了可用的实验类型。通过单击条形图的相应部分 (4.1) 来选择人类样品的“原位 Hi-C”数据集。将显示相关数据集的筛选列表;选择源自分化为心脏成肌细胞的 H9 细胞的 Hi-C 数据集 (4.2)。所选数据集 (4.3) 将通过 Explore Data 按钮 (4.4) 在 HiGlass 浏览器中打开。在坐标框 (4.5) 中输入感兴趣的基因组区域,并将接触矩阵呈现为彩色缩放热图。较深的颜色(深红色到黑色)表示较强的染色质接触频率,而较浅的颜色(白色到橙色)表示较弱的相互作用。在启动子坐标处放置水平尺,在三个经过实验验证的对照增强子 (4.6) 的位置绘制垂直尺。这些交集用于定义严格的交互阈值,由启动子-增强子触点中最强的可见信号(最深的颜色)设置 (4.7)。在候选 H3K27ac 和 H3K4me1 标记的增强子的位置绘制额外的垂直规则(从步骤 3.8 开始)。启动子-增强子交叉点等于或暗于阈值的候选子被保留,而信号较弱(较浅的色方块)的候选子被排除在外 (4.8)。手动提取保留的坐标并以 BED 格式保存以供下游分析。(a. 增强子 2、b. 增强子 9 和 c. 增强子 16) 请单击此处查看此图的较大版本。

图 5:预测的 TBX5 增强子与 VISTA 验证的心脏增强子和对照增强子的基因组浏览器视图。 基因组浏览器快照显示增强子搜索范围 (STEP 2),将基于 Web 的方案检索到的预测增强子(底部)与 VISTA 验证的增强子(顶部)和实验验证的对照增强子(中)进行比较。主面板显示了带有注释调节元件的完整基因组位点,包括心肌细胞特异性 H3K4me1(黄色)、H3K27ac(蓝色)和 H3K4me3 峰(橙色)。三个放大图捕获了方案检索到的 VISTA 和对照增强子之间的比对。与控制增强器的重叠用红色框勾勒出来。每个子区域的坐标都显示在浏览器下方的面板中。 请单击此处查看此图的较大版本。
表 1:研究中使用的数据。请点击此处下载此表格。
表 2:研究中使用的基于 Web 的工具。请点击此处下载此表格。
补充文件 1:EnsEMBL Genome Browser 的故障排除说明。 请点击此处下载此文件。
补充文件 2:GRCh38 格式的 BED 文件,经实验验证的 TBX5 心脏控制增强剂23。请点击这里下载此文件。
补充文件 3:GRCh38 格式的 BED 文件, 通过 STEP3 从 EnsEMBL 检索的 TBX5 心脏增强剂。请单击此处下载此文件。
补充文件 4:GRCh38 格式的 BED 文件, 通过 STEP4 从 EnsEMBL 检索的 TBX5 心脏增强剂。请单击此处下载此文件。
补充文件 5:GRCh38 格式的 BED 文件,从 VISTA 心脏增强器浏览器中检索的 TBX5 心脏增强器 25。请单击此处下载此文件。
此处描述的基于 Web 的协议用作增强子检索的工作流程,而不是增强子预测。通过利用公开可用的组织特异性数据集、组蛋白修饰(H3K4me1 和 H3K27ac)和 Hi-C 相互作用数据,它缩小了与 GoI 相关的潜在增强子的范围。与依赖于机器学习或基于序列模型的计算预测工具不同,我们的方法仅侧重于根据实验数据检索增强子。 补充文件 1 中提供了 EnsEMBL 基因组浏览器的故障排除说明。
该工作流程集成了组蛋白标记和染色质相互作用数据等关键元素,类似于 ChIA-PET 和 PLAC-seq 等高级方法,它们以更高的准确性绘制增强子-启动子相互作用10,11。然而,当高分辨率实验技术不可行时,这种方法是有益的,因为该方案的资源密集度较低,可以节省大量时间。
上述方法的主要限制是它依赖于现有数据集的可用性和质量,这可能会影响检索到的交互的精度。在 TBX5 相关基因突变的背景下研究心脏发育过程中的增强子活性需要组织特异性分辨率。对于此类分析,胚胎心脏组织将是最合适的,因为它与发育调节相关。然而,在分析时,没有包含组蛋白修饰数据的胚胎数据集是公开的。为了解决这个问题,集成 ENCODE、Enhancer Atlas 2.0 或 EnhancerFinder 等替代资源可以通过为增强子识别和验证提供额外的数据集来扩展管道的效用18,19。
在 TBX5 REPFIX 中,基于 H3K4me1 的分析揭示了 22 个假定的增强子,作为进一步研究的起点。随后的 Hi-C 分析显示,基于先前组蛋白修饰标记的 22 个候选增强子中有 21 个与心肌细胞中的 TBX5 启动子相互作用(图 5)。这支持了基于组蛋白修饰标记的方法在预测感兴趣区域方面的可靠性。
我们选择不优先考虑该方案中跨物种的序列保守,尽管这是识别增强子的通用标准。因为它之前被证明对组织或物种特异性增强子的效果较差,其中许多增强子在进化过程中并不强烈保守26。鉴于此,我们选择专注于更直接指示功能增强子活性的基于染色质的标记物。然而,序列守恒在特定情况下可能仍然具有价值,例如研究具有进化意义的增强子。在这种情况下,对于对保守监管元素感兴趣的用户,可以将其添加为可选步骤。
该方案被证明可有效从 EnsEMBL 网站检索 TBX5 基因的 21 个心脏特异性增强子,而现有平台 VISTA Cardiac Enhancers Browser 之前无法识别该增强子。尽管此处描述的方案无法检索到其中一种增强子 hs498,这表明检测可能存在限制,但该方法发现了一些 VISTA 心脏增强子浏览器未检测到的增强子。然而,对检索到的增强子的进一步验证是必要的,因为与精选的 VISTA 数据库相比,该协议产生了更多的推定区域。较高的计数会增加假阳性的风险,并且预测的增强子数量越多并不一定表示特异性或功能相关性更高。结合额外的实验数据集或功能分析,如基因表达分析、CRISPR 扰动或报告基因分析,对于确认 Smemo 等人所进行的这些候选基因的生物学有效性至关重要。S 研究23.
与三个实验验证的增强子的交叉比较显示与预测区域部分重叠(图 5)。预测的增强子 “e1” 比经过验证的增强子 2 定位更广泛,而 “e18” 与增强子 16 部分重叠(图 5)。这些结果表明,这种方法成功地识别了具有已知调节活性的区域,尽管预测的增强子的更广泛跨度可能反映了增强子边界的灵活性。增强子通常作为模块化元件发挥作用,其活性可能取决于染色质环境、细胞类型和发育时间 2,27。因此,预测区域可能包括核心活性位点和有助于调节功能的相邻序列,尽管需要实验验证以确定这些更广泛预测的哪些部分在组织特异性或发育环境中具有功能活性。虽然 VISTA 心脏增强子浏览器在定义范围内确定了四个区域,但只有一个增强子 mm370 与经过实验验证的增强子 Enhancer 16 部分重叠,特别是在 H3K4me1 标记的区域(图 5)23。mm370 的剩余部分与增强子 16 不重叠,可能表明增强子28,29 的非功能或无活性亚区。
作者报告没有利益冲突。
我们感谢 Ramialison 实验室(转录组学和生物信息学,reNEW 生物信息学中心)的成员进行有益的讨论。MR 和 HTN 由 NHMRC 创意资助 (APP1180905) 提供支持。我们感谢 Richard Saffery 的支持。MR 由心脏基金会未来领袖奖学金 (107328) 资助。澳大利亚政府国家卫生和医学研究委员会独立研究所基础设施支持计划为默多克儿童研究所提供了额外的基础设施资金。澳大利亚再生医学研究所得到了维多利亚州政府和澳大利亚政府的资助。诺和诺德基金会干细胞医学中心得到了诺和诺德基金会赠款 (NNF21CC0073729) 的支持。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 计算机工作站 | 不适用 | 不适用 | 支持网络浏览器的计算机,Windows/Mac/Linux作系统 |
| 4DN 数据门户 | 4DN 数据门户 | https://data.4dnucleome.org/ | |
| 星系 | 星系 | https://usegalaxy.org/published/history?id=aff5db4e07064445 | |
| Github的 | Github的 | https://github.com/Ramialison-Lab/EnhancerWorkflow | |
| VISTA | VISTA 心脏增强器浏览器 | https://portal.nersc.gov/dna/RD/heart/ |
申请许可以重复使用本 JoVE 文章的文本或图表
申请许可