Structure of mitochondrial DNA control region and genetic diversity of Misgurnus anguillicaudatus
-
摘要:
采用PCR和DNA测序技术对4个泥鳅(Misgurnus anguillicaudatus)群体的线粒体DNA控制区序列进行比较，研究其遗传多样性和控制区的结构。结果显示，用于分析的线粒体DNA控制区片段序列为918~920 bp。32个序列中共发现了56个多态位点，其中32个转换位点，19个颠换位点，5个转换与颠换同时存在的位点，定义了32种单倍型。同时对控制区结构进行分析，识别了其终止序列区(ETAS)、中央保守区(CD)和保守序列区(CSB)的关键序列。4个地理群体的单倍型多样性(Hd)、核苷酸多样性(Pi)和平均核苷酸差异数(K)分别为0.992、0.012和10.698。群体间的平均Kimura双参数遗传距离(Kimura 2-parameter distance，K2-P)、遗传分化指数(Fst)、基因交流值(Nm)和分子方差分析(AMOVA)均表明4个泥鳅群体具有较高的遗传分化，基因交流贫乏。利用核苷酸序列构建的NJ分子系统树揭示4个群体分为2个谱系，除范县群体外，其余各群体间相互交叉聚集。
Abstract:We compare and study the genetic diversity and structure of DNA control region of 4 Misgurnus anguillicaudatus populations by using PCR and DNA sequencing. The results show that the sequence of mitochondrial DNA control region is 918~920 bp.Among 32 sequences, a total of 56 polymorphic loci were detected, including 32 transitions, 19 transversions, 5 transitions-transversions and 36 haplotypes. The key sequences of extended terminal associated sequences (ETAS), central domains (CSB-A, B, C, D, E and F) and conserved sequence blocks (CSB-1, CSB-2 and CSB-3) were discovered. The haplotypic diversity (Hd), average genetic distance (Pi) and average number of pairwise differences (K) is 0.992, 0.012 and 10.698, respectively. The average Kimura 2-parameter distance, Fst value, Nm value and AMOVA reveal great genetic differentiation among the 4 M.anguillicaudatus populations. The NJ phylogenetic tree constructed on nucleotide sequence indicates that the 4 populations form into 2 monophyletic lineages, and cross and cluster among one another except Fanxian population.
-
鸢乌贼(Sthenoteuthis oualaniensis)隶属枪形目、开眼亚目、柔鱼科、鸢乌贼属，广泛分布于印度洋、太平洋的热带和亚热带海域，其中在南海和印度洋西北部海域数量较大[1-2]。鸢乌贼具有非常复杂的种群结构。NESIS[3]将其分为大型种群、单内壳叶柄侧轴中型种群、双内壳叶柄侧轴中型种群(典型的鸢乌贼种群)、小型种群、微型种群等5个种群，后3个种群主要分布于太平洋海域。微型种群背部没有发光器官，主要分布在10°N~10°S的赤道海域，分布区域与双内壳叶柄侧轴中型种群重叠[4]。国内外学者已对印度洋西北海域鸢乌贼的种群结构、年龄和生长等做了较为深入的研究[2-7]。
南海的鸢乌贼主要分布于南海中南部海盆区，已成为当前南海外海渔业开发的热点[8]。东南亚渔业开发中心(SEAFDEC)、中国台湾和大陆都曾开展过南海鸢乌贼资源调查，并对其资源量分布、年龄与生长等生物和生态学习性进行过研究[9-18]。SEAFDEC在菲律宾西部海域的调查发现，南海存在中型和微型2个鸢乌贼种群，微型种群早熟且无背发光器，胴长范围9.8~12 cm，中型种群具背发光器，胴长范围12~24 cm[15]。颜云榕等[13]研究了西中沙海域鸢乌贼夏季渔业生物学，发现雌性个体出现2个性腺发育高峰期，即101~120 mm(小型个体，胴部较硬，胴背无卵圆形光斑)与161 mm以上(中型个体，胴部柔软，胴背有卵圆形光斑)胴长组。种群问题是渔业生物学研究的基础，文章根据2011年春季灯光罩网渔船南沙探捕所获资料，研究了南沙海域鸢乌贼的种群组成、胴长分布、雌雄比例、性腺发育等种群结构特征，旨在为今后深入开展南海鸢乌贼生物学研究、合理开发利用资源提供参考依据。
1. 材料与方法
1.1 资料来源
资料来源于2011年春季南海中南部海域的灯光罩网调查。调查船为“桂防渔96886”，钢质，总长41.80 m，总吨位413 GT；渔船主机1台，功率382 kW，发电机2台，功率分别为200 kW和150 kW，配340盏(×1 kW)金属卤化物集鱼灯及HE-670型垂直探鱼仪，包括船长在内船员7人，渔获加冰冷藏保存。调查所用罩网网具沉子纲长280.0 m，拉直高度82.5 m，网衣最大网目3.5 mm，网囊最小网目2.0 mm。
1.2 调查方法与数据处理
整个调查为生产性探捕，不设具体调查站位。航次时间为2011年4月19日~5月6日，历时26 d，作业18晚，作业海域范围9°39′N~14°54′N、113°16′E~115°30′E，其中前16晚在南沙海域(12°N以南)作业，后2晚是返航途中在中沙海域生产。
渔船每晚作业10~16网次不等，记录每个网次的渔获组成，每晚选择部分网次随机抽取50尾鸢乌贼进行种群、性别和胴长测定，部分网次解剖观测性腺成熟度和摄食强度。根据胴背发光器官的有无和性腺发育的快慢[3, 15]，将鸢乌贼分成中型和微型2个种群；根据左侧第4腕是否茎化和性腺结构的差异[1, 5]，进行鸢乌贼雌雄的判别；胴长用精度为1 mm的量鱼板测定；参照刘必林[5]将鸢乌贼性腺成熟度划分为性未成熟(Ⅰ期和Ⅱ期)、性成熟(Ⅲ期和Ⅳ期)、产卵后(Ⅴ期)5期3个阶段。因为第1晚的鸢乌贼渔获未按种群分开，文章仅对12°N以南海域15个夜晚的数据进行分析，具体采样站位见图 1。
以10 mm组距对胴长进行分组，例如100 mm胴长组表示胴长范围96~105 mm；运用t检验比较2个种群和雌雄个体之间的胴长分布是否存在显著性差异；运用D检验(Kolmogorov-Smirnov检验)分析胴长组成是否符合正态分布；利用适合性卡方(χ2)测验，验证鸢乌贼的雌雄比例是否与1 : 1的理论比例相符；利用不同胴长组性成熟个体百分比拟合Logistic曲线，分种群和性别推算鸢乌贼的初次性成熟胴长[19]：
$$ {P_i} = \frac{1}{{1 + {{\rm{e}}^{ - \left( {a + b{L_l}} \right)}}}};\quad M{L_{50\% }} = - \frac{a}{b} $$ 式中Pi为成熟个体占组内样本的百分比；Li为胴长组，ML50%为初次性成熟胴长(mm)。
2. 结果
2.1 种群组成
调查船15个夜晚作业194网次(表 1)，平均每晚作业12.9网，合计捕捞鸢乌贼68.06 t，网次产量为28.5~1 282.5 kg，平均为350.8 kg [标准差(SD)=237.3 kg]，每晚产量为1.54~10.94 t，平均为4.54 t(SD=2.69 t)。鸢乌贼每个网次都有捕获，但网次和每晚产量的变幅较大。随机观测了119网次5 964尾鸢乌贼的种群、性别和胴长。中型种群占抽样总尾数的88.92%，胴背有发光器官，性成熟个体较大；微型种群占抽样总尾数的11.08%，胴背无发光器官，性成熟个体较小。中型种群是南沙海域鸢乌贼的优势种群。微型种群15个夜晚都有捕获，占每晚抽样尾数的2.75%~30.29%；在119个抽样网次中出现110次，占网次抽样尾数的0~72.00%。微型种群每晚和网次所占比例的变幅较大。
表 1 取样站位和鸢乌贼渔获概况Table 1. Sampling stations and catches of S.oualaniensis夜晚
No.日期
date作业地点 position 作业网次
fishing times产量/t
yield抽样网次
sampling times抽样尾数/尾
sampling number微型占比/%
drawf-form/totalLat.(N) Long.(E) 1 04-20 114°44′ 11°24′ 10 1.65 4 200 4.00 2 04-21 114°02′ 10°59′ 13 3.33 13 651 4.30 3 04-22 114°31′ 11°10′ 10 1.54 10 501 12.57 4 04-23 113°52′ 09°53′ 12 3.90 10 501 3.59 5 04-24 113°44′ 09°48′ 12 2.91 9 453 5.30 6 04-25 113°16′ 09°39′ 15 3.59 10 505 15.84 7 04-26 113°51′ 10°29′ 12 2.68 6 302 16.56 8 04-27 115°01′ 11°22′ 16 5.10 8 401 19.95 9 04-28 115°07′ 11°51′ 14 10.94 7 349 10.32 10 04-29 115°10′ 11°50′ 13 4.62 7 351 16.52 11 04-30 115°07′ 11°49′ 11 3.08 7 350 30.29 12 05-01 114°56′ 11°25′ 13 6.90 7 350 3.14 13 05-02 114°54′ 11°25′ 16 5.87 7 350 18.29 14 05-03 114°56′ 11°24′ 16 9.15 8 400 2.75 15 05-04 114°55′ 11°25′ 11 2.79 6 300 8.00 2.2 胴长分布
鸢乌贼中型种群渔获胴长范围为69~231 mm，平均胴长为121.4 mm(SD=16.4 mm)；微型种群渔获胴长范围为51~125 mm，平均胴长为88.6 mm(SD=13.4 mm)。2个种群的胴长分布有重叠，中型种群的胴长分布范围和平均胴长明显大于微型种群(图 2和表 2)。t检验表明，2个种群的胴长分布存在显著性差异(P＜0.05)。D检验表明，中型和微型种群的胴长组成均不呈正态分布(P＜0.05)。
![]()
图 2 中型和微型种群鸢乌贼胴长分布Figure 2. Distribution of mantle length of medium-form and dwarf-form S.oualaniensis表 2 鸢乌贼渔获胴长组成Table 2. Composition of mantle length of S.oualaniensis种群
population性别
sex尾数
number胴长范围/mm
mantle length range优势胴长 dominant mantle length 平均胴长/mm
average mantle length标准差/mm
SD范围/mm 比例/% 中型 medium-form 雌性 2 571 71~231 110~140 82.77 125.5 20.1 雄性 2 732 69~153 100~130 93.23 117.5 10.3 合计 5 303 69~231 100~140 90.01 121.4 16.4 微型 dwarf-form 雌性 409 56~125 80~110 86.31 93.5 14.1 雄性 252 51~104 70~90 95.24 80.8 6.9 合计 661 51~125 70~110 94.40 88.6 13.4 中型种群雌性平均胴长125.5 mm(SD=20.1 mm)，雄性平均胴长117.5 mm(SD=10.3 mm)；微型种群雌性平均胴长93.5 mm(SD=14.1 mm)，雄性平均胴长80.8 mm(SD=6.9 mm)。2个种群雌性的胴长分布范围和平均胴长均大于雄性。t检验表明，2个种群雌、雄个体之间的胴长分布均存在显著性差异(P＜0.05)。
2.3 雌雄比例
随机抽样的5 303尾中型种群鸢乌贼中，雌性个体占48.48%。卡方检验结果显示，总性比符合1 : 1的关系，雌雄组成无显著差异(P＞0.05)；70 mm、80 mm、90 mm和100 mm胴长组中，雌性个体分别占55.56%、46.88%、48.48%和45.53%，雌雄组成均无显著差异(P＞0.05)；110 mm、120 mm、130 mm、140 mm和150 mm体长组中，雌性个体分别占40.51%、43.64%、40.40%、73.58%和90.96%，雌雄组成均有显著差异(P＜0.05)。70~100 mm胴长组，雌雄比例接近；110~130 mm胴长组，雄性个体明显多于雌性；140 mm胴长组开始，雌性个体明显多于雄性；160 mm胴长组以上全为雌性(图 3-a)。D检验表明，中型种群雌雄胴长组成均不呈正态分布(P＜0.05)，雌性胴长分布的偏度系数(skewness)为1.16，分布为右偏，雄性胴长分布的偏度系数为-0.22，分布为左偏。
![]()
图 3 鸢乌贼中型(a)和微型(b)种群渔获性别组成Figure 3. Sex composition of medium-form (a) and dwarf-form (b) S.oualaniensis随机抽样的661尾微型种群鸢乌贼中，雌性个体占61.88%。卡方检验结果显示，总性比不符合1 : 1的关系，雌性个体明显多于雄性(P＜0.05)；60 mm和70 mm胴长组中，雌性个体分别占54.55%和37.88%，雌雄组成均无显著差异(P＞0.05)；80 mm、90 mm和100 mm体长组中，雌性个体分别占38.67%、71.43%和92.94%，雌雄组成均有显著差异(P＜0.05)。50~70 mm胴长组，雌雄比列接近；80 mm胴长组，雄性个体明显多于雌性；90 mm胴长组开始，雌性个体明显多于雄性，110 mm胴长组以上全为雌性(图 3-b)。D检验表明，微型种群雌、雄胴长组成均不呈正态分布(P＜0.05)，雌性胴长分布的偏度系数为0.10，分布为右偏，雄性胴长分布的偏度系数为-0.35，分布为左偏。
总体而言，鸢乌贼2个种群不同胴长组性比表现出相似的规律性：随着胴长的增加，雌雄比例由开始阶段的基本接近，到雄性个体占优，再到雌性个体占优，直至全部为雌性。2个种群雌性和雄性的胴长组成均不呈正态分布，雌性胴长分布都为右偏，雄性胴长分布都为左偏。
2.4 性腺发育
调查期间，合计解剖观测了3 833尾鸢乌贼(中型种群3 182尾，微型种群651尾)的性腺成熟度，其中随机抽样3 363尾(中型种群2 926尾，微型种群437尾，分属68个网次，抽样渔获全部解剖)，非随机抽样470尾(中型种群256尾，微型种群214尾，分属51个网次，抽样渔获少量解剖，小个体渔获为主，解剖主要是为了确定抽样个体的种群或雌雄)。
随机抽样的3 363尾鸢乌贼中，中型种群性成熟个体占15.07%，微型种群性成熟个体占42.11%，微型种群的性腺发育明显早于中型种群(表 3)。中型种群雌性和雄性中性成熟个体分别占5.12%和24.29%，微型种群雌性和雄性中性成熟个体分别占39.51%和45.36%，2个种群雄性个体的性腺发育都早于雌性个体。2个种群都未捕到产卵后(Ⅴ期)个体。
表 3 鸢乌贼性腺成熟度组成Table 3. Composition of sex maturity of S.oualaniensis种群
population性别
sex样本数/尾
number性腺成熟度/% sexual maturity 性未成熟
immature性成熟
mature产卵后
spawned中型 medium-form 雌性 1 407 94.88 5.12 0 雄性 1 519 75.71 24.29 0 合计 2 926 84.93 15.07 0 微型 dwarf-form 雌性 243 60.49 39.51 0 雄性 194 54.64 45.36 0 合计 437 57.89 42.11 0 解剖观测的3 833尾鸢乌贼中，中型种群雌、雄性成熟个体最小胴长分别为123 mm和111 mm，雌、雄性未成熟个体最大胴长分别为199 mm和133 mm；微型种群雌、雄性成熟个体最小胴长分别为80 mm和73 mm，雌、雄性未成熟个体最大胴长分别为123 mm和92 mm。根据不同胴长组的性成熟个体百分比，拟合推算了鸢乌贼初次性成熟胴长(图 4和图 5)。
![]()
图 4 中型种群鸢乌贼初次性成熟胴长Figure 4. Mantle length at first maturity of medium-form S.oualaniensis![]()
图 5 微型种群鸢乌贼初次性成熟胴长Figure 5. Mantle length at first maturity of dwarf-form S.oualaniensis中型种群雌、雄个体的性成熟度-胴长曲线分别为：
雌性：${P_i} = \frac{1}{{1 + {{\rm{e}}^{ - \left( { - 15.45\;2 + 0.086{L_i}} \right)}}}}\;\;\;\;\left( {{R^2} = 0.963, \;\;n = 576} \right) $；
雄性：${P_i} = \frac{1}{{1 + {{\rm{e}}^{ - \left( { - 29.89\;3 + 0.238{L_i}} \right)}}}}\;\;\;\;\left( {{R^2} = 0.998, \;\;n = 606} \right) $；
中型种群初次性成熟胴长(ML50%)，雌性为179.7 mm，雄性为125.6 mm。
微型种群雌、雄个体的性成熟度-胴长曲线分别为：
雌性：${P_i} = \frac{1}{{1 + {{\rm{e}}^{ - \left( { - 29.87\;7 + 0.315{L_i}} \right)}}}}\;\;\;\;\left( {{R^2} = 0.995, \;\;n = 401} \right) $；
雄性：${P_i} = \frac{1}{{1 + {{\rm{e}}^{ - \left( { - 18.16\;1 + 0.224{L_i}} \right)}}}}\;\;\;\;\left( {{R^2} = 0.998, \;\;n = 250} \right) $；
微型种群初次性成熟胴长(ML50%)，雌性为94.8 mm，雄性为81.1 mm。
3. 讨论
3.1 南沙海域鸢乌贼种群和雌雄判别
NESIS[3]将分布在太平洋和印度洋海域的鸢乌贼分为3个主要种群和2个次要种群。大型种群仅分布在印度洋北部的红海、亚丁湾和阿拉伯海，在阿拉伯海其胴长为400~500 mm，最大胴长为650 mm；中型种群存在于整个鸢乌贼种群分布区域内，成熟雄性胴长为120~150 mm，成熟雌性胴长为190~250 mm；微型种群分布在赤道附近，无背发光器，成熟雄性胴长为90~100 mm，成熟雌性胴长为90~120 mm，最大胴长为140~150 mm。中型种群根据内壳特征(双或者单内壳叶柄侧轴)又可以分为2个群体，其中单内壳叶柄侧轴群体仅出现在红海、阿丁湾和阿拉伯海15°N~17°N海域。小型种群形态近似中型种群，但性成熟个体更小，雌性胴长范围为90~160 mm，成熟胴长为120~140 mm，近似于微型种群，该种群主要分布在东热带太平洋和西印度洋。微型种群具有许多与中型种群相区别的形态学特征，如微型种群的胴背部没有发光器官、茎化腕的差异以及精囊结构和叶柄侧轴结构的差异等[4]。在仔稚鱼和幼鱼阶段，微型种群与中型种群在形态和外观上不存在差异，因此，以往的研究者并不赞同把微型种群划分为另外一个种类[4, 20]，但STAAF等[21]对中太平洋海域鸢乌贼分子系统学的研究显示，微型和中型种群之间的遗传差异较大，可能已达到物种水平。
该研究根据胴背发光器的有无，结合性腺发育的快慢，将南沙海域的鸢乌贼分为中型和微型2个种群。调查中发现，中型种群胴长大于120 mm的个体，胴背发光器官呈明显的黄色椭圆形斑块，胴长小于120 mm的个体，发光器官则不明显，容易误判为微型种群，但仔细观察胴背位置仍可发现密布的淡白色小斑点；解剖也可发现，微型种群的性腺发育要明显早于相同胴长的中型种群。鸢乌贼雄性个体左侧第4腕茎化，但未成熟个体的茎化腕有时并不明显，此时也需观察内部性腺结构判别雌雄；鸢乌贼雄性的精巢从小就清晰可辨，雌性的卵巢在性腺发育成熟时才充实成形。因此，根据外形特征结合解剖观测，可以较为准确的进行鸢乌贼2个种群与雌雄的判别。
该研究参照LABE[15]对菲律宾西部海域鸢乌贼种群的研究，将胴背具发光器官的鸢乌贼都归为中型种群。调查发现的中型种群雌、雄性成熟个体胴长范围分别为123~231 mm和111~153 mm，推算的ML50%胴长分别为179.7 mm和125.6 mm。相比NESIS[3]的分类标准，中型种群性成熟胴长变小，雌性还捕获10尾胴长140 mm以下的性成熟个体，这是由于中型种群对南海水文环境的适应性变化，还是由于混有小型种群等其他鸢乌贼种群，有待深入分析研究。
3.2 采样方法对鸢乌贼渔获种群和性比的影响
SEAFDEC调查船1998年4月~5月在菲律宾西部海域(11°N~20°N，117°E~120°E)的鱿钓调查显示，鸢乌贼是唯一渔获种类，产量合计2 538尾，393.1 kg，自动钓机平均CPUE(单位捕捞努力量渔获量)6尾· (线· h)-1；调查海域存在2个鸢乌贼种群，中型种群胴长范围为12~24 cm，占总尾数的86%(2 175尾，377.5 kg)；微型种群胴长范围为9.8~12 cm，占总尾数的14%(363尾，15.6 kg)；中型种群中雌性个体占83%，微型种群中雌性个体占70%，雌性明显占优；微型种群每个站位都有捕获(占比2.6%~31.2%)，渔获量少但遍布整个调查海域[15]。
与之比较，该次调查也发现南沙海域存在中型和微型2个鸢乌贼种群，分别占总渔获尾数的88.92%和11.08%，中型种群为优势种群，微型种群资源量少但也遍布整个调查海域。这些发现与菲律宾西部海域的调查基本一致。与之不同的是，该次调查中小个体渔获明显增多，98 mm以下个体占总渔获的11.05%，中型种群小于120 mm的个体占53.12%；渔获的雌雄比例相对均衡，中型和微型种群中雌性个体分别占48.48%和61.88%。采样方法的不同应该是造成这种差异的主要原因。SEAFDEC调查是鱿钓取样，而该调查是灯光罩网取样。实践证明，鱿钓作业选择性较强，渔获以大个体为主且雌性居多[6, 11, 17]。颜云榕等[13]2010年6月~7月在南沙北部海域(12°N~16°N，112°E~116°E)开展的自动鱿钓和灯光罩网联合调查也发现，鱿钓和罩网鸢乌贼渔获胴长范围分别为108~215 mm和65~231 mm，雌雄性比分别为32.0和1.4，钓具的总体雌雄比明显较高。ZAKARIA[17]认为，由于小个体鸢乌贼样品的缺失，鱿钓渔获结构无法反映真实的海域资源结构。相对而言，罩网网目尺寸小、选择性差、渔获量大，渔获组成更接近海域资源实际。
调查期间，鸢乌贼微型种群网次和每晚所占比例的变幅较大，即使同个夜晚的相邻网次中也经常出现这种现象。分析原因，一方面是受资源分布时空差异的影响，另一方面也与灯光罩网捕捞方法有关。罩网作业是夜晚通过集鱼灯将捕捞对象诱集至表层后再放网罩扣捕捞，资源分布水层越深逃逸率越高，调查船所用网具的有效捕捞水深约40 m。张引[11] 分析台湾调查船声学数据时发现，当某处存在大个体鸢乌贼时，小个体鸢乌贼倾向栖息于更深水层，当大个体鸢乌贼消失后，小个体鸢乌贼则上升到大个体鸢乌贼原先所在的水层；认为同类相残是小个体鸢乌贼栖息于更深水层的原因，因为它们要避开大个体鸢乌贼的捕食。微型种群鸢乌贼体型较小，栖息水层相对较深，放网时所处水层的不同和逃逸率的差异明显影响其在渔获中所占的比例。该次调查未捕到产卵后(Ⅴ期)个体，推测鸢乌贼产卵后期不再趋光并且栖息于深水层直至死亡，灯光罩网无法有效捕捞，具体还有待进一步验证。
3.3 鸢乌贼不同胴长组性比与性腺发育的关系
历史研究表明，鸢乌贼两性比例会随性成熟及繁殖情况而有变化，在未成熟群体中，性比大体为1 : 1；在繁殖初期，性成熟较早、行动较快的雄性个体先行到达繁殖场，此时雄性的比例甚高于雌性；在繁殖后期，完成交配活动的雄性个体又先行死去，从繁殖场大量消失，此时雌性的比例甚高于雄性[1]。
该研究发现，鸢乌贼中型和微型种群不同胴长组性比表现出相似的规律性：随着渔获胴长的增加，雌雄比例由开始阶段的基本接近，到雄性个体占优，再到雌性个体占优，直至全部为雌性，其中尤以雄性占优到雌性占优的转变过程最为明显(图 3)。陈新军[22]进行新西兰海域双柔鱼研究时有过类似发现，渔获物性别组成生长初期接近1 : 1，成长期雄性个体的比例略占优势，繁殖期雌性个体占优势。该次调查发现的鸢乌贼雄性最小成熟胴长，中型和微型种群分别为111 mm和73 mm，分属110 mm和70 mm胴长组，估算的雄性ML50%胴长，中型和微型种群分别为125.6 mm和81.1 mm，分属130 mm和80 mm胴长组。中型种群70~100 mm和微型种群50~60 mm胴长组，鸢乌贼处于未成熟阶段，此时的雌雄比例接近；中型种群110~130 mm和微型种群70~80 mm胴长组，为雄性鸢乌贼开始性成熟并大规模交配繁殖阶段，此时雄性比例要高于雌性；中型种群140 mm和微型种群90 mm胴长组开始，因为完成交配活动的雄性个体先行死去，雌性比例又开始高于雄性直至全部为雌性。因此，鸢乌贼不同胴长组雌雄比例的变化可能与雄性个体的生长发育密切相关。
该次探捕中每个站位都是成体和幼体混捕，中型种群110~130 mm和微型种群70~80 mm胴长组雄性比例偏高，无法用雄性个体先行到达繁殖场理论解释。ZAKARIA[17]研究发现：南海鸢乌贼的生命周期在1年之内，雌性寿命长于雄性，胴长达到100 mm后雌性个体的生长速度快于雄性。因此，中型种群110~130 mm胴长组雄性比例偏高，可能是由于雌性生长相对较快，导致相应胴长组雄性个体积压所致。微型种群70~80 mm胴长组应该可用同样的理论解释。2个种群雄性胴长分布都为左偏，雌性胴长分布都为右偏，应该也与雌雄个体寿命的长短、生长的快慢有关。
中型和微型种群抽样渔获中雌性个体分别占48.48%和61.88%，卡方检验显示，中型种群雌雄组成无显著差异(P＞0.05)，微型种群雌雄组成则有显著差异(P＜0.05)。南海鸢乌贼的主要产卵季节是夏季[11, 13]，该次探捕为春季调查，中型和微型种群中性成熟个体分别占15.07%和42.11%，微型种群的性腺成熟度明显高于中型种群，其中，雄性中性成熟个体所占比例，微型和中型种群分别为45.36%和24.29%。微型种群雌性比例偏高，应该是由于部分雄性个体已完成产卵活动先行死去，中型种群还未到产卵高峰期，雌雄比例相对均匀。因此，鸢乌贼总渔获中雌雄比例的变化也与雄性个体的生长发育密切相关。
-
![]()
图 1 基于线粒体控制区序列变异构建的4个泥鳅群体邻接树节点处数值为1 000次bootstrap检验的支持率
Figure 1. Neighbor-joining tree based on sequence variation of mtDNA control region of 4 M.anguillicaudatus populations Numbers at the nodes indicate bootstrap values with 1 000 replicate
表 1 泥鳅样本数量、体长与体质量
Table 1 Number, length and weight of M.anguillicaudatus
群体population 数量/尾number 体长/cm length 体质量/g weight 怀远Huaiyuan(HY) 9 12.56±2.14 15.43±5.43 范县Fanxian(FX) 9 11.78±1.34 13.98±5.31 舒城Shucheng(SC) 9 11.63±1.52 14.45±4.36 池州Chizhou(CZ) 9 12.63±3.26 16.86±6.29 
下载: 导出CSV
表 2 4个泥鳅群体内(对角线)与群体间的平均K2-P遗传距离
Table 2 Average genetic distance within (diagonal) and among 4 M.anguillicaudatus populations based on Kimura 2-parameter model
群体population HY FX SC CZ 怀远Huaiyuan(HY) 0.008 范县Fanxian(FX) 0.011 0.006 舒城Shucheng(SC) 0.011 0.014 0.004 池州Chizhou(CZ) 0.016 0.014 0.015 0.012
下载: 导出CSV
表 3 4个泥鳅群体内遗传多样性参数
Table 3 arameters of genetic diversity within 4 M.anguillicaudatus populations
群体population 样本数number of samples 单倍型数number of haplotypes 单倍型多样性haplotype diversity (Hd) 多态位点数polymorphic loci (s) 平均核苷酸差异数average number of pairwise difference (K) 核苷酸多样性nucleotide diversity (Pi) 怀远Huaiyuan(HY) 9 8 0.972 23 8.028 0.009 范县Fanxian(FX) 9 8 0.972 22 6.611 0.007 舒城Shucheng(SC) 9 7 0.917 14 3.278 0.004 池州Chizhou(CZ) 9 9 1.000 37 11.306 0.012 合计total 36 32 0.992 56 10.698 0.012
下载: 导出CSV
表 4 4个泥鳅群体间的遗传分化指数Fst(对角线下方)和基因交流值Nm(对角线上方)
Table 4 Genetic differentiation Fst value (below diagonal) and gene flow Nm value (above diagonal) among 4 M.anguillicaudatus populations
群体population HY FX SC CZ 怀远Huaiyuan(HY) 0.490 0.274 0.468 范县Fanxian(FX) 0.338 0.127 0.470 舒城Shucheng(SC) 0.477 0.663 0.238 池州Chizhou(CZ) 0.348 0.347 0.512
下载: 导出CSV
表 5 4个泥鳅群体的AMOVA分析
Table 5 AMOVA analysis of 4 M.anguillicaudatus populations
变异来源source of variation 自由度df 方差总和sum of squares 变异组分variance components 变异百分比/% percentage of variation 群体间among populations 3 97.809 3.137 38.980 群体内个体间within populations 32 157.110 4.910 61.020 总计total 35 254.917 8.046 -
下载: 导出CSV
-
[1] CANN R L, BROWN W M, WILSON A C. Polymorphic sites and the mechanism of evolution in human mitochondrial DNA[J]. Genetics, 1984, 106(2): 479-499. doi: 10.1093/genetics/106.3.479
[2] ROSEL P E, DIZON A E, HAYGOOD M G. Variability of the mitochondrial control region in populations of the harbour porpoise, Phocoena phocoena, on interoceanic and regional scales[J]. Can J Fish Aquat Sci, 1995, 52(3): 1421-1429. doi: 10.1139/F95-118
[3] ROSEL P E, HAYGOOD M G, PERRIN W F. Phylogenetic relationship among the true porpoise (Cetacean: Phocoenidae)[J]. Mol Phyl Evol, 1995, 4(4): 463-474. doi: 10.1006/mpev.1995.1043
[4] 孙玉华, 王伟, 刘思阳, 等. 中国胭脂鱼线粒体控制区遗传多样性分析[J]. 遗传学报, 2002, 29(9): 787-790. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsQ0hJTmV3UzIwMjQxMTA1MTcxMzA0Eg15Y3hiMjAwMjA5MDA3GghiazZiejkzZg%3D%3D SUN Yuhua, WANG Wei, LIU Siyang, et al. Genetic diversity research of mtDNA D-loop area from Chinese Myxocyprinus asiaticus[J]. Acta Genetica Sinica, 2002, 29(9): 787-790. (in Chinese) https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsQ0hJTmV3UzIwMjQxMTA1MTcxMzA0Eg15Y3hiMjAwMjA5MDA3GghiazZiejkzZg%3D%3D
[5] 杨金权, 胡雪莲, 唐文乔, 等. 长江口邻近水域刀鲚的线粒体控制区序列变异与遗传多样性[J]. 动物学杂志, 2008, 43(1): 8-15. doi: 10.3969/j.issn.0250-3263.2008.01.002 YANG Jinquan, HU Xuelian, TANG Wenqiao, et al. mtDNA control region sequence variation and genetic diversity of Coilia nasus in Yangtze River estuary and its adjacent waters[J]. Chin J Zool, 2008, 43(1): 8-15. (in Chinese) doi: 10.3969/j.issn.0250-3263.2008.01.002
[6] 彭珊, 代应贵. 濒危鱼类稀有白甲鱼清水江种群mtDNA D-loop序列多态性[J]. 水产学报, 2009, 33(2): 196-200. doi: 10.3321/j.issn:1000-0615.2009.02.004 PENG Shan, DAI Yinggui. Sequence polymorphism of mtDNA D-loop in the population of the endangered species Onychostom arara from the Qingshui River[J]. J Fish China, 2009, 33(2): 196-200. (in Chinese) doi: 10.3321/j.issn:1000-0615.2009.02.004
[7] 袁娟, 张其中, 李飞, 等. 铜鱼线粒体控制区的序列变异和遗传多样性[J]. 水生生物学报, 2010, 34(1): 9-19. doi: 10.3724/SP.J.1035.2010.00009 YUAN Juan, ZHANG Qizhong, LI Fei, et al. mtDNA control rigion sequence variation and genetic diversity of Coreius heterodon (Bleeker) in the upper and middle sections of the Yangtz River[J]. Acta Hydaobiologica Sinica, 2010, 34(1): 9-19. (in Chinese) doi: 10.3724/SP.J.1035.2010.00009
[8] 刘海林, 章群, 唐优良, 等. 黄渤海松江鲈鱼线粒体控制区结构与序列多态性分析[J]. 海洋通报, 2010, 29(3): 283-288. doi: 10.3969/j.issn.1001-6392.2010.03.008 LIU Hailin, ZHANG Qun, TANG Youliang, et al. Structure and genetic diversity of mtDNA D-Loop sequences among Trachidermus fasciatus stocks in Yellow Sea and Bohai Sea of China[J]. Mar Sci Bull, 2010, 29(3): 283-288. (in Chinese) doi: 10.3969/j.issn.1001-6392.2010.03.008
[9] 朱翠兰, 余为一, 陆依群, 等. 分布稻田的三种不同体色泥鳅组织LDH同工酶的研究[J]. 安徽农业大学报, 1996, 23(2): 155-159. doi: 10.1007/BF02009548 ZHU Cuilan, YU Weiyi, LU Yiqun, et al. Study on characterization of the electrophoretic patterns of LDH isoenzymes of Migurnus anguillicaudatus distributed over bog[J]. J Anhui Agric Univ, 1996, 23(2): 155-159. (in Chinese) doi: 10.1007/BF02009548
[10] 常重杰, 周家荣, 余其兴. 两种泥鳅不同群体遗传变异的RAPD分析[J]. 动物学报, 2001, 47(1): 89-93. doi: 10.3969/j.issn.1674-5507.2001.01.015 CHANG Chongjie, ZHOU Jiarong, YU Qixing. Genetic variation of two loach species revealed by RAPD analysis[J]. J Anim, 2001, 47(1): 89-93. (in Chinese) doi: 10.3969/j.issn.1674-5507.2001.01.015
[11] 曾柳根, 甘云飞, 王军花, 等. 鄱阳湖区泥鳅的微卫星DNA多态性分析[J]. 南昌大学学报: 理科版, 2008, 32(1): 84-88. doi: 10.3969/j.issn.1006-0464.2008.01.020 ZENG Liugen, GAN Yunfei, WANG Junhua, et al. Genetic diversity of loaches (Misgurnus anguillicaudatus) from Poyang Lake determined by microsatellites DNA[J]. J Nanchang Univ: Science Edition, 2008, 32(1): 84-88. (in Chinese) doi: 10.3969/j.issn.1006-0464.2008.01.020
[12] 单磊. 长江流域二倍体泥鳅野生群体遗传多样性的微卫星分析[D]. 武汉: 华中农业大学, 2009. 10.7666/d.y1812498 SHAN Lei. Population genetic diversity of natural diploid loach in the Yangtze River basin as revealed by microsatellite DNA analyses[D]. Wuhan: Huazhong Agricultural University, 2009. (in Chinese) 10.7666/d.y1812498
[13] PEREZ-ENRIQUZE R, TANIGUCHI N. Genetic structure of red sea bream (Pagrus major) population of Japan and the southwest Pacific using microsatellite DNA marks[J]. Fish Sci, 1999, 65(1): 23-30. doi: 10.2331/fishsci.65.23
[14] LIU Huanzhang, TZENG C S, TENG Huiyu. Sequence variations in the mitochondrial DNA control region and their implications for the phylogeny of the Cypriniformes[J]. Can J Zool, 2002, 80(6): 569-581 doi: 10.1139/z02-035
[15] SACCONE C, PESOLE G, SBISA E. Structral elements highly preserved during the evolution of the D-loop containing region in vertebrate mitochodrial DNA[J]. Mol Evol, 1987, 26(4): 205-211.
[16] 郑冰蓉, 张亚平, 肖蘅. 鲤属鱼类mtDNA控制区(D-环区)序列的变异性分析[J]. 水产学报, 2002, 26(4): 289-294. doi: 10.1007/BF02099853 ZHENG Bingrong, ZHANG Yaping, XIAO Heng. mtDNA control region of Cyprinus (D-loop region) sequence analysis of the variability[J]. J Fish China, 2002, 26(4): 289-294. (in Chinese) doi: 10.1007/BF02099853
[17] 刘焕章. 鱼类线粒体DNA控制区的结构与进化: 以鱼类为例[J]. 自然科学进展, 2002, 12(3): 266-270. doi: 10.3321/j.issn:1002-008X.2002.03.008 LIU Huanzhang. The structure and evolution of mitochondrial DNA control region of fish: a case study to bitterlings[J]. Prog Natural Sci, 2002, 12(3): 266-270. (in Chinese) doi: 10.3321/j.issn:1002-008X.2002.03.008
[18] BROWN G G. Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA[J]. J Mol Biol, 1986, 192(2): 503-511. doi: 10.1016/0022-2836(86)90272-x
[19] BROUGHTON R E, DOWLING T E. Length variation in mitochondrial DNA of the minnow Cyprinella spiloptera[J]. Genetics, 1994, 138(1): 179-190. doi: 10.1093/genetics/138.1.179
[20] LEE W J, CONROY J, HOWELL W H, et al. Structure and evolution of teleost mitochondrial control regions[J]. J Mol Evol, 1995, 41(1): 54-66. doi: 10.1007/BF00174041
[21] 张燕, 张鹗, 何舜平. 中国鲿科鱼类线粒体DNA控制区结构及其系统发育分析[J]. 水生生物学报, 2003, 27 (5): 463-467. doi: 10.3321/j.issn:1000-3207.2003.05.004 ZHANG Yan, ZHANG E, HE Shunping. Studies on the Structure of the control region of the Bagridae in China and its phylogentic significance[J]. Acta Hydaobiologica Sinica, 2003, 27(5): 463-467. (in Chinese) doi: 10.3321/j.issn:1000-3207.2003.05.004
[22] 唐琼英, 刘焕章, 杨秀平, 等. 沙鳅亚科鱼类线粒体DNA控制区结构分析及系统发育关系的研究[J]. 水生生物学报, 2009, 26(6): 645-653. http://ssswxb.ihb.ac.cn/article/id/6ab96820-95a7-44db-920c-235de6ab9a17 TANG Qiongying, LIU Huanzhang, YANG Xiuping, et al. Studies on the struction of the mitochondrial DNA control region and phylogenetic relationship of the subfamily Botiinae[J]. Acta Hydaobiologica Sinica, 2009, 26(6): 645-653. (in Chinese) http://ssswxb.ihb.ac.cn/article/id/6ab96820-95a7-44db-920c-235de6ab9a17
[23] SOUTHERN S O, SOUTHERN P J, DIZON A E. Molecular characterization of a cloned dolphin mitochondrial genome[J]. J Mol Evol, 1988/1989, 28(1/2): 32-42. doi: 10.1007/BF02143495
[24] SBISA E, TANZARIELLO F, REYES A, et al. Mammalian mitochondrial D-loop region structural analysis: identification of new conserved sequences and their functional and evolutionary implications[J]. Gene, 1997, 205(1/2): 125-140. doi: 10.1016/s0378-1119(97)00404-6
[25] SACCONE C, PESOLE G, SBISA E. The main regulatory region of mammalian mitochondrial DNA: structure-function model and evolutionary pattern[J]. J Mol Evol, 1991, 33(1): 83-91. doi: 10.1007/BF02100199
[26] 曾青兰, 刘焕章. 大口胭脂鱼线粒体DNA控制区序列的研究[J]. 湖北大学学报: 自然科学版, 2001, 23(3): 261-264. doi: 10.3969/j.issn.1000-2375.2001.03.021 ZENG Qinglan, LIU Huanzhang. Study on mitochondrial DNA control region of the Ictiobus cyprinellus[J]. J Hubei Univ: Natural Science, 2001, 23(3): 261-264. (in Chinese) doi: 10.3969/j.issn.1000-2375.2001.03.021
[27] 赵凯, 李军祥, 张亚平, 等. 青海湖裸鲤mtDNA遗传多样性的初步研究[J]. 遗传, 2001, 23(5): 445-448. https://xueshu.baidu.com/usercenter/paper/show?paperid=00d28978b80365611304219818059c47&site=xueshu_se ZHAO Kai, LI Junxiang, ZHANG Yaping, et al. Mitochondrial DNA diversity in naked carps in Qinghai-lake (Gymnocypris prezwalskii prezwalskii)[J]. Hereditas, 2001, 23(5): 445-448. (in Chinese) https://xueshu.baidu.com/usercenter/paper/show?paperid=00d28978b80365611304219818059c47&site=xueshu_se
[28] 赵刚, 周剑, 杜军, 等. 长薄鳅(Leptobotia elongata)线粒体DNA控制区遗传多样性研究[J]. 西南农业学报, 2010, 23(3): 930-937. doi: 10.3969/j.issn.1001-4829.2010.03.068 ZHAO Gang, ZHOU Jian, DU Jun, et al. Genetic diversity of elongate loach (Leptobotia elongata) inferred from mitochondrial DNA control region[J]. Southwest China J Agric Sci, 2010, 23(3): 930-937. (in Chinese) doi: 10.3969/j.issn.1001-4829.2010.03.068
[29] 杨博, 陈小勇, 杨君兴. 浪白鱼线粒体DNA控制区结构和种群遗传多样性分析[J]. 动物学研究, 2008, 29(4): 379-385. doi: 10.3321/j.issn:0254-5853.2008.04.005 YANG Bo, CHEN Xiaoyong, YANG Junxing. Structure of the mitochondrial DNA control region and population genetic diversity analysis of Anabarilius grahami (Regan)[J]. Zool Res, 2008, 29(4): 379-385. (in Chinese) doi: 10.3321/j.issn:0254-5853.2008.04.005
[30] 代应贵, 韩雪, 张晓杰. 小口白甲鱼都柳江种群mtDNA D环的序列变异及遗传多样性[J]. 动物学杂志, 2010, 45(2): 115-120. https://xueshu.baidu.com/usercenter/paper/show?paperid=a528705686b212f0d25207fde0a875cd&site=xueshu_se&hitarticle=1 DAI Yinggui, HAN Xue, ZHANG Xiaojie. Genetic diversity and sequence variation of mtDNA D-loop in the population of vulnerable species Onychostoma lini from the Duliu River[J]. China J Zool, 2010, 45(2): 115-120. (in Chinese) https://xueshu.baidu.com/usercenter/paper/show?paperid=a528705686b212f0d25207fde0a875cd&site=xueshu_se&hitarticle=1
[31] SONG Zhaobin, SONG Jun, YUE Bisong. Population genetic diversity of Prenant's schizothoracin, Schizothorax prenanti, inferred from the mitochondrial DNA control region[J]. Environ Biol Fish, 2008, 81(3): 247-252. doi: 10.1007/s10641-007-9197-6
-
期刊类型引用(26)
1. 温利红，张衡，方舟，陈新军. 鸢乌贼渔业资源研究进展. 水产科学. 2023(03): 527-537 . 
百度学术
2. 招春旭，吴文秀，邱星宇，周倍合，谢嘉仪，康斌，颜云榕. 南海不同海域鸢乌贼生长与死亡参数比较. 上海海洋大学学报. 2021(02): 294-300 . 百度学术
3. 谢嘉仪，张丽姿，吴文秀，周倍合，陈秋杰，招春旭，何雄波，徐军，颜云榕. 南沙群岛海域鸢乌贼摄食习性与营养生态位. 水产学报. 2021(12): 1993-2002 . 百度学术
4. 陈炫妤，陆化杰，王洪浩，何静茹，刘凯，陈新军. 西北印度洋鸢乌贼角质颚色素沉积特性分析. 动物学杂志. 2020(04): 468-476 . 百度学术
5. 朱凯，张立川，肖楚源，陈新军，林东明，朱俊磊. 南海鸢乌贼微型群雌性个体繁殖力研究. 渔业科学进展. 2020(06): 140-148 . 百度学术
6. 朱凯，姚吉祥，陈新军，刘维达，孙程婕，林东明. 南海鸢乌贼微型群肌肉和性腺组织能量积累及其分配. 上海海洋大学学报. 2020(06): 910-920 . 百度学术
7. 朱凯，王雪辉，杜飞雁，张鹏，邱永松. 南海中南部鸢乌贼中型群与微型群形态指标的分析比较. 中国海洋大学学报(自然科学版). 2019(01): 43-54 . 百度学术
8. 李敏，张鹏，张俊，张魁，陈作志. 南海鸢乌贼的遗传差异:种群分化还是种间分化. 中国水产科学. 2019(01): 133-140 . 百度学术
9. 朱凯，张立川，陈新军，姚吉祥，韩飞，林东明. 南海鸢乌贼中型群雄性个体肌肉和性腺组织能量积累及其分配. 热带海洋学报. 2019(04): 41-51 . 百度学术
10. 朱凯，张立川，陈新军，陆化杰，林东明，姚吉祥，马有成. 基于精荚数量对鸢乌贼中型群雄性个体有效繁殖力的研究. 动物学杂志. 2019(04): 517-528 . 百度学术
11. 黄佳兴，龚玉艳，徐姗楠，陈作志，张俊，于文明. 南海中西部海域鸢乌贼中型群和微型群的营养生态位. 应用生态学报. 2019(08): 2822-2828 . 百度学术
12. 刘玉，王雪辉，杜飞雁，刘必林，张鹏，刘梦娜，邱永松. 南海鸢乌贼耳石微量元素差异性分析. 南方水产科学. 2019(05): 15-24 . 本站查看
13. 朱凯，孙程婕，陈新军，林东明，刘海涛，古昊. 南海鸢乌贼中型群雌性个体肌肉和性腺组织能量积累及其分配. 海洋渔业. 2019(06): 641-651 . 百度学术
14. 刘玉，王雪辉，杜飞雁，刘必林，张鹏，刘梦娜，邱永松. 基于耳石微结构的南海鸢乌贼日龄和生长研究. 热带海洋学报. 2019(06): 62-73 . 百度学术
15. 龚玉艳，孔啸兰，杨玉滔，詹凤娉，张鹏，江艳娥，陈作志. 南海鸢乌贼微型群体的摄食习性研究. 海洋渔业. 2018(04): 394-403 . 百度学术
16. 徐亮，宁加佳，王雪辉，徐磊，刘玉，杜虹，杜飞雁. 脂类去除对南海鸢乌贼肌肉碳、氮稳定同位素分析的影响. 南方水产科学. 2018(04): 88-93 . 本站查看
17. 江淼，马胜伟，吴洽儿. 南海鸢乌贼资源探捕与开发. 中国渔业经济. 2018(02): 65-70 . 百度学术
18. 粟丽，陈作志，张鹏，李杰，王欢欢，黄佳兴. 2017年南海中南部渔场灯光罩网渔获物组成及渔获率时空分布. 南方水产科学. 2018(05): 11-20 . 本站查看
19. 王雪辉，邱永松，张鹏，杜飞雁. Natural mortality estimation and rational exploitation of purpleback flying squid Sthenoteuthis oualaniensis in the southern South China Sea. Chinese Journal of Oceanology and Limnology. 2017(04): 902-911 . 百度学术
20. 于鑫，单秀娟，李忠义，朱仁，陈云龙，金显仕. 渤海底拖网主要渔业生物类别时空分布的初步研究. 南方水产科学. 2017(02): 9-17 . 本站查看
21. 黄卉，杨丽芝，杨贤庆，李来好，郝淑贤，魏涯，王锦旭. 南海鸢乌贼墨汁多糖分离纯化及组分分析. 食品科学. 2017(24): 118-123 . 百度学术
22. 招春旭，陈昭澎，何雄波，邓玉淑，冯波，颜云榕. 基于耳石微结构的南海春季鸢乌贼日龄、生长与种群结构的研究. 水生生物学报. 2017(04): 884-890 . 百度学术
23. 朱凯，杜飞雁，王雪辉，邱永松，张鹏. 南海中南部鸢乌贼中型群群体结构. 应用生态学报. 2017(04): 1370-1376 . 百度学术
24. 粟丽，陈作志，张鹏. 南海中南部海域春秋季鸢乌贼繁殖生物学特征研究. 南方水产科学. 2016(04): 96-102 . 本站查看
25. 龚玉艳，詹凤娉，杨玉滔，张鹏，孔啸兰，江艳娥，陈作志. 南海鸢乌贼摄食习性的初步研究. 南方水产科学. 2016(04): 80-87 . 本站查看
26. 朱凯，王雪辉，张鹏，杜飞雁，邱永松. 南海南部鸢乌贼中型群与微型群形态学差异及其判别分析. 热带海洋学报. 2016(06): 82-88 . 百度学术
其他类型引用(8)
计量
- 文章访问数: 4695
- HTML全文浏览量: 161
- PDF下载量: 2727
- 被引次数: 34
粤公网安备 44010502001741号
