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

    凡納濱對(duì)蝦進(jìn)口親蝦子一代群體的遺傳變異分析

    2021-03-05 20:47:41徐煜黃小帥胡曉娟徐武杰蘇浩昌文國樑曹煜成楊鏗
    關(guān)鍵詞:遺傳變異

    徐煜 黃小帥 胡曉娟 徐武杰 蘇浩昌 文國樑 曹煜成 楊鏗

    摘要:【目的】明確凡納濱對(duì)蝦(Litopenaeus vannamei)進(jìn)口親蝦直接繁殖子一代種苗群體的遺傳變異情況,旨在評(píng)估遺傳特征與養(yǎng)殖生產(chǎn)性能的關(guān)聯(lián)性,為我國養(yǎng)殖對(duì)蝦優(yōu)質(zhì)種苗大數(shù)據(jù)建設(shè)及質(zhì)量評(píng)估體系的完善提供基礎(chǔ)數(shù)據(jù)支持?!痉椒ā恳詠碓从谡笊駷澄r苗場(chǎng)的Stock 1~Stock 5群體、正大漳浦蝦苗場(chǎng)的Stock 6群體、海尚種苗培育基地的Stock 7和Stock 8群體及中海水產(chǎn)種苗科技有限公司的Stock 9群體等9個(gè)不同批次凡納濱對(duì)蝦進(jìn)口親蝦子一代種苗為研究對(duì)象,利用11個(gè)微衛(wèi)星分子標(biāo)記(M1、Pvan1758、Pvan1815、HLJN-008、HLJN-023、TUMXLv7.121、TUMXLv8.256、TUMXLv9.43、TUMXLv10.312、TUYFLvL16.1a和TUDGLv1-3.224)對(duì)不同凡納濱對(duì)蝦種苗群體進(jìn)行遺傳變異分析?!窘Y(jié)果】9個(gè)凡納濱對(duì)蝦種苗群體的等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、觀測(cè)雜合度(Ho)、期望雜合度(He)、Shannon’s信息指數(shù)(I)和多態(tài)信息含量(PIC)的平均值分別在5.5455~8.8182、3.2795~5.4735、0.3769~0.5182、0.6189~0.6948、1.2823~1.6092和0.5863~0.6646,即各凡納濱對(duì)蝦種苗群體均具有較高的遺傳多樣性,具體排序?yàn)镾tock 3>Stock 2>Stock 4>Stock 5>Stock 7>Stock 6>Stock 9>Stock 8>Stock 1。在9個(gè)凡納濱對(duì)蝦種苗群體中有91.82%的遺傳變異源自群體內(nèi),僅有8.18%的變異源自群體間。Stock 2、Stock 3、Stock 4和Stock 5等4個(gè)群體間及Stock 7群體與Stock 8群體間的遺傳分化系數(shù)(Fst)均小于0.05(0.0079~0.0371),即群體間無遺傳分化,而其他群體間表現(xiàn)為輕度遺傳分化(0.05<Fst <0.15)?;贜ei’s遺傳距離,9個(gè)凡納濱對(duì)蝦種苗群體大致可劃分為三大分支,第一分支包括Stock 1~Stock 6群體,第二分支由Stock 7群體和Stock 8群體組成,而Stock 9群體獨(dú)立為第三分支?!窘Y(jié)論】不同批次凡納濱對(duì)蝦進(jìn)口親蝦子一代種苗群體的遺傳特征存在明顯差異,各群體間的遺傳距離隨著采樣時(shí)間和地點(diǎn)的改變而發(fā)生明顯變化,且這些遺傳背景差異可能會(huì)造成不同群體在養(yǎng)殖生產(chǎn)性能方面產(chǎn)生差異。

    關(guān)鍵詞: 凡納濱對(duì)蝦;進(jìn)口親蝦;子一代群體;遺傳變異;微衛(wèi)星分子標(biāo)記

    中圖分類號(hào): S945.49? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼: A 文章編號(hào):2095-1191(2021)11-3139-08

    Genetic variation analysis among first filial generation of introduced Litopenaeus vannamei stocks

    XU Yu1,2, HUANG Xiao-shuai1, HU Xiao-juan1,2, XU Wu-jie1,2, SU Hao-chang1,2,

    WEN Guo-liang1, CAO Yu-cheng1,2*, YANG Keng1*

    (1South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences/Key Laboratory of South China Sea Fishery Resources Exploitation & Utilization, Ministry of Agriculture and Rural Affairs/Guangdong Provincial Key Laboratory of Fishery Ecology and Environment, Guangzhou? 510300, China; 2 Shenzhen Base of South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Shenzhen, Guangdong? 518121, China)

    Abstract:【Objective】The main purpose of this study was to analyze and assess the genetic diversity of first filial ge-neration of Litopenaeus vannamei, in order to evaluate the correlation between its genetic characteristics and aquaculture production performance, and provide basic data support for the big data construction and establishment of quality evaluation system of high-quality post-larvae in China. 【Method】The objects were L. vannamei Stock 1-Stock 5 from Zhengda Shenwan shrimp larvae farm, Stock 6 from Zhengda Zhangpu shrimp larvae farm, Stock 7 and Stock 8 from Haishang shrimp larvae farm, Stock 9 from Zhonghai aquatic larvae culture company. The genetic variation of nine first filial gene-ration stocks were investigated using eleven microsatellite DNA markers(M1, Pvan1758, Pvan1815, HLJN-008, HLJN-023, TUMXLv7.121, TUMXLv8.256, TUMXLv9.43, TUMXLv10.312, TUYFLvL16.1a and TUDGLv1-3.224). 【Result】This paper showed that the average values of the number of alleles (Na), number of effective alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), Shannon’s diversity index (I), and polymorphism information content (PIC) were among 5.5455-8.8182, 3.2795-5.4735, 0.3769-0.5182, 0.6189-0.6948, 1.2823-1.6092 and 0.5863-0.6646, respectively. That presented a high genetic diversity of these filial generation of introduced L. vannamei stocks, and according to the richness, they ranked as follows: Stock 3>Stock 2>Stock 4>Stock 5>Stock 7>Stock 6>Stock 9>Stock 8>Stock 1. AMOVA analysis indicated that the genetic variation within stock was 91.82% and the remaining 8.18% of the variation occurred among stocks. The genetic differentiation among Stocks 2-5 and among Stocks 7-8 showed a low level based on their Fst values (0.0079-0.0371), which were lower than 0.05, whereas the other stocks had no genetic differentiation (0.05<Fst <0.15). Based on Nei’s genetic distance, the nine stocks could be roughly divided into three branches: the first branch including Stock 1 - Stock 6 stocks, the second branch being composed of Stock 7 and Stock 8, and the third branch only included Stock 9. 【Conclusion】These findings indicated that different batches of first filial generations have distinct genetic characteristics. Dramatically, genetic distance of first filial generations are clearly fluctuated with the breeding farm and breeding time. Changes of genetic background might directly affect the shrimp traits in subsequent farming.

    Key words: Litopenaeus vannamei; exotic parents; first filial generation; genetic variation; microsatellite molecular marker

    Foundation item: Central Public-interest Scientific Institution Basal Research Fund of Chinese Academy of Fishery Sciences(2020TD54); China Agriculture Research System of MOF and MARA(CARS-48); Guangdong Provincial Special Fund for Modern Agriculture Industry Technology Innovation Teams (2019KJ149)

    0 引言

    【研究意義】凡納濱對(duì)蝦(Litopenaeus vannamei)又稱南美白對(duì)蝦,隸屬于十足目(Decapoda)游泳亞目(Natantia)對(duì)蝦科(Penaeidae),是當(dāng)今最主要的對(duì)蝦養(yǎng)殖品種。凡納濱對(duì)蝦自1988年引入后逐步發(fā)展成為我國對(duì)蝦養(yǎng)殖產(chǎn)業(yè)的主導(dǎo)品種,近年來我國的凡納濱對(duì)蝦養(yǎng)殖產(chǎn)量穩(wěn)居世界首位(代平等,2018)。凡納濱對(duì)蝦為熱帶型種類,原產(chǎn)自中南美洲太平洋沿岸水域,其優(yōu)良種質(zhì)資源長期受到國外育種機(jī)構(gòu)的壟斷(王興強(qiáng)等,2004;代平等,2018;彭敏等,2020)。當(dāng)前全球凡納濱對(duì)蝦種蝦資源主要分布在美國的夏威夷、邁阿密和關(guān)島,以及泰國、新加坡和印尼,我國每年需從境外進(jìn)口大量親蝦,由進(jìn)口親蝦直接繁育子一代種苗具有生長速度快、養(yǎng)成群體個(gè)頭均勻等優(yōu)點(diǎn),在我國對(duì)蝦養(yǎng)殖業(yè)中占有較高的市場(chǎng)份額。但受種質(zhì)資源退化及養(yǎng)殖環(huán)境惡化等因素的影響,對(duì)蝦養(yǎng)殖業(yè)面臨病害頻發(fā)、養(yǎng)殖風(fēng)險(xiǎn)劇增等突出問題(王麗花等,2012;文國樑等,2015;Qiu et al.,2017)。從不同批次進(jìn)口親蝦子一代種苗的養(yǎng)殖實(shí)踐來看,蝦苗的成活率、抗逆性及生長周期等表現(xiàn)良莠不齊,嚴(yán)重影響了養(yǎng)殖經(jīng)濟(jì)效益。因此,針對(duì)凡納濱對(duì)蝦進(jìn)口親蝦直接繁殖子一代種苗群體開展遺傳變異分析,對(duì)評(píng)估其遺傳特征與養(yǎng)殖生產(chǎn)性能的關(guān)聯(lián)性具有重要意義?!厩叭搜芯窟M(jìn)展】微衛(wèi)星分子標(biāo)記分析具有簡便、快捷及穩(wěn)定性好等優(yōu)點(diǎn),已廣泛應(yīng)用于種質(zhì)資源評(píng)價(jià)和遺傳育種等研究領(lǐng)域(孫效文等,2008;趙志英等,2018;周大顏等,2019)。近年來,微衛(wèi)星分子標(biāo)記也常應(yīng)用于水生生物遺傳多樣性分析。劉連為等(2014)利用8個(gè)微衛(wèi)星位點(diǎn)對(duì)北太平洋柔魚(Ommastrephes bartramii)的6個(gè)群體進(jìn)行分析,結(jié)果發(fā)現(xiàn)群體間具有較高的遺傳多樣性。毛守康等(2016)利用14個(gè)微衛(wèi)星位點(diǎn)對(duì)4個(gè)野生梭魚(Liza haematocheila)地理群體進(jìn)行遺傳多樣性分析,結(jié)果表明,4個(gè)群體遺傳多樣性處于中等多態(tài)水平,且按其地理位置由北向南依次聚為一類。Wang等(2016)利用8個(gè)微衛(wèi)星位點(diǎn)分析中國對(duì)蝦親本和回捕群體的遺傳多樣性,發(fā)現(xiàn)這2個(gè)種群擁有豐富遺傳多樣性,但存在較高程度的近親繁殖,且在中國對(duì)蝦種群中呈現(xiàn)一定的近交衰退。黃小芳等(2020)利用8對(duì)微衛(wèi)星引物對(duì)廣西不同地區(qū)克氏原螯蝦群體的遺傳多樣性進(jìn)行分析,結(jié)果發(fā)現(xiàn)各克氏原螯蝦群體具有一定的遺傳多樣性,但其雜合子缺失現(xiàn)象普遍存在,應(yīng)通過適當(dāng)引種及加強(qiáng)不同地區(qū)群體間的基因交流,保護(hù)好優(yōu)良種質(zhì)資源。在凡納濱對(duì)蝦方面,童馨等(2007,2009)在研究不同世代凡納濱對(duì)蝦養(yǎng)殖群體生長性狀變異的基礎(chǔ)上,利用微衛(wèi)星分子標(biāo)記分析種群遺傳多樣性的變化,發(fā)現(xiàn)隨著繁育世代增加,凡納濱對(duì)蝦的生長性狀逐代分化,生長速度變緩,群體遺傳多樣性降低,種質(zhì)發(fā)生退化;吳怡迪等(2016)基于微衛(wèi)星分子標(biāo)記對(duì)凡納濱對(duì)蝦選育系間雜交的生長性狀及遺傳多樣性進(jìn)行分析,結(jié)果表明,凡納濱對(duì)蝦的生長性能與群體遺傳多樣性水平呈顯著正相關(guān),選育系間雜交可使雜交后代的生長性狀和遺傳多樣性水平得到改善;李強(qiáng)勇等(2020)基于第三代高通量測(cè)序技術(shù)——單分子實(shí)時(shí)(SMRT)測(cè)序開發(fā)出凡納濱對(duì)蝦微衛(wèi)星分子標(biāo)記,有助于開展凡納濱對(duì)蝦遺傳育種研究工作;彭敏等(2020)從GenBank中篩選11對(duì)能有效擴(kuò)增的微衛(wèi)星引物對(duì)2016—2018年桂海1號(hào)凡納濱對(duì)蝦連續(xù)3個(gè)世代選育群體進(jìn)行遺傳多樣性分析,監(jiān)測(cè)相鄰2個(gè)世代選育群體間的遺傳變異情況,結(jié)果顯示桂海1號(hào)凡納濱對(duì)蝦各世代選育群體尚具有較豐富的遺傳變異,但群體遺傳結(jié)構(gòu)處于不穩(wěn)定狀態(tài)。【本研究切入點(diǎn)】遺傳變異是影響對(duì)蝦養(yǎng)殖性能表現(xiàn)的重要因素,但近年來鮮見從遺傳變異水平研究進(jìn)口親蝦子一代種苗養(yǎng)殖性能差異的問題?!緮M解決的關(guān)鍵問題】選取不同批次凡納濱對(duì)蝦進(jìn)口親蝦子一代種苗群體,利用微衛(wèi)星分子標(biāo)記分析不同群體的遺傳變異及其多樣性水平,旨在評(píng)估遺傳特征與養(yǎng)殖生產(chǎn)性能的關(guān)聯(lián)性,為我國養(yǎng)殖對(duì)蝦優(yōu)質(zhì)種苗大數(shù)據(jù)建設(shè)及質(zhì)量評(píng)估體系的完善提供基礎(chǔ)數(shù)據(jù)支持。

    1 材料與方法

    1. 1 試驗(yàn)材料

    供試凡納濱對(duì)蝦種苗群體為2017—2018年收集于廣東中山、茂名、汕尾及福建漳州等地的蝦苗場(chǎng)或養(yǎng)殖場(chǎng),共9個(gè)群體,分別命名為Stock 1~Stock 9。各凡納濱對(duì)蝦種苗群體樣品具體信息如表1所示。

    1. 2 DNA提取

    取凡納濱對(duì)蝦種苗腹部組織10~30 mg,參照E.Z.N.A. Tissue DNA Kit (OMEGA Bio-tek)試劑盒說明提取各樣品基因組DNA。采用1.0%瓊脂糖凝膠電泳檢測(cè)DNA完整性,以NanoDrop 2000紫外分光光度計(jì)測(cè)定其濃度,-80 ℃保存?zhèn)溆谩?/p>

    1. 3 微衛(wèi)星位點(diǎn)

    根據(jù)相關(guān)文獻(xiàn)報(bào)道,選取凡納濱對(duì)蝦中多態(tài)性高且擴(kuò)增效果好的11個(gè)微衛(wèi)星位點(diǎn):M1(Garcia et al.,1996;Alcivar-Warren et al.,2007a)、Pvan1758(Cruz et al.,2002)、Pvan1815(Cruz et al.,2002)、HLJN-008(Jia et al.,2006)、HLJN-023(Jia et al.,2006)、TUMXLv7.121(Meehan et al.,2003)、TUMX Lv8.256(Meehan et al.,2003;Alcivar-Warren et al.,2007b)、TUMXLv9.43(Meehan et al.,2003;Alcivar-Warren et al.,2007b)、TUMXLv10.312(Meehan et al.,2003;Alcivar-Warren et al.,2007b)、TUYFLvL16.1a(Alcivar-Warren et al.,2007b)和TUDGLv1-3.224(Alcivar-Warren et al.,2007b;Garcia et al.,2007),在各位點(diǎn)上游引物5'端分別進(jìn)行6-FAM、ROX、HEX或TAMRA熒光標(biāo)記,引物序列及熒光標(biāo)記類型詳見表2,并委托北京睿博興科生物技術(shù)有限公司合成。

    1. 4 PCR擴(kuò)增及產(chǎn)物檢測(cè)

    PCR反應(yīng)體系15.0 μL:DNA模板(50 ng/μL)1.0 μL,Premix TaqTM Mix(TaKaRa)7.5 μL,上、下游引物(0.5 μmol/L)各0.5 μL,雙蒸水5.5 μL。擴(kuò)增程序:95 ℃預(yù)變性3 min;95 ℃ 30 s,50 ℃ 30 s,72 ℃ 20 s,進(jìn)行32個(gè)循環(huán);72 ℃延伸6 min。以1.0%瓊脂糖凝膠電泳檢測(cè)PCR擴(kuò)增產(chǎn)物,并委托北京睿博興科生物技術(shù)有限公司,根據(jù)產(chǎn)物熒光特性進(jìn)行毛細(xì)管電泳混檢,測(cè)序結(jié)果使用Gene Marker 2.2.0讀取。

    1. 5 數(shù)據(jù)分析

    根據(jù)各凡納濱對(duì)蝦種苗群體在各微衛(wèi)星位點(diǎn)的基因型,利用GENE Stock 1.32進(jìn)行統(tǒng)計(jì)分析(Yeh and Boyle,1997),計(jì)算各微衛(wèi)星位點(diǎn)的等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、觀測(cè)雜合度(Ho)、期望雜合度(He)、Shannon’s信息指數(shù)(I)及Nei’s遺傳距離;采用PIC_CALC 0.6(https://github.com/luansheng/PIC_CALC)計(jì)算多態(tài)信息含量(PIC);運(yùn)用Arlequin 3.1計(jì)算群體間遺傳分化系數(shù)(Fst)及進(jìn)行分子變異分析(AMOVA) (Excoffier et al.,2005);同時(shí)基于Nei’s遺傳距離,通過MEGA 5.0以UPGMA法進(jìn)行聚類分析(Tamura et al.,2011)。

    2 結(jié)果與分析

    2. 1 微衛(wèi)星位點(diǎn)多態(tài)性分析結(jié)果

    利用11對(duì)微衛(wèi)星引物對(duì)9個(gè)凡納濱對(duì)蝦種苗群體共計(jì)258個(gè)個(gè)體的基因組進(jìn)行擴(kuò)增,結(jié)果(圖1和圖2)顯示,11個(gè)微衛(wèi)星位點(diǎn)均能進(jìn)行有效擴(kuò)增,其Na、Ne、Ho、He、I及PIC等信息見表3。各微衛(wèi)星位點(diǎn)的Na為5~39,Ne為1.6582~14.5320,Ho為0.2023~0.6357,He為0.3969~0.9312,I為0.8357~3.0296,PIC為0.3789~0.9274。其中,TUMXLv8.256、TUMXLv10.312和TUYFLvL16.1a位點(diǎn)的PIC分別為0.3821、0.4605和0.3789,而其余8個(gè)微衛(wèi)星位點(diǎn)的PIC均大于0.6000。

    2. 2 凡納濱對(duì)蝦種苗群體遺傳多樣性分析結(jié)果

    9個(gè)凡納濱對(duì)蝦種苗群體的Na、Ne、Ho、He、I和PIC平均值分別在5.5455~8.8182、3.2795~5.4735、0.3769~0.5182、0.6189~0.6948、1.2823~1.6092和0.5863~ 0.6646(表4)。其中,Stock 3群體的Na、Ne和I平均值最高,Stock 2群體的He和PIC平均值最高,而Stock 1群體的Na、Ne、Ho、PIC和I平均值均最低。在9個(gè)凡納濱對(duì)蝦種苗群體中,He平均值以Stock 6群體最低,而Ho平均值以Stock 7群體最高??梢姡鞣布{濱對(duì)蝦種苗群體均具有較高的遺傳多樣性,遺傳多樣性水平排序?yàn)镾tock 3>Stock 2>Stock 4>Stock 5>Stock 7>Stock 6>Stock 9>Stock 8>Stock 1。

    2. 3 凡納濱對(duì)蝦種苗群體遺傳分化及遺傳距離分析結(jié)果

    9個(gè)凡納濱對(duì)蝦種苗群體的AMOVA分析結(jié)果如表5所示,有91.82%的變異來自于群體內(nèi),僅有8.18%的變異來自于群體間。兩兩群體間的Fst范圍在0.0079~0.1494(表6),即整體遺傳分化表現(xiàn)為中低水平。其中,Stock 6群體與Stock 9群體間的Fst最高(0.1494),而Stock 2群體與Stock 4群體間的Fst最低(0.0079)。Stock 2、Stock 3、Stock 4和Stock 5等4個(gè)群體間及Stock 7群體與Stock 8群體間的Fst均小于0.05(0.0079~0.0371),即群體間無遺傳分化,而其他群體間表現(xiàn)為輕度遺傳分化(0.05<Fst <0.15)。

    9個(gè)凡納濱對(duì)蝦種苗群體間的Nei’s遺傳距離為0.0568~0.5084(表6),以Stock 2群體與Stock 4群體的Nei’s遺傳距離最近,而Stock 8群體與Stock 9群體的Nei’s遺傳距離最遠(yuǎn)。根據(jù)Nei’s遺傳距離構(gòu)建的UPGMA聚類圖(圖3)顯示,9個(gè)凡納濱對(duì)蝦種苗群體大致可劃分為三大分支,第一分支包括Stock 1~Stock 6群體,其中Stock 2群體和Stock 4群體先聚為一簇,然后依次同Stock 3、Stock 5、Stock 1和Stock 6等4個(gè)群體聚為一支;第二分支由Stock 7群體和Stock 8群體組成,而Stock 9群體獨(dú)立為第三分支。

    3 討論

    微衛(wèi)星分子標(biāo)記是物種群體遺傳多樣性分析的重要手段,而基因型分析是微衛(wèi)星分子標(biāo)記研究的關(guān)鍵環(huán)節(jié)。本研究采用基于毛細(xì)管電泳的熒光標(biāo)記引物進(jìn)行基因分型,選取11個(gè)具有多態(tài)性的微衛(wèi)星位點(diǎn),針對(duì)不同批次的凡納濱對(duì)蝦進(jìn)口親蝦一代苗群體開展遺傳多樣性分析,結(jié)果顯示各微衛(wèi)星位點(diǎn)的PIC為0.3789~0.9274。PIC能反映等位基因的多態(tài)性水平,PIC在0.25~0.50為中度多態(tài),在0.50之上則為高度多態(tài)(Shete et al.,2000;Nagy et al.,2012)。說明本研究選取的11個(gè)微衛(wèi)星分子標(biāo)記均呈中高度多態(tài),是評(píng)估群體遺傳多樣性的有效分子標(biāo)記。

    遺傳多樣性是種質(zhì)資源評(píng)價(jià)的重要依據(jù),Na、Ho、He和PIC等均是衡量群體遺傳多樣性的重要指標(biāo)。Zhang等(2014)選用7個(gè)微衛(wèi)星分子標(biāo)記對(duì)由新加坡和美國引進(jìn)的7個(gè)親蝦群體遺傳多樣性進(jìn)行分析,結(jié)果顯示,各群體的Na平均值為4.286~11.143,Ne平均值為3.090~6.407,Ho平均值為0.526~0.754,He平均值為0.600~0.802。本研究中,由引進(jìn)親蝦直接繁育的9個(gè)子一代苗群體的Na、Ne和He平均值均在以上相應(yīng)范圍內(nèi),僅Ho平均值小于上述指標(biāo),推測(cè)是由于各凡納濱對(duì)蝦群體幼體來源相同,其遺傳多樣性差異相對(duì)較小。包秀鳳(2014)針對(duì)親蝦源于正大集團(tuán)的養(yǎng)殖種苗群體開展微衛(wèi)星分析,發(fā)現(xiàn)受試群體的Na、Ne、He和PIC平均值分別為3.571、2.215、0.509和0.448,各項(xiàng)指標(biāo)均低于本研究結(jié)果。此外,于2018年5月在正大神灣蝦苗場(chǎng)采集的4個(gè)批次種苗群體(Stock 2、Stock 3、Stock 4和Stock 5)的I分別為1.5908、1.6092、1.4919和1.4874,其遺傳多樣性相近,且均高于其余5個(gè)群體(Stock 1、Stock 6、Stock 7、Stock 8和Stock 9),故推測(cè)遺傳多樣性水平能反映后續(xù)對(duì)蝦養(yǎng)殖生產(chǎn)性能的差異。

    Fst是衡量群體間遺傳分化程度的重要參數(shù)(黃小芳等,2020;彭敏等,2020)。Fst為0~0.05時(shí)表示群體間無遺傳分化,F(xiàn)st為0.05~0.15時(shí)表示群體間存在輕度遺傳分化,F(xiàn)st為0.15~0.25表示群體間存在中度遺傳分化,F(xiàn)st大于0.25時(shí)表示群體間存在高度遺傳分化。謝麗等(2009)研究指出Molokai、OI、SIS和Kona Bay等4個(gè)養(yǎng)殖凡納濱對(duì)蝦群體間的遺傳分化占15.83%;包秀鳳(2014)對(duì)國內(nèi)外7個(gè)選育凡納濱對(duì)蝦群體的遺傳分析結(jié)果表明,群體間的Fst平均值為0.2155;Zhang等(2014)研究發(fā)現(xiàn)7個(gè)國外引進(jìn)凡納濱對(duì)蝦群體間的遺傳分化水平為11.30%。在本研究中,9個(gè)凡納濱對(duì)蝦種苗群體間的Fst平均值為0.0818,即8.18%的變異來自于群體間,表現(xiàn)為中低水平的遺傳分化;但兩兩群體間的遺傳分化水平存在明顯差異,Stock 2、Stock 3、Stock 4和Stock 5等4個(gè)群體間及Stock 7群體與Stock 8群體間的Fst均小于0.05(0.0079~0.0371),即群體間無遺傳分化,而其他群體間表現(xiàn)為輕度遺傳分化,與各群體間的Nei’s遺傳距離基本對(duì)應(yīng)。一般情況下,群體間的遺傳距離越小,其親緣關(guān)系越近。基于Nei’s遺傳距離的聚類分析結(jié)果也表明,2018年5月采集自正大神灣蝦苗場(chǎng)的4個(gè)批次群體(Stock 2、Stock 3、Stock 4和Stock 5)間及2017年8月采集自海尚種苗培育基地的2個(gè)批次群體(Stock 7和Stock 8)間的親緣關(guān)系較接近,而較早采集自正大神灣蝦苗場(chǎng)的Stock 1群體與Stock 2~Stock 5群體的親緣關(guān)系相對(duì)較遠(yuǎn),其遺傳分化達(dá)中等水平。

    本研究利用微衛(wèi)星分子標(biāo)記分析不同批次進(jìn)口親蝦子一代種苗的遺傳變異特征,結(jié)果發(fā)現(xiàn)不同群體間存在遺傳多樣性差異,且表現(xiàn)出明顯的遺傳分化;由同一種苗培育(基地)企業(yè)提供的蝦苗在短期內(nèi)各群體間的遺傳分化保持在較低水平,但隨著時(shí)間的推移可能會(huì)發(fā)生較大變化;在同一時(shí)期由不同種苗培育(基地)企業(yè)培育獲得的蝦苗即使其親蝦來源相同,也會(huì)存在明顯的遺傳分化。這可能與近年來進(jìn)口親蝦子一代種苗在生長性能和養(yǎng)殖成活率等方面的性能變化有關(guān)聯(lián)。因此,有必要進(jìn)一步系統(tǒng)分析凡納濱對(duì)蝦進(jìn)口親蝦子一代種苗的遺傳變異狀況,并通過與之對(duì)應(yīng)的養(yǎng)殖實(shí)踐以獲得不同群體的養(yǎng)殖生產(chǎn)性能信息,建立各群體遺傳特征與養(yǎng)殖生產(chǎn)性能的關(guān)聯(lián)數(shù)據(jù)庫,為今后我國養(yǎng)殖對(duì)蝦優(yōu)質(zhì)種苗的大數(shù)據(jù)建設(shè)及其推廣應(yīng)用提供基礎(chǔ)數(shù)據(jù)支撐。

    4 結(jié)論

    不同批次凡納濱對(duì)蝦進(jìn)口親蝦子一代種苗群體的遺傳特征存在明顯差異,各群體間的遺傳距離隨著采樣時(shí)間和地點(diǎn)的改變而發(fā)生明顯變化,且這些遺傳背景差異可能會(huì)造成不同群體在養(yǎng)殖生產(chǎn)性能方面產(chǎn)生差異。

    參考文獻(xiàn):

    包秀鳳. 2014. 凡納濱對(duì)蝦選育群體遺傳多樣性分析[D]. 湛江:廣東海洋大學(xué). [Bao X F. 2014. Analysis of genetic diversity of selected stocks of Litopenaeus vannamei[D]. Zhanjiang:Guangdong Ocean University.] doi:10.7666/d.D522247.

    代平,孔杰,欒生. 2018. 我國凡納濱對(duì)蝦種質(zhì)資源引進(jìn)與分析[J]. 科學(xué)養(yǎng)魚,(1):3-5. [Dai P,Kong J,Luan S. 2018. Introduction and analysis of germplasm resources of Litopenaeus vannamei in China[J]. Scientific Fish Farming,(1):3-5.] doi:10.14184/j.cnki.issn1004-843x.2018.01. 002.

    黃小芳,唐章生,劉俊丹,張宏燕,鐘一治,盧智發(fā),侯樹鑒,王大鵬,陸專靈. 2020. 廣西不同地區(qū)克氏原螯蝦群體遺傳多樣性微衛(wèi)星分析[J]. 南方農(nóng)業(yè)學(xué)報(bào),51(2):437-444. [Huang X F,Tang Z S,Liu J D,Zhang H Y,Zhong Y Z,Lu Z F,Hou S J,Wang D P,Lu Z L. 2020. Genetic diversity microsatellite analysis of Procambarus clarkii po-pulations in different regions of Guangxi[J]. Journal of Southern Agriculture,51(2):437-444.] doi:10.3969/j.issn. 2095-1191.2020.02.025.

    李強(qiáng)勇,李旻,曾地剛,朱威霖,彭敏,楊春玲,劉青云,趙永貞,陳秀荔,陳曉漢. 2020. 凡納濱對(duì)蝦微衛(wèi)星分子標(biāo)記的開發(fā)及不同養(yǎng)殖家系遺傳多態(tài)性分析[J]. 南方農(nóng)業(yè)學(xué)報(bào),51(2):429-436. [Li Q Y,Li M,Zeng D G,Zhu W L,Peng M,Yang C L,Liu Q Y,Zhao Y Z,Chen X L,Chen X H. 2020. Development of microsatellite markers of Litopenaeus vannamei and genetic polymorphism analy-sis of different cultured families[J]. Journal of Southern Agriculture,51(2):429-436.] doi:10.3969/j.issn.2095-1191. 2020.02.024.

    劉連為,陳新軍,許強(qiáng)華,李偉文. 2014. 北太平洋柔魚微衛(wèi)星標(biāo)記的篩選及遺傳多樣性[J]. 生態(tài)學(xué)報(bào),34(23):6848-6854. [Liu L W,Chen X J,Xu Q H,Li W W. 2014. Isolation and genetic diversity of microsatellite DNA of Ommastrephes bartramii in the North Pacific Ocean[J]. Acta Ecologica Sinica,34(23):6848-6854.] doi:10.5846/stxb 201303030334.

    毛守康,馬愛軍,丁福紅,閆喜武,羅海忠,李偉業(yè),徐寶榮,王懷忠,王寶義,王廣寧. 2016. 梭魚(Liza haematocheila)4個(gè)野生地理群體遺傳多樣性的微衛(wèi)星分析[J]. 漁業(yè)科學(xué)進(jìn)展,37(2):68-75. [Mao S K,Ma A J,Ding F H,Yan X W,Luo H Z,Li W Y,Xu B R,Wang H Z,Wang B Y,Wang G Y. 2016. Analysis of genetic structures of four wild geographic populations of mullet Liza haematocheila by using microsatellite marker technique[J]. Progress in Fishery Sciences,37(2):68-75.] doi:10.11758/yykxjz.20150407001.

    彭敏,陳慧芳,李強(qiáng)勇,楊春玲,曾地剛,劉青云,趙永貞,陳曉漢,林勇,陳秀荔. 2020. 凡納濱對(duì)蝦連續(xù) 3 個(gè)世代選育群體的遺傳多樣性分析[J]. 南方農(nóng)業(yè)學(xué)報(bào),51(6):1442-1450. [Peng M,Chen H F,Li Q Y,Yang C L,Zeng D G,Liu Q Y,Zhao Y Z,Chen X H,Lin Y,Chen X L. 2020. Genetic diversity of three consecutive generations of Litopenaeus vannamei[J]. Journal of Southern Agriculture,51(6):1442-1450.] doi:10.3969/j.issn.2095-1191. 2020.06.026.

    孫效文,張曉鋒,趙瑩瑩,張研,賈智英,常玉梅,魯翠云,梁利群. 2008. 水產(chǎn)生物微衛(wèi)星標(biāo)記技術(shù)研究進(jìn)展及其應(yīng)用[J]. 中國水產(chǎn)科學(xué),15(4):689-703. [Sun X W,Zhang X F,Zhao Y Y,Zhang Y,Jia Z Y,Chang Y M,Lu C Y,Liang L Q. 2008. Development and application of microsatellite markers in aquatic species[J]. Journal of Fishery Sciences of China,15(4):689-703.] doi:10.3321/j.issn:1005-8737.2008.04.022.

    童馨,龔世圓,喻達(dá)輝,杜博,黃桂菊,李莉好,郭奕惠,李色東. 2007. 凡納濱對(duì)蝦不同世代生長性狀的變異[J]. 南方水產(chǎn)科學(xué),3(6):30-33. [Tong X,Gong S Y,Yu D H,Du B,Huang G J,Li L H,Guo Y H,Li S D. 2007. Variation of growth traits at different generations of Pacific white shrimp(Litopenaeus vannamei)[J]. South China Fisheries Science,3(6):30-33.] doi:10.3969/j.issn.2095-0780.2007.06.005.

    童馨,龔世圓,喻達(dá)輝,黃桂菊,杜博,李色東. 2009. 凡納濱對(duì)蝦(Litopenaeus vannamei)不同世代養(yǎng)殖群體的遺傳多樣性分析[J]. 海洋與湖沼,40(2):214-220. [Tong X,Gong S Y,Yu D H,Huang G J,Du B,Li S D. 2009. Genetic diversity of cultured Pacific white shrimp(Litopenaeus vannamei) stocks of different generations in China[J]. Oceanologia et Limnologia Sinica,40(2):214-220.] doi:10.3321/j.issn:0029-814X.2009.02.017.

    王麗花,曹煜成,李卓佳. 2012. 溶藻細(xì)菌控藻作用及其在對(duì)蝦養(yǎng)殖池塘中的應(yīng)用前景[J]. 南方水產(chǎn)科學(xué),8(4):76-82. [Wang L H,Cao Y C,Li Z J. 2012. Algae-constrai-ning effect of algae-lysing bacteria and its application prospect in shrimp ponds[J]. South China Fisheries Science,8(4):76-82.] doi:10.3969/j.issn.2095-0780.2012.04.012.

    王興強(qiáng),馬甡,董雙林. 2004. 凡納濱對(duì)蝦生物學(xué)及養(yǎng)殖生態(tài)學(xué)研究進(jìn)展[J]. 海洋湖沼通報(bào),(4):94-100. [Wang X Q,Ma S,Dong S L. 2004. Studies on the biology and cultural ecology of Litopenaeus vannamei:A review[J]. Transactions of Oceanology and Limnology,(4):94-100.] doi:10.3969/j.issn.1003-6482.2004.04.016.

    文國樑,曹煜成,徐煜,胡曉娟,徐武杰,李卓佳. 2015. 養(yǎng)殖對(duì)蝦肝胰腺壞死綜合癥研究進(jìn)展[J]. 廣東農(nóng)業(yè)科學(xué),42(11):118-123. [Wen G L,Cao Y C,Xu Y,Hu X J,Xu W J,Li Z J. 2015. Review on hepatopancreas necrosis syndrome of shrimp[J]. Guangdong Agricultural Sciences,42(11):118-123.] doi:10.3969/j.issn.1004-874X.2015.11. 022.

    吳怡迪,駱軒,楊章武,黃永春,游偉偉. 2016. 凡納濱對(duì)蝦選育系間雜交的生長性狀及遺傳多樣性分析[J]. 廈門大學(xué)學(xué)報(bào)(自然科學(xué)版),55(5):646-653. [Wu Y D,Luo X,Yang Z W,Huang Y C,You W W. 2016. Growth performance and genetic diversity analysis for hybrids between selective lines of Litopenaeus vannamei[J]. Journal of Xiamen University (Natural Science),55(5):646-653.] doi:10.6043/j.issn.0438-0479.201511010.

    謝麗,陳國良,葉富良,栗志民. 2009. 凡納濱對(duì)蝦4個(gè)選育群體遺傳多樣性的SSR分析[J]. 廣東海洋大學(xué)學(xué)報(bào),29(4):5-9. [Xie L,Chen G L,Ye F L,Li Z M. 2009. Genetic diversity of four selected stocks of Litopenaeus vannamei as revealed by SSR marker[J]. Journal of Guangdong Ocean University,29(4):5-9.] doi:10.3969/j.issn. 1673-9159.2009.04.002.

    趙志英,梁麗運(yùn),白麗蓉. 2018. 斑節(jié)對(duì)蝦3個(gè)野生群體遺傳多樣性的微衛(wèi)星標(biāo)記分析[J]. 熱帶海洋學(xué)報(bào),37(3):65-72. [Zhao Z Y,Liang L Y,Bai L R. 2018. Analysis of genetic diversity among three wild populations of Penaeus monodon using microsatellite marker[J]. Journal of Tropical Oceanography,37(3):65-72.] doi:10.11978/2017084.

    周大顏,張志新,黃彩林,招志杰,莫飛龍. 2019. 3種白鯽雜交子代的轉(zhuǎn)錄組學(xué)分析[J]. 南方農(nóng)業(yè)學(xué)報(bào),50(6):1328-1338. [Zhou D Y,Zhang Z X,Huang C L,Zhao Z J,Mo F L. 2019. Transcriptomic analysis of three Carassius auratus cuvieri hybrids[J]. Journal of Southern Agriculture,50(6):1328-1338.] doi:10.3969/j.issn.2095-1191.2019. 06.24.

    Alcivar-Warren A,Meehan-Meola D,Park S W,Xu Z K,Dela-ney M,Zuniga G. 2007a. ShrimpMap:A low-density,microsatellite-based linkage map of the pacific white leg shrimp,Litopenaeus vannamei:Identification of sex-linked markers in linkage group 4[J]. Journal of Shellfish Research,26(4):1259-1277. doi:10.2983/0730-8000(2007)26[1259:SALMLM]2.0.CO;2.

    Alcivar-Warren A,Song L S,Meehan-Meola D,Xu Z K,Xiang J H,Warren W. 2007b. Characterization and mapping of expressed sequence tags isolated from a subtrac-ted cDNA library of Litopenaeus vannamei injected with white spot syndrome virus[J]. Journal of Shellfish Resear-ch,26(4):1247-1258. doi:10.2983/0730-8000(2007)26[1247:CAMOES]2.0.CO;2.

    Cruz P,Mejia-Ruiz C H,Perez-Enriquez R,Ibarra A M. 2002. Isolation and characterization of microsatellites in Pacific white shrimp Penaeus(Litopenaeus) vannamei[J]. Molecular Ecology Notes,2(3):239-241. doi:10.1046/j.1471- 8286.2002.00211.x-i2.

    Excoffier L,Laval G,Schneider S. 2005. Arlequin (version 3.0):An integrated software package for population genetics data analysis[J]. Evolutionary Bioinformatics,1:47-50. doi:10.1177/117693430500100003.

    Garcia D K,Alcivar-Warren A. 2007. Characterization of 35 new microsatellite genetic markers for the pacific whiteleg shrimp,Litopenaeus vannamei:Their usefulness for studying genetic diversity of wild and cultured stocks,tracing pedigree in breeding programs,and linkage mapping[J]. Journal of Shellfish Research,26(4):1203-1216. doi:10.2983/0730-8000(2007)26[1203:CONMGM]2.0.CO;2.

    Garcia D K,Dhar A K,Alcivar-Warren A. 1996. Molecular analysis of a RAPD marker(B20) reveals two microsatellites and differential mRNA expression in Penaeus vannamei[J]. Molecular Marine Biology and Biotechnology,5(1):71-83.

    Jia Z Y,Sun X W,Liang L Q,Li D Y,Lei Q Q. 2006. Isolation and characterization of microsatellite markers from Pacific white shrimp(Litopenaeus vannamei)[J]. Molecular Ecology Notes,6(4):1282-1284. doi:10.1111/j.1471-8286.2006.01515.x.

    Meehan D,Xu Z K,Zuniga G,Alcivar-Warren A. 2003. High frequency and large number of polymorphic microsatellites in cultured shrimp,Penaeus(Litopenaeus) vannamei [Crustacea:Decapoda][J]. Marine Biotechnology,5:311-330. doi:10.1007/s10126-002-0092-z.

    Nagy S,Poczai P,Cernák I,Gorji A M,Heged?s G,Taller J. 2012. PICcalc:An online program to calculate polymorphic information content for molecular genetic studies[J]. Biochemical Genetics,50:670-672. doi:10.1007/s10528-012-9509-1.

    Qiu L,Chen M M,Wan X Y,Li C,Zhang Q L,Wang R Y,Cheng D Y,Dong X,Yang B,Wang X H,Xiang J H,Huang J. 2017. Characterization of a new member of Iridoviridae,Shrimp hemocyte iridescent virus(SHIV),found in white leg shrimp(Litopenaeus vannamei)[J]. Scienti-fic Reports,7:11834. doi:10.1038/s41598-017-10738-8.

    Shete S,Tiwari H,Elston R C. 2000. On estimating the he-terozygosity and polymorphism information content value[J]. Theoretical Population Biology,57(3):265-271. doi:10.1006/tpbi.2000.1452.

    Tamura K,Peterson D,Peterson N,Stecher G,Nei M,Kumar S. 2011. MEGA5:Molecular evolutionary genetics analysis using maximum likelihood,evolutionary distance,and maximum parsimony methods[J]. Molecular Biology and Evolution,28(10):2731-2739. doi:10.1093/molbev/msr121.

    Wang M S,Wang W J,Xiao G X,Liu K F,Hu Y L,Tian T,Kong J,Jin X S. 2016. Genetic diversity analysis of spaw-ner and recaptured populations of Chinese shrimp(Fenneropenaeus chinensis) during stock enhancement in the Bohai Bay based on an SSR marker[J]. Acta Oceanologica Sinica,35(8):51-56. doi:10.1007/s13131-016-0830-0.

    Yeh F C,Boyle T J B. 1997. Population genetic analysis of co-dominant and dominant markers and quantitative traits[J]. Belgian Journal of Botany,129(2):157.

    Zhang K,Wang W J,Li W Y,Zhang Q Q,Kong J. 2014. Analysis of genetic diversity and differentiation of seven stocks of Litopenaeus vannamei using microsatellite mar-kers[J]. Journal of Ocean University of China,13(4):647-656. doi:10.1007/s11802-014-2208-2.

    收稿日期:2020-05-26

    基金項(xiàng)目:中國水產(chǎn)科學(xué)研究院基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2020TD54);國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-48);廣東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系創(chuàng)新團(tuán)隊(duì)建設(shè)專項(xiàng)(2019KJ149)

    通訊作者:曹煜成(1979-),https://orcid.org/0000-0002-0995-1131,博士,研究員,主要從事水產(chǎn)健康養(yǎng)殖及養(yǎng)殖生態(tài)環(huán)境調(diào)控與修復(fù)研究工作,E-mail:cyc_715@163.com;楊鏗(1975-),https://orcid.org/0000-0002-6970-0668,副研究員,主要從事水產(chǎn)健康養(yǎng)殖技術(shù)及產(chǎn)品推廣應(yīng)用研究工作,E-mail:yangkeng66@163.com

    第一作者:徐煜(1987-),https://orcid.org/0000-0003-2364-6659,主要從事對(duì)蝦分子生理及健康養(yǎng)殖研究工作,E-mail:xuyublq@163.com

    1443501186395

    猜你喜歡
    遺傳變異
    葡萄醇酰基轉(zhuǎn)移酶編碼基因遺傳變異研究
    武漢市園科院“蓮基因組遺傳變異與表達(dá)數(shù)據(jù)庫”數(shù)據(jù)信息被國家基因庫收錄
    園林科技(2021年2期)2022-01-19 02:22:38
    先導(dǎo)編輯技術(shù)可編輯近90%的人類遺傳變異
    基于改進(jìn)遺傳變異算子的海島算法
    電子制作(2019年24期)2019-02-23 13:22:18
    寬葉藍(lán)靛果葉片不同海拔和郁閉度的遺傳變異研究
    森林工程(2018年3期)2018-06-26 03:40:46
    豬偽狂犬病毒的分離鑒定及其gE、TK全基因序列遺傳變異分析(內(nèi)文第133~144頁)圖版
    火力楠子代遺傳變異分析及優(yōu)良家系選擇
    江南油杉優(yōu)樹子代家系苗期遺傳變異與早期選擇
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    廣西鴨圓環(huán)病毒流行毒株遺傳變異分析
    欧美高清成人免费视频www| 日本av手机在线免费观看| 亚洲欧美日韩卡通动漫| 热re99久久精品国产66热6| 99久久精品一区二区三区| 人体艺术视频欧美日本| 精品久久久久久电影网| 大片免费播放器 马上看| 国产精品一二三区在线看| av专区在线播放| 免费观看av网站的网址| 久久久久久伊人网av| 色播亚洲综合网| 三级国产精品片| 小蜜桃在线观看免费完整版高清| 天堂中文最新版在线下载 | 欧美性感艳星| 少妇 在线观看| 成年女人在线观看亚洲视频 | 乱系列少妇在线播放| 亚洲天堂国产精品一区在线| 欧美三级亚洲精品| 午夜激情久久久久久久| 国产精品99久久99久久久不卡 | 狂野欧美激情性xxxx在线观看| 日本欧美国产在线视频| 久久精品熟女亚洲av麻豆精品| 在线天堂最新版资源| 少妇人妻 视频| 亚洲色图综合在线观看| 熟妇人妻不卡中文字幕| 国产欧美亚洲国产| 亚洲怡红院男人天堂| 国产综合精华液| 国产亚洲一区二区精品| 伊人久久精品亚洲午夜| 日韩不卡一区二区三区视频在线| 最后的刺客免费高清国语| 日本三级黄在线观看| 国产欧美亚洲国产| 看非洲黑人一级黄片| 搡老乐熟女国产| 热99国产精品久久久久久7| 涩涩av久久男人的天堂| 永久网站在线| 国产男人的电影天堂91| 大片免费播放器 马上看| 日本三级黄在线观看| h日本视频在线播放| 最新中文字幕久久久久| 最后的刺客免费高清国语| 极品少妇高潮喷水抽搐| 人体艺术视频欧美日本| 高清av免费在线| 国产精品嫩草影院av在线观看| 丝袜喷水一区| 日韩一本色道免费dvd| 国产亚洲午夜精品一区二区久久 | 亚洲精品日韩av片在线观看| 嫩草影院入口| 免费少妇av软件| 99精国产麻豆久久婷婷| 久久人人爽人人爽人人片va| 日韩 亚洲 欧美在线| 国产伦精品一区二区三区四那| 中文天堂在线官网| 久久影院123| 免费在线观看成人毛片| 日本熟妇午夜| 美女主播在线视频| videossex国产| 国产在线一区二区三区精| 成人高潮视频无遮挡免费网站| 亚洲国产高清在线一区二区三| 日产精品乱码卡一卡2卡三| eeuss影院久久| 国产乱人视频| 精品久久久久久久久亚洲| 免费av毛片视频| 久久久久久久午夜电影| 超碰av人人做人人爽久久| 丰满少妇做爰视频| 久久久久国产精品人妻一区二区| 欧美丝袜亚洲另类| 男人狂女人下面高潮的视频| 午夜精品国产一区二区电影 | 欧美三级亚洲精品| 亚洲av中文av极速乱| 菩萨蛮人人尽说江南好唐韦庄| 欧美少妇被猛烈插入视频| 亚洲国产最新在线播放| 久久精品熟女亚洲av麻豆精品| 永久网站在线| 18禁在线播放成人免费| 亚洲美女搞黄在线观看| 全区人妻精品视频| 亚洲国产欧美人成| 亚洲欧美日韩另类电影网站 | 人人妻人人澡人人爽人人夜夜| 国产精品av视频在线免费观看| 天天躁日日操中文字幕| 男女无遮挡免费网站观看| 久久热精品热| 亚洲aⅴ乱码一区二区在线播放| 亚洲av在线观看美女高潮| 99视频精品全部免费 在线| 国产精品人妻久久久影院| 在线免费十八禁| 亚洲av中文字字幕乱码综合| 视频区图区小说| 日韩欧美精品免费久久| 国产精品人妻久久久久久| 2021天堂中文幕一二区在线观| 一级毛片久久久久久久久女| 国产亚洲av片在线观看秒播厂| 肉色欧美久久久久久久蜜桃 | 一个人看视频在线观看www免费| 亚洲精品日韩av片在线观看| 久久国内精品自在自线图片| 97超碰精品成人国产| 黄片无遮挡物在线观看| 丰满少妇做爰视频| 18+在线观看网站| 久久99蜜桃精品久久| 91在线精品国自产拍蜜月| 丝瓜视频免费看黄片| 一级毛片黄色毛片免费观看视频| 国产高清国产精品国产三级 | 看十八女毛片水多多多| 日本猛色少妇xxxxx猛交久久| a级毛色黄片| 国产欧美亚洲国产| 精品午夜福利在线看| 狂野欧美白嫩少妇大欣赏| 亚洲无线观看免费| 国产成人91sexporn| 国产精品嫩草影院av在线观看| 国产又色又爽无遮挡免| 97人妻精品一区二区三区麻豆| 少妇猛男粗大的猛烈进出视频 | 91久久精品国产一区二区三区| 国产亚洲一区二区精品| 精品熟女少妇av免费看| 免费大片18禁| 又黄又爽又刺激的免费视频.| 日韩欧美 国产精品| 嫩草影院精品99| 中文字幕亚洲精品专区| 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区 | 久热久热在线精品观看| 大话2 男鬼变身卡| 啦啦啦中文免费视频观看日本| 久久精品国产a三级三级三级| 麻豆久久精品国产亚洲av| 激情 狠狠 欧美| 中文字幕亚洲精品专区| 男女啪啪激烈高潮av片| 欧美zozozo另类| 丝袜喷水一区| 婷婷色综合大香蕉| 久久久久精品性色| 久久国产乱子免费精品| 91在线精品国自产拍蜜月| 成年av动漫网址| 亚洲av成人精品一区久久| 国产淫片久久久久久久久| 黄色视频在线播放观看不卡| 国产老妇伦熟女老妇高清| 国产高清不卡午夜福利| 精品国产乱码久久久久久小说| 久久久久久久精品精品| 久久久久久久大尺度免费视频| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区| av一本久久久久| 国产片特级美女逼逼视频| 九九爱精品视频在线观看| 久久久精品免费免费高清| 少妇高潮的动态图| 亚洲天堂av无毛| 啦啦啦在线观看免费高清www| 国模一区二区三区四区视频| 欧美3d第一页| 国产成人一区二区在线| 91aial.com中文字幕在线观看| 亚洲精品久久午夜乱码| 狠狠精品人妻久久久久久综合| 日韩电影二区| 免费av不卡在线播放| 亚洲欧美精品专区久久| 国产伦精品一区二区三区视频9| 校园人妻丝袜中文字幕| 久久久久精品性色| 蜜桃久久精品国产亚洲av| 秋霞伦理黄片| 国产一区二区三区综合在线观看 | av国产久精品久网站免费入址| 国产精品嫩草影院av在线观看| 蜜桃亚洲精品一区二区三区| 午夜视频国产福利| 老司机影院毛片| 亚洲av成人精品一二三区| 69av精品久久久久久| 亚洲国产成人一精品久久久| 国产伦在线观看视频一区| 国产爱豆传媒在线观看| 国产精品久久久久久精品电影| 欧美日韩视频高清一区二区三区二| 少妇人妻精品综合一区二区| 日韩视频在线欧美| 波野结衣二区三区在线| 自拍偷自拍亚洲精品老妇| 狂野欧美激情性xxxx在线观看| 最近2019中文字幕mv第一页| 欧美日韩亚洲高清精品| 成人鲁丝片一二三区免费| 亚洲精品乱码久久久v下载方式| 欧美变态另类bdsm刘玥| 精品国产一区二区三区久久久樱花 | 日本欧美国产在线视频| 大又大粗又爽又黄少妇毛片口| 国产免费一级a男人的天堂| 大码成人一级视频| 亚洲国产精品成人久久小说| 午夜免费男女啪啪视频观看| 人妻夜夜爽99麻豆av| 高清日韩中文字幕在线| 日韩欧美精品免费久久| 伊人久久国产一区二区| 久久久久国产网址| 久久99蜜桃精品久久| 丰满乱子伦码专区| 大香蕉97超碰在线| 免费人成在线观看视频色| 精品人妻偷拍中文字幕| 免费不卡的大黄色大毛片视频在线观看| 神马国产精品三级电影在线观看| 精品少妇久久久久久888优播| 伊人久久国产一区二区| 国产女主播在线喷水免费视频网站| 91久久精品国产一区二区三区| 香蕉精品网在线| 插逼视频在线观看| 有码 亚洲区| 热re99久久精品国产66热6| 亚洲成色77777| 国产精品人妻久久久影院| 久久99热6这里只有精品| 国产精品一区二区三区四区免费观看| 99热国产这里只有精品6| 春色校园在线视频观看| 一级毛片aaaaaa免费看小| 2021天堂中文幕一二区在线观| 六月丁香七月| 免费少妇av软件| 天天躁夜夜躁狠狠久久av| 超碰av人人做人人爽久久| 欧美成人午夜免费资源| 国产高清不卡午夜福利| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 免费看光身美女| 日韩伦理黄色片| 大香蕉97超碰在线| 免费少妇av软件| 日韩av在线免费看完整版不卡| 国产成人精品婷婷| 99久国产av精品国产电影| av在线app专区| 亚洲经典国产精华液单| 久久99精品国语久久久| 国产av国产精品国产| 亚洲精品中文字幕在线视频 | 边亲边吃奶的免费视频| 精品久久国产蜜桃| 久久精品夜色国产| 久久久久久久久久成人| 国产一区二区在线观看日韩| 久久午夜福利片| 国产色爽女视频免费观看| 国产精品一区二区性色av| 九色成人免费人妻av| 亚洲欧美日韩另类电影网站 | 在线观看一区二区三区激情| www.色视频.com| 蜜桃亚洲精品一区二区三区| 日韩欧美精品v在线| 最近的中文字幕免费完整| 寂寞人妻少妇视频99o| 亚洲高清免费不卡视频| 少妇高潮的动态图| 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 男人舔奶头视频| 国产片特级美女逼逼视频| 99精国产麻豆久久婷婷| 卡戴珊不雅视频在线播放| 亚洲无线观看免费| av在线天堂中文字幕| 久热这里只有精品99| 成人国产av品久久久| av免费观看日本| 少妇 在线观看| 国产女主播在线喷水免费视频网站| 午夜老司机福利剧场| 国产成人精品福利久久| 校园人妻丝袜中文字幕| 免费看光身美女| 日韩av免费高清视频| 亚洲欧美日韩另类电影网站 | 精品久久久久久久人妻蜜臀av| 亚洲成人中文字幕在线播放| 色婷婷久久久亚洲欧美| 九草在线视频观看| 久久久久久久久大av| 国产精品久久久久久精品电影小说 | 99久久人妻综合| 国内精品美女久久久久久| 有码 亚洲区| 国产精品久久久久久久久免| 夜夜看夜夜爽夜夜摸| 国精品久久久久久国模美| 欧美+日韩+精品| 男女国产视频网站| 国产精品人妻久久久久久| 日本猛色少妇xxxxx猛交久久| 国产成人a区在线观看| 国产成人免费观看mmmm| 亚洲人成网站在线播| 少妇人妻一区二区三区视频| 蜜桃亚洲精品一区二区三区| 欧美日韩精品成人综合77777| 成人鲁丝片一二三区免费| 2022亚洲国产成人精品| 青春草国产在线视频| 在线观看av片永久免费下载| 黄色视频在线播放观看不卡| 久久精品人妻少妇| 麻豆乱淫一区二区| 久久久午夜欧美精品| 久久精品夜色国产| av在线观看视频网站免费| 最近2019中文字幕mv第一页| 国产黄色免费在线视频| 色视频www国产| 亚洲精品自拍成人| 欧美xxxx性猛交bbbb| 欧美精品一区二区大全| 最近的中文字幕免费完整| 久久久久久九九精品二区国产| 国产片特级美女逼逼视频| 永久网站在线| 国产人妻一区二区三区在| 成人毛片60女人毛片免费| 91久久精品国产一区二区成人| 国产色婷婷99| 成人国产麻豆网| 亚洲成人久久爱视频| 看十八女毛片水多多多| 综合色av麻豆| 亚洲国产成人一精品久久久| 午夜福利网站1000一区二区三区| 18禁在线播放成人免费| 久久精品国产亚洲av天美| 国产爱豆传媒在线观看| 国产成人a∨麻豆精品| 18+在线观看网站| 美女国产视频在线观看| 国产亚洲一区二区精品| 免费观看av网站的网址| 最近中文字幕高清免费大全6| 成人午夜精彩视频在线观看| 青春草亚洲视频在线观看| 国产精品不卡视频一区二区| 18禁在线播放成人免费| 人妻制服诱惑在线中文字幕| 嫩草影院入口| 肉色欧美久久久久久久蜜桃 | 久久久午夜欧美精品| 少妇熟女欧美另类| 日产精品乱码卡一卡2卡三| 建设人人有责人人尽责人人享有的 | 欧美xxⅹ黑人| 国产精品99久久99久久久不卡 | 欧美高清性xxxxhd video| 少妇被粗大猛烈的视频| 国产乱人偷精品视频| 一级毛片 在线播放| 免费黄频网站在线观看国产| 欧美精品国产亚洲| av专区在线播放| 国国产精品蜜臀av免费| 午夜视频国产福利| 欧美日韩视频高清一区二区三区二| 亚洲在久久综合| 日韩成人av中文字幕在线观看| 在线精品无人区一区二区三 | 狠狠精品人妻久久久久久综合| 亚洲欧洲国产日韩| 国产综合懂色| 有码 亚洲区| 免费看不卡的av| 久久综合国产亚洲精品| 国产伦在线观看视频一区| 亚洲精品久久午夜乱码| .国产精品久久| 日韩免费高清中文字幕av| 国产午夜福利久久久久久| 好男人在线观看高清免费视频| 午夜福利高清视频| 国产精品国产三级国产专区5o| 欧美日韩在线观看h| 在线亚洲精品国产二区图片欧美 | 午夜免费观看性视频| 国产亚洲最大av| 欧美老熟妇乱子伦牲交| 久久久久久久精品精品| av天堂中文字幕网| 老司机影院成人| 91久久精品电影网| 少妇熟女欧美另类| 久久精品国产自在天天线| 直男gayav资源| 日韩大片免费观看网站| 人妻 亚洲 视频| 麻豆国产97在线/欧美| 亚洲精品一二三| 最近中文字幕高清免费大全6| 白带黄色成豆腐渣| 亚洲四区av| 国产爱豆传媒在线观看| 成人鲁丝片一二三区免费| 日韩伦理黄色片| 国产日韩欧美亚洲二区| 九色成人免费人妻av| 欧美高清成人免费视频www| 午夜亚洲福利在线播放| 成人亚洲精品一区在线观看 | 性色av一级| 哪个播放器可以免费观看大片| 丝瓜视频免费看黄片| 爱豆传媒免费全集在线观看| 看黄色毛片网站| 精品少妇久久久久久888优播| 三级经典国产精品| 好男人在线观看高清免费视频| 男人舔奶头视频| 91精品国产九色| 亚洲国产精品专区欧美| 久久99蜜桃精品久久| 综合色丁香网| 成年女人在线观看亚洲视频 | 五月天丁香电影| 九九在线视频观看精品| 国产成人免费无遮挡视频| 色5月婷婷丁香| 有码 亚洲区| tube8黄色片| 少妇人妻 视频| 久久久久久久大尺度免费视频| 一区二区三区精品91| 一本久久精品| 亚洲成人久久爱视频| 九草在线视频观看| 观看美女的网站| 91精品一卡2卡3卡4卡| 国产真实伦视频高清在线观看| 水蜜桃什么品种好| 日本黄大片高清| 亚洲最大成人av| 欧美zozozo另类| 熟女av电影| 91久久精品国产一区二区三区| 久久久久久久久久久丰满| 日本午夜av视频| 秋霞在线观看毛片| 22中文网久久字幕| 一本一本综合久久| 熟女电影av网| 18禁裸乳无遮挡免费网站照片| 国产男女超爽视频在线观看| 久久影院123| 人妻一区二区av| 深爱激情五月婷婷| 午夜免费男女啪啪视频观看| 亚洲精品国产成人久久av| 少妇的逼水好多| 99热网站在线观看| 最近最新中文字幕大全电影3| 我的女老师完整版在线观看| 97在线人人人人妻| 99热国产这里只有精品6| 特大巨黑吊av在线直播| av国产免费在线观看| 国产中年淑女户外野战色| 午夜日本视频在线| 日韩 亚洲 欧美在线| 男女啪啪激烈高潮av片| 狂野欧美白嫩少妇大欣赏| 神马国产精品三级电影在线观看| 秋霞在线观看毛片| 国产精品福利在线免费观看| 美女高潮的动态| 免费看a级黄色片| 韩国高清视频一区二区三区| 汤姆久久久久久久影院中文字幕| 少妇熟女欧美另类| 2021天堂中文幕一二区在线观| 男人狂女人下面高潮的视频| 亚洲成人中文字幕在线播放| 中文字幕免费在线视频6| 视频区图区小说| 中文字幕人妻熟人妻熟丝袜美| 三级经典国产精品| 国产一区二区亚洲精品在线观看| 大又大粗又爽又黄少妇毛片口| 高清av免费在线| 尾随美女入室| 亚洲aⅴ乱码一区二区在线播放| 精品国产一区二区三区久久久樱花 | 下体分泌物呈黄色| 亚洲激情五月婷婷啪啪| 最近的中文字幕免费完整| 国产精品偷伦视频观看了| 91精品国产九色| 国产精品秋霞免费鲁丝片| 国产成人精品一,二区| 久久久久久久久久久免费av| 日本一二三区视频观看| 在线观看av片永久免费下载| 国产综合精华液| 免费av不卡在线播放| 黄色怎么调成土黄色| 三级男女做爰猛烈吃奶摸视频| 亚洲精品国产av成人精品| 国产精品一区二区三区四区免费观看| av国产久精品久网站免费入址| 亚洲欧美精品自产自拍| 日韩亚洲欧美综合| 天天躁夜夜躁狠狠久久av| 欧美丝袜亚洲另类| 婷婷色综合大香蕉| 亚洲天堂国产精品一区在线| 精品国产三级普通话版| 中文在线观看免费www的网站| 一级毛片 在线播放| 国产成人精品一,二区| 日本猛色少妇xxxxx猛交久久| 国产色婷婷99| 亚洲自偷自拍三级| 80岁老熟妇乱子伦牲交| 激情 狠狠 欧美| 亚洲在久久综合| 国产精品久久久久久av不卡| 九九在线视频观看精品| 好男人在线观看高清免费视频| 中文字幕av成人在线电影| 成年免费大片在线观看| 亚洲精品日韩av片在线观看| 99久久人妻综合| 大码成人一级视频| a级一级毛片免费在线观看| 国产精品精品国产色婷婷| 国产成人freesex在线| 男男h啪啪无遮挡| 国产熟女欧美一区二区| 国产精品熟女久久久久浪| 亚洲av中文av极速乱| 国产一级毛片在线| 少妇人妻 视频| 国产免费又黄又爽又色| 听说在线观看完整版免费高清| 黄色配什么色好看| 好男人视频免费观看在线| 欧美xxⅹ黑人| 性色av一级| 久久久成人免费电影| 日韩电影二区| 超碰av人人做人人爽久久| 波野结衣二区三区在线| 下体分泌物呈黄色| 深夜a级毛片| 午夜精品国产一区二区电影 | 中文字幕免费在线视频6| 综合色av麻豆| 天堂中文最新版在线下载 | 街头女战士在线观看网站| 中文字幕久久专区| 亚洲婷婷狠狠爱综合网| 亚洲成人精品中文字幕电影| 成年女人在线观看亚洲视频 | 亚洲av中文av极速乱| 真实男女啪啪啪动态图| 欧美zozozo另类| 亚洲精品久久久久久婷婷小说| av在线老鸭窝| 91久久精品国产一区二区成人| 婷婷色麻豆天堂久久| 国产 一区 欧美 日韩| 啦啦啦中文免费视频观看日本| 日韩av在线免费看完整版不卡| 国产精品秋霞免费鲁丝片| 亚洲精华国产精华液的使用体验| 国产亚洲精品久久久com| 日韩精品有码人妻一区| 日韩在线高清观看一区二区三区| 日本一二三区视频观看| 又大又黄又爽视频免费| 91狼人影院| 中文乱码字字幕精品一区二区三区| 99久久精品热视频| 精品少妇黑人巨大在线播放| 狂野欧美白嫩少妇大欣赏| 国产一区亚洲一区在线观看|