華北地區中華蒙潮蟲(甲殼動物亞門：等足目)種群遺傳多樣性研究

張瑞， 王明曉， 安建梅

(山西師范大學生命科學學院，山西臨汾041000)

中華蒙潮蟲Mongoloniscussinensis(Dollfus，1901)隸屬于節肢動物門Arthropoda甲殼動物亞門Crustacea真軟甲亞綱Eumalacostrace軟甲綱Malacacostra囊蝦總目Peracarida等足目Isopoda潮蟲亞目Oniscidea緣潮蟲科Agnaridae蒙潮蟲屬Mongoloniscus(Boykoetal.，2008)，中國特有種。中華蒙潮蟲的環境適應能力非常強，主要分布于山西、北京、內蒙古、吉林等地(陳國孝，2000)，是我國古北界典型的溫帶特有種，本課題組在山東、陜西西安、遼寧、天津、河南、河北均發現該物種。師二燕(2015)基于線粒體細胞色素c氧化酶亞基Ⅰ基因(mtCOⅠ)對中華蒙潮蟲進行了譜系地理研究，共基于26個居群(跨越5省15市26縣區)進行了分析，推測該物種可能的進化路線為從北往南，且最有可能的發源地為我國東北平原。有研究表明，中華蒙潮蟲還可作為環境的重金屬指示物種(牛曉倩，2015；李夢雯，2016)。關于多基因聯合分析中華蒙潮蟲種群遺傳多樣性的相關研究未見報道。

分子標記以蛋白質、核酸的突變為基礎，檢測生物遺傳結構及其變異(白玉，2007)。動物線粒體基因組遵循嚴格的母系遺傳且分子的各部分均共享同一祖先(Wilsonetal.，2010)，且線粒體基因組進化速率快，是單拷貝核DNA的5～10倍，群體內變異大，群體遺傳進化常用此標記(陳復生等，2003；Itoetal.，2011；Zhouetal.，2011；Lietal.，2012)，也是物種生物地理研究的常用標記之一(Aviseetal.，1987；徐慶剛，花保禎，2001；王蘭萍等，2013)。

線粒體內不同基因的進化速率不同。NADH5氧化還原酶基因(ND5)進化速率相當快，是種群系統發育研究中較有效的基因之一(智妍等，2008)。COⅠ基因常用于種以下的系統發育研究，具有進化速率快、富含系統發育遺傳信息等特點，常用于探討物種分類及近緣種和種群遺傳結構多樣性等系統進化關系(鄧春興，2014；劉平等，2018)。本研究采用基因聯合(COⅠ+ND5)的方法對中國華北地區10個地理種群89只中華蒙潮蟲進行群體遺傳多樣性分析，以期了解華北地區中華蒙潮蟲種群間的遺傳分化和基因交流。

1 材料和方法

1.1 樣品采集

2011年4月—2018年6月采集了10個地理種群共89只中華蒙潮蟲(表1，圖1)，標本浸泡于無水乙醇。

外群選用海蟑螂Ligiaoceanica(Linnaeus，1767)(GenBank登錄號：DQ442914)和Cylisticusconvexus(De Geer，1778)(GenBank登錄號：KR013002)。

表1 中華蒙潮蟲采集信息Table 1 Sampling informations of Mongoloniscus sinensis

圖1 中華蒙潮蟲采樣點Fig. 1 The sample locations of Mongoloniscus sinensis

1.2 DNA的提取及測序

DNA的提取使用Chelex-100法(戴文申等，2007)，選取綠豆大小的中華蒙潮蟲胸足或腹足肌肉于0.5 mL離心管中，加純水清洗，隨后放入干凈的培養皿中晾干。將晾干的組織樣品放入200 μL 5%的Chelex-100離心管中，并加入5 μL 5 mg·mL-1的蛋白激酶，混勻振蕩，確保肌肉組織完全浸沒于Chelex-100的顆粒中。56 ℃ 3 h以上，取出振蕩；100 ℃ 8 min，振蕩；13 000 r·min-1離心3 min，取上清用于PCR擴增或-20 ℃保存備用。

COⅠ基因引物(Folmeretal.，1994)為COⅠF：5’-GGTCAACAAATCATAAAGATATTGG-3’；COⅠR：5’-TAAACTTCAGGGTGACCAAAAAATCA-3’。ND5基因引物為877F：5’-TTTATCTTTTGGGTTCGCTA-3’；1608R：5’-TAAAATTAAATCCTTGCCCTC-3’，是選取中華蒙潮蟲ND5基因全序列(GenBank登錄號：MG729627)的877～1 608 bp片段(編碼區序列)，用Oligo 7(Rychlik，2007)加上手工調整而設計的。引物合成委托蘇州金唯智生物科技有限公司完成。PCR的擴增體系為50 μL：2×HieffTMPCR Master Mix 25 μL，上、下游引物各3 μL，DNA模板7 μL，ddH2O 12 μL。

PCR在TCA0096擴增儀上進行：94 ℃ 2 min；94 ℃ 45 s，退火45 s(退火溫度為COⅠ基因：50 ℃；ND5基因：45.5 ℃)，68 ℃ 1 min，35個循環；72 ℃ 10 min。擴增產物用0.7%瓊脂糖凝膠進行電泳檢測，擴增效果良好且足量的樣品原液委托蘇州金唯智生物科技有限公司進行純化及雙向測序。

1.3 數據分析

測序所得序列峰圖用Geneious(Kearseetal.，2012)檢查和篩選，經校驗后進行序列雙向拼接，剪切掉引物部分，并分別進行BLAST搜索以確定每個基因片段的同源性。所有序列在Geneious中排序比對。序列拼接使用SequenceMatrix，將2個基因序列合并為1個完整序列數據集，并命名為COⅠ+ND5。利用PAUP 4.0對序列集的堿基組進行同質性檢驗。使用MEGA 7(Kumaretal.，2016)對獲得的序列數據集計算堿基組成、保守位點、變異位點以及序列間的轉換顛換比率；采用Kimura(Kimura，1980)雙參數模型計算種群遺傳距離。使用DnaSP 5.0(Librado & Rozas，2009)統計單倍型，計算各種群核苷酸多樣性(Pi)、單倍型多樣性(Hd)，并構建遺傳分化指數(Fst)，Fst反映2個種群之間的遺傳差異程度，其值與基因流(Nm)成反比，即Fst越大，基因交流程度越小；中性檢驗通過計算Tajima’sD值和Fu’sFs值來表示種群是否經歷過擴張，二者值接近0時，種群較為穩定，小于0(P<0.05)時，種群近期經歷過擴張，大于0時，種群可能出現了分化(Fu，1997)。構建單倍型岐點分布圖，通過可視化曲線觀察種群近期是否經歷過擴張。系統發育分析分別采用MEGA 7構建最大簡約(MP)樹和MrBayes 3.0(Ronquistetal.，2012)構建貝葉斯(BI)樹。BI樹中，由于每個基因、編碼基因中每個密碼子位點進化速率的不同，采用Partitionfinder 1.1.1(Lanfearetal.，2012)分別對COⅠ與ND5基因的第一、第二、第三密碼子進行分區計算。

采用Arlequin 3.5(Excoffer & Lischer，2010)進行分子變異分析(AMOVA)。但所有單倍型幾乎未按地理來源分支，因此按照系統發育劃分出支系進行區域分組，將分布零散的運城獨立為1個組。使用Network 4.1構建單倍型網絡圖。

2 結果與分析

2.1 序列特征

中華蒙潮蟲COⅠ部分基因長604 bp，ND5部分基因長615 bp，拼接后長度為1 219 bp。對數據集(COⅠ+ND5)進行統計分析，發現保守位點716個，變異位點503個(占總序列長度41.3%)，其中，單一多態位點235個，簡約信息位點268個。聯合序列的堿基A、T、C、G平均含量分別為30.8%、41.0%、11.2%、17.0%，A+T(71.8%)明顯高于C+G(28.2%)，符合節肢動物門高A+T含量特點(孫紅英等，2003)，序列間的轉換/顛換比值為2.8。

數據集(COⅠ+ND5)同質化檢驗結果顯示，P=0.053，表明二者為同質性樣品，可進行聯合分析。

2.2 單倍型多樣性和遺傳多樣性

89只中華蒙潮蟲線粒體聯合基因序列(COⅠ+ND5)共45種單倍型(表2)。其中，單倍型H1[北京東城區(DCQ)、河北石家莊(SJZ)、遼寧葫蘆島(HLD)]、H16[河南新鄉(XX)、山西臨汾(LF)、山西呂梁(LL)] 和H21[SJZ、LL、山東桃村(TC)]為3個種群間的共享單倍型；H15[LF、山西運城(YC)]、H41[XX、陜西西安未央區(WYQ)]為2個種群間的共享單倍型，其余為種群獨享單倍型，表現出地理種群之間明顯的遺傳分化。10個種群總的Hd為0.964(DCQ僅1種單倍型，無種群內多樣性體現)。其余9個地理種群中，SJZ、HLD、LL、山西大同(DT)、YC顯示出較高的單倍型多樣性(Hd>0.900)，但總體核苷酸多樣性(0.005 6)較低。

2.3 種群間的遺傳距離和遺傳分化

WYQ、LF與XX之間的遺傳距離最小(0.00～0.01)，這3個種群與其他種群的遺傳距離隨著地理位置增大，越往北延伸遺傳距離越大。DT和其他9個種群的遺傳距離均較大(0.08～0.11)，DCQ與HLD種群的遺傳距離為0。所有種群的平均遺傳距離為0.06。

Fst<0.150的種群主要有：TC與HLD、LL與SJZ、DCQ與HLD、YC與LL和SJZ，以及SJZ與LL。LF、WYQ和XX這3個種群分別與DCQ和HLD種群的遺傳分化系數較大(Fst>0.900)，基因交流水平貧乏，絕大多數種群之間的遺傳分化系數為0.150～0.500(表3)。平均Fst為0.513，平均Nm為0.24。

2.4 種群動態分析

僅TC種群顯示出明顯的種群擴張(P<0.05)，XX和HLD種群的Tajima’sD雖為負值，但均不顯著，未達到種群擴張(P>0.10)。SJZ種群的Tajima’sD值為不顯著負值(0.05

表2 中華蒙潮蟲(COⅠ+ND5)聯合基因單倍型分布Table 2 Haplotypes of mt DNA (COⅠ+ND5) genes of Mongoloniscus sinensis

表3 中華蒙潮蟲(COⅠ+ND5)序列的種群遺傳距離(對角線上方)和各種群間遺傳分化系數(對角線下方)Table 3 Genetic distance (above diagonal) and pairwise fixation indices of genetic variation (below diagonal) of Mongoloniscus sinensis among different populations based on COⅠ and ND5 genes

單倍型歧點分析結果見圖2。基于10個種群的所有單倍型的岐點分布曲線未成單一峰型，表明中華蒙潮蟲近期未曾經歷種群擴張事件，即種群之間存在一定程度的分化。這與種群之間較大的遺傳分化系數吻合。

2.5 分子系統發育分析

Partitionfinder計算所得COⅠ+ND5聯合序列的最佳分區方案和各區最適核苷酸替代模型見表4。

中華蒙潮蟲不同地理種群的單倍型系統進化關系幾乎得到了完全一致的系統發育樹，僅節點支持率略有差異，這里僅顯示BI樹(圖3)。中華蒙潮蟲個別地理種群的單倍型并沒有完全按照地理來源形成明顯的譜系地理結構。單倍型較好的DCQ、HLD、DT、WYQ、XX以及LF種群按地理分布進行劃分，其余種群如YC、SJZ、LL以及TC種群則出現了單倍型混雜的分布情況。按照進化的先后以及不同地理種群的單倍型聚合情況大致可以分為5個主要支系；C1為DT；隨之進化出的主要類群為TC和SJZ，劃分為C2；C3主要為DCQ和HLD；LL雖分布零散，但絕大部分個體位于C3之后，因此劃分為C4；最后演化出的是WYQ、LF及XX種群，劃分為C5；YC單倍型個體在各個支系均有涉及，主要分布于C4與C5。

圖2 所有單倍型總體錯配分布Fig. 2 Mismatch distribution analysis for all haplotypes

2.6 分子變異分析

用Arlequin對劃分后的5個組進行分子變異分析，中華蒙潮蟲的分子變異在各組間以及種群內差異不大，種群內略高(53.12%)(表5)。

表5 基于COⅠ與ND5組合的組間和群體間的變異分析Table 5 AMOVA analysis for Mongoloniscus sinensis among groups and populations based on COⅠ and ND5 genes

2.7 單倍型網絡圖分析

中華蒙潮蟲各個地理種群的單倍型幾乎未按地理來源進行劃分，這與分子系統樹所得結果一致。單倍型H1、H21及H16出現的頻率較高且屬于2～3個種群的共享單倍型，推測其為祖先類型。LL、SJZ種群較為分散，與系統樹類型一致；WYQ、LF及XX種群明顯聚成一個簇群。

3 討論

遺傳多樣性是生態系統多樣性、物種多樣性和景觀多樣性的基礎，通常遺傳多樣性最直接的表現形式就是決定進化潛力的遺傳變異水平。對于任何一個物種而言，個體生命雖短暫，但種群或種群系統在自然界卻具有特定的分布格局，這些分布格局會隨著外界條件不斷發生改變，因此，遺傳多樣性不僅包括變異水平的高低，而且包括變異的分布格局(Shen & Liu，2001)。

中華蒙潮蟲屬于生物圈中分布極廣的節肢動物門，雖不如昆蟲綱Insecta種類繁多，但也具備極強的繁殖生長能力，種群數量相對較大。該物種的研究報道十分稀少，遺傳多樣性的研究也非常薄弱。本研究首次聯合使用COⅠ及ND5基因序列對我國華北地區中華蒙潮蟲種群遺傳多樣性進行了探討。

圖3 基于COⅠ與ND5基因組合序列構建的中華蒙潮蟲貝葉斯樹Fig. 3 Bayesian tree of Mongoloniscus sinensis based on COⅠ and ND5 genes

圖4 中華蒙潮蟲(COⅠ+ND5)單倍型網絡圖Fig. 4 Haplotype network of Mongoloniscus sinensis based on COⅠ and ND5 genes

3.1 種群遺傳多樣性

衡量一個種群遺傳多樣性水平可通過Hd及Pi綜合體現。本研究中，DCQ種群僅1種單倍型，該地區的遺傳分化極為貧乏，結合目前城市化水平的升高導致物種棲息地破壞、生境破碎、基因交流斷裂，加之北京西部是太行山余脈的西山，北部是燕山山脈的軍都山，兩山在南口關溝相交形成的大山彎，山脈的阻隔降低種群交流的程度。SJZ種群的Hd(0.964)和Pi(0.058 0)均顯示出較高的遺傳多樣性；HLD、LL、DT和TC種群均具有較高的Hd(>0.900)，但Pi整體不高。一般認為，如果1個種群有較高的Hd，但核苷酸水平較低，推測可能是種群的建群者效應，即1個較小的有效種群通過內部的變異，積累了大量單倍型多態性，但卻未達到積累核苷酸序列多樣性的水平，由此種群迅速增長(Wright，1943)。所有種群顯示出較高的Hd以及較低的Pi，據此推斷華北地區中華蒙潮蟲種群遺傳多樣性水平中等。

3.2 種群遺傳結構

種群遺傳結構的差異是遺傳多樣性的重要體現，種群遺傳變異決定該種群內物種進化潛力以及抵御不良環境的能力(Grant，1991)。單倍型岐點圖以及中性檢驗顯示該種群近期未經歷擴張，但總體有平穩分化現象。基于AMOVA顯示，中華蒙潮蟲的變異主要來自種群內部，Fst(0.513)顯示各個種群之間遺傳分化程度明顯，平均Nm(0.24)很低，考慮中華蒙潮蟲是陸地爬行類物種，因此，生境限制其無法遠距離活動，造成種群之間交流稀缺，這也是中華蒙潮蟲種群內部不斷分化的原因。

3.3 系統發育分析

中華蒙潮蟲作為中國古北界特有種，BI樹和單倍型網絡圖顯示，不同種群之間的單倍型并未完全按地理來源進行分布。位于系統發育起始端主要為華北偏北地區的種群，如DT、SJZ；而最后演化出來的為WYQ、LF和XX的全部個體，這也在一定程度上與師二燕(2015)的研究結果相吻合，即中華蒙潮蟲的主要演化路線為從北向南，同理，與冰川影響下物種的遷移路線基本吻合。另外，TC及YC的單倍型個體較為分散。Hap36(TC)是除外群之外最早分化出的一種單倍型，考慮山東位于黃河下游，東臨渤海黃海，且中華蒙潮蟲的進化未經歷淡水階段，直接從海洋過渡到陸生(Brolyetal.，2013)，可以將Hap36看做中華蒙潮蟲的原始單倍型；其次Hap39/37/38(TC)與緯度較為接近的SJZ主要支系交叉分布；基于中性檢驗分析，TC數據顯著，近期可能經歷了擴張。YC樣本量雖偏少，但結合系統發育樹以及單倍型分布情況，主要劃分為C4與C5支系，其地理位置與演化情況幾乎和研究結果相一致，推測數據對路徑的演化影響不大。

本研究中，由于個別種群采集數量較少，單倍型分布零散，考慮小的樣本量可能造成實驗結果有一定偏差，后期將繼續加大樣本量，完善數據信息。近年來分子標記作為生物學領域一個強有力的工具，在遺傳多樣性以及生態領域均發揮重要作用，然而，當種群樣本數很少或不平均的情況下，利用分子標記對遺傳多樣性和遺傳分化進行解釋就需要謹慎。樣本數太大會花費很多時間、資源和經費，而樣本數太小往往會導致錯誤的結論(王麗等，2010)，因此，樣本量是一個非常重要的問題。后期應加大采集量，轉向南方未涉及到的區域，加大種群覆蓋地區，細化中華蒙潮蟲的系統演化路徑以及種群分布的動態情況。

4 結語

土壤中生活著豐富的生物類群，它們在自然生態系統中扮演著消費者和分解者的角色，是重要的地下生物寶庫，對全球物質循環和能量流動起著不可替代的作用。然而，由于人類活動的強烈干擾，生境的不斷喪失，土壤生物多樣性已日趨減少，這又加劇了土地的退化。由于土壤生物與地上動植物息息相關，只有將對二者的保護結合起來，才能實現完整的全球意義上的生物多樣性保護(章家恩，1999)。本研究基于線粒體COⅠ和ND5基因的聯合序列，對華北地區中華蒙潮蟲遺傳多樣性進行了研究，發現其遺傳多樣性處于中等水平，主要變異來自于種群內部，個別種群的單倍型沒有完全按照地理分布形成簇群；系統發育和單倍型網絡圖顯示，基于華北地區10個種群的中華蒙潮蟲的主要演化是從北向南，但存在個別單倍型的混雜分布。本研究基于多基因的聯合分析與師二燕(2015)基于單基因的研究結果基本一致，這也與Avise(2009)對系統地理學歸納的第六種模式相吻合，即“即使一個單一的線粒體基因樹也能提供重要的生物學見解”。后期，將擴大采集范圍(由古北界向東洋界過渡)，完善該物種在中國的分布地區，結合形態標記與核基因進行中華蒙潮蟲演化路徑的研究。