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