跳到主要内容
访问密钥 NCBI主页 MyNCBI主页 主要内容 主导航
自然细胞生物学。作者手稿;PMC 2019年1月18日提供。
以最终编辑形式发布为:
预防性维修识别码:项目经理6145444
EMSID:EMS78112标准
PMID:30022119

屏蔽蛋白复合物促进BRCA1-null细胞中DNA末端连接并对抗同源重组

关联数据

补充资料

摘要

BRCA1缺陷会导致乳腺癌、卵巢癌、前列腺癌和其他癌症,并使肿瘤对PARP抑制剂过敏。为了了解耐药机制,我们对使用PARP抑制剂治疗的BRCA1缺陷乳腺癌细胞进行了全基因组CRISPR-Cas9合成活性/耐药筛选。我们鉴定了两种以前没有特征化的蛋白质,C20orf196和FAM35A,它们的失活导致强烈的PARP抑制剂耐药性。从力学角度来看,我们发现C20orf196和FAM35A形成了一个复合物“Shieldin”(SHLD1/2),FAM35A-通过其C末端OB折叠区与单链DNA相互作用。我们确定,Shieldin作为53BP1/RIF1/MAD2L2的下游效应器,通过限制DSB切除来促进DNA双链断裂(DSB)末端连接,并通过拮抗BRCA1缺陷细胞中BRCA2/RAD51负载来对抗同源重组。值得注意的是,屏蔽蛋白失活进一步使BRCA1缺陷细胞对顺铂敏感,这表明确定BRCA1缺乏肿瘤的SHLD1/2状态可能有助于患者分层并产生新的治疗机会。为了突出这一潜力,我们记录了SHLD1/2在显示固有或获得性PARP抑制剂耐药性的人类乳腺癌中的表达降低。

介绍

DNA双链断裂(DSB)是一种高度细胞毒性的细胞损伤,必须进行有效和准确的修复,以保持基因组稳定性,防止过早衰老、神经退行性变、免疫缺陷、癌症和其他疾病1为了响应DSB检测,顶端激酶ATM、ATR和PRKDC(DNA-PKcs)被激活并磷酸化许多底物以启动细胞DNA损伤反应(DDR)4随后的分子DDR事件的级联,由各种翻译后修饰促进,包括蛋白质磷酸化、泛素化、sumoylation和聚(ADP-核糖基)化,影响无数的细胞成分,以及其他导致DDR因子在DNA损伤位点的组装的因素,细胞周期进展的阻滞或减慢,以及DNA修复机制的激活4,5DSB修复途径的两种主要类型是在整个细胞周期中活跃的非同源末端连接(NHEJ)和同源重组(HR),后者通常需要姐妹染色单体作为模板,因此仅在细胞周期的S期和G2期运行。DSB修复途径的选择部分由HR促进因子BRCA1和NHEJ促进蛋白如TP53BP1(53BP1)、RIF1和MAD2L2(REV7)之间的功能拮抗作用决定613.

遗传或获得性突变巴西航空公司1巴西航空公司2导致蛋白质丢失或突变BRCA1/2蛋白的基因会导致乳腺癌、卵巢癌、前列腺癌和其他癌症,并使肿瘤对PARP抑制剂药物(如奥拉帕林)过敏1417不幸的是,固有或获得性PARP抑制剂耐药性经常导致反应迟钝或患者复发和肿瘤再生15,18在临床上,迄今为止报道的最常见PARP抑制剂耐药机制是BRCA1/2表达或功能的恢复。值得注意的是,53BP1在各种三阴性乳腺癌中表达缺失7这可能是PARP抑制剂耐药的某些临床相关示例的原因。然而,在大部分BRCA1/2缺陷肿瘤中,导致PARP抑制剂耐药的机制仍不清楚18,19.

为了系统地研究PARP抑制剂耐药性的遗传机制,我们对使用PARP抑制剂治疗的人类BRCA1缺陷乳腺癌细胞进行了全基因组CRISPR-Cas9合成活性/耐药性筛选。除了识别已知的阻力系数,如53BP1、RIF1和MAD2L2损失613,我们鉴定了两个先前未经特征化的蛋白质,C20orf196和FAM35A,其失活使PARP抑制剂对BRCA1缺陷细胞产生耐药性。我们接下来的工作引导我们定义“Shieldin”(SHLD1C20或196/SHLD2型FAM35A系列)通过充当53BP1、RIF1和MAD2L2的下游效应器促进NHEJ的复合物,限制DSB切除,并通过拮抗BRCA2和RAD51对切除的单链DNA(ssDNA)上的复制蛋白A(RPA)的替换,对抗BRCA1缺陷细胞中的HR。最后,我们报告SHLD1C20或196/SHLD2型家庭35a丢失使DNA交联剂顺铂过敏,并降低SHLD1C20或196或SHLD2FAM35A系列在患者衍生的BRCA1缺陷乳腺癌异种移植模型和表现出固有PARP抑制剂耐药性的BRCA1-突变癌症中,PARP抑制剂的表达与PARP抑制剂抗性的演变有关。

结果

FAM35A或C20orf196丢失抑制PARP抑制剂的敏感性巴西航空公司1-突变细胞

为了系统地探讨PARP抑制剂耐药性的遗传机制,我们利用GeCKO文库进行了全基因组CRISPR-Cas9基因激活筛选20在中巴西航空公司1-突变型乳腺癌细胞株SUM149PT与PARP抑制剂olaparib、talazoparib(BMN673)或AZD2461平行治疗(图1a、b,补充图1a-c). 除了识别已知的抗性基因外第53页/页,RIF1级MAD2L2型其产品形成一个复合体21,我们确定了几个新的抑制物候选物(补充表1,补充图1d-e). 其中包括DYNL1,一个已知的53BP1交互合作伙伴22和TEN1,CST端粒包被复合物(CTC1/STN1/TEN1)的一种成分,也促进端粒DNA复制23然而,在我们随后的研究中,我们将重点放在非特异性蛋白质FAM35A和C20orf196上,这些蛋白质在我们的筛选中获得了最高的分数(图1b补充表1). 因此,通过在非转化的hTERT永生化人类RPE1细胞中执行短干扰RNA(siRNA)介导的mRNA沉默(补充图1f),我们确定,对于53BP1损失7FAM35A或C20orf196的缺失显著抑制了BRCA1失活引起的PARP抑制剂敏感性,而对BRCA1阳性细胞无明显影响(图1c,补充图1g). 这一结论通过从头开始的CRISPR-Cas9基因编辑得到独立证实,FAM35A或C20orf196失活可以通过重新导入野生型FAM35A或C20or f196来抵消,从而减轻BRCA1缺陷细胞的olaparib超敏反应(图1d-e;补充图1h; 如所示补充图1i,这些影响并没有反映出细胞周期的改变)。

保存图片、插图等的外部文件。对象名称为emss-78112-f001.jpg
CRISPR-Cas9筛查确定PARP-抑制剂敏感性的抑制物巴西航空公司1-突变细胞。

a、,屏幕程序示意图。b、,特定药物治疗后指南富集的MAGeCK分析;虚线表示的错误发现率(FDR)为0.1;n=每个药物治疗3个技术重复。c、,siRNA介导的克隆生存分析中hits的验证;下部面板显示曲线下面积(AUC);n=3个独立实验日期:,克隆形成生存分析中FAM35A的从头Cas9介导的敲除(ko)验证和补充(AUC中显示了多个ko克隆);n=4个独立实验,FAM35Ako(#14)(n=2)、FAM35Ako(#40)(n=3)、BRCA1ko/FAM35Ako(#34)(n=2)和BRCA1ko/FAM35Ako(#2)+FAM35A(n=3)除外。e、,如(d)所示,但适用于C20或196;n=3个独立实验,BRCA1ko/C20或196ko+C20或196除外(n=2)。c-e公司条形代表平均值±SEM,单向方差*p<0.05,**p<0.01,***p<0.001,***p<0.05,ns=不显著(p≥0.05)。各个数据点绘制在条形图上,包括精确p值在内的统计源数据可以在补充表5.

FAM35A/C20orf196复合物与53BP1/RIF1/MAD2L2相互作用并作用于其下游

序列分析表明,FAM35A和C20orf196在脊椎动物中具有很好的保守性。此外,结构预测建模(RaptorX;http://raptorx.uchicago.edu/)发现FAM35A含有一个无序的N末端和一个有序的C末端区域,其中包含三个OB折叠,最后一个C末端OB折叠/FAM结构域包含CXXC型锌指基序(图2a). 值得注意的是,该组织与ssDNA结合蛋白RPA的RPA1亚基和也结合ssDNA的CST复合物的CTC1亚基非常相似23在这方面,我们注意到,虽然C20或196 N末端(残基1-70)被预测为内在无序,但其C末端部分结构更为复杂,可能含有单翼或双翼螺旋(WH)结构域(图2a)与酵母CST亚单位Stn1中的类似23,24表明C20orf196和Stn1可能发挥类似或互补的作用。

保存图片、插图等的外部文件。对象名称为emss-78112-f002.jpg
FAM35A和C20或196域、相互作用和IRIF形成。

a、,FAM35A和C20或使用的196个预测域和变体、OB折叠(OB)、FAM域(OB3/FD)。b、,通过mCherry-LacR-C20orf196将FAM35A/衍生物GFP-融合物招募到染色体Lacooperator阵列。所示数据代表3个实验,量化如所示补充图2a.比例尺10µm。c、,(左侧和中间面板)纯化的重组GST-FAM35A直接与重组His-C20orf196相互作用。c、 (右图)通过联合免疫沉淀和免疫印迹分析表达GFP-FAM35A/衍生物和HA-C20orf196的细胞提取物。日期:,V5-FAM35A与GFP-MAD2L2共同免疫沉淀;与C20或196的相互作用如所示补充图2c.e、,IR(5Gy)处理后5 hγH2AX阳性细胞中诱导型GFP-FAM35A(左侧)和GFP-C20orf196(右侧)IRIF的定量。N=4个独立实验,除了(左面板)si53BP1(N=3)、siRIF1和siMAD2L2(N=2);以及(右侧面板)siCTRL(n=5)、siRIF1(n=3)、siFAM35A(n=3)。f、,如(e)所示,但对于可诱导的GFP-FAM35A N末端;n=4个独立实验,siRIF1除外(n=3)。克,与GFP-FAM35A N末端的内源性MAD2L2共同免疫沉淀。e-f、,条形代表平均值±SEM,单向方差*p<0.05,**p<0.01,***p<0.001,***p<0.05,ns=不显著(p≥0.05);在条形图上绘制的单个数据点。统计源数据(包括精确的p值)如所示补充表5所有免疫印迹都是两个独立实验的代表;免疫印迹的未处理扫描如所示补充图8.

通过结合细胞共定位和共免疫沉淀实验,我们确定FAM35A和C20orf196直接相互作用的方式主要是,但不完全是由FAM35AOB3/FAM结构域介导的(图2b-c,补充图2a-b). 由于FAM35A或C20orf196的缺失与BRCA1缺陷细胞中53BP1/RIF1/MAD2L2的缺失有相似的影响,我们测试了这些因素之间可能的相互作用。因此,通过联合免疫沉淀和质谱(MS)研究,我们发现C20orf196和FAM35A都与MAD2L2相互作用,后者是53BP1/RIF1/MAD2L2轴的最远端因子,介导BRCA1缺陷细胞中PARP抑制剂敏感性613(图2d,补充图2c).

在电离辐射(IR)诱导的核病灶(IRIF)中,许多DDR蛋白积聚在DSB位点5我们确定FAM35A和C20orf196均形成IRIF,并且通过活细胞成像研究发现,这些蛋白质也被招募到激光微辐射诱导的局部DNA损伤位点(补充图2d). 此外,我们通过共焦显微镜和超分辨显微镜确定FAM35A与已建立的DSB标记物磷酸化组蛋白H2AFX共定位25(γH2AX)和53BP126(图2e-f,补充图2e). 值得注意的是,siRNA/shRNA-缺失实验证明,虽然53BP1 IRIF和MAD2L2水平以及IRIF没有因FAM35A或C20或196缺失而显著受损(补充图2f-h)FAM35A和C20或196形成IRIF需要53BP1、RIF1和MAD2L2,但不需要PTIP(图2e补充图3a-c; 注意,GFP标记的FAM35A/C20orf196的总水平受53BP1/RIF1/MAD2L2耗竭的影响最小)。我们还确定,C20或196 IRIF几乎被FAM35A耗尽完全废除,而C20或1196耗尽减少,但没有废除FAM35AIRIF(图2e). 此外,当用DNA拓扑异构酶I抑制剂喜树碱(CPT)处理细胞时,FAM35A形成核病灶;补充图3d). 值得注意的是,FAM35A的N末端对其IRIF的形成是必要的和充分的,这些IRIF依赖于53BP1、RIF1、MAD2L2和C20orf196,并且该区域可以与MAD2共同免疫沉淀(图2f-g,补充图3e-f). 总之,这些发现表明FAM35A和C20orf196作为53BP1/RIF1/MAD2L2分子组装的下游组件21DSB现场。

FAM35A和C20或196促进NHEJ

自53BP1、RIF1和MAD2L2推动NHEJ以来613,我们测试了FAM35A和C20或196是否发挥了类似的作用。事实上,对于NHEJ因子XRCC4的缺失,53BP1、FAM35A或C20orf196的siRNA缺失通过质粒DNA随机整合到染色体中来检测NHEJ受损27(图3a). 此外,FAM35A或C20或196缺失使人和小鼠细胞对IR过敏(图3b补充图4a). 53BP1和相关因子促进NHEJ介导的免疫球蛋白重链位点的类开关重组(CSR),这是一个允许B淋巴细胞将抗体产生从一种类型改变为另一种类型的过程28通过CRISPR-Cas9基因编辑小鼠CH12F3(CH12)B淋巴细胞29我们确定,对于53BP1/RIF1/MAD2L2失活6,8,1012,30,31FAM35A或C20或196的丢失显著降低了CSR(图3c-d补充图4b-d). 此外,对这些细胞中期染色体扩散的分析表明,FAM35A或C20orf196失活导致染色体断裂和异常CSR的易位症状32 (图3e-f; 来自的备注补充图4e-fCSR效应与细胞增殖缺陷、Aid或种系Sα开关区转录缺陷无关)。

保存图片、插图等的外部文件。对象名称为emss-78112-f003.jpg
FAM35A和C20orf196促进NHEJ和免疫球蛋白类开关重组。

a、,随机质粒整合试验。b、,用IR处理FAM35Ako和C20orf196ko细胞,并分析其克隆存活率,右侧显示AUC。a-b中,条形代表平均值±SEM,单向方差;n=3个独立实验,b中的C20或196ko除外(n=4),单个数据点绘制在条形图上;统计源数据可以在中找到补充表5.c、,小鼠IgM类开关重组和染色体不稳定性的示意图+B细胞(带有C的生殖系配置μ转录)诱导表达AID并进行CSR至IgA(与C的开关配置α转录)。CSR水平以细胞因子刺激72小时后IgA阳性细胞的百分比来测量,DNA荧光原位杂交(FISH)使用12号染色体特异性涂料(灰色染色体)和伊格用于测量染色体不稳定性的位点特异性探针(红点和绿点)伊格CSR诱导位点。日期:,与野生型(WT)CH12-Cas9细胞相比,Fam35Ako和C20orf196ko CH12-Cas8细胞的CSR水平降低。条形表示平均值±SEM,单向Anova。3个克隆的N=4个独立实验,但53BP1ko+细胞因子除外,其中N=3个克隆,53BP1ko-细胞因子,其中N=2个克隆;在条形图上绘制单个数据点。e、,的代表性图像伊格异常中期的易位和断裂,在f中定量。f、,量化伊格所示CH12-Cas9细胞中期的断裂和移位。水平条代表平均值,Fisher精确检验;n=2个独立实验,Fam35ako和C20orf196ko除外,其中n=3。对于a、b、d和f,*p<0.05,**p<0.01,***p<0.001,****p<0.0001,ns=不显著(p≥0.05);统计源数据,包括这些面板的精确p值,可以在补充表5.

FAM35A和C20orf196拮抗DNA末端切除

为了进一步探索FAM35A和C20orf196的功能,我们在含有端粒相关因子TRF2(TRF2ts)的温度敏感等位基因的小鼠细胞中进行了检测。TRF2ts在高温下失活导致染色体末端失去保护,并导致NHEJ介导的端粒融合(图4a)6,33令人惊讶的是,通过使用短发夹RNA(shRNA)介导的mRNA沉默,我们发现FAM35A或C20或196缺失减少了这种染色体融合,就像MAD2L2缺失引起的那样(图4b,补充图5a-b).

保存图片、插图等的外部文件。对象名称为emss-78112-f004.jpg
FAM35A和C20orf196促进端粒介导的融合并限制DNA末端切除。

a、,TRF2ts实验装置示意图。b、,shRNA缺失FAM35A(左面板)或C20或196(右面板)可减少未标记端粒介导的染色体融合。横线表示平均值。实验进行了两次,每个条件下计数的染色体数≥1300条,并在条形图上绘制了个人数据点;源数据可以在中找到补充表5.c、,FAM35Ako和C20orf196ko RPE1细胞用BrdU(10μM)标记48小时,然后用1µM喜树碱(CPT)处理1小时,预提取、固定并在非变性条件下对BrdU进行染色,以显示ssDNA。中心线位于中线、箱限位于25的框线和胡须图第个/75第个百分位数和胡须±1.5xIQR;单向Anova;n=3个独立实验。日期:,由于Fam35a或C20或f196沉默,IR诱导的pRPA(S4/8)在MEF中增强。横线表示平均值。实验进行了两次,每个数据点绘制在条形图上;源数据可以在中找到补充表5.e、,RPE1-FAM35Ako或-C20orf196ko细胞显示超DNA末端切除(用1μM喜树碱处理细胞1h)。来自3个独立实验的代表性图像。比例尺10µm。(f)RPE1-FAM35Ako或-C20orf196ko细胞显示过度DNA末端切除的BLM和CtIP依赖性标记。中心线位于中线、箱限位于25的框线和胡须图第个/75第个百分位数和胡须±1.5xIQR;单向Anova;n=3个独立实验。克,与激光微辐射后2小时固定和染色的野生型(WT)RPE1细胞相比,FAM35Ako和C20orf196ko细胞BLM增加。左侧面板中显示的代表性图像和右侧面板中的量化图像。比例尺10µm。中心线位于中线、箱限位于25的框线和胡须图第个/75第个百分位数和胡须±1.5xIQR;单向Anova;n=3个独立实验。对于c、f和g,*p<0.05,**p<0.01,***p<0.001,****p<0.0001,ns=不显著(p≥0.05);包含精确p值的统计源数据可以在补充表5.

TRF2ts系统中53BP1、RIF1或MAD2L2缺失的影响与这些抵消DSB切除的因素有关613,34因此,我们探讨了FAM35A和C20orf196是否也具有此功能。事实上,对于53BP1/RIF1/MAD2L2失活613,35用RPA和喜树碱处理后预提取细胞核中ssDNA染色强度测定人类细胞中FAM35A或C20orf196的缺失增强了DSB切除(图4c-e; FAM35A、C20或196缺失不会改变激光微辐射诱导DNA损伤部位的RPA1动力学,补充图5c). 根据这种由典型通路介导的切除,它被切除促进因子RBBP8(CtIP)或BLM的缺失所减少(图4f). 此外,我们还确定,FAM35A或C20或196耗尽增强了BLM向激光微辐射位点的募集(图4g). 类似地,如小鼠细胞中53BP1缺失所示6IR处理后,这些细胞中FAM35A或C20orf196缺失导致DNA末端切除标记Ser4/8磷酸化RPA2水平升高(图4d). 结合我们的其他发现,这些数据确定FAM35A和C20orf196是53BP1/RIF1/MAD2L2介导的染色体NHEJ的关键成分,并表明其前NHEJ功能与限制DSB切除有关。

FAM35A OB折叠区与ssDNA相互作用并促进IR存活

与我们对FAM35A和RPA1之间结构相似性的预测一致,FAM35AC末端可以通过与ssDNA寡核苷酸的相互作用从细胞提取物中检索到(图5a). 根据RPA与ssDNA结合的关键类似残基,对RPA1的序列比对和FAM35A的结构建模确定了两个Trp(W)残基,预测它们位于蛋白-ssDNA界面(图5b,补充图5d). 根据这一预测,我们通过电泳凝胶迁移率分析(EMSA)发现,细菌表达的纯化的FAM35A C末端区域优先与ssDNA结合,而不是与双链DNA结合(图5c,补充图5e)当两个Trp残基突变为Ala(W489/W640A;图5c). 此外,全长FAM35A携带这些突变W489/W640A型)仍与C20或196交互(补充图5f)并且在细胞中形成IRIF,这些IRIF始终低于野生型FAM35A蛋白(图5d). 这表明,在通过其N末端区域招募IRIF后,FAM35A C末端ssDNA结合区域可能允许进一步招募、保留和/或稳定FAM35A。此外,与野生型蛋白质不同,FAM35AW489/W640A型当重新引入FAM35A无效细胞时,没有产生显著的IR抗性(图5e). 在平行研究中,FAM35A C末端的表达并不补充FAM35A-空细胞的IR超敏反应。此外,无论细胞是否表达内源性FAM35A,FAM35A N端的表达都会使细胞对IR过敏,这意味着FAM35A的N端IRIF形成结构域可能对NHEJ具有显性负面影响(图5f; 在野生型背景下,FAM35A衍生物的过度表达不会影响奥拉帕林的敏感性,补充图5g).

保存图片、插图等的外部文件。对象名称为emss-78112-f005.jpg
FAM35A OB折叠介导ssDNA相互作用,是IR抗性所必需的。

a、,残基W489/W640突变为A的FAM35A示意图(顶面板)。W489和W640位置野生型FAM35A的预测3D结构(左下面板)。FAM35A W489/W640促进细胞提取物中ssDNA的高效结合(右下角)。b、,yRPA1与FAM35A C端子对齐;对yRPA1 ssDNA结合至关重要的氨基酸和FAM35A中相应的氨基酸残基被装箱。c(c),EMSA在含有10nM ssDNA或dsDNA的天然(非变性)凝胶上,以及所示纯化、细菌表达的FAM35A C末端或W489/W640A突变的µM中。日期:,可诱导GFP-FAM35A W489/W640A未能有效形成IRIF(IR 5Gy后12小时)。比例尺10µm。来自3个独立实验的代表性图像。e、,在克隆存活试验中补充FAM35A衍生物的FAM35Ako RPE1细胞;右侧面板显示AUC。f、,在克隆形成分析中,FAM35A N端而非C端或全长的过度表达使野生型细胞对IR敏感;右侧面板显示AUC。e-f、,条形代表平均值±SEM,单向方差*p<0.05,**p<0.01,***p<0.001,***p<0.0001,ns=无显著性(p≥0.05);n=3个独立实验,e中第2组和第4组除外(n=2),单个数据点绘制在条形图上;包含精确p值的统计源数据可以在补充表5所有免疫印迹都是两个独立实验的代表;免疫印迹的未处理扫描如所示补充图8.

FAM35A和C20orf196拮抗BRCA1缺陷细胞的HR

PARP抑制剂产生复制相关的DNA损伤,需要BRCA1介导的HR才能有效修复1953BP1/RIF1/MAD2L2的缺失部分恢复了BRCA1缺陷细胞修复此类损伤的能力613这导致了一个模型,其中BRCA1和53BP1/RIF1/MAD2L2分别在引导DSB流向HR或NHEJ方面发挥拮抗作用。因此,我们推测BRCA1可能拮抗FAM35A/C20或196的作用。因此,FAM35A和C20或196 IRIF,而不是53BP1 IRIF在BRCA1耗竭后在数量和强度上显著增强,但在BRCA2耗竭后没有显著增强(图6a-b,补充图6a-b).

保存图片、插图等的外部文件。对象名称为emss-78112-f006.jpg
FAM35A或C20或196缺失可恢复BRCA1缺陷细胞的HR。

a、,BRCA1或BRCA2耗竭后(5Gy后5h)U2OS细胞中GFP-FAM35A(左侧面板)和GFP-C20或196(右侧面板)IRIF的定量。条形表示平均值±SEM,单向Anova;n=3个独立实验,FAM35A siCTRL(n=4)、FAM35AsiBRCA2(n=2)和C20或196 siCTRR(n=5)除外;在条形图上绘制单个数据点。b、,在有或无BRCA1缺失的U2OS细胞中对53BP1和可诱导GFP-FAM35A IRIF进行定量(5Gy,指示时间点)。条形代表平均值±SEM,单向方差;n=4个独立实验,除53BP1 1.5h siCTRL(n=2)、53BP1 1.5 h siBRCA1和53BP1 16h siCTRR(n=3)、FAM35A 1.5h siCTRL(n=5)外;在条形图上绘制单个数据点。c、,如图所示,在Cyclin A(CycA)阳性RPE1ko细胞系中RAD51 IRIF(5Gy后5.5小时)的代表性图像(左面板)和量化(右面板)。条形代表平均值±SEM,单向方差;n=3个独立实验,在条形图上绘制单个数据点。比例尺10µm。日期:,FAM35A/C20orf196在激光诱导DNA损伤部位2小时后恢复BRCA2募集巴西航空公司1-空单元格(有关量化,请参阅补充图6e). 比例尺10µm。e(电子),用指定的siRNAs处理的U2OS-TLR细胞的HR检测(有关门控策略,请参见补充图6f). 条形代表平均值±SEM,单向方差;n=4个独立实验,在条形图上绘制单个数据点。f、,FAM35A/C20orf196失活可减轻BRCA1ko细胞自发染色体畸变的形成。图中显示了中期扩散的代表性图像,并进行了量化;条形图表示平均值,n=2个独立实验,FAM35Ako和C20或196ko(n=1)除外,在条形图上绘制单个数据点。用指示的RPE1ko和补体细胞系进行Olaparib克隆存活试验。条形代表平均值±SEM,单向方差;n=4个独立实验,第4组和第5组(n=3)和第3组(n=2)除外;AUC如所示补充图6g对于a-c和e,*p<0.05,**p<0.01,***p<0.001,****p<0.0001,ns=不显著(p≥0.05);包含精确p值的统计源数据可以在补充表5.

总之,我们的结果表明,FAM35A/C20orf196作用于BRCA1和53BP1/RIF1/MAD2L2相反功能之间的接口,以调节DSB修复途径的选择。虽然这种作用至少可以部分通过控制DSB切除来实现,但DSB切除在BRCA1缺陷细胞中调节不当且动力学较慢11,36,我们推断FAM35A/C20orf196也可能导致BRCA1缺陷细胞DNA损伤部位BRCA2介导的RAD51负载的严重缺陷11,35,37事实上,对于53BP1失活,FAM35A或C20orf196的缺失恢复了RAD51 IRIF的形成巴西航空公司1-空单元格(图6c). 在探索这种效应的机制时,我们发现,在没有BRCA1的情况下,通过DNA损伤部位的RPA募集来测量,FAM35A和C20orf196敲除细胞的切除水平仍然较高(补充图6c-d). 此外,FAM35A/C20orf196的缺失也减轻了BRCA1缺陷细胞中BRCA2向DNA损伤位点募集的受损(图6d,补充图6e). 因此,研究基于细胞的染色体红绿灯报告(TLR)HR系统38,39确定BRCA1缺陷细胞中FAM35A或C20orf196的缺失使HR恢复到与在此情况下53BP1缺失时获得的HR相似的水平(图6e,补充6f). 此外,去除FAM35A或C20orf196可挽救BRCA1敲除细胞的自发基因组不稳定表型(图6f). 基于我们的发现,FAM35A N末端区域在很大程度上调节了其对IRIF的定位(图2f,补充图3e-f)引入FAM35A N端而非C端,增强了BRCA1/FAM35A-空细胞的PARP抑制剂敏感性(图6g,补充图6g). 此外,在BRCA1基因敲除细胞中,FAM35A的失活与53BP1的失活相关,这与PARP抑制剂的耐药性有关(补充图6h). 综合考虑我们的研究结果,我们建议将C20orf196和FAM35A分别命名为SHLD1和SHLD2,或统称为“Shieldin复合体”,因为它可以保护DSB免受不当活动的影响,并促进适当的DSB修复模式。

FAM35A/C20orf196缺失与癌症中PARP抑制剂耐药性相关

已确定SHLD1C20或f196和SHLD2FAM35A系列在调节BRCA1缺陷乳腺癌细胞系的PARP抑制剂敏感性时,我们推测这可能也适用于更多的生理环境。因此,我们采用了BRCA1缺陷型乳腺癌患者衍生异种移植(PDX)模型,该模型在存在或不存在奥拉帕林(队列2)的小鼠体内繁殖(图7a). 耐药肿瘤进一步连续传至新的宿主,这些宿主在有或无olaparib的情况下接受治疗,以确认并维持耐药性(图7a; 另请参阅补充图7a). 然后采集肿瘤并进行全基因组RNA测序。值得注意的是,我们的分析表明,与其他队列相比,几乎所有慢性奥拉帕林治疗后的耐药肿瘤(队列4)都与SHLD1 mRNA表达降低相关C20或196、SHLD2FAM35A系列、53BP1和/或PARP1(图7a; 每个heatmap列表示一个肿瘤/小鼠样本)。因为这个肿瘤模型是多克隆的40,我们的数据表明,奥拉帕林抵抗机制可能是通过平行的进化轨迹收敛于屏蔽蛋白活性的丧失而产生的。此外,当我们通过SHLD1/2的表达对BRCA1-deficient PDX肿瘤队列进行分层时,随后的分析表明低SHLD1C20或196转录水平与olaparib内在耐药相关(图7b). 其中一个抗奥拉帕林模型(PDX127)显示伴随的SHLD1丢失C20或196和SHLD2FAM35A系列而其他两种SHLD1/2转录水平正常的耐药模型则含有有害的53BP1突变。值得注意的是,其中一些抗药性PDX模型也显示BRCA1核病灶41这表明由于肿瘤异质性和/或机制合作,存在多种耐药机制。

保存图片、插图等的外部文件。对象名称为emss-78112-f007.jpg
FAM35A或C20orf196缺失与癌症中PARP抑制剂耐药性相关。

a、,体内PDX研究示意图(顶面板)。mRNA测序生成的热图显示了相应PDX样本中指示基因的缩放表达水平(下面板);第1-4组分别为6、5、7、8只小鼠。b条,C20orf196/FAM35A在来自BRCA1缺陷肿瘤的乳腺癌和卵巢癌PDX中的表达。y轴:每百万分之log2转录本。线代表平均值±SEM;对于SHLD1高、SHLD1低、SHLD2高、SHLD2低组,n分别为12、4、15、1;双尾非配对学生t检验***p=0.0003。可在中找到PDX的统计源数据补充表5方法.c-d、,IR(c)或顺铂治疗(d)后显示的具有AUC的RPE1ko细胞系的克隆存活率测定补充图7b和7c分别为。所示数据代表平均值±SEM(n=3个独立实验,但c组7和d组7除外,其中n=2)e(电子)FAM35A/C20orf196缺失导致顺铂诱导的FANCD2病灶增加。条形代表平均值±SEM,单向方差*p<0.05,**p<0.01,***p<0.001,***p<0.05,ns=不显著(p≥0.05);n=4个独立实验,在条形图上绘制单个数据点;统计源数据可以在中找到补充表5.比例尺10µm。f、,在存在或不存在功能性BRCA1的情况下,SHLD1/2在DSB修复中的作用的拟议模型。

最后,我们发现与53BP1缺乏相比42、SHLD1C20或f196或SHLD2FAM35A系列缺失增加了BRCA1-profective和BRCA1-null细胞对IR的敏感性,更显著地提高了它们对DNA交联剂顺铂的敏感性(图7c-d,补充图7b-c). 此外,增强SHLD1对顺铂的敏感性C20或196或SHLD2FAM35A系列失活与参与检测和修复DNA交联的FANCD2蛋白增加DNA损伤焦点形成有关(图7e,补充图7d). 因此,这些发现表明,如果患者中SHLD1/2的表达缺失/减少,可能会提供临床上可以利用的副作用治疗漏洞。

讨论

在过去二十年中,真核细胞已经进化出多种DNA DSB修复机制,这些机制以复杂复杂的方式进行调节,以优化基因组稳定性。特别是,人们关注的焦点是细胞如何战略性地使用DSB修复的两种主要模式——NHEJ和HR,这两种模式相互对抗,在不同的环境中以最佳方式运行,其相对使用受染色质结构和细胞周期阶段等因素的规定。除了学术兴趣外,这类课题的研究也具有临床意义,特别是在癌症治疗中,DSB诱导的化疗药物被频繁使用,而PARP抑制剂等分子靶向药物正越来越多地用于特定环境。患有以下疾病的患者体内固有或产生的PARP抑制剂耐药性巴西1/2然而,突变是一个日益严重的临床问题。使用全基因组CRISPR-Cas9合成活性/耐药性筛选,我们发现了两种以前没有特征的蛋白质——SHLD1C20或196和SHLD2FAM35A系列-其缺失介导PARP抑制剂耐药性,我们已经证明其是53BP1/RIF1/MAD2L2分子轴的最远端因子,在BRCA1缺乏的环境中促进NHEJ并限制HR。我们的筛选还确定了其他候选PARP抑制剂耐药因子,等待未来研究的验证。

机械方面,我们已经证明SHLD1C20或196和SHLD2FAM35A系列与SHLD1形成复合物,称为ShieldinC20或196通过与SHLD2的相互作用招募DNA损伤位点FAM35A系列和其他因素,以及SHLD2FAM35A系列通过其C-末端OB折叠区与单链DNA相互作用。此外,我们已经确定SHLD1C20或196和SHLD2FAM35A系列促进NHEJ的方式可能通过其对限制DNA末端切除的作用来介导,并通过拮抗BRCA2/RAD51在切除的ssDNA上用RPA替代BRCA1而成为HR的屏障。我们的工作与最近的报告一致43独立识别SHLD1C20或196和SHLD2FAM35A系列在BRCA1缺陷细胞中作为NHEJ的传播因子和HR的拮抗剂。本研究还确定了第三种成分,RINN1/SHLD3CTC-534A2.2型这被认为是RIF1到MAD2L2和SHLD1/2的分子桥梁。

虽然Shieldin的丢失可能通过影响切除和BRCA2/RAD51负荷来促进BRCA1缺陷细胞的HR恢复,但这些机制的相对重要性需要进一步研究。我们注意到,在没有Shieldin的情况下进行更广泛和更快的切除可能会增加BRCA2/RAD51的负荷。或者,或者另外,在BRCA1缺陷细胞中,Shieldin可以作为一个物理屏障,阻止BRCA2/RAD51在dsDNA/ssDNA连接处加载,可能是通过53BP1/RIF1/MAD2L2复合物将其连接到DSB侧翼染色质,从而使远端SHLD2的C末端相互作用FAM35A系列和ssDNA(参见图7f对于模型,以及补充图7e对于SHLD2家庭35a域功能摘要)。然而,我们发现SHLD2的过度表达FAM35A系列N-而非C-末端赋予奥拉帕林超敏反应BRCA1/FAM35A双敲除细胞,表明至少在这种情况下,SHLD2结合染色质FAM35A系列在限制人力资源方面起主导作用。相比之下,我们发现SHLD2FAM35A系列N端和C端对IR抗性很重要(在BRCA1-profective细胞中)。由于Shieldin缺陷细胞的IR敏感性可能反映NHEJ受损,我们推测Shieldin通过限制DSB切除以及与其他NHEJ生成因子组装以将DSB末端连接在一起以促进其并置和修复来增强NHEJ。

值得注意的是,我们发现SHLD1C20或196或SHLD2FAM35A系列失活增强了顺铂对BRCA1-null或BRCA1-profective细胞的敏感性。这种敏感性可能不能反映Shieldin在促进NHEJ、限制DNA末端切除或拮抗BRCA1介导的BRCA2/RAD51负载中的作用,因为在我们手中,53BP1缺失对顺铂敏感性没有显著影响。顺铂和其他化合物产生的序列内DNA交联(ICL)由范科尼贫血(FA)途径检测和修复,FA的关键蛋白是FANCD2,在这些损伤部位形成病灶44我们观察到,顺铂治疗后,FANCD2病灶在SHLD1的细胞中更明显C20或196或SHLD2FAM35A系列已停用。因此,我们有兴趣确定屏蔽蛋白样MAD2L2(与REV3L一起)是否是反义DNA合成(TLS)聚合酶Pol zeta的调节亚单位,其双等位基因失活是否与FA有关4547–还可以通过TLS机制促进ICL修复。

最后,根据我们的研究结果,评估患者肿瘤活检中SHLD1/2的表达,确定该信息是否可用于PARP抑制剂治疗的患者分层,将是一件有趣的事情,并确定BRCA1缺陷癌症患者在PARP抑制剂治疗后出现耐药性时是否出现SHLD1/2表达变化。在这方面,我们注意到,如果在BRCA1缺陷的癌症中下调Shieldin组分确实会产生临床耐药性,这可能允许替代治疗,例如基于铂化合物的治疗。

方法

CRISPR-Cas9屏幕

使用基因组规模(GeCKO)2.0版执行20SUM149PT细胞在感染倍数(MOI)为0.3倍和250倍的文库覆盖率下转导。然后用嘌呤霉素筛选细胞7天,然后用3种不同的PARPi再处理14天。使用的IC有:;Olaparib IC95-2μM,BMN673 IC95-5 nM,AZD2461 IC70-4µM。收集每个条件下的存活克隆,分离基因组DNA(gDNA)(血液和细胞培养DNA Midi Kit,Qiagen),并使用Illumina兼容引物进行PCR,然后进行Illuminia测序。通过软件包MAGeCK 0.5.5版(参见“代码可用性”).

细胞培养

将U2OS、U2OS衍生、HEK293、HEK2903T-LentiX细胞培养为39.RPE1 p53缺失FRT48在添加17 ml NaHCO的F-12(Ham’s F-12;Sigma-Aldrich)中培养RPE1 p53缺失FRT衍生细胞7.5%/500ml(Sigma-Aldrich)。所有培养基均添加10%(v/v)胎牛血清(FBS;BioSera)、100 U/ml青霉素、100µg/ml链霉素(Sigma-Aldrich)和2 mM L-谷氨酰胺。SUM149PT细胞在添加5%(v/v)FBS(BioSera)、10 mM HEPES、1μg/ml氢化可的松、5μg/ml胰岛素、抗生素的火腿F12营养混合物(ThermoFisher)中培养。为了维持和选择稳定表达GFP或GFP标记构建物的RPE1 FRT p53无效或U20S Trex细胞,使用2μg/ml速溶素(Sigma-Aldrich)和0.5 mg/ml G418(Invitrogen)。U2OS-TLR按39除RPE1 p53缺失FRT衍生细胞外,稳定表达诱导构建物的U20S Trex细胞与1µg/ml多西环素(Sigma-Aldrich)培养24-48h以诱导GFP构建物的表达。所有细胞最初都是从ATCC细胞库中获得的,并且常规测试为无支原体。U2OS和RPE1细胞系最近通过Affymetrix SNP6拷贝数分析进行了鉴定。Trf2-/-;p53-/-;如前所述的TRF2(Ile468Ala)MEF(TRF2ts MEF)33,49.CH12F3(CH12)29和CH12-Cas9细胞系在补充有10%FBS、100U/ml青霉素、100μg/ml链霉素、50µM 2-巯基乙醇、1xMEM非必需氨基酸、1mM丙酮酸钠和10mM HEPES的RPMI 1640中培养。

人类稳定细胞系的产生和敲除

通过转染含有GFP标记结构和pOG44(分别为1:4)的pcDNA5/FRT/TO-neo,产生稳定表达可诱导GFP标记构建物的U2OS Trex或RPE1 p53缺失FRT衍生细胞。在48小时时使用0.5mg/ml G418(Invitrogen)开始选择。通过转染“All-in-one”质粒在RPE1 p53阴性细胞中产生敲除48使用MoFlo(Beckman Coulter)通过GFP表达进行单细胞排序。扩增单个克隆,提取基因组DNA并通过PCR、TOPO克隆和测序进行筛选。Fam35a和C20orf196的已验证小鼠胚胎干细胞(mESC)敲除物来自Haplobank(www.haplobank.at网站).

质粒和克隆

请参见补充表3.

siRNA与质粒转染

siRNAs从MWG或IDT中获得,并根据制造商的协议使用Lipofectamine RNAiMAX(Invitrogen)转染。根据制造商的方案,使用TransIT-LT1(Mirus Bio)进行质粒转染。对于siRNA和DNA联合转染,质粒在siRNA处理后8小时转染。请参见补充表4.

随机质粒整合试验

按照前面的描述执行27.

使用化学试剂、电离辐射和激光微辐射诱导DNA损伤

按照前面的描述执行39.

FRAP与缔合动力学

按照前面的描述执行50.

TLR分析

交通灯报告器(TLR)分析和此处使用的构造已在前面详细描述38,39.

细胞周期分析

按照前面的描述执行39.

克隆原性存活试验

按照前面的描述执行27,39.

全细胞提取物和免疫印迹

按照之前的描述执行39为了检测磷酸-RPA(pS4/S8-RPA2),通过在2xSDS缓冲液中刮取细胞,然后使用4-12%双三凝胶(Invitrogen)进行SDS-PAGE,制备裂解物,并使用SuperSignal West Pico PLUS(Thermo Scientific)进行免疫印迹。使用IRDye800CW和IRDye680标记的二级抗体在奥德赛红外成像仪(LI-COR)上进行检测。使用ImageJ对斑点进行量化。所有蛋白质浓度均采用BCA分析(Pierce)测定。所有抗体都列在补充表2.

免疫沉淀

如前所述,所有免疫沉淀程序执行两次39。如图所示的联合免疫沉淀图2d(FAM35A),293T细胞与pMSCV-blas-eGFP-MAD2L2和pLX304-blast-V5-Empty或pLX304-V5-FAM35A。转染后72h的细胞暴露于25Gy IR,然后恢复3h。使用GFP-Trap_MA珠(ChromoTek),并根据制造商的方案进行免疫沉淀。对于C20或196(补充图2c),293T细胞用pMSCV-blas-3xFlag-hMAD2L2和pcDNA5.1-GFP或pcDNA5.1-GFP-C20或f196共转染。转染后72小时,细胞暴露于25Gy IR,然后恢复3小时。用冷PBS洗涤后,在补充有上述相同抑制剂的1ml裂解缓冲液(50mM Tris HCl pH7.4;150mM NaCl;1mM EDTA;1%Triton X-100)中裂解细胞。在冰上培养30分钟,然后离心(16000 g),将预先用TBS(50mM Tris HCl,150 mM NaCl pH7.4)清洗过的反标记M2磁珠(M8823,Sigma-Aldrich)添加到裂解液中,并在4°C下旋转过夜。免疫复合物经5min煮沸洗脱。

DNA下拉实验

有关程序的详细说明,请参见39使用具有以下序列的寡核苷酸:

5'BiosG/ATCGCATTGGCATTGGAATGCAGACGACTGATCGAGCTGATTCGATCGCTGATCCGATCGTATCCGTACATATCGAT-3'(IDT)

体外GST下拉

用冰镇PBS清洗葡硫酮琼脂糖珠(GE Healthcare),并用添加10%细菌裂解物(非诱导BL21细胞,用PBS/溶菌酶裂解)的PBS封闭30分钟,然后在结合缓冲液(10mM Tris pH7.5,150 mM NaCl,0.5%NP40,0.5 mM EDTA,0.5%BSA)中重新悬浮。将纯化的GST(细菌表达)、GST-FAM35A(Novus Biologicals)和His-C20orf196(Creative BioMart)以2 pmol添加到珠子中,并在4°C下培养30分钟。用10 mM Tris,pH 7.5,250 mM NaCl,0.5%NP40,0.5 mM EDTA洗涤珠子5次,并用100 mM Tris-pH 8,20 mM还原谷胱甘肽,120 mM NaCl在4°C下旋转15分钟洗脱。将洗脱液煮沸5分钟,装入4-12%Bis-Tris凝胶(Invitrogen)并进行western印迹。用指示的抗体探测斑点。

重组蛋白纯化和电泳迁移率变化分析(EMSA)

使用相同的方法纯化野生型和突变型FAM35A C末端结构域。收集的细胞在50 mM Tris pH 8.0、5%甘油、150 mM NaCl、2 mMβ-巯基乙醇、10mM咪唑、蛋白酶抑制剂(罗氏)和40 g/ml脱氧核糖核酸酶I(西格玛)中通过超声分解。在30000g离心30分钟后,将上清液加载到含有Ni-NTA亲和树脂(Qiagen)的重力柱上,该树脂与50 mM Tris pH 8.0、5%甘油、150 mM NaCl、2 mMβ-巯基乙醇和10 mM咪唑预平衡。在用相同的缓冲液洗涤小球10倍柱体积后,使用50 mM Tris pH 8.0、5%甘油、150 mM NaCl、2 mMβ-巯基乙醇和100 mM咪唑洗脱蛋白质。用Q柱(GE healthcare)缓冲液A(20 mM Tris pH 8.0,50 mM NaCl,5%甘油和2 mMβ-巯基乙醇)透析洗脱液,并加载到5 ml Q柱上。在缓冲液B(20mM Tris pH 8.0,1M NaCl,5%甘油和2mMβ-巯基乙醇)的梯度中洗脱蛋白质。收集含有FAM35A蛋白的组分,并通过Superdex 200 10/300色谱柱(GE Healthcare)进一步纯化,该色谱柱在缓冲液GF(20 mM Tris pH 8,150 mM NaCl,5%甘油和5 mM DTT)中平衡。纯化每个步骤中的蛋白质样品在4-12%双三凝胶(Invitrogen)上进行分析。将最终纯化样品浓缩并储存在-80°C下。正向和反向90-bases DNA寡核苷酸(IDT)(F:6-FAM(6-羧基荧光素)-ATCGCATTGGCATTGCAATGCGATACGACTGAGCTAGCTGAGCTGATTCCGGCTTATTCCGGTACATATCGAT);R: 将6-FAM-ATCCGATGTATGTAGTACAGGAATAGCCGGATATCAGCTAGTGTACCTCGATCAGTGTATCCATCATCATCATGCAATGCAT溶解在退火缓冲液(10 mM Tris pH 8.0、50 mM NaCl和1 mM EDTA)中,使最终浓度达到100μM。DNA寡核苷酸F被用作EMSA的ssDNA。混合等体积的DNA寡核苷酸F和R并退火(加热至95°C 2分钟,冷却至25°C超过45分钟)以生成dsDNA。每20μl EMSA反应包含10 nM的ssDNA/dsDNA,与不同浓度的蛋白质在20 mM Tris-HCl pH 7.5、50 mM KCl、5%(v/v)甘油、100μM DTT、10μg/ml BSA中孵育。样品在37°C下孵育15分钟,然后在0.5xTBE缓冲液中的5%聚丙烯酰胺天然凝胶上进行4°C电泳。台风9000(GE Healthcare)对DNA进行了可视化。

用于质谱分析的GFP-Trap下拉

HEK293T细胞在含有L-精氨酸和L-赖氨酸或L-精氨碱的SILAC培养基中培养[13C6,15N4]和L-赖氨酸[13C6,15N个2](剑桥同位素实验室)如前所述51转染后48小时,在补充有蛋白酶、磷酸酶抑制剂和N-乙基马来酰亚胺的改良RIPA缓冲液(50mM Tris pH 7.5、150mM NaCl、1mM EDTA、1%NP-40、0.1%脱氧胆酸钠)中裂解细胞。通过在4°C下16000×g离心15分钟来清除裂解液,并使用QuickStart-Bradford蛋白质分析(BioRad)估计蛋白质浓度。根据SILAC条件,将20µl预先平衡的GFP-Trap-A珠(ChromoTek)添加到2 mg裂解液中,并在4°C下旋转培养1小时,然后用改良的RIPA缓冲液进行3次洗涤。将结合的蛋白质在补充有1mM二硫代苏糖醇的NuPAGE®LDS样品缓冲液(Life Technologies)中洗脱,在70ºC下加热10分钟,并在室温下用5.5mM氯乙酰胺烷基化。将样品加载到4-12%梯度SDS-PAGE凝胶上,使用胶体蓝染色试剂盒(Life Technologies)对蛋白质进行染色用胰蛋白酶在凝胶中消化。从凝胶中提取肽并在反相C18 StageTips上脱盐52.

质谱分析

如前所述,在配备EASY-nLC 1000(Thermo Scientific)的四极轨道质谱仪(Q Exactive Plus,Thermo科学)上分析肽组分53将肽样品加载到C18反相色谱柱上,用含0.1%甲酸的8~40%乙腈线性梯度洗脱2小时。质谱仪以数据相关模式运行,在MS和MS2采集之间自动切换。Orbitrap中获得了测量全扫描MS光谱(m/z 300–1650)。十种最强烈的离子被HCD依次分离和破碎54在Orbitrap质量分析仪中获得碎片光谱。使用MaxQuant(开发版本1.5.2.8)分析原始数据文件55。使用Andromeda搜索引擎,在包含92578个人类蛋白质序列的数据库中搜索母离子和MS2光谱,该序列来自2016年12月发布的UniProtKB56在MS模式下以6 ppm的质量容限、HCD MS2模式下以20 ppm的质量容限、严格的胰蛋白酶特异性和最多允许三次误滤进行光谱搜索。半胱氨酸氨基甲酰化被搜索为固定修饰,而蛋白质N-末端乙酰化、蛋氨酸氧化、半胱氨酸的N-乙基马来酰亚胺修饰被搜索为可变修饰。根据后验错误概率对数据集进行过滤,以获得低于1%的错误发现率,该错误发现率使用目标恢复方法进行估计57.

免疫荧光和显微成像

对γH2AX、RAD51、RPA、ssDNA(BrdU)、BLM、BRCA2、FANCD2、Cyclin A和GFP(FAM35A和C20orf196)进行共焦成像,如39对于RAD51和Cyclin A,省略了预提取步骤,固定后将细胞在PBS中0.2%Triton X-100(Sigma)中渗透15分钟。使用Deltavision OMX 3D-SIM System V3 BLAZE(GE Healthcare公司的Applied Precision)采集超分辨率图像,该系统配备3台sCMOS摄像机、405、488、592.5nm二极管激光照明、Olympus Plan Apo N 60x 1.42NA油物镜以及标准激发和发射滤波器组。使用所述照明模式的三个角度和五个相移,依次对每个通道进行成像58。以0.125μm z步长采集截面。在SoftWoRx软件版本6.5.2(Applied Precision,GE Healthcare公司)中重建原始OMX数据并注册通道。使用斐济的自定义脚本(Butler R)测量GFP-FAM35A信号相对于53BP1信号的Voxelwise近邻距离(https://github.com/gurdon-institute/OMX-空间分析). 该脚本使用Kapur的最大熵阈值法映射信号量59和使用精确的签名3D欧氏距离变换测量距离,内部距离设置为零,以便在直方图上显示。对于所有图像,比例尺=10μm。

多重荧光原位杂交(M-FISH)

随后进行人24色多重FISH(M-FISH)探针制备和载玻片处理60对于每个人类细胞样本,根据M-FISH分类和DAPI带型对10-30个中期进行核型分析。如前所述,使用BAC探针对中期传播进行FISH61并手动计数。对于类开关重组(CSR)分析,如前所述在中期扩散上进行DNA FISH61并手动计数。每个基因型至少评估470个中期,每个条件至少使用2个独立克隆。对于使用Alexa488-(CCCTAA重复序列)肽核酸定制探针(PN-TC060-005,Panagene/Eurogentec)进行端粒去盖、细胞采集、中期扩散的制备和端粒FISH,如前所述进行中期染色体分析6这些数据代表2个独立的实验,每个条件下≥1300条染色体,在盲法基因型后手动计数。

端粒融合分析:MEFs病毒转导

细胞像以前一样转导49使用从MISSION shRNA文库(Sigma)获得的pLKO-puro shRNA慢病毒,对抗描述的或未翻译的小鼠基因(“空”)。端粒NHEJ的评估.TRF2ts MEF在39°C的非许可温度下生长24小时,以灭活TRF2并因端粒去盖而诱导NHEJ依赖的染色体端到端融合。

类别开关重组分析:CH12和CH12-Cas9细胞系

CH12-Cas9细胞系是通过转导CH12-Cas9细胞而产生的,CH12-Cas8细胞是用HEK293T中包装的慢病毒颗粒进行自旋感染制成的。质粒:pKLV2-EF1aBsd2Acas9-W、pxPAX2(添加基因#12260)、VSV-G和pMD2.G(添加基因#12259)。转导一周后48小时开始选择Blastidin(10μg/ml)。用于类开关重组分析的sgRNA表达质粒。gRNA被用于靶向Fam35a、C20orf196和Mad2l2/Rev7小鼠基因(每个靶基因2个sgRNAs,序列列于补充表4将.sgRNAs克隆到pKLV翻转的U6gRNA_CCDB_PB_BbsI_PGKpuro2ABFP载体中62.生成野生型和敲除型CH12-Cas9细胞克隆。53bp1空CH12细胞克隆(来自Fred Alt的礼物)如前所述63用2.2μg的sgRNA-1和sgRNA2以及0.6μg的piggyBac转座酶表达载体对1200万CH12-Cas9细胞进行核感染64,使用Amaxa Nucleofector、Nucleofector®Kit V解决方案(Lonza)和程序X-001。两天后,用3μg/ml嘌呤霉素筛选BFP阳性/耐嘌呤菌素的CH12-Cas9细胞,为期一周。然后将细胞单细胞稀释成96个平板,进一步培养并通过PCR进行筛选,并使用中列出的PCR引物进行Sanger测序补充表4.分类开关重组和细胞增殖分析。将CH12细胞以50000个细胞/ml的速度接种在添加了抗CD40抗体(1μg/ml,Miltenyi)、IL-4(20 ng/ml,Millenyi)和TGF-β(1 ng/ml)的完整RPMI中,以诱导IgM向IgA转换。3天后,使用IgA-PE抗体(eBiosciences)和Canto II分析仪(BD Bioscinces)通过流式细胞术分析细胞的类别转换。使用Casy细胞计数器(Roche)对活细胞进行计数。在三个独立的实验中,对3xwild-type(WT)、3xFam35a敲除(Fam35a)、3x20orf196敲除(C20或f196)、2x53bp1敲除(53bp1)和3xMad2l2敲除(Mad2l 2)进行了CSR和增殖测定。RT-PCR分析。 伊格,α种系转录本(αGLT)和援助如前所述,对mRNA进行定量12底漆列于补充表4.

患者源性肿瘤异种移植

如前所述,PDX是从同意的乳腺癌或卵巢癌患者样本中生成和建立的40该研究得到了剑桥郡2 REC国家研究伦理服务局(REC参考号:08/H0308/178)和瓦尔德赫布伦医院临床调查伦理委员会(PR(AG)183/2012)的适当批准。STG201是本研究中使用的PDX模型,是一种BRCA-null模型,以BRCA1启动子甲基化、BRCA1 mRNA和蛋白表达缺失为特征。我们之前已经证明其在体内和PDX衍生细胞中对PARP抑制剂(包括奥拉帕林)的敏感性。STG201也与深层分子和药物敏感性注释有关40http://caldaslab.cruk.cam.ac.uk/bcape/所有其他PDX均来自乳腺或卵巢肿瘤巴西航空公司1-突变携带者或BRCA1型启动子超甲基化引起的表观遗传沉默41.PDX127未显示BRCA1的任何联合表达,但在FAM35A和C20orf196表达中均较低。5个PARPi敏感的PDX中没有一个表现出低水平的C20orf196、FAM35A或53BP1缺失,也没有BRCA1亚型。该研究符合所有涉及动物使用和人类参与者的研究相关道德规范。

体内获得性耐药性的产生

如前所述,随机化后(50mg/kg,5天/周),将AZD2281(Olaparib/Lynparza)作为PARP抑制剂(50mg/kg,5IW)给予免疫功能低下的荷瘤小鼠40为了对皮下植入物的反应进行分类,我们修改了RECIST标准,以持续奥拉帕林治疗后肿瘤体积变化百分比为基础:完全反应(CR),最佳反应≤-95%;部分反应(PR),-95%<最佳反应≤-30%;稳定期疾病(SD),-30%<最佳反应≤+20%;进行性疾病(PD),治疗第21天肿瘤体积变化百分比>+20%。PARPi耐药PDX表现为PD,而PARPi敏感模型表现为SD、PR或CR。对于STG201,在治疗25天后收集时间匹配的溶媒和奥拉帕林治疗的样品(PARPi初始PDX),并进行RNA提取和测序处理。研究中的几只小鼠继续接触奥拉帕林,直到肿瘤再生。其中一个耐药肿瘤在治疗126天后连续传至新的宿主小鼠(PARPi耐药PDX),并用进一步的载体或olaparib进行治疗。治疗58天后,确认耐药表型,并收集样本并进行RNA测序,如下所述。生长曲线显示每个试验组至少5个独立肿瘤体积的平均和标准偏差。所有实验程序均由剑桥大学动物福利和伦理审查委员会以及瓦尔德赫布伦医院临床调查伦理委员会和动物使用委员会批准。RNA测序。按照制造商的说明,使用Qiagen-miRNeasy或RNeasy-Mini试剂盒(Cat ID,217004或74104)从所有样品中提取RNA。使用TruSeq Stranded mRNA HT试剂盒或Total RNA Library Prep试剂盒和Ribo-Zero Gold(Cat ID,RS-122-2103或RS-122-2301,Illumina)制备Illuminia测序文库。500克RNA完整性数(RIN)大于8的总RNA用于文库制备。样品按照制造商的HS(高样品)说明(零件号15031048 Rev.E,Illumina)进行处理,在DNA片段富集步骤中使用12或15个PCR循环。使用KAPA文库定量试剂盒Illumina ROX Low(Cat ID,KK4873,KAPA Biosystems)对所有文库进行定量并进行归一化。按照制造商的指示,将库集中在等体积中,并将池用于HiSeq4000测序流动池上的聚类。对双诱导文库使用150bp或100bp配对序列运行类型进行测序。在比对之前,读取的测序质量是使用Trim-Galore!(v0.4.2)http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/然后,如Callari等人或Ahdesmaki等人所述65使用STAR(v2.5.2b)将读取结果与人类(hg19)和小鼠(mm10)的联合参考基因组对齐66,67使用featureCounts(v1.5.2)将计数分配给基因组特征,从而使用比对分数准确区分来自人类和小鼠的读数68。对多个测序运行的计数进行合并,然后使用edgeR包进行归一化69,70.

代码可用性

OMX分析中使用的自定义FIJI脚本可以在以下位置找到https://github.com/gurdon-institute/OMX-Spatial-Analysis网站.

用于CRISPR-Cas9屏幕的MAGeCK命令为:

mageck test-k counts.csv-c DMSO-t WC_2461-n WC_2461

mageck test-k counts.csv-c DMSO-t WC-673-n WC-673

mageck test-k counts.csv-c DMSO-t WC-2281-n WC-2281

统计和再现性

除非另有说明,否则Prism v7.0b(GraphPad软件)用于生成图形、执行统计测试和计算p值。误差条、统计测试和独立重复次数(n)用图表图例表示,统计源数据包括中提供的精确p值补充表5统计检验包括双尾Student t检验、Fisher’s Exact检验和单向方差分析(ANOVA),后者均按照多重比较的建议进行校正。使用ImageJ/FIJI或Volocity 6.3(Perkin Elmer)进行显微镜图像分析。每个药物治疗均使用三个克隆进行CRISPR筛查。GFP-FAM35A和GFP-C20orf196的质谱在两个独立的实验中进行。由于跨多个通道进行测序(在进一步分析之前进行合并),每个试验组的RNA测序为三个重复。对以下数量的独立生物样本进行了研究:队列1中的6个PDX、队列2中的5个PDX,队列3中的7个PDX和队列4中的8个PDX。对于SHLD1高表达和低表达队列,分别评估了12个和4个独立PDX。除非另有规定,否则所有免疫荧光分析定量数据表示3个独立生物重复的平均值±SEM,每个条件下n≥30个细胞。所有免疫印迹都是两个独立实验的代表,免疫印迹的未处理扫描如补充图8.

补充材料

报告摘要

单击此处查看。(100K,pdf格式)

补充图1-8

单击此处查看。(1500万,文件)

补充表1

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

补充表2

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

补充表3

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

补充表4

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

补充表5

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

致谢

我们感谢所有SPJ实验室成员的支持和建议,感谢剑桥大学的同事N.Lawrence对OMX超分辨率显微镜的支持,感谢R.Butler对计算图像分析和编程的帮助。我们感谢S.Selivanova和S.Hough在质粒扩增、样品制备和组织培养维护方面的帮助,K.Dry的广泛编辑协助,F.Muñoz Martinez在CRISPR-Cas9敲除生成方面的协助,L.Radu在蛋白质纯化方面的协助,C.Lord(伦敦癌症研究所)针对SUM149PT细胞,D.Durocher(加拿大多伦多大学)针对U2OS LacSceIII细胞,F.Alt(美国哈佛大学)针对CH12F3细胞和P53结合蛋白1淘汰CH12F3细胞克隆,T.Honjo(日本京都大学)获得使用CH12F4细胞系的许可,以及雅各布斯实验室的J.Serrat获得技术援助。SPJ实验室主要由英国癌症研究(CRUK)计划拨款C6/A18796和威康信托(WT)研究员奖206388/Z/17/Z资助。核心基础设施资金由CRUK拨款C6946/A24843和WT拨款WT203144提供。SPJ从剑桥大学领取薪水。HD由WT Clinical Fellowship 206721/Z/17/Z资助。TWC由剑桥国际奖学金资助。D.P.由英国癌症研究生奖学金C6/A21454资助。PB实验室由德国研究基金会的Emmy Noether计划(BE 5342/1-1)和欧洲委员会的Marie Curie职业整合拨款(630763)支持。LP实验室由WT(研究员奖104641/Z/14/Z)和医学研究委员会(项目拨款MR/N000161/1)资助。CC实验室得到了CRUK的资助。JJ实验室得到了欧洲研究委员会拨款ERC-StG 311565、荷兰癌症协会拨款KWF 10999和荷兰科学研究组织(NWO)的支持,作为荷兰国家路线图大型研究设施的一部分,蛋白质@工作(荷兰癌症研究所蛋白质组学设施项目编号184.032.201)。LD实验室由巴斯德研究所、国家癌症研究所(#PLBIO16-181)和欧洲癌症研究委员会资助(启动拨款协议#310117)。WW是巴斯德-巴黎大学(PPU)国际博士项目的一部分,该项目得到了中国CNBG公司的资助。QW由Wellcome信托基金资助(200814/Z/16/Z to T.L.B)。VS实验室工作由西班牙经济和创新部的Salud Carlos III研究所(ISCIII)资助,该研究所部分由欧洲区域发展联邦储备基金(PI17-01080 to VS)、欧洲研究区域网、Transcan-2(AC15/00063)、与阿斯利康英国有限公司的非商业研究协议支持,以及来自Agència de Gestiód’Ajuts Universitaris i de Recerca(AGAUR,2017 SGR 540)和Orozco家族的结构性资金。VS获得了ISCIII(CP14/00228,MV15/00041)和FERO基金会对CC实验室的工资和差旅支持。

脚注

数据可用性

SUM149细胞中全基因组CRISPR-Cas9筛查的原始数据文件可通过注册号PRJNA471892在NIH序列读取档案(SRA)上获得。PDX RNA测序的原始数据文件可通过注册号PRJNA473981在NIH序列读取档案(SRA)上获得。质谱的原始数据文件可通过ProteomeXchange Consortium在PRIDE合作伙伴存储库中获得,数据集标识符为PXD009830。图形的源数据可以在中找到补充表5。应合理要求,可从相应作者处获得支持本研究结果的所有其他数据。

贡献者

作者贡献

TWC和SPJ构思了该项目,TWC通过执行CRISPR-Cas9筛查启动了该项目。MH进行生物信息分析。LD监督,CL和WW执行类开关重组伊格轨迹不稳定性实验。JJ监督,IK执行重组MAD2L2 co-IP、MAD2L2-IRIF/IB、pS4/8-RPA2 IB和端粒融合实验。MO、PB进行了质谱测定。JC进行并分析了克隆存活实验和随机质粒整合实验;MD和MS生成人类敲除细胞系,MS执行在体外下拉实验;JL开展了寡核苷酸相互作用研究;DP进行末端切除分析;TWC和RB产生的RPE1 p53ko、RPE1 p 53ko/BRCAko和p53ko/BRC A1ko/53BP1ko细胞系;YG和SPJ对上述情况进行了监督。GB在小鼠细胞中进行IR存活,小鼠sgRNA克隆,以及用FY和BF进行FISH。LP进行结构分析和FAM35A建模。QW对细菌表达的重组FAM35A蛋白和EMSA进行纯化。AM、AS、AB和CC对PARPi诱导的耐药性进行了患者源性异种移植实验。VS、MO'C和ZL在本征PARPi电阻实验中建立、执行、分析和表征了PDX。HD协助完成了上述许多实验,并设计和执行了所有其他实验。HD、YG和SPJ根据所有其他人的意见撰写了手稿。LD、JJ、YG和SPJ监督了该项目。

竞争性金融利益

SPJ获得阿斯利康的一些研究资助,是描述PARP抑制剂在癌症治疗中的使用的专利的指定发明人。VS的实验室得到阿斯利康的研究资金支持。MJO和ZL是阿斯利康的员工和股东。其他作者声明没有相互竞争的经济利益。

工具书类

1Goldstein M,Kastan MB。DNA损伤反应:对肿瘤放射和化疗反应的影响。医疗年度收入。2015;66:129–143.[公共医学][谷歌学者]
2Jackson SP,Bartek J.人类生物学和疾病中的DNA损伤反应。自然。2009;461:1071–1078. [PMC免费文章][公共医学][谷歌学者]
三。Schwertman P,Bekker-Jensen S,Mailand N。泛素和泛素样修饰物对DNA双链断裂修复的调节。Nat Rev Mol细胞生物学。2016;17:379–394.[公共医学][谷歌学者]
4Blackford AN、Jackson SP.ATM、ATR和DNA-PK:DNA损伤反应的核心三位一体。分子细胞。2017;66:801–817.[公共医学][谷歌学者]
5Bekker-Jensen S,Mailand N.哺乳动物细胞中DNA双链断裂修复灶的组装和功能。DNA修复(Amst)2010;9:1219–1228.[公共医学][谷歌学者]
6Boersma V等。MAD2L2通过抑制5'末端切除来控制端粒处的DNA修复和DNA断裂。自然。2015;521:537–540. [PMC免费文章][公共医学][谷歌学者]
7Bunting SF等,53BP1通过阻断DNA断裂的切除来抑制Brca1缺陷细胞中的同源重组。单元格。2010;141:243–254. [PMC免费文章][公共医学][谷歌学者]
8Chapman JR等。RIF1对于53BP1依赖的非同源末端连接和抑制DNA双链断裂切除至关重要。分子细胞。2013;49:858–871. [PMC免费文章][公共医学][谷歌学者]
9.Daley JM,Sung P.RIF1在DNA断裂修复途径选择中的应用。分子细胞。2013;49:840–841. [PMC免费文章][公共医学][谷歌学者]
10Di Virgilio M等人。Rif1防止DNA断裂的切除,并促进免疫球蛋白的类别转换。科学。2013;339:711–715. [PMC免费文章][公共医学][谷歌学者]
11Escribano-Diaz C等。由53BP1-RIF1和BRCA1-CtIP组成的细胞周期依赖性调节电路控制DNA修复途径的选择。分子细胞。2013;49:872–883.[公共医学][谷歌学者]
12Xu G等。REV7抵消DNA双链断裂切除并影响PARP抑制。自然。2015;521:541–544. [PMC免费文章][公共医学][谷歌学者]
13Zimmermann M,Lottersberger F,Buonomo SB,Sfeir A,de Lange T.53BP1使用Rif1调节DSB修复,以控制5’端切除。科学。2013;339:700–704. [PMC免费文章][公共医学][谷歌学者]
14Capoluongo E等人,《卵巢癌患者BRCA1/2肿瘤检测指导声明》。塞明·昂科尔。2017;44:187–197.[公共医学][谷歌学者]
15Lord CJ,Ashworth A.PARP抑制剂:临床中的合成致命性。科学。2017;355:1152–1158. [PMC免费文章][公共医学][谷歌学者]
16Bryant HE等。用聚ADP-核糖聚合酶抑制剂特异性杀死BRCA2缺陷肿瘤。自然。2005;434:913–917.[公共医学][谷歌学者]
17Farmer H等。针对BRCA突变细胞中的DNA修复缺陷作为治疗策略。自然。2005;434:917–921.[公共医学][谷歌学者]
18Lord CJ,Ashworth A.针对BRCA突变癌的治疗耐药机制。自然医学。2013;19:1381–1388.[公共医学][谷歌学者]
19Basu B、Sandhu SK、de Bono JS。PARP抑制剂:作用机制及其在癌症预防和治疗中的潜在作用。药物。2012;72:1579–1590.[公共医学][谷歌学者]
20Shalem O等人,《人类细胞中基因组尺度的CRISPR-Cas9基因敲除筛查》。科学。2014;343:84–87. [PMC免费文章][公共医学][谷歌学者]
21Simonetta M等人H4K20me2通过53BP1-RIF1-MAD2L2区分复制前和复制后染色质,以适当地指导DNA修复途径的选择。细胞周期。2018;17(1):124–136. [PMC免费文章][公共医学][谷歌学者]
22Lo KW等。8-kDa动力蛋白轻链与p53结合蛋白1结合,并介导DNA损伤诱导的p53核积累。生物化学杂志。2005;280:8172–8179.[公共医学][谷歌学者]
23Giraud-Panis MJ、Teixeira MT、Geli V、Gilson E.CST与遮蔽物会面以检查端粒。分子细胞。2010;39:665–676.[公共医学][谷歌学者]
24Sun J等。Stn1-Ten1是端粒上的一种Rpa2-Rpa3样复合物。基因发育。2009;23:2900–2914. [PMC免费文章][公共医学][谷歌学者]
25Scully R,Xie A.组蛋白H2AX的双链断裂修复功能。突变研究。2013;750:5–14. [PMC免费文章][公共医学][谷歌学者]
26Kakarougkas A,Jeggo PA。DNA DSB修复路径选择:协调移交机制。英国无线电杂志。2014;8720130685[PMC免费文章][公共医学][谷歌学者]
27Galanty Y等。哺乳动物SUMO E3-ligas PIAS1和PIAS4促进对DNA双链断裂的反应。自然。2009;462:935–939. [PMC免费文章][公共医学][谷歌学者]
28Feldman S等人。53BP1在类开关重组期间对Igh位点染色质拓扑有贡献。免疫学杂志。2017;198:2434–2444. [PMC免费文章][公共医学][谷歌学者]
29Nakamura M等。IgM+B淋巴瘤克隆CH12F3到IgA+细胞的高频类别转换。国际免疫学。1996;8:193–201.[公共医学][谷歌学者]
30Manis JP等。53BP1将DNA损伤反应途径与免疫球蛋白重链类开关重组联系起来。自然免疫学。2004;5:481–487.[公共医学][谷歌学者]
31Ward IM等人53BP1是类开关重组所必需的。细胞生物学杂志。2004;165:459–464. [PMC免费文章][公共医学][谷歌学者]
32Lieber MR人类淋巴染色体易位的机制。Nat Rev癌症。2016;16:387–398. [PMC免费文章][公共医学][谷歌学者]
33Konishi A,de Lange T。TRF2的DNA结合域中的ts突变揭示了端粒保护和NHEJ的细胞周期控制。基因发育。2008;22:1221–1230. [PMC免费文章][公共医学][谷歌学者]
34Dimitrova N,Chen YC,Spector DL,de Lange T.53BP1通过增加染色质流动性促进端粒的非同源末端连接。自然。2008;456:524–528. [PMC免费文章][公共医学][谷歌学者]
35Ochs F等人53BP1培养了同源定向DNA修复的保真度。自然结构分子生物学。2016;23:714–721.[公共医学][谷歌学者]
36Cruz-Garcia A、Lopez-Savedra A、Huertas P.BRCA1加速CtIP介导的DNA末端切除。单元格代表。2014;9:451–459.[公共医学][谷歌学者]
37Gudmundsdottir K,Ashworth A.BRCA1和BRCA2及其相关蛋白在维持基因组稳定性中的作用。致癌物。2006;25:5864–5874.[公共医学][谷歌学者]
38Certo MT等。在单个DNA断点处跟踪基因组工程结果。自然方法。2011;8:671–676. [PMC免费文章][公共医学][谷歌学者]
39Schmidt CK等。系统E2筛查揭示UBE2D-RNF138-CtIP轴促进DNA修复。自然细胞生物学。2015;17:1458–1470. [PMC免费文章][公共医学][谷歌学者]
40Bruna A等人,具有保留的肿瘤内异质性的乳腺癌移植物生物库,用于筛选抗癌化合物。单元格。2016;167:260–274 e222。 [PMC免费文章][公共医学][谷歌学者]
41Cruz C等。RAD51病灶是种系BRCA突变乳腺癌中同源重组修复和PARP抑制剂耐药性的功能性生物标记物。安·昂科尔。2018 [PMC免费文章][公共医学][谷歌学者]
42Bunting SF等。BRCA1在DNA链间交联修复中独立于同源重组发挥作用。分子细胞。2012;46:125–135. [PMC免费文章][公共医学][谷歌学者]
43Gupta R等。DNA修复网络分析揭示屏蔽蛋白是NHEJ和PARP抑制剂敏感性的关键调节器。单元格。2018;173:972–988 e923。 [PMC免费文章][公共医学][谷歌学者]
44Ceccaldi R、Sarangi P、D'Andrea AD。范科尼贫血途径:新参与者和新功能。Nat Rev Mol细胞生物学。2016;17:337–349.[公共医学][谷歌学者]
45Makarova AV,Burgers PM。真核DNA聚合酶zeta。DNA修复(Amst)2015;29:47–55. [PMC免费文章][公共医学][谷歌学者]
46Tomida J等。REV7通过哺乳动物DNA聚合酶zeta中的两个REV3L结合位点对DNA损伤耐受性至关重要。核酸研究。2015;43:1000–1011. [PMC免费文章][公共医学][谷歌学者]
47Bluteau D等。REV7的双等位基因失活与范科尼贫血相关。临床投资杂志。2016;126:3580–3584. [PMC免费文章][公共医学][谷歌学者]
48Chiang TW、le Sage C、Larrieu D、Demir M、Jackson SP.CRISPR-Cas9(D10A)基于尼克斯的基因型和表型筛查,以增强基因组编辑。科学代表。2016;6:24356. [PMC免费文章][公共医学][谷歌学者]
49Peuscher MH,Jacobs JJ。DNA损伤反应和未封盖端粒的修复活性取决于RNF8。自然细胞生物学。2011;13:1139–1145.[公共医学][谷歌学者]
50Galanty Y,Belotserkovskaya R,Coates J,Jackson SP。RNF4,一种SUMO靶向的泛素E3连接酶,促进DNA双链断裂修复。基因发育。2012;26:1179–1195. [PMC免费文章][公共医学][谷歌学者]
51Ong SE等。细胞培养中氨基酸的稳定同位素标记,SILAC,作为一种简单而准确的表达蛋白质组学方法。分子细胞蛋白质组学。2002;1:376–386.[公共医学][谷歌学者]
52Rappsilber J,Mann M,Ishihama Y。使用StageTips对蛋白质组学中的肽进行微纯化、富集、预分离和储存的方案。国家协议。2007;2:1896–1906.[公共医学][谷歌学者]
53Michalski A等。使用高性能台式四极轨道质谱仪Q Exactive的基于质谱的蛋白质组学。分子细胞蛋白质组学。2011;10M111 011015。[PMC免费文章][公共医学][谷歌学者]
54Olsen JV等。肽修饰分析的高能C-trap解离。自然方法。2007;4:709–712.[公共医学][谷歌学者]
55Cox J,Mann M.MaxQuant实现了高肽识别率、个体化p.p.b.范围质量准确度和蛋白质全蛋白定量。国家生物技术。2008;26:1367–1372.[公共医学][谷歌学者]
56.Cox J等人,Andromeda:一个集成到MaxQuant环境中的肽搜索引擎。蛋白质组研究杂志。2011;10:1794–1805.[公共医学][谷歌学者]
57Elias JE,Gygi SP。通过质谱法提高大规模蛋白质鉴定信心的目标经济搜索策略。自然方法。2007;4:207–214.[公共医学][谷歌学者]
58Gustafsson MG等人,通过结构照明在宽场荧光显微镜中实现三维分辨率加倍。生物物理学杂志。2008;94:4957–4970. [PMC免费文章][公共医学][谷歌学者]
59Kapur JN、Sahoo PK、Wong AKC。一种利用直方图熵进行灰度图像阈值化的新方法。计算视觉图。1985;29:273–285. [谷歌学者]
60Agu CA等。从室温下保存长达48小时的血液样本中成功生成人诱导多能干细胞系。干细胞报告。2015;5:660–671. [PMC免费文章][公共医学][谷歌学者]
61Lescale C等。RAG2和XLF/Cernunnos相互作用揭示了RAG复合体在DNA修复中的新作用。国家公社。2016;7:10529. [PMC免费文章][公共医学][谷歌学者]
62Metzakopian E等人。增强基因组编辑工具箱:全基因组CRISPR阵列库。科学代表。2017;7:2244. [PMC免费文章][公共医学][谷歌学者]
63Panchakshari RA等。DNA双链断裂反应因子影响IgH类开关和一般易位连接的末端连接特征。美国国家科学院院刊。2018;115:762–767. [PMC免费文章][公共医学][谷歌学者]
64Yusa K,Zhou L,Li MA,Bradley A,Craig NL.一种用于哺乳动物的超活性piggyBac转座酶。美国国家科学院院刊。2011;108:1531–1536. [PMC免费文章][公共医学][谷歌学者]
65Ahdesmaki MJ,Gray SR,Johnson JH,Lai Z.消除歧义:一个开源应用程序,用于消除嫁接样本下一代测序数据中两个物种的歧义。F1000研究。2016;5:2741. [PMC免费文章][公共医学][谷歌学者]
66Callari M等人。在患者源性肿瘤异种移植物中区分人类和小鼠序列的计算方法。BMC基因组学。2018;19:19. [PMC免费文章][公共医学][谷歌学者]
67Casper J等人。UCSC基因组浏览器数据库:2018年更新。核酸研究。2018;46:D762–D769。 [PMC免费文章][公共医学][谷歌学者]
68.Liao Y,Smyth GK,Shi W.featureCounts:一个有效的通用程序,用于将序列读取分配给基因组特征。生物信息学。2014;30:923–930.[公共医学][谷歌学者]
69McCarthy DJ、Chen Y、Smyth GK。生物变异方面多因子RNA-Seq实验的差异表达分析。核酸研究。2012;40:4288–4297. [PMC免费文章][公共医学][谷歌学者]
70Robinson医学博士、McCarthy DJ、Smyth GK。edgeR:用于数字基因表达数据差异表达分析的Bioconder软件包。生物信息学。2010;26:139–140. [PMC免费文章][公共医学][谷歌学者]