10X單核轉錄組 && 空間轉錄組揭示癌癥相關成纖維細胞是宮頸鱗狀細胞癌預后和治療的潛在靶點

hello氛雪,大家好济丘,這一篇我們要看一下單核轉錄組和空間轉錄組的聯(lián)合分析驮肉,參考文章在Transcriptomic Profiling Reveals Cancer-Associated Fibroblasts as Potential Targets for the Prognosis and Treatment of Cervical Squamous Cell Carcinoma戚长,其實最主要的目的就是在于了解單核轉錄組與空間轉錄組的聯(lián)合分析及運用田藐。

圖片.png

ABSTRACT

為了了解宮頸鱗狀細胞癌 (CSCC) 的病因、結構和免疫學特征兑徘,對來自 20 個人的宮頸樣本進行了單核 RNA 測序 (snRNA-seq) 和空間轉錄組學 (ST) 實驗。在探索影響 CSCC 個體內(nèi)免疫異質(zhì)性的可能因素時羡洛,確定了一組在某些腫瘤周圍富集的癌癥相關成纖維細胞 (CAF)挂脑,它們高度表達 ACTA2、POSTN欲侮、ITGB4 和 FAP崭闲。結果表明,CAFs可能通過抑制淋巴細胞浸潤和重塑腫瘤細胞外基質(zhì)來支持腫瘤的生長和轉移威蕉。此外刁俭,高 CAF 信號預測 CSCC 患者的臨床結果較差。數(shù)據(jù)還揭示了 HPV 在腫瘤中的感染情況韧涨、參與宮頸癌病變進展的關鍵因素以及腫瘤代謝與免疫反應強度之間的關聯(lián)牍戚。總的來說虑粥,研究結果可能會改善 CSCC 的預后和治療方法如孝。

INTRODUCTION

Cervical cancer is the fourth most common cancer affecting women’s health globally, especially in low- and middle-income regions. Since 2018, the World Health Organization (WHO) has called for the global elimination of cervical cancer, quantifying actions in vaccination, screening, and disease treatment/management,which require joint efforts from different parties for decades(這一段大家了解即可).

Currently, over 12 types of human papillomaviruses (HPVs) are carcinogenic 娩贷。 其中第晰,HPV16占宮頸癌病例的60%-70%,尤其是宮頸鱗狀細胞癌(CSCC)彬祖。 HPV 編碼三種癌蛋白茁瘦,包括 E5、E6 和 E7储笑,它們可以破壞細胞周期并導致不受控制的細胞分裂甜熔。 CSCC的發(fā)展一般要經(jīng)歷低級別鱗狀上皮內(nèi)病變(LSIL)、高級別鱗狀上皮內(nèi)病變(HSIL)和浸潤性宮頸癌南蓬,這可能需要幾十年的時間纺非。 CSCC 的腫瘤微環(huán)境 (TME) 是由病毒哑了、癌細胞、基質(zhì)細胞和免疫細胞之間復雜的相互作用形成的烧颖,但由于技術限制和缺乏研究弱左,尚未完全表征。 破譯推動疾病進展的因素可能會產(chǎn)生新的預后和干預策略炕淮。

早期宮頸癌行根治性子宮切除術雖然可以獲得良好的預后拆火,但晚期宮頸癌的5年總生存率或無病生存率均不盡如人意。目前涂圆,化療(如紫杉醇们镜、順鉑、貝伐珠單抗等)和放療仍是轉移或復發(fā)患者的主要姑息治療策略润歉,反應率低(48%)模狭,生存期短(17個月).免疫療法通過逆轉疲憊或受抑制的免疫活動,為治療無法治愈的宮頸癌帶來新希望踩衩。針對程序性細胞死亡 1 (PD1)嚼鹉、程序性細胞死亡配體 1 (PD-L1) 和細胞毒性 T 淋巴細胞抗原 4 (CTLA4) 的免疫檢查點阻斷 (ICB) 藥物目前正在試驗中用于治療復發(fā)性/轉移性宮頸癌。不幸的是驱富,ICB 治療的總體反應率很低锚赤,從 4% 到 26% 不等。闡明 CSCC 的免疫狀況褐鸥,尤其是 TME 中的免疫抑制狀態(tài)线脚,可能有助于我們更好地解決這一現(xiàn)象并調(diào)整我們對宮頸癌的治療策略

單細胞測序和空間轉錄組學 (ST) 是揭示腫瘤細胞異質(zhì)性和微環(huán)境的最先進工具叫榕,但此類技術在 CSCC 研究中的應用仍然很少浑侥。 在這項研究中,收集了 20 個人的宮頸樣本翠霍,并結合單核 RNA 測序 (snRNA-seq) 和 ST 來研究病因锭吨、結構 和 CSCC 的免疫學特征。 破譯 HPV 誘導的 CSCC 的轉錄組變化將為 CSCC 的診斷寒匙、預后和治療提供新的見解零如。

RESULTS

snRNA-seq data revealed the cellular composition of CSCC

為了充分表征宮頸組織的細胞組成,收集了 5 名患者的 CSCC 樣本進行 snRNA-seq锄弱。共有67,003個細胞和30,996個基因通過質(zhì)控考蕾,根據(jù)經(jīng)典細胞標志物從中鑒定出14種細胞類型,包括癌細胞(6,960)会宪、柱狀上皮細胞(CECs肖卧,22,396)、內(nèi)皮細胞(6,340)掸鹅、平滑肌細胞(4,502)塞帐、成纖維細胞 (9,836)拦赠、B 細胞 (689)、單核細胞 (5,281)葵姥、T 細胞 (4,930)荷鼠、調(diào)節(jié)性 T (Treg) 細胞 (1,081)、漿細胞 (3,236)榔幸、髓樣樹突細胞 (DC) (955 )允乐、漿細胞樣 DC (272)、肥大細胞 (384) 和自然殺傷 (NK) 細胞 (141)削咆。子宮頸包含兩種類型的細胞在其表面排列牍疏,子宮頸外層為復層鱗狀上皮細胞,子宮頸內(nèi)膜和隱窩上有簡單的柱狀上皮細胞拨齐。鱗狀上皮細胞的發(fā)育不良導致 CSCC鳞陨。因此,癌細胞主要表達已知的CSCC相關基因SERPINB3(Serpin Family B Member 3)瞻惋,鱗狀細胞的腫瘤基因TP63炊邦、CDKN2A和角蛋白基因KRT15。由于組織主要來自晚期癌癥患者(FIGO IB2-IIIC1期熟史;在化療前宮頸活檢期間收集沒有分期信息的組織),因此很少分離出正常的上皮鱗狀細胞窄俏。 snRNA-seq counts與高風險 HPV 參考基因組的映射揭示了癌細胞中病毒基因的存在蹂匹。還鑒定了一大群柱狀上皮細胞,它們高度表達 MUC5B 和 WFDC2凹蜈。這種細胞類型主要位于子宮頸內(nèi)皮上皮限寞,但也可以出現(xiàn)在成人子宮頸和一些腺體的鱗柱交界處。平滑肌細胞仰坦、成纖維細胞和內(nèi)皮細胞是構成宮頸間質(zhì)的主要細胞類型履植。平滑肌細胞高表達MYH11、MYLK悄晃、ACTG2玫霎、COL3A1和COL1A1,成纖維細胞除COL3A1和COL1A1外還高表達LAMA2妈橄,而內(nèi)皮細胞可通過EMCN庶近、FLT1和EGFLT的高表達來區(qū)分。除了子宮頸的結構細胞外眷蚓,還鑒定了多種免疫細胞類型鼻种,其中單核細胞(ITGAX、MX4A7)沙热、T細胞(CD3E叉钥、CD247)和漿細胞(MZB1罢缸、IGKC)最為豐富。簡而言之投队,snRNA-seq 數(shù)據(jù)提供了一個全面的圖譜枫疆,顯示了 CSCC 組織的結構和免疫細胞組成(細胞marker注意搜集一下)。

圖片.png

圖片.png

Spatial transcriptomic characterization of CSCC

空間信息對于理解組織中的細胞間相互作用至關重要蛾洛,不幸的是养铸,這在 snRNA-seq 數(shù)據(jù)中被忽略了。因此轧膘,利用 SpaTial Enhanced REsolution Omics-sequencing (Stereo-seq) 技術來獲取原位基因表達譜钞螟。獲得了來自 2 名非癌癥患者和 14 名 CSCC 患者的宮頸樣本,并將其嵌入 OCT谎碍。從每個 OCT 塊中解剖出 10 μm 厚的連續(xù)冷凍切片鳞滨,用于立體測序、蘇木精和伊紅 (HE) 染色以及免疫組織化學 (IHC) 染色蟆淀。最后拯啦,共成功獲得 18 張 ST 載玻片(非癌,n=3熔任;CSCC褒链,n=15)。 Stereo-seq 芯片中的捕獲點直徑為 220 nm疑苔,兩個相鄰點之間的中心距為 500 nm甫匹。捕獲spot被分組到箱中以包含足夠的基因以進行準確的聚類。初步分析顯示惦费,腫瘤區(qū)域的 RNA 豐度比基質(zhì)區(qū)域高得多兵迅。為了平衡腫瘤和基質(zhì)之間的表達差異,在 bin(100 x 100 點)處對 CSCC ST 載玻片進行了注釋薪贫,以充分展示組織成分恍箭,其覆蓋面積約為 49.72 x 49.72 μm。 CSCC ST 載玻片的每個 bin 的平均基因數(shù)介于 1,767 到 4,152 之間瞧省。由于來自非癌癥患者的三張 ST 載玻片的基因表達強度低于 CSCC 載玻片扯夭,因此它們在 bin200 (99.72 x 99.72 μm) 處進行了注釋。統(tǒng)一流形近似和投影 (UMAP) 分析表明,CSCC 和非癌癥的 bin cluster傾向于彼此分離,而 CSCC 的 bin cluster則顯示出一些收斂吨述。

考慮到癌組織的復雜結構,根據(jù)專業(yè)病理評估(基于HE和IHC染色結果)和標記基因表達模式手動進行了ST載玻片的初始注釋藕筋。在分析的 CSCC 樣本中,通常鑒定出六種類型的組織cluster,包括腫瘤隐圾、間質(zhì)(無明顯炎癥)伍掀、炎癥(間質(zhì)具有彌漫性炎癥或局灶性炎癥)、腺體暇藏、血管和壞死蜜笤。腫瘤、基質(zhì)和炎癥cluster廣泛分布在 CSCC 樣本的 ST 載玻片中盐碱,某些樣本包含壞死把兔、腺體和血管。根據(jù)基因表達譜瓮顽,組織cluster可以進一步分為帶有數(shù)字后綴的subcluster县好。應該記住,ST cluster實際上是一個定義區(qū)域內(nèi)的小區(qū)的混合物暖混,但由其主要特征指定缕贡。例如,腫瘤cluster不完全由癌細胞組成拣播,也可能包含其他細胞類型晾咪,但數(shù)量很少,如免疫細胞贮配、成纖維細胞等谍倦。組織特異性基因也顯示出空間模式。 KRT5泪勒、CDKN2A剂跟、SERPINB3等癌性鱗狀細胞相關基因主要富集在ST腫瘤區(qū),IGKC和IGLC2主要富集在炎癥區(qū)酣藻,VIM在間質(zhì)區(qū),ADRA2A在血管中鳍置,MUC5B在腺體中富集辽剧。一般而言,在 ST 腫瘤區(qū)域中觀察到的轉錄和翻譯活性税产、細胞增殖怕轿、氧化磷酸化和免疫反應高于其他區(qū)域。 ST 載玻片的初始手動注釋用于輔助 TME 的下游分析辟拷。

圖片.png

圖片.png

圖片.png

圖片.png

圖片.png

圖片.png

關于HPV的病因學作用撞羽,病毒基因在宮頸組織中的表達一直是人們感興趣的話題。 將 ST 測序讀數(shù)與 HPV 參考基因組作圖衫冻,發(fā)現(xiàn)捐獻 15 個樣本的所有 14 名 CSCC 患者均為 HPV 陽性诀紊,其中 85.7%(12/14)被 HPV16 感染,7.1%(1/14)被 HPV33 感染隅俘,并且 HPV58 為 7.1% (1/14)邻奠。 HPV 讀數(shù)覆蓋了病毒基因組的 8% - 100%(約 600 bp - 7905 bp)笤喳,主要在 ST 腫瘤區(qū)域被識別。 空間可視化展示了同一樣本不同腫瘤部位病毒基因的不同捕獲信號碌宴,包括 E5杀狡、E6、E7贰镣、和 L1 經(jīng)常觀察到呜象。 In contrast, only marginal HPV reads were identified in the non-CSCC samples 。

圖片.png

圖片.png

ST data resolved the evolutionary trajectory of CSCC(空間進化軌跡)

由于 ST 數(shù)據(jù)能夠揭示同一組織切片內(nèi)不同腫瘤區(qū)域的空間分布碑隆,這為研究 CSCC 的進化軌跡提供了機會恭陡。高級別上皮內(nèi)瘤變、原位癌和浸潤癌代表宮頸的不同病理狀態(tài)干跛,它們經(jīng)常同時發(fā)生在同一患者身上子姜。然而,由于技術限制楼入,同一個人內(nèi)這些病理過程的動態(tài)變化仍然模糊不清哥捕。在這里,首先通過計算其 CytoTRACE 評分(關于CytoTRACE嘉熊,大家可以參考我的文章10X單細胞軌跡分析(擬時分析)之cytotrace遥赚,非常棒的方法)來評估 ST 載玻片中每個腫瘤亞群的分化狀態(tài)。同時使用 HE 染色來識別宮頸上皮及其相鄰腫瘤區(qū)域阐肤。癌性上皮區(qū)被認為是浸潤前癌病灶凫佛,而鄰近腫瘤區(qū)域的CytoTRACE評分相對較高,表明分化程度相對較低孕惜,被認為是浸潤性癌病灶愧薛。由于遠離上皮的腫瘤的轉移路徑難以確定,因此該分析不包括此類腫瘤區(qū)域衫画。最后毫炉,確定了來自獨立患者的兩個樣本(TJH37 和 TJH90),其中包含侵襲前和侵襲性病變削罩。cluster 14 和 16 被識別為侵襲前病變瞄勾,而cluster 12 和 13 分別被確定為 TJH37 和 TJH90 的侵襲性病變。

基于侵襲前病變與侵襲性病變(Log 2FC > 0.5 或 Log 2FC < -0.5弥激,P < 0.001)的差異表達基因(DEGs)的功能分析表明进陡,上調(diào)基因主要參與皮膚發(fā)育、角質(zhì)形成細胞分化微服、 對傷口和細菌的反應趾疚、維生素 D 受體 (VDR) 通路、血管發(fā)育和中性粒細胞脫顆粒,主要與分化和對有害刺激(如 HPV 感染盗蟆、上皮損傷)的反應增強有關戈二。 相比之下,下調(diào)的基因與 RNA 代謝喳资、氧化磷酸化觉吭、線粒體組織、細胞呼吸和病毒 mRNA 翻譯顯著相關仆邓,表明侵襲性癌癥病變中的代謝增加鲜滩。 VDR 通路在侵襲前病變中的積極參與值得注意,這表明其在宮頸粘膜免疫中的重要作用节值。

然后徙硅,試圖查明參與癌性病變進展的關鍵基因。兩個樣品共有 16 個 DEG(Log2FC > 1.0 或 Log2FC < -1.0搞疗,P < 0.001)嗓蘑。偽時間分析表明,12個基因(ANXA1匿乃、AL136982.6桩皿、CRCT1、S100A7幢炸、S100A9泄隔、SPRR1A、SPRR1B宛徊、SPRR2A佛嬉、SPRR2D、SPRR2E闸天、KRTDAP和KRT17)被其他浸潤前癌病變基因過表達(而CD74暖呕、MT2A、KRT5 和 KRT15) 往往被侵襲性病變過度表達苞氮。由于早期預后對于宮頸上皮內(nèi)瘤變 (CIN) 患者的分類至關重要湾揽,我們很好奇這些 DEG 是否具有任何臨床價值。我們從基因表達綜合(GEO)數(shù)據(jù)庫中下載了 GSE63514 數(shù)據(jù)集葱淳,其中包含 66 個宮頸浸潤前(CIN)和 38 個浸潤性病變(宮頸癌)樣本,并進行 PCA 分析以驗證 DEG 的區(qū)分能力抛姑。結果表明赞厕,11個DEGs的組合(其余5個基因在該數(shù)據(jù)集中未鑒定)可以很好地區(qū)分浸潤前(CIN)和浸潤性病變(宮頸癌)。隨機森林 (RF) 測試顯示這 11 個基因(ANXA1定硝、S100A9皿桑、S100A7、KRT15、SPRR1A诲侮、SPRR1B镀虐、CD74、KRTDAP沟绪、MT2A刮便、CRCT1 和 KRT5)具有相似的重要性。 RF 分類模型也表現(xiàn)出良好的性能绽慈,在 10 次測試中達到了 0.848 的平均 AUC(ROC 曲線下面積)恨旱。總之坝疼,這些 DEG 可能為 CIN 和癌癥之間的鑒別診斷提供新的見解搜贤。

圖片.png

圖片.png

Variable immune inhibition in CSCC(免疫抑制)

關于 ICB 治療宮頸癌的低反應率,決定仔細檢查 CSCC 的免疫狀況以尋找線索钝凶。在 snRNA-seq 和 ST 數(shù)據(jù)中評估了具有不同免疫功能的三個基因組的表達譜仪芒,即共刺激、細胞毒性/效應子和共抑制/耗竭耕陷。在單細胞水平上掂名,共刺激基因在先天免疫和適應性免疫細胞中均有表達,尤其是在 Treg啃炸、T 和 NK 細胞中铆隘。發(fā)現(xiàn) Treg 細胞高度表達 CD27、CD28南用、CD40LG膀钠、ICOS、TNFRSF18裹虫、TNFRSF4 和 TNFRSF9肿嘲。雖然這些基因是 Treg 細胞成熟和正常抑制功能所必需的,但 Treg 細胞中 CD27 的過度表達可能會抑制抗腫瘤免疫反應筑公。在空間上雳窟,共刺激基因傾向于在腫瘤和炎癥區(qū)域富集。特別是匣屡,TNFRSF18(也稱為糖皮質(zhì)激素誘導的TNF受體封救,GITR)在炎癥和腫瘤區(qū)域中普遍表達。然而捣作,其在非癌癥樣本的上皮細胞中的表達也上調(diào)誉结。在snRNA-seq 數(shù)據(jù)中,該基因主要在 Treg券躁、NK惩坑、T 細胞和肥大細胞中檢測到掉盅。雖然TNFRSF18與腫瘤中Treg細胞的免疫抑制有關,但其高空間表達水平可能是由多種類型的免疫細胞共同貢獻的以舒。免疫細胞毒/效應基因主要由 T 和 NK 細胞表達趾痘,其中一些在 ST 載玻片的炎癥和腫瘤區(qū)域通常上調(diào),包括 GNLY蔓钟、GZMA永票、GZMB 和 NKG7。這些基因主要由 NK 細胞表達奋刽,表明 NK 細胞在針對 CSCC 的細胞毒性反應中起著不可或缺的作用瓦侮。對于共抑制/耗竭基因,我們未能在我們的 ST 數(shù)據(jù)中檢測到 CTLA4 和 PD-1 的任何顯著表達佣谐,盡管 Treg 細胞可以高度表達 CTLA4肚吏,Treg、T 和 NK 細胞可以高度表達 PD-1狭魂。 PD-L1 主要由漿細胞樣 DCs 表達罚攀,僅在一小部分 ST 樣本的腫瘤或炎癥區(qū)域過度表達。雖然 CD276雌澄、ENTPD1斋泄、IDO1、LGALS9 和 VSIR 通常在 ST 樣本中檢測到镐牺,但與非癌癥樣本相比炫掐,只有 IDO1 和 LGALS9 在 CSCC 樣本中似乎具有更高和更廣泛的表達。這兩個基因都由 DC 表達睬涧,并且可以下調(diào)細胞毒性 T 細胞的活性募胃。 IDO1 和 LGALS9 是否可以成為比 CTLA4 和 PD-L1 更好的 ICB 治療 CSCC 的靶點仍有待探索。此外畦浓,當我們放大檢查同一ST載玻片中的免疫基因時痹束,它們的表達在不同腫瘤區(qū)域之間可能會有很大差異。在樣本 TJH08 中讶请,腫瘤通常表達高 IDO1祷嘶、低 PD-L1 和非常低的 CTLA4。相比之下夺溢,樣本 TJH37 中只有一個腫瘤區(qū)域表達這些基因论巍。總的來說风响,盡管我們的 snRNA-seq 和 ST 數(shù)據(jù)都顯示了 CSCC 患者免疫衰竭的證據(jù)嘉汰,但患者之間和患者內(nèi)部的免疫微環(huán)境差異很大

圖片.png

Hypermetabolic tumors were associated with low immune response

代謝可以調(diào)節(jié)腫瘤的免疫微環(huán)境钞诡,這可能是癌癥治療的推定干預目標郑现。分別對六個途徑進行了基因集變異分析 (GSVA),包括缺氧荧降、乳酸接箫、糖酵解、脂質(zhì)代謝朵诫、磷酸戊糖和氧化磷酸化途徑辛友。然后將上述六種途徑的 GSVA 評分的平均值計算為每個腫瘤區(qū)域的代謝評分。根據(jù) GSVA 代謝評分剪返,排名前 20 位的 ST 腫瘤cluster被歸類為高代謝腫瘤废累,而那些排名最后 20 位的為低代謝腫瘤。一般來說脱盲,高代謝腫瘤在氧化磷酸化邑滨、糖酵解和乳酸途徑中表現(xiàn)出更高的活性,表明增殖的癌細胞中存在活躍的有氧糖酵解钱反,即 Warburg 效應掖看。此外,高代謝腫瘤還伴隨著嚴重的缺氧和活躍的脂質(zhì)代謝面哥,表明快速生長的腫瘤存在強烈的氧化和營養(yǎng)應激哎壳。進一步觀察到代謝和免疫反應之間的負相關。發(fā)現(xiàn)高代謝腫瘤經(jīng)歷了較低水平的適應性和先天免疫干擾尚卫、淋巴細胞浸潤和淋巴管生成归榕。根據(jù)每個 bin(100 x 100 個點)中的基因模塊表達評分進一步探討了腫瘤分化、缺氧和免疫之間的關系吱涉。較低的腫瘤分化水平(由較高的 CytoTRACE 評分表示)與較嚴重的缺氧呈正相關(r = 0.28刹泄,P < 0.001)。此外邑飒,缺氧與適應性免疫呈負相關(r = -0.19循签,P < 0.001),尤其是淋巴細胞浸潤(r = -0.22疙咸,P < 0.001)县匠。值得注意的是,Th1細胞(r = -0.17撒轮,P < 0.001)乞旦、B細胞(r = 0.30,P < 0.001)题山、未成熟B細胞(r = -0.18兰粉,P < 0.001)和肥大細胞(r = -0.22, P < 0.001) 在缺氧腫瘤中顯著降低。雖然缺氧對先天免疫評分沒有顯著影響顶瞳,但在缺氧的腫瘤中發(fā)現(xiàn) CD56-NK 細胞增加(r = 0.27玖姑,P < 0.001)和未成熟的樹突狀細胞(r = 0.23愕秫,P < 0.001),這可能無法有效對抗癌細胞焰络。進一步嘗試分別使用 TOP2A 和 PTPRC 作為增殖細胞和免疫細胞的標志物在空間上驗證代謝和免疫之間的相關性戴甩。在樣本 TJH34 和 TJH35 中,能夠在同一 ST 載玻片中識別高代謝和低代謝腫瘤區(qū)域闪彼。在低代謝腫瘤區(qū)域內(nèi)外檢測到更高的 PTPRC 表達甜孤,與 HE 圖像中的淋巴細胞分布模式一致。這些觀察結果表明畏腕,由活躍增殖的腫瘤引起的氧氣和營養(yǎng)缺乏環(huán)境可能會阻止 B 細胞缴川、肥大細胞以及功能性 DC、NK 細胞和輔助 T 細胞的積累描馅“芽洌總的來說,分化差的腫瘤往往代謝更活躍铭污,這可能導致缺氧和缺乏免疫干預扎即。

圖片.png

圖片.png

圖片.png

Identification of a cluster of cancer-associated fibroblasts (CAFs) in CSCC

在探索樣本 TJH34 中高代謝和低代謝腫瘤區(qū)域之間的免疫差異時,注意到高代謝腫瘤區(qū)域外有一個獨特的空間cluster况凉。這個cluster不同于大多數(shù)間質(zhì)cluster谚鄙,看起來像一條包裹著腫瘤的絲帶。分析認為這個特定的cluster可能對塑造腫瘤區(qū)域之間的免疫異質(zhì)性至關重要刁绒。由于該cluster是基質(zhì)的一部分闷营,仔細檢查了 snRNA-seq 數(shù)據(jù)中的成纖維細胞。幸運的是知市,鑒定了來自所有五個樣本的一小組成纖維細胞傻盟,它們高度表達了 CAF 的報告標記基因,包括 ACTA2嫂丙、POSTN娘赴、ITGB4 和 FAP。 CAFs 的干細胞評分低于癌細胞和成纖維細胞跟啤,并且處于不同的細胞周期階段诽表。基于標志基因集(MSigDB v7.4隅肥,https://www.gsea-msigdb.org/gsea/msigdb/)的功能富集表明CAF與成纖維細胞和癌細胞具有共同的活動竿奏。 CAF 參與成纖維細胞的類似途徑,包括紫外線反應減弱腥放、血管生成泛啸、肌生成和上皮-間質(zhì)轉化。對于包括 p53 通路秃症、KRAS 信號轉導候址、雌激素反應吕粹、有絲分裂紡錘體、G2/M 檢查點和 E2F 目標在內(nèi)的通路岗仑,CAF 顯示出與癌細胞相似的活性昂芜。在基因水平,CAFs不僅高表達成纖維細胞的標志基因赔蒲,如膠原蛋白家族(COL1A1、COL3A1良漱、COL4A1舞虱、COL5A2愈魏、COL6A3等)巨坊,還高表達惡性鱗狀細胞的標志基因,如KRT4和KRT13角蛋白家族的巩检。盡管癌癥中的 CAFs 可能有不同的起源患久,但上皮特征表明 CSCC 樣本中 CAFs 的起源可能與 EMT 相關椅寺。

為了定位 CSCC 組織中 CAF 的空間分布,采用了 Moncada 等人開發(fā)的多模態(tài)交叉分析 (MIA) 方法蒋失,整合snRNA-seq和ST數(shù)據(jù)(關于MIA的方法返帕,大家可以參考我的文章MIA用于單細胞和空間的聯(lián)合分析。簡而言之篙挽,該方法計算了由snRNA-seq數(shù)據(jù)識別的細胞類型特異性基因與由ST數(shù)據(jù)表征的區(qū)域特異性基因的表達水平的重疊程度荆萤。由此產(chǎn)生的 p 值越小,在后面的描述中被稱為 MIA 分數(shù)铣卡,定義的細胞類型和 ST 區(qū)域之間的相關性就越強链韭。初始 MIA 結果表明,ST 聚類結果符合相應區(qū)域的預期細胞組成煮落。不幸的是敞峭,單獨的 MIA 分數(shù)不能完全反映細胞的空間特異性,尤其是在 RNA 豐度低的區(qū)域蝉仇。因此旋讹,同時利用經(jīng)實驗證實與 CAF 相關的 POSTN 的高表達水平和 CAF 的高 MIA 分數(shù)來定義 CAF 的 ST cluster。結果表明轿衔,CAFs 在 ST 載玻片中的一些腫瘤區(qū)域周圍富集骗村,包括樣本 TJH34 的高代謝腫瘤區(qū)域。通過使用相同樣本的連續(xù)組織切片對 POSTN 進行 IHC 染色呀枢,進一步證實了 CSCC 中 CAF 的存在胚股。值得注意的是,并非所有腫瘤區(qū)域都被 CAF 包圍裙秋,這讓我們對與這些細胞存在相關的生物學差異感到好奇琅拌。

圖片.png

CAFs might facilitate the growth and metastasis of CSCC from diverse aspects

為了全面揭示 CAFs 在 CSCC 中的生物學功能缨伊,將 ST 載玻片中的腫瘤區(qū)域分為兩種類型:CAFs 包圍的腫瘤區(qū)域(CAFs+腫瘤)和未被 CAFs 包圍的腫瘤區(qū)域(CAFs-腫瘤)。 發(fā)現(xiàn)三張 ST 載玻片包含 CAFs+ 和 CAFs- 腫瘤區(qū)域进宝,并用于下游分析刻坊。 然后使用這3個樣本的CAFs+和CAFs-腫瘤區(qū)域之間上調(diào)和下調(diào)的DEG進行GO富集分析。 結果表明党晋,CAFs+腫瘤在能量使用谭胚、代謝、有絲分裂和細胞生長方面比CAFs-腫瘤更活躍未玻。 同時灾而,細胞粘附、凋亡和免疫反應在 CAFs+ 腫瘤中下調(diào)扳剿。 上述觀察結果與 CSCC 的免疫和代謝異質(zhì)性非常吻合旁趟,尤其是在樣品 TJH34 中。 這些表明CAFs的存在可以從不同方面支持腫瘤進展庇绽。

圖片.png

圖片.png

接下來锡搜,計算了 993 個單個細胞關于免疫基因集的基因模塊表達分數(shù),以評估免疫細胞豐度瞧掺。 結果顯示 CAFs+ 腫瘤中 B 細胞耕餐、CD4 T 細胞、CD8 T 細胞辟狈、中性粒細胞蛾方、DC、NK 細胞和 Th1 細胞的數(shù)量顯著減少上陕。 有趣的是桩砰,在 CAFs+ 腫瘤中發(fā)現(xiàn)了更多的巨噬細胞,可能是腫瘤相關巨噬細胞 (TAM)释簿,這可能促進癌細胞的增殖和遷移亚隅。 因此,CAFs 不僅可以作為物理屏障阻止促免疫細胞浸潤到腫瘤區(qū)域庶溶,還可以招募抗免疫巨噬細胞以促進腫瘤的生長煮纵。

作為基質(zhì)的一部分,CAF 必須與癌細胞偏螺、基質(zhì)細胞和免疫細胞密切相互作用行疏。事實上,對 snRNA-seq 數(shù)據(jù)的分析表明套像,CAF 與其他細胞之間在細胞外基質(zhì) (ECM) 形成和細胞-細胞接觸方面存在復雜的相互作用酿联。 CAFs 高表達膠原家族基因,特別是 COL1A1、COL1A2贞让、COL4A1周崭、COL4A2、COL4A5喳张、COL6A1续镇、COL6A2 和 COL6A3)與免疫細胞和平滑肌細胞表達的 CD44 相互作用,可能參與細胞粘附和遷移销部。膠原蛋白還與由癌細胞摸航、免疫細胞和基質(zhì)細胞表達的整合素家族的不同成員相互作用。同樣舅桩,由 CAF 表達的 FN1(纖連蛋白 1酱虎,一種可溶性糖蛋白)和層粘連蛋白(LAMA2、LAMA3江咳、LAMA4、LAMA5哥放、LAMB1歼指、LAMB2、LAMB3甥雕、LAMC1踩身、LAMC3)也通過整合素與其他細胞類型相互作用。整聯(lián)蛋白是由 α 和 β 亞基組成的膜受體蛋白社露,參與細胞粘附和識別挟阻。有趣的是,CAF 使用不同的異二聚體形式的整聯(lián)蛋白來接觸其他細胞類型峭弟。它們可能通過由α2β1附鸽、α3β1和αvβ8亞基組成的整合素與癌細胞相互作用,同時通過由α9β1瞒瘸、α6β1和α1β1亞基組成的整合素與內(nèi)皮細胞坷备、CEC、平滑肌細胞相互作用情臭,并通過與T細胞和NK細胞相互作用省撑。由亞基 α1β1 組成的整合素。重要的是俯在,CAF 可以利用 F11R(也稱為 JAM1竟秫,連接粘附分子 1)通過 F11R 和 JAM3 與癌細胞和基質(zhì)細胞形成緊密連接,這可能會阻止免疫細胞的浸潤跷乐。 CAF 還可能表達其他基質(zhì)蛋白肥败,包括 THBS1(血小板反應蛋白 1)、THBS2 和 TNC(生腱蛋白 C),以通過 CD44拙吉、整合素(α3β1)和 SDC1 與免疫細胞潮孽、平滑肌細胞、癌細胞和漿細胞進行交流筷黔。此外往史,CAF 過度表達了幾種組織重塑因子,包括 POSTN(骨膜蛋白佛舱,一種分泌的細胞外基質(zhì)蛋白)椎例、FAP(成纖維細胞活化蛋白,一種絲氨酸蛋白酶)请祖、MMP1(基質(zhì)金屬蛋白酶 1)订歪、TNC(肌腱蛋白-C,一種基質(zhì)蛋白)肆捕、和 LOXL1(賴氨酰氧化酶 1刷晋,催化膠原蛋白和彈性蛋白的交聯(lián))。這一證據(jù)表明 CAF 在塑造腫瘤細胞外環(huán)境中的關鍵作用慎陵。

除了在 ECM 構建中的作用外眼虱,CAF 還可以通過過度表達分泌因子(包括 SEMA3C、POSTN席纽、CXCL6)來增強癌細胞的干性和增殖捏悬。據(jù)報道,SEMA3C 可促進癌癥干細胞維持润梯、血管生成和侵襲 36,37过牙。 POSTN 可以通過整合素 αvβ3 激活 Akt/PKB 通路來增加癌細胞的存活率。它還可能通過 PTK7-Wnt/β-Catenin 信號通路促進癌癥生長纺铭。 CXCL6(C-X-C 基序趨化因子配體 6)雖然主要與免疫反應有關寇钉,但據(jù)報道可促進食管鱗狀細胞癌的生長和轉移。 CAF 中另一個高度表達的基因 SNAI2 (Slug) 是一種與蝸牛相關的鋅指轉錄因子舶赔,可抑制細胞凋亡并促進癌癥進展摧莽。 CAF 還表達了幾種常見的生長因子,例如 TGFB1(轉化生長因子β1)顿痪、EGF(表皮生長因子)和 VEGFA(血管內(nèi)皮生長因子 A)镊辕。此外,CAFs 中 LSD1(組蛋白賴氨酸去甲基化酶 1)的上調(diào)可能會抑制 IFN 激活以逃避免疫攻擊 42蚁袭。CAFs 產(chǎn)生的 Wnt5a 信號蛋白也可能抑制免疫反應以促進腫瘤轉移征懈。總之揩悄,CAFs 可能能夠增強 TME卖哎,促進腫瘤的進展。

The presence of CAFs was associated with poorer outcomes of CSCC patients

為了驗證 CAF 的促腫瘤作用,首先使用來自癌癥基因組圖譜 (TCGA)45 的數(shù)據(jù)集進行生存分析亏娜,其中包含 252 名 CSCC 患者焕窝。 使用標記基因組(ACTA2、POSTN维贺、ITGB4 和 FAP)計算每位 CSCC 患者的 CAF GSVA 評分它掂。 毫不奇怪,較高的 CAF 信號預測了不利的無進展生存期(HR = 1.66溯泣,95%CI = 1.03-2.67虐秋,P = 0.038)和總生存期(HR = 1.69,95%CI = 1.00-2.84垃沦,P = 0.05) 用于 CSCC 患者客给。 接下來,測量了 POSTN(一種 CAF 的生物標志物)在 61 個存檔的福爾馬林固定石蠟包埋 (FFPE) CSCC 樣本的癌癥病變附近的基質(zhì)中的表達水平肢簿。 相關性分析表明靶剑,POSTN 表達水平越高,分化越差池充、腫瘤分期越長桩引、淋巴結轉移越頻繁、外周血鱗狀細胞抗原 (SCC) 濃度越高纵菌,這進一步證實了 CAF 的促腫瘤發(fā)生能力阐污。

總的來說休涤,分析結果表明 CAFs 是 CSCC TME 的重要組成部分咱圆,形成屏障以保護癌細胞免受免疫監(jiān)視和清除,分泌細胞因子以刺激細胞增殖和血管生成功氨,抑制細胞凋亡序苏,并重塑 ECM 以促進腫瘤轉移 . 考慮到 CAFs 的惡化作用,它們可能成為 CSCC 預后和治療的潛在靶點捷凄。

圖片.png

圖片.png

圖片.png

DISCUSSION

盡管疫苗和根治性子宮切除術是預防和治療宮頸癌的有效措施忱详,但復發(fā)/轉移宮頸癌的治療仍然是實現(xiàn)消除宮頸癌目標的一大障礙。 在此跺涤,我們結合 snRNA-seq 和 ST 技術描繪了 CSCC 的高分辨率轉錄組圖譜匈睁,詳細描述了 CSCC 的病因、結構和免疫學特征桶错,這可能有助于 HPV 誘導的宮頸癌的管理和治療航唆。

數(shù)據(jù)顯示 HPV 在癌性宮頸鱗狀細胞中的轉錄或翻譯活躍,致癌基因(E5院刁、E6 和 E7)高度表達糯钙。由于 HPV 具有 DNA 基因組,從轉錄組數(shù)據(jù)中恢復幾乎完整的病毒基因組是表明晚期腫瘤中病毒活性活躍的生動證據(jù)。 HPV 基因組包含兩個多聚腺苷酸化位點任岸。早期多聚腺苷酸化 (AE) 位點位于 E5 和 L2 之間再榄,由 E6、E7享潜、E1困鸥、E2、E4 和 E5 的轉錄本共享米碰。晚期多聚腺苷酸化 (AL) 位點位于 L1 開放閱讀框之后窝革,L1、L2 甚至 E7 和 E4 的轉錄本都使用它吕座。由于 Stereo-seq 使用 polyT 捕獲 mRNA虐译,因此與聚腺苷酸化位點相鄰的基因可能有更高的測序機會,這反映在 E5 和 L1 的普遍高于其他基因的表達上吴趴。盡管如此漆诽,早期基因,尤其是 E5锣枝、E6 和 E7厢拭,往往比包括 L1 和 L2 在內(nèi)的晚期基因具有更高的表達。因為 HPV 在分化的鱗狀細胞中組裝和釋放子代病毒粒子撇叁,所以在增殖的癌細胞中早期基因的表達高于晚期基因是合理的供鸠。此外,雖然與 A E 位點不相鄰陨闹,但癌基因 E6 和 E7 在腫瘤中顯示出與 E5 相似甚至更高的表達水平楞捂,表明它們在宮頸鱗狀細胞的腫瘤發(fā)生中起關鍵作用。事實上趋厉,針對病毒 E6 和 E7 基因的治療性疫苗在治療宮頸癌前病變和癌癥方面已顯示出一些有希望的結果寨闹。

傳統(tǒng)的bulk RNA 測序和單細胞 RNA 測序通常不足以獲得連續(xù)病理狀態(tài)的基因表達譜。在此君账,我們使用包含侵襲前和侵襲性癌癥病變的組織樣本繁堡,通過尖端 ST 方法揭示了與發(fā)育不良進展相關的關鍵基因和途徑。發(fā)現(xiàn)在侵襲前狀態(tài)上調(diào)的 VDR 途徑對于維持粘膜屏障穩(wěn)態(tài)和預防微生物感染至關重要乡数。 VDR 通路可以被 Toll 樣受體 (TLR) 上調(diào)椭蹄,以誘導抗菌因子,如人體中的導管素净赴。此外绳矩,據(jù)報道,服用維生素 D 可以改善由結核病感染引起的肺部病變劫侧。因此埋酬,TLR-VDR 軸可能是對抗 HPV 感染的關鍵哨啃。 TLR 家族的激動劑,如咪唑喹啉和瑞喹莫特写妥,以及維生素 D拳球,可能是治療癌前病變和癌病變的可能療法。此外珍特,在我們的測試數(shù)據(jù)中祝峻,侵襲前和侵襲性病變之間的DEG 在區(qū)分 CIN 與宮頸癌方面表現(xiàn)出良好的性能。四種上調(diào)基因扎筒,包括 CD74莱找、MT2A、KRT5 和 KRT15嗜桌,與幾種類型的鱗狀細胞癌的腫瘤進展相關奥溺,它們在 CSCC 中的預后價值可能值得進一步研究。

如今骨宠,ICB 療法浮定,尤其是那些使用 PD-L1/PD-1 和 CTLA4 抑制劑的療法,是治療轉移性宮頸癌的新方法之一层亿。幾項研究報告了 PD-L1 在宮頸癌中的廣泛表達桦卒,陽性率從 34% 到 96% 不等。然而匿又,單獨的 PD-L1 表達與宮頸癌患者的疾病結果無關 61方灾。事實上,不同試驗對 PD-1/PD-L1 和 CTLA-4 抑制劑的反應率波動很大碌更,而且療效似乎獨立于相關檢查點基因的表達狀態(tài)裕偿。在我們的研究中,除 LGALS9 和 IDO1 外针贬,CSCC 腫瘤和炎癥區(qū)域中大多數(shù)免疫檢查點基因的表達水平并不顯著高于非癌樣本中的表達水平击费。 LGALS9(即半乳糖凝集素 9)通過與 T 細胞表面的 Tim-3 結合或抑制 DC 的抗原呈遞能力來下調(diào)效應 T 細胞免疫拢蛋。雖然已證明半乳糖凝集素 9 信號通路的破壞可誘導患有胰腺導管腺癌的小鼠的腫瘤消退桦他,但在肺轉移小鼠模型中報告了相反的效果。 IDO1 主要在 DC 中表達谆棱,有助于將色氨酸降解為犬尿氨酸快压,從而抑制 T 細胞功能。發(fā)現(xiàn)抑制IDO1可以增強HeLa和SiHa腫瘤球細胞的放射敏感性垃瞧,表明IDO1抑制劑作為放射增敏劑的潛在應用蔫劣。靶向 LGALS9 和 IDO1 是否可以改善 CSCC 的治療需要進一步探索。我們的研究還揭示了高氧和高營養(yǎng)壓力的腫瘤區(qū)域的免疫監(jiān)視不足个从,突出了代謝調(diào)節(jié)在 TME 中的關鍵作用脉幢。事實上歪沃,最近的臨床試驗將檢查點抑制劑與針對葡萄糖、氨基酸和核苷酸代謝的代謝藥物結合起來嫌松。更好地了解免疫反應和新陳代謝之間的串擾將進一步有益于癌癥治療沪曙。

CAFs 已在多種類型的癌癥中表征,并且可以分為不同的亞型萎羔。雖然細胞系研究表明 CAFs 在宮頸癌細胞增殖中的支持作用液走,但這是第一項系統(tǒng)描述 CSCC 臨床樣本中促腫瘤 CAFs 的空間分布和生物學特性的研究。由于組織異質(zhì)性贾陷,CAFs 的標記基因因癌癥而異缘眶。在本研究中,ACTA2髓废、POSTN巷懈、ITGB4 和 FAP 是鑒定 CSCC 中 CAF 的足夠標記基因。其他基因如 KRT4慌洪、ITGA1砸喻、COL24A1 和 COL7A1 可能作為 CSCC 中 CAF 的互補標記基因,顯示成纖維細胞和癌性鱗狀細胞的細胞特性蒋譬。這些 CAF 對 TME 的異質(zhì)性做出了顯著貢獻割岛,TME 通過促進腫瘤生長、轉移和免疫逃避而表現(xiàn)出促腫瘤發(fā)生表型犯助。 CAFs+腫瘤增殖活躍但缺乏淋巴細胞浸潤癣漆。將這些免疫逃避的腫瘤暴露于免疫系統(tǒng)對于根除癌細胞至關重要。因此剂买,建議在治療晚期 CSCC 患者時結合針對 CAF 的療法惠爽。研究人員已嘗試使用抗體或抑制劑干擾 CAF 的激活、作用和正乘埠撸化過程婚肆,并正在進行多項臨床試驗。由于 CAF 高度表達的基因對正常組織也是必不可少的坐慰,因此需要密切監(jiān)測其功效和副作用较性。

總之,數(shù)據(jù)系統(tǒng)地證明了 CSCC 中病毒基因表達结胀、免疫反應和代謝的高度異質(zhì)性赞咙,表明針對多個生物學過程的聯(lián)合藥物或療法將是治療 CSCC 的更好實踐。 除了治療性疫苗和 ICB 療法外糟港,我們的研究結果表明攀操,對 CAF 和腫瘤代謝的干預可以補充目前對 CSCC 的治療。 對這些生物學方面的進一步研究可能有助于開發(fā)針對 CSCC 和其他 HPV 誘導的鱗狀細胞癌的新藥物或療法秸抚。

Methods

Identification of viral RNA

圖片.png

Analysis of the progression trajectory of cervical cancerous lesions with ST data

圖片.png

圖片.png

Multimodal intersection analysis (MIA)

圖片.png

生活很好速和,有你更好

最后編輯于
?著作權歸作者所有,轉載或內(nèi)容合作請聯(lián)系作者
禁止轉載歹垫,如需轉載請通過簡信或評論聯(lián)系作者。
  • 序言:七十年代末颠放,一起剝皮案震驚了整個濱河市县钥,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌慈迈,老刑警劉巖若贮,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異痒留,居然都是意外死亡谴麦,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進店門伸头,熙熙樓的掌柜王于貴愁眉苦臉地迎上來匾效,“玉大人,你說我怎么就攤上這事恤磷∶婧撸” “怎么了?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵扫步,是天一觀的道長魔策。 經(jīng)常有香客問我,道長河胎,這世上最難降的妖魔是什么闯袒? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任,我火速辦了婚禮游岳,結果婚禮上政敢,老公的妹妹穿的比我還像新娘。我一直安慰自己胚迫,他們只是感情好喷户,可當我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著访锻,像睡著了一般褪尝。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上朗若,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天恼五,我揣著相機與錄音昌罩,去河邊找鬼哭懈。 笑死,一個胖子當著我的面吹牛茎用,可吹牛的內(nèi)容都是我干的遣总。 我是一名探鬼主播睬罗,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼旭斥!你這毒婦竟也來了容达?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤垂券,失蹤者是張志新(化名)和其女友劉穎花盐,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體菇爪,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡算芯,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了凳宙。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片飘哨。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡敛纲,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情孽拷,我是刑警寧澤,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布瞒御,位于F島的核電站钾腺,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏饺汹。R本人自食惡果不足惜昔字,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望首繁。 院中可真熱鬧作郭,春花似錦、人聲如沸弦疮。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽胁塞。三九已至咏尝,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間啸罢,已是汗流浹背编检。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留扰才,地道東北人允懂。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓,卻偏偏與公主長得像衩匣,于是被迫代替她去往敵國和親蕾总。 傳聞我的和親對象是個殘疾皇子粥航,可洞房花燭夜當晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內(nèi)容