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

图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, 祁连山。

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

本研究共采集扁桃亚属4种野生植物的37个种群, 每个种群采集8-12份样品。采样的扁桃株间至少相距10 m, 均采集新鲜、幼嫩、无病斑的叶片, 用硅胶干燥保存于自封袋, 带回实验室置于-20 ℃冷藏。利用Excel随机抽取公式随机挑选每个种群的4个个体, 共148份样品开展4种扁桃植物的系统发育和物种分化研究(图1)。

前人利用叶绿体序列和单拷贝核基因序列建树得出, 新疆野扁桃、长柄扁桃与欧洲李、黑刺李(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)。

本文的其它图/表

• 图2 扁桃亚属4种野生植物28个单倍型(R1-R28)的地理分布及其单倍型网络。图中种群采样点编码与图1一致, 饼图表示各种群的单倍型频率。A, 单倍型网络, 图中圆圈大小与单倍型频率成正比, 节点间的分支长度大致与单倍型的突变数成正比, 相应分支附近附有步长; 稠李和山桃作为外类群。
  
• 图3 基于ITS序列计算的扁桃亚属4种野生植物种间和种内的遗传距离。
  
• 图4 扁桃亚属4种野生植物基于种群水平的前3个坐标的主坐标分析(PCoA)。
  
• 图5 扁桃亚属4种野生植物基于ITS序列的单倍型系统发育树。A, 最大似然(ML)树。分支点上方的数字为大于等于70的自展支持率。B, 贝叶斯树。分支节点右侧的数字表示大于0.70的后验概率值, 两系统发育树右侧的中括号表示为4个物种相应的亚支。
  
• 图6 扁桃亚属4种野生植物生态位相似性与驱动因子贡献率。A-F, 生态位可视化。绿色表示第一个种的生态位空间, 红色表示第二个种的生态位空间, 蓝色表示重叠空间。D为生态位重叠分数。G, 驱动因子主成分分析。箭头表示相关的方向(同一方向表示高度相关)。Bio1, 年平均气温; Bio4, 温度季节性变化标准差; Bio5, 最暖月最高气温; Bio6, 最冷月最低气温; Bio15, 降水量的季节性变异系数; Bio17, 最干季降水量; Bio18, 最暖季降水量。环境因子后括号里的数字为贡献率排名。