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

    中部太平洋大眼金槍魚(yú)(Thunnus obesus)種群遺傳結(jié)構(gòu)的初步研究*

    2014-04-16 03:29:17吳智超許強(qiáng)華許柳雄戴小杰朱江峰
    海洋與湖沼 2014年1期
    關(guān)鍵詞:控制區(qū)金槍魚(yú)中西

    吳智超 許強(qiáng)華, 2, 3, 4 許柳雄, 2, 3, 4 戴小杰, 2, 3, 4 朱江峰, 2, 3, 4

    ?

    中部太平洋大眼金槍魚(yú)()種群遺傳結(jié)構(gòu)的初步研究*

    吳智超1許強(qiáng)華1, 2, 3, 4①許柳雄1, 2, 3, 4戴小杰1, 2, 3, 4朱江峰1, 2, 3, 4

    (1. 上海海洋大學(xué)海洋科學(xué)學(xué)院 上海 201306; 2. 大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn)室 上海 201306; 3. 農(nóng)業(yè)部大洋漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站 上海 201306; 4. 國(guó)家遠(yuǎn)洋漁業(yè)工程技術(shù)研究中心 上海 201306)

    大眼金槍魚(yú)()是世界最重要的漁業(yè)資源之一, 經(jīng)濟(jì)價(jià)值較高, 捕撈壓力的加大使中部太平洋大眼金槍魚(yú)過(guò)度捕撈程度愈發(fā)嚴(yán)重。開(kāi)展種群遺傳結(jié)構(gòu)研究將為大眼金槍魚(yú)資源的可持續(xù)利用與開(kāi)發(fā)提供科學(xué)依據(jù)。本文以取自中西太平洋和中東太平洋共182尾大眼金槍魚(yú)個(gè)體為研究對(duì)象, 分析研究其線粒體DNA控制區(qū)第一個(gè)高變區(qū)(the first hypervariable region, HVR-1)457bp的序列, 共發(fā)現(xiàn)115個(gè)變異位點(diǎn), 定義了178種單倍型。序列多樣性分析結(jié)果顯示: 8個(gè)采樣點(diǎn)核苷酸多樣性在0.03074—0.04362之間, 單倍型多樣性在各個(gè)采樣點(diǎn)都較高, 變動(dòng)范圍在0.996—1.000之間; 分子方差分析顯示99.97%的遺傳變異性來(lái)自于種群內(nèi); 群體間的高基因交流值m和低分化指數(shù)st揭示中東太平洋與中西太平洋群體間基因交流頻繁, 不存在顯著的遺傳分化。

    大眼金槍魚(yú)(); 種群遺傳結(jié)構(gòu); 線粒體DNA; 控制區(qū)

    大眼金槍魚(yú)()隸屬于鱸形目(Perciformes)、鯖亞目(Scombroidei)、鯖科(Scombridae)、金槍魚(yú)屬(), 是世界最重要的漁業(yè)資源之一, 經(jīng)濟(jì)價(jià)值較高, 分布在太平洋、大西洋、印度洋的熱帶、亞熱帶和溫帶廣闊水域(趙傳絪等, 1983; 苗振清等, 2003), 在我國(guó)主要分布在東海和南海(苗振清等, 2002)。WCPFC(中西太平洋漁業(yè)委員會(huì))的評(píng)估報(bào)告指出, 2008年中西太平洋大眼金槍魚(yú)的漁獲量雖然從2007年的14.3萬(wàn)噸增至15.7萬(wàn)噸, 但其捕撈壓力已超越可持續(xù)生產(chǎn)的水平, 資源正在進(jìn)一步惡化(Fonteneau, 2005)。另外, 在海洋漁業(yè)中, 漁獲物經(jīng)常由多個(gè)原始群組成, 盡管這些原始群在種群形態(tài)上十分相似, 但是在種群遺傳結(jié)構(gòu)上還可能進(jìn)行種群再分。不加區(qū)別的捕撈, 將會(huì)導(dǎo)致一些微弱原始群的衰退或滅絕。弄清大眼金槍魚(yú)的種群遺傳結(jié)構(gòu)將有助于掌握其生活史、估算種群數(shù)量及資源變動(dòng), 從而為大眼金槍魚(yú)資源的可持續(xù)利用和開(kāi)發(fā)提供科學(xué)依據(jù)(Excoffier, 1992)。

    對(duì)大眼金槍魚(yú)種群的分子遺傳分析始于1962年Suzuki(1962)的工作, 之后利用分子標(biāo)記對(duì)大眼金槍魚(yú)種群遺傳結(jié)構(gòu)的研究層出不窮。Grewe和Hampton(1998)同時(shí)利用線粒體DNA標(biāo)記和微衛(wèi)星標(biāo)記對(duì)太平洋內(nèi)的大眼金槍魚(yú)進(jìn)行了系統(tǒng)研究, 兩種標(biāo)記的分析結(jié)果都認(rèn)為太平洋內(nèi)大眼金槍魚(yú)不存在遺傳分化。Alvarado-Bremer等(1998)和Chow等(1993, 2000)在利用PCR-RFLP技術(shù)研究大洋間大眼金槍魚(yú)遺傳結(jié)構(gòu)時(shí), 也證明了太平洋大眼金槍魚(yú)呈現(xiàn)單一的遺傳結(jié)構(gòu)。因此, 這之后的幾年間, 學(xué)界都公推這一結(jié)論。由于線粒體DNA的控制區(qū)進(jìn)化速率最快(Brown, 1983), 且能夠有效地檢測(cè)出傳統(tǒng)生物學(xué)方法無(wú)法識(shí)別的種群分化, 線粒體DNA控制區(qū)已成為群體遺傳學(xué)中最有效的分子標(biāo)記(Ravaoarimanana, 2004)。而在這其中, 線粒體DNA控制區(qū)第一個(gè)高變區(qū)(the first hypervariable region, HVR-1)標(biāo)記越來(lái)越多地成為研究魚(yú)類(lèi)種群遺傳結(jié)構(gòu)的首選(Falk, 2003)。Martinez等(2006)和Chiang等(2006)先后將HVR-1分子標(biāo)記應(yīng)用到大眼金槍魚(yú)種群遺傳結(jié)構(gòu)的研究中。Chiang等(2006)以西太平洋大眼金槍魚(yú)為研究對(duì)象, 證實(shí)西太平洋大眼金槍魚(yú)是一個(gè)隨機(jī)交配的種群, 但鄰接樹(shù)和最小跨度網(wǎng)絡(luò)揭示了兩種線粒體世系的存在。同樣, 王中鐸等(2012)利用線粒體控制區(qū)標(biāo)記分析認(rèn)為, 南海大眼金槍魚(yú)群體與東太平洋群體間存在較顯著的分化, 而與西太平洋群體分化不顯著。

    作為連接?xùn)|太平洋與西太平洋的海體, 中部太平洋的大眼金槍魚(yú)群體遺傳結(jié)構(gòu)如何呢?鑒于中部太平洋擁有非常重要的金槍魚(yú)漁場(chǎng), 且目前還沒(méi)有關(guān)于該海域大眼金槍魚(yú)群體遺傳結(jié)構(gòu)的研究, 本研究將利用線粒體DNA控制區(qū)HVR-1分子標(biāo)記, 對(duì)中部太平洋大眼金槍魚(yú)野生群體進(jìn)行種群分析, 以期為大眼金槍魚(yú)的資源保護(hù)與合理利用提供重要的理論依據(jù)。

    1 材料與方法

    1.1 樣本采集

    2009年10月至2011年2月分別由3艘商業(yè)漁船作業(yè)采集到大眼金槍魚(yú)182尾(表1), 其中中西太平洋5個(gè)采樣點(diǎn), 分別命名為WP1—WP5, 中東太平洋3個(gè)采樣點(diǎn), 分別命名為EP1—EP3(圖1)。剪取大眼金槍魚(yú)的肌肉組織樣本, 將其固定在裝有95%酒精的離心管并凍存于-20°C冰箱, 備用。

    1.2 DNA提取、PCR擴(kuò)增和測(cè)序

    大眼金槍魚(yú)肌肉樣本DNA的提取采用標(biāo)準(zhǔn)的苯酚-氯仿抽提法。電泳后用紫外分光光度計(jì)檢測(cè)DNA濃度和純度。

    用于擴(kuò)增線粒體控制區(qū)的引物(Alvarado-Bremer, 1998)為L(zhǎng)15998-PRO(5′-TACCCCAAACTCCC-A-A-AGCTA-3′)以及CSBDH(5′-TGAATTAGGAACC-AGATGCCAG-3′)。PCR擴(kuò)增反應(yīng)體系為20μL:1.0× PCR反應(yīng)緩沖液2μL; 基因組模板DNA(10μg/mL) 1μL; Taq酶(5U/μL)0.12μL; dNTP(2.5mmol/L)1.6μL; 引物(10μmol/L)各0.5μL; 剩余用ddH2O補(bǔ)足。PCR反應(yīng)在Eppendorf熱循環(huán)儀上進(jìn)行, 反應(yīng)條件為:95°C預(yù)變性3min; 95°C 30s, 55°C 30s, 72°C 1min, 共30個(gè)循環(huán); 72°C延伸10min; 4°C保存。PCR產(chǎn)物純化后由上海生工生物工程技術(shù)服務(wù)有限公司進(jìn)行測(cè)序。

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

    利用ClustalX1.83軟件對(duì)控制區(qū)序列進(jìn)行多重比對(duì), 并輔以手工校正。采用MEGA4.0軟件統(tǒng)計(jì)DNA序列的堿基組成, 通過(guò)鄰接法基于Kimura雙參數(shù)模型構(gòu)建系統(tǒng)發(fā)育樹(shù), 并用bootstrap的方法重復(fù)1000次進(jìn)行檢驗(yàn)。利用TCS1.21軟件構(gòu)建單倍型間的系譜關(guān)系, 簡(jiǎn)約上限為默認(rèn)值95%。單倍型多樣性指數(shù)()與核苷酸多樣性指數(shù)()由軟件DnaSP version 4.00來(lái)進(jìn)行計(jì)算。采用ARLEQUIN version3.0軟件分析種群間分化指數(shù)st, 通過(guò)1000次重抽樣來(lái)檢驗(yàn)群體間st值的顯著性。另外, 對(duì)大眼金槍魚(yú)遺傳結(jié)構(gòu)進(jìn)行的分子方差分析也在ARLEQUIN version3.0中完成?;蛄鱩根據(jù)公式M=1/4(1/st?1)得出。

    2 結(jié)果

    2.1 基因序列多樣性

    經(jīng)PCR擴(kuò)增后, 獲得182個(gè)長(zhǎng)為457bp的線粒體控制區(qū)HVR-1序列, 發(fā)現(xiàn)115個(gè)變異位點(diǎn), 變異率為25.16%, 定義了178個(gè)單倍型。序列中各堿基在DNA雙鏈中平均所占的比例為:A(腺嘌呤)=35.4%, T(胸腺嘧啶)=28.9%, C(胞嘧啶)=22.0%, G(鳥(niǎo)嘌呤)=13.7%, A/T堿基對(duì)含量64.3%明顯高于C/G堿基對(duì)含量35.7%。群體內(nèi)線粒體DNA控制區(qū)HVR-1單倍型多樣性指數(shù)()和核苷酸多樣性指數(shù)()見(jiàn)表2。單倍型多樣性指數(shù)在8個(gè)采樣點(diǎn)都很高, 平均值達(dá)0.9998, 核苷酸多樣性介于0.03074—0.04362之間, 最高值出現(xiàn)在EP3采樣點(diǎn), 最低值出現(xiàn)在WP5采樣點(diǎn)。

    表1 大眼金槍魚(yú)采樣位點(diǎn), 采樣時(shí)間及采樣數(shù)目

    Tab.1 The sampling locations, date and number of Thunnus obesus

    圖1 大眼金槍魚(yú)采樣位點(diǎn)

    182尾大眼金槍魚(yú)個(gè)體中共檢測(cè)到178種單倍型。除了單倍型49在采樣點(diǎn)EP2和EP3同時(shí)出現(xiàn)外, 其他單倍型均為對(duì)應(yīng)各采樣點(diǎn)特有。另外, 還有其他3種單倍型也都不是為一個(gè)個(gè)體所獨(dú)有, 分別是單倍型39、62以及131, 對(duì)應(yīng)的取樣點(diǎn)為EP2、EP1、WP5(表3)。

    表2 大眼金槍魚(yú)線粒體DNA控制區(qū)HVR-1的遺傳多樣性參數(shù)

    表3 各采樣位點(diǎn)的大眼金槍魚(yú)單倍型分布

    Tab.3 Haplotype distribution of Thunnus obesus in different localities

    2.2 系統(tǒng)發(fā)育樹(shù)和單倍型網(wǎng)絡(luò)的構(gòu)建

    基于178種線粒體DNA控制區(qū)HVR-1單倍型構(gòu)建的NJ系統(tǒng)發(fā)育樹(shù)結(jié)果顯示:各單倍型廣泛分布在單倍型鄰接樹(shù)上, 未形成明顯的分支, 不存在顯著的譜系結(jié)構(gòu)?;贖VR-1構(gòu)建的單倍型間系譜關(guān)系顯示:178種單倍型只形成單一的一個(gè)網(wǎng)絡(luò), 各單倍型間呈現(xiàn)高水平的平行演化(圖略)。

    2.3 群體遺傳結(jié)構(gòu)分析

    采用AMOVA分析評(píng)估大眼金槍魚(yú)的群體遺傳結(jié)構(gòu), 分析結(jié)果表明99.97%的變異組分來(lái)自種群內(nèi); 而只有0.03%的變異組分出現(xiàn)在種群間(表4)。利用st進(jìn)一步檢驗(yàn)大眼金槍魚(yú)群體的群體遺傳結(jié)構(gòu)。根據(jù)兩兩群體間的st值大小可以發(fā)現(xiàn):所有的st值都較小, 并且除了WP3群體與EP3群體間統(tǒng)計(jì)檢驗(yàn)極其顯著(=0.00901)外, 其他任何兩個(gè)群體間的統(tǒng)計(jì)檢驗(yàn)都不顯著。結(jié)果暗示采樣區(qū)域內(nèi)大眼金槍魚(yú)不存在顯著的遺傳結(jié)構(gòu)(表5)。除此之外, 8個(gè)采樣點(diǎn)間大眼金槍魚(yú)的基因交流m也非常頻繁, WP3群體與EP3群體間基因流最小, 為6.75(表5)。

    表4 大眼金槍魚(yú)種群的分子方差分析

    Tab.4 AMOVA analysis on Thunnus obesus populations

    表5 采樣區(qū)域間大眼金槍魚(yú)種群的st分析及基因流m

    Tab.5 Fst analysis and Nm values of Thunnus obesus populations between sampling regions

    st值位于斜線下方;m值位于斜線上方; *差異顯著(<0.05); **差異極其顯著(<0.01)。值采用1000次重復(fù)來(lái)進(jìn)行估計(jì)。

    3 討論

    Fonteneau等(2005)指出, 除了西太平洋大眼金槍魚(yú)過(guò)度開(kāi)發(fā)程度較小外, 東太平洋、印度洋以及大西洋大眼金槍魚(yú)都處于嚴(yán)重過(guò)度捕撈的程度。為了更好地進(jìn)行大眼金槍魚(yú)的管理和保護(hù)工作, 業(yè)內(nèi)將大西洋、太平洋和印度洋大眼金槍魚(yú)分別作為各自單一的群體, 并由四大國(guó)際性委員會(huì):ICCAT(大西洋金槍魚(yú)類(lèi)保護(hù)委員會(huì))、IATTC(美洲熱帶金槍魚(yú)委員會(huì))、WCPFC(中西太平洋金槍魚(yú)養(yǎng)護(hù)委員會(huì))以及IOTC(印度洋金槍魚(yú)委員會(huì))來(lái)進(jìn)行管理。中西太平洋擁有太平洋最大的金槍魚(yú)漁場(chǎng), 對(duì)該范圍及附近海域大眼金槍魚(yú)的群體遺傳學(xué)研究將為合理保護(hù)與利用該漁業(yè)資源提供理論依據(jù)。

    本研究182尾大眼金槍魚(yú)樣本中, 132尾取自中西太平洋, 50尾取自中東太平洋。利用線粒體DNA控制區(qū)HVR-1標(biāo)記, 本研究共鑒定出178種單倍型, 中西太平洋和中東太平洋分別具有131和48種單倍型, 與樣本數(shù)量也基本相同。此項(xiàng)結(jié)果表明中部太平洋海域大眼金槍魚(yú)擁有極高的單倍型多樣性指數(shù), 這與許多學(xué)者對(duì)其他大洋性洄游魚(yú)類(lèi)線粒體控制區(qū)片段特征的描述相吻合(Grant, 1998; Carlsson, 2004)。另外, 對(duì)樣本線粒體DNA控制區(qū)HVR-1片段的序列分析顯示A+T含量遠(yuǎn)遠(yuǎn)高于C+G含量, 這與之前學(xué)者關(guān)于線粒體DNA控制區(qū)堿基組成的研究也相一致(Brown, 1986)。

    分子方差分析顯示99.97%的變異組分來(lái)自種群內(nèi); 而只有0.03%的變異組分出現(xiàn)在種群間, 種群內(nèi)差異遠(yuǎn)遠(yuǎn)大于種群間差異。對(duì)兩兩群體間進(jìn)行st分析顯示除了WP3群體與EP3群體間統(tǒng)計(jì)檢驗(yàn)極其顯著(=0.00901)外, 其他任何兩個(gè)群體間的統(tǒng)計(jì)檢驗(yàn)都不顯著。而基于線粒體DNA控制區(qū)HVR-1單倍型構(gòu)建的NJ系統(tǒng)發(fā)育樹(shù)和單倍型間系譜關(guān)系網(wǎng)絡(luò)圖均顯示不存在顯著的譜系結(jié)構(gòu), 這些分析結(jié)果暗示大眼金槍魚(yú)在中部太平洋不存在顯著的遺傳結(jié)構(gòu)。根據(jù)Liu等(2006)對(duì)鳀魚(yú)的研究發(fā)現(xiàn), 由于海洋作為一個(gè)大的水體, 魚(yú)類(lèi)在其內(nèi)能相對(duì)自由的運(yùn)動(dòng)并且海洋環(huán)流能夠幫助輸送海洋魚(yú)類(lèi)的卵和幼體, 因此海洋魚(yú)類(lèi)群體的遺傳分化一般都比較低, 不易形成顯著的遺傳結(jié)構(gòu)。所以我們可以推測(cè):本研究中所采集的中西太平洋和中東太平洋樣本間較小的遺傳分化可以歸因于個(gè)體間頻繁的基因交流。

    而事實(shí)上, 在赤道附近盛行東北信風(fēng)和東南信風(fēng), 兩種方向不同的信風(fēng)驅(qū)使海水向兩側(cè)分流, 使得在東太平洋的赤道地區(qū)的冷水上翻來(lái)補(bǔ)充被分流的海水, 此為東太平洋的冷水區(qū)。相反, 赤道太平洋表層暖的海水也被輸送至西太平洋, 海平面上升, 熱量不斷積累, 此為西太平洋的暖池區(qū)(陳雪冬等, 2006)。一般冷暖水團(tuán)的相遇區(qū)由于生產(chǎn)力較高, 常形成天然的漁場(chǎng)。中部太平洋集合了這一特點(diǎn)(周甦芳等, 2004), 其兩側(cè)的水團(tuán)驅(qū)使中西太平洋和中東太平洋大眼金槍魚(yú)群體在產(chǎn)卵季節(jié)不斷涌入該區(qū)域進(jìn)行繁殖, 使得兩個(gè)群體間的基因交流較高。在Chiang等(2006)的研究中, 根據(jù)HVR-1區(qū)域單倍型構(gòu)建的鄰接樹(shù)和最小跨度網(wǎng)絡(luò)揭示了西太平洋兩種線粒體世系的存在。而本研究基于HVR-1單倍型構(gòu)建的NJ系統(tǒng)發(fā)育樹(shù)和單倍型間系譜關(guān)系圖均未表現(xiàn)出遺傳分化的跡象, 猜測(cè)這可能與西太平洋和中部太平洋分別具有獨(dú)立的大眼金槍魚(yú)漁場(chǎng), 其間的基因交流相對(duì)較少有關(guān)。另外, 由于本研究中, 中東太平洋海域擁有的樣本總量較少, 甚至個(gè)別采樣點(diǎn)采樣數(shù)目較低(EP2=14, EP3=12), 這可能影響了中東太平洋內(nèi)的遺傳分化結(jié)果。但總的來(lái)看, 我們可以認(rèn)為, 中部太平洋只存在一種隨機(jī)交配的大眼金槍魚(yú)群體, 種群內(nèi)部不存在顯著的遺傳分化。

    作為太平洋金槍魚(yú)漁場(chǎng)最大的集中地, 中部太平洋大眼金槍魚(yú)的管理和保護(hù)工作由兩個(gè)金槍魚(yú)管理組織(IATTC和WCPFC)共同進(jìn)行。本研究初步確定中西太平洋與中東太平洋大眼金槍魚(yú)無(wú)明顯遺傳分化, 這一結(jié)果支持將中部太平洋大眼金槍魚(yú)作為單一管理單元進(jìn)行管理; 它同時(shí)也提示管理中部太平洋大眼金槍魚(yú)的兩個(gè)金槍魚(yú)管理組織(IATTC和WCPFC)應(yīng)該彼此協(xié)調(diào), 制訂相對(duì)比較統(tǒng)一的管理策略。值得一提的是:本研究中的中東太平洋海域采樣樣本數(shù)量相對(duì)不足(表1), 需要在今后的研究中進(jìn)一步增加采樣點(diǎn)和樣本數(shù), 以期獲得更可靠的實(shí)驗(yàn)結(jié)果。

    王中鐸, 郭昱嵩, 顏云榕等, 2012. 南海大眼金槍魚(yú)和黃鰭金槍魚(yú)的群體遺傳結(jié)構(gòu). 水產(chǎn)學(xué)報(bào), 36(2): 191—201

    陳雪冬, 崔雪森, 2006. 衛(wèi)星遙感在中東太平洋大眼金槍魚(yú)漁場(chǎng)與環(huán)境關(guān)系的應(yīng)用研究. 遙感信息, (1): 25—28

    苗振清, 黃錫昌, 2002. 世界金槍魚(yú)漁業(yè)現(xiàn)狀分析. 浙江海洋學(xué)院學(xué)報(bào)(自然科學(xué)版), (4): 307—313

    苗振清, 黃錫昌, 2003. 遠(yuǎn)洋金槍魚(yú)漁業(yè). 上海: 上海科學(xué)技術(shù)文獻(xiàn)出版社, 2—30

    周甦芳, 沈建華, 樊偉, 2004. ENSO現(xiàn)象對(duì)中西太平洋鰹魚(yú)圍網(wǎng)漁場(chǎng)的影響分析. 海洋漁業(yè), 26(3): 167—172

    趙傳絪, 陳思行, 1983. 金槍魚(yú)類(lèi)和金槍魚(yú)漁業(yè). 北京: 海洋出版社, 1—25

    Alvarado-Bremer J R, Stequert B, Robertson N W, 1998. Genetic evidence for inter-oceanic subdivision of bigeye tuna () populations. Marine Biology, 132: 547—557

    Brown G G, Gadaleta G, Pepe G, 1986. Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA. Molecular Biology, 192: 503—511

    Brown W M, 1983. Evolution of animal mitochondrial DNAs. In: Nei M, Koehn R K eds. Evolution of Genes and Proteins. Sunderland, M A. Sinauer Associates, UK: 62—88

    Carlsson J, McDowell J R, Diaz-Jaimes P, 2004. Microsatellite and mitochondrial DNA analyses of Atlantic bluefin tuna () population structure in the Mediterranean Sea. Molecular Ecology, 13: 3345—3356

    Chiang H C, Hsu C C, Lin H D, 2006. Population structure of bigeye tuna () in the South China Sea, Philippine Sea and western Pacific Ocean inferred from mitochondrial DNA. Fisheries Research, 79: 219—225

    Chow S, Inoue S, 1993. Intra-and inter-specific restriction fragment length polymorphism in mitochondrial genes oftuna species. Bulletin of the National Research Institute of Far Seas Fisheries, 30: 207—224

    Chow S, Okamoto H, Miyabe N, 2000. Genetic divergence between Atlantic and Indo-Pacific stocks of bigeye tuna () and admixture around South Africa. Molecular Ecology, 9: 221—227

    Excoffier L, Smouse P E, Quattro J M, 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics, 131: 479—491

    Falk T M, Teugels G G, Abban E K, 2003. Phylogeographic patterns in populations of the black chinned tilapia complex (,) from coastal areas in West Africa: Support for the refuge zone theory. Molecular Phylogenetics and Evolution, 27: 81— 92

    Fonteneau A, Ariz J, Delgado A, 2005. A comparison of bigeye () stocks and fisheries in the Atlantic, Indian and Pacific oceans. ICCAT, 57(2): 41—66

    Grant W S, Bowen B W, 1998. Population histories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation. Heredity, 89: 415—426

    Grewe P M, Hampton J, 1997. An assessment of bigeye () population structure in the Pacific Ocean based on mitochondrial DNA and DNA microsatellite analysis. The 7th meeting of the western yellowfin tuna research group, Fuji, Nadi: 98—320

    Liu J X, Gao T X, Zhuang Z M, 2006. Late Pleistocene divergence and subsequent population expansion of two closly related fish species, Japanese anchovy () and Australian anchovy (). Molecular Phylogenetics and Evolution, 40: 712—723

    Martinez P, Gonzalez E G, Castilho R, 2006. Genetic diversity and historical demography of Atlantic bigeye tuna (). Molecular Phylogenetic Evolotion, (39): 404—416

    Ravaoarimanana I B, Tiedemann R, Montagnon R, 2004. Molecular and cytogenetic evidence for cryptic speciation within a rare endemic Malagasy lemur, the Northern Sportive Lemur (). Molecular Phylogenetics and Evolution, 31: 440—448

    Suzuki A, 1962. On the blood types of yellowfin and bigeye tuna. American Naturalist, 96: 239—246

    POPULATION GENETIC STRUCTURE OFIN CENTRAL PACIFIC OCEAN

    WU Zhi-Chao1, XU Qiang-Hua1, 2, 3, 4, XU Liu-Xiong1, 2, 3, 4, DAI Xiao-Jie1, 2, 3, 4, ZHU Jiang-Feng1, 2, 3, 4

    (1. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China; 2. Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai 201306, China; 3. Scientific Observing and Experimental Station of Oceanic Fishery Resources, Ministry of Agriculture, Shanghai 201306, China; 4. National Distant-water Fisheries Engineering Research Center, Shanghai 201306, China)

    Overexploitation of bigeye tuna in the Central Pacific Ocean has become a serious concern. We investigated the population genetic structure aiming at providing scientific references for the sustainable development of the resources. In this study, 182 bigeye tuna specimens were sampled from the Central Pacific Ocean. Population genetic diversity was evaluated by using the sequence data of the first HVR-1 (hypervariable region, HVR-1) of mitochondrial control region. In total, 115 variable sites were observed and 178 haplotypes were identified. Analysis of mtDNA HVR-1 sequences of bigeye tuna from eight localities with the nucleotide diversities ranged 0.03074-0.04362 and high haplotype diversities ranged 0.996-1.000. AMOVA test of bigeye tuna revealed that 99.97% of the genetic variation occurred within populations. Lowst(fixation index) values and high gene flow (m) values revealed in this investigation indicate no significant genetic differentiation and high rates of gene flow between the Central Western Pacific bigeye tuna populations and the Central Eastern Pacific bigeye tuna populations.

    bigeye tuna; population genetic structure; mitochondrial DNA; control region

    10.11693/hyhz20121021001

    * 上海市捕撈學(xué)重點(diǎn)學(xué)科專(zhuān)項(xiàng)基金項(xiàng)目, S30702號(hào); 國(guó)家863計(jì)劃, 2012AA092303號(hào)。吳智超, E-mail: wuzhichao1988@126.com

    許強(qiáng)華, 教授, E-mail: qhxu@shou.edu.cn

    2012-10-21,

    2013-01-18

    Q75

    猜你喜歡
    控制區(qū)金槍魚(yú)中西
    金槍魚(yú)淹死事件
    可以加速的金槍魚(yú)機(jī)器人
    軍事文摘(2021年22期)2022-01-18 06:22:00
    金槍魚(yú)與海豚
    基于OMI的船舶排放控制區(qū)SO2減排效益分析
    影響中西的波希戰(zhàn)爭(zhēng)
    國(guó)內(nèi)新型遠(yuǎn)洋金槍魚(yú)圍網(wǎng)漁船首航
    商周刊(2018年19期)2018-12-06 09:49:41
    管好高速建筑控制區(qū)
    阿什河流域非點(diǎn)源污染優(yōu)先控制區(qū)識(shí)別
    山東日照劃定大氣污染物排放控制區(qū)
    “味其道”與“理其道”(一)——中西詩(shī)與思比較談片
    琪琪午夜伦伦电影理论片6080| 99国产极品粉嫩在线观看| 777久久人妻少妇嫩草av网站| 国语自产精品视频在线第100页| netflix在线观看网站| 国产成人av激情在线播放| 成人亚洲精品av一区二区| 狂野欧美激情性xxxx| 欧美黄色淫秽网站| 亚洲av电影不卡..在线观看| 国产精品永久免费网站| 久久影院123| 成人亚洲精品av一区二区| 亚洲国产精品久久男人天堂| 日韩欧美国产在线观看| 午夜免费鲁丝| 国产精品98久久久久久宅男小说| 婷婷六月久久综合丁香| 99riav亚洲国产免费| 亚洲第一欧美日韩一区二区三区| 国产熟女xx| 国产精品综合久久久久久久免费 | 国产亚洲精品久久久久5区| 国产精品美女特级片免费视频播放器 | 亚洲中文日韩欧美视频| 亚洲全国av大片| 亚洲av美国av| 在线观看午夜福利视频| 精品第一国产精品| 老司机靠b影院| 午夜福利,免费看| 在线观看午夜福利视频| 国产在线精品亚洲第一网站| 老司机靠b影院| 长腿黑丝高跟| 99久久国产精品久久久| 精品卡一卡二卡四卡免费| 欧美一区二区精品小视频在线| www.999成人在线观看| 黄片大片在线免费观看| 黄色a级毛片大全视频| 国产精品亚洲一级av第二区| 午夜福利,免费看| 桃红色精品国产亚洲av| 亚洲欧美日韩高清在线视频| 999精品在线视频| 一二三四社区在线视频社区8| 黄片大片在线免费观看| 亚洲 欧美一区二区三区| 免费久久久久久久精品成人欧美视频| 又黄又粗又硬又大视频| 午夜免费成人在线视频| 欧美日韩乱码在线| 女人被狂操c到高潮| 搡老熟女国产l中国老女人| 国产av在哪里看| 国产精品爽爽va在线观看网站 | 成年女人毛片免费观看观看9| 搡老熟女国产l中国老女人| av中文乱码字幕在线| 在线天堂中文资源库| 在线观看免费日韩欧美大片| 波多野结衣巨乳人妻| 伊人久久大香线蕉亚洲五| 女生性感内裤真人,穿戴方法视频| 精品欧美国产一区二区三| 人人澡人人妻人| 在线观看www视频免费| 亚洲第一青青草原| 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www| 午夜免费成人在线视频| 久久久久久久精品吃奶| 国内久久婷婷六月综合欲色啪| 又紧又爽又黄一区二区| 国产精品一区二区三区四区久久 | 无人区码免费观看不卡| 欧洲精品卡2卡3卡4卡5卡区| 怎么达到女性高潮| 久久 成人 亚洲| 看免费av毛片| 色综合婷婷激情| 好看av亚洲va欧美ⅴa在| 热re99久久国产66热| 久久久久久久久中文| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美日韩另类电影网站| 亚洲一区二区三区色噜噜| 亚洲自拍偷在线| 久久精品aⅴ一区二区三区四区| 亚洲中文av在线| 国产成人精品久久二区二区91| www.www免费av| 精品日产1卡2卡| 欧美不卡视频在线免费观看 | 亚洲美女黄片视频| 一级,二级,三级黄色视频| 亚洲第一电影网av| 精品国产超薄肉色丝袜足j| 天堂动漫精品| 国内精品久久久久久久电影| 少妇熟女aⅴ在线视频| 国产不卡一卡二| 一级,二级,三级黄色视频| 最近最新中文字幕大全电影3 | 亚洲美女黄片视频| 一个人观看的视频www高清免费观看 | 麻豆久久精品国产亚洲av| 99在线人妻在线中文字幕| 欧美日韩福利视频一区二区| 一本大道久久a久久精品| 美女免费视频网站| 日韩高清综合在线| 久久国产亚洲av麻豆专区| 亚洲精品一区av在线观看| 国产亚洲精品久久久久5区| 亚洲男人的天堂狠狠| 一二三四在线观看免费中文在| 国产欧美日韩一区二区三| 青草久久国产| 久久国产亚洲av麻豆专区| 男女午夜视频在线观看| av网站免费在线观看视频| 午夜久久久久精精品| 国产午夜福利久久久久久| 久久久久久久午夜电影| 波多野结衣一区麻豆| 国产精品一区二区免费欧美| 国产精品乱码一区二三区的特点 | 国产亚洲av嫩草精品影院| 夜夜夜夜夜久久久久| 夜夜躁狠狠躁天天躁| 精品不卡国产一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 色在线成人网| 色综合亚洲欧美另类图片| 黄色视频,在线免费观看| 欧美成狂野欧美在线观看| 亚洲成a人片在线一区二区| 欧美另类亚洲清纯唯美| 亚洲国产日韩欧美精品在线观看 | 男男h啪啪无遮挡| 国产亚洲av高清不卡| 亚洲国产精品999在线| 国产亚洲欧美在线一区二区| 国产成人欧美在线观看| 国产一级毛片七仙女欲春2 | 老司机福利观看| 精品不卡国产一区二区三区| 日韩视频一区二区在线观看| 黄色丝袜av网址大全| 老熟妇仑乱视频hdxx| av在线天堂中文字幕| 欧美日韩亚洲国产一区二区在线观看| 丝袜在线中文字幕| 亚洲欧美一区二区三区黑人| 看免费av毛片| 亚洲片人在线观看| 99香蕉大伊视频| 欧美乱码精品一区二区三区| 久热这里只有精品99| 成年人黄色毛片网站| 精品卡一卡二卡四卡免费| 国产不卡一卡二| 成年人黄色毛片网站| 在线观看舔阴道视频| 欧美黑人精品巨大| 三级毛片av免费| 曰老女人黄片| 视频区欧美日本亚洲| 久久香蕉国产精品| 国产97色在线日韩免费| 久久香蕉精品热| 19禁男女啪啪无遮挡网站| 国产欧美日韩综合在线一区二区| 亚洲色图av天堂| 久久精品影院6| 女人被狂操c到高潮| 国产亚洲av高清不卡| 视频区欧美日本亚洲| 亚洲精品国产一区二区精华液| 999精品在线视频| 少妇裸体淫交视频免费看高清 | 国产精品久久久久久人妻精品电影| 成人18禁高潮啪啪吃奶动态图| 黑人巨大精品欧美一区二区蜜桃| 麻豆av在线久日| 精品国内亚洲2022精品成人| 在线视频色国产色| 制服诱惑二区| 97人妻天天添夜夜摸| 日韩三级视频一区二区三区| 国产激情久久老熟女| 身体一侧抽搐| 亚洲五月天丁香| 日日夜夜操网爽| 亚洲国产毛片av蜜桃av| 久久九九热精品免费| 欧美乱码精品一区二区三区| 老汉色av国产亚洲站长工具| 法律面前人人平等表现在哪些方面| 麻豆久久精品国产亚洲av| 18禁美女被吸乳视频| 91老司机精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲无线在线观看| 日本vs欧美在线观看视频| 97超级碰碰碰精品色视频在线观看| 这个男人来自地球电影免费观看| 成熟少妇高潮喷水视频| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 日本a在线网址| 国产精品一区二区免费欧美| 欧美一区二区精品小视频在线| 成人亚洲精品一区在线观看| 性少妇av在线| 亚洲最大成人中文| 精品国产美女av久久久久小说| 亚洲av熟女| 女人精品久久久久毛片| 天堂影院成人在线观看| 在线观看www视频免费| 啦啦啦 在线观看视频| 黄色片一级片一级黄色片| 国产1区2区3区精品| 亚洲精品一区av在线观看| 天天添夜夜摸| 欧美另类亚洲清纯唯美| 超碰成人久久| 夜夜看夜夜爽夜夜摸| 国产不卡一卡二| 色播亚洲综合网| 亚洲精品国产色婷婷电影| 欧美日韩黄片免| 亚洲av美国av| 国产成人av教育| 国产熟女xx| 午夜福利18| 高清黄色对白视频在线免费看| www.精华液| 校园春色视频在线观看| 男男h啪啪无遮挡| 欧美+亚洲+日韩+国产| 国产精品精品国产色婷婷| 免费不卡黄色视频| 日韩国内少妇激情av| 国产精品,欧美在线| 国产av又大| 91大片在线观看| 久久人人97超碰香蕉20202| 久久中文字幕人妻熟女| 国产区一区二久久| 国产精品国产高清国产av| av视频免费观看在线观看| 日本 欧美在线| 亚洲精品国产一区二区精华液| av免费在线观看网站| 亚洲男人天堂网一区| 午夜精品久久久久久毛片777| 一级毛片高清免费大全| 亚洲一区二区三区色噜噜| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| 女同久久另类99精品国产91| 国产成人精品在线电影| 日韩欧美在线二视频| 精品国产亚洲在线| 午夜福利一区二区在线看| 午夜精品久久久久久毛片777| 亚洲av成人一区二区三| 精品久久久久久,| 99国产精品一区二区蜜桃av| 欧美中文综合在线视频| 一边摸一边做爽爽视频免费| 男女床上黄色一级片免费看| 女人精品久久久久毛片| 亚洲国产高清在线一区二区三 | 国产精品免费视频内射| 美女高潮喷水抽搐中文字幕| 久久这里只有精品19| 757午夜福利合集在线观看| 人人妻人人爽人人添夜夜欢视频| 日韩有码中文字幕| 91精品国产国语对白视频| 女人高潮潮喷娇喘18禁视频| 嫩草影院精品99| 国产成人欧美在线观看| 亚洲中文字幕日韩| 天堂√8在线中文| 午夜精品在线福利| 少妇熟女aⅴ在线视频| 欧美激情 高清一区二区三区| 大香蕉久久成人网| 丝袜美足系列| 国产亚洲精品第一综合不卡| 狠狠狠狠99中文字幕| 国产片内射在线| 欧美不卡视频在线免费观看 | √禁漫天堂资源中文www| 少妇熟女aⅴ在线视频| 国产精品影院久久| 亚洲七黄色美女视频| 国产精品久久电影中文字幕| 老司机深夜福利视频在线观看| 久久天堂一区二区三区四区| 精品无人区乱码1区二区| 亚洲三区欧美一区| 成人三级黄色视频| 亚洲国产精品成人综合色| 成人精品一区二区免费| 国产av一区二区精品久久| 操美女的视频在线观看| 日韩中文字幕欧美一区二区| 免费女性裸体啪啪无遮挡网站| 日日夜夜操网爽| 一进一出好大好爽视频| 国产单亲对白刺激| 天天躁夜夜躁狠狠躁躁| 99riav亚洲国产免费| 无限看片的www在线观看| 一二三四社区在线视频社区8| 中文字幕色久视频| 欧美日韩黄片免| 母亲3免费完整高清在线观看| 男男h啪啪无遮挡| 亚洲成人国产一区在线观看| 国产97色在线日韩免费| 女人被躁到高潮嗷嗷叫费观| 99国产精品一区二区蜜桃av| 国产亚洲欧美在线一区二区| 欧美乱妇无乱码| 亚洲av五月六月丁香网| 黄色视频,在线免费观看| 亚洲在线自拍视频| 国产不卡一卡二| aaaaa片日本免费| 国产一区二区三区视频了| 亚洲国产欧美一区二区综合| 在线观看舔阴道视频| 天堂√8在线中文| www.精华液| 日韩精品青青久久久久久| 1024视频免费在线观看| 窝窝影院91人妻| 在线观看免费视频网站a站| 99在线视频只有这里精品首页| 久久精品亚洲精品国产色婷小说| 岛国视频午夜一区免费看| 操美女的视频在线观看| 中文亚洲av片在线观看爽| 免费观看人在逋| 我的亚洲天堂| 99re在线观看精品视频| 69精品国产乱码久久久| 国产成+人综合+亚洲专区| 亚洲第一青青草原| 老司机午夜福利在线观看视频| 老司机靠b影院| 国产精品 国内视频| 亚洲国产欧美网| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 午夜视频精品福利| 欧美成人一区二区免费高清观看 | 91字幕亚洲| 久久亚洲真实| 99国产精品99久久久久| 亚洲视频免费观看视频| www国产在线视频色| 成人av一区二区三区在线看| 国产高清有码在线观看视频 | 久久香蕉国产精品| 午夜久久久久精精品| 欧美+亚洲+日韩+国产| 国产亚洲精品久久久久5区| 精品国产乱子伦一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 中亚洲国语对白在线视频| 欧美最黄视频在线播放免费| 神马国产精品三级电影在线观看 | 一级a爱视频在线免费观看| 88av欧美| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久人人人人人| 电影成人av| 久久久久久久久免费视频了| 久久人人精品亚洲av| 十八禁人妻一区二区| 一级毛片精品| 亚洲七黄色美女视频| 人人妻人人澡人人看| 欧美国产精品va在线观看不卡| 精品国产美女av久久久久小说| 欧美色视频一区免费| 又黄又爽又免费观看的视频| 岛国视频午夜一区免费看| 女人精品久久久久毛片| 妹子高潮喷水视频| 黑丝袜美女国产一区| 两个人看的免费小视频| 波多野结衣高清无吗| 亚洲精品久久国产高清桃花| 亚洲色图综合在线观看| 亚洲 国产 在线| 亚洲五月天丁香| 怎么达到女性高潮| 欧美在线黄色| 午夜精品在线福利| 亚洲午夜精品一区,二区,三区| 日本免费a在线| 99国产精品一区二区三区| 18禁美女被吸乳视频| 亚洲午夜精品一区,二区,三区| 最近最新免费中文字幕在线| 亚洲 欧美 日韩 在线 免费| 一进一出抽搐动态| 国产成人系列免费观看| 国产精品永久免费网站| 长腿黑丝高跟| 欧美乱码精品一区二区三区| 亚洲黑人精品在线| 日韩中文字幕欧美一区二区| 国产国语露脸激情在线看| 此物有八面人人有两片| 久久精品国产99精品国产亚洲性色 | 99在线人妻在线中文字幕| 一级毛片高清免费大全| а√天堂www在线а√下载| av免费在线观看网站| 午夜福利视频1000在线观看 | 国产精品亚洲av一区麻豆| 国产又色又爽无遮挡免费看| 久久午夜综合久久蜜桃| 久久久久久国产a免费观看| 国产熟女xx| 久久久久久久精品吃奶| 精品国产超薄肉色丝袜足j| 国产精品美女特级片免费视频播放器 | 欧美成狂野欧美在线观看| 国产高清视频在线播放一区| 国产精品亚洲av一区麻豆| 国产高清激情床上av| 午夜久久久久精精品| 亚洲精品美女久久av网站| 99久久99久久久精品蜜桃| 久久人妻福利社区极品人妻图片| 国产精品98久久久久久宅男小说| 麻豆国产av国片精品| 自线自在国产av| 国产区一区二久久| 国产精品久久久人人做人人爽| 亚洲欧美日韩另类电影网站| 亚洲五月天丁香| 亚洲人成网站在线播放欧美日韩| 级片在线观看| 一二三四社区在线视频社区8| 麻豆一二三区av精品| 国产伦人伦偷精品视频| 国产亚洲精品久久久久久毛片| 波多野结衣一区麻豆| 好男人电影高清在线观看| 非洲黑人性xxxx精品又粗又长| 搡老岳熟女国产| 天堂影院成人在线观看| 亚洲无线在线观看| 亚洲av五月六月丁香网| 欧美色视频一区免费| 涩涩av久久男人的天堂| 高清在线国产一区| 久久久久久久午夜电影| 亚洲一区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 99国产极品粉嫩在线观看| 成人三级做爰电影| 91精品三级在线观看| 国产熟女xx| 亚洲av第一区精品v没综合| 亚洲人成网站在线播放欧美日韩| 午夜老司机福利片| 国产不卡一卡二| 精品久久久久久久久久免费视频| 99国产综合亚洲精品| 国产精品1区2区在线观看.| 制服诱惑二区| 久久国产乱子伦精品免费另类| 国产精品永久免费网站| 国产精品日韩av在线免费观看 | 午夜精品在线福利| 国产野战对白在线观看| 在线播放国产精品三级| 国产精品日韩av在线免费观看 | 男人操女人黄网站| 欧美av亚洲av综合av国产av| 国产日韩一区二区三区精品不卡| 成在线人永久免费视频| 女同久久另类99精品国产91| 极品人妻少妇av视频| 欧美黑人精品巨大| 久久国产亚洲av麻豆专区| 国产高清有码在线观看视频 | 日日夜夜操网爽| 精品免费久久久久久久清纯| 久久精品91蜜桃| 非洲黑人性xxxx精品又粗又长| 一区福利在线观看| 精品国内亚洲2022精品成人| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2 | 一边摸一边抽搐一进一出视频| 亚洲五月婷婷丁香| 最好的美女福利视频网| 91大片在线观看| 国产精品精品国产色婷婷| 在线观看免费视频日本深夜| 久9热在线精品视频| 国产欧美日韩一区二区三区在线| 国产极品粉嫩免费观看在线| 老熟妇仑乱视频hdxx| 18禁黄网站禁片午夜丰满| 成人国产一区最新在线观看| 国产精品,欧美在线| 性欧美人与动物交配| 亚洲精品一区av在线观看| 久久久久国产一级毛片高清牌| 男女午夜视频在线观看| 女人被狂操c到高潮| 国产在线观看jvid| 欧美中文综合在线视频| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 午夜老司机福利片| 亚洲三区欧美一区| 制服诱惑二区| 婷婷精品国产亚洲av在线| 91麻豆av在线| 巨乳人妻的诱惑在线观看| 成人国语在线视频| 国产真人三级小视频在线观看| 天堂影院成人在线观看| 亚洲精品粉嫩美女一区| 咕卡用的链子| 国产精品久久久久久人妻精品电影| 亚洲中文字幕日韩| 制服诱惑二区| 麻豆久久精品国产亚洲av| 一区在线观看完整版| 嫩草影视91久久| 日韩一卡2卡3卡4卡2021年| 国产三级在线视频| 在线av久久热| av在线播放免费不卡| 高清在线国产一区| 国产成年人精品一区二区| 人妻丰满熟妇av一区二区三区| 天堂动漫精品| 窝窝影院91人妻| 免费高清在线观看日韩| 欧美 亚洲 国产 日韩一| 大型av网站在线播放| 久久国产精品影院| 制服诱惑二区| 久久欧美精品欧美久久欧美| 欧美黄色淫秽网站| 在线观看一区二区三区| 黄色丝袜av网址大全| 亚洲熟妇中文字幕五十中出| 一级毛片高清免费大全| 欧美日韩亚洲综合一区二区三区_| 国产精品乱码一区二三区的特点 | 亚洲 欧美一区二区三区| 欧美在线一区亚洲| 久久欧美精品欧美久久欧美| 国产日韩一区二区三区精品不卡| 性色av乱码一区二区三区2| 成人av一区二区三区在线看| 久久伊人香网站| 黄色视频,在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 国产精品免费视频内射| e午夜精品久久久久久久| 亚洲第一青青草原| 女性被躁到高潮视频| 国产成人av教育| 欧美老熟妇乱子伦牲交| 无限看片的www在线观看| 一级作爱视频免费观看| 丰满的人妻完整版| 香蕉久久夜色| 男人操女人黄网站| 两个人免费观看高清视频| 美国免费a级毛片| 久久久国产成人免费| 黄色女人牲交| 欧美乱色亚洲激情| 男人操女人黄网站| 国产免费av片在线观看野外av| 国产伦一二天堂av在线观看| 亚洲色图 男人天堂 中文字幕| 久久亚洲精品不卡| 黑人操中国人逼视频| 精品欧美国产一区二区三| 日韩欧美一区视频在线观看| av电影中文网址| 别揉我奶头~嗯~啊~动态视频| 色综合亚洲欧美另类图片| 老司机深夜福利视频在线观看| 女警被强在线播放| 亚洲激情在线av| 欧美性长视频在线观看| 精品国产一区二区久久| 国产精品一区二区三区四区久久 | 亚洲专区中文字幕在线| 亚洲精品国产色婷婷电影| 亚洲精品中文字幕一二三四区|