楊瑩 高珊 潘寶平 鄭津輝 高虹

摘要 利用轉錄組測序獲得了青蛤（Cyclina sinensis）白介素-1受體相關激酶-4（Interleukin-1 receptor-associated kinase 4，IRAK-4）基因的序列信息，利用巢式及熒光定量PCR克隆并分析了CsIRAK-4基因在青蛤各組織中的表達情況，進一步在鰻弧菌的脅迫下檢測了CsIRAK-4基因在青蛤血淋巴中的時序表達情況。結果顯示，CsIRAK-4基因的序列全長共2 687 bp，其開放閱讀框長1 926 bp，共編碼641個氨基酸，CsIRAK-4基因在青蛤的血淋巴、肝臟、外套膜、閉殼肌、性腺和鰓中均表達，其中在血淋巴中表達量最高，肝臟中表達量最低。青蛤在鰻弧菌脅迫下該基因的表達量在3 h即達到最大值，與對照組差異極顯著（P<0.01），證明該基因所表達的產物是青蛤重要的免疫信號通路蛋白。

關鍵詞 青蛤；轉錄組文庫；熒光定量PCR；CsIRAK-4

中圖分類號 S917 文獻標識碼 A 文章編號 0517-6611（2015）27-033-04

青蛤（Cyclina sinensis）是我國重要的海產經濟貝類，棲息于富于泥沙、泥漿的沿海潮汐地帶，廣泛分布于東海、黃海、南海等地[1]。青蛤的營養(yǎng)價值豐富，肉質鮮美，且其生長周期短，耐鹽性好，存活率高，在離水后可較長時間生存[2-3]，對于不良生存環(huán)境有較好的抵御能力[4-5]。近年來，青蛤增養(yǎng)殖產業(yè)的擴大及養(yǎng)殖環(huán)境的惡化，導致青蛤的病害連年發(fā)生，使該項產業(yè)受到嚴重威脅[6-8]。因此對于青蛤免疫抗病害研究具有重要的實踐意義。

貝類屬于非特異性免疫系統(tǒng)，其中Toll樣受體（Toll-like receptors，TLRs）是非特異性免疫反應中一類重要的識別受體，在對抗病原物的侵染過程中起著關鍵作用[9-14]。白介素-1受體相關激酶（Interleukin-1 receptor-associated kinases，IRAKs）是TLRs信號通路中的重要成員之一[9，15]。其中白介素-1受體相關激酶-4（Interleukin-1 receptor-associated kinase 4，IRAK-4）是IRAKs家族中最后發(fā)現(xiàn)的一個成員[16-17]。IRAK-4基因與IRAKs家族中其他基因結構相似，但缺少一個C端結構域，且行使功能時需要利用其自身激酶的活性[18-19]。研究表明，IRAK-4基因與細胞免疫過程中T細胞受體的識別過程密切相關[20]，這表明，先天免疫系統(tǒng)與適應性免疫系統(tǒng)之間有著錯綜復雜的聯(lián)系，而IRAK-4基因在兩者之間起著媒介作用。由此可見，IRAK-4基因在機體的免疫過程中至關重要。

目前，IRAK-4基因在脊椎動物中的相關研究較多，而在無脊椎動物中研究極少，僅在長牡蠣（Crassostrea gigas）和雜色鮑（Haliotis diversicolor）中有相關報道。在與青蛤同源性最高的雜色鮑中，IRAK-4基因在各個組織中均表達[9]。目前并未發(fā)現(xiàn)IRAK-4基因在青蛤中的相關研究，該研究在提高青蛤免疫能力及探索其抗病害機理方面具有重要的理論意義。

1 材料與方法

1.1 材料

青蛤樣品采集于天津大港灘涂，暫養(yǎng)于人工通氣海水中，海水密度1.02～1.04 g/cm3，水溫21～24 ℃，pH 7.0，投喂5‰小球藻。飼養(yǎng)7 d后，選取表面無破損，形態(tài)上無顯著差異的個體[平均殼寬（19.12±0.57）mm，平均殼長（29.14±1.23）mm，平均殼高（29.52±1.47）mm）]進行試驗。

1.2 方法

1.2.1 材料處理。

將鰻弧菌（Vibrio anguillarum）菌種在2216E培養(yǎng)基中培養(yǎng)24 h，后用無菌海水離心洗脫3次，重懸菌液，將濃度調為OD600=0.4。隨機分組，每次試驗設置10個平行組。將鰻弧菌菌液注入至試驗組青蛤體內，50 μl/只，對照組注射等量的滅菌海水。在注射前提取血淋巴、肝臟、外套膜、鰓、閉殼肌和性腺，分別在注射后3、6、12、24、48、96 h提取血淋巴。準確稱取各組織50 mg，后放入液氮中冷凍備用。

1.2.2 青蛤轉錄組文庫的構建。

利用TRIZOL提取青蛤各個組織的總RNA，用QIAGEN公司的Oligotex mRNA Kits試劑盒分離mRNA。采用第二代MiSeq測序儀，pair end 雙端模式完成青蛤轉錄組測序，利用De novo RNA-seq analysis技術，綜合分析相關基因功能注釋和代謝途徑，后通過篩選得到CsIRAK-4基因的類似序列。

1.2.3 生物信息學分析。

利用獲得的CsIRAK-4基因類似序列設計克隆引物CsIRAK-4-S：5′-AGATGGAGGCAGGGTGAT-3′，CsIRAK-4-A：5′-TGGCAGAGGTGGCGTATT-3′，并進行克隆，與GenBank中的核酸數據庫進行BLASTX比對，利用開放閱讀框（open reading frame，ORF）在線分析軟件分析其開放閱讀框，使用Sig-nalP 3.0分析該基因的信號肽序列，利用SMART查找其結構域，ProtParam工具在線預測此序列的分子量、分子式和等電點，使用Clustal W對氨基酸序列進行同源性分析和多重比對。

1.2.4 IRAK-4基因在青蛤各組織內的表達。

利用TRIZOL法提取青蛤血淋巴、肝臟、外套膜、閉殼肌、鰓和性腺的總RNA，反轉成cDNA，-20 ℃保存?zhèn)溆谩Ｒ驭?actin基因為內參基因，實時定量引物分別為βactin-S：5′-CACCACAACTGCCGAGAG-3′，βactin-A：5′-C-CGATAGTGATGACCTGACC-3′；CsIRAK-4-S：5′-AGATGGAGGCAGGGTGAT-3′，CsIRAK-4-A：5′-TGGCAGAGGTGGCGTATT-3′；反應在 Rotor-Gene 6000實時定量 PCR 儀上進行，擴增體系為 20 μl，反應程序：95 ℃預變性30 s，94 ℃變性5 s，56 ℃退火30 s，72 ℃ 延伸 30 s，40 個循環(huán)。采用2-ΔΔCT法進行數據處理，采用SPSS軟件進行數據分析。

1.2.5 在鰻弧菌刺激下青蛤IRAK-4基因在血淋巴內的時序性表達。

利用TRIZOL法提取鰻弧菌侵染后各時間點血淋巴的總RNA，反轉成cDNA。熒光定量PCR引物、反應體系及程序等參照“1.2.4”。采用2-ΔΔCT法進行數據處理，采用SPSS軟件進行數據分析。

2 結果與分析

2.1 青蛤IRAK-4基因的結構

CsIRAK-4基因開放閱讀框長1 926 bp，共編碼641個氨基酸（圖1）；其共有2個結構域：DEATH結構域（35～129）和S_Tkc結構域（380～638）；其理論分子量為72 552.8，等電點為5.89，分子式為C3206H5069N879O1005S17，氨基酸組成中亮氨酸（Leu）最高，占9.2%。此基因在GenBank的注冊號為KR732936。

2.2 青蛤IRAK-4基因的分子系統(tǒng)學分析

將CsIRAK-4基因與NCBI上已有的其他物種的IRAK-4基因進行同源性比對，使用MEGA4.1軟件構建CsIRAK-4基因的系統(tǒng)樹（圖2），利用bootstrap1000個循環(huán)檢驗拓撲結構的置信度。結果發(fā)現(xiàn)，雜色鮑（Haliotis diversicolor）、紅鮑螺（Haliotis rufescens）、霸王蓮花青螺（Lottia gigantea）、夏威夷短尾魷魚（Euprymna scolopes）與CsIRAK-4基因在同一分支，此基因與紅鮑螺（Haliotis rufescens）IRAK-4基因的同源性最高，達57%。

2.3 青蛤IRAK-4基因在不同組織間的表達

以β-actin基因在各組織中的表達量為內參對照，利用實時熒光定量PCR分析CsIRAK-4基因在其血淋巴、肝臟、外套膜、閉殼肌、性腺和鰓6個組織中的表達情況。結果表明，CsIRAK-4基因在此6個組織中普遍表達（圖3），但表達量明顯不同，在血淋巴中表達量最高，與其他組織相比差異極顯著（P<0.01），閉殼肌次之，肝臟中表達量最低。表明CsIRAK-4基因主要在血淋巴中表達。

2.4 在鰻弧菌刺激下青蛤IRAK-4基因在血淋巴中的時序性表達

以鰻弧菌侵染青蛤β-actin基因為內參基因，使用實時熒光定量PCR儀對CsIRAK-4基因在血淋巴中表達的時序性變化進行分析（圖4）。結果顯示，試驗組在刺激后3 h的表達量達到最大值，明顯高于對照組，約為對照組的17倍，并與對照組存在極顯著差異（P<0.01）；3 h后表達量開始下降。

3 討論

該研究通過構建青蛤轉錄組文庫的方法獲得了CsIRAK-4基因的類似序列，利用此序列設計引物并克隆，后進行測序，通過序列比對發(fā)現(xiàn)，該基因在進化中比較保守，與許多水生動物的IRAK-4基因具有明顯的相似性。通過SMART分析其蛋白序列，發(fā)現(xiàn)CsIRAK-4基因共有2個結構域，分別為死亡結構域（Death Domain）和激酶結構域（S_TKc Domain），死亡結構域位于第35～129個氨基酸之間，激酶結構域位于第380～638個氨基酸之間。死亡結構域作為一種中介將死亡受體和胞漿信號相連通，在控制機體的凋亡過程中有重要作用。在IRAK基因家族中，只有IRAK-1基因、IRAK-4基因具有包含激酶結構域的功能催化位點[21]，但只有IRAK-4基因在下游信號轉導通路中使用其激酶活性。同時，IRAK-4基因還具有一個高度保守的賴氨酸殘基，此殘基位于ATP結合位點，在ATP的水解過程中起關鍵作用[17，22]。IRAKs在TLRs的信號轉導過程中起著重要的媒介作用[20，23]，因此在先天免疫過程中IRAKs極為關鍵。

該研究結果顯示，CsIRAK-4基因在各個組織普遍表達，在血淋巴中表達量最高，閉殼肌次之，后依次為鰓、外套膜、性腺，在肝臟中表達量最低。貝類所具有的非特異性免疫體系中，血液的吞噬作用往往是抗病害的第一道防線，青蛤血液中含有大量的吞噬細胞、抗菌性物質及酶原物質，在病原物侵染過程中其免疫作用非常關鍵。血細胞中的CsIRAK-4基因上調最為明顯，說明血淋巴作為軟體動物的非特異性免疫防御的首要組織擔當著機體非特異性免疫防御的重任。而鰓和外套膜在軟體動物生存過程中最先接觸外界環(huán)境，是軟體動物抵御外界傷害的第一道防線，因此在其非特異性免疫過程中也極其重要。

CsIRAK-4所表達的產物是TLRs信號通路中重要的信號蛋白之一，當青蛤受到鰻弧菌脅迫后，其血淋巴中該基因在3 h即表達至最大值，說明病原物對TLRs免疫受體及下游 CsIRAK-4信號通路蛋白具有強烈的識別與誘導作用，可使其血液中快速合成相應的免疫信號蛋白，對抗外來脅迫物對機體自身的侵害。并且CsIRAK-4基因的表達量呈先升高后降低的趨勢，證明CsIRAK-4基因參與了非特異免疫的全部過程。該試驗結果為今后貝類增養(yǎng)殖中的病害防治提供了重要的試驗數據。

參考文獻

[1] YOON J M.Geographic variations and genetic distance of three geographic cyclina clam（Cyclina sinensis Gmelin）populations from the Yellow Sea[J].Development and reproduction，2012，16（4）：315-320.

[2] 王慧珍.浙江動物志：軟體動物[M].杭州：浙江科學技術出版社，1991：243-256.

[3] LIU W，MA Y，HU S，et al.Rearing Venus clam seeds Cycline sinensis on a commercial scale[J]. Acquaculture，2002，211（1）：109-114.

[4] WANG H Z.Zhejiang fauna[M].Molluscs articles，1991.

[5] WANG R C，WANG Z P，ZHANG J Z.Shellfish aquaculture science[M].Qingdao：Qingdao Ocean University Press，1993：208-209.

[6] 任毅鵬，高晶，潘寶平，等.青蛤（Cyclina sinensis）TLR2 基因的克隆與表達分析[J].海洋與湖沼，2014，45（5）：1037-1043.

[7] 孫國銘，萬夕和，劉培庭，等.通州海區(qū)灘涂青蛤死亡原因的初步分析[J].水產養(yǎng)殖，2004，25（2）：26-27.

[8] 曹華.沿海灘涂青蛤死亡原因初探及對策[J].科學養(yǎng)魚，2004（4）：47-48.

[9] GE H，WANG G D，ZHANG L L，et al.Molecular cloning and expression of interleukin-1 receptor-associated kinase 4，an important mediator of Toll-like receptor signal pathway，from small abalone Haliotis diversicolor[J].Fish & shellfish immunology，2011，30（4/5）：1138-1146.

[10] MEDZHITOV R，PRESTON-HURLBURT P，JANEWAY JR C A.A human homologue of the Drosophila Toll protein signals activation of adaptive immunity[J].Nature，1997，388（6640）：394-397.

[11] ZHANG D，ZHANG G，HAYDEN M S，et al.A toll-like receptor that prevents infection by uropathogenic bacteria[J].Science，2004，303（5663）：1522-1526.

[12] TAKEUCHI O，KAWAI T，SANJO H，et al.TLR6：A novel member of an expanding toll-like receptor family[J].Gene，1999，231（1/2）：59-65.

[13] CHUANG T H，ULEVITCH R J.Cloning and characterization of a sub-family of human toll-like receptors：hTLR7，hTLR8 and hTLR9[J].Eur Cytokine Netw，2000，11（3）：372-378.

[14] HEMMI H，TAKEUCHI O，KAWAI T，et al.A Toll-like receptor recognizes bacterial DNA[J].Nature，2000，408（6813）：740-745.

[15] JANSSENS S，BEYAERT R.Functional diversity and regulation of different inter-leukin-1 receptor-associated kinase（IRAK）family members[J].Mol Cell，2003，11：293-302.

[16] GOTTIPATI S，RAO N，F(xiàn)UNG-LEUNG W.IRAK1：A critical signaling mediator of innate immunity[J].Cell Signal，2008，20（2）：269-276

[17] LI S，STRELOW A，F(xiàn)ONTANA E，et al.IRAK-4：a novel member of the IRAK family with the properties of an IRAK-kinase[J].Proc Natl Acad Sci USA，2002，99（8）：5567-5572.

[18] EDWARDS D，TOWB P，WASSERMAN S.An activity-dependent network of inter-actions links the Rel protein Dorsal with its cytoplasmic regulators[J].Development，1997，124：38-55.

[19] TOWB P，BERGMANN A，WASSERMAN S.The protein kinase Pelle mediates feed-back regulation in the Drosophila Toll signaling pathway[J].Development，2001，128（23）：4729-4736.

[20] SUZUKI N，SUZUKI S，MILLAR D，et al.A critical role for the innate immune signaling molecule IRAK-4 in T cell activation[J].Science，2006，311（5769）：1927-1932.

[21] WIENS M，KORZHEV M，PEROVI -OTTSTADT S，et al.Toll-like receptors are part of the innate immune defense system of sponges（Demospongiae：Porifera）[J].Mol Biol Evol，2007，24（3）：792-804.

[22] ZHENG J，KNIGHTON D，TEN EYCK L，et al.Crystal structure of the catalytic subunit of cAMP-dependent protein kinase com- plexed with magnesium-ATP and peptide inhibitor[J]. Biochemistry，1993，32（9）：2154-2161.

[23] LASKER M，GAJJAR M，NAIR S.Cutting edge：Molecular structure of the IL-1R-associated kinase-4 death domain and its implications for TLR signaling[J].J Immunol，2005，175（7）：4175-4179.