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

    松江鱸魚野生群體遺傳多樣性的RAPD分析和SCAR標(biāo)記的轉(zhuǎn)化

    2012-12-25 06:41:04劉至治潘連德唐文喬耿云皓
    Zoological Research 2012年2期
    關(guān)鍵詞:種群特異性遺傳

    曾 珍, 劉至治, 潘連德, 唐文喬, 王 茜, 耿云皓

    (上海海洋大學(xué) 省部共建水產(chǎn)種質(zhì)資源發(fā)掘與利用教育部重點(diǎn)實(shí)驗(yàn)室, 上海 201306)

    松江鱸魚野生群體遺傳多樣性的RAPD分析和SCAR標(biāo)記的轉(zhuǎn)化

    曾 珍, 劉至治*, 潘連德*, 唐文喬, 王 茜#, 耿云皓#

    (上海海洋大學(xué) 省部共建水產(chǎn)種質(zhì)資源發(fā)掘與利用教育部重點(diǎn)實(shí)驗(yàn)室,上海201306)

    首先, 從294條10個(gè)堿基隨機(jī)引物中, 篩選出32條多態(tài)性引物, 對富春江、黃河、灤河和鴨綠江等4個(gè)松江鱸魚 (Trachidermus fasciatus)野生群體共120尾個(gè)體進(jìn)行RAPD分析。結(jié)果表明, 松江鱸魚野生群體的遺傳多樣性較豐富,其主要表現(xiàn)在:① 在擴(kuò)增得到的591個(gè)位點(diǎn)中, 有515個(gè)(87.14%)位點(diǎn)呈現(xiàn)多態(tài)性,群體間多態(tài)位點(diǎn)比率(P)的大小順序?yàn)椋焊淮航后w89.17%>黃河群體87.99%>鴨綠江群體86.63%>灤河群體83.25%。② 松江鱸魚群體間的Shannon信息指數(shù)(IT)和Nei’s遺傳多樣性指數(shù)(HT)分別在0.3393~0.3566和0.2157~0.2279間, 灤河群體的值較其他3群體稍低;若作為一個(gè)整體, 則總的 Shannon信息指數(shù)(IT)和Nei’s遺傳多樣性指數(shù)(HT)分別為0.3710±0.2153和0.2336±0.1643。③ 雖然群體間基因流值(Nm)在5.76103~19.84497間, 顯示各地理群體間存在程度不同的基因交流, 但分子方差分析(AMOVA)結(jié)果卻表明, 各群體間存在顯著(P<0.05)或極顯著(P<0.01)的遺傳分化。④ 聚類分析表明, 鴨綠江群體首先與黃河群體聚為一支, 再與富春江群體相聚, 最后與單獨(dú)一支的灤河群體聚類, 表明鴨綠江、黃河、富春江等3群體間的遺傳距離與彼此間的地理距離遠(yuǎn)近密切相關(guān), 而灤河群體與它們的遺傳距離較遠(yuǎn)。其次, 從獲得的S1225525bp、S1225605bp、S1225841bp、S1345695bp、S1345825bp等5個(gè)特異RAPD條帶中, 成功地由S1225605bp、S1225841bp條帶分別轉(zhuǎn)化出SCAR01560bp、SCAR02443bp的SCAR標(biāo)記。這兩個(gè)標(biāo)記的出現(xiàn)頻率, 在鴨綠江群體最高(96.67%和93.33%)、富春江群體其次(83.33%和90%)、黃河群體再其次(56.67%和66.67%)、灤河群體最低(13.33%和20%)。因此, SCAR01560bp、SCAR02443bp可作為鑒別松江鱸魚灤河群體與其他3群體的分子標(biāo)記。

    松江鱸魚; 遺傳多樣性; RAPD; SCAR標(biāo)記

    松江鱸魚 (Trachidermus fasciatus, Heckel)隸屬于鲉形目(Scorpaeniformes)杜父魚科(Cottidae)松江鱸魚屬(Trachidermus), 俗稱“四鰓鱸”等, 是一種降海溯河洄游性魚類(Hubei Institute of Hydrobiology Laboratory of Ichthyology, 1976)。它個(gè)體不大, 但自古被譽(yù)為佳品, 其分布曾被認(rèn)為僅局限于長江口咸淡水域及附近小河中(Hubei Institute of Hydrobiology Laboratory of Ichthyology, 1976);不過, 最近Wang (2006)報(bào)道, 該魚從日本南部、朝鮮西南部到中國東部沿海及相鄰淡水水域都有分布。在中國, 隨著經(jīng)濟(jì)的發(fā)展, 因氣候變化、環(huán)境污染及水利設(shè)施的修建等諸多因素的影響, 松江鱸魚資源日益枯竭、種群瀕危, 已于1988年被列為國家二級保護(hù)動物 (Wang & Liang,2008)。進(jìn)入21世紀(jì)以來, 先后有研究者從形態(tài)、生化遺傳到分子水平上對松江鱸魚的種質(zhì)資源進(jìn)行了分析(Liu et al, 2010; Wang et al, 2001, 2002; Wang & Liang, 2008; Xu et al, 2008, 2009), 但這些報(bào)道受限于野生個(gè)體難覓等原因, 研究結(jié)果尚難以全面、客觀地反映中國現(xiàn)有松江鱸魚野生種群的遺傳多樣性現(xiàn)狀; 另一方面, 關(guān)于松江鱸魚種群分化問題, 研究者的觀點(diǎn)也不盡一致(Liu et al, 2010; Wang et al, 2001, 2002; Wang & Liang,2008; Xu et al, 2008, 2009)。因此, 采用其他分子標(biāo)記方法, 從不同側(cè)面進(jìn)一步探討這些問題十分必要。

    隨機(jī)擴(kuò)增多態(tài)性DNA(random amplified polymorphic DNA, RAPD)技術(shù)是由Wiliams et al(1990)和Welsh et al(1990)提出的一種具有簡單、快速、高效等特點(diǎn)的分子標(biāo)記, 已被廣泛地應(yīng)用于生物的親緣關(guān)系、遺傳多樣性、群體遺傳、種質(zhì)鑒定等多方面的研究(Barman et al, 2003; Gao et al, 2007; Yan et al, 2005)。序列特定擴(kuò)增區(qū)域標(biāo)記(sequence characterized applied region, SCAR)是在RAPD技術(shù)基礎(chǔ)上發(fā)展起來的一種可靠、穩(wěn)定、可長期利用的標(biāo)記技術(shù), 因其具有穩(wěn)定性好、重現(xiàn)性高等特點(diǎn), 已被國內(nèi)外研究者廣泛應(yīng)用于大西洋鮭(Salmo salar)、虹鱒(Oncorhynchus mykiss)(Zhang et al, 2007)、團(tuán)頭魴(Megalobrama amblycephala) (Zou et al, 2005)、“新吉富”尼羅羅非魚(Oreochromis niloticus)(Li et al, 2010)、甌江彩鯉(Cyprinus carpiovar. color)(Lu et al, 2011)、奧利亞羅非魚(Oreochromis aureus)(Gao et al, 2008)等許多水產(chǎn)動物的群體間特異性分子標(biāo)記或品系鑒別。本研究采用RAPD-SCAR分析方法, 既從核基因組角度分析松江鱸魚的群體遺傳多樣性, 又力求尋找地理群體間的分子標(biāo)記, 為松江鱸魚種質(zhì)資源的保護(hù)與利用、種群鑒別提供基礎(chǔ)資料。

    1 材料與方法

    1.1 材料

    試驗(yàn)所用的野生松江鱸魚樣本, 分別采自浙江富春江(簡稱富春江群體, 用FR表示)、山東黃河入???簡稱黃河群體, 用YR表示)、河北灤河(簡稱灤河群體, 用LR表示)、遼寧鴨綠江(簡稱鴨綠江群體, 用YL表示)。每個(gè)群體均采集活魚30~40尾,剪取約0.2 g右胸鰭, 放入1.5 mL離心管中, 先用95%酒精保存48 h后, 換用75%酒精保存于4 ℃冰箱中備用。

    1.2 方法

    1.2.1 基因組DNA的提取 采用常規(guī)的“酚-氯仿”法, 參照文獻(xiàn)(Liu et al, 2009)的方法進(jìn)行。提取的DNA于?20 ℃冰箱保存?zhèn)溆谩?/p>

    1.2.2 RAPD-PCR擴(kuò)增 10 堿基隨機(jī)引物294條(引物編號為S11-S40、S122-S185、S1210-S1409), 購自上海生工生物工程技術(shù)服務(wù)有限公司。參照說明書配置濃度 5 μmol/L, ?20 ℃冰箱保存?zhèn)溆?。PCR反應(yīng)體系25 μL:包括含2.5 μL buffer、1 μL dNTPs (2.5 μmol/L), 引物2 μL、2μL基因組DNA、0.6 μLTaqDNA聚合酶(北京艾德萊生物科技有限公司, 5 U/μL)、16.9 μL去離子水(天根生化有限公司)。反應(yīng)條件為:94 ℃預(yù)變性3 min, 94 ℃變性45 s、36℃退火45 s、72 ℃延伸90 s, 45個(gè)循環(huán)后, 72 ℃延伸10 min。所有PCR 反應(yīng)均在Mastercycler ep gradient S (Eppendorf)型PCR儀上進(jìn)行。取10 μL PCR擴(kuò)增產(chǎn)物, 經(jīng)凝膠(含EB 0.5 μg/mL)電泳(5 V/cm)1.5 h后, 用BIO-RAD Gel Doc XR凝膠成像系統(tǒng)拍照和分析。

    1.2.3 數(shù)據(jù)分析 經(jīng)電泳獲得的RAPD 圖譜, 在同一電泳遷移位置上, 有 RAPD 擴(kuò)增條帶的計(jì)為1, 沒有的計(jì)為 0, 將每群體所有個(gè)體在全部位點(diǎn)上的譜帶式樣組成一個(gè)“0或 1”矩陣。在假定Hardy-Weinberg平衡下, 用POPGENE 32軟件進(jìn)行分析, 計(jì)算多態(tài)位點(diǎn)數(shù)、多態(tài)位點(diǎn)比例(P)、Nei’s 遺傳多樣性指數(shù)(HT)、Shannon 信息指數(shù)(IT)、群體間遺傳相似度及群體間遺傳距離等。根據(jù)群體間的遺傳距離, 用MEGA4.0 (Tamura et al, 2007)軟件包中的“UPGMA”程序?qū)?個(gè)群體進(jìn)行聚類分析, 用Arlequin 3.0(Excoffier et al, 2005) 軟件進(jìn)行分子方差分析(analysis of molecular variance, AMOVA), 計(jì)算遺傳分化指數(shù)(Fst)。

    1.2.4 RAPD產(chǎn)物的克隆 利用 DNA 快速回收試劑盒(北京艾德萊生物科技有限公司), 從1.8%瓊脂糖凝膠中回收純化RAPD特異片段。將純化片段連接到pMD19-T Vector載體(TaKaRa公司), 然后轉(zhuǎn)化入TOP10感受態(tài)細(xì)胞(天根生化有限公司)。根據(jù)藍(lán)白斑原理, 挑選白色菌落進(jìn)行搖菌, 用M13+/-引物做菌落PCR、篩選鑒定陽性克隆, 每個(gè)個(gè)體挑選3個(gè)陽性克隆進(jìn)行測序(上海美吉生物醫(yī)藥有限公司)。

    1.2.5 SCAR標(biāo)記的篩選 將RAPD擴(kuò)增產(chǎn)生的特異性條帶回收并克隆、測序。根據(jù)測序結(jié)果, 用Primer Premier3.0(Rozen & Skaletsky, 2000)軟件設(shè)計(jì)19~25 bp的正反引物, 先經(jīng)小樣本(每群體8個(gè)體)PCR篩選, 對仍能擴(kuò)增出目的條帶的引物再進(jìn)行大樣本(每群體30個(gè)體)驗(yàn)證。PCR 擴(kuò)增參數(shù)基本與RAPD分析保持一致, 但退火溫度稍做變動(表1)、循環(huán)次數(shù)減為35個(gè)。

    2 結(jié) 果

    2.1 RAPD擴(kuò)增結(jié)果

    本研究從294條10堿基隨機(jī)引物中, 篩選出32條擴(kuò)增結(jié)果穩(wěn)定、重復(fù)性好的引物(圖1), 對4個(gè)松江鱸魚群體進(jìn)行RAPD分析。獲得的總位點(diǎn)數(shù)為591個(gè), 其中多態(tài)位點(diǎn)數(shù)515個(gè)(87.14%), 平均每條引物產(chǎn)生18.46個(gè)位點(diǎn), 擴(kuò)增片段大小集中在200~2 500 bp間。松江鱸魚群體間的多態(tài)性位點(diǎn)比例(P)的大小順序?yàn)椋焊淮航后w(FR)89.17%>黃河群體(YR)87.99%>鴨綠江群體(YL)86.63%>灤河群體(LR)83.25%。

    表1 SCAR標(biāo)記的引物序列、退火溫度及擴(kuò)增片段大小Tab. 1 Primer sequence, annealing temperature and size of PCR band of SCAR markers

    圖1 引物S1255(A)和S1345(B)的擴(kuò)增結(jié)果Fig. 1 The amplification results of primer S1255 (A) and S1345(B)

    在4群體間, 共獲得5條特異性RAPD擴(kuò)增帶(圖1A,B, 箭頭所指)。其中, 引物S1255有3條特異性帶:條帶S1255525bp在FR、YL兩群體的出現(xiàn)率很低, 僅6.67%, 即只有2個(gè)體檢測到此位點(diǎn), 而在YR、LR群體卻較高, 分別為23.33%、76.67%; 與此相反, 條帶S1255605bp在FR、YL兩群體的出現(xiàn)率(90%、93.33%)最高, 而在YR群體為73.33%、在LR群體僅13.33%; 條帶S1255841bp在4群體間的出現(xiàn)頻率與S1255605bp類似, 大小順序?yàn)椋篩L (96.67%)>FR (90%)>YR (83.33%)>LR(13.33%)。引物S1345有2條特異性帶:S1345695 bp和S1345825 bp, 其出現(xiàn)率在YL群體均為100%, 在FR和YR兩群體均為93.33%, 而在LR群體中卻只有16.67%。

    2.2 群體遺傳多樣性與遺傳分化

    4個(gè)松江鱸魚群體間的Shannon信息指數(shù)(IT)、Nei’s遺傳多樣性指數(shù)(HT)分別在0.3393~0.3566和0.2157~0.2279間(表2)。比較而言, LR群體的值較其它3群體稍低。若將4群體看成一個(gè)整體, 則總?cè)后w的Shannon信息指數(shù)(IT)、Nei’s遺傳多樣性指數(shù)(HT)分別為0.3710±0.2153、0.2336±0.1643??梢?松江鱸魚群體的基因多樣性較高。

    表2 松江鱸魚群體的Shannon信息指數(shù)(IT)和Nei’ s遺傳多樣性指數(shù)(HT)Tab. 2 Shannon’s information index (IT)and the Nei’ s genetic diversity (HT) in the populations of Trachidermus fasciatus

    分子方差分析(AMOVA)結(jié)果表明(表3), 在遺傳變異的總方差中, 僅4.63%的差異來自群體間,而95.37%的差異存在于群體內(nèi)個(gè)體間。4群體總遺傳分化指數(shù)(Fst)為0.0463。由表4可知, 群體間的Fst值在0.02458~0.07986間, 其中LR群體與FR、YL、YR群體的Fst值都較大, YR、YL群體間的Fst值最小(0.02458), 群體間存在顯著(P<0.05)或極顯著(P<0.01)的遺傳分化。不過, 由基因流值(Nm)的大小可知, 各地理群體間均存在一定程度的基因交流,其中YR、YL群體間的Nm最強(qiáng)(19.84497), 而FR、LR群體Nm最弱(5.76103)(表4)。

    表3 松江鱸魚群體間遺傳差異的 AMOVA分析Tab. 3 AMOVA analyses in the populations of Trachidermus fasciatus

    表4 松江鱸魚群體間基因交流值(Nm, 對角線上方)與遺傳分化指數(shù) (Fst, 對角線下方)Tab. 4 The value of gene flow (Nm,above diagonal) and Genetic differentiation index (Fst,below diagonal) among the four populations of Trachidermus fasciatus

    2.3 聚類分析

    YR和YL群體間的遺傳相似度最大而遺傳距離最小, 分別為0.9919和0.0082; 相反, FR和LR群體間的遺傳相似度最小而遺傳距離最大, 分別為0.9757和0.0246(表5)。根據(jù)遺傳距離進(jìn)行聚類分析表明, YL群體與YR群體首先聚為一支, 然后與FR群體聚類, 最后與單獨(dú)為一支的LR群體相聚(圖2)。顯然, LR群體與其他3個(gè)群體間的遺傳距離較遠(yuǎn), 而YL、YR、FR群體間的親緣關(guān)系與彼此間地理距離的遠(yuǎn)近密切相關(guān)。

    表5 群體間的遺傳相似度(對角線上方)和遺傳距離(對角線下方)Tab. 5 Genetic identity (above diagonal) and genetic distance (below diagonal) within or among populations of Trachidermus fasciatus

    圖2 基于遺傳距離的4個(gè)松江鱸魚群體間的UPGMA聚類分析Fig. 2 Dendrogram o f the four populations of Trachidermus fasciatus by UPGMA based on genetic distance

    2.4 特異性RAPD條帶的克隆、測序及SCAR標(biāo)記的建立與驗(yàn)證

    對5個(gè)特異性RAPD條帶進(jìn)行測序后發(fā)現(xiàn), 它們的大小分別為520、600、840、690和820 bp左右, 與RAPD擴(kuò)增的條帶大小基本一致(圖1), 表明得到的序列是真實(shí)的目的片段。在這5個(gè)特異性條帶上成功設(shè)計(jì)了5對SCAR標(biāo)記引物, 但經(jīng)小樣本PCR驗(yàn)證后發(fā)現(xiàn), 只有S1255605bp和S1255841bp(GenBank序列號:JN714551、JN714552)仍能擴(kuò)增出特異性目的片段(表1, 圖3), 其余3個(gè)標(biāo)記(序列未上傳)的特異性消失。又因考慮引物中G、C堿基的含量, 特異性引物的上、下游位置未從序列的兩端開始, 而是移向了序列的中間部分, 故所得的SCAR標(biāo)記SCAR01560bp、SCAR02443bp的長度分別要比對應(yīng)的RAPD標(biāo)記短(表1, 圖3)。接著, 為進(jìn)一步檢測SCAR標(biāo)記的有效性, 在松江鱸魚4群體共120尾個(gè)體中進(jìn)行大樣本驗(yàn)證。結(jié)果表明, SCAR01560bp標(biāo)記在FR、YR和YL群體的出現(xiàn)頻率分別高達(dá)83.33%(25/30)、56.67%(17/30)及96.67% (29/30); 而在LR群體的30個(gè)體中, 僅4尾具有此標(biāo)記, 出現(xiàn)頻率很低(13.33%)。同樣, SCAR02443bp標(biāo)記在FR、YL兩群體的出現(xiàn)頻率也很高, 分別為90%(27/30)和93.33%(28/30); 在YR群體的出現(xiàn)頻率也較高, 為66.67%(20/30); 但在LR群體出現(xiàn)的頻率比較低, 僅20%(6/30)??梢? SCAR01560bp和SCAR02443bp可作為鑒別松江鱸魚LR群體與其他3群體的特異性分子標(biāo)記。

    圖3 S1255605bp和S1255841bp片段的 SCAR轉(zhuǎn)化Fig. 3 SCAR bands from S1255605bp and S1255841bp markersA: SCAR01560bp; B: SCAR02443bp。 M:Marker;1~4:FR群體(FR population); 5~8:YR群體(YR population); 9~12:LR群體(LR population); 13~16:YL群體(YL population)。

    3 討 論

    3.1 松江鱸魚群體遺傳多樣性與種群遺傳結(jié)構(gòu)

    現(xiàn)代遺傳學(xué)的觀點(diǎn)認(rèn)為, 遺傳多樣性的高低與一個(gè)物種的適應(yīng)能力、生存能力及進(jìn)化潛能密切相關(guān)。遺傳多樣性的下降, 可能導(dǎo)致物種對環(huán)境適應(yīng)能力的降低(Liu et al, 2009; Vrijenhoek,1994)。因此,了解瀕危物種的群體遺傳多樣性, 對其保護(hù)與管理策略的制定具有重要意義。關(guān)于松江鱸魚群體遺傳多樣性, 先后有過一些生化遺傳、分子遺傳等方面的報(bào)道。同工酶的分析結(jié)果表明, 鴨綠江流域天然種群的多態(tài)位點(diǎn)比例(10%)和平均雜合度(0.03), 在魚類中均低于平均值(Wang et al, 2002); Xu et al (2009)通過ISSR技術(shù), 發(fā)現(xiàn)松江鱸魚丹東、秦皇島群體的多態(tài)位點(diǎn)比例(44%)等遺傳指標(biāo)均較低。而Liu et al (2010)通過線粒體控制區(qū)序列的多態(tài)性分析發(fā)現(xiàn), 松江鱸魚丹東、盤錦、秦皇島群體的單倍型多樣指數(shù)(Hd)維持在一個(gè)較高水平(0.989~1.00)。這些結(jié)論的不一致, 可能與作者所用同工酶多態(tài)位點(diǎn) (2個(gè))、ISSR多態(tài)性引物數(shù)量(4個(gè))及群體樣本數(shù)量(6~15尾個(gè)體)偏少有關(guān)。本研究借助RAPD技術(shù), 使用32個(gè)多態(tài)性10堿基隨機(jī)引物, 從整個(gè)基因組水平上, 對分布于較大地理范圍里的4個(gè)松江鱸魚地理種群共120個(gè)體(每群體30尾樣本)進(jìn)行分析, 結(jié)果表明,鴨綠江群體(YL)與黃河群體(YR)、富春江群體(FR)的多態(tài)位點(diǎn)比例(P)在86.63%~89.17%間, 只灤河群體(LR)略低 (83.25%)。這些值雖低于野生大彈涂魚(Boleophthalmus pectinirostris, 93.02%~96.35%, Liu et al, 2009)、野生小黃魚(Larimichthys polyactis, 91.03%, Han et al, 2007), 但高于半滑舌鰨(Cynoglossus semilaevisGüntuer, 76%~80%, Han et al, 2007)。故此, 我們認(rèn)為野生松江鱸魚的遺傳多樣性尚較豐富。

    種群遺傳結(jié)構(gòu)是指遺傳多樣性在群內(nèi)、群間的分布, 即遺傳分化。種群遺傳結(jié)構(gòu)可反映種群大小、繁育系統(tǒng)、隔離程度和遷徙格局的時(shí)空動態(tài)。分化是種群偏離共同祖先所經(jīng)歷的進(jìn)化過程, 反映了種群間遺傳差異的大小和距離的遠(yuǎn)近。環(huán)境為種群分化提供了特定場所, 自然選擇使分化得以實(shí)現(xiàn)(Tang et al, 2011)。關(guān)于松江鱸魚種群分化問題, 先后有些零星的報(bào)道。從表型性狀上看, Wang et al (2008)發(fā)現(xiàn), 松江鱸魚鴨綠江、青龍河和富春江等3個(gè)群體的分化雖未達(dá)到亞種水平, 但鴨綠江群體與其它兩個(gè)群體的差異較大。在分子水平上, AFLP和ISSR(Xu et al, 2008; 2009)的分析結(jié)果表明, 松江鱸魚丹東群體、秦皇島群體間的分化不明顯。相反, 線粒體D-Loop區(qū)的序列分析結(jié)果則顯示, 這兩個(gè)群體間存在顯著的遺傳分化(Fst=0.1422,P<0.05), 而盤錦群體與它們間的分化卻不明顯(P>0.01) (Liu et al, 2010)。上述觀點(diǎn)的不盡一致, 應(yīng)與研究者所采用分析方法靈敏度的不一、樣本數(shù)量多少不等及采樣群體地理位置的遠(yuǎn)近等等有關(guān)。

    本研究發(fā)現(xiàn), 雖然松江鱸魚4個(gè)地理群體間的Fst值(0.02~0.08間)不是很大, 但這些群體間卻都存在顯著(P<0.05)或極顯著(P<0.01)的遺傳分化, 特別是LR與YL、LR與FR群體間的Fst值(0.05903和0.07986)較大, 揭示灤河群體與鴨綠江、富春江群體間的分化已達(dá)到Wright (1978)標(biāo)準(zhǔn)的中等程度。這一結(jié)果與Wang et al (2008)的形態(tài)研究結(jié)果相符, 也與Liu et al (2010)報(bào)道的線粒體分析結(jié)果部分的一致。不過需要指出的是, 本研究的4個(gè)群體間的基因交流值(Nm=5~20>4)都比較大, 特別是YL、YR、FR群體間的基因交流值更大(12~20), 表明它們間的基因交流較為頻繁, 這可能掩蓋了各群體間固有的部分遺傳差異, 從而影響了彼此間遺傳分化指數(shù)的大小。松江鱸魚繁殖期間及其幼魚在近海岸生長階段, 季節(jié)性海洋沿岸流增加了不同地理種群間的混雜與基因交流的機(jī)會, 特別是地理位置相對較近、同屬于黃海海區(qū)的YR和YL群體, 這種交流更易發(fā)生。

    在UPGMA聚類圖中(圖2), YL、YR、FR群體間的關(guān)系與彼此地理位置的遠(yuǎn)近相一致, 呈現(xiàn)出從北到南的分布格局, 而LR群體與它們的親緣關(guān)系都較遠(yuǎn), 顯示出明顯的遺傳差異。

    3.2 SCAR標(biāo)記的應(yīng)用

    由特異性的RAPD條帶轉(zhuǎn)化而來的SCAR標(biāo)記,克服了RAPD重復(fù)性欠佳的弱點(diǎn), 已被廣泛應(yīng)用于分子標(biāo)記輔助育種、種質(zhì)鑒定等方面(Li et al, 2010; Lu et al, 2011; Zou et al, 2005)。本研究由S1255605bp和S1255841bp兩條RAPD特異性條帶成功轉(zhuǎn)化出SCAR01560bp和SCAR02443bp的SCAR標(biāo)記。大樣本檢測表明, 它們的出現(xiàn)頻率在YL、FR群體都很高(83%~97%), 在YR群體比較高(56%~67%), 而在LR群體卻很低(13%~20%)。因此, SCAR01560bp和SCAR02443bp可作為鑒別灤河群體與其他3群體的分子標(biāo)記。至于這兩個(gè)標(biāo)記與什么性狀相關(guān)聯(lián),有待以后的研究來證實(shí)。有趣的是, 在每群體的30尾樣本中, YL、FR群體未出現(xiàn)這兩個(gè)標(biāo)記的樣本數(shù)(1~2和3~5個(gè)), 恰與LR群體出現(xiàn)這兩個(gè)標(biāo)記的樣本數(shù)(4~6個(gè))接近。若假設(shè)這些個(gè)體恰好在3群體間發(fā)生了互換, 則YL、FR群體出現(xiàn)與LR群體不出現(xiàn)標(biāo)記的頻率幾乎達(dá)到100%。LR群體的養(yǎng)殖性能較其他海區(qū)的差(Pan et al, 2010), 在4群體中的遺傳多樣性最低、UPGMA聚類樹中與其他群體的關(guān)系最遠(yuǎn)等, 都揭示該群體具有自身的遺傳結(jié)構(gòu)特點(diǎn), 推測渤海灣海區(qū)內(nèi)可能存在松江鱸魚的獨(dú)立產(chǎn)卵場, 值得深入調(diào)查。

    本研究的另外3個(gè)RAPD標(biāo)記, 經(jīng)克隆、測序后, 雖曾設(shè)計(jì)多對引物進(jìn)行驗(yàn)證, 但都未能成功轉(zhuǎn)化為SCAR標(biāo)記。類似的現(xiàn)象在鯉(Gao et al, 2007)、團(tuán)頭魴(Zou et al, 2005)、奧利亞羅非魚(Gao et al, 2008)的SCAR標(biāo)記轉(zhuǎn)化中都存在。究其原因可能是:(1)S1255525 bp、S1345695 bp和S1345825 bp片段在基因組中可能是高度重復(fù)序列, 在轉(zhuǎn)化為SCAR標(biāo)記時(shí)特異性消失; (2)若產(chǎn)生 RAPD 特異性位點(diǎn)發(fā)生在隨機(jī)引物的退火區(qū), 用隨機(jī)引物進(jìn)行 RAPD擴(kuò)增時(shí)發(fā)生錯(cuò)配, 退火不良而不能擴(kuò)增, 而隨機(jī)引物 3'端增加了 9~15個(gè)堿基后, 在非特異性區(qū)提供了足夠的同源性, 克服了最初的錯(cuò)配允許引物退火擴(kuò)增,使原有的特異性消失(Gao et al, 2007; Xu et al, 1995); (3)也可能是特異RAPD帶位于引物的 3'端或者與基因的遺傳距離較遠(yuǎn), 連鎖不緊密造成的(Gao et al, 2007)。

    接下來, 我們將著力尋找YL、YR、FR群體間的特異性分子標(biāo)記, 為松江鱸魚天然種質(zhì)資源庫的保護(hù)提供分子水平上的鑒定手段。

    致謝:遼寧省水產(chǎn)苗種管理局檢疫檢測中心唐作鵬研究員、丹東市水產(chǎn)技術(shù)推廣站梁天紅高級工程師、山東省墾利縣周輝先生和浙江省杭州市金良庭先生協(xié)助采集松江鱸魚樣品, 上海海洋大學(xué)王成輝教授和本校水產(chǎn)與生命學(xué)院的2009級李雪松、馬召騰幫助處理數(shù)據(jù), 在此一并致謝。

    Barman HK, Barat A, Yadav BM, Banerjee S, Meher PK, Reddy PVGK, Jana RK. 2003. Genetic variation between four species of Indian major carps as revealed by random amplified polymorphic DNA assay[J]. Aquaculture,213(1-4): 115-123.

    Excoffier L, Laval G, Schneider S.2005. Arlequin (version 3.0): An integrated software package for population genetics data analysis[J]. Evol Bioinform,1: 47-50.

    Gao FY, Ye X, Lu WX, Huang ZH, Zhang T. 2008. RAPD analysis and SCAR transformation of three populations of Oreochromis aureus[J]. J Shanghai Fish Univ,17(1): 22-27. [高風(fēng)英,葉星,盧邁新,黃樟翰,張庭. 2008. 3個(gè)奧利亞羅非魚群體遺傳多樣性的RAPD分析和SCAR標(biāo)記的轉(zhuǎn)化[J].上海水產(chǎn)大學(xué)學(xué)報(bào),17(1): 22-27.]

    Gao JS, Sun XW, Liang LQ. 2007. Conversion of RAPD markers linked to cold tolerance to SCAR markers in common carp[J]. J Dalian Fish Univ,22(3): 194-197. [高俊生,孫效文,梁利群. 2007.鯉抗寒性狀的RAPD標(biāo)記轉(zhuǎn)化為SCAR標(biāo)記的研究[J].大連水產(chǎn)學(xué)院學(xué)報(bào),22(3): 194-197.]

    Han ZQ, Zhuang ZM, Gao TX, Liu JX, Li YH, Wang ZY, Tang QS. 2007. Genetic diversity in Cynoglossus semilaevis by AFLP, RAPD and mtDNA markers[J]. J Fish Sci Chn,14(2): 192-200. [韓志強(qiáng),莊志猛,高天翔,劉進(jìn)賢,李玉暉,王志勇,唐啟升. 2007.半滑舌鰨DNA的群體遺傳變異[J].中國水產(chǎn)科學(xué),14(2): 192-200.]

    Hubei Institute of Hydrobiology Laboratory of Ichthyology. 1976. Fishes in Yangtze River[M]. Beijing:Science Press, 214-215. [湖北省水生生物研究所魚類研究室. 1976.長江魚類[M].北京:科學(xué)出版社, 214-215.]

    Li SF, Tang SJ, Cai WQ. 2010. RAPD-SCAR markers for genetically improved NEW GIFT Nile Tilapia (Oreochromis niloticus niloticus L.) and their application in strain identification[J]. Zool Res,31(2): 147-153. [李思發(fā),唐首杰,蔡完其. 2010.遺傳改良“新吉富” (NEW GIFT)尼羅羅非魚RAPD-SCAR標(biāo)記開發(fā)及其在品系鑒別中的應(yīng)用[J].動物學(xué)研究,31(2): 147-153.]

    Liu HL, Zhang Q, Tang YL, Yu FY, Zhou JY. 2010. Structure and genetic diversity of mtDNA D-Loop sequences among Trachidermus fasciatus stocks in Yellow Sea and Bohai Sea of China[J]. Mar Sci Bull,29(3): 283-288. [劉海林,章群,唐優(yōu)良,余帆洋,周佳怡. 2010.黃渤海松江鱸魚線粒體控制區(qū)結(jié)構(gòu)與序列多態(tài)性分析[J].海洋通報(bào),29(3): 283-288.]

    Liu ZZ, Yang JQ, Wang ZQ, Tang WQ. 2009. Genetic structure and population history of Beleophthalmus petinirostris in Yangtze River Estuary and its Southern Adjacent Regions[J]. Zool Res,30(1): 1-10. [劉至治,楊金權(quán),王正琦,唐文喬. 2009.長江口及其南部鄰近地區(qū)大彈涂魚種群遺傳結(jié)構(gòu)及種群歷史分析[J].動物學(xué)研究,30(1): 1-10.]

    Lu YP, Hu ZH, Wang CH, Xiang SP, Su XP, Yao ZL, Wang KW. 2011.SRAP and SCAR molecular marker related to body color of “whole red” patterns of Cyprinus carpio var. color[J]. Acta Hydrobiol Sin,35(1): 45-50. [呂耀平,胡則輝,王成輝,項(xiàng)松平,蘇小平,姚子亮,王凱偉. 2011.與“全紅”甌江彩鯉體色相關(guān)的SRAP及SCAR分子標(biāo)記[J].水生生物學(xué)報(bào),35(1): 45-50.]

    Pan LD, Cai F, Ma ZT, Jiang X. 2010. The resources conservation and population features of Trachidermus fasciatus of China[J]. Fish Sci Technol Inform,37(5): 211-215. [潘連德,蔡飛,馬召騰,蔣鑫. 2010.中國境內(nèi)松江鱸魚的種群特征以及資源保護(hù)[J].水產(chǎn)科技情報(bào),37(5): 211-215.]

    Rozen S, Skaletsky HJ. 2000. Primer3 on the WWW for generalusers and for biologist programmers[M] // Krawetz S, Misener S. Bioinformatics Methods and Protocols: Methods in Molecular Biology. New Jersey: Humana Press, 365-386.

    Tamura K, Dudley J, Nei M, Kumar S. 2007. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0[J]. Mol Biol Evol,24(8): 1596-1599.

    Tang SJ, Li SF, Cai WQ. 2011. Mitochondrial DNA analysis on genetic variation of wild, domesticated, and genetically selected populations of blunt snout bream (Megalobrama amblycephala)[J]. J Fish Sci Chn,18(3): 1-10. [唐首杰,李思發(fā),蔡完其. 2011.團(tuán)頭魴野生、馴養(yǎng)、選育3類遺傳生態(tài)群體遺傳變異的線粒體DNA分析[J].中國水產(chǎn)科學(xué),18(3): 1-10.]

    Vrijenhoek RC.1994. Genetic diversity and fitness in small populations[M] // Loeschcke V, Tomiuk J, Jian S K, eds. Conservation Genetics. Basel: Birkh?user, 37-53.

    Wang JQ, Cheng G, Tang ZP. 2001. The distribution of Roughskin Sculpin (Trachidermus fasciatus Heckel) in Yalu River Basin, China[J]. J Fudan Univ: Nat Sci,40(5): 471-475. [王金秋,成功,唐作鵬. 2001.鴨綠江流域中國境內(nèi)松江鱸的分布[J].復(fù)旦學(xué)報(bào):自然科學(xué)版,40(5): 471-475.]

    Wang JQ, Liang HZ. 2008. Populations identification in morphology of roughskin sculpin Trachidermus fasciatus in China[J]. Ocean Limnol Sin,39(4): 348-353. [王金秋,梁鴻子. 2008.基于形態(tài)性狀的松江鱸魚(Trachidermus fasciatus)種群鑒別[J].海洋與湖沼,39(4): 348-353.]

    Wang JQ, Shi C, Cheng G. 2002. Isozyme analysis for a wild population of Roughskin Sculpin, Trachidermus fasciatus, from the Yalu River Basin in China[J]. J Fudan Univ: Nat Sci,41(6): 688-693. [王金秋,石椿,成功. 2002.鴨綠江流域中國境內(nèi)松江鱸魚(Trachidermus fasciatus)天然種群的同工酶分析[J].復(fù)旦學(xué)報(bào):自然科學(xué)版,41(6): 688-693.]

    Wang YH. 2006. Notes on the scientific name, type-locality and geographical distribution of roughskin sculpin (Trachidermis fasciatus Heckel,1840)[J]. Mar Fish,28(4): 299-303. [王幼槐. 2006.關(guān)于淞江鱸學(xué)名和模式產(chǎn)地以及地理分布之探討[J].海洋漁業(yè),28(4):299-303.]

    Welsh J, McClelland M. 1990. Fingerprinting genomes using PCR with arbitrary primers[J]. Nucleic Acids Res,18(24): 7213-7218.

    Williams JGK, Kubelik AR, Livak KJ, Rafalski JA Tingey SV. 1990. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers[J]. Nucleic Acids Res,18(22): 6531-6535.

    Wright S. 1978. Evolution and the Genetics of Populations. Vol. 4: Variabilty Within and Among Natural Populations[M]. Chicago: University of Chicago Press.

    Xu H, Wilson DJ, Arulsekar S, Bakalinsky AT. 1995. Sequence-specific polymerase chain-reaction markers derived from randomly amplified polymorphic DNA markers for fingerprinting grape (Vitis) rootstocks[J]. J Amer Soc Hort Sci,120(5): 714-720.

    Xu JR, Han XL, Li N, Xu P, Bao ZM. 2008. Analysis of genetic diversity in roughskin sculpin Trachidermus fasciatus by AFLP markers[J]. J Dalian Fish Univ,23(6): 437-441.[徐建榮,韓曉磊,李寧,許璞,包振民. 2008.松江鱸群體遺傳多樣性的AFLP分析[J].大連水產(chǎn)學(xué)院學(xué)報(bào),23(6): 437-441.]

    Xu JR, Han XL, Yu JF, Bao F, Xu P. 2009. Analysis of genetic diversity of Trachidermus fasciatus by ISSR[J]. Freshwat Fish,39(1): 21-24. [徐建榮,韓曉磊,郁建鋒,鮑峰,許璞. 2009.松江鱸群體遺傳多樣性的ISSR分析[J].淡水漁業(yè),39(1): 21-24.]

    Yan JP, Liu SJ, Sun YD, Zhang C, Luo KK, Liu Y. 2005. RAPD and microsatellite analysis of diploid gynogens from allotetraploid hybrids of red crucian carp (Carassius auratus)× common carp ( Cyprinus carpio)[J]. Aquaculture,234(1-4): 49-60.

    Zhang JB, Wang HJ, Cai ZP. 2007. The application of DGGE and AFLP-derived SCAR for discrimination between Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss)[J]. Food Control,18(6): 672-676.

    Zou SM, Li SF, Cai WQ. 2005. SCAR transformation of a RAPD marker in blunt snout bream “Pujiang No. 1”[J]. J Fish Chn,29(3): 296-299. [鄒曙明,李思發(fā),蔡完其. 2005.團(tuán)頭魴“浦江1號”一個(gè)RAPD標(biāo)記的SCAR轉(zhuǎn)化[J].水產(chǎn)學(xué)報(bào),29(3): 296-299.]

    Analysis of genetic diversity in wild populations ofTrachidermus fasciatusby RAPD and the transformation of two SCAR markers

    ZENG Zhen, LIU Zhi-Zhi*, PAN Lian-De*, TANG Wen-Qiao, WANG Qian#, GENG Yun-Hao#

    (Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources Certificated by Ministry of Education, Shanghai Ocean University, Shanghai 201306, China)

    Firstly, RAPD was conducted to analyze genetic diversity ofTrachidermus fasciatusin the Fuchun River population (FR), Yellow River population (YR), Luan River population (LR), and Yalu River population (YL), with 32 polymorphic 10-bp random primers selected from 294 ones. Thirty wild individuals were detected in each population. The results indicated that the genetic diversity ofT.fasciatuswas relatively rich. The major results were as the following: 1) Altogether, 591 bands were detected and 515 of them were polymorphic, accounted for 87.14%. The range of proportion of polymorphic loci (P) was: FR(89.17%)>YR(87.99%)>YL(86.63%)>LR(83.25%). 2) The Shannon’s information index(IT) and Nei’s genetic diversity(HT) among populations were 0.3393?0.3566 and 0.2157?0.2279, respectively. Compare to other three populations, LR population had relative lower values. If took the populations as a whole, the total Nei’s genetic diversity(HT) and Shannon’s information index(IT) was 0.2336±0.1643 and 0.3710±0.2153, respectively. 3) The value of gene flow (Nm) (5.76103?19.84497) were high, indicating certain gene exchange existedamong the four populations. But the AMOVA results exhibited significantly differentiation (P<0.05 orP<0.01) among the populations. 4) In the UPGMA tree constructed according to genetic distance, YL and YR populations clustered firstly, then with FR population, and finally they joined to LR population. Obviously, the YL, YR and FR populations had relatively close relationship according to their geographic distance, whereas LR population showed clear divergence to the other three populations. Secondly, out of the five special RAPD bands (S1225525bp, S1225605bp, S1225841bp, S1345695bpand S1345825bp), SCAR maker SCAR01560bpand SCAR02443bpwere successfully transformed from S1255605bpand S1255841bp, respectively. After large samples examination of the two markers, we found the highest frequency (96.67% and 93.33%) in the YL population, higher frequency (83.33% and 90%) in the FR population, high frequency (56.67% and 66.67%) in the YR population, and the lowest frequency (13.33% and 20 %) in the LR population. Therefore, SCAR01560bpand SCAR02443bpcan be used as special molecular markers for the population identification between LR and other three populations.

    Trachidermus fasciatus; Genetic diversity; RAPD; SCAR marker

    Q959.483; Q343

    A

    0254-5853-(2012)02-0203-08

    10.3724/SP.J.1141.2012.02203

    2011-11-03;接受日期:2012-02-14

    上海市教育委員會重點(diǎn)學(xué)科(海洋生物學(xué))建設(shè)項(xiàng)目(J50701);上海市科技興農(nóng)項(xiàng)目(滬農(nóng)科推字2006-3-5);上海市水生生物學(xué)重點(diǎn)學(xué)科建設(shè)項(xiàng)目(S30701);上海市重點(diǎn)學(xué)科建設(shè)項(xiàng)目(Y1101)

    ?通信作者(Corresponding authors),E-mail: zzliu@shou.edu.cn; ldpan@shou.edu.cn

    #2007生物科學(xué)(海洋生物方向本科生)

    曾珍(1987-),女,江西撫州人,碩士研究生。專業(yè)方向?yàn)轸~類分子進(jìn)化與免疫遺傳學(xué)。E-mail:zengzhen2010@126.com

    猜你喜歡
    種群特異性遺傳
    邢氏水蕨成功繁衍并建立種群 等
    非遺傳承
    山西省發(fā)現(xiàn)刺五加種群分布
    還有什么會遺傳?
    還有什么會遺傳
    還有什么會遺傳?
    精確制導(dǎo) 特異性溶栓
    BOPIM-dma作為BSA Site Ⅰ特異性探針的研究及其應(yīng)用
    重復(fù)周圍磁刺激治療慢性非特異性下腰痛的臨床效果
    兒童非特異性ST-T改變
    亚洲欧美精品专区久久| 啦啦啦啦在线视频资源| 秋霞伦理黄片| 日本黄大片高清| 两个人的视频大全免费| 日韩电影二区| 国精品久久久久久国模美| 亚洲精品色激情综合| 丰满乱子伦码专区| 国产美女午夜福利| 欧美潮喷喷水| 大话2 男鬼变身卡| 亚洲精品久久久久久婷婷小说| 久久鲁丝午夜福利片| 免费黄色在线免费观看| 国产人妻一区二区三区在| 青春草国产在线视频| 极品教师在线视频| 久久精品国产鲁丝片午夜精品| 99热这里只有是精品50| 少妇的逼水好多| 精品久久久久久久人妻蜜臀av| 人人妻人人澡人人爽人人夜夜 | 久久6这里有精品| 色视频www国产| 国产色婷婷99| 高清午夜精品一区二区三区| 亚洲av福利一区| 亚洲va在线va天堂va国产| 精品久久久久久成人av| 国产精品1区2区在线观看.| 亚洲精品成人久久久久久| 中国国产av一级| 国产黄色视频一区二区在线观看| 国产又色又爽无遮挡免| 91aial.com中文字幕在线观看| 男女国产视频网站| 亚洲成人久久爱视频| 欧美97在线视频| 精品国产露脸久久av麻豆 | 在线免费观看不下载黄p国产| 精品一区二区三卡| 麻豆久久精品国产亚洲av| 国产精品av视频在线免费观看| 97人妻精品一区二区三区麻豆| 午夜精品国产一区二区电影 | 亚洲经典国产精华液单| 精华霜和精华液先用哪个| 国产综合精华液| 久久久色成人| 成人美女网站在线观看视频| 欧美区成人在线视频| 日本wwww免费看| 国精品久久久久久国模美| 亚洲国产成人一精品久久久| 国产男女超爽视频在线观看| 亚洲av男天堂| 久久久色成人| 国产在线男女| 日日啪夜夜撸| 91久久精品国产一区二区三区| 熟女人妻精品中文字幕| 色5月婷婷丁香| 亚洲真实伦在线观看| 精品久久国产蜜桃| 51国产日韩欧美| 又爽又黄a免费视频| 老师上课跳d突然被开到最大视频| 在线观看一区二区三区| 亚洲av成人精品一二三区| 午夜福利在线在线| 3wmmmm亚洲av在线观看| 少妇猛男粗大的猛烈进出视频 | 日韩电影二区| 亚洲乱码一区二区免费版| 中国国产av一级| 国产亚洲一区二区精品| 欧美不卡视频在线免费观看| 国产成人91sexporn| 性色avwww在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品日韩在线中文字幕| 综合色丁香网| 极品教师在线视频| 全区人妻精品视频| 国产亚洲91精品色在线| 久久精品夜色国产| 搞女人的毛片| 亚洲国产精品专区欧美| 2021少妇久久久久久久久久久| 男女那种视频在线观看| 麻豆久久精品国产亚洲av| 有码 亚洲区| 日韩 亚洲 欧美在线| 午夜视频国产福利| 日本色播在线视频| 最近手机中文字幕大全| 国产又色又爽无遮挡免| 亚洲精品久久久久久婷婷小说| 国产精品女同一区二区软件| 亚洲不卡免费看| 波野结衣二区三区在线| 精品久久久久久电影网| 国产国拍精品亚洲av在线观看| 男女边吃奶边做爰视频| 丝袜美腿在线中文| 大片免费播放器 马上看| 久久国内精品自在自线图片| 久久精品国产亚洲av天美| 久久这里有精品视频免费| 不卡视频在线观看欧美| 欧美日韩在线观看h| 日本爱情动作片www.在线观看| 特大巨黑吊av在线直播| 欧美日韩一区二区视频在线观看视频在线 | 免费看日本二区| 亚洲精品一区蜜桃| 中国美白少妇内射xxxbb| 日韩一本色道免费dvd| www.色视频.com| 麻豆国产97在线/欧美| 欧美成人一区二区免费高清观看| 成年免费大片在线观看| 观看免费一级毛片| 搞女人的毛片| 国产极品天堂在线| 日本猛色少妇xxxxx猛交久久| 校园人妻丝袜中文字幕| 国产中年淑女户外野战色| 国产91av在线免费观看| 黄色欧美视频在线观看| 欧美变态另类bdsm刘玥| 国产有黄有色有爽视频| 久久精品夜夜夜夜夜久久蜜豆| av又黄又爽大尺度在线免费看| 国产高清有码在线观看视频| 亚洲国产精品成人久久小说| 狠狠精品人妻久久久久久综合| 国产成人一区二区在线| 午夜福利在线观看免费完整高清在| 日韩欧美精品v在线| www.色视频.com| 又粗又硬又长又爽又黄的视频| 国产亚洲午夜精品一区二区久久 | 亚州av有码| 欧美精品国产亚洲| 亚洲自偷自拍三级| 国产高清有码在线观看视频| 老司机影院成人| 国产精品一二三区在线看| 99热这里只有是精品50| 黄色欧美视频在线观看| 日韩,欧美,国产一区二区三区| 国产高潮美女av| 国产 亚洲一区二区三区 | 国产av国产精品国产| 亚洲精品第二区| 精品一区二区三卡| 国产一级毛片七仙女欲春2| 亚洲怡红院男人天堂| 亚洲精品国产成人久久av| av女优亚洲男人天堂| 一级av片app| 亚洲美女搞黄在线观看| 尾随美女入室| 亚洲av中文av极速乱| 成人特级av手机在线观看| 能在线免费观看的黄片| 久久久久久久久久久丰满| 熟妇人妻久久中文字幕3abv| 精品不卡国产一区二区三区| 91午夜精品亚洲一区二区三区| 精品国产三级普通话版| 久久久精品94久久精品| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 18禁在线无遮挡免费观看视频| 午夜精品一区二区三区免费看| 国产黄色免费在线视频| 插阴视频在线观看视频| 欧美人与善性xxx| 色吧在线观看| 男人狂女人下面高潮的视频| 91久久精品电影网| 嫩草影院入口| 欧美人与善性xxx| 国产 亚洲一区二区三区 | 日韩亚洲欧美综合| 色综合色国产| 欧美另类一区| 亚洲天堂国产精品一区在线| 91aial.com中文字幕在线观看| 最后的刺客免费高清国语| 看黄色毛片网站| 男女下面进入的视频免费午夜| 精品一区二区三区视频在线| 熟妇人妻久久中文字幕3abv| 亚洲av成人av| 亚洲图色成人| 国产精品福利在线免费观看| 国产精品嫩草影院av在线观看| 国产亚洲精品av在线| 日韩伦理黄色片| 一本一本综合久久| 国产精品.久久久| 欧美精品国产亚洲| 国产乱人偷精品视频| 神马国产精品三级电影在线观看| 插阴视频在线观看视频| 国产探花极品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 在现免费观看毛片| 国产男女超爽视频在线观看| 日本爱情动作片www.在线观看| 欧美一级a爱片免费观看看| 久久国产乱子免费精品| 一本一本综合久久| 国产高清不卡午夜福利| 观看美女的网站| 十八禁国产超污无遮挡网站| av播播在线观看一区| 日本一二三区视频观看| 99re6热这里在线精品视频| 七月丁香在线播放| 国产91av在线免费观看| 亚洲成人精品中文字幕电影| 国产精品一区www在线观看| 免费大片黄手机在线观看| 国产国拍精品亚洲av在线观看| 欧美极品一区二区三区四区| 国产黄频视频在线观看| 国产成人aa在线观看| 人人妻人人澡人人爽人人夜夜 | 成人美女网站在线观看视频| 亚洲伊人久久精品综合| 亚洲精品456在线播放app| av.在线天堂| 久久99热这里只频精品6学生| 国产日韩欧美在线精品| 搡女人真爽免费视频火全软件| 精品人妻熟女av久视频| 国产精品久久久久久久电影| 毛片一级片免费看久久久久| 国产成人freesex在线| 中文精品一卡2卡3卡4更新| 91精品国产九色| 纵有疾风起免费观看全集完整版 | 一级毛片aaaaaa免费看小| 日韩伦理黄色片| 午夜爱爱视频在线播放| 国产精品1区2区在线观看.| 中国国产av一级| 国产成人精品一,二区| 26uuu在线亚洲综合色| 日韩 亚洲 欧美在线| 亚洲成人av在线免费| 国产精品精品国产色婷婷| 日本一本二区三区精品| 亚洲国产av新网站| 蜜桃亚洲精品一区二区三区| 亚洲精品日韩在线中文字幕| 2021少妇久久久久久久久久久| 黄色配什么色好看| 噜噜噜噜噜久久久久久91| 一级爰片在线观看| 一级毛片 在线播放| 麻豆av噜噜一区二区三区| 日本-黄色视频高清免费观看| 一夜夜www| 久久精品熟女亚洲av麻豆精品 | 晚上一个人看的免费电影| 一区二区三区高清视频在线| 性插视频无遮挡在线免费观看| 日韩成人伦理影院| 国产亚洲午夜精品一区二区久久 | 日本-黄色视频高清免费观看| 嫩草影院新地址| 成人欧美大片| 亚洲18禁久久av| 亚洲av在线观看美女高潮| 男女啪啪激烈高潮av片| 亚洲最大成人手机在线| 国产精品1区2区在线观看.| 偷拍熟女少妇极品色| 成年女人在线观看亚洲视频 | 亚洲精品国产av成人精品| 亚洲在久久综合| 日本猛色少妇xxxxx猛交久久| 午夜久久久久精精品| 国产亚洲5aaaaa淫片| 国精品久久久久久国模美| 啦啦啦韩国在线观看视频| 久久国产乱子免费精品| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 搡老乐熟女国产| 午夜精品一区二区三区免费看| 婷婷六月久久综合丁香| 欧美成人午夜免费资源| 黄片wwwwww| 亚洲国产色片| 亚洲成人av在线免费| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 久久综合国产亚洲精品| 麻豆成人av视频| 日本猛色少妇xxxxx猛交久久| av在线老鸭窝| 日日撸夜夜添| 国产大屁股一区二区在线视频| 亚洲欧洲国产日韩| 亚洲美女搞黄在线观看| 亚洲美女视频黄频| 精品99又大又爽又粗少妇毛片| 亚洲自拍偷在线| 少妇熟女欧美另类| 亚洲精品亚洲一区二区| 亚洲精品第二区| 91久久精品电影网| 人妻少妇偷人精品九色| 亚洲精品aⅴ在线观看| 最后的刺客免费高清国语| 网址你懂的国产日韩在线| 免费大片18禁| 欧美激情国产日韩精品一区| 久久综合国产亚洲精品| 免费av观看视频| 床上黄色一级片| 国产成人精品一,二区| 高清日韩中文字幕在线| 国产乱人偷精品视频| 最新中文字幕久久久久| 久久久久久久国产电影| 国产真实伦视频高清在线观看| 非洲黑人性xxxx精品又粗又长| 麻豆久久精品国产亚洲av| 禁无遮挡网站| 三级国产精品欧美在线观看| 亚洲国产色片| 国产乱人视频| 亚洲精品第二区| 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 亚洲精品aⅴ在线观看| 国产精品熟女久久久久浪| 精品国产一区二区三区久久久樱花 | 老司机影院毛片| 少妇猛男粗大的猛烈进出视频 | 精品国产三级普通话版| 直男gayav资源| 日韩中字成人| 欧美 日韩 精品 国产| 亚洲国产精品sss在线观看| 久久精品国产亚洲网站| 国产亚洲5aaaaa淫片| 久久精品国产亚洲网站| 亚洲自拍偷在线| 欧美3d第一页| 日韩欧美精品v在线| 成年人午夜在线观看视频 | 国产午夜福利久久久久久| 一级毛片电影观看| 国产亚洲av片在线观看秒播厂 | 在线观看免费高清a一片| 日韩欧美一区视频在线观看 | 在线观看免费高清a一片| 日韩一区二区视频免费看| 欧美日韩在线观看h| 亚洲四区av| 亚洲在久久综合| 日本猛色少妇xxxxx猛交久久| 亚洲三级黄色毛片| 欧美日韩一区二区视频在线观看视频在线 | 在线播放无遮挡| 搡女人真爽免费视频火全软件| 国产精品熟女久久久久浪| 中文资源天堂在线| 国产男人的电影天堂91| 日韩视频在线欧美| 亚洲欧洲日产国产| 成人综合一区亚洲| 男女国产视频网站| 亚洲,欧美,日韩| 国产精品女同一区二区软件| 看免费成人av毛片| av国产久精品久网站免费入址| 国产成人精品福利久久| 成人无遮挡网站| 国产高清不卡午夜福利| 男人和女人高潮做爰伦理| 永久网站在线| 亚洲人与动物交配视频| 街头女战士在线观看网站| 色视频www国产| 亚洲国产精品sss在线观看| 如何舔出高潮| 亚洲欧洲国产日韩| 国产色爽女视频免费观看| 五月玫瑰六月丁香| 成人午夜精彩视频在线观看| 久久久精品欧美日韩精品| 91狼人影院| 国产综合精华液| 看免费成人av毛片| 午夜老司机福利剧场| 九九在线视频观看精品| 在现免费观看毛片| freevideosex欧美| 特大巨黑吊av在线直播| 亚洲av.av天堂| 国产成人精品一,二区| 2018国产大陆天天弄谢| 97超碰精品成人国产| 听说在线观看完整版免费高清| 中国美白少妇内射xxxbb| 搡老妇女老女人老熟妇| 夫妻性生交免费视频一级片| 高清欧美精品videossex| 国产免费一级a男人的天堂| 最近2019中文字幕mv第一页| 亚洲精品,欧美精品| 女人被狂操c到高潮| av国产免费在线观看| 午夜精品国产一区二区电影 | 国产精品三级大全| 高清视频免费观看一区二区 | 欧美xxⅹ黑人| 亚洲精品成人av观看孕妇| 91aial.com中文字幕在线观看| 非洲黑人性xxxx精品又粗又长| 日本黄大片高清| 国产 一区精品| 精华霜和精华液先用哪个| 中文乱码字字幕精品一区二区三区 | 一级毛片aaaaaa免费看小| 99久久精品热视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产欧美在线一区| 97超碰精品成人国产| 久久99热这里只频精品6学生| 成人亚洲精品一区在线观看 | 久久久久久久久久黄片| av专区在线播放| 欧美性猛交╳xxx乱大交人| 看免费成人av毛片| 日本一二三区视频观看| 婷婷色麻豆天堂久久| 嘟嘟电影网在线观看| 九九久久精品国产亚洲av麻豆| 亚洲天堂国产精品一区在线| 偷拍熟女少妇极品色| 欧美激情国产日韩精品一区| 成人无遮挡网站| 两个人的视频大全免费| 欧美3d第一页| 久久久久久伊人网av| 成年女人在线观看亚洲视频 | 国产一级毛片在线| 亚洲国产精品专区欧美| 国产欧美日韩精品一区二区| 欧美一区二区亚洲| 人人妻人人澡人人爽人人夜夜 | 在线a可以看的网站| 99久国产av精品国产电影| 天天躁夜夜躁狠狠久久av| 成人二区视频| 亚洲一级一片aⅴ在线观看| 国产激情偷乱视频一区二区| 啦啦啦中文免费视频观看日本| 国产精品一区二区在线观看99 | 精品不卡国产一区二区三区| 熟女人妻精品中文字幕| 精品久久久久久久久亚洲| 伊人久久国产一区二区| 少妇的逼好多水| 熟妇人妻久久中文字幕3abv| 亚洲综合精品二区| 在线观看av片永久免费下载| 成人亚洲精品av一区二区| 日韩三级伦理在线观看| 中国美白少妇内射xxxbb| 一个人观看的视频www高清免费观看| 亚洲第一区二区三区不卡| 亚洲人与动物交配视频| 最新中文字幕久久久久| 国产探花在线观看一区二区| 国产男人的电影天堂91| 亚洲在线自拍视频| 国产精品福利在线免费观看| 最近最新中文字幕免费大全7| 极品教师在线视频| 久久久亚洲精品成人影院| 一本久久精品| 国产精品一区二区三区四区免费观看| 欧美性感艳星| av专区在线播放| 丰满乱子伦码专区| 久久鲁丝午夜福利片| 免费观看的影片在线观看| 国产亚洲5aaaaa淫片| 久久精品国产亚洲av天美| 日本猛色少妇xxxxx猛交久久| 一个人看视频在线观看www免费| 一本久久精品| 2021天堂中文幕一二区在线观| 精品久久久久久成人av| 国产精品三级大全| 成年版毛片免费区| 国产精品美女特级片免费视频播放器| 亚洲成人中文字幕在线播放| 日本一本二区三区精品| 综合色丁香网| 国产一级毛片七仙女欲春2| 天堂俺去俺来也www色官网 | 简卡轻食公司| 五月伊人婷婷丁香| 国产淫语在线视频| 最近中文字幕高清免费大全6| 成人亚洲精品一区在线观看 | 国产亚洲最大av| 搞女人的毛片| 高清欧美精品videossex| 只有这里有精品99| 久久久久久久国产电影| 又粗又硬又长又爽又黄的视频| 中文字幕久久专区| 中文精品一卡2卡3卡4更新| 免费黄频网站在线观看国产| 色吧在线观看| 国产乱人偷精品视频| 日韩强制内射视频| 精品不卡国产一区二区三区| 久久精品夜色国产| 亚洲国产精品专区欧美| 欧美日本视频| 日本一二三区视频观看| 欧美日韩在线观看h| 又黄又爽又刺激的免费视频.| 中国国产av一级| 国产男女超爽视频在线观看| 免费在线观看成人毛片| 久久久成人免费电影| 国产在线男女| 丰满少妇做爰视频| 国产精品福利在线免费观看| 成人美女网站在线观看视频| 日韩成人av中文字幕在线观看| 全区人妻精品视频| 日本黄色片子视频| 国产精品一区二区三区四区久久| 人妻制服诱惑在线中文字幕| 男人爽女人下面视频在线观看| 日韩伦理黄色片| 免费人成在线观看视频色| 一二三四中文在线观看免费高清| 国产精品综合久久久久久久免费| 国产久久久一区二区三区| 国产精品日韩av在线免费观看| 欧美不卡视频在线免费观看| 亚洲欧美一区二区三区国产| 久久这里只有精品中国| 人人妻人人看人人澡| 国产精品久久视频播放| 午夜福利在线在线| 午夜福利成人在线免费观看| 男女下面进入的视频免费午夜| 激情 狠狠 欧美| 精品国内亚洲2022精品成人| 色哟哟·www| 亚洲丝袜综合中文字幕| 久久久久久国产a免费观看| 一级毛片aaaaaa免费看小| 一边亲一边摸免费视频| 欧美97在线视频| 久久久久久国产a免费观看| 国产又色又爽无遮挡免| 国产黄色小视频在线观看| 九草在线视频观看| 禁无遮挡网站| 亚洲av中文av极速乱| a级毛色黄片| 亚洲欧美精品自产自拍| 亚洲欧美中文字幕日韩二区| 人人妻人人澡欧美一区二区| 一级黄片播放器| 内射极品少妇av片p| 国产精品久久久久久av不卡| 麻豆成人av视频| 国产精品久久视频播放| 国产白丝娇喘喷水9色精品| 大片免费播放器 马上看| 1000部很黄的大片| 日韩欧美精品免费久久| 国产免费福利视频在线观看| 久久6这里有精品| 高清视频免费观看一区二区 | 久久久久九九精品影院| 国产高清三级在线| 国产精品日韩av在线免费观看| 欧美高清成人免费视频www| 亚洲欧洲国产日韩| 久久精品国产亚洲av天美| 美女黄网站色视频| 久久久欧美国产精品| 一级毛片电影观看| 日韩欧美一区视频在线观看 | 国产黄色免费在线视频| 91久久精品国产一区二区成人| 日韩,欧美,国产一区二区三区| 一个人观看的视频www高清免费观看| 亚洲av中文字字幕乱码综合| 欧美性猛交╳xxx乱大交人| 亚洲av成人精品一区久久| 日韩av免费高清视频| 精品99又大又爽又粗少妇毛片|