《Phycology》：An Atlas of Macroalgal Epiphytic Bacterial Genomes Reveals Broad Biosynthetic and Metabolic Potential
编辑推荐：
宏藻表面栖息着附生细菌群落，它们在宿主代谢、胁迫耐受、化学防御和生态互作中具有潜在作用，然而在主要宏藻谱系中，这些细菌的基因组解析多样性仍未得到充分表征。研究人员从108个最初为宿主基因组测序而生成的培养来源宏藻鸟枪法数据集，以及另外三个具有配对高精度长读长和
宏藻表面栖息着附生细菌群落，它们在宿主代谢、胁迫耐受、化学防御和生态互作中具有潜在作用，然而在主要宏藻谱系中，这些细菌的基因组解析多样性仍未得到充分表征。研究人员从108个最初为宿主基因组测序而生成的培养来源宏藻鸟枪法数据集，以及另外三个具有配对高精度长读长和短读长数据的宏藻物种中，重建了附生细菌宏基因组组装基因组（metagenome-assembled genomes，MAGs）。对匹配MAGs的基准测试显示，不同测序策略之间的组装大小和基因组质量相当，支持基于短读长从宿主相关宏藻数据集中恢复MAG。在所有数据集中，研究人员恢复了3272个初步MAG，并保留1915个非冗余、质量受控的MAG，其中1543个无法通过基于平均核苷酸一致性（average nucleotide identity，ANI）的GTDB-Tk被分配到物种水平分类单元。这些基因组以假单胞菌门（Pseudomonadota）、拟杆菌门（Bacteroidota）和浮霉菌门（Planctomycetota）为主，三者合计占该集合的88.3%，并包含许多未在物种水平上分配的基因组。功能分析鉴定出12950个生物合成基因簇，归为7985个基因簇家族，20个生物合成丰富的MAG，以及222个KEGG模块。该图谱为研究宏藻全生物体功能和微生物天然产物潜力提供了基因组解析资源。
宏藻表面附生细菌具有影响宿主代谢、胁迫耐受、化学防御和生态互作的潜在功能，然而跨主要宏藻谱系的基因组水平多样性尚未被系统揭示。已有研究多聚焦于少数代表性宿主，缺乏广泛的多物种比较。为此，研究人员利用公共数据库中大量原本用于宿主基因组测序的宏藻宏基因组数据，系统重建了附生细菌基因组，旨在构建宏藻附生细菌基因组图谱并评估其生物合成与代谢潜力。该研究整合了108个来自藻类培养物保藏中心（CCAP/CCMP）的宏藻样本Illumina短读长数据集，以及红藻Kappaphycus striatus、绿藻Okellya marina和褐藻Saccharina japonica三个物种的配对高精度长读长和短读长数据。通过组装、细菌序列分类、分箱、去冗余和质量控制，最终获得1915个非冗余高质量MAG。比较分析表明，短读长数据可在宿主DNA背景下可靠恢复细菌MAG。分类学分析显示这些MAG以假单胞菌门（Pseudomonadota）、拟杆菌门（Bacteroidota）和浮霉菌门（Planctomycetota）为主，三者占88.3%，且1543个MAG无法在物种水平上归属，揭示大量未表征的细菌多样性。功能分析共鉴定出12950个生物合成基因簇（BGCs），归为7985个基因簇家族（GCFs），其中萜烯类BGC最为丰富；另有20个生物合成丰富的MAG，主要分布于Myxococcota_A等谱系。KEGG模块分析识别出222个模块，以中心碳代谢、脂肪酸合成、核苷酸合成等保守功能为主，而DMSP降解途径等模块的完整性在不同MAG间呈现差异，提示功能分化。该图谱为宏藻全生物体功能研究和微生物天然产物挖掘提供了重要基因组资源。
研究人员采用的关键技术方法包括：对短读长数据使用fastp进行质控，并用MEGAHIT组装；长读长数据使用metaMDBG组装。利用Tiara对组装contigs进行分类，提取细菌序列；通过bwa mem将短读比对至细菌contigs，并使用MetaBAT2和VAMB分别进行分箱，DAS_Tool整合分箱结果，CheckM2评估基因组质量。MAGs用MAGmax以95%完整性、50%一致性、90%纯度标准去冗余。分类学注释采用GTDB-Tk（GTDB release 232），基因预测使用Prodigal，功能注释基于eggNOG-mapper，KEGG模块完整性由kegg-pathways-completeness-tool计算，生物合成基因簇用antiSMASH预测并用BiG-SLiCE聚类。
研究结果如下：
3.1 测序数据与宏藻附生细菌的基因组解析重建。该部分分析了108个宏藻培养样本的测序数据，建立了从序列组装、细菌contigs分类、宏基因组分箱到基因组质量评估的完整工作流程，为后续MAG恢复提供了方法基础。
3.2 短读与长读宏藻相关数据集中MAG重建的比较。针对三个同时具有长读和短读数据的宏藻物种，独立重建MAG并比较。结果显示，匹配MAG的组装大小在两种策略间高度线性相关（R2=0.94），而基因组质量分数无显著线性对应（R2=0.003），长读数据未系统性优于短读数据。这表明，在宿主DNA背景下，短读数据足以支撑可靠的MAG重建。
3.3 重建宏藻附生MAG的系统发育多样性与质量。从所有数据集共获得3272个初步MAG，经去冗余和质量过滤后保留1915个非冗余MAG。基于bac120标记蛋白的系统发育树显示这些MAG覆盖广泛细菌多样性，且多数具有高完整性和低污染。门水平分布表明，Pseudomonadota、Bacteroidota和Planctomycetota合计占88.3%，其中1543个MAG无法基于ANI在物种水平分类，说明宏藻表面存在大量参考数据库中尚未涵盖的基因组多样性。
3.4 宏藻附生MAG的生物合成基因簇与代谢模块谱。antiSMASH预测共鉴定出12950个BGCs，约20%位于contig边缘，产物类型以萜烯类最多，其次为核糖体合成和翻译后修饰肽（RiPP）、非核糖体肽合成酶（NRPS）和聚酮合酶（PKS）等。BiG-SLiCE聚类得到7985个GCFs（默认参数）和5409个GCFs（--complete参数）。以每个基因组含20个以上BGC为阈值，识别出20个生物合成丰富的MAG，主要集中在Myxococcota_A，另有Bacteroidota、Actinomycetota等。KEGG模块分析识别出222个模块，分布最广的是脂肪酸生物合成（M00082、M00083）、糖酵解和中心碳代谢（M00002、M00009等）、核苷酸生物合成（M00048等）、能量代谢（M00155、M00157）以及氨基酸和辅因子代谢（M00621等），而DMSP降解脱甲基途径（M00990）的完整性在不同MAG中变化较大，提示该代谢功能可能局限于特定附生细菌谱系。
讨论部分指出，深度测序的宏藻宿主基因组数据集是尚未充分利用的附生细菌基因组资源。参考基因组数据库正在快速扩展，但未培养的宏藻附生细菌仍有大量物种水平未归属的基因组。宏藻表面蕴含丰富的生物合成潜力，与珊瑚相关微生物组中发现的多样性相当。该研究局限性在于未在组装前去除宿主reads，可能导致低丰度细菌基因组恢复不完整，因此所得MAG集合可能低估了实际细菌多样性。研究结论：该图谱为研究宏藻全生物体功能和微生物天然产物潜力提供了基因组解析资源。
