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

    基于SNP標(biāo)記的辣木群體遺傳分析

    2023-11-11 15:55:02普天磊韓學(xué)琴羅會(huì)英鄧紅山鄒枚伶金杰夏志強(qiáng)王文泉
    熱帶作物學(xué)報(bào) 2023年9期
    關(guān)鍵詞:辣木遺傳多樣性

    普天磊 韓學(xué)琴 羅會(huì)英 鄧紅山 鄒枚伶 金杰 夏志強(qiáng) 王文泉

    關(guān)鍵詞:辣木;SNP;雜合度;群體結(jié)構(gòu);遺傳多樣性

    中圖分類號(hào):S792.99 文獻(xiàn)標(biāo)識(shí)碼:A

    辣木(Moringa oleifera Lam.)屬于辣木科辣木屬的多年生落葉喬木,又被稱為鼓槌樹(shù),辣木原產(chǎn)于印度,是埃塞俄比亞、尼日利亞、菲律賓和蘇丹的重要農(nóng)作物,在非洲、美洲等熱帶亞熱帶地區(qū)均有分布,在我國(guó)的云南、海南、福建、廣州等地均有種植[1]。辣木有13 個(gè)種,其葉片中約含有20%~30%的蛋白質(zhì),葉片、花、果實(shí)含有豐富的維生素A、維生素B、維生素C 和鈣、鎂等礦物質(zhì),種子含高油酸,可用作化妝品、烹飪和機(jī)械潤(rùn)滑油,種子榨油后的剩余物可用于凈化污水、飼喂動(dòng)物[2-4]。同時(shí)辣木含有豐富的皂苷、生物堿、黃酮、酚類等次生代謝產(chǎn)物,具有抗氧化、抗炎、細(xì)胞保護(hù)、神經(jīng)保護(hù)、抗癌等藥理作用[5-6]。

    SNP 標(biāo)記相比較于RFLP、SSR 等傳統(tǒng)分子標(biāo)記而言,可檢測(cè)單個(gè)堿基的插入、缺失、轉(zhuǎn)換和顛換,具有變異數(shù)量多,分布廣,遺傳穩(wěn)定性高,檢測(cè)快、通量高的優(yōu)點(diǎn)[7]?;赟NP 標(biāo)記進(jìn)行的遺傳分析在植物學(xué)領(lǐng)域應(yīng)用較多,例如,高嵩等[8]利用SNP 芯片進(jìn)行玉米遺傳多樣性、群體遺傳結(jié)構(gòu)和類群間遺傳關(guān)系分析,選育并審定了玉米新品種;韓志剛等[9]基于SNP 標(biāo)記對(duì)148 份馬鈴薯種質(zhì)遺傳多樣性進(jìn)行分析,認(rèn)為馬鈴薯絕大部分栽培種遺傳相似性高,遺傳背景不夠豐富。目前,國(guó)內(nèi)并沒(méi)有利用SNP 分子標(biāo)記對(duì)辣木種群的遺傳學(xué)進(jìn)行分析的報(bào)道。AFSM 技術(shù)為簡(jiǎn)化基因組測(cè)序技術(shù),該法分別利用EcoR I-Msp I 和EcoR I-Hpa II 兩種酶對(duì)基因組DNA 進(jìn)行雙酶切,并在兩端加上區(qū)分不同樣本的標(biāo)簽和接頭,樣品混合后進(jìn)行雙端測(cè)序,測(cè)序后獲得的SNP 標(biāo)記數(shù)量多,比傳統(tǒng)分子標(biāo)記更好地代表全基因組的遺傳信息,具有成本低、準(zhǔn)確性和穩(wěn)定性高、易于操作的優(yōu)點(diǎn)[10]。

    雜合度分析有助于深入了解辣木的遺傳組成情況,確定繁育類型,合理規(guī)劃育種,傳統(tǒng)研究繁育類型的方法主要是基于對(duì)花器官的形態(tài)學(xué)分析,傳粉媒介的觀察以及溫室雜交試驗(yàn)展開(kāi),主要通過(guò)表型性狀進(jìn)行評(píng)估,易受環(huán)境、氣候、栽培措施等因素影響,不能準(zhǔn)確地反映植物基因型[11-13]。辣木群體結(jié)構(gòu)的研究對(duì)于辣木種質(zhì)資源的挖掘、利用和保護(hù)具有重要的理論和實(shí)踐意義,遺傳多樣性及群體分化分析是遺傳學(xué)研究的核心內(nèi)容,親本的遺傳關(guān)系很大程度決定子代種子的質(zhì)量,親本間存在差異的遺傳信息會(huì)隨著雜交或自交過(guò)程傳遞給子代,使之在單核苷酸水平上呈現(xiàn)出來(lái)。目前辣木的繁育類型頗受爭(zhēng)議,還沒(méi)有學(xué)者基于SNP 對(duì)辣木的繁育類型進(jìn)行研究,國(guó)內(nèi)辣木育種工作進(jìn)展緩慢,沒(méi)有自主產(chǎn)權(quán)的辣木品種,存在種子繁殖會(huì)發(fā)生性狀分離及種子管理不規(guī)范等問(wèn)題,造成辣木優(yōu)良品種缺乏、品種混亂的現(xiàn)象[14],辣木親本和子代群體的遺傳分析對(duì)于確定辣木繁育類型、分析親緣關(guān)系及選育優(yōu)良品種具有重要的意義。

    本研究以96 份辣木為研究材料,結(jié)合基因組AFSM 高通量測(cè)序技術(shù),與參考基因組進(jìn)行比對(duì)后,進(jìn)行基因型分析、雜合度分析、群體結(jié)構(gòu)、遺傳多樣性、群體分化及連鎖不平衡分析,以揭示辣木親本與子代間的遺傳關(guān)系,為辣木繁育類型和種質(zhì)親緣關(guān)系提供理論指導(dǎo),以及為發(fā)掘控制辣木種質(zhì)優(yōu)良性狀的優(yōu)異等位基因提供理論依據(jù)。

    1 材料與方法

    1.1 材料

    選取來(lái)源于同一母本通過(guò)自然授粉得到的YMLM002 辣木種子94 粒,該種質(zhì)是經(jīng)過(guò)連續(xù)3a的跟蹤觀測(cè)篩選出的果用型優(yōu)良單株材料,具有產(chǎn)量高、果型好、種子飽滿的特點(diǎn)[15]。辣木種子先用清水浸泡10 h,軟化種子硬殼,再用100 mg/L高錳酸鉀溶液浸泡0.5 h 消毒,清水洗凈后點(diǎn)播于穴盤中(紅土∶蛭石=1∶1),適時(shí)補(bǔ)充水分保證濕潤(rùn),待苗長(zhǎng)至15 cm 左右,收集94 份子代樣品、1 份母株和1 份扦插苗樣品備用。

    1.2 方法

    1.2.1 辣木基因組DNA 提取及建庫(kù) 采用CTAB 法提取辣木樣品DNA,用Nano Drop ND-1000 對(duì)DNA 樣品濃度進(jìn)行檢測(cè),并調(diào)節(jié)樣品濃度至100 ng/μL,置于–20 ℃保存。采用AFSM 技術(shù)進(jìn)行建庫(kù),利用EcoR I-Msp I 和EcoR I-Hpa II兩種酶對(duì)96 份辣木DNA 樣品進(jìn)行混合雙酶切,再將酶切產(chǎn)物連接加上用于區(qū)分不同樣品的接頭標(biāo)簽,純化后進(jìn)行PCR 擴(kuò)增,樣品混合后再用高通量測(cè)序平臺(tái)Illumina 進(jìn)行雙端測(cè)序,并計(jì)算GC含量和Q30 評(píng)估測(cè)序數(shù)據(jù)質(zhì)量。

    1.2.2 辣木群體基因型分析 利用Perl 腳本對(duì)原始測(cè)序數(shù)據(jù)進(jìn)行過(guò)濾,使用Bowtie 軟件將過(guò)濾數(shù)據(jù)比對(duì)到辣木參考基因組ASM980114v1,再使用VCFtools 和BCFtools 軟件檢測(cè)并統(tǒng)計(jì)SNP 和Indel 位點(diǎn)信息。

    1.2.3 辣木基因雜合度分析 使用AWK 語(yǔ)言分析96 份樣品的雜合位點(diǎn),并計(jì)算個(gè)體內(nèi)基因的雜合位點(diǎn)比率即為個(gè)體內(nèi)雜合度;同時(shí)通過(guò)將子代數(shù)據(jù)分別與親本進(jìn)行比對(duì),找出差異位點(diǎn),統(tǒng)計(jì)差異位點(diǎn)概率即為子代與親本比對(duì)雜合度,分別生成個(gè)體內(nèi)雜合度及子代與親本比對(duì)雜合度統(tǒng)計(jì)圖。

    1.2.4 群體結(jié)構(gòu)分析、遺傳多樣性及群體分化分析 利用Plink 對(duì)變異位點(diǎn)進(jìn)行過(guò)濾,過(guò)濾掉最小等位基因頻率低于0.05 及基因型缺失率小于5%的位點(diǎn),哈迪–溫伯格檢驗(yàn)顯著性P>0.0001,保留高質(zhì)量的變異位點(diǎn),再使用ADMIXTURE 軟件進(jìn)行群體結(jié)構(gòu)分析,將亞群數(shù)K 值范圍設(shè)置為1~10,根據(jù)得到的交叉驗(yàn)證錯(cuò)誤率(cross-validationerror, CV error)值選擇合適的亞群數(shù)K 值,以個(gè)體占亞群的遺傳成分系數(shù)確定個(gè)體歸屬的類群,用R 軟件繪制群體遺傳結(jié)構(gòu)矩陣圖。采用GCTA 軟件對(duì)過(guò)濾得到的高質(zhì)量文件進(jìn)行主成分分析,并用R 軟件繪圖;采用VCFtools 軟件計(jì)算辣木群體的遺傳多樣性指數(shù)(π)及群體分化指數(shù)(Fst)。

    1.2.5 連鎖不平衡分析 使用LDBlockShow 軟件,計(jì)算不同標(biāo)記距離下的D值,并生成單體型塊圖以展示位點(diǎn)間的連鎖不平衡程度。

    2 結(jié)果與分析

    2.1 辣木群體基因型分析

    96 份辣木樣品基因組DNA 經(jīng)過(guò)AFSM技術(shù)建庫(kù)、測(cè)序,將數(shù)據(jù)過(guò)濾并比對(duì)至辣木參考基因組ASM980114v1,參考基因組信息來(lái)源于NCBI 數(shù)據(jù)庫(kù)(https://www.ncbi.nlm.nih.gov/data-hub/taxonomy/3735/),基因組大小為253.9 Mb,測(cè)序得到1.8 G 數(shù)據(jù)文件,346 615 757 條reads,測(cè)序長(zhǎng)度為150 bp,平均GC 值為50.53%,平均Q30 為94.49%。采用VCFtools 和BCFtools 軟件處理樣品基因組數(shù)據(jù)后,得到1 187 831 個(gè)SNP 和150 861 個(gè)Indel位點(diǎn),以及11 158 個(gè)多等位基因位點(diǎn),4930 個(gè)多等位基因SNP 位點(diǎn)。SNP 同類型堿基之間的突變?yōu)檗D(zhuǎn)換,不同類型堿基之間的突變?yōu)轭崜Q,SNPs發(fā)生轉(zhuǎn)換概率與顛換概率的比值為2.08,單一序列發(fā)生轉(zhuǎn)換次數(shù)為804 031,單一序列發(fā)生顛換次數(shù)為383 471(圖1)。

    SNPs 發(fā)生堿基轉(zhuǎn)換和顛換概率隨著位點(diǎn)的增大呈現(xiàn)先逐漸增加后緩慢降低的趨勢(shì)(圖1)。辣木不同類型的突變位點(diǎn)中,堿基轉(zhuǎn)換的變異數(shù)量顯著大于顛換的數(shù)量,其中堿基G/A 和C/T 的替換都較高,分別為243 672 和241 616 次;堿基A/G和T/C 的替換次之,分別為161 258 和158 648 次,堿基插入/缺失發(fā)生的次數(shù)隨著堿基插入/缺失長(zhǎng)度的增加而呈現(xiàn)出迅速下降的趨勢(shì)。

    2.2 辣木自然結(jié)實(shí)子代基因雜合度分析

    采用1 187 831 個(gè)SNP 位點(diǎn)和150 861 個(gè)Indel位點(diǎn)分析96 份辣木樣品的雜合度(圖2)。辣木同源染色體上的SNP 位點(diǎn)為同一類型堿基,則該SNP 位點(diǎn)稱為純合SNP 位點(diǎn),若為不同類型堿基,則為SNP 雜合位點(diǎn)。由圖2A 可知,辣木個(gè)體內(nèi)雜合度在10.79%~0.36%之間,個(gè)體內(nèi)平均雜合度為4.89%,其中,母株雜合度為5.65%,扦插苗雜合度為5.34%。由圖2B 可知,子代與親本比對(duì)雜合度在21.22%~35.33%之間,子代與親本的比對(duì)平均雜合度為24.85%。由此可知,導(dǎo)致辣木子代雜合的基因中,約有4.89%的基因?yàn)樽陨黼s合基因,19.96%為外來(lái)遺傳物質(zhì)導(dǎo)致雜合的基因,基本表明辣木通過(guò)自花和異花2 種授粉方式繁衍后代。

    2.3 辣木群體結(jié)構(gòu)分析

    采用Plink 對(duì)變異位點(diǎn)進(jìn)行過(guò)濾后,得到141 323 個(gè)SNP 位點(diǎn),再通過(guò)軟件利用所有SNP和Indels 分子標(biāo)記對(duì)96 份辣木樣品進(jìn)行群體遺傳結(jié)構(gòu)分析,由圖3A 可知,當(dāng)K 值為3 時(shí),隨著K 值的增大,CV error 逐漸增大。由于K 值為2和3 時(shí),CV error 值均較小且較為接近,分別為0.401 和0.404,但當(dāng)K 值為2,即將96 份樣品分為2 個(gè)亞群時(shí),各亞群的個(gè)體呈現(xiàn)分布不集中的現(xiàn)象,故將96 份辣木樣品分為3 個(gè)亞群(subgroup1-3)。根據(jù)個(gè)體在3 個(gè)亞群的Q 值,將個(gè)體歸類到Q 值占比最大的亞群(圖3B),發(fā)現(xiàn)3個(gè)亞群中分別有47、31、18 份材料,其中母株和扦插苗屬于亞群1,亞群2、亞群3 均為子代樣品。

    主成分分析發(fā)現(xiàn)(圖3C),亞群1 和亞群2在PC1 軸上有分布差距,而亞群3 與亞群1、2在PC2 軸上有一定的分布差距。大部分亞群可以聚類在一起,表明聚類結(jié)果與群體結(jié)構(gòu)的劃分具有一致性。同時(shí),上述結(jié)果(辣木親本與子代樣本聚類為3 個(gè)亞群)再次論證了雜合度分析結(jié)果,即在生殖遺傳的過(guò)程中,辣木并非以自花授粉的方式繁衍后代,在一定程度上接受了外來(lái)的花粉,導(dǎo)致后代在不包含母株的另外2 個(gè)群體中有分布。

    由圖4 可知,群體進(jìn)化樹(shù)的聚類結(jié)果與群體結(jié)構(gòu)的劃分相一致,各亞群大致能聚在一起,且樣品間有一定的交叉。相比較而言,亞群1 的分枝長(zhǎng)度較短,有4 個(gè)個(gè)體分散在亞群2 中;亞群2 的分布總體集中,有7 個(gè)個(gè)體與亞群3 有交叉;同時(shí),亞群3 有3 個(gè)個(gè)體與亞群1 有交叉。

    2.4 辣木群體遺傳多樣性及分化分析

    3 個(gè)亞群的π 值差距較小且均較低,平均π值也較低,為0.0010,表明96 個(gè)辣木群體的遺傳多樣性水平低。各亞群的Fst 在0.0049~0.0110 之間,其中亞群1 和亞群2 間的Fst 值最小,亞群2和亞群3 間的Fst 值最大,各亞群間的Fst 值均小于0.05,表明各樣本之間存在較弱的遺傳分化(當(dāng)Fst 等于0 或1 時(shí),分別表明亞群間沒(méi)有分化或完全分化;當(dāng)Fst 為0~0.05 時(shí),表明亞群間的分化較弱;當(dāng)Fst 為0.05~0.15 時(shí),表明亞群間為中度分化;當(dāng)Fst 為0.15~0.25 時(shí),表明亞群間的分化較強(qiáng)[16]),各亞群間的親緣關(guān)系相對(duì)較近。

    2.5 連鎖不平衡分析

    結(jié)合多態(tài)性核心SNP 位點(diǎn)在辣木基因組上對(duì)應(yīng)位置分析,發(fā)現(xiàn)共有136 個(gè)Scaffold,主要Scaffold 統(tǒng)計(jì)情況見(jiàn)表1(SNP 數(shù)量前十的Scaffold)。其中,Scaffold 1 的SNP 數(shù)目最多,為62 225個(gè),Scaffold 122 的SNP 數(shù)目最少,為288 個(gè)。通過(guò)LDBlockShow 軟件對(duì)Scaffold 1 在6.748~6.749 Mb 區(qū)域內(nèi)的變異信息進(jìn)行連鎖不平衡分析,發(fā)現(xiàn)6 748 044~6 748 185 位點(diǎn)之間具有強(qiáng)連鎖不平衡關(guān)系,而6 748 040 與6 748 041、6 748 041與6 748 044 等位點(diǎn)間的連鎖關(guān)系弱(圖5)。

    3 討論

    AFSM 技術(shù)采用EcoR I-Hpa II 和EcoR I-MspI 兩組雙酶切體系簡(jiǎn)化基因組DNA 的復(fù)雜度,目前已發(fā)展得較為成熟,已用于檢測(cè)巴西木薯、澳洲堅(jiān)果、麻瘋樹(shù)等植物的SNP、Indel 及甲基化位點(diǎn)[17-19],該技術(shù)DNA 處理步驟和數(shù)據(jù)分析步驟相對(duì)簡(jiǎn)單,效率高,測(cè)定的位點(diǎn)穩(wěn)定,無(wú)需進(jìn)行超聲剪切或熒光標(biāo)記,試驗(yàn)成本低,適用于大量非模型物種的基因分型。本研究利用該技術(shù)得到1 187 831 個(gè)SNP 和150 861 個(gè)Indel 位點(diǎn),可實(shí)現(xiàn)辣木親本及子代遺傳分析的目的。

    國(guó)內(nèi)外相關(guān)學(xué)者從不同的角度對(duì)辣木的繁育系統(tǒng)進(jìn)行研究,呂亞等[20]發(fā)現(xiàn)狹瓣辣木在開(kāi)花第一天就有花粉活力和微弱的柱頭可授性,且開(kāi)花之初柱頭高于雄蕊,之后逐漸低于雄蕊。MULUVI等[21]利用AFLP 分子標(biāo)記研究肯尼亞種源辣木的繁育系統(tǒng),表明該種源辣木種子是自交和異交的混合產(chǎn)物。起國(guó)海[22]研究辣木對(duì)干熱河谷傳粉網(wǎng)絡(luò)的影響,并表明辣木單花能提供5~30 μL,含糖量高達(dá)60.5%的花蜜報(bào)酬物,屬于昆蟲(chóng)傳粉植物,主要傳粉者為蜂類。本研究中辣木個(gè)體內(nèi)平均雜合度4.89%,子代與親本的比對(duì)平均雜合度為24.85%,表明辣木繁殖方式為自交與異交同時(shí)存在。因此,在進(jìn)行辣木雜種優(yōu)勢(shì)利用時(shí),需要關(guān)注相關(guān)個(gè)體間的隨機(jī)化分布和最小距離,以最大限度地增加差異品種/系間的雜交受精,并盡量減少品種內(nèi)部的自交。

    植物的繁育系統(tǒng)、選擇、遺傳漂移、突變和遷移是影響植物群體遺傳結(jié)構(gòu)的進(jìn)化因子[23],本研究利用ADMIXTURE 軟件對(duì)辣木的群體結(jié)構(gòu)進(jìn)行分析,將96 個(gè)辣木群體劃分成3 個(gè)亞群,該結(jié)果與聚類分析和主成分分析的結(jié)果類似,3 種群體結(jié)構(gòu)分析方法相互補(bǔ)充印證,表明辣木群體的遺傳結(jié)構(gòu)劃分可靠。在群體結(jié)構(gòu)劃分中,大部分亞群可以聚類在一起,其中,1 亞群有親本及子代樣品,而2、3 亞群均為子代樣品,該結(jié)果表明辣木自然結(jié)實(shí)子代群體除了攜帶親本的遺傳信息外,還攜帶有外來(lái)的遺傳信息,即辣木繁衍后代的方式不僅為自花授粉,而且還存在異花授粉。子代樣品在3 個(gè)亞群中均有分布,可能是由于親本植株種植于保存有不同辣木種質(zhì)的資源圃內(nèi),不同來(lái)源材料的花粉傳播至親本植株所導(dǎo)致的。

    群體的π 值和Fst 值是衡量群體遺傳分化程度的重要參數(shù),RAJALAKSHMI 等[24]使用ISSR、SRAP 標(biāo)記研究印度種源辣木的遺傳多樣性,表明辣木的平均遺傳分化系數(shù)為0.15,總遺傳多樣性指數(shù)為0.17。本研究發(fā)現(xiàn)辣木群體的π 值為0.0010,F(xiàn)st 值在0.0049~0.0110 之間,表明本研究所用辣木群體遺傳分化弱,遺傳多樣性較低,該現(xiàn)象應(yīng)該是由于選取的辣木群體是親本與子代親緣關(guān)系較近造成的。同時(shí),也表明引起子代雜合的外來(lái)基因與親本的基因型差異不大,這可能是由于本研究文采用的是栽培種辣木資源,經(jīng)歷了多次的人工選擇育種,資源間豐度低造成的,后續(xù)可引進(jìn)印度、非洲種源的優(yōu)良辣木種質(zhì),以豐富資源圃內(nèi)的辣木種質(zhì)[25]。

    當(dāng)2 個(gè)距離較近的等位基因在同一單體型上同時(shí)出現(xiàn)的頻率高于隨機(jī)出現(xiàn)的頻率時(shí),表明它們處于連鎖不平衡狀態(tài)。在定位克隆中,通過(guò)連鎖可檢測(cè)到產(chǎn)生連鎖信號(hào)的變異,在關(guān)聯(lián)分析中,利用鄰近位點(diǎn)形成的強(qiáng)連鎖不平衡,有助于找到與性狀相關(guān)的位點(diǎn)[26-27]。本文對(duì)辣木的SNP 位點(diǎn)進(jìn)行了連鎖不平衡分析,并在單體型塊圖上發(fā)現(xiàn)了連鎖不平衡關(guān)系強(qiáng)的基因區(qū)域,可為研究多個(gè)處于連鎖不平衡的位點(diǎn)與重要性狀的關(guān)聯(lián)性提供參考依據(jù)。

    猜你喜歡
    辣木遺傳多樣性
    不同品種辣木全氮的動(dòng)態(tài)變化特征研究
    緬甸辣木受國(guó)外市場(chǎng)追捧
    眾耘牌辣木健康食品 人人喜愛(ài)成新寵
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    辣木之歌
    辣木的抗腫瘤活性及作用機(jī)制研究進(jìn)展
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    楊梅種質(zhì)資源遺傳多樣性研究進(jìn)展
    成年版毛片免费区| 麻豆久久精品国产亚洲av| 99国产极品粉嫩在线观看| 国产精品免费一区二区三区在线| 深夜精品福利| 午夜福利在线观看免费完整高清在 | 色哟哟哟哟哟哟| 天堂网av新在线| 国产亚洲av高清不卡| 成年版毛片免费区| 在线观看免费视频日本深夜| 午夜福利在线观看免费完整高清在 | 又粗又爽又猛毛片免费看| 精品不卡国产一区二区三区| 女同久久另类99精品国产91| 老司机在亚洲福利影院| 国产精品98久久久久久宅男小说| 小说图片视频综合网站| 不卡av一区二区三区| 国产亚洲精品综合一区在线观看| 中文字幕高清在线视频| 亚洲欧洲精品一区二区精品久久久| www.熟女人妻精品国产| 亚洲天堂国产精品一区在线| 精品一区二区三区av网在线观看| 99国产精品99久久久久| 国产一级毛片七仙女欲春2| 久久久成人免费电影| 亚洲七黄色美女视频| 国产日本99.免费观看| 成人一区二区视频在线观看| 精品电影一区二区在线| 亚洲国产欧美人成| 中文字幕精品亚洲无线码一区| 日本黄色片子视频| 亚洲va日本ⅴa欧美va伊人久久| 又粗又爽又猛毛片免费看| 老汉色av国产亚洲站长工具| 丝袜人妻中文字幕| 最近最新中文字幕大全电影3| 色老头精品视频在线观看| 成人三级黄色视频| 亚洲七黄色美女视频| 国产精品亚洲av一区麻豆| 亚洲av成人一区二区三| 亚洲国产精品sss在线观看| 久久久精品欧美日韩精品| 操出白浆在线播放| 亚洲国产欧美网| 97碰自拍视频| 香蕉丝袜av| 亚洲人成伊人成综合网2020| 亚洲国产日韩欧美精品在线观看 | 草草在线视频免费看| 国产精品av久久久久免费| 欧美日韩福利视频一区二区| 黑人操中国人逼视频| 18禁观看日本| 99精品欧美一区二区三区四区| 色精品久久人妻99蜜桃| 日本一本二区三区精品| 99久久久亚洲精品蜜臀av| 97超视频在线观看视频| 麻豆成人午夜福利视频| 一区二区三区国产精品乱码| 在线视频色国产色| 91久久精品国产一区二区成人 | 日韩欧美在线二视频| 岛国在线观看网站| www日本黄色视频网| 国产黄色小视频在线观看| 夜夜夜夜夜久久久久| 人人妻人人看人人澡| 国内少妇人妻偷人精品xxx网站 | 一个人看的www免费观看视频| 五月伊人婷婷丁香| 99riav亚洲国产免费| 99久久成人亚洲精品观看| 黄色日韩在线| 欧美黄色片欧美黄色片| 狂野欧美白嫩少妇大欣赏| 国产成人精品久久二区二区91| 黄色片一级片一级黄色片| 啦啦啦观看免费观看视频高清| 不卡一级毛片| 国产一区二区三区在线臀色熟女| 亚洲精品456在线播放app | 51午夜福利影视在线观看| 成人鲁丝片一二三区免费| 中文资源天堂在线| 欧美在线黄色| 长腿黑丝高跟| 一级作爱视频免费观看| 香蕉久久夜色| 中文字幕人成人乱码亚洲影| 狠狠狠狠99中文字幕| 午夜福利在线观看吧| 性欧美人与动物交配| 高清在线国产一区| 精品一区二区三区视频在线观看免费| 国产精品日韩av在线免费观看| 桃红色精品国产亚洲av| 国产精品亚洲av一区麻豆| 国产视频内射| 五月玫瑰六月丁香| 亚洲熟妇熟女久久| 极品教师在线免费播放| 亚洲五月天丁香| 禁无遮挡网站| 国产精品亚洲av一区麻豆| 桃色一区二区三区在线观看| 亚洲avbb在线观看| av黄色大香蕉| 国产综合懂色| 好看av亚洲va欧美ⅴa在| 黄色日韩在线| 男人舔女人的私密视频| 丁香欧美五月| 国产精品,欧美在线| 天堂√8在线中文| 日本a在线网址| 亚洲国产色片| 怎么达到女性高潮| 人人妻,人人澡人人爽秒播| 国产精品女同一区二区软件 | 在线a可以看的网站| 国产精品久久久人人做人人爽| 国产成人影院久久av| 日本一本二区三区精品| 亚洲真实伦在线观看| 一夜夜www| 熟女少妇亚洲综合色aaa.| 精品不卡国产一区二区三区| 视频区欧美日本亚洲| 日本黄色片子视频| 女同久久另类99精品国产91| a级毛片a级免费在线| 久久久国产成人精品二区| 日韩精品中文字幕看吧| 变态另类丝袜制服| 99久久精品国产亚洲精品| 亚洲欧美日韩东京热| 国内揄拍国产精品人妻在线| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 久久久成人免费电影| 欧美成狂野欧美在线观看| 丁香六月欧美| 五月伊人婷婷丁香| 久久久国产成人免费| 亚洲中文日韩欧美视频| 久久这里只有精品中国| 国产真实乱freesex| 91在线观看av| 一区福利在线观看| 午夜精品久久久久久毛片777| 亚洲国产欧美网| 亚洲成人中文字幕在线播放| 真人一进一出gif抽搐免费| 看黄色毛片网站| 巨乳人妻的诱惑在线观看| a级毛片在线看网站| 亚洲精品乱码久久久v下载方式 | 久久久水蜜桃国产精品网| 97超级碰碰碰精品色视频在线观看| 亚洲最大成人中文| 国产av不卡久久| e午夜精品久久久久久久| 日韩精品中文字幕看吧| 亚洲熟妇熟女久久| 亚洲国产中文字幕在线视频| ponron亚洲| av在线天堂中文字幕| 色综合欧美亚洲国产小说| 在线观看免费视频日本深夜| 叶爱在线成人免费视频播放| 国产精品一区二区免费欧美| 美女免费视频网站| 国产真实乱freesex| 国产成+人综合+亚洲专区| 九九久久精品国产亚洲av麻豆 | 老熟妇乱子伦视频在线观看| 久久久久久久午夜电影| 级片在线观看| 亚洲欧美日韩无卡精品| 国产真实乱freesex| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区三| 亚洲成人久久爱视频| 一个人看视频在线观看www免费 | 欧美乱妇无乱码| 亚洲国产欧美网| 久久精品夜夜夜夜夜久久蜜豆| 香蕉国产在线看| 亚洲av片天天在线观看| 少妇丰满av| 亚洲成av人片免费观看| 一本综合久久免费| 十八禁网站免费在线| 免费在线观看成人毛片| 午夜两性在线视频| 亚洲av第一区精品v没综合| 国产人伦9x9x在线观看| 在线视频色国产色| 国产真实乱freesex| 超碰成人久久| 琪琪午夜伦伦电影理论片6080| 成人三级做爰电影| 三级国产精品欧美在线观看 | 欧美激情久久久久久爽电影| 美女黄网站色视频| 99在线人妻在线中文字幕| 国产成人aa在线观看| 国产精品国产高清国产av| 美女高潮喷水抽搐中文字幕| 亚洲成a人片在线一区二区| 亚洲欧美精品综合一区二区三区| 国产黄a三级三级三级人| 首页视频小说图片口味搜索| 免费在线观看影片大全网站| 久久亚洲真实| 午夜日韩欧美国产| 国产精品久久久人人做人人爽| 黄色片一级片一级黄色片| 欧美日韩瑟瑟在线播放| 欧美性猛交黑人性爽| 亚洲精品乱码久久久v下载方式 | 听说在线观看完整版免费高清| 亚洲国产中文字幕在线视频| 此物有八面人人有两片| 给我免费播放毛片高清在线观看| 老熟妇乱子伦视频在线观看| 最近视频中文字幕2019在线8| av黄色大香蕉| 国产真人三级小视频在线观看| 后天国语完整版免费观看| 欧美绝顶高潮抽搐喷水| 国产一区二区激情短视频| 欧美黄色淫秽网站| 在线观看舔阴道视频| 人人妻人人澡欧美一区二区| 久久中文字幕一级| 两性夫妻黄色片| 色尼玛亚洲综合影院| 天天躁日日操中文字幕| 在线观看日韩欧美| 国产单亲对白刺激| 99国产精品99久久久久| 精品福利观看| 少妇人妻一区二区三区视频| 1024香蕉在线观看| 久久久久九九精品影院| 成人鲁丝片一二三区免费| 黑人操中国人逼视频| 国内精品美女久久久久久| 久久99热这里只有精品18| 在线观看美女被高潮喷水网站 | 久久久久免费精品人妻一区二区| ponron亚洲| 1024香蕉在线观看| 免费搜索国产男女视频| 两个人看的免费小视频| 男女做爰动态图高潮gif福利片| 午夜激情欧美在线| 动漫黄色视频在线观看| 国产精品 国内视频| 在线国产一区二区在线| 久久久水蜜桃国产精品网| 国产精品 欧美亚洲| 久久欧美精品欧美久久欧美| av女优亚洲男人天堂 | 欧美色视频一区免费| 免费无遮挡裸体视频| 90打野战视频偷拍视频| 亚洲 欧美一区二区三区| 男女做爰动态图高潮gif福利片| 在线免费观看不下载黄p国产 | 亚洲成av人片免费观看| 婷婷亚洲欧美| 51午夜福利影视在线观看| 国产精品野战在线观看| 一二三四社区在线视频社区8| 国产成人啪精品午夜网站| av中文乱码字幕在线| 两性夫妻黄色片| 老司机在亚洲福利影院| 国产高清有码在线观看视频| 久久久久久久久中文| 亚洲av免费在线观看| 校园春色视频在线观看| 精华霜和精华液先用哪个| 国产成人av教育| 精品国产超薄肉色丝袜足j| cao死你这个sao货| 亚洲一区二区三区不卡视频| 舔av片在线| 夜夜躁狠狠躁天天躁| 一区福利在线观看| 亚洲欧美精品综合久久99| 亚洲av成人不卡在线观看播放网| 欧美日韩福利视频一区二区| 日韩中文字幕欧美一区二区| 国产精品 国内视频| 国产免费男女视频| 亚洲性夜色夜夜综合| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三| 中出人妻视频一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成a人片在线一区二区| 每晚都被弄得嗷嗷叫到高潮| 香蕉国产在线看| 在线a可以看的网站| 两个人视频免费观看高清| 91在线观看av| 一个人免费在线观看的高清视频| 一级毛片高清免费大全| 美女黄网站色视频| 久久午夜综合久久蜜桃| 国产免费av片在线观看野外av| 国产成人欧美在线观看| 亚洲av成人不卡在线观看播放网| bbb黄色大片| 午夜福利在线观看免费完整高清在 | 免费av不卡在线播放| 18禁美女被吸乳视频| 亚洲av成人一区二区三| 99在线人妻在线中文字幕| 国产综合懂色| 一a级毛片在线观看| 国产成人欧美在线观看| 在线视频色国产色| 成年女人毛片免费观看观看9| 神马国产精品三级电影在线观看| 熟女人妻精品中文字幕| 日韩人妻高清精品专区| 国产麻豆成人av免费视频| 老汉色∧v一级毛片| 久久久久九九精品影院| 小说图片视频综合网站| 这个男人来自地球电影免费观看| 国产精品一区二区三区四区免费观看 | a级毛片a级免费在线| 免费在线观看日本一区| 久久精品影院6| 黑人欧美特级aaaaaa片| 国产私拍福利视频在线观看| 日韩欧美国产在线观看| 色吧在线观看| 国产精品一区二区免费欧美| 天堂动漫精品| 国产探花在线观看一区二区| 亚洲精品国产精品久久久不卡| 淫妇啪啪啪对白视频| 久久天堂一区二区三区四区| 国产久久久一区二区三区| 日韩 欧美 亚洲 中文字幕| 欧美3d第一页| 丁香欧美五月| 中出人妻视频一区二区| 国产一区在线观看成人免费| 欧美中文日本在线观看视频| 91av网一区二区| 国产精品亚洲av一区麻豆| 神马国产精品三级电影在线观看| 亚洲成a人片在线一区二区| 久久久久久人人人人人| 国产成人一区二区三区免费视频网站| 日韩免费av在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕一区二区三区有码在线看 | 成人精品一区二区免费| 美女高潮喷水抽搐中文字幕| 两个人看的免费小视频| 十八禁网站免费在线| 亚洲在线观看片| 男女床上黄色一级片免费看| 成人一区二区视频在线观看| 九九热线精品视视频播放| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 久久精品国产综合久久久| 国产91精品成人一区二区三区| 无人区码免费观看不卡| 欧美+亚洲+日韩+国产| 日韩有码中文字幕| 国产极品精品免费视频能看的| 桃红色精品国产亚洲av| 黄色日韩在线| 叶爱在线成人免费视频播放| 亚洲中文日韩欧美视频| 听说在线观看完整版免费高清| 伦理电影免费视频| 国产成人系列免费观看| 欧美另类亚洲清纯唯美| 九九久久精品国产亚洲av麻豆 | 亚洲电影在线观看av| 国产精品爽爽va在线观看网站| 国产97色在线日韩免费| 精品一区二区三区视频在线 | 婷婷丁香在线五月| 国产极品精品免费视频能看的| or卡值多少钱| 国产精品久久电影中文字幕| 欧美三级亚洲精品| 免费观看人在逋| 日本一本二区三区精品| 色在线成人网| 十八禁人妻一区二区| 国内精品一区二区在线观看| 99国产综合亚洲精品| 亚洲欧美日韩高清在线视频| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 搞女人的毛片| 国产欧美日韩精品亚洲av| 精品午夜福利视频在线观看一区| 黄色日韩在线| 国产精品自产拍在线观看55亚洲| 国产精品 欧美亚洲| 日本五十路高清| 一边摸一边抽搐一进一小说| 亚洲无线观看免费| 亚洲精品美女久久av网站| 99久久综合精品五月天人人| 日韩欧美国产一区二区入口| 青草久久国产| 午夜福利高清视频| 在线看三级毛片| 亚洲成人精品中文字幕电影| www.精华液| 国产欧美日韩一区二区三| 亚洲人成网站高清观看| 成年女人永久免费观看视频| 久久精品亚洲精品国产色婷小说| 亚洲黑人精品在线| 精品国产亚洲在线| 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 午夜久久久久精精品| 十八禁网站免费在线| 日日摸夜夜添夜夜添小说| 亚洲狠狠婷婷综合久久图片| 757午夜福利合集在线观看| 亚洲精华国产精华精| 亚洲 欧美 日韩 在线 免费| 欧美日韩福利视频一区二区| 人人妻人人澡欧美一区二区| 老司机福利观看| 一级毛片高清免费大全| 在线免费观看的www视频| 看免费av毛片| 天天添夜夜摸| 成人国产综合亚洲| 日本一二三区视频观看| 999久久久国产精品视频| 日本撒尿小便嘘嘘汇集6| 巨乳人妻的诱惑在线观看| 日本在线视频免费播放| 欧美黑人欧美精品刺激| 人人妻人人看人人澡| 精品免费久久久久久久清纯| 女同久久另类99精品国产91| 精品一区二区三区视频在线 | 成年女人永久免费观看视频| 人人妻人人看人人澡| 黄频高清免费视频| 精品不卡国产一区二区三区| 亚洲欧美日韩高清专用| 国产男靠女视频免费网站| 特大巨黑吊av在线直播| 黄片小视频在线播放| 国产激情偷乱视频一区二区| 国产黄色小视频在线观看| 免费av不卡在线播放| 给我免费播放毛片高清在线观看| 久久久久亚洲av毛片大全| 精品午夜福利视频在线观看一区| 亚洲,欧美精品.| 99久久精品一区二区三区| 制服丝袜大香蕉在线| 亚洲欧美精品综合一区二区三区| 一级作爱视频免费观看| 午夜福利18| 丰满人妻熟妇乱又伦精品不卡| 中文字幕人妻丝袜一区二区| 国产熟女xx| 亚洲黑人精品在线| 1000部很黄的大片| 久久久国产成人精品二区| 日韩欧美一区二区三区在线观看| 怎么达到女性高潮| 俄罗斯特黄特色一大片| av福利片在线观看| 99在线视频只有这里精品首页| 日本五十路高清| www.精华液| 十八禁网站免费在线| 在线观看舔阴道视频| 久久亚洲真实| 国产精品综合久久久久久久免费| 亚洲精品456在线播放app | 中文字幕久久专区| 观看美女的网站| 黑人巨大精品欧美一区二区mp4| 久久久国产成人免费| 村上凉子中文字幕在线| 精品久久久久久久久久久久久| 国产一区二区在线av高清观看| 两性午夜刺激爽爽歪歪视频在线观看| 俺也久久电影网| 男人舔奶头视频| 精品久久久久久成人av| 亚洲熟妇熟女久久| 亚洲欧美日韩东京热| 免费观看的影片在线观看| av天堂在线播放| svipshipincom国产片| 哪里可以看免费的av片| 亚洲最大成人中文| 男人和女人高潮做爰伦理| 色在线成人网| 最近最新中文字幕大全电影3| 淫秽高清视频在线观看| 免费无遮挡裸体视频| 真人一进一出gif抽搐免费| 久久久国产欧美日韩av| 国产av一区在线观看免费| 欧美一区二区国产精品久久精品| 色老头精品视频在线观看| 后天国语完整版免费观看| 欧美日本亚洲视频在线播放| 国产综合懂色| 国产精品自产拍在线观看55亚洲| 久久国产精品影院| 老司机午夜十八禁免费视频| 久久久成人免费电影| 成人精品一区二区免费| 男女视频在线观看网站免费| 男人舔女人的私密视频| 国产精品久久电影中文字幕| 国产三级黄色录像| 老司机在亚洲福利影院| 亚洲成人久久爱视频| 一a级毛片在线观看| 啦啦啦韩国在线观看视频| 欧美成人性av电影在线观看| 12—13女人毛片做爰片一| 91麻豆av在线| 床上黄色一级片| 日韩精品中文字幕看吧| 日韩有码中文字幕| avwww免费| 麻豆一二三区av精品| 欧美黑人巨大hd| 亚洲片人在线观看| 真人做人爱边吃奶动态| 久久天躁狠狠躁夜夜2o2o| 波多野结衣高清无吗| 丰满人妻熟妇乱又伦精品不卡| 90打野战视频偷拍视频| 在线a可以看的网站| 99在线视频只有这里精品首页| 五月伊人婷婷丁香| 美女被艹到高潮喷水动态| 亚洲国产精品成人综合色| 日韩人妻高清精品专区| 欧美国产日韩亚洲一区| 最近视频中文字幕2019在线8| 女同久久另类99精品国产91| 美女高潮的动态| 丰满人妻一区二区三区视频av | 国产精品1区2区在线观看.| 亚洲人成电影免费在线| 婷婷六月久久综合丁香| 欧美高清成人免费视频www| 一边摸一边抽搐一进一小说| 亚洲av成人不卡在线观看播放网| 法律面前人人平等表现在哪些方面| 国产一区二区激情短视频| 亚洲欧美精品综合一区二区三区| 99久久国产精品久久久| 免费av毛片视频| 久久精品亚洲精品国产色婷小说| tocl精华| 国内揄拍国产精品人妻在线| 看黄色毛片网站| 十八禁人妻一区二区| 欧美av亚洲av综合av国产av| 99热精品在线国产| 免费搜索国产男女视频| 97碰自拍视频| 久久国产精品人妻蜜桃| 亚洲七黄色美女视频| 99久久久亚洲精品蜜臀av| 法律面前人人平等表现在哪些方面| or卡值多少钱| 成年女人看的毛片在线观看| 色av中文字幕| 亚洲中文av在线| 熟女电影av网| 91麻豆av在线| 九色成人免费人妻av| 中文资源天堂在线| 成人午夜高清在线视频| 国产视频一区二区在线看| 99在线视频只有这里精品首页| 不卡一级毛片| 狂野欧美白嫩少妇大欣赏| 18禁观看日本| 国产97色在线日韩免费| 国产精品久久久av美女十八| 亚洲男人的天堂狠狠|