西白令海阿拉斯加狭鳕矢耳石形态特征研究

段咪1、2,魏联1、2,朱国平1、2、3

(1.上海海洋大学 海洋科学学院,上海 201306;2.上海海洋大学 大洋渔业资源可持续开发教育部重点实验室极地海洋生态系统研究室,上海 201306;3.国家远洋渔业工程研究中心,上海 201306)

摘要:为研究阿拉斯加狭鳕Gadus chalcogrammus耳石与鱼体生长的关系,对采自西白令海水域53尾阿拉斯加狭鳕样本的矢耳石外部形态进行了观察与测量。结果表明:阿拉斯加狭鳕矢耳石外形较大,呈近椭圆形,基叶和翼叶明显,主间沟不明显;体长为35~43 cm的个体耳石长和耳石宽分别为12.93~17.05 mm和5.70~7.78 mm;左、右矢耳石形态无显著性差异(P>0.05),随着个体体长的增加,耳石绝对尺寸增加,但相对尺寸减小,耳石的形态结构比例基本不变;通过对耳石外部形态9个指标主成分分析,可以用耳石面积、背长和耳石长来描述阿拉斯加狭鳕耳石的主要外部形态特征,耳石主要形态指标与体长、湿体质量的关系显示,耳石长、耳石宽、背长、基叶长、翼叶长、耳石面积和耳石周长与鱼体体长呈指数关系;而背宽与鱼体体长呈幂函数关系;耳石长、背宽、背长、基叶长、耳石面积和耳石周长与鱼体湿质量呈幂函数关系,耳石宽、翼叶长与鱼体湿质量呈指数关系。本研究结果可为阿拉斯加狭鳕年龄鉴定和种群鉴别提供基础数据。

关键词: 阿拉斯加狭鳕;耳石;形态特征;主成分分析

阿拉斯加狭鳕Gadus chalcogrammus也称狭鳕鱼、明太鱼,隶属于鳕形目Gadiformes、鳕鱼科Gadidae、狭鳕属Gadus,主要分布于北太平洋,为当今世界产量最高的鱼种之一[1]。该鱼种生长速度相对较快,为白令海生态系统的主要组成部分[2]。阿拉斯加狭鳕属于冷水性鱼类,幼鱼主要摄食桡足类,成鱼摄食磷虾类、端足类和狭鳕幼鱼等[3]。中国狭鳕每年加工量约50万t,主要用于生产冷冻鳕鱼片、鳕鱼块和鱼糜制品等[4]

耳石是存在于硬骨鱼类内耳膜迷路内的一种钙化组织,主要由碳酸钙组成,平时起到平衡和听觉作用[5]。耳石具有较高的物种特异性[6],其形态大小、微结构特征因种类不同而存在差异,因此,鱼类耳石不仅可作为鉴定鱼类年龄与分析生长的重要材料[7],还可用于种类或种群鉴定[8-10]

迄今为止,有关阿拉斯加狭鳕的研究多与生态及资源分布[11]和资源量评估[2]有关。Blood等[12]对阿拉斯加狭鳕早期胚胎发育进行了研究,Jung等[13]分析了20世纪70年代白令海峡东南岸海域狭鳕鱼卵和仔稚鱼的生态学特征。关于耳石仅见Yoklavich等[14]通过耳石对阿拉斯加狭鳕幼鱼年龄和生长进行了研究,目前,尚未见国内学者对北极鱼类耳石形态特征进行专门研究。为此,本研究中对西白令海域阿拉斯加狭鳕进行生物学测定,并对其耳石外部形态特征进行测定和分析,旨在揭示阿拉斯加狭鳕矢耳石的形态特征及各形态指标与阿拉斯加狭鳕体长和湿体质量的关系,为其年龄鉴定及种群鉴别提供数据资料。

1 材料与方法

1.1 材料

阿拉斯加狭鳕样本由大型拖网加工船“连兴海”轮采集于西白令海俄罗斯专属经济区水域,采集时间为2016年9月25日,采集地点为61°25′N、178°15′W。采集的样本经冷冻保藏运回实验室进行后续分析。

1.2 方法

1.2.1 生物学测定和耳石处理 样本带回实验室

解冻后进行常规生物学测量,包括体长(L,由吻端至尾鳍基部之间的距离)和湿体质量(W)等,其中体长测量精确到0.1 cm,湿体质量测量精确到0.1 g。阿拉斯加狭鳕体长为35~43 cm,湿体质量为296.0~613.1 g。使用镊子和解剖刀由头部平衡囊内摘取耳石,将摘取后的耳石进行编号,用酒精清洗耳石表面的黏膜和有机质,之后用超纯水冲洗干净,进行干燥处理。干燥后用信封袋保存,最后得到完整的耳石样本共45对。采用微量天平(CPA225D)称量耳石质量(精确到0.000 1 g)。

1.2.2 耳石图像采集与测量 将矢耳石凸面朝上、凹面向下,在OLYMPUS解剖镜下(10倍率)对耳石进行拍照,图像采集后对耳石进行打磨,以确定其核心(图1),核心确定后,采用Digimizer图像分析软件测量耳石形态学参数值,精确到0.01 mm。配对样本t检验显示,左、右矢耳石形态间无显著性差异(P>0.05),因此,试验中统一使用左矢耳石。结合阿拉斯加狭鳕矢耳石外部形态特征及参考李凤莹等[15]和魏联等[16]的研究方法,试验中使用了以下8组形态学测量指标,即耳石长(OL,mm),耳石宽(OW,mm),背长(DL,mm),背宽(DW,mm),翼叶长(ARL,mm),基叶长(RL,mm),周长(P,mm),面积(S,mm2)(图2)。

图1 阿拉斯加狭鳕矢耳石打磨后平面图
Fig.1 Sagittal plane of otolith with ground surface in Alaska pollock Gadus chalcogrammus

图2 阿拉斯加狭鳕矢耳石形态参数示意图
Fig.2 Sketch photo on morphological parameters of otolith in Alaska pollock Gadus chalcogrammus

1.2.3 耳石的形态描述 参照李凤莹等[15]、魏联等[16]、李史民等[17]、Gaemers[18]和Messiech[10]的相关研究,结合阿拉斯加狭鳕耳石外部形态特征,使用如下术语描述耳石形态:基叶,耳石腹部前端的突起;翼叶,耳石背部前端的突起;主间沟,基叶和翼叶间的凹陷部分;辐射状条纹,耳石侧表面的辐射状纹路;核心为耳石的中心区域。

1.3 数据处理

(1) 首先利用Shapiro-Wilk检验对数据序列进行正态性检验。若符合正态分布,则利用配对t检验分析两组数据是否存在显著性差异(P<0.05),若不符合正态分布,则利用配对Kolmogorov-Smirnov检验。

(2) 采用幂函数、指数函数、对数函数和线性函数等模型,拟合阿拉斯加狭鳕耳石形态特征参数与体长和湿体质量间的关系,并检验其显著性,利用赤池信息准则(AIC)选取拟合函数类型,计算公式如下:

AIC=n×ln(RSS/n)+2k

其中:n为样本数量;k为方程中参数常数的数量;RSS为残差平方和。多种生长模型中,取AIC值最小的模型为最适生长模型[19-20]

(3) 计算耳石各个形态学特征参数值与体长的比值,分析耳石各形态学参数的相对尺寸。

(4) 计算耳石各形态特征参数值与耳石长度之比,分析不同体长下耳石各部分的形态变化。

(5) 用SPSS 19.0软件对阿拉斯加狭鳕耳石形态特征参数进行主成分分析,得到用于描述阿拉斯加狭鳕耳石外型的主要特征参数。

2 结果与分析

2.1 耳石外型的测量

Shapiro-Wilk检验表明,耳石长(Z=0.977,P=0.496)和耳石宽(Z=0.987,P=0.882)均符合正态分布。

阿拉斯加狭鳕耳石外型较大,耳石长度范围为12.93~17.05 mm(15.46±1.03),耳石宽度范围为5.70~7.78 mm(6.52±0.46),即耳石长明显大于耳石宽(配对样本t检验,t=71.710,P<0.001),其整体上趋于近椭圆形,耳石狭长,表面光滑,中间厚,两端较薄,耳石呈拱形突起,内侧凹面,外侧凸面。背部和腹部边缘有许多不规则的锯齿状突起,背部突起较大且稀疏,而腹部突起较小且密集。耳石表面辐射状条纹清晰,具有明显的基叶和翼叶,但主间沟不明显,耳石整体结构凹凸不平。阿拉斯加狭鳕耳石各形态参数统计值如表1所示。

表1 阿拉斯加狭鳕耳石形态参数值

Tab.1 Morphological parameters of otolith in Alaska pollock Gadus chalcogrammus

形态参数 morphological parameter最小值minimal value最大值maximal value平均值±标准差mean±S.D.耳石长OL/mm12.9317.0515.46±1.03耳石宽OW/mm5.707.786.52±0.46背长DL/mm11.7316.0114.23±1.13背宽DW/mm2.913.893.40±0.25基叶长RL/mm6.509.958.36±0.77翼叶长ARL/mm4.948.476.81±0.91耳石周长P/mm36.5949.1843.01±3.54耳石面积S/mm255.8093.8873.20±8.93耳石质量WW/g0.15070.29480.2003±0.0327

2.2 主成分分析

试验中对阿拉斯加狭鳕耳石的9项形态特征参数进行主成分分析(表2),本研究中共选取了两个主成分,主成分1和主成分2的贡献率依次为66.08%和14.41%,累计贡献率达80.49%,这2个主成分的累计贡献率较高,足以代表整体指标,因此,可认为其保留了阿拉斯加狭鳕耳石形态的主要信息,故可用这两个形态因子描述阿拉斯加狭鳕耳石的形态特征。由表2可得出第1主成分与耳石面积、背长、耳石长存在密切的相关性,载荷系数均在0.9以上,故可以用上述参数代表9项形态指标来描述阿拉斯加狭鳕耳石的形态特征。

表2 阿拉加狭鳕9个形态参数中两个主要成分负荷值和贡献率

Tab.2 Loading and contribution of two principal components in nine morphological parameters of Alaska pollock Gadus chalcogrammus based on principal component analysis

形态参数morphological parameter主成分principal component12耳石面积S0.976-0.117背长DL0.9370.073耳石长OL0.9200.157耳石周长P0.880-0.095耳石宽OW0.877-0.389耳石质量WW0.811-0.058背宽DW0.664-0.510基叶长RL0.5760.691翼叶长ARL0.5460.588贡献率contribution rate/%66.0814.41

2.3 矢耳石的生长

2.3.1 耳石相对尺寸 耳石的相对尺寸(即各形态参数值与体长比)显示,阿拉斯加狭鳕耳石长、耳石宽、背宽、背长、基叶长、翼叶长和耳石周长与体长之比随着鱼体长的增加呈现下降趋势,这7个形态参数与体长的比值分别由体长为35 cm时的4.61%、1.80%、0.94%、4.36%、2.44%、2.18%和12.31%,下降到体长为43 cm时的4.05%、1.65%、0.81%、3.78%、1.91%、1.63%和11.21%。阿拉斯加狭鳕耳石生长速度比鱼体生长速度慢。

2.3.2 耳石主要形态特征与体长的关系 根据AIC值最小原则(表3),阿拉斯加狭鳕耳石长、耳石宽、背长、基叶长、翼叶长、面积和周长与鱼体体长(L)间以指数函数拟合为最佳,而背宽与鱼体体长以幂函数拟合为最佳(图3)。利用回归分析检验上述拟合的函数关系(P<0.001),说明耳石各形态指标与鱼体体长拟合的函数关系成立。

表3 耳石主要形态参数与体长关系拟合方程的AIC值比较

Tab.3 Comparison of main morphological parameters of otolith with AIC values of different growth equations for body length

形态参数morphological parameter 线性方程linear对数方程logarithmic幂函数方程power指数方程exponential耳石长OL-4.19-3.83-249.21-249.49耳石宽OW-89.23-192.41-257.08-257.41背宽DW-131.27-131.33-240.64-240.41背长DL4.555.00-232.85-233.24基叶长RL-27.13-27.16-215.84-215.85翼叶长ARL-7.67-7.61-180.07-180.13耳石周长P107.37107.54-229.50-229.68耳石面积S179.78180.41-204.55-205.07

2.3.3 耳石主要形态特征与湿体质量的关系 根据AIC最小原则(表4),阿拉斯加狭鳕耳石长、背宽、背长、基叶长、周长和面积与湿体质量(W)间以幂函数拟合为最佳,而耳石宽、翼叶长与湿体质量间以指数函数拟合为最佳(图4)。利用回归分析检验上述拟合的函数关系(P<0.001),说明耳石各形态指标与鱼体湿体质量(W)拟合的函数关系成立。

2.3.4 耳石各形态参数值与耳石长之比与体长的关系 阿拉斯加狭鳕耳石宽、背宽、背长、基叶长、翼叶长与耳石长之比的平均值基本保持不变,分别维持在42.01%(95%置信区间为41.31%~42.71%)、21.87%(95%置信区间为21.37%~22.36%)、91.88%(95%置信区间为90.95%~92.81%)、54.13%(95%置信区间为52.96%~55.29%)、44.16%(95%置信区间为42.59%~45.74%)。随着鱼体体长的增加,耳石各区的总体形态基本维持不变(图5),即耳石外部形态趋于较稳定状态。

图3 阿拉斯加狭鳕耳石各形态参数与体长的关系
Fig.3 Relationship between morphological parameters of otolith and body length of Alaska pollock Gadus chalcogrammus

表4 耳石主要形态参数与湿体质量关系拟合方程的AIC值比较

Tab.4 Comparison of main morphological parameters of otolith with AIC values of different growth equations for wet body weight

形态参数morphological parameter线性方程linear对数方程logarithmic幂函数方程power指数方程exponential耳石长OL-3.59-4.26-249.77-249.21耳石宽OW-92.65-92.29-260.16-260.52背宽DW-132.21-132.62-241.80-241.33背长DL2.061.91-235.64-235.44基叶长RL-27.19-27.80-216.25-215.71翼叶长ARL-7.83-7.76-180.13-180.19耳石周长P108.75108.09-229.32-228.62耳石面积S177.56177.44-207.31-206.88

3 讨论

3.1 耳石外型特征

阿拉斯加狭鳕体长范围为35~43 cm时,矢耳石相对较大,耳石平均长度和宽度分别为15.47、6.52 mm,耳石背部和腹部具有许多锯齿状的波纹状突起,耳石表面辐射状条纹清晰。耳石具有明显的基叶和翼叶,且基叶长大于翼叶长,主间沟不明显。耳石两端向上翘,前端稍尖,后端闭合呈钝圆形,整体结构凹凸不平,呈拱形状突起。阿拉斯加狭鳕耳石外型与Cardinale等[21]研究的同属中的大西洋鳕鱼Gadus morhua耳石外型存在较大差异,这说明同种、属鱼类耳石形态可能也会差异较大[15-16],但其相同点是均具有明显的基叶和翼叶,主间沟不明显。由主成分分析可得,阿拉斯加狭鳕的前2个主成分方差累计贡献率可达80.49%,耳石面积、背长、耳石长为阿拉斯加狭鳕耳石最具代表性的形态指标,因此,这些指标可代表阿拉斯加狭鳕其他形态指标来描述矢耳石外型特征,故在今后开展类似鱼类耳石形态的研究中可侧重于研究这几项形态指标。

图4 阿拉斯加狭鳕耳石各形态参数与湿体质量的关系
Fig.4 Relationship between morphological parameters of otolith and wet body weight of Alaska pollock Gadus chalcogrammus

3.2 耳石外型与鱼体之间的关系

不同鱼类的耳石与鱼体大小呈现不同的关系[22-23],因此,可用不同的生长模型拟合耳石各形态指标与鱼体体长和体质量的关系。本研究中,阿拉斯加狭鳕耳石形态指标与鱼体体长和湿体质量呈显著的指数函数和幂函数关系,这说明阿拉斯加狭鳕耳石生长与个体生长密切相关。Longenecker[24]曾利用耳石大小推算个体或群体的生长发育时期和生长规律特性,在今后阿拉斯加狭鳕的研究中,可利用本研究中建立的函数关系式,通过测量耳石形态指标推算阿拉斯加狭鳕的体长和湿体质量。

图5 耳石各参数值长度比与阿拉斯加狭鳕体长的关系
Fig.5 Relationship between ratio of different parameters to otolith length and body length of Alaska pollock Gadus chalcogrammus

随着阿拉斯加狭鳕体长的增加,耳石各形态指标也随之增加,但其相对尺寸(即耳石各形态参数与体长之比)却逐渐减少,这表明在阿拉斯加狭鳕生长发育过程中,体长在增加,耳石整体形态却不易发生变化,同时也说明在生长过程中鱼体生长要大于耳石生长。由耳石各形态指标与耳石长度之比也可看出,耳石宽、背宽、背长、基叶长和翼叶长与耳石长度之比基本维持稳定状态,分别为42.01%、21.87%、91.88%、54.13%和44.16%。耳石的长宽比例并不随体长及湿体质量的增加而发生改变,耳石的形状与鱼类体长和湿体质量的增长不存在显著的相关关系。Campana等[6]研究表明,成鱼的耳石外型特征一般比较稳定,但在幼鱼阶段的生长发育过程中,耳石的形态会发生较大变化,可推测本试验样本可能是阿拉斯加狭鳕成鱼样本。由于耳石具有物种特异性和种间差异性,因此,其耳石形态多用于种群鉴定和近缘种的鉴别[25],如郭弘艺等[26]通过对中国鲚属鱼类的矢耳石形态特征进行研究,鉴别了七丝鲚、凤鲚、刀鲚和短颌鲚这4种鲚属鱼类,但本研究中阿拉斯加狭鳕矢耳石形态特征不存在差异性,耳石外型基本相同,因此,可推测本试验样本为同一种属的阿拉斯加狭鳕。

由于本试验样本为商业捕捞获取,故采样时间和海域受到限制,采集样本数量较少,对研究的结果可能会有一定的影响,同时缺少早期幼鱼耳石的样本,故无法比较阿拉斯加狭鳕幼体和成体耳石外型差异,以及探讨幼体和成体鱼耳石大小与鱼体长和湿体质量之间关系的差异,这需要在以后的研究中加以补充和完善。

致谢:感谢辽宁远洋渔业有限公司提供阿拉斯加狭鳕样本!

参考文献

[1] 高天翔,张肖荣,王丹,等.几种鳕鱼的生物学初步研究[J].海洋湖沼通报,2003(1):35-42.

[2] Ianelli J N,Barbeaux S,Honkalehto T,et al.Assessment of Alaska pollock stock in the Eastern Bering Sea[R].Seattle:Alaska Fisheries Science Center,National Marine Fisheries Service,2005:32-124.

[3] 孙枫.东北太平洋狭鳕资源的开发利用[J].国外水产,1990(1):40-42

[4] 楼乔明,张问,刘连亮,等.狭鳕鱼皮脂肪酸组成分析及其营养评价[J].核农学报,2016,30(2):332-337.

[5] 张治国,王卫民.鱼类耳石研究综述[J].湛江海洋大学学报,2001,21(4):77-83.

[6] Campana S E,Casselman J M.Stock discrimination using otolith shape analysis[J].Canadian Journal of Fisheries and Aquatic Sciences,1993,50(5):1062-1083.

[7] 廖锐,区又君.鱼类耳石研究和应用进展[J].南方水产,2008,4(1):69-75.

[8] Arkhipkin A I.Towards identification of the ecological lifestyle in nektonic squid using statolith morphometry[J].Journal of Molluscan Studies,2003,69(3):171-178.

[9] Panella G.Fish otoliths:daily growth layers and periodical patterns[J].Science,1971,173(4002):1124-1127.

[10] Messieh S N.Use of otoliths in identifying herring stocks in the southern Gulf of St.Lawrence and adjacent waters[J].Journal of the Fisheries Research Board of Canada,1972,29(8):1113-1118.

[11] 林景祺.狭鳕等三种鳕鱼生态和资源[J].海洋科学,1994,18(2):25-29.

[12] Blood D M,Matarese A C,Yoklavich M M.Embryonic development of walleye pollock,Theragra chalcogramma,from Shelik of Strait,Gulf of Alaska[J].Fishery Bulletin,1994,92(2):207-222.

[13] Jung K M,Kang S,Kim S,et al.Ecological characteristics of walleye pollock eggs and larvae in the southeastern Bering Sea during the late 1970s[J].Journal of Oceanography,2006,62(6):859-871.

[14] Yoklavich M M,Bailey K M.Hatching period,growth and survival of young Walleye pollock Theragra chalcogrammaas determined from otolith analysis[J].Marine Ecology Progress Series,1990,64:13-23.

[15] 李凤莹,朱国平,许柳雄.西南大西洋拉式南美南极鱼耳石外型特征[J].生态科学,2014,33(1):7-11.

[16] 魏联,刘子俊,朱国平.阿根廷外海拉式南美南极鱼矢耳石形态特征分析[J].海洋渔业,2017,39(3):256-266.

[17] 李史民,朱国平.南极南设得兰群岛南极电灯鱼耳石形态学研究[J].大连海洋大学学报,2015,30(6):674-680.

[18] Gaemers P A M.Taxonomic position of the Cichlidae (Pisces,Perciformes) as demonstrated by the morphology of their otoliths[J].Netherlands Journal of Zoology,1983,34(4):566-595.

[19] Akaike H.Information theory and an extension of the maximumlikelihood principle[M]//Petrov B N,Csaki F.Second International Symposium on Information Theory.Budapest:Akademiai Kiado,1973:267-281.

[20] Burnhan K P,Anderson D R.Model selection and multi-model inference:a practical information-theoretic approach (2nd Edition)[M].New York:Springer,2003.

[21] Cardinale M,Doering-Arjes P,Kastowsky M,et al.Effects of sex,stock,and environment on the shape of known-age Atlantic cod (Gadus morhua) otoliths[J].Canadian Journal of Fisheries and Aquatic Sciences,2004,61(2):158-167.

[22] 韩青鹏,陆化杰,金岳,等.南海北部海域中国枪乌贼耳石的形态学分析[J].广东海洋大学学报,2017,37(1):1-6.

[23] 王晓杰,朱清澄,张阳,等.西北太平洋秋刀鱼耳石形态学分析[J].江苏农业科学,2013,41(4):210-213.

[24] Longenecker K.Relationships between otolith size and body size for Hawaiian reef fishes[J].Pacific Science,2008,62(4):533-540.

[25] Campana S E.Otolith science entering the 21 stcentury[J].Marine and Freshwater Research,2005,56(5):485-495.

[26] 郭弘艺,唐文乔,魏凯,等.中国鲚属鱼类的矢耳石形态特征[J].动物学杂志,2007,42(1):39-47.

Morphometric features of sagittal otolith for Alaska pollock Gadus chalcogrammus in the western Bering Sea

DUAN Mi1,2, WEI Lian1,2, ZHU Guo-ping1,2,3

(1.College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China; 2.Polar Marine Ecosystem Laboratory, The Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China;3.National Engineering Research Center for Oceanic Fisheries, Shanghai 201306, China)

Abstract The morphologic features of sagittal otolith were investigated in 53 Alaska pollock Gadus chalcogrammus samples collected from western Bering Sea to evaluate the relationship between otolith and growth of Alaska pollock. The Alaska pollock has large and nearly oval sagittal otolith, with obvious rostrum and antirostrum, without obvious excisura, and Alaska pollock with body length of 35-43 cm had otolith with length (OL) of 12.93-17.05 mm and width (OW)of 5.70-7.78 mm. There is no significant difference in otolith morphology between left and right otolith(P>0.05). The statolith was increased gradually in size with the increase in body length (L), but the ratio of each morphologic parameter to body length was decreased slowly, with almost constant proportion of morphological structure. The principal component analysis(PCA) of the nine indices describing the otolith morphology revealed that otolith acreage(A), dorsal length(DL) and otolith length(OL)was used as typical morphological characteristics of Alaska pollock otolith. The relationship between the main morphological indices of the otolith and body length and wet body weight(W) indicated that there was an exponential function between OL, OW, DL, rostrum length(RL), antirostrum length(ARL), and tolith area (S) and circumference (P) and wet body weight. There was a power function relationship between OL, and DW, and wet body weight. The findings provide basic data for age determination and population identification.

Key words Gadus chalcogrammus; otolith; morphologic feature; principal component analysis

DOI10.16535/j.cnki.dlhyxb.2018.04.013

中图分类号S912

文献标志码:A

文章编号:2095-1388(2018)04-0492-07

收稿日期 2017-11-03

基金项目国家自然科学基金资助项目(41776185);教育部留学回国人员科研启动基金资助项目

作者简介段咪(1995—),女,硕士研究生。E-mail:duanmi_1995@163.com

通信作者朱国平(1976—),男,博士,副教授。E-mail:gpzhu@shou.edu.cn