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

    基于生長和抗逆功能基因SNP分子標(biāo)記的凡納濱對蝦野生及選育群體遺傳多樣性分析

    2020-02-22 13:17申淑慧戴習(xí)林
    關(guān)鍵詞:遺傳多樣性

    申淑慧 戴習(xí)林

    摘要:【目的】分析凡納濱對蝦野生群體與選育群體的遺傳多樣性及遺傳差異,為進(jìn)一步改良優(yōu)化凡納濱對蝦種質(zhì)提供參考依據(jù)?!痉椒ā客ㄟ^凡納濱對蝦轉(zhuǎn)錄組測序數(shù)據(jù)挖掘候選SNP位點(diǎn),設(shè)計(jì)擴(kuò)增引物,采用10尾野生凡納濱對蝦進(jìn)行直接測序驗(yàn)證,篩選獲得11對有效引物及24個(gè)位于生長和抗逆功能基因上的SNPs位點(diǎn),并用于分析2個(gè)野生群體(Y1和Y2)和2個(gè)選育群體(F1和F2)的遺傳多樣性?!窘Y(jié)果】凡納濱對蝦野生群體的觀測雜合度(Ho)、期望雜合度(He)及多態(tài)信息含量(PIC)等遺傳參數(shù)平均值明顯高于選育群體。在2個(gè)野生群體中,Y2的各項(xiàng)遺傳參數(shù)平均值均高于Y1;在2個(gè)選育群體中,F(xiàn)2的各項(xiàng)遺傳參數(shù)平均值均高于F1,且在CATL、ANT、Tryptase、Crustin-like和Penaeidin4c等功能基因的某些SNP位點(diǎn)上,2個(gè)選育群體的各項(xiàng)遺傳參數(shù)平均值為0。4個(gè)凡納濱對蝦群體中均存在SNP位點(diǎn)不同程度偏離Hardy-Weinberg平衡現(xiàn)象。野生群體Y1與選育群體F1間的近交系數(shù)(Fis)最高(0.0875),且基因流(Nm)也較高(4.2728);野生群體Y2與選育群體F1間的遺傳分化指數(shù)(Fst)最高(0.0947),但Nm最低(2.3886),表明二者間的遺傳背景較遠(yuǎn)。4個(gè)凡納濱對蝦群體的遺傳變異主要來源于群體內(nèi)(85.15%),僅有14.85%變異來源于群體間。4個(gè)凡納濱對蝦群體的遺傳相似性系數(shù)為0.9351~0.9747,遺傳距離為0.0256~0.0671;基于凡納濱對蝦群體遺傳距離的UPGMA聚類分析結(jié)果顯示,Y1、F1和F2聚為一支,而Y2單獨(dú)聚為一支?!窘Y(jié)論】凡納濱對蝦野生群體在抗逆功能基因SNP分子標(biāo)記中的遺傳多樣性顯著高于選育群體,且野生群體與選育群體間的遺傳距離較遠(yuǎn),具有較高的遺傳改良潛力。在實(shí)際生產(chǎn)中,可通過雜交措施將野生群體的優(yōu)良抗逆基因引入選育群體,提高凡納濱對蝦選育群體的免疫及抗病能力。

    關(guān)鍵詞: 凡納濱對蝦;野生群體;選育群體;SNP分子標(biāo)記;生長功能基因;抗逆功能基因;遺傳多樣性

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

    Genetic diversity analysis in wild and selected populations of Penaeus vannamei based on SNP markers in growth

    and stress resistance genes

    SHEN Shu-hui1,2,3, DAI Xi-lin1,2,3*

    (1National Demonstration Center for Experimental Fisheries Science Education(Shanghai Ocean University), Shanghai 201306, China; 2Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture and Rural Affairs (Shanghai Ocean University), Shanghai ?201306, China; 3Shanghai Collaborative Innovation for Aquatic Animal

    Genetics and Breeding(Shanghai Ocean University), Shanghai ?201306, China)

    Abstract:【Objective】The genetic diversity and genetic differences between wild populations and selected populations of Litopenaeus vannamei were analyzed to provide reference for improvement and optimization of germplasm resources. 【Method】Single nucleotide polymorphisms(SNPs) found from amplification primers, which designed with high-throughput sequencing of L. vannamei, were verified using ten wild L. vannamei by sequencing. Finally, 11 primers and 24 identified SNPs in growth and stress resistance genes, were selected to analyze genetic diversity of Y1, Y2, which were wild populations, and F1, F2, which were selected populations. 【Result】The results showed that the genetic parameters such as ?observed heterozygosity(Ho), expected ?heterozygosity(He) and polymorphism information content(PIC) in Y1 and Y2 were higher than those in F1 and F2. Generally, all the genetic parameters in Y2 were higher than Y1. Similarly, in two selected families, genetic parameters of F2 were higher than F1. In addition, the genetic parameters in the two selected populations were 0 at some SNP loci of functional genes such as CATL, ANT, Tryptase, Crustin-like and Penaeidin4c. Deviation from Hardy-Weinberg equilibrium existed in all four populations differently. The coefficient of inbreeding (Fis) was the highest(0.0875) and the gene flow(Nm) was also high (4.2728) between wild population Y1 and selected population F1. The genetic differentiation index(Fst) between Y2 and F1 was the highest(0.0947), but Nm was the lowest (2.3886), indicating a distant genetic relationship between the two populations. AMOVA analysis found that 85.15% genetic variation of four populations was mainly within the population, and only 14.85% was found between populations.The genetic similarity coefficients of the four populations ranged from 0.9351 to 0.9747, and the genetic distance varied from 0.0256 to 0.0671. UPGMA clustering analysis showed that Y1, F1 and F2 clustered into a cluster, while Y2 clustered into a single cluster. 【Conclusion】Generally, in wild populations, genetic diversity of SNP markers in stress resistance function gene is significantly higher than that of selected populations. And wild populations are far from the selected populations, which shows a high potential for breeding. Therefore, it is necessary that the excellent stress resistance genes from wild populations should be introduced into the selected population through hybridization in production process, to improve the immunity and disease resistance of the selected population.

    Key words: Litopenaeus vannamei; wild population; ?selected populations; SNP markers; growth functional genes; ?stress resistance genes;genetic diversity

    Foundation item: Special Project of Shanghai Shrimp Industry Technology System Construction (HNKGZ〔2014〕7-1-11)

    0 引言

    【研究意義】凡納濱對蝦(Litopenaeus vannamei)又稱南美白對蝦,原產(chǎn)于太平洋沿岸,因具有含肉率高、抗逆性強(qiáng)、生長快、餌料需求低、便于活蝦運(yùn)輸及耐高密度和低鹽度養(yǎng)殖等優(yōu)點(diǎn),已發(fā)展成為世界上最主要的蝦類養(yǎng)殖品種,約占蝦類總產(chǎn)量的85%(楊銘等,2017;李強(qiáng)勇等,2020)。隨著對蝦養(yǎng)殖業(yè)的不斷發(fā)展,因病害頻發(fā)和種質(zhì)退化造成的對蝦抗病力下降、生產(chǎn)周期延長及產(chǎn)量下降等問題日益凸顯,給對蝦養(yǎng)殖帶來巨大經(jīng)濟(jì)損失,嚴(yán)重制約了凡納濱對蝦養(yǎng)殖業(yè)的持續(xù)健康發(fā)展(戴習(xí)林等,2018;蔣葛等,2019)。因此,培育抗逆性強(qiáng)且生長快的優(yōu)良對蝦養(yǎng)殖品種,是對蝦養(yǎng)殖業(yè)今后發(fā)展的重要任務(wù)(楊銘等,2017;葉寧等,2017)。其中,引進(jìn)原種或改良現(xiàn)有養(yǎng)殖品種是恢復(fù)和優(yōu)化凡納濱對蝦種質(zhì)的有效途徑,但在遺傳改良過程中必須掌握已有群體的種質(zhì)資源現(xiàn)狀及遺傳背景?!厩叭搜芯窟M(jìn)展】對蝦的遺傳育種工作起步相對較晚,但隨著PCR和分子標(biāo)記技術(shù)的不斷成熟及推廣應(yīng)用,近年來對蝦遺傳育種研究取得了長足進(jìn)展。于洋(2014)、郝敏(2016)、李東宇等(2016)、王全超(2017)、Yu等(2019)先后從凡納濱對蝦遺傳圖譜構(gòu)建、數(shù)量性狀定位、遺傳多樣性分析、親權(quán)鑒定及分子輔助育種等方面開展研究,研究結(jié)果為凡納濱對蝦的良種選育提供了參考依據(jù)。目前,針對凡納濱對蝦的遺傳學(xué)研究主要是利用等位酶、單拷貝核DNA、線粒體DNA、限制性片段長度多態(tài)性(RFLP)、擴(kuò)增片段長度多態(tài)性(AFLP)及微衛(wèi)星(SSR)等分子標(biāo)記進(jìn)行群體遺傳背景分析(張靈俠等,2008;葉寧等,2017;陳錦豪等,2019),其中又以SSR分子標(biāo)記在凡納濱對蝦遺傳多樣性研究中的應(yīng)用最廣泛。Rezaee等(2016)、劉洪濤等(2018)、方振朋等(2019)、黃小帥等(2019)分別研發(fā)出SSR分子標(biāo)記并應(yīng)用于凡納濱對蝦群體遺傳多樣性研究,但SSR分子標(biāo)記也存在一些缺陷,包括可能存在無效等位基因、無法存在于功能基因中及費(fèi)時(shí)費(fèi)力等。單核苷酸多態(tài)性標(biāo)記(Single nucleotide polymorphism,SNP)作為第三代分子標(biāo)記,屬于共顯性標(biāo)記,雖然每個(gè)標(biāo)記的基礎(chǔ)信息量較少,但能穩(wěn)定、均勻、廣泛地分布于整個(gè)基因組中,且可能出現(xiàn)在功能基因內(nèi)部,更有利于開展生長、抗逆及繁殖等相關(guān)功能基因的研究,是分析遺傳變異的有效分子工具(張曉萌等,2013;王全超,2017;陳靜等,2018)。有關(guān)SNP分子標(biāo)記在水產(chǎn)動(dòng)物群體遺傳學(xué)方面的應(yīng)用研究已有較多報(bào)道,Agarwal等(2016)通過高通量測序獲得大量羅氏沼蝦SNP分子標(biāo)記,并用于分析印度不同區(qū)域羅氏沼蝦群體的遺傳背景;于磊等(2018)基于酶切測序技術(shù)開發(fā)獲得38個(gè)多態(tài)性較高的SNP分子標(biāo)記,并以此分析遼寧丹東和福建三沙2個(gè)日本鰻鱺群體的遺傳結(jié)構(gòu)及其本地適應(yīng)性;Thongda等(2018)研究表明,開發(fā)獲得的58個(gè)東方牡蠣SNP分子標(biāo)記均可用于遺傳分化、群體結(jié)構(gòu)分析及親緣關(guān)系鑒定等相關(guān)研究?!颈狙芯壳腥朦c(diǎn)】SNP分子標(biāo)記是一種研究物種遺傳變異的理想分子標(biāo)記,但至今鮮見應(yīng)用于凡納濱對蝦遺傳多樣性分析的研究報(bào)道?!緮M解決的關(guān)鍵問題】基于凡納濱對蝦轉(zhuǎn)錄組數(shù)據(jù),挖掘與生長和抗逆相關(guān)候選基因中的SNP位點(diǎn),并以此分析凡納濱對蝦野生群體與選育群體的遺傳多樣性及遺傳差異,為進(jìn)一步改良優(yōu)化凡納濱對蝦種質(zhì)提供參考依據(jù)。

    1 材料與方法

    1. 1 樣品采集及DNA提取

    凡納濱對蝦樣本為2個(gè)厄瓜多爾野生群體和2個(gè)選育群體。其中,厄瓜多爾野生群體為2016和2019年從厄瓜多爾引進(jìn)的野生凡納濱對蝦,捕撈自瓜亞基樂灣水域,分別命名為Y1和Y2;選育群體是以美國關(guān)島大學(xué)引進(jìn)家系選育獲得的子二代,分別命名為F1和F2。從各群體中隨機(jī)取30尾樣本,無水酒精固定后以天根海洋動(dòng)物組織基因組DNA提取試劑盒提取其DNA,-20 ℃冰箱中保存?zhèn)溆谩?/p>

    1. 2 SNP位點(diǎn)發(fā)掘與引物篩選

    根據(jù)農(nóng)業(yè)農(nóng)村部淡水水產(chǎn)種質(zhì)資源重點(diǎn)實(shí)驗(yàn)室已有轉(zhuǎn)錄組數(shù)據(jù)中的基因注釋,挖掘與凡納濱對蝦生長和抗逆相關(guān)候選基因的預(yù)測SNP位點(diǎn),包括胰蛋白酶(Tryptase)、糜蛋白酶(Chymotrypsin)、保幼激素(JH)(杜育哲等,2001)、腺苷酸轉(zhuǎn)移酶(ANT)、性腺抑制激素(GIH)(邱慶慶,2019)、組織蛋白酶L(CATL)、有絲分裂原活化蛋白激酶(MAPK)、Alpha淀粉酶(AMY)(劉立鶴等,2006)、熱休克蛋白(HSP90)(夏西超等,2012)、ATP轉(zhuǎn)運(yùn)酶、對蝦素(Penaeidin)(杜志強(qiáng),2012,2013)、酚氧化酶(Phenoloxidase,PO)(樊廷俊,2002)、β-葡聚糖結(jié)合蛋白(BGBP) (胡世康等,2018)、Toll受體(黃旭雄等,2012)和Crustin-like(張艷艷等,2012)等功能基因。然后選用10尾野生凡納濱對蝦樣本進(jìn)行驗(yàn)證,在挑選基因序列中的預(yù)測SNP位點(diǎn)兩端設(shè)計(jì)引物,PCR擴(kuò)增結(jié)果送至上海邁浦生物科技有限公司測序。測序結(jié)果使用BioEdit進(jìn)行分析,讀取轉(zhuǎn)錄本預(yù)測位置的基因型,在10尾野生凡納濱對蝦中其雜合基因型概率達(dá)20%的SNP位點(diǎn)所對應(yīng)的引物即為有效引物,可用于后續(xù)的凡納濱對蝦群體基因型檢測。

    1. 3 統(tǒng)計(jì)分析

    獲得4個(gè)凡納濱對蝦群體的基因分型數(shù)據(jù)后,使用PopGen32計(jì)算群體的觀測雜合度(Ho)、期望雜合度(He)、Nei基因多樣性指數(shù)(H)、Hardy-Weinberg平衡偏離指數(shù)(D)、近交系數(shù)(Fis)、遺傳分化指數(shù)(Fst)、遺傳相似性系數(shù)及遺傳距離等(戴習(xí)林等,2017);利用Cervus 3.0計(jì)算SNP位點(diǎn)的多態(tài)信息含量(PIC),剔除多態(tài)性過低的SNP位點(diǎn);采用Arlequin V3.5分析Fst和基因流(Nm),并以MEGA 5.0進(jìn)行UPGMA聚類分析。

    2 結(jié)果與分析

    2. 1 SNP位點(diǎn)發(fā)掘及引物篩選結(jié)果

    從轉(zhuǎn)錄組中挖掘獲得70個(gè)SNPs位點(diǎn),設(shè)計(jì)對應(yīng)的引物并在10尾野生蝦中進(jìn)行PCR擴(kuò)增,結(jié)果顯示有46對引物能成功擴(kuò)增出目的片段。圖1為部分?jǐn)U增結(jié)果的凝膠電泳圖,將擴(kuò)增獲得且片段大小在預(yù)期范圍內(nèi)的單一明亮條帶(編號33、34、35和39)送至上海邁浦生物科技有限公司測序。根據(jù)測序結(jié)果讀取峰圖,并統(tǒng)計(jì)基因型,發(fā)現(xiàn)SNP雜合基因型概率達(dá)20%的有效引物共計(jì)15對,包含33個(gè)SNPs位點(diǎn)。采用15對篩選出的有效引物對4個(gè)凡納濱對蝦群體進(jìn)行PCR擴(kuò)增,分別統(tǒng)計(jì)33個(gè)SNPs位點(diǎn)的基因型,并剔除在4個(gè)凡納濱對蝦群體中突變率均低于15%的SNP位點(diǎn),最終獲得11對有效引物(表1),包含24個(gè)SNPs位點(diǎn)。其中,HSP90、GIH、CATL、ANT、Tryptase和Chymotrypsin為生長功能基因,含有9個(gè)SNPs位點(diǎn);Toll、Crustin-like、Penaeidin4c和PO為抗逆功能基因,含有15個(gè)SNPs位點(diǎn)。

    2. 2 遺傳多態(tài)性分析結(jié)果

    基于9個(gè)生長功能基因SNPs位點(diǎn)和15個(gè)抗逆功能基因SNPs位點(diǎn),研究4個(gè)凡納濱對蝦群體(2個(gè)野生群體和2個(gè)選育群體)的遺傳多樣性。結(jié)果(表2)顯示,24個(gè)SNP分子標(biāo)記在4個(gè)凡納濱對蝦群體中的PIC平均值分別為Y1:0.2587、Y2:0.2762、F1:0.1341和F2:0.1513,2個(gè)野生群體明顯高于2個(gè)選育群體。其中,Y1和Y2屬于中度多態(tài)水平(0.25

    總體上,2個(gè)野生群體(Y1和Y2)的各項(xiàng)遺傳參數(shù)平均值明顯高于2個(gè)選育群體(F1和F2)。在2個(gè)野生群體中,Y2的各項(xiàng)遺傳參數(shù)平均值均高于Y1,在生長功能基因中Y1的PIC和He平均值高于Y2,但在抗逆功能基因中Y2高于Y1。在2個(gè)選育群體中,F(xiàn)2的各項(xiàng)遺傳參數(shù)平均值均高于F1,在抗逆功能基因中F1和F2基本相同,但在生長功能基因中二者存在明顯差異,表現(xiàn)為F2高于F1。在CATL、ANT、Tryptase、Crustin-like和Penaeidin4c等功能基因的某些SNP位點(diǎn)上,2個(gè)選育群體(F1和F2)的各項(xiàng)遺傳參數(shù)平均值為0,可能是雜合子缺失所引起。4個(gè)凡納濱對蝦群體中均存在SNP位點(diǎn)不同程度偏離Hardy-Weinberg平衡現(xiàn)象(P<0.05為顯著偏離,P<0.01為極顯著偏離)。其中,Y1群體的SNP19和SNP20位點(diǎn)顯著偏離Hardy-Weinberg平衡,處于Penaeidin4c基因上;SNP12和SNP13位點(diǎn)極顯著偏離Hardy-Weinberg平衡,位于Crustin-like基因上;Y2群體的SNP3和SNP7位點(diǎn)顯著偏離Hardy-Weinberg平衡,分別位于GIH和Chymotrypsin基因上;F1群體中除突變?nèi)笔稽c(diǎn)外,僅SNP3位點(diǎn)顯著偏離Hardy-Weinberg平衡,位于GIH基因上;F2群體中除突變?nèi)笔稽c(diǎn)外,僅SNP6和SNP7位點(diǎn)顯著偏離Hardy-Weinberg平衡,均位于Chymotrypsin基因上。

    2. 3 遺傳分化分析結(jié)果

    4個(gè)凡納濱對蝦群體在24個(gè)SNPs位點(diǎn)上的Fis平均值為-0.0224,且負(fù)值較多(表3),表明不同群體間的近交程度較低。不同群體間各SNP位點(diǎn)的Fst變化范圍為0.0050~0.2754,說明4個(gè)凡納濱對蝦群體在24個(gè)SNPs位點(diǎn)上表現(xiàn)出明顯的遺傳分化。4個(gè)凡納濱對蝦群體相互間的Fis、Fst和Nm詳見表4。其中,Y1與F1間的Fis最高(0.0875),且Nm也較高(4.2728);Y2與F1間的Fst最高(0.0947),但其Nm最低(2.3886),表明Y2與F1間的遺傳背景較遠(yuǎn)?;赟NP分子標(biāo)記的凡納濱對蝦群體遺傳變異分子生物學(xué)方差分析結(jié)果(表5)顯示,4個(gè)凡納濱對蝦群體的遺傳變異大部分來源于群體內(nèi)(85.15%),僅有14.85%變異來源于群體間。

    2. 4 遺傳距離及聚類分析結(jié)果

    由表6可知,4個(gè)凡納濱對蝦群體的遺傳相似性系數(shù)為0.9351~0.9747,遺傳距離為0.0256~0.0671。其中,F(xiàn)1與F2的遺傳相似性系數(shù)最大(0.9747),Y2與F1的遺傳相似性系數(shù)最小(0.9351);F1與F2的遺傳距離最?。?.0256),Y2與F1的遺傳距離最大(0.0671)?;诜布{濱對蝦群體間的遺傳距離,采用MEGA 5.0進(jìn)行UPGMA聚類分析,結(jié)果(圖2)顯示,Y1、F1和F2聚為一支,而Y2單獨(dú)聚為一支。

    3 討論

    SNP可通過改變堿基引起所處基因的表達(dá)變化,導(dǎo)致其編碼蛋白發(fā)生結(jié)構(gòu)功能變異,進(jìn)而影響表型性狀。目前,開展水產(chǎn)動(dòng)物遺傳多樣性分析時(shí)有關(guān)分子標(biāo)記的選擇尚存在隨機(jī)性,且主要針對相關(guān)功能基因的分子標(biāo)記進(jìn)行研究(崔蕾等,2012;孫成飛等,2015;武菲等,2016),利用SNP分子標(biāo)記進(jìn)行遺傳多樣性的研究相對較少。本研究基于凡納濱對蝦的生長和抗逆相關(guān)功能基因,并通過轉(zhuǎn)錄組注釋篩選其中的候選SNP位點(diǎn),最終篩選出24個(gè)SNPs位點(diǎn)(9個(gè)SNPs位點(diǎn)位于生長功能基因上,15個(gè)SNPs位點(diǎn)位于抗逆功能基因上)用于凡納濱對蝦野生和選育群體遺傳多樣性分析。生長和抗逆相關(guān)功能基因中的SNP位點(diǎn)差異會(huì)影響其所處功能基因表達(dá),相對于其他中性位點(diǎn)更容易造成機(jī)體生長性狀及抗病能力的差異,若將其應(yīng)用到群體遺傳背景分析中,則有望從DNA水平了解群體生長性狀及抗病能力的遺傳差異,更有利于后續(xù)將分子標(biāo)記與性狀聯(lián)系起來,輔助良種選育工作的開展。

    本研究通過轉(zhuǎn)錄組數(shù)據(jù)結(jié)合直接測序挖掘出24個(gè)功能基因相關(guān)SNPs位點(diǎn)(9個(gè)生長功能基因SNPs位點(diǎn),15個(gè)抗逆功能基因SNPs位點(diǎn)),并應(yīng)用于凡納濱對蝦野生和選育群體遺傳多樣性分析,結(jié)果表明,野生群體Y1和Y2的遺傳多樣性水平顯著高于選育群體F1和F2,與董在杰等(2018)在鯉魚、王軍等(2018)在中國明對蝦、王沈同等(2018)在草魚中的研究結(jié)果相似,可能是由于選育過程中封閉群體的親本數(shù)量限制、近親交配、遺傳漂變及對養(yǎng)殖環(huán)境適應(yīng)等因素所引起。2個(gè)野生群體具有更豐富的遺傳多樣性,說明其進(jìn)化潛力更高,即具有較高的遺傳育種改良潛力和適應(yīng)生存能力。此外,在2個(gè)野生群體中,Y2的各項(xiàng)遺傳參數(shù)平均值均高于Y1,在生長功能基因中Y1的PIC和He平均值高于Y2,在抗逆功能基因中則Y2高于Y1,究其原因可能是不同野生種群在自然環(huán)境中遷移或變異所造成。在2個(gè)選育群體中,F(xiàn)2的各項(xiàng)遺傳參數(shù)平均值均高于F1,在抗逆功能基因中F1和F2基本相同,但在生長功能基因中二者存在明顯差異,表現(xiàn)為F2高于F1,說明在選擇育種中F2較F1更具優(yōu)勢。在CATL、ANT、Tryptase、Crustin-like和Penaeidin4c等功能基因的某些SNP位點(diǎn)上,2個(gè)選育群體(F1和F2)的各項(xiàng)遺傳參數(shù)平均值為0,推測F1和F2存在雜合子缺失現(xiàn)象,與頡曉勇等(2008)、葉寧等(2017)的研究結(jié)果相似,是由于人工選育過程中的近親交配現(xiàn)象或?qū)δ骋惶囟ɑ蛐蛡€(gè)體的篩選,造成某些基因從該基因庫中流失。凡納濱對蝦遺傳多樣性降低,極易造成種質(zhì)資源退化、生長滯緩,進(jìn)而嚴(yán)重影響其養(yǎng)殖業(yè)的健康發(fā)展(李強(qiáng)勇等,2020)。因此,在凡納濱對蝦選育過程中,需及時(shí)了解其遺傳多樣性水平,并在此基礎(chǔ)上進(jìn)行科學(xué)管理,才能在遺傳改良的同時(shí)維持選育群體具有較高的遺傳多樣性水平,或從野生群體引種,進(jìn)一步豐富凡納濱對蝦選育群體的遺傳多樣性。

    從本研究的遺傳分化及UPGMA聚類分析可看出,凡納濱對蝦野生群體Y1與2個(gè)選育群體(F1和F2)聚為一支,且Nm和Fis均較高,與三者間遺傳距離較近的結(jié)果相對應(yīng);而凡納濱對蝦野生群體Y2與2個(gè)選育群體(F1和F2)間的Fst較高,且Nm較低,與三者間遺傳距離最遠(yuǎn)的結(jié)果相對應(yīng),表明野生群體Y2在家系選擇育種中更具引種優(yōu)勢和選育潛力。分子生物學(xué)方差分析結(jié)果也顯示,4個(gè)凡納濱對蝦群體的遺傳變異主要來源于群體內(nèi)(85.15%),僅有14.85%變異來源于群體間,與童馨等(2009)、Rezaee等(2016)的研究結(jié)果類似,說明不同凡納濱對蝦群體間相似性較小且基因交流較少,仍具有較高的選育潛力。本研究還發(fā)現(xiàn),在生長功能基因SNP位點(diǎn)遺傳多樣性分析中,F(xiàn)2的遺傳多樣性水平明顯高于F1,且在相同養(yǎng)殖條件下F2的生長性能也優(yōu)于F1,說明F2無論是在遺傳上還是生長性狀上均優(yōu)于F1,在后續(xù)執(zhí)行選育計(jì)劃時(shí)應(yīng)優(yōu)先考慮。但從抗逆性狀的遺傳背景來看,2個(gè)選育群體(F1和F2)均不具備抗逆優(yōu)勢,在傳染病頻發(fā)的養(yǎng)殖現(xiàn)狀中極易被淘汰。因此,在實(shí)際生產(chǎn)中非常有必要通過選育群體與野生群體雜交以提高凡納濱對蝦后代的抗逆能力,可采取抗逆性狀高的野生群體Y2與選育群體相雜交的措施,將Y2的優(yōu)良抗逆基因引入選育群體,優(yōu)化其種質(zhì)資源,提高凡納濱對蝦選育群體的免疫及抗病能力。

    4 結(jié)論

    凡納濱對蝦野生群體在抗逆功能基因SNP分子標(biāo)記中的遺傳多樣性顯著高于選育群體,且野生群體與選育群體間的遺傳距離較遠(yuǎn),具有較高的遺傳改良潛力。在實(shí)際生產(chǎn)中,可通過雜交措施將野生群體的優(yōu)良抗逆基因引入家系選育群體,優(yōu)化其種質(zhì)資源,以提高凡納濱對蝦選育群體的免疫及抗病能力。

    參考文獻(xiàn):

    陳錦豪,鄭錦濱,王攀攀,李天驕,毛勇,蘇永全,王軍. 2019. 凡納濱對蝦養(yǎng)殖親本群體遺傳多樣性分析[J]. 漁業(yè)研究,41(1):1-10. [Chen J H,Zheng J B,Wang P P,Li T J,Mao Y,Su Y Q,Wang J. 2019. Analysis for genetic diversity of broodstock populations of breeding Litopenaeus vannamei[J]. Journal of Fisheries Research,41(1):1-10.]

    陳靜,何吉祥,樊佳佳,黃龍,吳本麗,宋光同,汪翔,武松. 2018. 草魚MyoD基因SNP和Indel標(biāo)記的篩選及其與生長性狀的關(guān)聯(lián)分析[J]. 江蘇農(nóng)業(yè)學(xué)報(bào),34(3):612-616. [Chen J,He J X,F(xiàn)an J J,Huang L,Wu B L,Song G T,Wang X,Wu S. 2018. Screening of SNP and Indel mar-ker of MyoD gene and its association with growth traits in Ctenopharyngodon idella[J]. Jiangsu Journal of Agricultural Sciences,34(3):612-616.]

    崔蕾,謝從新,李艷和,于永耀,劉肖蓮,張杰,周小云. 2012. 斑點(diǎn)叉尾鮰4個(gè)群體遺傳多樣性的微衛(wèi)星分析[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào),31(6):744-751. [Cui L,Xie C X,Li Y H,Yu Y Y,Liu X L,Zhang J,Zhou X Y. 2012. Analysis of genetic diversity among four different channel catfish populations by using microsatellite markers[J]. Journal of Huazhong Agricultural University,31(6):744-751.]

    戴習(xí)林,劉潔,李晶晶,江宗冰,明磊,蔣飛. 2017. 羅氏沼蝦種群SSR分析中樣本量及標(biāo)記量對遺傳多樣性指標(biāo)的影響[J]. 水產(chǎn)學(xué)報(bào),41(7):1083-1095. [Dai X L,Liu J,Li J J,Jiang Z B,Ming L,Jiang F. 2017. Effects of sample size and loci number on genetic diversity index in Macrobrachium rosenbergii with microsatellite markers[J]. Journal of Fisheries of China,41(7):1083-1095.]

    戴習(xí)林,楊展昆,朱其建. 2018. 凡納濱對蝦室內(nèi)養(yǎng)殖密度和簡易水質(zhì)調(diào)控措施對水質(zhì)及養(yǎng)殖效果的影響[J]. 上海海洋大學(xué)學(xué)報(bào),27(6):894-906. [Dai X L,Yang Z K,Zhu Q J. 2018. Effects of indoor aquaculture density and simple water quality control measures on water quality and culture efficiency of Litopenaeus vannamei[J]. Journal of Shanghai Ocean University,27(6):894-906.]

    董在杰,劉念,傅建軍,朱文彬,王蘭梅,蘇勝彥. 2018. 6個(gè)野生與選育鯉群體的微衛(wèi)星遺傳分析[J]. 南方水產(chǎn)科學(xué),14(4):48-57. [Dong Z J,Liu N,F(xiàn)u J J,Zhu W B,Wang L M,Su S Y. 2018. Genetic analysis for six wild and selection populations of common carp(Cyprinus carpio) using microsatellites[J]. South China Fisheries Science,14(4):48-57.]

    杜育哲,郭世宜,劉麗. 2001. 甲殼動(dòng)物保幼激素的研究[J]. 天津農(nóng)學(xué)院學(xué)報(bào),8(3):34-36. [Du Y Z,Guo S Y,Liu L. 2001. Juvenile hormone in crustacean[J]. Journal of Tianjin Agricultural College,8(3):34-36.]

    杜志強(qiáng). 2012. 先天免疫系統(tǒng)中抗菌肽的研究進(jìn)展[J]. 廣東農(nóng)業(yè)科學(xué),39(10):231-233. [Du Z Q. 2012. Research progress in antimicrobial peptides of innate immune system[J]. Guangdong Agricultural Sciences,39(10):231-233.]

    杜志強(qiáng). 2013. 凡納濱對蝦對蝦素3a基因的克隆與序列分析[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)),28(1):79-82. [Du Z Q. 2013. Gene cloning and sequence analysis of Penaeidin3a in Litopenaeus vannamei[J]. Journal of Yunnan Agricultural University(Natural Science),28(1):79-82.]

    樊廷俊. 2002. 對蝦非特異性免疫與對蝦疾病監(jiān)控的研究進(jìn)展[J]. 海洋科學(xué),26(4):26-31. [Fan T J. 2002. Resear-ching Progresses in shrimp non-specific immunity and di-sease control[J]. Marine Sciences,26(4):26-31.]

    方振朋,孟憲紅,李旭鵬,欒生,曹家旺,陳寶龍,孔杰,閆茂倉,胡利華. 2019. 基于微衛(wèi)星分子標(biāo)記的凡納濱對蝦商業(yè)苗種遺傳多樣性分析[J].漁業(yè)科學(xué)進(jìn)展,41(5):101-109. [Fang Z P,Meng X H,Li X P,Luan S,Cao J W,Chen B L,Kong J,Yan M C,Hu L H. 2019. Genetic diversity analysis of domestic commercial brands seedlings of Litopenaeus vannamei based on microsatellite molecular markers[J]. Progress in Fishery Sciences,41(5):101-109.]

    郝敏. 2016. 凡納濱對與生長繁殖相關(guān)的遺傳標(biāo)記的篩選[D]. 上海:上海海洋大學(xué). [Hao M. 2016. The genetic marker screening for growth and reproduction properties in Litopenaeus vannamei[D]. Shanghai:Shanghai Ocean University.]

    胡世康,王博,秦海鵬,廖栩崢,孫成波,劉永勝,曲朋. 2018. 抗菌肽對凡納濱對蝦抗病性和免疫指標(biāo)的影響[J]. 熱帶生物學(xué)報(bào),9(3):281-286. [Hu S K,Wang B,Qin H P,Liao X Z,Sun C B,Liu Y S,Qu P. 2018. Effect of microbial antimicrobial peptides on disease resistance and immune indices of Litopenaeus vannamei[J]. Journal of Tropical Biology,9(3):281-286.]

    黃小帥,徐煜,胡曉娟,徐武杰,蘇浩昌,文國樑,楊鏗,曹煜成. 2019. 利用微衛(wèi)星標(biāo)記分析7個(gè)凡納濱對蝦引進(jìn)群體子一代的遺傳多樣性[J]. 南方水產(chǎn)科學(xué),15(1):54-62. [Huang X S,Xu Y,Hu X J,Xu W J,Su H C,Wen G L,Yang K,Cao Y C. 2019. Genetic diversity analysis of first filial generation of seven introduced Litopenaeus vannamei populations using microsatellite DNA markers[J]. South China Fisheries Science,15(1):54-62.]

    黃旭雄,羅詞興,郭騰飛,溫文,華雪銘,周洪琪. 2012. 凡納濱對蝦Toll受體和溶菌酶mRNA組織表達(dá)的差異研究[J]. 上海海洋大學(xué)學(xué)報(bào),21(1):28-32. [Huang X X,Luo C X,Guo T F,Wen W,Hua X M,Zhou H Q. 2012. Expre-ssions of Toll receptor and lysozyme mRNA among ti-ssues of white shrimp,Litopenaeus vannamei[J]. Journal of Shanghai Ocean University,21(1):28-32.]

    蔣葛,沈輝,萬夕和,喬毅,黎慧,王李寶,史文軍. 2019. 凡納濱對蝦急性肝胰腺壞死綜合癥病蝦與健康蝦腸道優(yōu)勢菌群比較分析[J]. 江蘇農(nóng)業(yè)學(xué)報(bào),35(1):142-148. [Jiang G,Shen H,Wan X H,Qiao Y,Li H,Wang L B,Shi W J. 2019. Difference analysis of intestinal flora between healthy Litopenaeus vannamei and acute hepatopancreatic necrosis syndrome L. vannamei[J]. Jiangsu Journal of Agricultural Sciences,35(1):142-148.]

    李東宇,孔杰,孟憲紅,欒生,羅坤,盧霞,曹寶祥. 2016. 凡納濱對蝦(Litopenaeus vannamei)微衛(wèi)星多重PCR體系的建立及其在家系親權(quán)鑒定中的應(yīng)用[J]. 漁業(yè)科學(xué)進(jìn)展,37(3):58-67. [Li D Y,Kong J,Meng X H,Luan S,Luo K,Lu X,Cao B X. 2016. Development of multiplex PCR systems of microsatellite markers for Pacific white shrimp(Litopenaeus vannamei) and its application for parentage identification[J]. Progress in Fishery Sciences,37(3):58-67.]

    李強(qiáng)勇,李旻,曾地剛,朱威霖,彭敏,楊春玲,劉青云,趙永貞,陳秀荔,陳曉漢. 2020. 凡納濱對蝦微衛(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 ana-lysis of different cultured families[J]. Journal of Southern Agriculture,51(2):429-436.

    劉洪濤,楊明秋,何玉貴,唐賢明. 2018. 凡納濱對蝦八個(gè)地理群體遺傳多樣性的微衛(wèi)星分析[J]. 海南大學(xué)學(xué)報(bào)(自然科學(xué)版),36(2):146-152. [Liu H T,Yang M Q,He Y G,Tang X M. 2018. Microsatellite analysis of genetic diversity in eight geographical populations of Litopenaeus vannamei[J]. Natural Science Journal of Hainan University,36(2):146-152.]

    劉立鶴,陳立僑,周永奎,江洪波. 2006. 甲殼動(dòng)物消化酶的研究[J]. 飼料工業(yè),27(18):56-62. [Liu L H,Chen L Q,Zhou Y K,Jiang H B. 2006. Studies on digestive enzymes of crustaceans[J]. Feed Industry,27(18):56-62.]

    邱慶慶. 2019. 羅氏沼蝦生長性狀候選基因的克隆及mRNA表達(dá)規(guī)律研究[D]. 南寧:廣西大學(xué). [Qiu Q Q. 2019. Cloning and mRNA expression of candidate genes for growth performance of Macrobrachium rosenbergii[D]. Nanning:Guangxi University.]

    孫成飛,葉星,董浚鍵,田園園,梁健輝. 2015. 羅氏沼蝦6個(gè)養(yǎng)殖群體遺傳多樣性的微衛(wèi)星分析[J]. 南方水產(chǎn)科學(xué),11(2):20-26. [Sun C F,Ye X,Dong J J,Tian Y Y,Liang J H. 2015. Genetic diversity analysis of six cultured populations of Macrobrachium rosenbergii using microsate-llite markers[J]. South China Fisherise Science,11(2):20-26.]

    童馨,龔世圓,喻達(dá)輝,黃桂菊,杜博,李色東. 2009. 凡納濱對蝦(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]

    王軍,王清印,孔杰,孟憲紅,曹家旺,王明珠,馮亞萍,呂丁. 2018. 中國明對蝦人工選育群體與野生群體遺傳多樣性的SSR分析[J]. 漁業(yè)科學(xué)進(jìn)展,39(2):104-111. [Wang J,Wang Q Y,Kong J,Meng X H,Cao J W,Wang M Z,F(xiàn)eng Y P,Lü D. 2018. SSR analysis on genetic diversity in breeding and wild populations of Fenneropenaeus chinensis[J]. Progress in Fishery Sciences,39(2):104-111.]

    王全超. 2017. 凡納濱對蝦生長和抗病性狀的全基因組關(guān)聯(lián)分析與基因組選擇育種研究[D]. 北京:中國科學(xué)院大學(xué). [Wang Q C. 2017. Genome-wide association study and genomic selection of growth and disease resistance traits in Litopenaeus vannamei[D]. Beijing:University of Chinese Academy of Sciences.]

    王沈同,沈玉幫,孟新展,王榮泉,李家樂. 2018. 草魚野生與選育群體遺傳變異微衛(wèi)星分析[J]. 水產(chǎn)學(xué)報(bào),42(8):1273-1284. [Wang S T,Shen Y B,Meng X Z,Wang R Q,Li J L. 2018. Genetic variability in wild and selected populations of Ctenopharyngodon idella using microsatellite markers[J]. Journal of Fisheries of China,42(8):1273-1284.]

    王帥豪. 2012. 淅川烏骨雞群體遺傳多樣性及綠殼蛋性狀相關(guān)基因SNP分析[D]. 鄭州:河南農(nóng)業(yè)大學(xué). [Wang S H. 2012. SNP analysis about green egg shell correlation gene and population genetic diversity of Xichuan silky[D]. Zhengzhou:Henan Agricultural University.]

    武菲,胡文革,王翠華,郭飏,莫超. 2016. 河鱸養(yǎng)殖與野生群體遺傳多樣性比較分析[J]. 水生生物學(xué)報(bào),40(1):181-188. [Wu F,Hu W G,Wang C H,Guo Y,Mo C. 2016. Genetic diversity of the cultivated and natural Perca fluviatilis[J]. Acta Hydrobiologica Sinica,40(1):181-188.]

    夏西超,楊洪,王文鋒,李昕,寧黔冀. 2012. 凡納濱對蝦熱休克蛋白90基因cDNA全長克隆及表達(dá)分析[J]. 水產(chǎn)科學(xué),31(11):673-678. [Xia X C,Yang H,Wang W F,Li X,Ning Q J. 2012. Molecular cloning and expression of heat shock protein 90(HSP90) gene in Pacific white leg shrimp Litopenaeus vannamei[J]. Fisheries Science,31(11):673-678.]

    頡曉勇,蘇天鳳,陳文,張志,顏遠(yuǎn)義,江世貴. 2008. 凡納濱對蝦6個(gè)養(yǎng)殖群體遺傳多樣性的比較分析[J]. 南方水產(chǎn),4(6):42-49. [Xie X Y,Su T F,Chen W,Zhang Z,Yan Y Y,Jiang S G. 2008. Analysis on genetic diversity of six cultured stocks of Litopenaeus vannamei[J]. South China Fisheries Science,4(6):42-49.]

    楊銘,于洋,張曉軍,王全超,劉敬文,李富花,相建海. 2017. 基于轉(zhuǎn)錄組數(shù)據(jù)的凡納濱對蝦微衛(wèi)星標(biāo)記開發(fā)[J]. 海洋科學(xué),41(2):96-102. [Yang M,Yu Y,Zhang X J,Wang Q C,Liu J W,Li F H,Xiang J H. 2017. Development of microsatellite markers from the transcriptome sequences of Pacific white shrimp(Litopenaeus vannamei)[J]. Marine Sciences,41(2):96-102.]

    葉寧,包秀鳳,劉建勇. 2017. 凡納濱對蝦引進(jìn)群體與養(yǎng)殖群體遺傳多樣性的AFLP分析[J]. 水產(chǎn)學(xué)報(bào),41(3):339-346. [Ye N,Bao X F,Liu J Y. 2017. AFLP analysis of genetic diversity of imported Litopenaeus vannamei popu-lations and cultured stocks[J]. Journal of Fisheries of China,41(3):339-346.]

    于磊,劉炳艦,劉進(jìn)賢. 2018. 日本鰻鱺離散SNP標(biāo)記篩選及群體遺傳多樣性分析[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),48(4):37-47. [Yu L,Liu B J,Liu J X. 2018. Outlier SNPs development and preliminary genetic diversity ana-lysis of Japanese eel(Anguilla japonica)[J]. Periodical of Ocean University of China,48(4):37-47.]

    于洋. 2014. 凡納濱對蝦分子標(biāo)記的開發(fā)及其在遺傳育種中的應(yīng)用[D]. 北京:中國科學(xué)院大學(xué). [Yu Y. 2014. Deve-lopment of molecular markers and their applications in selective breeding of the Pacific white shrimp,Litopenaeus vannamei[D]. Beijing:University of Chinese Academy of Sciences.]

    張靈俠,沈琪,胡超群. 2008. 凡納濱對蝦的遺傳育種研究現(xiàn)狀[J]. 海洋科學(xué),32(2):91-95. [Zhang L X,Shen Q,Hu C Q. 2008. A review of genetics and breeding of Litopenaeus vannamei[J]. Marine Sciences,32(2):91-95.]

    張曉萌,馬普,王洪迪,王秀利. 2013. SNPs在水產(chǎn)動(dòng)物中的研究進(jìn)展[J]. 生物技術(shù)通報(bào),(8):7-11. [Zhang X M,Ma P,Wang H D,Wang X L. 2013. Progresses of SNPs studies in aquaculture animals[J]. Biotechnology Bulletin,(8):7-11.]

    張艷艷,劉小林,黃海洪,錢昭英,王憲宗,相建海,李義軍. 2012. 凡納濱對蝦Crustin-like基因的克隆及在副溶血弧菌感染條件下的表達(dá)分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),40(12):41-48. [Zhang Y Y,Liu X L,Huang H H,Qian Z Y,Wang X Z,Xiang J H,Li Y J. 2012. Molecular cloning of Crustin-like gene in the white shrimp(Litopenaeus vannamei) and its mRNA expression with Vibrio parahaemolyticus challenge[J]. Journal of Nor-thwest A & F University(Natural Science Edition),40(12):41-48.]

    Agarwal D,Aich N,Pavan-Kumar A,Kumar S,Sabnis S,Joshi C G,Koringa P,Pandya D,Patel D,Karnik T,Bhingarde R,Gireesh-Babu P,Chaudhari A. 2016. SNP mi-ning in transcripts and concomitant estimation of genetic variation in Macrobrachium rosenbergii stocks[J]. Conservation Genetics Resources,8(2):159-168.

    Rezaee S,F(xiàn)arahmand H,Nematollahi M A. 2016. Genetic diversity status of Pacific white shrimp(Litopenaeus vannamei) using SSR markers in Iran[J]. Aquaculture International,24(2):479-489.

    Thongda W,Zhao H G,Zhang D D,Jescovitch L N,Liu M,Guo X M,Schrandt M,Powers S P,Peatman E. 2018. Development of SNP panels as a new tool to assess the genetic diversity,population structure,and parentage ana-lysis of the eastern oyster(Crassostrea virginica)[J]. Marine Biotechnology(New York),20(3):385-395.

    Yu Y,Wang Q C,Zhang Q,Luo Z,Wang Y,Zhang X J,Huang H,Xiang J H,Li F H. 2019. Genome scan for genomic regions and genes associated with growth trait in Pacific white shrimp Litopenaeus vannamei[J]. Marine Biotechnology(New York),21(3):374-383.

    (責(zé)任編輯 蘭宗寶)

    收稿日期:2020-02-15

    基金項(xiàng)目:上海市蝦類產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(滬農(nóng)科攻字〔2014〕第7-1-11號)

    作者簡介:*為通訊作者,戴習(xí)林(1969-),教授,主要從事水產(chǎn)養(yǎng)殖與海洋生物學(xué)研究工作,E-mail:xldai@shou.edu.cn。申淑慧(1995-),研究方向?yàn)榱_氏沼蝦遺傳育種,E-mail:18994535472@163.com

    猜你喜歡
    遺傳多樣性
    從葉綠體DNA角度分析云南省砂梨地方品種遺傳多樣性
    寧夏外引水稻種質(zhì)資源表型性狀遺傳多樣性分析
    寧夏外引水稻種質(zhì)資源表型性狀遺傳多樣性分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    金魚起源及遺傳多樣性研究進(jìn)展
    金銀花SSR指紋圖譜的構(gòu)建及遺傳多樣性分析
    山西大豆自然群體遺傳多樣性的研究
    2個(gè)鯉魚群體遺傳多樣性的RAPD分析
    日韩高清综合在线| 成年版毛片免费区| 国产精品综合久久久久久久免费| 美女黄网站色视频| 青草久久国产| 黄色丝袜av网址大全| 赤兔流量卡办理| 色哟哟哟哟哟哟| 免费av毛片视频| 欧美成人性av电影在线观看| 亚洲av.av天堂| 深爱激情五月婷婷| 两人在一起打扑克的视频| 啦啦啦韩国在线观看视频| 欧美日韩综合久久久久久 | 欧美高清成人免费视频www| 身体一侧抽搐| 亚洲国产欧美人成| 最近在线观看免费完整版| 一本精品99久久精品77| 在线天堂最新版资源| 国产伦精品一区二区三区视频9| 91在线精品国自产拍蜜月| 免费在线观看成人毛片| 亚洲av中文字字幕乱码综合| 日本五十路高清| 国产三级黄色录像| 性色avwww在线观看| 婷婷六月久久综合丁香| 午夜福利18| 一个人免费在线观看的高清视频| 国产美女午夜福利| 免费观看精品视频网站| 欧美一区二区精品小视频在线| 亚洲国产欧洲综合997久久,| 真人一进一出gif抽搐免费| 最近最新中文字幕大全电影3| 波野结衣二区三区在线| 免费av毛片视频| 久久午夜福利片| 中文字幕免费在线视频6| 99国产精品一区二区三区| 2021天堂中文幕一二区在线观| 久久久久久大精品| 国产视频内射| 少妇高潮的动态图| 他把我摸到了高潮在线观看| 真人一进一出gif抽搐免费| 精品一区二区三区av网在线观看| 动漫黄色视频在线观看| 欧美精品国产亚洲| 久久久久性生活片| 日韩国内少妇激情av| 日韩成人在线观看一区二区三区| 超碰av人人做人人爽久久| 亚洲成a人片在线一区二区| 在线国产一区二区在线| 亚洲国产精品999在线| 制服丝袜大香蕉在线| 久久99热6这里只有精品| 9191精品国产免费久久| 婷婷精品国产亚洲av在线| 丰满乱子伦码专区| 久久精品国产99精品国产亚洲性色| 日日摸夜夜添夜夜添av毛片 | 精品午夜福利在线看| 男人的好看免费观看在线视频| 高清在线国产一区| 变态另类成人亚洲欧美熟女| 国产高清视频在线观看网站| 热99re8久久精品国产| 亚洲三级黄色毛片| 亚洲乱码一区二区免费版| 国产精品女同一区二区软件 | 99国产精品一区二区三区| 婷婷色综合大香蕉| 欧美成人一区二区免费高清观看| 看免费av毛片| 2021天堂中文幕一二区在线观| 制服丝袜大香蕉在线| 极品教师在线视频| 熟女电影av网| 草草在线视频免费看| 色精品久久人妻99蜜桃| 脱女人内裤的视频| 精品人妻1区二区| 俄罗斯特黄特色一大片| 老司机午夜福利在线观看视频| 国产精品99久久久久久久久| 亚洲精品成人久久久久久| 亚洲av熟女| 精品99又大又爽又粗少妇毛片 | 他把我摸到了高潮在线观看| 国产人妻一区二区三区在| 国产精品野战在线观看| 国产av在哪里看| 我要搜黄色片| 中文字幕人妻熟人妻熟丝袜美| 中文亚洲av片在线观看爽| 亚洲精品色激情综合| 真人一进一出gif抽搐免费| 亚洲av中文字字幕乱码综合| 少妇被粗大猛烈的视频| 亚洲五月婷婷丁香| 亚洲中文日韩欧美视频| 3wmmmm亚洲av在线观看| 亚洲最大成人手机在线| 国产精品亚洲美女久久久| 丁香欧美五月| 美女cb高潮喷水在线观看| 看免费av毛片| 亚洲国产精品sss在线观看| 亚洲欧美日韩卡通动漫| 国产高清三级在线| 老熟妇仑乱视频hdxx| or卡值多少钱| 亚洲不卡免费看| 小说图片视频综合网站| a级毛片免费高清观看在线播放| 在线看三级毛片| 亚洲欧美日韩东京热| 国产主播在线观看一区二区| 99在线视频只有这里精品首页| 国产激情偷乱视频一区二区| 午夜福利欧美成人| 国产亚洲精品久久久com| 美女被艹到高潮喷水动态| 国产精品久久久久久亚洲av鲁大| 夜夜爽天天搞| 可以在线观看毛片的网站| 久久精品人妻少妇| 女人十人毛片免费观看3o分钟| 女人被狂操c到高潮| 亚洲va日本ⅴa欧美va伊人久久| 国产主播在线观看一区二区| 亚洲成人久久爱视频| 国产亚洲精品综合一区在线观看| 日本 av在线| 十八禁人妻一区二区| 中文字幕久久专区| 91av网一区二区| 久久天躁狠狠躁夜夜2o2o| av在线老鸭窝| 久久人人精品亚洲av| 观看美女的网站| 免费黄网站久久成人精品 | 一级作爱视频免费观看| 久久久久久久久久成人| 香蕉av资源在线| 精品欧美国产一区二区三| 美女被艹到高潮喷水动态| 亚洲真实伦在线观看| 亚洲七黄色美女视频| 99热这里只有精品一区| 国产高清视频在线播放一区| 亚洲av熟女| 日韩欧美国产在线观看| 国产三级中文精品| 性插视频无遮挡在线免费观看| 乱码一卡2卡4卡精品| 男女之事视频高清在线观看| av在线老鸭窝| 色在线成人网| 在现免费观看毛片| 日本熟妇午夜| 免费av毛片视频| 真人一进一出gif抽搐免费| netflix在线观看网站| 精品久久久久久,| 色精品久久人妻99蜜桃| 91在线观看av| 亚洲精品成人久久久久久| 一区二区三区免费毛片| 中文资源天堂在线| 亚洲,欧美精品.| 51国产日韩欧美| 国产免费男女视频| 波多野结衣高清作品| 老司机午夜十八禁免费视频| 亚洲人成网站在线播放欧美日韩| 制服丝袜大香蕉在线| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区 | 午夜福利成人在线免费观看| 国内少妇人妻偷人精品xxx网站| 一个人免费在线观看电影| www.熟女人妻精品国产| av在线观看视频网站免费| 啪啪无遮挡十八禁网站| 99久久99久久久精品蜜桃| 亚洲成人中文字幕在线播放| 俺也久久电影网| 国产大屁股一区二区在线视频| 成年版毛片免费区| 久久精品国产自在天天线| 观看美女的网站| 日韩亚洲欧美综合| 又爽又黄a免费视频| 99久久99久久久精品蜜桃| 日韩中字成人| 黄色日韩在线| 一进一出抽搐gif免费好疼| 五月玫瑰六月丁香| www.www免费av| 日本撒尿小便嘘嘘汇集6| 国产伦在线观看视频一区| 午夜免费成人在线视频| 亚洲精品粉嫩美女一区| 嫁个100分男人电影在线观看| 特大巨黑吊av在线直播| 久久久久免费精品人妻一区二区| 人妻夜夜爽99麻豆av| 99在线人妻在线中文字幕| 日韩亚洲欧美综合| 国产高清三级在线| 免费av不卡在线播放| 身体一侧抽搐| 18禁黄网站禁片免费观看直播| 国产精品三级大全| 国产亚洲精品久久久com| 在线播放国产精品三级| 成人国产综合亚洲| 日本与韩国留学比较| 国产一区二区三区在线臀色熟女| 国产精品嫩草影院av在线观看 | 久久久久久久午夜电影| 99热6这里只有精品| 9191精品国产免费久久| 97人妻精品一区二区三区麻豆| 51午夜福利影视在线观看| 99久久精品国产亚洲精品| 国产午夜精品论理片| 国产精品国产高清国产av| 最近最新中文字幕大全电影3| 在线十欧美十亚洲十日本专区| 国产野战对白在线观看| 免费人成视频x8x8入口观看| 九色国产91popny在线| 久久热精品热| 婷婷六月久久综合丁香| 欧美色视频一区免费| 啦啦啦观看免费观看视频高清| 国产一区二区三区视频了| 露出奶头的视频| 蜜桃亚洲精品一区二区三区| 人妻制服诱惑在线中文字幕| 色在线成人网| 嫩草影院新地址| 成人亚洲精品av一区二区| h日本视频在线播放| 神马国产精品三级电影在线观看| 麻豆成人av在线观看| 亚洲av成人不卡在线观看播放网| 国产黄a三级三级三级人| 亚洲av免费高清在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产av麻豆久久久久久久| 人妻夜夜爽99麻豆av| 欧美一区二区国产精品久久精品| 毛片女人毛片| 中文字幕高清在线视频| 在线观看一区二区三区| 一区二区三区免费毛片| 久久久久精品国产欧美久久久| 亚洲成人免费电影在线观看| 日本 av在线| 国产精品一区二区免费欧美| 精品久久久久久久久久免费视频| 午夜影院日韩av| 99久久无色码亚洲精品果冻| 日韩欧美国产在线观看| 亚洲 国产 在线| 最近视频中文字幕2019在线8| 成人亚洲精品av一区二区| 日韩免费av在线播放| 性色av乱码一区二区三区2| 色av中文字幕| 黄色一级大片看看| 午夜福利18| 毛片一级片免费看久久久久 | 国产私拍福利视频在线观看| 搡老妇女老女人老熟妇| 少妇人妻精品综合一区二区 | 国产精品99久久久久久久久| 草草在线视频免费看| www.www免费av| 亚洲第一电影网av| 国产精品99久久久久久久久| 亚洲精品成人久久久久久| 精品免费久久久久久久清纯| 日本一二三区视频观看| 午夜福利成人在线免费观看| 日韩欧美在线乱码| 亚洲专区中文字幕在线| 国产精品三级大全| 美女高潮的动态| 内射极品少妇av片p| 精华霜和精华液先用哪个| 欧美在线一区亚洲| 欧美日韩亚洲国产一区二区在线观看| 99视频精品全部免费 在线| 日韩欧美三级三区| 日韩有码中文字幕| av专区在线播放| 亚洲国产精品999在线| 日本免费a在线| 黄色丝袜av网址大全| 欧美+亚洲+日韩+国产| 桃红色精品国产亚洲av| 国产精品亚洲一级av第二区| 中国美女看黄片| 亚洲最大成人中文| 最近中文字幕高清免费大全6 | .国产精品久久| 免费观看人在逋| 亚洲中文日韩欧美视频| 国产成人福利小说| 美女被艹到高潮喷水动态| 日本精品一区二区三区蜜桃| 十八禁国产超污无遮挡网站| 欧美黄色片欧美黄色片| 国产精品永久免费网站| 91午夜精品亚洲一区二区三区 | 俺也久久电影网| 欧美国产日韩亚洲一区| 一区二区三区激情视频| 长腿黑丝高跟| 精品人妻熟女av久视频| 小蜜桃在线观看免费完整版高清| 国产精品影院久久| 一个人免费在线观看的高清视频| 精品午夜福利视频在线观看一区| 一本综合久久免费| 色吧在线观看| 久久久久久大精品| 亚洲第一电影网av| 自拍偷自拍亚洲精品老妇| 亚洲成av人片免费观看| 91狼人影院| 亚州av有码| 99久久精品国产亚洲精品| 熟女电影av网| 91av网一区二区| 国产探花极品一区二区| 黄片小视频在线播放| 欧美丝袜亚洲另类 | 少妇裸体淫交视频免费看高清| 日本三级黄在线观看| 国产高潮美女av| 日本a在线网址| 欧美潮喷喷水| eeuss影院久久| 午夜精品久久久久久毛片777| 免费无遮挡裸体视频| 免费av毛片视频| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 国产黄色小视频在线观看| 国产成人啪精品午夜网站| 一区二区三区高清视频在线| 99国产综合亚洲精品| 中文在线观看免费www的网站| 日日摸夜夜添夜夜添小说| 亚洲av五月六月丁香网| 久久亚洲精品不卡| or卡值多少钱| 别揉我奶头 嗯啊视频| 一个人看的www免费观看视频| 成人特级av手机在线观看| 蜜桃久久精品国产亚洲av| 国产精品影院久久| 禁无遮挡网站| 好看av亚洲va欧美ⅴa在| 深夜a级毛片| 欧美xxxx黑人xx丫x性爽| 午夜视频国产福利| av国产免费在线观看| 亚洲av免费高清在线观看| 熟女电影av网| 欧美在线一区亚洲| 可以在线观看毛片的网站| 免费电影在线观看免费观看| 亚洲av电影不卡..在线观看| 神马国产精品三级电影在线观看| 麻豆成人午夜福利视频| av天堂中文字幕网| 精品久久国产蜜桃| 国产色爽女视频免费观看| 搡老岳熟女国产| 老司机午夜十八禁免费视频| 99久久久亚洲精品蜜臀av| 精品免费久久久久久久清纯| 国产av麻豆久久久久久久| 国产精品一及| 色综合欧美亚洲国产小说| 午夜福利高清视频| 中亚洲国语对白在线视频| 免费人成在线观看视频色| 久久久成人免费电影| 欧美色欧美亚洲另类二区| 欧美三级亚洲精品| av视频在线观看入口| 好看av亚洲va欧美ⅴa在| 亚洲人成网站高清观看| 熟妇人妻久久中文字幕3abv| 日本五十路高清| 欧美乱妇无乱码| 亚洲精品一卡2卡三卡4卡5卡| 有码 亚洲区| 性欧美人与动物交配| 两个人的视频大全免费| 91久久精品国产一区二区成人| 亚洲国产日韩欧美精品在线观看| 亚洲熟妇中文字幕五十中出| 最近最新中文字幕大全电影3| 精品国产三级普通话版| 婷婷丁香在线五月| 2021天堂中文幕一二区在线观| 国产亚洲av嫩草精品影院| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久午夜电影| 黄色丝袜av网址大全| 久久精品国产自在天天线| 色视频www国产| 中文字幕av成人在线电影| 国产成人影院久久av| 天天一区二区日本电影三级| 亚洲内射少妇av| 久久人妻av系列| 国产高清激情床上av| 午夜激情福利司机影院| 精品国产亚洲在线| 欧美激情在线99| 在线观看午夜福利视频| netflix在线观看网站| 日韩欧美三级三区| 久久伊人香网站| av欧美777| 极品教师在线视频| 亚洲乱码一区二区免费版| 国产视频内射| 午夜免费成人在线视频| 热99re8久久精品国产| 麻豆久久精品国产亚洲av| 十八禁人妻一区二区| 欧美区成人在线视频| 中国美女看黄片| 伦理电影大哥的女人| 一个人免费在线观看电影| 国产成人a区在线观看| 99久久九九国产精品国产免费| 国产精品久久久久久亚洲av鲁大| 亚洲自偷自拍三级| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲国产精品sss在线观看| 亚洲成人久久爱视频| 好男人电影高清在线观看| 国产美女午夜福利| 精品人妻偷拍中文字幕| 男人的好看免费观看在线视频| 国产精品1区2区在线观看.| 欧美激情在线99| 搡老熟女国产l中国老女人| 90打野战视频偷拍视频| 国产av麻豆久久久久久久| 久久久久久国产a免费观看| 日日摸夜夜添夜夜添小说| 国产成年人精品一区二区| 亚洲午夜理论影院| 亚洲av成人精品一区久久| 免费av毛片视频| 老女人水多毛片| 色吧在线观看| 尤物成人国产欧美一区二区三区| 国产一区二区在线观看日韩| 麻豆av噜噜一区二区三区| 日韩免费av在线播放| 18禁黄网站禁片午夜丰满| 国产三级中文精品| 99热6这里只有精品| 国产精品综合久久久久久久免费| 久久性视频一级片| 网址你懂的国产日韩在线| av视频在线观看入口| 搞女人的毛片| 嫩草影院精品99| 亚洲中文日韩欧美视频| 色综合亚洲欧美另类图片| 我要看日韩黄色一级片| 人妻丰满熟妇av一区二区三区| 国产欧美日韩精品亚洲av| 国产在线精品亚洲第一网站| 免费av毛片视频| 99热只有精品国产| 国产久久久一区二区三区| 国产精品久久久久久精品电影| 日本三级黄在线观看| 午夜影院日韩av| 午夜亚洲福利在线播放| 日日摸夜夜添夜夜添小说| 欧美性猛交╳xxx乱大交人| 成人精品一区二区免费| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添av毛片 | 久久精品影院6| 精品国内亚洲2022精品成人| 亚洲精品色激情综合| 日本一本二区三区精品| 悠悠久久av| 国产一区二区在线av高清观看| 国产精品电影一区二区三区| 亚洲av日韩精品久久久久久密| 韩国av一区二区三区四区| 午夜福利免费观看在线| 亚洲人成网站高清观看| 少妇的逼好多水| 国产亚洲精品综合一区在线观看| 露出奶头的视频| 在线播放无遮挡| 51午夜福利影视在线观看| 韩国av一区二区三区四区| 色综合站精品国产| 中文资源天堂在线| 狠狠狠狠99中文字幕| 黄色配什么色好看| 国产一区二区激情短视频| 亚洲无线观看免费| 免费大片18禁| 国产免费男女视频| 亚洲aⅴ乱码一区二区在线播放| av福利片在线观看| 国产亚洲精品久久久com| av在线蜜桃| 国产麻豆成人av免费视频| 99国产极品粉嫩在线观看| 性欧美人与动物交配| 一a级毛片在线观看| 国产精品永久免费网站| 99热只有精品国产| 小蜜桃在线观看免费完整版高清| 18禁在线播放成人免费| 好看av亚洲va欧美ⅴa在| 欧美激情在线99| 国产精品1区2区在线观看.| 日本一二三区视频观看| 久久久久久九九精品二区国产| 老司机午夜十八禁免费视频| 婷婷精品国产亚洲av在线| 69av精品久久久久久| 熟女人妻精品中文字幕| av在线蜜桃| 好男人在线观看高清免费视频| 国产免费一级a男人的天堂| 国产美女午夜福利| 少妇人妻一区二区三区视频| 欧美国产日韩亚洲一区| 丰满的人妻完整版| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 99久久成人亚洲精品观看| 国产高清视频在线观看网站| 亚洲欧美日韩高清在线视频| 久久久久久久午夜电影| 亚洲va日本ⅴa欧美va伊人久久| 久久天躁狠狠躁夜夜2o2o| 动漫黄色视频在线观看| 成人性生交大片免费视频hd| 亚洲精品在线观看二区| 性欧美人与动物交配| 亚洲片人在线观看| 天堂网av新在线| 露出奶头的视频| 在线十欧美十亚洲十日本专区| 亚洲五月天丁香| 国产一区二区三区视频了| 国产成人aa在线观看| 精品99又大又爽又粗少妇毛片 | 大型黄色视频在线免费观看| 欧美日韩福利视频一区二区| 一进一出抽搐gif免费好疼| 成人鲁丝片一二三区免费| 2021天堂中文幕一二区在线观| 99久国产av精品| 国内揄拍国产精品人妻在线| bbb黄色大片| 校园春色视频在线观看| av专区在线播放| 国产精品三级大全| 精品久久久久久久人妻蜜臀av| 中文字幕熟女人妻在线| 亚洲,欧美,日韩| 岛国在线免费视频观看| 又紧又爽又黄一区二区| 国产精品乱码一区二三区的特点| 最后的刺客免费高清国语| 日韩免费av在线播放| 特大巨黑吊av在线直播| 亚州av有码| 国产伦在线观看视频一区| 男女视频在线观看网站免费| 国产午夜福利久久久久久| 非洲黑人性xxxx精品又粗又长| 夜夜看夜夜爽夜夜摸| 18禁裸乳无遮挡免费网站照片| 在线观看av片永久免费下载| 国产v大片淫在线免费观看| 国产精品免费一区二区三区在线| 午夜日韩欧美国产| 久久久久久久精品吃奶| 三级国产精品欧美在线观看| 少妇人妻精品综合一区二区 | av在线蜜桃| 老司机午夜福利在线观看视频| 国产精品免费一区二区三区在线|