低氧—復(fù)氧脅迫下脊尾白蝦鰓組織差異表達(dá)基因趨勢分析

史文軍　王盼　胡潤豪　萬夕和　沈輝　黎慧　王李寶　喬毅　蔣葛　成婕

摘要：【目的】探究脊尾白蝦（Exopalaemon carinicauda）鰓組織在漸變式低氧—復(fù)氧脅迫下的分子調(diào)控機制，為今后開展脊尾白蝦耐低氧品系（種）的選育提供理論指導(dǎo)?！痉椒ā客ㄟ^模擬自然低氧環(huán)境的形成過程，分別于低氧處理0（對照）、3和6 h及復(fù)氧后1和8 h采集脊尾白蝦鰓組織，利用Illumina HiSeq^TM4000測序平臺進(jìn)行轉(zhuǎn)錄組測序分析，經(jīng)過濾和Trinity組裝獲得Unigenes，選取Nr、Swiss-Prot、KEGG和COG/KOG等數(shù)據(jù)庫進(jìn)行注釋分析，在Omicsmart平臺上完成差異表達(dá)基因篩選及其表達(dá)趨勢分析，然后進(jìn)行GO功能注釋分析和KEGG信號通路富集分析，并隨機選取5個差異表達(dá)基因進(jìn)行實時熒光定量PCR驗證。【結(jié)果】脊尾白蝦鰓組織樣本轉(zhuǎn)錄組測序數(shù)據(jù)經(jīng)過濾后的Clean reads進(jìn)行Trinity組裝共獲得93227條Unigenes，其長度范圍在201～35402 bp，平均長度為834 bp，N50長度為1352 bp。通過組間兩兩比較分析，共鑒定出4750個差異表達(dá)基因，其中上調(diào)差異表達(dá)基因3557個、下調(diào)差異表達(dá)基因2829個;超過50%的差異表達(dá)基因被顯著富集到6種基因表達(dá)趨勢模式中（P<0.01），具體表現(xiàn)為：Profile 0模式富集到415個基因，Profile 5模式富集到201個基因，Profile 11模式富集到371個基因，Profile 13模式富集到841個基因，Profile 17模式富集到387個基因，Profile 18模式富集到411個基因。6種基因表達(dá)趨勢模式中的差異表達(dá)基因被注釋到代謝進(jìn)程、細(xì)胞進(jìn)程、單一有機體進(jìn)程、細(xì)胞、細(xì)胞零件、大分子復(fù)合物及催化活性等GO功能條目上;而KEGG信號通路富集分析結(jié)果顯示，以Profile 13模式中的差異表達(dá)基因富集到最多信號通路（86條），其中呈顯著富集的有8條，分別為核糖體、碳代謝、氧化磷酸化、氨基酸生物合成、內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工、糖酵解/糖異生、谷胱甘肽代謝和蛋白輸出?！窘Y(jié)論】脊尾白蝦鰓組織在受低氧脅迫早期通過合成蛋白質(zhì)及提高代謝能力來抵御低氧環(huán)境，隨著低氧脅迫時間的延長，物質(zhì)合成和能量代謝活動均顯著下降;但在復(fù)氧后隨著復(fù)氧時間的延長，其蛋白質(zhì)合成和能量代謝水平又逐漸升高恢復(fù)至常氧水平。

關(guān)鍵詞：脊尾白蝦;低氧—復(fù)氧;鰓組織;轉(zhuǎn)錄組測序;趨勢分析

中圖分類號： S917;S945.49? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼： A 文章編號：2095-1191（2022）03-0735-13

Trend analysis of differentially expressed genes in gill of Exopalaemon carinicauda under hypoxia-reoxygenation stress

SHI Wen-jun WANG Pan HU Run-hao WAN Xi-he SHEN Hui LI Hui WANG Li-bao QIAO Yi JIANG Ge CHENG Jie

（1Institute of Oceanology and Marine Fisheries， Nantong， Jiangsu? 226007， China; 2Yantai Institute of Coastal Zone Research， Chinese Academy of Sciences， Yantai， Shandong? 264003， China;? 3National Demonstration Center for Experimental Fisheries Science Education（Shanghai Ocean University）， Shanghai? 201306， China）

Abstract：【Objective】To investigate the molecular regulation mechanisms in gill tissue of Exopalaemon carinicauda under gradually changing hypoxia-reoxygenation stress， so as to provide theoretical reference for the breeding of hypoxia tolerance strains （species） of E. carinicauda. 【Method】By simulating the formation process of hypoxia in the natural environment， gill tissues were collected at 0 （control）， 3 and 6 h， and 1 and 8 h after reoxygenation， respectively. Transcriptome sequencing analysis was performed using Illumina HiSeq^TM4000 sequencing platform. Unigenes were obtained by filtration and Trinity assembly. Nr， Swiss-Prot， KEGG and COG/KOG databases were selected for annotation analysis. The differentially expressed genes （DEGs） were screened and their trends were analyzed on the Omicsmart platform. Then GO functional annotation analysis and KEGG signal pathway enrichment analysis were performed， and five DEGs were randomly selected for real-time quantitative PCR（qRT-PCR） validation. 【Result】A total of 93227 uni genes whose range were from 201 to 35402 bp were obtained through transcriptome sequencing， with their average length of 834 bp and the N50 of 1352 bp. 4750 genes were identified as DEGs through pairwise comparison between groups， of which the up-regulated genes were 3557 and the down-regulated genes were 2829. More than fifty percent were significantly enriched in six gene expression trends （P<0.01）：Profile 0 was enriched in 415 genes;Profile 5 was enriched in 201 genes;Profile 11 was enriched in 371 genes;Profile 13 was enriched in 841 genes;Profile 17 was enriched in 387 genes and Profile 18 was enriched in 411 genes.The DEGs in six trends were annotated to metabolic processes， cellular processes， single organism process， cell， cell parts， macromolecular complex and catalytic activity by GO functional analysis. The KEGG signaling pathway enrichment analysis showed that DEGs in Profile 13 were enriched in 86 pathways at most and significantly enriched? in ribosome，carbon metabolism，oxidative phosphorylation，biosynthesis of amino acids，protein processing in endoplasmic reticulum，glycolysis or gluconeogenesis，glutathione metabolism and protein export. 【Conclusion】The gill tissues of E.carinicauda synthesizes protein and increases metabolic capacity to resist hypoxia in the early stage of hypoxia stress. But with the prolongation of hypoxia， both substance synthesis and energy metabolism were affected and decreased. After reoxygenation， with the prolongation ofreoxygenation， the protein synthesis and energy metabolism gradually recover to level under normal oxygen.038C597C-F774-40C1-9DC1-F6FE2E5B0876

Key words： Exopalaemon carinicauda; hypoxia-reoxygenation; gilltissue; transcriptome sequencing; trend analyses

Foundation items： Agricultural Major New Breed Creation of Jiangsu Province （PZCZ201747）;Agricultural Scien-ce and Technology Independent Innovation Project of Jiangsu Province （CX18-2010）; The Fifth “226 Project” Cultivation Fund Support Project of Nantong （NT226-202021）; Science and Technology Plan Project of Nantong （JC2019057）

0 引言

【研究意義】脊尾白蝦（Exopalaemon carinicauda）廣泛分布于我國黃海、渤海沿岸海域，是一種重要的海洋漁業(yè)和池塘養(yǎng)殖品種（Ma et al.，2020），因具有良好的口感及豐富的營養(yǎng)而深受消費者喜愛，加之生長速度快、繁殖能力強及環(huán)境適應(yīng)性廣等優(yōu)點，其養(yǎng)殖規(guī)模逐年擴(kuò)大，已發(fā)展成為池塘單養(yǎng)和魚蝦蟹貝混養(yǎng)的重要經(jīng)濟(jì)蝦類（李玉全，2014）。據(jù)統(tǒng)計，我國每年的脊尾白蝦養(yǎng)殖面積和產(chǎn)量分別為2萬ha和4.5萬t（劉俊杰等，2018;Shao et al.，2018;Li et al.，2019）。為了追求高產(chǎn)量，當(dāng)前的池塘養(yǎng)殖密度普遍較高，但高密度養(yǎng)殖易導(dǎo)致池底殘餌和排泄物過度累積，而脊尾白蝦大部分時間生活在池塘底層，由于呼吸耗氧作用及積累的有機物分解，極易出現(xiàn)低氧情況，尤其是在晚上和清晨，低氧嚴(yán)重時會威脅到脊尾白蝦的生命（Cheng et al.，2003）。除了依靠開放式曝氣維持水中的含氧量外，目前尚無其他更好的解決方法，但這種方法不僅增加養(yǎng)殖成本，還無法完全解決低氧的問題。因此，亟待明確脊尾白蝦對低氧的適應(yīng)機制，培育出耐低氧的新品系（種）以解決這一突出問題。【前人研究進(jìn)展】至今，有關(guān)低氧對蝦類影響的研究主要集中在生長發(fā)育（胡賢德和林北堃，1999）、生活習(xí)性（陳琴等，2001）、免疫能力（Cheng et al.，2002）、相關(guān)基因克隆（Li et al.，2009）、代謝過程（Dupont-Prinet et al.，2013）、行為變化（侯文杰等，2014）及組織損傷（楊明等，2019）等方面。近年來，隨著新一代轉(zhuǎn)錄組測序技術(shù)的迅速發(fā)展，轉(zhuǎn)錄組分析已廣泛用于揭示逆境脅迫下對蝦的分子生物學(xué)過程和生理功能變化機制（Wang et al.，2020）。Hu等（2015）對低鹽脅迫下凡納濱對蝦鰓組織進(jìn)行轉(zhuǎn)錄組分析，獲得了一些與離子和滲透調(diào)節(jié)相關(guān)的基因;Huang等（2017）對凡納濱對蝦低溫脅迫不同時間和低溫脅迫后恢復(fù)的肌肉組織進(jìn)行轉(zhuǎn)錄組測序分析，結(jié)果發(fā)現(xiàn)絲氨酸/蘇氨酸—蛋白激酶信號通路在冷適應(yīng)中可能發(fā)揮著更重要作用;Li等（2018）對氨脅迫后的斑節(jié)對蝦肝胰腺組織進(jìn)行轉(zhuǎn)錄組分析，證實氨應(yīng)激誘導(dǎo)氧化應(yīng)激，最終導(dǎo)致細(xì)胞凋亡。然而，通過轉(zhuǎn)錄組學(xué)挖掘分析蝦類低氧脅迫過程中差異表達(dá)基因的研究相對較少。Sun等（2014）通過比較分析低氧脅迫下日本沼蝦肝胰腺組織中基因的表達(dá)情況，發(fā)現(xiàn)差異表達(dá)基因主要與呼吸代謝、血藍(lán)蛋白合成、能量消耗、抗氧化能力和免疫應(yīng)答相關(guān)。曹梅等（2020）通過研究低氧脅迫下脊尾白蝦全蝦組織中相關(guān)基因的表達(dá)情況，證實缺氧能促使蝦體糖酵解等無氧代謝途徑增強?！颈狙芯壳腥朦c】現(xiàn)有的研究報道主要是將蝦類直接從常氧水體轉(zhuǎn)移至低氧水體中或分析常氧與低氧極限組的轉(zhuǎn)錄學(xué)，其低氧處理方式可能無法準(zhǔn)確反映實際養(yǎng)殖環(huán)境中脊尾白蝦經(jīng)歷的低氧形成過程：常氧—溶解氧降低—低氧—極限低溶解氧—死亡。鰓組織是蝦類呼吸的重要場所，但至今有關(guān)脊尾白蝦鰓組織在漸變式低氧—復(fù)氧脅迫下的適應(yīng)性調(diào)節(jié)機制研究鮮見報道?！緮M解決的關(guān)鍵問題】通過模擬自然低氧環(huán)境的形成過程，對漸變式低氧—復(fù)氧不同時間點的脊尾白蝦鰓組織進(jìn)行轉(zhuǎn)錄組測序分析，旨在探討脊尾白蝦對低氧—復(fù)氧響應(yīng)的分子機制，為今后開展脊尾白蝦耐低氧品系（種）的選育提供理論指導(dǎo)。

1 材料與方法

1. 1 試驗材料

試驗在江蘇省海洋水產(chǎn)研究所如東海水蝦類科研基地進(jìn)行，供試用蝦為取自基地外塘飼養(yǎng)的脊尾白蝦。選取體長5.3±0.1 cm，體重2.57±0.21 g，無外傷、活力好、健康的脊尾白蝦，置于1000 L的PVC圓桶中暫養(yǎng)7 d。暫養(yǎng)期間的海水鹽度為25‰、水溫20～21 ℃，持續(xù)充氧，每天投喂適量的商品配合飼料。正式試驗前1 d停止喂食;每天換水1次，換水量約為總水量的1/3。

1. 2 試驗設(shè)計及樣本采集

以預(yù)試驗中的脊尾白蝦存活量和水體溶解氧為參考，將密封低氧時間定為6 h，然后立即開始充氣復(fù)氧。正式試驗在鹽度25‰、水溫20～21 ℃的條件下進(jìn)行，將暫養(yǎng)穩(wěn)定后的脊尾白蝦隨機分裝到裝有5 L海水的三角瓶中，每瓶10尾，每3瓶為1個組，共5組，充氣適應(yīng)0.5 h后，停止充氣，并用帶有橡膠導(dǎo)管的橡膠塞密封，開始計時。分別于低氧處理0（對照）、3和6 h及復(fù)氧后1和8 h采集水樣和脊尾白蝦鰓組織。每瓶取1個水樣，每組3個重復(fù)，以堿性碘化鉀和氯化錳固定液固定，用于后續(xù)的水體溶解氧測定;取完水樣后立即從每瓶中選取5～6尾活蝦的完整鰓組織作為1個生物學(xué)重復(fù)，每組3個重復(fù)，樣品編號對應(yīng)為js-0、js-3、js-6、js-7和js-14，液氮速凍后置于 -80 ℃超低溫冰箱保存，用于總RNA提取。

1. 3 溶解氧測定

參照GB 17378.4—2007《海洋監(jiān)測規(guī)范 第4部分：海水分析》中的碘量法測定水體溶解氧含量。038C597C-F774-40C1-9DC1-F6FE2E5B0876

1. 4 RNA提取、cDNA文庫構(gòu)建及高通量測序

按照TRIzol試劑盒（美國Invitrogen公司）說明提取樣本總RNA，分別使用Agilent 2100生物分析儀評估RNA質(zhì)量、1.0%瓊脂糖凝膠電泳檢測其完整性、NanoDrop微量分光光度計測定其純度?？俁NA檢測合格后采用帶有Oligo（dT）的磁珠富集mRNA，然后利用片段緩沖液將富集的mRNA裂解成短片段。以片段化的mRNA為模板、隨機寡核苷酸為引物，在ProtoScript II逆轉(zhuǎn)錄酶體系中合成cDNA第一條鏈;隨后采用RNaseH降解RNA鏈，在DNA Polymerase I體系下以dNTPs為原料合成cDNA第二條鏈。純化后的雙鏈cDNA經(jīng)末端修復(fù)及加poly（A）尾，并連接測序接頭，以AMPure XP beads篩選200 bp左右的cDNA進(jìn)行PCR擴(kuò)增，使用AMPure XP beads純化PCR產(chǎn)物以構(gòu)建cDNA文庫，并通過Illumina HiSeq^TM4000測序平臺進(jìn)行高通量測序。其中，cDNA文庫構(gòu)建及高通量測序均委托廣州基迪奧生物科技有限公司完成。

1. 5 轉(zhuǎn)錄組組裝與Unigenes功能注釋

測序獲得的原始序列（Raw reads）利用fastp（0.18.0）進(jìn)行數(shù)據(jù)質(zhì)量控制（Chen et al.，2018），過濾低質(zhì)量數(shù)據(jù)，具體包括：（1）去除含Adapter的Reads;（2）去除含N比例大于10%的Reads;（3）去除全部都是A堿基的Reads;（4）去除低質(zhì)量Reads（Q≤20堿基數(shù)占整條Reads的50%以上），剩余的有效序列（Clean reads）使用Trinity進(jìn)行組裝（Grabherr et al.，2011）。使用BLASTx程序?qū)⑺薪M裝獲得的序列分別在Nr、Swiss-Prot、KEGG及COG/KOG等數(shù)據(jù)庫中進(jìn)行比對，設(shè)定閾值1e-5，根據(jù)最佳比對結(jié)果獲得序列注釋信息。

1. 6 差異表達(dá)基因鑒定及其表達(dá)趨勢分析

基因表達(dá)量采用RPKM法進(jìn)行計算（Mortazavi et al.，2008），以DESeq2進(jìn)行兩組間基因表達(dá)差異分析（Love et al.，2014），將錯誤發(fā)現(xiàn)率（FDR）<0.01和差異倍數(shù)絕對值（|Fold Change|）≥2作為差異基因篩選標(biāo)準(zhǔn)。然后對所有差異表達(dá)基因進(jìn)行趨勢分析（Ernst and Bar-Joseph，2006），設(shè)定模塊數(shù)量為20，P為0.01，以獲得顯著性和非顯著性基因表達(dá)模式。差異表達(dá)基因篩選及其趨勢分析均在Omicsmart平臺（https：//www.omicsmart.com）上完成。

1. 7 實時熒光定量PCR驗證

為驗證轉(zhuǎn)錄組測序結(jié)果的可靠性，隨機挑選5個差異表達(dá)基因進(jìn)行實時熒光定量PCR驗證，以18S為內(nèi)參基因，以樣本RNA反轉(zhuǎn)錄合成的cDNA為模板，按照TB Green試劑盒操作說明進(jìn)行實時熒光定量PCR驗證（引物見表1），設(shè)3次重復(fù)。反應(yīng)體系20.0 μL：TB Green Premix Ex Taq 10.0 μL，10 μmol/L正、反向引物各0.4 μL，ROX Reference Dye （50×） 0.4 μL，cDNA模板2.0 μL，DEPC水6.8 μL。在Step-One Plus熒光定量PCR儀上進(jìn)行擴(kuò)增，擴(kuò)增程序：95 ℃預(yù)變性30 s;95 ℃ 5 s，60 ℃ 30 s，進(jìn)行40個循環(huán);熔解曲線設(shè)為95 ℃ 15 s，60 ℃ 1 min，95 ℃ 15 s?；蛳鄬Ρ磉_(dá)量以2-ΔΔCt 法進(jìn)行換算。

2 結(jié)果與分析

2. 1 水樣溶解氧含量變化趨勢

低氧處理0 h（對照）的水樣溶解氧含量為8.04±0.08 mg/L，低氧處理3 h的水樣溶解氧含量為2.12±0.20 mg/L，即此時已達(dá)低氧環(huán)境;低氧處理6 h的水樣溶解氧含量為1.04±0.12 mg/L，溶解氧繼續(xù)降低，低氧脅迫更加劇烈;充氣復(fù)氧后1和8 h的水樣溶解氧含量分別為8.10±0.28和8.05±0.15 mg/L，已恢復(fù)至對照水樣溶解氧水平。說明脊尾白蝦在此試驗過程中已經(jīng)歷了由常氧環(huán)境逐漸到低氧環(huán)境再恢復(fù)到常氧環(huán)境的脅迫過程。

2. 2 轉(zhuǎn)錄組數(shù)據(jù)及其組裝情況

高通量測序結(jié)果（表2）顯示，每個脊尾白蝦鰓組織樣本獲得39381904～64573816條Raw reads，經(jīng)數(shù)據(jù)過濾處理后，平均在每個樣本獲得52394125條Clean reads，平均有效率為99.84%。Q20（測序堿基正確識別率在99.0%以上）大于98.00%，Q30（測序堿基正確識別率在99.9%以上）大于94.00%，且過濾后樣本的平均GC含量為47.26%，表明測序質(zhì)量較好。過濾后的Clean reads經(jīng)Trinity組裝后共獲得93227條Unigenes（表3），其長度范圍在201～35402 bp（圖1），平均長度為834 bp，N50長度為1352 bp，表明組裝效果較良好。

2. 3 Unigenes功能注釋分析結(jié)果

分別在Nr、Swiss-Prot、KEGG和COG/KOG數(shù)據(jù)庫中對Unigenes進(jìn)行注釋分析，結(jié)果（表4）顯示，在Nr數(shù)據(jù)庫中注釋到37346條Unigenes，在Swiss-Prot數(shù)據(jù)庫中注釋到29042條Unigenes，在KEGG數(shù)據(jù)庫中注釋到33955條Unigenes，在COG/KOG數(shù)據(jù)庫中注釋到24238條Unigenes，累計被注釋的Unigenes為40719條，注釋率為43.68%。

2. 4 樣本重復(fù)性分析結(jié)果

Pearson相關(guān)系數(shù)熱圖（圖2）顯示，除js-6-2樣本外，其他樣本在組內(nèi)的生物學(xué)重復(fù)相關(guān)性均在0.80以上（0.81～0.99），相關(guān)性較好。造成js-6-2樣本相關(guān)性偏低的原因可能是取樣、提取或測序等某一環(huán)節(jié)發(fā)生污染，因此后續(xù)分析與驗證均將該樣本測序數(shù)據(jù)剔除。038C597C-F774-40C1-9DC1-F6FE2E5B0876

2. 5 差異表達(dá)基因鑒定及趨勢分析結(jié)果

通過組間兩兩比較分析，共鑒定出4750個差異表達(dá)基因（圖3），其中上調(diào)差異表達(dá)基因3557個、下調(diào)差異表達(dá)基因2829個。由圖4可看出，與對照組相比，隨著低氧處理時間的延長，差異表達(dá)基因的數(shù)量呈增長趨勢。差異表達(dá)基因最多的出現(xiàn)在js-7 vs js-14，其中上調(diào)差異表達(dá)基因1940個、下調(diào)差異表達(dá)基因815個;而js-6 vs js-7鑒定出的差異表達(dá)基因最少（上調(diào)差異表達(dá)基因240個、下調(diào)差異表達(dá)基因360個），說明低氧處理后再復(fù)氧對脊尾白蝦鰓組織基因表達(dá)影響顯著。

對4750個差異表達(dá)基因進(jìn)行表達(dá)趨勢分析，結(jié)果（圖5）顯示在默認(rèn)的20種基因表達(dá)趨勢模式中有6種顯著的基因表達(dá)趨勢模式（P<0.01）。其中，Profile 0模式表示基因表達(dá)量在低氧和復(fù)氧時一直下降（415個）;Profile 5模式表示基因表達(dá)量在低氧時先下降后上升，復(fù)氧后再下降（201個）;Profile 11和Profile 17模式表示基因表達(dá)量在低氧時上升，復(fù)氧后呈先下降再上升的變化趨勢（371和387個）;Profile 13模式表示基因表達(dá)量在低氧時先上升后下降，復(fù)氧后再上升（841個）;Profile 18模式表示基因表達(dá)量在低氧時上升，復(fù)氧后下降（411個）。綜合分析發(fā)現(xiàn)，超過50%的差異表達(dá)基因被富集到這6種基因表達(dá)趨勢模式中。

2. 6 差異表達(dá)基因GO功能注釋分析結(jié)果

GO功能注釋分析結(jié)果（圖6）顯示，顯著富集在6種基因表達(dá)趨勢模式中的差異表達(dá)基因主要注釋在代謝進(jìn)程（Metabolic processes）、細(xì)胞進(jìn)程（Cellular processes）、單一有機體進(jìn)程（Single-organism process）、細(xì)胞（Cell）、細(xì)胞零件（Cell part）、大分子復(fù)合物（Macromolecular complex）、催化活性（Catalytic activity）及結(jié)合（Binding）等GO功能條目上。

2. 7 差異表達(dá)基因KEGG信號通路富集分析結(jié)果

KEGG信號通路富集分析結(jié)果（表5）顯示，在Profile 0模式中富集到42條信號通路，呈顯著富集的有6條，分別為：①細(xì)胞外基質(zhì)受體相互作用（ECM-receptor interaction）;②脂肪酸降解（Fatty acid degradation）;③溶酶體（Lysosome）;④纈氨酸，亮氨酸和異亮氨酸的生物合成（Valine，leucine and isoleucine biosynthesis）;⑤甘氨酸，絲氨酸和蘇氨酸代謝（Glycine，serine and threonine metabolism）;⑥類固醇生物合成（Steroid biosynthesis）。在Profile 5模式中富集到41條信號通路，呈顯著富集的有6條，分別為：①MAPK信號通路（MAPK signaling pathway-fly）;②泛醌及其他萜醌類生物合成（Ubiquinone and other terpenoid-quinone biosynthesis）;③凋亡（Apoptosis-fly）;④自噬（Autophagy-other eukaryotes）;⑤自噬（Autophagy-animal）;⑥生理節(jié)律（Circadian rhythm-fly）。在Profile 11模式中富集到4條信號通路，呈顯著富集的有3條，分別為：①花生四烯酸代謝（Arachidonic acid metabolism）;②ABC轉(zhuǎn)運蛋白（ABC transporters）;③蛋白酶體（Proteasome）。在Profile 13模式中富集到86條信號通路，呈顯著富集的有8條（圖7），分別為：①核糖體（Ribosome）;②碳代謝（Carbon metabolism）;③氧化磷酸化（Oxidative phosphorylation）;④氨基酸生物合成（Biosynthesis of amino acids）;⑤內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工（Protein processing in endoplasmic reticulum）;⑥糖酵解/糖異生（Glycolysis/gluconeogenesis）;⑦谷胱甘肽代謝（Glutathione metabolism）;⑧蛋白輸出（Protein export）。在Profile 17模式中富集到17條信號通路，呈顯著富集的有4條，分別為：①核糖體（Ribosome）;②RNA降解（RNA degradation）;③嘧啶代謝（Pyrimidine metabolism）;④壽命調(diào)節(jié)途徑（Longevity regulating pathway-worm）。在Profile 18模式中富集到65條信號通路，呈顯著富集的有3條，分別為：①TGF-β信號通路（TGF-beta signaling pathway）;②MAPK信號通路（MAPK signaling pathway）;③纈氨酸、亮氨酸和異亮氨酸的生物合成（Valine，leucine and isoleucine biosynthesis）。

2. 8 實時熒光定量PCR驗證結(jié)果

隨機挑選5個差異表達(dá)基因（Gene4007、Gene 29395、Gene40966、Gene56231和Gene82609）進(jìn)行實時熒光定量PCR驗證，結(jié)果（圖8）表明其相對表達(dá)量的變化趨勢與轉(zhuǎn)錄組測序獲得的基因表達(dá)量變化趨勢基本一致，表明本研究的轉(zhuǎn)錄組測序結(jié)果真實可靠。

3 討論

3. 1 低氧—復(fù)氧對脊尾白蝦呼吸耗氧及存活的影響

溶解氧是影響水生動物呼吸的重要環(huán)境因子，與耗氧量直接相關(guān)（曹祥德等，2014;Ulaje et al.，2019）。一般情況下，甲殼類動物可通過自身調(diào)節(jié)使機體處于最小耗氧量的狀態(tài)，以此耐受一定程度的低氧環(huán)境，但超出其耐受極限則會引起死亡（Larkin et al.，2008）。本研究發(fā)現(xiàn)，在脊尾白蝦自身呼吸耗氧的情況下，隨著低氧脅迫時間的延長，水體溶解氧含量逐漸降低，但低氧處理3～6 h的水樣溶解氧含量降低幅度明顯低于低氧處理0～3 h的水樣，表明隨著低氧脅迫時間的延長，脊尾白蝦的耗氧率在逐漸降低。王愛敏（1997）研究表明羅氏沼蝦蚤狀幼體的耗氧速率隨水體溶解氧含量的降低而下降;單保黨等（2018）也研究發(fā)現(xiàn)凡納濱對蝦p12～p13階段仔蝦的耗氧率隨水體溶解氧含量的降低而下降，故推測這種通過降低耗氧率的調(diào)節(jié)方式是蝦類普遍應(yīng)對低氧環(huán)境的一種策略。本研究還發(fā)現(xiàn)，低氧處理6 h時的水體溶解氧含量降至1.04 mg/L，此時脊尾白蝦的死亡率為47%。李明云等（1992）研究發(fā)現(xiàn)平均體長為5.02 cm的脊尾白蝦在水溫（18.0±0.5）℃中死亡50%時的溶解氧含量為0.78 mg/L;曹梅等（2020）研究證實平均體重為3.15±0.26 g的脊尾白蝦在水溫22～24 ℃、溶解氧含量為1.13±0.26 mg/L的條件下處于昏厥狀態(tài)?？梢?，蝦體規(guī)格和環(huán)境條件均會影響脊尾白蝦對低氧的耐受程度。038C597C-F774-40C1-9DC1-F6FE2E5B0876

3. 2 低氧—復(fù)氧對脊尾白蝦鰓組織基因表達(dá)的影響

鰓組織是甲殼動物進(jìn)行呼吸、排泄及調(diào)節(jié)滲透壓的重要器官，其直接與水體接觸，易受水環(huán)境影響（盧建平，2001）。低氧會影響蝦類鰓組織的正常生理功能，進(jìn)而引起能量代謝及其他物質(zhì)合成代謝紊亂、抗氧化應(yīng)激損傷等（Paschke et al.，2010;Sun et al.，2015），但脊尾白蝦鰓組織如何適應(yīng)低氧—復(fù)氧的分子調(diào)控機制尚未明確，因此利用轉(zhuǎn)錄組測序分析其鰓組織在漸變式低氧—復(fù)氧過程中的調(diào)節(jié)機制具有重要意義。趨勢分析可從聚類結(jié)果中挑選符合一定生物學(xué)特性的基因集，從而降低數(shù)據(jù)的復(fù)雜性和分析難度（郭睿等，2018）。為此，本研究基于轉(zhuǎn)錄組測序分析挖掘低氧—復(fù)氧脅迫下各時間點脊尾白蝦鰓組織的差異表達(dá)基因，并對所有差異表達(dá)基因進(jìn)行表達(dá)趨勢分析，GO功能注釋分析結(jié)果顯示，顯著富集在6種基因表達(dá)趨勢模式中的差異表達(dá)基因被注釋到生物學(xué)進(jìn)程（Biological process）、細(xì)胞組分（Cellular component）和分子功能（Molecular function）三大功能中，主要涉及代謝進(jìn)程、細(xì)胞進(jìn)程、單一有機體進(jìn)程、細(xì)胞、細(xì)胞零件、大分子復(fù)合物及催化活性等，說明脊尾白蝦鰓組織在低氧—復(fù)氧脅迫下的代謝和細(xì)胞活動較活躍。而KEGG信號通路富集分析結(jié)果顯示，以Profile 13模式中的差異表達(dá)基因富集到最多信號通路（86條），其中呈顯著富集的有8條，分別為核糖體、碳代謝、氧化磷酸化、氨基酸生物合成、內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工、糖酵解/糖異生、谷胱甘肽代謝和蛋白輸出。

3. 2. 1 對脊尾白蝦鰓組織蛋白合成相關(guān)基因的影響

核糖體主要由核糖體RNA（rRNA）、核糖體蛋白（RPs）和小核糖體RNA（snoRNAs）組成，是蛋白翻譯機制的重要組成部分（Nakhoul et al.，2014），其中RPs在rRNA加工過程中促進(jìn)rRNA折疊形成功能三維結(jié)構(gòu)及穩(wěn)定核糖體最終空間構(gòu)象方面發(fā)揮重要作用（Xie et al.，2018）。RPs除了在蛋白翻譯中發(fā)揮作用外，還參與核糖體外的DNA修復(fù)、細(xì)胞凋亡和細(xì)胞內(nèi)穩(wěn)態(tài)等功能（Shenoy et al.，2012）。內(nèi)質(zhì)網(wǎng)是負(fù)責(zé)跨膜蛋白和分泌蛋白合成與折疊的細(xì)胞器（Iurlaro and Munoz-Pinedo，2016）。本研究發(fā)現(xiàn)，有80個差異表達(dá)基因富集在核糖體通路上，包括核糖體大亞基蛋白RPLP0、RPLP1和RPLP2及核糖體小亞基蛋白RPS10e、RPS11e和RPS12e;有19個差異表達(dá)基因富集在內(nèi)質(zhì)網(wǎng)蛋白加工通路上，如蛋白質(zhì)二硫鍵異構(gòu)酶A3（PDIA3）、鈣網(wǎng)蛋白（CALR）和X-盒結(jié)合蛋白1（XBP1）。此外，有17個差異表達(dá)基因富集在氨基酸生物合成通路上，有5個差異表達(dá)基因富集在蛋白輸出通路上。Land等（1993）研究發(fā)現(xiàn)西部錦龜（Chrysemys picta Bellii）肝細(xì)胞蛋白質(zhì)合成速率在25 ℃缺氧12 h后下降92%，但能在復(fù)氧2 h后恢復(fù)正常。Smith等（1996）研究發(fā)現(xiàn)，鯽魚在缺氧狀態(tài)暴露48 h會導(dǎo)致其肝臟、心臟和肌肉的蛋白質(zhì)合成率下降，但腦組織的蛋白合成率未受影響，在復(fù)氧24 h后各組織的蛋白合成率基本恢復(fù)至常氧組水平。本研究結(jié)果表明，這些通路上的差異表達(dá)基因表達(dá)量在低氧時均隨溶解氧含量的降低而呈先升高后降低的變化趨勢，故推測脊尾白蝦鰓組織在低氧早期蛋白合成較活躍，但溶解氧含量降至低氧水平后蛋白合成能力減弱;在復(fù)氧后1 h內(nèi)蛋白合成的能力尚未恢復(fù)到常氧水平，隨著復(fù)氧時間的延長，蛋白合成能力逐漸恢復(fù)并有所升高。

3. 2. 2 對脊尾白蝦鰓組織能量代謝相關(guān)基因的影響

呼吸代謝是生物能量代謝研究的重要內(nèi)容之一，不僅能反映生物的生理狀態(tài)，還可表征環(huán)境條件對生物生理活動的影響（馮雪等，2012）。本研究在Profile 13表達(dá)模式中分別發(fā)現(xiàn)有24、19和11個差異表達(dá)基因富集在碳代謝、氧化磷酸化和糖酵解/糖異生通路上，包括蘋果酸酶（MDH）、異檸檬酸脫氫酶（IDH）、6-磷酸葡糖酸脫氫酶（PGDH）、甘油醛-3-磷酸脫氫酶（GAPDH）和果糖-二磷酸醛縮酶（ALDO）等重要呼吸代謝酶類。MDH催化蘋果酸鹽和草酰乙酸鹽的可逆氧化還原反應(yīng)及蘋果酸鹽氧化脫羧成丙酮酸鹽，且在有氧糖酵解過程中參與蘋果酸天冬氨酸穿梭，其電子受體為NAD+或NADP+（Soldatov et al.，2020）;IDH催化異檸檬酸氧化脫羧生成α-酮戊二酸，在此過程中也利用NAD+或NADP+生產(chǎn)NADH或NADPH（Kim and Park，2010）;PGD催化6-磷酸葡萄糖酸產(chǎn)生核酮糖-5-磷酸和NADPH;GAPDH催化D-甘油醛-3-磷酸（GAP）在有NAD+和無機磷酸鹽存在的條件下可逆氧化和磷酸化成1，3-二磷酸甘油酸;ALDO催化3-磷酸甘油醛和磷酸二羥丙酮的醇醛縮合反應(yīng)可逆生成果糖-1，6-二磷酸（王志剛等，2019）。這些呼吸代謝酶在糖酵解、三羧酸循環(huán)、磷酸戊糖途徑及糖異生途徑中發(fā)揮重要作用，為生物體物質(zhì)合成代謝提供能量與底物。已有較多研究證實，低氧—復(fù)氧對水生動物鰓能量代謝有影響，但不同水生動物的響應(yīng)模式可能存在差異（Crocker et al.，2013;Li et al.，2017）。

3. 2. 3 對脊尾白蝦鰓組織抗氧化酶相關(guān)基因的影響

本研究發(fā)現(xiàn)有11個差異表達(dá)基因富集到谷胱甘肽代謝通路上。谷胱甘肽（GSH）能維持機體正常的免疫系統(tǒng)功能，具有抗氧化作用，而NADPH也是抗氧化劑的關(guān)鍵電子供體（Saetan et al.，2020）。此外，本研究中的谷胱甘肽過氧化物酶（GPX）和谷胱甘肽-S-轉(zhuǎn)移酶（GST）基因表達(dá)變化趨勢與王盼等（2020）研究脊尾白蝦鰓組織在低氧—復(fù)氧條件下酶活性變化的趨勢相似，進(jìn)一步說明這2個基因的表達(dá)變化可能在低氧—復(fù)氧的抗氧化損傷中發(fā)揮重要作用。038C597C-F774-40C1-9DC1-F6FE2E5B0876

4 結(jié)論

脊尾白蝦鰓組織在受低氧脅迫早期通過合成蛋白質(zhì)及提高代謝能力來抵御低氧環(huán)境，隨著低氧脅迫時間的延長，物質(zhì)合成和能量代謝活動均顯著下降;但在復(fù)氧后隨著復(fù)氧時間的延長，其蛋白質(zhì)合成和能量代謝水平又逐漸升高恢復(fù)至常氧水平。

參考文獻(xiàn)：

曹祥德，張根玉，范武江，李雪松. 2014. 溫度與溶氧水平對日本沼蝦瞬時耗氧速率的影響[J]. 廣東農(nóng)業(yè)科學(xué)，41（12）：123-126. [Cao X D，Zhang G Y，F(xiàn)an W J，Li X S. 2014. Effects of temperature and dissolved oxygen level on instantaneous rate of oxygen consumption of Macrobrachium nipponense[J]. Guangdong Agricultural Scien-ces，41（12）：123-126.] doi：10.3969/j.issn.1004-874X. 2014.12.028.

曹梅，王興強，秦傳新，沈曄，張子楊，錢詩悅. 2020. 脊尾白蝦對低氧響應(yīng)的轉(zhuǎn)錄組學(xué)分析[J]. 漁業(yè)科學(xué)進(jìn)展，41（2）：112-123. [Cao M，Wang X Q，Qin C X，Shen Y，Zhang Z Y，Qian S Y. 2020. Transcriptome analysis of Palaemon carincauda subjected to hypoxic stress[J]. Progress in Fishery Sciences，41（2）：112-123.] doi：10.19663/j.issn2095-9869.20190924001.

陳琴，陳曉漢，羅永巨，黃鈞，李貴福，寧良坤. 2001. 南美白對蝦耗氧率和窒息點的初步測定[J]. 水利漁業(yè)，21（2）：14-15. [Chen Q，Chen X H，Luo Y J，Huang J，Li G F，Ning L K. 2001. Determination of the oxygen consumption rate and suffocation point of Litopenaeus vannamei[J]. Reservoir Fisheries，21（2）：14-15.] doi：10.15928/j.1674-3075. 2001.02.008.

馮雪，陳丕茂，秦傳新. 2012. 漁業(yè)資源生物呼吸代謝研究進(jìn)展[J]. 廣東農(nóng)業(yè)科學(xué)，39（4）：92-96. [Feng X，Chen P M，Qin C X. 2012. Advancement of the study on respiratory metabolism of fishery resources organism[J]. Guangdong Agricultural Sciences，39（4）：92-96.] doi：10.16768/j.issn.1004-874x.2012.04.033.

郭睿，解彥玲，熊翠玲，尹偉軒，鄭燕珍，付中民，陳大福. 2018. 意大利蜜蜂4、5和6日齡幼蟲腸道發(fā)育過程中差異表達(dá)基因的趨勢分析[J]. 上海交通大學(xué)學(xué)報（農(nóng)業(yè)科學(xué)版），36（4）：14-21. [Guo R，Xie Y L，Xiong C L，Yin W X，Zheng Y Z，F(xiàn)u Z M，Chen D F. 2018. Trend analysis for differential expressed genes in the developmental process of 4-，5- and 6- day-old larval guts of Apis mellifera ligustica[J]. Journal of Shanghai Jiaotong University （Agricultural Science），36（4）：14-21.] doi：10.3969/J.ISSN. 1671-9964.2018.04.003.

侯文杰，周文玉，潘桂平，畢浩，周裕華，劉本偉，張年國. 2014. 溶氧水平對日本對蝦耗氧能力與生理行為的影響[J]. 江蘇農(nóng)業(yè)科學(xué)，42（2）：190-192. [Hou W J，Zhou W Y，Pan G P，Bi H，Zhou Y H，Liu B W，Zhang N G. 2014. Effect of dissolved oxygen level on the oxygen consumption capacity and physiological behavior of Penaeus japonicus[J]. Jiangsu Agricultural Sciences，42（2）：190-192.] doi：10.3969/j.issn.1002-1302.2014.02.066.

胡賢德，林北堃. 1999. 水體溶氧量對對蝦消化道蛋白酶活力的影響[J]. 沈陽農(nóng)業(yè)大學(xué)學(xué)報，30（4）：457-459. [Hu X D，Lin B K. 1999. Effect of oxygen content in water on proteases vitality in prawn digestive system[J]. Journal of Shenynag Agricultural University，30（4）：457-459.]

李明云，包堅敏，吳春娥. 1992. 脊尾白蝦窒息點與耗氧率的試驗觀察[J]. 海洋漁業(yè)，（6）：251-253. [Li M Y，Bao J M，Wu C E. 1992. Experimental observations of asphy-xiation point and oxygen consumption rate in ridgetail white prawn，Exopalaemon carinicauda[J]. Marine Fishe-ries，（6）：251-253.]038C597C-F774-40C1-9DC1-F6FE2E5B0876

李玉全. 2014. pH脅迫對脊尾白蝦代謝酶活力的影響[J]. 南方農(nóng)業(yè)學(xué)報，45（6）：1098-1101. [Li Y Q. 2014. Effects of pH stress on activities of phosphatase in Exopalaemon carinicauda Holthuis[J]. Journal of Southern Agriculture，45（6）：1098-1101.] doi：10.3969/j：issn.2095-1191.2014. 6.1098.

劉俊杰，肖李霞，陶冶，孟玲玲，朱善央，陸錦天，沈新強，王云龍，馬凌波，顧邦祥. 2018. 脊尾白蝦引種及多品種混養(yǎng)集成與示范[J]. 科學(xué)養(yǎng)魚，（1）：33-34. [Liu J J，Xiao L X，Tao Y，Meng L L，Zhu S Y，Lu J T，Shen X Q，Wang Y L，Ma L B，Gu B X. 2018. Integration and introduction of Ridgetail white shrimp （Exopalaemon carinicauda） and multispecies mixed cultivation demonstration[J]. Scientific Fish Farming，（1）：33-34.] doi：10.14184/j.cnki.issn1004-843x.2018.01.022.

盧建平. 2001. 日本沼蝦鰓細(xì)胞的超微結(jié)構(gòu)[J]. 東海海洋，19（3）：19-24. [Lu J P. 2001. Ultramicrostructure of the gill cells of Macrobrachium nipponense[J]. Donghai Marine Science，19（3）：19-24.] doi：10.3969/j.issn.1001-909X. 2001.03.003.

單保黨，胡壯濤，蔡娟，丁玉龍，黎秋云. 2018. 利洋1號和利洋2號南美白對蝦仔蝦耐低氧能力研究[J]. 當(dāng)代水產(chǎn)，43（7）：87-89. [Shan B D，Hu Z T，Cai J，Ding Y L，Li Q Y. 2018. Study on the hypoxia tolerance of Liyang No.1 and No.2 whiteleg shrimp（Litopenaeus vannamei）[J]. Current Fisheries，43（7）：87-89.] doi：10.3969/j.issn.1674-9049.2018.07.027.

王愛敏. 1997. 羅氏沼蝦蚤狀幼體耗氧速率與窒息點測定[J]. 齊魯漁業(yè)，14（2）：37-43. [Wang A M. 1997. Experiment on oxygen consumption rates and asphyxiant points of zoeal larvae of Macrobrachium rosenbergii[J]. Shandong Fisheries，14（2）：37-43.]

王志剛，劉士力，李賢露，呂卓云，毛麗盈，唐超然，鄭榮泉. 2019. 翹嘴紅鲌（Erythroculter ilishaeformis）果糖-1，6-二磷酸醛縮酶（ALDO-C）基因定位、克隆及表達(dá)分析[J]. 海洋與湖沼，50（5）：1138-1145. [Wang Z G，Liu S L，Li X L，Lü Z Y，Mao L Y，Tang C R，Zheng R Q. 2019. Location，cloning and expression analysis of ALDO-C gene in Erythroculter ilishaeformis[J]. Oceanologia et Limnologia Sinica，50（5）：1138-1145.] doi：10.11693/hyhz2019 0300049.

王盼，史文軍，萬夕和，沈輝，沙士兵，黎慧，王李寶，孫瑞健，蔣葛，吳旭干. 2021. 低氧–復(fù)氧對脊尾白蝦呼吸代謝和抗氧化酶活力的影響[J]. 漁業(yè)科學(xué)進(jìn)展，42（4）：106-115. [Wang P，Shi W J，Wan X H，Shen H，Sha S B，Li H，Wang L B，Sun R J，Jiang G，Wu X G. 2021. Effects of hypoxia and reoxygenation on respiratory metabolism enzyme and antioxidant enzyme activities in Exopalaemon carinicauda[J]. Progress in Fishery Sciences，42（4）：106-115.] doi：10.19663/j.issn2095-9869.20200424 002.

楊明，孫盛明，傅洪拓，趙才源，王亞冰，胡宇寧，喬慧，張文宜，龔永生，蔣速飛. 2019. 低氧和復(fù)氧對日本沼蝦抗氧化酶活力及組織結(jié)構(gòu)的影響[J]. 中國水產(chǎn)科學(xué)，26（3）：493-503. [Yang M，Sun S M，F(xiàn)u H T，Zhao C Y，Wang Y B，Hu Y N，Qiao H，Zhang W Y，Gong Y S，Jiang S F. 2019. Hypoxia and reoxygenation on antioxidant enzyme activities and his-tological structure of Macrobrachium nipponense[J]. Journal of Fishery Sciences of China，26（3）：493-503.] doi：10.3724/SP.J.1118.2019.18285.038C597C-F774-40C1-9DC1-F6FE2E5B0876

Chen S F，Zhou Y Q，Chen Y R，Jia G. 2018. fastp：An ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics，34（17）：i884-i890. doi：10.1093/bioinformatics/bty560.

Cheng W，Liu C H，Hsu J P，Chen J C. 2002. Effect of hypo-xia on the immune response of giant freshwater prawn Macrobrachium rosenbergii and its susceptibility to pathogen Enterococcus[J]. Fish & Shellfish Immunology，13（5）：351-365. doi：10.1006/fsim.2001.0411.

Cheng W，Liu C H，Kuo C M. 2003. Effects of dissolved oxygen on hemolymph parameters of freshwater giant prawn，Macrobrachium rosenbergii （de Man）[J]. Aquaculture，220（1）：843-856. doi：10.1016/S0044-8486（02）00534-3.

Crocker C D，Chapman L J，Martínez M L. 2013. Hypoxia-induced plasticity in the metabolic response of a widespread cichlid[J]. Comparative Biochemistry and Phy-siology. Part B：Biochemistry & Molecular Biology，166（2）：141-147. doi：10.1016/j.cbpb.2013.08.002.

Dupont-Prinet A，Pillet M，Chabot D，Hansen T，Tremblay R，Audet C. 2013. Northern shrimp （Pandalus borealis） oxygen consumption and metabolic enzyme activities are severely constrained by hypoxia in the Estuary and Gulf of St. Lawrence[J]. Journal of Experimental Marine Biology and Eco-logy，448：298-307. doi：10.1016/j.jembe.2013.07.019.

Ernst J，Bar-Joseph Z. 2006. STEM：A tool for the analysis of short time series gene expression data[J]. BMC Bioinformatics，7（1）：191. doi：10.1186/1471-2105-7-191.

Grabherr M G，Haas B J，Yassour M，Levin J Z，Thompson D A，Amit I，Adiconis X，F(xiàn)an L，Raychowdhury R，Zeng Q D，Chen Z H，Mauceli E，Hacohen N，Gnirke A，Rhind N，di Palma F，Birren B W，Nusbaum C，Lindblad-Toh K，F(xiàn)riedman N，Regev A. 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nature Biotechnology，29（7）：644-652. doi：10.1038/nbt.1883.

Hu D X，Pan L Q，Zhao Q，Ren Q. 2015. Transcriptomic response to low salinity stress in gills of the Pacific white shrimp，Litopenaeus vannamei[J]. Marine Genomics，24（3）：297-304. doi：10.1016/j.margen.2015.07.003.

Huang W，Ren C H，Li H M，Huo D，Wang Y H，Jiang X，Tian Y S，Luo P，Chen T，Hu C Q. 2017. Transcriptomic analyses on muscle tissues of Litopenaeus vannamei provide the first profile insight into the response to low temperature stress[J]. PLoS One，12（6）：e0178604. doi：10. 1371/journal.pone.0178604.

Iurlaro R，Munoz-Pinedo C. 2016. Cell death induced by endoplasmic reticulum stress[J]. The FEBS Journal，283（14）：2640-2652. doi：10.1111/febs.13598.

Kim S Y，Park J W. 2010. Modulation of hypoxia-inducible factor-1α expression by mitochondrial NADP-dependent isocitrate dehydrogenase[J]. Biochimie，92（12）：1908-1913. doi：10.1016/j.biochi.2010.08.004.038C597C-F774-40C1-9DC1-F6FE2E5B0876

Land S C，Buck L T，Hochachka P W. 1993. Response of protein synthesis to anoxia and recovery in anoxia-tolerant hepatocytes[J]. The American Journal of Physiology，265（1）：R41-R48. doi：10.1152/ajpregu.1993.265.1.R41.

Larkin G，Closs G P，Peake B. 2008. Tolerance and behaviour of the mysid shrimp Tenagomysis novae-zealandiae to low dissolved oxygen[J]. New Zealand Journal of Marine and Freshwater Research，41（3）：317-323. doi：10. 1080/00288330709509919.

Li F H，Luan W，Zhang C S，Zhang J Q，Wang B，Xie Y S，Li S H，Xiang J H. 2009. Cloning of cytoplasmic heat shock protein 90 （FcHSP90） from Fenneropenaeus chinensis and its expression response to heat shock and hypoxia[J]. Cell Stress and Chaperones，14（2）：161-172. doi：10.1007/s12192-008-0069-6.

Li H L，Lin H R，Xia J H. 2017. Differential gene expression profiles and alternative isoform regulations in gill of Nile tilapia in response to acute hypoxia[J]. Marine Biotechnology，19（6）：551-562. doi：10.1007/s10126-017-9774-4.

Li J T，Lü J J，Liu P，Chen P，Wang J J，Li J. 2019. Genome survey and high-resolution backcross genetic linkage map construction of the ridgetail white prawn Exopalaemon carinicauda applications to QTL mapping of growth traits[J]. BMC Genomics，20（1）：598. doi：10.1186/s12864-019-5981-x.

Li Y D，Zhou F L，Huang J H，Yang L S，Jiang S，Yang Q B，He J G，Jiang S G. 2018. Transcriptome reveals involvement of immune defense，oxidative imbalance，and apoptosis in ammonia-stress response of the black tiger shrimp （Penaeus monodon）[J]. Fish & Shellfish Immunology，83：162-170. doi：10.1016/j.fsi.2018.09.026.

Love M I，Huber W，Anders S. 2014. Moderated estimation of fold change and dispersion for RNA-Seq data with DESeq2[J]. Genome Biology，15（12）：550. doi：10.1186/s13059-014-0550-8.

Ma H K，Sun J Q，Xu W Y，Gao W，Hu G W，Lai X F，Yan B L，Gao H. 2020. Cloning and functional study of lipocalin：Retinol-binding protein-like gene family of the ridgetail white prawn，Exopalaemon carinicauda[J]. Molecular Genetics & Genomics，295（2）：453-464. doi：10. 1007/s00438-019-01633-0.

Mortazavi A，Williams B A，McCue K，Schaefferet L，Wold B. 2008. Mapping and quantifying mammalian transcriptomes by RNA-Seq[J]. Nature Methods，5（7）：621-628. doi：10.1038/nmeth.1226.

Nakhoul H，Ke J W，Zhou X，Liao W J，Zeng S X，Hua L. 2014. Ribosomopathies：Mechanisms of disease[J]. Clinical Medicine Insights：Blood Disorders，7：7-16. doi：10. 4137/CMBD.S16952.

Paschke K，Cumillaf J P，Loyola S，Gebauer P，Urbina M，Chimal M E，Pascual C，Rosas C. 2010. Effect of dissolved oxygen level on respiratory metabolism，nutritional phy-siology，and immune condition of southern king crab Lithodes santolla （Molina，1782） （Decapoda，Lithodidae）[J]. Marine Biology，157（1）：7-18. doi：10.1007/s00227-009-1291-1.038C597C-F774-40C1-9DC1-F6FE2E5B0876

Saetan W，Tian C X，Yu J W，Lin X H，He F X，Huang Y，Shi H J，Zhang Y L，Li G L. 2020. Comparative transcriptome analysis of gill tissue in response to hypoxia in silver sillago （Sillago sihama）[J]. Animals，10（4）：628. doi：10.3390/ani10040628.

Shao C N，Zheng H N，Zhou Z X，Li J，Luo X W，Hui G H，Zhao Z D. 2018. Ridgetail white prawn （Exopalaemon carinicauda） K value predicting method by using electronic nose combined with non-linear data analysis model[J]. Food Analytical Methods，11：3121-3129. doi：10. 1007/s12161-018-1297-8.

Shenoy N，Kessel R，Bhagat T D，Bhattacharyya S，Yu Y T，Mcmahon C，Verma A. 2012. Alterations in the riboso-mal machinery in cancer and hematologic disorders[J]. Journal of Hematology & Oncology，5：32. doi：10.1186/1756-8722-5-32.

Smith R W，Houlihan D F，Nilsson G E，Brechin J G. 1996. Tissue-specific changes in protein synthesis rates in vivo during anoxia in crucian carp[J]. The American Journal of Physiology，271（4）：R897-R904. doi：10.1152/ajpregu.1996.271.4.R897.

Soldatov A A，Golovina I V，Kolesnikova E E，Sysoeva I V，Kladchenko E S. 2020. Activity of energy metabolism enzymes and ATP content in the brain and gills of the black sea scorpionfish Scorpaena porcus under short-term hypoxia[J]. Journal of Evolutionary Biochemistry & Physiology，56（3）：224-234. doi：10.1134/S002209302 0030059.

Sun S M，Xuan F J，Ge X P，F(xiàn)u H T，Zhu J，Zhang S Y. 2014. Identification of differentially expressed genes in hepatopancreas of oriental river prawn，Macrobrachium nipponense exposed to environmental hypoxia[J]. Gene，534（2）：298-306. doi：10.1016/j.gene.2013.10.036.

Sun S M，Xuan F J，F(xiàn) J，F(xiàn)u H T，Zhu J，Ge X P，Gu Z M. 2015. Transciptomic and histological analysis of hepatopancreas，muscle and gill tissues of oriental river prawn （Macrobrachium nipponense） in response to chronic hypoxia[J]. BMC Genomics，16（1）：491. doi：10.1186/s12864-015-1701-3.

Ulaje S A，Rojo-Arreola L，Lluch-Cota S E，Ascencio F，Cruz-Hernández P，Sicard M T. 2019. Gene expression and energetic metabolism changes in the whiteleg shrimp （Litopenaeus vannamei） in response to short-term hypoxia[J]. Aquaculture Research，50（4）：994-1004. doi：10.1111/are.13960.

Wang Z L，Zhou J，Li J Y，Lü W，F(xiàn)an L F. 2020. A new insight into the intestine of Pacific white shrimp：Regulation of intestinal homeostasis and regeneration in Litopenaeus vannamei during temperature fluctuation[J]. Comparative Biochemistry and Physiology. Part D：Genomics & Proteomics，35：100687. doi：10.1016/j.cbd.2020.100687.

Xie X，Guo P，Yu H，Wang Y，Chen G. 2018. Ribosomal proteins：Insight into molecular roles and functions in hepatocellular carcinoma[J]. Oncogene，37（3）：277-285. doi：10.1038/onc.2017.343.

（責(zé)任編輯 蘭宗寶）038C597C-F774-40C1-9DC1-F6FE2E5B0876