全 文 : 福建省教育厅基金项目(K20024 、K02060) 、福建省省长基金(2002)资助
收稿日期 :20050424 改回日期 :20050719
基于 RAPD标记的芸苔属蔬菜茎叶镉累积特性与耐性分析
王松良
(福建农林大学作物科学学院 福州 350002)
摘 要 应用 RAPD技术分析 13个小白菜和 11个结球甘蓝品种的 Cd累积特性与耐性的结果表明 ,8个随机引物
中有 5个引物扩增到 31条多态性片段 ,大小在 100 ~ 2000bp之间 。聚类分析表明 ,茎叶 Cd高累积的品种有聚集
成丛的趋势 ;茎叶 Cd含量与品种间相对遗传距离呈极显著的线性回归关系 ;遗传多样性则随着茎叶 Cd含量的增
加而升高 ;高 、低 Cd累积组间分子方差(AMOVA)也呈现显著差异 。初步确认了与供试材料茎叶 Cd高累积和耐性
相关的 2个 RPAD标记 :小白菜 1个(350bp) ,结球甘蓝 1个(410bp) 。
关键词 随机扩增多态性 DNA(RAPD) 芸苔属蔬菜 Cd累积 遗传多样性 分子方差(AMOVA)
A RAPD analysis of shoot cadmium accumulating characteristic and tolerance of vegetable cultivars from Brassica
genus .WANG SongLiang (College of Crop Sciences ,Fujian Agriculture and Forestry University ,Fuzhou 350002 ,Chi
na) ,CJEA ,2006 ,14(4) :141 ~ 145
Abstract 13 cultivars of Brassica Chinensis and 11 of B . oleracea were analyzed using Random Amplified Polymorphic
DNA (RAPD) markers to investigate the pattern of genetic variation related to Cd accumulation and tolerance of shoot un
der Cd stress .31 wellamplified and highly reproductive polymorphic bands by 5 of total 8 random primers were scored ,
with length between 100 ~ 2000bp .Two polygenic trees ,of B . Chinensis and B . oleracea , respectively ,computed us
ing microevolution method show that the cultivars with high shoot Cd accumulation present tight clustering tendence ,and
it is of highly significant linear regression between the shoot Cd content and their correspondingly genetic distance . The
genetic diversity is found higher in subgroup of high shoot Cd content than that in subgroup of low shoot Cd content .
Analysis of Molecular Variance (AMOVA) of the RAPD markers give highly statistical significance to the subgroupings of
populations according to the shoot Cd content ,indicating that RAPD markers detect a polymorphism related to the Cd ac
cumulation .Finally , two RAPD bands ,one for B . chinensis (350bp) ,another for B . oleracea (410bp) ,presented in
high shoot Cd content subgroup ,but exclusive to the low one are identified .
Key words RAPD , Brassica vegetable ,Cd accumulation ,Genetic diversity ,AMOVA
(Received April 24 ,2005 ; revised July 19 ,2005)
我国农用土壤和水体普遍受到重金属的污染[1] 。研究表明在自然条件下 ,不同植物对重金属污染物的
适应可引起植物微进化(Microevolution)[5] 。耐性植物与非耐性植物在酶多态性等方面存在明显的差异[6] 。
DNA随机扩增多态性(Random amplified polymorphic DNA ,RAPD) 分析提供一个基于 DNA 的标记[7] ,是
揭示植物的遗传变异和重金属耐性之间关系的可能途径 。 目前发现的重金属高累积 、高耐性植物在分类上
大多为芸苔属植物[8] 。理论上 ,芸苔属蔬菜可作为植物对重金属耐性的生态遗传学和系统生物学研究的模
式植物 ;在应用上 ,小白菜 、结球甘蓝 、芥菜等芸苔属蔬菜是全国各地广为种植的植物食品 。 研究芸苔属蔬
菜的重金属耐性有望获得潜在的土壤重金属植物修复材料和低重金属累积蔬菜品种 ,这对保障蔬菜食品安
全意义重大[2] 。 Cd是环境中最重要的重金属污染物 ,它被认为是人类的致癌物(Carcinogen)之一 。 同时 Cd
通过干扰生物 DNA的碱基配对引起生物包括人体 DNA的变异 ,对生物危害极大 。 本研究旨在揭示芸苔属
蔬菜茎叶 Cd累积量和耐性的差异及其遗传变异方式 ,同时通过种内 Cd耐性区分 ,寻求与植物 Cd 耐性连锁
的分子标记 ,洞见植物因重金属胁迫引起的遗传变异与微进化的关系 ,以进一步评价芸苔属植物作为植物
修复土壤重金属污染物的潜力 。
第 14卷第 4期 中 国 生 态 农 业 学 报 Vol .14 No .4
2 0 0 6年 1 0月 Chinese Journal of EcoAgriculture Oct ., 2006
1 实验材料与方法
供试材料为“日本冬妃” 、“虹桥矮青” 、“早生华京” 、“上海青” 、“福绿 1 号” 、“葵扇白菜” 、“94N1” 、“清江
白菜” 、“黄牙菜” 、“日本华冠” 、“黑叶白菜” 、“半青黄” 、“山西白”等 13 个小白菜品种(处理代码依次为 G11 ~
G113)和“Fieldfirce” 、“Morris” 、“ Surprice” 、“京丰 1号” 、“Matsumo” 、“Gideon” 、“Ducati” 、“Invinto” 、“夏秋 3
号” 、“中甘 11号” 、“Pacel”等 11个结球甘蓝品种(处理代码依次为 G214 ~ G224) ,实验在福建农林大学生
物技术中心温室进行 。植物材料的培育参考文献[2]并略做修改 :30d 蛭石培养苗在标准 Hoagland 营养液
预培养 10d后 ,移入含 50μmol/L Cd的 Hoagland营养液中 ,每品种设置 3 个重复 ,1 个对照 。 胁迫培养 10d
后 ,3个重复混合取样[3]用于 DNA提取 。茎叶 Cd含量测定按 3个重复分别取样 、测定 。
基因组 DNA 的提取按微量 CTAB 法[9]略加改进 :05g 新鲜叶片加液氮研磨成粉 ,粉末解冻前转入
2mL预热提取液[1 % CTAB (W/V) ,100 mmol/L TrisHCl(pH80) ,14mmol/L NaCl ,20mmol/L EDTA ,
01 % (V/V)巯基乙醇]中 ,65 ℃ 水浴 30min 后 ,搅拌混匀 ,加 15mL 氯仿/异戊醇(24∶1) ,反复倒置混匀
15s ,室温下以 5000r/min离心 5min ,收集上清液 ,加等体积冷异丙醇 ,混匀 10 ~ 20s 。 DNA沉淀后 ,用灭菌后
的牙签或枪头挑取 DNA ,沉淀物用 70 %乙醇漂洗 ,用滤纸吸干乙醇 。 溶于 250 ~ 500μL TE 缓冲液 ,65 ℃ 水
浴溶解 。加 1/10体积 3mmol/L NaAc (pH68)和 2 倍体积 95 % 乙醇 ,沉淀 DNA 。 挑出 DNA ,溶于 250μL
TE溶液中 ,分装待用 。
RAPD体系的建立 。 经预先优化后的 RAPD 体系为 25μL ,含 10mmol/L KCl ,2mmol/L TrisHCl
(pH75) ,001 % Triton X100 ,3mmol/L MgCl2 ,02mmol/L dNTP mix (上海生工) ,1U TaqDNA poly
merase (Promega) ,500 ng 引物 ,25 ng 模板 DNA 。 8 个 10pb 的随机引物引用 Mengoni 等[10]的类似研究 。
PCR在 MJ PTC200 型 PCR多功能扩增仪(基因公司生产)上进行 ,程序参见 Mengoni 等的方法[10] :94 ℃
3min ,经 40 个循环(94 ℃ 1min ,36 ℃ 2min ,72 ℃ 15min)后 ,72 ℃ 延伸 10min ,最后保持在 4 ℃ 。 扩增产物
在 BIORAD mini电泳仪上进行 ,2 %琼脂糖 ,电压 10V/cm ,电泳时间 15h ;EB染色 ,在 BIORAD凝胶成像
系统观察并紫外拍摄 。
RAPD扩增条带的检测与统计分析 。应用 Quality one软件(BIORAD)机检 ,并辅以手工检带除去因胶
板上亮点所引起的计算机误检 。清晰的电泳条带用于统计分析 ,按扩增带的有无记数 ,预先假设每个 RAPD
条带代表独立的遗传位点 ,当某一扩增带出现时 ,赋值为“1” ,不存在时赋值为“0” ,并记录在 Microsoft excel
表格上 。应用 Phylib软件包(version 36b)[11]的软件 ,计算出的遗传距离矩阵 ,导入 Mega21 软件[12] ,应用
微进化(ME)聚类 。根据聚类结果 ,计算每个品种的相对遗传距离 ,及其与植株茎叶 Cd 含量的回归关系 。
经分类的 RAPD标记数据导入 Arlequin 20 软件包 ,按欧氏平方距离法进行样品组之间遗传多样性(Genetic
diversity) 、分子方差分析(Analysis of molecular variance ,AMOVA) 及其显著测验(φtest)[13] 。
2 结果与分析
21 DNA提取和 RAPD扩增结果
按经改进后 CTAB法提取小白菜和结球甘蓝的基因组 DNA ,其电泳检测结果显示 ,获得了较为理想的
RAPD扩增模板 ,两类实验材料基因组长度约为 21000bp(图 1) 。 应用优化后的 PCR体系分别扩增 24 个处
图 1 供试材料基因组 DNA
Fig1 The molecular size of genome DNA
M ,DNA marker — λDNA/Eco .RI + Hind Ⅲ (上海生工 ,Sangon) ;
1 ,2 为小白菜基因组 DNA ,3 ,4 为结球甘蓝基因组 DNA 。
图 2 引物 1247对 24个供试 DNA样品的扩增条带
Fig2 Amplied bands of 24 DNA samples by primer 1247
M ,DNA marker — λDNA/Eco .RI + Hind Ⅲ (上海生工) ;1 ~ 13 ,14 ~ 24 分别
为引物 1247 对小白菜 G11 ~ G113 和结球甘蓝 G214 ~ G224
的基因组 DNA 样品的扩增条带 。
142 中 国 生 态 农 业 学 报 第 14卷
理 DNA样品 ,8 条引物中的 5 条扩增
出共 35 条清晰的多态性片段 ,其中 ,
具有高复制性的多态性条带 31 条 ,平
均每个引物可扩增出 6 个多态性条
带 ,多态性达 89 % 。片段长度在100 ~
2000bp之间(见表 1) ,可用作 RAPD位
点分析 。图 2为引物 1247对 24个处理
样品扩增条带 。
表 1 用于RAPD分析的随机引物序列及相应的扩增多态性条带数目和碱基程度
Tab1 Oligonucleotide primers used for RAPD analysis ,
number and size range of polymorphic bands scored
引物代码 序列 多态性条带数目 碱基长度范围/bp
Primer codes Sequence (5’ — 3’) No .of polymorphic bands Molecular size range
1253 GT TCCGCCCC 3 550 ~ 100
1247 AAGAGCCCGT 9 2000 ~ 110
SP1 CAGGCCTCCG 8 1300 ~ 260
SP2 GGGCATCGGC 6 560 ~ 150
OPP05 CAGGCCTCCG 5 1500 ~ 300
总计 31 100 ~ 2000
图 3 基于 RAPD标记的小白菜品种茎叶 Cd含量与品种间相对遗传距离的关系
Fig3 Microevolutionary tree and relationship between correspondingly genetic distance and
Cd tolerance of B .chinensis based on RAPD analysis
茎叶 Cd 含量数据为 3 个重复的平均值 ,其后英文字母表示 5 % 水平的差异显著性 ,图 4 同 。
22 RAPD标记与
Cd 累积特性
和耐性的关系
通过 Phylib软
件运算分别获得
13 个小白菜和 11
个结球甘蓝品种间
的相对遗传距离 。
经 Mega21 软 件
按 ME法聚类后分
别得到进化树图 ,
同时把对应品种相
对遗传距离和茎叶
Cd 含量列在进化
树图右边(见图 3 、
4) 。通过相关和线
性回归分析表明 ,
两类供试蔬菜品种
的相对遗传距离
图 4 基于 RAPD标记的结球甘蓝品种茎叶 Cd含量与品种间相对遗传距离的关系
Fig4 Microevolutionary tree and relationship between correspondingly genetic distance
and Cd tolerance of B .oleracea based on RAPD analysis
( X)与茎叶 Cd 含
量( Y)呈极显著线性回归
关系 ,回归方程小白菜为 :
Y = 852656 X + 86318
( R = 0946 ,p < 0001)
(1)
结球甘蓝为 :
Y = 669691 X + 173661
( R = 0937 ,p < 0001)
(2)
上述结果表明 ,基于
RAPD 标记的相对遗传
距离可作为反映植物的
重金属耐性指标 。 根据
上述结果 ,分别把小白
菜 、结球甘蓝按茎叶 Cd
含量分为高 、低累积型 2
个亚组(Subgroup) ,应用
第 4期 王松良 :基于 RAPD标记的芸苔属蔬菜茎叶镉累积特性与耐性分析 143
表 2 Cd耐性分组及其遗传多样性估计
Tab2 Classifications of Cd tolerance and their genetic diversity estimation
Cd 累积类型 Cd 耐性类型分组 遗传多样性
Cdaccumulating type Subgroups divided according to Cdaccumulation Genetic diversity
小白菜
高 G11 ,G12 ,G15 ,G16 ,G17 ,G18 ,G110 0465
低 G13 ,G14 ,G19 ,G111 ,G112 ,G113 0181
结球甘蓝
高 G214 ,G217 ,G219 ,G220 ,G222 0512
低 G213 ,G215 ,G216 ,G218 ,G221 ,G223 ,G224 0387
Arlequin 20软件运算各亚组的遗传多样
性(见表 2) 。 结果表明 ,高 Cd 累积亚组
的遗传多样性相应比低 Cd 累积亚组高 。
其中 ,小白菜高 Cd累积亚组的遗传多样性
为 0465 ,比低累积亚组高 1569 % ;结球
甘蓝高 Cd累积亚组则为 0512 ,比低累积
亚组高 323 % 。 从 ME 聚类结果也反映
出 ,在 Cd胁迫下 ,不管是小白菜还是结球
甘蓝 ,除了个别基因型外(如小白菜 G110
和结球甘蓝的 G219) ,不同植物基因型在微
进化中 ,有按对环境适应力聚集成丛(Cluster)
的趋势 。遗传多样性的降低可能意味着在进
化中对不良环境适应力处于劣势的地位[4] 。
用 AMOVA检验植物 Cd耐性与遗传变
异之间差异[12] ,分析结果表明(见表 3) ,小白菜
与结球甘蓝类间及 24 个品种间差异不显著
(φCT = - 0019 ,P < 0342 ;φST = 0059 ,P <
0055) ,而类内组间差异达显著水平( φSC =
表 3 31条RAPD多态性条带的分子方差分析
Tab3 Analysis of Molecular Variance(AMOVA) for 24
samples using 31 RAPD polymorphic bands
变异
Source of variation
自由度
D .f .
平方和
Sum of squares
方差
Variance
测验
statistics test P
类间 1 7493 - 0122 CT = - 0019 < 0342
类内组间 2 17855 0498 SC = 0077 < 0046
属内品种间 20 119152 5958 ST = 0059 < 0055
总计 23 144500 6334
0077 ,P < 0046) 。由此进一步确认 ,Cd 胁迫导致遗传变异 ,致使不同植物个体基因型在 Cd 耐性出现分
异 ,最终导致抗 Cd污染生态型的形成[4] 。
对与 Cd耐性连锁的 RAPD 标记的分析表明 ,在本试验条件下 ,31 条多态性扩增带中 ,350bp(引物
1247)仅出现在高 Cd累积亚组中 ,而排除在低 Cd累积亚组之外 ;甘蓝组中 ,410bp(引物 SP1)仅出现在高 Cd
累积亚组中 ,而排除在低 Cd累积亚组之外 。它们可被初步认定为与茎叶 Cd高累积连锁的 RAPD标记 。
3 小结与讨论
Cd是一种人体和植物体不必需的重金属元素 ,然而 ,由于自然和人为原因 ,农业土壤中大量存在 Cd 等
重金属污染物 ,经食物链浓缩后对人体和生物具有很强的危害性 。 研究蔬菜对 Cd 不同累积特性和耐性具
有现实意义 。植物种群遗传组成的差异是选择不同 Cd 累积类型蔬菜的理论前提 。 为了证实该理论假设 ,
RAPD标记检测是一条可能的技术途径 。近年来 ,越来越多的研究者把 RAPD 技术应用于环境监测方面的
研究[14] 。 Koch等 1999年的研究指出 ,重金属胁迫下的拟南芥与对照呈现不同的 RAPD 多态性 ,意味着应
用 RAPD可探寻环境胁迫下的植物遗传多样性变化[12] 。本实验结果表明 ,RAPD标记可以有效地展示与植
物 Cd耐性或抗性相关的分子多态性 ,进一步为人工探寻或室内筛选不同 Cd累积和耐性类型植物提供理论
依据 。自然环境的变异和人工干扰都可能导致自然物种遗传多态性的变化[15] 。本实验结果表明 ,重金属污
染可以促进植物物种遗传多样性丧失 ,影响其微进化进程 。 具体表现一是高 Cd 累积型的芸苔蔬菜类型与
其基于 RAPD的相对遗传距离成正相关 ,这个结果与 Ducousso等(1990)[16]对植物种群( A r rhenantherium)
基于形态和同工酶变异的研究一致 ;二是 AMOVA 和遗传多样性分析表明 ,具有高累积或耐性的品种群体
往往具有高的遗传多样性和分子方差 ,φ测验也表明其与低累积的群体之间呈现显著的差异 ;三是关于遗传
多样性的分布 ,ME聚类结果(见图 3 ,4)表明 ,除个别品种外 ,高累积的品种总是聚集在一起 ,低累积群体遗
传多样性丧失可能意味者其耐性的降低[17] 。 此外 ,本实验还发现与 Cd 高累积连锁的 2 个 RAPD 标记 ,即
350bp(引物 1247)和 410bp(引物 SP1) 。这 2个标记尚待进一步证实 。 总之 ,利用 RAPD标记分析不同重金
属污染背景下作物 DNA 多态性 ,辨认与重金属耐性基因连锁的 DNA 片段 ,并进一步克隆 、转化和表达 ,以
改善植物的重金属耐性 ,促进植物修复土壤重金属污染技术走向田间应用[18] 。
致谢 承蒙周元昌博士提供部分试验材料 ,福建农林大学农学 2000 级林颖慧 、范兆烽同学参与部分试验操
作 ,植株重金属的测定得到林梅讲师的热心帮助 。在此一并表示诚挚的谢意 。
144 中 国 生 态 农 业 学 报 第 14卷
参 考 文 献
1 陈怀满 ,郑春荣 ,涂 从等 .中国土壤重金属污染现状与防治对策 . Ambio ., 1999 ,28(2) :130 ~ 134
2 王松良 ,郑金贵 .芸苔属蔬菜的 Cd 富集特性及其修复土壤 Cd 污染的潜力 .福建农林大学学报(自然科学版) ,2004 ,33(1) :94 ~ 99
3 江树业 ,陈启锋 ,方宣钧 .利用 cDNARAPD 分析籼稻光敏核不育基因的差异表达 .科学通报 ,1998 ,43(23) :2521 ~ 2524
4 熊治廷 .植物抗污染进化及其遗传生态学代价 .生态学杂志 ,1997 ,16(1) :53 ~ 57
5 Lefebvre C .,Vernet P .Microevolutionary processes on contaminated deposits . In :Shaw J .(ed .) .Heavy Metal Tolerance in Plants : Evolu
tionary Aspects . Boca Raton ,FL : CRC , 1990 .286 ~ 297
6 Macnair M .R .The genetics of metal tolerance in vascular plants . New Phytologist ,1993 ,124 : 541 ~ 559
7 Willams J . G .K ., Kubelik A .R .K ., Livak J .,et al . DNA polymorphisms amplified by arbit rary primers are useful as genetic markers .
Nucleic . Acids . Res .,1990 ,18 :6531 ~ 6535
8 Zhu Y .L ., Elizabeth A .H ., Smit P .,et al . Cadmium tolerance and accumulation in Indian mustard is enhanced by overexpressing γglu
tamylcysteine synthase . Plant Physiology , 1999 ,121 : 1169 ~ 1177
9 Saghai M .M .A .,Soliman K .M .,Jorgensen R .A .,et al .Ribosomal DNA spacer length polymorphism in population dynamics .Proceedings
of the National Academy of Sciences ,USA ,1998 ,81 : 8014 ~ 8018
10 Mengoni A .,Gonnelli C .F .,Galardi R .,et al .Genetic diversit y and heavy metal tolerance in populations of Silene paradoxaL .(Caryophyl
laceae) : a random amplified polymorphic DNA analysis .Molecular Ecology , 2000 ,9 : 1319 ~ 1324
11 Felsentein J . Phylogenies from rest riction sites , a maximum likelihood approach . Evolution , 1992 ,46 :159 ~ 173
12 Koch M .,Bishop J ., Mitchell O .T .Molecular systematics and evolution of A rabidopsis and A rabis . Plant Biology , 1999 ,1 :529 ~ 537
13 Excoffier L ., Smouse P ., Quatt ro J . M . Analysis of molecular variance inferred from metric distances among DNA haplotypes :application
to human mitochondrial DNA rest riction data . Genetics , 1992 , 131 :479 ~ 491
14 Slbiger R .N ., Christ S .A ., Leonard A . C ., et al . Preliminary studies on the population genetics of the cent ral stoneroller( Cam postoma
anomalum) from the Great Miami River Basin , Ohio . Environ . Monitor . Assess ., 1998 , 51 : 481 ~ 495
15 Hedrik W ., Ginevan M . E ., Ewing E .P . Genetic polymorphism in heterogeneous environment . Annual Review of Ecology and Systemat
ics , 1976 ,7 :1 ~ 32
16 Ducousso A ., Petit D .,Valero M ., et al .Genetic variation bet ween and within populations of a perennial g rass :A r rhenantherium elatius .
Heredity , 1990 ,65 :179 ~ 188
17 Van den Koornhuyse P ., Leyval C .,Bonnin I . High genetic diversity in arbuscular mycorrhizal fungi : evidence for recombination events .
Heredity , 2001 , 87 : 243 ~ 253
18 Waugh R ., Russell J . Screening for ‘useful variation’ using RAPD markers . In : Karp A ., Isaac P .G , Ingram D .S .( eds .) .Molecular
Tools for Screening Biodiversity . London : Chapman & Hall ,1998 . 23 ~ 57
欢迎订阅 2007年枟农业现代化研究枠
枟农业现代化研究枠是由中国科学院农业研究委员会和中国科学院亚热带农业生态研究所主办的农业综合性学术刊物 ,
科学出版社出版 ,主要刊登农业发展战略和农业基础科学及其交叉学科的理论研究 、基础研究和应用研究方面的学术论文 、
科研报告 、研究简报和文献综述等 ,包括农业发展战略 、农业可持续发展 、区域农业 、生态农业 、农业生物工程 、信息农业 、农村
生态环境 、农业经济 、农业产业化 、农业系统工程 、农业机械化 、高新技术应用 、资源利用与保护和国外农业等 ,适于广大农业
科技工作者 、农业院校师生 、各级领导干部和管理人员 、农业工作者阅读 。本刊为双月刊 ,逢单月出版 ,国际标准大 16开本 ,
国内外公开发行 ,国内邮发代号 :4246 ,国外发行代号 :BM6665 ,每册定价 1200元 ,全年 7200元 ,全国各地邮局(所)均可订
阅 ,漏订者也可直接汇款至本刊补订 ,地址 :(410125)湖南省长沙市马坡岭中国科学院亚热带农业生态研究所枟农业现代化研
究枠编辑部 ;电话 :07314615231 ;Email :nyxdhyj @ isa .ac .cn 。
第 4期 王松良 :基于 RAPD标记的芸苔属蔬菜茎叶镉累积特性与耐性分析 145