陈雪萍^*, 荆凌云^*, 王嘉, 荐红举, 梅家琴, 徐新福, 李加纳, 刘列钊^*
西南大学农学与生物科技学院 / 重庆市油菜工程技术研究中心, 重庆 400715
^*通讯作者(Corresponding author): 刘列钊, E-mail:liezhao2003@126.com, Tel: 023-68251383 第一作者联系方式: E-mail:1473718700@qq.com
^**同等贡献(Contributed equally to this work)
收稿日期:2017-01-07 接受日期:2017-04-20网络出版日期:2017-05-08基金:本研究由国家自然科学基金项目(31371655)和重庆市科委(cstc2016shmszx80083)项目资助

摘要菌核病是一类非专一性的植物真菌病原菌, 寄主范围广泛, 严重危害农作物的生产。对高世代重组自交系群体(RIL)及F₂群体终花期茎秆进行菌核病抗性接种鉴定, 根据构建的近红外模型对接种鉴定的茎秆木质素含量、单体组分比例进行测定, 并进行相关性分析和QTL定位。结果表明在2013年和2014年RIL群体茎秆菌斑大小与木质素含量呈极显著负相关, 相关系数分别为-0.348和-0.286, 与单体G/S呈显著正相关, 相关系数分别为0.198和0.167。2014年F₂群体菌斑大小与木质素含量呈极显著负相关, 相关系数为-0.306, 与单体G/S相关性为0.142。F_2:3家系抗(感)植株茎部切片间苯三酚染色观察表明抗性较强的材料木质素含量高于抗性较弱的材料。根据已构建的重组自交系高密度SNP遗传图谱, 利用复合区间作图法对上述性状进行QTL分析, 共检测到18个QTL, 其中9个菌核病抗性相关QTL分布于A05、A06、C04和C06染色体, 单个QTL可解释的表型变异为2.38%~12.05%; 3个木质素含量QTL分别位于A04、A05和C01染色体, 单个QTL可解释表型变异的2.03%~13.75%。6个木质素单体G/S QTL分布于A08、C03和C07染色体, 单个QTL可解释表型变异的2.06%~8.66%。本文研究结果为油菜菌核病抗性育种提供了新的思路和理论基础。

关键词:甘蓝型油菜; 菌核病; 相关性分析; 数量性状位点; 木质素; 单体G/S
Correlation Analysis of SclerotiniaResistance with Lignin Content and Monomer G/S and Its QTL Mapping in Brassica napus L.
CHEN Xue-Ping^*, JING Ling-Yun^*, WANG Jia, JIAN Hong-Ju, MEI Jia-Qin, XU Xin-Fu, LI Jia-Na, LIU Lie-Zhao^*
College of Agronomy and Biotechnology, Southwest University / Chongqing Engineering Research Center for Rapeseed, Chongqing 400715, China
Fund:This study was supported by the National Natural Science Foundation of China (31371655) and the Science and Technology Committee of Chongqing (cstc2016shmszx80083)
Abstract Sclerotinia sclerotiorumis a fungal pathogen causing disease in a wide range of plants, resulting in serious damage in crop production. The detached stem inoculation assay of RIL and F₂ populations at final flowering stage was conducted, near infrared (NIR) spectroscopy was used to measure lignin content and monomer G/S in the stem, and correlation analysis and QTL mapping for these traits were performed. The lesion size of the RIL had a significantly negative correlation with lignin content, with a correlation coefficient at -0.348 and -0.286 in 2013 and 2014, respectively. The monomer G/S was significantly correlated with lesion size in the RIL population, and the correlation coefficient at 0.198 and 0.167 in 2013 and 2014, respectively. The lesion size of F₂ in 2014 was significantly and negatively correlated with lignin content in the stem, with a correlation coefficient at -0.306. The cross sections of resistant and susceptible plants from F_2:3family were stained with phloroglucinol-HCl, showing that the content of lignin was significantly lower in the less resistant materials than in more resistant plants. According to the high density SNP genetic maps and composite interval mapping, a total of 18 QTLs were identified, which were located on A04, A05, A06, A08, C01, C03, C04, C06, and C07 chromosomes, with the explained phenotypic variation by individual QTL ranging from 2.38% to 12.05% for nine QTLs of lesion size, from 2.03% to 13.75% for three QTLs of lignin content, and from 2.06% to 8.66% for six QTLs of monomer G/S. The research results provide some new insights for the Sclerotiniaresistance breeding in B. napus.

Keyword: Brassica napus; Sclerotinia sclerotiorum; Correlation; QTL; Lignin; Monomer G/S
Show Figures
Show Figures








油菜(Brassica napusL.)属于十字花科芸薹属, 是世界范围内广泛种植的油料作物之一, 为全世界提供15%食用植物油^[1]。油菜在中国占据重要的地位, 每年国内一半的食用植物油均为菜籽油^[2]。菌核病作为油菜三大病害(菌核病、根肿病和病毒病)之一, 严重影响我国油菜的产量和品质, 每年造成的损失高达5%~30%。长江流域的发病情况最为严重, 轻病年可使油菜减产10%~30%, 重病年甚至可达80% ^[3]。该病害在澳大利亚可导致油菜年产量减少24% ^[4], 在北美以及欧洲主要油菜种植区危害也越来越严重^[5]。核盘菌的寄主范围非常广泛, 产生的菌核以及子囊孢子能在土壤中存活若干年, 给该病的预防和治疗带来很大困难。尽管化学农药能在一定程度上抑制该病菌, 但植株发育阶段、外界环境以及施用时间等条件的不同都会影响预防效果^[6], 令生产成本增加^[7], 环境负担加重。选择和培育抗病品种是控制油菜菌核病最经济有效的途径, 但到目前为止, 仅发现一些对菌核病具一定生理抗性或避逆性的品种, 还未培育出对菌核病免疫的品种^[8]。
菌核病作为油菜生产的第一大病害, 倍受广大遗传育种家的关注。随着QTL定位、关联分析及转录组测序分析在油菜菌核病研究中的普及和应用, 对油菜菌核病的研究取得了很大进展。Zhao等^[9]在F_2:3家系中检测到3个与叶片菌核病抗性相关的QTL, 分布在LG3、LG12和LG17连锁群上, 3个与茎秆菌核病抗性相关的QTL, 分布于LG7、LG10和LG15连锁群; Wu等^[10]对DH群体苗期叶片和成熟期茎秆的QTL定位, 共检测到13个核盘菌抗性相关的QTL, 重复检测到2个主效QTL, 分别位于A09和C06染色体, 解释表型变异率为8.54%~32.61%; 周李鹏等^[11]利用DH群体共检测4个菌核病抗性相关的QTL, 分布于A02、A03和C02连锁群, 解释8.00%~13.80%的表型变异; 梅家琴^[12]在甘蓝F₂无性系群体中检测到8个与叶片菌核病抗性相关的QTL, 4个与茎秆菌核病抗性相关的QTL, 位于C01、C03、C06、C08和C09染色体, 解释表型变异率为4.9%~18.8%; Wei等^[13]对甘蓝型油菜全基因组关联分析, 检测到17个茎秆菌核病抗性相关的QTL, 分别位于A08和C06连锁群; Wu等^[14]对接有核盘菌的油菜抗感植株进行转录组分析, 筛选到9001个差异表达基因, 它们主要参与病原体识别、MAPK信号级联、WRKY转录调控、茉莉酸和乙烯信号转导以及防御蛋白生物合成等途径。
木质素作为植物机械组织中重要组成部分和结构抗性物质, 其代谢途径与其他次生物质代谢途径相互交叉, 在植物防御病原体的过程中发挥重要的作用^{[15, 16, 17]}。当植物受到病原体入侵或是机械损伤时便会大量产生木质素及其类似物^{[18, 19, 20]}, 而木质化的细胞壁可以增强草本和木本植物的抗菌能力^{[21, 22, 23]}。杨向东^[24]发现油菜木质素含量与抗菌核病之间存在显著正相关性。虽然疾病诱导木质素合成相关基因被激活的分子机理尚不清楚, 但在病原体入侵等生物胁迫下会改变木质素的组成即单体比例^{[25, 26, 27, 28]}。木质素单体之间的相对比例及不同类型之间的化学键决定着木质素的复杂性和特异性。油菜是双子叶植物, 木质素主要为S、G型, H型木质素平均仅占总木质素含量的10.96%。S型木质素(syringyl lignin)由紫丁香基丙烷结构单体聚合而成, G型木质素(guaiacyl lignin)由愈创木基丙烷结构单体聚合而成, 不同材料间木质素单体S/G有极显著差异。因此本研究主要从单体组分比例(G/S)入手, 探讨其与油菜菌核病抗性之间的关系。
目前对油菜木质素含量研究的主要目的是降低种子其含量以提高饼粕利用效率。Snowdon等^[29]利用SSR标记在甘蓝型油菜A09染色体上检测到1个与种皮木质素含量相关的主效QTL; Liu等^[30]利用BSA法对RIL和BC₁F₂群体进行SSR标记分析, 检测到1个与种子木质素含量相关的主效QTL, 位于A09连锁群。曲存民等^[31]在2个环境下共检测到12个种皮木质素含量QTL, 分布于4个连锁群, 可解释表型变异的4.50%~8.79%; Liu等^[32]利用高密度SNP遗传图谱, 对4个环境下RIL群体种子木质素含量进行QTL定位, 共检测到8个QTL, 在所有环境中均能稳定检测到2个主效QTL, 分别位于A09和C05染色体, 解释的表型变异范围分别为31.6%~ 42.8%和8.1%~14.1%。木质素作为苯基丙烷衍生物单体构成的聚合物, 极大地限制了生物燃料的开发和应用。Chen等^[33]在转基因苜蓿中通过下调木质素合成过程中相关基因的表达, 能大幅提高植物中糖化作用的效率以及纤维素乙醇的产量。黄杰恒^[34]分析表明, 干旱处理和浇水对照条件下, 木质素单体S/G与抗折性均极显著正相关。Eynck等^[35]以2个菌核病抗性不同的亚麻芥为材料, 发现在抗病植株中木质素单体S/G高于感病植株。在油菜中木质素单体G/S与菌核病抗性间的相关性以及木质素单体G/S相关的QTL尚未见报道。
1 材料与方法1.1 试验材料和田间试验高世代重组自交系群体亲本分别为GH06和P174, F₂群体以中双11 (菌核病抗性较强)和E54 (菌核病抗性较弱)为亲本, 选用甘蓝型油菜中双11为抗病对照。随机选取172份重组自交系进行SNP标记分析, 构建高密度SNP遗传连锁图谱。
2013年10月播种RIL10群体, 次年4月鉴定其离体茎秆抗性; 2014年10月播种RIL11、F₂群体, 次年4月鉴定其离体茎秆抗性。选取中双11为对照, 种植所有材料于重庆市油菜工程中心试验田, 随机区组设计3个重复, 每个小区2行, 每行10株。行距0.40 m, 株距0.24 m, 按常规生产方式进行田间管理, 确保所有样本的外部生长环境一致。
1.2 菌核病抗性鉴定将纯化的菌核于PDA培养基传代培养至第3代。终花期对172份RIL群体和311份F₂群体茎秆离体接种, 随机选取每份材料3株, 距离地面20 cm处截取长度约30 cm的主茎, 将茎秆两端用保鲜膜缠绕, 以防水分散失和伤到茎秆表面而影响试验结果。将茎秆排布在铺有湿毛巾和滤纸的铝合金框架上, 将PDA菌丝块接种于茎秆平滑处, 带菌面紧贴茎秆表面, 注意避开分枝结口部位, 每段茎秆接种2个菌丝块(两菌丝块间隔10 cm)。接种完毕, 将温湿度计与接种的茎秆一同密闭, 设置室内温度为22℃, 相对湿度为90%以上。接种后第3天统计菌斑长度。通过“ 相对感病度” (R)衡量材料抗性水平, R=V/V_对照, 即待测植株菌斑大小与对照中双11病斑大小的相对比值。R值越小表明材料抗性越强, 反之代表材料越易感病, 若R=1, 则表明材料的抗性与中双11相当。
1.3 菌斑大小、木质素含量和单体G/S测定及相关性分析将接种鉴定后的离体茎秆晒干, 于60℃烘箱中彻底烘干(约48 h), 用FW80-1粉碎机彻底粉碎, 收集至自封袋, 密闭于干燥器中保存。用构建的近红外模型测定2013年和2014年种植的RIL群体及2014年种植的F₂群体样品粉末木质素含量和单体G/S, 用SPSS18.0软件和Microsoft Excel 2013软件统计数据及分析表型相关性。
1.4 茎秆切片观察终花期从F_2:3家系中选取株高一致的菌核病抗(感)性材料各2株, 于距离植株顶端10 cm处徒手横切茎端组织进行间苯三酚染色。将Leica CM1800冷冻切片机机箱温度降至-20℃。截取油菜茎秆(长约0.1~0.3 cm), 置样品托头上, 用OCT包埋剂包埋, 待包埋剂变白后将其移至切片机的机头上, 设置切片厚度(25~100 μ m), 进行切片。挑选组织结构完整的薄片, 用毛笔将其放入盛水的培养皿中清洗包埋剂。将组织薄片转移到载玻片上, 吸干剩余水分, 滴加1%的间苯三酚染色, 静置5 min, 再加入1滴35%浓盐酸显色。将载玻片置Olympus DP80型光学显微镜下观察, 并采用cellSens Standard软件采集图片, NIS-Elements AR 3.2分析图片。
1.5 QTL定位分析所用的重组自交系高密度SNP遗传图谱, 包含2795个SNP多态性标记位点, 总长1832.9 cM, 相邻标记间平均距离为0.66 cM ^[32]。采用软件Windows QTL Cartographer 2.5及复合区间作图(Composite Interval Mapping, CIM)法对2年RIL群体的菌斑大小、木质素含量以及单体G/S进行QTL定位及效应检测。CIM分析时, 选取1 cM的步长(walking speed), 按照假定检测10和Zmapqtl模型6, 选取参数1000次回归, 显著水平为0.01, 检测区间可能存在的QTL, 软件分析的同时给出性状QTL的加性效应和解释的表型变异。当加性效应大于0时, 表明性状效应值来自父本P174, 加性效应小于0时, 性状效应值来自母本GH06。参照McCouch等^[36]的方法, 以“ q” 加年份再加染色体编号和性状命名检测到的QTL。

2 结果与分析2.1 甘蓝型油菜茎秆相对感病度、木质素含量和单体G/S的测定如图1所示, 两年里RIL群体和F₂群体中相对感病度、木质素含量及单体G/S均呈明显的正态分布, 变幅范围表现出双向超亲分离, 表明菌核病抗性、木质素含量和单体G/S是受多基因控制的典型数量性状。偏斜度和峰度值几乎都小于1, 适合于QTL分析(表1和表2)。
图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 RIL和F2群体茎秆相对感病度、木质素含量和单体G/S在2013年和2014年的频率分布Fig. 1 Frequency distribution of stem relative susceptibility, lignin content and monomer G/S fromB. napusRIL in 2013 and 2014 and F2in 2014

表1
Table 1
表1(Table 1)

表1 RIL群体茎秆相对感病度、木质素含量及单体G/S的表型分布特征 Table 1 Phenotype distribution feature of stem relative susceptibility, lignin content and monomer G/S for RIL populations 年份和性状 Year and trait 亲本Parent RIL群体RIL population GH06 P174 范围 Range 均值 Mean 标准差 SD 变异系数 CV (%) 方差 Variance 峰度 Kurtosis 偏度 Skewness 2013R 0.94 0.77 0.28-1.36 0.83 0.21 25.90 0.05 -0.41 0.15 2014R 1.11 0.85 0.34-1.33 0.85 0.20 23.60 0.04 -0.24 0.32 2013ADL 7.60 9.29 4.20-11.80 8.94 1.61 18.00 2.60 0.66 -0.88 2014ADL 8.47 9.60 6.02-11.63 9.55 0.94 10.00 0.89 0.41 -0.25 2013G/S 0.93 0.81 0.75-1.40 1.10 0.11 10.00 0.01 0.27 0.25 2014G/S 0.89 0.83 0.30-1.30 0.80 0.18 23.00 0.03 0.34 0.06 R: relative susceptibility; ADL: lignin content; G/S: monomer G/S. R: 相对感病度; ADL: 木质素含量; G/S: 单体G/S。

表1 RIL群体茎秆相对感病度、木质素含量及单体G/S的表型分布特征 Table 1 Phenotype distribution feature of stem relative susceptibility, lignin content and monomer G/S for RIL populations

表2
Table 2
表2(Table 2)

表2 F2群体茎秆相对感病度、木质素含量及单体G/S的表型分布特征 Table 2 Phenotype distribution feature of stem relative susceptibility, lignin content and monomer G/S for F2 populations 年份和性状 Year and trait 亲本Parent F2群体F2 population 中双11 E54 范围 Range 均值 Mean 标准差 SD 变异系数 CV (%) 方差 Variance 峰度 Kurtosis 偏度 Skewness 2014R 1.00 1.25 0.25-1.90 0.96 0.32 32.90 0.10 0.23 0.63 2014ADL 10.34 9.06 4.23-13.44 10.48 1.21 12.00 1.45 2.23 -0.72 2014G/S 0.84 0.90 0.45-1.25 0.93 0.15 16.00 0.02 -0.29 -0.31 R: relative susceptibility; ADL: lignin content; G/S: monomer G/S. R: 相对感病度; ADL: 木质素含量; G/S: 单体G/S。

表2 F2群体茎秆相对感病度、木质素含量及单体G/S的表型分布特征 Table 2 Phenotype distribution feature of stem relative susceptibility, lignin content and monomer G/S for F2 populations

2.2 油菜终花期菌斑大小与木质素含量、单体G/S的相关性分析如表3所示, 2013年RIL群体菌斑大小与木质素含量呈极显著负相关, 相关系数为-0.348, 与单体G/S呈显著正相关, 相关系数为0.198, 木质素含量与单体G/S之间无显著相关性; 2014年RIL群体菌斑大小与木质素含量呈极显著负相关, 相关系数为-0.286, 与单体G/S呈显著正相关, 相关系数为0.167, 木质素与单体G/S之间相关性未达到显著水平; 表4所示, 2014年F₂群体菌斑大小与木质素含量呈极显著负相关, 相关系数为-0.306, 与单体G/S相关系数为0.142, 木质素含量与单体G/S之间无显著相关性。这说明木质素和单体G/S在油菜菌核病防御过程中均发挥重要的作用。同一性状在2013和2014年RIL群体及2014年F₂群体中均表现为极显著正相关, 相关系数较高, 说明菌斑大小、木质素含量、单体G/S性状遗传稳定。
表3
Table 3
表3(Table 3)

表3 2013年和2014年RIL群体菌斑大小、木质素含量及单体G/S的相关性分析 Table 3 Correlation coefficients among lesion size , lignin content and monomer G/S for RIL population in 2013 and 2014 性状 Trait 13LS 14LS 13ADL 14ADL 13G/S 2014LS 0.865* * 2013ADL -0.348* * -0.280* * 2014ADL -0.236* * -0.286* * 0.219* * 2013G/S 0.198* 0.178* 0.064 -0.152 2014G/S 0 0.167* -0.073 0.088 0.206* * LS: lesion size; ADL: lignin content; G/S: monomer G/S. * and * * denote significant correlation at the 0.05 and 0.01 probability levels, respectively. LS: 菌斑大小; ADL: 木质素含量; G/S: 单体G/S。* 和* * 分别代表在0.05和0.01显著水平。

表3 2013年和2014年RIL群体菌斑大小、木质素含量及单体G/S的相关性分析 Table 3 Correlation coefficients among lesion size , lignin content and monomer G/S for RIL population in 2013 and 2014

表4
Table 4
表4(Table 4)

表4 2014年F2群体菌斑大小、木质素含量及单体G/S的相关性分析 Table 4 Correlation coefficients among lesion size, lignin content and monomer G/S for F2 population in 2014 性状 Trait 14LS 14ADL 2014ADL -0.306* * 2014G/S 0.142 0.079 LS: lesion size; ADL: lignin content; G/S: monomer G/S. * and * * denote significant correlation at the 0.05 and 0.01 probability levels, respectively. LS: 菌斑大小; ADL: 木质素含量; G/S: 单体G/S。* 和* * 分别代表在0.05和0.01显著水平。

表4 2014年F2群体菌斑大小、木质素含量及单体G/S的相关性分析 Table 4 Correlation coefficients among lesion size, lignin content and monomer G/S for F2 population in 2014

2.3 茎秆切片观察细胞产生的木质素会在细胞壁沉积, 间苯三酚-盐酸染液能将木质化的细胞壁染成樱桃红色或红紫色, 颜色越深, 木质素含量越高。A和C分别代表2个菌核病抗性相对较弱的材料(Z258和Z234), B和D分别代表2个菌核病抗性相对较强的材料(Z219和Z232)。图2为抗(感)植株茎秆局部冷冻横切图片, 在显微镜视野下观察到木质素含量以抗病较弱植株(A、C)低于抗病较强植株(B、D)。
图2
Fig. 2

	Figure OptionViewDownloadNew Window
	图2 抗病较弱材料(A, C)和抗病较强材料(B, D)在终花期的茎部横切面间苯三酚-盐酸染色观察以示木质素含量差异(10× 4)Fig. 2 Stem cross-section stained with phloroglucinol-HCl to show the difference of lignin content in the less resistant (A, C) and more resistant (B, D) lines by the end of flowering time (10× 4)

2.4 茎秆菌斑大小、木质素含量和单体G/S的QTL分析总共检测到18个QTL (表5和图3), 菌核病抗性相关的QTL有9个, 分布于A05、A06、C04和C06染色体, 单个QTL可解释的表型变异为2.38%~12.05%, LOD为2.78~7.45, 性状效应值均来自母本。两年内在A05染色体上均检测到1个有重叠置信区间的菌核病抗性QTL (q2013A05LS-2和q2014A05LS-2), 加性效应均为负, 说明其性状效应值来自母本GH06, 即抗性增效基因来自父本P174。对于木质素含量检测到3个QTL, 分别位于A04、A05和C01染色体, LOD为3.00~3.61, 解释表型变异为2.03%~7.35%, 其中q2013A05Lig与q2013A05LS-1的置信区间部分重叠。对于木质素单体G/S共检测到6个QTL, 分布于A08、C03和C07染色体, 单个QTL可解释表型变异在2.06%~8.66%之间, LOD值介于2.67~3.40之间, 在C07染色体上也检测到1个有重叠置信区间的木质素单体G/S的QTL (q2013C07LigG/S-1、q2013C07LigG/S-2和q2014C07LigG/S), 其加性效应为正, 说明增效基因均来自父本P174。检测到的数量性状位点中, 与菌斑大小相关的q2013A05LS-2和与木质素含量相关的q2014C01Lig, 其表型贡献率最大(R²> 12%)。
表5
Table 5
表5(Table 5)

表5 利用复合区间作图法检测到茎秆菌斑大小(LS)、木质素含量(Lig)及单体G/S (Lig G/S)在2个环境中的QTL Table 5 Putative QTL detected by CIM for stem lesion size (LS), lignin content (Lig), and monomer G/S (Lig G/S) in two environments 数量性状位点 QTL 染色体 Chromosome 标记区间 Marker interval (cM) 阈值 LOD 贡献率 R2 (％) 加性效应 Additive 置信区间 Confidence interval q2013A05LS-1 A05 SNP11632A05-SNP11708A05 3.40 9.18 -0.47 29.56-35.13 q2013A05LS-2 A05 SNP11775A05-SNP4033A02 4.84 12.05 -0.54 41.20-47.54 q2014A05LS-1 A05 SNP11166A05-SNP11242A05 3.55 5.17 -0.22 3.96-6.85 q2014A05LS-2 A05 SNP11775A05-SNP3643A02 3.39 7.40 -0.37 41.20-43.70 q2014A05LS-3 A05 SNP12064A05-SNP11984A05 5.38 4.39 -0.19 62.87-68.61 q2014A05LS-4 A05 SNP11825A05-SNP11305A05 2.78 2.38 -0.19 73.90-76.90 q2014A06LS A06 SNP14966A06-SNP11305A05 3.26 5.28 -0.26 57.00-76.90 q2014C04LS C04 SNP45588-SNP42480 7.45 5.74 -0.26 2.56-6.60 q2014C06LS C06 SNP29157C06-SNP29132 3.40 8.90 -0.65 99.75-102.11 q2013A05Lig A05 SNP11594A05-SNP11716A05 3.61 7.35 0.48 22.27-34.10 q2014A04Lig A04 SNP8405A04-SNP9618A04 3.00 2.03 0.14 61.18-83.63 q2014C01Lig C01 SNP27002C01-SNP34230C01 3.54 13.75 -0.83 0-14.33 q2013A08LigG/S A08 SNP12620A05-SNP12606A05 2.85 5.67 -0.03 11.56-13.52 q2013C03LigG/S-1 C03 SNP34516C03-SNP37474 2.67 5.04 0.02 15.89-28.43 q2013C03LigG/S-2 C03 SNP36062C03-SNP31881 3.20 4.99 0.02 34.10-41.00 q2013C07LigG/S-1 C07 SNP31994C07-SNP28392C07 3.10 8.66 0.04 18.60-34.33 q2013C07LigG/S-2 C07 SNP31994C07-SNP26531C07 3.40 7.41 0.03 18.60-43.00 q2014C07LigG/S C07 SNP26531C07-SNP25787C07 2.79 2.06 0.04 43.00-55.14

表5 利用复合区间作图法检测到茎秆菌斑大小(LS)、木质素含量(Lig)及单体G/S (Lig G/S)在2个环境中的QTL Table 5 Putative QTL detected by CIM for stem lesion size (LS), lignin content (Lig), and monomer G/S (Lig G/S) in two environments

图3
Fig. 3

	Figure OptionViewDownloadNew Window
	图3 甘蓝型油菜茎秆菌斑大小、木质素含量及单体G/S QTL在SNP连锁群上分布图 QTL的染色体上仅展示了QTL区段的标记及染色体两端各2个标记。Fig. 3 Putative QTL locations of lesion size, lignin content, and monomer G/S for stem on the SNP genetic map Showing the markers in the QTL confidence intervals, along with the terminal two markers at each end of chromosomes containing the QTL.

3 讨论近年来, 对甘蓝型油菜茎秆木质素含量的研究逐渐增多, 尤其在茎秆抗病抗倒伏方面。杨向东等^[24]分析了29个油菜材料茎秆木质素含量与抗菌核病间的相关关系, 相关系数为0.6971, 表明木质素含量与油菜抗菌核病之间存在显著正相关性。本研究表明, 2013年和2014年RIL群体菌斑大小和木质素含量相关系数分别为-0.348和-0.286, 2014年F₂群体菌斑大小和木质素含量相关系数为-0.306, 菌核病抗性与木质素含量均呈极显著正相关, 与杨向东的结果一致。对F_2:3家系菌核病抗(感)材料冷冻切片观察也表明, 抗性强的材料中木质素含量明显高于抗性弱的材料, 与上述结论一致。木质素单体之间的相对比例及不同类型之间的化学键决定着木质素的复杂性、特异性及稳定性。在病原体入侵等生物胁迫下木质素的组成即单体比例会发生改变^{[25, 26, 27, 28]}。黄杰恒^[34]分析表明, 干旱处理和浇水对照条件下, 油菜木质素单体S/G与抗折性均极显著正相关。Eynck等^[35]发现, 亚麻芥木质素单体S/G在抗性植株中(1.7)高于在感性植株中(1.2)。在我们的研究中发现2013年和2014年RIL群体菌斑大小和单体G/S均呈显著正相关, 相关系数分别为0.198和0.167, 2014年F₂群体中与单体G/S相关系数为0.142, 菌核病抗性与木质素单体G/S均呈负相关, 这与Eynck等^[35]的结果一致。上述研究结果表明木质素单体组分比例(G/S)不仅与油菜抗倒伏相关, 还影响油菜菌核病抗性。
Li等^[37]对油菜茎秆菌核病抗性相关QTL整合分析, 最终得到的保守QTL位于A09 (22.5~27.5 Mb)和C06 (29.5~36.1 Mb)染色体。尽管油菜菌核病抗性QTL定位的研究结果具有一定共性, 但是目前检测到的QTL分布相对比较分散, 很少能重复检测到QTL。目前油菜菌核病抗性相关QTL主要集中在A01、A02、A03、A05、A06、A08、A09、C02、C04、C06、C07、C08和C09染色体上^{[3, 9, 10, 11, 12, 13, 38, 39, 40]}, 鉴定出的油菜抗菌核病QTL数量较多但主效QTL较少, 且单个QTL能解释的表型变异有限。QTL检测的方法、被检测群体的种类和大小、接种的环境、方法和接种年份都会影响QTL的检测, 因此能检测出的主效QTL较少。本研究连续两年对RIL群体菌核病抗性进行QTL定位, 检测到9个与菌核病抗性相关的QTL, 分布于A05、A06、C04和C06染色体上, 与前人研究结果一致。同时其加性效应均为负, 说明抗病QTL的抗病性贡献都来自抗性较强的父本P174, 在抗性较弱亲本中没有检测到相关QTL, 这可能和遗传图谱的覆盖度有很大关系, 不能否定抗性较弱亲本中存在抗性QTL。在A05染色体上检测到1个有重叠置信区间的菌核病抗性QTL (q2013A05LS-2和q2014A05LS-2), 且抗性增效基因均来自父本P174, 可能为同一QTL。Tuberosa等^[41]认为性状相关表型在QTL分析结果上, 可能存在控制不同性状QTL之间的紧密连锁或同一QTL位点控制不同性状的现象即一因多效。Zhao等^{[3, 9]}发现油菜菌核病抗性QTL与种子硫甙含量及油菜开花时间QTL间有重叠区间。梅家琴^[12]在F₂无性系群体中发现油菜菌核病抗性QTL与开花期QTL区间存在部分区域重叠。本研究发现油菜菌核病抗性相关QTL (q2013A05LS-1)与茎秆木质素含量相关QTL (q2013A05Lig)的置信区间部分重叠, 佐证了菌斑大小与木质素含量极显著负相关性。上述研究结果表明在油菜抗菌核病QTL中也存在一因多效的现象。
甘蓝型油菜木质素的遗传基础报道相对较少, 多数研究认为木质素遗传基础较复杂, 表型易受外界环境影响, 属于多基因控制的数量性状。Snowdon等^[29]在甘蓝型油菜染色体A09上发现1个种皮木质素含量相关的主效QTL。Liu等^[30]利用BSA法对RIL和BC₁F₂群体进行SSR标记分析, 检测到1个与种子木质素含量相关的主效QTL, 位于A09连锁群上。Liu等^[32]以4个环境下RIL9为材料, 对甘蓝型油菜种皮木质素含量进行定位, 共检测到8个QTL, 所有环境中均能稳定检测到2个主效QTL, 分别位于A09和C05上, 解释的表型变异率分别为31.6%~42.8%和8.1%~14.1%。前人关于油菜种皮木质素含量数量性状位点研究结果表明, 相关QTL主要集中在A09和C05上, 且多次在A09上检测到主效QTL, 由于不同遗传图谱上利用的分子标记不同, 因此还无法确定上述各研究中的QTL是否相同。本研究中检测到3个木质素含量QTL, 分别位于A04、A05和C01染色体, 单个QTL可解释表型变异的2.03%~13.75%。与前人对种子或种皮木质素含量相关QTL定位结果不一致, 推测种子或种皮木质素含量和茎秆木质素含量可能受不同基因控制, 因此它们在遗传上可能相互独立。本试验共检测到6个与木质素单体G/S相关的QTL, 分布在A08、C03和C07染色体上, 单个QTL可解释表型变异的2.06%~ 8.66%。其中q2013C07LigG/S-1和q2013C07LigG/S-2的置信区间重叠, q2013C07LigG/S-2和q2014C07 LigG/S置信区间部分相同, 3个QTL的加性效应均来自父本, 因此可能是同1个QTL。此外, 这种影响同一性状的不同QTL聚集在C07染色体的现象, 也可能是同一位点的复等位基因在多倍体形成过程中发生染色体重排所致。目前与木质素单体G/S相关的QTL尚未见报道, 因此还不能与前人研究作比较。
4 结论茎秆菌斑大小与木质素含量呈极显著负相关, 与单体G/S呈显著正相关。木质素含量以菌核病抗性较强的材料高于抗性较弱材料。共获得18个QTL, 分别位于A04、A05、A06、A08、C01、C03、C04、C06和C07染色体, 含9个菌核病抗性相关的QTL, 3个与木质素含量相关的QTL, 6个与单体G/S相关的QTL。其中q2013A05LS-2和q2014A05LS-2置信区间基本重合, 且抗性增效基因均来自父本, 可能是同1个QTL。q2013A05Lig和q2013A05LS-1的置信区间部分重叠。q2013C07LigG/S-1和q2013C07 LigG/S-2的置信区间相重叠, q2013C07LigG/S-2和q2014C07LigG/S的置信区间有相同部分, 且3个QTL的加性效应均来自父本, 因此也可能是同1个QTL。
The authors have declared that no competing interests exist.


参考文献View Option
原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1]	Carr R A, McDonald B E. Rapeseed in a changing world: Processing and utilization. GCIRC Eighth International Rapeseed Congress. Saskatoon, 1991. pp 39-56[本文引用:1]
[2]	李加纳, 谌利, 张学昆. 甘蓝型黄籽油菜的研究与思考. 北京: 中国农业科技出版社, 2004. pp 29-39 Li J N, Chen L, Zhang X K. Research and thinking of yellow- seeded rapeseed (Brassica napus L. ). Beijing: China Science and Technology Press, 2004. pp 29-39(in Chinese)[本文引用:1]
[3]	Zhao J, Meng J. Genetic analysis of loci associated with partial resistance to Sclerotina sclerotiorum in rapeseed (Brassica napus L. ). Theor Appl Genet, 2003, 106: 759-764[本文引用:3]
[4]	Garg H, Sivasithamparam K, Banga S. Cotyledon assay as a rapid and reliable method of screening for resistance against Sclerotinia sclerotiorum in Brassica napus genotypes. Australas Plant Path, 2008, 37: 106-111[本文引用:1]
[5]	Bradley C, Hamey H. Canola disease situation in North Dakota, USA, 1993-2004. 14th Australian Research Assembly on Brassicas. Port Lincoln, 2005. pp 33-34[本文引用:1]
[6]	Turkington T K, Morrall R A A. Use of petal infestation to forecast Sclerotinia stem rot of canola: the influence of inoculum variation over the flowering period and canopy density. Phytopathology, 1993, 83: 682-689[本文引用:1]
[7]	del Río L E, Bradley C A, Henson R A, Endres G J, Hanson B K, McKay K, Halvorson M, Porter P M, Le Gare D G, Lamey H A. Impact of Sclerotinia stem rot on yield of canola. Plant Dis, 2007, 91: 191-194[本文引用:1]
[8]	王汉中, 刘贵华, 郑元本, 王新发, 杨庆. 抗菌核病双低油菜新品种中双9号选育及其重要防御酶活性变化规律的研究. 中国农业科学, 2004, 37: 23-28 Wang H Z, Liu G H, Zheng Y B, Wang X F, Yang Q. Breeding of the Brassica napus cultivar zhongshuang 9 with high-resistance to Sclerotinia sclerotiorum and dynamics of its important defense enzyme activity. Sci Agric Sin, 2004, 37: 23-28 (in Chinese with English abstract)[本文引用:1]
[9]	Zhao J, Udall J A, Quijada P A, Grau C R, Meng J, Osborn T C. Quantitative trait loci for resistance to Sclerotinia sclerotiorum and its association with a homeologous non-reciprocal transposition in Brassica napus L. Theor Appl Genet, 2006, 112: 509-516. [本文引用:3]
[10]	Wu J, Cai G, Tu J, Li L, Liu S, Luo X, Zhou L, Fan C, Zhou Y. Identification of QTLs for resistance to Sclerotinia stem rot and BnaC. IGMT5. a as a cand idate gene of the major resistant QTL SRC6 in Brassica napus. PloS One, 2013, 8: e67740[本文引用:2]
[11]	周李鹏. 甘蓝型油菜抗菌核病QTL定位. 华中农业大学硕士学位论文, 湖北武汉, 2014. pp 20-21 Zhou L P. QTL Mapping for Resistance to Sclerotinia Stem Rot in Brassica napus. MS Thesis of Huazhong Agricultural University, Wuhan, China, 2014. pp 20-21 (in Chinese with English abstract)[本文引用:2]
[12]	梅家琴. 甘蓝与甘蓝型油菜C亚基因组遗传关系调查及甘蓝抗菌核病QTL定位. 西南大学博士学位论文, 重庆, 2011. pp 62-73 Mei J Q. Genetic Investigation of Relationships between Brassica oleracea and C Subgenome of B. napus and Mapping QTL for Sclerotinia sclerotiorum Resistance in B. oleracea. PhD Dissertation of Southwest University, Chongqing, China, 2011. pp 62-73 (in Chinese with English abstract)[本文引用:3]
[13]	Wei L, Jian H, Lu K, Filardo F, Yin N, Liu L, Qu C, Wei L, Du H, Li J. Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus. Plant Biotechnol J, 2015, 14: 1368-1380[本文引用:2]
[14]	Wu J, Zhao Q, Yang Q, Liu H, Li Q, Yi X, Cheng Y, Guo L, Fan C, Zhou Y. Comparative transcriptomic analysis uncovers the complex genetic network for resistance to Sclerotinia sclerotiorum in Brassica napus. Sci Rep, 2016, 6: 19007[本文引用:1]
[15]	Hoffman D D, Diers B W, Hartman G L, Nickell C D, Nelson R L, Pedersen W L, Cober E R, Graef G L, Steadman J R, Grau C R, Nelson B D, del Rio L E, Helms T, Anderson T, Poysa V, Rajcan I, Stienstra W C. Selected soybean plant introductions with partial resistance to Sclerotinia sclerotiorum. Plant Dis, 2002, 86: 971-980[本文引用:1]
[16]	Kim H S, Diers B W. Inheritance of partial resistance to Sclerotinia stem rot in soybean. Crop Sci, 2000, 40: 55-61[本文引用:1]
[17]	余叔文, 汤章城. 植物生理与分子生物学. 北京: 科学出版社, 2002. pp 770-781 Yu S W, Tang Z C. Plant Physiology and Molecular Biology. Beijing: Science Press, 2002. pp 770-781(in Chinese)[本文引用:1]
[18]	Boudet A M, Lapierre C, Grima-Pettenati J. Biochemistry and molecular biology of lignification. New Phytologist, 1995, 129: 203-236[本文引用:1]
[19]	Nicholson R L, Hammerschmidt R. Phenolic compounds and their role in disease resistance. Annu Rev Phytopathol, 1992, 30: 369-389[本文引用:1]
[20]	Vance C P, Kirk T K, Sherwood R T. Lignification as a mechanism of disease resistance. Annu Rev Phytopathol, 1980, 18: 259-288[本文引用:1]
[21]	Dushnicky L G, Ballance G M, Sumner M J, MacGregor A W. The role of lignification as a resistance mechanism in wheat to a toxin-producing isolate of Pyrenophora tritici-repentis. Can J Plant Pathol, 1998, 20: 35-47[本文引用:1]
[22]	Hammerschmidt R, Kuć J. Lignification as a mechanism for induced systemic resistance in cucumber. Physiol Plant Pathol, 1982, 20: 61-71[本文引用:1]
[23]	Southerton S G, Deverall B J. Histochemical and chemical evidence for lignin accumulation during the expression of resistance to leaf rust fungi in wheat. Physiol Mol Plant Pathol, 1990, 36: 483-494[本文引用:1]
[24]	杨向东. 木质素合成调控及其与甘蓝型油菜抗菌核病和抗倒伏性关系研究. 华中农业大学博士学位论文, 湖北武汉, 2006. pp 28-30 Yang X D. The Study on Relationship between Lignin Biosynthesis Manipulation and Brassica napus’ Resistance to Sclerotinia sclerotiorum and Lodging. PhD Dissertation of Huazhong Agricultural University, Wuhan, China, 2006. pp 28-30 (in Chinese with English abstract)[本文引用:2]
[25]	Gayoso C, Pomar F, Novo-Uzal E, Merino F, de Ilárduya Ó M. The Ve-mediated resistance of the tomato to Verticillium dahliae involves H2O2, peroxidase and lignins and drives PAL gene expression. BMC Plant Biol, 2010, 10: 1[本文引用:2]
[26]	Hammerschmidt R, Bonnen A M, Bergstrom G C, Baker K K. Association of epidermal lignification with nonhost resistance ofcucurbits to fungi. Can J Bot, 1985, 63: 2393-2398[本文引用:2]
[27]	Pomar F, Novo M, Bernal M A, Merino F, Barceló A R. Changes in stem lignins (monomer composition and crosslinking) and peroxidase are related with the maintenance of leaf photosynthetic integrity during Verticillium wilt in Capsicum annuum. New Phytol, 2004, 163: 111-123[本文引用:2]
[28]	Barceló A R. Lignification in plant cell walls. Int Rev Cytol, 1997, 176: 87-132[本文引用:2]
[29]	Snowdon R J, Wittkop B, Rezaidad A, Hasan M, Lipsa F, Stein A, Friedt W. Regional association analysis delineates a sequenced chromosome region influencing antinutritive seed meal compounds in oilseed rape. Genome, 2010, 53: 917-928[本文引用:2]
[30]	Liu L, Stein A, Wittkop B, Sarvari P, Li J, Yan X, Dreyer F, Frauen M, Friedt W, Snowdon R J. A kockout mutation in the lignin biosynthesis gene CCR1 explains a major QTL for acid detergent lignin content in Brassica napus seeds. Theor Appl Genet, 2012, 124: 1573-1586[本文引用:2]
[31]	曲存民, 付福友, 卢坤, 谢景梅, 刘晓兰, 黄杰恒, 李波, 王瑞, 谌利, 唐章林, 李加纳. 不同环境中甘蓝型油菜种皮木质素含量的QTL定位. 作物学报, 2011, 37: 1398-1405 Qu C M, Fu F Y, Lu K, Xie J M, Liu X L, Huang J H, Li B, Wang R, Chen L, Tang Z L, Li J N. Identification of QTLs for lignin content of seed coat in Brassica napus L. in different environments. Acta Agron Sin, 2011, 37: 1398-1405 (in Chinese with English abstract)[本文引用:1]
[32]	Liu L Z, Qu C M, Wittkop B, Yi B, Xiao Y, He H Y, Snowdon R, Li J. A high-density SNP map for accurate mapping of seed fibre QTL in Brassica napus L. PLoS One, 2013, 8: e83052[本文引用:3]
[33]	Chen F, Dixon R A. Lignin modification improves fermentable sugar yields for biofuel production. Nat Biotechnol, 2007, 25: 759-761[本文引用:1]
[34]	黄杰恒. 干旱胁迫下油菜抗倒伏相关性状动态变化及木质素关键基因表达特性分析. 西南大学博士学位论文, 重庆, 2013. pp 42-43 Huang J H. Lodging Resistant Traits and Lignin Related Gene Analysis in B. napus under Drought Stress. PhD Dissertation of Southwest University, Chongqing, China, 2013. pp 42-43 (in Chinese with English abstract)[本文引用:2]
[35]	Eynck C, Séguin-Swartz G, Clarke W E, Parkin I A. Monolignol biosynthesis is associated with resistance to Sclerotinia sclerotiorum in Camelina sativa. Mol Plant Pathol, 2012, 13: 887-899[本文引用:3]
[36]	McCouch S R, Cho Y G, Yano M, Paul E, Blinstrub M, Morishima H, Kinoshita T. Report on QTL nomenclature. Rice Genet Newsl, 1997, 14[本文引用:1]
[37]	Li J, Zhao Z, Hayward A, Cheng H, Fu D. Integration analysis of quantitative trait loci for resistance to Sclerotinia sclerotiorum in Brassica napus. Euphytica, 2015, 205: 483-489[本文引用:1]
[38]	Behla R S, Fernand o W G D, Li G. Identification of quantitative trait loci for resistance against Sclerotinia stem rot in Brassica napus. Can J Plant Pathol, 2009, 31: 477-478[本文引用:1]
[39]	Yin X R, Yi B, Chen W, Zhang W J, Tu J X, Fernand o W G D, Fu T D. Mapping of QTLs detected in a Brassica napus DH population for resistance to Sclerotinia sclerotiorum in multiple environments. Euphytica, 2010, 173: 25-35[本文引用:1]
[40]	Wei D, Mei J, Fu Y, Disi J O, Li J, Qian W. Quantitative trait loci analyses for resistance to Sclerotinia sclerotiorum and flowering time in Brassica napus. Mol Breed, 2014, 34: 1797-1804[本文引用:1]
[41]	Tuberosa R, Salvi S, Sanguineti M C, Land i P, Maecaferri M, Conti S. Mapping QTLs regulating morpho-physiological traits and yield: case studies, shortcomings and perspectives in drought-stressed maize. Ann Bot, 2002, 89: 941-963[本文引用:1]