研究者利用DART-seq技術(shù),通過APOBEC1-YTH融合蛋白在活細(xì)胞中誘導(dǎo)m6A位點(diǎn)附近的C-to-U突變,以實(shí)現(xiàn)單核苷酸分辨率的m6A識(shí)別。研究發(fā)現(xiàn),這些C-to-U突變并非集中于m6A鄰近位置,而傾向于在距離m6A位點(diǎn)100 nt范圍內(nèi)聚集成簇。那么這些距離在10100nt范圍內(nèi)的突變能夠用于標(biāo)記m6A位點(diǎn),而不會(huì)干擾ONT DRS中m6A甲基化5-mer的電流信號(hào)?;?/span>這種內(nèi)源性標(biāo)記策略,研究團(tuán)隊(duì)利用在ONT DRS數(shù)據(jù)中提取了1,020,237個(gè)讀長(zhǎng)水平的5-mer?m6A信號(hào),并結(jié)合深度殘差神經(jīng)網(wǎng)絡(luò),成功開發(fā)了能夠在單條讀長(zhǎng)水平上識(shí)別m6A的工具m6Aiso。此外,研究者還采用了真邁生物的SURFSeq 5000平臺(tái)完成了A549細(xì)胞系的GLORI測(cè)序,為m6Aiso 在多個(gè)細(xì)胞系中的準(zhǔn)確性評(píng)估提供了有力的支撐,同時(shí)也反映了SURFSeq 5000平臺(tái)進(jìn)行二代高通量測(cè)序具有高度的可靠性。如下圖所示,研究團(tuán)隊(duì)在多個(gè)細(xì)胞系中比較了m6Aiso以及GLORI測(cè)序得到的m6A修飾水平,m6Aiso在HEK293T, HeLa, A549, mESC中都和GLORI高度一致,其中HEK293T, HeLa和 mESC的GLORI數(shù)據(jù)來源于原始的GLORI論文,而A549細(xì)胞的GLORI數(shù)據(jù)為團(tuán)隊(duì)自行建庫后采用SURFSeq 5000平臺(tái)完成的測(cè)序。

圖1 不同細(xì)胞系m6Aiso和GLORI測(cè)序得到的m6A修飾水平比較
利用工具m6Aiso,研究者闡述了相同m6A位點(diǎn)不同的異構(gòu)體上形成差異化m6A修飾的三種機(jī)制:1、EJC(外顯子連接復(fù)合體)對(duì)m6A抑制的差異,與之前的研究一致,發(fā)現(xiàn)EJC會(huì)抑制靠近外顯子邊界的m6A沉積,即在異構(gòu)體中外顯子邊界200 bp范圍內(nèi)的m6A甲基化水平較低;2、RNA結(jié)合蛋白對(duì)不同異構(gòu)體RNA的特異性結(jié)合,如TARBP2結(jié)合內(nèi)含子的鄰近外顯子上的m6A位點(diǎn)的甲基化水平顯著高于未結(jié)合TARBP2內(nèi)含子的鄰近外顯子;3、轉(zhuǎn)錄因子也會(huì)調(diào)控m6A,主要通過共轉(zhuǎn)錄的方式招募m6A 的甲基化轉(zhuǎn)移酶來修飾其轉(zhuǎn)錄的RNA,因此使用不同啟動(dòng)子的RNA由于轉(zhuǎn)錄因子不同而產(chǎn)生差異的m6A修飾。由于變化的m6A位點(diǎn)位于共同的3'UTR,而短讀長(zhǎng)測(cè)序無法確定每個(gè)亞型的具體5'區(qū)域,因此這一發(fā)現(xiàn)進(jìn)一步強(qiáng)調(diào)了在研究m6A動(dòng)態(tài)變化時(shí)使用長(zhǎng)讀長(zhǎng)測(cè)序的必要性。

圖2?m6A在單reads上的內(nèi)源性標(biāo)記和m6Aiso在單分子水平上的結(jié)果評(píng)價(jià)