• <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)性
    精品熟女少妇av免费看| 欧美激情在线99| 日韩av在线大香蕉| 美女cb高潮喷水在线观看| 国产精品一区www在线观看| 美女 人体艺术 gogo| 极品教师在线视频| 久久久久久久亚洲中文字幕| 联通29元200g的流量卡| 国产精品不卡视频一区二区| 赤兔流量卡办理| 看片在线看免费视频| 国产伦一二天堂av在线观看| 国产成人freesex在线| 美女脱内裤让男人舔精品视频 | 国产黄a三级三级三级人| 国产一区二区激情短视频| 免费大片18禁| 女人十人毛片免费观看3o分钟| 免费搜索国产男女视频| 狂野欧美白嫩少妇大欣赏| 在线免费观看不下载黄p国产| 欧美成人免费av一区二区三区| 天堂影院成人在线观看| 亚洲国产精品国产精品| 十八禁国产超污无遮挡网站| 免费人成在线观看视频色| 久久鲁丝午夜福利片| 国产精品久久久久久久电影| 我的老师免费观看完整版| 日韩精品青青久久久久久| 村上凉子中文字幕在线| 欧美bdsm另类| 少妇的逼好多水| 亚洲五月天丁香| 麻豆精品久久久久久蜜桃| 级片在线观看| 亚洲av电影不卡..在线观看| 成人综合一区亚洲| 亚洲人成网站在线观看播放| 欧美最新免费一区二区三区| 久久精品人妻少妇| 看非洲黑人一级黄片| 国产私拍福利视频在线观看| 丰满人妻一区二区三区视频av| 99久久精品热视频| av在线天堂中文字幕| 亚洲va在线va天堂va国产| 夜夜夜夜夜久久久久| 中文字幕制服av| 联通29元200g的流量卡| 老女人水多毛片| 久久亚洲国产成人精品v| 亚洲人成网站在线播放欧美日韩| av在线观看视频网站免费| 日韩欧美精品v在线| 欧美bdsm另类| 一边亲一边摸免费视频| 国产高清有码在线观看视频| 乱系列少妇在线播放| 亚洲天堂国产精品一区在线| 亚洲性久久影院| 婷婷亚洲欧美| 麻豆久久精品国产亚洲av| 人妻夜夜爽99麻豆av| www日本黄色视频网| 精品久久久久久久久av| 欧美成人一区二区免费高清观看| 99在线视频只有这里精品首页| 秋霞在线观看毛片| 高清毛片免费看| 色尼玛亚洲综合影院| 久久久国产成人精品二区| 尤物成人国产欧美一区二区三区| 看黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看| 精品不卡国产一区二区三区| 不卡一级毛片| 此物有八面人人有两片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 熟女电影av网| 91久久精品电影网| 亚洲av免费在线观看| 一夜夜www| 久99久视频精品免费| 亚洲图色成人| 亚洲成a人片在线一区二区| 久久久国产成人精品二区| 色尼玛亚洲综合影院| 久久久久国产网址| 国产三级中文精品| 91麻豆精品激情在线观看国产| 有码 亚洲区| 能在线免费看毛片的网站| 成人性生交大片免费视频hd| 亚洲第一电影网av| 久久久久久久久中文| 色哟哟哟哟哟哟| 一区二区三区高清视频在线| 日本五十路高清| 寂寞人妻少妇视频99o| 亚洲乱码一区二区免费版| 亚洲精品色激情综合| 国产av不卡久久| 97热精品久久久久久| 亚洲av电影不卡..在线观看| 精品少妇黑人巨大在线播放 | 少妇高潮的动态图| 最好的美女福利视频网| 少妇的逼好多水| 国产亚洲精品久久久com| 欧洲精品卡2卡3卡4卡5卡区| 免费av观看视频| 国产精品久久久久久久电影| 2022亚洲国产成人精品| 日韩欧美一区二区三区在线观看| 亚洲av一区综合| 一进一出抽搐动态| 一级毛片aaaaaa免费看小| 国产精品一区二区性色av| 欧美日韩在线观看h| 久久精品影院6| 成人午夜高清在线视频| 久久久久性生活片| 别揉我奶头 嗯啊视频| 日本成人三级电影网站| 国内揄拍国产精品人妻在线| 国产成人a∨麻豆精品| 内地一区二区视频在线| 黄色日韩在线| 一区福利在线观看| av卡一久久| 日本黄色片子视频| 尾随美女入室| 中文欧美无线码| 一级二级三级毛片免费看| 日韩亚洲欧美综合| 又黄又爽又刺激的免费视频.| 欧美最新免费一区二区三区| 亚洲自拍偷在线| 色哟哟·www| 麻豆国产av国片精品| 亚洲天堂国产精品一区在线| 中国美白少妇内射xxxbb| 2021天堂中文幕一二区在线观| 国产黄色小视频在线观看| 国产日本99.免费观看| 欧美变态另类bdsm刘玥| 亚洲成人久久性| 国产精品免费一区二区三区在线| 精品久久久久久久末码| 国产精品嫩草影院av在线观看| 国产精品电影一区二区三区| 欧美三级亚洲精品| 久久精品国产清高在天天线| 美女内射精品一级片tv| 亚洲不卡免费看| 不卡视频在线观看欧美| 欧美日本亚洲视频在线播放| 国产精品精品国产色婷婷| 非洲黑人性xxxx精品又粗又长| 麻豆一二三区av精品| 国内精品美女久久久久久| 国产伦精品一区二区三区视频9| 一级二级三级毛片免费看| 欧美zozozo另类| 一卡2卡三卡四卡精品乱码亚洲| 日本熟妇午夜| 国产真实伦视频高清在线观看| 激情 狠狠 欧美| 久久久久性生活片| 欧美变态另类bdsm刘玥| 两个人的视频大全免费| 91av网一区二区| 美女xxoo啪啪120秒动态图| 99热这里只有精品一区| 三级经典国产精品| 久久久久久伊人网av| 亚洲精品国产av成人精品| 老师上课跳d突然被开到最大视频| 伦理电影大哥的女人| 亚洲精品成人久久久久久| av视频在线观看入口| 麻豆乱淫一区二区| 国产精品不卡视频一区二区| 又粗又爽又猛毛片免费看| 亚洲无线在线观看| 午夜爱爱视频在线播放| 色视频www国产| 人妻制服诱惑在线中文字幕| 精品不卡国产一区二区三区| 男人舔奶头视频| 久久久久久国产a免费观看| 简卡轻食公司| 天堂av国产一区二区熟女人妻| 欧美三级亚洲精品| 亚洲一级一片aⅴ在线观看| 欧美成人一区二区免费高清观看| 桃色一区二区三区在线观看| 国产亚洲av片在线观看秒播厂 | 国产黄色小视频在线观看| 中文字幕免费在线视频6| 国产精品久久久久久久电影| 麻豆国产97在线/欧美| 成年免费大片在线观看| 欧美最新免费一区二区三区| 一边亲一边摸免费视频| 在线免费观看不下载黄p国产| 国产伦在线观看视频一区| 国产极品天堂在线| 少妇被粗大猛烈的视频| av卡一久久| 熟女人妻精品中文字幕| 国产日韩欧美在线精品| ponron亚洲| 色综合站精品国产| 99热这里只有是精品在线观看| 天天一区二区日本电影三级| 男女做爰动态图高潮gif福利片| 一本久久中文字幕| 啦啦啦韩国在线观看视频| 此物有八面人人有两片| 精品免费久久久久久久清纯| 女的被弄到高潮叫床怎么办| 国产白丝娇喘喷水9色精品| 美女国产视频在线观看| 给我免费播放毛片高清在线观看| 青春草视频在线免费观看| 日本免费一区二区三区高清不卡| 中文字幕熟女人妻在线| 美女高潮的动态| 国产白丝娇喘喷水9色精品| 99久久中文字幕三级久久日本| 久久久色成人| 欧美日本视频| 精华霜和精华液先用哪个| 午夜福利在线观看吧| 亚洲人成网站高清观看| 成人毛片a级毛片在线播放| 日韩强制内射视频| 欧美成人免费av一区二区三区| 夫妻性生交免费视频一级片| 在线免费观看的www视频| av女优亚洲男人天堂| 久久久久久久久中文| 国产精品福利在线免费观看| 好男人视频免费观看在线| a级一级毛片免费在线观看| eeuss影院久久| av女优亚洲男人天堂| 亚洲欧美成人精品一区二区| 少妇被粗大猛烈的视频| 亚洲av免费高清在线观看| 麻豆一二三区av精品| 黄片wwwwww| 99久国产av精品国产电影| 国产亚洲精品av在线| 久久亚洲精品不卡| 国产精品久久久久久精品电影小说 | 五月伊人婷婷丁香| av视频在线观看入口| 国产成人精品一,二区 | 热99re8久久精品国产| 国产精品99久久久久久久久| 日本色播在线视频| 黄色视频,在线免费观看| 欧美3d第一页| 少妇高潮的动态图| 国产精品免费一区二区三区在线| 国产精品一区www在线观看| 99久久精品一区二区三区| 欧美一区二区精品小视频在线| 色综合站精品国产| av女优亚洲男人天堂| 国产av麻豆久久久久久久| 国产精品久久久久久亚洲av鲁大| 国产高清有码在线观看视频| 晚上一个人看的免费电影| 日本熟妇午夜| 丰满的人妻完整版| 欧美成人一区二区免费高清观看| 午夜久久久久精精品| 91午夜精品亚洲一区二区三区| 亚洲av一区综合| 久久99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 国产精品福利在线免费观看| 国产av不卡久久| 最近2019中文字幕mv第一页| 亚洲精华国产精华液的使用体验 | 乱系列少妇在线播放| 老熟妇乱子伦视频在线观看| 国产亚洲91精品色在线| 美女cb高潮喷水在线观看| 国产探花在线观看一区二区| 91aial.com中文字幕在线观看| 成人无遮挡网站| 中文字幕av在线有码专区| 一区福利在线观看| 精品一区二区三区视频在线| 国产伦理片在线播放av一区 | 国产精品美女特级片免费视频播放器| 亚洲美女搞黄在线观看| av专区在线播放| 你懂的网址亚洲精品在线观看 | 日韩欧美精品v在线| 国产成人福利小说| 一级毛片我不卡| 午夜福利成人在线免费观看| 欧美一级a爱片免费观看看| 亚洲中文字幕一区二区三区有码在线看| 哪里可以看免费的av片| 在线免费十八禁| 麻豆成人av视频| 久久久久久久久大av| 非洲黑人性xxxx精品又粗又长| 夫妻性生交免费视频一级片| 26uuu在线亚洲综合色| 一个人免费在线观看电影| 桃色一区二区三区在线观看| 91精品一卡2卡3卡4卡| 日韩人妻高清精品专区| 精品一区二区免费观看| av视频在线观看入口| 欧美成人精品欧美一级黄| 中文资源天堂在线| 欧美又色又爽又黄视频| 久久精品影院6| 亚洲最大成人手机在线| 久久久久久久久久久丰满| .国产精品久久| 日韩制服骚丝袜av| 一区二区三区免费毛片| 亚洲第一电影网av| 日本成人三级电影网站| 久久精品国产亚洲av香蕉五月| 日本免费a在线| 能在线免费观看的黄片| 久久人妻av系列| av.在线天堂| 亚洲av成人av| 老司机影院成人| 久久这里有精品视频免费| av在线蜜桃| 国产精品国产三级国产av玫瑰| 18禁在线无遮挡免费观看视频| 特大巨黑吊av在线直播| 内射极品少妇av片p| av天堂在线播放| 国产一区二区激情短视频| 99久久精品国产国产毛片| 国产综合懂色| 成人特级黄色片久久久久久久| 在线免费观看的www视频| 老司机影院成人| 波多野结衣巨乳人妻| 22中文网久久字幕| 久久久成人免费电影| 校园人妻丝袜中文字幕| 欧美极品一区二区三区四区| 三级经典国产精品| 九九久久精品国产亚洲av麻豆| 日韩高清综合在线| 亚洲国产精品国产精品| 你懂的网址亚洲精品在线观看 | 日本与韩国留学比较| 搡老妇女老女人老熟妇| 精品99又大又爽又粗少妇毛片| 久久久久久久久久久免费av| 亚洲aⅴ乱码一区二区在线播放| 亚洲图色成人| 国产一级毛片七仙女欲春2| 亚洲欧洲日产国产| 又爽又黄a免费视频| 三级国产精品欧美在线观看| or卡值多少钱| 别揉我奶头 嗯啊视频| 国产精品综合久久久久久久免费| 嫩草影院新地址| 日本一本二区三区精品| 日韩av不卡免费在线播放| 亚洲欧美日韩卡通动漫| 黄色欧美视频在线观看| 亚洲精品日韩在线中文字幕 | 亚洲精华国产精华液的使用体验 | 99热6这里只有精品| 免费av观看视频| avwww免费| 久久国内精品自在自线图片| 赤兔流量卡办理| 五月伊人婷婷丁香| 一本久久精品| 午夜免费男女啪啪视频观看| 我的女老师完整版在线观看| 国产91av在线免费观看| 91午夜精品亚洲一区二区三区| 国产亚洲精品久久久com| 国产精品1区2区在线观看.| 中文字幕熟女人妻在线| 色尼玛亚洲综合影院| 婷婷色综合大香蕉| 日本黄色片子视频| 男的添女的下面高潮视频| 男插女下体视频免费在线播放| 婷婷亚洲欧美| 日本爱情动作片www.在线观看| 欧美高清成人免费视频www| 有码 亚洲区| 亚洲天堂国产精品一区在线| 国产 一区精品| 日韩 亚洲 欧美在线| 免费观看在线日韩| 精品免费久久久久久久清纯| 免费观看精品视频网站| 一区福利在线观看| 女人被狂操c到高潮| 日本免费一区二区三区高清不卡| 欧美成人一区二区免费高清观看| 国产成人freesex在线| 国产一区二区激情短视频| 久久精品国产亚洲av天美| 22中文网久久字幕| 好男人视频免费观看在线| 亚洲性久久影院| 国产男人的电影天堂91| 成熟少妇高潮喷水视频| 悠悠久久av| 毛片一级片免费看久久久久| 18+在线观看网站| 国产精品一区二区三区四区久久| 欧美高清性xxxxhd video| 久久99蜜桃精品久久| 国产人妻一区二区三区在| 国产亚洲精品久久久久久毛片| 亚洲aⅴ乱码一区二区在线播放| 国产精品不卡视频一区二区| 精品久久久久久久久av| 丝袜喷水一区| 久久久成人免费电影| 99久国产av精品| 最好的美女福利视频网| 国产精品美女特级片免费视频播放器| 亚洲电影在线观看av| 六月丁香七月| 欧美激情久久久久久爽电影| 亚洲av一区综合| 亚洲精华国产精华液的使用体验 | 看免费成人av毛片| 亚洲欧美清纯卡通| 亚洲av第一区精品v没综合| 国产av麻豆久久久久久久| 综合色av麻豆| 国产精品一区二区三区四区免费观看| 国产成人影院久久av| 91久久精品电影网| 一级av片app| 男女做爰动态图高潮gif福利片| av在线老鸭窝| 成人特级av手机在线观看| 一个人看视频在线观看www免费| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久电影| 国产精品美女特级片免费视频播放器| 哪个播放器可以免费观看大片| 亚洲国产精品国产精品| 国产一区二区激情短视频| 中国美白少妇内射xxxbb| 两个人视频免费观看高清| 99视频精品全部免费 在线| 午夜福利视频1000在线观看| 亚洲av中文av极速乱| 中文字幕制服av| 亚洲一级一片aⅴ在线观看| 看片在线看免费视频| 国产精品女同一区二区软件| 亚洲av二区三区四区| 日本欧美国产在线视频| 日产精品乱码卡一卡2卡三| 精品久久久久久成人av| 黑人高潮一二区| kizo精华| 99久久无色码亚洲精品果冻| 国产女主播在线喷水免费视频网站 | 国产蜜桃级精品一区二区三区| 国产成人精品一,二区 | 国产高清不卡午夜福利| 国产真实乱freesex| 2021天堂中文幕一二区在线观| 高清毛片免费看| av又黄又爽大尺度在线免费看 | 老熟妇乱子伦视频在线观看| 国内精品宾馆在线| 亚洲精品久久久久久婷婷小说 | 可以在线观看毛片的网站| 18禁在线播放成人免费| 欧美日韩综合久久久久久| 国产真实乱freesex| 美女大奶头视频| 亚洲第一电影网av| av免费在线看不卡| 乱码一卡2卡4卡精品| 美女黄网站色视频| 国产亚洲精品av在线| 一级二级三级毛片免费看| 亚洲无线观看免费| 一级毛片久久久久久久久女| 亚洲经典国产精华液单| 日韩 亚洲 欧美在线| 18禁黄网站禁片免费观看直播| 男女那种视频在线观看| 亚洲无线在线观看| 国语自产精品视频在线第100页| 男人和女人高潮做爰伦理| 黄片无遮挡物在线观看| 精品99又大又爽又粗少妇毛片| 欧美精品国产亚洲| 国内精品久久久久精免费| 少妇裸体淫交视频免费看高清| 小说图片视频综合网站| 少妇人妻精品综合一区二区 | 欧美xxxx性猛交bbbb| 久久国产乱子免费精品| 又爽又黄a免费视频| 久久久久久久午夜电影| 人妻制服诱惑在线中文字幕| 给我免费播放毛片高清在线观看| 亚洲成人精品中文字幕电影| 老司机福利观看| 久久精品影院6| 中文资源天堂在线| 午夜免费激情av| 又粗又爽又猛毛片免费看| 亚洲美女搞黄在线观看| 亚洲av成人av| 国产精品国产高清国产av| 99久久精品国产国产毛片| 亚洲欧美精品专区久久| 国产大屁股一区二区在线视频| 青春草亚洲视频在线观看| 夫妻性生交免费视频一级片| 免费在线观看成人毛片| 国产精品久久久久久精品电影| 人人妻人人看人人澡| av视频在线观看入口| 12—13女人毛片做爰片一| 国产v大片淫在线免费观看| 亚洲av成人av| 久久鲁丝午夜福利片| 午夜精品在线福利| 日韩中字成人| 国产精品电影一区二区三区| 久久6这里有精品| 麻豆久久精品国产亚洲av| 中文字幕制服av| 九色成人免费人妻av| 亚洲色图av天堂| 亚洲成人久久爱视频| 日韩制服骚丝袜av| 人人妻人人看人人澡| 能在线免费看毛片的网站| 午夜福利在线在线| 国产又黄又爽又无遮挡在线| 美女xxoo啪啪120秒动态图| 免费观看在线日韩| 国产单亲对白刺激| АⅤ资源中文在线天堂| 大香蕉久久网| 国产美女午夜福利| 日韩制服骚丝袜av| 啦啦啦观看免费观看视频高清| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 亚洲无线在线观看| 内射极品少妇av片p| 日本av手机在线免费观看| 亚洲国产精品合色在线| 一区二区三区免费毛片| 国产精品久久久久久久电影| 亚洲精品国产av成人精品| av.在线天堂| 亚洲最大成人手机在线| 亚洲一区二区三区色噜噜| 亚洲欧美日韩高清专用| 男的添女的下面高潮视频| 人妻少妇偷人精品九色| 神马国产精品三级电影在线观看| 欧美色欧美亚洲另类二区| av天堂在线播放| 一级av片app| 国产一区二区在线观看日韩| 亚洲av电影不卡..在线观看| 久久99热这里只有精品18| 老女人水多毛片| 精品午夜福利在线看| 欧美变态另类bdsm刘玥| 超碰av人人做人人爽久久| 熟妇人妻久久中文字幕3abv| 久久午夜福利片| 成年女人看的毛片在线观看| 国产精品久久久久久精品电影小说 | av女优亚洲男人天堂| av国产免费在线观看| 成人一区二区视频在线观看| 99国产精品一区二区蜜桃av| 午夜福利在线在线| 午夜久久久久精精品| 91在线精品国自产拍蜜月| 国产精品综合久久久久久久免费| 又爽又黄a免费视频| 欧美日韩乱码在线| 日韩人妻高清精品专区| av黄色大香蕉| 观看美女的网站| 久久精品人妻少妇| 尤物成人国产欧美一区二区三区|