• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于線粒體Cyt b基因的中國(guó)南方唇魚骨(Hermibarbuslabeo)和間魚骨(Hermibarbusmedius)種群分化及親緣地理研究

    2016-11-28 09:17:33藍(lán)昭軍范明君黃小林
    生態(tài)學(xué)報(bào) 2016年19期
    關(guān)鍵詞:魚骨親緣海南島

    藍(lán)昭軍,范明君,黃小林,趙 俊

    1 廣東省水產(chǎn)健康安全養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,廣東省高等學(xué)校生態(tài)與環(huán)境科學(xué)重點(diǎn)實(shí)驗(yàn)室,廣州市亞熱帶生物多樣性與環(huán)境生物監(jiān)測(cè)重點(diǎn)實(shí)驗(yàn)室,華南師范大學(xué)生命科學(xué)學(xué)院,廣州 510631 2 韶關(guān)市水產(chǎn)研究所, 韶關(guān) 512006 3 深圳華大水產(chǎn)科技有限公司, 深圳 518083 4 中國(guó)水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所, 廣州 510300

    ?

    基于線粒體Cytb基因的中國(guó)南方唇魚骨(Hermibarbuslabeo)和間魚骨(Hermibarbusmedius)種群分化及親緣地理研究

    藍(lán)昭軍1,2,范明君1,3,黃小林1,4,趙 俊1,*

    1 廣東省水產(chǎn)健康安全養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室,廣東省高等學(xué)校生態(tài)與環(huán)境科學(xué)重點(diǎn)實(shí)驗(yàn)室,廣州市亞熱帶生物多樣性與環(huán)境生物監(jiān)測(cè)重點(diǎn)實(shí)驗(yàn)室,華南師范大學(xué)生命科學(xué)學(xué)院,廣州 510631 2 韶關(guān)市水產(chǎn)研究所, 韶關(guān) 512006 3 深圳華大水產(chǎn)科技有限公司, 深圳 518083 4 中國(guó)水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所, 廣州 510300

    為了解中國(guó)南方唇魚骨(Hermibarbuslabeo)和間魚骨(Hermibarbusmedius)的種群分化、親緣地理格局及物種有效性,作者對(duì)唇魚骨 8個(gè)水系及間魚骨 9個(gè)水系共148尾樣本的Cytb基因全序列進(jìn)行了測(cè)定。在所有序列中,共有128個(gè)變異位點(diǎn),共檢測(cè)出了41個(gè)單倍型,單倍型多樣性為0.954,核苷酸多樣性為0.02153?;贑ytb基因全序列構(gòu)建的NJ樹(shù)顯示,中國(guó)南方的唇魚骨與間魚骨合聚為兩大支,其中韓江和九龍江的全部唇魚骨樣本組成了Ⅰ支,而其他水系的唇魚骨和間魚骨樣本則組成了Ⅱ支。兩支系間的遺傳距離為5.1%,而唇魚骨和間魚骨之間的遺傳距離為3.2%?,F(xiàn)有證據(jù)不支持間魚骨與唇魚骨達(dá)到種一級(jí)的分化。單倍型網(wǎng)絡(luò)圖顯示,韓江、九龍江種群和其他水系的種群分化較大;海南島三大水系種群和漠陽(yáng)江種群的單倍型分支與代表珠江水系單倍型的分支之間的親緣關(guān)系較近,與其他地理區(qū)間種群的親緣關(guān)系則相對(duì)較遠(yuǎn);西江可能為珠江水系、漠陽(yáng)江水系和海南島三大水系間魚骨種群的擴(kuò)散中心,其中一支往東向北江和東江擴(kuò)散;另一支往南向海南島擴(kuò)散,海南島三大水系種群爾后在冰期海退時(shí)向北擴(kuò)散至漠陽(yáng)江水系及珠江水系。AMOVA分析表明,唇魚骨和間魚骨地理區(qū)之間變異約占54.50%,地理區(qū)內(nèi)種群間變異約占18.64%,種群內(nèi)的變異占26.86%,這說(shuō)明,唇魚骨和間魚骨種群Cytb的遺傳分化主要是來(lái)自地理區(qū)之間。錯(cuò)配分析及中性檢驗(yàn)顯示,全部種群、唇魚骨種群、間魚骨種群在歷史上均沒(méi)有發(fā)生過(guò)明顯的擴(kuò)張,兩個(gè)mtDNA支系亦未發(fā)生過(guò)擴(kuò)張,而海南島的昌化江種群曾發(fā)生過(guò)種群擴(kuò)張。

    唇魚骨;間魚骨;種群分化;親緣地理;細(xì)胞色素b基因;中國(guó)南方

    親緣地理學(xué),又稱譜系生物地理學(xué),是研究有密切親緣關(guān)系的種間和種內(nèi)支系現(xiàn)有的地理分布格局的形成過(guò)程和形成機(jī)制的學(xué)科。生物在漫長(zhǎng)的演化過(guò)程中,受各種因素(地質(zhì)因素、氣候因素等)的影響,都會(huì)在其演化的過(guò)程中留下印跡,如影響生物的地理隔離、遷移、區(qū)域擴(kuò)張和滅絕等,而生物的這些活動(dòng)都會(huì)具有影響一定區(qū)域內(nèi)物種基因流的作用,改變遺傳變異的分布,從而產(chǎn)生新的物種進(jìn)化模式[1]。尤其是,地質(zhì)地貌的格局對(duì)于種群基因流的改變和分化往往起著重要的作用[2- 6]。而分子親緣地理學(xué)是采用分子生物技術(shù),在分子水平上探討物種的系統(tǒng)地理格局的形成機(jī)制,進(jìn)一步闡述其起源進(jìn)化歷史,分析區(qū)域類群在時(shí)間上和空間上的發(fā)展變化,從而重建生物區(qū)系的進(jìn)化歷史過(guò)程[7]。

    唇魚骨(Hemibarbuslabeo)和間魚骨(Hemibarbusmedius)均隸屬于鯉形目(Cypriniformes),鯉科(Cyprinidae),鮈亞科(Gobioninae),魚骨屬(Hemibarbus)[8- 9],兩者為近緣種。間魚骨是在原唇魚骨種群中劃分出的一個(gè)種[8],在此之前,一直被歸為唇魚骨[10- 14]。1995年,樂(lè)佩琦對(duì)我國(guó)魚骨屬魚類進(jìn)行分類整理后認(rèn)為,唇魚骨華南地區(qū)及西南部分地區(qū)的種群與唇魚骨其他地方種群的差別已經(jīng)達(dá)到了物種級(jí)的水平,并以海南島的12尾標(biāo)本作為模式標(biāo)本,將華南及西南部分地區(qū)的唇魚骨種群描述為一新種——間魚骨(H.mediusYue),并指出其與唇魚骨的鑒別特征主要為吻長(zhǎng)較短,唇稍薄,不甚發(fā)達(dá),下唇側(cè)葉稍狹窄,無(wú)皺褶,背鰭刺較細(xì)弱及鰓耙數(shù)較少等[8]。1998年,樂(lè)佩琦等較為詳細(xì)地記載唇魚骨和間魚骨的地理分布,其中唇魚骨分布于我國(guó)臺(tái)灣各水系和閩江、錢塘江、長(zhǎng)江、黃河等直至黑龍江水系;間魚骨分布于華南、西南部分地區(qū)各水系[9]。

    然而,目前不同學(xué)者對(duì)于兩者的分布區(qū)域與樂(lè)佩琦等的觀點(diǎn)不盡相同[15- 18]。唇魚骨和間魚骨之間的鑒別特征較為模糊,且某些性狀可能存在連續(xù)變化,不易區(qū)分。因此,造成對(duì)兩者分布區(qū)描述差異的原因可能為物種鑒定存在偏差。此外,兩者的分布均較為廣泛,是研究親緣地理的理想素材。有關(guān)唇魚骨和間魚骨的親緣地理的研究,目前已有林宗鍵[17]及Lin等[18]基于線粒體控制區(qū)序列對(duì)唇魚骨的報(bào)道,但該研究并未涉及珠江等水系的種群。因此,作者在中國(guó)南方唇魚骨和間魚骨的分布區(qū),采集不同地理種群的樣本,對(duì)其線粒體細(xì)胞色素b(Cytb)基因序列進(jìn)行了測(cè)定與分析,旨在澄清二者物種有效性的前提下,探討其種群分化和分布特征,推測(cè)其種群演化歷史及親緣地理模式。

    1 材料與方法

    1.1 樣品采集

    本研究所用樣品均為2009年9月—2010年8月采自貴州省、廣西壯族自治區(qū)、廣東省、海南省、福建省、江西省和浙江省等地,所有樣本根據(jù)樂(lè)佩琦等提出的分類鑒定特征[8- 9]進(jìn)行鑒定,其中間魚骨樣本共92尾;唇魚骨樣本56尾。將所有個(gè)體按所采樣的水系進(jìn)行種群劃分,即用水系名稱對(duì)種群進(jìn)行命名。具體的樣品信息及采樣點(diǎn)見(jiàn)表1。

    1.2 DNA的提取及測(cè)序

    基因組 DNA的提取參考慶寧等[19]的方法,即采用從上海生工生物工程技術(shù)服務(wù)有限公司購(gòu)買的“基因組 DNA抽提試劑盒”所示方法,略有改動(dòng)。PCR反應(yīng)體系:反應(yīng)體系約25μL,按順序分別將PCR Master Mix 12.5μL、Primer1 3μL、Primer2 3μL、模板DNA 4μL及滅菌雙蒸水2.5μL混合,每一次反應(yīng)都用ddH2O代替模板DNA做陰性對(duì)照。反應(yīng)在PTC100或PTC200型PCR儀上進(jìn)行,具體擴(kuò)增程序如下:A: 95℃預(yù)變性4min,B: 94℃變性40sec,C: 58℃退火45sec,D: 72℃延伸1min,E: 重復(fù)步驟B—D45次;F:最后72℃延伸10min。PCR及測(cè)序所用引物均為通用引物L(fēng)14724和H15915,由上海英駿生物技術(shù)有限公司合成,其序列分別為:L14724:5′-GACTTG AAA AACCAC CGTTG- 3′;H15915:5′-CTCCGATCTCCGGATTACAAGAC- 3′。PCR產(chǎn)物用1%瓊脂糖電泳檢驗(yàn)?zāi)康钠魏?將PCR產(chǎn)物送上海生工生物工程技術(shù)服務(wù)有限公司純化并測(cè)序。

    1.3 數(shù)據(jù)分析

    對(duì)測(cè)得的序列使用Clustal X 1.81 軟件進(jìn)行排序比對(duì)并輔以人工校對(duì)。用軟件MEGA 5.0[20]計(jì)算所分析序列的堿基組成。采用Kimura′s two-parameter模型計(jì)算遺傳距離,并采用該軟件中的鄰接法(Neighbor-joining method,NJ),以似鮈(Pseudogobiovaillanti)作為外類群,構(gòu)建分支系統(tǒng)樹(shù)。用軟件DNAsp 5.0[21]對(duì)各種群序列的多樣性進(jìn)行分析,統(tǒng)計(jì)核苷酸多樣性(Nucleotide diversity,π)與單倍型多樣性(Haplotype diversity,h),并以此估計(jì)不同水系種群的遺傳多樣性;并利用該軟件進(jìn)行錯(cuò)配分布(Mismatch-distribution)分析及Tajima′s D和Fu and Li′s檢驗(yàn),以檢驗(yàn)種群是否發(fā)生過(guò)擴(kuò)張。使用Network 4.6 軟件,以Median-joining法構(gòu)建各單倍型之間網(wǎng)絡(luò)關(guān)系圖,并對(duì)圖中各個(gè)單倍型進(jìn)行群體對(duì)應(yīng)關(guān)系分析。利用Arlequin Ver 2000[22]軟件進(jìn)行分子變異分析(AMOVA),測(cè)試種群間與地理區(qū)間遺傳變異程度。

    2 結(jié)果與分析

    2.1 Cyt b 基因的序列變異

    本研究共得到148尾唇魚骨與間魚骨的Cytb基因全序列,序列長(zhǎng)1140bp,無(wú)堿基的插入或缺失,Cytb基因所有序列中,共128個(gè)變異位點(diǎn),占全序列的11.23%,其中14個(gè)為單個(gè)多態(tài)位點(diǎn),114個(gè)為簡(jiǎn)約信息位點(diǎn)。大多數(shù)的變異發(fā)生在密碼子第3位(82.81%),第1位(14.02%)和第2位(3.12%)的變異較低。唇魚骨與間魚骨在Cytb序列中A、C、T、G平均含量分別為27.5%、28.4%、28.4%和15.6%,G的含量明顯低于其他3種堿基含量,表現(xiàn)出很強(qiáng)的堿基組成偏向性。其中A+T的含量55.9%,G+C的含量為44.0%,A+T的含量明顯高于G+C的含量。

    2.2 Cyt b 基因的單倍型多樣性和遺傳多樣性

    在唇魚骨56尾樣本中,檢測(cè)出19個(gè)單倍型,其中17個(gè)為特有單倍型,2個(gè)為與間魚骨共享的單倍型;間魚骨92尾樣本,檢測(cè)出24個(gè)單倍型,其中22個(gè)為獨(dú)有單倍型,2個(gè)為與唇魚骨共享的單倍型(表1)。在全部148尾樣本中,共檢測(cè)出了41個(gè)單倍型,單倍型多樣性為0.954。

    在一些不同的種群中,存在共享單倍型,表明這些種群的親緣關(guān)系十分密切。其中海南島三大水系之一南渡江(NAN)、廣東西部獨(dú)立入海水系漠陽(yáng)江(MOY)及珠江水系第二大支流北江(BEI)的種群共享單倍型NAN1-MOY1-BEI5;西江支流桂江(GUI)和柳江(LIU)的間魚骨種群與長(zhǎng)江水系支流湘江(XIA)的唇魚骨種群共享單倍型GUI1-LIU1-XIA1。

    單倍型多樣性(h)和核苷酸多樣性(π)分析表明(表1),148尾唇魚骨和間魚骨樣本的單倍型多樣性為0.954,核苷酸多樣性為0.02153,遺傳多樣性較高。核苷酸多樣性(π)分析表明(表1),湘江、桂江、北江、東江、九龍江、閩江種群的核苷酸遺傳多樣性較高,均在0.002以上。

    2.3 種群間遺傳距離及遺傳分化

    統(tǒng)計(jì)表明,唇魚骨和間魚骨之間的遺傳距離為3.2%,兩者遺傳差異不大。表2分別顯示了不同水系的唇魚骨與間魚骨基于Cytb基因的遺傳距離和及遺傳分化指數(shù)。由表2可知,在各地理區(qū)(如海南島地區(qū)、東南地區(qū)等)和水系內(nèi)部(如長(zhǎng)江水系和珠江水系)的遺傳距離、遺傳分化指數(shù)較低;珠江水系和長(zhǎng)江水系的遺傳距離、遺傳分化指數(shù)亦較低。在唇魚骨和間魚骨各個(gè)種群之間,韓江、九龍江和閩江種群與其他各種群的遺傳距離較大,遺傳分化指數(shù)較高,而余下各種群之間的遺傳距離和遺傳分化指數(shù)則較低。

    2.4 系統(tǒng)發(fā)育分析

    為更好地了解中國(guó)南方地區(qū)唇魚骨與間魚骨各種群之間的親緣關(guān)系,以似鮈為外類群,對(duì)148尾樣品的Cytb基因中的41個(gè)單倍型進(jìn)行系統(tǒng)發(fā)育分析,利用最大鄰接法(NJ)構(gòu)建分子系統(tǒng)發(fā)育樹(shù)(圖1)。

    在我國(guó)南方的唇魚骨與間魚骨種群間,聚成Ⅰ和Ⅱ兩支。韓江和九龍江的全部樣本組成了Ⅰ支,而其他水系的樣本則組成了Ⅱ支,支系Ⅰ和Ⅱ之間的遺傳距離為5.1%,這說(shuō)明韓江和九龍江種群已和其他種群發(fā)生了一定程度的遺傳分化。在Ⅱ支中,又可分為A和B個(gè)兩姐妹群。A群包含了海南島三大水系、漠陽(yáng)江、柳江、賀江、北江、錢塘江和甌江的全部樣本及桂江、東江、湘江和閩江的部分樣本;B群則包含了沅江的全部樣本及桂江、湘江和閩江的部分樣本。在支系Ⅱ的A群中,又可分為A1和A2兩個(gè)姐妹群。A1有兩個(gè)分支,分別為海南島三大水系、漠陽(yáng)江的全部樣本和北江水系一尾樣本組成一個(gè)分支,說(shuō)明海南島三大水系種群與漠陽(yáng)江水系種群的親緣關(guān)系較近;而賀江、柳江的全部標(biāo)本及湘江、北江和東江的部分樣本組成另一支,這則說(shuō)明珠江水系三大支流之間的遺傳差異不大,而珠江水系與湘江水系之間的遺傳差異亦較低。A2則包含了錢塘江、甌江的全部樣本及閩江部分樣本和東江的一尾樣本。

    就唇魚骨與間魚骨這兩個(gè)物種而言,除了支系Ⅰ全部為唇魚骨、A1其中的一個(gè)分支全部為間魚骨之外,其余均含有唇魚骨與間魚骨兩個(gè)物種,這也說(shuō)明兩者間遺傳差異不大。

    2.5 單倍型網(wǎng)絡(luò)分析

    圖2顯示了中國(guó)南方唇魚骨和間魚骨基于Cytb基因的單倍型網(wǎng)絡(luò)親緣關(guān)系,大多數(shù)單倍型之間的變異為1—7步,遺傳差異不大。但是桂江和湘江種群的共享單倍型GUI2-XIA2與韓江和九龍江種群之間的變異近40步;桂江與錢塘江、甌江與閩江種群的變異達(dá)約20步,表明這些種群的單倍型間存在一定的遺傳分化。

    物種Species水系Riversystems采集地Samplelocations樣品數(shù)Samplesize單倍型(個(gè)體數(shù))Haplotypes(numberofindividuals)變異位點(diǎn)數(shù)Numberofpolymorphicsites(S)單倍型數(shù)Numebrofhaplotypes單倍型多樣性Haplotypediversity(h)核苷酸多樣性Nucleotidediversity(π)間魚骨H.medius昌化江(CHJ)五指山10CHJ1(3)、CHJ2(2)、CHJ3(1)、CHJ4(1)、CHJ5(2)450.8670.00121南渡江(NAN)屯昌9NAN1-MOY1-BEI5*(9)010.0000.00000萬(wàn)泉河(WAN)瓊海9WAN1(7)、WAN2(2)120.3890.00034漠陽(yáng)江(MOY)陽(yáng)春11NAN1-MOY1-BEI5(10)、MOY2(1)120.1820.00016東江(DON)河源10DON1(1)、DON2(1)、DON3(5)、DON4-HEJ1(3)*、2940.7110.00532北江(BEI)連州9BEI1(1)、BEI2(10)、BEI3(1)、BEI4(1)、1750.5060.00224始興5NAN1-MOY1-BEI5(1)賀江(HEJ)賀州11DON4-HEJ1(4)、HEJ2(7)120.5090.00045桂江(GUI)桂林6GUI1-LIU1-XIA1(4)*、GUI2-XIA2(1)*、GUI3(2)、3050.8060.00614荔浦3GUI4-LIU2(1)*、GUI5(1)柳江(LIU)從江3GUI1-LIU1-XIA1(3)、GUI4-LIU2(2)、LIU3(4)430.6670.00137永福6唇魚骨H.labeo韓江(HAN)上杭10HAN1(8)、HAN2(1)、HAN3(1)330.3780.00053九龍江(JIU)漳平5JIU1(2)、JIU2(3)120.6000.03072閩江(MIN)邵武9MIN1(2)、MIN2(1)、MIN3(1)、MIN4(2)、MIN5(3)4250.8610.01652湘江(XIA)全州9GUI1-LIU1-XIA1(6)、GUI2-XIA2(3)2720.5000.01184沅江(YUA)錦屏7YUA1(7)010.0000.00000信江(XIN)上饒3XIN1(2)、XIN2(1)120.6670.00058甌江(OUJ)麗水9OUJI(4)、OUJ2(5)220.5000.00088錢塘江(QIA)開(kāi)化4QIA1(3)、QIA2(1)120.5000.00044合計(jì)Total148128410.9540.02153

    ﹡表示此單倍型為不同河流共享單倍型

    表2 不同水系的唇魚骨與間魚骨Cyt b基因遺傳距離(對(duì)角線下方)和遺傳分化指數(shù)(% 對(duì)角線上方)

    圖1 唇魚骨與間魚骨中國(guó)南方種群基于Cyt b基因的最大鄰接NJ樹(shù)Fig.1 Neighbor-joining tree of H.labeo and H.medius populations of south China based on Cyt b gene

    圖2 唇魚骨與間魚骨中國(guó)南方種群基于Cyt b基因的單倍型網(wǎng)絡(luò)圖Fig.2 Minimum spanning network of Cyt b sequences of H.labeo and H.medius populations of south China

    單倍型網(wǎng)絡(luò)圖還揭示,桂江及柳江水系種群與長(zhǎng)江水系支流湘江水系種群的親緣關(guān)系十分密切,兩者具有共享單倍型且各個(gè)單倍型之間的變異不大;海南島三大水系的種群與漠陽(yáng)江種群的親緣關(guān)系比較密切;錢塘江、甌江及閩江種群的親緣關(guān)系亦較為密切;韓江與九龍江種群的親緣關(guān)系較近;而我國(guó)東南部水系閩江的種群不僅與長(zhǎng)江水系支流湘江、沅江的關(guān)系較近,且與珠江水系支流東江、長(zhǎng)江水系支流信江種群的關(guān)系亦較密切。這與NJ樹(shù)的結(jié)果相似。從單倍型網(wǎng)絡(luò)圖推測(cè),西江可能為珠江水系、漠陽(yáng)江水系和海南島水系間魚骨種群的擴(kuò)散中心,其中一支往東向北江和東江擴(kuò)散;另一支往南向海南島擴(kuò)散,海南島種群爾后在第四紀(jì)冰期海退時(shí)向北擴(kuò)散至漠陽(yáng)江水系及珠江水系。

    2.6 分子變異分析(AMOVA)

    為了解唇魚骨和間魚骨種群分子變異的分布模式,根據(jù)系統(tǒng)發(fā)育分析及單倍型網(wǎng)絡(luò)分析的結(jié)果,將唇魚骨和間魚骨17個(gè)水系的種群分成5個(gè)地理區(qū),進(jìn)行分子變異分析(AMOVA)。地理區(qū)的劃分主要根據(jù)各種群親緣關(guān)系的親疏及地理位置,其中海南島的南渡江、萬(wàn)泉河、昌化江和漠陽(yáng)江為一地理區(qū);珠江水系為一地理區(qū),包括西江支流柳江、桂江、賀江、珠江另外兩條支流北江及東江;長(zhǎng)江水系為一地理區(qū),包括湘江、沅江及信江;韓江和九龍江為一地理區(qū);閩江、甌江和錢塘江為一個(gè)地理區(qū)。AMOVA分析分析表明,唇魚骨和間魚骨地理區(qū)之間變異約占54.50%,地理區(qū)內(nèi)種群間約占18.64%,種群內(nèi)的變異占26.86%(表3)。這說(shuō)明,唇魚骨和間魚骨種群Cytb的遺傳分化主要是來(lái)自地理區(qū)之間。

    2.7 種群動(dòng)態(tài)分析

    除了昌化江、桂江、北江、東江、閩江等種群外,其他種群的單倍型數(shù)目未達(dá)到4,故無(wú)法計(jì)算Fu and Li′s D和Tajima′s D值。對(duì)唇魚骨和間魚骨的全部種群、唇魚骨種群和間魚骨種群進(jìn)行錯(cuò)配分析發(fā)現(xiàn),全部種群、唇魚骨及間魚骨種群均呈現(xiàn)多峰分布,說(shuō)明其均沒(méi)有發(fā)生過(guò)明顯的擴(kuò)張。而中性檢驗(yàn)也顯示,全部種群、唇魚骨及間魚骨種群的Fu and Li′s D和Tajima′s D值均沒(méi)有呈現(xiàn)顯著的負(fù)值,表明其在歷史上均沒(méi)有發(fā)生過(guò)明顯的擴(kuò)張(表4)。表4還表明,兩個(gè)支系的Fu and Li′s D和Tajima′s D值均沒(méi)有呈現(xiàn)顯著的負(fù)值,這說(shuō)明兩個(gè)支系亦未發(fā)生過(guò)擴(kuò)張。

    表3 唇魚骨和間魚骨分子變異分析

    表4 唇魚骨與間魚骨種群基于Cyt b基因的Tajima′s D 和 Fu and Li′s D檢驗(yàn)

    就各個(gè)水系的種群而言,桂江、北江、東江、閩江等種群的Fu and Li′s D和Tajima′s D均不為顯著的負(fù)值,說(shuō)明這幾個(gè)種群亦相對(duì)穩(wěn)定,沒(méi)有發(fā)生過(guò)擴(kuò)張;而昌化江(CHJ)的Fu and Li′s D為顯著的負(fù)值,而Tajima′s D則不顯著,但其錯(cuò)配分析圖為單峰分布(圖3),所以可以推測(cè)昌化江種群曾發(fā)生過(guò)種群擴(kuò)張。

    圖3 昌化江種群的歧點(diǎn)分布分析圖(其中實(shí)線代表期望值,虛線代表觀察值)Fig.3 Mismatch-distribution analysis of H.medius of Changhuajiang population of mtDNA haplotype sequences(solid lines represent the distribution expeted under constant population size, and dotted lines represent the observe distribution)

    3 討論

    3.1 間魚骨物種有效性

    在本研究中,遺傳距離分析表明,唇魚骨和間魚骨之間的遺傳距離為3.2%,低于唇魚骨種內(nèi)各種群間的最大遺傳距離(韓江與信江,5.9%)。系統(tǒng)分析顯示,間魚骨并未形成單系群,而是與部分唇魚骨種群組成一支,而韓江和九龍江的唇魚骨組成了另一支。單倍型分布分析顯示,湘江唇魚骨種群的所有單倍型均與桂江或柳江的間魚骨種群共享,湘江、桂江和柳江種群的遺傳差異極小。系統(tǒng)發(fā)育分析及單倍型網(wǎng)絡(luò)圖還顯示,除了唇魚骨的湘江、沅江種群與間魚骨的桂江、柳江種群的親緣關(guān)系較近外,唇魚骨的信江和閩江種群與間魚骨的東江種群亦有著較近的親緣關(guān)系。此外,藍(lán)昭軍基于多變量形態(tài)度量學(xué)的研究亦表明,唇魚骨和間魚骨存在一定的形態(tài)差異,但各個(gè)種群間存在梯度變異,差異并不顯著,無(wú)法截然分開(kāi)[23];而林龍峰基于線粒體CO I和ND5基因分別對(duì)中國(guó)南方唇魚骨與間魚骨的種群遺傳分化進(jìn)行了研究,結(jié)果均表明間魚骨與唇魚骨之間的遺傳差異較小,兩者存在共享單倍型,且間魚骨不能單獨(dú)構(gòu)成單系群,而是與唇魚骨部分樣本構(gòu)成一個(gè)分支[24]。因此,綜合考慮基于線粒體Cytb、CO I和ND5基因和多變量形態(tài)度量學(xué)的分析結(jié)果,我們認(rèn)為唇魚骨和間魚骨之間的遺傳差異未達(dá)到種級(jí)水平,間魚骨應(yīng)為唇魚骨的次定同物異名。

    3.2 韓江、九龍江種群與其他種群間的遺傳分化

    系統(tǒng)發(fā)育分析顯示,在我國(guó)南方的唇魚骨與間魚骨種群間,可以分成Ⅰ和Ⅱ兩支。韓江和九龍江種群遺傳差異不大,兩條江的全部樣本組成了Ⅰ支,而其他種群的樣本則組成了Ⅱ支,兩支系之間的遺傳距離為5.1%,這說(shuō)明韓江和九龍江種群已和其他種群發(fā)生了一定程度的遺傳分化。林弘都[25]和梁曉旭等[26]對(duì)部分淡水魚類的研究亦表明,華南大陸東部的蓮花山脈是粵東地區(qū)的一個(gè)重要地理屏障,而慶寧等[27]對(duì)華南地區(qū)黑眶蟾蜍(BufoMelanostictus)的研究亦得出了相似的結(jié)果。

    韓江是廣東省東部地區(qū)的最大河流,其有東西兩源,東源稱梅江,發(fā)源于蓮花山山脈北麓;西源稱汀江,發(fā)源于武夷山南麓。蓮花山脈位于廣東省東北部,其山勢(shì)由粵北部的大埔向西南一直延伸到香港附近,最高峰為銅鼓障,高1559 m; 整個(gè)山脈高度基本在980—1 170 m之間[28]。九龍江為福建省第二大水系,發(fā)源于福建玳瑁山南麓。由此可知,韓江和九龍江被蓮花山、武夷山、玳瑁山等較高的山脈圍繞,將其與珠江水系、長(zhǎng)江水系及閩江水系隔離。據(jù)此,我們推測(cè):一方面,由于蓮花山、武夷山、玳瑁山等山脈的隔離,韓江和九龍江的種群無(wú)法與其他種群進(jìn)行基因交流,因而產(chǎn)生了相對(duì)較大的遺傳分化;另一方面,韓江和九龍江在冰期海退時(shí),可能曾發(fā)生過(guò)連接,使得兩種群之間得以基因交流,因而遺傳差異不顯著。

    3.3 珠江及長(zhǎng)江種群間的遺傳分化

    分子系統(tǒng)發(fā)育分析顯示,除湘江與桂江和柳江種群的親緣關(guān)系較近以及東江少數(shù)樣本與信江、閩江種群的親緣關(guān)系較近外,唇魚骨和間魚骨剩下的絕大部分單倍型則分別歸入了不同的分支之內(nèi)。這表明,在唇魚骨和間魚骨分布區(qū)之間,可能存在地理障壁,造成其不同的分布范圍,而這一地理障壁可能為長(zhǎng)江水系和珠江水系的分水嶺之一——南嶺山脈。南嶺山脈由越城嶺、都龐嶺、萌渚嶺、騎田嶺和大庾嶺5條主要山嶺所組成。其橫亙?cè)谙婀?、湘粵、贛粵之間,向東延伸至閩南,東西長(zhǎng)約600km,南北寬約200km。廣義的南嶺還包括苗兒山、海洋山、九嶷山、香花嶺、瑤山、九連山等。有研究表明,南嶺山脈在11.06—8.04百萬(wàn)年前已經(jīng)形成[29],但在隆起初期,其高度相對(duì)較低,可能不足以阻隔兩側(cè)的淡水魚類的基因交流。而在青藏高原隆起的時(shí)期,南嶺山脈也呈加速隆起之勢(shì)[25]。南嶺的加速隆起,可能使其成為了南嶺以北水系與南嶺以南水系間的地理障壁,因而限制了南北兩側(cè)水系種群的基因交流,使得其兩側(cè)的種群產(chǎn)生了遺傳分化,進(jìn)而形態(tài)上也發(fā)生了一定程度的分化。

    在本研究中,多種分析結(jié)果亦顯示,長(zhǎng)江水系支流湘江、沅江及信江的唇魚骨種群與珠江水系間魚骨種群的親緣關(guān)系較近,其中湘江與沅江的種群在系統(tǒng)發(fā)育樹(shù)中與桂江和柳江的親緣關(guān)系均較近;而東江的部分樣本則與信江種群的親緣關(guān)系較近。在Li等對(duì)馬口魚(Opsariichthysbidens)的研究[30]及Yang等對(duì)大鰭鳠(Hemibagrusmacropterus)的研究[31]中,其結(jié)果均顯示珠江水系的種群與長(zhǎng)江水系的種群親緣關(guān)系較近。

    我們推測(cè),造成這種情況的原因可能有二,一是由于地殼運(yùn)動(dòng)、洪水泛濫等原因,使長(zhǎng)江水系和珠江水系在歷史上曾發(fā)生過(guò)襲奪事件或連接事件,使長(zhǎng)江南部支流上游的部分支流被珠江襲奪或與珠江水系支流相連接,從而使兩水系種群得以基因交流;二是水利樞紐工程的影響,在湘江水系與桂江水系上游,有一條建成已有2200多年人工運(yùn)河——靈渠將兩江聯(lián)通[32],這也可能促進(jìn)長(zhǎng)江水系與珠江水系種群的基因交流。

    3.4 海南島種群的起源、擴(kuò)散及隔離

    從唇魚骨和間魚骨單倍型的分布可知,海南島的南渡江水系、廣東西部水系漠陽(yáng)江及北江水系共享一個(gè)單倍型,表明三種群間親緣關(guān)系較近。在單倍型網(wǎng)絡(luò)圖中,單倍型WAN1與桂江、柳江和湘江的共享單倍型GUI1-LIU1-XIA1及桂江的單倍型GUI3之間的變異均較小,約為10步,表明這些單倍型間有著較近的親緣關(guān)系;萬(wàn)泉河的單倍型WAN1與南渡江、漠陽(yáng)江和北江的共享單倍型NAN1-MOY1-BEI5只有3步,兩者之間亦有著很近的親緣關(guān)系。此外,單倍型網(wǎng)絡(luò)圖還顯示,海南島中部山脈以東的南渡江和萬(wàn)泉河種群與西邊的昌化江種群之間關(guān)系相對(duì)較遠(yuǎn)(16步)。

    有關(guān)華南沿海地質(zhì)構(gòu)造的研究顯示,在第四紀(jì)時(shí),瓊州海峽至少經(jīng)歷過(guò)4次海侵與海退[33];而在第四紀(jì)冰期,由于海平面下降,華南沿海的海岸線在海南島以南,海南島北部均為遼闊的濱海平原[34];在更新世,在海南島和越南之間是一片陸地,其間有一條聯(lián)系西江、廣西、越南沿海和海南島諸水系的古河道[35]。有鑒于此,結(jié)合本研究系統(tǒng)發(fā)育分析及單倍型網(wǎng)絡(luò)圖等研究結(jié)果,海南島的間魚骨種群有可能為西江水系的種群沿著上述這一古河道擴(kuò)散而來(lái),然后進(jìn)一步在海南島內(nèi)擴(kuò)散。隨后,由于海平面的上升,海南島與大陸隔離,因而使海南島種群與大陸種群的基因交流受到了限制,發(fā)生了遺傳分化。而到了晚更新世冰期時(shí),由于海退等原因,雷州半島與海南島連成一片山地丘陵[36],這使得海南島的種群得以向北擴(kuò)散,到達(dá)雷州半島,進(jìn)而擴(kuò)散至漠陽(yáng)江水系,并重新回到珠江水系。間魚骨的這一擴(kuò)散途徑,與海南紋胸鮡(Glyptothoraxhainanensis)[37]及中間黃顙魚(Pelteobagrusintermedius)[38]有相似之處。

    在林宗鍵和Lin等對(duì)唇魚骨的親緣地理研究中,基于線粒體控制區(qū)序列的分析結(jié)果顯示,海南島昌化江種群(實(shí)為間魚骨)與長(zhǎng)江中游種群的親緣關(guān)系較其與廣東東部水系榕江種群的關(guān)系近,而昌化江與長(zhǎng)江的距離遠(yuǎn)大于其與榕江的距離,他們認(rèn)為其原因是長(zhǎng)江和海南島的水系曾通過(guò)古紅河水系發(fā)生過(guò)連接[17- 18]。值得指出的是,在Lin等的研究中,并沒(méi)有采集珠江水系的樣本,而海南島也只有昌化江水系的樣本。在本研究中,也發(fā)現(xiàn)海南島三大水系種群與長(zhǎng)江水系種群的親緣關(guān)系近于其與韓江和九龍江水系種群的親緣關(guān)系。由于韓江與榕江的地理位置很近,故此該結(jié)果與Lin等的相似??梢哉J(rèn)為,造成本研究及Lin等所述之結(jié)果的另一重要原因是海南島的種群由西江水系種群沿著古河道擴(kuò)散而來(lái),而西江曾與長(zhǎng)江水系發(fā)生過(guò)襲奪或連接事件,其種群與長(zhǎng)江水系種群親緣關(guān)系較近。因此,這一現(xiàn)象進(jìn)一步印證了我們對(duì)海南島種群來(lái)源及擴(kuò)散之推測(cè)。

    另一方面,在海南島內(nèi)部,單倍型網(wǎng)絡(luò)圖顯示昌化江種群與萬(wàn)泉河及南渡江種群的變異達(dá)到16步,說(shuō)明昌化江種群與其他兩種群之間,發(fā)生了一定的遺傳分化。地質(zhì)研究表明,在第三紀(jì)末期,海南島中部山地已經(jīng)形成,在多雨山區(qū)發(fā)源的河流也四散分流,獨(dú)自入海,不相聯(lián)系匯合[39]。在間魚骨擴(kuò)散至海南島各主要水系后,可能由于海南島中部山地,如五指山等的隆起,將海南島西部的昌化江種群與東部的萬(wàn)泉河及南渡江種群隔離,進(jìn)而使昌化江種群產(chǎn)生了分化。而左艷玲等基于線粒體DNA控制區(qū)序列對(duì)擬平鰍(Liniparhomalopteradisparis)的研究也表明:擬平鰍的昌化江種群與南渡江和萬(wàn)泉河的種群亦存在較大的遺傳分化[40]。因此,本研究的結(jié)果進(jìn)一步說(shuō)明了海南島中部山脈可能對(duì)海南島的淡水魚類存在著較為廣泛的隔離。

    致謝:華南師范大學(xué)陳湘粦教授和慶寧教授幫助實(shí)驗(yàn)方案設(shè)計(jì)及修改,林弘都博士和張雄同學(xué)參與了部分樣品的采集,Radhakrishnan博士幫助寫作,特此致謝。

    [1] Hewitt G. The genetic legacy of the quaternary ice ages. Nature, 2000, 405(6789): 907- 913.

    [2] Angers B, Magnan P, Plante M, Bernatchez L. Canonical correspondence analysis for estimating spatial and environmental effects on microsatellite gene diversity in brook charr (Salvelinusfontinalis). Molecular Ecology, 1999, 8(6): 1043- 1053.

    [3] Keyghobadi N, Roland J, Strobeck C. Influence of landscape on the population genetic structure of the alpine butterflyParnassiussmintheus(Papilionidae). Molecular Ecology, 1999, 8(9): 1481- 1495.

    [4] Castric V, Bonney F, Bernatchez L. Landscape structure and hierarchical genetic diversity in the brook charr,Salvelinusfontinalis. Evolution, 2001, 55(5): 1016 - 1028.

    [5] Cassel A, Tammaru T. Allozyme variability in central, peripheral and isolated populations of the scarce heath (Coenonymphahero: Lepidoptera, Nymphalidae); implications for conservation. Conservation Genetics, 2003, 4(1): 83- 93.

    [6] Costello A B, Down T E, Pollard S M, Pacas C J, Taylor E B. The influence of history and contemporary stream hydrology on the evolution of genetic diversity within species: an examination of microsatellite DNA variation in bull trout,Salvelinusconfluentus(Pisces: Salmonidae). Evolution, 2003, 57(2): 328- 344.

    [7] Hewitt G M. Some genetic consequences of ice ages, and their role in divergence and speciation. Biological Journal of the Linnean Society, 1996, 58: 247- 276.

    [8] 樂(lè)佩琦. 魚骨和間魚骨屬魚類的分類整理(鯉形目: 鯉科). 動(dòng)物分類學(xué)報(bào), 1995, 20(1): 116- 123.

    [9] 陳宜瑜. 中國(guó)動(dòng)物志硬骨魚綱鯉形目(中卷). 北京: 科學(xué)出版社, 1998: 239- 241.

    [10] 伍獻(xiàn)文. 中國(guó)鯉科魚類志(下卷). 上海: 上海人民出版社, 1977: 439- 537.

    [11] 廣西壯族自治區(qū)水產(chǎn)研究所, 中國(guó)科學(xué)院動(dòng)物研究所. 廣西淡水魚類志. 南寧: 廣西人民出版社, 1981: 112- 113.

    [12] 中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所等. 海南島淡水及河口魚類志. 廣州: 廣東科學(xué)技術(shù)出版社, 1986: 93- 94.

    [13] 鄭慈英. 珠江魚類志. 北京: 科學(xué)出版社, 1989: 119- 121.

    [14] 中國(guó)水產(chǎn)科學(xué)研究院珠江水產(chǎn)研究所等. 廣東淡水魚類志. 廣州: 廣東科技出版社, 1991.

    [15] Yang J Q, He S P, Freyhof J, Witte K, Liu H Z. The Phylogenetic Relationships of the Gobioninae (Teleostei: Cyprinidae) Inferred from Mitochondrial CytochromebGene Sequences. Hydrobiologia, 2006, 553(1): 255- 266.

    [16] 廣西壯族自治區(qū)水產(chǎn)研究所, 中國(guó)科學(xué)院動(dòng)物研究所. 廣西淡水魚類志(第二版). 南寧: 廣西人民出版社, 2006, 193- 195.

    [17] 林宗鍵. 唇魚骨和間魚骨的親緣地理研究. 臺(tái)灣:國(guó)立成功大學(xué), 2009.

    [18] Lin C J, Lin H D, Wang J P, Chao S C, Chiang T Y. Phylogeography ofHemibarbuslabeo(Cyprinidae): secondary contact of ancient lineages of mtDNA. Zoologica Scripta, 2010, 39(1): 23- 35.

    [19] 慶寧, 丘城鋒, 廖偉群, 馬天峰, 梁曉旭, 列金妮. 華南沿海西部美麗小條鰍基于線粒體控制區(qū)的種群遺傳變異及親緣地理格局. 生態(tài)學(xué)報(bào), 2010, 30(1): 258- 264.

    [20] Kumar S, Nei M, Dudley J, Tamura K. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Briefings in Bioinformatics, 2008, 9(4): 299- 306.

    [21] Librado P, Rozas J. DNASP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 2009, 25(11): 1451- 1452.

    [22] Excoffier L, Smouse P E, Quattor J M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics, 1992, 131(2): 479- 491.

    [23] 藍(lán)昭軍. 中國(guó)南方唇魚骨和間魚骨的種群分化、親緣地理及間魚骨的生物學(xué)特征[D]. 廣州: 華南師范大學(xué), 2011.

    [24] 林龍峰. 基于COI和ND5序列的中國(guó)南方唇魚骨和間魚骨種群遺傳變異與親緣地理分析[D]. 廣州: 華南師范大學(xué), 2014.

    [25] 林弘都. 臺(tái)灣與中國(guó)大陸地區(qū)鯉科魚類之親緣地理研究[D]. 臺(tái)灣: 國(guó)立成功大學(xué), 2008.

    [26] 梁曉旭, 慶寧, 楊柯林, 萬(wàn)彩霞, 趙俊, 陳湘粦. 廣東地區(qū)寬鰭鱲種群遺傳變異和親緣地理. 水生生物學(xué)報(bào), 2010, 34(4): 806- 814.

    [27] 慶寧, 馬天峰, 梁曉旭, 林弘都, 盧文華, 鹽司櫓. 華南地區(qū)黒眶蟾蜍的遺傳變異和地理分化. 動(dòng)物分類學(xué)報(bào), 2011, 36(2): 356- 367.

    [28] 曾昭璇, 黃偉峰. 廣東自然地理. 廣州: 廣東人民出版社, 2001.

    [29] 地質(zhì)礦產(chǎn)部《南嶺項(xiàng)目》構(gòu)造專題組. 南嶺區(qū)域構(gòu)造特征及控巖控礦構(gòu)造研究. 北京: 地質(zhì)出版社, 1988.

    [30] Li G Y, Wang X Z, Zhao Y H, Zhang J, Zhang C G, He S P. Speciation and phylogeography ofOpsariichthysbidens(Pisces: Cypriniformes: Cyprinidae) in China: analysis of the cytochrome b gene of mtDNA from diverse populations. Zoological Studies, 2009, 48(4): 569- 583.

    [31] Yang L, Mayden R L, He S P. Population genetic structure and geographical differentiation of the Chinese catfishHemibagrusmacropterus(Siluriformes, Bagridae): evidence for altered drainage patterns. Molecular Phylogenetics and Evolution, 2009, 51(2): 405- 411.

    [32] 宗鳳書. 靈渠地質(zhì)地貌初步研究. 中國(guó)巖溶, 1990, 9(1): 88- 93.

    [33] 趙煥庭, 王麗榮, 袁家義. 瓊州海峽成因與時(shí)代. 海洋地質(zhì)與第四紀(jì)地質(zhì), 2007, 27(2): 33- 40.

    [34] 張虎男, 陳偉光. 華南沿海新構(gòu)造運(yùn)動(dòng)與地質(zhì)環(huán)境. 北京: 地震出版社, 1990.

    [35] Harold K V. Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations. Journal of Biogeography, 2000, 27(5): 1153- 1167.

    [36] 中國(guó)科學(xué)院南海海洋研究所海洋地質(zhì)研究室. 華南沿海第四紀(jì)地質(zhì). 北京: 科學(xué)出版社, 1978.

    [37] Chen X L, Chiang T Y, Lin H D, Zheng H S, Shao K T, Zhang Q, Hsu K C. Mitochondrial DNA phylogeography ofGlyptothoraxfokiensisandGlyptothoraxhainanensisin Asia. Journal of Fish Biology, 2007, 70(S): 75- 93.

    [38] 慶寧, 呂鳳義, 廖偉群, 趙俊, 丘城鋒. 中間黃顙魚群體遺傳變異與親緣生物地理. 動(dòng)物學(xué)報(bào), 2007, 53(5): 845- 852.

    [39] 曾昭璇, 曾憲中. 海南島自然地理. 北京: 科學(xué)出版社, 1989.

    [40] 左艷玲, 林岳光, 梁曉旭, 馬天峰, 慶寧. 基于mtDNA控制區(qū)序列的擬平鰍遺傳變異和種群分化. 水產(chǎn)學(xué)報(bào), 2009, 33(6): 925- 931.

    Population diversity and phylogeography ofHemibarbuslabeoandHemibarbusmediusin South China

    LAN Zhaojun1,2, FAN Mingjun1,3, HUANG Xiaolin1,4, ZHAO Jun1,*

    1GuangdongProvincialKeyLaboratoryforHealthyandSafeAquaculture,KeyLaboratoryofEcologyandEnvironmentScienceinGuangdongHigherEducation,GuangzhouKeyLaboratoryofSubtropicalBiodiversityandBiomonitor,SchoolofLifeScience,SouthChinaNormalUniversity,Guangzhou510631,China2ShaoguanFisheriesResearchInstitute,Shaoguan512006,China3ShenzhenBGIFisheriesScience&Technologycorporation,Ltd,Shenzhen518083,China4SouthChinaSeaFisheriesResearchInstitute,ChineseAcademyofFisherySciences,Guangzhou510300,China

    HemibarbuslabeoandHemibarbusmedius(Cypriniformes: Cyprinidae: Gobioninae:Hemibarbus) are primary freshwater fish species, and their wide distribution makes them ideal models for phylogeographical studies. In this study, variations in the population genetics and phylogeographical patterns of 148 specimens belonging to 8 populations ofH.labeoand 9 populations ofH.mediuscollected from 17 drainage systems in south China were investigated using the nucleotide sequences of the mtDNA Cytochromebgene (1140 bp). The results reveal 128 variable sites (11.23%) within the gene sequences ofH.labeoandH.medius. A total of 41 haplotypes were identified, with the haplotype diversity (h) and nucleotide diversity (π) being 0.954 and 0.02153, respectively. This indicates a high level of genetic diversity and evolutionary potential in both species. The results of the neighbor-joining tree demonstrate thatH.labeoandH.mediusindividuals fall into two major clades (clade I and clade II): clade I is composed of all specimens collected from the Hanjiang and Jiulongjiang Rivers, all of which were identified asH.labeo. Meanwhile, all remaining populations fell into clade II, which includes specimens of bothH.labeoandH.medius. The genetic distance between clade I and clade II was 5.1%, while that betweenH.labeoandH.mediuswas 3.2%. Our results indicate that the specimens collected from these sampling localities may not be differentiated into two or more subspecies. The haplotype network indicates that the populations of the Hanjiang and Jiulongjiang Rivers exhibit a relatively high level of genetic variation compared to that of the rest of the rivers in the region, and that the populations from Hainan Island and Moyangjiang River were genetically close to those from the Pearl River System, but not to those from the Yangtze River System,Hanjiang River, Jiulongjiang River, Minjiang River, Oujiang River and Qiantangjiang River. TheH.mediuspopulations of south China seem to have originated from the Xijiang River and reached south China by one of two routes: one route may be through the Xijiang River to the Beijiang and Dongjiang rivers; and the other may be though the southern Guangxi drainages to Hainan Island, and then to the Moyangjiang and Beijiang rivers through the Qiongzhou Strait. According to the neighbor-joining tree and the haplotype network, five regions of genetic distribution were defined: Hainan Island and the Moyangjiang River region; the Pearl River region; the Yangtze River region; the Hanjiang and Jiulongjiang river region; and the Minjiang, Oujiang, and Qiantangjiang river region. An analysis of molecular variance (AMOVA) showed that the genetic variation of populations among regions was 54.50%; conversely, the genetic variation among populations within regions was 18.64%, and that within the populations of different regions was 26.86%. This indicates that the greatest genetic variation is found among the populations of the different regions. Mismatch distribution and tests of neutrality taking populations ofH.labeoandH.mediusinto account both together and separately were all multimodal types, and the value of Fu, Li′s D, and Tajima′s D for all comparisons were all negative but non-significant or positive. The results reveal relatively stable populations ofH.labeoandH.medius. Mismatch distribution and tests of neutrality also showed that both clade I and clade II were relatively stable; however, the population of the Changhuajiang River has undergone an obvious population expansion.

    Hemibarbuslabeo;Hemibarbusmedius; population diversity; phylogeography;Cytochromebgene; South China

    國(guó)家自然科學(xué)基金項(xiàng)目(31372178);公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)經(jīng)費(fèi)資助(201303048);廣東省科技計(jì)劃項(xiàng)目(2015A070706009)

    2014- 11- 04;

    日期:2016- 01- 15

    10.5846/stxb201411042169

    *通訊作者Corresponding author.E-mail: zhaojun@scnu.edu.cn

    藍(lán)昭軍,范明君,黃小林,趙俊.基于線粒體Cytb基因的中國(guó)南方唇魚骨(Hermibarbuslabeo)和間魚骨(Hermibarbusmedius)種群分化及親緣地理研究.生態(tài)學(xué)報(bào),2016,36(19):6091- 6102.

    Lan Z J, Fan M J, Huang X L, Zhao J.Population diversity and phylogeography ofHemibarbuslabeoandHemibarbusmediusin South China.Acta Ecologica Sinica,2016,36(19):6091- 6102.

    猜你喜歡
    魚骨親緣海南島
    谷子近緣野生種的親緣關(guān)系及其利用研究
    中國(guó)醫(yī)學(xué)科學(xué)院藥用植物研究所藥用植物親緣學(xué)研究中心
    奶奶愛(ài)拼魚骨畫
    菊科藥用植物遺傳多樣性及親緣關(guān)系的ISSR分析
    一星期沒(méi)換水的夢(mèng)境
    詩(shī)潮(2018年3期)2018-03-26 12:29:30
    小白菜種質(zhì)遺傳多樣性與親緣關(guān)系的SRAP 和SSR分析
    魚骨千萬(wàn)別丟 它能幫你增壽
    工友(2016年4期)2016-09-18 05:57:49
    是海南省還是海南島?
    “魚骨圖”在項(xiàng)目教學(xué)中的應(yīng)用
    河南科技(2014年12期)2014-02-27 14:10:59
    關(guān)于海南島戰(zhàn)役作戰(zhàn)方針的考察
    軍事歷史(1994年6期)1994-08-15 08:56:38
    亚洲精品久久国产高清桃花| 久久久久久久午夜电影| 老司机福利观看| 欧美又色又爽又黄视频| 人人妻人人澡欧美一区二区| 久久精品91蜜桃| 国产午夜精品论理片| 午夜日韩欧美国产| 亚洲中文日韩欧美视频| 亚洲三级黄色毛片| 可以在线观看毛片的网站| 亚洲中文日韩欧美视频| 搡老熟女国产l中国老女人| 国产精品人妻久久久久久| 可以在线观看毛片的网站| 天堂av国产一区二区熟女人妻| 国产亚洲精品综合一区在线观看| 欧美日本视频| 成人永久免费在线观看视频| h日本视频在线播放| 老熟妇乱子伦视频在线观看| 亚洲专区国产一区二区| 白带黄色成豆腐渣| 色综合色国产| 人妻制服诱惑在线中文字幕| 国产精品乱码一区二三区的特点| 亚洲无线观看免费| 欧美+日韩+精品| 精品福利观看| 国产极品精品免费视频能看的| 成年女人看的毛片在线观看| 国产免费一级a男人的天堂| 2021天堂中文幕一二区在线观| 97人妻精品一区二区三区麻豆| 最新中文字幕久久久久| 草草在线视频免费看| 色综合亚洲欧美另类图片| 免费电影在线观看免费观看| 深夜精品福利| 久久久久久久午夜电影| 一个人免费在线观看电影| 日韩av在线大香蕉| 日日摸夜夜添夜夜爱| 高清毛片免费观看视频网站| 国产精品嫩草影院av在线观看| 午夜老司机福利剧场| 在线看三级毛片| 成年av动漫网址| 精品99又大又爽又粗少妇毛片| 欧美高清性xxxxhd video| 国产 一区 欧美 日韩| 啦啦啦韩国在线观看视频| 国产激情偷乱视频一区二区| 亚洲无线观看免费| 欧美日韩综合久久久久久| 日韩一本色道免费dvd| 亚洲精品日韩av片在线观看| 亚洲熟妇熟女久久| 免费观看人在逋| 1024手机看黄色片| 亚洲美女搞黄在线观看 | 亚洲精品日韩在线中文字幕 | 精品一区二区免费观看| 哪里可以看免费的av片| 丰满乱子伦码专区| av在线播放精品| 日韩 亚洲 欧美在线| 成年女人毛片免费观看观看9| 一级毛片aaaaaa免费看小| 国产免费男女视频| av在线蜜桃| av国产免费在线观看| 免费电影在线观看免费观看| 国产一区二区三区在线臀色熟女| 精品久久久久久久末码| 在线看三级毛片| 国产精品一区二区免费欧美| 欧美人与善性xxx| 免费在线观看影片大全网站| 亚洲真实伦在线观看| 精品国产三级普通话版| 波多野结衣高清作品| 在线观看美女被高潮喷水网站| 久久午夜福利片| 久久精品国产亚洲网站| 国内久久婷婷六月综合欲色啪| 国产精品不卡视频一区二区| 在线观看66精品国产| 蜜臀久久99精品久久宅男| 高清日韩中文字幕在线| 嫩草影院精品99| 男人和女人高潮做爰伦理| 女人被狂操c到高潮| 别揉我奶头 嗯啊视频| 大又大粗又爽又黄少妇毛片口| 不卡视频在线观看欧美| 午夜福利在线观看免费完整高清在 | 成人午夜高清在线视频| 老司机影院成人| 最新在线观看一区二区三区| 亚洲av成人精品一区久久| 国产免费男女视频| 午夜福利在线观看免费完整高清在 | 99热精品在线国产| 中文资源天堂在线| a级一级毛片免费在线观看| 天天一区二区日本电影三级| 国产一级毛片七仙女欲春2| 久久国产乱子免费精品| 日韩三级伦理在线观看| 日韩成人av中文字幕在线观看 | 欧美区成人在线视频| 国产精品久久视频播放| 美女 人体艺术 gogo| 国产精品永久免费网站| 欧美高清成人免费视频www| 老熟妇乱子伦视频在线观看| 国产精品综合久久久久久久免费| 午夜a级毛片| 淫秽高清视频在线观看| 精品久久久久久久末码| 成年免费大片在线观看| 精品久久久久久久久久免费视频| 日韩,欧美,国产一区二区三区 | 国产69精品久久久久777片| 老女人水多毛片| 国产午夜精品论理片| 一区二区三区高清视频在线| 亚洲欧美成人综合另类久久久 | 欧美激情久久久久久爽电影| 成人毛片a级毛片在线播放| 久久精品国产自在天天线| 97在线视频观看| 俺也久久电影网| 成年女人毛片免费观看观看9| 精品国内亚洲2022精品成人| 免费无遮挡裸体视频| 欧美一区二区精品小视频在线| 在线观看一区二区三区| 亚洲内射少妇av| av卡一久久| 成人性生交大片免费视频hd| 色尼玛亚洲综合影院| 成人精品一区二区免费| 永久网站在线| 一a级毛片在线观看| 九色成人免费人妻av| 少妇猛男粗大的猛烈进出视频 | 欧美精品国产亚洲| 18禁在线播放成人免费| 久久久久国内视频| 亚洲av电影不卡..在线观看| 看片在线看免费视频| 亚洲专区国产一区二区| 91av网一区二区| 欧美一区二区国产精品久久精品| 男人的好看免费观看在线视频| 国产单亲对白刺激| 赤兔流量卡办理| 熟女电影av网| 在线观看一区二区三区| а√天堂www在线а√下载| 简卡轻食公司| 老女人水多毛片| 久久久久国产精品人妻aⅴ院| 夜夜看夜夜爽夜夜摸| 免费电影在线观看免费观看| 嫩草影院新地址| 韩国av在线不卡| 中文字幕av在线有码专区| 亚洲中文字幕日韩| 久久婷婷人人爽人人干人人爱| 免费看美女性在线毛片视频| 99在线视频只有这里精品首页| 亚洲图色成人| 久久久久国产精品人妻aⅴ院| 成人av一区二区三区在线看| 国产私拍福利视频在线观看| 91久久精品国产一区二区成人| 日本黄色片子视频| 国产欧美日韩精品一区二区| 欧美区成人在线视频| 少妇被粗大猛烈的视频| 少妇的逼好多水| 久久韩国三级中文字幕| 三级毛片av免费| 久久人妻av系列| 一级黄色大片毛片| 男女下面进入的视频免费午夜| 久久久久久大精品| 1024手机看黄色片| 国产在视频线在精品| 日韩国内少妇激情av| 波多野结衣高清无吗| 国产精品一区二区三区四区久久| 成人午夜高清在线视频| 国产视频内射| 一级毛片aaaaaa免费看小| 欧美三级亚洲精品| 精品乱码久久久久久99久播| 国产成人a∨麻豆精品| 最新在线观看一区二区三区| 午夜日韩欧美国产| 亚洲精品亚洲一区二区| 九色成人免费人妻av| 国产欧美日韩一区二区精品| 天堂网av新在线| 亚洲精品一卡2卡三卡4卡5卡| 久久国产乱子免费精品| 久久久久久国产a免费观看| 亚洲欧美日韩高清专用| 老师上课跳d突然被开到最大视频| 少妇高潮的动态图| 五月伊人婷婷丁香| 午夜亚洲福利在线播放| 久久人人爽人人片av| 日日摸夜夜添夜夜添小说| 国产av一区在线观看免费| 免费看光身美女| 亚洲欧美日韩无卡精品| 嫩草影院精品99| 99热全是精品| 久久久欧美国产精品| 男女啪啪激烈高潮av片| 少妇的逼好多水| 最近在线观看免费完整版| 欧美一级a爱片免费观看看| 亚洲第一区二区三区不卡| 国产精品一及| 午夜福利在线观看吧| 亚洲美女搞黄在线观看 | 女生性感内裤真人,穿戴方法视频| 午夜爱爱视频在线播放| 在线播放无遮挡| 精品人妻视频免费看| 亚洲性夜色夜夜综合| 别揉我奶头~嗯~啊~动态视频| 日韩大尺度精品在线看网址| 日韩欧美 国产精品| 九色成人免费人妻av| 国产精品一区www在线观看| 天堂网av新在线| 亚洲第一区二区三区不卡| 亚洲熟妇中文字幕五十中出| 最近在线观看免费完整版| 成人欧美大片| av在线蜜桃| 亚洲18禁久久av| 99久国产av精品国产电影| 亚洲在线观看片| 1000部很黄的大片| 麻豆国产97在线/欧美| 国产黄色视频一区二区在线观看 | 少妇猛男粗大的猛烈进出视频 | 91在线精品国自产拍蜜月| 精品久久久久久成人av| 国产毛片a区久久久久| 午夜a级毛片| 亚洲aⅴ乱码一区二区在线播放| 中出人妻视频一区二区| 国产精品乱码一区二三区的特点| av在线观看视频网站免费| 男女边吃奶边做爰视频| 国产中年淑女户外野战色| av专区在线播放| 日韩人妻高清精品专区| 激情 狠狠 欧美| 精品欧美国产一区二区三| 精品国内亚洲2022精品成人| 亚洲av中文av极速乱| 亚洲av第一区精品v没综合| 能在线免费观看的黄片| 亚洲三级黄色毛片| 91午夜精品亚洲一区二区三区| 成人特级av手机在线观看| 一进一出抽搐动态| 国产欧美日韩精品一区二区| 国产免费一级a男人的天堂| 亚洲精品粉嫩美女一区| 美女 人体艺术 gogo| 国产伦一二天堂av在线观看| 精品人妻视频免费看| 狠狠狠狠99中文字幕| 亚洲av电影不卡..在线观看| 日韩在线高清观看一区二区三区| 搡老岳熟女国产| 两个人视频免费观看高清| 日本三级黄在线观看| 亚洲成人av在线免费| 日韩制服骚丝袜av| 春色校园在线视频观看| 成人综合一区亚洲| 久久欧美精品欧美久久欧美| 男女边吃奶边做爰视频| 老司机午夜福利在线观看视频| 亚洲国产精品合色在线| 国产精品av视频在线免费观看| 六月丁香七月| 99热这里只有是精品50| 全区人妻精品视频| 长腿黑丝高跟| 99热6这里只有精品| 欧美日本视频| 97人妻精品一区二区三区麻豆| 免费人成视频x8x8入口观看| 国产在视频线在精品| 亚洲成人久久性| 美女xxoo啪啪120秒动态图| 成人二区视频| 女同久久另类99精品国产91| 十八禁网站免费在线| 日本撒尿小便嘘嘘汇集6| 欧美绝顶高潮抽搐喷水| 亚洲欧美成人综合另类久久久 | 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 精品日产1卡2卡| 欧美区成人在线视频| 色吧在线观看| 日韩 亚洲 欧美在线| 欧美激情国产日韩精品一区| 婷婷六月久久综合丁香| 亚洲va在线va天堂va国产| 日韩一区二区视频免费看| 久久人人精品亚洲av| 一区二区三区四区激情视频 | 国产久久久一区二区三区| 国产熟女欧美一区二区| 国产黄a三级三级三级人| 欧美色欧美亚洲另类二区| 国产探花极品一区二区| 国内少妇人妻偷人精品xxx网站| 永久网站在线| 高清午夜精品一区二区三区 | 乱系列少妇在线播放| 亚洲av免费高清在线观看| 91久久精品国产一区二区成人| 成人鲁丝片一二三区免费| 中国美白少妇内射xxxbb| 少妇丰满av| 给我免费播放毛片高清在线观看| 亚洲精品粉嫩美女一区| 精品无人区乱码1区二区| 国产毛片a区久久久久| 岛国在线免费视频观看| 99久国产av精品| av.在线天堂| 久久久久久久午夜电影| 晚上一个人看的免费电影| 熟女电影av网| or卡值多少钱| 校园春色视频在线观看| 免费观看在线日韩| 中国国产av一级| 国产91av在线免费观看| 婷婷精品国产亚洲av| 天堂动漫精品| 99热这里只有是精品在线观看| 成年版毛片免费区| 国产精品一区二区三区四区久久| 97超级碰碰碰精品色视频在线观看| 黄色配什么色好看| 精品午夜福利在线看| 美女高潮的动态| 在线天堂最新版资源| 欧美一区二区国产精品久久精品| 色av中文字幕| 日本黄大片高清| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美| 日韩精品青青久久久久久| 91久久精品国产一区二区成人| 精品久久久噜噜| 日本与韩国留学比较| videossex国产| 日韩成人伦理影院| 一区福利在线观看| 久久久国产成人免费| 欧美激情国产日韩精品一区| 国产精品免费一区二区三区在线| 久久久久性生活片| 三级经典国产精品| 国产一区二区亚洲精品在线观看| 特大巨黑吊av在线直播| 国产精品一区二区性色av| 91在线精品国自产拍蜜月| 成人一区二区视频在线观看| 久久精品久久久久久噜噜老黄 | 国产三级在线视频| 最后的刺客免费高清国语| 亚洲,欧美,日韩| 久久久久久久久久黄片| 国产亚洲欧美98| 日本-黄色视频高清免费观看| 69人妻影院| 看黄色毛片网站| 亚洲一级一片aⅴ在线观看| 老司机午夜福利在线观看视频| 国产高潮美女av| 少妇猛男粗大的猛烈进出视频 | 国产黄a三级三级三级人| 午夜福利成人在线免费观看| 国产成人精品久久久久久| 两个人视频免费观看高清| 日韩一区二区视频免费看| 婷婷六月久久综合丁香| 韩国av在线不卡| 插逼视频在线观看| 免费大片18禁| 99久久精品一区二区三区| 亚洲欧美成人综合另类久久久 | 亚洲av.av天堂| 日本撒尿小便嘘嘘汇集6| 亚洲国产高清在线一区二区三| 69人妻影院| 亚洲丝袜综合中文字幕| 久久精品国产亚洲av香蕉五月| 久久人人爽人人片av| 熟女电影av网| 一个人免费在线观看电影| 久久久精品大字幕| 看非洲黑人一级黄片| 精品无人区乱码1区二区| av.在线天堂| 卡戴珊不雅视频在线播放| 亚洲精华国产精华液的使用体验 | 国产亚洲91精品色在线| 如何舔出高潮| 97碰自拍视频| 男人的好看免费观看在线视频| 熟妇人妻久久中文字幕3abv| 精品久久久久久久久av| 尤物成人国产欧美一区二区三区| 亚洲自拍偷在线| 精品人妻熟女av久视频| 亚洲电影在线观看av| 日日摸夜夜添夜夜爱| 色5月婷婷丁香| 国产一区二区亚洲精品在线观看| 天堂网av新在线| 亚洲综合色惰| 91精品国产九色| 日日啪夜夜撸| 国产片特级美女逼逼视频| 亚洲精品一卡2卡三卡4卡5卡| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 热99在线观看视频| 97超级碰碰碰精品色视频在线观看| 1000部很黄的大片| 18禁裸乳无遮挡免费网站照片| 91狼人影院| 精品久久国产蜜桃| 亚洲人成网站高清观看| 一a级毛片在线观看| 久久99热6这里只有精品| 日韩欧美 国产精品| 午夜免费激情av| 老熟妇乱子伦视频在线观看| 一夜夜www| 亚洲欧美成人精品一区二区| 天堂√8在线中文| 久久人妻av系列| 久久这里只有精品中国| 一级av片app| 午夜福利18| 秋霞在线观看毛片| 国产成人精品久久久久久| 99九九线精品视频在线观看视频| 欧美性猛交黑人性爽| 亚洲欧美日韩高清在线视频| 精品一区二区免费观看| 女人被狂操c到高潮| 日韩一本色道免费dvd| 久久午夜亚洲精品久久| 免费av不卡在线播放| а√天堂www在线а√下载| 亚洲欧美成人精品一区二区| 99在线人妻在线中文字幕| 99热这里只有是精品50| 日韩欧美一区二区三区在线观看| 特大巨黑吊av在线直播| 国产男靠女视频免费网站| 亚洲成人久久爱视频| 色哟哟哟哟哟哟| 亚洲欧美日韩高清专用| 久久精品91蜜桃| 最近手机中文字幕大全| 婷婷精品国产亚洲av| 免费人成视频x8x8入口观看| 欧美高清成人免费视频www| 哪里可以看免费的av片| 国产综合懂色| www.色视频.com| 国产精品久久电影中文字幕| 国产真实伦视频高清在线观看| 精品无人区乱码1区二区| 久久久久久九九精品二区国产| 日本一二三区视频观看| 三级男女做爰猛烈吃奶摸视频| 成人午夜高清在线视频| 日本熟妇午夜| 在线观看美女被高潮喷水网站| 久久久欧美国产精品| 久久久久久久午夜电影| 99精品在免费线老司机午夜| 丰满人妻一区二区三区视频av| 又爽又黄a免费视频| 午夜a级毛片| 亚洲av不卡在线观看| 搡老岳熟女国产| 亚洲av电影不卡..在线观看| 美女大奶头视频| 性欧美人与动物交配| 久久人人精品亚洲av| 国产精品一二三区在线看| 国产男靠女视频免费网站| 秋霞在线观看毛片| 亚洲精品日韩在线中文字幕 | 国产中年淑女户外野战色| 日韩av不卡免费在线播放| 久99久视频精品免费| 久久久久久九九精品二区国产| 国产人妻一区二区三区在| 寂寞人妻少妇视频99o| 好男人在线观看高清免费视频| 免费观看在线日韩| 俺也久久电影网| 你懂的网址亚洲精品在线观看 | 日韩人妻高清精品专区| 国产精品一区www在线观看| 国产精品综合久久久久久久免费| 国产真实乱freesex| 国产免费一级a男人的天堂| 好男人在线观看高清免费视频| 一区二区三区高清视频在线| 国产精品一区www在线观看| 91久久精品电影网| 在线免费十八禁| 国内少妇人妻偷人精品xxx网站| 能在线免费观看的黄片| 大型黄色视频在线免费观看| 狂野欧美白嫩少妇大欣赏| 午夜福利18| 精品午夜福利视频在线观看一区| 久久精品国产自在天天线| 99在线视频只有这里精品首页| 男女啪啪激烈高潮av片| 99热6这里只有精品| av在线观看视频网站免费| 色视频www国产| 精品久久久久久久久亚洲| 欧洲精品卡2卡3卡4卡5卡区| 久久99热这里只有精品18| 天堂影院成人在线观看| 丰满人妻一区二区三区视频av| 亚洲在线自拍视频| 丰满的人妻完整版| 色视频www国产| 亚洲精品日韩在线中文字幕 | 精品人妻熟女av久视频| 亚洲国产高清在线一区二区三| 美女xxoo啪啪120秒动态图| 国产伦精品一区二区三区四那| 少妇的逼好多水| 色尼玛亚洲综合影院| 国产 一区精品| 免费黄网站久久成人精品| 精品一区二区免费观看| 国产高清三级在线| 长腿黑丝高跟| 成人精品一区二区免费| 中国国产av一级| 91狼人影院| 丝袜美腿在线中文| 变态另类成人亚洲欧美熟女| 一本久久中文字幕| 两个人的视频大全免费| 久99久视频精品免费| 国产精品久久久久久亚洲av鲁大| 麻豆av噜噜一区二区三区| 欧美bdsm另类| 中文字幕人妻熟人妻熟丝袜美| 天堂动漫精品| 在线国产一区二区在线| 嫩草影院新地址| 人人妻人人澡人人爽人人夜夜 | 国产成人a∨麻豆精品| 一级黄色大片毛片| 午夜精品在线福利| 久久欧美精品欧美久久欧美| 97超级碰碰碰精品色视频在线观看| 国产精品乱码一区二三区的特点| 免费av毛片视频| 免费黄网站久久成人精品| 日本 av在线| 午夜激情福利司机影院| 亚洲av成人av| 久久久久久久久久久丰满| 亚洲av免费在线观看| 99精品在免费线老司机午夜| 国产一区二区三区在线臀色熟女| 国产精品乱码一区二三区的特点| 一区二区三区免费毛片| 国内精品美女久久久久久| 中文字幕精品亚洲无线码一区| 日本色播在线视频| 尤物成人国产欧美一区二区三区| 哪里可以看免费的av片| 久久久久免费精品人妻一区二区| 69av精品久久久久久| 免费观看在线日韩| 非洲黑人性xxxx精品又粗又长| 日本色播在线视频| 午夜久久久久精精品|