《Biology》：Microbial Functional Potentials Differ Among Monospecific and Mixed Moss Biocrusts in an Alpine Sandy Ecosystem
编辑推荐：
研究人员采用鸟枪法宏基因组测序（shotgun metagenomic sequencing）比较了青藏高原东北缘共和盆地中，与扭叶墙藓（Didymodon constrictus）单种结皮（mossC）、锈色墙藓（D. ferrugineus）单种结皮（mo
研究人员采用鸟枪法宏基因组测序（shotgun metagenomic sequencing）比较了青藏高原东北缘共和盆地中，与扭叶墙藓（Didymodon constrictus）单种结皮（mossC）、锈色墙藓（D. ferrugineus）单种结皮（mossF）及视觉共优势混生结皮（mossM）相关的结皮下土壤微生物功能谱。每类结皮设15个空间分离样方，合并为3个复合生物学重复（有效n=3）。KEGG、CAZy及靶向碳氮循环注释显示，各类结皮在碳水化合物加工、碳氮转化相关基因的相对丰度上存在类别关联差异。mossC中糖苷水解酶（glycoside hydrolases, GH）及若干复杂碳加工功能平均代表度更高；mossF中氮同化与乙酸相关功能代表度更高；mossM中部分碳降解、氮矿化及异化硝酸盐还原功能代表度更高。完整冗余分析（RDA）模型解释58.4%功能变异（校正R2=0.334；精确置换p=0.028340），总碳（TC）与土壤含水量（SWC）保留为单变量关联。属级分类差异在FDR校正后与KEGG、CAZy、碳循环及氮循环差异矩阵相关。mossM组均值相对两单种近似非加权中点呈正负描述性偏离，但未对这些偏离做推断检验。复合重复数少、藓类比例为视觉估算、土壤属性混杂且未测过程速率，限制了因果与确证性解释。
研究背景与立项依据
生物土壤结皮（biological soil crusts, biocrusts）是由蓝细菌、藻类、真菌、地衣、苔藓及异养微生物组成的地表薄层生命系统，在侵蚀控制、水分再分配及碳氮循环中具有重要作用。苔藓主导结皮较早期蓝细菌结皮具更高表面粗糙度、生物量与持水能力，但其下土壤微生物功能基因库在寒冷风蚀高山沙质生态系统中仍缺乏解析。扭叶墙藓（Didymodon constrictus）与锈色墙藓（D. ferrugineus）同属丛本藓科（Pottiaceae），分布于青藏高原高山草地。苔藓无真根，但通过假根形成密切接触界面，干湿循环释放可溶物、凋落组织及截获颗粒可为微生物提供底物并改变微生境水热条件。共和盆地属低温、少雨、高蒸发、强辐射与频繁风蚀的高山沙质区，约三分之一土地沙化，氮素尤为受限。在此背景下，比较单种与混生藓类结皮下微生物功能潜力，可探明藓类身份及共存是否关联不同的氮获取、保留与潜在损失能力。该研究发表于《Biology》。
主要技术方法概括
研究人员在青海共和盆地沙珠玉试验站布设三类藓类结皮（mossC、mossF、mossM），每类取15个10 cm×10 cm样方（相邻间隔≥10 m），5个合并为1个复合生物学重复，每类3重复（有效n=3）。去除可见结皮层后采集其下0–2 cm矿质土，?80 ℃保存。提取DNA后行Illumina NovaSeq 6000双端150 bp测序，fastp质控，MEGAHIT独立组装，Prodigal预测ORF，CD-HIT聚类非冗余基因，Bowtie2比对以TPM归一化。功能注释依托Kraken2（NCBI RefSeq）、KofamScan（KEGG）、DIAMOND（NCycDB、CAZyDB）、HMMER（dbCAN HMMdb）。统计以复合重复为单位，采用精确置换检验（one-way F）、BH校正、RDA（vegan）、Spearman相关及Mantel检验，混生结皮中点偏离按两单种均值近似非加权中点计算描述性偏离百分比。
研究结果
3.1 土壤理化性质
精确置换检验与BH校正后，砂、粉、黏、有机质（SOM）、总碳（TC）、总氮（TN）、土壤含水量（SWC）在三类间存在整体差异（pH无）。组均值上mossC砂最低而粉、黏、SOM、TC、TN、SWC最高，mossF砂最高，mossM居中；但24个两两比较BH校正后均无显著，仅作描述。
3.2 回收非冗余基因数
各样品组装所得非冗余基因数中值mossF<mossM<mossC；mossM组内跨度最大由1个复合样驱动，无推断比较，仅描述。
3.3 KEGG与CAZy功能谱
22个KEGG二级类中14类BH校正后组间差异显著。mossM在脂质代谢、辅因子维生素代谢、次生代谢物生物合成、萜类聚酮代谢、聚糖生物合成与代谢、膜转运、原核细胞群落、复制修复均值最高；mossF在碳水、能量、核苷酸代谢及转录最高；mossC在翻译及运输分解最高。CAZy中GH在mossC最高，CBM在mossF最高，CE与PL在mossM最高；GH、CBM、CE、PL组间差异显著（GT、AA不显著），数值为TPM相对组成非绝对拷贝。
3.4 碳循环功能与选定基因
三碳通路（碳降解、碳固定、还原型三羧酸循环rTCA）组间均显著。碳降解mossM最高，碳固定mossC最高，rTCA mossF最高。五个选定基因中乙酰辅酶A合成酶（ACSS1_2/acs）与异柠檬酸裂解酶（ICL/aceA）mossF最高；rfbB/rmlB/rffG与甘露糖苷酶（MAN）mossM最高；木糖调控相关（xylS/yicI）mossC最高。均为功能潜力非活性证明。
3.5 氮循环功能与选定基因
氮同化mossF最高，氮矿化与异化硝酸盐还原为铵（DNRA）mossM最高，反硝化整体占比低且mossF最高；四通路组间均显著。九选定基因中谷氨酸合酶大亚基（gltB）与谷氨酸脱氢酶（GDH2）mossC最高，gdhA与一氧化氮还原酶（norB）mossF最高且这四项显著；nasB、narB、nrfH、nosZ、napA不显著。
3.6 环境关联与分类-功能一致性
全RDA（含pH、TC、SWC）解释58.4%变异（校正R2=0.334，精确置换p=0.028340）；单变量中TC与SWC模型显著，pH不显著。128个Spearman相关BH校正后留10个：CBM与砂正相关；xylS/yicI与砂负相关、与粉、SOM、TC、TN、SWC正相关；GDH2与砂负相关、与粉、SWC正相关。属级分类Bray-Curtis差异与KEGG、CAZy、碳、氮功能差异矩阵Mantel检验FDR后均显著正相关（r=0.444–0.682），仅示矩阵协同变异不指认携基因 taxa。
3.7 混生结皮中点描述性偏离
mossM相对两单种非加权中点，CBM、CE、碳降解、氮矿化、DNRA、nosZ、nrfH为正偏离；GH、AA、rTCA、氮同化、反硝化、acs、ICL/aceA、gltB、GDH2、narB、norB为负偏离，norB负偏离超60%。无推断检验，仅为探索性描述。
讨论部分总结
研究人员指出结皮类别、伴生土壤条件与微生物功能基因谱在同点共变，但无法区分土壤驱藓、藓改土壤或互馈。RDA与相关将xylS/yicI、GDH2与SWC正向关联解释为湿态下碳水底物可及性与铵代谢的潜力模式，非直接调控证据。mossM非冗余基因跨度受单样品与技术因素混合影响。碳加工上mossC富GH与xylS/yicI指向水解潜力，mossF富acs与aceA指向乙酸激活与乙醛酸支路，mossM碳降解及CBM、CE正偏离或反映底物谱更广，但TPM相对归一致比例重分配。氮方面mossF氮同化最高不悖于低TN，因基因为获取潜力非池大小；mossM的DNRA高或借藓源可溶物C:N偏利还原，但缺组织C:N与微氧实测。混生中点正负偏离更似功能比例重分布而非均匀增强，norB大幅负偏与nosZ正偏说明基因非协同如通路。视觉估比、属性混杂、n=3、无表达与通量测度为四大限制。
结论部分翻译
在所调查高山沙质点，扭叶墙藓、锈色墙藓及视觉混生藓类结皮下微生物功能谱存在差异并与伴生土壤条件共变。mossC、mossF、mossM在碳水加工、碳代谢及氮转化相关基因相对代表度上对比鲜明。mossM组均值相对两单种近似非加权中点呈正负描述性偏离，但非统计确认的非加和效应。独立复合重复少、藓比视觉估、土壤混杂及缺过程速率测度限制因果推论与外推。验证需增生物学重复、定量混比、13C/15N示踪、宏转录组、酶活及温室气体通量测定。
