生态学报  2021, Vol. 41 Issue (18): 7322-7330

文章信息

杨安琪, 余海英, 黄化刚, 张锡洲, 李廷轩
YANG Anqi, YU Haiying, HUANG Huagang, ZHANG Xizhou, LI Tingxuan
籽粒镉低积累水稻地上部镉高积累遗传特性分析
Pivotal genetic mechanisms of the shoot cadmium accumulation in a low grain cadmium-accumulating rice
生态学报. 2021, 41(18): 7322-7330
Acta Ecologica Sinica. 2021, 41(18): 7322-7330
http://dx.doi.org/10.5846/stxb202006281668

文章历史

收稿日期: 2020-06-28
网络出版日期: 2021-06-15
籽粒镉低积累水稻地上部镉高积累遗传特性分析
杨安琪 , 余海英 , 黄化刚 , 张锡洲 , 李廷轩     
四川农业大学资源学院, 成都 611130
摘要: 农田土壤镉(Cd)污染日益严重,导致稻米Cd含量超标事件频繁出现,使粮食安全问题备受关注。因此,合理利用Cd污染农田、降低水稻籽粒Cd含量成为亟待解决的问题。籽粒Cd低积累水稻雅恢2816的地上部具有较强的Cd积累能力,研究旨在弄清其地上部Cd积累能力的遗传稳定性,进一步揭示控制该性状的遗传基础,为利用分子标记辅助选育地上部Cd富集能力强、籽粒Cd安全的水稻提供途径。以水稻雅恢2816和3个不同品种水稻分别组配获得的F1为研究对象,分析地上部Cd积累相关性状的杂种优势。进一步以优势组合C268A/雅恢2816构建F2作图群体,对地上部Cd积累相关性状进行QTL定位分析。结果表明:(1)F1地上部Cd积累相关性状杂种优势明显,遗传稳定性强。地上部Cd积累相关性状属数量性状,F2中/超亲分离现象明显。(2)在第4、6号染色体上共挖掘到1个控制水稻地上部生物量和3个控制地上部Cd积累量的QTL位点,分别为qSB-6、qSCdA-4、qSCdA-6-1和qSCdA-6-2,表型贡献率为10.6%—14.4%,且增效等位基因均来自雅恢2816。(3)地上部Cd积累量与地上部生物量、Cd含量,根、糙米的生物量、Cd积累量,根-地上部转移系数均呈极显著正相关,与地上部-籽粒转移系数呈极显著负相关,存在4个QTL集簇区Cl-4-1、Cl-6-1、Cl-6-2和Cl-6-3。(4)区间marker 04171-marker 04197控制着地上部生物量和Cd积累量,与控制糙米Cd含量的QTL不重叠。研究表明:籽粒Cd低积累水稻雅恢2816携带控制地上部Cd高积累的等位基因,可在世代间稳定遗传,QTL位点qSCdA-4、qSCdA-6-1、qSCdA-6-2是控制该性状的重要遗传基础,可为分子标记辅助选育地上部Cd高积累、籽粒Cd低积累水稻提供理论依据。
关键词: 水稻    Cd积累    杂种优势    数量性状    
Pivotal genetic mechanisms of the shoot cadmium accumulation in a low grain cadmium-accumulating rice
YANG Anqi , YU Haiying , HUANG Huagang , ZHANG Xizhou , LI Tingxuan     
College of Resources, Sichuan Agricultural University, Chengdu 611130, China
Abstract: Farmland soil is gradually polluted by cadmium (Cd), which leads to the excessive Cd concentration in grain. Food security has attracted much attention of people all over the world. Therefore, the utilization of Cd polluted farmland soil and reduction of the Cd concentration in grain have become an urgent issue to be resolved. Yahui2816, a low grain Cd-accumulating rice line with high shoot Cd-accumulation, is potential to achieve these goals simultaneously. The genetic stability of shoot Cd accumulation and the relevant genetic basis in Yahui 2816 was investigated in this study, providing an idea for molecular marker-assisted selection of low grain Cd-accumulating rice with high shoot Cd accumulation. Three types of F1 hybrids were constructed by rice restorer line Yahui 2816 and three male sterile lines (Lu98A, 5406A and C268A) to analyze the heterosis of shoot Cd accumulation. The F2 population was constructed by Yahui2816 and C268A. The F1 population (constructed by Yahui2816 and C268A) showed a distinct advantage and thus met the basic requirement of quantitative trait locus (QTL) analysis of phenotypic isolation. Furthermore, the F2 was used to analyze the QTLs of shoot Cd accumulation combined with the genetic linkage map. The results revealed that: (1) F1 hybrids showed a good heterosis phenomenon of shoot Cd accumulation, indicating Yahui2816 had strong genetic stability. Meanwhile, a mid-/super-parent segregation was observed in F2 population. (2) A QTL controlling shoot biomass (qSB-6) and three QTLs controlling shoot Cd accumulation (qSCdA-4, qSCdA-6-1 and qSCdA-6-2) were detected on chromosome 4 and 6, explaining 10.6%-14.4% of the phenotypic variance in the shoot Cd accumulation of the F2 population. (3) Shoot Cd accumulation showed significantly positive correlations with shoot biomass and Cd concentration, root and grain biomass and Cd accumulation, transfer factor of Cd to shoot, and a negative correlation with transfer factor of Cd to grain. In addition, there were 4 overlap regions (Cl-4-1, Cl-6-1, Cl-6-2 and Cl-6-3) controlling two or more traits at the same time. (4) Marker 04171-marker 04197 controlled shoot biomass and Cd accumulation simultaneously and did not overlap with the QTLs controlling grain Cd concentration. In conclusion, Yahui2816 carried an allele which can enhance the ability of shoot Cd accumulation, and its related traits could be stably inherited from generation to generation. QTLs (qSCdA-4, qSCdA-6-1 and qSCdA-6-2) controlling shoot Cd accumulation were beneficial for molecular marker-assisted breeding of rice with high shoot Cd accumulation and low grain Cd accumulation.
Key Words: Oryza sativa L.    Cd accumulation    heterosis    quantitative trait locus    

镉(Cd)是一种毒性极强的重金属元素, 其生物有效性高、迁移性强[1]。随着工业化和城镇化进程的加快, 农田土壤Cd污染日益严重, 且易随食物链进入人体、危害人类健康[2-3]。水稻是人类膳食Cd的主要来源, 其平均Cd含量显著高于其他谷类食物[4]。因此, 安全利用Cd污染农田土壤、降低糙米Cd含量成为亟待解决的粮食安全问题。选育籽粒Cd低积累水稻能保障Cd污染农田的安全生产, 同时一些生物量大、地上部Cd积累能力较强的水稻品种表现出较好的Cd修复潜力[5-6], 因而选育地上部Cd高积累、籽粒Cd低积累的水稻品种利于实现Cd污染农田的修复和安全生产。明晰水稻Cd积累的遗传机制是选育籽粒Cd低积累或地上部Cd高积累水稻的重要前提。不同水稻品种Cd积累存在明显的基因型差异, 利于寻找相关功能基因[7-8]。研究发现, 水稻Cd积累相关性状属数量性状, 遗传机制复杂, 同时受环境因素影响, 因此不能仅通过表型对基因型进行选择[9]。QTL定位通过建立数量性状表型值与DNA分子标记间的关系, 可以确定各个QTL位点在染色体上的位置、效应及其相关作用[10-11], 使复杂数量性状的遗传改良和分子操纵成为可能, 具有较大的应用前景[12-13]。近年来, 利用Cd积累能力差异显著的水稻品种构建的F2、RIL、CSSLs等群体进行QTL定位, 在水稻12条染色体上发现了控制水稻籽粒和地上部Cd积累相关的QTLs[14-16]。Xue等[17]利用DH群体, 在第7号染色体上得到控制地上部Cd含量的QTL qCDS7。Ueno等[18]和Tezuka等[19]利用不同亲本构建的F2群体, 也在第7号染色体上得到控制地上部Cd积累量的QTL, 其中QTL qCdT7包含的基因OsHMA3发生碱基突变, 以该基因作为分子标记, 筛选出了Cd污染农田修复效果显著的水稻品种Akita 110[19-22]。除第7号染色体外, Ueno等[18]和Yan等[23]分别还在第2、5、11和10号染色体上定位到控制水稻地上部Cd积累量的QTL。此外, 还有研究发现水稻Cd积累相关QTL如qCd-2、qCd-7、GCC7、qCd1-3、CAL1等影响籽粒、叶片、根等部位Cd含量[14-15, 24-25]。可见, 通过QTL定位挖掘水稻Cd积累的关键基因有助于阐明水稻Cd积累的遗传机制, 为利用分子标记辅助选择育种奠定基础。相同土壤环境条件下, 相较于其他籽粒Cd低积累水稻品种, 雅恢2816具有更强的地上部Cd积累能力, 可用于Cd污染农田边修复边生产。前期已对其籽粒Cd低积累特性进行了探讨[26-28], 但其地上部Cd积累的遗传稳定性和遗传机制尚不清楚。利用雅恢2816杂交后代F1和F2, 分析其地上部Cd积累相关性状的杂种优势, 挖掘控制地上部Cd积累相关性状的QTL, 明确其地上部Cd积累机制, 为分子标记辅助选择籽粒Cd低积累且地上部Cd高积累材料提供理论依据。

1 材料与方法 1.1 供试材料

水稻雅恢2816(R)、泸98A(S1)、5406A(S2)、C268A(S3)、泸98A/雅恢2816(H1)、5406A/雅恢2816(H2)、C268A/雅恢2816(H3)、F1群体(包含雅恢2816与3个不育系材料(泸98A、5406A、C268A)分别杂交得到的第一代群体H1、H2和H3)、F2群体(共120个株系, 为杂交组合H3的自交后代)。由四川农业大学农学院提供。雅恢2816是四川省常用杂交生产的籼稻品种, 具有地上部Cd高积累和籽粒Cd低积累特性;泸98A(S1)、5406A(S2)和C268A(S3)均为不育系水稻材料, 其地上部Cd含量与积累量在盆栽和大田试验中均显著低于雅恢2816[26, 28]

1.2 试验设计与处理

第一期试验共设21个试验小区(1 m×2 m), 种植亲本及F1, 每个材料设3次重复。第二期试验共设6个试验小区和1个试验大区(1 m×4 m), 分别种植亲本(雅恢2816、C268A)和F2。各小区随机区组排列, 各区间设20 cm缓冲带。每小区种植水稻32穴, 每大区种植水稻120穴, 每穴1株, 株行距为20 cm×25 cm。三叶期栽插, 采用旱地育秧, 于当年4月播种, 5月移栽, 按当地习惯进行水肥管理和病虫害防治。

试验区位于四川省成都平原某市的Cd污染农田, 属中纬度亚热带湿润气候区, 海拔507 m, 年均温15.9—16.1 ℃, 年均降水量900—1000 mm。试验土壤基本理化性质为pH 6.60、有机质34.61 g/kg、全氮2.43 g/kg、碱解氮100.31 mg/kg、有效磷19.03 mg/kg、速效钾69.25 mg/kg、全Cd 1.92 mg/kg、有效Cd 0.63 mg/kg。

1.3 样品采集与制备

于水稻抽穗期, 分单株采集雅恢2816、C268A和F2幼叶鲜样于2mL灭菌灭酶冻存管中, 经液氮冷冻后储存于-75℃超低温冰箱中, 用于DNA提取。于水稻成熟期采样, 将样品用自来水冲洗干净, 根部用20 mmol/L Na2-EDTA溶液浸泡15 min去除根表附着的Cd2+, 再用去离子水洗净, 最后用吸水纸擦干。将植株分为根、茎、叶、穗, 根、茎、叶于105 ℃下杀青30 min后, 75 ℃烘干至恒重, 称重后粉碎备用, 穗部经风干称重后, 脱粒制成糙米备用。

1.4 测定项目及方法

土壤基本理化性质采用常规分析方法[29];土壤全Cd含量采用HNO3-HClO4-HF(5∶1∶1, 体积分数)消化, 有效态Cd含量采用DTPA提取, 植株Cd含量采用HNO3-HClO4(5∶1, 体积分数)消化, 火焰原子吸收分光光度计(AA400P, Analytikjena, Germany)测定。

1.5 基因型分析

采用Takara试剂盒提取叶片DNA。利用简化基因组测序技术(Genotyping by sequencing)进行基因型鉴定。内切酶组合为EcorI/NIaIII(New England Biolabs, Ipswich, MA)。

1.6 QTL定位

利用已构建的遗传连锁图谱和表型数据[27], 采用WinQTL Cartographer v2.5程序包的复合区间作图法(Composite Interval Mapping, CIM)进行QTL定位分析。选择0.5 cM的扫描区间, 在整个遗传图谱上逐步扫描寻找QTL。同时, 选择距离待检测区间至少10 cM, 与性状相关性显著性最强的前15个标记作为协变量, 校正背景QTL效应对目标区间QTL分析的影响。以LOD值2.5作为阈值来判断QTL是否存在, 若标记区间中LOD>2.5, 则认为该区间LOD最高处所对应的位点为一个QTL, 并估计各个QTL对表型变异的贡献率。QTL的命名基本按照McCouch等[30]的方法并稍作改动:q+目标性状英文缩写名称(首字母大写)+“-” +连锁群的序号+“-”+同一连锁群上QTL数, 本文中QTL名称均用斜体表示。qSB指控制指地上部生物量的QTL;qSCdA指控制指地上部Cd积累量的QTL。

1.7 数据处理

地上部Cd积累量=茎生物量×茎Cd含量+叶生物量×叶Cd含量, 超亲优势=(F1表型值-高值亲本表型值)/高值亲本表型值×100%, 中亲优势=(F1表型值-双亲表型平均值)/双亲平均值×100%。

统计分析在SPSS 24.0中进行, 多重比较选择LSD法;图表制作采用Origin 8.0和Excel 2013。

2 结果与分析 2.1 F1地上部Cd积累相关性状的杂种优势

亲本及其F1地上部生物量、Cd含量和Cd积累量差异显著(表 1, 表 2)。总体而言, F1地上部生物量显著高于雅恢2816, 除H2, 杂交组合H1和H3与其不育系母本无显著差别。H1、H2和H3地上部生物量分别为雅恢2816的1.43、1.35和1.43倍。所有杂交组合地上部Cd含量显著低于雅恢2816, 显著高于不育系母本。所有杂交组合地上部Cd积累量显著高于各不育系母本, 与雅恢2816无显著差别。H1、H2和H3地上部Cd积累量分别为雅恢2816的1.17、1.16和1.16倍。F1地上部Cd积累相关性状表现出较好的中/超亲优势, Cd积累相关性状杂种优势受遗传效应的影响, C268A/雅恢2816后代表现较好, 可利用其F2进行遗传特性分析。

表 1 亲本及其F1地上部Cd积累相关性状 Table 1 Comparison of differences of traits related to shoot Cd accumulation between rice parents and their F1 hybrids
水稻材料
Rice lines
生物量
Biomass/(g/株)
Cd含量
Cd concentration/(mg/kg)
Cd积累量
Cd accumulation/(μg/株)
水稻材料
Rice lines
生物量
Biomass/(g/株)
Cd含量
Cd concentration/(mg/kg)
Cd积累量
Cd accumulation/(μg/株)
R 27.51±2.16c 8.88±0.22a 245.28±24.58ab H1 39.42±1.71ab 7.32±0.28bc 288.57±14.63a
S1 40.57±0.58a 5.02±0.35d 204.01±16.74b H2 37.22±0.49ab 7.64±0.17b 284.07±3.38a
S2 29.41±2.35c 6.77±0.11c 198.61±12.48bc H3 39.41±2.67ab 7.23±0.1bc 285.03±20.83a
S3 35.08±1.18b 4.29±0.25e 150.82±13.38c
字母代表不同水稻材料,即R: 雅恢2816, S1:泸98A, S2:5406A, S3:C268A, H1:泸98A/雅恢2816杂交第一代,H2:5406A/雅恢2816; 杂交第一代,H3:C268A/雅恢2816杂交第一代; 不同小写字母表示不同水稻间差异显著(P < 0.05)

表 2 F1地上部Cd积累相关性状的中/超亲优势 Table 2 Mid-/Super-parent heterosis of traits related to shoot Cd accumulation in F1 hybrids
杂交组合
Hybrids
中亲优势Mid-parent heterosis/% 超亲优势Super-parent heterosis/%
生物量
Biomass
Cd含量
Cd concentration
Cd积累量
Cd accumulation
生物量
Biomass
Cd含量
Cd concentration
Cd积累量
Cd accumulation
H1 15.83 18.67 28.46 -2.81 -17.56 17.65
H2 30.76 8.37 27.99 26.53 -14.04 15.82
H3 25.93 24.53 43.92 12.34 -18.64 16.21
2.2 F2地上部Cd积累相关性状的遗传特性分析

雅恢2816地上部Cd含量、Cd积累量显著高于C268A, 生物量无显著差异(表 3)。雅恢2816地上部Cd含量和Cd积累量分别为C268A的2.17和2.03倍。F2地上部生物量、Cd含量和Cd积累量表现为连续分布和中/超亲分离, 变异系数分别为51.14%、15.18%和57.79%, 其中Cd积累量变异较大。说明地上部Cd积累相关性状均为数量性状, 适合进行QTL定位。

表 3 亲本与F2地上部Cd积累相关性状的表型变异 Table 3 Phenotypic variation of traits related to shoot Cd accumulation in rice parents and F2 population
性状Trait 亲本Parents F2群体F2 population
雅恢2816 C268A 最大值
Maximum
最小值
Minimum
标准差
Standard deviation
平均值
Avgerage value
变异系数
Coefficient variation/%
生物量Biomass/(g/株) 40.46 43.18 93.70 7.56 19.06 37.27 51.14%
Cd含量Cd concentration/(mg/kg) 9.19* 4.23 11.37 5.93 1.26 8.28 15.18%
Cd积累量Cd accumulation/(μg/株) 371.77* 182.77 1007.00 57.95 179.55 310.69 57.79%
*表示两亲本间差异显著(P < 0.05)

对F2地上部Cd积累量与各性状间进行相关性分析(表 4), 除根和糙米Cd含量外, 地上部Cd积累量与地上部生物量、Cd含量, 根、糙米的生物量、Cd积累量和根-地上部转移系数呈极显著正相关, 与地上部-籽粒转移系数呈极显著负相关。说明水稻地上部Cd积累量与地上部生物量、Cd含量, 根、糙米的生物量、Cd积累量和根-地上部转移系数、变化趋势相同, 而与根-地上部转移系数变化趋势相反。

表 4 F2地上部Cd积累量与各器官Cd积累相关性状之间的相关性分析 Table 4 Correlation analysis between shoot Cd accumulation and traits related to organs Cd accumulation in F2 population
性状
Trait
地上部Cd积累量
Shoot Cd accumulation/(μg/株)
性状
Trait
地上部Cd积累量
Shoot Cd accumulation/(μg/株)
地上部生物量Shoot Biomass/(g/株) 0.953** 糙米生物量Grain Biomass/(g/株) 0.263**
地上部Cd含量Shoot Cd concentration/(mg/kg) 0.356** 糙米Cd含量Grain Cd concentration/(mg/kg) 0.067
根生物量Root Biomass/(g/株) 0.439** 糙米Cd积累量Grain Cd accumulation/(μg/株) 0.279**
根Cd含量Root Cd concentration/(mg/kg) 0.164 根-地上部转移系数Transfer factor of Cd to shoot 0.458**
根Cd积累量Root Cd accumulation/(μg/株) 0.458** 地上部-籽粒转移系数Transfer factor of Cd to grain -0.422**
*表示各性状相关性极显著(P < 0.01), *表示各性状相关性显著(P < 0.05)
2.3 F2地上部Cd积累相关性状QTL定位

利用C268A与雅恢2816构建的F2作图群体, 共定位到4个与地上部生物量和Cd积累量紧密连锁的主效QTL位点(表 5, 图 1)。其中, 地上部生物量紧密连锁的QTL qSB-6位于第6号染色体, 物理区间大小为4.28cM, 表型贡献率为14.4%。3个控制地上部Cd积累量的QTL qSCdA-4、qSCdA-6-1和qSCdA-6-2分别位于第4、6和6号染色体, 其物理区间大小分别是9.64、13.5和2.43cM, 表型贡献率分别为11.7%、14.4%和10.6%。4个QTL的加性效应均为负, 表明控制增加地上部生物量及Cd积累量的基因均来自父本雅恢2816。在第6号染色体上, 区间marker 04171-marker 04197同时控制地上部生物量与Cd积累量, 其物理区间大小为3.28cM。

表 5 F2地上部Cd积累相关性状QTL定位 Table 5 QTL mapping of traits related to shoot Cd accumulation in F2 population
性状
Trait
QTL 染色体编号
Chromosome
No.
标记区间
Maker interval
阈值
LOD
物理区间
Physical >interval
表型贡献率
R2 /%
加性效应
Additive effect
地上部生物量
Shoot Biomass
qSB-6 6 marker 04171-marker 04197 4.06 17.24-21.52 14.4 -8.94
地上部Cd积累量 qSCdA-4 4 marker 03334-marker 03367 3.25 92.98-102.62 11.7 -94.78
Shoot Cd Accumulation qSCdA-6-1 6 marker 04124-marker 04200 4.06 9.77-23.27 14.4 -91.97
qSCdA-6-2 6 marker 04225-marker 04239 2.92 26.26-28.69 10.6. -81.18

图 1 F2地上部Cd积累相关性状QTL的染色体定位 Fig. 1 Chromosome location of QTL related to shoot Cd accumulation in F2 population QTL qT-r-s-6-2、qGB-4和qRCdA-6来源于前期研究结果[27, 31]
3 讨论 3.1 水稻地上部Cd积累的杂种优势及遗传机理

Cd胁迫下, 植物的正常生长发育会受到影响, 首先体现在根、茎、叶生物量的变化[32]。本研究中, 杂交组合泸98A/雅恢2816、5406A/雅恢2816和C268A/雅恢2816的地上部生物量显著高于父本雅恢2816, 表现出中/超亲优势, 表明利用杂交育种, 可以明显提高水稻的生物量。同时, 各杂交组合地上部Cd含量和Cd积累量显著高于不育系母本, 表现出中/超亲优势。前期发现, 在Cd污染大田和不同Cd处理的盆栽试验中, 各杂交组合籽粒Cd含量显著低于双亲, 表现出明显的负向超亲优势[31, 33]。研究表明, 杂交水稻对Cd的积累主要受双亲遗传背景影响[34]。雅恢2816具有籽粒Cd低积累、地上部Cd高积累特性, 可稳定遗传给后代。因此, 以上杂交组合可应用于Cd污染农田的边生产边修复。具有杂种优势的性状一般属数量性状, 受多基因控制, 亲本材料的遗传背景差异可能导致杂交组合Cd积累量存在较大差异[34-35], 从分子层面深入挖掘水稻Cd积累的遗传机制可提高杂种优势利用效率。本研究利用地上部Cd积累特性差异显著的亲本材料, 在第4和6号染色体上得到控制地上部Cd积累量的QTL qSCdA-4、qSCdA-6-1和qSCdA-6-2, 与前人结果比较, 是新的控制水稻地上部Cd积累量的主效QTL。

3.2 水稻Cd积累相关性状遗传关系分析

了解水稻地上部与籽粒Cd积累的遗传关系对于实现Cd污染农田边生产边修复目标至关重要。Wang等[36]利用RIL群体, 在7号染色体上定位到同时控制水稻地上部Cd含量和糙米Cd含量的QTL qCd7.1。Ishikawa等[37]利用BILs群体, 在第2号染色体上控制籽粒Cd含量的QTL qGCd2与控制地上部Cd含量的QTL qSCd2紧密连锁, 在第7号染色体上控制籽粒Cd含量QTL qGCd7和地上部Cd含量的QTL qSCd7完全重叠。Hu等[38]与Abe等[39]在3号染色体上分别定位到控制糙米Cd含量的QTL qCCBR3和控制地上部Cd含量的QTL qlGCd3, 二者存在重叠区域。Ueno等[18]和Liu等[15]利用不同的作图群体, 以水稻地上部Cd积累量与糙米Cd含量为指标进行QTL定位, 在第2号染色体上, 两个QTL存在重叠部分。本研究中控制水稻地上部Cd积累量的QTL与前期得到的4个控制糙米Cd含量的QTL位点位于不同的染色体[27], 与Wang等[36]和Ishikawa等[37]的研究结果不同, 且与其他学者比较, 也并无重叠区域, 与Hu等[38]和Abe等[39]、Ueno等[18]和Liu等[15]结果不同。这可能与不同研究的双亲遗传差异和群体类型不同有关[40]。本研究中, 控制糙米Cd含量和地上部Cd积累量的QTL增效等位基因来源不同, 且定位于不同染色体, 表明源于雅恢2816的等位基因不会同时增加地上部Cd积累量和糙米Cd含量。

水稻地上部Cd积累量与各器官生物量、Cd积累量、转移系数之间存在明显的相关性, 且在4号和6号染色体出现QTL集簇区Cl4-1、Cl6-1、Cl6-2和Cl6-3。在4号染色体上, 集簇区Cl4-1发生qGB-4与qSCdA-4重叠, 且加性效应一致[31], 表明Cl4-1能增加地上部Cd积累量并提高稻米产量。经基因预测, 在集簇区Cl4-1内含有基因Os04g0613000(OsZIP3)。研究发现, OsZIP3在水稻节点表达较高, 负责扩大维管束的木质部中Zn的卸载, 从而调控Zn在地上部的分配[41]。Cd和Zn具有相似的化学性质, 一些Zn转运蛋白同时负责Cd的转运[42], 因而推测OsZIP3可能参与了Cd的分配, 进一步研究Os04g0613000的功能, 有助于解释雅恢2816地上部Cd积累机制。亲本雅恢2816具有较高的根-地上部转移系数, 且地上部-籽粒转移系数较低, 使其具有籽粒Cd低积累、地上部Cd高积累特性。相关研究表明, OsHMA2和OsHMA3在Cd从根向地上部的转运过程中起到重要作用, OsHMA2参与Cd向木质部的装载[43], 而OsHMA3能将Cd区隔化在根细胞的液泡中[20, 44]。但在本研究中同时控制根-地上部转移和地上部Cd积累量的集簇区Cl6-1不包含OsHMA2和OsHMA3, 有望在集簇区Cl6-1挖掘到类似功能的新基因。集簇区Cl6-2发生qSCdA-6-1、qSB-6、qRCdA-6重叠, Cl6-3发生qSCdA-6-1、qRCdA-6重叠[31], F2地上部Cd积累量与地上部生物量、根Cd积累量呈极显著正相关, 可见增加地上部生物量和根系Cd积累量, 能有效提高水稻对Cd的提取效率。控制数量性状的QTL在染色体上成簇分布的现象可见于绝大多数研究中, 这可能是同一基因控制不同性状的表达, 也可能是控制多个性状表达的基因紧密连锁的结果[39, 45], 后续可针对重叠区段进行深入研究, 以期挖掘控制地上部Cd积累的关键基因。

4 结论

籽粒Cd低积累亲本材料雅恢2816具有地上部Cd高积累特性, 且能稳定遗传。地上部Cd积累相关性状由多基因控制, F2代中/超亲优势明显。QTL qSB-6控制地上部生物量, QTL qSCdA-4、qSCdA-6-1和qSCdA-6-2控制地上部Cd积累量。区间marker 04171-marker 04197同时关联着水稻地上部生物量(qSB-6)与地上部Cd积累量(qSCdA-6-1), 且与控制籽粒Cd含量的QTL未发生重叠, 该区间为后期利用分子标记辅助育种同时实现籽粒Cd低积累、地上部Cd高积累提供可能。

参考文献
[1]
Touceda-González M, Brader G, Antonielli L, Ravindran V B, Waldner G, Friesl-Hanl W, Corretto E, Campisano A, Pancher M, Sessitsch A. Combined amendment of immobilizers and the plant growth-promoting strain Burkholderia phytofirmans PsJN favours plant growth and reduces heavy metal uptake. Soil Biology & Biochemistry, 2015, 91: 140-150.
[2]
Du B Y, Zhou J, Lu B X, Zhang C, Li D M, Zhou J, Jiao S J, Zhao K Q, Zhang H H. Environmental and human health risks from cadmium exposure near an active lead-zinc mine and a copper smelter, China. Science of the Total Environment, 2020, 720: 137585. DOI:10.1016/j.scitotenv.2020.137585
[3]
Tang L, Deng S H, Tan D, Long J M, Lei M. Heavy metal distribution, translocation, and human health risk assessment in the soil-rice system around Dongting Lake area, China. Environmental Science and Pollution Research, 2019, 26(17): 17655-17665. DOI:10.1007/s11356-019-05134-w
[4]
Song Y, Wang Y, Mao W F, Sui H X, Yong L, Yang D J, Jiang D G, Zhang L, Gong Y Y. Dietary cadmium exposure assessment among the Chinese population. PLoS One, 2017, 12(5): e0177978. DOI:10.1371/journal.pone.0177978
[5]
Hu Y A, Cheng H F, Tao S. The challenges and solutions for cadmium-contaminated rice in China: a critical review. Environment International, 2016, 92-93: 515-532. DOI:10.1016/j.envint.2016.04.042
[6]
Murakami M, AE N, Ishikawa S, Ibaraki T, Ito M. Phytoextraction by a high-Cd-accumulating rice: reduction of Cd content of soybean seeds. Environmental Science & Technology, 2008, 42(16): 6167-6172.
[7]
Zhao J L, Yang W, Zhang S H, Yang T F, Liu Q, Dong J F, Fu H, Mao X X, Liu B. Genome-wide association study and candidate gene analysis of rice cadmium accumulation in grain in a diverse rice collection. Rice, 2018, 11(1): 61. DOI:10.1186/s12284-018-0254-x
[8]
肖清铁, 王经源, 郑新宇, 戎红, 张国君, 王良华, 谢惠玲, 李艺, 陈珊, 林瑞余, 林文雄. 水稻根系响应镉胁迫的蛋白质差异表达. 生态学报, 2015, 35(24): 8276-8283.
[9]
Rao Y C, Li Y Y, Qian Q. Recent progress on molecular breeding of rice in China. Plant Cell Reports, 2014, 33(4): 551-564. DOI:10.1007/s00299-013-1551-x
[10]
Chen J G, Zou W L, Meng L J, Fan X R, Xu G H, Ye G Y. Advances in the uptake and transport mechanisms and QTLs mapping of cadmium in rice. International Journal of Molecular Sciences, 2019, 20(14): 3417. DOI:10.3390/ijms20143417
[11]
Sehgal D, Singh R, Rajpal V R. Quantitative trait loci mapping in plants: concepts and approaches//Rajpal V R, Rao S R, Raina S N, eds. Molecular Breeding for Sustainable Crop Improvement. Cham: Springer, 2016: 31-59.
[12]
邹文莉, 冷语佳, 李娇龙, 孟丽君, 贺浩华, 陈景光, 叶国友. 水稻镉积累的吸收与转运机制及镉胁迫相关QTL定位的研究. 分子植物育种, 2018, 16(24): 8128-8141.
[13]
Yoneyama T, Ishikawa S, Fujimaki S. Route and regulation of zinc, cadmium, and iron transport in rice plants (Oryza sativa L.) during vegetative growth and grain filling: metal transporters, metal speciation, grain Cd reduction and Zn and Fe biofortification. International Journal of Molecular Sciences, 2015, 16(8): 19111-19129. DOI:10.3390/ijms160819111
[14]
Luo J S, Huang J, Zeng D L, Peng J S, Zhang G B, Ma H L, Guan Y, Yi H Y, Fu Y L, Han B, Lin H X, Qian Q, Gong J M. A defensin-like protein drives cadmium efflux and allocation in rice. Nature Communications, 2018, 9(1): 645. DOI:10.1038/s41467-018-03088-0
[15]
Liu W Q, Pan X W, Li Y C, Duan Y H, Min J, Liu S X, Liu L C, Sheng X N, Li X X. Identification of QTLs and validation of qCd-2associated with grain cadmium concentrations in rice. Rice Science, 2019, 26(1): 42-49. DOI:10.1016/j.rsci.2018.12.003
[16]
Huang Y, Sun C X, Min J, Chen Y L, Tong C, Bao J S. Association mapping of quantitative trait loci for mineral element contents in whole grain rice (Oryza sativa L.). Journal of Agricultural and Food Chemistry, 2015, 63(50): 10885-10892. DOI:10.1021/acs.jafc.5b04932
[17]
Xue D W, Chen M C, Zhang G P. Mapping of QTLs associated with cadmium tolerance and accumulation during seedling stage in rice (Oryza sativa L.). Euphytica, 2009, 165(3): 587-596. DOI:10.1007/s10681-008-9785-3
[18]
Ueno D, Kono I, Yokosho K, Ando T, Yano M, Ma J F. A major quantitative trait locus controlling cadmium translocation in rice (Oryza sativa). New Phytologist, 2009, 182(3): 644-653. DOI:10.1111/j.1469-8137.2009.02784.x
[19]
Tezuka K, Miyadate H, Katou K, Kodama I, Matsumoto S, Kawamoto T, Masaki S, Satoh H, Yamaguchi M, Sakurai K, Takahashi H, Satoh-Nagasawa N, Watanabe A, Fujimura T, Akagi H. A single recessive gene controls cadmium translocation in the cadmium hyperaccumulating rice cultivar Cho-Ko-Koku. Theoretical and Applied Genetics, 2010, 120(6): 1175-1182. DOI:10.1007/s00122-009-1244-6
[20]
Miyadate H, Adachi S, Hiraizumi A, Tezuka K, Nakazawa N, Kawamoto T, Katou K, Kodama I, Sakurai K, Takahashi H, Satoh-Nagasawa N, Watanabe A, Fujimura T, Akagi H. OsHMA3, a P1B-type of ATPase affects root-to-shoot cadmium translocation in rice by mediating efflux into vacuoles. New Phytologist, 2011, 189(1): 190-199. DOI:10.1111/j.1469-8137.2010.03459.x
[21]
Murakami M, Nakagawa F, AE N, Ito M, Arao T. Phytoextraction by rice capable of accumulating Cd at high levels: reduction of Cd content of rice grain. Environmental Science & Technology, 2009, 43(15): 5878-5883.
[22]
Takahashi R, Ito M, Katou K, Sato K, Nakagawa S, Tezuka K, Akagi H, Kawamoto T. Breeding and characterization of the rice (Oryza sativa L.) line "Akita 110" for cadmium phytoremediation. Soil Science and Plant Nutrition, 2016, 62(4): 373-378. DOI:10.1080/00380768.2016.1150791
[23]
Yan Y F, Lestari P, Lee K J, Kim M Y, Lee S H, Lee B W. Identification of quantitative trait loci for cadmium accumulation and distribution in rice (Oryza sativa). Genome, 2013, 56(4): 227-232. DOI:10.1139/gen-2012-0106
[24]
Liu C L, Gao Z Y, Shang L G, Yang C H, Ruan B P, Zeng D L, Guo L B, Zhao F J, Huang C F, Qian Q. Natural variation in the promoter of OsHMA3 contributes to differential grain cadmium accumulation between Indica and Japonica rice. Journal of Integrative Plant Biology, 2020, 62(3): 314-329. DOI:10.1111/jipb.12794
[25]
Pan X W, Li Y C, Liu W Q, Liu S X, Min J, Xiong H B, Dong Z, Duan Y H, Yu Y Y, Li X X. QTL mapping and candidate gene analysis of cadmium accumulation in polished rice by genome-wide association study. Scientific Reports, 2020, 10(1): 11791. DOI:10.1038/s41598-020-68742-4
[26]
李堃, 余海英, 黄富, 陈光登, 张路. 镉低积累水稻亲本及其杂交组合镉积累特征分析. 农业环境科学学报, 2016, 35(8): 1444-1450.
[27]
Guo J Y, Li K, Zhang X Z, Huang H G, Huang F, Zhang L, Wang Y D, Li T X, Yu H Y. Genetic properties of cadmium translocation from straw to brown rice in low-grain cadmium rice (Oryza sativa L.) line. Ecotoxicology and Environmental Safety, 2019, 182: 109422. DOI:10.1016/j.ecoenv.2019.109422
[28]
Zhang H J, Zhang X Z, Li T X, Huang F. Variation of cadmium uptake, translocation among rice lines and detecting for potential cadmium-safe cultivars. Environmental Earth Sciences, 2014, 71(1): 277-286. DOI:10.1007/s12665-013-2431-y
[29]
鲁如坤. 土壤农业化学分析方法. 北京: 中国农业科技出版社, 1999.
[30]
McCouch S R, Cho Y G, Yano M, Paul E. Report on QTL nomenclature. Rice Genet Newsl, 1997, 14: 11-13.
[31]
李堃. 镉低积累水稻遗传分析及镉积累相关性状QTL定位[D]. 雅安: 四川农业大学, 2017.
[32]
顾翠花, 王懿祥, 白尚斌, 吴建强, 杨一. 四种园林植物对土壤镉污染的耐受性. 生态学报, 2015, 35(8): 2536-2544.
[33]
Li K, Yu H Y, Li T X, Chen G D, Huang F. Cadmium accumulation characteristics of low-cadmium rice (Oryza sativa L.) line and F1 hybrids grown in cadmium-contaminated soils. Environmental Science and Pollution Research, 2017, 24(21): 17566-17576. DOI:10.1007/s11356-017-9350-5
[34]
Sun L, Xu X X, Jiang Y R, Zhu Q H, Yang F, Zhou J Q, Yang Y Z, Huang Z Y, Li A H, Chen L H, Tang W B, Zhang G Y, Wang J R, Xiao G Y, Huang D Y, Chen C Y. Genetic diversity, rather than cultivar type, determines relative grain Cd accumulation in hybrid rice. Frontiers in Plant Science, 2016, 7: 1407.
[35]
Yu Y H, Zhu M M, Cui Y, Liu Y, Li Z Y, Jiang N, Xu Z J, Xu Q, Sui G M. Genome sequence and QTL analyses using backcross recombinant inbred lines (BILs) and BILF1 lines uncover multiple heterosis-related loci. International Journal of Molecular Sciences, 2020, 21(3): 780. DOI:10.3390/ijms21030780
[36]
Wang C C, Tang Z, Zhuang J Y, Tang Z, Huang X Y, Zhao F J. Genetic mapping of ionomic quantitative trait loci in rice grain and straw reveals OsMOT1;1 as the putative causal gene for a molybdenum QTL qMo8. Molecular Genetics and Genomics, 2020, 295(2): 391-407. DOI:10.1007/s00438-019-01632-1
[37]
Ishikawa S, Abe T, Kuramata M, Yamaguchi M, Ando T, Yamamoto T, Yano M. A major quantitative trait locus for increasing cadmium-specific concentration in rice grain is located on the short arm of chromosome 7. Journal of Experimental Botany, 2010, 61(3): 923-934. DOI:10.1093/jxb/erp360
[38]
Hu D W, Sheng Z H, Li Q L, Chen W, Wei X J, Xie L H, Jiao G A, Shao G N, Wang J L, Tang S Q, Hu P S. Identification of QTLs associated with cadmium concentration in rice grains. Journal of Integrative Agriculture, 2018, 17(7): 1563-1573. DOI:10.1016/S2095-3119(17)61847-1
[39]
Abe T, Nonoue Y, Ono N, Omoteno M, Kuramata M, Fukuoka S, Yamamoto T, Yano M, Ishikawa S. Detection of QTLs to reduce cadmium content in rice grains using LAC23/Koshihikari chromosome segment substitution lines. Breeding Science, 2013, 63(3): 284-291. DOI:10.1270/jsbbs.63.284
[40]
林晗, 徐江民, 胡瑚倩, 郑安, 徐婉璐, 漏平, 王跃星, 曾大力, 饶玉春. 水稻耐金属离子胁迫的QTL分析. 中国水稻科学, 2018, 32(1): 23-34.
[41]
Sasaki A, Yamaji N, Mitani-Ueno N, Kashino M, Ma J F. A node-localized transporter OsZIP3 is responsible for the preferential distribution of Zn to developing tissues in rice. The Plant Journal, 2015, 84(2): 374-384. DOI:10.1111/tpj.13005
[42]
Tan L T, Zhu Y X, Fan T, Peng C, Wang J R, Sun L, Chen C Y. OsZIP7 functions in xylem loading in roots and inter-vascular transfer in nodes to deliver Zn/Cd to grain in rice. Biochemical and Biophysical Research Communications, 2019, 512(1): 112-118. DOI:10.1016/j.bbrc.2019.03.024
[43]
Takahashi R, Ishimaru Y, Shimo H, Ogo Y, Senoura T, Nishizawa N K, Nakanishi H. The OsHMA2 transporter is involved in root-to-shoot translocation of Zn and Cd in rice. Plant, Cell & Environment, 2012, 35(11): 1948-1957.
[44]
Sasaki A, Yamaji N, Ma J F. Overexpression of OsHMA3 enhances Cd tolerance and expression of Zn transporter genes in rice. Journal of Experimental Botany, 2014, 65(20): 6013-6021. DOI:10.1093/jxb/eru340
[45]
Sato H, Shirasawa S, Maeda H, Nakagomi K, Kaji R, Ohta H, Yamaguchi M, Nishio T. Analysis of QTL for lowering cadmium concentration in rice grains from 'LAC23'. Breeding Science, 2011, 61(2): 196-200. DOI:10.1270/jsbbs.61.196