膽管癌栓(BDTT)是一種可能發(fā)生于晚期肝細(xì)胞癌(HCC)患者的嚴(yán)重并發(fā)癥,可導(dǎo)致生活質(zhì)量迅速下降,并與較差的臨床結(jié)局相關(guān)。然而,由于其細(xì)胞起源和發(fā)病機(jī)制不明,合并BDTT的HCC患者常被誤診誤治。因此,開發(fā)一種有效的肝癌合并膽管癌栓的診斷和治療方法是臨床上的迫切要求。既往研究主要集中在探索BDTT的HCC患者的臨床病理或治療措施,而對(duì)BDTT的發(fā)展過(guò)程知之甚少。根據(jù)初步的臨床研究,有研究推測(cè)高度侵襲性的腫瘤細(xì)胞可能通過(guò)鄰近膽管壁的上皮下侵襲,導(dǎo)致癌栓的形成。此外,他們注意到BMI1在有BDTT12的HCC患者中表達(dá)上調(diào)。BMI1在一系列細(xì)胞過(guò)程中發(fā)揮關(guān)鍵作用,包括干細(xì)胞自我更新、增殖和保存基因組穩(wěn)定性。此外,它涉及多種癌癥,包括血液系統(tǒng)惡性腫瘤和實(shí)體瘤,其過(guò)表達(dá)常與侵襲性腫瘤行為、腫瘤起始細(xì)胞(TICs)的存在和不良預(yù)后相關(guān)。miRNA是已知的小非編碼分子,參與調(diào)節(jié)基因表達(dá),其異常表達(dá)與TICs特征相關(guān)。然而,由于缺乏專門設(shè)計(jì)的BDTT動(dòng)物模型,BMI1及其下游效應(yīng)miRNA在BDTT中的潛在作用仍未被探索。溶酶體組織蛋白酶B(CTSB)是溶酶體半胱氨酸蛋白酶家族的一種成員,積極調(diào)節(jié)細(xì)胞外基質(zhì)成分的蛋白水解、破壞細(xì)胞內(nèi)通訊、降解基底膜和改變細(xì)胞-細(xì)胞之間的相互作用,因此在癌癥進(jìn)展中具有多重潛在作用。然而,BMI1是否通過(guò)表觀遺傳調(diào)控CTSB在TICs中的分泌從而通過(guò)miRNA調(diào)控BDTT的發(fā)育仍不清楚。該研究發(fā)表在《Nature Communications》,IF: 16.6。
技術(shù)路線
主要研究結(jié)果
1. 在HCC患者中,BMI1表達(dá)升高與TIC相關(guān)蛋白/基因標(biāo)簽表達(dá)升高、BDTT發(fā)生率升高和不良預(yù)后相關(guān)
由于BMI1表達(dá)與HCC患者預(yù)后之間的關(guān)系尚未得到充分探索,因此研究者對(duì)來(lái)自癌癥基因組圖譜(TCGA)數(shù)據(jù)庫(kù)的臨床和RNA測(cè)序數(shù)據(jù)進(jìn)行了生物信息學(xué)分析。此外,對(duì)收集的HCC患者隊(duì)列進(jìn)行了IHC檢查,發(fā)現(xiàn)在HCC患者中,BMI1的高表達(dá)水平與較差的總生存期、增加的腫瘤大小和升高的癌癥進(jìn)展相關(guān)(圖1a-c)。與之前的研究相比,研究者還收集了一個(gè)相對(duì)較大的有或無(wú)BDTT的HCC患者的樣本隊(duì)列,并使用他們對(duì)BMI1表達(dá)進(jìn)行了IHC分析,結(jié)果表明有BDTT的患者與較差的總生存率相關(guān)(圖1d),而BDTT和BMI1高表達(dá)的患者與較低的總生存率相關(guān)(圖1e)。此外,研究者發(fā)現(xiàn)BMI1表達(dá)的升高區(qū)分了那些BDTT發(fā)生率增加的HCC患者(圖1f)。
鑒于BMI1表達(dá)在控制TIC自我更新中的作用,進(jìn)一步研究了HCC患者中TIC相關(guān)標(biāo)簽與BMI1蛋白和mRNA表達(dá)之間的關(guān)系。在腫瘤區(qū)域中,TIC相關(guān)蛋白標(biāo)簽和BMI1的表達(dá)呈正相關(guān)(圖1g,h)。進(jìn)一步的臨床研究表明,在HCC患者中,TIC相關(guān)蛋白/基因標(biāo)簽的高表達(dá)與BDTT發(fā)病率的增加、較差的總生存率相關(guān)(圖1i-k)。最后,研究者發(fā)現(xiàn)與低表達(dá)TIC相關(guān)蛋白/基因標(biāo)簽和BMI1的HCC患者相比,高表達(dá)TIC相關(guān)蛋白/基因標(biāo)簽和BMI1的HCC患者的總生存率較差(圖1l,m)。綜上所述,該研究揭示了BMI1與TIC相關(guān)標(biāo)志物的表達(dá)之間的潛在關(guān)聯(lián),提示BMI1在BDTT的發(fā)生和HCC的癌癥進(jìn)展中發(fā)揮作用。這一發(fā)現(xiàn)促使研究者深入研究其在調(diào)節(jié)癌細(xì)胞特性中的功能。
圖1 在HCC中,BMI1高表達(dá)與TIC相關(guān)基因標(biāo)簽表達(dá)增強(qiáng)、BDTT發(fā)生率增加和總體生存率降低相關(guān)
2. BMI1high腫瘤起始細(xì)胞來(lái)源于惡性轉(zhuǎn)化的肝祖細(xì)胞
為了檢測(cè)BMI1在HCC發(fā)生中的細(xì)胞功能,研究者在大鼠LPC細(xì)胞系WBF344和人HCC細(xì)胞系中對(duì)BMI1的表達(dá)進(jìn)行了遺傳修飾(圖2a,b)。研究者的結(jié)果表明,與對(duì)照細(xì)胞相比,BMI1過(guò)表達(dá)增加了WBBMI1和PLCBMI1細(xì)胞的細(xì)胞生長(zhǎng)(圖2c,d)、集落形成(圖2e,f)、細(xì)胞遷移(圖2g,h)和細(xì)胞侵襲(圖2i,j)以及腫瘤球形成(圖2k,l)。
為了確定穩(wěn)定表達(dá)BMI1的細(xì)胞是否可以通過(guò)旁分泌效應(yīng)調(diào)節(jié)其侵襲能力,研究者在transwell matrigel侵襲系統(tǒng)中,將WBCtrl或PLCCtrl細(xì)胞暴露于來(lái)自WBBMI1/WBCtrl細(xì)胞或PLCBMI1/PLCCtrl細(xì)胞的條件培養(yǎng)基(CM),表明與分別暴露于來(lái)自WBCtrl和PLCCtrl細(xì)胞的CM的細(xì)胞相比,暴露于來(lái)自WBBMI1/PLCBMI1細(xì)胞的CM增加了WBCtrl和PLCCtrl細(xì)胞的侵襲(圖2m-o)。接下來(lái),研究者通過(guò)向裸鼠皮下注射少量WBBMI1或WBCtrl細(xì)胞來(lái)檢測(cè)它們?cè)隗w內(nèi)的成瘤性,表明注射5000個(gè)WBBMI1細(xì)胞足以在裸鼠體內(nèi)形成皮下腫瘤,并且100%有效。相反,WBCtrl細(xì)胞未能產(chǎn)生任何可觸及的腫瘤(圖2p,q)??傮w而言,這些發(fā)現(xiàn)表明,BMI1的過(guò)表達(dá)可誘導(dǎo)LPC惡性轉(zhuǎn)化為BMI1high TICs,稱為WBBMI1。
圖2 過(guò)表達(dá)BMI1可自發(fā)地將肝祖細(xì)胞轉(zhuǎn)化為腫瘤起始細(xì)胞
3. 小鼠肝/脾原位注射WBBMI1細(xì)胞可產(chǎn)生原發(fā)腫瘤/轉(zhuǎn)移灶,并發(fā)生自發(fā)性膽管瘤栓
為了檢查BMI1high TICs是否可以啟動(dòng)原位腫瘤形成和BDTT,研究者在裸鼠中將WBBMI1細(xì)胞原位植入肝臟的Glisson包膜(圖3a),表明40%植入WBBMI1細(xì)胞的小鼠在注射后4周出現(xiàn)黃疸(圖3b)。黃疸的出現(xiàn)是合并BDTT2的HCC患者的臨床表現(xiàn)之一。磁共振成像顯示,將WBBMI1細(xì)胞原位植入裸鼠肝臟Glisson包膜后4周內(nèi)產(chǎn)生了肝臟腫瘤。而WBCtrl移植細(xì)胞未能形成任何腫瘤(圖3c)。研究者還表明,WBBMI1植入細(xì)胞以100%的發(fā)生率產(chǎn)生了原發(fā)性肝臟腫瘤,而40%的WBBMI1荷瘤小鼠出現(xiàn)了膽囊和膽管增大(圖3d)。相比之下,WBCtrl植入小鼠均未顯示任何這些表型(圖3c,d)。對(duì)腫瘤連續(xù)切片的進(jìn)一步組織病理學(xué)研究表明,在黃疸的WBBMI1原位腫瘤荷瘤小鼠中存在膽管癌栓,但在植入WBCtrl細(xì)胞的小鼠中未見膽管癌栓(圖3e)。重要的是,研究者證明了瘤栓是由原發(fā)腫瘤部位的WBBMI1細(xì)胞侵襲膽管腔引起的(圖3e)。由于細(xì)胞角蛋白7(CK7)和細(xì)胞角蛋白19(CK19)一般表達(dá)于膽管上皮,因此研究者在WBCtrl或WBBMI1細(xì)胞移植的小鼠的肝組織/腫瘤切片中進(jìn)行了CK7和CK9的IHC染色,表明在WBBMI1移植的黃疸小鼠的癌栓表面存在CK7和CK19陽(yáng)性的膽管上皮碎片(圖3e)。
為了進(jìn)一步證實(shí)研究者的觀察結(jié)果,研究者在裸鼠中進(jìn)行了脾內(nèi)注射WBCtrl/WBBMI1細(xì)胞(圖3f),表明注射WBBMI1細(xì)胞的小鼠中約有17%在注射后4周出現(xiàn)黃疸,并且血清膽紅素水平升高(圖3g)。研究者的數(shù)據(jù)顯示,WBBMI1細(xì)胞脾內(nèi)種植形成原發(fā)性脾腫瘤的發(fā)生率為100%(圖3h),而67%的小鼠形成肝轉(zhuǎn)移,約17%的WBBMI1種植小鼠出現(xiàn)膽囊和膽管增大(圖3h)。再次,腫瘤連續(xù)切片的組織病理學(xué)研究顯示,在僅黃疸的WBBMI1移植小鼠中存在BDTT,同時(shí)還表明,癌栓源于肝轉(zhuǎn)移灶WBBMI1細(xì)胞侵襲膽管腔(圖3i)。隨后對(duì)這些組織切片中CK7和CK19表達(dá)的IHC分析支持了這一觀察結(jié)果(圖3i)。相比之下,WBCtrl移植小鼠沒有發(fā)生原發(fā)性脾腫瘤和肝轉(zhuǎn)移以及BDTT(圖3h,i)。總體而言,研究者的數(shù)據(jù)表明,原位肝/脾移植BMI1high TICs可以產(chǎn)生原發(fā)性腫瘤/肝轉(zhuǎn)移,并誘導(dǎo)自發(fā)性BDTT。更重要的是,本研究建立了肝細(xì)胞癌膽管癌栓的臨床前動(dòng)物模型,為探索膽管癌栓的分子機(jī)制提供了機(jī)會(huì)。
圖3 WBBMI1細(xì)胞原位移植入小鼠肝臟/脾臟可形成原發(fā)性肝臟腫瘤/轉(zhuǎn)移灶,并誘導(dǎo)自發(fā)性膽管瘤栓形成。
4. 比較蛋白質(zhì)組學(xué)分析顯示,CTSB介導(dǎo)的細(xì)胞成分和通路在BMI1high TICs中富集
接下來(lái),研究者試圖研究BMI1high TIC介導(dǎo)BDTT的分子機(jī)制。通過(guò)與WBCtrl和WBBMI1細(xì)胞進(jìn)行比較蛋白質(zhì)組學(xué)分析,與WBCtrl細(xì)胞相比,分泌因子CTSB在WBBMI1細(xì)胞中表達(dá)上調(diào),而進(jìn)一步的GO富集分析表明,在WBBMI1細(xì)胞中,CTSB相關(guān)的細(xì)胞成分和通路富集(圖4a,b)。能夠降解細(xì)胞外基質(zhì),并破壞細(xì)胞黏附和連接,從而刺激癌細(xì)胞侵襲。蛋白質(zhì)印跡分析表明,與對(duì)照細(xì)胞相比,WBBMI1/PLCBMI1細(xì)胞中CTSB的表達(dá)在翻譯水平上調(diào),但在mRNA水平?jīng)]有上調(diào)(圖4c,d)。此外,ELISA和CTSB活性測(cè)定實(shí)驗(yàn)均表明,與來(lái)自WBCtrl或PLCCtrl細(xì)胞的CM相比,來(lái)自WBBMI1或PLCBMI1細(xì)胞的CM中CTSB的表達(dá)水平和酶活性顯著增強(qiáng)(圖4e-h)。在臨床上,研究者的數(shù)據(jù)表明,在HCC患者中,CTSB表達(dá)升高與較差的總生存期和癌癥進(jìn)展增加相關(guān)(圖4i,k),但其表達(dá)與腫瘤大小無(wú)關(guān)(圖4l)。有趣的是,研究者還發(fā)現(xiàn)在HCC中,CTSB的高表達(dá)水平與BDTT的發(fā)生率增加相關(guān)(圖4m,n)。重要的是,患BDTT的HCC患者血清CTSB水平上調(diào)(圖4o)。此外,研究者發(fā)現(xiàn)在HCC患者中,BMI1和CTSB的高表達(dá)與較差的總生存期強(qiáng)相關(guān)(圖4p),而它們的表達(dá)也呈正相關(guān)(圖4q)。最后,研究者的數(shù)據(jù)顯示,在HCC患者中,TIC相關(guān)蛋白標(biāo)簽和CTSB的高表達(dá)與較差的總生存期強(qiáng)相關(guān)(圖4r),同時(shí)它們的表達(dá)呈正相關(guān)(圖4s)。綜上所述,該研究表明,CTSB可能作為BMI1介導(dǎo)的HCC進(jìn)展的下游效應(yīng)分子,為HCC合并BDTT患者提供了一個(gè)有價(jià)值的血清標(biāo)志物的前景。
圖4 比較蛋白質(zhì)組學(xué)分析揭示了BMI1high TICs中CTSB依賴的細(xì)胞成分和通路的富集
5. CTSB抑制劑治療禁止BDTT,并在不影響原發(fā)腫瘤生長(zhǎng)的情況下延長(zhǎng)小鼠生存期
為了確定CTSB在BMI1介導(dǎo)的癌細(xì)胞特性中的作用,研究者對(duì)WBBMI1/PLCBMI和BMI1敲除的MHCC97H細(xì)胞中的CTSB表達(dá)進(jìn)行了遺傳修飾。與非沉默對(duì)照siRNA(siNSC)轉(zhuǎn)染/安慰劑處理的細(xì)胞相比,通過(guò)siRNA或CTSB抑制劑CA-074敲除CTSB降低了WBBMI1和PLCBMI1細(xì)胞的侵襲能力(圖5a-d)。此外,與使用siNSC轉(zhuǎn)染的WBBMI1/PLCBMI1細(xì)胞的CM處理的細(xì)胞相比,使用來(lái)自CTSB敲除的WBBMI1/PLCBMI1細(xì)胞的CM處理的細(xì)胞相比,使用來(lái)自CTSB敲除的WBBMI1/PLCBMI1細(xì)胞的CM處理的細(xì)胞降低了其侵襲性(圖5e,f)。接下來(lái),研究者創(chuàng)建了體外跨肝內(nèi)膽管上皮遷移模型,將HCC細(xì)胞置于單層的人肝內(nèi)膽管上皮細(xì)胞上,在transwell小室的上室中,transwell小室涂有Matrigel基質(zhì),下室裝載DMEM + 20%胎牛血清(圖5g)。值得注意的是,研究者的數(shù)據(jù)表明,BMI1過(guò)表達(dá)增強(qiáng)了人PLC細(xì)胞的跨肝內(nèi)膽管上皮遷移,而通過(guò)siRNA或抑制劑去除CTSB降低了PLCBMI1細(xì)胞中觀察到的增強(qiáng)的跨肝內(nèi)膽管上皮遷移能力(圖5h,i)。此外,與scramble轉(zhuǎn)染細(xì)胞相比,在MHCC97H細(xì)胞中去除BMI1降低了其跨肝內(nèi)膽管上皮遷移能力(圖5j)。
為了進(jìn)一步驗(yàn)證體內(nèi)觀察結(jié)果,研究者將WBBMI1細(xì)胞注射到裸鼠皮下,然后用安慰劑或CTSB抑制劑CA-074處理裸鼠。正如預(yù)期,CA-074治療對(duì)WBBMI1皮下腫瘤的生長(zhǎng)沒有影響(圖5k,1)。相反,與安慰劑治療的小鼠相比,在WBBMI1原位荷瘤小鼠中,CA-074降低了黃疸、膽紅素水平升高、膽囊和膽管增大或BDTT的發(fā)生率(圖5m-o),而在這些小鼠中,它不影響腫瘤的形成和原發(fā)腫瘤的生長(zhǎng)(圖5n-p)。最后,與安慰劑治療組相比,CTSB抑制劑治療延長(zhǎng)了WBBMI1原位荷瘤小鼠的生存期(圖5q)。綜上所述,這些發(fā)現(xiàn)表明CTSB在BMI1依賴的膽管癌栓中發(fā)揮重要作用,并確定了肝癌合并膽管癌栓患者的潛在治療靶點(diǎn)。
圖5 CTSB抑制劑治療可抑制WBBMI1荷瘤小鼠的膽管腫瘤血栓形成,與原發(fā)腫瘤生長(zhǎng)無(wú)關(guān)
6. BMI1high TICs通過(guò)抑制miR-218-1-3p的表達(dá)促進(jìn)CTSB驅(qū)動(dòng)的膽管癌血栓形成
目前的研究數(shù)據(jù)表明BMI1在蛋白水平上調(diào)節(jié)CTSB的表達(dá)。研究者使用三種不同的在線搜索工具,使用CTSB mRNA序列預(yù)測(cè)其調(diào)控miRNA,而維恩圖分析表明,這些預(yù)測(cè)的miRNA中有41個(gè)被所有這些工具共同識(shí)別(圖6a)。接下來(lái),研究者在BMI1基因修飾的肝癌細(xì)胞中對(duì)這些miRNA進(jìn)行了RT-PCR定量。特別是,與對(duì)照細(xì)胞相比,miR-218-1-3p在BMI1穩(wěn)定表達(dá)的PLC細(xì)胞中表達(dá)下調(diào),而在BMI1敲除的MHCC97H細(xì)胞中表達(dá)顯著增加(圖6b)。維恩圖分析表明,miR-218-1-3p是only與BMI1表達(dá)相關(guān)的miRNA(圖6c)。接下來(lái),研究者使用之前發(fā)表的PRE序列分析了miR-218-1-3p的啟動(dòng)子,表明其啟動(dòng)子包含2個(gè)假定的PRE位點(diǎn)(圖6d),而染色質(zhì)免疫沉淀分析表明,與PLCCtrl細(xì)胞相比,PLCBMI1細(xì)胞中BMI1和三甲基化H3K27在這兩個(gè)PRE位點(diǎn)上的結(jié)合增加(圖6e)。BMI1可能與組蛋白甲基轉(zhuǎn)移酶EZH2協(xié)同作用,通過(guò)重塑染色質(zhì)結(jié)構(gòu)抑制miR-213-1-3p的表達(dá)。為了研究miR-218-1-3p對(duì)CTSB表達(dá)的調(diào)控作用,分別在WBBMI1, MHCC97H和PLCBMI1細(xì)胞中過(guò)表達(dá)miR-218-1-3p,結(jié)果顯示與陰性對(duì)照模擬物轉(zhuǎn)染細(xì)胞相比,在表達(dá)miR-218-1-3p的WBBMI1/MHCC97H細(xì)胞/PLCBMI1細(xì)胞中CTSB的表達(dá)降低(圖6f,g)。這表明溶酶體CTSB是潛在的靶點(diǎn)之一,因?yàn)槠銫TSB的3'UTR區(qū)域攜帶保守的miR-218-1-3p種子匹配序列(圖6h)。為了確定miR-218-1-3p是否與CTSB的3'UTR結(jié)合位點(diǎn)結(jié)合,將含有野生型或突變型CTSB 3'UTR結(jié)合位點(diǎn)的熒光素酶報(bào)告載體與miR-218-1-3p模擬物或NC模擬物共轉(zhuǎn)染MHCC97H和PLCBMI1細(xì)胞。值得注意的是,在轉(zhuǎn)染了攜帶WT-CTSB 3'UTR結(jié)合位點(diǎn)序列的報(bào)告載體的MHCC97H和PLCBMI1細(xì)胞中,miR-218-1-3p過(guò)表達(dá)抑制了熒光素酶活性,而在轉(zhuǎn)染了含有突變的CTSB 3'UTR結(jié)合位點(diǎn)序列的報(bào)告載體構(gòu)建的細(xì)胞中,miR-218-1-3p不影響熒光素酶活性(圖6i)。這些結(jié)果表明,BMI1和EZH2通過(guò)miR-218-1-3p協(xié)同調(diào)控CTSB的表達(dá)。
通過(guò)對(duì)HCC腫瘤切片進(jìn)行miR-218-1-3p的熒光原位雜交實(shí)驗(yàn),研究者發(fā)現(xiàn)在HCC腫瘤中,miR-218-1-3p的表達(dá)水平分別與BMI1/CTSB的表達(dá)呈負(fù)相關(guān)(圖6j)。為了檢測(cè)體內(nèi)miR-218-1-3p表達(dá)在BDTT中的重要性,將miR-218-1-3p或NC模擬物穩(wěn)定轉(zhuǎn)染W(wǎng)BBMI1細(xì)胞的裸鼠肝臟原位注射到裸鼠肝臟中,表明與NC模擬物轉(zhuǎn)染組相比,miR-219-1-3p過(guò)表達(dá)降低了WBBMI1種植小鼠的黃疸、膽紅素升高、膽囊和膽管增大或BDTT的發(fā)生率,而不影響原發(fā)性肝臟腫瘤的生長(zhǎng)(圖6k-n)。
總之,該研究揭示了BMI1過(guò)表達(dá)將LPCs惡性轉(zhuǎn)化為BMI1high TICs,而在小鼠原位肝移植中,這些細(xì)胞可以產(chǎn)生原發(fā)腫瘤和自發(fā)性BDTT。在機(jī)制上,BMI1通過(guò)PRC依賴的機(jī)制沉默miR-218-1-3p的表達(dá),然后表觀遺傳上調(diào)BMI1high TICs中的CTSB分泌,從而促進(jìn)其侵襲膽管形成BDTT。
圖6 BMI1通過(guò)抑制miR-218-1-3p的表達(dá)表觀遺傳上調(diào)CTSB的表達(dá)
結(jié)論
該研究結(jié)果解決了肝癌生物學(xué)領(lǐng)域中關(guān)于BDTT細(xì)胞起源和發(fā)病機(jī)制的謎題,表明BMI1high TICs誘導(dǎo)肝癌發(fā)生并經(jīng)歷CTSB驅(qū)動(dòng)的膽管侵犯形成癌栓。更重要的是,該研究為HCC合并BDTT患者提供了潛在的診斷標(biāo)志物和治療途徑,為未來(lái)的臨床應(yīng)用打開了大門。
機(jī)制圖
實(shí)驗(yàn)方法
細(xì)胞培養(yǎng)和轉(zhuǎn)染,體外腫瘤發(fā)生試驗(yàn),蛋白質(zhì)印跡,RT-PCR,免疫組化,免疫染色實(shí)驗(yàn),熒光素酶報(bào)告基因檢測(cè),體內(nèi)致瘤性測(cè)定,外科原位肝/脾癌模型,蛋白質(zhì)組學(xué)分析,經(jīng)肝內(nèi)膽道上皮遷移試驗(yàn),ELISA檢測(cè),CTSB活性測(cè)定法,染色質(zhì)免疫沉淀實(shí)驗(yàn),熒光原位雜交(FISH)實(shí)驗(yàn)
參考文獻(xiàn)
Xu LB, Qin YF, Su L, Huang C, Xu Q, Zhang R, et al. Cathepsin-facilitated invasion of BMI1-high hepatocellular carcinoma cells drives bile duct tumor thrombi formation. Nat Commun. 2023 Nov 3;14(1):7033. doi: 10.1038/s41467-023-42930-y.