《Hydrobiologia》：Shotgun metagenomics uncovers novel features of the cyanobacterial community in a large, shallow lake in Central Italy
编辑推荐：
本研究旨在通过鸟枪法宏基因组学分析，解析意大利中部一座大型浅水中营养—富营养湖泊——特拉西梅诺湖（Lake Trasimeno）中蓝细菌群落的分类结构与功能潜力。样品的宏基因组谱分析不仅表征了细菌群落结构，还识别出若干稀有物种；这些物种在后续组装与分箱过程中基
本研究旨在通过鸟枪法宏基因组学分析，解析意大利中部一座大型浅水中营养—富营养湖泊——特拉西梅诺湖（Lake Trasimeno）中蓝细菌群落的分类结构与功能潜力。样品的宏基因组谱分析不仅表征了细菌群落结构，还识别出若干稀有物种；这些物种在后续组装与分箱过程中基于宏基因组组装基因组（metagenome-assembled genomes, MAGs）的解析中未被检出。MAGs的分类鉴定基于基因组相似性指标的计算，即平均核苷酸一致性（average nucleotide identity, ANI）和平均氨基酸一致性（average amino acid identity, AAI），并结合系统基因组学分析，从而证实了监测调查中既往已记录的多种蓝细菌的存在，同时鉴定出地中海地区新记录类群，如荷兰原绿丝藻（Prochlorothrix hollandica）及一个尚未表征的原绿丝藻属（Prochlorothrix）类群。对原绿丝藻MAGs的基因组注释显示，其特征与产氧型、非固氮型、产生叶绿素b（chlorophyll-b, Chl-b）的光合蓝细菌一致。对蓝细菌MAGs的注释未发现编码蓝细菌毒素的操纵子。在保守分析策略下，抗菌药物抗性基因（antimicrobial resistance genes, ARGs）仅在霍乱弧菌（Vibrio cholerae，Enterobacterales）以及少数未分箱重叠群（unbinned contigs）中被识别，进一步表明ARGs在淡水环境中传播具有全球性重要意义。
该文发表于《Hydrobiologia》，围绕意大利中部特拉西梅诺湖（Lake Trasimeno）夏季蓝细菌群落的分类组成、功能潜力以及潜在健康风险展开系统研究。研究背景在于，传统蓝细菌分类依赖形态学、分子数据及生态生理特征相结合的多相分类策略，虽然在已分离培养藻株的鉴定中具有较高准确性，但其前提是必须获得纯培养或至少可分离材料，这使大量不可培养或难培养的蓝细菌难以被纳入研究。与此同时，显微镜方法对形态简单、缺乏足够鉴别性状的类群分辨能力有限，尤其难以识别微小蓝细菌及若干丝状类群，因此难以完整揭示天然水体中蓝细菌的真实多样性及其功能特征。近年来，高通量测序技术（HTS）推动环境DNA研究快速发展，但扩增子测序通常受限于标记基因片段较短、分类分辨率有限及功能推断能力不足。相较之下，鸟枪法宏基因组学可直接从环境样品中重建基因组草图，并结合基因注释识别代谢功能、次级代谢生物合成基因簇（biosynthetic gene clusters, BGCs）以及抗菌药物抗性基因，为蓝细菌生态学、毒素风险评估及微生物群落功能解析提供了更高分辨率的工具。
在这一背景下，研究人员选取特拉西梅诺湖两个滨岸站位的夏季环境样品开展鸟枪法宏基因组分析。该湖为大型浅水中营养—富营养湖泊，长期受蓝细菌优势影响，尤其以拉氏尖孢蓝细菌（Raphidiopsis raciborskii）占优。研究的核心目标包括：解析蓝细菌群落及主要细菌类群的分类组成；对蓝细菌MAGs进行分类与功能注释，重点评估蓝细菌毒素相关BGCs；对地中海区域新检出的两类原绿丝藻属类群开展系统基因组学、功能基因组学与泛基因组学比较；同时筛查细菌与蓝细菌基因组及重叠群中的ARGs。研究最终表明，特拉西梅诺湖蓝细菌群落不仅包含既往监测已知类群，还存在此前未在地中海地区记录的荷兰原绿丝藻及一个未定原绿丝藻类群；这些原绿丝藻基因组具有典型产氧光合、非固氮、含叶绿素b的特征。更重要的是，所有蓝细菌MAGs中均未发现编码经典蓝细菌毒素的操纵子，而ARGs仅在霍乱弧菌及少量未分箱重叠群中检出。该研究的意义在于，证明了宏基因组学能够弥补传统分类方法对隐蔽类群识别不足的问题，同时为淡水生态系统中蓝细菌毒性风险与抗性基因传播提供了更加客观和可重复的评估框架。
方法上，研究人员于2023年9月在特拉西梅诺湖Castiglione del Lago（CAST）与San Feliciano（SFEL）两个滨岸站位采集表层水样，并结合地方环境机构前一日采集的湖心理化数据。样品经0.2 μm滤膜过滤后提取环境DNA，采用Illumina NovaSeq 6000进行150 bp双端测序。生物信息学分析包括：读段质控、去除人源污染序列、基于Sylph的k-mer宏基因组谱分析；随后进行组装、分箱与MAGs构建，并借助GTDB-Tk、ANI、AAI和系统基因组树完成分类鉴定。针对原绿丝藻属进一步开展泛基因组学分析、KEGG代谢注释、antiSMASH次级代谢基因簇预测和AMRFinderPlus抗性基因筛查。样本来源明确为特拉西梅诺湖两个夏季滨岸环境样品。
研究结果部分首先报告了环境背景与毒素检测结果。采样期间两站位表层水温约24.6–25.0 °C，透明度仅0.3 m；高可溶性反应磷、总磷、叶绿素a与高pH、溶解氧饱和度共同表明湖泊处于中营养—富营养状态，并存在潜在氮限制。液相色谱-质谱联用（LC–MS）分析未在全部环境样品与网采浓缩水样中检测到可测定浓度的蓝细菌毒素。这一结果为后续基因组层面的毒素操纵子缺失提供了化学分析支撑。
在“Light microscopy”部分，显微观察显示两个滨岸站位的浮游植物群落均由R. raciborskii绝对优势主导，丰度很高。另有少量Anabaenopsis、Pseudanabaena、Planktothrix、Merismopedia及若干未定丝状蓝细菌，但由于高密度优势种造成的拥挤效应，稀有类群在显微镜下难以准确辨认。这一结果提示，仅靠显微形态学难以完整解析群落结构。
在“Read-based metagenomic profiling and identification of MAGs”部分，研究人员利用基于k-mer包含关系的宏基因组谱分析揭示，两站位细菌群落结构高度相似，主要类群包括放线菌门（Actinobacteria）、γ-变形菌纲（Gammaproteobacteria）、拟杆菌纲（Bacteroidia）和α-变形菌纲（Alphaproteobacteria）等。蓝细菌中，Raphidiopsis / R. raciborskii占绝对优势，其后依次为Nodosilinea、Prochlorothrix、Pseudanabaena、Anabaenopsis arnoldii、Sphaerospermopsis kisseleviana及Cyanobium等。谱分析还识别出部分稀有或非光合蓝细菌类群，如Vampirovibrionia。进一步的组装与分箱获得了一系列蓝细菌MAGs，且其覆盖度与Sylph估计的相对丰度高度吻合，说明宏基因组谱分析与MAG重建结果具有较强一致性。特别重要的是，在两个站位中均检测到两种不同的原绿丝藻类群，即P. hollandica和一个未定原绿丝藻属类群。
在“Taxonomy and main functional characteristics of the Lake Trasimeno Prochlorothrix assemblies”部分，研究人员详细解析了4个原绿丝藻MAGs。P. hollandica与未定原绿丝藻基因组大小约4.3–5.0 Mbp，完整度高、污染率低。标记基因分析支持P. hollandica的种水平鉴定，而未定类群与北美和亚洲若干环境序列仅呈中等一致性，提示其分类地位尚未完全澄清。代谢注释显示，这些基因组均包含产氧光合作用、Calvin循环、三羧酸循环与糖酵解等核心代谢模块，并具备叶绿素b合成关键酶叶绿素化物a加氧酶（chlorophyllide a oxygenase, CAO）以及叶绿素a/b结合捕光蛋白相关基因，符合原绿丝藻的特殊光合特征。与此同时，基因组中缺乏完整藻胆体（phycobilisome）生物合成基因簇，也未检测到固氮相关基因，说明其属于非固氮型光合蓝细菌。研究还识别出硝酸盐/亚硝酸盐转运、磷酸盐与膦酸盐转运、渗透保护物质转运及气泡囊泡（gas vesicle）相关基因，提示其具备适应富离子浅水环境和调节浮力的能力。
在“Comparison of the main functional characteristics with other Prochlorothrix assemblies analyzed in this work”部分，作者将湖泊中的原绿丝藻MAGs与公共数据库中的其他原绿丝藻基因组进行比较，发现P. hollandica和未定原绿丝藻在主要KEGG模块层面差异不大，均具有CAO和叶绿素a/b结合光捕获蛋白基因，但均缺乏完整藻胆体和固氮基因复合体。这表明这些类群在基本光合与氮代谢策略上具有共同性。
在“Phylogenomic and pangenomic analyses of Prochlorothrix”部分，系统基因组树明确将P. hollandica与未定原绿丝藻分为两个主分支，分支间ANI和AAI值较低，支持二者之间存在清晰的进化分化；而未定原绿丝藻主支内部分化较大，存在多个地理相关亚分支。泛基因组分析共识别出8772个基因簇，其中单拷贝核心基因簇占比较低，单例基因比例较高。基于Heaps定律的稀释分析表明，原绿丝藻属泛基因组属于开放型（open pangenome），意味着随着更多基因组加入，仍将不断发现新的基因家族。核心基因主要与基础代谢、翻译、能量生产和细胞结构维持相关，而壳层和云层基因则富集于信号转导、防御机制、可移动遗传元件及噬菌体相关功能，反映出该类群在不同生境中通过基因获得和基因组可塑性实现适应分化。
在“Secondary metabolites”部分，研究发现特拉西梅诺湖4个原绿丝藻基因组仅含少数常见次级代谢基因簇类型，如萜类（terpene）、NRPS/NRPS-like等，但与已知特定化合物基因簇无显著对应关系。扩展至全体蓝细菌MAGs后，虽然识别出多种BGC类型，并在部分Nostocales类群中检测到异形胞糖脂、hassallidin、mycosporine等相关簇，但所有蓝细菌基因组均未发现编码经典蓝细菌毒素的完整操纵子。结合LC–MS结果，表明该湖当次采样中未显示经典蓝细菌毒素生产潜力。
在“Antimicrobial resistance”部分，所有蓝细菌MAGs及全部原绿丝藻基因组均未检出ARGs。相反，在CAST站位的霍乱弧菌MAG中高置信度识别出与黏菌素（colistin）耐受相关的almEFG连续基因簇；对全部重叠群进一步分析后，还在未分箱序列中检出林可酰胺类、β-内酰胺类及汞抗性相关基因。这说明淡水环境中的抗性风险更多体现在伴生异养细菌而非蓝细菌本身。
讨论部分强调，宏基因组谱分析能够在无需培养的条件下高效识别稀有和难辨认类群，并补充MAG方法对低丰度类群重建不足的局限。研究证实，特拉西梅诺湖细菌群落呈现出适应高有机质、浅水、盐度波动环境的结构特征，而蓝细菌群落则远比显微观察揭示的更加复杂。对于原绿丝藻属，研究不仅首次在地中海区域获得其明确记录，还表明未定原绿丝藻类群具有较强的系统发育异质性和地理分散特征。功能上，P. hollandica和未定原绿丝藻均表现为含Chl-b、非固氮、产氧光合蓝细菌，但在气泡囊泡基因、信号转导、膜转运、噬菌体防御和可移动元件等方面存在差异，提示其生态适应策略并不完全相同。开放型泛基因组进一步说明该属仍具有较大的遗传多样性等待发掘。就应用层面而言，研究指出，宏基因组学对于娱乐性水体中潜在有毒蓝细菌及抗性基因风险的监测具有直接价值，但由于本研究仅基于夏季单次采样，尚不能代表全年动态。
研究结论可概括为：宏基因组分析在确认特拉西梅诺湖既有蓝细菌种类的同时，识别出此前未在地中海地区记录的P. hollandica和一个未表征的原绿丝藻类群；P. hollandica与未定原绿丝藻之间的进化关系仍需更多基因组与可能的表型证据加以厘清。全部蓝细菌MAGs中均未发现编码蓝细菌毒素的操纵子，这为环境监测提供了重要背景信息。基于保守标准，所有蓝细菌基因组中均未检出ARGs，同样具有重要意义；但霍乱弧菌及未分箱重叠群中检出的多类ARGs提示，淡水环境中的抗菌药物抗性传播仍需被谨慎评估。作者同时指出，未来仍有必要结合其他季节样品持续更新对该湖微生物群落结构与功能的认识。
