水稻香味的起源和进化

1、水稻香味的起源和进化在水稻作物中，香味是评介其品质的一个重要参数，虽然甜菜碱醛脱氢酶基因（BADH2）的起源和进化在这个特征下还不活楚。我们假定广泛存在8个BADH2勺非功能等位基因表明，这些基因具有不同的地理和遗传起源。尽管香味特征有多个起源，单一 等位基因badh2.1在几乎所有香稻品种占主要地位， 包括巴斯马蒂和茉莉花。单体型分 析，使我们能够建立一个粳稻品种群 badh2.1等位基因单一来源并使这个等位基因从粳 稻品质渗入釉稻品种中。巴斯马蒂样的种类，不论其香味基因的表型，都与粳稻原始种 类一样，跨越了 BADH观翼的5.5MB区域，展示了巴斯马蒂和粳稻基因库品种之间密切 的进化关系。

2、这些结果表明了水稻品种之间香味的关系，并对认为水稻香味特点起源丁 釉稻的传统假说提出了挑战。最近的研究表明，业洲栽培稻（Oryza sativa）是几个不同的遗传组组成，在这些遗 传组中，一些等位基因对关键驯化和粮食品质性状负责。由丁复杂的水稻进化史，这些 基因变异的起源是怎么样通过高度不同的水稻业种生存下来，在水稻生物进化中任然是个核心问题。香味被认为是最重要的粮食稻米品质性状之一，因为它是决定市场价格的关键因 素，同样关系到地方和国家的地位（1, 2）。对水稻香味基因的遗传基础进行调查研究， 发现了水稻8号染色体上的一个单一位点与水稻香味有关（3、4）。精细定位（5-7） 和随后的序列分析

3、鉴定甜菜碱醛脱氢酶基因，BADH2与香味表型有关命名如下（8）。把功能基因的突变创建隐性badh2.1等位基因描述为在第七外显子基因处三个核苜酸多 态性（SNP$和8 - bp的缺失，导致了密码子的提前终止以及假定存在缩短的BADH2S白（9）。其他序列比对已被用来描述这种复杂的突变（10, 11）。因此，badh2.1突变 被称为功能核苜酸多态性（FNP）。最近的调查各种不同的香味种质，支持 badh2.1与香 味有关（10, 12, 13）, 一个有关香味与非香味主要基因的转变实验，已经显示消除了香 味（14）,确认BADH2S水稻种是主要香味遗传决定基因。在香米品种中，100多种容易挥发

4、的物质已经被检测到，2-乙酰-1 -毗咯琳（2AP） （15, 16）是最主要的化合物对香米品质其主要作用。这种化合物在水稻植物除根以外 的所有器官中产生，含量非常低，人们能够立即在密集种植或者碎叶组织中发现，以及 在蒸煮过程或者以后或之前的谷物中发现（17）。然而2AP的化学合成途径还没有成功，BADH那化r-氨基丁醛的氧化（AB-ald; a 2AP的前身），丁是非功能等位基因在AB-ald 和其循环物中毗咯琳积累，导致 2AP的合成的增强。中国自古已经认识到（19, 20）业洲栽培稻（O.sativa ）分为两个主要品种群体，釉 稻和粳稻（谈到品种组时用大写）。用15种同工酶标记能 够进

5、一步的把这两个主要群体 分为遗传上的六个不同的业群，这和公认生态型是一致的（21）.随后用SSRs（22） andSNPs （23）标记分辨出了和上面的群体差不多 5种分界线活楚的遗传群组：即釉稻，澳 大利业稻，温带粳稻，热带粳稻和芳香稻（提到业群时小写）（图1）。系统进化分析法和Fst显示构成粳稻品种群的芳香稻，温带粳稻和热带粳业群之间有密切的进化关系； 然而釉稻，澳大利业稻业群具有不同的祖先，并作为釉稻品种组成员（22, 23）。尽管它的名字怎么样，芳香业群的表型多样，既包括香型也包括非香味型。为了避免混淆， 从今以后我们将用其同工酶的名字命名芳香业群，简称V组（24）。水稻香料新品种已经

6、确定至少包含三个不同的业群，包括v组（即巴斯马蒂和'沙德里品种），釉稻（即茉莉品种）与热带粳稻。在调查各种不同的新型香米 中，几乎都含有badh2.1等位基因，证明这个等位基因是在不同香稻中普遍遗传的（9, 10, 11, 12）。最近,BADH2S因的第二个突变位点 badh2.2,被证明与来自中国的种 质资源有有限的联系。有证据表明有额外的突变可能导致 2AP的合成来自一个严格的研 究，包括主要来自东南业的已经检测到几个新香料种质资源多样化小组，这已使2AP水平升高，但没有携带badh2.1等位基因（12）。考虑到单个基因主要是在负责水稻的香 味，这项研究的目的主要是在基因不同的水

7、稻业群之间研究这个等位基因的起源和探究 它的祖先。我们同样也把BADH舛主要负责独立遗传，香味基因定位的基因的额外功能 突变作为研究对象。结论badh2.1等位基因在不同的水稻种质资源中出现的概率。我们曾在280种野生稻（普通野生稻）中观察badh2.1等位基因的出现频率，发现它不在所有野生基因型中，除了 一个单一位点杂合的等位基因外（表 1）。这种野生插入位点显示了驯化水稻的几个特 征，包括白色果皮，表明这是栽培稻香味品种和野生稻近交的结果。共收集调查了 176 种不同的栽培稻种类，这些栽培稻的不同业种的身份发现是通过一套全基因组SSR和SNP分子标记确定的（22,23）。总之，badh2.

8、1等位基因在这些种类中被检测到17次（10%），香味基因在v组中出现的频率最大，在温带稻和澳大利业稻中出现的频率最小（表1 ） 0badh2.1等位基因的起源。考虑到在釉稻还是在粳稻品种群体中检测相同badh2.1等位基因，我们的目标是检测 badh2.1等位基因起源丁哪个种类。为了解决这个问题， 我们检测了一个包含242个水稻品种的的BADH鑫基因，包含起源组和拥有BADH2S因 的额外加入组（表s1）。在5kb的对齐基因系列中，我们发现106个SNPW记的系列插 入缺失和SSR多态性，其中54个已经显示出频率大丁 5%.。这些多态性被用丁构建 8 个基因单体型（GH ,以及这些单体基因型被

9、明确的分成两个不同的小组（图 2a ,表 s2）。在第一组中，粳稻品种组（Jap_GH）的所有品种都携带了野生类型的等位基因，然 而绝大多数（74%的种类在第二组中都是来自丁釉稻。尽管业种的种类不同，每一个 插入系列都在Jap_GH佥测到了 badh2.1基因，在Jap_GFfr,香味种类基因仅仅在编码 badh2.1等位基因的FNP中和原始非香味基因不同（在图2a中用黄色显示的部分）。这 些数据支持badh2.1等位基因在像粳稻型遗传背景的品种中有单一的起源。来自釉稻品种组携带有badh2.1等位基因的的所有香稻品种与 Jap-GH形成集群， 形成一个这些品种基因和基因组表型系统之间明显的矛

10、盾。如果釉稻品种能检测到具有釉稻遗传背景环绕BADH殊因一样的粳稻DN制定区域，这种不一致就能够得到很好的 理解。然后，在我们这个242个业洲栽培稻品种的BADH染因的跨度上游3.2MB和下游 2.1MB的基因组区域扩增24个位点系列，用以扩宽我们对单基因表型的分析。在这个 BADH观翼的5.3MB的范围内，总的有426个SNPW记的插入缺失，在一致的13KB区域 内确定了 SSR标记的多态性。其中271个多态性的出现率大丁 5%从这些数据中，78 个多态原始信息（AIPs）被找到并且用丁扩增单体型（见方法和材料，表 s3）。这些扩增单 体型被称为单体扩增一致系列显示出六类大的单体型（图2b）

11、。这些扩增的单体型和基因单体型是一致的，所以含有badh2.1基因插入系列都有BADH基因周围粳稻的DNAT 增区域。在釉稻品种中有 24个携带等位基因badh2.1，粳稻区域被上游650KB和下游 330KB的关键重组位点被扩展，其侧翼区和釉稻的原始扩展单基因表型是一致的（图2B）。 这支持了 badh2.1等位基因是经过基因渗入釉稻种类的假说。第v组种类携带badh2.1等位基因扩展单倍型和原始粳稻品种一样，含有一个扩展 的单基因型，除了 FNPft在14、30、69三个位点上的三个特殊多态型（图 2B中绿色标 注部分）。这三个FNPM翼的多态位点以及像粳稻品种那样的区域，都被检测到在所有

12、 的釉稻插入位点携带badh2.1等位基因。来自缅甸单一野生型种类基因 badh2.1在后代 是1:2:1的杂合表型，并且单基因表型分析确认了 携带有badh2.1等位基因的染色体同样含有上述v组的三个多态型。因此，我们能够在野生型杂合品种上追踪badh2.1等位基因的祖先和所有釉稻香味品种来自 v组的原始类型。有几个非香味粳稻品种和Ind_GH集群形成团体，这些是由釉稻非香味品种构成（图 2A）。这些来自粳稻的非香味品种包含了像釉稻基因组一样包围 BADH2S因的区域，并 且所有显示重组都回到粳稻原始类型的侧翼区（图 2B,粳稻重组）。这为一定数目的包 含了 BADH凛因的粳稻品种和Ind_

13、GH形成的集群提供了很好的解释。在badh2.1等位基因上，核酸多样性的减少和连锁不平衡的提高。为检查围绕badh2.1基因的选择证据，我们分析了这个 242个品种的BADH2S因的核苜酸多样性。 我们对比携带野生型基因和bdha2.1等位基因的5kb插入系列之间的BADH2S因的排列 系列的核酸多样性。携带香味基因的 badh2.1等位基因和非香味基因品种相比，在多样 性方面平均降低了 97.4% （表2）。扩展单倍型纯合子（EHH评估了这个可能性：两个随机选择的基因组区域的遗传 是相同的，允许功能突变导致遗传距离的增加的连锁不平衡（LD）衰变测量法（25）。在我们的种子资源组，通过比较香味

14、基因（badh2.1）和非香味基因品种的扩展单体型，我们为BADH2S因计算了 EHH携带香味基因（badh2.1）的品种显示了大范围围绕 FNP 的扩展LD,然而LD围绕的野生型等位基因减少很快（图 3）。当LD的程度在包含野生 型或者badh2.1等位基因的各个业种中对比时，观察到了相同的模式（图s1）。badh2.1, 核酸多样性显著减少，加上围绕这个等位基因的连锁不平衡区域的扩宽，为badh2.1 FNP. 的强阳性选择提供了证据。BADH2付加突变基因。BADH2付加突变基因有可能导致2AP的合成，BADH2S因在 26个香味基因品种中，这些系列没有任何先前所描述派生赋予香味等位基因

15、（12）。先前242个品种中的26个系列比对显示8个非同义替换的多态型，其中 4个移码诱导缺 失多态和其中一个是SNP成提前终止密码子，所有这些在缩短了的BADH2I白质中都 是公认了的（图4A,表s4）。其他三个潜在功能多态性，包括3 - bp的插入，编码区 两个不同SNPs.虽然其中的几个多态性仅仅在一个单一插入位点发现并确认，在香味基因插入位点找到4个，这些多态性看上去有某些地理学上的联系（图4B）。迄今为止，我们还不能确定两个插入位点候选功能突变能够为2AP的水平提高做出解释。讨论 水稻香味基因的起源。 从来自业洲各地不同的水稻种质资源大范围的遗传和地理 上，这项研究展现了 BADH基

16、因多样性的深度研究。携带badh2.1等位基因的香味插入基因和非香味基因插入系列相比，显示出核酸多样性的急剧减少以及连锁不平衡的提 升，这和badh2.1等位基因强阳性选择是一致的。badh2.1等位基因选择强度和这个报 道水稻基因控制谷物多态性特点的报道的相似（如wx, rc, gs3 ） （26-28）。单体型分析，使我们能够证明badh2.1等位基因在粳稻品种群体遗传背景产生。扩 展单体型分析显示粳稻基因组区域的明确基因渗入，在所有釉稻香味基因插入系列中够 包含badh2.1等位基因，包括茉莉花品种。所有香釉稻同样拥有被确定为v组业群的badh2.1等位基因侧翼的三个多态型，表明这些重要

17、的谷物性状等位基因起源丁v组业群（即Basmati巴斯马蒂）并被渗入釉稻品种（即 Jasmine，茉莉花）。在69位点多态性的起源以及香味插入 V组都是原始的和固定的，所有热带粳稻插 入位点都携带badh2.1等位基因。这些表明：（1） badh2.1等位基因可能来自v组和热 带粳稻的后代，第五组的两个特定基因多态性（14和30位点，图2B）来自第五组后代 的分离，并且基因组区域含有badh2.1等位基因和下游多态位点（69）可能是从v组渗 入到热带粳稻中。相同的序列延伸区域两侧BADH基因在V组和热带粳稻插入之间阻止 了断点检测重组，这将定义一个从 v组到热带粳稻的渗入体系（反之亦然）。从本

18、研究的证据表明，业洲栽培稻从原始祖先引种驯化后，badh2.1等位基因被选定为一个从头突变，也许粳稻业群以后的分化会引起等位基因频繁出现在v组以及温带粳稻上难以发现的低频率。在我们收集的280不同的插入系列的来自业洲的澳普通野生 稻/澳野生稻，我们发现仅仅一个来自缅甸的单一位点在badh2.1等位基因是杂合的。这些插入位点的系歹0分析显示，在v组香味基因品种染色体携带 badh2.1等位基因在釉 稻DNAff列边缘包含3个显著多态性固定位点。然而badh2.1等位基因被确定在其他来 自东南业的水稻种质中（29-31）,很可能这些有代表性的原始等位基因渗入野生品种 中并且在栽培稻的附近微弱的生长

19、。等位基因从业洲栽培稻侵入野生品种已经多次报道 和现在认为是一个普遍存在的现象（32.33）。这项研究改进了我们对业洲栽培稻种类之间基因组联系的理解和进一步明确业洲 栽培稻的进化历程。来自v组的巴斯马蒂品种经常被错误的认为是釉稻品种的一个成员 （34-36）。然而来自釉稻和v组的品种广泛的种植在业洲南部和中部，并有可能延长， 其形态是细长谷物粒，研究人员早就注意了这两个品种之间高度的杂种不育现象（37-40）。Glaszmann通过同工酶标记显示为杂交不亲合和提出了一个假说，V组和其他釉稻品种是明显不同的，分组和粳稻品种更接近（21.24）。业洲栽培稻的基因多样性附加检测是使用叶绿体标记，核酸

20、SSR函SNPs立决定，V组业群是一个和粳稻品种关系很近的独特的遗传实体(22,23,41 )。我们对这项研究的单体型分析表明V组品种，香型和非香型，和原始粳稻插入系列形成集群，两者都跨度整个BADHS因和8号染色体5.3MB的被研究区域。尽管它的形态相似和地理分布与釉稻重叠，在基因上是粳 稻种类的一个成员，这为V组提供了进一步的证据。有趣的是，在南业 V组和釉稻品种 的重叠分布区域可能为主要的香味基因进入釉稻品种提供了一个通道。香味水稻的独立起源。水稻品种显示的是提高了 2AP的位置，但是缺少任何已知 的BADH2勺非功能等位基因，为BADH*的额外香味等位基因的存在提供了可能性 (12)。

21、 对26个缺少香味品种的分析，在预计改变 BADH；®白的编码区域找到了 8个多态性。 就像前面所述，只有全长的BADH2羊本，才能产生完整的503-aa蛋白质，能够抑制2AP 的生产功能(140).在这个研究中确定的8个多态位点的4个(等位基因badh2.3 - 2.6 ), 被推测为引起转录提前终止，这将推定去除蛋白质功能和导致香味基因的不复存在。这 些所有突变导致 BADH2®白质在形成催化和/或者底物结合区域的关键残基之前缩短 (14)。此外，badh2.7等位基因，被6个澳大利业品种共有，同样被推测为消除蛋白 质聚集区域复制体的缩短。其他 3个位点(2.8-2.1

22、0 )，是额外的氨基酸位点(2.8) 或氨基酸替代位点(2.9-2.10 )。尽管这些香味基因和BADH2S因突变位点的联系，进 一步的工作是确定这些突变对 AP积累程度的影响。BADH笏位基因上品种之间携带相同等位基因突变点的地理连系表明香味基因的选 择在不同的地理区域、多种情况下是独立的。有趣的是，我们现在发现对高2AP生产率负责的主要等位基因，在水稻是来自丁粳稻不同种类，这看起来好像其他的香味基因是 来自丁其他的釉稻种类一样,比如 badh2.7等位基因，他是在几个澳大利业品种中找到 的。这些以及其他的结果表明，V组和澳大利业系列拥有有用的在水稻改良方面未得到 充分利用的等位基因(22,

23、 27, 47, 48)。26个香味品种缺少任何一个先前提到的香味识别基因，我们能够在编码或者启动区域找到两个品种没有突变，预计将改变 BADH2S白质或其表达。因此，这些插入位点可 能包括其他地方的遗传病变点在代谢途径上控制 2AP的合成；映射实验正在进行用以检 测这个假说。负责香味的额外基因的识别能够为改变水稻基因的品质适应不同文化的需 要提供机遇。水稻BADH基因更广范围的驯化和分化。在水稻种已经有记载一个在基因方面与水 稻驯化和谷物质量改进相关的多功能等位基因。在其中一些等位基因位点似乎是从直立 型野生种选择的，然而其他的是再次从栽培类型基因库中选取的(42.43 )。无论他们起源丁哪

24、里，这些等位基因在业洲栽培稻品种之间要么是独立的要么变得广泛播种。被釉 稻和粳稻同时享有派生等位基因许多特点在一个业群中显示单一等位基因来源，符合渗入过程。RC等位基因，在现种业洲栽培稻中98%!一个白色果皮基因，作为一个驯化相 关基因的例子，在粳稻品种遗传背景中传到其他种类时形成，即基因不同业群问形成(27)。在釉稻不同品种中，RC等位基因渗入的平均大小是 1MB和在这个研究中所述 的badh2.1等位基因渐渗大小差不多。同样地，在 GS3B因突变点拥有长粒谷物来源的 釉稻祖先，随后渗入到其他栽培稻中(26)。在这个研究中，我们发现了 badh2.1等位 基因的进化历史也届丁这个类型。 这些

25、等位基因通过渐渗杂交的过程能在不同的水稻业 群之间转移，在早期栽培稻高异原杂交率和物理距离的不同的栽培稻之间得到简化，结 果是人口的膨胀和人类的迁徙(42, 44-46)。值得注意的是，尽管明显的重要性杂交和 基因问的移动，反力量还是维持了不同业洲种之间的遗传分歧。解决这一悖论，需要对 业群隔离的催成关键因素作进一步的研究，同样更进一步对不同组之间的遗传交换动力 学的研究提供方法。材料和方法植物材料。我们的种质库乂 280种澳普通野生稻和来自38个国家的242中业洲栽 培稻组成。我们同样从先前的研究中获得了 26个缺少badh2.1等位基因香味品种(12)。 在表S1中是这个研究中用到的插入系

26、列的全部名单。所有未经报到过的业洲栽培稻都 是用SSRsS记检测从而确定他们的业群。DNA提取，PCR和测序。DNA提取是利用醋酸甲-SDS处理叶组织和处理精种子。 badh2.1的标记功能是用种质基因遗传型。对丁基因单体型分析，8个位点约700bp的扩增物测序是BADHS因的编码区域，形成5kb的排列系列。延长单体型分析，24个地 区进行测序，跨度BADHS因上游3.2MB和下游2.1MB范围，形成超过13kb的排列系 列。先前描述MITE多态性在在栽培稻同样是遗传的。表 5是32个引物序列的全部。康 奈尔大学生命科学中心的核心实验室 PCB物是在ABI Prism 3700/3100 DN

27、A分析仪上 纯化和测序的。序列排序是用 CodonCode®位仪，末端扩增产物被调整以除去低质量序 歹0。单独和不确定序歹U在必要的情况下被重新测序。在缺少任何已知的非功能等位基因 的BADH2®入位点中找到非同义多态，为确定这些他们被重新测序很多次。2AP表型。用改进的方法二氯甲烷萃取2AP (52).化学合成2AP是由日本Dr.T.Yoshihashi提供的(日本茨城，日本国际研究中心的农业科学)以及用来量化 2AP 样品。每个样品抽提和分析是在六个不同的场合以及至少三个复制的情况下进行。BADH密的单体型和遗传多样性分析。在242个业洲栽培稻插入位点BADH染因编 码

28、区域的八个扩增产物排序被导入 TASSE旺序提取所有多态型，频率在 5恕上，在构 建的单体型基因样本中（表si）.在20249280位点的高度多态性SSRgE没有延长也没有 缩短（TA11 - 14） （TA6- 8）,这和粳稻品种和釉稻品种分别是一致的，如果等位基因比8个重复序列多则标记“ 1”，少丁这标记“ 0”。总的8个基因单体型是从238个品种推导 出的，其中4个重组单基因型没有包含在图2A中。Bayesian集群结构被用来分析单基 因型和在2个集群中得到的最高的可能性，用丁标记粳稻单基因型集群（Jap_GH）和釉稻单基因集群（Ind_GH）对丁扩展单体型分析，我们在 242个业洲栽培

29、稻中的BADH殊因侧翼测定了 24区 域，如上所述（表S5）。TASSEIM用来提取5%以上抽样本216多态性的频率。为了减 少单倍型数，我们采用了以下标准选择 AIPs:粳稻品种和釉稻品种之间必须具有显著的 差异（P < 0.00003），在这个品种中出现频率高丁 20%来自16个扩展序列的78个多态 型中，FNP MITE多态性达到这个标准，见图 2B、和表S3中。用前述的方法（25）, 242个业洲栽培稻的BADH观翼的24个位点区域的多态性数 据被用来计算EHH图1 :业洲栽培稻业群结构。来自169个核SSRsB记的无根连接树。枝条颜色代 表叶绿体单基因表型。自助值（共100个）

30、是表示在分支点，这树活楚地说明了品种群 之间的两个主要分支（釉稻和粳稻），这进一步细分为5个水稻业群：釉稻，澳大利业， 热带粳稻，粳稻温带和 V组（芳香）。经许可复制和修改。图2 : BADH她区单倍型基因分析。（A） BADH基因单体型。序列读取242个业洲 栽培稻的整个BADH凛因，所有多态性（频率>5%）为申联，用丁创建8个基因的单倍 型。每个单倍型字母都代表 SNP位点的核苜酸选择性。数字表明了一个删除的大小 （0 为没有删除），每个多态性的相对位置表示沿 BADHS因模型，表型的单体型的编号1 至8和香味型是对应的以及来自拥有单体型的业群的插入数。两个基因单体型集群确定：粳稻基因单体型的群集（Jap_GH）和釉稻基因单体型集群（Ind_GH）.。蓝细胞代表Jap_GH多态性特点，红色代表Ind_GH细胞的特性。badh2.1 FNP描绘灰色/黄色，黄色代表香 味等位基因。（B）扩展单倍型。共延长了 17种单倍型为BADH2 6个一致单体型被每个显示表型 描述。字母或数字以及颜色所指都和 A图一样有（以下是简略缺失：1x_28,1y_12,1z_48）。问 断着色表明重组在哪里被发现。第51位，标有星号，代表前面所述的MITE多态位点（10）. badh2.1 FNP在第52位用深黄色标记出来，“ + “是指香味起源基因