• <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)展
    山西大豆自然群體遺傳多樣性的研究
    精品人妻偷拍中文字幕| 国产一区有黄有色的免费视频| 又爽又黄无遮挡网站| 热99国产精品久久久久久7| 久久久久九九精品影院| 男人和女人高潮做爰伦理| 亚洲成人久久爱视频| 久久久久性生活片| 如何舔出高潮| 不卡视频在线观看欧美| 美女主播在线视频| 午夜亚洲福利在线播放| 如何舔出高潮| 国产精品一区二区在线观看99| 欧美亚洲 丝袜 人妻 在线| 国产男人的电影天堂91| 一本一本综合久久| 欧美成人午夜免费资源| 久久精品国产亚洲av天美| 精品国产三级普通话版| 国产伦理片在线播放av一区| 新久久久久国产一级毛片| 亚洲国产色片| 婷婷色综合www| 一区二区三区免费毛片| 六月丁香七月| 亚洲色图综合在线观看| 久久97久久精品| 国产精品人妻久久久影院| 啦啦啦在线观看免费高清www| 亚洲av福利一区| 久久人人爽人人爽人人片va| 精品午夜福利在线看| 最新中文字幕久久久久| 黄片wwwwww| freevideosex欧美| 日韩av不卡免费在线播放| 日韩av不卡免费在线播放| 国产有黄有色有爽视频| 日韩国内少妇激情av| 狂野欧美激情性bbbbbb| 亚洲国产色片| 97精品久久久久久久久久精品| 免费大片黄手机在线观看| 日韩一区二区视频免费看| 80岁老熟妇乱子伦牲交| 日韩在线高清观看一区二区三区| 国产v大片淫在线免费观看| 99热这里只有精品一区| 老师上课跳d突然被开到最大视频| 久久99热6这里只有精品| 欧美成人精品欧美一级黄| 女人被狂操c到高潮| 中文字幕制服av| 一区二区三区四区激情视频| 看非洲黑人一级黄片| 亚洲欧美成人综合另类久久久| 欧美日韩视频精品一区| 尾随美女入室| 成年女人在线观看亚洲视频 | 又爽又黄无遮挡网站| 亚洲精品456在线播放app| 免费看av在线观看网站| 日本一二三区视频观看| 男人和女人高潮做爰伦理| 亚洲欧美中文字幕日韩二区| 亚洲欧美中文字幕日韩二区| 男人爽女人下面视频在线观看| 久久久久精品性色| 在线 av 中文字幕| 国产美女午夜福利| 少妇熟女欧美另类| 蜜桃亚洲精品一区二区三区| 别揉我奶头 嗯啊视频| 国产成人精品一,二区| 人妻一区二区av| 日韩欧美 国产精品| 亚洲自偷自拍三级| 麻豆乱淫一区二区| 岛国毛片在线播放| 97超视频在线观看视频| 51国产日韩欧美| 赤兔流量卡办理| 亚洲一区二区三区欧美精品 | 精品少妇黑人巨大在线播放| 少妇人妻 视频| 亚洲婷婷狠狠爱综合网| 秋霞伦理黄片| 精品人妻视频免费看| 日本猛色少妇xxxxx猛交久久| av在线播放精品| 深爱激情五月婷婷| 免费人成在线观看视频色| 如何舔出高潮| 欧美3d第一页| 狂野欧美激情性bbbbbb| 国产一级毛片在线| 亚洲av不卡在线观看| 国产v大片淫在线免费观看| 久久精品熟女亚洲av麻豆精品| 97超碰精品成人国产| 三级男女做爰猛烈吃奶摸视频| 肉色欧美久久久久久久蜜桃 | 七月丁香在线播放| 又粗又硬又长又爽又黄的视频| 嫩草影院精品99| 毛片一级片免费看久久久久| 久久久久九九精品影院| 免费黄色在线免费观看| 女的被弄到高潮叫床怎么办| 国产 一区 欧美 日韩| 免费av观看视频| 波野结衣二区三区在线| 国产永久视频网站| 国产黄频视频在线观看| 国产午夜精品一二区理论片| 亚洲三级黄色毛片| 国产精品99久久久久久久久| 99久久九九国产精品国产免费| 欧美老熟妇乱子伦牲交| 联通29元200g的流量卡| 精品酒店卫生间| 国产精品一区二区性色av| 国产精品无大码| 色吧在线观看| 狂野欧美激情性xxxx在线观看| 免费观看无遮挡的男女| 中文乱码字字幕精品一区二区三区| 在线免费观看不下载黄p国产| 中文字幕制服av| 精品久久久噜噜| 一个人看视频在线观看www免费| 日日摸夜夜添夜夜添av毛片| 欧美激情在线99| 69人妻影院| 国产视频首页在线观看| 少妇人妻久久综合中文| 亚洲av成人精品一二三区| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av免费高清在线观看| 亚洲精品第二区| 国产视频内射| 又爽又黄a免费视频| 色5月婷婷丁香| 国产av码专区亚洲av| 日韩av在线免费看完整版不卡| 99久久精品热视频| 亚洲精品国产色婷婷电影| 性色avwww在线观看| 久久午夜福利片| 黄色视频在线播放观看不卡| 国产欧美日韩精品一区二区| 寂寞人妻少妇视频99o| 亚洲av.av天堂| 午夜福利视频1000在线观看| av专区在线播放| 日韩av不卡免费在线播放| 精品一区在线观看国产| 三级国产精品片| 男女那种视频在线观看| 国产高潮美女av| 麻豆国产97在线/欧美| 熟女人妻精品中文字幕| 尾随美女入室| 街头女战士在线观看网站| 欧美日韩视频高清一区二区三区二| 干丝袜人妻中文字幕| 99re6热这里在线精品视频| av国产免费在线观看| 国精品久久久久久国模美| 亚洲美女搞黄在线观看| 99视频精品全部免费 在线| 亚洲av男天堂| 大香蕉97超碰在线| 赤兔流量卡办理| av播播在线观看一区| 日本爱情动作片www.在线观看| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 亚洲精品乱久久久久久| 亚洲精品成人av观看孕妇| 亚洲无线观看免费| eeuss影院久久| 欧美成人a在线观看| av黄色大香蕉| 欧美亚洲 丝袜 人妻 在线| 亚洲成人一二三区av| 日本午夜av视频| 国产成人91sexporn| 黄色怎么调成土黄色| 麻豆乱淫一区二区| 看黄色毛片网站| 欧美精品人与动牲交sv欧美| 久久久国产一区二区| 99久国产av精品国产电影| 老女人水多毛片| 久久女婷五月综合色啪小说 | 亚洲av福利一区| 国产老妇伦熟女老妇高清| 真实男女啪啪啪动态图| 高清午夜精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 女人十人毛片免费观看3o分钟| 欧美一级a爱片免费观看看| 亚洲欧美精品专区久久| 在线 av 中文字幕| 免费大片黄手机在线观看| 国产亚洲91精品色在线| 成年av动漫网址| 欧美精品人与动牲交sv欧美| 插阴视频在线观看视频| 欧美老熟妇乱子伦牲交| 国产黄a三级三级三级人| 女的被弄到高潮叫床怎么办| 国产亚洲5aaaaa淫片| 久久久久九九精品影院| 成年版毛片免费区| av在线天堂中文字幕| 国产亚洲av片在线观看秒播厂| 日本-黄色视频高清免费观看| 亚洲av免费在线观看| 好男人视频免费观看在线| 人妻 亚洲 视频| 午夜福利高清视频| 丰满人妻一区二区三区视频av| 能在线免费看毛片的网站| videos熟女内射| 欧美三级亚洲精品| 日韩大片免费观看网站| 一区二区三区乱码不卡18| 伊人久久国产一区二区| 波野结衣二区三区在线| av免费在线看不卡| 春色校园在线视频观看| 又黄又爽又刺激的免费视频.| 欧美xxⅹ黑人| 大片电影免费在线观看免费| 国产精品一区www在线观看| 国产精品人妻久久久久久| 亚洲av.av天堂| 欧美丝袜亚洲另类| 我的老师免费观看完整版| 日日摸夜夜添夜夜爱| 成人午夜精彩视频在线观看| 亚洲精品乱久久久久久| 亚洲精品乱码久久久v下载方式| 国产乱人视频| 一本一本综合久久| 精品国产一区二区三区久久久樱花 | 亚洲精品久久久久久婷婷小说| 能在线免费看毛片的网站| 人妻一区二区av| 国产亚洲91精品色在线| 国产综合懂色| 久久久成人免费电影| 人妻系列 视频| 九九在线视频观看精品| 亚洲国产精品成人综合色| 午夜视频国产福利| 日韩av免费高清视频| 亚洲av在线观看美女高潮| 中文欧美无线码| av女优亚洲男人天堂| 成年版毛片免费区| 久久精品久久精品一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 久久久亚洲精品成人影院| 欧美成人一区二区免费高清观看| 成年版毛片免费区| 亚洲,一卡二卡三卡| 人妻 亚洲 视频| 一本色道久久久久久精品综合| 欧美 日韩 精品 国产| 男人舔奶头视频| 亚洲四区av| 一个人看的www免费观看视频| 国产精品女同一区二区软件| 国产午夜精品一二区理论片| 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| 少妇人妻久久综合中文| 高清毛片免费看| 只有这里有精品99| 国产淫片久久久久久久久| 夫妻性生交免费视频一级片| 亚洲色图综合在线观看| 免费看光身美女| 噜噜噜噜噜久久久久久91| 又黄又爽又刺激的免费视频.| 日韩伦理黄色片| 久久久久性生活片| 神马国产精品三级电影在线观看| 午夜爱爱视频在线播放| 欧美xxⅹ黑人| 国产亚洲精品久久久com| 九草在线视频观看| av国产精品久久久久影院| 97超碰精品成人国产| 一边亲一边摸免费视频| 精品久久久噜噜| 日韩伦理黄色片| 热re99久久精品国产66热6| 纵有疾风起免费观看全集完整版| 国产亚洲精品久久久com| 人人妻人人澡人人爽人人夜夜| 黄色怎么调成土黄色| 亚洲精华国产精华液的使用体验| 日韩强制内射视频| 亚洲精华国产精华液的使用体验| 91久久精品国产一区二区三区| 国产色婷婷99| 国产成年人精品一区二区| 亚洲久久久久久中文字幕| 国产午夜精品久久久久久一区二区三区| 国产精品嫩草影院av在线观看| 网址你懂的国产日韩在线| 嫩草影院新地址| 欧美潮喷喷水| av免费观看日本| 男女边摸边吃奶| 在线观看一区二区三区| 日韩欧美精品免费久久| 99久久人妻综合| 麻豆久久精品国产亚洲av| 丰满乱子伦码专区| 在线免费观看不下载黄p国产| 一级毛片我不卡| 精品99又大又爽又粗少妇毛片| 99re6热这里在线精品视频| av国产免费在线观看| 老司机影院成人| 国产一区二区三区综合在线观看 | 可以在线观看毛片的网站| 男的添女的下面高潮视频| 久久精品综合一区二区三区| 午夜精品国产一区二区电影 | 麻豆国产97在线/欧美| 如何舔出高潮| 免费少妇av软件| 亚洲av二区三区四区| 国产男女内射视频| 男人狂女人下面高潮的视频| 高清午夜精品一区二区三区| 中国国产av一级| 免费黄频网站在线观看国产| 亚洲精品视频女| av国产精品久久久久影院| 夜夜爽夜夜爽视频| 国产精品国产三级专区第一集| 天天躁日日操中文字幕| 在线亚洲精品国产二区图片欧美 | 最近2019中文字幕mv第一页| 日韩在线高清观看一区二区三区| 日本午夜av视频| 国产黄片视频在线免费观看| 少妇 在线观看| 最近手机中文字幕大全| 国产一区有黄有色的免费视频| eeuss影院久久| 热99国产精品久久久久久7| 大又大粗又爽又黄少妇毛片口| 国产精品国产三级专区第一集| 国产一区有黄有色的免费视频| 成人午夜精彩视频在线观看| 一区二区三区精品91| 校园人妻丝袜中文字幕| 国产日韩欧美亚洲二区| 亚洲精品视频女| 久久久精品94久久精品| 久久99蜜桃精品久久| 亚洲av一区综合| 熟女人妻精品中文字幕| 欧美日本视频| 精品一区二区免费观看| 国产淫片久久久久久久久| 精品视频人人做人人爽| 少妇被粗大猛烈的视频| 欧美激情在线99| 国产午夜精品一二区理论片| 亚洲综合色惰| 欧美丝袜亚洲另类| 你懂的网址亚洲精品在线观看| 青青草视频在线视频观看| 国产精品麻豆人妻色哟哟久久| 亚洲成人中文字幕在线播放| 老司机影院成人| 三级经典国产精品| 久久久久久国产a免费观看| 熟女人妻精品中文字幕| 国产爱豆传媒在线观看| 一级黄片播放器| 国产爱豆传媒在线观看| 好男人视频免费观看在线| 国产成人一区二区在线| 九九久久精品国产亚洲av麻豆| 性色avwww在线观看| 男女无遮挡免费网站观看| 22中文网久久字幕| 嫩草影院新地址| 男女无遮挡免费网站观看| 欧美日韩视频精品一区| 国产精品久久久久久精品电影小说 | 在线 av 中文字幕| 亚洲经典国产精华液单| a级毛片免费高清观看在线播放| av在线app专区| 亚洲熟女精品中文字幕| 国产 精品1| 国产精品久久久久久久电影| 日韩成人av中文字幕在线观看| 午夜精品一区二区三区免费看| 精品久久久久久久久亚洲| 久久久久国产精品人妻一区二区| 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 久久热精品热| 亚洲成人一二三区av| 少妇被粗大猛烈的视频| 国产综合懂色| 亚洲国产欧美人成| 在线精品无人区一区二区三 | 丝瓜视频免费看黄片| 高清日韩中文字幕在线| 大又大粗又爽又黄少妇毛片口| 国产男女内射视频| 国产毛片在线视频| 久久97久久精品| 久久女婷五月综合色啪小说 | 国产黄片美女视频| 2021少妇久久久久久久久久久| 欧美精品一区二区大全| 日韩一区二区三区影片| 99热国产这里只有精品6| 丝瓜视频免费看黄片| 日韩大片免费观看网站| 欧美变态另类bdsm刘玥| 午夜亚洲福利在线播放| 欧美最新免费一区二区三区| 久久久久网色| 亚洲精品久久久久久婷婷小说| 美女视频免费永久观看网站| 99热网站在线观看| 美女被艹到高潮喷水动态| 综合色av麻豆| 看免费成人av毛片| 婷婷色综合www| 久久久久久九九精品二区国产| 亚洲精品第二区| 97人妻精品一区二区三区麻豆| 国产久久久一区二区三区| 亚州av有码| 赤兔流量卡办理| 久久ye,这里只有精品| 伊人久久精品亚洲午夜| 大片免费播放器 马上看| 亚洲欧美精品自产自拍| 青春草国产在线视频| 深夜a级毛片| 亚洲久久久久久中文字幕| 亚洲精品乱久久久久久| 好男人在线观看高清免费视频| 综合色丁香网| 一级毛片aaaaaa免费看小| 18禁裸乳无遮挡动漫免费视频 | 成人毛片a级毛片在线播放| 一本久久精品| 久久女婷五月综合色啪小说 | 国产日韩欧美亚洲二区| 日韩精品有码人妻一区| 特大巨黑吊av在线直播| 国产精品福利在线免费观看| 日韩一区二区三区影片| av网站免费在线观看视频| 日韩免费高清中文字幕av| 国产老妇女一区| 国产黄片视频在线免费观看| 91狼人影院| 国产亚洲91精品色在线| 校园人妻丝袜中文字幕| 色视频www国产| 在线天堂最新版资源| 国产精品爽爽va在线观看网站| a级一级毛片免费在线观看| 一级a做视频免费观看| 真实男女啪啪啪动态图| 精品久久久久久电影网| 极品教师在线视频| 亚洲熟女精品中文字幕| 亚洲av不卡在线观看| 亚洲精品成人久久久久久| 最近的中文字幕免费完整| 精品熟女少妇av免费看| 国产精品久久久久久av不卡| 国产成人精品福利久久| 亚洲,一卡二卡三卡| 欧美人与善性xxx| 国产精品女同一区二区软件| 欧美日本视频| 亚洲最大成人av| 日韩成人伦理影院| 成人漫画全彩无遮挡| 国产又色又爽无遮挡免| 久久久国产一区二区| 国产一区有黄有色的免费视频| 丰满人妻一区二区三区视频av| 免费黄网站久久成人精品| 亚洲精品日韩av片在线观看| 国产视频内射| 自拍偷自拍亚洲精品老妇| 韩国高清视频一区二区三区| 日韩av在线免费看完整版不卡| 最近2019中文字幕mv第一页| 日韩av在线免费看完整版不卡| 亚洲成人久久爱视频| 国产一区二区在线观看日韩| 午夜精品一区二区三区免费看| 久久综合国产亚洲精品| 欧美日韩一区二区视频在线观看视频在线 | 国产 一区 欧美 日韩| 久久ye,这里只有精品| 欧美激情久久久久久爽电影| av女优亚洲男人天堂| av国产免费在线观看| 欧美一级a爱片免费观看看| 综合色丁香网| 夫妻性生交免费视频一级片| 国产中年淑女户外野战色| 九九久久精品国产亚洲av麻豆| 80岁老熟妇乱子伦牲交| 一二三四中文在线观看免费高清| 国产精品一区二区性色av| 香蕉精品网在线| 黑人高潮一二区| 中国美白少妇内射xxxbb| 在线观看免费高清a一片| 国产永久视频网站| 久久久久久久国产电影| 国产一区二区亚洲精品在线观看| 亚洲经典国产精华液单| 秋霞在线观看毛片| 成人国产av品久久久| 国国产精品蜜臀av免费| 五月开心婷婷网| 国产成人aa在线观看| 国产伦精品一区二区三区视频9| 小蜜桃在线观看免费完整版高清| 国模一区二区三区四区视频| 亚洲精品,欧美精品| 国产色爽女视频免费观看| 欧美性猛交╳xxx乱大交人| 国产精品久久久久久久电影| 又爽又黄无遮挡网站| 国产亚洲91精品色在线| 精品酒店卫生间| 亚洲人与动物交配视频| 久久99热这里只有精品18| 欧美变态另类bdsm刘玥| 能在线免费看毛片的网站| 国产精品成人在线| 国产色爽女视频免费观看| 日韩成人av中文字幕在线观看| 国产精品一区二区性色av| 欧美极品一区二区三区四区| 一级毛片电影观看| 精品国产一区二区三区久久久樱花 | 午夜福利视频1000在线观看| 国产 一区精品| 成年女人在线观看亚洲视频 | 两个人的视频大全免费| 99热这里只有精品一区| 噜噜噜噜噜久久久久久91| 国产精品久久久久久久电影| 91精品国产九色| 男女边吃奶边做爰视频| 亚洲av二区三区四区| 久久99热6这里只有精品| xxx大片免费视频| 少妇高潮的动态图| 99热6这里只有精品| 亚洲精品乱久久久久久| 九草在线视频观看| 少妇人妻精品综合一区二区| 日韩大片免费观看网站| 久久久久国产网址| 在线免费观看不下载黄p国产| 成人国产av品久久久| 国产精品久久久久久av不卡| 2022亚洲国产成人精品| 国产成人精品一,二区| 国产成人免费无遮挡视频| 久久久午夜欧美精品| 国产成人a区在线观看| 91久久精品国产一区二区成人| 午夜福利在线在线| 欧美一级a爱片免费观看看| 国产成人a∨麻豆精品| 一区二区三区精品91| 男人添女人高潮全过程视频| 91久久精品国产一区二区成人| 2022亚洲国产成人精品| 夜夜爽夜夜爽视频| 大片免费播放器 马上看| 国产午夜精品久久久久久一区二区三区| 男人添女人高潮全过程视频| 精品久久久久久久末码| 国产精品福利在线免费观看| 久久精品熟女亚洲av麻豆精品| 亚洲国产最新在线播放| 亚洲成人精品中文字幕电影| 日韩一区二区视频免费看| 国产在视频线精品| 少妇丰满av| 91狼人影院| 熟女人妻精品中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 |