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

    SsC1qDC基因在毛蚶免疫防御中的作用

    2023-12-22 01:48:38王進京王春德劉相全劉桂龍寧軍號
    海洋科學(xué) 2023年9期
    關(guān)鍵詞:毛蚶弧菌補體

    彭 強, 王進京, 王春德, 盧 霞, 劉相全, 劉桂龍, 徐 鑫, 許 賀, 寧軍號

    基因在毛蚶免疫防御中的作用

    彭 強1, 2, 3, 王進京2, 王春德2, 盧 霞2, 劉相全3, 劉桂龍4, 徐 鑫4, 許 賀5, 6, 寧軍號2

    (1. 上海海洋大學(xué) 水產(chǎn)與生命學(xué)院, 上海 201306; 2. 中國科學(xué)院煙臺海岸帶研究所, 山東 煙臺 264003; 3. 山東海洋資源與環(huán)境研究所, 山東 煙臺 264006; 4. 煙臺海之春水產(chǎn)種業(yè)有限公司, 山東 煙臺 264006; 5. 江蘇寶源生物科技有限公司, 江蘇 連云港 222144; 6. 江蘇海泰海洋科技有限公司, 江蘇 連云港 222144)

    含C1q結(jié)構(gòu)域蛋白(C1qDCs)在無脊椎動物免疫應(yīng)答中發(fā)揮重要作用, 但目前有關(guān)毛蚶()C1qDCs免疫功能的研究鮮有報道。本研究從毛蚶中鑒定出一個含有典型C1q結(jié)構(gòu)域的基因(命名為), 并分析了其對主要病原感染的響應(yīng)特征。的cDNA全長序列為711 bp, 編碼含有N端信號肽和199個氨基酸的蛋白, 系統(tǒng)進化分析顯示SsC1qDC與雙殼貝類的C1qDCs聚在一起。在毛蚶卵細胞到眼點幼蟲時期的表達量隨發(fā)育階段顯著增加, 且該基因在成體毛蚶各組織中均有表達, 其中在肝胰腺和閉殼肌中表達水平較高。副溶血弧菌刺激可顯著誘導(dǎo)基因高表達, 敲降基因表達后毛蚶死亡率和血細胞凋亡水平顯著升高, 且導(dǎo)致毛蚶血細胞總數(shù)和吞噬活力顯著下降。同時, 敲降基因的表達可顯著抑制母源免疫相關(guān)基因、、、和的表達, 而、、和等基因的表達顯著上調(diào)。本研究結(jié)果表明,在毛蚶抗弧菌響應(yīng)中發(fā)揮重要作用, 對貝類免疫學(xué)和抗病新品系的選育有重要的指導(dǎo)意義。

    ; 毛蚶(); 表達圖譜; RNA干擾; 免疫防御

    補體C1q是補體C1復(fù)合蛋白的一個重要組分, 能夠激活補體經(jīng)典途徑, 該途徑是脊椎動物中連接獲得性免疫和先天免疫的重要樞紐[1, 2]。含C1q結(jié)構(gòu)域蛋白(C1q domain containing proteins, C1qDCs)是包含一個或多個球狀gC1q結(jié)構(gòu)域的蛋白超家族[1, 3], 廣泛存在于無脊椎動物和脊椎動物中[4, 5]。在脊椎動物中, C1qDCs不僅可以激活補體經(jīng)典通路, 還可以作為模式識別受體(pattern recognition receptors, PRRs), 通過結(jié)合病原表面的病原相關(guān)分子模式(pathogen-associated molecular patterns, PAMP), 快速啟動細胞吞噬等免疫過程, 有效遏制病原微生物的侵害[6]。例如, 鯉魚()的C1qDC蛋白具有細菌凝集活性及脂多糖和肽聚糖結(jié)合活性[7]。此外, 斑馬魚()C1qDC蛋白可以特異結(jié)合免疫球蛋白IgG和IgM, 并通過激活補體經(jīng)典途徑抑制血清的溶血活性[8]。

    研究發(fā)現(xiàn)文昌魚()C1q蛋白可以結(jié)合哺乳類的IgG, 但沒有直接證據(jù)表明無脊椎動物中存在補體的經(jīng)典激活途徑[9]。越來越多的研究表明, 無脊椎動物的C1qDCs參與了多種免疫應(yīng)答反應(yīng), 例如免疫識別、結(jié)合與凝集活性[10, 11], 抗菌活性[12, 13],促吞噬作用[14, 15]和細胞遷移[16-18]等。然而, C1qDCs在無脊椎動物特定免疫途徑中的主要作用鮮有報道, 需要進一步研究以豐富無脊椎動物的免疫機制。

    本研究以毛蚶為研究對象, 利用PCR和RACE技術(shù)克隆得到基因的cDNA序列。探究了基因在毛蚶不同組織、早期胚胎以及弧菌刺激后的表達特征, 并利用RNA干擾技術(shù)敲降毛蚶基因的表達, 分析其在弧菌脅迫后毛蚶血細胞吞噬和凋亡及抗弧菌感染中的功能特征。本研究為進一步闡明在毛蚶免疫防御中的潛在功能提供了理論基礎(chǔ)。

    1 材料與方法

    1.1 實驗材料

    實驗所用的毛蚶為合作公司的養(yǎng)殖群體, 均為殼形完整、健康無損傷的個體。實驗前隨機挑選30只毛蚶測量數(shù)據(jù)如下: 平均殼長為22.15±1.37 mm、平均殼高為24.17±1.08 mm、平均殼寬為19.56±0.44 mm、平均體質(zhì)量為9.63±0.71 g。在實驗室條件(pH 8.2±0.1、鹽度 23±1、溫度28±0.5 ℃)下暫養(yǎng)7 d進行實驗, 每日早晚各投喂1次螺旋藻餌料, 投喂后2 h更換全部海水。實驗所用的弧菌是副溶血弧菌()PL2菌株, 菌種由浙江省海洋水產(chǎn)養(yǎng)殖研究所清江基地實驗室提供。

    1.2 實驗方法

    1.2.1 毛蚶總RNA提取及cDNA合成

    采用Trizol提取試劑盒(Invitrogen, 美國)提取毛蚶各組織總RNA, 采用1.0%的瓊脂糖凝膠電泳檢測RNA的完整性, 利用Nanodrop 2000(Thermo Scientific, 美國)檢測RNA的純度和濃度。然后采用HiScript?III RT SuperMix反轉(zhuǎn)錄試劑盒(Vazyme, 中國)將毛蚶各組織的RNA反轉(zhuǎn)錄成cDNA, –20 ℃保存待用。

    1.2.2 毛蚶基因cDNA序列克隆

    以毛蚶的肝胰腺、血細胞和鰓等組織的cDNA為模板, 根據(jù)引物F、R和- F1,-F2(表1)擴增基因的ORF和cDNA全長, 利用PrimeSTAR?Max DNA Polymerase試劑盒(TaKaRa, 日本)進行PCR實驗, PCR反應(yīng)體系: 2.0 μL dNTP Mix、12.5 μL 10×PCR Buffer、0.5 μL正反向引物(10 μmol/L)、1.0 μL cDNA模板、0.2 μL Taq聚合酶和18.3 μL DEPC H2O(TaKaRa, 日本)。反應(yīng)程序: 94 ℃預(yù)變性 3 min, 35 個循環(huán)(94 ℃變性30 s, 54 ℃退火50 s, 72 ℃延伸60 s), 72 ℃延伸10 min。PCR產(chǎn)物經(jīng)1.0%瓊脂糖凝膠電泳切膠純化后, 將融合質(zhì)粒轉(zhuǎn)入大腸桿菌DH5α感受態(tài)細胞, 篩選陽性克隆質(zhì)粒送北京擎科生物有限公司測序。

    表1 本研究中所用的引物

    Tab.1 Primer sequences used in this study

    續(xù)表

    1.2.3 生物信息學(xué)分析

    在NCBI數(shù)據(jù)庫中BLAST模塊(http://www.ncbi. nlm.nih.gov/BLAST)進行核苷酸和氨基酸序列分析; 利用ORF Finder(http://www.ncbi.nlm.nih.gov/ gorf/gorf.htmL)在線預(yù)測基因的開放閱讀框和氨基酸序列; 使用 ExPASy(https://www.ncbi.nlm. nih.gov/orffinder/)進行理論等電點與相對分子量的預(yù)測; 使用 DNAMAN 8.0軟件對不同物種的C1qDC蛋白進行多序列比對; 利用MEGA-X軟件的鄰接法(neighbor-joining, NJ)對包括毛蚶在內(nèi)的15個C1qDC氨基酸序列構(gòu)建系統(tǒng)進化樹, bootstrap值設(shè)為1 000。

    1.2.4 不同組織、胚胎發(fā)育時期和弧菌刺激后基因表達情況

    解剖雌、雄毛蚶各5只, 分別取血細胞、性腺、閉殼肌、斧足、外套膜、肝胰腺和鰓組織等樣品, 其中毛蚶血細胞取樣具體操作如下: 先用移液槍吸取500 μL的抗凝劑(336 mmol/L NaCl, 27 mmol/L檸檬酸鈉, 115 mmol/L葡萄糖, 9 mmol/L EDTA-Na2, pH 7.0)加至無RNA酶的1.5 mL EP管中, 再用注射器或移液槍從毛蚶體內(nèi)吸取500 μL的血液, 注入提前加入抗凝劑的EP管中輕輕吹打均勻, 4 ℃、1 000 r/min離心5 min收集血細胞, 用PBS洗兩遍, 重懸于1 mmol/L的Trizol中。誘導(dǎo)毛蚶產(chǎn)卵后分別取卵細胞、卵裂期(2細胞到8細胞階段)、囊胚期、原腸期、D型幼蟲和眼點幼蟲等不同發(fā)育時期樣品, 每個發(fā)育時期取5個平行樣品。此外, 將120只健康無損的毛蚶隨機分成2組, 實驗組向海水加入副溶血弧菌溶液至終濃度為3×107CFU/mL, 對照組毛蚶在常規(guī)充氣的海水中飼養(yǎng)。在弧菌攻毒0、3、6、12、24和48 h, 隨機取5只毛蚶解剖取血細胞。上述所有樣品均在液氮中速凍后置于–80℃冰箱保存、待用。

    依據(jù)毛蚶基因cDNA序列設(shè)計熒光定量引物(表1), 利用 qRT-PCR技術(shù)檢測基因在毛蚶不同組織、胚胎和幼蟲發(fā)育時期及弧菌刺激后的表達模式, 實驗設(shè)置4個不同樣本重復(fù), 內(nèi)參基因為18S rRNA。利用Taq Pro Universal SYBR qPCR Master Mix試劑盒(Vazyme, 中國)進行qRT-PCR實驗, 反應(yīng)體系: 2×Taq Pro Universal SYBR qPCR Master Mix 5 μL、ddH2O 3.6 μL、cDNA模板1 μL、正、反向引物各0.2 μL。反應(yīng)程序: 95 ℃, 30 s; 40個循環(huán)(95 ℃, 5 s; 60 ℃, 35 s)。反應(yīng)在ABI QuantStudioTM 5 Real-Time PCR Instrument(Thermo Scientific, 美國)儀器中進行, 利用2–ΔΔCt法分析基因的相對表達水平[19]。

    1.2.5 RNA干擾分析

    根據(jù)基因的ORF序列設(shè)計其小干擾RNA序列為-siRNA(5′-CGACAAGGUAGUA AUGUAU-3′), 將該序列隨機打亂順序獲得- random-siRNA(5′-AGACUAAUGGAUACUAUGG-3′)序列, 利用T7啟動子體外轉(zhuǎn)錄試劑盒(Tiandz, 天恩澤, 中國)并按其說明書進行小RNA合成與純化。

    將120只健康無損的毛蚶隨機分為實驗組和對照組, 每組各60只, 利用微量注射器分別注射6 μL的-siRNAs和-random-siRNAs至毛蚶的體腔中。在注射后的0、12、24、48和72 h, 兩組分別隨機取5只毛蚶解剖取血細胞樣品, 液氮凍存后轉(zhuǎn)入–80℃冰箱保存。經(jīng)RNA提取和反轉(zhuǎn)錄成cDNA后, 通過qRT-PCR檢測實驗組和對照組毛蚶基因的表達量, 分析-siRNA的干擾效率。

    1.2.6 敲降毛蚶基因后母源性免疫基因的表達分析

    敲降基因的表達后, 通過qRT-PCR技術(shù)檢測母源免疫相關(guān)基因的表達情況, 這些基因包括卵黃蛋白原(vitellogenin,)、細胞黏附分子免疫球蛋白(down syndrome cell adhesion molecule,)、補體分子(complement,)、硫酯蛋白(thioester- containing protein,)、溶菌酶(lysozyme,)、酚氧化酶原激活因子3(phenoloxidase-activating factor 3,)以及同源基因和。用于檢測母源性免疫基因表達的特異性引物見表1, 18S rRNA作為內(nèi)參基因。

    1.2.7 敲降基因的毛蚶弧菌刺激后血細胞凋亡、吞噬活性、血細胞總數(shù)和存活率的變化

    將120只健康無損的毛蚶隨機分為2組, 分別注射-siRNA和-random-siRNA后12 h, 分別向兩組毛蚶養(yǎng)殖水體加入終濃度為3×107CFU/mL的副溶血弧菌, 弧菌攻毒后的0、3、6、12、24、48和72 h取樣, 每個時間點實驗組和對照組各隨機解剖5只毛蚶取血淋巴樣品。用預(yù)裝500 μL抗凝劑的無菌注射器收集500 μL血淋巴, 并用300目篩絹過濾。4 ℃、1 000 r/min離心5 min后倒去上清液, 重懸于1 mL的PBS緩沖液, 相同條件再次離心倒掉上清液, 重懸于1mL的PBS緩沖液中, 用制備好的血細胞工作液檢測血細胞凋亡、吞噬活性、血細胞總數(shù)。

    利用Annexin V-FITC凋亡試劑盒(碧云天, 中國)檢測血細胞凋亡水平, 操作如下: 取200 μL血細胞工作液避光條件下加入195 μL的FITC結(jié)合液、5 μL FITC和10 μL碘化丙染色液(PI), 輕輕混勻, 25℃避光孵育20 min后用流式細胞儀(CytoFLEX, 美國)檢測細胞凋亡率。設(shè)置僅添加FITC或PI的單一染料對照組和不添加染料的陰性對照組, 來確定散射圖十字門的位置(圖1)。毛蚶血細胞的凋亡率用FITC陽性和PI陰性細胞的數(shù)量占細胞總數(shù)的百分比來表示[20]。

    圖1 細胞凋亡散點圖中十字門的位置

    注: 未添加染料的陰性對照組(a)和僅添加FITC(b)或PI(c)的單一染料對照組, 來確定散點圖十字門的位置

    血細胞吞噬率檢測方法參考文獻[21]并進行相應(yīng)調(diào)整: 將200 mL血細胞工作液與20 μL 3.0%熒光微球(YG 1.0 μm, Polysciences, 美國)混勻, 在25 ℃黑暗環(huán)境中孵育1 h。后用10 μL 37%福爾馬林溶液(上海吉至生化科技有限公司)固定, 并在1 000 r/min, 4 ℃條件下離心10 min, 重懸于200 mL PBS緩沖液中并通過流式細胞儀(FL-1通道)檢測, 血細胞的吞噬率用吞噬微球的血細胞數(shù)目占血細胞總數(shù)的百分比表示(圖2)。

    使用 SYBR Green I 熒光染料(Biosharp, 中國)對血細胞染色, SYBR Green I可以通過和細胞核中的DNA結(jié)合, 產(chǎn)生綠色熒光。取180 μL血細胞工作液, 加入20 μL 10×的SYBR Green I 染料, 25 ℃避光孵育1 h, 孵育完成后用移液槍吸取50 μL反應(yīng)液, 利用流式細胞儀(FL-1通道)檢測綠色熒光的細胞并進行計數(shù)。

    180只健康的毛蚶被隨機分成-siRNA和-random-siRNA注射后加弧菌組以及-siRNA注射后加PBS 3組, 每組設(shè)置3個平行, 每個平行20只毛蚶。在siRNA注射后12 h,-siRNA和-random-siRNA注射組的毛蚶在28℃的含有副溶血弧菌(3×107CFU/ mL)的海水中養(yǎng)殖。-siRNA注射組中加入等量的 PBS緩沖液(Solarbio, 中國)作為空白對照。每12 h統(tǒng)計3組毛蚶的存活個數(shù)。

    1.2.8 數(shù)據(jù)分析

    實驗數(shù)據(jù)結(jié)果采用SPSS 26.0軟件進行統(tǒng)計分析, 數(shù)據(jù)以平均值±標(biāo)準(zhǔn)差(Mean±SD)表示, 單因素方差分析(One-Way ANOVA)比較不同數(shù)據(jù)組間的差異, 其中<0.05表示差異顯著,<0.01表示差異極顯著。

    2 結(jié)果與分析

    2.1 毛蚶SsC1qDC基因的克隆和序列特點

    利用PCR和RACE克隆技術(shù)得到毛蚶基因(基因序列號: OQ264768)的cDNA全長序列為711 bp(圖3), 其中5′-UTR 21 bp, 3′-UTR 80 bp, ORF 600 bp, 共編碼199個氨基酸。的C1q結(jié)構(gòu)域為第66~199個氨基酸, 預(yù)測相對分子量22.64 kDa, 理論等電點為8.60。

    圖2 通過流式細胞術(shù)檢測得到的毛蚶血細胞吞噬率的FITC直方圖

    a. 不加微球的陰性對照; b. 加微球的實驗組

    圖3 毛蚶SsC1qDC核苷酸及預(yù)測的氨基酸序列

    注: 起始密碼子(ATG)、信號肽、Clq結(jié)構(gòu)域和終止密碼子(TAG)分別用綠色、紅色和黃色以及灰色陰影標(biāo)出

    2.2 序列比對與系統(tǒng)進化關(guān)系

    將毛蚶SsC1qDC保守的C1q結(jié)構(gòu)域序列與其他物種的C1q結(jié)構(gòu)域序列對比, 結(jié)果如圖4所示, 毛蚶SsC1qDC與雙殼類SsC1qDCs的同源性高于脊椎動物。毛蚶SsC1qDC與美洲牡蠣(41.23%)CvC1qDC、硬殼蛤(39.47%)MmC1qDC的相似性較高, 與人()和家鼠()SsC1qDCs的相似性較低, 分別是28.83%、26.61%。

    利用15個不同動物門的C1qDCs結(jié)構(gòu)域序列構(gòu)建NJ 進化樹, 系統(tǒng)進化樹如圖5所示, 進化樹的拓?fù)浣Y(jié)構(gòu)明顯分無脊椎動物和脊椎動物2大分支。毛蚶SsC1qDC與紫貽貝()、硬殼蛤()C1qDCs聚為一支, 然后與其他貝類的C1qDCs形成無脊椎動物分支; 來自哺乳動物、鳥類、爬行動物、兩棲動物和魚類的C1qDCs聚為脊椎動物分支。

    圖4 C1qDCs蛋白質(zhì)的多序列比對

    注: C1qDCs的物種和登錄號如下: 人(NP_852100.3); 小鼠(NP_112150.1); 長牡蠣(XP_034333860.1); 美洲牡蠣(XP_022291269.1); 紅鮑(XP_046377716.2); 歐洲平牡蠣(CAG2203397.1); 亞馬遜花鳉(XP_007577129.1); 非洲爪蟾(XP_018122782.1); 硬殼蛤(XP_045158489.1); 貝氏隆頭魚(XP_020504487.1)

    2.3 毛蚶SsC1qDC基因在不同組織和胚胎時期的表達特征

    如圖6a所示,基因在毛蚶所有檢測組織中均有表達, 在肝胰腺、閉殼肌和鰓中表達水平較高,在雄性毛蚶性腺、斧足、外套膜和鰓中的表達量顯著高于雌性毛蚶, 而在肝胰腺中的表達量呈相反趨勢(圖6a所示)?;蛟谂咛グl(fā)育階段的表達模式如圖6b所示, 與卵細胞期相比,的表達從卵裂期到眼點幼蟲階段均顯著上調(diào)(<0.05)。

    2.4 毛蚶SsC1qDC基因在副溶血弧菌刺激后的表達模式

    副溶血弧菌刺激后, 毛蚶血細胞中基因的表達量如圖7所示, 與對照組相比,的表達量在3、6、24和48 h顯著升高(<0.01), 呈現(xiàn)先升高再下降而后再升高的變化趨勢, 并在6 h達到最大值(約為對照組的4.0倍)。隨后的表達水平顯著下降, 在24 h達到表達量最低, 顯著低于對照組的表達量(<0.01)。在48 h,的表達量再次顯著上調(diào)(<0.01)。

    注: 利用MEGA-X軟件, 采用鄰接法構(gòu)建系統(tǒng)進化樹, 重復(fù)1 000次。C1qDCs的物種和登錄號如下: 人(NP_852100.3); 小鼠(NP_112150.1); 蝙蝠(XP_0244 25520.1); 非洲爪蟾(XP_018122782.1); 斑馬魚(XP_ 009296972.1); 蝦夷扇貝(XP_021374682.1); 灘棲螺(KAG5686619.1); 太平洋牡蠣(XP_034333860.1); 長牡蠣(XP_022291269.1); 紅鮑(XP_046377716.2); 歐洲平牡蠣(CAG2203397.1); 紫貽貝(CAG2203397.1); 硬殼蛤(XP_045158489.1)

    圖7 弧菌刺激后SsC1qDC的表達情況(平均值±標(biāo)準(zhǔn)差, n=3)

    2.5 毛蚶SsC1qDC基因siRNA的干擾效率

    合成的-siRNA的干擾效率見圖8。與注射隨機打亂序列的siRNA對照組相比, 毛蚶注射-siRNA后12和24 h血細胞中基因干擾效果顯著(<0.01), 干擾效率分別為64.9%和59.9%。因此, 選擇注射-siRNA后12 h的血細胞用于后續(xù)實驗分析。

    2.6 敲降SsC1qDC基因后免疫相關(guān)基因的表達情況

    敲降基因的表達后, 母源性免疫基因的表達顯著下調(diào)(<0.01), 而、、和的表達顯著上調(diào)(<0.01), 其中和基因的表達量約是對照組的兩倍(圖9)。

    圖8 注射siRNA對SsC1qDC基因的干擾效率(平均值±標(biāo)準(zhǔn)差, n=3)

    圖9 敲降SsC1qDC后免疫相關(guān)基因的表達情況(平均值±標(biāo)準(zhǔn)差, n=3)

    2.7 敲降SsC1qDC基因表達對弧菌刺激后毛蚶的血細胞凋亡、吞噬活性、血細胞總數(shù)及存活率的影響

    如圖10a所示, 敲降基因的表達, 副溶血弧菌刺激后實驗組毛蚶血細胞的凋亡率在6 h顯著高于對照組(<0.05), 并維持較高的凋亡水平到72 h; 實驗組毛蚶血細胞的吞噬率在12和24 h顯著低于對照組(<0.05), 與弧菌刺激前相比, 48和72 h的實驗組和對照組毛蚶血細胞吞噬活力均明顯下降(圖10b); 實驗組毛蚶血細胞總數(shù)在3 h后顯著低于對照組, 與弧菌刺激前相比, 6 h后實驗組和對照組毛蚶血細胞總數(shù)均明顯下降(圖10c); 副溶血弧菌脅迫后實驗組的毛蚶的死亡率在24 h顯著高于對照組, 且其死亡率在36 h超過60%, 在84 h全部死亡, 而對照組毛蚶在36 h死亡率不到20%, 在84 h仍有超過20%的毛蚶存活(圖10d)。

    3 討論

    越來越多的研究發(fā)現(xiàn)雙殼貝類基因具有多種免疫功能, 在對不同病原的免疫防御反應(yīng)中起重要作用[7, 9, 11, 14, 22]。本研究從毛蚶中克隆獲得一個新的基因(命名為), 多序列比對顯示該基因編碼的蛋白不保守, 這與水生動物中C1qDC蛋白序列高度可變性的研究結(jié)果一致[14, 23]。系統(tǒng)進化分析顯示與紫貽貝、硬殼蛤的C1qDC蛋白聚為一支, 然后與其他軟體動物C1qDC聚成一大支, 表明SsC1qDC是軟體動物C1qDC家族中的新成員。

    基因在毛蚶所有檢測的組織中均有表達, 尤其在肝胰腺和鰓等免疫組織中明顯高表達, 提示該基因在毛蚶先天免疫中的潛在作用。已有的研究表明,基因在雙殼貝類不同組織中普遍表達[23-26], 通常在血細胞中高表達, 然而本研究發(fā)現(xiàn)在毛蚶血細胞中表達量較低, 表明基因在不同物種組織中存在表達差異。毛蚶SsC1qDC蛋白序列存在信號肽, 推測SsC1qDC是一種分泌型蛋白, 可通過開放式循環(huán)系統(tǒng)轉(zhuǎn)運至機體的不同組織, 以發(fā)揮其免疫功能[17]。有研究發(fā)現(xiàn)成體三疣梭子蟹()的母源性免疫基因在卵巢中的表達通常顯著高于精巢[27, 28]。本研究中,僅在雌性毛蚶肝胰腺中的表達水平顯著高于雄性毛蚶中的表達, 但在雄性鰓、性腺、足和外套膜等組織中的表達顯著高于雌性, 這一結(jié)果表明可能不是母源性免疫基因。因為母源的mRNAs通常在早期胚胎發(fā)育中因為消耗而顯著下調(diào)表達, 待合子基因組激活后通過轉(zhuǎn)錄而顯著增加其表達量[29-31],從受精卵到D形幼蟲階段的顯著上調(diào)表達也暗示不是母源基因。

    副溶血弧菌刺激后, 毛蚶的表達量在短時間內(nèi)顯著升高, 并且在6 h達到最高值, 提示SsC1qDC是一種可誘導(dǎo)的快速免疫效應(yīng)分子, 這與馬氏珠母貝()感染溶藻弧菌()后的表達模式一致[32]。通常在清除入侵的病原或外來刺激物時, C1qDCs的mRNA水平往往呈現(xiàn)先升高后降低的趨勢[23, 31-34]。在本研究中, 副溶血弧菌刺激后毛蚶基因表達呈先升高后降低而后再升高的變化趨勢,基因在6 h后表達降低, 可能是負(fù)反饋調(diào)節(jié)機制、抑制因子的產(chǎn)生以及清除被激活的細胞等調(diào)控作用造成基因表達下降[35, 36]。由于刺激仍然存在, 免疫細胞可能會形成免疫記憶[37], 以便更好地應(yīng)對未來的刺激, 造成48 h后基因表達的再次升高, 暗示該基因在毛蚶血液系統(tǒng)中起到免疫功能調(diào)節(jié)的作用。此外, 敲降的表達, 毛蚶在副溶血弧菌刺激后的死亡率顯著高于對照組, 這一結(jié)果支持上述參與毛蚶抗弧菌反應(yīng)的觀點。

    圖10 敲降SsC1qDC基因表達對弧菌刺激后毛蚶的血細胞凋亡率(a)、吞噬率(b)、血細胞總數(shù)(c)和毛蚶存活率(d)的影響(平均值±標(biāo)準(zhǔn)差, n=3)

    在脊椎動物中, 補體經(jīng)典途徑由補體C1q與抗體-抗原復(fù)合物的相互作用激活并觸發(fā)C3的蛋白水解, 進而清除入侵的病原[3, 11]。本研究中, 敲降基因可以影響補體C3樣基因(如和)的表達, 這與脊椎動物中補體經(jīng)典途徑的激活的調(diào)控過程類似, 提示補體的經(jīng)典激活途徑是動物由低等向高等進化過程中逐漸完善而來的。通常雙殼類動物缺乏免疫球蛋白, 因此缺少經(jīng)典的補體途徑, 但櫛孔扇貝()的CfC1qDC蛋白具有與免疫球蛋白IgG的結(jié)合活性[14], 推測無脊椎動物C1qDC蛋白可能通過凝集素激活途徑激活補體系統(tǒng)[38]。

    本研究中, 注射-siRNA后12 h對基因的干擾效果最佳, 12 h后干擾效果下降, 可能與-siRNA在毛蚶體內(nèi)的濃度降低有關(guān), 這是因為-siRNA的生物穩(wěn)定性較差, 可能被毛蚶體內(nèi)的核酸酶分解[39]。此外, 由于細胞更新和分裂, 注射的-siRNA會逐漸稀釋, 從而降低對基因的干擾效率[40]。因此注射- siRNA后的干擾效果只能在一定時間內(nèi)觀察到, 隨后干擾效果會逐漸降低。和是典型的母源性免疫基因, 在水生動物的早期發(fā)育和入侵病原快速清除中發(fā)揮重要作用[17, 28, 41, 42]。本研究中, 除補體外, 敲降基因表達后, 上述免疫基因的表達均顯著下調(diào), 表明基因可能參與調(diào)控毛蚶基礎(chǔ)免疫。這與在三疣梭子蟹中的研究結(jié)果相一致, 當(dāng)敲降基因表達后, 補體樣基因(和)的表達顯著上調(diào)[29]。

    吞噬作用作為免疫細胞的第一道防線, 具有識別、吞噬和清除外源異物或自身凋亡細胞的作用[43]。先前的研究表明, C1qDCs可以作為調(diào)理素促進后生動物的吞噬作用[12, 32, 44]。本研究中, 敲降毛蚶基因的表達后, 副溶血弧菌刺激后血細胞的吞噬率顯著下降, 提示可能參與調(diào)節(jié)毛蚶血細胞的吞噬作用。這與馬氏珠母貝中敲降-基因表達能明顯降低血細胞吞噬作用的結(jié)果一致[23]。同時, 已有文獻報道雙殼類的C1qDC蛋白可參與病原的識別, 從而直接增強其血細胞的吞噬作用[12, 17, 22]。

    在軟體動物的免疫反應(yīng)中, 細胞凋亡通過應(yīng)激信號或吞噬作用清除受損、衰老和受感染的細胞, 并引發(fā)凋亡細胞死亡[45]。本研究中, 敲降基因表達可顯著提升毛蚶血細胞的凋亡水平, 從而導(dǎo)致血細胞總數(shù)和吞噬活力的顯著降低。此外, 當(dāng)外界脅迫超過動物的耐受極限后, 隨著脅迫的持續(xù), 脅迫損傷變得不可修復(fù), 最終導(dǎo)致動物死亡[46]。本研究發(fā)現(xiàn), 敲降表達, 副溶血弧菌刺激后毛蚶的死亡率顯著高于對照組, 表明在毛蚶抗弧菌感染中具有關(guān)鍵作用。

    在三疣梭子蟹中的研究中發(fā)現(xiàn), PtgC1qR可在蛋白和轉(zhuǎn)錄水平上抑制其酚氧化酶原激活系統(tǒng)[29]。在本研究中, 敲降基因顯著增強了的表達, 提示基因可能通過酚氧化酶原激活系統(tǒng)參與毛蚶的體液免疫。有關(guān)毛蚶與下游調(diào)控基因及效應(yīng)分子之間的調(diào)控關(guān)系有待于進一步研究。

    4 結(jié)論

    綜上, 本研究克隆獲得了毛蚶基因的cDNA全長序列, 并分析了該基因在響應(yīng)主要病原感染后的功能特征。的表達量從毛蚶卵細胞到眼點幼蟲期隨發(fā)育階段顯著增加, 該基因在成體毛蚶各組織廣泛表達且可被副溶血弧菌誘導(dǎo)高表達。敲降基因表達后, 母源免疫相關(guān)基因的表達受到顯著影響。此外, 與對照組相比, 副溶血弧菌攻毒后基因干擾實驗組的毛蚶死亡率和血細胞凋亡水平顯著升高, 因此導(dǎo)致其血細胞總數(shù)和吞噬活力顯著降低。結(jié)果表明,基因可能參與了毛蚶抗弧菌感染免疫防御反應(yīng)。

    [1] KISHORE U, REID K B M. C1q: structure, function, and receptors[J]. Immunopharmacology, 2000, 49(1/2): 159-170.

    [2] GABORIAUD C, THIELENS N M, GREGORY L A, et al. Structure and activation of the C1 complex of complement: unraveling the puzzle[J]. Trends in Immunology, 2004, 25(7): 368-373.

    [3] GHAI R, WATERS P, ROUMENINA L T, et al. C1q and its growing family[J]. Immunobiology, 2007, 212(4-5): 253-266.

    [4] CARLAND T M, GERWICK L. The C1q domain containing proteins: Where do they come from and what do they do? [J]. Developmental & Comparative Immunology, 2010, 34(8): 785-790.

    [5] KISHORE U, GHAI R, GREENHOUGH T J, et al. Structural and functional anatomy of the globular domain of complement protein C1q[J]. Immunology Letters, 2004, 95(2): 113-128.

    [6] ADZIGBLI L, HAO R, JIAO Y, et al. Immune response of pearl oysters to stress and diseases[J]. Reviews in Aquaculture, 2020, 12(2): 513-523.

    [7] ZHU L, GAO Y Y, YUAN G L, et al. 2021. Identification and immune functional analysis of a novel C1q domain-containing protein in Yellow River Carp (L.) [J]. Aquaculture, 534: 736313.

    [8] HU Y L, PAN X M, XIANG L X, et al. Characterization of C1q in teleosts: insight into the molecular and functional evolution of C1q family and classical pathway[J]. The Journal of Biological Chemistry, 2010, 285(37): 28777-28786.

    [9] GAO Z, LI M Y, MA J, et al. AngC1q protein binds human IgG and initiates the classical pathway: implications for a C1q-mediated complement system in the basal chordate[J]. European Journal of Immunology, 2014, 44(12): 3680-3695.

    [10] MATSUSHITA M, MATSUSHITA A, ENDO Y, et al. Origin of the classical complement pathway: lamprey orthologue of mammalian C1q acts as a lectin[J]. Proceedings of the National Academy of Sciences, 2004, 101(27): 10127-10131.

    [11] WANG L L, WANG L L, ZHANG H, et al. A C1q domain containing protein from scallopserving as pattern recognition receptor with heat-aggregated IgG binding activity[J]. PLoS One, 2012, 7 (8): e43289.

    [12] LV Z, WANG L L, JIA Z H, et al. Hemolymph C1qDC promotes the phagocytosis of oysterhemocytes by interacting with the membrane receptor β-integrin[J]. Developmental & Comparative Immunology, 2019, 98: 42-53.

    [13] LIANG X R, XIONG X W, CAO Y F, et al. Globular C1q domain-containing protein frommartensii participates in the immune defense process[J]. Fish &Shellfish Immunology, 2022, 123: 521-527.

    [14] BOHLSON S S, RASER D A, TENNER A J. Complement proteins C1q and MBL are pattern recognition molecules that signal immediate and long-term protective immune functions[J]. Molecular Immunology, 2007, 44(1-3): 33-43.

    [15] JIANG S, Li H, ZHANG D X, et al. A C1q domain containing protein fromserves as pattern recognition receptor and opsonin with high binding affinity to LPS[J]. Fish & Shellfish Immunology, 2015, 45(2): 583-591.

    [16] TAHTOUH M, CROP F, VIZIOLI J, et al. Evidence for a novel chemotactic C1q domain-containing factor in the leech nerve cord[J]. Molecular Immunology, 2009, 46: 523-531.

    [17] LI D, WAN Z, LI X J, et al. Alternatively spliced down syndrome cell adhesion molecμLe (Dscam) controls innate immunity in crab[J]. Journal of Biological Chemistry, 2019, 294(44): 16440-16450.

    [18] MORISS K M, ADEN D P, KNOELES B B, et al. Complement biosynthesis by the human hepatoma- derived cell line HepG2[J]. The Journal of clinical investigation, 1982, 70(4): 906-913.

    [19] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2?ΔΔCTmethod[J]. Methods, 2001, 25(4): 402-408.

    [20] SCHUTTE B, NUYDENS R, GEERTS H, et al. Annexin V binding assay as a tool to measure apoptosis in differentiated neuronal cells[J]. Journal of Neuroscience Methods, 1998, 86(1): 63-69.

    [21] DELAPORTE M, SOUDANT P, MOAL J, et al. Effect of a mono-specific algal diet on immune functions in two bivalve species-and[J]. The Journal of Experimental Biology, 2003, 206(17): 3053-3064.

    [22] JIANG S, LI H, ZHANG D, et al. A C1q domain containing protein fromserves as pattern recognition receptor and opsonin with high binding affinity to LPS[J]. Fish & Shellfish Immunology, 2015, 45(2): 583-591.

    [23] WANG J, THAIMUANGPHOL W, CHEN Z, et al. A C1q domain-containing protein incontributes to the innate immune response and elimination of the pathogen[J]. Fish & Shellfish Immunology, 2022, 131: 582-589.

    [24] GERDOL M, MANFRIN C, DE MORO G, et al. The C1q domain containing proteins of the Mediterranean mussel: a widespread and diverse family of immune-related molecules[J]. Developmental and Comparative Immunology, 2011, 35(6): 635-643.

    [25] GERDOL M, GRECO S, PALLAVICINI A. Extensive tandem duplication events drive the expansion of the C1q- domain-containing gene family in bivalves[J]. Marine Drugs, 2019, 17(10): 583.

    [26] GERDOL M, VENIER P, PALLAVICINI A. The genome of the Pacific oysterbrings new insights on the massive expansion of the C1q gene family in bivalvia[J]. Developmental & Comparative Immunology, 2015, 49(1): 59-71.

    [27] NING J H, LIU Y, GAO F T, et al. 2019. Characterization and functional analysis of a novel gC1qR in the swimming crab[J]. Fish & Shellfish Immunology, 2019, 84: 970-978.

    [28] NING J H, LIU Y, CUI Z X. Identification and functional analysis of a thioester-containing protein fromreveals its involvement in the prophenoloxidase system, phagocytosis and AMP synthesis[J]. Aquaculture, 2019, 510: 9-21.

    [29] 寧軍號. 三疣梭子蟹補體樣分子結(jié)構(gòu)及免疫功能研究[D]. 青島: 中國科學(xué)院海洋研究所, 2019.NING Junhao. Study on complement-like molecular structure and immune function of[D]. Qingdao: Institute of Oceanography, Chinese Academy of Sciences, 2019.

    [30] TADROS W, LIPSHITZ H D. The maternal-to-zygotic transition: a play in two acts[J]. Development, 2009, 136: 3033-3042.

    [31] WALSER C B, LIPSHITZ H D. Transcript clearance during the maternal-to-zygotic transition[J]. Current Opinion in Genetics and Development, 2011, 21(4): 431-443.

    [32] WANG Z L, WANG B, CHEN G et al. An alpha-2 macroglobμLin in the pearl oyster: characterization and function in hemocyte phagocytosis of[J]. Fish & Shellfish Immunology, 2016, 55: 585-594.

    [33] ZHU L, GAO Y, YUAN G et al. Identification and immune functional Analysis of an novel C1q domain- containing protein in Yellow River carp (L.) [J]. Aquaculture, 2021, 534: 736313.

    [34] DU X, WANG G H, YUE B, et al. A novel C1q domain containing protein in black rockfish () serves as a pattern recognition receptor with immunore-gulatory properties and possesses binding activity to heat- aggregated IgG[J]. Fish & Shellfish Immunology, 2019, 87: 73-81.

    [35] PERRIMON N, MCMAHON A P. Negative feedback mechanisms and their roles during pattern formation[J]. Cell, 1999, 97(1): 13-16.

    [36] SANCHEZ-MARTIN A, SANCHON-SANCHEZ P, ROMERO M R, et al. Impact of tumor suppressor genes inactivation on the multidrug resistance phenotype of hepatocellular carcinoma cells[J]. Biomedicine & Pharmacotherapy, 2023, 165: 115209.

    [37] WERONIKA R, PAULINA N, BEATA T, et al. Immunological memory cells[J]. Central-European Journal of Immunology, 2018, 43(2): 194-203.

    [38] ENDO Y, TAKAHASHI M, FUJITA T. Lectin complement system and pattern recognition[J]. Immunobiolo-gy, 2006, 211(4): 283-293.

    [39] LIEBERMAN J, SONG E, LEE S, et al. Interfering with disease: opportunities and roadblocks to harness-sing RNA interference[J]. Trends in Molecular Medicine, 2003, 9(9): 397-403.

    [40] CASTANOTTO D, ROSSI J J. The promises and pitfalls of RNA-interference-based therapeutics[J]. Nature, 2009, 457(7228): 426-433.

    [41] MATSUMOTO T, YAMANO K, KITAMURA M, et al. Ovarian follicle cells are the site of vitellogenin synthesis in the Pacific abalone[J]. Comparative Biochemistry and Physiology. Part A, Molecular & Integrative Physiology, 2008, 149(3): 293-298.

    [42] WANG Z P, ZHANG S C, TONG Z, et al. Maternal transfer and protective role of the alternative complement components in zebrafish[J]. PLoS One, 2009, 4(2): e4498.

    [43] STUART L M, EZEKOWITZ R A. Phagocytosis and comparative innate immunity: learning on the fly[J]. Nature Reviews Immunology, 2008, 8(2): 131-141.

    [44] NEPOMUCENO R R, HENSCHEN-EDMAN A H, BURGESS W H, et al. cDNA cloning and primary structure analysis of C1qRP, the human C1q/MBL/SPA receptor that mediates enhanced phagocytosis in[J]. Immunity, 1997, 6(2): 119-129.

    [45] TERAHARA K, TAKAHASHI K G. Mechanisms and immunological roles of apoptosis in molluscs[J]. Current Pharmaceutical Design, 2008, 14(2): 131-137.

    [46] YANG C Y, GAO Q, LIU C, et al. The transcriptional response of the Pacific oysteragainst acute heat stress[J]. Fish & Shellfish Immunology, 2017, 68: 132-143.

    A novel C1q domain-containing gene () fromparticipates in the immune defense against pathogen infection

    PENG Qiang1, 2, 3, WANG Jin-jing2, WANG Chun-de2, LU Xia2, LIU Xiang-quan3, LIU Gui-long4, XU Xin4, XU He5, 6, NING Jun-hao2

    (1. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306, China; 2. Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai 264003, China; 3. Shandong Marine Resource and Environment Research Institute, Yantai 264006, China; 4. Yantai Spring-Sea AquaSeed Co., Ltd., Yantai 264006, China; 5. Jiangsu Baoyuan Biotechnology Co., Ltd., Lianyungang 222144, China; 6. Jiangsu Haitai MariTech Co., Ltd., Lianyungang 222144, China)

    The C1q domain-containing proteins (C1qDCs) have been demonstrated to play crucial roles in the immune responses of invertebrates, with very few reports in the ark shells. In the present study, a C1qDC gene with a typical C1q domain (designated as) was identified from the ark shelland its possible immune responses against pathogen infection was also characterized. The full-length cDNA sequence ofwas 711 bp, encoding a peptide of 199 amino acids with an N-terminal signal peptide. SsC1qDC was clustered with the C1qDCs from bivalves in the phylogenetic tree.was widely expressed in the tissues, with the high expression level in hepatopancreas and adductor muscles. The expression ofwas increased significantly from eggs to eye-point larvae stage and induced remarkably byinfection. Moreover, compared to the controls, knockdown ofresulted in a remarkable increase of mortality and apoptotic levels of ark shells afterchallenge, and thus a significant reduction of total number and phagocytosis of hemocytes. Meanwhile, knockdown ofcould significantly inhibit the expression of maternally derived immune-related genes (,,,and), while the expression of,,andwas up-regulated remarkably. Taken together, these results suggested thatperforms crucial roles in immune defense against pathogen infection in ark shells, which is of great guiding significance to the immunology of shellfish and the breeding of new disease-resistant strains.

    ;; expression profiles; RNA interference; immune defense

    May 18, 2023

    [the National Natural Science Foundation of China, No. 31972791; The Earmarked Fund for Agriculture Seed Improvement Project of Shandong Province, No. 2020LZGC016; The Scientific and Technological Project of Yantai, Shandong Province, No. 2022XCZX083; The Earmarked Fund for Shandong Modern Agro-Industry Technology Research System, No. SDAIT-14]

    S917

    A

    1000-3096(2023)9-0040-13

    10.11759/hykx20230518001

    2023-05-18;

    2023-08-06

    國家自然科學(xué)基金資助項目(31972791); 山東省農(nóng)業(yè)良種工程資助項目(2020LZGC016); 山東省煙臺市科技項目專項基金資助項目(2022XCZX083); 山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)項目(SDAIT-14)

    彭強(1995—), 男, 河北保定人, 碩士, 主要從事貝類遺傳育種研究, E-mail: pq16603173759@163.com; 寧軍號(1988—), 男, 通信作者, 助理研究員, 主要從事貝類遺傳育種研究, E-mail: jhning@yic.ac.cn

    (本文編輯: 譚雪靜)

    猜你喜歡
    毛蚶弧菌補體
    “補體法”在立體幾何解題中的妙用
    銷量增長200倍!“弧菌克星”風(fēng)靡行業(yè),3天殺滅98%弧菌
    水溫、底質(zhì)對不同規(guī)格毛蚶潛砂行為的影響
    毛蚶增殖放流技術(shù)
    副溶血弧菌檢測方法的研究進展
    補體因子H與心血管疾病的研究進展
    抗dsDNA抗體、補體C3及其他實驗室指標(biāo)對于診斷系統(tǒng)性紅斑狼瘡腎損傷的臨床意義
    2種升溫方式對毛蚶呼吸代謝及免疫相關(guān)酶活性的影響
    如何有效防控對蝦養(yǎng)殖中的弧菌病
    流感患兒血清免疫球蛋白及補體的檢測意義
    午夜福利在线观看免费完整高清在| 只有这里有精品99| 小蜜桃在线观看免费完整版高清| 国产91av在线免费观看| 亚洲国产成人一精品久久久| 人妻少妇偷人精品九色| 国产精品嫩草影院av在线观看| 91久久精品国产一区二区成人| 欧美一区二区亚洲| 变态另类丝袜制服| 成年女人永久免费观看视频| 长腿黑丝高跟| 精品久久国产蜜桃| 国产精品久久久久久精品电影| 久久国产乱子免费精品| 99久久成人亚洲精品观看| 热99re8久久精品国产| 亚洲欧美中文字幕日韩二区| 99久久精品国产国产毛片| 免费观看的影片在线观看| 亚洲四区av| 男人的好看免费观看在线视频| 国产精品久久久久久av不卡| 99久久无色码亚洲精品果冻| 两个人视频免费观看高清| 一区二区三区四区激情视频| 男女国产视频网站| a级毛片免费高清观看在线播放| 精品午夜福利在线看| 久久久久久久久久黄片| 长腿黑丝高跟| 国产亚洲91精品色在线| 国产亚洲午夜精品一区二区久久 | 国产精品国产三级国产av玫瑰| 亚洲国产精品专区欧美| 久久精品人妻少妇| 国产女主播在线喷水免费视频网站 | 国产成人福利小说| 99在线人妻在线中文字幕| 亚洲18禁久久av| 最近视频中文字幕2019在线8| 九色成人免费人妻av| 麻豆国产97在线/欧美| 欧美性猛交黑人性爽| 乱人视频在线观看| 国产精品女同一区二区软件| 九色成人免费人妻av| 我要看日韩黄色一级片| 久久久久久久午夜电影| 亚洲av电影在线观看一区二区三区 | 一区二区三区四区激情视频| 久久久精品94久久精品| 国产精品日韩av在线免费观看| a级毛片免费高清观看在线播放| 欧美一区二区亚洲| 亚洲高清免费不卡视频| 国产精品久久久久久av不卡| 欧美xxxx黑人xx丫x性爽| 色播亚洲综合网| 91在线精品国自产拍蜜月| av视频在线观看入口| 日日干狠狠操夜夜爽| 天天躁日日操中文字幕| 2021少妇久久久久久久久久久| 午夜福利在线观看吧| 男插女下体视频免费在线播放| 69av精品久久久久久| 天天躁日日操中文字幕| 国内精品宾馆在线| 99热这里只有是精品在线观看| 内地一区二区视频在线| 国产亚洲av片在线观看秒播厂 | 狂野欧美白嫩少妇大欣赏| 久久人人爽人人片av| 热99re8久久精品国产| 成人亚洲精品av一区二区| 精品久久久久久成人av| 国产 一区 欧美 日韩| 啦啦啦韩国在线观看视频| 啦啦啦韩国在线观看视频| 亚洲无线观看免费| 黄色一级大片看看| 亚洲成色77777| 精品一区二区三区人妻视频| 波多野结衣高清无吗| 最新中文字幕久久久久| 亚洲精品色激情综合| 最近最新中文字幕免费大全7| 国产av不卡久久| 国产一区亚洲一区在线观看| 国产精品爽爽va在线观看网站| 狠狠狠狠99中文字幕| 禁无遮挡网站| 国产淫片久久久久久久久| 99热全是精品| 国产在线一区二区三区精 | 91午夜精品亚洲一区二区三区| 午夜日本视频在线| 26uuu在线亚洲综合色| 精品久久久久久电影网 | 天堂影院成人在线观看| 欧美潮喷喷水| 尾随美女入室| 国产精品国产高清国产av| 亚洲av日韩在线播放| 91精品伊人久久大香线蕉| 久久久久久伊人网av| 七月丁香在线播放| 丰满乱子伦码专区| 久久精品久久精品一区二区三区| 毛片女人毛片| 午夜激情欧美在线| 国产精品伦人一区二区| 午夜精品国产一区二区电影 | 纵有疾风起免费观看全集完整版 | 不卡视频在线观看欧美| av福利片在线观看| 欧美一区二区精品小视频在线| 男人舔女人下体高潮全视频| 97超视频在线观看视频| 欧美性猛交黑人性爽| 欧美性感艳星| 男插女下体视频免费在线播放| 少妇熟女欧美另类| 尤物成人国产欧美一区二区三区| 亚洲精品乱码久久久v下载方式| 国产亚洲一区二区精品| 啦啦啦啦在线视频资源| 国产乱人偷精品视频| 蜜桃久久精品国产亚洲av| 天堂网av新在线| 中文字幕av在线有码专区| 亚洲性久久影院| 精品不卡国产一区二区三区| 成人毛片a级毛片在线播放| 亚洲精品成人久久久久久| 国产亚洲91精品色在线| 国产黄色小视频在线观看| 精品一区二区三区视频在线| 国产亚洲一区二区精品| 老女人水多毛片| 成人一区二区视频在线观看| 亚洲性久久影院| 日本黄色视频三级网站网址| 国产av在哪里看| 美女黄网站色视频| 久久99热这里只频精品6学生 | av线在线观看网站| 免费av观看视频| 人人妻人人澡人人爽人人夜夜 | 内射极品少妇av片p| 你懂的网址亚洲精品在线观看 | 午夜激情欧美在线| 日本三级黄在线观看| 麻豆成人av视频| 插逼视频在线观看| 久久久亚洲精品成人影院| 亚洲精品成人久久久久久| 国产男人的电影天堂91| 女人十人毛片免费观看3o分钟| 日韩 亚洲 欧美在线| 国产高清三级在线| 亚洲成av人片在线播放无| 夜夜看夜夜爽夜夜摸| 午夜福利在线在线| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美日韩卡通动漫| 1024手机看黄色片| 欧美高清性xxxxhd video| 看黄色毛片网站| 色网站视频免费| 久久热精品热| 菩萨蛮人人尽说江南好唐韦庄 | 秋霞在线观看毛片| 99久久人妻综合| 91精品一卡2卡3卡4卡| 亚洲国产高清在线一区二区三| 午夜福利在线在线| www.色视频.com| 九色成人免费人妻av| 六月丁香七月| 三级毛片av免费| 国产乱人偷精品视频| 高清视频免费观看一区二区 | 国产精品一区二区在线观看99 | 九草在线视频观看| 中文亚洲av片在线观看爽| 欧美日本亚洲视频在线播放| 精品久久久久久久久亚洲| 嫩草影院精品99| 99久久精品国产国产毛片| 久热久热在线精品观看| 国产伦精品一区二区三区四那| 女人久久www免费人成看片 | 少妇猛男粗大的猛烈进出视频 | 日本色播在线视频| 精品午夜福利在线看| 丰满少妇做爰视频| 免费黄网站久久成人精品| 少妇的逼水好多| 欧美激情久久久久久爽电影| 亚洲激情五月婷婷啪啪| 亚州av有码| 亚洲精品影视一区二区三区av| 精品无人区乱码1区二区| 国产亚洲精品久久久com| 日本欧美国产在线视频| 国产av一区在线观看免费| 麻豆成人av视频| 国产爱豆传媒在线观看| 97超视频在线观看视频| 一二三四中文在线观看免费高清| 久久精品久久久久久久性| 亚洲丝袜综合中文字幕| 天堂中文最新版在线下载 | 国产精品久久视频播放| 极品教师在线视频| 日韩在线高清观看一区二区三区| 边亲边吃奶的免费视频| 欧美精品国产亚洲| 国语对白做爰xxxⅹ性视频网站| 国产91av在线免费观看| 亚洲国产精品合色在线| 国产av在哪里看| 麻豆国产97在线/欧美| 午夜免费激情av| 国产在视频线在精品| 久久久久久久久久成人| 两个人的视频大全免费| 亚洲国产精品专区欧美| 嫩草影院新地址| 秋霞在线观看毛片| 99热精品在线国产| 国产老妇女一区| 国产一区二区在线观看日韩| 99九九线精品视频在线观看视频| 免费黄色在线免费观看| 啦啦啦韩国在线观看视频| 人体艺术视频欧美日本| 国产一区二区在线观看日韩| av在线播放精品| 嘟嘟电影网在线观看| 国产精品人妻久久久久久| 成人毛片a级毛片在线播放| 日本猛色少妇xxxxx猛交久久| 午夜福利视频1000在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品一区二区在线观看99 | 国产一级毛片七仙女欲春2| 丰满人妻一区二区三区视频av| 精品久久久久久久末码| 99久久人妻综合| 少妇熟女aⅴ在线视频| 久久久午夜欧美精品| 97在线视频观看| 最近的中文字幕免费完整| 亚洲精华国产精华液的使用体验| 成人毛片a级毛片在线播放| 亚洲内射少妇av| 久久久国产成人免费| 丝袜喷水一区| 白带黄色成豆腐渣| 搡女人真爽免费视频火全软件| 男女视频在线观看网站免费| 欧美日韩精品成人综合77777| 天美传媒精品一区二区| 三级毛片av免费| 久久久久久久久中文| 成人毛片a级毛片在线播放| 国产成人福利小说| 国产老妇伦熟女老妇高清| 亚洲人成网站高清观看| 亚洲国产高清在线一区二区三| 午夜a级毛片| av卡一久久| 精品久久久久久成人av| 伦理电影大哥的女人| 色噜噜av男人的天堂激情| 国产黄色小视频在线观看| 欧美成人精品欧美一级黄| www.av在线官网国产| 亚洲欧美成人综合另类久久久 | 波多野结衣巨乳人妻| 国产爱豆传媒在线观看| 国产黄a三级三级三级人| 欧美色视频一区免费| 精品久久久噜噜| 小说图片视频综合网站| 国产成人免费观看mmmm| 免费观看的影片在线观看| 十八禁国产超污无遮挡网站| 一区二区三区高清视频在线| 精品久久久久久久久av| 少妇人妻精品综合一区二区| 欧美高清成人免费视频www| 亚洲激情五月婷婷啪啪| 国产亚洲最大av| 三级国产精品片| 亚洲国产精品合色在线| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美精品免费久久| 国产av一区在线观看免费| 插逼视频在线观看| 国产成人福利小说| 亚洲人与动物交配视频| 欧美日韩国产亚洲二区| 成人三级黄色视频| 亚洲精品影视一区二区三区av| 十八禁国产超污无遮挡网站| 天堂中文最新版在线下载 | 久久精品国产亚洲网站| 啦啦啦韩国在线观看视频| 汤姆久久久久久久影院中文字幕 | 亚洲av成人精品一二三区| 国产精品久久视频播放| 男女那种视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 九色成人免费人妻av| 国产成人a∨麻豆精品| 91在线精品国自产拍蜜月| 九九久久精品国产亚洲av麻豆| 国产精品国产三级国产专区5o | 春色校园在线视频观看| 国产精品一区二区三区四区久久| 长腿黑丝高跟| 能在线免费观看的黄片| 国产伦在线观看视频一区| 午夜日本视频在线| 免费观看性生交大片5| 精品久久久久久久久久久久久| 国产高清不卡午夜福利| 天堂√8在线中文| 最近2019中文字幕mv第一页| 亚洲无线观看免费| 日日撸夜夜添| 国产黄色视频一区二区在线观看 | 噜噜噜噜噜久久久久久91| 亚洲av电影不卡..在线观看| 欧美最新免费一区二区三区| 久久精品影院6| 99热6这里只有精品| 97人妻精品一区二区三区麻豆| 精品一区二区免费观看| 狂野欧美白嫩少妇大欣赏| 国产探花极品一区二区| 亚洲国产精品sss在线观看| 嫩草影院新地址| 国产亚洲av片在线观看秒播厂 | 看片在线看免费视频| 亚洲成人精品中文字幕电影| 亚洲高清免费不卡视频| 美女被艹到高潮喷水动态| 国产午夜精品久久久久久一区二区三区| av卡一久久| 欧美潮喷喷水| 欧美变态另类bdsm刘玥| 麻豆成人av视频| 亚洲av男天堂| 三级毛片av免费| 亚洲,欧美,日韩| 亚洲综合色惰| 国产麻豆成人av免费视频| 免费无遮挡裸体视频| 婷婷色综合大香蕉| 91精品一卡2卡3卡4卡| 亚洲18禁久久av| 国产69精品久久久久777片| 国产av码专区亚洲av| 亚洲一级一片aⅴ在线观看| 亚洲国产成人一精品久久久| 22中文网久久字幕| 国产精华一区二区三区| 黄色配什么色好看| 淫秽高清视频在线观看| 国产单亲对白刺激| 国产高潮美女av| av在线天堂中文字幕| 秋霞伦理黄片| 国产成人免费观看mmmm| 久久99热6这里只有精品| 国产女主播在线喷水免费视频网站 | 联通29元200g的流量卡| 日本五十路高清| 亚洲欧美日韩卡通动漫| 成人欧美大片| 久久久久久久久大av| 三级国产精品片| 亚洲欧美日韩卡通动漫| 能在线免费观看的黄片| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 欧美3d第一页| 纵有疾风起免费观看全集完整版 | 日韩人妻高清精品专区| 久久99热这里只频精品6学生 | 老司机影院毛片| 九九热线精品视视频播放| 国产色爽女视频免费观看| 自拍偷自拍亚洲精品老妇| 岛国在线免费视频观看| 亚洲人成网站高清观看| 91aial.com中文字幕在线观看| 2021少妇久久久久久久久久久| 国产精品不卡视频一区二区| 久久婷婷人人爽人人干人人爱| 熟女人妻精品中文字幕| 黄色配什么色好看| 色哟哟·www| 级片在线观看| 日产精品乱码卡一卡2卡三| 国产精品久久久久久av不卡| 少妇的逼水好多| 国产成人91sexporn| 少妇熟女欧美另类| 国产熟女欧美一区二区| 中国国产av一级| 亚洲三级黄色毛片| 在线观看一区二区三区| a级毛色黄片| 身体一侧抽搐| 国产黄片视频在线免费观看| 国产伦精品一区二区三区四那| 亚洲国产精品专区欧美| 婷婷色av中文字幕| 身体一侧抽搐| 嫩草影院入口| 91狼人影院| 日本猛色少妇xxxxx猛交久久| 亚洲精品一区蜜桃| 深夜a级毛片| 2021少妇久久久久久久久久久| 亚洲av二区三区四区| 婷婷六月久久综合丁香| 亚洲最大成人手机在线| 亚洲最大成人中文| 久久国产乱子免费精品| 亚洲婷婷狠狠爱综合网| 男的添女的下面高潮视频| 噜噜噜噜噜久久久久久91| 国产精品国产三级国产专区5o | 晚上一个人看的免费电影| 高清午夜精品一区二区三区| 99九九线精品视频在线观看视频| av在线老鸭窝| 精华霜和精华液先用哪个| 亚洲人成网站高清观看| 国产精品电影一区二区三区| 国产伦精品一区二区三区视频9| 午夜激情欧美在线| av.在线天堂| 一区二区三区四区激情视频| 天堂av国产一区二区熟女人妻| 波多野结衣高清无吗| 中文乱码字字幕精品一区二区三区 | 欧美激情在线99| 亚洲精华国产精华液的使用体验| 长腿黑丝高跟| 亚洲精品乱码久久久v下载方式| 亚洲精品,欧美精品| 1024手机看黄色片| 你懂的网址亚洲精品在线观看 | 欧美不卡视频在线免费观看| 一边亲一边摸免费视频| 久久久a久久爽久久v久久| 看非洲黑人一级黄片| 韩国高清视频一区二区三区| ponron亚洲| 99久久精品国产国产毛片| 看十八女毛片水多多多| 亚洲av电影在线观看一区二区三区 | av播播在线观看一区| 能在线免费观看的黄片| 精品熟女少妇av免费看| 卡戴珊不雅视频在线播放| 国内精品一区二区在线观看| 欧美性感艳星| 波多野结衣巨乳人妻| 少妇的逼好多水| 午夜精品一区二区三区免费看| 日韩精品有码人妻一区| 国产亚洲一区二区精品| 国产乱人视频| 老司机影院成人| 国产一级毛片七仙女欲春2| 精品久久久久久久末码| 国产精品国产高清国产av| 国产黄色小视频在线观看| 日韩成人伦理影院| 丰满乱子伦码专区| 国产淫语在线视频| 国产男人的电影天堂91| 十八禁国产超污无遮挡网站| 日韩精品有码人妻一区| 国产一级毛片在线| 国产三级在线视频| 国产片特级美女逼逼视频| 久久精品91蜜桃| 午夜老司机福利剧场| 看十八女毛片水多多多| 岛国在线免费视频观看| 日韩国内少妇激情av| 成人午夜高清在线视频| 国产精品99久久久久久久久| 舔av片在线| 波多野结衣巨乳人妻| 亚洲精品色激情综合| 成人亚洲精品av一区二区| 少妇人妻一区二区三区视频| 韩国高清视频一区二区三区| 又粗又爽又猛毛片免费看| 国产精品嫩草影院av在线观看| 九九在线视频观看精品| 国语对白做爰xxxⅹ性视频网站| 老司机影院毛片| 成人无遮挡网站| 国产高清不卡午夜福利| 久久午夜福利片| 熟女电影av网| 欧美成人一区二区免费高清观看| 网址你懂的国产日韩在线| 亚洲av熟女| 精品不卡国产一区二区三区| 99久久无色码亚洲精品果冻| 黄色一级大片看看| 丝袜美腿在线中文| 免费无遮挡裸体视频| 性插视频无遮挡在线免费观看| 91狼人影院| 一级二级三级毛片免费看| 亚洲在久久综合| 一夜夜www| 亚洲在久久综合| 卡戴珊不雅视频在线播放| 国产精品1区2区在线观看.| 97热精品久久久久久| 久久久久国产网址| 91在线精品国自产拍蜜月| 国产成人freesex在线| 精品免费久久久久久久清纯| 国产探花极品一区二区| 岛国在线免费视频观看| 国产毛片a区久久久久| 久久韩国三级中文字幕| 99久久九九国产精品国产免费| 我要搜黄色片| 日本熟妇午夜| 丝袜美腿在线中文| 黄色日韩在线| 春色校园在线视频观看| 边亲边吃奶的免费视频| 久久人人爽人人片av| 欧美日本亚洲视频在线播放| 久久草成人影院| 久久久久精品久久久久真实原创| 国产乱人偷精品视频| 国产成年人精品一区二区| 成人午夜精彩视频在线观看| 一级二级三级毛片免费看| 日韩国内少妇激情av| 毛片一级片免费看久久久久| .国产精品久久| 国内精品宾馆在线| 精品人妻偷拍中文字幕| 国产乱来视频区| av国产免费在线观看| 欧美成人午夜免费资源| www.色视频.com| 99热这里只有是精品50| 一边亲一边摸免费视频| 国产午夜精品一二区理论片| 免费大片18禁| 国产精品一区二区三区四区免费观看| 亚洲婷婷狠狠爱综合网| 亚洲色图av天堂| 亚洲综合精品二区| 91久久精品电影网| 秋霞在线观看毛片| 亚洲国产精品sss在线观看| a级毛色黄片| eeuss影院久久| 国产乱人偷精品视频| 国产国拍精品亚洲av在线观看| 日韩欧美在线乱码| 观看美女的网站| 国产黄片美女视频| 一本一本综合久久| 日韩av不卡免费在线播放| 人妻系列 视频| 男人舔女人下体高潮全视频| videos熟女内射| 99在线人妻在线中文字幕| 亚洲欧美成人精品一区二区| 91在线精品国自产拍蜜月| 欧美区成人在线视频| 熟女人妻精品中文字幕| 国产精品无大码| 尤物成人国产欧美一区二区三区| av又黄又爽大尺度在线免费看 | 日韩av不卡免费在线播放| 又爽又黄无遮挡网站| 日本wwww免费看| 欧美一级a爱片免费观看看| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添av毛片| 伊人久久精品亚洲午夜| 乱人视频在线观看| 免费无遮挡裸体视频| 99久国产av精品国产电影| 日韩一本色道免费dvd| 国产一级毛片在线| 亚洲综合色惰| 99久久精品一区二区三区| 精品免费久久久久久久清纯| 精品不卡国产一区二区三区| 在线免费十八禁|