机械老化发展。作者手稿;PMC 2008年4月7日提供。
以最终编辑形式发布为:
PMCID公司:项目经理2288694
美国国立卫生研究院:美国国立卫生研究院27801
分析人-黑猩猩同源基因对以探索衰老的进化假说
马萨诸塞州波士顿哈佛医学院遗传学系02115
*通讯作者:Joáo Pedro de Magalháes。哈佛医学院遗传学系,路易斯·巴斯德大道77号,马萨诸塞州波士顿238室,邮编02115。电话:1-617-432-6512;传真:1-617-432-6513;电子邮件:最小ecnecsence@pj - 补充资料
1
GUID:76DA9867-7498-4C51-9755-E8ADA75708CD
摘要
与黑猩猩相比(黑猩猩)在人类中,衰老的开始似乎被推迟了。在此,我们研究了人类-黑猩猩同源基因对,以研究在不同模型系统中作用于与衰老相关基因的选择力,这使我们能够探索衰老的进化假说。我们的结果表明,衰老相关基因往往处于净化选择和强于平均水平的功能约束之下。我们几乎没有发现原始人类或人类谱系中与衰老相关的基因加速进化的证据,以前与衰老有关的途径在人类和黑猩猩之间基本上是保守的。特别是,与哺乳动物衰老相关的基因相比,非哺乳动物模型生物和细胞系统中与衰老相关的基因组似乎受到更强的功能约束。沃纳综合征基因可能在原始人类中经历了快速进化。总的来说,我们的发现为模型系统中作用于与衰老相关基因的进化力提供了新的见解。我们认为,模型生物中与衰老相关的基因可能是与衰老的多效性效应相关的保守路径的一部分,而这些多效性影响可能不会调节衰老中的物种差异。
关键词:比较基因组学、生命史、寿命、自然选择
2.材料和方法
2.1. 数据的来源和类型
我们使用了黑猩猩测序和分析协会收集的数据集(黑猩猩测序和分析协会,2005年),由13454个人类-黑猩猩同源基因对加上7043个人类、黑猩猩、小鼠和大鼠同源基因的4向基因比对组成。这些基因比对的直系学被认为是明确的。K一/K秒每个比对的比率之前由黑猩猩测序和分析协会(2005年)以及谱系特异性蛋白质变化的估计。然而,考虑到人类和黑猩猩基因之间的相似性,同义替换将很少,而且也经常是零。为了避免这个问题,使用了另一种对背景突变率的估计K我,这是根据局部基因间/内含子替换率计算得出的。K一/K我比率的解释与K一/K秒比率,以及K一/K我在查看单个基因时,比率为假阳性提供了额外的控制。平均值K一/K秒人-黑猩猩同源基因对的比率为0.23。平均值K一/K我比率也为0.23,4.4%的基因具有K一/K我>1和29%的人类和黑猩猩同源基因编码相同的蛋白质。对于四向比较,平均值K一/K秒原始人的比率为0.20,而鼠的比率为0.13(黑猩猩测序和分析协会,2005年). 我们还使用了Clark等人(2003)它由7645条人-黑猩猩-鼠排列组成,用于区分黑猩猩谱系中的选择和人类谱系中选择,并试图推断人类谱系的加速进化。第页-计算中模型2的值Clark等人(2003),给出了拒绝人类血统中无阳性选择的无效假设的似然比检验,与黑猩猩测序和分析协会(2005年).
恒河猴(猕猴)这是一种约2500万年前从猿类分化而来的旧大陆猴子,寿命长达40年,其衰老发生的年龄比人类或黑猩猩要早得多;老鼠(小家鼠)和老鼠(褐家鼠),大约8000万年前从灵长类进化而来的啮齿动物,寿命不到5年,比灵长类要快得多(芬奇,1990年;de Magalhaes,2006年). 作为一项附加测试,从ENSEMBL v.37下载了人-恒河猴、人-小鼠和人-大鼠同源基因对(Birney等人,2006年).K一/K秒只获得了我们有最高置信度的基因的比率,并且我们排除了伪基因和注释错误的基因对。(尽管如此,对于任何这些数据集,根本不消除基因对都不会显著影响我们的总体结果和结论。)总的来说,9857个人类恒河猴基因的平均值K一/K秒比率为0.20,12063人-小鼠基因的平均值为K一/K秒比率为0.13,对于11594个人鼠基因,平均值为K一/K秒比值为0.13。尽管可能,我们也雇佣了猩猩(蓬戈·皮格马乌斯)这些基因可以推断人类和黑猩猩蛋白质之间的差异是否是由于人类或黑猩猩的特定变化。
2.3. 计算平均值K一/K秒基因类别比率
这个K一和K秒GO类别的值取自黑猩猩测序和分析协会(2005年).计算K一和K秒其他基因类别或组的值,如GenAge或GenAge子集,我们采用了黑猩猩测序和分析协会(2005年)计算K一/K秒每个GO类别的比率。的值K一和K秒为给定类别中的所有基因串联C类计算方法如下:
存在一我和A类我非同义替换和位置的数量以及秒我和秒我每个基因的同义替换数和位点我在类别中C类.
2.4. 检测快速和缓慢进化的基因类别
我们使用了黑猩猩测序和分析协会(2005年)快速识别(即具有较高的平均值K一/K秒比率)和缓慢(即低于平均值K一/K秒比率)不断发展的GO类别。为了确定可能经历快速或缓慢进化的其他基因类别,我们还采用了这些作者使用的方法。二项式概率通常用于计算几个连续事件的概率,每个事件都有两个可能的结果。其中,事件是同义或非同义的替代。假设为二项式分布,我们计算二项式概率如下(黑猩猩测序和分析协会,2005年).
概率第页A类在给定类别中发生的非同义替换事件C类根据非同义替换的预期比例估算,使用:
哪里k个一和k个秒代表平均值K一和K秒所有人类黑猩猩基因对的值。因此,对于类别中的替换事件总数C类,概率第页C类我们观察到相等或更高数量的非同义替换计算如下:
存在一C类和秒C类范畴中非同义和同义替换的总数C类同样,我们通过计算观察到等量或更少数量的非同义替换的概率来确定缓慢发展的类别。
这种方法使我们能够识别具有K一/K秒比率明显高于或低于平均水平。二项式概率不允许直接拒绝特定类别中无加速度或无约束的零假设。相反,它是一种用于检测可能正在经历快速或缓慢演化的最极端异常值的度量。根据黑猩猩测序和分析协会(2005年),其截止值设置为0.001。测试统计值<0.001,发现98个GO类别的特征是升高K一/K秒但很可能只有30个。这意味着,与全基因组平均值相比,这些类别中的大多数(但不是全部)经历了快速进化。同样,251个GO类别的特征更低K一/K秒偶然预计只有32个(黑猩猩测序和分析协会,2005年).
2.5. 基因组覆盖率
因为黑猩猩测序和分析协会(2005年)源于黑猩猩基因组的草图序列,并非所有人类基因都与黑猩猩基因组中的同源基因相匹配,因此我们的结果不应被视为全面的。例如,在GenAge的242个蛋白质编码基因中,211个可以被正确分析。同样,在GenAge的242个蛋白编码基因中,我们获得了169个可以正确分析的人-恒河猴、216个人-小鼠和205个人-大鼠基因对。因此,我们的分析可能忽略了潜在感兴趣的基因,但即使有完整的结果,具有大量基因的群体(如GO或GenAge类别)的趋势也不太可能改变。
2.6. 计算和统计分析
为了检索感兴趣的基因并分析我们的结果,我们设计并实现了一个用Perl编写的自动化模块系统。该系统是老龄化研究计算工具(ARCT)的一部分,如de Magalhaes等人,2005年)并用于检索、解析和数据挖掘本工作中使用的不同数据集:黑猩猩测序和分析协会(2005年)、GenAge数据集和Clark等人(2003)。然后系统生成一组新的结果,例如K一/K秒每个GenAge类别的比率,可以进行统计分析。
如上所述,我们采用了黑猩猩测序和分析协会(2005年)它们也是在类似工作中使用的那些(例如。,尼尔森等人,2005年). 曼希特尼单位测试(MWU)用于比较不同组的基因。当将基因子集与整个数据集进行比较时,我们从数据集中排除了基因子集,以避免非独立性问题。尽管有必要,我们还是用随机选择的基因对进行了计算机模拟,以推断偶然观察我们结果的可能性。每次模拟涉及100000个随机排列的数据集。使用SPSS软件包11.5版(SPSS Inc.,芝加哥伊利诺伊州)进行统计分析。
3.结果
3.1. 作用于推测的衰老相关基因的选择力
为了研究模型系统中作用于与衰老相关基因的整体选择力,我们首先比较了K一/K秒这些基因和通路与K一/K秒所有人类黑猩猩同源基因对的比率。通常,净化选择在人类基因中占主导地位K一/K秒人-黑猩猩直系亲属的比率为0.23(黑猩猩测序和分析协会,2005年). 基因组显著高于K一/K秒比率是快速进化的假定情况,因此更可能与表型变化和副病毒有关。因此,我们试图寻找积极选择的证据,或者至少是与衰老相关的途径和基因快速或加速进化的证据,基于这样的假设,即这些基因和途径更有可能指定人类相对于黑猩猩寿命更长。如材料和方法所述,我们使用了由黑猩猩测序和分析协会(2005年)之前已经证明能够检测有意义的生物学发现。
3.1.1. 纯化选择在衰老相关基因中占主导地位
我们从GenAge数据库的构建12中获得了242个与衰老相关的蛋白编码人类基因,主要是在模型生物体中(de Magalhaes等人,2005年)其中211个基因可以正确分析。虽然尚不清楚这211个基因中控制人类衰老的比例(如果有的话),但我们有理由假设这一列表中包含了丰富的此类基因,这使我们能够估计影响推测的衰老相关基因的进化压力。我们惊讶地发现,GenAge数据集作为一个整体,即使用串联K一和K秒值(见材料和方法)-在统计学上显著降低K一/K秒比人类黑猩猩直系亲属的平均值(0.23)高0.16(第页< 0.001; MWU)。这表明,与偶然发生的情况相比,与衰老相关的基因往往受到更强的功能约束。同样,自从黑猩猩和人类分化以来,GenAge中的功能类别似乎经历了净化选择()我们几乎没有发现快速进化的证据(见下文)。
表1
GenAge类别 | n个 | K一/K秒 | 第页-值高一 | 第页-值低一 | n个四向比较b条 | K一/K秒在原始人中 | K一/K秒在murids | 人类AA差异。c(c) | 平移AA差异。c(c) |
---|
细胞凋亡 | 65 | 0.13 | 1 | 2.8电子-09 | 43 | 0.082 | 0.10 | 14 | 15 |
钙 | 6 | 0.18 | 0.75 | 0.25 | 三 | 0.29 | 0.17 | 2 | 0 |
代谢细胞周期控制 | 42 | 0.13 | 1 | 7.5E−07英寸 | 26 | 0.057 | 0.11 | 10 | 9 |
DNA缩合 | 28 | 0.22 | 0.67 | 0.33 | 8 | 0.32 | 0.16 | 6 | 9 |
DNA修复 | 38 | 0.25 | 0.24 | 0.79 | 16 | 0.17 | 0.17 | 9 | 11 |
DNA | 15 | 0.24 | 0.42 | 0.58 | 7 | 0.079 | 0.14 | 4 | 2 |
复制能量装置 | 19 | 0.12 | 1 | 5.1 E−05 | 12 | 0.13 | 0.082 | 7 | 6 |
增长和发展 | 66 | 0.14 | 1 | 5.9东至10 | 44 | 0.13 | 0.13 | 37 | 17 |
其他 | 47 | 0.17 | 1 | 2.9E−3 | 27 | 0.17 | 0.14 | 12 | 16 |
氧化还原和氧化 | 23 | 0.18 | 0.92 | 0.11 | 13 | 0.10 | 0.071 | 4 | 三 |
调节信号 | 52 | 0.13 | 1 | 1.9E−08年 | 27 | 0.13 | 0.097 | 15 | 16 |
应力响应 | 31 | 0.14 | 1 | 2.4欧-4 | 14 | 0.10 | 0.087 | 4 | 5 |
转录调控 | 39 | 0.15 | 1 | 2.6E−3级 | 22 | 0.062 | 0.084 | 三 | 7 |
功能受限的基因可能与谱系特异性表型变化有关,这可以被检测为谱系特异的加速进化,例如原始人类与鼠类相比(Dorus等人,2004年). 对于使用黑猩猩、人类、小鼠和大鼠直系图计算的所有7043个四向基因比对黑猩猩测序和分析协会(2005年),平均值K一/K秒原始人和鼠的比率分别为0.20和0.13。平均值K一/K秒GenAge中的基因比率在原始人中为0.13,在鼠类中为0.12,这并不表明原始人中与衰老相关的基因加速进化的一般模式。我们也未能在GenAge的功能类别中检测到加速进化(). 钙代谢等类别(n个=3)和DNA缩合(n个=8)有相当高的K一/K秒原始人类谱系中的比率,但鉴于这些类别中的基因数量很少,我们不能排除这些是假阳性。
我们还试图找到GenAge中人类谱系特异性蛋白变化的证据。总的来说,与55只黑猩猩的血统特异性变化相比,有66只人类的血统特异性变化。根据模拟,在一组氨基酸发生如此多变化的基因中,观察到人类谱系特异性变化的比例高于黑猩猩谱系特异性变化的概率为32%,这并不意味着GenAge中人类谱系特异性蛋白变化的比例高于预期。在GenAge类别中,生长和发育蛋白质的异常值最大,其中人类谱系发生的变化是黑猩猩谱系的两倍多(). 这可能意味着血统特异性选择,因为观察到如此高比例的人类血统特异变化的概率为1.2%。最后,我们使用了Clark等人(2003)以确定人类谱系特异性进化是否发生在GenAge或其类别中(参见材料和方法),但我们没有发现人类谱系在任何这些途径上加速进化的证据(未显示)。
由于人类和黑猩猩之间的相似性K一/K秒当关注基因组时,比率更准确。尽管如此,GenAge中单个基因的结果证实了我们的评估,即这些假定的衰老相关基因往往会经历净化选择。除非出现测序错误,而且已经排除了与黑猩猩序列一致的人类序列中少于95%的蛋白质,56种蛋白质在人类和黑猩猩之间是完全保守的。偶然会发现43个这样的蛋白质,即29%的GenAge蛋白质至少95%与黑猩猩的同源序列一致。模拟表明,观察到如此高比例的蛋白质在人类和黑猩猩之间完全保守的概率为8.1%。同样,尽管4.4%的基因具有K一/K我比率高于1,GenAge中211个基因中只有2个(APOC3型和PMCH公司)符合这一标准。根据模拟,观察到如此少数量的基因具有K一/K我1以上的比率为0.45%。除了巴西航空公司1之前的结果(包括连锁不平衡分析)表明,哪些是正选择(Huttley等人,2000年)、和警告(见下文),GenAge中其他具有正向选择潜力的基因包括CDKN2A型(K一/K我=0.90)和EMD公司(K一/K我=0.84和K一/K秒= 1.22).
我们还研究了K一/K秒使用人-恒河猴、人-小鼠和人-大鼠基因对的GenAge比率。尽管生物寿命短,进化性强,生物学上与人类距离较远(见材料和方法),但这些额外的物种为人类-黑猩猩基因对提供了一个补充模型。使用人类恒河猴基因对K一/K秒GenAge中的基因比率为0.16,显著低于所有人类-恒河猴基因对的平均值(0.20)(第页< 0.001; MWU)。然而,来自人-鼠和人-鼠基因对的结果并没有证实这一观点K一/K秒GenAge中的基因比率不低于偶然的预期(未显示)。因此,似乎在灵长类GenAge基因中观察到的功能约束强于平均水平,在啮齿动物谱系中没有观察到。我们对这一观察没有令人满意的解释。
3.1.2. 在页面关联的GO类别中进行选择
当然,可以说GenAge只包含一小部分(如果有的话)与人类衰老有关的基因,因为大多数基因尚待确定。为了进一步研究作用于推测的衰老相关通路的选择力,并探索正选择或快速进化是否在人类和黑猩猩分化后对修复和防御机制产生作用,我们检测了K一/K秒比率,从黑猩猩测序和分析协会(2005年)在GO类别中,可能与老化有关(). 大多数被假设在衰老和身体维持中起作用的途径似乎没有经历快速进化,但也有少数例外。例如,有证据表明DNA修复(GO:0006281;K一/K秒= 0.26;n个=121),在人类和黑猩猩分化后,与DNA修复相关的GO类别可能经历了快速进化(). 然而,尚不清楚这些GO类别中的加速进化是否发生在人类而非黑猩猩谱系中,因为这两个谱系中都发生了谱系特异性蛋白变化(未显示,但见讨论)。
表2
GO(开始)b条 | 姓名 | n个 | K一/K秒 | 第页-值高c(c) | 第页-值低c(c) |
---|
与DNA修复或处理相关的生物过程GO类别 |
0006974 | 对DNA损伤刺激的反应 | 148 | 0.27 | 2007年9月9日 | 1 |
0006281 | DNA修复 | 121 | 0.26 | 1.1E−04 | 1 |
0042770 | DNA损伤反应,信号转导 | 15 | 0.38 | 3.2东风-04 | 1 |
0006306 | DNA甲基化 | 11 | 0.087 | 1 | 1.8东风-04 |
0006471 | 蛋白质氨基酸ADP核糖基化 | 10 | 0.36 | 3.3 E−03 | 1 |
0006298 | 不匹配修复 | 9 | 0.32 | 5.8东风-03 | 1 |
0006282 | DNA修复调控 | 5 | 1.1 | 4.6电子-08 | 1 |
0030330 | p53类介体的DNA损伤反应和信号转导 | 三 | 1.4 | 1.3东经-13 | 1 |
其他生物工艺GO类别 |
0006950 | 对压力的反应 | 594 | 0.27 | 1.1E−17 | 1 |
0006118 | 电子输运 | 216 | 0.23 | 0.039 | 0.96 |
0006511 | 泛素依赖性蛋白质分解代谢 | 82 | 0.20 | 0.65 | 0.38 |
0006512 | 泛素循环 | 70 | 0.10 | 1 | 4.6东风-11 |
0006979 | 氧化应激反应 | 27 | 0.16 | 0.95 | 0.066 |
0007254 | JNK级联 | 23 | 0.12 | 1 | 2.6东经-04 |
0006120 | 线粒体电子传递,NADH到泛醌 | 12 | 0.12 | 0.97 | 0.054 |
0001666 | 缺氧反应 | 5 | 0.043 | 1 | 0.014 |
分子功能GO类别 |
0005179 | 激素活性 | 77 | 0.37 | 1.2E−08年 | 1 |
0004840 | 泛素结合酶活性 | 44 | 0.094 | 1 | 5.1E−06版 |
0003773 | 热休克蛋白活性 | 29 | 0.080 | 1 | 2008年4月4日 |
0003684 | 受损的DNA结合 | 24 | 0.37 | 8.5E−05 | 1 |
0005184 | 神经肽激素活性 | 15 | 0.39 | 4.8E−03型 | 1 |
0004364 | 谷胱甘肽转移酶活性 | 14 | 0.50 | 4.5东风-04 | 1 |
0004407 | 组蛋白脱乙酰酶活性 | 8 | 0.14 | 0.95 | 0.082 |
我们没有发现GO类别与氧化应激相关的证据,例如氧化应激反应(GO:0006979;K一/K秒= 0.16;n个=27),经历了快速演变(). 应激反应类别(GO:0006950)可能经历了快速进化,但这是一种误导,因为它的子类别,即DNA修复和对害虫/病原体/寄生虫的反应(GO:0009613)似乎经历了快速演变。当消除这些可能代表偏见来源的子类别的影响时,快速选择不再考虑压力反应(未显示)。同样,有证据表明进化速度很快(K一/K秒= 0.50;n个=14)谷胱甘肽转移酶活性(GO:0004364)。因为谷胱甘肽转移酶通常与衰老有关,例如GSTA4公司和GSTP1标准,没有发现正选择,并且由于来自该GO类别的一些蛋白质在应对通常与快速进化相关的环境应激源中发挥作用,这些发现不一定与衰老有关,即使很有趣。同样,激素活性(GO:0005179)似乎经历了快速进化(K一/K秒=0.37,带n个=77),但激素的变化可能与衰老以外的过程中的适应有关。事实上,来自GO类别的基因K一/K秒>1与不同的功能相关,包括通常与快速进化相关的过程,如免疫反应、生殖和发育,因此尚不清楚GO类别中的快速进化是否与衰老相关。
总的来说,之前发现人类和黑猩猩之间高度保守的与衰老相关的途径和机制包括拓扑异构酶,如TOP3B,JNK途径,包括MAPK8和MAPK9,IGF1,其受体IGF1R和IGFBP3,Sir2同源物SIRT1,TOR同源物FRAP1,TP53,热休克蛋白(GO:0003773),如HSPA1B,HSPA8,HSPA9B,HSP90AA1和HSPD1、组蛋白去乙酰化酶(GO:0004407),如HDAC1和HDAC3,泛素循环(GO:0006512),如SUMO、UBB、UBE2I和UCHL1。完整结果可在线获取(http://genomics.sensense.info/evolution/chimp.html).
3.2. 不同模型中与衰老相关的基因受到不同的选择力
如前所述,GenAge以不同模型系统中与衰老相关的基因为特征,包括来源于从酵母到啮齿动物、细胞模型和前体综合征等生物体的基因,尽管可能存在一些重叠。因此,我们研究了在不同的模型系统中,不同的进化压力是否作用于与衰老相关的基因(). 有趣的是,在非人类哺乳动物中,与衰老直接相关的基因(n个=21)和人类(n个=3)平均值最高K一/K秒比率(分别为0.33和0.47),而非哺乳动物模型(通常是无脊椎动物)中与衰老相关的基因的平均值最低K一/K秒比率(0.083)。较高的平均值K一/K秒非人类哺乳动物(通常是啮齿动物)中与衰老相关的基因比例可能是由一些基因造成的,例如警告和巴西航空公司1因为如果排除这两个基因的平均值K一/K秒这一类别的比率降至0.19。同样在非人类哺乳动物中,与过早衰老相关的基因也稍高K一/K秒比率(0.34;n个=16)比与寿命延长相关的基因(0.27;n个=6),但这似乎是由于巴西航空公司1和警告在前者中。与细胞衰老模型相关的基因也非常低K一/K秒比率(0.11)。最后,我们发现,与鼠相比,原始人类没有明显的加速进化迹象,除了哺乳动物中与衰老有关的基因,这些基因可能再次受到一些基因的影响,例如巴西航空公司1和警告.
表3
通过选择过程对GenAge类别进行操作的选择性部队。
GenAge类别一 | n个 | K一/K秒 | 第页-值高b条 | 第页-值低b条 | n个四向比较c(c) | K一/K秒原始人类 | K一/K秒在murids | 人类AA差异。d日 | 平移AA差异。d日 |
---|
人类 | 三 | 0.472 | 0.017 | 0.99 | 2 | 0.46 | 0.22 | 4 | 8 |
哺乳动物 | 21 | 0.331 | 0.0026 | 1 | 9 | 0.32 | 0.20 | 9 | 10 |
模型 | 22 | 0.083 | 1 | 4.6东至07 | 14 | 0.069 | 0.10 | 2 | 三 |
细胞 | 14 | 0.109 | 1 | 5.7E−4 | 6 | 0.068 | 0.19 | 2 | 1 |
功能性的 | 55 | 0.184 | 0.99 | 0.017 | 22 | 0.094 | 0.11 | 8 | 4 |
上游 | 25 | 0.138 | 1 | 1.1E−3 | 15 | 0.18 | 0.21 | 9 | 2 |
下游 | 46 | 0.097 | 1 | 7.6东至16 | 28 | 0.085 | 0.080 | 22 | 17 |
假定的 | 56 | 0.159 | 1 | 2.3东至4 | 38 | 0.15 | 0.13 | 16 | 19 |
我们还查看了K一/K秒基因比率(n个=14)与人类寿命或老年人生存有关,因此可能与健康老龄化有关。平均值K一/K秒比率为0.16,没有明显证据表明高于或低于平均值K一/K秒比率。这些基因都没有K一/K秒比率>1。最强的异常值是SIRT3型用一个K一/K秒=0.85(但K一/K我=0.42)和APOC3型用一个K一/K我=1.19(但K一/K秒= 0.35). 完整的结果可在线获取(http://genomics.sensecence.info/evolution/hypemp.html).
根据使用人类恒河猴、人鼠和人鼠基因对的选择聚类对GenAge数据集进行的分析与使用人类黑猩猩基因对获得的上述结果基本一致。例如,低等生物和细胞模型中与衰老相关的基因通常具有较低的K一/K秒哺乳动物中与衰老相关的基因的比例。
3.3. 节段性进展综合征相关基因的进化
如上所述,我们在GenAge基因中发现了积极选择或快速或加速进化的有限证据。除了巴西航空公司1目前尚不清楚自人类和黑猩猩分化以来,GenAge中是否有任何基因经历了正向选择。然而,我们想评估与人类早衰相关的基因是否也与人类谱系中的衰老进化有关。
突变导致与人类早衰最相似的三种基因是警告,ERCC8号机组、和LMNA公司,其中警告沃纳综合征基因是最引人注目的病例。其他可能被认为是节段性进展综合征的基因包括英国标准CL2,ERCC6号机组,NBN公司、和AGPAT2公司(马丁,1978年&2005;de Magalhaes,2005年). 我们未能在这些基因中的任何一种人或人中找到积极选择或快速进化的重要证据,除了警告.
在以下情况下警告,两个K一/K秒(0.91)和K一/K我(0.59)比率表明,自从人类和黑猩猩分化以来,该基因可能经历了快速进化。警告也有更高的K一/K秒原始人类的比率(0.75)比鼠类的比率(0.31)高,这可能表明原始人类进化加速,平均值更高K一/K秒人-恒河猴基因对(0.53)、人-小鼠基因对(0.32)和人-大鼠基因对(0.29)的比率。尽管这些结果可能是误报(见讨论),但它们表明警告人类谱系中的基因。以前,Clark等人(2003)报告说警告可能已经经历了人类血统的加速进化(第页= 0.028). 然而,通过人-黑猩猩-人蛋白质比对,我们发现了11种黑猩猩特异性蛋白质变化,但WRN中只有4种人类特异性蛋白质改变,加上2种不能归因于任何血统的变化。然而,由于黑猩猩基因组序列仍被认为是一个草案,其中的错误可能会人为地增加黑猩猩谱系特异性蛋白质变化的数量。
我们的结果和K一/K秒比率已纳入GenAge数据库(http://genomics.sensecence.info/genes/human.html)这可能有助于研究衰老的研究人员研究灵长类动物衰老相关基因的进化。他们可以在网上进一步获得(http://genomics.sensense.info/evolution/chimp.html).
4.讨论
为了理解导致人类衰老过程显著特征的遗传变化,需要多个灵长类基因组,因此我们的工作应被视为初步的或探索性的。我们用于识别导致物种寿命差异的基因的方法是新颖且有前途的,但其检测极限仍然局限于相对较强的信号(). 此外,到目前为止,这种大规模的全基因组分析仅在蛋白质编码区可行,但基因调控的差异可能也指明了物种在衰老过程中的差异,并且调控区、RNA基因甚至一些未知基因产物可能发生重要变化,而我们的方法无法检测到这些变化。尽管存在这些局限性,但我们的工作为模型系统中作用于与衰老相关基因的进化力提供了新的见解。
使用我们的测试统计量作为序列长度和K一/K秒比率。黑色直线代表K一/K秒一个给定的蛋白质编码序列将被检测为正在经历快速进化的比率阈值。序列长度可以是单个基因或一组基因的串联序列(人类黑猩猩的平均基因对,仅包括蛋白质编码序列,约为1.4 kb),以对数标度绘制。灰色虚线为平均值K一/K秒所有人类黑猩猩基因对的比率(0.23)。K一/K秒比率假定为常数K秒同义和非同义位点相对于序列长度的比例是恒定的,这两个位点都是根据所有人类黑猩猩基因对的平均值进行估计的。
虽然我们在一些与衰老相关的基因和途径中发现了选择的证据,但我们强烈建议谨慎解释这些结果。例如,GO类别中的基因可能与与衰老无关的功能有关,因此任何快速进化的证据都可能反映出对另一个性状的选择(见下文)。一个潜在的有趣发现是DNA修复蛋白的假定快速进化。其他研究报告了DNA代谢基因中正选择的证据,包括ERCC8号机组,在最近的人类进化中(Wang等人,2006年). 这些结果可能与预测DNA修复/代谢基因在物种衰老差异中起作用的衰老假说有关(参见de Magalhaes,2005年;Martin,2005年). 另一方面,DNA修复基因的快速进化并不是长寿人类或原始人独有的,事实上,一些结果表明,与原始人相比,鼠类DNA修复基因进化速度更快(黑猩猩测序和分析协会,2005年). 更详细的研究包括几个哺乳动物基因组,有必要研究DNA修复/代谢基因的变化是否与衰老的进化有关。
总的来说,我们使用GenAge数据集的结果表明,至少在灵长类谱系中,模型系统中与衰老相关的基因往往处于净化选择和强于平均水平的功能约束之下。尽管在强烈功能约束下的基因可能在谱系特异性性状的进化中发挥作用,但我们发现,在原始人类或人类谱系中,几乎没有证据表明与衰老相关的基因加速进化。具有这种选择模式的基因往往具有哺乳动物之间保守的关键生物功能,这些进化力量通常保留现有的表型。然而,因为我们的工作是第一次探索作用于衰老相关基因和途径的进化压力,不同的解释可以被概念化,我们的结果并不能证明这些基因可能与人类衰老有关。当然,我们的分析可能无法检测到一个或几个基因的表型重要变化,甚至可能是调控区域的变化。例如,尽管我们没有发现线粒体电子传递链(GO:0006118和GO:0006120)快速进化的证据,但在类人猿灵长类动物中已经报道了其中一些蛋白质的快速进化(Grossman等人,2004年),并且有人假设这代表了一种延长生命的适应(de Magalhaes,2005年b). 尽管是推测性的,但另一种假设是,与衰老有关的基因需要极其精确,因为它们可能具有其他功能,而且很微妙,迄今为止无法检测到导致人类寿命延长的变化。
除了GO类别的结果(可能会有不同的解释)外,我们发现几乎没有证据表明多种维护和维修机制得到了广泛优化,尽管这可能是由于我们分析中的局限性造成的。我们发现,在不同模型系统中,作用于与衰老相关基因的选择压力存在差异,例如观察到非哺乳动物模型生物中与衰老相关的基因比哺乳动物中与衰老有关的基因受到更强的功能约束,这是一种新颖的现象。这些结果可能反映出,在人类-黑猩猩基因对中,在非哺乳动物模型生物中,同源基因参与衰老的基因的整个进化过程中,具有较高的保守性。事实上,据报道,在模型生物中,人类AGIG相关基因同源物的比例高于预期(Budovsky等人,2007年). AGIG相关基因的高度保守性可能与其报道的高平均相互作用有关(Promislow,2004年;Budovsky等人,2007年). 我们的研究结果虽然只是探索性的,但似乎并不支持这样的假设,即大多数基因和机制都是模型生物研究的重点,它们很可能调节物种在衰老过程中的差异,至少在灵长类中是如此。即使迄今为止在模型系统中确定的衰老机制不能解释原始人之间衰老的差异,但这些机制在生物体之间可能是保守的,例如,作为将生殖和食物供应与衰老联系起来的保守机制的一部分(鹧鸪与宝石,2002). 换句话说,这些基因的生物学功能可能是与衰老的多效性效应相关的保守途径的一部分,这些多效性影响衰老的种内差异,但不调节衰老的种间差异。
正如引言中所提到的,我们在这项工作中的理论基础是现代进化论支持的概念(例如。,查诺夫,1993;Kirkwood和Austad,2000年)人类寿命延长和衰老延迟是直接选择压力的产物。然而,由于我们对灵长类动物衰老的遗传学知之甚少,我们不能排除其他选择。也许人类延缓衰老是选择作用于其他特征(如发育、生长或繁殖)的副产品。例如,一个假设是,在人类和黑猩猩分化后,自然选择可能主要作用于发展途径,以延长发育、生长和寿命(de Magalhaes and Church,2005年). 事实上,一种作用于发育、生长和衰老的机制是内分泌系统(参见麦金尼和麦克纳马拉,1991年;鹧鸪与宝石,2002;de Magalhaes,2005年;de Magalhaes and Church,2005年). 因此,激素的适应,甚至那些与发育和生殖有关的激素,可能会对衰老产生一些影响。与生长和发育相关的GenAge蛋白中人类特异性蛋白变化的比例较高,以及激素活性GO类别中推测的快速进化支持了这一观点。如果人类衰老延迟开始是人类发育延迟的副产品,那么模型系统中与衰老相关的基因可能会受到净化选择。
最后,我们的结果表明警告,一种突变导致表型类似于人类患者加速衰老的基因。我们的结果并不能证明警告经历了适应性的变化,诚然可能是偶然的结果,因为警告编码一种具有1432个氨基酸的相对较大的蛋白质K一/K秒与编码较小蛋白质的基因相比,这些比率不太可能被夸大。此外,一个有趣的巧合是,与最令人印象深刻的人类孕激素综合征相关的基因是GenAge中最高的基因之一K一/K秒比率。人们很容易猜测警告也参与了灵长类和原始人类的衰老进化。由于人类和黑猩猩的分化,WRN的大多数蛋白质变化都是黑猩猩谱系特异性的,这使得WRN的蛋白质变化不太可能对人类寿命的进化做出贡献。这可能意味着警告与在MCPH1型在谱系的早期部分选择最明显的基因(Vallender和Lahn,2004年). 灵长类和哺乳动物进化过程中WRN蛋白变化的分析值得进一步关注。
致谢
作者感谢David Reich和Joana Costa对之前草案的评论,感谢Steve Elledge和教会实验室成员的宝贵讨论。J.P.de Magalháes由NIH-NHGRI CEGS向G.M.Church拨款支持。
脚注
出版商免责声明:这是一份未经编辑的手稿的PDF文件,已被接受出版。作为对客户的服务,我们正在提供这份早期版本的手稿。手稿在以最终可引用的形式出版之前,将进行编辑、排版和校对。请注意,在制作过程中可能会发现可能影响内容的错误,适用于该期刊的所有法律免责声明都适用。
工具书类
- Ashburner M等人,《基因本体论:生物学统一的工具》。基因本体联盟。自然遗传学。2000;25:25–29. [PMC免费文章][公共医学][谷歌学者]
- Ayyadevara S、Engle MR、Singh SP、Dandapat A、Lichti CF、Benes H、Shmookler Reis RJ、Liebau E、Zimniak P.秀丽隐杆线虫的寿命和抗逆性是通过表达能够代谢脂质过氧化产物4-羟基壬醛的谷胱甘肽转移酶而提高的。老化细胞。2005;4:257–271.[公共医学][谷歌学者]
- Birney E等人,2006年合奏。核酸研究。2006;34:D556–561。 [PMC免费文章][公共医学][谷歌学者]
- Brown-Borg HM、Borg KE、Meliska CJ、Bartke A.矮小小鼠和衰老过程。自然。1996;384:33.[公共医学][谷歌学者]
- Budovsky A、Abramovich A、Cohen R、Chalifa-Caspi V、Fraifeld V。长寿网络:构建和含义。机械老化发展。2007;128:117–124.[公共医学][谷歌学者]
- Burkle A、Beneke S、Muiras ML.聚(ADP-核糖基)化和老化。Exp Gerontal公司。2004;39:1599–1601.[公共医学][谷歌学者]
- Cadenas E,Davies KJ。线粒体自由基生成、氧化应激和衰老。自由基生物医药。2000;29:222–230.[公共医学][谷歌学者]
- 查诺夫EL。生命史不变量:进化生态学中对称性的一些探索。牛津大学出版社;牛津:1993年。[谷歌学者]
- Clark AG、Glanowski S、Nielsen R、Thomas PD、Kejariwal A、Todd MA、Tanenbaum DM、Civello D、Lu F、Murphy B、Ferriera S、Wang G、Zheng X、White TJ、Sninsky JJ、Adams MD、Cargill M。从人-黑猩猩-鼠同源基因三组推断非中性进化。科学。2003;302:1960–1963.[公共医学][谷歌学者]
- de Magalhaes JP.哺乳动物衰老是由基因控制的吗?生物老年学。2003;4:119–120.[公共医学][谷歌学者]
- de Magalhaes JP。开放的怀疑论:从基因扰动推断人类衰老的因果机制。老化研究版本。2005年a;4:1–22.[公共医学][谷歌学者]
- de Magalhaes JP。非人类灵长类动物中与人类疾病相关的线粒体突变。分子进化杂志。2005年b;61:491–497.[公共医学][谷歌学者]
- de Magalhaes JP,Church GM。基因组优化生殖:发育计划导致的衰老。生理学(贝塞斯达)2005;20:252–259.[公共医学][谷歌学者]
- de Magalhaes JP、Costa J、Toussant O.HAGR:人类老龄化基因组资源。核酸研究。2005;33(数据库问题):D537–543。 [PMC免费文章][公共医学][谷歌学者]
- de Magalhaes JP。衰老和抗衰老研究比较研究中的物种选择。收件人:Conn PM,编辑。人类老龄化模型手册。爱思唯尔学术出版社;马萨诸塞州伯灵顿:2006年。第9-20页。[谷歌学者]
- Dorus S、Vallender EJ、Evans PD、Anderson JR、Gilbert SL、Mahowald M、Wyckoff GJ、Malcom CM、Lahn BT。智人起源中神经系统基因的加速进化。单元格。2004;119:1027–1040.[公共医学][谷歌学者]
- Erwin JM、Hof PR、Ely JJ、Perl DP。一个老年学:通过对类人猿和其他灵长类动物的研究,提高对衰老的理解。作者:Erwin JM,Hof PR,编辑。非人类灵长类的衰老。卡格尔;巴塞尔:2002年。第1-21页。[谷歌学者]
- 芬奇CE。长寿、衰老和基因组。芝加哥大学出版社;芝加哥和伦敦:1990年。[谷歌学者]
- 芬奇CE,斯坦福CB。肉类适应性基因和人类缓慢衰老的进化。Q生物版。2004;79:3–50.[公共医学][谷歌学者]
- 新泽西州霍尔布鲁克Finkel T。氧化剂、氧化应激和衰老生物学。自然。2000;408:239–247.[公共医学][谷歌学者]
- Flurkey K、Papaconstantinou J、Miller RA、Harrison DE。生长激素生成缺陷突变小鼠的寿命延长和免疫和胶原蛋白老化延迟。美国国家科学院院刊。2001;98:6736–6741. [PMC免费文章][公共医学][谷歌学者]
- Fraga CG,Shigenaga MK,Park JW,Degan P,Ames BN。衰老过程中对DNA的氧化损伤:大鼠器官DNA和尿液中的8-羟基-2′-脱氧鸟苷。美国国家科学院院刊。1990;87:4533–4537. [PMC免费文章][公共医学][谷歌学者]
- Gray DA、Tsirigotis M、Woulfe J.泛素、蛋白酶体和衰老的大脑。科学老化知识环境。2003年:RE6。2003. [公共医学][谷歌学者]
- Grossman LI、Wildman DE、Schmidt TR、Goodman M.加速了类人灵长类动物电子传输链的进化。趋势Genet。2004;20:578–585.[公共医学][谷歌学者]
- Hasty P、Campisi J、Hoeijmakers J、van Steeg H、Vijg J。衰老和基因组维护:小鼠的经验教训?科学。2003;299:1355–1359.[公共医学][谷歌学者]
- Hill K、Boesch C、Goodall J、Pusey A、Williams J、Wrangham R.野生黑猩猩死亡率。J Hum进化。2001;40:437–450.[公共医学][谷歌学者]
- Huttley GA、Easteal S、Southey MC、Tesoriero A、Giles GG、McCredie MR、Hopper JL、Venter DJ。人类和黑猩猩体内肿瘤抑制因子BRCA1的适应性进化。澳大利亚乳腺癌家族研究。自然遗传学。2000;25:410–413.[公共医学][谷歌学者]
- Katschinski DM。缺氧和衰老之间有分子联系吗?Exp Gerontal公司。2006;41:482–484.[公共医学][谷歌学者]
- 柯克伍德TB。了解衰老的奇怪科学。单元格。2005;120:437–447.[公共医学][谷歌学者]
- Kirkwood TB,Austad SN.我们为什么变老?自然。2000;408:233–238.[公共医学][谷歌学者]
- Lin YJ,Seroude L,Benzer S.果蝇突变株methuselah的寿命延长和抗逆性。科学。1998;282:943–946.[公共医学][谷歌学者]
- 伦巴第DB、Chua KF、Mostoslavsky R、Franco S、Gostissa M、Alt-FW。DNA修复、基因组稳定性和衰老。单元格。2005;120:497–512.[公共医学][谷歌学者]
- Martin GM。与衰老病理生物学潜在相关的男性遗传综合征。出生缺陷原始序列。1978;14:5–39.[公共医学][谷歌学者]
- Martin GM。智人衰老表型的遗传调控。单元格。2005;120:523–532.[公共医学][谷歌学者]
- Martin GM、Austad SN、Johnson TE。衰老的遗传分析:氧化损伤和环境应激的作用。自然遗传学。1996;13:25–34.[公共医学][谷歌学者]
- McKinney ML,McNamara KJ公司。异时性:个体发育的进化。增压室;纽约:1991年。[谷歌学者]
- Migliaccio E、Giorgio M、Mele S、Pelicci G、Reboldi P、Pandolfi PP、Lanfrancone L、Pelicsi PG。p66shc衔接蛋白控制哺乳动物的氧化应激反应和寿命。自然。1999;402:309–313.[公共医学][谷歌学者]
- 黑猩猩测序分析联盟。黑猩猩基因组的初始序列以及与人类基因组的比较。自然。2005;437:69–87.[公共医学][谷歌学者]
- Miller RA公司。克莱梅尔奖讲座:衰老有基因吗?老年医学杂志。1999;54:B297–307。[公共医学][谷歌学者]
- Morbeck ME,Galloway A,Sumner DR。在贡贝变老:野生黑猩猩的骨骼老化。作者:Erwin JM,Hof PR,编辑。非人类灵长类的衰老。卡格尔;巴塞尔:2002年。第48-62页。[谷歌学者]
- Moskovitz J、Bar-Noy S、Williams WM、Requena J、Berlett BS、Stadtman ER。蛋氨酸亚砜还原酶(MsrA)是哺乳动物抗氧化防御和寿命的调节器。美国国家科学院院刊。2001;98:12920–12925. [PMC免费文章][公共医学][谷歌学者]
- Nielsen R、Bustamante C、Clark AG、Glanowski S、Sackton TB、Hubisz MJ、Fledel-Alon A、Tanenbaum DM、Civello D、White TJ、J JS、Adams MD、Cargill M.对人类和黑猩猩基因组中的阳性选择基因进行扫描。《公共科学图书馆·生物》。2005;三:e170。 [PMC免费文章][公共医学][谷歌学者]
- Partridge L,Gems D.衰老机制:公共还是私人?自然资源部Genet。2002;三:165–175.[公共医学][谷歌学者]
- Partridge L,Gems D.超越衰老进化理论,从功能基因组学到进化基因组学。经济发展趋势。2006;21:334–340.[公共医学][谷歌学者]
- 蛋白质网络、多效性和衰老的进化。过程生物科学。2004;271:1225–1234. [PMC免费文章][公共医学][谷歌学者]
- Richardson B.衰老对DNA甲基化的影响。老化研究版本。2003;2:245–261.[公共医学][谷歌学者]
- Rogina B,Helfand SL,Frankel S.果蝇Rpd3脱乙酰酶和热量限制对长寿的调节。科学。2002;298:1745.[公共医学][谷歌学者]
- Rose先生。衰老进化生物学。牛津大学出版社;纽约:1991年。[谷歌学者]
- Sohal RS、Weindruch R.氧化应激、热量限制和衰老。科学。1996;273:59–63. [PMC免费文章][公共医学][谷歌学者]
- Tatar M,Bartke A,Antebi A.胰岛素样信号对衰老的内分泌调节。科学。2003;299:1346–1351.[公共医学][谷歌学者]
- Vallender EJ,Lahn BT.人类基因组的正选择。人类分子遗传学。2004;13(规范2):R245–254。[公共医学][谷歌学者]
- Varki A,Altheide TK。人类和黑猩猩基因组的比较:在干草堆中寻找针头。基因组研究。2005;15:1746–1758.[公共医学][谷歌学者]
- Wang ET,Kodama G,Baldi P,Moyzis RK。最近推断的智人达尔文选择的全球景观。美国国家科学院院刊。2006;103:135–140. [PMC免费文章][公共医学][谷歌学者]
- Wang MC、Bohmann D、Jasper H.JNK信号传递赋予果蝇对氧化应激的耐受性并延长寿命。开发单元。2003;5:811–816.[公共医学][谷歌学者]
- 渡边H等。黑猩猩22号染色体的DNA序列和比较分析。自然。2004;429:382–388.[公共医学][谷歌学者]
- 横山K、福本克K、村上春树T、原田S、细野R、瓦德瓦R、三井Y、大久马S。通过敲入额外的hsp70F(mot-2(mothlin)/mthsp70/Grp75的同源物)副本,延长秀丽隐杆线虫的寿命。FEBS信函。2002;516:53–57.[公共医学][谷歌学者]