第一作者联系方式: 李闯, E-mail: chuangli163@yeah.net; 刘承晨, E-mail: 731064949@qq.com
云南是世界公认的亚洲栽培稻遗传多样性和起源中心之一。利用51对SSR分子标记对云南哈尼梯田当前栽培的111份水稻材料进行多态性检测, 并分析其群体结构。聚类分析表明, 供试材料主要分为偏粳类群(group I)及偏籼类群(group II)两大类群, 其中以偏籼类群居多, 占83%, 这两类群并未完全按地理来源聚类。所选材料稻米的糊化温度变异较广, 可以分为低(<66℃)、中低(66~70℃)、中高(70~74℃)和高(>74℃) 4类。通过基因测序, 分析了稻米糊化温度控制基因 ALK的序列多样性, 发现其可以分为10种单倍型。关联分析结果表明, 所选水稻样品稻米的糊化温度主要由 ALK基因3个SNP位点组合控制。其中, A-GC和G-TT两种组合主要控制低(<66℃)和中低(66~70℃)糊化温度, G-GC组合主要控制中高(70~74℃)和高(>74℃)糊化温度。SNP3402 T类型与高、中高糊化温度有关。绝大多数偏籼类品种的 ALK基因含有G-GC组合, 因此都倾向于高糊化温度。上述研究结果为水稻分子育种和种质资源保护研究提供了一定的参考价值。
Yunnan province is one of the centers for genetic diversity of cultivated rice in Asia. A total of 111 rice landraces were collected from Hani’s rice terraces of Yunnan, and their population structure was analyzed by using 51 SSR markers. The 111 rice landraces were divided into japonicagroup (group I) and indicagroup (group II); however, each of the two groups was not completely clustered according to its geographical origin. Most collected landraces, accounting for 83%, belonged to indica group. Gelatinization temperature (GT) is an important parameter affecting rice cooking and eating quality. The GT of rice flour was measured, which was widely different among the collected rice varieties, showing four classes in very low ( <66°C), low (66-70°C), high (70-74°C) and very high ( >74°C) GT. Meanwhile, ALK gene, also known as SSSIIa and controlling GT, was sequenced and compared, showing ten haploid types of the ALKgene among tested samples. Association analysis, revealed that GT is mainly controlled by the combination of three SNPs of ALK gene. Among them, the A-GC or G-TT combination type controls the very low and low GT, while the combination of G-GC controls high and very high GT. SNP3402 type T is associated with high and very high GT. It was also noticed that most of the varieties from indicagroup tended to contain the G-GC mutation and thus showed a relative high GT compared with those from the japonicagroup. These results are meaningful for rice molecular breeding and germplasm resources protection.
水稻(Oryza sativa L.)是世界上最重要粮食作物, 也是适应性最广的栽培作物之一[1]。云南是世界上公认的亚洲栽培稻遗传多样性和起源中心之一, 拥有中国全部野生稻种、籼粳两个栽培稻亚种以及六大生态群, 是中国稻种资源最大的遗传和生态多样性中心及天然宝库[2, 3]。因此, 利用现代分子生物学技术手段研究云南当前栽培水稻品种, 对于保护稻种资源多样性以及栽培稻的遗传改良具有重要的意义[4]。云南哈尼梯田分布于红河州内哀牢山区红河南岸元阳、红河、绿春和金平四县, 总面积5.47万公顷。其中, 位于北纬23° 04° ~23° 12° 、东经102° 40° ~ 102° 48° 的红河州元阳县中部地区的梯田面积广大、分布集中, 代表性最强。
近年来, 利用SSR分子标记分析云南水稻群体遗传结构的报道较多。董超等[5]利用48个SSR标记分析了哈尼梯田的98份月亮谷种质, 将其分成5个居群, 这5个居群与其村寨间地理位置的距离相关。董树斌等[6]利用24个SSR标记对云南660份样本遗传多样性和遗传关系的分析, 表明云南水稻传统品种具有较高水平的遗传多样性。我们对收集自云南哈尼梯田的47份栽培水稻的分析证明, 可将其分为偏籼和偏粳两大类群; 偏籼类群又可分为籼稻亚群、偏籼糯稻亚群及红米亚群; 偏粳类群又可分为粳稻亚群及偏粳糯稻亚群[1]。
稻米的食用方法多种多样, 对其品质的要求随用途而异[7]。糊化温度(gelatinization temperature, GT)是评价稻米蒸煮与食味品质的重要指标之一[8]。已有研究表明, 糊化温度受胚乳淀粉结构如支链淀粉链长分布、不同生育期环境温度等因素的影响[9]。支链淀粉的合成受众多酶催化, 其中可溶性淀粉合酶(soluble starch synthase, SSS)是较重要的一种。水稻中的SSS可分为SSSI、SSSII、SSSIII和SSSIV 4类, 而SSSII又含有SSSIIa (SSSII-3)、SSSIIb (SSSII-2)和SSSIIc (SSSII-1) 3个同工型[10]。其中, SSSIIa主要负责延伸支链淀粉的短支链(A链+B1链), 同时还负责合成中等长度的葡聚糖进而形成中等长度的分支(B2+B3链), 从而影响淀粉颗粒的形态结构[11]。该酶的编码基因ALK (又称为SSSIIa)是控制稻米糊化温度的主效基因[12]。
本研究以近年来从云南哈尼梯田直接采集的111份栽培水稻品种为材料, 选取均匀覆盖水稻全基因组的51对SSR标记进行遗传多样性和群体结构分析。在此基础上, 对111份水稻种质中ALK基因全长序列测序并分析其核酸多态性, 以期探析云南哈尼梯田当前栽培水稻中ALK基因多态性及其对稻米糊化温度的影响。
分别于2012年9月底和2013年10月上旬分两次在云南省元阳县10个代表性村寨的自然梯田内采集水稻成熟种子, 共111份(附表1)。这10个共村寨分别是下新城、小新街、爱春村、全福庄、箐口村、黄草岭、大鱼塘、阿党寨、核桃寨及上主鲁村。为进行比较分析, 在试验同时选择8份常规水稻品种作为对照, 包括南粳46、武运29154、武运粳8号和日本晴4份粳稻品种, 以及明恢63、龙特甫、黄花占和9311共4份籼稻品种。2013年11月将所有种质种于海南省陵水县扬州大学南繁基地, 于2014年4月收获成熟种子供品质分析用。
参照CTAB法[13]从水稻叶片中提取基因组总DNA, 20° C下保存备用。参考刘承晨等[1]的标记, 选取均匀分布于水稻全基因组且在9311和日本晴均有多态性的51对SSR标记。按常规PCR方法进行SSR检测, 在扩增图谱上与9311迁移率相同的记为“ 1/1” , 与日本晴相同的记为“ 2/2” , 如果两者都不同则重新编号(如“ 3/3” ), 无带记为“ ?/?” 。利用PowerMarker V 3.25软件按非加权配对法(Unweighted Pair-Group Method With Arithmetic Means, UPGMA)及Nei 1983遗传距离对SSR检测结果聚类, 并在MEGA 6.0中编辑。
以粳稻品种日本晴ALK基因序列(GenBank登录号为NC_029261.1)为参考, 选取从起始密码子ATG到终止密码子TGA共4422 bp的序列。利用引物设计软件Primer 5.0设计ALK-a1/ALK-a2、ALK-b1/ALK-b2、ALK-c1/ALK-c2和ALK-d1/ALK-d2 (表1) 4对引物扩增目的基因, 把目的基因分为4段(a、b、c、d)分段扩增。PCR程序为95° C预变性2 min; 98° C, 变性12 s; 49~54° C (各个片段扩增时的温度稍有差异), 退火10 s; 72° C, 延伸2 min, 35个循环; 最后72° C延伸5 min; 18° C保温。PCR产物经1%的琼脂糖凝胶电泳分离, 然后在Bio-Rad Μ V 1000核酸成像仪上观察并切胶送南京金斯瑞生物科技有限公司(http://www.genscript.com.cn/)测序。
水稻种子完全成熟后混收, 自然晾干2个月后, 用砻谷机(Model SY88-TH, 韩国)出糙, Kett 精米机(Tokyo, 日本)磨成精米, FOSS旋风式磨粉机(FOSS, 瑞典)磨成粉。过100目筛后于40° C烘箱烘干2 d, 室温下平衡2 d后装入小塑封袋中密封。利用DSC 200 F3差示扫描量热仪(德国耐驰公司)测量供试米粉的热力学特性, 参照Yan等[14]的方法用配套的分析软件分析其特征参数米粉糊化时的热焓(Δ H)变化、起始温度(To)及峰值温度(Tp)等。
根据与8份籼粳稻对照品种聚类关系的远近, 可将供试材料分为偏粳(group I)和偏籼(group II)两大类群(图1)。其中偏粳类群有19份, 占总数的17%, 又可分为4个亚类。偏籼类群92份, 占供试材料的83%, 又可分为7个亚类。上述聚类结果说明111份采集自哈尼梯田的水稻材料尽管在地理位置上很集中, 但在遗传上存在很大的分化, 以籼型为主, 但有明显的籼粳稻分化。这可能与其独特的气候和纬度条件有关。
图2表明, 所测定材料稻米起始糊化温度To变异范围为62~77° C, 其中以71~73° C间的居多。按起始糊化温度To, 可分为4档, 即低糊化温度( < 66° C)、中低糊化温度(66~70° C)、中高糊化温度(70~74° C)和高糊化温度( > 74° C)。糊化的峰值温度Tp变异范围在66~81° C之间, 以75~78° C间的分布较多, 与To分布规律较为相似。热焓值Δ H在5.5~ 12.0 J g-1之间, 且呈近似正态分布, 平均值为8.21 J g-1。糊化温度特征值参数Δ H和Tp在0.05水平上呈正显著性相关, 相关系数为0.243; Tp和To在0.01水平上呈正显著性相关, 相关系数达0.947。
在111份云南稻中共发现32处变异位点, 包括2处InDel (Insert/Delet)类型和30个SNP类型(表2)。32个变异位点中, 9个SNP位于外显子区域, 其中6个为有义突变, 包括: (1)第1外显子185 G/A突变位点导致一个精氨酸到组氨酸的变化; (2)第1外显子264 C/G突变位点导致一个天冬氨酸到谷氨酸的变化; (3)第2外显子613 G/C突变位点导致一个天冬氨酸到组氨酸的变化; (4)第8外显子3799 A/G突变位点导致一个丝氨酸到甘氨酸的变化; (5)第8外显子4198 G/A突变位点导致蛋氨酸到缬氨酸的变化; (6)第8外显子4330 C/T突变位点导致一个亮氨酸到苯丙氨酸变化。
据多态性分析结果, 可将供试材料中ALK基因分为10种单倍型(Hp1~Hp10)(图3), 其单倍型多样性为0.748, 单倍型多样性方差为0.000 59, 单倍型多样性标准差为0.024, 平均核酸多态性π 为0.000 59, 多态性参值θ w为0.001 28。
对ALK基因不同单倍型的分析显示, 第8外显子内的3个SNP (4198 G/A、4329 G/T和4330 C/T)呈一定规律的组合形式, 在111份云南稻中共出现A-GC、G-GC和G-TT三种组合类型。其中, 主要以G-GC类型为主(占72%), G-TT与A-GC各占总材料的14%。含G-GC类型的水稻以籼型为主, 其中仅有4份材料属于粳稻(占G-GC类型总数的5%); G-TT类型以籼稻为主, 占75%, A-GC类型以粳稻为主, 约占73%。
供试云南栽培水稻稻米糊化温度的变异主要是与第8外显子中3个SNP (4198 G/A、4329 G/T和4330 C/T)有关, 其中与4198 G/A位点呈极显著关联(P < 0.0001, R2= 0.211), 4198处为碱基A的类型与低和中低糊化温度有关; 而其余2个SNP位点4329G/T和4330C/T也达显著关联(P < 0.01, R2= 0.089)。由这3个SNP组成的3种组合类型(A-GC、G-GC和G-TT)中, A-GC和G-TT两种组合主要控制低( < 66° C)和中低(66~70° C)糊化温度, 而G-GC组合则主要控制中高(70~74° C)和高( > 74° C)糊化温度。此外, 第7内含子3402 C/T对稻米糊化温度的变异也有显著关联(P < 0.01, R2= 0.093), 3402 T类型与高和中高糊化温度有关(表2和表3)。
由ALK基因的10种单倍型与糊化温度间的相关分析可知(表3), Hp2与对照粳稻品种日本晴属同一单倍型, 控制低糊化温度; 单倍型Hp1、Hp3、Hp8和Hp10主要控制中高和高档糊化温度; 其余5种单倍型主要控制低和中低档糊化温度。从群体结构分析发现, 偏籼类群以高和中高糊化温度为主(G-GC和G-TT类型), G-GC类型中95%的材料属于籼稻, G-TT类型中75%为籼稻。偏粳类群以低和中低糊化温度为主(A-GC类型), A-GC类型中约73%为粳稻(图3和表3)。
SSR标记分析表明, 本研究所检测的云南栽培水稻的遗传多样性低于前人报道的。供试材料主要分为偏粳(group I)及偏籼(group II)两大类群, 但并未完全按地理来源聚类, 其中以偏籼类群居多(占83%)。偏粳类群(group I)糊化温度中等, 偏籼类群(group II)糊化温度较高。说明在主要理化品质方面当地人对籼稻具有一定偏好性。分析其原因可能是: (1)由于自然条件等的影响, 对栽培水稻的选择具有一定偏向性; (2)本研究供试材料来源于元阳中部具有代表性的哈尼族人村寨, 因此所选材料可能与哈尼族人的喜好有一定的关系。总体而言, 元阳哈尼梯田现有栽培水稻的遗传多样性较之前的报道有所降低, 这些水稻的主要农艺性状和品质性状均具有趋同的趋势, 这可能是杂交水稻等育成品种的推广、当地农民的选择偏好、以及传统文化习俗的淡化所致。从自然资源保护角度出发, 建议进一步加强对云南元阳哈尼梯田地方水稻品种和传统稻作文化的保护力度。
稻米的蒸煮与食味品质主要由淀粉的物化特性决定[19]。已有研究表明, 糊化温度并非单基因控制, 其他种类淀粉合酶基因也会影响糊化温度的高低[20, 21, 23]。由群体籼粳分化与ALK基因的系统进化关系结果可知, 单倍型Hp1、Hp5、Hp8、Hp10全部来自籼稻品种; Hp6、Hp7、Hp9全部来自粳稻品种; Hp3、Hp4中多数为籼稻, 分别占95%和80%; Hp2多数为粳稻, 约占73%。从群体结构与品质比较来看, 偏籼类群稻米糊化温度总体上明显高于偏粳类群。由ALK基因10种单倍型与稻米糊化温度的关联分析可见, 单倍型Hp1、Hp3、Hp8、Hp10主要控制中高和高档糊化温度, 其余单倍型主要控制低和中低档糊化温度。Hp1单倍型的第7内含子3402C/T位点呈显著关联(P< 0.01), 3402T类型与高、中高糊化温度有关, 可能为新的单倍型, 未曾见有报道; 在Hp2 (NIP型)中检测到4198 A/G位点与糊化温度也呈极显著关联(P< 0.0001), 4198A类型与低和中低糊化温度有关。
本研究供试材料中的ALK基因共有32个序列变异位点, 其中外显子中有9个。Bao等就曾在30份水稻品种的ALK序列2705~4755位点间的2051 bp区间内就发现25个SNP, 其中9个SNP位于外显子区域。高振宇等[17]曾报道ALK基因在籼稻双科早和粳稻日本晴内含子中存在19个SNP位点, 外显子中发现4个SNP位点。本研究共发现有义突变6个, 其中SNP185和SNP613未曾有人报道, 其余位点均已被报道。可能功能位点SNP4198、SNP4329、SNP4330前人均有报道, 而SNP3402位点前人未曾报道。前人报道的可能功能位点SNP264在本研究中也检测到, 但与糊化温度并未达到显著关联, 与前人报道有所不同[7]。从本研究对所选种质ALK基因序列变异及其与糊化温度的关联分析来看, 糊化温度主要由该基因3个SNP位点控制。含A-GC和G-TT两种序列组合的稻米表现为低( < 66° C)或中低(66~70° C)糊化温度, 而含G-GC序列组合的水稻则具有中高(70~74° C)或高( > 74° C)糊化温度, 这些结果与前人报道相似。综上, 通过测序和关联分析, 本研究鉴定了ALK基因的10种单倍型, 为进一步揭示该基因的遗传变异及育种利用提供了有用的遗传信息, 但关于各单倍型的具体功能尚需通过构建近等基因系或转基因系等来具体地验证分析。
云南哈尼梯田当前栽培水稻可分为偏粳(group I)及偏籼(group II)两大类群, 其中以偏籼类群居多, 占83%; 稻米糊化温度变异较广, 其控制基因ALK可分为10种单倍型。所选水稻样品稻米糊化温度主要由ALK基因3个SNP位点组合控制, 其中A-GC和G-TT两种组合主要控制低(< 66° C)和中低(66~70° C)糊化温度, G-GC组合主要控制中高(70~74° C)和高(> 74° C)糊化温度。绝大多数偏籼类品种的ALK基因含有G-GC组合, 因此都倾向于高糊化温度。
The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|