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

    基于線粒體Cytb基因序列探討花斑鰍三亞種的分類地位

    2015-08-11 11:54:40陳詠霞梁娜李秀潘曉睿武大勇
    關鍵詞:花斑亞種西伯利亞

    陳詠霞,梁娜,李秀,潘曉睿,武大勇

    (1.河北大學生命科學學院,河北保定071002;2.衡水學院生命科學系,河北衡水053000)

    基于線粒體Cytb基因序列探討花斑鰍三亞種的分類地位

    陳詠霞1,梁娜1,李秀1,潘曉睿1,武大勇2

    (1.河北大學生命科學學院,河北保定071002;2.衡水學院生命科學系,河北衡水053000)

    以線粒體DNA細胞色素b(Cytb)基因的全序列作為遺傳標記,探討花斑鰍3亞種(指名亞種Cobitismelanoleuca melanoleuca、北方鰍亞種C.m.granoei和格氏鰍亞種C.m.gladkovi)的分類地位.基于Kimura雙參數(shù)法計算河北花斑鰍指名亞種個體間的遺傳距離為0.2%~1.1%,俄羅斯花斑鰍指名亞種、北方鰍亞種和格氏鰍亞種之間的平均遺傳距離為0.1%~1.5%,而河北花斑鰍指名亞種與俄羅斯花斑鰍3個亞種之間的平均遺傳距離為1.7%~2.6%.分子系統(tǒng)樹顯示,花斑鰍分為兩大分支:分支I為河北花斑鰍指名亞種8個單倍型類群,分支Ⅱ含俄羅斯花斑鰍指名亞種、北方鰍亞種和格氏鰍亞種以及遼寧北方鰍亞種.結合已有研究結果得出:花斑鰍指名亞種僅分布于中國,而俄羅斯也僅有北方鰍亞種1個物種,兩者的分化時間約為2.50~3.09百萬年.

    花斑鰍;線粒體細胞色素b基因;遺傳標記;遺傳距離;分類地位

    鰍屬(Cobitis)隸屬于硬骨綱(Osteichthyes)鯉形目(Cypriniformes)鰍科(Cobitidae),為淡水底層魚類,喜棲息于砂泥底質的清澈流水環(huán)境中,少數(shù)棲息于湖泊內,廣泛分布于古北區(qū)和東洋區(qū)[1-2].據(jù)統(tǒng)計,世界鰍屬魚類約有84種[3-4],是鰍科種類最為繁盛和分化最為明顯的屬之一.

    鰍屬魚類因其高度的表型可塑性、兩性異型和小生境適應性,給該類群的形態(tài)分類與物種鑒定帶來很大困難[5-6],導致其近緣種的分類和命名非?;靵y,同物異名現(xiàn)象時有發(fā)生[7-8].過去,鰍屬的模式種花鰍(C. taenia Linnaeus,1758)被認為是廣泛分布的多型種,歐亞大陸的許多鰍魚被記述為花鰍或其亞種.例如,Nichols[9]描記的新亞種花斑鰍(C.t.melanoleuca Nichols,1925,模式產地:中國山西省清徐縣附近的汾河),Rendahl[10]描記的新亞種北方鰍(C.t.granoei Rendahl,1935,模式產地:俄羅斯西伯利亞鄂木斯克附近的額爾齊斯河)以及Gladkov[11]描記的新亞種西伯利亞鰍(C.t. sibirica Gladkov,1935,模式產地:俄羅斯烏拉爾南部Turaoiak湖泊).

    Berg[12]認為北方鰍和西伯利亞鰍為同種,北方鰍是西伯利亞鰍的次定同物異名.Bǎcescu和Maier[13]將西伯利亞鰍提升為種.但是,Nalbant等[14]認為北方鰍是有效種名,并描記了蒙古鰍新亞種(C.granoeiolivai Nalbant,Holcík etPivnicka,1970,模式產地:蒙古后杭愛省色楞格水系鄂爾渾河支流Ogijn湖泊和Narijn河).Vasil’eva[15]也認為北方鰍是西伯利亞鰍的先定同物異名,但是根據(jù)形態(tài)特征及核型分析,他否定了蒙古鰍亞種的有效性.Nalbant[5]研究了花斑鰍的模式標本(1920年采自中國山西清徐)和俄羅斯西伯利亞庫達河(Kuda)與Ordinskoe湖泊北方鰍標本的形態(tài)特征,指出花斑鰍與北方鰍并沒有明顯差異,并根據(jù)命名優(yōu)先原則,認為花斑鰍是北方鰍的先定同物異名,并認為花斑鰍廣泛分布于黃河以北包括韓國、蒙古、俄羅斯和哈薩克斯坦等地區(qū).Kottelat[16]對采自蒙古克魯倫河(Kherlen)、色楞格河(Selenge)和布爾干河(Bulgan)的魚類進行體色斑紋研究,也認為蒙古鰍亞種、北方鰍和西伯利亞鰍均為花斑鰍的次定同物異名.之后,Prokofiev[17]檢視采自蒙古鄂爾渾河(Orkhon)鰍魚的體色斑紋,將蒙古鰍亞種(C.granoeiolivai)提升為種.Vasil’ev等[18]依據(jù)核型分析,認為蒙古鰍和西伯利亞鰍是北方鰍的次定同物異名,并認為北方鰍是花斑鰍的亞種(C.melanoleuca granoei),同時又描記新亞種格氏花鰍(C.melanoleuca gladkovi Vasil’ev et Vasil’eva,2008,模式產地:俄羅斯Staraya Stanitsa州,Starodubovka村,頓河支流北頓涅茨河).然而,Perdices等[8]依據(jù)Cytb、RAG1、S7和RHO等基因序列對采自俄羅斯的北方鰍亞種(C.m.granoei)、格氏花鰍亞種(C.m.gladkovi)和花斑鰍指名亞種(C.m. melanoleuca)進行研究后認為,北方鰍和格氏花鰍均為花斑鰍的次定同物異名.

    雖然對花斑鰍及其亞種的分類地位進行了很多研究報道,但中國的種群數(shù)據(jù)缺乏,其研究結果仍存在疑問.為了進一步確立花斑鰍及其亞種的分類地位,本文通過分析中國、韓國、俄羅斯等地區(qū)花斑鰍及其亞種的線粒體Cytb基因序列在種間、亞種間的變異,以探討花斑鰍及其亞種的分類地位.

    1 材料與方法

    1.1 實驗材料

    采用電網和蝦籠捕撈鰍屬魚類.于2013年4—10月份采集于河北省拒馬河(易縣紫荊關)、洋河(下花園區(qū)雞鳴驛)、湯河(豐寧鄧珊子村)、灤河(塞罕壩)和白河(赤城縣田鄉(xiāng)).標本浸泡于體積分數(shù)95%的乙醇中,-20℃保存?zhèn)溆?標本保存在河北大學博物館.

    1.2 DNA提取、序列擴增及測定

    采用上海捷瑞生物有限公司生產的細胞/組織基因組DNA提取試劑盒(離心柱型)提取純化花斑鰍基因組DNA.用于擴增線粒體Cytb的引物分別為L14724(5′-GACTTGAAAAACCACCGTTG-3′)和H15915(5′-CTCCGATCTCCGGATTACAAGAC-3′)[19].PCR反應體系總體積為50μL,其中:引物濃度1.0μmol/L各2μL,DNA模版4μL,2×Taq MasterMix Loading Dye 25μL,ddH2O 17μL.PCR反應條件為:95℃預變性3 min,95℃變性30 s,50℃退火30 s,72℃延伸30 s,循環(huán)30次;72°C延伸10min,4℃終止反應.

    PCR產物經質量分數(shù)為1.2%的瓊脂糖凝膠電泳檢測后,用膠回收試劑盒(天根生化科技有限公司)進行PCR產物的回收和純化.純化后的PCR產物由北京諾賽基因組研究中心進行雙向測通.

    1.3 基因序列信息及數(shù)據(jù)分析

    從GenBank下載了花斑鰍指名亞種(C.m. melanoleuca)、格氏花鰍亞種(C.m.gladkovi)、北方鰍亞種(C.m.granoei)、西伯利亞鰍亞種(C.m.sibirica)(=北方鰍亞種)、喬氏鰍(C.choii)、太平鰍(C.pacifica)、四線鰍(C.tetralineata)和韓鰍(C.hankugensis)等鰍屬魚類28個個體的Cytb基因同源序列,與本研究所檢測的12個個體(白河(HBBH)2個個體、湯河(HBTH)2個個體、拒馬河(HBJMH)3個個體、洋河(HBYH)3個個體和灤河(HBLH)2個個體)的Cytb基因序列,共計40個一并分析.另選取與鰍屬親緣關系較近的薩瓦納鰍屬(Sabanejewia)的南斯拉夫薩瓦納鰍(Sabanejewia balcanica)作為外群進行系統(tǒng)發(fā)育分析,相關信息見表1.

    使用BioEdit7.2.5軟件進行同源序列拼接、排列及手動矯正.用DnaSP 5.10軟件[20]確定所測序列的單倍型.利用MEGA 5.0軟件[21]對測定序列和GenBank下載的序列進行比對,排序后去除兩端冗余序列,統(tǒng)計序列的平均堿基組成、變異位點、簡約信息位點數(shù)、轉換/顛換比率(Ts/Tv)和遺傳距離,其中遺傳距離的計算采用Kimura雙參數(shù)進化模型.

    表1 本研究中所分析的鰍屬魚類和來源Tab.1 Speciesand sourcesof Cobitis studied

    采用最大似然法(ML)和貝葉斯法(BI)構建分子系統(tǒng)樹,并用jModelTest2.1.4軟件[22]選擇最適合的進化模型.運用MEGA 5.0軟件構建ML法的分子系統(tǒng)樹,替代模型為GTR,進行1 000次自展分支檢驗.應用MrBayes 3.0軟件[23]構建BI法分子系統(tǒng)樹,隨機選取起始樹,以后驗概率(posterior probability,PP)來表示各分支的可信度.替代模型根據(jù)hLRTs檢驗結果設定為“nst=6 rates=gamma”,計算出最適合的替代模型為GTR+G+I.馬爾科夫鏈的蒙特卡洛方法(Markov chain MonteCarlo process)設置為4條鏈同時運行2×106代,3條熱鏈1條冷鏈.計算2×106代,每100代對系統(tǒng)樹進行抽樣,最終得到20 001棵系統(tǒng)發(fā)育樹,重復一次以保證所獲得結果的準確度.將運行過程中所得的對數(shù)似然值(Log-likelihood scores)與相應代數(shù)進行作圖,找到對數(shù)似然值達到飽和的位置,飽和前的數(shù)據(jù)被作為老化樣本而舍棄.在舍棄老化樣本(burn-in=4 000)后,根據(jù)剩余的樣本構建一致樹(consensus tree)并計算相關參數(shù).

    2 結果與分析

    2.1 堿基組成及序列變異

    所測定的河北省5條河流12個個體的線粒體Cytb基因序列(1 140 bp)中有8個單倍型序列(表1),與GenBank下載的28個同源序列合并分析.結果表明,鰍屬8種(亞種)36個Cytb單倍型序列的平均堿基組成為T 32.1%、C 24.8%、A 28.0%、G 15.1%,A+T含量(60.1%)高于C+G含量(39.9%),堿基組成表現(xiàn)出明顯的T偏倚.其中,第1密碼子位點C+G含量(50.3%)顯著高于第2和第3密碼子位點(分別為38.7%和30.7%).

    所分析序列共含有保守位點816個,約占序列總長(下同)的71.6%;變異位點數(shù)324個,約占28.4%;簡約信息位點255個,約占22.4%.全部位點中不變位點有1 045個,轉換位點有70個,顛換位點有17個,序列中的轉換明顯多于顛換,轉換顛換比值(Ts/Tv)平均為4.13,序列突變未達到飽和.其中不變位點數(shù)第2密碼子位點最多,有377個;而轉換位點和顛換位點均是第3密碼子位點最多,分別有58個和16個.

    2.2 遺傳距離

    表2列出了基于Kimura雙參數(shù)進化模型計算的鰍屬種(亞種)間的平均遺傳距離.結果顯示,鰍屬6種間平均遺傳距離最大的是西伯利亞鰍(=北方鰍亞種)和黑龍江鰍(19.1%),四線鰍與韓鰍的最小(8.8%);種間遺傳距離明顯大于花斑鰍亞種間的平均遺傳距離(0.5%~2.6%).在花斑鰍3個亞種中,亞種間的平均遺傳距離最大值為河北花斑鰍指名亞種與西伯利亞鰍(=北方鰍亞種)(2.6%),其次是河北花斑鰍指名亞種與北方鰍亞種、俄羅斯花斑鰍指名亞種(均為1.9%),最小值為俄羅斯花斑鰍指名亞種與北方鰍亞種、格氏花鰍亞種(均為0.5%),該最小值接近于花斑鰍亞種內的平均遺傳距離(0.1%~0.6%).河北花斑鰍指名亞種與俄羅斯花斑鰍3個亞種間的平均遺傳距離(1.7%~2.6%)明顯大于俄羅斯花斑鰍3個亞種之間的平均遺傳距離(0.5%~1.4%).

    表2 基于Kimura雙參數(shù)進化模型的鰍屬8種(亞種)間Cytb基因單倍型的平均遺傳距離Tab.2 Averagegenetic distancesamong eight species(subspecies)of Cobitis for Cytb genehap lotypes according to Kimura-2-parametermodel%

    2.3 分子系統(tǒng)進化樹

    應用貝葉斯法(BI法)和最大似然法(ML法)分別構建了鰍屬魚類的分子系統(tǒng)進化樹(圖1和圖2),各分支的數(shù)字表示ML樹中該支的bootstrap支持率和BI樹中該支的后驗概率.兩系統(tǒng)樹的拓撲結構基本一致,鰍屬8種(亞種)36個個體分為2大類:第1大類(類群Ⅰ)由花斑鰍3個亞種和喬氏鰍以高置信度(BI 100%,ML 96%)聚在一起;黑龍江鰍、太平鰍、四線鰍和韓鰍4個種類聚為第2類群(類群Ⅱ)(BI85%,ML 82%).

    在類群Ⅰ中,分布于俄羅斯的北方鰍亞種、西伯利亞鰍亞種(=北方鰍亞種)、格氏花鰍亞種、花斑鰍指名亞種以及分布于我國遼寧的北方鰍亞種以低置信度(BI58%,ML 91%)聚為一支;分布于河北的花斑鰍指名亞種同樣以較低置信度(BI64%,ML 80%)聚為另一支,然后這2支以高置信度(BI100%,ML 100%)聚在一起構成一大分支.與此分支構成姐妹群關系的是由韓國喬氏鰍和分布于我國黑龍江的北方鰍亞種以高置信度(BI 100%,ML 100%)聚在一起構成的分支.

    圖1 基于8種(亞種)鰍屬魚類Cytb基因單倍型構建的貝葉斯法(BI)分子系統(tǒng)樹Fig.1 M olecular phylogenetic tree of 8(subspecies)speciesof Cobitis based on Cytb gene hap lotepes using Bayesian(BI)inferencemethod

    圖2 基于8種(亞種)鰍屬魚類Cytb基因單倍型構建的最大似然法(M L)分子系統(tǒng)樹Fig.2 M olecular phylogenetic tree of8(subspecies)species of Cobitis based on Cytb gene hap lotypesusing maximum likelihood inference(M L)method

    3 討論

    遺傳距離分析顯示,分布于河北的花斑鰍與分布于俄羅斯的花斑鰍3個亞種(花斑鰍指名亞種、北方鰍亞種和格氏鰍亞種)間的平均遺傳距離為1.7%~2.6%,大于河北花斑鰍個體間的遺傳距離(0.2%~1.1%),也大于俄羅斯花斑鰍3個亞種間的遺傳距離(0.1%~1.5%).遺傳距離雖然不是物種的鑒別特征,但可以為物種的分類地位提供佐證[8].Bohlen等[24]基于Cytb基因序列對小亞細亞和巴爾干半島鰍屬魚類的研究中,斑條鰍(C.punctilineata)和瘤鰍(C.strumicae)的遺傳距離為1.9%.Perdices等[8]基于Cytb基因序列分析俄羅斯境內的花斑鰍指名亞種、北方鰍亞種和格氏鰍亞種的遺傳距離(0.0%~1.3%),認為上述3個亞種的遺傳距離與地中海淡水魚類種內遺傳距離(<2%)[25]相近,并以此為依據(jù),認為俄羅斯地區(qū)僅有1個種,無亞種分化.由此可見,從遺傳距離來看,河北花斑鰍與俄羅斯花斑鰍之間的遺傳分化應屬于亞種級分化水平.

    在系統(tǒng)進化上,分布于中國河北、遼寧和俄羅斯的花斑鰍以高置信度(BI 100%,ML 100%)聚為一支單系,其中分布于中國河北的花斑鰍聚為一支,而分布于遼寧和俄羅斯的花斑鰍聚為另一支,雖然兩分支的置信度(分別為BI 64%和ML 80%;BI 58%和ML 91%)相對較低,但提示兩者已有分化,系統(tǒng)發(fā)育關系的結果與遺傳距離顯示的結果一致.從地理起源和演化上看,鰍屬魚類從中新世,約15百萬年前[26]由中國南方向北擴散,經中國東部、日本、朝鮮、蒙古到達歐洲和北非[27-28].花斑鰍分布廣泛,從中國黃河以北到韓國、蒙古、俄羅斯至哈薩克斯坦均有分布[4].遼寧北方鰍亞種在地理分布上,介于河北花斑鰍和俄羅斯花斑鰍之間.系統(tǒng)進化上,遼寧北方鰍亞種與俄羅斯花斑鰍3個亞種具有較近的親緣關系,遺傳分化程度也支持二者聚為一支;遼寧北方鰍亞種與俄羅斯花斑鰍3個亞種的遺傳分化(遺傳距離為0.9%~1.4%)明顯低于其與河北花斑鰍的遺傳分化(遺傳距離為1.7%~2.2%).由于本研究中遼寧北方鰍亞種Cytb序列僅1條,不排除因序列少導致該亞種在系統(tǒng)樹分析中位置有差異.但以上分析表明,河北花斑鰍與俄羅斯花斑鰍為亞種級分化水平.

    俄羅斯花斑鰍亞種名為北方鰍(C.m.granoei)還是西伯利亞鰍(C.m.sibirica),國內外學者存在明顯分歧.Berg[12]和Bǎcescu等[13]認為北方鰍是西伯利亞鰍的次定同物異名.但是,Nalbant等[14]、Vasil’eva[15]、Vasil’ev等[18]認為北方鰍是西伯利亞鰍的先定同物異名.Bogutskaya等[29]研究這2個種名的先定性時,認為北方鰍(C.t.granoei)刊載時間(于1935年10月16日刊載于《NaturaeNovitiates》)早于西伯利亞鰍(C.t.sibirica)(于1935年11月19日刊載于《Naturae Novitiates》).因此,根據(jù)命名優(yōu)先原則,俄羅斯花斑鰍亞種名為北方鰍亞種(C.m.granoei),分布于遼河以北至韓國、蒙古、俄羅斯和哈薩克斯坦;花斑鰍指名亞種(C.m. melanoleuca)分布于黃河以北至河北境內.按鰍屬魚類線粒體Cytb的進化速率為每百萬年變化0.68~0.84的模型計算[8],花斑鰍指名亞種和北方鰍亞種的分歧年代距今為2.50~3.09百萬年.

    本研究中,分布于中國黑龍江的北方鰍亞種與韓國喬氏鰍以高置信度(BI 100%,ML 100%)聚為一支單系,系統(tǒng)發(fā)育結果與遺傳分化程度一致;黑龍江北方鰍亞種與其他花斑鰍亞種的遺傳分化(平均遺傳距離為14.1%~15.2%)已達到種的分化水平,而與韓國喬氏鰍的遺傳分化(遺傳距離為0.3%~1.1%)屬于種內分化范圍.因此,根據(jù)該研究的成果,黑龍江的北方鰍亞種應與喬氏鰍是同種,為我國的新紀錄種,但需結合形態(tài)特征作進一步研究.

    [1]SAWADAY.Phylogenyand zoogeographyof thesuperfamily Cobitoidea(Cyprinoidei,Cypriniformes)[J].Memoirs of the Faculty of Fisheries HokkaidoUniversity,1982,28(2):65-223.

    [2]BAˇNAˇRESCU P.General Distribution and Dispersal of Freshwater Animals:Zoogeography of FreshWaters[M].Wiesbaden:Aula-Verlag. 1990.

    [3]FROESE R,PAULY D.FishBase[M/OL].[2014-05-26].http://www. fishbase.org.

    [4]CHEN Y X,SUIX Y,CHEN Y F,etal.Three new speciesof cobitid fish genus Cobitis(Teleostei,Cobitidae)from the River Pearl basin of China[J].Folia Zoologica,2015,64(1),1-16.

    [5]NALBANT T T.Some problems in the systematics of the genus Cobitis and its relatives(Pisces,Ostariophysi,Cobitidae)[J].RevueRoumaine deBiologie:SerieBiologie Animale,1993,38(2):101-110.

    [6]JANKOK,SˇFLAJHANSM,CHOLEVA L,etal.Diversity ofEuropean spined loaches(genus Cobitis L.):an update of the geographic distribution of the Cobitis taenia hybrid complex with a description of newmolecular tools forspeciesdetermination[J].JournalofFish Biology,2007,71(Supplement):387-408.

    [7]BORON A,OZOUF-COSTAZ C,COUTANCEAU JP,et al.Gene mapping of 28Sand 5S rDNA sites in the spined loach Cobitis taenia(Pisces,Cobitidae)from a diploid population and a diploic-tetraploid population[J].Genetica,2006,128(1/2/3):71-79.

    [8]PERDICESA,VASIL’EVA E,VASIL’EV V.From Asia to Europe across Siberia:phylogeographyof the Siberian spined loach(Teleostei,Cobitidae)[J].Zoologica Scripta,2015,44(1):29-40.

    [9]NICHOLSJT.The Two Chinese Loachesof the Genus Cobitis[M].New York:The AmericanMuseum ofNaturalHistory,1925.

    [10]RENDAHLH.Ein paarneueUnteraren von Cobitis taenia[J].Menoranda Societatispro FaunaetFlora Fennica,1935(10):330-332.

    [11]GLADKOV NA.Sur la variabilitégéographique de Cobitis taenia L.[J]. Sbornik Trudov Zoologicheskogo Muzeya Moskovskogo Universiteta,1935(2):69-74.

    [12]BERGLS.FreshwaterFishesoftheU.S.S.R.and AdjacentCountries[M]. Moscow-Leningrad:Izdatel’stvo Acad Nauk USSR,1949.

    [13]BAˇCESCUM,MAIERR.The Cobitis of the Don and Volga[J].Journal of Ichthyology,1969(9):38-44.

    [14]NALBANT T T,HOLCˇIK J,PIVNICˇKA K.A new loach,Cobitis granoei olivai,ssp.n.,from Mongolia,with some remarks on the Cobitis elongata-belseli-macrostigma group(Pisces,Ostariophysi,Cobitidae)[J].Véstnik CˇeskoslovenskéSpole Cˇnosti Zoologické,1970(34):121-128.

    [15]VASIL’EVAED.Redescription,morpho-ecologicalcharacteristicsand distribution of Cobitisgranoei(Teleostei,Cobitidae)[J].Zoologicheskii Zhurnal,1988(67):1025-1036.

    [16]KOTTELATM.FishesofMongolia.ACheck-Listof the FishesKnown to Accur inMongoliawith Commentson Systematicsand Nomenclature[M]. Washington:TheWorld Bank,2006.

    [17]PROKOFIEV AM.Cobitisolivai Nalbant,Hol ik&Pivni ka 1970,a valid species of spined loaches from Mongolia:redescription and comparison with C.melanoleuca Nichols1925(Osteichthyes,Cypriniformes,Cobitidae)[J].SenckenbergianaBiologica,2007(87):111-124. [18]VASIL’EV V P,VASIL’EVA ED.Comparativekaryologyofspeciesof thegenera Misgurnus and Cobitis(Cobitidae)from theAmur Riverbasin in connection with their taxonomic relations and the evolution of karyotypes[J].Journalof Ichthyology,2008,48(1):1-13.

    [19]XIAOW H,ZHANG Y P,LIUH Z.Molecular systematicsof Xenocyprinae(Teleoatei:Cyprinidae):taxonomy,biogengraphy,and coevolution of a special group restricted in East Asia[J].Molecular Phylogenetics and Evolution,2001,18(2):163-173.

    [20]LIBEADO P,ROZAS J.Dnasp v5:a software for comprehensive analysisofDNA polymorphism data[J].Bioinformatics,2009,25(11):1451-1452.

    [21]TAMURAK,PETERSOND,PETERSONN,etal.MEGA 5:molecular evolutionary genetics analysis usingmaximum likelihood,evolutionary distance,andmaximum parsimonymethods[J].Molecular Biology and Evolution,2011,28(10):2731-2739.

    [22]GUINDON S,GASCUELO.A simple,fast,and accurate algorithm to estimate largephylogeniesbymaximum likelihood[J].Systematic Biology,2003,52(5):696-704.

    [23]HUELSENBECK JP,RONQUIST F.MrBayes:bayesian inference of phylogeny trees[J].Bioinformatics,2001,17(8):754-755.

    [24]BOHLENJ,PERDICESA,DOADRIOI,etal.Vicariance,colonisation,and fast localspeciation in AsiaMinorand the Balkansas revealed from thephylogenyofspined loaches(Osteichthyes;Cobitidae)[J].Molecular Phylogeneticsand Evolution,2006,39(2):552-61.

    [25]GEIGERM F,HERDER F,MONAGHANM T,etal.Spatial heterogeneity in the Mediterranean Biodiversity hotspot affects barcoding accuracy of its freshwater fishes[J].Molecular Ecology Resources,2014,14(6):1210-1221.

    [26]周家健.山東山旺中中新世鰍科化石[J].古脊椎動物學報.1992,30(1):71-76.

    [27]陳景星,朱松泉.鰍科魚類亞科的劃分及其宗系發(fā)生的相互關系[J].動物分類學報,1984,9(2):201-207.

    [28]LECHTOV V,BOHLEN J,PERDICESA.Molecular phylogeny of the freshwater fish family Cobitidae(Cyriniformes:Teleostei):delimitation of genera,mitochondrial infrogression and evolution of sexual dimorphism[J].Molecular Phylogenetics and Evolution,2008,47(2):812-831.

    [29]BOGUTSKAYANG,NASEKAAM.CatalogueofAgnathansand Fishes ofFreshand BrackishWatersofRussiawith CommentsonNomenclature and Taxonomy[M].Moscow:KMK Scientific Press,2004.

    (責任編校 紀翠榮)

    Taxonom ic statusof three subspeciesof Cobitismelanoleuca based on sequencesofm itochondrial cytochrome b gene

    CHEN Yongxia1,LIANGNa1,LIXiu1,PANXiaorui1,WUDayong2
    (1. College of Life Sciences,Hebei University,Baoding 071002,Hebei Province,China;2. College of Life Sciences,Hengshui University,Hengshui 053000,Hebei Province,China)

    The taxonomic status of three subspecies(Cobitismelanoleucamelanoleuca,C.m.granoei,C.m.gladkovi)of Cobitismelanoleuca was discussed based on the whole sequences ofmitochondrial cytochrome b gene taken as hereditary marker.The genetic distance was calculated by themethod of Kimura-2-parameter.The results showed that the individual genetic distance of C.m.melanoleuca in Hebei Province was 0.2%-1.1%;the average genetic distance among the three subspeciesof Cobitismelanoleuca in Russiawas0.1%-1.5%.And the average genetic distance between C.m.melanoleuc and the three subspeciesof Cobitismelanoleuca in Russiawas1.7%-2.6%.Themolecular phylogenic trees demonstrated that the speciesof Cobitismelanoleuca wasdivided into two distinctclades.One clade included eighthyplogroupsof C.m.melanoleuca in HebeiProvince;the other clade included the three subspecies in Russia and C.m.granoei in Liaoning Province in China. Based on the existing researches and the presentstudy,the conclusion was that C.m.melanoleuca only distributed in China and therewasonly one subspecies,i.e.C.m.granoei in Russia.The divergence time between these two subspecieswasabout 2.50-3.09Mya.

    Cobitismelanoleuca;mitochondrial cytochrome b gene;hereditarymarker;genetic distance;taxonomic status

    Q959.4

    A

    2015-04-25

    國家自然科學基金資助項目(31372166).

    陳詠霞(1976—),女,副教授,主要從事魚類分類、系統(tǒng)發(fā)育與生物地理學方面的研究.

    1671-1114(2015)03-0001-06

    猜你喜歡
    花斑亞種西伯利亞
    A comprehensive review on the chemical constituents,sesquiterpenoid biosynthesis and biological activities of Sarcandra glabra
    亞沉茶漬亞洲亞種Lecanora subimmersa subsp. asiatica Zahlbr.的訂正研究
    盤羊新亞種
    ——和田盤羊
    花斑牛的學習秘訣
    西伯利亞的熊媽媽
    小讀者(2019年20期)2020-01-04 02:13:32
    花斑牛的學習秘訣
    西伯利亞博亞雷巖畫與歐亞大草原的青銅時代
    花斑熊
    大灰狼(2015年5期)2015-07-25 20:57:36
    大西洋鮭殺鮭氣單胞菌無色亞種的分離鑒定和致病性研究
    花斑豬的幸運日
    天堂中文最新版在线下载| 欧美老熟妇乱子伦牲交| 亚洲九九香蕉| 成人黄色视频免费在线看| 欧美成人免费av一区二区三区| 黄片播放在线免费| 亚洲三区欧美一区| 男女做爰动态图高潮gif福利片 | 99精品欧美一区二区三区四区| 国产亚洲精品久久久久5区| 在线观看免费日韩欧美大片| 在线播放国产精品三级| 97人妻天天添夜夜摸| 涩涩av久久男人的天堂| 欧美日韩乱码在线| aaaaa片日本免费| 看免费av毛片| 国产精品综合久久久久久久免费 | 黄色视频不卡| av福利片在线| 美女福利国产在线| 97碰自拍视频| 精品久久久久久久毛片微露脸| 成人黄色视频免费在线看| 最新美女视频免费是黄的| 亚洲午夜精品一区,二区,三区| 国产男靠女视频免费网站| 黄色视频不卡| 性少妇av在线| 国产伦一二天堂av在线观看| 久久久久亚洲av毛片大全| 长腿黑丝高跟| 性色av乱码一区二区三区2| 91国产中文字幕| 色尼玛亚洲综合影院| 99国产精品一区二区三区| 国产精品影院久久| 香蕉丝袜av| 久久这里只有精品19| 久久青草综合色| 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| a级毛片在线看网站| 老司机福利观看| av电影中文网址| 自线自在国产av| 亚洲av成人av| 国产麻豆69| 国产精品野战在线观看 | 夜夜夜夜夜久久久久| 在线国产一区二区在线| 伊人久久大香线蕉亚洲五| 久久精品影院6| 欧美成狂野欧美在线观看| 久久国产精品人妻蜜桃| 亚洲欧美激情在线| 亚洲久久久国产精品| 超色免费av| 久久午夜综合久久蜜桃| 身体一侧抽搐| 视频区欧美日本亚洲| 国产极品粉嫩免费观看在线| 久久精品国产综合久久久| 国产深夜福利视频在线观看| 亚洲av第一区精品v没综合| 97超级碰碰碰精品色视频在线观看| 国产成人免费无遮挡视频| 免费在线观看影片大全网站| 国产欧美日韩精品亚洲av| 久久性视频一级片| 一级毛片精品| 免费看a级黄色片| 久久久精品欧美日韩精品| 无人区码免费观看不卡| 成人永久免费在线观看视频| 精品国产一区二区三区四区第35| 日韩大码丰满熟妇| 黄片播放在线免费| 亚洲av成人不卡在线观看播放网| 香蕉丝袜av| 久久精品国产99精品国产亚洲性色 | 日韩av在线大香蕉| 精品第一国产精品| 日本免费一区二区三区高清不卡 | 精品熟女少妇八av免费久了| 色综合婷婷激情| 免费看a级黄色片| 日韩欧美免费精品| 性欧美人与动物交配| 好看av亚洲va欧美ⅴa在| 久久人人精品亚洲av| 国产av精品麻豆| 欧美黄色淫秽网站| 欧美日韩av久久| 丰满的人妻完整版| 久久午夜亚洲精品久久| 色播在线永久视频| 啦啦啦免费观看视频1| 久久精品国产亚洲av高清一级| 中出人妻视频一区二区| 在线观看免费视频日本深夜| 亚洲精品美女久久av网站| 午夜精品久久久久久毛片777| 日韩大尺度精品在线看网址 | 99在线人妻在线中文字幕| 免费在线观看黄色视频的| 精品高清国产在线一区| 亚洲第一av免费看| 在线观看午夜福利视频| 久久热在线av| 国产男靠女视频免费网站| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| 在线观看www视频免费| 日本精品一区二区三区蜜桃| 国产亚洲欧美98| 中文字幕精品免费在线观看视频| 精品国产美女av久久久久小说| 成人手机av| 一级毛片女人18水好多| 99久久人妻综合| 黄色视频,在线免费观看| 性欧美人与动物交配| 色老头精品视频在线观看| x7x7x7水蜜桃| 脱女人内裤的视频| 免费看a级黄色片| 久久99一区二区三区| 大码成人一级视频| 亚洲欧美日韩高清在线视频| 午夜日韩欧美国产| 桃红色精品国产亚洲av| 91麻豆av在线| 欧洲精品卡2卡3卡4卡5卡区| 免费观看精品视频网站| 国产精品电影一区二区三区| 一边摸一边抽搐一进一出视频| 自线自在国产av| 免费少妇av软件| 五月开心婷婷网| 天堂中文最新版在线下载| 国产精品美女特级片免费视频播放器 | 午夜福利影视在线免费观看| 色婷婷久久久亚洲欧美| 日韩免费av在线播放| 三上悠亚av全集在线观看| 欧美黑人欧美精品刺激| 99在线人妻在线中文字幕| 一进一出抽搐动态| 亚洲人成电影免费在线| 国产99久久九九免费精品| 超碰97精品在线观看| 国产成年人精品一区二区 | 亚洲国产欧美网| 国产精品综合久久久久久久免费 | 久久婷婷成人综合色麻豆| 精品国产亚洲在线| 国产精品一区二区在线不卡| 亚洲av片天天在线观看| 国产精品成人在线| 18禁美女被吸乳视频| 国产野战对白在线观看| 最好的美女福利视频网| 一级毛片精品| 久久伊人香网站| 精品久久久久久成人av| √禁漫天堂资源中文www| 高清av免费在线| 俄罗斯特黄特色一大片| 国产免费男女视频| 99久久精品国产亚洲精品| 国产精品 国内视频| 99在线人妻在线中文字幕| 一a级毛片在线观看| 国产亚洲精品久久久久久毛片| 搡老熟女国产l中国老女人| 国产在线精品亚洲第一网站| www.精华液| 久久这里只有精品19| 成年人黄色毛片网站| 久久人妻熟女aⅴ| 宅男免费午夜| av福利片在线| 咕卡用的链子| 精品国产美女av久久久久小说| 大香蕉久久成人网| 亚洲中文日韩欧美视频| 日韩成人在线观看一区二区三区| 18禁裸乳无遮挡免费网站照片 | 亚洲成国产人片在线观看| 色婷婷久久久亚洲欧美| 国产午夜精品久久久久久| 三级毛片av免费| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三卡| 欧美日韩福利视频一区二区| 岛国在线观看网站| 99精品欧美一区二区三区四区| 国产欧美日韩综合在线一区二区| 亚洲自偷自拍图片 自拍| 久久精品aⅴ一区二区三区四区| 91国产中文字幕| 99国产精品99久久久久| 亚洲精品久久成人aⅴ小说| 久久性视频一级片| 亚洲一区二区三区不卡视频| 亚洲精品一二三| av福利片在线| 亚洲欧美精品综合一区二区三区| 亚洲av电影在线进入| 婷婷丁香在线五月| 女人被躁到高潮嗷嗷叫费观| 国产成人精品久久二区二区91| 国产高清videossex| 成人av一区二区三区在线看| 欧美日韩乱码在线| 精品久久久久久久久久免费视频 | 日韩高清综合在线| 久久精品国产综合久久久| 在线观看一区二区三区| 国产欧美日韩综合在线一区二区| 午夜激情av网站| 久久亚洲精品不卡| 精品一区二区三区av网在线观看| 国产高清视频在线播放一区| 国产高清视频在线播放一区| 桃色一区二区三区在线观看| 国产一区二区在线av高清观看| 91大片在线观看| 中文字幕人妻丝袜一区二区| 精品一区二区三区视频在线观看免费 | 久久精品国产清高在天天线| 日韩精品青青久久久久久| 午夜激情av网站| 精品第一国产精品| 叶爱在线成人免费视频播放| 一级作爱视频免费观看| 在线十欧美十亚洲十日本专区| 他把我摸到了高潮在线观看| 久久久精品国产亚洲av高清涩受| 精品福利观看| 女性生殖器流出的白浆| 欧洲精品卡2卡3卡4卡5卡区| 日本免费一区二区三区高清不卡 | 久久这里只有精品19| 手机成人av网站| 天堂动漫精品| 91麻豆av在线| 一进一出抽搐gif免费好疼 | 国产精品一区二区免费欧美| 国产免费现黄频在线看| www国产在线视频色| 午夜精品久久久久久毛片777| 一夜夜www| 亚洲久久久国产精品| 亚洲av美国av| www.999成人在线观看| 精品国产美女av久久久久小说| 久久人人精品亚洲av| ponron亚洲| 999久久久精品免费观看国产| 又黄又爽又免费观看的视频| 午夜91福利影院| 久久人人爽av亚洲精品天堂| 午夜成年电影在线免费观看| 成年女人毛片免费观看观看9| 后天国语完整版免费观看| 在线永久观看黄色视频| 三上悠亚av全集在线观看| 亚洲专区字幕在线| 成人亚洲精品av一区二区 | 国产亚洲精品综合一区在线观看 | 一二三四社区在线视频社区8| 国产91精品成人一区二区三区| 日韩精品青青久久久久久| 97人妻天天添夜夜摸| 色综合欧美亚洲国产小说| 男女下面插进去视频免费观看| 亚洲精品在线美女| 高清av免费在线| 1024视频免费在线观看| 亚洲视频免费观看视频| 91字幕亚洲| 又紧又爽又黄一区二区| 超碰成人久久| 99久久99久久久精品蜜桃| 丰满迷人的少妇在线观看| www.精华液| 久久人人爽av亚洲精品天堂| 精品午夜福利视频在线观看一区| 成年人免费黄色播放视频| 中文字幕人妻熟女乱码| 99国产精品免费福利视频| 婷婷六月久久综合丁香| 国产三级黄色录像| 男女下面插进去视频免费观看| 成人亚洲精品一区在线观看| 日韩精品免费视频一区二区三区| 国产成年人精品一区二区 | 亚洲精品中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| av网站在线播放免费| 别揉我奶头~嗯~啊~动态视频| 久久精品91无色码中文字幕| 国内毛片毛片毛片毛片毛片| 极品教师在线免费播放| 国产亚洲欧美在线一区二区| 91国产中文字幕| 国产亚洲精品综合一区在线观看 | 国产有黄有色有爽视频| 久久午夜综合久久蜜桃| 黑人巨大精品欧美一区二区蜜桃| 久久99一区二区三区| 欧美另类亚洲清纯唯美| 国产精品久久视频播放| 少妇的丰满在线观看| 亚洲国产精品合色在线| 日韩欧美一区视频在线观看| 久久精品国产亚洲av高清一级| 色婷婷av一区二区三区视频| 久久性视频一级片| 正在播放国产对白刺激| 悠悠久久av| 亚洲九九香蕉| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 国产黄a三级三级三级人| 亚洲免费av在线视频| 老熟妇乱子伦视频在线观看| 欧美黑人欧美精品刺激| 免费不卡黄色视频| 亚洲av五月六月丁香网| 一级a爱片免费观看的视频| 一夜夜www| 丝袜美足系列| 国产一区二区激情短视频| xxx96com| 国产精品一区二区三区四区久久 | av在线天堂中文字幕 | 午夜福利一区二区在线看| 99在线视频只有这里精品首页| 少妇的丰满在线观看| a级毛片在线看网站| 女人被狂操c到高潮| 97人妻天天添夜夜摸| 久久精品国产清高在天天线| 99香蕉大伊视频| 99久久99久久久精品蜜桃| 久久天堂一区二区三区四区| 亚洲一区高清亚洲精品| 亚洲欧美激情在线| 亚洲性夜色夜夜综合| 首页视频小说图片口味搜索| 日韩国内少妇激情av| 国产又爽黄色视频| 日韩视频一区二区在线观看| 人妻久久中文字幕网| 亚洲色图av天堂| 他把我摸到了高潮在线观看| 久久久久国内视频| 成人精品一区二区免费| 国产一区在线观看成人免费| 麻豆国产av国片精品| 欧美人与性动交α欧美软件| 99热只有精品国产| 97超级碰碰碰精品色视频在线观看| 欧美日韩av久久| 成年人免费黄色播放视频| 成人三级做爰电影| 不卡一级毛片| 在线看a的网站| 免费在线观看日本一区| 美女大奶头视频| 午夜精品在线福利| 热99国产精品久久久久久7| 亚洲男人天堂网一区| 色综合欧美亚洲国产小说| 国产无遮挡羞羞视频在线观看| 99国产精品一区二区蜜桃av| 97人妻天天添夜夜摸| a级毛片在线看网站| 一二三四在线观看免费中文在| 久久亚洲真实| 亚洲精品成人av观看孕妇| 18美女黄网站色大片免费观看| 国产亚洲av高清不卡| 超色免费av| 一边摸一边做爽爽视频免费| 久久性视频一级片| 日本免费一区二区三区高清不卡 | www.精华液| 久热这里只有精品99| 亚洲avbb在线观看| 91麻豆精品激情在线观看国产 | 久久精品人人爽人人爽视色| 国产亚洲精品综合一区在线观看 | 亚洲五月色婷婷综合| 久久精品91无色码中文字幕| 999久久久精品免费观看国产| 欧美色视频一区免费| 午夜精品国产一区二区电影| 国产麻豆69| 91成人精品电影| 国产精品亚洲av一区麻豆| 亚洲欧美日韩无卡精品| 久久久国产一区二区| 精品一品国产午夜福利视频| 正在播放国产对白刺激| 一二三四社区在线视频社区8| 国产精品野战在线观看 | a级毛片黄视频| 午夜福利在线观看吧| 亚洲五月色婷婷综合| 久久九九热精品免费| 亚洲精品中文字幕在线视频| 国产麻豆69| 天天躁狠狠躁夜夜躁狠狠躁| svipshipincom国产片| 亚洲第一青青草原| 国产99久久九九免费精品| 国产精品九九99| 一级毛片高清免费大全| 丝袜美腿诱惑在线| 国产精品国产高清国产av| 一级片免费观看大全| 亚洲专区国产一区二区| 一区在线观看完整版| 黑人巨大精品欧美一区二区蜜桃| 欧美一级毛片孕妇| 国产精品久久久人人做人人爽| 欧美久久黑人一区二区| 午夜福利免费观看在线| 久久久久国产精品人妻aⅴ院| av网站免费在线观看视频| 在线观看免费视频日本深夜| 欧美日韩福利视频一区二区| 99久久人妻综合| 久久久久国产一级毛片高清牌| 男人舔女人下体高潮全视频| 久久午夜综合久久蜜桃| 久久久久久久午夜电影 | 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播放欧美日韩| 久9热在线精品视频| 国产高清视频在线播放一区| 亚洲成国产人片在线观看| 国产高清激情床上av| 伦理电影免费视频| 长腿黑丝高跟| 韩国av一区二区三区四区| 大型av网站在线播放| 亚洲专区国产一区二区| 欧美日韩一级在线毛片| 国产野战对白在线观看| 高清在线国产一区| 国产三级黄色录像| www日本在线高清视频| 一级片免费观看大全| 欧美日韩黄片免| 国产欧美日韩综合在线一区二区| 天天影视国产精品| www.熟女人妻精品国产| 国产欧美日韩精品亚洲av| 日本免费a在线| 又黄又粗又硬又大视频| 最近最新免费中文字幕在线| 国产无遮挡羞羞视频在线观看| 精品一区二区三区四区五区乱码| 巨乳人妻的诱惑在线观看| 久久伊人香网站| 长腿黑丝高跟| 老司机深夜福利视频在线观看| 啪啪无遮挡十八禁网站| 老司机靠b影院| x7x7x7水蜜桃| 婷婷六月久久综合丁香| 久久精品亚洲熟妇少妇任你| 亚洲五月天丁香| 国产99白浆流出| 极品人妻少妇av视频| 国产又爽黄色视频| 精品久久久久久久久久免费视频 | 91老司机精品| 亚洲av电影在线进入| 老熟妇乱子伦视频在线观看| 久久久国产成人免费| 久久久国产欧美日韩av| 制服人妻中文乱码| 久久精品人人爽人人爽视色| 欧美日韩精品网址| 男人操女人黄网站| 91成人精品电影| 久久久久精品国产欧美久久久| 黄色视频不卡| 超碰成人久久| 成年人黄色毛片网站| 最好的美女福利视频网| 级片在线观看| 国产成人精品在线电影| 一二三四在线观看免费中文在| 精品熟女少妇八av免费久了| 一a级毛片在线观看| 国产精品一区二区免费欧美| 国内毛片毛片毛片毛片毛片| 亚洲熟女毛片儿| 亚洲欧美激情综合另类| 国产精品国产av在线观看| 美女福利国产在线| 老鸭窝网址在线观看| 激情视频va一区二区三区| 国产精品秋霞免费鲁丝片| 一a级毛片在线观看| 亚洲精品久久午夜乱码| 日本a在线网址| 亚洲三区欧美一区| 精品久久久久久电影网| 欧美黑人精品巨大| 99国产精品99久久久久| av超薄肉色丝袜交足视频| 12—13女人毛片做爰片一| 首页视频小说图片口味搜索| 人人妻人人爽人人添夜夜欢视频| 搡老熟女国产l中国老女人| 日韩人妻精品一区2区三区| 黄色 视频免费看| 亚洲一区高清亚洲精品| 国产成年人精品一区二区 | www.www免费av| 99久久人妻综合| 黄色丝袜av网址大全| 亚洲中文av在线| 757午夜福利合集在线观看| 男人操女人黄网站| 欧美在线黄色| 窝窝影院91人妻| а√天堂www在线а√下载| 亚洲 欧美 日韩 在线 免费| 夜夜夜夜夜久久久久| 久久中文看片网| 午夜福利影视在线免费观看| 日韩高清综合在线| 久久人人精品亚洲av| 国产在线精品亚洲第一网站| 男女之事视频高清在线观看| 老司机靠b影院| 校园春色视频在线观看| 国产精品日韩av在线免费观看 | 国产精品国产av在线观看| 国产亚洲精品综合一区在线观看 | 亚洲av成人不卡在线观看播放网| 久久久国产一区二区| 国产精品日韩av在线免费观看 | 一进一出好大好爽视频| 欧美乱码精品一区二区三区| 国产成人影院久久av| 乱人伦中国视频| 男女高潮啪啪啪动态图| 国产av精品麻豆| 欧美成人午夜精品| 国产成人欧美在线观看| 欧美黑人欧美精品刺激| 欧美日韩瑟瑟在线播放| 日本一区二区免费在线视频| 丰满的人妻完整版| 国产亚洲精品综合一区在线观看 | 国产深夜福利视频在线观看| 免费av中文字幕在线| 黑丝袜美女国产一区| 激情视频va一区二区三区| 91麻豆av在线| 成人免费观看视频高清| 一边摸一边做爽爽视频免费| 国产单亲对白刺激| www.999成人在线观看| 欧美不卡视频在线免费观看 | 日韩精品免费视频一区二区三区| 一区二区三区国产精品乱码| 午夜福利在线观看吧| 国产野战对白在线观看| 精品国产乱子伦一区二区三区| 狠狠狠狠99中文字幕| 每晚都被弄得嗷嗷叫到高潮| 大码成人一级视频| 琪琪午夜伦伦电影理论片6080| 国产亚洲av高清不卡| 国产成人欧美在线观看| 国产一区二区在线av高清观看| 欧美 亚洲 国产 日韩一| 久久久久久久久久久久大奶| 激情视频va一区二区三区| 亚洲国产精品sss在线观看 | 精品久久久久久久久久免费视频 | 女警被强在线播放| 女人被躁到高潮嗷嗷叫费观| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 50天的宝宝边吃奶边哭怎么回事| 色综合站精品国产| 中文字幕精品免费在线观看视频| 日韩av在线大香蕉| 午夜两性在线视频| 亚洲色图 男人天堂 中文字幕| 成人国产一区最新在线观看| 无人区码免费观看不卡| 无遮挡黄片免费观看| 国产精品久久电影中文字幕| 黄片小视频在线播放| 中文字幕高清在线视频| 狠狠狠狠99中文字幕| 国产亚洲精品久久久久久毛片| 日韩有码中文字幕| 搡老岳熟女国产| 精品第一国产精品| 黑丝袜美女国产一区|