扩展功能
文章信息
- 李超, 索栋, 杨承忠, 赵元莙
- LI Chao, SUO Dong, YANG Chengzhong, ZHAO Yuanjun
- 楔形两极虫(黏体门,双壳目)的形态学重描述及其近缘种的分子系统学研究
- Redescription of Myxidium cuneiforme Fujita, 1924 (Myxosporea: Bivalvulida) and Molecular Phylogeny with its Relative Species
- 四川动物, 2016, 35(3): 384-390
- Sichuan Journal of Zoology, 2016, 35(3): 384-390
- 10.11984/j.issn.1000-7083.20150350
-
文章历史
- 收稿日期: 2015-11-06
- 接受日期: 2016-01-15
Bütschli(1882)从淡水的白斑狗鱼Esox lucius膀胱中检获并报道李氏两极虫Myxidium lieberkuehni,从而建立了两极虫属Myxidium。两极虫孢子一般是纺锤形、直线形或新月形,甚至是S形,末端稍尖;壳瓣表面光滑或具条纹,孢子缝线直,极囊梨形,位于孢子的两极端,极丝细长,极丝在极囊内横向盘曲;孢子质位于两极囊之间(Lom & Dyková,2006)。迄今,该属已记述232种,是粘孢子虫中继碘泡虫属Myxobolus Bütschli,1882之后最大的一属(Eiras et al.,2011)。其主要营腔寄生,少数可行组织寄生。两极虫多寄生于海、淡水鱼类,两栖动物、爬行动物和哺乳动物则鲜有报道(宋微波等,2003;Canning & Okamura,2004;Prunescu et al.,2007)。
因技术条件限制,过去两极虫的分类鉴定主要以成熟孢子的形态特征作为依据。但是,粘孢子虫种类繁多,种间相似度高,加之在形态判定上主观因素的干扰,其分类及系统地位研究虽历史悠久但混乱较多。两极虫属、楚克拉虫属Zschokklella和弧形虫属Sphaeromyxa物种形态较为相似,其中,两极虫与楚克拉虫形态尤为相似,用于物种鉴定的形态特征常常不能将两者有效区分,这点已有学者对其提出了不同的看法(Diamant et al.,1994;Whipps & Zhao,2015)。
18S rDNA分布于所有的细胞生物体内,并且具有相对缓慢的进化历程。其中,保守区可用于构建所有生命的进化树,而变异区可用来区分属或者种。18S rDNA长度适中,通常为1 200~1 900 bp,其信息量充足,但又不过长且容易获得。因此,在不同物种的研究领域中得到了广泛应用。以18S rDNA为分子标记,学者们对包括粘孢子虫在内的低等动物的分类学和系统学问题进行了大量研究,解决了许多悬疑问题(Wright et al.,1997;Kent et al.,2001;钟霞等,2007;Zhao et al.,2013)。众多研究表明,18S rDNA是一种粘孢子虫等低等动物分类和系统发生研究的良好分子标记(Michel et al.,2003;鲁义善,聂品,2004;Lom & Dyková,2006;Bartošová et al.,2009)。
本研究在对采自涪江重庆市潼南县江段的楔形两极虫Myxidium cuneiforme Fujita,1924进行形态学重描述的基础上,利用18S rDNA为分子标记对其分子系统学进行研究,以期进一步厘清两极虫与楚克拉虫的分类学关系。
1 材料与方法 1.1 标本采集与虫种鉴定宿主鲤Cypyinus carpio于2008年11月采自涪江重庆市潼南县江段。20条鱼体中12条被感染粘孢子虫,感染率为60%。寄生粘孢子虫(楔形两极虫)标本的检获、观察、处理、测量、鉴定方法及相关术语参照赵元莙等(2001)完成。
1.2 DNA的提取、扩增和测序 1.2.1 DNA的提取从鲤胆囊获得游离的粘孢子虫经冲洗收集、离心富集,DNA提取前采用低速离心的方法,用灭菌蒸馏水清洗悬浮液2~3次以除去杂质。DNA抽提采用Dneasy Tissue Kit(QIAGEN)试剂盒,操作步骤均按照厂家提供的说明书进行。提取的基因组DNA置于-20 ℃冰箱保存备用。
1.2.2 DNA扩增和测序选取粘孢子虫通用扩增引物:ERIB1(5’-ACCTGGTTGATCCTGCCAG-3’)和ERIB10(5’-CTTCCGCAGGTTCACCTACGG-3’)(Ba-rta et al.,1997),对样本进行18S rDNA扩增。50 μL的PCR反应体系包括5 μL 10×Ex Taq buffer,40 μmol·L-1 dNTPs,20 μmol·L-1引物各1 μL,100 ng的模板DNA,1.5 μL Ex Taq酶(5 U·μL-1)(TaKaRa),用灭菌双蒸水补足至终体积50 μL。反应程序为:95 ℃预变性5 min;95 ℃变性1 min,48 ℃退火1 min,72 ℃延伸2 min,循环35次;最后72 ℃延伸10 min,反应完成后12 ℃保温。取3 μL扩增产物于1%的琼脂糖凝胶电泳检测,然后将PCR产物用胶回收试剂盒Gel Extraction Kit(OMEGA)进行纯化回收。最后将纯化产物送上海英骏生物技术有限公司测序(测序仪为ABI,3730XL)。
1.3 序列选取和系统分析将本研究的楔形两极虫重庆种群18S rDNA片段在GenBank中通过BLAST进行序列同源性比对。根据比对结果,选取同源性较高的序列以及其他相关类群的11条DNA序列(表 1),序列所涉及物种有两极虫属、楚克拉虫属、弧形虫属和鼩粘虫属Soricimyxum等4个属的种类,另选取鲑四囊虫Tetracapusuloides bryosalmonae作为外群构建系统发育树。通过Clustal W程序按照缺省参数进行序列多重比对,采用PAUP 4.0(Swofford,2001)基于GTR模型按照缺省参数完成最大似然(ML)系统树的构建,ML树运算100代。贝叶斯(BI)树通过MrBayes3.0b4(Ronquist & Huelsenbeck,2003)基于GTR模型构建而成,BI树运算1 000 000代,由Treeview32和Photoshop完成系统树的绘制。
粘孢子虫物种名 | GenBank 登录号 |
Myxidium cuneiforme | DQ377709 |
Zschokkella parasilui | DQ377689 |
Z. soleae | JX271832 |
Z. nova | DQ377688 |
M. truttae | AF201374 |
M. chelonarum | DQ377694 |
Sphaeromyxa hellandi | DQ377693 |
S. zaharoni | AY538662 |
Soricimyxum fegati | EU232760 |
M. bergense | DQ377702 |
M. laticurvum | JN232760 |
对所选序列两两之间的遗传距离利用MEGA6.0(Tamura et al.,2013)进行计算。所选粘孢子虫之间的相似度用GenBank的BLAST(bl2seq)工具计算获得。变异位点分析与GC含量的计算借助Bioedit和MEGA 6.0完成。
2 结果 2.1 楔形两极虫重庆种群的形态学重描述楔形两极虫孢子壳面观长条形,中部稍凹陷或孢子一边突出,一边稍凹;壳瓣上有6~8条与缝嵴平行的条纹。缝嵴直,缝面观呈梭形。极囊2个,呈梨形,分布于孢子两极端;极丝细长且明显,盘曲5~6圈。孢子量度(n=20)为:孢子长12.5 μm±0.3 μm(12.0~13.1 μm),孢子宽5.4 μm±0.3 μm(4.8~6.1 μm),极囊长4.3 μm±0.2 μm(4.0~4.6 μm),极囊宽3.1 μm±0.2 μm(2.6~3.4 μm)(图 1)。
2.2 楔形两极虫18S rDNA分子特征及其系统发生本研究所得到的楔形两极虫18S rDNA部分序列片段长度为1 622 nt。楔形两极虫重庆种群比武汉种群GC含量略低(46.58% vs. 46.97%)。楔形两极虫的18S rDNA经BLAST比对的结果显示,在可比范围内,本研究的楔形两极虫重庆种群与武汉种群(Fiala,2006)相似度高达99%。本种群与鲇楚克拉虫Zschokklella parasiluri、新楚克拉虫Z. nova以及鲑两极虫M. truttae的相似度为91%、90%、90%。本研究对序列相似度最高的2条序列,即楔形两极虫武汉种群和重庆种群进行比对分析,结果表明武汉种群和重庆种群共有8个变异位点(图 2)。
基于18S rDNA序列片段的遗传距离分析表明,本研究的楔形两极虫重庆种群与武汉种群的遗传距离为0.001,楔形两极虫重庆种群与鲑两极虫之间的遗传距离为0.046,楔形两极虫重庆种群和新楚克拉虫之间的遗传距离为0.046,鲑两极虫和新楚克拉虫之间的遗传距离为0.006(表 2)。
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | |
1 | 85% | 88% | 99% | 86% | 90% | 89% | 91% | 90% | 90% | 91% | 90% | |
2 | 0.231 | 85% | 85% | 92% | 84% | 84% | 87% | 88% | 85% | 84% | 85% | |
3 | 0.058 | 0.228 | 88% | 86% | 89% | 89% | 91% | 90% | 90% | 88% | 89% | |
4 | 0.001 | 0.233 | 0.059 | 86% | 89% | 89% | 92% | 91% | 90% | 91% | 90% | |
5 | 0.219 | 0.043 | 0.212 | 0.220 | 87% | 85% | 88% | 87% | 87% | 88% | 87% | |
6 | 0.046 | 0.238 | 0.045 | 0.047 | 0.221 | 88% | 92% | 91% | 99% | 90% | 92% | |
7 | 0.094 | 0.241 | 0.074 | 0.096 | 0.234 | 0.087 | 90% | 89% | 90% | 89% | 88% | |
8 | 0.096 | 0.218 | 0.075 | 0.097 | 0.211 | 0.087 | 0.110 | 89% | 89% | 91% | 90% | |
9 | 0.103 | 0.238 | 0.082 | 0.105 | 0.222 | 0.095 | 0.117 | 0.058 | 90% | 91% | 89% | |
10 | 0.046 | 0.235 | 0.045 | 0.047 | 0.228 | 0.006 | 0.082 | 0.096 | 0.095 | 91% | 92% | |
11 | 0.043 | 0.233 | 0.049 | 0.045 | 0.211 | 0.036 | 0.093 | 0.094 | 0.099 | 0.035 | 91% | |
12 | 0.050 | 0.233 | 0.048 | 0.052 | 0.222 | 0.039 | 0.093 | 0.087 | 0.090 | 0.038 | 0.040 | |
注Notes: 1. Myxidium cuneiforme CQ, 2. M. bergense, 3. M. chelonarum, 4. M. cuneiforme WH, 5. M. laticurvum, 6. M. truttae, 7. Soricimyxum fegati, 8. Sphaeromyxa hellandi, 9. S. zaharoni, 10. Zschokkella nova, 11. Z. parasiluri, 12. Z. soleae. |
利用ML和BI法分别构建了楔形两极虫及其相关种的18S rDNA系统树,两者呈一致的拓扑结构(图 3)。粘孢子虫首先聚类为两大支系,即海水支系(图 3,B支)和淡水支系(图 3,A支)。海水支系由海水两极虫(M. bergense、M. laticurvum)组成,而淡水支系则由两极虫、弧形虫、楚克拉虫和肝鼩粘虫Soricimyxum fegati组成。两极虫属和楚克拉虫属均非单系发生,两属间物种相互交叉聚支(图 3)。本研究得到的楔形两极虫位于淡水支系,其首先和楔形两极虫武汉种群聚在一起,且显示两者具有较高支持率(PP=1.00,BS=100)支持的同源性,它们形成的一支与鲇楚克拉虫Z. parasilui形成姐妹群关系,弧形虫支则位于淡水两极虫-楚克拉虫支系的基部。
通过分析得到该12条序列共有保守区33个,较淡水种类而言,海水两极虫整体上具有更多的变异位点,弧形虫的变异位点个数介于海水两极虫和淡水种类之间,肝鼩粘虫序列保守区更相似于淡水两极虫(图 4)。
3 讨论楔形两极虫在日本从鲤胆囊中检获并首次描述(Fujita,1924),随后,在前苏联,日本和中国的辽宁、杭州西湖、四川、重庆、湖北武汉等地都有过相关报道(陈启鎏,马成伦,1998)。本研究从鲤胆囊中检获两极虫,用20个样品进行了形态学数据的测量,采用统计学方法对各种形态特征进行描述,经比较发现与楔形两极虫源报道种群(Fujita,1924)所描述的形态特征基本一致,但孢子和极囊大小与原始描述略有差异(表 3),同时与《中国动物志》(陈启鎏,马成伦,1998)记载描述种群相比较,整体偏小(表 3)。这可能与宿主种类和地理分布等有关。然而比较遗憾的是,除本研究外,唯一公布有18S rDNA分子数据的楔形两极虫报道却未给出其形态学数据,因此不能与该种群进行形态比较。本研究综合了该物种的形态和分子数据,补充和丰富了粘孢子虫物种基础资料。
项目 | 本研究种群 | 源报道种群(Fujita,1924) | 《中国动物志》记载种群(陈启鎏,马成伦,1998) |
孢子长/μm | 12.5(12.0~13.1) | 13 | 15.5 |
孢子宽/μm | 5.4(4.8~6.1) | 6 | 5.8 |
极囊长/μm | 4.3(4.0~4.6) | 3 | 5.4 |
极囊宽/μm | 3.1(2.6~3.4) | — | 3.9 |
虫体表面的条纹数 | 6~8 | — | 8~9 |
宿主 | 鲤Cypyinus carpio haematopterus | 鲫Carassius vulgaris syn. Carassius auratus | 鲫Carassius auratus或镜鲤Cyprinus carpio carpio |
地理分布 | 重庆 | 秋田(日本) | 辽宁、重庆(长寿)、浙江、湖北 |
较保守序列的变异位点的个数可以反映生物体之间的关系,变异位点越少,2种生物体之间的亲缘关系可能越近,反之亦然。序列变异位点分析结果显示,在1 622 nt长度的18S rDNA序列上,本研究所得的楔形两极虫重庆种群与武汉种群之间共有8个变异位点,其主要表现为碱基的缺失,如在531位点、570位点、771位点、1281位点、1282位点、1283位点分别缺失1个碱基。据Zhao等(2013)对粘孢子虫18S rDNA序列的分子鉴别标准,本研究所得的粘孢子虫应为楔形两极虫,这从分子水平进一步验证了从形态上对该种鉴定的准确性。
冉佼等(2014)基于18S rDNA为分子标记的遗传距离的粘孢子虫分类研究认为,种内水平的遗传距离范围为0~0.042,其中大部分物种的遗传距离集中在0~0.007范围内;种间水平的遗传距离范围为0.010~0.323,其中大多数物种的遗传距离集中在0.135 2~0.260 4范围内;属间水平的遗传距离范围为0.105~0.296,其中大多数物种的遗传距离集中在0.181 4~0.257 8范围内。本研究中,选取与楔形两极虫重庆种群相似度最高的11种粘孢子虫进行18S rDNA遗传距离计算,结果显示:本研究得到的楔形两极虫重庆种群与武汉种群之间的遗传距离为0.001,为种内遗传水平(Zhao et al.,2013;冉佼等,2014),进一步证实了本研究所得物种为楔形两极虫。同时,序列相似度的分析得到了同样的结论。
本研究所选取的11种粘孢子虫利用BI/ML法构建的系统树一致显示,来自海水呈单系发生的弧形虫却位于淡水支(两极虫-楚克拉虫)的基部,这表明弧形虫与淡水两极虫有较近的亲缘关系。同时,保守区变异位点分析显示,弧形虫的序列差异整体上介于海水两极虫和淡水两极虫之间。弧形虫宿主大多生活于海洋性环境,其中一些有洄游习性(如Sphaeroxyum hellandi的宿主黑线鳕Melanogrammus aeglefinus),这可能对弧形虫在海-淡水中的演化有重要影响。这些表明,弧形虫在两极虫类群中地位较为特殊,其物种的形成与分化规律有待今后深入研究。18S rDNA保守区的变异位点分析(图 4)和系统发育分析(图 3)同时表明,陆生哺乳动物肝脏寄生的肝鼩粘虫与淡水粘孢子虫具有更近的亲缘关系,这与前人的研究结果一致(Dyková et al.,2007),这可能与其宿主鼩鼱Sorex araneus淡水饮水环境(Michielsen,1966)有关。
目前,粘孢子虫基于分子数据的系统发生关系与依据孢子形态的分类具有很多差异(张婧等,2015a,2015b)。根据两极虫与楚克拉虫形态的相似性和差异性,前人把两者归入两极科下的2个平行属(陈启鎏,马成伦,1998;Lom & Dyková,2006)。然而,两属间诸多分类学特征却甚为模糊,如两极虫缝线多为直线形,楚克拉虫缝线弯曲成"S"形,两极虫属的极囊多为梨形,楚克拉虫属的极囊多近似圆形。但是,依据这些特征并不能解决其所有的形态分类问题,如当极囊既不是严格的梨形也不是严格的圆形时,极囊形态作为一种分类依据就显得很混乱。类似这样的鉴别性描述还有很多,如极囊开口方向,缝线形状等。但系统发育分析显示这2属的物种相互交叉聚支(图 3),并不能将两者有效地分开。同时,淡水支系的两极虫与楚克拉虫在保守区的变异位点上也显示出其较近的亲缘关系(图 4)。基于以上形态和分子数据的比较研究,本研究支持Diamant等(1994)将两极虫属和楚克拉虫属合并为1个属的观点。