青海湖裸鯉鹽堿耐受過程中的滲透、免疫、代謝相關基因篩選和分析

張海琛，馬清花，許?？桑⒘樟?，梁 健

(青海大學生態(tài)環(huán)境工程學院，省部共建三江源生態(tài)與高原農(nóng)牧業(yè)國家重點實驗室，青海 西寧 810016)

青海湖是我國面積最大的內(nèi)陸咸水湖泊，湖水pH 為9.1～9.4，26～32 mmol/L，鹽度為11.9～13.49 g/L，是世界上海拔最高的咸水湖之一，其獨特的自然生態(tài)環(huán)境造就了獨特的高寒生態(tài)系統(tǒng)[1-2]。青海湖裸鯉(Gymnocypris przewalskii)俗稱湟魚、無鱗魚等，屬鯉形目 (Cypriniformes)鯉科(Cyprinidae)裸鯉屬 (Gymnocypris)，可以在高海拔的缺氧環(huán)境中生長繁殖，生命力強，是青海湖重要的經(jīng)濟魚類，處于青海湖整個生態(tài)系統(tǒng)的核心地位[3-4]。青海湖裸鯉屬溯河性產(chǎn)卵魚類，每年的4—8 月為其產(chǎn)卵繁殖期，6 月中旬達到其產(chǎn)卵盛期。產(chǎn)卵時從青海湖溯游到淡水支流布哈河、泉吉河等水流緩慢的河灘道，完成產(chǎn)卵后又洄游到青海湖內(nèi)。同時青海湖裸鯉在洄游過程中先是面臨著鹽堿度的下降，在回到青海湖過程中又面臨著鹽堿度的升高[5]，表現(xiàn)出極強的鹽堿適應性和耐受性。我國鹽堿水資源豐富，但絕大部分的水生生物對于水中鹽堿度的變化比較敏感。高鹽度會改變水生生物的滲透壓，高堿度會損傷水生生物DNA[6]、消化酶及免疫相關酶活性下降[7]、耗氧率升高[8]等。同時，鹽度、堿度和pH 對水生生物的生長、代謝、繁殖等方面表現(xiàn)出交互協(xié)同作用[9]。轉(zhuǎn)錄組測序作為研究復雜轉(zhuǎn)錄組變化的強有力工具，可用于物種中尤其是那些基因組未知的物種中的基因發(fā)現(xiàn)和分子信號通路的鑒定，而目前關于青海湖裸鯉在鹽堿度耐受轉(zhuǎn)錄組學方面的研究仍然比較少。本研究中，對青海湖裸鯉在青海湖河口水域(有一定鹽堿度)和淡水下的鰓、腎組織進行轉(zhuǎn)錄組測序，并通過比較分析篩選出在鹽堿耐受過程中發(fā)揮作用的滲透、免疫、代謝相關基因，同時為青海湖裸鯉在鹽堿脅迫下的基因表達水平變化研究奠定基礎。

1 材料與方法

1.1 實驗材料

實驗用青海湖裸鯉取樣于青海湖及其周邊水系，從同一時期 (5—7 月) 的淡水(無脅迫處理，青海省海北藏族自治州剛察縣泉吉河，99.89° E，37.24° N)以及青海湖與哈爾蓋河口 (鹽堿環(huán)境脅迫，pH 9.58，鹽度為6，堿度3.22 mmol/L,100.42° E，37.12° N)處分別收集3 尾青海湖裸鯉(3～4 齡，體長200～300 mm，體重100～150 g)。獲得裸鯉樣品后將魚立即處死，快速分離其鰓組織和腎臟組織，并儲存于液氮中備用。本研究獲得了青海大學實驗動物管理和使用倫理委員會批準，實驗過程 中操作人員嚴格遵守倫理規(guī)范，并按照青海大學倫理委員會制定的規(guī)章制度執(zhí)行。

1.2 RNA 提取和文庫構建

使用TRIzol 法參照說明書分別從2 個采樣地的青海湖裸鯉鰓組織和腎臟組織中提取總RNA，然后用DNase I (TaKaRa，日本)去除基因組DNA。通過使用NanoPhotometer-NP80 (Implen，德國)對RNA 的濃度和純度進行檢測，1%瓊脂糖凝膠電泳用于測量RNA 的完整性。

轉(zhuǎn)錄組文庫的制備及測序由上海美吉生物醫(yī)藥科技有限公司完成。使用Oligo(dT)的磁珠從總RNA 中分離出具有poly(A)的mRNA，進行純化，用fragmentation buffer 將其片段化 (100～400 bp)，并使用隨機引物反轉(zhuǎn)錄成第一鏈cDNA，然后通過添加緩沖液、dNTPs、DNA 聚合酶I 和RNase H 合成雙鏈cDNA。使用AMPure XP beads再次純化雙鏈cDNA。首先對純化的雙鏈cDNA進行末端修復，加A 尾并連接到測序接頭，然后使用AMPure XPbeads 進行片段大小的選擇。最后進行PCR 擴增，并使用AMPure XP beads 純化PCR 產(chǎn)物以獲得最終的cDNA 文庫。構建好文庫后，使用Qubit 2.0 進行初步定量，并將文庫稀釋至1.5 ng/μL。然后使用Agilent 2 100 檢測文庫的插入片段大小 (insert size)。通過實時熒光定量PCR 技術 (qRT-PCR)方法準確定量文庫的有效濃度 (文庫有效濃度>2 nmol/L)以確保文庫質(zhì)量。最后基于 Illumina HiSeq 2000 測序平臺，以Pairedend 雙末端測序技術完成轉(zhuǎn)錄組測序。

1.3 數(shù)據(jù)拼接和功能注釋

由于原始測序數(shù)據(jù)中會包含測序接頭序列、低質(zhì)量讀段、N (N 表示不確定堿基信息)率較高序列及長度過短序列，這將嚴重影響后續(xù)分析的質(zhì)量。因此在分析之前需先對測序所得的原始數(shù)據(jù)進行統(tǒng)計和質(zhì)控，從而得到高質(zhì)量的質(zhì)控數(shù)據(jù)(clean data)以保證后續(xù)分析結(jié)果的準確性，然后將質(zhì)控后得到的有效數(shù)據(jù) (clean reads)使用Trinity 軟件[10]進行無參轉(zhuǎn)錄組從頭組裝。使用Blast X 程序?qū)nigenes 與六大數(shù)據(jù)庫 (NR、Swiss-Prot、Pfam、COG、GO 和KEGG 數(shù)據(jù)庫) 進行比對(E<10-5)，獲得在各數(shù)據(jù)庫的注釋信息，并對其在各數(shù)據(jù)庫的注釋情況進行統(tǒng)計。

1.4 差異表達基因(DEGs)分析

本實驗根據(jù) FPKM (每千個堿基的轉(zhuǎn)錄每百萬映射讀取的片段)方法評估基因的表達量。使用基于負二項分布的DESeq2 軟件對raw counts 進行統(tǒng)計分析，以P-adjust <0.05 &|log2FC| ≥2 作為篩選條件以獲得比較組間差異表達基因(differentially expressed genes，DEGs)，為了控制整體推斷結(jié)果發(fā)生錯誤的概率或頻次，使用BH (Benjamini hochberg)方法對統(tǒng)計檢驗獲得的P-value 進行矯正，然后對樣本間差異表達基因DEGs 進行GO功能注釋和KEGG 通路富集分析。分別將注釋、富集到GO term 和KEGG 通路中的基因進行篩選，獲得滲透、免疫、代謝相關基因。使用Cytoscape 軟件進行網(wǎng)絡的可視化分析，通過計算節(jié)點權重以及節(jié)點的連通度，篩選可能的關鍵基因。

1.5 qRT-PCR 分析

為驗證RNA-Seq 數(shù)據(jù)，使用qRT-PCR 分析隨機選擇的具有不同功能的差異基因的組織表達特異性。使用Primer premier 5 軟件和在線引物設計網(wǎng)站 (https://bioinfo.ut.ee/primer3-0.4.0/) 對所選差異表達基因設計引物 (表1)。使用TRIzol 試劑(Invitrogen，美國) 從青海湖水域和淡水泉吉河的青海湖裸鯉鰓、腎臟組織中提取總RNA。在測量提取的RNA 的完整性和濃度后，使用Prime ScriptTMRT-PCR 試劑盒 (TaKaRa，日本)反轉(zhuǎn)錄為cDNA。以β-actin 為內(nèi)參，使用SYBR?Premix ExTaqTMⅡ (TliRNase H Plus) 試劑盒 (TaKaRa，日本)，在Roche LightCycler 480 (Roche，瑞士) 上進行qRT-PCR 檢測。反應程序為95 °C 30 s,95 °C 5 s,60 °C 30 s (40 個循環(huán))。最后使用2-△△CT[11]方法計算基因的相對表達量，且CT值為各樣品3 個生物學重復的平均值。

表1 實驗所用引物Tab.1 Primers used in the experiments

2 結(jié)果

2.1 測序數(shù)據(jù)組裝及質(zhì)控

本研究分別提取了青海湖水域和淡水泉吉河中的青海湖裸鯉鰓、腎臟組織的總 RNA，獲得并構建了12 個相應的 cDNA 文庫，同時采用 Illumina HiSeq 2000 高通量測序技術對其轉(zhuǎn)錄組進行了測序。青海湖水域獲得的青海湖裸鯉鰓、腎臟組織分別編號為B 和E，淡水泉吉河中獲得的青海湖裸鯉鰓、腎臟組織編號為C、F。各樣品質(zhì)控數(shù)據(jù) (clean data)均達到 6.74 Gb 以上，Q30 堿基百分比在94.38%以上 (表2)。使用Trinity 對所有樣本clean data 進行從頭組裝，將并對組裝結(jié)果進行優(yōu)化評估，結(jié)果顯示，組裝得到的轉(zhuǎn)錄本 (transcript)個數(shù)為880 998，以聚類最長的轉(zhuǎn)錄本為unigene，其個數(shù)為541 429，N50 (組裝獲得的unigene 按長度排列，累加達到總unigene 長度的一半時所對應的unigene 的長度)平均長度為612 bp。其中長度為0～500 bp 的unigenes 共有404 000 個，占比最大，占總unigenes 的75% (圖1)。

圖1 青海湖裸鯉 unigenes 的長度分布Fig.1 Unigene Sequence Length distribution of G. przewalskii

表2 測序數(shù)據(jù)統(tǒng)計表Tab.2 Statistics table of sequencing data

2.2 轉(zhuǎn)錄組注釋

將經(jīng)質(zhì)控后獲得的unigenes 與6 大數(shù)據(jù)庫(NR、Swiss-Prot、Pfam、COG、GO 和KEGG 數(shù)據(jù)庫)進行比對，獲得在各數(shù)據(jù)庫的注釋信息，并對其在各數(shù)據(jù)庫的注釋情況進行統(tǒng)計。共有107 568個unigenes 得到了注釋，注釋率為19.87%。其中有97 816 個unigenes 注釋到了NR 庫，注釋率為18.07%；68 226 個 unigenes 注釋到了 Swiss-Prot 庫，注釋率為12.6%；46 459 個unigenes 注釋到了Pfam 庫，注釋率為8.58%；16 894個unigenes 注釋到了COG 庫，注釋率為3.12%；13 560個unigenes 注釋到了GO 庫，注釋率為2.5%；51 153 個unigenes 注釋到了KEGG 庫，注釋率為9.45% (圖2)。拼接所得的 unigene 與NR 蛋白數(shù)據(jù)庫進行同源性比對，相似序列所占比例最高的幾種物種分別是鯉 (Cyprinus carpio，21.42%)、犀角金線鲃 (Sinocyclocheilus rhinocerous，17.01%)、安水金線鲃 (S.anshuiensis，14.32%)、滇池金線鲃(S.grahami，13.26%)、中華菊頭蝠 (Rhinolophus sinicus，5.39%)、斑馬魚 (Danio rerio，4.71%) (圖3)。

圖2 不同數(shù)據(jù)庫注釋結(jié)果1.NR 數(shù)據(jù)庫,2.Swiss-Prot 數(shù)據(jù)庫,3.Pfam 數(shù)據(jù)庫,4.COG 數(shù)據(jù)庫,5.GO 數(shù)據(jù)庫,6.KEGG 數(shù)據(jù)庫。Fig.2 Annotated results of different databases1.annotated in NR,2.annotated in Swiss-Prot,3.annotated in Pfam,4.annotated in COG,5.annotated in GO,6.annotated in KEGG.

圖3 NR 注釋物種分布Fig.3 NR annotated species distribution of G. przewalskii

2.3 差異表達基因

利用 DEseq 2 軟件分析2 個區(qū)域的青海湖裸鯉鰓、腎臟中差異表達的基因。結(jié)果顯示，共有832 個基因在2 個區(qū)域的組織中共表達，3 269 個基因僅在鰓中 (Bvs.C)表達，其中表達上調(diào)的DEGs 共2 403 個，下調(diào)的1 698 個；4 281 個基因僅在腎臟 (Evs.F)中表達，其中表達上調(diào)的DEGs 共3 506 個，下調(diào)的1 607 個。

2.4 差異表達基因的 GO 功能分類分析

本研究中，GO 功能注釋分析結(jié)果表明共有4 101 (Bvs.C)和5 113 (Evs.F)個差異表達基因分別注釋到了GO 的3 大分類中，即生物過程、細胞組分和分子功能。其中，Bvs.C 和EvsF 中注釋到結(jié)合 (binding)的DEGs 數(shù)量最多，分別為226 個和261 個。其次為細胞過程 (cellular process)，注釋到該GO term 的DEGs 分別為239 個和232 個。注釋到代謝過程 (metabolic process)的DEGs 分別為226 個和194 個，而注釋到單一生物過程 (single-organism process)的DEGs 分別為190個和207 個(圖4)。此外，Bvs.C 和Evs.F 分別有20、13 個DEGs 注釋到了免疫系統(tǒng)過程 (immune system process)GO term。

圖4 差異表達基因的GO 功能注釋分類圖(a) B vs.C;(b) E vs.F 左邊縱坐標表示GO 的二級分類，橫坐標表示包含在該二級分類中的unigene/transcript 的數(shù)目，右邊三種顏色由上到下分別表示GO 的3 大分支 (即生物過程、細胞組分和分子功能)。Fig.4 Cassification of differentially expressed genes in three GO categories in G. przewalskii The ordinate on the left represents the secondary classification of GO,the abscissa represents the number of unigene/transcript included in the secondary classification,and the three colors on the right represent the three major branches of GO (ie BP,CC,MF) from top to bottom.

(圖4 Fig.4)

2.5 差異表達基因的 KEGG 通路分析

2 個區(qū)域的差異表達基因在KEGG 不同的通路中出現(xiàn)不同程度的富集。Bvs.C 中共富集到307 條通路，其中顯著富集的通路主要集中在癌癥 (pathways in cancer)、PI3K-Akt 信號通路 (PI3Kakt signaling pathway)、癌癥中的蛋白多糖 (proteoglycans in cancer)、黏著斑 (focal adhesion)、核糖體 (ribosome)、人T 細胞白血病病毒1 感染(human T-cell leukemia virus 1 infection)、甲型流感(influenza A)、阿米巴原蟲病 (amoebiasis)、癌癥中的MicroRNAs (MicroRNAs in cancer)、細胞因子和細胞因子受體相互作用 (cytokine-cytokine receptor interaction)(圖5-a)。部分基因富集到了吞噬體 (phagosome)、ABC 轉(zhuǎn)運蛋白 (ABC transporters)、IL-17 信號通路 (IL-17 signaling pathway)、補體和凝血級聯(lián) (complement and coagulation cascades)、造血細胞系 (hematopoietic cell lineage)、精氨酸和脯氨酸代謝 (arginine and proline metabolism)、色氨酸代謝 (tryptophan metabolism)、果糖和甘露糖代謝 (fructose and mannose metabolism)等與滲透、免疫和代謝相關的信號通路(表3)。

在產(chǎn)能建設效益、效果上，青海油田明確“產(chǎn)量優(yōu)先、進站優(yōu)先、壓裂優(yōu)先”的產(chǎn)建思路，強調(diào)地面與地下同步規(guī)劃、同步實施，做到新投一口、進站一口。通過有效銜接鉆井作業(yè)、新投施工、地面配套等環(huán)節(jié)，實現(xiàn)投產(chǎn)周期較往年縮短11天，進站周期縮短22天，保障新井產(chǎn)量貢獻。在11月中旬冬供峰值來臨之前，油田天然氣產(chǎn)能建設已完成工作的78%，日增氣超過200萬立方米。

圖5 差異表達基因的KEGG 富集分析Fig.5 KEGG classification of putative functions of differentially expressed genes(a) B vs.C;(b) E vs.F.

表3 鰓組織滲透、免疫、代謝相關信號通路Tab.3 Osmotic,immune,metabolism related pathway in gill

Evs.F 中共富集到303 條通路，顯著富集的通路主要包括PI3K-Akt 信號通路 (PI3K-Akt signaling pathway)、癌癥 (pathways in cancer)、癌癥中的蛋白多糖 (proteoglycans in cancer)、Rap1 信號通路 (Rap1 signaling pathway)、肌動蛋白細胞骨架的調(diào)節(jié) (regulation of actin cytoskeleton)、黏著斑(focal adhesion)、癌癥中的MicroRNAs (MicroRNAs in cancer)、吞噬體 (phagosome)、MAPK 信號通路 (MAPK signaling pathway)、軸突引導 (axon guidance)(圖5-b)。同樣，有部分基因富集到了礦物質(zhì)吸收 (mineral absorption)、近端小管碳酸氫鹽再 生 (proximal tubule bicarbonate reclamation)、cAMP 信號通路 (cAMP signaling pathway)、白細胞跨內(nèi)皮遷移 (leukocyte transendothelial migration)、補體和凝血級聯(lián) (complement and coagulation cascades)、造血細胞系 (hematopoietic cell lineage)、精氨酸合成 (arginine biosynthesis)、泛酸和輔酶A生物合成 (pantothenate and CoA biosynthesis)等與滲透、免疫、代謝相關的信號通路(表4)。

表4 腎組織滲透、免疫、代謝相關信號通路Tab.4 Osmotic,immune,metabolism related pathway in kidney

2.6 鰓組織滲透、免疫、代謝相關基因篩選

根據(jù)差異表達基因的GO 注釋和KEGG 信號通路富集分析，然后使用Cytoscape 軟件，進行網(wǎng)絡的可視化分析，通過計算節(jié)點權重以及節(jié)點的連通度，我們篩選到的青海湖裸鯉滲透相關的基因則主要有絲裂原活化蛋白激酶 (mitogen-activated protein kinase)、鈣/鈣調(diào)蛋白依賴性蛋白激酶(calcium/calmodulin-dependent protein kinase)、瞬時受體電位陽離子通道 (transient receptor potential cation channel subfamily)、囊性纖維化跨膜 (cystic fibrosis transmembrane)、氯通道蛋白(chloride channel protein) 等。免疫相關基因主要包括補體(complement)、整合素 (integrin)、白細胞介素(interleukin)、顆粒酶B (granzyme B)等。代謝相關的基因主要有乙醇脫氫酶 (aldehyde dehydrogenase)、1,25-二羥基維生素D(3)24-羥化酶 (1,25-dihydroxyvitamin D(3) 24-hydroxylase)、3-磷酸肌醇合酶 (inositol-3-phosphate synthase)、細胞色素 P450(cytochrome P450)、果糖二磷酸醛縮酶C (fructosebisphosphate aldolase C)等 (附表1)。

附表1 鰓組織滲透、免疫、代謝相關基因Tab.S1 Osmotic,immune,metabolism related genes in gill

2.7 腎組織滲透、免疫、代謝相關基因篩選

篩選得到的青海湖裸鯉腎臟組織滲透相關蛋白包括溶質(zhì)載體家族 (solute carrier family)、鈣/鈣調(diào)蛋白依賴性蛋白激酶 (calcium/calmodulindependent protein kinase)、絲裂原活化蛋白激酶(mitogen-activated protein kinase)、瞬時受體電位陽離子通道 (transient receptor potential cation channel subfamily)、氯通道蛋白 (chloride channel protein)、鈉/鉀轉(zhuǎn)運ATP 酶 (sodium/potassium-transporting ATPase)等。免疫相關基因主要有整合素 (integrin)、干擾素誘導蛋白VLIG 蛋白亞家族 (interferon-induced very large GTPase)、顆粒酶B (granzyme B)、白細胞介素 (interleukin)、補體 (complement)。代謝相關基因則包括一氧化氮合酶 (nitric oxide synthase)、1,25-二羥基維生素D(3)24-羥化酶 (1,25-dihydroxyvitamin D(3) 24-hydroxylase)、6-磷酸果糖-2-激酶/果糖-2,6-雙磷酸酶3 (6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3)、細胞色素P450 (cytochrome P450)、磷酸烯醇丙酮酸羧激酶 (phosphoenolpyruvate carboxykinase)、磷酸甘油酸變位酶 (phosphoglycerate mutase)等(附表2)。

附表2 腎組織滲透、免疫、代謝相關基因Tab.S2 Osmotic,immune,metabolism related genes in kidney

2.8 qRT-PCR 結(jié)果

為驗證轉(zhuǎn)錄組測序數(shù)據(jù)結(jié)果的可靠性，從上述2 個區(qū)域的實驗數(shù)據(jù)分析結(jié)果中隨機挑選16 個基因進行qRT-PCR，檢測差異表達基因的表達情況，并將該結(jié)果與轉(zhuǎn)錄組測序結(jié)果進行比較 (圖6)。結(jié)果表明，所挑選的基因表達量與轉(zhuǎn)錄組測序結(jié)果趨勢一致，從而驗證了轉(zhuǎn)錄組測序結(jié)果的可靠性。

圖6 差異表達基因的qRT-PCR 檢測RNA-Seq 結(jié)果1.溶質(zhì)載體家族 41 成員 1 基因，2.鈉/鉀轉(zhuǎn)運 ATP 酶亞基 α-1 基因，3.鈣依賴膜結(jié)合蛋白基因，4.酶原顆粒膜蛋白16 基因，5.GTPase IMAP 家族成員 7 基因，6.補體C3 基因，7.NADPH 氧化酶激活劑 1 基因，8.自然殺傷細胞受體2B4 基因，9.一氧化氮合酶基因，10.核因子7 基因，11.TNF 受體相關因子2 基因，12.環(huán)指蛋白 233 基因，13.黏蛋白-17 基因，14.S100 蛋白-A1 基因，15.細胞周期調(diào)控蛋白，16.ras 樣蛋白家族成員10B 基因。Fig.6 Detection of RNA-seq by qRT-PCR of differentially expressed genes1.solute carrier family 41 member 1,2.sodium/potassium-transporting ATPase subunit alpha-1,3.copine-2,4.zymogen granule membrane protein 16,5.GTPase IMAP family member 7,6.complement C3,7.NADPH oxidase activator 1,8.natural killer cell receptor 2B4,9.nitric oxide synthase,10.nuclear factor 7,11.TNF receptor-associated factor 2,12.RING finger protein 223,13.mucin-17,14.protein S100-A1,15.cell cycle control protein 50B,16.ras-like protein family member 10B.

3 討論

水體環(huán)境中鹽堿度含量的變化使得魚類通過改變自身的生理狀態(tài)以調(diào)節(jié)體內(nèi)滲透壓平衡來維持體內(nèi)環(huán)境的穩(wěn)態(tài)[12]。魚類滲透壓調(diào)節(jié)變化的過程中，同時還伴隨著魚體免疫功能的變化以及能量的大量消耗[13]。

鰓是魚類主要的呼吸器官，因具有表面積大和與外界水環(huán)境直接接觸的特點而被認為在水分運輸和離子轉(zhuǎn)運中起主要作用[14]，且對水體環(huán)境的變化高度敏感。此外，鰓組織所分泌的黏液也是魚體抵抗環(huán)境病原微生物入侵和擴散的有力方式。鰓上皮組織中含有大量的上皮組織細胞，內(nèi)含有豐富的Na+-K+-ATP 酶。魚類能通過鰓組織上的泌氯細胞來調(diào)節(jié)體內(nèi)Na+、Cl-以及體液滲透平衡。在高滲環(huán)境下魚類還可以通過腎臟的腎小管重吸收作用，增強對水和離子的重吸收能力[15]，從而維持體內(nèi)環(huán)境的穩(wěn)態(tài)。同時，腎臟作為魚類主要的免疫器官，在鹽堿耐受過程中同樣發(fā)揮著重要的作用。水環(huán)境中鹽堿度的增加，導致魚體抗病變能力下降易受到水體中病原微生物的侵入，在此過程中腎臟上的免疫相關因子的含量就會發(fā)生相應的變化以應對外界環(huán)境的刺激。

通過轉(zhuǎn)錄組測序分析發(fā)現(xiàn)，青海湖裸鯉鹽堿耐受的過程很復雜，從KEGG 通路分析來看，既涉及到癌癥相關的通路 (pathway in cancer)、又有信號轉(zhuǎn)導的過程，如PI3K-Akt 信號通路 (PI3KAkt signaling pathway)、還包括蛋白多糖 (proteoglycans in cancer)的變化等，故我們從滲透、免疫和代謝這3 個方面探討鹽堿耐受過程中發(fā)生作用的相關基因。

3.1 轉(zhuǎn)錄組測序分析

本研究通過無參轉(zhuǎn)錄組高通量測序技術，比較分析了鹽堿水環(huán)境和淡水環(huán)境下青海湖裸鯉鰓、腎臟2 個組織轉(zhuǎn)錄組以及差異基因的表達情況。轉(zhuǎn)錄組結(jié)果分析表明，所構建的12 個cDNA 文庫中，各樣品clean data 均達到 6.74 Gb 以上，且Q30 堿基百分比在94.38%以上。使用Trinity 對所有樣本clean data 進行從頭組裝，共得到541 429個unigene，其中共有107 568 個unigenes 得到了注釋，注釋率為19.87%。本實驗轉(zhuǎn)錄組測序結(jié)果注釋率不高，這可能是由于目前的幾大公共數(shù)據(jù)庫中所記錄到的與青海湖裸鯉相近的物種基因注釋信息較少的原因。NR 數(shù)據(jù)庫比對注釋的結(jié)果顯示，有97 816 個unigenes 注釋到了NR 庫，注釋率為18.07%，其中與鯉 (21.42%)的同源性比對最高，其次為犀角金線鲃 (17.01%)、安水金線鲃(14.32%) 和滇池金線鲃 (13.26%)，這可能是因為青海湖裸鯉與鯉、金線鲃屬魚類同為鯉科魚類，也反映了青海湖裸鯉與金線鲃屬魚類具有很近的親緣關系[3-4]。

3.2 滲透相關基因的篩選及分析

由于魚類生活在水體中，故其在整個生命過程中常常會面臨水環(huán)境滲透壓的改變。維持體液平衡是一個復雜的過程，涉及到了各種滲透調(diào)節(jié)器官，主要包括鰓、消化道、腎臟等[16]。在海水適應 (魚類生活環(huán)境從淡水向海水過渡)和淡水適應 (魚類生活環(huán)境從海水向淡水過渡) 過程中，魚類分別通過吞飲海水、減少尿量、排出離子和停止飲水，增加尿量，降低鰓上皮氯化鈉排出，從水環(huán)境中吸收Na+和Cl-等方式來調(diào)節(jié)滲透壓[17-19]。本研究中腎臟組織中篩選得到的滲透相關基因種類多于鰓組織，表明青海湖裸鯉腎組織在滲透脅迫中也發(fā)揮了重要作用。許多細胞的細胞膜表明都分布有Na+-K+-ATP 酶 (Na+-K+-ATPase,NKA)，這是離子轉(zhuǎn)運過程中起關鍵作用的酶，具有調(diào)節(jié)體內(nèi)Na+、Cl-以及體液滲透平衡的功能。Na+-K+-ATPase 可以通過水解產(chǎn)生能量進行Na+、K+之間的互換，同時也為其他離子轉(zhuǎn)運供能，以此來維持體內(nèi)滲透壓的穩(wěn)定，并為細胞運輸提供能量[12,20]。已有大量研究表明，魚類NKA 對外界環(huán)境滲透壓的變化會隨著其生態(tài)位的改變而有所改變。在對莫桑比克羅非魚(O.mossambicus)[21]、青鳉[22]、尼羅羅非魚[23]等的研究中均發(fā)現(xiàn)，處于高滲環(huán)境下的NKA 活性增加，而低滲環(huán)境下的NKA 的活性則最低。本研究中篩選所得的NKA 的表達同樣也符合這一規(guī)律，腎臟中篩選出的NKA 其表達量同樣為鹽堿脅迫環(huán)境高于淡水環(huán)境。而本研究未在鰓組織中篩選出NKA，這可能是由于鹽堿脅迫下青海湖裸鯉免疫能力的下降從而使得鰓組織NKA表達量過低而未被篩選到，也可能是NKA 在鰓和腎組織中發(fā)揮的作用并不相同，需要進一步研究。有研究表明[24]，胞內(nèi)第二信使cAMP 濃度可隨著鹽度的升高而增加，在激活細胞基底側(cè)Cl-通道的同時還能激活NKA，從外界吸收Na+，從而實現(xiàn)水生動物的滲透調(diào)節(jié)平衡。我們在2 個區(qū)域的腎組織中發(fā)現(xiàn)部分DEGs 富集到了cAMP 信號通路，且差異顯著，提示青海湖裸鯉腎臟也可依賴cAMP來調(diào)節(jié)體內(nèi)滲透平衡。此外，Péqueux 等[25]的研究證明鈣調(diào)蛋白參與離子轉(zhuǎn)運，可激活Cl-通道并對Na+的吸收產(chǎn)生影響。我們也在2 個區(qū)域的鰓、腎組織中均篩選到了作為鈣調(diào)蛋白下游靶基因的鈣/鈣調(diào)蛋白依賴性蛋白激酶的顯著表達。

3.3 免疫相關基因的篩選及分析

鹽堿度是影響水生生物生命活動最主要的外界因素之一，在改變魚體內(nèi)環(huán)境滲透壓穩(wěn)態(tài)的同時還伴隨著魚體免疫功能的變化[13]。腎臟不僅具有泌尿和滲透調(diào)節(jié)的作用，作為硬骨魚主要的免疫器官，在鹽堿耐受過程中發(fā)揮著重要的免疫防御作用。隨著水環(huán)境中鹽堿度的增加，魚體抗病變能力下降，導致易受到水體中病原微生物的侵入，在此過程中腎臟上的免疫相關因子的含量就會發(fā)生相應的變化以應對外界環(huán)境的刺激。

由于青海湖具有鹽堿度含量高、高海拔、紫外線輻射強、貧營養(yǎng)等極端的水文環(huán)境特點，使得生活在其中的青海湖裸鯉在生長較其他魚類相比較為緩慢[26]，并較其他魚表現(xiàn)出極強的鹽堿適應性和耐受性，且高鹽堿環(huán)境脅迫下會導致青海湖裸鯉的免疫能力下降。本研究結(jié)果顯示，免疫相關基因如白細胞介素、補體、免疫球蛋白、干擾素誘導蛋白、顆粒酶B 等的顯著差異表達，說明青海湖裸鯉通過表達這些基因以減緩鹽堿環(huán)境脅迫所導致的應激反應，且大部分免疫相關基因都顯著下調(diào)。這與郭雯翡等[27]的研究結(jié)果相符，其通過構建青海湖裸鯉在高鹽堿環(huán)境下的抑制消減cDNA 文庫，在轉(zhuǎn)錄水平上研究了基因表達差異，且發(fā)現(xiàn)免疫相關的基因表達下調(diào)。Tong 等[28]在提取青海湖裸鯉腎臟和鰓組織樣本進行轉(zhuǎn)錄組測序后，發(fā)現(xiàn)了大量免疫相關候選基因 (如Tolllike 樣受體、干擾素調(diào)節(jié)因子白細胞介素以及腫瘤壞死因子)。同樣，我們也在免疫相關基因中篩選獲得了白細胞介素的顯著差異表達，其他免疫相關基因未在本研究中篩選獲得。但在本研究中，我們篩選到整合素家族相關基因在鰓、腎臟組織中均顯著差異表達，表明除白細胞介素外，該基因也是青海湖裸鯉免疫防御系統(tǒng)中的重要組成因子。此外，本研究中還發(fā)現(xiàn)，青海湖水域中大部分免疫相關基因的下調(diào)可能是由于青海湖水域的青海湖裸鯉所面對的微生物種類較淡水少，故使得其免疫應激能力較淡水下調(diào)。

補體作為一種存在于脊椎動物血清和組織液中的血清蛋白質(zhì)，可在機體被病原體侵染時激活，通過吞噬或裂解病原微生物來介導免疫應答和炎癥反應。已有對哺乳動物的研究表明，補體和凝血系統(tǒng)涉及了炎癥的發(fā)病機制[29]。本研究中，鰓組織和腎臟中顯著富集到補體和凝血級聯(lián) (complement and coagulation cascades)信號通路的DEGs 也主要為補體 (complement)和補體成分(complement component)。腎臟組織中補體基因的表達顯著上調(diào)表明了在免疫防御和清除病原體功能方面，補體發(fā)揮了十分重要的作用，這與Encinas 等[30]在對斑馬魚不同組織器官進行VHSV 病毒感染后發(fā)現(xiàn)補體基因表達上調(diào)的結(jié)果一致。C3作為補體系統(tǒng)激活通路的中間樞紐，可通過激活末端補體通路來發(fā)揮清除病原體功能，其主要在肝臟中表達。且有研究表明魚類C3 表達具有組織差異性，如虹鱒 肝臟中C3 表達量最高，肝外組織表達量則較低[31]。在對牙鲆 (Paralichthys olivaceus)進行鰻弧菌感染前后的脾臟轉(zhuǎn)錄組研究中未檢測到C3 在脾臟中表達[32]。在對日本沼蝦(Macrobrachium nipponense)進行氨氮脅迫后的轉(zhuǎn)錄組分析中也發(fā)現(xiàn)肝胰腺、鰓和心臟組織中C3 基因的表達量相對肌肉和血液組織表達量高[33]，本研究也僅在鰓組織中篩選獲得C3 基因顯著表達。此外，鰓組織中補體基因的下調(diào)也表明不同組織在補體防御和清除病原體方面可能具有不同的響應機制。并且本研究2 個區(qū)域的青海湖裸鯉鰓組織中篩選獲得的免疫相關基因的數(shù)量多于腎臟，推測在裸鯉鹽堿耐受過程的免疫防御能力方面，鰓組織的免疫調(diào)節(jié)能力可能更強。

3.4 代謝相關基因的篩選及分析

外界水環(huán)境的變化通常會導致魚類能量的大量消耗，不同滲透壓環(huán)境下的魚類為適應滲透壓會通過改變自身的生理變化來調(diào)整機體的滲透壓平衡，使其能夠產(chǎn)生足夠的耗氧量以維持自身正常生命活動[34-35]。本研究中，青海湖裸鯉在從淡水組到河口組時生態(tài)位發(fā)生了一定的改變，導致其可能選擇不同的代謝策略以適應環(huán)境中鹽堿度的增加。

Su 等[36]在對3 種脅迫下雜交羅非魚的鰓轉(zhuǎn)錄組研究中發(fā)現(xiàn)，鹽堿脅迫下富集到的通路主要為生物合成和代謝相關，同樣本研究中篩選得到的代謝相關信號通路也主要與生物合成和代謝相關，主要集中在氨基酸代謝相關通路。Jiang 等[37]研究也證實了氨基酸代謝在滲透調(diào)節(jié)中起著關鍵作用。精氨酸是一種在有機體內(nèi)具有重要生理、代謝調(diào)控功能的半必須氨基酸，同時其作為內(nèi)源性一氧化氮合成的唯一前體，可在一氧化氮合成酶的作用下生成一氧化氮，生成的一氧化氮進而可以調(diào)節(jié)機體內(nèi)產(chǎn)能底物的代謝過程[38]。精氨酸激酶可催化生成ADP 和磷酸精氨酸，為細胞提供能量需求，同時生成的精氨酸則作為一氧化氮合成酶催化生成一氧化氮的唯一反應底物，結(jié)合生成一氧化氮從而其能夠發(fā)揮生理功能。本研究中，青海湖水域組和淡水組鰓組織在精氨酸和脯氨酸代謝途徑中表達上調(diào)的重要基因為一氧化氮合成酶基因，而在2 個區(qū)域的腎臟組織中該基因則顯著下調(diào)表達，說明在不同組織中一氧化氮合成酶可能是通過不同的調(diào)節(jié)方式參與機體的能量代謝過程。細胞色素P450 (cytochromeP450,CYP450)可參與生物活性分子如類固醇、維生素和脂肪酸的代謝。維生素D 的生理功能主要通過活性形式的1,25-(OH)D3來實現(xiàn)的[39]。細胞色素P450 將與維生素D3結(jié)合蛋白結(jié)合后的維生素D3 羥化為2,5-(OH)D3，轉(zhuǎn)運到腎臟后生成具有活性的1,25-(OH)D3。而腎臟分泌的1,25-二羥基維生素 D(3) 24-羥化酶可以在2,5-(OH)D3或1,25-(OH)D3濃度過高時通過與之相互作用來抑制其表達，在實現(xiàn)體內(nèi)維生素D 含量的動態(tài)平衡的同時發(fā)揮調(diào)節(jié)機體鈣磷代謝以及部分免疫調(diào)節(jié)功能[40]。我們在不同區(qū)域的鰓、腎臟組織中均篩選得到了1,25-二羥基維生素 D(3)24-羥化酶和細胞色素P450 基因的顯著表達，表明鹽堿耐受過程中青海湖裸鯉可通過細胞色素P450 和1,25-二羥基維生素 D(3) 24-羥化酶的上調(diào)來穩(wěn)定體內(nèi)的鈣磷代謝過程。

（作者聲明本文無實際或潛在的利益沖突）