CMMM 计算和数学方法在医学 1748 - 6718 1748 - 670 x Hindawi 10.1155 / 2021/2649123 2649123 研究文章 转录因子的识别签名可以预测乳腺癌的生存 https://orcid.org/0000 - 0002 - 8009 - 9607 方ydF4y2Ba Chunni 1 https://orcid.org/0000 - 0002 - 1521 - 9756 剑士 2 https://orcid.org/0000 - 0002 - 4433 - 290 x 1 Alexov 埃米尔 1 的乳房手术 吉林大学第三医院 长春 吉林130033 中国 jlu.edu.cn 2 部血管手术 吉林大学第三医院 长春 吉林130033 中国 jlu.edu.cn 2021年 20. 2 2021年 2021年 5 7 2020年 6 1 2021年 27 1 2021年 20. 2 2021年 2021年 版权©2021 Chunni风扇等。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

背景。转录因子的表达模式(TFs)可以用来开发潜在的癌症预后的生物标志物。在这项研究中,我们旨在识别和验证TF签名预测乳腺癌的无病生存期(DFS) (BRCA)患者。<我talic> 方法。套索和Cox回归分析应用于构建一个基于TF签名,TCGA的基因表达数据集。的预后价值TF签名TCGA调查数据库,进一步验证模型的可靠性3独立从基因表达数据集综合(GEO)。TF签名的预测性能比较与先前发表4基因签名。研究癌症的TF签名和特征之间的关系,基因集富集(GSEA)进行了分析。TF的相关性签名和免疫渗透的水平也被调查。<我talic> 结果。11-TF预后签名是由对BRCA病人生存预测性能好。通过使用风险评分模型基于11-TF签名,BRCA患者分层为低收入和高危人群和显示好,可怜的无病生存期(DFS),分别。分数是一个独立的风险预测指标调整其他临床病理的因素。此外,11-TF签名有一个更好的生存预测性能相比4以前公布的基因签名。此外,风险评分是一个癌症的标志。最后,一个高风险得分更高的渗透M0和M2巨噬细胞,降低渗透休息CD4记忆+T细胞、CD8+T细胞。<我talic> 结论。在这个研究结果确认和验证小说预后TF签名,这是一个独立的预测生物标志物DFS BRCA病人。

1。介绍

乳腺癌(BRCA)是死于癌症的主要原因之一,妇女和代表一个异构组肿瘤起源于上皮细胞衬乳导管( 1, 2]。传统上,肿瘤大小;淋巴结的转移状况;和地位的雌激素受体(ER)、孕激素受体(PR),和HER2作为有用的预后标志物的BRCA诊所( 3- - - - - - 5]。然而,这些预后生物标志物仍局限在准确预测生存由于BRCA患者的遗传异质性 6]。因此,建立一个精确的和健壮的签名来指导预后分层BRCA是相当重要的。

转录因子(TFs) dna结合蛋白质与基因的启动子序列和随后调节基因的表达 7]。先前的研究已经表明,TFs参与了几个关键的细胞过程,如细胞生长、分化、增殖和细胞死亡 8- - - - - - 10]。TFs BRCA患者往往异常表达,TFs的表达之间的关系,已经证明了病人生存BRCA [ 11]。然而,TFs的表达模式和预后能力BRCA只在一些研究调查。

在目前的研究中,我们确定了几个survival-related TFs的预测和开发了一个11-TF签名BRCA患者的无病生存期(DFS)通过分析基因表达和相应的临床数据的BRCA患者癌症基因组图谱(TCGA)。TF签名验证的预测性能3独立从基因表达数据集综合(GEO)并与4以前公布的基因签名。此外,协会之间的风险评分基于11-TF签名和癌症的标志,和免疫细胞渗透了。

2。材料和方法

所有数据分析使用R软件( http://www.r-project.org/version3.5.1)。研究过程的流程图如图 1

这项研究的流程图。

2.1。数据收集

2018年,兰伯特等人发表了一篇引人注目的评论覆盖≥1600可能人类TFs潜在的人类生理学和疾病( 7]。在这项研究中,我们收集了1639 TFs文学的进一步探索。从BRCA基因表达数据患者从TCGA通过下载加州大学圣克鲁兹齐娜浏览器(UCSC齐娜, http://xena.ucsc.edu/)[ 12]。TCGA的重叠TFs 1639 TFs被提取,和1554年得到了TFs TCGA的微分表达式分析队列(114配对正常和肿瘤样本)。

共有950名患者相应的基因组数据(读计数)和临床信息包括DFS TCGA数据是从数据库中,作为训练集。基因表达谱和相应的临床信息的BRCA患者选择在线公共基因表达综合(地理, http://www.ncbi.nlm.nih.gov/geo/)数据库,包括三个独立数据集的数据加入GSE20685 [ 13],GSE21653 [ 14, 15],GSE42568 [ 16]。GSE20685由327 BRCA病例,GSE21653由266 BRCA病人,GSE42568包含121个样本(104 BRCA 17患者和正常对照组,表 1)。共有697个肿瘤样本生成使用相同的芯片平台(GPL570 Affymetrix人类基因组U133 + 2.0数组),和这些数据集作为验证集。

Cox回归分析的临床病理因素之间的关系(包括风险分数)和BRCA患者的无病生存。

变量 集团 患者(<我nline-formula> N ) 单变量分析 患者(<我nline-formula> N ) 多变量分析
人力资源(95%置信区间) P 人力资源(95%置信区间) P
,TCGA
风险评分 低/高 475/475 4.18 (2.50 - -6.97) 4.37 E 08年 418/411 4.40 (2.31 - -8.39) 6.62 E 06
年龄 ≤50 / > 50 290/660 0.67 (0.43 - -1.03) 6.59 E 02 251/578 0.66 (0.39 - -1.11) 1.14 E 01
病理阶段 i ii / iii iv 723/216 3.34 (2.15 - -5.19) 8.81 E 08年 645/184 4.28 (2.49 - -7.34) 1.31 E 07年
ER状态 N / P 214/697 0.49 (0.31 - -0.77) 1.78 E 03 191/638 1.35 (0.32 - -5.61) 6.82 E 01
公共关系状态 N / P 304/604 0.46 (0.30 - -0.71) 5.34 E 04 274/555 0.56 (0.26 - -1.19) 1.29 E 01
HER2状态 N / P 680/160 0.93 (0.47 - -1.82) 8.22 E 01 671/158 0.72 (0.30 - -1.69) 4.47 E 01
三负 没有/是的 804/146 1.62 (0.97 - -2.71) 6.39 E 02 685/144 1.89 (0.43 - -8.35) 4.02 E 01
GSE20685
风险评分 低/高 164/163 2.41 (1.58 - -3.68) 4.16 E 05年 164/158 2.22 (1.43 - -3.43) 3.55 E 04
年龄 ≤50 / > 50 209/118 0.74 (0.48 - -1.14) 1.71 E 01 206/116 0.79 (0.49 - -1.27) 3.36 E 01
T台 i ii / iii iv 101/188 1.93 (1.50 - -2.49) 3.27 E 07年 100/186 1.22 (0.86 - -1.74) 2.57 E 01
N阶段 N0 / N1-3 137/190 3.73 (2.26 - -6.15) 2.49 E 07年 135/187 3.46 (1.99 - -6.01) 1.06 E 05年
M阶段 M0、M1 319/8 21.96 (10.49 - -45.99) 2.59 E 16 315/7 5.78 (2.02 - -16.58) 1.10 E 03
相邻CT 没有/是的 54/268 2.07 (1.07 - -3.99) 3.05 E 02 54/268 1.34 (0.65 - -2.76) 4.26 E 01
GSE21653
风险评分 低/高 126/126 2.83 (1.77 - -4.53) 1.32 E 05年 115/113 3.04 (1.80 - -5.11) 3.00 E 05年
年龄 ≤50 / > 50 92/160 1.14 (0.72 - -1.79) 5.83 E 01 81/147 1.31 (0.79 - -2.17) 2.95 E 01
T台 i ii / iii iv 178/66 1.70 (1.06 - -2.73) 2.75 E 02 164/64 1.20 (0.72 - -2.01) 4.87 E 01
N阶段 N0 / N1-3 116/133 1.54 (0.98 - -2.40) 5.84 E 02 105/123 1.61 (0.97 - -2.69) 6.81 E 02
ER状态 N / P 110/140 0.66 (0.43 - -1.02) 5.94 E 02 100/128 0.46 (0.17 - -1.26) 1.31 E 01
公共关系状态 N / P 124/126 0.84 (0.55 - -1.30) 4.35 E 01 114/114 1.43 (0.60 - -3.38) 4.15 E 01
HER2状态 N / P 207/26 1.59 (0.84 - -3.03) 1.58 E 01 204/24 1.09 (0.42 - -2.85) 8.60 E 01
三负 没有/是的 160/85 1.19 (0.74 - -1.89) 4.73 E 01 145/83 0.88 (0.26 - -2.90) 8.29 E 01
GSE42568
风险评分 低/高 52/52 3.68 (1.97 - -6.88) 4.55 E 05年 51/50 4.18 (2.13 - -8.17) 2.99 E 05年
年龄 ≤50 / > 50 27/77 0.66 (0.36 - -1.21) 1.79 E 01 25/76 0.72 (0.36 - -1.46) 3.63 E 01
N阶段 N0 / N1-3 45/59 4.35 (2.16 - -8.76) 3.88 E 05年 44/57 4.72 (2.24 - -9.95) 4.53 E 05年
ER状态 N / P 34/67 0.44 (0.24 - -0.79) 6.35 E 03 34/67 0.40 (0.21 - -0.77) 5.77 E 03

人力资源:风险比;置信区间:置信区间;呃:雌激素受体;公关:孕激素受体;HER2 / erbb2:表皮生长因子受体2;相邻CT:相邻的化疗;N:消极;P:积极的。

2.2。肿瘤样本

肿瘤样本和匹配相邻正常组织收集从20 BRCA病人在吉林大学第三医院(吉林、中国),2020年1月至2020年11月。口头同意了所有的病人。所有实验的伦理道德委员会批准吉林大学第三医院(吉林,中国)。所有样品都是匿名分析之前保证保护隐私。

2.3。差异表达TFs的识别

Limma包是用来识别差异表达TFs 114双的BRCA和正常控制样品( 17]。发现了差异表达TFs的截止条件<我nline-formula> 日志 2 褶皱 改变 足球俱乐部 > 1 和调整<我nline-formula> P 值≤0.05。所有<我nline-formula> P 值调整业务使用Benjamini和方法来控制错误发现率(罗斯福) 18]。火山地块创建基于R包中的ggplot2(3.2.0版本, https://cran.r-project.org/web/packages/ggplot2/index.html)[ 19]。层次聚类分析对微分TFs的表达进行了使用pheatmap包(1.0.12版本, https://cran.r-project.org/web/packages/pheatmap/index.html)[ 20.]。

2.4。预后特遣部队建设的签名

确定候选人预后TF,单变量回归分析进行评估的相关性差异表达TFs和DFS BRCA队列。只有TFs考克斯<我nline-formula> P < 0.05 和log-rank<我nline-formula> P < 0.05 被视为关联的生存。至少绝对收缩和选择操作符(套索)回归模型的敏感性和特异性,其具有相当大的优势及其系数可用于确定真正的药效监管组织( 21]。在这项研究中,套索回归分析来确定候选人是否有效的变量DFS,和那些<我nline-formula> 系数 = 0 被消除。基于签名的预后模型被用来计算每个病人使用的风险评分 (1) 风险 分数 = = 1 n E β , 在哪里<我nline-formula> n 是选择基因的数量,<我nline-formula> E 代表每个基因的表达,<我nline-formula> 和<我nline-formula> β 代表基因的系数<我nline-formula> 。值得注意的是,平均风险评分被定义为截止标准分类样本到高风险(<我nline-formula> 分数 中位数 )和低风险(<我nline-formula> 分数 中位数 )组。

为进一步验证的预测性能TF签名,时间的AUC接受者操作特征(ROC)曲线评估( 22]。定义点设置为1 - 2 -,3 -,4 - 5 -,10年时间,ROC曲线用来评估时间的风险评分结果的预测价值( 22]。更高的AUC表示更高的灵敏度。

2.5。基因集富集分析(GSEA) TF的签名和癌症的标志

癌症特征总结数据集的相关信息,从而减少变异和冗余。这提供了更精确和简洁投入GSEA分析( 23]。总共50标志基因集是从分子签名数据库(MSigDB,招募 http://software.broadinstitute.org/gsea/msigdb),导入到GSEA ( http://software.broadinstitute.org/gsea/index.jsp)。

2.6。实时定量PCR(存在)

验证mRNA表达水平之间的TFs TCGA肿瘤和正常乳腺组织样本数据集,这些TFs mRNA表达水平是衡量在收集新鲜冷冻肿瘤组织样本中存在。随后,存在实施从癌症和收集的RNA正常溶解产物。总之,总RNA分离通过FastPure细胞/组织总RNA隔离设备V2 (Vazyme生物技术有限公司,南京,中国)根据制造商的指导方针。接下来,RNA转录到cDNA使用HiScript II第一链cDNA合成装备。所有反应都是使用一个实现AceQ普遍SYBR qPCR大师混合(Vazyme生物技术有限公司,南京,中国)在QuantStudio™5实时PCR系统(ABI,生活技术,热费希尔科学、沃尔瑟姆,妈,美国)。表达式计算为δ的褶皱变化阈值周期(<我nline-formula> Δ Δ Ct 归一化后)一个内部参考。的引物对存在如表所示 S1

2.7。反褶积的渗透到肿瘤微环境中免疫细胞

细胞类型鉴定评估相对子集的RNA转录(CIBERSORT)是一种算法,用于描述复杂组织的细胞成分基于基因表达谱( 24]。在目前的研究中,CIBERSORT方法应用于研究免疫细胞浸润程度TCGA的肿瘤样本数据库。

2.8。统计分析

生存使用kaplan meier方法进行估计( 25]。DFS之间的时间被定义为诊断和疾病复发。统计生存时间差异评估使用log-rank测试( 26]。单变量和多变量回归分析进行使用Cox比例风险模型( 27]。<我nline-formula> P 0.05 被认为是具有统计学意义。

3所示。结果 3.1。预后TF签名的识别

在这项研究中,共有287个差异表达TFs被确定阈值以下的BRCA与正常对照组之间<我nline-formula> 日志 2 足球俱乐部 > 1 和<我nline-formula> P < 0.05 。287年TFs 85是调节和202是表达下调。分布的差异表达TFs描述使用火山的情节<我nline-formula> 日志 2 足球俱乐部 横坐标和吗<我nline-formula> - - - - - - 日志 10 罗斯福 是纵坐标(图 2(一个))。这些差异表达TFs的热图的可视化呈现在图 2 (b)和显示,BRCA(肿瘤)和正常对照组(正常)可能会被这些TFs(图 2 (b))。

识别的预后11-transcription因子(11-TF)签名。(一)火山差异表达TFs的情节。(b)差异表达TFs的热图。(c)无病生存期(DFS)曲线的12 survival-correlated TFs的高风险和低风险组。(d)套索回归分析的12 TF的候选人。(e)系数从套索算法获得12个TF的候选人。(f)预后11-TF签名的斯皮尔曼相关分析的预测乳腺癌的DFS (BRCA)。X:<我nline-formula> P > 0.05 。颜色酒吧:斯皮尔曼的系数(<我nline-formula> r )。

单变量回归分析的DFS导致提取12 survival-associated TFs候选人预后签名,包括E2F2 EGR3, EMX1, FOXD1, FOXJ1, NKX6-1, NR3C2, PAX7, STAT4, MYBL2 ZNF552, MSX1。如图 2 (c)、高表达E2F2 FOXD1 MSX1, MYBL2, NKX6-1, PAX7与贫穷DFS,虽然EGR3表达升高,EMX1, FOXJ1, NR3C2, STAT4和DFS ZNF552与改善。确定一个预后TF签名,套索回归分析进行的12 prognosis-associated TFs TCGA数据集(图 2 (d))。如图 2 (e)11 TFs选择构建预后TF签名,而MYBL2 0的回归系数是移除。一般信息的11个TFs表所示 2。斯皮尔曼的相关分析表明,表情识别TFs,组成特遣部队签名显示没有相关性,从而表明每个特遣部队(图的一个独立的预测能力 2 (f))。

细节的11个转录因子用于构造预后签名。

基因符号 基因稳定的ID 基因型 空空的 基因开始(bp) 基因(bp)结束
E2F2 ENSG00000007968 蛋白质编码 1 23506438 23531233
EGR3 ENSG00000179388 蛋白质编码 8 22687659 22693480
EMX1 ENSG00000135638 蛋白质编码 2 72916260 72936071
FOXD1 ENSG00000251493 蛋白质编码 5 73444827 73448777
FOXJ1 ENSG00000129654 蛋白质编码 17 76136333 76141245
NKX6-1 ENSG00000163623 蛋白质编码 4 84491987 84498450
MSX1 ENSG00000163132 蛋白质编码 4 4859665 4863936
NR3C2 ENSG00000151623 蛋白质编码 4 148078762 148444698
PAX7 ENSG00000009709 蛋白质编码 1 18630846 18748866
STAT4 ENSG00000138378 蛋白质编码 2 191029576 191151596
ZNF552 ENSG00000178935 蛋白质编码 19 57803841 57814913

缩写:装备=染色体。

3.2。基于11 TF签名预测的风险评分BRCA患者的生存

我们接下来的表达模式分析11 TFs TCGA的数据库。正常样本,表达E2F2相比,FOXJ1, EMX1, PAX7, NKX6-1 FOXD1, ZNF552显著提高(<我nline-formula> P < 0.05 )和MSX1的表达、STAT4 NR3C2, EGR3显著降低肿瘤样本(所有<我nline-formula> P < 0.05 )(图 3(一个))。与这些结果一致,存在证实这些TFs的信使rna表达水平显示类似的趋势TCGA表达式模式的数据库(图 3 (b))。随后,构建了风险评分模型,并用来预测BRCA患者的预后。根据风险评分中位数,BRCA患者分为低风险和高风险组(图 3 (c))。高危患者得分更容易复发患者相比,低风险的分数(图 3 (c))。分层集群性能表明,11预后TFs分类之间的低风险和高风险组(图 3 (d))。

TF表达水平、患者复发状态和集群的性能。(a)表达模式的11个TF预后签名BRCA患者之间(肿瘤)和正常对照组(正常),TCGA数据库。(b)的mRNA表达水平实验验证11 TF签名BRCA患者之间(肿瘤)和正常对照组(正常)。(c)复发BRCA患者在高风险和低风险组的状态。(d)的分层集群性能11预后TFs高风险和低风险组之间。

接下来,我们分析了生存预测能力的风险评分,TCGA数据库,这被认为是训练集,如图 4(一)低风险组中,病人已经大大延长DFS相比,病人在高危人群(<我nline-formula> 人力资源 = 4.44 ,<我nline-formula> 95年 % CI = 2.64 7.48 ,考克斯<我nline-formula> P = 2.13 e 08年 ,log-rank<我nline-formula> P = 8.77 e 10 )。接下来,预测风险评分的性能验证3地理数据和显示类似的预测价值(GSE20685:<我nline-formula> 人力资源 = 2.41 ,<我nline-formula> 95年 % CI = 1.58 3.68 ,考克斯<我nline-formula> P = 4.16 e 05年 ,log-rank<我nline-formula> P = 2.35 e 05年 ;GSE21653:<我nline-formula> 人力资源 = 2.83 ,<我nline-formula> 95年 % CI = 1.77 4.53 ,考克斯<我nline-formula> P = 1.32 e 05年 ,log-rank<我nline-formula> P = 5.45 e 06 ;和GSE42568:<我nline-formula> 人力资源 = 3.68 ,<我nline-formula> 95年 % CI = 1.97 6.88 ,考克斯<我nline-formula> P = 4.55 e 05年 ,log-rank<我nline-formula> P = 1.31 e 05年 )。

kaplan meier和接受者操作特征(ROC)曲线的BRCA患者分层风险评分。(a) kaplan meier和ROC曲线BRCA患者分层的风险评分在训练数据集(TCGA)和验证数据集(GSE20685, GSE21653, GSE42568)。(b) kaplan meier对先前发表的预测曲线和ROC曲线签名。

评估的敏感性和特异性11-TF签名,时间对DFS ROC曲线预测构造(图 4(一))。不包括1年期AUC值GSE20685 GSE21653,所有其他AUC值≥0.6。11-TF签名实现5年生存的AUC值为0.689,0.647,0.672,和0.721 TCGA GSE20685, GSE21653和GSE42568分别。这些结果表明,分数基于11-TF签名是有效的风险预测的DFS的BRCA病人和一个健壮的预测性能。

3.3。11-TF签名显示优越的预测性能相比以前公布的预后签名

在这项研究中,我们11-TF签名与4以前公布的预后签名面板:9-TF签名报道陈et al。 28周],12-lncRNA签名et al。 29日),13-gene后生签字保et al。 30.毛],12基因预后签字et al。 31日]。如图 4 (b)低风险组中,所有患者有一个长DFS相比,病人在高危人群(Chen等人:log-rank<我nline-formula> P = 1.24 e 04 ;包等。:log-rank<我nline-formula> P = 1.24 e 04 ;毛et al。: log-rank<我nline-formula> P = 3所示。4 e 03 ;和周et al。: log-rank<我nline-formula> P = 1.52 e 03 )。4以前公布的基因签名相比,我们获得一个更小的log-rank 11-TF签名<我nline-formula> P 价值。此外,我们11-TF签名1-10-year AUC值最高。因此,这些发现强调了我们11-TF签名的最佳性能。

3.4。风险评分是一个独立的预后因素

在这里,我们首先进行单变量回归分析探讨临床和病理因素之间的相关性和BRCA DFS的病人。如表所示 1风险评分,病理阶段,状态,和公关状态显著相关患者的DFS BRCA TCGA的数据库;风险评分,T台,N舞台,舞台,和相邻CT与GSE20685 DFS的患者;只有风险得分和T台的DFS BRCA GSE21653病人;风险评分,而N阶段,ER地位在GSE42568 DFS的病人有关。有趣的是,我们发现风险评分与整个BRCA DFS组相关。随后,多元回归分析探讨哪些因素是独立的预后因素。如表所示 1,只有风险评分在所有数据集是一个独立的预后因子。

我们进一步进行分层分析患者临床病理的变量,分子亚型和风险评分。如图 5(一个),BRCA病人低风险组相比更长的生存时间在所有分层的高危人群临床病理的变量组在每个数据集。具体来说,高,低风险组之间比较的曲线显示与log-rank显著差异<我nline-formula> P < 0.05 和考克斯<我nline-formula> P < 0.05 。此外,分子亚型的分层分析表明,高风险患者显示可怜的生存所有分子亚型(图的趋势 5 (b))。然而,没有观察到的差异为正常,DFS和基底亚型之间的高收入和低风险组。综上所述,这些数据表明,风险评分是一个独立的预后因素的预测的DFS BRCA病人。

kaplan meier和ROC曲线的BRCA患者临床病理因素分层,分子亚型和风险评分。(a) kaplan meier和ROC曲线BRCA患者分层的临床病理的因素和风险评分。(b) kaplan meier和ROC曲线的BRCA患者分层分子亚型和风险评分。

3.5。风险评分与免疫细胞浸润有关

使用CIBERSORT方法,我们估计22免疫细胞亚型的丰度调查之间的关系风险评分和肿瘤微环境中免疫细胞浸润。如图 6的分数M0和M2巨噬细胞,自然杀伤(NK)细胞,调节性T细胞,和记忆B细胞明显高于在高危人群,而本地B细胞,肥大细胞,单核细胞,内存和本地CD4休息+T细胞、CD8+T细胞,树突状细胞和巨噬细胞M1被明显高于低风险组。此外,两个免疫细胞比例最高的高危人群是M0和M2巨噬细胞。两种类型的免疫细胞,包括休息CD4记忆+T细胞、CD8+低风险患者的T细胞比例最高。

的作文分数高和低风险组之间的免疫细胞。

3.6。风险评分与癌症的特点相关联

调查如果11-TF签名与肿瘤生物学过程相关,基因集富集(GSEA)进行了分析。如图 7,共有6个特征基因集富集(<我nline-formula> P < 0.05 )。四个特征(包括MYC_TARGETS_V1 MYC_TARGETS_V2、糖酵解和DNA_REPAIR)与高风险有关分数,这表明这些生物过程的激活可能参与BRCA进展。相比之下,另外两个特点(ESTROGEN_RESPONSE_EARLY和ESTROGEN_RESPONSE_LATE)与低风险相关的成绩,表明他们的激活抑制肿瘤恶化,改善BRCA患者的生存。

(GSEA)的特征基因集富集分析。(一)四个特征丰富的高危人群。(b)两个特征丰富低风险组中。

4所示。讨论

因为BRCA异构的特点,传统的预测因素似乎并不是足够的BRCA患者的生存。因此,survival-specific预后有价值的生物标志物的发展是至关重要的。近年来,一些研究已经确定基因可能作为预后指标BRCA面板。例如,陈等人。 28]TCGA使用基因表达数据和临床数据和地理数据库识别9-TF签名,这可能发挥着重要在BRCA患者预后的作用。太阳et al。 32]分析了lncRNA表达谱的BRCA患者从GEO数据库来识别nine-lncRNA签名可以是一个有价值的预后生物标志物预测BRCA患者转移风险。在另一项研究中,BRCA患者的lncRNA表达谱分析了构建12-lncRNA签名来预测recurrence-free生存( 29日]。2019年,保et al。 30.)建立了一个13-gene后生签名通过结合mRNA表达和DNA甲基化数据集。12基因预后签名识别的基础上,结合独立BRCA数据库通过基因coexpression网络分析( 31日]。值得注意的是,所有这些研究集中在预后签名的识别;然而,他们并不认为免疫渗透分析和实施中存在的表达水平来验证这些签名使用冷冻组织样本。在我们的研究中,除了建立和验证预测模型的一个重要预测DFS 11-TF签名BRCA病人通过集成TF数据,基因表达数据集,也表现免疫渗透和临床信息,分析探讨相关的免疫细胞和TF签名。此外,我们验证了TFs使用中存在使用冷冻组织样本的表达水平。我们的研究结果显示,成功的风险评分是一个独立的预后因素分类肿瘤患者分为高风险和低风险组DFS的显著差异。因此,这些结果表明的重现性和可靠性对DFS的预测BRCA 11-TF签名。

11-TF预后模型是由E2F2 EGR3, EMX1, FOXD1, FOXJ1, NKX6-1, NR3C2, PAX7, STAT4 ZNF552, MSX1。在先前的研究中,陈等人报道,E2F家族共同或个别的TFs调节细胞增殖在癌症 33]。具体来说,BRCA e2f有预后价值,独立于临床参数( 34]。在先前的研究中,据报道,E2F2的激活与ER的耐抗雌激素治疗<我talic> α阳性的BRCA [ 35]。HER2-regulated E2F2表达式进一步影响cell-matrix附着力,与潜在后果转移殖民( 36]。此外,调节E2F2充当一个重要中间BRCA HER2-directed信号电路中( 37]。Malorni等人签名表明基因表达可能是建立基于与E2F1基因关联的调查和E2F2表达BRCA TCGA在数据库中( 38]。我们的结果是一致的结果在这些研究中,从而进一步验证在BRCA E2F2预后的关键功能。EGR3积极参与雌激素信号( 39]。井上等人证明EGR3中起关键作用的生理正常和恶性乳腺细胞,通过过量的感应基因( 40]。在这项研究中,GSEA表明雌激素反应的特点(ESTROGEN_RESPONSE_EARLY和ESTROGEN_RESPONSE_LATE)被浓缩在低风险病人组。我们推断,EGR3抑制TFs通过促进雌激素信号在BRCA进展的能力。狐狸家族的DNA结合蛋白质调节DNA转录和修复( 41]。有趣的是,在这项研究中,两个狐狸的家人(FOXD1和FOXJ1)被确定为在BRCA DFS预测患者预后签名。FOXD1函数作为癌基因在肺癌、乳腺癌、和脑癌 42, 43),是调节促进乳腺癌细胞增殖和药物抗性诱导G1年代过渡( 42]。FOXJ1 hypermethylated BRCA细胞系和临床组织样本,揭示其作为假定的肿瘤抑制基因( 44]。此外,DNA_REPAIR是癌症的一个重要标志,在高危人群丰富,这是按照福克斯家族的功能,从而验证我们的研究的一致性。

癌症是一种复杂的多级过程涉及遗传和表观遗传变化,导致致癌基因的激活信号和/或肿瘤抑制信号的失活( 45]。癌细胞获得一定数量的变化,促进肿瘤的生长和入侵 46]。GSEA MSigDB,最初开发,仍然是最大的和最受欢迎的存储库的基因数据集( 47]。GSEA数据库集中协调的微分表达式注释的基因或基因集和生产数据,可以更容易地解释有关的生物过程( 48]。我们利用GSEA从MSigDB选择DFS预测的重要标志和11-TF预后签名。MYC_TARGETS_V1, MYC_TARGETS_V2、糖酵解和DNA_REPAIR丰富的高危人群,而ESTROGEN_RESPONSE_EARLY和ESTROGEN_RESPONSE_LATE被显著富集低风险组。MYC促进细胞周期进程通过激活到Cdc25A, E2F1, E2F2 [ 49]。此外,减少扩散是由减少MYC enteroendocrine细胞/ E2F1活动( 50]。已经表明,MYC基因表达与GSEA高分数为MYC标志(MYC_TARGETS_V1和MYC_TARGETS_V2) ( 23]。2019年,玉等人表明,MYC在激进的扩散中扮演着一个关键的角色由BRCA细胞表型表现出表达ER<我talic> α突变( 51]。因此,针对MYC结合其他致癌途径提供了一种有前途的治疗策略BRCA [ 52]。在这项研究中,11 TFs被合并成一个面板,及其预后价值DFS BRCA成立。

最近,在几项研究,表明肿瘤浸润免疫细胞与癌症患者的预后相关( 53, 54]。进一步确定免疫浸润细胞之间的关系和TF签名,CIBERSORT分析揭示免疫浸润细胞的成分。在目前的研究中,高渗透的M2巨噬细胞与一个高风险的得分,而高渗透的休息CD4记忆+T细胞、CD8+T细胞是低风险评分。这些结果符合这些免疫细胞的生物学功能在癌症的发展过程中。在先前的研究中,据透露,M2巨噬细胞扮演了一个重要的角色在促进肿瘤的生长和转移 55]。此外,M2巨噬细胞已报告是与不利三阴性乳腺癌患者的预后( 56]。休息CD4记忆+T细胞、CD8+T细胞已被证明是与增加的总生存期(OS)和DFS的BRCA [ 57]。基于这些结果,我们相信这11-TF签名是可靠的风险评分预测BRCA患者的预后。

总之,我们已经确认和验证的预测预后11-TF签名DFS的BRCA患者。11-TF签名是一个独立的因素,可能作为一种补充临床病理因素的预后因素。

数据可用性

和/或使用的数据集分析在当前研究可从相应的作者以合理的要求。

同意

TCGA的知情同意和地理科目得到TCGA的地理协会,分别。

的利益冲突

我们国家没有利益冲突。

作者的贡献

CF负责收集文献和写作。JD和NL提出这项工作的想法,进行数据分析,实现实验验证和修订后的手稿。

补充材料

补充表S1:引物对qPCR。

DeSantis C。 J。 布莱恩 l Jemal 一个。 乳腺癌的统计数据,2013年 CA:临床医生的癌症杂志》上 2014年 64年 1 52 62年 10.3322 / caac.21203 2 - s2.0 - 84892795827 24114568 波里亚克 K。 在乳腺癌的异质性 临床研究杂志 2011年 121年 10 3786年 3788年 10.1172 / JCI60534 2 - s2.0 - 80053406583 21965334 我。 Q。 找到 K。 预测并为浸润性乳腺癌预后标记 病理学国际 2002年 52 3 186年 194年 10.1046 / j.1440-1827.2002.01335.x 2 - s2.0 - 0036961812 11972862 海斯 d F。 艾萨克斯 C。 斯登 V。 在乳腺癌预后因素:当前和转移的新预测 乳腺生物学和肿瘤杂志》上 2001年 6 4 375年 392年 10.1023 /:1014778713034 2 - s2.0 - 0035495969 12013528 Rodenhiser d . I。 安德鲁斯 j . D。 范登堡 t。 钱伯斯 答:F。 乳腺癌的进展和转移基因签名 乳腺癌研究 2011年 13 1 201年 10.1186 / bcr2791 2 - s2.0 - 84860389800 21345283 卡多佐 F。 凡不偏离 l . J。 博尔加特 J。 Slaets l Viale G。 Delaloge 年代。 Pierga j . Y。 大脑 E。 Causeret 年代。 DeLorenzi M。 格拉斯 a . M。 Golfinopoulos V。 Goulioti T。 诺克斯 年代。 马托斯 E。 Meulemans B。 Neijenhuis p。 Nitz U。 Passalacqua R。 Ravdin P。 卢比奥 i T。 Saghatchian M。 Smilde t·J。 Sotiriou C。 l Straehle C。 托马斯。 G。 汤普森 a . M。 van der胡芬 j . M。 Vuylsteke P。 伯纳德 R。 Tryfonidis K。 罗格斯大学 E。 Piccart M。 70 -基因签名作为援助在早期乳腺癌的治疗决定 新英格兰医学杂志》上 2016年 375年 8 717年 729年 10.1056 / nejmoa1602253 2 - s2.0 - 84984697640 兰伯特 美国一个。 Jolma 一个。 Campitelli l F。 达斯 p K。 Y。 Albu M。 程ydF4y2Ba X。 Taipale J。 休斯 t·R。 Weirauch m . T。 人类转录因子 细胞 2018年 172年 4 650年 665年 10.1016 / j.cell.2018.01.029 2 - s2.0 - 85041527619 29425488 阿伦特 D。 j . M。 贝克 c . v . H。 伯格曼 一个。 Cepko C。 欧文 d . H。 Pavlicev M。 Schlosser G。 威德 年代。 Laubichler m D。 瓦格纳 g . P。 细胞的起源和演化 自然遗传学评论 2016年 17 12 744年 757年 10.1038 / nrg.2016.127 2 - s2.0 - 84994079481 布沙尔 C。 “斯太尔 P。 Eilers M。 通过Myc基因控制的细胞增殖 细胞生物学的趋势 1998年 8 5 202年 206年 10.1016 / s0962 - 8924 (98) 01251 - 3 2 - s2.0 - 0032080864 9695840 Shaulian E。 卡琳 M。 生与死AP-1作为管理者的细胞 自然细胞生物学 2002年 4 5 E131 E136 10.1038 / ncb0502-e131 2 - s2.0 - 0036098552 11988758 C。 l . M。 阿尔维斯 P。 格斯坦 M。 系统的识别与病人相关的转录因子在癌症生存 BMC基因组学 2009年 10 1 225年 10.1186 / 1471-2164-10-225 2 - s2.0 - 66349112272 克莱因 m . S。 工艺 B。 Swatloski T。 高盛 M。 年代。 Haussler D。 J。 探索TCGA pan-cancer UCSC的癌症基因组浏览器的数据 科学报告 2013年 3 1,第2652条 10.1038 / srep02652 2 - s2.0 - 84885137302 花王 k·J。 k . M。 h . C。 a . T。 芯片乳腺癌分子亚型的相关性及临床结果:影响治疗的优化 BMC癌症 2011年 11 1 143年 10.1186 / 1471-2407-11-143 2 - s2.0 - 79955121426 Sabatier R。 Finetti P。 赛尔维拉 N。 Lambaudie E。 Esterni B。 Mamessier E。 Tallet 一个。 Chabannon C。 额外的 j . M。 Jacquemier J。 来吧 P。 伯恩鲍姆 D。 Bertucci F。 一个基因表达特征识别两种基底乳腺癌预后子组 乳腺癌研究和治疗 2011年 126年 2 407年 420年 10.1007 / s10549 - 010 - 0897 - 9 2 - s2.0 - 79958728723 20490655 Sabatier R。 Finetti P。 阿德莱德 J。 Guille 一个。 Borg j . P。 Chaffanet M。 车道 l 伯恩鲍姆 D。 Bertucci F。 ECRG4下调,候选肿瘤抑制基因,在人类乳腺癌 《公共科学图书馆•综合》 2011年 6 11条e27656 10.1371 / journal.pone.0027656 2 - s2.0 - 81155131661 22110708 克拉克 C。 马登 美国F。 Doolan P。 Aherne s T。 乔伊斯 H。 奥德利 l 加拉格尔 w·M。 轩尼诗 b . T。 莫里亚蒂 M。 皇冠 J。 肯尼迪 年代。 Clynes M。 乳腺癌相关转录网络生存:大规模coexpression分析 致癌作用 2013年 34 10 2300年 2308年 10.1093 / carcin / bgt208 2 - s2.0 - 84885131104 23740839 里奇 m E。 Phipson B。 D。 Y。 法律 c·W。 W。 史密斯 g·K。 limma权力分析RNA-sequencing与微阵列研究微分表达式 核酸的研究 2015年 43 7 e47 e47 10.1093 / nar / gkv007 2 - s2.0 - 84926507971 25605792 Brohawn d·G。 O ' brien l . C。 班尼特 j . P。 Jr。 RNAseq分析识别肿瘤坏死factor-mediated ALS脊髓炎症的主要异常 《公共科学图书馆•综合》 2016年 11 8条e0160520 10.1371 / journal.pone.0160520 2 - s2.0 - 84983071394 27487029 韦翰 H。 ggplot2:优雅的图形进行数据分析 2009年 斯普林格出版公司,注册 10.1007 / 978-0-387-98141-3 Kolde R。 pheatmap:漂亮的热图 2015年 Y。 Y。 W。 M。 C。 一个套索回归模型microRNA-target监管网络的建设 生物信息学 2011年 27 17 2406年 2413年 10.1093 /生物信息学/ btr410 2 - s2.0 - 80051931399 21743061 Heagerty p . J。 广告 T。 佩佩 m . S。 对审查的生存时间ROC曲线数据和诊断标记 生物识别技术 2000年 56 2 337年 344年 10.1111 / j.0006 - 341 x.2000.00337.x 2 - s2.0 - 0033936550 10877287 Liberzon 一个。 birge C。 Thorvaldsdottir H。 甘地 M。 Mesirov j . P。 Tamayo P。 的分子签名数据库特征基因集集合 电池系统 2015年 1 6 417年 425年 10.1016 / j.cels.2015.12.004 2 - s2.0 - 84955328286 26771021 纽曼 a . M。 c . L。 绿色 m·R。 叠成 a·J。 W。 Y。 黄平君 c, D。 Diehn M。 Alizadeh 答:一个。 健壮的枚举的细胞组织表达谱的子集 自然方法 2015年 12 5 453年 457年 10.1038 / nmeth.3337 2 - s2.0 - 84928927858 25822800 Y。 罗伊斯顿 P。 从发表kaplan meier曲线重建比较数据 占据杂志 2017年 17 4 786年 802年 10.1177 / 1536867 x1801700402 2 - s2.0 - 85041535161 Koletsi D。 Pandis N。 生存分析,第2部分:kaplan meier方法和生存率较 美国口腔正畸学杂志》和Dentofacial矫形手术 2017年 152年 4 569年 571年 10.1016 / j.ajodo.2017.07.008 2 - s2.0 - 85032182021 琼斯 m P。 克罗利 J。 生存分析的非参数测试的一般类 生物识别技术 1989年 45 1 157年 170年 10.2307 / 2532042 2 - s2.0 - 0024523040 2655728 程ydF4y2Ba H。 X。 M。 M。 l T。 转录因子分析来预测recurrence-free生存在乳腺癌:开发和验证的诺模图优化临床管理 遗传学前沿 2020年 11 10.3389 / fgene.2020.00333 M。 l W。 太阳 Y。 Z。 H。 l 太阳 J。 发现潜在的预后长非编码RNA生物标志物预测乳腺癌患者的肿瘤复发的风险 科学报告 2016年 6 1,第31038条 10.1038 / srep31038 2 - s2.0 - 84981206974 X。 Anastasov N。 Y。 Rosemann M。 小说表观遗传特征对整体的生存预测乳腺癌患者 转化医学杂志》 2019年 17 1 10.1186 / s12967 - 019 - 2126 - 6 x Y。 m·J。 J。 C。 法律 s M。 Snijders a . M。 全基因组屏幕识别小说对乳腺癌的生存预后的签名 Oncotarget 2017年 8 8 14003年 14016年 10.18632 / oncotarget.14776 2 - s2.0 - 85013377381 28122328 太阳 J。 程ydF4y2Ba X。 Z。 M。 H。 X。 l M。 一个潜在的预后长非编码RNA的签名来预测metastasis-free乳腺癌患者的生存 科学报告 2015年 5 1 16553年 10.1038 / srep16553 2 - s2.0 - 84946905311 程ydF4y2Ba h . Z。 s Y。 里昂 G。 新兴的角色e2f在癌症:退出细胞周期控制 自然评论癌症 2009年 9 11 785年 797年 10.1038 / nrc2696 2 - s2.0 - 70350518233 19851314 约翰逊 J。 Thijssen B。 麦克德莫特 U。 加内特 M。 欧洲鞋号 l·f·A。 伯纳德 R。 针对乳腺癌的RB-E2F通路 致癌基因 2016年 35 37 4829年 4835年 10.1038 / onc.2016.32 2 - s2.0 - 84959218324 26923330 米勒 t·W。 大多数 j . M。 福克斯 e . M。 Ghazoui Z。 Dunbier 一个。 安德森 H。 Dowsett M。 一个。 史密斯 r。 玛丽亚 s M。 曼宁 h . C。 Gonzalez-Angulo a . M。 米尔斯 g . B。 海厄姆 C。 Chanthaphaychith 年代。 库巴地毯 m·G。 米勒 w·R。 Shyr Y。 Arteaga c . L。 呃<我talic> α端依赖E2F转录调节抗雌激素可以剥夺人类乳腺癌 癌症的发现 2011年 1 4 338年 351年 2159 - 8290. - 10.1158 / cd - 11 - 0101 2 - s2.0 - 82055187255 22049316 Shibue T。 温伯格 r。 整合素<我talic> β1-focal黏附激酶信号指导肺部转移性癌细胞的扩散传播 美国国家科学院院刊》上的美利坚合众国 2009年 106年 25 10290年 10295年 10.1073 / pnas.0904227106 2 - s2.0 - 67649849871 Bollig-Fischer 一个。 l Mitrea C。 J。 克鲁格 一个。 Manca V。 Drăghici 年代。 建模时间HER2致癌基因转录的影响,发现在乳腺癌cell-matrix E2F2粘附作用 生物信息学 2014年 30. 21 3036年 3043年 10.1093 /生物信息学/ btu400 2 - s2.0 - 84911436760 25028721 Malorni l 广场 年代。 Y。 Guarducci C。 Bonechi M。 Biagioni C。 哈特 c, D。 Verardo R。 利奥 答:D。 阶段 我。 视网膜母细胞瘤的基因表达特征功能丧失是抵抗的预测生物标志物palbociclib在乳腺癌细胞系和ER阳性早期乳腺癌患者的预后 Oncotarget 2016年 7 42 68012年 68022年 10.18632 / oncotarget.12010 2 - s2.0 - 85015019288 27634906 井上 一个。 Omoto Y。 山口那津男 Y。 Kiyama R。 哈亚希 我美国。 转录因子EGR3参与estrogen-signaling通路在乳腺癌细胞 分子内分泌学杂志》 2004年 32 3 649年 661年 10.1677 / jme.0.0320649 2 - s2.0 - 3042717657 15171706 井上 一个。 吉田 N。 Omoto Y。 Oguchi 年代。 Yamori T。 Kiyama R。 哈亚希 年代。 发展为过量的基因表达谱芯片技术 分子内分泌学杂志》 2002年 29日 2 175年 192年 10.1677 / jme.0.0290175 2 - s2.0 - 0036803020 12370120 Katoh M。 Igarashi) M。 福田 H。 Nakagama H。 Katoh M。 癌症遗传学和基因组学人类福克斯家族的基因 癌症的信 2013年 328年 2 198年 206年 10.1016 / j.canlet.2012.09.017 2 - s2.0 - 84870376109 23022474 y F。 j . Y。 H。 k . S。 H。 z G。 x J。 FOXD1促进乳腺癌扩散,针对p27化疗耐药性 生物化学和生物物理研究通信 2015年 456年 1 232年 237年 10.1016 / j.bbrc.2014.11.064 2 - s2.0 - 84919698298 25462566 y F。 T。 x Y。 c . X。 l j . Y。 H . H。 z Q。 沉默Forkhead盒D1抑制胶质瘤细胞增殖和迁移 肿瘤的报道 2017年 37 2 1196年 1202年 10.3892 / or.2017.5344 2 - s2.0 - 85011933915 28075458 Siu m . k . Y。 e . s . Y。 香港 d . s . H。 常ydF4y2Ba h . Y。 l o·g·W。 e·w·F。 常ydF4y2Ba k . k . L。 h . y S。 x F。 a . N。 干细胞转录因子通过失调NANOG控制细胞迁移和入侵上皮和FoxJ1导致不良临床结果卵巢癌 致癌基因 2012年 32 30. 3500年 3509年 10.1038 / onc.2012.363 2 - s2.0 - 84880923339 22945654 特里帕西 b K。 Anderman m F。 朵罗修 j . H。 罗伊 d·R。 摘要4513年。针对DLC1 Rho-GAP肿瘤抑制蛋白治疗实体肿瘤 诉讼:AACR年会2019年 2019年 亚特兰大,乔治亚州 Sanchez-Tillo E。 Y。 德·巴里奥斯 O。 大雨如注 l Fanlo l Cuatrecasas M。 亲爱的 d S。 迪安 d . C。 中文版 一个。 Postigo 一个。 EMT-activating转录因子在癌症:除了EMT和肿瘤侵袭性 细胞和分子生命科学 2012年 69年 20. 3429年 3456年 10.1007 / s00018 - 012 - 1122 - 2 2 - s2.0 - 84867578091 22945800 Liberzon 一个。 萨勃拉曼尼亚 一个。 Pinchback R。 Thorvaldsdottir H。 Tamayo P。 Mesirov j . P。 3.0分子特征数据库(MSigDB) 生物信息学 2011年 27 12 1739年 1740年 10.1093 /生物信息学/ btr260 2 - s2.0 - 79958135926 21546393 诉K。 P。 米勒 K。 Bunkenborg J。 帝国 M。 Hjerrild M。 Delmonte T。 维伦纽夫 一个。 Sladek R。 F。 米切尔 g。 莫林 C。 曼恩 M。 哈德逊 t·J。 罗宾逊 B。 Rioux j . D。 着陆器 大肠。 识别的基因导致人类细胞色素c氧化酶缺乏整合基因组学 美国国家科学院院刊》上的美利坚合众国 2003年 One hundred. 2 605年 610年 10.1073 / pnas.242716699 2 - s2.0 - 0037458031 12529507 迈耶 N。 l Z。 反思与MYC 25年 自然评论癌症 2008年 8 12 976年 990年 10.1038 / nrc2231 2 - s2.0 - 56749184298 19029958 Basak O。 J。 Wiebrands K。 Seno H。 范Oudenaarden 一个。 聪明的 H。 诱导静止Lgr5就+干细胞在肠道瀑样使分化hormone-producing enteroendocrine细胞 细胞干细胞 2017年 20. 2 177年 e4 190. 10.1016 / j.stem.2016.11.001 2 - s2.0 - 85008153076 27939219 l l C。 Duraki D。 j·E。 R。 赫尔弗里希 w·G。 纳尔逊 e·R。 公园 b . H。 夏皮罗 d . J。 Estrogen-independent Myc基因超表达赋予内分泌治疗抵抗乳腺癌细胞表达ER<我talic> αY537S和ER<我talic> αD538G突变 癌症的信 2019年 442年 373年 382年 10.1016 / j.canlet.2018.10.041 2 - s2.0 - 85057072193 30419347 J。 程ydF4y2Ba Y。 Olopade o . I。 MYC和乳腺癌 基因与癌症 2010年 1 6 629年 640年 10.1177 / 1947601910378691 21779462 浓密冰 r D。 Sotiriou C。 Piccart-Gebhart m·J。 Haanen j·b·a·G。 van Vugt M a t M。 德弗里斯 E . g . E。 施罗德 c·P。 Fehrmann r . s . N。 肿瘤浸润免疫细胞组成和功能的相关性疾病结果乳腺癌 美国国家癌症研究所杂志》上 2016年 109年,文章djw192 1 10.1093 / jnci / djw192 2 - s2.0 - 85014873742 Y。 Q。 Y。 T。 J。 l G。 H。 谅解备忘录 T。 年代。 ImmunoScore签名:胃癌的预后和预测工具 年报的手术 2016年 1 M。 J。 C。 J。 J。 ICAM-1抑制肿瘤转移通过抑制巨噬细胞通过封锁efferocytosis M2极化 细胞死亡和疾病 2015年 6 6条e1780 10.1038 / cddis.2015.144 2 - s2.0 - 84989967720 X。 R。 T。 Y。 W。 综合分析单细胞RNA-seq和散装RNA-seq拆解肿瘤异质性+ M2-like TNBC肿瘤相关巨噬细胞浸润和侵略性 癌症免疫学和免疫疗法 2021年 70年 1 189年 202年 10.1007 / s00262 - 020 - 02669 - 7 研究所。 Z.-Q。 黄永发。 G.-M。 Y。 Y。 J。 欧阳 Y。 Z。 肿瘤浸润免疫细胞在乳腺癌的临床意义 癌症杂志》 2019年 10 24 6175年 6184年 10.7150 / jca.35901 31762828