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

    基于微衛(wèi)星標(biāo)記的不同種質(zhì)資源羅氏沼蝦遺傳多樣性研究

    2021-05-06 03:16:56許珊華章嘉淇盧婷
    江蘇農(nóng)業(yè)科學(xué) 2021年6期
    關(guān)鍵詞:遺傳多樣性種質(zhì)資源

    許珊華 章嘉淇 盧婷

    摘要:為探究江浙地區(qū)不同羅氏沼蝦(Macrobrachium rosenbergii)群體的遺傳多樣性和遺傳結(jié)構(gòu),利用16個(gè)微衛(wèi)星標(biāo)記對(duì)潮豐(CF)、藍(lán)天(LT)、源泉A(YQA)、源泉B(YQB)、嘉豐(JF)、南太湖(NTH)等6個(gè)江浙地區(qū)養(yǎng)殖群體和1個(gè)泰國(guó)正大(ZD)引進(jìn)群體進(jìn)行遺傳變異分析。結(jié)果顯示,16個(gè)微衛(wèi)星位點(diǎn)均為高度多態(tài)位點(diǎn);7個(gè)群體均為較高的遺傳多樣性,期望雜合度(He)和多態(tài)信息含量(PIC)均大于0.7,遺傳多樣性大小排序?yàn)閆D>YQA>YQB>JF>LT>NTH>CF;遺傳分化指數(shù)(Fst)分析結(jié)果顯示,泰國(guó)正大與江浙地區(qū)各群體間存在中等程度的遺傳分化,F(xiàn)st介于 0.101 45~0123 48之間;江浙地區(qū)群體除NTH與CF群體間為中等程度遺傳分化(Fst=0.050 98)外,其余群體間的遺傳分化程度較低,F(xiàn)st介于0.015 71~0.040 99之間。AMOVA分析結(jié)果顯示,遺傳變異主要發(fā)生在個(gè)體內(nèi)和群體內(nèi)個(gè)體間,群體間遺傳變異僅占2.10%。依據(jù)Neis遺傳距離構(gòu)建的非加權(quán)組平均法(UPGMA)系統(tǒng)進(jìn)化樹顯示,泰國(guó)正大群體獨(dú)占一支,江浙地區(qū)6個(gè)群體聚為另一支。Structure分析結(jié)果顯示,所有樣本被劃分為2個(gè)理論群,江浙地區(qū)6個(gè)群體為一個(gè)集群,泰國(guó)正大群體為一個(gè)集群。該研究不但揭示了江浙地區(qū)羅氏沼蝦養(yǎng)殖群體的遺傳多樣性現(xiàn)狀,而且也為羅氏沼蝦種質(zhì)資源的保護(hù)、利用以及優(yōu)良品種的選育提供了參考信息。

    關(guān)鍵詞:羅氏沼蝦;微衛(wèi)星標(biāo)記;遺傳多樣性;遺傳結(jié)構(gòu);種質(zhì)資源

    羅氏沼蝦(Macrobrachium rosenbergii)別稱馬來西亞大蝦、泰國(guó)蝦、淡水長(zhǎng)臂大蝦等,是世界上最大的淡水蝦,原產(chǎn)于東南亞地區(qū),具有重要的經(jīng)濟(jì)價(jià)值和營(yíng)養(yǎng)價(jià)值[1-2]。其食性廣、生長(zhǎng)快、個(gè)體大、殼薄體肥、肉質(zhì)鮮嫩、營(yíng)養(yǎng)豐富,因而成為廣泛養(yǎng)殖的淡水蝦,2016年全球產(chǎn)量達(dá)23.4 萬(wàn)t[3-4]。羅氏沼蝦于1976年由中國(guó)農(nóng)業(yè)科學(xué)院首次引入我國(guó)大陸[5],隨后在我國(guó)十多個(gè)省市自治區(qū)廣泛進(jìn)行養(yǎng)殖,且產(chǎn)量不斷增加,2018年產(chǎn)量達(dá)13.33萬(wàn)t[6]。目前,我國(guó)已成為全球羅氏沼蝦養(yǎng)殖量最多的國(guó)家,養(yǎng)殖產(chǎn)量占世界總產(chǎn)量的50%以上[7]。但在羅氏沼蝦產(chǎn)業(yè)迅速發(fā)展的同時(shí)也面臨諸多問題,如抗病力下降、生長(zhǎng)速度減緩、個(gè)體質(zhì)量減小、性成熟提早等,究其原因主要在于種質(zhì)資源遺傳多樣性的喪失。養(yǎng)殖的羅氏沼蝦在世代繁衍過程中,往往由于親本更新不及時(shí)和種群較小,而出現(xiàn)近親繁殖的現(xiàn)象,導(dǎo)致了其遺傳多樣性的逐漸喪失和種質(zhì)的退化[8-10]。Miao等報(bào)道了一些養(yǎng)殖物種因近交而產(chǎn)生了遺傳衰退,而這種遺傳衰退已經(jīng)是一個(gè)很嚴(yán)重的問題了,因?yàn)樗粌H影響產(chǎn)量和經(jīng)濟(jì)效益,還影響了衰退物種所在的自然生態(tài)系統(tǒng)和發(fā)展的可持續(xù)性[11]。遺傳多樣性是物種進(jìn)化依賴的基礎(chǔ),與種群的適應(yīng)能力、生存能力及進(jìn)化能力呈正相關(guān)[12]。準(zhǔn)確評(píng)價(jià)種質(zhì)資源的遺傳多樣性可以為親本選擇、后代遺傳變異和雜種優(yōu)勢(shì)預(yù)測(cè)提供預(yù)測(cè)指導(dǎo),提高育種效率[13]。因此,了解現(xiàn)有羅氏沼蝦種質(zhì)資源的遺傳多樣性和遺傳結(jié)構(gòu)對(duì)種質(zhì)資源的有效保護(hù)、高效利用,現(xiàn)有性狀的改良,新品種的選育等都至關(guān)重要。

    微衛(wèi)星(microsatellite)標(biāo)記憑借著其在基因組中分布廣泛、多態(tài)性豐富、具有共顯性、高突變率等優(yōu)點(diǎn)而被廣泛用于群體遺傳多樣性檢測(cè)和遺傳結(jié)構(gòu)的分析[14-15]。目前,已有一些采用微衛(wèi)星標(biāo)記對(duì)羅氏沼蝦不同養(yǎng)殖和野生群體進(jìn)行遺傳多樣性研究的報(bào)道,但這些報(bào)道多數(shù)是研究我國(guó)臺(tái)灣[16]和國(guó)外[17-21]羅氏沼蝦群體遺傳多樣性的,我國(guó)大陸報(bào)道的研究主要集中在廣東和廣西地區(qū)[19-22]的羅氏沼蝦養(yǎng)殖群體,而關(guān)于江蘇省和浙江省的羅氏沼蝦養(yǎng)殖群體的報(bào)道較少,僅見于文獻(xiàn)[19-21],且群體數(shù)量較少。江蘇省和浙江省是羅氏沼蝦的主要養(yǎng)殖區(qū),2018年這2個(gè)省的產(chǎn)量占全國(guó)羅氏沼蝦養(yǎng)殖總產(chǎn)量的62%,這2個(gè)省也均是羅氏沼蝦苗種的培育基地[6]。羅氏沼蝦引入我國(guó)已有44年,這期間盡管有多次小規(guī)模的重新引種,也有選育新品種(南太湖2號(hào)),但多數(shù)育苗場(chǎng)仍然缺乏科學(xué)的選育技術(shù)和保種措施[20]。因此,本研究選取江蘇和浙江地區(qū)的6個(gè)養(yǎng)殖群體和1個(gè)泰國(guó)引進(jìn)群體(作為參比群體),利用微衛(wèi)星標(biāo)記對(duì)江浙和泰國(guó)群體的遺傳變異進(jìn)行比較分析,旨在了解江浙地區(qū)現(xiàn)有羅氏沼蝦種質(zhì)資源遺傳多樣性的現(xiàn)狀,以期為我國(guó)羅氏沼蝦種質(zhì)的保護(hù)利用和優(yōu)良品種的選育提供參考信息。

    1 材料與方法

    1.1 試驗(yàn)材料與基因組DNA的提取

    試驗(yàn)樣品分別取材于不同的羅氏沼蝦養(yǎng)殖戶,其蝦苗來源于6家育苗企業(yè)(圖1)。江浙地區(qū)6個(gè)養(yǎng)殖群體分別為潮豐(CF,30尾,湖州市潮豐水產(chǎn)育苗公司)、藍(lán)天(LT,30尾,浙江藍(lán)天生態(tài)農(nóng)業(yè)開發(fā)有限公司)、源泉A(YQA,30尾,湖州源泉水產(chǎn)有限公司)、源泉B(YQB,30尾,湖州源泉水產(chǎn)有限公司)、嘉豐(JF,30尾,揚(yáng)州市嘉豐羅氏沼蝦良種繁殖有限公司)和南太湖(NTH,30尾,浙江南太湖淡水水產(chǎn)種業(yè)有限公司),1個(gè)泰國(guó)正大引進(jìn)群體(ZD,30尾,嘉興豐源農(nóng)業(yè)科技有限公司)。每尾蝦均剪取其背部肌肉組織,放置于無(wú)水乙醇中儲(chǔ)存。

    取少量樣品于1.5 mL的離心管中,盡量剪碎,于烘箱中烘干。加入HOM Buffer和蛋白酶K,搖勻后置于55 ℃搖床中消化至透明,再加入三氯甲烷抽提、離心,將上清液轉(zhuǎn)移,于-20 ℃沉淀,最后漂洗、干燥,加入雙蒸水溶解。采用1%瓊脂糖凝膠電泳檢測(cè)DNA的純度和完整性。采用紫外分光光度計(jì)檢測(cè)DNA質(zhì)量和濃度,部分稀釋至50 ng/μL,剩余DNA保存于-20 ℃冰箱中。

    1.2 微衛(wèi)星引物、PCR擴(kuò)增和PCR樣品測(cè)序

    16對(duì)引物來自文獻(xiàn)[20-22],由表1可知,引物均由武漢天一輝遠(yuǎn)生物科技有限公司合成。

    PCR擴(kuò)增體系為25 μL, 含MgCl2(25 mmol/L)1.5 μL、dNTP(2.5 mmol/L)1.0 μL、10×Buffer 2.5 μL、上下游引物各1.0 μL、ddH2O 15.9 μL、DNA 2.0 μL、rTaq 0.1 μL。PCR擴(kuò)增程序:95 ℃預(yù)變性3 min;95 ℃變性30 s,55~65 ℃退火40 s,72 ℃ 延伸30 s,35個(gè)循環(huán);最后72 ℃延伸10 min。

    合格的PCR樣品委托武漢天一輝遠(yuǎn)生物科技有限公司進(jìn)行測(cè)序分析。

    1.3 數(shù)據(jù)統(tǒng)計(jì)和分析

    利用PopGene 1.32軟件計(jì)算等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、Shannon信息指數(shù)(I)、觀測(cè)雜合度(Ho)、期望雜合度(He)、Nei氏遺傳距離、遺傳相似度、近交系數(shù)(Fis)、基因流(Nm)[23];采用Cervus 3.0軟件進(jìn)行多態(tài)信息含量(PIC)的計(jì)算[24]。根據(jù)Nei氏遺傳距離用MEGA 4.0軟件來構(gòu)建群體間聚類圖[25]。采用Arlequin 3.5軟件進(jìn)行群體間遺傳分化指數(shù)(Fst)計(jì)算和分子方差分析(AMOVA)[26]。采用Structure 231軟件分析群體遺傳結(jié)構(gòu)[27],利用在線軟件(http://clumpak.tau.ac.il/)得出最佳的理論群體數(shù)(K值),并繪制群體遺傳結(jié)構(gòu)圖。

    1.4 試驗(yàn)時(shí)間和地點(diǎn)

    試驗(yàn)于2018年7月開始,2019年12月結(jié)束。試驗(yàn)地點(diǎn)為中國(guó)水產(chǎn)科學(xué)研究院水生動(dòng)物繁育與營(yíng)養(yǎng)重點(diǎn)實(shí)驗(yàn)室。

    2 結(jié)果與分析

    2.1 微衛(wèi)星位點(diǎn)的多態(tài)性分析

    由微衛(wèi)星位點(diǎn)多樣性統(tǒng)計(jì)結(jié)果(表2)可知,16個(gè)位點(diǎn)共檢測(cè)出323個(gè)等位基因(Na),Na為10~29個(gè),每個(gè)位點(diǎn)平均有20個(gè)。Ne為2.611 5~12431 6個(gè),平均值為7.939 3個(gè)。Na與Ne數(shù)值相差較大,說明群體中等位基因分布不均勻。I為1376 0~2.738 2,平均值為2.276 5,I值較大,說明群體的遺傳多樣性高。Ho為0.329 7~1.000 0,平均值為0650 2。He為0.618 6~0.921 8,平均值為 0.834 7。16個(gè)位點(diǎn)中有13個(gè)位點(diǎn)的Ho低于He,說明群體內(nèi)自交率較高。PIC為0.539 7~0.872 4,平均值為0.776 1。根據(jù)Botstein等劃分標(biāo)準(zhǔn),16個(gè)位點(diǎn)的PIC均高于0.5,為高度多態(tài)性[28]。F檢驗(yàn)數(shù)據(jù)顯示,有4個(gè)位點(diǎn)的Fis為負(fù)值,其余12個(gè)為正值,說明近交程度高。根據(jù)Wright建議,有6個(gè)位點(diǎn)的遺傳分化很?。‵st<0.05),有10個(gè)位點(diǎn)的遺傳分化為中等(0.05 2.2 羅氏沼蝦7個(gè)群體多態(tài)性分析? ? 7個(gè)群體的遺傳多樣性。由表3可知,7個(gè)群體的Na介于9.625 0~12.187 5之間,Ne介于5.078 1~7261 2之間, I介于1.770 6~2.148 4之間,Ho介于0.621 4~0.702 8之間,He介于0.764 4~0.859 6 之間,PIC介于0.744 3~0.825 8之間。泰國(guó)正大群體的遺傳多樣性參數(shù)均高于江浙地區(qū)的6個(gè)羅氏沼蝦群體,CF群體的各參數(shù)均為最小。Ne、I和PIC在7個(gè)群體內(nèi)趨勢(shì)一致,表現(xiàn)為ZD>YQA>YQB> JF>LT>NTH>CF。 7個(gè)群體的多態(tài)信息含量均高于0.7,說明該研究的7個(gè)群體有豐富的遺傳多樣性。

    2.3 群體間遺傳分化與遺傳距離的分析

    由表4可知,群體間遺傳分化,泰國(guó)正大與所有群體間的Fst介于0.101 45~0.123 48之間,為中等程度的遺傳分化(0.05? ?由表6可知,7個(gè)羅氏沼蝦群體的遺傳距離介于0.095 0~1.027 4之間,CF和ZD群體遺傳距離最遠(yuǎn),為1.027 4,遺傳相似度最低,為0.357 9。LT和CF群體遺傳距離最小,為0.095 0,遺傳相似度最大,為0.909 4。由圖2可知,為根據(jù)Neis遺傳距離構(gòu)建的非加權(quán)組平均法(UPGMA)聚類樹,此樹分為2支:泰國(guó)正大群體獨(dú)占一支,江浙地區(qū)的群體聚為一支。江浙地區(qū)的群體這一支又分為2支:NTH群體為一支,其余群體為一支(此支又分為2支:CF與LT聚為一支,YQB、JF和YQA聚為一支)。

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

    將Burnin運(yùn)算長(zhǎng)度設(shè)置為100 000,K值預(yù)設(shè)為1~7,每個(gè)K值重復(fù)10次,在K=2時(shí),Delta K最大,即所有參試個(gè)體最佳分組為2個(gè)理論群。由圖3可知, 在K=2的情況下, 江浙地區(qū)的6個(gè)群體為一個(gè)集群,泰國(guó)正大群體為一個(gè)集群。在K=3的情況下,CF、LT、YQB、YQA和JF聚為一群,NTH和ZD分別獨(dú)立成群。在K=4的情況下,CF和LT為一個(gè)集群,YQB、JF和YQA為一個(gè)集群,NTH和ZD分別獨(dú)立成群。

    3 討論

    3.1 微衛(wèi)星位點(diǎn)的遺傳多樣性

    多態(tài)信息含量可體現(xiàn)微衛(wèi)星位點(diǎn)的遺傳變異程度,數(shù)值越大,遺傳變異程度越高,即多樣性越高。該研究所用的16個(gè)微衛(wèi)星位點(diǎn)均為高度多態(tài)性位點(diǎn)(PIC>0.5),每個(gè)位點(diǎn)平均等位基因數(shù)(Na)為20.19個(gè),與孫成飛等的研究結(jié)果[20,22,30-31]相比,均高于這些研究的平均等位基因數(shù),說明這些位點(diǎn)可提供豐富的遺傳信息,可有效地用于羅氏沼蝦群體遺傳變異分析。

    3.2 群體遺傳多樣性分析

    等位基因數(shù)和期望雜合度是衡量群體遺傳多樣性常用的2個(gè)參數(shù)[32]。本研究7個(gè)群體的Na介于9.625 0~12.187 5之間,He介于0.764 4~0859 6之間。Chareontawee等研究了臺(tái)灣5個(gè)養(yǎng)殖群體2個(gè)野生群體的遺傳多樣性,認(rèn)為養(yǎng)殖群體和野生群體相似,且所有群體都顯示出了相對(duì)較高的遺傳變異,Na介于7.50~10.67之間,He介于064~0.73之間[16]。與該報(bào)道相比,本研究中Na和He都較高。Schneider等研究了美國(guó)、印度、以色列和緬甸4個(gè)國(guó)家7個(gè)養(yǎng)殖群體和2個(gè)野生群體的遺傳多樣性,緬甸野生群體的Na和He均為最高,印度野生群體的Na和He介于7個(gè)養(yǎng)殖群體值之間,9個(gè)群體的Na介于3.96~20.45之間,He介于0.579 4~0935 6之間[17]。與該報(bào)道相比,本研究的Na和Ne都處于中等水平。Nguyen Thanh等研究了中國(guó)6個(gè)養(yǎng)殖群體和越南2個(gè)野生群體的遺傳多樣性,其中1個(gè)野生群體的Na和He均為最高,另1個(gè)野生群體的Na和He介于6個(gè)養(yǎng)殖群體值之間,8個(gè)群體的Na介于9.666~12.000之間,He介于0.785~0.835之間[19]。與該報(bào)道相比,本研究的Na、He與其相似。Khan等研究了孟加拉國(guó)3個(gè)野生群體的遺傳多樣性,3個(gè)群體的Na介于9.28~957之間,He介于0.804~0.827之間[18]。與該報(bào)道相比,本研究的Na高于其Na、He??梢姡c同類研究相比,本研究7個(gè)養(yǎng)殖群體的Na和He相對(duì)處于中等偏高的水平,與一些野生群體相當(dāng)。當(dāng)雜合度在0.5~0.8之間即可認(rèn)為該群體具有較高的多樣性[33],本研究7個(gè)群體的期望雜合度均在0.7以上,故屬較高的遺傳多樣性水平。7個(gè)群體的多態(tài)信息含量介于0.744 3~0.825 8之間,均高于0.5,也說明5個(gè)群體有豐富的遺傳多樣性。但是7個(gè)群體的觀測(cè)雜合度均低于期望雜合度,表明7個(gè)群體都存在一定程度的近交現(xiàn)象,因而為了避免群體因進(jìn)一步近交而出現(xiàn)衰退現(xiàn)象,要及時(shí)引入外來的優(yōu)勢(shì)群體、野生群體或選擇雜合子較多的近緣群體進(jìn)行育種。

    3.3 群體遺傳分化與遺傳結(jié)構(gòu)分析

    基因流(Nm)是不同群體間由于個(gè)體的交換而發(fā)生基因交流的過程[34-35],它最基本的作用是削弱群體間的遺傳差異。Wright認(rèn)為Nm大于1可抑制種群間的分化,本研究16個(gè)位點(diǎn)的Nm均大于1,均值高達(dá)3.792 1,表明群體間有較高的基因交流,遺傳分化程度較低[36]。遺傳分化指數(shù)分析結(jié)果表明,除泰國(guó)正大群體與江浙地區(qū)所有群體間存在中等程度的遺傳分化(Fst在0.101 45~0.123 48之間)外,江浙地區(qū)6個(gè)群體間遺傳分化很小,即使遺傳距離最遠(yuǎn)的NTH和CF群體間Fst也僅為0.050 98,剛超過0.05。采用Evanno等人的方法進(jìn)行分析計(jì)算,得到最佳K值為2,即所有樣本被劃分為2個(gè)理論群[37]。該Structure分析結(jié)果與根據(jù)Neis遺傳距離構(gòu)建的UPGMA系統(tǒng)樹相一致,均直觀表明江浙地區(qū)的6個(gè)群體為一個(gè)集群,泰國(guó)正大群體獨(dú)自為一個(gè)集群。所以江浙地區(qū)的6個(gè)群體從遺傳上講實(shí)為一個(gè)群體,可能是這6個(gè)群體所在的育苗公司地理位置較近,平時(shí)存在親蝦或苗種的互相供給,所以基因交流較多,遺傳差異較小。推測(cè)它們的親本來源均與選育新品種南太湖2號(hào)有關(guān)。系統(tǒng)聚類樹分支長(zhǎng)度表明江浙的養(yǎng)殖群體和引進(jìn)的泰國(guó)群體間親緣關(guān)系相對(duì)較遠(yuǎn),此結(jié)果與孫成飛等的研究結(jié)果[20]一致,據(jù)此可以考慮引入泰國(guó)群體與江浙地區(qū)的群體進(jìn)行雜交育種,以獲得雜種優(yōu)勢(shì)。

    綜上,江浙地區(qū)6個(gè)羅氏沼蝦養(yǎng)殖群體存在一定程度的近交,群體間遺傳分化很小,均擁有豐富的遺傳多態(tài)性,具有良好的選育潛力和前景,但須及時(shí)引種。與泰國(guó)正大群體相比,江浙地區(qū)養(yǎng)殖群體的遺傳多態(tài)性略低。同時(shí),ZD群體與江浙地區(qū)的6個(gè)羅氏沼蝦群體的親緣關(guān)系較遠(yuǎn),存在中等程度的遺傳分化,因此建議在今后育種時(shí),可考慮引入正大的羅氏沼蝦個(gè)體,以增加江浙地區(qū)種群的遺傳多態(tài)性,避免種質(zhì)衰退。該研究不但揭示了江浙一帶羅氏沼蝦種質(zhì)資源的遺傳變異現(xiàn)狀,而且也為羅氏沼蝦育種者和生產(chǎn)者在羅氏沼蝦種質(zhì)資源保護(hù)、優(yōu)良品種選育以及生產(chǎn)實(shí)踐上提供參考。

    參考文獻(xiàn):

    [1]New M B,Singholka S. Freshwater prawn farming. A manual for the culture of Macrobrachium rosenbergii[M]. FAO Fisheries Technical Paper,1982,225:116.

    [2]謝忠明,李增崇,趙明森. 淡水經(jīng)濟(jì)蝦類養(yǎng)殖技術(shù)[M]. 北京:中國(guó)農(nóng)業(yè)出版社,2002.

    [3]孫麗慧,潘 茜,陳雪峰,等. 羅氏沼蝦高效飼料轉(zhuǎn)化家系與原種后代肌肉營(yíng)養(yǎng)成分的比較[J]. 飼料工業(yè),2019,40(4):40-45.

    [4]FAO. The state of world fisheries and aquaculture 2018:meeting the sustainable development goals[R]. Rome:FAO,2018:23.

    [5]Yang G L,F(xiàn)rinsko M,Chen X F,et al. Current status of the giant freshwater prawn (Macrobrachium rosenbergii) industry in China,with special reference to live transportation[J]. Aquaculture Research,2012,43(7):1049-1055.

    [6]農(nóng)業(yè)農(nóng)村部漁業(yè)漁政管理局,全國(guó)水產(chǎn)技術(shù)推廣總站,中國(guó)水產(chǎn)學(xué)會(huì).中國(guó)漁業(yè)統(tǒng)計(jì)年鑒[M]. 北京:中國(guó)農(nóng)業(yè)出版社,2019:1-152.

    [7]唐瓊英,夏正龍,蔡繆熒,等. 羅氏沼蝦養(yǎng)殖群體表型性狀間的相關(guān)性及類群差異分析[J]. 中國(guó)水產(chǎn)科學(xué),2019,26(6):1075-1085.

    [8]呂 敏,黃光華,李 旻,等. 異型雄性羅氏沼蝦遺傳多樣性的微衛(wèi)星分析[J]. 水產(chǎn)科學(xué),2019,38(3):355-360.

    [9]韓耀全. 羅氏沼蝦迅速衰落的主要原因[C]//廣西水產(chǎn)研究所論文集,2006:196-198.

    [10]肖楚康,方 劉,阮國(guó)良,等. 羅氏沼蝦淡化養(yǎng)殖的現(xiàn)狀與展望[J]. 江蘇農(nóng)業(yè)科學(xué),2019,47(8):29-33.

    [11]Miao W M,Ge X P. Freshwater prawn culture in China:an overview[J]. Aquaculture Asia,2002,7(1):9-12.

    [12]Vandewoestijne S,Schtickzelle N,Baguette M. Positive correlation between genetic diversity and fitness in a large,well-connected metapopulation[J]. BMC Biology,2008,6:46.

    [13]Liu F,Qu Y K,Geng C,et al. Analysis of the population structure and genetic diversity of the red swamp crayfish (Procambarus clarkia) in China using SSR markers[J]. Electronic Journal of Biotechnology,2020,47:59-71.

    [14]楊子拓,劉 麗. 微衛(wèi)星分子標(biāo)記技術(shù)在水產(chǎn)養(yǎng)殖行業(yè)的應(yīng)用[J]. 廣東農(nóng)業(yè)科學(xué),2014,41(18):114-116,126.

    [15]陳萬(wàn)光,賈曉慧. 微衛(wèi)星標(biāo)記在水產(chǎn)動(dòng)物遺傳研究方面的應(yīng)用[J]. 安徽農(nóng)業(yè)科學(xué),2009,37(12):5659-5660.

    [16]Chareontawee K,Poompuang S,Na-Nakorn U,et al. Genetic diversity of hatchery stocks of giant freshwater prawn (Macrobrachium rosenbergii) in Thailand[J]. Aquaculture,2007,271(1/2/3/4):121-129.

    [17]Schneider K J,Tidwell J H,Gomelsky B,et al. Genetic diversity of cultured and wild populations of the giant freshwater prawn Macrobrachium rosenbergii(de Man,1879)based on microsatellite analysis[J]. Aquaculture Research,2012,44(9):1-13.

    [18]Khan S R,Akter H,Sultana N,et al. Genetic diversity in three river populations of the giant freshwater prawn (Macrobrachium rosenbergii) in Bangladesh assessed by microsatellite DNA markers[J]. International Journal of Agriculture&Biology,2014,16(1):195-200.

    [19]Nguyen Thanh H,Liu Q G,Zhao L J,et al. Genetic diversity of the cultured giant freshwater prawn (Macrobrachium rosenbergii) in China based on microsatellite markers[J]. Biochemical Systematics and Ecology,2015,59:144-154.

    [20]孫成飛,葉 星,董浚鍵,等. 羅氏沼蝦6個(gè)養(yǎng)殖群體遺傳多樣性的微衛(wèi)星分析[J]. 南方水產(chǎn)科學(xué),2015,11(2):20-26.

    [21]陳佳毅. 羅氏沼蝦親蝦群體遺傳結(jié)構(gòu)分析及育苗參數(shù)比較[D]. 揚(yáng)州:揚(yáng)州大學(xué),2016:12-13.

    [22]鐘丹丹,林 勇,賓石玉,等. 兩個(gè)羅氏沼蝦種群的遺傳多樣性研究[J]. 廣東農(nóng)業(yè)科學(xué),2015,42(24):140-145.

    [23]Yeh F C,Yang R C,Boyle T J,et al. PopGene32,microsoft windows-based freeware for population genetic analysis.Version 1.32[Z]. 2000.

    [24]Kalinowski S T,Taper M L,Marshall T C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment[J]. Molecular Ecology,2007,16(5):1099-1106.

    [25]Tamura K,Dudley J,Nei M,et al. MEGA4:molecular evolutionary genetics analysis (MEGA) software version 4.0[J]. Molecular Biology and Evolution,2007,24(8):1596-1599.

    [26]Excoffier L,Lischer H E. Arlequin suite ver 3.5:a new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular Ecology Resources,2010,10(3):564-567.

    [27]Pritchard J K,Stephens M,Donnelly P. Inference of population structure using multilocus genotype data[J]. Genetics,2000,155(2):945-959.

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

    [29]Wright S. The interpretation of population structure by F-statistics with special regard to systems of mating[J]. Evolution,1965,19(3):395-420.

    [30]朱其建,戴習(xí)林,鄒衛(wèi)麗,等. 羅氏沼蝦抗病選育群體的抗病性能及其遺傳多樣性分析[J]. 水產(chǎn)學(xué)報(bào),2013,37(10):1468-1478.

    [31]蔣 飛,戴習(xí)林,朱其建,等. 羅氏沼蝦5個(gè)專門化品系選擇系生長(zhǎng)比較及其遺傳結(jié)構(gòu)分析[J]. 上海海洋大學(xué)學(xué)報(bào),2014,23(3):329-337.

    [32]Kallinowski S T. HP-RARE 1.0:a computer program for performing rarefaction on measures of allelic richness[J]. Molecular Ecology Notes,2005,5(1):187-189.

    [33]張 智,俞 丹,劉 飛,等. 西昌華吸鰍的微衛(wèi)星引物篩選及赤水河四個(gè)地理種群的 遺傳多樣性分析[J]. 水生生物學(xué)報(bào),2019,43(6):1224-1230.

    [34]Husband B C,Barrett S H. Estimates of gene flow in Eichhornia paniculate(Pontederiaceae):effects of range substructure[J]. Heredity,1995,75(6):549-560.

    [35]曲若竹,侯 林,呂紅麗,等. 群體遺傳結(jié)構(gòu)中的基因流[J]. 遺傳,2004,26(3):377-382.

    [36]Wright S. Evolution in mendelian populations[J]. Bulletin of Mathematical Biology,1990,52 (1/2):241-295.

    [37]Evanno G,Regnaut S,Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE:a simulation study[J]. Molecular Ecology,2005,14(8):2611-2620.胡振陽(yáng),程 宏,盧 臣,等. 施氮量和植物生長(zhǎng)調(diào)節(jié)劑對(duì)優(yōu)質(zhì)稻抗倒能力及產(chǎn)量的調(diào)控效應(yīng)[J]. 江蘇農(nóng)業(yè)科學(xué),2021,49(6):52-60.

    猜你喜歡
    遺傳多樣性種質(zhì)資源
    云南德宏羊奶果資源收集和評(píng)價(jià)簡(jiǎn)報(bào)
    綠肥作物紫云英研究進(jìn)展
    大白菜種質(zhì)資源抗根腫病基因CRa和CRb的分子標(biāo)記鑒定與分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    玉米種質(zhì)資源抗旱性鑒定研究進(jìn)展
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    楊梅種質(zhì)資源遺傳多樣性研究進(jìn)展
    山西大豆自然群體遺傳多樣性的研究
    看片在线看免费视频| www日本在线高清视频| 国产精品乱码一区二三区的特点| 国产毛片a区久久久久| 国产精品一区二区精品视频观看| 日本成人三级电影网站| 麻豆久久精品国产亚洲av| 99国产极品粉嫩在线观看| 国产一区在线观看成人免费| 男女做爰动态图高潮gif福利片| 国产欧美日韩精品亚洲av| 久久中文字幕人妻熟女| 最近最新免费中文字幕在线| 麻豆成人午夜福利视频| 无人区码免费观看不卡| 少妇的丰满在线观看| 在线播放国产精品三级| 亚洲av中文字字幕乱码综合| 少妇熟女aⅴ在线视频| 老汉色∧v一级毛片| 国产三级在线视频| 午夜免费观看网址| 久久久久久久午夜电影| 91在线精品国自产拍蜜月 | 最近在线观看免费完整版| 欧美丝袜亚洲另类 | 夜夜夜夜夜久久久久| 欧美中文综合在线视频| 日日摸夜夜添夜夜添小说| 欧美精品啪啪一区二区三区| 不卡一级毛片| 成人三级黄色视频| 国产精品一区二区三区四区久久| 久久久精品欧美日韩精品| 欧美成人性av电影在线观看| 亚洲成av人片免费观看| 欧美丝袜亚洲另类 | 国产午夜福利久久久久久| 国产爱豆传媒在线观看| 女人高潮潮喷娇喘18禁视频| 国产私拍福利视频在线观看| 视频区欧美日本亚洲| 啪啪无遮挡十八禁网站| 在线观看免费午夜福利视频| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| 国产精品野战在线观看| 757午夜福利合集在线观看| 精品国产亚洲在线| 午夜福利18| 亚洲18禁久久av| 欧美zozozo另类| 欧美日韩福利视频一区二区| 成人亚洲精品av一区二区| 亚洲av中文字字幕乱码综合| 亚洲专区国产一区二区| 国产精品影院久久| 色吧在线观看| 国产精品精品国产色婷婷| 久久天躁狠狠躁夜夜2o2o| 免费在线观看亚洲国产| 国产在线精品亚洲第一网站| 欧美激情久久久久久爽电影| 国产午夜福利久久久久久| 国产一区在线观看成人免费| 国产三级黄色录像| 母亲3免费完整高清在线观看| 非洲黑人性xxxx精品又粗又长| 观看免费一级毛片| 熟妇人妻久久中文字幕3abv| 亚洲一区二区三区色噜噜| 精华霜和精华液先用哪个| 亚洲av成人av| 国产高清有码在线观看视频| av女优亚洲男人天堂 | 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区在线av高清观看| 久久精品91无色码中文字幕| e午夜精品久久久久久久| 国内精品久久久久精免费| 欧美日韩瑟瑟在线播放| 日韩免费av在线播放| 久久婷婷人人爽人人干人人爱| 精品久久久久久久毛片微露脸| 久久九九热精品免费| 欧美中文综合在线视频| 91在线观看av| 午夜福利在线观看免费完整高清在 | 亚洲精品粉嫩美女一区| 中国美女看黄片| 色视频www国产| 黑人欧美特级aaaaaa片| 日本熟妇午夜| ponron亚洲| 日韩欧美一区二区三区在线观看| 中文字幕熟女人妻在线| 后天国语完整版免费观看| 欧美极品一区二区三区四区| 看黄色毛片网站| 精品国产亚洲在线| 国内揄拍国产精品人妻在线| 听说在线观看完整版免费高清| 成年人黄色毛片网站| 精品99又大又爽又粗少妇毛片 | 日韩欧美国产一区二区入口| 欧美又色又爽又黄视频| 日韩成人在线观看一区二区三区| 精华霜和精华液先用哪个| 色吧在线观看| 夜夜爽天天搞| 欧美黄色淫秽网站| 亚洲 欧美一区二区三区| 黑人巨大精品欧美一区二区mp4| av国产免费在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲色图 男人天堂 中文字幕| 高清在线国产一区| 一个人看视频在线观看www免费 | 无遮挡黄片免费观看| 天堂动漫精品| 久久精品91无色码中文字幕| 久久久久国产一级毛片高清牌| 精品久久久久久久人妻蜜臀av| 精品国产亚洲在线| 在线观看一区二区三区| 99国产综合亚洲精品| 中文字幕精品亚洲无线码一区| 怎么达到女性高潮| 99视频精品全部免费 在线 | www.精华液| 曰老女人黄片| 婷婷精品国产亚洲av在线| 久久这里只有精品中国| 日韩欧美一区二区三区在线观看| 91在线精品国自产拍蜜月 | 丰满人妻熟妇乱又伦精品不卡| 99国产精品一区二区三区| 久久人人精品亚洲av| 麻豆国产97在线/欧美| 久久久色成人| 国产aⅴ精品一区二区三区波| 成人国产综合亚洲| 人人妻,人人澡人人爽秒播| 日韩三级视频一区二区三区| 99国产精品一区二区蜜桃av| 国产麻豆成人av免费视频| 最新在线观看一区二区三区| 美女大奶头视频| 不卡一级毛片| 午夜精品一区二区三区免费看| 欧美精品啪啪一区二区三区| 2021天堂中文幕一二区在线观| 五月玫瑰六月丁香| 成人无遮挡网站| 色综合欧美亚洲国产小说| 欧美乱码精品一区二区三区| 久久久久久九九精品二区国产| 变态另类丝袜制服| 好男人在线观看高清免费视频| 国产又黄又爽又无遮挡在线| 免费av不卡在线播放| avwww免费| а√天堂www在线а√下载| av片东京热男人的天堂| 亚洲av中文字字幕乱码综合| 伊人久久大香线蕉亚洲五| 欧美日韩中文字幕国产精品一区二区三区| 1024手机看黄色片| 黑人欧美特级aaaaaa片| 亚洲人成网站在线播放欧美日韩| 午夜a级毛片| 母亲3免费完整高清在线观看| 成人亚洲精品av一区二区| 日韩精品中文字幕看吧| 啦啦啦免费观看视频1| 欧美成人性av电影在线观看| 人妻久久中文字幕网| 亚洲最大成人中文| 精华霜和精华液先用哪个| 最近最新中文字幕大全免费视频| 亚洲av电影不卡..在线观看| 欧美乱色亚洲激情| 无遮挡黄片免费观看| 亚洲av第一区精品v没综合| 国产单亲对白刺激| 两性午夜刺激爽爽歪歪视频在线观看| 国产野战对白在线观看| 我要搜黄色片| 成人鲁丝片一二三区免费| 国产伦一二天堂av在线观看| 精品人妻1区二区| 噜噜噜噜噜久久久久久91| 久久这里只有精品19| 欧美日韩国产亚洲二区| 男女午夜视频在线观看| 午夜福利免费观看在线| 非洲黑人性xxxx精品又粗又长| 欧美高清成人免费视频www| 国产三级在线视频| 少妇裸体淫交视频免费看高清| 成人鲁丝片一二三区免费| 两性午夜刺激爽爽歪歪视频在线观看| 搡老岳熟女国产| 亚洲av熟女| 人妻久久中文字幕网| 成年版毛片免费区| 九九久久精品国产亚洲av麻豆 | 亚洲成a人片在线一区二区| 亚洲一区二区三区不卡视频| 欧美在线黄色| 免费观看人在逋| www国产在线视频色| 亚洲最大成人中文| 色精品久久人妻99蜜桃| 老司机午夜十八禁免费视频| 一区二区三区国产精品乱码| www.自偷自拍.com| 久久久精品欧美日韩精品| 桃色一区二区三区在线观看| 欧美日韩国产亚洲二区| 国产三级中文精品| 日韩欧美免费精品| 欧美中文日本在线观看视频| 国产亚洲精品av在线| 日本与韩国留学比较| 观看美女的网站| 欧美成人免费av一区二区三区| 日韩国内少妇激情av| 最新中文字幕久久久久 | xxx96com| 又粗又爽又猛毛片免费看| 亚洲av电影在线进入| 欧美黄色淫秽网站| 亚洲色图av天堂| 精品国产美女av久久久久小说| 美女cb高潮喷水在线观看 | 琪琪午夜伦伦电影理论片6080| 中文字幕精品亚洲无线码一区| 国产伦在线观看视频一区| 中文字幕人成人乱码亚洲影| 99久久无色码亚洲精品果冻| 一卡2卡三卡四卡精品乱码亚洲| 女同久久另类99精品国产91| 日韩欧美一区二区三区在线观看| 两性夫妻黄色片| 看片在线看免费视频| 色在线成人网| 亚洲国产欧美一区二区综合| 91老司机精品| 一进一出抽搐gif免费好疼| 亚洲熟妇熟女久久| 亚洲欧美激情综合另类| 国产精品乱码一区二三区的特点| 亚洲国产高清在线一区二区三| 最新中文字幕久久久久 | 欧美大码av| 每晚都被弄得嗷嗷叫到高潮| 99在线人妻在线中文字幕| 亚洲人成伊人成综合网2020| 国产成+人综合+亚洲专区| av国产免费在线观看| cao死你这个sao货| 黑人欧美特级aaaaaa片| ponron亚洲| 99国产综合亚洲精品| 亚洲一区高清亚洲精品| 不卡av一区二区三区| 久久这里只有精品中国| www.www免费av| 怎么达到女性高潮| 欧美+亚洲+日韩+国产| 俄罗斯特黄特色一大片| 成人三级做爰电影| 国语自产精品视频在线第100页| 国产主播在线观看一区二区| 青草久久国产| 精品久久久久久久人妻蜜臀av| 国产精品久久久av美女十八| 亚洲aⅴ乱码一区二区在线播放| 国产av一区在线观看免费| 两个人的视频大全免费| 日韩欧美免费精品| 伦理电影免费视频| 一区二区三区高清视频在线| 91av网站免费观看| 国产精品永久免费网站| 99久久久亚洲精品蜜臀av| 国产高清有码在线观看视频| 亚洲天堂国产精品一区在线| 亚洲专区中文字幕在线| 一个人免费在线观看电影 | 91麻豆精品激情在线观看国产| 成人性生交大片免费视频hd| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久人人做人人爽| 狂野欧美白嫩少妇大欣赏| 又黄又爽又免费观看的视频| 美女高潮喷水抽搐中文字幕| 99久久国产精品久久久| 亚洲欧美日韩东京热| 亚洲av美国av| 美女 人体艺术 gogo| 欧美日韩瑟瑟在线播放| 亚洲av成人av| 国产成人欧美在线观看| 亚洲成人免费电影在线观看| 成年女人永久免费观看视频| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久久电影 | 夜夜爽天天搞| 免费在线观看日本一区| 日韩成人在线观看一区二区三区| 日韩高清综合在线| 麻豆一二三区av精品| 18禁黄网站禁片免费观看直播| 美女cb高潮喷水在线观看 | 黄频高清免费视频| 国产高清videossex| 日本免费一区二区三区高清不卡| 亚洲欧美一区二区三区黑人| 国产 一区 欧美 日韩| 婷婷精品国产亚洲av在线| 色精品久久人妻99蜜桃| 99久久精品一区二区三区| 啦啦啦免费观看视频1| 香蕉国产在线看| 成人永久免费在线观看视频| 一进一出好大好爽视频| 黄色片一级片一级黄色片| 成人国产综合亚洲| 国产黄a三级三级三级人| 色综合欧美亚洲国产小说| 免费在线观看影片大全网站| 欧美日韩乱码在线| 亚洲av成人精品一区久久| 日韩成人在线观看一区二区三区| 国产精华一区二区三区| 一个人看视频在线观看www免费 | 一夜夜www| 黑人欧美特级aaaaaa片| 日韩成人在线观看一区二区三区| 精品欧美国产一区二区三| 亚洲专区字幕在线| 国产精品一区二区三区四区久久| 99久国产av精品| 免费在线观看影片大全网站| 亚洲欧美一区二区三区黑人| 怎么达到女性高潮| 色噜噜av男人的天堂激情| 99在线视频只有这里精品首页| 亚洲欧美日韩东京热| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式 | 男女下面进入的视频免费午夜| 一级作爱视频免费观看| 成年版毛片免费区| a级毛片a级免费在线| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 桃红色精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 老汉色av国产亚洲站长工具| 成人三级做爰电影| 国产乱人视频| 亚洲国产欧美一区二区综合| 亚洲国产中文字幕在线视频| 一区二区三区激情视频| 国产精品自产拍在线观看55亚洲| 真人一进一出gif抽搐免费| 亚洲va日本ⅴa欧美va伊人久久| 天堂√8在线中文| 亚洲精品乱码久久久v下载方式 | 国产私拍福利视频在线观看| 性欧美人与动物交配| 12—13女人毛片做爰片一| 日本 欧美在线| 欧美高清成人免费视频www| 岛国视频午夜一区免费看| 91麻豆av在线| 又粗又爽又猛毛片免费看| 这个男人来自地球电影免费观看| 日韩欧美在线乱码| 88av欧美| 法律面前人人平等表现在哪些方面| 国产精品影院久久| 久久午夜亚洲精品久久| 精品福利观看| 免费大片18禁| 国产探花在线观看一区二区| 国产一区在线观看成人免费| 国产成人一区二区三区免费视频网站| 在线观看免费视频日本深夜| 日本免费a在线| 一级毛片高清免费大全| 亚洲av美国av| www.精华液| 少妇人妻一区二区三区视频| 在线观看日韩欧美| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 国产av麻豆久久久久久久| 91麻豆精品激情在线观看国产| 午夜两性在线视频| 久久久久久久精品吃奶| 搡老熟女国产l中国老女人| 亚洲熟妇熟女久久| 日韩欧美国产在线观看| 欧美在线黄色| 免费在线观看亚洲国产| 我的老师免费观看完整版| 亚洲精华国产精华精| 国产精品亚洲美女久久久| 国产黄a三级三级三级人| 欧美大码av| 国内揄拍国产精品人妻在线| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩东京热| 亚洲九九香蕉| 看片在线看免费视频| 天堂动漫精品| 欧美大码av| 亚洲美女视频黄频| 精品无人区乱码1区二区| 老汉色∧v一级毛片| 亚洲精品久久国产高清桃花| 在线十欧美十亚洲十日本专区| 久久久精品欧美日韩精品| 99精品在免费线老司机午夜| 97超视频在线观看视频| 美女扒开内裤让男人捅视频| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| 精品电影一区二区在线| 在线永久观看黄色视频| 亚洲成av人片免费观看| 精品福利观看| 性色av乱码一区二区三区2| 一区二区三区国产精品乱码| 一二三四社区在线视频社区8| 国产 一区 欧美 日韩| 男女做爰动态图高潮gif福利片| 欧美成人免费av一区二区三区| 久久久国产欧美日韩av| 国产美女午夜福利| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| 变态另类成人亚洲欧美熟女| 99久久无色码亚洲精品果冻| 成人高潮视频无遮挡免费网站| 国产伦人伦偷精品视频| 黑人操中国人逼视频| 国产三级黄色录像| 窝窝影院91人妻| 一二三四社区在线视频社区8| 亚洲av成人精品一区久久| 国产精品免费一区二区三区在线| 国产高清videossex| 男女做爰动态图高潮gif福利片| 精品99又大又爽又粗少妇毛片 | 亚洲五月天丁香| 亚洲色图av天堂| 欧美在线一区亚洲| 激情在线观看视频在线高清| 欧美成人性av电影在线观看| 免费在线观看视频国产中文字幕亚洲| 舔av片在线| 亚洲国产欧美网| or卡值多少钱| 国产亚洲精品综合一区在线观看| 亚洲国产欧美人成| 99久久精品一区二区三区| 成年女人永久免费观看视频| 国产精品野战在线观看| 久久香蕉国产精品| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 啦啦啦韩国在线观看视频| 天堂√8在线中文| 免费观看精品视频网站| 深夜精品福利| 国产高清三级在线| 丰满人妻一区二区三区视频av | 久久久国产成人精品二区| 亚洲自拍偷在线| 国产成人欧美在线观看| 国产成人系列免费观看| 91麻豆av在线| 国产不卡一卡二| 国产成人啪精品午夜网站| 国产高清激情床上av| 欧美国产日韩亚洲一区| 婷婷六月久久综合丁香| 国产av一区在线观看免费| 午夜免费观看网址| 国产综合懂色| 99精品欧美一区二区三区四区| 欧美一级a爱片免费观看看| АⅤ资源中文在线天堂| 成人欧美大片| 欧美性猛交黑人性爽| 91在线观看av| 91av网站免费观看| 日韩国内少妇激情av| 99国产精品一区二区蜜桃av| 欧美日韩精品网址| 嫩草影院入口| 中文亚洲av片在线观看爽| 两性夫妻黄色片| 成人av一区二区三区在线看| 99精品久久久久人妻精品| 国产精品久久久久久亚洲av鲁大| 亚洲成a人片在线一区二区| 欧美av亚洲av综合av国产av| 特大巨黑吊av在线直播| 国产精品av视频在线免费观看| 久久久久精品国产欧美久久久| 91麻豆精品激情在线观看国产| 亚洲天堂国产精品一区在线| 97超视频在线观看视频| 精品国产乱子伦一区二区三区| 97超视频在线观看视频| 中文字幕久久专区| 久久久国产精品麻豆| 9191精品国产免费久久| 精品国产乱码久久久久久男人| 亚洲欧美日韩东京热| 老鸭窝网址在线观看| 中文亚洲av片在线观看爽| 亚洲色图 男人天堂 中文字幕| 日本免费一区二区三区高清不卡| 免费在线观看成人毛片| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久成人av| 1024香蕉在线观看| 免费在线观看日本一区| 天天添夜夜摸| 免费在线观看视频国产中文字幕亚洲| 久久亚洲真实| 欧美日韩黄片免| 午夜激情福利司机影院| 亚洲欧美日韩无卡精品| 成人18禁在线播放| 亚洲av免费在线观看| 免费一级毛片在线播放高清视频| 亚洲性夜色夜夜综合| 国内精品一区二区在线观看| 久久精品91无色码中文字幕| 国产99白浆流出| 一进一出抽搐动态| 国产成人系列免费观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产欧洲综合997久久,| 国产亚洲精品综合一区在线观看| а√天堂www在线а√下载| 精品福利观看| 久久午夜亚洲精品久久| 又粗又爽又猛毛片免费看| 欧美中文日本在线观看视频| 久久精品国产亚洲av香蕉五月| 欧美大码av| 成年女人永久免费观看视频| 亚洲自偷自拍图片 自拍| 亚洲欧美日韩卡通动漫| 一级黄色大片毛片| 久久九九热精品免费| 国产爱豆传媒在线观看| 国产精品亚洲av一区麻豆| 丰满的人妻完整版| 免费观看精品视频网站| 久久中文字幕一级| 黄色日韩在线| 国产淫片久久久久久久久 | 欧美3d第一页| 精品国产美女av久久久久小说| 亚洲 欧美一区二区三区| 免费观看人在逋| 在线观看66精品国产| 操出白浆在线播放| 午夜影院日韩av| 国产精品 欧美亚洲| 波多野结衣高清作品| 香蕉久久夜色| 在线观看午夜福利视频| 亚洲专区中文字幕在线| 性色avwww在线观看| 国产一区二区三区视频了| aaaaa片日本免费| 99久久精品热视频| 日日夜夜操网爽| 狠狠狠狠99中文字幕| 一区二区三区激情视频| 亚洲人成伊人成综合网2020| av欧美777| 怎么达到女性高潮| 制服丝袜大香蕉在线| 在线观看免费视频日本深夜| 99久久国产精品久久久| 国产不卡一卡二| 中文字幕久久专区| 一个人免费在线观看电影 | 男人和女人高潮做爰伦理| 亚洲乱码一区二区免费版| svipshipincom国产片| 久久香蕉国产精品| 成人高潮视频无遮挡免费网站| 久久性视频一级片| 国产男靠女视频免费网站| 黄色丝袜av网址大全| 成人特级av手机在线观看| 婷婷亚洲欧美| 成人av在线播放网站| 丁香欧美五月| 国产av一区在线观看免费|