癌变·畸变·突变 ›› 2021, Vol. 33 ›› Issue (6): 405-409.doi: 10.3969/j.issn.1004-616x.2021.06.001

• 论著 •    下一篇

甲基丙烯酸环氧丙酯诱导16HBE细胞恶性转化相关m6A修饰异常mRNA的分析

王苗1,2, 王全凯1,2, 李昕苇1,2, 马顺鹏1,2, 乌瀚宝栎尔1, 许建宁1,2   

  1. 1. 中国疾病预防控制中心职业卫生与中毒控制所, 北京 100050;
    2. 中国疾病预防控制中心化学污染与健康安全重点实验室, 北京 100050
  • 收稿日期:2021-08-26 修回日期:2021-11-03 出版日期:2021-11-30 发布日期:2021-12-04
  • 通讯作者: 许建宁,E-mail:jnx999@263.net E-mail:jnx999@263.net
  • 作者简介:王苗,E-mail:2091419301@qq.com。
  • 基金资助:
    国家自然科学基金(81673221)

m6A methylation of mRNAs in glycidyl methacrylate-induced malignant transformation of 16HBE cells

WANG Miao1,2, WANG Quankai1,2, LI Xinwei1,2, MA Shunpeng1,2, WUHAN Baolier1, XU Jianning1,2   

  1. 1. National Institute of Occupational Health and Poison Control, Chinese Center for Disease Control and Prevention, Beijing 100050;
    2. Key Laboratory of Chemical Safety and Health, Chinese Center for Disease Control and Prevention, Beijing 100050, China
  • Received:2021-08-26 Revised:2021-11-03 Online:2021-11-30 Published:2021-12-04

摘要: 目的: 检测甲基丙烯酸环氧丙酯(GMA)诱导16HBE细胞恶性转化相关mRNA的m6A甲基化修饰水平,筛选出m6A修饰的mRNA并进行功能注释。方法: GMA(8 μg/mL)重复染毒16HBE细胞后,收集GMA染毒组和DMSO对照组的第30代细胞,采用高通量人类表观转录组芯片检测16HBE细胞恶性转化相关mRNA的m6A修饰水平,应用Visual Studio Code软件进行m6A修饰mRNA的筛选,并利用Omicshare工具对这些mRNA进行GO富集分析和KEGG通路分析。结果: 与对照组细胞比较,GMA诱导的恶性转化16HBE细胞中m6A修饰水平异常的mRNA共有454个,其中高甲基化水平mRNA 334个,低甲基化水平mRNA 120个;表达水平异常的mRNA共有434个,其中高表达mRNA 236个,低表达mRNA 198个。m6A修饰的mRNA有45个,GO富集分析结果显示上述mRNA主要富集于SNAP受体活性、SNARE结合、SNARE复合体等生物学过程,KEGG通路分析主要富集于囊泡运输中的SNARE相互作用、非同源末端连接、苯丙氨酸代谢等通路。结论: m6A修饰mRNA可能在GMA诱导16HBE恶性转化过程发挥重要作用。

关键词: N6-甲基腺嘌呤, 甲基丙烯酸环氧丙酯, 高通量基因芯片, 恶性转化细胞

Abstract: OBJECTIVE: To investigate m6A methylation levels of mRNAs which are related to glycidyl methacrylate (GMA)-induced malignant transformation of 16HBE cells and their functional roles. METHODS: After 16HBE cells were repeatedly exposed to GMA(8 μg/mL),the 30th generation cells of GMA group and DMSO control group were collected. m6A methylation levels of mRNAs which were related to 16HBE malignant transformation were detected by high-throughput human apparent transcriptome chip. The differentially m6A-methylated mRNAs and differentially expressed mRNAs were screened by Visual Studio Code. The Omicshare tools were used for GO enrichment and KEGG pathway analyses of these mRNAs. RESULTS: A total of 454 lncRNAs with differential m6A methylation were identified in GMA-induced 16HBE malignant transformed cells,including 334 which were hypermethylated mRNAs and 120 which were hypomethylated mRNAs (|FC|>2.0,P<0.05). There was a total of 434 differentially expressed mRNAs,among which 236 were up-regulated and 198 were down-regulated (|FC|>2.0,P<0.05). There were 45 m6A-methylated mRNAs. The GO analyses showed that SNAP receptor activity,SNARE binding and SNARE complex were the main biological processes. KEGG pathway analyses showed that m6A-methylated mRNAs were enriched in SNARE interactions in vesicular transport,non-homologous end-joining,and phenylalanine metabolism. CONCLUSION: Our data indicate that m6A methylation played an important role in the process of GMA-induced 16HBE malignant transformation.

Key words: N6-methyladenosine, glycidyl methacrylate, high throughput microarray, malignant transformation

中图分类号: