SPECIES IDENTIFICATION AND POPULATION GENETIC DIVERSITY OF NON-NATIVE APPLE SNAILS (AMPULLARIIDEA: POMACEA) IN THE WEST LAKE
-
摘要: 为了解杭州西湖景区外来生物福寿螺的分布现状和群体遗传变异特征, 研究详细调研了西湖景区福寿螺的分布现状, 并针对连续六年采集的福寿螺样本开展了分子鉴定及种群遗传多样性研究。调查结果表明, 西湖福寿螺集中发生在苏堤以西, 其中西北部的西里湖、茅家埠、杭州花圃及岳湖水域的发生密度高, 在西南部的乌龟潭、浴鹄湾等水域发生密度较低。在602 bp长的线粒体COⅠ序列中, 共检测到79个变异位点, 生成了6个单倍型。序列相似性和系统发育关系分析表明, 西湖水域中Hap1-Hap5为小管福寿螺, Hap6为斑点福寿螺; 单倍型多样性、核苷酸多样性和核苷酸平均差异数分析表明, 不同年份的小管福寿螺单倍型多样性Hd为0.229—0.600, 核苷酸多样性π为0.011—0.027, 核苷酸平均差异数k为6.797—16.000, 两种福寿螺均有较低的种群遗传多样性, 显著性分析表明2018—2020年与2015—2017年间的种群遗传多样性无显著性差异。此外, 结合核基因EF1α序列分析表明西湖水域的福寿螺种群包括CC型(23.08%)、CM型(23.08%)和CB型(53.84%), 杂种所占比例高达76.92%以上, 群体广泛发生了种间杂交渐渗。由于引配水工程的连通性, 西湖水域的福寿螺向周边水系扩散风险较大, 研究为西湖景区福寿螺的防控提出了建议。Abstract: The exotic apple snails are big freshwater snails that have spread and damaged the West Lake for nearly 10 years. To understand the species and population genetic diversity of apple snails, we explored the distribution of the apple snails in the West Lake. Based on the samples collected in the past six years, identified the species, and analyzed their population genetic diversity. The results showed that the distribution of apple snails was concentrated in west of the Su Causeway, with a high density in the northwest area including Xilihu, Maojiabu, Hangzhou Flowerbed, and Yuehu, and low density in the southwest aera such as Wuguitan and Yuhuwan. We detected 79 variable sites, and defined 6 haplotypes based on mitochondrial COⅠ sequences. Sequences similarity and phylogenetic relationships analysis revealed that the Hap1-Hap5 were Pomacea canaliculata and the Hap6 was P. maculata in the West Lake. The analysis of haplotype diversity, nucleotide diversity, and average number of nucleotide difference indicated that the haplotype diversity (Hd) of P. canaliculata was 0.229—0.600, the nucleotide diversity (π) was 0.011—0.027, the average number of nucleotide difference (k) was 6.797—16.000, and both P. canaliculata and P. maculata in the West Lake have low genetic diversity. There was no significant difference in genetic diversity between 2018—2020 and 2015—2017. In addition, nuclear EF1α sequence analysis demonstrated that the apple snail populations of P. canaliculata in the West Lake included CC (23.08%), CM (23.08%), and CB (53.84%) types. Introgressive hybridization were widely observed in the apple snail populations in the West Lake with a rate over 76.92%. Due to the connectivity of the water diversion project, there is a risk of apple snails spreading to the surrounding wetlands in the West Lake. This study provides some suggestions for the establishment of effective prevention and control strategies of apple snails in the West Lake.
-
福寿螺隶属于腹足纲Gastropoda新进腹足目Caenogastropoda瓶螺科Ampullariidae福寿螺属Pomacea, 是原产于南美洲、中美洲和加勒比地区的一类大型淡水螺[1]。当前, 福寿螺入侵亚洲、欧洲、北美洲和大洋洲的许多国家和地区, 被列入世界100种危害最大的入侵生物名单[2], 也是我国首批公布的外来入侵有害生物之一。福寿螺在我国最早于1981作为食用螺引入到台湾地区, 随后引入广东省中山市, 并在之后掀起的养殖热中引至全国各地[3—5]。由于福寿螺口味差, 未能产生良好的市场效益而被大量弃养或逃逸到农田、湿地等自然生境。因具有繁殖力强、成熟周期短、适应性和耐受力强等特点, 加上缺乏本土天敌控制, 福寿螺种群迅速在我国南方地区扩散及暴发, 造成巨大经济损失[4, 6]。它们不仅严重为害水稻幼苗及多种水生蔬菜[7], 并且破坏水生生态系统平衡[8, 9]。此外, 福寿螺作为广州管圆线虫等多种人畜共患寄生虫的中间宿主, 严重威胁着人类身体健康[10]。
由于福寿螺种间形态近似, 种内形态特征多变, 基于形态特征鉴别种类十分困难[11]。最新研究表明, 福寿螺在入侵地广泛发生了种间渐渗杂交, 导致卵形态特征趋同, 基于形态特征鉴别种类难度进一步增大[12]。近年来, 国内外研究学者基于DNA条形码技术极大程度地解决了福寿螺种类鉴定的难题, 并基于形态与分子特征结合方法构建了福寿螺综合分类系统。Rawlings等[13]基于线粒体12S—16S rDNA及线粒体COⅠ基因序列分析, 从小管福寿螺P. canaliculata中鉴别出近缘种孤岛福寿螺P. insularum; 进一步地, Hayes等[11]通过形态学结合分子系统学研究将其归并为斑点福寿螺P. maculata的同物异名。Lv等[14]和Yang等[15]基于DNA条形码分析从我国福寿螺样品中鉴别出一个隐种; Yang和Yu等[16]利用系统的螺壳、软体解剖学、卵及幼螺形态特征描述, 将这一隐种命名为隐秘福寿螺P. occulta。当前, 小管福寿螺、斑点福寿螺和隐秘福寿螺是在我国湿地生态系统中严重危害的3个种。
近年来, 种间杂交渐渗的新发现成为了福寿螺系统进化研究的焦点。在原产地及入侵地均检测到小管福寿螺和斑点福寿螺的杂交渐渗种群[17, 18]。杂种的生物学特性(如耐低温水平)[19]、形态特征(如卵形态)[12]呈现小管福寿螺和斑点福寿螺的中间型。一般来说, 杂种可以理解为不同物种间交配形成的新生命体, 也可以更广泛地理解为形态可区分或遗传间断的群体间的个体交配产生的后代, 这个后代与父母本物种或者群体在形态或者遗传上是连续的, 并且没有形成自己独特的遗传特征[20, 21]。
当前, 福寿螺在西湖、西溪国家湿地公园、洱海等景区水域大量发生, 景观及湿地的生物多样性受到了严重的威胁。由于福寿螺在水面上方产卵, 卵块呈鲜红色, 在春夏繁殖高峰期, 大量卵块附着在水生植物、石壁及建筑物上, 与景区绿色的主体背景形成强烈对比, 引起游客不适, 引发社会的广泛关注[22]。然而有关福寿螺在景区的种类、分布特点及种群遗传多样性尚未有系统报告。本研究以杭州西湖景区为例, 在系统调研福寿螺分布情况的基础上, 结合连续6年采集的样本开展种类鉴定及种群遗传多样性分析, 以期为景区螺害的综合治理提供理论指导。
1. 材料与方法
1.1 西湖水域福寿螺样品采集与发生调研
2017年至2020年春夏季连续开展了杭州西湖水域福寿螺样品的采集, 包括成螺及卵块共计82个样品(表 1和图 1)。于2020年7月在福寿螺发生盛期, 对杭州西湖水域及杭州植物园和太子湾公园水域(30°13′51″—30°15′58″N, 120°07′34″—120°10′09″E)开展福寿螺发生调研。在西湖水域沿岸每100 m为一个调查点, 在杭州植物园和太子湾公园的水域沿岸每50 m为一个调查点。每个调查点仔细观察是否有福寿螺卵块及成螺, 记录每个福寿螺栖息地景观植物种类。
表 1 西湖水域福寿螺样品、序列及单倍型分布信息表Table 1. Samples, sequences, and haplotype distribution of the apple snails from the West Lake日期
Date数量(成螺/卵块)
No. (Adult/Egg clutch)GenBank序列号
Accession No.Hap1 Hap2 Hap3 Hap4 Hap5 Hap6 参考文献
Reference2020.05 5/14 MW051334-353 14 4 0 0 1 0 本研究This study 2019.05 23/17 MW051295-334 35 4 1 0 0 0 本研究This study 2018.05 13/0 MW051282-294 8 4 0 1 0 0 本研究This study 2017.07 0/10 MW051272-281 1 3 6 0 0 0 本研究This study 2016.11 26/2 MK189337-62, MH211409-10 17 9 0 0 0 2 [24] 2015.06 0/10 KT852718-20, 31-36, 61 4 6 0 0 0 0 [25] 总计Total 67/53 / 79 30 7 1 1 2 / ![]()
图 1 西湖水域的福寿螺形态图A—C. 雌螺; D. 雄螺, 均为小管福寿螺Figure 1. Morphology of apple snails from the West LakeA—C. female; D. male, four snails are all P. canaliculata1.2 种类分子鉴定
DNA提取、PCR及测序 从福寿螺成螺中取约10 mg腹足、卵块中取单粒卵, 使用血液/细胞/组织基因组DNA提取试剂盒(天根生化科技有限公司)提取样品基因组DNA, 具体步骤参考试剂盒说明书。每个样品获取的DNA溶液均溶解于200 μL的ddH2O中, 使用超微量紫外分光光度仪(NanoDrop, ND-1000)检测浓度及纯度后, 置于–20℃保存。
使用通用引物对LCO1490(5′-GGTCAACAAATCATAAAGATATTGG-3′)和HCO2198(5′-TAAACTTCAGGGTGACCAAAAAATCA-3′)[23]扩增线粒体COⅠ基因的DNA条形码区段。PCR扩增在25 µL的体系反应中进行, 包括为12.5 µL的TaKaRa Premix Taq, 9.5 µL的ddH2O, 1 µL的基因组DNA, 上下游引物(10 µmol/L)各1 µL。PCR循环条件为94℃预变性3min; 94℃变性30s, 50℃退火30s, 72℃延伸1min, 共34个循环; 72℃延伸10min, 4℃终止反应(PCR扩增体系及条件参考Yang等[5])。
使用1%的琼脂糖凝胶检测PCR扩增产物, 将有单一明亮条带的扩增产物送由生工生物工程(上海)股份有限公司纯化并测序。使用Geneious 11.1.5软件结合人工方式对获取的原始序列中的模糊碱基进行校正并剪切引物。
线粒体COⅠ数据集及单倍型分析 从GenBank中下载分别于2015年和2016年采集自杭州西湖的福寿螺样品的线粒体COⅠ序列共计38条(表 1)。结合本研究测定的序列, 使用Geneious 11.1.5软件进行全局比对并剪切齐末端后, 形成共计120条序列的602 bp的线粒体COⅠ数据集。使用DnaSP 5软件计算线粒体COⅠ数据集的单倍型, 并使用Geneious 11.1.5软件分析单倍型的多态位点数。
遗传距离及系统发育关系分析 由于已有研究表明, 我国野外种群的福寿螺主要包括小管福寿螺、斑点福寿螺和隐秘福寿螺3个种, 因而从GenBank下载了3种福寿螺的线粒体COⅠ序列各3条用于比较分析, 包括阿根廷种群的小管福寿螺(P. canaliculata_AR, GenBank序列号: EF514961), 巴西种群的斑点福寿螺(P. maculata_BR, EF514999), 日本种群的小管福寿螺(P. canaliculata_JP, AB433765)和斑点福寿螺(P. maculata_JP, AB433777), 中国种群的小管福寿螺(P. canaliculata_CN, KP310268)、斑点福寿螺(P. maculata_CN, KP310490)和隐秘福寿螺(P. occulta_CN1-3, KT852790, KT852789, KP310451)。
利用Geneious 11.1.5软件基于Kimura 2-parameter(K2P)遗传距离模型计算西湖种群福寿螺各单倍型序列及其与参考序列间的相似度。以神秘福寿螺(EF515067)为外群, 分别构建邻接法(Neighbor-joining, NJ)、最大似然法(Maximum likelihood, ML)和贝叶斯法(Bayesian inference, BI)系统发育树。其中, 分别使用MEGA X软件基于K2P模型构建NJ树, 置信值设置为1000次; 利用IQ-Tree 1.6.12构建ML树, 置信值设置为1000次, 其中由内置的ModelFinder软件计算的最佳模型为HKY+F+I; 使用MrBayes 3.2.5软件构建BI树, 运行4个马尔可夫(MCMC)链5000000代, 每1000代储存一次树, 利用burnin参数舍弃前25%后, 计算合一树。
种群遗传多样性分析 利用DnaSP 5软件计算种群遗传多样性参数, 包括单倍型多样性(Haplotype diversity, Hd)、核苷酸多样性(Nucleotide diversity, π)和核苷酸平均差异数(Average number of nucleotide difference, k), 并进行Tajima’s D 中性检验。使用SPSS软件对2018—2020年和2015—2017年的遗传多样性参数进行显著性分析。
1.3 杂交渐渗分析
已有报道表明, 我国分布福寿螺广泛发生种间遗传渐渗, 且同一卵块的不同个体因多父性或杂交分离而产生不同渐渗型[12]。因而, 本研究针对2015—2020年采集的65头小管福寿螺成螺样品, 开展基于核基因EF1α分型的种间杂交渐渗分析, 具体参考Yang等[12]建立的特异引物多重PCR方法。
2. 结果
2.1 西湖水域福寿螺发生情况
在西湖水域福寿螺的发生调研中, 通过目测福寿螺卵块的发生共43个调查点发现了福寿螺, 成螺仅在曲院风荷景区的一个样点发现。成螺个体大, 螺高5.36—7.60 cm, 螺宽4.04—5.74 cm, 螺壳较厚, 壳面光滑、呈绿色或橄榄绿色, 缝合线深、凹入, 螺旋尖的高度低, 壳口内侧颜色红褐色(图 1)。调研结果表明, 福寿螺仅分布于苏堤以西的水域, 其中西里湖西岸、茅家埠北岸、杭州花圃景区及岳湖的曲院风荷景区发生密度较高, 在茅家埠南岸、西里湖南部及乌龟潭发生密度低。在苏堤以东的西湖外湖、北里湖和三潭印月景区均未发现有福寿螺分布。北部的杭州植物园、南部的小南湖和太子湾公园的水域沿岸未发现有福寿螺分布(图 2)。
![]()
图 2 西湖水域福寿螺的分布情况坐标符号表示福寿螺的发生地点; 黑色箭头为进水口, 1—6分别为莲花峰、小南湖、浴鹄湾、乌龟潭、茅家埠1、茅家埠2进水口; 灰色箭头为出水口, 7—15分别为岳湖、北里湖、圣塘闸、五公园、一公园、大华饭店、涌金闸、涌金池、柳浪出水口Figure 2. Distribution of the apple snails in the West LakeThe coordinates indicate the location of the apple snails; the black arrow indicate the inlet of the West Lake, 1—6. Lianhuafeng, Xiaonanhu, Yuhuwan, Wuguitan, Maojiabu 1, Maojiabu 2 the grey arrow indicate the outlet of the West Lake, 7—15. Yuehu, Beilihu, Shengtangzha, Wugongyuan, Yigongyuan, Dahuafandian, Yongjinzha, Yongjinchi, Liulang西湖水域福寿螺的分布地点均有大型水生及半水生植物覆盖, 没有植物栽培的水域未观察到福寿螺成螺及卵块。其中, 卵块附着的大型挺水植物主要包括: 荷花Nelumbo nucifera(在19个栽培荷花的调查点中的14个点发现福寿螺卵, 14/19)、再力花Thalia dealbata(9/15)、芋Colocasia esculenta(4/16)、黄菖蒲Iris pseudacorus(2/8)、香蒲Typha orientalis(1/4)、萍蓬草Nuphar pumilum(1/4)、水葱Schoenoplectus tabernaemontani(1/5)和菖蒲Acorus calamus(1/5)等。其中, 在荷花和再力花种植区发现福寿螺卵块的比例较高。
2.2 西湖福寿螺种类鉴定
本研究测序获得了于2017—2020年采集自西湖水域的福寿螺线粒体COⅠ序列82条, 序列长度为602 bp, 均已提交至GenBank(序列号MW051272-MW051353)。采集自2015—2020年西湖水域福寿螺的120条线粒体COⅠ序列集共生成6个单倍型, 其中Hap1-Hap6分别包含79条、30条、7条、1条、1条和2条序列(表 1)。在6个单倍型序列中共发现了79个多态位点, 占分析位点总数的13.12%(图 3)。
![]()
图 3 西湖福寿螺线粒体COⅠ基因单倍型变异位点Figure 3. Variation sites of mitochondrial COⅠ gene haplotype in the apple snails from the West Lake根据序列相似度, 西湖水域的福寿螺线粒体COⅠ单倍型与GenBank序列可分为3个组。其中, Hap1-5与不同地理种群的小管福寿螺为一组, 相似度≥94.63%, 而与斑点福寿螺和隐秘福寿螺相似度≤89.60%, 其中Hap1和Hap4-Hap5与P. canaliculata_AR的相似度≥99.16%, Hap2-3与P. canaliculata_JP和P. canaliculata_CN的相似度≥97.82%; Hap6与不同地理种群的斑点福寿螺为一组, 序列相似度≥98.32%, 与小管福寿螺序列相似度≤89.93%, 与隐秘福寿螺相似度≤94.13%; 隐秘福寿螺序列间为一组(表 2)。
表 2 西湖福寿螺单倍型序列与GenBank下载序列的相似性及碱基差异数Table 2. Similarity and number of variation sites among the haplotypes and sequences retrieved from GenBank1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 P. canaliculata_AR 0 1 4 29 30 30 30 56 55 60 60 57 56 56 2 Hap1 100 1 4 29 30 30 30 56 55 60 60 57 56 56 3 Hap4 99.83 99.83 5 30 31 31 31 57 56 61 61 58 57 57 4 Hap5 99.33 99.33 99.16 31 32 32 32 56 56 60 60 58 57 57 5 P. canaliculata_JP 95.13 95.13 94.97 94.80 12 12 1 60 57 60 60 55 54 54 6 P. canaliculata_CN 94.97 94.97 94.80 94.63 97.99 0 13 62 59 58 58 56 55 55 7 Hap2 94.97 94.97 94.80 94.63 97.99 100 13 62 59 58 58 56 55 55 8 Hap3 94.97 94.97 94.80 94.63 99.83 97.82 97.82 59 56 59 59 56 55 55 9 P. maculata_BR 90.60 90.60 90.44 90.60 89.93 89.60 89.60 90.10 9 10 10 41 40 40 10 P. maculata_JP 90.77 90.77 90.60 90.60 90.44 90.10 90.10 90.60 98.49 9 9 36 35 35 11 P. maculata_CN 89.93 89.93 89.77 89.93 89.93 90.27 90.27 90.10 98.32 98.49 0 35 34 34 12 Hap6 89.93 89.93 89.77 89.93 89.93 90.27 90.27 90.10 98.32 98.49 100 35 34 34 13 P. occulta_CN1 90.44 90.44 90.27 90.27 90.77 90.60 90.60 90.60 93.12 93.96 94.13 94.13 1 3 14 P. occulta_CN2 90.60 90.60 90.44 90.44 90.94 90.77 90.77 90.77 93.29 94.13 94.30 94.30 99.83 2 15 P. occulta_CN3 90.60 90.60 90.44 90.44 90.94 90.77 90.77 90.77 93.29 94.13 94.30 94.30 99.50 99.66 注: 对角线以下为序列相似性, 对角线以上为碱基差异数。AR、BR、CN、JP分别表示阿根廷、巴西、中国和日本地理种群Note: Data below the diagonal are similarity and those above the diagonal are variable sites. AR, BR, CN, and JP indicate the populations of Argentina, Brazil, China, and Japan 尽管置信值不同, NJ、ML和BI系统发育树呈现一致的系统分枝结构。其中, Hap1-5与小管福寿螺聚类为一枝, Hap6与斑点福寿螺聚类为一枝, 且具有较高的置信值(图 4)。在小管福寿螺分枝内, Hap1、Hap4-5与P. canaliculata_AR聚为亚枝, 表明具有较近的系统发育关系, 而Hap2-3分别与P. canaliculata_CN和P. canaliculata_JP具有较近的系统发育关系(图 4)。
![]()
图 4 西湖水域福寿螺的单倍型及GenBank序列构建的系统发育树系统发育分支节点依次为NJ/ML/BI置信值, 仅显示>50%的数值Figure 4. The phylogenetic trees reconstructed by the haplotypes of apple snails from the West Lake and the GenBank sequencesThe branch nodes of phylogenetic trees denote NJ/ML/BI bootstrap supports, and only>50% values are displayed相似度分组与NJ系统发育分枝结果一致表明, Hap1—5与小管福寿螺的基因序列的一致性达94.63%—100%, Hap6与斑点福寿螺的基因序列的一致性达98.32%—100%。综合系统发育关系和遗传距离分析结果表明, 本研究中120个样本包括小管福寿螺118个和斑点福寿螺2个。
2.3 西湖福寿螺种群遗传多样性
西湖水域仅检测到一个斑点福寿螺单倍型, 单倍型多样性低; 小管福寿螺共形成5个单倍型。不同年份的小管福寿螺单倍型多样性Hd为0.229—0.600, 核苷酸多样性π为0.011—0.027, 核苷酸平均差异数k为6.797—16.000; 总体单倍型多样性Hd为0.488, 核苷酸多样性π为0.023, 核苷酸平均差异数k为13.494。统计分析表明, 2018—2020年和2015—2017年间的小管福寿螺种群遗传多样性无显著性差异。中性检验结果显示, 2017年和2019年小管福寿螺种群的Tajima’s D为负值, 其他年份均为正值, 总Tajima’s D为2.334, 表现出显著偏离中性(P<0.05); 2015年和2016年种群的Tajima’s D分别为2.443和2.956, 达到极显著性检验标准(P<0.01), 其他年份的Tajima’s D均无显著性(表 3)。
表 3 小管福寿螺种群内的遗传多样性分析Table 3. Population genetic diversity of P. canaliculata from the West Lake日期
Date单倍型
多样性
Haplotype diversity (Hd)核苷酸
多样性
Nucleotide diversity (π)核苷酸平均
差异数
Average number
of nucleotide difference (k)Tajima’s D中性检验
Tajima’s D neutrality test2020.05 0.433 0.018 10.901 0.616 2019.05 0.229 0.011 6.797 –0.764 2018.05 0.564 0.023 14.000 1.764 2017.07 0.600 0.019 11.200 –0.696 2016.11 0.471 0.024 14.123 2.956** 2015.06 0.533 0.027 16.000 2.443** 总计Total 0.488 0.023 13.494 2.334* 注: *表示显著偏离中性(P<0.05), **表示P<0.01Note: *indicate significantly deviated from neutrality (P<0.05), **represent P<0.01 2.4 杂交渐渗分析
基于特异引物多重PCR的方法针对核基因EF1α进一步检测到3种不同的杂交渐渗型, 包括单一125 bp扩增产物小管福寿螺型(C型), 单一151 bp扩增产物斑点福寿螺型(M型), 及同时生成125和151 bp两条扩增产物小管福寿螺-斑点福寿螺杂合型(B型; 图 5)。经线粒体COⅠ序列分析, 2015—2020年采集自西湖水域的65头小管福寿螺(C型), 结合核基因型可分为15头CC型(23.08%)、15头CM型(23.08%)和35头CB型(53.84%)。
![]()
图 5 基于EF1α基因多重PCR扩增的西湖福寿螺杂交渐渗型检测M为50 bp DNA ladder, 1—5为西湖水域采集的小管福寿螺样品, 分别为C, M, B, B, M型Figure 5. Introgressive hybridization of the apple snails from the West Lake based on EF1α gene multiplex PCR methodM. 50 bp DNA ladder, 1—5: samples of P. canaliculata collected in the West Lake, represent C, M, B, M type3. 讨论
福寿螺自19世纪80年代引入我国以来, 已在我国南方各省(自治区、直辖市)广泛扩散, 对湿地系统的生物多样性和功能产生重要影响, 例如危害水生植物、捕食本土淡水生物、改变蜉蝣生物组分变化等[25]。本研究调研了杭州西湖水域福寿螺的发生情况, 以苏堤为界, 福寿螺集中分布在苏堤以西的水域, 西北部发生较密集, 南部堤岸及东南部水域分散分布; 其他水域均未发现福寿螺。调研结果表明, 西湖水域的福寿螺分布范围总体得到了较好的控制。此外, 西湖水域的福寿螺主要分布在有大型水生植物生长的浅水区。水生及半水生植物可为福寿螺在新区定殖提供必要条件, 并对扩散有促进作用[26]。福寿螺取食范围广, 对不同的植物的取食程度差别较大。叶建人等[26]根据取食率与生长率的关系将福寿螺寄主植物分为多食快长型、少食快长型、少食慢长型与多食慢长型。在非农业湿地中, 福寿螺发生与危害程度受植物组成影响较大[27, 28]。例如, 福寿螺尤其喜欢营养成分高且干物质含量低的植物, 对苋Amaranthus gangeticus、节节草Commelina diffusa、蕹菜Ipomoea aquatica和浮萍Lemna minor等取食量大; 相反, 对干物质含量高的植物, 如许多常见大型挺水植物包括芦苇Phragmites australis、荷花、芋和萍蓬草等取食量较少[28]。西湖水域分布的大型挺水植物多为福寿螺不喜食或取食后生长速率缓慢的植物, 这可能是限制种群发展的重要因素之一。然而, 在苏堤以东有多处大面积景观植物种植区, 尚未有福寿螺发生, 具有较高的定殖风险, 需加强防控。
福寿螺属是所在瓶螺科种类鉴定最困难的类群[1]。基于DNA条形码技术的分子鉴定结果显示, 西湖水域分布的福寿螺包括小管福寿螺和斑点福寿螺2种, 其中以小管福寿螺为优势种群。虽然隐秘福寿螺已在我国广泛扩散, 在浙江省亦有发生报道[15, 16, 29], 尚未在西湖水域检测到该物种。已有研究表明, 斑点福寿螺在我国集中分布在四川盆地和浙江省, 在浙江仅局限于江干区的局部水域[15], 本研究在西湖水域发现斑点福寿螺进一步扩大了我们对其分布范围的认识。当前已报道的基于线粒体COⅠ序列分析, 共在全国检测到24个小管福寿螺单倍型和3个斑点福寿螺单倍型[5, 14], 我们在西湖分别检测到的两种福寿螺的其中5个和1个单倍型, 具有较低的种群遗传多样性。不同年份群体间的Taijima’s D中性检验结果表明近几年的福寿螺种群总体上达到一个动态平衡, 未经历明显的种群扩张[30]。已有研究表明, 因同一物种不同群体之间存在遗传差异, 引入外来种群可能与已定殖种群进行杂交, 群体之间发生基因交流, 进而影响已定殖种群的遗传多样性和遗传结构[31]。福寿螺已在杭州西湖景区入侵扩散十余年, 通过对近3年(2018—2020)与2015—2017这3年的福寿螺遗传多样性指标进行显著性分析, 包括单倍型多样性Hd、核苷酸多样性π和核苷酸平均差异数k, 发现西湖福寿螺的遗传多样性无显著差异, 表明近6年西湖水域的福寿螺的遗传多样性比较稳定, 推测没有新引入种群。
基因渐渗是提升种群遗传多样性的重要机制之一[12]。核基因EF1α序列中包含快速进化的内含子, 能够为近期进化提供系统发育信息[32, 33], 被用于检测福寿螺种间杂交渐渗的分子标记[17]。Yang等[12]基于核基因EF1α序列分析, 发现所检测的单个小管福寿螺的卵块中均能检测到C型、M型和B型3种基因型, 可能与多父性生殖策略有关, 同时幼螺和卵中杂种的比例显著高于成螺, 推测可能与不同杂交渐渗型的卵孵化率和幼螺存活率不同有关。本研究中基于核基因EF1α和线粒体COⅠ基因的种间杂交渐渗检测结果表明, 西湖水域的小管福寿螺发生了较高比例的种间遗传渐渗, 纯种比例为23.08%, 杂种包含53.84%的CB型和23.08%的CM型个体。Yang等[12]针对全国242头小管福寿螺开展种间基因渐渗研究, 发现CM型仅在南部沿海种群中检测到。与全国小管福寿螺相比, 西湖福寿螺中CM型所占比例明显增高。此外, 在日本、韩国野外种群中发现的CM杂交型也很少[12, 17], 认为可能与斑点福寿螺耐低温能力显著低于小管福寿螺有关[19, 34]。例如, 在日本低温北部仅有小管福寿螺分布, 斑点福寿螺分布于南部及中部且发生比例向高纬度逐渐降低, 杂种相较于斑点福寿螺更向北扩散[17]。西湖福寿螺CM型福寿螺较多, 可能是由于相对于浅水的稻田等湿地系统, 湖泊具有更稳定的生境, 且在淤泥底部冬季仍可以保持较高的温度[35]。
福寿螺幼螺和成螺虽具有一定的主动扩散能力, 但自然扩散以随水流传播的方式为主。西湖苏堤东西两侧湖水通过宽5—10 m的六个桥孔相连, 受风生流、风浪和进出水口等多种因素综合影响流场流态较复杂[36]。此外, 2003年实施的西湖引水工程(从钱塘江引水至西湖), 设有6个进水口和9个出水口, 福寿螺随水流在湖内扩散及向周边河流、公园和湿地等扩散的风险较大。鉴于福寿螺对景观及湿地生态系统的严重危害, 应加强对景区的发生监测, 并采取措施加强防控。例如可通过设置防螺网, 阻止福寿螺随水流扩散。此外, 应加强防控福寿螺的人为扩散风险, 如在景观植物引种时, 要仔细观察去除种苗茎干上附着的卵块和螺, 设定无福寿螺发生的码头及通行桥孔, 设立警示牌提醒游客请勿携带等。
4. 结论
当前, 西湖水域发生的福寿螺集中分布在苏堤以西, 岳湖水域发生密度最高。西湖分布有小管福寿螺和斑点福寿螺两种, 种群遗传多样性低。西湖福寿螺中较不耐低温的CM型比例较高, 可能与湖泊生境冬季淤泥底部可以保持较高的温度有关。西湖景区福寿螺随水流及人为携带在湖内及向湖外水域扩散的风险较高, 因而有必要加强对景区福寿螺的发生监测, 并采取措施防控福寿螺的扩散。
致谢:
感谢高跃、周焱、金钰婕、郭雨涵、刘浩文、刘苏汶、贺超、蒲佳佳在福寿螺发生调研及样品采集中的贡献, 感谢庞利铮在植物鉴定中予以的帮助!
-
![]()
图 1 西湖水域的福寿螺形态图
A—C. 雌螺; D. 雄螺, 均为小管福寿螺
Figure 1. Morphology of apple snails from the West Lake
A—C. female; D. male, four snails are all P. canaliculata
![]()
图 2 西湖水域福寿螺的分布情况
坐标符号表示福寿螺的发生地点; 黑色箭头为进水口, 1—6分别为莲花峰、小南湖、浴鹄湾、乌龟潭、茅家埠1、茅家埠2进水口; 灰色箭头为出水口, 7—15分别为岳湖、北里湖、圣塘闸、五公园、一公园、大华饭店、涌金闸、涌金池、柳浪出水口
Figure 2. Distribution of the apple snails in the West Lake
The coordinates indicate the location of the apple snails; the black arrow indicate the inlet of the West Lake, 1—6. Lianhuafeng, Xiaonanhu, Yuhuwan, Wuguitan, Maojiabu 1, Maojiabu 2 the grey arrow indicate the outlet of the West Lake, 7—15. Yuehu, Beilihu, Shengtangzha, Wugongyuan, Yigongyuan, Dahuafandian, Yongjinzha, Yongjinchi, Liulang
![]()
图 3 西湖福寿螺线粒体COⅠ基因单倍型变异位点
Figure 3. Variation sites of mitochondrial COⅠ gene haplotype in the apple snails from the West Lake
![]()
图 4 西湖水域福寿螺的单倍型及GenBank序列构建的系统发育树
系统发育分支节点依次为NJ/ML/BI置信值, 仅显示>50%的数值
Figure 4. The phylogenetic trees reconstructed by the haplotypes of apple snails from the West Lake and the GenBank sequences
The branch nodes of phylogenetic trees denote NJ/ML/BI bootstrap supports, and only>50% values are displayed
![]()
图 5 基于EF1α基因多重PCR扩增的西湖福寿螺杂交渐渗型检测
M为50 bp DNA ladder, 1—5为西湖水域采集的小管福寿螺样品, 分别为C, M, B, B, M型
Figure 5. Introgressive hybridization of the apple snails from the West Lake based on EF1α gene multiplex PCR method
M. 50 bp DNA ladder, 1—5: samples of P. canaliculata collected in the West Lake, represent C, M, B, M type
表 1 西湖水域福寿螺样品、序列及单倍型分布信息表
Table 1 Samples, sequences, and haplotype distribution of the apple snails from the West Lake
日期
Date数量(成螺/卵块)
No. (Adult/Egg clutch)GenBank序列号
Accession No.Hap1 Hap2 Hap3 Hap4 Hap5 Hap6 参考文献
Reference2020.05 5/14 MW051334-353 14 4 0 0 1 0 本研究This study 2019.05 23/17 MW051295-334 35 4 1 0 0 0 本研究This study 2018.05 13/0 MW051282-294 8 4 0 1 0 0 本研究This study 2017.07 0/10 MW051272-281 1 3 6 0 0 0 本研究This study 2016.11 26/2 MK189337-62, MH211409-10 17 9 0 0 0 2 [24] 2015.06 0/10 KT852718-20, 31-36, 61 4 6 0 0 0 0 [25] 总计Total 67/53 / 79 30 7 1 1 2 / 
下载: 导出CSV
表 2 西湖福寿螺单倍型序列与GenBank下载序列的相似性及碱基差异数
Table 2 Similarity and number of variation sites among the haplotypes and sequences retrieved from GenBank
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 P. canaliculata_AR 0 1 4 29 30 30 30 56 55 60 60 57 56 56 2 Hap1 100 1 4 29 30 30 30 56 55 60 60 57 56 56 3 Hap4 99.83 99.83 5 30 31 31 31 57 56 61 61 58 57 57 4 Hap5 99.33 99.33 99.16 31 32 32 32 56 56 60 60 58 57 57 5 P. canaliculata_JP 95.13 95.13 94.97 94.80 12 12 1 60 57 60 60 55 54 54 6 P. canaliculata_CN 94.97 94.97 94.80 94.63 97.99 0 13 62 59 58 58 56 55 55 7 Hap2 94.97 94.97 94.80 94.63 97.99 100 13 62 59 58 58 56 55 55 8 Hap3 94.97 94.97 94.80 94.63 99.83 97.82 97.82 59 56 59 59 56 55 55 9 P. maculata_BR 90.60 90.60 90.44 90.60 89.93 89.60 89.60 90.10 9 10 10 41 40 40 10 P. maculata_JP 90.77 90.77 90.60 90.60 90.44 90.10 90.10 90.60 98.49 9 9 36 35 35 11 P. maculata_CN 89.93 89.93 89.77 89.93 89.93 90.27 90.27 90.10 98.32 98.49 0 35 34 34 12 Hap6 89.93 89.93 89.77 89.93 89.93 90.27 90.27 90.10 98.32 98.49 100 35 34 34 13 P. occulta_CN1 90.44 90.44 90.27 90.27 90.77 90.60 90.60 90.60 93.12 93.96 94.13 94.13 1 3 14 P. occulta_CN2 90.60 90.60 90.44 90.44 90.94 90.77 90.77 90.77 93.29 94.13 94.30 94.30 99.83 2 15 P. occulta_CN3 90.60 90.60 90.44 90.44 90.94 90.77 90.77 90.77 93.29 94.13 94.30 94.30 99.50 99.66 注: 对角线以下为序列相似性, 对角线以上为碱基差异数。AR、BR、CN、JP分别表示阿根廷、巴西、中国和日本地理种群Note: Data below the diagonal are similarity and those above the diagonal are variable sites. AR, BR, CN, and JP indicate the populations of Argentina, Brazil, China, and Japan
下载: 导出CSV
表 3 小管福寿螺种群内的遗传多样性分析
Table 3 Population genetic diversity of P. canaliculata from the West Lake
日期
Date单倍型
多样性
Haplotype diversity (Hd)核苷酸
多样性
Nucleotide diversity (π)核苷酸平均
差异数
Average number
of nucleotide difference (k)Tajima’s D中性检验
Tajima’s D neutrality test2020.05 0.433 0.018 10.901 0.616 2019.05 0.229 0.011 6.797 –0.764 2018.05 0.564 0.023 14.000 1.764 2017.07 0.600 0.019 11.200 –0.696 2016.11 0.471 0.024 14.123 2.956** 2015.06 0.533 0.027 16.000 2.443** 总计Total 0.488 0.023 13.494 2.334* 注: *表示显著偏离中性(P<0.05), **表示P<0.01Note: *indicate significantly deviated from neutrality (P<0.05), **represent P<0.01
下载: 导出CSV
-
[1] Cowie R H, Hayes K A, Thiengo S C. What are Apple Snails? Confused Taxonomy and Some Preliminary Resolution [M]//Joshi R C, Sebastian L S (Eds.), Global Advances in Ecology and Management of Golden Apple Snails. Munñoz: Philippine Rice Research Institute, 2006: 3-23.
[2] Lowe S, Browne M, Boudjelas S, et al. 100 of the World’s Worst Invasive Alien Species, a Selection from the Global Invasive Species Database [M]. Auckland, New Zealand: The Invasive Species Specialist Group (ISSG) a Specialist of the Species Survival Commission (SSC) of the World Conservation Union (IUCN), 2004: 6-7.
[3] 尹绍武, 颜亨梅, 许芳. 金宝螺营养成份的初步研究 [J]. 水生生物学报, 1999, 23(6): 744-746. Yin S W, Yan H M, Xu F. Studies on nutrition of Ampullaria gigas [J]. Acta Hydrobiologica Sinica, 1999, 23(6): 744-746.
[4] 杨叶欣, 胡隐昌, 李小慧, 等. 福寿螺在中国的入侵历史、扩散规律和危害的调查分析 [J]. 中国农学通报, 2010, 26(5): 245-250. Yang Y X, Hu Y C, Li X H, et al. Historical invasion, expansion process and harm investigation of Pomacea canaliculata in China [J]. Chinese Agricultural Science Bulletin, 2010, 26(5): 245-250.
[5] Yang Q Q, Liu S W, He C, et al. Distribution and the origin of invasive apple snails, Pomacea canaliculata and P. maculata (Gastropoda: Ampullariidae) in China [J]. Scientific Reports, 2018, 8(1): 1185. doi: 10.1038/s41598-017-19000-7
[6] 俞晓平, 和田节, 李中方, 等. 稻田福寿螺的发生和治理 [J]. 浙江农业学报, 2001, 13(5): 247-252. doi: 10.3969/j.issn.1004-1524.2001.05.002 Yu X P, He T J, Li Z F, et al. Occurrence of golden apple snail, Pomacea canaliculata (Lamarck), in paddy fields and its management [J]. Acta Agriculture Zhejianggensis, 2001, 13(5): 247-252. doi: 10.3969/j.issn.1004-1524.2001.05.002
[7] 章家恩, 赵本良, 罗明珠, 等. 外来生物福寿螺入侵的生态风险及其评价探讨 [J]. 佛山科学技术学院学报 (自然科学版), 2010, 28(5): 1-6. Zhang J E, Zhao B L, Luo M Z, et al. The ecological risk analysis and assessment on the invasion of golden apple snails [J]. Journal of Foshan University (
Natural Science Edition ) , 2010, 28(5): 1-6. [8] Kwong K L, Chan R K Y, Qiu J W. The potential of the invasive snail Pomacea canaliculata as predator of various life-stages of five species of freshwater snails [J]. Malacologia, 2009, 51(2): 343-356. doi: 10.4002/040.051.0208
[9] Kwong K L, Dudgeon D, Wong P K, et al. Secondary production and diet of an invasive snail in freshwater wetlands: implications for resource utilization and competition [J]. Biological Invasions, 2010, 12(5): 1153-1164. doi: 10.1007/s10530-009-9537-x
[10] Lü S, Zhang Y, Steinmann P, et al. Emerging angiostrongyliasis in the mainland of China [J]. Emerging Infectious Diseases, 2008, 14(1): 161-164. doi: 10.3201/eid1401.061529
[11] Hayes K A, Cowie R H, Thiengo S C, et al. Comparing apples with apples: clarifying the identities of two highly invasive Neotropical Ampullariidae (Caenogastropoda) [J]. Zoological Journal of the Linnean Society, 2012, 166(4): 723-753. doi: 10.1111/j.1096-3642.2012.00867.x
[12] Yang Q Q, He C, Liu G F, et al. Introgressive hybridization between two non-native apple snails in China: widespread hybridization and homogenization in egg morphology [J]. Pest Management Science, 2020, 76(12): 4231-4239. doi: 10.1002/ps.5980
[13] Rawlings T A, Hayes K A, Cowie R H, et al. The identity, distribution, and impacts of non-native apple snails in the continental United States [J]. BMC Evolutionary Biology, 2007, 7(1): 97. doi: 10.1186/1471-2148-7-97
[14] Lü S, Zhang Y, Liu H X, et al. Phylogenetic evidence for multiple and secondary introductions of invasive snails: Pomacea species in the People’s Republic of China [J]. Diversity and Distributions, 2013, 19(2): 147-156. doi: 10.1111/j.1472-4642.2012.00924.x
[15] Yang Q Q, Liu S W, He C, et al. Invisible apple snail invasions: importance of continued vigilance and rigorous taxonomic assessments [J]. Pest Management Science, 2019, 75(5): 1277-1286. doi: 10.1002/ps.5241
[16] Yang Q Q, Yu X P. A new species of apple snail in the genus Pomacea (Gastropoda: Caenogastropoda: Ampullariidae) [J]. Zoological Studies, 2019(58): e13.
[17] Matsukura K, Okuda M, Cazzaniga N J, et al. Genetic exchange between two freshwater apple snails, Pomacea canaliculata and Pomacea maculata invading East and Southeast Asia [J]. Biological Invasions, 2013, 15(9): 2039-2048. doi: 10.1007/s10530-013-0431-1
[18] Glasheen P M, Burks R L, Campos S R, et al. First evidence of introgressive hybridization of apple snails (Pomacea spp. ) in their native range [J]. Journal of Molluscan Studies, 2020, 86(2): 96-103. doi: 10.1093/mollus/eyz035
[19] Matsukura K, Izumi Y, Yoshida K, et al. Cold tolerance of invasive freshwater snails, Pomacea canaliculata, P. maculata, and their hybrids helps explain their different distributions [J]. Freshwater Biology, 2016, 61(1): 80-87. doi: 10.1111/fwb.12681
[20] Leducq J B, Nielly-Thibault L, Charron G, et al. Speciation driven by hybridization and chromosomal plasticity in a wild yeast [J]. Nature Microbiology, 2016, 1(2): 15003.
[21] 田代科, 李春, 肖艳, 等. 中国秋海棠属植物的自然杂交发生及其特点 [J]. 生物多样性, 2017, 25(6): 654-674. doi: 10.17520/biods.2017050 Tian D K, Li C, Xiao Y, et al. Occurrence and characteristics of natural hybridization in Begonia in China [J]. Biodiversity Science, 2017, 25(6): 654-674. doi: 10.17520/biods.2017050
[22] 邓晓玲, 刘德清, 郭学飞, 等. 洱海流域湖泊湖滨带福寿螺时空分布 [J]. 生态毒理学报, 2018, 13(4): 134-142. doi: 10.7524/AJE.1673-5897.20180527003 Deng X L, Liu D Q, Guo X F, et al. Temporal and spatial distribution of Pomacea canaliculata in Er’hai Lake catchment [J]. Asian Journal of Ecotoxicology, 2018, 13(4): 134-142. doi: 10.7524/AJE.1673-5897.20180527003
[23] Folmer O, Black M, Hoeh W, et al. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates [J]. Molecular Marine Biology and Biotechnology, 1994(3): 294-299.
[24] Yang Q Q, Liu S W, Li J N, et al. Microsatellite evidence for multiple paternity in non-native populations of Pomacea canaliculata (Caenogastropoda: Ampullariidae) in China [J]. Aquatic Invasions, 2019, 14(4): 656-666. doi: 10.3391/ai.2019.14.4.06
[25] Fang L, Wong P K, Lin L, et al. Impact of invasive apple snails in Hong Kong on wetland macrophytes, nutrients, phytoplankton and filamentous algae [J]. Freshwater Biology, 2010, 55(6): 1191-1204. doi: 10.1111/j.1365-2427.2009.02343.x
[26] 叶建人, 冯永斌, 林贤文, 等. 福寿螺的寄主植物及其对福寿螺体重的影响 [J]. 生物安全学报, 2011, 20(2): 124-131. Ye J R, Feng Y B, Lin X W, et al. Host plants and their effects on body mass in the golden apple snail, Pomacea canaliculata Lamarck [J]. Journal of Biosafety, 2011, 20(2): 124-131.
[27] Qiu J W, Kwong K L. Effects of macrophytes on feeding and life-history traits of the invasive apple snail Pomacea canaliculata [J]. Freshwater Biology, 2009, 54(8): 1720-1730. doi: 10.1111/j.1365-2427.2009.02225.x
[28] Wong P K, Kwong K L, Qiu J W. Complex interactions among fish, snails and macrophytes: implications for biological control of an invasive snail [J]. Biological Invasions, 2009, 11(10): 2223-2232. doi: 10.1007/s10530-008-9378-z
[29] 杨倩倩, 刘苏汶, 茹炜岽, 等. 基于DNA条形码技术对浙江省外来入侵福寿螺进行分子鉴定 [J]. 生物多样性, 2016, 24(3): 341-350. Yang Q Q, Liu S W, Ru W D, et al. Molecular identification of invasive golden apple snails in Zhejiang Province based on DNA barcoding [J]. Biodiversity Science, 2016, 24(3): 341-350.
[30] Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J]. Genetics, 1989, 123(3): 585-595. doi: 10.1093/genetics/123.3.585
[31] Sekino M, Hara M, Taniguchi N. Loss of microsatellite and mitochondrial DNA variation in hatchery strains of Japanese flounder Paralichthys olivaceus [J]. Aquaculture, 2002, 213(1-4): 101-122. doi: 10.1016/S0044-8486(01)00885-7
[32] Sanchis A, Michelena J, Latorre A, et al. The phylogenetic analysis of variable-length sequence data: elongation factor-1α introns in European populations of the parasitoid wasp genus Pauesia (Hymenoptera: Braconidae: Aphidiinae) [J]. Molecular Biology and Evolution, 2001, 18(6): 1117-1131. doi: 10.1093/oxfordjournals.molbev.a003882
[33] Kim S, Lee W, Lee E. Estimation of a new molecular marker of the genus Stathmopoda (Lepidoptera: Stathmopodidae): comparing EF1α and COI sequences [J]. Journal of Asia-Pacific Entomology, 2017, 20(1): 269-280. doi: 10.1016/j.aspen.2016.12.002
[34] Yoshida K, Matsukura K, Cazzaniga N J, et al. Tolerance to low temperature and desiccation in two invasive apple snails, Pomacea canaliculata and P. maculata (Caenogastropoda: Ampullariidae), collected in their original distribution area (northern and central Argentina) [J]. Journal of Molluscan Studies, 2014, 80(1): 62-66. doi: 10.1093/mollus/eyt042
[35] Ito K. Environmental factors influencing overwintering success of the golden apple snail, Pomacea canaliculata (Gastropoda: Ampullariidae), in the northernmost population of Japan [J]. Applied Entomology and Zoology, 2002, 37(4): 655-661. doi: 10.1303/aez.2002.655
[36] 尤爱菊, 滑磊, 韩曾萃, 等. 杭州西湖三维流场特性的现场观测及模拟研究 [J]. 水力发电学报, 2017, 36(9): 111-120. doi: 10.11660/slfdxb.20170912 You A J, Hua L, Han Z C, et al. Field observation and simulation study of three-dimensional flows in West Lake, Hangzhou [J]. Journal of Hydroelectric Engineering, 2017, 36(9): 111-120. doi: 10.11660/slfdxb.20170912
-
期刊类型引用(6)
1. 李锦，杨玲玲，赖利平，周昭敏，罗怡，靳联娟. 南充市福寿螺种类与遗传多样性分析. 西华师范大学学报(自然科学版). 2024(04): 343-349 . 
百度学术
2. 蒋玲，赵永波，李天美. 我国外来入侵生物福寿螺的物种鉴定. 中国寄生虫学与寄生虫病杂志. 2024(03): 384-388+398 . 百度学术
3. 高雪，赵志淼，鲁仙，高雪晴，吴泓权，张饮江. 中华天胡荽、菖蒲水提液对福寿螺的毒性影响及损伤机理. 水生生物学报. 2024(11): 1862-1870 . 本站查看
4. 郭靖，李娟，张少斌，彭红元，贺爱兰. 广东省稻田福寿螺的发生与防治情况——基于全省597份问卷的调查分析. 中南农业科技. 2023(08): 133-138 . 百度学术
5. 罗鑫，张君倩，刘洋，张多鹏，白雪，杨佳莉，谢志才，刘金兰. 桂江流域福寿螺种群结构时空分布特征及其驱动机制. 生态学报. 2023(21): 9048-9060 . 百度学术
6. Wang Xin，Zhao Congcong，Liu Yang. Review of Restoration Technology and Multi-environmental Factors Analysis of West Lake Ecosystem in Hangzhou City. Journal of Northeast Agricultural University(English Edition). 2022(03): 37-47 . 必应学术
其他类型引用(6)
计量
- 文章访问数: 2166
- HTML全文浏览量: 605
- PDF下载量: 141
- 被引次数: 12