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

    基于SNP標記的504份番茄種質(zhì)資源遺傳多樣性分析

    2022-10-30 08:56:02宋麗娜王曉敏劉文娟付金軍胡新華高艷明李建設(shè)
    核農(nóng)學報 2022年12期
    關(guān)鍵詞:種質(zhì)寧夏多態(tài)性

    宋麗娜 王曉敏,2,3,4,* 劉文娟 楊 宏 付金軍 胡新華 高艷明,2,3,4 李建設(shè),2,3,4

    (1 寧夏大學農(nóng)學院,寧夏 銀川 750021;2 寧夏現(xiàn)代設(shè)施園藝工程技術(shù)研究中心,寧夏 銀川 750021;3 寧夏優(yōu)勢特色作物現(xiàn)代分子育種重點實驗室,寧夏 銀川 750021;4 寧夏設(shè)施園藝(寧夏大學)技術(shù)創(chuàng)新中心,寧夏 銀川 750021;5 寧夏巨豐種苗有限責任公司,寧夏,銀川 750021)

    番茄(Solanumlycopersicom),原產(chǎn)于南美洲,屬于茄科草本植物[1]。番茄果實富含多種營養(yǎng)元素,具有抗氧化能力,能夠預防多種疾病[2]。截止2020年,我國番茄的年產(chǎn)量已達6 515萬噸,種植面積達百萬公頃,目前我國已成為世界上最大的番茄生產(chǎn)國[3]。寧夏地處黃河中上游的氣候干爽,日照充足,為本地番茄栽培創(chuàng)造了有利的環(huán)境條件,尤其在全國大部分地區(qū)遇到持續(xù)高溫天氣時,寧夏地區(qū)因氣候冷涼、晝夜溫差大,種植的番茄果實品質(zhì)佳,使得該地區(qū)成為我國越夏番茄的新產(chǎn)區(qū)[3]。番茄作為寧夏重要農(nóng)產(chǎn)品之一,種質(zhì)資源數(shù)量龐大,因此,開展番茄種質(zhì)資源的遺傳多樣性相關(guān)研究,對番茄品種改良和雜種優(yōu)勢利用,了解物種起源、預測物種適應性、預測遺傳資源分布、開發(fā)與利用種質(zhì)資源等具有重要意義[4-5]。

    迄今為止,已有前人開展了關(guān)于番茄種質(zhì)遺傳多樣性的研究,但是對于種質(zhì)資源遺傳多樣性的評價主要依賴于表型。例如何潤銘等[6]對95份番茄種質(zhì)資源的23個表型進行了遺傳多樣性分析,發(fā)現(xiàn)當前番茄育種在一定程度上趨于同質(zhì)化,應該重視對不同遺傳背景材料的保護;李艷紅等[7]利用表型性狀對華南地區(qū)的櫻桃番茄的種質(zhì)資源進行評價,結(jié)果表明,數(shù)量性狀在不同類型番茄的種質(zhì)個體間差異較大,供試材料的表型遺傳變異豐富。由于不同材料對蕃茄種質(zhì)的多樣性影響有差異,且表型性狀易受基因型和環(huán)境條件的影響,分子標記的發(fā)展為作物遺傳多樣性分析提供了更有力、更精準的工具[8]。目前,已有多種分子標記檢測方法被開發(fā)并應用于番茄育種及遺傳多樣性分析,包括隨機擴增多態(tài)性DNA標記(random amplified polymorphic DNA, RAPD)、簡單重復序列區(qū)間(inter-simple sequence repeat,ISSR)、限制性片段長度多態(tài)性(restriction fragment length polymorphisms,RFLP)、單核苷酸多態(tài)性(single nucleotide polymorphism,SNP)、簡單重復序列(simple sequence repeat,SSR)等。宋建等[9]利用400對SSR標記分析了國內(nèi)外具有代表性的幾大類番茄材料,結(jié)果表明,番茄材料間存在不同程度的遺傳分化,栽培番茄種內(nèi)親緣關(guān)系相對野生種種間親緣關(guān)系較近,遺傳背景十分狹窄,充分反映了我國栽培各類番茄的實際情況。SNP由于具有遺傳穩(wěn)定性和通量高、成本低等優(yōu)點[10],已被廣泛用于國內(nèi)外番茄種質(zhì)資源的遺傳多樣性分析[11]。Sim等[12]利用從整個基因組挖掘出的7 720個SNP標記對收集到的426份番茄遺傳資源進行基因分型,將其劃分為7個不同的番茄類型。馮晶晶等[13]利用59個高質(zhì)量SNP標記、競爭性等位基因特異性PCR(kompetitive allele specific PCR,KASP)技術(shù)對收集到的來源廣泛的433份醋栗番茄種質(zhì)資源進行遺傳多樣性分析,結(jié)果表明,醋栗番茄能夠為栽培番茄提供遺傳變異資源庫。從2015年開始至今,寧夏的蕃茄育種在種質(zhì)資源遺傳多樣性方面開展了一些研究工作[14-16],但是基于分子標記尤其是SNP標記對番茄種質(zhì)資源遺傳多樣性的研究還處于初步發(fā)展階段。

    鑒于此,本研究以60對SNP引物對來自不同地區(qū)504份番茄種質(zhì)資源的遺傳多樣性、聚類分析、主成分分析及群體遺傳結(jié)構(gòu)進行研究,以期揭示不同材料在DNA水平上的遺傳變異,為提高寧夏地區(qū)種質(zhì)資源的利用效率及核心種質(zhì)構(gòu)建提供理論依據(jù)和材料基礎(chǔ)。

    1 材料與方法

    1.1 試驗材料

    本試驗所用504份番茄種質(zhì)資源由寧夏大學農(nóng)學院蔬菜課題組(本課題組)和寧夏巨豐種苗有限公司共同提供,編號為S001~S504,詳細信息見電子附表1。

    1.2 試驗方法

    1.2.1 番茄葉片基因組DNA的提取與檢測 參照芮文婧[17]的方法,采用十六烷基三甲基溴化銨(cetyltrimethylammonium,CTAB)法提取番茄葉片基因組DNA。

    1.2.2 引物篩選與PCR擴增反應 本研究利用寧夏大學瓜菜遺傳育種實驗室篩選的均勻覆蓋12條染色體的60對引物(電子附表2)對504份番茄種質(zhì)資源進行遺傳多樣性分析,PCR擴增反應總體系為5 μL:20~50 ng·μL-1DNA模板2.20 μL、KASP Master Mix 2.50 μL、KBD Assay mix 0.07 μL、ddH2O 0.23 μL。PCR反應程序:95℃預變性15 min;95℃變性20 s,55~61℃延伸60 s,每循環(huán)退火延伸溫度為0.6℃,10個循環(huán);95℃變性20 s,57℃延伸60 s,32個循環(huán)。設(shè)置4個空白對照,模板用ddH2O代替。PCR反應在Pro-flexTMPCR擴增儀(ABI美國)上進行,在Omega F SNP熒光檢測儀(LGC,英國)讀取數(shù)據(jù),結(jié)果利用軟件Kluster Caller獲得待測樣品的SNP位點的數(shù)據(jù)。

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

    根據(jù)SNP分型結(jié)果,利用PowerMarker軟件計算各材料間的遺傳距離,在此基礎(chǔ)上采用加權(quán)平均法(unweighted pair-group method with arithmetic means,UPGMA)構(gòu)建聚類圖,使用MEGA 5.0軟件繪制聚類圖;并計算每個SNP位點的等位變異數(shù)、等位基因頻率、基因多樣性、雜合率、多態(tài)性信息含量(polymorphism information content,PIC)。利用Structure 2.3.4軟件分析群體遺傳結(jié)構(gòu),估算最佳群體組群數(shù)K,并計算Q值。利用GenALex 6.5軟件進行主成分分析。

    2 結(jié)果與分析

    2.1 番茄種質(zhì)資源遺傳多樣性分析

    利用60對SNP標記對504份番茄材料的遺傳多樣性進行分析(表1),結(jié)果表明,60對SNP引物共檢測到181個等位基因,平均每個位點有3個等位基因,NDP7582、NDP7560、NDP7783檢測到的等位基因數(shù)最多,為4個;60個SNP標記的主要等位基因頻率(major allele frequency, MAF)平均值為0.659,變幅為0.375~0.905,其中NDP7560最高,NDP7781最低;期望雜合率(Heterozygosity)平均值為0.069;基因多樣性在0.178~0.658之間,平均值為0.450,其中DNP7781的基因多樣性最高,NDP7560的基因多樣性最低,基因多樣性在0.600以上的有NDP7723(0.625)、NDP7529(0.623)、NDP7625(0.652)、DNP7596(0.643)、DNP7661(0.612)、DNP7781(0.658),表明番茄種質(zhì)資源有一定的多樣性。

    表1 60個SNP位點的多態(tài)性分析Table 1 Genetic diversity analysis of 60 SNP loci

    表1(續(xù))

    由表1可知,多態(tài)性信息值(PIC)變幅為0.171~0.583,平均值為0.381,PIC在0.5以上的有NDP7723(0.554)、NDP7589(0.501)、NDP7529(0.551)、NDP7710(0.508)、DNP7625(0.577)、NDP7596(0.569)、DNP7661(0.537)、NDP7368(0.512)、DNP7781(0.583),共9個SNP標記。Botstein等[18]指出當PIC值>0.5時,表明其具有高度多態(tài)性;當0.25

    圖1 SNP標記多態(tài)性信息值(PIC)柱形圖Fig.1 Histogram of SNP marker polymorphism information value (PIC)

    2.2 番茄種質(zhì)資源的聚類分析與主成分分析

    基于遺傳距離構(gòu)建504份番茄種質(zhì)資源的UPGMA聚類圖,結(jié)果如圖2所示。504份番茄種質(zhì)資源的遺傳距離范圍在0.02~2.64之間,平均遺傳距離為0.62,在遺傳距離為0.36時被分為7類。

    圖2 基于SNP標記的504份番茄種質(zhì)資源聚類圖Fig.2 Cluster map of 504 tomato germplasm resources based on SNP marker

    第Ⅰ類包含S460和S458共2份材料,均為來自甘肅的紅果櫻桃番茄,植株長勢中等,早熟;S458的生長習性為有限生長型,果形為橢圓,果實較硬;S460的生長習性為無限生長型,果形為圓形,果實較軟。

    第Ⅱ類包含120份材料,其中來自寧夏27份,甘肅73份,四川和云南各3份,遼寧10份,陜西4份;這一類群的材料中有15份材料為櫻桃番茄,其余均為大果番茄;大多數(shù)材料為無限生長型、生長勢較強的中晚熟番茄,果形主要為圓形,成熟果色有紅色、粉色和黃色三類。

    第Ⅲ類包含8份材料,其中S055來自甘肅,其余均來自寧夏;S055、S056、S057材料為櫻桃番茄,S058、 S059、 S060、 S061和S062材料為口感型的大果番茄,這8份材料均為植株長勢強的無限生長類型早熟番茄,果形有圓形和扁圓形,果實較軟,成熟果色有粉色和紫色。

    第Ⅳ類別包含72份材料,均為櫻桃番茄,其中來自甘肅24份,寧夏43份,陜西5份;這一類群材料主要為無限生長型、生長勢較強的早熟和中早熟番茄,果形主要為長圓形,果實較硬,成熟果色有紅、粉、黃色三類。

    第Ⅴ類和第Ⅵ類均包含2份材料,其中第Ⅴ類番茄S064和S179來自甘肅,均為無限生長型、生長勢強的早熟番茄品種,果形為圓形,果實較軟;S064成熟果色為粉色,S179成熟果色為紅色。第Ⅵ類2份番茄材料來自寧夏,均為長勢強的無限生長型早熟紅果番茄,果形為圓形,果實較硬。

    第Ⅶ類包含298份材料,其中來自寧夏118份、甘肅144份、四川8份、遼寧10份、陜西14份和山東5份;這類材料大多數(shù)為無限生長型、生長勢較強的早熟大果番茄,果形主要是圓形和扁圓形,成熟果色有紅色和粉色。

    第Ⅱ類的甘肅粉果材料S435與第Ⅶ類的遼寧紅果材料S372間的遺傳距離最大,為2.64,說明這兩個材料之間的親緣關(guān)系較遠,遺傳背景復雜且遺傳多樣性較高;第Ⅶ類S212與S214材料間遺傳距離最小,為0.02,二者均來源于甘肅地區(qū)且果實不容易轉(zhuǎn)色,說明這兩個材料之間親緣關(guān)系較近。

    由主成分分析二維圖可知(圖3),第1主成分與第2主成分的貢獻率分別為11.4%和10.01%。二維圖中的位置距離能夠直觀地代表番茄種質(zhì)資源的親緣關(guān)系的遠近,group1與group2主要分布在第1、第2、第4象限,且兩個群體之間有重疊現(xiàn)象,group3集中分布在第2、第3象限,且與其他群體之間基本沒有重疊,證明group3包含的材料與其他兩組的親緣關(guān)系較遠。

    圖3 基于SNP標記的504個番茄種質(zhì)資源主成分分析二維圖Fig.3 Two-dimensional map of principal component analysis of 504 tomato germplasm based on SNP markers

    2.3 群體遺傳結(jié)構(gòu)分析

    基于SNP分子標記數(shù)據(jù),運用Structure 2.3.4軟件對504份番茄種質(zhì)資源進行群體遺傳結(jié)構(gòu)分析,將分析得到的dimensional K繪制折線圖(圖4),結(jié)果表明在K=3時,ΔK達到最大值。因此,將504份番茄種質(zhì)資源劃分為3個類群(圖5),第一個群體(紅色)包含79份材料,第二個群體(綠色)包含302份材料,第三個群體(藍色)包含123份材料。參考Evanno等[19]研究發(fā)現(xiàn),當Q>0.6時,則表示該材料遺傳背景簡單,反之則表示遺傳背景復雜。從各類群的Q值分布表發(fā)現(xiàn)(表2),95.44%的番茄Q值大于等于0.6,有70.24%的番茄Q值≥0.9,表明大多數(shù)參試材料的遺傳背景簡單,大數(shù)為純合或者高代材料;4.56%的參試材料的Q值小于0.6,表明這23份番茄種質(zhì)資源可能為雜合材料。

    圖4 K值與ΔK的變化折線圖Fig.4 The change of K value and ΔK line diagram

    表2 各類群Q值分布Table 2 Q value distribution of various groups

    3 討論

    優(yōu)異的種質(zhì)資源是新品種的選育基礎(chǔ),種質(zhì)資源遺傳多樣性評價與鑒定是挖掘優(yōu)良育種材料、提高育種效率的有效工具[20]。分子標記技術(shù)是鑒定和評價番茄種質(zhì)資源遺傳多樣性的有效手段。鄧學斌等[21]利用46對SNP標記對收集的不同類型的加工番茄種質(zhì)資源進行遺傳多樣性評價,發(fā)現(xiàn)3 026份資源的期望雜合度均值為0.435,平均PIC為0.338。劉希艷等[22]利用44對SNP標記對615份番茄材料進行研究,表明番茄材料的基因多樣性均值為0.470,平均PIC為0.392 1;王濤等[20]利用篩選的103對SNP標記對210份番茄資源進行基因分型,結(jié)果表明平均期望雜合度為0.029,基因多樣性均值為0.346,平均PIC為0.308。在本研究中,利用篩選的60對SNP分子標記運用KASP技術(shù)對收集的504份番茄種質(zhì)資源進行了遺傳多樣性分析,共檢測到181個等位基因,期望雜合度均值為0.069,基因多樣性均值為0.450,平均PIC為0.381。比較發(fā)現(xiàn),本研究收集的番茄種質(zhì)資源的期望雜合度(0.069)與王濤等[20]的研究結(jié)果基本一致(0.029),遠低于鄧學斌等[21]的研究結(jié)果(0.435),原因可能是本研究參試的番茄材料多數(shù)為高代或純合材料,雜合率較低,這與群體遺傳結(jié)構(gòu)分析結(jié)果基本相同,鄧學斌等[21]利用的番茄材料全部為單株多代自交材料,而非高代或純合材料,因此期望雜合度較高。基因多樣性的大小(0.450)與劉希艷等[22]的研究結(jié)果(0.470)相當,高于王濤等[20]的研究結(jié)果(0.389),可能的原因是本研究的試材數(shù)量較多,種類多樣,而王濤等[20]的參試材料主要為鮮食番茄,類型較少。本研究篩選的60對SNP標記的PIC平均值較高(0.381),表明本研究篩選的SNP分子標記表現(xiàn)為中度偏高的多態(tài)性,這些SNP標記可以用于番茄種質(zhì)資源的遺傳多樣性分析和后期其他遺傳關(guān)系的研究。

    根據(jù)聚類分析結(jié)果將504份番茄材料分為7類,第Ⅰ類材料均為植株長勢中等、早熟的紅果櫻桃番茄,可用于早熟紅果櫻桃番茄的培育;第Ⅱ類材料大多數(shù)為植株長勢較強、無限生長的中晚熟番茄,果色豐富,遺傳潛力較大,可用于培育各類型中晚熟大果或櫻桃番茄品種,滿足周年生產(chǎn)的需求;第Ⅲ類材料均為植株長勢強、無限生長型的早熟番茄,包括5份口感型大果番茄材料和2份紫色櫻桃番茄材料,因此,這類材料可用于培育口感型粉色大果或紫色櫻桃番茄等特色品種;第Ⅳ類材料均為果實較硬、果色豐富、早熟的長圓形櫻桃番茄,這類材料可用于培育耐貯運的長圓形櫻桃番茄品種;第Ⅴ、第Ⅵ和第Ⅶ類材料多數(shù)為生長勢強、無限生長型的圓形或扁圓形的大果早熟番茄,可用于培育早熟大果鮮食番茄品種。遺傳距離可判斷作物種質(zhì)資源的親緣關(guān)系,遺傳距離越大親緣關(guān)系越遠,反之亦然。本研究的504份參試材料有490份材料均集中在第Ⅱ、第Ⅳ、第Ⅶ類;材料間的遺傳距離的變化范圍0.02~2.64,平均值為0.62,說明本研究參試材料的具有一定遺傳差異和遺傳多樣性,但豐富性仍不高,可能的原因是供試材料均為地方選育的栽培種。Sacco等[23]和Sim等[24]的研究結(jié)果表明,地方栽培品種的遺傳變異水平較低,推測是由于育種者通常從最好的果實中收集種子,而不是收集更優(yōu)異的基因型[25]。因此,種質(zhì)資源的收集至關(guān)重要。本課題組已經(jīng)開始了對國外番茄種質(zhì)的收集,包括野生種質(zhì)與栽培品種,以此豐富本課題組的番茄種質(zhì)資源。

    主成分分析與群體結(jié)構(gòu)分析均劃分為三類,與聚類分析劃分類群數(shù)不同,類群Ⅳ與group1和群體Ⅰ、類群Ⅶ與group2和群體Ⅱ、類群Ⅱ與group3和群體Ⅲ包含的番茄材料基本一致。從三種不同的聚類方法能夠發(fā)現(xiàn)參試的504份種質(zhì)資源的聚類結(jié)果與地理來源沒有明顯的關(guān)系,原因在于:一方面可能是種質(zhì)資源的遺傳背景不夠明確,無法細分,番茄種質(zhì)通過自交,導致種內(nèi)出現(xiàn)基因交流,從群體結(jié)構(gòu)也可發(fā)現(xiàn)種群內(nèi)存在基因交流現(xiàn)象;另一方面可能是引物數(shù)量較少,檢測位點未能覆蓋所有變異,這與周艷朝等[26]和胡永超等[27]的推斷一致。

    4 結(jié)論

    本研究利用60對SNP引物,對收集的504份番茄種質(zhì)資源進行遺傳多樣性分析,結(jié)果表明,每個SNP標記平均有3個等位基因,多態(tài)性信息值(PIC)平均值為0.381;同時發(fā)現(xiàn)有NDP7723(0.554)、NDP7589(0.501)、NDP7529(0.551)、NDP7710(0.508)、DNP7625(0.577)、NDP7596(0.569)、DNP7661(0.537)、NDP7368(0.512)、DNP7781(0.583)共9個具有高度多態(tài)性較高的SNP標記?;赟NP標記對參試材料進行UPGMA聚類分析、主成分分析和群體遺傳結(jié)構(gòu)分析,了解到番茄種質(zhì)資源之間存在一定的親緣關(guān)系,反映了不同番茄材料之間存在基因交流。

    猜你喜歡
    種質(zhì)寧夏多態(tài)性
    華南地區(qū)最大農(nóng)作物種質(zhì)資源保護庫建成
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進展
    寧夏
    寧夏畫報(2019年5期)2019-09-19 14:42:26
    亞麻抗白粉病種質(zhì)資源的鑒定與篩選
    寧夏煤電博弈
    能源(2017年5期)2017-07-06 09:25:50
    寧夏
    寧夏畫報(2016年5期)2016-06-28 05:48:08
    貴州玉米種質(zhì)資源遺傳多樣性及核心種質(zhì)庫構(gòu)建
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    紅錐種質(zhì)早期生長表現(xiàn)
    GlobalFiler~? PCR擴增試劑盒驗證及其STR遺傳多態(tài)性
    亚洲精品乱码久久久久久按摩| 日韩欧美三级三区| 欧美zozozo另类| 国产精品一区www在线观看| 久久久久久久久久久丰满| 日日摸夜夜添夜夜添av毛片| 日本黄色视频三级网站网址| 国产欧美另类精品又又久久亚洲欧美| 国产在视频线精品| 国产精品久久电影中文字幕| 熟妇人妻久久中文字幕3abv| 久久国产乱子免费精品| 床上黄色一级片| 日韩欧美精品免费久久| 久久精品国产鲁丝片午夜精品| 特级一级黄色大片| 欧美日本亚洲视频在线播放| 久久久精品欧美日韩精品| 国产淫片久久久久久久久| av黄色大香蕉| 国产熟女欧美一区二区| 欧美bdsm另类| 色播亚洲综合网| 成人av在线播放网站| 超碰97精品在线观看| 欧美3d第一页| 最近最新中文字幕大全电影3| 国产淫语在线视频| 欧美日本视频| 国产精品一区二区性色av| 亚洲国产精品专区欧美| 丰满少妇做爰视频| 尤物成人国产欧美一区二区三区| 精品国产一区二区三区久久久樱花 | videos熟女内射| 国产精品国产高清国产av| 亚洲精品,欧美精品| 免费av毛片视频| 女人十人毛片免费观看3o分钟| 成人鲁丝片一二三区免费| 亚洲人成网站在线观看播放| 波多野结衣高清无吗| АⅤ资源中文在线天堂| 人妻少妇偷人精品九色| 可以在线观看毛片的网站| av.在线天堂| 99久久中文字幕三级久久日本| 免费搜索国产男女视频| 直男gayav资源| 亚洲欧美日韩东京热| 国产精品国产三级专区第一集| 国产精品人妻久久久影院| videossex国产| 欧美另类亚洲清纯唯美| 人人妻人人澡人人爽人人夜夜 | 日韩精品有码人妻一区| 亚洲国产精品专区欧美| 国产精品女同一区二区软件| 一级黄片播放器| 麻豆一二三区av精品| 日韩欧美 国产精品| 97超视频在线观看视频| 1024手机看黄色片| 午夜视频国产福利| 神马国产精品三级电影在线观看| 国产亚洲最大av| 午夜爱爱视频在线播放| 91精品一卡2卡3卡4卡| 汤姆久久久久久久影院中文字幕 | 噜噜噜噜噜久久久久久91| 麻豆成人午夜福利视频| 国产伦精品一区二区三区四那| 午夜福利网站1000一区二区三区| 国产亚洲最大av| 一级黄片播放器| 麻豆一二三区av精品| 中文字幕av成人在线电影| 天堂√8在线中文| 欧美色视频一区免费| 久久久精品大字幕| 天天躁夜夜躁狠狠久久av| 床上黄色一级片| 亚洲成人久久爱视频| 免费人成在线观看视频色| 国产精品一区二区性色av| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩高清专用| 啦啦啦观看免费观看视频高清| 国产成人精品久久久久久| 国产欧美另类精品又又久久亚洲欧美| 人妻少妇偷人精品九色| 久久久久久久久久成人| 日韩精品有码人妻一区| videossex国产| videos熟女内射| 亚洲第一区二区三区不卡| 亚洲精品乱久久久久久| 51国产日韩欧美| 小蜜桃在线观看免费完整版高清| 日本色播在线视频| 岛国在线免费视频观看| 久热久热在线精品观看| 国产毛片a区久久久久| 成人国产麻豆网| 毛片一级片免费看久久久久| 亚洲av成人精品一区久久| 麻豆成人av视频| av免费观看日本| 丝袜美腿在线中文| 老师上课跳d突然被开到最大视频| 精品国产露脸久久av麻豆 | 最近最新中文字幕免费大全7| 一区二区三区乱码不卡18| 国产大屁股一区二区在线视频| 国产老妇女一区| 日韩大片免费观看网站 | 国产在线男女| 国产精品国产三级专区第一集| 一级二级三级毛片免费看| 成年免费大片在线观看| 国产毛片a区久久久久| 永久免费av网站大全| 久久久久久久午夜电影| 亚洲国产成人一精品久久久| 国产av码专区亚洲av| 日韩三级伦理在线观看| 秋霞在线观看毛片| 一个人看的www免费观看视频| av免费在线看不卡| 亚洲成人久久爱视频| 亚洲国产日韩欧美精品在线观看| 26uuu在线亚洲综合色| 性插视频无遮挡在线免费观看| 尤物成人国产欧美一区二区三区| 久久久亚洲精品成人影院| 综合色av麻豆| 高清毛片免费看| 色网站视频免费| 国产成人91sexporn| 亚洲国产最新在线播放| 午夜视频国产福利| 97人妻精品一区二区三区麻豆| 91在线精品国自产拍蜜月| 中文字幕制服av| 在线观看66精品国产| 国产成人免费观看mmmm| 精品欧美国产一区二区三| 天堂网av新在线| 久久久亚洲精品成人影院| 综合色av麻豆| 色尼玛亚洲综合影院| 全区人妻精品视频| a级毛片免费高清观看在线播放| 秋霞在线观看毛片| 色综合色国产| 又粗又硬又长又爽又黄的视频| 午夜激情福利司机影院| 久久久久久久久久成人| 色综合亚洲欧美另类图片| 午夜福利视频1000在线观看| 大又大粗又爽又黄少妇毛片口| 午夜福利高清视频| 日本猛色少妇xxxxx猛交久久| 中文字幕精品亚洲无线码一区| 中文资源天堂在线| 国产精品一及| 麻豆精品久久久久久蜜桃| 亚洲无线观看免费| 日本与韩国留学比较| 狠狠狠狠99中文字幕| 高清毛片免费看| 国产午夜福利久久久久久| 99在线人妻在线中文字幕| 国产麻豆成人av免费视频| 欧美日本视频| 日本wwww免费看| 午夜精品在线福利| 色5月婷婷丁香| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 久久久色成人| 最近最新中文字幕大全电影3| 久久久久久九九精品二区国产| 一级二级三级毛片免费看| 亚洲美女搞黄在线观看| 寂寞人妻少妇视频99o| 97在线视频观看| 婷婷色av中文字幕| 最近中文字幕2019免费版| АⅤ资源中文在线天堂| 老司机福利观看| 国产av码专区亚洲av| 97在线视频观看| 老司机福利观看| 亚洲电影在线观看av| 舔av片在线| 别揉我奶头 嗯啊视频| 成人毛片a级毛片在线播放| 免费看美女性在线毛片视频| 国产成人a区在线观看| 亚洲精品久久久久久婷婷小说 | 日日摸夜夜添夜夜添av毛片| 美女黄网站色视频| 色视频www国产| 国产黄片美女视频| 噜噜噜噜噜久久久久久91| 内地一区二区视频在线| 日韩精品青青久久久久久| 91久久精品国产一区二区成人| 国产高清三级在线| 免费观看a级毛片全部| 亚洲四区av| 麻豆成人av视频| 国产成年人精品一区二区| 国产69精品久久久久777片| 欧美最新免费一区二区三区| 午夜亚洲福利在线播放| 国产欧美日韩精品一区二区| 色综合亚洲欧美另类图片| 少妇人妻精品综合一区二区| 亚洲中文字幕日韩| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| 久久99热这里只有精品18| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 国产又黄又爽又无遮挡在线| 日韩中字成人| 亚洲欧美中文字幕日韩二区| 人人妻人人澡人人爽人人夜夜 | 精品免费久久久久久久清纯| 国产在线男女| 亚洲国产欧美在线一区| 高清毛片免费看| 神马国产精品三级电影在线观看| 久久久久网色| 美女脱内裤让男人舔精品视频| 亚洲av电影在线观看一区二区三区 | av在线播放精品| 国产日韩欧美在线精品| 日本wwww免费看| 中文字幕精品亚洲无线码一区| 黄色日韩在线| 亚洲av一区综合| 亚洲在线自拍视频| 久久久亚洲精品成人影院| 亚洲国产精品成人综合色| 亚洲丝袜综合中文字幕| 日韩欧美三级三区| 久久精品国产鲁丝片午夜精品| 亚洲av中文av极速乱| a级毛片免费高清观看在线播放| 国产美女午夜福利| 亚洲国产精品合色在线| 亚洲国产欧洲综合997久久,| 亚洲国产色片| 波野结衣二区三区在线| 亚洲av熟女| a级一级毛片免费在线观看| 亚洲欧美日韩无卡精品| 久久综合国产亚洲精品| 国产亚洲av片在线观看秒播厂 | 又爽又黄a免费视频| 国产精品福利在线免费观看| 精品欧美国产一区二区三| 91狼人影院| 91av网一区二区| 久久久久网色| 看十八女毛片水多多多| 久久久久久大精品| 中文字幕人妻熟人妻熟丝袜美| 久久久亚洲精品成人影院| 99国产精品一区二区蜜桃av| 联通29元200g的流量卡| 中文欧美无线码| 欧美xxxx黑人xx丫x性爽| 一级毛片aaaaaa免费看小| 久久久久久伊人网av| 国产美女午夜福利| 成人亚洲欧美一区二区av| 国产男人的电影天堂91| 亚洲人成网站在线播| 国内精品美女久久久久久| 乱人视频在线观看| 十八禁国产超污无遮挡网站| 免费播放大片免费观看视频在线观看 | 一边摸一边抽搐一进一小说| 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看 | 床上黄色一级片| 日本黄色片子视频| 精品国产露脸久久av麻豆 | 中文精品一卡2卡3卡4更新| 一级黄片播放器| av福利片在线观看| 一级毛片电影观看 | 91久久精品电影网| 99九九线精品视频在线观看视频| 别揉我奶头 嗯啊视频| 少妇人妻一区二区三区视频| av女优亚洲男人天堂| 日韩国内少妇激情av| 久久精品综合一区二区三区| 99久久人妻综合| 在线观看av片永久免费下载| 好男人视频免费观看在线| 国产淫语在线视频| 成人国产麻豆网| 51国产日韩欧美| 夫妻性生交免费视频一级片| 我要搜黄色片| a级一级毛片免费在线观看| 久久精品国产自在天天线| 亚洲国产色片| 欧美精品一区二区大全| 亚洲婷婷狠狠爱综合网| 91狼人影院| 欧美一区二区国产精品久久精品| 久久精品国产亚洲av天美| 亚洲精品自拍成人| 2021少妇久久久久久久久久久| 三级经典国产精品| 国产黄a三级三级三级人| 亚洲av电影不卡..在线观看| 汤姆久久久久久久影院中文字幕 | 日日摸夜夜添夜夜爱| 国产精品不卡视频一区二区| 丰满人妻一区二区三区视频av| 美女高潮的动态| av在线蜜桃| 亚洲综合精品二区| 日本熟妇午夜| a级毛片免费高清观看在线播放| 精品国内亚洲2022精品成人| 免费av毛片视频| 色视频www国产| 亚洲成人av在线免费| 久久韩国三级中文字幕| 中文资源天堂在线| 久久人人爽人人爽人人片va| 日本一本二区三区精品| 亚洲怡红院男人天堂| 在线免费观看不下载黄p国产| 女人久久www免费人成看片 | 久久久久久久亚洲中文字幕| 国产不卡一卡二| 午夜激情福利司机影院| 两个人的视频大全免费| 岛国在线免费视频观看| 免费一级毛片在线播放高清视频| 日韩三级伦理在线观看| 嘟嘟电影网在线观看| 国产精品综合久久久久久久免费| 国产黄色视频一区二区在线观看 | 91久久精品国产一区二区成人| 国产探花极品一区二区| 亚洲18禁久久av| 精品99又大又爽又粗少妇毛片| 少妇人妻一区二区三区视频| 2021少妇久久久久久久久久久| 国产av码专区亚洲av| 亚洲国产最新在线播放| 亚洲色图av天堂| 少妇猛男粗大的猛烈进出视频 | 简卡轻食公司| 在线天堂最新版资源| av国产久精品久网站免费入址| 日韩在线高清观看一区二区三区| 亚洲精品色激情综合| 老司机福利观看| 最近中文字幕2019免费版| 午夜老司机福利剧场| 黄色配什么色好看| 狠狠狠狠99中文字幕| 欧美变态另类bdsm刘玥| 国产黄a三级三级三级人| 观看免费一级毛片| 美女国产视频在线观看| 级片在线观看| av免费观看日本| 国产老妇伦熟女老妇高清| 岛国毛片在线播放| 久久99热这里只频精品6学生 | 日韩精品有码人妻一区| 国产中年淑女户外野战色| 国产精品国产高清国产av| 久久精品综合一区二区三区| 成人性生交大片免费视频hd| 国产一级毛片七仙女欲春2| 午夜福利在线观看吧| 免费av不卡在线播放| 最近中文字幕高清免费大全6| 天堂av国产一区二区熟女人妻| 久久精品综合一区二区三区| 成人一区二区视频在线观看| 国产精品av视频在线免费观看| 女人被狂操c到高潮| 精品无人区乱码1区二区| 久久久亚洲精品成人影院| 免费黄色在线免费观看| 男人舔奶头视频| 国产老妇伦熟女老妇高清| 久久久久久九九精品二区国产| 日日摸夜夜添夜夜爱| 天堂网av新在线| 高清av免费在线| 干丝袜人妻中文字幕| 精华霜和精华液先用哪个| 国产亚洲一区二区精品| 精品久久久久久久久亚洲| 人妻系列 视频| 村上凉子中文字幕在线| 韩国av在线不卡| 99视频精品全部免费 在线| 美女内射精品一级片tv| 免费看光身美女| 久久精品久久精品一区二区三区| 国产探花极品一区二区| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 成人午夜高清在线视频| 99热精品在线国产| 又粗又硬又长又爽又黄的视频| 国产精品蜜桃在线观看| 亚洲精品乱码久久久久久按摩| 国产视频内射| 国产一区二区在线观看日韩| 亚洲国产成人一精品久久久| 看非洲黑人一级黄片| 欧美成人午夜免费资源| 亚洲精品成人久久久久久| 欧美极品一区二区三区四区| 久久久久精品久久久久真实原创| av线在线观看网站| 亚洲国产欧洲综合997久久,| 亚洲在线观看片| 久久精品影院6| 国产精品永久免费网站| 亚洲av成人精品一二三区| 一级av片app| 欧美bdsm另类| 日韩一区二区三区影片| 久久久久精品久久久久真实原创| 男女那种视频在线观看| 国产久久久一区二区三区| 人人妻人人澡欧美一区二区| 国产美女午夜福利| 亚洲精品亚洲一区二区| 如何舔出高潮| 色噜噜av男人的天堂激情| 级片在线观看| 国产精品人妻久久久影院| 麻豆一二三区av精品| 精品久久久久久久久av| 秋霞伦理黄片| 国产一区二区亚洲精品在线观看| 久久亚洲国产成人精品v| 舔av片在线| 91狼人影院| 丰满人妻一区二区三区视频av| 亚洲丝袜综合中文字幕| 18禁在线播放成人免费| 国产探花极品一区二区| 欧美xxxx黑人xx丫x性爽| 亚洲成色77777| 一本一本综合久久| 欧美精品一区二区大全| 一级av片app| 日韩欧美精品v在线| 直男gayav资源| 亚洲av成人av| 国产探花在线观看一区二区| 少妇裸体淫交视频免费看高清| 麻豆成人av视频| 亚洲精品影视一区二区三区av| 日韩国内少妇激情av| 一边摸一边抽搐一进一小说| 伦理电影大哥的女人| 久久久久九九精品影院| 91午夜精品亚洲一区二区三区| 国产麻豆成人av免费视频| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 97超视频在线观看视频| 啦啦啦韩国在线观看视频| 中国国产av一级| 久久久久免费精品人妻一区二区| 舔av片在线| 人人妻人人看人人澡| 国产精品一及| 国产亚洲av片在线观看秒播厂 | 亚洲欧美日韩东京热| 男女边吃奶边做爰视频| 精品久久久噜噜| 久久亚洲精品不卡| 乱码一卡2卡4卡精品| 九九热线精品视视频播放| 欧美日本视频| 久久午夜福利片| 精品国产露脸久久av麻豆 | 性色avwww在线观看| 全区人妻精品视频| 精品无人区乱码1区二区| 亚洲真实伦在线观看| 久久久国产成人免费| 国产 一区精品| 国产精品久久久久久久久免| 五月玫瑰六月丁香| 午夜视频国产福利| 乱人视频在线观看| 一级黄色大片毛片| 2021少妇久久久久久久久久久| 夜夜看夜夜爽夜夜摸| 99在线视频只有这里精品首页| 人妻夜夜爽99麻豆av| 亚洲精品456在线播放app| 26uuu在线亚洲综合色| 亚洲精品乱码久久久v下载方式| 亚洲国产欧美在线一区| 少妇的逼水好多| 最近视频中文字幕2019在线8| 白带黄色成豆腐渣| 日韩欧美精品免费久久| 人妻夜夜爽99麻豆av| 久久草成人影院| 少妇丰满av| 91精品国产九色| 亚洲自偷自拍三级| 免费人成在线观看视频色| 最近中文字幕高清免费大全6| 国产毛片a区久久久久| 国产精品久久视频播放| 亚洲精品自拍成人| 免费看av在线观看网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色综合亚洲欧美另类图片| 又爽又黄a免费视频| 激情 狠狠 欧美| 国产精品乱码一区二三区的特点| 天天躁日日操中文字幕| 22中文网久久字幕| 少妇人妻精品综合一区二区| 精品人妻偷拍中文字幕| 精品久久久久久成人av| 亚洲精品色激情综合| 91狼人影院| av在线播放精品| 色综合色国产| 在现免费观看毛片| 美女大奶头视频| 亚洲美女搞黄在线观看| 丰满乱子伦码专区| 高清视频免费观看一区二区 | 深爱激情五月婷婷| 日韩国内少妇激情av| 亚洲最大成人手机在线| 国产高潮美女av| 麻豆成人午夜福利视频| av在线观看视频网站免费| 99热全是精品| 美女国产视频在线观看| 免费人成在线观看视频色| 国产精品三级大全| 18+在线观看网站| 欧美一区二区国产精品久久精品| av天堂中文字幕网| 午夜精品在线福利| 亚洲自拍偷在线| 亚洲国产精品久久男人天堂| 草草在线视频免费看| 成人美女网站在线观看视频| 综合色av麻豆| 在线a可以看的网站| 日韩亚洲欧美综合| av播播在线观看一区| 久久久色成人| 寂寞人妻少妇视频99o| 国产美女午夜福利| 日韩制服骚丝袜av| 亚洲丝袜综合中文字幕| 国产免费视频播放在线视频 | 丰满人妻一区二区三区视频av| 啦啦啦观看免费观看视频高清| 婷婷色综合大香蕉| 久久人人爽人人片av| 在线观看66精品国产| 偷拍熟女少妇极品色| 国产欧美日韩精品一区二区| 国产黄片视频在线免费观看| 国产高清不卡午夜福利| 免费观看的影片在线观看| 午夜日本视频在线| 如何舔出高潮| 国产精品三级大全| 日韩国内少妇激情av| av女优亚洲男人天堂| 亚洲精品,欧美精品| 久久久久久久亚洲中文字幕| 国产精品爽爽va在线观看网站| 一级黄色大片毛片| av女优亚洲男人天堂| 亚洲精品,欧美精品| 亚洲久久久久久中文字幕| 色综合站精品国产| 日日摸夜夜添夜夜添av毛片| 六月丁香七月| 久久久午夜欧美精品| 欧美成人免费av一区二区三区| 成年女人永久免费观看视频| 天天一区二区日本电影三级| 嫩草影院入口| 欧美日本亚洲视频在线播放| 九九爱精品视频在线观看| 最近视频中文字幕2019在线8|