• <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
    亚洲自偷自拍三级| 欧美3d第一页| 日韩一区二区视频免费看| 国产男女超爽视频在线观看| 免费播放大片免费观看视频在线观看| 女人十人毛片免费观看3o分钟| 亚洲精品成人av观看孕妇| 99热这里只有是精品50| 色网站视频免费| 久久国产乱子免费精品| 日韩成人av中文字幕在线观看| 久久久久网色| 黄色欧美视频在线观看| 亚洲欧美一区二区三区黑人 | 久久久久精品性色| 日韩视频在线欧美| 插逼视频在线观看| 日本wwww免费看| 啦啦啦中文免费视频观看日本| 亚洲欧美一区二区三区国产| a级一级毛片免费在线观看| 精品熟女少妇av免费看| 免费av毛片视频| 91久久精品国产一区二区成人| 蜜桃亚洲精品一区二区三区| 免费大片黄手机在线观看| 黄片无遮挡物在线观看| 日韩强制内射视频| 国产视频内射| 欧美+日韩+精品| 国产亚洲5aaaaa淫片| 久久精品综合一区二区三区| 最近的中文字幕免费完整| 亚洲丝袜综合中文字幕| 国产女主播在线喷水免费视频网站 | 欧美xxxx性猛交bbbb| 成人一区二区视频在线观看| 久久久亚洲精品成人影院| 国产成人精品婷婷| 国产亚洲av嫩草精品影院| av黄色大香蕉| 欧美三级亚洲精品| 亚洲熟妇中文字幕五十中出| 亚洲国产精品成人久久小说| 真实男女啪啪啪动态图| 国产成人精品久久久久久| 亚洲一级一片aⅴ在线观看| 一区二区三区免费毛片| 极品少妇高潮喷水抽搐| 三级男女做爰猛烈吃奶摸视频| 久久久久精品久久久久真实原创| 女人十人毛片免费观看3o分钟| 久热久热在线精品观看| 少妇高潮的动态图| 美女主播在线视频| 毛片一级片免费看久久久久| 日日啪夜夜爽| 国产午夜福利久久久久久| 97超视频在线观看视频| 国产淫片久久久久久久久| 成人综合一区亚洲| 日韩一本色道免费dvd| 一级a做视频免费观看| 亚洲欧美中文字幕日韩二区| 看黄色毛片网站| 人妻夜夜爽99麻豆av| 亚洲av免费在线观看| 国内精品宾馆在线| a级一级毛片免费在线观看| 男女边摸边吃奶| 亚洲国产欧美人成| 三级国产精品片| 赤兔流量卡办理| 婷婷色麻豆天堂久久| 青春草亚洲视频在线观看| 99热全是精品| 色5月婷婷丁香| 亚洲自拍偷在线| 少妇裸体淫交视频免费看高清| 欧美丝袜亚洲另类| 99re6热这里在线精品视频| 少妇人妻一区二区三区视频| 卡戴珊不雅视频在线播放| 亚洲久久久久久中文字幕| 国产精品一区二区在线观看99 | 日本黄色片子视频| 少妇高潮的动态图| 日韩三级伦理在线观看| 亚洲欧美一区二区三区国产| 一个人看视频在线观看www免费| 久久国产乱子免费精品| 能在线免费看毛片的网站| av在线老鸭窝| 国产永久视频网站| 国产国拍精品亚洲av在线观看| 在线免费十八禁| 久久久精品免费免费高清| 日韩精品青青久久久久久| 嫩草影院精品99| 免费看av在线观看网站| 成人高潮视频无遮挡免费网站| 少妇被粗大猛烈的视频| 性色avwww在线观看| 国产成人a区在线观看| 国产女主播在线喷水免费视频网站 | 别揉我奶头 嗯啊视频| 99视频精品全部免费 在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩一区二区三区影片| 大香蕉久久网| 久久99热这里只有精品18| 亚洲综合精品二区| 赤兔流量卡办理| 69av精品久久久久久| 精品99又大又爽又粗少妇毛片| h日本视频在线播放| 色吧在线观看| 亚洲成人久久爱视频| 亚洲精品一二三| 街头女战士在线观看网站| 最近中文字幕高清免费大全6| 国产精品一二三区在线看| 国产探花极品一区二区| 免费av毛片视频| 在线 av 中文字幕| 特级一级黄色大片| 美女国产视频在线观看| 国产成人aa在线观看| 国产伦理片在线播放av一区| 亚洲国产av新网站| 国内精品宾馆在线| 国产伦在线观看视频一区| 国产一区二区在线观看日韩| 一区二区三区免费毛片| 99久久精品一区二区三区| 自拍偷自拍亚洲精品老妇| 啦啦啦中文免费视频观看日本| 国产乱人偷精品视频| 大话2 男鬼变身卡| 日本av手机在线免费观看| 亚洲av在线观看美女高潮| 亚洲第一区二区三区不卡| 一级a做视频免费观看| 日本黄色片子视频| 国产老妇女一区| 高清av免费在线| 亚洲精品第二区| 十八禁网站网址无遮挡 | 大又大粗又爽又黄少妇毛片口| 男人舔奶头视频| 欧美性感艳星| 国产亚洲av片在线观看秒播厂 | 精品久久久精品久久久| 人人妻人人看人人澡| 日本色播在线视频| 男人和女人高潮做爰伦理| 赤兔流量卡办理| 3wmmmm亚洲av在线观看| 男人爽女人下面视频在线观看| kizo精华| 男女下面进入的视频免费午夜| 最近手机中文字幕大全| 91狼人影院| 99久久精品热视频| 少妇人妻精品综合一区二区| 国产黄片美女视频| 亚洲国产精品专区欧美| 欧美激情在线99| videos熟女内射| 精品酒店卫生间| 久久国产乱子免费精品| 男人和女人高潮做爰伦理| 亚洲国产日韩欧美精品在线观看| 国产成人免费观看mmmm| 狂野欧美白嫩少妇大欣赏| 少妇熟女aⅴ在线视频| 国产成人91sexporn| 国产精品久久久久久精品电影小说 | a级毛色黄片| 国产精品一区二区在线观看99 | 在线观看人妻少妇| 日韩欧美精品v在线| 亚洲国产最新在线播放| 国产精品不卡视频一区二区| 爱豆传媒免费全集在线观看| 免费看日本二区| 精品熟女少妇av免费看| 两个人视频免费观看高清| 一级爰片在线观看| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 99热这里只有是精品在线观看| 免费av观看视频| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 熟妇人妻久久中文字幕3abv| 一级片'在线观看视频| 超碰97精品在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品久久久久久精品电影| 久久久久精品性色| 伦理电影大哥的女人| 午夜激情久久久久久久| 成年av动漫网址| 国产永久视频网站| 六月丁香七月| 国产伦一二天堂av在线观看| 国产亚洲91精品色在线| 亚洲图色成人| 亚洲精品乱码久久久久久按摩| 女人十人毛片免费观看3o分钟| 大又大粗又爽又黄少妇毛片口| 国产 一区精品| 国产在视频线精品| 国产毛片a区久久久久| 国产成人精品婷婷| 伊人久久国产一区二区| 亚洲av中文av极速乱| 国产成人a∨麻豆精品| a级一级毛片免费在线观看| 又大又黄又爽视频免费| 建设人人有责人人尽责人人享有的 | 精品久久久久久久末码| 尤物成人国产欧美一区二区三区| 亚洲成人av在线免费| av国产免费在线观看| 成人特级av手机在线观看| 最近2019中文字幕mv第一页| 欧美变态另类bdsm刘玥| 久久久午夜欧美精品| 少妇裸体淫交视频免费看高清| 18禁动态无遮挡网站| 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 久久久久久久国产电影| 九草在线视频观看| 亚洲综合色惰| 国产亚洲精品久久久com| 在线观看人妻少妇| 1000部很黄的大片| 精品久久国产蜜桃| 男人爽女人下面视频在线观看| 真实男女啪啪啪动态图| 国产av码专区亚洲av| 国产成人免费观看mmmm| 禁无遮挡网站| 啦啦啦韩国在线观看视频| 午夜福利高清视频| 亚洲精品国产成人久久av| 天堂俺去俺来也www色官网 | av又黄又爽大尺度在线免费看| 春色校园在线视频观看| 又爽又黄a免费视频| 寂寞人妻少妇视频99o| av免费在线看不卡| 久久精品夜色国产| 亚洲欧美一区二区三区黑人 | 女人久久www免费人成看片| 一级毛片电影观看| 成人综合一区亚洲| 97超视频在线观看视频| 亚洲av在线观看美女高潮| 午夜福利在线在线| 国产美女午夜福利| 一级毛片黄色毛片免费观看视频| 大片免费播放器 马上看| 欧美高清成人免费视频www| 丰满人妻一区二区三区视频av| 女人被狂操c到高潮| 国产一级毛片七仙女欲春2| 少妇猛男粗大的猛烈进出视频 | 亚洲av中文av极速乱| 春色校园在线视频观看| 免费大片黄手机在线观看| 色综合站精品国产| 欧美一级a爱片免费观看看| 国产一区二区在线观看日韩| 欧美3d第一页| 国产精品国产三级专区第一集| 日韩亚洲欧美综合| 日本-黄色视频高清免费观看| 日本免费a在线| 狠狠精品人妻久久久久久综合| 国产成人freesex在线| av在线观看视频网站免费| 亚洲精品久久午夜乱码| 中文字幕免费在线视频6| 国产高清不卡午夜福利| 国产黄片视频在线免费观看| 伊人久久国产一区二区| av在线蜜桃| av在线天堂中文字幕| 黄色欧美视频在线观看| 最新中文字幕久久久久| 欧美成人精品欧美一级黄| 亚洲激情五月婷婷啪啪| 国产精品一区www在线观看| 日本免费在线观看一区| videossex国产| 十八禁网站网址无遮挡 | 国内精品一区二区在线观看| 国产综合精华液| 国产午夜精品一二区理论片| 国产精品一区二区三区四区久久| 久久久久久久亚洲中文字幕| 免费播放大片免费观看视频在线观看| 男人爽女人下面视频在线观看| 少妇熟女aⅴ在线视频| xxx大片免费视频| 日韩 亚洲 欧美在线| 亚洲,欧美,日韩| 国产 亚洲一区二区三区 | 日本黄色片子视频| 人人妻人人澡欧美一区二区| 久久人人爽人人片av| 国产女主播在线喷水免费视频网站 | 亚洲国产日韩欧美精品在线观看| 久久这里只有精品中国| 日本免费a在线| 久久久久久久久久黄片| 国产亚洲精品久久久com| 国产在线男女| 精品久久久久久电影网| 人妻少妇偷人精品九色| 九草在线视频观看| 久久久久久国产a免费观看| 国产精品伦人一区二区| 日本免费在线观看一区| 久久精品久久精品一区二区三区| 又黄又爽又刺激的免费视频.| 日日啪夜夜撸| 亚洲av免费高清在线观看| 久久久久久久久久久免费av| 777米奇影视久久| 国内精品宾馆在线| 日日摸夜夜添夜夜添av毛片| 黑人高潮一二区| 亚洲国产成人一精品久久久| 成年女人看的毛片在线观看| 丰满乱子伦码专区| 精品亚洲乱码少妇综合久久| 亚洲av中文av极速乱| 在线观看免费高清a一片| 91精品一卡2卡3卡4卡| 菩萨蛮人人尽说江南好唐韦庄| 久久精品夜色国产| 永久网站在线| 在线观看美女被高潮喷水网站| 99久久人妻综合| 国产极品天堂在线| 日韩亚洲欧美综合| 久久久久久久久久人人人人人人| or卡值多少钱| 99久久人妻综合| 久久久亚洲精品成人影院| 午夜免费激情av| 黄色一级大片看看| 欧美一区二区亚洲| 啦啦啦韩国在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 少妇熟女欧美另类| 男的添女的下面高潮视频| 最近中文字幕高清免费大全6| 80岁老熟妇乱子伦牲交| 国产成人精品婷婷| 免费看av在线观看网站| 亚洲精品自拍成人| 精品久久久精品久久久| 91狼人影院| 久久久久久久午夜电影| 99热这里只有是精品50| 看非洲黑人一级黄片| 免费黄色在线免费观看| 肉色欧美久久久久久久蜜桃 | 两个人视频免费观看高清| 亚洲国产日韩欧美精品在线观看| 一级a做视频免费观看| 99热这里只有是精品50| 国产黄频视频在线观看| 午夜激情福利司机影院| 国产一级毛片七仙女欲春2| 免费观看性生交大片5| 晚上一个人看的免费电影| 免费高清在线观看视频在线观看| 熟女人妻精品中文字幕| 少妇高潮的动态图| 嫩草影院精品99| 天堂影院成人在线观看| 亚洲精品aⅴ在线观看| 最近中文字幕高清免费大全6| 久久久精品免费免费高清| 国产精品人妻久久久影院| 高清在线视频一区二区三区| 搡女人真爽免费视频火全软件| 日韩一区二区三区影片| 在线a可以看的网站| 亚洲av电影在线观看一区二区三区 | 99九九线精品视频在线观看视频| 直男gayav资源| 国产精品国产三级专区第一集| 欧美日韩视频高清一区二区三区二| 99久久九九国产精品国产免费| 国产亚洲午夜精品一区二区久久 | 日日啪夜夜爽| 日本与韩国留学比较| 亚洲最大成人手机在线| 日韩国内少妇激情av| 床上黄色一级片| 精品人妻一区二区三区麻豆| 国产成人91sexporn| 男女那种视频在线观看| 国产黄a三级三级三级人| 久久久精品免费免费高清| 久久久精品欧美日韩精品| 99re6热这里在线精品视频| 欧美激情在线99| 国产黄片美女视频| 亚洲最大成人手机在线| 日韩精品青青久久久久久| 美女高潮的动态| 欧美bdsm另类| 国产色爽女视频免费观看| 日韩av免费高清视频| 国产在线男女| 免费观看性生交大片5| 成人亚洲精品一区在线观看 | 亚洲av福利一区| 久久这里有精品视频免费| or卡值多少钱| 精品酒店卫生间| 一区二区三区乱码不卡18| 亚洲欧美一区二区三区黑人 | 亚洲伊人久久精品综合| 成人美女网站在线观看视频| 淫秽高清视频在线观看| 伊人久久精品亚洲午夜| 一区二区三区乱码不卡18| av在线蜜桃| 一级黄片播放器| 国产真实伦视频高清在线观看| 18禁动态无遮挡网站| 亚洲综合色惰| 亚洲人成网站在线播| 天天躁夜夜躁狠狠久久av| 亚洲国产最新在线播放| 麻豆成人av视频| www.色视频.com| 亚洲一区高清亚洲精品| 成年版毛片免费区| 一级片'在线观看视频| 国产大屁股一区二区在线视频| 日本午夜av视频| 免费在线观看成人毛片| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 免费看不卡的av| 麻豆国产97在线/欧美| 可以在线观看毛片的网站| 日韩三级伦理在线观看| 超碰97精品在线观看| 久久热精品热| 久久久久精品性色| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产高清在线一区二区三| 18禁在线无遮挡免费观看视频| 精品少妇黑人巨大在线播放| 国产视频首页在线观看| 可以在线观看毛片的网站| 亚州av有码| 亚洲天堂国产精品一区在线| 男女下面进入的视频免费午夜| 亚洲精品成人av观看孕妇| 中文天堂在线官网| 国产亚洲精品久久久com| 大陆偷拍与自拍| 女人十人毛片免费观看3o分钟| 欧美人与善性xxx| 美女xxoo啪啪120秒动态图| 18+在线观看网站| 一级毛片电影观看| 99久久精品一区二区三区| 天堂√8在线中文| 成人美女网站在线观看视频| 日韩欧美三级三区| av黄色大香蕉| www.色视频.com| 99久久人妻综合| 久久这里只有精品中国| 欧美日韩亚洲高清精品| 深夜a级毛片| 欧美丝袜亚洲另类| av播播在线观看一区| 99热这里只有精品一区| 床上黄色一级片| 伊人久久精品亚洲午夜| 2021少妇久久久久久久久久久| 日产精品乱码卡一卡2卡三| 亚洲精品乱码久久久v下载方式| 国产av不卡久久| 免费看不卡的av| 69av精品久久久久久| 亚洲av日韩在线播放| 干丝袜人妻中文字幕| 最近最新中文字幕免费大全7| 久久久欧美国产精品| 女人十人毛片免费观看3o分钟| 国产免费福利视频在线观看| 男人狂女人下面高潮的视频| 亚洲av男天堂| 日韩在线高清观看一区二区三区| videossex国产| 男人舔奶头视频| 亚洲不卡免费看| 2018国产大陆天天弄谢| 麻豆成人av视频| 能在线免费看毛片的网站| 亚洲精品乱码久久久久久按摩| 99热这里只有精品一区| 男女边吃奶边做爰视频| 女人久久www免费人成看片| 男女视频在线观看网站免费| 黄色日韩在线| 欧美xxxx黑人xx丫x性爽| 亚洲经典国产精华液单| 亚洲国产精品成人久久小说| 精品久久久久久久久亚洲| 晚上一个人看的免费电影| 深夜a级毛片| 欧美极品一区二区三区四区| 国内揄拍国产精品人妻在线| 欧美高清性xxxxhd video| 国产精品久久久久久精品电影| 亚洲精品一区蜜桃| 欧美另类一区| 国产精品1区2区在线观看.| 国产精品美女特级片免费视频播放器| 亚洲欧洲国产日韩| 天堂av国产一区二区熟女人妻| 国产免费视频播放在线视频 | 天堂俺去俺来也www色官网 | 九九爱精品视频在线观看| 精品久久国产蜜桃| 网址你懂的国产日韩在线| 国产成人aa在线观看| 日韩精品有码人妻一区| 精品久久久久久久人妻蜜臀av| 国产成人a∨麻豆精品| 久久精品人妻少妇| 狠狠精品人妻久久久久久综合| 国产黄a三级三级三级人| 女人久久www免费人成看片| 久99久视频精品免费| 国内揄拍国产精品人妻在线| 国产在线男女| 一级毛片我不卡| 日本与韩国留学比较| 国产成人a区在线观看| 日本av手机在线免费观看| 我的老师免费观看完整版| 欧美 日韩 精品 国产| 国产黄a三级三级三级人| 大香蕉久久网| 少妇人妻精品综合一区二区| 久久久久久伊人网av| 麻豆成人av视频| 午夜精品在线福利| 日韩中字成人| a级一级毛片免费在线观看| 亚洲婷婷狠狠爱综合网| 日日啪夜夜撸| 久久久久久久亚洲中文字幕| 啦啦啦中文免费视频观看日本| 国产精品综合久久久久久久免费| 亚洲美女搞黄在线观看| or卡值多少钱| 免费观看在线日韩| 午夜精品国产一区二区电影 | 丝袜美腿在线中文| 91午夜精品亚洲一区二区三区| 亚洲国产精品国产精品| 久久综合国产亚洲精品| 亚洲国产高清在线一区二区三| 极品教师在线视频| 一级a做视频免费观看| 三级男女做爰猛烈吃奶摸视频| 国产乱来视频区| 国产精品国产三级国产av玫瑰| 特级一级黄色大片| 人妻少妇偷人精品九色| 九九久久精品国产亚洲av麻豆| 91狼人影院| 国产老妇伦熟女老妇高清| 亚洲av成人精品一二三区| 熟妇人妻久久中文字幕3abv| 五月天丁香电影| 26uuu在线亚洲综合色| 亚洲人成网站在线观看播放| 国模一区二区三区四区视频| 亚洲成人精品中文字幕电影| 亚洲三级黄色毛片| 国产精品久久视频播放| 国产精品无大码| 成年女人在线观看亚洲视频 | 熟妇人妻不卡中文字幕| 淫秽高清视频在线观看| 精品国产一区二区三区久久久樱花 | 久久久久久伊人网av| 国产精品99久久久久久久久| av女优亚洲男人天堂| 草草在线视频免费看| 高清视频免费观看一区二区 | 亚洲熟女精品中文字幕| 99久久精品一区二区三区| 国语对白做爰xxxⅹ性视频网站| 精华霜和精华液先用哪个|