Study on ITS 2 molecular markers of six pearl oyster species in the genus Pinctada
-
摘要:
对珠母贝属的大珠母贝、珠母贝、白珠母贝、黑珠母贝、长耳珠母贝、黑珠母贝和合浦珠母贝6个种的内部转录间隔区2(ITS 2)序列及其两侧的5.8S和28S的部分序列进行了比较分析。其中黑珠母贝的序列来自GenBank。PCR扩增片段大小为600 bp左右，测序结果表明，ITS 2长211~254 bp，两端的5.8S和28S分别长84 bp和272 bp(均含引物)。序列比对分析结果表明，5.8S和28S序列高度保守，不适合于种类鉴定，而ITS 2序列高度变异，270个比对位点中有146个位点发生突变，其中72个位点发生插入/缺失突变。除白珠母贝和黑珠母贝之间的遗传距离较小外，其余种类之间的遗传距离远远大于种内遗传距离。基因型分析表明，每个种具有各自特有的基因型。基因型和序列变异分析表明ITS 2序列可作为珍珠贝种类鉴定的分子标记。可用于种间、杂交育种、幼体和珍珠贝肉等材料的种类鉴定与遗传分析。
Abstract:The ITS 2 complete sequences and flanked 5.8S and 28S rDNA partial sequences of six pearl oysters, namely, Pinctada maxima, P.margaritifera, P.albina, P.nigra, P.chemnitzi and P.fucata in the genus Pinctada were analyzed, among which the sequences of P. nigra were obtained from GenBank. The PCR product is about 600 bp, which contains 84 bp of 5.8S and 272 bp of 28S rDNA fragments(both including primer sequence). The full length of ITS 2 ranges from 211 to 254 bp. The sequence variation analysis indicated that 5.8S and 28S rDNA fragments are highly conservative and not suitable for species identification whereas ITS 2 is highly variable with 146 mutation sites (including 72 insertion/deletions) within 270 alignment positions. The interspecific genetic distances are much greater than the intraspecific distances except that between P.albina and P. nigra. The genotype analysis also demonstrated that each species has its own genotypes. These findings suggested that ITS 2 is a high-resolution molecular marker useful in identification of species, hybrids, larvae and tissue samples of the pearl oysters.
-
Keywords:
- species identification /
- genotype /
- sequence variation /
- ribosomal DNA
-
硫酸盐还原菌(sulfate-reducing bacteria，SRB)在生态系统中是属于土著微生物类群，是兼性厌氧菌，广泛存在于各种环境中，是自然界硫循环中SO42-异化还原的专门承担者，将SO42-还原同时降解有机物获取能量，产物为对周围环境有害的有毒气体H2S。目前，国内外对海洋沉积物中SRB研究多集中在海洋工程、石油勘探、管道铺设等金属防腐方面[1-2]。近年来因病害使国内外海水养殖产量出现大幅度的滑坡，其主要原因是养殖环境恶化，尤其是底质的污染严重，大部分养殖海区投喂天然饵料，饵料系数较高，养殖年限较长，网箱底部堆积大量富含有机物的沉积物，为SRB大量繁殖提供了有利条件，SRB繁殖会产生大量有毒气体H2S，造成水体污染并毒害养殖生物，使养殖水体环境状况进入恶性循环状态。目前，国内外对网箱养殖环境中SRB的相关研究报道不多[3]，因此，开展网箱养殖海区沉积物中SRB生态特性和相关防治是很有必要的。本文通过研究大鹏澳网箱养殖海域的SRB，了解此环境中SRB的生态特性、SRB数量变化，以及SRB与其他环境因子之间的关系，为以后修复养殖水体环境，提高养殖生物产量，实现海水网箱养殖业的可持续发展提供理论依据。
大鹏澳是大亚湾西南部的一个小内湾，水面面积约1 400 km2。20世纪90年代以来，海水网箱养殖得到了迅速发展。目前，网箱养殖区水面面积约有20 km2，养殖网箱约4 200箱，养殖区平均水深4.5 m，污染相对比较严重[4-5]。
1. 材料与方法
1.1 样品采集
2006年11月，在大亚湾大鹏澳网箱养殖海域利用GPS定位5个站点，其中1、2、3号站位于鱼类网箱养殖区，4号站位于浮筏贝类养殖区，5号站作为对照站，采样站点如图 1所示。在每个站用柱状采泥器采集柱状沉积物样品2管，将沉积物置于内径为5 cm的PVC管内并用橡胶塞密封好，其中1管用于测定硫化物、氧化还原电位(oxidation-reduction potential，ORP)和pH，另1管用于测定SRB含量，并清晰标记样品垂直方向，其中柱状沉积物SRB计数按照1 cm层、5 cm层、10 cm层深度进行取样。
![]()
图 1 大鹏澳网箱养殖海域采样站位Figure 1. Sampling stations of fish cage aquaculture area in Dapeng′ao Cove1.2 SRB培养基配方和计数方法
培养基配方为K2HPO4 0.5 g · L-1，NH4Cl 1.0 g · L-1，MgSO4 2.0 g · L-1，Na2SO4 0.5 g · L-1，CaCI2 0.1 g · L-1，酵母膏1.0 g · L-1，乳酸钠4 mL ·L-1依次加入净化陈海水中溶解。调节pH值约为7.0~7.5[6]。密封后用蒸气压力灭菌器121~125℃灭菌15~20 min后冷却至室温。称取硫酸亚铁铵1.2 g，抗坏血酸0.4 g，在无菌箱(室)内均匀地摊在离紫外线灯30 cm处灭菌30 min。在无菌操作下，把硫酸亚铁铵和维生素C溶解于事先准备好的40 mL无菌水中。按每100 mL培养基各加入1.0 mL硫酸亚铁铵溶液和1.0 mL维生素C溶液。
硫酸盐还原菌的计数(本标准适用于工业循环冷却水中硫酸盐还原菌的测定，也适用于原水、生活用水及粘泥中硫酸盐还原菌的测定)采用我国国家标准《工业循环冷却水中硫酸盐还原菌的测定》(GB/T14643.5-1993)规定的MPN法[7]。
1.3 沉积物硫化物、ORP、pH的测定
1.3.1 硫化物的测定
沉积物中硫化物的测定采用《海洋监测规范》[8]中的硫离子选择电极法。取5.00 g混匀的湿样置于50 mL烧杯中，加入20 mL抗氧化络合剂使用液，充分搅拌5 min，静置，待沉积物完全沉淀后，用倾斜法将上清液倒入50 mL量瓶中；再重复2次用15 mL抗氧化络合剂使用液按上述步骤浸取残留的沉积物，合并3次上清液，再用抗氧化络合剂使用液定容至刻度。测定时将上清液移入50 mL烧杯中，插入参比电极和活化后的硫离子选择电极，一边搅拌至读数稳定后记录结果。每次测定后均需用蒸馏水洗涤电极并用滤纸吸干，以备下一样品测定。
1.3.2 ORP的测定
沉积物ORP的测定参照《海洋监测规范》[8]中的电位计法。
1.3.3 pH的测定
沉积物pH的测定按照《海洋监测规范》[8]规定采用复合电极方法。
上述分析仪器为上海雷磁PHSJ-4A实验室pH/ORP测定仪，配有硫离子选择电极、参比电极、氧化还原电极、感温电极和pH复合电极。
2. 数据处理
利用Microsoft Excel对SRB的检出率、含量、变化幅度等按照不同站位、不同层次进行分析，利用数理统计软件SPSS对SRB与硫化物、ORP、pH做相关分析。
3. 结果与讨论
3.1 沉积物中的SRB的数量
2006年11月大亚湾大鹏澳网箱养殖海域沉积物中SRB的数量如表 1所示，沉积物中SRB数量波动范围在900~110 000 ind · g-1之间。
表 1 大鹏澳网箱养殖海域沉积物中SRB的数量Table 1. Content of SRB in cage culture area sediment in Dapeng′ao Coveind · g-1 垂直深度/cm
vertical depth1号站
station 12号站
station 23号站
station 34号站
station 4对照站
control station1 45 000 45 000 110 000 25 000 30 000 5 30 000 4 500 45 000 15 000 15 000 10 11 500 2 500 4 500 900 4 500 图 2显示了各个站位表层沉积物SRB的分布情况，SRB平面分布特征是鱼类网箱区＞对照区＞浮筏贝类区，网箱区的数量明显高于对照区和浮筏贝类区，这可能是由于网箱区内的网箱设置比较密集，饵料主要是冰鲜小杂鱼，残饵和粪便形成的生物沉积，为SRB的生长提供了有利的环境；浮筏贝类区SRB数量最少，可能是由于此区域的养殖生物为太平洋牡蛎，无需投饵，故底部富含有机质的沉积物较少，SRB含量也较少。
![]()
图 2 采样站点表层沉积物SRB的平面分布Figure 2. Plane distribution of SRB in surface layer of sediment各个站位沉积物中SRB的垂直分布如图 3所示，可以明显看出SRB数量的分布特征为1 cm层＞5 cm层＞10 cm层，呈现由上到下逐渐减少的趋势。验证了海洋沉积物中硫酸盐还原菌数量是由上向下逐渐降低的结论[9]，这也与高爱国等[10]在北极楚科奇海调查得到的SRB分布规律相似。
![]()
图 3 沉积物中的SRB的垂直分布Figure 3. Vertical distribution of sediment SRB in the five sampling stations3.2 沉积物中的硫化物
2006年11月大亚湾大鹏澳网箱养殖海域5个站位沉积物中硫化物的含量如表 2所示，其硫化物含量波动范围为372.48~703.98 μg · g-1。
表 2 表层沉积物中SRB与其他环境因子Table 2. The amount of SRB and other factors of surface sediment站位
stationSRB数量
content of SRB硫化物含量/μg·g-1
sulfide content氧化还原电位
ORPpH 1号站 station 1 45 000 529.44 -226.5 8.12 2号站 station 2 45 000 527.76 -270.3 8.18 3号站 station 3 110 000 703.98 -384.5 8.16 4号站 station 4 25 000 372.48 -171.5 8.21 对照站 control station 30 000 445.18 -199.5 8.17 由图 4可以明显看出，鱼类网箱养殖区(1、2、3号站)沉积物中硫化物含量要高于贝类区和对照区，其主要原因是鱼类网箱区需要大量投饵，底部富含有机质，这一环境利于SRB大量繁殖，异化还原硫酸盐形成硫化物；浮筏贝类养殖区的硫化物含量最低，这是因为贝类养殖过程中不需投饵，底部沉积物有机质含量较少；另外，该海域贝类区养殖年限较短，水体底部的溶解氧和pH值均较高，不利于SRB的生长，因此，硫化物含量较低。而对照区的5号站，其硫化物含量相对较高于贝类养殖区，这可能是因为5号站位于大鹏澳湾口，在海流作用下，受网箱区沉积物往湾口外迁移的影响要大于贝类区有关。
![]()
图 4 各采样站点沉积物中硫化物的含量Figure 4. Content of sediment sulfide in the five sampling stations3.3 沉积物中的ORP
各采样站位沉积物的氧化还原电位如表 2所示，其波动范围在-171.5~-384.5 mV之间。图 5显示了2006年11月大亚湾大鹏澳网箱养殖海域5个站位沉积物的ORP平面分布特征。可以看出，鱼类网箱区的ORP要低于贝类区和对照区，其主要原因是由于网箱区残饵和排泄粪便长期积累，网箱区底部富含有机质，微生物分解大量有机质，产生还原性的物质(如：S2-、NH4+、Fe2+等)，从而导致氧化还原电位较低。
3.4 沉积物中的pH
图 6显示了2006年11月大亚湾大鹏澳网箱养殖海域各个站位沉积物的pH值，其特征为网箱区低于对照区和贝类区，影响沉积物pH变化的主要因素是呼吸作用[11]，因为网箱区沉积物的有机质含量高于对照区和贝类区，所以网箱区沉积物中的细菌呼吸作用和有机质在细菌参与下的分解作用相对活跃，导致pH下降。
3.5 SRB与其他因子的相关关系
2006年11月大亚湾大鹏澳网箱养殖海域表层沉积物中SRB含量、硫化物含量、ORP、pH值如表 2所示。利用SPSS对SRB与其他环境因子所做的相关分析结果显示，SRB与硫化物的含量为显著性正相关(R=0.96, P < 0.05，表 3)，这与李培英等[12]研究的浙江至闽北陆架积物硫酸盐还原菌与硫化物的相关关系的结果相一致。由于SRB在沉积物环境中主要电子受体是硫酸盐，当SRB大量存在时就会有更多的还原态的硫化物生成。SRB含量与Eh的关系显示极显著性负相关关系(R=-0.97，P < 0.01，表 3)，这与陈皓文等[13]对北部湾东侧沉积物SRB研究结果相一致。ORP的高低主要取决于沉积物中的DO含量，当沉积物处于严重缺氧时，有机物最终分解为还原态物质，Eh为负值，而厌氧的低氧化还原状态环境适合SRB的大量繁殖，所以这一分析结果与理论相一致。而SRB与pH的关系系数为R=-2.71(P>0.05，表 3)，两者不具有显著性相关关系，由此可以看出，硫化物含量和ORP是反映SRB数量的2个非常重要的指标。
表 3 SRB与其他因子的相关系数Table 3. Correlation between SRB and other chemical factors对子 parallelism 相关系数R correlation coefficient n number of sample P significance factor SRB-硫化物 SRB-sulfide 0.96 5 0.011 SRB-氧化还原电位 SRB-ORP -0.97 5 0.006 SRB-pH -2.71 5 0.611 致谢: 本研究承蒙澳大利亚的Dr. Wayne O′Connor提供澳大利亚的珍珠贝样品，广西海洋研究所阎冰先生提供广西北海的样品，广东省珍珠养殖场黄碧光场长提供广东大亚湾的样品，谨此致谢！ -
![]()
图 1 种间ITS 2序列比对分析结果(显示单核苷酸多态(SNP)位点)
Figure 1. The alignment of six species ITS 2 sequences(showing the single nucleotide polymorphic sites).
![]()
图 2 6种珠母贝的UPGMA无根系统发育树
注：图中数字为bootstrap检测值，小于50%的未给出
Figure 2. Unrooted UPGMA tree showing the genetic relationships of the six Pinctada species
Note: The number near the branch is bootstrap test value, value less than 50% are not shown.
表 1 样品种类、采集地点以及ITS 2的扩增片段长度、基因型及其序列号
Table 1 Species, sampling localities as well as amplified fragment length, genotype and accession numbers of ITS 2
种类及采样地点
species and locality代码
code基因型
genotypeITS 2(bp) 总长(bp)
length(bp)GenBank序列号
accession number大珠母贝P.maxima
(海南三亚/ Sanya, Hainan)pmax pmax1 211 525 AY877505 pmax2 211 525 AY883851 pmax3 211 525 AY877504 珠母贝P.margaritifera
(海南三亚/ Sanya, Hainan)pmar pmar1 215 529 AY877507 pmar2 214 528 AY883850 pmar3 215 529 AY877506 白珠母贝P.albina
(澳大利亚/ Port Stephens)palb palb1 251 565 AY877508 palb2 251 565 AY883846 长耳珠母贝P.chemnitzi
(海南三亚/ Sanya, Hainan)pche pche1 251 564 AY877511 pche2 251 564 AY883848 pche3 251 564 AY877510 pche4 252 565 AY877509 pche5 252 565 AY883847 黑珠母贝P.nigra pnig pnig0 254 - AY192714 pnig1 254 - AY282728 pnig2 254 - AY282729 pnig3 254 - AY282730 合浦珠母贝P.fucata
(海南三亚/ Sanya, Hainan(hn)，广东大亚湾/ Daya Bay (db), Guangdong和广西北海/Beihai (bh), Guangxi)db db1 235 548 AY877581 db2 231 544 AY877604 bh bh1 233 546 AY877583 bh2 231 544 AY877605 hn hn1 237 550 AY877592 hn2 230 543 AY877597 
下载: 导出CSV
表 2 种间单核苷酸突变
Table 2 Interspecific single nucleotide mutation
种类
species大珠母贝 珠母贝 白珠母贝 黑珠母贝 长耳珠母贝 合浦珠母贝 大珠母贝 P.maxima - 8 21 21 23 20 珠母贝 P.margaritifera 6 - 25 26 23 20 白珠母贝 P.albina 31 26 - 1 19 18 黑珠母贝 P.nigra 33 31 4 - 19 17 长耳珠母贝 P.chemnitzi 22 22 34 34 - 14 合浦珠母贝 P.fucata 21 19 31 32 24 - 注：对角线下为颠换突变，对角线上为转换突变
Note：transiversion-lower diagonal, transition-upper diagonal
下载: 导出CSV
表 3 珠母贝属种内与种间的遗传距离
Table 3 Intraspecific and interspecific genetic divergences of pearl oysters in Pinctada
种类species 1 2 3 4 5 6 1max 0.010 2mar 0.088 0.011 3alb 0.271 0.280 0.004 4nig 0.287 0.301 0.034 0.014 5che 0.318 0.313 0.096 0.115 0.006 6pfuc 0.268 0.280 0.233 0.238 0.223 0.014 注：对角线为种内遗传距离，对角线下为种间遗传距离
Note：intraspecific-diagonal, interspecific-belon diagonal
下载: 导出CSV
-
[1] 王祯瑞. 中国近海珍珠贝科的研究[J]. 海洋科学集刊, 1978, (14): 101-117. [2] 王祯瑞. 软体动物门双壳纲珍珠贝亚目. 中国动物志, 无脊椎动物, 第31卷[M]. 北京: 科学出版社, 2002.68-98. [3] Hynd J S. A revision of the Australian pearl-shells, genus Pinctada (Lamellibranchia)[J]. Aust J Mar Fresh Res, 1955, 6(1): 98-137. doi: 10.1071/MF9550098
[4] 李刚, 姜卫国, 魏贻尧. 合浦珠母贝、长耳珠母贝和大珠母贝种间人工杂交的研究: Ⅲ同工酶谱的比较研究[J]. 热带海洋, 1983, 2(4): 321-328. https://d.wanfangdata.com.cn/periodical/Ch9QZXJpb2RpY2FsQ0hJTmV3UzIwMjQxMTA1MTcxMzA0EhVDQVMyMDEzMDMwNDAwMDAyMDIwMjkaCDh2ZGF1MTJp [5] Hillis D M, Dixon M T. Ribosomal DNA: Molecular evolution and phylogenetic inference[J]. Q Rev Biol, 1991, 66(4): 411-453. doi: 10.1086/417338
[6] Hillis D M, Davis S K. Ribosomal DNA: Intraspecific polymorphism, concerted evolution, and phylogeny reconstruction[J]. Syst Zool, 1988, 37(1): 63-66. doi: 10.2307/2413191
[7] Sanderson M J, Doyle J J. Reconstruction of organismal and gene phylogenies from data on multigene families: concerted evolution, homoplasy, and confidence[J]. Syst Biol, 1992, 41(1): 4-17. doi: 10.1093/sysbio/41.1.4
[8] Vogler A P, DeSalle R. Evolution and phylogenetic information content of the ITS-1 region in the tiger beetle, Cicindela dorsalis[J]. Mol Biol Evol, 1994, 11(3): 393-405. doi: 10.1093/oxfordjournals.molbev.a040121
[9] Anderson T J, Adlard R D. Nucleotide qequence of a rDNA internal transcribed spacer supports synonymy of Saccostrea commercialis and S. glomerata[J]. J Moll Stud, 1994, 60(2): 196-197. doi: 10.1093/mollus/60.2.196
[10] Miller B R, Crabtree M B, Savage H M. Phylogeny of fourteen Culex mosquito species, including the Culex pipiens complex, inferred from the internal transcribed spacers of ribosoma DNA[J]. Insect Mol Biol, 1996, 5(1): 93-107. doi: 10.1111/j.1365-2583.1996.tb00044.x
[11] Karvonen P, Szmidt A E, Savolainen O. ITS 2 length variation in the internal transcribed spacers of ribosomal DNA in Picea abies and related species[J]. Theor Appl Genet, 1994, 89(6): 969-974.
[12] López-Piňón M J, Insua A, Méndez J. Identification of four scallop species using PCR and restriction analysis of the ribosomal DNA internal transcribed spacer region[J]. Mar Biotech, 2002, 4(5): 495-502. doi: 10.1007/s10126-002-0030-0
[13] Remigio E A, Blair D. Relationships among problematic North American stagnicoline snails (Pulmonata: Lymnaeidae) reinvestigated using nuclear ribosomal DNA internal transcribed spacer sequences[J]. Can J Zool, 1997, 75(11): 1540-1545. doi: 10.1139/z97-779
[14] Beauchamp K A, Powers D A. Sequence variation of the first internal transcribed spacer (ITS-1) of ribosomal DNA in ahermatypic corals from California[J]. Mol Mar Biol Biotech, 1996, 5(3): 357-362. https://pubmed.ncbi.nlm.nih.gov/8983201/
[15] Chen C A, Chen C P, Fan T Y, et al. Nucleotide sequences of ribosomal internal transcribed spacers and their utility in distinguishing closely related Perinereis polychaetes (Annelida; Polychaeta; Nereididae)[J]. Mar Biotech, 2002, 4(1): 17-29. doi: 10.1007/s10126-001-0069-3
[16] Harris D J, Crandall K A. Introgenomic variation within ITS 1 and ITS 2 of freshwater crayfishes (Decapoda: Cambaridae): Implications for phylogenetic and microsatellite studies[J]. Mol Biol Evol, 2000, 17(2): 284-291. doi: 10.1093/oxfordjournals.molbev.a026308
[17] He M, Huang L, Shi J, et al. Variability of ribosomal DNA ITS-2 and its utility in detecting genetic relatedness of pearl oyster[J]. Mar Biotech, 2005(In press). doi: 10.1007/s10126-004-0003-6
[18] 喻达辉, 李有宁, 吴开畅. 中国、日本和澳大利亚珍珠贝的ITS 2序列特征分析[J]. 南方水产, 2005, 1(2): 1-5. doi: 10.3969/j.issn.2095-0780.2005.02.001 [19] Thompson J D, Gibson T J, Plewniak F, et al. The CLUSTAL＿X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucl Acids Res, 1997, 25(24): 4876-4882. doi: 10.1093/nar/25.24.4876
[20] Kumar S, Tamura K, Nei M. MEGA3: molecular evolutionary genetics analysis software[M]. Tempe, Arizona: Arizona State University, USA, 2003.
[21] Rozas J, Sánchez-DelBarrio J C, Messeguer X, et al. DnaSP, DNA polymorphism analyses by the coalescent and other methods[J]. Bioinformatics, 2003, 19(18): 2496-2497. doi: 10.1093/bioinformatics/btg359
[22] 何毛贤, 黄良民. 长耳珠母贝核rRNA基因ITS-2序列分析[J]. 热带海洋学报, 2004, 23(5): 81-84. doi: 10.3969/j.issn.1009-5470.2004.05.011
计量
- 文章访问数: 4981
- HTML全文浏览量: 208
- PDF下载量: 3732
粤公网安备 44010502001741号
