高脂饮食诱导的小鼠NAFLD模型肝组织中m6A甲基化修饰表达谱分析
1.
2.
Analysis of m6A methylation expression profiles in liver tissue of high-fat diet-induced mouse models of NAFLD
1.
2.
通讯作者: 马万山,电子信箱:mwsqianyi@163.com。
编委: 吴洋
收稿日期: 2023-04-14 接受日期: 2023-09-12 网络出版日期: 2023-10-28
基金资助: |
|
Corresponding authors: MA Wanshan, E-mail:mwsqianyi@163.com.
Received: 2023-04-14 Accepted: 2023-09-12 Online: 2023-10-28
目的·利用微阵列芯片技术检测高脂饮食诱导的小鼠非酒精性脂肪性肝病(non-alcoholic fatty liver disease,NAFLD)模型中肝组织mRNA的m6A甲基化修饰和基因表达的变化。方法·以6~8周龄雄性C57BL/6J小鼠为实验动物,高脂饲料喂养16周诱导建立NAFLD模型(高脂组,n=10);另设基础组(n=10)为对照,给予含10%脂肪的基础饲料喂养。苏木精-伊红(hematoxylin-eosin,H-E)染色评估小鼠肝组织病理改变,判断NAFLD模型是否构建成功。运用甲基化RNA免疫共沉淀技术(methylated RNA immunoprecipitation,MeRIP)和微阵列测序技术(microarray表达谱分析)检测2组小鼠肝组织中mRNA的m6A甲基化和表达水平变化。结果·基础组小鼠肝脏呈鲜红色,少见脂肪沉积;高脂组小鼠肝脏边界黄润,H-E染色可见肝细胞中脂滴弥漫浸润且相互融合,提示高脂饲料诱导的NAFLD模型构建成功。MeRIP-微阵列芯片检测结果显示,与基础组相比,高脂组小鼠肝脏中共有320个基因m6A甲基化修饰水平变化显著(P<0.05且变化倍数>1.5),其中有108个上调基因和212个下调基因。将组间m6A甲基化水平差异显著的基因与mRNA差异表达基因取交集,发现有163个基因m6A甲基化水平和mRNA表达水平均差异显著。结论·高脂饮食诱导的小鼠NAFLD模型中肝组织mRNA的m6A修饰变化显著,且该变化与mRNA的基因表达有关。
关键词:
Objective ·To detect the differences in m6A methylation modification and gene expression of liver tissue mRNA in high-fat diet-induced mouse models of non-alcoholic fatty liver disease (NAFLD) using microarray technology. Methods ·The NAFLD models were established in 6-8 weeks old male C57BL/6J mice fed with high-fat chow for 16 weeks (high-fat group, n=10). The basal group (n=10) was given 10% fat diet. Hematoxylin-eosin (H-E) staining was used to assess the histopathological changes in liver tissue and to determine the success of the NAFLD models. The changes of mRNA m6A methylation and expression levels in the liver tissues of the two groups were detected by using methylated RNA immunoprecipitation (MeRIP) and microarray expression profiling. Results ·The livers of the mice in the basal group were bright red with few fat deposits, while the livers of the mice in the high-fat group were yellowish with diffuse infiltration and fusion of lipid droplets in the hepatocytes by H-E staining, suggesting that the high-fat diet-induced NAFLD models were successfully constructed. The results of the MeRIP-microarray showed that the m6A methylation levels of 320 genes in the livers of mice in the high-fat group were significantly altered compared with those in the basal group (P<0.05 and fold change>1.5), of which 108 genes were up-regulated and 212 genes were down-regulated. Genes with significant differences in m6A methylation levels between the two groups were intersected with those with differentially expressed mRNAs, and 163 genes were found to have significant differences in both m6A methylation level and mRNA expression level. Conclusion ·The change in m6A modification of liver tissue mRNA in the high-fat diet-induced mouse models of NAFLD is significant and the change is associated with the gene expression of mRNA.
Keywords:
本文引用格式
刘君君, 逯素梅, 张炳杨, 李永清, 马万山.
LIU Junjun, LU Sumei, ZHANG Bingyang, LI Yongqing, MA Wanshan.
随着社会经济的进步和人们物质生活的改善,非酒精性脂肪性肝病(non-alcoholic fatty liver disease,NAFLD)的发病率激增,已成为我国最常见的慢性肝病。有数据显示我国成年人NAFLD发病率为12%~15%[1],预计到2030年我国NAFLD患病人数将达到3.14亿[2]。NAFLD主要表现为肝细胞脂质代谢异常,肝细胞受损后引起炎症、胰岛素抵抗、纤维化等反应[3]。既往研究[4-7]表明,肥胖、2型糖尿病、高脂血等代谢综合征与NAFLD的发生和发展关系密切,两者之间互为因果,并行发展[8-9]。NAFLD可引起失代偿期肝硬化、肝细胞癌,还参与2型糖尿病和动脉粥样硬化的发病[10-12]。NAFLD使患者生活质量降低,寿命缩短,还加重了国家的卫生经济负担[13]。
本研究借助高脂饲料喂养构建小鼠NAFLD模型,利用甲基化RNA免疫共沉淀技术(methylated RNA immunoprecipitation,MeRIP)和微阵列测序技术,检测差异m6A甲基化修饰和差异基因表达,对数据进行联合分析,为NAFLD发病的基因调控机制研究提供理论依据。
1 对象与方法
1.1 实验动物
SPF级C57B/6J雄性小鼠,6~8周龄,体质量18~20 g,购自北京维通利华实验动物技术有限公司。动物生产许可证号为SCXK(京)2016-0006,使用许可证号为SYXK(鲁)2018-0009。饲养环境采用12 h明和12 h暗交替,温度为20~26 ℃、日温差1 ℃,相对湿度达40%~70%,新鲜空气换气次数15次/h,气流速度<0.18 m/s,压差25 Pa,氨浓度≤15 mg/m3,噪音<55 dB,照度150~300 Lux。
1.2 仪器与试剂
高速冷冻离心机(美国Thermo Fisher Scientific),摇床(海门市其林贝尔仪器制造有限公司),FluorChem-E凝胶成像分析平台(美国Cell Biosciences),生物安全柜(中国济南鑫贝西生物技术有限公司),Nano drop 2000分光光度计(美国Thermo Fisher Scientific),微量移液器(德国Eppendorf),TRIzol Reagent(美国Invitrogen),总RNA提取试剂(中国上海奕杉生物科技有限公司),异丙醇(上海国药试剂有限公司),三氯甲烷(上海国药试剂有限公司),去离子水(北京天根生化科技有限公司),Western一抗稀释液(上海碧云天生物技术有限公司),m6A抗体(德国Synaptic Systems),辣根过氧化物酶HRP标记山羊抗兔二抗(北京中杉金桥生物技术有限公司)。
1.3 方法
1.3.1 小鼠NAFLD模型建立及分组
将20只小鼠随机分为基础组和高脂组,每组10只,分别给予10%脂肪的基础饲料(美国Research Diets)和60%脂肪的高脂饲料(美国Research Diets)喂养。喂养期间每3 d更换1次水和饲料,每周称量小鼠体质量,监测精神状态等。持续喂养小鼠16周。
1.3.2 肝组织病理观察
小鼠造模成功后,将小鼠深度麻醉,取其肝脏组织。先观察肝脏的大小、颜色、质地等,再将生理盐水清洗干净的肝脏用10%的甲醛溶液固定,常规处理后制备石蜡包埋切片。将石蜡切片脱蜡水化后进行苏木精-伊红(hematoxylin-eosin,H-E)染色,光学显微镜下观察肝脏脂肪变性情况。
1.3.3 RNA的提取与质量检测
用分析天平称重小鼠肝脏组织100 mg,加入1 mL TRIzol试剂研磨组织后进行总RNA的提取,并根据试剂说明书用DNase处理以去除潜在的DNA污染物。样品经主成分分析(principal component analysis,PCA)合格后,使用Nano drop 2000分光光度计检测总RNA的浓度与质量,RNA均为单峰且D(260 nm)/D(280 nm)的比值在1.8~2.0视为质量合格。将RNA片段化处理后均分为2个部分,一部分使用能特异识别m6A修饰的抗体进行免疫共沉淀(immunoprecipitation,IP)作为IP样品,另一部分不进行IP留作Input样品。设计有m6A修饰的阳性对照引物和非m6A修饰的阴性对照引物,对上述样品进行反转录和qPCR。若阳性对照的 %(MeRIP/Input)值是阴性对照2倍以上,即认为IP富集实验成功,RNA甲基化免疫共沉淀(methylated RNA immunoprecipitation,MeRIP)质检合格。将质量合格的样品送至上海康成生物科技有限公司进行mRNA表观转录组学芯片表达谱分析。同时将提取好的RNA稀释成3个浓度进行斑点杂交(Dot-blot)。使用m6A抗体(1∶1 000,1 ng/μL)孵育过夜,再加入HRP二抗(1∶1 000)孵育1 h,加入发光液后放置于凝胶成像分析仪中,观察基础组和高脂组小鼠肝组织中RNA的m6A表达情况。
1.4 统计学方法
通过表达差异倍数(Fold change>1.5)和统计学显著性(P<0.05),筛选2组之间差异m6A甲基化基因和差异表达基因。
2 结果
2.1 小鼠NAFLD肝组织中m6A修饰水平变化
高脂饲料喂养16周后,基础组和高脂组的小鼠体质量出现明显差异,高脂组小鼠平均体质量为(45±5)g,而基础组小鼠为(25±5)g。高脂组的小鼠肝脏边界黄润,表面有颗粒感,且有明显的弥漫性肿大(图1A)。H-E染色见高脂组小鼠肝组织细胞质中
图1
图1
高脂饮食诱导的小鼠NAFLD肝组织中m6A修饰水平变化
Note: A. General observation of the livers. B.H-E staining showed the cytoplasm of liver tissue in the high-fat group was filled with lipid droplets. C.Dot-blot showed the increase of m6A methylation levels in the high-fat group.
Fig 1
Changes of m6A modification levels in high-fat diet mouse livers
2.2 基因微阵列芯片实验流程的质量控制
图2
图2
样品质控和MeRIP-qPCR质量检测
Note: A. PCA-plot. B.MeRIP-qPCR quality inspection bar chart showed negative and positive control.
Fig 2
Sample quality control and MeRIP-qPCR quality testing
2.3 小鼠肝脏mRNA差异m6A修饰基因与差异表达基因的联合分析
2组小鼠肝组织基因芯片差异甲基化分析发现,与基础组相比,高脂组有320个mRNA发生m6A甲基化差异修饰(differentially methylated m6A site,DMMS)(变化倍数>1.5,P<0.05),其中有108个上调基因、212个下调基因(图3A)。将2组间差异表达的基因进行排序筛选,将m6A上调差异的前10位基因和m6A下调差异的前10位基因作热图分析,见图3B。上调和下调前10位基因的具体信息见表1。进一步将含有m6A甲基化水平与mRNA表达水平差异均1.5倍及以上的基因取交集,共得到163个差异基因,包括25个m6A甲基化水平升高且mRNA表达上调的基因,46个m6A甲基化水平升高且mRNA表达下调的基因,58个m6A甲基化水平降低且mRNA表达上调的基因,34个m6A甲基化水平降低且mRNA表达下调的基因(图3C)。
图3
图3
小鼠肝脏mRNA差异m6A修饰基因与差异表达基因的联合分析
Note: A. Volcano plot of differential m6A methylation mRNAs. B.Cluster diagram showed top 10 differentially m6A-methylated mRNAs. C.Four quadrant diagrams showed differential methylation associated with gene expression in the high-fat group and the control group.
Fig 3
Conjoint analysis of differential m6A methylation mRNAs and differentially expressed genes in the livers of mice
表1 前10位m6A差异修饰的mRNA信息
Tab 1
Gene | mRNA_ID | Locus | RNA length/nt | Fold change | P value |
---|---|---|---|---|---|
Up-regulated gene | |||||
Micalc1 | ENSMUST00000033033 | chr7: 112374345-112395355: + | 2 259 | 2.840 | 0.004 |
Cyp2b10 | ENSMUST00000005477 | chr7: 25897676-25926559: + | 1 822 | 2.721 | 0.011 |
Ksr2 | ENSMUST00000180430 | chr5: 117414000-117775003: + | 13 125 | 2.484 | 0.001 |
Ybx1 | ENSMUST00000079644 | chr4: 119277981-119294604: - | 1 628 | 2.359 | 0.023 |
Paqr7 | ENSMUST00000095074 | chr4: 134497004-134510235: + | 3 550 | 2.274 | 0.002 |
Casc5 | ENSMUST00000099542 | chr2: 119047119-119104121: + | 6 525 | 2.270 | 0.000 |
Ubr3 | ENSMUST00000131553 | chr2: 69897303-69938659: + | 1 801 | 2.113 | 0.005 |
Dpm3 | ENSMUST00000107462 | chr3: 89259358-89267077: + | 827 | 1.978 | 0.025 |
Cd36 | ENSMUST00000082367 | chr5: 17782016-17835696: - | 3 016 | 1.965 | 0.021 |
Tsc22d1 | ENSMUST00000142683 | chr14: 76487759-76506998: + | 797 | 1.928 | 0.028 |
Down-regulated gene | |||||
Esd | ENSMUST00000176957 | chr14: 74732384-74749848: + | 1 117 | 0.356 | 0.010 |
Slc22a26 | ENSMUST00000120522 | chr19: 7781041-7802578: - | 2 854 | 0.421 | 0.010 |
Folr1 | NM_001252553 | chr7: 101858331-101870788: - | 1 337 | 0.437 | 0.000 |
Slc46a3 | ENSMUST00000138244 | chr5: 147893881-147894815: - | 388 | 0.437 | 0.001 |
Inpp5d | ENSMUST00000167032 | chr1: 87676239-87720502: + | 4 125 | 0.462 | 0.010 |
Hoxd9 | ENSMUST00000059272 | chr2: 74697727-74700208: + | 2 136 | 0.471 | 0.012 |
Romo1 | ENSMUST00000088610 | chr2: 156144039-156145797: + | 569 | 0.484 | 0.005 |
Litaf | ENSMUST00000117360 | chr16: 10960824-10975579: - | 763 | 0.492 | 0.014 |
Mdh2 | ENSMUST00000019323 | chr5: 135778480-135790391: + | 1 456 | 0.509 | 0.032 |
Sycp3 | ENSMUST00000125612 | chr10: 88459664-88473236: + | 1 169 | 0.513 | 0.044 |
2.4 小鼠肝组织差异m6A修饰与差异表达的基因信息
进一步对小鼠肝组织m6A修饰和表达差异均1.5倍及以上的基因与m6A修饰差异显著的基因(变化倍数>1.5,且P<0.05)取交集。将高甲基化高表达的基因和m6A甲基化修饰显著上调的基因取交集后,筛选出19个基因(图4A);将高甲基化低表达的基因和m6A
图4
图4
小鼠肝组织差异显著的m6A修饰且差异表达的基因信息
Fig 4
Differential m6A modification and differentially expressed genes in mouse liver tissue
3 讨论
NAFLD的病因及发病机制比较复杂,目前仍未完全阐明。现阶段NAFLD的治疗主要通过改善生活方式以稳定机体调节机制,积极寻找并去除可能的病因和诱因,以及减肥等,尚未发现可以根治该病的药物[16]。因此,NAFLD的治疗也是全球面临的一大难题。
近年来,微阵列芯片作为一种快速、并行、高效的分子生物学研究工具,对探究基因表达的调控机制有重要的辅助作用。本研究通过高脂饲料喂养小鼠,构建了稳定可靠的小鼠NAFLD模型。研究发现NAFLD小鼠总m6A甲基化水平升高,在RNA样本质检合格的前提下,进一步进行真核生物转录组基因m6A微阵列芯片表达谱分析,基于检测结果建立了高脂饮食诱导下肝组织m6A修饰差异表达谱。其中,有320个基因m6A甲基化修饰水平变化显著(P<0.05且变化倍数>1.5),且发现有163个基因m6A甲基化水平和mRNA表达水平均差异显著。该结果提示m6A甲基化水平可能影响转录组基因表达的变化,且与肝脏脂质代谢的调控密切相关,这163个差异基因在NAFLD发生和发展中可能发挥重要的调控作用。
现有文献证实,筛选出的部分差异表达基因与NAFLD的发生和发展有不同程度的相关性。例如,有研究[19]表明小鼠肝脏组织Hacd2缺失可通过减少血清中长链脂肪酸的生成增加以产热为主的能量消耗,改善高脂饮食诱导的NAFLD等。Cd36缺失的小鼠肝脏中Srebp1、Pparγ和Fasn等成脂相关基因表达降低,在HepG2细胞中也证实CD36通过与INSIG2结合,促进肝脏SREBP1表达和脂质从头合成途径,从而促进NAFLD的发生和发展[20]。外周血PPARγ和CIDEA的mRNA分子表达水平与血浆三酰甘油和总胆固醇水平呈显著负相关[21]。SANS等[22]的研究表明:在肥胖患者中,肝脏CIDEA的表达与NAFLD活性评分和肝损伤密切相关;肝脏CIDEA的表达与肥胖程度呈正相关,且会加重NAFLD的进展。TANG等[23]报道,高脂饮食小鼠肝脏中CIDEA的m6A甲基化水平升高与mRNA表达水平呈显著正相关,经实验证实mettl16介导的m6A甲基化通过上调CIDEA表达水平促进NAFLD进展。有研究[24]显示,Ces1/Es-x在小鼠肝脏脂肪代谢中起调节作用,Ces1/Es-x基因敲除小鼠出现肥胖、脂肪肝、高胰岛素血症等NAFLD表现。高脂诱导的肥胖小鼠中,MECR基因的表达水平与2型糖尿病密切相关,ATP降低MECR蛋白水平可能参与胰岛素抵抗的发生[25]。高脂饮食可诱导小鼠肝脏发生严重的脂肪样变,而纤维化相关基因Pdgfrb、Tgfb在小鼠体内表达升高[26]。以上研究都提示,本次研究中筛选出的差异基因对NAFLD的发病机制研究有重要意义。
本研究未对筛选出的差异基因进行进一步验证。后续我们将围绕差异基因进行验证,并通过干扰或过表达差异基因,开展进一步的功能分析,为NAFLD的预防和治疗提供新的靶点。
综上所述,本研究通过比较正常小鼠与NAFLD模型小鼠中肝脏组织的m6A甲基化及mRNA表达差异,发现NAFLD小鼠肝组织中m6A修饰水平发生了显著变化,且mRNA表达水平变化显著,表明m6A甲基化可能通过影响mRNA的表达调控NAFLD的发生和发展。本研究建立了主要差异基因表达谱,为筛选更多的NAFLD治疗靶点提供了一定的理论基础。
作者贡献声明
马万山、逯素梅、刘君君参与了实验设计以及论文的写作和修改,刘君君、张炳杨、李永清负责实验操作与数据分析。所有作者均阅读并同意了最终稿件的提交。
The study was designed by MA Wanshan, LU Sumei and LIU Junjun. The manuscript was drafted and revised by MA Wanshan, LU Sumei and LIU Junjun. The experiments were performed by LIU Junjun, ZHANG Bingyang and LI Yongqing. The data were analyzed by LIU Junjun, ZHANG Bingyang and LI Yongqing. All the authors have read the last version of paper and consented for submission.
利益冲突声明
所有作者声明不存在利益冲突。
All authors disclose no relevant conflict of interests.
参考文献
/
〈 |
|
〉 |
![](../../../richhtml/1674-8115/richHtml_jats1_1/images/1.png)