孫麗冉 ，齊迎菊 ，田曉軒

蜈蚣是多足亞門（Myriapoda）唇足綱（Chilopoda）最大最古老的陸生節(jié)肢動(dòng)物，大約有3 300～3 500種蜈蚣分布在除南極洲以外的每個(gè)洲[1]。蜈蚣主要以昆蟲、蚯蚓、蝸牛甚至小脊椎動(dòng)物為食；因?yàn)槿狈ο炠|(zhì)的昆蟲表皮，容易失水，通常棲息于森林的葉、樹皮、木材和土壤中[2]。蜈蚣的第一對(duì)四肢從身體向前伸展覆蓋口，四肢末端有鋒利的爪子和毒腺。毒腺分泌的毒液不僅可以殺死獵物，也可防御天敵，保護(hù)自己的安全。被毒液傷害的人會(huì)有局部疼痛，發(fā)紅，腫脹，表面壞死，寒戰(zhàn)，發(fā)熱，虛弱等癥狀，嚴(yán)重的繼發(fā)感染甚至?xí)?dǎo)致人的死亡[3－6]。

蜈蚣作為傳統(tǒng)中藥，主治半身不遂、破傷風(fēng)、小兒驚風(fēng)、中風(fēng)等病[7]，其藥效物質(zhì)主要是脂肪酸、蛋白質(zhì)、酶及膽甾醇等物質(zhì)[8]。其中，神經(jīng)毒素、離子通道作用的成分及毒液過敏原是哈氏蜈蚣毒液的主要成分。純化的哈氏蜈蚣毒液蛋白或多肽具有血小板聚集活性、抗凝活性磷脂酶A2活性、胰蛋白酶抑制活性、電壓門控鉀通道活性等藥理特性[9]?，F(xiàn)代藥理學(xué)研究發(fā)現(xiàn)，蜈蚣可調(diào)節(jié)機(jī)體脂代謝，降低血脂，保護(hù)心肌細(xì)胞及血管內(nèi)皮細(xì)胞免受脂質(zhì)過氧化及缺血再灌損傷，增加冠脈血流量。同時(shí)，蜈蚣能夠提高人體免疫力，具有改善機(jī)體免疫功能和抗衰老功能。除此以外，蜈蚣水提液可增強(qiáng)腸胃功能，且具有一定的鎮(zhèn)靜、鎮(zhèn)痛作用。體外實(shí)驗(yàn)結(jié)果表明，蜈蚣水提取物對(duì)肝癌、腎癌、肺癌、胃癌、結(jié)腸癌、卵巢癌等多種癌癥細(xì)胞有顯著的生長(zhǎng)抑制作用[10－14]。

中國(guó)供藥用的蜈蚣除藥典收載的少棘巨蜈蚣外，還有哈氏蜈蚣、墨江蜈蚣、多棘蜈蚣等6個(gè)物種在中國(guó)不同地區(qū)作藥用蜈蚣使用。哈氏蜈蚣（Scolopendra subspinipes dehaani）在分類上屬于蜈蚣目（Scolopendromorpha）蜈蚣科（Scolopendridae），是一種體形較長(zhǎng)的（通常為6～20 cm）、毒液毒性較大的蜈蚣[15－16]。哈氏蜈蚣主要分布在廣東、廣西、海南及云南西雙版納等南部熱帶地區(qū)，數(shù)量較少。不同蜈蚣品種間所含藥效物質(zhì)及毒性物質(zhì)不盡相同，對(duì)蜈蚣品種的準(zhǔn)確鑒定可保證用藥安全，而線粒體基因組序列可為其鑒定提供依據(jù)。

動(dòng)物線粒體基因組中含有從分子序列到基因結(jié)構(gòu)的所有信息，其基因的組成排列相對(duì)保守，可以對(duì)許多特征進(jìn)行綜合比較[17]。且與核基因組相比，線粒體基因組相對(duì)較小，易于分析，利于進(jìn)行基于基因組的進(jìn)化與生物類群系統(tǒng)發(fā)生的研究[18]，被廣泛應(yīng)用于動(dòng)物的系統(tǒng)發(fā)育以及種群保護(hù)等領(lǐng)域的研究[19]。本文采用Illumina HiSeq X Ten測(cè)序技術(shù)，獲得了哈氏蜈蚣線粒體基因組的完整序列，是國(guó)內(nèi)外首個(gè)哈氏蜈蚣線粒體基因組的公開報(bào)道。

1 實(shí)驗(yàn)材料

1.1 實(shí)驗(yàn)樣品 哈氏蜈蚣樣品購(gòu)于廣東佛山，并經(jīng)形態(tài)學(xué)鑒定和COI基因分子鑒定確定購(gòu)買樣品為該物種。Voucher ID為hlwg－1。

1.2 實(shí)驗(yàn)方法

1.2.1 線粒體DNA富集 取新鮮肌肉組織立即使用。利用差速離心的方法富集線粒體細(xì)胞器，然后采用酚氯仿抽提法，獲得哈氏蜈蚣總基因組DNA，并使用Nanodrop 2000粗略測(cè)定DNA濃度及純度,然后利用 REPLI－g（R）Mitochondrial DNA Kit試劑盒對(duì)總DNA進(jìn)行滾環(huán)擴(kuò)增，提升線粒體DNA相對(duì)于基因組DNA的比例，得到更多線粒體DNA。

1.2.2 Illumina HiSeq X Ten高通量測(cè)序

1.2.2.1 建庫(kù) 利用超聲波將滾環(huán)擴(kuò)增產(chǎn)物打斷為特定長(zhǎng)度的DNA片段，并在這些片段兩端連接上序列已知的通用接頭構(gòu)建出文庫(kù)。

1.2.2.2 Flowcell 構(gòu)建好的文庫(kù)通過flowcell時(shí)會(huì)隨機(jī)附著在flowcell表面的channel上，每個(gè)channel表面有很多和建庫(kù)過程中連接在DNA片段兩端的接頭互補(bǔ)配對(duì)的接頭。

1.2.2.3 橋式擴(kuò)增 橋式PCR以flowcell表面固定的DNA鏈為模板，進(jìn)行橋式擴(kuò)增，經(jīng)過變性、退火、延伸的不斷循環(huán)，最終每個(gè)DNA片段都在各自的位置上形成簇。

1.2.2.4 上機(jī)測(cè)序 向反應(yīng)體系中同時(shí)添加包含DNA聚合酶、接頭引物及帶有保護(hù)基團(tuán)和特異熒光的4種dNTP的混合物。在測(cè)序過程中，當(dāng)dNTP被結(jié)合到合成鏈上之后，所有其他未使用的dNTP及DNA聚合酶都被沖洗掉，再加入所需的緩沖液，激發(fā)熒光，記錄熒光信號(hào)并轉(zhuǎn)化為堿基信號(hào)，然后加入化學(xué)試劑淬滅熒光基團(tuán)及去掉保護(hù)基團(tuán)以便進(jìn)行下一輪的測(cè)序反應(yīng)。

1.2.3 序列拼接 將測(cè)序的原始數(shù)據(jù)從測(cè)序服務(wù)器上下載至本地，將低質(zhì)量序列、過短序列以及錯(cuò)誤序列過濾掉得到用于下游分析的序列，根據(jù)目的基因組大小，調(diào)整合適的覆蓋倍數(shù)，使用MIRA 4.0等軟件進(jìn)行de novo拼接，得到片段較長(zhǎng)的contigs，再結(jié)合近緣物種的序列進(jìn)行mapping組裝，取得質(zhì)量較好的consensus替代近源物種序列，繼續(xù)進(jìn)行mapping組裝，經(jīng)反復(fù)迭代后，進(jìn)行最高質(zhì)量堿基合意，合意后得到完整線粒體全基因組序列。

1.2.4 注釋分析和提交 通過與近緣種序列進(jìn)行比較以及MITOS server線粒體在線注釋系統(tǒng)[20]，對(duì)全序列進(jìn)行注釋。以盲蜈蚣科Scolopocryptops sp.（GenBank KC200076）的線粒體基因組全序列為參照，對(duì)蛋白質(zhì)編碼基因、tRNA、rRNA及控制區(qū)等進(jìn)行注釋，tRNA二級(jí)結(jié)構(gòu)由MITOS[12]折疊生成。使用軟件MEGA 4.0及Bioedit等統(tǒng)計(jì)哈氏蜈蚣線粒體基因組的序列總長(zhǎng)、堿基組成以及GC含量等。完整的哈氏蜈蚣線粒體基因組使用Sequin對(duì)各基因進(jìn)行詳細(xì)注釋后，提交到GenBank，登錄號(hào)為KY947341。

圖1 哈氏蜈蚣線粒體基因組結(jié)構(gòu)圖Fig.1 The mitochondrial genome structure of Scolopendra subspinipes dehaani

2 實(shí)驗(yàn)結(jié)果

2.1 序列拼接結(jié)果 組裝得到哈氏蜈蚣線粒體基因組后，經(jīng)bowtie2比對(duì)[21]，確認(rèn)共2 557 168條序列（reads）被用于拼接該基因組，這些序列的平均長(zhǎng)度為150 bp，測(cè)序深度達(dá)39 293倍。

2.2 整體結(jié)構(gòu) 哈氏蜈蚣線粒體基因組結(jié)構(gòu)見于圖1與表1。哈氏蜈蚣共編碼35個(gè)基因，包括了20個(gè)tRNA，13個(gè)蛋白編碼基因，2個(gè) rRNA（16SrRNA、12SrRNA）并有1個(gè)控制區(qū)。其中除ND5、ND1、ND4、ND4L、16SrRNA、12SrRNA 和 8 個(gè) tRNA（tRNA－Leu、tRNA－Val、tRNA－Gln、tRNA－Pro、tRNA－Cys、tRNAPhe、tRNA－Tyr和 tRNA－His）由輕鏈（L 鏈）編碼，其余21個(gè)基因均由重鏈（H鏈）編碼。

2.3 堿基組成 哈氏蜈蚣線粒體基因組全長(zhǎng)14 538 bp，其中H鏈的堿基含量為：A=6 085 bp（41.9%），C=2 519 bp（17.3%），G=1 248 bp（8.6%），T=4 686 bp（32.2%），GC=3767 bp（25.9%），A＞T＞C＞G，AT%＞GC%。

2.4 蛋白編碼基因 哈氏蜈蚣線粒體基因組中蛋白編碼基因全10 810 bp，約占全基因組的74.36%。A=4 542 bp（42.0%），C=1 924 bp（17.8%），G=958 bp

（8.9%），T=3 386 bp（31.3%），GC=2 882 bp（26.7%），A＞T＞C＞G，AT%＞GC%。所有蛋白編碼區(qū)均無內(nèi)含子。除ND1以TTG為起始密碼子，ND5和ND4以ATT為起始密碼子，CO1以ACG為起始密碼子,ND2、ATP8、ND3、ND4L、ND6 以 ATA 起始外，其余均使用ATG為起始密碼子。除ND3以TAG為終止密碼子，ND2、COX2、ND1和 ATP8 以 TAA 為終止密碼子，其余蛋白均以T堿基終止，并通過轉(zhuǎn)錄過程中3’端的poly A序列補(bǔ)全，形成完整的TAA終止密碼子。氨基酸使用頻率的分析見于表2。在3 598 個(gè)氨基酸中，氨基酸 Cys（1.1%）、Gln（1.8%）和Arg（1.6%）出現(xiàn)頻率較低，而頻率最高的為L(zhǎng)eu（14.0%）。

表1 哈氏蜈蚣線粒體基因組表Tab.1 The table of Scolopendra subspinipes dehaani mitochondrial genome

表2 哈氏蜈蚣線粒體基因組蛋白編碼基因氨基酸組成分析Tab.2 The amino acid composition in the protein coding genes on Scolopendra subspinipes dehaani mitochondrial genome

2.5 rRNA和tRNA基因 哈氏蜈蚣的12S rRNA和16S rRNA長(zhǎng)度分別為765 bp和1 220 bp，GC含量分別為25.2%和24.9%。tRNA由MITOS[12]注釋生成。20個(gè)tRNA基因中，最長(zhǎng)的為 tRNA－Ala（anticodon=TGC，80 bp）， 最 短 的 為 tRNA －Asp（anticodon=GTC，42 bp）。

2.6 非編碼區(qū) 哈氏蜈蚣線粒體基因組的控制區(qū)（control region）位于tRNA－Ile與12S rRNA之間，為732 bp。除控制區(qū)外，哈氏蜈蚣線粒體基因組中有6個(gè)非編碼間隔區(qū)域，其中最長(zhǎng)的間隔區(qū)長(zhǎng)46 bp，位于12S rRNA和16S rRNA之間。

3 討論

本實(shí)驗(yàn)中采用差速離心以及滾環(huán)擴(kuò)增的方法，得到純度較高的線粒體DNA，打斷后進(jìn)行測(cè)序。用這種方法明顯提升測(cè)序樣品中線粒體DNA的比例，降低了細(xì)胞核基因的影響。相比于從總DNA測(cè)序結(jié)果中篩選拼接線粒體的方法[22－24]，該方法對(duì)于線粒體DNA的測(cè)序深度顯著增加，能夠盡可能降低測(cè)序錯(cuò)誤及假陽(yáng)性結(jié)果的發(fā)生，得到更準(zhǔn)確的線粒體基因組信息。

線粒體基因組具母系遺傳特征，基因結(jié)構(gòu)相對(duì)保守，其蛋白編碼區(qū)及tRNA區(qū)結(jié)構(gòu)變化在系統(tǒng)發(fā)育研究中應(yīng)用較廣[19]。本研究中僅發(fā)現(xiàn)20個(gè)tRNA,與蜈蚣目已發(fā)表的其他物種相比[25]，分別在tRNASer與 tRNA－Phe之間缺少一個(gè) tRNA－Glu、在ND1與tRNA－Leu之間缺少一個(gè)tRNA－Leu。tRNA的缺失可能是由于tRNA的截?cái)嗷蚱涠?jí)結(jié)構(gòu)不對(duì)稱造成的[26]。這種現(xiàn)象的出現(xiàn)有多種可能的原因，如：缺少完全配對(duì)的受體莖、tRNA區(qū)與蛋白編碼區(qū)或tRNA區(qū)的重疊、不對(duì)稱的tRNA序列無法折疊成標(biāo)準(zhǔn)的三葉草結(jié)構(gòu)等。在許多后生動(dòng)物的線粒體基因組中都有類似現(xiàn)象出現(xiàn)[27－28]。

目前，蜈蚣目中僅有一個(gè)物種已發(fā)表線粒體基因組全序列，其近源物種的線粒體基因組研究未見報(bào)道。對(duì)哈氏蜈蚣線粒體基因組的研究，填補(bǔ)了蜈蚣科線粒體基因組數(shù)據(jù)的空白。未來更多近源物種數(shù)據(jù)的發(fā)表，將為藥用蜈蚣的資源普查及保護(hù)提供可靠的依據(jù)。

[1] Undheim EA,King GF.On the venom system of centipedes(Chilopoda),a neglected group of venomous animals[J].Toxicon,2011,57(4):512－524.

[2] Edgecombe G,Giribet G.Evolutionary biology of centipedes(Myriapoda:Chilopoda)[J].Annual Review of Entomology,2007,52(1):151－170.

[3] Uzel AP,Steinmann G,Bertino R,et al.Necrotizing fasciitis and cellulitis of the upper limb resulting from centipede bite:two case reports[J].Chirurgie De La Main,2009,28(5):322－325.

[4] Knysak I,Martins R,Bertim C R.Epidemiological aspects of centipede (Scolopendromorphae:Chilopoda)bites registered in greater S.Paulo,SP,Brazil[J].Revista De Saúde Pública,1998,32(6):514.

[5] Mohri S,Sugiyama A,Saito K,et al.Centipede bites in Japan[J].Cutis,1991,47(3):189－190.

[6]Fung HT,Lam SK,Wong OF.Centipede bite victims:a review of patients presenting to two emergency departments in Hong Kong[J].Hong Kong medical journal,2011,17(5):381－385.

[7] 方秀桐,莫可元.蜈蚣的藥理研究進(jìn)展[J].中國(guó)醫(yī)藥指南,2015,13(18):32－34.

[8] 梁潔,黃華營(yíng).藥用蜈蚣的化學(xué)成分及藥理活性研究近況[J].廣西中醫(yī)藥,2005,28(4):6－7.

[9] Liu ZC,Zhang R,Zhao F,et al.Venomic and transcriptomic analysis of centipede Scolopendra subspinipes dehaani[J].Journal of Proteome Research,2012,11(12):6197－6212.

[10]陳武忠,曹海燕.中藥蜈蚣的研究進(jìn)展[J].中國(guó)現(xiàn)代中藥,2011,13(7):54－57.

[11]黃曉蒂,陳 園,田 莎,等.蜈蚣提取物對(duì)人肝癌細(xì)胞甲基化譜的影響[A].第八次全國(guó)中西醫(yī)結(jié)合中青年學(xué)術(shù)論壇論文集[C].北京:2016.

[12]姜璐璐,曹艷華,王詩(shī)敏,等.蜈蚣毒素的研究進(jìn)展[J].西北藥學(xué)雜志,2009,24(6):517－520.

[13]羅 雷,楊仕隆,賴 仞.蜈蚣毒素研究進(jìn)展[J].生命科學(xué),2016(1):27－32.

[14]周永芹,韓 莉.中藥蜈蚣的研究進(jìn)展[J].中藥材,2008,(02):315－319.

[15]康四和,鄧海英,江珍玉,等.我國(guó)藥用蜈蚣分類鑒定及資源研究[J].中藥材,2016,39(4):727－731.

[16]康四和,鄧海英,劉義梅,等.我國(guó)藥用蜈蚣資源與品質(zhì)研究現(xiàn)狀及對(duì)策[J].中南民族大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,35(2):61－66.

[17]養(yǎng)小蟲賺大錢本刊讀者服務(wù)部舉辦中華土元、紅龍蜈蚣、黃粉蟲與竹鼠養(yǎng)殖技術(shù)培訓(xùn)班[J].農(nóng)村百事通,2013,32(11):84-85.

[18]Liu ZC,Zhang R,Zhao F,et al.Venomic and transcriptomic analysis of centipede Scolopendra subspinipes dehaani[J].Journal of Proteome Research,2012,11(12):6197-6212.

[19]左寶平.少棘蜈蚣線粒體基因組全序列和平耳孔蜈蚣的地理種群遺傳變異分析[D].南京：南京師范大學(xué),2007.

[20]龔美蓉.多足動(dòng)物的系統(tǒng)發(fā)生關(guān)系及少棘蜈蚣線粒體基因組研究[D].南京：南京師范大學(xué),2003.

[21]劉 杰,田曉軒,崔 英,等.藥用動(dòng)物烏梢蛇線粒體基因組全序列分析[J].天津中醫(yī)藥大學(xué)學(xué)報(bào),2016,35(3):187-191.

[22]Bernt M,Al E.MITOS:improved de novo metazoan mitochondrial genome annotation.[J].Molecular Phylogenetics&Evolution,2013,69(2):313-319.

[23]Langmead B,Salzberg SL.Fast gapped-read alignment with Bowtie 2[J].Nature Methods,2012,9(4):357-359.

[24]Hahn C,Bachmann L,Chevreux B.Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—a baiting and iterative mapping approach[J].Nucleic Acids Research,2013,41(13):e129.

[25]Gillett CPDT,Cramptonplatt A,Timmermans MJTN,et al.Bulk De Novo Mitogenome Assembly from Pooled Total DNA Elucidates the Phylogeny of Weevils(Coleoptera:Curculionoidea)[J].Molecular Biology&Evolution,2014,31(8):2223-2237.

[26]Tang M,Tan M,Meng G,et al.Multiplex sequencing of pooled mitochondrial genomes—a crucial step toward biodiversity analysis using mito-metagenomics[J].Nucleic Acids Research,2014,42(22):e166.

[27]Gai Y,Ma H,Ma J,et al.The complete mitochondrial genome of Scolopocryptops sp.(Chilopoda:Scolopendromorpha:Scolopocryptopidae)[J].Mitochondrial Dna,2014,25(3):192-193.

[28]Robertson HE,Lapraz F,Rhodes AC,et al.The Complete Mitochondrial Genome of the Geophilomorph Centipede Strigamia maritima[J].Plos One,2015,10(3):e0121369.

[29]Lavrov DV,Brown WM,Boore JL.A novel type of RNA editing occurs in the mitochondrial tRNAs of the centipede Lithobius forficatus[J].Proceedings of the National Academy of Sciences of the United States of America,2000,97(25):13738-13742.

[30]Masta SE,Boore JL.The complete mitochondrial genome sequence of the spider Habronattus oregonensis reveals rearranged and extremely truncated tRNAs[J].Molecular Biology&Evolution,2004,21(5):893-902.