肾阳虚体质大鼠肝细胞全基因组DNA甲基化的测序与分析

欧阳双; 李晗; 蔡海洋; 吴晶; 汤朝晖; 李炜弘

doi:10.13288/j.11-2166/r.2022.09.015

您当前的位置：

首页 >

文章列表页 >

肾阳虚体质大鼠肝细胞全基因组DNA甲基化的测序与分析

实验研究 | 更新时间：2022-04-25

- 肾阳虚体质大鼠肝细胞全基因组DNA甲基化的测序与分析
- Whole-genome DNA Methylation Sequencing and Analysis of Hepatocytes in Rats with Kidney Yang Deficiency Constitution
- 中医杂志 2022年63卷第9期页码：875-881
- 作者机构：
  
  成都中医药大学，成都市温江区柳台大道1166号，611137
- 作者简介：
  
  lwh@cdutcm. edu. cn
- 基金信息：
  
  国家自然科学基金(81873204)
- DOI：10.13288/j.11-2166/r.2022.09.015
  中图分类号：
- 收稿：2021-11-09，
  
  修回：2022-02-25，
  
  纸质出版：2022-05-02
- 稿件说明：
移动端阅览
欧阳双,李晗,蔡海洋等.肾阳虚体质大鼠肝细胞全基因组DNA甲基化的测序与分析[J].中医杂志,2022,63(09):875-881.

OUYANG Shuang,LI Han,CAI Haiyang,et al.Whole-genome DNA Methylation Sequencing and Analysis of Hepatocytes in Rats with Kidney Yang Deficiency Constitution[J].Journal of Traditional Chinese Medicine,2022,63(09):875-881.
欧阳双,李晗,蔡海洋等.肾阳虚体质大鼠肝细胞全基因组DNA甲基化的测序与分析[J].中医杂志,2022,63(09):875-881. DOI： 10.13288/j.11-2166/r.2022.09.015.

OUYANG Shuang,LI Han,CAI Haiyang,et al.Whole-genome DNA Methylation Sequencing and Analysis of Hepatocytes in Rats with Kidney Yang Deficiency Constitution[J].Journal of Traditional Chinese Medicine,2022,63(09):875-881. DOI： 10.13288/j.11-2166/r.2022.09.015.

摘要

目的

探讨肾阳虚体质大鼠肝细胞DNA甲基化特征。

方法

选取本课题组前期研发的自发性肾阳虚体质模型大鼠子一代、子四代雄性各3只分别组成子一代组和子四代组，3只正常SD大鼠设为空白组。取各组大鼠肝脏组织，进行DNA抽取与质检，构建DNA文库，进行甲基化测序，与参考基因组进行唯一比对，定位富集区域。将子一代组、子四代组富集区域分别与空白组比较，取差异富集区域。子一代组、子四代组差异富集区域取交集，将富集区域关联到基因，并根据组间富集区域富集样本个数差异筛选组间差异的甲基化富集区域。针对差异关联基因集进行GO、KEGG功能分析。

结果

子一代组共获得276 091个富集区域，子四代组共获得309 430个富集区域，空白组共获得275 413个富集区域。与空白组大鼠相比，子一代组肾阳虚体质大鼠差异富集区域共372个，其中甲基化水平增高89个，甲基化水平降低283个；子四代组肾阳虚体质大鼠差异富集区域共51条，其中甲基化水平增高30个，甲基化水平降低21个。取子一代组与子四代组差异富集区域的交集，初步筛出差异富集区域12个。对差异富集区域进行GO、KEGG功能分析，密切关联到基因11个，其中甲基化水平增高有Gdap1、Cdkal1、Cd300lb，甲基化水平降低有Zmiz1、Fam111a、Sh3bp4、As3mt、Col16a1、LOC103690271、Wdpcp、LOC102553861，基因功能涉及线粒体膜生成、代谢过程、细胞增殖、细胞吞噬、信号转导、生物调节、蛋白定位等。

结论

肾阳虚体质大鼠存在肝细胞DNA甲基化的差异表达，差异甲基化位点其关联基因集中于肿瘤发生和线粒体能量代谢、糖脂代谢等生物学功能紊乱。

Abstract

Objective

To explore the characteristics of DNA methylation of hepatocytes in rats with kidney

yang

deficiency constitution.

Methods

Three first-generation and three fourth-generation male offspring of spontaneous kidney

yang

deficiency model rats developed by our previous research were selected as the first-generation offspring group and fourth-generation offspring group， respectively， and three normal SD rats were set as the blank group. The liver tissues of rats in each group were collected for DNA extraction and quality inspection， and DNA libraries were then constructed. Methylation sequencing was performed， and the results were uniquely compared with the reference genome to identify the location of enrichment regions. The enriched region of the first-generation offspring group and the fourth-generation offspring group were compared to that of the blank group， respectively， to identify the differentially enriched regions. The differentially enriched regions of first-generation offspring group and the fourth-generation offspring group were intersected， and the enrichment region was associated with genes. The methylation-enriched regions that differ between groups were then identified according to the differences in the number of enriched samples in enriched regions between groups. GO and KEGG functional analysis was performed on differentially expressed genes.

Results

A total of 276 091 enrichment regions were obtained in the first-generation offspring group， while 309 430 enrichment regions and 275 413 enrichment regions were identified in the fourth-generation offspring group and the blank group， respectively. Compared to the rats in the blank group， there were 372 differentially enriched regions in the kidney

yang

deficiency rats in the first-generation offspring group， among which 89 regions had increased methylation levels and 283 had decreased methylation levels； there were 51 differentially enriched regions in the kidney

yang

deficiency rats in the fourth-generation offspring group， among which 30 regions had increased methylation levels and 21 had decreased methylation levels. The intersection of the differentially enriched areas of the first-generation offspring group and the fourth-generation offspring group preliminarily identified 12 differentially enriched regions. GO and KEGG functional analysis of differentially enriched regions showed that 11 genes were closely related. Gdap1， Cdkal1， and Cd300lb had an increase in the methylation level， while Zmiz1， Fam111a， Sh3bp4， As3mt， Col16a1， LOC103690271， Wdpcp， and LOC102553861 had a decrease in the methylation level. Gene function was involved with mitochondrial membrane formation， metabolic process， cell proliferation， cell phagocytosis， signal transduction， biological regulation， protein localization， etc.

Cnclusion

DNA methylation of hepatocytes is differentially expressed in rats with kidney

yang

deficiency constitution， and the genes associated with differential methylation sites are mainly involved with tumorigenesis and biological function disorders such as mitochondrial energy metabolism and glucose and lipid metabolism.

关键词

Keywords

references

BIRD A ． Perceptions of epigenetics ［J］． Nature ， 2007 ， 447 （ 7143 ）： 396 - 398 ．

王米渠，丁维俊，李炜弘，等 . 基因表达谱芯片发现热药疗寒的相关七类基因初报［J］. 中医杂志， 2003 ， 44 （ 2 ）： 119 - 122 .

曾跃琴，李炜弘，张天娥，等 . 肾阳虚证免疫相关基因CPG岛调控机制研究［J］. 时珍国医国药， 2013 ， 24 （ 6 ）： 1515 - 1517 .

和靖 . 自发性肾阳虚体质动物模型的构建及评价研究［D］. 成都：成都中医药大学， 2020 .

陈博武，贾茹，钱晨，等 . 基于中医四诊方法对哮喘大鼠肾阳虚证模型建立的评价研究［J］. 中国中医药科技， 2017 ， 24 （ 5 ）： 537 - 539，543 .

卢德赵，沃兴德，沃立科，等 . 温补肾阳药对甲状腺切除的肾阳虚大鼠肝线粒体蛋白质组的影响［J］. 中国中药杂志， 2009 ， 34 （ 10 ）： 1251 - 1256 .

李炜弘，雍小嘉，范怀昌，等 . 老龄肾阳虚证的差异表达基因谱分析［J］. 时珍国医国药， 2009 ， 20 （ 5 ）： 1210 - 1212 .

李炜弘，范怀昌，丁维俊，等 . 运用基因芯片技术分析老龄肾阳虚证与Wnt/B-catenin信号通路的相关性［J］. 四川中医， 2009 ， 27 （ 12 ）： 15 - 17 .

李炜弘，许嗣立，张天娥，等 . 老龄肾阳虚证差异表达基因的筛选与鉴定［J］. 辽宁中医杂志， 2012 ， 39 （ 4 ）： 590 - 592 .

乔卫龙 . 基于线粒体母系遗传研究肾阳虚体质能量代谢低下的分子机制［D］. 成都：成都中医药大学， 2017 .

卢德赵，沃兴德，施孟如，等 . 激素型肾阳虚动物肝线粒体蛋白质组与能量代谢相关性［J］. 中国生物化学与分子生物学报， 2005 ， 21 （ 6 ）： 807 - 814 .

NIEMANN A ， RUEGG M ， PADULA VLA ， et al . Ganglioside-induced differentiation associated protein 1 is a regulator of the mitochondrial network： new implications for Charcot-Marie-Tooth disease ［J］. J Cell Biol ， 2005 ， 170 （ 7 ）： 1067 - 1078 .

PEDROLA L ， ESPERT A ， WU X ， et al . GDAP1， the protein causing Charcot-Marie-Tooth disease type 4A， is expressed in neurons and is associated with mitochondria ［J］. Hum Mol Genet ， 2005 ， 14 （ 8 ）： 1087 - 1094 .

LACKNER LL . Shaping the dynamic mitochondrial network ［J］. BMC Biol ， 2014 ， 12 ： 35 . doi： 10.1186/1741-7007-12-35 http://dx.doi.org/10.1186/1741-7007-12-35 .

SAFIULINA D ， KAASIK A . Energetic and dynamic： how mitochondria meet neuronal energy demands ［J］. PLoS Biol ， 2013 ， 11 （ 12 ）： e1001755 .

卢德赵，沃兴德，沃立科，等 . 甲状腺切除的肾阳虚大鼠肝线粒体蛋白质组的研究［J］. 中医药学报， 2008 ， 36 （ 4 ）： 11 - 15 .

OHARA-IMAIZUMI M ， YOSHIDA M ， AOYAGI K ， et al . Deletion of CDKAL1 affects mitochondrial ATP generation and first-phase insulin exocytosis ［J］. PLoS One ， 2010 ， 5 （ 12 ）： e15553 .

陈继承，庄耀东，李灿东，等 . 代谢综合征痰证人群不同症候群与CDKAL1、CDKN2A的相关性［J］. 中华中医药杂志， 2017 ， 32 （ 8 ）： 3426 - 3430 .

KAWAI VK ， LEVINSON RT ， ADEFURIN A ， et al . A genetic risk score that includes common type 2 diabetes risk variants is associated with gestational diabetes ［J］. Clin Endocrinol （Oxf）， 2017 ， 87 （ 2 ）： 149 - 155 .

TARNOWSKI M ， MALINOWSKI D ， PAWLAK K ， et al . GCK， GCKR， FADSI， DGKB/TMEM195 and CDKALL1 gene polymorphisms in women with gestational diabetes ［J］. Can J Diabetes ， 2017 ， 41 （ 4 ）： 372 - 379 .

佘茁萃，吕俊，蔡昂 . 癌症恶病质中医证治探讨［J］. 中国民间疗法， 2020 ， 28 （ 21 ）： 11 - 13 .

郭威堂，汪鸣，韩琦，等 . 论阳虚体质与肿瘤的发病关系［J］. 中医药信息， 2016 ， 33 （ 6 ）： 17 - 19 .

AKAMATSU S ， TAKATA R ， HAIMAN CA ， et al . Common variants at 11q12， 10q26 and 3p11.2 are associated with prostate cancer susceptibility in Japanese ［J］. Nat Genet ， 2012 ， 44 （ 4 ）： 426 - 429 .

FERNANDEZ-RETANA J ， ZAMUDIO-MEZA H ， RODRIGUEZ-MORALES M ， et al . Gene signature based on degradome-related genes can predict distal metastasis in cervical cancer patients ［J］. Tumour Biol ， 2017 ， 39 （ 6 ）： 1010428317711895 .

LI X ， THYSSEN G ， BELIAKOFF J ， et al . The novel PIAS-like protein hZimp10 enhances Smad transcriptional activity ［J］. J Biol Chem ， 2006 ， 281 （ 33 ）： 23748 - 23756 .

LEE J ， BELIAKOFF J ， SUN Z . The novel PIAS-like protein hZimp10 is a transcriptional co-activator of the p53 tumor suppressor ［J］. Nucleic Acids Res 2007 ， 35 （ 13 ）： 4523 - 4534 .

MATHIOS D ， HWANG T ， XIA Y ， et al . Genome-wide investigation of intragenic DNA methylation identifies ZMIZ1 gene as a prognostic marker in glioblastoma and multiple cancer types ［J］. Int J Cancer ， 2019 ， 145 （ 12 ）： 3425 - 3435 .

宋晓丽，耿志荣，王志林．无机砷甲基化研究进展［J］．无机化学学报， 2012 ， 28 （ 10 ）： 2025 - 2035 .

SMITH AH ， LOPIPERO PA ， BATES MN ， et al . Public health Arsenic epidemiology and drinking water standards ［J］. Science ， 2009 ， 296 （ 5576 ）： 2145 - 2146 .

SUN M ， TAN J ， WANG M ， et al . Inorganic arsenic-mediated upregulation of AS3MT promotes proliferation of nonsmall cell lung cancer cells by regulating cell cycle genes ［J］. Environ Toxicol ， 2021 ， 36 （ 2 ）： 204 - 212 .

RECIO-VEGA R ， HERNANDEZ-GONZALEZ S ， MICHEL-RAMIREZ G ， et al . Association between the polymorphism of three genes involved in the methylation and efflux of arsenic （As3MT， MRP1， and P-gp） with lung cancer in a Mexican cohort ［J］. J Appl Toxicol ， 2021 ， 41 （ 9 ）： 1357 - 1366 .

SIMONS M ， GLOY J ， GANNER A ， et al . Inversin， the gene product mutated in nephronophthisis type Ⅱ， functions as a molecular switch between Wnt signaling pathways ［J］. Nat Genet ， 2005 ， 37 （ 5 ）： 537 - 543 .

LI Z ， MAO S ， JIN J . Activation of hedgehog pathway in acute myeloid leukemia patients ［J］. Int J Clin Exp Pathol ， 2017 ， 10 （ 8 ）： 8605 - 8609 .

MERCHANT AA ， MATSUI W . Targeting Hedgehog：a cancer stem cell pathway ［J］. Clin Cancer Res ， 2010 ， 16 （ 12 ）： 3130 - 3140 .

姚柳，赵芳，秦玉婷，等 . Shh信号通路在骨髓增生异常综合征中的表达及SMO抑制剂对MUTZ-1细胞作用的研究［J］. 中国现代医学杂志， 2022 ， 32 （ 1 ）： 32 - 39 .

张香红，王荟荟，黄小燕，等 . Wnt信号通路在慢性阻塞性肺疾病中的研究进展［J］. 临床肺科杂志， 2021 ， 26 （ 4 ）： 628 - 631，645 .

BEROUKHIM R ， GETZ G ， NGHIEMPHU L ， et al . Assessing the significance of chromosomal aberrations in cancer：methodology and application to glioma ［J］. Proc Natl Acad Sci USA ， 2007 ， 104 （ 50 ）： 20007 - 20012 .

BEROUKHIM R ， MERMEL CH ， PORTER D ， et al . The landscape of somatic copy number alteration across human cancers ［J］. Nature ， 2010 ， 463 （ 7283 ）： 899 - 905 .

FRANCAVILLA C ， RIGBOLT KT ， EMDAL KB ， et al . Functional proteomics defines the molecular switch underlying FGF receptor trafficking and cellular outputs ［J］. Mol Cell ， 2013 ， 51 （ 6 ）： 707 - 722 .

KIM YM ， KIM DH . dRAGging amino acid-mTORC1 signaling by SH3BP4 ［J］. Mol Cells ， 2013 ， 35 （ 1 ）： 1 - 6 .

郑亚秋，潘琦，梁利，等 . SH3BP4基因对肝癌细胞迁移和增殖作用的实验研究［J］. 重庆医科大学学报， 2018 ， 43 （ 5 ）： 617 - 623 .

查看原文

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

从“体质-证候-病位相关”探讨高血压共病的中医辨治思路

“天癸”与女性生殖遗传疾病相关性探究

外湿因素诱导下银屑病脾虚湿阻证病证结合大鼠模型的建立

冠心病稳定型心绞痛患者再入院中西医风险预测模型构建与应用评估

支气管哮喘不同临床表型的体-病-证特征及三辨诊疗模式的应用