MORPHOLOGICAL ANALYSIS USED BY GEOMETRIC MORPHOMETRICS COMBINED WITH MICRO CT AMONG GYMNOCYPRIS ECKLONI IN TWO DRAINAGE (TELEOSTEI: CYPRINIDAE)
-
摘要: 研究采用现代几何形态测量法并结合Micro CT扫描技术对两水系(黄河上游与柴达木水系格尔木河)花斑裸鲤(Gymnocypris eckloni Herzenstein)的形态进行定量分析, 讨论两水系花斑裸鲤是否存在先前研究未能发现的形态差异。结果显示: 花斑裸鲤的两个地理种群在基于其头部形态轮廓二维坐标的主成分分析和判别分析的结果均聚在一起, 这与之前文献记载的两水系花斑裸鲤的形态学特征相一致, 并且对两水系花斑裸鲤进行的典型变量分析也同样支持这一结果; 对花斑裸鲤两个地理种群背鳍支鳍骨插入椎骨间的相对位置这一形态特征进行卡方检验, 结果表明在骨骼这一更具分类学意义的形态特征上表现出显著性差异。通过对花斑裸鲤两个地理种群的几何形态与Micro CT扫描, 准确定量地分析了两水系花斑裸鲤形态上的异同, 为正确评价两水系花斑裸鲤真实的演化关系奠定了形态学基础。
-
关键词:
- 花斑裸鲤 /
- 现代几何形态测量法 /
- Micro CT扫描技术 /
- 定量分析 /
- 形态鉴定
Abstract: Gymnocypris eckloni is widely distributed in the Yellow River and Qaidam River in the northeast of the Qinghai-Tibet Plateau, China. Inconsistence between morphology and phylogenetics was observed in two different drainage G. eckloni. On this basis, we speculated that the first was the morphological convergence caused by the similar ecological pressure of the two water systems, and the second was the traditional measurement method, which makes some subtle morphological features undiscovered. Hence, the combination of geometric morphometrics and Micro CT revealed the phenotypic difference among G. eckloni of two water systems. The structure of G. eckloni’s head was almost identical between two drainages, demonstrating a significant difference in their bone structure. The divergence between head morphology and molecular data were presumably attributed to one reason. It may be that the morphological convergence caused by the similar ecological environment of the two river systems misled the taxonomic status of Gymnocypris eckloni of two water systems. The current study has laid a morphology foundation for correctly evaluating the real relationship between Gymnocypris eckloni of two water systems, and suggest that bone is an accurate standard for taxonomic classification in schizothoracin fishes. More molecular markers are required for further research to uncover the evolutionary pattern of G. eckloni from different geographic distribution. -
花斑裸鲤(Gymnocypris eckloni Herzenstein), 隶属鲤形目(Cypriniformes), 鲤科(Cyprinidae), 裂腹鱼亚科(Schizothoracinae), 主要分布在我国柴达木水系格尔木河和黄河水系上游干支流[1]。花斑裸鲤物种最早是由俄国鱼类学家Herzenstein于1891年提出, 他运用尽可能多的标本研究花斑裸鲤的形态性状, 指出花斑裸鲤存在两个地理种群, 主要分布于黄河水系上游及柴达木水系格尔木河, 并将格尔木河的花斑裸鲤定为模式标本[2]。在随后的近一百年里, 我国鱼类学家相继对花斑裸鲤的形态和分类进行过大量研究。武云飞等[1]对两水系花斑裸鲤的形态进行大量统计分析, 结果显示除扎陵湖与鄂陵湖种群的鳃耙数偏高, 呈现湖泊种所具有的地区性差异外, 其余特征, 包括口裂、背鳍距尾鳍基部距离、下颌形态等在两种群间均没有差异。此外, 陈毅峰和曹文宣[3]经过大量形态学研究也得到了类似结果, 认为两水系花斑裸鲤形态相似, 是同一物种的不同地理种群。同样, 赵铁桥[4]在对格尔木河、黄河水系上游与河西水系不同地理种群花斑裸鲤的形态学研究中提出, 花斑裸鲤具有中间形态, 其形态上的差异难以构成种级分化标准。以上学者基于形态学研究均认为黄河与格尔木河的花斑裸鲤是同一物种的不同地理种群。
然而, 自从20世纪80年代随着分子生物学的发展, 分子水平上对两水系裸鲤系统关系的研究却呈现出与形态学研究完全不一致的结果。赵凯等[5]利用线粒体标记Cyt b基因对黄河水系、柴达木水系格尔木河花斑裸鲤的系统发育关系进行了研究, 发现两水系花斑裸鲤经过地理隔离已经发生了独立进化。随后, Zhao等[6]在增加了祁连裸鲤的样品后, 发现花斑裸鲤的两个地理种群没有形成单系群, 而显示为相互独立起源的并系, 再次肯定了上述结果的有效性。这一研究结果对花斑裸鲤的分类学及其两个地理种群的分类地位提出了不同的看法。
长期以来, 对花斑裸鲤的形态学研究多采用传统形态测量方法, 这一方法容易受主观因素的影响, 也容易忽视一些肉眼难以观察到的性状[7, 8], 同时传统的形态分类标准, 如口裂、鳃耙等多为环境适应性特征, 易受环境影响而发生改变, 不能客观反映物种的形态特征。随着科学技术的发展, 现代几何形态测量法(Modern geometric morphometrics) 和Micro CT扫描技术(Micro computed tomography, 微计算机断层扫描技术) 被广泛应用于物种形态鉴定。几何形态测量法是于20世纪90年代Bookstein[9]和Rohlf[10]以数学为基础所创, 将统计学方法应用于形态分析, 通过标志点(Landmark) 与滑动半标志点(Semi-landmark) 相互补充, 对整体轮廓进行描绘[11], 并最终通过多元统计分析对物种形态进行相似性比较和差异分析[12－14], 进而定量判定形态特征的异同。经过近30年的发展, 几何形态测量法已渗透到多个领域中[15－17], 且在裂腹鱼类的形态鉴定中也已应用成熟[18－20]。Micro CT扫描是基于X射线对活体小动物进行断层扫描的一项技术, 是利用不同组织对X射线吸收和透过率不同的原理实现的[21]。Micro CT扫描技术因其超高的分辨率更容易发现骨骼结构上的差异, 现也广泛应用于小鼠骨密度研究、硬骨鱼类脊椎研究等多种小动物的骨骼形态研究方面[22－24]。
因此, 为清楚阐释两水系裸鲤是否因为传统测量方法的局限性, 使其分子与形态学研究不一致这一问题, 本研究运用新型的几何形态测量法并结合Micro CT扫描技术, 对两水系花斑裸鲤的头部形态及骨骼特征进行了定量分析, 这对正确评价两水系花斑裸鲤真实的系统发育关系具有重要意义。
1. 材料与方法
1.1 实验材料
本研究所用花斑裸鲤两个地理种群的样品分别采集于青藏高原东北部的黄河水系上游河段和柴达木水系格尔木河(图 1), 共采集116个样本(表 1), 同时作为本次实验对照组的青海湖裸鲤(Gymnocypris przewalskii przewalskii) 样品均采集于青海湖水系(图 1), 共19个样本(表 1)。采集时间分别为2010年、2012年和2014年。采集时使用流刺网、地笼捕捞, 并保存于95%酒精中, 部分保存于福尔马林溶液中, 现均存于中国科学院西北高原生物研究所高原鱼类进化与功能基因组学实验室和中国科学院西北高原生物研究所标本馆内。
![]()
图 1 采样地点示意图图中圆圈代表采样地点Figure 1. Map of the head-water area of the Yellow River, Golmud River, and Qinghai Lake showing the sampling sitesThe circle represents the sampling site表 1 样品的采集水系、地点及数量Table 1. Populations, numbers, river systems and sampling locations used in this study编号No. 物种名Species 水系Drainage 采集地Sampling site 数量Number of specimens 1 花斑裸鲤 格尔木河 清水河 19 2 花斑裸鲤 格尔木河 泽格斯特河 23 3 花斑裸鲤 格尔木河 郭勒木德镇傲色河 10 4 花斑裸鲤 格尔木河 温泉水库 10 5 花斑裸鲤 黄河 扎陵湖 6 6 花斑裸鲤 黄河 鄂陵湖 13 7 花斑裸鲤 黄河 共和曲沟 5 8 花斑裸鲤 黄河 贵德 30 9 青海湖裸鲤 青海湖水系 哈尔盖河入湖口 19 总计Total 135 1.2 图片信息获取
选取头部侧面轮廓和下颌轮廓对黄河与格尔木河花斑裸鲤(N=116) 进行形态比较。此外, 补充了19个青海湖裸鲤样本作为外类群进行对照。具体方法如下。
所有样品均使用锡恩帝HX700数码相机拍摄, 每个样品采集两张照片, 分别为鱼头部左侧和腹部图片, 拍摄时用镊子或大头针固定鱼的形态使其保持自然状态[19]。为减少系统误差, 使用翻拍架固定相机位置, 布置蓝色背景布, 放置刻度尺, 统一样品摆放位置, 以确保样品与采集面保持平行, 与相机高度距离相等。为减少拍摄时的人为误差与偶然误差, 所有拍摄均由一人完成。
1.3 标志点与半标志点选取
标志点需选取既有同源性又能准确反映物种形态信息的位点, 根据类似的研究[19, 20], 并结合花斑裸鲤实际形态特点, 本次实验选取了花斑裸鲤头部侧面5个标志点: ①头骨背缘末端, ②吻上唇顶点, ③吻下唇顶点, ④口裂末端, ⑤腹中线上腮盖条膜后缘。为使其能更加精确地描绘头部轮廓, 结合文献采用滑动技术的方法选取了17个半标志点。半标志点的确定是在IMP软件包MakeFan8 (canisius.edu/~sheets/IMP%208.htm) 软件中实现的, 首先确定同源性位点, 然后以两者之间的连线做11条距离相等的垂线, 其垂线与头部轮廓的交点即为半标志点(图 2a)。花斑裸鲤下颌轮廓选取了2个标志点: ①下颌骨上端, ②下颌骨下端。下颌骨轮廓半标志点的确定同样在MakeFan8中实现, 以确定的两个标志点作连线为直径画半圆, 再以一组射线将其分成11份, 则射线与下颌轮廓的交点即为半标志点(图 2b)。随后利用TpsDig 2.04软件记录标志点与半标志点。
![]()
图 2 花斑裸鲤标志点、半标志点的选取a. 头部侧面轮廓特征点; b. 下颌轮廓特征点Figure 2. Landmarks and semi-landmarks of G. ecklonia. Landmarks in the leftside of the head; b. Landmarks of the Mandibular contour1.4 图像数字化处理与统计分析
采用TPS系列软件对图片进行数字化处理[25]。利用TpsUtil (tpsutil.updatestar.com) 修改样本的ID, 把本研究所涉及到的照片数据“矢量化”, 最终生成扩展名为“.NTS”的二维坐标矩阵文件, 用tpsDig 2.17 (life.bio.sunysb.edu/morph) 统一所有照片的比例尺以及按上述标志点对样本头部侧面和下颌进行标定。随后将生成的“.NTS”文件导入MorphoJ[26](flywings.org.uk/MorphoJ_page.htm) 软件中, 进行New Procrustes fit消除非形变因素对数据分析的影响, 即对样本进行普氏印叠(Generalized Procrustes Analysis, GPA)。GPA实际上是对数据标准化的过程, 即通过对每个样品的标志点进行平移、缩放、旋转的过程来消除非形状因素产生的误差, 使样本的几何形态得到精准的客观复原[9－13]。然后运用MorphoJ对普氏坐标数据进行主成分分析(Principal components analysis, PCA), 用以分析不同水系花斑裸鲤几何形态上的差异; 之后将二维坐标数据导入到PAST3.0[27]中进行辨别分析(Discriminant analysis, DLA), 根据其不同地理种群二维坐标的特征值判断两水系花斑裸鲤的种群归属问题; 进行典型变量分析(Canonical variable analysis, CVA), 讨论花斑裸鲤两个地理种群间的差异, 并计算马氏距离(Mahalanobis distance) 和普氏距离(Procrustes distance), 来定量分析两种群间的形态差异量。
1.5 Micro CT扫描
本次实验使用QuantumGX Micro CT扫描系统, 分别选取黄河水系花斑裸鲤21条, 格尔木河花斑裸鲤22条进行骨骼扫描, 分析两水系花斑裸鲤在骨骼方面是否存在未被发现的细微特征。参数设置: 电压70 kV, 电流114 μA(该电流为此电压下的最大电流), 视场86 mm, X-射线渗透Cu 0.1 mm, 扫描模式采用High Resolution 4min, 其余保持默认参数。
2. 结果
2.1 两水系花斑裸鲤头部形状的主成分分析
几何形态学分析显示, 花斑裸鲤的两个地理种群在主成分分析上并没有形成两个形态组, 而是聚在一起 (图 3)。头部轮廓的前两个PCA共解释了85.15%的变异量, 其中PC1的贡献率为70.42%, PC2的贡献率为14.73%。在PCA二维分布散点图中, 黄河与柴达木水系格尔木河的花斑裸鲤聚在一起, 有很大程度的重叠, 并没有明显的形态差异。沿PC1轴的正方向, 两水系花斑裸鲤最明显的变化表现在1、2、3号特征点上, 即口型位置发生了变化, 逐渐表现为口亚下位的特点, 符合花斑裸鲤的口型特征。
![]()
图 3 黄河与格尔木河水系花斑裸鲤头部侧面轮廓(a) 与下颌轮廓(b) 的主成分分析Figure 3. Principal component analysis of the head profile (a) and mandibular profile (b) of the G. eckloni in Yellow River and Golmud River下颌轮廓的前两个PCA共解释了87.9%的变异量, 其中PC1的贡献率为80.98%, PC2的贡献率为6.92%。在下颌轮廓的PCA散点图中, 两水系花斑裸鲤与青海湖裸鲤都有一定程度的重叠, 彼此之间均没有明显差异。下颌形态的网格图同样显示两水系花斑裸鲤在下颌形态上没有差异, 均表现为吻钝圆, 下颌弧形。
2.2 两水系花斑裸鲤头部形态的判别分析
除通过PCA分析计算研究样本几何形态上的差异外, 本研究采用判别分析探讨两水系花斑裸鲤在形态上是否属于同一类别(图 4和图 5), 结果显示: 在侧面轮廓的判别分析中, Axis1和Axis2分别占总体变异的55.5%和44.5%, 判别分析的结果显示黄河水系与柴达木水系格尔木河的花斑裸鲤在侧面轮廓有很大程度的重叠, 属于同一类别, 与青海湖裸鲤存在明显差异(图 4); 图 5为下颌骨形态的判别分析结果, 其中Axis1和Axis2分别代表总变异的61.24%和38.76%, 花斑裸鲤2个地理种群与青海湖裸鲤之间均有不同程度的重叠, 三者在下颌形态轮廓特征上均无明显差异。判别分析结果与PCA结果基本一致。
![]()
图 4 花斑裸鲤头部侧面轮廓的判别分析Figure 4. The scatter diagram based on the LDA for the cephalic side contour of the G. eckloni![]()
图 5 花斑裸鲤下颌骨轮廓的判别分析Figure 5. The mandible morphological variations of the G. eckloni based on LDA2.3 两水系花斑裸鲤头部形态的典型变量分析
为量化花斑裸鲤两个地理种群间的差异, 本文对两水系花斑裸鲤分别进行了头部形态的典型变量分析(图 6), 并计算马氏距离与普氏距离(表 2—表 5)。头部侧面轮廓的典型变量分析中CV1和CV2分别代表总体变异的59.42%和40.58%, 在CV1轴上两水系花斑裸鲤没有差异, 沿CV1轴正方向表现出口亚下位的特点, 与PCA分析一致。在CV2轴上两水系花斑裸鲤有部分重叠, 其差异主要表现在头部长宽比上, 沿正值方向主要表现为头部形状略圆一些, 正值处的特征在黄河水系花斑裸鲤种群中更明显一些。
表 2 花斑裸鲤不同种群头部侧面轮廓的马氏距离比较Table 2. Comparison of Mahalanobis distance of head profile of the G. eckloni in Yellow River and Golmud River表 5 花斑裸鲤两个地理种群下颌轮廓的普氏距离比较Table 5. Procrustes distance comparison of mandibular contour between two geographic populations of G. eckloni花斑裸鲤(格尔木河水系) G. eckloni in Golmud River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 0.0748* 青海湖裸鲤G. przewalskii 0.0899* 0.0216 ![]()
图 6 黄河与格尔木河水系花斑裸鲤头部侧面轮廓(a) 与下颌轮廓(b) 的典型变量分析Figure 6. Canonical variable analysis of the head profile (a) and mandibular profile (b) of the G. eckloni in Yellow River and Golmud River下颌形态的典型变量分析中CV1代表总体变异量的75.74%, CV2代表总体变异量的24.26%。CV1展示的差异形态主要表现在下颌骨的形状上, CV1正值处格尔木水系花斑裸鲤的特征比较明显。两水系花斑裸鲤在CV2轴上存在重叠, 没有明显差异。
此外通过计算花斑裸鲤两个地理种群与青海湖裸鲤头部形态的普氏距离与马氏距离, 得到了与上述相似的结果, 即三者之间没有明显的形态变异。特别是在头部侧面轮廓计算中相较于青海湖裸鲤, 花斑裸鲤的两个地理种群间的变异较小, 无明显形态特征。
表 3 花斑裸鲤不同种群头部侧面轮廓的普氏距离比较Table 3. Comparison of Procrustes distance of head profile of the G. eckloni in Yellow River and Golmud River花斑裸鲤(格尔木河水系) G. eckloni in Golmud River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 0.0329* 青海湖裸鲤G. przewalskii 0.0739* 0.0852* 表 4 花斑裸鲤两个地理种群下颌轮廓的马氏距离比较Table 4. Mahalanobis distance of mandibular contour between two geographic populations of G. eckloni花斑裸鲤(格尔木河水系) G. eckloni in Golmud River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 2.0805* 青海湖裸鲤G. przewalskii 2.9728* 1.9776* 2.4 Micro CT扫描结果
本研究分别对两水系花斑裸鲤的骨骼进行Micro CT扫描, 并补充13条青海湖裸鲤的CT扫描结果作为对照, 结果显示在背鳍支鳍骨插入椎骨间的位置上有所不同。背鳍支鳍骨插入椎骨的相对位置如表 6所示: 在黄河水系花斑裸鲤中, 背鳍第一支鳍骨位于11—12椎骨间, 第二支鳍骨位于13—14椎骨间(图 7a)的样本占到总研究样本的61.9%; 在格尔木河水系花斑裸鲤的研究中, 背鳍第一支鳍骨位于9—10椎骨间, 第二支鳍骨位于11—12椎骨间(图 7b)的样本占到总研究样本的59.1%; 青海湖裸鲤在该特征上存在显著区别于两水系花斑裸鲤的形态特征, 即第一支鳍骨与第二支鳍骨插入椎骨无间隔的性状(图 7c), 且占到研究样本的61.5%。此外三者均存在背鳍第一支鳍骨位于10—11椎骨间, 第二支鳍骨位于12—13椎骨间的样品, 但用卡方检验进行统计分析, 结果显示**P<0.001, 说明三者在背鳍支鳍骨插入椎骨间的相对位置上存在极显著差异。
表 6 不同水系样品背鳍支鳍骨与椎骨相对位置统计表Table 6. The relative position of dorsal fin bone and vertebra in different drainage of sample物种Species 水系Drainage 背鳍支鳍骨与椎骨的相对位置The relative position of
the dorsal fin bone to the vertebra样本量Number of specimens 花斑裸鲤 G. eckloni 黄河水系 第一支鳍骨位于10—11; 第二支鳍骨位于12—13 7 花斑裸鲤G. eckloni 黄河水系 第一支鳍骨位于11—12; 第二支鳍骨位于13—14 13 花斑裸鲤G. eckloni 黄河水系 第一支鳍骨位于12—13; 第二支鳍骨位于14—15 1 花斑裸鲤G. eckloni 格尔木河水系 第一支鳍骨位于9—10; 第二支鳍骨位于11—12 13 花斑裸鲤G. eckloni 格尔木河水系 第一支鳍骨位于10—11; 第二支鳍骨位于12—13 9 青海湖裸鲤G. przewalskii 青海湖水系 第一支鳍骨位于10—11; 第二支鳍骨位于12—13 5 青海湖裸鲤G. przewalskii 青海湖水系 第一支鳍骨与第二支鳍骨插入椎骨无间隔 8 总计Total 56 ![]()
图 7 不同水系样品背鳍支鳍骨与椎骨相对位置的部分展示图a. 黄河水系花斑裸鲤样图; b. 柴达木水系格尔木河花斑裸鲤样图; c. 青海湖裸鲤样图Figure 7. The relative position of the dorsal fin branch bone and vertebra in different drainage of samplea. The sample of G. eckloni in Yellow River; b. The sample of G. eckloni in Golmud River; c. The sample of G. przewalskii3. 讨论
3.1 不同水系花斑裸鲤的形态分析
武云飞等[1]、曹文宣和陈毅峰[3]根据形态特征对黄河与柴达木水系格尔木河花斑裸鲤的形态学特征进行了研究, 他们基于20多项形态分析得出了一致结论, 即两水系花斑裸鲤的形态学特征没有差异, 应属同一物种的不同地理种群。此外, 赵铁桥在对134尾格尔木河、黄河上游水系及河西水系不同地理种群花斑裸鲤的形态学研究中发现, 其下颌角质内缘、第一鳃弓鳃耙数等七项鉴别特征中均存在变异性大、性状不稳定等特点, 难以作为种级鉴别特征, 其形态上的细微差异不能作为物种分类标准[4]。本文通过采用更为先进的形态测量方法, 对两水系花斑裸鲤头部形态及骨骼特征进行精确的定量分析, 得到了比之前的研究更为丰富的结果。对花斑裸鲤两个地理种群头部轮廓定量分析的结果显示与上述前人的研究结果一致, 即对头部轮廓特征点做PCA分析、判别分析和典型变量分析(图 3—图 6)均没有发现二者在头部轮廓形态上存在差异, 马氏距离与普氏距离也同样显示两水系花斑裸鲤头部轮廓的变异量很小。然而值得关注的是, 基于骨骼3D模型的比较, 发现花斑裸鲤两个地理种群在“背鳍第一支鳍骨插入椎骨位置”这一性状上有着显著性差异(图 7), 该性状在虾虎鱼属中是一个重要的鉴别特征[28]。本研究是首次发现该性状在裂腹鱼类种群间存在显著性差异, 但是否可以作为裂腹鱼类的一个鉴别特征还需今后广泛取样, 在整个裂腹鱼亚科上进行实验验证。
3.2 不同水系花斑裸鲤形态演化的生态学机制
物种形态的演化深受生态压力的影响, 特别是地质地貌的改变对物种的形态分化往往起着重要作用, 鱼类口型下颌等环境适应性特征更是与其所处的流水环境和食物组成息息相关。高原河流生态位少, 食物供给相对有限, 在高原鱼类中发生形态趋同的现象极为普遍, 如汤永涛等[20]在对高度特化等级裂腹鱼类基于简化基因组水平的研究中就发现裸鲤属与裸裂尻鱼属存在的大量形态趋同导致了分类学混乱。然而骨骼是构成脊椎动物的基本结构, 在长期的形态演化中受环境影响较小, 不仅起到支撑身体形态的作用, 也是物种分类的重要依据[29]。早在20世纪90年代, 孔晓瑜等[30]根据鱼类的骨骼特征, 头骨、椎骨及附肢骨骼的差异将中国鳜亚科分为了三个属。在虾虎鱼的分类鉴定中, Birdsong等[31]将“背鳍支鳍骨插入椎骨间的相对位置”作为一个重要的分类依据, 对200个属500种虾虎鱼进行了分类鉴定。
就本文而言, 花斑裸鲤的两个地理种群(黄河和格尔木河)基于前人在形态学上的研究应属同一物种, 但对两者分子水平上的研究显示他们是由不同祖先进化而来的并系。本研究所用的几何形态测量法, 与传统测量方法相比, 更能精确定量物种的形态特征, 但采用几何形态测量法, 分析两水系花斑裸鲤的环境适应性特征, 头部形态及下颌骨轮廓, 仍未发现两者存在显著差异, 几何形态测量法支持了传统测量结果。但通过Micro CT扫描技术对花斑裸鲤两地理种群的骨骼结构进行分析, 却观察到在“背鳍第一支鳍骨插入椎骨位置”这一更具分类学意义的特征上具有显著性差异。本文认为这很可能是由于两水系相似的生态环境造成的适应性形态趋同。青藏高原东部的黄河上游与柴达木水系的格尔木河地理相邻, 同时发源于昆仑山脉[32]。作为花斑裸鲤模式生物产地的格尔木河, 河水补给主要来源于高山融雪, 食物组成相对有限, 硅藻类是其水体生态系统的重要组成部分[33], 底栖动物仅有寡毛类和钩虾。黄河上游同处高海拔, 河水多来源于高山冰雪融水, 且食物组成中同为硅藻类占据优势[1, 34], 兼有少量摇蚊幼虫等。花斑裸鲤的两个地理种群均为杂食性鱼类, 以硅藻、摇蚊幼虫为主要食物, 两个地理种群食性结构相似, 在长期自然选择中很可能导致与食性相关的形态特征, 如口裂、下颌形状等发生了趋同进化。此外, 另一值得关注的结果是花斑裸鲤两个地理种群与青海湖裸鲤的“背鳍支鳍骨相对于椎骨位置”这一性状, 在三者的种群内部均出现了分化, 且均存在“背鳍第一支鳍骨插入第10—11椎骨间, 第二支鳍骨插入12—13椎骨间”性状, 这一结果又对两水系花斑裸鲤形态上的异同做出了另一种解释, 可能暗示了三者在“背鳍支鳍骨插入椎骨相对位置”这一性状上尚未完全分化。青海湖裸鲤与花斑裸鲤属近缘种, 且遗传分化时间短[5], 很可能在青海湖裸鲤快速成种后, 其骨骼特征分化速度较慢, 还处于分化阶段。而花斑裸鲤一直以来被认为是广布种, 具有广泛的形态多态性和不同水系的遗传分化[4, 6], 因而在青海湖裸鲤快速成种后, 两水系花斑裸鲤在头部保留了原来的特征, 支鳍骨相对于椎骨的位置却呈现一定的多态性。
总之, 鱼类的表型性状深受环境压力的影响, 在长期自然选择条件下, 一些环境适应性特征往往会发生改变, 而骨骼系统, 如椎骨数、鳍条数等则受环境压力相对很小。因而采用几何形态测量法与Micro CT扫描技术对鱼类进行分类鉴定具有十分重要的意义, 它比传统测量方法更能精确分析形态特征。近年来, 现代几何形态测量法已在裂腹鱼上得到应用, 但Micro CT扫描技术却是在花斑裸鲤为代表的裂腹鱼上的首次尝试, 并发现了“背鳍支鳍骨相对于椎骨位置”这一性状在裂腹鱼亚科不同种群间存在显著差异的特性。综上所述, 本研究丰富了花斑裸鲤不同地理种群间的形态学数据, 并在形态学角度为两水系花斑裸鲤的分类地位及种群生态学提供了理论支持。但若要清楚阐释两水系花斑裸鲤的进化关系, 真正厘清两者形态异同的生态学机制还需结合更广泛的分子数据, 挖掘更多有意义的核基因或简化基因组数据进行深入研究, 同时“背鳍支鳍骨相对于椎骨位置”这一性状是首次在裂腹鱼类上发现, Micro CT也是首次应用在裂腹鱼类的骨骼结构研究中, 今后可作为一种模式在裂腹鱼类上进行广泛的实验验证, 建立裂腹鱼类特征数据库, 判断该性状是否可作为裂腹鱼类的鉴别特征。
-
![]()
图 1 采样地点示意图
图中圆圈代表采样地点
Figure 1. Map of the head-water area of the Yellow River, Golmud River, and Qinghai Lake showing the sampling sites
The circle represents the sampling site
![]()
图 2 花斑裸鲤标志点、半标志点的选取
a. 头部侧面轮廓特征点; b. 下颌轮廓特征点
Figure 2. Landmarks and semi-landmarks of G. eckloni
a. Landmarks in the leftside of the head; b. Landmarks of the Mandibular contour
![]()
图 3 黄河与格尔木河水系花斑裸鲤头部侧面轮廓(a) 与下颌轮廓(b) 的主成分分析
Figure 3. Principal component analysis of the head profile (a) and mandibular profile (b) of the G. eckloni in Yellow River and Golmud River
![]()
图 4 花斑裸鲤头部侧面轮廓的判别分析
Figure 4. The scatter diagram based on the LDA for the cephalic side contour of the G. eckloni
![]()
图 5 花斑裸鲤下颌骨轮廓的判别分析
Figure 5. The mandible morphological variations of the G. eckloni based on LDA
![]()
图 6 黄河与格尔木河水系花斑裸鲤头部侧面轮廓(a) 与下颌轮廓(b) 的典型变量分析
Figure 6. Canonical variable analysis of the head profile (a) and mandibular profile (b) of the G. eckloni in Yellow River and Golmud River
![]()
图 7 不同水系样品背鳍支鳍骨与椎骨相对位置的部分展示图
a. 黄河水系花斑裸鲤样图; b. 柴达木水系格尔木河花斑裸鲤样图; c. 青海湖裸鲤样图
Figure 7. The relative position of the dorsal fin branch bone and vertebra in different drainage of sample
a. The sample of G. eckloni in Yellow River; b. The sample of G. eckloni in Golmud River; c. The sample of G. przewalskii
表 1 样品的采集水系、地点及数量
Table 1 Populations, numbers, river systems and sampling locations used in this study
编号No. 物种名Species 水系Drainage 采集地Sampling site 数量Number of specimens 1 花斑裸鲤 格尔木河 清水河 19 2 花斑裸鲤 格尔木河 泽格斯特河 23 3 花斑裸鲤 格尔木河 郭勒木德镇傲色河 10 4 花斑裸鲤 格尔木河 温泉水库 10 5 花斑裸鲤 黄河 扎陵湖 6 6 花斑裸鲤 黄河 鄂陵湖 13 7 花斑裸鲤 黄河 共和曲沟 5 8 花斑裸鲤 黄河 贵德 30 9 青海湖裸鲤 青海湖水系 哈尔盖河入湖口 19 总计Total 135 
下载: 导出CSV
表 2 花斑裸鲤不同种群头部侧面轮廓的马氏距离比较
Table 2 Comparison of Mahalanobis distance of head profile of the G. eckloni in Yellow River and Golmud River
下载: 导出CSV
表 5 花斑裸鲤两个地理种群下颌轮廓的普氏距离比较
Table 5 Procrustes distance comparison of mandibular contour between two geographic populations of G. eckloni
花斑裸鲤(格尔木河水系) G. eckloni in Golmud River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 0.0748* 青海湖裸鲤G. przewalskii 0.0899* 0.0216
下载: 导出CSV
表 3 花斑裸鲤不同种群头部侧面轮廓的普氏距离比较
Table 3 Comparison of Procrustes distance of head profile of the G. eckloni in Yellow River and Golmud River
花斑裸鲤(格尔木河水系) G. eckloni in Golmud River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 0.0329* 青海湖裸鲤G. przewalskii 0.0739* 0.0852*
下载: 导出CSV
表 4 花斑裸鲤两个地理种群下颌轮廓的马氏距离比较
Table 4 Mahalanobis distance of mandibular contour between two geographic populations of G. eckloni
花斑裸鲤(格尔木河水系) G. eckloni in Golmud River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 花斑裸鲤(黄河水系) G. eckloni in Yellow River 2.0805* 青海湖裸鲤G. przewalskii 2.9728* 1.9776*
下载: 导出CSV
表 6 不同水系样品背鳍支鳍骨与椎骨相对位置统计表
Table 6 The relative position of dorsal fin bone and vertebra in different drainage of sample
物种Species 水系Drainage 背鳍支鳍骨与椎骨的相对位置The relative position of
the dorsal fin bone to the vertebra样本量Number of specimens 花斑裸鲤 G. eckloni 黄河水系 第一支鳍骨位于10—11; 第二支鳍骨位于12—13 7 花斑裸鲤G. eckloni 黄河水系 第一支鳍骨位于11—12; 第二支鳍骨位于13—14 13 花斑裸鲤G. eckloni 黄河水系 第一支鳍骨位于12—13; 第二支鳍骨位于14—15 1 花斑裸鲤G. eckloni 格尔木河水系 第一支鳍骨位于9—10; 第二支鳍骨位于11—12 13 花斑裸鲤G. eckloni 格尔木河水系 第一支鳍骨位于10—11; 第二支鳍骨位于12—13 9 青海湖裸鲤G. przewalskii 青海湖水系 第一支鳍骨位于10—11; 第二支鳍骨位于12—13 5 青海湖裸鲤G. przewalskii 青海湖水系 第一支鳍骨与第二支鳍骨插入椎骨无间隔 8 总计Total 56
下载: 导出CSV
-
[1] 武云飞, 吴翠珍. 青藏高原鱼类 [M]. 成都: 四川科学技术出版社, 1992: 437-441 Wu Y F, Wu C Z. The Fishes of the Qinghai-Xizang Plateau [M]. Chengdu: Sichuan Publishing House of Science & Technology, 1992: 437-441
[2] Herzenstein S M. Wissenschaftliche Resultat der von N. M. Przewalski nach Central-Asien [J]. Zool. Theil. III, 1891, 2(3): 181-262.
[3] 陈毅峰, 曹文宣. 裂腹鱼亚科 [A]//乐佩琦. 中国动物志硬骨鱼纲鲤形目(下卷) [C]. 北京: 科学出版社, 2000: 361-362 Chen Y F, Cao W X. Schinzothoracinae [A]//Yue P Q (Eds.), Fauna Sinica Osteichthyes Cypriniformes Ⅲ [C]. Beijing: Science Press, 2000: 361-362
[4] 赵铁桥. 关于花斑裸鲤(Cymnocypris eckloni Herzenstein)的分类问题 [J]. 动物世界, 1986, 3(4): 49-55. Zhao T Q. On the taxonomic problem of Cymnocypris eckloni Herzenstein (Cyprinidae) [J]. La Animala Mondo, 1986, 3(4): 49-55.
[5] 赵凯, 李俊兵, 杨公社, 等. 青海湖及其相邻水系特有裸鲤属鱼类的分子系统发育 [J]. 科学通报, 2005, 50(13): 1348-1355. doi: 10.3321/j.issn:0023-074X.2005.13.011 Zhao K, Li J B, Yang G S, et al. Molecular phylogeny of endemic Gymnocypris in Qinghai Lake and its adjacent water [J]. Chinese Science Bulletin, 2005, 50(13): 1348-1355. doi: 10.3321/j.issn:0023-074X.2005.13.011
[6] Zhao K, Duan Z Y, Peng Z G, et al. The youngest split in sympatric schizothoracine fish (Cyprinidae) is shaped by ecological adaptations in a Tibetan Plateau glacier lake [J]. Molecular Ecology, 2009, 18(17): 3616-3628.
[7] 白明, 杨星科. 几何形态测量法在生物形态学研究中的应用 [J]. 昆虫知识, 2007, 44(1): 143-147. Bai M, Yang X K. Application of geometric morphometrics in biological researches [J]. Entomological Knowledge, 2007, 44(1): 143-147.
[8] 张颖奇, 严亚玲, 刘毅弘, 等. 几何形态测量学方法在小哺乳动物化石分类鉴定中的应用——4种䶄类化石大样本的个案研究 [J]. 古脊椎动物学报, 2012, 50(4): 361-372. Zhang Y Q, Yan Y L, Liu Y H, et al. Applying geometric morphometrics to the classification and identification of small mammals—A case study of large samples of four Arvicoline species [J]. Vertebrata Palasiatica, 2012, 50(4): 361-372.
[9] Bookstein F L. Thin-Plate Splines and the Atlas Problem for Biomedical Images [A]. Biennial International Conference on Information Processing in Medical Imaging [C]. Springer Berlin Heidelberg, 1991(511): 326-342
[10] Rohlf F J, Marcus L F. A Revolution in Morphometrics [J]. Trends in Ecology & Evolution, 1993, 8(4): 129-132.
[11] Mitteroecker P, Bookstein F L. Linear discrimination, ordination, and the visualization of selection gradients in modern morphometrics [J]. Evolutionary Biology, 2001, 38(1): 100-114.
[12] Perez S I, Bernal V, Gonzalez P N. Differences between sliding semi-landmark methods in geometric morphometrics, with an application to human craniofacial and dental variation [J]. Journal of Anatomy, 2006, 208(6): 769-784. doi: 10.1111/j.1469-7580.2006.00576.x
[13] Harda J J, Berejikianb B A, Tezakb E P, et al. Evidence for morphometric differentiation of wild and captively reared adult coho salmon: a geometric analysis [J]. Environmental Biology of Fishes, 2000, 58: 61-73.
[14] Bookstein F L. Biometrics, biomathematics and the morphometric synthesis [J]. Bulletin of Mathematical Biology, 1996, 58(2): 313-365. doi: 10.1007/BF02458311
[15] Ozpicak M, Polat N. Morphological comparison of six coastal stream populations of crimean barbel (Barbus tauricus Kessler, 1877) from the Southern Black Sea basin [J]. Turkish Journal of Fisheries and Aquatic Sciences, 2019, 19(6): 485-494.
[16] Bookstein F L, Gunz P, Mitterœcker P, et al. Cranial integration in Homo: singular warps analysis of the midsagittal plane in ontogeny and evolution [J]. Journal of Human Evolution, 2003, 44(2): 167-187.
[17] Martino A L, Dehling J M, Sinsch U. Integrative taxonomic reassessment of Odontophrynus populations in Argentina and phylogenetic relationships within Odontophrynidae (Anura) [J]. Biodiversity and Conservation, 2019, 7(82): 1-32.
[18] Tang Y T, Feng C G, Wang H K Y, et al. Taxonomic status of a population of Gymoncypris waddelli Regan, 1905 (Cypriniformes: Schizothoracinae) distributed in Pengqu River, Tibet, China [J]. Zootaxa, 2016, 4126(1): 123-137. doi: 10.11646/zootaxa.4126.1.7
[19] 王贺崐元, 汤永涛, 李国刚, 等. 高原裸裂尻鱼头部轮廓形状及其变异的几何形态测量分析 [J]. 水生生物学报, 2017, 41(1): 182-193. doi: 10.7541/2017.23 Wang H K Y, Tang Y T, Li G G, et al. Geometric morphometrics of the cephalic contour and its morphological variations among Shizopygopsis stoliczkai (Teleostei: Cyprinidae) [J]. Acta Hydrobiologica Sinica, 2017, 41(1): 182-193. doi: 10.7541/2017.23
[20] Tang Y T, Li C H, Wang H K Y, et al. Convergent evolution misled taxonomy in schizothoracine fishes (Cypriniformes: Cyprinidae) [J]. Molecular Phylogenetics and Evolution, 2019, 134: 323-337. doi: 10.1016/j.ympev.2019.01.008
[21] Wang G, Jiang M, Lee S W, et al. Cone-beam reconstruction for Micro-CT [J]. Biomedical Imaging, 2002, 11(3): 384-387.
[22] Oliviero S, Lu Y, Viceconti M, et al. Effect of integration time on the morphometric, densitometric and mechanical properties of the mouse tibia [J]. Journal of Biomechanics, 2017, 65: 203-211.
[23] Curtis A A, Simmons N B. Unique turbinal morphology in horseshoe bats (Chiroptera: Rhinolophidae) [J]. The Anatomical Record, 2017, 300(2): 309-325.
[24] Sakashita M, Sato M, Kondo S. Comparative morphological examination of vertebral bodies of teleost fish using high-resolution micro-CT scans [J]. Journal of morphology, 2019, 280(8): 1-18.
[25] Rohlf F J. The tps series of software [J]. Hystrix, the Italian Journal of Mammalogy, 2015, 26(1): 9-12.
[26] Klingenberg C P. MorphoJ: an integrated software package for geometric morphometrics [J]. Molecular Ecology Resources, 2011, 11(2): 353-357.
[27] Hammer O, Harper D A T, Ryan P D, et al. PAST-Palaeontological statistics, ver. 1.89 [J]. Palaeontologia Electronica, 2009, 4(1): 1-92.
[28] 伍汉霖, 钟俊生. 中国动物志硬骨鱼纲鲈形目(五) [M]. 上海: 上海水产大学, 2007: 43-45 Wu H L, Zhong J S. Zoology of China Osteichthyes Perciformes (V) [M]. Shanghai: Shanghai Ocean University, 2007: 43-45
[29] 王映, 赵金良, 柯巧珍, 等. 大黄鱼养殖群体和野生群体骨骼micro CT扫描显微结构比较 [J]. 研究报告, 2016, 40(12): 36-40. Wang Y, Zhao J L, Ke Q Z, et al. Comparison of bone microstructures in cultured and wild stocks of large yellow croaker (Larimichthys crocea) using micro CT scanning [J]. Reports, 2016, 40(12): 36-40.
[30] 孔晓瑜, 周才武. 中国鳜亚科(Sinipercinae)七种鱼类骨骼形态特征的比较研究 [J]. 青岛海洋大学学报, 1993, 23(3): 116-123. Kong X Y, Zhou C W. Characteristics of seven Sinipercinae fishes of china [J]. Journal of Ocean University of Qingdao, 1993, 23(3): 116-123.
[31] Birdsong R S, Murdy E O, Pezold F L. A Study of the vertebral column and median fin osteology in gobioid fishes with comments on gobioid relationships [J]. Bulletin of Marine Science, 1988, 42(2): 174-214.
[32] 赵凯, 杨公社, 李俊兵, 等. 黄河裸裂尻鱼群体遗传结构和Cyt b序列变异 [J]. 水生生物学报, 2006, 30(2): 129-133. doi: 10.3321/j.issn:1000-3207.2006.02.001 Zhao K, Yang G S, Li J B, et al. Phylogenetic structure of Schizopygopsis pylzovi populations from mitochondrial cytochrome b gene sequence variations [J]. Acta Hydrobiologica Sinica, 2006, 30(2): 129-133. doi: 10.3321/j.issn:1000-3207.2006.02.001
[33] 陈燕琴, 申志新, 刘玉婷. 柴达木水系格尔木河秋季浮游植物群落结构特征初步研究 [J]. 青海农林科技, 2012(2): 10-16. doi: 10.3969/j.issn.1004-9967.2012.02.005 Chen Y Q, Shen Z Z, Liu Y T. Study on community structure of phytoplankton in Geermu River of Chaidamu River system in autumn [J]. Science and Technology of Qinghai Agriculture and Forestry, 2012(2): 10-16. doi: 10.3969/j.issn.1004-9967.2012.02.005
[34] 张军燕, 张建军, 杨兴中, 等. 黄河上游玛曲段春季浮游生物群落结构特征 [J]. 生态学杂志, 2009, 28(5): 983-987. Zhang J Y, Zhang J J, Yang X Z, et al. Characteristics of plankton community structure in Maqu section at upper reaches of Yellow River in spring [J]. Chinese Journal of Ecology, 2009, 28(5): 983-987.
-
期刊类型引用(9)
1. 刘莹，于超勇，赵文溪，王英俊，侯富晟，邹琰. 基于几何形态测量学方法的大菱鲆不同群体形态变异分析. 中国海洋大学学报(自然科学版). 2025(01): 66-75 . 
百度学术
2. 寇若彬，刘丹，高强，晁燕，张存芳，聂苗苗，谭瑾，郭守全，祁得林. 花斑裸鲤红细胞发育相关基因KLF17启动子区域互作蛋白筛选. 水生生物学报. 2024(05): 808-819 . 本站查看
3. 黄青峰，邓楚可，夏静怡，严红娟，夏继刚. 秦岭细鳞鲑与同域物种拉氏鱥几何形态学特征：生活史阶段效应及种间差异. 生态学杂志. 2024(04): 922-929 . 百度学术
4. 李永项，李兆雨，杨怀龙，谢坤. 青海共和盆地发现鱼类化石. 第四纪研究. 2024(06): 1713-1716 . 百度学术
5. 牛依萌，刘思嘉，田菲，冯晨光，赵凯. 渭河高原鳅线粒体基因组特征及其系统发育分析. 基因组学与应用生物学. 2024(Z2): 1813-1825 . 百度学术
6. 刘鸿，牛建功，胡江伟，海萨·艾也力汗，咸玉兰，贺疆滔，蔡林钢，张人铭. 新疆叶尔羌河流域4种裂腹鱼多变量形态学对比分析. 水产学报. 2023(04): 77-90 . 百度学术
7. 周启苓，王刘永，杨云生，马骞，吴雨薇，陈刚. 不同发育期虹鳟脊椎骨的显微结构及钙、磷元素含量分析. 渔业科学进展. 2023(05): 115-124 . 百度学术
8. 唐静，薛竣仁，刘洪波，姜涛，陈修报，杨健. 中华绒螯蟹螯足长节产地差异的几何形态测量学研究. 安徽农业大学学报. 2023(05): 802-808 . 百度学术
9. 朱国平，刘芳沁. 几何形态测量学及其在鱼类生态学研究中的应用进展. 上海海洋大学学报. 2022(05): 1180-1189 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 3227
- HTML全文浏览量: 947
- PDF下载量: 106
- 被引次数: 14