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

    長江刀鱭選育群體轉(zhuǎn)錄組EST-SSR的分布特征分析*

    2019-09-27 01:42:12于愛清施永海徐嘉波陸根海張海明謝永德劉永士
    漁業(yè)科學(xué)進(jìn)展 2019年5期
    關(guān)鍵詞:基序微衛(wèi)星核苷酸

    于愛清 施永海 徐嘉波 陸根海 張海明 謝永德 劉永士

    長江刀鱭選育群體轉(zhuǎn)錄組EST-SSR的分布特征分析*

    于愛清 施永海①徐嘉波 陸根海 張海明 謝永德 劉永士

    (上海市水產(chǎn)研究所 上海市水產(chǎn)技術(shù)推廣站 上海 200433)

    本研究利用MISA軟件挖掘長江刀鱭()肌肉和肝臟轉(zhuǎn)錄組中的微衛(wèi)星標(biāo)記,為刀鱭選育群體的種質(zhì)資源評(píng)估和分子標(biāo)記輔助育種奠定基礎(chǔ)。結(jié)果顯示,從71869條Unigenes中共獲得33896條重復(fù)單元長度為1~6堿基的微衛(wèi)星序列;刀鱭轉(zhuǎn)錄組中不同類型微衛(wèi)星的重復(fù)基序具有不同的分布特征,其中,單核苷酸重復(fù)、二核苷酸重復(fù)和三核苷酸重復(fù)為主要的微衛(wèi)星重復(fù)類型,分別占總微衛(wèi)星數(shù)目的34.94%、49.47%和13.34%;不同微衛(wèi)星重復(fù)類型的優(yōu)勢(shì)重復(fù)基序亦有所不同,其中,A/T為單核苷酸重復(fù)基序的優(yōu)勢(shì)重復(fù)基序占86.25%,AC/GT為二核苷酸重復(fù)基序的為優(yōu)勢(shì)重復(fù)基序占75.25%,AGG/CCT為三核苷酸重復(fù)基序的優(yōu)勢(shì)重復(fù)基序占28.57%;不同微衛(wèi)星重復(fù)基序核苷酸的數(shù)量和重復(fù)次數(shù)亦有所不同,重復(fù)次數(shù)伴隨著重復(fù)單元中核苷酸數(shù)量的增加而呈現(xiàn)降低的趨勢(shì);從100對(duì)四核苷酸重復(fù)的SSR引物中篩選獲得了16對(duì)多態(tài)性微衛(wèi)星標(biāo)記,并以此為基礎(chǔ),對(duì)長江刀鱭選育群體(F3)的遺傳學(xué)特征進(jìn)行了初步評(píng)估,結(jié)果顯示,長江刀鱭選育群體F3的平均有效等位基因數(shù)(N)、平均觀測(cè)雜合度(H)、平均期望雜合度(H)和Shannon多樣性指數(shù)分別為1.7580、0.3414、0.3977和0.6278。以上結(jié)果表明,基于刀鱭轉(zhuǎn)錄組數(shù)據(jù)批量開發(fā)微衛(wèi)星是切實(shí)可行的,所開發(fā)的多態(tài)性微衛(wèi)星標(biāo)記能夠應(yīng)用于長江刀鱭選育群體的遺傳背景評(píng)估和進(jìn)一步的遺傳育種研究。

    刀鱭;轉(zhuǎn)錄組;微衛(wèi)星;分子標(biāo)記

    刀鱭()隸屬于鯡形目(Clupeiformes)、鳀科(Engraulidae)、鱭屬(),又名長頜鱭,俗稱刀魚、毛刀魚,主要分布于我國黃渤海、東海海域及通海的江河,其中,又以長江下游的產(chǎn)量最高,因其魚體豐腴肥滿、肉質(zhì)細(xì)嫩鮮美而享有“長江三鮮”之首的美譽(yù)(施永海等, 2014、2015; 鄧平平等, 2014)。歷史上,我國刀鱭的種質(zhì)資源極為豐富。然而,近年來,由于酷捕濫漁、水域環(huán)境污染、水文條件改變、生態(tài)環(huán)境惡化及海岸工程建設(shè)等諸多因素,長江刀鱭自然種質(zhì)資源急劇衰退,產(chǎn)量波動(dòng)很大,且呈現(xiàn)出逐年降低且個(gè)體日趨小型化的趨勢(shì),甚至不能形成優(yōu)勢(shì)種群(施永海等, 2014、2015; 魏廣蓮等, 2013; 袁傳宓, 1988),亟需對(duì)刀鱭的種質(zhì)資源進(jìn)行保護(hù)和恢復(fù)。而刀鱭基因組信息匱乏、遺傳背景不清晰、繁殖周期較長等問題亦在一定程度上制約了其種質(zhì)資源的開發(fā)和恢復(fù)。刀鱭種質(zhì)資源的有效保護(hù)和恢復(fù),除了政府相關(guān)部門采取相關(guān)行政干預(yù)(建立禁漁期、加強(qiáng)捕撈監(jiān)管、適時(shí)增殖放流等)之外,還需相關(guān)部門適時(shí)地監(jiān)測(cè)洄游型刀鱭的遺傳多樣性,從整體上把握刀鱭野生群體的遺傳多樣性水平和遺傳背景,同時(shí)加快刀鱭良種的選育進(jìn)程,特別是分子標(biāo)記輔助良種選育的研究。

    分子標(biāo)記由于能夠?qū)?dòng)物整個(gè)基因組的遺傳變異進(jìn)行探究,且相對(duì)穩(wěn)定,而被應(yīng)用于水產(chǎn)動(dòng)物的遺傳育種研究,其中,微衛(wèi)星標(biāo)記(Microsatellite)的應(yīng)用最為廣泛。微衛(wèi)星標(biāo)記,又稱為簡(jiǎn)單重復(fù)序列(Simple sequence repeat, SSR),是指由1~6個(gè)核苷酸組成的簡(jiǎn)單串聯(lián)重復(fù)序列,這些序列的重復(fù)單位和重復(fù)次數(shù)的差異使其具有高度變異性和多態(tài)性,被廣泛應(yīng)用于水產(chǎn)動(dòng)物種質(zhì)資源評(píng)估和鑒定、遺傳圖譜的構(gòu)建及分子標(biāo)記輔助育種等方面(楊銘等, 2017; 管奧等, 2018; 白翠翠等, 2016)。近年來,學(xué)者利用線粒體基因D-loop(張燕萍等, 2017; Cheng, 2011)、線粒體Cytb (魏廣蓮等, 2012; Ma, 2010)、RAPD (馬春艷等, 2004)、AFLP(葛家春等, 2008)等方法研究了刀鱭不同群體的遺傳結(jié)構(gòu)和遺傳差異,而對(duì)刀鱭微衛(wèi)星標(biāo)記開發(fā)的報(bào)道相對(duì)較少。Rong等(2013)利用FIASO法構(gòu)建了微衛(wèi)星富集文庫方法,從89對(duì)刀鱭微衛(wèi)星引物中獲得了20對(duì)多態(tài)性微衛(wèi)星標(biāo)記。馬春燕等(2011)利用FIASO法構(gòu)建微衛(wèi)星富集文庫的方法,從70對(duì)刀鱭微衛(wèi)星引物中獲得了12對(duì)多態(tài)性微衛(wèi)星標(biāo)記。陳芳等(2012)利用磁珠富集法,從198對(duì)刀鱭微衛(wèi)星引物中獲得了34對(duì)多態(tài)性微衛(wèi)星標(biāo)記。鄧平平等(2014)利用磁珠富集法,從59對(duì)刀鱭微衛(wèi)星引物中獲得了9對(duì)多態(tài)性微衛(wèi)星標(biāo)記。利用磁珠富集法從水產(chǎn)動(dòng)物基因組中開發(fā)的微衛(wèi)星標(biāo)記,多態(tài)性相對(duì)較高,但相對(duì)費(fèi)時(shí)費(fèi)力,效率相對(duì)較低,且標(biāo)記大規(guī)模的開發(fā)成本相對(duì)較高,遠(yuǎn)遠(yuǎn)不能滿足刀鱭高密度遺傳圖譜的構(gòu)建、QTL精細(xì)定位和分子標(biāo)記輔助育種等對(duì)分子標(biāo)記需求量較大的研究,亟需大規(guī)模開發(fā)分子標(biāo)記。而高通量轉(zhuǎn)錄組測(cè)序技術(shù)所開發(fā)的微衛(wèi)星,又被稱為表達(dá)序列標(biāo)簽SSR (Expressed sequence tag SSR, EST-SSR),其多態(tài)性雖然相對(duì)低于基因組所開發(fā)的微衛(wèi)星(G-SSR),但由于數(shù)量多、通量大、周期短、成本低,且直接與某些功能基因相關(guān)聯(lián)(李東明等, 2017),而被廣泛應(yīng)用于水產(chǎn)動(dòng)物分子標(biāo)記大規(guī)模的開發(fā)研究。

    本研究基于刀鱭肌肉和肝臟的高通量轉(zhuǎn)錄組數(shù)據(jù),對(duì)長江刀鱭轉(zhuǎn)錄組中的微衛(wèi)星位點(diǎn)進(jìn)行篩選,分析刀鱭EST-SSR的序列特征、分布和組成等信息,并利用所篩選的多態(tài)性微衛(wèi)星標(biāo)記對(duì)長江刀鱭選育群體(F3)的遺傳學(xué)特征進(jìn)行初步評(píng)估,以期能夠?yàn)榈恩q遺傳多樣性的分析、種質(zhì)資源的保護(hù)和利用、高密度遺傳連鎖圖譜的構(gòu)建及分子標(biāo)記輔助良種選育奠定基礎(chǔ)。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料與數(shù)據(jù)來源

    實(shí)驗(yàn)所用長江刀鱭來自于上海市水產(chǎn)研究所核心選育群體(F3)的2齡刀鱭,提取刀鱭肌肉和肝臟的RNA后,構(gòu)建高通量轉(zhuǎn)錄組測(cè)序文庫,利用Illumina HiSeq4000進(jìn)行高通量轉(zhuǎn)錄組測(cè)序。刀鱭肌肉和肝臟的轉(zhuǎn)錄組測(cè)序總共獲得28.20 Gb的數(shù)據(jù)量,原始測(cè)序數(shù)據(jù)經(jīng)Trinity軟件(http://trinityrnaseq.sourceforge. net/)組裝并去冗余后,獲得平均長度為932.49 bp的Unigenes 71869條,N50達(dá)到1364 bp,質(zhì)量不低于20的堿基比例(Q20)為97.31%,表明測(cè)序結(jié)果良好(未發(fā)表)。

    1.2 EST-SSR的篩查

    以刀鱭轉(zhuǎn)錄組中的Unigenes作為參考序列,利用MISA(MIcroSAtellite identification tool)軟件(http: //pgrc.ipk-gatersleben.de/misa/misa.html)搜索刀鱭肝臟和肌肉轉(zhuǎn)錄組中的SSR位點(diǎn),其中,單核苷酸、二核苷酸、三核苷酸、四核苷酸、五核苷酸和六核苷酸的最少重復(fù)次數(shù)分別設(shè)置10、6、5、5、5和5。復(fù)合SSR兩個(gè)位點(diǎn)間最大間隔堿基數(shù)設(shè)置為100。

    1.3 EST-SSR多態(tài)性引物的設(shè)計(jì)、篩選和初步驗(yàn)證

    利用Primer Premier 3.0軟件對(duì)篩選出的微衛(wèi)星序列批量設(shè)計(jì)引物,主要設(shè)置參數(shù)為:引物長度和PCR產(chǎn)物分別為18~25 bp和100~500 bp,55℃≤退火溫度(T)≤65℃,40%≤GC含量≤60%,上下游引物的T≤2℃??紤]到四堿基重復(fù)的微衛(wèi)星標(biāo)記相對(duì)二堿基重復(fù)和三堿基重復(fù)具有較長的演化歷史和更加穩(wěn)定的突變積累,本研究隨機(jī)選取100對(duì)基于四核苷酸重復(fù)微衛(wèi)星序列所設(shè)計(jì)的PCR擴(kuò)增引物交由生工生物工程(上海)股份有限公司合成,并以刀鱭核心選育群體(F3)為實(shí)驗(yàn)材料篩選具有多態(tài)性的微衛(wèi)星引物。PCR擴(kuò)增產(chǎn)物經(jīng)10%非變性聚丙烯酰胺凝膠電泳檢測(cè)后,利用銀染法顯色定影,條帶顯色清晰后用數(shù)碼相機(jī)拍照保存,進(jìn)行微衛(wèi)星引物后續(xù)的多態(tài)性分析。利用GenAlEx 6.0軟件(Peakall, 2006)計(jì)算長江刀鱭選育群體的有效等位基因數(shù)(N)、觀測(cè)雜合度(H)、期望雜合度(H)、Shannon多樣性指數(shù)()等遺傳多樣性參數(shù)。參照Bostein等(1980)的計(jì)算公式計(jì)算多態(tài)信息含量(Polymorphism information content, PIC):

    式中,PP分別為群體中第、個(gè)等位基因頻率,為等位基因數(shù)。

    2 結(jié)果

    2.1 刀鱭轉(zhuǎn)錄組中EST-SSR位點(diǎn)的數(shù)量與分布

    利用MISA軟件挖掘刀鱭轉(zhuǎn)錄組中的EST-SSR標(biāo)記(搜索結(jié)果見表1)。從刀鱭71869條Unigenes中共獲得33896個(gè)SSR位點(diǎn),占所評(píng)估序列總數(shù)目的47.16%;22463條Unigenes包含SSR位點(diǎn),占總Unigenes的31.25%;含有1個(gè)以上的SSR位點(diǎn)的Unigenes有7701個(gè),約占總Unigenes的10.71%;具有復(fù)合型SSR位點(diǎn)的Unigenes有2989個(gè),約占總Unigenes的4.16%(表1)。

    2.2 刀鱭EST-SSR的基序類型和頻率特征

    基于微衛(wèi)星重復(fù)單元的類型和重復(fù)次數(shù),將不同的SSR位點(diǎn)進(jìn)行分類(表2)。結(jié)果顯示,從單核苷酸重復(fù)SSR到六核苷酸重復(fù)SSR在刀鱭轉(zhuǎn)錄組中均有分布,但數(shù)量分布上具有不同程度的差異。其中,二核苷酸重復(fù)SSR數(shù)量最多,占SSR位點(diǎn)總數(shù)的49.47%;單核苷酸重復(fù)SSR數(shù)量次之,占SSR位點(diǎn)總數(shù)的34.94%;三核苷酸重復(fù)SSR的數(shù)量亦較為豐富,占SSR總位點(diǎn)數(shù)的13.34%;四核苷酸重復(fù)SSR的數(shù)量相對(duì)較少,占SSR總位點(diǎn)數(shù)的2.11%;五核苷酸和六核苷酸重復(fù)SSR的數(shù)量最少,分別占SSR總位點(diǎn)數(shù)的0.09%和0.06%。

    表1 刀鱭轉(zhuǎn)錄組中EST-SSR位點(diǎn)的搜索結(jié)果

    Tab.1 EST-SSR search results in the transcriptome of C. ectenes

    此外,刀鱭轉(zhuǎn)錄組EST-SSR中共有94種核苷酸重復(fù)基序,主要核苷酸重復(fù)基序亦有不同類型的優(yōu)勢(shì)重復(fù)基序。其中,單核苷酸重復(fù)SSR的優(yōu)勢(shì)重復(fù)基序?yàn)锳/T,占86.25%,這可能與Poly A被納入計(jì)算結(jié)果有關(guān);二核苷酸重復(fù)SSR的優(yōu)勢(shì)重復(fù)基序?yàn)锳C/GT,占75.25%;三核苷酸重復(fù)SSR的優(yōu)勢(shì)重復(fù)基序?yàn)锳GC/CTG和AGG/CCT,分別占27.16%和28.57%;四核苷酸重復(fù)SSR的優(yōu)勢(shì)重復(fù)基序?yàn)锳AAC/GTTT和ACAG/CTGT,分別占17.51%和16.81%。

    2.3 刀鱭EST-SSR不同基序重復(fù)次數(shù)分布

    刀鱭SSR主要重復(fù)基序的重復(fù)次數(shù)的分布情況見表3。刀鱭SSR基序的重復(fù)次數(shù)大部分為5~24次,其他重復(fù)次數(shù)相對(duì)較少,僅ATC/ATG出現(xiàn)了45次重復(fù)。低重復(fù)次數(shù)(5~9次重復(fù))的SSR位點(diǎn)最多,約占61.21%;中等重復(fù)次數(shù)(10~11次重復(fù))的SSR位點(diǎn)次之,約占22.84%;而>11次重復(fù)的SSR位點(diǎn)僅占15.95%。其中,單核苷酸基序重復(fù)次數(shù)均≥10,且重復(fù)11次以上的最多,占45.60%;二核苷酸基序以低重復(fù)次數(shù)為主,占92.26%,重復(fù)11次以上的較少,僅占0.02%;三核苷酸至六核苷酸基序均以重復(fù)次數(shù)為5的最多,其余重復(fù)次數(shù)的重復(fù)基序相對(duì)較少。

    2.4 刀鱭EST-SSR多態(tài)性引物的初步篩選和分析

    以本單位核心選育群體F3隨機(jī)選擇的30尾魚苗的基因組DNA為模板,針對(duì)100對(duì)四核苷酸重復(fù)SSR位點(diǎn)的引物進(jìn)行多態(tài)性檢驗(yàn),最終篩選到16對(duì)具有多態(tài)性的微衛(wèi)星引物(表4)。

    基于16對(duì)經(jīng)篩選獲得的多態(tài)性微衛(wèi)星標(biāo)記,對(duì)長江刀鱭核心選育群體(F3)進(jìn)行遺傳多樣性評(píng)估 (表5)。結(jié)果顯示,長江刀鱭核心選育群體(F3)的平均有效等位基因數(shù)(N)、平均觀測(cè)雜合度(H)、平均期望雜合度(H)和Shannon信息指數(shù)()分別為1.7580、0.3414、0.3977和0.6278?;诘任换蝾l率計(jì)算不同位點(diǎn)的多態(tài)信息含量,結(jié)果顯示,16對(duì)多態(tài)性微衛(wèi)星標(biāo)記中,有2對(duì)屬于高度多態(tài),10對(duì)屬于中度多態(tài),4對(duì)屬于低度多態(tài)??梢姡诘恩q轉(zhuǎn)錄組開發(fā)的多態(tài)性EST-SSR能夠應(yīng)用于刀鱭選育群體的種質(zhì)資源評(píng)估,具有較高實(shí)用性。

    表2 刀鱭轉(zhuǎn)錄組中EST-SSR的主要重復(fù)基序的類型和數(shù)量

    Tab.2 The type and number of main repeat motif of EST-SSR in the transcriptome of C. ectenes

    表3 刀鱭轉(zhuǎn)錄組主要EST-SSR重復(fù)基序的重復(fù)次數(shù)

    Tab.3 The repeat number of main repeat motif of EST-SSR in the transcriptome of C. ectenes

    表4 刀鱭多態(tài)性微衛(wèi)星引物

    Tab.4 Polymorphic SSR primers of C. ectenes

    3 討論

    微衛(wèi)星分子標(biāo)記在動(dòng)物進(jìn)化過程中,由于DNA分子的復(fù)制、缺失、滑移或錯(cuò)配及姐妹染色單體的不均交換等原因而具有高度變異性,被廣泛應(yīng)用于水產(chǎn)動(dòng)物群體遺傳結(jié)構(gòu)分析(滕爽爽等, 2018)、親子鑒定(張丹等, 2018)、遺傳圖譜構(gòu)建(郭香等, 2013)、數(shù)量性狀定位(顧穎等, 2016)和分子標(biāo)記輔助育種(孫際佳等, 2017; 顧穎等, 2016)等方面的研究。對(duì)于具有基因組信息的物種而言,微衛(wèi)星標(biāo)記的大規(guī)模開發(fā)相對(duì)方便;而對(duì)于目前還沒有基因組信息的物種,則可利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)大規(guī)模開發(fā)微衛(wèi)星標(biāo)記。刀鱭作為我國具有較高經(jīng)濟(jì)、營養(yǎng)和養(yǎng)殖價(jià)值的江海洄游性魚類,其基因組信息、轉(zhuǎn)錄組信息及可用分子標(biāo)記資源均極為匱乏,有鑒于全基因組測(cè)序的成本相對(duì)較高,因此,亟需開展刀鱭的轉(zhuǎn)錄學(xué)研究以獲得大量的可用分子標(biāo)記。此外,我國長江流域的刀鱭產(chǎn)量最為豐富,相應(yīng)的野生群體遺傳多樣性較高和資源分布較廣,是開展刀鱭遺傳育種研究的理想種質(zhì)資源庫。而本研究的刀鱭最初的繁育親本為長江流域的野生群體,歷經(jīng)3代高強(qiáng)度的人工選育,亟需適時(shí)地利用合適的分子標(biāo)記檢測(cè)其遺傳變異情況,以制定科學(xué)有效的育種策略。

    表5 16個(gè)微衛(wèi)星標(biāo)記的遺傳參數(shù)

    Tab.5 Genetic parameters at 16 microsatellite markers

    本研究基于刀鱭轉(zhuǎn)錄組數(shù)據(jù),利用生物信息學(xué)方法挖掘刀鱭EST-SSR,獲得了33896條刀鱭EST-SSR,這些EST-SSR的種類較為豐富,且1~6核苷酸重復(fù)基序均有不同程度的分布。其中,以出現(xiàn)頻率為49.47%的二核苷酸重復(fù)基序?yàn)橹?,其次是單核苷酸重?fù)基序(34.94%),三核苷酸重復(fù)基序亦較為豐富,占13.34%,其他類型的核苷酸重復(fù)基序相對(duì)較少。蔡磊等(2015)研究表明,多數(shù)物種EST-SSR中的重復(fù)基序均以二核苷酸重復(fù)基序?yàn)橹?,這與本研究結(jié)果一致;而其余類型的優(yōu)勢(shì)重復(fù)基序則有較大的差異,這可能與不同物種轉(zhuǎn)錄組文庫構(gòu)建的組織來源、不同物種的種間特異性、位點(diǎn)突變頻率及選擇性進(jìn)化機(jī)制有關(guān)。在刀鱭所有EST-SSR中,不同的核苷酸重復(fù)基序具有不同的重復(fù)次數(shù),其中,單核苷酸重復(fù)基序大部分屬于中高等重復(fù)次數(shù),其余類型重復(fù)基序的大部分重復(fù)次數(shù)在5~9之間,相對(duì)低于其他魚類基因組中的微衛(wèi)星重復(fù)次數(shù),可能與編碼區(qū)和非編碼區(qū)所受到的選擇壓力差異有關(guān)。不同重復(fù)基序的堿基優(yōu)勢(shì)重復(fù)類型表現(xiàn)出不同程度的偏倚性,單核苷酸重復(fù)基序以A/T類型為主,可能是由于富含A/T的微衛(wèi)星退火溫度相對(duì)較低,利于雙鏈DNA的解鏈,通過DNA復(fù)制的重組和滑動(dòng)機(jī)制而增加A/T重復(fù)基序的機(jī)率(倪守勝等, 2018)。而二核苷酸重復(fù)基序以AC/GT為主,這與魚類基因組微衛(wèi)星中以AC/GT重復(fù)為主相符;其余重復(fù)基序的優(yōu)勢(shì)重復(fù)堿基類型在不同物種間差異極大,可能與不同物種的種間差異性有關(guān)。此外,本研究發(fā)現(xiàn),刀鱭EST-SSR重復(fù)單元的拷貝次數(shù)隨著重復(fù)單元長度的增加而下降,微衛(wèi)星的數(shù)量亦呈現(xiàn)出逐步減少的趨勢(shì)。

    近年來,相關(guān)學(xué)者在魚類中利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)陸續(xù)開發(fā)出了多態(tài)性較高的微衛(wèi)星標(biāo)記。蔡磊等(2015)利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)從諸氏鯔蝦虎魚()肝臟轉(zhuǎn)錄組中獲得6225個(gè)微衛(wèi)星位點(diǎn),并從76對(duì)引物中篩選獲得32對(duì)具有多態(tài)性的微衛(wèi)星標(biāo)記。龔詩琦等(2016)利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)從黃姑魚()轉(zhuǎn)錄組獲得12254個(gè)微衛(wèi)星標(biāo)記,并從80對(duì)引物中篩選獲得18個(gè)具有多態(tài)性的微衛(wèi)星標(biāo)記。岳華梅等(2016)利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)從興國紅鯉()的垂體和性腺轉(zhuǎn)錄組中獲得13652個(gè)微衛(wèi)星標(biāo)記,并從30對(duì)引物中篩選獲得了9個(gè)具有多態(tài)性的微衛(wèi)星標(biāo)記。本研究從100對(duì)四堿基重復(fù)的微衛(wèi)星引物中篩選獲得了16對(duì)具有多態(tài)性的微衛(wèi)星標(biāo)記,引物多態(tài)率相對(duì)低于其他魚類的研究結(jié)果。造成多態(tài)性較低的原因,一方面是由于本研究選取的四堿基重復(fù)的微衛(wèi)星相對(duì)于二堿基和三堿基重復(fù)的微衛(wèi)星比較保守,且變異率低;另一方面是由于基于轉(zhuǎn)錄所獲得的EST-SSR擴(kuò)增產(chǎn)物可能跨越內(nèi)含子或者PCR引物位于內(nèi)含子和外顯子的結(jié)合處,而造成擴(kuò)增失敗導(dǎo)致多態(tài)率較低。此外,篩選多態(tài)性引物所用樣本的種類、數(shù)量、遺傳差異程度以及不同物種的DNA序列的保守性對(duì)于微衛(wèi)星位點(diǎn)多態(tài)性的篩選亦有不同程度的影響。

    基于Botstein等(1980)對(duì)于多態(tài)信息含量(PIC)的評(píng)估標(biāo)準(zhǔn),當(dāng)PIC>0.5、0.250.5、PIC<0.25時(shí)分別代表高度多態(tài)、中度多態(tài)和低度多態(tài)性。Rong等(2013)利用FIASCO法所開發(fā)的微衛(wèi)星標(biāo)記大部分為高度多態(tài)(70%),中度多態(tài)占25%,而本研究所獲得的多態(tài)性微衛(wèi)星標(biāo)記則大部分屬于中度多態(tài)(62.5%),高度多態(tài)僅占12.5%,表明不同的微衛(wèi)星開發(fā)方法所獲得的位點(diǎn)多態(tài)性有所不同,而實(shí)驗(yàn)樣本的來源、數(shù)量和位點(diǎn)的核苷酸的重復(fù)基序、次數(shù)的不同位點(diǎn)多態(tài)性亦有不同程度的影響。雜合度作為度量群體遺傳變異的重要參考指標(biāo),能夠有效反映出群體遺傳變異程度的高低。在本研究中,基于刀鱭多態(tài)性EST-SSR所獲得的刀鱭選育群體(F3)的平均觀測(cè)雜合度(H)和平均期望雜合度(H)分別為0.3414和0.3977,低于Rong等(2013)(H=0.49,H=0.68)、Ma等(2011)(H=0.675,H=0.650)、Chen等(2012)(H=0.53,H=0.77)利用磁珠富集法所獲的G-SSR的研究結(jié)果。究其原因,一方面是由于本研究選用的四核苷酸重復(fù)微衛(wèi)星序列在動(dòng)物進(jìn)化過程中相對(duì)保守,相對(duì)磁珠富集法所獲得的二核苷酸多態(tài)性和變異率較低,另一方面,本研究所選用的實(shí)驗(yàn)樣本經(jīng)過了3代高強(qiáng)度的人工選育,而其他研究樣本大部分為野生群體??傮w而言,基于轉(zhuǎn)錄組所獲得的EST-SSR多態(tài)性雖然相對(duì)低于基于基因組所獲得的G-SSR,但仍然能夠有效區(qū)分個(gè)體或者群體的遺傳差異,加上所開發(fā)的微衛(wèi)星數(shù)量相對(duì)較多,具有相對(duì)較高的實(shí)用性。

    本研究利用生物信息學(xué)方法篩選刀鱭高通量轉(zhuǎn)錄組Unigenes中的微衛(wèi)星標(biāo)記,分析刀鱭EST-SSR分布和組成特征,并利用所開發(fā)的多態(tài)性微衛(wèi)星標(biāo)記初步分析了刀鱭核心選育群體(F3)的遺傳特征。研究結(jié)果不僅豐富了刀鱭分子標(biāo)記的數(shù)量,而且對(duì)刀鱭后續(xù)選育群體的育種效果評(píng)估及分子標(biāo)記輔助育種實(shí)踐具有重要的指導(dǎo)意義。

    Bai CC, Liu SF, Zhuang ZM. Characteristic analysis of microsatellite DNA in the genome of Gobiidae. Progress in Fishery Sciences, 2016, 37(5): 9–15 [白翠翠, 柳淑芳, 莊志猛. 蝦虎魚科(Gobiidae)基因組微衛(wèi)星DNA的分布特征. 漁業(yè)科學(xué)進(jìn)展, 2016, 37(5): 9–15]

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

    Cai L, Yu LJ, Chen XQ,. A preliminary screening and characteristic analysis of microsatellite markers from transcriptome sequences in. Biotechnology Bulletin, 2015(9): 146–151 [蔡磊, 余露軍, 陳小曲, 等. 諸氏鯔蝦虎魚轉(zhuǎn)錄組序列中微衛(wèi)星標(biāo)記的初步篩選及特征分析. 生物技術(shù)通報(bào), 2015(9): 146–151]

    Chen F, Zeng LY, Cheng QQ. Development of thirty-four novel polymorphic microsatellite markers in(Clupeiformes: Engraulidae) and cross-species amplification in two closely related taxa. Journal of Genetics, 2012, 91(1): 1–7

    Cheng QQ, Zhang QY, Ma CY,. Genetic structure and differentiation of four lake populations of(Clupeiformes: Engraulidae) based on mtDNA control region sequences. Biochemical Systematics and Ecology, 2011, 39(4–6): 544–552

    Deng PP, Shi YH, Zhang GY,. Isolation of microsatellite inby magnetic beads. Biotechnology Bulletin, 2014(6): 145–149 [鄧平平, 施永海, 張根玉, 等. 磁珠富集法分離刀鱭微衛(wèi)星標(biāo)記. 生物技術(shù)通報(bào), 2014(6): 145–149]

    Ge GC, Cao T, Chen CJ,. Analysis of genetic diversity based on amplified fragment length polymorphism fingerprint offrom Yangtze River. Journal of Nanjing University (Natural Sciences), 2008, 44(3): 332–338 [葛家春, 曹廷, 陳嬋娟, 等. 利用擴(kuò)增片斷長度多態(tài)性技術(shù)分析長江刀鱭的遺傳多樣性. 南京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008, 44(3): 332–338]

    Gong SQ, Wang ZY, Xiao SJ,. Development and verification of SSR based on transcriptome of yellow drum,. Journal of Jimei University (Natural Science), 2016, 21(4): 241–246 [龔詩琦, 王志勇, 肖世俊, 等. 黃姑魚轉(zhuǎn)錄組SSR的開發(fā)與驗(yàn)證. 集美大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 21(4): 241–246]

    Gu Y, Lu CY, Zhang Q,. Analysis of genetic diversity and growth traits inusing microsatellite QTL markers. Freshwater Fisheries, 2016, 46(4): 9–18 [顧穎, 魯翠云, 張芹, 等. 微衛(wèi)星QTL標(biāo)記分析豫選黃河鯉群體遺傳結(jié)構(gòu)及生長性狀相關(guān)性. 淡水漁業(yè), 2016, 46(4): 9–18]

    Guan A, Wu YT, Chen Y,. Deep sequence-based transcriptome analysis of microsatellites in the cuttlefish (). Progress in Fishery Sciences, 2018, 39(3): 144–151 [管奧, 毋玉婷, 陳宇, 等. 曼氏無針烏賊轉(zhuǎn)錄組微衛(wèi)星特征分析. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(3): 144–151]

    Guo X, Li Q, Kong LF,. Microsatellite-based consensus map of the Pacific oyster (). Journal of Fisheries of China, 2013, 37(6): 823–829 [郭香, 李琪, 孔令鋒, 等. 基于微衛(wèi)星標(biāo)記整合長牡蠣遺傳圖譜. 水產(chǎn)學(xué)報(bào), 2013, 37(6): 823–829]

    Li DM, Yang AG, Wu B,. Development and application of the EST-SSR markers in. Progress in Fishery Sciences, 2017, 38(2): 137–142 [李東明, 楊愛國, 吳彪, 等. 櫛江珧()EST-SSR標(biāo)記的開發(fā)與應(yīng)用. 漁業(yè)科學(xué)進(jìn)展, 2017, 38(2): 137–142]

    Ma CY, Cheng QQ, Zhang QY. Genetic variation of(Clupeiformes: Engraulidae) revealed by the complete cytochrome b sequences of mitochondrial DNA. Journal of Experimental Marine Biology and Ecology, 2010, 385(1): 14–19

    Ma CY, Cheng QQ, Zhang QY. Development of 12 polymorphic microsatellite markers inJordan and Seale, 1905 (Clupeiformes: Engraulidae) and cross-species amplification inLinnaeus, 1758. Environmental Biology of Fishes, 2011, 91(2): 243–249

    Ma CY, Liu M, Ma LB,. Genetic diversity inby RAPD analysis. Marine Fisheries Research, 2004, 25(5): 19–24 [馬春艷, 劉敏, 馬凌波, 等. 長江口刀鱭遺傳多樣性的隨機(jī)擴(kuò)增多態(tài)DNA(RAPD)分析. 海洋水產(chǎn)研究, 2004, 25(5): 19–24]

    Ni SS, Yang Y, Liu SF,. Microsatellite analysis ofusing next-generation sequencing method. Progress in Fishery Sciences, 2018, 39(1): 107–113 [倪守勝, 楊鈺, 柳淑芳, 等. 基于高通量測(cè)序的蝦夷扇貝基因組微衛(wèi)星特征分析. 漁業(yè)科學(xué)進(jìn)展, 2018, 39(1): 107–113]

    Peakall R, Smouse PE. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Resources, 2006, 6(1): 288–295

    Rong XJ, Xu YJ, Wang QY,. Isolation and characterization of polymorphic microsatellite markers from. Genetics and Molecular Research, 2013, 12(4): 6011–6017

    Shi YH, Zhang GY, Zhang HM,. Analysis and comparison of nutritional quality of muscle fromFED formulated feed and live feed. Chinese Journal of Animal Nutrition, 2014, 26(2): 427–436 [施永海, 張根玉, 張海明, 等. 配合飼料和活餌料喂養(yǎng)刀鱭肌肉營養(yǎng)品質(zhì)分析與比較. 動(dòng)物營養(yǎng)學(xué)報(bào), 2014, 26(2): 427–436]

    Shi YH, Zhang GY, Zhang HM,. Complete artificial propagation and embryonic development of. Journal of Shanghai Ocean University, 2015, 24(1): 36–43 [施永海, 張根玉, 張海明, 等. 刀鱭的全人工繁殖及胚胎發(fā)育. 上海海洋大學(xué)學(xué)報(bào), 2015, 24(1): 36–43]

    Sun JJ, Li GF, Liu L,. Analysis of microsatellite DNA markers emphasis on correlation with some economically important traits in Mandarin fish. Chinese Journal of Fisheries, 2017, 30(1): 11–18 [孫際佳, 李桂峰, 劉麗, 等. 翹嘴鱖微衛(wèi)星標(biāo)記及其與主要經(jīng)濟(jì)性狀的相關(guān)分析. 水產(chǎn)學(xué)雜志, 2017, 30(1): 11–18]

    Teng SS, Fang J, Shao YQ,. Microsatellite analysis of genetic variation in the fast growth families of the second generation ofand correlation with growth traits. Acta Hydrobiologica Sinica, 2018, 42(4): 681–689 [滕爽爽, 方軍, 邵艷卿, 等. 泥蚶G2代快速生長家系遺傳結(jié)構(gòu)的微衛(wèi)星分析及其與生長性狀的關(guān)聯(lián). 水生生物學(xué)報(bào), 2018, 42(4): 681–689]

    Wei GL, Xu GC, Gu RB,. Studies on the genetic diversity of farmed and wild populations ofby analysing mitochondrial DNA Cytb genes. Acta Agricultural Universitatis Jiangxiensis, 2012, 34(6): 1216–1221 [魏廣蓮, 徐鋼春, 顧若波, 等. 基于mtDNA Cytb序列分析養(yǎng)殖與野生刀鱭群體的遺傳多樣性. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào), 2012, 34(6): 1216–1221]

    Wei GL, Xu GC, Gu RB,. Biology and artificial breeding progress in. Journal of Yangtze University (Natural Science), 2013, 9(7): 31–35 [魏廣蓮, 徐鋼春, 顧若波, 等. 刀鱭的生物學(xué)及人工養(yǎng)殖研究進(jìn)展. 長江大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 9(7): 31–35]

    Yang M, Yu Y, Zhang XJ,. 2017. Development of microsatellite markers from the transcriptome sequences of Pacific white shrimp (). Marine Sciences, 41(2): 96–102 [楊銘, 于洋, 張曉軍, 等. 基于轉(zhuǎn)錄組數(shù)據(jù)的凡納濱對(duì)蝦微衛(wèi)星標(biāo)記開發(fā). 海洋科學(xué), 2017, 41(2): 96–102]

    Yuan CM. The resources ofin middle and lower Yangtze River, and the status and the causation of the population formation changes. Chinese Journal of Zoology, 1988, 23(3): 12–15 [袁傳宓. 長江中下游刀鱭資源和種群組成變動(dòng)狀況及其原因. 動(dòng)物學(xué)雜志, 1988, 23(3): 12–15]

    Yue HM, Zhai Q, Song MY,. Development of microsatellite markers inusing next- generation sequencing. Freshwater Fisheries, 2016, 46(1): 24–28 [岳華梅, 翟晴, 宋明月, 等. 基于轉(zhuǎn)錄組測(cè)序的興國紅鯉微衛(wèi)星標(biāo)記篩選. 淡水漁業(yè), 2016, 46(1): 24–28]

    Zhang D, Fu JJ, Zhang LD,. The parentage analysis of bighead carp () based on ten microsatellite markers. Genomics and Applied Biology, 2018, 1–11 [張丹, 傅建軍, 張利德, 等. 鳙基于10個(gè)微衛(wèi)星標(biāo)記的親子鑒定分析. 基因組學(xué)與應(yīng)用生物學(xué), 2018, 1–11]

    Zhang YP, He G, Wang S,. Analysis of genetic diversity ofandbased on mitochondrial DNA D-loop sequence. Jiangsu Agricultural Sciences, 2017, 45(20): 59–64 [張燕萍, 賀剛, 王生, 等. 基于線粒體DNA D-loop序列分析刀鱭與短頜鱭的遺傳多樣性. 江蘇農(nóng)業(yè)科學(xué), 2017, 45(20): 59–64]

    Characteristic Analysis of Microsatellites in SelectedUsing a Transcriptome Dataset

    YU Aiqing, SHI Yonghai①, XU Jiabo, LU Genhai, ZHANG Haiming, XIE Yongde, LIU Yongshi

    (Shanghai Fisheries Research Institute, Shanghai Fisheries Technology Extension Station, Shanghai 200433)

    The objective of this study was to identify microsatellites from the transcriptome sequences ofusing MISA software. These microsatellites could create very useful resources in the evaluation of germplasm resource and marker-assisted breeding of. The results showed that a total of 33,896 microsatellites were identified with repeating units, at lengths of 1~6 bases, from 71,869 unigenes. Different types of repeat SSRs had considerably different distribution characteristics. The majority of the microsatellite loci consisted of mono-, di-, and tri-nucleotide motifs (34.94%, 49.47%, and 13.34%, respectively). Dinucleotide microsatellite repeating units were the most abundant in thetranscriptome, and the AC/GT repeating units were the most ascendant repeating unit (75.25%). The dominant repeating units for the mononucleotide and trinucleotide motifs were A/T (86.25%) and AGG/CCT (28.57%), respectively. Different nucleotide repeat motifs had different repetitions that had a reducing trend, with the increase in the number of nucleotides in repeat motifs. Among the 100 designed primer pairs, 16 pairs proved to be polymorphic microsatellite markers. In the present study, 16 polymorphic microsatellite loci were characterized to evaluate the genetic diversity of the selected(F3). The results showed that the average number of effective alleles (N), the average observed heterozygosity (H), the average expected heterozygosity (H),and the average Shannon′s information index () of the selected(F3) were 1.7580, 0.3414, 0.3977, and 0.6278, respectively. These results indicate that it is feasible to develop microsatellite markers based on thetranscriptome and polymorphic microsatellite loci obtained in this study will facilitate further studies on the population genetic management and conservation of.

    ; Transcriptome; Microsatellites; Molecular markers

    SHI Yonghai, E-mail: yonghais@163.com

    * 上海市科技興農(nóng)重點(diǎn)攻關(guān)項(xiàng)目[滬農(nóng)科攻字(2016)第6-2-2號(hào)]、上海市科學(xué)技術(shù)委員會(huì)重點(diǎn)科技攻關(guān)項(xiàng)目(17391900300)和上海市科學(xué)技術(shù)委員會(huì)重點(diǎn)科技攻關(guān)項(xiàng)目(11391901300)共同資助 [This work was supported by the Shanghai Key Scientific and Technological Project on Agriculture from Shanghai Municipal Agricultural Commission (2016-2-2), Shanghai Key Scientific Technological Project from Shanghai Science and Technology Commission (17391900300), and Shanghai Key Scientific Technological Project from Shanghai Science and Technology Commission (11391901300)]. 于愛清,E-mail: aiqingyu0125@163.com

    施永海,教授級(jí)高級(jí)工程師,E-mail: yonghais@163.com

    2018-08-02,

    2018-08-29

    S937.3

    A

    2095-9869(2019)05-0101-09

    10.19663/j.issn2095-9869.20180802001

    http://www.yykxjz.cn/

    于愛清, 施永海, 徐嘉波, 陸根海, 張海明, 謝永德. 劉永士. 長江刀鱭選育群體轉(zhuǎn)錄組EST-SSR的分布特征分析. 漁業(yè)科學(xué)進(jìn)展, 2019, 40(5): 101–109

    Yu AQ, Shi YG, Xu JB, Lu GH, Zhang HM, Xie YD, Liu YS. Characteristic analysis of microsatellites in selectedusing a transcriptome dataset. Progress in Fishery Sciences, 2019, 40(5): 101–109

    (編輯 馮小花)

    猜你喜歡
    基序微衛(wèi)星核苷酸
    EPIYA 基序與幽門螺桿菌感染相關(guān)胃病關(guān)系的研究進(jìn)展
    梅花鹿基因組微衛(wèi)星分布特征研究
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進(jìn)展
    徐長風(fēng):核苷酸類似物的副作用
    肝博士(2022年3期)2022-06-30 02:48:28
    帶TRS基序突變的新型冠狀病毒威脅更大
    芥藍(lán)Aux/IAA家族基因生物信息學(xué)與表達(dá)分析
    Acknowledgment to reviewers—November 2018 to September 2019
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國“一箭雙星”成功將“遙感衛(wèi)星二十一號(hào)”與“天拓二號(hào)視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    韩国av在线不卡| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美一区二区三区黑人 | 99热国产这里只有精品6| tube8黄色片| 丝袜在线中文字幕| 老司机影院成人| 在线观看三级黄色| 亚洲精品,欧美精品| 一级毛片 在线播放| 亚洲av.av天堂| 人人妻人人添人人爽欧美一区卜| 看免费av毛片| 蜜桃国产av成人99| 不卡视频在线观看欧美| 一级毛片黄色毛片免费观看视频| 18+在线观看网站| 亚洲国产最新在线播放| 久久ye,这里只有精品| 两个人看的免费小视频| 制服人妻中文乱码| 国产精品亚洲av一区麻豆 | 一级,二级,三级黄色视频| 成人黄色视频免费在线看| 成人国语在线视频| 夫妻午夜视频| 国产片特级美女逼逼视频| 老女人水多毛片| 欧美日本中文国产一区发布| 国产欧美日韩综合在线一区二区| 国产免费一区二区三区四区乱码| 91成人精品电影| 99精国产麻豆久久婷婷| 久久久久国产一级毛片高清牌| 午夜91福利影院| 久久久久久人妻| 黄色配什么色好看| 亚洲伊人久久精品综合| 在线观看免费高清a一片| 免费看av在线观看网站| 国产成人午夜福利电影在线观看| 美女主播在线视频| 少妇猛男粗大的猛烈进出视频| 国产精品久久久久久精品古装| 男女无遮挡免费网站观看| 亚洲欧美精品自产自拍| 欧美另类一区| 国产黄色视频一区二区在线观看| 亚洲精华国产精华液的使用体验| 一区二区三区精品91| 日本猛色少妇xxxxx猛交久久| 久久国内精品自在自线图片| 久久婷婷青草| 免费人妻精品一区二区三区视频| 久久久久久久精品精品| 国产一级毛片在线| 在线精品无人区一区二区三| 欧美bdsm另类| 咕卡用的链子| 欧美在线黄色| 精品国产一区二区三区久久久樱花| 精品一区二区免费观看| 久久午夜福利片| 国产在线一区二区三区精| 亚洲人成77777在线视频| 日韩不卡一区二区三区视频在线| 又黄又粗又硬又大视频| 精品人妻在线不人妻| 男的添女的下面高潮视频| 一边亲一边摸免费视频| 久久精品夜色国产| 欧美日本中文国产一区发布| 极品人妻少妇av视频| 好男人视频免费观看在线| 久久精品国产鲁丝片午夜精品| 色94色欧美一区二区| 亚洲国产av新网站| 国产片特级美女逼逼视频| 久久婷婷青草| 国产 一区精品| 我的亚洲天堂| 久久精品国产鲁丝片午夜精品| 人体艺术视频欧美日本| 亚洲精品国产一区二区精华液| 欧美日韩精品网址| 亚洲人成网站在线观看播放| 美女主播在线视频| 又黄又粗又硬又大视频| 精品视频人人做人人爽| 亚洲天堂av无毛| 最近2019中文字幕mv第一页| 18禁观看日本| 中文字幕另类日韩欧美亚洲嫩草| 国产精品欧美亚洲77777| 侵犯人妻中文字幕一二三四区| 赤兔流量卡办理| 亚洲人成电影观看| 美女福利国产在线| 亚洲第一av免费看| av在线观看视频网站免费| 国产高清国产精品国产三级| 日韩 亚洲 欧美在线| 国产日韩欧美亚洲二区| 9热在线视频观看99| 边亲边吃奶的免费视频| av福利片在线| 91久久精品国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 亚洲国产成人一精品久久久| 亚洲第一青青草原| 老司机影院毛片| 只有这里有精品99| 亚洲精品国产一区二区精华液| 久久 成人 亚洲| 中文乱码字字幕精品一区二区三区| 亚洲伊人色综图| 国产野战对白在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 女的被弄到高潮叫床怎么办| 99re6热这里在线精品视频| 国产极品天堂在线| 亚洲激情五月婷婷啪啪| 日本色播在线视频| 婷婷色麻豆天堂久久| 爱豆传媒免费全集在线观看| 精品少妇久久久久久888优播| 亚洲精品日本国产第一区| 天天躁狠狠躁夜夜躁狠狠躁| 高清视频免费观看一区二区| 午夜福利网站1000一区二区三区| 在线看a的网站| av免费在线看不卡| 国产一区二区 视频在线| 春色校园在线视频观看| 欧美国产精品va在线观看不卡| 老司机影院成人| 乱人伦中国视频| 亚洲,欧美,日韩| 精品福利永久在线观看| 亚洲欧洲国产日韩| 国产男女内射视频| 国产精品二区激情视频| av有码第一页| 欧美亚洲 丝袜 人妻 在线| 最新中文字幕久久久久| 免费久久久久久久精品成人欧美视频| 久久97久久精品| 亚洲综合精品二区| 中文字幕制服av| 一区二区日韩欧美中文字幕| 久久人妻熟女aⅴ| 亚洲五月色婷婷综合| 伊人久久国产一区二区| 色婷婷久久久亚洲欧美| 观看美女的网站| 国产男人的电影天堂91| 国产精品av久久久久免费| 91精品三级在线观看| 久久久久久久亚洲中文字幕| 亚洲成人一二三区av| 丰满少妇做爰视频| 成人18禁高潮啪啪吃奶动态图| 色吧在线观看| 久久久精品国产亚洲av高清涩受| a级毛片在线看网站| 久久久久久久久久人人人人人人| 大香蕉久久网| 久久久久精品人妻al黑| 黄色配什么色好看| 狠狠精品人妻久久久久久综合| 亚洲欧美日韩另类电影网站| 午夜福利影视在线免费观看| 青青草视频在线视频观看| 人人妻人人添人人爽欧美一区卜| 国产人伦9x9x在线观看 | 午夜影院在线不卡| av在线老鸭窝| 日日啪夜夜爽| 极品少妇高潮喷水抽搐| 韩国精品一区二区三区| 韩国高清视频一区二区三区| 国产免费又黄又爽又色| 晚上一个人看的免费电影| 日本欧美国产在线视频| 久久午夜福利片| 高清不卡的av网站| 亚洲情色 制服丝袜| 免费久久久久久久精品成人欧美视频| 精品久久蜜臀av无| 各种免费的搞黄视频| 少妇的丰满在线观看| 国产97色在线日韩免费| 成人国产av品久久久| 一区二区三区乱码不卡18| 久久午夜福利片| 国产精品久久久久成人av| 亚洲国产最新在线播放| av国产久精品久网站免费入址| 日韩成人av中文字幕在线观看| 国产1区2区3区精品| 嫩草影院入口| 天天躁日日躁夜夜躁夜夜| 日韩免费高清中文字幕av| 热99国产精品久久久久久7| 老司机影院成人| 国产 精品1| 精品少妇一区二区三区视频日本电影 | 蜜桃国产av成人99| 纯流量卡能插随身wifi吗| 国产成人欧美| 成年av动漫网址| 丝袜喷水一区| 久久久国产欧美日韩av| 国产成人精品一,二区| 久久狼人影院| 免费看不卡的av| 中文字幕人妻丝袜一区二区 | 成人午夜精彩视频在线观看| 寂寞人妻少妇视频99o| 亚洲欧美精品自产自拍| av.在线天堂| 男女下面插进去视频免费观看| 天天躁夜夜躁狠狠久久av| 两性夫妻黄色片| 亚洲国产av影院在线观看| 少妇人妻久久综合中文| 精品久久久久久电影网| 亚洲精品乱久久久久久| 国产精品免费大片| 99久久中文字幕三级久久日本| 国产熟女午夜一区二区三区| 黄频高清免费视频| 亚洲一码二码三码区别大吗| 精品人妻偷拍中文字幕| 久久99蜜桃精品久久| 国产国语露脸激情在线看| 久久久久视频综合| 有码 亚洲区| 黄片播放在线免费| 美女福利国产在线| 2021少妇久久久久久久久久久| 丝袜脚勾引网站| 激情视频va一区二区三区| 乱人伦中国视频| h视频一区二区三区| 高清视频免费观看一区二区| videosex国产| 交换朋友夫妻互换小说| 国产黄频视频在线观看| 日韩精品有码人妻一区| 久久久久视频综合| 美女xxoo啪啪120秒动态图| 久久韩国三级中文字幕| 2021少妇久久久久久久久久久| 久久人人爽人人片av| 韩国精品一区二区三区| 成年女人在线观看亚洲视频| 亚洲精品一二三| 久久久久久久久久人人人人人人| 精品福利永久在线观看| 欧美日韩精品成人综合77777| 久久久久久伊人网av| 夫妻午夜视频| 亚洲精品成人av观看孕妇| 人妻少妇偷人精品九色| 少妇熟女欧美另类| 九九爱精品视频在线观看| 女人高潮潮喷娇喘18禁视频| 日本欧美视频一区| 捣出白浆h1v1| 成人毛片60女人毛片免费| 哪个播放器可以免费观看大片| av免费在线看不卡| 黑人欧美特级aaaaaa片| 午夜福利影视在线免费观看| 丰满饥渴人妻一区二区三| 亚洲av在线观看美女高潮| 国产午夜精品一二区理论片| 黄色配什么色好看| 国产福利在线免费观看视频| 中文欧美无线码| 久久这里只有精品19| 久久ye,这里只有精品| videossex国产| 青春草国产在线视频| tube8黄色片| 99国产精品免费福利视频| 最近最新中文字幕大全免费视频 | 久久国内精品自在自线图片| 另类精品久久| 在线观看一区二区三区激情| 亚洲精品美女久久av网站| 亚洲精品aⅴ在线观看| 欧美日韩精品成人综合77777| 国产精品免费大片| 99九九在线精品视频| 亚洲精品成人av观看孕妇| 天天躁狠狠躁夜夜躁狠狠躁| 欧美人与善性xxx| 成年美女黄网站色视频大全免费| 91国产中文字幕| 2018国产大陆天天弄谢| 少妇猛男粗大的猛烈进出视频| 亚洲欧美精品综合一区二区三区 | 大香蕉久久网| 国产黄色免费在线视频| 少妇的丰满在线观看| 免费观看无遮挡的男女| 亚洲欧美精品综合一区二区三区 | 午夜福利,免费看| 久久99一区二区三区| 日本wwww免费看| 在线观看国产h片| 春色校园在线视频观看| 亚洲精品国产av蜜桃| 欧美变态另类bdsm刘玥| 国产av精品麻豆| 观看av在线不卡| 欧美另类一区| 婷婷色综合www| 九九爱精品视频在线观看| 午夜免费男女啪啪视频观看| 成年美女黄网站色视频大全免费| 免费在线观看黄色视频的| 青春草亚洲视频在线观看| 欧美人与性动交α欧美软件| 18禁观看日本| 亚洲,欧美精品.| 母亲3免费完整高清在线观看 | 黑人欧美特级aaaaaa片| 久久ye,这里只有精品| 午夜福利视频精品| 一区二区三区四区激情视频| 欧美日韩精品网址| 精品亚洲成a人片在线观看| 亚洲在久久综合| av在线观看视频网站免费| 男女边吃奶边做爰视频| 午夜免费鲁丝| 亚洲三区欧美一区| 男人添女人高潮全过程视频| 精品国产露脸久久av麻豆| 美女主播在线视频| 午夜激情久久久久久久| 美国免费a级毛片| 亚洲av综合色区一区| 国产不卡av网站在线观看| 午夜福利在线观看免费完整高清在| 久久久久精品久久久久真实原创| 亚洲情色 制服丝袜| 久久人人爽人人片av| 亚洲精品久久成人aⅴ小说| 肉色欧美久久久久久久蜜桃| 日韩不卡一区二区三区视频在线| 黑人欧美特级aaaaaa片| 亚洲欧美日韩另类电影网站| 一区福利在线观看| 国产女主播在线喷水免费视频网站| 七月丁香在线播放| 午夜日韩欧美国产| 久久久久久久国产电影| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 婷婷色综合大香蕉| 天美传媒精品一区二区| 在线观看www视频免费| 亚洲精品一二三| 成人毛片a级毛片在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 大片免费播放器 马上看| 熟妇人妻不卡中文字幕| a级片在线免费高清观看视频| av网站免费在线观看视频| 最近最新中文字幕免费大全7| 在线观看美女被高潮喷水网站| 国产亚洲一区二区精品| 黄片小视频在线播放| 午夜福利在线免费观看网站| 如日韩欧美国产精品一区二区三区| 美女xxoo啪啪120秒动态图| 秋霞在线观看毛片| av视频免费观看在线观看| 成人毛片60女人毛片免费| 日韩中文字幕视频在线看片| 国产片特级美女逼逼视频| 国产精品麻豆人妻色哟哟久久| 久久精品人人爽人人爽视色| 欧美人与性动交α欧美精品济南到 | 婷婷色综合大香蕉| 久久久久久久精品精品| 亚洲三区欧美一区| 国产精品不卡视频一区二区| 这个男人来自地球电影免费观看 | 亚洲av福利一区| 精品久久久精品久久久| 三上悠亚av全集在线观看| 看免费成人av毛片| 青草久久国产| 亚洲国产精品一区三区| 少妇 在线观看| 欧美变态另类bdsm刘玥| 精品人妻偷拍中文字幕| 国产男女内射视频| 国产精品熟女久久久久浪| av在线老鸭窝| 国产成人精品久久久久久| 18禁裸乳无遮挡动漫免费视频| 男男h啪啪无遮挡| 久久人人97超碰香蕉20202| 女人久久www免费人成看片| 黄片无遮挡物在线观看| 成人漫画全彩无遮挡| 午夜激情av网站| 日韩欧美一区视频在线观看| 国产成人精品无人区| 欧美精品国产亚洲| 国产老妇伦熟女老妇高清| 99久久人妻综合| 国产在线视频一区二区| 69精品国产乱码久久久| 中文字幕最新亚洲高清| 国产成人aa在线观看| 国产色婷婷99| 亚洲精品aⅴ在线观看| 国产精品免费大片| 97人妻天天添夜夜摸| 亚洲国产看品久久| 男女免费视频国产| 99久久中文字幕三级久久日本| 嫩草影院入口| 国产精品麻豆人妻色哟哟久久| 免费观看在线日韩| 精品久久蜜臀av无| 亚洲一区中文字幕在线| 日韩一区二区视频免费看| 美女大奶头黄色视频| 制服人妻中文乱码| 在线免费观看不下载黄p国产| 国产精品国产av在线观看| 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 女人被躁到高潮嗷嗷叫费观| videosex国产| 一级毛片 在线播放| 男人爽女人下面视频在线观看| 一区福利在线观看| 视频在线观看一区二区三区| 久久热在线av| 丝袜在线中文字幕| 成人亚洲精品一区在线观看| 日日撸夜夜添| 日韩欧美精品免费久久| av在线老鸭窝| 日韩一本色道免费dvd| 中国国产av一级| 亚洲五月色婷婷综合| 欧美人与善性xxx| 国产欧美日韩综合在线一区二区| 国产av码专区亚洲av| 成人国语在线视频| av女优亚洲男人天堂| 亚洲视频免费观看视频| 亚洲熟女精品中文字幕| 国产精品久久久久久av不卡| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 亚洲欧洲精品一区二区精品久久久 | 美女高潮到喷水免费观看| 啦啦啦在线观看免费高清www| a 毛片基地| 人妻一区二区av| 亚洲欧洲国产日韩| 日韩免费高清中文字幕av| 精品一区二区三卡| 欧美日韩成人在线一区二区| 中文字幕av电影在线播放| 国产亚洲精品第一综合不卡| 寂寞人妻少妇视频99o| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲五月色婷婷综合| 亚洲国产日韩一区二区| 99久久中文字幕三级久久日本| 老司机影院成人| 亚洲在久久综合| 国产成人91sexporn| 亚洲国产色片| 久久亚洲国产成人精品v| 亚洲精品日韩在线中文字幕| 午夜免费鲁丝| videos熟女内射| 99国产综合亚洲精品| 亚洲国产欧美在线一区| 国产黄色视频一区二区在线观看| 精品少妇黑人巨大在线播放| 国产成人免费观看mmmm| 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放| 成人国语在线视频| 麻豆乱淫一区二区| 欧美日韩一区二区视频在线观看视频在线| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 精品国产一区二区三区久久久樱花| 亚洲综合色网址| 亚洲精品中文字幕在线视频| 精品国产露脸久久av麻豆| 亚洲精品国产色婷婷电影| 国产精品久久久久久久久免| 人体艺术视频欧美日本| 国产片内射在线| 韩国精品一区二区三区| 寂寞人妻少妇视频99o| 看免费av毛片| 巨乳人妻的诱惑在线观看| 国产精品国产三级专区第一集| 亚洲婷婷狠狠爱综合网| 一区福利在线观看| 亚洲欧美成人综合另类久久久| 国产精品嫩草影院av在线观看| 黄片小视频在线播放| av免费在线看不卡| 天堂俺去俺来也www色官网| 国产av精品麻豆| 制服人妻中文乱码| 亚洲伊人色综图| 中文天堂在线官网| 久久青草综合色| 久久狼人影院| 久久久久精品久久久久真实原创| 国产片特级美女逼逼视频| 国产一区二区三区综合在线观看| 中文字幕色久视频| 亚洲三区欧美一区| 免费观看在线日韩| 91在线精品国自产拍蜜月| 巨乳人妻的诱惑在线观看| 青草久久国产| 女性被躁到高潮视频| 久久久久国产网址| 亚洲国产欧美在线一区| 中文字幕制服av| 亚洲精品自拍成人| 最近手机中文字幕大全| 免费大片黄手机在线观看| 色哟哟·www| 在线观看一区二区三区激情| 欧美日韩精品成人综合77777| 国产精品免费视频内射| 一级毛片黄色毛片免费观看视频| 我的亚洲天堂| 成年女人在线观看亚洲视频| 亚洲成国产人片在线观看| 国产爽快片一区二区三区| 午夜日本视频在线| 热99国产精品久久久久久7| 日韩熟女老妇一区二区性免费视频| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三卡| 久久这里只有精品19| 香蕉精品网在线| 亚洲人成电影观看| 色婷婷久久久亚洲欧美| 一级毛片黄色毛片免费观看视频| 中文字幕av电影在线播放| 亚洲视频免费观看视频| www.精华液| 国产日韩欧美在线精品| av福利片在线| 美女视频免费永久观看网站| 亚洲精品日韩在线中文字幕| 日本欧美国产在线视频| 久久国产精品大桥未久av| 久久人人爽av亚洲精品天堂| 99久久人妻综合| 日本午夜av视频| 美女视频免费永久观看网站| 日产精品乱码卡一卡2卡三| 9热在线视频观看99| 久久人妻熟女aⅴ| 丰满少妇做爰视频| videosex国产| 免费观看无遮挡的男女| 99re6热这里在线精品视频| 激情五月婷婷亚洲| 国产免费又黄又爽又色| 免费高清在线观看视频在线观看| a级毛片黄视频| 中文字幕人妻丝袜制服| 久久精品国产综合久久久| 久久女婷五月综合色啪小说| 久久97久久精品| 人妻一区二区av| 精品第一国产精品| 亚洲精华国产精华液的使用体验| 国产精品 欧美亚洲| 巨乳人妻的诱惑在线观看| 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 少妇被粗大猛烈的视频| 久热久热在线精品观看| 欧美日韩综合久久久久久| 国产免费一区二区三区四区乱码| 看十八女毛片水多多多| 国产1区2区3区精品| 久久精品国产亚洲av涩爱| 美女xxoo啪啪120秒动态图| 1024香蕉在线观看| 97在线视频观看| 国产又色又爽无遮挡免| 男女国产视频网站| 久久久欧美国产精品| 亚洲欧洲日产国产| 亚洲av在线观看美女高潮| 搡女人真爽免费视频火全软件| 亚洲av综合色区一区| 男人操女人黄网站|