哪个网站上做自媒体最好,做配资网站多少钱,android开发和网站开发,wordpress 获取指定文章表观基因组学分析是解释非编码基因组功能的重要方法之一。表观基因组学特征作为肿瘤细胞的基本特征#xff0c;对发病机制、临床行为和治疗具有影响。在所有表观基因组标记中#xff0c;组蛋白修饰和 DNA 甲基化已得到最广泛的研究#xff0c;以深入了解表观基因组失调。
染…表观基因组学分析是解释非编码基因组功能的重要方法之一。表观基因组学特征作为肿瘤细胞的基本特征对发病机制、临床行为和治疗具有影响。在所有表观基因组标记中组蛋白修饰和 DNA 甲基化已得到最广泛的研究以深入了解表观基因组失调。
染色质可及性是 DNA 调控元件的一个标志新出现的证据表明它在癌症中发挥着重要作用。使用测序 (ATAC-seq) 检测转座酶可及染色质的方法的出现和优化使得在原发性癌症的全基因组范围内分析染色质可及性成为可能。利用这项技术最近的一项研究表明淋巴细胞特异性开放染色质区域预先决定了 ALL 的糖皮质激素耐药性这表明染色质可及性特征在 B-ALL 耐药和复发中的潜在作用。然而关于原发性儿童 B-ALL 的染色质可及性概况以及复发期间发生的可及性变化的认知仍然缺乏。
2023年10月25日国家儿童医学中心/上海儿童医学中心刘宇和沈树红研究团队在 期刊Nature Communications发表题为 Chromatin accessibility landscape of relapsed pediatric B-lineage acute lymphoblastic leukemia 的研究论文该研究从染色质可及性accessibility的角度对复发B-ALL基因组特征进行了分析。 研究单位国家儿童医学中心/上海儿童医学中心 发表时间2023年10月25日 发表杂志Nature Communications16.6/Q1 研究技术ATAC-seq、RNA-seq和ChIP-seq等爱基百客可提供
01 研究背景
急性淋巴细胞白血病ALL是最常见的儿童癌症。B 系急性淋巴细胞白血病 (B-ALL) 约占儿童 ALL 病例的 80%。对大型队列的基因组分析已识别出 20 多种具有不同遗传改变的 B-ALL 亚型这使得风险分层和精准治疗成为可能。这与其他治疗进展相结合已将患者生存率提高到 90% 以上。然而难治性和复发性B-ALL患者的预后较差5年生存率50% 对于大约一半的儿童 B 系急性淋巴细胞白血病 (B-ALL) 患者来说复发的分子机制仍不清楚。
02 研究思路 03 研究结果
1.儿童B-ALL的染色质可及性景观
从上海儿童医学中心治疗的61名复发B-ALL患者收集的79个儿童B-ALL肿瘤中生成了总共144个染色质可及性图谱。先前的研究也提供了多个基因组测序数据包括来自32名患者的诊断-缓解-复发三组的WGS数据、来自57名B-ALL患者的89个肿瘤的RNA-seq数据以及来自11名 B-ALL 患者的12个肿瘤的H3K27ac ChIP-seq数据图 1a。通过整合WGS的驱动基因组易位、RNA-seq的融合和基因表达特征以及临床测试的核型和FISH 结果来确定每位 B-ALL 患者的分子亚型。
本分析共纳入以下11种B-ALL亚型分别为超二倍体 (n 20)、ETV6::RUNX1 ( n 11)、TCF3::PBX1 ( n 5)、KMT2A 重排 (n 5)、BCR ::ABL1 ( n 3)、BCR::ABL1-like ( n 4)、ZNF384 ( n 3)、PAX5alt ( n 2)、TCF3::HLF ( n 1)、亚二倍体 ( n 1) ), MEF2D (n 1)还有 5 例未分类亚型命名为B-other。
经过质量控制和重复组合后获得了59名患者的75份高质量ATAC-seq 图谱。每个B-ALL样本中鉴定出的ACRs中位数为138,366范围为57,941 至204,563。这些ACRs进一步合并为代表儿童B-ALL队列ACR 758,738 个 ACRsc-ACR。根据表观基因组路线图项目将ACRs注释到八个功能基因组区域。通过在Blueprint epigenomics Consortiumwith ChromHMM中分析从原代B-ALL细胞收集的全基因组组蛋白修饰获得B-ALL基因组的功能划分。
作者观察到在B-ALL肿瘤基因组中ACRs的功能分布具有可比性图 1b。就 ACRs 的数量和开放性而言与活跃基因转录相关的基因组区域显示出更高的染色质可及性图 1b。与增强子区域Enh相关的ACRs中位数占基因组中所有ACRs的31.30%其次是活性转录位点TssA20.38%、转录相关区域Tx6.38%和二价Tss/Enh BivR4.24%图 1b。
另一方面转录抑制相关区域的可及性较低包括PolyComb区域(ReprPC, 8.00%)、异染色质区域(Het, 0.87%)和ZNF基因和重复序列(ZNF/Rpts, 0.03%)(图1b)。此外与活跃转录区域的ACRs相比抑制区域的ACRs更具异质性。令人惊讶的是代表没有明确组蛋白修饰的基因组区域的Quies区域在 B-ALL 中也显示出染色质可及性占所有 ACR 的中位数为 27.95%图1b 。超过一半 (54.94%) 的 c-ACR 与Quies区域重叠。为了进一步表征Quies 区域作者对12个B-ALL肿瘤样本进行了H3K27ac ChIP-seq分析。结果显示64.83%的Quies ACR位于基因区域。在这些基因中中位数为70.59%的基因也显示出 H3K27ac 信号并且基因转录增加。这些数据表明Quies区域的 ACR 参与主动转录的调节。 Fig 1a-b.儿童B-ALL患者染色质可及性模式
2. 与白血病发生相关的B-ALL特异性染色质可及区域
目前认为B-ALL起源于B细胞前体作者将上述鉴定的c-ACRs与先前发表的从胎儿骨髓中分选的pre-pro B和pro B细胞的染色质可及性谱进行了比较这些细胞代表了B细胞祖细胞中可接近的染色质状态。作者发现与pre-pro B和pro B细胞相比B-ALL在整个基因组中染色质可及性的数量上没有显著差异。在pre-pro B细胞(98.57%)和pro B细胞(98.35%)中检测到的大部分ACRs在B-ALL中仍然是开放的(图1c)。这些数据支持B-ALL来源于pre-pro/pro B细胞。另一方面585,248个 (78.39%)ACRs是B-ALL特异性的。进一步分析发现与B-ALL和B细胞祖细胞之间重叠的ACRs相比B-ALL特异性ACRs的异质性明显更高(图1d)与上述B-ALL肿瘤细胞染色质可及性的异质性一致。
接下来作者比较了B-ALL和B细胞祖细胞之间的差异 ACRs。总共 252,028 个ACRs在B-ALL中显示出更高的可及性。这些ACRs位于2332个蛋白质编码基因的启动子区域 (TSS ±1 kb) 内。富集分析表明这些基因与肿瘤相关的生物过程相关包括增殖和分化、信号转导、免疫过程、细胞反应和代谢过程 图1e。其中有61个潜在癌基因包括IL7R、TCL1A、TCF3、RHOA和ELL。如图 1f所示在这些癌基因的启动子区域观察到染色质可及性增加表明这些 ACRs 具有潜在的调节功能。此外在远端增强子区域也观察到染色质可及性增加的 ACRs。图1g中给出了一个例子是远端血液增强子簇BENC据报道它是激活MYC转录的超级增强子在该区域的多个增强子中观察到染色质可及性增加这与B-ALL中MYC活性的增加一致。这些发现表明B-ALL中可及性增加的ACRs参与了疾病的发展。 Fig 1c-g.B-ALL与B细胞祖细胞染色质可及性对比
3. B-ALL 中染色质可及性与亚型特异性转录调控相关
除了B-ALL和B细胞祖细胞pre-pro B和pro B细胞之间的 ACRs 存在差异外作者还观察到分子亚型之间染色质开放程度的显著差异。如图2a所示将B-ALL与复发性c-ACRs进行无监督聚类时将B-ALL按亚型分组。这通过计算 B-ALL 样本之间ACRs的成对相关性得到了支持图2b。所有功能基因组区域均观察到与亚型相关的可及性图 2c包括Enh和 BivR在内的远端调控区显示出最显著的亚型特异性图 2c这与之前报道的远端染色质开放区域的组织特异性一致。接下来作者分析了不同亚型间 ACRs的差异该分析仅包括具有三个以上亚型的病例在亚型之间观察到不同的染色质可及性。作者发现 ETV6::RUNX1和ZNF384 B-ALL样本的 ACRs明显低于其他B-ALL样本。在 625,287个复发性c-ACRs中17,981 个被确定为亚型特异性ACRs每个亚型的ACRs中位数为3083个范围708-5288。分层聚类热图显示这些ACRs显示出很强的亚型特异性可及性图 2d。
通过将转录因子(TF)基序分析与使用RNA-seq数据的基因转录分析相结合作者鉴定了109个与这些亚型特异性ACRs相关的TF。这些转录因子根据它们在每个亚型中的富集程度分为九个簇图 2e。除了针对特定B-ALL亚型富集的TF 之外作者还观察到某些亚型之间TF富集高度相似。这包括TCF3::PBX1和ETV6::RUNX1亚型KMT2A和ZNF384亚型以及BCR::ABL1\BCR::ABL1样亚型和超二倍体亚型的共享TF。其中一些观察结果得到了之前报道的支持。例如KMT2A重排亚型和 ZNF384 B-ALL亚型之间的重叠与这两种亚型都显示出骨髓转录趋势的事实一致。虽然机制仍有待进一步研究但共享的转录调控表明亚型之间细胞分化状态相似。
为了进一步探讨这些TF的调节作用作者对富集亚组的肿瘤样本与其他亚组的肿瘤样本之间的TF及其潜在靶基因进行了表达分析。在富集亚型中发现14个TF的转录显著增加图 2f表明与TF表达直接相关的转录调控。对于靶基因分析作者重点分析了53个在基因启动子区(TSS±1 kb)具有结合基序的TF并将这些基因作为每个转录因子的靶基因进行了分析。进一步组合单个靶基因分析的结果来代表富集亚型中转录因子的调节功能。如图2g所示该分析中包含的53个转录因子中有13个被发现靶基因的表达显著升高支持在富集亚型中增加转录调控活性。作者注意到靶基因分析中的13个TF中有12个没有显示TF本身的表达变化这表明B-ALL亚型之间存在上下文依赖的转录调控。在13个TF 中E2F6被鉴定为在 ETV6::RUNX1 亚型中特异性富集。该基因在细胞周期的控制中起着至关重要的作用并且与多种肿瘤的肿瘤生长或化疗敏感性相关。与此同时与其他亚型相比ETV6::RUNX1 B-ALL亚型中E2F6及其靶基因的转录明显更高(图2f, g)。这些结果进一步揭示了B-ALL亚型特异性转录调控。 Fig2.B-ALL亚型特异性染色质可及性
4. B-ALL 中等位基因特异性开放染色质与白血病相关
对来自32名B-ALL患者的总共44份样本包括13份诊断样本和31份复发样本以及配对的 ATAC-seq 和 WGS 数据进行了等位基因特异性开放染色质 (ASOC) 分析。每个样本平均识别出 3616个ASOC 区域。ASOC区域占全基因组ACRs的中位数14.39%显著低于双等位基因开放染色质 (BiOC) 区域图3a。此外与BiOCs相比在活性转录相关区域TssA、Tx 和 Enh中发现的ASOC区域较少图3b。进一步分析表明与 BiOC ACRs相比ASOC ACRs之间的距离更近并且更有可能被归为单个拓扑相关域TAD。这些数据表明等位基因之间染色质可及性的调节符合三维基因组结构的调节。
接下来作者使用EpiMap 32的白血病相关单核苷酸多态性 (SNP) 研究了染色质可及性。在ASOC区域中总共发现了46个与白血病相关的等位基因之间染色质可及性不平衡的SNP其中至少5个样本中存在的7个SNP图3c, d。在这些最常见的SNP中rs7090445是先前预测的通过破坏RUNX3与c等位基因的结合来减少ARID5B的转录。有趣的是作者观察到在21个具有杂合C/T等位基因的B-ALL样本中有14个样本 (66.67%) 中该位点的 T 等位基因的染色质可及性显着高于C等位基因的染色质可及性这与C等位基因在白血病中的作用一致。另一个复发性SNP是 rs13401811G 等位基因先前被报道为慢性淋巴细胞白血病的风险等位基因。ATAC-seq数据显示9个G/A基因型B-ALL样本中有8个(88.89%)G等位基因的染色质可及性显著高于A等位基因。值得注意的是rs13401811位于BCL2L11上游约262 kb, BCL2L11编码一种参与ALL耐药的促凋亡蛋白20。这些结果表明染色质可及性与这些疾病相关SNP的功能有关 Fig3.B-ALL等位基因特异性开放染色质(ASOC)区
5. B-ALL 治疗后染色质可及性发生变化
作者分析了每个亚型的诊断和复发之间具有差异可及性的ACRs。本分析仅纳入具有超过五个诊断或复发样本的亚型。每个亚型中鉴定出平均有945个差异ACRs范围从268到4072图 4a。值得注意的是在复发样本高复发中只有 1.54%5911例中有91例具有较高可及性的ACRs 中在复发样本低复发中只有0.14%1423例中有2例具有较低可及性的ACRs在两个或多个亚型之间共享这表明不同亚型在治疗期间染色质可及性变化具有显著的异质性图4b。
为了获得这些在复发期间失调的差异ACRs的靶基因作者使用配对的 ATAC-seq和RNA-seq数据对总共52个B-ALL样本进行了ACRs到基因预测并定义了116,307个ACRs-gene相关性。由此总共1259个基因被确定为复发期间ACRs失调的潜在靶标。正如预期的那样观察到显著的异质性在任意两个亚型之间只有13.98%1259个中有176个靶基因是共享的图4cd。富集分析表明这些靶基因与细胞粘附相关的生物过程相关图 4e表明骨髓微环境中白血病细胞和间充质基质细胞之间的相互作用可能失调这在之前已被证明与白血病细胞的化疗耐药有关。
在此基础上作者整合了CTD^236的药物敏感性数据以研究在复发B-ALL患者中观察到的基因失调是否与临床治疗相关。首先作者确定了来自着11个B-ALL细胞系从CCLE项目收集的基因转录数据与细胞对8种常用B-ALL治疗药物(包括阿糖胞苷和甲氨蝶呤)的反应之间的关系。该分析总共产生了14,680对drug-gene对代表与药物反应相关的转录改变。有趣的是由复发相关的ACRs调控的潜在靶基因与药物治疗显著相关包括伊马替尼和依托泊苷图 4f。当作者对个体B-ALL亚型中鉴定的药物基因对进行分析时观察到类似的结果。BCR::ABL1\BCR::ABL1 样亚型差异ACRs的靶基因与伊马替尼和达沙替尼显著相关图4f两者都是用于BCR::ABL1\BCR::ABL1样B-ALL治疗的酪氨酸激酶抑制剂而 ETV6::RUNX1 亚型则观察到与多柔比星显著相关如图 4f。这些结果表明治疗可以重塑染色质的可及性从而影响B-ALL复发期间的基因转录调控。 Fig4. 复发B-ALL患者染色质可及性的改变
6. 染色质可及性特征影响患者预后
作者分析了无复发生存率 (RFS)来研究染色质可及性与B-ALL预后的关系。对42例接受CCCG-ALL-2009 ( n 37) 和 CCCG-ALL-2015 ( n 5)方案治疗的复发B-ALL患者的ATAC-seq数据进行了分析。两种方案之间没有观察到患者预后的显著差异。在 625,287 个复发性 c-ACRs 中总共鉴定出了 70,573 个 (11.29%) RFS相关 ACRs(FDR 0.05)。根据ACRs与基因的联系预测这些RFS相关ACRs的潜在靶标在细胞周期和白细胞分化相关的生物过程中富集图 5a表明这些ACRs对B-ALL母细胞的增殖和分化有调节作用。使用RFS相关ACRs对42例复发 B-ALL 患者进行无监督聚类结果得出两个B-ALL 组A 组和 B 组的复发时间和预后不同图 5b,c。使用不同的聚类方法观察到一致的结果。为了验证这一观察结果作者用来自治疗适用研究以产生有效治疗 (TARGET) 项目的数据进行分析结果表明染色质可及性的异常反映了患者的预后。
与之前的报道一致A组患者富含ETV6::RUNX1和超二倍体亚型预后相对较好而B组患者主要为KMT2A和BCR::ABL1\BCR:ABL1样亚型预后较差图 5b,c。值得注意的是尽管TCF3::PBX1 B-ALL 患者通常被认为是低风险在本分析中所有4例TCF3::PBX1病例与kmt2a和BCR::ABL1\BCR::ABL1样病例在B组进行分组并在诊断后22个月内复发图 5b。有趣的是超二倍体B-ALL病例被分为两个不同的组图 5b。3例复发的超二倍体病例A118R、A174R 和 A233R聚集在B组具有KMT2A和BCR-ABL1\BCR::ABL1样并且显示较差的RFS图5b, c。在两个超二倍体亚组之间总共鉴定出 7566 个差异ACRs图5d。在这3156个ACRs中B 组RFS较差的超二倍体病例中染色质可及性增加。作者从ACR-gene的结果中获得了这些ACRs的603个潜在靶基因。
进一步分析表明这些靶基因模仿干细胞和骨髓祖细胞的表达特征包括巨核细胞-红细胞祖细胞和粒细胞-巨噬细胞祖细胞图5e。从功能上讲目标基因被发现富集在与耐药性相关的迁移/粘附/运动相关类别。这表明与其他超二倍体 B-ALL 患者相比3例超二倍体 B-ALL 患者表现出更高的谱系可塑性具有更多的骨髓样细胞和较差的治疗反应。作者进一步分析了与TARGET队列高二倍体b - all患者不良预后相关的ACRs模式。43例高二倍体样本(7例为诊断-复发配对)根据上述603个靶基因聚类为3组图5f。聚类3由3例患者的5个样本组成其靶基因表达最高代表高危人群。因此簇3患者在EFS和OS方面的预后最差图5g跟作者的观察结果一致。然而该分析的患者数量有限异常ACRs仍有待进一步的大队列研究。 Fig5.染色质可及性特征与B-ALL预后相关
04 总 结
本文的研究展示了儿童B-ALL染色质可及性的概况并描述了B-ALL患者和不同分子亚型中特异富集的ACRs。更重要的是该研究展示了药物治疗下染色质重塑的发生并确定了与早期复发相关的染色质可及性特征。这些结果扩展了研究者对B-ALL背后的基因组畸变的理解并强调了表观基因组特征对于这种恶性肿瘤风险分层的重要性。
