侯林涛^**, 王腾岳^**, 荐红举, 王嘉, 李加纳, 刘列钊^*
西南大学农学与生物科技学院, 重庆 400715
^*通讯作者(Corresponding author): 刘列钊, E-mail: liezhao2003@126.com, Tel: 023-68250701 ^**同等贡献(Contributed equally to this work) 第一作者联系方式: E-mail: 594632283@qq.com
收稿日期:2016-04-20 接受日期:2016-09-18网络出版日期:2016-09-27基金:本研究由国家自然科学基金项目(31371655)资助

摘要盐胁迫是非生物胁迫中影响作物产量的一个主要因素, 利用分子标记方法选育油菜耐盐品种对提高油菜产量具有重要意义。选用来自GH06与P174杂交后通过单粒传法连续自交获得的高世代重组自交系群体, 以含16 g L^-1 NaCl的Hoagland溶液培养幼苗进行盐胁迫处理25 d后, 分别测定叶和根的鲜重及干重, 根据已构建的高密度SNP遗传连锁图谱进行QTL定位, 在QTL物理区间筛选耐盐相关基因并以极端表型材料进行qRT-PCR分析。采用复合区间作图法(CIM), 在对照和盐胁迫处理中共检测到19个QTL, 其中与盐胁迫相关的有6个, 可解释的表型变异7.16%~16.15%, 分布在A02、A04和C03染色体上, 将QTL置信区间序列和拟南芥中与盐胁迫相关的基因比对分析, 共找到8个候选基因。对其中4个候选基因在极端表型材料中的表达分析表明, BnaA02g14680D与 BnaA02g14490D基因在盐胁迫处理后的48 h或72 h表达量均高于对照组, 即基因的表达由盐胁迫引起, 而 BnaC03g64030D在敏感型材料中的相对表达量高于在耐盐型材料中, BnaC03g62830D在敏感型材料中没有明显变化, 但在耐盐型材料中呈现先升高后降低的表达特征, 其表达可能会增强植株对盐胁迫的耐受力。本研究为油菜耐盐基因功能挖掘和油菜耐盐品种选育奠定基础。

关键词:甘蓝型油菜; 盐胁迫; 数量性状位点; 单核苷酸多态性; 候选基因
QTL Mapping for Seedling Dry Weight and Fresh Weight under Salt Stress and Candidate Genes Analysis in Brassica napusL
HOU Lin-Tao^**, WANG Teng-Yue^**, JIAN Hong-Ju, WANG Jia, LI Jia-Na, LIU Lie-Zhao^*
College of Agronomy and Biotechnology, Southwest University, Chongqing 400715, China
Fund:The study was supported by the National Natural Science Foundation of China (31371655)
AbstractSalt stress is one of the main abiotic stresses affecting crop yield and it would be very important by using the salt tolerance related markers in rapeseed breeding to improve the oilseed production. In this research, the Brassica napus L. high generation recombinant inbred lines (RIL) population derived from the cross of GH06 and P174 via single seed descent propagation was used for QTL mapping and candidate gene analysis. The fresh and dry weight of leaf and root were measured at 25 days after the seedlings were grown in Hoagland solution with 16 g L^-1NaCl. Composite interval mapping (CIM) was used to identify the related QTLs according to the high density SNP genetic map, and the candidate gene expression in the extreme lines tested by qRT-PCR. A total of 19 QTLs were identified in the control and salt stress treatment, and six QTLs were mapped on chromosomes A02, A04, and C03 under salt stress, with contribution rate ranged from 7.16% to 16.15%. Eight genes were identified according to the BLAST of genes in the QTL confidence intervals and the salt stress related genes in Arabidopsis. The expression of four candidate genes in the extreme lines showed that BnaA02g14680D and BnaA02g14490D under salt stress treatment for 48 or 72 hours had higher expression than the control, which indicates that the expressions are induced by salt stress. The relative expressions of gene BnaC03g64030Din sensitive extreme lines were higher than those in tolerant extreme lines. There were no changed in expression for gene BnaC03g62830D in sensitive extreme lines but increased expression at 48 hours and reduced expression at 72 hours after salt treatment in tolerant extreme lines, showing the enhance of plant salt tolerance possibly. Our research laid a foundation for the function research of salt tolerant gene in rapeseed and the breeding of salt tolerant rapeseed.

Keyword: Brassica napusL; Salt stress; QTL; SNP; Candidate gene
Show Figures
Show Figures










土壤盐渍化是危害农业生产的主要因素之一, 中国盐渍土面积大、分布广、类型多, 总面积达9913万公顷, 约占国土面积的1.03% ^[1]。土壤中过量的盐分会引起土壤物理和化学性质的改变, 导致农作物生长环境的恶化^[2]。在我国, 油菜是继水稻、玉米、小麦、大豆之后的第五大作物, 也是第一大油料作物, 常年种植面积约670万公顷, 总产超过1100万吨, 年产菜籽油约450万吨, 占自产植物油总量的40%以上, 占草本植物油总量的57%以上^[3]。甘蓝型油菜在三大栽培类型油菜中种植面积最大, 被认为是比较耐盐的作物之一^[4]。因此我们期望通过甘蓝型油菜耐盐性定位研究为其耐盐性和品种选育提供研究基础, 进而实现盐碱地的农业高效利用, 这对我国耕地农业生产能力的提升, 耕地数量的增加, 国家粮食安全的保障具有重要意义^[5]。
随着对油菜研究的不断深入及SNP芯片和重测序的普及, 利用SNP分子标记技术建立高密度的遗传图谱, 对油菜重要性状的QTL分析已经成为新常态, 并且已经定位了与油菜含油量、千粒重、株高、产量等相关的QTL ^{[6, 7, 8, 9, 10]}, 对于油菜在盐胁迫条件下的研究大多集中在生理反应的测定^{[11, 12, 13]}, 而疏于对油菜耐盐性状基因的定位。Li等^[14]利用甘蓝型油菜DH系群体在正常和干旱胁迫条件下的5个指标(株高、根长、叶干重、根干重和总重)数据的QTL分析, 分别检测到28个和31个QTL, 其中在正常环境下, 以油菜幼苗干重为指标检测到18个QTL, 单个位点可解释的表型变异为6.43%~18.97%; 旱胁迫环境下, 检测到与幼苗干重相关的20个QTL, 单个位点可解释的表型变异为5.87%~15.97%。荐红举等^[15]利用本研究的RIL群体对盐胁迫下种子发芽率的QTL分析, 检测到11个相应的QTL, 位于A03和A09染色体上, 可解释的表型变异为4.9%~10.9%。Moursi^[16]利用DH群体在正常和盐胁迫环境下的干重和鲜重数据进行QTL分析, 正常环境下找到4个QTL, 分别位于A06、A05、C02和C03染色体上, 贡献率为36%; 盐胁迫环境下找到3个QTL, 分别位于C03和C06染色体上, 共同可解释的遗传变异25%。
本试验利用已经构建的高密度SNP遗传连锁图谱对盐胁迫环境下RIL群体叶和根的鲜重和干重进行QTL定位, 进而筛选出与盐胁迫相关的候选基因并分析其在极端表型材料中的表达量, 为油菜耐盐基因的进一步研究提供理论基础。
1 材料与方法1.1 试验材料以母本GH06与父本P174杂交, 采用单粒传法连续自交10代, 构建高世代重组自交系群体(recombinant inbred line, RIL), 之后随机选取其中的172个重组自交系进行SNP标记分析, 构建高密度SNP遗传连锁图谱^[17]。
1.2 盐胁迫处理挑选健康饱满、大小均一的油菜种子, 置铺2层定性滤纸的9 cm培养皿内, 每皿20粒, 加入适量的1/4 Hoagland营养液培养7~8 d。待种子发芽长成幼苗时, 从每株系各取生长健康一致的1植株, 3个重复, 分别置装有质量浓度为0 g L^-1和16 g L^-1NaCl溶液^[15](以Hoagland营养液为母液配制)的塑料盆中, 于温室培养25 d, 设定温度25℃, 光照16 h, 黑暗8 h。
1.3 幼苗鲜重、干重测定用滤纸吸干每株系3个植株水分后分别用电子天平测量幼苗的叶鲜重、根鲜重, 取平均值。将幼苗放入105℃烘箱杀青30 min, 70℃烘箱烘干48 h至恒重, 然后用电子天平称量幼苗的叶干重、根干重, 取平均值。
1.4 图谱构建及QTL分析所用SNP分子标记构建的遗传连锁图谱^[17], 包含2795个SNP多态性标记位点, 总长1832.9 cM, 相邻标记间平均距离为0.66 cM。采用QTL分析软件Windows QTL Cartographer 2.5^[18]及复合区间作图(composite interval mapping, CIM)法对RIL群体幼苗的叶鲜重、根鲜重与叶干重、根干重进行QTL定位及效应检测^[19]。CIM分析时, 选取1 cM的步长(walking speed), 按照假定检测10和Zmapqtl模型3, 选取参数1000次回归, 显著水平为0.01。LOD≥ 3时, 即认为该区间可能存在一个QTL。运行软件同时给出性状QTL的加性效应和解释的表型变异。按照McCouch等^[20]的方法命名检测到的QTL, 以“ q” 加相对应的性状再加染色体编号表示, 字体为斜体。
1.5 候选基因的筛选将检测到的QTL的置信区间在甘蓝型油菜基因组^[21]上查询到对应的序列, 然后与拟南芥基因组序列进行BlastN, E值设定为E-20, 最后筛选出每个QTL置信区间内匹配E值小于阈值的盐胁迫相关基因, 共筛选出8个与盐胁迫相关的基因。
1.6 RNA提取根据RIL群体叶、根鲜重在盐胁迫环境与正常环境下的比值, 筛选出对盐敏感型材料(F215的比值分别为0.78和0.42, F328分别为0.76和0.37)和耐盐型(F285的比值为分别1.02和1.50, F293分别为1.00和1.38)的极端表型材料, 在盐胁迫处理24、48、72和96 h后取叶片, 将样品迅速于-80℃保存。按照EASYspin植物RNA快速提取试剂盒(RN09)说明提取样品的总RNA。采用Bio-Rad iScript cDNA Synthesis Kit试剂盒反转录, 反转录体系为5× iScript reaction mix 4 μ L, iScript reverse transcriptase 1 μ L, RNA模板(1 μ g) 1~3 μ L, ddH₂O补充体积至20 μ L。反应条件为25℃ 5 min, 42℃ 30 min, 85℃ 5 min, 4℃保存。
1.7 实时荧光定量PCR(qRT-PCR)在甘蓝型油菜基因组^[21]中查找到候选基因的cDNA序列, 利用软件Primer Premier 5.0设计特异引物, 由上海生工生物工程股份有限公司合成(表1)。按照Bio-Rad公司的SsoAdvanced Universal SYBR Green Supermix 试剂盒使用说明在Bio-Rad CFX96上运行实时荧光定量PCR。反应体系为SsoAdvanced universal SYBR Green supermix (2× ) 10 μ L, 正向引物0.4 μ L, 反向引物0.4 μ L, cDNA 2 μ L, ddH₂O补充体积至20 μ L。反应初始变性温度和时间为95℃/30 s, 然后95℃/10 s, 退火时间为30 s, 温度因引物而异, 40个循环。熔解曲线为65~95℃之间以0.5℃递增, 5 s。所有反应3个重复。用2^{-Δ Δ CT}方法以内参(油菜Actin基因-AF111812)为标准计算每个基因的相对表达量。
表1
Table 1
表1(Table 1)

表1 荧光定量PCR引物 Table 1 Primers used for RT-PCR 基因名称 Gene name 同源基因 Homologous gene 上游引物序列 Forward primer sequence (5'-3') 下游引物序列 Revers primer sequence (5'-3') 熔解温度 Tm (℃) 扩增片段大小 Fragment size (bp) BnaA02g14680D ATRPK1/RPK1 CTCGTCCTCGTTGTCTTG TCGTAGGTCAGAGGGATA 54.0 118 BnaA02g14490D ATHVA22C TTCCGATATGGCCGTTCT TTCGCTCAAATGCGTCAG 60.2 245 BnaC03g64030D ATCPL1/SHI4 TCCTCTTCATCCTACAGTCTCCAAA TAGTCCACCTTAGTTCCA 60.2 263 BnaC03g62830D ATCESA8/LEW2 AAAGCGATGAGTTTCAACAAGACCA GTCCGTATTTCTCCACAA 60.2 138

表1 荧光定量PCR引物 Table 1 Primers used for RT-PCR

1.8 统计数据利用SPSS 17.0软件完成相关性分析。利用MapChart 2.2作图软件分析染色体和QTL之间的遗传距离, 在Microsoft Excel 2003和WPS 2016中处理其他数据。

2 结果与分析2.1 油菜幼苗在盐胁迫条件下鲜重及干重的变化RIL群体的鲜重、干重, 各数值均呈连续分布和双向超亲分离(表2、图1和图2), 表明这些性状均为多基因控制的数量性状, 其偏斜度和峰值均小1, 适合于QTL分析。除了盐胁迫处理下根干重均值高于对照组外, 其余叶、根的鲜重、干重均值均低于对照组, 这表明盐胁迫环境对幼苗的生长和发育产生了影响, 而根干重高于对照组可能由于幼苗在盐胁迫环境下促进根系的生长。盐胁迫下的叶鲜重均值要远远低于对照组, 而叶干重却相差不大, 可能是由于盐胁迫条件下离子浓度过高而影响植株的水分吸收。
表2
Table 2
表2(Table 2)

表2 油菜幼苗在盐胁迫处理和对照处理下干重、鲜重的分布特征 Table 2 Phenotype variation of dry weight and fresh weight under salt stress and control conditions 重量 Weight 亲本 Parent RIL群体 RIL population GH06 (mg) P174 (mg) 范围 Range (mg) 均值 Mean (mg) 标准差 SD 变异系数 CV 偏度 Skewness 峰度 Kurtosis 对照叶鲜重 CK-Leaf-FW 536.93 118.90 113.8-921.0 404.92 164.81 0.41 0.78 0.76 盐胁迫叶鲜重 Salt-Leaf-FW 340.10 607.30 82.1-607.3 293.99 89.52 0.30 0.62 0.90 对照根鲜重 CK-Root-FW 27.53 11.50 8.3-47.7 25.71 8.55 0.33 0.27 -0.44 盐胁迫根鲜重 Salt-Root-FW — 51.00 8.7-51.0 24.58 7.99 0.33 0.79 1.35 对照叶干重 CK-Leaf-DW 25.23 48.00 2.2-48.6 19.21 8.30 0.43 0.90 1.22 盐胁迫叶干重 Salt-Leaf-DW 23.85 36.70 6.4-39.3 18.70 5.76 0.31 0.86 1.15 对照根干重 CK-Root-DW 2.93 1.55 0.66-4.10 2.18 0.66 0.30 0.37 -0.17 盐胁迫根干重 Salt-Root-DW 2.20 4.50 0.5-4.5 2.39 0.70 0.29 0.34 0.56

表2 油菜幼苗在盐胁迫处理和对照处理下干重、鲜重的分布特征 Table 2 Phenotype variation of dry weight and fresh weight under salt stress and control conditions

图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 油菜幼苗叶片和根鲜重在盐胁迫和对照下的频率分布图Fig. 1 Frequency distribution of fresh weight of B. napusunder control and salt stress conditions

图2
Fig. 2

	Figure OptionViewDownloadNew Window
	图2 油菜幼苗干重在盐胁迫处理和对照处理下的频率分布图Fig. 2 Frequency distribution of dry weight of B. napus under salt stress and control conditions

2.2 油菜幼苗在2种环境下鲜重、干重的相关性分析从表3可以看出2种处理下的叶干重显著正相关, 相关系数为0.205。盐胁迫环境下叶干重和根干重极显著正相关, 相关系数为0.771, 盐胁迫环境下叶鲜重和根鲜重极显著相关, 相关系数为0.380。相关性分析表明, 相同的处理下, 相同部位干重和鲜重均极显著相关。在对照和处理间, 只有叶片的干重和鲜重表现显著和极显著相关, 相关系数分别为0.205和0.246, 而根部在对照和处理间相关性不显著。
表3
Table 3
表3(Table 3)

表3 油菜幼苗盐胁迫处理和对照干重、鲜重的相关系数 Table 3 Correlation coefficients between indication of dry weight and fresh weight of B. napus in salt stress treatment and control 指标 Indicator 对照叶干重 CK-Leaf-DW 对照根干重 CK-Root-DW 盐胁迫叶干重 Salt-Leaf-DW 盐胁迫根干重 Salt-Root-DW 对照叶鲜重 CK-Leaf-FW 盐胁迫叶鲜重 Salt-Leaf-FW 对照根鲜重 CK-Root-FW 对照根干重 CK-Root-DW 0.749* * 叶干重 Salt-Leaf-DW 0.205* 0.165 根干重 Salt-Root-DW 0.052 0.041 0.771* * 对照叶鲜重 CK-Leaf-FW 0.957* * 0.740* * 0.246* * 0.057 叶鲜重 Salt-Leaf-FW 0.120 0.027 0.487* * 0.429* * 0.129 对照根鲜重 CK-Root-FW 0.858* * 0.789* * 0.162 0.037 0.843* * 0.080 根鲜重 Salt-Root-FW 0.087 0.067 0.870* * 0.852* * 0.095 0.380* * 0.078 * * and * denote significant correlation at the 0.01 and 0.05 probability levels, respectively. * * 和* 分别代表在0.01和0.05显著水平。

表3 油菜幼苗盐胁迫处理和对照干重、鲜重的相关系数 Table 3 Correlation coefficients between indication of dry weight and fresh weight of B. napus in salt stress treatment and control

2.3 正常环境下油菜幼苗干重、鲜重的QTL分析总共检测到13个QTL (表4和图3)。油菜叶和根干重相关的QTL有4个, 分布在染色体A01、A10和C08上, 单个QTL可解释的表型变异为8.43%~ 12.19%。其中C08染色体上的qCKDLC08贡献率为12.19%, LOD值为6.28。油菜幼苗叶和根鲜重的QTL有9个, 单个QTL可解释的表型变异为7.08%~ 15.08%, 分布在A01、A08、A10和C08染色体上。其中分布在C08染色体上的QTL有4个, 贡献率为7.56%~15.08%, 置信区间有部分重叠。检测到A08染色体上3个QTL, 贡献率为8.00%~8.84%, 置信区间相互靠近。在A01染色体上的QTL qCKWLA01和qCKDLA01置信区间完全重合, 推测在此区间内可能有控制这一性状的基因。
表4
Table 4
表4(Table 4)

表4 利用复合区间作图法检测油菜幼苗干重、鲜重在正常环境下的QTL Table 4 Putative QTL detected by composite interval mapping for dry weight and fresh weight of B. napus in normal environment 性状 Trait 数量性状位点 QTL 染色体 Chr. 标记区间 Marker interval (cM) 加性效应 Additive 贡献率 R2 (%) 阈值 LOD 置信区间 Confidence interval 对照叶干重CK-Leaf-DW qCKDLA01 A01 SNP2668A01-SNP2079A01 -27.55 9.44 4.67 92.50-96.20 qCKDLA10 A10 SNP23037A10-SNP23245A10 -25.91 8.43 4.49 60.00-65.60 qCKDLC08 C08 SNP30078C08-SNP30068C08 31.24 12.19 6.28 35.50-36.70 对照根干重CK-Root-DW qCKDRA01 A01 SNP2668A01-SNP1348A01 -2.13 9.13 3.34 92.50-94.60 对照叶鲜重CK-Leaf-FW qCKWLA01 A01 SNP2668A01-SNP2079A01 -664.61 13.38 6.34 92.50-96.20 qCKWLA10 A10 SNP22233A10-SNP23011A10 -493.58 7.08 3.81 57.50-60.00 qCKWLC08-1 C08 SNP30088C08-SNP30043C08 715.82 15.08 6.58 35.30-36.80 qCKWLC08-2 C08 SNP30060-SNP30034C08 701.67 14.91 6.67 37.50-38.10 qCKWLC08-3 C08 SNP31434-SNP30196 599.27 11.66 5.42 26.70-31.80 qCKWLC08-4 C08 SNP30229C08-SNP30122C08 531.73 7.56 3.36 42.30-47.20 对照根鲜重CK-Root-FW qCKWRA08-1 A08 SNP18210A08-SNP18201A08 32.33 8.84 3.28 40.40-41.70 qCKWRA08-2 A08 SNP12599A05-SNP18516A08 -28.81 8.62 3.59 12.50-13.50 qCKWRA08-3 A08 SNP18508A08-SNP18387A08 -29.86 8.00 4.48 14.90-26.80

表4 利用复合区间作图法检测油菜幼苗干重、鲜重在正常环境下的QTL Table 4 Putative QTL detected by composite interval mapping for dry weight and fresh weight of B. napus in normal environment

图3
Fig. 3

	Figure OptionViewDownloadNew Window
	图3 甘蓝型油菜幼苗盐胁迫处理下干重、鲜重QTL在连锁群上的分布情况 QTL的染色体上仅展示了QTL区段的标记及染色体两端各2个标记。Fig. 3 Putative QTL locations of dry weight and fresh weight in salt stress treatment on the SNP genetic map Showing the markers in the QTL confidence intervals, along with the terminal two markers at each end of the QTL containing chromosomes.

2.4 盐胁迫环境下油菜幼苗干重、鲜重的QTL定位检测到与盐胁迫相关的QTL 6个(表5和图4), 阈值为3.21~4.59, 单个QTL可解释的表型变异为7.16%~16.15%, 分布在染色体A02、A04和C03上。分布在A02染色体上的QTL有3个, 单个QTL可解释的表型变异为7.16%~11.75%。A04染色体上1个QTL, LOD值为4.66, 贡献率为11.60%, 置信区间为58.00~77.40。C03染色体上2个QTL, 为qSTWRC03和qSTWRC03, LOD值分别为4.47和4.59, 共同可解释的表型变异是31.75%, 置信区间有部分重叠。
表5
Table 5
表5(Table 5)

表5 利用复合区间作图法检测到油菜幼苗干重、鲜重在盐胁迫环境下的QTL Table 5 Putative QTL detected by composite interval mapping for dry weight and fresh weight of B. napus in salt stress environment 性状 Trait 数量性状位点 QTL 染色体 Chr. 标记区间 Marker interval (cM) 加性效应 Additive 贡献率 R2 (%) 阈值 LOD 置信区间 Confidence interval 盐胁迫叶干重 Salt-Leaf-DW qSTDLA02-1 A02 SNP2849A02-SNP2787A02 -19.97 11.75 4.58 26.30-27.80 qSTDLA02-2 A02 SNP4846A02-SNP4791A02 -15.65 7.16 3.31 33.70-37.30 qSTDLA04 A04 SNP8445A04-SNP8405A04 -26.14 11.60 4.66 58.00-77.40 盐胁迫根鲜重 Salt-Root-FW qSTWRA02 A02 SNP4046A02-SNP11915A05 -26.77 7.16 3.21 48.50-55.00 qSTWRC03-1 C03 SNP39078C03-SNP37474 39.53 16.15 4.47 20.90-28.40 qSTWRC03-2 C03 SNP36113C03-SNP36062C03 38.67 15.60 4.59 29.90-35.30

表5 利用复合区间作图法检测到油菜幼苗干重、鲜重在盐胁迫环境下的QTL Table 5 Putative QTL detected by composite interval mapping for dry weight and fresh weight of B. napus in salt stress environment

图4
Fig. 4

	Figure OptionViewDownloadNew Window
	图4 甘蓝型油菜4个基因在盐胁迫下表达量的比值Fig. 4 Expression ratio of four genes under salt stress (ST) and normal environment (CK) in B. napus

2.5 候选基因筛选在盐胁迫环境下共检测到6个QTL, 但由于qSTWRA02和qSTWRC03-2在甘蓝型油菜基因组中没有找到相应的序列, 故只能将4个QTL置信区间序列与拟南芥中盐胁迫相关基因比对, 共检测到8个候选基因(表6和表7)。其中在qSTDLA02-1区间内1个, 匹配E值为E-137。qSTDLA02-2区间内3个, 匹配E值介于0~2E-52之间。qSTDLA04区间内只1个, 匹配E值为5E-25。qSTWRC03-1区间内3个, 匹配E值为4E-35~E-167。
表6
Table 6
表6(Table 6)

表6 甘蓝型油菜基因组中QTL置信区间候选基因与拟南芥盐胁迫相关基因的比对 Table 6 Alignment of candidate genes in QTL confidence interval in B. napus with the salt related genes in A. thaliana 数量性状座位 QTL 盐胁迫相关基因 Genes related to salt stress inA. thaliana 性状 Trait 名称 Name 物理区间 Physical interval 预测基因 Gene prediction 基因名 Gene name 登录号 Accessions number E值 E-value 盐胁迫叶干重 Salt-Leaf-DW qSTDLA02-1 11237912-10815211 BnaA02g18520D SKIP/SNW AT1G77180 E-137 qSTDLA02-2 8466214-8090196 BnaA02g14680D ATRPK1/RPK1 AT1G69270 0 qSTDLA02-2 8466214-8090196 BnaA02g14490D ATHVA22C AT1G69700 2E-52 qSTDLA02-2 8466214-8090196 BnaA02g14430D ATHB29/ZFHD1 AT1G69600 E-138 qSTDLA04 11329756-10914108 Not found ATGRDP2 AT4G37900 5E-25 盐胁迫根鲜重 Salt-Root-FW qSTWRC03-1 53543616-51523070 BnaC03g64020D BnaC03g64030D ATCPL1/SHI4 AT4G21670 E-141 qSTWRC03-1 53543616-51523070 BnaC03g63430D BnaC03g63440D CHIC AT4G19810 E-167 qSTWRC03-1 53543616-51523070 BnaC03g62830D BnaC03g62840D ATCESA8/LEW2 AT4G18780 4E-35

表6 甘蓝型油菜基因组中QTL置信区间候选基因与拟南芥盐胁迫相关基因的比对 Table 6 Alignment of candidate genes in QTL confidence interval in B. napus with the salt related genes in A. thaliana

表7
Table 7
表7(Table 7)

表7 筛选的在拟南芥中盐胁迫相关候选基因功能 Table 7 Functions of the screened candidate genes related to salt stress in A. thaliana 名称 Name 基因 Gene 基因功能 Gene function 参考文献 Reference qSTDLA02-1 SKIP/SNW 编码一个假定的转录因子, 在酵母中表现转录活性, 参与到对脱落酸、盐等渗透性胁迫反应中。 Encodes a putative transcriptional factor. Shows transcriptional activator activity in yeast. Involved in response to abscisic acid, salt and osmotic stress. Feng et al.[22] qSTDLA02-2 ATRPK1/RPK1 RPK1是一个富含亮氨酸的受体酶, 位于离子膜中, 由脱落酸、脱水、高盐、低温等引起, 和其他植物激素无关。 RPK1 is a leucine-rich receptor-like kinase located in the plasma membrane which is upregulated by abscisic acid, dehydration, high salt, low temperature, but not by other plant hormones. Hong et al.[23] ATHVA22C 它是AtHVA22 家族中的一员, 蛋白表达是由脱落酸和胁迫引起的。 Part of the AtHVA22 family. Protein expression is ABA- and stress-inducible. Guo et al.[24] ATHB29/ZFHD1 编码 ZFHD1, 它是锌指同源转录因子家族的一个成员, 其表达是由干旱、高盐和脱落酸引起的。 Encodes ZFHD1, a member of the zinc finger homeodomain transcriptional factor family. Expression of ZFHD1 is induced by drought, high salinity and abscisic acid. Elfving et al.[25] qSTDLA04 ATGRDP2 编码未知功能的一种蛋白质, 参与植株生长和发育的许多方面以及应对盐胁迫。 Protein of unknown function that expressed in many plant tissues and is involved in many aspects of plant growth and development as well as response to salt stress. Ortega et al.[26] qSTWRC03-1 ATCPL1/SHI4 编码的蛋白质在体外培养具有先天性磷酸酶活性。突变体对脱落酸, 低温和氯化钠十分敏感。 Encodes a protein that exhibits innate phosphatase activity in vitro. Mutants exhibit hyper responsiveness to ABA, cold, and NaCl. Seifert et al.[27] CHIC 反转录PCR实验表明, ChiC的转录水平在应对脱落酸、茉莉酸和氯化钠压力时增加。 RT-PCR experiments demonstrate that ChiC transcript levels are increased in response to abscisisc acid, jasmonic acid, and NaCl stress. Ohnuma et al.[28] ATCESA8/LEW2 编码一个纤维素合酶家族的成员, 其参与次生细胞壁的生物合成。突变体有异常木质部形成, 会降低纤维素含量、增强干旱和渗透压力胁迫。 Encodes a member of the cellulose synthase family involved in secondary cell wall biosynthesis. Mutants have abnormal xylem formation, reduced cellulose content, and enhanced drought and osmotic stress tolerance. Glass et al.[29]

表7 筛选的在拟南芥中盐胁迫相关候选基因功能 Table 7 Functions of the screened candidate genes related to salt stress in A. thaliana

2.6 候选基因在2种环境下的表达水平由图4可知, 在敏感型材料中, BnaA02g14680D基因盐处理后的48 h或72 h表达量较高, 之后96 h表达量降低; 而在耐盐型材料中, BnaA02g14680D基因在盐处理后的表达量先升高后降低之后在96 h又升高。BnaA02g14490D基因与BnaA02g14680D基因的表达水平特征相似。而BnaC03g64030D基因在敏感型材料中盐处理后的表达水平都有急剧升高的时期, 之后下降; 而其在耐盐型材料中的表达量则没有明显变化。BnaC03g62830D基因则在敏感型材料中没有明显变化, 在耐盐型材料中呈现先升高后降低的表达特征。

3 讨论作物生长会受到多种因素的影响, 包括生物胁迫和非生物胁迫, 非生物胁迫中盐胁迫成为限制作物产量的一个重要因素^[30]。油菜的营养价值和经济价值极高, 已经成为植物油的主要来源, 然而其生长、产量和出油量会因为土壤盐渍化而有所降低, 特别是油菜种子发芽和幼苗早期的生长发育更容易受到盐害的影响^[31]。从油菜幼苗在正常环境下和盐胁迫环境下鲜重和干重的分布特征可以看出, 盐胁迫下叶鲜重、干重的分布范围较窄, 均值也低于对照组, 尤其是叶鲜重, 盐胁迫环境下的叶鲜重均值几乎为对照组的一半, 根干重和鲜重则变化不大。推测可能是由于外界环境盐浓度高而影响了幼苗对水分的吸收, 从而导致植株叶鲜重有所下降, 这一结果与Tunç tü rk等^[32]的研究结果一致。
植物对盐胁迫的响应是一个复杂的生理过程^{[33, 34]}, 它会引起很多相关性状的改变, 如种子发芽率、植株鲜重干重等。我们采用SNP高密度遗传连锁图谱对正常环境和盐胁迫环境下的幼苗重量QTL分析, 发现两种环境下检测到的QTL有很大差异, 且分布的染色体也不同, 可能是由于油菜幼苗在不同处理下差异基因适应环境而产生不同的响应。盐胁迫环境下共检测到6个QTL, 其中分布在染色体A02上的有3个, 置信区间十分靠近且有部分重叠, A04染色体上有1个, 加性效应均来自于父本P174; C03染色体上有2个, qSTWRC03-1与qSTWRC03-2共同可解释的表型变异为31.75%, 是盐胁迫环境下找到的贡献率最高的两个QTL, 而在正常环境下并没有检测到, 因此推测它们可能是与盐胁迫环境下幼苗根鲜重相关的主效QTL, 加性效应来自于母本GH06。Moursi^[16]利用DH群体对正常和盐胁迫环境下干重和鲜重得到的数据进行QTL分析, 盐胁迫环境下找到3个QTL, 共同可解释的遗传变异为25%, 分别位于C03和C06染色体上, 但与qSTWRC03-1和qSTWRC03-2位置相距较远。荐红举等^[15]利用本研究的RIL群体对盐胁迫下种子发芽率进行QTL分析, 检测到11个相应的QTL, 可解释的表型变异为4.9%~10.9%, 位于A03和A09染色体上, 与我们检测到的QTL所在染色体不一致, 推测可能是由于油菜在盐胁迫环境下种子发芽和植株生长响应机制不同。
本研究将盐胁迫环境下检测到的4个QTL置信区间与拟南芥盐胁迫相关基因比对, 在每个区间内都找到了相应的候选基因, 共8个, 并对其中4个进行实时荧光定量PCR分析。研究结果表明, BnaA02g14680D基因在拟南芥中对应基因ATRPK1/RPK1的表达是由脱落酸、脱水、高盐和低温引起的, 本实验中该基因在盐胁迫环境下的表达量明显高于对照组, 与Hong等^[23]的研究结果一致。Osakabe等^[35]在实验中也发现在拟南芥中过量表达RPK1基因会增强对盐胁迫的耐受力。但Shi等^[36]认为ATRPK1/RPK1基因的过量表达会使拟南芥对盐胁迫的耐受力降低, 当ATRPK1基因的表达量被限制时, 拟南芥对盐胁迫的耐受力明显升高。BnaA02g14490D基因在拟南芥中对应基因ATHVA22C的表达也是由脱落酸和盐胁迫引起的, 在敏感型和耐盐型材料中的表达特征与BnaA02g14680D基因相似, 在耐盐型材料中两个基因在盐处理后96 h表达量明显升高(图4), 推测可能是由于基因在耐盐型材料中的表达量较高而使其对盐胁迫环境的耐受力增强。BnaC03g64030D基因在拟南芥中对应ATCPL1/SHI4基因功能为编码的蛋白质在体外培养具有先天性磷酸酶活性, 突变体对脱落酸, 低温和氯化钠十分敏感^[26], 其在敏感型材料中的表达量在48 h或72 h高于耐盐型材料(图4), 可能是由于其对盐胁迫环境更为敏感而高量表达以响应盐胁迫。BnaC03g62830D基因在拟南芥中对应ATCESA8/LEW2基因功能为编码一个纤维素合酶家族的成员, 其参与次生细胞壁的生物合成^[28], 与对照组相比, 该基因在耐盐型材料中的盐处理后的表达量在48 h明显升高(图4), 之后降低, 推测是由于受到盐胁迫环境影响而高量表达以增强对盐胁迫环境的耐受力, 而敏感性材料则没有, 所以对逆境的耐受力较差。Long等^[37]通过RNA测序和PCR荧光定量方法在油菜中找到31个与盐胁迫相关的基因, 但与我们发现的基因位点不一致, 推测可能是由于取材时间和部位不同, 由基因表达的时间和组织特异性引起的。
4 结论盐胁迫环境下油菜幼苗的鲜重和干重均低于对照组, 表明盐胁迫对油菜生长有一定的影响。在染色体A02、A04和C03定位了6个盐胁迫环境下油菜幼苗干重鲜重的QTL, 贡献率超过10%的QTL有4个, 表明油菜耐盐性是由多个微效基因调控的数量性状。共找到8个候选基因, 对其中4个基因在极端表型材料中的qRT-PCR分析表明, BnaA02g14680D和BnaA02g14490D基因的表达由盐胁迫引起, BnaC03g62830D基因相对高量表达可能会增强植株对盐胁迫环境的耐受力。今后的研究应对筛选出的候选基因进行功能验证等深入分析, 为耐盐基因定位和盐胁迫相关机制的研究提供理论依据。
The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.


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

[1]	杨真, 王宝山. 中国盐渍土资源现状及改良利用对策. 山东农业科学, 2015, 47(4): 125-130 Yang Z, Wang B S. Present status of saline soil resources and countermeasures for improvement and utilization in China. Shand ong Agric Sci, 2015, 47(4): 125-130 (in Chinese with English abstract)[本文引用:1]
[2]	Qadir M, Ghafoor A, Murtaza G. Amelioration strategies for saline soils: a review. Land Degrad Devel, 2000, 11: 501-521[本文引用:1]
[3]	沈金雄, 傅廷栋. 我国油菜生产、改良与食用油供给安全. 中国农业科技导报, 2011, 13(1): 1-8 Shen J X, Fu T D. Rapeseed production, improvement and edible oil supply in China. J Agric Sci Technol, 2011, 13(1): 1-8 (in Chinese with English abstract)[本文引用:1]
[4]	陈宗金, 蔡士宾, 杨继书, 张巧凤, 吴纪中, 蒋彦婕, 颜伟, 吴小有. 主要农作物芽期耐盐性比较研究. 农业科学, 2012, 2(4): 59-65 Chen Z J, Cai S B, Yang J S, Zhang Q F, Wu J Z, Jiang Y J, Yan W, Wu X Y. Comparison of salinity tolerance among main crops at germination stage. Hans J Agric Sci, 2012, 2(4): 59-65 (in Chinese with English abstract)[本文引用:1]
[5]	王佳丽, 黄贤金, 钟太洋, 陈志刚. 盐碱地可持续利用研究综述. 地理学报, 2011, 66: 673-684 Wang J L, Huang X J, Zhong T Y, Chen Z G. Review on Saline-alkali land sustainable utilization research. Acta Geograph Sin, 2011, 66: 673-684 (in Chinese with English abstract)[本文引用:1]
[6]	易斌, 陈伟, 马朝芝, 傅廷栋, 涂金星. 甘蓝型油菜产量及相关性状的QTL分析. 作物学报, 2006, 32: 676-682 Yi B, Chen W, Ma C Z, Fu T D, Tu J X. Mapping of quantitative trait loci for yield and yield components in Brassica napus L. Acta Agron Sin, 2006, 32: 676-682 (in Chinese with English abstract)[本文引用:1]
[7]	Yang P, Shu C, Chen L, Xu J S, Wu J S, Liu K D. Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L. ). Theor Appl Genet, 2012, 125: 285-296[本文引用:1]
[8]	孙美玉. 甘蓝型油菜含油量QTLs定位及候选基因筛选. 中国农业科学院博士论文, 湖北武汉, 2012. pp 51-54 Sun M Y. Mapping of QTLs and Screening of Cand idate Genes for Oil Content in Brassica napus. PhD Dissertation of Chinese Academy of Agricultural Sciences, Wuhan, China, 2012. pp 51-54 (in Chinese with English abstract)[本文引用:1]
[9]	Yang M, Ding G D, Shi L, Feng J, Xu F S, Meng J L. Quantitative trait loci for root morphology in response to low phosphorus stress in Brassica napus. Theor Appl Genet, 2010, 121: 181-193[本文引用:1]
[10]	张凤启, 刘越英, 程晓辉, 童超波, 董彩华, 唐敏强, 黄军艳, 刘胜毅. 利用高密度SNP标记定位甘蓝型油菜株高QTL. 中国油料作物学报, 2014, 36: 695-700 Zhang F Q, Liu Y Y, Cheng X H, Tong C B, Dong C H, Tang M Q, Huang J Y, Liu S Y. QTL mapping of plant height using high density SNP markers in Brassica napus. Chin J Oil Crop Sci, 2014, 36: 695-700 (in Chinese with English abstract)[本文引用:1]
[11]	Yıldız M, Akçalı N, Terzi H. Proteomic and biochemical responses of canola (Brassica napus L. ) exposed to salinity stress and exogenous lipoic acid. J Plant Physiol, 2015, 179: 90-99[本文引用:1]
[12]	刘国红, 姜超强, 刘兆普, 梁明祥, 殷祥贞, 郑青松. 盐胁迫对油菜幼苗生长和光合特征的影响. 生态与农村环境学报, 2012, 28(2): 157-164 Liu G H, Jiang C Q, Liu Z P, Liang M X, Yin X Z, Zheng Q S. Effects of salt stress on growth and photosynthetic traits of Canola seedlings. J Ecol Rural Environ, 2012, 28(2): 157-164 (in Chinese with English abstract)[本文引用:1]
[13]	郑青松, 刘海燕, 隆小华, 刘兆普, 牛丹丹, 高影影. 盐胁迫对油菜幼苗离子吸收和分配的影响. 中国油料作物学报, 2010, 32: 65-70 Zheng Q S, Liu H Y, Long X H, Liu Z P, Niu D D, Gao Y Y. Effects of salt stress on ionic absorption and distribution of rapeseed seedlings. Chin J Oil Crop Sci, 2010, 32: 65-70 (in Chinese with English abstract)[本文引用:1]
[14]	Li Z, Mei S F, Mei Z, Liu X L, Fu T D, Zhou G S, Tu J X. Mapping of QTL associated with waterlogging tolerance and drought resistance during the seedling stage in oilseed rape (Brassica napus). Euphytica, 2014, 197: 341-353[本文引用:1]
[15]	荐红举, 肖阳, 李加纳, 马珍珍, 魏丽娟, 刘列钊. 利用SNP遗传图谱定位盐、旱胁迫下甘蓝型油菜种子发芽率的QTL. 作物学报, 2014, 40: 629-635 Jian H J, Xiao Y, Li J N, Ma Z Z, Wei L J, Liu L Z. QTL mapping for germination percentage under salinity and drought stresses in Brassica napus L. using a SNP genetic map. Acta Agron Sin, 2014, 40: 629-635 (in Chinese with English abstract)[本文引用:3]
[16]	Moursi Y S S. Genetic Mapping of QTL Controlling Salt Tolerance and Glucosinolates in Brassica napus and Brassica oleracea. PhD Dissertation of Georg-August-University, Germany, Göttingen, 2014. pp 73-78[本文引用:2]
[17]	Liu L Z, Qu C M, Wittkop B, Yi B, Xiao Y, He Y J, Snowdon R J, Li J N. High-Density SNP Map for Accurate Mapping of Seed Fibre QTL in Brassica napus L. PLoS One, 2013, 8: e83052. DOI: DOI:10.1371/journal.pone.0083052[本文引用:2]
[18]	Wang S, Basten C J, Zeng Z B. Windows QTL Cartographer. Ver. 2. 5 [computer program] Department of Statistics, North Carolina State University, Raleigh, NC, 2006. http://statgen.ncsu.edu/qtlcart/WQTLCart.htm[本文引用:1]
[19]	Land er E S, Botstein D. Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics, 1989, 121: 185-199[本文引用:1]
[20]	Mccouch S R, Cho Y G, Yano M, Paul E, Blinstrub M, Morishima H, Kinoshita T. Report on QTL nomenclature. Rice Genet Newslett, 1997, 14: 11-13[本文引用:1]
[21]	Chalhoub B, Denoeud F, Liu S Y, Isobe I A, Tang H B, Wang X Y, Chiquet J, Belcram H, Tong C B, Samans B, Corréa M, Silva C D, Just J, Falentin C, Koh S H, Clainche I L, Bernard M, Bento P, Noel B, Labadie K, Alberti A, Charles M, Arnaud D, Guo H, Daviaud C, Alamery S, Jabbari K, Zhao M X, Edger P P, Chelaifa H, Tack D, Lassalle G, Mestiri I, Schnel N, Paslier M L, Fan G Y, Renault V, Bayer P E, Golicz A A, Manoli S, Lee T H, Thi V H D, Chalabi S, Hu Q, Fan C C, Tollenaere R, Lu Y H, Battail C, Shen J X, Sidebottom C H. D, Wang X F, Canaguier A, Chauveau A, Bérard A, Deniot G, Guan M, Liu Z S, Sun F M, Lim Y P, Lyons E, Christopher D, Town, Ian B, Wang X W, Meng J L, Ma J X, Pires J C, Graham J, King, Brunel D, Delourme R, Renard M, Aury J M, Adams K L, Batley J, Snowdon R J, Tost J, Edwards D, Zhou Y M, Hua W, Sharpe A G, Paterson A H, Guan C Y, Wincker P. Early allopolyploid evolution in the post-neolithic Brassica napus oilseed genome. Science, 2014, 345: 950-953[本文引用:2]
[22]	Feng J L, Li J J, Gao Z X, Lu Y R, Yu J Y, Zheng Q, Yan S N, Zhang W J, He H, Ma L G, Zhu Z G. SKIP confers osmotic tolerance during salt stress by controlling alternative gene splicing in Arabidopsis. Mol Plant, 2015, 87: 1038-1052[本文引用:1]
[23]	Hong S W, Jon J H, Kwak J M, Nam H G. Identification of a receptor-like protein kinase gene rapidly induced by abscisic acid, dehydration, high salt, and cold treatments in Arabidopsis thaliana. Plant Physiol, 1997, 113: 1203-1212[本文引用:1]
[24]	Woei-Jiun G, Tuan-Hua D H. An abscisic acid-induced protein, HVA22, inhibits gibberellin-mediated programmed cell death in cereal aleurone cells. Plant Physiol, 2008, 148: 1182-1182[本文引用:1]
[25]	Elfving N, Davoinea C, Benllochb R, Blomberga J, Brännströma K, Müllerc D, Nilssond A, Ulfstedtd M, Ronned H, Wingsleb G, Nilssonb O, Björklunda S. The Arabidopsis thaliana Med25 mediator subunit integrates environmental cues to control plant development. Proc Natl Acad Sci USA, 2011, 108: 8245-8250[本文引用:1]
[26]	Ortega-Amaro M A, Rodriguez-Hernand ez A A, Rodriguez- Kessler M, Hernand ez-Lucero E, Rosales-Mendoza S, Ibañez- Salazar A, Delgado-Sanchez P, Jimenez-Bremont J F. Overexpression of AtGRDP2, a novel glycine-rich domain protein, accelerates plant growth and improves stress tolerance. Front Plant Sci, 2015, 5: 782-783[本文引用:1]
[27]	Seifert G J, Xue H, Acet T. The Arabidopsis thaliana Fasciclin like arabinogalactan protein 4 gene acts synergistically with abscisic acid signalling to control root growth. Ann Bot, 2014, 114: 1125-1133[本文引用:1]
[28]	Ohnuma T, Numata T, Osawa T, Mizuhara M, Lampela O, Juffer A H, Skriver K, Fukamizo T. A class V chitinase from Arabidopsis thaliana: gene responses, enzymatic properties, and crystallographic analysis. Planta, 2011, 234: 123-137[本文引用:1]
[29]	Glass M, Barkwill S, Unda F, Mansfield S D. Endo-β-1, 4- glucanases impact plant cell wall development by influencing cellulose crystallization. J Integr Plant Biol, 2015, 57: 396-410[本文引用:1]
[30]	Munns R, Tester M. Mechanisms of salinity tolerance. Annu Rev Plant Biol, 2008, 59: 651-681[本文引用:1]
[31]	Ashraf M, McNeilly T. Salinity tolerance in Brassica oilseeds. Crit Rev Plant Sci, 2004, 23: 157-174[本文引用:1]
[32]	Tunçtürk M, Tunçtürk R, Yildirim B, Çiftçi V. Changes of micronutrients, dry weight and plant development in canola (Brassica napus L. ) cultivars under salt stress. Afr J Biotechnol, 2011, 10: 3726-3730[本文引用:1]
[33]	Yeo A R. Molecular biology of salt tolerance in the context of whole-plant physiology. J Exp Bot, 1998, 49: 915-929[本文引用:1]
[34]	Flowers T J, Yeo A R. Breeding for salinity resistance in crop plants: where next?Aust J Plant Physiol, 1995, 22: 875-884[本文引用:1]
[35]	Osakabe Y, Mizuno S, Tanaka H, Maruyama K, Osakabe K, Todaka D, Fujita Y, Kobayashi M, Shinozaki K, Yamaguchi- Shinozaki K. Overproduction of the membrane-bound receptor- like protein kinase 1, RPK1, enhances abiotic stress tolerance in Arabidopsis. J Biol Chem, 2010, 285: 9190-9201[本文引用:1]
[36]	Shi C C, Feng C C, Yang M M, Li J L, Li X X, Zhao B C, Huang Z J, Ge R C. Overexpression of the receptor-like protein kinase genes AtRPK1 and OsRPK1 reduces the salt tolerance of Arabidopsis thaliana. Plant Sci, 2014, 217: 63-70[本文引用:1]
[37]	Long W H, Zou X L, Zhang X K. Transcriptome analysis of canola (Brassica napus) under salt stress at the germination stage. PLoS One, 2015, 10(2): e0116217[本文引用:1]