《Environmental Microbiology Reports》：Metagenomic Bioprospecting of Bioremediation-Related Genes in Amazonian Aquatic Ecosystems
编辑推荐：
受污染淡水系统的环境监测仍缺乏针对生物修复导向筛选的精选基因组资源。本研究开发了一个包含来自80个微生物物种的546个基因的精选目录，这些基因与重金属抗性、烃类降解、解毒及相关功能相关，并将其应用于巴西亚马逊帕拉州四个市镇（阿尔塔米拉、布拉干萨、马拉巴和帕拉戈
受污染淡水系统的环境监测仍缺乏针对生物修复导向筛选的精选基因组资源。本研究开发了一个包含来自80个微生物物种的546个基因的精选目录，这些基因与重金属抗性、烃类降解、解毒及相关功能相关，并将其应用于巴西亚马逊帕拉州四个市镇（阿尔塔米拉、布拉干萨、马拉巴和帕拉戈米纳斯）河流宏基因组的功能注释。功能注释采用eggNOG-mapper进行，宏基因组读数则利用Bowtie2与该精选目录进行比对。在所有河流系统中，汞和铜相关基因是最常检测到的功能类别，这一模式可能既反映了宏基因组中的序列代表性，也反映了精选参考目录的不均衡组成。布拉干萨显示出最大的检测基因库，而多样性分析、杰卡德指数、UpSet图和基于Dunn检验的分析揭示了显著的空间异质性和局部功能富集位点。反复关联的细菌分类群包括不动杆菌属（Acinetobacter）、肠杆菌属（Enterobacter）、假单胞菌属（Pseudomonas）、寡养单胞菌属（Stenotrophomonas）、贪铜菌属（Cupriavidus）和鞘脂菌属（Sphingopyxis）。这些发现表明，该目录为环境筛查和生物勘探提供了有针对性的基因组参考，但宏基因组检测代表的是功能潜力，而非经实验证实的修复活性。
研究背景与科学问题
淡水生态系统正日益受到城市化、采矿、农业和工业活动产生的持久性污染物（如烃类、重金属、农药等）的影响。这些污染物会破坏水生群落结构并威胁水质安全。生物修复利用微生物的代谢与适应能力转化、固定或解毒污染物，宏基因组学的发展极大地提升了识别相关基因与通路的能力。然而，现有的通用基因组数据库虽广泛可用，但针对热带及生物多样性丰富地区（如亚马逊流域）的、整合了微生物分类群、功能注释及污染物相关基因的精选目录仍然十分有限。因此，亟需建立一套针对性强的基因组资源，用于环境分子监测与生物勘探。
研究设计与主要结论
为填补这一空白，研究人员开发了一个包含546个生物修复相关基因的精选目录，涵盖80个微生物物种，涉及重金属抗性、污染物降解、解毒、生物膜形成等功能。随后，将该目录应用于巴西亚马逊帕拉州四个市镇——阿尔塔米拉（Igarapé Altamira河）、布拉干萨（Cereja河）、马拉巴（Itacaiúnas河）和帕拉戈米纳斯（Uraim河）的河流宏基因组数据中。研究假设，不同地点的人类活动压力差异会导致生物修复功能基因分布的空间变异。结果表明，汞和铜相关基因是所有河流系统中最常检测到的功能类别，这既反映了实际宏基因组中的序列特征，也与目录本身不均衡的组成有关。其中，布拉干萨拥有最丰富的检测基因库，而杰卡德距离、UpSet图和Dunn检验等多种分析共同揭示了显著的空间异质性与局部功能富集位点。反复关联的细菌类群主要包括鲍曼不动杆菌（Acinetobacter baumannii）、路德维希肠杆菌（Enterobacter ludwigii）、铜绿假单胞菌（Pseudomonas aeruginosa）、嗜麦芽窄食单胞菌（Stenotrophomonas maltophilia）、耐金属贪铜菌（Cupriavidus metallidurans）等。该研究证实，亚马逊水生生态系统中存在着高度空间结构化的生物修复基因库，所开发的精选目录可作为环境分子监测与生物勘探的有效工具，但其检测结果反映的是功能潜力，并非经实验验证的修复活性。
关键技术方法概述
研究人员首先通过结构化文献筛选，从130个候选分类群中筛选出80个满足条件的微生物物种，其基因组来源于NCBI的RefSeq数据库。功能注释采用eggNOG-mapper（基于DIAMOND搜索模式）结合eggNOG v5数据库完成，并通过BLAST（基本局部比对搜索工具）验证，仅保留相似度≥95%的序列。最终构建了包含546个基因的精选目录。宏基因组DNA提取后经NovaSeq 6000平台进行150 bp双端测序，原始数据经FastQC评估质量、Trimmomatic过滤（去除短于75 bp及Phred质量值<20的序列），随后使用MEGAHIT进行组装。宏基因组读数通过Bowtie2与精选目录进行比对，允许种子区最多一个错配，比对结果以每百万映射读段计数（CPM）进行标准化处理。所有统计分析均在R环境中完成，包括Pearson相关性、Shannon多样性指数、Jaccard相异性及Dunn事后检验。
研究结果分述
3.1 精选生物修复基因目录
通过结构化文献筛选与基因组筛选流程，研究人员最终获得了一个包含80个微生物物种、546个生物修复相关基因的精选目录。目录涵盖15个功能类别，但分布极不均衡：汞相关基因占主导地位（225条，41.21%），其次为铜相关基因（115条，21.06%）和碲/亚碲酸盐相关基因（56条，10.26%），汞和铜两类合计占比达62.27%。其他类别如生物膜形成、砷抗性、有机过氧化物抗性等频率较低。这种不均衡性意味着下游解读中，汞和铜相关序列的主导地位可能同时反映了生物学真实模式和目录自身的偏差。
3.2 四条河流中生物修复基因的空间格局
热图分析显示，每条河流沿程五个采样点的基因丰度均呈现明显的空间结构化。阿尔塔米拉的下游位点（AP3–AP5）相对丰度较高；布拉干萨的BP2和BP4最为富集；马拉巴的MP1呈现出独特的丰度模式；帕拉戈米纳斯的PP5相对富集最强。这些结果表明，尽管生物修复相关基因在整个河段均有检出，但其相对丰度存在显著的局部变化，可能与人类活动压力的空间差异有关。
3.3 污染物相关基因类别的分布
在所有四个市镇中，汞和铜相关基因均为最常检出的功能类别。例如，布拉干萨的汞相关基因超过90条，铜相关基因亦居次席。砷、镉、铬/铬酸盐、酚类、有机过氧化物等基因的检出频率较低且在不同市镇间存在差异。研究人员强调，汞和铜基因的主导地位不应被直接解读为这两种金属是水中浓度最高的污染物，而是序列代表性与目录组成的综合体现。
3.4 微生物多样性、组成相似性与统计支持
Pearson相关性分析表明，同一河流内部分采样点之间的基因丰度谱较为相似，但不同河流间的相似性格局各异。Shannon多样性指数显示，各河流内部存在局部多样性和丰富度的差异，例如阿尔塔米拉的AP5、布拉干萨的BP5、马拉巴的MP1以及帕拉戈米纳斯的PP4–PP5分别表现出较高的功能多样性。Jaccard相异性分析进一步揭示，即使整体丰度模式相似的位点，其基因组成也可能存在较大差异。Dunn事后检验确认了显著的空间分化位点：阿尔塔米拉的下游位点（AP4–AP5）在TerC和Chrome_Resist基因上分化明显；布拉干萨的BP4–BP5在TerD、TehB和DUF446上呈基因依赖性分化；马拉巴的MP1在CutA1和DUF446上反复分化；帕拉戈米纳斯的PP5则在DUF386和TerD上表现出一致的分化模式。
3.5 共享与位点特异性生物修复基因
UpSet图分析表明，每条河流都存在一个有限的共有功能组分（核心基因），同时伴随大量的局部变异。布拉干萨的总基因数最大（213条），核心基因为16条；马拉巴的MP1拥有最多的位点特有基因（37条）；阿尔塔米拉和帕拉戈米纳斯的核心基因较少（分别为2条和5条），但各自也有位点特异性的基因富集（AP3和PP4）。这些结果突出了局部功能富集位点作为未来生物勘探和功能验证目标的价值。
3.6 关键微生物类群与生物勘探意义
反复检出的细菌类群包括鲍曼不动杆菌（A. baumannii）、路德维希肠杆菌（E. ludwigii）、铜绿假单胞菌（P. aeruginosa）、嗜麦芽窄食单胞菌（S. maltophilia）和耐金属贪铜菌（C. metallidurans），它们主要与汞和铜抗性功能相关，同时也涉及砷转化、镉/铅抗性、苯酚降解等。研究人员特别指出，这些检测结果仅代表功能潜力，不能等同于实际的修复活性，且部分菌种（如A. baumannii和P. aeruginosa）属于条件致病菌，在后续生物勘探中需谨慎评估生物安全问题。
结论讨论与总结
研究人员在讨论部分指出，本研究展示的是一种基于宏基因组序列关联的功能潜力推断，而非直接的污染物浓度测量或基因表达证据。所有空间格局均应理解为序列检测层面的关联，而非因果生态关系。未来研究应聚焦于实验验证、微生物分离培养和功能测试，以确认最具应用潜力的候选菌株与基因。结论部分明确总结道：亚马逊水生生态系统拥有一个空间结构化的生物修复基因库，其中汞和铜相关功能占据主导地位。通过将精选目录与宏基因组筛查相结合，研究人员在帕拉州的四条河流中识别出了共享的抗性组组分和位点特异性的热点区域。该精选目录可作为环境分子监测与生物勘探的有价值的基因组参考，但研究本身基于精选参考数据和宏基因组推断，并不证明实际的修复性能。
