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

    斯氏副柔線蟲CAMK/TSSK蛋白激酶基因的克隆及其蛋白質結構與抗原表位預測分析

    2018-04-04 01:04:02王文龍趙學亮馮陳晨白麗艷王姝懿王夢雅劉春霞
    畜牧獸醫(yī)學報 2018年3期
    關鍵詞:斯氏表位親水性

    王文龍,趙學亮,孫 柯,馮陳晨,白麗艷,王姝懿,3,王夢雅,劉春霞

    (1. 內蒙古農(nóng)業(yè)大學獸醫(yī)學院,農(nóng)業(yè)部動物疾病臨床診療技術重點實驗室,呼和浩特 010018;2. 內蒙古農(nóng)業(yè)大學動物科學學院,呼和浩特 010018;3. 內蒙古自治區(qū)綜合疾病預防控制中心,呼和浩特 010018; 4. 內蒙古農(nóng)業(yè)大學生命科學學院,呼和浩特010018)

    駱駝斯氏副柔線蟲(Parabronemaskrjabini)屬于旋尾目,副柔屬[1-2],是一種寄生在駱駝、牛、羊等反芻動物真胃引起疾病的寄生蟲,其中駱駝是其最適宜的終末宿主[3-4]。當駱駝大量感染斯氏副柔線蟲時會導致腹瀉、生長遲緩,真胃黏膜潰瘍、出血,淋巴細胞浸潤等現(xiàn)象,嚴重者會導致貧血死亡。國內自楊曉野[5]最早發(fā)現(xiàn)斯氏副柔線蟲病以來,經(jīng)流行病學調查發(fā)現(xiàn)斯氏副柔線蟲已廣泛存在,并且在不同地區(qū)的感染率為72%~100%。據(jù)報道,內蒙古地區(qū)雙峰駝斯氏副柔線蟲感染率高達100%,感染強度最高可達7 611條[6]。目前,斯氏副柔線蟲病仍然嚴重危害我國牧區(qū)駱駝的健康,制約著當?shù)仞B(yǎng)駝業(yè)的發(fā)展,給農(nóng)牧民造成巨大的經(jīng)濟損失。

    由于駱駝主要分布在經(jīng)濟欠發(fā)達國家或地區(qū),導致斯氏副柔線蟲病的相關研究資料較少,尤其是駱駝斯氏副柔線蟲病致病作用的資料極度缺乏。目前,國內外對于駱駝斯氏副柔線蟲病的研究主要在于形態(tài)學鑒定、流行病學調查等方面[7],尚無有效疫苗用于該病的預防,也未見報道該病有效的生前診斷方法。目前,在生產(chǎn)實踐中,往往是在沒有診斷依據(jù)的情況下盲目給藥驅蟲,不僅造成經(jīng)濟損失,而且也會導致藥物殘留超標和耐藥蟲株的產(chǎn)生。因此,對該病的研究應致力于在探索斯氏副柔線蟲病生物學特征的基礎上,尋找有效的生前診斷方法和制定新的防治措施。

    本試驗選取駱駝斯氏副柔線蟲為研究對象,參考已報道的線蟲免疫學診斷抗原的特點,在斯氏副柔線蟲轉錄組測序分析的基礎上[8-9],篩選出與斯氏副柔線蟲免疫相關的基因,通過Blastx比對確定駱駝斯氏副柔線蟲免疫相關基因CTPK(CAMK/TSSK protein kinase)。通過克隆斯氏副柔線蟲CTPK基因的CDS區(qū),利用生物信息學方法對該基因結構、理化性質及抗原表位等進行分析和預測,以期找到優(yōu)勢抗原表位,為駱駝斯氏副柔線蟲病生前診斷方法iELISA的建立和DNA疫苗的研究奠定了理論基礎。

    1 材料與方法

    1.1 試驗材料

    1.1.1蟲體和質粒駱駝斯氏副柔線蟲樣本采自于巴彥淖爾市烏拉特后旗雙峰駝真胃。在體視鏡下根據(jù)蟲體不同形態(tài)特征鑒定雌雄后分裝到2 mL凍存管中,標記后放入液氮中由本實驗室保存,用于RNA的提取。

    試驗所用菌株:E.coliDH5α,北京全式金(TransGen Biotech)生物技術有限公司;質粒載體:克隆載體pMD19-T simple vector,大連寶生物工程有限公司。

    1.1.2主要試劑與儀器Trizol 通用型RNA快速提取試劑盒、高保真酶Prime STAR?GXL DNA Polymerase、標準相對分子質量DNA Marker購自TaKaRa公司;DNA凝膠回收試劑盒,購自AXYGEN公司,DEPC處理水。NanoDROP核酸濃度測定儀,ND2000美國Thermo公司;凝膠影像儀,北京賽智創(chuàng)業(yè)科技有限公司等。

    1.2 CTPK的擴增

    1.2.1總RNA提取與反轉錄按照Trizol[10]抽提總RNA說明書,提取斯氏副柔線蟲總RNA。將獲得的總RNA進行濃度、純度檢測后,利用反轉錄試劑盒合成cDNA。

    1.2.2目的基因的擴增CTPK序列來自本實驗室駱駝斯氏副柔線蟲高通量轉錄組測序序列,應用Primer5.0和Oligo6.24軟件設計引物,同時在引物序列中加入酶切位點EcoRⅠ和XhoⅠ,并由華大基因有限公司合成,預期擴增的DNA片段約為519 bp。引物序列,上游引物:5′-GCAGAATTCATGTCCACAAAGACCAC-3′(下劃線為EcoRⅠ酶切位點);下游引物:5′-CCACTCGAGACCGTTTTGTTGCACG-3′(下劃線為XhoⅠ酶切位點)。以合成的cDNA為模板進行PCR擴增,PCR反應體系50 μL:cDNA 4 μL;Up-Primer 1 μL;Down-Primer 1 μL;Prime Star GXL DNA Polymerase 2 μL;5×Prime Star GXL Buffer(Mg2+plus) 10 μL;dNTP Mixture 4 μL;ddH2O 補齊至 50 μL。PCR擴增條件: 95 ℃預變性5 min; 94 ℃變性30 s, 55 ℃退火30 s, 72 ℃延伸1 min,共31個循環(huán);72 ℃延伸10 min,4 ℃保存。將PCR產(chǎn)物用1%瓊脂糖凝膠電泳進行檢測。

    1.3 CTPK基因的克隆和測序

    將PCR產(chǎn)物目的片段按照Axygen DNA膠回收試劑盒說明書進行膠回收,利用寶生物公司的DNA A-Tailing Kit對CTPK基因PCR產(chǎn)物進行平末端DNA片段的3′末端加一個“A-尾”,從而提高了平末端DNA片段的TA克隆效率。具體步驟:10×A-Tailing Buffer 5 μL;dNTP Mixture 4 μL;A-Tailing Enzyme 0.5 μL;平末端DNA片段0.5~5 μg;補ddH2O至50 μL;緩慢混勻后,72 ℃反應20 min;在冰浴中靜置2 min。然后按照pMD19-T說明書,將加尾后的PCR產(chǎn)物與克隆載體 pMD19-T 16 ℃過夜連接,將連接產(chǎn)物轉化至Trans1-T1感受態(tài)細胞中:將100 μL 感受態(tài)細胞放入冰盒中融化,加入過夜連接產(chǎn)物10 μL,混勻后冰浴30 min,42 ℃水浴45 s,冰浴2 min,然后加入無抗性LB培養(yǎng)液500 μL,37 ℃ 190 r·min-1搖菌2~3 h,最后取100 μL涂布在氨芐青霉素抗性的LB固體培養(yǎng)基上,37 ℃倒置過夜培養(yǎng)。挑取單克隆菌落放入含有氨芐抗性的LB液體培養(yǎng)基中,37 ℃搖菌一段時間后,觀察菌液渾濁程度,如有渾濁可進行菌液PCR鑒定,陽性菌株送華大基因公司進行雙向序列測定。將測序結果正確的質粒命名為pMD19T-CTPK。

    1.4 CTPK基因生物信息學分析

    使用DNAStar中的EditSeq推導斯氏副柔線蟲CTPK基因的氨基酸序列,然后利用Prot-Param在線軟件分析CTPK蛋白的氨基酸序列組成和理化性質;利用ProtScale軟件分析氨基酸殘基的親疏水性;利用TMHMM Server和SignaIP 4.1 Server[11]分析工具預測該蛋白質的跨膜區(qū)及信號肽;利用SOPMA在線軟件分析蛋白質的二級結構;根據(jù)在線軟件Phyre2.0[12]預測CTPK蛋白的三級結構;應用NetPhos3.0 Server軟件預測磷酸化位點;運用SecretomeP軟件預測非典型分泌蛋白;利用DNAStar分析CTPK蛋白質表面可行性,B細胞、T細胞抗原表位等[13]。

    2 結 果

    2.1 目的基因的PCR擴增

    以cDNA為模板擴增CTPK基因,目的片段為519 bp,用1%瓊脂糖凝膠電泳檢測結果與預期目的基因片段大小一致(圖1)。

    M. DNA相對分子質量標準DL2000;1~4. CTPK基因擴增產(chǎn)物M. DNA marker DL2000;1-4. Products of CTPK gene圖1 CTPK基因PCR擴增結果Fig.1 PCR amplification of CTPK gene

    2.2 CTPK基因的生物信息學分析

    2.2.1CTPK基因序列分析與測定利用NCBI中的ORF Finder工具,對CTPK基因cDNA全長序列進行分析,結果表明,該基因cDNA序列全長共519 bp,含有一個編碼173個氨基酸的完整ORF。重組基因測序后序列與轉錄組測序序列相似性為98.46%,發(fā)現(xiàn)8個突變位點:C139T、G165A、C189T、G199A、A213G、G321A、A348T、G408A。利用在線軟件Rare Codon Caltor對斯氏副柔線蟲重組基因測序序列和轉錄組序列進行稀有密碼子分析,發(fā)現(xiàn)是同義突變,且突變前后稀有密碼子種類、數(shù)量沒有改變。

    2.2.2CTPK蛋白理化性質分析使用Port-Param在線軟件分析CTPK蛋白的理化性質,結果表明CTPK蛋白相對分子質量為20.264 ku,理論等電點(theoretical PI)為9.53,原子總數(shù)共2 831個,分子式(Formula)為C895H1403N263O264S6,包含20種氨基酸,其中含量最高的是Lys(占8.7%),含量最低的是Trp占0.6%;帶正電荷數(shù)(Arg+Lys)為27個,帶負電荷數(shù)(Asp+Glu)為16個;蛋白質不穩(wěn)定指數(shù)(instability index)為53.22,說明該蛋白質是不穩(wěn)定蛋白(注:不穩(wěn)定系數(shù)小于40,說明是穩(wěn)定蛋白,反之則是不穩(wěn)定蛋白);脂肪系數(shù)(Aliphatic index)為64.74;總平均親水指數(shù)(Grand average of hydropathicity)為-0.925,說明該蛋白質是親水性蛋白;假定所有的半胱氨酸都形成二硫鍵,CTPK在280 nm處的摩爾消光系數(shù)為23 630 (mol·cm),在0.1%(=1 g·L-1)溶液中的摩爾消光系數(shù)為1.154;而當所有的二硫鍵打開時280 nm處摩爾消光系數(shù)為23 380 (mol·cm),在0.1%(=1 g·L-1)溶液中摩爾消光系數(shù)1.517。當成熟肽N端的1個氨基酸為Met時,CTPK蛋白在體外哺乳動物網(wǎng)狀紅細胞中的半衰期為30 h,在酵母體內的半衰期大于20 h,在大腸埃希菌體內的半衰期大于10 h。

    2.2.3CTPK蛋白疏水性/親水性預測分析氨基酸殘基根據(jù)其親疏水性可分為親水性氨基酸和疏水性氨基酸。親水性氨基酸位于蛋白質的表面,連續(xù)在一起形成抗原表位,同時也是判斷B細胞抗原表位的主要參數(shù)[14]。而疏水性氨基酸位于蛋白內部。應用ProtScale在線軟件對斯氏副柔線蟲CTPK基因編碼蛋白親水性/疏水性進行分析(圖2)。由輸出的文本可知其疏水性在143~144位點有最大值2.222,在27位點有最小值-3.267。根據(jù)氨基酸分值越高疏水性越強或氨基酸分值越低親水性越強的規(guī)律分析得出,CTPK基因編碼蛋白多肽鏈大部分區(qū)域分值為負值,說明該蛋白質是親水性蛋白,與Port-Param預測結果相一致。

    2.2.4CTPK蛋白跨膜區(qū)和信號肽預測膜蛋白是鑲嵌在生物膜中的蛋白質,在細胞中起著至關重要的作用,它們構成了各種激素、受體和神經(jīng)信號,同時也是離子跨膜的通道。應用TMHMM預測CTPK蛋白無跨膜區(qū),不是跨膜蛋白,且整個蛋白全部位于膜外(圖3),跨膜區(qū)中的氨基酸期望值約為0.176(該數(shù)值大于18時表明有跨膜區(qū))。

    圖2 CTPK蛋白疏水性預測Fig.2 Prediction of hydrophobicity of CTPK protein

    信號肽是由蛋白質N端15~30個氨基酸殘基組成的一段疏水肽段,決定蛋白質在細胞中的最終定位。信號肽不但決定蛋白質是否為分泌蛋白,而且與新生肽鏈在細胞內全方位定位有關。應用在線軟件SignalP對CTPK蛋白的信號肽進行預測,綜合圖形和輸出文本分析發(fā)現(xiàn)CTPK不存在信號肽(圖4),但是在運用SecretomeP軟件預測非典型分泌蛋白發(fā)現(xiàn),NN-Score=0.912(若NN-Score>0.5,判定為非典型分泌蛋白),所以該蛋白質屬于分泌蛋白。

    圖3 CTPK 蛋白跨膜結構預測Fig.3 Prediction of CTPK protein transmembarne domain

    圖4 CTPK 蛋白信號肽序列預測Fig.4 Prediction of signal peptide sequence of CTPK protein

    2.2.5CTPK蛋白磷酸化位點的預測磷酸化是最重要的蛋白質翻譯后修飾之一,與蛋白質的功能密切相關,蛋白質的磷酸化與去磷酸化為真核細胞提供了調節(jié)機制。隨著生物信息學的發(fā)展,對磷酸化位點的預測和發(fā)現(xiàn)已成為生物信息學焦點問題之一。通過在線軟件NetPhos2.0 Server預測CTPK蛋白磷酸化位點發(fā)現(xiàn)(圖5),當潛在磷酸化位點的閾值為0.5時,CTPK蛋白存在24個潛在的磷酸化位點,其中包括10個絲氨酸位點,11個蘇氨酸位點,3個酪氨酸位點。

    圖5 CTPK蛋白磷酸化位點預測Fig.5 Prediction of CTPK protein phosphorylation sites

    2.2.6CTPK蛋白二級結構預測通過在線軟件SOPMA對斯氏副柔線蟲CTPK基因編碼蛋白的二級結構進行預測,可見CTPK蛋白中參與α螺旋形成的氨基酸有34個,占19.65%;參與β轉角形成的氨基酸有25個,占14.45%;參與延伸鏈形成的氨基酸有32個,占18.50%;參與無規(guī)則卷曲形成的氨基酸為82個,占47.40%。

    2.2.7CTPK蛋白三級結構預測應用在線軟件 Phyre2.0 對CTPK蛋白的三級結構預測,結果表明:CTPK蛋白由α螺旋、β轉角、延伸鏈和無規(guī)卷曲組合形成超二級結構,經(jīng)折疊、彎曲等一系列復雜的過程形成了CTPK蛋白的三級結構(圖6)。

    圖6 CTPK三級結構的預測Fig.6 Prediction of tertiary structure of CTPK

    2.2.8CTPK蛋白抗原表位的分析利用DNAStar軟件對CTPK蛋白抗原表位進行分析,由圖7可知其親水性位點主要在1—69、72—74、78—81、84—95、99—103、113—138、158—162、170—173位氨基酸處,與之前ProtScale預測的結果相符。根據(jù)Plot-Emini方法預測其表面可行性(圖8),在1—16、19—33、36—60、67—68、82—85、89—91、101—104、111—121、128—137、159—162、171—173位的氨基酸處有較高的可行性。根據(jù)Jameson-Wolf 方法預測CTPK蛋白的抗原性系數(shù),該蛋白質含有較多潛在的B細胞抗原表位(圖9),主要位于1—63、67—73、82—96、98—99、101—107、111—126、133—140、148—151、158—164、171—173位氨基酸殘基或其附近,提示這些區(qū)段含有潛在優(yōu)勢抗原表位。然后再結合蛋白質二級結構、親疏水性和表面可能性,綜合分析CTPK蛋白潛在的蛋白抗原表位,在1—60、82—95、111—126、133—137、158—162、171—173等位氨基酸殘基區(qū)域β轉角和無規(guī)則卷曲結構較多、親水性和表面可能性指數(shù)較高,暴露于表面的概率較大,提示作為B細胞抗原表位的可能性較大。T細胞表位是可以被TCR識別的抗原表位。應用Rothbard-Tayloy方法預測含有特定基序(motif)的潛在T淋巴細胞抗原表位,利用AMPHI方法預測免疫優(yōu)勢輔助性T淋巴細胞抗原位點,綜合兩種方法預測T細胞抗原表位(圖10),CTPK蛋白的T細胞抗原表位可能有9個,主要位于4—10、42—52、75—77、90—98、105—113、115—118、123—128、130—140、158—161位氨基酸殘基或其附近,且分布較均勻。

    3 討 論

    近年來,斯氏副柔線蟲病在駱駝養(yǎng)殖中仍處于較高的感染水平,由于駱駝主要分布在經(jīng)濟欠發(fā)達國家或地區(qū),導致斯氏副柔線蟲病的相關研究資料較少,到目前為止,其確切的生活史仍未完全研究清楚。在斯氏副柔線蟲病的防治上,尚缺乏有效的生前診斷方法和防治新措施。本研究通過克隆CTPK的編碼區(qū)序列,利用生物信息學技術對分泌蛋白CTPK進行分析,采用不同方法和多種參數(shù)聯(lián)合預測,以提高預測的準確性,為駱駝斯氏副柔線蟲病的生前診斷建立和疫苗研制提供參考。

    生物信息學是近年來發(fā)展起來的一種高效的研究手段,促進了生物學研究領域的快速發(fā)展[15],在寄生蟲的研究上,也得到了廣泛應用[16]。生物信息學技術已被廣泛應用于預測和篩選有效的疫苗候選基因[17]。寄生蟲的一些分泌性蛋白是蟲體的表皮、腸上皮或分泌排泄器官不斷分泌的抗原物質。這些抗原暴露于宿主免疫系統(tǒng),是刺激宿主產(chǎn)生特異性免疫應答的重要靶抗原,被認為是與寄生蟲免疫相關的抗原蛋白,在寄生蟲建立寄生生活,免疫調節(jié)和免疫逃避中起關鍵作用[18-19]。這些功能性抗原可以為篩選疫苗抗原、發(fā)掘藥物靶點和確定診斷抗原提供依據(jù)[20-21]。

    圖7 CTPK蛋白親水性預測Fig.7 Prediction of hydrophilicity of CTPK protein

    圖8 CTPK蛋白表面可及性預測結果Fig.8 Prediction of surface probability of CTPK protein

    圖9 CTPK蛋白B細胞抗原表位預測結果Fig.9 Prediction of B-cell antigenic epitopes of CTPK protein

    圖10 CTPK蛋白T細胞抗原表位預測結果Fig.10 Prediction for T-cell antigenic epitopes of CTPK protein

    本試驗克隆了CTPK的編碼區(qū)序列,采用生物信息學方法對分泌蛋白CTPK的理化性質、二級結構、B/T細胞抗原表位進行預測分析,結果表明:蛋白質理論相對分子質量為20.264 ku,理論等電點為9.53,在組成CTPK蛋白的20種氨基酸中,含量最高的是Lys占8.7%,含量最低的是Trp占0.6%;蛋白質不穩(wěn)定指數(shù)為53.22,說明該蛋白質是不穩(wěn)定蛋白;總平均親水指數(shù)為-0.925,說明該蛋白質是親水性蛋白,親水區(qū)的存在有利于蛋白的高效表達;CTPK蛋白是一種不具有信號肽序列的非典型分泌蛋白(NN-Score=0.912),蛋白跨膜區(qū)預測分析表明該蛋白質不是跨膜蛋白,主要位于胞外區(qū)。利用DNAStar對斯氏副柔線蟲CTPK蛋白的T/B細胞抗原表位進行預測,從該蛋白質的親水性與表面可行性等方面分析,該蛋白質親水性多肽較多,其合成的多肽很可能易溶于水;表面可行性較高區(qū)域可能位于CTPK蛋白分子表面。綜合分析發(fā)現(xiàn)在1—60、82—95、111—126、133—137、158—162、171—173位氨基酸殘基存在潛在的B細胞抗原表位;在4—10、42—52、75—77、90—98、105—113、115—118、123—128、130—140、158—161位氨基酸殘基存在潛在的T細胞抗原表位,具有較好的抗原性,符合寄生蟲免疫學診斷和免疫防治候選基因的特點。本試驗為駱駝斯氏副柔線蟲病早期快速診斷試劑盒的建立奠定基礎,同時也為其DNA疫苗的研制提供理論依據(jù)。

    4 結 論

    擴增斯氏副柔線蟲免疫相關基因CTPK的CDS區(qū),其cDNA序列全長共519 bp,編碼173個氨基酸,推導CTPK的相對分子質量為20.26 ku,理論等電點為9.53。結構分析發(fā)現(xiàn),蛋白質無跨膜區(qū),無規(guī)則卷曲構成二級結構的主要組分,親水性氨基酸比例超過60%;抗原表位預測表明,CTPK蛋白是一種抗原性較高的親水性蛋白。推測CTPK有望用作斯氏副柔線蟲的免疫診斷抗原和疫苗候選抗原。

    參考文獻(References):

    [1]盧俊杰, 靳家聲. 人和動物寄生線蟲圖譜[M]. 北京: 中國農(nóng)業(yè)科學技術出版社, 2002: 405-459.

    LU J J, JIN J S. The map of human and animal parasitic nematodes[M]. Beijing: China Agricultural Science and Technology Press, 2002: 405-459. (in Chinese)

    [2]沈杰, 黃兵. 中國家畜家禽寄生蟲名錄[M]. 北京: 中國農(nóng)業(yè)科學技術出版社, 2004: 105-111.

    SHEN J, HUANG B. A list of parasites for livestock and poultry in China[M]. Beijing: China Agricultural Science and Technology Press, 2004: 105-111. (in Chinese)

    [3]HASHEMINASAB S S. Molecular characterization of the first internal transcribed spacer of rDNA ofParabronemaskrjabinifor the first time in sheep[J].AnnParasitol, 2015, 61(4): 241-246.

    [4]HASHEMINASAB S S, JALOUSIAN F, MESHGI B. Molecular and morphological characterization ofParabronemaskrjabiniof sheep and goats at three different geographical zones in Iran[J].AnnParasitol, 2016, 62(1): 55-61.

    [5]楊曉野. 內蒙古自治區(qū)駱駝寄生蟲種類調查[D]. 呼和浩特: 內蒙古農(nóng)牧學院, 1985: 40-41.

    YANG X Y. Camel parasitic species investigation in Inner Mongolia Autonomous Region[D]. Hohhot: Inner Mongolia Agricultural University, 1985: 40-41. (in Chinese)

    [6]楊蓮茹, 楊曉野, 劉珍蓮, 等. 內蒙古地區(qū)駱駝斯氏副柔線蟲(Parabronemaskrjabini)病調查[J]. 內蒙古農(nóng)業(yè)大學學報, 2004, 25(1): 43-45.

    YANG L R, YANG X Y, LIU Z L, et al. Investigation ofParabronemaSkrjabinidisease of camels in Inner Mongolia Region[J].JournalofInnerMongoliaAgriculturalUniversity, 2004, 25(1): 43-45. (in Chinese)

    [7]趙治國. 我國駱駝斯氏副柔線蟲病傳播媒介的研究[D]. 呼和浩特: 內蒙古農(nóng)業(yè)大學, 2010.

    ZHAO Z G. Study on the vector of camelParabronemosisin China[D]. Hohhot: Inner Mongolia Agricultural University, 2010. (in Chinese)

    [8]張曉東, 楊曉野, 楊蓮茹, 等. 斯氏副柔線蟲rDNA-ITS片段的克隆及序列分析[J]. 畜牧獸醫(yī)學報, 2009, 40(5): 775-779.

    ZHANG X D, YANG X Y, YANG L R, et al. Cloning and sequence analysis of the rDNA-ITS ofParabronemaskrjabini[J].ActaVeterinariaetZootechnicaSinica, 2009, 40(5): 775-779. (in Chinese)

    [9]馮陳晨, 王文龍, 呼和巴特爾. 駱駝斯氏副柔線蟲轉錄組的高通量測序及分析[J]. 中國獸醫(yī)學報, 2017, 37(4): 671-675.

    FENG C C, WANG W L, HUHR B T E. Sequencing and analysis of the transcriptome ofParabronemaskrjabiniin camle[J].ChineseJournalofVeterinaryScience, 2017, 37(4): 671-675. (in Chinese)

    [11]PETERSEN T N, BRUNAK S, VON HEIJNE G, et al. SignalP 4.0: discriminating signal peptides from transmembrane regions[J].NatMethods, 2011, 8(10): 785-786.

    [12]KELLEY L A, MEZULIS S, YATES C M, et al. The Phyre2 web portal for protein modeling, prediction and analysis[J].NatProtoc, 2015, 10(6): 845-858.

    [13]李玲玲, 朱珊麗, 李文姝, 等. EB病毒核蛋白-1 B細胞表位的預測及其同源性分析[J]. 生物醫(yī)學工程學雜志, 2011, 28(2): 371-375.

    LI L L, ZHU S L, LI W S, et al. Prediction and research on homology of B-cell epitopes of Epstein-Barr virus nuclear antigen-1[J].JournalofBiomedicalEngineering, 2011, 28(2): 371-375. (in Chinese)

    [14]SIVALINGAM G N, SHEPHERD A J. An analysis of B-cell epitope discontinuity[J].MolImmunol, 2012, 51(3-4): 304-309.

    [15]BLAZIE S M, BABB C, WILKY H, et al. Comparative RNA-Seq analysis reveals pervasive tissue-specific alternative polyadenylation inCaenorhabditiselegansintestine and muscles[J].BMCBiol, 2015, 13: 4.

    [16]LIU G H, XU M J, CHANG Q C, et al.Denovotranscriptomic analysis of the female and male adults of the blood flukeSchistosomaturkestanicum[J].ParasitVectors, 2016, 9: 143.

    [17]成文敏, 潘偉榮, 尹俊芳, 等. 版納微型豬近交系DLDH基因編碼區(qū)序列擴增及生物信息學分析[J]. 畜牧獸醫(yī)學報, 2015, 46(10): 1706-1712.

    CHENG W M, PAN W R, YIN J F, et al. Amplification and bioinformatics analysis of coding region of theDLDHgene in Banna Mini-pig inbred line (BMI)[J].ActaVeterinariaetZootechnicaSinica, 2015, 46(10): 1706-1712. (in Chinese)

    [18]TRITTEN L, CLARKE D, TIMMINS S, et al.Dirofilariaimmitisexhibits sex-and stage-specific differences in excretory/secretory miRNA and protein profiles[J].VetParasitol, 2016, 232: 1-7.

    [19]EBERLE R, BRATTIG N W, TRUSCH M, et al. Isolation, identification and functional profile of excretory-secretory peptides fromOnchocercaochengi[J].ActaTrop, 2015, 142: 156-166.

    [20]CAO X D, FU Z Q, ZHANG M, et al. iTRAQ-based comparative proteomic analysis of excretory-secretory proteins of schistosomula and adult worms ofSchistosomajaponicum[J].JProteomics, 2016, 138: 30-39.

    [21]CAO X D, FU Z Q, ZHANG M, et al. Excretory/secretory proteome of 14-day schistosomula,Schistosomajaponicum[J].JProteomics, 2016, 130: 221-230.

    猜你喜歡
    斯氏表位親水性
    H3N2流感病毒HA保守Th表位對CD4+T細胞活化及分化的影響
    雙負載抗生素親水性聚氨酯泡沫的制備與表征
    中蜂斯氏蜜蜂繭蜂的生物防治
    斯氏鈍綏螨對西花薊馬的捕食功能反應
    聯(lián)合T、B細胞表位設計多肽疫苗的研究進展①
    空氣中納秒脈沖均勻DBD增加聚合物的表面親水性
    小反芻獸疫病毒化學合成表位多肽對小鼠的免疫效果研究
    結核分枝桿菌抗原Lppx和MT0322人T細胞抗原表位的多態(tài)性研究
    關于“斯氏體系”的系統(tǒng)化研究
    戲劇之家(2015年2期)2015-03-13 00:27:27
    水刺型空氣加濕器濾材的親水性改性研究
    欧美日韩综合久久久久久 | 国产伦在线观看视频一区| 嫁个100分男人电影在线观看| 国产精品一及| 国产精品久久久久久久电影 | 国产精品98久久久久久宅男小说| 国产综合懂色| 美女 人体艺术 gogo| 国产单亲对白刺激| 欧美中文综合在线视频| 校园春色视频在线观看| 日韩欧美精品v在线| 精品久久久久久久久久久久久| 亚洲人成电影免费在线| 香蕉丝袜av| 国产高清激情床上av| 国产乱人视频| 欧美日韩精品网址| 日日摸夜夜添夜夜添小说| a级毛片a级免费在线| 国产免费av片在线观看野外av| 搡老岳熟女国产| 久久精品人妻少妇| 91久久精品电影网| 激情在线观看视频在线高清| 亚洲久久久久久中文字幕| 在线看三级毛片| 99riav亚洲国产免费| 丰满乱子伦码专区| 日本 av在线| 淫妇啪啪啪对白视频| 精品免费久久久久久久清纯| 成人三级黄色视频| 日本黄色视频三级网站网址| 舔av片在线| 神马国产精品三级电影在线观看| 欧美日韩精品网址| 99热这里只有是精品50| 99热6这里只有精品| 丰满乱子伦码专区| 婷婷六月久久综合丁香| 男人舔奶头视频| 淫妇啪啪啪对白视频| 成人av在线播放网站| 观看美女的网站| 女人高潮潮喷娇喘18禁视频| 国产探花极品一区二区| 国产在视频线在精品| 国产成人影院久久av| 国产一区二区在线av高清观看| 97人妻精品一区二区三区麻豆| 悠悠久久av| 亚洲av五月六月丁香网| 久久久久久久久中文| 熟女少妇亚洲综合色aaa.| svipshipincom国产片| 成人亚洲精品av一区二区| 欧洲精品卡2卡3卡4卡5卡区| 欧美中文综合在线视频| 蜜桃久久精品国产亚洲av| 国产美女午夜福利| 亚洲五月天丁香| or卡值多少钱| 国产欧美日韩精品亚洲av| 少妇丰满av| 亚洲一区二区三区色噜噜| 无限看片的www在线观看| 亚洲第一电影网av| 三级毛片av免费| 亚洲成av人片免费观看| 日本a在线网址| 热99在线观看视频| 一本综合久久免费| 夜夜躁狠狠躁天天躁| 免费人成在线观看视频色| 中文字幕av在线有码专区| 国产精品av视频在线免费观看| 久久亚洲真实| 亚洲电影在线观看av| 特大巨黑吊av在线直播| 99国产精品一区二区蜜桃av| 欧美日韩瑟瑟在线播放| 黄片大片在线免费观看| 国产高清激情床上av| 特级一级黄色大片| 午夜福利免费观看在线| 天堂av国产一区二区熟女人妻| 老汉色av国产亚洲站长工具| 麻豆一二三区av精品| 久久九九热精品免费| 美女大奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 老鸭窝网址在线观看| 高清日韩中文字幕在线| 国产精品亚洲一级av第二区| 国产高清videossex| 国产成人欧美在线观看| 1024手机看黄色片| 一夜夜www| 97超视频在线观看视频| 2021天堂中文幕一二区在线观| 国产精品久久久久久久电影 | 岛国在线观看网站| 最近最新中文字幕大全免费视频| 精品人妻偷拍中文字幕| 亚洲久久久久久中文字幕| 国产一区二区三区在线臀色熟女| 国内少妇人妻偷人精品xxx网站| 91在线观看av| aaaaa片日本免费| 亚洲国产精品成人综合色| 2021天堂中文幕一二区在线观| 欧美+日韩+精品| 又爽又黄无遮挡网站| 欧美日韩乱码在线| 黄色女人牲交| 91九色精品人成在线观看| 久久久国产成人精品二区| 成人午夜高清在线视频| 婷婷精品国产亚洲av在线| 国产伦精品一区二区三区四那| 午夜激情福利司机影院| 亚洲国产欧美人成| 全区人妻精品视频| 亚洲不卡免费看| 美女cb高潮喷水在线观看| 草草在线视频免费看| 手机成人av网站| 亚洲不卡免费看| 欧美日韩福利视频一区二区| 搡老妇女老女人老熟妇| 久久中文看片网| 午夜福利在线观看免费完整高清在 | 久久精品国产综合久久久| 亚洲成av人片免费观看| 久久久国产精品麻豆| 搡老岳熟女国产| 丁香六月欧美| 欧美黄色淫秽网站| 亚洲av熟女| 黄片小视频在线播放| 91久久精品电影网| 老司机在亚洲福利影院| 欧美乱妇无乱码| 国产野战对白在线观看| 99热精品在线国产| x7x7x7水蜜桃| 女人十人毛片免费观看3o分钟| 国产高清视频在线观看网站| 国产主播在线观看一区二区| 色在线成人网| 亚洲 欧美 日韩 在线 免费| 一本精品99久久精品77| xxxwww97欧美| 一本久久中文字幕| 国产精品 国内视频| 偷拍熟女少妇极品色| 欧美黄色片欧美黄色片| 久久国产乱子伦精品免费另类| 身体一侧抽搐| 看黄色毛片网站| 高清在线国产一区| 三级毛片av免费| 国产精华一区二区三区| 亚洲熟妇熟女久久| 国产三级中文精品| 99国产精品一区二区三区| 天堂影院成人在线观看| 男女之事视频高清在线观看| 男女床上黄色一级片免费看| 亚洲av五月六月丁香网| 免费一级毛片在线播放高清视频| 18禁裸乳无遮挡免费网站照片| 国产一区二区亚洲精品在线观看| 亚洲成人精品中文字幕电影| 婷婷亚洲欧美| 熟女少妇亚洲综合色aaa.| 亚洲五月婷婷丁香| 亚洲狠狠婷婷综合久久图片| 99久久精品热视频| 日本免费一区二区三区高清不卡| 免费av不卡在线播放| 美女高潮的动态| 一个人免费在线观看的高清视频| 国产三级黄色录像| 精品久久久久久久毛片微露脸| 99国产综合亚洲精品| 精品免费久久久久久久清纯| 一级黄片播放器| 99久久精品一区二区三区| 99视频精品全部免费 在线| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看成人毛片| 美女 人体艺术 gogo| 一进一出抽搐动态| 色视频www国产| 欧美日韩亚洲国产一区二区在线观看| 男女那种视频在线观看| 麻豆一二三区av精品| 少妇高潮的动态图| 每晚都被弄得嗷嗷叫到高潮| 成人三级黄色视频| 五月伊人婷婷丁香| 国产免费av片在线观看野外av| АⅤ资源中文在线天堂| 叶爱在线成人免费视频播放| 色播亚洲综合网| 午夜激情欧美在线| 欧美另类亚洲清纯唯美| 99热精品在线国产| 亚洲国产精品合色在线| 国产一级毛片七仙女欲春2| 久久精品人妻少妇| 在线视频色国产色| 亚洲,欧美精品.| 一二三四社区在线视频社区8| 欧美日韩乱码在线| 日韩高清综合在线| 午夜激情福利司机影院| 别揉我奶头~嗯~啊~动态视频| 禁无遮挡网站| 桃色一区二区三区在线观看| 日韩欧美在线二视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 身体一侧抽搐| 国产野战对白在线观看| 国产av在哪里看| 久久精品夜夜夜夜夜久久蜜豆| 美女高潮喷水抽搐中文字幕| 老汉色av国产亚洲站长工具| 一个人免费在线观看电影| 国产成人福利小说| 久久久久久国产a免费观看| 在线观看美女被高潮喷水网站 | 国产老妇女一区| 免费在线观看影片大全网站| 国内少妇人妻偷人精品xxx网站| 国产蜜桃级精品一区二区三区| 亚洲av免费在线观看| 国产伦人伦偷精品视频| 久久性视频一级片| 亚洲人成伊人成综合网2020| 亚洲国产精品久久男人天堂| 亚洲av中文字字幕乱码综合| 亚洲国产精品合色在线| 变态另类成人亚洲欧美熟女| 日本成人三级电影网站| 99久久久亚洲精品蜜臀av| 一级a爱片免费观看的视频| 精品久久久久久久毛片微露脸| 最后的刺客免费高清国语| 三级毛片av免费| 国产av一区在线观看免费| 亚洲人成网站在线播| 草草在线视频免费看| 天天添夜夜摸| 狂野欧美白嫩少妇大欣赏| 麻豆一二三区av精品| 亚洲av电影不卡..在线观看| 中文字幕人妻熟人妻熟丝袜美 | 九九久久精品国产亚洲av麻豆| 久久精品夜夜夜夜夜久久蜜豆| 色噜噜av男人的天堂激情| 免费电影在线观看免费观看| 国产一区二区三区在线臀色熟女| 国产高清有码在线观看视频| 麻豆国产av国片精品| 欧洲精品卡2卡3卡4卡5卡区| 真人做人爱边吃奶动态| 午夜两性在线视频| 亚洲自拍偷在线| 亚洲精品美女久久久久99蜜臀| 欧美zozozo另类| 亚洲欧美激情综合另类| 最近最新中文字幕大全电影3| 欧美绝顶高潮抽搐喷水| 国产高清视频在线播放一区| av国产免费在线观看| 女生性感内裤真人,穿戴方法视频| 在线观看午夜福利视频| 欧美丝袜亚洲另类 | 村上凉子中文字幕在线| 国产av不卡久久| 18美女黄网站色大片免费观看| 美女高潮的动态| 99精品在免费线老司机午夜| 黄色日韩在线| 国产三级在线视频| 国产一区二区亚洲精品在线观看| 亚洲av免费高清在线观看| 久久6这里有精品| 日本免费a在线| 国产麻豆成人av免费视频| 亚洲avbb在线观看| 少妇的逼水好多| 久久国产精品影院| 国产不卡一卡二| 色吧在线观看| 欧美黑人巨大hd| 在线观看免费午夜福利视频| 观看免费一级毛片| 国产精品乱码一区二三区的特点| 大型黄色视频在线免费观看| 国产精品av视频在线免费观看| 观看美女的网站| 很黄的视频免费| 两个人看的免费小视频| 国产乱人伦免费视频| 久久国产乱子伦精品免费另类| 国产国拍精品亚洲av在线观看 | 亚洲人成伊人成综合网2020| 两个人看的免费小视频| 亚洲色图av天堂| av天堂在线播放| 久久精品综合一区二区三区| 国产成人影院久久av| 午夜福利成人在线免费观看| 久久久久免费精品人妻一区二区| 老熟妇仑乱视频hdxx| 丰满人妻一区二区三区视频av | 久久久国产成人免费| 国产 一区 欧美 日韩| 免费无遮挡裸体视频| 最好的美女福利视频网| 国产亚洲欧美在线一区二区| 欧美日韩乱码在线| 免费看十八禁软件| 老司机福利观看| 亚洲国产欧美网| 精品福利观看| 99久久久亚洲精品蜜臀av| 久久亚洲精品不卡| 亚洲美女视频黄频| 老鸭窝网址在线观看| 欧美极品一区二区三区四区| 欧美绝顶高潮抽搐喷水| 在线天堂最新版资源| 亚洲欧美日韩无卡精品| 九九热线精品视视频播放| 亚洲七黄色美女视频| 精品无人区乱码1区二区| 精品久久久久久久久久免费视频| 三级男女做爰猛烈吃奶摸视频| av欧美777| 久久6这里有精品| 特大巨黑吊av在线直播| 国产精品99久久久久久久久| ponron亚洲| 亚洲av美国av| 俺也久久电影网| 亚洲不卡免费看| 免费av不卡在线播放| 国产野战对白在线观看| 亚洲性夜色夜夜综合| 久久午夜亚洲精品久久| 国产成人a区在线观看| 色哟哟哟哟哟哟| 少妇的逼好多水| 黄色片一级片一级黄色片| 熟妇人妻久久中文字幕3abv| 日本五十路高清| 国产精品三级大全| 午夜精品在线福利| 夜夜躁狠狠躁天天躁| 两个人的视频大全免费| 欧美一区二区国产精品久久精品| 好男人在线观看高清免费视频| 国产在线精品亚洲第一网站| 欧美日韩精品网址| av国产免费在线观看| 岛国在线免费视频观看| 日本三级黄在线观看| 日本成人三级电影网站| 美女高潮喷水抽搐中文字幕| 熟女人妻精品中文字幕| 搡老熟女国产l中国老女人| 久久欧美精品欧美久久欧美| 在线观看美女被高潮喷水网站 | 91麻豆精品激情在线观看国产| 黄片大片在线免费观看| 怎么达到女性高潮| 身体一侧抽搐| 成年女人永久免费观看视频| 日本撒尿小便嘘嘘汇集6| 韩国av一区二区三区四区| 国产精品香港三级国产av潘金莲| 欧美xxxx黑人xx丫x性爽| 国产av一区在线观看免费| 国产亚洲精品久久久久久毛片| 男女床上黄色一级片免费看| 国产美女午夜福利| 搡老熟女国产l中国老女人| 欧美av亚洲av综合av国产av| 色播亚洲综合网| 中文字幕人妻熟人妻熟丝袜美 | 亚洲aⅴ乱码一区二区在线播放| 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 久久久成人免费电影| 亚洲成av人片免费观看| 波多野结衣巨乳人妻| 最好的美女福利视频网| 午夜亚洲福利在线播放| 欧洲精品卡2卡3卡4卡5卡区| 亚洲 欧美 日韩 在线 免费| 噜噜噜噜噜久久久久久91| 色综合欧美亚洲国产小说| 黄色成人免费大全| 亚洲国产欧美人成| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲美女久久久| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费激情av| 久久久久久久久中文| 亚洲国产精品sss在线观看| 色播亚洲综合网| 一级毛片女人18水好多| 亚洲av成人不卡在线观看播放网| 一个人免费在线观看电影| 级片在线观看| 亚洲专区中文字幕在线| 最近最新中文字幕大全免费视频| 级片在线观看| 一区二区三区高清视频在线| 18禁在线播放成人免费| 欧美国产日韩亚洲一区| 动漫黄色视频在线观看| 露出奶头的视频| 国产精品久久久久久久久免 | www国产在线视频色| 国产视频一区二区在线看| 亚洲精品国产精品久久久不卡| 两人在一起打扑克的视频| 日本熟妇午夜| 色吧在线观看| 精品无人区乱码1区二区| 欧美日韩精品网址| 99精品欧美一区二区三区四区| 两人在一起打扑克的视频| 人人妻人人看人人澡| 日韩欧美国产一区二区入口| 国产伦人伦偷精品视频| 国内久久婷婷六月综合欲色啪| 亚洲精品一区av在线观看| 嫩草影院精品99| 亚洲天堂国产精品一区在线| 日韩欧美 国产精品| 深爱激情五月婷婷| 热99re8久久精品国产| 成人高潮视频无遮挡免费网站| avwww免费| 亚洲人成电影免费在线| 国产毛片a区久久久久| 99在线人妻在线中文字幕| 欧美3d第一页| 精品国产超薄肉色丝袜足j| 色综合站精品国产| 国产精品久久久久久久久免 | 国产成人啪精品午夜网站| 午夜免费成人在线视频| 天美传媒精品一区二区| 亚洲在线自拍视频| 国产成人系列免费观看| 毛片女人毛片| 欧美激情久久久久久爽电影| 国产成年人精品一区二区| 51国产日韩欧美| 女同久久另类99精品国产91| 欧美日韩中文字幕国产精品一区二区三区| www国产在线视频色| 人妻夜夜爽99麻豆av| 免费人成视频x8x8入口观看| 日日夜夜操网爽| 国产高清激情床上av| 国产中年淑女户外野战色| 亚洲精品一卡2卡三卡4卡5卡| 国产三级黄色录像| 一卡2卡三卡四卡精品乱码亚洲| 精品一区二区三区视频在线观看免费| 日韩欧美免费精品| 最近最新中文字幕大全电影3| 亚洲精品在线观看二区| a在线观看视频网站| 国产精品99久久久久久久久| 美女高潮喷水抽搐中文字幕| 日本黄色片子视频| 一进一出抽搐gif免费好疼| 日韩高清综合在线| 亚洲国产日韩欧美精品在线观看 | 国产成人欧美在线观看| 国产一区二区三区视频了| 久久6这里有精品| 国产午夜精品论理片| 亚洲国产欧洲综合997久久,| 亚洲国产中文字幕在线视频| av中文乱码字幕在线| 好男人在线观看高清免费视频| 欧美日韩瑟瑟在线播放| 一个人免费在线观看的高清视频| 国产综合懂色| 日韩 欧美 亚洲 中文字幕| 欧美色欧美亚洲另类二区| 老汉色∧v一级毛片| 日本 欧美在线| 免费看光身美女| 国产精品精品国产色婷婷| 啦啦啦韩国在线观看视频| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 日本成人三级电影网站| 超碰av人人做人人爽久久 | 欧美三级亚洲精品| 男女视频在线观看网站免费| 国产午夜精品论理片| or卡值多少钱| 美女被艹到高潮喷水动态| 欧美黄色片欧美黄色片| 精品久久久久久成人av| 色在线成人网| 国产视频内射| 久久久久久人人人人人| 性色avwww在线观看| 国产精品国产高清国产av| 亚洲狠狠婷婷综合久久图片| 日日夜夜操网爽| 88av欧美| 在线播放国产精品三级| 啦啦啦免费观看视频1| 免费在线观看亚洲国产| 久久久久久人人人人人| 亚洲最大成人手机在线| 精品福利观看| 亚洲欧美一区二区三区黑人| 国产亚洲精品综合一区在线观看| 天堂动漫精品| 精品一区二区三区av网在线观看| 久久中文看片网| 黄色丝袜av网址大全| av欧美777| 色综合欧美亚洲国产小说| 2021天堂中文幕一二区在线观| 国产成人av教育| 91久久精品电影网| 日韩亚洲欧美综合| 在线观看免费视频日本深夜| 母亲3免费完整高清在线观看| 天天一区二区日本电影三级| 亚洲av一区综合| 在线观看舔阴道视频| 久久久久久大精品| 欧美高清成人免费视频www| 久久久国产精品麻豆| 久久久精品大字幕| 亚洲av免费在线观看| 亚洲第一电影网av| 日本 av在线| 午夜福利成人在线免费观看| 美女免费视频网站| 精品无人区乱码1区二区| 老熟妇仑乱视频hdxx| 精品熟女少妇八av免费久了| 又紧又爽又黄一区二区| 天天躁日日操中文字幕| 免费看光身美女| 制服丝袜大香蕉在线| 久久精品国产综合久久久| 中文字幕熟女人妻在线| 成人精品一区二区免费| 国产三级黄色录像| 99精品在免费线老司机午夜| 中文字幕人妻丝袜一区二区| xxx96com| 黄片小视频在线播放| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三| 色综合欧美亚洲国产小说| 成人三级黄色视频| 又爽又黄无遮挡网站| 欧美av亚洲av综合av国产av| 狂野欧美激情性xxxx| 精品久久久久久,| 成年免费大片在线观看| 亚洲真实伦在线观看| 97超级碰碰碰精品色视频在线观看| 日韩欧美国产一区二区入口| 18禁黄网站禁片午夜丰满| ponron亚洲| 在线免费观看不下载黄p国产 | 日本一二三区视频观看| 午夜福利在线观看免费完整高清在 | 国产亚洲欧美98| av中文乱码字幕在线| 成人鲁丝片一二三区免费| 精华霜和精华液先用哪个| 桃色一区二区三区在线观看| 日本与韩国留学比较| 人妻久久中文字幕网| 国产主播在线观看一区二区| 日本三级黄在线观看| 国产精品免费一区二区三区在线| 欧美国产日韩亚洲一区| 一区福利在线观看| 99热6这里只有精品| 99视频精品全部免费 在线| 噜噜噜噜噜久久久久久91| 日本在线视频免费播放| 日本五十路高清| a在线观看视频网站| 网址你懂的国产日韩在线| 欧美性感艳星| 天堂影院成人在线观看| 免费看a级黄色片| 亚洲欧美激情综合另类| 无遮挡黄片免费观看| 亚洲男人的天堂狠狠|