拟轮枝镰孢与玉米籽粒互作的差异基因筛选及代谢通路分析  

Screening of Differential Genes and Analysis of Metabolic Pathways in the Interaction Between Fusarium verticillioides and Maize Kernels

作　　者：渠清 刘宁 邹金鹏[2,3] 张雅璇 贾慧 孙蔓莉[2,3] 曹志艳 董金皋 QU Qing;LIU Ning;ZOU JinPeng;ZHANG YaXuan;JIA Hui;SUN ManLi;CAO ZhiYan;DONG JinGao(College of Life Sciences,Hebei Agricultural University,Baoding 071001,Hebei;Hebei Key Laboratory of Plant Physiology and Molecular Pathology/State Key Laboratory of North China Crop Improvement and Regulation,Baoding 071001,Hebei;College of Plant Protection,Hebei Agricultural University,Baoding 071001,Hebei)

机构地区：[1]河北农业大学生命科学学院,河北保定071001 [2]河北省植物生理与分子病理学重点实验室/华北作物改良与调控国家重点实验室,河北保定071001 [3]河北农业大学植物保护学院,河北保定071001

出　　处：《中国农业科学》2023年第6期1086-1101,共16页Scientia Agricultura Sinica

基　　金：国家玉米产业技术体系(CARS-02);华北作物改良与调控国家重点实验室自主科研项目(NCCIR2020ZZ-12);河北省重点研发计划(20326502D)。

摘　　要：【目的】拟轮枝镰孢(Fusarium verticillioides)引起的玉米穗腐病是我国玉米产区发生严重的病害之一,论文旨在了解病原菌与玉米籽粒互作过程中的基因表达差异,为揭示病原菌的致病机制和玉米抗病机制提供依据。【方法】对拟轮枝镰孢侵染玉米籽粒0、4、12和72 h的样品进行转录组测序,之后采用生物信息学分析,分别以玉米和拟轮枝镰孢基因组为参考,以|log_(2)FC|≥1,P-adjust<0.05为阈值筛选互作过程中玉米和拟轮枝镰孢的差异表达基因,利用GO和KEGG对其进行功能注释及富集分析。Goatools软件分析植物-病原互作、MAPK途径和植物激素信号转导通路相关差异基因的表达变化,采用实时荧光定量PCR(qRT-PCR)方法对测序差异基因进行验证。【结果】在互作4、12和72 h后拟轮枝镰孢分别有140、400和1945个基因上调表达,有9、302和1784个基因下调表达;玉米分别有293、692和1426个基因上调表达,320、482和153个基因下调表达。GO和KEGG富集分析显示,侵染早期拟轮枝镰孢在细胞间隙生长,差异基因主要富集在RNA生物合成、细胞壁结构成分、脂肪酸生物合成、蛋白质代谢、碳水化合物代谢、生物过程和代谢过程等通路中。拟轮枝镰孢早期侵染触发了玉米活性氧(ROS)爆发,差异基因主要富集在对活性氧、过氧化氢的反应,几丁质酶、单加氧酶活性,木质素代谢过程等相关通路中。侵染后期拟轮枝镰孢继续在籽粒中定殖及扩展,差异基因主要富集在碳水化合物和细胞壁多糖分解代谢过程、跨膜转运、氧化还原酶活性等功能通路中。玉米主要通过苯丙素、木质素、类黄酮生物合成,MAPK信号通路、植物-病原互作和植物激素信号转导等途径差异基因大量表达响应病原菌侵染。随机选取6个玉米和6个拟轮枝镰孢的差异表达基因进行qRT-PCR分析,基因表达规律与转录组测序结果一致,证实了RNA-seq的准�【Objective】Maize ear rot caused by Fusarium verticillioides is one of the most serious diseases in maize producing areas in China.The objective of this study is to understand the differences in gene expression during the plant-pathogen interaction at different stages,and to provide a basis for pathogenic mechanism of the pathogen infection and resistance mechanism of maize.【Method】Illumina platform was used to sequence the transcriptome of maize kernels infected with F.verticillioides at 0,4,12,and 72 h.The differentially expressed genes(DEGs)of maize and F.verticillioides were screened with|log_(2)FC|≥1,P-adjust<0.05 as threshold and clean reads were compared with genome of maize and F.verticillioides,separately.Functional annotation and enrichment analysis of DEGs were carried out by using GO and KEGG databases.Goatools software was used to analyze the expression changes of genes related to plant-pathogen interaction,MAPK signaling pathway and plant hormone signal transduction pathway.Sequencing results were verified by quantitative real-time PCR(qRT-PCR).【Result】A total of 140,400 and 1945 DEGs were up-regulated and 9,302,and 1784 DEGs were down-regulated in F.verticillioides after 4,12 and 72 h interaction,respectively.A total of 293,692,and 1426 DEGs were up-regulated and 320,482,and 153 DEGs were down-regulated in maize after 4,12 and 72 h interaction,respectively.GO and KEGG enrichment analysis of DEGs showed that F.verticillioides grew in intercellular space at the early stage of pathogen infection.The DEGs were enriched in RNA biosynthesis,cell wall structural component,fatty acid biosynthesis,protein metabolism,carbohydrate metabolism,biological process,and metabolic process.Reactive oxygen species(ROS)was triggered in maize at the early stage of infection.The DEGs were enriched in response to ROS,hydrogen peroxide,chitinase activity,monooxygenase activity,lignin metabolism.At the later stage of infection,F.verticillioides colonized and expanded in maize,and the DEGs were enriched in carbo

关 键 词：拟轮枝镰孢 玉米穗腐病 转录组 植物-病原互作 基因表达 差异表达基因 

分 类 号：S435.131[农业科学—农业昆虫与害虫防治]