编辑推荐：
摘要翻译
Apoidea是拥有超过20,000个已识别物种的最重要传粉昆虫类群之一，其生存环境复杂多样。不同地理区域中Apoidea的蜜源植物组成存在差异，促使该昆虫类群在形态和生理水平上演化出多种适应性特征以增强其生存能力。消化酶的适应性进化对于Apoi
摘要翻译
Apoidea是拥有超过20,000个已识别物种的最重要传粉昆虫类群之一，其生存环境复杂多样。不同地理区域中Apoidea的蜜源植物组成存在差异，促使该昆虫类群在形态和生理水平上演化出多种适应性特征以增强其生存能力。消化酶的适应性进化对于Apoidea高效吸收营养至关重要；然而，蜜蜂消化酶多样化的分子基础仍不清楚。本研究分析了来自Apoidea六个科的GH13家族中代表Hbg1、Hbg2和Hbg3三个亚家族的63条α-葡萄糖苷酶序列。系统发育分析表明，α-葡萄糖苷酶基因家族在Apoidea总科的早期演化过程中经历了两次分歧事件，产生了GH13家族成员Hbg1、Hbg2和Hbg3；与祖先谱系相比，这三个成员在表达模式和酶功能上均表现出显著分歧。蜜蜂中的GH13基因具有显著的结构多样性。此外，亚家族之间观察到明显的功能分化，并在特定谱系和密码子位点伴随正选择信号。位于α-葡萄糖苷酶配体结合区附近的28个位点受到正选择和功能分化，这可能与α-葡萄糖苷酶的配体结合及水解特性有关。总体而言，α-葡萄糖苷酶基因家族在Apoidea中经历了由正选择驱动的适应性进化，从而导致了该基因家族的扩张和多样化。
论文解读文章
**研究背景与科学问题**
Apoidea总科隶属于膜翅目，是陆地生态系统中重要的传粉者，对维持生态平衡和被子植物多样性具有关键作用。蜜蜂类昆虫的起源时间（约140–110百万年前）与早期被子植物花粉的出现（约139–131百万年前）及其快速多样化时期高度吻合，为研究昆虫与植物的协同进化提供了重要证据。在长期演化过程中，Apoidea在形态、生理和分子水平上形成了多层次的适应机制。碳水化合物是昆虫的主要能量来源，其中α-葡萄糖苷酶是糖苷水解酶（GH）超家族中的一种保守水解酶，在昆虫中主要分布于中肠等消化组织，催化淀粉和寡糖的水解以提供能量。根据催化结构域的氨基酸序列相似性，α-葡萄糖苷酶被划分为GH4、GH13、GH31、GH97和GH122五个家族，昆虫α-葡萄糖苷酶主要属于GH13家族。目前，对Apoideaα-葡萄糖苷酶的研究主要集中在西方蜜蜂（*Apis mellifera*），已鉴定的主要编码基因为*Hbg1*、*Hbg2*和*Hbg3*，它们在组织表达模式和酶学特性上存在显著差异。然而，这些功能差异是否代表整个Apoidea的普遍进化模式，以及该基因家族的起源、重复历史、功能分化和选择压力，尚未被系统研究。因此，本研究旨在从比较分子进化角度，系统分析Apoidea的GH13α-葡萄糖苷酶基因家族的起源、分歧和选择压力，揭示*Hbg1*、*Hbg2*和*Hbg3*的不同进化模式。
**研究内容与主要结论**
研究人员收集了来自Apoidea六个科（Apidae、Andrenidae、Colletidae、Halictidae、Melittidae和Megachilidae）的63条GH13家族α-葡萄糖苷酶序列，综合运用系统发育分析、基因结构比较、染色体共线性分析、功能分化检测和选择压力分析等方法，阐明了该基因家族的进化历史。结果表明，α-葡萄糖苷酶基因家族在Apoidea早期演化中经历了两次分歧事件，产生了*Hbg1*、*Hbg2*和*Hbg3*三个亚家族；*Hbg1*与*Hbg3*位于同一染色体上呈串联排列，而*Hbg2*位于另一染色体上。基因结构分析显示，*Hbg2*仅含一个外显子，而*Hbg1*和*Hbg3*各含八个外显子，表明*Hbg2*可能通过逆转录转座机制产生。功能分化分析显示各亚家族间存在显著的功能分歧，其中*Hbg1*与*Hbg3*之间的功能分化系数最高。选择压力分析检测到*Hbg1*、*Hbg3*以及一些谱系特异性重复分支受到正选择，而*Hbg2*相对保守。在配体结合区附近有28个位点同时受到正选择和功能分化，可能参与底物识别和催化特性的改变。该研究揭示了GH13α-葡萄糖苷酶在基因家族多样化后经历了不对称的选择压力和功能分化，为理解Apoidea消化酶基因家族的适应性进化提供了新的分子证据。
**关键技术方法**
本研究使用的样本包括数据库中的基因组/转录组序列以及野外采集样本。对于数据库缺失的物种，研究人员通过PCR扩增和Sanger测序获得全长序列（如*Macropis* sp.和*Hylaeus* sp.），并从NCBI SRA数据库下载五个物种（*Andrena cineraria*、*Andrena camellia*、*Colletes gigas*、*Melitta haemorrhoidalis*、*Dasypoda hirtipes*）的转录组数据进行组装。主要技术方法包括：使用TBLASTN同源搜索和Trinity组装获得候选序列；采用MAFFT进行多序列比对，Gblocks提取保守位点；利用RAxML和MrBayes构建最大似然（ML）和贝叶斯（BI）系统发育树；使用MEGA7中的RelTime方法估算分歧时间；利用DIVERGE 3.0进行Type I功能分化分析；使用PAML 4.8中的CodeML进行位点模型、分支模型和分支-位点模型的选择压力分析；通过Phyre2进行同源建模，3DLigandSite预测配体结合位点。
**研究结果**
**3.1 序列获取与比对**
通过整合数据库检索、RT-PCR扩增和转录组组装，共获得来自Apoidea六个科的63条保守GH13家族核苷酸序列。多序列比对显示，GH13家族α-葡萄糖苷酶序列相似性高，催化相关的α-Amy结构域高度保守，表明该家族在扩张和分化后仍保留了基本催化功能。相比之下，*Hbg1*、*Hbg2*和*Hbg3*在α-Amy结构域的非必需区域存在长度差异，可能导致底物和活性的差异。
**3.2 系统发育关系**
基于贝叶斯法构建的系统发育树显示，膜翅目GH13家族分为*Hbg1*、*Hbg2*和*Hbg3*三个亚家族，其中*Hbg1*和*Hbg3*聚在一起，*Hbg2*单独成支，表明*Hbg2*分歧更早。此外，在部分类群中发生了谱系特异性重复，如切叶蜂科*Osmia*属中*Hbg1*发生重复形成*Hbg11*和*Hbg12/13*，蜜蜂科中*Hbg3*重复为两个拷贝，反映了食物需求驱动的适应性进化。
**3.3 染色体共线性分析**
以*A. mellifera*的*Hbg1*、*Hbg2*和*Hbg3*为参照，发现*Hbg1*和*Hbg3*在同一染色体上串联排列，而*Hbg2*位于不同染色体上，表明两者起源于不同的基因重复机制：前者为局部串联重复，后者为染色体间重复或易位。鳞翅目和双翅目昆虫的GH13基因也分布在两条染色体上，提示类似重复模式在远缘昆虫谱系中保守。
**3.4 基因结构分析**
所有膜翅目GH13家族的*Hbg2*基因仅含一个外显子，而*Hbg1*和*Hbg3*各含八个外显子，支持*Hbg2*可能由反转录转座产生。*Hbg1*亚家族外显子组成变异较大，仅三个外显子长度保守；*Hbg3*亚家族外显子组成相对保守，多数成员具有六个保守外显子。蜜蜂科中*Hbg3*重复产生的两个拷贝（*Hbg31*和*Hbg32*）基因结构高度相似。
**3.5 进化分歧时间**
RelTime分析显示，Apoidea GH13家族的早期分歧约在145.08百万年前，与被子植物早期辐射时间重叠；*Hbg2*在约123百万年前分歧，*Hbg1*/*Hbg3*重复约在107.33百万年前，对应晚白垩世被子植物爆发期。后续谱系特异性重复事件，如*Apis*和*Bombus*属中的重复约在17–20百万年前，可能与中新世气候变迁和寄主植物多样化有关。
**3.6 功能分化分析**
DIVERGE 3.0的Type I功能分化分析显示，*Hbg1*、*Hbg2*和*Hbg3*两两之间的功能分化系数θ均显著大于0（p<0.01），表明各亚家族间存在显著功能分化。其中*Hbg1*与*Hbg3*的功能分化系数最高（θ=0.4096），反映出这两个亚家族在重复后经历了显著的功能进化。多数潜在功能分化位点位于GH13家族保守催化结构域内，可能影响α-葡萄糖苷酶的催化活性和酶学特性。
**3.7 选择压力分析**
位点模型表明，整个GH13家族整体受到强烈纯化选择（M0模型ω=0.23878），但M8模型检测到15个潜在正选择位点，且M8显著优于M7。分支模型显示，*Hbg1*谱系ω值高达33.7879，检测到20个正选择位点；*Hbg2*（ω=0.026）、*Hbg3*（ω=0.0025）及*Hbg1*-*Hbg3*祖先谱系（ω=0.0024）均受纯化选择。*Osmia*属特异性重复形成的*Hbg12*分支（ω=2.13336）也受正选择。分支-位点模型在*Hbg1*和*Hbg3*谱系分别检测到20个和31个正选择位点，而*Hbg2*未检测到正选择。广食性社会蜜蜂的*Hbg1*和*Hbg3*部分位点受正选择，其中257H、295V、480G和510S同时是潜在功能分化位点，提示食性变化驱动了α-葡萄糖苷酶的功能改变。
**3.8 同源建模及关键位点空间定位**
以家蚕蔗糖水解酶BmSUH（PDB ID: 6LGH）为模板构建了AmHbg1、AmHbg2和AmHbg3的三维结构模型，模板覆盖率分别为91%、95%和96%，建模置信度均为100%。3DLigandSite预测AmHbg1含有一个镁离子结合域和两个糖结合域。将28个同时受正选择和功能分化的位点映射到结构上，发现它们主要分布在糖和镁离子结合区域附近，表明这些位点可能影响配体结合和蛋白质空间构象，进而参与底物特异性与催化特性的分化。
**讨论与结论**
讨论部分指出，Apoidea GH13家族的祖先基因产生于早白垩世（约144百万年前），随后经历两次主要分歧和多次谱系特异性重复，同时伴随基因丢失事件（如切叶蜂科丢失*Hbg3*）。不同亚家族表现出明显的功能分化和差异化的选择压力，正选择位点集中在配体结合区附近，提示其与食物资源利用策略相关。然而，高ω值可能受有限替换数和比对不确定性影响，这些位点应视为候选功能残基，需进一步通过结合亲和力和配体相互作用实验验证。研究结论表明，α-葡萄糖苷酶基因家族在Apoidea早期进化中经历了两轮多样化，产生*Hbg1*、*Hbg2*和*Hbg3*，它们具有不同的表达模式和功能特化；现存GH13家族分布在两条染色体的不同区域，至少发生五次基因重复事件，包含大量促进功能分化的特定位点；*Hbg1*、*Hbg3*及某些亚支受正选择影响，而不同位点承受不同选择压力；28个氨基酸位点同时受正选择并导致功能分化，且集中在蛋白-配体结合区附近，影响底物特异性、催化活性和三维构象。谱系特异性重复基因在酶活性上的进一步分化，推动了GH13家族的扩展和功能多样化。
