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

    田螺科4種代表性貝類基因組調(diào)研及SSR特征分析

    2023-05-08 06:23:22季秋婷王俊杰唐伯平
    江蘇農(nóng)業(yè)科學(xué) 2023年6期
    關(guān)鍵詞:田螺貝類雜合

    朱 穎, 陳 瑜, 唐 俊, 季秋婷, 王俊杰, 唐伯平, 王 剛

    (鹽城師范學(xué)院濕地學(xué)院/江蘇省鹽土生物資源研究重點(diǎn)實(shí)驗(yàn)室,江蘇鹽城 224007)

    田螺科(Viviparidae)貝類是我國淡水水域常見的軟體動(dòng)物,其屬于腹足綱、中腹足目,廣泛分布于湖泊、沼澤、水庫、池塘及溪流等處[1],能對(duì)水體健康起到維持作用,凈化水體,是重要的生態(tài)指示物種[2]。田螺科物種多為經(jīng)濟(jì)貝類,其營養(yǎng)價(jià)值高,螺肉鮮美,是低脂肪、高蛋白食品,在我國廣受歡迎,同時(shí)也是家禽,水產(chǎn)養(yǎng)殖魚類、蝦蟹類的天然飼料[3]。田螺科貝類屬于世界物種,我國包含9個(gè)屬,70多個(gè)種,其中河螺屬(Rivularia)和螺螄屬(Margarya)為特有屬,環(huán)棱螺屬和圓田螺屬物種最為常見[4]。隨著資源環(huán)境的破壞和對(duì)田螺的需求日益增加,加上外來淡水螺的不斷入侵,我國本土田螺科物種分布和數(shù)量受到嚴(yán)重威脅,很多種群在數(shù)量上明顯衰退,嚴(yán)重的已瀕臨滅絕[5]?;谔锫菘莆锓N的重要生態(tài)價(jià)值和經(jīng)濟(jì)價(jià)值,保護(hù)和恢復(fù)田螺科物種的多樣性及種群數(shù)量的問題顯得尤為重要。

    隨著組學(xué)時(shí)代的到來和發(fā)展,測序成本下降,高通量測序被廣泛用于軟體動(dòng)物基因組測序中。在軟體動(dòng)物中,首個(gè)貝類——太平洋牡蠣(Pacific oyster)基因組圖譜于2012年完成,標(biāo)志著軟體動(dòng)物組學(xué)時(shí)代的到來,填補(bǔ)了軟體動(dòng)物家族基因組圖譜的空白[6]。隨后,霸王蓮花青螺(Lottiagigantea)[7]、章魚(Octopusbimaculoides)[8]、牡蠣(Crassostreavirginica)[9]等其他軟體動(dòng)物的基因組圖譜陸續(xù)發(fā)表,這為其他軟體動(dòng)物的基因組圖譜的構(gòu)建提供了基礎(chǔ)。截至2022年6月,在NCBI數(shù)據(jù)庫中已公布92個(gè)軟體動(dòng)物基因組圖譜。基因組的解析可為貝類優(yōu)良性狀基因挖掘、遺傳機(jī)制解析提供研究基礎(chǔ)。然而,目前淡水的田螺科貝類分子生物學(xué)研究基礎(chǔ)相對(duì)薄弱,關(guān)于田螺科的研究大多集中在初步調(diào)查、系統(tǒng)發(fā)育和分類地位等方面的分析[10-12]。在NCBI數(shù)據(jù)庫中,提供的田螺科分子數(shù)據(jù)十分有限,缺乏田螺科貝類參考基因組信息。在構(gòu)建基因組精細(xì)譜圖之前進(jìn)行基因組特征調(diào)研,可為后續(xù)的測序方案提供參考依據(jù)。目前,主要采用K-mer分析方法對(duì)基因組進(jìn)行調(diào)研評(píng)估獲取物種雜合率、重復(fù)度等信息[13-15]。本研究選取淡水田螺科的中華園田螺(Cipangopaludinacathayensis)、銅銹環(huán)棱螺(Bellamyaaeruginosa)、多棱角螺(Angulyagrapolyzonata)、湄公螺(Mekongiarivularia)等4種代表性貝類進(jìn)行基因組調(diào)研,通過K-mer分析,進(jìn)行物種雜合率、基因組重復(fù)序列比例及基因組大小的評(píng)估,以期為田螺科貝類基因組學(xué)研究提供參考。同時(shí),利用Illumina 測序數(shù)據(jù)對(duì)4個(gè)田螺科貝類進(jìn)行初步組裝和初步比較分析,采用MISA工具對(duì)簡單重復(fù)序列(SSR)進(jìn)行鑒定和分析,為進(jìn)一步分子標(biāo)記開發(fā)和資源鑒定提供研究基礎(chǔ)。

    1 材料與方法

    1.1 材料

    中華園田螺、銅銹環(huán)棱螺、多棱角螺和湄公螺,于2021年5—6月分別采集自四川(30°23′36.82″N,104°4′50.46″E)、江蘇鹽城(33°22′33.4″N,120°12′7.19″E)、湖南(27°49′36.25″N,113°7′59.75″E)、江西(28°30′7.34″N,115°48′22.93″E),采集的田螺科樣本保存于江蘇鹽城師范學(xué)院沿海灘涂重點(diǎn)實(shí)驗(yàn)室。

    1.2 樣品DNA提取

    在江蘇省鹽城師范學(xué)院沿海灘涂重點(diǎn)實(shí)驗(yàn)室,使用Omega公司的試劑盒D3373-01 Mollusc DNA Kit(50)提取田螺科4種代表性螺的基因組DNA(取螺的部分樣品),用NanoDrop 光譜儀測定DNA 的濃度和純度。DNA樣本標(biāo)準(zhǔn)為:DNA總量≥20 μg 且濃度≥300 ng/μL;樣本純度要求為:D260 nm/280 nm≥1.8,D260 nm/230 nm≥1.8。

    1.3 4種貝類樣本測序和基因組評(píng)估

    DNA樣本送到武漢菲沙基因信息有限公司進(jìn)行文庫制備,基于高通量測序儀(MGI-SEQ2000)平臺(tái)采取雙端150 bp(PE150)技術(shù)進(jìn)行測序,測得的數(shù)據(jù)經(jīng)過質(zhì)控后用來進(jìn)行下一步數(shù)據(jù)分析。4種貝類Illumina測序數(shù)據(jù)存儲(chǔ)在國家基因組科學(xué)數(shù)據(jù)中心數(shù)據(jù)庫(https://ngdc.cncb.ac.cn/),項(xiàng)目編號(hào):PRJCA013294。

    文庫構(gòu)建參照Illumina標(biāo)準(zhǔn)規(guī)程[16],構(gòu)建并插入片段大小為270 bp的文庫[17]并測序。完成構(gòu)建的文庫在Illumina HiSeq 2500平臺(tái)進(jìn)行測序。選用AdapterRemoval(version 2.1.7)軟件對(duì)測序數(shù)據(jù)去接頭,用FastQC軟件進(jìn)行質(zhì)量評(píng)估,并過濾質(zhì)粒污染。采用基于K-mer分析進(jìn)行估算,K值取21,在測序reads均勻分布的前提下,計(jì)算基因組大小(基因組大小=總堿基數(shù)/平均測序深度=總K-mer數(shù)平均值/平均K-mer深度)。利用Jellyfish和 GenomeScope 軟件對(duì)K-mer頻數(shù)分布數(shù)據(jù)進(jìn)行統(tǒng)計(jì)并擬合作圖,得到K-mer分布圖。

    1.4 基因組組裝和GC深度分析

    對(duì)過濾數(shù)據(jù)利用SOAPdenovo軟件進(jìn)行基因組組裝,將讀長(reads)分成100 bp長度的K-mer,再基于 K-mer 數(shù)據(jù)構(gòu)建de Bruijn圖。利用perl腳本對(duì)組裝的基因組進(jìn)行滑窗統(tǒng)計(jì),設(shè)置窗口為10 kb,根據(jù)GC分布和覆蓋深度統(tǒng)計(jì)結(jié)果,應(yīng)用R腳本繪制散點(diǎn)圖。

    1.5 簡單重復(fù)序列(SSR)特征分析

    對(duì)初步完成組裝的基因組進(jìn)行簡單重復(fù)序列分析。使用微衛(wèi)星識(shí)別工具(microsatellite identification tool,MISA)搜索SSR位點(diǎn)對(duì)于150 bp以上的序列片段(Scaffolds)[18]。同時(shí)運(yùn)行MISA軟件程序的配置文件misa.ini,使用MISA默認(rèn)參數(shù)。

    2 結(jié)果與分析

    2.1 田螺科4種貝類基因組測序質(zhì)量和基因組大小評(píng)估

    對(duì)Illumina Hiseq PE測序平臺(tái)測得的中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺等4個(gè)代表性貝類的原始數(shù)據(jù),經(jīng)過濾和校正等質(zhì)控步驟后,分別獲得52.59、58.64、55.71、54.90 G的有效數(shù)據(jù),有效數(shù)據(jù)均大于96%(表1)。測序質(zhì)量評(píng)估結(jié)果顯示,4種螺的Q20分別為 96.5%、96.1%、96.4%、96.6%,Q30分別為88.5%、87.5%、88.3%、88.8%,Q20均大于95%,Q30均大于85%,表明4種貝類的測序質(zhì)量較好,滿足分析要求。4種田螺科貝類GC 含量分別為33.4%、33.7%、33.1%、33.8%,GC含量基本持平。

    表1 田螺科4種螺基因組調(diào)查測序結(jié)果統(tǒng)計(jì)

    本研究采用K-mer 的分析方法,取K為21來進(jìn)行物種基因組特征分析,K-mer分布圖如圖1,實(shí)際觀測的K-mer分布為藍(lán)色區(qū)域;測序頻數(shù)較低的K-mer為紅色線條下方區(qū)域(測序錯(cuò)誤);黑色線條下方是可靠的K-mer數(shù)據(jù),用于4種貝類基因組大小的評(píng)估;黃色線條下方區(qū)域是非重復(fù)區(qū)域。由圖1可知,中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺的主峰分別在測序深度18×、18×、45×、23×左右,均只有1個(gè)主峰,說明4種螺均為2倍體。在測序數(shù)據(jù)中去除深度異常的K-mer后,評(píng)估得到4種螺的基因組大小分別為1 353、896、991、1 196 Mb;雜合度分別為0.77%、1.37%、2.41%、0.91%。4種貝類雜合率均大于0.5%,為高雜合基因組。4種螺的重復(fù)序列比例分別為26.2%、26.7%、27.8%、28.8%,比例基本一致(表2)。

    表2 田螺科4種螺基因組評(píng)估結(jié)果

    2.2 田螺科4種貝類基因組序列初步組裝

    利用SOAPdenovo軟件分別對(duì)4種貝類經(jīng)過質(zhì)控處理后得到的高質(zhì)量read,進(jìn)行初步組裝,中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺的基因組大小分別為 1 598、1 845、1 202、1 322 Mb,其中銅銹環(huán)棱螺基因組長度最大,多棱角螺基因組長度最小。4種貝類的Contigs N50分別為690、200、1 025、889 bp(表3)。初步組裝得到的基因組大小與K-mer評(píng)估得到基因組大小有差異,其組裝的基因組大小偏大,這可能雜合度有關(guān),也有可能與二代測序數(shù)據(jù)讀長較短和覆蓋度較低等因素有關(guān)。

    2.3 田螺科4種貝類初步組裝的GC含量與覆蓋度分布

    對(duì)4種田螺科貝類初步組裝得到的基因組,過濾覆蓋深度10×以下的contig序列,基于覆蓋深度和GC含量繪制成散點(diǎn)圖。由圖2可知,中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺4種螺的GC含量較為一致,在27%~43%范圍內(nèi),但區(qū)域測序深度分布有所差異。其中中華園田螺測序深度較高,在20~40×區(qū)域深度分布比較集中,其次分別是湄公螺(18~25×)、多棱角螺(10~38×)、銅銹環(huán)棱螺(10~20×)。4種田螺科貝類均為1個(gè)主峰,說明無外源污染。

    2.4 田螺科4種貝類SSR特征分析

    對(duì)中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺等4種貝類初步組裝的基因組,通過微衛(wèi)星識(shí)別工具M(jìn)ISA軟件對(duì)組裝基因組的Scaffolds進(jìn)行SSR分析。4種貝類總共分別搜索到 928 562、824 342、739 299、739 407個(gè)SSR。其中,4種貝類分別有 196 158、127 288、158 667、160 694條序列包含1個(gè)以上堿基重復(fù);分別有13 360、127 897、118 293、118 258 個(gè)復(fù)合形式存在。4種貝類SSR位點(diǎn)的核苷酸重復(fù)類型可劃分6種類型。其中,單核苷酸重復(fù)的SSR基元類型最豐富,五核苷酸重復(fù)、六核苷酸重復(fù)基元類型和分布相對(duì)較少。具體而言,中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺的單核苷酸重復(fù)、雙核核苷酸重復(fù)序列占比最高,分別達(dá)79.2%、79.9%、78.7%、80.9%(表4)。對(duì)檢測出的重復(fù)基元進(jìn)行排序,其中,4種貝類優(yōu)勢重復(fù)基元分別是A/T、AC/GT、AAT/ATT、AGAT/ATCT、ATAAT/ATTAT、CTAACC/GGTTAA。對(duì)雙核苷酸重復(fù)、三核苷酸重復(fù)、四核苷酸重復(fù)中基元進(jìn)行比較發(fā)現(xiàn),4種貝類的優(yōu)勢重復(fù)基元大體相同(圖3)。

    表4 田螺科4種螺基因組SSR類型

    3 結(jié)論與討論

    近年來,隨著測序技術(shù)的發(fā)展,越來越多的軟體動(dòng)物基因組被解析,為研究軟體動(dòng)物起源、進(jìn)化、生殖、發(fā)育、性別調(diào)控和免疫等問題提供分子生物學(xué)基礎(chǔ)。目前基因組調(diào)研的主要手段是高通量二代測序技術(shù)。本研究通過Illumina數(shù)據(jù),利用生物信息學(xué)手段,首次對(duì)腹足綱淡水軟體動(dòng)物田螺科貝類中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺等代表性物種的基因組大小、雜合度進(jìn)行評(píng)估分析。4種貝類的基因組大小分別為1 353、896、991、1 196 Mb,說明田螺科物種基因組大小存在差異,但與已公布大多數(shù)腹足綱物種的基因組大小接近[19-21]。腹足綱軟體動(dòng)物基因組大小存在豐富的多樣性,如在已測腹足綱中,基因組最大的物種是地中海芋螺(Conusventricosus),基因組大小為 3.6 G;霸王蓮花青螺(Lottiagigantea)基因組最小(359.5 Mb)[8],這2個(gè)物種基因組大小差異約為10倍?;蚪M是反映生物物種遺傳信息的重要指標(biāo),在物種進(jìn)化和適應(yīng)中起著重要作用,其大小與生物參數(shù)有關(guān),如細(xì)胞大小、抗逆性和個(gè)體生長等。在本研究中,田螺科貝類生長個(gè)體大小與基因組大小沒有直接的相關(guān)性。

    前人將一個(gè)物種雜合度高于0.5%定義為高雜合度物種[22],在本研究中,中華園田螺、銅銹環(huán)棱螺、多棱角螺、湄公螺4種貝類的雜合度分別約0.77%、1.35%、2.41%、0.91%,均大于0.5%,說明這4種田螺科貝類是高雜合率物種,這可能與田螺科物種的特殊性有關(guān)。田螺科貝類屬于雌雄異體,且生境基本相同,不同種群間基因交流可能是造成雜合率高的原因之一。這在Wang等應(yīng)用線粒體基因組構(gòu)建系統(tǒng)發(fā)育樹的研究中有所顯示[23],說明田螺科存在著種間基因交流,而部分種間是否存在生殖隔離還存在爭議,需要進(jìn)一步研究。實(shí)踐證明,高雜合度基因組會(huì)導(dǎo)致組裝過程中同源區(qū)域的read無法合并,分支結(jié)構(gòu)過多,連續(xù)性降低,從而造成組裝的基因組偏大。在本研究中也存在這樣的問題,應(yīng)用二代測序數(shù)據(jù)從頭組裝結(jié)果大小均高于應(yīng)用K-mer評(píng)估的結(jié)果。另外,高雜合還導(dǎo)致組裝結(jié)果中Contig N50較小,造成GC平均深度及分布較小。以上結(jié)果說明高雜合度提高了對(duì)田螺科物種精細(xì)基因圖譜的構(gòu)建難度。在本研究中4個(gè)田螺科物種重復(fù)序列占全基因組較小,分別為26.2%、26.7%、27.8%、28.8%,與已發(fā)表的其他腹足綱物種接近[24-26],為低重復(fù)序列基因組(重復(fù)序列比例高于50%是高重復(fù)基因組)[27],但非洲大蝸牛(Achatinafulica)例外,其重復(fù)序列比例較大,占整個(gè)基因組的70%,前人研究發(fā)現(xiàn)其發(fā)生基因組復(fù)制事件[28]。本研究4個(gè)田螺科貝類基因組重復(fù)序列占比較小,推測較少發(fā)生全基因組復(fù)制事件。

    物種GC比例是評(píng)估、分析物種調(diào)研圖準(zhǔn)確度,評(píng)估后續(xù)基因組圖譜組裝難度的重點(diǎn)衡量標(biāo)準(zhǔn)之一[29]。在本研究中,田螺科4種貝類的GC含量分別為33.4%、33.7%、33.1%、33.8%,腹足綱其他物種GC含量在26.8%~45.5%的區(qū)間內(nèi),其中綠唇鮑(Haliotislaevigata)GC含量最小(26.8%)[30],印度洋熱液噴口蝸牛(Alviniconchamarisindica)最大(GC含量為45.5%)[31]。事實(shí)上腹足綱物種線粒體基因組中也有存在GC含量較低的規(guī)律,如福壽螺(Pomaceacanaliculata)、斑點(diǎn)果瓶螺(P.maculat)和皺紋盤鮑(Haliotisdiscushannai)的線粒體GC含量在33%~40%范圍內(nèi)[32],說明腹足綱物種核基因組和線粒體基因組均具有AT偏好性。研究發(fā)現(xiàn)物種在基因組圖譜構(gòu)建中,過高或過低的GC含量將導(dǎo)致測序錯(cuò)誤,GC貧乏或富集區(qū)通常會(huì)引起擴(kuò)增效果較差,并影響組裝數(shù)據(jù)的準(zhǔn)確性[33],且對(duì)準(zhǔn)確度的干擾高于完整性。這也是本研究4個(gè)貝類物種中基因組初步組裝質(zhì)量較差,以及其他貝類基因組序列組裝難度較高的因素之一[34]。

    SSR分子標(biāo)記具有多態(tài)性高、等位差異顯著的優(yōu)點(diǎn)[35],是遺傳研究中應(yīng)用廣泛的分子標(biāo)記[36-37]。本研究利用MISA分析方法對(duì)田螺科代表性4種螺進(jìn)行SSR特征分析,與植物的SSR主要是雙核苷酸和三核苷酸不同[22],本研究4種貝類中單、雙核核苷酸重復(fù)序列占比最高,分別達(dá)79.2%、79.9%、78.7%、80.9%,這與扁玉螺(Neveritadidyma)[35]、泥東風(fēng)螺(Babylonialutosa)[38]的研究結(jié)果相一致。從4種田螺科貝類的SSR重復(fù)基元來看,單核苷酸中SSR位點(diǎn)占比最多的為A/T,雙核苷酸重復(fù)基元中,AC/GT占比均最高,這與前人研究的結(jié)果[39]相一致。SSR的長度會(huì)影響其本身多態(tài)性,而4種螺所含SSR的長度絕大部分在10~24 bp 之間的,分別占71.4%、64.9%、68.4%、70.2%,推測田螺科物種可能受到強(qiáng)烈的趨同選擇壓力。本研究為田螺科物種基因圖譜構(gòu)建、物種分子鑒定提供了前期基礎(chǔ)。

    猜你喜歡
    田螺貝類雜合
    我國海水貝類養(yǎng)殖低碳效應(yīng)評(píng)價(jià)
    田螺姑娘的夏天
    甘藍(lán)型油菜隱性上位互作核不育系統(tǒng)不育系材料選育中常見的育性分離及基因型判斷
    種子(2021年3期)2021-04-12 01:42:22
    QuEChERS-液相色譜-高分辨質(zhì)譜法測定貝類中6種親脂性貝類毒素
    鮮美貝類可能暗藏毒素
    食品與生活(2019年8期)2019-10-30 12:13:09
    摸田螺
    田螺
    快樂語文(2017年33期)2017-12-06 09:19:46
    吃田螺
    小布老虎(2017年4期)2017-08-10 08:22:40
    從翻譯到文化雜合——“譯創(chuàng)”理論的虛涵數(shù)意
    雄激素可調(diào)節(jié)的腎臟近端腎小管上皮細(xì)胞靶向雜合啟動(dòng)子的優(yōu)化
    日本五十路高清| 久久久国产精品麻豆| 亚洲av熟女| 欧美日韩福利视频一区二区| 国产99久久九九免费精品| 在线观看66精品国产| 亚洲乱码一区二区免费版| 91麻豆精品激情在线观看国产| 黄频高清免费视频| 女人被狂操c到高潮| 真人做人爱边吃奶动态| 大型黄色视频在线免费观看| 国产高清视频在线观看网站| 岛国在线免费视频观看| 男插女下体视频免费在线播放| 国产激情久久老熟女| 91在线观看av| 在线观看www视频免费| 悠悠久久av| 在线永久观看黄色视频| 99久久精品热视频| 日本a在线网址| 美女 人体艺术 gogo| 亚洲人成网站高清观看| 久久久久久免费高清国产稀缺| 美女大奶头视频| 最近视频中文字幕2019在线8| 欧美人与性动交α欧美精品济南到| www.熟女人妻精品国产| 在线a可以看的网站| 日本 av在线| 久久精品国产亚洲av高清一级| 最新在线观看一区二区三区| 国产亚洲精品综合一区在线观看 | 毛片女人毛片| 久久热在线av| 国产午夜福利久久久久久| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 婷婷丁香在线五月| 国产精品国产高清国产av| 精品久久久久久,| 嫩草影院精品99| 三级毛片av免费| 国产激情欧美一区二区| 日韩欧美国产一区二区入口| 亚洲人成77777在线视频| 国产精品美女特级片免费视频播放器 | 久久久国产成人精品二区| 看黄色毛片网站| 国产一区二区三区视频了| 宅男免费午夜| 亚洲五月婷婷丁香| 窝窝影院91人妻| 老司机靠b影院| 亚洲男人天堂网一区| 又大又爽又粗| 黑人欧美特级aaaaaa片| 我的老师免费观看完整版| 国产成年人精品一区二区| 黄片大片在线免费观看| 久久午夜综合久久蜜桃| 麻豆久久精品国产亚洲av| 亚洲男人天堂网一区| 成人手机av| 91大片在线观看| 熟女电影av网| 长腿黑丝高跟| 99久久国产精品久久久| 99国产精品一区二区三区| 亚洲精品av麻豆狂野| 91麻豆av在线| 99在线人妻在线中文字幕| 国产亚洲精品久久久久5区| 岛国视频午夜一区免费看| 久久精品国产亚洲av香蕉五月| 午夜免费成人在线视频| 日本成人三级电影网站| 又黄又粗又硬又大视频| 精品国产乱子伦一区二区三区| 亚洲人成网站高清观看| 午夜精品一区二区三区免费看| 黄色丝袜av网址大全| 伦理电影免费视频| 国产成人欧美在线观看| 亚洲五月天丁香| 久久久精品欧美日韩精品| 亚洲乱码一区二区免费版| 手机成人av网站| 50天的宝宝边吃奶边哭怎么回事| 少妇的丰满在线观看| 18禁裸乳无遮挡免费网站照片| 女生性感内裤真人,穿戴方法视频| 我要搜黄色片| 欧美zozozo另类| 真人一进一出gif抽搐免费| 欧美色欧美亚洲另类二区| 久久久久精品国产欧美久久久| 国产1区2区3区精品| 国内精品一区二区在线观看| 日本免费一区二区三区高清不卡| 啦啦啦免费观看视频1| 国产精品日韩av在线免费观看| 成年人黄色毛片网站| 国产又黄又爽又无遮挡在线| 亚洲专区字幕在线| 亚洲 欧美 日韩 在线 免费| 日韩 欧美 亚洲 中文字幕| 久久久精品欧美日韩精品| 久久久精品国产亚洲av高清涩受| 精品第一国产精品| 少妇熟女aⅴ在线视频| 999精品在线视频| 美女免费视频网站| 国产av不卡久久| 久99久视频精品免费| 日日摸夜夜添夜夜添小说| 好男人电影高清在线观看| 国产精品九九99| x7x7x7水蜜桃| 国语自产精品视频在线第100页| 国产精品av视频在线免费观看| 中国美女看黄片| 欧美久久黑人一区二区| 国内少妇人妻偷人精品xxx网站 | 男女视频在线观看网站免费 | 亚洲av成人av| 一进一出抽搐动态| netflix在线观看网站| 亚洲自拍偷在线| av视频在线观看入口| 日韩欧美在线二视频| 女同久久另类99精品国产91| 淫秽高清视频在线观看| 美女扒开内裤让男人捅视频| 三级国产精品欧美在线观看 | 99精品久久久久人妻精品| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 久久精品91无色码中文字幕| 可以在线观看毛片的网站| 国产精品久久久久久亚洲av鲁大| cao死你这个sao货| 欧美久久黑人一区二区| 女人高潮潮喷娇喘18禁视频| av超薄肉色丝袜交足视频| 男女那种视频在线观看| 日本五十路高清| 最近最新中文字幕大全免费视频| 成人18禁高潮啪啪吃奶动态图| 亚洲一区二区三区色噜噜| 久久香蕉国产精品| 精品欧美一区二区三区在线| 精品国产亚洲在线| 久久久久久大精品| 伦理电影免费视频| 欧美又色又爽又黄视频| 国内毛片毛片毛片毛片毛片| 婷婷丁香在线五月| 久久精品91蜜桃| 国产精品九九99| 首页视频小说图片口味搜索| 黑人巨大精品欧美一区二区mp4| netflix在线观看网站| 男女午夜视频在线观看| 日韩大码丰满熟妇| 久久久久久久久中文| 天天躁狠狠躁夜夜躁狠狠躁| 毛片女人毛片| 老司机深夜福利视频在线观看| 午夜影院日韩av| av超薄肉色丝袜交足视频| 怎么达到女性高潮| 欧美在线黄色| 亚洲国产欧美一区二区综合| 男男h啪啪无遮挡| www日本黄色视频网| 99riav亚洲国产免费| 欧美最黄视频在线播放免费| 97人妻精品一区二区三区麻豆| 亚洲精华国产精华精| 一级毛片女人18水好多| 国产又色又爽无遮挡免费看| 国产熟女午夜一区二区三区| 757午夜福利合集在线观看| 欧美最黄视频在线播放免费| av欧美777| 制服诱惑二区| 成人高潮视频无遮挡免费网站| 亚洲激情在线av| 国产真实乱freesex| 亚洲国产日韩欧美精品在线观看 | 亚洲精品中文字幕在线视频| 欧美在线一区亚洲| 波多野结衣巨乳人妻| 国产一区二区在线观看日韩 | 色综合亚洲欧美另类图片| 欧美色视频一区免费| 欧美日韩国产亚洲二区| 十八禁网站免费在线| 国产日本99.免费观看| 色老头精品视频在线观看| 丰满人妻一区二区三区视频av | 国产亚洲欧美在线一区二区| 欧美成人免费av一区二区三区| 亚洲 欧美 日韩 在线 免费| 丝袜美腿诱惑在线| 高清在线国产一区| 在线视频色国产色| 可以在线观看的亚洲视频| 岛国在线观看网站| 变态另类丝袜制服| 亚洲国产精品久久男人天堂| 免费在线观看影片大全网站| 制服丝袜大香蕉在线| 国产人伦9x9x在线观看| 日韩成人在线观看一区二区三区| 真人一进一出gif抽搐免费| 欧美 亚洲 国产 日韩一| 国产91精品成人一区二区三区| 国产成人一区二区三区免费视频网站| 一边摸一边做爽爽视频免费| 午夜亚洲福利在线播放| 国产精品野战在线观看| 草草在线视频免费看| 18禁黄网站禁片午夜丰满| 亚洲国产精品合色在线| 久久久久久久午夜电影| 欧美大码av| 最好的美女福利视频网| 后天国语完整版免费观看| 久久精品亚洲精品国产色婷小说| 女生性感内裤真人,穿戴方法视频| 成熟少妇高潮喷水视频| 人人妻,人人澡人人爽秒播| 欧美成人一区二区免费高清观看 | 天天躁夜夜躁狠狠躁躁| 一本久久中文字幕| 少妇熟女aⅴ在线视频| 国产欧美日韩精品亚洲av| 婷婷丁香在线五月| 国产av在哪里看| 叶爱在线成人免费视频播放| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品一区二区www| 波多野结衣高清无吗| 国产午夜精品久久久久久| 91国产中文字幕| 久久午夜综合久久蜜桃| 一进一出抽搐gif免费好疼| 国产高清视频在线观看网站| 午夜激情福利司机影院| 熟妇人妻久久中文字幕3abv| 搡老岳熟女国产| 欧美日韩亚洲综合一区二区三区_| 国产麻豆成人av免费视频| 99精品在免费线老司机午夜| 亚洲 国产 在线| 视频区欧美日本亚洲| 成人av在线播放网站| 午夜精品一区二区三区免费看| 国产99白浆流出| 亚洲精品美女久久久久99蜜臀| 亚洲五月婷婷丁香| 一区二区三区国产精品乱码| 午夜福利在线在线| 狠狠狠狠99中文字幕| 欧美中文综合在线视频| 不卡一级毛片| 久久久水蜜桃国产精品网| 别揉我奶头~嗯~啊~动态视频| 国产激情欧美一区二区| 在线免费观看的www视频| 久久草成人影院| 久久精品人妻少妇| 好男人在线观看高清免费视频| 曰老女人黄片| 日本a在线网址| 日日爽夜夜爽网站| 亚洲国产欧洲综合997久久,| 午夜成年电影在线免费观看| 精品久久久久久,| 操出白浆在线播放| 熟女电影av网| 在线观看午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 久久人妻福利社区极品人妻图片| 精品高清国产在线一区| 88av欧美| 一本精品99久久精品77| 男女做爰动态图高潮gif福利片| 久久久久性生活片| 亚洲av成人不卡在线观看播放网| 欧美日韩福利视频一区二区| 亚洲中文字幕一区二区三区有码在线看 | 狠狠狠狠99中文字幕| 精品高清国产在线一区| 亚洲中文av在线| 九色国产91popny在线| 一a级毛片在线观看| 最新在线观看一区二区三区| 人人妻,人人澡人人爽秒播| 最近在线观看免费完整版| 久久国产精品影院| 午夜福利在线观看吧| 岛国在线观看网站| 岛国在线免费视频观看| 日本熟妇午夜| 性色av乱码一区二区三区2| 精品久久久久久久人妻蜜臀av| 在线免费观看的www视频| 亚洲国产精品久久男人天堂| 亚洲国产精品999在线| 三级毛片av免费| 久久婷婷成人综合色麻豆| www.www免费av| 亚洲av五月六月丁香网| 国产精品自产拍在线观看55亚洲| 麻豆一二三区av精品| 国产成年人精品一区二区| 黄色丝袜av网址大全| 男人舔女人下体高潮全视频| 熟女少妇亚洲综合色aaa.| 日韩欧美精品v在线| 国产真实乱freesex| 欧美日韩国产亚洲二区| 日本黄色视频三级网站网址| 亚洲av五月六月丁香网| 久久精品国产清高在天天线| 在线观看www视频免费| 国产亚洲av嫩草精品影院| 香蕉丝袜av| 一本久久中文字幕| 国产单亲对白刺激| 久久伊人香网站| 国产片内射在线| 国产免费男女视频| 日本熟妇午夜| 中文资源天堂在线| 午夜两性在线视频| 午夜福利高清视频| 久久久久久大精品| 亚洲av电影不卡..在线观看| 99久久99久久久精品蜜桃| 嫩草影视91久久| 九九热线精品视视频播放| 又爽又黄无遮挡网站| 一区二区三区激情视频| 90打野战视频偷拍视频| www国产在线视频色| 亚洲精品在线美女| 亚洲人成77777在线视频| 午夜成年电影在线免费观看| 变态另类成人亚洲欧美熟女| 国产精品久久久久久精品电影| 久久精品国产亚洲av高清一级| 国产精品亚洲av一区麻豆| 又紧又爽又黄一区二区| 岛国视频午夜一区免费看| 欧美黄色淫秽网站| 日韩中文字幕欧美一区二区| 99精品欧美一区二区三区四区| 18禁黄网站禁片午夜丰满| 久久久精品大字幕| 久久久久亚洲av毛片大全| 国产精品九九99| 老司机午夜福利在线观看视频| 老司机在亚洲福利影院| av有码第一页| www.熟女人妻精品国产| 麻豆久久精品国产亚洲av| 99国产精品一区二区蜜桃av| 波多野结衣巨乳人妻| 少妇被粗大的猛进出69影院| av片东京热男人的天堂| 黑人操中国人逼视频| 亚洲国产看品久久| 可以在线观看的亚洲视频| 久热爱精品视频在线9| 国内揄拍国产精品人妻在线| 欧美色视频一区免费| 亚洲成av人片在线播放无| 国产精品一及| 成人国语在线视频| 色哟哟哟哟哟哟| 精品久久久久久久毛片微露脸| 精品日产1卡2卡| 女同久久另类99精品国产91| 亚洲九九香蕉| 高潮久久久久久久久久久不卡| 成人特级黄色片久久久久久久| 国产一区二区激情短视频| 欧美色欧美亚洲另类二区| 我的老师免费观看完整版| 亚洲国产欧美一区二区综合| 亚洲欧洲精品一区二区精品久久久| 午夜两性在线视频| 制服丝袜大香蕉在线| 蜜桃久久精品国产亚洲av| 久久欧美精品欧美久久欧美| 禁无遮挡网站| 中文资源天堂在线| 国产精品久久久人人做人人爽| 可以免费在线观看a视频的电影网站| 欧美日韩黄片免| 国产三级黄色录像| 国产激情久久老熟女| 午夜激情av网站| 少妇被粗大的猛进出69影院| 香蕉国产在线看| 成人国产综合亚洲| 成在线人永久免费视频| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久亚洲av鲁大| 成人手机av| 女人爽到高潮嗷嗷叫在线视频| 久久天堂一区二区三区四区| 少妇的丰满在线观看| 国产高清视频在线播放一区| 午夜亚洲福利在线播放| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 夜夜夜夜夜久久久久| 中文亚洲av片在线观看爽| 国产亚洲欧美在线一区二区| 欧美另类亚洲清纯唯美| av天堂在线播放| 在线观看免费日韩欧美大片| 久久久久久久久中文| 91麻豆精品激情在线观看国产| 婷婷丁香在线五月| 夜夜看夜夜爽夜夜摸| 大型av网站在线播放| 99热这里只有是精品50| 一级a爱片免费观看的视频| 给我免费播放毛片高清在线观看| 色综合婷婷激情| 日韩三级视频一区二区三区| 精品第一国产精品| 国产三级在线视频| 黄片大片在线免费观看| 熟妇人妻久久中文字幕3abv| 亚洲熟妇熟女久久| 亚洲avbb在线观看| 免费在线观看影片大全网站| 午夜福利视频1000在线观看| 免费无遮挡裸体视频| 97人妻精品一区二区三区麻豆| 亚洲欧美一区二区三区黑人| 久久国产精品人妻蜜桃| 麻豆久久精品国产亚洲av| 男人的好看免费观看在线视频 | 国产三级在线视频| 一进一出抽搐动态| 老熟妇乱子伦视频在线观看| 欧美一区二区国产精品久久精品 | 少妇粗大呻吟视频| 久久久久久久午夜电影| 久久草成人影院| 熟女电影av网| 国产精品九九99| 可以免费在线观看a视频的电影网站| 午夜激情av网站| 1024香蕉在线观看| xxxwww97欧美| 欧美又色又爽又黄视频| 精品无人区乱码1区二区| 又紧又爽又黄一区二区| 欧美黑人巨大hd| 亚洲国产精品合色在线| 叶爱在线成人免费视频播放| 级片在线观看| 国产在线观看jvid| 99久久精品国产亚洲精品| 12—13女人毛片做爰片一| а√天堂www在线а√下载| 级片在线观看| 黄色女人牲交| 亚洲avbb在线观看| 少妇人妻一区二区三区视频| 99热6这里只有精品| 99久久久亚洲精品蜜臀av| 久久久久久久精品吃奶| 不卡av一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲av中文字字幕乱码综合| 日本精品一区二区三区蜜桃| 欧美国产日韩亚洲一区| 日韩成人在线观看一区二区三区| 免费在线观看亚洲国产| 欧美一区二区精品小视频在线| 给我免费播放毛片高清在线观看| 国产精品98久久久久久宅男小说| 香蕉丝袜av| 中文字幕av在线有码专区| 亚洲精品美女久久av网站| 欧美高清成人免费视频www| 亚洲精品中文字幕一二三四区| 国产亚洲av高清不卡| 久99久视频精品免费| 2021天堂中文幕一二区在线观| 国产激情久久老熟女| 深夜精品福利| 麻豆一二三区av精品| 狠狠狠狠99中文字幕| 亚洲国产高清在线一区二区三| 中文字幕精品亚洲无线码一区| 精品久久久久久久毛片微露脸| 国产精品 欧美亚洲| 少妇粗大呻吟视频| 国产单亲对白刺激| 亚洲人成网站在线播放欧美日韩| 久久久久九九精品影院| av免费在线观看网站| 舔av片在线| 久久人妻av系列| 亚洲国产精品sss在线观看| 午夜福利视频1000在线观看| 超碰成人久久| 久久精品影院6| 国产精品亚洲av一区麻豆| 不卡一级毛片| 久久人妻福利社区极品人妻图片| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 国产成人aa在线观看| 亚洲专区国产一区二区| 欧美性猛交黑人性爽| 精品久久久久久久久久免费视频| 日韩欧美精品v在线| 久99久视频精品免费| 国产视频一区二区在线看| 亚洲精品中文字幕一二三四区| 丁香六月欧美| 香蕉av资源在线| 两人在一起打扑克的视频| 亚洲 欧美 日韩 在线 免费| 一个人免费在线观看电影 | 一二三四社区在线视频社区8| 可以免费在线观看a视频的电影网站| 久久婷婷成人综合色麻豆| 制服丝袜大香蕉在线| 啦啦啦观看免费观看视频高清| 97人妻精品一区二区三区麻豆| 久久精品国产99精品国产亚洲性色| 亚洲国产精品sss在线观看| 成年女人毛片免费观看观看9| 国产一区在线观看成人免费| 婷婷丁香在线五月| 亚洲欧美日韩无卡精品| 可以免费在线观看a视频的电影网站| 欧美在线黄色| 丰满人妻熟妇乱又伦精品不卡| 久久天堂一区二区三区四区| 日本一本二区三区精品| 波多野结衣高清作品| 天堂av国产一区二区熟女人妻 | 三级毛片av免费| 久久亚洲精品不卡| 久久久久久久久中文| 又大又爽又粗| 久久国产乱子伦精品免费另类| 亚洲人成网站高清观看| 日日干狠狠操夜夜爽| 国产精品综合久久久久久久免费| 99热这里只有精品一区 | 免费搜索国产男女视频| 此物有八面人人有两片| 一本精品99久久精品77| 国产黄色小视频在线观看| 两个人的视频大全免费| 美女扒开内裤让男人捅视频| 久久久久久大精品| 亚洲中文字幕日韩| www国产在线视频色| 变态另类丝袜制服| 亚洲精品一卡2卡三卡4卡5卡| 国产成人一区二区三区免费视频网站| 91麻豆精品激情在线观看国产| 一本大道久久a久久精品| 真人一进一出gif抽搐免费| 一a级毛片在线观看| 人人妻人人看人人澡| 高清毛片免费观看视频网站| 亚洲自拍偷在线| 国产探花在线观看一区二区| 欧美在线一区亚洲| 午夜福利在线在线| 淫妇啪啪啪对白视频| 日韩av在线大香蕉| 午夜免费成人在线视频| 99热这里只有是精品50| 人成视频在线观看免费观看| 亚洲精品在线观看二区| 亚洲在线自拍视频| 男女视频在线观看网站免费 | 欧洲精品卡2卡3卡4卡5卡区| 成在线人永久免费视频| 三级国产精品欧美在线观看 | 国产亚洲av高清不卡| 精品久久蜜臀av无| 亚洲无线在线观看| 午夜两性在线视频| 色av中文字幕| 亚洲人成77777在线视频| 妹子高潮喷水视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品爽爽va在线观看网站| 手机成人av网站| 99热6这里只有精品| 免费电影在线观看免费观看| 亚洲一区二区三区不卡视频| 男插女下体视频免费在线播放| 午夜福利欧美成人| 日韩欧美在线乱码|