研究样本由重庆师范大学媒介昆虫重点实验室从重庆市沙坪坝区歌乐山国家森林公园垃圾堆放点采集，后续已完成前期驯化、饲养多代呈现稳定性状的蝇虫。

将昆虫经实验室饲养，所得到的卵(100粒)、幼虫(1龄、2龄、3龄各50条)、蛹(50只)、成虫(雌、雄各25只)分别制作成标本，采用KEYENCE-VHX-5000超景拍照系统(日本，基恩士)拍照，并采用OLYMPUS SZ61昆虫解剖镜进行观察。

根据通用型柱式基因组提取试剂盒说明提取DNA，将所得的DNA样本放在-20 ℃保存，以备后续使用。

PCR扩增引物信息如表 1所示。参照表 1文献中的引物序列合成引物，扩增COI、16S rDNA基因序列。

PCR反应体系(50 μL)：模板DNA 1 μL，正反引物(10 μmol/L)各2 μL，2X San Taq PCR Mix 25 μL，Sterilized ddH₂O补充到50 μL。PCR反应程序：95 ℃预变性5 min；95 ℃变性20 s，50 ℃退火30 s，72 ℃延伸30 s；35个循环后，72 ℃延伸10 min，于4 ℃保存。PCR扩增完成后，取1 μL扩增产物进行2%琼脂糖凝胶电泳检测，将凝胶板置于紫外灯下观察，并利用凝胶成像系统进行拍照分析。扩增产物由北京擎科生物科技股份有限公司测序。

将获得的COI和16S rDNA碱基序列片段通过DNAMAN 9.0软件进行拼接处理，并各截取250 bp和278 bp长度的核酸序列。将获取的COI和16S rDNA碱基序列片段提交至BLAST平台(网址：https://blast.ncbi.nlm.nih.gov/Blast.cgi)，并上传至GenBank数据库完成登录注册，进而获得序列登录号。借助MEGA 11.0软件内置的ClustalW工具执行序列比对操作，采用Kimura 2-parameter模型来计算遗传距离，并通过Bootstrap检验方法进行验证，设置1 000次重复抽样^[26]。借助PhyloSuite V 1.2软件整合样品序列与下载序列，构建序列矩阵，通过MAFFT(Multiple Alignment using Fast Fourier Transform)序列比对工具进行矩阵序列的排列对齐，并结合BioEdit进行人工微调，随后采用最大似然法(Maximum Likelihood)和贝叶斯法(Bayesian Inference)构建系统发育树。在最大似然法构建系统发育树过程中，通过IQ-TREE筛选出最优模型TIM2+R5+F，并基于1 000次重复序列，采用逐步添加法计算分支支持值^[27]。对于贝叶斯法构建系统发育树，采用MrBayes软件筛选出最佳适配模型GTR+I+G4+F，设置双链并行运算，每条链独立生成随机树，共运行200万代，每100代进行抽样^[28]，最终确保平均标准偏差值(ASDSF)低于0.01的收敛标准。

利用配制的饲料吸引蝇虫进行产卵。实验室培养的卵为乳白色，呈椭圆形米粒状，100粒卵的长度均在1~1.5 mm范围，宽度在0.25~0.35 mm范围(图 1)。在实验室环境中，当温度控制在30 ℃，光照条件为12 L(光亮)∶12 D(黑暗)，湿度保持在70%时，卵期约13 h。卵孔位于卵一端，呈圆形。

幼虫呈蠕虫状，体色介于淡黄色与乳白色之间，体表光滑，身体近透明，消化道清晰可见，虫体有12节，前尖后粗。在第二节段和第三节段之间观察到前螺旋瓣，有10~12个螺旋瓣分支和孔，幼虫第三节的刺呈单尖排列(图 2a)。在腹面终末节的水平面上有两个棕色圆形区域是后气门。幼虫尾节的边缘被突起包围，称为结节。不同的幼虫期通过后气孔裂口数量来确定，分为3龄。

1龄期幼虫体壁极薄，通体透明，50条幼虫平均体长约为3 mm。头咽部骨架发育未完全，呈浅黄色半透明状，口钩细长且呈特征性“7”字形弯曲，基部狭长，无附属骨片结构。呼吸器官未发育，无前气门，仅见极小后气门单裂结构(图 2b)。在30 ℃，12 L∶12 D光照周期及70%湿度条件下，该龄期持续约23 h。

2龄期幼虫体壁增厚至半透明，50条幼虫体长均增至4~6 mm(图 2c)，头咽部骨架几乎完整呈不规则“H”形，在不完整的周膜上有色素沉着，口钩增粗并形成“C”形弯曲。呼吸器官开始分化，前气门出现淡黄色细长色素管(图 2d)。后气门扩大为双裂结构(图 2e)，龄期后期可见腹缘气门环开口。该龄期持续约53 h。

3龄幼虫进入预蛹阶段，体色转为黄白，体壁高度几丁质化接近不透明，长度达到最大值，50条幼虫平均体长约为9 mm(图 2f)，虫体前端较细，头咽部骨架结实，背角比腹角长，附口骨点状出现，口钩退化为倒“L”形，基部矩形扩展并出现腹角突(图 2g)。呼吸器官完全成熟，前气门具11个指状突起且色素沉积加深(图 2h)。后气门显著增宽呈三裂结构，下方伴生凹陷区(图 2i)。该龄期最长持续72 h，后期幼虫停止进食，呈现扩散行为并寻找化蛹场所。

该物种蛹属围蛹类型，体呈纺锤形，50只蛹长度在(8±1) mm范围，宽度在(3±0.5) mm范围。蛹的颜色随着蛹的发育逐渐加深，由乳白色(早期蛹)变为淡黄褐色，最后变为黑褐色(图 3a)。呼吸器官形成粗短的牛角状呼吸角，后气门孔缘呈卵形，包裹典型的三裂式气门结构(图 3b)。蛹期持续4 d，占整个未成熟发育阶段的50%，是大头金蝇完成变态发育的关键节点。

成虫属于中型类群，50只成虫体长均在9~10 mm范围，体呈金属绿色，具强光泽。头部可自由转动，口器为典型舔吸式，口器周围鬃毛丛发达。颊部呈现出杏黄至橙色的渐变色调，颊毛呈黄色。胸部具铜绿金属光泽，前盾片覆薄层灰白色毛，侧片覆黑毛，后小盾片退化。腹部呈蓝绿铜色，各背板(除第5节)具紫黑后缘带，第1腹板具黄毛，其余腹板及背板侧缘呈现出黄黑毛混生足基节紧密排列，后基节与触角密布刚毛。翅膜透明，翅脉及腋瓣呈棕色调(图 3c，雄蝇)。

性别二态性特征显著：雄性复眼具大小眼面分区(上部2/3大型眼面与下部1/3小型眼面分界清晰)(图 3d)，第5腹节后缘具裂口，尾毛显著长于雌性；雌性额区在眼前缘处内凹(图 3e，雌蝇)，复眼间距较宽(图 3f)，产卵器呈套筒状伸缩结构。

将实验室待鉴定蝇种的DNA提取后，使用核酸检测仪测得蝇种胸部DNA体积质量比为132.666 ng/μL，腿部DNA体积质量比为17.08 ng/μL。经PCR扩增，扩增产物通过琼脂凝胶电泳可见两条条带，分别为该物种的COI和16S rDNA扩增条带。

PCR扩增蝇种COI基因和16S rDNA基因经测序，获得COI基因片段序列为278 bp，碱基A占30.2%，碱基T占37.8%，碱基C占17.6%，碱基G占14.4%，其中A+T(68.0%)远高于C+G(32.0%)。所采集样本16S rDNA基因片段序列为250 bp，A、T、C、G这4种碱基的平均组成为：碱基A占38.8%，碱基T占43.6%，碱基C占6.0%，碱基G占11.6%，其中A+T(82.4%)远高于C+G(17.6%)(表 2)。

将所测得的大头金蝇COI基因和16S rDNA基因片段序列在GenBank上进行Blast比对，结果显示两个基因片段均与大头金蝇参考序列呈现出高度同源性(核苷酸一致性大于99%)。其中，COI基因遗传距离为0.000~0.000，16S rDNA基因遗传距离为0.000~0.004，均显著低于与丽蝇科其他物种的遗传距离(p＜0.01)。该序列已上传至GenBank数据库，登录号分别为COI：PQ680654和16S rDNA：PQ676744。

在美国国家生物技术信息中心(NCBI)上下载蝇科家蝇(Musca domestica)COI序列和16S rDNA序列作为外群，根据形态及基因遗传距离分析待鉴定物种为丽蝇科金蝇属的大头金蝇(Chrysomya megacephala)。在麻蝇科中下载棕尾别麻蝇(Sarcophaga peregrina)16S rDNA序列；在丽蝇科中分别下载与鉴定物种相似及相关种的COI序列和16S rDNA序列：绿蝇属的丝光绿蝇(Lucilia sericata)、叉叶绿蝇(Lucilia caesar)，阿丽蝇属的巨尾阿丽蝇(Aldrichina graham)，金蝇属的肥躯金蝇(Chrysomya pinguis)、裸金蝇属的绯颜裸金蝇(Chrysomya rufifacies)及6个大头金蝇(Chrysomya megacephala)数据(确保待鉴定物种序列的准确性)，采用PhyloSuite V 1.2软件对以上12个COI序列、11个16S rDNA序列(肥躯金蝇、绯颜裸金蝇无16S rDNA序列)进行单基因建树(图 4a、b，图 5a、b)和多基因建树(图 4c、图 5c)分析。IQ-TREE和MrBayes两种方法构建的系统发育树均显示：以家蝇(Musca domestica)为外群时，待鉴定样本与大头金蝇在3种基因树上均形成高置信度分支(COI：83%、76.5%，16S rDNA：80%、82.7%，多基因：88%、61.2%)，而与丽蝇科其他物种显著分离(图 4、图 5)。在丽蝇属中，丝光绿蝇(Lucilia sericata)与叉叶绿蝇(Lucilia caesar)在部分进化树中聚为一支，获得66%、89%等Bootstrap支持率，体现了较近的亲缘关系；金蝇属的大头金蝇(Chrysomya megacephala)多个序列形成明显聚类分支，其中待测序列PQ680654与大头金蝇序列分支支持率达83%，且在图 5带置信值的进化树中分支数值进一步验证了聚类可靠性；阿丽蝇属的巨尾阿丽蝇(Aldrichina graham)亦呈现出独立分支特征。进化树结果、分子系统学分析结果与BLAST比对结论一致，明确支持待测样本为大头金蝇的物种鉴定。

系统发育树的结果进一步证实了待鉴定蝇种与大头金蝇的亲缘关系最近，从分子进化角度为大头金蝇鉴定提供了有力的支持。通过对COI基因和16S rDNA基因的分析，不仅能够准确鉴定大头金蝇，还能够揭示其与其他蝇类的进化关系，为大头金蝇分类和系统发育研究提供了重要的分子依据。

本研究运用形态学与分子生物学相结合的手段对未知蝇种进行鉴定，明确了从室外诱捕经实验室多代培养的蝇种为大头金蝇Chrysomya megacephala。在鉴定大头金蝇时要基于形态学特征和多个基因分析，才能得到准确的结果。本研究从卵、幼虫(3个龄期)、蛹、成虫(雌、雄)各阶段分别对大头金蝇进行了形态特征描述，提供了各虫态的特征彩图。同时，扩增了COI基因、16S rDNA基因，发现基因序列片段小且基因两端序列较为保守是大头金蝇鉴定和系统进化分析常用的分子靶标。

研究证实，16S rDNA与COI基因片段的种内遗传距离(均值＜0.01)显著低于种间距离(＞0.03)，且两者分布区间无重叠，符合分子鉴定标记的有效性标准^[29-30]。以大头金蝇为例，其16S rDNA基因种内距离为0.000~0.004，最大种内距离小于与其他蝇种的种间距离；COI基因种内距离为0.000~0.000，与丽蝇科其他物种的种间距离形成显著差异。双基因标记体系(16S rDNA+COI)在蝇类种属鉴定中展现出了高效的鉴定能力^[31-33]。

本研究采用PhyloSuite V 1.2软件构建IQ-TREE和MrBayes系统发育树，通过统计距离法推断亲缘关系，分子鉴定结果与形态学特征高度一致。不同蝇种在系统发育树上形成显著分离的分支，证明该方法可实现嗜尸性蝇类的种级鉴定。实验获得的COI和16S rDNA序列与GenBank大头金蝇序列相似度达99.85%~100%，在系统发育树中稳定聚于同一分支，明确支持物种鉴定结论。该方法通过检测mtDNA中16S rDNA(250 bp)和COI(278 bp)短片段即可实现种属鉴别，具有操作简便、成本低廉的优势，对实验设备及操作人员要求较低，适合在基础实验室推广，为传统形态学鉴定提供了高效辅助手段^[34-36]。尽管小片段分子标记存在局限性，但其快速经济的特点在基层应用中具有重要价值。未来结合表皮碳氢化合物等生化特征的多维鉴定体系，有望突破近缘种鉴定难题^[37-39]，推动嗜尸性蝇类鉴定技术的精准化发展。

明确大头金蝇的形态及生物学特征，有助于开展后续实验。虽然已有关于大头金蝇鉴定及形态特征的报道，但反映该昆虫各阶段详细且清晰的彩图较少，本研究展示的大头金蝇不同虫态图片可增加对该虫的直观认识。基于IQ-TREE和MrBayes系统发育树的分子系统发育分析，同时结合形态特征进行鉴定和分类，为蝇种分类和鉴定提供了新的思路和方案。