1. 從“表觀”到“本質(zhì)”為什么DNA甲基化研究如此重要如果你在生物醫(yī)學領(lǐng)域待過一陣子無論是做腫瘤研究、發(fā)育生物學還是探索衰老與神經(jīng)退行性疾病大概率都繞不開“DNA甲基化”這個詞。它就像一個無處不在的“化學標簽”在不改變DNA序列本身的情況下深刻地調(diào)控著基因的表達。簡單來說你可以把它想象成一本寫在基因上的“使用說明書”哪些章節(jié)基因可以大聲朗讀高表達哪些章節(jié)需要保持沉默低表達或不表達很大程度上就由這些甲基化“標簽”來決定。我剛開始接觸這個領(lǐng)域時也常常被各種縮寫和概念搞得暈頭轉(zhuǎn)向CpG島、全基因組甲基化、啟動子區(qū)高甲基化、印記基因……感覺每個詞都懂但連起來就不知道研究該從哪里下手。更讓人頭疼的是實驗做出來了數(shù)據(jù)拿到了但面對海量的甲基化位點信息如何解讀其背后的生物學意義如何將它與具體的疾病表型聯(lián)系起來往往是比實驗本身更大的挑戰(zhàn)。這篇文章我就想結(jié)合自己這些年踩過的坑和積累的經(jīng)驗和你系統(tǒng)地聊聊DNA甲基化研究中最核心的十個問題。這不僅僅是概念梳理更是實戰(zhàn)指南。我們會從“為什么要研究它”開始一步步深入到“如何設計實驗”、“如何分析數(shù)據(jù)”最后探討“如何將發(fā)現(xiàn)轉(zhuǎn)化為真正的生物學洞見”。無論你是剛?cè)腴T的研究生還是正在設計課題的博士后希望這些內(nèi)容能幫你理清思路避開那些我當年走過的彎路。2. 核心概念掃盲甲基化到底“標記”了什么在深入具體問題之前我們必須先統(tǒng)一“語言”。DNA甲基化最經(jīng)典的形式就是在胞嘧啶C的第5位碳原子上添加一個甲基基團-CH3形成5-甲基胞嘧啶5mC。這個過程主要發(fā)生在DNA序列中“CpG”二核苷酸位點上也就是胞嘧啶后面緊跟著一個鳥嘌呤G的情況。注意這里說的“主要”非常關(guān)鍵。隨著技術(shù)發(fā)展我們現(xiàn)在知道甲基化也存在于非CpG上下文如CHH, CHG其中H代表A、T或C尤其在植物和哺乳動物的胚胎干細胞、神經(jīng)元中較為顯著。但在大多數(shù)哺乳動物體細胞研究中CpG位點的甲基化仍是焦點。2.1 CpG島基因調(diào)控的“戰(zhàn)略要地”大約60%的人類基因啟動子區(qū)域富含CpG序列這些區(qū)域被稱為“CpG島”。你可以把它們想象成基因開關(guān)附近的“控制面板”。在正常細胞中大多數(shù)CpG島是保持低甲基化或未甲基化狀態(tài)的這為轉(zhuǎn)錄因子結(jié)合和基因轉(zhuǎn)錄開啟提供了可能。相反基因體Gene Body區(qū)域的CpG位點則常常呈現(xiàn)高甲基化這種甲基化被認為與轉(zhuǎn)錄延伸、抑制基因內(nèi)部異常轉(zhuǎn)錄起始有關(guān)其功能更為復雜。2.2 甲基化的“書寫者”、“擦除者”與“閱讀者”理解甲基化的動態(tài)調(diào)控離不開這三類功能蛋白書寫者WritersDNA甲基轉(zhuǎn)移酶DNMTs。DNMT3A和DNMT3B負責建立全新的甲基化模式從頭甲基化而DNMT1則更傾向于在DNA復制后根據(jù)母鏈的甲基化模式在子鏈對應的CpG位點上添加甲基維持甲基化模式的遺傳維持甲基化。擦除者ErasersTET家族蛋白TET1/2/3。它們能將5mC逐步氧化成5-羥甲基胞嘧啶5hmC、5-甲?；奏?fC和5-羧基胞嘧啶5caC最終通過堿基切除修復途徑實現(xiàn)主動去甲基化。5hmC本身也被視為一種重要的表觀標記尤其在腦組織中含量豐富。閱讀者Readers含有甲基化CpG結(jié)合域MBD的蛋白如MeCP2、MBD1-4等。它們能特異性識別并結(jié)合甲基化的DNA進而招募染色質(zhì)重塑復合物、組蛋白修飾酶等共同構(gòu)成抑制性的染色質(zhì)環(huán)境導致基因沉默。這個“寫-讀-擦”的動態(tài)平衡構(gòu)成了表觀遺傳調(diào)控的基礎框架。任何一個環(huán)節(jié)的失調(diào)都可能引發(fā)基因表達程序的混亂進而導致疾病。3. 技術(shù)選型困境我該用哪種方法檢測甲基化這是實操中第一個攔路虎。方法太多各有優(yōu)劣選錯了不僅浪費經(jīng)費更可能得到有偏差甚至錯誤的結(jié)論。下面這個表格對比了最主流的幾種技術(shù)你可以根據(jù)你的科學問題和樣本情況對號入座。方法原理簡述分辨率覆蓋范圍優(yōu)點缺點適用場景全基因組亞硫酸氫鹽測序WGBS用亞硫酸氫鹽處理DNA將未甲基化的C轉(zhuǎn)化為U測序為T甲基化的C不變通過測序比對計算每個CpG位點的甲基化水平。單堿基全基因組幾乎所有CpG位點黃金標準無偏倚可發(fā)現(xiàn)新的甲基化區(qū)域。成本高數(shù)據(jù)量大分析復雜對DNA質(zhì)量要求高無法區(qū)分5mC和5hmC。探索性研究需要全基因組甲基化圖譜如胚胎發(fā)育、細胞分化過程。簡化代表性亞硫酸氫鹽測序RRBS使用限制性內(nèi)切酶如MspI富集CpG密集區(qū)域如啟動子區(qū)再進行亞硫酸氫鹽處理和測序。單堿基覆蓋約1-3百萬個CpG位點側(cè)重CpG島和啟動子。成本低于WGBS針對性強數(shù)據(jù)量更易管理。覆蓋不全會丟失CpG島外的調(diào)控區(qū)域酶切效率影響結(jié)果。大樣本量的病例-對照研究專注于啟動子/CpG島甲基化差異。甲基化DNA免疫共沉淀測序MeDIP-seq用抗5mC抗體富集甲基化的DNA片段然后進行測序。約100-300 bp全基因組范圍但偏向高甲基化區(qū)域。成本相對較低對DNA量要求不高適合研究高甲基化區(qū)域。分辨率低無法精確定位到單個CpG定量是半定量的抗體特異性可能有問題?？焖俸Y查基因組范圍的高甲基化變化如腫瘤樣本的異常高甲基化區(qū)域。甲基化芯片如Illumina EPIC基于探針雜交檢測預先設計好的約85萬個CpG位點。單堿基覆蓋約85萬個精心挑選的CpG位點包括ENCODE、FANTOM5等數(shù)據(jù)庫中的調(diào)控區(qū)域。通量高成本低技術(shù)成熟穩(wěn)定數(shù)據(jù)分析流程標準化。只能檢測預設位點無法發(fā)現(xiàn)新位點動態(tài)范圍可能受限。大規(guī)模流行病學隊列研究、臨床樣本回顧性分析、需要快速獲得可靠數(shù)據(jù)的項目。靶向甲基化測序如靶向Panel通過多重PCR或探針捕獲富集特定基因或區(qū)域如某個信號通路的所有基因啟動子再進行亞硫酸氫鹽測序。單堿基自定義通常幾十到幾百個基因區(qū)域。成本最低深度極高適合低頻變異檢測可用于FFPE等降解樣本。目標區(qū)域外一無所知panel設計依賴前期知識。臨床驗證、液體活檢ctDNA甲基化標志物、已知候選基因的深度研究。我的選型心得“探索未知”用WGBS “驗證已知”用芯片或靶向。如果你的研究是開創(chuàng)性的比如研究一個全新細胞類型的甲基化圖譜WGBS是不二之選。如果你是基于已有線索如GWAS發(fā)現(xiàn)的位點、已知的腫瘤抑制基因進行大樣本驗證那么850K芯片性價比最高。樣本量和預算決定上限。對于動輒上千樣本的隊列研究芯片幾乎是唯一可行的選擇。RRBS是一個不錯的折中方案但務必確認你關(guān)心的區(qū)域能被MspI酶有效覆蓋。別忘了樣本本身。如果是福爾馬林固定石蠟包埋FFPE的組織樣本DNA降解嚴重WGBS和RRBS可能失敗率很高此時靶向測序或經(jīng)過特殊優(yōu)化的芯片方案是更穩(wěn)妥的選擇。考慮下游分析需求。如果你需要做甲基化數(shù)量性狀位點mQTL分析或者與染色質(zhì)可及性ATAC-seq、組蛋白修飾ChIP-seq等多組學數(shù)據(jù)整合那么單堿基分辨率的WGBS或RRBS數(shù)據(jù)能提供更大的靈活性。4. 實驗設計與樣本準備中的“隱形陷阱”確定了方法下一步就是設計實驗和準備樣本。這里面的坑往往在數(shù)據(jù)分析階段才會暴露到時追悔莫及。4.1 對照組的設立不僅僅是“健康 vs 疾病”很多研究簡單地用健康人組織作為疾病的對照這存在很大問題。DNA甲基化具有高度的細胞類型特異性和組織特異性。比如從血液中檢測到的阿爾茨海默病相關(guān)甲基化變化可能僅僅反映了疾病狀態(tài)下免疫細胞組成比例的改變?nèi)缌馨图毎麥p少、中性粒細胞增多而非神經(jīng)元本身的表觀遺傳改變。這就是“細胞異質(zhì)性”造成的混淆。實操建議盡可能使用同一種細胞類型進行比較。如果使用混合組織如全血、大腦皮層強烈建議在分析時進行細胞類型比例校正。可以使用像“Houseman”算法這樣的參考數(shù)據(jù)集通過甲基化數(shù)據(jù)反推樣本中各種細胞的比例并將其作為協(xié)變量納入統(tǒng)計模型?；蛘咴诮?jīng)費和樣本允許的情況下通過流式分選或單細胞測序獲取純化的細胞群體。4.2 批次效應數(shù)據(jù)質(zhì)量的“頭號殺手”甲基化檢測尤其是芯片和測序?qū)嶒灢僮鞣浅Ｃ舾小２煌瑫r間點提取的DNA、不同批次處理的試劑、甚至不同操作員都可能引入系統(tǒng)性誤差這種誤差的強度有時會遠超真實的生物學差異。如果你把病例和對照分別放在兩個批次里處理那么你發(fā)現(xiàn)的任何“差異”都很可能是批次效應。實操建議隨機化在實驗流程中務必打亂病例和對照樣本的順序確保每個處理批次中都包含兩類樣本。技術(shù)重復在預算中留出空間對部分樣本如5-10%進行技術(shù)重復從同一樣本DNA開始獨立走完全部流程。這是評估數(shù)據(jù)重現(xiàn)性和批次效應大小的最直接方式。陽性對照使用商業(yè)化的標準甲基化DNA如來自不同細胞系的混合DNA在每批次中運行監(jiān)控批次間的一致性。數(shù)據(jù)分析校正在生物信息學分析中必須使用如ComBat在svaR包中等工具對批次效應進行校正。但切記這不能替代良好的實驗設計它只是補救措施。4.3 DNA質(zhì)量與量被低估的關(guān)鍵亞硫酸氫鹽處理對DNA是破壞性的會導致DNA片段化。因此輸入DNA的質(zhì)量至關(guān)重要。降解的DNAOD260/280異常電泳呈拖尾會嚴重影響轉(zhuǎn)化效率導致數(shù)據(jù)覆蓋度不均和偏差。量對于850K芯片通常要求至少250ng DNAWGBS要求更高通常需要微克級。一定要在樣本收集階段就規(guī)劃好寧多勿少。質(zhì)使用Qubit進行精確定量比Nanodrop準確并用凝膠電泳或安捷倫生物分析儀評估完整性RIN值或DV200值。5. 生物信息學分析流水線從原始數(shù)據(jù)到差異位點拿到原始數(shù)據(jù)芯片的IDAT文件或測序的fastq文件后真正的挑戰(zhàn)才開始。下面我以最常用的850K芯片和WGBS為例梳理核心分析步驟。5.1 850K芯片數(shù)據(jù)分析流程數(shù)據(jù)導入與質(zhì)控使用minfi或sesame等R包讀取IDAT文件。第一步就是質(zhì)控檢查檢測p值每個探針在每個樣本中的檢出信號是否顯著高于背景通常我們會剔除在任何樣本中檢測p值 0.01的探針。同時檢查樣本的甲基化β值密度分布圖異常樣本分布明顯偏離其他樣本需要排查原因。預處理與歸一化這是關(guān)鍵步驟旨在消除技術(shù)變異。包括背景校正扣除熒光背景信號。染料偏倚校正850K芯片使用兩種熒光染料需要校正其效率差異。亞硫酸氫鹽轉(zhuǎn)化效率校正通過內(nèi)參探針評估。探針類型校正Infinium技術(shù)有I型和II型兩種探針其化學性質(zhì)不同信號強度范圍不同必須校正如BMIQ方法。歸一化使樣本間數(shù)據(jù)分布具有可比性常用ssNoob、FunctionalNormalization等方法。批次效應校正如前所述使用ComBat。差異甲基化分析計算每個CpG位點在兩組間的甲基化水平差異。對于芯片數(shù)據(jù)每個位點的甲基化水平通常用β值表示甲基化信號強度 / 總信號強度范圍0-1。常用的差異分析工具包括limma適用于芯片非常強大、DSS或methylKit對測序和芯片都適用。你需要設定閾值通常包括差異倍數(shù)如Δβ 0.1或0.2和統(tǒng)計學顯著性如FDR校正后的p值 0.05。差異甲基化區(qū)域DMR分析單個CpG位點的差異可能不穩(wěn)定生物學功能往往由基因組上相鄰的多個差異位點構(gòu)成的區(qū)域共同行使。使用DMRcate、bumphunter等工具可以找出這些連續(xù)的DMRs這通常比單個位點更有生物學意義。5.2 WGBS數(shù)據(jù)分析流程WGBS分析更為復雜計算量巨大。原始數(shù)據(jù)質(zhì)控與比對使用FastQC檢查測序質(zhì)量用Trim Galore或Cutadapt進行接頭修剪和質(zhì)量修剪。然后使用專門為亞硫酸氫鹽處理數(shù)據(jù)設計的比對工具如Bismark或BS-Seeker2將 reads 比對到參考基因組。這些工具會進行“三重比對”以區(qū)分轉(zhuǎn)化后的序列。甲基化位點提取比對后使用Bismark的methylation_extractor功能逐個計算每個CpG位點的甲基化率。輸出文件通常包含每個C位點的覆蓋深度和甲基化C的計數(shù)。樣本間差異分析由于WGBS數(shù)據(jù)是計數(shù)型數(shù)據(jù)甲基化C數(shù) vs 總覆蓋C數(shù)更適合使用基于負二項分布的模型。DSS和methylKit是常用選擇。它們能考慮覆蓋深度的波動更準確地識別差異位點。平滑與DMR鑒定WGBS覆蓋度可能不均勻直接進行單點分析噪音較大。BSmooth或DSS的平滑功能可以先對每個樣本的甲基化水平在基因組上進行平滑處理然后再進行差異區(qū)域?qū)ふ疫@樣找到的DMR更可靠。一個常見的坑過度依賴默認參數(shù)。無論是質(zhì)控閾值、歸一化方法還是DMR尋找的窗口大小、閾值都沒有“放之四海而皆準”的標準。你需要根據(jù)自己數(shù)據(jù)的特點覆蓋深度、樣本數(shù)、生物學變異大小進行調(diào)試和驗證。例如在樣本量較小的情況下尋找DMR的統(tǒng)計效力不足可能更適合先關(guān)注差異最顯著的單個位點再通過功能注釋判斷其是否位于有意義的調(diào)控區(qū)域。6. 功能注釋與生物學意義挖掘從數(shù)據(jù)列表到故事找到了成百上千個差異甲基化位點或區(qū)域報告一個p值和Δβ值列表是遠遠不夠的。如何解讀它們這是將數(shù)據(jù)轉(zhuǎn)化為生物學發(fā)現(xiàn)的關(guān)鍵。6.1 位置位置還是位置首先看DMR落在基因組的什么位置。使用annotatr、ChIPseeker等R包或在線工具如GREAT將DMR注釋到最近的基因、以及相關(guān)的基因組功能元件上。啟動子區(qū)特別是CpG島這里的甲基化通常與基因沉默直接相關(guān)是最經(jīng)典的調(diào)控模式。高甲基化往往意味著基因表達下調(diào)。增強子區(qū)增強子的甲基化狀態(tài)同樣調(diào)控其活性。但增強子的鑒定本身更復雜需要結(jié)合組蛋白修飾如H3K27ac和染色質(zhì)可及性數(shù)據(jù)來精確定義?；蝮wGene Body功能復雜可能與轉(zhuǎn)錄延伸、選擇性剪接相關(guān)不能簡單地與表達抑制劃等號。絕緣子/邊界元件如CTCF結(jié)合位點。這些區(qū)域的甲基化可能影響染色質(zhì)拓撲結(jié)構(gòu)從而遠距離調(diào)控基因表達。6.2 通路與網(wǎng)絡分析將差異基因與DMR相關(guān)的基因列表導入通路富集分析工具如clusterProfiler進行GO和KEGG富集分析。這能告訴你這些甲基化變化的基因是否集中在某個特定的生物學過程、分子功能或信號通路上。例如在癌癥研究中你很可能看到“Wnt signaling pathway”、“Cell cycle”、“PI3K-Akt signaling pathway”等經(jīng)典通路被富集。更進一步可以構(gòu)建蛋白質(zhì)-蛋白質(zhì)相互作用PPI網(wǎng)絡使用STRING數(shù)據(jù)庫、Cytoscape軟件找出網(wǎng)絡中的核心樞紐Hub基因。這些基因可能扮演更關(guān)鍵的角色。6.3 整合多組學數(shù)據(jù)建立“甲基化-表達”關(guān)聯(lián)這是提升研究層次的關(guān)鍵。單獨看甲基化變化是孤立的如果能證明甲基化變化導致了基因表達的改變你的故事就完整多了。如果有配對的轉(zhuǎn)錄組數(shù)據(jù)RNA-seq這是最理想的情況。直接計算差異甲基化基因的mRNA表達變化進行相關(guān)性分析。注意啟動子高甲基化通常與表達負相關(guān)但并非絕對存在延遲效應或間接調(diào)控。如果沒有配對數(shù)據(jù)可以利用公共數(shù)據(jù)庫如TCGA中同類疾病的數(shù)據(jù)驗證在你發(fā)現(xiàn)的DMR相關(guān)基因上是否普遍存在表達變化。雖然證據(jù)力度稍弱但能提供重要的支持。共定位分析將你的DMR與公開的染色質(zhì)狀態(tài)圖譜如ENCODE的ChIP-seq數(shù)據(jù)、染色質(zhì)可及性數(shù)據(jù)ATAC-seq進行疊加看DMR是否落在活躍的啟動子或增強子區(qū)域這能增強其功能重要性。7. 因果推斷的挑戰(zhàn)是“因”還是“果”這是表觀遺傳學研究尤其是觀察性研究如病例-對照中最棘手的問題。我們檢測到的疾病相關(guān)甲基化差異究竟是導致疾病的“原因”驅(qū)動因素還是疾病發(fā)生發(fā)展過程中產(chǎn)生的“結(jié)果”伴隨現(xiàn)象或者是第三方因素如環(huán)境暴露、年齡、生活方式導致的“旁觀者”7.1 利用自然實驗和縱向研究設計孟德爾隨機化如果某個甲基化位點受遺傳變異甲基化數(shù)量性狀位點mQTL強烈影響我們可以利用這些遺傳變異作為工具變量來推斷甲基化對疾病的潛在因果效應。這需要大型的基因組和表觀基因組關(guān)聯(lián)研究數(shù)據(jù)?？v向隊列在疾病發(fā)生前多年收集生物樣本如血液發(fā)病后再進行檢測。如果某些甲基化變化在疾病臨床診斷前就已出現(xiàn)那么它作為“原因”或早期標志物的可能性就大大增加。例如很多癌癥早期診斷的研究都在朝這個方向努力。干預性研究在細胞或動物模型中通過基因編輯如CRISPR-dCas9-TET1/DNMT3A特異性改變目標位點的甲基化狀態(tài)然后觀察表型變化。這是證明因果關(guān)系的“金標準”實驗。7.2 區(qū)分驅(qū)動因素與乘客事件在腫瘤研究中這個概念尤為重要。腫瘤基因組中存在大量異常甲基化但只有少數(shù)位于關(guān)鍵基因調(diào)控區(qū)、并能促進細胞生長優(yōu)勢的甲基化事件才是“驅(qū)動因素”其余大部分可能是基因組不穩(wěn)定導致的隨機“乘客事件”。如何區(qū)分一個常用的思路是看該甲基化事件是否在多種腫瘤類型或大量獨立樣本中重復出現(xiàn)復現(xiàn)性以及它是否影響已知的癌癥相關(guān)基因或通路。8. 甲基化“時鐘”如何準確解讀生物學年齡近年來基于DNA甲基化水平構(gòu)建的“表觀遺傳時鐘”如Horvath時鐘、Hannum時鐘、PhenoAge、GrimAge非?；馃?。它們能高精度地預測個體的生物學年齡而生物學年齡與衰老相關(guān)疾病、死亡率密切相關(guān)。8.1 時鐘的構(gòu)建與應用這些時鐘本質(zhì)上是一組特定的CpG位點通過機器學習算法通常是彈性網(wǎng)絡回歸訓練使其甲基化水平加權(quán)組合后與實際年齡高度相關(guān)。使用時將你的樣本數(shù)據(jù)經(jīng)過與訓練集相同的歸一化處理代入公式即可計算出“表觀遺傳年齡”。年齡加速如果一個人的表觀遺傳年齡大于其實際年齡即年齡加速通常意味著其生物學衰老更快健康風險更高。臨床應用評估衰老干預措施如飲食、運動、藥物的效果預測年齡相關(guān)疾病如阿爾茨海默病、心血管疾病的風險甚至在法醫(yī)學中用于年齡推斷。8.2 使用時鐘的注意事項訓練集匹配不同的時鐘在不同人群、組織中的表現(xiàn)不同。Horvath時鐘是多組織時鐘適用性廣Hannum時鐘是基于血液開發(fā)的。選擇與你樣本組織最匹配的時鐘。技術(shù)平臺一致性時鐘的CpG位點必須在你使用的技術(shù)平臺如850K芯片上存在。確保你的數(shù)據(jù)預處理流程特別是歸一化方法與時鐘原始論文中描述的一致否則計算結(jié)果會偏差很大。解讀需謹慎年齡加速是一個綜合指標它反映了多種內(nèi)源遺傳和外源環(huán)境、生活方式因素對衰老的累積影響。它指向風險但本身不是特定疾病的診斷標志物。需要結(jié)合其他臨床指標和生物學知識進行解讀。9. 從實驗室到臨床甲基化作為生物標志物的前景與挑戰(zhàn)基于甲基化的液體活檢如檢測循環(huán)腫瘤DNA的甲基化標志物是當前精準醫(yī)療的熱點。相比基因突變甲基化信號在早期腫瘤中可能更早、更普遍地出現(xiàn)且具有組織溯源特性不同組織有特異的甲基化模式。9.1 標志物開發(fā)的流程發(fā)現(xiàn)階段使用WGBS或甲基化芯片在少量組織樣本如癌與癌旁中篩選出差異巨大的甲基化區(qū)域。驗證與優(yōu)化在獨立的、更大的樣本集中使用靶向甲基化測序如甲基化特異性PCR、數(shù)字PCR驗證這些區(qū)域。進一步縮小范圍找到最敏感、特異的標志物組合Panel。臨床驗證在前瞻性隊列中使用血液等液體活檢樣本評估該Panel對疾病的早期檢測、預后預測或療效監(jiān)測的效能。這是最耗時耗力但決定成敗的階段。9.2 面臨的挑戰(zhàn)信號微弱早期腫瘤或微小殘留病灶釋放到血液中的ctDNA量極少甲基化信號容易被背景噪音淹沒。需要超高深度的測序或極其靈敏的檢測技術(shù)。異質(zhì)性腫瘤內(nèi)部存在甲基化異質(zhì)性單個活檢組織發(fā)現(xiàn)的標志物可能無法代表所有腫瘤細胞。液體活檢理論上能捕獲這種異質(zhì)性但分析更復雜。標準化與合規(guī)要將一個甲基化檢測推向臨床需要經(jīng)過嚴格的分析驗證檢測性能和臨床驗證臨床效用并最終獲得監(jiān)管機構(gòu)如FDA、NMPA的批準。這個過程漫長且昂貴。10. 單細胞甲基化測序看見細胞群體的“暗物質(zhì)”傳統(tǒng)批量測序得到的是細胞群體的平均信號會掩蓋細胞間的異質(zhì)性。而細胞異質(zhì)性正是發(fā)育、免疫、腫瘤微環(huán)境等核心生物學過程的關(guān)鍵。scBS-seq、scRRBS、以及基于液滴的scNOME-seq等技術(shù)讓我們能在單細胞分辨率下觀察DNA甲基化。10.1 技術(shù)特點與數(shù)據(jù)特性單細胞甲基化數(shù)據(jù)極其稀疏。一個細胞僅能檢測到全基因組中很小一部分CpG位點覆蓋率約5-20%且不同細胞覆蓋的位點不同。這導致無法直接進行細胞間的位點對位點比較。分析的核心是降維、聚類和細胞分型將甲基化模式相似的細胞聚在一起推斷細胞狀態(tài)或類型。需要整合其他模態(tài)數(shù)據(jù)單獨的單細胞甲基化數(shù)據(jù)解釋性有限。目前的主流是與單細胞轉(zhuǎn)錄組scRNA-seq或染色質(zhì)可及性scATAC-seq數(shù)據(jù)進行多組學整合分析例如通過“加權(quán)最近鄰”WNN方法綜合兩種模態(tài)的信息來更準確地定義細胞類型和狀態(tài)。10.2 能回答的新問題細胞命運決定的動態(tài)圖譜在胚胎發(fā)育或器官類分化過程中追蹤單個細胞甲基化模式的連續(xù)變化揭示表觀遺傳重編程的軌跡。腫瘤微環(huán)境的表觀遺傳景觀解析腫瘤中癌細胞、免疫細胞、成纖維細胞等不同組分各自獨特的甲基化特征以及它們之間的相互作用。罕見細胞類型的鑒定發(fā)現(xiàn)那些在批量測序中被淹沒的、具有獨特甲基化模式的稀有細胞亞群。我的體會是單細胞甲基化研究目前仍處于方法學和資源積累階段計算分析門檻很高。對于大多數(shù)研究者更現(xiàn)實的路徑可能是先利用已發(fā)表的單細胞甲基化圖譜數(shù)據(jù)庫作為參考來幫助解構(gòu)自己批量樣本數(shù)據(jù)中隱含的細胞組成異質(zhì)性。當你真正需要探索一個完全未知的、高度異質(zhì)的細胞群體內(nèi)部動態(tài)時再考慮投入這項昂貴而復雜的技術(shù)。回顧這十個問題它們貫穿了DNA甲基化研究從概念到設計、從實驗到分析、從數(shù)據(jù)到生物學、從基礎到臨床轉(zhuǎn)化的全鏈條。這個領(lǐng)域技術(shù)迭代飛快新的修飾如6mA、新的檢測方法層出不窮。但萬變不離其宗清晰的研究邏輯、嚴謹?shù)膶嶒炘O計、對數(shù)據(jù)本質(zhì)的深刻理解以及將分子變化與宏觀表型緊密聯(lián)系的生物學思維始終是做出扎實工作的基石。希望這篇長文能成為你手邊的一份實用地圖當你在甲基化研究的森林中探索時能幫你辨明方向少走些彎路。