《Journal of Fungi》：The Complete Mitogenome of Alternaria brassicae Causing Black Spot Disease on Rapeseed-Mustard
编辑推荐：
芸薹生链格孢（Alternaria brassicae）是一种破坏性真菌病原菌，引起油菜-芥菜及其他十字花科作物黑斑病，导致全球范围内显著产量损失。尽管其具有经济重要性，A. brassicae的线粒体基因组此前未见报道。研究人员利用HiMT工具包对公开可用的
芸薹生链格孢（Alternaria brassicae）是一种破坏性真菌病原菌，引起油菜-芥菜及其他十字花科作物黑斑病，导致全球范围内显著产量损失。尽管其具有经济重要性，A. brassicae的线粒体基因组此前未见报道。研究人员利用HiMT工具包对公开可用的PacBio HiFi读段进行de novo组装，获得一条环状线粒体DNA分子，全长50,713 bp，平均测序深度309.0×。该基因组编码19个蛋白质编码基因（protein-coding genes, PCGs），包括12个具名线粒体基因（atp6、cox1、cox3、cob、nad1、nad2、nad3、nad4、nad4L、nad5、nad6和rps3）及7个内含子开放阅读框（intronic open reading frames, iORFs）（lagli_0–lagli_5和giy），2个rRNA基因（rrnL、rrnS）和37个tRNA基因，总体GC含量29.3%。在nad4、nad4L、nad5、cox3、cob和cox1基因中预测到含内含子结构；这些预测未经转录本数据验证。共检测33个简单序列重复（simple sequence repeat, SSR）位点，以A/T单核苷酸重复为主，密度为一处SSR每1.52 kb。串联重复主要聚集于33–40 kb和45–50 kb两区；散在重复分析显示倒位重复占45.83%、互补重复占30.28%、正向重复占12.22%、回文重复占11.67%。基于22个线粒体基因组的极大似然系统发育分析将A. brassicae（CNSA: CNP0008733）置于A. burnsii的姐妹位，自举支持率65%；其与A. gossypina更深层聚类的支持率为25%。相较于已发表的A. tenuissima线粒体基因组（57,475 bp、34个tRNA、2个rRNA），本次组装短6762 bp且多3个tRNA。该首个人A. brassicae完整线粒体基因组资源为后续群体遗传学、分子诊断与进化动态研究奠定基础，助力十字花科作物黑斑病改良管理策略开发。
研究背景方面，十字花科作物（含油菜Brassica napus、白菜Brassica rapa ssp. pekinensis、甘蓝Brassica oleracea var. capitata及萝卜Raphanus sativus）是全球食用油、膳食纤维与必需营养素的重要来源。芸薹生链格孢（Alternaria brassicae (Berk.) Sacc.，座囊菌纲Dothideomycetes、格孢腔菌目Pleosporales、格孢腔菌科Pleosporaceae）作为坏死营养型真菌，是黑斑病（亦称链格孢叶斑病）首要病原，可造成油菜-芥菜产量损失15%至50%以上。该菌经气传或种传分生孢子萌发、附着胞形成及气孔或直接穿透侵染，致叶片、茎秆与角果出现褐黑坏死斑伴黄晕，引发早衰、落叶、光合下降与角果开裂。当前化学防治（QoI类、DMI类杀菌剂）与耕作措施难以遏制病害，商用种质缺乏高水平抗病性且面临杀菌剂抗药性进化风险。线粒体为真核细胞呼吸产ATP、ROS稳态与程序性死亡关键细胞器，真菌线粒体基因组（mitogenome）编码氧化磷酸化（OXPHOS）电子传递链核心亚基及rRNA/tRNA，其架构可塑性（大小、共线性、内含子、重复）使其成为进化与种群标记；且Complex III为QoI类靶标，解析线粒体基因结构对监测靶标突变具直接意义。Alternaria属核基因组研究虽进展，但A. brassicae完整线粒体基因组长期缺位，构成细胞器进化与呼吸链认知缺口。研究人员据此利用三代PacBio HiFi长读段，首次完成A. brassicae线粒体基因组de novo组装、注释、重复解析及基于22分类单元的系统发育重建，论文发表于《Journal of Fungi》。
关键技术方法上，研究人员从NCBI SRA获取公开PacBio HiFi数据（SRR32051149，源自油菜黑斑病斑分离物），经FastQC质检与Filtlong过滤（Q≥20、长度≥1000 bp）；以HiMT v1.0基于OLC算法从全基因组长读段中分选线粒体读段并环化组装；MITOS2与MFannot初注，BLASTp比对curated 12个PCGs，tRNAscan-SE识tRNA，RNAmmer联合BLASTn定rRNA，按同源比对划内含子边界并CDD注LAGLIDADG/GIY-YIG；MISA与TRF析SSR/串联重复，BIBIServ2 RepeatFinder析散在重复；取12个共有PCGs经MAFFT(L-INS-i)氨基酸比对、TrimAl修剪、回译密码子比对，ModelTest-NG选GTR+G+I模型，MEGA11极大似然法（1000自举）建树。
研究结果部分，3.1节“A. brassicae线粒体基因组的测序、组装与基本特征”中，HiFi干净读段1,385,133条回贴确认环状50,713 bp序列（CNP0008733），平均覆盖309.0×无断裂；注释得19 PCGs（12具名+7 iORFs: lagli_0–lagli_5、giy）、2 rRNA（rrnL、rrnS）、37 tRNA，GC 29.3%；33 SSR（1157 bp，密度1/1.52 kb，单核苷酸占57.6%且A/T占单核苷酸94.7%），最长单体重复140 bp（45,171–45,310），串联重复富集33–40 kb与45–50 kb，23,604–23,873存132 bp二倍重复；散在重复中倒位165、互补109、正向44、回文42；nad4、nad4L、nad5、cox3、cob、cox1预测含内含子（未获转录验证），Cob第143位为甘氨酸且其后无Group I内含子。3.2节“系统发育与比较线粒体基因组分析”中，22分类单元（7 Alternaria、13 Pleosporales、2 Ophiostomatales外类群，Acrocalymma sp. PP151286生根）12 PCGs串联矩阵ML树显示A. brassicae姐妹于A. burnsii（PV822416，自举65%），再与A. gossypina聚合并获25%支持；较A. tenuissima（57,475 bp、34 tRNA）短6762 bp且多3 tRNA，rRNA数同，本研究未做基因序比较故不议共线性。
讨论部分指出，A. brassicae线粒体基因组（50,713 bp，GC 29.3%）较A. alternata（60.0 kb）、A. solani（58.1 kb）更紧凑，大小差异由非编码间隔与内含子动态扩张驱动而非核心编码改变；33 SSR与360散在重复为滑链错配多态性候选标记，其是否介导分子内重组待验；六呼吸基因内含子预测符真菌移动元件特征但未转录验证；Cob G143及cox1后无Group I内含子为靶标突变监测基线而非表型敏感佐证；线粒体ML树A. brassicae姐妹A. burnsii仅中等支持（65%），可与核基因组/多位点分类互参。结论部分译文如下：本研究报道首例芸薹生链格孢完整线粒体基因组。环状基因组长50,713 bp，编码19个蛋白质编码基因、37个tRNA与2个rRNA。系统发育分析厘清该物种在链格孢属内分类位置并揭示其线粒体基因组结构特征。上述结果可为A. brassicae群体遗传学、分子诊断及杀菌剂抗性相关变异研究提供基础基因组资源。
