• <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
    水刺型空氣加濕器濾材的親水性改性研究
    国产深夜福利视频在线观看| 中文字幕色久视频| 亚洲久久久国产精品| 久久精品aⅴ一区二区三区四区| 国产蜜桃级精品一区二区三区 | 国产成人系列免费观看| 中文字幕精品免费在线观看视频| 亚洲国产精品合色在线| 又黄又粗又硬又大视频| 欧美日韩av久久| 欧美性长视频在线观看| 九色亚洲精品在线播放| 亚洲熟女毛片儿| 伊人久久大香线蕉亚洲五| 久久婷婷成人综合色麻豆| 一进一出抽搐gif免费好疼 | 不卡一级毛片| 中文字幕制服av| 91字幕亚洲| videos熟女内射| 免费看a级黄色片| 精品国产一区二区三区四区第35| 黄色女人牲交| 久久香蕉国产精品| 欧美日韩乱码在线| 男女下面插进去视频免费观看| 超色免费av| 国产极品粉嫩免费观看在线| 欧美人与性动交α欧美精品济南到| 99精品久久久久人妻精品| 国产精品美女特级片免费视频播放器 | 亚洲av美国av| 五月开心婷婷网| 亚洲国产欧美日韩在线播放| 又黄又爽又免费观看的视频| 色婷婷久久久亚洲欧美| 在线国产一区二区在线| 美国免费a级毛片| 欧美激情久久久久久爽电影 | 日韩大码丰满熟妇| 午夜精品在线福利| 两人在一起打扑克的视频| 国产精品免费大片| 人妻丰满熟妇av一区二区三区 | 色在线成人网| 他把我摸到了高潮在线观看| 亚洲精品国产精品久久久不卡| 叶爱在线成人免费视频播放| 国产精品免费一区二区三区在线 | 一夜夜www| 嫁个100分男人电影在线观看| 久久性视频一级片| 久99久视频精品免费| 欧美黄色淫秽网站| 精品午夜福利视频在线观看一区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲情色 制服丝袜| 亚洲av片天天在线观看| 久热爱精品视频在线9| 久久久久久人人人人人| 天天躁狠狠躁夜夜躁狠狠躁| 两个人免费观看高清视频| 在线免费观看的www视频| 亚洲av电影在线进入| 亚洲av第一区精品v没综合| 麻豆av在线久日| 成年人黄色毛片网站| 亚洲免费av在线视频| 久久久久国产一级毛片高清牌| 91大片在线观看| 午夜日韩欧美国产| 亚洲精品成人av观看孕妇| 免费看a级黄色片| 精品乱码久久久久久99久播| 女性生殖器流出的白浆| 99国产精品一区二区蜜桃av | 亚洲va日本ⅴa欧美va伊人久久| 亚洲专区国产一区二区| 大型黄色视频在线免费观看| 日韩欧美一区视频在线观看| 天堂√8在线中文| 亚洲精品中文字幕一二三四区| 午夜福利在线观看吧| 精品久久久久久电影网| 国产精品永久免费网站| 窝窝影院91人妻| 欧美精品啪啪一区二区三区| 激情视频va一区二区三区| 十八禁高潮呻吟视频| 热99久久久久精品小说推荐| 亚洲av电影在线进入| 桃红色精品国产亚洲av| 国产精华一区二区三区| 国产成人欧美在线观看 | 美女高潮喷水抽搐中文字幕| 亚洲九九香蕉| 久久精品亚洲av国产电影网| 成人国语在线视频| 久久久久久久久免费视频了| 色老头精品视频在线观看| 国产精品综合久久久久久久免费 | 桃红色精品国产亚洲av| av片东京热男人的天堂| 国产精品久久电影中文字幕 | 自线自在国产av| 国产xxxxx性猛交| 亚洲久久久国产精品| 天天操日日干夜夜撸| 精品国产一区二区久久| 国产高清videossex| 又黄又粗又硬又大视频| 欧美日韩黄片免| 午夜福利一区二区在线看| 午夜福利在线免费观看网站| 母亲3免费完整高清在线观看| 黄色片一级片一级黄色片| 校园春色视频在线观看| 美女福利国产在线| 久久香蕉精品热| 国产真人三级小视频在线观看| 啦啦啦 在线观看视频| 夜夜爽天天搞| 精品无人区乱码1区二区| 国产人伦9x9x在线观看| 热99国产精品久久久久久7| 亚洲av成人不卡在线观看播放网| 中文字幕人妻丝袜一区二区| 午夜福利在线观看吧| 最新在线观看一区二区三区| 日日夜夜操网爽| 下体分泌物呈黄色| 久久人人97超碰香蕉20202| 欧美在线黄色| 国产成人精品在线电影| 嫁个100分男人电影在线观看| 狂野欧美激情性xxxx| 精品人妻在线不人妻| 国产精品香港三级国产av潘金莲| 亚洲中文av在线| 视频在线观看一区二区三区| 热99re8久久精品国产| 亚洲情色 制服丝袜| www日本在线高清视频| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩一区二区三| 91av网站免费观看| 成人手机av| 777米奇影视久久| 亚洲国产欧美网| 999久久久精品免费观看国产| 久久久国产欧美日韩av| 午夜福利,免费看| 欧美精品人与动牲交sv欧美| 国产精品美女特级片免费视频播放器 | 国产成人精品在线电影| 久久精品91无色码中文字幕| 精品一区二区三区四区五区乱码| 丰满迷人的少妇在线观看| 色老头精品视频在线观看| 久久久久久人人人人人| 国产又爽黄色视频| 欧美激情极品国产一区二区三区| 99riav亚洲国产免费| 亚洲国产中文字幕在线视频| 女人爽到高潮嗷嗷叫在线视频| 动漫黄色视频在线观看| 精品久久久久久,| 人人妻人人澡人人爽人人夜夜| 99riav亚洲国产免费| 一级片'在线观看视频| 亚洲国产欧美网| videos熟女内射| 露出奶头的视频| 免费看十八禁软件| 欧美日韩亚洲国产一区二区在线观看 | 一级毛片女人18水好多| 成年人午夜在线观看视频| 在线视频色国产色| 女人久久www免费人成看片| 久久性视频一级片| 成人影院久久| 亚洲五月色婷婷综合| avwww免费| 亚洲av熟女| 精品午夜福利视频在线观看一区| 好看av亚洲va欧美ⅴa在| 亚洲专区中文字幕在线| a级片在线免费高清观看视频| 免费女性裸体啪啪无遮挡网站| 久99久视频精品免费| 精品少妇久久久久久888优播| 一本大道久久a久久精品| 在线观看免费视频日本深夜| 国产又爽黄色视频| 51午夜福利影视在线观看| 精品无人区乱码1区二区| 午夜91福利影院| 午夜精品国产一区二区电影| 亚洲av成人av| 中文欧美无线码| 搡老乐熟女国产| 又紧又爽又黄一区二区| 一区福利在线观看| 在线看a的网站| 99热网站在线观看| 韩国av一区二区三区四区| 韩国精品一区二区三区| 亚洲精品粉嫩美女一区| 香蕉国产在线看| 制服人妻中文乱码| 在线观看66精品国产| av欧美777| 电影成人av| www.精华液| 岛国在线观看网站| videos熟女内射| 天天添夜夜摸| x7x7x7水蜜桃| 窝窝影院91人妻| 丝袜在线中文字幕| 亚洲熟妇中文字幕五十中出 | 亚洲色图 男人天堂 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 又黄又爽又免费观看的视频| 欧美在线黄色| 亚洲国产精品一区二区三区在线| 久久久国产成人精品二区 | 国产精品免费视频内射| 日韩欧美一区视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久久精品吃奶| 女人精品久久久久毛片| 精品高清国产在线一区| 亚洲伊人色综图| 日韩视频一区二区在线观看| 精品国产一区二区久久| 免费一级毛片在线播放高清视频 | 成熟少妇高潮喷水视频| 色播在线永久视频| 日韩制服丝袜自拍偷拍| 视频区欧美日本亚洲| 国产精品美女特级片免费视频播放器 | 亚洲精品av麻豆狂野| 悠悠久久av| 在线天堂中文资源库| 午夜老司机福利片| 男男h啪啪无遮挡| 看黄色毛片网站| 成人av一区二区三区在线看| 久久久久久久午夜电影 | 国产成人免费观看mmmm| 一本综合久久免费| 不卡一级毛片| 91成年电影在线观看| 国产精品美女特级片免费视频播放器 | 高潮久久久久久久久久久不卡| 亚洲精品在线美女| 高清视频免费观看一区二区| 99精品久久久久人妻精品| 久久影院123| 国产亚洲精品久久久久久毛片 | 91老司机精品| 欧美+亚洲+日韩+国产| 久久国产精品大桥未久av| 欧美精品一区二区免费开放| 国内久久婷婷六月综合欲色啪| 国产精品99久久99久久久不卡| 丰满的人妻完整版| 欧美精品人与动牲交sv欧美| 欧美成人免费av一区二区三区 | 午夜福利一区二区在线看| 国产成人影院久久av| 午夜免费鲁丝| www.自偷自拍.com| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 中文字幕av电影在线播放| 村上凉子中文字幕在线| 中国美女看黄片| 日日夜夜操网爽| 国产男靠女视频免费网站| 母亲3免费完整高清在线观看| xxxhd国产人妻xxx| 极品少妇高潮喷水抽搐| 久久天躁狠狠躁夜夜2o2o| 国产99久久九九免费精品| 在线观看一区二区三区激情| 中文欧美无线码| 国产亚洲精品一区二区www | 麻豆av在线久日| xxx96com| 亚洲国产欧美日韩在线播放| 视频区欧美日本亚洲| 777米奇影视久久| 久久久精品区二区三区| 国产精品九九99| 亚洲成a人片在线一区二区| 在线十欧美十亚洲十日本专区| 国产欧美日韩一区二区三区在线| 黄色视频不卡| 欧美色视频一区免费| 成人精品一区二区免费| 亚洲av成人不卡在线观看播放网| 国产男女超爽视频在线观看| 777米奇影视久久| 国产乱人伦免费视频| 99re6热这里在线精品视频| 精品久久久久久电影网| 亚洲一区中文字幕在线| 久久精品国产清高在天天线| bbb黄色大片| 色婷婷久久久亚洲欧美| 久久亚洲精品不卡| 在线观看午夜福利视频| av国产精品久久久久影院| av天堂在线播放| 精品免费久久久久久久清纯 | 国产男靠女视频免费网站| 亚洲一码二码三码区别大吗| 黄色丝袜av网址大全| 国产一区二区三区视频了| 国产精品一区二区在线观看99| 人妻久久中文字幕网| 黄色成人免费大全| 国产午夜精品久久久久久| 亚洲一区高清亚洲精品| 久久香蕉精品热| 亚洲人成77777在线视频| 一级片'在线观看视频| 999久久久精品免费观看国产| 中文字幕另类日韩欧美亚洲嫩草| 国产成人啪精品午夜网站| 老汉色∧v一级毛片| 久久中文字幕人妻熟女| 熟女少妇亚洲综合色aaa.| 国产精品亚洲一级av第二区| 一本一本久久a久久精品综合妖精| 99久久国产精品久久久| 777久久人妻少妇嫩草av网站| 黄片小视频在线播放| 极品人妻少妇av视频| 亚洲成人手机| 欧美av亚洲av综合av国产av| 一边摸一边抽搐一进一小说 | 国产欧美亚洲国产| 大码成人一级视频| 精品电影一区二区在线| 人成视频在线观看免费观看| 在线观看免费高清a一片| 久久天堂一区二区三区四区| 少妇的丰满在线观看| 成人av一区二区三区在线看| 电影成人av| 亚洲av成人一区二区三| 一a级毛片在线观看| 99热国产这里只有精品6| 亚洲专区国产一区二区| 精品国产一区二区三区久久久樱花| av一本久久久久| 18禁观看日本| 久久久久久久久久久久大奶| 欧美最黄视频在线播放免费 | 黑人巨大精品欧美一区二区mp4| 久久香蕉国产精品| 免费av中文字幕在线| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 淫妇啪啪啪对白视频| 亚洲欧美色中文字幕在线| 超色免费av| 91大片在线观看| 日韩欧美国产一区二区入口| 国产精品 国内视频| 91在线观看av| 久久久精品免费免费高清| 99香蕉大伊视频| 国产乱人伦免费视频| 日韩三级视频一区二区三区| 91麻豆av在线| 悠悠久久av| 天堂中文最新版在线下载| 免费av中文字幕在线| 国产三级黄色录像| 制服人妻中文乱码| 91成人精品电影| 男女免费视频国产| 怎么达到女性高潮| 国产一区在线观看成人免费| 男人操女人黄网站| www日本在线高清视频| 亚洲精品国产精品久久久不卡| 高清av免费在线| 亚洲全国av大片| 中国美女看黄片| 欧美乱码精品一区二区三区| 国产一卡二卡三卡精品| 精品久久久久久电影网| 亚洲三区欧美一区| 99国产综合亚洲精品| 我的亚洲天堂| 日本欧美视频一区| 99精品久久久久人妻精品| 亚洲精品一卡2卡三卡4卡5卡| videosex国产| 国内久久婷婷六月综合欲色啪| 久久婷婷成人综合色麻豆| 亚洲av熟女| 在线国产一区二区在线| 黄片小视频在线播放| 国产亚洲欧美98| 国产99白浆流出| 狠狠狠狠99中文字幕| 久久国产亚洲av麻豆专区| 国产av一区二区精品久久| 久久久久久亚洲精品国产蜜桃av| 国产亚洲一区二区精品| 亚洲国产精品一区二区三区在线| 少妇猛男粗大的猛烈进出视频| 亚洲精品一二三| 亚洲第一欧美日韩一区二区三区| 欧美日韩亚洲高清精品| 国产精品 国内视频| 在线观看舔阴道视频| 欧美日本中文国产一区发布| 国产精品 欧美亚洲| 日本五十路高清| 超色免费av| 亚洲一区中文字幕在线| 国产一区二区三区在线臀色熟女 | 丰满人妻熟妇乱又伦精品不卡| 日本五十路高清| 亚洲熟女精品中文字幕| 老司机在亚洲福利影院| 大陆偷拍与自拍| 91老司机精品| 午夜福利欧美成人| 亚洲色图 男人天堂 中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 大陆偷拍与自拍| 男男h啪啪无遮挡| 亚洲男人天堂网一区| 国产精品永久免费网站| 国产一区有黄有色的免费视频| 久久久精品国产亚洲av高清涩受| 国产精品久久久久久精品古装| 亚洲五月天丁香| 欧美精品啪啪一区二区三区| 精品久久久久久久毛片微露脸| 亚洲精品久久午夜乱码| 国产欧美日韩一区二区三| 大码成人一级视频| 黄色视频,在线免费观看| 色精品久久人妻99蜜桃| 亚洲精品av麻豆狂野| 午夜福利一区二区在线看| 18禁裸乳无遮挡动漫免费视频| 欧美+亚洲+日韩+国产| 亚洲熟妇熟女久久| 国产精品永久免费网站| 美女午夜性视频免费| 黄色片一级片一级黄色片| 国产99白浆流出| 亚洲国产欧美一区二区综合| 亚洲自偷自拍图片 自拍| 久久久水蜜桃国产精品网| av天堂久久9| 国产又爽黄色视频| 亚洲成a人片在线一区二区| 黄色 视频免费看| 久久九九热精品免费| 757午夜福利合集在线观看| 黄色视频,在线免费观看| 脱女人内裤的视频| 欧美老熟妇乱子伦牲交| 成人国语在线视频| 久久人人97超碰香蕉20202| 欧美色视频一区免费| 九色亚洲精品在线播放| 丝袜美腿诱惑在线| 成人精品一区二区免费| 天堂动漫精品| 在线播放国产精品三级| 精品国产一区二区三区四区第35| 老司机午夜十八禁免费视频| ponron亚洲| 18禁裸乳无遮挡免费网站照片 | 欧美大码av| 国产不卡一卡二| 欧美黄色淫秽网站| 亚洲 国产 在线| 欧美人与性动交α欧美精品济南到| 黄色成人免费大全| 久久精品国产亚洲av高清一级| 日韩免费高清中文字幕av| 亚洲人成77777在线视频| 中文字幕高清在线视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久人人做人人爽| 老司机福利观看| 亚洲国产精品sss在线观看 | 亚洲国产欧美一区二区综合| 免费观看精品视频网站| 亚洲男人天堂网一区| 欧美久久黑人一区二区| avwww免费| 狂野欧美激情性xxxx| 亚洲少妇的诱惑av| 亚洲人成伊人成综合网2020| 国产在线一区二区三区精| 美女 人体艺术 gogo| 女人被躁到高潮嗷嗷叫费观| 男人操女人黄网站| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩成人在线一区二区| 亚洲精品成人av观看孕妇| 一级片免费观看大全| 久久久久国产精品人妻aⅴ院 | av不卡在线播放| 在线看a的网站| 亚洲国产精品合色在线| 啪啪无遮挡十八禁网站| 水蜜桃什么品种好| 青草久久国产| 精品卡一卡二卡四卡免费| 国产男女超爽视频在线观看| 精品卡一卡二卡四卡免费| 午夜福利在线免费观看网站| 精品国产国语对白av| 欧美激情高清一区二区三区| 午夜福利影视在线免费观看| 欧美日韩亚洲综合一区二区三区_| 91九色精品人成在线观看| 精品久久久久久久毛片微露脸| 亚洲中文日韩欧美视频| 亚洲人成77777在线视频| 王馨瑶露胸无遮挡在线观看| 女人久久www免费人成看片| 美女福利国产在线| 老汉色av国产亚洲站长工具| 99精品在免费线老司机午夜| 国产亚洲欧美在线一区二区| 热re99久久国产66热| 99精国产麻豆久久婷婷| av电影中文网址| 黄频高清免费视频| 又黄又粗又硬又大视频| 亚洲av日韩在线播放| 精品福利观看| 大型av网站在线播放| a级毛片黄视频| 午夜福利在线观看吧| 欧美成人午夜精品| 午夜精品在线福利| 国产欧美日韩一区二区三区在线| 精品国产一区二区久久| 久久人人97超碰香蕉20202| 国产欧美日韩综合在线一区二区| 久久精品国产a三级三级三级| 香蕉国产在线看| 国产一区二区三区视频了| 两性夫妻黄色片| 欧美 日韩 精品 国产| 亚洲人成电影免费在线| 日韩欧美国产一区二区入口| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 熟女少妇亚洲综合色aaa.| 午夜福利乱码中文字幕| 国产精品永久免费网站| 18禁国产床啪视频网站| 欧美精品av麻豆av| 国产精品二区激情视频| 国产激情欧美一区二区| 露出奶头的视频| 欧美 亚洲 国产 日韩一| av网站免费在线观看视频| 国产视频一区二区在线看| 波多野结衣一区麻豆| 成人国产一区最新在线观看| 成人国语在线视频| avwww免费| 成人18禁高潮啪啪吃奶动态图| 欧美日本中文国产一区发布| 真人做人爱边吃奶动态| 一个人免费在线观看的高清视频| 亚洲自偷自拍图片 自拍| 最近最新中文字幕大全免费视频| 九色亚洲精品在线播放| 久久九九热精品免费| 嫩草影视91久久| 在线av久久热| 一区二区日韩欧美中文字幕| 欧美久久黑人一区二区| 成人亚洲精品一区在线观看| 国产免费男女视频| 老鸭窝网址在线观看| 久久久久久免费高清国产稀缺| 国产成人av教育| 国产又爽黄色视频| 伦理电影免费视频| 亚洲精品粉嫩美女一区| 欧美人与性动交α欧美软件| 一级,二级,三级黄色视频| 妹子高潮喷水视频| 精品无人区乱码1区二区| 精品一区二区三卡| 搡老乐熟女国产| 欧美激情极品国产一区二区三区| 欧美乱码精品一区二区三区| 久久久久久久久久久久大奶| 妹子高潮喷水视频| 免费不卡黄色视频| 午夜福利乱码中文字幕| 精品少妇久久久久久888优播|