999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于SNP標(biāo)記的玉米自交系S155、PHA458和A01的遺傳多樣性分析

2023-04-29 00:00:00孫佩張培風(fēng)張玉紅周聯(lián)東王蕊王文潔張瑞平李祥馬朝陽李合順王學(xué)軍
陜西農(nóng)業(yè)科學(xué) 2023年9期

摘要: 為了分析玉米品種新科891的遺傳結(jié)構(gòu)和遺傳多樣性,本研究利用Affymetrix Axiom 10K SNP芯片對(duì)玉米品種新科891的雙親自交系(母本S155、父本PHA458)及自交系A(chǔ)01進(jìn)行全基因組掃描,獲得多態(tài)性標(biāo)記4 926個(gè)。結(jié)果表明,群體聚類分析和群體遺傳結(jié)構(gòu)分析的雜種優(yōu)勢(shì)類群劃分結(jié)果基本一致。母本S155為SS類種質(zhì),父本PHA458含有NSS和塘四平頭類血緣,自交系A(chǔ)01為塘四平頭類血緣,新科891的雜種優(yōu)勢(shì)模式為SS×NSS,這為新科891雙親的后續(xù)改良與創(chuàng)新奠定基礎(chǔ),并為雜交組配的育種工作提供依據(jù)。

關(guān)鍵詞:玉米;新科891;SNP標(biāo)記;類群劃分;雜種優(yōu)勢(shì)模式

中圖分類號(hào):S513文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):0488-5368(2023)09-0001-07

收稿日期:2022-11-23修回日期:2023-01-25

基金項(xiàng)目:河南省重大科技專項(xiàng)“機(jī)收玉米新品種選育及配套生產(chǎn)技術(shù)研究”(212300111100);河南省農(nóng)業(yè)良種聯(lián)合攻關(guān)項(xiàng)目“高產(chǎn)抗病抗逆宜機(jī)收玉米種質(zhì)創(chuàng)新及新品種選育”(2022010201)。

第一作者簡(jiǎn)介:孫佩(1981-),女,河南新鄉(xiāng)人,副研究員,博士,主要從事玉米遺傳育種研究。

通信作者:周聯(lián)東。Genetic Diversity of Maize Inbred Lines S155, PHA458

and A01 Based on SNP Markers

SUN Pei, ZHANG Peifeng, ZHANG Yuhong, ZHOU Liandong, WANG Rui, WANG Wenjie,

ZHANG Ruiping, LI Xiang, MA Chaoyang, LI Heshun, WANG Xuejun

(Xinxiang Academy of Agricultural Sciences,Xinxiang, Henan" 453002," China)

Abstract:" In order to analyze the genetic structure and diversity of Xinke 891 maize, in this study, affymetrix Axiom 10K SNP chip was used to scan the whole genome of its two-parent inbred lines (female S155, male PHA458) and inbred line A01, and 4 926 polymorphic markers were obtained. The results showed that cluster analysis and genetic structure analysis for classification were basically consistent. The blood relationship between S155 and PHA458 was complicated, S155 belonged to SS germplasm, PHA458 contained NSS and TSPT, and A01 belonged to TSPT. The heterosis mode of Xinke 891 was SS×NSS. This study lays the foundation for the follow-up improvement and innovation of Xinke 891 parents, and provides a basis for breeding of crossbreeding.

Key words: Maize; Xinke 891; SNP markers; Genetic group classification; Heterosis model

優(yōu)良玉米品種的選育依賴于優(yōu)異豐富的玉米種質(zhì)資源。然而,當(dāng)前我國(guó)玉米種質(zhì)資源正面臨著日益狹窄的問題[1, 2]。通過引進(jìn)國(guó)內(nèi)地方種質(zhì)和國(guó)外種質(zhì)材料,對(duì)我國(guó)玉米自交系進(jìn)行不斷的改良和創(chuàng)新,可以形成很多骨干自交系,并在育種實(shí)踐中得以廣泛應(yīng)用[3]。與此同時(shí),對(duì)玉米種質(zhì)資源的遺傳多樣性進(jìn)行分析,來揭示種質(zhì)資源的遺傳特性,可以為玉米新品種的選育提供信息,并為品種的改良和選育提供參考依據(jù)[4, 5]。

合理的雜種優(yōu)勢(shì)類群劃分可以降低雜交組合組配的盲目性,可以提高玉米育種效率,加速玉米育種進(jìn)程,對(duì)玉米育種工作具有重要指導(dǎo)意義[6, 7]。隨著分子標(biāo)記技術(shù)的發(fā)展,利用分子標(biāo)記對(duì)玉米雜種優(yōu)勢(shì)類群劃分已成為新的手段[8]。作為目前最具有發(fā)展?jié)摿Φ姆肿訕?biāo)記,SNP(單核苷酸多態(tài)性)標(biāo)記因具有遺傳穩(wěn)定性更高、可靠性更強(qiáng)等特點(diǎn),更容易實(shí)現(xiàn)高通量、大規(guī)模的基因型檢測(cè)[9, 10],目前,利用SNP標(biāo)記對(duì)玉米雜種優(yōu)勢(shì)類群劃分也得以廣泛應(yīng)用[11-13]。

本研究采用SNP分子標(biāo)記對(duì)玉米品種新科891雙親的遺傳結(jié)構(gòu)和遺傳多樣性進(jìn)行分析,以期對(duì)其雙親自交系的雜種優(yōu)勢(shì)類群進(jìn)行劃分,進(jìn)而了解其雜種優(yōu)勢(shì)模式,為玉米種質(zhì)創(chuàng)新和遺傳改良提供理論依據(jù)。1材料與方法1.1試驗(yàn)材料

供試材料為玉米自交系S155、PHA458、A01(S155和PHA458分別為玉米品種新科891的母本和父本,PHA458是由PHB1M/A01雜交選育而成),新科891審定編號(hào)為國(guó)審玉20186115。以45份國(guó)內(nèi)常用玉米骨干自交系作為劃分雜種優(yōu)勢(shì)類群的參照(CK),見表1。

1.2SNP分子標(biāo)記分析

每份材料選擇10株葉片等面積混合取樣,采用改良的CTAB法提取DNA并進(jìn)行純化[14]。SNP分子標(biāo)記檢測(cè)基于北京中玉金標(biāo)記生物技術(shù)股份有限公司Affymetrix Axiom平臺(tái),采用10K玉米育種專用芯片—中玉芯1號(hào),對(duì)玉米自交系S155、PHA458、A01及45份國(guó)內(nèi)常用玉米骨干自交系進(jìn)行SNP分子標(biāo)記分析,該芯片有9 573個(gè)SNP標(biāo)記,標(biāo)記均勻分布10條染色體上,標(biāo)記間的物理距離平均都在200 kb左右。1.3數(shù)據(jù)統(tǒng)計(jì)分析

利用Axiom Analysis Suite V5.2進(jìn)行SNP分型分析,PowerMarker 3.25進(jìn)行最小等位基因頻率、樣品雜合率、缺失率和基因多樣性分析,Treebest 1.9.2構(gòu)建進(jìn)化樹,GCTA 1.26.0進(jìn)行主成分分析(PCA, Principal Component Analysis),Admixture 1.3.0進(jìn)行群體結(jié)構(gòu)分析。使用https://www.bioinformatics.com.cn在線網(wǎng)址繪制SNP密度圖。2結(jié)果與分析2.1SNP質(zhì)控

利用Maize 10K芯片對(duì)供試自交系進(jìn)行SNP基因型分析,位點(diǎn)質(zhì)控剩余8 932個(gè)標(biāo)記,占93.3%(dish QCgt;0.82、call rategt;97)。然后進(jìn)行過濾,缺失率(miss)lt;0.1,最小等位基因頻率(maf)gt;0.05,最終獲得4 926個(gè)可用標(biāo)記(圖1),占51.5%。48份玉米自交系的雜合率統(tǒng)計(jì)見表2,雜合率的計(jì)算方式是雜合標(biāo)記數(shù)除以總標(biāo)記數(shù)。由于48份玉米自交系是高代自交系,純合程度高,雜合率為0.003~0.009。

2.2群體聚類分析

利用Treebest的nj-tree模型構(gòu)建系統(tǒng)進(jìn)化樹,以45份國(guó)內(nèi)常用的玉米骨干自交系作為劃分雜種優(yōu)勢(shì)類群的參照,如圖2所示。48份自交系材料明顯分為六大類群,Ⅰ類群為歐美種質(zhì)群,包括PHA458、先玉335父本、S155、德美亞1號(hào)母本、德美亞2號(hào)母本、KWS9384母本、德美亞3號(hào)母本;Ⅱ類群為塘四平頭類血緣,包括德美亞1號(hào)父本、德美亞2號(hào)父本、掖81515、A01、吉單27母本、LX9801、K12、昌7-2、吉444、黃野四、黃早四;Ⅲ類群為蘭卡斯特類血緣,包括KWS9384父本、德美亞3號(hào)父本、獲唐黃、合344、吉002、Mo17、Mo17ht;Ⅳ類群為自330和旅大紅骨血緣,包括吉846、綜3、綜31、鄭58、吉853、C052、丹340、OH43、莫群17、自330;Ⅴ類群為改良瑞德類血緣,包括先玉335母本、掖478、丹9046、C8605-2、鐵7922、U8112、B73、鄭32;Ⅵ類群為PB群種質(zhì),包括豫自87-1、沈3336、X178、齊319、沈137。

2.3主成分分析

根據(jù)篩選的SNP結(jié)果,對(duì)48份玉米自交系材料進(jìn)行主成分分析(Principal component analysis, PCA)。PC1和PC2的二維主成分聚類圖見圖2,主成分PC1和PC2貢獻(xiàn)率分別為25.09%和18.37%,累積貢獻(xiàn)率43.46%(圖3)。根據(jù)48份自交系的集中程度,可清晰地鑒定出六大類群,這與群體聚類分析結(jié)果基本一致。

S155與先玉335母本、掖478、德美亞1號(hào)母本、德美亞2號(hào)母本、KWS9384母本、德美亞3號(hào)母本、X178、豫自87-1、齊319、沈137距離較近;PHA458與德美亞1號(hào)父本、德美亞2號(hào)父本、掖81515、X178、豫自87-1、齊319、沈137、A01距離較近;A01掖81515、吉單27母本距離較近。由于PHA458是由PHB1M/A01雜交選育而成,PHA458與A01距離較近。

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

利用Admixture軟件分析48份玉米自交系材料的群體遺傳結(jié)構(gòu),當(dāng)K=6時(shí)最適合,可分為六大類群(圖4),不同顏色表示不同的類群。第Ⅰ類群為自330和旅大紅骨血緣,包括鄭58、綜3、綜31、OH43、吉853、C052、丹340、自330、莫群17;第Ⅱ類群為改良瑞德類血緣,包括先玉335母本、C8605-2、掖478、丹9046、鄭32、U8112、B73、鐵7922;第Ⅲ類群為PB群,包括吉846、沈3336、沈137、豫自87-1、齊319、X178;第Ⅳ類群為蘭卡斯特類血緣,包括德美亞2號(hào)父本、KWS9384父本、德美亞3號(hào)父本、獲唐黃、Mo17ht、吉002、合344、Mo17;第Ⅴ類群為塘四平頭類血緣,包括德美亞1號(hào)父本、掖81515、A01、吉單27母本、黃野四、LX9801、黃早四、K12、吉444、昌7-2;第Ⅵ類群為歐美種質(zhì)群,包括S155、PHA458、先玉335父本、德美亞1號(hào)母本、KWS9384母本、德美亞3號(hào)母本、德美亞2號(hào)母本。類群劃分結(jié)果與Treebest軟件構(gòu)建系統(tǒng)進(jìn)化樹結(jié)果基本一致。

由不同顏色所占的比例(Q值)可以分析出S155、PHA458和A01在相應(yīng)類群內(nèi)所占的遺傳背景(表3)。根據(jù)劉麗華等[15]的研究結(jié)果,Q值gt;0.6視為血緣相對(duì)單一,Q值lt;0.6視為具有混合來源。S155和PHA458的歐美種質(zhì)血緣Q值較低,分別為0.318和0.354,遺傳背景比較復(fù)雜,A01的塘四平頭血緣Q值為0.574,雖然Q值接近0.6,但遺傳背景也為混合來源。

3結(jié)論與討論

本研究以45份國(guó)內(nèi)常用的玉米骨干自交系作為劃分雜種優(yōu)勢(shì)類群的參照,利用SNP芯片技術(shù),對(duì)玉米品種新科891的雙親及自交系A(chǔ)01進(jìn)行群體聚類分析和群體遺傳結(jié)構(gòu)分析。結(jié)果表明,48份自交系材料明顯分為六大類群,群體聚類分析和群體遺傳結(jié)構(gòu)分析結(jié)果基本一致。新科891的父本S155和母本PHA458被劃分到歐美種質(zhì)群,該群還包括先玉335父本、德美亞1號(hào)母本、KWS9384母本、德美亞3號(hào)母本、德美亞2號(hào)母本。德國(guó)KWS公司的德美亞系列早熟品種具有早熟、耐密、抗倒、籽粒收獲時(shí)含水量低等特點(diǎn),如德美亞1號(hào)、德美亞3號(hào)[16]。新科891的雙親與德國(guó)KWS公司系列品種的母本劃分到同一類群,可能是因?yàn)樾驴?91是早熟易機(jī)收品種,屬于易機(jī)收種質(zhì)。

群體改良是種質(zhì)創(chuàng)新的重要方式,也是主要方式,有利于群體遺傳多樣性的保持[17, 18],但是由于群體中含有不同類型種質(zhì)資源,遺傳背景就變得更加復(fù)雜。PN78599種質(zhì)自交系與Reid類群組成基礎(chǔ)選系材料,選育的二環(huán)系能提高Reid類群自交系的抗性和子粒品質(zhì)[19]。新科891的母本S155是由PN78599、B73、PH6WC等自交系組建的SS基礎(chǔ)群體選育而成,具有復(fù)雜的遺傳背景,這與本文群體遺傳結(jié)構(gòu)分析結(jié)果相符合。

自交系A(chǔ)01被劃分為塘四平頭類血緣,該群還包括掖81515、吉單27母本、LX9801、K12、昌7-2、吉444、黃野四、黃早四等。玉米自交系新01A3是以優(yōu)良玉米雜交種為基礎(chǔ)材料選育而成,具有中早熟、高抗、配合力高、品質(zhì)優(yōu)等突出優(yōu)點(diǎn)[20]。A01作為新01A3的姊妹系,屬于塘四平頭類種質(zhì)[21],這與本文分析結(jié)果相吻合。A01群體遺傳結(jié)構(gòu)分析Q值為0.574,略低于0.6,可能是由于A01是二環(huán)系選育而成,雖然被劃分為塘四平頭類血緣,但是遺傳背景比較復(fù)雜。父本PHA458是由PHB1M/A01雜交選育而成。PHB1M是先玉696的父本,為NSS類種質(zhì);因此PHA458含有NSS和塘四平頭類血緣,遺傳背景較為復(fù)雜,這與本文群體遺傳結(jié)構(gòu)分析結(jié)果相吻合。

美國(guó)的一些公益性機(jī)構(gòu)如艾奧瓦州立大學(xué)、北達(dá)科他州立大學(xué)等依據(jù)雜種優(yōu)勢(shì)模式,用不同來源的自交系或地方品種組建并改良了多個(gè)育種群體,如SS群體、NSS群體等,SS種質(zhì)主要來源于B73,PHB1M為Non-Reid類(NSS)種質(zhì)[22~24]。雖然本研究結(jié)果新科891的父本S155和母本PHA458被劃分到同一類群(歐美種質(zhì)群),但是根據(jù)品種選育來源,S155為SS類種質(zhì),PHA458含有NSS和塘四平頭類血緣,雜種優(yōu)勢(shì)模式為SS×NSS。“PN78599群×塘四平頭群”是黃淮海夏玉米區(qū)主要雜優(yōu)模式之一[25]。以PN78599和A01為基礎(chǔ)種質(zhì),選育出的玉米品種新科891,因具有產(chǎn)量高、抗逆性強(qiáng)、品質(zhì)優(yōu)良、綜合性狀好等突出特點(diǎn),已在生產(chǎn)上大面積推廣應(yīng)用[26]。通過SNP芯片技術(shù)對(duì)新科891的雙親進(jìn)行群體聚類分析和群體遺傳結(jié)構(gòu)分析,將為新科891雙親的后續(xù)改良與創(chuàng)新奠定基礎(chǔ),并為雜交組配的育種工作提供依據(jù)。參考文獻(xiàn):

[1]鄒成林, 譚華, 黃開健, 等. 基于SNP芯片劃分CIMMYT和廣西玉米自交系雜種優(yōu)勢(shì)群[J/OL].分子植物育種. https://kns.cnki.net/kcms/detail/46.1068.s.20220303.1730.022.html.

[2]段燦星, 董懷玉, 李曉, 等. 玉米種質(zhì)資源大規(guī)模多年多點(diǎn)多病害的自然發(fā)病抗性鑒定[J]. 作物學(xué)報(bào), 2020, 46(8): 1 135-1 145.

[3]黎裕, 王天宇. 我國(guó)玉米育種種質(zhì)基礎(chǔ)與骨干親本的形成[J]. 玉米科學(xué), 2010, 18(5): 1-8.

[4]趙久然, 李春輝, 宋偉, 等. 基于SNP芯片揭示中國(guó)玉米育種種質(zhì)的遺傳多樣性與群體遺傳結(jié)構(gòu)[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51(4): 626-634.

[5]張亞菲, 劉松濤, 曹雯梅, 等.黃淮海夏玉米品種主要性狀遺傳多樣性研究[J].種子, 2021, 40(4): 96-100.

[6]吳金鳳, 宋 偉, 王 蕊, 等. 利用SNP標(biāo)記對(duì)51份玉米自交系進(jìn)行類群劃分[J]. 玉米科學(xué), 2014, 22(5): 29-34.

[7]姜龍, 陳殿元, 周嵐, 等. 早熟類Iodent種質(zhì)與我國(guó)主要玉米類群的雜種優(yōu)勢(shì)分析[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020, 48(3):75-81.

[8]姚文華, 韓學(xué)莉, 汪燕芬, 等. 我國(guó)甜玉米育種研究現(xiàn)狀與發(fā)展對(duì)策[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào), 2011, 13(2): 1-8.

[9]唐立群, 肖層林, 王偉平. SNP分子標(biāo)記的研究及其應(yīng)用進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2012, 28(12): 154-158.

[10]劉傳光, 張桂權(quán). 水稻單核苷酸多態(tài)性及其應(yīng)用[J]. 遺傳, 2006(28): 737-744.

[11]盧柏山, 史亞興, 宋偉, 等. 利用SNP標(biāo)記劃分甜玉米自交系的雜種優(yōu)勢(shì)類群[J]. 玉米科學(xué), 2015, 23(1): 58-62, 68 .

[12]高嵩, 劉宏偉, 何歡, 等. 利用SNP芯片進(jìn)行玉米遺傳多樣性和群體遺傳結(jié)構(gòu)分析及新品種選育[J]. 玉米科學(xué), 2021, 29(1): 39-45.

[13]張鵬, 管俊嬌, 黃清梅, 等. 基于SNP芯片的云南玉米自交系遺傳多樣性和群體遺傳結(jié)構(gòu)分析[J]. 南方農(nóng)業(yè)學(xué)報(bào), 2020, 51(9): 2 082-2 089.

[14]藍(lán)碧秀, 王凜, 吳子愷, 等. 利用改良CTAB法快速小量提取微胚乳玉米基因組DNA[J]. 基因組學(xué)與應(yīng)用生物學(xué), 2015, 34(1): 190-194.

[15]劉麗華, 王立新, 趙昌平, 等. 光溫敏二系雜交小麥恢復(fù)系遺傳多樣性和群體結(jié)構(gòu)分析[J]. 中國(guó)生物化學(xué)與分子生物學(xué)報(bào), 2009, 25(9): 867-875.

[16]王元東, 趙久然, 張華生, 等. 玉米宜粒收品種的選育及種質(zhì)創(chuàng)新策略[J]. 分子植物育種, 2020, 18(10): 3 455-3 460.

[17]黎裕, 王天宇. 玉米種質(zhì)創(chuàng)新——進(jìn)展與展望[J]. 玉米科學(xué), 2017, 25(3): 11-18.

[18]李蘆江, 陳文生, 蘭海,等. 不同輪回選擇方法對(duì)玉米窄基群體遺傳多樣性的影響[J].核農(nóng)學(xué)報(bào), 2015, 29(1): 21-28.

[19]王秀鳳, 景希強(qiáng), 王孝杰, 等. PN78599種質(zhì)在我國(guó)玉米育種和生產(chǎn)中的應(yīng)用[J]. 玉米科學(xué), 2012, 20(4): 50-52.

[20]洪德峰, 張學(xué)舜, 馬毅, 等. 優(yōu)良玉米自交系新01A3及其改良系選育與應(yīng)用[J]. 玉米科學(xué), 2021, 29(1): 15-19.

[21]劉經(jīng)緯, 周聯(lián)東, 魏鋒, 等. 玉米自交系01A3和姊妹系A(chǔ)01、新6的選育及其衍生系在玉米雜交種選育中的利用[J]. 湖北農(nóng)業(yè)科學(xué), 2021, 60(17): 16-18.

[22]雍洪軍, 王建軍, 張德貴, 等. 美洲地區(qū)主要玉米群體特征及其利用途徑分析[J]. 遺傳, 2013, 35(6): 703-713.

[23]石雷. 現(xiàn)代美國(guó)馬齒型玉米商業(yè)育種的種質(zhì)基礎(chǔ)[J]. 玉米科學(xué), 2011, 19(5): 1-5.

[24]張艷輝, 欒奕, 王麗萍,等. 不同遺傳基礎(chǔ)玉米自交后代的誘導(dǎo)率和加倍率[J]. 吉林農(nóng)業(yè)大學(xué)學(xué)報(bào), 2017, 39(1): 10-14.

[25]李文才, 劉治先, 孫琦, 等. “PN78599群×塘四平頭群”玉米雜優(yōu)模式的利用及連續(xù)改良[J]. 山東農(nóng)業(yè)科學(xué), 2019, 51(6): 31-34.

[26]劉經(jīng)緯, 周聯(lián)東, 孫佩, 等. 多抗機(jī)收玉米新品種新科891的選育與應(yīng)用[J]. 中國(guó)種業(yè), 2019(3): 79-81.

主站蜘蛛池模板: 美女国产在线| 最新亚洲人成网站在线观看| 四虎永久在线视频| 欧美在线导航| 欧美人与动牲交a欧美精品| 一级一级一片免费| 久青草免费在线视频| 国产69精品久久| 9久久伊人精品综合| 国产不卡网| 欧美v在线| 成人精品免费视频| 欧美亚洲欧美| 精品国产毛片| 看看一级毛片| 国产尤物视频网址导航| 女人18毛片水真多国产| 国产青榴视频| 国产成人无码AV在线播放动漫 | 国产视频久久久久| 日韩一区精品视频一区二区| 怡红院美国分院一区二区| 自偷自拍三级全三级视频| 在线日韩一区二区| 四虎精品免费久久| 男女男免费视频网站国产| 国产第二十一页| 伊人大杳蕉中文无码| 日韩一区二区三免费高清| 亚洲永久免费网站| 欧美a级完整在线观看| 欧美区在线播放| 免费在线a视频| yy6080理论大片一级久久| 无码久看视频| 亚洲成a人片| 国产精品免费久久久久影院无码| 成人毛片免费在线观看| 亚洲色图另类| 国产精欧美一区二区三区| 午夜三级在线| 欧美 国产 人人视频| 欧美无专区| 国产aⅴ无码专区亚洲av综合网| 精品久久人人爽人人玩人人妻| 夜夜拍夜夜爽| 2021国产精品自产拍在线| 久久77777| 亚洲福利一区二区三区| 久久久久免费精品国产| 91久久精品国产| 亚洲欧美另类日本| 久久这里只有精品8| 尤物视频一区| 黄片在线永久| 看国产一级毛片| 91网在线| 97久久精品人人| 国产精品99久久久久久董美香| 91视频青青草| 国产一区二区三区夜色| 日韩专区欧美| 久久亚洲精少妇毛片午夜无码| 久久久久久午夜精品| 色婷婷综合激情视频免费看 | 欧美精品成人一区二区视频一| 国产欧美日韩专区发布| AV熟女乱| 在线精品欧美日韩| 91福利片| 伊人久久青草青青综合| 九九热在线视频| 日本手机在线视频| 色婷婷亚洲综合五月| 91人妻在线视频| 天天综合网在线| 91热爆在线| 999福利激情视频| 亚洲综合欧美在线一区在线播放| 2019国产在线| 亚洲有无码中文网| 国产综合色在线视频播放线视|