跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
乳腺癌研究治疗。作者手稿;可在2018年4月10日的PMC中获得。
以最终编辑形式发布为:
预防性维修识别码:项目经理5892182
尼姆斯:美国国立卫生研究院189761
PMID:20012355

预测新佐剂反应的基因表达途径分析多西他赛联合卡培他滨治疗乳腺癌

关联数据

补充资料

摘要

新辅助化疗已被证明相当于乳腺癌术后治疗,并允许评估化疗反应。在多西他赛(T)和卡培他滨(X)的先导试验中对于第II/III期BC的新辅助化疗,我们评估了基线基因表达和肿瘤对治疗的反应,并检查变化与治疗相关的基因表达。患者接受四个周期的TX.从Mammotome获得的肿瘤组织核心活检TX的预处理(BL)和后循环1(C1)被快速冷冻并储存在加工前为-70°C。利用基因表达分析Affymetrix HG-U133 Plus 2.0基因芯片阵列。进行了统计分析在RMA规范化后使用BRB阵列工具。基因本体(GO)途径随机方差分析t吨具有显著性水平的测试属于P(P)< 0.005. 对于GO确定的基因类别通路分析作为重要的,单个基因的表达水平使用单变量在不同类别之间比较这些路径t吨测验;那些具有显著性水平的基因P(P)< 0.05已报告。对肿瘤样本进行PAM50分析以进行研究生物亚型和复发风险(ROR)。使用GO通路分析,39个基因应答者和非应答者之间的类别区别,最显著的是参与微管组装和调控的基因。比较pre和化疗后标本,我们鉴定出71个差异表达基因包括DNA修复和细胞增殖调控。当时有化疗一个周期后表达变化的45条GO途径应答者和非应答者之间存在显著差异。大多数肿瘤样本分为基底样和管腔B类。ROR得分化疗反应减少;这种变化在样本中更加明显根据临床标准分类为应答者的患者。GO通路分析确定了一些与治疗反应相关的基因类别,以及这可能是一种信息丰富的方法,用于识别对化疗。使用此处描述的方法进行更大的研究是必要的充分评估化疗反应中的基因表达变化。

关键词:新佐剂、基因表达、化疗反应、微管、DNA修复、PAM50

介绍

过去几十年来,乳腺癌死亡率可能有所下降由于乳房X光检查和系统治疗的改进[1]. 实质性的在美国,被诊断患有乳腺癌的女性比例辅助化疗[2],但是我们能够先验地预测哪些患者最有可能受益化疗是有限的。肿瘤组织分子图谱研究进展建议可能采用更个性化的方法治疗乳腺癌[,4]. 使用术前或新佐剂的试验化疗在评估以下因素的分子相关性方面具有独特的地位化疗敏感性,因为肿瘤对治疗的反应可以实时评估。这些研究还有一个额外的优点,即允许评估成对的化疗前后肿瘤标本,以便更好地了解生物学以及乳腺癌的遗传异质性。这反过来会导致更多个体化治疗方法——允许特定选择一种或多种特定肿瘤的化疗药物易受影响。

表达改变的基因很可能是对化疗在那些对治疗有反应的人(即基因)中是不同的赋予敏感性)和不赋予敏感性(即赋予抗性的基因)。在这项研究中,我们试图确定与化疗敏感性相关的基因和阻力。我们利用了第二阶段研究中的治疗前和治疗后活检新辅助多西紫杉醇和卡培他滨化疗的相关性评估在基线基因表达和肿瘤治疗反应之间与治疗相关的基因表达变化。我们还试图评估其他研究中确定的基因表达谱预测数据集中的响应。

材料和方法

这项II期新辅助治疗研究的设计之前描述[5]. 简单地说,女人们新诊断的2或3期乳腺癌,肿瘤大小≥2cm,正常的实验室参数和Zubrod性能状态从0到2为符合条件。如果患者有出血障碍、心脏病,或怀孕或哺乳。对方案进行了审查并经国家癌症研究所机构审查委员会批准。

所有患者均接受了四个周期的多西他赛和卡培他滨治疗每21天。患者最初接受多西他赛(75 mg/m2静脉注射)和卡培他滨(1000 mg/m2p.o.)每天两次每21天2–15次,共四个周期。由于第一阶段毒性过大10名按方案治疗的患者,两种药物的剂量均减少(多西紫杉醇降至60毫克/米2卡培他滨至937.5 mg/m2每天两次)。之后在完成四个周期的多西他赛/卡培他滨(TX)治疗后,所有患者都接受了四次阿霉素循环(60 mg/m2)和环磷酰胺(600毫克/米2)第1天和每21天;6名反应不佳的患者首次化疗前接受阿霉素和环磷酰胺(AC)手术后,其余患者在手术后接受了额外的治疗。所有肿瘤本分析中使用的样本是在AC化疗前采集的用于评估反应的测量也在AC放射治疗之前进行根据手术类型、肿瘤大小和淋巴结受累情况进行个体化治疗。激素受体阳性疾病患者接受他莫昔芬和/或芳香化酶抑制剂完成化疗、手术和放射治疗后。本研究中没有患者接受佐剂或新佐剂曲妥珠单抗。

我们使用双向和最长维度测量记录肿瘤大小在基线检查时、完成TX后和手术前进行体检。我们根据以下方面的变化将患者分为有应答者和无应答者临床检查和病理反应的肿瘤大小(参见表1). 对于术前接受AC治疗的患者,在完成TX时评估反应。简言之完全缓解(三名患者)、手术时的微创疾病(两名患者),或在四个TX周期(三名患者)后出现临床完全反应被视为响应者。所有归类为应答者的患者均接受AC治疗术后。13名患者被视为无应答者;11,带部分缓解(TX四个周期后17–75%残留疾病)四年后,两人患有进展性疾病(基线的110-154%)TX周期。其中,进展性疾病患者和4名患者部分缓解的患者接受了阿霉素的额外术前治疗然而,肿瘤测量用于确定反应在进行额外治疗之前进行分类。包括所有21名患者基线分析中,14名患者的基线和周期1后配对成对分析中包括可评估样本(详细描述(见下文)。

表1

肿瘤特征、临床反应、PAM50的生物亚型和分析中患者肿瘤标本的复发(ROR)评分

零件编号。诊断阶段受体状态基线时的临床肿瘤测量(厘米)四个周期后的临床肿瘤测量TX(厘米)手术病理结果纠正临床反应响应类别基线PAM50学龄后PAM50基线ROR学龄后ROR
响应者
18b条T2N1M0型ER−/PgR+/HER2−  7  01.5厘米R(右)基础正常46−2
27b条T3N1M0型ER+/PgR+/HER2型+  7  1IDC微焦点公共关系R(右)赫兹2微安46  1
25b条T3N0M0型ER+/PgR+/HER2−  5.3  1.51.2厘米c(c)公共关系R(右)亮度A微安29  8
12b条T3N1M0型ER+/PgR+/HER2−120个无残留疾病R(右)亮度B微安49−1
14b条T2N0M0型ER+/PgR−/HER2+  4.5  0无残留疾病R(右)亮度B微安4828
6b条T2N0M0型ER+/PgR+/HER2−  3.5  0残留病灶大小不文件R(右)亮度B微安82  8
T2N1M0型ER+/PgR−/HER2+  5  0.8DCIS;少焦点IDC公共关系R(右)流明B不适用43不适用
21T2N0M0型ER−/PgR+/HER2−  3.50个1.2厘米R(右)基础不适用48不适用
无响应者
8b条T3N1M0型ER−/PgR−/HER2−  8.5  33厘米公共关系尼泊尔卢比玄武岩基础7570
29b条T3N1M0型ER−/PgR+/HER2+  7.5  3.53.5厘米c(c)公共关系尼泊尔卢比基础她的26969
20b条T3N0M0型ER−/PgR−/HER2−12  92.0厘米c(c)公共关系尼泊尔卢比基础基础5648
13b条T3N1M0型ER−/PgR+/HER2−  5.5  8.53.5厘米c(c)PD公司尼泊尔卢比基础基础7520
19b条T3N3M0型ER−/PgR−/HER2−16115.5厘米c(c)公共关系尼泊尔卢比卢马正常−540
26b条T4N0M0型ER+/PgR+/HER2−  53.53.5厘米c(c)公共关系尼泊尔卢比亮度B微安3932
15b条T3N1M0型ER+/PgR+/HER2型+16  58厘米公共关系尼泊尔卢比亮度B卢马4817
17b条T3N1M0型ER−/PgR−/HER2−  8  2.21.2厘米公共关系尼泊尔卢比正常基础  452
4T3N0M0型ER+/PgR−/HER2−  6  11厘米公共关系尼泊尔卢比正常不适用24不适用
24T4N0M0型ER−/PgR−/HER2−13.6  42.5厘米公共关系尼泊尔卢比基础不适用55不适用
16T3N1M0型ER−/PgR−/HER2−10111厘米c(c)标准偏差尼泊尔卢比基础不适用72不适用
7T3N0M0型ER+/PgR−/HER2−  62.6在整个乳房公共关系尼泊尔卢比微安不适用18不适用
28T3N0M0型ER+/PgR+/HER2+  7.5  3.9浸润小叶的散在区域c(c)公共关系尼泊尔卢比微安不适用17不适用
如果雌激素受体和/或孕激素受体阳性>IHC染色阳性的肿瘤细胞占1%;肿瘤如果IHC为3+或FISH检测HER2:CEP17比值>2.0
b条基线检查和配对化疗前后患者样本分析
c(c)患者术前接受AC化疗

评估肿瘤的雌激素(ER)和孕激素受体(PgR)免疫组化6F11抗体对ER和克隆的阳性反应PgR为16(美国新墨西哥州)。如果出现以下情况,则认为肿瘤ER或PgR阳性>1%的细胞呈阳性。如果出现以下情况,则认为肿瘤HER2阳性他们或者被阅读病理学家归类为3+免疫组织化学(DAKO HercepTest)或如果HER2:CEP17比值通过荧光原位杂交(FISH)>2.0。

组织收集和处理

肿瘤组织由Mammotome获得活检前在多西他赛/卡培他滨化疗一个周期后最终手术标本(OR)。本报告仅使用基线和周期1后的样本,以及基因表达结果分析不受阿霉素/环磷酰胺摄入的影响。用液氮快速冷冻组织并在−70°C下保存直到处理。分析中包含的所有样本均为阳性福尔马林固定液上苏木精和伊红染色法检测恶性细胞活检标本(29%)、触摸准备(6%)或两者兼而有之(65%)。

使用以下方法从冷冻组织中提取总RNA三唑®试剂(Life Technologies,Inc.)依据标准协议。RNA的数量和质量由安捷伦科技2100生物分析仪(生物测定软件版本A.02.01。,安捷伦科技公司)和使用纳米液滴®(纳米滴技术)。样本包含使用RNeasy迷你套件(Qiagen)或RNeasy-微型套件(Qaagen),取决于样品大小。作为此操作的一部分,执行了“on column”摘要使用RNase-Free DNA酶组(Qiagen)进行纯化步骤。

用该制剂合成了双标记cDNA基因芯片®表达3′-双周期扩增符合Affymetrix协议的cDNA合成试剂盒基因芯片®真核靶点制备。产生的结果使用基因芯片纯化双链cDNA®样品清理模块(Qiagen)。生物素标记cRNA的合成通过in获得纯化模板cDNA的体外转录基因芯片®表达3′扩增试剂用于IVT标记(Affymetrix)。通过使用基因芯片®样本清理模块(Qiagen)。A 1:10稀释用1μl纯化的IVT反应产物制成,然后继续运行安捷伦科技2100生物分析仪(生物测定软件版本A.02.01。,安捷伦科技),以确定样品质量和纳米液滴®用于定量。然后调整cRNA根据Affymetrix协议,在使用基因芯片®样本清理模块(Qiagen)。一微升从裂解反应中移除,并在安捷伦上进行检查Technologies 2100生物分析仪(生物测定软件版本A.02.01.,安捷伦技术),以在杂交之前确认适当的片段。

制备了每个裂解反应的杂交鸡尾酒根据Affymetrix基因芯片®协议。这个杂交鸡尾酒应用于Affymetrix HU-U133 Plus 2.0基因芯片®数组,在昂飞®Fluidics Station 400随后在Agilent GeneArray扫描仪,配备Affymetrix Microarray Suite 5.0.0.032版软件。

所有分析都是在不知道任何类别分配的情况下进行的(应答者或非应答者)或获取肿瘤样本的时间点。

统计分析

对探针水平的数据进行标准化,并对基因表达总结进行使用稳健多芯片分析(RMA)计算每个探针集。统计学使用BRB阵列工具3.4.1版进行分析。所有基因表达摘要强度(基因摘要)低于50的阈值为50个基因的变异性明显小于中间基因筛选出变异性(基于方差的筛选标准)(P(P)<0.01,基于χ2测试为在BRB数组工具中实现[6]).

使用凝聚层次聚类和使用以平均值为中心的24224个基因通过基于方差的筛选中心相关作为距离度量,平均关联作为链接功能。

BRB阵列工具中实现的基因集表达用于确定具有更多差异表达基因的基因本体(GO)组意外地超过了预期[1]应答者与非应答者[2]基线与周期1后的样本。基因Set Expression方法确定哪些基因集包含更多差异使用两种统计数据(LS和堪萨斯州)。我们认为GO类别非常重要P(P)LS或KS置换测试的值为0.005或更少。通过基于方差的筛选标准的24224个基因中的每一个可以映射到多个GO类别。总的来说,1666个GO类别考虑过的。使用以下方法比较周期1后和基线标本配对数据。

此外,我们确定了GO类别和改变的基因从基线检查到周期1后,in的表达在响应者与非响应者。为了进行此分析我们计算了周期1后基因和基线基因之间的对数比率摘要(如上所述设置阈值)并比较应答者和非应答者;19165个基因通过了基于方差的筛选(适用于对数比率和这个简化的样本集)和1514个唯一GO考虑了类别。

为了评估其他研究中确定的基因是否可以预测样本中的反应,我们匹配了基因图谱中的92个基因由Chang等人出版[7]基于Locuslink ID的阵列探测集,识别HG-U133 Plus 2.0芯片上的251个探头组。类预测器反应是根据我们的数据及其性能使用这组基因开发的因为预测模型必须在我们的数据(只有先前研究的基因列表保持不变),其预测为了避免偏差,使用leave-on-out交叉验证来评估准确性。

最后,使用PAM50算法对患者样本进行评估。对于该分析,来自基线和周期1后样本的原始数据为使用RMA归一化进行处理,每个基因以训练集为中心中位数用于解释阵列平台差异。然后将探针映射到使用Affymetrix注释文件输入基因名称。映射到的探针相同的Entrez基因名称被折叠到样本平均值。50个基因提取了兴趣(列在{Parker,2009#30}中)每个样本和每个亚型之间的Spearman等级相关性计算质心。每个样本都被分配到相关质心。Parker等人概述的算法[8]被跟踪分配生物亚型和复发风险(ROR)评分。

结果

患者和肿瘤特征

患者在2001年1月至2003年8月期间累积到研究中。30名患者参与了本研究;人口统计学和肿瘤特征之前已经描述了整个研究人群[5]. 一名患者在第一个化疗周期,并从所有分析中排除。21岁患者进行了含有恶性细胞的基线肿瘤活检纳入应答者与对照组的基线基因表达分析。无响应者(参见表1). 其中,14名患者在基线检查时和一个周期后接受了可评估的肿瘤样本化疗,并纳入后续分析。

获得的所有肿瘤样本的聚类树状图如所示图1一般来说,来自同一个不同时间点的患者聚集在一起。此外,三阴性肿瘤倾向于聚集在一起,激素受体阳性肿瘤也聚集在一起。

保存图片、插图等的外部文件。对象名为nihms189761f1.jpg

患者标本的聚类树状图。使用对样本进行聚类凝聚层次聚类与以平均值为中心的24224个基因传递以中心相关为距离度量的基于方差的筛选平均值作为联动函数。“ID”列中的数字对应患者ID号;每个病人最多有三个样本

基因本体(GO)类别分析

我们比较了应答者和无应答者,以及基线与周期1后样本,使用GO分类,其根据其相关的生物过程来描述基因产物,细胞成分和分子功能。基线基因分析表达包括来自21名患者的基线样本:8名应答者和13名无应答者。包括配对基线的基因表达变化分析14例患者的1个周期后标本;这六名患者中有应答者和8名患者为无应答者。

当比较有应答者和无应答者时,1666个GO类别中有39个测试的类别在显著性水平上区分了两组属于P(P)< 0.005. 平均来说,我们预计会找到8个如果GO类别与类。差异显著的GO类别包括在中列出表2.单个基因在这些表现出差异表达的类别中(P(P)<0.05)列在附录1这些基因包括许多参与调控的基因微管解聚(微管相关蛋白-2;地图-2)和细胞周期阻滞(微管-肌动蛋白交联因子1;MACF1型). 其他生物相关基因差异表达的包括血管内皮生长因子-B(VEGF-B(血管内皮生长因子-B))和表皮生长因子受体(表皮生长因子受体).

表2

GO类别在应答者和基线肿瘤样本中的无应答者

GO描述GO术语基因数量LS排列P(P)价值KS置换P(P)价值
溶酶体膜科科斯群岛  51.00电子–050.14641
细胞色素-b5还原酶活性MF公司  5东经1.00–040.0161
钙调素结合MF公司520.000190.00129
铜离子结合MF公司170.000420.00034
氧化还原酶活性,作用于NADH或NADPH公司MF公司330.000440.00444
L-氨基酸转运体活性MF公司  90.000820.09359
铜离子转运活性MF公司  60.001460.00015
氨基酸转运体活性MF公司400.001710.0191
组织形态发生英国石油公司120.002110.0094
胺转运蛋白活性MF公司440.002340.01258
蛋白质自动处理/蛋白质氨基酸磷孢子化英国石油公司100.002740.00069
脂溶性维生素代谢英国石油公司  50.002970.04321
细胞周期阻滞英国石油公司320.002970.03979
二酰甘油结合MF公司160.003440.01502
谷胱甘肽生物合成英国石油公司100.003520.00481
角化作用英国石油公司  60.003740.08781
磷酸肌醇介导的信号转导英国石油公司230.004170.00126
L-氨基酸转运英国石油公司  60.004180.20218
微管聚合或解聚英国石油公司  90.004310.09298
1-酰基甘油-3-磷酸-酰基转移酶活性MF公司140.00440.18002
G蛋白信号耦合到IP3秒信使(磷脂酶C激活)英国石油公司210.004540.00332
细胞连接科科斯群岛870.004820.03406
表皮形态发生英国石油公司7个0.004880.09918
酰基甘油-酰基转移酶活性MF公司160.004950.06529
蛋白质过氧化物酶体靶向英国石油公司  50.005330.00151
表皮生长因子受体信号通路英国石油公司220.006220.0018
中性脂质代谢英国石油公司  80.012930.00068
消化英国石油公司160.013450.00182
单加氧酶活性MF公司100.015450.00112
钙和钙调素依赖蛋白激酶复合物科科斯群岛  90.016080.00274
滑面内质网科科斯群岛  80.020460.00191
肝素结合MF公司290.025930.00443
乙醇胺磷酸转移酶活性MF公司  50.028610.00132
核定位序列绑定MF公司  90.036310.00399
水运输活动MF公司  50.07180.00183
乙醇代谢英国石油公司  50.096960.00482
肽激素分泌英国石油公司  50.111540.00211
核染色质科科斯群岛  90.231620.0036
髓样细胞分化英国石油公司  80.337650.00446

科科斯群岛蜂窝组件,MF公司分子的功能,英国石油公司生物过程

当比较基线肿瘤标本和一次后获得的标本时多西他赛和卡培他滨化疗周期,1666个GO类别中的71个P(P)< 0.005 (表3). 平均来说,我们希望偶然发现8个差异表达的GO类别。个人与这些途径表达显著不同的基因有在中列出补充附录2.显著差异表达的基因包括参与DNA修复的人[乳腺癌1例早期发病;巴西航空公司1,聚(ADP-核糖)聚合酶家族,成员2和3;第2部分第3部分)]、和单元格增殖调节[成纤维细胞生长因子受体癌基因合伙人(FGFR1OP公司)和表皮生长因子受体(表皮生长因子受体)]. 此外,一些纺锤体相关检查点基因表现出差异表达,包括苯并咪唑蛋白1和3不抑制芽接(BUB1(业务单元1)业务单元3),有丝分裂停滞缺陷蛋白2L1(MAD2L1型)和微管蛋白γ亚型。

表3

GO类别在一个周期后显著差异表达化疗与基线比较

GO描述GO术语GO类别中的基因数LS置换P(P)价值KS置换P(P)价值
染色体中心周围区科科斯群岛  26东经1.00–05东经1.00–05
动粒科科斯群岛  16东经1.00–05东经1.00–05
浓缩染色体科科斯群岛  18东经1.00–050.00175
主轴科科斯群岛  67东经1.00–050.00265
染色体部分科科斯群岛  33东经1.00–052005年9月
层粘连蛋白结合MF公司    5东经1.00–05东经1.00–05
有丝分裂细胞周期M期英国石油公司  69东经1.00–05东经1.00–05
有丝分裂细胞周期英国石油公司88个东经1.00–05东经1.00–05
M相英国石油公司  79东经1.00–05东经1.00–04
葡萄糖分解代谢英国石油公司  80东经1.00–050.01684
DNA复制英国石油公司  98东经1.00–05东经1.00–05
DNA修复英国石油公司  92东经1.00–05东经1.00–05
有丝分裂英国石油公司  64东经1.00–05东经1.00–05
单糖分解代谢英国石油公司  83东经1.00–050.01815
复制体科科斯群岛15个2005年9月0.00108
依赖DNA的DNA复制英国石油公司44个2005年9月0.00924
复制分叉科科斯群岛  16东经1.00–040.00287
糖酵解英国石油公司  63东经1.00–040.0016
DNA复制起始英国石油公司    6东经1.00–040.00016
浓缩核染色体科科斯群岛  140.000130.01502
Rho GTPase激活剂活性MF公司  140.000170.00562
碳水化合物分解代谢英国石油公司1000.000310.04332
线粒体腔科科斯群岛  160.000310.03546
酒精分解代谢英国石油公司  860.000310.03545
己糖分解代谢英国石油公司  820.000320.0261
纺锤组织和生物成因英国石油公司  110.000360.00177
有丝分裂的调节英国石油公司  170.000360.00089
细胞增殖的正向调节英国石油公司  490.000381.00电子–05
染色质结合MF公司  44东经4.00–040.01305
染色体分离英国石油公司  120.000480.00792
核仁科科斯群岛  990.000650.05974
细胞周期的调节英国石油公司  190.000860.0036
线粒体基质科科斯群岛  640.000910.01216
对辐射的反应英国石油公司  110.000950.00032
泛醌代谢英国石油公司6个0.001020.03387
对温度的响应英国石油公司9个0.001490.00069
对寒冷的反应英国石油公司    60.001540.00086
生物聚合物甲基化英国石油公司  220.001680.00038
氧化还原辅酶代谢英国石油公司  250.001940.25704
有丝分裂纺锤体组织和生物成因英国石油公司  100.001960.00456
多聚体科科斯群岛    90.00198东经1.00–05
细胞的正向调节新陈代谢英国石油公司  550.002110.03108
辅酶生物合成英国石油公司  610.002510.00018
转换启动英国石油公司  630.002720.00303
线粒体部分科科斯群岛  220.003030.20335
RNA聚合酶II转录因子活性增强剂结合MF公司    70.003620.08965
核染色体科科斯群岛  420.003760.05191
葡萄糖输入调节英国石油公司    70.003980.01164
细胞分裂英国石油公司  800.00410.0135
蛋白质氨基酸甲基化/烷基化英国石油公司  150.004610.00476
微管细胞骨架组织和生物成因英国石油公司  510.004790.0234
细胞命运承诺英国石油公司  200.006320.00421
细胞定位英国石油公司  380.007020.00355
微管运动英国石油公司84个0.009280.00048
谷胱甘肽生物合成英国石油公司10个0.00970.00019
翻译起始的调节英国石油公司  290.015620.00011
磷酸丙酮酸水合酶复合物科科斯群岛    50.017820.00395
周心材料科科斯群岛    70.017890.00261
细胞骨架依赖性细胞内运输英国石油公司  900.019190.00122
谷胱甘肽代谢英国石油公司  130.02492东经6.00–04
mRNA 3′-UTR结合MF公司    50.033750.00172
RNA代谢的调节英国石油公司  130.034790.00275
氧转运蛋白活性MF公司    70.05040.00431
蛋白质聚合英国石油公司46个0.054020.0012
辅因子转运体活性MF公司7个0.058270.00454
应激激活蛋白激酶信号通路/JNK级联英国石油公司  160.09780.00475
中胚层发育英国石油公司  140.104020.00199
细胞发育英国石油公司  240.130610.00421
辅因子运输英国石油公司    70.14560.00454
Ras GTPase激活剂活性MF公司  100.171160.00448
不溶部分科科斯群岛  130.227030.00299

我们寻找哪些基因在有应答者和无应答者从基线到第1周期后的表达。在这项分析中,我们发现1514个GO类别中有45个类别(表4)其中一个之后的表达式变化两类患者的化疗周期不同(8个为只是偶然的)。这些类别中的单个基因在这项分析中有统计学上的显著差异,如补充附录3.在应答者发生较大变化的基因中,有两种热休克蛋白(热休克蛋白A9B热休克蛋白5),抑制素β-A,A泛素相关基因数量、肿瘤坏死因子α诱导蛋白6(TNFαIP6)和骨膜炎(邮政).表达变化显著大于无应答者比应答者包括转化生长因子β受体III(转化生长因子βR3),VEGFB公司、和成纤维细胞生长因子受体-1(FGFR1号机组)和相关基因在活性氧代谢中,如谷胱甘肽过氧化物酶-3(GPX3系列)和醛氧化酶-1(AOX1型).

表4

GO类别,其中基因表达从基线到周期后的变化1个化疗方案在应答者和无应答者

GO描述GO术语基因数量LS置换P(P)价值KS置换P(P)价值
内质网腔科科斯群岛19东经1.00–050.02961
图案装订MF公司45东经1.00–05东经1.00–05
电子转运体活性MF公司73东经1.00–050.00539
糖胺聚糖结合MF公司44东经1.00–05东经1.00–05
小蛋白激活酶活性MF公司41东经1.00–050.0392
泛素样激活酶活性MF公司40东经1.00–050.02938
多糖结合MF公司44东经1.00–05东经1.00–05
信号序列绑定MF公司16东经1.00–040.00115
肌动蛋白丝的调节聚合英国石油公司  60.000180.00195
肌动蛋白纤维聚合英国石油公司120.000360.007
对温度的响应英国石油公司  80.000380.00889
蛋白质转锁酶活性MF公司140.000630.02231
蛋白质-线粒体靶向英国石油公司  5东经9.00–040.08933
氧化还原酶活性作用于单个掺入分子氧的供体MF公司  50.001950.0092
蛋白质聚合英国石油公司320.001950.00354
蛋白质载体活性MF公司170.001960.03679
细胞表面科科斯群岛390.002040.03331
中性氨基酸转运体活性MF公司  60.002070.01723
钠:二羧酸同向转运体活性MF公司  80.002160.02185
细胞周期阻滞英国石油公司200.004210.00081
肝素结合MF公司250.00428东经1.00–05
蛋白质进口英国石油公司320.004410.26526
辅因子转运体活性MF公司  60.004460.00871
蛋白质靶向英国石油公司650.004710.19948
caspase活性的调节英国石油公司120.004840.0657
tRNA修饰英国石油公司130.005070.00369
RNA修饰英国石油公司180.005140.00182
循环英国石油公司340.00570.0038
氧气和活性氧种类新陈代谢英国石油公司380.008810.00133
氧化还原酶活性作用于CH–OH捐赠者组MF公司350.009250.00241
蛋白质酪氨酸磷酸酶活性MF公司130.009580.00064
STAT蛋白核移位英国石油公司  70.010290.00432
tRNA代谢英国石油公司470.018820.00227
磷蛋白磷酸酶活性MF公司290.019540.00143
肌肉肌球蛋白科科斯群岛100.020720.0026
肌球蛋白科科斯群岛100.020720.0026
成纤维细胞生长因子受体活动MF公司  50.022930.00435
蛋白激酶A结合MF公司110.030010.00386
COPI囊泡涂层科科斯群岛  70.030330.00301
FACIT胶原蛋白科科斯群岛  80.03366东经1.00–04
神经肌肉接头发育英国石油公司  90.036670.00087
锚定胶原蛋白科科斯群岛  90.044560.00216
G蛋白偶联受体结合MF公司310.04840.0033
心脏收缩率的调节英国石油公司100.051040.00422
骨骼肌纤维发育英国石油公司100.056990.00422

班级预测

为了评估其他研究中确定的基因是否为化疗反应的重要决定因素在我们的样本,我们在数据集上使用Chang等人鉴定的基因[7]并使用leave-on-out评估其绩效交叉验证。基因的起始集是由Chang等人[7]作为应答者和非应答者之间在他们的数据集。在开发预测工具时,我们从Chang中选择基因,即在数据集中的显著性水平为0.1。为了避免在由于分类器使用我们的数据而导致的预测精度估计构建过程中,我们使用了leave-on-out交叉验证。具体来说对每个人重复基因选择和模型拟合的整个过程离场交叉验证培训集,确保公平代表预测精度。

用于响应的复合协变量分类器的灵敏度为通过交叉验证估计50%和61.5%的特异性。阳性和阴性预测值分别为44%和67%。分类正确率的交叉验证估计为57%。使用其他预测方法或进一步的方法获得了类似的结果选择要包含在分类器中的特征。具体来说,我们将分析限制在Chang列表中显示我们数据集中的一致“方向”,即指向那些应答者和非应答者之间的折叠变化均大于或小于1研究,无论数据集的统计意义如何。这个仅使用一致基因构建的分类器的预测准确性为与上述内容类似。

生物亚型和ROR评分

生物亚型分析结果和ROR评分如所示表1基线检查时,21人中的38%分析的肿瘤分类为基底样,29%为管腔B,19%为管腔A,9%为正常,5%为富含HER2。基线检查时,内脏A肿瘤的ROR得分最低。平均值有应答者和无应答者的ROR相似,但在在无应答者中记录到基线ROR(应答者的平均ROR=47,四分位范围=3;无应答者的平均ROR=48,四分位范围=51)。在配对分析中,来自患者的肿瘤根据临床标准分类为应答者往往有很大的从基线检查到第1周期后ROR得分下降[第1周期平均值6名患者的ROR=4.5(四分位间距=8.5)包含在配对分析中];相反,较小的减少或无应答者的ROR评分增加[第1周期后平均值八个样本中的ROR=44(四分位范围=27.25)配对分析中的无应答者]。每个项目的ROR变化单个患者如所示图2.有趣的是,在配对分析中,11/14名患者的生物亚型化疗一个周期后发生变化;最常见的变化是从管腔B至管腔A亚型。在应答者中,所有化疗后标本都映射到管腔A组或正常样组。在无应答者中化疗一个周期后,亚型变化较大。热图显示对所有患者样本(A组)和仅基线样本(B组)进行聚类如所示图3.

保存图片、插图等的外部文件。对象名为nihms189761f2.jpg

有应答者化疗后ROR评分的变化(面板)和无应答者(面板b)

保存图片、插图等的外部文件。对象名为nihms189761f3a.jpg
保存图片、插图等的外部文件。对象名为nihms189761f3b.jpg

来自PAM50的热图使用(a)所有患者的层次聚类分析标本(b)仅基线标本。示例密钥:蓝色管腔的A中,浅蓝色管腔B,粉红色富含HER2,红色类基底,绿色像正常人一样

讨论

在本研究中,通过cDNA微阵列分析乳腺中的基因表达接受新辅助多西他赛和卡培他滨治疗的癌症患者,我们发现在表征反应方面似乎很重要的生物途径的数量对化疗,包括纺锤体调节和微管解聚,DNA修复和细胞增殖。此外,我们发现之前发表的与治疗反应相关的基因集对我们的样品。

多西紫杉醇是一种紫杉烷,用于治疗转移性和辅助设置。紫杉烷与微管蛋白结合并抑制微管解聚,破坏有丝分裂细胞分裂从而导致细胞死亡[9]. 微管相关蛋白,MAP2、MAP4和MAP-tau通过与微管蛋白结合来稳定微管,并增加微管动力学与MAP表达改变有关紫杉烷抗性[10,11]. 在我们的研究中,无应答者显著提高地图2而不是回应者,建议MAP2蛋白升高可能与患者的耐药性有关。地图4也是GO途径中显示在我们的分析中倾向于差异表达,但在应答者和非应答者之间的基因表达<2,并且P(P)单个基因探针集的值对应于该基因无统计学意义。

除了这些基线差异外,我们还发现比较肿瘤标本时微管和纺锤体相关基因的数量化疗一个周期前后。这些包括主轴组件检查点基因BUB1(业务单元1)业务单元3MAD2L1型以及一些微管蛋白亚型。这些基因可能起作用在确定对含有紫杉烷的反应中也起着重要作用化疗。

我们还注意到一些基因的表达显著增加参与一个化疗周期后的DNA修复途径,包括巴西航空公司1,第2部分、和CHEK1(检查1).DNA损伤是许多具有活性的化疗药物的作用机制乳腺癌和DNA修复途径的改变在乳腺癌的发病机制巴西航空公司1巴西航空公司2突变[12]. 此外,一类新的靶向药物,PARP目前正在对抑制一种参与DNA修复的关键酶的抑制剂进行评估两者都是针对BRCA公司-相关癌症及其合并症DNA破坏剂[13,14]. 已知两者之间的重叠BRCA1突变和三阴性肿瘤,这些药物也在三阴性疾病的设置(http://www.cancer.gov/search/ResultsClinicalTrialsAdvanced.aspx?protocolsearchid=5122624).大约三分之一的患者患有三重阴性肿瘤在我们的研究中,应答者的百分比患有三阴性疾病。数据来自大型新辅助化疗研究表明,患有三阴性肿瘤的女性与荷尔蒙组相比,更有可能实现病理性完全反应敏感肿瘤,但那些没有反应的患者预后较差[15]. 因此,如果DNA修复在这些肿瘤的发病机制以及对如我们的数据所示,治疗后这些基因在诊断时间和化疗早期可能都是预后因素和预测。

FGF途径中的基因最近被显示为全基因组关联与乳腺癌风险相关的研究[16,17],FGF的肿瘤表达被认为是导致化疗耐药的原因,特别是对紫杉醇的耐药性[18,19]. 在我们的分析中,的表达式FGFR1OP公司一个周期后显著增加化疗。该基因的产物是一种中心体蛋白,涉及在微管的锚定中,并被认为是一种预后指标非小细胞肺癌[20]. 也与这个假设一致,FGFR1号机组在我们的分析中,应答者的基因表达下降。

被认为参与决定对卡培他滨的反应性的基因,例如胸苷合成酶、胸苷磷酸化酶和二氢嘧啶脱氢酶[21]没有在我们的研究中发挥着关键作用。Pusztai等人最近建议第二阶段的研究不太可能识别出显著的单个基因与对个别治疗药物的反应有关[22]. 一个明显的挑战是缺乏权力在小型研究中,再加上多重比较的问题很难区分真正差异表达的基因和噪声。另外,在这些小型研究中,肿瘤之间的显著异质性可能导致稀释可能明显的基因表达的真正差异如果对更均匀的肿瘤进行研究。我们的患者人数相当多肿瘤的分子和组织学特征的异质性。越来越明显的是,乳腺癌的不同亚类具有非常不同的预测,并且可能对特定类别有不同的反应治疗药物。因此,对DNA修复表达变化的分析仅限于三阴性肿瘤的基因可能与确定结果比涉及混合肿瘤组的结果要好。

对我们基线标本的PAM50分析显示ER和PgR表达和亚型分类之间的一致性IHC表达ER或PgR的样本往往属于管腔类型,而三阴性肿瘤表达基本表型。此外,尽管根据PAM50对化疗前后患者的分析,这些数字非常小标本引发了几个有趣的问题。首先,我们发现ROR得分化疗一个周期后总体下降。这可能反映了细胞周期因化疗而被捕。此外,我们还看到,经过一个循环化疗后,我们大多数标本的指定表型发生了变化最常见的变化是从管腔B型变为管腔A型。这些数据可以反映了肿瘤细胞增殖的减缓,使其与管腔A亚型。或者,我们的数据可以反映选择性杀戮对化疗敏感的肿瘤细胞,会留下亚群激素敏感、化疗反应较低的肿瘤细胞。虽然我们的结果是从生物学角度来看,这些数据必须谨慎解释,因为样本量和确定亚型成员的方法学挑战新数据集。在解释分析数据时必须谨慎,这意味着使用基因中心化,特别是在小样本和以下情况下研究人群中的亚型流行率与最初用于创建算法的示例[23]. 重申我们使用RT-PCR的样品可能已经产生了不同的结果,并可能对未来的研究有所帮助。因此,更大需要对可用于RT-PCR的样品进行研究,以确定化疗的分子亚型是临床上的真实现象启示。

如前所述,使用高通量分析的小型II期研究技术可能受到两个重要问题的限制:高错误发现率对大量基因的询问,以及基因表达导致识别重要个体的能力不足基因。在我们的分析中,我们试图通过寻找与生物相关的基因路径和类别,而不是单个基因,并使用之前设计的基于已证明与结果相关的多个基因[8]. 我们的GO途径分析确定39个GO类别区分应答者和非应答者,71个一个化疗周期后表达差异显著的类别在1666个类别中,包含符合阈值的基因条件。使用P(P)显著性小于0.005,我们预计只有八个类别会对每个类别产生显著差异这些分析纯属偶然。

对小型第二阶段研究的识别能力持怀疑态度化疗反应的药物基因组标记[22]当然是有保证的。然而,我们相信使用生物相关途径而非个体进行分析基因是一种可行的策略,尤其是对于有限的研究检测单个基因表达显著差异的敏感性。因此,通路方法代表了基因领域的一个重要进展表达分析。

总的来说,我们的研究有几个明显的优势。新鲜的可用性在治疗过程中的多个时间点收集的冷冻组织是一种稀有而宝贵的资源。虽然在所有研究中我们没有配对样本参与者,我们不认为我们的研究样本有偏见。似乎没有化疗前后标本与肿瘤是否相关大小、肿瘤特征或治疗反应。

然而,应注意几个重要的限制,包括样本量和我们的发现缺乏验证集。因此,数据这里的介绍必须被认为是探索性的。当我们进入个性化时代肿瘤治疗,获得高质量肿瘤样本的研究,以及储存用于研究,特别是用于分子诊断至关重要。此外,在新辅助化疗研究中化疗后可以使用新的分子技术进行检查至关重要促进我们对乳腺癌生物学和个体肿瘤对治疗的反应。

总之,我们得出结论,一些生物学上似乎合理的途径,包括微管和细胞纺锤体的结构和功能,DNA修复,以及生长因子受体可能在确定敏感性和对以紫杉醇为基础的化疗。我们肿瘤标本中的PAM50分析显示合理IHC分类和分子亚型之间的一致性,如其他研究。此外,我们的数据表明,肿瘤ROR评分随着化疗,尤其是对反应性肿瘤的化疗,尽管有更大的研究有必要确认这一发现。使用GO类别进行基因通路分析应该作为限制基因错误发现率的一种手段表达谱研究。我们的发现与来自其他研究人员检查了基因表达与对治疗的反应对于确定哪些特定基因是最有价值的重要。在现有肿瘤的大型研究中证实了这些发现标本有望使我们更接近实现个性化的目标乳腺癌的治疗。

补充材料

1

单击此处查看。(224K,pdf)

致谢

这项研究得到了NIH和国家癌症研究所。我们要感谢马修·埃利斯博士和查尔斯·E。佩罗感谢他们对PAM50的分析和解释提供的帮助。我们也会感谢参与本研究的患者的参与。

脚注

电子辅助材料本文的在线版本(doi:10.1007/s10549-009-0651-3)包含补充材料,可供授权人员使用用户。

参与者信息

拉里萨·科尔德,美国医科大学肿瘤医学系华盛顿/西雅图癌症护理联盟,G3639,825 Eastlake Ave,E,西雅图,WA 98109,美国。

劳拉·卢萨,美国大学生物统计与医学信息研究所卢布尔雅那,卢布尔雅那,斯洛文尼亚。

丽莎·麦克谢恩,癌症治疗和美国贝塞斯达NCI诊断(DCTD)。

Peter F.Lebowitz,葛兰素史克公司(Glaxo Smith Kline),美国Collegeville。

LuAnne Lukes,美国贝塞斯达NCI人口遗传学实验室。

凯文·坎普豪森,NCI癌症研究中心放射肿瘤科,美国贝塞斯达。

乔尔·帕克,Lineberger综合癌症中心遗传学系,美国北卡罗来纳大学教堂山分校。

桑德拉·斯温,华盛顿医院中心,美国华盛顿。

肯特·亨特,美国贝塞斯达NCI人口遗传学实验室。

乔·安妮·祖耶夫斯基,癌症治疗评估项目临床调查处,美国贝塞斯达NCI DCTD。

工具书类

1Berry DA、Cronin KA、Plevritis SK、Fryback DG、Clarke L、Zelen M等。筛查和辅助治疗对乳腺死亡率的影响癌症。N英格兰医学杂志。2005;353(17):1784–1792.[公共医学][谷歌学者]
2Harlan LC、Clegg LX、Abrams J、Stevens JL、Ballard-Barbash R.化疗和激素治疗的社区应用早期乳腺癌:1987年至2000年。临床肿瘤学杂志。2006;24(6):872–877.[公共医学][谷歌学者]
三。Andre F、Mazouni C、Hortobagyi GN、Pusztai L.DNA阵列作为佐剂/新佐剂疗效的预测因子乳腺癌患者的化疗:当前数据和研究问题设计。Biochim生物物理学报。2006;1766(2):197–204.[公共医学][谷歌学者]
4Chuthapitish S,Eremin JM,Eremin O.预测乳腺癌新辅助化疗的反应:分子影像学、全身生物标志物与肿瘤代谢组(审查)Oncol代表。2008;20(4):699–703.[公共医学][谷歌学者]
5Lebowitz PF、Eng-Wong J、Swain SM、Berman A、Merino MJ、Chow CK等。新佐剂多西紫杉醇和卡培他滨治疗局部晚期乳腺癌。临床癌症研究。2004;10(20):6764–6769.[公共医学][谷歌学者]
6Simon R,Lam A,Li MC,Ngan M,Menenzes S,Zhao Y.使用BRB阵列分析基因表达数据工具。癌症信息。2007;2:11–17. [PMC免费文章][公共医学][谷歌学者]
7Chang JC、Wooten EC、Tsimelzon A、Hilsenbeck SG、Gutierrez MC、Elledge R等。预测治疗效果的基因表达谱乳腺癌患者对多西紫杉醇的反应。柳叶刀。2003;362(9381):362–369.[公共医学][谷歌学者]
8Parker JS、Mullins M、Cheang MC、Leung S、Voduc D、Vickery T等。基于内在因素的乳腺癌风险预测值监测子类型。临床肿瘤学杂志。2009;27(8):1160–1167. [PMC免费文章][公共医学][谷歌学者]
9Eisenhauer EA、Vermorken JB。分类单元。临床药理学与治疗学的比较潜力。药物。1998;55(1):5–30.[公共医学][谷歌学者]
10Goncalves A、Braguer D、Kamath K、Martello L、Briand C、Horwitz S等。肺癌细胞对紫杉醇的耐药性与微管动力学增强。美国国家科学院程序。2001;98(20):11737–11742. [PMC免费文章][公共医学][谷歌学者]
11McGragan BT、Gilmartin B、Carney DN、McCann A.紫杉烷、微管和耐药乳腺癌症。Biochim生物物理学报。2008;1785(2):96–132.[公共医学][谷歌学者]
12Welsh PL、King MC、BRCA1和BRCA2与乳腺和卵巢的遗传学癌症。人类分子遗传学。2001;10(7):705–713.[公共医学][谷歌学者]
13Farmer H、McCabe N、Lord CJ、Tutt AN、Johnson DA、Richardson TB等。以BRCA突变细胞中的DNA修复缺陷为靶点治疗策略。自然。2005;434(7035):917–921.[公共医学][谷歌学者]
14Martin SA、Lord CJ、Ashworth A.DNA修复缺陷是癌症。当前操作基因开发。2008;18(1):80–86.[公共医学][谷歌学者]
15Wolff AC、Berry D、Carey LA、Colleoni M、Dowsett M、Ellis M等。影响术前系统治疗的研究问题可手术的乳腺癌。临床肿瘤学杂志。2008;26(5):806–813.[公共医学][谷歌学者]
16Easton DF、Pooley KA、Dunning AM、Pharoah PD、Thompson D、Ballinger DG等。全基因组关联研究确定新的乳腺癌易感基因座。自然。2007;447(7148):1087–1093. [PMC免费文章][公共医学][谷歌学者]
17Hunter DJ、Kraft P、Jacobs KB、Cox DG、Yeager M、Hankinson SE等。全基因组关联研究确定FGFR2的等位基因与散发性绝经后乳腺风险相关癌症。自然遗传学。2007;39(7):870–874. [PMC免费文章][公共医学][谷歌学者]
18Gan Y、Wientjes MG、Au JL。碱性成纤维细胞生长因子的表达与人类肿瘤患者对紫杉醇的耐药性。药物研究。2006;23(6):1324–1331.[公共医学][谷歌学者]
19Song S、Wientjes MG、Gan Y、Au JL。成纤维细胞生长因子:一种广泛的表观遗传机制对抗癌药物的光谱抗性。美国国家科学院程序。2000;97(15):8658–8663. [PMC免费文章][公共医学][谷歌学者]
20Mano Y、Takahashi K、Ishikawa N、Takano A、Yasui W、Inai K等。成纤维细胞生长因子受体1癌基因伴侣作为新的肺癌的预后生物标志物和治疗靶点。癌症科学。2007;98(12):1902–1913. [PMC免费文章][公共医学][谷歌学者]
21Vallbohmer D、Yang DY、Kuramochi H、Shimizu D、Danenberg KD、Lindebjerg J等。DPD是卡培他滨在结直肠癌。国际癌症杂志。2007;31(2):413–418.[公共医学][谷歌学者]
22Pusztai L、Anderson K、Hess KR。二期临床试验中药物基因组预测因子的发现治疗乳腺癌。临床癌症研究。2007;13(20):6080–6086.[公共医学][谷歌学者]
23Lusa L,McShane LM,Reid JF,De Cecco L,Ambrogi F,Biganzoli E,et al.跨基因预测聚类结果的挑战表达式分析数据集。美国国家癌症研究所杂志。2007;99(22):1715–1723.[公共医学][谷歌学者]