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

    漢灘病毒糖蛋白Gn的T細(xì)胞表位預(yù)測(cè)及其跨種屬驗(yàn)證

    2023-10-18 12:56:38孫報(bào)增秦聰聰張俊琦王永凱張文標(biāo)劉瑞波白天原張志輝張宇絲楊琨姜東伯
    解放軍醫(yī)學(xué)雜志 2023年9期
    關(guān)鍵詞:免疫原性表位親和力

    孫報(bào)增,秦聰聰,張俊琦,王永凱,張文標(biāo),劉瑞波,白天原,張志輝,張宇絲,楊琨,2,3,姜東伯,2,4*

    1空軍軍醫(yī)大學(xué)基礎(chǔ)醫(yī)學(xué)院免疫學(xué)教研室,陜西 西安 710032;2空軍軍醫(yī)大學(xué)基礎(chǔ)醫(yī)學(xué)院生物武器損傷防治藥物重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710032;3空軍醫(yī)科大學(xué)唐都醫(yī)院風(fēng)濕科,陜西 西安 710032;4空軍軍醫(yī)大學(xué)基礎(chǔ)醫(yī)學(xué)院微生物學(xué)教研室,陜西 西安 710032

    漢坦病毒屬(Hantavirus)屬于布尼亞病毒科,包括漢灘病毒(Hantaan virus,HTNV)、多布拉-貝爾格萊德病毒、普馬拉病毒和首爾病毒等,其在歐亞大陸主要引起腎綜合征出血熱(hemorrhagic fever with renal syndrome,HFRS)[1-2]。HFRS 以鼠類等嚙齒類動(dòng)物為主要傳染源,表現(xiàn)為發(fā)熱、出血、急性腎功能損傷等。我國是HTNV 所致HFRS 的主要疫源地,占全球感染病例的70%~90%,發(fā)病率和病死率均較高[3]。HTNV 為三股單負(fù)鏈RNA 病毒,由S、M 和L片段組成,分別編碼RNA 依賴性RNA 聚合酶(RdRp)、糖蛋白(GP)及核衣殼蛋白(NP)[4]。其中GP構(gòu)成HTNV 的主要抗原,同時(shí)可被加工分解成包膜糖蛋白-N 末端(Gn)及包膜糖蛋白-C 末端(Gc),構(gòu)成病毒吸附蛋白,能夠識(shí)別并結(jié)合細(xì)胞表面受體,刺激機(jī)體產(chǎn)生中和抗體,在病毒對(duì)機(jī)體細(xì)胞的侵襲及感染過程中發(fā)揮重要作用,因此可作為HFRS 疫苗研究的潛在靶點(diǎn)[5]。然而,在現(xiàn)有的表位研究中,篩選鑒定的限制性表位不全面,針對(duì)表位的研究比較單一,不能全面揭示表位特性。同時(shí),尋找人類的表位較困難,難以開展相關(guān)的實(shí)驗(yàn),而前期研究發(fā)現(xiàn),采用計(jì)算機(jī)預(yù)測(cè)的方法具有全面準(zhǔn)確、經(jīng)濟(jì)效益高的優(yōu)勢(shì)[6-8]。因此,本研究對(duì)Gn 進(jìn)行了系統(tǒng)分析,以探索主要組織相容性復(fù)合體(MHC)-Ⅰ、Ⅱ類表位免疫反應(yīng)性的相關(guān)規(guī)律,指導(dǎo)新型HFRS 表位疫苗的研發(fā),預(yù)防漢坦病毒的跨物種傳播風(fēng)險(xiǎn)。

    1 材料與方法

    1.1 Gn 蛋白序列獲取 從UniProt 數(shù)據(jù)庫(https://www.uniprot.org/)檢索HTNV 76-118 株Gn 蛋白序列(PRO_0000036816)[9],氨基酸序列從aa1至aa649。

    1.2 Gn蛋白表位親和力預(yù)測(cè)

    1.2.1 預(yù)測(cè)工具的選取 利用5 種預(yù)測(cè)算法,包括IEDB(http://tools.iedb.org/mhci/)[10]、 SMMPMBEC(https://github.com/ykimbiology/smmpmbec)[11]、NetMHCpan 4.1(https://services.healthtech.dtu.dk/service.php?NetMHCpan-4.1)[12]、SYFPEITHI(http://www.syfpeithi.de/bin/MHCServer.dll/EpitopePrediction.htm)[13]、Rankpep(http://imed.med.ucm.es/Tools/rankpep.html)[14]全長(zhǎng)無偏倚地獲得高親和力的表位候選肽,并對(duì)HTNV Gc 蛋白序列進(jìn)行序列寡肽分割,計(jì)算MHC-Ⅰ類分子與9 肽表位之間的親和力;利用4 種預(yù)測(cè)算法,包括IEDB、NetMHC Ⅱpan3.2(https://services.healthtech.dtu.dk/service.php? NetMHCIIpan-3.2)[12]、SYFPEITHI、Rankpep來預(yù)測(cè)MHC-Ⅱ類分子與15肽表位之間的親和力。

    1.2.2 MHC-Ⅰ類限制性表位的預(yù)測(cè)與篩選 對(duì)IEDB推薦的人類白細(xì)胞抗原(HLA)-Ⅰ類亞型等位基因及Allel Frequency[15]網(wǎng)站(http://www.allelefrequencies.net/)中各地區(qū)基因頻率>10%的HLA基因型進(jìn)行親和力分析,包括HLA-A1(-A*0101、HLA-A*2601、-A*3001、-A*3002)、HLA-A2(HLA-A*0201、-A*0203、-A*0206、-A*0207、-A*0212、-A*6802)、HLA-A3(HLA-A*0301、 -A*1101、 -A*3001、 -A*3101、 -A*3301、 -A*6801)、 HLA-A24(HLA-A*2301、 -A*2402、 -A*2407、-A*3201)、HLA-A*3303、HLA-A*3108、HLA-A*3401、HLA-B7(HLA-B*0702、-B*3501、-B*5101、-B*5301)、HLA-B8(HLA-B*0801)、 HLA-B*1301、 HLA-B*1405、HLA-B15(HLA-B*1501、-B*1502、-B*1506、-B*1512、-B*1521)、 HLA-B*1807、 HLA-B*3508、 HLA-B*3903、HLA-B44(HLA-B*4001、-B*4002、-B*4006、-B*4402、-B*4403、 -B*4404)、 HLA-B46(HLA-B*4601)、 HLAB50(HLA-B*5001)、 HLA-B*5110、 HLA-B*5201、HLA-B*5502、 HLA-B*5601、 HLA-B58(HLA-B*5701、-B*5801)。上述基因型覆蓋了超過99.3% 人群的MHC-Ⅰ類亞型。以每種算法結(jié)果中前2%的預(yù)測(cè)表位為候選優(yōu)勢(shì)表位,至少在3 種預(yù)測(cè)算法結(jié)果中出現(xiàn)的表位作為優(yōu)勢(shì)表位[16]。

    1.2.3 MHC-Ⅱ類限制性表位的預(yù)測(cè)與篩選 MHC-Ⅱ類分子包括4 個(gè)HLA-Ⅱ超家族:DRB1(DRB1*0101、DRB1*0301、DRB1*0304、DRB1*0401、DRB1*0405、DRB1*0407、DRB1*0412、DRB1*0701、DRB1*0704、DRB1*0802、DRB1*0803、DRB1*0901、DRB1*1101、DRB1*1104、DRB1*1201、DRB1*1202、DRB1*1301、DRB1*1302、DRB1*1304、DRB1*1401、DRB1*1402、DRB1*1408、DRB1*1409、DRB1*1501、DRB1*1502、DRB1*1503、DRB1*1602),DRB3/4/5(DRB3*0101、 DRB3*0202、 DRB4*0101、 DRB5*0101),DQA1/DQB1(DQA1*0501/DQB1*0201、 DQA1*0501/DQB1*0301、DQA1*0301/DQB1*0302、DQA1*0401/DQB1*0402、DQA1*0101/DQB1*0501、DQA1*0102/DQB1*0602),DPB1(DPB1*0101、DPB1*0201、DPB1*0401、DPB1*0402、DPB1*0501、DPB1*1401)。上述基因型覆蓋了超過98%人群的MHC-Ⅱ類亞型,且在Allel Frequency 網(wǎng)站(http://www.allelefrequencies.net/)中各地區(qū)的基因型頻率>10%。以每種算法結(jié)果中前2%的預(yù)測(cè)表位為候選優(yōu)勢(shì)表位,至少在2種預(yù)測(cè)算法結(jié)果中出現(xiàn)的表位作為優(yōu)勢(shì)表位。

    1.3 HTNV Gn MHC-Ⅰ/Ⅱ類限制性表位的免疫原性分析 高親和力肽不能充分誘導(dǎo)免疫應(yīng)答,因此需對(duì)優(yōu)勢(shì)表位進(jìn)行免疫原性分析[17]。免疫原性由氨基酸序列所決定。本研究采用VaxiJen v2.0(http://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html)[18]計(jì)算9肽、15肽表位的免疫原性,以評(píng)分>0.5為陽性標(biāo)準(zhǔn),即該表位被視為強(qiáng)免疫原性,否則為不具有免疫原性。

    為進(jìn)一步探究高親和力、強(qiáng)免疫原性的優(yōu)勢(shì)表位跨人群、物種的交叉免疫反應(yīng)性,將各個(gè)基因型篩選出來的高親和力、強(qiáng)免疫原性的優(yōu)勢(shì)表位進(jìn)行對(duì)比分析,15 肽表位用其核心9 肽來表示,并采用Excel繪制三維柱狀圖。

    1.4 HTNV Gn MHC-Ⅰ/Ⅱ類限制性表位的保守性預(yù)測(cè) 為了確定HTNV 76-118毒株Gn序列中優(yōu)勢(shì)表位的進(jìn)化保守性,采用Blastp(https://blast.ncbi.nlm.nih.gov/Blast.cgi)[19]工具預(yù)測(cè)高親和力、強(qiáng)免疫原性的9肽和15肽表位的種間種內(nèi)保守性。其中,種間保守的評(píng)價(jià)標(biāo)準(zhǔn)為除漢坦病毒(HTNV,taxid:1980471)、人(Human,taxid:9606)、鼠(Mouse,taxid:10088)外的正漢坦病毒(Orthohantavirus,taxid:1980442),而種內(nèi)保守評(píng)價(jià)標(biāo)準(zhǔn)為除漢灘病毒76-118 株(Hantaan virus 76-118,taxid:11602)、人(Human,taxid:9606)、鼠(Mouse,taxid:10088)外的漢灘病毒(Hantaan hantavirus,taxid:1980471)。在分析結(jié)果中,若E值<10-5,則認(rèn)為保守。本實(shí)驗(yàn)選擇種間、種內(nèi)均保守的表位作為候選靶標(biāo)。

    1.5 雙向聚類分析 MHC 分子多態(tài)性和表位肽多樣性構(gòu)成了pMHC 相互作用的差異性。為了將其中的關(guān)系可視化,利用TBtools[20]對(duì)MHC-Ⅰ、Ⅱ基因型與HTNV Gn 表位親和力的排名進(jìn)行雙向?qū)哟尉垲?。親和力排序數(shù)據(jù)通過以2為底的對(duì)數(shù)和Z得分標(biāo)準(zhǔn)化處理之后,用Complete 方法進(jìn)行歐氏距離的雙向?qū)哟尉垲?。分析包?7 個(gè)MHC-Ⅰ類分子與HTNV Gn 9 肽的相互作用及51 個(gè)MHC-Ⅱ類分子與HTNV Gn 15肽表位的相互作用。

    1.6 模擬pMHC分子對(duì)接

    1.6.1 HTNV Gn 9肽與MHC-Ⅰ類分子模擬對(duì)接HPEPDOCK(http://huanglab.phys.hust.edu.cn/hpepdock/)[21]是一種新型的網(wǎng)絡(luò)服務(wù)器,通過輸入9 肽優(yōu)勢(shì)表位序列與MHC-Ⅰ類分子PDB 格式文件即可模擬pMHC對(duì)接并獲得對(duì)接模型。使用RCSB PDB數(shù)據(jù)庫(https://www.rcsb.org/)即可獲得相關(guān)的MHC-Ⅰ類分子3D 結(jié)構(gòu)數(shù)據(jù)[HLA-A0206(3OXR)、 HLA-A2402(5WWJ)、 HLA-A0101(4NQX)、 HLA-B5301(1A1N)、HLA-B4001(6IEX)、 HLA-B5801(5V5L)、 HLA-B0801(3X13)、H2Kb(6JQ3)、H2-Ld(6L9M)]。對(duì)高親和力、強(qiáng)免疫原性和種間種內(nèi)均保守的9 肽優(yōu)勢(shì)表位與上述MHC-Ⅰ類分子模擬對(duì)接,觀察對(duì)接形式及相應(yīng)的對(duì)接分?jǐn)?shù)。每個(gè)9肽與MHC-Ⅰ類分子對(duì)接會(huì)產(chǎn)生100 個(gè)模擬結(jié)構(gòu),前10 個(gè)被選為高置信度對(duì)接結(jié)果。將高置信度對(duì)接結(jié)果數(shù)據(jù)標(biāo)準(zhǔn)化后,絕對(duì)值高于總體中位數(shù)的數(shù)量稱為可結(jié)合核心數(shù);對(duì)接分?jǐn)?shù)絕對(duì)值最高者稱為最優(yōu)核心結(jié)合力。將上述兩種數(shù)據(jù)利用RStudio制作出氣泡圖。

    1.6.2 HTNV Gn 15肽與MHC-Ⅱ類分子模擬對(duì)接使用EpiDOCK(http://www.ddg-pharmfac.net/epidock/EpiDockPage.html)[22]工具對(duì)具有高親和力、強(qiáng)免疫原性和種間種內(nèi)保守性的人鼠共同的15肽優(yōu)勢(shì)表位與MHC-Ⅱ類分子進(jìn)行模擬對(duì)接。15肽表位被作為7個(gè)9肽核心進(jìn)行分子模擬對(duì)接,將得到的預(yù)測(cè)分?jǐn)?shù)減去閾值即為對(duì)接分?jǐn)?shù),對(duì)接分?jǐn)?shù)>0則認(rèn)為是真正的結(jié)合物。每種情況的結(jié)合分?jǐn)?shù)最大值稱為最優(yōu)核心結(jié)合力,絕對(duì)值高于總體中位數(shù)的數(shù)量稱為可結(jié)合核心數(shù)。將上述兩種數(shù)據(jù)利用RStudio制作出氣泡圖。

    1.7 質(zhì)粒疫苗和滅活疫苗 二價(jià)HFRS 滅活疫苗(HANPUWEI?)由中國長(zhǎng)春生物制品研究所有限責(zé)任公司生產(chǎn)。該疫苗由倉鼠腎細(xì)胞中滅活和純化的HTNV 和SEOV 制成(中國無單價(jià)HFRS 疫苗)。質(zhì)粒pVAX-Gn疫苗為本課題組在先前的研究中合成[23-24]。

    1.8 動(dòng)物分組及免疫 20只無特定病原體的近交系8 周齡雌性BALB/c 小鼠購自空軍軍醫(yī)大學(xué)實(shí)驗(yàn)動(dòng)物中心。BALB/c小鼠位于MHC單倍型H2-d內(nèi)。將小鼠隨機(jī)分為滅活疫苗組(n=10)與pVAX-Gn 疫苗組(n=10)。滅活疫苗組小鼠接種50 μl 1:5 稀釋的二價(jià)滅活疫苗(長(zhǎng)春生物制品研究所有限責(zé)任公司,參照說明書確定免疫劑量),pVAX-Gn 疫苗組小鼠注射50 μl 質(zhì)粒pVAX-Gn 疫苗,均為每3 周注射1 次,共注射2 次。第1 次免疫后5 周通過吸入CO2處死小鼠,采用ELISpot實(shí)驗(yàn)評(píng)估細(xì)胞免疫反應(yīng)性。所有動(dòng)物均飼養(yǎng)在無病原體、安靜的動(dòng)物實(shí)驗(yàn)室內(nèi),用標(biāo)準(zhǔn)飼料和過濾水喂養(yǎng),每周更換墊料,小鼠均正常進(jìn)食、飲水,確保10 h以上的連續(xù)夜間時(shí)間以維持晝夜節(jié)律。本動(dòng)物實(shí)驗(yàn)獲空軍軍醫(yī)大學(xué)動(dòng)物中心倫理委員會(huì)審批(2021-1201),實(shí)驗(yàn)過程符合國家及單位有關(guān)動(dòng)物管理和使用的規(guī)定。

    1.9 ELISpot實(shí)驗(yàn)驗(yàn)證表位的細(xì)胞免疫反應(yīng)性 由于其生產(chǎn)中固有的困難,僅對(duì)HTNV Gn上所預(yù)測(cè)出的5 個(gè)15 肽表位進(jìn)行合成(ChinaPeptides,中國上海),用20 μg/ml 的PBS 稀釋單個(gè)肽用于γ 干擾素(INF-γ)ELISpot 分析,4 ℃條件下在包被有IFN-γ 特異性捕獲抗體的ELISpot板上過夜。通過吸入CO2處死小鼠并解剖脾臟組織,研磨制成(0.1~2.0)×106個(gè)/ml 濃度的單細(xì)胞懸液。紅細(xì)胞裂解后,沖洗脾細(xì)胞并重新懸浮。在室溫下用含10%胎牛血清的RPMI 1640培養(yǎng)基封閉ELISpot板2 h后,將1×106個(gè)脾細(xì)胞添加至每個(gè)孔中,用終濃度為5 μg/ml 的合成Gn 肽刺激脾細(xì)胞,并將平板在37 ℃、5% CO2培養(yǎng)箱中培養(yǎng)24 h,其中,陰性對(duì)照為完整的培養(yǎng)基,陽性對(duì)照為10 μg/ml 的刀豆蛋白A(Con A)。孵育后,將培養(yǎng)板用去離子水和PBST(含0.05% Tween-20的PBS)分別洗2~3 次,孔內(nèi)加入100 μl 終濃度為2 μg/ml 生物素化抗IFN-γ抗體,室溫孵育2 h。用PBST洗滌后加入鏈霉親和素-辣根過氧化物酶(streptavidin-HRP)室溫孵育1 h。PBST 洗4 次后加入AEC(DAKEWEI,China)底物顯色,并用去離子水洗滌終止反應(yīng)。待96孔板風(fēng)干后用CTL ELISpot讀取器(美國CTL公司)對(duì)板中的斑點(diǎn)進(jìn)行計(jì)數(shù)。所有結(jié)果均顯示為每106個(gè)脾細(xì)胞的斑點(diǎn)形成單位(SFU)平均值。采用GraphPad Prism 9繪制條形圖。

    1.10 統(tǒng)計(jì)學(xué)處理 采用GraphPad Prism 9 進(jìn)行統(tǒng)計(jì)分析。本研究中僅ELISpot 實(shí)驗(yàn)涉及到統(tǒng)計(jì)學(xué)檢驗(yàn),各組數(shù)據(jù)均以±s 表示,組間差異比較采用t檢驗(yàn)。P<0.05為差異有統(tǒng)計(jì)學(xué)意義。

    2 結(jié) 果

    2.1 HTNV Gn上小鼠H-2和人HLA-Ⅰ、Ⅱ類限制性表位的親和力、免疫原性及保守性分析 表1 列出了預(yù)測(cè)工具SMMPMBEC、IEDB、NetMHCⅡpan3.2、NetMHCpan4.1、SYFPEITHI、Rankpep 預(yù)測(cè)的Gn 優(yōu)勢(shì)表位數(shù)目。重新計(jì)算排除重復(fù)9 肽和15 肽后,在H-2亞型中獲得了61個(gè)優(yōu)勢(shì)表位,在HLA-Ⅰ亞型中得到234個(gè)優(yōu)勢(shì)表位,HLA-Ⅱ亞型中得到212個(gè)優(yōu)勢(shì)表位。在H-2 亞型中,H-2-Kb 優(yōu)勢(shì)表位數(shù)目最多,有12 個(gè)9 肽。在HLA-Ⅰ類亞型中,HLA-B7 優(yōu)勢(shì)表位數(shù)目最多,有27 個(gè)9 肽。在同等條件下,免疫原性強(qiáng)的表位更有可能充分誘導(dǎo)免疫應(yīng)答。表1 列出了小鼠和人對(duì)Gn的強(qiáng)免疫原性表位數(shù)目。排除重復(fù)的表位后,在H-2 亞型中獲得了23 個(gè)強(qiáng)免疫原性優(yōu)勢(shì)表位,在HLA-Ⅰ亞型中獲得了110 個(gè)強(qiáng)免疫原性優(yōu)勢(shì)表位,HLA-Ⅱ亞型中獲得了42個(gè)強(qiáng)免疫原性優(yōu)勢(shì)表位。

    表1 小鼠H-2 和人HLA 分子對(duì)Gn 表位的親和力、免疫原性和保守性分析結(jié)果Tab.1 Analysis results of affinity, immunogenicity and conservation of mouse H-2 and human HLA molecules to Gn epitopes

    采用Blastp工具預(yù)測(cè)上述高親和力、強(qiáng)免疫原性9肽和15肽優(yōu)勢(shì)表位的種間種內(nèi)保守性,結(jié)果如表1所示。篩選獲得MHC-Ⅰ類限制性表位6 個(gè),MHC-Ⅱ類表位81 個(gè)。其中,上述MHC-Ⅰ、Ⅱ類限制性表位中存在3個(gè)交叉表位。

    各個(gè)基因型篩選出來的高親和力、強(qiáng)免疫原性優(yōu)勢(shì)表位比對(duì)結(jié)果顯示,DRB1基因型與各個(gè)基因型的高親和力、強(qiáng)免疫原性表位均有重疊,與H2-A、H2-E基因型的重疊數(shù)目最多(圖1),而H2-Ld、DPB1等基因型與其他基因型重疊表位較少。

    圖1 HTNV Gn上各MHC超家族高親和力、強(qiáng)免疫原性表位的交叉反應(yīng)性Fig.1 Cross-reactivity of MHC superfamilies with high affinity and strong immunogenicity on HTNV Gn

    2.2 HTNV Gn 9 肽、15 肽表位與MHC-Ⅰ、Ⅱ類分子的雙向?qū)哟尉垲惤沂究缁蛐烷g交叉免疫反應(yīng)性

    2.2.1 HTNV Gn 9 肽與MHC-Ⅰ類分子的相互作用關(guān)系 如圖2所示,57個(gè)MHC-Ⅰ Gn類分子被分為3組:HLA-Ⅰ亞型組、交叉反應(yīng)組(HLA 主要)及交叉反應(yīng)組(H2 主要)。在HLA-Ⅰ亞型組中,HLA-A3201與HLA-B58(-B5701、-B5801)得分相似,表明HLAA3201 基因型在提呈HTNV Gn 抗原過程中表現(xiàn)出HLA-B58 樣特征;HLA-A2601 與HLA-A3401 歸屬于同一類別,兩個(gè)基因型得分相似;HLA-B4601 與HLA-B15(-1502,-1513,-1521)得分相似,表明HLAB4601 基因型在提呈HTNV Gn 抗原過程中表現(xiàn)出HLA-B15 樣特征;HLA-A6801、HLA-A3301、HLAA0301、HLA-A1101 得分相似,歸屬于同一類別;HLA-B0801、HLA-B1405、HLA-B3903 得分歸屬于同一類;HLA-B5301 與HLA-B3501、HLA-B3508歸屬于同一類別。對(duì)于交叉反應(yīng)組(H2 主要),H2-Ld 得分與HLA-B0702相似,歸屬于同一類;H2-Db、H2-Dd和H2-Kb得分體現(xiàn)出在鼠種內(nèi)提呈HTNV Gn的相似性。在交叉反應(yīng)組(HLA 主要)中,H2-Kd 得分與HLA-A24(-A2301、-A2402、-A2407)相似,表明MHC在提呈9肽表位上的跨基因型相似性。

    2.2.2 HTNV Gn 15 肽與MHC-Ⅱ類分子的相互作用關(guān)系 圖3顯示了不同MHC-Ⅱ亞型間的差異。51個(gè)MHC-Ⅱ類分子被分為3組:HLA-Ⅱ亞型組、交叉反應(yīng)組(HLA主要)及交叉反應(yīng)組(H2主要)。在HLA-Ⅱ亞型組中, DRB1(-0407、 -0412、 -0803、 -1402、-1409、-0304、-0301、-0101、-0202、-1202、-1401、-1408、 -1104、 -1301、 -1304) 和 HLA-DQA10101-DQB10501、 HLA-DPA10103(-DPB10101、 -DPB10201、-DPB10401、-DPB10402)歸屬于同一類別,上述基因型得分相似;DRB1(-1302、-1201、-1501)和HLADQA10501-DQB10201、 HLA-DQA10301-DQB10301、HLA-DQA10401-DQB10402、HLA-DQA10501-DQB10301 基因型相近,表明其在提呈HTNV Gn 抗原過程中表現(xiàn)出相同的特征;DRB1(-0704、-1602、-0701、-0101、-0901、-0401、-0405、-0802)和HLADPA10103-DPB11401 歸屬于同一類別,基因型得分相似。對(duì)于交叉反應(yīng)組(H2 主要),H2-IAd 得分與HLA-DQA10102-DQB10602相似,歸屬于同一類,具有相似的抗原呈遞結(jié)果。在交叉反應(yīng)組(HLA 主要)中,H2-IEK與DRB0101基因型得分相似。

    2.3 MHC分子與表位的模擬分子對(duì)接驗(yàn)證親和力

    2.3.1 HTNV Gn 9 肽表位與MHC-Ⅰ類分子對(duì)接圖4A 所示為每個(gè)肽-MHC 對(duì)接后的可結(jié)合核心數(shù)及最優(yōu)核心結(jié)合力的氣泡圖,3個(gè)9肽優(yōu)勢(shì)表位均與多個(gè)MHC 分子對(duì)接較好。表位RYKSRCYIF 與HLAB5801、HLA-A0206、HLA-B5301 等基因型對(duì)接的最優(yōu)核心結(jié)合力較為良好,表明該表位與上述基因型分子結(jié)合的牢固程度更高,同時(shí),該表位與HLAB5801、HLA-A0206、HLA-B5301、HLA-B0801 等基因型分子對(duì)接可結(jié)合核心數(shù)較多,表明其與基因型分子對(duì)接的概率較大; 表位VVYERTYCM、QVVYERTYC與H2-Kb和HLA-B4001基因型對(duì)接情況較弱,表明其與上述基因型分子對(duì)接的概率及對(duì)接牢固程度均較低。此外,HLA-B5801、HLA-A0206、HLA-B5301、HLA-B0801 基因型對(duì)3 個(gè)優(yōu)勢(shì)表位的可結(jié)合核心數(shù)多,表明上述基因型分子對(duì)優(yōu)勢(shì)表位的結(jié)合概率較大??傮w而言,3 個(gè)優(yōu)勢(shì)表位與人HLA基因型對(duì)接時(shí)的最優(yōu)核心結(jié)合力和可結(jié)合核心數(shù)均優(yōu)于小鼠H-2基因型。如圖4B中b、d所示,分別選取表位RYKSRCYIF 與H2-Kb、QVVYERTYC 與H2-Ld作為對(duì)接優(yōu)良的例子,通過相應(yīng)的對(duì)接模擬結(jié)構(gòu)可以看出,在上述兩個(gè)對(duì)接過程的10種對(duì)接模擬結(jié)構(gòu)中,9 肽均位于MHC 分子抗原結(jié)合槽內(nèi)。選取表位QVVYERTYC 與H2-Kb、表位VVYERTYCM 與H2-Kb作為對(duì)接較差的例子,如圖4B中a、c所示,上述兩種對(duì)接模擬結(jié)果中,多數(shù)9 肽模擬狀態(tài)均位于MHC分子抗原結(jié)合槽外,其中有一個(gè)游離于結(jié)合槽;同樣的情況見于表位VVYERTYCM 與H2-Kb 對(duì)接結(jié)果中。通過比較上述分子對(duì)接模型,從側(cè)面驗(yàn)證了分子模擬對(duì)接存在分子結(jié)構(gòu)上的基礎(chǔ),與分子對(duì)接結(jié)合可保持一致。

    2.3.2 HTNV Gn 15 肽優(yōu)勢(shì)表位與MHC-Ⅱ類分子模擬對(duì)接 如圖5 所示,5 個(gè)MHC-Ⅱ類限制性優(yōu)勢(shì)表位均與多個(gè)MHC基因型分子對(duì)接良好,其中,與表位GGIFNITSSMCLVSK 對(duì)接良好的基因型最多。此外,DPA1*0201/DPB1*0101 基因型與優(yōu)勢(shì)表位對(duì)接結(jié)果均較差。

    圖5 HTNV Gn 15肽表位與MHC-Ⅱ類分子對(duì)接分?jǐn)?shù)氣泡圖Fig.5 The bubble map of HTNV Gn 15-mer epitope and MHC-Ⅰ molecule docking

    2.4 ELISpot實(shí)驗(yàn)評(píng)估15肽表位與MHC-Ⅱ類分子對(duì)接的細(xì)胞免疫反應(yīng)性 ELISpot實(shí)驗(yàn)結(jié)果顯示,在滅活疫苗組及pVAX-Gn 組中,5 個(gè)表位均可引起脾細(xì)胞分泌IFN-γ(圖6),提示通過表位篩選促進(jìn)細(xì)胞免疫應(yīng)答確實(shí)行之有效。同時(shí),兩組的SFU 平均值接近,提示Gn 上的3 個(gè)人鼠共優(yōu)勢(shì)15 肽表位由Gn 和GP全長(zhǎng)所引起的免疫應(yīng)答效果相近。

    圖6 5個(gè)15肽表位ELISpot實(shí)驗(yàn)結(jié)果Fig.6 ELISpot experiment results of 5 15-mer epitopes

    3 討 論

    對(duì)于HTNV 的預(yù)防與控制,現(xiàn)有雙價(jià)滅活疫苗已被證實(shí)可誘導(dǎo)抗病毒應(yīng)答,但在中和抗體效價(jià)、細(xì)胞免疫應(yīng)答和長(zhǎng)效免疫記憶方面尚存在不足,同時(shí)在T 細(xì)胞免疫應(yīng)答水平方面也有待提高。本研究以比較免疫學(xué)和結(jié)構(gòu)生物學(xué)原理為基礎(chǔ),選取可特異性激活抗HTNV 免疫應(yīng)答的Gn 分子為研究對(duì)象,對(duì)9 肽和15 肽表位免疫反應(yīng)性進(jìn)行評(píng)估,最終獲得高親和力、強(qiáng)免疫原性、種間種內(nèi)保守的MHC-Ⅰ類限制性優(yōu)勢(shì)表位6個(gè),MHC-Ⅱ類限制性優(yōu)勢(shì)表位81個(gè)。使用上述方法可精準(zhǔn)預(yù)測(cè)出HTNV Gn上細(xì)胞免疫反應(yīng)性表位,在短時(shí)間內(nèi)找到具有潛在疫苗作用靶點(diǎn)的抗原表位,極大地節(jié)省了實(shí)驗(yàn)成本。同時(shí),該方法已用于近年流行的新型冠狀病毒等病毒的研究[25-26]。

    人類遺傳因素也可影響對(duì)漢坦病毒的易感性及其所致疾病的嚴(yán)重性[27]。有研究發(fā)現(xiàn)了HTNV 感染患者中HLA基因型的連鎖性,HLA-B*46-DRB1*09和HLA-B*51-DRB1*09 基因型在中國漢族危重HFRS 患者中更常見[28-29]。本研究比較了HLA-Ⅰ類基因型和HLA-Ⅱ類基因型之間HTNV GP 免疫反應(yīng)性表位的相關(guān)性,發(fā)現(xiàn)DRB1基因型與HLA-B15、-B46基因型相同的免疫反應(yīng)性表位較多,提示DRB1-B15、DRB1-B46的化學(xué)結(jié)構(gòu)基礎(chǔ)相似,具有上述基因的不同地區(qū)、不同種族的人群可能對(duì)漢灘病毒感染具有相近的免疫應(yīng)答反應(yīng),同一表位引起的效應(yīng)可能相近。同時(shí),本研究對(duì)小鼠基因型與人類基因型的表位進(jìn)行對(duì)比發(fā)現(xiàn),人類DRB1 基因型與小鼠H2-A、H2-E 基因型部分表位更為相近,HLA-A02 基因型與H2-Kd 基因型在提呈HTNV GP 上9 肽表位時(shí)親和力數(shù)據(jù)呈現(xiàn)相似的趨勢(shì),提示在不同物種存在對(duì)GP表位相似的反應(yīng),具有上述基因型的小鼠與人類在漢灘病毒感染及表位免疫刺激后可能產(chǎn)生相近的免疫應(yīng)答反應(yīng),提示在缺乏人源化轉(zhuǎn)基因小鼠的情況下,使用特定基因型的小鼠可能是實(shí)驗(yàn)?zāi)P偷挠杏锰娲颷16,30]。

    近年來,隨著生物信息學(xué)的迅速發(fā)展,工具與算法不斷更新。不同工具和算法的原理大相徑庭,但目的均是對(duì)MHC表位的親和力進(jìn)行模擬分析。本研究將多個(gè)算法進(jìn)行整合運(yùn)用,保證了分析的高準(zhǔn)確率和泛MHC覆蓋率。然而,正因如此,本研究更注重于結(jié)果的獲取,對(duì)其中的算法過程考慮較少,因而衍生出以下問題:(1)不同運(yùn)算過程將賦予表位一些意料之外的特殊內(nèi)涵;(2)多種算法多個(gè)工具同時(shí)使用可造成多重檢驗(yàn)問題;(3)由于多重檢驗(yàn)造成的假陰性、假陽性和過擬合等;(4)生物信息學(xué)算法和數(shù)據(jù)庫的更新及對(duì)表位結(jié)果的影響;(5)調(diào)整工具優(yōu)先級(jí)、參數(shù)設(shè)置以提高預(yù)測(cè)準(zhǔn)確率的問題。這些問題背后的含義,將是后續(xù)深入挖掘探索生物信息學(xué)的方向。

    綜上所述,本研究預(yù)測(cè)并驗(yàn)證了HTNV Gn上可激活細(xì)胞免疫反應(yīng)性的表位,證實(shí)在MHC提呈中病毒抗原存在跨基因、種屬、物種的交叉免疫反應(yīng)性。然而,本研究?jī)H在小鼠脾細(xì)胞中進(jìn)行了表位驗(yàn)證,其在人體內(nèi)的免疫反應(yīng)性仍需進(jìn)一步探討。與Dawood[25]、Araújo 等[26]通過數(shù)據(jù)庫預(yù)測(cè)獲得SARSCoV-2 結(jié)構(gòu)蛋白相關(guān)的強(qiáng)免疫原性、種間種內(nèi)保守的候選表位相比,本研究創(chuàng)造性地進(jìn)行種屬交叉分析和實(shí)驗(yàn)驗(yàn)證,提高了抗原優(yōu)勢(shì)表位篩選的系統(tǒng)性和可靠性,為在免疫遺傳學(xué)背景下開發(fā)基于HTNV Gn分子特異性群體免疫保護(hù)性表位疫苗提供了理論基礎(chǔ)。

    猜你喜歡
    免疫原性表位親和力
    H3N2流感病毒HA保守Th表位對(duì)CD4+T細(xì)胞活化及分化的影響
    高端訪談節(jié)目如何提升親和力
    新聞傳播(2018年11期)2018-08-29 08:15:30
    高端訪談節(jié)目如何提升親和力探索
    新聞傳播(2018年13期)2018-08-29 01:06:52
    豬乙型腦炎PrM-E重組腺病毒的構(gòu)建及免疫原性
    聯(lián)合T、B細(xì)胞表位設(shè)計(jì)多肽疫苗的研究進(jìn)展①
    親和力在播音主持中的作用探究
    新聞傳播(2016年9期)2016-09-26 12:20:34
    牛傳染性鼻氣管炎IBRV/JZ06-3基礎(chǔ)毒株的安全性和免疫原性研究
    布魯菌缺失疫苗株M5-ΔznuA和M5-Δbp26-ΔznuA的構(gòu)建及毒力和免疫原性的評(píng)估
    將親和力應(yīng)用于播音主持中的方法探討
    新聞傳播(2015年7期)2015-07-18 11:09:57
    小反芻獸疫病毒化學(xué)合成表位多肽對(duì)小鼠的免疫效果研究
    久久97久久精品| 亚州av有码| 人体艺术视频欧美日本| 久久久久久久大尺度免费视频| 久久久久人妻精品一区果冻| 黄色一级大片看看| 另类精品久久| 国产精品久久久久久av不卡| 久久久国产一区二区| 成人国产麻豆网| 亚洲美女黄色视频免费看| 中文字幕人妻丝袜制服| .国产精品久久| 伦理电影大哥的女人| 成年女人在线观看亚洲视频| 欧美丝袜亚洲另类| 中文字幕精品免费在线观看视频 | 欧美xxxx性猛交bbbb| 六月丁香七月| 一区在线观看完整版| 国产欧美亚洲国产| 久久久久久久亚洲中文字幕| 午夜激情福利司机影院| 插阴视频在线观看视频| 日本爱情动作片www.在线观看| 国产亚洲最大av| 啦啦啦在线观看免费高清www| 国产男人的电影天堂91| 菩萨蛮人人尽说江南好唐韦庄| 日韩人妻高清精品专区| 女人久久www免费人成看片| 亚洲av成人精品一区久久| 嘟嘟电影网在线观看| 精品卡一卡二卡四卡免费| av卡一久久| 精品国产一区二区三区久久久樱花| 晚上一个人看的免费电影| 国产美女午夜福利| 午夜老司机福利剧场| 午夜91福利影院| 久久久久久久精品精品| 多毛熟女@视频| 国产精品成人在线| 精品国产露脸久久av麻豆| 久久狼人影院| 自线自在国产av| 下体分泌物呈黄色| 国产在视频线精品| 亚洲欧美中文字幕日韩二区| 99热全是精品| 三级国产精品片| 精品少妇内射三级| 啦啦啦中文免费视频观看日本| av天堂久久9| a级毛片免费高清观看在线播放| 日本wwww免费看| 日韩人妻高清精品专区| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜添av毛片| 久久精品国产鲁丝片午夜精品| 丝瓜视频免费看黄片| 国产中年淑女户外野战色| 六月丁香七月| 18禁在线播放成人免费| 男女边吃奶边做爰视频| 日韩大片免费观看网站| 亚洲成色77777| 欧美变态另类bdsm刘玥| 久久精品国产鲁丝片午夜精品| 色婷婷久久久亚洲欧美| 午夜福利视频精品| 久久97久久精品| 丰满乱子伦码专区| 91aial.com中文字幕在线观看| 免费观看av网站的网址| 久久久久久久久久久久大奶| 人人妻人人澡人人爽人人夜夜| 色5月婷婷丁香| 一级爰片在线观看| 亚洲成人一二三区av| 丝袜脚勾引网站| 夜夜骑夜夜射夜夜干| 亚洲,一卡二卡三卡| 青春草亚洲视频在线观看| 一边亲一边摸免费视频| 啦啦啦在线观看免费高清www| xxx大片免费视频| 爱豆传媒免费全集在线观看| 日韩中字成人| 亚洲av成人精品一二三区| 久久久久久久久大av| 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美在线一区| 日韩中字成人| 欧美日韩视频精品一区| 午夜日本视频在线| 日韩,欧美,国产一区二区三区| 97在线人人人人妻| 三上悠亚av全集在线观看 | 久久精品熟女亚洲av麻豆精品| 国产高清三级在线| 久久久久久久久久久丰满| 十八禁高潮呻吟视频 | 少妇被粗大猛烈的视频| 国产精品无大码| 亚洲精品,欧美精品| 在线观看av片永久免费下载| 国产成人精品久久久久久| 99精国产麻豆久久婷婷| 成人亚洲精品一区在线观看| 性高湖久久久久久久久免费观看| 久久99一区二区三区| 婷婷色综合大香蕉| 色婷婷av一区二区三区视频| 亚洲综合色惰| 丰满人妻一区二区三区视频av| 综合色丁香网| 午夜老司机福利剧场| 性色avwww在线观看| 国产精品熟女久久久久浪| 免费观看的影片在线观看| 亚洲精品久久午夜乱码| 成人特级av手机在线观看| 日韩中文字幕视频在线看片| 五月玫瑰六月丁香| 精品久久久久久电影网| 国产一区亚洲一区在线观看| 一边亲一边摸免费视频| 成人漫画全彩无遮挡| 啦啦啦中文免费视频观看日本| 黄色怎么调成土黄色| 欧美国产精品一级二级三级 | 亚洲欧美一区二区三区国产| 麻豆成人午夜福利视频| 国产精品嫩草影院av在线观看| av女优亚洲男人天堂| 国产亚洲一区二区精品| 国产美女午夜福利| 久久精品久久久久久久性| 蜜桃久久精品国产亚洲av| 麻豆精品久久久久久蜜桃| 纯流量卡能插随身wifi吗| 国产欧美日韩一区二区三区在线 | 在现免费观看毛片| 嫩草影院入口| 一本一本综合久久| 熟女人妻精品中文字幕| 在线 av 中文字幕| 免费黄色在线免费观看| 香蕉精品网在线| 精品人妻一区二区三区麻豆| 国产免费又黄又爽又色| 熟女av电影| 亚洲成人av在线免费| 在线观看av片永久免费下载| 男女边摸边吃奶| 亚洲精品乱码久久久久久按摩| 22中文网久久字幕| av网站免费在线观看视频| 这个男人来自地球电影免费观看 | 久久99热这里只频精品6学生| 在线观看国产h片| 国产高清国产精品国产三级| 成人特级av手机在线观看| 激情五月婷婷亚洲| 国产日韩一区二区三区精品不卡 | 久久久亚洲精品成人影院| 91午夜精品亚洲一区二区三区| 91午夜精品亚洲一区二区三区| 亚洲av中文av极速乱| 亚洲精品乱码久久久久久按摩| 亚洲真实伦在线观看| 观看av在线不卡| 看非洲黑人一级黄片| 日韩精品有码人妻一区| 亚洲怡红院男人天堂| 久久久久久伊人网av| 爱豆传媒免费全集在线观看| 黑人猛操日本美女一级片| 国产综合精华液| 亚洲欧美日韩卡通动漫| 18禁在线无遮挡免费观看视频| 一区二区三区免费毛片| 2021少妇久久久久久久久久久| 久久精品夜色国产| av网站免费在线观看视频| 色94色欧美一区二区| 最近中文字幕2019免费版| 天堂中文最新版在线下载| 永久免费av网站大全| 久久99一区二区三区| 亚洲精品亚洲一区二区| 男女边摸边吃奶| 深夜a级毛片| 最黄视频免费看| www.色视频.com| 精品久久久久久久久av| 成人漫画全彩无遮挡| 精品少妇久久久久久888优播| 精品国产国语对白av| 精品卡一卡二卡四卡免费| 热99国产精品久久久久久7| 欧美xxxx性猛交bbbb| 免费看av在线观看网站| 国产在线男女| 亚洲精品视频女| 制服丝袜香蕉在线| 这个男人来自地球电影免费观看 | 91成人精品电影| 国产精品人妻久久久久久| 欧美三级亚洲精品| 视频区图区小说| 综合色丁香网| 日本黄色片子视频| 精品少妇黑人巨大在线播放| 边亲边吃奶的免费视频| 一区二区av电影网| 久久精品久久精品一区二区三区| 下体分泌物呈黄色| 老女人水多毛片| 成年人免费黄色播放视频 | 久久久久久久久大av| 亚洲国产精品一区三区| 亚洲精品国产成人久久av| 岛国毛片在线播放| 我的老师免费观看完整版| 亚洲一级一片aⅴ在线观看| 18禁在线无遮挡免费观看视频| 成年av动漫网址| 亚洲在久久综合| 永久免费av网站大全| 免费播放大片免费观看视频在线观看| 亚洲国产精品一区三区| 亚洲国产色片| 国产亚洲av片在线观看秒播厂| 精品少妇黑人巨大在线播放| 男女边摸边吃奶| 多毛熟女@视频| 亚洲欧美中文字幕日韩二区| 免费少妇av软件| 国产精品一区二区性色av| 国产伦精品一区二区三区视频9| 欧美日韩视频高清一区二区三区二| 日韩 亚洲 欧美在线| 秋霞伦理黄片| 伊人亚洲综合成人网| 各种免费的搞黄视频| 一二三四中文在线观看免费高清| 免费看日本二区| 人妻一区二区av| 久久99精品国语久久久| 亚洲精品色激情综合| 丰满少妇做爰视频| 高清在线视频一区二区三区| 亚洲欧美日韩东京热| 一区二区三区精品91| 在线观看免费日韩欧美大片 | 美女国产视频在线观看| 日本av手机在线免费观看| 建设人人有责人人尽责人人享有的| 大又大粗又爽又黄少妇毛片口| 国产男女超爽视频在线观看| 亚洲熟女精品中文字幕| 久久综合国产亚洲精品| 久久99热6这里只有精品| 赤兔流量卡办理| 久久久久久久久久人人人人人人| 亚洲精品日韩av片在线观看| 成人18禁高潮啪啪吃奶动态图 | 国产男女内射视频| 国产亚洲最大av| 精品一区二区免费观看| 高清不卡的av网站| 啦啦啦啦在线视频资源| 蜜桃在线观看..| 91久久精品电影网| 亚洲精品色激情综合| av卡一久久| 国产乱来视频区| 在线观看人妻少妇| 99国产精品免费福利视频| 伦精品一区二区三区| 久久久国产精品麻豆| 日韩av不卡免费在线播放| 日韩大片免费观看网站| 精品一区在线观看国产| 国产片特级美女逼逼视频| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久久性| 国产精品99久久久久久久久| 日本黄大片高清| 99久久精品热视频| 黄色毛片三级朝国网站 | 成年美女黄网站色视频大全免费 | av免费观看日本| 日本爱情动作片www.在线观看| 久久国产亚洲av麻豆专区| 国产黄色视频一区二区在线观看| 国产91av在线免费观看| 日韩欧美一区视频在线观看 | 中文字幕精品免费在线观看视频 | 精品久久久久久电影网| 91久久精品国产一区二区成人| 韩国av在线不卡| 精品一品国产午夜福利视频| av.在线天堂| 免费不卡的大黄色大毛片视频在线观看| 美女福利国产在线| 中文字幕亚洲精品专区| 国产成人freesex在线| 亚洲美女黄色视频免费看| 国产黄片视频在线免费观看| 青春草国产在线视频| 熟女av电影| 国产黄频视频在线观看| 日韩强制内射视频| 三级国产精品片| 国产在线视频一区二区| 国产真实伦视频高清在线观看| 国产高清国产精品国产三级| 国产精品久久久久久精品电影小说| 老熟女久久久| 女性生殖器流出的白浆| 26uuu在线亚洲综合色| 国产日韩一区二区三区精品不卡 | 国产成人91sexporn| 校园人妻丝袜中文字幕| 一级黄片播放器| 80岁老熟妇乱子伦牲交| av播播在线观看一区| 午夜91福利影院| 一个人免费看片子| 免费观看在线日韩| 女人精品久久久久毛片| 少妇猛男粗大的猛烈进出视频| 啦啦啦中文免费视频观看日本| 国产精品久久久久成人av| 国产日韩欧美亚洲二区| 性高湖久久久久久久久免费观看| 欧美日韩精品成人综合77777| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲网站| 国产伦精品一区二区三区四那| 日本-黄色视频高清免费观看| 亚洲伊人久久精品综合| 黄色毛片三级朝国网站 | 内射极品少妇av片p| 国产永久视频网站| 黑人高潮一二区| 久久97久久精品| 三上悠亚av全集在线观看 | 一本色道久久久久久精品综合| 偷拍熟女少妇极品色| 久久国内精品自在自线图片| 免费看日本二区| 国产精品国产av在线观看| 永久免费av网站大全| 色哟哟·www| 欧美丝袜亚洲另类| 色94色欧美一区二区| 午夜久久久在线观看| 久久久久久久大尺度免费视频| 国产白丝娇喘喷水9色精品| 我的女老师完整版在线观看| 97精品久久久久久久久久精品| av有码第一页| 日本猛色少妇xxxxx猛交久久| 亚洲精品中文字幕在线视频 | 亚洲第一区二区三区不卡| 精品久久久久久久久av| 老司机亚洲免费影院| 国产色婷婷99| 女人久久www免费人成看片| 一区在线观看完整版| 在线观看国产h片| 国产日韩欧美视频二区| 在线观看人妻少妇| 国产精品一区二区在线观看99| 日本午夜av视频| 国产成人免费无遮挡视频| 精品亚洲成a人片在线观看| 永久网站在线| 我要看日韩黄色一级片| 精品99又大又爽又粗少妇毛片| 亚洲av二区三区四区| 中文精品一卡2卡3卡4更新| 水蜜桃什么品种好| 欧美高清成人免费视频www| 国产一区二区三区av在线| 一级毛片电影观看| 亚洲av不卡在线观看| 两个人的视频大全免费| 丰满迷人的少妇在线观看| 久久精品国产亚洲av天美| 亚洲国产精品一区二区三区在线| 啦啦啦在线观看免费高清www| 三级国产精品片| 色94色欧美一区二区| 97超视频在线观看视频| 另类精品久久| 亚洲欧美精品自产自拍| 中国国产av一级| 久久久久久久久久久丰满| 纵有疾风起免费观看全集完整版| 久久久久人妻精品一区果冻| 中文字幕制服av| 永久网站在线| 欧美日韩视频高清一区二区三区二| 中文字幕亚洲精品专区| 肉色欧美久久久久久久蜜桃| 青春草亚洲视频在线观看| 亚洲精品中文字幕在线视频 | 不卡视频在线观看欧美| 美女国产视频在线观看| 99热6这里只有精品| 性高湖久久久久久久久免费观看| 国产亚洲欧美精品永久| 国产精品国产三级国产专区5o| 精品国产一区二区三区久久久樱花| 成人特级av手机在线观看| 美女视频免费永久观看网站| 日本-黄色视频高清免费观看| 日韩强制内射视频| 另类精品久久| 日韩人妻高清精品专区| 欧美日韩国产mv在线观看视频| 少妇的逼水好多| 十分钟在线观看高清视频www | 久久精品国产亚洲av天美| 在线看a的网站| 欧美日韩视频高清一区二区三区二| 老女人水多毛片| 这个男人来自地球电影免费观看 | 久久精品国产亚洲av涩爱| 亚洲一级一片aⅴ在线观看| 国产69精品久久久久777片| 国产高清不卡午夜福利| √禁漫天堂资源中文www| 嘟嘟电影网在线观看| 日韩欧美 国产精品| 毛片一级片免费看久久久久| 蜜桃在线观看..| 亚洲国产精品专区欧美| a 毛片基地| 亚洲中文av在线| 国产欧美另类精品又又久久亚洲欧美| 三级国产精品片| 中文资源天堂在线| 国产成人91sexporn| 青春草国产在线视频| 天天操日日干夜夜撸| 欧美精品国产亚洲| 99热这里只有是精品在线观看| 精品国产一区二区三区久久久樱花| 亚洲精品国产色婷婷电影| 精品亚洲乱码少妇综合久久| 久久热精品热| 黄色毛片三级朝国网站 | 69精品国产乱码久久久| 午夜激情福利司机影院| 日韩av免费高清视频| 一本一本综合久久| 麻豆成人av视频| 一级毛片 在线播放| 麻豆乱淫一区二区| 天堂中文最新版在线下载| 男男h啪啪无遮挡| 日日摸夜夜添夜夜添av毛片| 各种免费的搞黄视频| 九九爱精品视频在线观看| 国产精品成人在线| 乱系列少妇在线播放| 熟女电影av网| 久久 成人 亚洲| 高清视频免费观看一区二区| 啦啦啦在线观看免费高清www| 香蕉精品网在线| 2022亚洲国产成人精品| 春色校园在线视频观看| .国产精品久久| 一区二区三区四区激情视频| 久久毛片免费看一区二区三区| av免费在线看不卡| 国产中年淑女户外野战色| 九九在线视频观看精品| 色视频在线一区二区三区| 成人午夜精彩视频在线观看| 久久精品久久久久久噜噜老黄| 亚洲人与动物交配视频| 久久久久久久久久人人人人人人| 久久人人爽av亚洲精品天堂| 国产一区有黄有色的免费视频| 最黄视频免费看| 国产精品久久久久久精品电影小说| 亚洲成色77777| 丰满迷人的少妇在线观看| 夜夜爽夜夜爽视频| 国产伦精品一区二区三区视频9| 色视频www国产| 亚洲精品亚洲一区二区| 国产美女午夜福利| 高清黄色对白视频在线免费看 | 亚洲人成网站在线播| 纵有疾风起免费观看全集完整版| 嫩草影院入口| 在线精品无人区一区二区三| 国产国拍精品亚洲av在线观看| 自线自在国产av| 国产精品福利在线免费观看| 在线观看免费高清a一片| 国产高清三级在线| 亚洲国产精品一区二区三区在线| 午夜激情久久久久久久| 亚洲国产色片| 久久毛片免费看一区二区三区| 日日啪夜夜撸| 一级二级三级毛片免费看| 日本黄色片子视频| 亚洲欧美日韩东京热| 少妇丰满av| 精品国产国语对白av| 丁香六月天网| 国产成人精品无人区| 欧美国产精品一级二级三级 | 熟女av电影| av免费观看日本| 精品卡一卡二卡四卡免费| 新久久久久国产一级毛片| 在线天堂最新版资源| 亚洲精品色激情综合| 久久久久网色| 国产精品人妻久久久久久| 免费观看a级毛片全部| 91精品伊人久久大香线蕉| 新久久久久国产一级毛片| 大话2 男鬼变身卡| 黑人猛操日本美女一级片| 韩国高清视频一区二区三区| 欧美日本中文国产一区发布| 亚洲伊人久久精品综合| 在现免费观看毛片| 国产精品99久久久久久久久| 中国三级夫妇交换| 午夜激情福利司机影院| 成年人免费黄色播放视频 | 最后的刺客免费高清国语| 国产一区二区三区综合在线观看 | 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 国产精品99久久久久久久久| 美女xxoo啪啪120秒动态图| 久久久欧美国产精品| 日韩av免费高清视频| 精品亚洲成a人片在线观看| 亚洲伊人久久精品综合| 久久久久久久国产电影| av福利片在线| 成年av动漫网址| 国产极品粉嫩免费观看在线 | 精品久久久久久久久av| 大又大粗又爽又黄少妇毛片口| 久久久久久久久大av| 国产色爽女视频免费观看| 桃花免费在线播放| 久久久国产欧美日韩av| 久久久久视频综合| 日本91视频免费播放| 97超视频在线观看视频| 伊人久久国产一区二区| 女人精品久久久久毛片| 亚洲国产毛片av蜜桃av| 97超视频在线观看视频| 精品亚洲成a人片在线观看| 99久久精品热视频| 内地一区二区视频在线| 日韩制服骚丝袜av| 黄色怎么调成土黄色| 亚洲无线观看免费| 美女xxoo啪啪120秒动态图| 高清视频免费观看一区二区| 美女xxoo啪啪120秒动态图| 亚洲无线观看免费| 在线观看国产h片| 中文资源天堂在线| 美女xxoo啪啪120秒动态图| 日韩成人伦理影院| 久久午夜综合久久蜜桃| 国产成人午夜福利电影在线观看| 晚上一个人看的免费电影| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 日本与韩国留学比较| 日本黄色片子视频| 街头女战士在线观看网站| 免费黄色在线免费观看| 制服丝袜香蕉在线| 国产精品熟女久久久久浪| 精品熟女少妇av免费看| 精品卡一卡二卡四卡免费| av天堂中文字幕网| 我要看黄色一级片免费的| 国产精品久久久久久久电影| 久久久久精品性色| 91aial.com中文字幕在线观看| 免费看日本二区| 最新中文字幕久久久久| 久久久久久人妻| 亚洲美女视频黄频| 国产在线一区二区三区精| 爱豆传媒免费全集在线观看| 高清在线视频一区二区三区| 国产成人91sexporn| 观看免费一级毛片| 777米奇影视久久| 赤兔流量卡办理|