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

    聯(lián)合二代和三代測序技術(shù)鑒定新冠病毒

    2023-12-21 08:49:06黃枝妙曾志偉陳明建翁育偉
    中國人獸共患病學(xué)報(bào) 2023年11期
    關(guān)鍵詞:變異基因組位點(diǎn)

    林 琦,黃枝妙,曾志偉,陳 煒,陳明建,翁育偉

    新型冠狀病毒感染(Coronavirus Disease 2019, COVID-19)是由新型冠狀病毒(Severe acute respiratory syndrome coronavirus 2,SARS-CoV-2)引起的急性呼吸道傳染病。SARS-CoV-2屬于β屬冠狀病毒,主要通過呼吸道飛沫和密切接觸等途徑傳播[1-3]。隨著全球COVID-19持續(xù)流行,SARS-CoV-2不斷發(fā)生變異并導(dǎo)致其傳播力和致病性發(fā)生改變。目前,Omicron變異株已成為國內(nèi)外主要流行毒株,雖然在流行過程中其致病力明顯減弱,但其潛伏期縮短(2~4 d)、傳播力變強(qiáng)、傳播速度變快,免疫逃逸能力增強(qiáng),對全球COVID-19防控提出了新挑戰(zhàn)[4-5]。

    及時(shí)發(fā)現(xiàn)各種SARS-CoV-2的新型變異株,了解其生物學(xué)及流行病學(xué)特征,是當(dāng)前COVID-19防控的重要手段,因此高通量測序技術(shù)大量應(yīng)用于新冠疫情防控工作中。目前,SARS-CoV-2高通量測序主要包括Illumina[6]、Ion Genestudio S5[7]和MGI[8]等二代測序平臺,以及Oxford nanopore[9]和Pacific bioscience[10]等三代測序平臺。其中Illumina和Oxford nanopore使用的更為廣泛。不同測序平臺的測序結(jié)果各有優(yōu)缺點(diǎn),如以O(shè)xford nanopore為代表的三代測序平臺的測序速度快,但其測序的精度尚有不足;而以Illumina為代表的二代測序平臺的測序精度高,但耗時(shí)較長,測序過程相對復(fù)雜[11-12]。在新冠防控工作具體實(shí)踐特別是在新冠疫情處置中,對測序要求既要“快速”又要“精準(zhǔn)”,因此需要兼顧測序的實(shí)效性和精確性。本研究采用三代和二代測序聯(lián)用的方式,測定了福建省2例BQ.1進(jìn)化分支病毒感染的輸入性COVID-19的病毒全基因組序列,并對其進(jìn)行了變異和進(jìn)化分析。在及時(shí)準(zhǔn)確判斷病毒變異類型的前提下,為研判和預(yù)測疫情提供了準(zhǔn)確的病原學(xué)信息。

    1 材料與方法

    1.1 病例基本信息 病例A和B均為男性,波蘭人,年齡分別為65歲和63歲,同為某船務(wù)代理公司船員。病例A已接種4劑次輝瑞mRNA新冠疫苗,病例B疫苗接種史不詳。兩病例于2022年10月6日從波蘭抵達(dá)上海并隔離,10月16日從上海乘機(jī)抵達(dá)福州市,隨后由專車接往寧德市,入住酒店后均未外出。兩病例10月18日核酸檢測結(jié)果為陽性,均無明顯發(fā)熱、咳嗽、流涕等癥狀。兩病例的鼻咽拭子標(biāo)本被送至福建省疾病預(yù)防控制中心進(jìn)行全基因組測序。

    1.2 核酸提取及real-time RT-PCR復(fù)核 使用磁珠法核酸提取試劑盒(貨號:T014,西安天隆)提取陽性病例鼻咽拭子標(biāo)本中的病毒RNA。利用SARS-CoV-2核酸檢測試劑盒(貨號:DA0932,廣州達(dá)安)對提取的RNA進(jìn)行real-time RT-PCR復(fù)核。留取復(fù)核后的陽性核酸于-80 ℃保存?zhèn)溆谩?/p>

    1.3 病毒全基因組測序 使用二步法長片段新冠病毒全基因組捕獲試劑盒(貨號:BK-WCoV024IITS,杭州柏熠)進(jìn)行病毒全基因組擴(kuò)增,按照Rapid Barcoding Kit(貨號:SQK-RBK004,Oxford Nanopore)操作步驟構(gòu)建三代測序文庫,使用Oxford Nanopore MinION MK1C三代測序儀進(jìn)行測序,采用柏熠新冠病毒全基因組分析系統(tǒng)對測序數(shù)據(jù)進(jìn)行質(zhì)控和病毒全基因組序列拼接。二代測序采用ULSEN超靈敏度新型冠狀病毒全基因組捕獲試劑盒(貨號:V-090418-1,北京微未來)進(jìn)行病毒全基因組擴(kuò)增,按照Nextera XT DNA Library Preparation Kit(貨號:FC-131-1096,Illumina)操作步驟構(gòu)建二代測序文庫,使用Illumina Miseq二代測序儀進(jìn)行測序,以Wuhan-Hu-1(NC_045512.2)基因組作為參考序列,使用CLC Genomics Workbench軟件(ver.23.0)對測序數(shù)據(jù)進(jìn)行病毒全基因組序列組裝拼接。

    1.4 全基因組序列分析 使用SARS-CoV-2分型系統(tǒng)Pangolin(https://pangolin.cog-uk.io/)對測定的病毒全基因序列進(jìn)行分型[13-14];通過在線工具Nextclade(https://clades.nextstrain.org/)對序列拼接結(jié)果進(jìn)行質(zhì)量檢查;從GISAID(https://www.gisaid.org)上篩選和下載同源序列以及不同Omicron亞分支變異株的代表序列作為參比序列,采用MAFFT軟件(ver.7.0)(https://mafft.cbrc.jp/alignment/server/)進(jìn)行多序列比對,使用MEGA軟件(ver.6.0)的NJ法(Bootstrap=1000)構(gòu)建系統(tǒng)進(jìn)化樹;使用Interactive Tree Of Life (https://itol.embl.de/)對進(jìn)化樹進(jìn)行注釋;使用Nextclade的SARS-CoV-2變異分析模塊對病毒全基因組序列進(jìn)行核苷酸與氨基酸突變分析、ACE2受體親和力及免疫逃逸能力預(yù)測[15]。

    2 結(jié) 果

    2.1 核酸復(fù)核結(jié)果 Real-time RT-PCR檢測病例A和B的鼻咽拭子核酸提取物,結(jié)果均為SARS-CoV-2核酸陽性,復(fù)核后的Ct值(ORF1ab/N)分別為:27.51/25.41、19.75/18.86,均滿足COVID-19診療方案(第九版)[16]要求的病毒全基因組測序標(biāo)準(zhǔn)。

    2.2 全基因組測序 采用三代測序技術(shù)從病例A的鼻咽拭子標(biāo)本中測得1株SARS-CoV-2全基因組序列,序列長度為28 910 bp,基因組覆蓋度為89.8%(<99.0%),平均測序深度為1 129×,Pangolin分型結(jié)果顯示病例A感染病毒屬于Omicron變異株BQ.1亞分支,但質(zhì)控分析顯示三代測序所獲序列存在基因組覆蓋度不足、因重復(fù)性堿基測序缺失引起的氨基酸翻譯錯(cuò)誤(移碼)等問題(詳見表1)。采用二代測序技術(shù)從兩病例標(biāo)本中測得2株SARS-CoV-2全基因組序列,序列長度分別為29 665 bp和29 682 bp,基因組覆蓋度分別99.3%和99.6%(均>99.0%),平均測序深度分別為2 973×和3 695×,Pangolin分型結(jié)果均為Omicron BQ.1變異株,與三代測序病毒分型結(jié)果一致;質(zhì)控分析顯示二代測序所獲2株SARS-CoV-2全基因組序列的完整度、基因組覆蓋度、測序深度、刺突蛋白(S蛋白)完整性等均較好,氨基酸翻譯正確無移碼,表明二代測序獲得的2株病毒序列均屬于高質(zhì)量序列,詳見表1。

    表1 SARS-CoV-2全基因組測序質(zhì)量分析結(jié)果

    2.3 全基因組進(jìn)化分析 根據(jù)病毒的分型結(jié)果,從GISAID數(shù)據(jù)庫中篩選并下載同源序列及不同Omicron亞分支的代表性序列作為參比序列,包括:BQ.1、BA.3、BA.4、BA.5、BA.5.1、BA.5.2、BA.5.3、BA.5.3.1、BF.7、XBB等10種Omicron亞分支。以Wuhan-Hu-1(NC_045512.2)作為“樹根”(Root),將篩選出的68條參比序列與2株SARS-CoV-2基因組序列構(gòu)建系統(tǒng)進(jìn)化樹。結(jié)果顯示,2株SARS-CoV-2基因組序列(2022XG-2599、2022XG-2600)共同處于Omicron變異株BQ.1進(jìn)化亞分支上,且位于同一進(jìn)化簇,表明2株SARS-CoV-2高度同源,見圖1。

    注:●表示武漢參考株序列Wuhan-Hu-1(NC_045512.2);▲表示本研究2株SARS-CoV-2全基因組序列;10種不同的Omicron亞分支序列用不同背景色顯示,亞分支名稱標(biāo)注在環(huán)形進(jìn)化樹外側(cè)。

    2.4 全基因組突變分析 用二代測序結(jié)果進(jìn)行SARS-CoV-2全基因組突變分析。結(jié)果顯示,同Wuhan-Hu-1相比,2022XG-2599(病例A)全基因組序列共有77個(gè)核苷酸突變位點(diǎn)(不計(jì)缺失/插入突變),其中1個(gè)突變位點(diǎn)位于5′UTR區(qū)域,76個(gè)突變位點(diǎn)位于蛋白編碼區(qū)。蛋白編碼區(qū)核苷酸突變涉及8種病毒蛋白編碼序列,包括2種非結(jié)構(gòu)蛋白、4種結(jié)構(gòu)蛋白和2種輔助蛋白。蛋白編碼區(qū)核苷酸突變導(dǎo)致57個(gè)錯(cuò)義突變和19個(gè)同義突變,主要突變發(fā)生在S蛋白(32個(gè))、ORF1a蛋白(18個(gè))、ORF1b蛋白(10個(gè))和 N蛋白(7個(gè))編碼區(qū)。在2022XG-2599全基因組氨基酸突變中,除S蛋白區(qū)A27S突變和第24-26位氨基酸缺失外,其余突變均屬于Omicron變異株BQ.1亞分支的特征性氨基酸突變(75%以上該型別變異株共有的突變位點(diǎn));同時(shí)S蛋白中的氨基酸突變中有19個(gè)位于受體結(jié)合域(receptor binding domain,RBD)(詳見表2)。2022XG-2600(病例B)序列與2022XG-2599雖然在全基因組長度及覆蓋度上有微小差異(序列長度相差17 bp,覆蓋度相0.3%),但二者共享全部核苷酸與氨基酸突變位點(diǎn)(詳見表2)。

    2.5 受體親和力及免疫逃逸能力預(yù)測 通過Nextclade在線預(yù)測本研究2株SARS-CoV-2的RBD氨基酸突變對人體細(xì)胞表面血管緊張素轉(zhuǎn)換酶2(angiotensinconverting enzyme 2,ACE2)受體的親和力及逃避中和抗體識別能力(免疫逃逸)的影響。結(jié)果顯示,以O(shè)micron變異株BA.2亞型為參照毒株,2株SARS-CoV-2的ACE2受體親和力評估值均為0.66,即ACE2受體的親和力比BA.2亞型高出約4.6倍(10∧0.66≈4.6);以O(shè)micron變異株BA.2亞型為參照毒株, 2株SARS-CoV-2的中和抗體逃逸能力評估值均為0.63,即只保留了約53.3%(e∧0.63≈0.533)BA.2亞型的抗體中和活性。預(yù)測結(jié)果表明本研究2株SARS-CoV-2較Omicron變異株BA.2亞型可能在ACE2受體親和力和免疫逃逸能力方面均有明顯提升。

    3 討 論

    本研究聯(lián)合采用Illumina Miseq(二代測序)和Nanopore MinION(三代測序)兩種高通量測序技術(shù),成功自2例COVID-19病例鼻咽拭子標(biāo)本中獲得SARS-CoV-2全基因組序列,分型結(jié)果顯示2株病毒均屬于Omicron BQ.1變異株,為福建省內(nèi)首次鑒定到的BQ.1病例;同時(shí)結(jié)合流行病學(xué)資料,推測兩例病例可能自境外輸入,并暴露于同一傳染源。本研究采用三代測序?qū)?biāo)本進(jìn)行測序,從核酸提取到獲得病毒全基因組序列僅用7 h,能夠初步鑒定出病毒型別屬于Omicron變異株BQ.1進(jìn)化分支,但僅成功測得1例病例病毒序列,且基因組覆蓋度低于99.0%,不能準(zhǔn)確分析病毒全基因組突變情況;采用二代測序則成功獲得兩例病例的病毒基因組序列,且基因組覆蓋度和測序深度均優(yōu)于三代測序,分型結(jié)果與三代測序一致,為后續(xù)基于突變的病毒精準(zhǔn)溯源提供可靠的數(shù)據(jù)基礎(chǔ),但耗時(shí)相對較長(20~22 h)。盡管Nanopore三代測序偶爾會出現(xiàn)連續(xù)堿基區(qū)域測序缺失的問題,使其測序準(zhǔn)確度受影響,但由于其耗時(shí)相對較短(6~8 h),因此能夠在疫情暴發(fā)初期迅速鑒別新型變異株,有利于第一時(shí)間發(fā)現(xiàn)感染源[17]。Illumina二代測序雖然耗時(shí)較長,但其測序錯(cuò)誤率相對較低,能夠更準(zhǔn)確分析病毒全基因組的突變情況,有利于精準(zhǔn)分析病毒傳播鏈和研判疫情傳播態(tài)勢。本研究結(jié)合不同測序技術(shù)的優(yōu)點(diǎn),將三代和二代測序技術(shù)聯(lián)用,既滿足了測序的時(shí)效性又解決了測序的準(zhǔn)確性問題,這種測序模式能夠從不同應(yīng)用層面更好地為疫情風(fēng)險(xiǎn)評估、流行病學(xué)監(jiān)測及公共衛(wèi)生決策提供技術(shù)支持。

    2022年7月,尼日利亞首次報(bào)道了Omicron變異株BQ.1亞型感染病例,之后BQ.1便在歐洲和北美迅速擴(kuò)散,BQ.1是Omicron變異株BA.5分支在人群流行過程中變異產(chǎn)生的亞分支,該變異株在S蛋白的一些關(guān)鍵抗原位點(diǎn)上發(fā)生突變,這些突變賦予其免疫逃逸優(yōu)勢,使其具有較高的再感染和疫苗突破潛力,可能會導(dǎo)致全球范圍內(nèi)的感染和死亡病例數(shù)量激增,將對公共衛(wèi)生安全構(gòu)成新的威脅[18-19]。WHO將BQ.1變異株定義為VUMs(Currently circulating variants under monitoring)[20],WHO的統(tǒng)計(jì)數(shù)據(jù)顯示[21],2023年第3~7周,BQ.1系列變異株在所有VUMs中的周流行率為20.03%~36.40%,已累計(jì)在133個(gè)國家檢測到。BQ.1變異株的S蛋白在BA.5變異株基礎(chǔ)上新增K444T和N460K特征性突變位點(diǎn),本研究2株BQ.1變異株在RBD區(qū)均含有這2個(gè)突變位點(diǎn),有研究顯示[22],K444T和N460K的組合突變拓展了氫鍵相互作用的網(wǎng)絡(luò),可能有利于對中和抗體的抵抗,從而增強(qiáng)了其免疫逃避能力。RBD是介導(dǎo)S蛋白與ACE2受體結(jié)合和誘導(dǎo)宿主免疫系統(tǒng)產(chǎn)生中和抗體的重要結(jié)構(gòu)域,該結(jié)構(gòu)域突變可能會導(dǎo)致病毒免疫逃逸能力增強(qiáng)和影響疫苗免疫效果[23-24]。本研究2株BQ.1變異株在RBD區(qū)包含19個(gè)氨基酸突變,其中K417N、L452R、S477N、N501Y等是WHO定義的MOI(Mutation of Interest)位點(diǎn),對其人群適應(yīng)性具有重要作用[25-26]。相關(guān)研究發(fā)現(xiàn)[27],Omicron變異株可能通過RBD區(qū)氨基酸殘基組合的趨同突變來逃避抗體的中和作用,同時(shí)又保持足夠的ACE2結(jié)合能力,如N460K、F486V與N460K組合、L452R與F486S/V組合、K444N/T與F486V組合等。本研究2株BQ.1變異株在RBD區(qū)出現(xiàn)了K444T、L452R、N460K、F486V等熱點(diǎn)突變,這些突變的不同組合可能對變異株的免疫逃避能力和ACE2受體親和力產(chǎn)生趨同作用。這與在線軟件預(yù)測的2株BQ.1變異株的ACE2受體親和力及免疫逃逸能力結(jié)果一致,表明變異株的生物學(xué)特性可能發(fā)生較明顯變化,需給予持續(xù)監(jiān)測。當(dāng)前SARS-CoV-2疫苗保護(hù)效力及BQ.1變異株可能導(dǎo)致突破性感染和再感染風(fēng)險(xiǎn)還需要進(jìn)一步研究。

    當(dāng)前我國疫情防控已進(jìn)入新階段,隨著入境政策的優(yōu)化調(diào)整,更多SARS-CoV-2新型變異株將持續(xù)輸入國內(nèi)并引起本土流行傳播,通過對具有代表性的新冠病毒感染者樣本進(jìn)行病毒全基因組測序和突變分析,可以動態(tài)監(jiān)測本土流行的變異株組成和分布,及時(shí)發(fā)現(xiàn)國內(nèi)外重點(diǎn)關(guān)注變異株,評估變異株的傳播力、致病力和免疫逃逸能力等變化和風(fēng)險(xiǎn)研判,為本土新增SARS-CoV-2變異株的大規(guī)模傳播提供早期預(yù)警。本研究采用二代和三代測序技術(shù)聯(lián)用進(jìn)行病毒基因組測序,兼顧了準(zhǔn)確性和時(shí)效性,可為SARS-CoV-2變異監(jiān)測工作及制定疫情防控策略提供技術(shù)參考與支持。

    利益沖突:無

    猜你喜歡
    變異基因組位點(diǎn)
    鎳基單晶高溫合金多組元置換的第一性原理研究
    上海金屬(2021年6期)2021-12-02 10:47:20
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    CLOCK基因rs4580704多態(tài)性位點(diǎn)與2型糖尿病和睡眠質(zhì)量的相關(guān)性
    變異危機(jī)
    變異
    二項(xiàng)式通項(xiàng)公式在遺傳學(xué)計(jì)算中的運(yùn)用*
    變異的蚊子
    百科知識(2015年18期)2015-09-10 07:22:44
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    含內(nèi)含子的核糖體蛋白基因轉(zhuǎn)錄起始位點(diǎn)情況分析
    菩萨蛮人人尽说江南好唐韦庄| 亚洲国产色片| 久久久精品国产亚洲av高清涩受| 一级毛片黄色毛片免费观看视频| 美女午夜性视频免费| 国产黄色视频一区二区在线观看| 精品一区在线观看国产| 国产 一区精品| 久久狼人影院| 亚洲四区av| 免费观看性生交大片5| 国产成人精品婷婷| 久久精品久久久久久噜噜老黄| 永久网站在线| 成年女人在线观看亚洲视频| 久久久久精品性色| 国产综合精华液| 日日撸夜夜添| 日韩av不卡免费在线播放| 日韩熟女老妇一区二区性免费视频| 午夜福利视频在线观看免费| 少妇人妻精品综合一区二区| 免费观看性生交大片5| 国产精品亚洲av一区麻豆 | 成人亚洲欧美一区二区av| 成年女人在线观看亚洲视频| 九九爱精品视频在线观看| 女人被躁到高潮嗷嗷叫费观| 老司机亚洲免费影院| 亚洲精品av麻豆狂野| 亚洲内射少妇av| 久久精品国产自在天天线| 美女大奶头黄色视频| 人人妻人人澡人人爽人人夜夜| 亚洲第一av免费看| 在线亚洲精品国产二区图片欧美| 日韩免费高清中文字幕av| 最近中文字幕高清免费大全6| 亚洲综合色网址| 精品国产一区二区久久| freevideosex欧美| 熟女电影av网| 亚洲在久久综合| 婷婷色av中文字幕| 熟女少妇亚洲综合色aaa.| 美女脱内裤让男人舔精品视频| 午夜福利视频在线观看免费| 成年av动漫网址| 免费看av在线观看网站| 男人爽女人下面视频在线观看| 国产在视频线精品| 久久热在线av| √禁漫天堂资源中文www| 可以免费在线观看a视频的电影网站 | 一二三四中文在线观看免费高清| 国产精品 欧美亚洲| 日韩一卡2卡3卡4卡2021年| 亚洲熟女精品中文字幕| 成人毛片60女人毛片免费| 国产白丝娇喘喷水9色精品| 91精品三级在线观看| 日韩一区二区视频免费看| 亚洲伊人色综图| 女人被躁到高潮嗷嗷叫费观| 日韩一区二区三区影片| 亚洲一区中文字幕在线| 日韩免费高清中文字幕av| 激情五月婷婷亚洲| 日日撸夜夜添| 最新中文字幕久久久久| 日韩av在线免费看完整版不卡| 亚洲精品aⅴ在线观看| 精品国产乱码久久久久久男人| 人人妻人人爽人人添夜夜欢视频| 午夜福利一区二区在线看| 老熟女久久久| 欧美成人午夜精品| av不卡在线播放| 亚洲国产欧美网| 久久青草综合色| 国产精品欧美亚洲77777| 亚洲一级一片aⅴ在线观看| 男人操女人黄网站| 亚洲精品日韩在线中文字幕| 日韩av不卡免费在线播放| 寂寞人妻少妇视频99o| 久久国产精品大桥未久av| 90打野战视频偷拍视频| 捣出白浆h1v1| 考比视频在线观看| 考比视频在线观看| 日韩一区二区视频免费看| 男人添女人高潮全过程视频| 精品人妻一区二区三区麻豆| 中文字幕亚洲精品专区| 一区二区三区激情视频| 少妇 在线观看| 啦啦啦啦在线视频资源| 人妻系列 视频| 97精品久久久久久久久久精品| 亚洲一区中文字幕在线| 丝袜脚勾引网站| 黑人巨大精品欧美一区二区蜜桃| 成人午夜精彩视频在线观看| 亚洲第一av免费看| 岛国毛片在线播放| 高清黄色对白视频在线免费看| 国产精品成人在线| 日韩中文字幕视频在线看片| 热re99久久国产66热| 欧美日韩精品网址| 亚洲美女搞黄在线观看| 赤兔流量卡办理| 黑丝袜美女国产一区| 亚洲一区二区三区欧美精品| 香蕉国产在线看| 最近最新中文字幕免费大全7| 看免费av毛片| 午夜福利影视在线免费观看| 久热久热在线精品观看| 欧美精品高潮呻吟av久久| 性高湖久久久久久久久免费观看| 久久精品国产亚洲av高清一级| 满18在线观看网站| 久久久精品国产亚洲av高清涩受| 国产高清国产精品国产三级| 国产精品熟女久久久久浪| 18禁动态无遮挡网站| 在线观看人妻少妇| 久久久久国产网址| 交换朋友夫妻互换小说| av在线观看视频网站免费| 亚洲av福利一区| 高清不卡的av网站| 亚洲av免费高清在线观看| 成年av动漫网址| 国产av一区二区精品久久| 国产成人精品在线电影| 久久久亚洲精品成人影院| 中国三级夫妇交换| 亚洲 欧美一区二区三区| 国产精品久久久久久精品电影小说| 精品99又大又爽又粗少妇毛片| 午夜免费男女啪啪视频观看| 国产熟女欧美一区二区| 成人手机av| 亚洲av.av天堂| 久久热在线av| 免费日韩欧美在线观看| 美女视频免费永久观看网站| 在线观看人妻少妇| 日韩,欧美,国产一区二区三区| 边亲边吃奶的免费视频| 亚洲伊人久久精品综合| 黄片播放在线免费| 王馨瑶露胸无遮挡在线观看| 美女视频免费永久观看网站| 欧美激情 高清一区二区三区| 国产精品一区二区在线观看99| 18在线观看网站| 久久久久视频综合| 国产成人免费无遮挡视频| 国产精品久久久久成人av| 久久久久久久国产电影| 国产成人午夜福利电影在线观看| 国产免费福利视频在线观看| 成人国产av品久久久| 久久狼人影院| 夫妻性生交免费视频一级片| 观看av在线不卡| 两个人看的免费小视频| 热99国产精品久久久久久7| 波多野结衣av一区二区av| 久久久久久久久久久免费av| 男人操女人黄网站| av天堂久久9| 国产老妇伦熟女老妇高清| 久久久久久免费高清国产稀缺| 久久精品熟女亚洲av麻豆精品| 成人亚洲精品一区在线观看| 在线观看免费日韩欧美大片| 丰满乱子伦码专区| 日韩 亚洲 欧美在线| 国产乱人偷精品视频| 天天影视国产精品| 精品一品国产午夜福利视频| 美女高潮到喷水免费观看| 日本-黄色视频高清免费观看| 在线观看免费日韩欧美大片| 最近最新中文字幕免费大全7| 国产精品 欧美亚洲| 啦啦啦中文免费视频观看日本| 高清欧美精品videossex| 亚洲三区欧美一区| 国产日韩欧美在线精品| 日韩,欧美,国产一区二区三区| 亚洲av中文av极速乱| 免费播放大片免费观看视频在线观看| 久久久久国产一级毛片高清牌| 天堂8中文在线网| 亚洲成人手机| 最近2019中文字幕mv第一页| 黑人巨大精品欧美一区二区蜜桃| 成年动漫av网址| 91精品三级在线观看| 叶爱在线成人免费视频播放| 在线观看免费日韩欧美大片| 午夜久久久在线观看| 美女福利国产在线| 不卡av一区二区三区| 不卡av一区二区三区| 啦啦啦啦在线视频资源| 男女午夜视频在线观看| 热99久久久久精品小说推荐| 亚洲成人av在线免费| 大话2 男鬼变身卡| 亚洲色图 男人天堂 中文字幕| 大片电影免费在线观看免费| 精品国产国语对白av| 成人二区视频| 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看| 成人免费观看视频高清| 黄片无遮挡物在线观看| 2018国产大陆天天弄谢| 91精品伊人久久大香线蕉| 国产 精品1| 黑人巨大精品欧美一区二区蜜桃| 久久久精品免费免费高清| 免费女性裸体啪啪无遮挡网站| 青草久久国产| 日韩av免费高清视频| 久久97久久精品| 中文字幕制服av| 国产精品.久久久| 国产精品一国产av| 美女国产视频在线观看| 久久久国产一区二区| 久久热在线av| 卡戴珊不雅视频在线播放| 夜夜骑夜夜射夜夜干| 午夜91福利影院| 中文乱码字字幕精品一区二区三区| 在线免费观看不下载黄p国产| 一区二区三区激情视频| 男人舔女人的私密视频| 欧美av亚洲av综合av国产av | 亚洲五月色婷婷综合| 国产无遮挡羞羞视频在线观看| 黄色 视频免费看| 自线自在国产av| 国产深夜福利视频在线观看| 国产免费又黄又爽又色| 午夜免费男女啪啪视频观看| 在线天堂中文资源库| 久久人妻熟女aⅴ| 这个男人来自地球电影免费观看 | 国产白丝娇喘喷水9色精品| 一级毛片黄色毛片免费观看视频| 日韩av在线免费看完整版不卡| 69精品国产乱码久久久| 爱豆传媒免费全集在线观看| 国产精品三级大全| 天美传媒精品一区二区| 国产免费一区二区三区四区乱码| 久久久久网色| 免费看av在线观看网站| 在线观看人妻少妇| 欧美老熟妇乱子伦牲交| 九草在线视频观看| 永久网站在线| 成人毛片60女人毛片免费| 91国产中文字幕| 亚洲av.av天堂| 亚洲精品国产色婷婷电影| 美女福利国产在线| 国产 精品1| 如何舔出高潮| av在线老鸭窝| 日韩在线高清观看一区二区三区| 亚洲人成77777在线视频| 亚洲av欧美aⅴ国产| 看非洲黑人一级黄片| 大片电影免费在线观看免费| 亚洲欧洲精品一区二区精品久久久 | 9热在线视频观看99| 国产乱来视频区| 精品亚洲成a人片在线观看| 成人午夜精彩视频在线观看| 深夜精品福利| 成年美女黄网站色视频大全免费| 国产一区二区三区av在线| 亚洲av福利一区| 91在线精品国自产拍蜜月| 2021少妇久久久久久久久久久| 激情视频va一区二区三区| 精品亚洲成国产av| 人妻一区二区av| 肉色欧美久久久久久久蜜桃| 青草久久国产| 国产男人的电影天堂91| 亚洲,一卡二卡三卡| 日韩 亚洲 欧美在线| 久久精品国产a三级三级三级| av视频免费观看在线观看| 精品少妇黑人巨大在线播放| 亚洲一码二码三码区别大吗| 午夜av观看不卡| 成年女人毛片免费观看观看9 | 久久精品人人爽人人爽视色| 天天操日日干夜夜撸| 日本vs欧美在线观看视频| 久久久久人妻精品一区果冻| 咕卡用的链子| 亚洲欧美日韩另类电影网站| 午夜精品国产一区二区电影| 日韩电影二区| 成人免费观看视频高清| 超色免费av| 国产男人的电影天堂91| 日日摸夜夜添夜夜爱| 国产成人欧美| 大陆偷拍与自拍| 成年女人毛片免费观看观看9 | 日韩大片免费观看网站| 性色avwww在线观看| 国产日韩欧美在线精品| 制服丝袜香蕉在线| 国产av国产精品国产| 日本欧美视频一区| 永久免费av网站大全| √禁漫天堂资源中文www| 青春草国产在线视频| 亚洲精品美女久久av网站| 亚洲国产av影院在线观看| a级毛片黄视频| 日韩在线高清观看一区二区三区| 精品国产一区二区久久| 久久国内精品自在自线图片| 久久婷婷青草| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 99精国产麻豆久久婷婷| 人人妻人人添人人爽欧美一区卜| 亚洲欧美色中文字幕在线| 精品国产乱码久久久久久小说| av女优亚洲男人天堂| 免费大片黄手机在线观看| 日本vs欧美在线观看视频| 一级黄片播放器| 高清黄色对白视频在线免费看| 9色porny在线观看| 精品亚洲成国产av| 在线观看免费日韩欧美大片| 精品卡一卡二卡四卡免费| 免费在线观看完整版高清| 中文精品一卡2卡3卡4更新| 国产日韩欧美视频二区| 久久这里有精品视频免费| 母亲3免费完整高清在线观看 | a级毛片黄视频| 欧美老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 国产免费视频播放在线视频| 男人添女人高潮全过程视频| 一区在线观看完整版| 精品卡一卡二卡四卡免费| 日本欧美视频一区| 精品福利永久在线观看| 亚洲情色 制服丝袜| 美女中出高潮动态图| 26uuu在线亚洲综合色| 日韩 亚洲 欧美在线| 老司机影院毛片| 天美传媒精品一区二区| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美成人精品一区二区| 日本午夜av视频| 久久精品久久精品一区二区三区| 国产精品免费视频内射| 亚洲人成网站在线观看播放| 中文字幕人妻熟女乱码| 国产日韩欧美视频二区| 久久国产精品大桥未久av| 免费久久久久久久精品成人欧美视频| 下体分泌物呈黄色| 亚洲精品美女久久av网站| 丰满乱子伦码专区| 欧美97在线视频| 十分钟在线观看高清视频www| 久久久久久人人人人人| 1024视频免费在线观看| 不卡av一区二区三区| 国产在线视频一区二区| 在线观看国产h片| 伊人久久国产一区二区| 亚洲欧美一区二区三区黑人 | 免费黄频网站在线观看国产| 制服诱惑二区| 精品少妇一区二区三区视频日本电影 | 制服人妻中文乱码| 天天影视国产精品| 2022亚洲国产成人精品| 热99国产精品久久久久久7| 国产精品不卡视频一区二区| 少妇人妻 视频| 久久久久久免费高清国产稀缺| 日韩人妻精品一区2区三区| 一区二区三区乱码不卡18| 欧美黄色片欧美黄色片| 亚洲欧洲国产日韩| 最近的中文字幕免费完整| 18+在线观看网站| 欧美人与性动交α欧美精品济南到 | 自拍欧美九色日韩亚洲蝌蚪91| 一级黄片播放器| 国产免费现黄频在线看| 国产精品秋霞免费鲁丝片| 99久国产av精品国产电影| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件| 啦啦啦中文免费视频观看日本| tube8黄色片| 欧美亚洲 丝袜 人妻 在线| 亚洲色图综合在线观看| 少妇精品久久久久久久| 中文天堂在线官网| 久久国产亚洲av麻豆专区| 久久国产精品男人的天堂亚洲| 亚洲人成电影观看| 久久青草综合色| 国产日韩欧美在线精品| 欧美精品一区二区免费开放| 丝袜美足系列| 一个人免费看片子| 亚洲av中文av极速乱| 在线观看免费高清a一片| 国产一区有黄有色的免费视频| 免费人妻精品一区二区三区视频| 丰满少妇做爰视频| 两个人免费观看高清视频| 免费播放大片免费观看视频在线观看| 夜夜骑夜夜射夜夜干| 国产不卡av网站在线观看| 日本vs欧美在线观看视频| 侵犯人妻中文字幕一二三四区| 亚洲国产精品成人久久小说| 国产一级毛片在线| 亚洲成人av在线免费| 成人国产av品久久久| 国产毛片在线视频| 亚洲一级一片aⅴ在线观看| 建设人人有责人人尽责人人享有的| 伦理电影大哥的女人| av网站在线播放免费| 自线自在国产av| 亚洲精品视频女| 咕卡用的链子| 熟女av电影| 大码成人一级视频| 亚洲av在线观看美女高潮| 老鸭窝网址在线观看| 好男人视频免费观看在线| 久久av网站| av又黄又爽大尺度在线免费看| 蜜桃国产av成人99| 欧美精品一区二区大全| 亚洲精品av麻豆狂野| 久热这里只有精品99| 色吧在线观看| 美女视频免费永久观看网站| 亚洲一码二码三码区别大吗| xxxhd国产人妻xxx| 美女中出高潮动态图| 日本wwww免费看| 建设人人有责人人尽责人人享有的| 99精国产麻豆久久婷婷| 人人澡人人妻人| 欧美97在线视频| 菩萨蛮人人尽说江南好唐韦庄| 一级a爱视频在线免费观看| 老司机影院毛片| 日本爱情动作片www.在线观看| 一级爰片在线观看| 欧美日韩精品成人综合77777| 亚洲欧美精品综合一区二区三区 | 国产精品偷伦视频观看了| 热re99久久国产66热| kizo精华| 国产精品人妻久久久影院| 纯流量卡能插随身wifi吗| 丝袜人妻中文字幕| 性色avwww在线观看| 欧美激情极品国产一区二区三区| 色哟哟·www| 国产亚洲av片在线观看秒播厂| 日本色播在线视频| 日韩 亚洲 欧美在线| 丝袜美腿诱惑在线| 青青草视频在线视频观看| 男女无遮挡免费网站观看| 亚洲美女视频黄频| 天堂俺去俺来也www色官网| 高清黄色对白视频在线免费看| 啦啦啦在线免费观看视频4| 蜜桃在线观看..| 久久久久国产一级毛片高清牌| 午夜精品国产一区二区电影| 国产黄色免费在线视频| 在线 av 中文字幕| 黑丝袜美女国产一区| 校园人妻丝袜中文字幕| 国产色婷婷99| 日本wwww免费看| 飞空精品影院首页| 日韩一本色道免费dvd| 欧美人与性动交α欧美软件| 亚洲中文av在线| 成人漫画全彩无遮挡| 日韩电影二区| 七月丁香在线播放| 午夜精品国产一区二区电影| 日产精品乱码卡一卡2卡三| 在现免费观看毛片| 亚洲成人av在线免费| 毛片一级片免费看久久久久| 深夜精品福利| 男女国产视频网站| 亚洲三级黄色毛片| 国产精品亚洲av一区麻豆 | 啦啦啦啦在线视频资源| 日韩一本色道免费dvd| 国产精品熟女久久久久浪| 欧美日韩综合久久久久久| 黄色怎么调成土黄色| 午夜福利,免费看| 男女国产视频网站| 精品一区二区三卡| 啦啦啦在线观看免费高清www| 精品亚洲成国产av| 少妇人妻 视频| 精品99又大又爽又粗少妇毛片| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 美女大奶头黄色视频| 久久国内精品自在自线图片| 人人妻人人澡人人爽人人夜夜| 国产日韩一区二区三区精品不卡| 亚洲国产看品久久| 日韩三级伦理在线观看| 新久久久久国产一级毛片| 成人国产麻豆网| 亚洲国产成人一精品久久久| 免费高清在线观看视频在线观看| av又黄又爽大尺度在线免费看| 久久精品国产亚洲av天美| av福利片在线| 香蕉精品网在线| 亚洲av国产av综合av卡| 人妻少妇偷人精品九色| 免费女性裸体啪啪无遮挡网站| av国产久精品久网站免费入址| 久久久国产精品麻豆| 久久久久久久精品精品| 国产伦理片在线播放av一区| www.熟女人妻精品国产| 一区二区日韩欧美中文字幕| 国产男人的电影天堂91| 一级爰片在线观看| 可以免费在线观看a视频的电影网站 | 丰满乱子伦码专区| 成人午夜精彩视频在线观看| av一本久久久久| 国产不卡av网站在线观看| 久久久精品94久久精品| av视频免费观看在线观看| 97在线视频观看| 欧美变态另类bdsm刘玥| 亚洲精品国产色婷婷电影| 少妇被粗大的猛进出69影院| 午夜久久久在线观看| 男人添女人高潮全过程视频| 亚洲综合色网址| 一区福利在线观看| 男女高潮啪啪啪动态图| 亚洲精华国产精华液的使用体验| 男女国产视频网站| a级片在线免费高清观看视频| 搡老乐熟女国产| 黄片小视频在线播放| 欧美日韩视频高清一区二区三区二| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91精品国产国语对白视频| 国产精品.久久久| 丝袜美足系列| 一级毛片黄色毛片免费观看视频| 国产在视频线精品| 黄频高清免费视频| 丝袜脚勾引网站| 国产黄色免费在线视频| 免费黄频网站在线观看国产| 日本wwww免费看| 丝瓜视频免费看黄片| 亚洲国产av影院在线观看| 涩涩av久久男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 日日啪夜夜爽| 日韩三级伦理在线观看| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 免费看不卡的av| 久久久久久久久免费视频了| 亚洲久久久国产精品| 高清不卡的av网站| 国产熟女欧美一区二区|