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

    新型牡蠣相關(guān)圓環(huán)病毒基因組的鑒定

    2022-08-23 07:36:32楊李玲郭迎香位紅穎方藝菲姜敬哲
    南方水產(chǎn)科學 2022年4期
    關(guān)鍵詞:進化樹牡蠣圓環(huán)

    楊李玲,郭迎香,位紅穎,王 萌,方藝菲,朱 鵬,姜敬哲

    1. 天津農(nóng)學院,天津 300384

    2. 中國水產(chǎn)科學研究院南海水產(chǎn)研究所/農(nóng)業(yè)農(nóng)村部南海漁業(yè)資源開發(fā)利用重點實驗室,廣東 廣州 510300

    3. 廣州美格基因科技有限公司,廣東 廣州 510005

    4. 從化區(qū)農(nóng)業(yè)農(nóng)村局,廣東 廣州 510925

    5. 上海美吉生物醫(yī)藥科技有限公司,上海 201321

    6. 上海海洋大學,上海 201306

    牡蠣屬于軟體動物門、雙殼綱、珍珠貝目,廣泛分布于世界各地,其營養(yǎng)價值高,是全球養(yǎng)殖范圍最廣、產(chǎn)量最高的海水養(yǎng)殖品種。中國作為最大的牡蠣生產(chǎn)國,2019 年牡蠣產(chǎn)量達到 8 259×104t,占世界總產(chǎn)量的85.3% (數(shù)據(jù)來源:http://www.fao.org)。然而,受環(huán)境變化、陸源污染及種質(zhì)退化等因素影響,牡蠣病害的發(fā)生日趨頻繁,對其養(yǎng)殖業(yè)的健康可持續(xù)發(fā)展造成負面影響[1]。

    1972年Farley等[2]首次報告了牡蠣病毒性疾病病例,證明在高溫條件下牡蠣死亡事件與皰疹病毒存在關(guān)聯(lián),并將這種病毒命名為1型牡蠣皰疹病毒 (Ostreid Herpesvirus-1, OsHV-1)。該病毒在長牡蠣 (Crassostrea gigas) 中具有高致病性,同時對其他多種雙殼類也造成嚴重影響[3-4]。除OsHV-1外,其他報道過的與牡蠣相關(guān)的病毒還包括:乳多空病毒科 (Papovaviridae) 病毒,可能導致牡蠣“卵囊炎”,引起卵子和配子細胞肥大;鰓壞死病毒,屬虹彩病毒科 (Iridoviridae),可能是導致20世紀60 年代后期葡萄牙牡蠣 (C. angulata) 種群大規(guī)模死亡的主要原因;另外,披膜病毒 (Togaviridae)、呼腸孤病毒 (Reoviridae) 和小 RNA 病毒 (Picornavridae) 在貝類宿主中也有多次報道[5]。牡蠣等濾食性貝類在過濾大量海水獲取有機顆粒物營養(yǎng)的同時,也會將水體中的其他物質(zhì)如病毒、污染物等帶入體內(nèi),并通過外套膜進行濃縮。因此,如果牡蠣養(yǎng)殖水域受到人類活動影響,牡蠣體內(nèi)可能還會富集諸如諾如病毒 (Norovirus, NV)、甲型肝炎病毒 (Hepatitis A Virus, HAV) 和星狀病毒 (Astrovirus, AV) 等引起人類消化道疾病的食源性病毒,但這些病毒并不是牡蠣的病原[6]。病毒基因組是疾病診斷、流行病檢測和病原生物學研究的基礎(chǔ)。由于無脊椎動物沒有能夠穩(wěn)定傳代的細胞系,除了OsHV-1和一些食源性病毒外,其他牡蠣相關(guān)病毒的鑒定證據(jù)主要基于病理學和電子顯微鏡觀察,很少能夠拿到基因組序列[7]。因此,基因組數(shù)據(jù)的匱乏是牡蠣病害研究進展緩慢的一個主要原因。高通量測序和宏病毒組技術(shù)的運用,克服了傳統(tǒng)病毒學研究對宿主細胞培養(yǎng)的依賴,極大提高了新病毒的發(fā)現(xiàn)和鑒定效率。近些年,已在海洋環(huán)境、無脊椎動物、低等脊椎動物及人類等相關(guān)研究中得到了廣泛應(yīng)用,鑒定到大量的新型病毒[8-11],極大地擴展了人們對病毒世界的認知。

    本研究基于前期工作所獲得的華南沿海多地養(yǎng)殖牡蠣宏病毒組數(shù)據(jù),進一步運用生物信息學工具對其進行深入挖掘,發(fā)現(xiàn)并鑒定到5條新型牡蠣相關(guān)圓環(huán)病毒基因組,為促進牡蠣病毒及貝類病害相關(guān)研究提供參考。

    1 材料與方法

    1.1 數(shù)據(jù)來源與篩選

    課題組前期工作通過對中國華南沿海多地養(yǎng)殖牡蠣進行宏病毒組測序,得到約25億條測序原始序列 (reads)[12]。對測序數(shù)據(jù)進行 Fastp (V0.20.0)[13]質(zhì)控去除低質(zhì)量和接頭序列,用Megahit (V1.2.9)[14-15]將reads組裝成重疊群 (contig) 片段,以NCBI非冗余蛋白質(zhì)數(shù)據(jù)庫為參考,用Diamond[16](V0.9.24.125) 進行比對注釋,并用 Megan 6[17]進一步對注釋結(jié)果進行分類,篩選其中5條基因組完整的、疑似為圓環(huán)病毒科 (Circoviridae) 的病毒序列 (Contig ID:ML2-k141_11943、ML2-k141_27933、QZd1-k141_2132、T5S3-k141_267932、ZHd1-k141_220676)進一步深入分析。

    1.2 開放閱讀框預(yù)測與比對

    基于 NCBI ORF finder (https://www.ncbi.nlm.nih.gov/orffinder/) 進行開放閱讀框 (Open reading frame, ORF) 預(yù)測,其中 Minimal ORF length (nt) 選取 75,Genetic code選取 standard,ORF start codon to use 選擇 Any sense codon,忽略 nested ORFs。用預(yù)測到的ORF進行Blastp比對,其中數(shù)據(jù)庫選擇nr蛋白數(shù)據(jù)庫,e為0.000 001。對于有比對結(jié)果的ORF,取一致性最高的蛋白用tblastn方法反向比對本研究中病毒的基因組序列,驗證5條ORF序列是否完整,并用DNAstar進行翻譯,得到最終的蛋白序列。

    1.3 基因組進化樹的構(gòu)建及比較基因組分析

    將5條病毒全基因組序列提交至VIPTree(V1.9)[18](https://www.genome.jp/viptree/),選擇核酸類型為ssDNA,選擇宿主病毒類型為Eukaryote,選擇基因預(yù)測的遺傳密碼為Eukaryotic,構(gòu)建基因組進化樹。利用MAFFT[19]將5條病毒基因組序列與近源的病毒基因組序列進行比對,并使用TrimAI[20]進一步修剪以去除對齊不明確的區(qū)域,利用 MEGA[21]構(gòu)建 Neighbor-joining Tree,設(shè)置 Bootstrap Replications為 1 000,選擇 substitution model為 p-distance,設(shè)置 Site Coverage Cutoff為50%。

    根據(jù)與5條病毒復(fù)制酶蛋白一致性最高的參考蛋白序列 (YP_006281010.1),計算參考蛋白與5條病毒序列之間的 Amino Acid Identity (AAI)。下載YP_006281010.1對應(yīng)的基因組序列 (NC_017843.2),利用 Orthologous Average Nucleotide Identity Tool (OAT) 對 NC_017843.2 和 5 條病毒序列計算兩兩序列的 Average Nucleotide Identity (ANI),得到兩兩序列的ANI,使用R的ggcorrplot包對數(shù)值進行可視化。同時在VIPtree中進行基因組兩兩比較分析:設(shè)置 positioning of each sequence 為auto,colors為 default,ticks of genomes為 shown,vertical dashed lines為 shown,representation of tBLASTx hits為 normal。

    1.4 復(fù)制酶蛋白進化樹分析

    因5條病毒復(fù)制酶蛋白序列Blatsp比對的結(jié)果存在重疊,因此僅以其中1條序列(ML2-k141_11 943)比對結(jié)果前50的序列和5條序列一起構(gòu)建系統(tǒng)進化樹,建樹過程參考1.3。

    1.5 復(fù)制酶蛋白結(jié)構(gòu)域分析

    在NCBI上下載ARI44308.1[22],用Jalview對ARI44308.1、YP_006281010.1和這5條病毒復(fù)制酶蛋白序列進行結(jié)構(gòu)位點分析,使用ClustalW進行比對,下載保存SVG格式文件。用SMART (http://smart.embl-heidelberg.de/) 軟件進行結(jié)構(gòu)域預(yù)測,同時采用Swiss model進行三維結(jié)構(gòu)預(yù)測。

    1.6 病毒豐度分析

    為了計算每個病毒的相對豐度值,首先合并5 條病毒基因組序列用 salmon (V0.13.1) index[23]命令生成參考基因組數(shù)據(jù)集,然后再利用salmon quant[23]命令將全部牡蠣病毒組文庫的clean reads逐一映射 (Mapping) 到參考基因組上,統(tǒng)計映射上的reads數(shù),根據(jù)調(diào)整后的FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) 計算公式計算各病毒的相對豐度值。

    式中:G為比對的reads數(shù);T為測序reads總數(shù)(×106);L為基因組長度 (kb)。

    2 結(jié)果

    2.1 病毒注釋結(jié)果

    通過對中國華南沿海多地養(yǎng)殖牡蠣樣品進行宏病毒組測序,消除低質(zhì)量序列并組裝,共獲得3 375 091 條非冗余 contig序列,進一步與 NCBI的NR庫進行blastp比對和分類注釋最終獲得731 536條疑似病毒來源的contig,篩選其中5條基因組相對完整的、被注釋為圓環(huán)病毒的基因組序列進一步分析,具體信息見表1。

    表1 牡蠣相關(guān)圓環(huán)病毒基因組信息Table 1 Genome information of oyster-related circoviruses

    2.2 基因組水平進化樹分析

    為了明確病毒的分類,用5條基因組序列與全部單鏈DNA病毒一起構(gòu)建基因組水平的進化樹(VIPTree,圖1)?;蚪M樹中內(nèi)圈顏色代表病毒科(Family),主要包括圓環(huán)病毒科、雙生病毒科(Geminiviridae)、類雙生病毒科 (Genomovirdae) 和細項圈病毒科 (Anellovridae)。從科水平分類看,圓環(huán)病毒科在整個單鏈DNA基因組樹中占據(jù)絕大部分。候選的5條序列均聚在了一個小的分支 (紅色五角星標識,圖1),且在圓環(huán)病毒大分支內(nèi),說明它們屬于圓環(huán)病毒科成員。圖1的最外圈代表病毒的宿主,包括脊索動物門、節(jié)肢動物門、綠藻門和軟體動物門。從宿主來源看,候選病毒分支周圍主要是節(jié)肢動物宿主和其他來源宿主,而節(jié)肢動物與軟體動物同屬無脊椎動物。

    圖1 單鏈DNA病毒基因組進化樹 (VIPtree)Fig. 1 Phylogenetic tree of single-stranded DNA virus genome (VIPtree)

    為了解5條病毒序列的進化關(guān)系,將5條病毒基因組序列與近源的病毒基因組序列進行比對。5 條病毒序列與 NC 017843.2、NC 027797.1、NC 030466.1均聚在了一支,其中除了NC 030468.1來源于海水樣品外,NC 017843.2、NC 027797.1、NC 030466.1分別來源于橈足動物、珊瑚、大猩猩樣品(圖2)。該聚類結(jié)果與5條候選序列的來源背景基本一致 (海洋、動物)。值得注意的是,QZd1-k141_2132與T5S3-k141_267932在全基因組水平關(guān)系最接近。

    圖2 牡蠣相關(guān)圓環(huán)病毒基因組與近源基因組序列構(gòu)建的系統(tǒng)發(fā)育樹Fig. 2 Phylogenetic tree of five oyster-related circoviruses and relative viral genome sequences

    2.3 全基因組序列和復(fù)制酶蛋白序列一致性分析

    用候選序列中預(yù)測到的復(fù)制酶蛋白序列在NCBI上進行Blastp比對,發(fā)現(xiàn)5條病毒序列均與YP_006281010.1的蛋白序列相似性較高。該序列注釋為圓環(huán)病毒復(fù)制酶蛋白,該蛋白來源于一種海洋節(jié)肢動物——夏唇角水蚤 (Labidocera aestiva)。用其基因組序列 (NC_017843.2) 和5條序列進行兩兩比較,NC_017843.2基因組序列與5條候選序列的ANI介于0%~72%,其中最高的是QZd1-k141_2132 (63%),最低的是 ML2-k141_27933 (0%,圖3)。5條病毒序列兩兩間的一致性為0%~94%,其中QZd1-k141_2132與T5S3-k141_267932的最高(94%),ML2-k141_27933和其余4條病毒序列均為0%。根據(jù)ICTV規(guī)定的圓環(huán)病毒科在種水平劃分閾值80%來看,5條序列可能屬于4個不同的新種,且均與NC_017843.2不同,即QZd1-k141_2132與T5S3-k141_267932可能屬于同一種,而其余3條序列屬于另外3個種。復(fù)制酶蛋白序列一致性分析結(jié)果一致。

    圖3 圓環(huán)病毒基因組序列間的ANI (a) 和復(fù)制酶蛋白序列間的AAI (b)Fig. 3 ANI value among genomic sequences (a) and AAI value among replication proteins of circoviruses (b)

    2.4 比較基因組分析

    比較基因組分析表明,夏唇角水蚤圓環(huán)病毒基因組序列 (NC_017843.2) 與5條候選基因組序列均在復(fù)制酶蛋白區(qū)域 (即YP_006281010.1) 具有相似性 (圖4),但在另一個ORF,即疑似衣殼蛋白區(qū)域,均未見明顯的相似性??梢姾蜻x病毒的衣殼蛋白屬于新型的衣殼蛋白,與牡蠣宿主細胞受體的特異性識別相關(guān)。

    圖4 牡蠣相關(guān)圓環(huán)病毒比較基因組分析圖Fig. 4 Comparative genome analysis of oyster-related circoviruses

    2.5 Neighbor-joining 進化樹構(gòu)建

    因5條病毒復(fù)制酶蛋白序列Blatsp比對的結(jié)果存在重疊,因此僅以其中1條序列(ML2-k141_11943) 比對結(jié)果的前50的序列和5條序列一起構(gòu)建系統(tǒng)進化樹。ML2-k141_11943、ZHd1-k141_220676、QZd1-k141_2132、T5S3-k141_267932 這4個病毒和YP_006281010.1、GAC77785.1、QNG41123.1、AXQ66182.1形成了一個分支,其中YP_006281010.1、QNG41123.1和 AXQ66182.1的來源是動物宏基因組 (具體物種不明),GAC77785.1來源于海洋沉積物 (圖5)。ML2-k141_27933序列與AXH76045.1、AXH74760.1形成了一個分支,該序列來源同樣是動物宏基因組 (具體物種不明)。

    圖5 牡蠣相關(guān)圓環(huán)病毒復(fù)制酶蛋白進化樹Fig. 5 Phylogenetic tree of replication proteins of oyster-related circoviruses

    2.6 功能分析

    為查明復(fù)制酶ORF中是否有相關(guān)功能位點和結(jié)構(gòu)域,用ARI44308.1和YP_006281010.1作為參考序列與5條候選序列進行全局比對 (圖6)。雖然5條候選序列在典型環(huán)狀病毒的復(fù)制酶結(jié)構(gòu)域(HUH phage replication sub-family) 處與參考蛋白序列差異較大。但SMART分析發(fā)現(xiàn),除ZHd1-k141_220676外,余下的4條候選序列與2條參考序列都鑒定到明確的HUH結(jié)構(gòu)域。復(fù)制酶蛋白結(jié)構(gòu)域三維結(jié)構(gòu)預(yù)測表明,5條候選序列,尤其是ML2-k141_11943、ML2-k141_27933和T5S3-k141_267932,與YP_006281010.1復(fù)制酶蛋白的空間結(jié)構(gòu)存在較高的相似性 (圖7)。

    圖6 牡蠣相關(guān)圓環(huán)病毒復(fù)制酶蛋白序列比對Fig. 6 Sequence alignment of oyster-related circoviruses

    圖7 牡蠣相關(guān)圓環(huán)病毒復(fù)制酶蛋白三維結(jié)構(gòu)預(yù)測圖Fig. 7 Three-dimensional structure prediction of replication protein of oyster-related circoviruses

    2.7 病毒豐度分析

    為了計算每個病毒的相對豐度值,用salmon統(tǒng)計映射上的reads數(shù)。在牡蠣測序文庫中均檢測到了5條備選病毒序列的reads,F(xiàn)PKM豐度各不相同 (表2)。整體來看,T5S3-k141_267932的分布最廣泛,ML2-k141_11943在特定樣品中的豐度最高。具體來看,ML2-k141_11943和ML2-k141_27933 在 ML1 中豐度最高 (148.67 和 2 080.33);QZd1-k141_2132和T5S3-k141_267932在ChYJ1509Da中豐度最高 (542.87和76.93);ZHd1-k141_220676在 ChZH1511Da 豐度最高 (55.17)。

    表2 牡蠣相關(guān)圓環(huán)病毒豐度分析Table 2 Analysis of abundance of oyster-related circovirus virus

    3 討論

    近年來,我國沿海地區(qū)養(yǎng)殖貝類大規(guī)模死亡事件頻發(fā),涉及的貝類品種、養(yǎng)殖面積和區(qū)域范圍不斷擴大,使貝類養(yǎng)殖業(yè)遭受重大經(jīng)濟損失[24]。傳統(tǒng)病毒學研究嚴重依賴宿主細胞的培養(yǎng),由于缺乏可穩(wěn)定傳代的無脊椎動物細胞培養(yǎng)技術(shù),加之相關(guān)研究投入有限,貝類病毒研究進展十分緩慢[25]。由于高通量測序和病毒組技術(shù)的應(yīng)用,近年幾項海洋無脊椎動物病毒研究成果的發(fā)表,加深了人們對于海洋動物病毒的了解。2016年,Shi等[26]報道了關(guān)于無脊椎動物RNA病毒組的重要研究成果,對包括多種貝類生物在內(nèi)的約400種無脊椎動物 (包括牡蠣) 組織混合后的約200個文庫進行了測序,發(fā)現(xiàn)了1 400多種新型RNA病毒,極大地豐富了RNA病毒的多樣性和種類。Rosani和Gerdol[27]從長牡蠣 (C. gigas)和地中海貽貝 (Mytilus galloprovincialis) 宿主轉(zhuǎn)錄組數(shù)據(jù)中鑒定了7個幾乎完整的RNA病毒基因組,屬于小核糖核酸病毒目 (Picornavirales),其中有6種是全新的病毒。Ian等[28]對健康的與被感染的海星進行了比較研究,鑒別出1種細小病毒科 (Parvoviridae) 的疑似致病病原,并且明確其在浮游生物和海洋沉積物中也廣泛存在。Cui[29]從不同海域共采集了3門6綱58種的海洋無脊椎動物樣品 (包括牡蠣),并利用宏轉(zhuǎn)錄組的測序方法對不同物種的RNA病毒組進行研究,鑒定出涵蓋9個病毒家族 (Durnavirales、Totiviridae、Bunyavirales、Chuviridae、Picornavirales、Flaviviridae、Hepelivirales、Solemoviridae、Tombusviridae) 的363個RNA病毒,其中的315個與已知的RNA病毒相似度較低,被認為是全新的RNA病毒,其中牡蠣來源病毒4個。牡蠣等海洋無脊椎動物體內(nèi)可能存在著非常豐富的病毒種類。但前述研究多基于轉(zhuǎn)錄組測序數(shù)據(jù),發(fā)現(xiàn)的新病毒大部分為RNA病毒,對DNA病毒的相關(guān)報道則十分有限。

    圓環(huán)病毒科病毒是已知基因組最小的一類動物致病性DNA病毒,呈二十面體對稱球形,無囊膜,直徑為17~22 nm,基因組由大約1.7~2.3 kb的單鏈環(huán)狀DNA組成。該家族成員種類繁多,可感染各類動物宿主,包括蝙蝠、嚙齒動物、畜禽動物和人類等,其中了解比較清楚的是感染畜禽動物和人類的種類,如雞傳染性貧血病毒 (Chicken Infectious anemia)、豬圓環(huán)病毒 (Porcine circovirus)、鴿圓環(huán)病毒 (Pigeon circovirus)等[30-31]。截止到本研究前,在牡蠣中尚未有圓環(huán)病毒的報道。經(jīng)查詢Virus-Host數(shù)據(jù)庫確認 (數(shù)據(jù)截止至2021年5月),目前已知的所有圓環(huán)病毒宿主都是動物界內(nèi)的兩側(cè)對稱動物,加之基因組和蛋白序列進化樹表明,它們與夏唇角水蚤圓環(huán)病毒最為接近,同屬無脊椎動物宿主范疇。這意味著這些圓環(huán)病毒不可能是牡蠣體內(nèi)的某些共生細菌或原生動物的病毒。另一方面,圓環(huán)病毒的宿主特異性是由其衣殼蛋白序列決定的[32]。而本研究發(fā)現(xiàn)5條病毒基因組序列中,僅在ZHd1-k141_220676序列中比對到了疑似衣殼蛋白序列,且與水蚤圓環(huán)病毒衣殼蛋白無相似性。因此,這5個新型圓環(huán)病毒的宿主是水蚤等節(jié)肢動物的可能性不高。當然,現(xiàn)有證據(jù)尚不足以確認牡蠣是其真正宿主,還需要結(jié)合水體浮游蚤類篩查、擴大養(yǎng)殖牡蠣調(diào)查范圍及牡蠣人工感染實驗來進一步明確其真正的宿主。

    根據(jù)ICTV標準,同一圓環(huán)病毒科物種的成員在其整個基因組中應(yīng)該具有大于80%的核苷酸同一性[33]。本研究中的5條候選序列與已知圓環(huán)病毒在整個基因組水平的相似性低于80%,在衣殼蛋白方面未觀察到明顯的相似性。因此,本研究鑒定的5個病毒可能與現(xiàn)有圓環(huán)病毒科家族成員有較大差異,具體的分類界定還有待發(fā)現(xiàn)更多家族成員序列后的分析。

    病毒是海洋中最豐富的生命體,以貝類為代表的軟體動物門也是海洋中最大的動物類群。但目前對于兩者的交叉領(lǐng)域——貝類病毒的了解仍十分有限。宏病毒測序方法已經(jīng)在眾多生物和環(huán)境樣品中得到廣泛應(yīng)用,突顯了高通量測序技術(shù)在發(fā)現(xiàn)新病毒方面的巨大潛力[12]。然而,病毒組在未知序列解析和新病毒分類方面還有很多不足尚待解決。本研究發(fā)現(xiàn)的5個牡蠣圓環(huán)病毒的分類地位、家族成員以及宿主范圍、致病性、流行情況等方面還有很多工作要做。后續(xù)研究應(yīng)該在發(fā)揮病毒組技術(shù)優(yōu)勢的基礎(chǔ)上,同時配合以經(jīng)典的病毒學研究方法,獲得更為全面和深入的數(shù)據(jù),進而為提升為牡蠣等貝類病害的防控水平、支撐相關(guān)產(chǎn)業(yè)發(fā)展發(fā)揮作用。

    猜你喜歡
    進化樹牡蠣圓環(huán)
    加權(quán)全能量最小的圓環(huán)形變
    基于心理旋轉(zhuǎn)的小學生物進化樹教學實驗報告
    告別自汗用牡蠣,四季都輕松
    豬圓環(huán)病毒病的發(fā)生、診斷和防治
    一例鴨圓環(huán)病毒病的診斷
    常見的進化樹錯誤概念及其辨析*
    生物學通報(2021年9期)2021-07-01 03:24:44
    圓環(huán)上的覆蓋曲面不等式及其應(yīng)用
    艾草白粉病的病原菌鑒定
    曇石山文化的牡蠣器
    大眾考古(2015年6期)2015-06-26 08:27:16
    層次聚類在進化樹構(gòu)建中的應(yīng)用
    在线免费观看的www视频| 亚洲av成人不卡在线观看播放网| 欧美日本亚洲视频在线播放| 在线天堂中文资源库| а√天堂www在线а√下载| 美女 人体艺术 gogo| 成人影院久久| 一级黄色大片毛片| 日韩欧美一区二区三区在线观看| 性色av乱码一区二区三区2| 视频区图区小说| 久久久久久久久免费视频了| 久热爱精品视频在线9| 美女午夜性视频免费| 国产精品爽爽va在线观看网站 | 久久精品国产清高在天天线| 免费观看人在逋| 无人区码免费观看不卡| 午夜日韩欧美国产| 日本五十路高清| 亚洲久久久国产精品| 在线观看午夜福利视频| 精品第一国产精品| 久热爱精品视频在线9| 午夜a级毛片| 欧美日韩乱码在线| 色综合欧美亚洲国产小说| 香蕉久久夜色| 亚洲少妇的诱惑av| 男女之事视频高清在线观看| 少妇的丰满在线观看| 麻豆一二三区av精品| 999精品在线视频| 精品乱码久久久久久99久播| 欧美日韩视频精品一区| 一级毛片精品| 国产欧美日韩一区二区三区在线| 国产成人精品在线电影| e午夜精品久久久久久久| 国产成人免费无遮挡视频| 亚洲国产欧美日韩在线播放| 久久久久久大精品| 最新在线观看一区二区三区| 757午夜福利合集在线观看| 亚洲欧美日韩无卡精品| 国产精品 国内视频| 国产亚洲欧美精品永久| 欧美日韩av久久| 日韩精品免费视频一区二区三区| www.www免费av| 欧美大码av| 亚洲欧美激情在线| 亚洲久久久国产精品| 亚洲色图 男人天堂 中文字幕| 久久人人爽av亚洲精品天堂| 午夜精品国产一区二区电影| 俄罗斯特黄特色一大片| 欧美黄色淫秽网站| 国产精品一区二区精品视频观看| 欧美黄色片欧美黄色片| 美女 人体艺术 gogo| 在线国产一区二区在线| 亚洲第一av免费看| 国产黄色免费在线视频| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 成年人免费黄色播放视频| 99精品久久久久人妻精品| av国产精品久久久久影院| 亚洲国产欧美日韩在线播放| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 好看av亚洲va欧美ⅴa在| 国产男靠女视频免费网站| 99精国产麻豆久久婷婷| 欧美日韩亚洲国产一区二区在线观看| 日本黄色视频三级网站网址| 国产精品爽爽va在线观看网站 | 91av网站免费观看| 国产精品 欧美亚洲| 欧美在线黄色| 久久久精品欧美日韩精品| 成年人黄色毛片网站| 人人妻人人添人人爽欧美一区卜| 国产深夜福利视频在线观看| 亚洲七黄色美女视频| 老司机福利观看| 激情在线观看视频在线高清| 丰满饥渴人妻一区二区三| 免费看a级黄色片| 久久香蕉精品热| 18禁观看日本| 一区福利在线观看| 欧美激情极品国产一区二区三区| 午夜亚洲福利在线播放| 国产熟女午夜一区二区三区| 亚洲avbb在线观看| 色综合婷婷激情| 一边摸一边做爽爽视频免费| 亚洲自拍偷在线| 欧美日本中文国产一区发布| 啪啪无遮挡十八禁网站| 99久久精品国产亚洲精品| xxx96com| 岛国在线观看网站| 国产高清激情床上av| 这个男人来自地球电影免费观看| 国产精品野战在线观看 | 熟女少妇亚洲综合色aaa.| 国产精品免费视频内射| 国产午夜精品久久久久久| 久久这里只有精品19| 国产伦一二天堂av在线观看| 91麻豆精品激情在线观看国产 | 午夜影院日韩av| 51午夜福利影视在线观看| www.精华液| 亚洲专区国产一区二区| 成人免费观看视频高清| 欧洲精品卡2卡3卡4卡5卡区| 天天躁夜夜躁狠狠躁躁| 俄罗斯特黄特色一大片| www.精华液| 热re99久久精品国产66热6| 黄色视频不卡| av网站在线播放免费| 一二三四在线观看免费中文在| 18禁国产床啪视频网站| 老司机午夜十八禁免费视频| 亚洲欧美日韩另类电影网站| 亚洲人成77777在线视频| www国产在线视频色| 曰老女人黄片| 在线天堂中文资源库| 夜夜夜夜夜久久久久| 精品福利观看| 国产亚洲欧美在线一区二区| 亚洲avbb在线观看| 久久人人97超碰香蕉20202| 99国产综合亚洲精品| www国产在线视频色| 亚洲五月色婷婷综合| 天堂√8在线中文| 精品一区二区三区四区五区乱码| 亚洲欧美精品综合一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 香蕉丝袜av| 色综合欧美亚洲国产小说| 国产亚洲av高清不卡| 国产三级黄色录像| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 亚洲av日韩精品久久久久久密| 久久久国产成人精品二区 | 波多野结衣一区麻豆| 午夜日韩欧美国产| 亚洲,欧美精品.| 精品熟女少妇八av免费久了| 亚洲精品一卡2卡三卡4卡5卡| 久久精品影院6| 两人在一起打扑克的视频| 国产深夜福利视频在线观看| 日韩大码丰满熟妇| av电影中文网址| 精品一品国产午夜福利视频| 欧美人与性动交α欧美精品济南到| 亚洲一区二区三区不卡视频| 国产成年人精品一区二区 | 久久亚洲精品不卡| 国产成人影院久久av| 国产精品野战在线观看 | 黄色怎么调成土黄色| 久久人人精品亚洲av| 亚洲五月天丁香| 免费观看精品视频网站| 男女下面进入的视频免费午夜 | 黄色视频不卡| 正在播放国产对白刺激| 波多野结衣高清无吗| 亚洲黑人精品在线| 亚洲国产欧美网| 亚洲精品久久成人aⅴ小说| 欧美乱妇无乱码| 韩国av一区二区三区四区| 他把我摸到了高潮在线观看| 少妇的丰满在线观看| 一区二区三区激情视频| 99精国产麻豆久久婷婷| 日韩欧美国产一区二区入口| 亚洲第一av免费看| 操出白浆在线播放| 亚洲欧美日韩另类电影网站| 18禁美女被吸乳视频| xxx96com| 欧美黑人精品巨大| 亚洲欧美激情综合另类| 国产成人欧美| 最近最新免费中文字幕在线| 看黄色毛片网站| 精品一品国产午夜福利视频| e午夜精品久久久久久久| 国产视频一区二区在线看| 日韩成人在线观看一区二区三区| 欧美色视频一区免费| 国产亚洲精品综合一区在线观看 | 国产成人欧美在线观看| 亚洲专区中文字幕在线| 欧美在线一区亚洲| 成人亚洲精品av一区二区 | 亚洲精品国产精品久久久不卡| 精品免费久久久久久久清纯| 一进一出抽搐动态| 久久人人爽av亚洲精品天堂| 一区二区三区激情视频| 黄片小视频在线播放| 91国产中文字幕| 国产日韩一区二区三区精品不卡| 亚洲熟女毛片儿| 午夜福利,免费看| 一区二区三区精品91| 最新美女视频免费是黄的| 成人特级黄色片久久久久久久| 日本精品一区二区三区蜜桃| 成年人黄色毛片网站| 一级毛片精品| 亚洲少妇的诱惑av| 精品免费久久久久久久清纯| 精品午夜福利视频在线观看一区| 国产成人欧美| 久久久久九九精品影院| 欧美激情久久久久久爽电影 | 99热国产这里只有精品6| 亚洲午夜精品一区,二区,三区| 天堂√8在线中文| 高清av免费在线| 亚洲专区中文字幕在线| 久久九九热精品免费| 久久精品亚洲熟妇少妇任你| 国产精品爽爽va在线观看网站 | 制服诱惑二区| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 日日摸夜夜添夜夜添小说| 久热这里只有精品99| x7x7x7水蜜桃| 女人被躁到高潮嗷嗷叫费观| 超色免费av| 欧美乱色亚洲激情| 嫩草影院精品99| 人人澡人人妻人| 亚洲国产精品sss在线观看 | 欧美日韩亚洲高清精品| 亚洲五月婷婷丁香| 精品福利观看| 亚洲av熟女| 国产成人精品久久二区二区免费| 亚洲av成人av| 国产精品一区二区免费欧美| 亚洲视频免费观看视频| 老司机福利观看| 日韩精品青青久久久久久| 一边摸一边抽搐一进一小说| 亚洲午夜精品一区,二区,三区| 国产成人精品在线电影| 国产精品野战在线观看 | 亚洲午夜理论影院| 国产精品九九99| 国产高清videossex| 一级a爱片免费观看的视频| 亚洲人成电影免费在线| 妹子高潮喷水视频| 亚洲人成电影观看| 又紧又爽又黄一区二区| 久久久久久久精品吃奶| 亚洲自偷自拍图片 自拍| 国产午夜精品久久久久久| 国产高清国产精品国产三级| 两个人免费观看高清视频| 女人精品久久久久毛片| av天堂在线播放| 成年人黄色毛片网站| 九色亚洲精品在线播放| 亚洲一码二码三码区别大吗| 国产成+人综合+亚洲专区| 亚洲av成人av| 搡老乐熟女国产| 欧美在线黄色| 99久久综合精品五月天人人| 精品久久久久久,| 黄色片一级片一级黄色片| 亚洲三区欧美一区| 免费看十八禁软件| 一级黄色大片毛片| 视频在线观看一区二区三区| 手机成人av网站| 丝袜美腿诱惑在线| 91九色精品人成在线观看| 久久香蕉国产精品| 国产激情欧美一区二区| 亚洲熟妇中文字幕五十中出 | 成人免费观看视频高清| 国产1区2区3区精品| 免费在线观看完整版高清| 欧美另类亚洲清纯唯美| 日韩 欧美 亚洲 中文字幕| e午夜精品久久久久久久| 18禁观看日本| 国产国语露脸激情在线看| 色哟哟哟哟哟哟| 麻豆av在线久日| 国产精品国产高清国产av| 精品国产乱码久久久久久男人| 在线观看免费高清a一片| 深夜精品福利| 国产精品一区二区免费欧美| 国产成年人精品一区二区 | 90打野战视频偷拍视频| 午夜日韩欧美国产| 久久久精品欧美日韩精品| 国产1区2区3区精品| 国产一区二区三区综合在线观看| 脱女人内裤的视频| 国产三级黄色录像| 热re99久久国产66热| 国产一区二区三区综合在线观看| 亚洲 欧美一区二区三区| 久久精品国产清高在天天线| 最新美女视频免费是黄的| 国产蜜桃级精品一区二区三区| 大码成人一级视频| 美女午夜性视频免费| 国产成人欧美在线观看| 国产又爽黄色视频| 婷婷丁香在线五月| av在线播放免费不卡| 91麻豆av在线| 亚洲国产看品久久| 欧美另类亚洲清纯唯美| 亚洲国产看品久久| 免费久久久久久久精品成人欧美视频| 久久国产精品男人的天堂亚洲| 美女 人体艺术 gogo| 日韩免费av在线播放| 久久香蕉激情| 国产成人欧美| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 国产成人av教育| 18禁观看日本| 精品一区二区三区四区五区乱码| 免费看a级黄色片| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 中文字幕人妻熟女乱码| 大型黄色视频在线免费观看| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区三区在线| 精品国产一区二区久久| av天堂在线播放| 精品一区二区三区四区五区乱码| 天堂俺去俺来也www色官网| 国产精品爽爽va在线观看网站 | 国产av精品麻豆| 亚洲国产欧美网| 久久久精品欧美日韩精品| 亚洲av片天天在线观看| 韩国av一区二区三区四区| 亚洲精品国产精品久久久不卡| 日日夜夜操网爽| 国产一区二区在线av高清观看| 免费观看精品视频网站| 视频区欧美日本亚洲| 国产精品野战在线观看 | 欧美不卡视频在线免费观看 | 777久久人妻少妇嫩草av网站| 亚洲精品粉嫩美女一区| 首页视频小说图片口味搜索| 亚洲专区国产一区二区| 日本黄色日本黄色录像| 日本撒尿小便嘘嘘汇集6| 精品国产亚洲在线| 亚洲九九香蕉| 50天的宝宝边吃奶边哭怎么回事| 一级a爱片免费观看的视频| 亚洲色图 男人天堂 中文字幕| 无限看片的www在线观看| 成年人免费黄色播放视频| 高清在线国产一区| 日本五十路高清| 欧美乱妇无乱码| 9热在线视频观看99| 国产成人精品在线电影| 制服人妻中文乱码| 在线视频色国产色| av福利片在线| 一级a爱视频在线免费观看| 国产精品自产拍在线观看55亚洲| 人人澡人人妻人| 精品一品国产午夜福利视频| 国产亚洲欧美98| 自线自在国产av| a在线观看视频网站| 在线观看66精品国产| 超色免费av| 欧美亚洲日本最大视频资源| 欧美日韩视频精品一区| 国产深夜福利视频在线观看| 久久久久久亚洲精品国产蜜桃av| 国产又爽黄色视频| 搡老乐熟女国产| 国产成年人精品一区二区 | 午夜亚洲福利在线播放| 国产精品99久久99久久久不卡| 一边摸一边抽搐一进一小说| 悠悠久久av| 人人妻人人澡人人看| 女人高潮潮喷娇喘18禁视频| 大码成人一级视频| 午夜免费观看网址| 啦啦啦在线免费观看视频4| 国产精品日韩av在线免费观看 | 亚洲自拍偷在线| av网站在线播放免费| 丰满的人妻完整版| 免费观看精品视频网站| 国产av精品麻豆| 色婷婷久久久亚洲欧美| 国产成人影院久久av| 欧美亚洲日本最大视频资源| 日本vs欧美在线观看视频| 国产在线观看jvid| 男女高潮啪啪啪动态图| 99在线人妻在线中文字幕| 中文字幕高清在线视频| 欧美成人性av电影在线观看| 午夜两性在线视频| 黄色片一级片一级黄色片| 村上凉子中文字幕在线| 午夜久久久在线观看| 国产精品免费一区二区三区在线| 国产精品野战在线观看 | 亚洲欧美激情在线| 国产极品粉嫩免费观看在线| 丁香欧美五月| 久久影院123| 国产一区二区在线av高清观看| 咕卡用的链子| 亚洲国产精品一区二区三区在线| 午夜两性在线视频| 日本wwww免费看| 亚洲国产精品sss在线观看 | 国产又色又爽无遮挡免费看| 精品第一国产精品| 午夜福利在线免费观看网站| 老熟妇乱子伦视频在线观看| 如日韩欧美国产精品一区二区三区| 久久人人97超碰香蕉20202| 少妇的丰满在线观看| 黄色女人牲交| 亚洲国产毛片av蜜桃av| 女人精品久久久久毛片| av视频免费观看在线观看| 黄片播放在线免费| 淫秽高清视频在线观看| 欧美国产精品va在线观看不卡| 国产黄色免费在线视频| 日韩免费av在线播放| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 一区二区三区国产精品乱码| 久久精品国产亚洲av高清一级| 黑人巨大精品欧美一区二区mp4| 一级作爱视频免费观看| 美女大奶头视频| www.熟女人妻精品国产| 亚洲熟妇熟女久久| 老司机深夜福利视频在线观看| 麻豆av在线久日| 老熟妇仑乱视频hdxx| 中文字幕色久视频| 免费在线观看完整版高清| 在线观看免费高清a一片| 欧美成人性av电影在线观看| 不卡一级毛片| 久久久国产成人精品二区 | 中文字幕高清在线视频| 免费人成视频x8x8入口观看| 丝袜在线中文字幕| 成年女人毛片免费观看观看9| 国产精品国产av在线观看| 一级毛片精品| 69av精品久久久久久| 在线观看www视频免费| 国产成人一区二区三区免费视频网站| 国产又色又爽无遮挡免费看| 亚洲精品久久午夜乱码| 性色av乱码一区二区三区2| 久久香蕉精品热| 国产成人一区二区三区免费视频网站| 一区福利在线观看| 国产精品野战在线观看 | 老汉色∧v一级毛片| 免费不卡黄色视频| 91麻豆精品激情在线观看国产 | 一级片免费观看大全| 99久久久亚洲精品蜜臀av| 国产单亲对白刺激| 99精品欧美一区二区三区四区| 午夜日韩欧美国产| 色精品久久人妻99蜜桃| 亚洲全国av大片| 欧美日韩亚洲国产一区二区在线观看| 天堂√8在线中文| 不卡av一区二区三区| 超色免费av| 女人高潮潮喷娇喘18禁视频| av天堂久久9| 国产精品爽爽va在线观看网站 | 成人18禁高潮啪啪吃奶动态图| 午夜福利在线免费观看网站| 久久国产精品影院| 久久这里只有精品19| 男女做爰动态图高潮gif福利片 | 国产乱人伦免费视频| 久9热在线精品视频| 亚洲国产精品sss在线观看 | 一进一出好大好爽视频| 久久精品国产综合久久久| 一边摸一边抽搐一进一小说| 国产av在哪里看| 在线播放国产精品三级| 国产精品自产拍在线观看55亚洲| 一边摸一边做爽爽视频免费| 色婷婷av一区二区三区视频| 国产精品免费一区二区三区在线| 国产亚洲精品一区二区www| 高清毛片免费观看视频网站 | 午夜精品久久久久久毛片777| 久久国产精品男人的天堂亚洲| 日韩大尺度精品在线看网址 | 国产精品久久电影中文字幕| 成人免费观看视频高清| 久热爱精品视频在线9| 亚洲人成伊人成综合网2020| 母亲3免费完整高清在线观看| 国产av又大| 日韩欧美国产一区二区入口| 正在播放国产对白刺激| 久久中文字幕一级| 中文亚洲av片在线观看爽| 欧美丝袜亚洲另类 | 男女之事视频高清在线观看| 丝袜在线中文字幕| 亚洲国产中文字幕在线视频| 国产免费现黄频在线看| 一边摸一边做爽爽视频免费| 久久久国产成人精品二区 | 超碰成人久久| 麻豆国产av国片精品| 亚洲国产看品久久| 美女福利国产在线| 麻豆av在线久日| 最好的美女福利视频网| 淫妇啪啪啪对白视频| 亚洲一码二码三码区别大吗| 伊人久久大香线蕉亚洲五| 黄色丝袜av网址大全| 香蕉丝袜av| 每晚都被弄得嗷嗷叫到高潮| 夫妻午夜视频| 国产精品av久久久久免费| 久久精品91蜜桃| 成人三级做爰电影| 色综合站精品国产| 亚洲成人免费电影在线观看| 亚洲欧美精品综合久久99| 亚洲精品国产精品久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 操美女的视频在线观看| 亚洲熟妇中文字幕五十中出 | 亚洲av美国av| 久久影院123| 午夜亚洲福利在线播放| 悠悠久久av| 国产精品永久免费网站| 成在线人永久免费视频| 久久婷婷成人综合色麻豆| 日韩欧美在线二视频| 国产三级黄色录像| 午夜两性在线视频| 精品午夜福利视频在线观看一区| 欧美日韩中文字幕国产精品一区二区三区 | 精品国产乱子伦一区二区三区| 高清黄色对白视频在线免费看| 三上悠亚av全集在线观看| 法律面前人人平等表现在哪些方面| 免费高清在线观看日韩| 亚洲精品中文字幕在线视频| 熟女少妇亚洲综合色aaa.| 日韩视频一区二区在线观看| 亚洲精品在线美女| 亚洲男人的天堂狠狠| 极品教师在线免费播放| 亚洲精品国产区一区二| 日本三级黄在线观看| 中文字幕人妻丝袜制服| ponron亚洲| 国产人伦9x9x在线观看| 最近最新中文字幕大全免费视频| 黑人巨大精品欧美一区二区蜜桃| 久久久久久久久免费视频了| 电影成人av| 亚洲欧美日韩无卡精品| 99国产精品99久久久久| 嫁个100分男人电影在线观看|