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