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

    芒果轉(zhuǎn)錄組中SSR位點信息分析與引物篩選

    2015-05-30 05:13:09羅純武紅霞姚全勝王松標許文天馬小衛(wèi)
    熱帶作物學報 2015年7期
    關(guān)鍵詞:轉(zhuǎn)錄組芒果多態(tài)性

    羅純 武紅霞 姚全勝 王松標 許文天 馬小衛(wèi)

    摘 要 利用MicroSAtellite軟件分析篩選芒果轉(zhuǎn)錄組中的SSR位點。結(jié)果表明,在54 207條Unigene中共搜索得到4 103個SSR位點,出現(xiàn)頻率為7.57%。其中,三核苷酸重復(fù)為主導(dǎo)重復(fù)類型,占SSR總數(shù)的38.19%,其次是單核苷酸重復(fù)(26.91%)。AG/CT和AAG/CTT分別是二核苷酸和三核苷酸重復(fù)的優(yōu)勢基元。根據(jù)SSR側(cè)翼序列共設(shè)計6 915對SSR引物,隨機挑選了230對進行PCR檢測,獲得93對多態(tài)性引物。由此可見,芒果轉(zhuǎn)錄組數(shù)據(jù)可以作為大量開發(fā)SSR標記的資源,這些SSR標記也將有助于芒果遺傳多樣性和種質(zhì)資源鑒定的研究。

    關(guān)鍵詞 芒果;轉(zhuǎn)錄組;SSR信息;多態(tài)性

    中圖分類號 S667.7 文獻標識碼 A

    Abstract A MicroSAtellite software was used to find the SSR loci in the mango transcriptome. As a result, a total of 4 103 SSRs were found from 54 207 unigenes, with a frequency of 7.57%. The trinucleotide repeats were the most frequent SSRs among them(38.19%). AG/CT and AAG/CTT were the most common motifs in the dinucleotide and trinucleotide repeats, respectively. Based on the flanking sequences of these SSR loci, 6 915 SSR primer pairs were designed. A total of 230 primer pairs were randomly selected to be verified by PCR, and 93 pairs(40.43%)of them exhibited to be polymorphism. The results showed that transcriptome data would be an effective source for the exploration of SSR markers for mango and these SSR markers would contribute greatly to the research into genetic diversity and germplasm characterization in mango.

    Key words Mango; Transcriptome; SSR information; Polymorphism

    doi 10.3969/j.issn.1000-2561.2015.07.013

    Simple sequence repeat(SSR,簡單重復(fù)序列)是指以1~6個核苷酸為單位多次串聯(lián)重復(fù)的DNA序列,在人類和其他動植物基因組中普遍存在且隨機分布[1]。SSR重復(fù)單位的數(shù)目高度變異,但是某一特定SSR序列的側(cè)翼序列通常為高度保守的單一序列,根據(jù)SSR側(cè)翼序列設(shè)計引物通過PCR技術(shù)即可分析SSR的多態(tài)性。與其它分子標記技術(shù)相比,SSR標記具有多態(tài)性高、共顯性、易用PCR檢測、重復(fù)性高、數(shù)量豐富和對基因組有很好的覆蓋性等特點[2]?;谝陨线@些特點,SSR已成為在生物上應(yīng)用廣泛和重要的DNA分子標記。在植物遺傳育種中,SSR標記已經(jīng)在植物遺傳多樣性研究、品種鑒定、分子標記輔助選擇育種(MAS)、遺傳連鎖圖譜的構(gòu)建和外源染色體片段或目的基因鑒定和定位等方面得到應(yīng)用。根據(jù)SSR的來源可將其分為基因組SSR和EST-SSR。傳統(tǒng)的基因組SSR標記引物設(shè)計過程較慢,需要構(gòu)建基因組DNA文庫,探針雜交,重復(fù)序列克隆的識別和篩選,測序等許多繁瑣的過程,不但耗時費力,而且開發(fā)成本也比較高,效率較低[3]。另外,利用該方法所獲得的SSR還會受開發(fā)過程所采用探針的SSR重復(fù)單元種類限制[4]。而EST-SSR分子標記與基因組SSR標記相比,在不同的物種之間具有較好通用性,其開發(fā)方法較簡單而且成本又較低。

    芒果(Mangifera indica L.)是世界五大水果之一,且因其風味獨特、營養(yǎng)豐富,經(jīng)濟效益高等優(yōu)點,種植面積不斷擴大,產(chǎn)量也逐年上升。中國是世界芒果主要生產(chǎn)國之一,芒果種植在國民經(jīng)濟上具有重要作用。目前,國內(nèi)外學者利用多種DNA分子標記技術(shù)如RAPD[5-7]、ISSR[8]、AFLP[9-10]等對芒果進行了品種或真實雜種鑒定,并且對野生、砧木或栽培品種等資源進行了親緣關(guān)系及遺傳多樣性分析。芒果SSR標記也有開發(fā)[11-13],但是目前基于EST數(shù)據(jù)開發(fā)的SSR標記非常少,限制了芒果分子水平的研究發(fā)展。因此,本研究對芒果轉(zhuǎn)錄組測序產(chǎn)生的海量EST序列進行全轉(zhuǎn)錄組SSR鑒定及引物設(shè)計,揭示芒果轉(zhuǎn)錄組SSR的分布規(guī)律和特性,為芒果遺傳圖譜構(gòu)建、系譜分析和MAS育種提供高效的分子標記資源。

    1 材料與方法

    1.1 材料

    以‘Zill4個發(fā)育時期的果皮和果肉混合樣品為材料進行轉(zhuǎn)錄組測序,共獲得54 207條Unigene[14],總長度為45 425 223 bp。

    1.2 方法

    1.2.1 芒果轉(zhuǎn)錄組SSR的篩選 利用SSR軟件 MicroSAtellite(MISA,http://pgrc.ipk -gatersleben.de/misa/)對‘Zill轉(zhuǎn)錄組Unigene序列進行SSR查找,具體篩選標準為:SSR重復(fù)基元的長度為1~6 bp,單核苷酸重復(fù)次數(shù)≥12,二核苷酸重復(fù)次數(shù)≥6,三、四核苷酸重復(fù)次數(shù)≥5,五、六核苷酸重復(fù)次數(shù)≥4。只保留SSR重復(fù)單元在Unigene上前后序列均不小于150 bp的SSR,使用Primer3軟件設(shè)計引物。

    1.2.2 芒果SSR-PCR擴增及產(chǎn)物檢測 PCR反應(yīng)體系20 μL,其中包括:1×buffer,40 ng DNA模板,1.8 mmol/L MgCl2,0.4 mmol/L dNTPs,0.2 μmol/L引物,0.75 U Taq酶。擴增反應(yīng)是在TaKaRa PCR Thermal Cycler Dice儀器上進行,擴增程序如下:94 ℃預(yù)變性2.5 min;94 ℃變性30 s,58 ℃復(fù)性45 s,72 ℃延伸1 min,35個循環(huán);72 ℃延伸5 min。擴增產(chǎn)物在6%變性聚丙烯酰胺凝膠上電泳分離。

    2 結(jié)果與分析

    2.1 芒果轉(zhuǎn)錄組中SSR位點的數(shù)量與分布

    利用軟件MISA對芒果轉(zhuǎn)錄組54 207條Unigene進行分析查找,共3 683條Unigene含有SSR,發(fā)生頻率(含有SSR的Unigene數(shù)目/總Unigene數(shù)目)為6.79%。其中,3 301條Unigene含有單個SSR位點,占含有SSR序列條數(shù)的89.63%,其余382條Unigene含有多個SSR位點,占含有SSR序列條數(shù)的10.37%。芒果轉(zhuǎn)錄組中共查找到SSR位點4 103個,出現(xiàn)頻率為(檢出的SSR個數(shù)/總Unigene數(shù)目)7.57%,平均距離(芒果轉(zhuǎn)錄組Unigene總長度/SSR數(shù)目)為11.07 kb,即芒果轉(zhuǎn)錄組中平均每11.07 kb出現(xiàn)1個SSR(表1)。

    芒果轉(zhuǎn)錄組SSR數(shù)量較為豐富,單核苷酸重復(fù)至六核苷酸重復(fù)均有分布,但是類型SSR出現(xiàn)頻率存在差異(表1),主要是單、二和三核苷酸重復(fù)類型。三核苷酸重復(fù)類型數(shù)量最多,占全部SSR的38.19%;其次是單核苷酸和二核苷酸重復(fù)類型,分別占26.91%和25.83%,四核苷酸重復(fù)類型數(shù)量最少,占2.10%。

    芒果轉(zhuǎn)錄組中SSR位點重復(fù)次數(shù)以5次(1 017)最多,占總SSR的24.79%;其次為重復(fù)次數(shù)6、7、4、8、9和10次,SSR位點個數(shù)在97~859。統(tǒng)計4~10次重復(fù)的SSR位點共有2 932個,占71.46%;11~23次重復(fù)的SSR位點共有1 171個,占總SSR數(shù)量的28.54%(圖1)。

    2.2 芒果轉(zhuǎn)錄組中SSR特性

    芒果轉(zhuǎn)錄組中共查找到153種重復(fù)基元,其中以單核苷酸重復(fù)基元A/T最多,占總SSR位點的26.81%,其次是二核苷酸重復(fù)基元AG/CT和三核苷酸重復(fù)基元AAG/CTT,分別占總SSR數(shù)目的19.79%和13.21%。具體來說,二核苷酸重復(fù)基元共查找到3種,AG/CT占主導(dǎo)地位,占該重復(fù)類型的76.60%,其它2種AC/GT和AT/AT數(shù)量相當;三核苷酸重復(fù)基元共10種,出現(xiàn)次數(shù)最多的是AAG/CTT,占該重復(fù)類型的34.59%,其次是ATC/ATG和ACC/GGT,分別占16.46%和13.08%,出現(xiàn)次數(shù)最少的是ACT/AGT,僅占該重復(fù)類型的1.21%;四、五和六核苷酸重復(fù)基元最多的分別是AAAG/CTTT、AAGAG/CTCTT和ACCCTG/AGGGTC,其它重復(fù)基元出現(xiàn)頻率較低(表2)。

    2.3 芒果轉(zhuǎn)錄組SSR的可用性評價

    SSR多態(tài)性的高低是評價其有用性的主要標準,而SSR的長度影響其多態(tài)性。當SSR長度在20 bp以上時,在不同品種間顯示出較高的多態(tài)性,長度在12~20 bp的SSR多態(tài)性低些,而長度在12 bp以下的SSR多態(tài)性很低[15],因此本研究在查找SSR時將長度小于12 bp的已經(jīng)除去。芒果轉(zhuǎn)錄組SSR的長度為12~60 bp,長度在12~20 bp的SSR數(shù)量最多,共3 532個,占全部SSR的86.08%。其次是21~30 bp的SSR,共561個(13.67%)(圖2)。依據(jù)高級基元SSR多態(tài)性普遍比低級基元的低[16],而在長度大于20 bp的芒果SSR中,主要是低級基元SSR,包含單、二和三核苷酸SSR,共378個,可以預(yù)計這部分多態(tài)性潛能高的SSR在芒果上應(yīng)具有較高的利用價值。

    2.4 芒果轉(zhuǎn)錄組SSR引物設(shè)計與篩選

    使用Primer 3軟件對上下游序列均不小于150 bp的SSR設(shè)計引物,每條序列產(chǎn)生5對引物,共有1 383條芒果SSR序列獲得6 915對引物。為了驗證這些引物的有效性和多態(tài)性,隨機挑選230對引物對7個不同的芒果品種進行PCR擴增(圖3)。結(jié)果表明除12對引物未獲得目的產(chǎn)物外,其余218對引物均獲得穩(wěn)定和可重復(fù)的目的產(chǎn)物,表現(xiàn)多態(tài)性的引物共93對,多態(tài)性比例為40.43%(表3)。

    3 討論與結(jié)論

    本研究在芒果轉(zhuǎn)錄組54 207條Unigene序列中搜索到SSR位點4 103個,分布于3 683條Unigene上,出現(xiàn)頻率為7.57%,平均距離為11.07 kb。與其他植物相比,芒果轉(zhuǎn)錄組SSR出現(xiàn)頻率低于枇杷(22.02%)、油棕(22.6%)和橡膠樹(35.58%)[17-19],但高于大麥(2.8%)、棉花(5.29%)和小麥(7.4%)[20-22],與蓮花(7.6%)和木豆(7.6%)相當[23-24]。SSR的出現(xiàn)頻率在不同物種、不同研究結(jié)果中存在差異主要是搜索重復(fù)類型和長度等標準不同以及分析數(shù)據(jù)量的大小不同所造成的[25]。

    已報道的大多數(shù)植物中SSR以二核苷酸和三核苷酸重復(fù)類型最為豐富,與芝麻、獼猴桃、咖啡、橡膠樹等二核苷酸SSR為主導(dǎo)的植物不同[26-29],芒果轉(zhuǎn)錄組中三核苷酸SSR數(shù)量最多,占全部SSR的38.19%(表1),與甘蔗、小麥、紅薯和柑橘等植物類似[30-33]。芒果中二核苷酸優(yōu)勢重復(fù)基元為AG/CT,與多數(shù)植物中報道的情況相同[26,29]。前人研究發(fā)現(xiàn)雙子葉植物中三核苷酸重復(fù)以AAG/CTT重復(fù)為主,CCG/CGG出現(xiàn)頻率很低[34-35],芒果三核苷酸SSR中AAG/CTT數(shù)量最多,占總SSR的13.21%,而CCG/CGG僅占0.95%,驗證了前人研究。

    在芒果上,國內(nèi)外已報道的SSR標記數(shù)量并不多,主要是通過構(gòu)建基因組文庫,利用探針雜交富集篩選得到,如Duval等[36]、Viruel等[37]、Honsho等[38]和Chiang等[39]通過該方法分別開發(fā)了28、16、6和20個SSR標記。國內(nèi)黃啟星等[40]利用公共數(shù)據(jù)庫中1 665條芒果EST序列分析查找到EST-SSR位點32個,EST-SSR平均距離為20.2 kb,大于本研究中利用轉(zhuǎn)錄組數(shù)據(jù)查找的SSR的平均距離,可能是由于轉(zhuǎn)錄組測序數(shù)據(jù)較豐富,從而含有SSR位點也更豐富。同時黃啟星等[40]還發(fā)現(xiàn)芒果EST-SSR以三核苷酸重復(fù)類型數(shù)量最多,其中又以AAG/CTT重復(fù)為主,與本研究得到的結(jié)果一致。

    總的來說,芒果轉(zhuǎn)錄組中SSR的出現(xiàn)頻率雖然不高,但是類型比較豐富,長度大于20 bp的主要是低基元重復(fù),多態(tài)性潛能高,設(shè)計的引物經(jīng)實驗驗證其有效性也高。本研究結(jié)果為進一步構(gòu)建芒果遺傳圖譜、系譜分析和MAS育種提供了高效的分子標記資源。

    參考文獻

    [1] Tautz D, Renz M. Simple sequence repeats are ubiquitous repetitive components of eukaryotic genomes[J]. Nucl Acids Res, 1984, 12(10): 4 127-4 138.

    [2] Powell W, Machray G C, Provan J. Polymorphism revealed by simple sequence repeat[J]. Trends Plant Sci, 1996, 1(7): 215-222.

    [3] Roder M S, Korzun V, Wendehake K, et al. A microsatellite map of wheat[J]. Genetics, 1998, 149(4): 2 007-2 023.

    [4] Chen C, Zhou P, Choi Y A, et al. Mining and characterizing microsatellites from citrus ESTs[J]. Theor Appl Genet, 2006, 112(7): 1 248-1 257.

    [5] Karihaloo J L, Dwivedi Y K, Sunil A, et al. Analysis of genetic diversity of Indian mango cultivars using RAPD markers[J]. J Hortic Sci Biotech, 2003, 78: 285-289.

    [6] Shukla M, Babu R, Mathur V K, et al. Diverse genetic bases of Indian polyembryonic and monoembryonic mango(Mangifera indica L)cultivars[J]. Curr Sci, 2004, 87: 870-871.

    [7] Rahman M L, Rabbani M G, Siddique M N A, et al. Molecular characterization of 28 mango germplasm using RAPD[J]. Plant Tiss Culture Biotech, 2007, 17(1): 71-77.

    [8] Rocha A, Carlos L, Salomao C, et al. Genetic diversity of ‘Ubá mango tree using ISSR markers[J]. Mol Biotechnol, 2012, 50(2): 108-113.

    [9] 房經(jīng)貴, 章 鎮(zhèn), 馬正強, 等. AFLP標記在兩個芒果品種間雜交F1代的多態(tài)性及分離方式[J]. 中國農(nóng)業(yè)科學, 2000, 33(3): 19-24.

    [10] Kashkush K, Fang J G, Tomer E, et al. Cultivar identification and genetic map of mango(Mangifera indica)[J]. Euphytica, 2001, 122: 129-136.

    [11] Hirano R, Ishii H, Oo T H, et al. Propagation management methods have altered the genetic variability of two traditional mango varieties in Myanmar, as revealed by SSR[J]. Plant Genet Resour, 2011, 9(3): 404-410.

    [12] Tsai C C, Chen Y H, Chen C H, et al. Cultivar identification and genetic relationship of mango(Mangifera indica)in Taiwan using 37 SSR markers[J]. Sci Hortic, 2013, 164: 196-201.

    [13] Surapaneni M, Vemireddy L R, Begum H, et al. Population structure and genetic analysis of different utility types of mango(Mangifera indica L.)germplasm of Andhra Pradesh state of India using microsatellite markers[J]. Plant Syst Evol, 2013, 299(7): 1 215-1 229.

    [14] Wu H X, Jia H M, Ma X W, et al. Transcriptome and proteomic analysis of mango(Mangifera indica Linn)fruits[J]. J Proteomics, 2014, 105: 19-30.

    [15] Temnykh S, De Clerck G, Lukashova L, et al. Computational andexperimental analysis of microsatellites in rice(Oryzasativa L.)[J]. Genome Res, 2001, 11: 1 441-1 452.

    [16] Dreisigacker S, Zhang P, Warburton M L, et al. SSR and pedigree analyses of genetic diversity among CIMMYT wheat lines targeted to different mega environments[J]. Crop Sci, 2004, 44: 381-388.

    [17] Li X Y, Xu H X, Feng J J, et al. Development and application of genic simple sequence repeat markers from the transcriptome of loquat[J]. J Amer Soc Hort Sci, 2014, 139(5): 507-517.

    [18] 周麗霞, 肖 勇, 楊耀東. 油棕轉(zhuǎn)錄組SSR標記開發(fā)研究[J]. 廣東農(nóng)業(yè)科學, 2014(14): 136-138.

    [19] Mantello C C, Cardoso-Silva C B, da Silva C C, et al. De Novo assembly and transcriptome analysis of the rubber tree (Hevea brasiliensis)and SNP markers development for rubber biosynthesis pathways[J]. PLoS ONE, 2014, 9: e102 665.

    [20] Varshney R, Grosse I, Hahnel U, et al. Genetic mapping and BAC assignment of EST-derived SSR markers shows non-uniform distribution of genes in the barley genome[J]. Theor Appl Genet, 2006, 113(2): 239-250.

    [21] Zhang X W, Ye Z W a, Wang T K, et al. Characterization of the global transcriptome for cotton(Gossypium hirsutum L.)anther and development of SSR marker[J]. Gene, 2014, 551(2): 206-213.

    [22] Peng J, Lapitan N L V. Characterization of EST-derived microsatellites in the wheat genome and development of eSSR markers[J]. Funct Integr Genomic, 2005, 5(2): 80-96.

    [23] Zhang W, Tian D, Huang X, et al. Characterization of flower-bud transcriptome and development of genic SSR markers in asian lotus(Nelumbo nucifera Gaertn.)[J]. PLoS ONE, 2014, 9(11): e112 223.

    [24] Dutta S, Kumawat G, Singh B P, et al. Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan(L.)Millspaugh][J]. BMC Plant Biol, 2011, 11: 17.

    [25] Varshney R K, Graner A, Sorrells M E. Genic microsatellite markers in plants: features and applications[J]. Trends Biotechnol, 2005, 23(1): 48-55.

    [26] Wei W, Qi X, Wang L, et al. Characterization of the sesame (Sesamum indicum L.)global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers[J]. BMC Genomics, 2011, 12: 451.

    [27] Fraser L G, Harvey C F, Crowhurst R N, et al. EST-derived microsatellites from Actinidia species and their potential for mapping[J]. Theor Appl Genet, 2004, 108(6): 1 010-1 016.

    [28] Aggarwal R K, Hendre P S, Varshney R K, et al. Identification, characterization and utilization of EST-derived genic microsatellite markers for genome analyses of coffee and related species[J]. Theor Appl Genet, 2007, 114(2): 359-372.

    [29] Li D, Deng Z, Qin B, et al. De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree(Hevea brasiliensis Muell. Arg.)[J]. BMC Genomics, 2012, 13: 192.

    [30] Marconi T G, Costa E, Miranda H R, et al. Functional markers for gene mapping and genetic diversity studies in sugarcane[J]. BMC Res Notes, 2011, 4: 264.

    [31] La Rota M, Kantety R V, Yu J K, et al. Nonrandom distribution and frequencies of genomic and EST-derived microsatellite markers in rice, wheat, and barley[J]. BMC Genomics, 2005, 6: 23.

    [32] Wang Z, Fang B, Chen J, et al. De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweet potato(Ipomoea batatas)[J]. BMC Genomics, 2010, 11: 726.

    [33] Chen C, Zhou P, Choi Y, et al. Mining and characterizing microsatellites from citrus ESTs[J]. Theor Appl Genet, 2006, 112(7): 1 248-1 257.

    [34] Kumpatla S, Mukhopadhyay S. Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species[J]. Genome, 2005, 48(6): 985-998.

    [35] Liang X, Chen X, Hong Y, et al. Utility of EST-derived SSR in cultivated peanut(Arachis hypogaea L.)and Arachis wild species[J]. BMC Plant Biol, 2009, 9: 35.

    [36] Duval M F, Bunel J, Sitbon C, et al. Development of microsatellite markers for mango(Mangifera indica L.)[J]. Mol Ecol Notes, 2005, 5(4): 824-826.

    [37] Viruel M A, Escribano P, Barbieri M, et al. Fingerprinting, embryo type and geographic differentiation in mango(Mangifera indica L., Anacardiaceae)with microsatellites[J]. Mol Breeding, 2005, 15(4): 383-393.

    [38] Honsho C, Nishiyama K, Eiadthong W, et al. Isolation and characterization of new microsatellite markers in mango(Mangifera indica)[J]. Mol Ecol Notes, 2005, 5(1): 152-154.

    [39] Chiang Y C, Tsai C M, Chen Y K H, et al. Development and characterization of 20 new polymorphic microsatellite markers from Mangifera indica(Anacardiaceae)[J]. Am J Bot, 2012, 99(3): e117-e119.

    [40] 黃啟星, 左 嬌, 孔 華, 等. 11種熱帶植物EST-SSR標記的開發(fā)和多樣性分析[J]. 熱帶作物學報, 2012, 33(7): 1 208-1 214.

    猜你喜歡
    轉(zhuǎn)錄組芒果多態(tài)性
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進展
    我才不要穿
    小洞會“咬”人
    小洞會“咬”人
    基于轉(zhuǎn)錄組測序的山茱萸次生代謝生物合成相關(guān)基因的挖掘
    金釵石斛轉(zhuǎn)錄組SSR位點信息分析
    人參屬藥用植物轉(zhuǎn)錄組研究進展
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    GlobalFiler~? PCR擴增試劑盒驗證及其STR遺傳多態(tài)性
    大豆轉(zhuǎn)錄組測序研究進展綜述
    国产成人免费观看mmmm| 日日摸夜夜添夜夜添av毛片| 三上悠亚av全集在线观看 | 麻豆精品久久久久久蜜桃| 熟女人妻精品中文字幕| 日韩免费高清中文字幕av| 亚洲欧美精品专区久久| 成年美女黄网站色视频大全免费 | 一级毛片久久久久久久久女| www.色视频.com| 69精品国产乱码久久久| 91aial.com中文字幕在线观看| 亚洲美女黄色视频免费看| 亚洲欧美成人综合另类久久久| www.色视频.com| 黑人猛操日本美女一级片| 人人妻人人添人人爽欧美一区卜| 久久韩国三级中文字幕| 久久国内精品自在自线图片| 黄色视频在线播放观看不卡| 久久99热6这里只有精品| 亚洲av在线观看美女高潮| 国产成人精品一,二区| 91精品一卡2卡3卡4卡| 国产av码专区亚洲av| 国产一区二区在线观看av| 黄色毛片三级朝国网站 | 一区二区av电影网| 精品国产国语对白av| 国产91av在线免费观看| h日本视频在线播放| 免费播放大片免费观看视频在线观看| 伦精品一区二区三区| 成人亚洲精品一区在线观看| 午夜福利网站1000一区二区三区| 国产精品三级大全| 97精品久久久久久久久久精品| 十八禁高潮呻吟视频 | 亚洲国产精品成人久久小说| 日韩视频在线欧美| 亚洲激情五月婷婷啪啪| 久久午夜综合久久蜜桃| 亚洲精品中文字幕在线视频 | 最近最新中文字幕免费大全7| 国产免费一区二区三区四区乱码| 一级爰片在线观看| 22中文网久久字幕| av天堂久久9| 99热6这里只有精品| 两个人免费观看高清视频 | 国产黄色视频一区二区在线观看| 秋霞伦理黄片| 交换朋友夫妻互换小说| 男女边吃奶边做爰视频| 国产黄色视频一区二区在线观看| 日韩成人伦理影院| 中文欧美无线码| 91久久精品国产一区二区三区| 国产美女午夜福利| 日韩欧美一区视频在线观看 | 免费播放大片免费观看视频在线观看| 成人亚洲欧美一区二区av| 欧美xxⅹ黑人| 国产精品国产三级专区第一集| 亚洲欧美精品专区久久| 桃花免费在线播放| 国产精品秋霞免费鲁丝片| 少妇人妻一区二区三区视频| 日本黄大片高清| av天堂中文字幕网| 国产探花极品一区二区| 国产伦精品一区二区三区四那| 国产亚洲av片在线观看秒播厂| 免费看av在线观看网站| 最近的中文字幕免费完整| 国产黄片视频在线免费观看| √禁漫天堂资源中文www| 亚洲怡红院男人天堂| 中文欧美无线码| 婷婷色麻豆天堂久久| 国产亚洲5aaaaa淫片| 久久99一区二区三区| 人体艺术视频欧美日本| 久久女婷五月综合色啪小说| 99re6热这里在线精品视频| 在线观看人妻少妇| 高清毛片免费看| 久久99热6这里只有精品| 狂野欧美白嫩少妇大欣赏| 亚洲综合精品二区| 老司机影院成人| 国产在视频线精品| 99热6这里只有精品| 综合色丁香网| 国产精品人妻久久久影院| 日本vs欧美在线观看视频 | 久久ye,这里只有精品| 中文精品一卡2卡3卡4更新| 国产亚洲91精品色在线| 亚洲av.av天堂| 国产一区二区在线观看av| 国产一区二区三区av在线| 日韩av在线免费看完整版不卡| av播播在线观看一区| 日韩中字成人| 欧美bdsm另类| 久久精品久久久久久噜噜老黄| 菩萨蛮人人尽说江南好唐韦庄| 少妇人妻久久综合中文| 久久久久久久久大av| 精品久久久久久电影网| 国内揄拍国产精品人妻在线| 哪个播放器可以免费观看大片| 2022亚洲国产成人精品| 国产成人精品无人区| 国产一区二区三区综合在线观看 | 综合色丁香网| 人人澡人人妻人| 日日摸夜夜添夜夜爱| 在线观看美女被高潮喷水网站| 精品一区在线观看国产| 黑人高潮一二区| 一级毛片 在线播放| 免费大片18禁| 99久国产av精品国产电影| 国产亚洲一区二区精品| 中文在线观看免费www的网站| a级一级毛片免费在线观看| 两个人免费观看高清视频 | 国产一区二区在线观看av| 女人精品久久久久毛片| av天堂久久9| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 精品熟女少妇av免费看| 黑丝袜美女国产一区| 寂寞人妻少妇视频99o| 亚洲精品自拍成人| 欧美日韩在线观看h| 国产精品秋霞免费鲁丝片| av一本久久久久| 亚洲av电影在线观看一区二区三区| 噜噜噜噜噜久久久久久91| 青青草视频在线视频观看| 免费观看性生交大片5| 免费大片18禁| 色视频在线一区二区三区| 久久精品久久精品一区二区三区| 久久韩国三级中文字幕| 热re99久久国产66热| 亚洲国产最新在线播放| 黄色怎么调成土黄色| 国产成人a∨麻豆精品| 日韩亚洲欧美综合| 少妇熟女欧美另类| 啦啦啦啦在线视频资源| 欧美成人午夜免费资源| 少妇精品久久久久久久| 国产免费又黄又爽又色| 亚洲国产毛片av蜜桃av| 久热这里只有精品99| 免费看不卡的av| 久久鲁丝午夜福利片| 国产男女超爽视频在线观看| 日本免费在线观看一区| 久久热精品热| 街头女战士在线观看网站| 精华霜和精华液先用哪个| 爱豆传媒免费全集在线观看| 2018国产大陆天天弄谢| 亚洲av电影在线观看一区二区三区| 亚洲,一卡二卡三卡| 午夜激情福利司机影院| 老女人水多毛片| 在线观看免费视频网站a站| 亚洲性久久影院| 99精国产麻豆久久婷婷| 两个人的视频大全免费| 国产精品福利在线免费观看| 99久久精品热视频| 成人影院久久| 免费观看无遮挡的男女| 国产中年淑女户外野战色| 三上悠亚av全集在线观看 | 久久精品国产自在天天线| 国产精品免费大片| 啦啦啦视频在线资源免费观看| 久久久久网色| 国产亚洲一区二区精品| 老女人水多毛片| 大香蕉久久网| 国产极品天堂在线| 国产欧美日韩综合在线一区二区 | 亚洲欧美中文字幕日韩二区| 乱人伦中国视频| 蜜桃在线观看..| 久久久久久久久久人人人人人人| 寂寞人妻少妇视频99o| 日本午夜av视频| 青青草视频在线视频观看| 韩国高清视频一区二区三区| 精品一区二区三区视频在线| 黄色视频在线播放观看不卡| 免费看日本二区| 亚洲三级黄色毛片| 国产探花极品一区二区| 亚洲欧美日韩东京热| 国产男女超爽视频在线观看| 亚洲欧美精品自产自拍| 午夜福利,免费看| 高清视频免费观看一区二区| 99re6热这里在线精品视频| 天天躁夜夜躁狠狠久久av| 9色porny在线观看| 久久精品国产亚洲网站| av专区在线播放| 国产精品久久久久久久电影| 亚洲自偷自拍三级| 国产 精品1| 久久国产精品男人的天堂亚洲 | 男人爽女人下面视频在线观看| 韩国高清视频一区二区三区| 国产精品人妻久久久久久| 日韩熟女老妇一区二区性免费视频| 婷婷色综合大香蕉| 久久久国产一区二区| 只有这里有精品99| 亚洲内射少妇av| 2021少妇久久久久久久久久久| 麻豆乱淫一区二区| 欧美成人午夜免费资源| 成年人免费黄色播放视频 | 欧美bdsm另类| 午夜久久久在线观看| 免费观看无遮挡的男女| 精品一区二区三区视频在线| 久久精品国产亚洲网站| 91精品国产九色| av天堂久久9| 国产视频首页在线观看| 97超视频在线观看视频| 国产精品欧美亚洲77777| 午夜福利视频精品| 日本欧美国产在线视频| 国产亚洲一区二区精品| 伊人久久国产一区二区| 少妇被粗大的猛进出69影院 | 国内少妇人妻偷人精品xxx网站| 亚洲美女视频黄频| 少妇的逼好多水| 美女xxoo啪啪120秒动态图| 性色avwww在线观看| 秋霞在线观看毛片| 妹子高潮喷水视频| 新久久久久国产一级毛片| h视频一区二区三区| 亚洲成色77777| 伊人亚洲综合成人网| 久久99精品国语久久久| 日韩制服骚丝袜av| 国产亚洲午夜精品一区二区久久| 在线观看美女被高潮喷水网站| 免费黄色在线免费观看| 97超碰精品成人国产| 国产熟女欧美一区二区| 亚洲欧洲国产日韩| 18禁在线播放成人免费| 国产精品一区www在线观看| 97在线人人人人妻| 国产在线视频一区二区| 日日啪夜夜撸| 高清av免费在线| 国产在视频线精品| videos熟女内射| 久久鲁丝午夜福利片| 国产黄色免费在线视频| 七月丁香在线播放| 观看美女的网站| 国产亚洲精品久久久com| 成人二区视频| 99re6热这里在线精品视频| 哪个播放器可以免费观看大片| 爱豆传媒免费全集在线观看| 国产成人a∨麻豆精品| 欧美日韩亚洲高清精品| 高清视频免费观看一区二区| 伊人亚洲综合成人网| 国产精品秋霞免费鲁丝片| 91在线精品国自产拍蜜月| 国产片特级美女逼逼视频| 久久久久久人妻| 欧美变态另类bdsm刘玥| 国产成人精品无人区| 国产伦精品一区二区三区视频9| 久久人人爽av亚洲精品天堂| 亚洲欧美中文字幕日韩二区| 美女中出高潮动态图| 亚洲欧美精品自产自拍| 又粗又硬又长又爽又黄的视频| 99久国产av精品国产电影| 亚洲精品乱码久久久v下载方式| 岛国毛片在线播放| 亚洲欧洲国产日韩| 日韩成人av中文字幕在线观看| 丝袜脚勾引网站| 两个人免费观看高清视频 | 精品久久久精品久久久| 欧美一级a爱片免费观看看| 国产女主播在线喷水免费视频网站| 夫妻午夜视频| 夜夜骑夜夜射夜夜干| 熟妇人妻不卡中文字幕| 另类亚洲欧美激情| 日本黄大片高清| 日韩不卡一区二区三区视频在线| 丰满乱子伦码专区| 麻豆成人午夜福利视频| 久久久久久久久久久丰满| 亚洲在久久综合| 国内少妇人妻偷人精品xxx网站| 美女国产视频在线观看| 天美传媒精品一区二区| av女优亚洲男人天堂| 日本91视频免费播放| 在线精品无人区一区二区三| 高清黄色对白视频在线免费看 | av女优亚洲男人天堂| 丁香六月天网| 久久女婷五月综合色啪小说| 午夜av观看不卡| 欧美性感艳星| 最近的中文字幕免费完整| 国产色婷婷99| av福利片在线观看| 日韩中字成人| 18禁裸乳无遮挡动漫免费视频| 久久久精品94久久精品| 精品卡一卡二卡四卡免费| 欧美日韩av久久| 日本91视频免费播放| .国产精品久久| 日本wwww免费看| 国产高清有码在线观看视频| 在线观看www视频免费| 97超碰精品成人国产| 最新的欧美精品一区二区| 亚洲精品国产成人久久av| 亚洲欧美清纯卡通| 亚洲精品乱码久久久v下载方式| 黄色一级大片看看| 国产精品偷伦视频观看了| 欧美国产精品一级二级三级 | 99国产精品免费福利视频| a级毛片在线看网站| 国产日韩欧美在线精品| 亚洲成色77777| 精华霜和精华液先用哪个| 偷拍熟女少妇极品色| 99热全是精品| 久久久欧美国产精品| 交换朋友夫妻互换小说| 一边亲一边摸免费视频| 一本大道久久a久久精品| 免费黄网站久久成人精品| 国产探花极品一区二区| 成人毛片a级毛片在线播放| 国产精品国产三级国产av玫瑰| 久久 成人 亚洲| 中文乱码字字幕精品一区二区三区| 久久精品国产a三级三级三级| 国产一区二区三区综合在线观看 | 久久国内精品自在自线图片| 丁香六月天网| 夜夜看夜夜爽夜夜摸| 日韩不卡一区二区三区视频在线| 午夜免费鲁丝| 亚洲欧美日韩另类电影网站| 久久狼人影院| 制服丝袜香蕉在线| 亚洲伊人久久精品综合| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩另类电影网站| 青春草视频在线免费观看| 制服丝袜香蕉在线| 亚洲av福利一区| 国产在视频线精品| 18禁动态无遮挡网站| 成年人免费黄色播放视频 | 国产成人精品婷婷| 国产老妇伦熟女老妇高清| 国产一区有黄有色的免费视频| 黄色日韩在线| 一区二区av电影网| 国产成人免费无遮挡视频| 国产爽快片一区二区三区| 少妇高潮的动态图| 最新中文字幕久久久久| 日日爽夜夜爽网站| 欧美 日韩 精品 国产| 伊人久久精品亚洲午夜| 精品国产乱码久久久久久小说| 亚洲欧洲精品一区二区精品久久久 | 日韩精品有码人妻一区| 久久久久久久精品精品| 国产精品一二三区在线看| 在线精品无人区一区二区三| 如日韩欧美国产精品一区二区三区 | 精品人妻一区二区三区麻豆| .国产精品久久| 国内揄拍国产精品人妻在线| 视频中文字幕在线观看| 嘟嘟电影网在线观看| 日本黄色日本黄色录像| 人人妻人人澡人人看| 一本久久精品| 午夜影院在线不卡| 哪个播放器可以免费观看大片| 美女内射精品一级片tv| 亚洲精品一区蜜桃| 黄色一级大片看看| 99久久精品热视频| 六月丁香七月| 乱码一卡2卡4卡精品| 三级国产精品欧美在线观看| xxx大片免费视频| 国产精品嫩草影院av在线观看| 久久婷婷青草| 日韩av在线免费看完整版不卡| 精品一区在线观看国产| 99热6这里只有精品| 久久av网站| av免费观看日本| 在线观看免费高清a一片| 高清毛片免费看| 寂寞人妻少妇视频99o| 日韩av免费高清视频| 久久人人爽人人片av| 亚洲无线观看免费| 最新中文字幕久久久久| 毛片一级片免费看久久久久| 97在线人人人人妻| 午夜福利网站1000一区二区三区| 一级毛片电影观看| 麻豆精品久久久久久蜜桃| 黄色一级大片看看| 爱豆传媒免费全集在线观看| 国产 精品1| 亚洲性久久影院| 国产一区亚洲一区在线观看| 日韩一区二区三区影片| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 国产精品一区www在线观看| 蜜桃在线观看..| 亚洲精品日韩av片在线观看| 一区二区三区免费毛片| 日韩视频在线欧美| 国产一区二区三区av在线| 欧美日韩国产mv在线观看视频| 99热这里只有精品一区| 久久鲁丝午夜福利片| 亚洲一级一片aⅴ在线观看| 韩国av在线不卡| 久久女婷五月综合色啪小说| 2018国产大陆天天弄谢| 热99国产精品久久久久久7| 午夜福利,免费看| 国产午夜精品久久久久久一区二区三区| 欧美少妇被猛烈插入视频| 丝袜在线中文字幕| 久久久精品免费免费高清| 久久免费观看电影| 黄色毛片三级朝国网站 | 久久久久国产网址| 99热6这里只有精品| 国产精品一区二区在线观看99| 我要看日韩黄色一级片| 日日啪夜夜爽| 亚洲精品久久久久久婷婷小说| av卡一久久| 国产淫语在线视频| 国内揄拍国产精品人妻在线| 国产黄色视频一区二区在线观看| 久久久久久久久久久免费av| 少妇的逼好多水| 极品教师在线视频| 久久热精品热| 国语对白做爰xxxⅹ性视频网站| 91精品一卡2卡3卡4卡| 99热网站在线观看| 丰满饥渴人妻一区二区三| 亚洲中文av在线| 一本色道久久久久久精品综合| 亚洲av综合色区一区| 国产免费福利视频在线观看| 国产免费视频播放在线视频| 久久毛片免费看一区二区三区| 精品久久久久久电影网| 丝袜在线中文字幕| 久久99蜜桃精品久久| 中文乱码字字幕精品一区二区三区| 免费大片黄手机在线观看| 精品久久国产蜜桃| 欧美成人午夜免费资源| 成人毛片60女人毛片免费| 国产成人免费无遮挡视频| a级一级毛片免费在线观看| 久久久久久久久大av| 麻豆精品久久久久久蜜桃| 国产女主播在线喷水免费视频网站| 插阴视频在线观看视频| 不卡视频在线观看欧美| 国产高清国产精品国产三级| 国产色婷婷99| 桃花免费在线播放| 欧美日韩亚洲高清精品| 亚洲精品视频女| 欧美 亚洲 国产 日韩一| 亚洲经典国产精华液单| 在线观看美女被高潮喷水网站| 一级毛片aaaaaa免费看小| 色吧在线观看| 免费黄网站久久成人精品| 亚洲精品,欧美精品| 国产精品久久久久久精品电影小说| 日韩熟女老妇一区二区性免费视频| 一级毛片久久久久久久久女| 日韩一区二区视频免费看| 一本色道久久久久久精品综合| 免费黄网站久久成人精品| 亚洲激情五月婷婷啪啪| a 毛片基地| 国产免费福利视频在线观看| 国产成人a∨麻豆精品| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 99久久精品一区二区三区| 女人精品久久久久毛片| 高清毛片免费看| 少妇丰满av| 亚洲成人一二三区av| 2022亚洲国产成人精品| 精品人妻熟女av久视频| 国产有黄有色有爽视频| 午夜福利视频精品| 内射极品少妇av片p| 一级毛片电影观看| 国产伦精品一区二区三区四那| 中国三级夫妇交换| 久久精品国产鲁丝片午夜精品| av福利片在线| 欧美日韩亚洲高清精品| 欧美日韩在线观看h| 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 中文乱码字字幕精品一区二区三区| 在线看a的网站| 久久久久精品性色| 男女免费视频国产| 亚洲人与动物交配视频| 婷婷色麻豆天堂久久| 在线观看www视频免费| 国产真实伦视频高清在线观看| 丝袜在线中文字幕| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 亚洲国产精品一区二区三区在线| 一本久久精品| 亚洲av二区三区四区| 日日爽夜夜爽网站| 国产精品99久久久久久久久| 亚洲精品国产色婷婷电影| 人人妻人人看人人澡| 男女国产视频网站| 免费看不卡的av| 搡老乐熟女国产| 99热这里只有是精品50| 久久99精品国语久久久| 99热这里只有是精品在线观看| 男女啪啪激烈高潮av片| 精品国产国语对白av| 99久久精品一区二区三区| 国产亚洲午夜精品一区二区久久| 一级黄片播放器| 永久网站在线| 嫩草影院入口| 尾随美女入室| 国产成人免费观看mmmm| 99热这里只有是精品在线观看| freevideosex欧美| 国产黄色视频一区二区在线观看| 精品一区二区三区视频在线| 你懂的网址亚洲精品在线观看| 美女cb高潮喷水在线观看| 亚洲在久久综合| 人人妻人人添人人爽欧美一区卜| 国产精品麻豆人妻色哟哟久久| 色婷婷av一区二区三区视频| 欧美 亚洲 国产 日韩一| 99视频精品全部免费 在线| 午夜老司机福利剧场| 人人妻人人看人人澡| 男女国产视频网站| 久久久久久久久久久丰满| 草草在线视频免费看| 曰老女人黄片| 男的添女的下面高潮视频| 亚洲av男天堂| 久久人人爽人人片av| 亚洲在久久综合| 美女福利国产在线| 波野结衣二区三区在线| 视频中文字幕在线观看| 伊人久久精品亚洲午夜|