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

    花鱸微衛(wèi)星標(biāo)記分離及其多態(tài)性分析

    2011-12-25 01:02:46季相山曾勇慶楊萍萍
    Zoological Research 2011年5期
    關(guān)鍵詞:花鱸微衛(wèi)星雜合

    趙 燕, 季相山, 曾勇慶, 丁 雷, 楊萍萍, 王 慧

    (山東農(nóng)業(yè)大學(xué) 動(dòng)物科技學(xué)院, 山東 泰安 271018)

    花鱸微衛(wèi)星標(biāo)記分離及其多態(tài)性分析

    趙 燕, 季相山, 曾勇慶, 丁 雷, 楊萍萍, 王 慧*

    (山東農(nóng)業(yè)大學(xué) 動(dòng)物科技學(xué)院, 山東 泰安 271018)

    該文利用FIASCO法(fast isolation by AFLP of sequences containing repeats)和GenBank數(shù)據(jù)庫(kù)搜索法開(kāi)發(fā)花鱸微衛(wèi)星標(biāo)記, 并對(duì)篩選的標(biāo)記進(jìn)行多態(tài)性檢測(cè)。兩種方法共獲得54條能夠設(shè)計(jì)引物的序列, 擴(kuò)增結(jié)果顯示15對(duì)引物具有多態(tài)性, 多態(tài)性微衛(wèi)星位點(diǎn)的等位基因數(shù)為 2~10個(gè)。15個(gè)多態(tài)性位點(diǎn)中, 4個(gè)位點(diǎn)偏離了Hardy-Weinberg平衡; 各位點(diǎn)間沒(méi)有連鎖不平衡現(xiàn)象; 僅位點(diǎn)SP52可能存在無(wú)效等位基因; 除SP17和SP468外,其余引物的PIC值均在0.5以上, 可用于花鱸群體遺傳分析等研究。

    FIASCO法; 數(shù)據(jù)庫(kù)查詢法; 花鱸; 微衛(wèi)星

    微衛(wèi)星具有多態(tài)性豐富、易于快速檢測(cè)等特點(diǎn),在種質(zhì)資源保存利用、遺傳結(jié)構(gòu)和遺傳變異及遺傳連鎖圖譜構(gòu)建等方面的研究中發(fā)揮著非常重要的作用, 已經(jīng)發(fā)展成為當(dāng)前主流的分子標(biāo)記之一(Wang et al, 2009)。FIASCO法(fast isolation by AFLP of sequences containing repeats)是基于AFLP技術(shù)的一種非常有效的分離微衛(wèi)星位點(diǎn)的方法(Zane et al, 2002)。它將AFLP與雜交選擇相結(jié)合, 能夠簡(jiǎn)單快速提高陽(yáng)性率, 發(fā)展前景較大, 已經(jīng)在黃鱔(Monopterus albus) (Li et al, 2007)、鱖魚(yú)(Siniperca chuatsi) (Kuang et al, 2009)等多種魚(yú)類微衛(wèi)星的開(kāi)發(fā)上使用。數(shù)據(jù)庫(kù)搜索法是從公用的DNA 序列數(shù)據(jù)庫(kù)(GenBank、EMBL和DDBJ等)中查找所研究物種的序列信息, 并用于開(kāi)發(fā)微衛(wèi)星標(biāo)記的一種簡(jiǎn)單經(jīng)濟(jì)有效的方法(Brown et al, 1996)。

    當(dāng)前已開(kāi)發(fā)的花鱸微衛(wèi)星標(biāo)記較少, 僅 Jiang et al (2007, 2008)發(fā)表了22個(gè)微衛(wèi)星標(biāo)記, 其中有13個(gè)微衛(wèi)星標(biāo)記為混合微衛(wèi)星, 不適于種群遺傳學(xué)研究, 遠(yuǎn)遠(yuǎn)滿足不了花鱸遺傳育種等研究的需要。在微衛(wèi)星開(kāi)發(fā)的幾種方法中, 從基因組文庫(kù)中篩選,耗時(shí)、費(fèi)力, 效率較低; 從近緣中擴(kuò)增, 成功率較低。本研究選擇目前應(yīng)用較廣泛的FIASCO法構(gòu)建花鱸部分基因組文庫(kù)及參照 GenBank上花鱸的相關(guān)序列設(shè)計(jì)微衛(wèi)星引物兩種方法, 對(duì)花鱸微衛(wèi)星位點(diǎn)的多態(tài)性進(jìn)行研究, 這可為花鱸群體遺傳結(jié)構(gòu)的分析、遺傳資源的保護(hù)利用和分子標(biāo)記輔助育種等打下基礎(chǔ)。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料

    實(shí)驗(yàn)所用花鱸活體樣本捕自山東煙臺(tái)近海, 共30個(gè)個(gè)體, 取其尾鰭置于75%酒精中, 保存于-20 ℃冰箱中備用。

    1.2 微衛(wèi)星富集文庫(kù)的構(gòu)建與序列測(cè)定

    用蛋白酶K消化和酚/氯仿抽提法從花鱸的鰭條中分離、純化基因組DNA (Sambrook & Russell, 2002)。

    基因組DNA用限制性內(nèi)切酶MseⅠ酶切后回收250~1 000 bp的片段, 回收產(chǎn)物與接頭oligo-MseⅠA (HO 5'-TACTCAGGACTCAT-3' OH)和oligo-MseⅠB (HO 3'-GAGTCCTGAGTAGCAG-5' OH) 連接。建立17.5 μL連接體系(接頭20 μmol/L, 2.5 μL; T4 DNA ligase 6 U, 10×buffer 1.75 μL, 純化的酶切產(chǎn)物 10 μL), 置16 ℃水浴中過(guò)夜連接。利用預(yù)擴(kuò)引物MseⅠ- N(5'-ATGAGTCCTGAGTAAN-3')對(duì)連接了接頭的DNA片段進(jìn)行PCR擴(kuò)增, 程序?yàn)?4 ℃預(yù)變性5 min, 94 ℃變性30 s, 53 ℃退火30 s, 72 ℃延伸 1 min, 共進(jìn)行 30個(gè)循環(huán), 最后 72 ℃延伸10 min?;厥?50~1 000 bp的預(yù)擴(kuò)增片段用于雜交分析。

    雜交體系包含生物素標(biāo)記的(GT)13的探針 6.0 μL(40 μmol/L), 雜交液(6×SSC + 0.1% SDS)70 μL,預(yù)擴(kuò)增回收產(chǎn)物 20 μL。磁珠經(jīng) 250 μL的TEN100(10 mmol/L Tris-HCl, 1 mmol/L EDTA, 100 mmol/L NaCl, pH 7.5)洗滌3次, 每次5 min后, 將雜交產(chǎn)物和300 μL TEN100加入到預(yù)處理過(guò)的磁珠中, 室溫反應(yīng)30 min。分別用400 μL TEN1000(0.2 ×SSC + 0.1% SDS)和洗脫液(0.2×SSC + 0.1% SDS)各洗滌磁珠3次, 每次5 min, 其間不時(shí)用移液器輕輕吹打, 重懸磁珠。最后用400 μL 1×SSC洗滌磁珠5 min。加50 μL TE, 95 ℃水浴中變性5 min, 洗脫含有微衛(wèi)星序列的單鏈 DNA。以洗脫的含有微衛(wèi)星的單鏈 DNA為模板再次用引物MseI-N 進(jìn)行擴(kuò)增, 擴(kuò)增程序同上, 回收 250~1 000 bp的片段。

    將回收的目的片段連接到pMD18-T載體中并轉(zhuǎn)化至大腸桿菌 DH5α中, 涂平板后挑取單克隆,用引物MseI-N進(jìn)行PCR檢測(cè), 挑選陽(yáng)性克隆送上海生物工程技術(shù)服務(wù)有限責(zé)任公司測(cè)序。利用DNAMAN Version 4.0 和GenBank中的Blastn軟件對(duì)測(cè)得序列進(jìn)行比對(duì), 去掉冗余序列。

    1.3 GenBank數(shù)據(jù)庫(kù)搜索

    輸入“Laterolabrax japonicus”作為搜索關(guān)鍵詞,在 NCBI公共數(shù)據(jù)庫(kù)(http: // www.ncbi.nlm.nih.gov)搜索核苷酸序列, 下載所有相關(guān)序列以 FSATA格式保存?zhèn)溆谩?/p>

    1.4 序列分析及微衛(wèi)星引物設(shè)計(jì)

    利用 SSRHUNTER 1.3(http://en.bio-soft.net/ dna/SSRHunter.html)軟件, 對(duì)以上兩種方法得到的序列進(jìn)行分析, 含二堿基、三堿基、四堿基重復(fù) 5次以上的為含微衛(wèi)星的序列。根據(jù)微衛(wèi)星的側(cè)翼序列, 利用軟件Premier 5.0設(shè)計(jì)引物, 并交由上海生物工程技術(shù)服務(wù)有限責(zé)任公司合成。

    1.5 微衛(wèi)星標(biāo)記的多態(tài)性檢測(cè)

    利用合成的引物分別在 30個(gè)花鱸個(gè)體中進(jìn)行PCR擴(kuò)增用于分析其多態(tài)性。PCR反應(yīng)總體積為20μL , 內(nèi)含20 ng的模板DNA, 0.2 μmol/L的引物, 200 μmol/L的dNTPs, 1.5 mmol/L 的Mg2+, 1 X PCR 反應(yīng)緩沖液, 1 U 的Taq DNA 聚合酶。PCR反應(yīng)為 30個(gè)循環(huán), 每個(gè)循環(huán)包括: 94℃變性 30 s,退火30 s, 72 ℃延伸1 min; 首次循環(huán)前預(yù)變性5 min; 最后一次循環(huán)結(jié)束后 72 ℃再延伸 10 min。PCR 擴(kuò)增產(chǎn)物經(jīng) 6%變性聚丙烯酰胺凝膠電泳, 硝酸銀染色后數(shù)碼相機(jī)拍照保存。

    1.6 數(shù)據(jù)分析

    微衛(wèi)星標(biāo)記呈共顯性, 可直接根據(jù)每個(gè)個(gè)體產(chǎn)生的條帶的位置確定基因型, 利用Popgene32 處理數(shù)據(jù), 計(jì)算多態(tài)性位點(diǎn)的等位基因數(shù)目、觀測(cè)雜合度(HO)、期望雜合度(HE), 并按下列公式計(jì)算多態(tài)性信息含量(PIC):

    其中,Pi和Pj分別為第i和第j個(gè)等位基因的頻率,n為等位基因個(gè)數(shù)。

    利用 MICRO-CHECKER Version 2.2.3軟件估計(jì)位點(diǎn)的無(wú)效等位基因頻率。所有多重比較均經(jīng)Bonferroni校正。

    2 結(jié) 果

    2.1 FIASCO法篩選的微衛(wèi)星標(biāo)記

    從構(gòu)建的文庫(kù)中隨機(jī)挑取150個(gè)克隆, 利用預(yù)擴(kuò)引物進(jìn)行PCR鑒定, 共得到陽(yáng)性克隆94個(gè), 陽(yáng)性率為62.7%; 對(duì)陽(yáng)性克隆測(cè)序得到66個(gè)含微衛(wèi)星序列的克隆, 陽(yáng)性序列比例為70.2%。66條微衛(wèi)星序列中共含微衛(wèi)星座位108個(gè), 其中CA/GT重復(fù)序列97個(gè), 占89.8%, GA/CT重復(fù)4個(gè), 占3.7%, 其他重復(fù)序列7個(gè), 占6.5%。在這些含微衛(wèi)星重復(fù)的序列中, 微衛(wèi)星重復(fù)數(shù)在10次以上的, 占70.3%。在這些位點(diǎn)的重復(fù)序列中, 共3個(gè)復(fù)合型(sp52、sp416、sp468), 其他為完美型。利用 DNAMAN 4.0 和GenBank中的Blastn軟件對(duì)這66條微衛(wèi)星序列進(jìn)行比對(duì), 去掉冗余序列后全部提交 GenBank, 并利用 Premier 5.0軟件共設(shè)計(jì)引物42對(duì), 部分微衛(wèi)星序列及引物信息見(jiàn)表1。

    2.2 GenBank 數(shù)據(jù)庫(kù)搜索結(jié)果

    表1 FIASCO法獲得的部分花鱸微衛(wèi)星序列Tab. 1 Microsatellite sequences obtained by FIASCO

    截至2010年2月, 在GenBank數(shù)據(jù)庫(kù)中共搜索到540條花鱸核苷酸序列及132條EST序列, 經(jīng)SSRHUNTER 1.3軟件分析, 含微衛(wèi)星的序列為22條, 其中8條為基因的全長(zhǎng)cDNA序列; 1條為基因的部分cDNA序列; 13條為EST序列。在22條微衛(wèi)星序列中含微衛(wèi)星座位 29個(gè), 其中CA/GT重復(fù)序列14個(gè), 占48.3%,GA/CT重復(fù)10個(gè), 占34.5%,三堿基重復(fù)序列5個(gè), 占17.2%。在這些含微衛(wèi)星重復(fù)的序列中, 微衛(wèi)星重復(fù)數(shù)在10次以上的僅6個(gè),重復(fù)數(shù)在5~7次間的為22個(gè)(75.8%)。在8條基因的全長(zhǎng) cDNA序列中, 微衛(wèi)星座位多數(shù)位于內(nèi)含子、3′-UTR、5′-UTR等區(qū)域, 僅有1個(gè)位于外顯子區(qū)域(表2)。

    2.3 微衛(wèi)星標(biāo)記多態(tài)性分析

    兩種方法篩選的微衛(wèi)星序列共設(shè)計(jì)引物54對(duì),利用8尾花鱸的DNA對(duì)各引物擴(kuò)增情況作初步篩選, 發(fā)現(xiàn)32對(duì)能擴(kuò)增出與期望PCR產(chǎn)物大小相等或相近的片段, 其中 15對(duì)引物的擴(kuò)增片段有穩(wěn)定的多態(tài)性(表3), 其余的都為單態(tài)。應(yīng)用這15對(duì)引物對(duì)30個(gè)個(gè)體進(jìn)行了PCR擴(kuò)增, 擴(kuò)增條帶清晰、穩(wěn)定, 其中微衛(wèi)星引物SPd10擴(kuò)增產(chǎn)物的聚丙烯酰胺凝膠電泳圖譜見(jiàn)圖1。

    表2 GenBank數(shù)據(jù)庫(kù)查詢得到的部分微衛(wèi)星序列Tab. 2 Microsatellite sequences obtained by bioinformatic screening in GenBank

    在野生花鱸30個(gè)個(gè)體中, 15個(gè)花鱸多態(tài)性微衛(wèi)星位點(diǎn)共檢測(cè)到89個(gè)等位基因, 平均每個(gè)位點(diǎn)5.9個(gè)等位基因, 表現(xiàn)出高度多態(tài)性, 各等位基因大小在141~349之間(表3)。Popgene32 分析結(jié)果顯示, 15個(gè)多態(tài)性微衛(wèi)星位點(diǎn)的觀測(cè)雜合度在0.6000~序的微衛(wèi)星序列中, CA/GT重復(fù)序列97個(gè), 占1.0000之間, 期望雜合度在 0.5079~0.8890之間。經(jīng)Bonferroni 校正后(P< 0.003), 各位點(diǎn)間沒(méi)有連鎖不平衡現(xiàn)象, 4個(gè)微衛(wèi)星位點(diǎn)(SP17、SP52、 SP94 和SP468)偏離了 Hardy-Weinberg平衡。MICROCHECKER 2.2.3軟件分析表明僅位點(diǎn)SP52可能存在無(wú)效等位基因。同時(shí), 對(duì) 15對(duì)微衛(wèi)星引物的多態(tài)性信息含量(PIC)進(jìn)行了評(píng)估, SP17位點(diǎn)提供的信息含量為 0.4189; SP468為 0.3747, 其余引物的PIC值均在0.5以上, 可用于花鱸群體遺傳分析等研究。

    圖1 微衛(wèi)星標(biāo)記SPd10在花鱸30個(gè)個(gè)體中的擴(kuò)增圖譜Fig. 1 Amplification results of microsatellite locus SPd10 in 30 Lateolabrax japonicus individuals

    表3 15個(gè)多態(tài)性花鱸微衛(wèi)星標(biāo)記分析結(jié)果Tab. 3 Characteristics of 15 polymorphic microsatellite markers in Lateolabrax japonicus

    3 討 論

    目前開(kāi)發(fā)微衛(wèi)星標(biāo)記的方法很多, 概括起來(lái)有3種: (1)從基因組文庫(kù)中篩選微衛(wèi)星; (2)從數(shù)據(jù)庫(kù)或有關(guān)文章中查詢; (3)從近緣物種中篩選微衛(wèi)星(Zhan et al, 2008)。本文研究發(fā)現(xiàn), 數(shù)據(jù)庫(kù)搜索法得到的花鱸微衛(wèi)星位點(diǎn)的多態(tài)率為 41.7%, 多態(tài)性較高。其他的研究也得到類似的結(jié)果, 如Ji et al(2009)搜索 GenBank數(shù)據(jù)庫(kù)中的牙鲆(Paralichthys olivaceus)全cDNA序列, 得到15個(gè)具有高多態(tài)性的微衛(wèi)星位點(diǎn), 多態(tài)率為 46.8%; Furukawa et al (2004)用同樣的方法搜索到11個(gè)具有多態(tài)性的河豚微衛(wèi)星位點(diǎn), 多態(tài)率為73%。從數(shù)據(jù)庫(kù)中查找微衛(wèi)星位點(diǎn), 雖然省時(shí)省力, 但受限于數(shù)據(jù)庫(kù)中提交序列的數(shù)量, 多數(shù)物種因提交的序列較少只能開(kāi)發(fā)出有限的多態(tài)性微衛(wèi)星標(biāo)記。FIASCO法是Zane et al (2002)年提出的一種高效的篩選微衛(wèi)星標(biāo)記的方法。此法一經(jīng)提出, 便被迅速?gòu)V泛地應(yīng)用于許多物種的微衛(wèi)星位點(diǎn)的篩選(Reed et al, 2001)。本文利用FIASCO法構(gòu)建了花鱸(GT)13的微衛(wèi)星富集文庫(kù),在所有測(cè)序的微衛(wèi)星序列中,CA/GT重復(fù)序列 97個(gè),占89.8%。而通過(guò)傳統(tǒng)的篩選小插入片段基因文庫(kù)的方法分離微衛(wèi)星標(biāo)記, 得到的陽(yáng)性克隆比例極低(0.1%~6.4%, 平均1.96%) (Zane et al., 2002)。

    雖然微衛(wèi)星可分布于基因組的任何地方, 但是在真核生物中, 微衛(wèi)星位于非編碼區(qū)的頻率遠(yuǎn)高于位于編碼區(qū)的(Toth et al, 2000)。本文研究發(fā)現(xiàn), 在花鱸8條基因的全長(zhǎng)cDNA序列中, 微衛(wèi)星座位多數(shù)位于內(nèi)含子、3'-UTR、5'-UTR等區(qū)域, 僅有1個(gè)位于外顯子區(qū)域, 占 12.5%, 這與其它脊椎動(dòng)物中微衛(wèi)星在外顯子區(qū)的分布頻率(9%~15%)(Moran et al, 1993; Jurka et al, 1995)基本一致。這種微衛(wèi)星位點(diǎn)在編碼區(qū)的低頻現(xiàn)象可能是為了防止蛋白質(zhì)翻譯時(shí)發(fā)生移碼突變的一種保護(hù)機(jī)制(Metzgar et al, 2000; Li et al, 2007)。

    本研究得到的15個(gè)多態(tài)性微衛(wèi)星位點(diǎn)中有4個(gè)(SP17、SP52、SP94和SP468)偏離了Hardy-Weinberg平衡。偏離HWE的原因有很多, 例如, 無(wú)效等位基因、檢測(cè)個(gè)體不具有代表性、交配群體等。軟件分析表明, 僅位點(diǎn)SP52可能存在無(wú)效等位基因, 造成其他3個(gè)位點(diǎn)偏離HWE的原因可能是檢測(cè)個(gè)體數(shù)目過(guò)小, 檢測(cè)的個(gè)體不具代表性等所致。

    微衛(wèi)星位點(diǎn)的多態(tài)性一般用3個(gè)指標(biāo)評(píng)價(jià): 等位基因的數(shù)目、雜合度和多態(tài)性信息含量值(PIC)。其中等位基因數(shù)目和雜合度能較好地反映群體中等位基因的豐富程度和均勻程度, 而多態(tài)性信息含量值反映了微衛(wèi)星位點(diǎn)能夠提供的遺傳信息容量。本研究得到的花鱸 15個(gè)多態(tài)性微衛(wèi)星位點(diǎn)共檢測(cè)到 89個(gè)等位基因, 平均每個(gè)位點(diǎn) 5.9個(gè)等位基因,表現(xiàn)出高度多態(tài)性, 且 15個(gè)多態(tài)性微衛(wèi)星位點(diǎn)的觀測(cè)雜合度在 0.6000~1.0000之間, 期望雜合度在0.5079~0.8890之間; 而Jiang et al(2007)對(duì)22個(gè)微衛(wèi)星位點(diǎn)的多態(tài)性分析顯示, 其觀測(cè)雜合度在0.114~0.914之間, 期望雜合度在0.111~0.938之間,這表明花鱸野生群體具有較高的遺傳多樣性。據(jù)Gu et al (2008), 當(dāng)PIC>0.5 時(shí), 表明該遺傳標(biāo)記可提供豐富的遺傳信息。本文研究結(jié)果表明, 15個(gè)花鱸多態(tài)性微衛(wèi)星位點(diǎn)中的13個(gè)位點(diǎn)PIC值均大于0.5, 說(shuō)明本文新開(kāi)發(fā)的微衛(wèi)星位點(diǎn)提供的遺傳信息容量大,能為群體遺傳性分析提供有力的工具。

    本文利用 FIASCO法和數(shù)據(jù)庫(kù)查詢法兩種方法開(kāi)發(fā)了花鱸微衛(wèi)星位點(diǎn) 15個(gè), 并對(duì)其多態(tài)性進(jìn)行了研究, 結(jié)果表明數(shù)據(jù)庫(kù)查詢法簡(jiǎn)便快捷, 多態(tài)性好; FIASCO法富集效率高, 適合用于花鱸微衛(wèi)星標(biāo)記大規(guī)模的開(kāi)發(fā)。開(kāi)發(fā)的多態(tài)性微衛(wèi)星標(biāo)記可為花鱸群體遺傳學(xué)研究、連鎖圖譜構(gòu)建等提供有力的支持。

    Brown SM, Hopkins MS, Mitchell SE. 1996. Multiple methods for the identification of polymorphic simple sequence repeats (SSRs) in sorghum [sorghum bicolor(L.) Moench] [J].Theol Appl Genet, 93: 190-198.

    Furukawa S, Takeshima H, Otaka T, Misuboshi T, Shirasu K, Ikeda D, Kaneko G, Nishida M, Watabe S. 2004. Isolation of microsatellite markers by in silico screening implicated for genetic linkage mapping in Japanese pufferfishTakifugu rubripe s[J].Fish Sci, 70: 620-628.

    Gu XY, Zeng QG, You ZJ, Lin ZH.2008. Isolation and characteristics of six microsatellite primers oftegillarca granosa[J]. Oceanol Limnol Sin, 39: 661-664. [顧曉英, 曾慶國(guó), 尤仲杰. 2008. 泥蚶(Tegillarca granosa)6個(gè)微衛(wèi)星引物的分離和鑒定. 海洋與湖沼. 39: 661-664.]

    Ji XS, Chen SL, Ma HY, Xu TJ, Liao XL, Jiang YL. 2009. Isolation and characterization of 19 EST-linked polymorphic microsatellite loci for olive flounder (Paralichthys olivaceus) [J].Aquac Res, 40: 980-983.

    Jiang X, Liao MJ, Liu YJ, Gao T, Yang GP. 2007. Isolation and characterization of 22 polymorphic microsatellite DNA markers of Japanese sea bass (Laterolabrax japonicus). [J].Mol Ecol Notes, 7: 492-494.

    Jiang X, Yang G, Liao M, Liu Y, Gao T, Wang D, Wei Q, Zou G. 2008. Microsatellite DNA polymorphism of Japanese sea bass (Laterolabrax japonicus) inhabiting Chinese and Japanese coasts[J].J Appl Ichthy, 24: 180-186.

    Jurka J, Pethiyagoda C. 1995. Simple repetitive DNA sequences from primates: compilation and analysis[J].J Mol Evol, 40: 120-126.

    Kuang GQ, Lu SQ, Zheng SM, Wu Q. 2009. Isolation and evaluation of 18 microsatellite markers inSiniperca chuatsi(Basilewsky) [J].Mol Ecol Resour, 9: 1473-1475.

    Li WT, Liao XL, Yu XM, Cheng L, Tong J. 2007. Isolation and characterization of polymorphic microsatellites in a sex-reversal fish, rice field eel (Monopterus albus) [J].Mol Ecol Notes, 7: 705-707.

    Metzgar D, Wills C. 2000. Evolutionary changes in mutation rates and spectra and their influence on the adaptation of pathogens[J].Microbes Infect, 2: 1513-1522.

    Moran C. 1993. Microsatellite repeats in pig (Sus domestica) and chicken (Gallus domesticus) genomes[J].J Hered, 84: 274-280.

    Reed K M, Beanie C W. 2001. Isolation of 105 microsatellite loci from an ovine genomic library enriched for microsatellites[J].Anim Biotech, 12: 77- 86.

    Toth G, Gaspari Z, Jurka J. 2000. Microsatellites in different eukaryotic genomes: survey and analysis[J].Genome Res, 10: 967-981.

    Wang JP, Yu XW, Shen PY, 2009. Microsatellite loci of portunus trituberculatus: isolation and characterization[J].Oceanol Limnol Sin, 40: 83-87. [王建平, 余曉巍, 沈龐幼, 2009. 三疣梭子蟹(Portunus trituberculatus)微衛(wèi)星位點(diǎn)的分離和序列分析。海洋與湖沼, 40: 83-87.]

    Zane L, Bargelloni L, Patarnello T. 2002. Strategies for microsatellite isolation: a review[J].Mol Ecol, 11: 1-16.

    Zhan AB, Bao ZM, Hu XL, Wang S, Peng W, Wang ML, Hu JJ, Liang CZ, Yue ZQ. 2008. Characterization of 95 novel microsatellite markers for Zhikong scallopChlamys farreriusing FIASCO-colony hybridization and EST database mining[J].Fish Sci, 74: 516-526.

    Zhang YW, Zhang YP, Oliver AR. 2001. Microsatellites and its application[J].Zool Res,22(4): 315-320.[張?jiān)莆? 張亞平. 2001. 微衛(wèi)星及其應(yīng)用. 動(dòng)物學(xué)研究, 22(4): 315-320.]

    Isolation of microsatellite markers forLateolabrax japonicusand polymorphic analysis

    ZHAO Yan, JI Xiang-Shan, ZENG Yong-Qing, DING Lei, YANG Ping-Ping, WANG Hui*

    (College of Animal Science and Technology,Shandong Agriculture University,Taian271018,China)

    To investigate population structure and marker assisted breeding, fast isolation by AFLP of sequences containing repeats (FIASCO) and GenBank database mining were used to develop novel microsatellite markers for sea perch (Lateolabrax japonicus). Genomic DNA fragments containing SSR sequences were captured by hybridization to (GT)13biotin-labeled probe and were ligated to PMD18-T vector. Among 150 randomly chosen clones from the SSR-enriched library, 66 sequences contained microsatellite motif over five repeats. In addition, 540 cDNA sequences and 132 ESTs ofLateolabrax japonicuswere downloaded from GenBank and screened for di-, tri- and tetra-nucleotide repeats, while 22 sequences were found to contain microsatellites. As a result, 15 microsatellite loci were shown to be polymorphic in 30Lateolabrax japonicusindividuals, with the alleles ranging from two to ten, the observed heterozygosities from 0.6000 - 1.0000, and the expected heterozygosities from 0.5079 - 0.8890. Four loci (SP17, SP52, SP94 and SP468) were deviated from HWE in the sampled population after Bonferroni’s correction, and no linkage disequilibrium was found among all loci (P<0.003), whereas null alleles were detected at locus SP52 (P<0.05). Among 15 polymorphic loci, the PIC values, which can be used for related population genetics analysis, were all above 0.5, with the exception of SP17 and SP468.

    FIASCO; Database mining;Lateolabrax japonicus; Microsatellite

    Q95-33; Q959.483

    A

    0254-5853-(2011)05-0515-06

    10.3724/SP.J.1141.2011.05515

    2011-02-21; 接受日期: 2011-07-08

    國(guó)家“863”計(jì)劃項(xiàng)目(2008AA101010); 國(guó)家發(fā)改委生物育種高技術(shù)專項(xiàng)(2007249051); 山東省農(nóng)業(yè)良種工程重大項(xiàng)目(2007LZ013); 博士后創(chuàng)新基金項(xiàng)目(200703044)

    ?通訊作者(Corresponding author),王慧(1963 - ), 女, 教授, 博導(dǎo),Tel: 13001771684; E-mail: wanghui2328@sdau.edu.cn

    猜你喜歡
    花鱸微衛(wèi)星雜合
    梅花鹿基因組微衛(wèi)星分布特征研究
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見(jiàn)的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    我國(guó)花鱸養(yǎng)殖產(chǎn)業(yè)現(xiàn)狀與種子工程研究進(jìn)展
    花鱸運(yùn)輸密度對(duì)其成活率的影響
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    花鱸攝食調(diào)控相關(guān)因子全長(zhǎng)cDNA的克隆和序列分析
    基于轉(zhuǎn)錄組測(cè)序的波紋巴非蛤微衛(wèi)星標(biāo)記研究
    中國(guó)“一箭雙星”成功將“遙感衛(wèi)星二十一號(hào)”與“天拓二號(hào)視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    海產(chǎn)花鱸IGFBP-1、2基因的克隆及表達(dá)分析
    免费观看人在逋| 国语自产精品视频在线第100页| av在线蜜桃| 午夜a级毛片| 精品国产乱码久久久久久男人| 亚洲欧美日韩高清专用| 成人欧美大片| 好男人电影高清在线观看| 深夜精品福利| 欧美色视频一区免费| 日韩欧美国产在线观看| 全区人妻精品视频| 波多野结衣高清作品| 国产精品久久久久久人妻精品电影| 午夜激情福利司机影院| 中文字幕久久专区| 麻豆成人av在线观看| 在线a可以看的网站| 日韩成人在线观看一区二区三区| 免费av毛片视频| 日韩欧美免费精品| 欧美一级毛片孕妇| 免费无遮挡裸体视频| 精华霜和精华液先用哪个| 99国产综合亚洲精品| av女优亚洲男人天堂 | 国产一区二区三区视频了| 久久精品夜夜夜夜夜久久蜜豆| 精品国产乱码久久久久久男人| 97超级碰碰碰精品色视频在线观看| 看片在线看免费视频| 人妻夜夜爽99麻豆av| 在线免费观看的www视频| 国产精品久久视频播放| 欧美国产日韩亚洲一区| 精品国产亚洲在线| 美女黄网站色视频| 一级作爱视频免费观看| 在线国产一区二区在线| 97人妻精品一区二区三区麻豆| 成在线人永久免费视频| 成人三级做爰电影| 黄色日韩在线| 国产精品av视频在线免费观看| 啪啪无遮挡十八禁网站| 国产人伦9x9x在线观看| 老汉色∧v一级毛片| 五月玫瑰六月丁香| 女警被强在线播放| 欧美乱色亚洲激情| 少妇丰满av| 亚洲av熟女| 99久久综合精品五月天人人| 免费av不卡在线播放| 一区福利在线观看| 色精品久久人妻99蜜桃| 人妻丰满熟妇av一区二区三区| 成年女人看的毛片在线观看| 久久精品影院6| 1024手机看黄色片| 9191精品国产免费久久| 日本与韩国留学比较| 叶爱在线成人免费视频播放| 99久久国产精品久久久| 亚洲自偷自拍图片 自拍| 亚洲美女黄片视频| 黄色视频,在线免费观看| 精品无人区乱码1区二区| 成人精品一区二区免费| 美女大奶头视频| 日日干狠狠操夜夜爽| 精品午夜福利视频在线观看一区| 白带黄色成豆腐渣| 亚洲精品色激情综合| 久久午夜亚洲精品久久| 亚洲国产欧美网| 久久精品人妻少妇| 精品免费久久久久久久清纯| 亚洲欧洲精品一区二区精品久久久| 一二三四在线观看免费中文在| 国产精品久久久人人做人人爽| 麻豆av在线久日| av黄色大香蕉| 精品国产超薄肉色丝袜足j| 欧美黑人欧美精品刺激| 国产视频内射| 免费av不卡在线播放| 色吧在线观看| 亚洲国产色片| 熟妇人妻久久中文字幕3abv| 婷婷六月久久综合丁香| 男人和女人高潮做爰伦理| 啦啦啦观看免费观看视频高清| 日韩欧美 国产精品| 亚洲精品国产精品久久久不卡| 精品熟女少妇八av免费久了| 搞女人的毛片| 在线永久观看黄色视频| 久9热在线精品视频| 天堂av国产一区二区熟女人妻| 99国产极品粉嫩在线观看| 黄色视频,在线免费观看| 久久草成人影院| 国产精品国产高清国产av| 久久亚洲精品不卡| 好男人在线观看高清免费视频| 国产1区2区3区精品| 此物有八面人人有两片| cao死你这个sao货| 国内久久婷婷六月综合欲色啪| 国产亚洲av高清不卡| 亚洲熟妇中文字幕五十中出| 又爽又黄无遮挡网站| 特级一级黄色大片| 欧美午夜高清在线| 亚洲人与动物交配视频| 少妇熟女aⅴ在线视频| 少妇丰满av| 99riav亚洲国产免费| 天堂网av新在线| 又粗又爽又猛毛片免费看| 天天添夜夜摸| 99久久成人亚洲精品观看| 亚洲精品在线美女| 久久中文字幕人妻熟女| 亚洲第一电影网av| 精品电影一区二区在线| 欧美成人性av电影在线观看| 午夜精品一区二区三区免费看| 12—13女人毛片做爰片一| 欧美性猛交黑人性爽| 国产69精品久久久久777片 | 在线观看日韩欧美| 香蕉久久夜色| 国产精品久久久久久精品电影| 亚洲五月婷婷丁香| 中文资源天堂在线| 18禁国产床啪视频网站| 两性夫妻黄色片| 看片在线看免费视频| 亚洲国产精品999在线| cao死你这个sao货| 看片在线看免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 男插女下体视频免费在线播放| 99久久精品一区二区三区| 人人妻人人澡欧美一区二区| 亚洲五月婷婷丁香| 免费在线观看视频国产中文字幕亚洲| 97超视频在线观看视频| 亚洲色图 男人天堂 中文字幕| 精品久久久久久久久久免费视频| www.自偷自拍.com| 99久久无色码亚洲精品果冻| 99久久精品一区二区三区| 成人三级黄色视频| 给我免费播放毛片高清在线观看| 色综合婷婷激情| 欧美在线黄色| 禁无遮挡网站| 成人高潮视频无遮挡免费网站| 国产欧美日韩精品亚洲av| 精品无人区乱码1区二区| 色老头精品视频在线观看| 欧美黄色片欧美黄色片| 国产97色在线日韩免费| 国产99白浆流出| 亚洲狠狠婷婷综合久久图片| 黄色丝袜av网址大全| 国产精品综合久久久久久久免费| 免费人成视频x8x8入口观看| 色综合婷婷激情| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美日韩一区二区精品| 欧美不卡视频在线免费观看| 九九热线精品视视频播放| 亚洲av电影不卡..在线观看| 午夜福利高清视频| 嫩草影院入口| 91九色精品人成在线观看| 亚洲美女视频黄频| 很黄的视频免费| 亚洲午夜精品一区,二区,三区| 亚洲18禁久久av| av片东京热男人的天堂| 精品国产超薄肉色丝袜足j| 黄色视频,在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 成人av一区二区三区在线看| 久久亚洲真实| 黑人操中国人逼视频| 男女做爰动态图高潮gif福利片| 国产欧美日韩一区二区精品| 黄色成人免费大全| 成在线人永久免费视频| 制服人妻中文乱码| 午夜免费成人在线视频| 欧美色视频一区免费| 免费看光身美女| 亚洲av免费在线观看| e午夜精品久久久久久久| 99视频精品全部免费 在线 | 99视频精品全部免费 在线 | 国产一级毛片七仙女欲春2| 最近视频中文字幕2019在线8| 日本在线视频免费播放| 18禁观看日本| 久久精品91无色码中文字幕| 久久精品亚洲精品国产色婷小说| 天堂√8在线中文| 中出人妻视频一区二区| 人妻夜夜爽99麻豆av| 岛国在线免费视频观看| 麻豆一二三区av精品| 最新中文字幕久久久久 | 三级男女做爰猛烈吃奶摸视频| 国产精品九九99| 熟女少妇亚洲综合色aaa.| 身体一侧抽搐| 欧美成人性av电影在线观看| 99久久99久久久精品蜜桃| 午夜免费观看网址| 午夜福利成人在线免费观看| 久久久久亚洲av毛片大全| 最近最新免费中文字幕在线| 特级一级黄色大片| 久久久国产精品麻豆| svipshipincom国产片| 又粗又爽又猛毛片免费看| 精品久久久久久久末码| 一本综合久久免费| 亚洲中文av在线| 亚洲精品乱码久久久v下载方式 | 日本a在线网址| 国产在线精品亚洲第一网站| 搡老岳熟女国产| 欧美色视频一区免费| 曰老女人黄片| 亚洲狠狠婷婷综合久久图片| 日日摸夜夜添夜夜添小说| 亚洲美女视频黄频| 久久久久国产一级毛片高清牌| 999久久久精品免费观看国产| 精品免费久久久久久久清纯| 久久精品国产综合久久久| 免费电影在线观看免费观看| 精品久久蜜臀av无| 国产精品久久久久久亚洲av鲁大| 一级作爱视频免费观看| 亚洲电影在线观看av| 我要搜黄色片| 在线国产一区二区在线| 深夜精品福利| 国产 一区 欧美 日韩| 日韩欧美精品v在线| 成人av在线播放网站| 黄色日韩在线| 91久久精品国产一区二区成人 | 国语自产精品视频在线第100页| 国产一区二区三区在线臀色熟女| 高潮久久久久久久久久久不卡| 国产不卡一卡二| 三级毛片av免费| 久久精品91蜜桃| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产欧美网| 五月玫瑰六月丁香| 九色国产91popny在线| 亚洲国产精品999在线| 女警被强在线播放| 99久久久亚洲精品蜜臀av| 母亲3免费完整高清在线观看| 9191精品国产免费久久| 午夜福利免费观看在线| 亚洲在线自拍视频| 99热只有精品国产| 高清在线国产一区| 国产伦精品一区二区三区视频9 | 精品国产乱子伦一区二区三区| 色噜噜av男人的天堂激情| 三级毛片av免费| 不卡一级毛片| 天天躁日日操中文字幕| 免费在线观看视频国产中文字幕亚洲| 亚洲九九香蕉| 亚洲午夜理论影院| 国产精品日韩av在线免费观看| 观看免费一级毛片| 丰满的人妻完整版| 午夜福利成人在线免费观看| 国产三级在线视频| 精品一区二区三区四区五区乱码| 国内少妇人妻偷人精品xxx网站 | 人人妻人人澡欧美一区二区| 免费在线观看成人毛片| 91字幕亚洲| 香蕉国产在线看| 黑人巨大精品欧美一区二区mp4| 欧美激情在线99| 欧美日本视频| 国产黄a三级三级三级人| 日韩欧美国产在线观看| 久久精品国产清高在天天线| 久久伊人香网站| 国产成人系列免费观看| 免费看日本二区| 日韩欧美国产一区二区入口| 在线免费观看不下载黄p国产 | 首页视频小说图片口味搜索| 最好的美女福利视频网| 两人在一起打扑克的视频| 一区二区三区高清视频在线| 亚洲色图 男人天堂 中文字幕| 九九在线视频观看精品| 18美女黄网站色大片免费观看| 99热这里只有是精品50| www日本黄色视频网| 日本黄色片子视频| 成年女人看的毛片在线观看| 性色av乱码一区二区三区2| 亚洲av中文字字幕乱码综合| 一二三四社区在线视频社区8| av黄色大香蕉| 丰满人妻一区二区三区视频av | 最近最新中文字幕大全电影3| av中文乱码字幕在线| 熟女人妻精品中文字幕| 老熟妇仑乱视频hdxx| 午夜福利视频1000在线观看| 动漫黄色视频在线观看| 毛片女人毛片| 成年女人看的毛片在线观看| 俺也久久电影网| 黄色成人免费大全| 又大又爽又粗| 亚洲成人免费电影在线观看| 午夜福利视频1000在线观看| 亚洲中文字幕一区二区三区有码在线看 | 两性夫妻黄色片| 国产视频一区二区在线看| 99国产精品一区二区蜜桃av| 色老头精品视频在线观看| 亚洲精品中文字幕一二三四区| 韩国av一区二区三区四区| 欧美一区二区国产精品久久精品| 亚洲国产欧美网| 九九热线精品视视频播放| 看片在线看免费视频| 床上黄色一级片| 免费观看精品视频网站| 久久精品影院6| 国内精品久久久久精免费| 99精品欧美一区二区三区四区| 久久人妻av系列| 久久亚洲精品不卡| 国产人伦9x9x在线观看| 欧美黑人欧美精品刺激| 国产激情欧美一区二区| 国产三级黄色录像| 中文字幕av在线有码专区| 女生性感内裤真人,穿戴方法视频| 欧美最黄视频在线播放免费| 国产三级黄色录像| 黑人操中国人逼视频| 亚洲狠狠婷婷综合久久图片| 欧美乱色亚洲激情| 亚洲aⅴ乱码一区二区在线播放| 特级一级黄色大片| 69av精品久久久久久| www日本在线高清视频| 日本a在线网址| www日本在线高清视频| av视频在线观看入口| 日韩精品中文字幕看吧| 夜夜爽天天搞| 在线a可以看的网站| 丝袜人妻中文字幕| 亚洲一区二区三区色噜噜| 男人和女人高潮做爰伦理| 亚洲一区二区三区色噜噜| 真人一进一出gif抽搐免费| 三级国产精品欧美在线观看 | 久久国产乱子伦精品免费另类| 夜夜夜夜夜久久久久| 日韩精品青青久久久久久| 美女 人体艺术 gogo| 日韩人妻高清精品专区| xxxwww97欧美| 国产一区二区在线观看日韩 | 成人欧美大片| 天天添夜夜摸| 亚洲欧美日韩东京热| 亚洲七黄色美女视频| 岛国在线观看网站| 麻豆av在线久日| 国产精品久久久久久久电影 | 99精品在免费线老司机午夜| 久久久国产精品麻豆| 欧美日本视频| 一个人免费在线观看电影 | 两性夫妻黄色片| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 久久国产乱子伦精品免费另类| 久久久久国产一级毛片高清牌| 亚洲国产看品久久| 国产av一区在线观看免费| 国产精品一及| 欧美一级a爱片免费观看看| 国产av不卡久久| 国产av麻豆久久久久久久| 淫秽高清视频在线观看| 一个人看视频在线观看www免费 | 国产综合懂色| 19禁男女啪啪无遮挡网站| 人人妻人人看人人澡| 香蕉丝袜av| 岛国在线观看网站| 俺也久久电影网| 1000部很黄的大片| 欧美色欧美亚洲另类二区| 久久国产乱子伦精品免费另类| 成人亚洲精品av一区二区| 黄色女人牲交| 99热精品在线国产| 精品久久久久久久毛片微露脸| 欧美乱色亚洲激情| 国产精品久久视频播放| 我的老师免费观看完整版| 三级毛片av免费| 亚洲美女黄片视频| 色尼玛亚洲综合影院| 国产成年人精品一区二区| 国产av麻豆久久久久久久| av在线天堂中文字幕| 免费在线观看视频国产中文字幕亚洲| 麻豆av在线久日| 色哟哟哟哟哟哟| 美女 人体艺术 gogo| 天天躁狠狠躁夜夜躁狠狠躁| 99久久成人亚洲精品观看| 国产高清视频在线观看网站| 国产亚洲精品av在线| 99国产精品一区二区蜜桃av| 亚洲电影在线观看av| 国产精品美女特级片免费视频播放器 | 俺也久久电影网| 一进一出抽搐gif免费好疼| 精品一区二区三区四区五区乱码| 国产不卡一卡二| 久久午夜综合久久蜜桃| 国产午夜精品久久久久久| 亚洲av片天天在线观看| 亚洲国产精品成人综合色| 日韩中文字幕欧美一区二区| 非洲黑人性xxxx精品又粗又长| 97人妻精品一区二区三区麻豆| 国内久久婷婷六月综合欲色啪| 又紧又爽又黄一区二区| 国产淫片久久久久久久久 | 亚洲一区二区三区不卡视频| a级毛片在线看网站| 成人特级黄色片久久久久久久| 亚洲人成电影免费在线| 久久国产乱子伦精品免费另类| 国产久久久一区二区三区| 亚洲av美国av| 91九色精品人成在线观看| 亚洲黑人精品在线| 午夜福利在线观看免费完整高清在 | 国产黄a三级三级三级人| 亚洲无线观看免费| 老鸭窝网址在线观看| 成年女人看的毛片在线观看| 欧美色欧美亚洲另类二区| 午夜a级毛片| 最新美女视频免费是黄的| 19禁男女啪啪无遮挡网站| 无限看片的www在线观看| 欧美成人一区二区免费高清观看 | 少妇的丰满在线观看| 国产免费av片在线观看野外av| 天堂网av新在线| 欧美性猛交╳xxx乱大交人| 国产乱人伦免费视频| 国产亚洲av嫩草精品影院| 大型黄色视频在线免费观看| 色av中文字幕| 深夜精品福利| 国产精品亚洲美女久久久| 中国美女看黄片| 999久久久国产精品视频| 欧美中文日本在线观看视频| 精品乱码久久久久久99久播| 中文资源天堂在线| 这个男人来自地球电影免费观看| 此物有八面人人有两片| 日本三级黄在线观看| 国产成人福利小说| 岛国在线观看网站| 欧美色视频一区免费| 无遮挡黄片免费观看| 丰满人妻熟妇乱又伦精品不卡| 麻豆国产97在线/欧美| 熟妇人妻久久中文字幕3abv| 国产亚洲av嫩草精品影院| 人人妻人人澡欧美一区二区| 欧美大码av| 九九热线精品视视频播放| 欧美日韩精品网址| 国产真实乱freesex| 91av网一区二区| 99热这里只有是精品50| 观看美女的网站| 黑人巨大精品欧美一区二区mp4| 无人区码免费观看不卡| 国产精品香港三级国产av潘金莲| 窝窝影院91人妻| 欧美成人性av电影在线观看| 亚洲成人精品中文字幕电影| av黄色大香蕉| 黄色视频,在线免费观看| av福利片在线观看| 亚洲人与动物交配视频| 男女做爰动态图高潮gif福利片| 99久久精品一区二区三区| 亚洲黑人精品在线| 男人舔女人的私密视频| 欧美中文综合在线视频| 久久久国产精品麻豆| 国语自产精品视频在线第100页| 亚洲成人久久爱视频| 日韩 欧美 亚洲 中文字幕| 高潮久久久久久久久久久不卡| 午夜亚洲福利在线播放| 999精品在线视频| 国产成人一区二区三区免费视频网站| 国产主播在线观看一区二区| 少妇的逼水好多| 麻豆av在线久日| 婷婷精品国产亚洲av在线| 成人18禁在线播放| 国产精品九九99| 精品福利观看| 亚洲人与动物交配视频| 国产蜜桃级精品一区二区三区| 人人妻人人看人人澡| 一边摸一边抽搐一进一小说| 亚洲精品在线观看二区| 制服人妻中文乱码| 日本免费a在线| 欧美日韩黄片免| www国产在线视频色| 老汉色av国产亚洲站长工具| 成年版毛片免费区| 99热6这里只有精品| 国产精品久久久人人做人人爽| 9191精品国产免费久久| 男插女下体视频免费在线播放| 欧美av亚洲av综合av国产av| 亚洲在线观看片| 国产黄色小视频在线观看| 两个人的视频大全免费| 校园春色视频在线观看| 亚洲精品国产精品久久久不卡| 国产精品一区二区精品视频观看| 午夜激情福利司机影院| 久久中文看片网| 特级一级黄色大片| 国产真实乱freesex| 男女做爰动态图高潮gif福利片| 99久久99久久久精品蜜桃| 亚洲五月婷婷丁香| 搡老岳熟女国产| 99久久精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 国产男靠女视频免费网站| 91老司机精品| 亚洲欧美日韩卡通动漫| 国产三级中文精品| 久久热在线av| а√天堂www在线а√下载| 成人av一区二区三区在线看| 日本在线视频免费播放| 国产精品女同一区二区软件 | 亚洲一区高清亚洲精品| 亚洲av成人av| 免费无遮挡裸体视频| aaaaa片日本免费| 久久天堂一区二区三区四区| netflix在线观看网站| www.www免费av| 日韩欧美精品v在线| 怎么达到女性高潮| 狠狠狠狠99中文字幕| 亚洲成人久久性| 欧美另类亚洲清纯唯美| 成年女人永久免费观看视频| 曰老女人黄片| tocl精华| 男女午夜视频在线观看| 真实男女啪啪啪动态图| 制服丝袜大香蕉在线| 黄色 视频免费看| 亚洲第一欧美日韩一区二区三区| 国产成年人精品一区二区| 免费看光身美女| 国产真实乱freesex| 日本 av在线| 一级毛片高清免费大全| 又大又爽又粗| 成人三级黄色视频| 看片在线看免费视频| 国产精品久久视频播放| 国产野战对白在线观看|