倒刺鲃两极虫的宿主新记录及其分子系统学研究

向 瑶 张金叶 赵元莙

(重庆师范大学生命科学学院, 动物生物学重庆市重点实验室, 重庆 401331)

黏孢子虫是一类个体微小结构简单的后生动物寄生虫, 主要寄生于鱼类, 有些种类能引起鱼类宿主的病害[1—3], 严重的甚至会导致死亡[4—7]。其中两极虫科(Myxidiidae Thélohan), 1892黏孢子虫与大多数黏孢子虫类群不同, 具有最广泛的宿主类群,该科的种类不仅能寄生于变温动物(海、淡水鱼类、两栖类、爬行类)[8—12], 还能寄生于恒温动物(鸟类、哺乳类)[13,14]。有些物种还能广泛寄生于多种不同的宿主, 如李氏肠黏虫(Enteromyxum leeiDiamant, Lom and Dykova), 1994的宿主已经超过40种[15]。两极虫属(MyxidiumButschli, 1882)隶属黏体门(Myxozoa)、黏孢子虫纲(Myxosporea)、双壳目(Bivalvulida)、两极虫科(Myxidiidae)[16], 是两极虫科中种类最多的属, 迄今已记录的种类超过230种[9]。两极虫属种类具如下特征: 2个孢子壳瓣,孢子呈纺锤形、梨形或椭圆形, 末端稍尖; 壳瓣表面光滑或具条纹, 孢子缝线直; 孢子的两极端各具一个梨形或球形的极囊且极囊向两端开口, 细长的极丝在极囊内横向盘曲; 孢子质位于两极囊之间。两极虫为典型的腔寄生类群, 绝大多数种类寄生于宿主的胆囊或泌尿系统, 极少数可营组织寄生[2,17],个别组织寄生种类甚至能导致寄生部位出现空泡影响宿主发育和繁殖[18]。

胭脂鱼(Myxocyprinus asiaticusBleeker, 1864),隶属鲤形目(Cypriniformes)胭脂鱼科(Catostomidae)胭脂鱼属(Myxocyprinus), 是胭脂鱼科唯一分布在我国的特有种[19,20], 现已在我国驯化为养殖鱼类,成为我国长江流域重要经济鱼类之一[21,22]。胭脂鱼是杂食性的底栖鱼类, 在人工养殖过程中常与其他不同食性和不同水层栖息的鱼类混养以充分利用不同水层的资源[23]。但混养易增加胭脂鱼与其他不同鱼类的密切接触机会, 加大寄生虫或其他病原传播的风险, 进而增高传染性造成经济损失[24]。因此开展对胭脂鱼的寄生虫及其病害的研究也是水产养殖中鱼病防治的需要。

倒刺鲃两极虫(Myxidium spinibarbaChen,et al.), 2020年最初从中华倒刺鲃[Spinibarbus sinensis(Bleeker, 1871)]的胆囊中检获[25]。本研究检获的倒刺鲃两极虫分离自胭脂鱼胆囊, 在对其重描述的基础上, 结合其形态学与18S rDNA分子数据, 比较分析寄生于不同宿主的倒刺鲃两极虫及其他相似种之间的分子系统学关系, 旨在为寄生胭脂鱼的黏孢子虫种类组成及两极虫的宿主多样性积累基础数据。

1 材料与方法

1.1 黏孢子虫样本采集和鉴定

宿主鱼于2019年10月11日从重庆北碚区某养殖基地(该养殖基地由8个8.35 m×5.66 m的养殖池相连, 且池水互通; 多种鱼类以不同年龄、不同规格同池或者单独轮养轮放)捕获。共获26尾无明显病症的胭脂鱼带回实验室解剖并进行黏孢子虫检测, 其中4尾鱼的胆囊检获到黏孢子虫(感染率为15.4%)。黏孢子虫收集、处理、测量和鉴定方法参照文献[26]处理。

1.2 DNA提取、PCR扩增和测序

从上述同一批次不同胭脂鱼个体的胆囊内分离获得游离黏孢子虫, 加入95%的乙醇溶液置于1.5 mL离心管保存。从离心管底部吸取10 μL黏孢子虫液体, 用超纯水清洗2—3次, 离心(2000 ×g)沉淀后用于DNA提取。基因组DNA的提取严格参照DNeasy Blood & Tissue试剂盒(QIAGEN, Hilden, Germany)说明书执行, 将提取的基因组DNA置于-20℃保存、备用。黏孢子虫18S rDNA序列通过PCR技术利用引物18e: 5′-CTGGTTGATCCTGCC-3′[27]和18R: 5′-CTACGGAAACCTTGTTACG-3′[28]扩增。反应体系和反应条件参照文献[29]。取3 μL扩增产物用于1%琼脂糖凝胶进行电泳检测, 目的产物经胶回收试剂盒(OMEGA, 美国)回收并送到上海英潍捷基贸易有限公司测序。

1.3 系统发育分析

将本研究获得的序列, 在GenBank中用BLAST进行同源序列检索比对。根据比对结果, 选取12条高同源性序列(序列相似度高于89%), 24条先前相关研究所用序列[25]及本研究所获的2条序列进行系统发育分析。另以Tetracapsuloides bryosalmonae(KF731712)和Buddenbrockia plumatellae(KF731691)为外群, 利用Mr Bayes 3.2.6[30]选用最优拟合模型(GTR)运行1000000代构建Bayes(BI)树, 利用在线软件CIPRES Science Gateway V.3.3(http://www.phylo.org/)选用RAxML-HPC2 XSEDE(8.2.12)模型构建ML(Maximum likelihood)树。最后用FigTree v1.4和Adobe Photoshop CC软件编辑注释系统树。将相关18S rDNA序列用MEGA 6.0[31]设置执行K2P模型多序列比对计算出遗传距离, 利用在线序列比对工具(http://www.ebi.ac.uk/Tools/psa/)计算得出两两序列相似度, 结合MEGA和Bioedit软件完成变异位点分析。

1.4 主成分分析

使用PAST 3软件进行主成分分析(Principal Component Analysis, PCA), 对不同宿主寄生的倒刺鲃两极虫株系的孢子长、孢子宽、极囊长、极囊宽和极囊间距5项形态学特征进行形态计量学比较[32]。利用变异协方差矩阵生成置信度为95%的散射图。

2 结果

2.1 倒刺鲃两极虫的形态学重描述

取宿主鱼的适量胆汁制片镜检, 未见营养体和未成熟孢子, 但可见大量成熟孢子。成熟孢子壳面呈梭形, 两端稍尖。孢子缝面呈卵圆形, 缝脊平直,孢子壳瓣上有7—10条与缝脊平行的条纹(图 1)。孢子长为(12.59±0.56) μm (11.55—13.36 μm), 孢子宽为(5.92±0.40) μm (5.11—6.86 μm)。两个极囊呈梨形位于孢子两端, 大小相等, 其开口倾向孢子两端。极囊长为(3.48±0.36) μm (2.63—4.22 μm), 极囊宽为(2.83±0.26) μm (2.29—3.21 μm)。两极囊间隔(3.66±0.29) μm (3.12—4.20 μm), 极丝盘曲4—5圈(表 1)。

表1 寄生不同宿主的倒刺鲃两极虫的特征比较Tab. 1 Comparison among Myxidium spinibarba from different hosts

图1 寄生于胭脂鱼的倒刺鲃两极虫孢子显微图片Fig. 1 Photomicrographs of Myxidium spinibarba from Myxocyprinus asiaticus

2.2 18S rDNA分子特征和系统发育分析

本研究从不同个体的胭脂鱼样本中获得2条18S rDNA序列, 长度分别为1785(MT775466)和1913 nt (MT775467)。比对结果显示, 所获两条18S rDNA序列之间相似度为100.00%, 在可比对长度(1785 nt)上无变异位点。两条序列在NCBI中BLAST分析显示: 与寄生于中华倒刺鲃的倒刺鲃两极虫(MH766654)的相似度最高, 为99.9%—100.0%,MT775467与Myxidium spinibarba(MH766654)在可比对长度范围内(1913 nt)仅有1个变异位点; 其次与Myxidiumsp.(MK913426)和Zschokkellasp.(KM 401441)具有较高相似度(表 2)。

表2 黏孢子虫基于18S rDNA序列遗传距离与相似度Tab. 2 Genetic distances and similarities of myxosporeans based on 18S rDNA sequences

基于18S rDNA序列构建的ML和BI树显示出一致的拓扑结构。两棵系统树均显示: 所选取黏孢子虫序列被分成A、B两个大分支(图 2)。A支是一混合支系, 包括寄生于爬行动物的两极虫、海水鱼类的两极虫和楚克拉虫及淡水鱼类的两极虫和楚克拉虫; 而B支则完全由寄生海水鱼类的两极虫组成。本研究株系位于A支内, 其中株系(MT775466)与倒刺鲃两极虫(MH766654)聚支后, 再与株系(MT775467)聚支, 形成倒刺鲃两极虫支系; 然后, 与Zschokkellasp. (KM401441)和Myxidiumsp. (MK913426)亚支再形成了姐妹支。

图2 基于18S rDNA序列构建的BI/ML系统发育树Fig. 2 Phylogenetic tree generated by BI/ML based on 18S rDNA sequences

2.3 基于寄生不同宿主的倒刺鲃两极虫形态学数据的主成分分析

主成分分析结果显示本研究两株系(MT775466和MT775467)与倒刺鲃两极虫(MH766654; 该株系原始数据来自本实验室)在5项形态特征结合生成的椭圆散点图上均有重合, 其中●株系(MT775466)散点图基本覆盖*株系(MH766654)散点图, ▲株系(MT775467)的散点图仅与前两者的散点图部分重叠(图 3)。

3 讨论

本研究两株系形态学和生态学特性均显示一致, 应属于同一物种。本研究获得的两极虫, 两端稍尖的梭形成熟孢子、平直的缝脊和壳瓣上的条纹数、两个大小相等的梨形极囊、极丝圈数及极丝开口方向等形态学特征均与倒刺鲃两极虫的形态学特征一致[25], 且两者的感染部位均为胆囊。同时, 形态学数据的主成分分析结果也显示: 尽管株系(MT775467)有略微的差异, 但本研究的2个株系的形态测量数据基本都与倒刺鲃两极虫的形态测量数据有大幅重叠(图 3), 仅存在少许差异, 这样的形态学量度差异在种内变化范围内[33]。综合形态学与主成分分析结果, 本研究获得的两极虫与倒刺鲃两极虫当属同种。

图3 寄生不同宿主倒刺鲃两极虫的孢子形态差异主成分分析Fig. 3 Principal component analysis on morphological differences of Myxidium spinibarba from different hosts

鉴于黏孢子虫结构简单, 学者们在以往的研究中认为在对黏孢子虫进行物种鉴定和描述时, 应当综合考虑该黏孢子虫的形态学特征、生态学特征(宿主、寄生部位、地理分布)和分子特征[34,35]。基于此, 尽管本研究株系的宿主与原始描述宿主不同,但其18S rDNA序列比对结果显示寄生胭脂鱼的两株系与中华倒刺鲃株系之间仅有1个变异位点, 遗传距离和相似度分别为0.000—0.001和99.9%—100.0%, 符合已有研究给出的黏孢子虫变异位点不超过10个, 遗传距离在0.000—0.007, 相似度在98.6%—100%属于种内水平的判别标准[36—39]。分子特征分析也表明本研究所获寄生胭脂鱼的两极虫株系(MT77546和MT775467)与寄生中华倒刺鲃的倒刺鲃两极虫(MH766654)为同一物种。

系统发育分析显示: 寄生不同宿主的倒刺鲃两极虫株系嵌合聚枝, 寄生胭脂鱼的倒刺鲃两极虫株系(MT775467)先分化, 并与前述的主成分分析和18S rDNA特征分析结果表现一致。同种两极虫可在亲缘关系较远的鱼类宿主寄生, 这与已有的关于弧形虫和肠黏虫的研究有着明显的一致性[33,40]。在本研究所显示的倒刺鲃两极虫姐妹支中, 寄生于露斯塔野鲮[Labeo rohitaHamilton, 1822(鲤形目)]的Zschokkellasp.(KM401441)与寄生于长臀[Cranoglanis bouderiusRichardson, 1846(鲇形目)]的Myxidiumsp.(MK913426)形成高支持率(0.99/84)的系统发生关系也得到了旁证: 尽管Zschokkellasp.(KM401441)与Myxidiumsp.(MK 913426)为不同种类且宿主亲缘关系较远, 但显示形态学特征较为相似, 都享有梭形的孢子并具有梨形极囊(图 2)。两极虫的系统发生关系与孢子形态、腔寄生的寄生部位趋向性等可能更相关。在通常情况下, 内寄生虫的进化在很大程度上受其与宿主相互作用的影响, 基于黏孢子虫SSU rDNA的大数据集系统发育分析结果认为组织寄生物种是由类似祖先型的腔寄生黏孢子虫进化而来[41,42]。寄生鱼类的碘泡虫类黏孢子虫多为组织寄生, 是黏孢子虫中进化较晚的类群[42]。绝大多数碘泡虫类的黏孢子虫寄生于淡水鱼且具有明显的宿主特异性乃至寄生部位特异性, 往往只感染一种宿主甚至只感染一种宿主的某一特定器官或组织[35]。而腔寄生的两极虫类是相对组织寄生的碘泡虫类较为祖先的类群, 进化过程中分化形成海水支系和淡水支系, 其变异程度要比碘泡虫等组织寄生类群更大,更易于对环境及宿主的变化形成适应[41,42]。此外,先前基于鱼体混养、受感染鱼宿主的生活水体、经口饲喂和肛门插管等途径对李氏肠黏虫传播实验研究表明, 李氏肠黏虫可感染多种不同且亲缘关系较远的鱼类宿主[43], 也说明腔寄生的两极虫类黏孢子虫对宿主具有更广泛的适应性。鉴于此, 推测本研究中倒刺鲃两极虫也可能具有多宿主特征, 能寄生多种不同宿主, 但这还需要更多研究予以证实。

综上, 本研究首次在胭脂鱼体内检获黏孢子虫,胭脂鱼是倒刺鲃两极虫的新宿主; 寄生胭脂鱼的倒刺鲃两极虫株系(MT775467)先分化。本研究也为寄生胭脂鱼的黏孢子虫种类组成积累了新的数据。