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

    基于線粒體控制區(qū)的中國南海海域卵形鯧鲹遺傳多樣性

    2017-07-12 16:19:40呂金磊楊喜書宮亞運
    海洋漁業(yè) 2017年3期
    關(guān)鍵詞:卵形控制區(qū)譜系

    呂金磊,章 群,楊喜書,宮亞運,曹 艷

    (暨南大學水生生物研究中心,教育部熱帶亞熱帶水生態(tài)工程研究中心,廣州 510632)

    基于線粒體控制區(qū)的中國南海海域卵形鯧鲹遺傳多樣性

    呂金磊,章 群,楊喜書,宮亞運,曹 艷

    (暨南大學水生生物研究中心,教育部熱帶亞熱帶水生態(tài)工程研究中心,廣州 510632)

    為探討中國南海重要經(jīng)濟魚類卵形鯧鲹(Trachinotus ovatus)的遺傳多樣性,測定了廣東閘坡、烏石、安鋪和廣西東興以及海南新盈等5個地理群體97 ind樣品的線粒體控制區(qū)5′端359 bp序列,發(fā)現(xiàn)47個變異位點,32個單倍型,總體呈現(xiàn)高單倍型多樣性(h=0.951)和高核苷酸(π=0.020 9)多樣性的特點。在鄰接樹和單倍型網(wǎng)絡(luò)圖中出現(xiàn)2個分化顯著但不存在明顯地理聚群的分支,推測二者的分化時間約為60~18萬年前(中更新世),可能是中更新世冰期海平面下降形成邊緣海而導(dǎo)致隔離,間冰期海平面上升后出現(xiàn)二次接觸。不同地理群體間的遺傳分化不顯著(Fst=-0.022 4~0.045 3),AMOVA分析也顯示97%以上的遺傳變異來源于群體內(nèi)個體間。卵形鯧鲹2個譜系及總體的核苷酸錯配圖呈現(xiàn)多峰,中性檢驗均為負值不顯著(P>0.05),表明都未經(jīng)歷過大規(guī)模的種群擴張,處于相對穩(wěn)定的狀態(tài)。

    卵形鯧鲹;中國南海;線粒體控制區(qū);遺傳多樣性

    卵形鯧鲹(Trachinotus ovatus)隸屬鱸形目(Perciformes)鲹科(Carangidae),俗稱金鯧[1];因其體表無鱗,無肌間刺且肉質(zhì)細嫩、味道鮮美,具有較高的食用和營養(yǎng)價值[2],而成為中國和東南亞等產(chǎn)地廣受消費者歡迎的名貴魚類,也是中國華南地區(qū)養(yǎng)殖規(guī)模最大、產(chǎn)量最高的主導(dǎo)品種之一[3-4]。目前養(yǎng)殖群體未經(jīng)選育,就已出現(xiàn)生長慢、成活率低、品質(zhì)差等種質(zhì)衰退的現(xiàn)象[4]。同時,隨著海洋環(huán)境的惡化和過度捕撈,南海漁業(yè)資源在不斷的衰減,卵形鯧鲹也不例外[5]。因此,亟待開展種質(zhì)資源的保護和良種選育研究,而了解經(jīng)濟魚類的遺傳背景,是良種選育和資源保護的基礎(chǔ)[6]。

    目前,國內(nèi)外卵形鯧鲹的研究報道多集中在人工養(yǎng)殖育種、胚胎發(fā)育、疾病診斷和治療等方面,種群遺傳方面的研究較少,僅有孫立元等[7]采用微衛(wèi)星標記分析海南、廣東2個育種群體的遺傳多樣性;彭敏等[8]利用AFLP技術(shù)分析比較了北部灣野生群體和養(yǎng)殖群體的遺傳多樣性;趙永貞等[9]采用微衛(wèi)星分子標記技術(shù)分析了南海4個群體的遺傳多樣性;吉磊等[10]分析了海南、廣東、福建3個養(yǎng)殖群體的微衛(wèi)星多態(tài)性等。與雙親遺傳的核基因組相比,線粒體DNA為母系遺傳,結(jié)構(gòu)簡單、無基因重組、進化速率快,有效群體大小僅為前者的1/4,是種群遺傳學研究的首選分子標記,也是核基因組遺傳分析的重要補充[11-12]。其中控制區(qū)(control region,CR)是線粒體基因組中進化速率最快的區(qū)段,常被用于分析魚類遺傳背景[13]。本研究測定了廣東陽江閘坡、廉江安鋪鎮(zhèn)、雷州烏石鎮(zhèn),海南臨高新盈鎮(zhèn),廣西東興等5個地區(qū)野生群體的線粒體控制區(qū)序列,以分析遺傳多樣性、種群遺傳結(jié)構(gòu)和動態(tài),旨在了解南海卵形鯧鲹的遺傳背景,為種質(zhì)資源的管理保護和合理開發(fā)利用提供科學依據(jù)。

    1 材料與方法

    實驗用的野生卵形鯧鲹樣品分別采集于廣東陽江閘坡鎮(zhèn)(YZP1~20)、廣東廉江安鋪鎮(zhèn)(YAP1~18)、廣東烏石鎮(zhèn)(YWS1~19)、海南臨高新盈鎮(zhèn)(QXY1~20)和廣西東興(GDX1~20)(圖1),保存于95%乙醇備用。DNA提取和PCR擴增條件參照樂小亮等[15]的方法,PCR擴增引物為本實驗室自行設(shè)計的L:5′-tggcttgaaaaaccaccgttg,R:5′-cctgaartaggaaccaaatgccag,將電泳檢測為條帶清晰明亮的擴增產(chǎn)物送至華大基因切膠純化,通過ABI-3730 DNA自動測序儀測序。

    圖1 卵形鯧鲹采樣分布圖(引自孫湘平[14]。A:中國沿岸流,B:南海暖流)Fig.1 Sam p ling sites for T.Ovatus(cited from SUN Xiang-ping[14].A:the coastal current of China;B:Nanhaiwarm current)

    利用MEGA 6.0[16]對測定序列進行校對與對位排列,計算堿基組成、多態(tài)位點、簡約信息位點、轉(zhuǎn)換與顛換比以及基于Kimura 2-Parameter模型的遺傳距離,構(gòu)建鄰接樹。利用Network 5.0[17]構(gòu)建單倍型網(wǎng)絡(luò)分布圖。通過DnaSP 5.1[18]計算單倍型數(shù)、單倍型多樣性(h)、核苷酸多樣性(π),以及遺傳分化系數(shù)(Fst)和基因流(Nm)值。通過Arlequin 3.5[19]進行AMOVA分析(analysis ofmolecular variance),評估遺傳變異在群體內(nèi)、群體間與組群間的分布;計算Tajima’s D和Fu’s Fs值與核苷酸不配對分析,分析種群歷史動態(tài)[20]。使用SAMOVA 2.0[21]進行空間分子變異分析(spatial analysis ofmolecular variance,SAMOVA),通過群體的地理位置及序列差異分析群體間的分化。

    2 結(jié)果與分析

    2.1 序列特征及遺傳多樣性

    在5個群體97個樣本的線粒體控制區(qū)5′端高變區(qū)359 bp序列中,沒有堿基的插入與缺失;T、C、A、G平均含量分別是30.6%、15.2%、37.4%、16.8%,表現(xiàn)出典型的反G偏倚,吻合脊椎動物線粒體DNA特征[22]。轉(zhuǎn)換與顛換比為14.15,表明序列變異未飽合,可用于種群遺傳分析[23]。47個變異位點(13.9%)中39個為簡約信息位點,定義32個單倍型,其中16個為多個群體共享。卵形鯧鲹整體呈現(xiàn)出高單倍型多樣性(h=0.951)與高核苷酸多樣性(π=0.020 9)的特點,其中新盈(h=0.953,π=0.021 8)和安鋪群體(h=0.961,π=0.021 9)的遺傳多樣性最高;烏石群體最低(h=0.895,π=0.018 4),詳見表1。

    表1 各采樣點的采樣數(shù)量及遺傳多樣性情況Tab.1 Samp le sizes and genetic diversity indices of T.ovatus in each sampling site

    2.2 卵形鯧鲹的遺傳分化

    鄰接樹與單倍型網(wǎng)絡(luò)圖拓撲結(jié)構(gòu)大體一致,都出現(xiàn)由不同地理來源個體混雜分布的2個譜系:Lineage A和Lineage B(單倍型網(wǎng)絡(luò)圖見圖2)。譜系間分化顯著(Fst=0.569,P<0.01),以控制區(qū)變異速率為3%~10%/百萬年[24]推算的分化時間約為60~18萬年前(中更新世)。5個地理群體的群體間和群體內(nèi)遺傳距離相近(0.019~0.023)(表2),遺傳分化系數(shù)Fst較?。?0.022 4~0.045 3,P>0.05)。將5個地理群體按照瓊州海峽東/西兩側(cè),大陸沿海/海南島,以及根據(jù)SAMOVA結(jié)果(新盈、烏石、安鋪/閘坡、東興)分組進行分子方差(AMOVA)分析,發(fā)現(xiàn)群體內(nèi)的變異比例均達到97%以上(表3),表明群體間遺傳分化不明顯,遺傳變異主要源于個體之間,地理分布的影響相對較小。

    圖2 基于線粒體控制區(qū)基因序列構(gòu)建的卵形鯧鲹單倍型網(wǎng)絡(luò)圖Fig.2 Parsimony network of hap lotypes of T.ovatus based on m tDNA control region sequences

    表3 不同分組的分子方差分析Tab.3 AMOVA analysis of various groups of different geographic populatons

    表2 群體內(nèi)遺傳距離(對角線)、群體間遺傳距離(對角線下)及群體間遺傳分化程度(Fst,對角線上)Tab.2 Genetic distancesw ithin(along diagonal)and among(below diagonal)populations and fixation index(Fst,above diagonal)among populations of T.ovatus

    2.3 卵形鯧鲹種群動態(tài)分析

    中性檢驗及吻合度檢驗結(jié)果顯示,Lineage A的Tajima’s D值和Fu’s Fs的值都為負值(-01.041,-3.824)不顯著(P>0.05);Lineage B的Tajima’s D值為正值(0.172)不顯著(P>0.05)、Fu’s Fs為負值(-1.739)不顯著(P>0.05);兩譜系的SSD值(0.014,0.038)和Rg值(0.029,0.072)小且不顯著(P>0.05)(表4);核苷酸不配對分析圖均為多峰(圖3)。卵形鯧鲹總體的Tajima’s D值(-0.616)和Fu’s Fs值(-6.562)都為負值且不顯著(P>0.05),SSD值(0.010)和Rg值(0.015)小且不顯著(P>0.05)。表明南海卵形鯧鲹總體及2個譜系都沒有出現(xiàn)種群擴張的情況。

    表4 南海卵形鯧鲹群體的中性檢驗及吻合度檢驗Tab.4 Neutrality test and goodness of fit test among T.ovatus populations in the South China Sea

    圖3 南海卵形鯧鲹不同譜系的核苷酸不配對分布圖Fig.3 M ismatch-distribution analysis of different lineage of T.ovatus in the South China Sea

    3 討論

    3.1 卵形鯧鲹的遺傳多樣性

    就線粒體序列分析而言,一般常用單倍型多樣性(h)和核苷酸多樣性(π)來表示遺傳多樣性的高低程度[25]。南海卵形鯧鲹線粒體控制區(qū)總體呈現(xiàn)出高單倍型多樣性和高核苷酸多樣性的特點(h=0.951,π=0.020 9),單倍型多樣性與同海域的真鯛(Pagrusmajor)(h=0.994)、黑鯛(Acanthopagrus schlegeli)(h=0.994)[26]以及多齒蛇鯔(Saurida tumbil)(h=0.957)[27]等相近;核苷酸多樣性與同海域的真鯛(π=0.022)相近,高于黑鯛(π=0.007 5)和多齒蛇鯔(π=0.007 4),表明目前南海卵形鯧鲹線粒體序列多樣性相對較高。一般認為海洋魚類出現(xiàn)高單倍型多樣性和高核苷酸多樣性是因為種群處于持續(xù)穩(wěn)定增長狀態(tài)或者2個出現(xiàn)相對分化的群體(或亞種)發(fā)生二次接觸[28]。由于南海海區(qū)卵形鯧鲹存在2個譜系,歷史上種群數(shù)量總體穩(wěn)定,推測較高的遺傳多樣性是譜系間二次接觸所導(dǎo)致。

    就南海海區(qū)卵形鯧鲹微衛(wèi)星分析而言,孫立元等[7]、吉磊等[10]認為育種、養(yǎng)殖群體的遺傳多樣性較高,而趙永貞等[9]則認為野生群體多樣性較低;彭敏等[8]利用AFLP分析,認為北部灣養(yǎng)殖群體和野生群體的多樣性都較低。本研究結(jié)果與已報道的不完全一致,可能原因有:1)遺傳標記的不同:線粒體控制區(qū)為母系遺傳;微衛(wèi)星和AFLP是核基因標記,雙親遺傳;二者進化速率不同[29],受趨同進化、選擇作用等的影響不同[30];同時,由于不同報道分析使用的核基因標記位點也不盡相同,難以直接比較。2)采樣群體與采樣點不同:已報道的研究樣品大多為養(yǎng)殖群體,部分來源或親本不明確。

    3.2 卵形鯧鲹的遺傳分化和種群動態(tài)分析

    南海卵形鯧鲹中存在2個分化時間約在60~18萬年前(中更新世)的譜系。歷史資料顯示,更新世冰期與間冰期旋回,對生物種群的數(shù)量和分布有著巨大的影響[31]。冰期海平面下降使得中國近海出現(xiàn)多個封閉的內(nèi)陸海[32],導(dǎo)致群體間基因交流中斷而出現(xiàn)分化。盡管如此,群體間Fst=-0.022 4~0.045 3(P>0.05);97%以上的變異來源于群體內(nèi)部,表明群體間無明顯分化,推測原因如下:1)間冰期海平面上升,不同海域得以連通。2)南海存在復(fù)雜洋流(季風漂流、黑潮分支、南海暖流等)及中國沿岸流,瓊州海峽的出現(xiàn)也為東西兩側(cè)的群體交流提供了通道,卵形鯧鲹的受精卵及剛孵化的仔稚魚不具備游泳能力,易被海流攜帶而擴散。3)成魚游動性強,有聚群索餌和產(chǎn)卵習性[33]。

    單倍型網(wǎng)絡(luò)圖呈現(xiàn)不規(guī)則圖形,沒有中心單倍型;中性檢驗和核苷酸不配對分析表明2個譜系及總體都沒有出現(xiàn)擴張的情況,與同一海域分布的鯔(Mugil cephalus)和梭魚(Chelon haematocheilus)情況相似[34-35],可能原因是南海地處冰期氣溫變化相對溫和的低緯度地區(qū),受冰期的影響相對較小。

    3.3 卵形鯧鲹種質(zhì)資源的保護和利用

    物種遺傳多樣性越高,對環(huán)境變化的適應(yīng)能力就越強,也可為遺傳育種提供更豐富的材料[6]。南海卵形鯧鲹雖較同海域的部分魚類表現(xiàn)出較高的線粒體遺傳多樣性,但作為南海重要的經(jīng)濟魚類,受過度捕撈、環(huán)境污染等因素影響,捕獲量日漸稀少[8];隨著卵形鯧鲹人工種苗生產(chǎn)的規(guī)?;嚓P(guān)研究的開展[36-37],養(yǎng)殖業(yè)快速興起,但養(yǎng)殖苗種已出現(xiàn)種質(zhì)衰退跡象[4],養(yǎng)殖逃逸也可能嚴重影響野生群體的多樣性[38]。因此,在充分利用其資源的同時不應(yīng)忽視卵形鯧鲹的保護工作。本文僅研究了5個地理群體,沒有涉及福建沿海等地理群體,采樣點的局限可能沒有完全反映出卵形鯧鲹的遺傳多樣性現(xiàn)狀;另外因群體中存在2個譜系,使得隨機采集得到的每個譜系的樣品數(shù)量不足。在今后的工作中,將增加樣品采集的地理范圍和每個地理群體的樣本量,綜合運用線粒體序列分析與核基因指紋技術(shù),以期全面了解南海卵形鯧鲹的遺傳背景,為種質(zhì)資源的保護和開發(fā)提供更為準確的科學依據(jù)。

    [1] 朱元鼎,張春霖,張有為,等.南海魚類志[M].北京:北京科學出版社,1962:392-394.ZHU Y D,ZHANG C L,ZHANG Y W,et al.Fishes in the South China Sea[M].Beijing:Beijing Science Press,1962:392-394.

    [2] 楊欣怡.網(wǎng)箱海養(yǎng)卵形鯧鲹肌肉營養(yǎng)品質(zhì)評價和風味物質(zhì)研究[D].上海:上海海洋大學,2016:1 -13.YANG X Y.Study on nutritional composition and flavor compounds of sea-cage cultured Trachinotus ovatu muscles[D].Shanghai:Shanghai Ocean University,2016:1-13.

    [3] 陳偉洲,徐鼎盛,王德強,等.卵形鯧鲹人工繁殖及育苗技術(shù)研究[J].臺灣海峽,2007,26(3):435-441.CHENW Z,XU D S,WANG D Q,et al.Study on the spawning and hatching technique for Trachinotus ovatus[J].Journal of Oceanography in Taiwan Strait,2007,26(3):435-441.

    [4] 區(qū)又君,吉 磊,李加兒,等.卵形鯧不同月齡選育群體主要形態(tài)性狀與體質(zhì)量的相關(guān)性分析[J].水產(chǎn)學報,2013,37(7):961-969.OU Y J,JIL,LIJE,et al.Correlation analysis of major morphometric traits and body weight of selective group at differentmonth ages of Trachinotus ovatus[J].Journal of Fisheries of China,2013,37(7):961-969.

    [5] 鞠海龍.南海漁業(yè)資源衰減相關(guān)問題研究[J].東南亞研究,2012(6):51-55.JU H L.Studies on the decrementof fishery resource in the South China Sea[J].Southeast Asian Studies,2012(6):51-55.

    [6] WARD R D.Genetics in fisheriesmanagement[J].Hydrobiologia,2000,420(1):191-201.

    [7] 孫立元,郭華陽,朱彩艷,等.卵形鯧鲹育種群體遺傳多樣性分析[J].南方水產(chǎn)科學,2014,10(2):67-71.SUN L Y,GUO H Y,ZHU C Y,et al.Genetic polymorphism of breeding populations of golden pompano(Trachinotus ovatus)[J].South China Fisheries Science,2014,10(2):67-71.

    [8] 彭 敏,陳曉漢,陳秀荔,等.卵形鯧鲹養(yǎng)殖群體與野生群體遺傳多樣性的AFLP分析[J].西南農(nóng)業(yè)學報,2011,24(5):1987-1991.PENG M,CHEN X H,CHEN X L,et al.Genetic diversity of wild and cultured Trachinotus ovatus populations by AFLP markers[J].Southwest ChinaJournal of Agricultural Sciences,2011,24(5):1987 -1991.

    [9] 趙永貞,陳秀荔,李詠梅,等.南海區(qū)卵形鯧鲹遺傳多樣性的研究[J].西南農(nóng)業(yè)學報,2014,27(4):1786-1790.ZHAO Y Z,CHEN X L,LI Y M,et al.Genetic polymorphism of Trachinotusovatus in Nanhaidistrict[J].Southwest China Journal of Agricultural Sciences,2014,27(4):1786-1790.

    [10] 吉 磊,區(qū)又君,李加兒.卵形鯧鲹3個養(yǎng)殖群體的微衛(wèi)星多態(tài)性分析[J].熱帶海洋學報,2011,30(3):62-68.JI L,OU Y J,LI J E.Genetic polymorphism of three cultured populations of golden pompano Trachinotus ovatus as revealed bymicrosatellites[J].Journal of Tropical Oceanography,2011,30(3):62 -68.

    [11] AVISE JC.Phylogeography:The history and formation of species[M].London Harvard University Press.2000:9-32.

    [12] CANINO M F,SPIES I B,LOWE S A,et al.Highly discordant nuclear and mitochondrial DNA diversity in Atka Mackerel[J].Marine and Coastal Fisheries:Dynamics,Management,and Ecosystem Science,2010(2):375-387.

    [13] RAVAOARIMANANA I B,TIEDEMANN R,MONTAGNON D,et al.Molecular and cytogenetic evidence for cryptic speciation within a rare endemic Malagasy lemur,the Northern sportive lemur(Lepilemur septentrionalis)[J].Molecular Phylogenetics&Evolution,2004,31(2):440-448.

    [14] 孫湘平.中國近海區(qū)域海洋[M].北京:海洋出版社,2006:97-98.SUN X P.China offshore area[M].Beijing:Ocean Press,2006:97-98.

    [15] 樂小亮,章 群,趙 爽,等.一種高效快速的魚類標本基因組DNA提取方法[J].生物技術(shù)通報,2010(2):202-204.YUE X L,ZHANG Q,ZHAO S,et al.A fast and efficientmethod for isolation of genomic DNA from fish specimens[J].Biotechnology Bulletin,2010(2):202-204.

    [16] TAMURA K,STECHER G,PETERSON D,et al.MEGA6:Molecular evolutionary genetics analysis version 6.0[J].Molecular Biology and Evolution,2013,30(12):2725-2729.

    [17] BANDELT H J,F(xiàn)ORSTER P,ROHL A.Medianjoining networks for inferring intraspecific phylogenies[J].Molecular Biology and Evolution,1999,16(1):37-48.

    [18] LIBRADO P,ROZAS J.DnaSP v5:A software for comprehensive analysis of DNA polymorphism data[J].Bioinformatics,2009,25(11):1451-1452.

    [19] EXCOFFIER L,LISCHER H E L.Arlequin v 3.5:A new series of programs to perform population genetics analyses under Linux and Windows[J].Molecular Ecology Resources,2010,10(3):564-567.

    [20] NEIM,TAJIMA F.DNA polymorphism detectable by restriction endonucleases[J].Genetics,1981,97(1):145-163.

    [21] DUPANLOUP I,SCHNEIDER S,EXCOFFIER L.A simulated annealing approach to define the genetic structure of populations[J].Molecular Ecology,2002,11(12):2571-2581.

    [22] KNIGHT A,MINDELL D P.Substations bias,weighting of DNA sequence evolution,and the phylogenetic positions of fea′s viper[J].Systematic Biology,1993,42(1):18-31.

    [23] KOICHIRO T,GLEN S,DANIEL P,etal.MEGA6:Molecular evolutionary genetics analysis version 6.0[J].Molecular Biology and Evolution,2013(30):2725-2729.

    [24] HOCHACHKA PW,MOMMSEN T.Biochemistry and molecular biology of fishes:Environmental and ecological biochemistry[J].Science,1993(2):1 -38.

    [25] BONIN A,NICOLE F,POMPANON F,et al.Population adaptive index:A new method to help measure intraspecific genetic diversity and prioritize populations for conservation[J].Conservation Biology 2007,21(3):697-708.

    [26] 曹 艷.基于線粒體控制區(qū)序列的中國沿海3種鯛科魚類遺傳多樣性分析[D].廣州:暨南大學,2016:22-60.CAO Y.Genetic diversity of 3 sparid species in coastal waters of China based on mitochondrial control region sequences[D].Guangzhou:Jinan University,2016:22-60.

    [27] 孫冬芳,董麗娜,李永振,等.南海北部灣海域多齒蛇鯔的種群分析[J].水產(chǎn)學報,2010,34(9):1389 -1393.SUN D F,DONG L N,LI Y Z,et al.Population analysis of Saurida tumbil in the Northern South China Sea[J].Journalof Fisheries of China,2010,34(9):1389-1393.

    [28] GRANT W S,BOWEN B W.Shallow population histories in deep evolutionary linages of marine fishes:Insights from sardines and anchovies and lessons for conservation[J].Journal of Heredity,1998,89(5):415-426.

    [29] UPHOLTW B,DAUID IB.Mapping ofmitochondrial DNA of individual sheep and goat:Rapid evolution in the D-loop region[J].Cell,1997(11):57-583.

    [30] LARMUSEAU M H,RAEYMAEKERS J A,HELLEMANSB,et al.Mito-nuclear discordance in the degree of population differentiation in a marine goby[J].Heredity,2010(105):532-542.

    [31] DYNESIUSM,JANSSON R.Evolutionary consequences of changes in species'geographical distributions driven by Milankovitch climate oscillations.[J].Proceedings of the National Academy of Sciences,2000,97(16):9115-9120.

    [32] WANG P.Response of western pacificmarginal seas to glacial cycles:Paleoceanographic and sedimentological features 1[J].Marine Geology,1999,156(1-4):5-39.

    [33] 陳傅曉,唐賢明,譚 圖,等.卵形鯧鲹深水網(wǎng)箱養(yǎng)殖風險對策分析[J].中國漁業(yè)經(jīng)濟,2011,29(4):145-150.CHEN F X,TANG X M,TAN T,et al.The waveresisted deep-water net-cage breeding of Trachinotus ovatus for risk prevention countermeasures and economic benefit analysis[J].Chinese Fisheries Economics,2011,29(4):145-150.

    [34] JAMANDRE B W,DURAND JD,TZENG W N. Phylogeography of the flathead mullet Mugil cephalus in the north-west Pacific as inferred from themtDNA control region.[J].Journal of Fish Biology,2009,75(2):393-407.

    [35] LIU J X,GAO T X,WU S F,et al.Pleistocene isolation in the North-Western Pacific marginal seas and limited dispersal in a marine fish,Chelon haematocheilus(Temminck&Schlegel,1845)[J].Molecular Ecology,2007(16):275-288.

    [36] 區(qū)又君,李加兒,蔡文超.卵形鯧鲹腎臟的顯微鏡和超顯微結(jié)構(gòu)的觀察[J].海洋漁業(yè),2015,37(5):434-441.OU Y J,LI J E,CAI W C.Observation on microstructure and ultrastructure of kidney in Trachinotus ovatus[J].Marine Fisheries,2015,37(5):434-441.

    [37] 程大川,郭華陽,馬振華,等.3月齡卵形鯧鲹形態(tài)性狀對體質(zhì)量的影響分析[J].海洋漁業(yè),2016,38(1):26-34.CHENG D C,GUO H Y,MA Z H,et al.Mathematical analysis of morphometric attribute effects on body weight for three-month-old Trachinotus ovatus[J].Marine Fisheries,2016,38(1):26-34.

    [38] 張全啟,徐曉斐,齊 潔,等.牙鲆野生群體與養(yǎng)殖群體的遺傳多樣性分析[J].中國海洋大學學報(自然科學報),2004,34(5):816-820.ZHANG Q Q,XU X F,QI J,et al.The genetic diversity of wild and farmed Japanese flounder populations[J].Periodical of Ocean University of China,2004,34(5):816-820.

    Genetic diversity of Trachinotus ovatus in the South China Sea inferred from M itochondrial DNA control region sequences

    LV Jin-lei,ZHANG Qun,YANG Xi-shu,GONG Ya-yun,CAO Yan
    (Institute of Hydrobiology,Jinan University;Engineering Research Center of Tropical and Subtropical Aquatic Engineering,Ministry of Education,Guangzhou 510632,China)

    Due to the delicious taste,high nutritional value and increasingmarket demand,Trachinotusovatus has become one of the most economically important maricultured fish in the South China Sea.With the dwindling wild resources,the development of large-scale cultivation and the degraded germplasm of cultivars etc.,the accurate definition of population structure should be assessed for bettermanagement and sustainable exploitation of T.Ovatus resource in the South China Sea.In present study,359 bp mtDNA control region sequences of 97 individuals of T.Ovatus from Zhapo,Anpu and Wushi in Guangdong Province,Xinying in Hainan Province and Dongxing in Guangxi Provincewere analyzed to study the genetic diversity of this species in the South China Sea.A total of 47 polymorphic loci were found,of which 39 sites were parsimonyinformative and 32 haplotypes were defined.The resultant haplotype diversity(h)and nucleotide diversity(π)were 0.951±0.008 and 0.020 9±0.000 8 respectively,indicating the high level of haplotype diversity and nucleotide diversity of T.Ovatus populations in north of the South China Sea.In Neighborjoining tree and TCS network,two lineages were detected,and individuals from various geographical sites scattered in both lineages,thus no geographical cluster was found.The deduced differentiation time was approximately 600~180 ka BP.Despite the presence of two co-existing lineages,the absence of genetic heterogeneity among the 5 sampling siteswas noted,as the genetic fixation was low and not significant(Fst=-0.022 4~0.045 3,P>0.05),and the genetic distances were small(0.020-0.023).The analysis of molecular variance(AMOVA)also revealed that most of genetic variation resided within populations(>97%).This pattern of harboring obvious lineage structure but no geographic structure might have been shaped by glaciation cycles in Pleistocene.During the ice age,the closure ofmarginal seas caused by deeply lowered sea levelmighthave prevented gene flow among seperate seas,resulting in lineages differentiation.In interglacial stage,the rising of sea level might have reconnected marginal seas,promoting the secondary contact of once-isolated lineages,as the warm current of the South China Sea and the China Coastal current might have carried away fish eggs and larvae of T.Ovatus,and the formation of feeding and spawning schools also promoted the gene flow.Both nucleotide mismatch distribution and neutrality tests revealed that the lineage A,lineage B,and all populations as a whole in the South China Sea hadn’t experienced population expansion.In conclusion,although the genetic diversity of T.Ovatus in the South China Sea is relatively high yet,in the contextof over exploitation ofwild resources,degraded the germplastof cultivars,and the potential negative effects of breeding escape on their genetic diversity,propermanagement and sustainable exploitation of this species are urgently deserved.In further study,more populations across the species’range should be analyzed to determine the potential population genetic structure and demographic history with both mitochondrial and nuclear information.

    Trachinotus ovatus;South China Sea;mtDNA control region;genetic diversity

    Q 349

    A

    1004-2490(2017)03-0241-08

    2017-03-28

    國家自然科學基金項目(41071034);中央高?;究蒲袠I(yè)務(wù)費專項資金項目(21613105)

    呂金磊(1990-),男,碩士研究生,研究方向為分子生態(tài)學及生物多樣性保護。

    章 群(1968-),副研究員。E-mail:tqzhang@jnu.edu.cn

    猜你喜歡
    卵形控制區(qū)譜系
    卵形的變化
    ——奇妙的蛋
    神族譜系
    百年大黨精神譜系的賡續(xù)與文化自信
    王錫良陶瓷世家譜系
    基于OMI的船舶排放控制區(qū)SO2減排效益分析
    海洋通報(2020年5期)2021-01-14 09:27:06
    管好高速建筑控制區(qū)
    中國公路(2017年18期)2018-01-23 03:00:42
    再論東周時期銅簠的譜系和源流
    東方考古(2017年0期)2017-07-11 01:37:50
    阿什河流域非點源污染優(yōu)先控制區(qū)識別
    山東日照劃定大氣污染物排放控制區(qū)
    汽車縱橫(2017年3期)2017-03-18 23:19:33
    卵形鯧鲹PPARα基因cDNA序列的克隆、組織表達及生物信息學分析
    成人欧美大片| 亚洲 欧美 日韩 在线 免费| 97碰自拍视频| 麻豆国产97在线/欧美 | 成人一区二区视频在线观看| 一区二区三区激情视频| 欧美黄色片欧美黄色片| 嫁个100分男人电影在线观看| 99热这里只有精品一区 | av福利片在线观看| 国产精品久久久av美女十八| 亚洲欧美一区二区三区黑人| 色av中文字幕| 欧美乱色亚洲激情| 久久国产乱子伦精品免费另类| 视频区欧美日本亚洲| 欧美色视频一区免费| 正在播放国产对白刺激| 精品不卡国产一区二区三区| 18禁国产床啪视频网站| 18禁裸乳无遮挡免费网站照片| 成人亚洲精品av一区二区| 国产成+人综合+亚洲专区| 欧美又色又爽又黄视频| or卡值多少钱| 99久久精品国产亚洲精品| 国产精品野战在线观看| 日韩欧美国产一区二区入口| 日韩成人在线观看一区二区三区| 国产私拍福利视频在线观看| 免费人成视频x8x8入口观看| 国产精品久久电影中文字幕| 亚洲国产精品合色在线| 久久亚洲精品不卡| a级毛片在线看网站| 亚洲欧美精品综合久久99| 好看av亚洲va欧美ⅴa在| 此物有八面人人有两片| 丝袜美腿诱惑在线| 欧美中文日本在线观看视频| 久久久久久久精品吃奶| 国产v大片淫在线免费观看| 9191精品国产免费久久| 精品久久久久久,| 亚洲性夜色夜夜综合| 热99re8久久精品国产| 国产精品电影一区二区三区| 97碰自拍视频| 亚洲av第一区精品v没综合| 国产精品亚洲av一区麻豆| 国产亚洲精品av在线| 97碰自拍视频| 亚洲av第一区精品v没综合| 久久草成人影院| 特级一级黄色大片| 欧美av亚洲av综合av国产av| 69av精品久久久久久| 99热这里只有精品一区 | 亚洲专区国产一区二区| 88av欧美| 国产麻豆成人av免费视频| 久久 成人 亚洲| 午夜福利高清视频| 精品国产乱码久久久久久男人| av片东京热男人的天堂| 又黄又粗又硬又大视频| av欧美777| 母亲3免费完整高清在线观看| 国产激情久久老熟女| 老熟妇乱子伦视频在线观看| 国产精品香港三级国产av潘金莲| 久久久久国内视频| 看免费av毛片| 婷婷精品国产亚洲av在线| 90打野战视频偷拍视频| 欧美一级毛片孕妇| 婷婷亚洲欧美| 99热只有精品国产| 精品一区二区三区四区五区乱码| 天堂av国产一区二区熟女人妻 | 最近视频中文字幕2019在线8| 最近在线观看免费完整版| 中文字幕久久专区| 女人高潮潮喷娇喘18禁视频| 成在线人永久免费视频| 美女扒开内裤让男人捅视频| 亚洲中文字幕日韩| 国产一区二区在线av高清观看| 久99久视频精品免费| 国产高清videossex| 精品不卡国产一区二区三区| av有码第一页| 老司机午夜福利在线观看视频| 欧美乱色亚洲激情| 国产精品一及| 一级作爱视频免费观看| 国产亚洲精品综合一区在线观看 | 色av中文字幕| 麻豆成人av在线观看| 五月伊人婷婷丁香| 久久性视频一级片| 国产1区2区3区精品| 此物有八面人人有两片| 69av精品久久久久久| 国产av一区在线观看免费| 亚洲人成网站在线播放欧美日韩| 国产精品一区二区三区四区久久| 最近在线观看免费完整版| 看片在线看免费视频| 国产精品久久久久久亚洲av鲁大| 性色av乱码一区二区三区2| 精品日产1卡2卡| 欧美日韩精品网址| 一二三四社区在线视频社区8| 操出白浆在线播放| 国产熟女xx| 国产在线观看jvid| 香蕉国产在线看| 中文字幕高清在线视频| 一级a爱片免费观看的视频| 国产高清有码在线观看视频 | 亚洲欧美日韩无卡精品| 岛国在线免费视频观看| 国产99久久九九免费精品| 伊人久久大香线蕉亚洲五| 少妇裸体淫交视频免费看高清 | 最好的美女福利视频网| 久久精品国产亚洲av高清一级| 亚洲第一电影网av| 999久久久国产精品视频| 嫩草影院精品99| 国产v大片淫在线免费观看| 久久国产精品影院| 女生性感内裤真人,穿戴方法视频| 99热6这里只有精品| 国产成+人综合+亚洲专区| 亚洲美女视频黄频| 午夜精品久久久久久毛片777| 人人妻人人看人人澡| 免费看a级黄色片| 香蕉丝袜av| 国产亚洲欧美98| 欧美日韩精品网址| 麻豆国产97在线/欧美 | 大型黄色视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放 | а√天堂www在线а√下载| 日韩欧美在线二视频| 在线观看www视频免费| 久久香蕉国产精品| 久久热在线av| 婷婷亚洲欧美| 精品欧美国产一区二区三| 一进一出抽搐动态| 欧美 亚洲 国产 日韩一| 两个人免费观看高清视频| 毛片女人毛片| 欧美精品啪啪一区二区三区| 国产精品一区二区三区四区久久| 国内精品久久久久久久电影| 又紧又爽又黄一区二区| a级毛片a级免费在线| 午夜福利成人在线免费观看| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| 黄色毛片三级朝国网站| 国产蜜桃级精品一区二区三区| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 国产91精品成人一区二区三区| 午夜精品久久久久久毛片777| 熟妇人妻久久中文字幕3abv| 亚洲精品色激情综合| e午夜精品久久久久久久| 国产精品99久久99久久久不卡| 国产成年人精品一区二区| 免费在线观看成人毛片| 中文亚洲av片在线观看爽| 久久久国产精品麻豆| 婷婷丁香在线五月| 亚洲av成人不卡在线观看播放网| 久久久久国内视频| netflix在线观看网站| 亚洲精品久久成人aⅴ小说| 制服诱惑二区| 亚洲一区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 97碰自拍视频| 免费看十八禁软件| 亚洲欧美日韩东京热| 国产成年人精品一区二区| 成年免费大片在线观看| 亚洲国产高清在线一区二区三| 欧美性猛交黑人性爽| 一卡2卡三卡四卡精品乱码亚洲| 在线观看www视频免费| 日韩欧美免费精品| 天堂av国产一区二区熟女人妻 | 18禁黄网站禁片午夜丰满| 国产高清videossex| 日本成人三级电影网站| 国产精品av久久久久免费| 久久伊人香网站| 男人舔奶头视频| 国产精品久久视频播放| 亚洲国产日韩欧美精品在线观看 | 91字幕亚洲| 大型av网站在线播放| 高潮久久久久久久久久久不卡| 久久人人精品亚洲av| 在线观看免费日韩欧美大片| 亚洲最大成人中文| 麻豆国产97在线/欧美 | 啦啦啦免费观看视频1| 国产熟女午夜一区二区三区| 国产伦在线观看视频一区| 亚洲精品色激情综合| 老司机午夜福利在线观看视频| 欧美日韩黄片免| 亚洲九九香蕉| 波多野结衣巨乳人妻| 国产主播在线观看一区二区| 午夜老司机福利片| 我的老师免费观看完整版| 亚洲精品中文字幕一二三四区| 欧美最黄视频在线播放免费| 一进一出好大好爽视频| 久久久久国产一级毛片高清牌| 国产一级毛片七仙女欲春2| 男女之事视频高清在线观看| 成人精品一区二区免费| 88av欧美| 亚洲天堂国产精品一区在线| 欧美在线黄色| 精品一区二区三区四区五区乱码| 亚洲成人久久性| 一个人观看的视频www高清免费观看 | 91麻豆av在线| 久久久久久九九精品二区国产 | 久久久久久人人人人人| 国产成人影院久久av| 一本大道久久a久久精品| 男女床上黄色一级片免费看| 校园春色视频在线观看| 国产精品影院久久| 日韩有码中文字幕| 亚洲欧美精品综合一区二区三区| 一个人免费在线观看电影 | 成人国产一区最新在线观看| 国产熟女午夜一区二区三区| 动漫黄色视频在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜理论影院| 色综合站精品国产| 精品国产乱子伦一区二区三区| 一区福利在线观看| 国产单亲对白刺激| 国产精品影院久久| tocl精华| 欧美绝顶高潮抽搐喷水| 久久精品综合一区二区三区| 91大片在线观看| 国产一区在线观看成人免费| 欧美日韩福利视频一区二区| 人妻夜夜爽99麻豆av| 国产精品 国内视频| 黑人操中国人逼视频| 在线免费观看的www视频| 国产单亲对白刺激| 嫩草影视91久久| 给我免费播放毛片高清在线观看| 国产精品免费视频内射| 99久久无色码亚洲精品果冻| 欧美激情久久久久久爽电影| 亚洲精品中文字幕一二三四区| 国产1区2区3区精品| 免费在线观看完整版高清| 91麻豆av在线| 精华霜和精华液先用哪个| 国产精品美女特级片免费视频播放器 | 老汉色∧v一级毛片| 麻豆av在线久日| 日本熟妇午夜| 亚洲av第一区精品v没综合| 人人妻人人看人人澡| 国产成年人精品一区二区| 日韩欧美在线二视频| 亚洲国产精品999在线| 18禁黄网站禁片午夜丰满| 精品国产超薄肉色丝袜足j| 夜夜爽天天搞| 国产黄片美女视频| 欧美黑人巨大hd| 欧美乱色亚洲激情| 制服丝袜大香蕉在线| 男男h啪啪无遮挡| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 午夜免费成人在线视频| 无人区码免费观看不卡| ponron亚洲| 亚洲成av人片免费观看| 亚洲熟女毛片儿| 欧美在线一区亚洲| 在线十欧美十亚洲十日本专区| 听说在线观看完整版免费高清| 成人18禁高潮啪啪吃奶动态图| 人妻夜夜爽99麻豆av| 最近最新中文字幕大全免费视频| 特级一级黄色大片| 99精品欧美一区二区三区四区| 免费在线观看成人毛片| 日本五十路高清| 在线观看日韩欧美| 亚洲中文字幕一区二区三区有码在线看 | 国产熟女午夜一区二区三区| 久久婷婷成人综合色麻豆| 久久精品成人免费网站| av福利片在线观看| 久久久久国内视频| 全区人妻精品视频| 俺也久久电影网| 一级黄色大片毛片| 麻豆av在线久日| 国产黄片美女视频| 中国美女看黄片| 午夜激情av网站| 国产成人精品无人区| 国产亚洲精品一区二区www| 琪琪午夜伦伦电影理论片6080| 久久久精品大字幕| 男女床上黄色一级片免费看| 啦啦啦观看免费观看视频高清| 国产精品一及| 精品一区二区三区av网在线观看| 亚洲精品美女久久久久99蜜臀| 两个人看的免费小视频| 亚洲,欧美精品.| 制服人妻中文乱码| 国产91精品成人一区二区三区| 国产精品永久免费网站| 国产av一区二区精品久久| 欧美另类亚洲清纯唯美| 国产精品香港三级国产av潘金莲| 日本三级黄在线观看| av欧美777| 欧美另类亚洲清纯唯美| 99久久精品热视频| 国产片内射在线| 久久精品夜夜夜夜夜久久蜜豆 | 欧美日本亚洲视频在线播放| 亚洲成a人片在线一区二区| 日韩欧美国产一区二区入口| 夜夜夜夜夜久久久久| 国产三级中文精品| 在线看三级毛片| 又粗又爽又猛毛片免费看| 久久香蕉激情| 波多野结衣高清无吗| 日韩有码中文字幕| 午夜福利在线在线| 夜夜夜夜夜久久久久| 国产三级在线视频| 欧美成人免费av一区二区三区| 精品久久久久久,| 50天的宝宝边吃奶边哭怎么回事| 19禁男女啪啪无遮挡网站| 国产av麻豆久久久久久久| 久久香蕉国产精品| 亚洲午夜理论影院| 一夜夜www| 无限看片的www在线观看| 少妇裸体淫交视频免费看高清 | 90打野战视频偷拍视频| 亚洲一区二区三区不卡视频| av视频在线观看入口| 亚洲国产精品合色在线| 最近最新免费中文字幕在线| 日本成人三级电影网站| 久久午夜亚洲精品久久| 久久久久久人人人人人| 在线看三级毛片| 午夜影院日韩av| 可以免费在线观看a视频的电影网站| 搡老岳熟女国产| 老司机福利观看| 国产三级在线视频| 欧美极品一区二区三区四区| 亚洲激情在线av| 国产在线精品亚洲第一网站| 国产av不卡久久| 国产1区2区3区精品| 一个人免费在线观看的高清视频| netflix在线观看网站| 久久人妻福利社区极品人妻图片| 亚洲自偷自拍图片 自拍| 欧美成人性av电影在线观看| 亚洲一区中文字幕在线| a在线观看视频网站| 精品福利观看| 18美女黄网站色大片免费观看| 亚洲熟女毛片儿| 日本一二三区视频观看| 欧美三级亚洲精品| 国产成人av教育| 最好的美女福利视频网| 精品欧美一区二区三区在线| 午夜a级毛片| 精品久久久久久久久久久久久| 久久久久久久久久黄片| 天堂av国产一区二区熟女人妻 | 成人欧美大片| 日韩欧美免费精品| 婷婷精品国产亚洲av| 麻豆av在线久日| 精品国产亚洲在线| 久久久精品国产亚洲av高清涩受| 搡老岳熟女国产| √禁漫天堂资源中文www| 他把我摸到了高潮在线观看| 淫妇啪啪啪对白视频| 国产一区二区在线av高清观看| 欧美高清成人免费视频www| 欧美日韩乱码在线| 在线a可以看的网站| 嫁个100分男人电影在线观看| 老汉色∧v一级毛片| 国产激情欧美一区二区| 视频区欧美日本亚洲| 国产伦一二天堂av在线观看| 亚洲国产欧美一区二区综合| 色哟哟哟哟哟哟| 午夜激情福利司机影院| 天天躁夜夜躁狠狠躁躁| 好男人在线观看高清免费视频| 亚洲熟女毛片儿| 色综合亚洲欧美另类图片| 黄色视频不卡| svipshipincom国产片| 亚洲美女黄片视频| 免费看十八禁软件| 亚洲最大成人中文| 国产一区二区三区在线臀色熟女| 国产成人av激情在线播放| 丰满的人妻完整版| 老司机福利观看| 色在线成人网| 亚洲欧美激情综合另类| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频| 国产亚洲av高清不卡| 女生性感内裤真人,穿戴方法视频| 在线播放国产精品三级| 午夜a级毛片| 欧美一级毛片孕妇| 在线视频色国产色| 成人国语在线视频| 国产三级中文精品| bbb黄色大片| 麻豆久久精品国产亚洲av| 日韩高清综合在线| 老司机在亚洲福利影院| 精品免费久久久久久久清纯| 啦啦啦韩国在线观看视频| 国产一区在线观看成人免费| 久99久视频精品免费| 一边摸一边做爽爽视频免费| 国产精品免费一区二区三区在线| 国产成人aa在线观看| 国产午夜福利久久久久久| 免费看十八禁软件| 国产乱人伦免费视频| 99国产精品一区二区蜜桃av| 国产亚洲精品一区二区www| or卡值多少钱| 看黄色毛片网站| 亚洲中文av在线| 亚洲 欧美一区二区三区| www.999成人在线观看| 久久久久精品国产欧美久久久| 亚洲片人在线观看| 亚洲男人的天堂狠狠| 在线观看一区二区三区| 午夜福利成人在线免费观看| 国产一区在线观看成人免费| 成人国语在线视频| avwww免费| 99热这里只有是精品50| 听说在线观看完整版免费高清| 国产乱人伦免费视频| 久久久国产成人免费| 免费看a级黄色片| www.熟女人妻精品国产| 视频区欧美日本亚洲| 国产精品国产高清国产av| 国产精华一区二区三区| 久久人人精品亚洲av| 国产高清videossex| 琪琪午夜伦伦电影理论片6080| 亚洲国产精品合色在线| 最近在线观看免费完整版| 日韩成人在线观看一区二区三区| 九色国产91popny在线| 天天添夜夜摸| 在线观看免费视频日本深夜| 狂野欧美白嫩少妇大欣赏| 欧美在线黄色| 成人一区二区视频在线观看| 欧美zozozo另类| 亚洲成人久久性| 91麻豆av在线| 亚洲熟妇熟女久久| 亚洲一区中文字幕在线| 在线观看日韩欧美| 国产三级中文精品| 成年女人毛片免费观看观看9| 久久精品国产清高在天天线| 欧美乱码精品一区二区三区| 国模一区二区三区四区视频 | 美女午夜性视频免费| 国内精品久久久久精免费| 男女之事视频高清在线观看| 不卡一级毛片| 欧美成人免费av一区二区三区| 亚洲熟妇中文字幕五十中出| 熟妇人妻久久中文字幕3abv| 制服诱惑二区| 中文字幕人成人乱码亚洲影| 高潮久久久久久久久久久不卡| 搡老妇女老女人老熟妇| 免费在线观看影片大全网站| 亚洲av电影不卡..在线观看| 亚洲男人天堂网一区| 欧美绝顶高潮抽搐喷水| 久久九九热精品免费| 亚洲精品一区av在线观看| 国产av一区在线观看免费| 国产一区二区三区在线臀色熟女| 久久香蕉精品热| 最新在线观看一区二区三区| 丰满人妻一区二区三区视频av | 亚洲人成伊人成综合网2020| 他把我摸到了高潮在线观看| 中文资源天堂在线| 青草久久国产| 国产真人三级小视频在线观看| 欧美最黄视频在线播放免费| 一个人观看的视频www高清免费观看 | 国产成人精品久久二区二区免费| 99久久无色码亚洲精品果冻| 波多野结衣高清无吗| av免费在线观看网站| 黄色片一级片一级黄色片| 国产蜜桃级精品一区二区三区| 亚洲av成人一区二区三| 黄片大片在线免费观看| 美女大奶头视频| 成人一区二区视频在线观看| 最近最新免费中文字幕在线| 日本五十路高清| 亚洲人成网站高清观看| 国产伦人伦偷精品视频| 少妇粗大呻吟视频| 99热这里只有精品一区 | 国产伦人伦偷精品视频| 12—13女人毛片做爰片一| 蜜桃久久精品国产亚洲av| 国产精品爽爽va在线观看网站| 亚洲熟女毛片儿| 可以免费在线观看a视频的电影网站| 长腿黑丝高跟| 色老头精品视频在线观看| 欧美日韩福利视频一区二区| 日韩欧美一区二区三区在线观看| 国产午夜精品论理片| 1024手机看黄色片| 久久精品aⅴ一区二区三区四区| 国产亚洲精品综合一区在线观看 | 亚洲精品一区av在线观看| 黄色a级毛片大全视频| 欧美在线一区亚洲| 99国产极品粉嫩在线观看| 高清在线国产一区| 我要搜黄色片| 我的老师免费观看完整版| 亚洲国产精品久久男人天堂| 99国产综合亚洲精品| 免费在线观看影片大全网站| 欧美国产日韩亚洲一区| 最近视频中文字幕2019在线8| 国产午夜福利久久久久久| 亚洲精品国产一区二区精华液| 亚洲人成77777在线视频| 日韩中文字幕欧美一区二区| 色av中文字幕| 天天添夜夜摸| 悠悠久久av| 中出人妻视频一区二区| 亚洲无线在线观看| 妹子高潮喷水视频| 亚洲国产精品久久男人天堂| 日韩欧美免费精品| 日韩免费av在线播放| 观看免费一级毛片| cao死你这个sao货| 人人妻人人澡欧美一区二区| 国产成人av激情在线播放| 亚洲国产欧美一区二区综合| 色在线成人网| 久久 成人 亚洲| 日本一二三区视频观看| 一区福利在线观看| 精品久久久久久久人妻蜜臀av| 99riav亚洲国产免费|