《Molecular Oncology》：Network divergence analysis identifies adaptive gene modules and two orthogonal vulnerability axes in pancreatic cancer
編輯推薦：
本研究運(yùn)用創(chuàng)新的基因-基因共表達(dá)發(fā)散框架，首次在胰腺導(dǎo)管腺癌（PDAC）單細(xì)胞轉(zhuǎn)錄組中系統(tǒng)性地揭示了驅(qū)動(dòng)腫瘤細(xì)胞可塑性的四種動(dòng)態(tài)調(diào)控模塊，并進(jìn)一步整合功能基因組學(xué)數(shù)據(jù)將其歸納為兩個(gè)正交的適應(yīng)性脆弱性軸，為理解PDAC異質(zhì)性、生存機(jī)制及精準(zhǔn)治療提供了全新的網(wǎng)絡(luò)生物學(xué)視角。
在復(fù)雜的腫瘤微環(huán)境中，癌細(xì)胞如何通過重塑其內(nèi)部的基因調(diào)控網(wǎng)絡(luò)，動(dòng)態(tài)適應(yīng)外界壓力并產(chǎn)生治療抵抗，是癌癥生物學(xué)的一大核心問題。尤其是在兇險(xiǎn)的胰腺導(dǎo)管腺癌（PDAC）中，腫瘤細(xì)胞表現(xiàn)出廣泛的異質(zhì)性和顯著的轉(zhuǎn)錄可塑性，這不僅使其極具侵襲性，也導(dǎo)致了極差的治療效果。傳統(tǒng)研究方法，如差異表達(dá)分析和穩(wěn)態(tài)共表達(dá)網(wǎng)絡(luò)，通常只能捕捉穩(wěn)定的轉(zhuǎn)錄狀態(tài)或平均的連接模式，難以揭示那些在不同細(xì)胞狀態(tài)間動(dòng)態(tài)重構(gòu)、驅(qū)動(dòng)適應(yīng)性行為的基因協(xié)調(diào)程序。為了突破這一局限，一項(xiàng)開創(chuàng)性研究采用了新穎的“網(wǎng)絡(luò)發(fā)散分析”框架，旨在解碼PDAC中隱藏的適應(yīng)性邏輯。
這項(xiàng)研究分析了來自42個(gè)PDAC腫瘤的77,155個(gè)惡性細(xì)胞的單細(xì)胞RNA測(cè)序（scRNA-seq）數(shù)據(jù)。通過整合與批次校正后，研究者識(shí)別出10個(gè)主要的惡性細(xì)胞亞群（Leiden聚類）。研究的關(guān)鍵創(chuàng)新在于其“發(fā)散邊”算法：它不是簡(jiǎn)單地比較基因的表達(dá)高低，而是計(jì)算每個(gè)基因?qū)Γ╣ene pair）在這些不同細(xì)胞亞群中的皮爾遜相關(guān)性，并量化這些相關(guān)性在亞群間變化的“發(fā)散”程度。那些相關(guān)性在不同背景下差異最大的基因?qū)Γ础鞍l(fā)散邊”，被認(rèn)為代表了動(dòng)態(tài)重連的、背景依賴性的基因協(xié)調(diào)關(guān)系，是細(xì)胞適應(yīng)性狀態(tài)轉(zhuǎn)換的潛在驅(qū)動(dòng)力。
通過構(gòu)建和可視化這些高發(fā)散邊構(gòu)成的網(wǎng)絡(luò)，并應(yīng)用社區(qū)檢測(cè)，研究驚人地、可重復(fù)地識(shí)別出了四個(gè)清晰且穩(wěn)定的適應(yīng)性轉(zhuǎn)錄模塊。這四個(gè)模塊各自代表了一套緊密協(xié)調(diào)的基因程序，它們突破了經(jīng)典PDAC狀態(tài)（如經(jīng)典型/基底樣）的界限，揭示了更深層次的調(diào)控結(jié)構(gòu)：
- 1.
集成生長(zhǎng)-能量模塊：這個(gè)模塊代表了一個(gè)向增強(qiáng)的生物合成和代謝活性轉(zhuǎn)變的轉(zhuǎn)錄程序。其核心基因富含核糖體蛋白和線粒體氧化磷酸化組件，功能富集于翻譯、細(xì)胞呼吸、內(nèi)質(zhì)網(wǎng)蛋白加工和mRNA剪接等過程。該模塊的啟動(dòng)與MYC和E2F轉(zhuǎn)錄因子網(wǎng)絡(luò)密切相關(guān)，其標(biāo)志性靶點(diǎn)如CDK4和LDHA也顯示出與核心基因的強(qiáng)共表達(dá)。有趣的是，研究者還發(fā)現(xiàn)染色質(zhì)重塑因子CHD9和生長(zhǎng)因子NRG4作為負(fù)相關(guān)的外周調(diào)節(jié)因子，提示該生長(zhǎng)程序受到內(nèi)在反饋回路的約束。
- 2.
應(yīng)激適應(yīng)性轉(zhuǎn)錄模塊：該模塊定義了一個(gè)在經(jīng)受代謝、蛋白質(zhì)毒性或環(huán)境壓力的PDAC細(xì)胞亞群中激活的協(xié)調(diào)基因表達(dá)程序。其核心特征是與整合應(yīng)激反應(yīng)通路緊密相關(guān)，包含ATF4、DDIT3等關(guān)鍵效應(yīng)因子。該模塊對(duì)氨基酸饑餓、自噬、缺氧等特定應(yīng)激類型敏感，但與非葡萄糖饑餓通路關(guān)聯(lián)較弱。模塊中包含大量長(zhǎng)鏈非編碼RNA和一組與DNA損傷反應(yīng)而非DNA復(fù)制更相關(guān)的組蛋白基因，提示了染色質(zhì)水平的應(yīng)激調(diào)控機(jī)制。Notch信號(hào)組件PSENEN作為負(fù)調(diào)節(jié)因子出現(xiàn)，暗示了Notch信號(hào)可能抑制過度的應(yīng)激反應(yīng)激活。
- 3.
IL-2關(guān)聯(lián)的免疫逃逸模塊：此模塊代表了一種獨(dú)特的、與細(xì)胞因子白介素-2相關(guān)的程序，在惡性PDAC細(xì)胞內(nèi)部協(xié)調(diào)表達(dá)多種免疫調(diào)節(jié)基因。深入分析顯示，其富集于IL-2信號(hào)通路，并伴隨IL2受體亞基及下游JAK3/STAT5A/B的共表達(dá)。更重要的是，該模塊核心基因的表達(dá)與多種免疫檢查點(diǎn)分子的表達(dá)呈強(qiáng)正相關(guān)，包括CTLA4、LAG3、TIGIT等，同時(shí)還與免疫抑制細(xì)胞因子TGFB1相關(guān)。這表明該模塊可能通過協(xié)調(diào)表達(dá)多個(gè)免疫抑制分子，幫助腫瘤細(xì)胞逃避免疫系統(tǒng)的攻擊。
- 4.
多通路集體侵襲模塊：這是四個(gè)模塊中最大的一個(gè)，其程序支持一種保留細(xì)胞連接的協(xié)同侵襲模式。功能富集分析指向肌動(dòng)蛋白絲組織、膜運(yùn)輸、細(xì)胞遷移、黏著斑形成等過程。該模塊整合了多種信號(hào)通路，包括Rho GTP酶、受體酪氨酸激酶、PI3K、TGF-β和Notch，其激活劑和抑制劑呈現(xiàn)協(xié)調(diào)的正相關(guān)，表明其維持的是一種精細(xì)的多通路信號(hào)平衡，而非單一通路的過度激活。與單細(xì)胞/間質(zhì)侵襲標(biāo)志物呈負(fù)相關(guān)，而與集體侵襲/上皮標(biāo)志物呈正相關(guān)，進(jìn)一步證實(shí)了其功能特性。自分泌的EGFR-EREG/TGFA配體-受體對(duì)的發(fā)現(xiàn)，為自我強(qiáng)化的EGFR信號(hào)環(huán)路維持定向集體侵襲提供了機(jī)制支持。
僅僅識(shí)別這些模塊并未完整揭示其生物學(xué)意義。研究的另一大突破在于，將這些從單細(xì)胞數(shù)據(jù)中衍生出的模塊活動(dòng)評(píng)分，與功能基因組學(xué)數(shù)據(jù)庫（DepMap）中的CRISPR-Cas9基因依賴性數(shù)據(jù)和藥物敏感性數(shù)據(jù)進(jìn)行整合分析。這一整合從根本上重塑了對(duì)這些模塊的理解，揭示出一個(gè)更高層次的組織架構(gòu)：這四個(gè)適應(yīng)性模塊并非四個(gè)獨(dú)立的平行狀態(tài)，而是坍縮為兩個(gè)正交的脆弱性軸。
- •
生物合成-代謝軸：主要由IGE模塊構(gòu)成。高IGE活性的細(xì)胞系強(qiáng)烈依賴于線粒體、核糖體和RNA加工相關(guān)基因，并顯示出對(duì)PARP抑制劑（如他拉唑帕尼）的選擇性敏感，這與MYC/E2F驅(qū)動(dòng)的、依賴持續(xù)翻譯和氧化磷酸化的生物合成狀態(tài)所伴隨的復(fù)制壓力和DNA修復(fù)需求相吻合。
- •
應(yīng)激-免疫-侵襲軸：由SAT、IL2和MPC三個(gè)模塊匯聚而成。盡管它們?cè)趩渭?xì)胞層面代表不同的生物學(xué)過程，但在功能依賴性上卻共享特征。SAT模塊依賴內(nèi)質(zhì)網(wǎng)蛋白折疊和蛋白質(zhì)穩(wěn)態(tài)調(diào)節(jié)因子；IL2模塊關(guān)聯(lián)自噬和細(xì)胞因子信號(hào)相關(guān)因子；MPC模塊則依賴細(xì)胞極性、細(xì)胞骨架重塑和NF-κB信號(hào)相關(guān)因子。在藥物反應(yīng)上，這三個(gè)模塊均表現(xiàn)出對(duì)加劇蛋白質(zhì)毒性負(fù)荷或干擾細(xì)胞骨架/信號(hào)傳導(dǎo)的化合物的相對(duì)敏感性，而對(duì)多種DNA損傷劑顯示出相對(duì)抵抗。
這種“兩軸”架構(gòu)僅在將發(fā)散邊模塊映射到功能基因組學(xué)空間時(shí)才顯現(xiàn)出來，凸顯了整合多維度數(shù)據(jù)在理解腫瘤生物學(xué)中的強(qiáng)大力量。
這些適應(yīng)性程序同樣具有臨床相關(guān)性。在TCGA-PAAD隊(duì)列的生存分析中，研究者發(fā)現(xiàn)SAT模塊的高活性與較差的總體生存期相關(guān)，MPC模塊也顯示出相似的不良趨勢(shì)。當(dāng)將SAT和MPC模塊的活動(dòng)合并考慮時(shí)，定義了一個(gè)明確的“應(yīng)激-侵襲”不良預(yù)后軸。相比之下，IGE模塊活動(dòng)與總體風(fēng)險(xiǎn)無顯著關(guān)聯(lián)，但一個(gè)由最優(yōu)切點(diǎn)定義的IGE高活性亞組卻顯示出生存期適度改善的趨勢(shì)。IL2模塊在整體隊(duì)列中未顯示出獨(dú)立的預(yù)后分層能力，這可能與其僅局限于少數(shù)惡性細(xì)胞亞群有關(guān)。在調(diào)整了病理分期等臨床病理學(xué)協(xié)變量的多變量模型中，模塊評(píng)分未能保持獨(dú)立的預(yù)后意義，表明這些適應(yīng)性程序可能是影響疾病進(jìn)展階段的上游調(diào)控狀態(tài)，而非獨(dú)立的生物標(biāo)志物。
這項(xiàng)研究通過聚焦于基因-基因協(xié)調(diào)關(guān)系的動(dòng)態(tài)重連，而非靜態(tài)的表達(dá)水平，成功地繪制了PDAC細(xì)胞適應(yīng)性的核心調(diào)控圖譜。所識(shí)別的四個(gè)適應(yīng)性模塊及其衍生的兩個(gè)脆弱性軸，不僅加深了我們對(duì)PDAC異質(zhì)性和可塑性的理解，更重要的是，它們將轉(zhuǎn)錄程序與實(shí)驗(yàn)可驗(yàn)證的基因依賴性和藥物敏感性直接聯(lián)系起來，為基于腫瘤特定適應(yīng)性狀態(tài)的精準(zhǔn)治療策略提供了全新的理論框架和潛在的靶點(diǎn)線索。未來的研究需要進(jìn)一步在空間維度、實(shí)驗(yàn)擾動(dòng)和縱向治療背景下驗(yàn)證這些模塊的因果角色和動(dòng)態(tài)演變，以最終實(shí)現(xiàn)將這些網(wǎng)絡(luò)生物學(xué)見解轉(zhuǎn)化為改善患者預(yù)后的有效療法。
