Comparative analysis of shell frame characteristics and mitochondrial 16S rRNA gene between wild and cultured mussels (Mytilus coruscus)
-
摘要:
了解厚壳贻贝 (Mytilus coruscus) 的生活环境对其生长的影响，对贻贝养殖业的高质量发展极为重要。以枸杞岛潮间带野生和浮筏养殖厚壳贻贝为研究对象，构建以壳长为基准的13个体框特征指标 (包括壳宽、壳高等)，比较两种厚壳贻贝的主要形态特征，并结合线粒体16S rRNA 基因的序列信息，分析其遗传关系与形态特征的关联性。结果表明，野生和养殖贻贝有10个体框特征指标存在显著性差异 (P<0.05)，野生的贝壳整体更宽厚、质量更大。应用逐步判别法筛选出的壳宽 (L1)、壳高 (L2)、壳顶至铰合部上端的距离 (L10)、壳顶至足丝孔的距离 (L12) 4个变量对野生和养殖厚壳贻贝的综合判别准确率达到94.9%。16S rRNA基因测序结果显示，野生厚壳贻贝的核苷酸多样性 (π) 为0.089，单倍型多样性 (Hd) 为0.894；养殖厚壳贻贝的π为0.087，Hd为0.682。野生与养殖厚壳贻贝的遗传多样性水平均较高，但系统发育树和单倍型网络图证明二者不存在显著的遗传分化，推测其形态差异可能主要受栖息环境 (波浪暴露强度、营养条件) 影响。
Abstract:Understanding the impact of the environment of Mytilus unguiculatus on its growth is extremely important for promoting high-quality development of mussel breeding industry. Taking the wild (From Gouqi intertidal area) and cultured (From mussel culture rafts) M. coruscus as research objects, we constructed 13 shell frame characteristics indexes based on shell length, including shell width and shell height to compare their main morphological differences, and sequenced their 16S rRNA gene to analyze the genetic and morphological associations between them. The results show that there were significant differences in ten shell frame characteristics between wild and cultured mussels (P<0.05), and the shells of wild mussels were generally wider, thicker and heavier. Four variables were selected by stepwise discrimination method: L1 (Shell width), L2 (Shell height), L10 (Distance from the top of shell to the upper end of hinge), L12 (Distance from the top of shell to the byssus orifice). The comprehensive discrimination accuracy of wild and cultured mussels was 94.9%. The 16S rRNA gene sequencing shows that the wild mussels nucleotide diversity (π) was 0.089, and haplotype diversity (Hd) was 0.894; cultured mussels' π was 0.087 and Hd was 0.682. The genetic diversity level of wild and cultured mussles were relatively high, but the phylogenetic tree and haplotype network prove that there was no significant genetic differentiation between them. It is speculated that their morphological differences might be mainly influenced by their habitat environments (Wave exposure, growth density, etc.).
-
Keywords:
- Mytilus coruscus /
- Morphological differences /
- 16S rRNA /
- Shell frame characteristics
-
大亚湾位于珠江口东侧，三面环山，为半封闭海湾，区域面积约600 km2，是我国南海北部重要的生物种质资源自然保护区之一。湾内岛屿众多，局部海域有海藻、珊瑚及人工鱼礁分布，生境复杂多样。目前，有关大亚湾鱼类资源与多样性评价方面的研究已有诸多报道。张敬怀等[1]基于大亚湾2004—2015年底拖网调查数据，研究分析了该海域物种组成及多样性时空变化；陈应华等[2]基于大亚湾大辣甲南人工鱼礁投放前后鱼类资源密度与优势种类组成的变化，初步评价了该海域人工鱼礁区建设效果；李娜娜等[3]运用声学方法对大亚湾杨梅坑水域生物资源量进行了科学评估；陈丕茂等[4]研究了大亚湾杨梅坑人工鱼礁区的渔业资源变动规律。上述研究主要围绕鱼类物种组成、资源动态、时空分布及多样性差异等方面，且研究区域主要集中在湾内人工鱼礁区及其毗邻水域等特殊生境。对于游泳生物群落格局及其生态效应方面的研究鲜有报道。
群落是实现海洋生态系统物质循环与能量传递的重要结构与功能单元[5]。基于群落时空格局及其变动规律的研究，能够很好地反映其所处生态系统的状态与变化趋势[6-8]。如王雪辉等[9]运用聚类与排序方法研究了大亚湾鱼类空间格局，并进一步探讨了鱼类群落与水环境因子的响应关系。大亚湾湾口连通湾内与湾外，平均水深约20 m。该区域流速较快，水交换较为频繁。涨潮时，形成逆时针环流；退潮时，潮流经中央列岛东侧流出，形成顺时针环流。此外，该区域营养盐和饵料生物较为丰富，是诸多海洋生物索饵、产卵和育肥的优良场所[3]。然而，有关湾口游泳生物组成、时空变动规律及群落结构特征的研究尚未见报道。
本研究基于2015年大亚湾湾口渔业资源与生态环境调查资料，分析了该海域游泳生物的组成、资源变动与群落结构特征。在此基础上，利用生态学多元统计方法，探究了导致该海域游泳生物群落时空异质结构的主要生物/非生物因子，以期增加对大亚湾游泳生物资源现状与群落结构特征的认识，为改善该水域生态状况、恢复和发展渔业资源提供基础资料。
1. 材料与方法
1.1 数据采集与处理
本研究所用游泳生物群落和水环境数据主要来源于大亚湾湾口区域，取样位点分布见图1。调查时间为2015年4月、8月、10月、12月，分别对应春、夏、秋、冬四季。调查渔船功率为135 kW，船长15.5 m，宽5.6 m，排水量50 t；通过虾拖网进行游泳生物采样，网长4.8 m，网口宽2.6 m，网目3 cm × 3 cm。各位点拖网1次，拖网时长为20~40 min，拖网时速为5.4~7.2 km·h–1。受地形及海况等条件的限制，不同航次和位点渔获强度存在一定差异，但各航次采样位点均分布在大亚湾湾口区域，分别代表不同季节调查区域游泳生物的资源组成特征，因此各季节不同位点渔获数据可用于大亚湾湾口游泳生物资源与群落季节异质特征的研究。为消除差异化拖网时长和航速的影响，各采样位点游泳生物资源密度均以0.5 h时长和6 km·h–1航速进行标准化处理。渔获样品均现场分类处理，主要分为鱼类、头足类、虾类、蟹类和虾蛄类，其中虾类、蟹类和虾蛄类统称甲壳类，难以鉴定的种类低温冷冻保存，带回实验室鉴定。对渔获个体数少于50的种类全部取样；渔获个体数大于50的种类，则根据个体体长范围进行分组，各体长组按数量百分比随机取样，取够50尾。测量并记录渔获样本体长/叉长/胴长和体质量，体长精确至1 cm，体质量精确至1 g。
各位点于拖网前，分表、底层进行水样采集。主要测定参数包括水温(WT)、溶解氧(DO)、盐度(Sal)、叶绿素a (Chl a)、水体透明度(Tra)和pH。其中Tra用塞氏盘测定，同时记录水深；WT、DO、Chl a、pH和Sal用便携式多参数水质分析仪YSI6600测定。
1.2 数据分析
1.2.1 物种相似性
本研究采用Jaccard物种相似性指数[10]，分析调查水域不同季度物种组成的变化。计算公式为：Sj = j/(a+b–j)。式中a为第一季度物种数；b为第二季度物种数；j为不同季度共有物种数。Jaccard指数1.0~0.75说明鱼类群落组成极为相似；0.75~0.5群落组成中等相似；0.5~0.25群落中等不相似；0.25~0群落极不相似。
1.2.2 优势度
利用相对重要性指数(IRI)[11]对大亚湾湾口各季度游泳生物优势种类组成进行分析。计算公式为：IRI = (N+W)×F。式中N为某一种类的尾数占渔获总尾数的百分比；W为某一种类的质量占渔获总质量的百分比；F为某一种类出现的位点数占总调查位点数的百分比(即出现频率)。将IRI≥500的物种定为优势种，100≤IRI<500为常见种，10≤IRI<100为一般种，IRI<10为少见种。
1.2.3 资源密度估算
游泳生物资源量密度D (kg·km–2) 根据扫海面积法估算，计算公式为：D = Y/[A×(1–E)]。式中Y为平均渔获率(kg·h–1)；A为每小时扫海面积(km2·h–1)；E为逃逸率，本研究中取0.5。调查区域各个航次不同游泳生物类群资源密度为各调查位点资源密度的算术平均值。
1.2.4 多样性评价
α多样性指数是度量生物多样性高低及其空间分布状态的主要数值指标，被广泛用于群落结构的变化与水生态系统环境质量评价。当调查水域不同种类及同种类个体间差异较大时，用生物量表示的多样性更接近种类间能量的分布[12]。因此，本研究中调查海域生物多样性的评价主要基于渔获生物量数据。该海域生物多样性的评价主要基于Shannon-Wiener多样性指数(H′)和Simpson优势集中度指数(C)，其计算公式为：
$$ \begin{aligned} H' = - \mathop \sum \nolimits_{i = 1}^s {P_i}{\rm{ln}}{P_i} \\ \;\;C = \mathop \sum \nolimits_{{i}}^{{s}} {\left( {{P_i}} \right)^2}\quad\;\, \end{aligned} $$ 式中Pi为і种所占渔获中生物量的百分比，S为样方内物种数。
1.2.5 群落结构与稳定性
依据调查海域各位点游泳生物聚类(Cluster)与非度量多维标度排序(NMDS)分析结果，了解不同季度游泳生物群落结构特征；并通过SIMPROF检验其结构差异的显著性。在此过程中，为平衡优势种和稀有种权重的影响，首先对各位点游泳生物丰度进行log (x+1)转化，再根据转化结果计算位点间Bray-Curtis相似性系数。在此基础上，通过相似性百分比(SIMPER)分析，判别对群落组间差异性和组内相似性贡献率较大的关键生物。ABC曲线通过生物丰度与生物量的位置关系，可分析生物群落在不同干扰程度下的生态效应[13]。
1.2.6 群落生态效应
基于游泳生物群落去趋势对应分析(DCA)结果，确定调查海域游泳生物群落生态效应研究的最适方法。根据DCA分析结果，四轴中最大值超过4，则选择单峰模型(典范对应分析CCA)较合适；小于3则线性模型(冗余分析RDA)较合适；两者之间则2种模型均适应。研究发现，调查海域游泳生物DCA分析结果，四轴中最大值为2.442，因此选择线性模型(RDA)分析环境变量与游泳生物之间的效应关系较为适宜。
应用RDA方法探究调查海域游泳生物群落与环境因子间的关系。首先对游泳生物群落数据进行log (x+1)转换和中心化处理，并剔除游泳生物群落中相对重要性指数(IRI)小于100的种类，以降低极端值的权重[14]。进而运用Monte Carlo偏置换检验，通过逐步迭代的方式，筛选对游泳生物群落空间特征的解释有显著贡献的环境变量，最终分析得出关键物种对特定生境的效应关系。
1.3 统计分析与作图
游泳生物资源密度、多样性指数及环境因子的季节差异，通过SPSS 21.0单因素方差分析(One-Way ANOVA)进行检验。多样性评价、群落结构与生态效应分析分别通过Primer 6.0和Canoco 4.5实现。本研究中绘图软件主要包括Primer 6.0、Canoco 4.5、ArcGIS 10.2.2及Origin 9.0。
2. 结果
2.1 物种组成
调查海域2015年春、夏、秋、冬季4个航次共捕获各类游泳生物177种，其中鱼类118种、蟹类32种、虾类16种、头足类7种、虾蛄类4种(附录1，详见http://dx.doi.org/10.12131/20180246的资源附件)。春、夏、秋、冬季游泳生物物种数分别为66、100、69、71种，可见调查水域夏季游泳生物物种数明显高于其他季度。调查海域各季度游泳生物组成见图2，鱼类在各类游泳生物组成中所占比重最高，其次为蟹类和虾类，其他种类所占比重相对较低。
![]()
图 2 大亚湾湾口2015年游泳生物百分比组成Figure 2. Percentage of nekton species at Daya Bay mouth in 2015本研究发现，调查海域2015年春、夏、秋、冬四季常见种类仅19种。根据Jaccard物种相似性研究结果(表1)，调查海域物种组成季节差异较大，季间物种相似性均在40%以下，其中秋季与夏季、冬季物种相似度相对较高(均为0.36)。整体上，调查海域游泳生物季间物种相似性处于中等不相似(0.5~0.25)和极不相似水平(0.25~0)。
表 1 大亚湾湾口2015年不同季节游泳生物Jaccard相似性指数Table 1. Jaccard similarity index of nekton at Daya Bay mouth in different seasons of 2015季节
season春季
spring夏季
summer秋季
autumn冬季
winter春季 spring 夏季 summer 0.27 秋季 autumn 0.29 0.36 冬季 winter 0.27 0.22 0.36 2.2 游泳生物资源密度季节变化
根据扫海面积法计算结果，调查海域春、夏、秋、冬季游泳生物资源量密度依次为1 536.4 kg·km–2、4 341.8 kg·km–2、1 634.1 kg·km–2、1 143.2 kg·km–2。各类游泳生物资源量密度详见表2。总体上鱼类资源密度相对较高，年均值为977.6 kg·km–2；其次为虾类和蟹类，年均质量密度分别为502.78 kg·km–2和439 kg·km–2；虾蛄类和头足类资源密度相对较低，年均值分别为175.63 kg·km–2和68.88 kg·km–2。单因素方差分析结果表明，夏季游泳生物资源量密度显著高于其他季节(F = 4.556, P = 0.01)。进一步分析发现，夏季鱼类资源密度显著高于其他季节(F = 10.323, P = 0.001)；而头足类(F = 0.346, P = 0.792)、虾蛄类(F = 2.559, P = 0.076)、虾类(F = 2.926, P = 0.052)和蟹类(F = 0.056, P = 0.982)资源密度季节差异不显著。
表 2 大亚湾湾口2015年游泳生物资源密度Table 2. Biomass density of nekton at Daya Bay mouth in 2015kg·km–2 季节
season头足类
Cephalopoda虾蛄类
Squilloidea虾类
Decapod蟹类
Crab鱼类
Fish春季 spring 41.2 225.5 67.6 399.4 802.7 夏季 summer 78.4 59.4 1 702.3 480.1 2 021.6 秋季 autumn 93.8 101.0 211.7 453.5 774.1 冬季 winter 62.1 316.6 29.5 423.0 312.0 2.3 优势度物种
基于相对重要性指数IRI研究结果(表3)，2015年春季，调查海域游泳生物优势种类主要包括二长棘鲷(Paerargyrops edita)等7种，累积数量和质量百分比分别为59.41%和40.61%。夏季优势种类主要有黄鳍马面鲀(Thamnaconus hypargyreus)等6种，累积数量和质量百分比分别为79.96%和49.42%。秋季优势种类主要有短吻鲾(Leiognathus brevirostris)等8种，累积数量和质量百分比分别为65.14%和62.26%。冬季优势种类主要包括拥剑梭子蟹(Portunus gladiator)和伪装关公蟹(Dorippe facchino)等7种，累积数量和质量百分比分别为59.93%和50.59%。综上可知：1)调查海域2015年各季节游泳生物优势种类存在较大差异。春季、夏季、秋季优势种类主要以鱼类为主，而冬季优势种类主要为甲壳类(蟹类、虾类、虾蛄类)。2)调查海域2015年各季节游泳生物个体差异较大。其中，夏季和春季优势种类数量百分比明显高于质量百分比，说明小型个体所占比重相对较高。
表 3 大亚湾湾口2015年游泳生物优势种类组成Table 3. Composition of dominant species of nekton at Daya Bay mouth in 2015季节
season种名
species数量比
N%质量比
W%频率
F%相对重要性指数
IRI春季 spring 二长棘鲷 Paerargyrops edita 31.6 15.6 100 4 717 口虾蛄 Oratosquilla oratoria 4.6 6.2 87.5 942 直额蟳 Charybdis truncatus 6.1 5.0 75 831 拟矛尾鰕虎鱼 Parachaeturichthys polynema 4.7 3.4 100 813 短脊鼓虾 Alpheus brevicristatus 6.8 2.1 87.5 780 隆线强蟹 Eucrate crenata 3.3 4.6 87.5 691 断脊口虾蛄 Oratosquillina interrupta 2.3 3.7 87.5 531 夏季 summer 黄鳍马面鲀 Thamnaconus hypargyreus 43.4 19.6 100 6 301 宽突赤虾 Metapenaeopsis palmensis 17.9 9.1 100 2 697 墨吉对虾 Penaeus merguiensis 6.3 13.7 57.1 1 142 侧线天竺鲷 Apogon lateralis 4.2 2.2 100 635 拥剑梭子蟹 Portunus gladiator 5.6 1.6 85.7 611 短吻鲾 Leiognathus brevirostris 2.6 3.4 100 593 秋季 autumn 短吻鲾 Leiognathus brevirostris 26.8 19.2 100 4 594 中线天竺鲷 Apogon kallopterus 13.1 3.8 87.5 1 472 红星梭子蟹 Portunus sanguinolentus 3.5 12.6 87.5 1 411 细条天竺鱼 Apogonichthys lineatus 10.3 3.7 87.5 1 228 墨吉对虾 Penaeus merguiensis 2.9 8.5 100 1 142 锈斑蟳 Charybdis feriatus 2.0 7.7 87.5 851 猛虾蛄 Harpiosquilla harpax 3.2 3.4 87.5 575 断脊口虾蛄 Oratosquillina interrupta 3.3 3.4 75 507 冬季 winter 拥剑梭子蟹 Portunus gladiator 19.0 5.0 100 2 394 伪装关公蟹 Dorippe facchino 15.2 11.8 87.5 2 364 断脊口虾蛄 Oratosquillina interrupta 7.4 9.4 100 1 682 杜氏枪乌贼 Uroteuthis duvauceli 6.7 10.4 75 1 280 猛虾蛄 Harpiosquilla harpax 4.2 6.5 75 804 口虾蛄 Oratosquilla oratoria 4.6 5.3 62.5 616 香港蟳 Charybdis hongkongensis 2.8 2.3 100 510 2.4 游泳生物群落特征
2.4.1 多样性水平
Shannon-Wiener多样性指数与Simpson优势集中度指数研究表明，调查海域不同季节生物多样性差异较小(图3)，但位点间差异明显(误差轴)。春、夏、秋、冬四季，Shannon-Wiener多样性指数分别为2.5、2.66、2.5、2.35；Simpson优势集中度指数分别为0.16、0.14、0.15、0.16。单因素方差分析结果表明，调查海域游泳生物Shannon-Wiener多样性指数(F = 0.597, P = 0.623)与Simpson优势集中度指数(F = 0.065, P = 0.978)季节差异均不显著。
![]()
图 3 大亚湾湾口2015年游泳生物多样性差异误差线表示地理空间差异Figure 3. Diversity variation of nekton at Daya Bay mouth in 2015Error bar represents geographical spatial difference.2.4.2 稳定性评价
湾口海域2015年各季节游泳生物丰度和生物量ABC曲线见图4。春、夏两季游泳生物丰度曲线在生物量曲线之上，w值分别为–0.112和–0.145，且第一丰度优势种起点较高(约30%~40%)，说明调查海域游泳生物主要以小型个体为主，且群落在春、夏两季受外界干扰强度较大，群落结构稳定性较差，应多予以保护。秋季和冬季渔获丰度和生物量累积优势度曲线基本处于交错状态，其w值基本为0，说明秋、冬两季调查海域游泳生物个体相对较大，群落受外界干扰相对较小，群落结构相对稳定。
![]()
图 4 大亚湾湾口2015年游泳生物群落ABC曲线Figure 4. ABC curve of nekton community at Daya Bay mouth in 20152.4.3 群落时空格局
Cluster (图5)与SIMPROF相似性检验显示，大亚湾湾口游泳生物群落结构季节差异显著(Pi = 9.3, P = 0.1%)；同一季节游泳生物群落空间差异较小，未表现出明显的空间格局。根据NMDS排序结果(图6)，2D Stress值为0.15，说明二维排序结果能较好地反映春、夏、秋、冬四季游泳生物群落结构特征。
![]()
图 5 大亚湾湾口2015年游泳生物群落聚类分析Figure 5. Cluster analysis of nekton community at Daya Bay mouth in 2015![]()
图 6 大亚湾湾口2015年游泳生物群落NMDS分析Figure 6. NMDS analysis of nekton community at Daya Bay mouth in 20152.5 游泳生物群落生态效应
2.5.1 生物效应因子
根据SIMPER相似性百分比分析结果(表4)，湾口海域游泳生物群落空间相似性为43.4%~51%，季间特异性为69.13%~80.77%。群落空间相似性贡献率前三位的种类，春季为二长棘鲷、拟矛尾鰕虎鱼(Parachaeturichthys polynema)和口虾蛄(Oratosquilla oratoria)；夏季为黄鳍马面鲀、宽突赤虾(Metapenaeopsis palmensis)和远海梭子蟹(Portunus pelagicus)；秋季为短吻鲾、墨吉对虾(Penaeus merguiensis)和红星梭子蟹(Portunus sanguinolentus)；冬季为断脊口虾蛄(Oratosquillina interrupta)、拥剑梭子蟹和香港蟳 (Charybdis hongkongensis)。群落季间特异性贡献率位于前三位的种类，春夏为黄鳍马面鲀、侧线天竺鲷(Apogon lateralis)和隆线强蟹(Eucrate crenata)；春秋为二长棘鲷、墨吉对虾和红星梭子蟹；春冬为二长棘鲷、拥剑梭子蟹和棘突猛虾蛄(Harpiosquilla raphidea)；夏秋为二长棘鲷、宽突赤虾和侧线天竺鲷；夏冬为二长棘鲷、黄鳍马面鲀和侧线天竺鲷；秋冬为短吻鲾、墨吉对虾和中线天竺鲷(Apogon kallopterus)。
表 4 大亚湾湾口2015年游泳生物群落空间相似性与季节差异性贡献率位于前3位的种类Table 4. Top three species contributing to spatial similarity and seasonal dissimilarity of nekton community at Daya Bay mouth in 2015% 种类
species空间相似性贡献率
similarity contribution rate in space季节差异性贡献率
dissimilarity contribution rate among seasonsA B C D A&B A&C A&D B&C B&D C&D 51 48.21 48.06 43.40 73.50 76.05 75.75 69.13 80.77 69.96 二长棘鲷 Parargyrops edita 8.68 5.64 0 0 0.55 3.28 3.57 2.53 2.3 0 拟矛尾鰕虎鱼 Parachaeturichthys polynema 6.62 1.67 2.5 0.67 1.02 1.24 2.12 0.96 0.93 1.39 口虾蛄 Oratosquilla oratoria 6.13 1.46 3.23 3.59 1.34 1.47 1.77 1.27 1.3 1.85 黄鳍马面鲀 Thamnaconus hypargyreus 0 8.04 3.93 0 3.37 1.87 0 1.91 3.19 2.13 宽突赤虾 Metapenaeopsis palmensis 5.3 6.23 0 1.79 1.33 1.76 1.23 2.81 1.91 1.25 远海梭子蟹 Portunus pelagicus 4.54 5.99 0.31 0.33 0.94 2 2.17 2.16 1.99 1.32 短吻鲾 Leiognathus brevirostris 2.53 5.45 9.75 1.48 1.38 1.97 1.69 0.76 1.74 3.13 墨吉对虾 Penaeus merguiensis 0 1.81 8.91 0.77 1.82 3.29 0.83 1.87 1.69 2.89 红星梭子蟹 Portunus sanguinolentus 0 0.26 6.52 2.39 0.74 3 1.63 2.13 1.26 2.38 断脊口虾蛄 Oratosquillina interrupta 5.18 0.37 3.66 9.64 1.38 1.42 1.1 1.54 1.71 1.67 拥剑梭子蟹 Portunus gladiator 0 2.36 2.21 9.25 1.59 1.23 2.58 1.3 1.03 1.6 香港蟳 Charybdis hongkongensis 3.51 0 0.14 7.04 1.48 1.63 1.13 0.32 1.67 2.06 中线天竺鲷 Apogon kiensis 0 0 5.92 0 0 2.39 0 2.02 0 2.73 侧线天竺鲷 Apogon lateralis 0 5.3 0 0 2.37 0 0 2.46 2.24 0 隆线强蟹 Eucrate crenata 4.88 0 0 1.33 1.92 2.33 1.9 0 0.77 1.13 棘突猛虾蛄 Harpiosquilla raphidea 4.29 0 0 0 1.74 2.11 2.3 0 0 0 注：A、B、C、D分别为春、夏、秋、冬季
Note: A, B, C and D represent spring, summer, autumn and winter, respectively.2.5.2 非生物效应因子
大亚湾湾口主要水环境要素2015年各季节调查结果见表5。区域年均WT约23.76 ℃、Sal为34.96、Tra为3.85 m、pH为8.32、DO质量浓度为5.79 mg·L–1、Chl a质量浓度为1.96 μg·L–1。根据单因素方差分析结果，各水环境因子季节差异显著(P < 0.01)。
表 5 大亚湾湾口2015年主要水环境因子季节差异Table 5. Seasonal variation of main water environment factors at Daya Bay mouth in 2015环境因子
environment factor春季
spring夏季
summer秋季
autumn冬季
winterF P 温度/℃ WT 21.56 26.72 27.39 19.35 1 047.81 0.001 盐度 Sal 34.89 35.74 34.81 34.42 134.89 0.001 透明度/m Tra 5.69 3.43 3.98 2.31 5.76 0.004 pH 8.22 8.25 8.38 8.43 21.57 0.001 溶解氧/mg·L –1 DO 6.91 4.58 5.42 6.23 53.67 0.001 叶绿素a/μg·L–1 Chl a 2.39 2.82 1.75 0.89 4.94 0.007 RDA研究发现(图7)，轴1和轴2能较好地解释研究海域游泳生物群落季节异质特征及其与环境因子的对应关系。蒙特卡洛统计结果，进一步检验了其相关显著性(F = 5.823，P = 0.001)。轴1和轴2对游泳生物群落季节异质性的解释率分别为0.302和0.166；游泳生物群落与环境对应关系的解释率分别为0.51和0.279。夏季，游泳生物群落主要沿轴1正方向分布，代表生物主要包括黄鳍马面鲀和侧线天竺鲷等；冬季，游泳生物群落主要沿轴1负方向分布，代表生物主要包括断脊口虾蛄和猛虾蛄等；秋季，游泳生物群落主要沿轴2正方向分布，代表生物主要包括中线天竺鲷、细条天竺鱼和侧线天竺鲷等，主要与WT和pH正相关；春季，游泳生物群落主要沿轴2负方向分布，代表生物主要包括二长棘鲷和拟矛尾鰕虎鱼等。环境因子中Sal、DO、WT和Chl a主要与轴1相关，相关系数分别为0.87、–0.77、0.63和0.41；pH和Tra主要与轴2相关，相关系数分别为0.67和–0.31。
![]()
图 7 大亚湾湾口2015年游泳生物群落与水环境因子RDA关系图sp1~sp21. 游泳生物种类；A1~D31. 各季度取样位点Figure 7. RDA ordination of nekton community and water parameters at Daya Bay mouth in 2015sp1–sp21. nekton species; A1–D31. sampling sites in each season3. 讨论
3.1 资源动态与多样性评价
渔业资源评估与动态监测旨在为渔业资源的保护和合理开发利用提供理论支持。本研究发现，大亚湾湾口夏季(8月)游泳生物密度明显高于其他季节，其中以鱼类和虾类资源较为丰富，而其他游泳生物类群(头足类等)季节差异不显著。这可能与我国南海伏季休渔政策相关(2015年调查海域休渔期为5月16日至8月1日)。诸多研究显示，伏季休渔政策能有效地缓解渔业资源衰退的局面，为鱼类等水生生物提供较稳定的繁殖、索饵与生长发育环境，进而达到资源保护与生态修复的目的。刘爱霞[15]研究发现，大亚湾2016年休渔期后渔业资源CPUE较休渔前明显上升。余景等[16]基于遥感数据定量评估了大亚湾2015年伏季休渔效果，研究发现休渔期后总渔获物CPUE从休渔前的0~10 kg·h–1上升至0~40 kg·h–1，且个体平均体长和体质量明显增加。此外，各季节游泳生物优势种类存在较大差异，该结果与张文超等[17]对北部湾洋浦海域鱼类群落结构的研究结果一致。春、夏、秋季第一优势种类分别为二长棘鲷、黄鳍马面鲀和短吻鲾，而冬季优势种类主要为拥剑梭子蟹等甲壳类生物(蟹类、虾类、虾蛄类)，该结果可能与调查海域游泳生物的繁殖、索饵及季度迁移习性相关[18]。
α多样性指数作为生物群落研究的重要指标之一[19]，常用来对比分析不同时期或水域生物群落多样性状况，探讨人类活动对鱼类群落多样性的影响及群落稳定性等[20]。如袁华荣等[21]通过Shannon-Wiener多样性指数、Pielou均匀度指数等对南海柘林湾鱼类多样性状况进行了评价，并分析了鱼类多样性与优势种类季节差异的相互关系。本研究中，夏季Shannon-Wiener多样性指数相对较高，但不同季节差异并不显著。湾口位置生境复杂，水体交换亦相对频繁，导致多样性指数空间差异较大(图3误差轴)。该差异经算术平均后被弱化，因此研究区域多样性指数季节差异并不显著。然而，当群落中的某一物种被与其生态功能相近的物种所代替时，这种群落结构的变化往往不能通过资源量或生物多样性指数等指标反映出来[22]。因此，单一的α多样性指数并不能有效地描述特定区域范围内生物群落的组成与结构特征[23]。
3.2 群落结构与稳定性评价
β多样性应用于生物群落时空格局研究，较α多样性指数能反映更多的群落信息，不仅描述生境内物种数量，同时也考虑到不同种类彼此之间的位置关系[24]。本研究采用生态学多元统计的方法，旨在研究大亚湾湾口海域游泳生物群落季节变化特征。由于湾口区域面积较小，调查过程中受天气、船舶通航及渔民捕捞作业等客观条件的影响，不同季节调查位点难以一一对应。各季节生物学采样位点均按照随机抽样的方式，均匀分布在湾口海域，整体上能反映该海域各个季节游泳生物组成特点，且各个季节不同调查位点渔获组成均进行了标准化处理，以保证其具有可比性。因此，本研究各个季节均采用随机抽样调查的方式，理论上是可行的。Cluster与NMDS研究表明，调查水域游泳生物群落季节差异显著(Pi = 9.3，P = 0.1%)，群落季间相似性均小于35%。诸多研究发现，群落中优势种类繁殖和索饵习性的差异(如多数鱼类的繁殖主要集中在春夏两季)，是导致群落季节异质特征的主要影响因子[25-26]。根据相对性指数计算结果，调查海域游泳生物优势种类季节差异显著。SIMPER相似性百分比进一步分析发现(表4)，二长棘鲷、黄鳍马面鲀和侧线天竺鲷为夏季和冬季游泳生物群落中特异性优势物种；二长棘鲷、墨吉对虾和红星梭子蟹为春季和秋季游泳生物群落中特异性优势物种。相关研究指出，优势种为海洋生态系统中的关键因子，对于维持生态系统的稳定以及渔业资源的可持续发展具有十分重要的意义[27]。根据上述研究结果，建议针对大亚湾湾口游泳生物群落中季节特异性优势物种，制定科学合理的保护措施(如限额捕捞、人工鱼礁建设等)，以促进该海域游泳生物资源的发展与群落结构的稳定。
ABC曲线中丰度和生物量优势度曲线的相对位置能反映游泳生物群落受外界环境的干扰程度[28]。有关研究显示，游泳生物群落稳定性表现出明显的季节变化特征。受水温、水位、水流等环境因子的影响，夏季游泳生物生理机能明显恢复，摄食增加、生长快速、繁殖旺盛且群落中小型个体所占比重大幅上升。近年来，随着大亚湾周边工业、港口建设、滨海旅游和规模化养殖等的迅速发展，该海域生态环境发生了较大变化，对鱼类等水生生物群落产生了不小的扰动，具体表现为渔业资源量下降、生物多样性降低、群落稳定性下降等[29]。本研究中，春季和夏季游泳生物丰度优势度曲线均在生物量优势度曲线之上，说明游泳生物群落以个体较小、生长快速、性成熟较早的r选择物种为主。另有研究指出，ABC曲线对优势群体的丰度和生物量变化比较敏感，优势种的丰度和生物量决定了生物量优势度曲线和丰度优势度曲线的位置[30]。本研究发现，夏季和春季第一优势种类数量百分比明显高于质量百分比，而秋季和冬季优势种类数量和质量百分比差异相对较小(表3)。综上可知，春、夏时期调查海域游泳生物群落受外界环境扰动的影响较大，群落结构处于不稳定状态，应采取有效措施，如伏季休渔等，为游泳生物提供良好的繁育环境。
3.3 游泳生物群落结构与水环境效应关系
有关研究指出，生境的变迁在一定时空范围内能导致群落异质结构的形成[31]。生态学多元统计分析(RDA、CCA等)应用于生物群落生态效应研究，能同时实现定量与定性分析的要求，因而逐渐成为群落生态研究的重要方法[32]。本研究中根据游泳生物群落与重要环境参数冗余分析结果，第一、二轴对游泳生物群落与环境关系的解释率分别为0.51和0.279。夏季和冬季游泳生物群落沿轴1分布，主要环境效应因子为Sal、DO、WT和Chl a；秋季和春季游泳生物群落沿轴2分布，主要环境效应因子为pH和Tra。由此可知，夏季游泳生物群落中二长棘鲷、黄鳍马面鲀和侧线天竺鲷的繁殖生长，对Sal、WT和Chl a等环境因子需求较高；冬季和秋季游泳生物群落中甲壳类(墨吉对虾、红星梭子蟹等)对pH耐受性较强，而不利于二长棘鲷等鱼类的种群增长。国外相关研究指出，pH水平与鱼类群落的组成和分布关系十分密切[33]。Ribeiro等[34]研究表明，水温在时间尺度上对鱼类群落影响最为显著，水温的季节变动直接导致了鱼类群落的季节差异[35]。另有研究指出，透明度在鱼类群落的时空分布中也起着十分重要的作用，其影响主要体现在索饵与避敌2个方面[36]。然而，生物群落与生态环境之间的效应关系因时因地而异，如Lunven等[37]研究指出，春季斯开湾浮游生物群落主要受磷酸盐调节，而夏季受硝酸盐影响较大。张文超等[17]研究发现，洋浦湾口海水的频繁交换，使得该海域水环境复杂，进而对该海域鱼类群落结构产生影响。因此，需尽可能系统全面地考虑潜在环境因子(如水深、海流、底质、人类活动等)的影响，以提高环境因子对生物群落时空异质特征的解释率。
本研究综合大亚湾湾口游泳生物资源动态、群落结构及其生态效应等多方面调查结果发现，调查海域主要水环境参数(Sal、WT、DO、pH、Chl a和Tra)及游泳生物中优势种类(二长棘鲷、黄鳍马面鲀、短吻鲾等)生理活动的季节性差异，是导致该海域游泳生物群落季节异质特征的主要影响因子。
-
![]()
图 2 野生与养殖厚壳贻贝壳长和壳质量的关系
Figure 2. Relationship between shell length and shell mass of wild and cultured M. coruscus
![]()
图 3 野生与养殖厚壳贻贝体框特征差异
Figure 3. Differences in shell frame characteristics between wild and cultured M. coruscus
![]()
图 5 野生与养殖厚壳贻贝主成分分析得分图
Figure 5. Principal component scores of wild and cultured M. coruscus
![]()
图 6 野生与养殖厚壳贻贝16S rRNA基因的单倍型网络图
Figure 6. Haploid network diagram of 16S rRNA gene in wild and cultured M. coruscus
![]()
图 7 邻接距离法构建的厚壳贻贝16S rRNA 系统发育树
注：YS代表野生厚壳贻贝，YZ代表养殖厚壳贻贝。
Figure 7. Neighbour-joining tree based on 16S rRNA gene of M. coruscus
Note: YS represents wild M. coruscus, and YZ represents cultured M. coruscus.
表 1 野生与养殖厚壳贻贝体框特征参数
Table 1 Shell frame characteristics parameters of wild and cultured M. coruscus
体框架特征
Characteristics代号
Code野生群体 Wild population 养殖群体 Cultured population Mean±SD CV/% Mean±SD CV/% W/L L1 0.323±0.021 6.52 0.374±0.025 6.76 H/L L2 0.490±0.035 7.10 0.535±0.034 6.36 AB/L L3 0.531±0.025 4.79 0.553±0.037 6.63 BC/L L4 0.198±0.034 16.95 0.199±0.042 21.22 CD/L L5 0.446±0.027 6.08 0.473±0.039 8.27 DE/L L6 0.447±0.047 10.49 0.460±0.049 10.64 AE/L L7 0.577±0.057 9.94 0.539±0.051 9.42 OA/L L8 0.683±0.039 5.67 0.729±0.046 6.37 OB/L L9 0.428±0.041 9.60 0.482±0.034 7.04 OC/L L10 0.306±0.023 7.67 0.337±0.028 8.37 OD/L L11 0.341±0.0.36 10.46 0.339±0.039 11.65 OE/L L12 0.266±0.044 16.55 0.299±0.037 12.55 OF/L L13 0.134±0.013 10.07 0.138±0.014 10.52 
下载: 导出CSV
表 2 野生与养殖厚壳贻贝体框特征之间的主成分分析
Table 2 Principal component analysis of shell frame characteristics between wild and cultured M. coruscus
方差解释及载荷
Variance explantion and loading主成分 Principal component PC1 PC2 PC3 PC4 总方差解释 Total variance explanation 特征根 Eigenvalue 4.217 1.823 1.473 1.262 差异解释率 Proportion of variance/% 32.441 14.022 11.328 9.707 累积解释率 Cumulative proportion of variance/% 32.441 46.463 57.791 67.498 因子载荷 Loading L1 0.680 0.322 0.281 0.083 L2 0.763 0.001 0.142 0.404 L3 0.502 0.167 −0.255 0.416 L4 0.308 −0.548 0.165 0.497 L5 0.292 0.633 0.078 −0.482 L6 0.544 −0.533 0.130 −0.453 L7 −0.640 0.417 −0.180 0.314 L8 0.675 0.186 −0.438 0.025 L9 0.858 0.055 −0.066 0.042 L10 0.545 0.585 0.004 −0.017 L11 −0.184 0.188 0.776 −0.037 L12 0.702 −0.329 −0.043 −0.309 L13 0.142 0.092 0.652 0.192
下载: 导出CSV
-
[1] 常抗美, 吴剑锋. 厚壳贻贝人工繁殖技术的研究[J]. 南方水产, 2007, 3(3): 26-30. [2] 张义浩, 赵盛龙. 嵊山列岛贻贝养殖种类生长发育调查[J]. 浙江海洋学院学报 (自然科学版), 2003, 22(1): 67-73. [3] 吴剑锋, 顾庆庭, 常抗美, 等. 厚壳贻贝人工育苗中幼虫合理培养密度的探讨[J]. 动物学杂志, 1982, 17(1): 19-22. [4] SÉBASTIEN L G, HAMILTON D J, BARBEAU M A. Effect of mussel density and size on the morphology of blue mussels (Mytilus edulis) grown in suspended culture in Prince Edward Island, Canada[J]. Aquaculture, 2005, 249(1/2/3/4): 265-274.
[5] HOBDAY A. Body-size variation exhibited by an intertidal limpet: influence of wave exposure, tidal height and migratory behavior[J]. J Exp Mar Biol Ecol, 1995, 189(1/2): 29-45.
[6] BOULDING E G, HOLST M, PILON V. Changes in selection on gastropod shell size and thickness with wave-exposure on northeastern pacific shores[J]. J Exp Mar Biol Ecol, 1999, 232(2): 217-239. doi: 10.1016/S0022-0981(98)00117-8
[7] 刘硕博, 唐祖蓉, 申望, 等. 基于F型mtDNA D-Loop的厚壳贻贝 (Mytilus coruscus) 群体遗传多样性研究[J]. 海洋与湖沼, 2019, 50(2): 355-364. [8] 管成林, 李继姬, 郭宝英. 基于线粒体CO I基因比较分析野生与养殖厚壳贻贝种群的遗传多样性[J]. 安徽农业科学, 2012, 40(13): 7660-7662. [9] 刘二田, 陈一, 战欣, 等. 企鹅珍珠贝微卫星标记与生长性状的相关性研究[J]. 水生生物学报, 2023, 47(2): 332-338. [10] YUAN Y, KANNO M, KIJIMA A. Genetic diversity of wild populations of Chlamys farreri in Japan and their genetic relationship with cultured stocks in China[J]. Aquaculture, 2012, 370: 109-122.
[11] HE X, WU F C, QI H G, et al. Whole-genome resequencing reveals the single nucleotide polymorphisms associated with shell shape in Crassostrea gigas[J]. Aquaculture, 2022, 547: 737502. doi: 10.1016/j.aquaculture.2021.737502
[12] 徐义平, 许会宾, 金凯, 等. 浙江乐清湾缢蛏的形态和遗传多样性[J]. 上海海洋大学学报, 2017, 26(1): 31-37. doi: 10.12024/jsou.20151001575 [13] 程亮, 徐善良, 刘飞, 等. 厚壳贻贝性腺不同发育时期肥满度与生化成分分析[J]. 海洋学研究, 2013, 31(4): 68-73. [14] 白晓倩, 杨阳, 邹李昶, 等. 嵊泗列岛海域三种贻贝贝体框架特征的差异[J]. 海洋与湖沼, 2014, 45(5): 1078-1084. [15] 李思发, 李晨虹, 李家乐. 尼罗罗非鱼品系间形态差异分析[J]. 动物学报, 1998(4): 75-82. [16] SUDHIR K, GLEN S, KOICHIRO T. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Mol Biol Evol, 2016, 33(7): 1870-1874. doi: 10.1093/molbev/msw054
[17] GREFSRUD E S, OIVIND S. Comparison of shell strength in wild and cultured scallops (Pecten maximus)[J]. Aquaculture, 2006, 251(2/3/4): 306-313.
[18] 顾志峰, 王嫣, 石耀华, 等. 马氏珠母贝两个不同地理种群的形态性状和贝壳珍珠质颜色比较分析[J]. 渔业科学进展, 2009, 30(1): 79-86. [19] 张义浩. 浙江沿海贻贝种类形态比较研究[J]. 渔业经济研究, 2009(2): 14-20. [20] FRANK E. ANDERSON. Phylogeny and historical biogeography of the loliginid squids (Mollusca: Cephalopoda) based on mitochondrial DNA sequence data[J]. Mol Phylogenet Evol, 2000, 15(2): 191-214. doi: 10.1006/mpev.1999.0753
[21] 孙明媛, 朱锐, 孙晓晴, 等. 基于CO I和16S rRNA的库页岛马珂蛤遗传多样性和分子系统进化研究[J]. 中国水产科学, 2018, 25(4): 891-901. [22] 闫永斌, 葛玉双, 程起群, 等. 基于线粒体16S rRNA和CO I序列的江苏启东海域4种贝类遗传多样性分析[J]. 海洋渔业, 2022, 44(2): 201-208. doi: 10.3969/j.issn.1004-2490.2022.02.007 [23] 谭月, 邱明, 李金龙, 等. 基于线粒体16S rRNA基因序列的中日朝菲律宾蛤仔野生群体遗传多样性分析[J]. 大连海洋大学学报, 2021, 36(1): 127-134. [24] SHEN Y B, LI J L, FENG B B. Genetic analysis of cultured and wild populations of Mytilus coruscus based on mitochondrial DNA[J]. Zool Res, 2009, 30(3): 240-246. doi: 10.3724/SP.J.1141.2009.03240
[25] 叶莹莹, 徐梅英, 郭宝英, 等. 我国东南沿海厚壳贻贝 (Mytilus coruscus) 4个群体线粒体16S rRNA序列及遗传结构分析[J]. 海洋与湖沼, 2012, 43(2): 376-381. doi: 10.11693/hyhz201202028028 [26] MONTAUDOUIN X D, LUCIA M, BINIAS C, et al. Why is Asari (=Manila) clam Ruditapes philippinarum fitness poor in Arcachon Bay: a meta-analysis to answer?[J]. Estuar Coast Shelf Sci, 2016, 179(20): 226-235.
[27] HAN Z Q, LI Q, LIU S K, et al. Genetic variability of an orange-shell line of the Pacific oyster Crassostrea gigas during artificial selection inferred from microsatellites and mitochondrial CO I sequences[J]. Aquaculture, 2019, 508: 159-166. doi: 10.1016/j.aquaculture.2019.04.074
[28] SMKA B, AMSB C. Molecules and morphology reveal 'new' widespread North American freshwater mussel species (Bivalvia: Unionidae)[J]. Mol Phylogenet Evol, 2019(138): 182-192.
[29] 祁剑飞, 宁岳, 巫旗生, 等. 菲律宾蛤仔5个地理群体形态差异分析[J]. 应用海洋学学报, 2020, 39(2): 189-194. doi: 10.3969/J.ISSN.2095-4972.2020.02.005 [30] 董建宇, 胡成业, 杨晓龙, 等. 基于GLM模型Logistic回归的菲律宾蛤仔地理群体判别[J]. 海洋与湖沼, 2018, 49(6): 1318-1324. doi: 10.11693/hyhz20180100002 [31] 仉佳文, 唐君玮, 张远帆, 等. 皱纹盘鲍 (Haliotis discus hannai) 不同群体贝壳形态和相对生长变异研究[J]. 海洋与湖沼, 2022, 53(6): 1539-1547. doi: 10.11693/hyhz20220300050 [32] MOESER G M, CARRINGTON E. Seasonal variation in mussel byssal thread mechanics[J]. J Exp Biol, 2006, 209(10): 1996-2003. doi: 10.1242/jeb.02234
[33] STEFFANI C N, BRANCH G M. Growth rate, condition, and shell shape of Mytilus galloprovincialis: responses to wave exposure[J]. Mar Ecol Prog Ser, 2003, 246(6): 197-209.
[34] AKESTER R J, MARTEL A L. Shell shape, dysodont tooth morphology, and hinge-ligament thickness in the bay mussel Mytilus trossulus correlate with wave exposure[J]. Can J Zool, 2000, 78(2): 240-253. doi: 10.1139/z99-215
[35] FRCHETTE M, DESPLAND E. Impaired shell gaping and food depletion as mechanisms of asymmetric competition in mussels[J]. Ecoscience, 1999, 6(1): 1-11. doi: 10.1080/11956860.1999.11952200
[36] 杨冠林. 典型贻贝浮筏式养殖海域物理-生态-贝类生长耦合数值模型的构建及应用[D]. 上海: 上海海洋大学, 2022: 2-8. [37] 温英. 贻贝浮筏养殖设施水流阻力的量化及其生态效应研究[D]. 上海: 上海海洋大学, 2021: 16-20. [38] FRÉCHETTE M, BOURGET E. Energy flow between the pelagic and benthic zones: factors controlling particulate organic matter available to an intertidal mussel bed[J]. Can J Fish Aquat Sci, 1985, 42(6): 1158-1165. doi: 10.1139/f85-143
[39] WATANABE S, KATAYAMA S. Relationships among shell shape, shell growth rate, and nutritional condition in the Manila clam (Ruditapes philippinarum) in Japan[J]. J Shellfish Res, 2014, 29(2): 353-359.
[40] 李杰. 日本镜蛤 (Dosinorbis japonica) 数量性状、营养成分及幼体生态学研究[D]. 大连: 大连海洋大学, 2018: 23-30. [41] 张学开, 范超, 霍忠明, 等. 不同群体长牡蛎壳型相关性状的比较分析[J]. 中国海洋大学学报 (自然科学版), 2021, 51(1): 1-9. [42] 江勇波. 企鹅珍珠贝摄食生理能量学研究[D]. 海口: 海南大学, 2014: 12-14. -
期刊类型引用(6)
1. 陈志，翁朝红，邹丽珍，周洪磊，蒋奕雄，方垂弘，李裕红，关杰耀. 中国鲎人工驯养繁殖技术初步研究. 渔业研究. 2023(03): 246-253 . 
百度学术
2. 钟赟民，翁朝红，王志婷，李佳怡，谢仰杰. 中国鲎室内人工繁育和鲎卵、幼体常见病原鉴定及其防治. 中国水产科学. 2023(06): 799-811 . 百度学术
3. 陈晓海，李银康，应紫薇，郭清扬，颉晓勇. 人为干扰栖息地退化对潮间带中华鲎幼鲎种群数量的影响. 生态学报. 2023(16): 6824-6832 . 百度学术
4. YING Ziwei，BAO Yuyuan，LI Yinkang，YE Guoling，ZHANG Shuhuan，XU Peng，ZHU Junhua，XIE Xiaoyong. Impact of Different Diets on Adult Tri-Spine Horseshoe Crab, Tachypleus tridentatus. Journal of Ocean University of China. 2022(03): 541-548 . 必应学术
5. 郭清扬，左娅，林吴颖，肖晓波，赵心忆，李银康，应紫薇，周传江，颉晓勇. 北部湾中华鲎（Tachypleus tridentatus）幼鲎生态位及种间竞争分析. 海洋学报. 2022(10): 109-118 . 百度学术
6. 应紫薇，李银康，颉晓勇. 中华鲎(Tachypleus tridentatus)稚鲎对不同比例泥沙类型选择行为特征. 海洋与湖沼. 2022(05): 1242-1249 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 315
- HTML全文浏览量: 99
- PDF下载量: 71
- 被引次数: 7
粤公网安备 44010502001741号
