日本晴/中嘉早17重组自交系产量性状QTL定位

删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

本站小编 Free考研考试/2021-12-26

张应洲, 罗荣剑, 圣忠华, 焦桂爱, 唐绍清, 胡培松, 魏祥进. 日本晴/中嘉早17重组自交系产量性状QTL定位[J]. , 2017, 50(19): 3640-3651 https://doi.org/10.3864/j.issn.0578-1752.2017.19.002
ZHANG YingZhou, LUO RongJian, SHENG ZhongHua, JIAO GuiAi, TANG ShaoQing, HU PeiSong, WEI XiangJin. QTL Mapping of Yield Associated Traits of Nipponbare ×Zhongjiazao 17 RIL Population[J]. Scientia Acricultura Sinica, 2017, 50(19): 3640-3651 https://doi.org/10.3864/j.issn.0578-1752.2017.19.002

0 引言

【研究意义】水稻是世界重要粮食作物,世界上1/3以上的人口以水稻为主食,在中国水稻约占粮食播种面积与总产量的27%和35%^[1]。作为中国第一大口粮作物,水稻承担了中国约60%人口的口粮。因此,提高水稻产量一直都是水稻育种的重要目标^[2]。水稻产量主要由每株穗数、每穗粒数和千粒重等因素构成,这些性状都属于多基因控制的数量性状。因此,发掘控制水稻产量性状的优异资源并对水稻产量性状的遗传机制进行解析,将为高产水稻品种的选育提供重要的指导意义。【前人研究进展】随着分子遗传学的快速发展,水稻产量性状的数量性状位点（QTL）定位与克隆取得了巨大进步。据不完全统计^[3-4],目前,单株产量共检测到至少80个QTL,分布于全部12条染色体上,其中第1、7染色体上检测频次较高;有效穗数共检测到207个QTL,分布于全部12条染色体上,其中第1、2、4、7染色体上检测的频次较高;每穗粒数共检测到90多个QTL,分布于除第5染色体外的11条染色体上,其中第2、3、6、7、11染色体上检测的频次较高;结实率共检测到50多个QTL,分布于全部12条染色体上,其中第2、6染色体上检测的频次较高;千粒重共检测到160多个QTL,分布于全部12条染色体上,其中在第1、2、3、5、6、10、11染色体上检测的频次较高。水稻粒型很大程度上直接决定着千粒重,也是研究者们关注的重点。目前,已经定位了100多个粒长QTL、90多个粒宽QTL,这些QTL分布于12条染色体上,其中第1、3、5染色体上检测的频次较高^[5-6]。虽然水稻产量性状相关QTL数量众多,但是只有一些主效QTL被克隆。其中,主要包括调控水稻穗粒数的QTL/基因,如GN1A^[7]、TAW1^[8]、DEP1^[9]、OsSPL14^[10]、Ghd7^[11]、DTH8^[12]等;控制水稻粒型或粒重的QTL/基因,如GW2^[13]、qSW5^[14]、GW5^[15]、GS3^[16]、GS5^[17]、GS6^[18]、GW8^[19]、TGW6^[20]、qGL3^[21]和GL7^[22] 等。【本研究切入点】尽管在水稻产量相关性状分子机理方面已取得长足的进步,然而,由于产量相关性状的复杂性、受多基因调控以及易受环境影响的特点,水稻产量形成的遗传及调控机制仍不清楚,一些精细定位或已克隆的基因在育种实践中难以得到广泛运用。因此,更多的水稻产量相关QTL/基因的定位与克隆将为水稻高产分子育种提供更完备理论基础与基因资源。中嘉早17是2009年通过国家审定的超级早籼稻品种,由于其高产、稳产、适合加工等优点,2010—2017年连续8年被农业部推荐为全国水稻主导品种。目前,年推广面积超过66.67万公顷,已经成为近年来中国南方稻区推广面积最大的水稻品种。但是中嘉早17高产、稳产的遗传机理尚未明确。【拟解决的关键问题】以粳稻品种日本晴和中嘉早17为亲本衍生了一套RIL群体,为RIL群体构建了分子连锁图谱,同时考察RIL群体多年多环境的产量相关性状,检测产量相关性状QTL。以期阐明中嘉早17的高产机理,挖掘水稻新的产量相关QTL/基因,为超高产品种育种提供理论基础和基因资源。

1 材料与方法

1.1 研究材料

以粳稻日本晴为母本,超级稻中嘉早17为父本杂交产生F₁,F₁自交后得到F₂群体,以后各株系均利用单粒传法连续自交直至第10代,最终获得由200个株系组成的重组自交系（RIL）群体。RIL群体2015年7月种植于杭州富阳（2015HZ）;2016年1月种植于海南陵水（2016HN）;2016年6月种植于杭州富阳（2016HZ）。材料的种植均按亲本和群体每个株系种植64株（8行×8列）,株行间距为20 cm×20 cm,区块间距25 cm,单本栽插,田间管理同大田生产。

1.2 连锁图谱的构建

在分蘖盛期采集亲本及RIL群体各家系叶片,按CTAB法提取DNA^[23]。针对全基因组选用629对SSR引物（序列来自Gramene网站,由上海英骏生物公司合成）对亲本日本晴和中嘉早17进行多态性筛选。根据多态标记在物理图谱中的位置,挑选出均匀覆盖水稻12条染色体的多态性SSR标记,对RIL群体200个家系进行基因型检测和分析。通过PCR和聚丙烯酰胺凝胶电泳获得SSR标记的基因型分型数据。分别将中嘉早17和日本晴相似的带型标记为“1”和“2”、杂合带和模糊带型分别记为“3”和“0”。应用Mapmaker/EXP3.0对筛选的标记进行连锁分析^[24],以LOD值为3.0确定连锁群内标记顺序,采用Kosambi函数将重组值转换成遗传图距离（cM）,利用mapdraw 2.1绘制染色体,构建连锁图谱。

1.3 产量相关农艺性状考查

于成熟期,每个家系选取表型基本一致的5个单株（剔除边行）,参考MONCADA等^[25]提出的方法考查8个性状,分别为有效穗数（effective panicles,EP）、每穗粒数（number of grain per panicle,NGPP）、结实率（seed setting rate,SSR）、单株产量（yield per plant,YPP）、千粒重（1 000-grain-weight,TGW）、粒长（grain length,GL）、粒宽（grain width,GW）和粒厚（grain thickness,GT）。以5个单株平均数作为分析的基本数值。

1.4 QTL定位

采用Windows QTL Cartographer 2.5软件^[26],以复合区间作图法（CIM）,在全基因组范围内检测产量相关性状QTL。以Permutation 1 000次运算确定LOD阈值。当某个位置检测到的LOD值大于LOD阈值时,就认为该位置存在1个QTL。QTL命名法则遵循BLAIR等^[27]的原则。同时,基于MICM法,采用QTLNetwork2.2^[28]联合检测3个环境下产量相关性状QTL与环境的互作效应。

2 结果

2.1 日本晴×中嘉早17重组自交系的分子遗传连锁图谱

从629对SSR标记中共检测到246对在亲本间表现多态性,挑选出均匀覆盖水稻12条染色体的163对多态性标记,检测群体家系的基因型,获得163个标记的基因型数据。经卡方检验,73%标记符合1:1单基因理论分离比,27%标记显著偏离1:1分离比,发生偏分离的标记大部分偏向父本中嘉早17。利用163个SSR标记分析RIL群体基因型,其中有效带型（1,2）占总数的99.25%,杂合带和模糊带型（3,0）仅占0.75%。父本中嘉早17基因型分布频率见图1。在200个株系中,父本基因型频率变幅13%—83%,平均值为56%;在163个标记中,父本基因型频率为37%—85%,平均值为56%。不同个体和标记的基因型频率分布基本符合正态分布,表明该研究群体具有较好的遗传构成。应用Mapmaker/EXP3.0对163个多态性标记间进行连锁分析,构建覆盖12条染色体的连锁图谱（表1）。每条染色体连锁群标记数为9—22个,遗传距离变幅78.0—190.2 cM。遗传连锁图谱总图距为1 479.40 cM,覆盖整个基因组的93.85%,标记间平均图距为9.08 cM,变幅为0.8—26.7 cM,达到QTL作图的要求。
Table 1
表1
表1遗传图谱信息
Table 1The genetic linkage map information

染色体 Chr.	总长度 Total length (Mb)	图谱长度 Length of map		标记数 Number of markers	覆盖率 Coverage (%)	标记平均图距 Average interval (cM)
染色体 Chr.	总长度 Total length (Mb)	物理距离 Physical distance (Mb)	遗传距离 Genetic distance (cM)	标记数 Number of markers	覆盖率 Coverage (%)	标记平均图距 Average interval (cM)
1	45.05	43.20	177.30	18	95.89	9.85
2	36.78	35.71	134.80	18	97.09	7.49
3	37.37	36.02	190.20	22	96.39	8.65
4	36.15	34.92	125.00	14	96.60	8.93
5	30.00	29.43	160.00	16	98.10	10.00
6	31.60	30.77	84.60	11	97.37	7.69
7	30.28	28.91	113.40	12	95.48	9.45
8	28.57	27.67	127.00	14	96.85	9.07
9	30.53	22.29	78.00	9	73.01	8.67
10	23.96	22.75	89.00	10	94.95	8.90
11	30.76	28.02	95.40	9	91.09	10.60
12	27.77	25.22	104.70	10	90.82	10.47
总计Total	388.82	364.91	1479.40	163	93.64	9.08

新窗口打开

显示原图|下载原图ZIP|生成PPT
图1父本基因型频率的分布
-->Fig. 1Distributions of paternal genotype frequency in the RILs population
-->

2.2 群体及其亲本产量相关性状表型

对RIL群体及其亲本的8个产量相关农艺性状在3种环境下的表型值统计分析发现,RIL群体各性状表现明显分离,连续分布且有超亲分离现象（表2、图2）。产量相关性状在不同环境下有一定的差异,其中2015HZ和2016HZ的有效穗数平均值少于2016HN;2015HZ和2016HZ的每穗粒数、单株产量和结实率3个性状平均值明显高于2016HN;千粒重、粒长、粒宽、粒厚在3个环境下基本无差异,认为这4个性状受环境的影响较小。RIL群体各性状的峰值和偏值都较小,其中粒长和粒宽3个环境下峰值均>1,每穗粒数和结实率在2015HZ、2016HZ峰值和偏值均>1,而其他性状的峰值介于-0.50—0.55,偏值介于-0.64—0.99,基本符合数量性状遗传的正态分布特点。
Table 2
表2
表2RIL群体与亲本间水稻产量性状的表型分析
Table 2Agronomic traits performance data of RILs population and parents

性状 Trait	环境 Environment	重组自交系 RIL populations				亲本 Parent
性状 Trait	环境 Environment	变幅 Variation	平均值±标准差 Mean±SD	峰值 Kurtosis	偏值 Skewness	中嘉早17 Zhongjiazao 17	日本晴 Nipponbare
有效穗数EP	2015HZ	5.80—21.00	10.49±2.33	0.35	0.76	11.00	10.00
	2016HN	5.67—26.67	12.28 ±3.06	-0.40	0.42	11.33	11.67
	2016HZ	6.00—18.04	10.76±2.31	0.55	0.71	10.19	9.98
每穗粒数NGPP	2015HZ	35.00—254.90	117.21±44.35	-0.50	0.44	154.67	63.58
	2016HN	22.78—251.10	96.87±36.87	0.16	0.57	140.96	70.18
	2016HZ	37.91—262.61	122.99±44.49	-0.16	0.55	164.94	78.37
结实率SSR (%)	2015HZ	0.12—0.93	0.67±0.17	1.10	-1.06	0.75	0.74
	2016HN	0.08—0.91	0.61±0.20	-0.50	-0.64	0.80	0.52
	2016HZ	0.15—0.91	0.66±0.13	1.40	-1.05	0.79	0.66
单株产量YPP (g)	2015HZ	3.47—37.78	18.32±6.21	0.29	0.42	27.32	11.86
	2016HN	0.35—39.45	17.05±7.23	-0.35	0.41	26.12	16.34
	2016HZ	1.36—44.05	17.82±9.52	0.14	0.62	27.29	10.61
千粒重TGW (g)	2015HZ	16.92—34.20	23.71±2.74	0.99	0.73	26.57	25.26
	2016HN	15.92—32.95	23.47±2.57	-0.32	0.22	26.50	24.12
	2016HZ	13.38—32.04	22.73±2.90	0.47	0.15	26.18	24.92
粒长GL (mm)	2015HZ	7.00—9.98	8.11±0.47	2.30	0.99	7.99	7.64
	2016HN	6.96—9.76	7.92±0.49	1.51	0.82	8.03	7.47
	2016HZ	6.66—9.82	7.76±0.50	1.34	0.85	7.89	7.50
粒宽GW (mm)	2015HZ	2.47—3.78	2.99±0.21	1.33	0.42	2.99	3.09
	2016HN	2.25—3.79	2.98±0.21	1.25	0.33	3.09	3.02
	2016HZ	2.24—3.50	2.82±0.21	1.51	0.67	2.99	2.90
粒厚GT (mm)	2015HZ	1.83—2.50	2.12±0.11	-0.14	-0.19	2.24	2.10
	2016HN	1.81—2.37	2.08±0.11	-0.37	-0.18	2.15	2.10
	2016HZ	1.80—2.40	2.02±0.12	-0.17	-0.09	2.14	2.09

新窗口打开

显示原图|下载原图ZIP|生成PPT
图2RIL群体水稻产量相关性状的表型频率分布
虚线箭头指示日本晴表型平均值;实线箭头指示中嘉早17表型平均值
-->Fig. 2Phenotypic frequency distribution of yield related traits
Dotted arrows indicate average phenotype of Nipponbare; Solid arrows indicate average phenotype of Zhongjiazao 17
-->

2.3 产量相关性状之间的相关性

相关性分析表明,RILs群体大部分产量性状之间存在显著的相关性（表3）。部分性状间相关性在3个环境下方向一致。其中,有效穗数与每穗粒数在3个环境中均呈现显著负相关（r≤-0.493 ）;单株产量和每穗粒数及结实率呈显著正相关（r≥0.451;r≥0.211）;千粒重与粒长及粒宽呈显著正相关（r≥0.427;r≥0.387）。以上结果表明,产量性状是复杂的性状,各性状间关系紧密,相互作用共同影响水稻产量。同时性状间相关系数受环境的影响,不同环境间会出现变化。
Table 3
表3
表3不同产量相关性状间的相关系数
Table 3Correlation coefficients among yield relate traits in RILs population

性状 Trait	地点 Station	有效穗数 EP	每穗粒数 NGPP	结实率 SSR	千粒重 TGW	粒长 GL	粒宽 GW	粒厚 GT
有效穗数EP	2015HZ
	2016HN
	2016HZ
每穗粒数NGPP	2015HZ	-0.493**
	2016HN	-0.595**
	2016HZ	-0.558**
结实率SSR	2015HZ	0.206**	0.150*
	2016HN	-0.047	-0.072
	2016HZ	-0.042	-0.464
千粒重TGW	2015HZ	-0.079	-0.143*	0.188**
	2016HN	-0.065	-0.082	-0.037
	2016HZ	-0.131	0.145*	0.288**
粒长GL	2015HZ	-0.244**	0.108	0.037	0.455**
	2016HN	-0.194**	0.050	0.039	0.427**
	2016HZ	-0.239**	0.191**	0.133	0.427**
粒宽GW	2015HZ	-0.071	0.038	-0.167*	0.387**	-0.072
	2016HN	-0.056	0.004	-0.088	0.608**	-0.146*
	2016HZ	-0.064	-0.128*	-0.007	0.430**	0.123
粒厚GT	2015HZ	-0.088	0.039	0.006	0.017	-0.008	0.101*
	2016HN	-0.039	-0.160*	-0.218**	0.634**	0.045	0.503**
	2016HZ	-0.202**	0.068	0.140	0.590**	0.154	0.438**
单株产量YPP	2015HZ	0.095	0.451**	0.378**	0.250**	0.224**	-0.070	0.025
	2016HN	-0.094	0.591**	0.738**	0.021	0.012	0.028	-0.107
	2016HZ	0.030	0.475**	0.211**	0.032	0.047	0.001	0.021

* and** mean significant levels at P<0.05 and P<0.01, respectively. The same as below*、**分别表示P<0.05、P<0.01显著水平。下同

新窗口打开

2.4 产量相关性状的QTL定位

对3个环境下产量相关性状QTL定位结果见表4和图3。在3个环境下共检测到46个QTL,分布于除第11染色体外的其他染色体不同区间,LOD值介于2.50—15.47,贡献率介于3.78%—25.45%。其中,10个QTL在3个环境下能被重复检测到,分别为qEP1、qEP2、qEP4a、qSRT1、qSRT7、qNGPE7、qTGW2、qGL3、qGL9和qGW2b。
各产量相关性状QTL具体分布为有效穗数：2015HZ、2016HN和2016HZ 3个环境下共检测到5个QTL,分布于第1、2、4和7染色体,LOD值介于2.79—8.57,贡献率介于3.78%—14.96%。所检测到的QTL中,qEP1、qEP2、qEP4a在 3个环境重复检测到,3个环境平均LOD值分别为4.20、5.42和5.19,平均贡献率分别为7.14%、9.20%和8.62%,平均加性效应值分别为-0.70、-0.77和0.84。qEP1和qEP2增效等位来自日本晴,qEP4a增效等位基因来自中嘉早17。
每穗粒数：2015HZ、2016HN和2016HZ 3个环境下共检测到5个每穗粒数QTL,分布于第1、6、7、8和10染色体,LOD值介于2.93—15.47,贡献率介于3.90%—25.45%。其中qNGPE1、qNGPE7在3个环境下重复检测到,平均LOD值分别为5.77和11.83,平均贡献率分别为10.96%和20.09%,平均加性效应值分别为13.62、19.05。qNGPE1增效等位来自中嘉早17,qNGPE7增效等位基因来自日本晴。
结实率：2015HZ、2016HN和2016HZ 3个环境下共检测到7个影响结实率的QTL,分布在第1、5、6、7、8和12染色体上,LOD值介于2.57—7.81,贡献率介于4.84%—14.90%,增效等位基因全部来自中嘉早17。其中qSRT7在3个环境下重复检测到,平均LOD值为4.71,平均贡献率为9.18%,平均加性效应值为3.01。
Table 4
表4
表4检测到的产量相关性状QTL分布情况
Table 4The distribution of yield relate traits QTLs which detected in RILs population

性状 Trait	QTL	染色体 Chr.	区间 Interval	2015HZ			2016HN			2016HZ
性状 Trait	QTL	染色体 Chr.	区间 Interval	LOD	a0	R0 (%)	LOD	a0	R0 (%)	LOD	a0	R0 (%)
有效穗数EP	qEP1	1	RM1329—RM5359	4.07	-0.62	6.67	3.53	-0.82	6.98	5.00	-0.66	7.76
	qEP2	2	RM318—RM6139	4.64	-0.64	7.08	3.04	-0.74	5.57	8.57	-0.92	14.96
	qEP4a	4	RM16720—RM17143	6.68	0.85	10.39	3.35	0.91	7.36	5.55	0.76	8.11
	qEP4b	4	RM17377—RM349	3.53	0.62	6.22				2.68	0.49	3.78
	qEP7	7	RM5875—RM6403							2.79	0.48	4.14
每穗粒数NGPP	qNGPE1	1	RM1329—RM600	4.70	12.64	8.06	8.08	14.94	16.04	4.54	13.28	8.77
	qNGPE6	6	RM5963—RM30	3.87	12.30	5.86				2.70	10.09	3.90
	qNGPE7	7	RM21242—RM11	15.06	-22.80	25.24	4.95	-11.52	9.59	15.47	-22.83	25.45
	qNGPE8	8	RM408—RM1376				3.19	8.65	5.41	2.93	9.66	4.63
	qNGPE10	10	RM7492—RM24957	2.64	9.01	3.95
结实率SSR	qSRT1a	1	RM1329—RM600	4.72	0.09	12.86	6.80	0.12	13.71
	qSRT1b	1	RM6738—RM3285							2.63	2.14	4.92
	qSRT5	5	RM163—RM3790							3.84	2.91	9.26
	qSRT6	6	RM5963—RM30							2.57	2.35	4.75
	qSRT7	7	RM3484—RM11	7.81	6.43	14.90	3.69	0.06	7.81	2.62	2.54	4.84
	qSRT8	8	RM1376—RM22694	5.17	5.33	10.14
	qSRT12	12	RM247—RM28102	2.74	-3.79	5.16	4.39	-0.07	10.28
千粒重TGW	qTGW1	1	RM128—RM8236				3.53	0.69	7.02
	qTGW2	2	RM6233—RM5812	2.66	-0.66	5.70	4.91	-0.82	9.75	4.73	-0.93	9.72
	qTGW3	3	RM14574—RM14791				2.50	0.73	6.35
	qTGW6	6	RM7434—RM5463				3.28	0.76	8.01
粒长GL	qGL3	3	RM6849—RM14791	3.76	0.16	8.85	2.54	0.12	4.39	2.68	0.16	7.46
	qGL5	5	RM6054—RM3664	2.84	-0.12	5.66
	qGL6	6	RM7158—RM3414	3.40	-0.12	6.02
	qGL7	7	RM6403—RM172	2.75	0.12	6.85	2.86	0.11	5.06
	qGL9	9	RM6839—RM205	5.80	0.17	12.74	4.05	0.14	7.96	2.53	0.12	5.06
	qGL12	12	RM5927—RM5746				2.77	0.12	6.18
粒宽GW	qGW1	1	RM8068—RM3746	2.50	0.05	6.11	6.09	0.07	11.73
	qGW2a	2	RM6367—RM12363	2.97	-0.05	5.10				2.76	-0.05	5.06
	qGW2b	2	RM6233—RM5812	4.96	-0.07	9.51	3.70	-0.06	6.01	2.80	-0.05	6.19
	qGW3	3	RM282—RM3646	6.01	-0.08	12.81	3.04	-0.05	5.57
	qGW6	6	RM7434—RM5463				2.84	0.05	5.50
	qGW9	9	RM23662—RM5526	2.59	0.04	4.33	3.51	0.05	5.71
粒厚GT	qGT1a	1	RM8068—RM3746							2.98	0.03	6.51
	qGT2b	1	RM8004—RM6738				2.86	0.03	4.61
	qGT2a	2	RM6233—RM5812				3.57	-0.03	7.38	4.12	-0.03	7.87
	qGT2b	2	RM3763—RM3542				3.03	0.03	5.25
	qGT3a	3	RM3202—RM6849	3.40	0.04	7.40	4.36	0.04	7.68
	qGT3b	3	RM3204—RM15694							3.94	-0.03	8.33
	qGT7	7	RM3484—RM214				3.03	-0.03	5.88
	qGT8	8	RM408—RM3702				2.59	-0.02	4.47	3.41	-0.03	6.23
单株产量YPP	qYPP1	1	RM5359—RM10782				4.60	2.48	10.43
	qYPP4	4	RM6659—RM16720	3.70	1.85	8.50
	qYPP7	7	RM3394—RM20897				2.50	1.69	4.82
	qYPP8a	8	RM408—RM310	4.52	2.12	11.34				2.84	0.04	6.70
	qYPP8b	8	RM23175—RM3845				5.15	2.68	11.35

新窗口打开

显示原图|下载原图ZIP|生成PPT
图3分子标记连锁图谱及重复检测到的QTL和已克隆产量相关基因在染色体上的分布
E1：2015杭州;E2：2016海南;E3：2016杭州;“+”和“–”分表代表增效等位基因来自中嘉早17和日本晴
-->Fig. 3The genetic linkage map and the location of the repeatability QTLs or cloned gene for yield related traits on marker interval
E1: 2015HZ; E2: 2015HN; E3: 2016HZ; “+” and “–” represent the addictive effects of the alleles are contributed by Zhongjiazao 17 and Nipponbare, respectively
-->

千粒重：2015HZ、2016HN和2016HZ 3个环境下共检测到4个千粒重QTL,分别位于第1、2、3和6染色体上,LOD值介于2.50—4.91,贡献率介于6.35%—7.95%。其中只有qTGW2在3个环境下重复检测到,平均LOD值为4.10,平均贡献率为8.39%,平均加性效应值为-0.80,增效等位基因来自日本晴,而其他千粒重QTL增效等位基因多数来自中嘉早17。
粒长：2015HZ、2016HN和2016HZ 3个环境下共检测到6个粒长QTL,分布在第3、5、6、7、9和12染色体上,LOD值介于2.53—5.80,贡献率介于4.39%—12.74%。其中3个环境下重复检测到qGL3和qGL9,平均LOD值分别为2.99和4.29,平均贡献率分别为6.90%和8.59%,平均加性效应值分别为0.99和0.14,这两个QTL增效等位基因均来自中嘉早17。
粒宽：2015HZ、2016HN和2016HZ 3个环境下共检测到6个粒宽QTL,分布在第1、2、3、6和9染色体上,LOD值介于2.59—6.09,贡献率介于4.33%—12.81%。其中qGW2b在3个环境下重复检测到,平均LOD值分别为3.82,平均贡献率分别为7.24%,平均加性效应值分别为-0.06,增效等位基因均来自日本晴。
粒厚：2015HZ、2016HN和2016HZ 3个环境下共检测到8个粒厚QTL,分布于第1、2、3、7和8染色体上,LOD值介于2.59—4.36,贡献率介于4.47%—8.33%。无3个环境下重复检测到的QTL,但是多数QTL增效等位基因来自中嘉早17。
单株产量：2015HZ、2016HN和2016HZ 3个环境下共检测到5个单株产量QTL,分别位于第1、4、7和8染色体上,LOD值介于2.50—5.15,贡献率介于4.80%—11.35%。所检测到的QTL增效等位基因均来自中嘉早17。

2.5 产量性状QTL环境互作

应用QTLNetwork2.2的MCIM分析方法联合检测产量相关性状与环境互作效应（表5）,共检测到6个QTL与环境互作效应显著,分别为控制每穗粒数的qNGPE1和qNGPE7;控制结实率的qSRT1a和qSRT7;控制单株产量的qYPP1和qYPP7。与环境互作效应QTL主要集中在第1染色体RM1329—RM600和第7染色体RM214—RM11区间。表明每穗粒数、结实率的和单株产量性状易受环境影响。
Table 5
表5
表5产量性状QTL与环境的互作效应
Table 5QTLs with addictive/environment (AE) interaction for yield related traits

性状 Trait	QTL	染色体 Chr.	区间 Interval	环境互作效应 Environment interaction effect
性状 Trait	QTL	染色体 Chr.	区间 Interval	2015HZ	2016HN	2016HZ
每穗粒数NGPP	qNGPE1	1	RM1329—RM600	5.6051**	4.2916*	-9.7202**
	qNGPE7	7	RM21242—RM11	-9.3251**	0.5841	8.7623**
结实率SSR	qSRT1a	1	RM1329—RM600	-7.1407**	-3.7795*	10.7944**
	qSRT7	7	RM3484—RM11	11.3593**	5.2059**	-16.0430**
单株产量YPP	qYPP1	1	RM5359—RM10782	-0.9237	-0.1243	1.0495*
	qYPP7	7	RM3394—RM20897	0.8749	0.3000	-1.1921*

新窗口打开

3 讨论

3.1 籼粳交分子连锁图谱

籼稻与粳稻是具有一定生殖隔离的2个主要栽培稻亚种,它们之间在形态性状、农艺性状和生理生化反应等方面存在明显的差异。本研究构建RIL群体选用的亲本为粳稻品种日本晴和籼稻品种中嘉早17。其中中嘉早17连续多年被中国农业部推荐为南方主栽水稻品种,具有高产、稳产等优良特征;而另一亲本日本晴具有最为完整的水稻基因组参考序列。同时中嘉早17的每穗粒数、结实率、单株产量、千粒重和粒长等产量相关显著明显优于日本晴（表2）。另外籼粳亚种间杂交后代基因型具有较大的变幅,因此,利用日本晴和籼稻品种中嘉早17衍生的RIL群体非常适合构建分子连锁遗传图谱开展水稻重要农艺性状QTL定位研究。而本研究结果也显示利用该重组自交系构建的连锁图谱具有丰富的多态性标记（共有163个多态性标记）,覆盖水稻基因组的93.64%,可较好地满足产量性状QTL定位需求。
在遗传作图研究中,标记基因型偏分离现象普遍存在。董少玲等^[29]以TD70和Kasalath构建的重组自交系,51.8%的标记显著偏分离,偏向父本TD70;沈希宏等^[30]用协优9308和协青早B构建的重组自交系48.5%的标记显著偏分离,大部分偏向父本协优9308;徐小飒等^[31]以培矮64和93-11构建的重组自交系,50.4%标记显著偏分离,大部分偏向父本培矮64。本研究构建的连锁图谱中有44个标记也出现了显著偏分离,占标记总标记的27.0%,其中大部分偏向父本中嘉早17。目前对偏分离的遗传机制尚不清楚,推测这可能亲本间的亲缘关系、雌雄配子体的选择性、群体构建过程中人工抽样有关。在实际运用中少量的偏分离标记对QTL检测的影响是有限的,会随着QTL间连锁距离的增大而消失^[32]。

3.2 QTL定位比较

在不同生态环境下检测农艺性状QTL,在理论上可以减少环境因素的影响,有利于发现在多种环境下稳定表达的主效QTL。本研究在两年三季3个环境下开展了水稻产量相关性状的QTL分析,共检测到了46个QTL。分析发现,qEP2位置存在已克隆的控制有效穗数基因MADS57^[33];而qEP4a在3个环境下均能被检测到,平均贡献率为7.62%,目前尚未发现有关报道。qEP4b与占小登等^[34]发现的有效穗数QTL位于相同的区域。qNGPE1与CHEN^[35]精细定位的抽穗期和穗粒数QTL qHD1在相近的区间。qNGPE6、qNGPE7和qNGPE8区间分别包含已经克隆同时控制抽穗期、株高、穗粒数的多功能基因位点Hd1^[36]、Ghd7^[11]和DTH8^[12]。qGL3和qGL9区间分别包含已经克隆的控制粒长的基因PGL1^[37]和SG1^[38];qGW2b和qGW3区间分别覆盖已被克隆、报道控制粒宽的功能基因GW2^[13]和GS3^[16]。上述结果表明,本研究检测到的主效QTL大部分与前人定位或克隆的产量相关基因位置一致,同时本研究也发现了一些新的主效QTL,如控制有效穗数的qEP4a等（图3）。另外,本研究发现的多个产量相关主效QTL其有利于增加产量的等位基因均来自中嘉早17,如控制有效穗数的qEP4a,控制每穗粒数的qNGPE1,控制粒长的qGL3和qGL9,控制粒长的qGL3。单株产量与结实率相关QTL的增效等位基因也大部分来自中嘉早17。该研究结果与中嘉早17的每穗粒数、结实率、单株产量、千粒重和粒长等产量相关明显优于日本晴的结果一致,这些产量增效QTL可能是保证中嘉早17高产稳产的遗传基础。

3.3 产量QTL与环境互作

水稻产量相关性状多属于遗传机制较为复杂的数量性状,其表现型不仅受主效基因和微效基因多基因调控、同时易受环境条件的影响。在不同的群体或者不同环境条件下控制同一性状的QTL表达情况有一定差异,即使同一群体在不同环境下检测到的QTL也不尽一致,因此,环境是影响产量相关性QTL表达的重要因素之一^[39]。占小登等^[34]在2个环境下共检测到产量相关性状30个,其中10个能在不同环境下重复检测到,5个与环境互作显著;苏相文等^[3]在2个环境下共检测到产量相关性状31个,8个能重复检测到,21个对环境敏感;杨占烈等^[40]在3个环境条件下检测到26个QTL,4个能重复检测到;本研究利用RIL群体进行产量相关性状QTL定位,3个同环境下检测到46个QTL,其中10个QTL在3个环境下都能被重复检测到。这些能在不同环境中都能检测到的QTL可能与水稻品种的广适性有关。本研究也发现了6个产量性状QTL与环境互作效应显著（表5）,这6个QTL主要是调控水稻结实率、每穗粒数及单株产量的QTL,说明这些性状对环境较为敏感。

4 结论

以日本晴×中嘉早17构建的重组自交系群体连锁图谱具有丰富的多态性标记,覆盖水稻基因组的93.64%,可较好地满足水稻重要农艺性状QTL定位要求。利用该套群体检测到多个产量相关性状QTL,其中,多数控制每穗粒数、结实率、粒长和单株产量的QTL的增效等位基因均来自中嘉早17。该结果与中嘉早17的每穗粒数、结实率、单株产量、千粒重和粒长等产量相关显著明显优于日本晴的结果一致,这些产量增效QTL可能是中嘉早17高产、稳产的遗传基础。
The authors have declared that no competing interests exist.