多组织转录组揭示荒漠草地蝗虫意大利蝗解毒代谢基因的系统鉴定与组织表达特征

《Insects》：Multi-Tissue Transcriptomics for Systematic Identification and Expression Analysis of Detoxification Genes in Calliptamus barbarus

　　

意大利蝗（Calliptamus barbarus）是中国新疆荒漠及半荒漠草原的优势害虫，但其解毒代谢基因库及组织水平上的功能布局此前尚不清楚。研究人员对四个关键解毒组织（胃盲囊、中肠、马氏管和脂肪体；各三个生物学重复）进行了转录组测序，并基于黑腹果蝇（Drosophila melanogaster）和飞蝗（Locusta migratoria）的参考蛋白序列，采用BLAST和HMMER联合搜索，共鉴定出五个家族1602个解毒代谢基因（羧酸酯酶COE 522个、细胞色素P450 470个、UDP-糖基转移酶UGT 276个、ABC转运蛋白233个、谷胱甘肽S-转移酶GST 101个）。与飞蝗和沙漠蝗（Schistocerca gregaria）的系统发育比较显示，外源物代谢相关亚家族（CYP3、CYP4、α-酯酶、Sigma-GST、UGT305及ABCB/ABCC）在意大利蝗中组成比例显著偏高且呈物种特异性聚集，而承担基础生理功能的亚家族则较为保守。在188个组织表达基因中，表达呈现跨组织的功能分工：胃盲囊主导II相缀合反应（GST、UGT）和水解反应（COE），而马氏管主导I相氧化反应（P450）和III相外排（ABC）。对覆盖全部五个家族的11个组织偏向基因进行RT-qPCR验证，11/11个基因均重现了转录组方向（双侧符号检验p = 0.001），且与FPKM倍数变化定量一致（Pearson r = 0.98），确认了Ugt49B1和MRP-PO（胃盲囊和中肠，超过25倍）具有最强的组织偏向性。这些发现为理解意大利蝗的外源物代谢和环境适应提供了解毒基因资源及组织水平功能图谱。

意大利蝗解毒代谢基因家族的系统鉴定、进化特征与组织功能分工研究解读

一、研究背景与科学问题

昆虫代谢解毒外源有毒化合物（包括合成杀虫剂和天然植物毒素）的能力，是其环境适应性和抗药性进化潜力的核心分子基础。昆虫解毒代谢系统通常包含三个功能阶段：I相反应由细胞色素P450单加氧酶（P450）和羧酸酯酶（COE）对外源化合物进行氧化或水解修饰；II相反应由谷胱甘肽S-转移酶（GST）和UDP-糖基转移酶（UGT）催化内源极性分子与底物的缀合，提高产物水溶性并降低生物活性；III相反应由ATP结合盒（ABC）转运蛋白作为外排泵，利用ATP水解释放的能量将终代谢产物跨膜运出细胞。这五个解毒基因家族的组成多样性和表达可塑性直接决定昆虫对外源化合物的代谢清除能力，也是代谢性杀虫剂抗性产生的关键分子基础。

近年来，高通量测序技术推动了许多昆虫物种解毒基因家族的系统鉴定与比较分析。在鳞翅目中，草地贪夜蛾（Spodoptera frugiperda）玉米品系与水稻品系的比较基因组分析揭示了CYP6AE亚家族的品系特异性扩张及UGT33B/UGT33F的品系间收缩；核桃举肢蛾（Atrijuglans aristata）作为单食性昆虫，其基因组解毒基因总数显著少于多食性近缘种。在膜翅目中，西方蜜蜂（Apis mellifera）全解毒基因组的表达谱分析确认马氏管、脂肪体和中肠为解毒代谢核心组织。在直翅目中，飞蝗（Locusta migratoria）作为经典解毒代谢研究模型，其P450基因CYP408B1和CYP409A1在溴氰菊酯诱导下显著上调，中肠特异性高表达的羧酸酯酶LmCesA1基因簇通过物理螯合和催化水解双重机制高效降解有机磷制剂。

意大利蝗（Calliptamus barbarus）隶属于直翅目斑腿蝗科，分布于从温带荒漠到草甸草原的多种草地类型，是新疆平原荒漠草原的优势害虫。该虫为多食性蝗虫，主要取食蒿属（菊科）和藜科植物，这些寄主植物普遍含有结构多样的次生代谢物（萜类、生物碱、黄酮类等）。近年来，在全球气候变化影响下，新疆草地蝗虫适生区格局持续变迁，长期高频率化学防治也导致草地蝗虫产生不同程度的抗药性。然而，与飞蝗和沙漠蝗相比，意大利蝗解毒代谢基因家族的组成、进化模式和组织表达特征完全未知。作为栖息于植物多样性较高的荒漠草原环境的多食性蝗虫，其解毒基因家族的组成规模是否适应其广食性取食范围，五个解毒基因家族在关键代谢组织中是否存在功能分工，阐明这些问题对于理解蝗虫解毒代谢系统的进化规律和制定针对性防控策略具有重要意义。

二、研究设计与技术方法

研究人员于2025年7月在中国新疆博乐市采集意大利蝗成虫，实验室条件下饥饿24小时排空肠道内容物后，在冰上分别解剖收集胃盲囊、中肠、马氏管和脂肪体四个关键解毒代谢组织。共解剖15头成虫，其中9头随机分配为三个生物学重复，每个重复由三头蝗虫的对应组织混合而成。总RNA提取后构建cDNA文库，在Illumina NovaSeq X Plus平台进行150 bp双端测序。测序数据经fastp质控后，用Trinity对各样本分别进行从头组装，CD-HIT去冗余构建非冗余unigene参考集，BUSCO评估组装完整性，TransDecoder预测开放阅读框。解毒基因家族鉴定采用BLASTp和HMMER联合策略，以飞蝗和黑腹果蝇的参考蛋白序列进行同源搜索，并以Pfam数据库结构域模型验证。系统发育树用MAFFT比对、trimAl修剪、IQ-TREE最大似然法构建，1000次自举评估节点置信度。组织表达谱通过RSEM定量计算FPKM值，单因素方差分析（ANOVA）结合Tukey HSD事后检验鉴定组织特异性表达基因，DESeq2进行两两组织间差异表达分析。RT-qPCR以β-actin为内参基因，对覆盖五个家族的11个组织偏向基因进行独立验证。

三、主要研究结果

3.1 转录组测序与组装质量评估

四个组织共12个样本的转录组测序结果显示，所有组织样本的完整BUSCO比例均高于88%，马氏管样本最高（94.0%），胃盲囊样本最低（88.5%）。单拷贝基因比例为55%-70%，片段化基因比例不超过4.1%，缺失基因比例不超过7.4%。所有样本Q20大于98.6%，Q30大于93.4%，GC含量为40.26%-44.70%。这些指标表明转录组组装具有较好的完整性和连续性。

3.2 解毒代谢基因家族的鉴定与种间比较

通过BLAST和HMMER联合鉴定及低置信度序列过滤，三个物种共鉴定出3017个解毒代谢基因。意大利蝗共获得1602个解毒基因，包括COE 522个（32.6%）、P450 470个（29.3%）、UGT 276个（17.2%）、ABC 233个（14.5%）和GST 101个（6.3%）。飞蝗和沙漠蝗分别鉴定出582个和833个解毒基因。P450和COE在三个物种中均为最大的两个家族，ABC比例在三个物种间相对稳定（14.5%-20.1%）。Fisher精确检验显示，与飞蝗相比，意大利蝗COE α-酯酶亚家族比例显著升高（OR = 2.05, FDR = 0.006），UGT305亚家族比例显著升高（OR = 1.99, FDR = 0.046），而表皮酯酶亚家族比例显著降低（OR = 0.41, FDR = 0.023）；与沙漠蝗相比，GST Sigma亚家族比例显著升高（OR = 3.62, FDR = 0.007），CYP2亚家族比例显著降低（OR = 0.32, FDR = 0.001）。

3.3 系统发育分析

ABC基因家族系统发育树包含493条蛋白序列，所有基因聚为ABCA-ABCH八个亚家族，其中ABCB（70个）和ABCC（72个）在意大利蝗中成员数最多，ABCB亚家族形成显著物种特异性簇，最大簇含39个ABCC成员。ABCE、ABCD和ABCH亚家族在三个物种间成员数有限且相似，序列高度保守。

COE基因家族系统发育树包含897条蛋白序列，分为胞内催化类、分泌催化类和神经发育类三大功能类群。α-酯酶构成最大单系群，意大利蝗含364个基因（占COE总数的69.7%），远超飞蝗（80个）和沙漠蝗（157个），形成209个物种特异性簇。β-酯酶呈现相似模式。乙酰胆碱酯酶（AChE）和神经连接蛋白（Neuroligin）等非催化亚家族在物种间数量相近且保守。

GST基因家族系统发育树包含196条蛋白序列，分为Delta、Sigma、Epsilon、Omega、Theta和Zeta六个胞质亚家族。Sigma亚家族在意大利蝗中最丰富（50个，占GST总数的49.5%），远超飞蝗（12个）和沙漠蝗（13个）。Delta亚家族含30个基因，Epsilon形成独立小分支，Omega、Theta和Zeta位于树基部且物种间差异小。

P450基因家族系统发育树包含900条蛋白序列，分为CYP2、CYP3、CYP4和线粒体P450四个单系群。CYP3群构成最大分支，意大利蝗含273个基因（占P450总数的58.1%），远超飞蝗（101个）和沙漠蝗（137个）。CYP4群在意大利蝗中成员数也较多（147个）。线粒体P450群形成独立且高度保守的小分支。CYP2群在意大利蝗中仅含25个基因，显著少于沙漠蝗（37个）。

UGT基因家族系统发育树包含531条蛋白序列，聚为UGT50、UGT36和UGT305三个单系群。UGT305亚家族在意大利蝗中构成最大分支（171个，占UGT总数的62.0%），远超飞蝗（45个）和沙漠蝗（75个）。UGT50含69个基因，UGT36在三个物种间差异小。

综合五个家族的系统发育分析，解毒相关亚家族（ABCB、ABCC、CYP3、CYP4、α-酯酶、β-酯酶、Sigma-GST、UGT305）在意大利蝗中倾向于紧密聚集形成物种特异性簇，而保守亚家族（ABCE、ABCD、Mito-P450、AChE、UGT36）成员少、分支短、序列高度保守。

3.4 解毒代谢基因的组织表达谱

对188个具有完整表达数据的解毒代谢基因进行组织表达谱分析。GST家族在胃盲囊（平均FPKM = 56.51）和马氏管（平均FPKM = 55.34）中表达最高，脂肪体最低（平均FPKM = 6.71）。GstT2在马氏管中表达最高（ANOVA p = 0.009），GST-se在中肠显著富集（ANOVA p = 0.002）。

COE家族在胃盲囊（平均FPKM = 27.81）和脂肪体（平均FPKM = 17.24）中表达相对较高。CG1128aE9_alpha_esterase在胃盲囊中相对于脂肪体、中肠和马氏管均显著上调，CG4757_integument_esterase在马氏管中相对于胃盲囊显著上调。

P450家族整体在马氏管中表达最高（平均FPKM = 8.63），其次为胃盲囊（平均FPKM = 7.10）。CYP3_Cyp6g1在马氏管中显著富集（ANOVA p = 0.036），CYP3_Cyp9f2在中肠显著上调（ANOVA p = 0.050）。

ABC家族在马氏管中表达最高（平均FPKM = 14.33），其次为中肠（平均FPKM = 11.04）。两两DESeq2比较鉴定出10个ABCB、ABCC和ABCG亚家族成员具有组织偏好性表达。

UGT家族主要在胃盲囊（平均FPKM = 9.07）和中肠（平均FPKM = 5.95）表达。Ugt35E2在马氏管中表达极显著高于其他三个组织（ANOVA p < 0.001），Ugt50B3在胃盲囊中显著富集（ANOVA p = 0.045）。

3.5 解毒代谢基因的差异表达分析

以Padj < 0.05和|log₂FC| > 1为阈值，在六个两两组织比较中筛选差异表达解毒基因。胃盲囊vs脂肪体比较中，CG1128aE9_alpha_esterase、Mito_shd、CYP3_Cyp6a20、GstS1和GstD1显著上调。胃盲囊vs马氏管比较中，Ugt49B1、CG1128aE9_alpha_esterase、ABCC_CG31792和ABCG_CG32091显著上调。中肠vs马氏管比较中，ABCC_MRP-PO和ABCC_Cftr显著上调。马氏管vs脂肪体比较中，Mito_shd、ABCG_CG11069、Ugt35A1、ABCG_CG31121和ABCG_CG5853显著上调。差异表达基因的组织分布模式表明，COE和部分GST基因在胃盲囊特异性上调，ABC和线粒体P450在马氏管中高度富集，特定UGT成员主要在中肠表达。

3.6 RT-qPCR验证组织偏向性表达

对11个覆盖全部五个家族的解毒基因进行RT-qPCR验证，所有基因在转录组鉴定的目标组织中的平均相对表达量均高于对照组织（方向一致性11/11，双侧符号检验p = 0.001）。Welch t检验显示10/11个基因的对照-目标组织差异达到显著水平（p = 0.0006-0.015），唯一例外为CG4757（p = 0.113）。RT-qPCR倍数变化范围为2.85倍（CG4757）至27.27倍（Ugt49B1），与RNA-seq FPKM倍数变化（3.49倍和24.73倍）高度一致，两种方法log₂倍数变化的中位绝对偏差为0.14（约1.11倍），Pearson r = 0.98（p < 0.0001），Spearman ρ = 0.95（p < 0.0001），回归斜率0.97，R² = 0.96。Ugt49B1（胃盲囊27.27倍）和ABCC_MRP-PO（中肠25.06倍）的最强组织偏向性得到确认。

四、讨论与结论

讨论部分指出，意大利蝗解毒基因的组成差异不仅体现在总数上，更集中于特定解毒相关亚家族：CYP3（273 vs 101 vs 137）、CYP4（147 vs 45 vs 60）、α-酯酶（364 vs 80 vs 157）、Sigma-GST（50 vs 12 vs 13）和UGT305（171 vs 45 vs 75）。从生态适应角度，解毒基因家族的组成规模与昆虫取食范围广度密切相关。意大利蝗作为典型多食性蝗虫，其寄主植物涵盖多个科属，产生的次生代谢物化学结构高度多样，对解毒系统的底物广度提出更高要求。解毒相关亚家族在意大利蝗谱系内的组成比例偏移可能代表其适应多样化植物化学环境的重要遗传基础。

系统发育分析揭示了一致的组成模式：直接参与解毒功能的亚家族（CYP3、CYP4、α-酯酶、β-酯酶、Sigma-GST、UGT305、ABCB、ABCC）在意大利蝗中成员数较多、比例偏高且多为物种特异性簇；而承担基础生理功能的保守亚家族（Mito-P450、AChE、ABCE、ABCD、UGT36）在三个物种间成员数相似、分支紧凑、组成比例稳定。α-酯酶在COE家族中的极高比例（69.7%）是意大利蝗最显著的组成特征之一，可能反映了长期适应寄主植物多样化化学防御过程中通过基因复制和功能分化积累了广谱酯水解能力。Sigma-GST亚家族比例的显著升高可能增强其II相缀合解毒能力。

组织表达谱揭示了意大利蝗解毒代谢系统的空间功能分工特征：胃盲囊是II相缀合反应（GST、UGT）和I相水解代谢（COE）的核心场所，马氏管在I相氧化代谢（P450）和III相外排（ABC）中发挥主导作用。单基因水平上，GstT2（马氏管）、GST-se（中肠）、Ugt35E2（马氏管）、CYP3_Cyp6g1（马氏管）和CG1128aE9_alpha_esterase（胃盲囊）等基因表现出显著组织特异性。ABC转运蛋白的差异表达模式进一步证实了组织间功能分化。

RT-qPCR独立验证了转录组推断的组织偏向性表达的方向和幅度，两种方法高度一致，支持了基于FPKM的表达谱分析的可靠性，并确认了本研究鉴定的组织偏向候选基因可作为后续功能研究的稳健靶标。

研究结论：本研究较为完整地解析了意大利蝗解毒代谢基因的家族组成，共鉴定1602个解毒代谢基因。五个家族中构成其外源物代谢主要基础的解毒亚家族为CYP3、CYP4、α-酯酶、β-酯酶、Sigma、Delta、UGT305、UGT50、ABCC和ABCB。组织表达谱阐明了四个关键解毒组织的功能分工：胃盲囊是II相缀合反应（GST、UGT）和I相水解代谢（COE）的核心场所，马氏管在I相氧化代谢（P450）和III相外排（ABC）中发挥主导作用。RT-qPCR验证的11个组织偏向基因全部重现转录组方向（Pearson r = 0.98），确认了表达谱分析的可靠性，并确立了这些基因作为后续RNAi和杀虫剂诱导功能研究的稳健候选。本研究为阐明意大利蝗环境适应的分子基础提供了重要资源，也为蝗虫解毒代谢系统的比较进化研究提供了新的数据支持。