《Fermentation》：Diversity and Functional Genome Analysis of Lactic Acid Bacteria Isolated from Trás-os-Montes Artisanal Alheira
编辑推荐：
阿雷拉（alheira）是葡萄牙特拉斯-奥斯-蒙特斯（Trás-os-Montes）地区生产的一种传统烟熏、发酵、非即食肉香肠，其中乳酸菌（Lactic Acid Bacteria, LAB）是决定产品安全性与品质的主要微生物。研究人员采用全基因组测序对来自六
阿雷拉（alheira）是葡萄牙特拉斯-奥斯-蒙特斯（Trás-os-Montes）地区生产的一种传统烟熏、发酵、非即食肉香肠，其中乳酸菌（Lactic Acid Bacteria, LAB）是决定产品安全性与品质的主要微生物。研究人员采用全基因组测序对来自六个市镇手工阿雷拉中分离的59株乳酸菌进行了表征。共鉴定出隶属于6个属的8个物种，平均核苷酸一致性（Average Nucleotide Identity, ANI）分析归并为24个非冗余菌株组。功能注释显示出广泛的碳水化合物活性酶（Carbohydrate-Active Enzymes, CAZymes）、蛋白酶、脂肪酶、转运系统，以及酸、氧化、胆汁、热和盐胁迫耐受相关基因组合。研究人员以22个由无质量标记基因组组装体代表的菌株组进行基因型-表型推断。经错误发现率（False Discovery Rate, FDR）校正后，碳代谢、转运、能量代谢、总功能基因、酸胁迫、胆汁胁迫和总胁迫基因计数均与总pH下降呈正相关。混合效应模型显示标准化总功能基因计数与酸化呈正相关（边际R2，条件R2），但组内-组间分解表明该关联是在物种之间而非物种内可检测。正则化回归、随机森林（Random Forest）、AICc模型比较与偏最小二乘回归均指向碳代谢、转运、能量代谢及酸或胆汁胁迫功能为主要探索性信号；然而用于变量选择的稀疏LASSO回归不稳定，对未观察物种的预测能力有限。靶向的乙醛、柠檬酸、二乙酰、胞外多糖和乳糖代谢五条通路经多重检验校正后均无显著意义。这些发现为遴选本土发酵剂候选菌株提供了基因组与表型基础，同时强调需进行独立的功能与安全性验证。
研究背景方面，阿雷拉是葡萄牙东北部特拉斯-奥斯-蒙特斯具有文化与经济重要性的传统烟熏发酵肉香肠，由熟肉、面包、香料灌肠后经冷烟熏和干燥制成。手工生产不添加发酵剂，依赖原料与环境中土著微生物，乳酸菌为优势菌，影响安全、保质期和感官特性。该产品水分活度常高于0.98，零售生品易受单核增生李斯特菌（Listeria monocytogenes）、沙门菌（Salmonella enterica）和金黄色葡萄球菌（Staphylococcus aureus）污染。行业传统依赖冷冻、真空包装和亚硝酸盐等化学防腐剂，但加工肉被国际癌症研究机构列为1类致癌物，亚硝胺贡献结直肠癌风险，消费者与监管机构推动清洁标签（clean-label）防腐策略。乳酸菌因发酵能力与长期安全使用史成为天然替代，但功能基因组内容能否定量解释传统发酵食品土著乳酸菌的工艺表型，尤其在校正细菌基因组强分类结构后，仍不清楚。研究人员由此在既往表型工作基础上，用全基因组测序刻画阿雷拉乳酸菌功能与安全基因组，并将基因组预测与乳酸生成、蛋白水解、抗菌和pH动力学等表型整合，评估与发酵剂遴选相关的基因型-表型关系。
关键技术方法上，样本为来自特拉斯-奥斯-蒙特斯六市镇手工阿雷拉分离的59株乳酸菌，其中米兰德拉萨（Mirandela）20株、瓦尔帕苏什（Valpa?os）16株、维尼艾什（Vinhais）12株、维米奥索（Vimioso）6株、布拉干萨（Bragan?a）4株、莫加多鲁（Mogadouro）1株。研究人员以Illumina NextSeq 6000Dx平台测序，用Shovill/SPAdes流程de novo组装，以fastANI算ANI并按≥99.9%归并菌株组，用GTDB-Tk做分类与系统发育，kbDRAM做KEGG代谢注释，Prokka做技术与胁迫性状挖掘，AMRFinderPlus和Abricate做抗性与毒力筛查，并以线性混合效应模型、Random Forest、LASSO/Elastic Net、AICc与偏最小二乘回归做基因型-表型建模。
研究结果部分，3.1 生态与分类多样性：研究人员通过GTDB-Tk和ANI分析将59株定为乳杆菌目（Lactobacillales），8种6属；肠球菌（Enterococcus lactis）在米兰德拉萨高度克隆，明串珠菌（Leuconostoc mesenteroides）、清酒乳杆菌（Latilactobacillus sakei）、植物乳植杆菌（Lactiplantibacillus plantarum）等呈现明显地理分布，说明市镇间物种组成差异但采样不平衡，地理模式仅为描述性。
3.2 功能基因组表征：3.2.1 碳利用，研究人员发现EMP、PPP普遍保留，ED途径部分保留，lac操纵子在去冗余后仅肠球菌两组有lacY，且与物种完全混淆；3.2.2 碳水化合物活性酶，GH13近普遍，GH43、GH5有属种限制，GT2近普遍提示胞外多糖潜力，PL分布狭窄；3.2.3 胁迫耐受，F0F1-ATP酶普遍，精氨酸脱亚胺酶系统（arcA）与谷氨酸脱羧酶（gadB）有种级分布，氧化防御在ahpC与katA间互补，胆盐水解酶（bsh）和胆酰甘氨酸水解酶（cbh）种级分布，热休克与甜菜碱转运普遍；3.2.4 转运体，血红素、钴/镍、相容性溶质转运均呈种级结构；3.2.5 蛋白酶与3.2.6 脂肪酶，广泛存在但作为技术潜力而非酸化直接预测因子；3.2.7 安全谱，bceAB、tetA、vraDE等抗性转运系统、acm/sgrA黏附基因、hdcA/tdcA/ldcA/odc生物胺基因在部分菌株出现，提示非QPS（Qualified Presumption of Safety）菌株需表型MIC等验证。
3.3 基因组与表型整合：研究人员在菌株组水平发现E. lactis平均总pH下降最大但靶向通路基因最少，Lb. paracasei酸化较强且靶向通路基因最多；抗菌活性随病原、温度、物种变化，L. plantarum和Ln. mesenteroides对李斯特菌抑制较强。
3.4 基因型-表型关联与模型稳健性：3.4.1 菌株水平相关，运输与酸胁迫等计数与总pH下降FDR显著；3.4.2 主要混合模型与组内-组间分解，总功能基因计数正关酸化，但ICC约0.430，关联主要来自物种间而非物种内；3.4.3 正则化、集成与潜变量分析，LASSO不稳定，Ridge较优，Random Forest指碳代谢、能量代谢、酸胁迫重要，PLS一成分解释相当方差；3.4.4 AICc比较，碳+转运模型与胁迫模型几乎等权；3.4.5 模块与糖降解筛查，无单KEGG模块稳健，半乳糖醛酸降解具探索性相关；3.4.6 补充结果与敏感性，总功能基因不与乳酸生成或蛋白水解指数显著相关，删单株物种不改变结论，物种在多项基因组计数上差异显著但总pH下降分布无显著种间差。
讨论部分总结，研究人员指出面包淀粉基质下GH13与多糖系统反映基质适应，但途径完整性因种而异；蛋白水解/脂解基因为潜力而非表达证据；安全筛查标明E. lactis、Weissella viridescens、E. durans非QPS，S11、S12、S17、S19因缺乏抗性/毒力/生物胺基因成初步候选，但需标准安全、发酵与挑战试验。宽泛功能类别比单通路更能解释酸化，碳代谢/转运与胁迫是互补组分；变量选择不稳定源于小样本与高相关；混合模型信号主要为物种间协方差，不代表物种内可预测。未来需更大种内面板、长读长/杂交组装、转录组与在制品重复试验。
结论部分翻译，对59株阿雷拉源乳酸菌的全基因组测序解析出跨8个物种的24个非冗余菌株组，揭示碳水化合物代谢、转运、胁迫耐受、蛋白水解、脂解和安全相关基因的显著多样性。酸化更一致地关联宽泛功能，尤其是碳代谢、转运、能量代谢及酸与胆汁胁迫响应，而非任何单一靶向通路。主要混合模型识别出总功能基因计数与酸化的正相关，但组内-组间分析表明信号主要在物种之间，留一物种验证显示对未观察物种泛化能力有限。LASSO、Ridge、Random Forest、PLS与AICc比较广泛支持分布式多变量信号，同时也表明现有样本无法声称存在稳定稀疏预测集。候选菌株S11、S12、S17与S19是从相对有利基因组安全谱与技术特征平衡中识出的初步候选；而一些酸化更强的菌株如E. lactis具有需追加安全评价的基因组特征，这一权衡凸显需将技术性能与安全联合考量，而非仅按酸化能力选发酵剂。应用前这些菌株需抗菌药敏与生物胺表型确认、阿雷拉基质重复发酵试验、感官评估及病原挑战验证。该整合基因组、表型与统计框架为下一阶段实验提供了可辩护基础。论文发表于《Fermentation》。
