急性高溫低氧脅迫下刺參體腔液細胞miRNA表達譜的構(gòu)建和靶標基因的功能分析?

尹英超，劉富祥，張玉欽，陳慕雁??

(1.中國海洋大學(xué)水產(chǎn)學(xué)院，山東 青島 266003；2.海水養(yǎng)殖教育部重點實驗室(中國海洋大學(xué))，山東 青島 266003；3.煙臺中集藍海洋科技有限公司，山東 煙臺 264000)

刺參(Apostichopusjaponicus)屬棘皮動物門(Echinodermata)海參綱(Holothuroidea)辛那參目(Synallactida)刺參科(Stichopodidae)，廣泛分布于中國、日本、韓國和俄羅斯等東北亞沿海[1]。刺參作為底棲生物以底泥、碎屑為食，在維持生態(tài)系統(tǒng)平衡中發(fā)揮著重要的作用[2]。此外，刺參具有很高的營養(yǎng)價值和藥用價值，長期以來在亞洲市場上作為一種食品和民間藥物備受歡迎，是中國重要的海水養(yǎng)殖經(jīng)濟品種之一，目前已經(jīng)成為北方水產(chǎn)養(yǎng)殖中的支柱性產(chǎn)業(yè)[1,3]。

全球氣候的變化導(dǎo)致海水溫度升高和缺氧現(xiàn)象已成為一個世界性環(huán)境問題，嚴重影響到海洋生物的生存?！岸驙柲嶂Z”是導(dǎo)致太平洋海水表面溫度異常升高的重要因素之一，可以引發(fā)嚴重而廣泛的天氣變化，如高溫和暴雨[4]。2013年夏末以來，中國山東、河北兩省因“厄爾尼諾”現(xiàn)象引發(fā)的連續(xù)高溫暴雨天氣，導(dǎo)致海水溫度高達30 ℃，溶解氧水平降至2 mg/L以下；由于刺參長期營底棲生活，活動緩慢，受環(huán)境因素影響非常顯著，極端高溫、低氧脅迫可導(dǎo)致刺參大量死亡，嚴重制約刺參養(yǎng)殖產(chǎn)業(yè)的可持續(xù)發(fā)展。為了給刺參應(yīng)對逆境關(guān)鍵指標篩選和抗逆品種選育提供理論支持，刺參在急性高溫、低氧脅迫下內(nèi)在分子調(diào)控機制的科學(xué)研究亟待開展。

棘皮動物通常缺乏特異性免疫，體腔液是棘皮動物進行機體防御的重要免疫組織。刺參體腔液中懸浮著各種體腔液細胞，同時也含有凝集素、水解酶、酚氧化酶、超氧化物歧化酶類白介素和類補體樣因子等體液免疫因子，兩者共同作用負責(zé)刺參的天然免疫防御[5-6]。已有大量研究表明棘皮動物在響應(yīng)外界溫度、鹽度、溶解氧和pH等環(huán)境脅迫因子時，體腔液及其細胞的相關(guān)免疫抗病能力會發(fā)生變化[7-10]。

miRNA是一類長度為21～24 nt的內(nèi)源性非編碼RNA，miRNA能夠通過種子序列與特定的靶基因位點相結(jié)合，引起mRNA降解或抑制mRNA翻譯，從而對靶基因的表達起到調(diào)節(jié)作用。目前，關(guān)于miRNA介導(dǎo)高溫、低氧脅迫的調(diào)節(jié)機制在水生生物中進行了廣泛研究。Zheng等在日本對蝦(Marsupenaeusjaponicus)體內(nèi)共發(fā)掘了15個miRNA參與高溫脅迫應(yīng)答，靶標基因的預(yù)測和功能分析表明這些miRNA參與了信號傳導(dǎo)、轉(zhuǎn)錄、抗脅迫應(yīng)激等細胞過程[11]。Zhou等研究發(fā)現(xiàn)虹鱒魚(Oncorhynchusmykiss)在致死溫度22.5 ℃脅迫48 h下，頭、腎組織中共發(fā)掘出12個上調(diào)和5個下調(diào)的miRNA，6個差異表達的miRNA共預(yù)測到22個靶基因，注釋結(jié)果表明在熱應(yīng)激下miRNA主要參與了虹鱒魚的代謝和免疫反應(yīng)[12]。Hadj-Moussa等發(fā)現(xiàn)miRNA調(diào)控的代謝抑制是保障洪堡魷魚(Dosidicusgigas)能夠在缺氧水域生存的重要機制，在缺氧誘導(dǎo)下魷魚的大腦、外套膜肌肉和心臟中表達量變化的miRNA主要參與了關(guān)鍵細胞功能維持和細胞保護的過程[13]。Wang等在南美白對蝦(Litopenaeusvannamei)應(yīng)對缺氧的研究中，鑒定出許多與缺氧耐受相關(guān)的miRNA，其差異表達的miRNA可能通過調(diào)節(jié)應(yīng)激反應(yīng)基因(如HIF1a和p53)的表達來調(diào)節(jié)宿主對缺氧的反應(yīng)，并影響與缺氧適應(yīng)相關(guān)的關(guān)鍵細胞事件[14]。Sun等在黒鱸(Micropterussalmoides)缺氧脅迫的研究中篩選出的13個差異表達的miRNAs的靶基因在血管內(nèi)皮生長因子(VEGF)信號通路、絲裂原活化蛋白激酶(MAPK)信號通路和磷脂酰肌醇信號系統(tǒng)中顯著富集[15]。關(guān)于刺參miRNA研究起步較晚，最早于2012年Li等從刺參的體腔液細胞中鑒定到miRNAs[16]，之后的研究表明miRNA廣泛參與了刺參補體系統(tǒng)、脂多糖以及Toll樣受體的調(diào)控[17-22]。近年來，關(guān)于miRNA在刺參響應(yīng)環(huán)境脅迫中的調(diào)控作用研究也相繼開展。Li等構(gòu)建了刺參在高溫脅迫下的miRNA表達譜，發(fā)現(xiàn)大量差異表達的miRNA參與免疫過程相關(guān)的調(diào)節(jié)[23]；Huo等在刺參缺氧脅迫研究中發(fā)現(xiàn)miRNA可能在氧化還原、運輸以及耐缺氧脅迫中發(fā)揮作用[24]。Tian等研究發(fā)現(xiàn)miR-10可能通過靶標TBC1D5基因參與低鹽脅迫下刺參的能量調(diào)節(jié)，并在維持細胞穩(wěn)態(tài)方面起到了關(guān)鍵作用[25]。Huo等以刺參呼吸樹為研究對象，構(gòu)建了在高溫、低氧脅迫下刺參的miRNA表達譜，并在高溫(26 ℃)、低氧(2 mg/L)和高溫低氧復(fù)合脅迫條件下，分別鑒定出21、26和22個差異表達miRNA，GO和KEGG通路分析顯示差異miRNA的潛在靶基因參與了生物合成、代謝、免疫、細胞生長和死亡、翻譯及信號轉(zhuǎn)導(dǎo)等過程[26]。

本研究以刺參高溫、低氧脅迫下的體腔液細胞為研究對象(由于高溫低氧協(xié)同組的刺參全部死亡，無法進行后續(xù)的研究)，采用Illumina測序技術(shù)，結(jié)合生物信息學(xué)分析，構(gòu)建了刺參在高溫、低氧脅迫下的miRNA表達譜，使用miRanda軟件，以構(gòu)建的三代全長轉(zhuǎn)錄組為參考文庫，預(yù)測了差異表達miRNA的靶標基因，并對候選靶標基因進行GO和KEGG富集分析，發(fā)掘了刺參在高溫、低氧脅迫下的差異miRNA參與調(diào)控的關(guān)鍵通路，探討了miRNA在刺參體腔液細胞響應(yīng)高溫、低氧脅迫中的潛在調(diào)控機制，為篩選刺參抗逆關(guān)鍵指標提供理論支持。

1 材料與方法

1.1 材料

1.1.1 實驗動物暫養(yǎng) 實驗所用的成熟刺參于2019年4月采集于威海刺參養(yǎng)殖場，低溫運到實驗室。挑選120頭健康無病、發(fā)育良好的刺參(體質(zhì)量(120±20)g)分別放入12個大小、形狀、材質(zhì)一致的藍色避光的養(yǎng)殖水箱中，每個水箱中放置10頭刺參，為了消除環(huán)境因素的差異，在海邊收集大小不一的碎石經(jīng)消毒清洗后，放入水箱中供刺參附著，箱內(nèi)水溫(T)為(15±1)℃，鹽度為30±1，pH為8.0±0.5，充氣保持溶解氧(DO)為8 mg/L，在此環(huán)境中馴化2周。期間每日上午定時清除水箱底部的糞便和剩余的飼料，飼料與底泥按1∶3比例混合投喂。

1.1.2 實驗處理 對照組的條件設(shè)置為T=15 ℃、DO=8 mg/L，高溫處理組T=30 ℃，DO=8 mg/L，低氧處理組T=15 ℃，DO=2 mg/L，高溫低氧處理組T=30 ℃，DO=2 mg/L。高溫設(shè)置方法為使用功率1 000 W的加熱棒加熱來達到實驗所需的水溫，并通過控溫器來保持水溫，水溫升高速率為2 ℃/h。低氧設(shè)置方法為向水體中充入氮氣，控制氮氣充入速率，使溶解氧下降速率保持在0.8 mg/L。達到實驗所需溶氧條件后，使用塑料薄膜覆蓋水箱的水面，在低氧脅迫期間使用便攜式溶氧儀實時監(jiān)測溶解氧水平。對照組和實驗組分別設(shè)置3個平行，對照組和實驗組在達到所需溫度溶氧條件后48 h進行取樣。

1.1.3 樣品收集 每個平行組中取4頭刺參，設(shè)3個平行組，每個處理組共取12頭刺參。將每個刺參的體腔液分別取3份，每份體積為1 mL，在4 ℃，1 000 r/min的冷凍離心機中離心30 s,去掉上清液后將剩余細胞混合，加入1 mL的Trizol試劑，立即放入液氮中速凍后轉(zhuǎn)入-80 ℃超低溫冰箱中保存?zhèn)溆?。由于刺參無法耐受48 h的高溫低氧(T=30 ℃, DO=2 mg/L)的協(xié)同作用而全部死亡，無法取得該實驗組樣品進行后續(xù)分析。

1.2 RNA文庫構(gòu)建

采用Trizol法提取細胞總RNA，具體步驟見李興可等[27]。使用Nanodrop檢測RNA的純度(OD260/280比值)；使用Qubit對RNA濃度進行精確定量；使用Agilent 2100精確檢測RNA的完整性。使用NEBNext? Multiplex Small RNA Library Prep Set for Illumina?(NEB, USA)生成測序文庫，文庫中對照組為CO1、CO2、CO3，高溫組為HS1、HS2、HS3，低氧組為HY1、HY2、HY3，建庫流程如圖1所示。

圖1 Illumina樣品檢測、建庫和測序的流程圖

1.3 miRNA測序與數(shù)據(jù)分析

采用Illumina Hiseq 2500進行miRNA測序，去除低質(zhì)量的reads(質(zhì)量值小于等于20 的堿基占整個reads的30%以上)；去除 N(N表示無法確定堿基信息)的比例大于10%的 reads；去除接頭污染的reads；去除3’接頭序列；去除polyA/T/G/C的reads。選擇一定長度的clean reads用于下游的分析，miRNA的長度集中在21～22 nt。將小RNA標簽映射到RepeatMasker和Rfam數(shù)據(jù)庫，除去源自蛋白質(zhì)編碼基因、重復(fù)序列、rRNA、tRNA、snRNA和snoRNA等。用bowtie將長度篩選后的sRNA定位到刺參基因組(NCBI BioProject ID: 29023486)序列上，分析small RNA在參考序列上的分布情況。

1.4 miRNA的鑒定和差異表達分析

使用miRBase20.0作為參考，鑒定已知的miRNA；結(jié)合miRNA前體發(fā)夾結(jié)構(gòu)特征，采用軟件miREvo和mirdeep2預(yù)測新的miRNA。將所有sRNA標簽與成熟miRNA進行比對，檢測可能有堿基編輯的miRNA，已知的miRNA使用miFam.dat(http://www.mirbase.org/ftp.shtml)查找家族。將新的miRNA前體提交到Rfam庫(http://rfam.sanger.ac.uk/search/)，以尋找Rfam家族。根據(jù)下面算法，應(yīng)用TPM(每百萬轉(zhuǎn)錄本)估算miRNA表達水平。

規(guī)范化表達式=(映射的讀數(shù)/總讀數(shù))×1 000 000。

使用DESeq R軟件包(1.8.3)對兩個實驗組進行差異表達分析。使用Benjamini&Hochberg方法調(diào)整P值。默認情況下，P<0.05為顯著差異表達。

1.5 靶標基因的預(yù)測和分析

利用miRanda分析軟件，以前期構(gòu)建的三代全長轉(zhuǎn)錄組數(shù)據(jù)作為參考文庫，從中提取3’非編碼區(qū)域，預(yù)測靶標基因。將篩選參數(shù)設(shè)置為自由能小于-83.68 kJ/mol和得分閾值大于140。

1.6 GO和KEGG富集分析

基于GOseq的Wallenius非中心超幾何分布調(diào)整基因長度偏差，用于GO富集分析。Pathway顯著性富集分析以KEGG Pathway為單位，應(yīng)用超幾何檢驗，該分析計算公式如下：

得到的P值經(jīng)Bonferroni方法校正，以校正后的P值小于0.05定義為顯著富集的pathway。

1.7 qRT-PCR驗證

應(yīng)用SYBR?Premix Ex TaqTMII Kit(Takara，Japan，Cat No.RR820)試劑盒，分別取刺參對照組、高溫組、低氧組的樣品各5個平行，每個平行樣品設(shè)置3個技術(shù)重復(fù)，將成熟miRNA(21～23 nt)的完整序列用作miRNA特異性5’端引物，引物序列如表1所示。

表1 miRNA qRT-PCR的引物序列

qPCR的3’端引物使用試劑盒提供的mRQ 3’通用引物，以U6作為內(nèi)參基因，其引物序列為：F:GGAACGATACAGAGAAGATTAGC, R: TGGAACG-CTTCACGAATTTGCG。變性的反應(yīng)條件：95 ℃ 10 s，95 ℃ 5 s，60 ℃ 20 s，40個循環(huán)。融解曲線的反應(yīng)條件：95 ℃ 60 s，55 ℃ 30 s，95 ℃ 30 s。使用2-ΔΔCt法計算 miRNA 的相對表達水平，使用單因素方差分析miRNA的組間差異顯著性，以P<0.05為差異顯著。

2 結(jié)果與分析

2.1 測序結(jié)果分析

建庫之后，利用Illumina HiSeqTM2500測序平臺將原始圖像數(shù)據(jù)文件經(jīng)過堿基識別轉(zhuǎn)化為測序序列，得到的Raw Data以FASTQ格式儲存，原始數(shù)據(jù)已上傳至NCBI(BioProject ID: PRJNA723369, SubmissionID: SUB9504899)。檢查測序錯誤率異常的堿基位置分布，發(fā)現(xiàn)9個文庫測序結(jié)果中每個位置堿基的錯誤率都低于0.5%，詳細結(jié)果見附表1。

2.2 sRNA的測序分析結(jié)果

為了保證信息分析的質(zhì)量，對得到的Raw Data數(shù)據(jù)進行過濾，得到Lean Data，篩選出sRNA進行后續(xù)的分析，詳細結(jié)果見附表2。

用bowtie將長度篩選后的sRNA比對到刺參參考基因組序列上，發(fā)現(xiàn)每組中90%以上sRNA能夠比對上參考基因組，其中HY3組能比對到參考基因組的數(shù)量最少(91.96%)，HY1組能比對上參考基因組的數(shù)量最多(98.27%)，與參考基因組方向相同鏈的reads數(shù)占全部clean reads數(shù)的83%以上，方向相反鏈的reads數(shù)占全部clean reads數(shù)的10%以下，結(jié)果如表2所示。

表2 sRNA與參考基因組比對結(jié)果統(tǒng)計

2.3 已知miRNA與ncRNA分析

所有樣本中共比對到33個已知的miRNA，比對上的sRNA的種類有3 242種。對Novel miRNA成熟體、前體、種類和總數(shù)進行統(tǒng)計，詳細結(jié)果見附表3。總共匹配到的前體和成熟體的數(shù)目分別為149和159，詳細結(jié)果見附表4。

2.4 miRNA家族分析

對檢測到的已知miRNA和新預(yù)測到的miRNA進行家族分析，結(jié)果顯示，已知miRNA和新miRNA總共屬于31個不同的miRNA家族，其中l(wèi)et-7在物種中的分布最豐富。大多數(shù)miRNA只檢測到一個前體，miR-10、miR-124和miR-25檢測到3個前體。miRNA家族分析部分結(jié)果見附表5。

2.5 miRNA差異表達分析

在HS組中共發(fā)掘出20個差異表達的miRNA，其中11個上調(diào)，9個下調(diào)；在HY組中共發(fā)掘出10個差異表達的miRNA，其中6個上調(diào)、4個下調(diào)(見表3)。兩組中共同差異表達的miRNA有6個：novel-23、novel-43、novel-45、novel-49、miR-200-3p和miR-71。為驗證結(jié)果的準確性，分別從高溫組和低氧組中隨機選取了10個差異miRNA，其中高溫組選取miR-133、miR-71、miR-10、miR-2012、miR-200-3p和Novel-32，低氧組選取miR-2006、miR-71、miR-200-3p和Novel-43，利用熒光定量PCR進行驗證。結(jié)果表明，高通量測序結(jié)果和qRT-PCR結(jié)果趨勢相同(見圖2)。

圖2 差異表達miRNAs的qRT-PCR驗證結(jié)果

表3 高溫組和低氧組差異表達miRNA的詳細結(jié)果

2.6 靶基因預(yù)測和功能富集分析

我們以刺參的三代全長轉(zhuǎn)錄組中的191 141個轉(zhuǎn)錄本作為參考文庫，預(yù)測差異表達miRNA的靶基因。結(jié)果表明，高溫組中20個差異表達miRNA共預(yù)測到14 827個靶基因，低氧組中10個差異表達的miRNA共預(yù)測到8 577個靶基因。

我們對每組的差異表達miRNA的靶基因分別進行了GO和KEGG富集分析。富集分析結(jié)果顯示，在高溫脅迫組中，差異表達miRNA的靶向基因在細胞組分中聚類最明顯的是細胞成分，在分子功能項中聚類最明顯的是轉(zhuǎn)運活性和跨膜轉(zhuǎn)運體活性(Transmembrane transporter activity)，在生物學(xué)進程中跨膜運輸(Transmembrane transport)聚類最為明顯。在低氧脅迫組中，差異表達miRNA的靶向基因在細胞組分中聚類最顯著的是細胞成分和膜部分(Membrane part)；在分子功能項中聚類最顯著的是轉(zhuǎn)運活性，在生物學(xué)進程中內(nèi)肽酶活性(Endopeptidase activity)聚類最顯著，詳細結(jié)果如圖3所示。

圖3 差異表達miRNA靶基因的GO富集分析

KEGG富集分析結(jié)果顯示，高溫組中顯著富集的通路有36個，其中富集到基因數(shù)目較多的通路主要集中在磷脂酰肌醇激酶-絲氨酸/蘇氨酸激酶信號通路(PI3K-Akt signaling pathway)、癌癥通路(Pathways in cancer)、內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工(Protein processing in endoplasmic reticulum)和黏著斑(Focal adhesion)。在低氧組中顯著富集的通路有37個，富集靶基因數(shù)目較多的通路主要集中在黏著斑(Focal adhesion)、癌癥通路、磷脂酰肌醇激酶-絲氨酸/蘇氨酸激酶信號通路和絲裂原活化蛋白激酶信號通路(MAPK signaling pathway-fly)，結(jié)果如圖4所示。

((a)高溫組；(b)低氧組。(a)Heat stress group；(b)Hypoxia group.)

3 討論

隨著高通量測序技術(shù)的快速發(fā)展，刺參不同組織中(腸道、肌肉、體壁、管足以及體腔液等)miRNA表達譜得到廣泛的構(gòu)建[22]，表達譜分析表明miRNA在刺參生物學(xué)過程中具有潛在的關(guān)鍵調(diào)控作用，但在刺參體腔液細胞中miRNA響應(yīng)高溫、低氧脅迫的調(diào)控機制的研究還少有報道。

體腔液作為刺參內(nèi)外環(huán)境的滲透介質(zhì)，外界環(huán)境的變化可能影響體腔液細胞的活性和功能，高溫和低氧是目前影響刺參生存最重要的環(huán)境因子。本研究共構(gòu)建了9個刺參體腔液細胞的miRNA文庫，測序結(jié)果顯示大部分成熟體miRNA長度分布在21～23 nt，峰值大小為22 nt，這與大部分海洋無脊椎動物的研究結(jié)果相似，略微的差異可能是由物種差異和組織特異性導(dǎo)致的。

本研究從刺參的體腔液細胞中共鑒定到33個已知的miRNA和149個Novel miRNA。多數(shù)miRNA在各物種中是高度保守的，發(fā)揮了廣泛而重要的調(diào)控作用。其中高溫組中miR-133、miR-10、miR-2012、miR-200-3p、miR-71和低氧組中miR-1、miR-2006、miR-200-3p、miR-71都有顯著差異表達。有研究表明，miR-10是刺參體腔液細胞中最豐富的miRNA[16]，且在細胞凋亡[28]、蛋白合成[29]、胚胎發(fā)育和分化[30-31]、造血干細胞調(diào)動[32]和炎癥[33]等生物過程中發(fā)揮著重要的調(diào)控作用，這表明miR-10在刺參響應(yīng)高溫脅迫中同樣可能參與調(diào)控多個生理學(xué)過程。De Lencastreetal等研究發(fā)現(xiàn)miR-71與細胞的衰老和DNA損傷反應(yīng)有關(guān)[34]，另有研究發(fā)現(xiàn)miR-71同樣可以增強對熱休克和氧化應(yīng)激的抗性[35]。Huo等最近研究發(fā)現(xiàn)miR-71家族成員在高溫脅迫中是顯著上調(diào)的[26]，這和本文的結(jié)果一致，也暗示了本研究中miR-71在刺參響應(yīng)高溫、低氧脅迫下的潛在功能。miR-1已被證實與凋亡相關(guān)基因如熱休克蛋白基因相關(guān)，并間接調(diào)節(jié)eNOs以及細胞凋亡[36]，以往研究也發(fā)現(xiàn)刺參呼吸樹中miR-1在低氧脅迫下顯著上調(diào)[24]，這和本文的研究結(jié)果一致。miR-133家族因為在抑制人類癌癥發(fā)生中的關(guān)鍵作用被充分研究[37]，在心肌細胞中miR-133通過靶向Hsp70抑制細胞凋亡[38]，相關(guān)研究還直接證明miR-133通過靶向AjIRAK-1參與Toll樣受體(TLR)級聯(lián)調(diào)節(jié)，從而促進刺參體腔細胞吞噬功能[39]。本研究中miR-133在高溫脅迫下顯著上調(diào)，這暗示了miR-133可能參與了高溫脅迫下體腔液細胞抗凋亡過程和免疫應(yīng)激。miR-200家族是哺乳動物中最著名的miRNA家族之一。作為該家族的重要一員，miR-200-3p被報道參與細胞增殖、凋亡和腫瘤發(fā)生等生物學(xué)過程[40-41]。本課題組在刺參夏眠的研究中也發(fā)現(xiàn)了miR-200-3p能夠通過靶向AjEHHADH和AjCA參與刺參夏眠期間氧化損傷和細胞周期調(diào)控[42-43]。在本研究中，miR-200-3p在高溫和低氧脅迫下均顯著上調(diào)，揭示了miR-200-3p可能在刺參響應(yīng)高溫、低氧脅迫下參與了廣泛的生物學(xué)功能調(diào)節(jié)。目前關(guān)于miR-2012的功能研究還有待開展。

為進一步發(fā)掘miRNA潛在的調(diào)控機制，我們以三代轉(zhuǎn)錄組數(shù)據(jù)為參考文庫，預(yù)測了差異表達miRNA的候選靶基因。單個miRNA可以靶向多個靶基因，每個靶基因也可以被多個miRNA調(diào)控，因此，miRNA通過介導(dǎo)復(fù)雜的網(wǎng)絡(luò)結(jié)構(gòu)發(fā)揮調(diào)控功能。靶基因GO分析結(jié)果顯示，大多數(shù)靶基因與細胞組分、轉(zhuǎn)運以及跨膜轉(zhuǎn)運活性有關(guān)。靶基因KEGG通路富集分析發(fā)現(xiàn)多數(shù)基因在絲裂原活化蛋白激酶信號通路、磷脂酰肌醇激酶-絲氨酸/蘇氨酸激酶信號通路、寡聚化核苷酸結(jié)合結(jié)構(gòu)域樣受體信號通路(NOD-like receptor signaling pathway)、黏著斑等信號通路中顯著富集。MAPK級聯(lián)協(xié)同傳遞各種細胞外信號，從而控制大量不同甚至相反的細胞過程，如增殖、分化、生存、發(fā)育、應(yīng)激反應(yīng)和凋亡[44]，有研究發(fā)現(xiàn)MAPK在調(diào)節(jié)凡納濱對蝦(Litopenaeusvannamei)、香港牡蠣(Crassostreahongkongensis)等多種水生生物的先天免疫過程中發(fā)揮著重要作用[45-46]。同樣在刺參體腔液高溫脅迫的研究中發(fā)現(xiàn)p44/42MAPK可能通過靶向p90RSK參與溫度誘導(dǎo)的代謝抑制[47]。PKB/Akt通路同樣能夠介導(dǎo)多種細胞功能，包括血管生成、代謝、生長、增殖、存活、免疫、轉(zhuǎn)錄和細胞凋亡[48]，在水產(chǎn)動物的研究中發(fā)現(xiàn)Akt基因可能參與介導(dǎo)水產(chǎn)動物機體的免疫防御。如香港牡蠣(Crassostreahongkongensis)在受到溶藻弧菌等病原菌感染后, 其體內(nèi)Akt基因表達上調(diào)[49]。研究發(fā)現(xiàn)海膽中第二大免疫受體家族是NOD-like受體(NLRs)，而在哺乳動物中的NLRs通過充當(dāng)支架蛋白來組裝蛋白復(fù)合物，從而激活NFkB和MAPK信號通路，在免疫應(yīng)答中發(fā)揮作用[6]。在刺參研究中同樣發(fā)現(xiàn)AjNLRP10蛋白參與了對細菌感染的先天免疫應(yīng)答[50]。這些都表明了本研究中miRNA通過靶向這些免疫相關(guān)通路參與刺參響應(yīng)高溫、低氧脅迫下的信號傳導(dǎo)、氧化應(yīng)激以及免疫反應(yīng)等過程。

4 結(jié)語

本研究以高溫組、低氧組和對照組的體腔液細胞為研究對象，利用Illumina測序平臺構(gòu)建了9個miRNA文庫，有效補充了刺參體腔液細胞的miRNA表達譜，并篩選了在高溫脅迫和低氧脅迫下差異表達的miRNAs。基于差異表達的miRNAs，以三代全長轉(zhuǎn)錄組為參考文庫預(yù)測了其候選靶基因，KEGG富集分析表明靶基因主要參與了刺參應(yīng)激、存活和免疫等重要生物學(xué)過程，進一步揭示了差異表達的miRNAs在刺參應(yīng)對高溫脅迫和低氧脅迫過程中的重要調(diào)控作用。后期將利用功能實驗驗證關(guān)鍵miRNAs及其靶基因的互作關(guān)系，進一步探討關(guān)鍵差異miRNA在刺參應(yīng)對高溫、低氧脅迫過程中的調(diào)控作用，以期為刺參應(yīng)對高溫、低氧逆境關(guān)鍵分子標記的篩選和刺參抗逆品種與品系的選育提供理論指導(dǎo)。