免费文献传递   相关文献

Preliminary Analysis on the Genetic Diversity of Myriophyllum spicatum from China

中国穗状狐尾藻遗传多样性的初步分析



全 文 :植物科学学报  2016ꎬ 34(1): 109~116
Plant Science Journal http: / / www.plantscience.cn
DOI:10􀆰 11913 / PSJ􀆰 2095 ̄0837􀆰 2016􀆰 10109
吕天锋ꎬ 蒿飞ꎬ 熊争ꎬ 刘星. 中国穗状狐尾藻遗传多样性的初步分析[J] . 植物科学学报ꎬ 2016ꎬ 34(1): 109-116
Lü TFꎬ Hao Fꎬ Xiong Zꎬ Liu X. Preliminary analysis on the genetic diversity of Myriophyllum spicatum from China[J] . Plant Science Jour ̄
nalꎬ 2016ꎬ 34(1): 109-116
中国穗状狐尾藻遗传多样性的初步分析
吕天锋ꎬ 蒿 飞ꎬ 熊 争ꎬ 刘 星∗
(武汉大学生命科学学院ꎬ 植物系统学与进化生物学实验室ꎬ 武汉 430072)
摘  要: 利用叶绿体 DNA三个片段( trnK ̄matK、 trnL ̄trnF、 rpl 32 ̄trnL)对中国大陆广布的水生植物穗状狐尾藻
(Myriophyllum spicatum L.)的遗传多样性进行了初步分析ꎬ 以探讨其自然居群的遗传结构及具有广泛分布格局
的可能机制ꎮ AMOVA 分析显示ꎬ 穗状狐尾藻 8 个居群间的遗传变异为 84􀆰97%ꎬ 而居群内的遗传变异为
15􀆰03%ꎬ 居群遗传分化系数(Fst)为 0􀆰85ꎬ 表明穗状狐尾藻具有较高的遗传多样性(Hd = 0􀆰83)且主要存在于居
群间ꎬ 奠基者效应可能导致了最初的遗传差异ꎬ 而隔离障碍(Nm =0􀆰09)又进一步导致了居群间的遗传分化ꎮ 基
于 17个单倍型构建的系统发育树和网络关系图均显示ꎬ 单倍型 H5 和 H6 在居群中的分布范围最广且出现频率
最高ꎬ 表明 H5和 H6可能为最古老的祖先单倍型ꎮ Mantel检验表明居群间的遗传距离与地理距离之间不存在显
著的相关性ꎬ 失配分布检测结果显示穗状狐尾藻在历史上曾发生过扩张事件ꎬ 而 Tajima’s、 Fu & Li’s D∗和 F∗
检测发现ꎬ 该物种不存在明显的谱系地理格局ꎬ 这可能与穗状狐尾藻种子的长距离扩散有关ꎮ
关键词: 穗状狐尾藻ꎻ 叶绿体 DNA片段ꎻ 遗传多样性
中图分类号: Q943ꎻ Q949􀆰762􀆰6          文献标识码: A          文章编号: 2095 ̄0837(2016)01 ̄0109 ̄08
      收稿日期: 2015 ̄09 ̄02ꎬ 退修日期: 2015 ̄10 ̄11ꎮ
  基金项目: 国家自然科学基金面上项目(31170203)ꎮ
This work was supported by a grant from the National Natural Science Foundation of China (31170203) .
  作者简介: 吕天锋(1989-)ꎬ 男ꎬ 硕士研究生ꎬ 研究方向为系统与进化植物学(E ̄mail: lvtianfeng@whu􀆰 edu􀆰 cn)ꎮ
  ∗通讯作者(Author for correspondence􀆰 E ̄mail: xingliu@whu􀆰 edu􀆰 cn)ꎮ
Preliminary Analysis on the Genetic Diversity of
Myriophyllum spicatum from China
LÜ Tian ̄Fengꎬ HAO Feiꎬ XIONG Zhengꎬ LIU Xing∗
(Laboratory of Plant Systematics and Evolutionary Biologyꎬ College of Life Sciencesꎬ Wuhan Universityꎬ Wuhan 430072ꎬ China)
Abstract: To determine the genetic structure of natural Myriophyllum spicatum populations in
mainland China and the possible mechanism that led to the widely distributed pattern of the
speciesꎬ eight populations were examined for preliminary genetic analysis based on three
chloroplast DNA segments ( trnK ̄matKꎬ trnL ̄trnFꎬ rpl 32 ̄trnL) . The AMOVA analysis results
showed that the genetic diversity among the eight populations was 84􀆰97% and within the
populations was 15􀆰 03%. Furthermoreꎬ the population genetic differentiation coefficient (Fst)
was 0􀆰85ꎬ indicating high genetic diversity in M. spicatum (Hd = 0􀆰83) among populations.
The founder effect may have led to the initial genetic differences and segregation barriers
(Nm = 0􀆰09) may have further resulted in genetic differentiation among populations. From the
constructed phylogenetic tree and haplotype network based on 17 haplotypesꎬ haplotypes H5
and H6 had the highest frequency and widest distributionꎬ indicating they were possibly the
oldest ancestral haplotypes. The Mantel test showed no significant correlation between genetic
and geographic distances among populationsꎬ and the mismatch distribution test showed that
historic M. spicatum populations likely experienced expansion events. Tajima’sꎬ Fu & Li’s D∗
and F∗ tests showed that this species did not have an obvious phylogeographic patternꎬ which
may have resulted from long distance dispersal of seeds.
Key words: Myriophyllum spicatumꎻ Chloroplast DNA segmentsꎻ Genetic diversity
    穗状狐尾藻(Myriophyllum spicatum L.)为小
二仙草科(Haloragidaceae)狐尾藻属多年生水生
草本植物ꎬ 其自然居群主要分布于欧亚大陆和北
非ꎮ 近年来ꎬ 国内外有关穗状狐尾藻的研究主要
集中在生态与生理学方面ꎬ 内容涉及形态可塑
性[1] 、 环境因子对生长发育的影响[2-5] 、 对水环
境污染的响应及繁殖对策[6ꎬ7]等ꎮ 目前ꎬ 穗状狐
尾藻在北美地区的生态入侵问题受到广泛关注ꎬ
作为入侵种它可通过与当地的土著种杂交进一步
强化其生态入侵能力[8-11] ꎮ 然而ꎬ 有关穗状狐尾
藻的遗传多样性评估及其世界性分布格局的成因
一直研究较少ꎮ
DNA测序分析即对 DNA片段的核苷酸排列顺
序进行测定ꎬ 是最能够充分揭示 DNA 多样性的一
种研究方法ꎮ 叶绿体 DNA 在大多数的被子植物中
具有稳定的母系遗传特性ꎬ 其基因组中存在很多相
对保守的序列[12]并为通用引物的设计提供了依据ꎬ
因此对叶绿体 DNA的特定片段进行 PCR扩增和测
序越来越被广泛应用于植物遗传多样性分析[13]、
系统发育重建[14]以及谱系地理学研究中[15-17]ꎮ
穗状狐尾藻广泛分布在我国南北各省区ꎬ 为土
著种ꎬ 常见于河流、 池塘、 小溪、 湖泊、 沼泽等水
生环境中[18]ꎮ 本研究拟以采自中国大陆的穗状狐
尾藻 8个居群为研究对象ꎬ 运用叶绿体 DNA 片段
( trnK ̄matK、 trnL ̄trnF、 rpl32 ̄trnL)对其遗传多样
性进行初步分析ꎬ 旨在探讨穗状狐尾藻自然居群的
遗传结构及其广泛分布格局形成的可能机制ꎮ
1  材料与方法
1􀆰 1  材料
穗状狐尾藻的 8个自然居群分别采自中国大陆
8个省、 自治区(黑龙江省、 内蒙古自治区、 湖北
省、 甘肃省、 西藏自治区、 云南省、 山东省、 广东
省)ꎬ 每个居群随机选取 5~6株个体ꎬ 个体间间距
至少 3 mꎬ 共取样 47份(表 1)ꎮ 将每株个体新鲜、
幼嫩的叶片洗干净后ꎬ 用吸水纸尽量吸干其表面水
分ꎬ 然后放入置有变色硅胶的取样袋中迅速干燥ꎬ
带回实验室后置于-20℃冰箱中长期保存、 备用ꎮ
表 1  穗状狐尾藻居群的采集信息
Table 1  Detailed information on the Myriophyllum spicatum populations examined in this study
居群
Populations
纬度
Latitude
经度
Longitude
个体数
No. of
individuals
生境
Habitat
伴生种
Companion species
黑龙江省五大连池(WD)
Wudalianchiꎬ Heilongjiang 48°39′12.26″ N 126°08′47.48″ E 6
湖泊
Lake 无
内蒙古自治区多伦县(DL)
Duolunꎬ Inner Mongolia 42°18′38.10″ N 116°29′32.70″ E 6
河流
River
金鱼藻 (Ceratophyllum demersum)、
篦齿眼子菜 ( Potamogeton pectina ̄
tus)
湖北省鄂州市梁子湖(EZ)
Liangzi Lakeꎬ Hubei 30°14′21.90″ N 114°38′56.10″ E 6
湖泊
Lake
菱(Trapa bispinosa)、 竹叶眼子菜(P.
malaianus)、 狐尾藻(M. verticillatum)
甘肃省兰州市小西湖(LZ)
Lanzhouꎬ Gansu 36°03′50.58″ N 103°47′48.31″ E 6
池塘
Pond
竹叶眼子菜(P. malaianus)、光叶眼子
菜(P. lucens)
西藏自治区米林县(ML)
Milinꎬ Tibet 29°11′10.05″ N 94°11′36.40″ E 6
池塘
Pond 无
云南省澄江县抚仙湖(CJ)
Chengjiangꎬ Yunnan 24°37′17.10″ N 102°56′39.42″ E 5
湖泊
Lake
竹叶眼子菜(P. malaianus)、篦齿眼子
菜(P. pectinatus)
山东省枣庄市(ZZ)
Zaozhuangꎬ Shandong 34°34′34.00″ N 117°29′52.53″ E 6
河流
River
竹叶眼子菜( P. malaianus)、 金鱼藻
(C. demersum)、 黑藻(Hydrilla verti ̄
cillata)
广东省英德市(YD)
Yingdeꎬ Guangdong 24°14′46.50″ N 113°38′55.50″ E 6
溪流
Stream 无
011 植 物 科 学 学 报 第 34卷 
1􀆰 2  基因组 DNA的提取
参考文献[19]的方法提取硅胶干燥叶片的基
因组 DNAꎬ 并利用 1%琼脂糖凝胶电泳进行 DNA
质量检测ꎮ
1􀆰 3  叶绿体 DNA片段的 PCR扩增与测序
叶绿体 DNA三个片段的 PCR扩增引物(表 2)
由生工生物工程(上海)有限公司合成ꎮ PCR扩增在
Eppendorf AG 22331 Hamburg PCR仪(德国)上进
行ꎮ PCR 反应体系为 25 μLꎬ 主要包含 25 ng 模板
DNA、 2􀆰5 mmol / L dNTP 2 μL (Takara)、 2􀆰 5 μL
10 × PCR Taq Buffer[含 100 mmol / L Tris ̄HCl (pH
8􀆰3)、 500 mmol / L KCl、 15 mmol / L MgCl2 ]
(Takara)、 5 U / μL Taq酶(Takara)、 10 μmol / L扩
增引物各 0􀆰5 μLꎮ
叶绿体 DNA三个片段的 PCR 扩增产物经 1%
琼脂糖凝胶电泳检测ꎬ 均得到清晰、 明亮的单一条
带ꎬ 其中 trnK ̄matK扩增片段约为 2400 bpꎬ trnL ̄
trnF约为 900 bpꎬ rpl32 ̄trnL约为 1200 bpꎮ PCR
产物直接送往生工生物工程(上海)有限公司进行
测序ꎮ
1􀆰 4  数据分析
采用 Contig Express Application 软件对叶
绿体 DNA 片段的测序序列进行拼接后ꎬ 用
ClustalX1􀆰83 进行序列比对 [23] ꎬ 并辅以人工校
正ꎮ 由于叶绿体基因不存在重组现象ꎬ 因此我
们对本研究中的叶绿体 DNA 的 3 个片段序列进
行联合分析ꎮ 利用 DNasp version 5 软件计算核
苷酸多样性指数(π) 、 基因流(Nm)及居群间的
遗传分化系数 ( Fst)等 [24] ꎬ 同时通过 Tajima’s
D[25] 、 Fu & Li’s D∗和 F∗检测 [26ꎬ27]其是否符合
中性进化模型ꎻ 利用失配分布对居群是否经历
扩张或者瓶颈效应进行判断ꎮ 通过 Arlequin ver
3􀆰1 软件的 AMOVA 程序进行居群间和居群内的
遗传变异组成分析 [28] ꎮ 运用 GeneAlEx 6􀆰501
软件对遗传距离与地理距离之间的相关性进行
mantel检验 [29] ꎮ 采用 PAUP∗4􀆰0b10 软件并通
过 1000 个重复 ( replicates) 进行启发式搜索
( heuristic search ) 构 建 单 倍 型 最 大 简 约 树
(maximum parsimonyꎬ MP) [30] ꎮ 单倍型之间
的网络关系图用 Network 软件中 median ̄joining
方法 [31]进行构建ꎮ
2  结果与分析
2􀆰 1  穗状狐尾藻单倍型多样性及其分布
测序结果表明ꎬ 叶绿体 DNA 三个片段的序列
长度分别为 2318 ~ 2319 bp( trnK ̄matK)、 829 ~
831 bp( trnL ̄trnF)、 1289~1316 bp( rpl32 ̄trnL)ꎮ
合并分析并进行序列比对后ꎬ 叶绿体 DNA 三个片
段的序列长度为 4384 bpꎮ 用 DNasp 分析获得了
17个单倍型(H1~H17ꎬ 表 3)ꎬ 其中单倍型 H5和
H6分布频率最高ꎬ 其所在的居群分布范围最广ꎮ
对穗状狐尾藻 8个自然居群的单倍型多样性及核苷
酸多样性分析表明ꎬ 总体表现出较高的遗传多样性
(Hd = 0􀆰83ꎬ π =4×10-4)ꎬ 居群间遗传多样性差
异较大(Hd = 0 ~ 0􀆰87ꎬ π = 0~1 ×10-4)ꎬ 其中甘
肃兰州居群(LZ)的遗传多样性相对较高 (Hd =
0􀆰87ꎬ π = 1 ×10-4)ꎮ 西藏米林居群(ML)只有一
个单倍型 H5ꎬ 可能与其在高原特殊环境下的无性
繁殖有关ꎬ 而其余居群均具有不同程度的多态性
(包含 2 ~ 5个单倍型)ꎮ
表 2  叶绿体 DNA 3个片段的 PCR扩增引物
Table 2  PCR primers for three chloroplast DNA segments of Myriophyllum spicatum
DNA片段
Region
引物序列(5′-3′)
Primer sequence
退火温度(℃)
Annealing temp
循环数
Cycles
参考文献
Reference
trnK ̄matK
trnK ̄3914F: TGGGTTGCTAACTCA ATGG
trnK ̄2R: AACTAGTCGGATGGAGTAG
trnK ̄AR:CGCTTCTTT CTT TACGAGG
trnK ̄2W1F:TCGGTCGATGCGTTTACCTG
50 30 Johnson等[20]
trnL ̄trnF
c: CGAAATCGGTAGACGCTACG
f: ATTTGAACTGGTGACACGAG
55 30 Taberlet等[21]
rpl32 ̄trnL
trnLUAG: CTGCTTCTTAGGAGCAGCGT
rpl32: CAG TTCCAA AAA AACGTACTTC
50 40 Shaw等[22]
111  第 1期                        吕天锋等: 中国穗状狐尾藻遗传多样性的初步分析
2􀆰 2  穗状狐尾藻居群遗传结构及谱系地理历史
利用 Arlequin ver 3􀆰1 软件的 AMOVA 程序对
穗状狐尾藻居群间和居群内的遗传变异进行分析ꎬ
发现 8个居群间的遗传变异为 84􀆰97%ꎬ 而居群内
的遗传变异为 15􀆰03%ꎬ 基因流(Nm)为 0􀆰09ꎬ 居
群遗传分化系数(Fst)为 0􀆰85ꎬ 表明 8 个居群间的
遗传变异水平很高ꎬ 而居群内的遗传变异水平较
低ꎮ Mantel检验结果显示(图 1) ꎬ 居群间的遗传
!
"
#
$
(
)
ge
ne
tic
d
is
ta
nc
e
cM
%&#$
( )Geographic distance km
R Pxy = -0.23, rxy-rand rxy-data = 0.18( ≥ )
0.700
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0.0
00
50
0.0
00
10
00
.00
0
15
00
.00
0
20
00
.00
0
25
00
.00
0
30
00
.00
0
35
00
.00
0
40
00
.00
0
图 1  穗状狐尾藻居群间的遗传距离与
地理距离之间的 mantel检验
Fig􀆰 1  Mantel test between genetic and geographical
distances among Myriophyllum spicatum populations
距离与地理距离之间不存在显著的相关性(Rxy =
-0􀆰 23ꎬ P( rxy ̄rand ≥ rxy ̄data) = 0􀆰18)ꎮ 同时ꎬ
AMOVA 分析显示ꎬ Gst = 0􀆰43ꎬ Nst = 0􀆰41ꎬ 即
Gst >Nstꎬ 表明该物种不存在明显的谱系地理格局ꎮ
穗状狐尾藻 8个居群的 Tajima’s test分析结果为负
值ꎬ 且具显著水平(D = -2􀆰09ꎬ P < 0􀆰05)ꎬ 表明
其不遵循中性进化模型ꎻ Fu & Li􀆳s D∗和 F∗检测结
果均为负值ꎬ 且具极显著水平(D∗ = -3􀆰68ꎬ P <
0􀆰02ꎻ F∗ = -3􀆰71ꎬ P <0􀆰02)ꎬ 也表明其不遵循中
性进化模型ꎮ 失配分布检测结果显示(图 2)ꎬ 穗状
0
0 05.
0 10.
0 15.
0 20.
0 25.
0.30
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30
&()*
Pairwise difference
+
,
-
.
Fr
eq
ue
nc
y
!"#
$%#
Observed value
Expected value
图 2  穗状狐尾藻叶绿体 DNA序列的失配分布分析
Fig􀆰 2  Analysis of mismatch distributions based on
chloroplast DNA sequences of Myriophyllum spicatum
狐尾藻历史上曾发生过扩张事件ꎮ
2􀆰 3  穗状狐尾藻居群单倍型的系统发育关系
对居群单倍型进行系统发育关系分析ꎬ 结果显
示(图 3)各单倍型之间不存在明显的地理分布格
局ꎬ 除单倍型 H5 和 H6 出现频率较高外ꎬ 其余单
倍型均出现一次ꎮ 单倍型网络关系图也显示出与系
统发育树相似的结构(图 4)ꎬ 其中单倍型 H5 和
H6为频率分布最高的单倍型(表 3)ꎬ 推测这两个
单倍型可能为最古老的祖先单倍型ꎬ 其余单倍型为
后来分化并扩散或定居到其他居群中的ꎮ
3  讨论
水生植物相对于陆生植物分布广泛的一个重要
原因是其有着很强的扩散能力[32]ꎬ Sculthorpe[33]
认为水生植物的扩散中ꎬ 断枝扩散和种子流扩散占
据主要地位ꎮ 穗状狐尾藻可以进行快速的克隆繁
殖ꎬ 并随水流扩散到其他水域ꎮ 本研究中单倍型
H6为穗状狐尾藻的广东英德居群(YD)、 云南澄
江居群(CJ)、 山东枣庄居群(ZZ)和黑龙江五大连
H6 WD CJ ZZ YD( )
H16( )YD
H YD( )15
H YD( )14
H WD( )13
H LZ( )11
H CJ( )17
H DL( )3
H DL( )9
H DL( )8
H DL( )10
H1( )EZ
H LZ( )12
H EZ( )7
H ML ZZ LZ( )5
H LZ( )4
H2 EZ( )
80
括号中的字母组合同表 1 中的居群编号ꎻ 分支上数值代表简约
分析中大于 50的自展值ꎮ
Population codes in brackets are the same as those in Table 1ꎻ
Numbers above the branches are bootstrap values (> 50) from
parsimony analysis.
图 3  基于穗状狐尾藻叶绿体 DNA片段序列的 17个
单倍型构建的 MP系统发育树
Fig􀆰 3  Maximum parsimony tree constructed by 17
haplotypes based on chloroplast DNA sequences
of Myriophyllum spicatum
211 植 物 科 学 学 报 第 34卷 
H2
H1
H7
H12 H3
H9
H4 H5
H8
H11
H10
H13
H14
H16
H15
H6
H17
图 4  基于中国穗状狐尾藻叶绿体 DNA 3个片段
序列的 17个单倍型构建的网络关系图
Fig􀆰 4  Network of 17 haplotypes constructed from
chloroplast DNA sequences of Myriophyllum
spicatum from mainland China
池居群(WD)的共有单倍型(表 3ꎬ 图 3)ꎬ 表明克
隆繁殖和长距离扩散可能导致了穗状狐尾藻的广泛
分布ꎻ Santamaría[34]发现狐尾藻属植物的种子可
被水鸟食用ꎬ Chen等[35]对小二仙草科生物地理历
史的研究也表明ꎬ 在长距离的扩散中ꎬ 水鸟携带种
子的扩散起着最重要的作用ꎮ 本研究中除单倍型
H5在 3 个居群(ZZ、 ML、 LZ)、 H6 在 4 个居群
(YD、 CJ、 ZZ、 WD)中共有外ꎬ 其余 15 个单倍
型都是居群的特有单倍型(表 3ꎬ 图 3)ꎬ 因此种子
扩散产生的奠基者效应可能也对穗状狐尾藻的遗传
结构产生了重要影响ꎮ
有研究显示ꎬ 一些水生植物在居群间具有很高
的遗传变异ꎬ 在居群内的变异则比较低[36-38]ꎬ 如
毛茛属植物 Ranunculus nipponicus(居群间的遗
传变异为 84􀆰1%ꎬ 居群内变异为 15􀆰9%) [39]和野
生稻 Oryza granulata (居群间变异为 84􀆰9%ꎬ 居
群内变异为 15􀆰1% ) [40]等ꎻ 本研究中穗状狐尾藻
自然居群间的遗传多样性很高(84􀆰97%)ꎬ 居群内
的遗传变异则很低(15􀆰03%)ꎮ 对这些水生植物而
言ꎬ 隔离障碍可能限制了居群间的基因流ꎬ 从而导
致居群的进一步分化[34]ꎮ 例如ꎬ 穗状狐尾藻的 8
个自然居群之间存在较远的地理距离ꎬ 其间的山脉
等地理环境可能成为基因流的障碍(Nm = 0􀆰09)ꎬ
因此隔离分化可能是导致穗状狐尾藻遗传多样性主
要存在于居群间的重要因素ꎮ
单倍型 H5和 H6 在穗状狐尾藻的多个居群中
出现(表 3ꎬ 图 3)ꎬ 推测其可能是祖先单倍型ꎻ 失
配分布和中性检验也表明穗状狐尾藻在历史上曾经
有基因交流并经历过扩张事件(图 2)ꎬ 但居群间的
遗传距离和地理距离之间并不存在显著相关性(图
1)ꎬ 同时ꎬ 该物种也不存在显著的谱系地理结构
(Gst(0􀆰43) > Nst(0􀆰41))ꎮ 内蒙古多伦居群(DL)
表 3  基于叶绿体 DNA 3个片段联合序列分析的穗状狐尾藻遗传多样性指数
Table 3  Genetic diversity index of different Myriophyllum spicatum populations
based on chloroplast DNA sequence analysis
居群编号
Population
code
单倍型多样性
Haplotype
diversity (Hd)
核苷酸多样性
Nucleotide
diversity (π)
单倍型数量
Haplotypes
单倍型分布
Haplotype
distribution
WD 0.33 1.5 ×10-4 2 H6 (5)ꎬ H13 (1)
DL 0.80 2.2 ×10-4 4 H3 (1)ꎬ H8 (1)ꎬ H9 (3)ꎬ H10 (1)
EZ 0.60 0.0 3 H1 (1)ꎬ H2 (1)ꎬ H7 (4)
LZ 0.87 1.0 ×10-3 4 H4 (2)ꎬ H5 (2)ꎬ H11 (1)ꎬ H12 (1)
ML 0.00 0.0 1 H5 (6)
CJ 0.40 2.7 ×10-4 2 H6 (4)ꎬ H17 (1)
ZZ 0.33 7.0 ×10-5 2 H5 (1)ꎬ H6 (5)
YD 0.80 2.2 ×10-4 4 H6 (3)ꎬ H14 (1)ꎬ H15 (1)ꎬ H16 (1)
总计 Total 0.83 4.0 ×10-4 17
    注: 居群编号同表 1ꎻ 括号中数字表示个体数ꎮ
Note: Population codes are the same as those in Table 1ꎻ Numbers in brackets represent No. of individuals.
311  第 1期                        吕天锋等: 中国穗状狐尾藻遗传多样性的初步分析
有 4个特有单倍型(Hd = 0􀆰80)ꎬ 甘肃兰州居群
(LZ)除共有单倍型 H5 外还具有 3 个特有单倍型
(Hd = 0􀆰87)ꎬ 广东英德居群(YD)除共有单倍型
H6外也具有 3 个特有单倍型(Hd = 0􀆰80)(表 3ꎬ
图 3)ꎬ 即穗状狐尾藻这 3 个居群的单倍型数量最
多ꎻ 山东枣庄居群 (ZZ)同时拥有单倍型 H5 和
H6ꎬ 这是否显示穗状狐尾藻存在多个避难所和重
复扩张还需要进一步研究ꎮ 根据 Chen 等[35]对小
二仙草科植物生物地理历史的研究推测ꎬ 穗状狐尾
藻在 4􀆰9 Mya 已经分化形成ꎬ 而青藏高原开始强
烈隆升发生在约 3􀆰6 Mya[41ꎬ42]ꎬ 本研究中穗状狐
尾藻的西藏米林居群(ML)仅具有一个单倍型 H5ꎬ
揭示其应在青藏高原隆升之前已扩散到米林等青藏
高原地区ꎮ 然而ꎬ 穗状狐尾藻这一扩散的途径与机
制还不清楚ꎮ
一个物种的遗传多样性越高ꎬ 其环境适应能力
就会越强ꎬ 而且更容易向其他地区扩张[43]ꎮ 穗状
狐属藻为世界性分布ꎬ 本研究通过对来自中国大陆
地区穗状狐尾藻 8个居群的叶绿体基因组片段序列
的遗传多样性初步分析发现ꎬ 该物种具有很高的遗
传多样性(Hd = 0􀆰83)ꎬ 这可能是导致其具有广泛
分布格局(从黑龙江到广东)及极强的环境适应能
力(西藏的高原环境等)的一个原因ꎮ 同时ꎬ 细胞
学证据表明穗状狐尾藻是六倍体物种[44ꎬ45]ꎬ 相对
于二倍体来说该物种更能够成功的进行长距离扩
散[46]ꎬ 这也可能是其具有广泛分布格局的一个重
要原因ꎮ 但穗状狐尾藻可能的遗传多样性中心及第
四纪冰期是否对其多倍化物种的形成和扩张产生影
响等依然不清楚ꎮ 因此ꎬ 需更广泛的取样ꎬ 并结合
其多倍化物种的形成过程来深入探讨穗状狐尾藻的
遗传结构、 谱系地理历史及扩张途径ꎮ
参考文献:
[ 1 ]   Strand JAꎬ Weisner SEB. Morphological plastic responses
to water depth and wave exposure in an aquatic plant
(Myriophyllum spicatum)[J] . J Ecolꎬ 2001ꎬ 89(2): 166-
175.
[ 2 ]   范国兰ꎬ 李伟. 穗状狐尾藻(Myriophyllum spicatum L.)在不
同程度富营养化水体中的营养积累特点及营养分配对策
[J] . 武汉植物学研究ꎬ 2005ꎬ 23(3): 267-271.
Fan GLꎬ Li W. Response of nutrient accumulation charac ̄
teristics and nutrient strategy of Myriophyllum spicatum L.
under different eutrophication conditions [ J] . Journal of
Wuhan Botanical Researchꎬ 2005ꎬ 23(3): 267-271.
[ 3 ]   钟爱文ꎬ 曹特ꎬ 张萌ꎬ 倪乐意ꎬ 谢平. 光照和黑暗条件下苦草
(Vallisneria natans)和穗状狐尾藻 (Myriophyllum spica ̄
tum)对铵态氮的吸收[J] . 湖泊科学ꎬ 2013ꎬ 25(2): 289-
294.
Zhong AWꎬ Cao Tꎬ Zhang Mꎬ Ni LYꎬ Xie P. Uptake of am ̄
monium by Vallisneria natans and Myriophyllum spicatum
under light and dark regimes [ J] . Journal of Lake Sci ̄
enceꎬ 2013ꎬ 25(2): 289-294.
[ 4 ]   Li GXꎬ Zhang DDꎬ Li QSꎬ Chen GY. Effects of pH on iso ̄
therm modeling and cation competition for Cd(Ⅱ) and Cu
(Ⅱ) biosorption on Myriophyllum spicatum from aqueous
solutions[J] . Environ Earth Sciꎬ 2014ꎬ 72 (11): 4237 -
4247.
[ 5 ]   Heine Sꎬ Schmitt Wꎬ Görlitz Gꎬ Schäffer Aꎬ Preuss TG.
Effects of light and temperature fluctuations on the growth
of Myriophyllum spicatum in toxicity tests—a model ̄based
analysis[ J ] . Environ Sci Pollut Resꎬ 2014ꎬ 21 ( 16 ):
9644-9654.
[ 6 ]   袁桂香ꎬ 符辉ꎬ 钟家有ꎬ 倪乐意ꎬ 朱天顺ꎬ 李威ꎬ 宋鑫. 铵
胁迫对狐尾藻(Myriophyllum spicatum)和金鱼藻(Cerato ̄
phyllum demersum)生物量分配和形态的影响[J] . 湖泊科
学ꎬ 2013ꎬ 25(5): 729-734.
Yuan GXꎬ Fu Hꎬ Zhong JYꎬ Ni LYꎬ Zhu TSꎬ Li Wꎬ Song X.
Changes in biomass allocation and morphology of Myrio ̄
phyllum spicatum and Ceratophyllum demersum under the
ammonium stress[J] . Journal of Lake Scienceꎬ 2013ꎬ 25
(5): 729-734.
[ 7 ]   Heine Sꎬ Schmitt Wꎬ Schaffer Aꎬ Gorlitz Gꎬ Buresová Hꎬ
Arts Gꎬ Preuss TG. Mechanistic modelling of toxicokinetic
processes within Myriophyllum spicatum [ J ] . Chemo ̄
sphereꎬ 2015ꎬ 120: 292-298.
[ 8 ]   LaRue EAꎬ Grimm Dꎬ Thum RA. Laboratory crosses and
genetic analysis of natural populations demonstrate sexual
viability of invasive hybrid watermilfoils (Myriophyllum spi ̄
catum ×M. sibiricum)[J] . Aquat Botꎬ 2013ꎬ 109: 49-53.
[ 9 ]   Weyl PSꎬ Coetzee JA. The invasion status of Myriophyllum
spicatum L. in southern Africa[J] . Managementꎬ 2014ꎬ 5
(1): 31-37.
[10]   Borrowman KRꎬ Sager EPSꎬ Thum RA. Distribution of bio ̄
types and hybrids of Myriophyllum spicatum and associa ̄
ted Euhrychiopsis lecontei in lakes of Central Ontarioꎬ
Canada[ J] . Lake Reserv Manageꎬ 2014ꎬ 30 ( 1): 94 -
104.
[11]   Moody MLꎬ Les DH. Evidence of hybridity in invasive wa ̄
termilfoil (Myriophyllum) populations[ J] . Proc Natl Acad
Sci USAꎬ 2002ꎬ 99(23): 14867-14871.
[12]   Clegg MTꎬ Gaut BSꎬ Learn GHꎬ Morton BR. Rates and
patterns of chloroplast DNA evolution[J] . Proc Natl Acad
Sci USAꎬ 1994ꎬ 91(15): 6795-6801.
[13]   李丹丹ꎬ 郭水良ꎬ 于晶ꎬ 李莎ꎬ 曹同. 基于 4 个叶绿体基因
411 植 物 科 学 学 报 第 34卷 
识别蓑藓属 (Macromitrium)植物的可行性研究[J] . 植物科
学学报ꎬ 2013ꎬ 31(1): 23-33.
Li DDꎬ Guo SLꎬ Yu Jꎬ Li Sꎬ Cao T. Feasibility study on the
identification of genus Macromitrium based on four chloro ̄
plast genes [ J] . Plant Sceince Journalꎬ 2013ꎬ 31 ( 1):
23-33.
[14]   赵晶华ꎬ 贾倩ꎬ 金洪涛ꎬ 盛彦敏. 木兰叶绿体 atpB 和 rbcL
基因的系统进化分析[J] . 长春师范大学学报: 自然科学版ꎬ
2014ꎬ 33 (6): 71-74.
Zhao JHꎬ Jia Qꎬ Jin HTꎬ Sheng YM. Evolution analysis of
the chloroplast genes atpB and rbcL of magnolia [ J ] .
Journal of Changchun Normal Universityꎬ 2014ꎬ 33 (6):
71-74.
[15]   宋敏舒ꎬ 乐霁培ꎬ 孙航ꎬ 李志敏. 横断山地区海仙报春的谱
系地理学研究[ J] . 植物分类与资源学报ꎬ 2011ꎬ 33(1):
91-100.
Song MSꎬ Yue JPꎬ Sun Hꎬ Li ZM. Phylogeographical
study on Primula poissonii (Primulaceae) from Hengduan
Mountains[ J] . Plant Diversity and Resourcesꎬ 2011ꎬ 33
(1): 91-100.
[16]   Chen JMꎬ Du ZYꎬ Sun SSꎬ Gituru RWꎬ Wang QF. Chloro ̄
plast DNA phylogeography reveals repeated range expan ̄
sion in a widespread aquatic herb Hippuris vulgaris in the
Qinghai ̄Tibetan Plateau and adjacent areas [ J ] . PloS
ONEꎬ 2013ꎬ 8(4): e60948.
[17]   Zhang Qꎬ Chiang TYꎬ George Mꎬ Liu JQꎬ Abbott RJ. Phy ̄
logeography of the Qinghai ̄Tibetan Plateau endemic Juni ̄
perus przewalskii (Cupressaceae) inferred from chloro ̄
plast DNA sequence variation [ J] . Mol Ecolꎬ 2005ꎬ 14
(11): 3513-3524.
[18]   陈家瑞. 中国植物志: 第 53 卷[M] . 北京: 科学出版社ꎬ
2000: 136.
Chen JR. Flora Republicae Popularis Sinicae: Vol. 53[M] .
Beijing: Science Pressꎬ 2000: 136.
[19]   Doyle JJ. A rapid DNA isolation procedure for small quan ̄
tities of fresh leaf tissue[ J] . Phytochem Bullꎬ 1987ꎬ 19:
11-15.
[20]   Johnson LAꎬ Soltis DE. matK DNA sequences and phylo ̄
genetic reconstruction in Saxifragaceae s. str. [ J] . Syst
Botꎬ 1994: 143-156.
[21]   Taberlet Pꎬ Gielly Lꎬ Pautou Gꎬ Bouvet J. Universal pri ̄
mers for amplification of three non ̄coding regions of chlo ̄
roplast DNA [ J] . Plant Mol Biolꎬ 1991ꎬ 17 ( 5): 1105 -
1109.
[22]   Shaw Jꎬ Lickey EBꎬ Schilling EEꎬ Small RL. Comparison
of whole chloroplast genome sequences to choose non ̄
coding regions for phylogenetic studies in angiosperms:
the tortoise and the hare Ⅲ[J] . Am J Botꎬ 2007ꎬ 94(3):
275-288.
[23]   Chenna Rꎬ Sugawara Hꎬ Koike Tꎬ Lopez Rꎬ Gibson TJꎬ
Higgins DGꎬ Thompson JD. Multiple sequence alignment
with the Clustal series of programs [ J ] . Nucleic Acids
Resꎬ 2003ꎬ 31(13): 3497-3500.
[24]   Rozas Jꎬ Rozas R. DnaSPꎬ DNA sequence polymor ̄
phism: an interactive program for estimating population
genetics parameters from DNA sequence data[J] . Com ̄
put Applic Biosciꎬ 1995ꎬ 11(6): 621-625.
[25]   Tajima F. Statistical method for testing the neutral mutation
hypothesis by DNA polymorphism [ J] . Geneticsꎬ 1989ꎬ
123(3): 585-595.
[26]   Fu YXꎬ Li WH. Statistical tests of neutrality of mutations
[J] . Geneticsꎬ 1993ꎬ 133(3): 693-709.
[27]   Fu YX. Statistical tests of neutrality of mutations against
population growthꎬ hitchhiking and background selection
[J] . Geneticsꎬ 1997ꎬ 147(2): 915-925.
[28]   Excoffier Lꎬ Laval Gꎬ Schneider S. Arlequin ( version
3􀆰0): an integrated software package for population ge ̄
netics data analysis[J] . Evol Bioinform Onlineꎬ 2007ꎬ 1:
47-50.
[29]   Peakall Rꎬ Smouse PE. GenAlEx 6􀆰5: genetic analysis in
Excel. Population genetic software for teaching and re ̄
search— an update[ J] . Bioinformaticsꎬ 2012ꎬ 28(19):
2537-2539.
[30]   Swofford DL. PAUP∗: Phylogenetic Analysis using Parsi ̄
mony (∗and other methods) [M] . Version 4􀆰0b10. Sun ̄
derland: Sinauer Associatesꎬ 2003.
[31]   Bandelt HJꎬ Forster Pꎬ Röhl A. Median ̄joining networks
for inferring intraspecific phylogenies [ J] . Mol Biol Evolꎬ
1999ꎬ 16(1): 37-48.
[32]   Barrett SCHꎬ Eckert CGꎬ Husband BC. Evolutionary
processes in aquatic plant populations [ J] . Aquat Botꎬ
1993ꎬ 44(2): 105-145.
[33]   Sculthorpe CD. Biology of Aquatic Vascular Plants [M] .
London: Edward Arnoldꎬ 1967.
[34]   Santamaría L. Why are most aquatic plants widely distribu ̄
ted? Dispersalꎬ clonal growth and small ̄scale heterogene ̄
ity in a stressful environment[ J] . Acta Oecolꎬ 2002ꎬ 23
(3): 137-154.
[35]   Chen LYꎬ Zhao SYꎬ Mao KSꎬ Les DHꎬ Wang QFꎬ Moody
ML. Historical biogeography of Haloragaceae: An out ̄of ̄
Australia hypothesis with multiple intercontinental disper ̄
sals[J] . Mol Phylogenet Evolꎬ 2014ꎬ 78: 87-95.
[36]   Laushman RH. Population genetics of hydrophilous angio ̄
sperms[J] . Aquat Botꎬ 1993ꎬ 44(2): 147-158.
[37]   Djebrouni M. Variabilité morphologiqueꎬ caryologique et
enzymatique chez quelques populations de Phragmites
australis (Cav.) Trin. Ex Steud[ J] . Folia Geobotanica et
Phytotaxonomicaꎬ 1992ꎬ 27(1): 49-59.
[38]   Hollingsworth PMꎬ Preston CDꎬ Gornall RJ. Genetic varia ̄
bility in two hydrophilous species of Potamogetonꎬ P. pec ̄
511  第 1期                        吕天锋等: 中国穗状狐尾藻遗传多样性的初步分析
tinatus and P. filiformis ( Potamogetonaceae) [ J] . Plant
Syst Evolꎬ 1996ꎬ 202(3 / 4): 233-254.
[39]   Koga Kꎬ Kadono Yꎬ Setoguchi H. The genetic structure of
populations of the vulnerable aquatic macrophyte Ranun ̄
culus nipponicus ( Ranunculaceae ) [ J ] . J Plant Resꎬ
2007ꎬ 120(2): 167-174.
[40]   Gao Lꎬ Ge Sꎬ Hong D. Low levels of genetic diversity
within populations and high differentiation among popula ̄
tions of a wild riceꎬ Oryza granulata Nees et Arn. ex Watt.ꎬ
from China[J] . Int J Plant Sciꎬ 2000ꎬ 161(4): 691-697.
[41]   Zhang MLꎬ Fritsch PW. Evolutionary response of Caraga ̄
na ( Fabaceae ) to Qinghai ̄Tibetan Plateau uplift and
Asian interior aridification[ J] . Plant Syst Evolꎬ 2010ꎬ 288
(3-4): 191-199.
[42]   Li Jꎬ Fang X. Uplift of the Tibetan Plateau and environmen ̄
tal changes[J] . Chinese Sci Bullꎬ 1999ꎬ 44(23): 2117-
2124.
[43]   Hamrick JLꎬ Godt MJW. Allozyme diversity in plant spe ̄
cies[M] / / Brown AHDꎬ Clegg MTꎬ Kahler ALꎬ Weir BSꎬ
eds. Plant Population Geneticsꎬ Breedingꎬ and Genetic
Resources. Sunderlandꎬ MA: Sinauer Associates Inc.ꎬ
1990: 43-63.
[44]   Aiken SG. Counts on Haloragaceae[J] . Taxonꎬ 1978ꎬ 27:
519-535.
[45]   Aiken SG. A conspectus of Myriophyllum (Haloragaceae)
in North America[J] . Brittoniaꎬ 1981ꎬ 33(1): 57-69.
[46]   Linder HPꎬ Barker NP. Does polyploidy facilitate long ̄dis ̄
tance dispersal? [ J] . Ann Botꎬ 2014ꎬ 113 ( 7): 1175 -
1183.
(责任编辑: 刘艳玲)
611 植 物 科 学 学 报 第 34卷