《Frontiers in Microbiology》：Exploring the Tenebrio molitor gut microbiota response to LDPE and PET: putative genetic indicators and methodological insights
编辑推荐：
昆虫肠道微生物组被认为是有可能携带与塑料代谢相关酶活性的潜在储库。研究人员使用16S rRNA测序和鸟枪法宏基因组学（shotgun metagenomics），考察了膳食暴露于低密度聚乙烯（low-density polyethylene, LDPE）和聚对
昆虫肠道微生物组被认为是有可能携带与塑料代谢相关酶活性的潜在储库。研究人员使用16S rRNA测序和鸟枪法宏基因组学（shotgun metagenomics），考察了膳食暴露于低密度聚乙烯（low-density polyethylene, LDPE）和聚对苯二甲酸乙二醇酯（polyethylene terephthalate, PET）条件下黄粉甲（Tenebrio molitor）肠道微生物群的分类学及功能动态。在扩增子序列变体（Amplicon Sequence Variant, ASV）水平上检测到显著的群落组成偏移，摄食塑料组富集了参与外源毒物（xenobiotic）代谢的分类单元。预测的功能变化提示芳香族化合物加工及氧化还原稳态相关通路丰度改变。研究人员通过可重复的开源流程进行宏基因组组装与功能注释，发现若干推定蛋白与邻苯二甲酸酯双加氧酶（phthalate dioxygenases）、脲烷酶（urethanases）及聚羟基脂肪酸酯解聚酶（polyhydroxyalkanoate depolymerases）具远缘同源性。归属于肠球菌属（Enterococcus）的一个宏基因组组装基因组（Metagenome-Assembled Genome, MAG）包含了大部分回收到的蛋白编码序列。尽管基因水平比较未显示统计学显著差异，但基因集富集分析（Gene Set Enrichment Analysis, GSEA）突出显示了塑料暴露条件下ABC转运蛋白特征及应激响应ATP酶的上调。总体而言，本探索性研究揭示了T. molitor肠道内的微生物偏移及代谢潜能的推定遗传指示因子，并为未来探究微生物在塑料生物转化中作用的研究提供了可重复的分析框架。
论文解读：黄粉甲肠道微生物对LDPE与PET的响应及宏基因组方法学评估
本研究发表于《Frontiers in Microbiology》。塑料污染是当前严峻的环境问题，传统回收与处理方式存在局限，微生物降解成为有前景的补充策略。黄粉甲等昆虫被证实可摄入并部分矿化PE和PS，其肠道微生物被认为参与塑料降解，但以往研究多依赖16S rRNA扩增子测序，缺乏高分辨率的鸟枪法宏基因组（shotgun metagenomics）研究，且缺乏对宿主DNA污染等技术瓶颈的系统评估，也缺少在同一流程下对比脂肪链聚合物（LDPE）与含芳香族酯键聚合物（PET）的研究。为此，研究人员以黄粉甲幼虫为模型，设置对照（麦麸）、麦麸+粉碎LDPE、麦麸+粉碎PET三组饲料，饲养30天后取肠道进行16S rRNA测序与鸟枪法宏基因组测序，结合可重复的生物信息学流程探讨微生物分类与功能响应，并评估宿主背景对昆虫关联宏基因组的功能信号检测影响。
主要关键技术方法：
研究人员使用西班牙Insectalia农场提供的Tenebrio molitor幼虫，设对照组（C，仅麦麸）、LDPE+麦麸组（LD）、PET+麦麸组（PET），各100条幼虫，饲喂30天（n=4生物学重复/组，每重复取30条幼虫 pooled肠道）。DNA提取采用QIAamp PowerFecal Pro Kit并纯化。16S rRNA扩增子靶向V4区（引物515F-Y/806R-B），Illumina MiSeq PE300测序，Qiime2(DADA2去噪→ASV；SILVA 138数据库注释；α/β多样性、指示种分析、PICRUSt2功能预测、STAMP差异检验)。鸟枪法宏基因组文库用Illumina DNA Prep制备，NovaSeq PE150测序；原始读段先用Bowtie2比对Tenebrio molitor参考基因组(icTenMoli1.1)去除宿主序列，经nf-core/mag流程(v3.2.1)：Fastp质控→再次去除PhiX及人基因组→MEGAHIT共组装→MetaBAT2与MaxBin2分箱(binning)→DAS Tool精炼MAG→CheckM评估质量→CAT/BAT进行MAG分类学注释→Salmon估算MAG相对丰度。功能注释采用Prokka(EC号)、OmicBox(GO)、GhostKOALA(KO)、InterProScan(蛋白结构域)。差异基因/蛋白采用limma(voom转换)，过代表征分析(Over-Representation Analysis, ORA)与基因集富集分析(Gene Set Enrichment Analysis, GSEA, fgsea包)进行功能富集。推定塑料降解酶通过本地BLASTp比对自建PlasticsDBPlus库（整合PlasticDB、PAZy、PMBD）筛选（E-value≤1×10?5，identity≥30%，query coverage≥20%）。
研究结果
3.1 Physiological outcomes and survival
所有组幼虫体长与体重均增加，LD组存活率97%，PET组99%，对照组100%，表明添加LDPE或PET未产生明显致死效应，黄粉甲对含塑膳食具生理耐受性。
3.2 Microbial community and inferred functional shifts (16S rRNA)
稀释曲线趋于饱和；PERMANOVA显示饮食组间微生物群落组成存在显著差异(F=4.19, R2=0.48, p=0.001)，NMDS与层次聚类可观察到分组趋势。指示种分析在ASV水平鉴定到与塑料摄食相关的分类单元。PICRuSt2推断功能显示塑料组中外源毒物生物降解与代谢、芳香族化合物处理等通路预测丰度发生变化，提示潜在代谢潜能调整，但该结果为基于分类组成的推断而非直接宏基因组证据。
3.3 Metagenomic sequencing and assembly quality control
尽管原始测序深度达每样本6400万–9100万条读段，经质控及宿主移除后仅保留6.2%–8.3%读段（典型昆虫肠道宏基因组应≥40%–65%），高读段重复率（>50%）与宿主DNA污染致有效微生物数据严重损失。共组装N50仅1079 bp（预期>5 kb），≥10 kb的contig仅6条，最长contig 19.5 kb，表明组装高度碎片化。分箱得到一株中质量MAG（MetaBAT2Refined-0.2，88%完整度/0.28%污染）归为Enterococcus（44%支持），另一bin被判定为宿主Tenebrio molitor DNA予以丢弃。说明宿主DNA是本研究所用昆虫肠道宏基因组的重大技术限制因素。
3.4 Metagenome-wide functional analysis and gene detection
功能基因稀释曲线未达平台，暗示测序深度不足以完全捕获功能多样性。基于Bray–Curtis的PCoA与NMDS显示饮食条件有一定分离趋势（PERMANOVA R2=0.285, p=0.008；ANOSIM R=0.181, p=0.051=临界），但组内变异大，残余方差占71.47%。limma单个基因差异表达分析、LEfSe分析及ORA均未获显著富集条目（adj. p>0.05）。功能注释显示777个非冗余推定蛋白中81.1%获Pfam注释，64.3%获GO注释，61.8%获KO注释，30.6%获EC号。
3.5 Targeted functional annotation and degradation potential
BLASTp比对NCBI NR库显示777个推定蛋白中约87%（677个）最佳命中为Enterococcus属（主要为E. villorum），其余少量来自Listeria、Streptococcus、Microbacterium（含M. tenebrionis）等。比对自建PlasticsDBPlus库初筛得14个具塑料/外源毒物代谢相关酶同源性的推定蛋白（邻苯二甲酸酯双加氧酶、聚羟基脂肪酸酯解聚酶poly(3-hydroxyalkanoate) depolymerase、聚乙烯醇脱氢酶polyvinyl alcohol dehydrogenase、脲烷酶urethanase、己内酰胺环状二聚体水解酶NylA等），严格阈值（identity≥50%）筛选保留3个PHB/PHA解聚酶同源蛋白。GSEA分析发现塑料暴露组中ABC转运蛋白家族特征（ProSitePatterns PS00211）及应激响应ATP酶显著富集；核心富集基因含聚羟基脂肪酸酯解聚酶同源基因、邻苯二甲酸酯双加氧酶及还原酶、4,5-二羟基邻苯二甲酸脱羧酶、邻苯二甲酸二氢二醇脱氢酶（与GSEA富集的醛-醇脱氢酶家族对应）、聚乙烯醇脱氢酶及脲烷酶同源序列；同时GO条目"ATP结合(GO:0005524)"、"ATP水解活性(GO:0016887)"及AAA+ATP酶、BRCT结构域、分子伴侣（蛋白解聚 chaperone）在塑料组富集，提示底物转运、氧化应激及DNA修复/应激响应激活。
讨论与结论翻译总结
讨论指出16S rRNA揭示的群落组成改变不能直接等同于Shotgun宏基因组可检出的功能改变，宿主DNA污染与低有效深度致基因水平差异不显著，传统ORA与LEfSe失效，而GSEA能捕捉协同微效的通路级信号——这是昆虫关联宏基因组在低微生物回收率下的推荐统计策略。Enterococcus的主导地位反映其核心共生菌地位而非特异塑料降解适应，但检出与塑料/外源毒物代谢相关酶的远缘同源序列提示肠道具一定遗传潜能。研究受限于短期暴露、无转录/酶活验证、宿主去除不充分及读段重复率高，未来应优化宿主DNA去除、加深测序（≥10–20 Gb/样本）、延长暴露期并做生化验证。
研究结论（翻译）：
综合来看，研究人员通过对摄食LDPE或PET饲料的黄粉甲幼虫肠道微生物组进行整合分析，发现了可能与塑料代谢相关的 tentative功能特征，包括与外源毒物降解相关的推定基因（如邻苯二甲酸酯双加氧酶、脲烷酶及聚羟基脂肪酸酯解聚酶同源基因）。基因集富集分析进一步提示ABC型底物转运蛋白及氧化应激相关酶的可能上调。这些初步观察支持昆虫关联微生物组在暴露于难降解碳源时可表现一定代谢可塑性。但受宿主DNA污染、读段重复及测序深度限制导致基因组与功能注释分辨率不足，本研究应定位为探索性及假设生成型调查而非确定性功能表征。研究同时提供了一套可重复的生物信息学流程以供后续宏基因组研究参考。
