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

    利用重測(cè)序技術(shù)開(kāi)發(fā)高粱多態(tài)性SSR分子標(biāo)記

    2019-11-21 11:09:26王平王春語(yǔ)張麗霞叢玲朱振興陸曉春
    生物技術(shù)通報(bào) 2019年11期
    關(guān)鍵詞:多態(tài)拷貝堿基

    王平 王春語(yǔ) 張麗霞 叢玲 朱振興 陸曉春

    (1遼寧省農(nóng)業(yè)科學(xué)院 高粱研究所,沈陽(yáng) 110161;2沈陽(yáng)農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,沈陽(yáng) 110161)

    高粱是世界上重要的糧食作物,是五大作物之一,在人類發(fā)展史上扮演了重要角色[1]。隨著分子生物學(xué)以及生物信息學(xué)的發(fā)展,越來(lái)越多的重要基因被克隆出來(lái),使得分子育種技術(shù)迅速發(fā)展。分子育種技術(shù)以及基因圖位克隆的深入發(fā)展,都離不開(kāi)分子標(biāo)記的開(kāi)發(fā)。分子標(biāo)記技術(shù)發(fā)展至今,已經(jīng)經(jīng)歷過(guò)幾代分子標(biāo)記的歷程,從AFLP,SSR標(biāo)記發(fā)展到今天的SNP標(biāo)記[2-4]。盡管SNP標(biāo)記在很多地方已經(jīng)利用起來(lái),但是SNP標(biāo)記利用需要用到一些特殊儀器設(shè)備,這些設(shè)備價(jià)格不菲,如利用KASP系統(tǒng),儀器昂貴[5];或是利用酶切方法開(kāi)發(fā)CAPS標(biāo)記或是dCAPS標(biāo)記,酶切的方法成本高,分辨率低,酶切不完全等問(wèn)題,經(jīng)常會(huì)對(duì)實(shí)驗(yàn)造成一些困擾[6]。雖然還有其他檢測(cè)SNP的方法,但總體操作相對(duì)復(fù)雜,核心問(wèn)題是SNP標(biāo)記僅有兩種多態(tài)型,以上這些因素都限制了SNP的應(yīng)用。而SSR標(biāo)記即短序列重復(fù)微衛(wèi)星標(biāo)記,因其開(kāi)發(fā)價(jià)格低廉,使用門檻很低,多態(tài)性好而且穩(wěn)定等原因在現(xiàn)階段的實(shí)踐中有很大的運(yùn)用空間[3]。測(cè)序技術(shù)的飛躍發(fā)展,使得高通量測(cè)序成本日益降低,在科研中發(fā)揮越來(lái)越重要的作用[7-8]。利用高通量測(cè)序技術(shù)在很多物種中進(jìn)行了SSR分子標(biāo)記的開(kāi)發(fā)和應(yīng)用,大多使用轉(zhuǎn)錄組測(cè)序的數(shù)據(jù),如在芝麻、豌豆、蘇丹草及馬尾松等作物中都用到RNA-Seq數(shù)據(jù)開(kāi)發(fā)SSR標(biāo)記[9-12]。因轉(zhuǎn)錄組測(cè)序通量、基因時(shí)空表達(dá)以及僅限基因區(qū)域序列等因素影響造成開(kāi)發(fā)的SSR標(biāo)記產(chǎn)量小,質(zhì)量不高。盡管美國(guó)高粱品種BTX623在2009年就完成了基因組測(cè)序[13],也預(yù)測(cè)了很多SSR標(biāo)記[14-15],但越來(lái)越多的研究表明很多參考基因組僅能代表一部分生態(tài)型的多樣性,這也是泛基因組學(xué)研究興起的重要原因[16-17]。來(lái)源于美國(guó)種質(zhì)的BTX623僅能代表部分生態(tài)型品種的序列,而且在國(guó)內(nèi)高粱育種中使用的很多資源與BTX623有較大的差異,造成SSR標(biāo)記的開(kāi)發(fā)在基因組一些區(qū)域的多態(tài)性標(biāo)記開(kāi)發(fā)效率低下,極大限制了高粱SSR分子標(biāo)記的應(yīng)用。

    為了有效開(kāi)發(fā)多態(tài)性SSR標(biāo)記,本研究利用26個(gè)包括不育系和恢復(fù)系的材料,然后進(jìn)行重測(cè)序。測(cè)序完成后,利用生物信息分析工具開(kāi)發(fā)SSR標(biāo)記,進(jìn)行一系列分析后最終開(kāi)發(fā)了兩萬(wàn)多個(gè)SSR標(biāo)記,并部分驗(yàn)證了SSR 標(biāo)記的多態(tài)性,這些標(biāo)記的成功開(kāi)發(fā),為將來(lái)基因定位、克隆、分子育種等實(shí)踐提供了重要的分子標(biāo)記信息。

    1 材料與方法

    1.1 材料

    重測(cè)序?qū)嶒?yàn)材料:高粱生產(chǎn)骨干親本:12份恢復(fù)系和14份不育系(表1)。所有實(shí)驗(yàn)材料于5月上旬在沈陽(yáng)實(shí)驗(yàn)基地種植,行長(zhǎng)4 m,株距15 cm。

    表112份恢復(fù)系和14份不育系樣品信息

    1.2 方法

    1.2.1 分子標(biāo)記開(kāi)發(fā)過(guò)程 生物信息學(xué)分析流程如下:對(duì)12份恢復(fù)系和14份不育系高粱材料進(jìn)行了10×重測(cè)序,然后從兩個(gè)途徑同時(shí)進(jìn)行數(shù)據(jù)分析。利用高通量測(cè)序數(shù)據(jù)中分析變異信息的軟件GATK(Genome Analysis ToolKit)算法分析全基因組indel差異,另一方面利用BTX623參考基因組和簡(jiǎn)單重復(fù)序列鑒定軟件MISA(MIcroSAtellite identification tool)分析全基因組SSR位點(diǎn)信息,然后綜合比較得到全基因組SSR差異位點(diǎn)信息,而后進(jìn)行單拷貝基因的分析,綜合SSR差異位點(diǎn)分析,依據(jù)SSR標(biāo)記重復(fù)單元堿基數(shù)進(jìn)行分析,然后依據(jù)兩端保守序列進(jìn)行引物的設(shè)計(jì)開(kāi)發(fā)。引物設(shè)計(jì)利用Primer 3軟件,遵循引物設(shè)計(jì)原則設(shè)計(jì)正反向擴(kuò)增引物,引物的GC含量在50%-60%之間,長(zhǎng)度在18-23 bp左右,退火溫度設(shè)在57-63℃之間,以58℃為最佳。設(shè)計(jì)時(shí)避免3'末端堿基為A,以及潛在二級(jí)結(jié)構(gòu)的產(chǎn)生。將挑選要進(jìn)行合成的引物,送往蘇州金唯智公司合成(https://www.genewiz.com.cn)。

    1.2.2 DNA提取 實(shí)驗(yàn)材料開(kāi)始抽穗時(shí)剪取葉片1-2 cm到2 mL離心管中,然后加入 600 μL TPS(100 mmol/L Tris-HCl(pH 8.0),10 mmol/L EDTA(pH8.0),1 mol/L KC,高壓濕熱滅菌后使用)提取液;放入一個(gè)不銹鋼珠(直徑約5 mm),在組織磨樣機(jī)(Tissuelyser Ⅱ,QIAGEN,德國(guó))上研磨,頻率設(shè)定為20/s,研磨 1-2 min;利用強(qiáng)力磁鐵取出鋼珠,將離心管放入 65℃水浴中保溫 30 min;12 000r/min 離心 5 min,小心吸出上清液轉(zhuǎn)移至新的 1.5 mL 離心管,加入等體積的異丙醇,輕輕顛倒混勻,clss= r/min 離心10 min,棄上清。加入 1 mL 70%乙醇洗滌沉淀,12000 r/min 離心 5 min,棄上清,留沉淀,打開(kāi)管蓋,待沉淀完全干燥后加入 50 μL包含0.2 mg/mL Rnase A(生工生物工程有限公司,上海,中國(guó))的滅菌水或是1×TE buffer溶解沉淀。提取的高粱葉片基因組DNA在0.8%-1.0%瓊脂糖電泳檢測(cè)合格后保存于-20℃冰箱備用。

    1.2.3 SSR電泳分析 PCR及聚丙烯酰胺凝膠電泳:高粱基因組 DNA 1 μL,F(xiàn)/R-primer(10 μmol/L)0.2μL,2×Taq PCR Master Mix 5 μL(北京艾德萊生物),ddH2O 3.8 μL總體積10 μL;反應(yīng)在PCR儀上(T100 Thermal cycler,BIO-RAD,美國(guó))進(jìn)行擴(kuò)增,PCR擴(kuò)增程序如下:95℃預(yù)變性3 min;94℃ 30 s、X℃30 s、72℃ 25 s,34個(gè)循環(huán),最后72℃ 延伸5 min;其中退火溫度X隨引物設(shè)計(jì)適度調(diào)整。PCR產(chǎn)物取2-2.5 μL用6%非變性聚丙烯酰胺凝膠電泳1.5-2.0 h左右(根據(jù)差異大小來(lái)決定電泳時(shí)間)。電泳結(jié)束后,關(guān)閉電源,從電泳槽中將膠板取下,用鏟子將兩板分開(kāi),緩慢將膠取下,并記錄膠的順序編號(hào),經(jīng)銀染、顯色、水洗等程序獲得擴(kuò)增條帶,然后在燈箱板上進(jìn)行膠片的拍攝(EOS 5D,Canon,日本)。

    2 結(jié)果

    2.126份重測(cè)序材料多態(tài)性SSR標(biāo)記篩選

    利用高粱中26份具有代表性的不育系和恢復(fù)系材料,進(jìn)行重測(cè)序序列分析。重測(cè)序完成后,進(jìn)行SSR標(biāo)記的篩選,為考慮將來(lái)SSR分子標(biāo)記檢測(cè)的便利性,本研究制定了比較嚴(yán)格的篩選標(biāo)準(zhǔn)。SSR篩選標(biāo)準(zhǔn)如下所示:(1)兩堿基單元重復(fù)SSR標(biāo)記的重復(fù)次數(shù)需≥6次重復(fù),即這個(gè)SSR標(biāo)記長(zhǎng)度必須在12 bp以上;三堿基重復(fù)序列SSR重復(fù)次數(shù)≥5次;四堿基重復(fù)序列SSR的重復(fù)次數(shù)≥4次;五堿基重復(fù)序列SSR的重復(fù)次數(shù)≥3;六堿基重復(fù)序列SSR的重復(fù)次數(shù)≥3次。(2)在所有26個(gè)高粱品種中要有多態(tài)性,至少是兩種。經(jīng)過(guò)了本研究的標(biāo)準(zhǔn)篩選,共得到24441個(gè)多態(tài)性SSR,以2或是3個(gè)基序重復(fù)的SSR為主,兩堿基重復(fù)SSR標(biāo)記11051個(gè);三堿基重復(fù)6082個(gè);四堿基重復(fù)2600個(gè);五堿基重復(fù)2449個(gè);六堿基重復(fù)2259個(gè)(圖1-A)。對(duì)篩選到的SSR標(biāo)記在26個(gè)品種間進(jìn)行基因型比較,基因型多態(tài)類型在2-7種之間,最多能達(dá)到7種多態(tài)型,大部分的標(biāo)記僅有兩種,然后多態(tài)性數(shù)目呈斷崖式下降,以2種多態(tài)性的為主,高達(dá)68%;其次是3種多態(tài)型的占21%;而其它幾種多態(tài)型就相對(duì)較少分布,7種多態(tài)型的標(biāo)記僅有77個(gè)(圖1-B)。

    圖1 不同重復(fù)單元堿基數(shù)的SSR 標(biāo)記的數(shù)量(A)及其在26個(gè)品種間的不同多態(tài)型標(biāo)記數(shù)量(B)

    2.2 單拷貝基因處多態(tài)性SSR分子標(biāo)記的篩選

    在分子育種實(shí)踐中,功能性標(biāo)記起到至關(guān)重要的作用。功能性標(biāo)記,主要是位于基因位置上的標(biāo)記,而一般這樣的基因,幾乎都是單拷貝基因。所以本研究的策略是篩選在單拷貝基因位置的SSR標(biāo)記。本研究篩選首先是進(jìn)行單拷貝基因的鑒定,然后分析這些單拷貝基因(包含該基因上下游)序列中所包含的SSR序列,篩選位于單拷貝基因(含上下游2 k)內(nèi)的多態(tài)性SSR標(biāo)記。相對(duì)所有的SSR標(biāo)記數(shù)量,位于基因區(qū)標(biāo)記的SSR數(shù)量,僅為總量的1/4左右,總計(jì)篩選到6733個(gè)標(biāo)記。2-3個(gè)基序重復(fù)的SSR標(biāo)記,僅有2000左右,而4-6重復(fù)的SSR標(biāo)記,約700-800個(gè)(圖2-A)。大部分的單拷貝基因上SSR標(biāo)記在26個(gè)品種僅有2種基因型,占了總SSR的69%;其次為3種基因型占總數(shù)的20%,4-7種基因型的SSR標(biāo)記比較少,7種基因型的SSR標(biāo)記僅為28個(gè)(圖2-B)。

    圖2 單拷貝基因處不同重復(fù)單元堿基數(shù)的SSR標(biāo)記數(shù)量(A)和26個(gè)品種間不同多態(tài)型標(biāo)記數(shù)量(B)

    2.3 單拷貝基因處多態(tài)性SSR分子標(biāo)記在基因上位置分析

    通過(guò)分析這些SSR標(biāo)記在基因上的位置,6733個(gè)SSR標(biāo)記中,位于編碼區(qū)的僅有820個(gè),UTR區(qū)1276個(gè),內(nèi)含子725個(gè),大多在上下游區(qū)域共3932個(gè)。如在2個(gè)基序重復(fù)的SSR中,占有63%,而分布在編碼區(qū)多為3堿基或6堿基重復(fù)(圖3-A)。這些SSR標(biāo)記在26個(gè)品種間的基因型數(shù),可以看出在基因上下游區(qū)域的SSR標(biāo)記居多,位于基因編碼區(qū)的最少,約占所有在編碼區(qū)附近的12%。大多數(shù)SSR在26個(gè)高粱品種有兩種基因型,依次遞減到7種多態(tài)型(圖3-B)。

    2.4 SSR標(biāo)記有效性評(píng)價(jià)

    為了檢測(cè)開(kāi)發(fā)SSR標(biāo)記的有效性,本研究在位于編碼區(qū)約800個(gè)SSR標(biāo)記引物中,在每條染色體上挑選了5對(duì)相對(duì)均一覆蓋整條染色體的SSR標(biāo)記,在10條染色體上進(jìn)行了引物的挑選、合成,共計(jì)50對(duì)引物,引物相關(guān)信息如表2所示。為了檢驗(yàn)引物的實(shí)用性,利用了參考基因組高粱材料BTX623和雜交種遼糯3號(hào),因?yàn)楸狙芯吭谟N實(shí)踐中很多時(shí)候利用的是雜交種材料。在50對(duì)引物中,兩個(gè)材料中僅有一對(duì)引物未擴(kuò)增出產(chǎn)物(Chr6-3兩個(gè)親本間無(wú)擴(kuò)增),引物成功率高達(dá)98%,另有兩對(duì)引物僅能在一個(gè)材料中進(jìn)行擴(kuò)增;另一個(gè)材料顯示無(wú)擴(kuò)增產(chǎn)物(Chr2-3,Chr2-5都是在BTX623中無(wú)擴(kuò)增,在遼糯3中能正常擴(kuò)增),最終顯示有18對(duì)引物在這兩個(gè)親本間存在多態(tài)性。

    為進(jìn)一步檢驗(yàn)該批引物的質(zhì)量,本研究同時(shí)挑選了引物Chr1-1,分析其在74份核心種質(zhì)資源中的多態(tài)性豐富程度。該引物在BTX623中擴(kuò)增的產(chǎn)物的大小,應(yīng)該是266 bp,擴(kuò)增產(chǎn)物中較濃的條帶,大小基本符合預(yù)期,該條帶即是目標(biāo)條帶。條帶大小比對(duì)分析表明,該標(biāo)記在74份核心種質(zhì)資源中至少有8種多態(tài)型,標(biāo)記Chr1-1具有比較高的分辨能力(圖4)。

    圖3 位于單拷貝基因上不同重復(fù)單元堿基數(shù)的SSR標(biāo)記位于基因不同位置的數(shù)量(A)及其在26個(gè)品種間的多態(tài)型數(shù)目(B)

    3 討論

    SSR分子標(biāo)記因其價(jià)廉物美,儀器設(shè)備要求不高、多態(tài)性好、實(shí)踐中可操作性強(qiáng)等特點(diǎn),一直是分子標(biāo)記領(lǐng)域的熱門開(kāi)發(fā)方向[3,18-19]。近年來(lái),由于高通量測(cè)序技術(shù)的發(fā)展,成本日益降級(jí),利用高通量測(cè)序開(kāi)發(fā)分子標(biāo)記,越來(lái)越高效和便捷。在很多研究中利用轉(zhuǎn)錄組數(shù)據(jù)開(kāi)發(fā)SSR標(biāo)記,由于轉(zhuǎn)錄組數(shù)據(jù)本身的缺陷,如測(cè)序主要是基因轉(zhuǎn)錄區(qū)段,

    很多基因具有時(shí)空表達(dá),于是會(huì)漏掉很多基因組序列信息,造成開(kāi)發(fā)多態(tài)性SSR標(biāo)記效率并不高[9-12]。

    表2 分布于10條染色體上的50對(duì)SSR引物信息

    圖4 引物Chr1-1在74份高粱微核心種質(zhì)資源中多態(tài)性分析

    在高粱作物中盡管BTX623的基因組序列在2009 年已經(jīng)完成[13],Yonemaru 等[15]利用該基因組序列開(kāi)發(fā)多態(tài)性SSR標(biāo)記,共獲得5599個(gè)非冗余SSR標(biāo)記,其中(AT/TA)n占所有SSR的26.1%,其次是(AG/TC)n 占20.5%、(AC/TG)n占13.7%和(CG/GC)n占11.8%。其中5012個(gè)SSR標(biāo)記染色體位置是通過(guò)e-PCR技術(shù)與34008個(gè)基因座的預(yù)測(cè)位置進(jìn)行比較來(lái)確定的。在通過(guò)片段分析驗(yàn)證的970個(gè)標(biāo)記中,67.8%(970個(gè)標(biāo)記中的658個(gè))的標(biāo)記成功地在高粱品系BTx623中擴(kuò)增,在11個(gè)高粱品系和一個(gè)蘇丹草品系的組合中,所有SSR基因座的平均多態(tài)性率為45.1%(658個(gè)標(biāo)記中的297個(gè))。BTX623來(lái)源于美國(guó),其基因組序列僅能代表部分生態(tài)型品種的序列,國(guó)內(nèi)高粱育種中使用的很多資源與BTX623有較大的差異,且遺傳范圍狹窄,造成SSR標(biāo)記的開(kāi)發(fā)在基因組一些區(qū)域的多態(tài)性標(biāo)記開(kāi)發(fā)效率低下,極大限制了高粱SSR分子標(biāo)記在國(guó)內(nèi)高粱育種上的應(yīng)用。

    為了能開(kāi)發(fā)多態(tài)性高的SSR標(biāo)記,本研究利用實(shí)驗(yàn)室先前進(jìn)行了26份不育系和恢復(fù)系材料重測(cè)序,利用這26份材料的重測(cè)序數(shù)據(jù)進(jìn)行SSR標(biāo)記的開(kāi)發(fā),有效提高了多態(tài)性標(biāo)記開(kāi)發(fā)的力度和效率。本研究設(shè)計(jì)開(kāi)發(fā)的50對(duì)SSR引物在兩個(gè)親本BTX623和雜交種遼糯3中進(jìn)行擴(kuò)增,有49對(duì)成功擴(kuò)增出條帶,其中有兩對(duì)引物僅在遼糯3號(hào)中能夠擴(kuò)增,表明很有可能某些基因組區(qū)域在BTX623中是缺失的,僅存在國(guó)內(nèi)某些高粱材料中。50對(duì)引物在這兩個(gè)品種中有多態(tài)性的SSR標(biāo)記為36%,并挑選了SSR標(biāo)記Chr1-1在74份微核心種質(zhì)資源中進(jìn)行分析,初步分析該標(biāo)記在這些種質(zhì)資源中有8種多態(tài)型之多,表明從開(kāi)發(fā)的SSR標(biāo)記中容易篩選出多態(tài)性好的標(biāo)記。

    本研究開(kāi)發(fā)的SSR分子標(biāo)記是以26份材料中必須至少有兩種多態(tài)型作為篩選標(biāo)準(zhǔn)。在這些入選的SSR標(biāo)記中,大部分入選的SSR標(biāo)記都是以2個(gè)堿基基序?yàn)橹貜?fù)單元的分子標(biāo)記,3個(gè)以上堿基重復(fù)的達(dá)到6082,26個(gè)品種中多態(tài)性超過(guò)2種的有24441,分布在單拷貝基因上SSR標(biāo)記占總標(biāo)記數(shù)1/4,分布單拷貝基因編碼區(qū)的標(biāo)記占落在單拷貝基因上的12%。這些標(biāo)記類型和特點(diǎn)為高粱基因定位、克隆以及分子育種等分子標(biāo)記的開(kāi)發(fā)提供了可靠的參考數(shù)據(jù)。

    4 結(jié)論

    利用14個(gè)不育系和12個(gè)恢復(fù)系作為重測(cè)序材料,開(kāi)發(fā)了在26份材料中至少含有2種多態(tài)型的SSR標(biāo)記24441個(gè)單拷貝基因處的多態(tài)性SSR 6733個(gè)。隨機(jī)挑選均勻覆蓋10條染色體的單拷貝基因處SSR標(biāo)記50對(duì),利用遼寧高粱雜交種遼糯3號(hào)進(jìn)行測(cè)試,其中49對(duì)能擴(kuò)增出產(chǎn)物,成功率高達(dá)98%,利用高粱雜交種遼糯3號(hào)、BTX623和74份微核心種質(zhì)資源測(cè)試表明,50對(duì)SSR標(biāo)記在2個(gè)品種中有18對(duì)表現(xiàn)出多態(tài)性,挑選了一對(duì)引物在74份微核心種質(zhì)中可見(jiàn)8種多態(tài)型。本研究表明利用不育系和恢復(fù)系材料進(jìn)行重測(cè)序能有效開(kāi)發(fā)多態(tài)性高的SSR標(biāo)記。

    猜你喜歡
    多態(tài)拷貝堿基
    分層多態(tài)加權(quán)k/n系統(tǒng)的可用性建模與設(shè)計(jì)優(yōu)化
    應(yīng)用思維進(jìn)階構(gòu)建模型 例談培養(yǎng)學(xué)生創(chuàng)造性思維
    參差多態(tài)而功不唐捐
    中國(guó)科學(xué)家創(chuàng)建出新型糖基化酶堿基編輯器
    生命“字母表”迎來(lái)4名新成員
    生命“字母表”迎來(lái)4名新成員
    中國(guó)生殖健康(2018年1期)2018-11-06 07:14:38
    人多巴胺D2基因啟動(dòng)子區(qū)—350A/G多態(tài)位點(diǎn)熒光素酶表達(dá)載體的構(gòu)建與鑒定及活性檢測(cè)
    煙堿型乙酰膽堿受體基因多態(tài)與早發(fā)性精神分裂癥的關(guān)聯(lián)研究
    文件拷貝誰(shuí)最“給力”
    精品第一国产精品| 女人高潮潮喷娇喘18禁视频| 成年人免费黄色播放视频| 一级,二级,三级黄色视频| 女人精品久久久久毛片| av电影中文网址| 黄色 视频免费看| 久久亚洲国产成人精品v| 国产福利在线免费观看视频| 亚洲国产中文字幕在线视频| 亚洲欧美一区二区三区国产| 成人午夜精彩视频在线观看| 老鸭窝网址在线观看| 777米奇影视久久| 色播在线永久视频| 91字幕亚洲| 这个男人来自地球电影免费观看| 多毛熟女@视频| 99久久综合免费| 欧美成人精品欧美一级黄| 人人妻人人爽人人添夜夜欢视频| 国产精品一国产av| 欧美日韩视频高清一区二区三区二| 亚洲欧美精品综合一区二区三区| 久久亚洲精品不卡| 这个男人来自地球电影免费观看| 黄色a级毛片大全视频| 一边摸一边做爽爽视频免费| 91精品三级在线观看| 免费人妻精品一区二区三区视频| 人人澡人人妻人| 秋霞在线观看毛片| 另类亚洲欧美激情| 久久九九热精品免费| 中文字幕最新亚洲高清| 欧美日韩黄片免| 赤兔流量卡办理| 精品福利永久在线观看| 日韩av免费高清视频| 欧美性长视频在线观看| 午夜福利乱码中文字幕| 国产精品 国内视频| 狂野欧美激情性xxxx| 精品久久久精品久久久| 欧美在线黄色| 中国美女看黄片| 黑人猛操日本美女一级片| 热99国产精品久久久久久7| 久久久国产欧美日韩av| 国产成人精品久久二区二区91| xxx大片免费视频| av不卡在线播放| av在线老鸭窝| 国产成人精品无人区| 91麻豆av在线| 欧美成人午夜精品| 精品少妇黑人巨大在线播放| 午夜老司机福利片| 咕卡用的链子| 波多野结衣av一区二区av| 国产伦人伦偷精品视频| 久久久久国产一级毛片高清牌| 一本大道久久a久久精品| 精品亚洲成a人片在线观看| 99精国产麻豆久久婷婷| 一边亲一边摸免费视频| 欧美日韩精品网址| 亚洲av美国av| 操出白浆在线播放| 色播在线永久视频| 精品视频人人做人人爽| 最新的欧美精品一区二区| 国产日韩一区二区三区精品不卡| 欧美日韩亚洲高清精品| 久久中文字幕一级| 色综合欧美亚洲国产小说| 国产精品熟女久久久久浪| 欧美精品高潮呻吟av久久| 狠狠婷婷综合久久久久久88av| 国产成人91sexporn| 欧美日韩国产mv在线观看视频| 97在线人人人人妻| 男女下面插进去视频免费观看| 日韩av免费高清视频| 777久久人妻少妇嫩草av网站| 久久亚洲精品不卡| 麻豆国产av国片精品| 一区二区三区四区激情视频| 亚洲 国产 在线| 亚洲成色77777| 国产一区二区在线观看av| 久久精品国产综合久久久| 亚洲精品国产av蜜桃| 国产91精品成人一区二区三区 | 亚洲欧美日韩高清在线视频 | 男女午夜视频在线观看| 亚洲av成人不卡在线观看播放网 | 又大又黄又爽视频免费| 又黄又粗又硬又大视频| 少妇被粗大的猛进出69影院| 18禁黄网站禁片午夜丰满| 亚洲国产最新在线播放| 精品一品国产午夜福利视频| 免费人妻精品一区二区三区视频| 丰满少妇做爰视频| 欧美黑人精品巨大| 精品熟女少妇八av免费久了| 超碰成人久久| 国产伦理片在线播放av一区| 老司机午夜十八禁免费视频| 一级a爱视频在线免费观看| 亚洲国产精品999| 成年人黄色毛片网站| 免费一级毛片在线播放高清视频 | 久久久久久久国产电影| 男女下面插进去视频免费观看| 精品一区在线观看国产| 美女视频免费永久观看网站| 波野结衣二区三区在线| av一本久久久久| 妹子高潮喷水视频| 亚洲九九香蕉| 高清av免费在线| 一边摸一边做爽爽视频免费| 欧美日韩综合久久久久久| 国产成人欧美在线观看 | 观看av在线不卡| 国产免费一区二区三区四区乱码| 日韩精品免费视频一区二区三区| 一区在线观看完整版| 中文字幕高清在线视频| 日韩视频在线欧美| 制服诱惑二区| 午夜91福利影院| 精品一区二区三卡| 香蕉国产在线看| 在线观看免费高清a一片| 新久久久久国产一级毛片| 男女下面插进去视频免费观看| 亚洲国产精品一区三区| 国产欧美日韩一区二区三区在线| 国产福利在线免费观看视频| 国产99久久九九免费精品| 国产在线一区二区三区精| 国产成人精品久久久久久| 亚洲av国产av综合av卡| 91精品国产国语对白视频| 精品亚洲成a人片在线观看| 久久久久久久久久久久大奶| 热re99久久国产66热| 国产真人三级小视频在线观看| 国产黄色免费在线视频| 国产成人a∨麻豆精品| 国产又爽黄色视频| 18在线观看网站| 久久 成人 亚洲| 欧美av亚洲av综合av国产av| 波多野结衣一区麻豆| 亚洲中文日韩欧美视频| 日本91视频免费播放| 叶爱在线成人免费视频播放| av一本久久久久| 久久久久久免费高清国产稀缺| 黄色视频不卡| av又黄又爽大尺度在线免费看| 老司机在亚洲福利影院| 赤兔流量卡办理| 在线观看免费日韩欧美大片| 99久久99久久久精品蜜桃| 国产日韩欧美在线精品| 国产色视频综合| 久久青草综合色| 少妇 在线观看| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲国产一区二区在线观看 | 18在线观看网站| 国产一区二区三区av在线| 久久久精品国产亚洲av高清涩受| xxx大片免费视频| www日本在线高清视频| 97精品久久久久久久久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 久久毛片免费看一区二区三区| 久久久久精品人妻al黑| 老鸭窝网址在线观看| 精品国产乱码久久久久久小说| 一边亲一边摸免费视频| 一级毛片我不卡| 国产一级毛片在线| 免费少妇av软件| 黄色一级大片看看| 久久精品国产亚洲av高清一级| 国产国语露脸激情在线看| 久久久久精品国产欧美久久久 | 欧美 亚洲 国产 日韩一| 久久精品aⅴ一区二区三区四区| 麻豆国产av国片精品| 又黄又粗又硬又大视频| 美女高潮到喷水免费观看| 久久精品国产亚洲av高清一级| 亚洲黑人精品在线| 男人操女人黄网站| 一级毛片电影观看| 一区二区三区乱码不卡18| av在线app专区| 在线天堂中文资源库| 成在线人永久免费视频| 麻豆av在线久日| 欧美日韩福利视频一区二区| www.熟女人妻精品国产| 99国产精品一区二区三区| 久久久久久久国产电影| 国产真人三级小视频在线观看| 国产一区二区三区综合在线观看| 国产熟女欧美一区二区| 欧美日韩视频精品一区| 天堂8中文在线网| 久久国产精品影院| 青春草视频在线免费观看| 午夜福利在线免费观看网站| 免费观看a级毛片全部| 多毛熟女@视频| 999精品在线视频| 天天影视国产精品| 美女国产高潮福利片在线看| 成年av动漫网址| 欧美日本中文国产一区发布| 99热国产这里只有精品6| 久久国产亚洲av麻豆专区| 色婷婷久久久亚洲欧美| 亚洲熟女毛片儿| 巨乳人妻的诱惑在线观看| 捣出白浆h1v1| 纵有疾风起免费观看全集完整版| 精品国产一区二区久久| 熟女少妇亚洲综合色aaa.| 脱女人内裤的视频| 久久久久网色| 涩涩av久久男人的天堂| 久久久欧美国产精品| 美女高潮到喷水免费观看| 18在线观看网站| 搡老乐熟女国产| 久久国产精品影院| 国产真人三级小视频在线观看| 精品国产一区二区三区久久久樱花| av又黄又爽大尺度在线免费看| 午夜日韩欧美国产| 少妇人妻 视频| 中文字幕av电影在线播放| 国产高清视频在线播放一区 | 妹子高潮喷水视频| 午夜福利一区二区在线看| 国产精品偷伦视频观看了| 国产成人a∨麻豆精品| www.自偷自拍.com| 国产亚洲精品久久久久5区| 啦啦啦在线观看免费高清www| 欧美精品啪啪一区二区三区 | 超碰成人久久| 在线精品无人区一区二区三| 日韩av在线免费看完整版不卡| 午夜视频精品福利| 日韩,欧美,国产一区二区三区| av福利片在线| 99久久99久久久精品蜜桃| 涩涩av久久男人的天堂| 少妇裸体淫交视频免费看高清 | av天堂久久9| 久久久久久免费高清国产稀缺| 日韩熟女老妇一区二区性免费视频| 99热国产这里只有精品6| 王馨瑶露胸无遮挡在线观看| 丁香六月天网| 久久精品人人爽人人爽视色| 日韩一卡2卡3卡4卡2021年| 欧美日韩精品网址| 亚洲精品一二三| 男人操女人黄网站| 丁香六月欧美| 国产淫语在线视频| 日韩大码丰满熟妇| 精品一区二区三区av网在线观看 | 大话2 男鬼变身卡| 免费久久久久久久精品成人欧美视频| 99热全是精品| 在线天堂中文资源库| 一级毛片 在线播放| 韩国高清视频一区二区三区| 国产精品二区激情视频| 欧美精品一区二区免费开放| 免费在线观看日本一区| 每晚都被弄得嗷嗷叫到高潮| 亚洲成色77777| av福利片在线| 9热在线视频观看99| 在线天堂中文资源库| 菩萨蛮人人尽说江南好唐韦庄| 王馨瑶露胸无遮挡在线观看| 97在线人人人人妻| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区在线不卡| a级毛片黄视频| 国产精品二区激情视频| 国产精品一国产av| 午夜91福利影院| avwww免费| 久久久精品国产亚洲av高清涩受| 2021少妇久久久久久久久久久| 捣出白浆h1v1| 桃花免费在线播放| 在线观看www视频免费| 99国产精品免费福利视频| 高清av免费在线| 亚洲精品第二区| www.av在线官网国产| 精品久久久久久电影网| 国语对白做爰xxxⅹ性视频网站| 老司机靠b影院| 亚洲男人天堂网一区| 天天操日日干夜夜撸| e午夜精品久久久久久久| 高清不卡的av网站| 亚洲精品国产av蜜桃| 男女之事视频高清在线观看 | 精品少妇内射三级| 天天躁狠狠躁夜夜躁狠狠躁| 肉色欧美久久久久久久蜜桃| 成人国产av品久久久| 午夜福利免费观看在线| 亚洲伊人久久精品综合| 午夜福利视频精品| 久久久国产欧美日韩av| 人人妻人人爽人人添夜夜欢视频| 亚洲国产av新网站| 久久久久国产精品人妻一区二区| cao死你这个sao货| 亚洲精品一卡2卡三卡4卡5卡 | 免费人妻精品一区二区三区视频| 一二三四在线观看免费中文在| 99国产精品一区二区蜜桃av | 成人手机av| 亚洲精品第二区| 久久久亚洲精品成人影院| 五月天丁香电影| 日本a在线网址| 人妻人人澡人人爽人人| 香蕉丝袜av| av国产久精品久网站免费入址| 午夜福利免费观看在线| 人妻人人澡人人爽人人| 日本a在线网址| 高清欧美精品videossex| 久久国产精品大桥未久av| 精品少妇内射三级| 叶爱在线成人免费视频播放| 99国产精品一区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲国产欧美网| 亚洲国产精品999| 一区在线观看完整版| 成人免费观看视频高清| 伊人亚洲综合成人网| 日韩中文字幕视频在线看片| 亚洲专区国产一区二区| 国产亚洲欧美在线一区二区| 麻豆国产av国片精品| 99国产精品99久久久久| 精品一品国产午夜福利视频| 国产精品三级大全| 国产一卡二卡三卡精品| 高潮久久久久久久久久久不卡| 国产成人一区二区在线| 日本av手机在线免费观看| 十八禁网站网址无遮挡| 精品高清国产在线一区| 久久人妻福利社区极品人妻图片 | 国产男女超爽视频在线观看| 欧美人与性动交α欧美软件| 欧美在线一区亚洲| 久久人人爽人人片av| 国产成人影院久久av| 亚洲精品第二区| 国产精品国产三级专区第一集| 久热爱精品视频在线9| 中国美女看黄片| 国产精品国产三级国产专区5o| 大香蕉久久成人网| 午夜视频精品福利| 色婷婷av一区二区三区视频| 国产免费又黄又爽又色| 丰满人妻熟妇乱又伦精品不卡| 久久久精品免费免费高清| 免费av中文字幕在线| 精品人妻在线不人妻| 欧美国产精品一级二级三级| 国产一卡二卡三卡精品| 久久青草综合色| 五月开心婷婷网| 国产成人a∨麻豆精品| 夫妻性生交免费视频一级片| 夜夜骑夜夜射夜夜干| 亚洲七黄色美女视频| 女性被躁到高潮视频| 纵有疾风起免费观看全集完整版| 一边摸一边做爽爽视频免费| 丝袜在线中文字幕| 午夜老司机福利片| 天天影视国产精品| 久久久久精品国产欧美久久久 | 国产高清不卡午夜福利| 老司机亚洲免费影院| 一级片免费观看大全| 亚洲精品成人av观看孕妇| 国产色视频综合| 在线亚洲精品国产二区图片欧美| 极品人妻少妇av视频| 少妇精品久久久久久久| 国产精品免费大片| 人人妻人人爽人人添夜夜欢视频| 99国产精品99久久久久| 久久人妻熟女aⅴ| 黄色视频在线播放观看不卡| 免费人妻精品一区二区三区视频| 狠狠精品人妻久久久久久综合| 日韩制服骚丝袜av| 美女午夜性视频免费| 亚洲欧洲日产国产| 久久精品国产亚洲av涩爱| 啦啦啦在线免费观看视频4| 亚洲人成77777在线视频| 欧美人与性动交α欧美软件| 中国美女看黄片| 国产精品久久久久久精品古装| 嫁个100分男人电影在线观看 | 少妇的丰满在线观看| 日韩 欧美 亚洲 中文字幕| 国产高清videossex| 亚洲久久久国产精品| 亚洲情色 制服丝袜| 九草在线视频观看| 欧美人与性动交α欧美软件| 亚洲一区二区三区欧美精品| 每晚都被弄得嗷嗷叫到高潮| 欧美人与性动交α欧美精品济南到| 97精品久久久久久久久久精品| 国产欧美日韩精品亚洲av| 亚洲国产看品久久| 午夜久久久在线观看| 国产主播在线观看一区二区 | 免费看不卡的av| 黄片播放在线免费| 久久女婷五月综合色啪小说| 国产老妇伦熟女老妇高清| 汤姆久久久久久久影院中文字幕| 少妇粗大呻吟视频| e午夜精品久久久久久久| 精品视频人人做人人爽| 我的亚洲天堂| 国产熟女欧美一区二区| 侵犯人妻中文字幕一二三四区| 久久精品亚洲av国产电影网| 老熟女久久久| 国产精品免费大片| 久久鲁丝午夜福利片| 午夜福利影视在线免费观看| 日韩人妻精品一区2区三区| 国产精品久久久久久精品古装| 国产麻豆69| 国产淫语在线视频| 看免费av毛片| 性色av乱码一区二区三区2| 久热爱精品视频在线9| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品麻豆人妻色哟哟久久| 精品一区在线观看国产| 国产精品.久久久| 老司机深夜福利视频在线观看 | 青春草亚洲视频在线观看| 成年人免费黄色播放视频| 亚洲精品国产av成人精品| 男女边摸边吃奶| 亚洲国产成人一精品久久久| 久久久久视频综合| 亚洲av美国av| av电影中文网址| 亚洲精品日韩在线中文字幕| 亚洲精品国产av蜜桃| 91麻豆av在线| 在线精品无人区一区二区三| 又黄又粗又硬又大视频| 亚洲五月婷婷丁香| 精品人妻一区二区三区麻豆| 激情五月婷婷亚洲| 色94色欧美一区二区| 91九色精品人成在线观看| 人成视频在线观看免费观看| 蜜桃国产av成人99| 国产精品亚洲av一区麻豆| 18禁国产床啪视频网站| 咕卡用的链子| 十八禁人妻一区二区| av线在线观看网站| 男女国产视频网站| 波野结衣二区三区在线| 一级a爱视频在线免费观看| 日韩伦理黄色片| 国产精品久久久av美女十八| 50天的宝宝边吃奶边哭怎么回事| 肉色欧美久久久久久久蜜桃| 在线观看www视频免费| 免费不卡黄色视频| 久久毛片免费看一区二区三区| 午夜福利视频精品| 色精品久久人妻99蜜桃| 亚洲图色成人| 欧美精品av麻豆av| 免费高清在线观看日韩| 久久人人爽人人片av| 极品人妻少妇av视频| 国产午夜精品一二区理论片| 精品卡一卡二卡四卡免费| 电影成人av| 亚洲 国产 在线| 国产精品三级大全| 视频区欧美日本亚洲| 国产精品99久久99久久久不卡| 免费在线观看黄色视频的| 国产精品香港三级国产av潘金莲 | 亚洲国产精品一区二区三区在线| 国产成人欧美| 2018国产大陆天天弄谢| 女人久久www免费人成看片| a 毛片基地| 视频区图区小说| 午夜福利一区二区在线看| 中文字幕高清在线视频| 视频在线观看一区二区三区| 国产爽快片一区二区三区| 亚洲欧洲日产国产| 久久女婷五月综合色啪小说| 欧美成人午夜精品| 首页视频小说图片口味搜索 | 亚洲国产最新在线播放| 只有这里有精品99| 波多野结衣av一区二区av| 黄色a级毛片大全视频| 51午夜福利影视在线观看| 国产无遮挡羞羞视频在线观看| 可以免费在线观看a视频的电影网站| 中文字幕人妻丝袜一区二区| 国产精品熟女久久久久浪| 免费观看a级毛片全部| 美女中出高潮动态图| 国产免费视频播放在线视频| 母亲3免费完整高清在线观看| 自线自在国产av| 久久ye,这里只有精品| 欧美日韩亚洲综合一区二区三区_| 久久久精品区二区三区| 国产男人的电影天堂91| 女人久久www免费人成看片| 久久国产精品影院| 亚洲免费av在线视频| 一本一本久久a久久精品综合妖精| 亚洲精品国产av成人精品| 欧美日韩视频精品一区| 一区福利在线观看| 99久久综合免费| 久久99精品国语久久久| 少妇精品久久久久久久| 亚洲,一卡二卡三卡| 色精品久久人妻99蜜桃| 91成人精品电影| 天堂8中文在线网| 纵有疾风起免费观看全集完整版| www.自偷自拍.com| 国产日韩一区二区三区精品不卡| 成年美女黄网站色视频大全免费| 亚洲av综合色区一区| 啦啦啦在线观看免费高清www| 女人久久www免费人成看片| 久久人人爽人人片av| 久久久久视频综合| 美女高潮到喷水免费观看| 国产亚洲av高清不卡| 性少妇av在线| 亚洲精品美女久久久久99蜜臀 | 亚洲熟女精品中文字幕| 在线观看免费视频网站a站| 99九九在线精品视频| 国产一区二区三区av在线| 一区福利在线观看| 国产成人a∨麻豆精品| 大香蕉久久成人网| 男女无遮挡免费网站观看| 99久久99久久久精品蜜桃| cao死你这个sao货| 丝瓜视频免费看黄片| 欧美久久黑人一区二区| 色网站视频免费| 飞空精品影院首页| 国产视频首页在线观看| 亚洲欧美一区二区三区黑人| 国产欧美日韩一区二区三区在线| 中文字幕最新亚洲高清| 女性生殖器流出的白浆| 侵犯人妻中文字幕一二三四区| 可以免费在线观看a视频的电影网站| 女人久久www免费人成看片| 黄网站色视频无遮挡免费观看| av又黄又爽大尺度在线免费看| 成人亚洲精品一区在线观看| 啦啦啦啦在线视频资源| 一边亲一边摸免费视频|