Name: 用于选择基因和组织特异性增强子的基于 Web 的工作流程 (视频) | JoVE Journal
Uploaded: 2025-07-18T13:10:10.000+00:00
Duration: 8 min 12 s
Description: 学习使用基于浏览器的无代码工作流程识别组织特异性基因增强子。该方案利用公共组蛋白标记和Hi-C数据，以TBX5为案例，帮助生物学家无需编程技能即可进行增强子发现。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

感谢Ramialison实验室（转录组学与生物信息学，reNEW生物信息学中心）成员的有益讨论。MR和HTN的研究由澳大利亚国家健康与医学研究委员会（NHMRC） Ideas项目基金（APP1180905）资助。感谢Richard Saffery提供的支持。MR由心脏基金会未来领袖奖学金（107328）资助。默多克儿童研究所的额外基础设施经费由澳大利亚政府国家健康与医学研究委员会独立研究机构基础设施支持计划提供。澳大利亚再生医学研究所由维多利亚州政府和澳大利亚政府的资助项目支持。诺和诺德基金会干细胞医学中心由诺和诺德基金会基金（NNF21CC0073729）支持。

注意：分步操作指南可在 https://github.com/Ramialison-Lab/EnhancerWorkflow 获取。本方案中使用的数据总结于表1和表2中。故障排除方法参见补充文件1。

1. 定位目的基因（图1）

<ol>
	<li>打开 EnsEMBL 基因组浏览器（https://www.ensembl.org）。</li>
	<li>选择与物种及版本匹配的适当基因组组装。</li>
	<li>在搜索框中输入目的基因（GoI），然后点击 Go。</li>
	<li>选择指向相应 EnsEMBL 基因 ID 的链接。</li>
	<li>点击摘要部分下方的 Region in details 标签页查看器链接，以浏览目的基因周围的区域。</li>
</ol>

2. 定义增强子检测区域（图2）

<ol>
	<li>&#160;确定目的基因（GoI）的两个相邻基因以定义检测区域。使用基因图例（Gene Legend）轨迹查找这些基因，它们被注释为与视觉相关的元件 合并的 EnsEMBL/Havana 在 来自GENCODE的基础基因注释 追踪。通过以下方式确定目的基因（GoI）的方向性 &#60; 和 &#62; 在GENCODE轨道上相应基因名称处的标记</li>
	<li>点击并拖动以选中这些基因之间的基因间区域，然后点击 跳转到区域 在弹出框中可视化所选区域。通过重复此步骤，可随时重新定义感兴趣区域。</li>
	<li>（可选）通过选择自定义显示方式 添加/移除轨道 位于轨迹查看器顶部。使用轨迹查看器上方的缩放/导航控件调节缩放级别。 &#160;</li>
</ol>

3. 组蛋白修饰分析（图3）

<ol>
	<li>在区域详情标签页查看器的侧边栏中，选择配置此页面。</li>
	<li>在配置区域图像标签页显示的侧边栏中，于调控下拉菜单下选择按细胞/组织的活性。</li>
	<li>使用细胞/组织搜索栏搜索并选择感兴趣的组织。或者，使用搜索栏下方的字母导航栏查找感兴趣的组织。</li>
	<li>选择位于细胞/组织旁边的实验标签页。</li>
	<li>选择H3K4me1和H3K27ac作为增强子的标记，选择H3K4me3作为启动子的标记。</li>
	<li>选择配置轨道显示。</li>
	<li>选择查看轨道，以可视化增强子检测区域内被H3K4me1标记的区域以及目的基因（GoI）上游被H3K4me3标记的区域。</li>
	<li>对于增强子候选区域，通过点击新添加的H3K4me1轨道中的彩色可视化/方框元素，获取在定义的检测区域内被H3K4me1标记的基因组区域的坐标。这将弹出一个包含该元素基因组位置（以碱基对 bp 表示）信息的&#34;Hists &amp;#38; Pols&#34;窗口。
	<ol>
		<li>通过选择H3K4me1和H3K27ac的彩色可视化/方框元素重叠的H3K4me1区域，获取活性增强子候选区域。缩放级别的精细程度可能影响该轨道中显示的增强子候选数量。</li>
		<li>或者，通过在轨道上点击并拖动以框选H3K4me1/H3K27ac轨道下方的图峰，手动定义每个基因组特征的兴趣区域。将基因组位置坐标复制到文本文件中，并保存为.bed格式。</li>
	</ol>
	</li>
	<li>同样地，对于启动子区域，使用H3K4me3轨道重复步骤3.8，重点关注目的基因（GoI）上游的区域。</li>
</ol>

4. 染色质构象捕获（Hi-C）分析（图4）

<ol>
	<li>访问 4DN 数据门户（https://data.4dnucleome.org/）。在 4DN 首页的主堆叠条形图中（图 5），确保 Y 轴选择为实验集，X 轴选择为实验类型，并按生物体进行分组。</li>
	<li>在主条形图的 X 轴上找到原位 Hi-C 条形。点击该条形中对应于人类实验集的部分，然后点击弹出窗口中的浏览按钮。使用左侧面板筛选相关数据集。</li>
	<li>点击目标组织相关生物样本的标题列中的链接。</li>
	<li>在已处理文件选项卡中点击探索数据，以更详细地查看 Hi-C 数据集。&#160;</li>
	<li>输入在目标组织中鉴定出的启动子坐标（来自步骤 3.9），通过右键点击热图并在相应位置标记水平线（图 4）。添加的线条可确保在热图中直观追踪启动子区域。若需删除误操作的线条，可右键点击该线条，选择水平/垂直标线，然后选择关闭系列。</li>
	<li>输入所有经实验验证的对照增强子23的基因组坐标，基于其最小值计算相互作用阈值。这些对照增强子此前已被证实与启动子区域存在相互作用，可作为定义最小显著相互作用的基准。</li>
	<li>根据矩阵右侧颜色图例中的最低非零相互作用得分，利用对照增强子确定启动子-增强子相互作用的阈值。</li>
	<li>输入所有与 H3K4me1 相关区域的基因组坐标（来自步骤 3.8），并在 Hi-C 热图上进行垂直标记。这将形成交叉核对模式，便于检查 H3K4me1 标记的增强子（垂直方向）与启动子区域（水平方向）之间的交点是否存在显著相互作用。
	<ol>
		<li>随后，通过将 H3K4me3 标记区域的相互作用得分与相互作用阈值（步骤 4.6）进行比较，过滤掉相互作用较弱的区域。&#160;</li>
		<li>在 4DN 门户热图中选择相互作用频率高于定义阈值的基因组坐标区域。这些区域在热图中表现为信号更集中（颜色更深），以 BED 格式保存。</li>
	</ol>
	</li>
</ol>

<ol><li>Barral, A., Déjardin, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+chromatin+signatures+of+enhancers+and+their+dynamic+regulation%5BTitle%5D)+AND+%22Nucleus%22%5BJournal%5D)">The chromatin signatures of enhancers and their dynamic regulation.</a> Nucleus. 14 (1), 2160551(2023).</li>
<li>Pennacchio, L. A., Bickmore, W., Dean, A., Nobrega, M. A., Bejerano, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhancers%3A+five+essential+questions%5BTitle%5D)+AND+%22Nat+Rev+Genet%22%5BJournal%5D)">Enhancers: five essential questions.</a> Nat Rev Genet. 14 (4), 288-295 (2013).</li>
<li>Kleinjan, D. A., van Heyningen, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Long-range+control+of+gene+expression%3A+emerging+mechanisms+and+disruption+in+disease%5BTitle%5D)+AND+%22Am+J+Hum+Genet%22%5BJournal%5D)">Long-range control of gene expression: emerging mechanisms and disruption in disease.</a> Am J Hum Genet. 76 (1), 8-32 (2005).</li>
<li>Kvon, E. Z., Waymack, R., Gad, M., Wunderlich, Z. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhancer+redundancy+in+development+and+disease%5BTitle%5D)+AND+%22Nat+Rev+Genet%22%5BJournal%5D)">Enhancer redundancy in development and disease.</a> Nat Rev Genet. 22 (5), 324-336 (2021).</li>
<li>Smith, E., Shilatifard, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhancer+biology+and+enhanceropathies%5BTitle%5D)+AND+%22Nat+Struct+Mol+Biol%22%5BJournal%5D)">Enhancer biology and enhanceropathies.</a> Nat Struct Mol Biol. 21 (3), 210-219 (2014).</li>
<li>Sur, I., Taipale, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+role+of+enhancers+in+cancer%5BTitle%5D)+AND+%22Nat+Rev+Cancer%22%5BJournal%5D)">The role of enhancers in cancer.</a> Nat Rev Cancer. 16 (8), 483-493 (2016).</li>
<li>Panigrahi, A., O'Malley, B. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mechanisms+of+enhancer+action%3A+the+known+and+the+unknown%5BTitle%5D)+AND+%22Genome+Biol%22%5BJournal%5D)">Mechanisms of enhancer action: the known and the unknown.</a> Genome Biol. 22 (1), 108(2021).</li>
<li>Bulger, M., Groudine, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+and+mechanistic+diversity+of+distal+transcription+enhancers%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Functional and mechanistic diversity of distal transcription enhancers.</a> Cell. 144 (3), 327-339 (2011).</li>
<li>Krivega, I., Dean, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhancer+and+promoter+interactions-long+distance+calls%5BTitle%5D)+AND+%22Curr+Opin+Genet+Dev%22%5BJournal%5D)">Enhancer and promoter interactions-long distance calls.</a> Curr Opin Genet Dev. 22 (2), 79-85 (2012).</li>
<li>Zhang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chromatin+connectivity+maps+reveal+dynamic+promoter-enhancer+long-range+associations%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Chromatin connectivity maps reveal dynamic promoter-enhancer long-range associations.</a> Nature. 504 (7479), 306-310 (2013).</li>
<li>Fang, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mapping+of+long-range+chromatin+interactions+by+proximity+ligation-assisted+ChIP-seq%5BTitle%5D)+AND+%22Cell+Res%22%5BJournal%5D)">Mapping of long-range chromatin interactions by proximity ligation-assisted ChIP-seq.</a> Cell Res. 26 (12), 1345-1348 (2016).</li>
<li>Park, P. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ChIP-seq%3A+advantages+and+challenges+of+a+maturing+technology%5BTitle%5D)+AND+%22Nat+Rev+Genet%22%5BJournal%5D)">ChIP-seq: advantages and challenges of a maturing technology.</a> Nat Rev Genet. 10 (10), 669-680 (2009).</li>
<li>Thurman, R. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+accessible+chromatin+landscape+of+the+human+genome%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">The accessible chromatin landscape of the human genome.</a> Nature. 489 (7414), 75-82 (2012).</li>
<li>Heintzman, N. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distinct+and+predictive+chromatin+signatures+of+transcriptional+promoters+and+enhancers+in+the+human+genome%5BTitle%5D)+AND+%22Nat+Genet%22%5BJournal%5D)">Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome.</a> Nat Genet. 39 (3), 311-318 (2007).</li>
<li>Visel, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ChIP-seq+accurately+predicts+tissue-specific+activity+of+enhancers%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">ChIP-seq accurately predicts tissue-specific activity of enhancers.</a> Nature. 457 (7231), 854-858 (2009).</li>
<li>Belton, J. -M., McCord, R. P., Gibcus, J. H., Naumova, N., Zhan, Y., Dekker, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Hi-C%3A+a+comprehensive+technique+to+capture+the+conformation+of+genomes%5BTitle%5D)+AND+%22Methods%22%5BJournal%5D)">Hi-C: a comprehensive technique to capture the conformation of genomes.</a> Methods. 58 (3), 268-276 (2012).</li>
<li>de Wit, E., de Laat, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+decade+of+3C+technologies%3A+insights+into+nuclear+organization%5BTitle%5D)+AND+%22Genes+Dev%22%5BJournal%5D)">A decade of 3C technologies: insights into nuclear organization.</a> Genes Dev. 26 (1), 11-24 (2012).</li>
<li>Gao, T., Qian, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((EnhancerAtlas+2.0%3A+an+updated+resource+with+enhancer+annotation+in+586+tissue%2Fcell+types+across+nine+species%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">EnhancerAtlas 2.0: an updated resource with enhancer annotation in 586 tissue/cell types across nine species.</a> Nucleic Acids Res. 47 (D1), D117-D123 (2019).</li>
<li>Erwin, G. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Integrating+diverse+datasets+improves+developmental+enhancer+prediction%5BTitle%5D)+AND+%22PLoS+Comput+Biol%22%5BJournal%5D)">Integrating diverse datasets improves developmental enhancer prediction.</a> PLoS Comput Biol. 10 (6), e1003677(2014).</li>
<li>Heintzman, N. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Histone+modifications+at+human+enhancers+reflect+global+cell-type-specific+gene+expression%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Histone modifications at human enhancers reflect global cell-type-specific gene expression.</a> Nature. 459 (7243), 108-112 (2009).</li>
<li>Spicuglia, S., Vanhille, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chromatin+signatures+of+active+enhancers%5BTitle%5D)+AND+%22Nucleus%22%5BJournal%5D)">Chromatin signatures of active enhancers.</a> Nucleus. 3 (2), 126-131 (2012).</li>
<li>Borsari, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enhancers+with+tissue-specific+activity+are+enriched+in+intronic+regions%5BTitle%5D)+AND+%22Genome+Res%22%5BJournal%5D)">Enhancers with tissue-specific activity are enriched in intronic regions.</a> Genome Res. 31 (8), 1325-1336 (2021).</li>
<li>Smemo, S., Campos, L. C., Moskowitz, I. P., Krieger, J. E., Pereira, A. C., Nobrega, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Regulatory+variation+in+a+TBX5+enhancer+leads+to+isolated+congenital+heart+disease%5BTitle%5D)+AND+%22Hum+Mol+Genet%22%5BJournal%5D)">Regulatory variation in a TBX5 enhancer leads to isolated congenital heart disease.</a> Hum Mol Genet. 21 (14), 3255-3263 (2012).</li>
<li>Boogerd, C. J., Evans, S. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((TBX5+and+NuRD+divide+the+heart%5BTitle%5D)+AND+%22Dev+Cell%22%5BJournal%5D)">TBX5 and NuRD divide the heart.</a> Dev Cell. 36 (3), 242-244 (2016).</li>
<li>Dickel, D. E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genome-wide+compendium+and+functional+assessment+of+in+vivo+heart+enhancers%5BTitle%5D)+AND+%22Nat+Commun%22%5BJournal%5D)">Genome-wide compendium and functional assessment of in vivo heart enhancers.</a> Nat Commun. 7 (1), 12923(2016).</li>
<li>Pennacchio, L. A., Visel, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Limits+of+sequence+and+functional+conservation%5BTitle%5D)+AND+%22Nat+Genet%22%5BJournal%5D)">Limits of sequence and functional conservation.</a> Nat Genet. 42 (7), 557-558 (2010).</li>
<li>Local, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Identification+of+H3K4me1-associated+proteins+at+mammalian+enhancers%5BTitle%5D)+AND+%22Nat+Genet%22%5BJournal%5D)">Identification of H3K4me1-associated proteins at mammalian enhancers.</a> Nat Genet. 50 (1), 73-82 (2018).</li>
<li>Creyghton, M. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Histone+H3K27ac+separates+active+from+poised+enhancers+and+predicts+developmental+state%5BTitle%5D)+AND+%22Proc+Natl+Acad+Sci+U+S+A%22%5BJournal%5D)">Histone H3K27ac separates active from poised enhancers and predicts developmental state.</a> Proc Natl Acad Sci U S A. 107 (50), 21931-21936 (2010).</li>
<li>Visel, A., Minovitsky, S., Dubchak, I., Pennacchio, L. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((VISTA+Enhancer+Browser-a+database+of+tissue-specific+human+enhancers%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">VISTA Enhancer Browser-a database of tissue-specific human enhancers.</a> Nucleic Acids Res. 35 (Database), D88-D92 (2007).</li>
</ol>

作者声明无竞争利益。

本文所述的基于网络的方案旨在实现增强子的检索，而非增强子的预测。通过利用公开可用的组织特异性数据集、组蛋白修饰（H3K4me1 和 H3K27ac）以及 Hi-C 相互作用数据，该方法可缩小与目的基因（GoI）相关的潜在增强子范围。与依赖机器学习或基于序列模型的计算预测工具不同，我们的方法仅基于实验数据进行增强子的检索。EnsEMBL 基因组浏览器的故障排除说明详见补充文件 1。
该工作流程整合了组蛋白修饰和染色质相互作用数据等关键要素，类似于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网站检索到21个心脏特异性增强子，这些增强子针对TBX5基因，而此前在现有平台VISTA Cardiac Enhancers Browser中未能被识别。尽管本方案未能检索到其中一个增强子hs498，提示其检测能力可能存在一定局限性，但该方法仍发现了一些VISTA Cardiac Enhancers Browser未能检测到的增强子。然而，对所获得的增强子仍需进一步验证，因为该方案产生的候选区域数量较经人工审编的VISTA数据库更多。较高的数量增加了假阳性风险，而预测增强子数量的增加并不一定意味着特异性或功能相关性的提升。正如Smemoet al.的研究中所采用的方法23，整合额外的实验数据集或功能检测手段，如基因表达分析、CRISPR扰动实验或报告基因检测，对于确认这些候选增强子的生物学有效性至关重要。
与三个经过实验验证的增强子进行交叉比较，结果显示预测区域与这些增强子存在部分重叠（图5）。预测的增强子&#34;e1&#34;比已验证的增强子2定位范围更广，而&#34;e18&#34;与增强子16存在部分重叠（图5）。这些结果表明，该方法能够成功识别具有已知调控活性的区域，尽管预测增强子范围较广，可能反映了增强子边界的灵活性。增强子通常作为模块化元件发挥作用，其活性可能依赖于染色质环境、细胞类型以及发育时机2,27。因此，预测区域可能包含核心活性位点以及对调控功能有贡献的邻近序列，但仍需通过实验验证来确定这些较宽泛预测区域中哪些部分在特定组织或发育背景下具有功能性活性。尽管VISTA心脏增强子浏览器在定义范围内识别出四个区域，但仅有一个增强子mm370与实验验证的增强子16存在部分重叠，且该重叠仅发生在H3K4me1标记的区域（图5）23。增强子mm370中未与增强子16重叠的其余部分，可能代表该增强子中的非功能性或失活亚区28,29。

增强子是一类调控基因转录、发育及细胞分化的非编码DNA区域1,2。增强子的突变可导致多种疾病，包括发育障碍、癌症及其他遗传性疾病3,4,5,6。因此，理解增强子对于阐明基因表达、突变及疾病机制至关重要。
为了理解增强子如何与其靶基因相互作用，确定其在基因组中的位置至关重要。然而，识别增强子的位置并不总是直接的，因为增强子既可以位于转录起始位点（TSS）附近，也可以远在数十至数百千碱基之外2,7,8,9。
尽管增强子的基因组位置难以预测，但它们具有独特的生物化学和结构特征，因此可以进行系统性追踪。通常，增强子富集于基因间区和内含子区域，仅有少数位于外显子内2,8。它们通常带有特定的组蛋白修饰和转录因子结合位点，这些特征决定了其调控功能，并在不同发育阶段和组织中决定其时空活性10,11。
ChIP-seq 用于鉴定转录因子结合位点（TFBSs）以及组蛋白修饰标记，例如增强子的标志性修饰 H3K4me1、活性增强子标记 H3K27ac 和富集于启动子区域的 H3K4me3 标记1,12,13,14,15。染色质构象捕获技术（3C）及其衍生技术，如 4C、5C、Hi-C 和 ChIA-PET，可用于绘制远端基因组区域之间的物理相互作用图谱。3C 技术针对特定组织中的特定相互作用，而 Hi-C 则可提供跨细胞类型的全基因组架构信息16,17。
除了现有方法外，研究人员还开发了专门用于表征增强子的方法，包括整合型增强子数据库，如 EnhancerAtlas 或 EnhancerFinder18,19。然而，这些工具通常要求研究者整合多个数据集，以在多种组织中研究多个增强子，这对于缺乏生物信息学和数据挖掘经验的生物学家而言可能难以应对。
本文介绍了一种完全基于现有网络工具的用户友好型增强子筛选方案。该方案允许研究人员查询目标基因（GoI），并获取相应的增强子。本方案根据一组特定标准筛选增强子：组蛋白修饰、染色质相互作用以及组织特异性1,12,13,14,15,16,17,20,21。位于内含子内的增强子相比基因间区增强子（位于基因间基因组区域）更可能表现出组织特异性活性22。为确保尽可能全面地覆盖潜在的活性增强子，我们将搜索范围定义在两个相邻目标基因之间，以提高捕获基因体外调控元件的可能性。我们利用增强子特异性的表观遗传标志H3K4me1和活性增强子标志H3K27ac来列出增强子候选区域。随后，基于Hi-C数据对这些候选区域进行进一步筛选，仅保留与相应启动子存在物理相互作用的增强子。本方案旨在指导生物学家仅使用公开可用的基于网络的工具完成增强子鉴定过程。通过整合表观遗传学和染色质相互作用数据，本文所述方法为提出有关潜在增强子的假设提供了实用框架，可用于后续实验验证。

<table><tbody><tr><td>计算机工作站</td><td>N/A</td><td>N/A</td><td>具备网页浏览功能的计算机，支持 Windows/Mac/Linux 操作系统</td></tr><tr><td>4DN 数据门户</td><td>4DN 数据门户</td><td>https://data.4dnucleome.org/</td><td></td></tr><tr><td>Galaxy</td><td>Galaxy</td><td>https://usegalaxy.org/published/history?id=aff5db4e07064445</td><td></td></tr><tr><td>Github</td><td>Github</td><td>https://github.com/Ramialison-Lab/EnhancerWorkflow</td><td></td></tr><tr><td>VISTA&amp;nbsp;</td><td>VISTA 心脏增强子浏览器</td><td>https://portal.nersc.gov/dna/RD/heart/</td><td></td></tr></tbody></table>

一种用于筛选基因和组织特异性增强子的基于网络的分析流程

A Web-Based Workflow for Selecting Gene- and Tissue-Specific Enhancers

增强子是调控基因表达的DNA区域。增强子内部的突变可能导致基因调控异常，从而引发疾病。因此，鉴定在特定组织中调控基因活性的增强子，对于理解疾病的遗传基础至关重要。然而，由于增强子不编码蛋白质，其识别较为困难。尽管目前已存在大量增强子数据库及识别工具，但这些工具的复杂性可能为生物学家带来使用上的挑战。为了帮助生物学家更便捷地利用这些资源，我们提供了一套面向生物学家的友好型实验流程（https://github.com/Ramialison-Lab/EnhancerWorkflow），该流程利用现有的基于网页的基因组学数据，如H3K4me1和H3K27ac组蛋白修饰标记以及染色质构象分析（Hi-C）数据，用于发现与目标组织中处于活性状态的感兴趣基因（GoI）相关的增强子。本流程完全基于网页操作，终端用户无需具备编程技能。我们通过鉴定调控&#160;TBX5（一个对心脏发育至关重要的基因）的候选增强子，验证了该方法的有效性。该流程有助于识别在左心室中与该基因相关的增强子。

为了说明本方案的应用，我们研究了人类心脏中的 TBX5 基因，利用包含 H3K4me1、H3K27ac 和 Hi-C 数据的综合工作流程，探索与 TBX5 相关的增强子。TBX5 是一个参与肢体和心脏发育的基因，包括心脏四个腔室的形成以及室间隔的分隔24。该基因的突变是导致霍尔特-奥拉姆综合征（Holt-Oram syndrome, HOS）的主要原因，该综合征可引起肢体异常和先天性心脏病（congenital heart disease, CHD），包括室间隔缺损24。与 TBX5 相关的心脏增强子发生突变可能对先天性心脏病产生关键影响24。先前一项研究在人类心脏特异性组织中发现了三个已知的 TBX5 增强子——分别为“Enhancer 2”、“Enhancer 9”和“Enhancer 16”（补充文件 2），并在转基因小鼠中证实其具有相似的表型23。
我们在人类中位于RBM19和TBX3这两个分别位于TBX5下游和上游的邻近基因之间，研究了H3K4me1和H3K27ac富集的区域，以获取TBX5位点的潜在增强子（图1和图2）。为了鉴定心脏特异性增强子，我们选择了心肌细胞。潜在的TBX5心脏增强子区域以坐标形式（chr12：起始-终止）提取，共鉴定出22个同时与H3K4me1和H3K27ac相关的区域（图3和补充文件3）。从EnsEMBL基因组数据库中提取潜在的TBX5心脏增强子区域，并与4DNucleome数据库中的Hi-C数据进行交叉比对（图4），以评估这些潜在增强子与TBX5心脏启动子之间可能存在的相互作用。根据本文所述的实验流程，在心肌细胞中确认了22个基因组区域中有21个与TBX5启动子（chr12: 114400143-114410103）存在相互作用（补充文件4）。有一个区域未检测到与启动子的物理相互作用（图4，步骤4.8）。最后，我们将本方案鉴定出的结果与已通过生物学验证的增强子以及当前心脏增强子的金标准数据库VISTA Cardiac Enhancers Browser进行比较，发现了一些该数据库尚未收录的额外增强子25。
我们对21进行了交叉比较 TBX5 通过本方案从现有数据库中检索到的增强子。我们检索到了4 TBX5 来自VISTA心脏增强子浏览器的增强子（补充文件 5)25在VISTA鉴定出的4个心脏增强子中，有3个增强子hs2329、mm1282和m370与本基于网络的增强子检测方案所鉴定出的区域重叠图5）。每个预测的增强子还与 Smemo 等人先前通过实验验证的增强子共享基因组区域。23增强子2（chr12:114025907-114026275，GRCh38）和增强子16（chr12:114415466-114420433，GRCh38），但它们与增强子9（chr12:114263402-114266886，GRCh38）无重叠。VISTA鉴定出的增强子之一， hs498 未与本方案预测的任何增强子或 Smemo 等人实验验证的增强子重叠23 (图5)，尽管该区域与H3K4me1标记显示出部分重叠图5）。类似地，增强子9并未与该分析流程预测的增强子区域重叠，但与H3K4me1标记相关联（图5).
<img alt="TBX5基因的Ensembl数据库搜索工作流程；基因组数据检索示意图。" class="xfigimg" src="/files/ftp_upload/66840/66840fig1.jpg" data-altupdatedat="1762425657000" data-alt="图1"> 
图1：在Ensembl基因组浏览器中定位目的基因（GoI）的分步指南。用户首先打开Ensembl主页（1.1），选择物种（人），并在搜索栏中输入目标基因（1.2–1.3）。在结果列表中选择正确的基因ID（1.4），进入该基因的摘要页面。随后，用户点击“区域详情”超链接（1.5），以可视化目的基因周围区域的基因组结构，包括邻近元件和调控特征。<a href="https://www.jove.com/files/ftp_upload/66840/66840fig1large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="Ensembl基因组浏览器界面；基因比对图；染色体区域分析。" class="xfigimg" src="/files/ftp_upload/66840/66840fig2.jpg" data-altupdatedat="1762425657000" data-alt="图2"> 
图2：利用Ensembl基因组浏览器定义目标基因（GoI）周围的增强子检测区域。为定义增强子检测区域，需使用GENCODE基础基因注释轨迹，识别位于目标基因两侧的相邻基因，这些基因以深黄色方块表示，并标注有合并的Ensembl/Havana注释。每个基因的转录方向由其基因名称旁的箭头符号（&lt; 或 &gt;）指示（2.1）。为选择相邻基因之间的基因间区域，点击并拖动选中目标区域，然后在弹出框中选择跳转至该区域以放大显示（2.2）。为添加调控或增强子相关注释，点击添加/移除轨迹（2.3）。<a href="https://www.jove.com/files/ftp_upload/66840/66840fig2large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="基因组浏览器界面步骤，突出显示设置、轨迹选项和基因活性可视化。" class="xfigimg" src="/files/ftp_upload/66840/66840fig3v2.jpg" data-altupdatedat="1762425657000" data-alt="图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 &amp; Pols”弹出窗口包含该区域以碱基对表示的基因组坐标（chr:start-end），可复制并保存用于下游分析。点击轨迹中的彩色元素后，将弹出“Hists &amp; Pols”窗口，其中包含该区域的基因组坐标（例如，启动子区域为 chr12:11443450-114451611），可复制并保存用于下游分析（3.8）。同样，为提取候选增强子，应优先选择 H3K4me1 和 H3K27ac 峰重叠的区域，如各轨迹中峰与方框的垂直对齐所示（3.9）。可通过点击其方框或手动点击并拖拽跨越对齐的峰来直接选择重叠区域（例如，活性候选区域为 chr12:114400143-114410103）。弹出窗口中显示的坐标应以 BED 格式保存，用于后续验证或可视化分析。 <a href="https://www.jove.com/files/ftp_upload/66840/66840fig3large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="4D核组数据门户界面，包含交互式图表、数据可视化和分析功能。" class="xfigimg" src="/files/ftp_upload/66840/66840fig4v2.jpg" data-altupdatedat="1762425657000" data-alt="图 4"> 
图 4：利用来自4D核组数据门户的Hi-C热图可视化启动子-增强子染色质相互作用。 4D核组数据门户主页显示了一个堆叠条形图，按生物体总结了可用的实验类型。“人类样本”的“原位 Hi-C”数据集通过点击条形图的相应部分被选中（4.1）。系统显示相关数据集的筛选列表；选择一个来源于H9细胞分化为心肌母细胞的Hi-C数据集（4.2）。通过“探索数据”按钮（4.4）在HiGlass浏览器中打开所选数据集（4.3）。在坐标输入框中输入感兴趣的基因组区域（4.5），接触矩阵以颜色标度的热图形式呈现。较深的颜色（深红至黑色）表示更强的染色质接触频率，而较浅的颜色（白色至橙色）代表较弱的相互作用。在启动子坐标处放置一条水平线，并在三个经实验验证的对照增强子位置绘制垂直线（4.6）。这些交点用于定义严格的相互作用阈值，该阈值由启动子-增强子接触中最强的可见信号（最深颜色）设定（4.7）。在候选H3K27ac和H3K4me1标记的增强子位置（来自步骤3.8）绘制额外的垂直线。保留其启动子-增强子交点颜色等于或深于阈值的候选区域，而信号较弱（颜色较浅的方块）的候选区域则被排除（4.8）。保留的坐标手动提取并以BED格式保存，用于后续分析。（a. 增强子2，b. 增强子9，和c. 增强子16） <a href="https://www.jove.com/files/ftp_upload/66840/66840fig4large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
<img alt="染色体12基因组序列分析图，包含详细区域注释和基因位置。" class="xfigimg" src="/files/ftp_upload/66840/66840fig5.jpg" data-altupdatedat="1762425657000" data-alt="图5"> 
图5：预测的TBX5增强子与VISTA验证的心脏增强子及对照增强子的基因组浏览器视图比较。基因组浏览器截图展示了增强子搜索范围（步骤2），将基于网络协议检索到的预测增强子（底部）与VISTA验证的增强子（顶部）以及实验验证的对照增强子（中部）进行比较。主图显示了完整的基因组位点，包含已注释的调控元件，如心肌细胞特异性的H3K4me1（黄色）、H3K27ac（蓝色）和H3K4me3峰（橙色）。三个放大图展示了由本协议检索到的增强子、VISTA增强子和对照增强子之间的比对情况。与对照增强子的重叠区域用红色框标出。各子区域的坐标显示在下方的浏览器面板中。<a href="https://www.jove.com/files/ftp_upload/66840/66840fig5large.jpg" target="_blank">请点击此处查看该图的高清版本。</a>
表1：本研究中使用数据。 <a href="https://www.jove.com/files/ftp_upload/66840/Table%201.xlsx" target="_blank">请点击此处下载本表格。</a>
表2：本研究中使用的基于网络的工具。 <a href="https://www.jove.com/files/ftp_upload/66840/Table%202.xlsx" target="_blank">请点击此处下载该表格。</a>
补充文件1：EnsEMBL基因组浏览器故障排除指南。<a href="https://www.jove.com/files/ftp_upload/66840/Supplementary%20File%201.docx" target="_blank">请点击此处下载该文件</a>。
补充文件 2：以 GRCh38 格式提供的 BED 文件，包含经实验验证的 TBX5 心脏调控增强子23。 <a href="https://www.jove.com/files/ftp_upload/66840/Supplementary%20File%202.bed" target="_blank">请点击此处下载该文件。</a>
补充文件 3：以 GRCh38 格式提供的 BED 文件，TBX5 心脏增强子由 STEP3 从 EnsEMBL 中检索获得。 <a href="https://www.jove.com/files/ftp_upload/66840/Supplementary%20File%203.bed" target="_blank">请点击此处下载该文件。</a>
补充文件4：以GRCh38格式提供的BED文件，TBX5心脏增强子由STEP4从EnsEMBL中检索获得。 <a href="https://www.jove.com/files/ftp_upload/66840/Supplementary%20File%204.bed" target="_blank">请点击此处下载该文件。</a>
补充文件5：以GRCh38格式提供的BED文件，TBX5心脏增强子数据来自VISTA心脏增强子浏览器25。 <a href="https://www.jove.com/files/ftp_upload/66840/Supplementary%20File%205.bed" target="_blank">请点击此处下载该文件。</a>

我们提供一种无需编程的生物学家工作流程，仅使用基于浏览器的工具即可识别组织特异性基因增强子。本方案利用公共的H3K4me1/H3K27ac组蛋白修饰标记和Hi-C数据，使不具备编程经验的研究人员也能够获取、分析并识别与其目标基因相关的潜在调控元件。

增强子是调控基因表达的DNA区域。增强子内部的突变可能导致基因调控异常，从而引发疾病。因此，……

我们致力于破译非编码基因组，以理解基因如何在正确的时间和正确的位点被调控表达。我们利用基于网络的基因组学工具开发了易于使用的实验方案，使所有生物学家都能便捷地进行深入发现。多模态数据的广泛可及性有助于缩小增强子位置的范围，而通过开源网络界面可轻松访问这些数据。这使得研究人员无需具备编程专业知识即可进行全面的增强子鉴定。本方案为任何没有增强子生物学背景的研究人员提供了基础，使其能够开始利用公开可用的数据集，检索目标基因的增强子。我们以在心脏生物学中具有重要作用的 TBX5 为例进行案例研究。我们的工作流程鉴定了21个增强子，可能为心脏发育机制及先天性心脏病的研究提供线索。我们将结合包括空间基因组学在内的新数据来源扩展本方案，同时开发更多便于生物学家使用的工具。首先，打开 Ensembl 基因组浏览器。选择与目标物种及版本匹配的基因组组装。在搜索框中输入目标基因，然后点击“Go”。在结果中，点击相应的 Ensembl 基因 ID。然后滚动至“Summary”部分，点击“Region in Detail”链接，以查看目标基因周围区域的可视化图谱。接下来，利用“Gene Legend”轨迹，搜索位于目标基因两侧的两个邻近基因。这些基因以视觉元素的形式呈现，代表了GENCODE轨道中基本基因注释所包含的Ensembl与Havana注释的合并结果。通过观察基因名称旁的大于号或小于号来确定目标基因的方向性。点击并拖动光标跨越这些基因之间的基因间区域，然后在弹出框中点击“Jump to Region”以查看选定的区域。通过点击轨道查看器顶部的“添加或移除轨道”来自定义显示。使用缩放和导航控件调整视图，以增强该区域的可视化效果。在“详细区域”标签查看器中，点击侧边栏中的“配置此页面”。在“配置区域图像”侧边栏中，于“调控”部分下选择“按细胞或组织划分的活性”。使用细胞或组织搜索栏查找并选择所需的组织，或使用搜索栏下方的字母索引。单击与“细胞或组织”选项相邻的“实验”标签页。选择 H3K4me1 和 H3K27ac 作为增强子的标志物，H3K4me3 作为启动子的标志物，然后点击“配置轨迹显示”。接着，选择“查看轨迹”，以可视化增强子检测区域内由 H3K4me1 标记的区域，以及目的基因上游由 H3K4me3 标记的区域。点击 H3K4me1 轨迹中的彩色可视化元素或方框，以获取在定义的检测区域内标记区域的基因组坐标。这将打开“Hists and Pols”弹窗，显示该元件以碱基对表示的基因组位置。或者，通过在轨迹图上点击并拖动，以框选H3K4me1或H3K27ac轨迹下图表峰所在的区域，从而手动定义每个基因组特征的兴趣区域。然后将基因组位置坐标复制到文本文件中，并将文件保存为BED格式。访问4DN数据门户。在主页上，确保将“实验集”选为堆叠柱状图的Y轴，“实验类型”选为X轴，并按生物体对图表进行分组。在X轴上找到原位Hi-C对应的柱状图部分，然后点击代表人类实验集的区段。在弹出窗口中，点击“浏览”按钮，选择心肌成肌细胞以及我们关注的与心脏器官发生相关的基因 TBX5。然后，点击相应生物样本“标题”列中的链接，该样本对应于目标组织。接着，在“已处理文件”选项卡中，点击“探索数据”以更详细地查看 Hi-C 数据集。输入在目标组织中鉴定出的启动子的坐标，然后右键点击热图以水平标记该区域。这些线条可实现启动子区域在热图上的可视化追踪。如有需要，可通过垂直拖动视图，直至搜索栏中Y坐标的上下边界与目标区域的边界对齐。要删除任何不需要的线条，右键单击该线条并选择“水平”或“垂直标尺”，然后点击关闭系列。输入所有经实验验证的对照增强子的坐标，以基于其最小相互作用值计算相互作用阈值。使用预设视图，在Ensembl提供的启动子区域旁边显示文献中的三个对照增强子。然后利用这些对照增强子，通过选择热图矩阵右侧颜色图例所示的最低非零相互作用得分，确定启动子-增强子的相互作用阈值。输入所有与H3K4me1相关的增强子区域的基因组坐标，并在Hi-C热图上垂直标记。通过将H3K4me3标记区域的相互作用得分与该阈值进行比较，过滤掉相互作用较弱的区域。选择在热图上显示为更密集或颜色更深、且相互作用频率高于此阈值的基因组坐标区域，并以 BED 格式保存。在 VISTA 心脏增强子浏览器中鉴定出的四个心脏增强子中有三个（HS2329、mm1282 和 mm370）与基于网络的增强子检测方案所预测的区域重叠。增强子二则同时与增强子检测方案预测的区域以及一个经过实验验证的增强子区域重叠。增强子16与预测的增强子区域以及先前的实验数据均显示出重叠。增强子Nine未与分析流程中预测的任何增强子区域重叠，但显示出部分H3K4me1信号富集。

a-web-based-workflow-for-selecting-gene-and-tissue-specific-enhancers

Research

科研

JoVE Journal

Genetics

遗传学

用于选择基因和组织特异性增强子的基于 Web 的工作流程

在本研究中，我们提供了一种实验流程，用于鉴定并实验验证活性增强子。...

视黄酸诱导胚胎干细胞分化过程中激活的基因调控元件的预测与验证

与DNA序列数据不同，表观基因组数据无法直接进行基于文本的搜索。本文介绍了使用GeNemo升级版的操作流程，GeNemo是一种基于网络的生物信息学工具，可用于在表观基因组数据中开展基于模式的相似性搜索，将用户数据与在线数据库（如DNA元件百科全书）中的数据进行比对。...

基于模式的表观基因组数据搜索：使用 GeNemo

本方案描述了使用失活融合蛋白SID4X-dCas9-KRAB（也称为增强子干扰，Enhancer-i）对增强子进行多重靶向设计与操作的步骤。该方案可用于鉴定调控基因表达的增强子，并有助于解析调控同一靶基因的多个增强子之间的相互关系。...

利用细胞系中基于多重CRISPR的增强子干扰解析增强子功能

本方案概述了构建一种模型系统所需的步骤，该系统可利用增强型lac阻遏蛋白和/或tet激活蛋白系统，在活体动物或细胞中条件性地调控目的内源基因的转录。...

体内应用REMOTE-控制系统调控内源性基因表达

我们报道了在由胚胎干细胞生成的拟胚体中应用定量染色体构象捕获技术并结合高通量测序的方法。该技术可用于鉴定和定量分析在胚胎干细胞分化过程中特定基因的潜在增强子与启动子区域之间的相互作用。...

通过定量染色体构象捕获技术（4C）鉴定类胚体中的增强子-启动子相互作用

本方案描述了一种基于重组腺相关病毒（rAAV）的瞬时增强子报告基因检测方法。该检测可用于诱导由增强子驱动的基因表达 体内...

基于增强子的表达载体的AAV递送 在体 小鼠脑中

在此，我们提出一种用于 体外...

体外筛选靶向表观遗传沉默的工程化转录抑制因子

基因表达受到基因启动子与远端调控元件相互作用的调控。本文介绍了低起始量的捕获染色体构象捕获技术（liCHi-C）如何实现对稀有细胞类型中这些相互作用的鉴定，而这些细胞类型在此前是无法测量的。...

一种用于在稀有细胞群体中研究以启动子为中心的时空基因组结构的整合工作流程

本方案提供了一条简化的计算流程，用于定量新生增强子转录本。通过整合染色质可及性、染色质特征和转录数据，该方法能够在复杂的基因内区域实现对增强子活性的准确检测和链特异性分析，同时适用于不具备深入生物信息学训练背景的研究人员。...

小鼠胚胎干细胞中基因间/基因内增强子RNA定量的计算分析流程

我们介绍一种用于分选单个哺乳动物细胞并量化每个细胞中多达96个目标基因表达水平的方法。该方法包括使用内部qPCR标准品，以实现对转录本绝对拷贝数的估算。...

使用流式细胞分选与含内参标准品的qPCR进行单细胞基因表达谱分析

学习使用基于浏览器的无代码工作流程识别组织特异性基因增强子。该方案利用公共组蛋白标记和Hi-C数据，以TBX5为案例，帮助生物学家无需编程技能即可进行增强子发现。

一种用于筛选基因和组织特异性增强子的基于网络的分析流程

探索更多视频