Heterogeneity of mantle length-body mass relationship in different Ommastrephes bartramii populations based on linear mixed model
-
摘要: 为研究不同群体柔鱼 (Ommastrephes bartramii) 生长的异质性,文章根据2016—2018年在北太平洋采集的柔鱼样本,对该物种胴长和体质量的性别和月间差异进行了分析。构建了广义线性模型 (Generalized linear model, GLM) 和9个线性混合模型 (Linear mixed model, LMM),用于研究不同群体柔鱼的胴长和体质量关系 (m=aLb) 在性别及时间上的差异。结果显示,东、西部群体不同性别柔鱼胴长、体质量均存在极显著差异 (P<0.01);各月份间也均存在极显著差异 (P<0.01)。东部群体异速生长参数 (b) 的固定值 (2.763 4) 小于3,表示其为负异速增长,肥满度随胴长增加而减少,胴体趋于细长;西部群体b的固定值 (3.117 19) 大于3,表示其为正异速增长,肥满度随胴长增加而增大,胴体趋于宽厚。研究表明,性别和月份对不同群体柔鱼胴长和体质量关系均有显著影响,线性混合模型能将性别和月份的异质性通过随机效应在单个模型中更准确、方便地体现出来,进一步证实了此模型在个体生长异质性研究中的优势。Abstract: We had collected neon flying squid (Ommastrephes bartramii) samples from 2016 to 2018 in North Pacific, and studied their gender and monthly differences of mantle length-body mass relationship, so as to understand the heterogeneity in the growth of different O. bartramii populations. One generalized linear (GLM) model and nine linear mixed models (LMM) were applied to study the difference in gender and time for mantle length-body mass relationship (m=aLb). The results show that there was a significant difference in the mantle length and body mass between female and male O. bartramii of eastern and western populations (P<0.01). Moreover, there was a significant difference between different months (P<0.01). The fixed value of b (2.763 4) was less than 3 in the eastern population, which indicates that O. bartramii was growing at negative allometry, and the degree of fullness was negatively related to mantle length, and the body tended to be slender. The fixed value of b (3.117 19) was more than 3 in the western population, which indicates that O. bartramii was growing at positive allometry, and the degree of fullness was positively related to the mantle length, and the body tended to be thick. The study indicates that month and gender have a significant effect on the mantle length-body mass relationship for different O. bartramii populations. The LMM reflects the heterogeneity of gender and month more accurately and conveniently in a single model through random effects, thereby further confirming the advantage of this model in the study of heterogeneity in individual growth.
-
图 4 柔鱼不同月份胴长体质量分布图
a. 东部群体胴长分布图;b. 东部群体体质量分布图;c. 西部群体胴长分布图;d. 西部群体体质量分布图
Figure 4. Distribution of mantle length and body mass of O. bartramii in different months
a. Distribution of mantle length for eastern population; b. Distribution of body mass for eastern population; c. Distribution of mantle length for western population t; d. Distribution of body mass for western population
图 5 不同群体柔鱼体质量-胴长关系
a. 东部群体雌性体质量-胴长关系;b. 东部群体雄性体质量-胴长关系;c. 西部群体雌性体质量-胴长关系;d. 西部群体雄性性体质量-胴长关系
Figure 5. Mantle length-body mass relationship for different O. bartramii populations
a. Relationship between BM and ML of female for eastern population; b. Relationship between BM and ML of male for eastern population; c. Relationship between BM and ML of female for western population; d. Relationship between BM and ML of male for western population
表 1 不同月份柔鱼样本基本信息
Table 1. Basic information of O. bartramii samples in different months
群体
Population年份
Year月份
Month样本数量
Number of sample胴长范围
Range of mantle length/mm东部群体 Eastern population 2016—2017 5 394 163~500 6 456 212~516 西部群体 Western population 2016—2018 7 124 187~290 8 283 195~357 9 331 199~386 10 165 190~377 11 105 215~421 表 2 不同群体柔鱼胴长-体质量关系模型的拟合效果指标
Table 2. Fitting effect index (AIC) of mantle length-body mass relationship models for different O. bartramii populations
种群
Population模型缩写
Model abbreviation模型
Model对数转化
Log-transformed赤池信息量准则
AIC均方根误差
RMSE东部群体
Eastern populationLMR m = a×Lb ln(m)=ln(a)+b×ln(L) −466.402 0.183 265 RSS m=a×L(b+RSS) ln(m)=ln(a)+(b+RSS)×ln(L) −472.338 0.181 617 RSI m=a×exp(RSI)×Lb ln(m)=[ln(a)+RSI]+b×ln(L) −472.27 0.181 641 RSI&S m=a×exp(RSI)×L(b+RSS) ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L) −469.057 0.181 428 RMS m=a×L(b+RMS) ln(m)=ln(a)+(b+RMS)×ln(L) −489.107 0.179 784 RMI m=a×exp(RMI)×Lb ln(m)=[ln(a)+RMI]+b×ln(L) −492.772 0.179367 RMI&S m=a×exp(RMI)×L(b+RSI) ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L) −595.137 0.166 468 RS&MS m=a×L(b+RSS+RMS) ln(m)=ln(a)+(b+RSS + RMS)×ln(L) −494.452 0.178 558 RS&MI m=a×exp(RSI)×exp(RMI)×Lb ln(m)=[ln(a)+RSI + RMI]+b×ln(L) −497.694 0.178 211 RS&MI&S m=a×exp(RSI)×exp(RMS)×L(b+RSI+RSI) ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L) −596.138 0.165 394 西部群体
Western populationLWR m = a×Lb ln(m)=ln(a)+b×ln(L) −1 795.63 0.099 003 RSS m=a×L(b+RSS) ln(m)=ln(a)+(b+RSS)×ln(L) −1 797.15 0.098 614 RSI m=a×exp(RSI)×Lb ln(m)=[ln(a)+RSI]+b×ln(L) −1 797.31 0.098 603 RSI&S m=a×exp(RSI)×L(b+RSS) ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L) −1 796.72 0.098 512 RMS m=a×L(b+RMS) ln(m)=ln(a)+(b+RMS)×ln(L) −1 807.04 0.097 626 RMI m=a×exp(RMI)×Lb ln(m)=[ln(a)+RMI]+b×ln(L) −1 807.41 0.097 615 RMI&S m=a×exp(RMI)×L(b+RSI) ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L) −1 807.59 0.097 215 RS&MS m=a×L(b+RSS+RMS) ln(m)=ln(a)+(b+RSS + RMS)×ln(L) −1 808.14 0.097 347 RS&MI m=a×exp(RSI)×exp(RMI)×Lb ln(m)=[ln(a)+RSI + RMI]+b×ln(L) −1 808.56 0.097 333 RS&MI&S m=a×exp(RSI)×exp(RMS)×L(b+RSI+RSI) ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L) −1 808.09 0.096 976 注:第一列为柔鱼的不同群体,第二列为第三、四列模型的缩写 RSS、RSI、RSI&S、RMS、RMI、RMI&S、RS&MS、RS&MI 和 RS&MI&S 分别为斜率 b 的性别随机效应、截距 ln(a) 的性别随机效应、截距 ln(a) 与斜率 b的性别随机效应、斜率 b 的月份随机效应、截距 ln(a) 的月份随机效应、截距 ln(a) 与斜率 b 的月份随机效应、斜率 b的性别与月份随机效应、截距 ln(a) 的性别与月份随机效应,以及截距 ln(a) 与斜率 b 的性别与月份随机效应Note: The first column shows different O. bartramii populations. The second column shows the abbreviations of models detailed in the third and fourth columns. RSS, RSI, RSI&S, RMS, RMI, RMI&S, RS&MS, RS&MI和RS&MI&S indicate random effects on slope (b) in sexual difference, random effects on intercept [ln (a)] in sexual difference, random effects on intercept [ln (a)] and slope (b) in sexual difference, random effects on slope (b) in months, random effects on intercept [ln (a)] in months, random effects on intercept [ln (a)] and slope (b) in months, random effects on slope (b) in sexual difference and months, random effects on intercept [ln (a)] fin sexual difference and months, and random effects on intercept [ln (a)] and slope (b) in sexual difference and months, respectively. 表 3 不同月份东西部群体柔鱼胴长与体质量关系
Table 3. Mantle length-body mass relationship for different O. bartramii populations in different months
群体
Population月份
Month胴长与体质量关系
Mantle length-mass relationship样本数量
Number of samples判定系数R2
Coefficient of determinationP 东部群体
Eastern population5 m =1.9×10−3L2.295 7 394 0.670 6 0.001 6 m =1×10−5L3.160 1 456 0.973 8 0.001 西部群体
Western population7 m = 9×10−6L3.205 3 124 0.889 5 0.001 8 m = 2×10−5L3.090 9 283 0.934 5 0.001 9 m = 2×10−5L3.059 7 331 0.952 0 0.001 10 m = 1×10−5L3.178 3 165 0.972 8 0.001 11 m = 5×10−5L2.902 8 105 0.947 7 0.001 表 4 不同性别东西部群体柔鱼胴长与体质量关系
Table 4. Mantle length-body mass relationship of O. bartramii of different sex for different populations
群体
Population性别
Sex胴长与体质量关系
Mantle length-mass relationship样本数量
Number of samples判定系数R2
Coefficient of determinationP 东部群体
Eastern population雌性 m=3×10−5L2.981 8 744 0.892 9 0.001 雄性 m=2×10−4L2.621 0 106 0.794 3 0.001 西部群体
Western population雌性 m=1×10−5L3.159 9 512 0.968 1 0.001 雄性 m =2×10−5L3.103 5 496 0.941 5 0.001 -
[1] 王尧耕, 陈新军. 世界大洋性经济柔鱼类资源及其渔业[M]. 北京: 海洋出版社, 2005: 124-137. [2] MURATA M. Oceanic resources of squids[J]. Mar Freshw Behav Phy, 1990, 18(1): 19-71. doi: 10.1080/10236249009378779 [3] BOWER J R, ICHII T. The red flying squid (Ommastrephes bartramii): a review of recent research and the fishery in Japan[J]. Fish Res, 2005, 76(1): 39-55. doi: 10.1016/j.fishres.2005.05.009 [4] 陈新军. 渔业资源与渔场学[M]. 北京: 海洋出版社, 2004: 29-212. [5] YATSU A, TANAKA H, MORI J. Population structure of the neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Jpn Mar Fish Resour Res Cent, 1998: 31-48. [6] FANG Z, LIU B, LI J, et al. Stock identification of neon flying squid (Ommastrephes bartramii) in the North Pacific Ocean on the basis of beak and statolith morphology[J]. Sci Mar, 2014, 78(2): 239-248. doi: 10.3989/scimar.03991.06A [7] 方舟, 陈新军, 瞿俊跃, 等. 北太平洋柔鱼角质颚形态及生长年间差异[J]. 上海海洋大学学报, 2020, 29(1): 109-120. doi: 10.12024/jsou.20181202485 [8] WATANABE H, KUBODERA T, ICHII T, et al. Feeding habits of neon flying squid Ommastrephes bartramii in the transitional region of the central North Pacific[J]. Mar Ecol Prog Ser, 2004, 266(1): 173-184. [9] CHEN C S, CHIU T S. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific[J]. Fish Res, 2003, 63(3): 349-366. doi: 10.1016/S0165-7836(03)00101-2 [10] YATSU A. Distribution of flying squid, Ommastrephes bartramii, in the north Pacific based on a jigging survey, 1976−1983[J]. Bull Natl Res Inst Far Seas Fish, 1992(29): 13-37. [11] YATSU A, MORI J. Early growth of the autumn cohort of neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Fish Res, 2000, 45(2): 189-194. doi: 10.1016/S0165-7836(99)00112-5 [12] SUN M C, ZHANG S H, QIAN W G. Approach on biological characteristics of Ommastrephes bartramii on the central and eastern squid fishing grounds in the North Pacific[J]. Mar Tieher, 2001, 23(1): 21-24. [13] YATSU A, MIDORIKAWA S, SHIMADA T, et al. Age and growth of the neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Fish Res, 1997, 29(3): 257-270. doi: 10.1016/S0165-7836(96)00541-3 [14] MCINTYRE A D. Cephalopod life cycles. Vol. 1. Species accounts[J]. J Anim Ecol, 1986, 55(1): 395-396. [15] GRANADOS-AMORES J, SALINAS-ZAVALA C A, FLORES-ORTEGA J R, et al. Length-weight relationship and condition factor for 7 loliginid squid species in Mexican waters[J]. Cienc Mar, 2019, 45(4): 175-180. doi: 10.7773/cm.v45i4.3015 [16] VASCONCELOS P, PEREIRA F, CARVALHO A N, et al. Weight-length relationships and relative growth of the cuttlefish (Sepia officinalis): causes and effects of hypoallometry[J]. Thalassas: Int J Mar Sci, 2018, 32(4): 323-331. [17] 衷思剑, 麻秋云, 刘淑德, 等. 基于线性混合效应模型的黄鮟体长体重关系的时空差异[J]. 中国水产科学, 2018, 25(6): 152-160. [18] XU H, SUN Y J, WANG X J, et al. Linear mixed-effects models to describe individual tree crown width for China-fir in Fujian province, southeast China[J]. PLoS One, 2015, 10(4): e0122257. doi: 10.1371/journal.pone.0122257 [19] MA Q Y, JIAO Y, REN Y P. Linear mixed-effects models to describe length-weight relationships for yellow croaker (Larimichthy spolyactis) along the north coast of China[J]. PLoS One, 2017, 12(2): e0171811. doi: 10.1371/journal.pone.0171811 [20] 瞿俊跃, 方舟, 陈新军. 基于线性混合模型的蓝点马鲛叉长体重生长关系的月间及性别差异[J]. 中国水产科学, 2020, 8: 953-960. [21] 廖为耕. 浅谈北太平洋之海洋环境及赤鱿之洄游[J]. 中国水产(台刊), 1980, 457: 39-47. [22] LIPINSKI M R, UNDERHILL L G. Sexual maturation in squid: quantum or continuum[J]. S Afr J Mar Sci, 1995, 15: 207-223. doi: 10.2989/02577619509504844 [23] 李春喜, 邵云, 姜丽娜. 生物统计学[M]. 北京: 科学出版社, 2008: 85-110. [24] GUEORGUIEVA R, KRYSTAL J H. Move over anova: progress in analyzing repeated-measures data and its reflection in papers published in the archives of general psychiatry[J]. Arch Gen Psychiat, 2004, 61(3): 310-317. doi: 10.1001/archpsyc.61.3.310 [25] 王斌会. 多元统计分析及R语言建模[M]. 广州: 暨南大学出版社, 2010: 157-158. [26] KEYS A B. The weight-length relationship in fishes[J]. P Natl A Sci USA, 1928, 14(12): 922-925. doi: 10.1073/pnas.14.12.922 [27] BURNHAM K P, ANDERSON D R. Model selection and multi-model inference: a practical information theoretic approach[M]. 2nd ed. New York: Springer Press, 2002: 484-488. [28] BERNAL-RUSIEL J L, GREVE D N, REUTER M, et al. Statistical analysis of longitudinal neuroimage data with Linear Mixed Effects models[J]. NeuroImage, 2013, 66(2): 249-260. [29] 黄洪亮, 郑元甲, 程家骅. 北太平洋海区柔鱼生物学特征研究[J]. 海洋渔业, 2003, 25(3): 21-24. [30] ICHII T, MAHAPATRA K, OKAMURA H, et al. Stock assessment of the autumn cohort of neon flying squid (Ommastrephes bartramii) in the North Pacific based on past large-scale high seas driftnet fishery data[J]. Fish Res, 2006, 78(2): 286-297. [31] KATUGIN O N. Patterns of genetic variability and population structure in the North Pacific squids Ommastrephes bartramii, Todarodes pacificus, and Berryteuthis magister[J]. Bull Mar Sci, 2002, 71(1): 383-420. [32] HILE R. Age and growth of the cisco, Leucichthys artedi (le Sueur), in the lakes of the Northeastern highlands, Wisconsin[J]. Bull Bur Fish, 1936, 48(11): 211-317. [33] KEITH A, BIGELOW, KEVIN C, et al. Hatch dates and growth of Ommastrephes bartramii paralarvae from Hawaiian waters as determined from statolith analysis[J]. Rec Adv Fish Biol, 1993: 15-24. [34] 唐峰华, 伍玉梅, 樊伟. 2010年北太平洋海区巴特柔鱼(Ommastrephe bartrami)的生物学特征[J]. 海洋科学进展, 2012, 30(2): 257-265. doi: 10.3969/j.issn.1671-6647.2012.02.013 [35] FORSYTHE J W. Accounting for the effect of temperature on squid growth in nature: from hypothesis to practice[J]. Mar Freshw Res, 2004, 55(4): 331-339. doi: 10.1071/MF03146 [36] ICHII T, MAHAPATRA K, SAKAI M, et al. Life history of the neon flying squid: effect of the oceanographic regime in the North Pacific Ocean[J]. Mar Ecol Prog Ser, 2009, 378: 1-11. doi: 10.3354/meps07873 [37] 章寒, 郑基, 虞聪达, 等. 东太平洋赤道海域茎柔鱼主要生物学特性比较研究[J]. 中国水产科学, 2019, 26(4): 745-755. [38] DOMINGUES P M, SYKES A, ANDRADE J P. The effects of temperature in the life cycle of two consecutive generations of the cuttlefish Sepia officinalis (Linnaeus, 1758), cultured in the Algarve (South Portugal)[J]. Aquacult Int, 2002, 10(3): 207-220. doi: 10.1023/A:1022148802078 [39] DOMINGUES P M, POIRIER R, DICKEL L, et al. Effects of culture density and live prey on growth and survival of juvenile cuttlefish, Sepia officinalis[J]. Aquacult Int, 2003, 11(3): 225-242. doi: 10.1023/A:1024803802486 [40] 陈锦辉, 王学昉, 田思泉, 等. 基于混合效应模型的印度洋大眼金枪鱼生长特征异质性分析[J]. 中国水产科学, 2020, 27(5): 570-578. -