坦桑尼亚小莫梅拉湖嗜极端细菌的分类多样性与功能预测

《Discover Life》：Taxonomic diversity and functional prediction of extremophilic bacteria of small Momela Lake, Tanzania

　　

苏打湖代表了地球上一类最极端的水生生态系统，然而其微生物多样性与功能潜力仍缺乏充分表征，尤其是在东非地区。本研究报道了对坦桑尼亚北部一处苏打湖——小莫梅拉湖——水体与沉积物中细菌群落开展的首次鸟枪法宏基因组（shotgun metagenomics）分析。样本采用Illumina NovaSeq测序，并结合MetaPhlAn4与HUMAnN3进行分析，以解析其分类组成、多样性格局及预测代谢功能。结果显示，群落中共检出27个不同物种，分别归属于Pseudomonadota、Bacteroidota、Firmicutes、Cyanobacteria和Verrucomicrobia，其中以参与硫循环的脱氮硫碱弧菌（Thioalkalivibrio denitrificans）及未培养的候选类群（Candidatus）谱系占优势。值得注意的是，基于当前可获得文献资料，研究检测到多个在全球苏打湖环境、非洲及东非苏打湖环境中首次报道的类群，从而拓展了嗜极端生物的已知生物地理分布。功能谱分析显示，核苷酸与氨基酸生物合成通路显著富集，同时还检测到与pH稳态、渗透胁迫耐受、硫代谢及碳固定相关的标志基因。这些特征提示该类群具有显著的代谢灵活性，从而支持微生物在极端碱性和高盐条件下持续存活。上述发现强调，小莫梅拉湖是新型功能性细菌资源的重要储库，在生物技术与生物勘探（bioprospecting）方面具有良好应用前景。

该论文发表于《Discover Life》，围绕坦桑尼亚阿鲁沙国家公园内小莫梅拉湖这一典型苏打湖生态系统，系统解析了其水体与沉积物细菌群落的分类多样性及潜在功能特征。研究背景在于，苏打湖属于高碱、高盐的极端水生环境，pH通常超过9，并富含碳酸钠等溶解盐类，仅有具特殊适应机制的微生物能够长期定殖和繁衍。尽管全球部分苏打湖已开展微生物生态研究，但坦桑尼亚苏打湖，尤其是小莫梅拉湖的微生物分类组成、遗传背景和代谢功能潜力仍缺乏系统性宏基因组证据。既往关于该湖的研究主要集中于可培养放线菌门（Actinomycetota）成员，尚未从群落整体层面揭示其细菌多样性及功能基因储备。因此，开展该研究的核心意义在于填补东非欠研究苏打湖在群落宏基因组学层面的基础数据空白，并为嗜碱、嗜盐微生物资源挖掘及生物技术开发提供依据。

研究人员以小莫梅拉湖岸线不同位置采集的水体与沉积物样本为对象，利用鸟枪法宏基因组测序从环境总DNA层面直接描绘微生物群落结构，进而结合分类学与功能学分析框架，评估该湖细菌群落的优势类群、群落多样性、代谢通路构成以及关键适应性标志基因。结果表明，小莫梅拉湖尽管处于极端理化环境下，仍蕴藏较丰富的细菌多样性。群落主要由Pseudomonadota、Bacteroidota、Firmicutes、Cyanobacteria和Verrucomicrobia构成，其中脱氮硫碱弧菌（Thioalkalivibrio denitrificans）和多个未培养候选谱系（Candidatus）占据优势。功能层面上，群落富集核苷酸与氨基酸合成相关通路，并携带与Na⁺/H⁺逆向转运、pH稳态维持、兼容溶质合成、硫氧化/硫还原及RuBisCO介导碳固定相关的标志基因，说明这些微生物可能通过较强的代谢可塑性适应强碱和高盐胁迫。论文的重要意义在于首次以宏基因组方式确认了该湖为功能性新型细菌资源库，并指出其在极端酶、代谢产物及生物勘探方面的潜在价值。

研究人员采用的关键技术方法主要包括以下几类：首先，在坦桑尼亚小莫梅拉湖岸线10个随机点采集三重复水体和沉积物样本，并同步原位测定温度、pH、盐度、电导率（EC）和总溶解固体（TDS）等理化参数；其次，分别提取水体与沉积物基因组DNA（gDNA），因经费限制按等比例混合形成水样复合样本SMW1和沉积物复合样本SMS1；随后使用Illumina NovaSeq平台和TruSeq DNA PCR-Free建库策略开展鸟枪法宏基因组测序；最后，采用Kneaddata、Trimmomatic和Bowtie2进行质控与污染序列去除，使用MetaPhlAn4进行物种水平分类注释和多样性分析，以HUMAnN3、ChocoPhlAn和UniRef90数据库进行功能基因家族及MetaCyc代谢通路预测分析。

在研究结果部分，论文首先报告了“3.1 Physicochemical parameters of Small Momella Lake’s water”。该部分通过现场理化测定证实，小莫梅拉湖具有典型极端环境属性。水温约为21.05 ± 0.62 °C至23.91 ± 0.08 °C，pH为9.75 ± 0.05至9.83 ± 0.03，盐度为11.97 ± 1.62 psu至14.01 ± 0.31 psu，显示该湖为显著碱性且中高盐的生态系统。研究据此指出，能够在该环境中存活的细菌必须具备独特的生物分子适应机制，这为后续嗜极端微生物的生物勘探提供了环境基础。

“3.2 Shotgun metagenomic analysis of bacterial species from Small Momela Lake”部分展示了该研究的核心分类学发现。经质控后，沉积物复合样本SMS1与水体复合样本SMW1分别获得约16,091,721和22,067,610条高质量序列。MetaPhlAn4分析显示，可识别序列全部归属于细菌域（Bacteria）。研究人员指出，尽管苏打湖通常也包含多样化古菌（Archaea），但本研究在严格标志基因覆盖阈值下未检出满足判定标准的古菌类群，原因可能与数据库参考基因组不足或相关类群分布于未分类部分有关。总体上共检测到27个物种，其中沉积物12个、水体25个，提示在本研究框架下水体复合群落显示出更高的可检出物种丰富度。分类上，这些物种分属5个门、12个纲、20个科和23个属。优势种包括Thioalkalivibrio denitrificans、Ca. GGB45747 SGB63468和Ca. GGB59200 SGB80967。论文特别强调，Thioalkalivibrio denitrificans在坦桑尼亚苏打湖中的检出属于首次报道，体现了其在东非苏打湖硫氮循环中的潜在生态地位。与此同时，多个Candidatus类群的大量存在说明未培养微生物可能是该极端生态系统中的关键成员。研究还指出，Glycocaulis profundi、Alkalicaulis satelles等若干物种可能构成苏打湖环境中的首次或区域首次报道，但作者对这些基于标志基因注释得出的物种级判定保持谨慎态度，认为仍需进一步验证。

“3.3 Bacterial diversity indices”部分对复合样本进行了α多样性（alpha diversity）和β多样性（beta diversity）描述。Shannon指数、Simpson指数及观测物种数均显示水体样本SMW1高于沉积物样本SMS1，提示在复合样本层面水体群落具有更高的丰富度与均匀度。研究人员指出，这可能与水体受光照、温度波动、风浪及营养不均匀分布影响更大，从而形成更多生态位有关；但同时强调，由于样本在测序前已混合，无法进行采样点间统计推断，因此这些差异只能作为探索性趋势。β多样性分析中，Bray–Curtis距离在水体中更高，提示其组成差异性更强；未加权UniFrac在沉积物中较高，暗示沉积物中可能包含系统发育上较独特但丰度较低的谱系；加权UniFrac在水体中较高，则表明水体优势类群在系统发育上相对更分散。

“3.4 Functional genetic potential of bacterial communities in small momela soda lake”部分揭示了群落功能潜力。研究共检测到350,147个UniRef90功能基因家族，其中84,674个为沉积物特异，167,607个为水体特异，另有相当部分为两类生境共享，反映出核心代谢功能与生境分化并存的格局。基于MetaCyc注释，最丰富的通路主要集中于核苷酸从头合成和氨基酸生物合成。高丰度通路包括PWY-7220和PWY-7222，分别对应dAMP和dGMP从头合成超通路，说明群落具备较强的DNA合成与修复前体生成潜力。PWY-7228、PWY-7221及其他与嘌呤/嘧啶代谢相关通路的广泛存在，进一步表明其在极端碱盐环境下维持核酸代谢的潜在能力。另一方面，VALSYN-PWY以及PWY-5103、PWY-7977、ARGSYN-PWY、HISTSYN-PWY、BRANCHED-CHAIN-AA-SYN-PWY、SERGLYSN-PWY等氨基酸合成通路在样本中也高度代表，尤其在水体中更为丰富，提示群落具有支持蛋白质合成及环境胁迫适应的广泛合成能力。

“3.5 Functional marker genes in small Momela Soda Lake”部分进一步从关键标志基因层面验证了功能推断结果。研究发现，两类样本均富含与pH稳态维持相关的阳离子/质子逆向转运蛋白基因家族，如mrp和nha，其功能是通过Na⁺/H⁺等离子交换维持胞质pH稳定，且在水体中的丰度高于沉积物，提示水体细菌为抵御碱性外环境需付出更高代谢代价。兼容溶质合成方面，L-ectoine synthase编码基因ectC在水体中约为沉积物的2倍，说明细菌可能通过合成甜菜碱样兼容溶质如ectoine来缓解渗透胁迫。硫循环方面，异化型硫酸盐/亚硫酸盐还原酶基因dsr在水体和沉积物中均高丰度存在，而硫氧化相关sox酶复合体基因在水体中更丰富，表明小莫梅拉湖可能存在由硫还原和硫氧化共同驱动的双向氧化还原硫循环。结合检出的Desulfohalovibrio alkalitolerans、Desulfonatronovibrio hydrogenovorans以及Thioalkalivibrio denitrificans、T. sulfidiphilus、T. thiocyanodenitrificans等类群，研究认为这些硫循环微生物很可能是上述关键功能基因的重要载体。碳固定方面，RuBisCO基因在两类样本中均有显著丰度，提示在富碳酸氢盐的苏打湖环境中，Calvin–Benson–Bassham循环可能是重要的初级生产途径，且可能主要由硫氧化细菌与蓝细菌共同承担。

论文讨论与结论部分强调，本研究首次从宏基因组层面对小莫梅拉湖进行了细菌分类多样性与推断功能潜力的系统描绘。结果表明，该湖虽然具有显著碱性和高盐特征，但仍栖息着多样化细菌群落，其中Pseudomonadota、Bacteroidota和Bacillota为主要门类。部分检出物种与既往其他苏打湖研究结果一致，另有若干物种可能构成全球、非洲或东非苏打湖中的首次观察，但这些判定仍需通过基因组解析、靶向分子验证或分离培养加以确认。功能上，群落表现出与核苷酸生物合成、氨基酸生物合成、pH稳态、硫代谢和碳固定相关的遗传特征，提示其具有适应极端碱盐条件的代谢灵活性。与此同时，作者也明确指出研究局限：样本DNA混合策略可能导致特定类群的低估或高估，且限制了空间异质性和统计学比较；采样仅限岸边和表层，深水及底栖微生物可能未被充分覆盖。基于此，研究建议未来应对各采样点独立测序，并开展旱季与雨季比较研究，同时通过宏转录组（metatranscriptomics）或定向培养等手段深化对Candidatus类群及关键功能的验证。

研究结论部分可译为：本研究首次对小莫梅拉湖进行了宏基因组学表征，提供了该湖细菌分类多样性及其与碱性、高盐环境相关推断功能潜力的探索性概览。尽管该湖具有极端理化特征，仍栖息着多样化细菌群落，其中以Pseudomonadota、Bacteroidota和Bacillota成员为优势。许多检出类群与其他苏打湖生态系统既往报道一致。此外，依据当前可获得文献和宏基因组参考数据库，部分物种水平注释可能构成苏打湖环境中的首次推定观察，但这些发现基于标志基因分型，应谨慎解读，仍需基因组解析宏基因组学、靶向分子验证或培养研究确认。小莫梅拉湖的宏基因组谱表明，该湖细菌群落具有与碱性和高盐条件下代谢灵活性相一致的遗传特征。与核苷酸和氨基酸生物合成相关的通路代表性较高，同时Na⁺/H⁺逆向转运蛋白、硫氧化（sox）、硫还原（dsr）及RuBisCO基因的共存，提示由Thioalkalivibrio、Desulfonatronum和Alkalilimnicola等类群主导的嗜盐嗜碱硫循环联合体，可能是该苏打湖生物地球化学过程的重要驱动者。这些观察结果为理解微生物群落如何在苏打湖严酷环境中持续存在提供了基础，但相关机制仍需时序性或功能性验证。