狗的衰老：分子生物学

机器翻译2026年04月09日约 28 分钟

宠物狗衰老过程中的血液生物标志物与品种遗传学

Vista Sohrab¹,²,³*, Michelle E. White³, Benjamin R. Harrison⁴, Rob Bierman⁵, Abbey Marye⁶, Kathleen Morrill Pirovich¹,²,³,⁷, Diane P. Genereux³, Kate Megquier³, Xue Li¹,²,³,⁸, Brittney Kenney¹,³, Cindy Reichel⁹, Dog Aging Project Consortium‡, Noah Snyder-Mackler¹⁰, Joshua M. Akey⁵, Daniel E. L. Promislow¹¹, Frances L. Chen¹,³†, Elinor K. Karlsson¹,³,¹²*†

摘要

宠物犬与人类共享相似的生活环境，且衰老速度更快，是衰老研究中极具价值的转化模型。本研究利用“犬类衰老项目”（Dog Aging Project）中7627 只犬的基因组和表型数据（其中976 只犬完成了159 种血液代谢物和临床分析物的检测），构建了首个犬类全基因组关联研究（GWAS）目录。犬类与人类的血液性状定位到同源基因位点，表明相关通路具有高度保守性。品种血统可解释血液性状的显著变异，且对被毛类型等外观特征的选择具有多效性代谢效应。借助混种犬的嵌合血统特征和纵向死亡率数据，本研究识别出短寿品种中升高且可预测个体死亡风险的血液性状（包括球蛋白和钾），以及长寿品种中富集的保护性性状（如乙醇胺）。尽管部分与衰老相关的性状与生长激素通路有关，但多数性状独立于该通路，表明犬类衰老具有多因素特性。这些发现确立了犬类作为转化系统的地位，为识别与延长健康寿命相关的衰老遗传决定因素和生物标志物提供了重要依据。

正文

衰老的特征包括生理和认知功能的渐进性衰退、广泛的代谢重塑，以及对慢性疾病的易感性增加，给个体、家庭和医疗系统带来沉重负担[1-5]。这些变化的速度和轨迹在个体间存在巨大差异，反映了遗传、环境和生活方式因素的复杂相互作用[6,7]。包括全基因组关联研究（GWAS）在内的大规模遗传研究已识别出数千个与极端长寿和年龄相关疾病相关的基因位点[8-12]。然而，将这些位点与功能衰退联系起来的生物学机制仍知之甚少，限制了其向有效治疗策略的转化。

与仅针对长寿的研究相比，血液性状的遗传研究为理解衰老提供了互补且更具机制性的方法。许多循环血液性状与衰老相关通路相关。血清化学图谱和全血细胞计数等临床分析物既受遗传变异影响[13-16]，也会随年龄相关的生理变化而改变[17-20]，并可作为疾病风险的生物标志物，反映与病理相关的过程[21-23]。血浆代谢组包含数千种反映代谢活动的小分子，可预测犬类和人类的死亡率、衰弱程度及功能衰退[24-28]。识别血液性状水平和纵向轨迹个体差异背后的遗传变异，有望为衰老相关疾病开发新的干预手段—— 例如低密度脂蛋白胆固醇（LDL-C）与他汀类药物治疗心血管疾病的研发[29]。

在人类研究中，识别影响衰老的可修饰代谢或环境因素受限于人类较长的寿命[30]。宠物犬（又称伴侣犬）具有许多与人类衰老相似的特征：其衰老轨迹受遗传、环境和生活方式因素影响而具有异质性[31-33]，但衰老速度更快，寿命约为9-15 年[34-36]。作为家庭动物，犬类与主人经历相似的饮食、社会和环境暴露，且通常接受相当水平的医疗护理[37,38]。这种紧密的人犬关系使得收集高分辨率纵向数据成为可能，包括兽医临床记录、血浆代谢组学分析以及主人报告的生理和认知功能评估[39-43]。犬类与人类中，与衰老相关的指标（如血浆代谢物死亡率风险比）具有强相关性，证实了该系统的转化价值[44]。

与人类相似，犬类的血液性状随年龄、性别和体型而变化，但同时也受品种影响[43,45-47]。现代犬品种是过去175 年内通过严重的种群瓶颈和针对物理及生理性状的持续选择性育种形成的近交种群[48-51]。在美国，大多数犬为纯种犬（single-breed）或具有多个品种血统的混种犬（mixed-breed）[52]。这种特殊的种群历史改变了等位基因频率谱，使得在规模小于人类研究的队列中也能开展高效力的GWAS [53,54]。

纯种犬的寿命通常短于混种犬[34]，部分大体型品种的寿命尤其短暂。品种间成年体型的变异与生长激素/ 胰岛素样生长因子1（GH/IGF-1）通路中的基因密切相关[55-57]。在啮齿动物模型中，该通路的调控会对衰老产生深远影响[58-61]，这引发了一种推测：对GH/IGF-1 变异的选择可能导致部分犬品种寿命缩短[62]。然而，仅体型无法完全解释品种间的寿命变异[63,64]。美国大多数混种犬具有四个或更多品种的血统，体型跨度广泛[52]，这形成了一个自然实验，可将体型的影响与其他品种相关遗传因素对衰老的作用区分开。

本研究利用“犬类衰老项目”（Dog Aging Project）中犬类独特的品种血统、环境多样性和深度纵向表型数据，识别与寿命变异相关的血液性状，并评估其与人类衰老的相关性。

犬类衰老项目

“犬类衰老项目”是规模最大的研究犬类衰老相关生物学、遗传和环境因素的项目[65]。截至2025年12月，该项目已招募52705只犬。本研究分析了2023年底前招募的47444只犬中具有遗传数据的7627只犬的子集。主人通过唾液拭子收集DNA样本，进行测序（平均测序深度0.79x±0.41x；范围0.14x至5.36x），并与CanFam4.0参考基因组比对[66]。基因型的估算来自包含1929只深度测序（12.95x至31.06x）犬类和相关犬科动物[67]。最终数据集包含29089701个双等位基因单核苷酸多态性（SNP），其中9871985个为常见SNP（次要等位基因频率≥0.01）。此前研究已表明，低覆盖度测序结合填充技术可达到与芯片相当的基因型准确性，同时产生10-50倍更高的标记密度，足以捕获大多数常见变异[52]。

遗传决定的品种血统

测序犬的入组年龄从数周至21岁不等（均值±标准差，4.6±3.5岁），比整个入组队列更年轻（6.8±4.3岁）（图1A、B；图S1A）。约一半为雌性（48.8%），82%已绝育，与美国宠物犬种群特征相似（图S1B）[68]。我们通过监督混合分析（supervised admixture analysis），以包含全球公认的360个品种中筛选出的109个常见品种、乡村犬和狼组成的参照组为依据，估算品种血统[69]。根据主人报告的信息将犬分为纯种或混种——此前研究已证实该分类方式可靠[52,69]，且可将参照组中未包含的品种纳入数据集。

图1 研究种群的人口统计学和遗传结构

（A）具有血液性状数据的犬比整个入组种群更年轻；（B）体重无差异。（C）近一半的犬（45%）具有可检测到的多个品种血统（>5%），包括简单混种（任一品种血统占比> 40% 的犬）和更复杂的混种。（D）按血统复杂性分层的体重分布显示，纯种犬和所有混种犬的体重呈双峰分布，而复杂混种犬不呈双峰分布（采用哈蒂根双峰检验验证单峰性，100,000 次bootstrap 重复）。（E）在纯种犬中，体重与品种标准体型相关。（F）在复杂混种犬中，体型与品种血统的关系不太明确—— 大型品种血统可能出现在小型犬中，反之亦然，且极少出现体型极小型或极大型的犬。

混种犬的品种血统比例与纯种犬大致相似。总体而言，45.7% 的测序犬为混种犬，其中一半（51%）为“复杂混种”（任一品种的血统占比不超过40%）（图1C）。纯种犬中常见的品种在混种犬的血统中占比也更高（斯皮尔曼相关系数ρ=0.70）。在最常见的品种中，金毛寻回犬（golden retrievers）和拉布拉多寻回犬（labrador retrievers）的血统在纯种犬中更普遍，而美国比特斗牛梗（American pit pull terriers,）、贵宾犬（poodles）、玩具贵宾犬（toy poodles）和吉娃娃（Chihuahuas）的血统在混种犬中更常见（图S2）。

品种相关的种群结构影响成年体重分布，导致其偏离随机交配种群中预期的近似正态分布[70]。在纯种犬中，极小型和极大型品种的受欢迎程度使得体重分布呈双峰分布（n=4163；哈蒂根双峰检验，p<1×10⁻⁵）（图1D、E）。相比之下，复杂混种犬的体重分布未否定单峰性（n=1760；双峰检验p=0.089），且混种犬极少出现极端体型（图1D、F）。品种血统对混种犬体重变异的解释度（R²=0.75）低于纯种犬（R²=0.84），在复杂混种犬中更低（R²=0.68）；所有关联均具有统计学显著性（置换检验p≤1×10⁻⁵；表S1）。这一模式支持我们的理论框架：由于纯种犬中体型与遗传血统的紧密关联[71] 在混种犬中被打破，混种犬为区分体型与其他影响衰老的品种相关因素提供了机会。

表型

所有主人在入组时完成了关于犬的环境、生活方式、健康状况和行为的基线调查（959 个问题）。我们从基线调查中筛选出154个问题进行遗传分析：这些问题均具有标准化应答类型、个体间存在有意义的变异、应答数量充足且可能具有遗传性（数据S1；图S3）。其中包括26个关于身体活动的问题、40个关于行为和认知健康的问题，以及29个关于环境、病史、饮食和体重的问题[65,72]。我们还纳入了项目开展一年后新增的两项调查（应答数量较少）：犬类肥胖风险与食欲（DORA）调查（34个问题，4895只犬）[73]和莫纳什犬主关系量表（Monash Dog Owner Relationship Scale（MDORS；25个问题，4915只犬）[74]（表S2）。

我们收集了部分犬的血液分子表型（血液性状）（图S1）。具有血液性状数据的犬（3.8±3.1标准差岁）比整个测序队列更年轻（4.6±3.5标准差岁），但其他特征相似（图1A、B）。排除分布受罕见极端值主导、不适合GWAS的性状后，我们分析了159种血液性状：36种临床分析物（最多971只犬检测）（全血细胞计数[CBC]和血清化学图谱[SCP]）和123种血浆代谢物（937只犬检测）（表S2）。

血液性状的遗传力高于调查性状

血液性状的SNP水平遗传力（h²_SNP）高于调查衍生性状（中位h²_SNP=0.39 vs. 0.15；置换检验p<1×10⁻⁵）（图2A）。我们采用限制最大似然法（校正连锁不平衡）估算h²_SNP [75]，并调整年龄、体重、性别（绝育状态纳入性别模型），对于血液性状还调整了空腹时间。与宠物犬身高的高遗传力一致[52]，体重的遗传力极高（h²_SNP=0.81±0.02；N=7378）。两种血液性状的遗传力更高：血清肌酐（h²_SNP=0.88±0.13；N=953）和同型半胱氨酸（h²_SNP=0.87±0.12；N=924）（图2A）[76,77]。总计29种血液性状的遗传力估计值超过50%，包括9种临床分析物和20种血浆代谢物（图S4）。

除体重外，所有调查衍生性状的遗传力均未超过50%（图2A）。这一模式可能反映了这些性状受非遗传因素的强烈影响。遗传力最高的调查性状涉及进食行为、体型相关行为以及影响饲养便利性的行为（如恐惧程度和顺从性）（图S4、数据S2）。相比之下，与犬环境相关的调查问题遗传力较低，而聚焦于犬主- 犬关系的调查遗传力最低。

图2 犬类血液性状和调查问题的遗传力及GWAS 分析

（A）血液性状和调查问题的SNP 遗传力显示，与犬相关的问题比与犬主- 犬关系（“主人”；[74]）或犬环境相关的问题遗传力更高。括号表示采用双侧置换检验（100,000 次置换）比较中位遗传力。箱线图显示中位数、四分位距和whiskers（1.5× 四分位距），并标注遗传力最高的性状。（B）SNP 遗传力越高的性状，其全基因组关联信号越强（斯皮尔曼相关系数ρ；双侧置换检验，100,000 次置换）。（C）血液性状和调查问题的最强关联。（D）超过指定显著性阈值的犬类GWAS 区域中，与相同或功能相关性状的人类全基因组显著性GWAS 信号重叠的比例。N 表示超过阈值的犬类GWAS 区域数量。（E）苏氨酸水平的GWAS 分析（924 只犬）；（F）同型半胱氨酸水平的GWAS 分析（924 只犬）。

构建首个犬类GWAS目录

为系统探究犬类性状的遗传结构，我们对每个可用表型进行了GWAS，并利用适用于犬类的PheWeb 软件[78]，生成了首个全面的犬类GWAS 目录（broad.io/dogPheWeb）。该交互式平台目前提供154 种调查性状和159 种血液性状的变异- 表型关联数据（图S5）。

为减轻复杂种群结构的影响，我们采用GCTA 软件中的留一染色体交叉混合线性模型（leave-one-chromosome-out mixed linear model，MLMA-LOCO）进行GWAS [52,79]，并纳入遗传力分析中使用的相同协变量。考虑到队列中品种血统的多样性以及相应较短的连锁不平衡[50,52,80]，我们采用标准的人类GWAS 阈值（全基因组显著性p<5.0×10⁻⁸；提示性关联p<1.0×10⁻⁶）[81,82]（数据S3）。

总计识别出298个全基因组显著性位点，平均大小为283kb（标准差= 378kb），定义为聚类邻近SNP（<250kb）且连锁不平衡（r²>0.2）[83,84]。遗传力越高的性状，越可能存在至少一个全基因组显著性位点（逻辑回归OR=9.5，p=0.0036），且总体显著性位点数量更多（拟泊松回归β=3.1，p=2.6×10⁻⁸）（图2B）。体重的结果清晰体现了这一模式，共产生33个不同的关联（图2C），包括已报道的IGF1（p=9.0×10⁻⁵⁰）、HMGA2（p=2.35×10⁻⁴⁹）和LCORL（1.33×10⁻⁴⁵）基因关联[52,55,56,77,85]。我们还识别出一个此前未报道的与microRNA-497邻近的关联（p=1.3×10⁻¹⁴），该microRNA是胰岛素信号基因邻近的IGF1R调控因子[86,87]。

其余153 个调查问题产生76 个显著性关联，其中67% 无显著性关联位点，22% 有1 个显著性关联位点，11% 有2 个及以上。除体重外，调查性状的关联显著性通常低于血液性状（图2C）。最强关联来自调查问题“与熟悉的犬玩耍玩具时表现出攻击性”（GWAS p=4.1×10⁻¹²），该位点在人类中与循环性激素结合球蛋白水平相关（rs79391862-A，p=2×10⁻¹²¹）[88]—— 这一表型主要反映代谢和内分泌状态，但可能间接影响行为。

相比之下，血液性状比调查性状产生更多、更显著的关联（图2C）。在159种血液性状中，识别出189个全基因组显著性关联（平均每种性状1.19±2.35个），其中82种血液性状（52%）至少有一个关联。最强关联来自血浆代谢物，且定位到具有明确生化关联的基因。例如，我们验证了此前报道的GPT基因与循环丙氨酸转氨酶水平的关联[54]。最显著的关联是代谢物3 -羟基异戊酸与MCCC2基因中的SNP——MCCC2 基因编码一种防止该代谢物积累的酶。

血液性状的位点异质性和多效性

犬类血液性状的遗传结构显示出广泛的位点异质性和多效性，与人类代谢组学和血液学GWAS中观察到的模式一致[29,89,90]。在82种至少有一个显著性关联的性状中，37种映射到多个独立位点，其中代谢物3 -羟基异戊酸（17个位点）、5 -甲基胞苷（14个位点）和S -甲基半胱氨酸（13个位点）的位点数量最多。189个全基因组显著性位点中，13 个与多种血液性状相关。

多效性最强的位点是13号染色体上包含ALB（编码白蛋白）的95kb区域。该位点的变异与色氨酸水平（p=1.8×10⁻⁶⁰）、白蛋白/球蛋白比值（p=1.1×10⁻¹⁷）、白蛋白水平（p=1.9×10⁻¹³）和球蛋白水平（p=2.7×10⁻⁸）相关。这一模式具有生物学合理性：约90%的循环色氨酸与白蛋白结合[91,92]，因此改变白蛋白浓度或结合能力的遗传变异会同时影响色氨酸的可利用性。与这一机制一致，ALB位点还与已知结合白蛋白的代谢物相关（图S6），包括L -犬尿氨酸（p=5.6×10⁻⁴⁹）、吲哚- 3 -乙酸（p=4.1×10⁻²⁹）、吲哚- 3 -乳酸（p=4.2×10⁻²⁷）和吲哚- 3 -丙酸（p=2.3×10⁻¹³）[93,94]。

犬类与人类血液性状位点的重叠

犬类血液性状最显著相关的位点与人类研究中涉及的位点同源，表明存在共享的生物学机制，突显了犬类作为转化模型的价值。我们将犬类关联与25 项研究中同源人类GWAS 位点（±250kb）的报道进行了比较[95-120]（数据S4）。在所有181 个犬类位点中，27% 与人类中相同或直接相关性状的关联重叠。当仅关注犬类最显著的信号时，这一重叠率大幅上升：在犬类中p<1×10⁻²⁰的14 个位点中，79% 与人类中相同或直接相关的性状相关（图2D）。

尽管存在这些广泛的相似性，部分代谢通路仍显示出明显的物种特异性差异。在犬类中，最强的GWAS信号之一是苏氨酸水平与编码L -苏氨酸脱氢酶的TDH基因的关联（p=3.9×10⁻⁵⁵）（图2E）。而在人类中，TDH是无功能的假基因，与苏氨酸或任何相关性状无关[120,121]。人类的苏氨酸水平反而与GCKR等位点相关，该基因调节肝脏代谢通量[104,122]。此外，犬类缺乏人类中观察到的色氨酸水平与TDO2基因的强关联——TDO2编码犬尿氨酸通路的限速酶，是系统性色氨酸周转的主要调节因子[123]。在犬类中，最强信号反而映射到编码白蛋白的ALB基因（p=1.8×10⁻⁵⁹），该基因可通过改变结合态与游离态的比例影响色氨酸水平。

品种塑造纯种犬的血液性状

品种血统是犬类血浆代谢物和临床分析物变异的主要决定因素。为评估品种与年龄、体重和性别等个体性状对变异的解释程度，我们对纯种犬（每种血液性状涉及81-87个品种，每个品种357-475只犬）进行了方差分析（II型ANOVA）（数据S5）。品种平均解释25.1±5.1%的变异，高于年龄（1.9±2.5%）、性别（1.8±1.8%）或体重（0.5±0.7%）（图3A）。正如预期，遗传力越高的性状，品种对其变异的解释度越高（斜率= 0.22，p<0.0001），但年龄、性别或体重无此规律（斜率分别为- 0.003、0.016和0.001；p值分别为0.83、0.22和0.95）（图3B）。

图3 品种解释纯种犬血液性状的显著变异

（A）品种对变异的解释度（广义η²，II 型方差分析）高于年龄、性别或体重。采用以表型为随机效应的混合效应模型检验差异；p 值来自品种与各效应的事后对比。（B）品种对变异的解释度随性状遗传力增加而升高（GREML–LDS 约束）。直线表示线性拟合；斜率和p 值来自检验遗传力- 效应交互作用的混合效应模型。（C）纳入品种后，部分性状的体重归因变异降低，与品种- 体重的共享结构一致。展示了包含（+ 品种）和不含（- 品种）品种的模型中体重的变异解释度；突出显示FDR 校正后p<0.05 的性状。（D）年龄、性别和体重对变异的解释度通常低于品种。灰色点表示非显著效应（FDR p>0.05）；标注10 个最强效应。（E）20 种受品种影响最大的性状在年龄、品种、性别和体重之间的变异分配（堆叠条形图）；均无显著的体重效应。（F）临床分析物的犬类与人类效应比较。左侧：已发表的人类风险比；右侧：犬类中品种对变异的解释度。

体重仅能解释所有血液性状变异的一小部分，部分原因是品种与体型的强相关性（图3A）[71]。多重共线性诊断显示，这种相关性通常不会掩盖模型中体重的独立效应（方差膨胀因子VIF=2.72）。然而，对于部分血液性状，纳入品种后会显著减弱体重的表观效应。例如，在不含品种的模型中，体重解释14.5% 的肌酐水平变异（肌酐是反映肌肉总量的代谢物）[124]；而纳入品种后，体重的解释度几乎为零（0.4%），品种的解释度为38.3%。在不含品种的模型中，体重对14 种血液性状的变异解释度超过5%，纳入品种后，体重效应平均降低94.1%±6.7%，剩余解释度仅为0.01%-1.4%（图3C）。

总体而言，96 种血液性状的变异可由品种显著解释（FDR p<0.05），多于年龄（N=68）、性别（N=28）或体重（N=3）（图3D）。例如，年龄对白蛋白- 球蛋白比值（13.9%；FDR p=2.7×10⁻¹¹）和血小板（11.0%，FDR p=8.5×10⁻⁷）的变异解释度较高，这两种性状均与人类衰老相关[125-127]。性别对肉碱水平（11.3%；p=7.5×10⁻⁸）有解释作用—— 这与人类中的性别二态性相关[128]，同时也影响犬尿氨酸（效应量= 9.2%；p=3.1×10⁻⁶）—— 犬尿氨酸是一种受睾酮和雌激素调节的色氨酸代谢物[129]。

品种对被毛类型的选择影响与长寿相关的血液性状

品种对犬类血液性状的影响可能部分源于对美学和功能特征的历史选择。品种对肌酐水平的强效应与对体型和肌肉量的选择一致[130-132]。同样，品种解释32.1% 的淀粉酶水平变异（FDR p=2.0×10⁻⁵）—— 此前研究表明，该性状在品种间的差异源于AMY2B 基因的拷贝数变异，这与对富含淀粉饮食的适应相关[133,134]。

至少有一种代谢物的品种效应可能是对看似无关的美学性状选择的多效性结果。同型半胱氨酸水平是受品种影响第二强的血液性状（40% 的变异可解释，FDR p=2.1×10⁻¹⁰），其关联的两个位点均位于决定被毛特征的基因中：毛发生长模式基因RSPO2（p=2.6×10⁻¹²）和毛发长度基因FGF5（p=1.2×10⁻⁹）[135]（图2F、图4A、B）。这与人类不同——人类的同型半胱氨酸水平与编码同型半胱氨酸γ-裂解酶的CTH基因变异相关[104]。在50只同时具有代谢物测量数据和照片的混种犬与纯种犬中，具有RSPO2相关装饰性被毛特征（胡须和浓密眉毛）的犬以及被毛较长的犬，同型半胱氨酸水平显著更高（图4C；表S3）。此外，与同型半胱氨酸水平最相关的变异也与相应的被毛性状相关（图S7；数据S6）。最后，同型半胱氨酸相关的SNP在不同被毛特征的品种间存在差异（图4D；数据S7），与这两个位点存在选择信号的证据一致[136]。

与衰老相关的血液性状并未受到品种的不成比例影响

仅针对纯种犬的分析显示，受品种影响最强的血液性状并未不成比例地与衰老或长寿相关。我们将与人类两种验证过的血液衰老时钟（Aging.AI [22] 和PhenoAge [137]）相关的32 种性状定义为衰老相关血液性状（表S4）。在控制表型类别（临床分析物vs 血浆代谢物）的线性模型中，这些性状受品种的影响并不比其他血液性状更强（β=0.019，FDR p=0.40）。在这组衰老时钟性状中，品种对其变异的解释度与人类中死亡率风险比更高的性状无相关性（斯皮尔曼相关系数ρ=-0.08，FDR p=0.91）。

图4 血统组成将血液性状变异与品种差异关联起来

（A、B）区域关联图显示，两个同型半胱氨酸位点与装饰性被毛（13 号染色体）和被毛长度（32 号染色体）的经典位点共定位。方框表示在50 只同时具有代谢物测量数据和表型记录的犬中，与同型半胱氨酸水平和指定被毛性状均显著相关的SNP。（C）50 只犬中，按被毛表型分组的标准化同型半胱氨酸水平。点表示个体；箱线图显示中位数和四分位距。p 值来自威尔科克森秩和检验（装饰性被毛）和肯德尔秩相关检验（被毛长度）；红色表示组均值。（D）按被毛表型分层的品种水平上，与同型半胱氨酸最相关的变异的次要等位基因频率；方框总结每组品种的分布。（E）采用线性混合效应回归估算的品种血统对肌酐的效应。血统评分为REML t 统计量；虚线表示无品种效应（t=0）。点表示品种（开放灰色：非显著；实心红色：FDR p<0.05）。直线表示品种效应量与品种平均体重的关系。（F）具有装饰性被毛和长毛的品种，其血统对同型半胱氨酸的效应更强。标注效应最大的品种；文本显示双侧t 检验p 值。

品种血统比例分离品种效应

为识别犬类品种间寿命差异背后的血液性状，我们利用宠物犬种群的遗传多样性，将混种犬纳入分析。这增加了样本量和统计效力，并能够分离品种效应与体型效应。对于每种血液性状，我们采用线性混合效应模型估算单个品种的血统对种群均值偏差的预测程度，得到“品种血统评分”（Breed Ancestry Score），量化品种血统对性状水平的影响[52]。随后，我们检验这些品种血统评分的变异是否可由品种寿命解释（数据S8）。

我们利用已知与品种表型差异相关的血液性状验证了该方法。例如，品种体重解释25.4% 的肌酐品种血统评分变异（肌酐creatinine,是肌肉量的生物标志物）[124]，体型较大品种的血统与较高的肌酐水平相关（n=82；图4E）。同样，具有装饰性被毛和长毛品种的血统与较高的同型半胱氨酸水平相关，与两个关联位点的多效性一致（图4F）。

血液性状变异与品种寿命

总体而言，19 种血液性状的品种血统评分变异可通过体型（17 种）或寿命（3 种，其中1 种同时与体型相关）在包含这两个变量的模型中显著解释（FDR p<0.05）（图5A、图S8、数据S9）。与纯种犬分析不同，体型的效应不再被其与品种的强相关性掩盖。长寿品种的血统与较低的钾水平、较高的钠钾比和较高的氯水平相关（β=-0.51 至0.54；所有FDR p=0.03；图5B）。

图5 品种血统评分区分可预测死亡率的血液性状

（A）品种血统评分变异可通过品种寿命（红色）和/ 或品种体重（蓝色）显著解释的血液性状（FDR p<0.05）。堆叠条形图显示每个协变量解释的变异（表示为线性混合效应模型的广义η²：品种血统评分～品种体重+ 寿命）。（B）（A）中性状的效应估计值，显示为同一模型的β±1.96 标准误，按协变量分面；星号表示FDR p≤0.05。（C）个体犬的死亡率风险比；星号表示FDR p≤0.05。（D）按性状类型分面的个体犬死亡率风险（圈出并标注的点；联合FDR p≤0.05）。（E）个体水平死亡率风险比与品种水平效应的一致性。点表示品种（圈出：FDR p≤0.05；标注：纳入（A））；蓝色直线表示线性拟合，阴影表示95% 置信区间。

图6 短寿品种中水平较高的两种血液性状与较低的存活概率相关

品种寿命效应与死亡率风险比的相关性

我们对957 只犬进行了三年的纵向随访（期间11% 的犬死亡），估算了所有159 种血液性状（包括血浆代谢物和临床分析物）的全因死亡率风险比。采用我们此前开发的方法[44]，在调整入组年龄、体重和性别的基础上，将死亡率风险建模为每种性状的函数（数据S10）。34 种性状与死亡率风险显著相关（FDR p<0.05）（图S9），其中部分性状的水平随品种血统和寿命而变化（图5C、D）。

品种寿命效应与死亡率风险比的相关性（斯皮尔曼相关系数ρ=-0.31；p=5.5×10⁻⁵）强于品种体重效应与死亡率风险比的相关性（斯皮尔曼相关系数ρ=-0.2，p=0.013）（图5E）。换句话说，短寿品种血统预测水平较高的性状，往往也会增加个体犬的死亡风险（“风险性状”），反之则为“保护性性状”。为确认这种关系并非源于Cox 分析中犬的品种血统，我们针对5 种显著的血液性状（FDR p<0.05），通过调整9 个影响最大品种的血统，评估品种血统是否能解释死亡率风险。调整血统后，所有5 种血液性状仍保持显著性（图S10；表S5）。

生长激素相关代谢物并无特殊性

已知受生长激素/ 胰岛素样生长因子1（GH/IGF-1）信号通路影响的血浆代谢物，在寿命相关性方面并不比其他代谢物更显著。在分析的123 种血浆代谢物中，24 种已被证明在啮齿动物模型中会随GH/IGF-1 通路的扰动而变化[138-142]（表S6）。GH 相关代谢物的绝对死亡率风险比更低（均值| z|=0.96 vs. 1.35；威尔科克森检验p=0.0092）。GH 相关代谢物和其他代谢物的品种血统评分与品种寿命的相关性相似（图S11）。

识别风险和保护性血液性状

为识别在品种和个体水平均与寿命相关的性状，我们结合了品种水平和个体水平分析的Z 分数（数据S11），得到13 种FDR p<0.05 的性状（6 种保护性性状，7 种风险相关性状）。例如，球蛋白和钾均为风险相关性状：短寿品种的血统预测其水平较高，且个体犬中水平升高与死亡风险增加相关（图6）。这与人类衰老一致—— 人类中球蛋白升高（HR>1.9）和钾升高（HR>1.4）均会增加死亡风险[143,144]。

8 种血浆代谢物（4 种风险相关，4 种保护性）在联合分析中被证实相关（图5E、图S12）。风险相关代谢物包括：辛二酸suberic acid（一种指示线粒体应激的二元羧酸，随人类年龄增长而积累[145]）、二甲基甘氨酸dimethylglycine,（肾功能障碍的标志物[146]，与全因死亡率增加相关[147,148]），以及乙基丙二酸ethylmalonic acid和琥珀酸succinate（均为线粒体功能障碍的标志物）。保护性代谢物包括：乙醇胺ethanolamine,（通过增强自噬延长酵母、果蝇和哺乳动物细胞系的寿命[149]）、脱氧肉碱deoxycarnitine（肉碱通路通量的生物标志物，其水平降低与衰弱、胰岛素抵抗、心血管疾病和肌肉无力相关[150,151]），以及胍乙酸guanidinoacetate（肌酸的直接前体，肌酸水平降低与肌肉减少症、疲劳和认知衰退相关[152-154]）。

最强的保护性候选物之一是羟脯氨酸hydroxyproline —— 它并非人类中典型的衰老生物标志物。羟脯氨酸主要源于胶原蛋白降解，可能反映富含胶原蛋白组织的周转[155]。长寿品种血统占比更高的犬，羟脯氨酸水平更高，且个体犬中羟脯氨酸水平与较低的死亡风险相关。

并非所有受品种寿命影响的性状都与个体犬的死亡率相关，反之亦然。例如，尽管短寿品种血统占比更高的犬氯水平更低（β=0.62，p=6.5×10⁻⁴，FDR p=0.038；图5B）—— 这与人类研究中氯水平低预测死亡风险高一致[156]—— 但个体犬的氯水平与本研究Cox 分析中的存活率无关联。黄素腺嘌呤二核苷酸flavin adenine dinucleotide（FAD，一种关键的线粒体氧化还原辅因子，其耗竭会导致代谢和细胞衰老[157]）与个体犬的存活率相关（风险比= 1.37；FDR p=0.0036），但品种血统评分与品种寿命无相关性（β=-0.041；p=0.83）。假尿苷Pseudouridine也呈现类似模式：它能强烈预测个体犬的存活率（与人类中死亡风险相关[26] 一致），但品种血统评分与品种寿命无相关性（β=0.172；p=0.29）。

讨论

通过整合数千只犬的基因组、代谢组和临床数据，我们确立了家犬作为解析衰老遗传和生理结构的独特强大系统。我们呈现了首个犬类GWAS 目录，以及迄今为止最大规模的非人类物种血液性状全基因组研究（包括临床分析物和血浆代谢物）。许多血液性状显示出显著的遗传力（与人类相关研究相似），且非常适合GWAS。尽管研究规模远小于典型的人类GWAS，但超过一半的血液性状至少有一个显著性关联。

在犬类和人类中，许多血液性状均受常见的强效变异影响，这使得即使在中等规模队列中也能获得稳健的遗传关联[99,111]。在犬类研究中，我们仅用不到1000 个个体就检测到强信号—— 这一规模比人类研究通常所需的样本量少近一个数量级[14]，这反映了犬类品种的人口统计学历史。严重的瓶颈效应、高近交系数和强烈的人工选择改变了等位基因频率[50,51,158]，显著提升了GWAS 效力[52,80]。这种遗传结构使得血液性状成为有效的内表型，有助于将遗传变异与犬类的疾病风险和衰老等复杂表型联系起来。

犬类血液性状最强的GWAS 信号几乎均与人类中相同或密切相关性状的关联位点重叠，突显了保守的生物学机制，并重申了犬类的转化价值。少数例外情况揭示了犬类与人类之间可能存在分歧的通路，这可能反映了物种特异性的代谢适应。例如，苏氨酸水平与犬类的TDH 基因遗传变异相关，但与人类无关—— 人类的TDH 已进化为无功能的假基因。人类的色氨酸水平与TDO2 基因遗传变异显著相关[122]，但犬类并非如此，这表明该通路存在功能差异。

对犬类品种美学性状的选择和遗传漂变重塑了代谢通路。品种解释了血液性状的显著变异，对89% 的临床分析物和60% 的血浆代谢物有显著影响。对于两种遗传力最高的性状（肌酐h²_SNP=0.88 和同型半胱氨酸h²_SNP=0.87），仅品种就占总变异的近40%。对于肌酐，这可能反映了品种在体型和肌肉量上的差异—— 体型较大品种的血统与较高的肌酐水平相关。然而，个体犬的肌酐水平与全因死亡风险无关联（图5E），表明该生物标志物的品种驱动差异并不会直接转化为寿命差异。与同型半胱氨酸相关的两个位点映射到决定品种特征性被毛性状的基因，且与本队列中的被毛类型相关，这体现了对可见性状选择的多效性效应。然而，同型半胱氨酸水平同样与死亡风险无关联（图5E）。这与人类不同—— 人类中同型半胱氨酸升高与心血管疾病和死亡相关[159]，这可能反映了犬类中心血管疾病的患病率较低[160]。

同型半胱氨酸代谢与被毛性状之间的多效性具有生物学合理性。同型半胱氨酸是甲硫氨酸代谢与半胱氨酸生物合成之间转硫通路的关键中间产物，同型半胱氨酸代谢的遗传疾病包括毛发断裂和脱发症状[161]。半胱氨酸的可利用性可能通过上调角蛋白合成[162] 或影响赋予毛发强度和弹性的二硫键形成[163]，从而影响毛发结构。犬类中与同型半胱氨酸水平和毛发性状均相关的两个位点，映射到已知在毛发生物学中起作用的基因：FGF5（调节毛囊周期[164]）和RSPO2（增强毛囊微环境中的Wnt/β- 连环蛋白信号传导[165]）。

通过分析所有犬类（包括纯种和混种），我们发现了纯种犬分析中无法检测到的体型相关血液性状影响—— 在纯种犬中，体型与品种密切相关。混种犬携带多个品种的嵌合血统，形成了一个自然实验，可评估单个品种血统的增量贡献。利用这一框架，我们发现肌酐与品种体型强相关（与其实作为肌肉量生物标志物的作用一致[124]）。我们还发现，大型品种血统占比更高的犬，胆红素水平往往更高—— 值得注意的是，我们同时发现胆红素水平升高预测犬类全因死亡风险降低（图S9），这与人类情况一致[166]。相比之下，仅针对纯种犬的分析显示，在模型中纳入品种后，体重对肌酐、胆红素或大多数其他血液性状无影响（图3C）。

利用“犬类衰老项目” 规模大、表型深入的纵向队列，我们发现品种血统与衰老之间的关系比“大型犬寿命短” 这一简单观察更为复杂[167]。巨型品种（如大丹犬great Danes）和小型品种（如哈巴狗pugs）的血统均预测球蛋白水平升高（炎症或感染的标志物）和存活率降低[143]（图6A）。不同体型的短头品种（如哈巴狗pugs和斗牛犬bulldogs）的血统与较高的钾水平相关（图6B），而钾水平升高反过来预测存活率降低（图6D）。这可能反映了先天性呼吸障碍的后果—— 慢性气道阻塞可导致呼吸性酸中毒和继发性血钾升高[168,169]。然而，非短头小型品种（如博美犬pomeranians Airedale terriers）和较大型品种（如艾尔谷梗）的血统也与较高的钾水平相关（图6B），表明可能存在其他作用机制。

图6 短寿品种中水平较高的两种血液性状与较低的存活概率相关

（A）球蛋白水平的品种血统评分（REML t）与品种寿命的关系；（B）钾水平的品种血统评分与品种寿命的关系。点表示品种（红色：FDR p≤0.05；灰色：FDR p>0.05），标注效应最强的品种；红色直线表示线性拟合，斯皮尔曼相关系数（p）标注于图中。灰色虚线表示无血统效应。图片展示具有最强效应的代表性品种。（C、D）根据入组年龄（老年：≥7 岁；非老年：<7 岁）和球蛋白或钾的标准化丰度（高：上四分位数）分层的犬类Kaplan-Meier 生存曲线，以入组后的年数为横轴[44]。

这些结果突显了混合映射（admixture mapping）在增强犬类遗传研究效力方面的潜力—— 前提是能够解决高维血统带来的计算挑战。混合映射利用混血个体的局部血统识别与性状相关的位点，但现有方法仅适用于建模数量相对较少的来源种群。在“犬类衰老项目” 中，对80% 的犬进行血统建模需要纳入至少70 个不同的来源种群（品种）（图S13）。能够适应这种高维血统空间的新方法，可能进一步提升犬类GWAS 的效力（尤其是对具有强品种效应的血液性状），并识别出与寿命相关变异的来源种群。

总之，我们的结果表明，可利用犬类的遗传多样性识别衰老生物标志物、评估其在个体中的预测价值，并明确保守的生理通路。这种双水平方法为解析衰老生物学机制和优先选择干预靶点提供了框架。家犬具有与人类相似的生理复杂性，且寿命较短，是衰老研究中兼具发现和转化价值的卓越模型。