MeRIP-seq揭示m6A修飾在肺動脈高壓發病機制中的潛在作用和新治療靶點
肺動脈高壓（pulmonary arterial hypertension， PAH）是一種不可逆的心血管疾病，其特征是肺血管阻力進行性增加，最終導致肺動脈高壓和右心衰竭。PAH主要成因包括血管收縮、肺血管重塑和細胞外基質（ECM）沉積。肺血管重塑主要由肺動脈平滑肌細胞（PASMCs）的過度增殖和抗凋亡能力引起。PAH是一種涉及遺傳、環境和表觀遺傳等多種因素的復雜疾病。表觀遺傳變化在血管重塑中的發病機制已被證實，但關于逆轉PAH的特定表觀遺傳靶點的了解仍然有限。
N6-甲基腺苷（m6A）是一種豐富的內源性化學修飾，對RNA的剪接、翻譯、定位和穩定性起著關鍵作用。m6A的動態和可逆性由甲基轉移酶（writers）、去甲基化酶（erasers）和不同的m6A結合蛋白（readers）調控。已有研究表明，m6A調控因子的表達失衡與多種生物過程功能障礙有關，包括凋亡/增殖、發育異常、腫瘤進展、自我更新能力降低和免疫系統異常。盡管有證據表明m6A酶失調參與了PAH發病機制，但m6A調控因子在PAH中的具體作用尚未明確。
2024年2月4日，中南大學湘雅二醫院鳳一路博士為第一作者、宋潔為通訊作者在《Biomedicines》雜志發表題為“Transcriptome-Wide N6-Methyladenosine Alternations in Pulmonary Arteries of Monocrotaline-Induced Pulmonary Arterial Hypertension in Rats and Novel Therapeutic Targets”的研究論文，文章通過甲基化RNA免疫沉淀測序（MeRIP-seq）和RNA轉錄組測序（RNA-seq）方法，分析單克隆素（MCT）誘導的PAH大鼠模型和對照組肺動脈組織中的mRNA表達和m6A位點變化，通過GO和KEGG通路富集分析進一步分析差異表達基因的功能。探討了m6A調控因子表達及其在調控PASMC生物學行為中的作用，以揭示潛在的分子機制。
研究摘要：
本研究利用MeRIP-seq和RNA-seq研究對照組和野百合堿（monocrotaline）誘導的PAH大鼠肺動脈組織中m6A和RNA表達差異。使用GO和KEGG進行功能富集分析，研究了差異表達功能。為了篩選候選m6A相關基因，使用STRING和Metascape數據庫構建了蛋白質-蛋白質互作（protein-protein interaction，PPI）網絡，并通過對候選相關基因表達進行實時PCR驗證。使用免疫組化染色、免疫熒光和Western blot技術進一步研究了m6A調控因子的表達水平，并對大鼠肺動脈平滑肌細胞（PASMC）進行增殖實驗。研究共鑒定出42個差異表達基因，這些基因的m6A甲基化水平表現出增加或減少（高甲基化或低甲基化）。主要與細胞外基質結構、MAPK和PI3K/AKT通路相關。在PAH組中檢測到候選基因著絲粒蛋白F（CENPF）表達增加。研究首次鑒定出一種富含亮氨酸的五肽重復序列（LRPPRC）的m6A reader，該reader在PAH大鼠模型中表達下調。體外實驗中，siRNA介導的Lrpprc下調導致大鼠PASMC增殖增強和Cenpf mRNA表達升高。本研究結果揭示了PAH大鼠肺動脈中轉錄組范圍的m6A修飾圖譜變化和相關調控機制，可能為未來治療策略提供新的靶點。
研究思路：
研究結果
（1）PAH大鼠模型的m6A基本特征
圖1： PAH大鼠模型的m6A基本特征。組織來自肺動脈，每組n=3。
A. 維恩圖顯示了兩組中m6A peaks（左）和基因（右）的重疊。
B. 餅圖顯示了轉錄本四個不重疊區域（（5' UTR、CDS、3' UTR和其他區域）的m6A peaks百分比。
C. 密度曲線顯示m6A peaks在轉錄本中的分布。轉錄本分為三部分，即5'UTR、CDS和3'UTR。
D. 兩組不同數量m6A peaks的基因比例。
E. 兩組RRACH motif情況。
m6A：N6-甲基腺苷；PAH：肺動脈高壓；MCT：野百合堿；UTR：非翻譯區；CDS：編碼序列。
（2）PAH中的差異m6A分析
B. 差異m6A peaks中富集的前五個motif。
C. 具有差異m6A的上調m6A標記的轉錄本GO分析。
D. 具有差異m6A的下調m6A標記的轉錄本GO分析。
E. 具有差異m6A的上調m6A標記的轉錄本KEGG分析。
F. 具有差異m6A的下調m6A標記的轉錄本KEGG分析。
FC：倍數變化；GO：基因本體論；KEGG：京都基因與基因組百科全書。
（3）m6A修飾基因的轉錄譜
B. 具有顯著變化的peak差異表達基因的四象限圖。
C. 總體m6A與mRNA表達水平呈正相關（r=0.1517793，p＜0.05）。
D-E. GO和KEGG分析分別具有顯著變化peak的差異表達基因。
F-G. 兩個數據集的DEG分布的火山圖。紅點代表上調DEG，藍點代表下調DEG，黑點代表非DEG。
H. 在兩個數據集中同時上調的DEG。
I. 在兩個數據集中同時下調的DEG。DEG：差異表達基因。
（4）PPI網絡建立和候選基因鑒定
D. MCT組和對照組（每組n=6）肺組織中Cenpf的相對mRNA表達，數據代表平均值±SEM，并使用Student t檢驗比較兩組組間差異。
E. 血小板衍生生長因子BB（PDGF-BB）刺激（30 ng/mL）下，大鼠PASMC中Cenpf mRNA相對表達水平（每組n=3），以α-Tublin作為內部對照歸一化。
F. 人類PASMCs在PDGF-BB刺激下的CENPF mRNA相對表達水平（每組n=3），以GAPDH作為內部對照歸一化。
G. 用siRNA轉染人PASMC 48h后，經PDGF-BB（30 ng/mL）再處理24h的CENPF mRNA相對表達水平（每組n=3）。
H-I. 以α-Tublin作為內部對照的PCNA蛋白水平的代表性Western blot和定量分析（每組n=3）。
J. 通過IGV觀察到Cenpf mRNA轉錄本上的m6A水平。
PPI：蛋白質-蛋白質互作；SEM：平均值的標準誤差；PASMCs：肺動脈平滑肌細胞；PDGF-BB：血小板衍生生長因子BB。
（5）肺動脈高壓（PAH）的肺動脈中，leucine-rich pentatricopeptide repeat-containing (LRPPRC) 蛋白表達水平下降
圖5： PAH患者的肺動脈LRPPRC表達降低。
A. 兩組（PAH和對照組）中23個m6A調控因子的表達情況熱圖。
B. 肺組織中LRPPRC的免疫組化染色（每組n=6），條形圖=100µm（左圖），放大的肺動脈（右圖）。
C. 肺組織中LRPPRC的免疫熒光染色（每組n=6），血管染色圖像（綠色）、目標蛋白染色圖像（紅色）、細胞核染色圖像（藍色），合并圖像（橙色），條形圖=50µm。
D-E. MCT和對照組肺組織中LRPPRC和α-Tubulin的代表性Western blot和定量分析（每組n=6）；數據表示為均值±標準誤（SEM），使用學生t檢驗比較兩組。
F-G. 大鼠PASMCs在PDGF-BB刺激下的LRPPRC水平的代表性Western blot和定量分析（30 ng/mL），以α-Tubulin作為內部對照（每組n=3）。
H. 在PDGF-BB刺激下（30 ng/mL），人類PASMCs中LRPPRC相對mRNA表達水平，以GAPDH作為內部對照（每組n=3）。
（6）LRPPRC下調促進PASMC增殖和Cenpf表達
B-C. 代表性的Western blot和PCNA水平的定量分析，以α-Tubulin作為內部對照歸一化（每組n=3）。
D. EdU摻入實驗評估si-Lrpprc處理的PASMCs在PDGF-BB刺激下的細胞增殖能力。條形圖=100µm。目標細胞為紅色。
E.計算EdU染色的細胞比率。
F. Cenpf水平的相對mRNA表達（每組n=3）。
參考文獻：
Feng Y, Yu Z, Tang M, Li J, Peng B, Juaiti M, Tang Y, Liang B, Ouyang M, Liu Q, Song J. Transcriptome-Wide N6-Methyladenosine Alternations in Pulmonary Arteries of Monocrotaline-Induced Pulmonary Arterial Hypertension in Rats and Novel Therapeutic Targets. Biomedicines. 2024 Feb 4;12(2) pii: biomedicines12020364. doi: 10.3390/biomedicines12020364. PubMed PMID: 38397966.
N6-甲基腺苷（m6A）是一種豐富的內源性化學修飾，對RNA的剪接、翻譯、定位和穩定性起著關鍵作用。m6A的動態和可逆性由甲基轉移酶（writers）、去甲基化酶（erasers）和不同的m6A結合蛋白（readers）調控。已有研究表明，m6A調控因子的表達失衡與多種生物過程功能障礙有關，包括凋亡/增殖、發育異常、腫瘤進展、自我更新能力降低和免疫系統異常。盡管有證據表明m6A酶失調參與了PAH發病機制，但m6A調控因子在PAH中的具體作用尚未明確。
2024年2月4日，中南大學湘雅二醫院鳳一路博士為第一作者、宋潔為通訊作者在《Biomedicines》雜志發表題為“Transcriptome-Wide N6-Methyladenosine Alternations in Pulmonary Arteries of Monocrotaline-Induced Pulmonary Arterial Hypertension in Rats and Novel Therapeutic Targets”的研究論文，文章通過甲基化RNA免疫沉淀測序（MeRIP-seq）和RNA轉錄組測序（RNA-seq）方法，分析單克隆素（MCT）誘導的PAH大鼠模型和對照組肺動脈組織中的mRNA表達和m6A位點變化，通過GO和KEGG通路富集分析進一步分析差異表達基因的功能。探討了m6A調控因子表達及其在調控PASMC生物學行為中的作用，以揭示潛在的分子機制。
本研究利用MeRIP-seq和RNA-seq研究對照組和野百合堿（monocrotaline）誘導的PAH大鼠肺動脈組織中m6A和RNA表達差異。使用GO和KEGG進行功能富集分析，研究了差異表達功能。為了篩選候選m6A相關基因，使用STRING和Metascape數據庫構建了蛋白質-蛋白質互作（protein-protein interaction，PPI）網絡，并通過對候選相關基因表達進行實時PCR驗證。使用免疫組化染色、免疫熒光和Western blot技術進一步研究了m6A調控因子的表達水平，并對大鼠肺動脈平滑肌細胞（PASMC）進行增殖實驗。研究共鑒定出42個差異表達基因，這些基因的m6A甲基化水平表現出增加或減少（高甲基化或低甲基化）。主要與細胞外基質結構、MAPK和PI3K/AKT通路相關。在PAH組中檢測到候選基因著絲粒蛋白F（CENPF）表達增加。研究首次鑒定出一種富含亮氨酸的五肽重復序列（LRPPRC）的m6A reader，該reader在PAH大鼠模型中表達下調。體外實驗中，siRNA介導的Lrpprc下調導致大鼠PASMC增殖增強和Cenpf mRNA表達升高。本研究結果揭示了PAH大鼠肺動脈中轉錄組范圍的m6A修飾圖譜變化和相關調控機制，可能為未來治療策略提供新的靶點。
研究思路：
研究結果
圖1： PAH大鼠模型的m6A基本特征。組織來自肺動脈，每組n=3。
B. 餅圖顯示了轉錄本四個不重疊區域（（5' UTR、CDS、3' UTR和其他區域）的m6A peaks百分比。
C. 密度曲線顯示m6A peaks在轉錄本中的分布。轉錄本分為三部分，即5'UTR、CDS和3'UTR。
D. 兩組不同數量m6A peaks的基因比例。
E. 兩組RRACH motif情況。
m6A：N6-甲基腺苷；PAH：肺動脈高壓；MCT：野百合堿；UTR：非翻譯區；CDS：編碼序列。
（2）PAH中的差異m6A分析
圖2：PAH中的差異m6A分析
A. 相對于對照組，MCT組中差異表達的m6A peaks火山圖（| log2 FC |>0.585，p<0.05）。B. 差異m6A peaks中富集的前五個motif。
C. 具有差異m6A的上調m6A標記的轉錄本GO分析。
D. 具有差異m6A的下調m6A標記的轉錄本GO分析。
E. 具有差異m6A的上調m6A標記的轉錄本KEGG分析。
F. 具有差異m6A的下調m6A標記的轉錄本KEGG分析。
FC：倍數變化；GO：基因本體論；KEGG：京都基因與基因組百科全書。
表1：差異甲基化m6A peaks和相關基因數
（3）m6A修飾基因的轉錄譜
圖3：m6A修飾基因的轉錄譜
A. 差異表達基因的火山圖。B. 具有顯著變化的peak差異表達基因的四象限圖。
C. 總體m6A與mRNA表達水平呈正相關（r=0.1517793，p＜0.05）。
D-E. GO和KEGG分析分別具有顯著變化peak的差異表達基因。
F-G. 兩個數據集的DEG分布的火山圖。紅點代表上調DEG，藍點代表下調DEG，黑點代表非DEG。
H. 在兩個數據集中同時上調的DEG。
I. 在兩個數據集中同時下調的DEG。DEG：差異表達基因。
表2：數據庫中基因的mRNA表達水平（Log2（FC）
（4）PPI網絡建立和候選基因鑒定
圖4：蛋白質-蛋白質互作（PPI）網絡建立和候選基因鑒定
A–C. 基于上述基因構建的中樞網絡。D. MCT組和對照組（每組n=6）肺組織中Cenpf的相對mRNA表達，數據代表平均值±SEM，并使用Student t檢驗比較兩組組間差異。
E. 血小板衍生生長因子BB（PDGF-BB）刺激（30 ng/mL）下，大鼠PASMC中Cenpf mRNA相對表達水平（每組n=3），以α-Tublin作為內部對照歸一化。
F. 人類PASMCs在PDGF-BB刺激下的CENPF mRNA相對表達水平（每組n=3），以GAPDH作為內部對照歸一化。
G. 用siRNA轉染人PASMC 48h后，經PDGF-BB（30 ng/mL）再處理24h的CENPF mRNA相對表達水平（每組n=3）。
H-I. 以α-Tublin作為內部對照的PCNA蛋白水平的代表性Western blot和定量分析（每組n=3）。
J. 通過IGV觀察到Cenpf mRNA轉錄本上的m6A水平。
PPI：蛋白質-蛋白質互作；SEM：平均值的標準誤差；PASMCs：肺動脈平滑肌細胞；PDGF-BB：血小板衍生生長因子BB。
（5）肺動脈高壓（PAH）的肺動脈中，leucine-rich pentatricopeptide repeat-containing (LRPPRC) 蛋白表達水平下降
圖5： PAH患者的肺動脈LRPPRC表達降低。
B. 肺組織中LRPPRC的免疫組化染色（每組n=6），條形圖=100µm（左圖），放大的肺動脈（右圖）。
C. 肺組織中LRPPRC的免疫熒光染色（每組n=6），血管染色圖像（綠色）、目標蛋白染色圖像（紅色）、細胞核染色圖像（藍色），合并圖像（橙色），條形圖=50µm。
D-E. MCT和對照組肺組織中LRPPRC和α-Tubulin的代表性Western blot和定量分析（每組n=6）；數據表示為均值±標準誤（SEM），使用學生t檢驗比較兩組。
F-G. 大鼠PASMCs在PDGF-BB刺激下的LRPPRC水平的代表性Western blot和定量分析（30 ng/mL），以α-Tubulin作為內部對照（每組n=3）。
H. 在PDGF-BB刺激下（30 ng/mL），人類PASMCs中LRPPRC相對mRNA表達水平，以GAPDH作為內部對照（每組n=3）。
表3：m6A調控因子的mRNA表達水平
（6）LRPPRC下調促進PASMC增殖和Cenpf表達
圖6：LRPPRC下調促進肺動脈平滑肌細胞（PASMCs）的增殖和Cenpf表達。
A.在PASMCs中，用siRNA轉染PASMC 48小時，然后用PDGF-BB（30 ng/mL）再處理24小時后，Lrpprc水平的相對mRNA表達（每組n=3）。B-C. 代表性的Western blot和PCNA水平的定量分析，以α-Tubulin作為內部對照歸一化（每組n=3）。
D. EdU摻入實驗評估si-Lrpprc處理的PASMCs在PDGF-BB刺激下的細胞增殖能力。條形圖=100µm。目標細胞為紅色。
E.計算EdU染色的細胞比率。
F. Cenpf水平的相對mRNA表達（每組n=3）。
參考文獻：
Feng Y, Yu Z, Tang M, Li J, Peng B, Juaiti M, Tang Y, Liang B, Ouyang M, Liu Q, Song J. Transcriptome-Wide N6-Methyladenosine Alternations in Pulmonary Arteries of Monocrotaline-Induced Pulmonary Arterial Hypertension in Rats and Novel Therapeutic Targets. Biomedicines. 2024 Feb 4;12(2) pii: biomedicines12020364. doi: 10.3390/biomedicines12020364. PubMed PMID: 38397966.
標簽：
RNA甲基化
Copyright(C) 1998-2025 生物器材網 電話：021-64166852;13621656896 E-mail：info@bio-equip.com