我国新记录绿潮物种Ulva laetevirens 的比较叶绿体基因组学研究*

王 静 王毅超,4 王洪淑 刘 峰① 陈楠生①

(1. 中国科学院海洋研究所 海洋生态与环境科学重点实验室 青岛 266071; 2. 青岛海洋科学与技术试点国家实验室海洋生态与环境科学功能实验室 青岛 266237; 3. 中国科学院海洋大科学研究中心 青岛 266071; 4. 中国科学院大学 北京 100049)

绿潮(green tide)在全球范围内频繁暴发, 欧洲、北美、澳洲和亚洲等海域均有报道。欧洲北亚得里亚海(Runcaet al, 1996)、波罗的海(Bäcket al, 2000)、德国瓦登海(Kolbeet al, 1995)等, 美国罗得岛州(Guidoneet al, 2013)、华盛顿州(Nelsonet al, 2003),亚洲中国(Liuet al, 2013)、日本(Shimadaet al, 2003)、菲律宾(Largoet al, 2004)等国的海域, 均发生过大规模的绿潮事件。近些年, 全球范围内绿潮呈现出常态化的暴发趋势, 对经济和社会活动造成了严重的负面影响, 对致灾海域的生态环境、群落结构、生物多样性以及水产养殖也造成极大的破坏性影响(Lyonset al, 2014; Liuet al, 2020)。绿潮暴发的原因种主要包括绿藻门石莼科的石莼属(Ulva)(Smetaceket al,2013)、硬毛藻属(Chaetomorpha)(De Paula Silvaet al,2008)和刚毛藻属(Cladophora)(De Paula Silvaet al,2013), 其中石莼属绿藻占主要部分(Fletcher, 1996;Smetaceket al, 2013; Fortet al, 2020)。石莼属绿潮物种中, 部分物种能够形成跨境侵袭的大规模绿潮灾害, 譬如中国黄海绿潮; 部分仅在局部海域暴发, 譬如美国纳拉甘塞特海湾连年发生的扁浒苔(Ulva compressa)-硬石莼(U. rigida)绿潮(Guidoneet al,2013)、韩国济州岛东海岸绿潮(Leeet al, 2019), 以及我国渤海自2015 年连续暴发的浒苔绿潮(Songet al,2019)。我国黄海浒苔绿潮在2007 年暴发之前, 浒苔物种或绿潮事件在我国并未引起广泛关注, 暴发之后投入大量的科研工作致力于浒苔绿潮的相关研究,这恰恰警示我们对潜在的绿潮物种和局地规模的绿潮事件应给予高度关注。调查和了解石莼属绿潮物种生物地理分布和局地规模的绿潮事件是预防和监测大规模绿潮暴发的首要工作之一。

我国在局地绿潮方面的研究报道相对较少。马元庆等(2010)报道了2007 年烟台金沙滩海域石莼(U.lactuca)绿潮事件相关的理化条件; Song 等(2019)研究发现秦皇岛海域自2015 年连年暴发的绿潮原因种是羽藻(Bryopsis plumosa)和浒苔(U. prolifera), 经巡航观测和对沿岸定生样本的分子标记rbcL分析发现,这两个原因种只在海草床中被发现, 最终确定海草床中的羽藻和浒苔是渤海绿潮的原因种。总体而言,我国在局地绿潮的相关研究工作正处于起步阶段。

我国海域石莼属物种的记录自20 世纪60 年代由曾呈奎等生态学前辈系统整理开始, 一直处于不断地补充和修正中。1963 年董美龄(1963)报道我国海域的“浒苔属”(Enteromorpha)绿藻有11 个种。2015年丁兰平等(2015)整理的中国海洋绿藻门新分类系统中“浒苔属”有10 个种, “石莼属”有14 个种。2020年Xie 等(2020)报道了我国海域的6 个新记录物种U.simplex、U. splitiana、U. partita、U. meridionalis、U. tepida和U. chaugulli, 后三者为绿潮物种。此外,还有一株淡水石莼物种U. shanxiensis(Chenet al,2015) 和一株浒苔亚种U. proliferasubsp.Qingdaoensis(Cuiet al, 2018)被鉴定。这说明我国海域可能还有更多的石莼属物种(包括石莼属绿潮物种)有待发现。由于石莼属部分近缘物种种间形态十分相近, 并且具有形态可塑性, 基于形态进行石莼属物种鉴定难度较大。随着分子生物技术的发展, 分子标记的开发和利用对石莼属物种的鉴定起了重要的推动作用(Haydenet al, 2003)。目前, 应用于石莼属物种鉴定的通用单基因分子标记是ITS、rbcL和tufA(Garyet al, 2010; Mudassaret al, 2016; Cuiet al, 2018)。然而,尽管这些通用分子标记使用广泛, 但是由于它们的分辨率比较有限, 在部分情况下无法做到准确定种。相对单基因分子标记, 细胞器基因组可以作为串联基因标记提供更可靠的分类依据。

本研究针对2020 年11 月初烟台第一海水浴场局地绿潮事件, 采用分子标记和叶绿体基因组比较分析对此次绿潮事件的原因种进行物种鉴定和分析,确定该局地绿潮原因种为U. laetevirens, 系该物种在我国海域的首次绿潮记录。同时, 探索了三者之间的遗传进化关系。

1 材料与方法

1.1 样品采集与处理

2020 年11 月2 日于山东省烟台市第一海水浴场的石莼绿潮事件(121°24′59″E, 37°32′6″N, 图1a)中随机采集了11 株石莼属绿藻, 分别编号CNS00533、CNS00592—CNS00601。现场采集的藻体置于内含无菌海水的离心管, 无菌海水中含有f/2 培养基、0.5‰GeO2、50 µg/mL 敌百虫和五种复合抗生素(50 µg/mL链霉素, 66.6 µg/mL 庆大霉素, 20 µg/mL 环丙杀星,2.2 µg/mL 氯霉素, 100 µg/mL 氨苄青霉素)(Shiblet al,2020)。采样后低温保存, 24 h 内运输回实验室。藻体样本到达实验室后, 立即用无菌海水冲洗藻体表面。最后, 将藻体置于含有30 mL 培养液的培养皿中, 内含f/2 培养基、0.5‰ GeO2、50 µg/mL 敌百虫和上述五种复合抗生素, 置于光照培养箱进行培养, 培养条件18 °C, 100—120 μmol photons/(m2·s), L:D = 12 h:12 h。待镜检无其他微藻、浮游动物等杂质, 即可用于后续实验。

1.2 DNA 文库制备与测序

藻体DNA 采用植物基因组DNA 提取试剂盒(DP305, 北京天根)进行提取, 所提取的 DNA 经Covaris S220 超声破碎仪进行片段化, 获得读长为350 bp 的DNA 序列; 经末端修复、加A 尾、加测序接头、纯化、PCR 及产物富集。文库质量分别采用NGS3K/Caliper 和 real-time PCR (Qubit 3.0 fluorometer, 美国Invitrogen)进行文库大小分布评估和定量分析。质量合格的文库采用美国 Illumina NovaSeq PE150 (北京诺禾致源科技股份有限公司)进行测序。

1.3 分子标记tufA 的序列获取及系统发育树构建

烟台绿潮11 个藻体样本的分子标记tufA序列信息通过软件BWA v0.7.17 (Liet al, 2010)、samtools v1.9 (Liet al, 2009)和SPAdes v3.14.0 (Bankevichet al,2012)从上述全基因组测序数据中筛选组装获得, 所采用的参考序列来自于U. laetevirens叶绿体基因组MT179351 注释的tufA基因全长序列。将11 个样品和来自GenBank 的tufA序列通过Clustal W 进行多重比对分析(Thompsonet al, 1997), 并辅以手工校正。通过Mega 7.0 软件的最大似然法构建进化树(Kumaret al, 2016)。Bootstrap 值为1 000。

1.4 叶绿体基因组组装与注释

采用GetOrganelle 软件(Jinet al, 2020)从头组装获得叶绿体基因组序列, 其中使用的组装软件为SPAdes (3.10.1)(Bankevichet al, 2012), 分析过程中使用的seed 参考序列为NCBI 登录序列MT179351(Ulva rigida叶绿体基因组)(Fortet al, 2021)。组装后的叶绿体基因组序列通过BWA v0.7.17 软件的MEM运算法对所构建的叶绿体基因组序列进行质检(Liet al, 2010), 并通过IGV v2.8.12 进行可视化比对分析(Robinsonet al, 2011)。本研究采用在线软件MFannot(https://megasun.bch.umontreal.ca/RNAweasel/)和ORF finder (https://www.ncbi.nlm.nih.gov/orffinder)对叶绿体基因组进行注释。

1.5 共线性分析

将CNS00533 分别与U. laetevirensMT179351、U. rigidaMT179352 和U. australisMT179348 进行共线性分析, 采用circos-0.69 进行作图(Krzywinskiet al, 2009)。

1.6 遗传进化分析

对19 株石莼属绿藻叶绿体基因组的71 个蛋白编码基因的氨基酸序列进行遗传进化分析。氨基酸序列采用MAFFT(默认参数)比对(Katohet al, 2013), 然后通过PhyloSuite 进行串联(Smithet al, 2008), 以基因的位置对串联后的序列进行分区, 每个分区的最适氨基酸替代模型通过 ModelFinder 确定(Kalyaanamoorthyet al, 2017), 以IQ-TREE(默认参数)构建基于氨基酸序列的遗传进化树(Trifinopouloset al, 2016), bootstrap 值设置为1 000。硅藻Lithodesmium undulatum(KC509525)和Odontella sinensis(Z67753)的叶绿体基因组作为外类群。

2 结果

2.1 2020 年11 月烟台局地绿潮原因种的鉴定

2020 年11 月2 日于山东省烟台市第一海水浴场暴发了局地绿潮, 浴场沿岸水体充斥绿色藻体碎片,在海浪的作用下大量的藻体积留海滩(图1b)。该绿藻藻体呈叶片状, 多孔、双层基质; 因藻体破碎无法辨别不同部位且未见假根, 难以从形态水平鉴别物种(Maoet al, 2014)(图1c)。

为鉴定此次烟台绿潮物种, 随机采集了12 株藻体样本, 对其构建了基于tufA的单基因遗传进化树(图 1d)。结果显示, 其中 11 株 CNS00533、CNS00592—CNS00601 序列完全一致, 且与U. rigida和U. laetevirens的参考序列聚类在一个clade 中; 在NCBI 数据库中, 11 株藻的tufA全长序列(1 224 bp)与U. laetevirens(MT179351)(Fortet al, 2021)的tufA序列(1 224 bp) 100%相同。另外1 株CNS00602 与U.australis参考序列聚类在一个 clade 中, 且与U.australisMT160640 (Fortet al, 2021)完全一致。综上所述, 除CNS00602 株系以外的11 株可以初步鉴定为同一物种——U. laetevirens, 在烟台绿潮的物种组成中占绝对优势, 是造成此次局地绿潮的原因种; 而CNS00602 可以初步鉴定为U. australis。

图1 烟台局地绿潮及其致灾物种的形态和分子鉴定Fig.1 A green tide in Yantai coast and the morphological and molecular analysis of the causative species

2.2 烟台绿潮物种U. laetevirens CNS00533 叶绿体基因组特征

相较于单基因分子标记, 集合多基因的叶绿体基因组信息可能提供更高的分辨率, 能够更准确地反映烟台绿潮原因种的分类地位, 进一步确定U.laetevirens、U. rigida以及U. australis的遗传进化关系。为此, 本研究构建并深入解析了烟台绿潮物种的叶绿体基因组。本研究共完成了4 个株系CNS00533、CNS00392、CNS00396 和CNS00398 叶绿体基因组的组装, 4 个叶绿体基因组碱基组成完全一致。以CNS00533 叶绿体基因组(图2)为参考, 比较分析了11 个株系之间的遗传多样性, 仅发现一个单碱基变异位点, 进一步表明从烟台绿潮中分离得到的这11个株系是同一个株系克隆繁殖的结果。后续分析仅以CNS00533 为代表。

图2 烟台绿潮物种U. laetevirens CNS00533 叶绿体基因组Fig.2 The chloroplast genome of U. laetevirens CNS00533 in the local green tide in Yantai

CNS00533 的叶绿体基因组为110 889 bp, 与U.laetevirens(MT179351, 103 444 bp)、U. rigida(MT179352, 96 673 bp)和U. australis(MT179348,99 820 bp)三个叶绿体基因组相比略长(表 1)。CNS00533 叶绿体基因组的AT 含量与U. laetevirens(MT179351)比较相近, 74.4%—74.6%, 比U. australis与U. rigida的AT 含量(分别为74.8%和75.4%)略低。四个叶绿体基因组共有71 个蛋白编码基因、3 个rRNA 基因(rrl、rrs和rrn5)和26 个tRNA 基因。CNS00533 与U. laetevirens(MT179351)的叶绿体基因组的tRNA 基因完全相同, 它们比U. australis多了一个tRNA 基因trnF(aaa)。同时,U. rigida相比其他三个叶绿体基因组多了一个trnR(ucu)。总体而言,CNS00533 与U. laetevirens(MT179351)的叶绿体基因组的基因组成完全一致, 进一步支持上述对CNS00533的物种鉴定。

表1 烟台绿潮物种U. laetevirens CNS00533 与近缘物种U. laetevirens(MT179351)、U. rigida(MT179352)、U.australis(MT179348)的叶绿体基因组特征Tab.1 Characteristics of chloroplast genomes of U. laetevirens CNS00533 (MW531676), U. laetevirens (MT179351), U. rigida(MT179352) and U. australis (MT179348)

2.3 烟台绿潮物种U. laetevirens CNS00533 叶绿体基因组共线性分析

为比较分析CNS00533 与其他3 个近缘物种叶绿体基因组的相似度, 本研究将烟台U. laetevirens株系CNS00533 (MW531676)分别与U. laetevirens(MT179351)、U. rigida(MT179352)、U. australis(MT179348)进行基因组共线性分析(图3)。结果显示,CNS00533 与U. laetevirens(MT179351)的叶绿体基因组基因数目和基因排列顺序完全一致(图 3a)。CNS00533 与U. rigida的叶绿体基因组相比,trnR基因为U. rigida的叶绿体基因组特有(图 3b)。CNS00533 与U. australis的叶绿体基因组存在显著的差异, 二者之间存在一个包括45 个基因的基因组反转(inversion)事件(图3c)。同时, 相比于U. australis的叶绿体基因组, CNS00533 还多一个基因trnF(aaa)。因此, 共线性分析也支持CNS00533 作为U.laetevirens的鉴定。

图3 U. laetevirens CNS00533 ( MW531676)与U. laetevirens (MT179351)(a)、U. rigida (MT179352)(b)、U. australis(MT179348)(c)叶绿体基因组共线性分析Fig.3 Synteny analysis of chloroplast genomes of U. laetevirens CNS00533 MW531676 vs. U. laetevirens MT179351 (a), U.laetevirens CNS00533 MW531676 vs. U. rigida MT179352 (b), and U. laetevirens CNS00533 MW531676 vs. U. australis MT179348 (c)

2.4 烟台绿潮物种U. laetevirens CNS00533 的进化分析

为了探究烟台绿潮物种U. laetevirensCNS00533的遗传进化关系, 本研究针对已发表的19 株石莼绿藻和CNS00533 的叶绿体基因组, 构建了基于蛋白质编码基因氨基酸序列的遗传进化树(图4)。结果显示,石莼属物种聚类为 4 个大的 clade, 分别是Ulva lactuca-ohnoi-laetevirens-rigida-gigantea、Ulva linzaprolifera-flexuosa、Ulva australis-fenestratapseudorotundata和Ulva compressa。烟台绿潮原因种CNS00533 聚类到Ulva lactuca-ohnoi-laetevirensrigida-gigantea, 且与U. laetevirens遗传进化距离更近。同时, 该结果也表明U. laetevirens与Ulva australis遗传进化关系较远, 不是同一物种, 与Hughey 等(2021)的结论相反。

图4 针对U. laetevirens CNS00533 和其他石莼属叶绿体蛋白质编码基因氨基酸序列的遗传进化树Fig.4 The phylogenetic tree based on concatenated amino acid sequences of chloroplast protein-coding genes of U. laetevirens CNS00533 and other Ulva species.

2.5 烟台绿潮物种U. laetevirens CNS00533 与近缘物种大片段差异序列和内含子组成分析

上述关于分子标记(tufA)、基因组成、共线性分析和遗传进化关系等多方面的分析表明CNS00533 是U. laetevirens的一个株系。比较分析CNS00533 和U. laetevirens(MT179351)的叶绿体基因组DNA 序列发现, 相较于U. laetevirensMT179351, CNS00533 叶绿体基因组(MW531676)多出6 个较长的插入序列(表2, 序号1—6), 长度范围1—2 kb。这6 个插入序列分别对应6 个完整的内含子(表3, 5 个I 型和1 个II 型内含子)。比较分析CNS00533 和U. rigida(MT179352)的叶绿体基因组发现, CNS00533 叶绿体基因组多了10 个插入序列(表2)。其中1—6 号插入序列与上面的6 个插入序列相同, 其他4 个插入序列分别是基因组atpB、rrl基因内的1 个II 型和1 个I 型内含子。为了验证不同物种间非插入序列的相似性, 本研究分析发现 CNS00533 MW531676 与U. laetevirensMT179351 插入序列之外的序列相似性达到99.8%。与此相比, CNS00533 MW531676 与U. rigidaMT179352 的DNA 序列相似性只有93.1%。

表2 烟台绿潮物种U. laetevirens CNS00533 叶绿体基因组中相比于U. laetevirens MT179351 和U. rigida MT179352 的大片段插入序列统计Tab.2 The long sequences in the chloroplast genome of U. laetevirens CNS00533 but not in U. laetevirens MT179351 or U. rigida MT179352

本研究进一步统计分析了U. l a e t e v i re n sCNS00533、U. laetevirens(MT179351)和U. rigida(MT179352)叶绿体基因组内含子的分布情况(表3)。结果显示,U. rigida(MT179352)只有1 个I 型内含子,位于其rrl基因内;U. laetevirens(MT179351)有1 个位于atpB基因内的II 型内含子和2 个位于rrl基因内的I 型内含子; 相对而言, CNS00533 含有7 个I 型内含子(位于psbD、psbC、psbB和rrl基因内)和2 个II型内含子(位于atpB基因内), 远远多于前两者。

表3 烟台绿潮物种U. laetevirens CNS00533、U. laetevirens (MT179351)和U. rigida (MT179352)叶绿体基因组内含子分布情况Tab.3 The intron contents in chloroplast genomes of U. laetevirens CNS00533, U. laetevirens (MT179351) and U. rigida (MT179352)

综上, CNS00533 与U. rigida(MT179352)叶绿体基因组的差异体现在内含子组成和基因编码区, 差异较大; 而CNS00533 与U. laetevirens(MT179351)叶绿体基因组基因编码区差异较小, 内含子组成是其种内差异的主要方面。

3 讨论

3.1 U. laetevirens 的形态和分子鉴定

根据文献报道,U. laetevirens具有掌状叶片、浅裂、不规则分叉、不规则的边缘突起和二层细胞基质等等, 这些形态特征与其他石莼属叶状藻体有不同程度的相似性, 特别是与其近缘的U. rigida极易混淆(Steinhagenet al, 2019b)。Sfriso 发现U. rigida的根状体细胞相比U. laetevirens更厚、更暗(Sfriso, 2010)。Sfriso (2010)的研究还发现U. laetevirens的假根和基部边缘没有突起。此外, 石莼属物种具有形态可塑性,不同的营养环境、光照以及共附生菌均会影响其形态(Tayloret al, 2001; Fuet al, 2008)。总体而言, 仅仅基于形态学的描述, 难以实现对包括U. laetevirens在内的石莼属物种进行准确的鉴定。

分子学方法已成为物种分类的一个有效手段,可以在形态学鉴定的基础上, 更加准确地鉴定石莼属物种。但是, 分子学方法的准确鉴定仍然依赖于分子标记的选择。在已报道的研究中,U. laetevirens、U. rigida和U. australis分类存在争议。Batista 等(2018)基于分子标记 ITS 和rbcL的分析结果认为U.laetevirens和U. rigida区分不开。Steinhagen 等(2019b)甚至基于此在其研究中将U. laetevirens和U. rigida作为同一物种。Hughey 等(2021)根据分子标记rbcL分析结果确定U. laetevirens与U. australis为同一物种。Lee 等(2019)在分析韩国济州岛的样品时, 尽管采用tufA和rbcL联合建树分析, 但由于其参考序列中没有U. rigida相关序列作为参考, 所以分析结果并没有说明其建树方式能有效区分U. laetevirens和U. rigida。Fort 等(2021)的研究中对欧洲海域采集的110 株叶状石莼进行单核苷酸多态性(single nucleotide polymorphism, SNP)分析, 结果明确聚类为6 个clade, 基于分子标记ITS 和rbcL的结果, 6 个SNP-clade分别对应6 个物种, 其中包括U. laetevirens和U. rigida, 其后续研究均视U. laetevirens和U.rigida为两个物种, 且细胞器基因组均显示二者种间差异极大, 种内差异极小。

此外, Fort 等(2021)比较了tufA、rbcL的单基因进化树以及线粒体、叶绿体的全部蛋白编码序列的进化树, 结果显示这四种方式的分类效果基本一致,tufA的bootstrap 支持率相对高于rbcL, 认为单分子标记tufA足以区分叶状石莼属物种。本研究为了鉴定烟台株系CNS00533, 对其tufA进行分析, 结果(图1d)显示CNS00533 与U. laetevirens和U. rigida的参考序列聚类在一个clade 中, NCBI 比对发现仅与U.laetevirens(MT179351)(Fortet al, 2021)的tufA序列(1 224 bp) 100%一致。另外, 根据石莼属属内物种的叶绿体基因组蛋白编码基因的进化树(图4), CNS00533与U. laetevirens有更近的遗传进化关系。此外, 为了排除叶绿体基因组中可能的插入序列的影响, 本研究剔除了CNS00533 与U. laetevirens(MT179351)、U. rigida(MT179352)参考叶绿体基因组彼此间的长片段插入序列后发现, CNS00533 与U. laetevirens相似度达到99.8%, 与U. rigida相似度只有93.1%。Fort等(2021)在进行6 个石莼属物种种间、种内叶绿体基因组多样性分析时, 发现种间蛋白编码区单核苷酸变异(single nucleotide variant, SNV)位点至少391 个,而U. laetevirens种内只有36 个SNV 位点。本研究比较分析了CNS00533 与爱尔兰U. laetevirens的叶绿体基因组全序列, 发现仅有51 个SNV 位点, 进一步说明二者为种内差异。以上分子标记、叶绿体蛋白编码基因遗传进化分析和全序列比对的结果, 均表明烟台株系CNS00533 为U. laetevirens。

3.2 U. laetevirens 烟台株系与爱尔兰株系间的异同

同一物种的不同株系之间可能存在遗传变异,包括大尺度的遗传变异。截至目前已公布的记录中,仅有爱尔兰海域和中国烟台海域的U. laetevirens株系的叶绿体基因组[MT179351 (Fortet al, 2021)和MW531676]被构建和解析。经比较发现, 二者的叶绿体基因组的基因种类、数量和排列顺序完全一致, 包括71 个蛋白编码基因、27 个tRNA 基因和3 个rRNA基因基因(图2, 3)。Fort 等(2021)的研究结果显示欧洲海域U. laetevirens种内差异极低, 主要体现在种内仅0.6 SNVs/kb 编码序列。本研究的结果显示烟台株系与爱尔兰株系之间SNV 也极低。同时, 烟台株系的叶绿体基因组比爱尔兰株系长7 445 bp, 这个长度差异在于二者内含子的组成差异, 前者比后者多了6 个完整的内含子序列(表2, 3)。种内内含子差异显著的现象并不例外, 本研究在整理NCBI 数据库的所有石莼属物种内含子分布情况时发现,U. compressa不同株系的叶绿体基因组中内含子相差较大, MK069584[U. mutabilis更名为U. compressa(Steinhagenet al,2019a)]含有10 个内含子, 而KX595275 只有4 个内含子。 另外, 在已报道的两株Metschnikowia santaceciliae线粒体中, 二者之间仅cob基因含有的内含子个数就相差6 个(Leeet al, 2020)。对于U.laetevirens烟台株系和爱尔兰株系间的内含子差异,可能是不同地理环境压力下自然选择的结果。

3.3 U. laetevirens 是我国新记录石莼属(绿潮)物种

截至2021 年1 月, 根据已有记载信息, 我国海域共记录29 个石莼属物种(丁兰平等, 2015; Xieet al,2020), 其中11 个物种已在全球有引发绿潮的报道,包括U. australis、U. chaugulli、Ulva clathrata、U.compress、U. flexuosa、U. intestinalis、U. linza、U.meridionalis、U. prolifera、Ulva ridiga、U. tepida(Guidoneet al, 2013; Smetaceket al, 2013; Gaoet al,2018; Golubkovet al, 2018; Xieet al, 2020)。本研究报道的烟台U. laetevirens是我国绿潮事件中该物种的首次记录。该物种于1854 年在澳大利亚维多利亚的菲利普港首次采集并记录, 现已在北美长岛湾、加拿大、爱尔兰、法国、荷兰、葡萄牙、突尼斯和韩国均有采集记录(Kirkendaleet al, 2013; Maoet al, 2014;Miladiet al, 2018; Leeet al, 2019; Fortet al, 2021)。此外, 本研究发现在 NCBI 数据库中U. laetevirensCNS00533 的ITS 序列与我国研究者Miao 等提交的8个标记“Ulvasp.”的ITS 序列(KY446828 等, 未发表)一致, 说明该物种在我国其他地区也有分布。U.laetevirens在澳洲、美洲、欧洲和亚洲海域均已有绿潮暴发记录, 该物种快速增殖和分布广泛的特征需要我们持续关注和监管。

4 结论

本研究聚焦了2020 年11 月初在烟台第一海水浴场发生的绿潮事件, 通过比较叶绿体基因组学分析证明该石莼属绿藻为U. laetevirens, 该物种在我国海域尚属首次记录。本研究构建了烟台株系U.laetevirensCNS00533 的叶绿体基因组(MW531676),相比于爱尔兰株系U. laetevirens(MT179351)二者之间蛋白编码基因、tRNA 基因和rRNA 基因的种类和数目一致, 仅有51 个SNV 位点, 序列组成差异在于烟台株系的叶绿体基因组的4 个基因psbD、psbC、psbB和atpB中共多了6 个内含子。U. laetevirens烟台株系与爱尔兰株系可能存在其他潜在的差异, 包括核基因组和线粒体基因组的组成差异等, 这些需要进一步解析以期为理解我国U. laetevirens绿潮暴发机制提供分子理论支持。对我们海域石莼属绿潮物种(包括新记录种)的全面调查和生物地理分布研究,将是对我国石莼绿潮实施有效监管和可靠预测的必要工作之一。