环境Aspergillus分离株的表型与基因组整合表征及其柠檬酸生产潜力

《MicrobiologyOpen》：Integrated Phenotypic and Genome-Informed Characterization of Environmental Aspergillus Isolates With Potential for Citric Acid Production

　　

Aspergillus（曲霉属）物种在工业生物技术中被广泛应用，尤其用于柠檬酸生产。本研究采用表型与基因组信息整合的方法，对来自科特迪瓦（C?te d'Ivoire）的具有柠檬酸生产潜力的环境Aspergillus分离株进行了表征。从香蕉根际土壤及腐烂果蔬中回收的分离株，先通过氨蒸气试验（ammonia vapor assay）筛查黄曲霉毒素相关反应，并在添加溴甲酚绿（bromocresol green）的Czapek-Dox琼脂上筛查酸化能力。入选分离株经全基因组测序（whole-genome sequencing）、Kaiju分类学分析和EggNOG-mapper功能注释进行分析。柠檬酸产量在固态发酵（solid-state fermentation）条件下采用校准比色法（calibrated colorimetric assay）定量。在回收的110株真菌分离株中，24株Aspergillus分离株被归为7个形态型（morphotypes）。4株分离株（E1-E4）在试验条件下表现出酸化活性且未检出黄曲霉毒素相关反应。形态学与全基因组分类学分析保守地将E1-E3归入Aspergillus niger species complex（黑曲霉复合群），而E4在分类学上仍未确定。E2分配到A. niger参考序列的reads比例最高（92%），且与柠檬酸代谢和碳水化合物降解相关的预测功能代表性最强。在固态发酵过程中，E2在测试条件下表现出最高的实测柠檬酸浓度。比较基因组分析亦显示E1-E3之间的相似性高于其与E4的相似性。总体而言，E2成为进一步开展柠檬酸生产研究的最有前景的候选菌株，支持将定量表型评估与基因组信息表征相结合用于环境真菌分离株研究的价值。

研究背景与问题：丝状真菌（filamentous fungi）是工业生物技术中最重要的微生物之一，可产生有机酸、酶、抗生素及生物活性化合物等。其中Aspergillus niger（黑曲霉）被认为是最高效的柠檬酸微生物生产者之一，全球柠檬酸年产量超过200万吨，A. niger因高产率、遗传稳定性和碳源利用广泛性而占据工业生产的绝大部分。持续寻找新的Aspergillus菌株仍是工业微生物学的重要目标，农业土壤、根际环境和腐烂植物材料是Aspergillus的丰富生态位。然而，环境Aspergillus分离株的开发利用面临挑战：其一，Aspergillus niger species complex内部存在显著的表型与遗传多样性，准确分类鉴定困难；其二，部分Aspergillus物种可产生霉菌毒素（mycotoxins），可能影响其工业应用安全性。传统筛选依赖形态学观察和生化筛查（产酸试验、霉菌毒素检测），但难以反映代谢能力的遗传基础，也可能无法准确反映近缘种间的分类关系。基因组测序与生物信息学工具的进步，使分类学分析和功能注释成为可能。因此，将传统微生物学筛查与基因组信息分析相结合，可为具有生物技术价值的环境真菌分离株的初步筛选提供更全面的框架。

研究开展与结论：研究人员从科特迪瓦采集香蕉根际土壤及腐烂果蔬样本，共100份独立环境样本，回收110株丝状真菌分离株，其中24株鉴定为Aspergillus属并归为7个形态型。基于酸化能力和氨蒸气试验未检出黄曲霉毒素相关反应，筛选出E1-E4四株进行全基因组测序、分类学分析和功能注释，并以菠萝废弃物为底物进行固态发酵定量柠檬酸产量。结果显示，E2在酸化表型、基因组分类学归属、柠檬酸代谢与碳水化合物降解相关预测功能以及实测柠檬酸产量方面均表现最优，被确定为最有前景的候选菌株。研究意义在于证明了定量表型筛选与基因组信息分析相结合在环境真菌分离株表征与优选中的价值，同时明确指出色谱确认、发酵工艺优化、全面霉菌毒素评估和多基因座分类学验证是工业应用前的必要步骤。

关键技术方法：研究人员采用氨蒸气试验在椰子提取物琼脂（Coconut Extract Agar，CEA）上筛查黄曲霉毒素相关反应；采用添加溴甲酚绿的改良Czapek-Dox琼脂进行酸化能力定性筛查并定量黄色晕圈直径；采用固态发酵以工业菠萝废弃物为底物，结合Marier与Boulet比色法（经Dhillon等改良）定量柠檬酸；采用改良CTAB法提取基因组DNA，经Illumina NextSeq平台双端测序（2 × 150 bp）；采用FastQC与Trim Galore进行质控，SPAdes组装、QUAST评估，Kaiju对照RefSeq Fungi数据库进行分类学分析，EggNOG-mapper进行功能注释并交叉参照KEGG柠檬酸循环（TCA循环，map00020）通路；采用方差分析（ANOVA）与Tukey's HSD检验进行统计比较。样本来源为科特迪瓦Yopougon、Koumassi、Grand-Bassam的柠檬、洋葱、胡萝卜及Bonoua的香蕉根际土壤。

研究结果：

3.1 真菌分离株的多样性与形态学表征：从环境样本回收110株丝状真菌，其中腐烂果蔬70株、香蕉根际土壤40株，24株初步鉴定为Aspergillus并分为7个形态型。各形态型在菌落颜色、质地、产孢模式和分生孢子形态上存在差异。E1-E3具双列（biseriate）瓶梗、辐射状头状体和粗糙壁分生孢子，符合Aspergillus section Nigri（黑曲霉组）特征；E4为单列（uniseriate）瓶梗、囊泡较小（20.0–35.0 μm）、疏松柱状头，形态学上明显区别于E1-E3。

3.2 黄曲霉毒素相关反应与酸化能力筛查：A2（A. candidus）与O5（A. flavus）呈极强阳性反应（++++），E1-E4及G3在所有重复中均未检出反应（ND）。酸化晕圈直径在四株间差异显著（p < 0.05），E2最大（57.50 ± 1.15 mm），E3与E1次之（44.50 ± 0.70 mm、41.40 ± 0.82 mm），E4最小（28.50 ± 0.60 mm）。酸化活性最高的分离株恰为未检出黄曲霉毒素相关反应者。

3.3 入选分离株的分类学表征：E1-E3在门、纲、目、属水平分别以Ascomycota、Eurotiomycetes、Eurotiales、Aspergillus为主，E4分类学构成更为异质。E1-E3组装基因组大小为29.95–35.76 Mb，E2组装连续性最高（scaffold N50为8211 bp）、BUSCO完整度最高（88.2%）；E4组装仅4.56 Mb且未能确定BUSCO完整度。E1、E2、E3中分配至A. niger参考序列的reads比例分别为64%、92%、65%，另有A. welwitschiae与A. luchuensis的分配。综合形态学、分类学谱、物种参考分配和组装特征，E1-E3被保守归入A. niger species complex，E4保留为分类学未确定。

3.4 分离株间比较基因组关系：基于共享基因组特征的分析显示E1-E3相似性高，E4相似性较低。基于注释基因含量的层次聚类中，E1与E3距离最短，E2形成近缘分支，E4单独分支；该聚类反映描述性基因组相似性，不构成系统发育证据。

3.5 柠檬酸代谢相关功能潜力：功能注释鉴定出柠檬酸合酶（citrate synthase，EC 2.3.3.1）、乌头酸酶（aconitase，EC 4.2.1.3）、异柠檬酸脱氢酶（isocitrate dehydrogenase，EC 1.1.1.41/42）和异柠檬酸裂解酶（isocitrate lyase，EC 4.1.3.1）等预测功能，E2的柠檬酸代谢相关预测功能代表性最高。但柠檬酸积累受碳通量、酶调控、底物可利用性、微量元素、pH等条件控制，注释仅代表预测代谢潜力。

3.6 碳水化合物降解功能潜力：注释显示α-淀粉酶（α-amylases）、β-葡萄糖苷酶（β-glucosidases）、葡聚糖酶（glucanases）和糖基水解酶（glycosyl hydrolases）等预测功能，E2代表性最高，E1与E3居中，提示分离株动员复杂碳源的预测能力存在差异。

3.7 固态发酵过程中的柠檬酸定量产量：分离株间差异显著（单因素ANOVA，p < 0.001）。E2全程保持最高实测浓度，第12天达最大值0.302 g/kg；E1、E3、E4最大值分别为0.240、0.180、0.140 g/kg，E2显著高于其他菌株（Tukey's HSD，p < 0.001）。产量在第12–15天趋于平台。实测产量顺序（E2 > E1 > E3 > E4）与基因组预测功能顺序一致，但该一致性不构成因果关系。

3.8 研究局限性：未进行HPLC等色谱确认；氨蒸气试验仅为推定筛查，赭曲霉毒素A（ochratoxin A）和伏马菌素B₂（fumonisin B₂）未定量；E1-E3未达物种级确定鉴定，E4分类学未确定；功能分析基于基因组注释，未验证基因表达、酶活性或代谢通量。

讨论与结论：讨论部分指出，腐烂果蔬中真菌分离株占优势可能与简单碳水化合物丰富和水分条件适宜有关；七个形态型的表型差异反映了环境Aspergillus种群的表型多样性及对生态位的适应；E1-E3与A. niger species complex的归属具有生物技术相关性，但分类归属本身不能预测柠檬酸生产能力；比较基因组相似性与既往报道的黑曲霉近缘种基因组多样性一致；E2的基因组预测模式与发酵表型的一致性需基因表达、酶活性、蛋白质组或代谢通量分析方能确立因果关系；本研究获得的柠檬酸浓度（0.140–0.302 g/kg）远低于优化工业发酵体系，结果应作为菌株筛选与后续工艺优化的基础而非工业性能指标。结论译文如下：本研究结合传统微生物学筛查、定量表型评估、固态发酵和基因组信息分析，表征了具有柠檬酸生产潜力的环境Aspergillus分离株。在从环境样本回收的24株Aspergillus分离株中，基于酸化能力及氨蒸气试验条件下未检出黄曲霉毒素相关反应，筛选出4株分离株（E1-E4）。形态学表征与全基因组分类学分析支持E1、E2和E3归属Aspergillus niger species complex，而E4在分类学上仍未确定。ITS区域和β-微管蛋白（benA）位点已从E1-E3基因组组装中经计算机（in silico）成功获取，提供了补充性分类学序列信息，但确定性的物种级鉴定仍需专门的多基因座系统发育分析。定量酸化测定与固态发酵显示分离株间存在显著差异，E2在测试条件下表现出最强的酸化表型和最高的实测柠檬酸浓度。功能注释亦显示E2具有最高代表性的柠檬酸代谢和碳水化合物降解相关预测功能，但这些基因组衍生注释仅代表预测代谢潜力，不能证明基因表达、酶活性或代谢通量。总体而言，定量表型数据、发酵结果与基因组信息分析的汇聚表明，E2是所评估分离株中进一步开展柠檬酸生产研究最有前景的候选菌株。然而，在考虑任何潜在工业应用之前，仍需进行柠檬酸的色谱确认、发酵工艺优化、全面霉菌毒素评估和多基因座分类学验证。本研究证明了将定量表型筛查与基因组信息分析相结合，用于具有生物技术价值的环境真菌分离株表征与优选的价值。