染色体尺度基因组分析揭示腐乳来源红曲霉Monascus ruber菌株BC20中橘霉素相关区域位点特异性破坏

《Foods》：Chromosome-Scale Genome Analysis Reveals Locus-Specific Disruption of the Citrinin-Associated Region in a Furu-Derived Monascus ruber Strain BC20

　　

红曲霉（Monascus）属物种广泛用于传统发酵食品中的色素与风味形成，但橘霉素（citrinin）污染仍是限制其更广泛应用的主要安全问题。因此，本研究旨在通过整合跨食品相关基质的表型筛选与基因组解析分析，评估一株腐乳来源的M. ruber菌株BC20的橘霉素风险。在包括真菌培养基及乳基、谷物基、麸皮基等8种基质中培养14天后，采用免疫亲和柱净化结合高效液相色谱-荧光检测法（HPLC–FLD，检出限4 μg/kg，定量限12 μg/kg）未检出橘霉素。为探究该表型的遗传基础，研究人员构建了BC20的染色体尺度基因组组装，并在总计19个红曲霉基因组间进行比对分析。平均核苷酸一致性（ANI）分析与系统基因组推断一致将BC20归入ruber–pilosus支系。比较共线性分析显示，BC20中橘霉素相关基因座不再保留完整簇构型，而呈现残基基因座（remnant-locus）架构，同支系数个相关基因组亦见类似模式。相比之下，莫纳可林K（monacolin K, mk）基因座在BC20中保持共线性保守，支持位点特异性结构扰动而非组装伪缺失。此外测定了其抗真菌敏感性。综上，BC20代表一株橘霉素未检出的M. ruber候选菌株，本研究为食品相关红曲霉分离株的橘霉素风险筛查提供了实用分析框架。

研究背景方面，红曲霉（Monascus）属作为东亚传统发酵食品（如红曲米、腐乳）中重要的丝状真菌，承担色素形成、风味发育及次生代谢物产出功能，并逐步拓展至谷物与乳基体系。然而橘霉素（citrinin）这一聚酮类真菌毒素在部分红曲发酵产品中被频繁检出，中国对红曲红设定0.04 mg/kg上限，欧盟对红曲米膳食补充剂设100 μg/kg上限，安全合规构成产业化障碍。既往研究指出橘霉素生物合成由以ctnA等为核心的基因簇（biosynthetic gene cluster, BGC）调控，工艺优化或转基因沉默虽可降低产量，但食品端受法规与消费者接受度制约。因此从环境源筛选天然稳定橘霉素阴性菌株并以可重复遗传证据支撑，成为更务实路径。本研究核心科学问题在于：天然分离株在跨食品相关条件下橘霉素未检出表型是否稳定，且该表型究竟源于橘霉素相关基因座的结构性破坏，还是完整基因簇的条件性沉默。研究人员以武汉农家腐乳分离所得M. ruber BC20为对象，提出假设——跨食品条件一致未检出橘霉素表型应伴随该基因座结构性破坏而非条件沉默，并设计表型—基因组联合验证。

关键技术方法上，研究人员从武汉腐乳样品经MEA与PDA纯化获得BC20等10株红曲霉分离株；采用ITS与β-微管蛋白（BenA）序列作分子鉴定并建邻位连接树。表型端设8套培养体系（YES、CY20、RC、RM、CM、WBM、RCW及市售灭菌乳FM），30℃培养14天，每条件3生物学重复，免疫亲和柱净化-HPLC–FLD定性筛查橘霉素（LOD 4 μg/kg，LOQ 12 μg/kg）。基因组端采用PacBio长读长加Illumina短读长混合组装成染色体尺度，AUGUSTUS与GeneWise联合预测6720个蛋白编码基因，BUSCO评估完整性99.3%。比较基因组纳入18个公共红曲霉基因组共计19株，用skani算ANI、OrthoFinder提单拷贝同源基因建系统基因组树、minimap2做组装间共线性、以M. ruber M7（KT781075.1/BGC0001338）橘霉素BGC为参考做序列级比对与变异注释，并以mk基因座作内对照排除组装伪缺失。抗真菌药敏按CLSI M38微量液稀释法测7种药剂MIC/MEC。

研究结果部分，3.1节通过在8种培养体系中HPLC–FLD色谱均无8.97 min特征峰，证明BC20在真菌培养基、乳基液体、谷物麸皮固体基质中14天培养后橘霉素均未检出（ND），表型跨基质稳定。3.2节形态学呈橙红色绒状菌落、分隔菌丝与亚球状分生孢子，ITS与BenA系统树均归M. ruber；混合组装得26.19 Mb基因组、GC 48.85%、6720编码基因、染色体尺度连续性高，降低把组装断裂误判为生物缺失的风险。3.3节BUSCO完整度99.3%，ANI热图与单拷贝同源蛋白系统基因组树一致将BC20定位ruber–pilosus支系，与M. purpureus明显分离。3.4节全基因组共线性显示BC20与该支系参考基因组线性保守，核心同源群占比大、变异限于附属组份，说明基因组骨架保守、便于靶点解读。3.5节橘霉素相关基因座共线性比对揭示：M. purpureus参考株保留含pksCT及修饰调控基因的完整簇，BC20对应共线性区呈残基架构（部分丢失/打断），同支系数株有相似样式；断点附近共享DDE转座酶插入位点，LINE样元件位置株间可变，BC20该区检出3个候选反向重复（IR）元件，保留编码区无start_lost/stop_gained变异，转座元件共定位仅作结构关联而非因果机制直证。3.6节mk基因座含mkA/mkC–mkH等核心基因且共线性保守、无大片段缺失，排除橘霉素区“消失”来自组装崩坏或次级代谢区普遍退化，证位点特异性扰动。3.7节COG显示碳水转运与代谢类富集，CAZyme以糖苷水解酶（GH）与糖基转移酶（GT）为主，分泌组跨属保守，提示底物利用潜力。3.8节CLSI M38药敏：卡泊芬净MEC 0.0156 μg/mL、两性霉素B MIC 0.5 μg/mL、伊曲康唑2 μg/mL、酮康唑8 μg/mL、氟康唑32 μg/mL、伏立康唑>32 μg/mL、5-氟胞嘧啶16 μg/mL。

讨论部分指出，跨8基质（营养、水分活度、氧传质差异大）一致ND支持表型在食品相关条件下稳定，且单菌接种排除竞争群落干扰；染色体尺度组装使BGC区真缺失与碎片假缺失可分辨；在ruber–pilosus支系保守骨架背景下，BC20橘霉素区残基架构与同支系共享DDE断点样式，暗示支系水平历史重排，LINE可变插入为株后事件，转座子共定位仅作结构关联待断点跨越读长/PCR验证；mk保守对照强化位点特异性论点；CAZyme/分泌组提示发酵底物利用潜质但需实基质验证产色产香与水解活性。结论部分翻译如下：本研究在红曲霉M. ruber BC20中鉴定出橘霉素相关区域位点特异性结构破坏。在8种培养基与食品相关基质中培养14天后均未检出橘霉素，与染色体尺度比较基因组所揭示的橘霉素相关基因座残基架构一致；相反mk基因座保持结构保守，进一步支持该基因组改变之位点特异性。表型与基因组证据相符，支持BC20作为橘霉素未检出M. ruber候选菌株，供谷物发酵食品与乳发酵体系进一步评估。鉴于目前缺基质特异性分析验证，后续开发应聚焦加标回收与基质效应验证、扩展安全评价及中试发酵试验。该文发表于《Foods》。