Science丨揭示L1逆轉座如何在腫瘤早期制造染色體互易易位

撰文｜雪月

L1逆轉座（ LINE-1 L1 retrotransposition ） 在腫瘤基因組中普遍存在，并可通過 “ 復制 - 粘貼 ” 式插入引發插入突變及大片段缺失等結構變異。然而，由于短讀長測序對重復序列與復雜斷點的解析能力有限， L1 是否以及如何促進平衡型結構變異（如互易易位）仍缺乏系統性證據。

來自西班牙圣地亞哥－德孔波斯特拉大學Jose M. C.Tubio團隊在 Science 上發表題為 Concurrent L1 retrotransposition events promote reciprocal translocations in human tumorigenesis 的文章。 本研究以“逆轉座高活性腫瘤”作為富集模型，結合長讀長測序與專門算法，旨在在序列層面直接重建L1相關斷點結構，并據此推斷其導致的結構重排類型、頻率與發生時間。

作者對 137 例腫瘤進行淺層 Illumina WGS 初篩以評估體細胞逆轉座負荷，并選取其中 10 例逆轉座事件數 >100 的腫瘤作為 “golden set” 。隨后對這 10 例腫瘤及配對鄰近正常組織進行 ONT 長讀長全基因組測序 。 為解析長讀長數據中的逆轉座插入與結構后果，作者開發 MEIGA ，可識別并重建 7 類逆轉座插入，并檢測 “ 由新獲得逆轉座插入橋接兩個遠端斷點 ” 的逆轉座介導基因組重排 。

方法學可靠性通過模擬數據基準測試進行評估，同時在插入負荷最高的樣本中引入 PacBio HiFi 與 Illumina 數據進行交叉驗證，并對部分 “MEIGA 獨有 ” 事件做 IGV 目檢確認。對關鍵互易易位，作者進一步采用 long-range PCR 排查斷點處是否存在 “ 預先存在的 L1 中間體 ” ，并在細胞系 NCI-H2009 中結合 FISH 與 Micro-C 等正交手段驗證跨染色體結構。分子計時方面，作者結合 VAF 、局部拷貝數與 WGD 定時框架，將插入 /RT-RG 劃分為克隆早期與晚期等類別，并以 “ 時鐘樣 ” 突變推斷 WGD 與診斷時間的相對關系。

在 10 例高活性腫瘤中， MEIGA 共鑒定 6,418 個體細胞逆轉座事件（中位數 622.5/ 例；范圍 296–1,311 ）。值得注意的是，與長讀長結果相比，短讀長流程漏檢了 43% 的經典插入事件（ 2,697/6,266 ），提示在高重復序列背景下，短讀長對插入負荷的估計可能存在系統性低估。在多平臺驗證中，插入最多的樣本中 92.5% 的典型插入可被至少一種獨立平臺 / 流程支持，同時對 “MEIGA 獨有 ” 事件的目檢也支持其為真實體細胞插入，增強了后續結構推斷的可信度。

在 golden set 中，作者識別出 152 個 RT-RG （約占全部體細胞逆轉座事件的 ~2% ）。 RT-RG 的類型譜顯示，缺失仍為最常見類型（ 43% ），但跨染色體重排占 30% ，倒位與重復分別占 18% 與 9% ，提示 L1 活動不僅造成局部插入或缺失，也可參與更大尺度的基因組結構重塑。在跨染色體事件中，約 44% （ 20/45 ）參與形成 13 個互易易位，這一 “ 平衡型交換 ” 在以往短讀長研究中不易捕捉，因而可能長期被低估。

作者對互易易位斷點進行序列層面精細解析，發現多例易位的兩條衍生染色體連接處分別包含不同的 L1 插入片段，并呈現 poly(A) 尾、 L1 內切酶切割基序及短微同源等特征；此外，斷點附近靶位點特征在兩條非同源染色體間呈 “ 互換 ” 樣式，符合 “ 兩次獨立插入在時間上并發 / 緊鄰發生后被共同解析 ” 的模型。

為檢驗替代解釋（例如斷點處存在預先形成的 L1 中間體并隨后發生重組），作者采用 long-range PCR 檢測斷點區域，未發現支持 “ 預先存在中間體 ” 的證據，從而進一步反對 “ 先插入后同源重組 ” 的路徑。

統計上， 13 個互易易位中有 10 個更符合 “ 兩次插入共同參與 ” 的結構模式，并且 long-range PCR 未檢測到 pre-existing 元件，整體支持 “ 并發或緊密連續發生 ” 的機制解釋。

為提供獨立體系的直接證據，作者進一步在肺腺癌細胞系 NCI-H2009 中進行 ONT+MEIGA 分析并發現 chr3q–chr6q 的互易易位由兩個獨立 L1 插入橋接；隨后 FISH 觀察到符合預期的易位構型， Micro-C 也在預測斷點處顯示跨染色體互作信號，從細胞遺傳學與三維基因組層面共同驗證了該互易易位的真實存在及其斷點結構。

基于納米孔數據的甲基化分析顯示，體細胞逆轉座事件數與 L1 啟動子甲基化水平呈顯著負相關（ r = ?0.537 ， P = 0.015 ），并且新獲得插入的 L1 往往也呈低甲基化，提示 L1 去甲基化 / 解抑制可能與逆轉座爆發相關。在 20 個腫瘤樣本的分析中， RT-RG 數與插入數呈正相關（ r = 0.66 ， P = 0.001 ），支持 “ 逆轉座活躍度越高，發生結構重排的機會越大 ” 的總體趨勢。同時，在 RT-RG 負擔最高的樣本中觀察到 Fanconi anemia （ FA ）通路關鍵基因的雙等位突變，為復制叉保護 / 修復通路影響 L1 相關結構重排風險提供了線索性關聯。在可定時的典型插入中，約 64.1% 被歸類為克隆早期事件，提示逆轉座活動在腫瘤克隆擴張早期已廣泛發生。

結合 “ 時鐘樣 ” 突變推斷， WGD 的發生中位時間約為診斷前 4.77 年，且所有分析腫瘤在 WGD 前均已發生大量逆轉座事件。對 RT-RG 而言，可置信定時的事件中約 3/4 為早期克 隆事件；尤其是 L1 介導的互易易位中，有 9 個發生在 WGD 之前，最早可推斷至診斷前 8.87 年以上，強調 L1 活動可能在腫瘤演化早期就參與推動染色體結構重排與不穩定。

這項工作將L1的角色從“插入突變源”拓展為“結構重排生成機制”，為理解腫瘤染色體不穩定的起源與早期演化提供了新的可檢驗框架。

原文鏈接：https://www.science.org/doi/10.1126/science.aee4513

制版人： 十一

學術合作組織

（*排名不分先后）

戰略合作伙伴

（*排名不分先后）

轉載須知

【非原創文章】本文著作權歸文章作者所有，歡迎個人轉發分享，未經作者的允許禁止轉載，作者擁有所有法定權利，違者必究。

BioArt

Med

Plants

人才招聘

近期直播推薦

點擊主頁推薦活動

關注更多最新活動！