二甲双胍使用史和全基因组DNA甲基化谱:衰老和长寿的潜在分子机制
摘要
背景:二甲双胍是一种常用的抗糖尿病药物,在临床前模型中已反复证明可以延缓衰老,并与降低人类死亡率相关。然而,从生物学的角度来看,二甲双胍如何能潜在地延长寿命尚不清楚。我们假设二甲双胍对长寿的潜在作用机制是通过其表观遗传修饰。
方法:为了验证我们的假设,我们对从有和无二甲双胍使用史的住院患者采集的全血中获得的可用全基因组DNA甲基化(DNAm)数据进行了事后分析。我们评估了171名患者(第一轮)和63名糖尿病患者(第二轮)的甲基化情况,并比较了二甲双胍使用者和非使用者的DNAm比率。
结果:《京都基因与基因组百科全书》(KEGG)的富集分析显示了与二甲双胍作用机制相关的途径,如长寿、AMPK和炎症途径。我们还确定了几种与谵妄相关的途径,其危险因素是衰老。此外,基因本体论(GO)的热门研究包括HIF-1α途径。然而,没有单个CpG位点显示全基因组统计意义(第页<5E-08)。
结论:本研究可能通过表观遗传学修饰和其他可能的作用机制阐明二甲双胍在长寿中的潜在作用。
介绍
我们生活在一个老龄化社会。根据美国人口普查局(U.S.Census Bureau)2017年全国人口预测,到2030年,每5名居民中就有1人将达到退休年龄[1]. 随后,有更大比例的人口将比以往任何时候都能忍受与年龄有关的疾病的挑战。针对这些疾病(如痴呆症或癌症)的治疗最多只能“延缓”疾病进程,但“治愈”能力有限。因此,人们对将衰老本身视为一种疾病越来越感兴趣[2].
来自基础和临床前模型的大量证据表明,一些干预措施,如运动、间歇性禁食,甚至摄入某些化合物,都可以延长寿命。这些有前途的化合物包括雷帕霉素[三,4],白藜芦醇[5–7],北美[8]和二甲双胍[9–11]. 我们小组还证实,与非二甲双胍使用者相比,使用二甲双胍的住院患者三年生存率有所提高[12]. 此外,我们的数据还表明,服用二甲双胍的患者的谵妄患病率低于未服用二甲双酮的患者[12].
已经对这些干预措施延长寿命或潜在延缓衰老的作用机制进行了深入研究。然而,没有确切的过程是清楚的,尤其是二甲双胍。人们认为表观遗传学是动植物衰老最重要的分子机制之一;因此,许多干预措施的“延长寿命”效应可能是通过改变表观遗传过程实现的。例如,一些报告显示运动后的表观遗传变化[13],禁食[14],雷帕霉素[三],白藜芦醇[5]和NAD[8]. 然而,只有少数研究调查二甲双胍对表观遗传变化的直接影响[15–17]这表明,目前关于二甲双胍对人类表观遗传特征影响的信息有限。
为了填补这一知识空白,我们通过检测有和无二甲双胍使用史的住院患者的全血样本中全基因组DNA甲基化(DNAm)来研究二甲双胍的潜在影响。
结果
人口统计学的
173名受试者参与了这项研究,但只有171人参与了下游数据分析。患者的平均年龄为74.4岁(标准偏差=9.8)。58名(33.9%)受试者是女性,而根据自我报告,几乎所有受试者都是白人(n个= 167; 97.7%). 108名患者为非糖尿病(non-DM),63名患者为糖尿病(DM)。在糖尿病组中,37名糖尿病患者有二甲双胍处方糖尿病(+)Met病史,26名糖尿病患者没有二甲双胍-二甲双胍的处方糖尿病(-)Met史。此外,43名(68.3%)糖尿病受试者有胰岛素使用史。查尔森共病指数(CCI)和体重指数(BMI)信息也包括在表1没有变量显示DM(−)Met和DM(+)Met之间的统计学显著差异。然而,正如预期的那样,与非糖尿病组相比,糖尿病组的CCI、BMI和胰岛素使用显著升高。
表1。患者特征。
分类 | 所有受试者 | 糖尿病 | 第页 | 统计检验 | 糖尿病受试者 | 第页 | 统计检验 |
非糖尿病患者 | DM公司 | DM(−)已达到 | DM(+)符合 |
N个 | 171 | 108 | 63 | | | 26 | 37 | | |
年龄-年 | 74.4 | 74.6 | 74.1 | 0.77 | t吨= 1.98 | 73.8 | 74.3 | 0.833 | t吨= 2.01 |
标准偏差 | 9.8 | 9.7 | 10 | | | 10.6 | 9.7 | | |
女性(n个) | 58 | 36 | 22 | 0.81 | χ2= 0.10 | 11 | 11 | 0.303 | χ2= 1.06 |
% | 33.9 | 33.6 | 34.9 | | | 42.3 | 29.7 | | |
种族,白人(n个) | 167 | 105 | 62 | 0.63 | χ2= 0.23 | 25 | 37 | 0.229 | χ2= 1.45 |
% | 97.7 | 97.2 | 98.4 | | | 96.2 | 100 | | |
CCI公司 | 3.8 | 3.1 | 4.9 | 7.5E-06页* | t吨= 1.98 | 4.8 | 5 | 0.756 | t吨= 2.00 |
标准偏差 | 2.7 | 2.7 | 2.4 | | | 2.4 | 2.5 | | |
体重指数 | 29.7 | 28.3 | 32.2 | 0.002* | t吨= 1.98 | 30 | 33.8 | 0.64 | t吨= 2.00 |
标准偏差 | 7.6 | 6.3 | 8.8 | | | 5 | 10.5 | | |
胰岛素使用史 | 43 | 0 | 43 | 第3.3页至第23页* | χ2= 98.48 | 15 | 28 | 0.131 | χ2= 2.28 |
% | 25.1 | 0 | 68.3 | | | 57.7 | 75.7 | | |
年龄、性别和种族在非糖尿病组和糖尿病组之间没有显著差异,而CCI、BMI和胰岛素使用则有显著差异。糖尿病组中二甲双胍非使用者DM(−)Met和二甲双胍者DM(+)Met之间的患者特征无统计学意义。缩写:SD:标准偏差;CCI:查尔森共病指数;BMI:体重指数。*第页< 0.05. |
Met与非Met(包括所有患者,无论糖尿病状态如何):热门、KEGG、GO
表2显示了在有和无二甲双胍使用史的患者(171名受试者)之间甲基化率差异最大的基因。没有一个位点符合全基因组统计显著性的标准(第页<5E-8)。
表2。所有患者中二甲双胍使用者和非使用者之间差异最大的20个CpG位点。
基因名称 | CpG站点 | 染色体 | 未达标(%) | 已达到(%) | %平均差(Δβ) | 第页-价值 |
PSME3型 | cg22769787号 | chr17 | 15.6% | 14.3% | 1.3% | 3.37E-07号 |
EPHA8型 | cg27136384号 | chr1号机组 | 83.2% | −2.7% | −2.7% | 4.84E-07 |
| cg22163972号 | chr17 | 92.1% | 4.2% | 4.2% | 4.89电子-07 |
| cg23047680号 | 氯甲烷 | 0.8% | −0.2% | −0.2% | 9.08E-07年9月 |
NEDD4公司 | cg11341892号 | chr15 | 4.7% | 0.6% | 0.6% | 2.82E-06号机组 |
PRKCG公司 | cg11293016 | chr19系列 | 52.9% | 4.0% | 4.0% | 4.68E-06页 |
SRSF11型 | cg12923877号 | chr1号机组 | 97.5% | −0.3% | −0.3% | 4.94E-06号机组 |
第15页 | cg24353272号 | chr1号机组 | 95.3% | −0.8% | −0.8% | 5.16E-06页 |
科威特航空1688 | cg07969649号 | chr8系列 | 91.1% | −1.6% | −1.6% | 5.22E-06日 |
修剪27 | cg02525926号 | 铬6 | 97.4% | 0.8% | 0.8% | 2006年6月9日 |
| cg23067796号 | chr12 | 93.7% | 1.7% | 1.7% | 2006年7月29日 |
瑞尔2 | cg04573831号 | chr1号机组 | 96.6% | −0.6% | −0.6% | 2006年11月8日 |
| cg15180899号 | chr18 | 93.9% | 1.7% | 1.7% | 8.67E-06号 |
| cg12222244号 | 氯甲烷 | 94.1% | 2.1% | 2.1% | 1.27E-05 |
C1或125 | cg20746459号 | chr1号机组 | 90.6% | 3.5% | 3.5% | 1.52E-05型 |
SERPINH1系列 | cg19586851号 | chr11号机组 | 97.2% | −0.5% | −0.5% | 1.55E-05型 |
PPL(公私合营) | cg12991522号 | chr16系列 | 1.8% | −0.5% | −0.5% | 2005年5月15日 |
ACO1公司 | cg13567378号 | chr9公司 | 89% | −1.3% | −1.3% | 2005年1月17日 |
| cg24525630 | chr17 | 1.6% | −0.3% | −0.3% | 1.72E-05 |
TCF7L1型 | cg20116596号 | chr2型 | 95.7% | −0.5% | −0.5% | 1.76E-05 |
接下来,我们根据二甲双胍使用者(Met)和非使用者(非Met)甲基化水平(β值)的绝对差异大于4%,以及第页-值小于0.01。KEGG顶部信号的富集分析显示了二甲双胍可能作用的相关途径,如“寿命调节途径”、“寿命调节通路-多物种”和“AMPK信号通路”(表3). 此外,还发现了其他通路,如“mTOR信号通路”、“胰岛素分泌”、“谷氨酸能突触”和“昼夜节律夹带”(表3). GO分析还揭示了相关的通路,如“低氧诱导因子-1α信号通路的调节”、“低氧诱发因子-1α的正调节信号通路”和“典型Wnt信号通路”(表4),尽管KEGG或GO中的路径均未达到假发现率(FDR)显著性水平(FDR<0.05)(表3和4).
表3。基于二甲双胍使用者和非使用者之间不同甲基化率的前30种KEGG途径。
通路 | N个 | 判定元件 | 第页-价值 | 财务总监 |
松弛素信号通路 | 129 | 6 | 0.007 | 1 |
长寿调节途径 | 89 | 5 | 0.008 | 1 |
谷氨酸能突触 | 114 | 6 | 0.008 | 1 |
库欣综合征 | 155 | 6 | 0.018 | 1 |
甲状旁腺激素的合成、分泌和作用 | 106 | 5 | 0.019 | 1 |
AMPK信号通路 | 119 | 5 | 0.021 | 1 |
调节干细胞多能性的信号通路 | 142 | 5 | 0.028 | 1 |
缝隙连接 | 88 | 4 | 0.033 | 1 |
胰岛素分泌 | 86 | 4 | 0.034 | 1 |
黑色素生成 | 101 | 4 | 0.043 | 1 |
长寿调节途径-多物种 | 62 | 三 | 0.051 | 1 |
醛固酮的合成和分泌 | 98 | 4 | 0.055 | 1 |
化学致癌作用-DNA加合物 | 69 | 2 | 0.056 | 1 |
昼夜节律夹带 | 97 | 4 | 0.058 | 1 |
类固醇激素生物合成 | 61 | 2 | 0.062 | 1 |
生热作用 | 219 | 5 | 0.063 | 1 |
胆汁分泌 | 89 | 三 | 0.063 | 1 |
细胞色素P450对外源物质的代谢 | 76 | 2 | 0.068 | 1 |
皮质醇合成和分泌 | 65 | 三 | 0.069 | 1 |
甲状腺激素合成 | 75 | 三 | 0.071 | 1 |
Wnt信号通路 | 166 | 5 | 0.071 | 1 |
加压素调节的水重吸收 | 44 | 2 | 0.081 | 1 |
胆碱能突触 | 113 | 4 | 0.090 | 1 |
逆行内源性大麻素信号 | 141 | 4 | 0.089 | 1 |
雌激素信号通路 | 137 | 4 | 0.091 | 1 |
矿物质吸收 | 60 | 2 | 0.111 | 1 |
胃癌 | 149 | 4 | 0.121 | 1 |
mTOR信号通路 | 155 | 4 | 0.123 | 1 |
蛋白质消化吸收 | 102 | 三 | 0.123 | 1 |
卵巢类固醇生成 | 51 | 2 | 0.123 | 1 |
甲状腺激素合成 | 75 | 三 | 0.071 | 1 |
KEGG的相关途径[58]突出显示。缩写:N个:每条通路中包含的基因数量;DE:差异表达基因的数量,即来自CpG顶部位点的基因数量;FDR:错误发现率。 |
表4。基于二甲双胍使用者和非使用者之间不同甲基化率的前30种GO途径。
通路 | 安大略省 | N个 | 判定元件 | 第页-价值 | 财务总监 |
通过质膜粘附分子实现嗜同性细胞粘附 | 英国石油公司 | 168 | 8 | 2004年7月20日 | 1 |
长期突触抑制 | 英国石油公司 | 31 | 4 | 9.77E-04号 | 1 |
运动行为 | 英国石油公司 | 198 | 9 | 0.001 | 1 |
中脑多巴胺能神经元分化 | 英国石油公司 | 17 | 三 | 0.002 | 1 |
细胞表面受体信号通路参与细胞间信号传导 | 英国石油公司 | 622 | 17 | 0.002 | 1 |
突触传递的负调控 | 英国石油公司 | 71 | 5 | 0.002 | 1 |
典型Wnt信号通路 | 英国石油公司 | 335 | 11 | 0.002 | 1 |
钙离子结合 | MF公司 | 698 | 17 | 0.003 | 1 |
己糖介导的信号传导 | 英国石油公司 | 6 | 2 | 0.003 | 1 |
糖介导的信号通路 | 英国石油公司 | 6 | 2 | 0.003 | 1 |
葡萄糖介导的信号通路 | 英国石油公司 | 6 | 2 | 0.003 | 1 |
细胞对酸性化学物质的反应 | 英国石油公司 | 209 | 8 | 0.003 | 1 |
细胞间信号 | 英国石油公司 | 1847 | 345 | 0.003 | 1 |
离子跨膜转运蛋白活性的调节 | 英国石油公司 | 256 | 9 | 0.004 | 1 |
中胚层发育 | 英国石油公司 | 133 | 6 | 0.004 | 1 |
神经细胞体膜 | 科科斯群岛 | 27 | 三 | 0.004 | 1 |
细胞体膜 | 科科斯群岛 | 28 | 三 | 0.004 | 1 |
跨膜转运蛋白活性的调节 | 英国石油公司 | 264 | 9 | 0.005 | 1 |
低氧诱导因子-1α信号通路的调节 | 英国石油公司 | 1 | 1 | 0.005 | 1 |
低氧诱导因子-1α信号通路的正调控 | 英国石油公司 | 1 | 1 | 0.005 | 1 |
细胞对维生素K的反应 | 英国石油公司 | 1 | 1 | 0.005 | 1 |
细胞对胰高血糖素刺激的反应 | 英国石油公司 | 25 | 三 | 0.005 | 1 |
碳水化合物介导的信号 | 英国石油公司 | 8 | 2 | 0.005 | 1 |
精囊形态发生 | 英国石油公司 | 1 | 1 | 0.005 | 1 |
胰高血糖素样肽1受体活性 | MF公司 | 1 | 1 | 0.005 | 1 |
行为 | 英国石油公司 | 593 | 15 | 0.005 | 1 |
烟酰胺磷酸核糖转移酶活性 | MF公司 | 1 | 1 | 0.006 | 1 |
对D-半乳糖的反应 | 英国石油公司 | 1 | 1 | 0.006 | 1 |
胚胎骨骼系统发育 | 英国石油公司 | 125 | 6 | 0.006 | 1 |
运输活动的监管 | 英国石油公司 | 279 | 9 | 0.006 | 1 |
重点介绍了相关途径。缩写:Ont:本体;BP:生物过程;CC:细胞成分;MF:分子功能;N个:每条通路中包含的基因数量;DE:差异表达基因的数量,即来自CpG顶部位点的基因数量;FDR:错误发现率。 |
Met与非Met(仅包括2型糖尿病患者):热门、KEGG、GO
表5显示了糖尿病组(63名受试者)中二甲双胍使用者和非使用者甲基化率差异最大的基因。与之前的分析类似,没有基因达到全基因组的统计显著性(第页<5E-8)。
表5。糖尿病组二甲双胍使用者和非使用者之间差异最大的20个CpG位点。
基因名称 | CpG站点 | 染色体 | 未达标(%) | 已达到(%) | 平均差(Δβ) | 第页-价值 |
| cg19873536号 | 铬10 | 78.3% | 67.9% | 10.4% | 1.28E-06 |
| cg13596208号 | chr9公司 | 1.9% | 2.7% | −0.9% | 2006年2月29日 |
HBA1型 | cg01704105号 | chr16系列 | 40.5% | 33.7% | 6.8% | 2006年5月42日 |
DUOX2型 | cg02550961号 | 第15章 | 1.5% | 1.9% | −0.4% | 2006年6月6日 |
近地天体1 | cg12516231号 | 第15章 | 2.2% | 3.2% | −0.9% | 6.97E-06号 |
C7或46 | cg06685724号 | chr7系列 | 2.1% | 2.9% | −0.8% | 1.28E-05 |
NAT15型 | cg00484396号 | chr16系列 | 9.8% | 4.9% | 4.8% | 1.56E-05型 |
| cg14685975 | 铬5 | 89.9% | 92.1% | −2.2% | 1.64E-05 |
CTSL公司 | cg02104500型 | chr9公司 | 3.6% | 4.9% | −1.4% | 1.66E-05 |
| cg12584257 | chr9公司 | 67.6% | 77.2% | −9.6% | 1.69E-05 |
NAT15型 | cg22508957号 | chr16系列 | 10.9% | 6.3% | 4.6% | 1.84E-05型 |
区域1 | cg11034672号 | chr14系列 | 11.6% | 15% | −3.3% | 1.86E-05型 |
| cg24651265号 | 铬10 | 1.1% | 1.7% | −0.5% | 2.12E-05日 |
CMBL公司 | cg17467873 | 铬5 | 1.7% | 2.1% | −0.4% | 2.21E-05年 |
EBF4类 | cg05857996号 | 铬20 | 77.6% | 63.6% | 13.9% | 2.23E-05日 |
| cg18482666号 | chr2型 | 95.8% | 94.8% | 1.0% | 2.39E-05号 |
HRASLS5型 | cg00489394号 | chr11号机组 | 6.6% | 7.1% | −0.5% | 2.40至05 |
AKAP13公司 | cg21530087号 | chr15 | 2.2% | 2.6% | −0.4% | 第2.59页至第05页 |
| cg15864571号 | 氯甲烷 | 93.4% | 95.0% | −1.6% | 2.67E-05年 |
FLJ35024型 | cg15981195 | chr9公司 | 2.3% | 3.5% | −1.1% | 2.91E-05年 |
使用甲基化水平差异的一致参数(β>4%)和第页-值(<0.01)。然而,目前的分析包括1283个CpG。KEGG显示了许多从先前分析中发现的相同信号,包括“寿命调节通路”、“谷氨酸能突触”、“胰岛素分泌”、“昼夜节律携带”和“胆碱能突触“(表6). 与第一次分析相比,GO也显示出重叠的通路,包括“低氧诱导因子-1α信号通路”,还有新的通路,如“白细胞介导的信号通路”、“白细胞凋亡过程的负调控”、“中性粒细胞内稳态”和“神经元投射”,尽管这些通路没有达到FDR显著性水平(FDR<0.05)(表7).
表6。糖尿病组中二甲双胍使用者和非使用者之间存在差异的前30种KEGG途径。
通路 | N个 | 判定元件 | 第页-价值 | 财务总监 |
醛固酮的合成和分泌 | 98 | 14 | 0.001 | 0.219 |
昼夜节律夹带 | 97 | 14 | 0.001 | 0.219 |
皮质醇合成和分泌 | 65 | 10 | 0.003 | 0.303 |
甲状腺激素合成 | 75 | 10 | 0.004 | 0.303 |
脂肪细胞脂肪分解的调节 | 55 | 8 | 0.006 | 0.330 |
甲状旁腺激素的合成、分泌和作用 | 106 | 13 | 0.006 | 0.330 |
胰岛素分泌 | 86 | 11 | 0.007 | 0.330 |
钙信号通路 | 238 | 21.5 | 0.009 | 0.388 |
cAMP信号通路 | 221 | 19 | 0.010 | 0.388 |
胆碱能突触 | 113 | 13 | 0.012 | 0.420 |
化学致癌-受体激活 | 212 | 16 | 0.015 | 0.435 |
谷氨酸能突触 | 114 | 13 | 0.016 | 0.435 |
Rap1信号通路 | 210 | 19 | 0.016 | 0.435 |
生热作用 | 219 | 15 | 0.020 | 0.468 |
苯丙胺成瘾 | 69 | 8 | 0.021 | 0.468 |
神经活性配体-受体相互作用 | 349 | 19.5 | 0.022 | 0.468 |
胰腺分泌物 | 101 | 9 | 0.029 | 0.552 |
长期增强 | 67 | 8 | 0.029 | 0.552 |
可卡因成瘾 | 49 | 6 | 0.036 | 0.552 |
磷脂酶D信号通路 | 147 | 14.5 | 0.038 | 0.552 |
cGMP-PKG信号通路 | 166 | 14.5 | 0.039 | 0.555 |
Apelin信号通路 | 139 | 12 | 0.040 | 0.555 |
尼古丁成瘾 | 40 | 5 | 0.040 | 0.555 |
EGFR酪氨酸激酶抑制剂耐药性 | 78 | 9 | 0.041 | 0.555 |
TRP通道的炎症介质调节 | 98 | 10 | 0.042 | 0.555 |
缝隙连接 | 88 | 9 | 0.042 | 0.555 |
II型糖尿病 | 46 | 6 | 0.043 | 0.558 |
长寿调节途径 | 89 | 9 | 0.045 | 0.561 |
唾液分泌 | 92 | 8 | 0.048 | 0.578 |
膀胱癌 | 41 | 5 | 0.054 | 0.610 |
KEGG的相关途径[58]突出显示。缩写:N个:每条通路中包含的基因数量;DE:差异表达基因的数量,即来自CpG顶部位点的基因数量;FDR:错误发现率。 |
表7。糖尿病组中二甲双胍使用者和非使用者之间存在差异的前30种GO途径。
通路 | Ont公司 | N个 | 判定元件 | 第页-价值 | 财务总监 |
神经元投射 | 科科斯群岛 | 1304 | 97.1 | 1.29E-05 | 0.286 |
第二信使介导信号 | 英国石油公司 | 438 | 38.5 | 第2.53页至第05页 | 0.286 |
中性粒细胞稳态 | 英国石油公司 | 16 | 6 | 3.77E-05日 | 0.286 |
突触信号 | 英国石油公司 | 725 | 59.5 | 2005年9月26日 | 0.505 |
跨突触信号 | 英国石油公司 | 708 | 57.5 | 1.67E-04页 | 0.505 |
白细胞凋亡过程的负调控 | 英国石油公司 | 46 | 8 | 1.93E-04号机组 | 0.505 |
钙介导信号 | 英国石油公司 | 218 | 22.5 | 1.94E-04号机组 | 0.505 |
化学突触传递 | 英国石油公司 | 700 | 56.5 | 2.00E-04年 | 0.505 |
顺行跨突触信号 | 英国石油公司 | 700 | 56.5 | 2.00E-04年 | 0.505 |
细胞基质粘附的正向调节 | 英国石油公司 | 51 | 10 | 2004年3月1日 | 0.682 |
多细胞生物过程的正向调节 | 英国石油公司 | 1802 | 106 | 3.39E-04号 | 0.682 |
质膜结合细胞投影 | 科科斯群岛 | 2093 | 130.1 | 4.05E-04号机组 | 0.682 |
白细胞介素-8受体活性 | MF公司 | 2 | 2 | 4.44E-04号机组 | 0.682 |
白细胞介素-8介导的信号通路 | 英国石油公司 | 2 | 2 | 4.44E-04号机组 | 0.682 |
成人行为 | 英国石油公司 | 144 | 17.5 | 4.73E-04号机组 | 0.682 |
细胞连接 | 科科斯群岛 | 1858 | 123.8 | 5.16E-04页 | 0.682 |
突触 | 科科斯群岛 | 1168 | 85.5 | 5.27E-04号 | 0.682 |
NMDA谷氨酸受体活性 | MF公司 | 7 | 4 | 6.05E-04(电子版) | 0.682 |
低氧诱导因子-1α信号通路 | 英国石油公司 | 6 | 三 | 6.39电子04 | 0.682 |
枝晶发育的调控 | 英国石油公司 | 148 | 19 | 6.44E-04号 | 0.682 |
阿克森 | 科科斯群岛 | 606 | 50.6 | 6.46E-04号 | 0.682 |
低电压门控钙通道活性 | MF公司 | 三 | 三 | 04年7月18日 | 0.682 |
树枝状岩发育 | 英国石油公司 | 232 | 26.5 | 2014年7月28日 | 0.682 |
前庭耳蜗神经发育 | 英国石油公司 | 10 | 4 | 7.64E-04页 | 0.682 |
离子型谷氨酸受体信号通路 | 英国石油公司 | 25 | 7 | 7.69E-04页 | 0.682 |
神经元投射发育 | 英国石油公司 | 976 | 74 | 8.82E-04号 | 0.682 |
细胞对葡萄糖刺激的反应 | 英国石油公司 | 132 | 15 | 2004年9月26日 | 0.682 |
运动行为 | 英国石油公司 | 198 | 21.5 | 9月36日-04 | 0.682 |
细胞对己糖刺激的反应 | 英国石油公司 | 134 | 15 | 1.08E-03版 | 0.682 |
细胞成分生物发生的正调控 | 英国石油公司 | 533 | 41 | 1.12E-03号机组 | 0.682 |
强调了相关途径。缩写:Ont:本体;BP:生物过程;CC:细胞成分;MF:分子功能;N个:每条通路中包含的基因数量;DE:差异表达基因的数量,即来自CpG顶部位点的基因数量;FDR:错误发现率。 |
DNA甲基化年龄加速
在糖尿病组中,非二甲双胍使用者的平均年龄加速为−8.07,而二甲双胍者的平均年龄加快为−4.47(第页= 0.11) (图1). 无论糖尿病状态如何,所有受试者之间的差异都较小(未使用二甲双胍的受试者为-5.92,使用二甲双双胍者为-4.47;第页= 0.34) (图2). 两项分析均未达到统计显著性。
图1。 糖尿病组二甲双胍使用者和非使用者之间的年龄加速。使用Horvath表观遗传时钟作为DNAm年龄-年代学年龄计算年龄加速。二甲双胍=0:无二甲双胍症使用史,二甲双胍=1:有二甲双胍血症使用史。第页= 0.11.
图2。 二甲双胍使用者和非使用者之间的年龄加速。使用Horvath表观遗传时钟作为DNAm年龄-年代学年龄计算年龄加速。二甲双胍=0:无二甲双胍症使用史,二甲双胍=1:有二甲双胍血症使用史。第页= 0.34.
讨论
在本研究中,我们比较了二甲双胍使用者和非使用者的全基因组DNA甲基化率,以研究二甲双胍的潜在表观遗传效应。采用富集分析阐明二甲双胍诱导的可能作用机制。我们的KEGG分析揭示了与“长寿”有关的表观遗传图谱差异的证据,例如“长寿调节途径”和“长寿调节路径-多物种”(表3和6). 虽然没有统计学意义,但KEGG分析中顶层信号中这些通路的出现证明了表观遗传过程在显示二甲双胍对寿命的影响方面的潜在作用。同样的KEGG分析也显示“AMPK信号通路”(表3). AMP-activated protein kinase(AMPK)是一种调节新陈代谢的能量传感器,通常被认为是二甲双胍假设作用机制的靶点之一[18,19],尽管也有证据表明二甲双胍的作用在一定程度上与AMPK无关[20]. 此外,AMPK激活与随后低氧诱导因子的激活有关[21]在我们的GO分析中也出现了“低氧诱导因子-1α信号通路的调节”和“低氧诱发因子-1α的信号通路的正向调节”(表4)以及“低氧诱导因子-1α信号通路”(表7). 低氧诱导因子-1α(HIF-1α)是一种在有核细胞中表达并由氧水平介导的转录因子。HIF-1α与年龄相关疾病、内皮细胞衰老进程、AMPK和许多其他途径有关[22]. 除了二甲双胍潜在的与长寿相关的表观遗传药物外,还发现了一些与谵妄相关的途径,如“昼夜节律夹带”、“胆碱能突触”和“谷氨酸能突触(表3和6). 从二甲双胍可能的“抗衰老”观点来看,这些途径很有趣,因为年龄是谵妄的主要危险因素。
二甲双胍对寿命的有益影响已被广泛研究。先前的研究报告称二甲双胍增加了秀丽线虫与…共同培养大肠杆菌超过35%[9,23],并延长了小鼠的寿命[10]. 服用二甲双胍的心血管疾病和癌症等年龄相关疾病患者的死亡率也较低[24,25]. 我们最近对1400多名住院患者进行的研究还表明,有二甲双胍使用史的糖尿病患者的3年死亡率明显低于未服用二甲双胍的糖尿病患者[12]. 然而,也有相互矛盾的报道。例如,在果蝇[26]. 此外,还报告了年龄依赖性、剂量依赖性和性别依赖性变量对小鼠寿命的影响[27,28]. 尽管这些先前的研究结果并不一致,但我们上面提到的队列(目前的数据是对其亚组的分析)清楚地表明二甲双胍的使用对糖尿病住院患者的生存率有积极影响[12].
我们在此提供的表观遗传学数据支持二甲双胍的广泛潜在作用,正如富集分析确定的途径所表明的那样。KEGG分析(表7)显示了一些与炎症和免疫系统相关的信号,如“白细胞介素-8受体活性”和“白细胞凋亡过程的负调控”。炎症相关通路的出现令人感兴趣,因为有强有力的证据表明老年人存在低度慢性炎症[29]. 我们研究中发现的这些信号可能支持我们的假设,即二甲双胍可以通过表观遗传修饰改变炎症过程并影响生存可能性。与我们的数据一致,Barath等人还报道,二甲双胍通过纠正自噬和线粒体生物能量学中与年龄相关的变化,抑制Th17产生细胞因子,表明其药物促进健康衰老的潜力[30]. 在支持二甲双胍抑制炎症作用的文献中,包括糖尿病预防计划(DPP)在内的临床试验[31]和旁路血管成形术再血管化研究2型糖尿病(BARI 2D)[32]进一步证明了二甲双胍在改变糖尿病患者炎症生物标志物水平方面的作用,而其他临床试验,如Lantus for C-反应蛋白减少早期治疗2型糖尿病(LANCET)[33],发现了相反的证据。虽然这里提到的几项研究已经调查了二甲双胍及其潜在抗炎作用之间的关系,但一项旨在证实二甲双双胍在衰老中作用的临床试验尚未见到[2,34]. 尽管如此,值得一提的是,一项小型临床研究表明,通过注射重组人生长激素(rhGH)、脱氢表雄酮(DHEA)和二甲双胍,患者的表观遗传年龄会下降[15]. 由于研究团队同时给受试者服用三种药物,因此无法区分仅由二甲双胍引起的表观遗传变化。值得一提的是,Horvath表观遗传时钟的意外结果,因为有二甲双胍使用史的受试者的年龄加速率相对高于无二甲双胍的受试对象。尽管如此,两者都没有达到统计显著性(第页< 0.05). 未来的前瞻性研究需要比较二甲双胍使用前后的表观遗传学标记,以更好地了解药物的直接作用。
仅在DM受试者中,A激酶锚定蛋白13(AKAP13公司)基因被发现(表5). 最近的一项研究表明,AKAP13抑制雷帕霉素复合物1(mTORC1)的哺乳动物靶点,该靶点在我们的富集分析中被称为“mTOR信号通路”(表3). 此外AKAP13公司肺腺癌细胞系中的表达与mTORC1活性相关[35]. 二甲双胍的抗炎作用是通过最终激活AMPK来实现的,AMPK也抑制mTOR信号通路[18]. 二甲双胍与AKAP13的关系尚待充分了解,值得进一步调查。
据我们所知,我们的研究是同类研究中规模最大的。此前的一项规模较小的研究也调查了二甲双胍对人类外周血全基因组DNA甲基化的影响,尽管他们的研究能力仅限于32名男性受试者的样本量[36]. 本研究中的丰富分析从无假设的方法揭示了长寿途径,进一步强化了我们的假设,即二甲双胍通过表观遗传过程显示其对长寿的潜在益处。我们还确定了与二甲双胍作用机制相关的其他相关通路,如AMPK信号通路和HIF-1α信号通路[37].
我们的研究有几个局限性。虽然在本研究中对171名受试者进行了回顾性分析,但更大样本量的对照前瞻性研究将更好地了解二甲双胍对长寿的表观遗传机制。此外,没有一个单独的CpG位点达到全基因组意义(第页<5E-08)。因此,我们的发现应该被解释为探索性和假设生成性。然而,我们发现它们与二甲双胍作用的生物学相关性仍然值得注意。由于受试者的糖尿病和二甲双胍使用状况是根据电子病历的回顾性图表审查确定的,尽管我们仍然能够找到与二甲双胍的作用相关的多个相关途径和相关基因,但仍有可能出现错误分类。此外,二甲双胍使用的持续时间没有得到精确评估,这使得我们对“二甲双胍史使用”的定义很宽泛,因为它可能包括仅服用二甲双胍几个月的患者和服用二甲双胍数年的患者。此外,其他类型的糖尿病药物,如磺酰脲类和格列奈类药物,也没有进行研究,因为我们使用了我们之前工作中已经完成的研究数据集。我们没有调查其他糖尿病药物的影响的理由是基于过去的文献,这些文献表明,二甲双胍以外的糖尿病药物对生存没有益处。事实上,有时它们会导致更严重的死亡率[38–40].
总之,这里提供的数据支持我们的假设,即表观遗传学,特别是DNA甲基化,可能会因二甲双胍的使用而改变,并且这种表观遗传学过程可能有助于导致长寿的分子机制。有必要对更大的样本量进行进一步的仔细调查。
方法
研究参与者和招募
我们之前在爱荷华州大学医院和诊所(UIHC)招募了患者,从2016年1月至2020年3月进行了一项与谵妄相关的单独研究[41–44]. 其中,我们使用了2017年11月至2020年3月招募的一组患者的数据,这些患者采集并处理了血样进行表观遗传学分析[45–47]. 与急诊科、骨科、普通医学科或重症监护室收治的18岁或以上患者进行接触。只有那些同意或其合法授权代表同意的人才被纳入研究。获得了所有参与者的书面知情同意。排除标准包括仅以舒适措施为护理目标的受试者、囚犯或采取防滴/接触措施的个人。有关研究对象和注册过程的更多详细信息,请参见上文[41–44].
我们测试了173名受试者的全基因组DNA甲基化(DNAm)状态,然后对可用数据进行事后分析,以评估二甲双胍的影响。该研究由爱荷华大学医院和诊所机构审查委员会批准,所有程序均符合赫尔辛基宣言。
临床信息
通过电子病历审查、患者访谈和家庭成员的附带信息收集临床变量[41–44]. 分别使用搜索词“二甲双胍”、“胰岛素”和“DM”或“糖尿病”获得二甲双胍的使用情况、胰岛素使用情况和2型糖尿病(DM)病史[12]. 只包括2型糖尿病(DM),不包括1型糖尿病或妊娠期糖尿病。如果在研究登记前有二甲双胍处方史,则患者被归类为二甲双胍(Met)使用者。参与后服用二甲双胍的患者未被归类为二甲双酮使用者(非二甲双酮),因为血液是在服用该处方之前获得的。
样品采集
在患者住院期间,在EDTA管中采集血样。样品被运送到研究实验室,并在−80°C下储存,直到下游分析时作为一批。
样品分析
使用MasterPure™DNA纯化试剂盒(Epicenter,MCD 85201)从全血中提取DNA。选择基于NanoDrop光谱法并通过Qubit dsDNA宽范围检测试剂盒(ThermoFischer Scientific,Q32850)进行足够数量的DNA通过质量控制,以分析全基因组DNAm状态。使用EZ DNA甲基化™试剂盒(Zymo Research,D5002)对每个样品中的500 ng基因组DNA进行亚硫酸氢盐转化,并使用Infinium HumanMethylationEPICBeadChip™试剂盒进行分析(Illumina,WG-317-1002)。Illumina iScan平台扫描了阵列。
统计学和生物信息学分析
所有分析均使用R进行。R包ChAMP[48]和minfi[49]用于处理数据。来自173名受试者的175份样本的数据被纳入统计和生物信息学分析。首先比较有和没有二甲双胍处方史的患者的每个CpG位点的DNAm水平(首次运行;补充表1). 然后,在仅有糖尿病患者和无二甲双胍处方史患者之间进行有限的比较,以避免糖尿病对DNAm状态的潜在影响(第二次运行;补充表2).
在质量控制过程中,根据密度分析图(作为我们质量控制管道的一部分),首次运行的2个样品和第二次运行的无样品被排除在外。2个样本也被排除在外,因为两名患者的血液被采集了两次。将第一次收集的样本纳入进一步分析,而将第二次样本排除在外,以保持所有受试者样本之间的一致性。因此,第一次试验的171名受试者和第二次试验的63名受试者仍有待分析。此外,在数据加载过程中,如果(i)检测到探针,则会将其过滤掉第页-值>0.01,(ii)每个探针至少5%的样本中有<3个珠子,(iii)非CpG、SNP相关或多点探针,或(iv)位于X或Y染色体上[50]用于对样本进行归一化,而战斗归一化方法用于在第一次运行中校正批量效应[51,52]. 第二次试验只包括糖尿病患者,但没有对批次效应进行校正,因为有个别患者不属于任何批次。
使用limma方法通过RnBeads包获得基于每个CpG位点差异的热门点击数[53,54]并将年龄、性别、胰岛素使用、BMI和细胞类型比例(CD8 T细胞、CD4 T细胞、自然杀伤细胞、B细胞和单核细胞)作为协变量。DNAm年龄计算器在线提供[55]通过前面报告的方法计算细胞类型比例[56].
在获得顶部CpG位点后,使用MismMethyl进行富集分析[57]使用EPIC阵列控制每个基因上CpG位点的不平衡数量。基因本体论(GO)和京都基因和基因组百科全书(KEGG)[58]进行了分析。分析中包含的CpG位点的数量由以下因素组合确定第页-每个CpG位点甲基化率的值和β值截止值(第页<.01和β>0.04)。全基因组重要性设定为第页-值小于=5.0E-08。
齐方检验比较了两组之间的分类数据(性别、种族和胰岛素使用),而韦尔奇的t吨-测试比较了两组之间的数字数据(年龄、BMI和CCI)。
DNA甲基化老化时钟分析
为了研究有二甲双胍使用史的受试者是否比没有二甲双胍的受试对象“年龄加速”较慢,我们将原始DNA甲基化β值提交给一个公开可用的工具,其中包括Horvath[55]方法。计算的结果是DNA甲基化年龄和时间年龄之间的差异。
数据材料的可用性
本研究期间分析的数据集可根据合理要求从相应作者处获得。
缩写
DM:糖尿病;非Met:无二甲双胍使用史的患者;Met:有二甲双胍使用史的患者;DM(−)Met:无二甲双胍使用史的糖尿病患者;糖尿病(+)Met:有二甲双胍使用史的糖尿病患者。
作者贡献
P.S.M.收集、分析数据并撰写手稿。T.Y.整理了临床数据集并编辑了手稿。K.J.C、Z.E.M.A.、M.M.、G.C.和T.T.收集临床数据和生物样本,并进行处理。N.E.W、K.J.C.、M.I.和H.R.C.对手稿进行了批判性审查。G.S.构思了该研究的想法,计划了其设计和协调,并编写和编辑了手稿。
致谢
作者感谢参与这项研究的患者。
利益冲突
Gen Shinozaki是Predelix Medical LLC的联合创始人,拥有以下未决专利:“用于预测和筛查谵妄的非侵入性设备”、PCT申请号PCT/US2016/064937和美国临时专利号62/263325;“利用新型脑电图设备预测患者预后”,美国临时专利第62/829411号;PCT申请号PCT/US19/51276和美国临时专利号62/731599中的“谵妄风险的表观遗传生物标记物”。Pedro S.Marra、Takehiko Yamanashi、Kaitlyn J.Crutchley、Nadia E.Wahba、Zoe-Ella M.Anderson、Manisha Modukuri、Gloria Chang、Tammy Tran、Masaaki Iwata和Hyunkeun Ryan Cho已声明不存在利益冲突。
道德声明和同意
该研究由爱荷华大学医院和诊所机构审查委员会批准,所有程序均符合赫尔辛基宣言。获得了所有参与者的书面知情同意。
工具书类