🧬 biology

Identification of Maize Fusarium Ear Rot Resistance Genes via Integrated BSA-Seq and RNA-Seq Analysis

本研究整合了 BSA-Seq、RNA-Seq 和 WGCNA 分析，以阐明玉米对禾谷镰刀菌（*Fusarium graminearum*）抗性的分子机制，鉴定了 8 个抗性相关位点和 6 个高置信度候选基因，这些基因涉及氧化应激、激素信号传导和细胞壁修饰，从而区分了抗性系 K0743 与感病系 K0742。

原作者： Jiaqin Liu, Wei Zhao, Fangguo Chang, Wei Wang, Fang Wang, Chun Wang, Yaqiong Shi

发布于 2026-07-01

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Jiaqin Liu, Wei Zhao, Fangguo Chang, Wei Wang, Fang Wang, Chun Wang, Yaqiong Shi

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ⚕️ 这是一篇未经同行评审的预印本的AI生成解释。这不是医疗建议。请勿根据此内容做出健康决定。阅读完整免责声明

核心概览：玉米地里的遭遇战

想象一下，一片玉米地就像一座繁忙的城市。有时，一个被称为“镰刀菌”（Fusarium graminearum）的邪恶入侵者（一种真菌）会袭击玉米穗。这不仅仅是令人烦恼的问题，它还会毁掉玉米，降低产量，并产生对人类和动物食用危险的隐形“毒气云”（霉菌毒素）。

科学家们想要弄清楚，为什么有些玉米植株在真菌攻击时像坚固的城堡（具有抗性），而另一些则像敞开的大门（易感）。他们研究了两种特定的玉米品系：

K0743（英雄）： 一种强韧、具有抗性的玉米品系。
K0742（受害者）： 一种脆弱、易感的玉米品系。

研究人员使用了一套“侦探工具包”，将传统生物学（观察细胞和化学物质）与高科技计算机科学（扫描植物的整个指令手册，即 DNA，以及其活跃的信息，即 RNA）结合在一起。

1. 物理战斗：崩塌的墙壁 vs. 坚固的城堡

类比： 把玉米粒想象成一座由淀粉砖块和蛋白质砂浆组成的房子。

受害者 (K0742)： 当真菌袭击时，房子迅速瓦解。到第3天时，“砖块”（淀粉）被侵蚀，“砂浆”消失了，真菌像征服者一样在空荡荡的房间里横冲直撞。到第7天时，房子已成废墟。
英雄 (K0743)： 这座房子守住了阵地。即使在真菌试图破门而入后，墙壁仍基本保持完整。“砖块”依然坚固，真菌无法立足。英雄玉米保持了更长时间的结构完整，有效地锁住了抵御入侵者的门。

2. 化学防御：警报系统与消防员

类比： 当小偷闯入时，一部智能房屋不会坐以待毙；它会触发警报并派出消防员。

研究人员测量了玉米体内的“防御化学物质”：

警报（酶）： 英雄玉米 (K0374) 发出的警报更快、更响亮。它产生了高水平的“防御酶”（如 PAL、PPO 和 SOD）。这些酶就像消防员，负责加固墙壁（使其更难被吞噬）并清理战斗中产生的有毒废物（氧化应激）。
燃料（渗透调节物质）： 英雄玉米还迅速储备了像脯氨酸和糖类这样的“应急燃料”。这让细胞保持水分和能量，以便持续战斗；相比之下，受害者玉米耗尽了能量并崩溃了。

3. 遗传侦探工作：寻找“抗性开关”

为了找出英雄玉米为何如此强韧，科学家使用了两种高科技方法：

方法 A：BSA-Seq（“寻找位置”的研究）

想象你有两座图书库（DNA）。一座属于英雄，另一座属于受害者。科学家将来自“超级英雄”组和“超级受害者”组的书籍混合成两大堆。

他们扫描这些书堆，查看哪些页面（基因）存在差异。
结果： 他们在玉米的指令手册（染色体）上找到了 8 个特定的位置，这些位置似乎控制着抗性。其中两个位置是“重磅选手”：
- qGER5（位于第 5 号染色体）
- qGER6-2（位于第 6 号染色体）
- 可以将它们视为玉米免疫系统的主要电路断路器。

方法 B：RNA-Seq（“寻找功能”的研究）

如果说 BSA-Seq 找到了开关的位置，那么 RNA-Seq 则告诉了他们这些开关究竟在做什么。他们观察了在感染期间哪些基因被“开启”（激活）了。

英雄的反应： 英雄玉米非常迅速地开启了数千个防御基因（在 72 小时时达到了 6,030 个！）。这就像是一场全面的军事动员。
受害者的反应： 受害者玉米反应迟钝，开启的基因数量极少（在 72 小时时仅为 805 个）。当它试图反击时，真菌早已取得了胜利。

4. 整合分析：“超级团队”

科学家将“位置”（BSA-Seq）和“功能”（RNA-Seq）数据结合起来。他们还使用了网络分析 (WGCNA) 来观察哪些基因在协同工作。

他们发现了 6 个“高置信度候选基因”，这些基因正好位于那两个主要的抗性区域（qGER5 和 qGER6-2）内。这些基因形成了两个特殊的团队：

“氧化应激”团队 (MEfloralwhite)： 这些基因充当清理小组。它们管理战斗中产生的有毒废物和压力，确保细胞不会因过度消耗而枯竭。
“信号与结构”团队 (MEblue4)： 这些基因充当将军和建筑队。它们发送“攻击指令”（激素信号）并加固玉米粒的物理墙壁。

结论

论文得出结论，英雄玉米 (K0743) 之所以能赢得战斗，是因为：

它具有更好的结构完整性（它的墙壁不会崩塌）。
它能通过大规模的防御基因进行更快速的反应。
它拥有 6 个特定的基因，通过完美的协调来管理压力并加固墙壁。

研究人员现在已经识别出了这 6 个特定基因以及它们在玉米 DNA 上的两个主要位置。这为育种家提供了一张“地图”，帮助他们寻找或培育未来能够抵御这种真菌、且无需依赖大量化学喷雾的玉米品种。

注：本文完全侧重于识别这些基因并理解其生物学机制。它并未声称这些基因已在商业作物中经过测试，也并未声称目前正用于种植用于销售的食物。

技术摘要：通过集成 BSA-Seq 与 RNA-Seq 分析鉴定玉米赤霉烯醇诱导的镰刀菌耳腐病抗性基因

问题陈述
玉米耳腐病主要由禾谷镰刀菌（Fusarium graminearum）引起，是一种破坏性极强的真菌病害，它会显著降低谷物产量和品质，同时引入威胁食品安全和动物健康的霉毒素（如 DON 和 ZEN）。虽然遗传抗性是最可持续的控制策略，但对禾谷镰刀菌的抗性是一种复杂的数量性状，受多个基因和环境因素的影响。以往的研究已经鉴定了许多数量性状位点（QTLs），但这些位点通常仅能解释一小部分表型变异，且在不同遗传背景下的稳定性不足。此外，玉米抗禾谷镰刀菌的具体分子机制和候选基因仍很大程度上不明，这阻碍了抗病品种的培育。

研究方法
本研究采用集成多组学方法，对比了人工接种禾谷镰刀菌后的抗病玉米自交系（K0743）与感病玉米自交系（K0742）。实验设计包括：

表型与细胞学分析： 通过内核伤口注射进行人工接种。根据 1–9 级评分进行发病程度评估。使用扫描电子显微镜（SEM）观察接种后 0、24、48、72、96 小时及 7 天时籽粒的超微结构变化。
生理与生化测定： 在多个时间点测量防御酶活性（PAL、PPO、SOD）和渗透调节物质（脯氨酸、可溶性糖、可溶性蛋白），以评估生理防御反应。
BSA-Seq（群体分离分析测序）： 使用由 K0743 和 K0742 杂交产生的 F₂ 群体。将极端抗性和感病群体（各 30 株）进行混合并进行 DNA 测序。利用深度学习（DL）和 $\Delta$ SNP-index 算法对抗性相关位点进行定位。
RNA-Seq（转录组分析）： 对两种品系的籽粒在接种后 0、72、96 小时及 7 天进行测序。鉴定差异表达基因（DEGs），随后进行基因本体（GO）和 KEGG 通路富集分析。
加权基因共表达网络分析（WGCNA）： 用于构建基因共表达模块并鉴定与抗病性状相关的枢纽基因（hub genes）。
集成候选基因鉴定： 通过将 BSA-Seq 的 QTL 区间与 RNA-Seq 的 DEGs 以及 WGNCA 的枢纽基因进行重叠分析，从而优先筛选高置信度的候选基因。
验证： 对九个选定的 DEGs 进行 RT-qPCR 实验，以验证转录组数据。

主要结果

生理与细胞学响应： 与感病系 K0742 相比，抗性系 K0743 保持了更完整的细胞结构和淀粉粒；而 K0742 在 7 天时表现出严重的组织降解和菌丝定殖。K0743 的防御酶（PAL、PPO、SOD）活性显著更高，且渗透调节剂（脯氨酸、可溶性糖、可溶性蛋白）的积累更快、更强劲。
BSA-Seq 定位： 在第 5、6、7、9 和 10 号染色体上鉴定了八个抗性相关位点。基于信号强度和重复性，两个主要位点 qGER5（第 5 染色体）和 qGER6-2（第 6 染色体）被优先确定。qGER5 同时被 DL 和 $\Delta$ SNP-index 方法检测到。
转录组动态： 抗性系 K0743 表现出更动态的转录响应，在所有时间点的 DEGs 数量均显著高于 K0742（例如，在 72 小时时，K0743 有 6,030 个 DEGs，而 K0742 仅有 805 个）。GO 和 KEGG 分析显示，K0743 特异性地富集了防御相关通路，包括“对生物刺激的反应”、“细胞壁组织”、“植物-病原体相互作用”和“MAPK 信号传导”，而 K0742 则主要富集于初级代谢过程。
WGCNA 与候选基因： 鉴定了两个关键的共表达模块：MEfloralwhite（富集于氧化应激反应和泛素介导的蛋白质降解）和 MEblue4（富集于植物激素信号转导、MAPK 信号传导和细胞壁修饰）。
集成候选基因： 通过取 BSA-Seq QTL、RNA-Seq DEGs 和 WGCNA 枢纽基因的交集，鉴定了六个高置信度候选基因：
- Zm00001eb242770 和 Zm00001eb244910（位于 qGER5 中）。
- Zm00001eb273840、Zm00001eb276160、Zm00001eb269250 和 Zm00001eb269150（位于 qGER6-2 中）。
  这些基因被注释为参与铁硫簇生物合成、生长素响应、生长调节、半乳糖变构酶活性、ARM 重复超家族功能以及 $\alpha/\beta$ -水解酶活性。其表达模式已通过 RT-qPCR 得到验证。

意义与主张
作者声称，本研究通过桥接细胞学、生理学和基因组学数据，为理解玉米抗镰刀菌耳腐病的分子机制提供了新见解。BSA-Seq 与 RNA-Seq 的集成结合 WGCNA，成功将抗性基因的搜索范围缩小至位于主要 QTL 区间内的特定六个高置信度候选基因。研究结果表明，K0743 的抗性是由一个涉及氧化应激管理（MEfloralwhite 模块）和结构/信号防御（MEblue4 模块）的协调网络所驱动的。作者认为，这些鉴定的位点和基因为抗性育种和功能特征化提供了宝贵的资源，解决了识别复杂数量性状中稳定抗性位点的挑战。该研究并未声称已完全阐明完整的调控网络，而是提供了一个用于未来验证的优先候选列表。