本文为一项关于黑曲霉组（Aspergillus section *Nigri*）中柠檬酸高产菌株Aspergillus lacticoffeatus WU-2020的基因组测序研究报告。该研究由Isato Yoshioka, Hiroki Takahashi, Yoko Kusuya, Takashi Yaguchi, Akira Shibata, Kohtaro Kirimura等作者共同完成，研究机构包括早稻田大学（Waseda University）的研究与工程研究所、科学与工程学院，以及千叶大学（Chiba University）的医学真菌学研究中心、分子手性研究中心和植物分子科学中心。该研究成果以题为“Draft genome sequence of Aspergillus lacticoffeatus WU-2020, a citric acid producer suitable for solid culture that belongs to Aspergillus section *Nigri*”的论文形式，于2023年1月4日发表在Microbiology Resource Announcements期刊（2023年，第12卷，第1期）。
研究的学术背景 柠檬酸（Citric Acid, CA）是一种用途广泛的有机酸，主要通过丝状真菌的工业发酵生产，其中黑曲霉组的菌株是主要的生产者。研究团队此前已报道了适用于深层发酵（submerged culture）的柠檬酸高产菌株Aspergillus tubingensis WU-2223L的基因组序列。然而，固态发酵（solid culture）在特定工业应用中具有优势，例如利用农业废弃物（如甘蔗渣）作为载体进行生产。为了发掘适用于固态培养的柠檬酸高产菌株，研究团队从日本土壤中分离筛选真菌。其研究目的在于获得一株适用于固态培养的柠檬酸高产菌株，并解析其基因组序列，以期为理解柠檬酸高产机制、评估菌株的工业应用潜力（包括生物安全性）奠定遗传学基础。
详细的研究流程 本研究是一个典型的微生物资源基因组测序与表征项目，主要包括菌株分离鉴定、基因组测序与组装、基因组注释与分析、以及初步的产酸能力与安全性评估四个主要环节。
第一环节：菌株的分离、筛选与鉴定 研究团队从日本采集的土壤样本开始。具体流程如下： 1. 样品处理与初筛：将少量土壤样品悬浮于生理盐水中，取0.1毫升悬液涂布于以葡萄糖为唯一碳源的察氏琼脂（Czapek-Dox agar）平板。利用溴甲酚绿（bromocresol green）作为酸碱指示剂，通过观察平板颜色变化（产酸区域变黄）来初步筛选能够产生柠檬酸的真菌菌株。 2. 高产菌株筛选：从初筛的菌株中，通过定量分析，选择柠檬酸产量最高的菌株，命名为WU-2020。该菌株在马铃薯葡萄糖琼脂（PDA）上于25°C培养7天，进行菌落形态学观察。 3. 分子生物学鉴定：为了准确鉴定WU-2020的分类学地位，研究人员提取了其基因组DNA，并测定了多个基因的部分序列进行系统发育分析。所使用的基因标记包括核糖体DNA内转录间隔区（Internal Transcribed Spacer, ITS）、钙调蛋白（calmodulin）基因和β-微管蛋白（β-tubulin）基因。通过邻接法（Neighbor-Joining method）构建系统发育树（使用MEGA11软件），将WU-2020的序列与黑曲霉组内其他已知菌株的序列进行比较。结果显示，WU-2020与*Aspergillus lacticoffeatus*的模式菌株CBS 101883聚为一支，从而将其鉴定为Aspergillus lacticoffeatus。 4. 产酸性能初步验证：为了确认WU-2020在固态培养下的高产特性，研究团队将其与*Aspergillus lacticoffeatus*的模式菌株CBS 101883进行了平行比较实验。实验采用典型的半固态培养（semisolid culture）体系，以甘蔗渣（bagasse）作为载体培养基。培养4天后，测定葡萄糖消耗量和柠檬酸产量。结果显示，WU-2020消耗了120 g/L葡萄糖，生产了85.2 g/L柠檬酸，产酸得率（柠檬酸产量/葡萄糖消耗量）高达71%（重量比）。而模式菌株CBS 101883消耗了94.6 g/L葡萄糖，生产了40.4 g/L柠檬酸，得率为43%。WU-2020的柠檬酸产量是模式菌株的2.1倍，证明了其作为固态培养柠檬酸高产菌株的潜力。
第二环节：基因组DNA提取、文库构建与测序 在确定WU-2020为高产菌株后，研究进入基因组测序阶段。 1. 培养与DNA提取：将WU-2020菌株在马铃薯葡萄糖肉汤（Potato Dextrose Broth）中于25°C培养1天，收集菌丝体。使用酚-氯仿提取法结合Genomic-tip 20/G kit（Qiagen）从菌丝体中提取高质量的基因组DNA。 2. 测序文库构建与平台：为了获得高质量的基因组草图，研究采用了混合测序策略（Hybrid Sequencing），结合了短读长（Illumina）和长读长（Oxford Nanopore Technologies, ONT）技术。 * Illumina文库：使用NEB Next Ultra II FS DNA Library Prep Kit构建，在Illumina HiSeq X平台上进行双末端测序（Paired-End sequencing），读长为2×150 bp。测序服务由BGI（日本神户）完成。 * ONT文库：使用Short-read Eliminator Kit（Nippon Genetics）制备，在Oxford Nanopore平台上进行长读长测序。 3. 原始数据预处理： * Illumina短读长数据：使用fastp v0.20.1软件进行质控过滤，去除质量值（Q score）低于20或长度小于40个碱基的读段。最终获得20,713,137对高质量双末端读段，估算的基因组覆盖深度约为169倍。 * ONT长读长数据：使用ONT Guppy v4.0.14进行碱基识别（base calling），随后使用Porechop v0.2.4去除接头，并使用NanoFilt v2.7.1进行过滤。
第三环节：基因组组装、抛光与注释 这是本研究方法学的核心部分，涉及多个生物信息学工具和流程。 1. 混合组装：研究团队参考了Saud等人（2021）报道的方法，并进行了一些修改。具体流程如下： * 长读长校正与修剪：首先使用fmlrc软件对95,487条ONT长读长（总计2.05 Gb，平均长度21,419 bp）进行错误校正，然后使用Canu v1.9进行修剪。 * 初步组装：使用校正和修剪后的长读长，通过Flye v2.7软件进行从头组装（de novo assembly），生成初步的重叠群（contigs）。 * 多轮抛光：为了提升组装序列的准确性，进行了多轮抛光（polishing）： * 基于长读长的抛光：使用minimap2 v2.17将原始ONT读段比对到初步组装序列上，然后使用Racon v1.4.13进行一致性校正，最后使用Medaka v0.11.5进行进一步抛光。 * 基于短读长的抛光：使用bwa-mem2 v2.0将高质量的Illumina短读段比对到经过长读长抛光后的组装序列上，然后使用Pilon v1.23进行迭代校正，以修正剩余的碱基错误和小规模的插入/缺失。 2. 线粒体基因组组装：使用NOVOPlasty v4.2软件，从全基因组测序数据中独立组装出线粒体基因组。 3. 基因组注释：使用funannotate v1.7.4流程对组装好的核基因组进行注释。 * 基因预测：采用从头预测（ab initio prediction）方法，结合了多个预测工具：Augustus v3.3.2、GeneMark-ES v4.38、GlimmerHMM v3.0.2和SNAP。为了提高预测准确性，使用了来自黑曲霉（*Aspergillus niger*）和米曲霉（*Aspergillus oryzae*）的蛋白质序列以及UniProt数据库的信息作为外显子提示（exon hints）。 * 非编码RNA预测：使用tRNAscan-SE v2.0.5预测转运RNA（tRNA），使用Barrnap v0.9预测核糖体RNA（rRNA）。 * 线粒体基因组注释：使用GeSeq v2.03对组装的线粒体基因组进行注释。 * 基因完整性评估：使用BUSCO（Benchmarking Universal Single-Copy Orthologs）工具，基于散囊菌目（Eurotiales）数据库评估组装和注释的完整性。
第四环节：基因组特征分析与生物安全初步评估 基于组装和注释结果，研究人员对WU-2020的基因组进行了详细描述，并对其潜在的生物安全问题进行了初步分析。
主要研究结果 1. 菌株鉴定与表型：系统发育分析明确将WU-2020鉴定为*Aspergillus lacticoffeatus*。在半固态培养实验中，WU-2020表现出卓越的柠檬酸生产能力，产量（85.2 g/L）和得率（71%）均显著高于其模式菌株CBS 101883（40.4 g/L， 43%），证实了其作为固态培养柠檬酸高产菌株的应用潜力。 2. 核基因组特征：Aspergillus lacticoffeatus WU-2020的基因组草图（draft genome）由11个支架（scaffolds）组成，总长度为35.9 Mb，GC含量为49.51%。共预测到11,490个蛋白质编码基因、273个tRNA基因和80个rRNA基因。使用BUSCO基于散囊菌目数据库评估，完整性得分高达98.6%，表明该基因组组装和注释的质量非常高。 3. 线粒体基因组特征：线粒体基因组为环状DNA，长度为31,324 bp，包含15个蛋白编码基因、26个tRNA基因和2个rRNA基因。 4. 生物合成基因簇与生物安全警示：通过对基因组序列的分析，研究人员发现WU-2020菌株含有与赭曲霉毒素A（Ochratoxin A）和伏马菌素B2（Fumonisin B2）生物合成相关的推定基因簇（putative biosynthetic gene clusters）。这与该物种的模式菌株Aspergillus lacticoffeatus CBS 101883的报道情况相似。这一发现至关重要，它提示尽管WU-2020是柠檬酸高产菌株，但其具有产生这些有害霉菌毒素（mycotoxin）的遗传潜力。因此，论文明确指出，在考虑将其用于工业生产之前，必须进行进一步的实验来评估其实际产毒能力和生物安全性。
研究的结论与价值 本研究成功分离、鉴定并公布了适用于固态培养的柠檬酸高产菌株Aspergillus lacticoffeatus WU-2020的高质量基因组草图。该基因组数据为深入研究黑曲霉组真菌柠檬酸高产（尤其是固态发酵条件下）的分子机制提供了宝贵的遗传资源。同时，研究揭示了该菌株基因组中潜在的真菌毒素生物合成基因簇，为工业菌株的生物安全评估提供了关键的遗传学预警信息。这项工作不仅贡献了一个有价值的微生物基因组资源，也体现了在开发工业微生物资源时，结合功能性能（高产）与安全评估（毒素基因筛查）的现代研究范式。
研究的亮点 1. 目标菌株的特殊性：聚焦于适用于固态培养的柠檬酸高产菌株，填补了此前基因组研究多集中于深层发酵菌株的空白，对开发生物质废弃物固态发酵生产柠檬酸的工艺具有指导意义。 2. 高质量的混合组装：采用Illumina短读长与Oxford Nanopore长读长相结合的混合测序与组装策略，并经过多轮精细抛光，最终获得了高度连续（仅11个支架）且完整度高（BUSCO 98.6%）的基因组草图，数据质量可靠。 3. 完整的资源报告：研究不仅报告了核基因组，还提供了完整的线粒体基因组序列，并进行了详尽的基因注释，为后续比较基因组学和功能基因组学研究奠定了坚实基础。 4. 强调生物安全意识：研究并未只宣扬菌株的高产特性，而是基于基因组数据，主动揭示了其携带潜在霉菌毒素合成基因簇的风险，并明确指出了后续进行生物安全评估的必要性，体现了负责任的研究态度。
其他有价值的信息 * 菌种保藏：菌株WU-2020已保藏于两个公共资源中心：日本国家技术与评价研究所生物资源中心（NBRC，保藏号NBRC 115414）和千叶大学医学真菌学研究中心（通过国家生物资源项目，保藏号IFM 65579），方便其他研究者获取。 * 数据公开：本研究产生的所有原始测序数据和组装后的基因组序列均已公开。全基因组序列提交至DDBJ/ENA/GenBank，登录号为BQMC01000001至BQMC01000011；线粒体基因组登录号为LC670769；原始测序读段提交至DDBJ序列读段档案，登录号为DRA013207。这完全符合开放科学的原则。 * 方法学细节：论文详细描述了从DNA提取、文库制备、测序、数据过滤、混合组装、抛光到注释的完整生物信息学流程，并列出了所用每个软件的版本号，具有很高的可重复性和参考价值。