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

    基于元分析的大豆胞囊線蟲病3號生理小種抗性基因挖掘

    2020-05-11 03:11:04韓英鵬田利崢姜海鵬包冬芳
    東北農(nóng)業(yè)大學學報 2020年4期
    關鍵詞:一致性大豆

    韓英鵬,田利崢,姜海鵬,包冬芳,王 俊,趙 雪

    (東北農(nóng)業(yè)大學大豆研究所,大豆生物學教育部重點實驗室,農(nóng)業(yè)農(nóng)村部東北大豆生物學與遺傳育種重點實驗室,哈爾濱150030)

    大豆胞囊線蟲?。⊿oybean cyst nematode,SCN,Heterodera glycines Ichinohe)是大豆生產(chǎn)重要病害,主要分布于中國、美國、巴西、日本等國家。大豆胞囊線蟲病害分布廣泛,可造成20~30%產(chǎn)量損失,病害嚴重時高達80%,甚至絕產(chǎn),嚴重危害大豆生產(chǎn)[1]。

    大豆胞囊線蟲病生理小種分化明顯,Riggs 等利用大豆4個鑒別寄主和1個對照感病寄主調(diào)查胞囊線蟲反應,將其劃分為16 個生理小種[2],根據(jù)7個鑒別寄主區(qū)劃分為HG 類型[3]。目前,我國大豆主產(chǎn)區(qū)發(fā)現(xiàn)1~7 號、9 號、11 號、13 號、14 號等11個SCN生理小種,其中3號生理小種在我國分布最廣,主要分布在東北三省和內(nèi)蒙古;4號生理小種浸染能力更強,主要分布在黃淮和西北地區(qū);6號生理小種主要發(fā)現(xiàn)于黑龍江省內(nèi)[4]??共∑贩N(系)選育是防治SCN 最有效方法,挖掘SCN 抗性基因是分子輔助育種基礎,而QTL 定位是挖掘抗性基因最有效方法之一。

    Weismann 等首次將抗性基因Rhg4 定位在pBLT24 與pBLT56 區(qū)間[5]。Concibido 定位多個SCN抗性QTL,且整合1994~2004年發(fā)表文獻報道的60多個SCN 抗性QTL,分布在除2、9、10 等以外17個染色體[6]。目前,國內(nèi)外文獻已報道大量不同SCN生理小種抗性QTL,但因遺傳抗性材料和接種方法不同,且不同地理環(huán)境和遺傳背景對QTL 穩(wěn)定性影響較大,導致QTL 基因組區(qū)域與實際所在基因組區(qū)域不完全相同。

    Goffinet 建立可整合不同遺傳背景和不同材料一致性QTL 元分析方法,基于最大鄰接法整合不同來源QTL,縮小置信區(qū)間,提高QTL 定位結果有效性和準確性[7]。通過meta分析在QTL區(qū)段發(fā)現(xiàn)相關性狀候選基因在大豆和其他作物中被廣泛利用。李長育等利用BioMercatorV2.1軟件通過meta分析獲得2個一致性QTL位點,得到與大豆根瘤相關的8 個候選基因[8]。在玉米中,通過meta 分析方法,在847 個有關玉米籽粒蛋白質(zhì)含量候選基因中,發(fā)現(xiàn)13個影響該性狀的關鍵基因[9]。

    大豆SCN抗性研究已定位和鑒定160多個QTL(http://www.soybase.or),其中,一些QTL未經(jīng)重復驗證,不同試驗條件下,鑒定結果不一致。目前,已克隆和鑒定2 個主要SCN 抗性位點rhg1 和Rhg4,并揭示新的大豆抗病機制。目前已有栽培品種對SCN 抗 性 主 要 來 自PI88788、PI437654 和Peking,SCN 在不斷適應過程中,易克服大豆寄主單一抗性,因此需鑒定新的SCN抗性位點和基因。

    本研究使用圖譜整合軟件BioMercatorV4.2 將不同來源SCN Race 3 抗性QTL 整合到遺傳圖譜GmComposite 2003,并通過元分析方法縮小QTL置信區(qū)間,進一步提高QTL 準確性并獲得一致性QTL,選取區(qū)間內(nèi)所有基因并作基因注釋,進一步篩選與SCN Race 3 相關抗性基因,分析其抗病機理,最終為分子輔助育種尤其選育抗SCN 品種(系)提供種質(zhì)資源和理論依據(jù)。

    1 材料與方法

    1.1 大豆抗胞囊線蟲QTL信息收集與整理

    收集已發(fā)表的SCN Race 3 抗性QTL 信息,包括QTL 位置、置信區(qū)間、LOD 值、遺傳貢獻率、作圖群體類型、群體規(guī)模和臨近標記等信息。其中,我國品種(系)有中品03-5373 和L-10,中品03-5373 對SCN Race 3 抗性來自Peking、PI437654和灰皮支黑豆等3 個抗源,L-10 對9 個胞囊線蟲小種均有抗性;美國品種(系)有PI437655、PI438489B、PI567516C、 PI89772、 Bell、 Msoy8001 和PI437654,均具有多種胞囊線蟲小種抗性,其中PI437655、 PI438489B、 PI567516C、 PI89772 和PI437654等均為野生型大豆。

    1.2 SCN Race 3抗性QTL映射與meta-QTL分析

    利用BioMercatorV4.2 軟件完成各作圖群體SCN Race 3 抗性QTL 映射及meta-QTL 分析?;诿總€QTL 圖譜名稱(Map name)、QTL 所在染色體(Chromosome)、群 體 規(guī) 模(Size)、群 體 類 型(Type)、QTL 名稱(QTL name)、QTL 位置(Position)、遺傳貢獻率(R2)、LOD 值和各原始圖譜與公共圖譜上相關標記,構建一致性遺傳圖譜(Consensus map),隨后利用ConsMap 和QTLProj 功能將每個原始圖譜中QTL 映射到GmComposite2003 上。利用Veyrieras元分析作meta-QTL分析。

    meta-QTL 分析后得到5 個模型,每個模型均按照最大似然函數(shù)比方法,利用高斯定理推導出在染色體上最可能排列位置和置信區(qū)間。根據(jù)AIC(Akaike-type-criteria-values)最小值原則,當在每條染色體上選取meta-QTL 數(shù)目時,首先考慮AIC 最小值對應的meta-QTL 數(shù)目模型,參考一致性圖譜上原始QTL 分布及數(shù)目,選取合適數(shù)目,獲取一致性置信區(qū)間。

    1.3 SCN Race3 抗性meta-QTL 區(qū)段候選基因發(fā)掘及分析

    利用已獲得的QTL 一致性置信區(qū)間,根據(jù)基因組預測注釋信息,挖掘SCN Race 3 抗性相關基因。查找已克隆的與大豆胞囊線蟲病抗性相關基因,下載其表達序列,分析其結構域特征。將本研究所獲置信區(qū)間內(nèi)基因作為候選基因,同時根據(jù)其包含的不同功能結構域,將抗性基因分類,將所得抗性基因與大豆EST 數(shù)據(jù)庫比對、注釋。

    1.4 SCN Race 3抗性相關基因蛋白結構域分析

    在Phytozome網(wǎng)站(https://phytozome.jgi.doe.gov/)利用Glycine max Williams82 基因組數(shù)據(jù)庫,獲得基因參考序列,通過網(wǎng)站Functional Annotation 工具分析蛋白結構域,并使用繪圖軟件IBS繪圖。

    2 結果與分析

    2.1 SCN Race 3抗性QTL信息收集與整理

    分析收集的SCN Race 3 抗性QTL 信息可見,共涉及7 個作圖群體,獲得33 個SCN Race 3 相關QTL 信息(見表1),分布于11 個連鎖群,其中G(chr18)分布較多,為5 個;F(chr13)、I(chr20)較少,均為1個。

    表1 大豆胞囊線蟲病3號小種抗性QTL相關信息Table 1 Information related to QTL resistance of soybean cyst nematode race 3 published

    2.2 SCN Race 3抗性QTL映射與meta-QTL分析

    共有21 個報道成功的QTL 映射到大豆參考圖譜,利用BioMercatorV4.2 軟件中的QTL meta-analyses程序分析計算,每次運算后以AIC最小值為原則,根據(jù)QTL模型得到3個一致性QTL,分別位于M(chr7)、E(chr15)和G(chr18)(見表2、圖1)等3個連鎖群中,其中置信區(qū)間指數(shù)(CI)為95%。這些一致性QTL 可為分子標記輔助育種、QTL 精細定位、挖掘抗性候選提供重要位置信息。

    2.3 SCN Race 3 抗性meta-QTL 區(qū)段候選基因發(fā)掘及分析

    挖掘利用BioMercatorV4.2 獲得的一致性置信區(qū)間內(nèi)候選基因,根據(jù)基因組測序注釋信息,共發(fā)現(xiàn)18 個SCN Race 3 抗性相關基因(見表3),17個蛋白激酶(PK)型基因和1 個亮氨酸重復序列類受體蛋白激酶(LRR)型基因。

    表2 候選QTL位點Table 2 Candidate QTL loci

    表3 抗病基因類型統(tǒng)計Table 3 Types statistics of resistance genes

    2.4 SCN Race 3相關抗性基因蛋白結構域分析

    通過Functional Annotation工具分析基因結構域,發(fā)現(xiàn)17個PK型基因分為絲/蘇氨酸蛋白激酶、半胱氨酸蛋白激酶、組/精/賴氨酸蛋白激酶3類(見圖2)。

    絲/蘇氨酸蛋白激酶包括Glyma.18G074200、Glyma.18G096500、Glyma.18G026700、Glyma.18G0 26800、Glyma.18G026900、Glyma.18G062500、Glyma.07G003200、 Glyma.18G060900、 Glyma.18G042200、Glyma.18G038500 和Glyma.07G021100, 其 中Glyma.18G074200和Glyma.18G096500具有鈣依賴蛋白激酶活性;Glyma.07G003200 和Glyma.18G060900具有分裂原活化蛋白激酶活性,該結構域參與大豆PTI和ETI兩種防御機制,從而影響大豆對胞囊線蟲抗性。半胱氨酸蛋白激酶包括Glyma.18G0 46100、Glyma.18G046200、 Glyma.18G046300、 Glyma.8G0 46600 和Glyma.18G046400,該結構域通過誘導大豆PTI防御機制從而影響大豆對胞囊線蟲抗性;組/精/賴氨酸蛋白激酶包括Glyma.18G054100;Glyma.15G155600 為LRR 型基因,具有亮氨酸重復序列類受體蛋白激酶活性,該結構域通過參與大豆ETI防御機制從而影響大豆對胞囊線蟲的抗性。

    3 討 論

    利用精細定位方法可從已定位QTL 中準確確定候選基因,但耗時耗力。利用元分析方法獲得QTL一個較小置信區(qū)間,將目的基因鎖定在較小置信區(qū)間內(nèi),再根據(jù)抗病基因保守性挖掘目標基因,最后利用EST 注釋信息進一步確認所獲基因,該方法省時省力且費用極低。?,|等利用BioMercator2.1 軟件,通過已知SCN 抗性QTL 和參考圖譜GmComposite2003獲得一致性位點,根據(jù)大豆20個連鎖標記(SSR、SNP 和RFLP)序列和大豆全基因組數(shù)據(jù)構建物理圖譜,與已獲得一致性QTL 位點,根據(jù)hmmsearch 搜索結果獲得77 個抗性基因[20]。在77個候選基因中,多數(shù)基因已證明與SCN抗性有關。因此,通過元分析計算確定一致性QTL位點從而準確篩選候選基因,為大豆分子輔助育種尤其是選育抗SCN 品種(系)提供一種新的思路和途徑。

    Li等研究中品03-5373(抗性)與中黃13(感性)及其后代所構成群體,根據(jù)其基因型和表型,經(jīng)全基因組關聯(lián)分析(GWAS),發(fā)現(xiàn)基因Glyma.18G026900具有SCN 抗性連鎖的單核苷酸多態(tài)性(SNP)位點(C/T)[21]。因此,可推斷其同源基因Glyma.18G026700、Glyma.18G026800 和Glyma.18G062500 等也可能具有SCN抗性連鎖的單核苷酸多態(tài)性(SNP)位點。但其余15個有關SCN抗性基因尚無報道。

    Glyma.18G026900及其同源基因具有絲氨酸/蘇氨酸蛋白激酶活性,屬于PK 酶(蛋白激酶)的一種。植物蛋白激酶在逆境脅迫信號傳導中發(fā)揮關鍵作用,維持植物體正常生命活動,根據(jù)底物特異性,植物蛋白激酶分為絲/蘇氨酸蛋白激酶、組/精/賴氨酸蛋白激酶、酪氨酸蛋白激酶、半胱氨酸蛋白激酶、天冬酰胺基/谷氨酰胺基蛋白激酶等5類[22]。其中絲/蘇氨酸蛋白激酶包括分裂原活化蛋白激酶(MAPK)和鈣依賴蛋白激酶(CDPK):當大豆受體蛋白參與植物體PTI 防御機制時,MAPK 在其下游的防御機制中發(fā)揮重要作用,植物體特異抗病蛋白,用于識別效應蛋白,即防御機制ETI,大豆通過MAPK在大豆胞囊線蟲防御機制ETI中發(fā)揮重要作用[23];CDPK 主要調(diào)節(jié)植物體中鹽信號傳導,其在大豆胞囊線蟲中抗性機理還需進一步研究。半胱氨酸類受體激酶(Cysteine-rich receptorlike-kinases,CRKs)屬于蛋白激酶,在植物體受病原菌脅迫時,CRK 誘導分泌甲基SA 誘導植物體PTI防御機制,使植物局部過敏性壞死,保護植物體[24]。其余蛋白激酶在植物非生物脅迫中發(fā)揮作用,但在生物脅迫中,尤其在SCN抗性中尚無報道。

    LRR 型蛋白因其結構域中富含亮氨酸重復序列而得名。參與植物體防御反應的LRR 型蛋白可分為類受體蛋白激酶、抗病基因編碼蛋白質(zhì)、多聚半乳糖醛酸酶抑制蛋白和伸展蛋白等4類[25]。類受體蛋白激酶主要參與植物體PTI防御過程;植物抗病基因編碼蛋白質(zhì)主要參與植物體ETI 過程[26];多聚半乳糖醛酸酶抑制蛋白(Polygalac-turonase-inhibiting proteins,PGIPs)通過保護植物細胞壁抵御病原菌入侵;伸展蛋白含有組成植物體細胞壁的結構蛋白,是防止病原物侵入和擴散的屏障。Kanazin 等在16 號染色體上定位到一個具有LRRTM 結構域的基因,已證實其與大豆胞囊線蟲抗性相關[27]。

    本研究在整合前人QTL 信息基礎上,采用圖譜整合軟件BioMercatorV4.2 挖掘SCN Race 3 抗性基因,將收集和整理的33 個QTL 整合在遺傳圖譜GmComposite2003上,并通過元分析方法得到一致性QTL,獲得3 個抗性候選區(qū)段,發(fā)現(xiàn)18 個與SCN Race 3 抗性相關基因。本研究選取的SCN Race 3抗性QTL信息及所構建遺傳圖譜描述較為詳細,因此可使用meta 分析方法發(fā)現(xiàn)與SCN Race 3抗性相關基因,對發(fā)現(xiàn)其他SCN 生理小種抗性基因具有借鑒意義。該方法為分子輔助育種提供新的思路和途徑,也為后期篩選候選基因,分析其抗病機理,最終為選育抗SCN 品種(系)提供種質(zhì)資源和理論基礎。

    4 結 論

    采集美國、中國和巴西等發(fā)表的33 個與SCN Race 3抗性相關QTL及其相關信息,利用QTL-meta分析得到3 個一致性QTL,從中發(fā)現(xiàn)18 個與SCN Race 3 抗性相關基因,其中17 個PK 型基因,1 個LRR 型基因。17 個PK 型基因位于連鎖群M(chr7)和連鎖群G(chr18),1個LRR型基因位于連鎖群E(chr15)。通過蛋白結構域分析,在17 個PK 基因中,1個基因具有組/精/賴氨酸蛋白激酶活性、5個基因具有半胱氨酸蛋白激酶活性、11 個基因具有絲/蘇氨酸蛋白激酶活性。

    猜你喜歡
    一致性大豆
    注意防治大豆點蜂緣蝽
    大豆帶狀種植技術如何提升我國大豆產(chǎn)量
    從大豆種植面積增長看我國糧食安全
    關注減污降碳協(xié)同的一致性和整體性
    公民與法治(2022年5期)2022-07-29 00:47:28
    巴西大豆播種順利
    注重教、學、評一致性 提高一輪復習效率
    對歷史課堂教、學、評一體化(一致性)的幾點探討
    大豆的營養(yǎng)成分及其保健作用
    IOl-master 700和Pentacam測量Kappa角一致性分析
    用大豆“炸沉”軍艦
    亚洲精品日本国产第一区| 久久久欧美国产精品| 在线亚洲精品国产二区图片欧美| 日本av免费视频播放| 少妇人妻久久综合中文| 久久综合国产亚洲精品| 侵犯人妻中文字幕一二三四区| 超碰97精品在线观看| 国产在线观看jvid| 国产欧美亚洲国产| 亚洲人成网站在线观看播放| 亚洲熟女精品中文字幕| 黑丝袜美女国产一区| 国产亚洲午夜精品一区二区久久| 18禁国产床啪视频网站| av线在线观看网站| 热99国产精品久久久久久7| 国产免费视频播放在线视频| 每晚都被弄得嗷嗷叫到高潮| 国产免费又黄又爽又色| cao死你这个sao货| 国产精品三级大全| 777久久人妻少妇嫩草av网站| 男女午夜视频在线观看| 高清视频免费观看一区二区| 在线观看免费视频网站a站| 男女午夜视频在线观看| 美女高潮到喷水免费观看| 丝袜人妻中文字幕| 久久综合国产亚洲精品| 欧美日韩亚洲高清精品| 大香蕉久久网| 少妇人妻久久综合中文| 国产精品久久久久成人av| 99久久人妻综合| 看免费成人av毛片| 视频区图区小说| 国产亚洲精品第一综合不卡| 啦啦啦 在线观看视频| 欧美日韩成人在线一区二区| 色视频在线一区二区三区| 亚洲成人免费电影在线观看 | 久久人妻福利社区极品人妻图片 | 国产日韩欧美在线精品| 国产女主播在线喷水免费视频网站| 日韩中文字幕视频在线看片| 男人爽女人下面视频在线观看| tube8黄色片| 成人亚洲欧美一区二区av| 亚洲人成电影免费在线| 99久久综合免费| 久久人人97超碰香蕉20202| 亚洲精品美女久久av网站| 久久99精品国语久久久| 看免费av毛片| 欧美成人精品欧美一级黄| 人人妻人人爽人人添夜夜欢视频| 午夜av观看不卡| 久久影院123| 亚洲精品第二区| 欧美黑人精品巨大| 只有这里有精品99| 免费在线观看日本一区| 亚洲av在线观看美女高潮| 性高湖久久久久久久久免费观看| 日本vs欧美在线观看视频| 精品高清国产在线一区| 丁香六月天网| 日本a在线网址| 精品久久蜜臀av无| 九色亚洲精品在线播放| 人妻一区二区av| 一区二区三区精品91| 制服人妻中文乱码| 女人被躁到高潮嗷嗷叫费观| 天天添夜夜摸| 成人三级做爰电影| 欧美成人午夜精品| 国产一区二区三区综合在线观看| 日本wwww免费看| 热re99久久国产66热| 91字幕亚洲| 国产亚洲欧美在线一区二区| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产av成人精品| av在线老鸭窝| 1024视频免费在线观看| 亚洲欧美日韩高清在线视频 | 亚洲精品第二区| 午夜免费男女啪啪视频观看| 午夜av观看不卡| av在线播放精品| 亚洲精品国产av蜜桃| 国产亚洲欧美在线一区二区| 久久久久视频综合| 国产又色又爽无遮挡免| 中国美女看黄片| 久久久精品94久久精品| 免费在线观看影片大全网站 | av网站在线播放免费| 在线 av 中文字幕| 亚洲精品一二三| 91国产中文字幕| 亚洲欧美日韩高清在线视频 | 日韩中文字幕欧美一区二区 | av电影中文网址| 国产精品久久久av美女十八| 欧美xxⅹ黑人| 亚洲欧洲精品一区二区精品久久久| 欧美成人午夜精品| 午夜福利免费观看在线| 午夜激情久久久久久久| 建设人人有责人人尽责人人享有的| 亚洲男人天堂网一区| 观看av在线不卡| 91精品伊人久久大香线蕉| 极品少妇高潮喷水抽搐| 久久鲁丝午夜福利片| 男女免费视频国产| 麻豆av在线久日| 中文字幕亚洲精品专区| www日本在线高清视频| 午夜精品国产一区二区电影| 蜜桃在线观看..| 国产av精品麻豆| 视频区图区小说| 丰满迷人的少妇在线观看| 婷婷色综合www| 久久久久久免费高清国产稀缺| 国产一区二区三区综合在线观看| xxx大片免费视频| 在现免费观看毛片| 韩国精品一区二区三区| 亚洲伊人久久精品综合| 色婷婷av一区二区三区视频| 最黄视频免费看| 美女扒开内裤让男人捅视频| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三区在线| 成年美女黄网站色视频大全免费| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 久久天躁狠狠躁夜夜2o2o | 两个人免费观看高清视频| 国产精品偷伦视频观看了| 国产1区2区3区精品| videos熟女内射| 美女大奶头黄色视频| 黄色怎么调成土黄色| 日本五十路高清| 永久免费av网站大全| 中国国产av一级| 自线自在国产av| 成年美女黄网站色视频大全免费| 看免费av毛片| 在线观看人妻少妇| 亚洲欧美一区二区三区黑人| 国产极品粉嫩免费观看在线| 97精品久久久久久久久久精品| 老司机亚洲免费影院| 丝袜在线中文字幕| 亚洲,欧美,日韩| 亚洲美女黄色视频免费看| 久久久国产一区二区| 人妻人人澡人人爽人人| 三上悠亚av全集在线观看| 久久av网站| 麻豆av在线久日| 色综合欧美亚洲国产小说| 一级毛片我不卡| 18禁观看日本| 人人妻人人爽人人添夜夜欢视频| 久久精品成人免费网站| 国产高清国产精品国产三级| 欧美亚洲日本最大视频资源| 国产在线视频一区二区| 日韩大码丰满熟妇| 男女高潮啪啪啪动态图| 天堂俺去俺来也www色官网| 日韩电影二区| 国产精品一区二区在线观看99| 在线观看免费日韩欧美大片| 欧美中文综合在线视频| 老司机亚洲免费影院| av在线app专区| 欧美激情高清一区二区三区| e午夜精品久久久久久久| 一本综合久久免费| 秋霞在线观看毛片| 精品欧美一区二区三区在线| 丰满迷人的少妇在线观看| 97人妻天天添夜夜摸| 精品人妻1区二区| 久久女婷五月综合色啪小说| 18禁黄网站禁片午夜丰满| 亚洲色图综合在线观看| 悠悠久久av| 精品人妻熟女毛片av久久网站| 91国产中文字幕| 色视频在线一区二区三区| 欧美亚洲日本最大视频资源| 少妇被粗大的猛进出69影院| 国产精品久久久久久人妻精品电影 | 国产女主播在线喷水免费视频网站| 成年女人毛片免费观看观看9 | 亚洲精品第二区| 久久久欧美国产精品| 久久性视频一级片| 青春草亚洲视频在线观看| 亚洲一区中文字幕在线| 日韩视频在线欧美| 成年美女黄网站色视频大全免费| 国产一区亚洲一区在线观看| 超碰成人久久| 国产日韩欧美在线精品| 国产成人精品无人区| 2018国产大陆天天弄谢| www日本在线高清视频| 日韩电影二区| 中国美女看黄片| 成人影院久久| 一区二区三区激情视频| 国产欧美日韩精品亚洲av| 一区在线观看完整版| 91精品国产国语对白视频| 天堂俺去俺来也www色官网| av欧美777| 久久鲁丝午夜福利片| 香蕉国产在线看| 夜夜骑夜夜射夜夜干| 亚洲成av片中文字幕在线观看| 一本色道久久久久久精品综合| 国产亚洲av高清不卡| 成人亚洲欧美一区二区av| 日韩制服骚丝袜av| 亚洲欧美一区二区三区久久| 国产在线视频一区二区| 99热网站在线观看| 亚洲欧美一区二区三区国产| 国产成人精品久久久久久| 日日夜夜操网爽| avwww免费| 国产成人一区二区三区免费视频网站 | videos熟女内射| 香蕉国产在线看| 亚洲国产精品成人久久小说| 香蕉丝袜av| 一级毛片女人18水好多 | 国产人伦9x9x在线观看| 99久久综合免费| 午夜激情久久久久久久| 一级毛片女人18水好多 | 国产成人精品在线电影| e午夜精品久久久久久久| 欧美人与性动交α欧美软件| 精品国产一区二区三区四区第35| 满18在线观看网站| 国产精品 欧美亚洲| 韩国高清视频一区二区三区| 亚洲精品美女久久av网站| 国产精品秋霞免费鲁丝片| 午夜福利影视在线免费观看| 日本欧美视频一区| 一级片免费观看大全| 久久午夜综合久久蜜桃| 国产欧美日韩一区二区三区在线| 人体艺术视频欧美日本| 国产亚洲欧美精品永久| 97精品久久久久久久久久精品| 精品人妻熟女毛片av久久网站| 亚洲国产av影院在线观看| 香蕉国产在线看| 中文精品一卡2卡3卡4更新| 高清不卡的av网站| 热re99久久精品国产66热6| 精品视频人人做人人爽| 香蕉国产在线看| 久久国产精品人妻蜜桃| 久久久久国产精品人妻一区二区| 精品一区二区三卡| 久久精品国产亚洲av高清一级| 中文精品一卡2卡3卡4更新| 亚洲av美国av| 丝袜人妻中文字幕| 亚洲成国产人片在线观看| 亚洲中文字幕日韩| 欧美黑人欧美精品刺激| 日本av免费视频播放| 久久综合国产亚洲精品| 国产精品久久久久久精品古装| 亚洲精品一二三| 最新的欧美精品一区二区| av一本久久久久| 啦啦啦啦在线视频资源| 女性生殖器流出的白浆| 最新在线观看一区二区三区 | 亚洲图色成人| 国产亚洲精品久久久久5区| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 国产亚洲欧美在线一区二区| 日韩中文字幕视频在线看片| 少妇 在线观看| 国产一区二区三区综合在线观看| av片东京热男人的天堂| 水蜜桃什么品种好| 成年人免费黄色播放视频| 国产成人91sexporn| 亚洲午夜精品一区,二区,三区| a级片在线免费高清观看视频| 啦啦啦在线观看免费高清www| 亚洲伊人色综图| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频| 久久精品久久久久久久性| 久久综合国产亚洲精品| 欧美国产精品va在线观看不卡| av欧美777| 十八禁人妻一区二区| 久久ye,这里只有精品| 欧美精品一区二区免费开放| 欧美国产精品一级二级三级| 亚洲三区欧美一区| 啦啦啦 在线观看视频| 亚洲免费av在线视频| 午夜精品国产一区二区电影| 又大又黄又爽视频免费| 国产爽快片一区二区三区| 欧美成狂野欧美在线观看| 侵犯人妻中文字幕一二三四区| 天天操日日干夜夜撸| 中文字幕精品免费在线观看视频| 性色av乱码一区二区三区2| 9191精品国产免费久久| 欧美日韩综合久久久久久| 热re99久久国产66热| 黄色一级大片看看| 热re99久久国产66热| 美女高潮到喷水免费观看| 脱女人内裤的视频| 欧美激情极品国产一区二区三区| 丁香六月欧美| 视频区欧美日本亚洲| 国产av国产精品国产| 一二三四社区在线视频社区8| 黑人欧美特级aaaaaa片| 久久久久精品人妻al黑| 中文字幕av电影在线播放| 国产成人一区二区在线| 高潮久久久久久久久久久不卡| 久久亚洲国产成人精品v| 国产高清国产精品国产三级| 一本一本久久a久久精品综合妖精| 国产精品一区二区在线不卡| 久久av网站| 欧美久久黑人一区二区| 热99久久久久精品小说推荐| 两个人看的免费小视频| 深夜精品福利| 免费久久久久久久精品成人欧美视频| 久热这里只有精品99| 免费久久久久久久精品成人欧美视频| 少妇被粗大的猛进出69影院| 两个人看的免费小视频| 日韩,欧美,国产一区二区三区| 脱女人内裤的视频| 人人妻人人澡人人爽人人夜夜| 日本av手机在线免费观看| 一本一本久久a久久精品综合妖精| 欧美久久黑人一区二区| 啦啦啦 在线观看视频| 亚洲 欧美一区二区三区| 看免费成人av毛片| 在线观看免费高清a一片| videosex国产| 一区二区三区乱码不卡18| 国产熟女欧美一区二区| 成年人午夜在线观看视频| 永久免费av网站大全| 亚洲欧美中文字幕日韩二区| 18在线观看网站| 黄网站色视频无遮挡免费观看| 久久精品亚洲av国产电影网| 中文字幕亚洲精品专区| 亚洲伊人久久精品综合| 国产又爽黄色视频| 国产精品 欧美亚洲| 热99久久久久精品小说推荐| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲欧美在线一区二区| 91九色精品人成在线观看| 国产av精品麻豆| 我要看黄色一级片免费的| 男女免费视频国产| 精品少妇一区二区三区视频日本电影| 制服诱惑二区| 纯流量卡能插随身wifi吗| 久久精品久久久久久噜噜老黄| 亚洲 国产 在线| 一区二区av电影网| 日本91视频免费播放| 亚洲精品乱久久久久久| 人妻 亚洲 视频| 中文字幕最新亚洲高清| 日本黄色日本黄色录像| netflix在线观看网站| 国产亚洲精品久久久久5区| 国产老妇伦熟女老妇高清| 90打野战视频偷拍视频| 99久久精品国产亚洲精品| 欧美日韩视频精品一区| 热re99久久国产66热| 久久亚洲精品不卡| 亚洲熟女精品中文字幕| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产最新在线播放| 看免费成人av毛片| 国产精品三级大全| videos熟女内射| 日本色播在线视频| 又紧又爽又黄一区二区| 亚洲欧美精品自产自拍| 国产成人精品久久二区二区免费| 亚洲国产精品999| 欧美中文综合在线视频| 午夜福利视频在线观看免费| 免费观看a级毛片全部| 亚洲五月色婷婷综合| 纵有疾风起免费观看全集完整版| 午夜福利影视在线免费观看| 亚洲av日韩精品久久久久久密 | 亚洲国产欧美日韩在线播放| 极品人妻少妇av视频| 欧美乱码精品一区二区三区| 黑丝袜美女国产一区| 久久精品亚洲熟妇少妇任你| 又大又爽又粗| 少妇 在线观看| 久久天躁狠狠躁夜夜2o2o | 久久国产精品大桥未久av| 久久久久久久精品精品| av欧美777| 国产人伦9x9x在线观看| 国产一区二区激情短视频 | 超碰97精品在线观看| 日本色播在线视频| 日本91视频免费播放| 午夜久久久在线观看| 午夜福利,免费看| 久久 成人 亚洲| 欧美激情高清一区二区三区| 一二三四在线观看免费中文在| 午夜av观看不卡| 久久精品国产a三级三级三级| 十分钟在线观看高清视频www| 黄色毛片三级朝国网站| 精品少妇黑人巨大在线播放| 99久久99久久久精品蜜桃| 男女免费视频国产| 亚洲精品在线美女| 久久这里只有精品19| 国产真人三级小视频在线观看| 国产男人的电影天堂91| 宅男免费午夜| 精品卡一卡二卡四卡免费| 高清黄色对白视频在线免费看| 亚洲欧美激情在线| 18禁黄网站禁片午夜丰满| 在线观看一区二区三区激情| 久久影院123| 女性生殖器流出的白浆| 国产片特级美女逼逼视频| 在线观看免费日韩欧美大片| 亚洲国产精品国产精品| 欧美另类一区| 在线天堂中文资源库| 婷婷色综合www| 亚洲一区二区三区欧美精品| 男人舔女人的私密视频| 亚洲国产欧美在线一区| 人妻一区二区av| netflix在线观看网站| 伊人亚洲综合成人网| 亚洲精品中文字幕在线视频| 日韩大片免费观看网站| a级片在线免费高清观看视频| 真人做人爱边吃奶动态| 国产精品久久久久久人妻精品电影 | 亚洲精品国产一区二区精华液| 高清视频免费观看一区二区| 国产在视频线精品| 免费高清在线观看视频在线观看| 精品亚洲乱码少妇综合久久| 欧美另类一区| 成人亚洲精品一区在线观看| 午夜激情av网站| 成人黄色视频免费在线看| 亚洲成人免费av在线播放| 精品一区在线观看国产| 人人妻人人澡人人爽人人夜夜| 手机成人av网站| 久9热在线精品视频| 欧美av亚洲av综合av国产av| 久久人人爽av亚洲精品天堂| 欧美激情高清一区二区三区| 激情视频va一区二区三区| 老司机午夜十八禁免费视频| 日韩欧美一区视频在线观看| 日日摸夜夜添夜夜爱| 亚洲一码二码三码区别大吗| 精品国产超薄肉色丝袜足j| 日本vs欧美在线观看视频| 午夜福利一区二区在线看| 亚洲欧美日韩高清在线视频 | videosex国产| 亚洲精品一卡2卡三卡4卡5卡 | 国产男女内射视频| 大香蕉久久网| 午夜激情av网站| 又大又黄又爽视频免费| 日韩精品免费视频一区二区三区| 日日爽夜夜爽网站| 欧美精品人与动牲交sv欧美| 亚洲专区中文字幕在线| 国产精品 欧美亚洲| 性高湖久久久久久久久免费观看| 国产精品99久久99久久久不卡| 久久精品久久久久久噜噜老黄| 国产淫语在线视频| 免费一级毛片在线播放高清视频 | 性色av乱码一区二区三区2| 性色av一级| 一本色道久久久久久精品综合| 亚洲,欧美精品.| 一级毛片我不卡| 免费高清在线观看视频在线观看| 精品亚洲成国产av| 久久国产精品人妻蜜桃| 少妇裸体淫交视频免费看高清 | 欧美在线黄色| 久久精品久久精品一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲专区中文字幕在线| 激情五月婷婷亚洲| 欧美日韩综合久久久久久| 午夜两性在线视频| 一本—道久久a久久精品蜜桃钙片| 人妻一区二区av| 99国产精品一区二区三区| 免费观看a级毛片全部| 久久狼人影院| 真人做人爱边吃奶动态| 亚洲一区二区三区欧美精品| 国产日韩欧美亚洲二区| 欧美日韩视频高清一区二区三区二| 精品免费久久久久久久清纯 | 各种免费的搞黄视频| www.999成人在线观看| 老司机影院成人| 69精品国产乱码久久久| av欧美777| 亚洲av片天天在线观看| 大香蕉久久成人网| 国产xxxxx性猛交| 亚洲,一卡二卡三卡| 国产精品 国内视频| 精品人妻1区二区| 亚洲情色 制服丝袜| 在线看a的网站| 国产欧美日韩一区二区三区在线| 欧美另类一区| 欧美日韩亚洲综合一区二区三区_| 欧美日韩黄片免| 国产男女内射视频| 免费高清在线观看日韩| 国产精品一二三区在线看| 精品久久蜜臀av无| 国产三级黄色录像| 亚洲精品成人av观看孕妇| 国产成人免费观看mmmm| 精品人妻熟女毛片av久久网站| 曰老女人黄片| 少妇人妻久久综合中文| 欧美人与善性xxx| 国产精品久久久久久精品电影小说| 国产精品秋霞免费鲁丝片| 熟女av电影| 女性被躁到高潮视频| 国产精品秋霞免费鲁丝片| 欧美黄色片欧美黄色片| 男人爽女人下面视频在线观看| 精品国产一区二区久久| 亚洲男人天堂网一区| 曰老女人黄片| 亚洲国产av影院在线观看| 美女午夜性视频免费| 激情五月婷婷亚洲| 国产一区二区三区综合在线观看| 亚洲av片天天在线观看| 老司机影院毛片| 777久久人妻少妇嫩草av网站| 欧美日韩亚洲综合一区二区三区_| 老司机影院成人| 这个男人来自地球电影免费观看| 悠悠久久av| 精品久久久久久电影网| 美女午夜性视频免费| 搡老乐熟女国产| 五月开心婷婷网| 久久精品亚洲熟妇少妇任你| 十八禁人妻一区二区| 嫩草影视91久久| 国产亚洲欧美精品永久| 深夜精品福利| 99九九在线精品视频|