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

    中華絨螯蟹野生群體和不同水系人工選育群體的遺傳多樣性分析*

    2019-12-06 01:23:02李晶晶陳麗梅耿緒云董學(xué)旺孫金生
    漁業(yè)科學(xué)進(jìn)展 2019年6期
    關(guān)鍵詞:里海微衛(wèi)星雜合

    李晶晶 陳麗梅 耿緒云 董學(xué)旺 孫金生

    中華絨螯蟹野生群體和不同水系人工選育群體的遺傳多樣性分析*

    李晶晶1,2陳麗梅3耿緒云2董學(xué)旺2孫金生1①

    (1. 天津師范大學(xué)生命科學(xué)學(xué)院 天津 300387;2. 天津市水生動(dòng)物疫病預(yù)防控制中心 天津 300221;3. 天津農(nóng)學(xué)院水產(chǎn)學(xué)院 天津市水產(chǎn)生態(tài)及養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室 天津 300384)

    本研究利用10個(gè)微衛(wèi)星分子標(biāo)記分析了中華絨螯蟹() 3個(gè)不同水系人工選育群體(“長(zhǎng)江1號(hào)”、“光合1號(hào)”和七里海河蟹)和1個(gè)海河流域自然群體的遺傳多樣性和遺傳分化水平。結(jié)果顯示,10個(gè)位點(diǎn)在4個(gè)群體中的等位基因數(shù)()為3~17,平均等位基因數(shù)為8.5~9.7,平均期望雜合度為0.720~0.745,平均觀測(cè)雜合度為0.566~0.661,平均多態(tài)信息含量為0.687~0.716,近交系數(shù)(is)范圍為–0.080~0.827。在40個(gè)群體–位點(diǎn)組合中,有13個(gè)群體–位點(diǎn)組合顯著偏離哈迪–溫伯格平衡(<0.05)。遺傳多樣性分析結(jié)果顯示,與海河自然群體相比,3個(gè)人工選育群體遺傳多樣性水平略有降低,但仍保持在較高水平,具有較大的選育潛力。遺傳分化分析結(jié)果顯示,群體間遺傳分化指數(shù)(st)范圍為0.015~0.075,遺傳相似度為0.7702~0.9401,遺傳距離為0.0617~0.2611?;贜ei’s遺傳距離構(gòu)建了群體UPGMA系統(tǒng)進(jìn)化樹,自然群體和“光合1號(hào)”聚為一支,而七里海河蟹群體單獨(dú)聚為一支。綜上所述,4個(gè)中華絨螯蟹群體間的遺傳分化水平較低,群體遺傳多樣性較高。本研究將為中華絨螯蟹選育繁育和種質(zhì)資源利用與管理等提供理論基礎(chǔ)。

    中華絨螯蟹;遺傳多樣性;遺傳分化;微衛(wèi)星

    中華絨螯蟹()俗稱河蟹,自然分布于我國(guó)長(zhǎng)江、黃河、遼河和甌江等河川及附近地區(qū)(Ma, 2007; 肖起珍等, 2017)。自20世紀(jì)70年代末突破人工育苗技術(shù)難關(guān),經(jīng)過30余年的發(fā)展,中華絨螯蟹養(yǎng)殖遍布我國(guó)遼寧、河北、天津、山東、江蘇、浙江等沿海省份及廣大內(nèi)陸地區(qū),成為我國(guó)重要經(jīng)濟(jì)甲殼動(dòng)物之一(李晶晶等, 2016)。隨著中華絨螯蟹養(yǎng)殖產(chǎn)業(yè)的規(guī)?;l(fā)展,過度捕撈野生自然資源,無序引進(jìn)異地苗種,大大加劇了中華絨螯蟹種質(zhì)退化和混雜,已成為當(dāng)前養(yǎng)殖生產(chǎn)發(fā)展的瓶頸(王成輝等, 2002)。研究和保護(hù)中華絨螯蟹種質(zhì)資源,對(duì)河蟹養(yǎng)殖業(yè)的持續(xù)穩(wěn)定發(fā)展具有重要意義。近年來,科研院所和養(yǎng)殖企業(yè)積極開展中華絨螯蟹良種人工選育工作,已培育出5個(gè)國(guó)家級(jí)新品種(“光合1號(hào)”、“長(zhǎng)江1號(hào)”、“長(zhǎng)江2號(hào)”、“江海21”和“諾亞1號(hào)”)。

    天津市地處環(huán)渤海中部,具有得天獨(dú)厚的海、淡水養(yǎng)殖環(huán)境條件,中華絨螯蟹養(yǎng)殖歷史悠久。轄區(qū)內(nèi)的七里海是我國(guó)重要的古海岸與古濕地保護(hù)區(qū),獨(dú)特的生態(tài)環(huán)境創(chuàng)出了味道鮮美、膏滿黃肥的天津市特色農(nóng)產(chǎn)品品牌——七里海河蟹(朱廣奇, 2012; 杜紅梅, 2014)。然而,50年代末期,天津海河入海口處建設(shè)了防潮閘,切斷了中華絨螯蟹的生殖洄游通道,導(dǎo)致本地區(qū)自然資源逐年減少。2000年,天津市水生動(dòng)物疫病預(yù)防控制中心采集七里海核心區(qū)中華絨螯蟹野生群體作為選育基礎(chǔ)群體,利用群體選育技術(shù)連續(xù)選育6代,獲得快速生長(zhǎng)的七里海河蟹群體。

    在傳統(tǒng)人工選育中,由于長(zhǎng)期不引進(jìn)外來親本,經(jīng)過連續(xù)多代選育獲得的新品種或新品系,其種群遺傳多樣性是否有所降低,有待進(jìn)行科學(xué)評(píng)價(jià)。一般認(rèn)為,遺傳多樣性的大小及其群體遺傳結(jié)構(gòu)跟一個(gè)物種的進(jìn)化潛力和抵御不良環(huán)境的能力密切相關(guān)(李博, 2000)。因此,分析檢測(cè)育種群體的遺傳特征,對(duì)于中華絨螯蟹選育效果預(yù)測(cè)、種質(zhì)資源保護(hù)以及養(yǎng)殖產(chǎn)業(yè)健康持續(xù)發(fā)展意義重大。

    目前,天津地區(qū)中華絨螯蟹資源主要有人工選育群體“七里海河蟹”、“長(zhǎng)江1號(hào)”、“光合1號(hào)”以及海河流域自然群體。本研究利用10個(gè)多態(tài)性豐富的微衛(wèi)星分子標(biāo)記,比較分析了天津地區(qū)4個(gè)中華絨螯蟹群體的遺傳多樣性和群體間的遺傳分化,旨在為開展中華絨螯蟹選育和種質(zhì)資源利用與管理等工作提供理論基礎(chǔ)。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料

    本實(shí)驗(yàn)共采集4個(gè)中華絨螯蟹群體,每個(gè)群體樣品數(shù)量為60只(表1)。七里海河蟹群體于2012年采自天津市蟹源水產(chǎn)養(yǎng)殖有限公司,為天津地區(qū)海河流域人工選育的第6代成蟹群體;自然群體于2012年采自海河入???,是天津地區(qū)海河流域野生成蟹種群;“長(zhǎng)江1號(hào)”和“光合1號(hào)”于2013年分別采自江蘇省淡水水產(chǎn)研究所和盤錦光合蟹業(yè)有限公司,為長(zhǎng)江水系和遼河水系人工選育幼蟹群體。采集的4個(gè)中華絨螯蟹群體同屬奇數(shù)年群體。分別剪取每只中華絨螯蟹步足肌肉組織,置于–20℃冰箱保存?zhèn)溆谩?/p>

    1.2 微衛(wèi)星引物

    本研究選用了10個(gè)已發(fā)表的多態(tài)性較高的微衛(wèi)星位點(diǎn),用于中華絨螯蟹4個(gè)群體的遺傳分析(Li, 2016)。引物由生工生物工程(上海)股份有限公司合成,引物序列及特征見表2。

    1.3 DNA提取與檢測(cè)

    采用天根生化科技有限公司的海洋動(dòng)物組織基因組DNA提取試劑盒,提取中華絨螯蟹肌肉組織的總DNA,并用1%瓊脂糖凝膠電泳檢測(cè)DNA質(zhì)量,紫外分光光度法測(cè)定DNA濃度和純度。用滅菌雙蒸水將DNA稀釋成100 μg/ml后,–20℃保存?zhèn)溆谩?/p>

    1.4 PCR反應(yīng)體系、反應(yīng)程序及擴(kuò)增產(chǎn)物的檢測(cè)

    PCR反應(yīng)在10 μl體系中進(jìn)行,包括1×buffer、2 mmol/L MgCl2、0.2 mmol/L dNTP、1 μmol/L引物、0.25 UDNA聚合酶(TaKaRa)和100 ng DNA模板。PCR程序:94℃預(yù)變性3 min;94℃變性1 min,退火1 min,72℃延伸1 min,35個(gè)循環(huán);最后72℃延伸5 min。PCR擴(kuò)增產(chǎn)物用6%的變性聚丙烯酰胺凝膠電泳檢測(cè),硝酸銀法染色,10 bp DNA ladder (Invitrogen)作為Marker檢測(cè)擴(kuò)增片段大小。

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

    使用Microsatellite AnalyseR (MSA)軟件計(jì)算各微衛(wèi)星位點(diǎn)在4個(gè)群體中的等位基因數(shù)()、期望雜合度(e)和觀測(cè)雜合度(o) (Dieringer, 2003)。用公式e=1/∑i2計(jì)算有效等位基因數(shù)(e) (i表示每個(gè)位點(diǎn)第個(gè)等位基因的頻率) (Kimura, 1964)。FSTAT 2.9.3軟件(Goudet, 2001)計(jì)算每個(gè)群體每個(gè)位點(diǎn)的近交系數(shù)(is)。應(yīng)用PIC-CALC 0.6軟件計(jì)算多態(tài)信息含量(PIC)。Genepop 4.2軟件檢測(cè)10個(gè)位點(diǎn)的哈迪-溫伯格平衡(HWE)偏離情況,并比較分析各位點(diǎn)在群體間等位基因頻率(Raymond, 1995)。MSA軟件計(jì)算st值(Weir, 1984)。應(yīng)用SPAGeDi 1.1軟件(Hardy, 2002)對(duì)無效等位基因進(jìn)行校正,Bonferroni檢驗(yàn)對(duì)顯著性水平的多點(diǎn)檢測(cè)進(jìn)行校正(Rice, 1989)。使用Popgene Version 1.31計(jì)算Nei’s遺傳距離(Nei, 1978),利用MEGA 5軟件構(gòu)建群體間UPGMA系統(tǒng)進(jìn)化樹(Tamura, 2011)。

    表1 中華絨螯蟹群體的取樣時(shí)間、地點(diǎn)、群體類型和樣本數(shù)

    Tab.1 The sampling time, location, population and size of E. sinensis

    表2 中華絨螯蟹10個(gè)微衛(wèi)星位點(diǎn)特征

    Tab.2 Characteristics of ten microsatellite loci of E. sinensis

    2 結(jié)果與分析

    2.1 微衛(wèi)星位點(diǎn)多態(tài)性及群體遺傳多樣性

    應(yīng)用選取的10對(duì)微衛(wèi)星引物,對(duì)4個(gè)中華絨螯蟹群體共240個(gè)個(gè)體進(jìn)行擴(kuò)增分析,遺傳多樣性參數(shù)見表3。等位基因數(shù)范圍為3~17,其中,位點(diǎn)SSR352和SSR1289等位基因數(shù)最多(=17),分別出現(xiàn)在七里海河蟹群體和自然群體中。七里海河蟹、自然群體、“長(zhǎng)江1號(hào)”和“光合1號(hào)”的平均等位基因數(shù)分別為9.4、9.7、9.0和8.5;平均觀測(cè)雜合度為0.658、0.661、0.566和0.618;平均期望雜合度為0.720、0.745、0.732和0.720;平均多態(tài)信息含量為0.691、0.716、0.694和0.687,所有參數(shù)的最高平均值均出現(xiàn)在自然群體中。所有位點(diǎn)的平均觀測(cè)雜合度均略低于平均期望雜合度,4個(gè)群體均具有較高的平均期望雜合度(>0.5)和平均多態(tài)信息含量(>0.5)。

    10個(gè)位點(diǎn)在4個(gè)群體中的近交系數(shù)is范圍為–0.080~0.827。在七里海河蟹群體中,位點(diǎn)SSR1240、SSR1053和SSR709的近交系數(shù)為負(fù)值;在“光合1號(hào)”群體中,位點(diǎn)SSR1289、SSR1126和SSR709的近交系數(shù)為負(fù)值;在自然群體中,位點(diǎn)SSR709的近交系數(shù)為負(fù)值,其余均為正值。經(jīng)過Bonferroni校正后,10個(gè)位點(diǎn)在4個(gè)群體中,有13個(gè)顯著偏離哈迪-溫伯格平衡(<0.05)。其中,“長(zhǎng)江1號(hào)”中偏離HWE平衡的位點(diǎn)最多,為5個(gè),七里海河蟹和海河自然群體均有3個(gè)位點(diǎn)出現(xiàn)顯著偏離,“光合1號(hào)”偏離位點(diǎn)最少,為2個(gè)。

    2.2 群體間的遺傳分化

    4個(gè)群體間的st值存在顯著差異(<0.05),結(jié)果如表4所示,七里海河蟹群體與“長(zhǎng)江1號(hào)”群體、“光合1號(hào)”群體的遺傳分化程度較高,達(dá)到中等水平(st>0.05);七里海河蟹群體與自然群體、自然群體與“長(zhǎng)江1號(hào)”群體、自然群體與“光合1號(hào)”群體、“長(zhǎng)江1號(hào)”群體與“光合1號(hào)”群體的遺傳分化程度處于低水平(st<0.05)。分別對(duì)每個(gè)位點(diǎn)進(jìn)行等位基因頻率的群體間兩兩比較,60個(gè)比較中有28個(gè)差異顯著(<0.05) (表6)。

    表3 4個(gè)中華絨螯蟹群體10個(gè)微衛(wèi)星位點(diǎn)的遺傳多樣性

    Tab.3 Genetic diversities of ten microsatellite loci in four populations of E. sinensis

    注:*表示顯著偏離哈迪-溫伯格平衡(<0.05)

    Note: * means significant deviation from Hardy-Weinberg equilibrium (<0.05)

    4個(gè)群體兩兩間的遺傳相似度范圍為0.7702~ 0.9401,遺傳距離(c)范圍為0.0617~0.2611,群體間遺傳距離越大的遺傳相似性越低(表5)?;贜ei’s遺傳距離構(gòu)建的群體UPGMA系統(tǒng)進(jìn)化樹顯示,自然群體與“光合1號(hào)”聚為一支后,再與“長(zhǎng)江1號(hào)”聚為一支,而七里海河蟹群體單獨(dú)聚為一支(圖1)。

    表4 群體兩兩間的遺傳分化指數(shù)(st)

    Tab.4 Genetic differentiation index (Fst) between any two populations

    注:*表示差異顯著(<0.05),下同

    Note: * means significant difference (<0.05). The same as below

    表5 4個(gè)中華絨螯蟹群體的遺傳相似度和遺傳距離

    Tab.5 Nei’s genetic identity and genetic distance of the four populations of E. sinensis

    注:對(duì)角線以上為遺傳相似度;對(duì)角線以下為遺傳距離

    Note: The diagonal above was genetic identity; the diagonal below was genetic distance

    圖1 基于Nei’s (1978)遺傳距離構(gòu)建的4個(gè)中華絨螯蟹群體的UPGMA系統(tǒng)進(jìn)化樹

    表6 10個(gè)微衛(wèi)星位點(diǎn)等位基因頻率的群體間兩兩比較

    Tab.6 Comparisons of allele frequency in ten microsatellite loci between any two populations

    續(xù)表6

    3 討論

    3.1 中華絨螯蟹自然和人工選育群體遺傳多樣性

    遺傳多樣性是生物多樣性的基礎(chǔ),對(duì)物種的生存、繁衍和進(jìn)化至關(guān)重要(張榮良等, 2016)。一個(gè)遺傳多樣性豐富的種群,往往具有較高的適應(yīng)能力、生存能力和進(jìn)化潛能,在育種和遺傳改良上也更具有潛力,所以對(duì)遺傳多樣性的準(zhǔn)確評(píng)估是生物多樣性保護(hù)和可持續(xù)利用的基礎(chǔ)和前提。人工選擇育種就是在保持選育群體具有一定遺傳多樣性的基礎(chǔ)上,篩選獲得具有目標(biāo)性狀的選育新品種或新品系(張榮良等, 2016)。大量相關(guān)研究表明,養(yǎng)殖群體因人工選擇作用和瓶頸效應(yīng)等原因,其遺傳多樣性往往低于未經(jīng)選育的野生群體。張輝等(2010)和王軍等(2018)分別應(yīng)用mtDNA控制區(qū)序列和微衛(wèi)星標(biāo)記分析了中國(guó)對(duì)蝦()養(yǎng)殖群體與野生群體,顯示野生群體的基因多樣性高于養(yǎng)殖群體。肖起珍等(2017)利用29對(duì)SSR標(biāo)記對(duì)中華絨螯蟹長(zhǎng)江水系野生和人工繁育大眼幼體的遺傳多樣性進(jìn)行了分析,發(fā)現(xiàn)野生群體的平均等位基因數(shù)、平均有效等位基因數(shù)和平均期望雜合度均高于人工繁育群體,香農(nóng)信息指數(shù)體現(xiàn)的群體遺傳多樣性野生群體也高于人工繁育群體。在魚類群體遺傳研究中也獲得了相似的結(jié)果,曾慶凱等(2017)應(yīng)用SSR標(biāo)記比較翹嘴鱖()湖南洞庭湖野生群體和廣東、安徽養(yǎng)殖群體的遺傳多樣性,分析結(jié)果也顯示野生群體的遺傳多樣性略高于養(yǎng)殖群體。

    遺傳多樣性主要用多態(tài)信息含量(PIC)和雜合度()來評(píng)估。其中,多態(tài)信息含量用于評(píng)價(jià)標(biāo)志基因(或標(biāo)志序列)的多態(tài)性大小。PIC>0.5的位點(diǎn)為高多態(tài)性位點(diǎn),0.25

    由于連續(xù)多代人工選育,不引進(jìn)外來親本,所以造成了中華絨螯蟹選育群體遺傳多樣性的降低,使選育群體的遺傳多樣性略低于自然群體。但在人工選育過程中,增加有效繁育親本數(shù)量,嚴(yán)格控制選育群體的近交衰退,在一定程度上保證了選育群體的遺傳穩(wěn)定性,使中華絨螯蟹的人工選育群體經(jīng)過累代繁殖多年之后,依然保持著較高的遺傳多樣性。

    3.2 群體間的遺傳分化

    遺傳分化指數(shù)(st)和遺傳距離(c)是衡量群體間遺傳差異的重要參數(shù)。Wright (1978)認(rèn)為,st值在0~0.05之間,說明群體間的遺傳差異很?。籹t值在0.05~0.15之間,說明群體間發(fā)生了中等水平的遺傳分化;在0.15~0.25之間,說明群體間存在較大的遺傳差異;當(dāng)st>0.25時(shí),說明群體間高度遺傳分化。本研究中,4個(gè)群體兩兩間比較,七里海河蟹群體與“長(zhǎng)江1號(hào)”群體、“光合1號(hào)”群體間的遺傳分化達(dá)到中等水平(st>0.05);其他比較組間的遺傳分化程度處于低水平(st<0.05)。兩兩群體間的Nei’s無偏遺傳距離為0.0617~0.2611,群體間遺傳距離較小。這與肖起珍等(2017)對(duì)長(zhǎng)江水系中華絨螯蟹野生與人工繁殖群體遺傳研究的結(jié)果相似。

    群體系統(tǒng)進(jìn)化樹顯示,自然群體和“光合1號(hào)”聚為一支后,再與“長(zhǎng)江1號(hào)”聚為一支,而七里海河蟹群體單獨(dú)聚為一支。這可能是由天津地區(qū)中華絨螯蟹養(yǎng)殖產(chǎn)業(yè)現(xiàn)狀導(dǎo)致的。自20世紀(jì)中后期,天津海河入??诮ㄔO(shè)了防潮閘后,海河流域的中華絨螯蟹自然資源逐年減少,同時(shí)養(yǎng)殖產(chǎn)業(yè)迅速發(fā)展。但天津地區(qū)遺傳育種研究起步晚,育種技術(shù)相對(duì)落后,本地人工繁育蟹苗量少,遠(yuǎn)遠(yuǎn)無法滿足當(dāng)?shù)仞B(yǎng)殖生產(chǎn)需求。所以,養(yǎng)殖企業(yè)大量引進(jìn)遼河、長(zhǎng)江等其他水系苗種,以“光合1號(hào)”和“長(zhǎng)江1號(hào)”居多。在養(yǎng)殖生產(chǎn)過程中,由于養(yǎng)殖群體個(gè)體逃逸,可能會(huì)與自然群體發(fā)生基因交流,從而導(dǎo)致了異地引進(jìn)群體與自然群體間的遺傳距離變小,打破了各水系中華絨螯蟹的原始分布格局。

    本研究以中華絨螯蟹3個(gè)水系人工選育群體和海河流域自然群體為對(duì)象,應(yīng)用微衛(wèi)星分子標(biāo)記技術(shù)分析各群體的遺傳多樣性和群體間的遺傳分化。結(jié)果表明,4個(gè)群體的遺傳分化程度處于中低水平,3個(gè)人工選育群體的遺傳多樣性與自然群體相比有一定程度的降低,但仍保持在較高水平上,具有較大的選育潛力。

    Botstein D, White RL, Skolnick M,. Construction of genetic linkage map in man using restriction fragment length polymorphisms. American Journal of Human Genetics, 1980, 32: 314–331

    Dieringer D, Schl?tterer C. Microsatellite analyser (MSA): A platform independent analysis tool for large microsatellite data sets. Molecular Ecology Notes, 2003, 3(1): 167–169

    Du HM. Brief discussion of Qilihai crab ecological cultivation. Rural Science and Technology, 2014(3): 29 [杜紅梅. 七里海河蟹生態(tài)養(yǎng)殖淺談. 鄉(xiāng)村科技, 2014(3): 29]

    Goudet J. FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). 2001, Available from http://www.unil.ch/izea/softwares/fstat.html

    Hardy OJ, Vekemans X. SPAGeDi: A versatile computer program toanalyze spatial genetic structure at the individual or population levels. Molecular Ecology Notes, 2002, 2(4): 618–620

    Kimura M, Crow JF. The number of alleles that can be maintained in a finite population. Genetics, 1964, 49(4): 725–738

    Li B. Ecology. Beijing: Higher Education Press, 2000, 338–341 [李博. 生態(tài)學(xué). 北京: 高等教育出版社, 2000, 338–341]

    Li JJ, Geng XY, Chen LM,. Identification and characterizationof thirty novel microsatellite DNA markers from the Chinesemitten crabexpressed sequence tags. BMC Research Notes, 2016, 105(9): 1–5

    Li JJ, Geng XY, Dong XW,. Evaluation of visible implant elastomer tags in Chinese mitten crab. Marine Fisheries, 2016, 38(4): 400–405 [李晶晶, 耿緒云, 董學(xué)旺, 等. 中華絨螯蟹可視熒光標(biāo)記應(yīng)用研究. 海洋漁業(yè), 2016, 38(4): 400–405]

    Ma HT, Chang YM, Yu DM,. Microsatellite variations among four populations of. Zoological Research, 2007, 28(2): 126–133

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

    Nei M. Molecular evolutionary genetics. New York: Columbia University Press, 1987

    Qin Y, Shi G, Sun Y. Evaluation of genetic diversity inusing SSR markers. Genetics and Molecular Research, 2013, 12(4): 5833–5841

    Raymond M, Rousset F. Population genetics software for exact tests and ecumenicism. Journal of Heredity, 1995, 86(3): 248–249

    Rice WR. Analyzing tables of statistical tests. Evolution, 1989, 43(1): 223–225

    Takezaki N, Nei M. Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA. Genetics, 1996, 144(1): 389–399

    Tamura K, Peterson D, Peterson N,. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution, 2011, 28(10): 2731–2739

    Wang CH, Li SF. Advances in studies on germplasm in Chinese mitten crab,. Journal of Fishery Sciences of China, 2002, 9(1): 82–86 [王成輝, 李思發(fā). 中華絨螯蟹種質(zhì)研究進(jìn)展. 中國(guó)水產(chǎn)科學(xué), 2002, 9(1): 82–86]

    Wang J, Wang QY, Kong J,. SSR Analysis on genetic diversity in breeding and wild populations of. Progress in Fishery Sciences, 2018, 39(2): 104– 111 [王軍, 王清印, 孔杰, 等. 中國(guó)明對(duì)蝦人工選育群體與野生群體遺傳多樣性的SSR分析. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(2): 104–111]

    Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution, 1984, 38(6): 1358–1370

    Wright S. Evolution and the genetics of population variability within and among natural population. University of Chicago Press, Chicago, 1978, 121–124

    Xiao QZ, Liu Q, Wu XG,. Genetic diversity analysis of wild and cultured megalopa population offrom Yangtze River. Genomics and Applied Biology, 2017, 36(5): 1935–1945 [肖起珍, 劉青, 吳旭干, 等. 長(zhǎng)江水系中華絨螯蟹野生和人工繁殖大眼幼體的遺傳多樣性分析. 基因組學(xué)與應(yīng)用生物學(xué), 2017, 36(5): 1935–1945]

    Zeng QK, Sun CF, Dong JJ,. Analysis of genetic diversity in three different populations of. Genomics and Applied Biology, 2017, 36(8): 3241–3250 [曾慶凱, 孫成飛, 董浚鍵, 等. 翹嘴鱖3個(gè)不同群體的遺傳多樣性分析. 基因組學(xué)與應(yīng)用生物學(xué), 2017, 36(8): 3241–3250]

    Zhang H, Gao TX, Zhuang ZM,. Comparative analysis of the mitochondrial control region between the cultured and wild population of Chinese shrimp (). Journal of Fisheries of China, 2010, 34(8): 1149–1155 [張輝, 高天翔, 莊志猛, 等. 中國(guó)對(duì)蝦養(yǎng)殖群體與野生群體線粒體控制區(qū)序列的比較. 水產(chǎn)學(xué)報(bào), 2010, 34(8): 1149–1155]

    Zhang RL, Wang WJ, Feng YW,. Assessment of genetic variability and microsatellite analysis of Pacific oyster () after artificial selection of the shell width. Progress in Fishery Sciences, 2016, 37(4): 90–96 [張榮良, 王衛(wèi)軍, 馮艷微, 等. 長(zhǎng)牡蠣()殼寬快速生長(zhǎng)選育群體遺傳多樣性及遺傳結(jié)構(gòu)的微衛(wèi)星標(biāo)記分析. 漁業(yè)科學(xué)進(jìn)展, 2016, 37(4): 90–96]

    Zhu GQ. Integration and application of key technologies for ecological cultivation of Qilihai crab. Tianjin Fisheries, 2012(1): 35–39 [朱廣奇. 七里海河蟹生態(tài)養(yǎng)殖關(guān)鍵技術(shù)集成與應(yīng)用. 天津水產(chǎn), 2012(1): 35–39]

    Genetic Diversity Analysis of Wild and Cultured Population offrom Different Water Systems

    LI Jingjing1,2, CHEN Limei3, GENG Xuyun2, DONG Xuewang2, SUN Jinsheng1①

    (1. College of Life Sciences, Tianjin Normal University, Tianjin 300387; 2. Tianjin Diseases Prevention and Control Center of Aquatic Animals, Tianjin 300221; 3. Tianjin Key Laboratory of Aqua-Ecology and Aquaculture, Fisheries College, Tianjin Agricultural University, Tianjin 300384)

    The Chinese mitten crabis one of the most economically important aquaculture species in China due to its taste and nutritional value. In this study, the genetic status of both wild and cultivated stocks ofwas assessed using simple sequence repeat (SSR) markers. Ten microsatellite markers reported in previous studies were selected to estimate the level of genetic diversity within one Haihe natural population and three selected artificial populations (Yangtze No.1, Guanghe No.1, and Qilihai crab) and to compare the degree of genetic differentiation among them. Unique PCR products and high levels of polymorphism were observed for all loci. For the four populations of, the number of alleles per locus ranged from 3 to 17, with a mean allele number of 8.5 to 9.7. The mean observed heterozygosities ranged from 0.566 to 0.661, while the mean expected heterozygosities ranged from 0.720 to 0.745. Additionally, the mean polymorphism information content ranged from 0.687 to 0.716. Furthermore, the inbreeding coefficient (F) ranged from-0.080 to 0.827, and there were 13 cases that significantly deviated from the Hardy-Weinberg equilibrium (<0.05) among the 40 population-locus cases (4 populations × 10 loci). The results of the microsatellite survey indicated that the genetic diversity of the three artificially selected populations was slightly lower than that of the natural population, but still maintained at a high level for further breeding. The analysis of genetic differentiation showed thatFranged from 0.015 to 0.075, genetic similarity ranged from 0.7702 to 0.9401, and genetic distance ranged from 0.0617 to 0.2611. UPGMA phylogenetic analysis divided these populations into two groups. Natural population and Guanghe No.1 stocks were allocated to the same cluster, and the Qilihai crab was allocated to another cluster. In conclusion, the four populations ofpossessed high genetic diversity but low to moderate levels of genetic differentiation. The information on the genetic variation and differentiation obtained in this study will provide a theoretical basis for further breeding and utilization of germplasm resources of.

    ; Genetic diversity; Genetic differentiation; Microsatellite

    S968.2

    A

    2095-9869(2019)06-0105-09

    10.19663/j.issn2095-9869.20180831002

    http://www.yykxjz.cn/

    李晶晶, 陳麗梅, 耿緒云, 董學(xué)旺, 孫金生. 中華絨螯蟹野生群體和不同水系人工選育群體的遺傳多樣性分析. 漁業(yè)科學(xué)進(jìn)展, 2019, 40(6): 105–113

    Li JJ, Chen LM, Geng XY, Dong XW, Sun JS. Genetic diversity analysis of wild and cultured population offrom different water systems. Progress in Fishery Sciences, 2019, 40(6): 105–113

    * 天津市財(cái)政資金項(xiàng)目(七里海優(yōu)質(zhì)河蟹繁養(yǎng)技術(shù)示范應(yīng)用)、天津市水產(chǎn)局科研推廣項(xiàng)目(J2013-21)、天津市水產(chǎn)產(chǎn)業(yè)技術(shù)體系專項(xiàng)(ITTFRS2017002)、天津市種業(yè)科技重大專項(xiàng)(15ZXZYNC00050)和天津市自然科學(xué)基金項(xiàng)目(16JCQNJC09500)共同資助 [This work was supported by Project of Tianjin Financial Funds (Demonstration and Application of Qilihai Crab Breeding and Culture Technology), Research and Extension Project of Tianjin Aquatic Products Bureau (J2013-21), Tianjin Aquaculture Research System (ITTFRS2017002), Major Project for Tianjin Seed Technology (15ZXZYNC00050), and Tianjin Natural Science Foundation (16JCQNJC09500)]. 李晶晶, E-mail: jingjingli206@163.com

    孫金生,研究員,E-mail: jinshsun@163.com

    2018-08-31,

    2018-10-20

    SUN Jinsheng, E-mail: jinshsun@163.com

    (編輯 馮小花)

    猜你喜歡
    里海微衛(wèi)星雜合
    梅花鹿基因組微衛(wèi)星分布特征研究
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    這才是擱置爭(zhēng)議的正確打開方式
    里海適用《聯(lián)合國(guó)海洋法公約》嗎?(答讀者問)
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國(guó)“一箭雙星”成功將“遙感衛(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)化
    個(gè)體基因雜合度對(duì)山西白豬體重和體尺性狀的影響
    色婷婷av一区二区三区视频| 亚洲免费av在线视频| 久久久久久久午夜电影 | 久久国产亚洲av麻豆专区| 别揉我奶头~嗯~啊~动态视频| 亚洲精品美女久久av网站| 在线观看66精品国产| 久久午夜综合久久蜜桃| 久久久国产成人免费| 深夜精品福利| 999精品在线视频| 淫妇啪啪啪对白视频| 99香蕉大伊视频| 露出奶头的视频| 国产一区在线观看成人免费| 黄色成人免费大全| 夫妻午夜视频| 青草久久国产| 在线天堂中文资源库| 自线自在国产av| 欧美人与性动交α欧美软件| 多毛熟女@视频| 最新的欧美精品一区二区| 无限看片的www在线观看| 香蕉丝袜av| 国产一区在线观看成人免费| 亚洲欧美激情在线| 国产激情欧美一区二区| 色婷婷av一区二区三区视频| 精品免费久久久久久久清纯 | 午夜老司机福利片| 精品久久久久久久毛片微露脸| 欧洲精品卡2卡3卡4卡5卡区| 一边摸一边抽搐一进一小说 | 欧美黄色淫秽网站| 中文字幕人妻丝袜制服| 久久精品国产亚洲av香蕉五月 | 久久99一区二区三区| 高清视频免费观看一区二区| 一进一出抽搐gif免费好疼 | 国内久久婷婷六月综合欲色啪| 国产成人影院久久av| 午夜福利免费观看在线| 亚洲片人在线观看| 黄色 视频免费看| 搡老熟女国产l中国老女人| 丰满人妻熟妇乱又伦精品不卡| 黄色片一级片一级黄色片| 国内毛片毛片毛片毛片毛片| 成年动漫av网址| 飞空精品影院首页| www日本在线高清视频| 欧美大码av| 欧美精品人与动牲交sv欧美| 天堂√8在线中文| 欧美精品高潮呻吟av久久| 老熟妇乱子伦视频在线观看| 亚洲 欧美一区二区三区| 欧美一级毛片孕妇| 超色免费av| 黄色视频,在线免费观看| 高清视频免费观看一区二区| 波多野结衣av一区二区av| 一级黄色大片毛片| 亚洲av片天天在线观看| 国产人伦9x9x在线观看| 视频区图区小说| 亚洲精华国产精华精| 精品久久久久久久毛片微露脸| 久久精品国产亚洲av香蕉五月 | 99国产精品一区二区蜜桃av | 国产区一区二久久| 18在线观看网站| 国产亚洲精品久久久久5区| 一级毛片高清免费大全| 午夜老司机福利片| 国产亚洲精品久久久久5区| 国产不卡av网站在线观看| 亚洲美女黄片视频| 国产亚洲av高清不卡| 久久国产亚洲av麻豆专区| 久久午夜综合久久蜜桃| 精品久久蜜臀av无| 99香蕉大伊视频| 久久久久久久久久久久大奶| 五月开心婷婷网| 久久精品aⅴ一区二区三区四区| 国产深夜福利视频在线观看| 国产亚洲精品久久久久久毛片 | 国产精品免费一区二区三区在线 | 宅男免费午夜| 欧美av亚洲av综合av国产av| 一区二区三区激情视频| 免费日韩欧美在线观看| 99精品在免费线老司机午夜| 99久久综合精品五月天人人| 99香蕉大伊视频| 亚洲精品自拍成人| 极品教师在线免费播放| 十八禁网站免费在线| 国产精品香港三级国产av潘金莲| 精品电影一区二区在线| 精品少妇一区二区三区视频日本电影| 国产免费男女视频| 免费在线观看日本一区| 欧美精品av麻豆av| 狂野欧美激情性xxxx| 久9热在线精品视频| 1024视频免费在线观看| 国产伦人伦偷精品视频| 99re6热这里在线精品视频| avwww免费| av中文乱码字幕在线| 国产无遮挡羞羞视频在线观看| 欧美不卡视频在线免费观看 | 99re在线观看精品视频| 国产男靠女视频免费网站| 在线播放国产精品三级| 久久久久久免费高清国产稀缺| 欧美日韩瑟瑟在线播放| 啦啦啦视频在线资源免费观看| 久久热在线av| 国产男女内射视频| 精品欧美一区二区三区在线| 欧美日韩福利视频一区二区| 这个男人来自地球电影免费观看| 国产有黄有色有爽视频| 国产一区二区激情短视频| 99国产精品免费福利视频| 精品人妻在线不人妻| 巨乳人妻的诱惑在线观看| av中文乱码字幕在线| 天天躁夜夜躁狠狠躁躁| 精品一区二区三区视频在线观看免费 | 脱女人内裤的视频| 中文字幕最新亚洲高清| 国产男女内射视频| 免费久久久久久久精品成人欧美视频| e午夜精品久久久久久久| 老司机午夜十八禁免费视频| 飞空精品影院首页| 欧美 日韩 精品 国产| 777久久人妻少妇嫩草av网站| www日本在线高清视频| 香蕉久久夜色| 亚洲av片天天在线观看| 男男h啪啪无遮挡| 看免费av毛片| 亚洲精品av麻豆狂野| 亚洲精品国产色婷婷电影| 午夜精品在线福利| 国产成人影院久久av| 色精品久久人妻99蜜桃| 国产成人精品无人区| 中文字幕另类日韩欧美亚洲嫩草| 国产一区二区三区在线臀色熟女 | 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费成人在线视频| 天天躁夜夜躁狠狠躁躁| 18禁观看日本| 老司机靠b影院| 欧美日韩黄片免| 日本黄色视频三级网站网址 | 一夜夜www| 国内毛片毛片毛片毛片毛片| 国产精品av久久久久免费| 久久精品人人爽人人爽视色| 国产一区二区激情短视频| 久久久水蜜桃国产精品网| 捣出白浆h1v1| 国产一区二区激情短视频| 大香蕉久久网| videosex国产| 亚洲性夜色夜夜综合| 国产精品自产拍在线观看55亚洲 | 丰满迷人的少妇在线观看| 91老司机精品| 精品少妇久久久久久888优播| 丁香六月欧美| 成人精品一区二区免费| 美国免费a级毛片| 亚洲午夜精品一区,二区,三区| 欧美日韩av久久| 妹子高潮喷水视频| 黄色 视频免费看| 久久久久久久国产电影| 中文字幕制服av| 丝瓜视频免费看黄片| 国产亚洲精品第一综合不卡| 久久国产乱子伦精品免费另类| 亚洲国产看品久久| 80岁老熟妇乱子伦牲交| 免费在线观看完整版高清| 一本一本久久a久久精品综合妖精| 天堂√8在线中文| videosex国产| 国产亚洲精品久久久久5区| 久久99一区二区三区| 最近最新免费中文字幕在线| 久久九九热精品免费| 亚洲一卡2卡3卡4卡5卡精品中文| 久久这里只有精品19| 国产97色在线日韩免费| 老司机福利观看| 黄色丝袜av网址大全| 高潮久久久久久久久久久不卡| 日日夜夜操网爽| 午夜91福利影院| 97人妻天天添夜夜摸| 国产视频一区二区在线看| 999久久久精品免费观看国产| 丝袜在线中文字幕| 99久久人妻综合| 99re在线观看精品视频| 少妇被粗大的猛进出69影院| 精品人妻熟女毛片av久久网站| 亚洲中文av在线| 免费观看人在逋| 色精品久久人妻99蜜桃| 午夜精品在线福利| 亚洲人成伊人成综合网2020| 成年女人毛片免费观看观看9 | 亚洲午夜理论影院| 欧美精品啪啪一区二区三区| 亚洲欧美日韩另类电影网站| 悠悠久久av| 亚洲美女黄片视频| 91麻豆精品激情在线观看国产 | 国产主播在线观看一区二区| 99热只有精品国产| 日韩视频一区二区在线观看| 精品午夜福利视频在线观看一区| 法律面前人人平等表现在哪些方面| 精品一区二区三区av网在线观看| 亚洲av片天天在线观看| 又黄又爽又免费观看的视频| 欧美日韩黄片免| 新久久久久国产一级毛片| 国产精品香港三级国产av潘金莲| 99re6热这里在线精品视频| 亚洲性夜色夜夜综合| 中亚洲国语对白在线视频| 99热国产这里只有精品6| 大码成人一级视频| 亚洲欧洲精品一区二区精品久久久| 欧美最黄视频在线播放免费 | 麻豆国产av国片精品| 亚洲成国产人片在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲av成人一区二区三| 亚洲精品一卡2卡三卡4卡5卡| 欧美 日韩 精品 国产| 99re在线观看精品视频| 中文字幕人妻丝袜一区二区| 中文字幕人妻丝袜制服| 久久九九热精品免费| 在线av久久热| 精品一区二区三区av网在线观看| 国产在视频线精品| 亚洲三区欧美一区| 精品亚洲成国产av| 18禁国产床啪视频网站| 老司机影院毛片| 麻豆av在线久日| 国产成人av教育| www.熟女人妻精品国产| 精品卡一卡二卡四卡免费| 亚洲欧美日韩另类电影网站| 免费看a级黄色片| 日韩欧美三级三区| 久久人妻熟女aⅴ| 国产熟女午夜一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 国产男女内射视频| 国产成人av激情在线播放| 久久精品亚洲精品国产色婷小说| 中文亚洲av片在线观看爽 | 一级作爱视频免费观看| 欧美成人免费av一区二区三区 | 精品少妇一区二区三区视频日本电影| 乱人伦中国视频| 热99久久久久精品小说推荐| 天天添夜夜摸| 日韩大码丰满熟妇| 女人被狂操c到高潮| 人妻丰满熟妇av一区二区三区 | 精品亚洲成a人片在线观看| 久久中文看片网| 黄频高清免费视频| 久久香蕉激情| 色尼玛亚洲综合影院| 在线国产一区二区在线| 中文亚洲av片在线观看爽 | 国产精品美女特级片免费视频播放器 | 免费人成视频x8x8入口观看| 无人区码免费观看不卡| 搡老乐熟女国产| 不卡av一区二区三区| 久久久久久人人人人人| 亚洲精品久久成人aⅴ小说| 久久久水蜜桃国产精品网| 日韩人妻精品一区2区三区| 老司机午夜十八禁免费视频| 国产精品一区二区免费欧美| 亚洲一区中文字幕在线| 欧美+亚洲+日韩+国产| 黄片小视频在线播放| 岛国毛片在线播放| 国产成人精品久久二区二区91| 色综合欧美亚洲国产小说| 精品一区二区三区av网在线观看| 又黄又爽又免费观看的视频| 午夜免费成人在线视频| 欧美午夜高清在线| 超色免费av| 亚洲成a人片在线一区二区| 国产成人免费无遮挡视频| 成年版毛片免费区| 一进一出抽搐gif免费好疼 | 91成人精品电影| 中文字幕人妻丝袜一区二区| 精品国产乱码久久久久久男人| bbb黄色大片| 亚洲欧洲精品一区二区精品久久久| 美女福利国产在线| 大型黄色视频在线免费观看| 亚洲视频免费观看视频| 99国产综合亚洲精品| 精品无人区乱码1区二区| 最新美女视频免费是黄的| 午夜福利影视在线免费观看| 国产亚洲精品第一综合不卡| 99re6热这里在线精品视频| www.自偷自拍.com| 12—13女人毛片做爰片一| 亚洲欧美日韩高清在线视频| 黄色 视频免费看| 久久精品aⅴ一区二区三区四区| 久久久久久亚洲精品国产蜜桃av| 久久久久视频综合| 一级毛片高清免费大全| av免费在线观看网站| 精品一品国产午夜福利视频| 久久中文字幕一级| 美女 人体艺术 gogo| 亚洲欧美日韩高清在线视频| 亚洲一区中文字幕在线| 国产一区在线观看成人免费| 黑丝袜美女国产一区| 人成视频在线观看免费观看| 久久久国产成人精品二区 | 亚洲国产中文字幕在线视频| 王馨瑶露胸无遮挡在线观看| 激情视频va一区二区三区| 99国产精品免费福利视频| 欧美在线一区亚洲| 99久久精品国产亚洲精品| 99久久99久久久精品蜜桃| 超碰97精品在线观看| 久久亚洲精品不卡| 久久天躁狠狠躁夜夜2o2o| 欧美黄色淫秽网站| 99国产精品99久久久久| 嫩草影视91久久| 99re6热这里在线精品视频| 母亲3免费完整高清在线观看| 免费av中文字幕在线| 无限看片的www在线观看| 好男人电影高清在线观看| 免费高清在线观看日韩| 如日韩欧美国产精品一区二区三区| 免费久久久久久久精品成人欧美视频| 侵犯人妻中文字幕一二三四区| 国产精品久久久av美女十八| 久久国产精品大桥未久av| 18禁裸乳无遮挡动漫免费视频| 窝窝影院91人妻| 精品国产超薄肉色丝袜足j| 亚洲熟妇熟女久久| 亚洲精品自拍成人| 成年女人毛片免费观看观看9 | 999精品在线视频| 国产精品国产av在线观看| 桃红色精品国产亚洲av| 国产精品 欧美亚洲| 精品福利观看| 国产成人系列免费观看| 久久99一区二区三区| 欧美一级毛片孕妇| 好看av亚洲va欧美ⅴa在| av在线播放免费不卡| 性色av乱码一区二区三区2| 久久久久久亚洲精品国产蜜桃av| 婷婷成人精品国产| 精品一品国产午夜福利视频| a在线观看视频网站| 国产在线精品亚洲第一网站| 国产精品综合久久久久久久免费 | 天天添夜夜摸| 国产精品秋霞免费鲁丝片| 啦啦啦在线免费观看视频4| 涩涩av久久男人的天堂| 亚洲七黄色美女视频| 大码成人一级视频| 久久香蕉精品热| 他把我摸到了高潮在线观看| 精品福利永久在线观看| 两个人看的免费小视频| 一级黄色大片毛片| 国产精品二区激情视频| 欧美日韩视频精品一区| 一边摸一边抽搐一进一出视频| 国产一区二区激情短视频| 国产aⅴ精品一区二区三区波| 国产主播在线观看一区二区| 欧美日韩视频精品一区| 欧美丝袜亚洲另类 | 一本一本久久a久久精品综合妖精| 成年人黄色毛片网站| 精品国产乱子伦一区二区三区| 国产99白浆流出| av有码第一页| 久久久精品国产亚洲av高清涩受| 91精品国产国语对白视频| 日韩 欧美 亚洲 中文字幕| 国产亚洲精品久久久久5区| 久久久久国产一级毛片高清牌| 侵犯人妻中文字幕一二三四区| 久久精品国产亚洲av香蕉五月 | 免费日韩欧美在线观看| 亚洲国产精品sss在线观看 | 国产成人系列免费观看| av超薄肉色丝袜交足视频| 精品少妇久久久久久888优播| 亚洲精品国产区一区二| 黑丝袜美女国产一区| 亚洲国产欧美日韩在线播放| 欧美老熟妇乱子伦牲交| 女警被强在线播放| av免费在线观看网站| 亚洲第一青青草原| 人人妻人人澡人人看| 欧美性长视频在线观看| 国产又色又爽无遮挡免费看| 免费在线观看视频国产中文字幕亚洲| 老熟妇仑乱视频hdxx| 99国产综合亚洲精品| 桃红色精品国产亚洲av| 欧美最黄视频在线播放免费 | 亚洲精品国产一区二区精华液| 激情视频va一区二区三区| 99精品欧美一区二区三区四区| 国产麻豆69| 欧美激情久久久久久爽电影 | aaaaa片日本免费| 国产精品久久久久成人av| 亚洲精品久久成人aⅴ小说| 欧美日韩福利视频一区二区| 一级毛片高清免费大全| 久久精品熟女亚洲av麻豆精品| 大香蕉久久成人网| 亚洲色图综合在线观看| 久久国产乱子伦精品免费另类| 久久亚洲真实| 久久人人爽av亚洲精品天堂| 搡老岳熟女国产| 日韩人妻精品一区2区三区| 国产精品电影一区二区三区 | 在线观看午夜福利视频| 国产成人av激情在线播放| 午夜精品国产一区二区电影| 国产av精品麻豆| 亚洲情色 制服丝袜| 老熟女久久久| 视频在线观看一区二区三区| 99久久99久久久精品蜜桃| 国产xxxxx性猛交| 桃红色精品国产亚洲av| 美国免费a级毛片| 亚洲片人在线观看| 亚洲人成伊人成综合网2020| 欧美日韩瑟瑟在线播放| 宅男免费午夜| 久久久久国产一级毛片高清牌| 水蜜桃什么品种好| 国产精品自产拍在线观看55亚洲 | 国产精品综合久久久久久久免费 | 久久久国产成人免费| 国产精品乱码一区二三区的特点 | 亚洲免费av在线视频| 亚洲五月色婷婷综合| 国产三级黄色录像| av天堂久久9| 自线自在国产av| 人成视频在线观看免费观看| 亚洲成人手机| 精品一品国产午夜福利视频| 99热网站在线观看| 久久国产精品影院| 久久久国产精品麻豆| 久久热在线av| 十八禁高潮呻吟视频| 久久性视频一级片| 中文字幕精品免费在线观看视频| 午夜免费观看网址| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片 | 国产激情欧美一区二区| 一区福利在线观看| 女人精品久久久久毛片| 丰满饥渴人妻一区二区三| 久久亚洲真实| 国产1区2区3区精品| 久久久精品免费免费高清| 很黄的视频免费| 久热爱精品视频在线9| 性色av乱码一区二区三区2| 人人澡人人妻人| 精品国产一区二区三区四区第35| 中文字幕人妻熟女乱码| 老司机亚洲免费影院| 国产一区二区三区在线臀色熟女 | 久久这里只有精品19| 黄色成人免费大全| 国产成人精品无人区| 一级毛片女人18水好多| 欧美在线黄色| 国产成人av教育| 侵犯人妻中文字幕一二三四区| 99国产精品一区二区蜜桃av | 一区福利在线观看| 亚洲人成伊人成综合网2020| 在线观看免费视频日本深夜| av网站免费在线观看视频| 一级毛片精品| 亚洲第一青青草原| 成熟少妇高潮喷水视频| 大型黄色视频在线免费观看| 亚洲色图 男人天堂 中文字幕| 不卡av一区二区三区| 最新在线观看一区二区三区| 午夜两性在线视频| 一夜夜www| 久久精品亚洲精品国产色婷小说| 亚洲第一青青草原| 成人18禁在线播放| 色在线成人网| 在线观看日韩欧美| 国产一区二区三区视频了| 满18在线观看网站| 国产精品乱码一区二三区的特点 | 亚洲精品国产一区二区精华液| 黑人巨大精品欧美一区二区蜜桃| 99riav亚洲国产免费| 少妇 在线观看| 国产精品av久久久久免费| 丝袜美足系列| 99国产精品免费福利视频| 欧美日韩亚洲国产一区二区在线观看 | 日韩免费高清中文字幕av| 飞空精品影院首页| 国产精品久久视频播放| 午夜激情av网站| 国产有黄有色有爽视频| 热re99久久国产66热| 黄片小视频在线播放| 最近最新中文字幕大全电影3 | 亚洲 欧美一区二区三区| 80岁老熟妇乱子伦牲交| 婷婷丁香在线五月| 免费看a级黄色片| 欧美日韩一级在线毛片| 久久久久久人人人人人| 日韩大码丰满熟妇| 一边摸一边做爽爽视频免费| 亚洲精品中文字幕一二三四区| 欧美人与性动交α欧美软件| 成人免费观看视频高清| 国产成人av激情在线播放| 欧美黑人欧美精品刺激| 高清在线国产一区| 丰满的人妻完整版| 校园春色视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲免费av在线视频| 精品亚洲成国产av| 久久久精品国产亚洲av高清涩受| 丝瓜视频免费看黄片| 国产av又大| 国产亚洲一区二区精品| 国产精华一区二区三区| 亚洲专区国产一区二区| 99riav亚洲国产免费| 97人妻天天添夜夜摸| 午夜福利在线观看吧| 成年人午夜在线观看视频| av网站免费在线观看视频| 手机成人av网站| 高清黄色对白视频在线免费看| 激情在线观看视频在线高清 | 精品国产乱子伦一区二区三区| 丰满的人妻完整版| 夫妻午夜视频| 一进一出好大好爽视频| 欧美中文综合在线视频| 久久精品国产亚洲av香蕉五月 | 久久人人97超碰香蕉20202| 香蕉国产在线看| 国产精品综合久久久久久久免费 | 狂野欧美激情性xxxx| 村上凉子中文字幕在线| 中出人妻视频一区二区| 亚洲欧美色中文字幕在线| 欧美另类亚洲清纯唯美| 日本欧美视频一区|