南麂列岛中华蛸(Octopus sinensis)形态与遗传多样性分析*

李嘉华 陈 舜 陈万东 谢尚微 倪孝品 郑小东①

南麂列岛中华蛸()形态与遗传多样性分析*

李嘉华1, 2陈 舜3陈万东3谢尚微3倪孝品3郑小东1, 2①

(1. 中国海洋大学海洋生物多样性与进化研究所 山东青岛 266003; 2. 中国海洋大学海水养殖教育部重点实验室 山东青岛 266003; 3. 南麂列岛国家海洋自然保护区管理局 浙江温州 325400)

2020年6～12月在南麂列岛采捕一种中大型章鱼44只, 描述了形态特征, 采用多元分析方法分析了形态多样性, 并比较了与真蛸的差异, 利用线粒体细胞色素氧化酶亚基I ()基因序列构建了系统发育树, 并计算了Kimura-2-Parameter (K2P)遗传距离。结果表明, 此种章鱼形态特征与中华蛸基本一致, 而与真蛸存在显著的形态差异: 其第二、三对腕的扩大吸盘位于第12～15个吸盘之间, 而真蛸则位于第15～19个吸盘之间; 茎化腕吸盘数少于真蛸; 舌叶指数小于真蛸。主成分分析与判别分析能将真蛸群体显著区分开, 判别准确率为100%, 而其余四个群体(南麂列岛与三个中华蛸群体)间存在重叠; 聚类分析表明南麂列岛群体与中华蛸更近, 二者均与真蛸存在较大距离。在遗传多样性分析中, 南麂列岛群体单倍型4个, 单倍型多样性水平为0.320±0.121, 多态位点6个, 核苷酸多样性指数为0.001 11; 系统发育树表明, 该群体与中华蛸亲缘关系最近, K2P遗传距离0.14%, 而与真蛸复合体其他类型为2.96%～12.11%。因此, 从形态和分子水平鉴定南麂列岛章鱼为中华蛸。

中华蛸; 南麂列岛; 形态多样性; 多元分析;基因

中华蛸(d’Orbigny, 1841)隶属软体动物门(Mollusca)、头足纲(Cephalopoda)、八腕目(Octopoda)、蛸科(Octopodidae)、蛸属(), 广泛分布于日本、韩国以及中国沿海的岩礁、砂底, 是重要的经济蛸类。1834年, d’Orbigny整理了现存头足类图谱, 根据当时日本百科全书《Wakansansaizue》添加了中华蛸的相关插图, 于1841年描述并命名新种—中华蛸(d’Orbigny, 1835-1848)。Sasaki(1929)指出中华蛸为真蛸()的同物异名, 此后, 日本采集的样品均被称为“真蛸”。21世纪初, 日本学者仍沿用种名真蛸(Sakaguchi, 2000, 2005), Warnke等(2004)通过线粒体基因与基因联合分析, 也认为在日本和中国台湾沿海采集的野生样本与真蛸为同一物种。然而, 由于存在地理隔离, Norman(2000)认为日本分布的所谓“真蛸”与真蛸模式种很可能不是同一物种。Norman等(2005)在蛸科分类修订中将“中华蛸”暂时归类于不确定种名。近年来, 基于形态学与分子标记的研究结果表明, 真蛸存在多个隐存种, 是一个大的复合物种(complex species), 如、、cf.、(Amor, 2017a, 2017b; Van Nieuwenhove, 2019; Avendaño, 2020)。Gleadall(2016)对日本九州与本州沿海采集的“真蛸”进行了重描述, 发现与真蛸存在显著形态差异, 重新确定d’Orbigny, 1841为有效种名。Amor等(2017b)采用形态学与分子标记法得到的结果也支持上述结论。

国内已发表的论文专著多沿用物种名“真蛸”, 指出该物种为我国东南沿海重要经济蛸类, 其肉嫩味美, 营养价值高, 是妇女生乳的滋补品(董正之, 1988), 有关其亲体培育、胚胎发育、幼体生长、生理生态以及基因组学、转录组学等方面研究已有报道(蔡厚才等, 2007, 2009; 刘兆胜等, 2011; 郑小东等, 2011; 孙田田等, 2012; 冯雪等, 2013; 肖懿哲等, 2019; Li, 2020)。然而, 该物种个体大, 不易保存, 我国在其分类、遗传多样性方面鲜有文献记载。本文采用形态参数与分子标记相结合的分析方法研究南麂列岛重要经济贝类“真蛸”(叶鹏等, 2006)以确定其正确种名, 旨在为南麂列岛种质资源保护、合理开发利用以及丰富保护区物种基因库提供理论依据。

1 材料与方法

1.1 实验材料

2020年6～12月, 从浙江省南麂列岛海域(大檑岛、马祖岙、上马鞍、下马鞍、三盘尾、柴屿)采捕活体样品共44只, 于–30 °C冷冻保存。

1.2 形态学指标测量

选取其中胴体完整、腕无残缺的样品37只(15雄、22雌)进行形态学测量。形态学测量指标及方法参见Roper等(1983)、Norman等(1997b)及Liao等(2009), 测量指标如下: 全长(TL, total length)、体重(TW, total weight)、胴背长(ML, mantle length)、胴背宽(MW, mantle width)、胴腹长(VML, ventral mantle length)、头宽(HW, head width)、外侧漏斗长(FL, funnel length)、内侧漏斗长(FFL, free funnel length)、腕间膜深(WD, web depth)、腕间膜式(WF, web formula)、腕长(AL, arm length)、腕式(AF, arm formula)、腕宽(AW, arm width)、右三腕吸盘数(SNR3, sucker number on 3rd arm of right side)、左三腕吸盘数(SNL3, sucker number on 3rd arm of left side)、吸盘直径(SD, sucker diameter)、舌叶长(LL, ligula length)、交接基长(CaL, calamus length)。

1.3 形态学数据处理与分析

用于形态学分析的5个群体数据见表1。除本文的南麂列岛数据, 其他中华蛸与真蛸数据引自Amor等(2017b)。

为消除不同规格个体对特征指标的影响, 对数据进行标准化处理, 标准化指数定义如下:

表1 形态多元分析的样品信息

Tab.1 Sample information for morphological multivariate analysis

注: 平均胴背长±标准差的单位为mm

ALLI: 左腕长/胴背长×100; ALRI: 右腕长/胴背长×100; AWI: 腕宽/胴背长×100; CaLI: 交接基长/舌叶长×100; FLI: 外侧漏斗长/胴背长×100; FFLI: 内侧漏斗长/外侧漏斗长×100; HWI: 头宽/胴背长×100; LLI: 舌叶长/茎化腕长×100; MWI: 胴背宽/胴背长×100; OAI: 茎化腕长/左三腕长×100; SDIn: 普通吸盘直径/胴背长×100; SDIe: 扩大吸盘直径/胴背长×100; WDI: 腕间膜深/最长腕长×100。

主成分分析 根据雌性19个指标ML、VML、MWI、HWI、FLI、FFLI、ALLI1、ALLI2、ALLI3、ALLI4、ALRI1、ALRI2、ALRI3、ALRI4、AWI、SDIn、SNL3、SNR3、WDI, 雄性22个指标ML、VML、MWI、HWI、FLI、FFLI、WDI、ALLI1、ALLI2、ALLI3、ALLI4、ALRI1、ALRI2、ALRI3、ALRI4、AWI、SDIn、SDIe、SNL3、SNR3、LLI、CaLI, 计算出互不关联的主成分。主成分贡献率和累计贡献率的计算方法参照Brzeski等(1988)。

判别分析 采用逐步判别法(SPSS 25.0)对雌性和雄性指标(同主成分)进行判别分析, 建立5个群体的判别函数。利用判别函数得分绘制前两个判别函数系数散点图, 判别准确率的计算公式为:

判别准确率1(%)=判别正确的个体数/该种群个体数×100%, (1)

判别准确率2(%)=判别正确的个体数/判入该种

群个体数×100%, (2)

式中,A和B分别为第个群体中被判别正确的个体数和实际判别的个体数,为群体数。

聚类分析 取雌性和雄性各群体指标(同主成分)的平均值, 采用欧式距离最短系统聚类法进行聚类分析(SPSS 25.0)。

1.4 COI序列分析

分别从采集样品的胴体部剪取少许肌肉组织保存于无水乙醇中, 采用CTAB法提取总DNA, 用引物Primer-F: 5′-GGTCAACAAATCATAAAGATATTGG-3′,Primer-R: 5′-TAAACTTCAGGGTGACCAAAAAATCA-3′ (Folmer, 1994)扩增片段, PCR反应体系: H2O 19 μL, Mix酶25 μL, F与R引物各2.5 μL, DNA 1 μL。PCR反应程序: 94 °C 4 min, 然后94 °C 10 s, 55 °C 20 s, 72 °C 30 s, 共32个循环, 最后72 °C 5 min延伸。扩增产物经琼脂糖凝胶电泳检测, 获得23个体扩增产物用于双向测序。

测序完成后, 用DNASTAR软件包中的SeqMan对测序结果进行拼接, 使用MEGA X (Kumar, 2018)中的Clustal W将全部序列进行多重比对分析, 序列已上传至GenBank (序列号: OK001740～ OK001762)。利用DnaSP v5 (Librado, 2009)计算群体的遗传多样性参数, 并将23条序列与GenBank下载的真蛸()、中华蛸()、、、、cf.的序列(表2)构建系统发育树, 幽灵蛸(AB385880)作为外群。建树前, 用ModelFinder (Kalyaanamoorthy, 2017)选择构建系统发育树的最佳模型, 用IQ-TREE (Nguyen, 2015)构建最大似然系统发育树。在计算真蛸复合体各物种间的遗传距离时, 利用MEGA X (Kumar, 2018)软件选择Kimura-2-Parameter (K2P)模型, 其他参数设置为默认值。

2 结果

2.1 形态学鉴别特征

南麂列岛样品为中等至大型章鱼, 成熟个体胴背长117.5～133.9 mm, 全长593.1～998.0 mm, 最高体重达2 140.0 g; 活体皮肤呈灰色或红棕色, 皮肤表面有纹理; 胴背长大于胴体宽(MWI 60.3～75.0), 头部窄于胴体部(HWI 30.1～65.7); 眼上具乳头状突起; 漏斗中等长(FLI 30.8～58.2), 内侧漏斗长多为外侧漏斗长的70%左右(FFLI 51.5～80.6), 腕长中等(ALI 218.1～531.1), 腕式通常为2>3>4>1, 最长腕长约为胴背长的4倍, 最短腕长约为胴背长的3倍; 雄性右三腕为茎化腕, 成熟雄性茎化腕多为左三腕的80%(OAI 78.1～96.4), 茎化腕吸盘数范围123～151; 腕宽中等(AWI 14.6～29.5), 腕间膜中等深度(WDI 15.9～26.2), 腕间膜式通常为C>D>B>E>A; 吸盘中等大小(SDIn 6.0～13.9), 雄性在第二、第三对腕的第12～15吸盘位置具有1～2个扩大吸盘; 舌叶小, 呈锥形, 中央具有沟槽, 多为茎化腕长的1%左右(LLI 0.6～1.2), 交接基多为舌叶的50% (CaLI 34.2～65.5)。

表2 用于构建最大似然系统发育树序列

Tab.2 COI sequences used for the construction of maximum likelihood phylogenetic tree

2.2 主成分分析

雌性主成分分析共构建5个主成分, 主成分的负荷值和贡献率见表3, 主成分1的贡献率为46.565%, 主成分2的贡献率为13.152%, 主成分3贡献率为9.005%, 主成分4贡献率为6.164%, 主成分5贡献率为5.800%, 累计贡献率为80.686%。在主成分1中, 影响群体形态差异的指标为HWI、ALLI1、ALLI2、ALLI3、ALLI4、ALRI1、ALRI2、ALRI3、ALRI4、SDIn、SNL3、SNR3, 主成分2中影响群体差异的指标为ML、VML, 主成分3中影响群体差异的指标为FFLI和WDI, 主成分4中影响群体差异的指标为MWI,主成分5中影响群体差异的指标为FLI和AWI。

雄性主成分分析共构建4个主成分, 主成分的负荷值和贡献率见表4, 主成分1贡献率为47.060%, 主成分2贡献率为13.121%, 主成分3贡献率为9.461%, 主成分4贡献率为6.151%, 累计贡献率为75.794%。在主成分1中, 影响群体形态差异的指标为HWI、FLI、ALLI1、ALLI2、ALLI3、ALLI4、ALRI1、ALRI2、ALRI3、ALRI4、AWI、SDIn、SDIe、SNL3、SNR3, 主成分2中影响群体差异的指标为ML、VML、LLI, 主成分3中影响群体差异的指标为MWI、FFLI、WDI, 主成分4中影响群体差异指标为CaLI。

表3 雌性群体形态特征主成分的负荷值和贡献率

Tab.3 Contribution and load of principal components on morphological characteristics of female populations

雌、雄各群体主成分1和主成分2的散点图如图1所示。在雌性各群体中, 西班牙真蛸(SP)均分布在轴右侧, 可以跟其他群体分开, 而南麂列岛(NJLD)与中华蛸各群体(ND、YL和JZ)存在重叠; 在雄性各群体中, 西班牙真蛸(SP)分布于轴右侧, 其他四个群体(NJLD、ND、YL和JZ)集中分布在轴左侧并存在部分重叠。

2.3 判别分析

利用逐步判别分析法分别对雌、雄群体进行费歇尔判别分析。

雌性各群体判别公式如下:

南麂列岛:1=4.0021+0.5052–78.379

宁德:2=3.3941+0.4802–66.822

宜兰:3=2.7891+0.4892–63.418

九州:4=4.1231+0.5722–95.204

表4 雄性群体形态特征主成分的负荷值和贡献率

Tab.4 Contribution and load of principal components on morphological characteristics of male populations

西班牙:5=5.6051+0.6192–126.995

判别式中,1为SDIn、2为SNR3。

雄性各群体判别公式如下:

南麂列岛:1=–0.0601+6.0932+0.7163+1.7274– 151.038

宁德:2= –0.1211+4.6342+2.428

3+1.6284–129.550

宜兰:3=–0.0091+2.6602+1.3203+1.5144–113.413

九州:4=–0.0731+4.7372+1.5933+1.6784–137.916

西班牙:5=–0.0811+3.2552+3.2663+2.1774–225.186

判别式中,1为ALRI1、2为SDIn、3为SDIe、4为SNR3。

判别分析散点图见图2, 西班牙真蛸(SP)可与其余四个群体明显区分开, 而南麂列岛(NJLD)与中华蛸各群体(ND、YL和JZ)存在重叠。为检验判别结果, 对所有雌、雄群体进行预测分类(见表5和表6), 雌、雄综合判别分析结果分别为66.0%与84.2%, 其中, 雌、雄真蛸(SP)均被100%准确判别。

图1 主成分分析散点图

注: 图例中字母表示各群体。NJLD: 南麂列岛, ND: 宁德, YL: 宜兰, JZ: 九州, SP: 西班牙

图2 判别分析散点图

注: 图例中字母表示各群体。NJLD: 南麂列岛, ND: 宁德, YL: 宜兰, JZ: 九州, SP: 西班牙

表5 雌性群体判别分析结果

2.4 聚类分析

对雌、雄5个群体的所有样本校正值进行聚类, 图3结果显示5个群体聚为2大类, 西班牙真蛸单独成为一支, 其余四个群体聚成一支。在雌性聚类分析中, 南麂列岛(NJLD)、宁德(ND)、宜兰群体(YL)先聚成一支, 再与九州群体(JZ)群体聚在一起。在雄性聚类分析中, NJLD、YL与JZ群体先后聚在一起, 再与ND群体聚成一支。结果表明, NJLD与中华蛸亲缘关系更近, 与真蛸存在明显差异。

表6 雄性群体判别分析结果

Tab.6 Discriminant results of male populations

图3 聚类分析图

注: 图例中字母表示各群体地理位置。NJLD: 南麂列岛, ND: 宁德, YL: 宜兰, JZ: 九州, SP: 西班牙

2.5 COI序列分析

南麂列岛群体的单倍型数(hap) 4个, 单倍型多样性(d)为0.320±0.121, 多态位点()为6个, 平均核苷酸差异数()为0.672, 核苷酸多样性指数()为0.001 11。南麂列岛群体与中华蛸的K2P遗传距离为0.14%, 而与真蛸复合种其他类型的遗传距离为2.96%～12.11% (表7), 这一结果支持南麂列岛采集样品为中华蛸。

系统发育树(图4)显示南麂列岛中华蛸与其他中华蛸序列聚成一支(序列OK001749除外), 这一支包括日本与中国沿海样本, 但未见明显遗传分化。其余真蛸物种复合体各类型均单独聚成一支,、cf.与、各聚为一支, 互为姊妹群, 然后与中华蛸聚在一起,亲缘关系较远, 是其他五个真蛸复合种的姊妹群。

表7 K2P模型下基于基因的种间遗传距离

Tab.7 Interspecies genetic distance under K2P model based on COI gene

图4 基于COI基因构建的ML系统发育树

注: 基于TIM2+F+G4进化模型, 节点数字为自举值

3 讨论

3.1 形态学分析

物种形态多元分析是划分不同群体、判断亲缘关系的有效方法, 物种形态学指标越多, 包含的鉴别特征信息就越多。头足类形态结构复杂, 可测量的形态学指标多达几十个。高晓蕾等(2019)利用形态多元分析方法表明中国沿海11个长蛸群体在形态上存在一定的地域差异, 陈唯(2018)研究了我国沿海短蛸6个群体的形态学特征, 结果显示群体间存在明显地理分化。在本研究中, 主成分分析、判别分析和聚类分析均将真蛸(SP)群体分离开来, 南麂列岛(NJLD)与其他三个中华蛸群体(ND、YL和JZ)存在重叠, 但雌性九州(JZ)与雄性宁德群体(ND)仍与其他中华蛸群体分离(图1～图3)。九州(JZ)中华蛸群体与中国沿海群体相距较远, 不同的生态与水文条件是导致其与中国沿海群体存在形态学差异的主要因素, 而宁德近海多港湾, 湾内风浪较小, 海水流速适中, 适宜海洋生物的生长, 这可能是造成该地区雄性中华蛸与其他群体存在形态差异的原因。

在头足类的物种鉴定中, 性别特征是重要的分类依据(Bello, 1995; Brakoniecki, 1996; Von Byern, 2010), 相较于其他可测量性状, 与性别相关的形态特征具有更加多样化的特点(Pomiankowski, 1995; Norman, 1997a; O’Dor, 1998), 雄性往往具有茎化腕、舌叶与扩大吸盘等特征, 因此, 基于形态学的物种多样性分析, 雄性的形态学差异分析结果可能更加可靠。本研究判别分析中, 雄性的综合判别分析准确率高于雌性(表5和表6), 雄性判别式共有4个参数, 分别为ALRI1、SDIn、SDIe与SNR3, 其中SDIe与SNR3为与雄性性别密切相关的形态指标。

Gleadall(2016)比较了中华蛸与真蛸的形态学异同, 重描述了日本九州与本州沿海采集的所谓“真蛸”, 恢复了中华蛸d’Orbigny, 1841种名的有效性。本文总结了南麂列岛章鱼的形态特征, 并将其与日本中华蛸标本(Gleadall, 2016)、真蛸(Norman, 2014; Gleadall, 2016)进行了比较(表8), 指出中华蛸主要鉴别特征有三点: (1) 成熟雄性在第二对腕和第三对腕的第12～15个吸盘之间存在1～2个扩大吸盘(通常在第13和14个吸盘); (2) 茎化腕长约为左三腕的80%, 吸盘数量范围119～152个; (3) 末端舌叶较真蛸小。Toll(1988)认为, 茎化腕吸盘数是鉴定蛸类的重要指标, 雄性茎化腕吸盘数量在同一物种中相对稳定, 不同种之间多存在差异, 南麂列岛样品的茎化腕吸盘数量与真蛸明显不同, 其数量更少。根据以上形态学数据分析, 南麂列岛样品与Gleadall(2016)重描述的中华蛸为同一物种, 二者与真蛸存在明显差异。因此, 从形态上可以确定南麂列岛章鱼为中华蛸。任静等(2021)在浙江南麂列岛与福建连江样品的肠道与盲囊中发现了一种新的寄生虫——多刺丛集球虫, 该寄生虫与真蛸中发现的寄生虫在形态与分子上均存在很大差异, 由于蛸类丛集球虫具有宿主特异性, 因此, 从寄生虫角度也可佐证南麂列岛章鱼与真蛸为不同种。

表8 主要形态鉴别特征比较

Tab.8 Comparison of morphological identification characteristics

3.2 COI序列分析

利用线粒体序列发现南麂列岛章鱼的hap为4个,d为0.320±0.121,为0.001 11显示其遗传多样性水平偏低。对于此结果, 需通过基因、微卫星等分子标记进一步验证, 同时建议当地有关部门对南麂列岛物种采取一定程度的保护措施, 如每年的5～7月为中华蛸繁殖期, 可设置禁捕或限捕期。

系统发育树显示序列OK001749未与其他中华蛸聚成一支(图4), Avendaño等(2020)用序列构建了真蛸物种复合体系统发育树, 发现墨西哥采集的7号样本与真蛸聚成一支,而用构建的系统发育树显示所有聚成一支, 两种分子标记存在差异, 这可能也是序列OK001749单独分离的原因, 因此, 今后可用多基因联合构建更可靠的系统发育树。

南麂列岛章鱼和其他中华蛸的K2P遗传距离仅为0.14%, 与真蛸复合体其他类型的遗传距离均大于2% (表7), 支持当地采集样品为中华蛸而非真蛸。此外, 中华蛸并未出现遗传分化(图4), 说明不同地理群体间存在一定的基因交流, 可能与其浮游幼虫阶段随黑潮等洋流远距离扩散有关。

4 结论

本研究表明南麂列岛采集的章鱼与真蛸在形态特征和遗传上均存在显著差异, 而与中华蛸鉴别特征一致, 支持该海域采集样品的正确种名为中华蛸, 为南麂列岛世界生物圈保护区种质资源保护提供了精准物种信息。

致谢 澳大利亚维多利亚博物馆卢重成教授对文章撰写提供了宝贵意见, 中国海洋大学贝类遗传育种研究室硕士研究生胡元海协助实验, 谨致谢忱。

叶鹏, 蔡厚才, 庄定根, 等, 2006. 南麂海区野生贝类增养殖种类初步筛选[J]. 渔业现代化(4): 26-28.

冯雪, 陈丕茂, 秦传新, 等, 2013. 不同温度和体质量对南海野生真蛸呼吸和排泄的影响[J]. 中国水产科学, 20(5): 968-974.

任静, 何卫彤, 王丽华, 等, 2021. 中华蛸寄生丛集球虫(顶复亚门: 丛集球虫科)新种的形态学与分子生物学研究[J]. 海洋与湖沼, 52(5): 1323-1331.

刘兆胜, 刘永胜, 郑小东, 等, 2011. 不同饵料对真蛸亲体产卵量、受精卵孵化率及初孵幼体大小的影响[J]. 海洋科学, 35(10): 81-85.

孙田田, 苏永全, 洪婧妮, 等, 2012. 真蛸热休克蛋白90基因(HSP90)的克隆及表达[J]. 水产学报, 36(9): 1367-1375.

肖懿哲, 姚成杰, 朱友芳, 等, 2019. 真蛸FAXDC2基因的克隆及其表达分析[J]. 海洋科学, 43(8): 56-63.

陈唯, 2018. 中国沿海短蛸()谱系地理格局与适应性分化研究[D]. 舟山: 浙江海洋大学: 14-19.

郑小东, 刘兆胜, 赵娜, 等, 2011. 真蛸()胚胎发育及浮游期幼体生长研究[J]. 海洋与湖沼, 42(2): 317-323.

高晓蕾, 许然, 张志新, 等, 2019. 中国沿海长蛸群体形态性状的差异[J]. 水产学报, 43(7): 1593-1602.

董正之, 1988. 中国动物志[M]. 北京: 科学出版社: 182-184.

蔡厚才, 庄定根, 叶鹏, 等, 2007. 浙江南麂岛真蛸网箱和水泥池养殖试验[J]. 南方水产, 3(2): 66-70.

蔡厚才, 庄定根, 叶鹏, 等, 2009. 真蛸亲体培育、产卵及孵化试验[J]. 海洋渔业, 31(1): 58-65.

AMOR M D, LAPTIKHOVSKY V, NORMAN M D,, 2017a. Genetic evidence extends the known distribution ofto the mid-Atlantic islands Ascension and St Helena [J]. Journal of the Marine Biological Association of the United Kingdom, 97(4): 753-758.

AMOR M D, NORMAN M D, ROURA A,, 2017b. Morphological assessment of thespecies complex evaluated in light of molecular-based phylogenetic inferences [J]. Zoologica Scripta, 46(3): 275-288.

AVENDAÑO O, ROURA Á, CEDILLO-ROBLES C E,, 2020.: a cryptic species of thespecies complex redescribed from the Caribbean [J]. Aquatic Ecology, 54(4): 909-925.

党的十八大就扎实推进社会主义文化强国建设做出了全面部署，明确指出“让人民享有健康丰富的精神文化生活，是全面建成小康社会的重要内容。要坚持以人民为中心的创作导向，提高文化产品质量，为人民提供更好更多精神食粮”。由此可见，“精神富有”赋予文艺界以新的使命，因此，承前启后，创新手段，扩大覆盖；提升标准，出人出作品；持久努力，建树品牌，是文艺创新的题中之意，是我们应当为之努力做出的一种姿态。

BELLO G, 1995. A key for the identification of the Mediterranean sepiolids (Mollusca: Cephalopoda) [J]. Bulletin de l’Institut Océanographique, Monaco, n° spécial 16: 41-55.

BRAKONIECKI T F, 1996. A revision of the genusVoss, 1953 (Cephalopoda; Myopsida) [J]. Bulletin of Marine Science, 58(1): 9-28.

BRZESKI V J, DOYLE R W, 1988. A morphometric criterion for sex discrimination in tilapia [C] // The Second International Symposium on Tilapia in Aquaculture. ICLARM Conference Proceedings 15. Department of Fisheries, Bangkok, Thailand, and International Center of Living Aquatic Resources Management, Manila, Philippines: 439-444.

D’ORBIGNY A, 1835-1848. Histoire naturelle générale et particulière des céphalopodes acétabulifères [M] // DE FÉRUSSAC A E, D’ORBIGNY A. Histoire Naturelle Générale et Particulière des Céphalopodes Acétabulifères Vivants et Fossiles, Paris, Lacour.

FOLMER O, BLACK M, HOEH W,, 1994. DNA primers for amplification of mitochondrial cytochromeoxidase subunit I from diverse metazoan invertebrates [J]. Molecular Marine Biology and Biotechnology, 3(5): 294-299.

GLEADALL I G, 2016.d’Orbigny, 1841 (Cephalopoda: Octopodidae): valid species name for the commercially valuable East Asian common octopus [J]. Species Diversity, 21(1): 31-42.

KALYAANAMOORTHY S, MINH B Q, WONG T K F,, 2017. ModelFinder: fast model selection for accurate phylogenetic estimates [J]. Nature Methods, 14(6): 587-589.

LI F H, BIAN L, GE J L,, 2020. Chromosome-level genome assembly of the East Asian common octopus () using PacBio sequencing and Hi-C technology [J]. Molecular Ecology Resources, 20(6): 1572-1582.

LIAO J X, Lu C C, 2009. A new species of(Cephalopoda: Octopodidae) from Taiwan and morphology of mucous pouches [J]. Journal of Molluscan Studies, 75(3): 502-509.

LIBRADO P, ROZAS J, 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data [J]. Bioinformatics, 25(11): 1451-1452.

NGUYEN L T, SCHMIDT H A, VON HAESELER A,, 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies [J]. Molecular Biology and Evolution, 32(1): 268-274.

NORMAN M D, 2000. Cephalopods: A World Guide [M]. Hackenheim, Germany: ConchBooks Press: 268.

NORMAN M D, FINN J K, HOCHBERG F G, 2014. Family octopodidae [M] // JEREB P, ROPER C F E, NORMAN M D,. Cephalopods of the World. An Annotated and Illustrated Catalogue of Cephalopod Species Known to Date. Vol 3 Octopods and Vampire Squids. FAO Species Catalogue for Fishery Purposes. Rome, Italy: FAO Press: 370.

NORMAN M D, HOCHBERG F G, 2005. The current state of octopus taxonomy [J]. Phuket Marine Biological Center Research Bulletin, 66: 127-151.

NORMAN M D, LU C C, 1997a. Redescription of the southern dumpling squidand a revision of the genus(Cephalopoda: Sepiolidae) [J]. Journal of the Marine Biological Association of the United Kingdom, 77(4): 1109-1137.

NORMAN M D, SWEENEY M J, 1997b. The shallow-water octopuses (Cephalopoda: Octopodidae) of the Philippines [J]. Invertebrate Taxonomy, 11(1): 89-140.

O'DOR R K, LIPINSKI M R, 1998. The genus(Cephalopoda; Ommastrephidae): characteristics, distribution and fisheries [M] // RODHOUSE P G, DAWE E G, O'DOR R K. Squid Recruitment Dynamics. Rome: FAO: 273.

POMIANKOWSKI A, MOLLER A P, 1995. A resolution of the lek paradox [J]. Proceedings of the Royal Society B: Biological Sciences, 260(1357): 21-29.

ROPER C F E, VOSS G L, 1983. Guidelines for taxonomic descriptions of cephalopod species [J]. Memoirs of Museum Victoria, 44: 48-63.

SAKAGUCHI H, HAMANO T, NAKAZONO A, 2000. Population structure ofestimated from catch size composition in northeastern iyo-nada of the seto Inland Sea, Japan [J]. Bulletin of the Japanese Society of Fisheries Oceanography, 64(4): 224-234.

SAKAGUCHI H, HAMANO T, NAKAZONO A, 2005. Growth ofin the northeastern iyo-nada of the seto Inland Sea, Japan [J]. Bulletin of the Japanese Society of Fisheries Oceanography, 66(1): 11-15.

SASAKI M, 1929. A monograph of the dibranchiate cephalopods of the Japanese and adjacent waters [J]. Journal of the College of Agriculture, Hokkaido Imperial University, 20(supplement): 1-357.

TOLL R B, 1988. The use of arm sucker number in Octopodid systematics (Cephalopoda, Octopoda) [J]. American Malacological Bulletin, 6(2): 207-211.

VAN NIEUWENHOVE A H M, RATSIMBAZAFY H A, KOCHZIUS M, 2019. Cryptic diversity and limited connectivity in octopuses: recommendations for fisheries management [J]. PLoS One, 14(5): e0214748.

VON BYERN J, KLEPAL W, 2010. Re-evaluation of taxonomic characters of(Cephalopoda, Mollusca) [J]. Malacologia, 52(1): 43-65.

WARNKE K, SÖLLER R, BLOHM D,, 2004. A new look at geographic and phylogenetic relationships within the species group surrounding(Mollusca, Cephalopoda): indications of very wide distribution from mitochondrial DNA sequences [J]. Journal of Zoological Systematics and Evolutionary Research, 42(4): 306-312.

MORPHOLOGICAL AND GENETIC DIVERSITY ANALYSIS OFIN NANJI ISLANDS

LI Jia-Hua1, 2, CHEN Shun3, CHEN Wan-Dong3, XIE Shang-Wei3, NI Xiao-Pin3, ZHENG Xiao-Dong1, 2

(1. Institute of Evolution & Marine Biodiversity, Ocean University of China, Qingdao 266003, China; 2. Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao 266003, China; 3.Nanji Islands National Marine Nature Reserve Administration, Wenzhou 325400, China)

Forty-four medium- to large-sized octopuses were collected from Nanji Islands in Zhejiang Province from June to December 2020, and the main morphological characteristics were described. The morphological diversity of this octopus and its differences betweenwere investigated by multivariate analysis. The cytochrome c oxidase subunit I () sequences were used to calculate the genetic diversity parameters and to construct the maximum likelihood phylogenetic tree. In addition, the Kimura-2-Parameter (K2P) genetic distance was also calculated. Results showed thatthis octopus was in line with the redescribedin morphological features, while they had some significant differences with. First, the enlarged suckers on arms 2 and 3 of males were between the 12thand the 15thin population from Nanji Islands, andwerebetween the 15thand the 19thinstead. Secondly, the number of suckers on the hectocotylus of the Nanji Islands group was less than that in. Thirdly, the male samples from Nanji Islands had a shorter ligula than. In the principal component and discriminant analysis,can be clearly distinguished from the other 4 populations with the discrimination ratio ofreaching 100%, while there was overlap among the other 4 groups. And the cluster analysis showed that the population of Nanji Islands was closely related to, both of which had a large distance from.In genetic diversity analysis, there were four haplotypes in Nanji Islands population, and the diversity of haplotypes (d) is 0.320±0.121. Besides, six polymorphic sites were examined, and the nucleotide diversity was 0.001 11. The phylogenetic tree showed that the octopus of Nanji Islands had the closest relationship with. Moreover, the K2P genetic distance between Nanji Islands population andwas 0.14%, while among Nanji Islands population andspecies complex ranged from 2.96%～12.11%. Therefore, samples from Nanji Islands can be determined asby the morphological and molecular analysis.

; Nanji Islands; morphological diversity; multivariate analysis;gene

*南麂列岛国家级海洋自然保护区院士专家工作站研究项目, NJKJ-2019-003号; 国家自然科学基金项目, 32170536号。李嘉华, 硕士研究生, E-mail: jhli123@163.com

郑小东, 教授, E-mail: xdzheng@ouc.edu.cn

2021-09-29,

2021-11-16

Q959.216

10.11693/hyhz20210900225