陳 郁， 龔 亮， 羅勇軍

(陸軍軍醫(yī)大學(xué)陸軍衛(wèi)勤訓(xùn)練基地軍事醫(yī)學(xué)地理學(xué)教研室， 重慶 400038)

當(dāng)人體進(jìn)入海拔2 500 m及其以上區(qū)域時，由于環(huán)境變化尤其是氧分壓的下降，會對機(jī)體造成明顯的生理學(xué)效應(yīng)，因而稱之為醫(yī)學(xué)高原[1]。在移居高原的部分人群中，由于個體對高原環(huán)境刺激反應(yīng)強(qiáng)烈，將發(fā)生包括急性高原反應(yīng)、高原肺水腫(high altitude pulmonary edema, HAPE)、高原腦水腫(high altitude cerebral edema, HACE)等急性高原病[2]。目前，全世界有1.4億人生活在海拔3 000 m以上地區(qū)[3]，而2 500 m以上地區(qū)生活的人群更多。近年來,到高原的平原人群顯著增多，如何維護(hù)其生命健康，保證其在高原正常生活是高原醫(yī)學(xué)研究的熱點(diǎn)和前沿。伴隨著高通量技術(shù)尤其是大規(guī)模測序和芯片的應(yīng)用，對人體在面臨高原環(huán)境時的基因表達(dá)調(diào)控和作用網(wǎng)絡(luò)的研究也越來越深入。其中，長鏈非編碼RNA(long non-coding RNA, lncRNA)及其調(diào)控網(wǎng)絡(luò)在低氧反應(yīng)中的作用得到了越來越多的關(guān)注。

lncRNA是一類長度大于200個核苷酸的非編碼RNA。研究發(fā)現(xiàn)，部分lncRNA具有編碼肽類的能力[4]。隨著研究的不斷深入，發(fā)現(xiàn)lncRNA在轉(zhuǎn)錄、轉(zhuǎn)錄后修飾，以及RNA翻譯等多個途徑參與基因的表達(dá)調(diào)控，進(jìn)而影響細(xì)胞的代謝、增殖和凋亡等[5]。按照來源和功能定位的不同，lncRNA也分為核基因來源和線粒體來源，兩者也在發(fā)揮著不同的生物學(xué)效應(yīng)[6]。研究發(fā)現(xiàn)，mtDNA3010位點(diǎn)的A基因型是漢族人群HAPE和急性高原反應(yīng)的保護(hù)因素，而G基因型是危險因素[7, 8]。同時也發(fā)現(xiàn)，該位點(diǎn)的A基因型是高原移居漢族高原紅細(xì)胞增多癥的危險因素，而G基因型是保護(hù)因素[9]。mtDNA3010位點(diǎn)位于編碼16s rRNA，難以從常規(guī)分子生物學(xué)途徑解釋其在缺氧過程中因變異引起的機(jī)制變化。而該位點(diǎn)A/G基因型在急性高原病和慢性高原病上相互矛盾的結(jié)果提示，mtDNA變異在缺氧反應(yīng)過程中的作用更為復(fù)雜，需要從多維度探討其可能的作用。

為進(jìn)一步探討mtDNA3010位點(diǎn)A/G基因型變異在缺氧反應(yīng)中的作用，本研究在構(gòu)建mtDNA3010位點(diǎn)A/G兩種基因型融合細(xì)胞的基礎(chǔ)上，采用表達(dá)譜芯片和qRT-PCR驗(yàn)證急性缺氧處理后lncRNA和mRNA表達(dá)譜的變化。通過生物信息學(xué)方法構(gòu)建共表達(dá)網(wǎng)絡(luò)，預(yù)測其潛在功能，為分析mtDNA變異在低氧反應(yīng)中的作用提供理論依據(jù)。

1 材料與方法

1.1 細(xì)胞來源與處理

ρ0206細(xì)胞是由骨肉瘤細(xì)胞143BTK-，經(jīng)溴化乙錠(EB)處理而形成一種特殊的無線粒體細(xì)胞，由浙江大學(xué)管敏鑫教授惠贈，作為供體用于構(gòu)建mtDNA3010A和mtDNA3010G兩種基因型融合細(xì)胞。

1.2 主要試劑

實(shí)驗(yàn)所需主要試劑如下：細(xì)胞總RNA提取試劑Trizol購于美國Invitrogen公司，RNA逆轉(zhuǎn)錄試劑盒、SYBR Green購于Takara公司；DMEM高糖培養(yǎng)基購于美國康寧公司，胎牛血清購于德國PAN公司；胰酶購于武漢博士德生物技術(shù)有限公司。

1.3 融合細(xì)胞缺氧處理

以ρ0206細(xì)胞為基礎(chǔ)，本課題組前期成功構(gòu)建了mtDNA3010A和mtDNA3010G兩種基因型融合細(xì)胞，具體見參考文獻(xiàn)[10]。細(xì)胞在37℃，5%CO2和常氧條件下培養(yǎng)，取對數(shù)期生長的兩種融合細(xì)胞，計(jì)數(shù)后取20×104的細(xì)胞接種于培養(yǎng)皿中，經(jīng)1%O2缺氧培養(yǎng)24 h后，建立細(xì)胞急性缺氧模型。缺氧結(jié)束后，棄培養(yǎng)上清，PBS清洗后加入Trizol，反復(fù)輕柔吹打以獲取所有細(xì)胞。每種融合細(xì)胞各送樣3份。

1.4 lncRNA表達(dá)譜和mRNA表達(dá)譜的篩選

采用北京博奧晶典生物技術(shù)有限公司提供的Human LncRNA Array v4(4x180K)芯片，包括了40916個lncRNA探針和34235個mRNA探針，涵蓋了人體現(xiàn)有的絕大部分lncRNA和mRNA?？俁NA的提取，標(biāo)本質(zhì)控以及后續(xù)芯片的處理、圖像采集、數(shù)據(jù)讀取在北京博奧晶典生物技術(shù)有限公司完成。

1.5 差異mRNA驗(yàn)證和功能預(yù)測

對兩種融合細(xì)胞經(jīng)缺氧處理后，采取Trizol法提取總RNA。根據(jù)mRNA芯片結(jié)果，以β-actin為內(nèi)參，使用qRT-PCR法驗(yàn)證差異表達(dá)的mRNA，每組實(shí)驗(yàn)重復(fù)3次。采用Gene Ontology數(shù)據(jù)庫和京都基因與基因組大百科全書(Kyoto Encyclopedia of Genes and Genomes, KEGG)數(shù)據(jù)庫對驗(yàn)證的差異基因進(jìn)行功能注釋，并進(jìn)行一級、二級信號通路和功能預(yù)測。

1.6 lncRNA-mRNA共表達(dá)網(wǎng)絡(luò)構(gòu)建和靶基因預(yù)測

根據(jù)上述GO和KEGG的預(yù)測結(jié)果，包括其可能涉及到的基因、蛋白之間的相互作用關(guān)系，采用Cytoscape軟件繪制lncRNA-mRNA共表達(dá)網(wǎng)絡(luò)圖，分析lncRNA和mRNA之間的信號轉(zhuǎn)導(dǎo)網(wǎng)絡(luò)。根據(jù)差異lncRNA預(yù)測其可能作用的mRNA，并通過GO和KEGG數(shù)據(jù)庫預(yù)測靶基因的相關(guān)功能變化。

1.7 統(tǒng)計(jì)學(xué)處理

芯片數(shù)據(jù)分析由北京博奧晶典生物技術(shù)有限公司完成，其他數(shù)據(jù)分析采用SPSS20.0完成。芯片雜交掃描后的tiff格式圖片數(shù)據(jù)采用Feature Extraction提數(shù)軟件進(jìn)行預(yù)處理分析，然后采用GeneSpring GX軟件計(jì)算基因表達(dá)差異和統(tǒng)計(jì)學(xué)差異顯著性P值。數(shù)據(jù)處理過程主要包括：芯片掃描tiff圖用Feature Extraction軟件處理得到原始數(shù)據(jù)文件 (.txt)。將原始數(shù)據(jù)文件(.txt)導(dǎo)入GeneSpring軟件，寫入分組等參數(shù)信息。進(jìn)行每個樣品的數(shù)據(jù)歸一化和QC分析。采用Cluster3.0軟件進(jìn)行Cluster分析和圖形化展示。根據(jù)分組信息，進(jìn)行差異比較，得到差異基因。所得數(shù)據(jù)結(jié)果以mtDNA3010G基因型融合細(xì)胞為對照，以|FC|>1.2、P<0.05、FRD<0.05為標(biāo)準(zhǔn)，篩選差異表達(dá)的lncRNA和mRNA。在驗(yàn)證芯片結(jié)果時，按照2-△△Ct法計(jì)算分析mtDNA3010A和mtDNA3010G基因型融合細(xì)胞缺氧后基因表達(dá)差異。

2 結(jié)果

2.1 兩種融合細(xì)胞缺氧后的lncRNA和mRNA表達(dá)譜差異

經(jīng)1%O2缺氧培養(yǎng)24 h后，根據(jù)芯片檢測結(jié)果，對不同樣本中差異表達(dá)的lncRNA和mRNA進(jìn)行聚類分析，系統(tǒng)熱圖分析結(jié)果見圖1。其中，每一個縱列代表1個樣本，AR和GR分別代表mtDNA3010A/G基因型融合細(xì)胞。紅色代表表達(dá)上調(diào)，綠色代表表達(dá)下調(diào)。和mtDNA3010G基因型融合細(xì)胞相比：mtDNA3010A基因型融合細(xì)胞表達(dá)上調(diào)的lncRNA共有668個，上調(diào)倍數(shù)最大的為ENST00000553679.1倍數(shù)=2.89)，表達(dá)下調(diào)的lncRNA共有1098個，下調(diào)倍數(shù)最大的為HIT000092395_03(倍數(shù)=2.80)(圖2)；表達(dá)上調(diào)的mRNA共有1151個，上調(diào)倍數(shù)最大的為成纖維細(xì)胞生長因子1(fibroblast growth factor 1, FGF1) (倍數(shù)=3.60)，表達(dá)下調(diào)的mRNA共有539個，下調(diào)倍數(shù)最大的為G蛋白耦聯(lián)受體151(G protein-coupled receptor 151, GPR151 )(倍數(shù)=2.07)(圖2，表1，表2)。上述結(jié)果表明，即便線粒體僅有一個位點(diǎn)變異，在急性缺氧時也可能誘導(dǎo)顯著的基因變化并被基因芯片檢測到。

Fig. 1 Heat map of differential expressions of lncRNA and mRNA in mtDNA3010A/G fusion cells after acute hypoxia treatment

Fig. 2 Scatter plots of differential expressions of lncRNA and mRNA in mtDNA3010A/G fusion cells after acute hypoxia treatment

橫坐標(biāo)、縱坐標(biāo)分別為mtDNA3010G基因型融合細(xì)胞和mtDNA3010A基因型融合細(xì)胞表達(dá)的lncRNA(圖2左)、mRNA(圖2右)經(jīng)log2轉(zhuǎn)化后的值。紅色直線為基因上調(diào)表達(dá)閾值線，綠色直線為基因下調(diào)表達(dá)閾值線，中間灰色直線為總體表達(dá)量的擬合線，方程為該擬合線的方程，R為mtDNA3010A基因型和mtDNA3010G基因型的兩組樣本的相關(guān)系數(shù)。

Tab. 1 Top 10 lncRNAs showed significant differences in expression changes

Tab. 2 Top 10 lncRNAs showed significant differences in expression changes

2.2 差異mRNA驗(yàn)證和功能預(yù)測

根據(jù)mRNA芯片結(jié)果，篩選了部分基因進(jìn)行驗(yàn)證，結(jié)果見表3。以FGF1為例，芯片結(jié)果為正 3.60，提示經(jīng)急性缺氧處理后，與mtDNA3010A融合細(xì)胞相比，mtDNA3010G融合細(xì)胞的該基因表達(dá)升高3.60倍。然而經(jīng)qRT-PCR驗(yàn)證后，發(fā)現(xiàn)該差異表達(dá)基因?yàn)榧訇栃?。上述結(jié)果提示，表達(dá)譜芯片檢測出的差異表達(dá)基因有許多假陽性的結(jié)果，必須通過qRT-PCR對結(jié)果驗(yàn)證。唯有驗(yàn)證成功的差異表達(dá)基因方能完成后續(xù)的生物信息學(xué)分析。

Tab. 3 Results of chips and verified by qRT-PCR

通過GO聚類分析，和mtDNA3010G基因型融合細(xì)胞相比，mtDNA3010A基因型融合細(xì)胞差異表達(dá)基因涉及的生物學(xué)過程包括纖毛合成、上皮纖毛運(yùn)動、細(xì)胞微管運(yùn)動、纖毛組裝、糖脂類生物合成以及心血管系統(tǒng)發(fā)生；涉及的細(xì)胞組分包括DNA復(fù)制復(fù)合物、核小體、微管等細(xì)胞骨架、細(xì)胞器、蛋白-DNA復(fù)合體等；差異表達(dá)基因涉及的分子功能包括芳香胺乙酰轉(zhuǎn)移酶活性，蛋白二聚體活性，磷脂酶A2活化因子活性，醌類輔助因子甲基轉(zhuǎn)移酶活性，脂蛋白轉(zhuǎn)運(yùn)體活性以及C-甲基轉(zhuǎn)移酶活性等。

KEGG通路分析提示，差異表達(dá)的基因涉及的信號通路包括到酒精中毒，系統(tǒng)性紅斑狼瘡，咖啡因代謝，化學(xué)物致癌，腫瘤轉(zhuǎn)錄失調(diào)節(jié)以及p53信號通路等(圖3)。

Fig. 3 KEGG-PATHWAY enrichment analysis of differentially expressed genesin mtDNA3010A/G fusion cells after acute hypoxia treatment

2.3 lncRNA-mRNA共表達(dá)網(wǎng)絡(luò)的構(gòu)建

按照相關(guān)性系數(shù)大于0.90，并且P<0.05的標(biāo)準(zhǔn)，篩選出顯著表達(dá)的lncRNA-mRNA前1000個共表達(dá)基因?qū)Γ捎肅ytoscape軟件生成并分析新構(gòu)建的lncRNA-mRNA共表達(dá)網(wǎng)絡(luò)。正相關(guān)表明隨著該lncRNA的表達(dá)增加，mRNA表達(dá)水平亦顯著上升，負(fù)相關(guān)即顯著下降，相關(guān)系數(shù)越大表明變化的幅度越顯著。本研究中發(fā)現(xiàn)，在急性缺氧處理后，與mtDNA3010G基因型融合細(xì)胞相比，mtDNA3010A基因型融合細(xì)胞差異表達(dá)lncRNA關(guān)聯(lián)程度最高的基因包括，錨蛋白3(Ankyrin 3, ANK3)，B細(xì)胞淋巴瘤樣因子13(BCL2 Like 13, BCL2L13)和發(fā)動蛋白結(jié)合蛋白(dynamin binding protein, DNMBP)等，其中最為顯著的是ANK3基因，預(yù)測有1695個lncRNA可能與其有相互作用，提示上述基因也可能參與了mtDNA3010A基因型對急性高原反應(yīng)和HAPE的保護(hù)作用。

2.4 lncRNA的靶基因預(yù)測

在構(gòu)建lncRNA-mRNA共表達(dá)網(wǎng)絡(luò)圖的基礎(chǔ)上，根據(jù)lncRNA在染色體上所屬區(qū)域，預(yù)測其在10kb范圍內(nèi)可能有作用的mRNA(即順式預(yù)測)，同時利用blast工具比對lncRNA序列和mRNA的3`-UTR區(qū)域(即反式預(yù)測)，并綜合不同lncRNA預(yù)測的同一靶基因數(shù)目，共篩選了10個與差異lncRNA存在高度相互作用的靶基因，具體見表4。生物信息學(xué)分析提示mtDNA3010位點(diǎn)基因型變異后，缺氧能夠誘導(dǎo)細(xì)胞出現(xiàn)一系列的動態(tài)調(diào)控網(wǎng)絡(luò)變化，有可能參與了HAPE的發(fā)生過程。

Tab. 4 Predicted target genes by differently expressed lncRNAs

2.5 預(yù)測靶基因的GO和KEGG分析

參照2.2的方法，對lncRNA預(yù)測的10個靶基因進(jìn)行GO和KEGG分析。結(jié)果提示，mtDNA3010G基因型在缺氧后可能會影響細(xì)胞在能量代謝和炎癥反應(yīng)相關(guān)基因的表達(dá)，能量代謝異常和炎癥反應(yīng)過度可能是HAPE發(fā)生中的重要環(huán)節(jié)。GO結(jié)果表明，和mtDNA3010G基因型融合細(xì)胞相比，mtDNA3010A基因型融合細(xì)胞經(jīng)急性缺氧處理后，差異lncRNA的靶基因涉及的生物學(xué)過程包括糖原分解代謝、細(xì)胞對白介素6(interleukin 6, IL-6)的反應(yīng)、內(nèi)皮細(xì)胞活化、碳水化合物的合成代謝等。在涉及到的細(xì)胞成分方面，主要包括糖原顆粒、蛋白磷酸酶1型復(fù)合體、Bcl3-Bcl10復(fù)合物、Bcl3-NFκB復(fù)合物和IL-6受體復(fù)合物等。在分子功能方面，主要為與IL-11結(jié)合、NADP轉(zhuǎn)氫酶活性、黃嘌呤氧化酶活性、IL-11受體活性、氧化還原酶活性等密切相關(guān)。KEGG結(jié)果提示，靶基因的功能富集于煙酰胺代謝、酒精代謝、咖啡因代謝、細(xì)胞連接、腎素-血管緊張素系統(tǒng)、胰島素抵抗和過氧化物酶活性等。

3 討論

人體在經(jīng)受高原環(huán)境刺激時，機(jī)體會出現(xiàn)多方面代償變化。良好的代償有助于機(jī)體習(xí)服高原環(huán)境，保持健康狀態(tài)，而失代償則會導(dǎo)致各種高原病的發(fā)生。大量的證據(jù)提示，在不同的人群中HAPE的發(fā)生與線粒體變異有關(guān)，表明線粒體變異及其誘導(dǎo)的功能變化是HAPE重要的發(fā)病環(huán)節(jié)[11]。本研究表明，mtDNA3010位點(diǎn)由A基因型突變?yōu)镚基因型后，在缺氧作用下，會導(dǎo)致lncRNA-mRNA調(diào)控網(wǎng)絡(luò)的顯著變化，提示lncRNA和mRNA之間存在著復(fù)雜而又緊密的相互作用。本研究還提示，線粒體與核基因之間有聯(lián)系，不僅僅是體現(xiàn)在mRNA表達(dá)差異，lncRNA同樣介導(dǎo)了線粒體與細(xì)胞核之間的相互作用，進(jìn)而影響細(xì)胞的代謝過程和機(jī)體功能的正常發(fā)揮。上述結(jié)果也表明，mtDNA3010位點(diǎn)由G基因型突變?yōu)锳基因型時，在缺氧條件下會引起細(xì)胞基因表達(dá)改變，然而每個差異表達(dá)基因的功能并不能完全解釋細(xì)胞在缺氧刺激時的功能變化，需要相互作用形成級聯(lián)放大效應(yīng)，最終導(dǎo)致疾病的發(fā)生。

lncRNA參與了動脈粥樣硬化、肺動脈高壓，尤其是缺氧性肺動脈高壓的發(fā)生，而缺氧性肺動脈高壓是HAPE發(fā)生的中心環(huán)節(jié)。研究發(fā)現(xiàn)，lncRNA LINC-PINT和LINC00599的多態(tài)性與漢族HAPE的發(fā)生風(fēng)險有關(guān)[12]。LINC-PINTrs157928位點(diǎn)的CT基因型是中國漢族人群HAPE的保護(hù)因素，LINC00599rs2272026位點(diǎn)的CC基因型是32歲以下人群HAPE的保護(hù)因素，而在32歲以上則為危險因素。目前，lncRNA在HAPE發(fā)生中的作用研究報(bào)道較少，尚未涉及具體的分子機(jī)制。本次研究中雖然也發(fā)現(xiàn)了大量差異表達(dá)的lncRNA，其中ENST00000553679.1和ENST00000453008.2差異最為明顯，但是有關(guān)這兩種lncRNA的具體功能的研究尚鮮有報(bào)告，其機(jī)制還需要進(jìn)一步探討。

在本研究中，采用qRT-PCR驗(yàn)證芯片結(jié)果后發(fā)現(xiàn)，A基因型的ATP結(jié)合盒亞家族C成員6(ATP binding cassette subfamily C member 6, ABCC6)、細(xì)胞色素P450家族1亞家族B成員1(cytochrome P450 family 1 subfamily B member 1, CYP1B1)和細(xì)胞色素C氧化酶Ⅰ(cytochrome C oxidase I, COX1)表達(dá)顯著上調(diào)，而脂氧合酶同源結(jié)構(gòu)域1(lipoxygenase homology domains 1, LOXHD1)表達(dá)顯著下調(diào)。ABCC6編碼的是一類跨膜蛋白，能夠轉(zhuǎn)運(yùn)細(xì)胞內(nèi)外多種物質(zhì)，并同時具備ATP酶活性，其突變與皮膚彈性假黃瘤有關(guān)。CYP1B1編碼重要的代謝酶類，參與了包括脂肪酸，類固醇激素和維生素的代謝，是重要的炎癥調(diào)控基因[13]。COX1是線粒體氧化呼吸鏈復(fù)合體Ⅳ的重要組成部分，在形成電化學(xué)梯度驅(qū)動中ATP合成發(fā)揮著關(guān)鍵作用。LOXHD1編碼一種高度保守的蛋白，主要表達(dá)在內(nèi)耳處的毛細(xì)胞，其突變與聽力缺失有關(guān)。目前有關(guān)上述基因在缺氧反應(yīng)及高原病發(fā)生中作用的研究報(bào)道。結(jié)合相關(guān)基因的功能預(yù)測，表達(dá)上調(diào)的ABCC6、CYP1B1和COX1可能會改善細(xì)胞的能量代謝，抑制炎癥反應(yīng)，促進(jìn)能量生成，有助于機(jī)體應(yīng)對高原環(huán)境的習(xí)服。表達(dá)下調(diào)的LOXHD1作用尚不清楚，根據(jù)其誘發(fā)聽力障礙的同類型基因潛在功能提示，可能會影響鈣泵的功能。鈣是收縮血管的重要因素，而鈣拮抗劑是降低血壓的主要藥物，而降低肺動脈壓是治療HAPE的重要方法。上述結(jié)果提示，本研究中發(fā)現(xiàn)的差異表達(dá)基因可能通過多種途徑參與了HAPE的發(fā)生過程。

在多個不同人群的HAPE遺傳易感性研究中，均發(fā)現(xiàn)HAPE的發(fā)生與mtDNA變異有關(guān)，進(jìn)一步分析提示，這些變異在缺氧條件下誘發(fā)的能量代謝異常有可能是導(dǎo)致HAPE發(fā)生的重要機(jī)制。氧化磷酸化是細(xì)胞的主要供能方式，在缺氧條件下這一方式受到明顯的抑制，難以提供足夠的能量供細(xì)胞代謝，導(dǎo)致細(xì)胞的某些代謝過程發(fā)生異常，在缺氧刺激下不能發(fā)揮有效的應(yīng)對，導(dǎo)致細(xì)胞功能失常，進(jìn)而通過一系列級聯(lián)反應(yīng)導(dǎo)致了HAPE的發(fā)生，而mtDNA的變異可能加劇這一病理反應(yīng)。此外，據(jù)報(bào)道，HAPE患者血清中的IL-6和腫瘤壞死因子(tumor necrosis factor-α, TNF-α)等炎癥因子含量明顯升高[14]，提示炎癥反應(yīng)與HAPE的發(fā)生密切相關(guān)[15]。在本實(shí)驗(yàn)中，GO分析和KEGG富集均提示，差異lncRNA預(yù)測的靶基因，以及差異表達(dá)基因參與了炎癥反應(yīng)和細(xì)胞能量代謝，表明mtDNA3010位點(diǎn)由G基因型突變?yōu)锳基因型時，也可導(dǎo)致與能量代謝和炎癥反應(yīng)相關(guān)因子的變化，為闡述HAPE的發(fā)生提供了新的理論依據(jù)。

綜上所述，在急性缺氧條件下，mtDNA3010位點(diǎn)由A基因型突變?yōu)镚基因型后，會導(dǎo)致相關(guān)lncRNA和mRNA的表達(dá)差異，并在生物學(xué)過程、細(xì)胞成分和分子功能，以及代謝途徑等方面影響了細(xì)胞的炎癥反應(yīng)、能量代謝等核心進(jìn)程，可能在誘導(dǎo)急性高原反應(yīng)和HAPE的發(fā)生、發(fā)展中發(fā)揮了重要作用。上述差異lncRNA和mRNA有望成為HAPE的易感標(biāo)志物，并為探索線粒體介導(dǎo)的細(xì)胞調(diào)控網(wǎng)絡(luò)變化在高原病發(fā)生中的機(jī)制提供了重要的研究方向，也為從線粒體角度分析能量代謝和炎癥反應(yīng)在缺氧反應(yīng)中的機(jī)制提供有益資料。