跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
癌细胞。作者手稿;2019年9月10日PMC发布。
以最终编辑形式发布为:
PMCID公司:PMC6158019型
NIHMSID公司:美国国家卫生研究院1504174
PMID:30146332

通过整合基因组、表观基因组和蛋白质组分析确定横纹肌肉瘤的治疗靶点

关联数据

补充资料

总结

针对患者肿瘤中的体细胞突变的个性化癌症治疗越来越多地被纳入实践。由于肿瘤特异性表观遗传扰动导致的基因表达变化导致的其他治疗脆弱性正在逐渐被认识到。这些基因组和表观基因组的变化最终在肿瘤蛋白质组和磷酸蛋白质组中表现出来。我们整合了转录组学、表观基因组学和蛋白质组/磷酸蛋白质组学数据,以阐明横纹肌肉瘤(RMS)的细胞起源和治疗脆弱性。我们发现肺泡RMS在发育过程中比胚胎RMS发生得更远。我们还发现RAS/MEK/ERK/CDK4/6,G2/并通过我们的综合分析揭示了蛋白质反应途径。综合临床前测试表明,靶向G区的WEE1激酶2/M路径是最有效的方法体内用于高风险RMS。

图形摘要

保存图片、插图等的外部文件。对象名为nihms-1504174-f0001.jpg

简介

横纹肌肉瘤(RMS)是一种以骨骼肌为特征的儿童实体瘤。约75%的局限性疾病患者通过传统的多模式治疗治愈,但复发或转移性疾病患者的总生存率分别为17%和30%(Pappo等人,1999年). 历史上,RMS分为胚胎RMS(ERMS)和肺泡RMS(ARMS),前者约占所有RMS病例的60%,后者占25%(牛顿等人,1988年). 大多数组织病理学分类的ARMS包含FOXO1系列在染色体13q14上PAX3系列染色体2q35或PAX7在染色体1p36上(Raney等人,2001年). 然而,ARMS的一个子集为易位阴性,其分子特征与ERMS更为相似(Williamson等人,2010年). ERMS患者预后较好(Ognjanovic等人,2009年)和当前临床试验使用PAX3/7-FOXO1系列易位状态用于风险分配而非组织学亚组。先前的基因组特征表明,易位阴性ERMS具有较高的结构、拷贝数和单核苷酸变异率,而易位阳性ARMS在癌症共识基因中的复发突变相对较少(Chen等人,2013;Shern等人,2014年).

为了提高对复发性RMS的认识并提供更多的临床前模型,我们建立了一组儿童实体瘤原位患者衍生异种移植(O-PDX)(Stewart等人,2017年). 我们所有的O-PDX及其相关数据都通过儿童实体肿瘤网络(CSTN)自由共享,没有义务合作(www.stjude.org/CSTN网站/)由于RMS患者的总生存率在过去15年中没有明显提高,因此这种简化的方法对于阐明有效的RMS治疗方法尤为重要(Smith等人,2014年).

结果

横纹肌肉瘤的表观遗传学特征

此前,我们为17个RMS生成O-PDX,并对其及其相应的患者肿瘤进行了表征(Stewart等人,2017年),根据目前的成人和儿童NCI-MATCH临床试验,其中11个具有潜在的可操作突变(表S1). 在这里,我们扩展了我们的分析,包括这些样本的全基因组亚硫酸氢盐测序(WGBS)。我们还对原始人类成肌细胞、肌管和胎儿骨骼肌(妊娠21周)进行了WGBS和RNA测序(RNA-seq)。大多数患者肿瘤和匹配的O-PDX的基因表达(r=0.926;范围0.76–0.99)和DNA甲基化(r=0.874;范围0.57–0.97)相关,肿瘤类型之间的差异与患者肿瘤中浸润的非肿瘤细胞有关(图1A和1B;表S2). 接下来,我们确定了ERMS和ARMS的DNA甲基化与基因表达之间的相关性差异。在ERMS中相对于ARMS上调的1782个基因中,427个(24%)与DNA甲基化呈负相关,而352个(20%)与279个(16%)呈正相关具有正相关和负相关区域(图1C;表S2). 在ERMS相对于ARMS下调的1576个基因中,134个(9%)与DNA甲基化呈负相关,386个(24%)呈正相关,368个(23%)均呈正相关和负相关(图1C;表S2).

保存图片、插图等的外部文件。对象名为nihms-1504174-f0002.jpg
DNA甲基化与基因表达的差异

(A和B)O-PDX样本和相应RMS肿瘤的RNA-seq数据代表性散点图(log2 FPKM,A)和WGBS分析的β值(B)。

(C) 描述ERMS中与低甲基化(蓝色)或高甲基化(红色)ARMS肿瘤相关的上调和下调基因(绿色)重叠的维恩图。记录各组的基因数。

(D-G)ARMS中相对于ERMS(装箱区)启动子下游高度甲基化并降低基因表达的基因示例(MFAP2型,D),即ARMS中转录起始位点侧翼区域相对于ERMS(装箱区域)发生低甲基化,基因表达增加(NHLH1型,E),即相对于ERMS,ARMS中的高甲基化,ERMS中的基因表达增加(基质金属蛋白酶16,F),该基因在ARMS中相对于ERMS被高甲基化,并增加了ARMS中的基因表达(NOS1号,G)。ERMS(蓝色)和ARMS(红色)绘制了由虚线表示的转录边界和每个CpG的β值。RNA-seq中每个样本的表达如右侧的柱状图(FPKM)所示。FPKM,碎片/千基百万。

另请参见表S1S2系列.

许多DNA甲基化和mRNA表达呈负相关的基因在启动子下游区域被发现,即下游相关区域(Hovestadt等人,2014年) (图1D和1E). 相反,在整个基因座上发现了许多与DNA甲基化和mRNA表达呈正相关的基因(图1F和1G)并含有部分甲基化结构域(李斯特等人,2009年). 差异表达和甲基化基因的通路分析表明,负相关基因主要是胚胎骨骼肌发育的介体,而正相关基因则主要是神经系统发育的中介体(表S2).

为了确定非DNA甲基化的表观遗传变化是否伴随RMS肿瘤发生期间基因表达的变化,我们对八种组蛋白修饰(即H3K4me1、H3K4me2、H3K4me3、H3K17me3、H3K27Ac、H3K 36me3、H 3K9/14Ac和H3K9me3)和CTCF、BRD4、,所有O-PDX、成肌细胞、肌管和胎儿骨骼肌中的RNAPII。我们总共使用了19种抗体,并对756个ChIP-seq文库进行了测序/分析,这些文库的数据用于进行染色质隐马尔可夫模型(chromHMM)(安永会计师事务所(Ernst and Kellis),2012年). 通过结合计算方法和人工审查,我们选择了18个chromHMM状态进行分析(图2A). 状态1至7具有与增强子和活性启动子一致的表观遗传学标记;状态8包含二价(H3K27me3/H3K4Me3)平衡启动子/增强子的表观遗传学标记;状态9至状态12包括基因体的表观遗传标记。状态13与基因组中具有H3K27Ac、H3K9me3和H3K27 me3修饰的大类锌指基因一致(O'Geen等人,2007年). 状态14包含多梳抑制(H3K27me3)基因组区域,在我们的分析中,状态15和状态16为安静/空白状态。状态17由异染色质组成,状态18代表受抑制的染色质(H3K27me3和H3K9me3)。ERMS和ARMS的这18个chromHMM状态在基因组中的分布非常相似(图2B). chromHMM状态的差异与这些肿瘤之间基因表达的差异密切相关(图2C)表明ChIP-seq和chromHMM分析的蛋白质和组蛋白修饰提供了比单独DNA甲基化更完整的表观基因组景观视图。

保存图片、插图等的外部文件。对象名为nihms-1504174-f0003.jpg
与启动子/增强子活性相关的表观遗传谱差异

(A) 此分析中使用的18个chromHMM状态的热图。

(B) 热图描绘了ERMS和ARMS中基因组注释区域的18个chromHMM状态的比例。单个标记或区域的比例与蓝色的强度成正比。

(C) ChIP-seq峰值和chromHMMMYOG公司。右侧显示刻度,基因边界用虚线标记。

(D-F)代表性chromHMM气体安全标准2在ERMS(D)中选择性表达,NOS1号在ARMS(E)中选择性表达AXL公司在成肌细胞(F)中选择性表达。编号1包含基因上游的肿瘤类型特异性超增强子(装箱的H3K27Ac区域)。基因边界用虚线标记,每个状态的颜色如(A)所示。比例最高的两个州是全高酒吧,其余州是半高酒吧。每个条形图的强度与该基因在所有阶段的每个状态的百分比成正比。对于半高杆,强度从最大强度的50%开始缩放。

(G) ERMS中相对于ARMS的上调或下调(绿色)基因的文氏图。与在chromHMM状态(蓝色)或肿瘤类型特异性超增强子(红色)上有差异的基因重叠显示。

另请参见表S3.

我们将表观基因组分析重点放在三大类结合DNA甲基化的研究上:(1)ERMS和ARMS之间差异表达的基因,其chromHMM状态发生相应变化(图2D;表S3), (2)与转录起始点500kb范围内的超增强子相关的肿瘤类型之间差异表达的基因(图2E;表S3)、和()H3K27me3差异基因(chromHMM状态8、13、14、18)与基因表达负相关(图2F;表S3). 相对于ARMS,chromHMM状态的差异导致ERMS中1782个上调基因中的749个(42%)和1576个下调基因中的522个(33%)(图2G). 只有36(2%)个上调基因和90(6%)个下调基因含有H3K27me3修饰(表S3). 重要的是,在ERMS中上调的98个基因包含ERMS特异性超增强子,在ARMS中上调了的174个基因包含ARMS特异性超级增强子(图2G;表S3)其中许多与肌肉发生有关。

为了将ChIP-seq数据与转录组和DNA甲基化数据整合,我们将ERMS和ARMS特异性H3K27me3、chromHMM和超增强信息添加到表S2在ERMS中相对ARMS上调并低甲基化的427个基因中,有4个基因H3K27me3相对降低,47个基因具有ERMS特异性增强子,174个基因具有与表达相关的chromHMM差异(表S2). 例如,KCNK5公司在ERMS中相对于ARMS上调,与差异基因表达一致的chromHMM变化与ERMS特异性超增强子相关,并与H3K27me3降低低甲基化(表S2).KCNK5公司编码与肌肉分化和肌细胞融合有关的钾通道(Afzali等人,2016年)ERMS肿瘤的表观遗传上调基因参与细胞外基质组织和胚胎形态发生,特别是四肢和骨骼系统的胚胎形态发生(表S3). 相反,ARMS肿瘤表观基因上调与肌肉分化有关的基因,表明其在肌肉发育的后期被阻止(表S3).

我们使用核心调节电路(CRC)映射(Saint-Andre等人,2016年)在我们的ChIP-seq和RNA-seq数据中确定骨骼肌发育中的核心转录因子,并确定在肌发生中重要的MYOD1、PAX7、FOXO1、MYF6、MEF2A和MYOG(表S3). 其他研究较少的核心转录因子包括FOSL2,这是最常见的CRC。结合SMAD3和TGIF1的CRC,提示TGF-β信号可能在RMS中起重要作用(Carthy等人,2015年) (表S3). 同样,刺猬(HH)、notch和维甲酸信号也与RMS有关(表S3).

横纹肌肉瘤蛋白质组和磷蛋白组的定量分析

我们使用了具有广泛肽分离能力和高质量分辨率的10-plex等压标记MS管道(Bai等人,2017年)在重复实验中量化12种OPDX和人类成肌细胞和肌管的蛋白质组和磷酸蛋白质组(表S4). 对三批10份样品进行溶解、消化,并用10个不同的串联质量标签(TMT)进行标记(图3A),每批中包含一份技术复制品作为对照(表S4). 用5%的样品池对均匀混合的样品进行分馏,以进行全利润分析(15.9×106MS/MS光谱)和其余95%的磷酸蛋白质组(8.0×106MS/MS光谱)(图3A). 我们分析了16523个人类蛋白质(12926个基因)和7177个蛋白质的62965个磷酸化事件。我们定量了来自10333个基因和12653个磷酸化位点的分批蛋白质产物,并通过免疫印迹验证了几种蛋白质和磷酸化蛋白(图3B和3C). 整个蛋白质组和磷酸蛋白质组的组内再现性均为0.98(图3D和3E),主成分分析显示,每种样品类型都有明显的分离(图3F和3G). 加权基因共表达网络分析(Langfelder和Horvath,2008年)鉴定了在ARMS、ERMS和发育中的肌肉中差异表达的六组蛋白质和磷蛋白(表S4).

保存图片、插图等的外部文件。对象名为nihms-1504174-f0004.jpg
横纹肌肉瘤蛋白质组和磷蛋白组的分析

(A) RMS肿瘤、肌管和成肌细胞的蛋白质组和磷酸蛋白质组分析的工作流程。液相色谱法;质谱法。

(B) 代表性免疫印迹验证NOS1、RPS6的磷酸化Ser235/236和β-肌动蛋白的蛋白和磷酸蛋白表达(负荷控制)。

(C) (B)和HMGA2中蛋白质的标准化相对折叠表达条形图。mb,成肌细胞;mt,肌管。

(D和E)相同批次和不同批次中SJRHB000026_X1重复样品的蛋白质(D)或磷酸蛋白(E)定量散点图。

(F和G)整个蛋白质组(F)和磷酸蛋白质组(G)的成肌细胞(浅绿色)、肌管(深绿色)、ERMS(红色)和ARMS(紫色)的主成分分析。RNA表达(FPKM)的(H和I)条形图MYOG公司(H) 和第五个多年电价(一) 发生肌肉、ARMS和ERMS肿瘤。红线表示人类成肌细胞的表达水平。ERMS高组是MYOG水平最高的样本,ERMS低组是其余样本。FQ21,妊娠21周时胎儿股四头肌。

(J) 原木散点图10来自MYOG和MYF5蛋白质组分析的蛋白质表达强度。线性回归(虚线)中排除了两个异常值样本(黑箱)。

(K) MYOG和MYF5样本的代表性免疫染色,以及每个O-PDX肿瘤的MYOG和MYF5免疫阳性细胞百分比散点图,具有最佳拟合线性回归(虚线)。比例尺=20μm。

另请参见图S1-S5表S4第5章.

我们对成肌细胞、肌管和RMS之间差异表达的蛋白质和磷酸蛋白进行了通路分析(Cossu和Borello,1999年),小时(Borycki等人,1999年),PI3K(Jiang等人,1998年),第38页/MAPK(Bennett和Tonks,1997年)、BMP(Reshef等人,1998年),和腺苷酸环化酶(Chen等人,2005年) (图S1-S5).

MYOG是一种与肌肉分化相关的转录因子,在ERMS和ARMS中表达(Hostein等人,2004年) (图S1B-C). ARMS表达的MYOG mRNA和蛋白水平高于ERMS,ERMS呈分级表达(图(图3H3小时S1C系列). 因为第五个多年电价在肌源性测定期间,由增殖的成肌细胞瞬时表达,我们测定了单个肿瘤中的表达变异性及其转录靶点。所有ARMS均显示较低水平第五个多年电价与肌肉相比,ERMS的一个子集表现出较高水平的MYF5 mRNA和蛋白质(图(图3I第三章S1C系列). 在大多数RMS中,MYF5和MYOG蛋白表达呈负相关(图3J)和MYOG转录靶基因以相同模式激活(图S1D第1版E;表S5). 我们对17例RMS肿瘤及其匹配的O-PDX进行了MYF5和MYOG免疫组化染色,发现染色、基因表达和蛋白质组数据之间存在一致性(图3K,表S5).

使用多个平台的集成分析

为了跨多个平台集成我们的数据,我们使用了两种不同的方法。首先,我们鉴定了与人类成肌细胞相关的肿瘤中差异表达的蛋白质(7546/10333)和mRNA(7317/12362),并对整合后跨两个平台保守的基因/蛋白质进行共表达聚类,包括四个集成蛋白簇(IPC)(图4A和4B;表S6). IPC1包含WEE1、CHEK2和RB1等基因/蛋白质,这些基因/蛋白质在肿瘤中相对于成肌细胞上调。IPC2包含在肿瘤中上调的基因/蛋白质,如MYC、GLI1/2和ATM,但ERMS中上调的水平高于ARMS。IPC3包含在ARMS中上调特异性的基因/蛋白,如NOS1、MYOG和CDK4。IPC4包含诸如AXL、MYH9、,与成肌细胞相关的肿瘤中下调的MURC。每个IPC用于通路分析,通路分析上覆盖有蛋白质相互作用(PPI)网络数据(图4C-E).

保存图片、插图等的外部文件。对象名为nihms-1504174-f0005.jpg
转录组、表观遗传和蛋白质组/磷酸蛋白质组数据的整合

(A) 基于差异表达聚类和PPI网络模块的蛋白质组和转录组数据整合方案。

(B) 综合分析中每个IPC中单个蛋白质归一化相对丰度的方框图。矩形表示第一到第三个四分位的四分位范围(IQR),方框内的水平线是中间值。在矩形上方和下方延伸的垂直线表示最大值和最小值,不包括IQR的3倍以上的异常值。将数据标准化为成肌细胞,并显示每个IPC中的基因/蛋白质数量。每个IPC中的代表性基因/蛋白质以灰色显示。

(C) 日志的热图10富集FDR用于每个IPC中基因/蛋白质的通路分析。对于IPC3,该集群中没有路径达到统计显著性(p值<0.05)。(D) IPC1的代表性PPI模块,显示E2F模块(蓝色圆圈)和G2/M个检查点模块(虚线圆圈)。

(E) G的扩展视图2/M检查点模块显示代表性基因/蛋白质。

另请参见表S6.

在我们的第二种方法中,我们将转录组、表观基因组、蛋白质组和磷酸蛋白质组数据与顺序统计相结合。对于表观遗传学分析,计算每个样本中每个基因的活跃状态(1、2、3、4、5、6和10)和抑制/不活跃状态(14、17和18)之间的分数差异,并根据对数进行排序2每个肿瘤相对于成肌细胞的折叠变化。对于蛋白质组和磷酸蛋白质组数据,计算每个蛋白质的总磷酸化水平或相对磷酸化水平,并根据对数对单个蛋白质进行排序2每个肿瘤与成肌细胞的比率(表S6).

MYOG是ARMS中排名最高的基因/蛋白质之一,MYF5在ERMS中排名第一(图3小时–3公里;表S6). ARMS相对于ERMS的上调途径为肌发生(表S6). 我们观察到ARMS和ERMS之间NOS1 mRNA/蛋白的差异表达(图(图2E第二版和3B)。3B公司). NOS1在ARMS中排名最高(表S6). 我们还发现了上调的G2/ARMS和ERMS中与成肌细胞相关的M和E2F靶点(表S6).

横纹肌肉瘤中RAS-CDK4/6通路的靶向性

我们分析了患者生殖系组织、RMS和匹配O-PDX的WGS和WES数据。我们在RMS中发现18个具有致病性癌症突变的基因,其中10个包含在NCI-CG儿科MATCH试验所用的致癌物综合分析小组中(表S1第7部分). NOTCH1、NOTCH3和p53的有限药物敏感性使这些靶点无法用于儿童或成人NCI-MATCH试验。然而,剩下的七个基因(NRAS、HRAS、NF1、CDKN2A、PIK3CA、FGFR1、和FGFR4公司)进入通过Rb/E2F途径调节增殖的信号转导级联(图5A). 虽然在我们的分析中,该途径中大多数单独突变的基因排名不高,但相对于成肌细胞,RMS中的RAS、mTORC1和PI3K-AKTmTOR信号通路被解除调控(表S7). 此外,我们发现这些途径下游的E2F靶点受到强烈的放松管制(图4C、4D、和第6a条;表S6).

保存图片、插图等的外部文件。对象名称为nihms-1504174-f0006.jpg
横纹肌肉瘤靶向RAS和CDK4/6

(A) RAS和CDK4/6通路的简化路径图。编码以紫色和粗体显示的蛋白质的基因的复发突变是以小分子抑制剂为靶点的。(B) 代表性ChIP-seq轨道和chromHMMCDK1型。标尺显示在右侧。为了简单起见,只显示了一种拼接变体。

(C) 成肌细胞、肌管和RMS肿瘤的标准化相对CDK1蛋白表达条形图。

(D) ERMS和ARMS肿瘤切片的代表性显微照片,苏木精和伊红(H&E)染色或细胞周期蛋白D2(ERMS)、细胞周期蛋白D3(ARMS)和磷酸化RB1(pRB1)的免疫组织化学染色。比例尺=完整图像中的50μm,放大插图中的10μm。

(E) RD细胞中10种浓度帕博西利和4种浓度曲美替尼的代表性联合药物敏感性研究。存活率是相对于未经处理的细胞(100%存活率)和使用阳性对照细胞毒性化合物的完全消融(0%存活率)绘制的。

(F) O-PDX肿瘤(蓝色)和血浆(红色)中帕尔博西利的PK图,每个治疗组3只小鼠。绘制每个时间点的平均偏差和标准偏差。

(G和H)SJRHB013759_X1 O-PDX肿瘤负荷的代表图,通过随时间变化的生物发光(G)和安慰剂对照组和帕博西利+曲美替尼治疗组PD小鼠的图像进行测量(H)。

另请参见图S6表S7.

表观遗传学分析证实E2F转录靶点如CDK1在RMS中被激活(图(图5B第5页第6页). 蛋白质组学/磷酸蛋白质组学和免疫组织化学染色数据表明,细胞周期蛋白D/CDK4/6-RB/E2F通路在O-PDX中是活跃的(图(图5C,5摄氏度,、5D、,第五天,S6B系列第6d页;表S7).

将CDK4/6抑制剂与MEK抑制剂结合是NCIMATCH试验中使用的既定策略(Vora等人,2014年)正在考虑参加NCI-CG儿科MATCH试验(Allen等人,2017年). 我们测试了MEK抑制剂(曲美替尼)与CDK4/6抑制剂(帕尔博西利或阿贝马西利)的联合使用,发现这些联合使用会协同降低RAS系统-突变型RMS(RD)和肺癌(A549)细胞株,但不影响皇家银行-缺陷型骨肉瘤细胞系(SAOS-2)(图5E;表S7数据未显示)。然而,RAS途径突变与对MEK抑制剂的敏感性或CDKN2A型CDK4/6抑制剂作为单一药物的突变和敏感性(https://braid.stjude.org/mastour网站/).

为了建立临床相关的小鼠等效剂量(MED),我们对每种药物在荷瘤小鼠中进行了药代动力学(PK)研究。我们还测量了每种药物对肿瘤的渗透性(图5F)联系体内暴露于在体外使用先前描述的数学模型的灵敏度(Stewart等人,2014年). 我们预测,将达到足够水平的帕尔博西卜和曲美替尼或阿贝马西卜和曲美替尼,以产生抗肿瘤作用体内(表S7). 如前所述,进行临床前1期耐受性研究(Stewart等人,2014年).

接下来,我们建立了两个ERMS(SJRHB000026_X1和SJRHB 012_Y)和一个ARMS(SJ RHB 013759_X1)的荧光素酶标记O-PDX(表S1). SJRHB000026_X1致癌人力资源管理系统PIK3CA公司突变,而SJRHB012_Y在该途径中没有突变(Stewart等人,2017年). 如前所述,我们在一项随机临床前2期试验中测试了帕尔博西利+曲美替尼和阿贝马西利+曲美替尼,并将其与复发性RMS标准治疗方案-利诺替康(IRN)+长春新碱(VCR)进行了比较(Stewart等人,2014年). 用IRN+VCR治疗的小鼠有稳定的疾病(SD),而用帕尔博西利+曲美替尼或阿贝马西利+曲美替尼治疗的小鼠表现为进行性疾病(PD;图5G和5H;表S7).

瞄准G2/横纹肌肉瘤中的M和未折叠蛋白反应(UPR)通路

因为以RAS-CDK4/6途径为目标体内是无效的,并且RMS中几乎不存在药物突变,我们使用我们的集成数据集来确定治疗相关的肿瘤脆弱性。我们重点研究了综合分析中阐明的途径,这些途径也通过O-PDX和细胞系的原代培养物表现出RMS选择性药物敏感性(https://braid.stjude.org/masttour网站/) (Stewart等人,2017年). 在RMS O-PDX中总共分析了1767种药物/肿瘤组合。最有效的药物是HDAC和蛋白酶体抑制剂(表S8). 在我们的综合分析中,这些途径没有受到明显干扰,以前使用这些药物(帕诺司他和硼替佐米)进行的临床前2期研究也没有显示疗效体内(Stewart等人,2017年). 针对UPR和G的HSP90抑制剂ganetespib(GSP)和WEE1抑制剂AZD17752/我们的综合分析分别发现M个检查点路径也显示出显著性在体外活动(表S8).

许多HSP90客户是信号转导途径的关键调节者(Karagoz和Rudiger,2015年) (图6A). 由于热休克转录因子HSF1的过度磷酸化是细胞应激的标志(Akerfelt等人,2010年),我们发现HSF1的蛋白质水平和磷酸化增加(图6B;表S4)以及HSF1靶基因表达,在RMS肿瘤中(图6C–E段;表S8). RMS细胞系和O-PDX细胞的原代培养物对GSP也比其他儿科实体瘤更敏感(图6F). GSP快速诱导RD细胞表达热休克蛋白70和其他HSF1靶基因(图6G). 我们筛选了156种药物和几种浓度的GSP(100 nM、33 nM、10 nM和3 nM)。我们使用了对加性相互作用剂量模型的双变量响应(Twarog等人,2016年)分析每种药物组合的效果,发现GSP增强了化疗诱导的RD和O-PDX细胞死亡(表S8).

保存图片、插图等的外部文件。对象名为nihms-1504174-f0007.jpg
HSP90在横纹肌肉瘤中的靶向作用

(A) RMS中解除调控并被HSP90调节的六条肌源性通路的图示。(B) 成肌细胞、肌管和RMS肿瘤中HSF1和S368磷酸化HSF1的条形图。数据标准化为成肌细胞(红线)。

通过RNA-seq测量,RMS肿瘤中396个HSF1靶基因(C)定量表达的(C和D)条形图和HSF1靶蛋白基因表达的代表性示例(D)相对于成肌细胞的表达进行了归一化。条形图是多个肿瘤的平均值和标准偏差。

(E) HSF1靶基因的ChromHMM(热休克蛋白A1)通过RNA-seq(FPKM)进行相应的基因表达。(F) 原木热图10测试药物子集的50%有效浓度。粗体字表示O-PDX的原代培养,而非粗体字则表示已建立的细胞系。(G) 单独使用GSP或IRN+VCR治疗后RD细胞HSP70表达的免疫印迹。对照组表示未经处理的细胞。在计算处理和未处理样品之间的折叠差异之前,对HSP70强度进行量化并归一化为β-肌动蛋白。

另请参见表S8.

AZD1775的靶点WEE1在RMS中的表达水平高于其他儿童癌症,但Wilms肿瘤和视网膜母细胞瘤除外(https://pecan.stjude.org网站) (周等,2016) (图7A). WEE1调节G2/细胞周期S期M检查点和DNA复制(Do等人,2013年) (图7B). RMS相对于成肌细胞的表观遗传学、基因表达和蛋白质组/磷酸蛋白质组数据揭示了该途径的上调(图7C和7D;表S6). 用AZD1775处理RMS细胞系导致G2/M期阻滞、有丝分裂突变和核分裂(图7E–7G). AZD1775与IRN或VCR联合使用可导致高DNA损伤率(图7G)提示WEE1和HSP90可能是RMS的药物靶点。

保存图片、插图等的外部文件。对象名为nihms-1504174-f0008.jpg
横纹肌肉瘤靶向WEE1

(A) 的方框图WEE1号机组儿童肿瘤的基因表达https://pecan.stjude.org/home网站矩形表示第一到第三个四分位的四分位范围(IQR),方框内的垂直线是中间值。从矩形延伸的水平线表示最大值和最小值,不包括IQR 3倍以上的异常值(圆)。脉络丛癌;EPD,室管膜瘤;HGG,高级别胶质瘤;LGG,低度胶质瘤;MB,髓母细胞瘤;急性髓细胞白血病;B细胞急性淋巴细胞白血病;MLL,混合系白血病;T细胞急性淋巴细胞白血病;肾上腺皮质肿瘤;MEL,黑色素瘤;NBL,神经母细胞瘤;骨肉瘤;视网膜母细胞瘤;WT、Wilms肿瘤。

(B) WEE1在G处规定的路径的简化图2/M检查点。箭头表示激活,条形表示压制。

(C) ChromHMM用于WEE1号机组样本中有相应的RNA-seq表达。(D) 成肌细胞、肌管和RMS肿瘤中归一化为成肌细胞的WEE1蛋白水平条形图(红线)。

(E-G)用AZD1775(E)处理前(黑色)和处理后(红色)RD细胞的DNA含量和细胞周期分布,用DAPI染色的这些细胞细胞核的代表性显微照片,突出显示细胞核碎片/有丝分裂突变(F),以及处理后具有与有丝分裂灾难一致核特征的细胞比例(G).F中的比例尺=2μm。

G中的横线是重复计分的平均值和标准偏差。

以G为目标2/M和UPR路径在体内

为了建立GSP的临床相关MED,我们对荷瘤小鼠进行了PK研究(图8A); AZD1775的PK先前已完成(Stewart等人,2017年). 肿瘤中GSP和AZD1775的持续水平高于培养中增强细胞毒性作用所需的水平(图8A和8B). 我们对临床相关的GSP和AZD1775组合进行了临床前1期研究(图8C)如前所述(Stewart等人,2017年;Stewart等人,2014年). 两种常用IRN方案治疗儿童RMS和其他实体瘤(Furman等人,1999年)包括在内。所有组对GSP(10 mg/kg)和AZD1775(60 mg/kg,每天两次)的MED耐受性良好(数据未显示)。

保存图片、插图等的外部文件。对象名为nihms-1504174-f0009.jpg
AZD1775和Ganetespib的临床前测试

(A) 使用SJRHB000026_X1 O-PDX原代培养物的GSP和AZD1775的剂量-反应曲线。

(B) O-PDX肿瘤(蓝色)和血浆(红色)中GSP的PK图,每个治疗组3只小鼠。绘制每个时间点的平均偏差和标准偏差。(C) 为每个联合治疗方案选择药物方案。显示了IRN的每日×5×2低剂量延长时间表,但我们也测试了替代的每日×6时间表。(D) O-PDX肿瘤(SJRHB012_Y)小鼠经IRN+VCR或GSP+IRN+VC治疗后的异种基因图像,这些小鼠患有PD或CR。

(E) 在一项临床前第2阶段研究中,通过生物发光随时间变化测量的肿瘤负担的代表性图,显示了三种O-PDX肿瘤中的每一种。每一行都是不同的鼠标。

(F) 安慰剂对照组、AZD1775组和GSP+IRN+VCR组中PD小鼠及其相应肿瘤的代表性图像。IRN+VCR和AZD1775+IRN+VC治疗组分别显示SD和CR小鼠。(G) 在SJRHB012_Y O-PDX模型的临床前3期研究中,通过生物发光随时间变化测量的肿瘤负荷线图。

(H) 所示四个治疗组的SJRHB012_Y O-PDX生存曲线。(一) 所有14个治疗组的四种O-PDX模型的CR和PR百分比直方图。对AZD1775和GSP的两种临床相关剂量(高(h)和低(l))以及小儿癌症患者使用的几种不同的IRN时间表和给药方案进行了测试。*与IRN+VCR相比,表示p值<0.05。

另请参见表S9.

接下来,我们使用荧光素酶标记的O-PDX进行了临床前2期试验,如前所述(Stewart等人,2017年;Stewart等人,2014年). 小鼠被随机分为八个治疗组:AZD1775+VCR+IRN、AZD11775+IRN,GSP+VCR+IRN、GSP+IRN,VCR+IRN、GSP、IRN或安慰剂,并接受四至六个疗程的治疗(每个疗程21天)(表S9). 每周通过生物发光测量肿瘤负荷,一些用GSP+IRN+VCR或AZD1775+IRN+VCR治疗的小鼠表现出部分反应(PR)或完全反应(CR)(Stewart等人,2017年)(图8D和8E;表S9).

为了进一步探讨GSP和AZD1775联合治疗的疗效,我们进行了一项盲法、随机、安慰剂对照的临床前3期试验。为了概括RMS患者的异质性,我们使用了两个来自复发患者的ERMS O-PDX(SJRHB000026_X1,SJRHB 012_Y)和两个ARMS O-PDXs,一个来自初始诊断(SJRHB013757_X2),另一个来自复发(表S1). 我们招募了499只小鼠,并将其随机分为14个治疗组(表S9). 我们包括三个不同的标准对照组(VCR+IRN),它们具有不同的IRN剂量和时间表(表S9). 四个治疗组加入了GSP,另外四个组加入了AZD1775和标准药物(表S9). 这些数据的子集(四个治疗组,140只小鼠)已在其他地方发表(Stewart等人,2017年)但由于它们是在499只小鼠的队列中随机分组的,因此被纳入此处进行直接比较。单剂GSP和AZD1775组小鼠的总体存活率与安慰剂组无差异(表S9). 虽然一些单独O-PDX的CR和PR发生在GSP+VCR+IRN中,但对于四个O-PDX中的任何一个,GSP+VC+IRN组的总体生存率与VCR+IR组没有差异(表S9).

与VCR+IRN、AZD1775单用或安慰剂相比,所有四种O-PDX在AZD11775+VCR+IR治疗时的疗效都有所改善(图8F–8I;表S9). 当我们在与NCI-CG儿科研究最相似的剂量组中测试AZD1775+IRN时(NCT02095132型)我们发现,与IRN+VCR相比,只有一种O-PDX(SJRHB012_Y)的反应明显更好(p=0.02);然而,在这些给药组中添加VCR显著提高了两种ERMS O-PDX受试者的生存率(图8F–8I;表S9). 四个O-PDX的汇编回复显示,AZD1775+VCR+IRN产生70%的CR+PR(26/37),而AZD11775+IRN则产生39%的CR+RR(14/36)(表S9). 个体O-PDX对VCR+IRN、GSP+VCR+IR和AZD1775+VCR+IRN的反应存在显著差异,这表明使用多种OPDX和相关标准方案进行比较的重要性(图8I;表S9).

讨论

尽管一些试验已经证明了精确药物治疗癌症儿童的可行性(Harris等人,2016年;Parsons等人,2016)这些试验对儿科癌症提出了一些独特的挑战。首先,由于患者人数通常较少,分子靶向治疗是根据基因损伤而非癌症类型来制定的。如果突变的细胞环境影响治疗的结果,那么反应将很难解释(Schneider等人,2017). 其次,目前儿童肿瘤中很少有可操作的突变可通过药物治疗,这使得难以招募足够的患者来达到统计显著性(Allen等人,2017年;Harris等人,2016年;Parsons等人,2016年). 这些障碍可能会妨碍人们仅从靶向临床试验中辨别精确药物对癌症儿童的巨大益处。

先前的基因组研究表明,大约50%至75%的中危和高危RMS患者的RAS途径发生突变(Chen等人,2013). 为了涵盖最大的复发性RMS患者群体,我们使用RAS通路的广义定义来包括NRAS、HRAS、KRAS、NF1、FGFR1/4PIKC3A项目.在一些成年人中RAS系统-以RAS和PI3K通路中的激酶为靶点的突变型癌症已显示出前景(Alagesan等人,2015年;Jokinen和Koivunen,2015年). 我们的综合分析表明,儿科RMS也通过RAS通路表现出强大的信号传导。尽管联合使用CDK4/6和MEK抑制剂在培养中协同诱导RMS细胞死亡,但它们未能达到疗效体内我们的PK研究表明,这并不是因为药物未能穿透肿瘤。其他信号转导途径的代偿性改变可能使肿瘤继续生长(Renshaw等人,2013年).

为了超越基因组分析并确定RMS中可能解除调控的基本途径,我们进行了综合转录组学、表观遗传学和蛋白质组学/磷酸蛋白质组学分析。我们确定了UPR和G2/M-有丝分裂纺锤体检查点通路是RMS综合分析中变化最大的。其他人最近证实,蛋白质稳态途径是RMS的基本脆弱性(Sabnis等人,2016年). 通过用GSP靶向HSP90,我们发现在RMS肿瘤中UPR通路很脆弱,部分原因是HSP90维持6条肌源性信号转导通路的负担增加。HSP90还调节G2/M检查点和有丝分裂纺锤体组装途径通过WEE1。由于HSP90水平几乎不足以控制RMS肿瘤中致癌转化所施加的细胞应激,因此化疗药物,尤其是针对HSP90的化疗药物,可能会增加细胞应激,从而达到合成致死性。事实上,我们观察到RMS细胞系和O-PDX对GSP的敏感性高于其他儿童实体瘤,WEE1抑制也是如此。此外,将GSP与亚治疗剂量的化疗药物相结合,这些化疗药物对自身协同消融培养的RMS细胞没有影响。我们发现GSP+VCR+IRN在临床前1期研究中耐受性良好,在2期研究中部分有效,但在盲法、随机和安慰剂对照的3期研究中无效。这突出了我们的1/2/3期方法的价值,该方法在对RMS儿童进行临床试验之前验证了治疗组合。尽管以UPR通路为靶点并不有效,但以该通路的其他成分为靶点,如HSP70,可能会在未来的研究中提高疗效(Sabnis等人,2016年).

我们研究的一个局限性是偏向于侵袭性难治性RMS肿瘤。我们的综合分析并不全面,因为它不包括低风险RMS,并且大多数O-PDX来源于复发肿瘤。这种局限性可能延伸到假定的可操作突变和药物敏感性之间的相关性。在我们的高通量筛选中测试的任何药物,包括AZD1775和MEK抑制剂,其可操作突变与药物敏感性之间没有相关性,这些药物被认为与RAS途径突变的肿瘤有关(https://braid.stjude.org/mastour网站/).

我们临床前研究的另一个局限性是缺乏药效学。虽然我们观察到有足够的药物渗透以引发抗肿瘤作用,但我们并没有直接证明这一点体内通过有效的药效学分析。因此,HSP90和UPR途径的其他成分可能仍然是RMS与未来几代药物的合理靶点。然而,AZD1775+VCR+IRN在我们的临床前研究中是有效的,我们的发现可能会为高危RMS患者的未来临床试验提供信息。根据我们广泛的第3阶段临床前研究,强调需要将VCR与AZD11775和IRN结合使用,以提高RMS患者反应的可能性,这一点很重要。

STAR方法

试剂和资源共享联系人

有关资源和试剂的更多信息和请求,请直接联系首席联系人Michael A.Dyer,并由其完成(gro.edujts@reyd.leahcim公司).

实验模型和受试者详细信息

动物和原位患者衍生异种移植物

本研究严格按照国家卫生研究所《实验动物护理和使用指南》中的建议进行。该方案得到圣犹德儿童研究医院动物护理和使用机构委员会的批准。所有的努力都是为了减少痛苦。所有的老鼠都是按照批准的IACUC协议饲养的。动物被安置在12-12光圈(早上6点照明,下午6点关闭)内,并提供食物和水随意.从查尔斯·里弗实验室(Charles River Laboratories)购买裸鼠(品系代码553)。NOD/SCID/伽马雌性小鼠购自Jackson Laboratories(菌株代码005557)。如前所述,生成并扩展了人O-PDX(Stewart等人,2017年). 此处描述的所有O-PDX样品可通过CSTN免费获得(http://www.stjude.org/CSTN网站/)没有合作的义务。

方法详细信息

基因组DNA提取

按照制造商说明(Qiagen 69504),使用DNeasy Blood&tissue试剂盒从15-20 mg组织中提取DNA。

染色质免疫沉淀和测序(ChIP-seq)

细胞固定

双份O-PDX肿瘤用于ChIP-seq。每个异种移植瘤(2-3克)直接从免疫功能低下小鼠后腿的注射部位采集。安乐死后立即采集组织。只有单侧肿瘤被用于生物复制。肿瘤切除后30分钟内,将组织切碎(1–3 mm)室温下,使用无菌手术刀片或剃须刀,在不含钙镁(PBS)的磷酸盐缓冲盐水(Lonza 17–516F)中。用新鲜的1%甲醛(Thermo scientific 28906)固定组织,轻轻混合不超过10分钟,包括离心时间。固定8分钟后,在4°C下以1500 RPM的转速离心组织2分钟。去除甲醛,用新制的1X甘氨酸(0.125 M)清洗组织5分钟,然后再次离心。用冷PBS洗涤组织一次,并用在冷PBS中以1:200稀释的蛋白酶抑制剂(Diagenode C12010011)重悬。将组织置于浸泡的组织研磨机中,在冰上浸泡20-40次,将组织分解成单细胞悬浮液。将细胞通过40μm筛网过滤器滤入冰冷的50 ml锥形管中。细胞在10×10中快速冷冻6等分样品,并在−80°C下储存。

细胞核提取和剪切

按照制造商说明,使用TruChIP染色质剪切试剂盒(Covaris 520154)进行细胞核提取和染色质剪切。简言之,10×10的原子核6使用Covaris E系列聚焦超声仪在0.1%SDS剪切缓冲液中剪切细胞12分钟(目标碱基对=350 bp,占空比=5%,强度=4200个周期/脉冲)。使用20μl经0.4 mg/ml蛋白酶K处理的染色质,在65°C下放置20分钟,分析染色质剪切效率。然后将剪切后的染色质在1%琼脂糖凝胶和Agilent 2200 TapeStation上运行,以检查最佳片段大小。

抗体优化、ChIP程序和实时PCR

CTCF、H3K4me1、H3K4me2、H3K 4me3、H3K17AC、H3K-9/14AC、H3 K36me3、H 3K27me3和H3K9me3的抗体优化和标准操作程序可通过儿童实体肿瘤网络获得(www.stjude.org/CSTN网站). 用Active Motif对所有O-PDX肿瘤和对照细胞进行RNA poll和BRD4的ChIP和文库准备。对于其余抗体,按照制造商说明,使用iDeal ChIP-seq试剂盒(Diagenode C01010051)执行ChIP程序。简言之,染色质为1×106细胞与抗体和蛋白A包裹的磁珠孵育过夜,并在4°C下旋转。珠子在4°C下清洗4次,然后用蛋白酶K(最终浓度为0.4 mg/ml)重新悬浮在Diagenode洗脱缓冲液中,然后在65°C下培养过夜,同时在900 RPM下摇晃。使用MinElute PCR纯化试剂盒(Qiagen 28606),用20μl不含RNAse/DNA酶的ChIP颗粒水洗脱DNA。使用Agilent 2200 TapeStation分析DNA片段大小,并使用Quant-IT dsDNA HS分析试剂盒(Thermo Fisher2010年第3312季度).

简单地说,使用Sybr Green Select Master Mix(Applied Biosystems 4472908)和定量实时PCR对2μl ChIP DNA进行测试,以富集每个抗体的阳性和阴性对照靶基因。以下列出了来自St.Jude的每个靶基因的引物序列。关键资源表中提供了商业来源的引物信息。循环条件:95°C下5 min;40次循环:95°C下30 s,60°C下30s,72°C下30ms;添加熔化曲线步骤。

阳性对照靶基因
蛋白质靶基因或
最近基因
正向底漆反向底漆
H3K4me1型5R3日历年CTGACCTAATGCCC公司
阿卡
GACTCTGGCACCTGT公司
GTC公司
H3K27me3型锌f333tgaagaacatct公司
GCGAACC公司
TCGCGCACTCATACG公司
TTTC公司
H3K9米3ZNF333型TGAAGACACATCT公司
GCGAACC公司
TCGCGCACTCATACG公司
TTTC公司
CTCF公司TRABD2B型ACTGTGGCAGCG公司
ACAGATT公司
TCCTGTCCCTCTCGCTG公司
TTC公司
阴性对照靶基因
蛋白质靶基因或
最近基因
正向底漆反向底漆
H3K4me1型MAG1(刀库1)CGCACACGTTT公司
ACTCC公司
AGGGCGGTTAGGAGGC公司
TTAG公司
H3K36me3型MAG1(刀库1)CGCACACGTTT ACTCCAGGGCGGTTAGGAGGC公司
TTAG公司

文库准备和排序

ChIPSeq文库是根据制造商的说明,使用NEB下一个ChIP-Seq文库Illumina(NEB)制备试剂集制备的,并进行了以下修改:免疫沉淀染色质的量在1到10 ng之间,具体取决于样品。接合器连接DNA的PCR扩增在10-15个周期之间变化,1 ng染色质输入被扩增15个周期,10 ng输入被扩增10个周期。使用0.8X胎圈比执行单安培胎圈清理步骤。未执行辅助尺寸选择。使用单读50周期协议在HiSeq 2000上生成数据。

全基因组测序和分析

如前所述进行全基因组测序和文库构建(Chen等人,2013)具有以下修改;使用Illumina兼容适配器输入250–500 ng基因组DNA用于文库构建,使用Kapa HiFi Hotstart ReadyMix(Kapa Biosystems)进行4–6个周期的扩增。如前所述,对SNV、CNV、SVs和Indels进行鉴定(Chen等人,2013).

RNA测序和分析

通过酚氯仿萃取从单个TRIzol(生命技术)制剂中分离出RNA。样品首先用组织匀浆器(Polytron,PT10–35GT)在17000 RPM下匀浆30 s。然后将1:4体积的氯仿(Sigma-Aldrich)添加到每个样品中,并在室温下培养3分钟,然后在12000下离心然后将水层转移到1.5 mL硅化管中,然后添加2.0μL糖原(罗氏应用科学公司)和500μL异丙醇(费希尔科学公司)。样品在室温下培养10分钟,并在12000下离心然后用冰镇80%乙醇(Fisher Scientific)将样品在4°C下清洗两次,以除去盐分,再悬浮在含有焦碳酸二乙酯的水中,并用NanoDrop分光光度计(Thermo Fisher科学)进行定量。

使用RNeasy Plus Mini Kit(Qiagen)从切除的肿瘤、O-PDX、原始人类成肌细胞和肌管以及胎儿骨骼肌中提取RNA。根据制造商说明(Illumina),使用TruSeq Straded Total RNA Library Prep Kit从~500 ng RNA制备文库。根据制造商说明(Illumina),使用HiSeq 2000或HiSeq 2500测序器进行配对、100周期测序。

使用TruSeq Straded Total RNA Library Prep Kit(Illumina)进行RNA-Seq,输入250 ng-1μg总RNA。如前所述,在RNA-Seq中进行基于基因的FPKM量化(Chen等人,2013).

DNA甲基化和分析

全基因组亚硫酸氢盐测序(WGBS)数据映射到GRCh37-lite(人类)或BSMAP2.74,参数如下:-z 33-f 5-g 3-r 0-m 17-x 600–u,原始甲基化状态由Bis-SNP调用(版本0.82.2)。我们开发了一种基于回归树的方法,使用m值识别样本的区域甲基化模式(=日志2#M(M)e(电子)t吨小时t吨e(电子)d日C类+0.5#U型n个e(电子)t吨小时t吨e(电子)d日C类+0.5). 相邻节段的平均β值差异不显著(p>1E-6,学生t检验)或很小(β=#M(M)e(电子)t吨小时t吨e(电子)d日C类#M(M)e(电子)t吨小时t吨e(电子)d日C类+#U型n个e(电子)t吨小时t吨e(电子)d日C类)(<0.1)递归合并。

收集单个样本的分段断点,生成一个区域的初始列表,用于测试甲基化和基因表达之间的相关性。计算可变甲基化区域(最高和最低平均m值之间的差异>=0.3)与可变表达相邻基因(甲基化区域重叠TSS到TES上游10 kb,最高和最低FPKM之间的差异至少为2倍,最高FPKM>=1)之间的相关性。保留FDR<=0.5的基因/区域对进行合并,合并具有相同方向(正相关或负相关)的相同基因的相邻区域(相距不超过1 kb)。重新评估相关性,并将最终FDR<=0.1的配对用作与基因表达显著相关的差异甲基化区域。

蛋白质组学和磷酸蛋白质组学

多重TMT-LC/LC-MS/MS蛋白质组定量分析

等压标记,如iTRAQ和TMT,正在成为一种强大的策略,用于高通量和重复性的深度多重蛋白质组学分析(Bai等人,2017年). 与该方法相关的一个局限性是,在LC-MS/MS分析过程中,目标肽离子通常与其他共洗脱肽离子共同分离,这会导致高噪声水平压缩定量比率并降低测量准确度。幸运的是,比率压缩效应也缓解了实验变化,因此对蛋白质差异表达分析的影响很小(Tan等人,2017年). 此外,通过广泛的肽分馏、狭窄的分离窗口和MS后校正可以减少比率压缩(Niu等人,2017年). 或者,MS3方法几乎可以消除这种比率压缩,但它需要更长的占空比、特定的MS设置以及使用低分辨率MS2进行鉴定,这通常会影响肽/蛋白质鉴定(McAlister等人,2014年). 为了平衡与等压标记相关的利弊,我们通过长梯度高分辨率LC/LC-MS/MS实现了广泛的分离,以实现深层蛋白质组覆盖,并减少量化过程中的离子压缩。此外,我们使用样本复制来促进差异表达蛋白质的统计推断,这是蛋白质组学分析中广泛采用的策略(Bai等人,2017年;Tan等人,2017年).

RMS O-PDX肿瘤组织、成肌细胞和肌管细胞的蛋白质组和磷酸蛋白质组分析108收集每个样本的成肌细胞和肌管细胞,用10ml冰镇PBS清洗两次,然后在液氮中快速冷冻。这些过程在5分钟内完成。将O-PDX移植的小鼠麻醉并在处死前用10ml PBS灌注,收集O-PDX肿瘤组织的中心切片,匀浆,并在液氮中快速冷冻。这些步骤在10分钟内完成。

蛋白质提取、消化、标记和汇集

蛋白质提取、消化、标记和汇集的操作与前面描述的类似(Bai等人,2017年;Tan等人,2017年). 在新制备的裂解缓冲液(50 mM HEPES,pH 8.5,8 M尿素,0.5%脱氧胆酸钠和磷酸酶抑制剂混合物(PhosSTOP,Roche))中裂解RMS O-PDX组织、成肌细胞和肌管细胞。样品裂解物的蛋白质浓度通过BCA蛋白质分析(赛默飞世尔科学公司)以滴定的BSA为标准进行定量~每个样品1 mg蛋白质首先在室温下用Lys-C(Wako,1:100 w/w)消化2小时,用50 mM HEPES(pH 8.5)稀释4倍,然后在室温下再用胰蛋白酶(Promega,1:100 w/w)进一步消化过夜。添加1%三氟乙酸以冷却消化反应,然后使用Sep-Pak C18筒(Waters)脱盐,并使用speedvac干燥脱盐的肽。然后,将样品重新悬浮在50 mM HEPES中,pH值为8.5,并按照制造商的说明,用10-plex TMT试剂进行标记。最后,将10个等压标记样品以等量汇集在一起,用Sep-Pak C18筒再次脱盐,然后进行speedvac干燥。

离线碱性pH反相液相色谱法

如前所述,在安捷伦1220 LC系统上进行碱性pH反相液相色谱肽预分离(Bai等人,2017年;Tan等人,2017年). 将合并的TMT标记样品溶解在缓冲液A(10 mM甲酸铵,pH值8)中,并在两个XBridge C18柱(3.5μm粒径,4.6 mM×25 cm,Waters)上分离成约180个组分,从15%到65%的缓冲液B(95%乙腈,10 mM甲酸铵,pH 8,流速:0.4 ml/min),梯度为220 min。将每个组合部分的5%进行干燥以进行全蛋白组分析,其余95%通过speedvac进行干燥以用于磷酸蛋白组分析。

TiO富集精制磷酸肽2

按照前面描述的改进方案进行磷酸肽富集(Tan等人,2015年). 简单地说,添加肽来清洁TiO2在结合缓冲液(65%乙腈、2%TFA和0.5 mM KH)中肽与珠重量比为1:4的珠(GL科学)2人事军官4)并培养20分钟。将浓缩的磷酸肽洗涤、洗脱、干燥并溶解在5%甲酸中,以进行LC-MS/MS分析。

长梯度酸性pH反相LC-MS/MS

分析是基于我们的优化平台进行的,如前所述(Bai等人,2017年;Tan等人,2017年). 干燥的肽部分在加载缓冲液(5%甲酸)中进行重组,加载到反相柱(75μm×50 cm,1.9μm C18树脂(德国Maisch GmbH博士))上,并与FUSION或Q Exactive HF质谱仪(Thermo Fisher Scientific)接口。用15–65%梯度的缓冲液B洗脱肽6小时(缓冲液a:0.2%甲酸,5%二甲基亚砜;缓冲液B:a缓冲液加65%乙腈,流速:0.25μl/min)。采用蝶形组合加热器(Phoenix S&T)在65°C下加热色谱柱,以降低背压。质谱仪以数据相关模式运行,MS1设置为60000分辨率,1×106AGC目标和最大50 ms离子时间,前20 ms/ms高分辨率扫描,MS2设置为1米/z隔离窗偏移0.2,分辨率60000,最大离子时间100ms,1×105AGC靶点、HCD、33归一化碰撞能量和40秒动态排斥(磷酸蛋白质组的35归一化冲突能量和20秒动态排斥)。

基于标签的混合搜索引擎JUMP识别肽

使用我们最近开发的JUMP搜索引擎进行肽鉴定,提高了灵敏度和特异性(Wang等人,2014年). 商用数据库搜索引擎可分为两类:基于标记的从头开始测序(例如灵敏度有限的峰值)和基于模式的数据库搜索(例如SEQUEST、MASCOT)。JUMP软件集成了这两种方法来对推测的肽进行评分,与这些商用工具相比,显示出了显著的改进(Wang等人,2014年). JUMP软件以前曾在许多出版物中使用过(Gong等人,2017). 分析与前面描述的类似(Tan等人,2017年)首先将MS/MS原始文件转换为mzXML格式,并根据复合目标/诱饵数据库FDR估计进行搜索。靶蛋白数据库由Uniprot小鼠和人类数据库编译而成(人类数据库:88965个蛋白质条目;小鼠数据库:52738个蛋白质条目,2015年2月下载),诱饵数据库通过反转靶蛋白序列生成。以前体离子和产物离子的±10ppm质量耐受性搜索光谱,具有完全胰蛋白酶限制,N-末端TMT标签和赖氨酸的静态修饰(+229.16293),丝氨酸、苏氨酸和酪氨酸的动态修饰(+79.996633,用于磷酸蛋白质组分析),三个最大修饰位点,两个最大遗漏劈理,以及,b条、和离子。肽谱匹配(PSM)首先通过MS质量精度(~2ppm,±4个标准偏差)进行过滤。在过滤之前,将JUMP Jscore大于30的双电荷肽PSM用于整体质量重新校准。存活的PSM首先按前体离子电荷状态分组,然后再按Jscore和dJn值进一步过滤。对这些值应用截止值,并进行调整,直到达到蛋白质FDR<1%。如果一个肽由多个蛋白质共享,则根据简约规则,PSM最高的蛋白质将代表该肽(Nesvizhskii和Aebersold,2005年).

Lscore从JUMP软件套件中分配磷灰石

为了确定磷酸化位点在肽上定位的可靠性,我们采用了磷酸化RS算法的概念(Taus等人,2011年)与前面描述的相同,计算每个PSM中的磷脂酶定位得分(Lscore,0–100%)(Tan等人,2017年). 将磷脂酶与蛋白质序列对齐,以生成除PSM Lscore之外的蛋白质水平Lscore。如果针对一个特定的磷酸位点鉴定出多个PSM,则蛋白质Lscore由最高的PSM Lscore表示。由于随机分配含有不明确磷酸位点的PSM通常会导致蛋白质上不可靠磷酸位点的数量过多,因此我们实施了一系列规则来缓解这个问题:(i)如果一个PSM中单个磷酸化肽的第一个和第二个位点之间的PSM L分数差距>10%,则选择L分数最高的位点;(ii)否则,我们检查相应蛋白质中的磷酸化位点,而不是选择具有最高蛋白质Lscore的位点。这允许低质量PSM借用从高质量PSM推断出的信息;(iii)如果PSM和蛋白质水平Lscore都无法区分,则根据发生顺序为磷酸化位点分配启发式优先级:SP-motif、S、T和Y;(iv)如果PSM不满足上述任何规则,则首先通过JUMP Jscores对这些PSM进行分类,然后我们选择由其他高Jscoress的PSM确定的蛋白质磷酸位点。

使用JUMP软件套件进行基于TMT的蛋白质和磷酸化位点定量

该分析按照与之前报告类似的步骤进行(Tan等人,2017年). (i) 提取每个PSM的TMT报告离子强度;(ii)根据每个标记试剂的同位素分布校正原始强度(Tan等人,2017年); (iii)排除了报告离子强度极低的PSM(例如,最小强度<1000,中间强度<5000);(iv)通过归一化所有PSM的修剪中值强度,校正样本加载偏差;(v) 计算了样本中以平均值为中心的强度(Tan等人,2017年); (vi)通过对相关PSM求平均值来总结蛋白质或磷脂酶的相对强度;(vii)通过将相对强度乘以前三个最丰富的PSM的grand-mean强度,得出蛋白质或磷脂酶的绝对强度。为了生成多批次的组合定量表,每个批次中都包含一个普通样品(SJRHB10468_X1)作为内部标准。根据该内部标准,对不同批次的MS强度进行标准化。

蛋白质组和磷酸蛋白质组的差异表达分析

使用LIMMA R软件包对整个蛋白质组和磷酸蛋白质组进行差异表达分析,该软件用于分析基因表达,同时比较多个基因靶点。LIMMA通过拟合线性模型跨基因借用信息,以克服小样本和复杂实验设计的问题,因此它是基于TMT的深层蛋白质组数据差异表达分析的理想工具。简言之,(i)线性模型适用于每个蛋白质或磷脂酶的表达数据;(ii)使用经验贝叶斯方法跨基因借用信息;(iii)通过Benjamin Hochberg方法调整p值;(iv)然后应用调整后的p值截止值0.05;(v) 在蛋白质组和磷酸蛋白质组的至少一组比较中,通过1.5倍和2.0倍的变化(相当于肿瘤生物复制品的3倍标准偏差)进一步过滤剩余蛋白质。

加权基因共表达网络分析(WGCNA)和通路注释

分析由WGCNA R软件包完成,与前面描述的类似(Tan等人,2017年). 只有DE蛋白和磷酸化位点分别用于定义蛋白质组和磷酸蛋白质组共表达簇(即WPC和PPC)。(i) Pearson相关矩阵是通过计算蛋白质之间的相关性生成的(只考虑了正相关性),并使用无标度拓扑原理计算邻接矩阵,将其进一步提高到16次方。(ii)然后通过混合动态树切割法确定共表达簇(Langfelder等人,2008年)具有用于合并模块的高度截断(例如0.2);(iii)基于每个共表达簇的第一主成分计算一致趋势(特征基因)。然后将蛋白质分配到相关性最高的共表达簇;(iv)然后使用从MsigDB下载的Hallmark通路数据库对每个共表达簇进行注释(Liberzon等人,2011年)(Hallmark数据库中未注释的肌源性调控途径是从KEGG数据库中手动提取并添加的),通过Fisher精确测试。在每个共表达簇中,选择B.H.调节p值小于0.05的途径作为解除调控的途径。

蛋白质组学的免疫印迹验证

异种组织在含有蛋白酶抑制剂的RIPA裂解缓冲液(Sigma,目录R0278)中进行裂解(Halt Inhibitor Cocktail,目录87785)。蛋白质在4–15%SDS-PAGE凝胶中溶解并转移到PVDF膜(Millipore)。在Odyssey阻断缓冲液(Licor,目录927-4000)中阻断膜,然后探测以下主要抗体:NOS1(1:1000,细胞信号,目录4234)、HMGA2(1:1000;细胞信号,类别5269)、RPS6的荧光-Ser235/236(1:1000、细胞信号,目2211)、CDC25,磷酸CDC25C(Ser216)(1:1000,细胞信号,目录9528)和肌动蛋白(1:5000,Sigma,AC-15)。次要抗体为抗兔抗体(Licor,目录926–32221);防鼠(Licor,目录926–32210)。

转录靶点分析

转录靶点来源于文献,补充信息中提供了单个靶点列表。然后对每个样本中每个目标的FPKM值进行平均,并将其归一化为成肌细胞。为了纳入分析,至少一个正常样品的FPKM必须高于1.0。

免疫组织化学

通过一些组织学特征比较患者H&E染色玻片和O-PDX衍生玻片的组织病理学特征的相似性。这些患者-异种移植物对通过以下属性进行评估:总体肿瘤细胞数、生长模式、细胞形态(包括多形性程度)和有丝分裂活性。

使用福尔马林固定石蜡包埋组织(厚度4μm)进行染色,并评估适当的阳性和阴性对照。有关所用抗体的信息如下:

通过免疫组织化学方法对所有患者样本和异种移植组织进行MYF5、CDK4、CDK6、CCND1、CCND2、CCND3、pRB、p16和肌生成素的评估。总共回顾了34张苏木精和伊红染色切片和306张免疫组化染色切片。抗体信息可在下文和关键资源表中找到。

抗体在免疫组织化学分析中的应用
抗体验证染色剂主机稀释检索放大/
舞台调度
试剂
供应商/
猫#
MYF5型鲁奥过激的马来西亚令吉1:25第2页–第92页
最小值
abcam;
ab125078号
肌生成素体外诊断过激的MP公司1:500CC1–64号机组
最小值
Dako;
M3559型
第16页体外诊断过激的MM(毫米)远程终端设备CC1–48号机组
最小值
罗氏;
705–4713
pRB(pRB)鲁奥过激的RP公司1:100CC1–64号机组
最小值
单元格
信号;
9301秒
CCND3号机组鲁奥奥姆尼斯MM(毫米)1:400高pH值
TR 30型
最小值
鼠标链接器-
10分钟
单元格
信号传导;
2936
CCND2号机组鲁奥债券-
RP公司1:100ER1–20
最小值
圣诞老人
克鲁兹;供应链-
593
CCND1号机组体外诊断过激的马来西亚令吉远程终端设备CC1–36号机组
最小值
单元格
品牌:
241R-18型
川东北6过激的MM(毫米)1:25CC1–64分钟圣诞老人
克鲁兹;
sc7961型
川东北4鲁奥过激的MM(毫米)1:50CC1–64号机组
最小值
AB公司
布洛克(罗氏)
圣诞老人
克鲁兹;
sc260标准

体外诊断-体外试验诊断,RUO-仅供研究使用,Bond-III-徕卡,Omnis-Dako,Ultra-Ventana BenchMark Ultra(Roche),MP-小鼠多克隆,RM-兔单克隆,MM-小鼠单克隆,RP-兔多克隆,CC1-细胞调节溶液1,TR-热回收,RTU-随时可用

评估抗体的染色强度,从无反应到3+。记录具有预期染色定位的肿瘤细胞百分比。这仅限于MYF5、CCND1、CCND2、CCND3、pRB、p16和肌生成素的核染色。CDK4和CDK6的细胞质或细胞核染色均可接受。根据预期的染色定位,评估患者和异种移植物样本之间的不一致性。当出现完全差异(一个样本阳性,另一个样本阴性或染色<5%)时,或当两个样本均阳性但样本之间存在显著差异(>50%)时,出现不一致。如果染色完全失败,则认为比较不适用。

普通药物筛选

对实体瘤异种移植物和细胞系进行药物筛选,并从儿童实体瘤网络获得数据(www.stjude.org/CSTN网站) (Stewart等人,2017年). 对所有异种移植物系进行最佳生长条件分析,包括电镀密度、DMS敏感性和阳性对照选择,如前所述(Stewart等人,2017年)在进行药物敏感性验证和分析之前。将所有RMS异种移植物系在SkBM-2细胞培养基(Lonza Cat#CC-3246)的涂胶板上进行体外培养,并补充SingleQuots补充剂(Lonsa Cat#CC3244)。在验证了分析条件后,可以对异种移植细胞系进行活性筛选。将分离的异种移植物细胞以先前确定的最佳密度在384孔板中进行电镀。电镀后20小时,将分析板与感兴趣的化合物分三份进行药物处理。对于剂量反应给药,每种化合物在每个板上以不同浓度对多个孔进行给药。给药72小时后,RLU由CellTiter-Glo(CTG)/EnVision系统确定。然后对数据进行分析,以确定单点用药情况下的活性化合物以及剂量-反应曲线。

ganetespib联合药物筛选

原始数据处理

从重复实验中汇集每个浓度下每个化合物的原始发光RLU(相对光单位)值,进行log2转换,并使用以下方程进行归一化以获得%活性:100*[(negctrl–化合物)/(negctrl–posctrl)]。

基于活性百分比的剂量反应分析

使用数字版权所有(丽兹酒店和斯特里比格酒店,2005年)R中的包。四参数S形函数LL2.4型使用的约束条件如下:−10≤山坡≤0;0≤y0≤最大中值归一化%活度;0≤yFin≤最大中值归一化%放射性;10−10≤ec50≤ 10−4(相当于本实验中测试的药物浓度范围)。对每种化合物测试的最高浓度值进行10%的加权,以减少曲线拟合伪影。

如果无法拟合S形剂量反应曲线平滑样条线R中的选项用于拟合曲线,该曲线可用于确定曲线下的面积。

联合Ganetespib药物的AUC分析

获得了每种药物与二甲基亚砜和4种浓度的Ganetespib(3 nM、10 nM、33 nM和100 nM)组合的AUC值。δAUC值是Ganetespib当前浓度下的AUC与仅含DMSO的AUC之差。较高的值意味着在Ganetespib的存在下对增殖的抑制作用更强。

选择药物组合的响应面模型分析

使用BRAID联合药物作用模型对七种药物进行联合研究(Twarog等人,2016年). 如前所述,在96孔板中以每孔1000个细胞的速度电镀RD细胞,并在8天后使用CTG和envision进行读取(Stewart等人,2017年). 在不同的日子进行生物复制。从重复实验中汇集每种浓度下每种化合物的原始发光RLU(相对光单位)值,转换log10,并通过减去DMSO对照孔的值进行归一化。

产生95%置信区间的药物组合的BRAID kappa参数用于指示相互作用的类型:如果kappa_l和kappa_ h结合0=相加性;如果kappa_l和kappa_ h均为正=协同作用;如果kappa_l和kappa_ h都为负=对抗。90%疗效时的可达到疗效指数(IAE90)与二维治疗窗口相当,可以直接比较药物组合的疗效:IAE90越高,联合疗效越高(Twarog等人,2016年).

药动学试验

帕尔博西利布

在单次口服16 mg/kg帕尔博西布林游离碱当量后,评估雌性裸鼠(Charles River Laboratories,12-16周龄)的帕尔博西布林总血浆和肿瘤PK。将盐酸帕尔博西利(MedKoo,批号SSC40317,纯度99%)悬浮于1%甲基纤维素(MC,类型400 cPs)和1%吐温80(1%MCT)超纯水中,最终标称浓度为1.6 mg/mL游离碱当量,灌胃体积为10 mL/kg。使用IACUC批准的方法在给药后10分钟、1小时、4小时、8小时和24小时处死小鼠,每个时间点3只小鼠。通过心脏穿刺将全血收集到含有10μL肝素钠的塑料微型离心管中,旋涡混合抗凝剂,立即离心至血浆中,并保存在干冰上,以备研究剩余时间之用。然后通过右心室用无钙和无镁的PBS灌注小鼠,切除原位异种移植物,用PBS冲洗,称重,分成等分试样,并置于干冰上。在体内程序结束时,从干冰中转移所有样品,并将其置于−80°C的温度下,直至分析。

使用灵敏且特异的液相色谱法和串联质谱法评估血浆和肿瘤匀浆中帕波西卜的总浓度。首先,用5:1体积的超纯水稀释肿瘤样本,并用基于珠的技术进行均质(Liang等人,2011年)在FastPrep-24系统上(MP Biomedicals,加州圣安娜)。将1.4 mm陶瓷球(MP Biomedicals,Lysing matrix D,10 mg/mg肿瘤)添加到含有样品的微离心管中。然后在FastPrep-24系统上对样品进行三次6.0 M/S振动循环,每次1分钟。为了防止摩擦导致过热,在每个循环之间将样品放置在湿冰上5分钟。然后将匀浆储存在−80°C下,直至分析。

在乙腈中制备盐酸帕尔博西利(MedKoo,批号SSC40317,纯度99%)储备溶液,并用于加标基质校准品和质量控制。血浆和肿瘤匀浆样品(各25μL)用75μL的10 ng/mL LEE011(Cayman Chemical Co.,Lot 0467691–2,纯度≥95%)在乙腈中进行蛋白质沉淀,作为内部标准。通过LEAP CTC PAL自动取样器将2μL等分提取的上清液注入岛津LC-20ADXR高效液相色谱系统。使用Phenomenex Gemini C6 Phenyl(3.0μm,30 mm x 2.0 mm)色谱柱进行液相色谱分离,色谱柱保持在60°C,梯度洗脱,流速为0.50 mL/min。二元流动相由储液罐a中pH为6.0(90:10 v/v)的水200 mm醋酸铵和pH为6.0的乙腈水-200 mm醋酸铵盐组成初始流动相由25%的B组成,在三分钟内线性增加到65%的B。然后在100%的B下冲洗柱两分钟,然后在初始条件下平衡两分钟,总运行时间为七分钟。在这些条件下,分析物和IS分别在1.36和1.04分钟洗脱。

使用SCIEX API 5500 Q-TRAP在正ESI模式下用串联质谱法检测分析物和IS,监测以下质量变化:Palbociclib 448.25->380.10,LEE011 435.26->322.20。

方法鉴定和生物分析运行均通过了P-PKSR的非GLP分析性能验收标准。二次模型(1/X2称重)在1至500 ng/mL范围内拟合校准品,相关系数(R)≥0.9951。定量下限(LLOQ)定义为峰面积信噪比为5或更大,与IS的基质空白相比为1 ng/mL。批内精密度和准确度分别为<14.9%CV和103%-114%。

根据矩阵和时间点对所得的帕尔博西布林浓度-时间(Ct)数据进行分组,手动插补低于定量下限(BLOQ)的数据如下:如果在任何时间点≥2/3rds的Ct结果高于LLOQ,则BLOQ数据替换为½LLOQ值,ELSE将整个时间点的数据视为缺失。然后,使用Phoenix WinNonlin 6.4(Certara USA,Inc.,Princeton,NJ)生成Ct数据汇总统计数据(算术平均值、标准偏差、%CV、最小值、中位数、最大值),并对每个基质的palbociclib算术平均Ct数据进行非部分药代动力学分析(NCA)。应用血管外模型(模型202),并使用“线性上升-下降”梯形规则估计Ct曲线下面积(AUC)值。终相被定义为Ct曲线末端的三个时间点,使用终相的未加权对数线性回归估计消除速率常数(Ke)。终端消除半衰期(T1/2)估计为0.693/Ke,从时间0到无穷大的AUC(AUCinf)估计为到最后一个时间点的AUC+预测的碎屑/Ke。

估计的其他NCA参数包括观察到的最大浓度(Cmax)、Cmax时间(Tmax),最后一个观察到的时间点的浓度(Clast)、Clast时间(Tlast)、表观清除率(CL/F=剂量/AUCinf)和表观终末分布体积(Vz/F)。给药间隔(Cavg)的平均浓度估计为AUCinf/给药间隔小时数。palbociclib从血浆到相关组织的表观分配系数(Kp,组织)估计为可用时AUCinf、组织与AUCinf血浆的比值。为了估计与临床相关的小鼠剂量,将合成的小鼠血浆AUCinf与报告的人类血浆PK值进行比较,假设单药帕尔博西利批准剂量为125 mg PO QD。所有的推断都是在小鼠和人类的非时间依赖性、线性和剂量比例PK假设下进行的。

考虑到系列牺牲设计,Palbociclib在研究中表现出中度PK变异性。血浆浓度变异系数(CV)为13.5%至58.4%,肿瘤浓度变异系数为25.3%至68.5%。最大的变化发生在第一个采样时间点10分钟,即在吸收阶段。表观血浆清除率(CL/F)估计值为5.47 L/hr/kg,与之前发表的文献中得出的结果相似(4.0至5.9 L/hr/kg])(Parrish等人,2015年). 与文献相比,表观终末半衰期(T1/2)稍短(2.61小时vs.约3.4小时)。Palbociclib在原位RMS异种移植物中分布良好,达到了比血浆高约6倍的浓度。肿瘤中的T1/2比血浆长1.8倍。NCA参数估计值和按矩阵分组的平均值(SD)palbociclib Ct图如下所示:

美国食品药品监督管理局批准的帕博西立普普通剂量为125 mg PO QD,持续21天,然后休息7天。以下PK和ADME信息来自帕博西立普美国食品药品监督管理局NDA审查。帕博昔单抗的群体PK分析显示,典型患者的平均血浆CL/F为60.2 L/小时,导致稳态时估计的总血浆AUCtau为2080小时ng/mL。据报道,患者体内帕尔博西利的血浆半衰期为22至29小时。此外,在两项临床研究中多次服用帕尔博西利125 mg PO QD后,平均血浆总Cmax范围为94.9至116 ng/mL,平均Ctrough为61 ng/mL。通过标准计算,患者的估计血浆总Cavg,ss约为85 ng/mL。由于小鼠和人类之间的血浆蛋白结合没有明显差异(Fu分别为0.159和0.147),因此可以使用血浆总浓度来估算小鼠等效剂量(MED)。假设小鼠的剂量和时间线性PK和最小研究间PK变异性,并基于人类血浆总AUCinf和AUCtau,MED估计为每天10 mg/kg帕博西利游离碱当量,该配方在类似的研究条件下。同样,小鼠血浆Cavg支持这一剂量,后估计的Cavg为76.9 ng/mL,与观察到的患者Cavg ss为85 ng/mL相似。

阿贝马奇利布

单次口服50 mg/kg阿贝马昔布游离碱当量后,评估雌性裸鼠(Charles River Laboratories,8–16周龄)的阿贝马昔布总血浆和肿瘤PK。将甲磺酸阿贝马西利(LY2835219,Abmole,M2112,Lot 2,纯度>98%)悬浮于1%羟乙基纤维素(HEC)、0.25%吐温80和~0.05%西米硫醇的超纯水中,最终标称浓度为5 mg/mL游离碱当量,灌胃体积为10 mL/kg。使用IACUC批准的方法在给药后10分钟、1小时、4小时、8小时和24小时处死小鼠,每个时间点3只小鼠。通过心脏穿刺收集全血至Sarstedt Microvette®500μl K3 EDTA微型离心管,旋涡混合抗凝剂,立即离心至血浆中,并储存在干冰上,以备剩余研究之用。然后通过主动脉用PBS灌注小鼠,切除原位异种移植物,用PBS冲洗,并放在干冰上。在体内程序结束时,从干冰中转移所有样品,并将其置于−80°C的温度下,直至分析。

使用灵敏且特异的液相色谱法和串联质谱法评估血浆和肿瘤匀浆Abemacilib的总浓度。首先,用体积为5:1的超纯水稀释肿瘤样本,并用基于珠的技术进行均质(Liang等人,2011年)在FastPrep-24系统上(MP Biomedicals,加州圣安娜)。将1.4 mm陶瓷球(MP Biomedicals,Lysing matrix D,10 mg/mg肿瘤)添加到含有样品的微离心管中。然后在FastPrep-24系统上对样品进行三次60 M/S的振动循环,每次1分钟。为了防止摩擦导致过热,在每个循环之间将样品放置在湿冰上5分钟。然后将匀浆储存在−80°C下,直至分析。

在乙腈中制备甲磺酸阿贝马西布(LY2835219,Abmole,M2112,Lot 2,纯度>98%)储备溶液,校正盐含量,并用于加标基质校准品和质量控制。血浆和肿瘤匀浆样品(各25μL)用100μL 150 ng/mL帕尔博西利(LC Labs,P-7788,Lot PLH-103,纯度>99%)在乙腈中沉淀蛋白作为内标。通过LEAP CTC PAL自动取样器将2μL等分提取的上清液注入岛津LC-20ADXR高效液相色谱系统。使用Phenomenex Kinetex 2.6μm EVO C18(100μm,50×2.1 mm)在50°C下进行LC分离,梯度洗脱,流速为0.5 mL/min。二元流动相由20 mm醋酸铵在H中组成2储液罐A中的O和储液罐B中的甲醇:乙腈:异丙醇(40:30:30 v/v)。初始流动相保持在25%B下0.5 min,然后在2.5 min内线性增加到100%B。然后在100%B下冲洗柱2 min,然后再在初始条件下平衡2 min,总运行时间为7 min。在这些条件下,分析物和IS分别在1.37和0.98 min洗脱。

使用SCIEX API 5500 Q-TRAP在正ESI模式下用串联质谱法检测分析物和IS,并监测以下质量转变:abemaciclib 507.28>393.20,palbociclib 448.25->380.10。

实验生物分析运行均被发现可用于单一非GLP临床前PK评估。线性模型(1/X2称重)在5.00至500 ng/mL范围内拟合校准品,相关系数(R)≥0.9967。定量下限(LLOQ)定义为峰面积信噪比为5或更大,与IS的基质空白相比为5.00 ng/mL。整个基质的批内精密度和准确度分别为<15.6%CV和90.7%至112%。

生成的阿贝马西布林浓度-时间(Ct)数据按矩阵和时间点分组,低于定量下限(BLOQ)的数据手动插补如下:如果在任何时间点≥2/3rds的Ct结果高于LLOQ,则BLOQ数据替换为½LLOQ值,ELSE将整个时间点的数据视为缺失。然后,使用Phoenix WinNonlin 6.4(Certara USA,Inc.,Princeton,NJ),生成Ct数据汇总统计数据(算术平均值、标准差、%CV、最小值、中值、最大值),并对每个矩阵的阿贝马昔单抗算术平均Ct数据进行非部门药代动力学分析(NCA)。应用血管外模型(模型202),并使用“线性上升-下降”梯形规则估计Ct曲线下面积(AUC)值。终相被定义为Ct曲线末端的三个时间点,使用终相的未加权对数线性回归估计消除速率常数(Ke)。终末消除半衰期(T1/2)估计为0.693/Ke,从时间0到无穷大的AUC(AUCinf)估计为到最后一个时间点的AUC(AUClast)+预测的Clast/Ke。

估计的其他NCA参数包括观察到的最大浓度(Cmax)、Cmax的时间(Tmax)、最后一个观察到的时间点的浓度(Clast)、Clast的时间(Tlast)、表观清除率(CL/F=剂量/AUCinf)和表观终末分布体积(Vz/F)。给药间隔(Cavg)的平均浓度估计为AUCinf/给药间隔小时数。阿贝马西利从血浆到相关组织的表观分配系数(Kp,组织)估计为可用时AUCinf、组织与AUCinf血浆的比值。为了估计与临床相关的小鼠剂量,将合成的小鼠血浆AUCinf与报告的人类血浆PK值进行比较,假设单药阿贝马西利最大耐受剂量(MTD)为200 mg PO BID(Patnaik等人,2016年). 所有的推断都是在小鼠和人类的非时间依赖性、线性和剂量比例PK假设下进行的。

Abemacilib浓度在小鼠血浆中表现出较低的变异性,在采样时间点上的变异系数为4.1%至40.1%,在血浆和肿瘤中的0.167小时时间点的变异性最高。所有血浆和肿瘤结果均高于LLOQ。肿瘤穿透似乎相对较快且广泛,基于AUCinf,Kp肿瘤值为6.20。阿贝马西利在血浆中的T1/2为5.2小时,在肿瘤中观察到类似的半衰期(5.98小时)。值得注意的是,我们在50 mg/kg PO条件下的阿贝马西利血浆PK与Tate报告的不同(Tate等人,2014年). 我们的小鼠在给药后4小时内显示出类似的血浆Ct曲线,但随后浓度迅速下降,导致AUC降低约3.8倍。尽管使用了甲磺酸盐形式和相同的混悬剂配方,但我们的小鼠并未表现出这些研究人员所描述的长期吸收和低CL/F。NCA参数估计值和平均值(SD)abemacilib Ct图(按矩阵分组)如下所示:

在一项持续服用阿贝马西利单药的第1阶段研究中,MTD为200 mg PO BID。稳态下人体总血浆AUCtau的几何平均值为3000 hr-ng/mL,相应的几何平均谷浓度为197 ng/mL。假设小鼠的剂量比例和线性PK过程,小鼠的等效剂量(MED)为提供相似的总血浆AUC为14.4 mg/kg。测定小鼠和人类血浆中未结合的平均体外分数(Fu,p),得出值分别为0.054和0.027(Raub等人,2015年). 由于这些值处于两倍阈值,定义了临床前研究的显著差异,因此尚不清楚是否应根据两个物种之间的血浆蛋白结合调整阿贝马西布剂量。因此,根据未结合和总血浆AUC,四舍五入和推荐的MED范围分别为7.5至15 mg/kg PO BID。此范围内的剂量还应提供与200 mg PO BID时人类类似的未结合稳态血浆浓度。可能需要进行更多的PK研究,以确定我们的血浆PK与Tate报告的差异的来源,或进一步确认我们的PK发现。

曲美替尼

在单次口服0.1 mg/kg曲美替尼游离碱当量后,评估雌性裸鼠(Charles River Laboratories,12-16周龄)的曲美替尼总血浆和肿瘤PK。将曲美替尼游离碱(Abmole,批号M1759,纯度>99%)悬浮于5%二甲基亚砜(DMSO)、10%克雷莫波尔RH40和10%PEG400中75%超纯水中,最终标称浓度为0.02 mg/mL游离碱当量,灌胃体积为5 mL/kg。使用IACUC批准的方法在给药后30分钟、1小时、4小时、8小时和16小时处死小鼠,每个时间点3只小鼠。通过心脏穿刺收集全血至Sarstedt Microvette®500μl K3 EDTA微型离心管,旋涡混合抗凝剂,立即离心至血浆中,并储存在干冰上,以备剩余研究之用。然后通过右心室向小鼠灌注无钙和无镁PBS,切除原位异种移植物,用PBS冲洗,称重,分成小份,放在干冰上。在体内程序结束时,从干冰中转移所有样品,并将其置于−80°C的温度下,直至分析。

采用灵敏、特异的液相色谱法和串联质谱法评估曲美替尼血浆和肿瘤匀浆的总浓度。首先,用体积为5:1的超纯水稀释肿瘤样本,并用基于珠的技术进行均质(Liang等人,2011年)在FastPrep-24系统上(MP Biomedicals,加州圣安娜)。将1.4 mm陶瓷球(MP Biomedicals,Lysing matrix D,10 mg/mg肿瘤)添加到含有样品的微离心管中。然后在FastPrep-24系统上对样品进行三次6.0M/S振动循环,每次1分钟。为了防止摩擦导致过热,在每个循环之间将样品放置在湿冰上5分钟。然后将匀浆储存在−80°C下,直至分析。

在乙腈中制备曲美替尼游离碱(Abmole,Lot M1759,纯度>99%)储备溶液,并用于刺入基质校准器和质量控制。血浆和肿瘤匀浆样品各25μL,用75μL的75 ng/mL雷美替尼(开曼化学公司,批号0464258-2,纯度≥95%)在乙腈中作为内标进行蛋白质沉淀。通过LEAP CTC PAL自动取样器将2μL等分提取的上清液注入岛津LC-20ADXR高效液相色谱系统。使用Phenomenex Gemini C6 Phenyl(3.0μm,30 mm x 2.0 mm)色谱柱进行液相色谱分离,色谱柱保持在60°C,梯度洗脱,流速为0.50 mL/min。二元流动相由储液罐a中pH为6.0(90:10 v/v)的水200 mm醋酸铵和pH为6.0的乙腈水-200 mm醋酸铵盐组成初始流动相由25%的B组成,在三分钟内线性增加到65%的B。然后在100%的B下冲洗柱两分钟,然后在初始条件下平衡两分钟,总运行时间为七分钟。在这些条件下,分析物和IS分别在2.91和2.8分钟洗脱。

使用SCIEX API 5500 Q-TRAP在正ESI模式下用串联质谱法检测分析物和IS,并监测以下质量转换:曲美替尼616.04->491.10,雷美替尼572.91->394.00。

实验生物分析运行均被发现可用于单一非GLP临床前PK评估。线性模型(1/X2称重)在1.00至500 ng/mL范围内拟合校准品,相关系数(R)≥0.9985。定量下限(LLOQ)定义为峰面积信噪比为5或更大,与IS的基质空白相比为1.00 ng/mL。整个基质的批内精密度和准确度分别为<5.73%CV和95.8%至103%。

将生成的曲美替尼浓度-时间(Ct)数据按矩阵和时间点分组,手动插补低于定量下限(BLOQ)的数据如下:如果在任何时间点≥2/3的Ct结果高于LLOQ,则BLOQ数据替换为½LLOQ值,ELSE将整个时间点的数据视为缺失。然后,使用Phoenix WinNonlin 6.4(Certara USA,Inc.,Princeton,NJ)生成Ct数据汇总统计数据(算术平均值、标准偏差、%CV、最小值、中位数、最大值),并对每个基质的曲美替尼算术平均Ct数据进行非部门药代动力学分析(NCA)。应用血管外模型(模型202),并使用“线性上升-下降”梯形规则估计Ct曲线下面积(AUC)值。终相被定义为Ct曲线末端的三个时间点,使用终相的未加权对数线性回归估计消除速率常数(Ke)。终端消除半衰期(T1/2)估计为0.693/Ke,从时间0到无穷大的AUC(AUCinf)估计为到最后一个时间点的AUC+预测的碎屑/Ke。

估计的其他NCA参数包括观察到的最大浓度(Cmax)、Cmax时间(Tmax),最后一个观察到的时间点的浓度(Clast)、Clast时间(Tlast)、表观清除率(CL/F=剂量/AUCinf)和表观终末分布体积(Vz/F)。给药间隔(Cavg)的平均浓度估计为AUCinf/给药间隔小时数。曲美替尼从血浆到相关组织的表观分配系数(Kp,组织)估计为可用时AUCinf、组织与AUCinf血浆的比值。为了估计与临床相关的小鼠剂量,在假定的单药曲美替尼最大耐受剂量(MTD)为2 mg PO QD时,将生成的小鼠血浆AUCinf与报告的人类血浆PK值进行比较。所有的推断都是在小鼠和人类的非时间依赖性、线性和剂量比例PK假设下进行的。

曲美替尼在研究中的PK变异性较低,浓度变异系数(CV)在5.57%至26.4%之间。在采样期间,肿瘤中没有足够的终相,因此肿瘤终相参数无法估计。利用AUClast、肿瘤与AUClast血浆a Kp的比值,估计肿瘤为5.49。NCA参数估计值和平均值(SD)曲美替尼Ct图(按矩阵分组)如下所示:

FDA批准的单药曲美替尼的常用剂量为每天一次性服用2毫克PO。曲美替尼的群体PK分析(Ouellet等人,2016年)结果显示,典型女性患者的平均血浆CL/F为4.91 L/hr,典型男性患者为6.19 L/hr,导致不同性别的平均AUC为365 hr-ng/mL。据报道,患者中曲美替尼的血浆半衰期为92.5至115小时。在每天一次性多次给药2 mg曲美替尼可稳定后,患者的平均血浆总Cmax范围为17.9至21.2 ng/mL,平均Ctrough范围为11.7至15.0 ng/mL。根据标准计算,患者的估计血浆总Cavg,ss约为15.1 ng/mL。小鼠和人类血浆中未结合的曲美替尼组分(Fu,p)分别为0.05和0.03。通过Fu,p对每个物种的AUC和Cavg的关键PK参数进行调整,在指定的共溶剂溶液中,0.1 mg/kg PO QD的MED在临床上与人类的2 mg PO QD相当。

AZD1775型

之前完成了AZD1775的小鼠等效剂量,小鼠等效剂量确定为60 mg/kg,每天两次给药(Stewart等人,2017).

甘奈特斯比

通过尾静脉缓慢静脉注射(IV)150 mg/kg后,评估雌性裸鼠(Charles River Laboratories,年龄8-16周)中ganetespib的总血浆和肿瘤PK。将Ganetespib(MedChem Express,HY-15205,批号08249,纯度>98%,摩尔Wt 364.40)溶于二甲基亚砜中,并用Kolliphor RH 40(Sigma)在5%注射用葡萄糖中进一步稀释,USP(D5W,Baxter)在10%DMSO/18%Kolliphor RH 40/72%D5W和10mL/kg注射体积中最终标称浓度为15mg/mL。使用IACUC批准的方法在给药后10分钟、1小时、4小时、8小时和16小时处死小鼠,每个时间点处死3只小鼠。通过心脏穿刺用肝素钠收集全血,立即离心至血浆中,并储存在干冰上,以供研究剩余时间使用。然后通过主动脉用PBS灌注小鼠,切除原位异种移植物,用PBS冲洗,并放在干冰上。在体内程序结束时,从干冰中转移所有样品,并将其置于−80°C的温度下,直至分析。

使用灵敏和特异的液相色谱、串联质谱(LC-MS/MS)分析法评估了甘奈特氏菌血浆和肿瘤匀浆的总浓度。首先,组织样品被浸渍,用5:1体积的超纯水稀释,并用基于珠的技术均质(Liang等人,2011年)在FastPrep-24系统上(MP Biomedicals,加州圣安娜)。将陶瓷溶珠基质(MP Biomedicals,Ceramic Bead lysing matrix,3mg/mg组织)添加到含有肿瘤样本的微离心管中。然后在FastPrep-24系统上对样品进行三次60 M/S的振动循环,每次1分钟。为了防止摩擦导致过热,在每个循环之间将样品放置在湿冰上5分钟。然后将匀浆储存在−80°C下,直至分析。

Ganetespib储备溶液在100%甲醇中制备,用于添加基质校准品和质量控制。血浆和组织匀浆样品各50μL,用10μL内标物(SLV320,Tocris Bioscience,纯度>99%)100 ng/mL加标溶液处理,然后使用甲基叔丁基醚进行液-液萃取程序,涡旋并在4°C下离心,上清液在CentriVap离心真空浓缩器(Labconco)中蒸发至干燥。然后用100μL 100%乙腈重新配制提取物,并通过岛津SIL-20ACXR自动取样器将5μL等分的重新配制提取物注入岛津LC20ADXR高效液相色谱系统。使用Phenomenex Luna C8 80 Au LC柱(4.0μm,30 mm x 2 mm)在室温下进行LC分离,梯度洗脱,流速为0.25 mL/min。二元流动相由储液罐a中超纯水中的0.1%甲酸和储液罐B中乙腈中0.1%甲酸组成。初始流动相由5%B组成,在2.5分钟内增加到90%B。在90%B下维持4分钟后,然后在100%B下冲洗柱1分钟,然后在初始条件下平衡1.5分钟,总运行时间为9分钟。在这些条件下,分析物和IS分别在3.52和2.86分钟洗脱。

使用SCIEX API 4000在正ESI模式下使用串联质谱法检测分析物和IS,并监测以下质量转变:ganetespib 365.2->131.1,SLV320 309.2->211.1。

实验生物分析运行均被发现可用于单一非GLP临床前PK评估。线性模型(1/X2称重)在5至200 ng/mL范围内拟合校准品,相关系数(R)≥0.99。在校准范围以上,以足够的精度和准确度稀释质控样品。定量下限(LLOQ)定义为峰面积信噪比为5或更大,与IS的基质空白相比,由于稀释因子,血浆为5 ng/mL,组织为30 ng/mL。整个矩阵的批内精密度和准确度分别≤9.62%CV和98.5%至111%。

根据矩阵和时间点对合成的ganetespib浓度-时间(Ct)数据进行分组,手动插补低于定量下限(BLOQ)的数据如下:如果在任何时间点≥2/3rds的Ct结果高于LLOQ,则BLOQ数据替换为½LLOQ值,ELSE将整个时间点的数据视为缺失。然后,使用Phoenix WinNonlin 6.4(Certara USA,Inc.,Princeton,NJ)生成Ct数据汇总统计数据(算术平均值、标准偏差、%CV、最小值、中位数、最大值),并对每个基质的ganetespib算术平均Ct数据进行非部门药代动力学分析(NCA)。采用静脉推注模型(模型201),并使用“线性向上-向下”梯形规则估计Ct曲线(AUC)值下的面积。终相被定义为Ct曲线末端的三个时间点,使用终相的未加权对数线性回归估计消除速率常数(Ke)。终端消除半衰期(T1/2)估计为0.693/Ke,从时间0到无穷大的AUC(AUCinf)估计为到最后一个时间点的AUC+预测的碎屑/Ke。

估计的其他NCA参数包括反萃初始浓度(C0)、观察到的最大浓度(Cmax)、Cmax时间(Tmax),最后观察到的时间点的浓度(Clast)、Clast时间(Tlast)、清除率(CL=剂量/AUCinf)和稳态分布体积(Vss)。ganetespib从血浆到相关组织的表观分配系数(Kp,组织)估计为可用时AUCinf、组织与AUCinf血浆的比值。为了估计临床相关的小鼠剂量,将得到的小鼠血浆AUCinf与在假定的单剂甘尼特司匹布(200 mg/m)下报告的人血浆PK值进行比较2这是建议的第2阶段剂量(RP2D)(Goldman等人,2013年). 所有推断都是在小鼠和人类中时间无关、线性和剂量比例PK的假设下进行的。

Ganetespib浓度在血浆中表现出中等的变异性,表明在整个采样时间点的变异系数为2.9%至65.6%。所有血浆和肿瘤结果均高于LLOQ。肿瘤穿透似乎很快但不太严重,基于AUCinf,Kp肿瘤值为0.639。ganetespib在血浆中的T1/2为2.03小时。肿瘤的表观半衰期比在血浆中观察到的半衰期长(8.42小时),表明ganetespab对原位肿瘤的亲和力很高。肿瘤的明显Vss也表明ganetespib具有较高的原位肿瘤保留率。值得注意的是,我们的血浆和肿瘤PK与Shimamura等人发表的小鼠静脉注射125 mg/kg的PK一致(Shimamura等人,2012年). NCA参数估计值和按矩阵分组的平均(SD)ganetespib-Ct图如下所示:

在单剂ganetespib的第1阶段研究中,每周静脉注射一次,为期3周,为期4周,200mg/m2是RP2D(Goldman等人,2013年). ganetespib的主要不良反应是腹泻和疲劳;然而,视网膜毒性一直备受关注。在人类研究的剂量范围内(7至259 mg/m),Cmax和AUC值似乎以剂量比例方式增加2). 200 mg/m时人体总血浆AUCinf的估计值2约为7000 hr-ng/mL。假设小鼠的剂量比例和线性PK过程,提供类似总血浆AUCinf的小鼠等效剂量(MED)将为12.8 mg/kg。ganetespib的血浆蛋白结合仅被报告为“高”,缺乏关于物种间差异程度的数据。因此,在小鼠和人类血浆蛋白结合相似的假设下,MED只能通过血浆总浓度得出。

值得注意的是,有数据表明,高亲和力HSP90抑制剂(如ganetespib)可能不会表现出剂量比例PK,尤其是在小鼠中(Yamazaki等人,2013年). 这被认为是由于靶点介导的药物处置,其中化合物对高度丰富和特异性靶点的亲和力驱动非线性组织分布。我们在小鼠体内静脉注射50 mg/kg ganetespib后产生了有限的血浆PK,这表明可能是这种情况(数据未显示)。在这种情况下,剂量从50毫克/千克增加到150毫克/千克的3倍会导致注射后10分钟血浆Cmax增加约17倍。因此,当我们从150 mg/kg剂量的AUC中反萃MED时,很难假设成比例。考虑到这一点,我们在小鼠药效模型中评估了ganetespib的两个剂量水平:10 mg/kg(可能接近MED)和50 mg/kg。

治疗RD细胞的蛋白质印迹

RD细胞系取自美国型培养物收集(ATCC,CCL-136)。

用0.1μM Ganetespib、0.25μM伊立替康和3 nM长春新碱或载体(DMSO)处理细胞,并在处理24或48小时后收获。细胞在含有蛋白酶抑制剂的RIPA裂解缓冲液(Sigma,目录R0278)中进行裂解(Halt Inhibitor Cocktail,Thermo Fisher,目录87785)。蛋白质在4–15%SDS-PAGE凝胶中溶解并转移到PVDF膜(Millipore)。在Odyssey阻断缓冲液(Licor,目录927-4000)中阻断膜,然后探测以下抗体:HSP70(1:1000,细胞信号,目录4872)和肌动蛋白(1:5000,Sigma,AC-15)。次要抗体为抗兔抗体(Licor,目录926–32221);防鼠(Licor,目录926–32210)。

流式细胞术

用1μM AZD1775、0.25μM伊立替康、3 nM长春新碱或载体(DMSO)处理RD细胞24或48小时。用碘化丙锭溶液(0.05 mg/mL碘化丙啶、0.1%柠檬酸钠、0.1%(v/v)Triton X-100、,2 ug/mL RNA酶),并在Becton-Dickson FACS系统(美国马里兰州Rockville)上通过流式细胞术分析细胞周期阶段。

QPCR(质量控制报告)

RD细胞用0.1μM-0.2μM Ganetespib、0.25μM-0.5μM伊立替康、3-6 nM长春新碱或载体(DMSO)处理,并在处理后24或48小时收获。使用Trizol提取RNA,并使用Superscript III系统(Invitrogen,目录18080093)合成cDNA。CDNA样品使用SYBR Select Master Mix for CFX(Applied Biosystems,目录4472942)进行分析,一式两份,并归一化为GAPDH(CCAGCAAGACAGGAA,GCCCCCCCTCTTCAAG)。对HSF-1靶基因进行QPCR(DARS:CCTGTAGTGCACAGGCTG,CGAACTTATACTGAATTACTCC;HSP70:CAGTTAGTCTTACTGAGCTT,GCTGTAGCTATAGCTCAACT)。对来自异种移植物组织的HSV-1靶基因的cDNA进行QPCR(ENY2:GAACCACTTCTGAATATAGCA,CTTAGGCAACAGTGTGGA;HNRNPA2B1:GGATGAGCCCAGAGGTA,TCTCTGCATGCTCTGTGGT;LTBP4:CATCCTGAGCGCTTAT,GTAGGTCCCTCTCTCTCCAGGT;USPL1:GAACTTACTGTGTTCAGGTa,GATTACGCTTCATGGGCA;ZNF467:CAGAGATCAGCGAACGTT,CACTCTTTCTCGCCCTGGTT)。

显微镜

用1μM AZD1775、0.25μM伊立替康、3 nM长春新碱或载体(DMSO)处理RD细胞24或48小时。用EdU处理细胞1小时,然后在4%PFA中培养以固定细胞。按照制造商的说明进行EdU标记(Click-iT EdU成像套件、Invitrogen、目录10083元)DNA用0.2μg/ml DAPI(Sigma-Aldrich)染色。细胞在蔡司Axioplan显微镜上进行评估。与治疗组相比,细胞核破碎的细胞以盲法评分为有丝分裂突变。多个随机字段重复计分。

临床前测试

如前所述,使用第一、二、三阶段试验设计进行临床前测试(Stewart等人,2014年).

临床前阶段I

使用非荷瘤裸鼠进行以下组合的耐受性试验:AZD1775+VCR+IRN和Ganetespib+VCR+IRN。给小鼠进行4个周期的化疗,每个周期21天。伊立替康在两个不同的时间表中进行测试(低剂量延长时间表以及与儿科患者相似的5天高剂量时间表)。所使用的化疗药物组合和时间表旨在模拟潜在的人类临床试验。每个治疗组包含3只小鼠。在第1天、第8天和第15天,每周腹腔注射一次长春新碱。伊立替康在第1-5天和第8-12天每天腹腔注射一次,或仅在第1-5.天给药。AZD1775在第1-5天每天口服两次。在第1天、第8天和第15天通过尾静脉每周给药一次Ganetespib。对每个治疗组的小鼠体重和全血细胞计数进行监测。在整个治疗过程中,每天监测动物的健康状况。所有治疗组的所有小鼠均存活了4个疗程。在所有测试组中均未发现明显的体重减轻。

临床前阶段II

采用上述肌肉注射技术,将来自SJRHB000026_X1(ERMS)、SJRHB12_Y(ERMS(Chen等人,2013). 每周用Xenogen筛选小鼠,并测量生物发光。小鼠在达到目标生物发光信号10后被纳入研究6−107光子/秒/厘米2或2周以上,或可触及肿瘤,并于下周一开始化疗。使用3种不同的异种移植物(SJRHB000026_X1、SJRHB012_Y、SJRHB013759_X1)进行以下临床前II期试验:

SJRHB000026_X1(ERMS):

将小鼠随机分为以下治疗组:

对照组、Ganetespib(GSP)、伊立替康(IRN)、GSP+IRN、AZD1775+IRN,长春新碱(VCR)+IRN和AZD11775+VCR+IRN。

每种药物的剂量和时间表如下:

普惠制50 mg/kg尾静脉,每周一次,第1、8、15天
AZD1775型60 mg/kg口服灌胃,每日两次,第1-5天
爱尔兰共和国1.25 mg/kg每天腹腔注射一次,第1-5天和第8-12天
录像机0.38 mg/kg每天腹腔注射一次,第1、8、15天
帕尔博西利布第1-21天每天口服10 mg/kg
阿贝马奇利布第1-21天每天口服15 mg/kg
曲美替尼第1-21天每天口服0.15 mg/kg
控制车辆控制(无化疗)

小鼠接受4个疗程的化疗(每个疗程3周),每周和治疗结束时监测生物发光。根据生物发光信号对疾病反应进行分类。信号为10的老鼠5光子/秒/厘米2或更少(与背景相似)被归类为完全应答,105-106光子/秒/厘米2作为部分响应,107最多10个8光子/秒/厘米2(与登记信号相似)为稳定疾病,且大于108光子/秒/厘米2作为进行性疾病。任何时候肿瘤负担超过体重20%的小鼠也被归类为进行性疾病。在接受化疗的同时,每天对小鼠进行监测。

SJRHB012_Y(ERMS):

将小鼠随机分为以下治疗组:

控制、GSP+IRN、VCR+IRN,AZD1775+VCR+IR,AZD2775+VC+IRN*,GSP+VCR+IRN,Palbociclib+Trametinib

每种药物的剂量和时间表如下:

普惠制每周第1、8、15天尾静脉注射50 mg/kg
AZD1775型60 mg/kg口服灌胃,每日两次,第1-5天
爱尔兰共和国1.25 mg/kg腹膜内注射,每天一次,第1-5天和第8-12天
爱尔兰共和国*3.125 mg/kg每天腹腔注射一次,第1-5天
录像机0.38 mg/kg腹膜内注射,每日一次,第1、8、15天
帕尔博西利布第1-21天每天口服10 mg/kg
曲美替尼第1-21天每天口服0.15 mg/kg
控制车辆控制(无化疗)

小鼠接受6个疗程的化疗,每周监测生物发光。根据生物发光信号对疾病反应进行分类。信号为10的老鼠5光子/秒/厘米2或更少(与背景相似)被归类为完全反应,105-106光子/秒/厘米2作为部分响应,107最多10个8光子/秒/厘米2(与登记信号相似)为稳定疾病,且大于108光子/秒/厘米2作为进行性疾病。任何时候肿瘤负担超过体重20%的小鼠也被归类为进行性疾病。在接受化疗的同时,每天对小鼠进行监测。

SJRHB013759_X1(臂):

将小鼠随机分为以下治疗组:

对照组、GSP+IRN、VCR+IRN,AZD1775+VCR+IR,GSP+VCR+IRN,Palbociclib+Trametinib

每种药物的剂量和时间表如下:

普惠制每周第1、8、15天尾静脉注射50 mg/kg
AZD1775型60 mg/kg口服灌胃,每日两次,第1-5天
爱尔兰共和国1.25 mg/kg每天腹腔注射一次,第1-5天和第8-12天
录像机0.38 mg/kg每天腹腔注射一次,第1、8、15天
帕尔博西利布第1-21天每天口服10 mg/kg
曲美替尼第1-21天每天口服0.15 mg/kg
控制车辆控制(无化疗)

小鼠接受6个疗程的化疗,并每周监测生物发光。根据生物发光信号对疾病反应进行分类。信号为10的小鼠5光子/秒/厘米2或更少(类似于背景)被归类为完全应答,105–106光子/秒/厘米2作为部分响应,107最多10个8光子/秒/厘米2(与登记信号相似)为稳定疾病,且大于108光子/秒/厘米2作为进行性疾病。任何时候肿瘤负担超过体重20%的小鼠也被归类为进行性疾病。在接受化疗的同时,每天对小鼠进行监测。

临床前III期

我们对来自4个不同O-PDX系的499只小鼠进行了随机、双盲安慰剂对照临床前III期试验。采用上述肌肉注射技术,将来自SJRHB000026_X1(ERMS)、SJRHB012_Y(ERMS,SJRHB13757_X2(ARMS)和SJRHB 013759_X1的荧光素酶标记细胞注射到受体裸鼠体内,建立RMS原位异种移植物。每周用Xenogen筛选小鼠,并测量生物发光。在目标生物发光信号达到10后的11周内,小鼠被纳入研究6−107光子/秒/厘米2或2周以上,或可触及肿瘤,并于下周一开始化疗。接受化疗的患者和进行生物发光成像的患者对所服用的药物以及特定的O-PDX肿瘤细胞系一无所知。

同时对以下14个治疗组的所有四种O-PDX进行了测试:在临床前III期试验中使用了以下药物和剂量,该试验由AUC引导给药,药物动力学研究获得的剂量与临床相关的人体剂量密切匹配:

控制车辆/未处理
AZD1775型AZD1775在第1-5天每天两次以60 mg/kg剂量给药。
普惠制在第1天、第8天和第15天,GSP每周剂量为50 mg/kg。
GSP+IRN在第1天、第8天、第15天,GSP每周剂量为50 mg/kg,IRN剂量为1.25
第1-5天和第8-12天每天mg/kg
GSP+IRNb条在第1、8、15天,GSP每周剂量为10 mg/kg,IRN剂量为1.25
第1-5天和第8-12天每天mg/kg
GSP+VCR+IRN在第1、8、15天,GSP每周剂量为50 mg/kg,VCR剂量为0.38
第1天、第8天、第15天每周服用mg/kg,每日服用1.25 mg/kg IRN
在第1-5天和第8-12天
GSP+VCR+IRNb条在第1、8、15天,GSP每周剂量为10 mg/kg,VCR剂量为0.38
第1天、第8天、第15天每周服用mg/kg,每日服用1.25 mg/kg IRN
第1-5天和第8-12天
AZD1775+免疫核糖核酸AZD1775在第1-5天每天两次以48 mg/kg的剂量给药,IRN给药
第1天至第5天每天1.875 mg/kg(剂量根据建议
ADVL1312临床试验第2阶段剂量)
AZD1775+录像机+IRNAZD1775在第1-5天每天两次以48 mg/kg的剂量给药,VCR的剂量为
第1、8、15天每周0.38 mg/kg,IRN剂量为1.875 mg/kg
每天第1-5天
AZD1775+录像机+IRNb条AZD1775在第1-5天每天两次以48 mg/kg的剂量给药,VCR的剂量为
第1、8、15天每周0.38 mg/kg,IRN剂量为1.25 mg/kg
每天第1-5天和第8-12天
AZD1775+录像机+IRNcAZD1775在第1-5天每天两次以60 mg/kg的剂量给药,VCR的剂量为
第1、8、15天每周0.38 mg/kg,IRN剂量为1.25 mg/kg
每天第1-5天和第8-12天
VCR+IRN在第1天、第8天、第15天,VCR每周剂量为0.38 mg/kg,IRN每周剂量为
第1-5天和第8-12天每天1.25 mg/kg(相当于低剂量
儿童IRN剂量为20 mg/m2的长期计划)
录像机+IRNb条在第1天、第8天、第15天,VCR每周剂量为0.38 mg/kg,IRN每周剂量为
第1-5天每天1.875 mg/kg(相当于5天IRN剂量30
儿童mg/m2)。
录像机+IRNc在第1天、第8天、第15天,VCR每周剂量为0.38 mg/kg,IRN每周剂量为
第1-5天每天3.125 mg/kg(相当于5天IRN剂量50
儿童mg/m2)

小鼠接受6个疗程的化疗(每个疗程3周),并在每个疗程的第3周和治疗结束时监测生物发光。根据上述临床前阶段II的生物发光信号对疾病反应进行分类。如前所述,对所有患有进展性疾病的小鼠进行了研究结束时的检查(Stewart等人,2014年).

临床前测试用化疗药物

下表显示了用于混合体内各化疗药物的载体:

化疗车辆
AZD1775型0.5%甲基纤维素
爱尔兰共和国生理盐水
录像机生理盐水
甘尼特斯匹布10%二甲基亚砜/18%科尔利弗RH 40/72%D5W
帕尔博西利布1%甲基纤维素(400 cPs型)和1%吐温80
阿贝马奇利布1%羟乙基纤维素(HEC)、0.25%吐温80和
超纯水中约0.05%硅氧烷
曲美替尼1%甲基纤维素(400 cPs型)和1%吐温80

氙气成像和量化

小鼠腹腔注射萤火虫D-荧光素(Caliper Life Sciences 3 mg/小鼠)。五分钟后,使用IVIS®200成像系统拍摄生物发光图像。在整个图像采集过程中进行麻醉(O中异氟醚1.5%2以2 L/min的速度交付)。使用Living Image 4.3软件(Caliper Life Sciences)生成包含最大生物发光信号下最大肿瘤的标准感兴趣区域(ROI)。使用相同的ROI来确定平均辐射(光子/s/cm2/sr)用于所有异种移植物。

量化和统计分析

ChIP测序分析

我们首先使用BWA(版本0.5.9-r26-dev,默认参数)将ChIP-Seq读数与人类基因组hg19(GRCh37-lite)对齐。然后使用Picard(版本1.65(1160))标记重复读取。然后samtools只保留非重复读取(参数“-q 1 F 1024”版本0.1.18(r982:295))。我们遵循ENCODE标准进行质量控制(QC),然后使用SPP的非重复版本(版本1.11)的数据绘制互相关,并在R(版本2.14.0)支持下使用软件包caTools(版本1.17)和bitops(版本1.0-6)计算相对链相关值(RSC),并估计片段大小。我们需要对点源因子(H3K4me2/3、H3K9/14Ac、H3K27Ac、CTCF、RNAPolII、BRD4)进行>10M的唯一映射读取,RSC>1。我们需要20 M宽标记(H3K9me3、H3K27me3和H3K36me3)的唯一映射读取。对于INPUT和RSC<1,我们需要10 M唯一映射读取。我们注意到H3K4me1在某些阶段是一个点源因子,而在其他阶段是广泛的,所以我们将H3K4me1作为广泛的标记物进行质控。然后,在手动检查SPP生成的互相关图后,使用估计的最佳片段大小(SPP在所有情况下估计的最小片段大小)来扩展每次读取并生成一个bigwig文件以在IGV上查看(版本2.3.40)。手动检查所有剖面是否有清晰的峰值和良好的信噪分离。

隐马尔可夫模型分析

非重复对齐读取通过上文定义的片段大小进行扩展,ChromHMM(1.10版,BinarizeBed使用“-colfields 0,1,2,5-center”)用于染色质状态建模。为了选择状态编号,我们首先对所有鼠标开发阶段(从7个状态到33个状态)进行建模,并在手动检查后选择16个状态的模型。为了更好地显示HMM状态在各个阶段的动态,我们通过基因和侧翼区域所覆盖状态的最大总百分比对颜色强度进行了标准化。为了确定代表基因的最佳区域,我们首先通过TSS在任何发育阶段的任何H3K4me3/H3K4me2/H3K27Ac/H3K9–14Ac峰值的2 kb范围内过滤注释亚型,然后我们选择在任何发育时期由cuffdiff估计的最高表达亚型,或者如果没有由cuffdiff估计的表达水平,则选择最长的亚型。最后,我们将HMM状态的间隔减少到半巴,如果它不在基因的前2个HMM状态中,则强度减少到标准化强度的一半。由于HMM状态可以由多个基因分配,因此跨基因的最大总百分比已用于这些标准化。

核心调节电路(CRC)分析

CRC映射器按前面所述运行(Aldiri等人,2017)除了我们进行了两次单独的分析来比较有启动子区和没有启动子区的数据。我们使用我们的H3K4me3 ChIP-seq数据来删除那些与启动子重叠的假定超增强子。这是因为一些非常活跃表达的基因具有被称为超增强子的H3K27Ac的大区域,并且很难区分这些假定的超增强器和基因上的强启动子。

DNA甲基化和分析

全基因组亚硫酸氢盐测序(WGBS)数据映射到GRCh37-lite(人类)或BSMAP2.74,参数如下:-z 33-f 5-g 3-r 0-m 17-x 600–u,原始甲基化状态由Bis-SNP调用(版本0.82.2)。我们开发了一种基于回归树的方法,使用m值识别样本的区域甲基化模式(=日志2#M(M)e(电子)t吨小时t吨e(电子)d日C类+0.5#U型n个e(电子)t吨小时t吨e(电子)d日C类+0.5). 相邻节段的平均β值差异不显著(p>1E-6,学生t检验)或很小(β=#M(M)e(电子)t吨小时t吨e(电子)d日C类#M(M)e(电子)t吨小时t吨e(电子)d日C类+#U型n个e(电子)t吨小时t吨e(电子)d日C类)(<0.1)递归合并。

收集单个样本的分段断点,生成一个区域的初始列表,用于测试甲基化和基因表达之间的相关性。计算可变甲基化区域(最高和最低平均m值之间的差异>=0.3)与可变表达相邻基因(甲基化区域重叠TSS到TES上游10 kb,最高和最低FPKM之间的差异至少为2倍,最高FPKM>=1)之间的相关性。保留FDR<=0.5的基因/区域对进行合并,合并具有相同方向(正相关或负相关)的相同基因的相邻区域(相距不超过1 kb)。重新评估相关性,并将最终FDR<=0.1的配对用作与基因表达显著相关的差异甲基化区域。

综合分析

为了优先考虑肿瘤中可能失调的基因,我们基于顺序统计整合了多种类型的数据源(Zhang等人,2012年). 具体来说,此方法将作为输入N个不同的基因排序集(详细信息如下所述),并通过组合所有N个数据源,基于公式:

(第页1,第页2,...,第页N个)=N个0第页11第页2...N个1第页N个d日N个d日N个1...d日1
(1)

哪里第页是数据源的秩比,N个是数据源的总数,并且第页0= 0. 实际上使用更有效的算法计算每个基因的值(Aerts等人,2006年):

(第页1,第页2,...第页N个)=N个V(V)N个
(2)

V(V)k个==1k个(1)1V(V)k个第页N个k个+1
(3)

统计数据不是均匀分布的(Stuart等人,2003年),我们按照公布的程序进行排列,以估计每个基因的经验p值(Zhang等人,2012年). 在每个排列中,缺失值被随机分配给k个基因(k由实际数据确定),并且每个数据源的基因等级被排列,并且根据排列数据矩阵计算值。使用至少具有3种数据类型(n=~12k)的所有基因构建零分布。经验p值由估计的零分布导出,然后使用BH方法进行多次测试校正(本杰米尼和霍奇伯格,1995年).

考虑四种类型的数据源(即全蛋白质组、磷酸蛋白质组、转录组和表观基因组)进行整合。对于基于RNA(来自转录组)和蛋白质(来自整个蛋白质组)丰度的基因排名,基因是根据每个肿瘤样本相对于对照(即成肌细胞样本)的对数2倍变化进行排名的。请注意,基因可以按降序排列(即顶部上调/激活程度最高的基因)或按升序排列(例如顶部下调/抑制程度最深的基因),这两种情况都在我们的分析中得到了考虑。对于表观基因组,我们计算了每个样本中每个基因(定义为TSS到TES加上5'端和3'端的2kb)的活性状态(1、2、3、4、5、6和10)和抑制/非活性状态(14、17和18)之间的分数差异,并根据每个肿瘤与对照的派别差异对基因进行排序。对于磷酸化蛋白质组,每种蛋白质的磷酸化水平通过所有相关磷酸肽的相对强度比的基于光谱计数的加权平均值来总结(Bai等人,未发表的数据),并通过log2比率对基因进行排序。在为每个数据源定义了基因排名后,通过上述顺序统计方法为每个肿瘤样本生成一个综合排名。为了总结每个肿瘤亚型(A/E-RMS)的基因排名,对于每个数据源,首先对每个亚型样本的log2比值进行平均,并对多个数据源的基因排名进行类似的整合。为了生成肿瘤亚型特异性基因排序,我们通过比较两种亚型的平均值(即ERMS与ARMS)来计算每个数据源的log2比率,然后进行整合。

最后,为了将基因水平的综合排名总结为通路水平,我们进行了基因集富集分析(Subramanian等人,2005年)在javaGSEA应用程序中使用预先排序的函数(http://software.broadinstitute.org/gsea/index.jsp). 简言之,基因通过综合蛋白质排序(针对每个肿瘤样本、亚型或亚型特异性)进行排序,FDR值通过排列路径数据库中的基因集得出。FDR<0.05的路径被认为是显著的。

临床前III期研究的统计分析

临床前III期研究的研究前咨询由圣犹德儿童研究医院生物统计部完成。数据包括来自4个O-PDX模型的14个不同治疗组共499只小鼠。该研究的目的是评估和比较各治疗组的肿瘤反应和无肿瘤进展生存率。按照生物统计学部提供的随机化代码,将小鼠随机分为每个治疗组。如上所述,肿瘤反应通过治疗结束时的生物发光来定义。如果一只小鼠在登记后的任何时间因肿瘤大小超过小鼠体重的20%而退出研究,则会自动分配给它们进行性疾病反应。采用Kaplan-Meier方法绘制肿瘤进展时间的生存曲线。采用对数秩检验比较各组的生存曲线。

数据和软件可用性

基因组、转录组、表观遗传和蛋白质组数据文件可在公共在线门户网站上获得(https://pecan.stjude.org/proteinpaint/study/RHB2018). 所有序列数据保存在欧洲分子生物学实验室核苷酸序列数据库中(登录号:EGAS00001002967)。有关软件可用性的详细信息,请参阅关键资源表。

集锦

  • 综合分析可深入了解RMS的发展起源
  • 与ERMS相比,ARMS在发展计划中发生得更远
  • RAS/MEK/ERK/CDK4/6,G2/RMS中M和UPR通路被解除管制
  • 靶向G中的WEE1激酶2/M通路在高危RMS中最有效

通过综合分析,Stewart等人揭示了横纹肌肉瘤(RMS)的发育起源,并确定了常见的非调控途径。全面的临床前测试有助于为未来的临床试验优先考虑药物,并表明WEE1是体内高风险RMS的高效治疗靶点。

重要性

本研究证明了整合转录组学、表观基因组学和蛋白质组/磷酸蛋白质组学数据的价值,以进一步完善发育性肿瘤的细胞起源,并确定超出基因组体细胞突变的肿瘤脆弱性。我们的结果证明了使用多个原位患者异种移植物进行全面临床前测试的重要性,以验证并优先考虑通过基因组或综合分析确定的脆弱性。根据这些数据,高危复发性RMS患者可能受益于包括AZD1775、伊立替康和长春新碱在内的联合化疗。为了促进这些数据的传播和未来的研究,我们开发了一个在线门户(https://pecan.stjude.org/proteinpaint/study/RHB2018)供公众使用。

补充材料

1

单击此处查看。(13K,xlsx)

9

单击此处查看。(114K,xlsx)

10

单击此处查看。(56M,pdf格式)

2

单击此处查看。(130万,xlsx)

单击此处查看。(88K,xlsx)

4

单击此处查看。(10万,xlsx)

5

单击此处查看。(26K,xlsx)

6

单击此处查看。(13M,xlsx)

7

单击此处查看。(94K,xlsx)

8

单击此处查看。(228K,xlsx)

致谢

我们感谢尼莎·巴德斯编辑手稿。这项工作部分得到了NCI癌症中心支持(CA21765)、NIH M.A.D(EY014867和EY018599和CA168875)和ALSAC的资助。医学博士还得到了亚历克斯柠檬水儿童癌症基金会、塔利家庭基金会和彼得森基金会的支持。E.S.得到了国家综合癌症网络(National Comprehensive Cancer Network)、国家儿科癌症基金会(National Pediatric Cancer Foundation)的支持,是圣·鲍德里克学院(St.Baldrick’S Scholar)的学者,得到了Invictus Fund的慷慨支持。J.P.由NIH资助(AG047928)。这项研究的大部分得到了霍华德·休斯医学研究所的支持。WGS、WGBS和RNAseq是圣犹德儿童研究医院/华盛顿大学儿科癌症基因组项目的一部分。临床前研究在动物影像共享资源的协助下进行;在细胞和组织成像共享资源的帮助下进行电子显微镜检查;药代动力学在临床前药代动力学共享资源的帮助下进行;MS在蛋白质组学共享资源中进行,组织病理学分析在St.Jude兽医病理学共享资源的协助下进行。

脚注

权益声明

作者声明没有相互竞争的利益。

出版商免责声明:这是一份未经编辑的手稿的PDF文件,已被接受出版。作为对客户的服务,我们正在提供这份早期版本的手稿。在以最终可引用的形式出版之前,手稿将经过文案编辑、排版和校对。请注意,在制作过程中可能会发现可能影响内容的错误,适用于期刊的所有法律免责声明都适用。

参考文献

  • Aerts S、Lambrechts D、Maity S、Van Loo P、Coessens B、De Smet F、Tranchevent LC、De Moor B、Marynen P、Hassan B等人(2006年)。通过基因组数据融合实现基因优先排序.Nat生物技术 24, 537–544. [公共医学][谷歌学者]
  • Afzali AM、Ruck T、Herrmann AM、Iking J、Sommer C、Kleinschnitz C、Preubetae C、Stenzel W、Budde T、Wiendl H等人(2016年)。钾通道TASK2和TREK1调节小鼠骨骼肌细胞的功能分化.美国生理学杂志细胞生理学 311,C583–C595。[公共医学][谷歌学者]
  • Akerfelt M、Morimoto RI和Sistonen L(2010年)。热休克因素:细胞应激、发育和寿命的集合体.《自然》杂志综述分子细胞生物学 11, 545–555.[PMC免费文章][公共医学][谷歌学者]
  • Alagesan B、Contino G、Guimaraes AR、Corcoran RB、Deshpande V、Wojtkiewicz GR、Hezel AF、Wong KK、Loda M、Weissleder R等(2015)。MEK和PI3K联合抑制小鼠胰腺癌模型.临床癌症研究 21, 396–404.[PMC免费文章][公共医学][谷歌学者]
  • Aldiri I、Xu B、Wang L、Chen X、Hiler D、Griffiths L、Valentine M、Shirinifard A、Thiagarajan S、Sablauer A等(2017)。视网膜在发育、重编程和肿瘤发生过程中的动态表观遗传景观.神经元 94,550-568 e510。[PMC免费文章][公共医学][谷歌学者]
  • Allen CE、Laetsch TW、Mody R、Irwin MS、Lim MS、Adamson PC、Seibel NL、Parsons DW、Cho YJ、Janeway K等人(2017年)。儿童肿瘤组的靶点和药物优先顺序-国家癌症研究所儿科MATCH试验.美国国家癌症研究所 109.[PMC免费文章][公共医学][谷歌学者]
  • Bai B、Tan H、Pagala VR、High AA、Ichaporia VP、Hendershot L和Peng J(2017)。通过等压标记、广泛液相色谱和质谱法对蛋白质组和磷蛋白组进行深度剖析.酶学方法 585,377–395页。[PMC免费文章][公共医学][谷歌学者]
  • Benjamini Y和Hochberg Y(1995)。控制错误发现率:一种实用而有效的多重测试方法. .J R统计Soc B 57, 289–300.[谷歌学者]
  • Bennett AM和Tonks NK(1997年)。丝裂原活化蛋白激酶对骨骼肌分化不同阶段的调控.科学类 278, 1288–1291. [公共医学][谷歌学者]
  • Borycki AG、Brunk B、Tajbakhsh S、Buckingham M、Chiang C和Emerson CP Jr.(1999)。声波刺猬通过Myf5激活控制轴旁肌的测定.发展 126, 4053–4063. [公共医学][谷歌学者]
  • Carthy JM、Sundqvist A、Heldin A、van Dam H、Kletsas D、Heldin-CH和Moustakas A(2015)。三苯氧胺通过ERK1/2阻断非MAD信号通路抑制TGF-β介导的肌成纤维细胞活化.J细胞生理学 230, 3084–3092. [公共医学][谷歌学者]
  • Chen AE、Ginty DD和Fan CM(2005年)。通过CREB的蛋白激酶A信号传导控制Wnt蛋白诱导的肌生成.自然 433, 317–322. [公共医学][谷歌学者]
  • Chen X、Stewart E、Shelat AA、Qu C、Bahrami A、Hatley M、Wu G、Bradley C、McEvoy J、Pappo A等(2013)。针对胚胎性横纹肌肉瘤的氧化应激.癌细胞 24, 710–724.[PMC免费文章][公共医学][谷歌学者]
  • Cossu G和Borello U(1999年)。Wnt信号转导与哺乳动物肌肉发生的激活.Embo J公司 18, 6867–6872.[PMC免费文章][公共医学][谷歌学者]
  • Do K、Doroshow JH和Kummar S(2013)。Wee1激酶作为癌症治疗的靶点.细胞周期 12, 3159–3164.[PMC免费文章][公共医学][谷歌学者]
  • Ernst J和Kellis M(2012年)。ChromHMM:自动化染色质状态发现和表征.自然方法 9, 215–216.[PMC免费文章][公共医学][谷歌学者]
  • Furman WL、Stewart CF、Poquette CA、Pratt CB、Santana VM、Zamboni WC、Bowman LC、Ma MK、Hoffer FA、Meyer WH等(1999年)。将长期伊立替康时间表从异种移植模型直接转换为儿童I期试验.临床肿瘤学杂志:美国临床肿瘤学会官方杂志 17, 1815–1824. [公共医学][谷歌学者]
  • Goldman JW、Raju RN、Gordon GA、El-Hariry I、Teofilivic F、Vukovic VM、Bradley R、Karol MD、Chen Y、Guo W等(2013)。固体恶性肿瘤患者每周服用一次Hsp90抑制剂ganetespib(STA-9090)的人体、安全性、药代动力学和临床活性第一阶段研究.BMC癌症 13, 152.[PMC免费文章][公共医学][谷歌学者]
  • 龚J、袁Y、病房A、康L、张B、吴Z、彭J、冯Z、刘J和徐XZS(2017)。秀丽隐杆线虫味觉受体同源物LITE-1是一种光受体.单元格 168, 325. [公共医学][谷歌学者]
  • Harris MH、DuBois SG、Glade Bender JL、Kim A、Crompton BD、Parker E、Dumont IP、Hong AL、Guo D、Church A等人(2016)。肿瘤分子谱分析为晚期儿童实体瘤的治疗决策提供信息的多中心可行性研究:个体化癌症治疗(iCat)研究.JAMA Oncol公司. [公共医学][谷歌学者]
  • Hostein I、Andraud-Fregeville M、Guillou L、Terrier-Lacombe MJ、Deminiere C、Ranchere D、Lussan C、Longavenne E、Bui NB、Delatter O和Coindre JM(2004)。横纹肌肉瘤:肌生成素表达分析和分子检测在肺泡亚型诊断中的价值:109例石蜡包埋标本的分析.癌症 101, 2817–2824. [公共医学][谷歌学者]
  • Hovestadt V、Jones DT、Picelli S、Wang W、Kool M、Northcott PA、Sultan M、Stachurski K、Ryzhova M、Warnatz HJ等人(2014年)。利用DNA甲基化测序解码髓母细胞瘤的调控景观.自然 510, 537–541. [公共医学][谷歌学者]
  • 姜博华,郑JZ,Vogt PK(1998)。磷脂酰肌醇3激酶在肌源性分化中的重要作用.美国国家科学院程序 95, 14179–14183.[PMC免费文章][公共医学][谷歌学者]
  • Jokinen E和Koivunen JP(2015)。实体瘤中MEK和PI3K抑制:理论基础和证据.Ther Adv Med Oncol公司 7, 170–180.[PMC免费文章][公共医学][谷歌学者]
  • Karagoz GE和Rudiger SG(2015年)。Hsp90与客户互动.生物化学趋势 40, 117–125. [公共医学][谷歌学者]
  • Langfelder P和Horvath S(2008年)。WGCNA:加权相关网络分析的R包.BMC生物信息学 9, 559.[PMC免费文章][公共医学][谷歌学者]
  • Langfelder P、Zhang B和Horvath S(2008)。从层次化集群树定义集群:R的Dynamic tree Cut包.生物信息学 24, 719–720. [公共医学][谷歌学者]
  • Liang X、Ubhayakar S、Liederer BM、Dean B、Ran-Ran Qin A、Shahidi-Latham S和Deng Y(2011)。用于制备小鼠组织样品以支持药物发现的均质技术评估.生物分析 , 1923–1933. [公共医学][谷歌学者]
  • Liberzon A、Subramanian A、Pinchback R、Thorvaldsdottir H、Tamayo P和Mesirov JP(2011年)。分子签名数据库(MSigDB)3.0.生物信息学 27, 1739–1740.[PMC免费文章][公共医学][谷歌学者]
  • Lister R、Pelizzola M、Dowen RH、Hawkins RD、Hon G、Tonti-Filippini J、Nery JR、Lee L、Ye Z、Ngo QM等(2009年)。碱基分辨率的人类DNA甲基体显示出广泛的表观基因组差异.自然 462, 315–322.[PMC免费文章][公共医学][谷歌学者]
  • McAlister GC、Nusinow DP、Jedrychowski MP、Wur M、Huttlin EL、Erickson BK、Rad R、Haas W和Gygi SP(2014年)。MultiNotch MS3能够准确、灵敏和多路检测癌细胞系蛋白质组的差异表达.分析化学 86, 7150–7158.[PMC免费文章][公共医学][谷歌学者]
  • Nesvizhskii AI和Aeberseld R(2005年)。霰弹枪蛋白质组学数据的解释:蛋白质推断问题。分子和细胞蛋白质组学:MCP公司 4, 1419–1440. [公共医学][谷歌学者]
  • Newton WA Jr.、Soule EH、Hamoudi AB、Reiman HM、Shimada H、Beltangady M和Maurer H(1988)。儿童肉瘤的组织病理学,组间横纹肌肉瘤研究I和II:临床病理相关性.临床肿瘤学杂志:美国临床肿瘤学会官方杂志 6, 67–75. [公共医学][谷歌学者]
  • Niu M、Cho JH、Kodali K、Pagala V、High AA、Wang H、Wu Z、Li Y、Bi W、Zhang H等(2017)。广泛肽分馏和基于y1离子的干扰检测方法,通过等压标记和质谱实现准确定量.分析化学 89, 2956–2963.[PMC免费文章][公共医学][谷歌学者]
  • O'Geen H、Squazzo SL、Iyengar S、Blahnik K、Rinn JL、Chang HY、Green R和Farnham PJ(2007年)。KAP1结合的全基因组分析提示KRAB-ZNFs的自身调节.公共科学图书馆-基因 ,e89。[PMC免费文章][公共医学][谷歌学者]
  • Ognjanovic S、Linabery AM、Charbonneau B和Ross JA(2009年)。1975年至2005年美国儿童横纹肌肉瘤发病率和生存率趋势.癌症 115, 4218–4226.[PMC免费文章][公共医学][谷歌学者]
  • Ouellet D、Kassir N、Chiu J、Mouksassi MS、Leonowens C、Cox D、DeMarini DJ、Gardner O、Crist W和Patel K(2016)。MEK抑制剂曲美替尼在BRAF V600突变阳性黑色素瘤患者中的人群药代动力学和暴露反应.癌症化疗与药理学 77, 807–817. [公共医学][谷歌学者]
  • Pappo AS、Anderson JR、Crist WM、Wharam MD、Breitfeld PP、Hawkins D、Raney RB、Womer RB、Parham DM、Qualman SJ和Grier HE(1999)。儿童和青少年横纹肌肉瘤复发后的生存率:来自组间横纹肌肉瘤研究小组的报告.临床肿瘤学杂志:美国临床肿瘤学会官方杂志 17, 3487–3493. [公共医学][谷歌学者]
  • Parrish KE、Cen L、Murray J、Calligaris D、Kizilbash S、Mittapalli RK、Carlson BL、Schroeder MA、Sludden J、Boddy AV等(2015)。PARP抑制剂Rucaparib在原位胶质母细胞瘤异种移植中的疗效受限于药物对中枢神经系统的无效渗透.摩尔癌症治疗 14, 2735–2743.[PMC免费文章][公共医学][谷歌学者]
  • Parsons DW、Roy A、Yang Y、Wang T、Scollon S、Bergstrom K、Kerstein RA、Gutierrez S、Petersen AK、Bavle A等(2016)。儿童实体瘤临床肿瘤诊断率和种系全外显子序列测定.JAMA Oncol公司.[PMC免费文章][公共医学][谷歌学者]
  • Patnaik A、Rosen LS、Tolaney SM、Tolcher AW、Goldman JW、Gandhi L、Papadopoulos KP、Beeram M、Rasco DW、Hilton JF等(2016年)。CDK4和CDK6抑制剂Abemacilib治疗乳腺癌、非小细胞肺癌和其他实体瘤患者的疗效和安全性.癌症发现 6, 740–753. [公共医学][谷歌学者]
  • Raney RB、Anderson JR、Barr FG、Donaldson SS、Pappo AS、Qualman SJ、Wiener ES、Maurer HM和Crist WM(2001)。前二十年的横纹肌肉瘤和未分化肉瘤:组间横纹肌肉瘤研究组经验和组间横纹肌肉瘤研究V的理论基础的选择性回顾.儿科血液学/肿瘤学杂志 23,215–220。[公共医学][谷歌学者]
  • Raub TJ、Wishart GN、Kulandhaivel P、Staton BA、Ajamie RT、Sawada GA、Gelbert LM、Shannon HE、Sanchez Martinez C和De Dios A(2015)。颅内胶质母细胞瘤异种移植中两种选择性双重CDK4和CDK6抑制剂的脑暴露及联合替莫唑胺抑制CDK4、CDK6的抗肿瘤活性.药物代谢和处置:化学品的生物命运 43, 1360–1371. [公共医学][谷歌学者]
  • Renshaw J、Taylor KR、Bishop R、Valenti M、De Haven Brandon A、Gowan S、Eccles SA、Ruddle RR、Johnson LD、Raynaud FI等人(2013年)。双重阻断PI3K/AKT/mTOR(AZD8055)和RAS/MEK/ERK(AZD 6244)通路在体内外协同抑制横纹肌肉瘤细胞生长.临床癌症研究 19, 5940–5951.[PMC免费文章][公共医学][谷歌学者]
  • Reshef R、Maroto M和Lassar AB(1998年)。背躯体细胞命运的调节:BMPs和Noggin控制肌源性调节因子表达的时间和模式.基因开发 12, 290303.[PMC免费文章][公共医学][谷歌学者]
  • Ritz C和Streibig JC(2005年)。使用R进行生物分析.J Statist软件 12.[谷歌学者]
  • Sabnis AJ、Guerriro CJ、Olivas V、Sayana A、Shue J、Flanagan J、Asthana S、Paton AW、Paton JC、Gestwicki JE等(2016)。联合化学遗传学方法鉴定横纹肌肉瘤细胞溶质HSP70依赖性.美国国家科学院程序 113, 9015–9020.[PMC免费文章][公共医学][谷歌学者]
  • Saint-Andre V、Federation AJ、Lin CY、Abraham BJ、Reddy J、Lee TI、Bradner JE和Young RA(2016年)。人类核心转录调控回路模型.基因组研究 26, 385–396.[PMC免费文章][公共医学][谷歌学者]
  • Schneider G、Schmidt-Supprian M、Rad R和Saur D(2017年)。组织特异性肿瘤发生:背景因素.Nat Rev癌症 17, 239–253.[PMC免费文章][公共医学][谷歌学者]
  • Shern JF、Chen L、Chmielecki J、Wei JS、Patidar R、Rosenberg M、Ambrogio L、Auclair D、Wang J、Song YK等(2014)。横纹肌肉瘤的综合基因组分析揭示了影响融合阳性和融合阴性肿瘤共同遗传轴的改变.癌症发现 4, 216–231.[PMC免费文章][公共医学][谷歌学者]
  • Shimamura T、Perera SA、Foley KP、Sang J、Rodig SJ、Inoue T、Chen L、Li D、Carretro J、Li YC等(2012)。Ganetespib(STA-9090)是一种非长效霉素HSP90抑制剂,在非小细胞肺癌的体内外模型中具有强大的抗肿瘤活性.临床癌症研究 18, 4973–4985.[PMC免费文章][公共医学][谷歌学者]
  • Smith MA、Altekruse SF、Adamson PC、Reaman GH和Seibel NL(2014年)。儿童和青少年癌症死亡率下降.癌症 120, 2497–2506.[PMC免费文章][公共医学][谷歌学者]
  • Stewart E、Federico SM、Chen X、Shelat AA、Bradley C、Gordon B、Karlstrom A、Twarog NR、Clay MR、Bahrami A等(2017年)。儿童实体瘤的原位患者来源异种移植物.自然 549, 96–100.[PMC免费文章][公共医学][谷歌学者]
  • Stewart E、Goshorn R、Bradley C、Griffiths LM、Benavente C、Twarog NR、Miller GM、Caufield W、Freeman BB 3rd、Bahrami A等人(2014年)。靶向尤因肉瘤的DNA修复途径.单元格代表 9,829–841页。[PMC免费文章][公共医学][谷歌学者]
  • Stuart JM、Segal E、Koller D和Kim SK(2003年)。全球发现保守遗传模块的基因表达网络.科学类 302, 249–255. [公共医学][谷歌学者]
  • Subramanian A、Tamayo P、Mootha VK、Mukherjee S、Ebert BL、Gillette MA、Paulovich A、Pomeroy SL、Golub TR、Lander ES和Mesirov JP(2005)。基因集富集分析:基于知识的全基因组表达谱解释方法.美国国家科学院程序 102, 15545–15550.[PMC免费文章][公共医学][谷歌学者]
  • Tan H、Wu Z、Wang H、Bai B、Li Y、Wang X、Zhai B、Beach TG和Peng J(2015)。磷酸盐添加剂富集精制磷酸肽及人脑磷酸蛋白质组分析.蛋白质组学 15, 500–507.[PMC免费文章][公共医学][谷歌学者]
  • Tan H、Yang K、Li Y、Shaw TI、Wang Y、Blanco DB、Wang X、Cho JH、Wang H、Rankin S等(2017)。综合蛋白质组学和磷蛋白组学分析揭示了T细胞激活的动态信号网络和生物能量学途径.免疫 46, 488–503.[PMC免费文章][公共医学][谷歌学者]
  • Tate SC、Cai S、Ajamie RT、Burke T、Beckmann RP、Chan EM、De Dios A、Wishart GN、Gelbert LM和Cronier DM(2014)。新型细胞周期素依赖性激酶4/6抑制剂LY2835219在荷人肿瘤异种移植小鼠中抗肿瘤活性的半机械药代动力学/药效学模型.临床癌症研究 20, 3763–3774. [公共医学][谷歌学者]
  • Taus T、Kocher T、Pichler P、Paschke C、Schmidt A、Henrich C和Mechtler K(2011年)。使用phosphorRS进行普遍且可靠的磷酸化位点定位.蛋白质组研究杂志 10, 5354–5362. [公共医学][谷歌学者]
  • Twarog NR、Stewart E、Hammill CV和Shelat AA(2016年)。BRAID:联合药物作用分析的统一范式.科学报告.[PMC免费文章][公共医学][谷歌学者]
  • Vora SR、Juric D、Kim N、Mino-Kenudson M、Huynh T、Costa C、Lockerman EL、Pollack SF、Liu M、Li X等人(2014)。CDK 4/6抑制剂使PIK3CA突变型乳腺癌对PI3K抑制剂增敏.癌细胞 26, 136–149.[PMC免费文章][公共医学][谷歌学者]
  • 王X、李毅、吴Z、王赫、谭赫和彭杰(2014)。JUMP:一种基于标签的数据库搜索工具,用于高灵敏度和准确性的肽鉴定。分子和细胞蛋白质组学:MCP公司 13, 3663–3673.[PMC免费文章][公共医学][谷歌学者]
  • Williamson D、Missiaglia E、de Reynies A、Pierron G、Thuille B、Palenzuela G、Thway K、Orbach D、Lae M、Freneaux P等人(2010年)。融合基因阴性的肺泡型横纹肌肉瘤在临床和分子上与胚胎型横纹肌肉瘤没有区别.临床肿瘤学杂志:美国临床肿瘤学会官方杂志 28, 2151–2158. [公共医学][谷歌学者]
  • Yamazaki S、Shen Z、Jiang Y、Smith BJ和Vicini P(2013)。靶向药物配置模型在小分子热休克蛋白90抑制剂中的应用.药物代谢和处置:化学品的生物命运 41,1285–1294。[公共医学][谷歌学者]
  • Zhang J、Benavente CA、McEvoy J、Flores Otero J、Ding L、Chen X、Ulyanov A、Wu G、Wilson M、Wang J等人(2012)。基于基因组和表观遗传学分析的新型视网膜母细胞瘤治疗.自然 481, 329–334.[PMC免费文章][公共医学][谷歌学者]
  • Zhou X、Edmonson MN、Wilkinson MR、Patel A、Wu G、Liu Y、Li Y、Zhang Z、Rusch MC、Parker M等(2016)。使用ProteinPaint探索小儿癌症的基因组改变.自然遗传学 48, 4–6.[PMC免费文章][公共医学][谷歌学者]