黑丽蝇不同蛹内发育阶段的基因调控机制的综合形态学和转录组学分析

《Insects》：Integrated Morphological and Transcriptomic Analyses of Gene Regulatory Mechanisms in Different Intra-Puparial Developmental Stages of Phormia regina

　　

黑丽蝇（Phormia regina (Meigen, 1826) (双翅目: 丽蝇科)）是一种具有重要法医学价值的食尸性丽蝇，在死后间隔（PMI）推断中发挥关键作用。该物种的蛹内时期约占其整个未成熟发育时长的一半。因此，精确确定蛹内样本的日龄对于提高PMI推断的准确性至关重要。本研究在恒温25?°C下进行，检查了不同发育阶段的蛹内形态学、组织切片和体重。此外，通过结合PacBio单分子实时（SMRT）测序和Illumina短读长下一代测序（NGS）的混合测序策略，进行了全长转录组分析。形态学研究揭示了蛹内形态特征、组织和器官的重建过程以及蛹期内体重的持续下降。转录组研究共鉴定出425,349条全长非嵌合（FLNC）序列。去除冗余后，获得84,852条转录本序列，其中46,325条转录本获得注释。以第0日蛹（D0）为对照，对后续发育阶段（D1–D5）的样本进行差异表达基因（DEG）分析，分别鉴定出4242、7964、9509、10,526和10,011个DEG。GO富集和KEGG通路分析结果为不同发育阶段的行为特征提供了合理解释。筛选出8个蛹内发育候选基因，其中大部分在蛹内时期的不同时间点高表达。实时荧光定量PCR（qRT-PCR）结果显示，基因表达趋势与RNA-seq数据一致，证实了RNA-seq结果的可靠性。本研究为法医学PMI推断中黑丽蝇的关键形态学和分子标记提供了依据。

### 论文解读：黑丽蝇蛹内发育的形态与转录组整合分析

#### 研究背景与问题

黑丽蝇（*Phormia regina*）是法医学中用于死后间隔（PMI）推断的重要指示昆虫，其蛹内发育期占整个未成熟期的近一半时间。然而，蛹内外部形态在化蛹数小时后即不再显著变化，仅凭外观难以精确区分后续发育阶段；传统形态学观察虽直观，但易受研究者主观判断影响。组织切片虽可揭示内部重建过程，但单独使用仍存在分辨率不足。转录组分析可提供高客观性的分子数据，但此前黑丽蝇缺乏覆盖连续蛹内发育阶段的全长转录组资源。因此，整合形态学与转录组分析，建立多指标年龄推断体系，是提升PMI准确性的关键。

#### 研究内容与结论

研究人员在恒温25?°C下，对黑丽蝇蛹内发育的6个时间点（D0–D5）系统开展了形态学观察、组织学切片、体重测量，以及结合PacBio单分子实时（SMRT）测序与Illumina短读长下一代测序（NGS）的混合全长转录组测序。主要结论包括：成功构建了内部形态与组织重建的详细图谱，发现体重在蛹内期间持续下降，尤以D0–D1阶段最为显著（占总减重的约70%）；全长转录组获得84,852条非冗余转录本，其中46,325条获得功能注释；以D0为对照，后续各阶段鉴定出4242至10,526个差异表达基因（DEG）；GO和KEGG富集分析揭示不同阶段分别富集于氨基酸代谢、黑色素合成、氧化磷酸化、三羧酸（TCA）循环等通路；时间序列聚类识别出持续上调（Cluster 9，富集于肌肉与能量代谢）和持续下调（Cluster 3，富集于自噬与蛋白降解）的基因模块；8个候选基因（如肌球蛋白调节轻链2、几丁质酶10、视蛋白Rh1等）的qRT-PCR验证与RNA-seq趋势一致。该研究发表在《Insects》期刊，为黑丽蝇蛹内日龄的高精度推断提供了关键形态学与分子标记，显著提升了法医学PMI估计的可靠性和分辨率。

#### 主要关键技术与方法

样本来源：黑丽蝇原始采集于中国宁夏石嘴山（38°98′ N, 106°52′ E），经形态鉴定与分子验证后，在实验室饲养建立稳定种群，实验所用个体均为同一亲代群体的后代。关键技术方法包括：（1）**形态学解剖与组织学切片**：每24小时取样8个蛹，经热水固定、乙醇保存后，在体视显微镜下剥除蛹壳观察内部结构；并行石蜡包埋、HE染色制备5–10 μm切片，记录组织重建进程。（2）**体重批量测量**：每24小时对77个活蛹进行三次独立称重取平均值，反映群体代谢特征。（3）**全长转录组测序**：将各发育阶段（D0–D5）的RNA样本等量混合，构建PacBio SMRT全长cDNA文库（插入1–6 kb）；同时每个时间点3个生物学重复独立构建Illumina文库（插入200–300 bp，PE150测序），每个样本平均6 Gb clean data。（4）**差异表达与时间序列分析**：以PacBio非冗余转录本为参考，使用STAR比对、Kallisto定量、DESeq2筛选DEG（|log₂FC|≥1，校正P≤0.05）；Mfuzz模糊聚类（k=9）分析18,637个独特DEG的时间表达模式。

#### 主要研究结果

**3.1 体重**：77个活蛹的总体重在蛹内发育期呈现持续下降趋势。D0至D1下降最剧烈（从3.2650 g降至2.9337 g，降幅0.3313 g），占D0–D5总减重（0.4742 g）的约70%；此后下降速度减缓，D5时达2.7908 g，反映出早期广泛组织重塑的高能耗与后期成熟阶段的稳定代谢。

**3.2 蛹内形态特征**：恒温25?°C下，D0蛹壳呈浅黄至浅褐色，内部组织与蛹壳紧密粘连。D1出现头、胸、腹分化但界限不清，可见粗厚的腿和翅膀。D2头胸腹分界清晰，触角、口器可见，腿翅变薄。D3翅膀完全折叠，胸部出现白色鬃毛，腹部可见分节。D4复眼呈紫红色，头部刚毛、上唇、胸鬃毛、翅膀及腿缘呈棕色。D5触角灰黑色，复眼棕红色，上唇深棕色，翅膀深灰色，腿完全变黑，为羽化做好准备。

**3.3 组织切片**：HE染色切片系统揭示了蛹内组织重建过程。D0细胞密集均质分布，可见幼虫头咽骨骼和盘绕的中肠。D1头胸腹开始分化但界限不清，中肠膨大成囊状占据体腔中央。D2头胸腹分界清晰，中肠局部收缩分化为后肠原基，出现初级肌纤维束。D3背部肌群呈现横纹结构，中肠前部管腔收缩，头部和胸部咽肌出现。D4咽扩张肌形成完整肌束，腹直肠囊发育至最大容积，肌纤维直径显著增大、深嗜伊红染色。D5肌纤维束成熟、横纹清晰致密，中肠转化为成虫型，体腔形成完整功能分区。

**3.4 转录组测序结果与数据质量**：18个RNA样本RIN值均>7.0。PacBio测序产生47.58 Gb数据，获得537,183个环形一致性（CCS）读段，其中425,349个为非嵌合全长（FLNC）序列，聚类后获得129,155个一致性序列，冗余去除后得到84,852个转录本。鉴定出38,379个开放阅读框（ORF），27,168个简单重复序列（SSR），50,283个长链非编码RNA（lncRNA）。46,325个转录本获功能注释。Illumina数据Q30 > 85%。

**3.5 主成分分析（PCA）**：PCA显示6个发育阶段样本在降维空间中明显分离，Dim1和Dim2分别解释总变异的21.2%和11.5%；Pearson相关系数热图证实各阶段内生物学重复一致性高。

**3.6 差异基因表达分析**：以D0为对照，D1–D5分别鉴定出4242、7964、9509、10,526和10,011个DEG，其中上调基因数逐渐增加（2092→6502），下调基因数先增后稳（2150→3509）。Venn图显示433个共有上调基因和647个共有下调基因。共有上调基因GO富集于结合、催化活性、膜系统、代谢过程及行为节律调控；共有下调基因富集于合酶/裂解酶活性、内质网/高尔基体等内膜系统及核苷酸/氨基酸合成代谢。各比较组GO分析表明：D1阶段富集于芳香族氨基酸代谢、血淋巴凝固；D2富集于IMP从头合成、黑色素合成；D3富集于氧化应激响应、肌肉系统过程；D4富集于ATP合成耦合质子转运、肌肉收缩调控；D5富集于TCA循环、肌肉收缩。KEGG分析显示，酪氨酸代谢、光转导通路在各阶段均富集，且氧化磷酸化、TCA循环等能量代谢通路在D4–D5显著增强。

**3.7 时间序列分析**：对18,637个独特DEG进行Mfuzz聚类（k=9），获得9个具有不同时间表达模式的簇。研究人员选取持续下调的Cluster 3和持续上调的Cluster 9进行富集分析：Cluster 3核心基因富集于泛素介导的蛋白水解、细胞蛋白分解、蛋白酶体/溶酶体等通路，提示早期幼虫组织降解；Cluster 9核心基因富集于线粒体能量代谢、TCA循环、肌肉细胞分化、氧化磷酸化等通路，反映晚期成虫肌肉形成与能量供应。

**3.8 基因表达水平验证**：选取8个与视觉和肌肉发育相关的候选基因（如视蛋白Rh1、肌球蛋白调节轻链2、神经钙黏蛋白等），以RPL23为内参，通过qRT-PCR检测D0–D5表达趋势，结果与RNA-seq数据高度一致，证实测序数据可靠。

#### 讨论总结与研究结论翻译

**讨论**：研究人员指出，蛹内体重快速下降与早期广泛幼虫组织溶解及器官重塑所致的高能耗吻合，此后下降减缓对应于主要组织重建完成和肌肉/内分泌系统成熟。形态学特征中，复眼最先着色（D4），随后刚毛、鬃毛及腿缘着色，最后触角、翅等变为深色，其顺序与已报道的丽蝇科物种一致，构成有价值的表型标记。但色素沉着速率受温度和地理种群影响，需在不同温度及种群中验证。转录组层次上，DEG数量随发育持续增加后趋于稳定；时间序列中持续上调的骨骼肌蛋白基因表达与间接飞行肌（IFM）发育进程吻合，持续下调的自噬相关基因体现从组织降解向重建的转变。视蛋白Rh1在D4–D5显著上调，与复眼着色同步；形态发生相关激活子在D3达峰值，对应翅形态发生关键期。这些稳定的表达趋势为区分早、中、晚期蛹提供了候选分子标记。qRT-PCR验证了数据可靠性。形态学与分子标记的整合可显著提高PMI推断的分辨率。

**研究结论**（翻译）：本研究通过在25?°C下整合形态学与全长转录组分析，系统探究了黑丽蝇蛹内的动态发育过程。在形态学层面，成功构建了详细的内部与组织发育图谱；同时观察到蛹内体重呈整体下降趋势，以化蛹初期下降最为显著。在转录组层面，首次将第二代与第三代测序技术联合应用于黑丽蝇蛹内6个发育阶段（0–5天），不仅鉴定出相邻阶段间的DEG，还阐明了骨骼肌蛋白等关键基因的表达模式，qRT-PCR结果进一步证实了测序数据的可靠性。综上所述，本研究清晰阐明了黑丽蝇蛹内发育的形态特征和分子调控机制，为法医学领域的PMI推断和昆虫证据分析提供了更精确的科学依据。然而，所有实验在恒温25?°C下完成，实际犯罪现场温度波动，因此本数据仅作为基线参考，需在变温条件下检验；本研究仅使用单一实验室种群，候选标记的跨种群稳健性需进一步验证；筛选出的DEG的体内外调控机制尚未经功能实验验证；所提出的形态与分子标记应视为候选标记，其预测精度尚未在独立数据集或真实场景中测试。未来研究将结合比色、形态测量等定量成像方法，并在变温及不同种群中验证标记，结合多组学数据与机器学习算法，构建高精度年龄预测模型，推动法医昆虫学从描述性研究迈向预测性研究。