青贮贮藏期间高氮与低氮含量大麦籽粒微生物组的时间与功能特征分析

《Food Science & Nutrition》：Temporal and Functional Profiling of the Microbiome of High and Low Nitrogen Content Barley Seed in Silo Storage

　　

大麦籽粒品质受氮含量与贮藏条件影响，然而二者对籽粒微生物组组成与功能的影响尚不清楚。研究人员结合宏分类学（16S rRNA与ITS）分析、宏基因组测序与宏蛋白质组分析，表征大麦籽粒微生物组的结构与功能。自单一大麦品种（Kadie）收获的高氮（>1.5%）与低氮（<1.5%）籽粒，在收获时及青贮贮藏3、6、9个月后取样。扩增子测序揭示群落以变形菌门（Proteobacteria）、厚壁菌门（Firmicutes）和子囊菌门（Ascomycota）为主，宏基因组学证实Erwinia、Pantoea和Pseudomonas等属丰度较高，与既往大麦内生菌报道一致。虽鉴定出一致的核心微生物属集，其相对丰度在贮藏期持续变化。宏基因组分析显示高氮籽粒微生物组具快速代谢潜能但收获后下降，低氮籽粒群落维持长期代谢潜能。宏蛋白质组学确认这些功能转变揭示从活跃生长到胁迫耐受的时间转换。本研究增进对长期贮藏期间大麦籽粒微生物组的认识，为优化麦芽与酿造贮藏提供见解。

研究背景方面，大麦（Hordeum vulgare L.）作为酿造啤酒的核心谷物，其籽粒品质直接决定麦芽与啤酒风味，而收获后青贮贮藏（silo storage）是打破种子休眠的标准工业环节。现有研究多聚焦单时间点或培养法，对贮藏期微生物组的时间演替与功能转换认识不足；同时氮含量影响啤酒品质已有定论，但其与贮藏微生物动态的关系尚不明确。为此，研究人员以南非西开普省Anheuser-Busch InBev设施中单一品种Kadie大麦为材料，按籽粒总氮（Dumas燃烧法测定）分为高氮（>1.5%）与低氮（<1.5%）两组，在收获及贮藏3、6、9月四个时点各取三重复，共24份样品，整合多组学解析微生物组演替。该研究发表于《Food Science & Nutrition》。

关键技术方法上，研究人员构建收获与贮藏3/6/9月×高/低氮的双因素样本队列（每时点每氮组3重复，宏基因组与宏蛋白质组各取8个子集样）。分类学采用16S rRNA V1–V3区与ITS1扩增子DADA2流程（SILVA v138.2、UNITE v10.0注释）；宏基因组用Illumina NovaSeq 6000双端测序，Kraken2/Bracken分类、MetaSPAdes组装、Prokka预测、eggNOG-mapper与KEGG/COG/CAZy功能注释；宏蛋白质组经PBS超声裂解、SDS-PAGE分带、胰酶消化、LC–MS/MS（AB Sciex 5600 TripleTOF）采集，Protein Pilot比对自建数据库（含大麦宿主蛋白与宏基因组注释蛋白）并以靶-诱饵法控假发现率（false discovery rate, FDR≤1%），Unipept做肽段分类与功能注释。

研究结果部分，3.1细菌与真菌α多样性显示，贮藏时长单独不显著改变细菌丰富度，但在各氮组内显著；高氮组收获时多样性最高，6月降至最低（Observed ASVs 14–28，Shannon 1.99–2.47），9月回升与低氮组趋同；低氮组3月内低位，6月反弹至收获水平。真菌α多样性随贮藏显著波动，3月峰值，6月回落，9月再升，氮含量主效应不显著。3.2细菌与真菌β多样性中，细菌结构受“贮藏时长×氮含量”交互作用主导（R²=26.48%，p=0.003），真菌结构仅由贮藏时长驱动（R²=55.66%，p=0.001）；PCoA显示9月样品聚集成簇。3.3分类组成与演替表明，细菌以变形菌门（均值75.5%）、厚壁菌门、蓝菌门为主，真菌91.1%为子囊菌门且Alternaria占90.7%；宏基因组子集印证上述 phylum 级轮廓。3.4属级多组学演替指出，核心七属（Pseudomonas、Pantoea、Chryseobacterium、Pedobacter、Stenotrophomonas、Paenibacillus、Massilia）跨三组学一致检出；收获时高氮组呈Pseudomonas(31.6%)+Pantoea(16.3%)+Chryseobacterium(10.9%)+Stenotrophomonas(10.2%)组合，低氮组几仅Pantoea(51.0%)+Pseudomonas(35.0%)；9月两组合流于Pantoea~37%，宏蛋白质组见Pseudomonas肽段占~70%。3.5宏基因组组装与COG注释获274,778 contigs、952,951 CDS，功能以转录(K)、氨基酸转运代谢(E)、无机离子转运(P)、碳水转运(G)为主，但COG与CAZyme剖面在贮藏/氮间无统计差异（PERMANOVA p>0.1）。3.6 CAZyme谱中GT2/GH13/GT4居首，6月低氮样GT1/GH18/GH23簇升高，9月大部家族z值转负。3.7宏蛋白质组功能分析检出2694 InterPro蛋白，外膜蛋白（outer membrane proteins, OMPs）、过氧化氢酶（catalase）、孔蛋白（porin）、TonB依赖受体为主；高氮组6月catalase达26.1%，低氮组3月TonB受体24.8%；6–9月高氮组注释蛋白覆盖率由85%降至25%，仅宏蛋白质组层PERMANOVA显示贮藏时长显著（R²=53.1%，p=0.033）。

讨论部分总结，贮藏时长是细菌与真菌群落结构的主驱因子，氮含量仅通过与贮藏交互影响细菌（真菌无此交互）；核心七属跨组学一致提示持久定殖而非污染；收获期高氮组群落异质、低氮组极简，后期向Pantoea主导及Paenibacillus等胁迫耐受型过渡，符合富营养向寡营养转换；非核心属（Lactococcus、Sphingobacterium于9月低氮样）反映局域微环境但受子集单重复限制不能泛化；真菌Alternaria单一优势（>90%）构成霉菌毒素风险，宏基因组/蛋白质组低频系细菌/真菌生物量比与技术局限所致；COG/CAZyme潜能稳定而宏蛋白质组表达随贮藏转向营养 scavenging 与氧化胁迫应答（catalase升、TonB受体早期高），9月蛋白量衰减暗示活而非培养态或休眠。结论部分译文如下：本研究结果显示贮藏时长是各数据集微生物群落结构的主导驱动因子，籽粒氮含量仅在细菌群落的“贮藏时长×氮含量”交互中具可检测效应。细菌群落响应交互作用，真菌群落仅随贮藏时长演替。组学间功能分异明显：编码功能潜能相对稳定，而表达的宏蛋白质组追踪贮藏时长并呈现从收获期表面定殖、到营养 scavenging、直至9月代谢活性降低的生理进程。产毒Alternaria在所有样品中持续存在构成食品安全隐患，提示需优化贮藏策略；Fusarium、Aspergillus、Penicillium低频或反映该贮藏系统竞争结构。DNA法不能区分生理状态，故以宏蛋白质组辅助解读。宏基因组/蛋白质组子集规模有限虽限制家族级统计力，但趋势表明多组学设计对捕获贮藏谷物微生物组功能响应不可或缺，单组学会低估之。