长读长基因组测序面包小麦促进抗病基因克隆

Long-read genome sequencing of bread wheat facilitates disease resistance gene cloning

作者信息Naveenkumar Athiyannan, Michael Abrouk, Willem H P Boshoff, Stéphane Cauet, Nathalie Rodde, David Kudrna, Nahed Mohammed, Jan Bettgenhaeuser, Kirsty S Botha, Shannon S Derman, Rod A Wing, Renée Prins, Simon G Krattinger
PMID35288708
期刊Nat Genet
发布时间2022-03
DOI10.1038/s41588-022-01022-1

实验完整度

实验包括基因组组装、遗传定位、EMS突变体分析、VIGS沉默验证、表达分析、单倍型分析和蛋白建模等多个独立证据层级。

主要模型

面包小麦品种 Kariega 近等基因系 Avocet 2B Avocet S × Avocet 2B F2群体

重点核对

基因组组装:HiFi reads覆盖率约34倍,使用hifiasm组装,光学图谱和Omni-C挂载 遗传定位:Avocet 2B × Avocet S F2群体(n=345),定位到2BS染色体1.4 cM区间 EMS诱变:0.6% EMS处理Avocet 2B种子,获得10个独立突变体 VIGS沉默:使用BSMV载体,三个沉默片段(Si1, Si2, Si3)靶向Yr27 序列验证:Sanger测序确认EMS突变,QYr.sgi-2B.1编码序列在Kariega和Avocet 2B中一致

摘要

克隆大型复杂作物基因组中的农艺重要基因仍然具有挑战性。这里,我们结合高保真长读长、光学作图和染色体构象捕获,生成了南非面包小麦(Triticum aestivum)品种Kariega的14.7千兆碱基染色体规模组装。所得组装比先前的小麦组装连续性高一个数量级。Kariega对毁灭性的真菌条锈病表现出持久抗性¹。我们确定了小种特异性抗病基因Yr27,其编码一种细胞内免疫受体,是该抗性的主要贡献者。Yr27与叶锈病抗性基因Lr13等位;Yr27和Lr13蛋白显示出97%的序列同一性²,³。我们的结果证明了生成染色体规模小麦组装以克隆基因的可行性,并例证了单拷贝基因的高度相似等位基因可以赋予对不同病原体的抗性,这可能为将来工程化具有多种识别特异性的Yr27等位基因提供基础。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
如何克隆面包小麦品种Kariega中持久抗条锈病数量性状位点QYr.sgi-2B.1的基因?
核心机制
Yr27编码的NLR免疫受体识别小种特异性效应因子,触发抗性反应。Yr27与叶锈病抗性基因Lr13为等位基因,两者编码蛋白高度相似,但差异氨基酸集中在对识别特异性重要的区域。
主要证据
遗传定位将QYr.sgi-2B.1定位到2BS染色体,候选区间内有两个NLR基因;10个独立EMS突变体在该基因上均有突变;VIGS沉默该基因降低了抗性。
研究意义
研究证明了从任何小麦品系生成染色体规模组装以指导基因克隆的可行性,并暗示高度相似的等位基因可以赋予对不同病原体的抗性,为将来工程化具有多重识别特异性的NLR受体提供了基础。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

基因组组装与注释

获得Kariega品种的高质量参考基因组,用于基因克隆。

利用PacBio HiFi长读长、Bionano光学图谱和Omni-C Hi-C数据进行从头组装,并结合RNA-seq和Iso-seq进行基因注释。

2

遗传定位

确定QYr.sgi-2B.1在染色体上的位置。

利用Avocet 2B × Avocet S的F2群体进行连锁分析,使用侧翼标记和KASP标记定位QTL区间。

3

候选基因鉴定

在物理区间内筛选引起抗性的候选基因。

锚定标记后提取区间内的注释基因和NLR预测基因,通过序列分析和表达数据筛选候选。

4

EMS突变体筛选验证

通过突变体确认候选基因的功能。

对Avocet 2B种子进行EMS诱变,筛选丧失抗性的突变体,并测序验证候选基因是否发生突变。

5

VIGS沉默验证

通过基因沉默验证候选基因的功能。

构建BSMV载体沉默候选基因,观察抗性表型变化。

6

基因结构和表达分析

确定Yr27基因的结构和组织表达模式。

利用Iso-seq数据确认基因结构;通过qRT-PCR分析接种后的表达变化。

7

序列对比和单倍型分析

比较Yr27与Lr13及其他小麦品系中的序列差异。

比较Yr27和Lr13蛋白序列,分析不同品系中的Yr27等位基因差异,计算Shannon熵。

8

蛋白结构预测与比较

预测Yr27蛋白三维结构,分析关键氨基酸的作用。

使用Phyre2进行同源建模,并与RPP1和ZAR1结构进行叠合。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
PacBio Sequel II系统----
SMRTbell Express模板制备试剂盒2.0Pacific Biosciences100-938-900
Saphyr系统Bionano Genomics--
Saphyr芯片G1.2Bionano Genomics--
Dovetail Omni-C试剂盒Dovetail--
Maxwell RSC植物RNA试剂盒PromegaAS1500
NEBNext单细胞/低起始量cDNA合成与扩增模块New England BiolabsE6421S
ProNex尺寸选择纯化系统PromegaNG2001
Sequel II结合试剂盒2.1Pacific Biosciences--
内部对照试剂盒1.0Pacific Biosciences101-843-000
测序试剂盒2.0Pacific Biosciences101-820-200
SMRT Cell 8MPacific Biosciences101-389-001
FC-43油3MFluorinert
EMSSigma-AldrichM0880
GoTaq Green预混液PromegaM7122
Phusion高保真DNA聚合酶New England BiolabsM0530
KASP预混液LGCKBS-1016-016
High-Capacity cDNA逆转录试剂盒Applied Biosystems4368814
PowerUp SYBR Green预混液Applied BiosystemsAS25741
QuantStudio 6 Flex实时荧光定量PCR系统ABI--
Femto Pulse系统Agilent--
Qubit dsDNA高灵敏度检测Invitrogen--
NanoDrop NP-1000NanoDrop Technologies--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
基因组组装
HiFi reads覆盖度(约34倍)、组装软件(hifiasm v0.11)及参数、光学图谱构建步骤、Omni-C文库制备试剂盒及测序平台。
阅读提示:Methods中'HiFi library preparation and sequencing'、'Optical map production'、'Omni-C library preparation and sequencing'、'Genome assembly and validation'
遗传定位
F2群体大小(n=345)、侧翼标记(Xbarc55和Xwmc344)、KASP标记、连锁分析参数(LOD≥3.0)、Kosambi作图函数。
阅读提示:Methods中'Genetic map'
EMS诱变
EMS浓度(0.6%)、处理时间(16小时)、种子数量(约2,000)、筛选世代(M1-M3)。
阅读提示:Methods中'EMS mutagenesis'
VIGS沉默
载体(pBS-BSMV-γ)、沉默片段序列及长度(Si1, Si2, Si3)、对照设计(GFP和GUS)、农杆菌菌株(GV3101)、接种方法及表型观察时间。
阅读提示:Methods中'VIGS'
表达分析
qRT-PCR引物(2BNLR5F17/R11)、内参基因(Ta.6863)、生物学重复数(5)、接种后取样时间点。
阅读提示:Methods中'Quantitative real-time PCR'