中国扁桃亚属四种野生扁桃的系统发育与物种分化 植物生态学报    2021, 45 (9): 987-995.   DOI: 10.17521/cjpe.2020.0366

图4 扁桃亚属4种野生植物基于种群水平的前3个坐标的主坐标分析(PCoA)。

正文中引用本图/表的段落

4种扁桃的核序列ITS1-ITS4多重比对后的总长为634 bp, 包括27个核苷酸变异位点, 共定义了28个单倍型(图2)。单倍型网络表明, 28个单倍型聚为4支, 分别对应新疆野扁桃、蒙古扁桃、长柄扁桃和西康扁桃, 种间无共享单倍型(图2)。4种植物的遗传距离分析显示: 种内和种间的遗传分化距离分别为0-0.009和0.014-0.066, 种内和种间的距离间隙为0.005, 表明4个物种种间存在显著的遗传分化(图3)。4种扁桃的PCoA显示, 新疆野扁桃、蒙古扁桃、长柄扁桃和西康扁桃总体上各聚为一组, 其中新疆野扁桃和蒙古扁桃及西康扁桃之间的遗传距离相对较小, 而三者与长柄扁桃的遗传距离相对较大(图4)。

目前, 对于李亚科的分属存在两种划分观点, 一类为广义李属(Benthan & Hooker, 1865), 即: 桃、李、梅、杏等核果统归于李属, 再下设7个组; 另一类为狭义李属(俞德浚, 1979), 采用小属的划分方法, 共7个属: 桃属、杏属(Armeniaca)、樱属(Cerasus)、桂樱属(Laurocerasus)、稠李属(Padus)和李属(Prunus), 每个属下又分亚属(如桃属分扁桃亚属和桃亚属(subgen. Persica)), 该方法同时也是《中国植物志》的分类方法。李亚科内所有植物均为单心皮核果, 关于其属下亚属或组, 以及种的划分, 几百年来均无统一结果。本研究基于ITS序列, 最大似然树、贝叶斯树、单倍型网络及PCoA均表明4种扁桃的单倍型呈4个亚支, 种间无共享单倍型(图2, 图4, 图5)。形态上, 4种扁桃的枝、叶、果实和种子等也均存在明显差异(赵一之, 1995; 毛祖梅等, 1996)。

前人利用叶绿体序列和单拷贝核基因序列建树得出, 新疆野扁桃、长柄扁桃与欧洲李、黑刺李(Prunus spinosa)和稠李(Padus avium)等李属、稠李属植物, 蒙古扁桃和西康扁桃与桃(Prunus persica)、紫叶桃(Prunus persica ‘Zi Ye Tao')、山桃等桃属植物分别具有较近的亲缘关系(Yazbek & Oh, 2010, 2013; 董山平, 2015)。本研究表明4种野生扁桃与两个外类群山桃和稠李间的亲缘及聚类关系均较远(图1, 图3, 图4)。4种扁桃果实成熟后开裂无汁液与山桃和稠李果实成熟后不开裂多汁也具有明显的形态差异(俞德浚, 1986)。4种扁桃的PCoA和系统发育分析均表明, 新疆野扁桃、蒙古扁桃、西康扁桃之间亲缘关系相对较近, 而与长柄扁桃的亲缘关系均相对较远(图4, 图5)。基于核标记和形态特征对扁桃亚属系统发育的研究也揭示了长柄扁桃和其他3种扁桃的亲缘关系相对较远(Delplancke et al., 2016)。从染色体核型上, 新疆野扁桃、蒙古扁桃、西康扁桃均为二倍体16条染色体, 长柄扁桃为十二倍体96条染色体(董山平, 2015)。根据植物由二倍体向多倍体演化的规律, 分布区为4种植物分布最东界的长柄扁桃很有可能为4种扁桃中最年轻的物种(李雅轩, 2001; Yazbek & Oh, 2013)。

本文的其它图/表

• 图1 4种扁桃亚属野生植物的野外采样点。BEC, 布尔津县; BTE, 大青山2; BTQ, 大青山1; DBL, 迭部县1; DBZ, 迭部县2; DMB, 达茂旗3; DMD, 达茂旗1; DMN, 达茂旗2; GQD, 固阳1; GQX, 固阳3; GX, 固阳2; HBK, 哈巴河县; MZS, 马鬃山; SPQ, 松潘县; TBB, 塔城2; TBT, 塔城1; TLB1, 托里县1; TLB2, 托里县2; WD, 乌拉山; WG, 乌盖苏木乡; WLB, 大坝口; WLJ, 乌加禾镇; WSH, 苏海图; WSM, 大麦力沟; YCQ, 龙首山; YLD1, 榆林1; YLD2, 榆林2; YMB1, 裕民县1; YMB2, 裕民县2; YMF, 阴山2; YMM, 阴山1; ZQK, 图克木嘎查; ZQL, 贺兰山; ZQX, 舟曲县; ZTL, 特莫乌拉嘎查; ZTM, 苏力图嘎查; ZYS, 祁连山。
  
• 图2 扁桃亚属4种野生植物28个单倍型(R1-R28)的地理分布及其单倍型网络。图中种群采样点编码与图1一致, 饼图表示各种群的单倍型频率。A, 单倍型网络, 图中圆圈大小与单倍型频率成正比, 节点间的分支长度大致与单倍型的突变数成正比, 相应分支附近附有步长; 稠李和山桃作为外类群。
  
• 图3 基于ITS序列计算的扁桃亚属4种野生植物种间和种内的遗传距离。
  
• 图5 扁桃亚属4种野生植物基于ITS序列的单倍型系统发育树。A, 最大似然(ML)树。分支点上方的数字为大于等于70的自展支持率。B, 贝叶斯树。分支节点右侧的数字表示大于0.70的后验概率值, 两系统发育树右侧的中括号表示为4个物种相应的亚支。
  
• 图6 扁桃亚属4种野生植物生态位相似性与驱动因子贡献率。A-F, 生态位可视化。绿色表示第一个种的生态位空间, 红色表示第二个种的生态位空间, 蓝色表示重叠空间。D为生态位重叠分数。G, 驱动因子主成分分析。箭头表示相关的方向(同一方向表示高度相关)。Bio1, 年平均气温; Bio4, 温度季节性变化标准差; Bio5, 最暖月最高气温; Bio6, 最冷月最低气温; Bio15, 降水量的季节性变异系数; Bio17, 最干季降水量; Bio18, 最暖季降水量。环境因子后括号里的数字为贡献率排名。