跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
Mol Cancer研究。作者手稿;2019年12月1日PMC提供。
以最终编辑形式发布为:
2019年3月1日在线发布。 数字对象标识:10.1158/1541-7786.MCR-18-1054
预防性维修识别码:PMC6548593型
NIHMSID公司:美国国家卫生研究院1523157
PMID:30824526

胆固醇硫转移酶SULT2B1b调节去势耐药前列腺癌对死亡受体配体TNF-α的敏感性

关联数据

补充材料

摘要

胆固醇硫转移酶SULT2B1b已被证明可以调节雄激素受体活性和细胞生长特性。然而,SULT2B1b改变前列腺癌细胞内这些特性的机制尚未描述。此外,尚不清楚SULT2B1b在前列腺癌细胞中表达的具体优势。在这些研究中,进行了单细胞信使核糖核酸测序(scRNA-seq),以比较SULT2B1b敲除(KD)的转录组控制KD LNCaP细胞。超过2000个差异表达(DE)基因被鉴定,以及许多典型通路的改变,包括死亡受体信号通路。本文的研究表明,SULT2B1b KD增加前列腺癌细胞中肿瘤坏死因子α(TNF)的表达,并以TNF依赖的方式导致NF-κB活化。更重要的是,SULT2B1b KD显著增强TNF-敏感LNCaP细胞和TNF-耐药C4-2细胞中TNF-介导的凋亡。LNCaP细胞中SULT2B1b的过度表达也降低了对TNF介导的细胞死亡的敏感性,这表明SULT2B1调节了决定前列腺癌细胞TNF敏感性能力的途径。对人类前列腺癌患者数据集的探测进一步支持了这项工作,提供了SULT2B1b表达与TNF相关基因呈负相关的证据,包括TNF公司,CD40长度,FADD公司、和NFKB1型总之,这些数据提供了证据,证明前列腺癌细胞中SULT2B1b的表达增强了对TNF的抵抗力,并可能提供生长优势。此外,靶向SULT2B1b可能会增强晚期前列腺癌对TNF治疗的疗效。

关键词:SULT2B1b、前列腺癌、TNFα、死亡受体信号、scRNA-seq

简介

前列腺癌是美国男性非皮肤癌死亡的第二大原因,也是最常见的诊断。(1)几十年来,死亡受体靶向疗法,特别是Fas、肿瘤坏死因子α(TNF)和TNF相关的凋亡诱导配体(TRAIL),已被研究为该疾病的潜在疗法。(2——4)TNF在40多年前被发现,并被证明可以诱导肿瘤细胞坏死。(5)在初步观察之后,TNF被定义为在前列腺癌细胞中具有看似矛盾的作用。TNF已被证明在一些前列腺癌细胞中诱导凋亡,如LNCaP细胞,而其他细胞,如去势抵抗前列腺癌(CRPC)细胞系C4-2,对TNF介导的细胞死亡具有抵抗力。(6,7)前列腺癌中的TNF研究为前列腺癌细胞中通过TNF-受体(TNFR1和TNFR2)的TNF信号转导产生的促肿瘤和抗肿瘤反应提供了证据。(8——14)然而,细胞内环境因素对这种自相矛盾的观察结果的作用机制尚待阐明。前列腺癌中细胞内环境可能影响TNF反应性的机制是一个尚未研究的领域,可能会影响治疗方法,尤其是对TNF介导的凋亡耐药的前列腺癌细胞。

前列腺癌细胞也表现出胆固醇的失调。(15,16)我们最近报道了胆固醇硫转移酶SULT2B1b的RNA干扰(RNAi)介导的敲除(KD)在多种前列腺癌细胞系中诱导凋亡,包括CRPC模型C4-2细胞。(17)为了更好地了解SULT2B1b诱导前列腺癌细胞死亡的机制,我们进行了单细胞mRNA测序(scRNA-seq)分析,比较了SULT2B1 bKD和非靶向(对照)KD LNCaP细胞。选择单细胞方法进行此分析,因为它可以在细胞水平上测定siRNA KD的质量,以及鉴定SULT2B1b调节通路的特定基因表达的相关性。使用LNCaP细胞作为模型限制了这些实验中的异质性,但也发现了SULT2B1b介导调控的许多可能途径。了解SULT2B1b在前列腺癌细胞中的作用以及该酶如何影响细胞生长特性,将有助于更深入地了解前列腺癌细胞的生物学特性,并有助于确定CRPC的替代治疗方案。在此,我们提供证据表明,SULT2B1b调节TNF介导的细胞信号传导,并且通过敲除靶向SULT2B1可增强TNF诱导的LNCaP细胞和TNF-耐药的CRPC C4-2细胞的细胞死亡。

材料和方法

单元格行

LNCaP细胞购自美国组织培养收藏中心(ATCC;Manassas,VA),并保存在与ATCC推荐的培养基条件相同的培养基中。C4-2细胞购自MD Anderson(德克萨斯州休斯顿),保存在补充有10%胎牛血清(FBS)和1%链霉素/青霉素的T培养基(Invitrogen)中。如前所述,LNCaP-SULT2B1b细胞通过连续慢病毒转导产生,通过添加强力霉素诱导SULT2B1b过度表达。(17)在前列腺癌基金会的慷慨支持下,所有使用LNCaP和C4-2细胞的实验都是在DDC Medical(俄亥俄州费尔菲尔德)之前鉴定的细胞的20代内完成的。(17)

抗体和试剂

针对以下蛋白质的主要抗体用于这些研究:TNFα(Thermo,P300A)和TRAIL/TNFSF10(R&D Systems,MAB375)用于中和;β-actin(Cell Signaling,3700)、phospho-IkB(Cell信号,9246)、IkBα(Cell信令,4814)、TRADD(BD,610572)和SULT2B1(Abcam,ab88085)用于蛋白质印迹;以及用于流式细胞术的TNFR1(Abcam,ab19139)和TNFR2(Abcam,ab109322)。次级抗体包括用于免疫印迹的山羊抗鼠IgG-HRP(Santa Cruz,sc-2005)和用于流式细胞术的驴抗兔IgG-AlexaFluor647(Biolegend,406414)。

使用控制/非靶向siRNA复合物(Dharmacon)和siRNA靶向完成RNAi研究SULT2B1系列(IDT,HSC.RNAI.N004605.12.2),TNF公司(IDT,hs.Ri.TNF.13.2),TL1(IDT、hs.Ri.TNFSF15.13.3),或DAXX公司(IDT,hs.Ri.DAXX.13.1)根据制造商的说明,使用Lipofectamine RNAiMax(ThermoFisher)。如前所述完成RNA分离和cDNA合成(17),源自EZNA Total RNA Kit I(Omega Bio-tek)。用于这些研究的PrimeTime®qRT-PCR基因探针(IDT)包括:AR公司(Hs.PT.56a.38770693),KLK3号机组(PSA)(Hs.PT.56a.38546086),SULT2B1系列(Hs.PT.56a.25562421.g),TNF公司(Hs.PT.58.45380900),以及SULT2B1b型(Hs.PT.58.22608626)。

此外,用于细胞裂解物的重组人TNFα(Peprotech,300–01A)、人TRADD细胞ELISA试剂盒(Abnova,KA3564)和人TNF-αELISA试剂瓶(Sigma-Aldrich,RAB1089)也用于这些研究。按照制造商的说明使用TRADD ELISA试剂盒,方法是在siRNA转染72小时后固定粘附细胞,并将其归一化为结晶紫吸光度,而TNFαELISA则分别使用在LNCaP和C4-2细胞中siRNA转导72小时或96小时后样品中的细胞裂解物进行。为了评估细胞活力,根据制造商的说明,在这些研究中使用了细胞计数试剂盒-8(Dojindo Molecular Technologies,Inc.)。

萤光素酶检测

通过转染pNF-κB-荧光素酶报告质粒(Stratagene)和雷尼利亚使用FuGENE HD转染试剂(Promega,E2311)的荧光素酶质粒(pRL-TK),然后使用双荧光素素酶报告试剂盒(Promeca,E1910)评估荧光素酶活性。(18)相对荧光素酶活性(RLU=萤火虫/雷尼利亚)NF-κB报告子结构的平均值显示为至少三个独立实验(一式三份)的折叠变化(在没有siRNA转染的样品上)+/-SEM的平均值。有关更多详细信息,请参阅补充方法.

流式细胞术

对于细胞周期分析,根据制造商的说明,用碘化丙啶(Biolegend#421301)对固定细胞进行染色,并在FACSCanto II(BD Biosciences)上进行分析。使用FlowJo软件版本9(TreeStar Inc.)对DNA含量进行进一步分析。对于蛋白质表达分析,在BD Fortessa LSR(BD Biosciences)上对细胞进行分析,并使用FlowJo软件将实验样品的平均荧光强度(MFI)与二级对照样品进行比较。

细胞制备和装入C1系统

siRNA转染48小时后,用胰蛋白酶采集LNCaP细胞,并用僵尸紫活性染料(Biolegend,423114)染色,然后在FACS ARIA III(BD Biosciences)上分选活细胞。细胞在17–25µm IFC板(Fluidigm,100–5761)中加载并染色,然后排除含有双重细胞或死细胞的捕获位点。对对照组和SULT2B1b KD细胞连续几天完成该程序,并对生物复制重复三次。总共使用209个对照siRNA和190个SULT2B1b siRNA处理的细胞进行分析。有关C1加载和后续库准备的更多详细信息,请参阅补充方法.

cDNA采集、定量、文库制备和测序

使用带有dsDNA HS分析试剂盒(Life Technologies,问题32851)并归一化为0.1–0.3 ng/µL的浓度,用于文库制备。使用Nextera®XT DNA样品制备试剂盒和索引试剂盒(Illumina,FC-131-1096和FC-131-1002)制备文库,然后根据Fluidgm C1 mRNA测序协议进行量化。在快速运行模式下,使用HiSeq 2500系统(Illumina)对未经训练的2×100bp读数进行测序。进行了三次单独的测序运行,每一次都在一条通道中使用共用的对照KD细胞库,在另一条通道上使用共用的SULT2B1b KD细胞图书馆。添加了cBot Duo(Illumina),以便在单个快速流动细胞中添加不同的样品池。

序列修剪和对齐

通过FastX-Toolkit版本0.0.13.2对序列进行质量调整并删除适配器(19). 使用了30的FastX修剪分数和50的修剪长度。最大长度设置为151个底座。FastQC版本0.11.2(20)运行以检查质量修整和适配器移除前后的数据质量。

使用Tophat2将读数与合集基因组参考联盟人类构建38(GRCh38.p3)一致(21,22)Tophat2 2.1.0版本使用默认设置运行,但有两个例外:一个不匹配是允许的,并且由于测序库的非链特异性,库类型设置为“fr unstranded”。HTSeq v.0.6.1(23)在“intersection-nonempty”模式下用于计算映射到特征的读取数;本分析中使用了Biopython v.2.7.3。

统计分析和路径分析

的统计分析在体外学生使用的研究t吨-使用GraphPad Prism 5软件进行方差测试或分析(ANOVA),根据需要进行或不进行多重比较测试。统计显著值用数字表示,如*p<0.05、**p<0.01、**p<0.001或****p<0.0001(无显著性以n.s.表示)。所有scRNA-seq数据都可以在GEO中获得,并带有登录号GSE117410共对36135个基因进行了测序,其中许多基因表现出低表达水平。移除所有样本中平均计数小于5的基因,然后使用Bioconder R软件包测试差异表达边缘R第3.2.2节,。使用Benjamini-Hochberg程序控制5%的错误发现率(FDR),获得2029个差异表达(DE)基因。DE基因、FDR、log(fold-change)和log(counts per million)被上传到Ingenuity Pathway Analysis(IPA)软件(Qiagen),并进行了典型通路分析和上游调节器分析。根据输入的DE基因在IPA中预测上游调控因子,并使用单侧Fisher精确检验确定p值。

人类前列腺癌数据库相关性

从cBioPortal数据库中检索到Robinson等人的RNA-seq数据,包括20例未经治疗的骨髓源性转移性前列腺癌样本和16例经紫杉烷和阿比特龙或苯扎鲁胺治疗的淋巴结源性CRPC样本。(24)数据通过log(RPKM+1)归一化。基因共表达相关性SULT2B1系列使用Pearson相关系数计算55个肿瘤坏死因子和受体相关基因,并使用10000轮随机模拟的排列测试评估其统计显著性。

结果

scRNA-seq分析成功鉴定了改变的通路和DE基因。

在这些研究中,进行了scRNA-seq以确定SULT2B1b KD中显著改变的基因和通路控制KD前列腺癌细胞。我们之前的研究表明,SULT2B1b KD在72小时内诱导LNCaP细胞凋亡。由于scRNA-seq需要活细胞,因此在非靶向或SULT2B1 siRNA(分别为对照KD或SULT2D1b KD)转染48小时后采集LNCaP,然后在Fluidigm C1单细胞自动制备系统上进行单细胞分离之前进行活细胞分选(补充图1A). 活细胞分选不影响SULT2B1b KD的效率(补充图1B). 分别为对照组或SULT2B1b KD完成三个独立实验(第1-3批),并将测序结果数据汇总用于质量控制和分析,分别得到209个对照KD和190个SULT2Bb KD细胞(图1A). 序列读取被确定为高质量,并且在分析过程中确定了最小的批量影响(补充图1C-D).

保存图片、插图等的外部文件。对象名为nihms-1523157-f0001.jpg
SULT2B1b KD的单链核糖核酸序列对照KD细胞证实AR活性降低。

(A) 测序的单个细胞数量概述。(B) 多维标度(MDS)图强调了单个细胞组之间因治疗而产生的差异。“C”表示控制KD,“KD”表示SULT2B1b KD。(C) 数量SULT2B1系列对于对照KD(对照组)或SULT2B1b KD(敲除组)组,指示单细胞水平上SULT2B1 b表达的读数。右上角表示每组中读取计数为零的单元格的百分比。(D) 左:小提琴图,表示scRNA-seq组AR的表达(读取计数);右:qRT-PCR表达AR公司从转染后48小时的大量样品中提取。(E) 显示AR靶基因的小提琴图。每个基因的表达(读取计数)用调整后的p值表示。

多维标度图显示了源于Control和SULT2B1b KD条件的细胞之间的分离程度(图1B). 采用负二项分布进行统计分析,然后采用Benjamini-Hochberg程序,在对照组和SULT2B1b KD组之间鉴定出2029个DE基因。我们的分析确定,与对照KD处理的细胞相比,SULT2B1b KD处理细胞中的SULT2B1 b读取计数有所减少,但并不显著(0个读取计数的数量较高;图1C)AR活性显著降低,如AR靶基因表达减少(图1D-E). 需要注意的是,0读计数表示转录水平低于此检测的检测限,而不是没有表达。值得注意的是,与scRNA-seq数据集中的对照KD细胞相比,SULT2B1b KD细胞在48小时时AR基因的表达保持不变,而正如我们之前报道的那样,在SULT2B1 KD细胞72小时后,通过qRT-PCR获得的大量AR表达减少(图1D,(17)).

使用Ingenuity Pathway Analysis(IPA)对scRNA-seq数据进行分析,确定了许多显著改变的典型路径(补充图2)和推测的DE基因上游调控因子。尽管最显著改变的途径显示了应激的细胞内环境(细胞凋亡在(17))死亡受体信号通路的上调被认为是SULT2B1b KD介导的细胞凋亡的一种令人惊讶的可能机制(补充图2). 死亡受体激活在前列腺癌中具有已知的意义,死亡受体已被证明定位于富含胆固醇的结构域,但死亡受体激活尚未被证明与胆固醇修饰相关的基因的敲除激活,例如磺化作用。(25,26)由于该通路与细胞死亡直接相关,我们旨在验证SULT2B1b是否参与死亡受体信号通路的调节,并确定SULT2B1 bKD可能导致死亡受体信号激活的机制。

NF-κB活化由SULT2B1b KD处理的LNCaP细胞内的TNF表达诱导。

作为死亡受体信号的贡献者和参与者,在scRNA-seq数据集中研究了TNF和NF-κB。TNF和NF-κB复合物均被预测激活(z-score=4.552,p=1.68E-08和z-score=4.383,p=2.75E-06)许多DE基因的上游调控因子。此外,四个TNF诱导的NF-κB调节基因在SULT2B1b KD细胞中显著上调,包括MnSOD、A20、cIAP2和Bcl-XL(图2). 利用LNCaP和C4-2细胞进行了进一步的研究,以验证SULT2B1b对该途径的调节。LNCaP细胞中SULT2B1b KD诱导TNF公司NF-κB的表达和激活,支持scRNA-seq数据(图3A-D). 在患有SULT2B1b KD的LNCaP细胞中,TNF蛋白的表达没有显著升高,但与对照KD细胞相比,其表达有上升趋势(图3B). 在携带SULT2B1b KD的C4-2细胞中,TNF表达和NF-κB活性均未显著增加,但SULT2B1 B KD确实激活了IκB的磷酸化(图3A-D).

保存图片、插图等的外部文件。对象名称为nihms-1523157-f0002.jpg
预计SULT2B1b KD可激活scRNA-seq中的TNF和NF-κB。

紫罗兰图表示对照KD中指定NF-κB靶基因的表达(读取计数)SULT2B1b KD组的scRNA-seq分析,调整后的p值显示在左上角。

保存图片、插图等的外部文件。对象名为nihms-1523157-f0003.jpg
SULT2B1b KD诱导LNCaP中TNF表达并激活NF-κB,但不激活C4-2细胞。

用对照或SULT2B1b siRNA转染LNCaP和C4-2细胞,并评估(A-B)TNF表达或(C-D)NF-κB活性。(A) qRT-PCR用于TNF公司siRNA KD后72小时表达。(B) siRNA KD后72小时,用基于细胞裂解物的ELISA定量TNF蛋白表达。(C) 荧光素酶检测在siRNA转染48小时后完成。NF-κB荧光素酶活性标准化为雷尼利亚荧光素酶。样品进一步标准化,以控制孔,无siRNA转染。横线表示使用Student t检验的至少四个独立实验的平均+/-SEM。(D) siRNA KD后72小时对指示蛋白进行Western blot。

为了确定SULT2B1b KD激活NF-κB是否是由于在LNCaP细胞中诱导TNF表达所致,对SULT2B1 B KD细胞中的TNF进行连续KD。这些数据表明,通过荧光素酶分析,SULT2B1b KD细胞中TNF表达的缺失消除了NF-κB活性的诱导(图4A-B). 了解到LNCaP细胞对TNF介导的细胞死亡敏感图3和4B4B类质疑SULT2B1b KD产生的TNF是否是这些细胞凋亡的原因。TNF中和实验是为了解决这个问题;这些研究表明,通过分析亚G1期细胞核,SULT2B1b-KD细胞中基于抗体的TNF中和作用对细胞死亡没有影响,而外源性TNF可以在这些细胞中有效中和(图4C). 因此,LNCaP细胞中的SULT2B1b KD产生TNF,导致NF-κB活化,但这些细胞中的TNF中和不能阻止SULT2B1 B KD诱导的凋亡。

保存图片、插图等的外部文件。对象名为nihms-1523157-f0004.jpg
SULT2B1b KD细胞中TNF的表达是NF-κB活化的原因,但不诱导细胞死亡。

(A) pRT-PCR指示SULT2B1b系列TNF公司连续转染SULT2B1b siRNA和TNF siRNA的LNCaP细胞中的表达。通过单因素方差分析和多重比较试验确定显著性。(B) NF-κB荧光素酶检测在siRNA转染48小时后完成,并显示相对荧光素酶单位(RLU)。条形图表示仅考虑对照KD和SULT2B1b KD组的双向方差分析得出的至少三份样品的平均+/-SEM,该图代表两个独立实验。(C) 用TNF中和抗体(40µg/mL)预处理LNCaP细胞,然后用对照或SULT2B1b siRNA转染或用指定浓度的外源性TNFα处理。在siRNA转染后约72小时,收获细胞并制备用于通过流式细胞术进行细胞周期分析。条形图表示三份样品的平均G1亚核+/-SEM百分比。该图代表了至少两个独立的实验。

SULT2B1b的表达改变了LNCaP和C4-2细胞对TNF介导的细胞死亡的敏感性。

虽然SULT2B1b KD通过TNF的表达激活NF-κB,但这似乎不是这些细胞凋亡的直接机制(图3——4). 数据显示于图3也证明SULT2B1b KD对LNCaP和C4-2细胞系中TNF表达和NF-κB活性有不同的影响。我们假设与死亡受体信号相关的基因上调会使SULT2B1b KD细胞对死亡配体刺激敏感。对SULT2B1b KD在LNCaP和C4–2细胞中的影响的进一步研究表明,SULT2B1b KD在二者都通过细胞周期分析确定LNCaP和C4-2细胞与对照组的比较(图5A-B). 当增加TNF浓度时,在两种细胞系的SULT2B1b KD细胞中也可以看到细胞死亡增强的剂量反应(图5A-B).

保存图片、插图等的外部文件。对象名为nihms-1523157-f0005.jpg
SULT2B1b调节LNCaP和C4-2细胞对TNF介导的细胞死亡的敏感性。

在添加指定浓度的TNF之前,用对照或SULT2B1b siRNA转染LNCaP(A)或C4-2(B)细胞。72小时后收集所有细胞,并通过流式细胞仪进行细胞周期分析。条形图表示三次孔的亚G1细胞核+/−SEM的平均百分比。对于(A)和(B)中的所有TNF处理细胞,通过双向方差分析,SULT2B1b KD显著提高LNCaP和C4-2细胞中G1以下细胞核的百分比(p<0.0001)。图是三个独立实验的代表。通过(C)qRT-PCR和(D)western blot评估LNCaP-SULT2B1b细胞中(C-D)SULT2B1b的表达。(E) 在添加指定浓度的TNF之前,用强力霉素处理LNCaP-SULT2B1b细胞。72小时后,收集所有细胞并通过流式细胞仪进行细胞周期分析。条形图表示通过双向方差分析和多重比较测试评估的三个孔的%亚G1核+/-SEM平均值。该图代表重复实验。

根据我们之前的报告,SULT2B1b KD降低了LNCaP细胞中AR的表达和活性,我们旨在确定TNF对LNCaP或C4-2细胞中SULT2Bb KD细胞的治疗是否会加剧AR活性的影响。在患有SULT2B1b KD的LNCaP细胞中,低剂量TNF处理导致SULT2B1 bKD和TNF处理在减少两者之间的显著交互作用AR公司(p=0.0361)和PSA(p=0.0421)表达,而通过双向方差分析,在C4-2细胞中SULT2B1b KD和TNF处理之间没有显著的相互作用(补充图3). 然而,与对照KD细胞相比,经SULT2B1b KD和TNF处理的C4-2细胞显著降低了PSA表达(补充图3B). 因此,尽管TNF和SULT2B1b KD能够独立抑制AR活性(6,17)虽然分子间的相互作用只发生在LNCaP细胞中,而不发生在CRPC模型中,但这种结合导致AR活性的最大降低。

由于流式细胞术未观察到TNFR1或TNFR2的显著增加,因此对TNF诱导的凋亡的敏感性增强似乎不是由于TNF受体表达的改变(补充图4). 先前的报告表明,TNF受体相关衔接蛋白(TRADD)表达的缺失是晚期前列腺癌细胞(如C4-2)对TNF介导的细胞死亡敏感性降低的原因。(6)SULT2B1b KD中TRADD表达的Western blot和ELISA分析对照KD细胞表明,SULT2B1b KD处理后LNCaP细胞中TRADD表达增加,但C4-2细胞中TRAND表达没有增加(补充图5). 因此,SULT2B1b可能调节TRADD水平,并可能解释LNCaP细胞对TNF敏感性的改变,但不能解释C4-2细胞对TNF-治疗敏感性的增强。

完成了额外的研究,以确定SULT2B1b KD介导的对TNF诱导的细胞死亡的敏感性增加是否仅仅是先前证明正在经历凋亡的细胞的伪影,或者SULT2B1 b是否是这种死亡受体途径的调节器。为了解决这个问题,通过添加多西环素,产生LNCaP-SULT2B1b细胞来诱导稳定的SULT2B1b过度表达。LNCaP-SULT2B1b的实验表明,SULT2B1b调节TNF介导的细胞死亡,因为与不诱导强力霉素相比,强力霉素的添加导致对外源性TNF治疗的敏感性显著降低(图5C-E). 这些研究证实了SULT2B1b作为前列腺癌细胞对TNF介导的细胞死亡敏感性的调节器的作用。

人类患者中SULT2B1b表达水平与肿瘤坏死因子相关基因呈负相关。

本文提供的数据表明,SULT2B1b KD诱导前列腺癌细胞的许多细胞内变化,特别是增加死亡受体信号的激活和对TNF的敏感性。为了研究SULT2B1b在临床环境中是否也可能在TNF敏感性或死亡受体信号传导中发挥作用,我们从cBioPortal获得了转移性前列腺癌的基因表达相关性数据。SULT2B1系列与骨转移样本中的多个TNF相关基因呈负相关,这些骨转移样本之前未接受过阿比特龙或恩扎鲁胺治疗,包括死亡受体配体TNF公司(r=−0.364)和CD40长度(r=-0.442),以及TNFR-相关基因贸易(r=-0.369)和TRAF2型(r=−0.439)(图6A). 此外,SULT2B1b与先前治疗的淋巴结衍生CRPC样本中的TNF相关基因呈负相关,包括NFKB1型(r=-0.533),凋亡调节因子FADD公司(r=−0.500)和死亡受体4(DR4,r=−0.544)(图6B). 这些数据表明,如本文所述的细胞系数据所示,SULT2B1b的表达与转移性前列腺癌和CRPC患者样本中的TNF和相关基因呈负相关。以下各项之间显著负相关的完整列表SULT2B1系列骨源性样本中的肿瘤坏死因子相关基因列于补充表1有趣的是,一个基因,TNFSF11型或RANKL与SULT2B1系列前列腺癌骨转移标本中r=0.746(图6C). 已知RANKL的高表达通过促进破骨细胞介导的骨吸收与骨转移生态位形成和癌症进展相关。(27,28)总之,这些数据表明,前列腺癌细胞中SULT2B1b的表达与TNF相关基因呈负相关,但与骨中破骨细胞激活和癌细胞存活相关的基因呈正相关。

保存图片、插图等的外部文件。对象名为nihms-1523157-f0006.jpg
在人类患者中,SULT2B1b的表达与许多TNF相关基因相关。

(A) 基因共表达表明SULT2B1系列与…负相关TNF公司,CD40长度,贸易、和TRAF2型来自未接受过治疗的人的骨髓源性转移性前列腺癌样本。(B) 共同表达SULT2B1系列具有NFKB1型,FADD公司和DR4确定了与淋巴结衍生CRPC样本中每个基因的显著负相关。(C) 基因共表达SULT2B1系列与RANKL在骨髓源性前列腺癌样本中显示出强烈的正相关。皮尔逊相关系数(r)和通过置换测试(p)确定的p值表示为每个相关图。圆点表示从单个患者中获得的样本。

讨论

我们在此提供证据证明,如scRNA-seq分析所检测到的,siRNA-mediated knowdown of SULT2B1b in LNCaP cells导致多种信号通路发生巨大变化。一个模型显示了SULT2B1b对死亡受体信号的拟议作用图7其中,SULT2B1b通过NF-κB和AR被证明是细胞死亡的关键调节剂。

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

SULT2B1b对晚期激素性前列腺癌和CRPC细胞中死亡受体信号和AR活性的影响示意图。

该分析的一个目的是了解各治疗组SULT2B1b表达的异质性。有趣的是,SULT2B1系列两组细胞内的表达水平均较低且不均匀(图1C)因此,在单细胞分析中,该基因的表达不能与其他感兴趣的基因直接相关。分析SULT2B1系列使用GTEx Portal在多种组织类型中的表达支持前列腺细胞中该基因的转录水平较低(补充图6)作为我们的数据图1这表明,但大量细胞零表达也可能是由于单细胞分析观察到的高脱落率。(29——31)即使表达SULT2B1系列相对较低,值得注意的是,SULT2B1b的丢失仍然很严重,有2000多个DE基因。

我们之前证明,SULT2B1b KD诱导LNCaP细胞凋亡,正如预期的那样,细胞内应激和凋亡相关通路显著改变,包括蛋白泛素化、线粒体功能障碍、BRCA1在DNA损伤反应中的作用、p53信号传导和HIV1诱导凋亡(补充图2). 令人惊讶的是,与对照KD细胞相比,SULT2B1b KD细胞在经典死亡受体信号通路的基因上产生了显著的改变,但在该通路中鉴定的基因发生了改变,包括TNFSF15型(TL1),TNFSF10型(TRAIL),以及DAXX公司基因与SULT2B1b KD诱导的凋亡无关(补充图7).

进一步研究在体外确定SULT2B1b-KD诱导TNF表达,导致LNCaP细胞中的NF-κB活化,并且这些作用由scRNA-seq数据预测。令人惊讶的是,SULT2B1b KD显著增加了四个抗凋亡基因的表达,说明SULT2B1 b的复杂细胞内效应(图2). SULT2B1b KD诱导LNCaP细胞表达TNF似乎不是SULT2B1b KD介导的凋亡的原因,因为TNF中和并不能消除凋亡(图4C). 因此,我们得出结论,SULT2B1b KD诱导TNF在蛋白水平上可能不会通过死亡受体直接诱导凋亡,但TNF表达增加确实会影响SULT2B1 bKD细胞的下游信号传导。SULT2B1b KD和过度表达研究(图5)表明SULT2B1b是前列腺癌对TNF介导的细胞死亡敏感性的新调节因子。中的数据进一步支持了这一点图6补充表1这表明SULT2B1b的表达与转移性前列腺癌中众多TNF-相关基因的表达呈负相关。

已经证明TNF信号可以降低AR的表达和活性。(32)我们之前显示SULT2B1b KD降低AR活性,但在CRPC中发生的程度较小(补充图3,(17)). SULT2B1b KD诱导的TNF表达不太可能是AR抑制的原因,因为没有检测到TNF蛋白的显著升高,TNF的中和作用对凋亡没有影响(图3——4). 因此,尚不清楚SULT2B1b调节AR的确切机制,但SULT2B1 bKD确实有助于降低TNF-处理的C4-2细胞中PSA的表达。也许SULT2B1b KD引起细胞内变化,使TNF治疗对CRPC中AR活性有更大的影响。

虽然近年来肿瘤坏死因子的治疗应用有所减少,但它仍然是肿瘤微环境中的一个重要分子,可能会影响治疗反应。与局限性前列腺癌患者相比,CRPC患者的循环TNF水平升高。(33)如果SULT2B1b的表达/活性如本文数据所示保持前列腺癌细胞对TNF的耐药性,那么靶向SULT2B1可能会使晚期前列腺癌细胞在微环境中对TNF敏感。然而,值得注意的是,TNF本身可能会影响SULT2B1b的表达,因为低剂量TNF治疗会降低C4-2细胞中SULT2B1的表达(补充图3B,第五天). 因此,微环境中升高的TNF可能抑制CRPC中SULT2B1b的表达。了解这种疾病患者的炎症状态对于确定这种靶向方法是否会改变癌细胞的生存能力/细胞毒性平衡至关重要,因为SULT2B1b和TNF的高表达可能对治疗效果是必要的。

SULT2B1b可能通过其作为LXR调节器在胆固醇代谢中的功能影响死亡受体信号通路。(17,34)众所周知,富含胆固醇的脂筏可以稳定细胞内信号,但也可能是死亡受体信号的起始位点。(16,25)虽然在LNCaP细胞中,SULT2B1b KD不会显著改变总胆固醇水平(数据未显示),但SULT2B1 KD确实会降低细胞内硫酸胆固醇水平。(17)尚不清楚硫酸胆固醇(未修饰胆固醇)对脂筏稳定性很重要,但脂筏中的扰动肯定会导致许多典型信号通路的大量中断,包括TNFR复合物的形成和信号传导。(35)

本文所证明的数据也提供了证据,表明SULT2B1b的表达与多个TNF相关基因的表达呈负相关,表明SULT 2B1b在转移性激素血症和去势抵抗前列腺癌中具有保护作用,但SULT2B1也与破骨细胞调节剂RANKL呈强正相关(图6C). 这些数据表明,SULT2B1b可能参与躲避死亡配体诱导的前列腺癌细胞凋亡,以及在疾病进展过程中建立前列腺癌骨转移生态位。这支持了SULT2B1b在激素性前列腺癌患者骨内转移病灶中的作用。靶向SULT2B1b可能会限制破骨细胞的形成,但该策略可能需要与其他靶向策略结合使用,以抵消转移性疾病中激活的信号通路,例如使用他汀类药物的脂筏稳定性。(36——38)在TNF水平较高的患者中,如CRPC患者,SULT2B1b靶向可能有助于使CRPC细胞对外源性TNF敏感和/或抑制AR活性,从而导致细胞死亡。SULT2B1b应被视为晚期前列腺癌的治疗靶点,其表达状态可与其他信息(即TNF水平、炎症状态、他汀类药物使用或其他治疗史等)一起用于指导患者的治疗选择。

启示:这些数据表明,SULT2B1b的表达增强了对TNF的抵抗力,并可能促进前列腺癌的发生。

补充材料

1

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

2

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

致谢

这些研究得到了国防部前列腺癌研究项目拨款#W81XWH-14–1-0588(TL Ratliff、AD Mesecar、SA Crist和RE Vickman)和沃尔特癌症基金会(TL Ratliff)的支持。作者感谢前列腺癌基金会的慷慨支持和DDC Medical对我们的细胞系进行认证。我们感谢普渡基因组核心设施(特别是菲利普·桑米格尔和保罗·帕克)、普渡流式细胞术和细胞分离设施、生物信息学共享资源、普渡大学癌症研究中心(NIH拨款P30 CA023168)、IU西蒙癌症中心(NIH-拨款P30 CA 082709,NA Lanman)、,以及沃尔瑟癌症基金会。我们也非常感谢Bennett Elzey博士和Ratliff实验室成员对该项目的持续建议和支持。

缩写列表:

方差分析方差分析
AR公司雄激素受体
中国铁建股份有限公司去势抵抗前列腺癌
DAXX公司死亡结构域相关蛋白
判定元件差异表达
FBS公司胎牛血清
财务总监错误发现率
投资促进机构创意路径分析
杜兰特击倒
LXR公司肝X受体
PSA公司前列腺特异性抗原
RANKL公司TNF配体超家族成员11
RNA干扰RNA干扰
每公里公里数每千字节读取数/百万映射读取数
scRNA-seq基因单细胞mRNA测序
SULT2B1b系列硫转移酶2B1b;胆固醇硫转移酶
TL1级TNF超级家族成员15
TNF公司肿瘤坏死因子
TNFR公司肿瘤坏死因子受体
贸易TNF受体相关死亡域
TRAIL(跟踪)肿瘤坏死因子相关凋亡诱导配体

脚注

利益冲突声明:这份手稿的作者声明没有潜在的利益冲突。

参考文献

1AC学会(2018)2018年癌症事实与数据(亚特兰大:美国癌症学会;)。[谷歌学者]
2An J等人(2003)蛋白酶体抑制剂硼替佐米与前列腺癌细胞毒性化疗、肿瘤坏死因子α和肿瘤坏死因子相关凋亡诱导配体的药物相互作用.临床癌症研究 9(12):4537–4545. [公共医学][谷歌学者]
三。Norris JS等人(2001年)Fas配体、TRAIL和Bax在前列腺癌基因治疗中的应用.Curr基因治疗 1(1):123–136. [公共医学][谷歌学者]
4O'Kane HF等人(2006年)针对膀胱癌、前列腺癌和肾癌中的死亡受体.泌尿学杂志 175(2):432–438. [公共医学][谷歌学者]
5Carswell EA等人(1975年)一种内毒素诱导的导致肿瘤坏死的血清因子.美国国家科学院程序 72(9):3666–3670.[PMC免费文章][公共医学][谷歌学者]
6王德等(2009)肿瘤坏死因子受体相关死亡域表达降低与前列腺癌进展相关.癌症研究 69(24):9448–9456.[PMC免费文章][公共医学][谷歌学者]
7Chopra DP、Menard RE、Januszewski J和Mattingly RR(2004)TNF-α介导的正常前列腺上皮细胞和肿瘤细胞凋亡.癌症快报 203(2):145–154. [公共医学][谷歌学者]
8Kimura K、Bowen C、Spiegel S和Gelmann EP(1999)肿瘤坏死因子α使前列腺癌细胞对γ射线诱导的凋亡敏感.癌症研究 59(7):1606–1614。[公共医学][谷歌学者]
9Chung TD等人(1998)肿瘤坏死因子α基因治疗增强前列腺癌的放射细胞毒性.癌症基因治疗 5(6):344–349. [公共医学][谷歌学者]
10Pirtskhalaishvili G、Shurin GV、Esche C、Trump DL和Shurin MR(2001)TNF-α保护树突状细胞免受前列腺癌诱导的凋亡.前列腺癌前列腺疾病 4(4):221–227. [公共医学][谷歌学者]
11Harada S等人(2001年)长期接触肿瘤坏死因子α导致LNCaP前列腺癌细胞对雄激素过敏和抗雄激素戒断现象.前列腺 46(4):319–326. [公共医学][谷歌学者]
12Mizokami A、Gotoh A、Yamada H、Keller ET和Matsumoto T(2000)肿瘤坏死因子α抑制LNCaP前列腺癌细胞系雄激素敏感性.泌尿学杂志 164(第3部分第1部分):800–805. [公共医学][谷歌学者]
13卢力等(2012)甘草酸通过PI3K/Akt和NF-kappaB信号通路抑制TNF-α诱导的人前列腺癌PC3细胞体外侵袭.中国药理学报 33(4):531–541.[PMC免费文章][公共医学][谷歌学者]
14Radhakrishnan P等人(2011年)TNFalpha通过刺激参与选择素配体合成的选择性糖基和硫转移酶基因的表达,增强前列腺癌细胞的运动性和侵袭性.生物化学-生物物理研究委员会 409():436–441.[PMC免费文章][公共医学][谷歌学者]
15岳S等(2014)PTEN缺失和PI3K/AKT激活诱导的胆固醇积累是前列腺癌侵袭性的基础.单元格元 19():393–406.[PMC免费文章][公共医学][谷歌学者]
16弗里曼MR和所罗门KR(2004)胆固醇与前列腺癌.J细胞生物化学 91(1):54–69. [公共医学][谷歌学者]
17Vickman RE等人(2016)胆固醇磺化酶SULT2B1b调节前列腺癌AR和细胞生长特性.摩尔癌症研究 14(9):776–786。[PMC免费文章][公共医学][谷歌学者]
18Hsu CC和Hu CD(2012)N-末端定位核输出信号在调节激活转录因子2(ATF2)亚细胞定位和转录活性中的关键作用.生物化学杂志 287(11):8621–8632.[PMC免费文章][公共医学][谷歌学者]
19Gordon AaH G(2010年)FASTX工具包.[谷歌学者]
20安德鲁斯·S(2010)快速质量控制.[谷歌学者]
21Kim D等人(2013)TopHat2:存在插入、缺失和基因融合时转录组的精确比对.基因组生物学 14(4):R36。[PMC免费文章][公共医学][谷歌学者]
22Trapnell C、Pachter L和Salzberg SL(2009年)TopHat:使用RNA-Seq发现拼接接头.生物信息学 25(9):1105–1111.[PMC免费文章][公共医学][谷歌学者]
23Anders S、Pyl PT和Huber W(2015)HTSeq——一个处理高通量测序数据的Python框架.生物信息学 31(2):166–169.[PMC免费文章][公共医学][谷歌学者]
24Robinson D等人(2015)晚期前列腺癌的综合临床基因组学.单元格 161(5):1215–1228.[PMC免费文章][公共医学][谷歌学者]
25Legler DF、Micheau O、Doucey MA、Tschopp J和Bron C(2003)向脂筏招募TNF受体1对TNFα介导的NF-kappaB活化至关重要.免疫 18(5):655–664. [公共医学][谷歌学者]
26Lewis AK等人(2016)死亡受体5网络需要膜胆固醇才能实现正确的结构和功能.分子生物学杂志 428(24点A):4843–4855.[PMC免费文章][公共医学][谷歌学者]
27Jones DH等人(2006)RANKL对癌细胞迁移和骨转移的调控.自然 440(7084):692–696. [公共医学][谷歌学者]
28李X等(2014)OPG/RANK/RANKL轴在前列腺癌侵袭和骨转移中的潜在作用.及国际医学权威期刊 32(6):2605–2611. [公共医学][谷歌学者]
29李伟伟、李JJ(2018)一种精确稳健的单细胞RNA-seq数据插补方法scImpute.国家公社 9(1):997.[PMC免费文章][公共医学][谷歌学者]
30Ziegenhain C等人(2017)单细胞RNA测序方法的比较分析.摩尔细胞 65(4):631–643 e634。[公共医学][谷歌学者]
31Kharchenko PV、Silberstein L和Scadden DT(2014)单细胞差异表达分析的贝叶斯方法.Nat方法 11(7):740–742.[PMC免费文章][公共医学][谷歌学者]
32Ko S、Shi L、Kim S、Song CS和Chatterjee B(2008)核因子-kappaB和B-myb在肿瘤坏死因子α负调控雄激素受体表达中的相互作用.分子内分泌学 22(2):273–286。[PMC免费文章][公共医学][谷歌学者]
33Chadha KC等人(2014)前列腺癌诊断的新血清生物标志物.临床癌症研究杂志 (1):72–79.[PMC免费文章][公共医学][谷歌学者]
34Chen W、Chen G、DL主管、Mangelsdorf DJ和Russell DW(2007)酶法降低氧甾醇对培养细胞和小鼠肝脏LXR信号转导的影响.单元格元 5(1):73–79.[PMC免费文章][公共医学][谷歌学者]
35Muppidi JR、Tschopp J和Siegel RM(2004)生死决定:TNF受体家族信号转导中的二级复合物和脂筏.免疫 21(4):461–465. [公共医学][谷歌学者]
36庄L、金J、亚当RM、所罗门KR和弗里曼MR(2005)胆固醇靶向改变前列腺癌细胞和异种移植物的脂筏组成和细胞存活率.临床研究杂志 115(4):959–968.[PMC免费文章][公共医学][谷歌学者]
37庄L、林J、卢ML、所罗门KR和弗里曼MR(2002)富含胆固醇的脂筏介导akt调节的前列腺癌细胞存活.癌症研究 62(8):2227–2231. [公共医学][谷歌学者]
38Srivatsav AT、Mishra M和Kapoor S(2018年)脂依赖性细胞过程的小分子调控抗癌:枪尖上的脂肪.生物识别识别 2018:6437371.[PMC免费文章][公共医学][谷歌学者]