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

    蝦夷扇貝養(yǎng)殖群體及其子代的遺傳多樣性分析

    2012-01-12 12:03:40韓斐斐張繼彪張國范閆喜武
    海洋科學(xué) 2012年9期
    關(guān)鍵詞:扇貝微衛(wèi)星雜合

    韓斐斐, 張繼彪, 李 莉, 張國范, 閆喜武

    (1. 中國科學(xué)院 海洋研究所, 山東 青島266071; 2. 大連海洋大學(xué), 遼寧 大連116023)

    蝦夷扇貝養(yǎng)殖群體及其子代的遺傳多樣性分析

    韓斐斐1,2, 張繼彪1, 李 莉1, 張國范1, 閆喜武2

    (1. 中國科學(xué)院 海洋研究所, 山東 青島266071; 2. 大連海洋大學(xué), 遼寧 大連116023)

    蝦夷扇貝(Patinopecten(Mizuhopecten)yessoensis)于1982年從日本引進(jìn), 在中國已經(jīng)成為一個(gè)很重要的產(chǎn)業(yè)。本研究利用微衛(wèi)星標(biāo)記, 對(duì)獐子島海域2個(gè)養(yǎng)殖群體蝦夷扇貝(大耗島、褡褳島), 1個(gè)日本野生群體遺傳多樣性進(jìn)行了評(píng)估, 并對(duì) 2個(gè)養(yǎng)殖群體的群體內(nèi)繁育和群體間雜交后代進(jìn)行了遺傳和生長檢測(cè)。結(jié)果表明:養(yǎng)殖群體與日本野生群體相比, 其遺傳多樣性沒有顯著降低, 養(yǎng)殖群體遺傳多樣性較高, 養(yǎng)殖狀況暫時(shí)良好。TL群體和DH群體遺傳距離很近, DH、TL群體和RB群體的遺傳距離遠(yuǎn)大于TL群體和DH群體間的遺傳距離。說明2個(gè)養(yǎng)殖蝦夷扇貝群體遺傳分化較小。分別對(duì)4個(gè)子代群體的殼長、殼高在幼苗培育第1、5、10、15、20、100、190天進(jìn)行測(cè)量, 結(jié)果表明生長性狀差異不顯著。4個(gè)子代群體的遺傳多樣性差異也不顯著, 說明養(yǎng)殖群體的遺傳背景不清楚, 種質(zhì)資源混交較嚴(yán)重。本研究結(jié)果反映了我國蝦夷扇貝種質(zhì)資源評(píng)估狀況, 并為蝦夷扇貝的健康養(yǎng)殖提供基礎(chǔ)數(shù)據(jù)。

    蝦夷扇貝(Patinopecten(Mizuhopecten)yessoensis); 微衛(wèi)星; 遺傳多樣性; 雜交; 生長檢測(cè)

    蝦夷扇貝(Patinopecten(Mizuhopecten)yessoensis)原產(chǎn)于日本海周邊海域如日本、朝鮮、俄羅斯的海參崴等, 為寒溫帶貝類。1982年從日本引進(jìn), 并人工育苗獲得成功。我國適宜養(yǎng)殖的海區(qū)僅限于黃海北部冷水團(tuán)能夠覆蓋到的海域, 如遼寧的長山群島,大連沿岸海域, 山東廟島列島的北部諸島及榮成桑溝灣等。其中以大連長??h養(yǎng)殖規(guī)模最大。據(jù)漁業(yè)年鑒中扇貝產(chǎn)量和蝦夷扇貝所占比例估算, 目前蝦夷扇貝年產(chǎn)量已達(dá)30萬t, 年產(chǎn)值達(dá)數(shù)十億元, 成為中國淺海養(yǎng)殖的支柱性產(chǎn)業(yè), 也是國際上貝類引種養(yǎng)殖最成功的例證之一。

    蝦夷扇貝主要有筏式養(yǎng)殖和底播增殖, 所用苗種幾乎全部依賴育苗室的人工培育。自引進(jìn)后, 蝦夷扇貝以其獨(dú)有的性狀優(yōu)勢(shì)為市場(chǎng)所歡迎, 使其養(yǎng)殖規(guī)模不斷擴(kuò)大, 現(xiàn)已成為北方海域海水貝類的主要養(yǎng)殖對(duì)象之一[1]。自 1990年以來, 蝦夷扇貝筏式養(yǎng)殖生產(chǎn)規(guī)模不斷擴(kuò)大, 2001~2003年, 蝦夷扇貝增養(yǎng)殖規(guī)模進(jìn)入了鼎盛時(shí)期[2]。但近幾年養(yǎng)殖蝦夷扇貝死亡率急劇增加, 在有些海區(qū), 據(jù)估算其死亡率已達(dá)70%~80%, 與此同時(shí), 育苗穩(wěn)定性和成功率急劇下降。

    自蝦夷扇貝引進(jìn)中國后, 主要采取人工培育方式進(jìn)行蝦夷扇貝苗種生產(chǎn), 近年來在黃海北部到渤海海峽雖然也發(fā)現(xiàn)一些天然苗種, 但數(shù)量有限且產(chǎn)量不穩(wěn)定, 養(yǎng)殖生產(chǎn)所用苗種還是主要依靠人工苗種[3]。隨著養(yǎng)殖持續(xù)時(shí)間的延展, 養(yǎng)殖業(yè)規(guī)模的不斷擴(kuò)大, 養(yǎng)殖所面臨的問題也不斷增加, 如育苗難度越來越大, 成功率越來越低, 這不但增加了生產(chǎn)成本, 也影響到產(chǎn)業(yè)的健康發(fā)展。目前中國蝦夷扇貝人工苗種生產(chǎn)所用親本的來源主要是從自然海區(qū)挑選,挑選的標(biāo)準(zhǔn)主要是由有經(jīng)驗(yàn)的技術(shù)人員檢查待選親本的規(guī)格、殼型、肥滿度、性腺等表型情況[4]。另外,為了增加育苗的成功率, 往往有目的隨機(jī)挑選兩個(gè)不同地點(diǎn)的親貝, 通過群體間雜交的方式進(jìn)行苗種生產(chǎn)[5]。但是由于對(duì)被挑選群體的遺傳背景知之甚少,故這種隨機(jī)挑選有很大的盲目性。

    目前, 對(duì)于群體的遺傳多樣性檢測(cè)主要利用微衛(wèi)星標(biāo)記。微衛(wèi)星的開發(fā)和利用在蝦夷扇貝中的應(yīng)用日趨增多。常亞青等[6]用 8對(duì)微衛(wèi)星標(biāo)記對(duì)不同區(qū)域的 5個(gè)蝦夷扇貝群體進(jìn)行群體遺傳學(xué)分析, 從5個(gè)群體的雜合子狀態(tài)表明遺傳變異主要存在群體內(nèi)個(gè)體之間。趙瑩瑩等[7]利用磁珠富集法篩選蝦夷扇貝微衛(wèi)星分子標(biāo)記, 并在俄羅斯和大連 2個(gè)蝦夷扇貝群體中進(jìn)行遺傳多樣性分析。李春艷等[8]采用新型生物素-磁珠吸附微衛(wèi)星與同位素雜交相結(jié)合的方法篩選微衛(wèi)星引物, 并對(duì)大連地區(qū) 3個(gè)蝦夷扇貝養(yǎng)殖群體的遺傳結(jié)構(gòu)進(jìn)行分析。楊璞等[9]用微衛(wèi)星引物的通用性對(duì)櫛孔扇貝(♀)×蝦夷扇貝(♂)進(jìn)行雜種鑒定。

    本研究利用微衛(wèi)星分子標(biāo)記, 對(duì)某育苗企業(yè)用于苗種生產(chǎn)的 2個(gè)養(yǎng)殖群體進(jìn)行雜交和自繁, 以一個(gè)日本群體作為參照, 研究了 2個(gè)親本群體和子代的遺傳多樣性。并且跟蹤檢測(cè)了雜交和自繁后代的生長情況, 以期從遺傳管理的角度, 對(duì)蝦夷扇貝的苗種生產(chǎn)給予相應(yīng)的指導(dǎo)和借鑒。

    1 材料與方法

    1.1 材料

    本實(shí)驗(yàn)用于進(jìn)行遺傳分析的7個(gè)蝦夷扇貝群體,分別是在大連獐子島的特定區(qū)域養(yǎng)殖的蝦夷扇貝,大耗島貝(簡稱DH)、褡褳島貝(簡稱TL)、從日本引種的日本貝(簡稱 RB)和大耗島貝與褡褳島貝雙列雜交子代:DH自繁 (WDH), TL♀×DH♂(WDT),DH♀×TL♂(WTD), TL 自繁(WTL)。其中, WDH、WDT、WTD、WTL四個(gè)子代用于遺傳分析的是稚貝(生長到190天測(cè)量生長指標(biāo)后取樣)。

    其中, 大耗島貝和褡褳島貝分別進(jìn)行群體內(nèi)自繁和群體間雜交繁育。從每個(gè)群體中選擇雌雄比例為3 :2的性腺發(fā)育完全成熟的個(gè)體。陰干1 h后將每一個(gè)扇貝單獨(dú)放入一個(gè) 100L的塑料桶中, 加滿14℃的新鮮海水充氣催產(chǎn)。配子期開始計(jì)時(shí), 2 h后開始逐一鏡檢排放的配子防止雌雄同體貝影響實(shí)驗(yàn)結(jié)果。當(dāng)每桶卵排放密度達(dá)到50個(gè)/mL時(shí)將親貝取出(保證實(shí)驗(yàn)用量), 分別等量混合每個(gè)群體的精和卵, 再把兩個(gè)群體的精和卵分別等份分開, 每一群體的卵分別與群體內(nèi)或群體間雄貝的精液授精, 得到四個(gè)交配組合:WDH、WDT、WTD、WTL。群體隨機(jī)各取48個(gè)樣本。

    1.2 基因組DNA的提取

    取蝦夷扇貝肌肉, 參照Sambrook等[10]采用常規(guī)的酚/氯仿抽提的方法提取DNA。

    1.3 微衛(wèi)星引物來源

    使用的10對(duì)蝦夷扇貝微衛(wèi)星引物的序列來源于已發(fā)表的引物序列[8,11-14], 詳細(xì)信息見表1。引物序列由上海生工生物工程有限公司合成。

    表1 蝦夷扇貝微衛(wèi)星基因座引物序列Tab. 1 Microsatellite primer sequence of Yesso scallop

    1.4 生長數(shù)據(jù)測(cè)量

    取蝦夷扇貝兩個(gè)養(yǎng)殖群體(大耗島 、褡褳島)的雙列雜交子代WDH、WTL、WDT、WTD個(gè)體各48個(gè), 在蝦夷扇貝幼蟲培育期的第 1、5、10、15、20天和稚貝期第100、190天, 分別使用顯微鏡和游標(biāo)卡尺測(cè)定其殼長、殼高。其中, 殼長指貝殼前后緣的最大距離, 殼高指殼頂?shù)礁咕壍淖畲缶嚯x[15]。獲得的數(shù)據(jù), 用Excel和SPSS13.0軟件分析和處理, 差異性顯著分析采用最小顯著差法(Least-Significant difference, LSD)。

    1.5 PCR反應(yīng)和數(shù)據(jù)分析

    PCR 反應(yīng)體系(15μL):20 ng的基因組 DNA、10×buffer緩沖液(20 mmol/L Mg2+)、0.15 μL Tag 酶、dNTP各0.2 mmol/L、引物各0.3 μmol/L。PCR反應(yīng)條件:95℃變性5 min; 94℃30 s, 退火30 s, 72℃30 s,30個(gè)循環(huán); 72℃延伸10 min, 4℃保存。各對(duì)引物的退火溫度見表1。用12%的非變性聚丙烯酰胺(丙烯酰胺和甲叉雙丙烯酰胺體積比為29 :1)凝膠電泳分離PCR反應(yīng)產(chǎn)物, 電泳液為1×TBE緩沖液, 電壓300 V,電泳 1~2 h(北京六一儀器廠 DYY-Ⅱ型電泳儀,DYCZ-30型電泳槽), Gel-red染色, 利用UVP凝膠成像儀成像。

    統(tǒng)計(jì)條帶以UVP公司Doc-It Ls Image Analysisi Software軟件為主, 以手工矯正統(tǒng)計(jì)輔助。利用POPGENEGENE (VERSION 1.31) 軟件統(tǒng)計(jì)微衛(wèi)星基因座的等位基因數(shù) (Observed number of alleles,Na)、有效等位基因數(shù) (Effective number of alleles,Nea)、觀測(cè)雜合度 (Observed heterozygosity,Ho)、期望雜合度 (Expected heterozygosity,He)、遺傳相似系數(shù) (Genetic similarity index,I)、群體間遺傳距離(Genetic Distance,Ds), 并進(jìn)行多態(tài)性信息含量(Polymorphism Information Content, PIC) 分析。

    參照 Botstein等[16]的方法計(jì)算多態(tài)性信息含量(Polymorphism Information Content, PIC):

    式中,Pi、Pj分別為群體中第i和第j個(gè)等位基因頻率,n為等位基因數(shù)。

    親本數(shù)目是影響蝦夷扇貝養(yǎng)殖群體遺傳多樣性的重要指標(biāo)。根據(jù)雜合度的下降來估計(jì)有效群體大小 (effective population size,Ne)[17], 公式為:

    公式中,Ne是估計(jì)有效群體個(gè)數(shù);t是世代數(shù);Ht是t代的平均期望雜合度;Ho是親代的平均期望雜合度。在本實(shí)驗(yàn)中,t= 1,Ht是F1群體的雜合度,Ho是親本群體的平均雜合度。

    近交系數(shù)(inbreeding coefficient,F)是指根據(jù)近親交配的世代數(shù), 基因的純化程度。近交率 (ΔF) 公式為:

    第t代的近交系數(shù) (Ft)[18]公式為:

    式中,Ne是估計(jì)有效群體個(gè)數(shù);t是世代數(shù)。

    2 結(jié)果

    2.1 生長數(shù)據(jù)分析

    蝦夷扇貝兩個(gè)養(yǎng)殖群體(大耗島、褡褳島)的雙列雜交子代WDH、WTL、WDT、WTD的殼長、殼高的幼蟲期第 1、5、10、15、20天和稚貝期第 100、190天測(cè)量結(jié)果見表2, 結(jié)果表明:4個(gè)實(shí)驗(yàn)組中, 在幼蟲期1、5日齡時(shí)雜交組WDT的殼長、殼高生長快于自繁組WDH和WTL, 但差異不顯著(P>0.05),雜交組WTD殼長、殼高生長慢于兩自繁組; 在幼蟲期10日齡時(shí), 自繁組 WTL的殼長生長快于其他 3個(gè)實(shí)驗(yàn)組, 但差異不顯著(P>0.05), WTD殼高生長最快, 與自繁組 WTL差異不顯著(P>0.05), 與其他兩個(gè)實(shí)驗(yàn)組差異顯著(P<0.05), 但不明顯。幼蟲期 15和 20日齡時(shí), 雜交組 WDT的殼長、殼高生長快于其他3個(gè)實(shí)驗(yàn)組, 15日齡時(shí)WDT的生長與其他 3組差異顯著(P<0.05), 20日齡時(shí) WDT與WDH差異不顯著(P>0.05)與其他兩個(gè)實(shí)驗(yàn)組差異顯著(P<0.05), 并不十分明顯。在稚貝期100日齡時(shí), 雜交組 WTD生長組快于自繁組差異不顯著(P>0.05), WDT生長最慢與自繁組WTL差異不顯著(P>0.05)與其他兩個(gè)實(shí)驗(yàn)組差異顯著(P<0.05); 190日齡時(shí)與 100日齡的生長趨勢(shì)相近, 即雜交組 WTD生長組快于自繁組差異顯著(P<0.05), WDT生長最慢與自繁組 WTL差異不顯著(P>0.05)與其他兩個(gè)實(shí)驗(yàn)組差異顯著(P<0.05); 稚貝期估算的兩群體并無顯著差異。

    2.2 PCR擴(kuò)增反應(yīng)結(jié)果

    實(shí)驗(yàn)合成了10對(duì)引物, 用每個(gè)群體中的8個(gè)樣本篩選引物, 引物均能顯示多態(tài)性。利用微衛(wèi)星引物在蝦夷扇貝群體 DH、TL、RB和子代群體 WDH、WTL、WTD、WDT進(jìn)行 PCR擴(kuò)增, 都能擴(kuò)增出特異性條帶(圖1)。其中, 除XY28的多態(tài)性稍低外, 其他位點(diǎn)均表現(xiàn)出較高的多態(tài)性。

    2.3 群體遺傳多樣性分析

    應(yīng)用 POPGENE(VERSION 1.31)軟件對(duì)蝦夷扇貝群體進(jìn)行群體遺傳學(xué)分析。10個(gè)微衛(wèi)星位點(diǎn)在蝦夷扇貝群體大耗島貝(DH)、褡褳島貝(TL)和日本貝(RB)中擴(kuò)增出的等位基因數(shù)、有效等位基因數(shù)、觀測(cè)雜合度、期望雜合度、雜交系數(shù)、多態(tài)性信息含量(表3)。在群體遺傳學(xué)中,Na、He和Ho等都是反映群體遺傳多樣性的度量, 其數(shù)值越大, 說明基因豐富度越高[19]。本研究中蝦夷扇貝養(yǎng)殖群體都表現(xiàn)出較高的遺傳多樣性水平。

    在蝦夷扇貝2個(gè)養(yǎng)殖群體DH、TL和野生群體RB中, TL群體的平均等位基因數(shù)最多, 為4.2613個(gè), DH群體最少, 為 3.8535個(gè), 與觀測(cè)到的等位基因數(shù)相比,每個(gè)位點(diǎn)的有效等位基因數(shù)較少。多態(tài)信息含量是衡量位點(diǎn)多樣性的較好指標(biāo)。在 10個(gè)多態(tài)位點(diǎn)中, 除XY28、XY35在DH、TL、RB群體表現(xiàn)為中度多態(tài), RB在 XY59中表現(xiàn)為中度多態(tài)外, 大部分位點(diǎn)多態(tài)信息含量較高。觀測(cè)雜合度為0.4467~0.6047, 期望雜合度為0.6436~0.6983, 多態(tài)性信息含量為0.5599~0.6457,這些位點(diǎn)在蝦夷扇貝中都具有較高的多態(tài)信息含量,可以應(yīng)用于蝦夷扇貝群體遺傳的研究。

    表2 蝦夷扇貝殼長、殼高的平均值和標(biāo)準(zhǔn)差及中親雜種優(yōu)勢(shì)Tab. 2 The mean value of shell length, height, standard devitation and heterosis in Yesso scallop

    圖1 XY4引物在蝦夷扇貝TL群體1-24號(hào)個(gè)體中的擴(kuò)增情況Fig. 1 The polymorphism at locus XY4 in Yesso scallop (TL population 1-24 individuals)

    表3 2個(gè)蝦夷扇貝養(yǎng)殖群體及其子代、1個(gè)野生群體和4個(gè)F1子代群體的遺傳多態(tài)性(均值)Tab. 3 The polymorphism information between two cultured populations, one natural population and four F1 populations of Yesso scallop (mean value)

    群體 DH 和 TL的等位基因數(shù)分別是 6.1~6.2個(gè), 略大于群體WDH 、WTL、WTD、WDT等位基因數(shù), 群體 DH 和 TL 與 WDH 、WTL、WTD、WDT的有效等位基因數(shù)相差不大。群體DH 和TL的觀測(cè)雜合度和期望雜合度略大于子代群體WDH 、WTL、WTD、WDT, 說明群體內(nèi)交配出現(xiàn)輕微的雜合度缺失。4個(gè)子代群體在 10對(duì)微衛(wèi)星引物中, 除有一對(duì)引物的多態(tài)性信息含量(PIC)值在0.19~0.31之間外,其余引物在各群體中的PIC值為0.36~0.85, 說明蝦夷扇貝目前的養(yǎng)殖狀況良好。

    2.4 遺傳分化的計(jì)算

    通過比較兩兩群體Fst值(表4)可以看出, DH群體和TL群體間的遺傳分化指數(shù)低于0.05, 說明兩個(gè)群體間的遺傳分化程度較低, 而TL群體和RB群體、DH群體和RB群體間的遺傳分化程度較高。

    表4 蝦夷扇貝DH、TL、RB群體的遺傳分化系數(shù)FstTab. 4 Pairwise estimates of Fst within three populations of Yesso scallop

    通過比較兩兩群體Fst值(表5)可以看出, 4個(gè)群體間的遺傳分化指數(shù)均低于0.05, 說明4個(gè)群體間的遺傳分化程度較低。

    表5 蝦夷扇貝F1代4個(gè)群體的Fst比較遺傳分化系數(shù)FstTab. 5 Pairwise estimates of Fst within four F1 populations of Yesso scallop

    2.5 Hardy-Weinberg平衡檢驗(yàn)

    在理想群體中, 群體的基因頻率和基因型頻率是穩(wěn)定的。應(yīng)用POPGENE(VERSION 1.31)數(shù)據(jù)分析軟件, 基于馬可夫鏈模型(Markov chain method)對(duì)各群體每個(gè)位點(diǎn)的Hardy-Weinberg P值進(jìn)行無偏估計(jì),4個(gè) F1子代群體的檢驗(yàn)結(jié)果表現(xiàn)出不同程度的平衡偏離, 其中9個(gè)位點(diǎn)沒有偏離平衡, 31(77.5%)位點(diǎn)偏離了平衡, 說明這些位點(diǎn)在各群體中基因頻率和基因型頻率穩(wěn)定性較差。

    2.6 群體間遺傳距離

    根據(jù)Nei[20]分別計(jì)算了DH、TL、RB 3個(gè)群體和子代4個(gè)群體的遺傳距離和遺傳相似性指數(shù)。TL群體和 DH群體遺傳距離最近, 遺傳相似性指數(shù)最高; 和RB群體的遺傳距離次之, DH群體和RB群體的遺傳距離最遠(yuǎn), 相似性指數(shù)最低(表6)。這說明從日本當(dāng)年引種的日本貝(RB)與大耗島貝(DH)、褡褳島貝(TL)間的親緣關(guān)系較遠(yuǎn)。WTL群體和 WTD群體遺傳距離較近, 遺傳相似性指數(shù)較高; WDH群體和 WDT群體遺傳距離較近, 遺傳相似性指數(shù)較高(表7)。

    表6 蝦夷扇貝 3個(gè)群體的 Nei遺傳距離(表左下角)及相似性指數(shù)(表右上角)Tab. 6 Pairwise estimates of Nei’s unbiased genetic distance (below diagonal) and similarity index(above diagonal) within three populations of Yesso scallop

    表7 蝦夷扇貝4個(gè)F1代群體的Nei遺傳距離(表左下角)及相似性指數(shù)(表右上角)Tab. 7 Pairwise estimates of Nei’s unbiased genetic distance (below diagonal) and similarity index(above diagonal) within four F1 populations of Yesso scallop

    2.7 有效群體大小和近交系數(shù)

    WDH、WTL、WTD、WDT實(shí)際使用的親本總數(shù)大于最大有效群體數(shù)和估計(jì)有效群體數(shù)(表8)。相對(duì)較低的Ne能導(dǎo)致近交率的增加, 這是近交衰退的原因之一。WDH、WTL、WTD、WDT的近交系數(shù)(Ft)分別為 0.0313、0.0263、0.0152、0.0217。

    3 討論

    微衛(wèi)星分子標(biāo)記不受生物發(fā)育時(shí)期和環(huán)境的影響, 為中性標(biāo)記。利用微衛(wèi)星分子標(biāo)記分析蝦夷扇貝親代與子代遺傳多樣性, 能夠更加客觀地反映生物之間的本質(zhì)差異, 從而能夠更加客觀地反映生物的種群結(jié)構(gòu)和進(jìn)化歷史。在群體遺傳學(xué)中,Na、He和Ho等都是反映群體遺傳多樣性的度量, 其數(shù)值越大,說明基因豐富度越高[19]。群體的遺傳多樣性越高, 其適應(yīng)環(huán)境變化的能力越強(qiáng)。與野生群體相比, 遺傳多樣性在很多情況下都會(huì)降低, 其中在一些海洋生物中都有報(bào)道[21-22]。本研究中養(yǎng)殖蝦夷扇貝群體的遺傳多樣性比野生群體沒有顯著降低, 可能是由于養(yǎng)殖群體的養(yǎng)殖環(huán)境適宜, 使自然條件下應(yīng)該被淘汰的個(gè)體保存下來, 保持了養(yǎng)殖群體的遺傳多樣性。本實(shí)驗(yàn)中所檢測(cè)的大部分位點(diǎn)偏離了 Hardy-Weinberg平衡, 說明蝦夷扇貝的種質(zhì)情況不容樂觀,混雜現(xiàn)象較為嚴(yán)重。這與常亞青[6]研究結(jié)果相一致。

    表8 4個(gè)子代群體的有效群體數(shù)(Ne)Tab. 8 The effective population size (Ne) within four F1 populations of Yesso scallop

    本研究中國內(nèi)養(yǎng)殖群體和日本群體的群體間遺傳距離大于 0.4, 這很可能是由于蝦夷扇貝作為一個(gè)引進(jìn)種, 有效群體較小, 加之在苗種生產(chǎn)過程中生產(chǎn)單位往往用較少的親貝尤其是較少的父本, 奠基者效應(yīng)使得養(yǎng)殖群體和日本野生群體發(fā)生較大的遺傳分化。常亞青等[6]研究了5個(gè)蝦夷扇貝的遺傳多樣性, 蝦夷扇貝群體內(nèi)遺傳距離大于0.5。

    在本研究中, 4個(gè)實(shí)驗(yàn)組子代在各日齡, 殼長、殼高生長雖然有一定程度的差異, 但兩養(yǎng)殖群體的雜交子代并未表現(xiàn)出明顯的差異。分析產(chǎn)生這種現(xiàn)象的原因是由于兩群體蝦夷扇貝生長性狀遺傳差異不顯著而引起的, 生產(chǎn)用親貝種質(zhì)混雜現(xiàn)象比較嚴(yán)重。張立冬等[23]研究了大長山群體、日本陸奧灣群體和獐子島群體 3個(gè)群體的遺傳多樣性, 發(fā)現(xiàn)陸奧灣群體和獐子島群體遺傳距離最近, 遺傳相似性指數(shù)最高; 陸奧灣群體和大長山群體的遺傳距離最遠(yuǎn),相似性指數(shù)最低, 并指出造成這種現(xiàn)象原因是由于獐子島在全人工育苗過程中多采用日本扇貝或利用日本扇貝與當(dāng)?shù)貎?yōu)良種貝進(jìn)行雜交, 故在遺傳多樣性水平上與日本蝦夷扇貝群體相差不大; 而大長山全人工苗種多采自當(dāng)?shù)? 累代繁殖及近交等因素影響了其遺傳多樣性水平, 造成了當(dāng)?shù)厝后w與其他兩個(gè)地區(qū)群體的差異。本研究所選用的兩養(yǎng)殖群體的親本可能是由于生產(chǎn)用親貝種質(zhì)混雜現(xiàn)象比較嚴(yán)重,從而在生長性狀并未獲得顯著的雜種優(yōu)勢(shì)。

    目前中國蝦夷扇貝人工苗種生產(chǎn)所用親本的來源主要是從自然海區(qū)挑選, 挑選的標(biāo)準(zhǔn)主要是由有經(jīng)驗(yàn)的技術(shù)人員檢查待選親本的規(guī)格、殼型、肥滿度、性腺等表型情況, 另外也要對(duì)既往用其作為親本的繁育成敗情況等做社會(huì)調(diào)查。但這些檢查和調(diào)查主要還是憑經(jīng)驗(yàn), 也有運(yùn)氣成分。通過本研究可以發(fā)現(xiàn), 蝦夷扇貝隨機(jī)雜交配易造成遺傳分化的降低,因此, 生產(chǎn)上這種隨機(jī)采用用兩個(gè)地理位置不同的群體雜交是盲目的, 造成蝦夷扇貝種質(zhì)衰退, 對(duì)蝦夷扇貝養(yǎng)殖產(chǎn)業(yè)造成沖擊。因此, 從分子水平來研究兩群體的自繁子代和雜交子代的遺傳多樣性, 從根本上解釋自繁群體和雜交群體的遺傳結(jié)構(gòu), 采用分子育種方法指導(dǎo)蝦夷扇貝的種質(zhì)改良是非常有意義的。

    鑒于蝦夷扇貝在養(yǎng)殖過程中很容易出現(xiàn)不同群體混雜的現(xiàn)象, 我們建議實(shí)行“種貝”和“菜貝”分區(qū)養(yǎng)殖的模式。所謂的“菜貝”去向是作為商品貝流向市場(chǎng)被消費(fèi)者購買和食用, 所謂的“種貝”是用于進(jìn)行苗種生產(chǎn)的貝類。用于“種貝”的親貝必須有系統(tǒng)的譜系分析和良好的保種措施。保種的過程中不但要確保各個(gè)種群之間具有嚴(yán)格的隔離, 而且在世代之間能夠確保遺傳結(jié)構(gòu)得到最大程度的保留。

    [1]李文姬, 薛真福. 持續(xù)發(fā)展蝦夷扇貝的健康養(yǎng)殖[J].水產(chǎn)科學(xué), 2005, 24(9):49-51.

    [2]張明明, 趙文. 我國蝦夷扇貝死亡原因的探討和控制對(duì)策[J]. 魚病診治, 2008, 2:65-74.

    [3]常亞青, 王慶志, 宋堅(jiān), 等. 黃海北部大連沿岸蝦夷扇貝天然苗采集技術(shù)研究[J]. 海洋水產(chǎn)研究, 2007,28(6):39-44.

    [4]王海濤, 王世黨, 姜啟平, 等. 蝦夷扇貝育苗綜合技術(shù)[J]. 中國水產(chǎn), 2008, 4:55-56.

    [5]楊愛國, 王清印, 劉志鴻, 等. 櫛孔扇貝與蝦夷扇貝雜交及子一代的遺傳性狀[J]. 海洋水產(chǎn)研究, 2004,25(5):1-5.

    [6]常亞青, 陳曉霞, 丁君, 等. 蝦夷扇貝(Patinopecten yessoensis)5 個(gè)群體的遺傳多樣性[J]. 生態(tài)學(xué)報(bào),2007, 27(3):1145-1152.

    [7]趙瑩瑩, 朱曉琛, 孫效文, 等. 蝦夷扇貝的多態(tài)性微衛(wèi)星座位[J]. 動(dòng)物學(xué)報(bào), 2006, 52(1):229-233.

    [8]李春艷, 丁君, 常亞青, 等. 蝦夷扇貝微衛(wèi)星標(biāo)記的分離及其養(yǎng)殖群體的遺傳結(jié)構(gòu)分析[J]. 中國水產(chǎn)科學(xué), 2009, 16(1):39-46.

    [9]楊璞, 楊愛國, 劉志鴻, 等. 櫛孔扇貝和蝦夷扇貝通用微衛(wèi)星引物的篩選及其在雜種鑒定中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué), 2008, 36(19):8287-8289.

    [10]Sambrook J, Fritsch E F, Maniatis T. Molecular Cloning:a Laboratory Manual, 2nd edn[M]. New York:Cold Spring Harbor Laboratory Press, 1989.

    [11]慧敏. 扇貝微衛(wèi)星標(biāo)記的篩選及其在物種鑒定中的應(yīng)用[D]. 青島:中國海洋大學(xué), 2008:28-29.

    [12]高祥剛, 曹潔, 劉瑩, 等. SSR 不對(duì)稱 PCR 法分析蝦夷扇貝遺傳多樣性[J]. 生物技術(shù)通報(bào), 2009, 2:118-123.

    [13]李云峰, 劉衛(wèi)東, 高祥剛, 等. 蝦夷扇貝外套膜和腎臟組織cDNA文庫構(gòu)建以及 EST的初步分析[J]. 中國水產(chǎn)科學(xué), 2010, 17(3):578-585.

    [14]趙海波. 蝦夷扇貝遺傳連鎖圖譜的構(gòu)建[D]. 青島:中國海洋大學(xué), 2009:34.

    [15]張存善, 常亞青, 曹學(xué)彬, 等. 蝦夷扇貝體形性狀對(duì)軟體重和閉殼肌重的影響效果分析[J]. 水產(chǎn)學(xué)報(bào),2009, 33(1):87-94.

    [16]Botstein D, White R L, Skolnick M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980,32(3):314-331.

    [17]Saavedra C, Guerra A. Allozyme heterozygosity, founder effect and fitness traits in a cultivated population of the European oyster, Ostrea edulis[J]. Aquaculture,1996, 139(3-4):203-224.

    [18]張國范, 鄭懷平. 海灣扇貝養(yǎng)殖遺傳學(xué)[M]. 北京:科學(xué)出版社, 2009:36-50.

    [19]陳蒙, 常亞青, 孫謙, 等. 蝦夷扇貝群體的遺傳結(jié)構(gòu)及微衛(wèi)星標(biāo)記與體尺、體重的相關(guān)性分析[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2009, 24(4):311-316.

    [20]Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals[J].Genetics, 1978, 89(3):583.

    [21]Hedgecock D, Sly F. Genetic drift and effective population sizes of hatchery-propagated stocks of the Pacific oysterCrassostrea gigas[J]. Aquaculture. 1990, 88:21-38.

    [22]Porta J, Porta J M, Martinez-Rodriguez G, et al. Genetic structure and genetic relatedness of a hatchery stock of Senegal sole (Solea senegalensis) inferred by microsatellites[J]. Aquaculture, 2006, 251:46-55.

    [23]張立冬, 王蕾, 高悅勉. 蝦夷扇貝遺傳結(jié)構(gòu)及微衛(wèi)星標(biāo)記與經(jīng)濟(jì)性狀相關(guān)分析[J]. 生物技術(shù)通報(bào), 2010,11:115-112.

    Genetic diversity in two-generation of cultured Japanese scallop (Patinopecten(Mizuhopecten)yessoensis) populations

    HAN Fei-fei1,2, ZHANG Ji-biao1, LI Li1, ZHANG Guo-fan1, YAN Xi-wu2
    (1. Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. Dalian Ocean University,Dalian 116023, China)

    Dec.,12,2011

    Patinopecten(Mizuhopecten)yessoensis; microsatellites; genetic diversity; hybridization; growth detection

    The Japanese scallop (Patinopecten(Mizuhopecten)yessoensis) was introduced into China in 1982. Since then, the culture of Japanese scallop has become an important industry in China. In this study microsatellite was used to investigate the genetic variations of two cultured populations [one from Dahao island (DH) and another from Talian island (TL)]and one wild Japanese scallop population (RB). The four F1progeny populations were the one hatchery population which from Dahao island (WDH), one population from Talian island (WTL), one outcrossed population (WTD population) from Dahao island (♀) and Talian island (♂), one outcrossed population(WDT population) from Dahao island (♂) and Talian island (♀). Forty-eight individuals from each population were analyzed. Amplification was performed in a PCR reaction and its products were separated by 12% non-denaturing polyacrylamide gel. The difference in genetic diversity among two cultural populations (DH and TL population) and one natural population (RB population) was not significant. The two cultural populations were in good condition temporarily. According to their genetic similarity (I) and genetic distance (D), we found that the genetic distance between DH and TL was the least. The phenotype characters (shell length and shell height) of four F1progeny populations scallop were measured at 1, 5, 10, 15, 20, 100, 190 days old. Non-significant differences were found between different characters in four F1progeny populations at different growth stages. The difference in genetic diversity among the four F1progeny populations was not significant, either. The results showed the genetic background of cultured populations was not clear and cultured populations were mixed with each other seriously. The results of this study reflect resources evaluation situation for Japanese scallop in our country and provide the basis of data for health breeding of Japanese scallop.

    Q751 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1000-3096(2012)09-0001-08

    2011-12-12;

    2012-04-22

    中國科學(xué)院海洋研究所獐子島漁業(yè)海洋生態(tài)養(yǎng)殖聯(lián)合實(shí)驗(yàn)室開放研究基金; 中國科學(xué)院重要方向項(xiàng)目群(KZCX2-EW-Q21)

    韓斐斐(1984-), 女, 碩士研究生, 從事水產(chǎn)養(yǎng)殖分子生物學(xué)研究, E-mail:hanfeifei0325310406@126.com; 張國范, 通信作者, 研究員, E-mail:guofanzhang2005@163.com

    致謝:感謝中國科學(xué)院海洋研究所李娟、王家豐、王威和大連海洋大學(xué)霍忠明、楊霏、王琳楠對(duì)實(shí)驗(yàn)的幫助。

    梁德海)

    猜你喜歡
    扇貝微衛(wèi)星雜合
    梅花鹿基因組微衛(wèi)星分布特征研究
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    扇貝的眼睛在哪里
    烤扇貝:貝殼反復(fù)利用,殼比肉還貴
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國“一箭雙星”成功將“遙感衛(wèi)星二十一號(hào)”與“天拓二號(hào)視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    扇貝吃了一頭大象
    雄激素可調(diào)節(jié)的腎臟近端腎小管上皮細(xì)胞靶向雜合啟動(dòng)子的優(yōu)化
    精品午夜福利在线看| 亚洲第一区二区三区不卡| 丝袜喷水一区| 欧美少妇被猛烈插入视频| 三级国产精品欧美在线观看| 国产精品一区二区在线观看99| 老女人水多毛片| 欧美高清性xxxxhd video| 91精品国产国语对白视频| 男女边吃奶边做爰视频| 亚洲av电影在线观看一区二区三区| 成人黄色视频免费在线看| 插逼视频在线观看| 国产又色又爽无遮挡免| 天天躁日日操中文字幕| 色5月婷婷丁香| 涩涩av久久男人的天堂| 久久精品夜色国产| 国产精品伦人一区二区| 亚洲综合精品二区| 日本av手机在线免费观看| 国内揄拍国产精品人妻在线| 女性生殖器流出的白浆| 亚洲不卡免费看| 丰满少妇做爰视频| 一级二级三级毛片免费看| 一边亲一边摸免费视频| 中文精品一卡2卡3卡4更新| 日本与韩国留学比较| 亚洲三级黄色毛片| av不卡在线播放| 亚洲色图综合在线观看| 亚州av有码| 久久99热这里只有精品18| 寂寞人妻少妇视频99o| 欧美激情极品国产一区二区三区 | 最近2019中文字幕mv第一页| 青青草视频在线视频观看| 色哟哟·www| 黄片wwwwww| 色网站视频免费| 99久久人妻综合| 久久久久久久大尺度免费视频| videossex国产| 男女下面进入的视频免费午夜| 久久久久久久久久人人人人人人| 韩国av在线不卡| 国产精品久久久久久av不卡| 婷婷色综合www| 亚洲欧美清纯卡通| 少妇裸体淫交视频免费看高清| 亚洲av成人精品一二三区| 午夜福利在线在线| av国产免费在线观看| 国产成人91sexporn| 高清日韩中文字幕在线| 亚洲精品aⅴ在线观看| 亚洲国产高清在线一区二区三| 插逼视频在线观看| 久久久亚洲精品成人影院| 黄片无遮挡物在线观看| 国产精品秋霞免费鲁丝片| 免费不卡的大黄色大毛片视频在线观看| 高清av免费在线| 久久人人爽av亚洲精品天堂 | av视频免费观看在线观看| 又黄又爽又刺激的免费视频.| 亚洲人成网站高清观看| 久久久久久久国产电影| 国产爽快片一区二区三区| 国产精品熟女久久久久浪| 一区二区三区四区激情视频| 少妇裸体淫交视频免费看高清| 黄色怎么调成土黄色| 18禁在线播放成人免费| 综合色丁香网| 啦啦啦视频在线资源免费观看| 久久国产亚洲av麻豆专区| 美女内射精品一级片tv| 97热精品久久久久久| 一本久久精品| 亚洲欧美中文字幕日韩二区| 国产 精品1| 色视频www国产| 欧美高清成人免费视频www| 国产一区二区三区综合在线观看 | 免费看光身美女| kizo精华| 亚洲中文av在线| 插逼视频在线观看| 亚洲成色77777| 国产在线视频一区二区| 日韩视频在线欧美| 久久精品国产亚洲av天美| 九九爱精品视频在线观看| 亚洲精品乱码久久久v下载方式| 欧美极品一区二区三区四区| 久久久久久久久久久丰满| 新久久久久国产一级毛片| 中文字幕久久专区| 赤兔流量卡办理| 免费大片黄手机在线观看| 夫妻性生交免费视频一级片| 91久久精品电影网| 久久99蜜桃精品久久| 欧美日韩精品成人综合77777| 精品视频人人做人人爽| 亚洲av中文av极速乱| 激情 狠狠 欧美| 色视频www国产| 欧美+日韩+精品| 国产真实伦视频高清在线观看| 天堂8中文在线网| 婷婷色av中文字幕| 99精国产麻豆久久婷婷| 国产亚洲91精品色在线| 久久99热6这里只有精品| 国产精品不卡视频一区二区| 少妇被粗大猛烈的视频| 日韩免费高清中文字幕av| 国产高清三级在线| 汤姆久久久久久久影院中文字幕| 亚洲精品久久午夜乱码| 久久精品久久久久久久性| 高清视频免费观看一区二区| 99久久精品热视频| 国产美女午夜福利| 少妇高潮的动态图| 日韩大片免费观看网站| 少妇 在线观看| 国产精品不卡视频一区二区| 最新中文字幕久久久久| 麻豆成人午夜福利视频| 久久精品国产亚洲av涩爱| 亚洲成人中文字幕在线播放| 日韩免费高清中文字幕av| 国产色婷婷99| 2022亚洲国产成人精品| 午夜福利视频精品| 一级av片app| 在线播放无遮挡| 国产一区二区三区av在线| 91久久精品国产一区二区成人| av线在线观看网站| 日韩av在线免费看完整版不卡| 噜噜噜噜噜久久久久久91| 欧美3d第一页| 黑丝袜美女国产一区| 一本色道久久久久久精品综合| 久久久久网色| 国产伦理片在线播放av一区| 七月丁香在线播放| 99热这里只有是精品在线观看| 久久鲁丝午夜福利片| 日本vs欧美在线观看视频 | 极品少妇高潮喷水抽搐| 蜜臀久久99精品久久宅男| 欧美zozozo另类| 日韩av不卡免费在线播放| 久久99精品国语久久久| 高清毛片免费看| 熟女av电影| 精品人妻偷拍中文字幕| 午夜激情久久久久久久| 国产精品久久久久成人av| 国产深夜福利视频在线观看| 国产成人精品福利久久| 欧美日韩在线观看h| 舔av片在线| 久久人人爽人人片av| 成人综合一区亚洲| www.av在线官网国产| 人妻夜夜爽99麻豆av| 观看av在线不卡| 18禁裸乳无遮挡动漫免费视频| 国精品久久久久久国模美| 久久精品久久久久久久性| 亚洲人成网站在线观看播放| 大话2 男鬼变身卡| av免费在线看不卡| 六月丁香七月| 国产精品蜜桃在线观看| 人妻一区二区av| 日韩av不卡免费在线播放| 少妇的逼好多水| 热99国产精品久久久久久7| 成人美女网站在线观看视频| 国产精品久久久久久久电影| 欧美精品国产亚洲| 国产精品久久久久久精品电影小说 | 秋霞在线观看毛片| 在线天堂最新版资源| 一级毛片电影观看| 九九爱精品视频在线观看| 中文天堂在线官网| 91久久精品国产一区二区成人| 欧美97在线视频| av播播在线观看一区| 亚洲欧美日韩卡通动漫| 久久久欧美国产精品| 亚洲熟女精品中文字幕| 中文字幕精品免费在线观看视频 | 嫩草影院新地址| 国产亚洲av片在线观看秒播厂| 亚洲精品亚洲一区二区| 日本欧美国产在线视频| h日本视频在线播放| 久久久国产一区二区| 欧美 日韩 精品 国产| 国产免费福利视频在线观看| 91精品伊人久久大香线蕉| 一级毛片aaaaaa免费看小| 精华霜和精华液先用哪个| 久久国内精品自在自线图片| 日本黄色日本黄色录像| 国产亚洲91精品色在线| av黄色大香蕉| 亚洲精品日韩在线中文字幕| 99久久精品国产国产毛片| 内射极品少妇av片p| 国产一区亚洲一区在线观看| 老熟女久久久| 欧美日韩在线观看h| 一区二区三区乱码不卡18| 国产又色又爽无遮挡免| 一级片'在线观看视频| 极品教师在线视频| 2022亚洲国产成人精品| 男人爽女人下面视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 日韩强制内射视频| 老师上课跳d突然被开到最大视频| 久久99精品国语久久久| 精品人妻熟女av久视频| 夜夜看夜夜爽夜夜摸| 一级毛片久久久久久久久女| 亚洲国产毛片av蜜桃av| 日韩不卡一区二区三区视频在线| 麻豆乱淫一区二区| 纵有疾风起免费观看全集完整版| 午夜老司机福利剧场| 国产成人午夜福利电影在线观看| 美女国产视频在线观看| 久久久欧美国产精品| 日韩中字成人| 水蜜桃什么品种好| 亚洲欧美一区二区三区黑人 | 国产亚洲一区二区精品| 91久久精品电影网| 在现免费观看毛片| 大陆偷拍与自拍| 丰满迷人的少妇在线观看| 国产精品麻豆人妻色哟哟久久| 国产真实伦视频高清在线观看| 国产老妇伦熟女老妇高清| 亚洲中文av在线| 国产精品一区二区在线不卡| 99久久精品一区二区三区| 日韩制服骚丝袜av| 女的被弄到高潮叫床怎么办| 只有这里有精品99| 在线免费观看不下载黄p国产| 久久女婷五月综合色啪小说| 一区二区三区四区激情视频| 人人妻人人爽人人添夜夜欢视频 | 国产视频首页在线观看| 成年女人在线观看亚洲视频| 国产精品人妻久久久影院| 欧美3d第一页| 欧美亚洲 丝袜 人妻 在线| 高清欧美精品videossex| 人妻制服诱惑在线中文字幕| 色视频在线一区二区三区| 大又大粗又爽又黄少妇毛片口| 亚洲欧美成人综合另类久久久| 久久久久性生活片| 在线观看av片永久免费下载| 春色校园在线视频观看| 欧美日韩综合久久久久久| 久久99热这里只频精品6学生| 人妻 亚洲 视频| 中国三级夫妇交换| av黄色大香蕉| 精品一区二区免费观看| 国产爽快片一区二区三区| 国产精品成人在线| 人妻少妇偷人精品九色| av在线蜜桃| 欧美成人一区二区免费高清观看| 国产精品一区二区三区四区免费观看| 青青草视频在线视频观看| a级一级毛片免费在线观看| 成人国产av品久久久| 日韩成人av中文字幕在线观看| 成人一区二区视频在线观看| 一级片'在线观看视频| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 蜜桃久久精品国产亚洲av| 日本与韩国留学比较| 亚洲精品日本国产第一区| 亚洲成色77777| 校园人妻丝袜中文字幕| 一个人看视频在线观看www免费| 精品一区在线观看国产| 另类亚洲欧美激情| 多毛熟女@视频| 欧美一级a爱片免费观看看| 美女国产视频在线观看| 国产乱人偷精品视频| 欧美bdsm另类| 在线亚洲精品国产二区图片欧美 | 99精国产麻豆久久婷婷| 色视频在线一区二区三区| 国产欧美日韩精品一区二区| 久久人人爽人人爽人人片va| 大又大粗又爽又黄少妇毛片口| 中国三级夫妇交换| 观看免费一级毛片| 日韩大片免费观看网站| 乱码一卡2卡4卡精品| 国产精品秋霞免费鲁丝片| 欧美日韩亚洲高清精品| 26uuu在线亚洲综合色| 干丝袜人妻中文字幕| 亚洲精品乱码久久久久久按摩| 在现免费观看毛片| 国产有黄有色有爽视频| 中文字幕精品免费在线观看视频 | 国产老妇伦熟女老妇高清| 国产有黄有色有爽视频| 美女cb高潮喷水在线观看| 欧美最新免费一区二区三区| 王馨瑶露胸无遮挡在线观看| 最近中文字幕2019免费版| av卡一久久| 国产成人freesex在线| 久久人人爽av亚洲精品天堂 | 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 国产精品女同一区二区软件| 91狼人影院| 国产精品一区二区在线不卡| 97超碰精品成人国产| 男女国产视频网站| 日日摸夜夜添夜夜添av毛片| 成人影院久久| 亚洲性久久影院| 亚洲国产精品一区三区| 人妻一区二区av| 亚洲av电影在线观看一区二区三区| 欧美变态另类bdsm刘玥| 搡老乐熟女国产| 99久久精品一区二区三区| 麻豆国产97在线/欧美| 成人一区二区视频在线观看| 久久久精品94久久精品| kizo精华| 人人妻人人添人人爽欧美一区卜 | 欧美国产精品一级二级三级 | 国产爱豆传媒在线观看| 亚洲精品第二区| 精品人妻偷拍中文字幕| 免费在线观看成人毛片| 亚洲精品,欧美精品| 青春草国产在线视频| 亚洲av中文av极速乱| 在线免费观看不下载黄p国产| 亚洲精品亚洲一区二区| 国产69精品久久久久777片| 汤姆久久久久久久影院中文字幕| 欧美丝袜亚洲另类| 只有这里有精品99| 如何舔出高潮| 99久久精品热视频| 国产一区有黄有色的免费视频| 我的女老师完整版在线观看| 纯流量卡能插随身wifi吗| 蜜臀久久99精品久久宅男| 99热全是精品| 精品人妻视频免费看| 婷婷色av中文字幕| 免费播放大片免费观看视频在线观看| 亚洲精品一区蜜桃| 亚洲av.av天堂| 又爽又黄a免费视频| 七月丁香在线播放| 中文字幕久久专区| 伦理电影免费视频| 欧美老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 内射极品少妇av片p| 欧美国产精品一级二级三级 | 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 简卡轻食公司| 国产精品一区二区在线观看99| 在线观看国产h片| 亚洲av综合色区一区| 国产av精品麻豆| 国产综合精华液| 日韩一区二区视频免费看| 精品一区二区免费观看| 干丝袜人妻中文字幕| 成年美女黄网站色视频大全免费 | 身体一侧抽搐| av卡一久久| 久久久久久久久久成人| 久久久久精品性色| 国产白丝娇喘喷水9色精品| 久久久久人妻精品一区果冻| 日本欧美视频一区| 99热网站在线观看| 亚洲国产成人一精品久久久| 国内少妇人妻偷人精品xxx网站| 天天躁日日操中文字幕| 人人妻人人澡人人爽人人夜夜| av国产久精品久网站免费入址| 欧美日韩综合久久久久久| 大香蕉97超碰在线| 亚洲精品视频女| 麻豆乱淫一区二区| 欧美bdsm另类| 观看免费一级毛片| 黑人猛操日本美女一级片| 中文资源天堂在线| av国产精品久久久久影院| 亚洲av二区三区四区| 最近手机中文字幕大全| 天天躁夜夜躁狠狠久久av| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 国产一级毛片在线| 美女国产视频在线观看| av在线老鸭窝| 亚洲精品国产成人久久av| 亚洲,欧美,日韩| 91精品国产国语对白视频| 99视频精品全部免费 在线| 男人和女人高潮做爰伦理| 男的添女的下面高潮视频| 久久精品国产亚洲网站| 在线观看一区二区三区| 亚洲av成人精品一区久久| 亚洲内射少妇av| 久久亚洲国产成人精品v| 综合色丁香网| 三级经典国产精品| 国产高清有码在线观看视频| 老司机影院毛片| 国产亚洲精品久久久com| 精品人妻偷拍中文字幕| 蜜桃亚洲精品一区二区三区| 久久久久久伊人网av| 免费人成在线观看视频色| 一级毛片 在线播放| videos熟女内射| 亚洲内射少妇av| 国产成人午夜福利电影在线观看| 国产av国产精品国产| 99视频精品全部免费 在线| 麻豆成人午夜福利视频| 国产真实伦视频高清在线观看| 高清毛片免费看| 五月玫瑰六月丁香| 色5月婷婷丁香| 亚洲成人手机| 国产免费一级a男人的天堂| 欧美精品亚洲一区二区| 美女主播在线视频| 国产大屁股一区二区在线视频| 九色成人免费人妻av| 亚洲性久久影院| 美女内射精品一级片tv| 午夜福利高清视频| 久久久久久久大尺度免费视频| 99久久中文字幕三级久久日本| av在线app专区| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片我不卡| 亚洲不卡免费看| 国产精品福利在线免费观看| 午夜激情福利司机影院| 欧美区成人在线视频| 观看美女的网站| 亚洲国产精品999| 免费av中文字幕在线| 男女无遮挡免费网站观看| 五月玫瑰六月丁香| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 久久精品国产鲁丝片午夜精品| 亚洲欧洲国产日韩| 欧美极品一区二区三区四区| 欧美日韩精品成人综合77777| 超碰97精品在线观看| 亚洲欧美精品专区久久| 性色av一级| 只有这里有精品99| 麻豆精品久久久久久蜜桃| 午夜日本视频在线| a 毛片基地| 国产女主播在线喷水免费视频网站| 久久久色成人| 身体一侧抽搐| 亚洲四区av| 国产成人精品久久久久久| 亚洲四区av| 国产精品一区二区性色av| 在现免费观看毛片| 午夜福利网站1000一区二区三区| 国产精品熟女久久久久浪| 国产黄色免费在线视频| 直男gayav资源| 多毛熟女@视频| 精品酒店卫生间| 男男h啪啪无遮挡| 免费av不卡在线播放| 最近手机中文字幕大全| 免费高清在线观看视频在线观看| 极品教师在线视频| 亚洲av成人精品一二三区| 亚洲精品一区蜜桃| 噜噜噜噜噜久久久久久91| av福利片在线观看| 一级黄片播放器| 久久99热这里只有精品18| 赤兔流量卡办理| 99久久综合免费| 亚洲精品视频女| 亚洲av在线观看美女高潮| 丰满乱子伦码专区| 久久国产精品大桥未久av | 中文字幕免费在线视频6| 一边亲一边摸免费视频| 精品酒店卫生间| 国产精品一区www在线观看| 欧美bdsm另类| 久久这里有精品视频免费| 免费观看无遮挡的男女| av网站免费在线观看视频| 国产精品蜜桃在线观看| 啦啦啦在线观看免费高清www| 人人妻人人看人人澡| av国产久精品久网站免费入址| 这个男人来自地球电影免费观看 | 免费看av在线观看网站| 一级毛片 在线播放| 尾随美女入室| 免费观看的影片在线观看| 国产一区二区三区av在线| 国产日韩欧美在线精品| 一级毛片久久久久久久久女| 久久人妻熟女aⅴ| 久久久久性生活片| 亚洲av.av天堂| 亚洲性久久影院| 国产精品熟女久久久久浪| 日本午夜av视频| 夜夜爽夜夜爽视频| 久久 成人 亚洲| 日韩av不卡免费在线播放| 精品亚洲乱码少妇综合久久| 在线观看国产h片| 亚洲电影在线观看av| 欧美3d第一页| 国产亚洲欧美精品永久| 国产黄片美女视频| 麻豆乱淫一区二区| 亚州av有码| 人妻制服诱惑在线中文字幕| 人妻少妇偷人精品九色| 久久久成人免费电影| 尾随美女入室| 亚洲精品日本国产第一区| 国内揄拍国产精品人妻在线| 免费看不卡的av| 美女cb高潮喷水在线观看| av免费观看日本| 人妻 亚洲 视频| 一级毛片我不卡| 亚洲熟女精品中文字幕| 国产欧美日韩精品一区二区| videossex国产| 国产午夜精品久久久久久一区二区三区| 成人美女网站在线观看视频| 一二三四中文在线观看免费高清| 日韩一本色道免费dvd| 乱码一卡2卡4卡精品| 久久久久久久精品精品| 我的老师免费观看完整版| 亚洲精品乱码久久久v下载方式| 成人18禁高潮啪啪吃奶动态图 | 久久久a久久爽久久v久久| 嫩草影院新地址| 日韩av不卡免费在线播放| 国产亚洲av片在线观看秒播厂| 精品少妇久久久久久888优播| 男的添女的下面高潮视频| av黄色大香蕉| 亚洲精品久久午夜乱码| 91午夜精品亚洲一区二区三区| 高清在线视频一区二区三区| 人人妻人人添人人爽欧美一区卜 | 18禁裸乳无遮挡免费网站照片| 久久国产乱子免费精品| 亚洲av中文字字幕乱码综合| 国产在视频线精品| 久久精品国产自在天天线| 91精品一卡2卡3卡4卡| 老司机影院成人| 少妇人妻久久综合中文| 国产毛片在线视频| 国产淫片久久久久久久久| 99热网站在线观看|