新生物信息技術(shù)-EISA：解析轉(zhuǎn)錄和轉(zhuǎn)錄后過(guò)程對(duì)基因表達(dá)變化所起的作用

分裂分析(EISA)是一種生物信息技術(shù)，通過(guò)對(duì)RNA-Seq的分析，將轉(zhuǎn)錄和轉(zhuǎn)錄后效應(yīng)分離開(kāi)來(lái)，并將其映射到外顯子和內(nèi)含子?？蓽y(cè)量在不同條件下mRNA和前體RNA水平的變化，可以區(qū)分轉(zhuǎn)錄水平上調(diào)控的RNA或在轉(zhuǎn)錄后水平——miRNA成熟或正在降解。上皮間充質(zhì)轉(zhuǎn)化（EMT)是一個(gè)備受關(guān)注的課題，認(rèn)識(shí)其機(jī)制有助于研究腫瘤轉(zhuǎn)移和藥物敏感性等問(wèn)題。雖然EMT受許多miRNAs和轉(zhuǎn)錄因子的調(diào)節(jié)，然而miRNAs和轉(zhuǎn)錄因子的功能很難區(qū)分，一部分原因是從測(cè)序數(shù)據(jù)集中難以直接識(shí)別miRNA靶基因，以及難以確定miRNA效應(yīng)，本文主要闡述miRNAs參與EMT過(guò)程中如何在轉(zhuǎn)錄組中發(fā)揮作用，而運(yùn)用了一種生物信息技術(shù)-EISA（exon-intron split analysis）。（Nucleic Acids Research，IF=11.147 ，Extensive transcriptional responses are co-ordinated by microRNAs as revealed by Exon-Intron Split Analysis (EISA)）
研究背景：miR-200c是一個(gè)促間質(zhì)轉(zhuǎn)化上皮的 microRNA。
一、EISA成功區(qū)分轉(zhuǎn)錄和轉(zhuǎn)錄后基因調(diào)控

      我們初步比較人乳腺上皮細(xì)胞(HMLEs)的內(nèi)含子豐度的變化與活性和非活性基因的標(biāo)志的一致性。
圖A：一種通過(guò)TGF進(jìn)入間充質(zhì)狀態(tài)EMT細(xì)胞模型(MesHMLE)， miR-200C驅(qū)動(dòng)間質(zhì)向上皮轉(zhuǎn)化后，前5%的轉(zhuǎn)錄上調(diào)和下調(diào)基因 (ΔI)，由ChIP-Seq獲得活性染色質(zhì)(H3K4me3，H3K9-14ac，H3K27c-綠色)和非活性染色質(zhì)(H3K27me3-紅色)的標(biāo)記物的相對(duì)富集。在一個(gè)跨越25個(gè)基因滑動(dòng)窗口中紅色/綠色陰影顯示了測(cè)量的均值。(每第5個(gè)窗口繪圖)。
圖B：由EISA分析所有組蛋白標(biāo)記與轉(zhuǎn)錄水平上的基因表達(dá)相關(guān)(大ΔI)。相反，活性基因和非活性基因的標(biāo)記與轉(zhuǎn)錄后調(diào)控的基因(大ΔE-ΔI，小ΔI，圖1B)的相關(guān)性小.因此，EISA度量ΔI表明了基因轉(zhuǎn)錄的影響?；趫DA的ΔE-ΔI的評(píng)分。
圖C：對(duì)轉(zhuǎn)染miR-200C的mesHMLE細(xì)胞的RNA-seq數(shù)據(jù)進(jìn)行EISA分析，以評(píng)價(jià)miR-200C相對(duì)轉(zhuǎn)錄和轉(zhuǎn)錄后效應(yīng)。基于ΔE:ΔI， 轉(zhuǎn)錄后miR-200 C靶點(diǎn)的測(cè)量值負(fù)值。
圖D：由ΔE-ΔI、ΔI、ΔE確定基因表達(dá)的累積分布（fold change），來(lái)預(yù)測(cè)響應(yīng)miR-200C的靶位點(diǎn)。
二、miR-200c轉(zhuǎn)錄后調(diào)節(jié)EGF信號(hào)網(wǎng)絡(luò)

      鑒于EISA能有效區(qū)分轉(zhuǎn)錄和轉(zhuǎn)錄后調(diào)節(jié)，我們關(guān)注在轉(zhuǎn)錄后下調(diào)基因。對(duì)miR-200C和預(yù)測(cè)的靶點(diǎn)進(jìn)行分析，以及確定在HMLE系統(tǒng)中由miR-200C直接調(diào)控的基因和通路。
圖A：通過(guò)TargetScan或MicroT-CDS預(yù)測(cè)miR-200c的靶基因，（ΔE-ΔI，前1000位轉(zhuǎn)錄后負(fù)調(diào)控基因）。通過(guò)基因本體論研究靶基因之間的關(guān)系，發(fā)現(xiàn)在EMT中EGFR信號(hào)通路的foid enrichment是最高的。 
圖B：轉(zhuǎn)錄后下調(diào)基因的大多數(shù)是預(yù)測(cè)的高重疊和預(yù)測(cè)目標(biāo)的靶基因。通過(guò)the targetscan and microT-CDS algorithms預(yù)測(cè)，miR200c和miR200c目標(biāo)反應(yīng)的1000、500和200個(gè)轉(zhuǎn)錄后下調(diào)的基因。
圖C：在EGF信號(hào)網(wǎng)絡(luò)3’ UTR 熒光素酶報(bào)告基因測(cè)定miR-200c靶標(biāo)。
圖D：在EGFR下游的信號(hào)網(wǎng)絡(luò)中miR-200C的靶基因(藍(lán)色陰影)。KEGG pathway顯示Akt和ERK被激活。
圖E：用EGF處理Meshmle細(xì)胞，Western Blot檢測(cè)在miR-200c表達(dá)或不表達(dá)的情況下的EGFR、Akt、MEK和ERK(或總蛋白)的磷酸化。
三、在EMT/MET過(guò)程中Mir-200 c和TGFβ共同調(diào)節(jié)的轉(zhuǎn)錄反應(yīng)。

圖A：ΔE-ΔI (轉(zhuǎn)錄后),ΔI (轉(zhuǎn)錄)表示EISA定義的基因調(diào)控效應(yīng). 雖然miRNAs的直接作用是轉(zhuǎn)錄后產(chǎn)生的，例如miR-200 c對(duì)EGF信號(hào)的抑制，圖中分別顯示miR-200 c與TGFβ在EMT/MET過(guò)程中發(fā)生的基因表達(dá)變化，其發(fā)生ΔE-ΔI (轉(zhuǎn)錄后)與ΔI (轉(zhuǎn)錄)。
圖B：雖然已經(jīng)確定miRNAs和轉(zhuǎn)錄因子(TGFβ)之間的相互作用，但miRNA研究的一個(gè)主要焦點(diǎn)是通過(guò)直接作用于靶標(biāo)來(lái)理解其功能。然而，EISA表明，許多對(duì)miR-200C的反應(yīng)是間接的(ΔI)。為了探討這些轉(zhuǎn)錄變化的意義，基因本體論(分子功能)分析了在MEHMLE細(xì)胞對(duì)miR-200反應(yīng)和在HMGLE細(xì)胞中對(duì)TGF反應(yīng)的轉(zhuǎn)錄水平上調(diào)和下調(diào)基因集。
圖C：在轉(zhuǎn)錄上調(diào)的基因中，例如ZEB 1等，間質(zhì)促進(jìn)轉(zhuǎn)錄因子的結(jié)合位點(diǎn)富集在ENCODE ChIP- Seq數(shù)據(jù)集。每個(gè)條形表示單個(gè)ChIP- Seq數(shù)據(jù)集。
圖D：相對(duì)于miR-200c(ΔE-ΔI)轉(zhuǎn)錄后下調(diào)的相對(duì)程度，顯示了TFs介導(dǎo)下游轉(zhuǎn)錄變化。轉(zhuǎn)錄后被miR-200 c下調(diào)且在它們的3’端內(nèi)具有miR-200 c結(jié)合位點(diǎn)的TFs。
圖E：TF靶基因是從 ChIP- Seq數(shù)據(jù)中獲取的，然然后與根據(jù)miR-200轉(zhuǎn)錄上調(diào)或下調(diào)的基因交叉引用。
四、功能顯著的轉(zhuǎn)錄調(diào)控是miRNAs的一個(gè)廣泛特征。

圖A-B：為了研究miRNA干擾中是否與轉(zhuǎn)錄反應(yīng)有關(guān)，我們對(duì)EISA進(jìn)行了多個(gè)數(shù)據(jù)集的分析，發(fā)現(xiàn)在所檢測(cè)的8個(gè)轉(zhuǎn)錄反應(yīng)中，27%到73%的基因是共同的。表達(dá)ZEB 1或多個(gè)miRNAs的轉(zhuǎn)錄和轉(zhuǎn)錄后機(jī)制的相對(duì)貢獻(xiàn)用散點(diǎn)圖(B)( ΔE-ΔI：ΔI)
圖C：無(wú)論這是轉(zhuǎn)錄上調(diào)基因還是下調(diào)基因，為每個(gè)選定的miRNA進(jìn)行10個(gè)基因 fold enrichment。