徐 敏，鄧道貴，張海軍，王文平，晁秋杰

(淮北師范大學(xué)生命科學(xué)學(xué)院，資源植物生物學(xué)安徽省重點(diǎn)實(shí)驗(yàn)室，淮北 235000)

三種裸腹溞的16S rDNA序列分析與系統(tǒng)分類探討*

徐 敏，鄧道貴**，張海軍，王文平，晁秋杰

(淮北師范大學(xué)生命科學(xué)學(xué)院，資源植物生物學(xué)安徽省重點(diǎn)實(shí)驗(yàn)室，淮北 235000)

為了檢驗(yàn)已記錄的3種裸腹溞(發(fā)頭裸腹溞(Moinairrasa)、短型裸腹溞(Moinabrachiata)、微型裸腹溞(Moinamicrura))的系統(tǒng)分類，用試劑盒法分別提取3種裸腹溞的基因組DNA.利用特異性引物，通過(guò)PCR擴(kuò)增了3種裸腹溞的16S rDNA部分序列，并與來(lái)自GenBank中每個(gè)種類相似度較高的裸腹溞屬序列進(jìn)行分析.結(jié)果表明，3種裸腹溞的平均種間相似度為88.7%，堿基中A+T含量均明顯高于G+C含量.本研究的發(fā)頭裸腹溞的16S rDNA序列與GenBank所下載的多刺裸腹溞(Moinamacrocopa)的16S rDNA序列相似度為99%，遺傳距離(K2P雙參考模型)為0.5%，屬種內(nèi)范圍；兩個(gè)地區(qū)的短型裸腹溞測(cè)得的16S rDNA序列與GenBank下載的歐洲短型裸腹溞的16S rDNA序列序列相似度相對(duì)較低(88%～90%)，遺傳距離較大(13.2%～13.5%左右)，已達(dá)到屬內(nèi)種間分化水平.基于16S rDNA構(gòu)建的NJ樹和貝葉斯樹也支持以上結(jié)論.結(jié)果表明，本研究的發(fā)頭裸腹溞可能為多刺裸腹溞，本研究用的短型裸腹溞與GenBank下載的歐洲短型裸腹溞已經(jīng)達(dá)到種間分化的標(biāo)準(zhǔn).由于缺乏物種形態(tài)資料和其他分子標(biāo)記的對(duì)比，3種裸腹溞的分類地位還需進(jìn)行更深入的探討.

裸腹溞； 16S rDNA； 系統(tǒng)分類

裸腹溞屬(Moina)是枝角類的重要類群，通稱水蚤.本屬的描述種類有幾十種，我國(guó)記載的有十種，存在同物異名現(xiàn)象[1-2].早期的枝角類種類鑒定主要依據(jù)其形態(tài)特征，由于枝角類種類的分類特征有限以及形態(tài)變異十分豐富，加上近源個(gè)體之間可以雜交[3]，至今只有40%～50%的種類被準(zhǔn)確描述，部分種類分類比較含糊，其中一些種類可能以種復(fù)合體的形式存在[4].隨著分子生物學(xué)技術(shù)的進(jìn)步，利用基因標(biāo)記對(duì)枝角類進(jìn)行分子系統(tǒng)學(xué)研究已逐漸成為趨勢(shì)，Petrusek等以線粒體12S rRNA基因作為分子標(biāo)記驗(yàn)證了微型裸腹溞在兩個(gè)大陸間(歐洲中部和澳大利亞)的種類不統(tǒng)一[5]；Prosser等用不同引物在不同條件下進(jìn)行對(duì)照，擴(kuò)增了包括微型裸腹溞在內(nèi)的不同小型甲殼動(dòng)物的線粒體COI基因[6]，但沒有進(jìn)行系統(tǒng)分化方面的分析.國(guó)內(nèi)學(xué)者對(duì)薄皮溞屬(Leptodora)、西藏?cái)M溞(Daphniopsistibetana)、擬同形溞(Daphniasimiloides)的系統(tǒng)分類進(jìn)行了研究[7-9]，有關(guān)裸腹溞屬系統(tǒng)分化的研究尚未見報(bào)道.本研究通過(guò)對(duì)安徽巢湖、江蘇太湖、洪澤湖、淮北乾隆湖4個(gè)湖泊中已發(fā)現(xiàn)的短型裸腹溞、發(fā)頭裸腹溞和微型裸腹溞的16S rDNA序列進(jìn)行分析，通過(guò)與GenBank上下載的與本研究序列相似度最大的裸腹溞屬種類的16S rDNA進(jìn)行比對(duì)，探討他們之間的分子系統(tǒng)關(guān)系，提出了一些與傳統(tǒng)分類上不一致的問(wèn)題，以期為裸腹溞屬的遺傳多樣性和物種歸屬問(wèn)題提供基礎(chǔ)資料.

1 材料與方法

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

研究用巢湖的短型裸腹溞、發(fā)頭裸腹溞和微型裸腹溞來(lái)自2012年3月巢湖及其支流(南淝河)沉積物中休眠卵的孵化，并分離鑒定[2，10](表1)；洪澤湖發(fā)頭裸腹溞于2013年6月采集于洪澤湖；太湖微型裸腹溞于2013年5月采自太湖；乾隆湖短型裸腹溞于2014年5月采自淮北乾隆湖.幾種裸腹溞種類均在智能光照培養(yǎng)箱中用斜生柵藻(Scenedesmusobliquus)培養(yǎng)，培養(yǎng)溫度為25℃±1℃、光照時(shí)間L∶D=12h∶12h.實(shí)驗(yàn)時(shí)用吸管吸取單只個(gè)體饑餓一晝夜后用雙蒸水沖洗，離心管收集備用.

1.2 基因組DNA提取

采用微量樣品基因組DNA試劑盒(TIANamp Micro DNA Kit)提取裸腹溞屬DNA.由于裸腹溞屬種類屬于小型甲殼動(dòng)物，外被甲殼，在提取前用滅菌的牙簽搗碎溞體，有利于蛋白酶K對(duì)組織進(jìn)行消化.

表1 本研究中樣本核苷酸序列的來(lái)源及編號(hào)Tab.1 The origin and number of nucleotide sequence of samples in this study

1.3 PCR擴(kuò)增及序列測(cè)定

用于擴(kuò)增16S rDNA基因引物[11]：L2510 5′ to 3′(CGCCTGTTTAACAAAAACAT) H3059 5′ to 3′(CCGGTCTGAACTCAGATCATGT)

PCR反應(yīng)體系25μl，每個(gè)反應(yīng)體系內(nèi)含10×LA-Taq BufferⅡ 2.5μl、dNTP(2.5 mol)4.0μl、Mg離子(2.5 mol)0.5μl、模板DNA(100 ng/μl左右)1.0μl、Taq酶(5U/μl)0.25μl、上下游引物各(10μm)1.0μl，雙蒸水補(bǔ)足到25μl.

擴(kuò)增條件為95℃預(yù)變性3min，95℃變性45s，50℃退火45s，72℃延伸45s，共35個(gè)循環(huán)，72℃ 10min充分延伸，4℃結(jié)束.

擴(kuò)增產(chǎn)物經(jīng)電泳檢測(cè)后，選擇目的帶清晰且明亮的PCR反應(yīng)產(chǎn)物送交生工生物工程公司(上海)純化并測(cè)序.所有樣品均采用雙向測(cè)序，正反向?qū)Ρ群蟠_定分析用序列，用DNAStar中的SeqMan比對(duì)，并對(duì)比峰圖修除兩端不可靠的堿基，經(jīng)人工校正后獲得CDL有效堿基417個(gè)，CFL有效堿基411個(gè)，CWL有效堿基417個(gè)，HFL有效堿基381個(gè)，TWL有效堿基406個(gè)，QDL有效堿基491個(gè).其它從GenBank上下載的用于分析的序列見表1.

1.4 系統(tǒng)發(fā)生分析

用DNAStar對(duì)測(cè)序結(jié)果進(jìn)行人工校正，并去除兩端冗余序列.使用DNAStar 軟件比較序列差異百分比；用Clustalx 1.81軟件多重序列比對(duì)；用DNAspV5分析變異位點(diǎn)；用Mega 4.1計(jì)算不同序列間轉(zhuǎn)換/顛換數(shù)，生成序列差異位點(diǎn)表，并計(jì)算遺傳距離和構(gòu)建NJ進(jìn)化樹.采用K2P雙參模型計(jì)算兩兩遺傳距離，系統(tǒng)樹設(shè)置采用P-distance模型，以大型溞(Daphniamagna)為外群，Bootstrap Replications1000次；用Mrbayes 3.2構(gòu)建貝葉斯樹，替換模型設(shè)置參數(shù)nst=6(GTR模型)，位點(diǎn)速率變異設(shè)置invgamma，其余設(shè)置為默認(rèn)值；采用mcmc法運(yùn)算500000代，每100代取樣一次，在舍棄25%老化樣本后，用剩余樣本構(gòu)建一致樹.

2 結(jié)果

CDL、QDL、CFL、CWL、HFL、TWL在DNAStar中經(jīng)過(guò)兩兩序列比對(duì)，TWL與CWL相似度為100%，HFL與CFL相似度為99%，有兩個(gè)位置發(fā)生轉(zhuǎn)換；CDL與QDL相似度為99%，有一個(gè)位置發(fā)生轉(zhuǎn)換，另一個(gè)位置發(fā)生顛換.堿基中A+T含量(67.6%)明顯高于G+C含量(32.4%)；3個(gè)物種間平均種間相似度為88.7%，用DNAsp v5進(jìn)行序列變異位點(diǎn)分析，上述6個(gè)裸腹溞屬個(gè)體可識(shí)別位點(diǎn)(包括缺失與丟失的位點(diǎn))381個(gè)，其中保守位點(diǎn)313個(gè)，變異位點(diǎn)68個(gè)，單一位點(diǎn)3個(gè)，簡(jiǎn)約信息點(diǎn)65個(gè)；用MEGA軟件分析平均轉(zhuǎn)換/顛換(si/sv)數(shù)為0.78，顛換數(shù)明顯高于轉(zhuǎn)換數(shù).與來(lái)自GenBank上Moina的16S rDNA進(jìn)行序列比對(duì)，計(jì)算遺傳距離見表3；就構(gòu)建的NJ樹(圖1a)和貝葉斯樹(圖1b)來(lái)看，兩者主要分枝基本一致.

圖1 基于9個(gè)裸腹溞屬個(gè)體16S rDNA序列構(gòu)建的系統(tǒng)發(fā)生樹(NJ樹)(a)和貝葉斯樹(b)Fig.1 NJ tree(a) and Bayes tree(b) based on the sequence of 16S rDNA in nine Moina individuals

3 討論

3.1 3種裸腹溞的16S rDNA序列分析

線粒體rRNA基因是非編碼蛋白質(zhì)的基因，進(jìn)化過(guò)程中未受到密碼子編碼的選擇壓力影響，進(jìn)化速率相對(duì)較低且比較保守[12].本研究所測(cè)3個(gè)裸腹溞屬種類的個(gè)體中，種間的16S rDNA序列相似度百分比在87%～90%之間(表2)，這與臂尾輪蟲屬的種間序列相似度較一致[13].從已測(cè)出的3種裸腹溞的16S rDNA部分堿基序列比對(duì)表可以看出，3個(gè)裸腹溞屬種類的變異區(qū)相對(duì)比較集中，以序列最短的洪澤湖發(fā)頭裸腹溞序列作為基準(zhǔn)定位，在第1位到第163位堿基之間僅有12處發(fā)生堿基的改變，而從164到282之間卻有51處發(fā)生變異；堿基中A+T(67.6%)與Petrusek等研究的多刺裸腹溞和微型裸腹溞12S rRNA基因的A+T含量(69.4%～72.5%)接近[5]，較高的A+T含量是無(wú)脊椎動(dòng)物線粒體DNA中的普遍現(xiàn)象[14-16].

3.2 3個(gè)裸腹溞屬種類的系統(tǒng)分類探討

經(jīng)過(guò)DNAStar序列比對(duì)，CFL與HFL的16S rDNA序列相似度為99%，遺傳距離(K2P雙參模型，下同)為0.5%，這可能是地理隔離造成的種內(nèi)遺傳差異.而CFL和HFL的16S rDNA序列與從GenBank下載的NCL的遺傳距離也較小(0.5%～1.1%)，遠(yuǎn)低于枝角類線粒體序列種內(nèi)差異上限5%的分歧度[17].NJ樹(P-distance模型)也支持NCL與CFL和HFL聚為一類.因此，就16S rDNA分析看，我國(guó)的CFL和HFL應(yīng)與來(lái)自歐洲的NCL為同一個(gè)物種.通過(guò)查閱文獻(xiàn)發(fā)現(xiàn)[18-21]，發(fā)頭裸腹溞這一名稱僅限國(guó)內(nèi)學(xué)者使用，國(guó)外學(xué)者較少使用，而Goulden把發(fā)頭裸腹溞歸到M.affinis一類[22].本研究的發(fā)頭裸腹溞雌性體長(zhǎng)一般大于1.00mm，孵育囊飽滿時(shí)，殼高常大于殼長(zhǎng).頭部膨大，背側(cè)叢生剛毛，無(wú)眼上凹，卵鞍長(zhǎng)卵形，內(nèi)儲(chǔ)2個(gè)冬卵.后腹部寬大，各側(cè)有10個(gè)羽狀肛刺.這與《中國(guó)動(dòng)物志·淡水枝角類》[2]中對(duì)發(fā)頭裸腹溞的描述一致.通過(guò)形態(tài)觀察發(fā)現(xiàn)，本研究的發(fā)頭裸腹溞在肛刺數(shù)量、卵鞍形狀及卵鞍內(nèi)冬卵的數(shù)量等方面與M.affinis存在較大的區(qū)別，與多刺裸腹溞的卵鞍形狀、卵鞍內(nèi)冬卵的數(shù)量也不同[2，22].Elmoor-Loureiro等報(bào)道了一種新發(fā)現(xiàn)的多刺裸腹溞，其外形、卵鞍特征和冬卵數(shù)量與本研究的發(fā)頭裸腹溞相似[23].結(jié)合16S rDNA序列的分析結(jié)果，推測(cè)本研究的發(fā)頭裸腹溞應(yīng)為多刺裸腹溞.由于缺乏更多的形態(tài)和分子對(duì)比資料，有關(guān)國(guó)內(nèi)發(fā)頭裸腹溞的歸屬問(wèn)題仍需進(jìn)一步的探討.

表2 9個(gè)裸腹溞屬個(gè)體16S rDNA的遺傳距離(對(duì)角線下)和序列相似度(對(duì)角線上)Tab.2 The genetic distances(below diagonal) and sequence similarity(above diagonal) of 16S rDNA in nine Moina individuals

本研究的微型裸腹溞雌性體長(zhǎng)在0.65～0.83mm之間，復(fù)眼很大，位于頭頂；頭部很大，背側(cè)有較深的眼上凹；后腹部短，側(cè)面有羽狀肛刺6個(gè)，叉狀肛刺1個(gè)；卵鞍近乎卵圓型，內(nèi)儲(chǔ)冬卵1個(gè).這些形態(tài)特征與《中國(guó)動(dòng)物志.淡水枝角類》[2]中關(guān)于微型裸腹溞的描述較一致.來(lái)自巢湖與太湖的微型裸溞相似度為100%，但在NCBI的BLAST中未找到與之相似度較高的物種，只有M.brachiata(JN651425)與之有92%的相似度；在GenBank中也沒有找到微型裸腹溞16S rDNA序列的數(shù)據(jù)與本研究進(jìn)行比對(duì).因此，關(guān)于微型裸腹溞的分子系統(tǒng)學(xué)仍需做更多的研究.

本研究的短型裸腹溞雌性體長(zhǎng)1.00～1.50mm，淺青灰色，殼瓣背緣均勻弓起，頭部無(wú)毛發(fā)，長(zhǎng)度超過(guò)殼長(zhǎng)的一半，向下傾斜，背側(cè)有大而深的眼上凹，后腹部較長(zhǎng)，各側(cè)有羽狀肛刺8個(gè)，叉狀肛刺一個(gè)，這與《中國(guó)動(dòng)物志·淡水枝角類》[2]中的短型裸腹溞描述較一致.卵鞍內(nèi)儲(chǔ)冬卵一個(gè)，這與文獻(xiàn)[22]中對(duì)短型裸腹溞的描述一致.本研究中，巢湖(CDL)和淮北乾隆湖(QDL)的短型裸腹溞相似度為99%，遺傳距離為0.5%.他們與GenBank下載的歐洲短型裸腹溞(NDL1、NDL2)的相似度相對(duì)較小(88%～90%)，遺傳距離(13.2%～13.5%)較大.其遺傳距離遠(yuǎn)遠(yuǎn)大于枝角類小于5%的種內(nèi)范圍[17]，達(dá)到端足類鉤蝦屬的種間分化水平[24].就NJ樹和貝葉斯樹分析，CDL、QDL與下載的歐洲短型裸腹溞(NDL1和NDL2)的分枝較遠(yuǎn)，也已經(jīng)不在一個(gè)種的界限內(nèi).造成這種差距的原因可能有兩點(diǎn)：一是長(zhǎng)期遠(yuǎn)距離的地理隔離，不同環(huán)境發(fā)生的自然選擇不同導(dǎo)致遺傳上的較大差異；二是存在形態(tài)相似而基因有差別的物種或稱隱種[25]，這種現(xiàn)象已經(jīng)發(fā)生在微型裸腹溞的跨洲(歐洲中部和澳大利亞)分化上[5].圓形盤腸溞曾被看作是枝角類中一個(gè)分布較廣、生活習(xí)性多樣化的種，但后來(lái)發(fā)現(xiàn)它實(shí)際上是以一種復(fù)合體的形式存在[26-27].其他廣泛分布的短鈍溞(D.obtusa)、長(zhǎng)額象鼻溞(Bosminalongirostris)也依然存在分類上的誤區(qū)[4].關(guān)于枝角類物種的統(tǒng)一性問(wèn)題一直受到質(zhì)疑[28-29].近年來(lái)，隨著分子工具的使用，大量的形態(tài)結(jié)構(gòu)和分子研究都證明了非統(tǒng)一性物種的存在，這將會(huì)改變我們對(duì)枝角類的發(fā)展史、生物地理學(xué)和物種多樣性的認(rèn)識(shí)[30-33].

[1] 趙 文.水生生物學(xué).北京：中國(guó)農(nóng)業(yè)大學(xué)出版社，2005：181-185.

[2] 蔣燮治，堵南山.中國(guó)動(dòng)物志·淡水枝角類.北京：科學(xué)出版社，1979：1-160.

[3] Billiones R， Brehm M， Klee Jetal. Genetic identification ofHyalodaphniaspecies and interspecific hybrids.Hydrobiologia， 2004， 526：43-53.

[4] Korovchinsky NM. How many species of Cladocera are there？Hydrobiologia， 1996， 321：191-204.

[5] Petrusek A， Cerny M， Audenaert E. Large intercontinental differentiation ofMoinamicrura(Crustacea： Anomopoda)： one less cosmopolitan cladoceran？Hydrobiologia， 2004， 526：73-81.

[6] Prosser S， Martinez-Arce A， Elias-Gutierrez M. A new set of primers for COI amplification from freshwater microcrustaceans.MolecularEcologyResources， 2013， 13(6)：1151-1155.

[7] Xu L， Han BP， Damme KVetal. Biogeography and evolution of the Holarctic zooplankton genusLeptodora(Crustacea： Branchiopoda： Haplopoda).JournalofBiogeography， 2011， 38：359-370.

[8] 趙 文，李 睿.4個(gè)品系西藏?cái)M溞12S rRNA基因序列及其分子進(jìn)化研究.大連海洋大學(xué)學(xué)報(bào)，2012，27(4)：300-305.

[9] 顧楊亮，徐 磊，林秋奇等.新亞種——中華擬同形溞(Daphniasimiloidessinensis).生態(tài)科學(xué)，2013，32(3)：308-312.

[10] Deng DG， Xie P， Zhou Qetal. Field and experimental studies on the combined impacts of cyanobacterial blooms and small algae on crustacean zooplankton in a large， eutrophic subtropical Chinese lake.Limnology， 2008， 9：1-11.

[11] Bouchon D， Souty-Grosset C， Raimond R. Mitochondrial DNA variation and markers of species identity in two penaeid shrimp species：PenaeusmonodonFabricius andP.japonicusBate.Aquaculture， 1994， 127(2)：131-144.

[12] 竇紅霞，張彤晴，許志強(qiáng)等.長(zhǎng)臂蝦亞科9個(gè)種系統(tǒng)發(fā)育關(guān)系的16S rDNA序列分析.淡水漁業(yè)，2008，38(1)：40-43.

[13] 程雙懷，席貽龍.基于16S rDNA基因序列探討十種臂尾輪蟲的系統(tǒng)關(guān)系和分類地位.動(dòng)物分類學(xué)報(bào)，2009，34(4)：934-940.

[14] Anderson FE. Phylogeny and historical biogeography of theloliginidsquids(Mollusca： Cephalopoda) based on mitochondrial DNA sequence data.MolecularPhylogeneticsandEvolution， 2000， 15(2)：191-214.

[15] 高天翔，張秀梅，吉崎悟朗等.日本絨螯蟹線粒體DNA序列研究I.12S rRNA.青島海洋大學(xué)學(xué)報(bào)，2000，30(1)：43-47.

[16] 劉 君，李 琪，孔令鋒等.基于線粒體COI的DNA條形碼技術(shù)在貽貝科種類鑒定中的應(yīng)用.水生生物學(xué)報(bào)，2011，35(5)：874-881.

[17] Adamowicz SJ， Petrusek A， Colbourne JKetal. The scale of divergence： a phylogenetic appraisal of intercontinental allopatric speciation in a passively dispersed freshwater zooplankton genus.MolecularPhylogeneticsandEvolution， 2009， 50(3)：423-436.

[18] 鄭又雄.發(fā)頭裸腹溞卵和胚胎發(fā)育形態(tài)研究.動(dòng)物學(xué)雜志，1993，28(6)：5-8.

[19] 高露嬌，賴 偉，陳立僑.發(fā)頭裸腹溞(MoinairrasaBrehm)種群增長(zhǎng)與食物濃度及光照的關(guān)系.水產(chǎn)科技情報(bào)，1997，24(2)：61-64.

[20] 饒小珍，許友勤，陳寅山等.發(fā)頭裸腹溞——一種新的溞類毒性試驗(yàn)材料.生物技術(shù)，2000，10(1)：17-19.

[21] 陳麗娜，李玉穎，鄧道貴等.種群密度和培養(yǎng)體積對(duì)發(fā)頭裸腹溞生長(zhǎng)和繁殖的影響.應(yīng)用生態(tài)學(xué)報(bào)，2012，23(7)：1959-1964.

[22] Goulden CE. The systematics and evolution of the Moinidae.TransactionofAmericanPhilosopnicalSociety， 1968， 58(6)：15-41.

[23] Elmoor-Loureiro LMA， Santangelo JM， Lopes PMetal. A new report ofMoinamacrocopa(Straus，1820)(Cladocera， Anomopoda) in South America.BrazilianJournalofMedicalandBiologicalResearch， 2010， 70(1)：225-226.

[24] Costa FO， de Waard J， Boutillier Jetal. Biological identifications through DNA barcodes： the case of the Crustacea.CanadianJournalofFisheriesandAquaticSciences， 2007， 64：272-295.

[25] 項(xiàng)賢領(lǐng)，席貽龍，溫新利等.萼花臂尾輪蟲種復(fù)合體遺傳分化的空間格局.動(dòng)物學(xué)研究，2010，3(31)：205-220.

[26] Frey DG. Questions concerning cosmopolitanism in Cladocera.ArchivfürHydrobiologia， 1982， 145：5-17.

[27] Frey DG. The taxonomy and biogeography of the Cladocera.Hydrobiologia， 1987， 145：19-28.

[28] Dumont HJ， Negrea SV. Introduction to the class Branchiopoda. In： Dumont HJ ed. Guides to the identification of the microinvertebrates of the continental waters of the world 19. Leiden： Backhuys Publishers， 2002.

[29] Frey DG. The taxonomic and phylogenetic significance of the head pores of the Chydoridae(Cladocera).InternationaleRevuedergesamtenHydrobiologie， 1959， 44：27-50.

[30] Adamowicz SJ， Hebert PDN， Marinone MC. Species diversity and endemism in theDaphniaof Argentina： A genetic investigation.ZoologicalJournaloftheLinneanSociety， 2004， 140：171-205.

[31] Cox AJ， Hebert PDN. Colonization， extinction， and phylogeographic patterning in a freshwater crustacean.MolecularEcology， 2001， 10：371-386.

[32] Schwenk KN， Posada D， Hebert PDN. Molecular systematics of European Hyalodaphnia： the role of contemporary hybridization in ancient species.ProceedingsoftheRoyalSocietyofLondon， 2000， B267：1833-1842.

[33] Taylor DJ， Ishikane CR， Haney RA. The systematics of Holarctic bosminids and a revision that reconciles molecular and morphological evolution.LimnologyandOceanography， 2002， 47：1486-1495.

Molecular phylogeny of threeMoinaspecies based on 16S rDNA gene sequences

XU Min， DENG Daogui， ZHANG Haijun， WANG Wenping & CHAO Qiujie

(AnhuiKeyLaboratoryofResourceandPlantBiology，SchoolofLifeScience，HuaibeiNormalUniversity，Huaibei235000，P.R.China)

In order to examine the systematic classification of threeMoinaspecies，M.irrasa，M.brachiataandM.micrurarecorded previously， genomic DNA of threeMoinaspecies were extracted through the kit method. PCR was used to amplify partial sequences of 16S rDNA of threeMoinaspecies by employing specific primers， and same sequence with higher matching rate of eachMoinaspecies from the GenBank were analysed. Based on 16S rDNA， the average interspecific similarity rate among threeMoinaspecies was 88.7%， and the content of A+T was obviously higher than that of G+C.The similarity ofM.irrasain this study andM.macrocopafrom the GenBank was 99%， and the divergence rate was 0.5%， within a range of same species. The similarity ofM.brachiatain this study and the GenBank was relatively low(88%-90%)， and the divergence rate was 13.2%-13.5%， without a range of same species. The NJ tree and Bayes tree based on 16S rDNA did also support the above results.The results suggested thatM.irrasain this study andM.macrocopafrom the GenBank might be the same species， andM.brachiatafrom this study and the GenBank had reached the standard of species divergence. Due to the lack of morphology and other genes analysis， the classification status of threeMoinaspecies in this study mentioned need to be further investigated.

Moina； 16S rDNA； systematic classification

*國(guó)家自然科學(xué)基金項(xiàng)目(31370470)資助.2014-03-18收稿；2014-07-02收修改稿.徐敏(1983～)，女，碩士研究生；E-mail：minxuzheng@126.com.

*通信作者；E-mail：dengdg@263.net.