猫科（Felidae）体型进化的基因组学解析

《Ecology and Evolution》：Genomic Insights Into the Body Size Evolution in Felidae (Mammalia: Carnivora)

　　

体型是哺乳动物最具基础性的适应性特征之一，但控制猫科（Felidae）中所观察到的大范围体型多样化的遗传架构仍未得到充分解析。在本研究中，研究人员利用来自17个物种的高质量参考基因组（包括两个新重新注释的基因组），首次对猫科体型进化进行了全面的比较基因组学分析。利用系统发育广义最小二乘（phylogenetic generalized least squares, PGLS）分析，研究人员鉴定了35个体型相关基因（body-size-associated genes, BAGs），其中大部分在功能上与骨骼发育、甲状腺激素合成和脂质代谢有关。功能富集分析显示，大体型猫科中的正选择基因（positively selected genes, PSGs）主要富集于DNA修复通路，而小体型猫科中的PSGs则主要参与能量代谢和脂质代谢过程。总之，核心信号通路和关键功能基因的适应性进化共同驱动了猫科体型多样化。这项研究为哺乳动物体型如何进化带来了新的分子学见解，同时生成了一个有价值的基因组数据集，供后续功能研究使用。

**一、研究背景与科学问题**
体型是哺乳动物最具基础性的适应性特征之一，与物种的生活史、生态位以及灭绝风险密切相关。关于哺乳动物体型进化方向，学术界存在多种假说，如“柯普法则”认为新生代哺乳动物总体趋向体型增大，另有研究主张哺乳动物存在单一或多个最适体型。猫科（Felidae）隶属于食肉目，现包含2个亚科、14属46种，除南极洲外广泛分布。猫科动物在体型上表现出惊人的差异：最小的锈斑豹猫和黑足猫（*Felis nigripes*）体重仅1至2千克，而最大的虎（*Panthera tigris*）某些亚种体重可达近300千克，两者相差超过200倍。因此，猫科是研究哺乳动物体型进化的理想类群。此前基于化石记录的研究推断猫科存在小体型和大体型两个最适演化峰值，但体型差异背后的遗传机制仍未被系统阐释。近年来高质量基因组资源的快速增长，使得从分子层面解析体型进化成为可能。Porto等人基于化石和贝叶斯框架推测体型增大与捕食大型猎物有关，但基因组层面的证据仍属空白。为此，研究人员利用17个高质量猫科基因组展开比较基因组学分析，以期揭示体型进化的分子机制。

**二、研究内容与总体结论**
研究人员对17个猫科物种的高质量参考基因组进行分析，其中对黑足猫和丛林猫（*Felis chaus*）两个物种的基因组进行了重新注释，并以小家鼠（*Mus musculus*）基因组作为外群。通过系统发育广义最小二乘（phylogenetic generalized least squares, PGLS）分析，鉴定出35个体型相关基因（body-size-associated genes, BAGs），多数与骨骼发育、甲状腺激素合成和脂质代谢有关。对大体型猫科（狮*Panthera leo*和虎*P. tigris*）和小体型猫科（黑足猫、豹猫*Prionailurus bengalensis*、家猫*Felis catus*和乔氏猫*Leopardus geoffroyi*）分别进行正选择分析，发现大体型猫科中的正选择基因（positively selected genes, PSGs）显著富集于DNA修复通路，而小体型猫科中的PSGs主要富集于能量代谢和脂质代谢。此外，研究人员在大型猫科分支鉴定出324个快速进化基因（rapidly evolving genes, REGs），富集于G蛋白偶联受体信号通路等。该研究表明核心信号通路和关键功能基因的适应性进化共同驱动了猫科体型分化，为哺乳动物体型进化提供了新见解，同时构建了有价值的基因组数据集。

**三、主要研究方法**
研究人员从NCBI数据库获取17个猫科物种基因组序列，并以小家鼠基因组作为外群；体重数据来自PanTHERIA数据库及《Fowler's Zoo and Wild Animal Medicine》，依据体重和猎物体型将物种划分小、中、大三种体型类别。利用RepeatModeler、RepeatMasker、Augustus等工具进行基因组注释，整合从头预测与同源预测证据构建非冗余基因集。采用OrthoFinder结合Diamond鉴定单拷贝直系同源基因，使用MUSCLE进行蛋白比对、PAL2NAL转换为密码子比对，基于RAxML构建系统发育树。运用PAML codeml计算基因进化速率（ω），在R环境中进行PGLS关联分析，并通过分支位点模型和分支模型检测PSGs和REGs，最后进行GO和KEGG富集分析。

**四、研究结果**

**4.1 两个猫科物种的基因组注释结果**
研究人员对*F. nigripes*和*F. chaus*基因组进行从头注释，重复序列含量分别为40.31%和40.67%，蛋白质编码基因数约为16,640和16,577个。BUSCO评估显示所有物种完整基因比例均超过85%，满足下游分析要求。

**4.2 直系同源基因筛选与系统发育树重建**
利用OrthoFinder在包括外群在内的18个基因组中鉴定出20,704个直系同源组，其中7875个为严格单拷贝基因。基于这些基因构建的系统发育树与既有外显子数据系统关系一致。

**4.3 体型相关基因的基因组扫描**
通过PGLS分析，以FDR校正p<0.05为阈值鉴定出35个BAGs，其中ABI3BP、CNR2、ERCC6、HERC1、IL1A、PPARGC1B、SLC30A10和TPO等既往已被证实参与哺乳动物生长、骨骼发育和脂质代谢调控。GO富集显示BAGs显著富集于MLL1复合体（GO:0071339）、ATP依赖性染色质重塑活性（GO:0140658）、甲基化（GO:0032259）、钙通道活性（GO:0005262）和钙离子转运（GO:0006816）等功能类群；KEGG富集显示其显著富集于破骨细胞分化（ko04380）、甲状腺激素合成（ko04918）和雌激素信号通路（ko04915）等13条通路。

**4.4 正选择分析**
在狮和虎中鉴定出92个PSGs，GO富集于DNA修复（GO:0006281）、Ras蛋白信号转导负调控（GO:0046580）及周期蛋白依赖性丝氨酸/苏氨酸激酶调节活性（GO:0016538）；KEGG富集于DNA修复和重组蛋白（ko03400）、JAK–STAT信号通路（ko04630）、矿物质吸收（ko04978）和消化系统（ko09154）。在小体型猫科中鉴定出207个PSGs，GO富集于脂肪酸代谢过程（GO:0006631）、脂质分解代谢过程（GO:0016042）、NAD⁺结合（GO:0070403）及Smo信号通路（GO:0007224）；KEGG富集于脂质代谢（ko09103）、类固醇生物合成（ko00100）和Hedgehog信号通路（ko04340）。

**4.5 快速进化基因分析**
在大型猫科分支中，利用分支模型鉴定出324个REGs，富集于G蛋白偶联受体信号通路（GO:0007186）、G蛋白偶联受体活性（GO:0004930）和肌动蛋白结合（GO:0003779）。

**五、讨论与结论**
讨论指出，猫科体型差异通常被理解为通过性状替代和营养生态位分化减少种间竞争的适应性策略。BAGs富集于MLL1复合体和甲状腺激素合成等通路，表明骨骼发育和甲状腺激素调控是猫科体型分化的核心因素。ERCC6、ABI3BP、HERC1、CNR2、IL1A、PPARGC1B、SLC30A10和TPO的功能证据进一步支持这些基因在生长、骨骼稳态和脂质代谢中的关键作用。大体型猫科PSGs在DNA修复通路中的富集为佩托悖论（Peto's paradox）提供了分子诠释，即大体型长寿物种可能通过增强DNA修复能力降低突变累积和癌症风险。REGs富集于GPCR信号通路，与Hippo通路协同调控器官和机体大小。小体型猫科PSGs富集于能量和脂质代谢通路，与小体型哺乳动物较高基础代谢率的需求相符。本研究也存在局限性，如大体型物种仅包含狮和虎、体型分组具有一定主观性，以及部分BAGs缺乏功能实验验证。

翻译研究结论部分：本研究对两个猫科物种的基因组进行了重新注释，并首次对猫科体型变异和进化进行了基因组学分析。结果表明，猫科体型进化受核心信号通路和个体关键基因的适应性变化共同驱动。比较基因组学分析显示，猫科体型进化与骨骼发育、脂质调控和生长激素调节密切相关。在通路水平上，这种关联体现在MLL1复合体和甲状腺激素合成等富集功能类群中。在基因水平上，研究人员鉴定了多个BAGs，包括ABI3BP、CNR2、ERCC6、PPARGC1B和SLC30A10。此外，部分结果与其他哺乳动物体型进化研究结论一致，表明调控体型的分子机制在哺乳动物谱系间可能具有共性，为今后解析复杂表型变异的分子机制提供参考。