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

    半滑舌鰨甘露糖結(jié)合凝集素相關(guān)絲氨酸蛋白酶1基因的克隆及表達(dá)分析*

    2020-06-09 11:32:40陳亞?wèn)|沙珍霞
    漁業(yè)科學(xué)進(jìn)展 2020年3期

    王 倩 張 雪 陳亞?wèn)| 沙珍霞,2

    半滑舌鰨甘露糖結(jié)合凝集素相關(guān)絲氨酸蛋白酶1基因的克隆及表達(dá)分析*

    王 倩1張 雪1陳亞?wèn)|2,3沙珍霞1,2①

    (1. 青島大學(xué)生命科學(xué)學(xué)院 青島 266071; 2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室海洋漁業(yè)科學(xué)與食物產(chǎn)出過(guò)程功能實(shí)驗(yàn)室 青島 266071;3. 中國(guó)水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所 農(nóng)業(yè)農(nóng)村部海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室 青島 266071)

    甘露糖結(jié)合凝集素相關(guān)絲氨酸蛋白酶1 (Mannan-binding lectin associated serine protease 1, MASP1)是補(bǔ)體凝集素途徑中起重要作用的激活蛋白。本研究以半滑舌鰨()為研究對(duì)象,應(yīng)用RACE技術(shù)和實(shí)時(shí)熒光定量qRT-PCR技術(shù)對(duì)半滑舌鰨1基因(11)進(jìn)行了克隆和表達(dá)模式分析。結(jié)果顯示,1基因cDNA序列全長(zhǎng)為2507 bp,5′非編碼區(qū)長(zhǎng)82 bp,3′非編碼區(qū)長(zhǎng)142 bp,開(kāi)放閱讀框長(zhǎng)2283 bp,共編碼760個(gè)氨基酸;1基因包含13個(gè)外顯子和12個(gè)內(nèi)含子,與多數(shù)已知魚(yú)類的1基因結(jié)構(gòu)一致;SMART分析顯示,CsMASP1包含6個(gè)結(jié)構(gòu)域,與哺乳動(dòng)物、鳥(niǎo)類、其他魚(yú)類的結(jié)構(gòu)域一致;同源比對(duì)發(fā)現(xiàn),CsMASP1和魚(yú)類的相似度較高,與金目鱸()的相似性最高,為76%。1基因在11種健康組織(血液、腦、鰓、性腺、心臟、頭腎、腸、肝、皮膚、脾和后腎)中均有表達(dá),其中,在肝、脾和頭腎中表達(dá)量較高;鰻弧菌()感染后,1基因在6種免疫組織(血液、鰓、頭腎、腸、肝和脾)中呈現(xiàn)不同的表達(dá)模式,在6種免疫組織中呈現(xiàn)明顯的上調(diào)表達(dá),肝的表達(dá)峰值出現(xiàn)在感染后24 h;脾和鰓的表達(dá)峰值出現(xiàn)在感染后6 h;腸、頭腎和血液的表達(dá)峰值均出現(xiàn)在感染后12 h;隨著病原菌感染時(shí)間增加,基因表達(dá)量逐漸降低并恢復(fù)至正常水平。研究表明,1基因參與了半滑舌鰨的免疫應(yīng)答過(guò)程,本結(jié)果為半滑舌鰨的免疫機(jī)理研究奠定了基礎(chǔ)。

    半滑舌鰨;1基因;基因克??;基因表達(dá);免疫應(yīng)答

    甘露糖結(jié)合凝集素相關(guān)絲氨酸蛋白酶1(MASP1)屬于絲氨酸蛋白酶(Serine proteases, SP)超家族成員,包括MASP1(Takahashi, 2013)、MASP2(Swierzko, 2014)和MASP3(Dahl, 2001)。MASP1是補(bǔ)體凝集素途徑中起重要作用的激活蛋白,MASP1可被甘露糖結(jié)合凝集素(Mannan-binding lectin, MBL)激活以響應(yīng)病原體感染,MBL-MASP復(fù)合物共同激活凝集素途徑(Sunyer, 1998; Schwaebl, 2002)。

    目前,關(guān)于MASP1的研究相對(duì)較少,2000年首次在人中克隆鑒定出1(Chen, 2000),隨后構(gòu)建了人的MASP1 N端片段原核表達(dá)載體并在大腸桿菌中進(jìn)行了表達(dá)(Cai, 2008)。迄今僅在文昌魚(yú)()(Endo, 2003)、七鰓鰻()(Endo, 2003)、草魚(yú)()(Dang, 2017)和菊海鞘()(Nicola, 2017)等水生生物中開(kāi)展了1的克隆表達(dá)及免疫應(yīng)答相關(guān)研究。硬骨魚(yú)MASP1的研究還極為有限。

    半滑舌鰨()是我國(guó)重要的海水養(yǎng)殖魚(yú)類。近年來(lái),隨著水產(chǎn)養(yǎng)殖行業(yè)的快速發(fā)展,病害問(wèn)題特別是弧菌病對(duì)半滑舌鰨養(yǎng)殖業(yè)造成嚴(yán)重的經(jīng)濟(jì)損失。有關(guān)半滑舌鰨免疫機(jī)制的研究和有效的病害防控手段還有待深入研究(Chen, 2015)。本研究通過(guò)PCR技術(shù)克隆了1基因的cDNA,并分析了其基因結(jié)構(gòu)及其進(jìn)化關(guān)系,還研究了1基因在半滑舌鰨11種健康組織中的表達(dá)特征,及鰻弧菌感染后該基因在6種免疫組織中的時(shí)空表達(dá)規(guī)律,為進(jìn)一步深入了解1基因在半滑舌鰨免疫防御中的作用奠定基礎(chǔ)。

    1 材料與方法

    1.1 實(shí)驗(yàn)材料

    半滑舌鰨購(gòu)自山東海陽(yáng)黃海水產(chǎn)有限公司,魚(yú)齡為2齡左右,平均體重為(400±5) g,平均體長(zhǎng)為(40.3±0.1) cm。健康實(shí)驗(yàn)魚(yú)在24℃水族箱中暫養(yǎng)7 d,每天換新鮮滅菌海水,盡量消除環(huán)境脅迫。鰻弧菌為本實(shí)驗(yàn)室保存菌種,經(jīng)16S rRNA測(cè)序進(jìn)行菌種鑒定。

    1.2 樣品處理和采集

    隨機(jī)選取3條健康的半滑舌鰨,分別采集每條魚(yú)的血液、腦、鰓、性腺、心臟、頭腎、腸、肝、皮膚、脾和后腎共11種組織,并立即投入液氮中冷凍,然后轉(zhuǎn)移至–80℃冰箱保存,用于總RNA提取。鰻弧菌感染參照(Sha, 2017)實(shí)驗(yàn)方法略作修改,采用腹腔注射方法對(duì)實(shí)驗(yàn)組進(jìn)行鰻弧菌感染,感染劑量(半致死劑量)為3.18×105CFU/g魚(yú)體,對(duì)照組注射相同體積的無(wú)菌PBS溶液。在注射后0、6、12、24、48和72 h共6個(gè)時(shí)間點(diǎn),每個(gè)時(shí)間點(diǎn)取3條魚(yú),分別收集實(shí)驗(yàn)組和對(duì)照組的血液、鰓、頭腎、腸、肝臟和脾臟6種組織,保存方式同上所述。

    1.3 RNA提取和cDNA制備

    采用液氮研磨收集的組織,按照RNA提取試劑盒(天根,北京)的說(shuō)明書(shū)進(jìn)行組織總RNA提取。提取的RNA用瓊脂糖凝膠電泳(1.00%)檢測(cè)其完整性。用Agilent 2100 biological analyzer (Applied Biosystems, 美國(guó))檢測(cè)其純度和濃度。按照PrimeScript Ⅱ 1st Strand cDNA Synthesis Kit (TaKaRa, 大連)試劑盒反轉(zhuǎn)錄成cDNA,–20℃冰箱保存?zhèn)溆谩?/p>

    1.4 CsMASP1基因的cDNA序列的擴(kuò)增

    通過(guò)對(duì)半滑舌鰨轉(zhuǎn)錄組和基因組數(shù)據(jù)進(jìn)行生物信息學(xué)分析,篩選了1基因的部分cDNA序列,并對(duì)該序列設(shè)計(jì)特異性引物P1和P2(表1),以半滑舌鰨11個(gè)組織混合cDNA為模板,進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系(20 μl):Trans?HiFi Polymerase (5 U/μl) 0.2 μl, 10×TransHiFi BufferⅡ2 μl, 2.5 mmol/L dNTPs 1.6 μl, P1 1 μl, P2 1 μl, cDNA模板2 μl, ddH2O 12.2 μl。反應(yīng)程序:95℃ 5 min;95℃ 30 s,58℃ 30 s,72℃ 30 s,38個(gè)循環(huán);72℃ 10 min;4℃保存。1.00%瓊脂糖凝膠電泳檢測(cè)后,用DNA膠回收試劑盒(諾唯贊, 南京)純化PCR擴(kuò)增特異性產(chǎn)物,連接到pEASY-T1(全式金, 北京)載體上,重組載體轉(zhuǎn)化到1感受態(tài)細(xì)胞(全式金, 北京)中,挑取陽(yáng)性克隆進(jìn)行菌液PCR驗(yàn)證后測(cè)序(上海派森諾生物科技有限公司)。在克隆獲得的序列5′端和3′端設(shè)計(jì)RACE (Rapid-Amplification of cDNA Ends)引物P3和P5(表1)和巢式特異性引物P4和P6(表1),使用SMARTerTMRACE cDNA Amplification Kit(TaKaRa, 大連)獲得該基因的5′和3′端序列,將克隆的序列進(jìn)行拼接獲得1基因的cDNA全長(zhǎng)序列。

    表1 半滑舌鰨1基因克隆和表達(dá)分析所用引物

    Tab.1 PCR primers used for cloning and expression analysis of MASP1

    1.5 CsMASP1基因序列分析

    將1的全長(zhǎng)cDNA序列提交到NCBI網(wǎng)站(http://blast.ncbi.nlm.nih.gov/blast.cgi),利用BLAST程序進(jìn)行核苷酸同源性比對(duì),利用ORF Finder (http://blast.ncbi.nlm.nih.gov/gorf/gorf.html)查找開(kāi)放閱讀框,并推導(dǎo)其相應(yīng)的氨基酸序列;采用SMART(Simple Modular Architecture Research Tool)在線軟件分析CsMASP1蛋白的結(jié)構(gòu)域,利用ExPASy軟件(http://www.expasy.org)預(yù)測(cè)其相對(duì)分子質(zhì)量和等電點(diǎn);用SignalP 4.0 Server (http://www.cbs.dtu. dk/services/SignalP/)軟件分析CsMASP1蛋白的信號(hào)肽。應(yīng)用Splign(http://www.ncbi.nlm.nih.gov/IEB/Research/Ostell/Spidey/)在線軟件確定1基因的內(nèi)含子與外顯子的數(shù)目,分析不同物種1的基因結(jié)構(gòu)以及其所在染色體的位置。

    1.6 MASP1同源比對(duì)系統(tǒng)進(jìn)化樹(shù)構(gòu)建

    從GenBank中查詢不同物種1基因的氨基酸序列,采用DNAMAN軟件進(jìn)行不同物種的MASP1蛋白的氨基酸序列多重比對(duì),用MEGA 6.0軟件中的Neighbor-joining法構(gòu)建系統(tǒng)進(jìn)化樹(shù)。

    1.7 CsMASP1基因表達(dá)的實(shí)時(shí)定量分析

    根據(jù)1基因設(shè)計(jì)擴(kuò)增引物P7和P8(表1),以半滑舌鰨18S rRNA基因?yàn)閮?nèi)參設(shè)計(jì)特異引物Cs18s-F和Cs18s-R (表1)。基因?qū)崟r(shí)定量表達(dá)參照已報(bào)道的方法進(jìn)行(Lu, 2013; Yu, 2017)。使用ChamQTMSYBR?Color qPCR Master Mix試劑(諾唯贊, 南京),按照說(shuō)明書(shū)于Applied Biosystems 7500 Real Time PCR儀上進(jìn)行該基因的實(shí)時(shí)定量分析(qRT-PCR)。反應(yīng)體系為20 μl:10 μl 2×ChamQ SYBR Color qPCR Master Mix,1 μl模板cDNA,0.4 μl引物P7,0.4 μl引物P8,0.4 μl 50×ROX Reference Dye 2,7.8 μl ddH2O。反應(yīng)程序:預(yù)變性95℃ 3 min;95℃ 10 s,60℃ 30 s,40個(gè)循環(huán);熔解曲線95℃ 15 s;60℃ 1 min;95℃ 15 s。每個(gè)樣品設(shè)置3個(gè)技術(shù)重復(fù),采用2–ΔΔCt方法計(jì)算目的基因的相對(duì)表達(dá)量,值取3個(gè)平行樣品的平均值。qRT-PCR檢測(cè)結(jié)果使用SPSS19.0統(tǒng)計(jì)軟件進(jìn)行單因素方差分析(One-way ANOVA)的Duncan法對(duì)多組樣本平均值進(jìn)行兩兩比較分析,<0.05時(shí),認(rèn)為存在顯著性差異。

    2 結(jié)果

    2.1 CsMASP1基因序列分析和結(jié)構(gòu)特征

    克隆獲得的1基因cDNA序列全長(zhǎng)為2507 bp (登錄號(hào):XM_008309208.2),5′非編碼區(qū)(Untranslated Region, UTR, 5′ UTR)長(zhǎng)82 bp,3′ UTR長(zhǎng)142 bp,開(kāi)放閱讀框(Open reading frame, ORF)長(zhǎng)2283 bp (圖1)。推測(cè)CsMASP1編碼氨基酸為760個(gè),預(yù)測(cè)其分子量為84.95 kDa,理論等電點(diǎn)(pI)為5.33。SignalP 4.0 Server分析顯示,此氨基酸序列無(wú)信號(hào)肽。SMART軟件分析顯示,CsMASP1蛋白含6個(gè)結(jié)構(gòu)域:分別是2個(gè)類Clr/Cls蛋白(like Clr/Cls Protein domain, CUB)結(jié)構(gòu)域(47~168 aa, 215~327 aa),1個(gè)表皮細(xì)胞生長(zhǎng)因子樣(Calcium-binding epidermal growth factor-like domain, EGF-CA)結(jié)構(gòu)域(169~212 aa),2個(gè)補(bǔ)體控制蛋白區(qū)(Complement control protein domain, CCP)結(jié)構(gòu)域(331~392 aa, 397~461 aa)和1個(gè)絲氨酸蛋白激酶區(qū)(Serine protease domain, SP)Tryp-SPc結(jié)構(gòu)域(478~752 aa) (圖2)。

    利用Splign軟件分析了7個(gè)物種的1基因結(jié)構(gòu)并進(jìn)行了比較(圖3),它們是半滑舌鰨(XM_ 008309208.2)、劍尾魚(yú)()(XM_ 023351221.1)、大黃魚(yú)()(XM_ 027280128.1)、斑馬魚(yú)()(XM_001341900. 6)、馬() (XM_001499629.5)、人()(NM_001879.5)和中國(guó)倉(cāng)鼠() (XM_003495604.2)。結(jié)果顯示,1含13個(gè)外顯子,位于4號(hào)染色體上;劍尾魚(yú)1含13個(gè)外顯子,位于18號(hào)染色體上;大黃魚(yú)1含13個(gè)外顯子,位于7號(hào)染色體上;斑馬魚(yú)1含12個(gè)外顯子,位于15號(hào)染色體上;馬1含16個(gè)外顯子,位于19號(hào)染色體上;人1含16個(gè)外顯子,位于3號(hào)染色體上;中國(guó)倉(cāng)鼠1含15個(gè)外顯子,位于4號(hào)染色體上。半滑舌鰨1基因結(jié)構(gòu)與多數(shù)已知魚(yú)類(劍尾魚(yú)和大黃魚(yú))一致,與哺乳類(人、馬和中國(guó)倉(cāng)鼠)在外顯子數(shù)目上有差異,說(shuō)明該基因在魚(yú)類保守。

    2.2 半滑舌鰨與其他物種MASP1同源比對(duì)和系統(tǒng)進(jìn)化分析

    通過(guò)BLAST比對(duì),將CsMASP1氨基酸序列(XP_008307430.1)與金目鱸()(XP_ 018560455.1)、牙鲆()(XP_ 019966501.1)、黃鱔()(XP_020477488.1)、伯氏樸麗魚(yú)()(XP_014184888.1)、非洲爪蟾()(AAI70221.1)、海龜()(EMP37646.1)、紅喉潛鳥(niǎo)() (XP_009815055.1)、皇帝企鵝() (KFM09600.1)、白鷺()(XP_009633907.1)、西伯利亞虎() (XP_015393589.1)、羊駝()(XP_ 006201039.1)、馬(XP_ 005601936.1)、中國(guó)倉(cāng)鼠(XP_007606939.1)、人(XP_ 016862358.1)的MASP1氨基酸序列進(jìn)行同源比對(duì),相似性范圍為76%~42%。其中,CsMASP1與金目鱸的相似性最高,為76%;而與中國(guó)倉(cāng)鼠的相似性最低,為42%。將半滑舌鰨和其他13個(gè)物種MASP1的氨基酸序列進(jìn)行多重比對(duì),可以看出不同物種MASP1存在高度保守的氨基酸殘基,如圖4所示,保守氨基酸出現(xiàn)在絲氨酸蛋白酶活性中心處:組氨酸His526、天冬氨酸Asp582和絲氨酸Ser705,這3個(gè)氨基酸的保守定位對(duì)于維持絲氨酸蛋白酶的活性是十分重要的。

    推導(dǎo)的氨基酸序列顯示在核苷酸序列下方,用大寫(xiě)字母表示。方框所示為起始密碼子(ATG)、星號(hào)(*)所示為終止密碼子(TGA);實(shí)線()所示為CUB結(jié)構(gòu)域;圓點(diǎn)線()所示為EGF-CA結(jié)構(gòu)域;短劃線()所示為CCP結(jié)構(gòu)域;雙劃線()所示為T(mén)ryp-SPc結(jié)構(gòu)域;波浪線()所示為poly A結(jié)構(gòu)

    Translated amino acid sequence was shown under the nucleotide sequence as uppercase. The iniation codon and termination codon were marked by a box and an asterisk, respectively; CUBdomain was shown in solid line(); EGF-CA domain was shown in dot line();CCPdomain was shown in dashed line();Tryp-SPcdomain was shown in double solid line(); poly A was shown in wavy line()

    利用MEGA 6.0軟件對(duì)上述各個(gè)物種的MASP1氨基酸序列進(jìn)行比對(duì)構(gòu)建了系統(tǒng)進(jìn)化樹(shù)(圖5)。結(jié)果顯示,半滑舌鰨與伯氏樸麗魚(yú)、金目鱸和黃鱔等魚(yú)類的MASP1聚為一支。兩棲類海龜和非洲爪蟾、鳥(niǎo)類白鷺、皇帝企鵝和紅喉潛鳥(niǎo),哺乳類馬、西伯利亞虎、羊駝、人和中國(guó)倉(cāng)鼠等動(dòng)物的MASP1聚為一支。但兩棲類、鳥(niǎo)類和哺乳類又按親緣關(guān)系遠(yuǎn)近聚成相應(yīng)的亞類。MASP1在進(jìn)化上與伯氏樸麗魚(yú)的親緣關(guān)系最近,與哺乳類、爬行類和鳥(niǎo)類等動(dòng)物的MASP1的親緣關(guān)系較遠(yuǎn)。從進(jìn)化的結(jié)果看,幾個(gè)關(guān)鍵節(jié)點(diǎn)的支持率較高,均在90%以上。例如:金目鱸和黃鱔,半滑舌鰨和伯氏樸麗魚(yú)這幾個(gè)物種聚類為一個(gè)分支,支持率分別為98%和100%。

    圖2 半滑舌鰨 MASP1 蛋白結(jié)構(gòu)域預(yù)測(cè)

    SMART軟件在線預(yù)測(cè)CsMASP1蛋白的結(jié)構(gòu)域,CsMASP1包含6個(gè)結(jié)構(gòu)域,從C端到N端分別是CUB結(jié)構(gòu)域(47~168 aa), EGF-CA結(jié)構(gòu)域(169~212aa), CUB結(jié)構(gòu)域(215~327 aa),CCP結(jié)構(gòu)域(331~392 aa),CCP結(jié)構(gòu)域(397~461 aa), Tryp-SPc結(jié)構(gòu)域(478~752 aa)

    SMART software predicts the domain of CsMASP1 protein online. CsMASP1 contains six domains:CUB domain (47~168 aa), EGF-CA domain (169~212 aa), CUB domain (215~327 aa), CCP domain (331~392 aa),CCP domain (397~461 aa), Tryp-SPc domain (478~752 aa)

    圖3 不同物種的MASP1基因結(jié)構(gòu)比較

    7個(gè)物種的1基因結(jié)構(gòu)圖,各個(gè)物種的1基因序列號(hào)與所在染色體基因組序列號(hào)分別為:人NM_001879.5, NG_029440.1;馬XM_001499629.5, NC_009162.3; 中國(guó)倉(cāng)鼠M_003495604.2, NW_003613596.1; 劍尾魚(yú)XM_023351221.1, NC_036460.1; 大黃魚(yú)XM_027280128.1, NC_040017.1; 半滑舌鰨XM_008309208.2, NC_024310.1; 斑馬魚(yú)XM_001341900.6, NC_007126.7

    The sequence numbers of the1 gene and the genome numbers of the chromosomes in seven species were as follows:NM_001879.5, NG_029440.1;XM_001499629.5, NC_009162.3;M_003495604.2, NW_003613596.1;XM_023351221.1, NC_036460.1;XM_027280128.1, NC_040017.1;XM_008309208.2, NC_024310.1;XM_001341900.6, NC_007126.7

    圖4 半滑舌鰨和其他物種MASP1氨基酸序列的多重比對(duì)

    MASP1氨基酸序列比對(duì),絲氨酸蛋白酶結(jié)構(gòu)域中保守氨基酸用黑框標(biāo)出:組氨酸His526、天冬氨酸Asp582和絲氨酸Ser705

    Amino acid sequence alignment of MASP1, conservative amino acids in serine protease domain were marked by the black box: His526, Asp582and Ser705

    圖5 MEGA 6.0軟件構(gòu)建的基于MASP1氨基酸序列的系統(tǒng)進(jìn)化樹(shù)

    2.3 CsMASP1基因在健康組織中的表達(dá)分析

    利用qRT-PCR分析1基因在11種健康組織中的表達(dá)特征。結(jié)果顯示,1基因在血液、腦、鰓、性腺、心臟、頭腎、腸、肝、皮膚、脾和后腎中均有表達(dá);其中,肝中相對(duì)表達(dá)量最高(185.4),其次為脾(138.1)和頭腎(76.2),在腦、心臟、血液、腸、后腎、性腺、皮膚和鰓中表達(dá)量較低且表達(dá)量依次降低(圖6)。

    2.4 鰻弧菌感染后CsMASP1基因在免疫組織中的表達(dá)分析

    鰻弧菌感染半滑舌鰨后可誘導(dǎo)1在免疫組織中的表達(dá)變化,在感染0、6、12、24、48、72 h后,1在PBS對(duì)照組和感染組半滑舌鰨的血液、鰓、頭腎、腸、肝和脾6種免疫組織中的表達(dá)水平如圖7所示。1基因表達(dá)總體趨勢(shì)表現(xiàn)為先升高后降低。在脾和鰓中,1基因的表達(dá)呈現(xiàn)相同的上調(diào)趨勢(shì),其表達(dá)峰值出現(xiàn)在鰻弧菌感染后6 h,在脾中為對(duì)照組表達(dá)量的15.4倍;在鰓中則是對(duì)照組的4.6倍;在肝中1基因表達(dá)上調(diào),最大表達(dá)量出現(xiàn)在感染后24 h,是對(duì)照組的12.4倍;在腸中1基因表達(dá)上調(diào),最大表達(dá)量出現(xiàn)在感染后12 h,是對(duì)照組的7.7倍,72 h時(shí)依然保持較高的表達(dá)水平;在頭腎和血液中,1基因的表達(dá)峰值均出現(xiàn)在鰻弧菌感染后12 h,在頭腎中為對(duì)照組表達(dá)量的9.7倍;在血液中則是對(duì)照組的6.8倍;隨后基因表達(dá)量逐漸下降至略低于正常水平值。

    3 討論

    3.1 CsMASP1基因cDNA的克隆及序列分析

    甘露聚糖結(jié)合凝集素相關(guān)絲氨酸蛋白酶是機(jī)體重要的天然免疫防御分子,通過(guò)激活補(bǔ)體和調(diào)理吞噬2種方式來(lái)清除病原體(Chen, 2000)。MASP1是形成C4bC2a轉(zhuǎn)化酶的必要因子,是凝集素激活途徑中的特異性啟動(dòng)酶(Oroszlan, 2016)。MASP最初是作為與Complement component 1s(Cls)功能相似可以裂解C4和C2的一種新型絲氨酸蛋白酶被發(fā)現(xiàn)的(Matsushita, 1992)。鑒于MASP與Cls相似,推測(cè)MBL-MASP1-MASP2與C1q-C1r2-C1s2相似,因此,MASP1和MASP2可能分別與類Complement component 1r (Clr)和類C1s絲氨酸蛋白酶有相似的功能(Matsushita, 1998、2000)。雖然在硬骨魚(yú)的免疫系統(tǒng)中已經(jīng)鑒定出了1,但其對(duì)細(xì)菌感染的免疫功能尚不清楚(Endo, 2003; Dang, 2017)。本研究成功克隆了1基因,預(yù)測(cè)其編碼蛋白的分子量為84.95 kDa,與已報(bào)道的草魚(yú)MASP1 (Dang, 2017)的分子量和人MASP1 (Matsushita, 1998)的分子量類似。MASP1編碼的氨基酸序列經(jīng)SignalP軟件預(yù)測(cè)無(wú)信號(hào)肽,而草魚(yú)MASP1蛋白序列中包含一段19 aa的信號(hào)肽(Dang, 2017)。研究表明,MASP的SP區(qū)(Serine protease domain)是催化功能區(qū),含有所有絲氨酸蛋白酶活性中心高度保守的組氨酸(His)、丙氨酸(Asp)和絲氨酸(Ser)等3個(gè)氨基酸殘基,與MASP的絲氨酸蛋白酶活性有關(guān)(Wallis, 2000、2007)。在本研究中,SMART結(jié)果顯示,MASP1含有6個(gè)結(jié)構(gòu)域,分別是2個(gè)CUB結(jié)構(gòu)域(47~168 aa, 215~327 aa)、1個(gè)EGF結(jié)構(gòu)域(169~212 aa)、2個(gè)CCP結(jié)構(gòu)域(331~392 aa, 397~461 aa)和1個(gè)Tryp-SPc結(jié)構(gòu)域(478~752 aa)。與文獻(xiàn)報(bào)道的人MASP1(Matsushita, 1998; Stover, 2001; Dobo, 2016)結(jié)構(gòu)域相同,與草魚(yú)MASP1(Dang, 2017)結(jié)構(gòu)域也相同。絲氨酸蛋白酶結(jié)構(gòu)域(SP)活性中心有1個(gè)典型的組氨酸環(huán)。氨基酸序列比對(duì)結(jié)果顯示,絲氨酸蛋白酶結(jié)構(gòu)域完全相似,組氨酸His526、天冬氨酸Asp582和絲氨酸Ser705定位相當(dāng)保守,與文獻(xiàn)報(bào)道的文昌魚(yú)和七鰓鰻(Endo, 2003)的比對(duì)結(jié)果一致,說(shuō)明該活性中心在進(jìn)化上相對(duì)保守。這些特性與MASP/C1r/C1s家族的進(jìn)化相關(guān)(Jia, 2003)。系統(tǒng)進(jìn)化分析結(jié)果顯示,MASP1與伯氏樸麗魚(yú)MASP1親緣關(guān)系最近,與金目鱸和黃鱔等魚(yú)類的MASP1親緣關(guān)系次之。

    圖6 CsMASP1基因在半滑舌鰨健康組織中的相對(duì)表達(dá)量

    Bl: 血液;Br:腦;Gi:鰓;Go:性腺;He:心臟;HK:頭腎;In:腸;Li:肝;Sk:皮膚;Sp:脾;TK:后腎。字母“a, b, c”代表SPSS多重分析的不同分組,有相同字母表示差異不顯著(>0.05),無(wú)相同字母表示差異顯著(<0.05);下同

    Bl: Blood; Br: Brain; Gi: Gills; Go: Gonad; He: Heart; Hk: Head kidney; In: Intestine; Li: Liver; Sk: Skin; Sp: Spleen; Tk: Trunk kidney. The letters of ‘a(chǎn), b, and c’ indicated the Duncan grouping in SPSS 19.0 software. The same letters indicate no significant difference (>0.05), the different letters indicate significant difference (<0.05); the same applies below

    圖7 鰻弧菌感染半滑舌鰨后不同時(shí)間點(diǎn)6種組織(血液、鰓、頭腎、腸、肝和脾)中CsMASP1基因的相對(duì)表達(dá)水平

    3.2 CsMASP1在健康組織和鰻弧菌感染后的免疫組織中的表達(dá)分析

    本研究利用qRT-PCR檢測(cè)1在半滑舌鰨11種健康組織和感染后6種免疫組織中的表達(dá)譜,結(jié)果顯示,在所檢測(cè)的11種組織中分布廣泛且具有組織特異性,在肝中表達(dá)水平最高,脾中表達(dá)水平次之,在頭腎中的表達(dá)相對(duì)其他組織也較高,推測(cè)其可能在免疫調(diào)控中發(fā)揮作用。與草魚(yú)的分布結(jié)果不相同,草魚(yú)1在心臟中表達(dá)水平最高,在肝和腦中表達(dá)水平次之(Dang, 2017)。

    MASP1是構(gòu)成先天性免疫防御體系的重要成員之一(Ammitzboll, 2013)。本研究報(bào)道了1參與鰻弧菌感染后機(jī)體反應(yīng)的免疫應(yīng)答。鰻弧菌感染半滑舌鰨72 h內(nèi),1在血液、鰓、頭腎、腸、肝和脾組織中均出現(xiàn)上調(diào)表達(dá)的情況,總體趨勢(shì)表現(xiàn)為表達(dá)量先升高,達(dá)到一個(gè)峰值后,隨感染時(shí)間的延長(zhǎng)表達(dá)量逐漸降低恢復(fù)至正常水平,說(shuō)明1在病原入侵后參與了機(jī)體的免疫應(yīng)答。鰓和腸中的1表達(dá)量在鰻弧菌感染后6 h大幅上調(diào),鰓在6 h達(dá)到最大值,腸比鰓稍微滯后,腸在12 h達(dá)到最高峰,分別是對(duì)照組的7.7倍和4.6倍,隨后下降到與對(duì)照組無(wú)顯著差異。腸是最大的免疫器官,腸上分布著大量的腸相關(guān)淋巴組織,腸道組織和鰓是腹腔注射鰻弧菌后直接接觸的部位,所以腸和鰓1的表達(dá)量均先于其他免疫組織上調(diào)。隨后,病原菌進(jìn)入其他免疫組織,激活其他免疫組織1基因的表達(dá)。與Dang等(2017)研究發(fā)現(xiàn)草魚(yú)感染嗜水氣單胞菌后,草魚(yú)1在鰓和腸中表達(dá)量短時(shí)快速增加結(jié)果相一致;但草魚(yú)1在腸中感染后4、24、72 h顯示3個(gè)表達(dá)高峰與本研究的結(jié)果不太一致,推測(cè)草魚(yú)1在腸中感染后出現(xiàn)3個(gè)峰值可能是由于感染菌不同,也可能是由于魚(yú)體差異等。脾1的表達(dá)量隨著時(shí)間增加而增加,在感染6 h后達(dá)到最大值,約是對(duì)照組的15倍,是所有免疫組織中增量最高的組織;這與腎、脾是魚(yú)類免疫細(xì)胞生成的主要場(chǎng)所,也是產(chǎn)生免疫應(yīng)答的主要部位這一觀點(diǎn)是一致的(Wang, 2018)。6種免疫組織表現(xiàn)出不同的表達(dá)模式可能與不同時(shí)間點(diǎn)的信號(hào)轉(zhuǎn)導(dǎo)有關(guān),在我們之前的研究中,感染后,補(bǔ)體激活途徑、凝血級(jí)聯(lián)激活途徑和弧菌感染相關(guān)途徑均被激活來(lái)抵御病原菌的入侵(Zhang, 2015)。1在免疫組織中的高表達(dá)暗示了1在機(jī)體免疫反應(yīng)中可能發(fā)揮著重要作用。本結(jié)果為深入研究1基因在半滑舌鰨免疫防御系統(tǒng)中的作用奠定了基礎(chǔ)。

    Ammitzboll CG, Steffensen R, Nielsen HJ,Polymorphisms in the1 gene are associated with serum levels of MASP-1, MASP-3, and MAp44. PLoS One, 2013, 8(9): e73317

    Cai XM, Zhao N, Zuo DM,. Construction and expression of prokaryotic expression vector of human MASP1 N-terminal fragment. Cell and Molecular Immunology, 2008, 24(6): 546–549 [蔡學(xué)敏, 趙娜, 左大明, 等. 人MASPlN端片段原核表達(dá)載體的構(gòu)建及其表達(dá). 細(xì)胞與分子免疫學(xué)雜志, 2008, 24(6): 546–549]

    Chen ZL, Lu X, Zhang LY,. Cloning and identification of human MBL-related serine protease 1 gene. Cellular and Molecular Immunology, 2000, 16(5): 417–419 [陳政良, 盧曉, 張麗蕓, 等. 人MBL相關(guān)絲氨酸蛋白酶1 cDNA的克隆與鑒定. 細(xì)胞與分子免疫學(xué)雜, 2000, 16(5): 417–419]

    Chen Y, Fan L, Liu TT,. Cloning and expression analysis of6 gene and1 gene of half-smooth tongue sole. Chinese Aquatic Science, 2015, 22(5): 867–876 [陳燕, 樊琳, 劉田田, 等. 半滑舌鰨6基因和1基因的克隆及表達(dá)分析. 中國(guó)水產(chǎn)科學(xué), 2015, 22(5): 867–876]

    Dahl MR, Thiel S, Matsushita M,. MASP-3 and its association with distinct complexes of the mannan-binding lectin complement activation pathway. Immunity, 2001, 15(1): 127–135

    Dang YF, Shen YB, Xu XY,Mannan-binding lectin- associated serine protease-1 (MASP-1) mediates immune responses againstandin grass carp. Fish and Shellfish Immunology, 2017(66): 93–102

    Dobo J, Pal G, Cervenak L,. The emerging roles of mannose-binding lectin-associated serine proteases (MASPs) in the lectin pathway of complement and beyond. Immunological Reviews, 2016, 274(1): 98–111

    Endo Y, Nonaka M, Saiga H,Origin of mannose-binding lectin-associated serine protease (MASP)-1 and MASP-3 involved in the lectin complement pathway traced back to the invertebrate, amphioxus. Journal of Immunology, 2003, 170(9): 4701–4707

    Jia TJ, Li P, Zhang SM. Research status of structure and function of MASPs. Medical Review, 2003, 9(4): 193–195 [賈天軍, 李萍, 張庶民. MASPs的結(jié)構(gòu)與功能研究現(xiàn)狀. 醫(yī)學(xué)綜述, 2003, 9(4): 193–195]

    Matsushita M, Fujita T. Activation of the classical complement pathway by mannose-binding protein in association with a novel Cls-like serine protease. Journal of Experimental Medicine, 1992, 176(6): 1497–1502

    Matsushita M, Endo Y, Nonaka M,. Complement-related serine proteases in tunicates and vertebrates. Current Opinion in Immunology, 1998, 10(1): 29–35

    Matsushita M, Thiel S, Jensenius JC,Proteolytic activities of two types of mannose-binding lectin-associated serine protease. Journal of Immunology, 2000, 165(5): 2637–2642

    Matsushita M, Endo Y, Fujita T. MASPl (MBL-associated serine protease 1). Immunobiology, 1998, 199(2): 340–347

    Nicola F, Loriano B. Morula cells as key hemocytes of the lectin pathway of complement activation in the colonial tunicate. Fish and Shellfish Immunology, 2017, 63: 157–164

    Oroszlan G, Kortvely E, Szakacs D,MASP-1 and MASP-2 do not activate pro-factor D in resting human blood, whereas MASP-3 is a potential activator: kinetic analysis involving specific MASP-1 and MASP-2 inhibitors. Journal of Immunology, 2016, 196(2): 857–865

    Schwaebl W, Dalai MR, Thiel S,. The mannan-binding lectin-associated serine proteases (MASPs) and MAp19: Four components of the lectin pathway activation complex encoded by two genes. Immunobiology, 2002, 205(4–5): 455–466

    Sha ZX, Wang LQ, Sun LM,Isolation and characterization of monocyte/macrophage from peripheral blood of half- smooth tongue sole (). Fish and Shellfish Immunology, 2017, (65): 256–266

    Stover C, Endo Y, Takahashi M,. The human gene for mannan binding lectin associated serine protease-2, the effector component of the lectin route of complement activation, is part of a tightly linked gene cluster on chromosome 1p36.2-3. Genes and Immunity, 2001, 2(3): 119–127

    Sunyer JO, Lambris JD. Evolution and diversity of the complement system of poikilothermic vertebrates. Immunological Reviews, 1998, 166(1): 39–57

    Swierzko AS, Szala A, Sawicki S,. Mannose-binding lectin (MBL) and MBL-associated serineprotease-2 (MASP-2) in women with malignant and benign ovarian tumours. Cancer Immunology Immunotherapy, 2014, 63(11): 1129–1140

    Takahashi M, Sekine H, Endo Y,. Comment on “Mannan-binding lectin-associated serine protease (MASP)-1 is crucial for lectin pathway activation in human serum, whereas neither MASP-1 nor MASP-3 is required for alternative pathway function”. Journal of Immunology, 2013, 190(6): 2477

    Wang LQ, Wang H, Chen YD,Isolation and identification of T lymphocyte from peripheral blood of half-smooth tongue sole and analysis ofgene immune response. Acta Hydrobiologica Sinica, 2018, 42(3): 480–487 [汪林慶, 王航, 陳亞?wèn)|, 等. 半滑舌鰨外周血T淋巴細(xì)胞的分離、鑒定及基因免疫應(yīng)答分析. 水生生物學(xué)報(bào), 2018, 42(3): 480–487]

    Wallis R, Dodd RB. Interaction of mannose-binding protein with associated serine proteases: Effects of naturally occurring mutations. Journal of Biological Chemistry, 2000, 275(40): 30962–30969

    Wallis R, Dodds AW, Mitchell DA,. Molecular interactions between MASP-2, C4, and C2 and their activation fragments leading to complement activation via the lectin pathway. Journal of Biological Chemistry, 2007, 282(11): 7844–7851

    Yu MJ, Sun LM, Chen YD,. Cloning and expression of thegene in immune response to pathogen challenge in half-smooth tongue sole (). Progress in Fisheries Science, 2017, 38(2): 31–39 [于孟君, 孫璐明, 陳亞?wèn)|, 等. 半滑舌鰨()基因的克隆及免疫應(yīng)答分析. 漁業(yè)科學(xué)進(jìn)展, 2017, 38(2): 31– 39]

    Zhang X, Wang SL, Chen SL,. Transcriptome analysis revealed changes of multiple genes involved in immunity induringinfection. Fish and Shellfish Immunology, 2015, 43(1): 209–218

    Cloning and Expression Analysis of the1 Gene from the Half-Smooth Tongue Sole ()

    WANG Qian1, ZHANG Xue1, CHEN Yadong2,3, SHA Zhenxia1,2①

    (1. College of Life Science, Qingdao University, Qingdao 266071; 2. Laboratory for Marine Fisheries Science and Food Production Processes, PilotNational Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071; 3. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Key Laboratory for Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Qingdao 266071)

    Mannan-binding lectin-associated serine protease 1 (MBL associated serine protease 1, MASP1) is an important activator in the complement lectin pathway. In this study, the cDNA of the1 gene from(1) was cloned using a RACE method based on the partial sequence of1, and gene expression was performed by qRT-PCR. The results showed that the full length of cDNA was 2507 bp in size, including a 5′-untranslated region (UTR) of 82 bp, 3′-UTR of 142 bp, and a complete open reading frame (ORF) of 2283 bp, encoding 760 amino acids. The theoretical isoelectric point (pI) of the predicted protein was 5.33 and the molecular weight was 84.95 kDa. Homologous alignment showed that the amino acids sequence of1 had a high identity with those of other species, approximately 42%~76%. The1 gene was expressed in all tested tissues (liver, intestine, spleen, head-kidney, gill, blood, brain, skin, heart, trunk-kidney, and gonad) in the healthyand the highest expression was in the liver (185.4). To study the expression patterns of the1 gene in an immune response, the specific expression of1 was performed afterinfection. The results showed that the expression of the1 gene was up-regulated in the liver, gill, blood, intestine, head-kidney, and spleen afterinfection. The most significantly up-regulated expression and the peak level at 6 h reached 15.4 times baseline in the spleen. The results indicate that the1 genes are involved in the immune response.

    ;1; Gene cloning; Gene expression; Immune response

    S971.4

    A

    2095-9869(2020)03-0049-11

    10.19663/j.issn2095-9869.20190122002

    沙珍霞,教授,E-mail: shazhenxia@163.com

    2019-01-22,

    2019-03-09

    * 國(guó)家自然科學(xué)基金項(xiàng)目(31572644)資助 [This work was supported by National Natural Science Foundation of China (31572644)]. 王 倩,E-mail: qianwangqd@163.com

    http://www.yykxjz.cn/

    王倩, 張雪, 陳亞?wèn)|, 沙珍霞. 半滑舌鰨甘露糖結(jié)合凝集素相關(guān)絲氨酸蛋白酶1基因的克隆及表達(dá)分析. 漁業(yè)科學(xué)進(jìn)展, 2020, 41(3): 49–59

    Wang Q, Zhang X, Chen YD, Sha ZX. Cloning and expression analysis of the1 gene from the half-smooth tongue sole (). Progress in Fishery Sciences, 2020, 41(3): 49–59

    SHA Zhenxia, E-mail: shazhenxia@163.com

    (編輯 馮小花)

    日本黄色片子视频| 久久这里只有精品中国| 99久久九九国产精品国产免费| 九草在线视频观看| 日日摸夜夜添夜夜添av毛片| 又爽又黄a免费视频| 亚洲av在线观看美女高潮| 国产久久久一区二区三区| 亚洲精品乱码久久久v下载方式| 久久草成人影院| 特大巨黑吊av在线直播| 蜜桃亚洲精品一区二区三区| 国产乱来视频区| 欧美3d第一页| 极品教师在线视频| 国产视频内射| 国产精品一二三区在线看| 精品人妻熟女av久视频| 国产人妻一区二区三区在| 国产精品美女特级片免费视频播放器| 99久久精品国产国产毛片| 久久午夜福利片| 欧美成人一区二区免费高清观看| 一个人观看的视频www高清免费观看| av在线天堂中文字幕| 国产高清国产精品国产三级 | 午夜视频国产福利| 国产伦精品一区二区三区四那| 久久久久久久久中文| av专区在线播放| 国产 一区 欧美 日韩| 国产精品麻豆人妻色哟哟久久 | 日韩精品有码人妻一区| 精品人妻视频免费看| 黄片wwwwww| 少妇丰满av| 亚洲乱码一区二区免费版| 22中文网久久字幕| 免费看a级黄色片| 少妇被粗大猛烈的视频| 婷婷六月久久综合丁香| 精品国产三级普通话版| 亚洲国产av新网站| 亚洲av成人精品一二三区| 亚洲精品成人久久久久久| 精品国产三级普通话版| 亚洲av男天堂| 久久99热这里只有精品18| 国产一区亚洲一区在线观看| 日韩大片免费观看网站| 日本wwww免费看| 少妇丰满av| 久久久成人免费电影| 国产 一区精品| 亚洲精品视频女| 国产成人精品久久久久久| 寂寞人妻少妇视频99o| 一个人看的www免费观看视频| 久久久久国产网址| 国产伦在线观看视频一区| 18+在线观看网站| 中文字幕久久专区| 亚洲三级黄色毛片| 免费少妇av软件| 夜夜爽夜夜爽视频| 国产淫语在线视频| 一级毛片aaaaaa免费看小| 成人一区二区视频在线观看| 国产亚洲精品av在线| 国产精品99久久久久久久久| 熟妇人妻久久中文字幕3abv| 看免费成人av毛片| 黄色一级大片看看| 亚洲成色77777| 日本免费在线观看一区| 精品酒店卫生间| 欧美激情国产日韩精品一区| 亚洲无线观看免费| 91狼人影院| 九草在线视频观看| 你懂的网址亚洲精品在线观看| 亚洲国产高清在线一区二区三| 国产免费一级a男人的天堂| av网站免费在线观看视频 | 美女内射精品一级片tv| 青青草视频在线视频观看| or卡值多少钱| 免费不卡的大黄色大毛片视频在线观看 | 男人舔奶头视频| 男女边摸边吃奶| 日韩人妻高清精品专区| 日韩精品青青久久久久久| 婷婷色综合大香蕉| 五月玫瑰六月丁香| 丰满乱子伦码专区| 91久久精品国产一区二区成人| 18禁在线无遮挡免费观看视频| 美女cb高潮喷水在线观看| 免费黄色在线免费观看| 欧美极品一区二区三区四区| 一级爰片在线观看| 免费看不卡的av| 免费大片黄手机在线观看| 欧美xxⅹ黑人| 婷婷色麻豆天堂久久| 深夜a级毛片| 精品久久久久久电影网| 亚洲18禁久久av| 亚洲国产欧美人成| 在线观看免费高清a一片| 男插女下体视频免费在线播放| 久99久视频精品免费| 麻豆av噜噜一区二区三区| av黄色大香蕉| 免费人成在线观看视频色| 久久鲁丝午夜福利片| av福利片在线观看| 三级男女做爰猛烈吃奶摸视频| 久久精品国产亚洲av涩爱| 亚洲三级黄色毛片| 精品国内亚洲2022精品成人| 国产精品99久久久久久久久| av天堂中文字幕网| 边亲边吃奶的免费视频| eeuss影院久久| 亚洲熟妇中文字幕五十中出| 一级毛片久久久久久久久女| 欧美成人一区二区免费高清观看| 欧美3d第一页| 18禁动态无遮挡网站| 大又大粗又爽又黄少妇毛片口| 日日摸夜夜添夜夜添av毛片| 国产精品1区2区在线观看.| 亚洲欧美精品自产自拍| 亚洲真实伦在线观看| 久久这里只有精品中国| 成人特级av手机在线观看| 欧美zozozo另类| 成人性生交大片免费视频hd| 成人鲁丝片一二三区免费| 老女人水多毛片| 亚洲最大成人中文| 日本-黄色视频高清免费观看| 亚洲乱码一区二区免费版| www.色视频.com| 搡老乐熟女国产| 好男人视频免费观看在线| 99热6这里只有精品| 久久草成人影院| 亚洲精品影视一区二区三区av| 老司机影院毛片| 黑人高潮一二区| 亚洲精品一区蜜桃| 国产精品一区二区在线观看99 | 国产成人freesex在线| 特级一级黄色大片| 成人亚洲欧美一区二区av| 免费高清在线观看视频在线观看| av在线天堂中文字幕| 久久精品夜色国产| 国产国拍精品亚洲av在线观看| 色尼玛亚洲综合影院| 大片免费播放器 马上看| 最近最新中文字幕免费大全7| 亚洲精品自拍成人| 国产成人精品福利久久| 一级爰片在线观看| 欧美 日韩 精品 国产| 毛片女人毛片| 欧美一区二区亚洲| 亚洲精品乱久久久久久| 少妇人妻一区二区三区视频| 国产精品美女特级片免费视频播放器| 又粗又硬又长又爽又黄的视频| 亚洲精品aⅴ在线观看| 岛国毛片在线播放| 国产精品人妻久久久久久| 免费黄网站久久成人精品| 日韩精品有码人妻一区| 大陆偷拍与自拍| 赤兔流量卡办理| 伊人久久国产一区二区| 欧美3d第一页| 国产黄色视频一区二区在线观看| 国模一区二区三区四区视频| 午夜亚洲福利在线播放| 内射极品少妇av片p| 日韩一区二区视频免费看| 熟妇人妻不卡中文字幕| 麻豆精品久久久久久蜜桃| 亚洲av免费高清在线观看| 国产免费视频播放在线视频 | 男人舔奶头视频| 亚洲欧美精品自产自拍| 岛国毛片在线播放| 亚洲人成网站在线观看播放| 亚洲人与动物交配视频| 久久精品综合一区二区三区| 婷婷色麻豆天堂久久| 我的老师免费观看完整版| 国产精品一区二区三区四区免费观看| 国产在线男女| 成人美女网站在线观看视频| 中文字幕免费在线视频6| 97精品久久久久久久久久精品| 99热全是精品| 久久精品国产鲁丝片午夜精品| 联通29元200g的流量卡| 黄色欧美视频在线观看| 久久鲁丝午夜福利片| 午夜福利高清视频| 又大又黄又爽视频免费| 国产精品日韩av在线免费观看| 卡戴珊不雅视频在线播放| 在线免费十八禁| 狠狠精品人妻久久久久久综合| 日韩视频在线欧美| 国产毛片a区久久久久| 三级国产精品欧美在线观看| 久久99蜜桃精品久久| 亚洲精品色激情综合| 日本免费a在线| 国产老妇伦熟女老妇高清| 如何舔出高潮| 免费av毛片视频| 久久精品久久久久久噜噜老黄| 久久久久久伊人网av| 中国国产av一级| 免费观看的影片在线观看| 欧美日韩亚洲高清精品| 国产伦精品一区二区三区四那| 黄色配什么色好看| 美女cb高潮喷水在线观看| 超碰97精品在线观看| 亚洲在线自拍视频| 国产成人精品福利久久| freevideosex欧美| 亚州av有码| 丝袜美腿在线中文| 欧美3d第一页| 午夜视频国产福利| 国产av码专区亚洲av| av国产免费在线观看| 国产精品三级大全| 最后的刺客免费高清国语| a级毛色黄片| 夜夜爽夜夜爽视频| 亚洲av二区三区四区| 只有这里有精品99| 亚洲国产欧美人成| 美女cb高潮喷水在线观看| 午夜精品国产一区二区电影 | 国产精品一二三区在线看| 免费少妇av软件| 亚洲av二区三区四区| 亚洲熟女精品中文字幕| 国产精品三级大全| 在现免费观看毛片| 午夜久久久久精精品| 人妻制服诱惑在线中文字幕| 日日啪夜夜撸| 国产国拍精品亚洲av在线观看| 一级毛片我不卡| 18禁动态无遮挡网站| 看黄色毛片网站| 精品熟女少妇av免费看| 久久精品久久久久久噜噜老黄| 国产亚洲最大av| 毛片女人毛片| 真实男女啪啪啪动态图| 网址你懂的国产日韩在线| 我的老师免费观看完整版| 亚洲国产精品专区欧美| 一夜夜www| 人妻少妇偷人精品九色| 能在线免费观看的黄片| 观看免费一级毛片| 欧美最新免费一区二区三区| 3wmmmm亚洲av在线观看| 国产精品女同一区二区软件| 亚洲成人久久爱视频| 午夜精品一区二区三区免费看| 黄色日韩在线| 国产一区二区在线观看日韩| 美女xxoo啪啪120秒动态图| 国产精品一区www在线观看| 最近最新中文字幕大全电影3| 黄色配什么色好看| 国产一区二区三区av在线| 床上黄色一级片| 国产亚洲精品久久久com| 色尼玛亚洲综合影院| 网址你懂的国产日韩在线| 日韩伦理黄色片| 啦啦啦中文免费视频观看日本| a级毛色黄片| 中文资源天堂在线| 亚洲精品日韩在线中文字幕| 亚洲内射少妇av| 国产亚洲av嫩草精品影院| 久久久久久久久久久丰满| 国产精品一区二区三区四区免费观看| 免费观看a级毛片全部| 噜噜噜噜噜久久久久久91| 久99久视频精品免费| 国产色爽女视频免费观看| 国产男女超爽视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 少妇裸体淫交视频免费看高清| 三级国产精品片| 国产伦精品一区二区三区四那| 成人毛片a级毛片在线播放| 国产男人的电影天堂91| 大香蕉97超碰在线| 乱系列少妇在线播放| 在线观看人妻少妇| 91精品一卡2卡3卡4卡| 精品国产露脸久久av麻豆 | 美女内射精品一级片tv| 国产91av在线免费观看| 女人久久www免费人成看片| 草草在线视频免费看| 亚洲精品影视一区二区三区av| 亚洲久久久久久中文字幕| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 一区二区三区免费毛片| 亚洲在线观看片| 男人舔奶头视频| 91精品伊人久久大香线蕉| 久久这里有精品视频免费| 国产精品无大码| 最近中文字幕高清免费大全6| 熟妇人妻不卡中文字幕| 午夜福利在线观看免费完整高清在| 一边亲一边摸免费视频| 男人狂女人下面高潮的视频| 亚洲经典国产精华液单| 一级a做视频免费观看| 如何舔出高潮| 国产免费福利视频在线观看| 18禁在线播放成人免费| 久久草成人影院| 亚洲国产高清在线一区二区三| 日韩av免费高清视频| 干丝袜人妻中文字幕| 国产男女超爽视频在线观看| 日韩av免费高清视频| 特大巨黑吊av在线直播| 国内精品一区二区在线观看| 丰满人妻一区二区三区视频av| 国产黄色视频一区二区在线观看| 免费大片黄手机在线观看| 26uuu在线亚洲综合色| 亚洲av免费高清在线观看| av在线亚洲专区| 99热网站在线观看| av国产久精品久网站免费入址| 国产高清三级在线| 中文欧美无线码| 免费观看的影片在线观看| 国产伦在线观看视频一区| 久久99热这里只有精品18| 国产精品一区www在线观看| 国产精品一区二区三区四区免费观看| 久久人人爽人人爽人人片va| 国产v大片淫在线免费观看| 亚洲第一区二区三区不卡| 淫秽高清视频在线观看| 欧美成人午夜免费资源| 午夜福利高清视频| 久久韩国三级中文字幕| 欧美潮喷喷水| 日本wwww免费看| 最近的中文字幕免费完整| 精品酒店卫生间| 精品久久久久久久久av| 日本三级黄在线观看| 一区二区三区高清视频在线| 美女内射精品一级片tv| 久久久欧美国产精品| 精品一区二区免费观看| 99热全是精品| 夜夜看夜夜爽夜夜摸| 亚洲aⅴ乱码一区二区在线播放| 国产成人精品一,二区| 亚洲乱码一区二区免费版| 2021少妇久久久久久久久久久| av在线亚洲专区| 亚洲熟妇中文字幕五十中出| 亚洲精品国产av蜜桃| 亚洲av国产av综合av卡| 一区二区三区免费毛片| 国产午夜精品论理片| 高清日韩中文字幕在线| 日韩 亚洲 欧美在线| 人人妻人人澡人人爽人人夜夜 | 人妻系列 视频| 永久网站在线| 一二三四中文在线观看免费高清| 国产激情偷乱视频一区二区| 久久久久久伊人网av| 一级av片app| 日本黄大片高清| 亚洲欧美日韩卡通动漫| 少妇的逼水好多| 亚洲av男天堂| 免费看日本二区| 欧美精品一区二区大全| 国产精品久久久久久av不卡| 亚洲国产日韩欧美精品在线观看| 久久久久久久久久黄片| 嘟嘟电影网在线观看| 欧美精品一区二区大全| 精品久久久久久久末码| 日韩成人伦理影院| 亚洲国产欧美在线一区| 啦啦啦啦在线视频资源| 国产有黄有色有爽视频| 国产亚洲最大av| 欧美一级a爱片免费观看看| 亚洲av成人精品一二三区| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 欧美激情在线99| 少妇的逼水好多| 国产一区亚洲一区在线观看| 亚洲美女搞黄在线观看| 91aial.com中文字幕在线观看| 简卡轻食公司| 全区人妻精品视频| 一级毛片久久久久久久久女| 欧美97在线视频| 99热网站在线观看| 亚洲国产成人一精品久久久| 亚洲精品久久久久久婷婷小说| 2021少妇久久久久久久久久久| 网址你懂的国产日韩在线| av在线天堂中文字幕| 成年女人在线观看亚洲视频 | 国产一区二区亚洲精品在线观看| av黄色大香蕉| 国产极品天堂在线| 亚洲综合精品二区| 国产成人精品福利久久| 欧美一区二区亚洲| 赤兔流量卡办理| 亚洲av福利一区| 伊人久久国产一区二区| 男人舔女人下体高潮全视频| 少妇熟女欧美另类| 成人毛片a级毛片在线播放| 国产久久久一区二区三区| 边亲边吃奶的免费视频| 禁无遮挡网站| 黄色配什么色好看| 亚洲精品国产成人久久av| 国产色婷婷99| 亚洲伊人久久精品综合| 亚洲成人一二三区av| 男的添女的下面高潮视频| 97热精品久久久久久| 日韩国内少妇激情av| 久久久久久久午夜电影| videos熟女内射| 久久精品久久久久久噜噜老黄| 三级毛片av免费| 亚洲第一区二区三区不卡| 在线观看免费高清a一片| av在线亚洲专区| 91久久精品国产一区二区三区| 成人美女网站在线观看视频| 国产精品久久久久久精品电影| 成人鲁丝片一二三区免费| 1000部很黄的大片| 国语对白做爰xxxⅹ性视频网站| 99久久九九国产精品国产免费| 免费看美女性在线毛片视频| 午夜免费男女啪啪视频观看| 中文精品一卡2卡3卡4更新| 亚洲天堂国产精品一区在线| 哪个播放器可以免费观看大片| 人妻系列 视频| 亚洲精品自拍成人| 久99久视频精品免费| 午夜福利网站1000一区二区三区| 丰满少妇做爰视频| 久久久久久久久久人人人人人人| 天天躁日日操中文字幕| 少妇熟女aⅴ在线视频| 一个人免费在线观看电影| 免费av观看视频| 国产视频首页在线观看| 波多野结衣巨乳人妻| a级毛片免费高清观看在线播放| 日韩电影二区| 天堂√8在线中文| 亚洲av成人精品一区久久| 亚洲av.av天堂| 十八禁国产超污无遮挡网站| 精品久久久久久电影网| 久久久久久久久久久丰满| 亚洲怡红院男人天堂| 国产91av在线免费观看| 26uuu在线亚洲综合色| 久久久久久久大尺度免费视频| 91aial.com中文字幕在线观看| 寂寞人妻少妇视频99o| 91久久精品国产一区二区成人| 搞女人的毛片| 一个人免费在线观看电影| 久久久久久久久久成人| 国产精品精品国产色婷婷| 看黄色毛片网站| 久久久精品欧美日韩精品| 亚洲性久久影院| av一本久久久久| 小蜜桃在线观看免费完整版高清| 久久精品综合一区二区三区| 久久久精品免费免费高清| 国产精品av视频在线免费观看| 国产国拍精品亚洲av在线观看| av在线老鸭窝| 高清日韩中文字幕在线| 十八禁国产超污无遮挡网站| 美女国产视频在线观看| 日本三级黄在线观看| 国产一区二区在线观看日韩| 精品99又大又爽又粗少妇毛片| 你懂的网址亚洲精品在线观看| 欧美日韩精品成人综合77777| 久久久久久久大尺度免费视频| 日韩精品有码人妻一区| 99久久人妻综合| 99热网站在线观看| 国产精品麻豆人妻色哟哟久久 | 成人鲁丝片一二三区免费| 久久国产乱子免费精品| 蜜桃亚洲精品一区二区三区| 日韩伦理黄色片| 又爽又黄a免费视频| 免费不卡的大黄色大毛片视频在线观看 | 欧美一级a爱片免费观看看| 精品久久久久久久久av| 国产又色又爽无遮挡免| 超碰av人人做人人爽久久| 人妻制服诱惑在线中文字幕| 欧美成人a在线观看| 你懂的网址亚洲精品在线观看| 五月天丁香电影| 国产乱人视频| 美女高潮的动态| 又大又黄又爽视频免费| 天堂影院成人在线观看| av又黄又爽大尺度在线免费看| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美| 看黄色毛片网站| 日本av手机在线免费观看| 国产黄片美女视频| 精品久久久精品久久久| 永久网站在线| 亚洲精品第二区| 日韩欧美精品v在线| 亚洲av电影不卡..在线观看| 亚洲av男天堂| 可以在线观看毛片的网站| 成人漫画全彩无遮挡| 日本熟妇午夜| 一级毛片我不卡| 国产永久视频网站| a级毛色黄片| 在线观看一区二区三区| 久久久a久久爽久久v久久| 国产精品日韩av在线免费观看| 国产伦一二天堂av在线观看| 观看美女的网站| 久久鲁丝午夜福利片| 男女视频在线观看网站免费| 精品国产露脸久久av麻豆 | 国产在线男女| 国产黄色小视频在线观看| av福利片在线观看| 老司机影院成人| 日韩av不卡免费在线播放| 国产乱来视频区| 国产真实伦视频高清在线观看| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久影院| 80岁老熟妇乱子伦牲交| 少妇裸体淫交视频免费看高清| 一个人看的www免费观看视频| 99久国产av精品国产电影| 中文天堂在线官网| 啦啦啦啦在线视频资源| 五月伊人婷婷丁香| 精品人妻一区二区三区麻豆| av又黄又爽大尺度在线免费看| 看黄色毛片网站| 女人十人毛片免费观看3o分钟| 亚洲成色77777| 99久国产av精品国产电影| 男女边摸边吃奶| 嘟嘟电影网在线观看| 一级毛片 在线播放| 免费人成在线观看视频色| 久久综合国产亚洲精品| 看十八女毛片水多多多| 国产在视频线精品| ponron亚洲| 韩国av在线不卡| 欧美bdsm另类| 丝袜喷水一区| av在线天堂中文字幕| .国产精品久久|