Bioscience Reports
Integrative analysis of miRNA–mRNA network in high altitude retinopathy by bioinformatics analysis
DOI 10.1042/BSR20200776 , Volume: 41 , Issue: 1
Article Type: research-article, Article History

High-altitude retinopathy (HAR) is an ocular manifestation of acute oxygen deficiency at high altitudes. Although the pathophysiology of HAR has been revealed by many studies in recent years, the molecular mechanism is not yet clear. Our study aimed to systematically identify the genes and microRNA (miRNA) and explore the potential biomarkers associated with HAR by integrated bioinformatics analysis. The mRNA and miRNA expression profiles were obtained from the Gene Expression Omnibus database. We performed Gene Ontology functional annotations and Kyoto Encyclopedia of Genes and Genomes pathway analysis. Potential target gene analysis and miRNA–mRNA network analysis were also conducted. Quantitative RT-PCR (qRT-PCR) was used to validate the results of the bioinformatics analysis. Through a series of bioinformatics analyses and experiments, we selected 16 differentially expressed miRNAs (DE-miRNAs) and 157 differentially expressed genes related to acute mountain sickness (AMS) and constructed a miRNA–mRNA network containing 240 relationship pairs. The hub genes were filtered from the protein-protein interaction network: IL7R, FOS, IL10, FCGR2A, DDX3X, CDK1, BCL11B and HNRNPH1, which were all down-regulated in the AMS group. Then, nine up-regulated DE-miRNAs and eight hub genes were verified by qRT-PCR in our hypoxia-induced HAR cell model. The expression of miR-3177-3p, miR-369-3p, miR-603, miR-495, miR-4791, miR-424-5p, FOS, IL10 and IL7R was consistent with our bioinformatics results. In conclusion, FOS, IL10, IL-7R and 7 DE-miRNAs may participate in the development of HAR. Our findings will contribute to the identification of biomarkers and promote the effective prevention and treatment of HAR in the future.

Su, Gu, Draga, Zhou, Lhamo, Zheng, and Qiu: Integrative analysis of miRNA–mRNA network in high altitude retinopathy by bioinformatics analysis


High-altitude retinopathy (HAR) is one clinical entity of acute high-altitude illness (AAI), which also includes acute mountain sickness (AMS), high-altitude pulmonary edema (HAPE) and high-altitude cerebral edema (HACE) [1,2]. As retinopathy is caused by the inability to adapt to acute oxygen deficiency at high altitudes, HAR mainly manifests as leakage of the peripheral retinal vessels, retinal hemorrhages, tortuous retinal vessels, optic disc edema and macular edema [3,4]. With the development of the economy and tourism in the plateau area, the incidence of HAR is up to 79% and is increasing [5]. However, the pathogenesis of HAR remains unclear, and there is a lack of effective prevention and treatment measures.

Current studies suggest that hypoxia is the main reason for the pathophysiological change in high-altitude situations [6,7]. However, under the same hypoxic conditions, some high-altitude residents will remain in good condition, while others may have AAI [8]. Several genes, including ACE, EDN1, ACYP2, RTEL1, and VEGF, are related to hypoxia [9–11]. In hypoxic environments, transcription of various genes, such as endothelial PAS domain-containing protein 1 (EPAS1) and prolyl hydroxylase domain-containing protein 2 (PHD2), is initiated by hypoxia-related pathways [12]. Also, genetic susceptibility has been reported as one of the major determinants of HAPE and AMS. The oxidative stress-related genes CYBA and GSTP1, contributing to endothelial damage under hypobaric hypoxia, are potential candidate genes for HAPE [13–17]. However, there are a few genetic studies on HAR.

MicroRNAs (miRNAs) can inhibit mRNA translation or induce degradation of target mRNAs by binding to complementary sequences in target regions [18,19]. Huang et al. revealed that saliva miR-134-3p and miR-15b-5p could be used as noninvasive biomarkers to predict AMS individuals exposed to high altitude in advance [19]. In a hypoxic atmosphere, the expression levels of miR-16, 20b, 22, 206, and 17/92 are reduced, which can inhibit ion channels, increase pulmonary artery pressure, and cause vascular dysfunction and the loss of cell integrity, promoting the occurrence of HAPE, and miR-155, 23b, 26a can inhibit TGF signals and promote pulmonary pressure, miR-210 can inhibit mitochondrial function. These may promote the occurrence of HAPE [20]. However, the differential expression of miRNAs has not been reported in HAR.

In the present study, we aimed to explore HAR-related miRNAs and mRNAs via analyzing the GSE90500 and GSE75665 datasets from the Gene Expression Omnibus (GEO). Gene Ontology (GO) functional annotations, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, potential target gene analysis, and construction of potential miRNA–mRNA network were also conducted. Besides, the HAR cell model was established to validate the results of the bioinformatics analysis via quantitative RT-PCR (qRT-PCR). The present study may help elucidate the pathogenic mechanism and identify valuable diagnostic biomarkers for HAR.

Materials and methods

Data source

The miRNA expression profile GSE90500 and the mRNA expression profile GSE75665 were obtained from the GEO database. The RNA-seq data of GSE90500 were based on the platform of GPL18058 (Exiqon miRCURY LNA microRNA array, 7th generation) and contained 13 patients with AMS and 9 non-AMS volunteers. The microarray data of GSE75665 were based on the platform GPL11154 (Illumina HiSeq 2000) and contained 5 patients with AMS and 5 non-AMS volunteers. Figure 1 shows the workflow for the study.

The workflow of bioinformatics analysis
Figure 1
The workflow of bioinformatics analysis

Identification of the differentially expressed miRNA and differentially expressed mRNA

The GEO2R online tool ( based on the R software ‘LIMMA’ package was used to identify differentially expressed miRNA (DE-miRNA) between AMS and non-AMS blood samples of GSE90500 with thresholds of |logFC| > 2 and P -value <0.01. A volcano plot was generated to visualize the significant miRNA expression changes using GraphPad Prism 8.0 (GraphPad, San Diego, CA, USA). For GSE75665, differentially expressed mRNAs (DE-mRNAs) related to AMS were obtained from the study of Liu et al. [21].

The prediction of target genes of DE-miRNA

We used three prediction tools to predict the target mRNAs of DE-miRNA: TargetScan (, miRDB (, and DIANA-microT ( The genes predicted by at least two programs were chosen as the targets of DE-miRNA. We selected the overlapping genes between targets and DE-mRNAs as target differentially expressed genes (DEGs). A heatmap was performed using the R package ‘Pheatmap’. Then, using Cytoscape software [22], the miRNA–mRNA regulatory network was constructed.

Functional enrichment analysis of miRNA-target DEG

GO annotation and KEGG pathway enrichment analyses were performed by utilizing the Database for Annotation, Visualization, and Integrated Discovery (DAVID, online software [23]. P<0.05 was set as the cut-off value. The GO analysis included three categories: biological process (BP), cellular component (CC), and molecular function (MF). Heatmaps of DEGs significantly enriched KEGG pathways was performed by the R package ‘Pheatmap’.

Construction of the protein–protein interaction network

By putting the target DEGs in the Search Tool for the Retrieval of Interacting Genes (STRING, database [24], the interaction relationships among DEGs at the protein level were screened. Then, the protein–protein interaction (PPI) networks for up- and down-regulated target DEGs were constructed under the criterion of a combined score > 0.4. Cytoscape software was used to visualize the networks. The connectivity degrees were calculated through network statistical methods.

Establishment of HAR cell model

Human retinal microvascular endothelial cells (HRMECs) were purchased from Cell Systems Corporations (catalog no. ACBRI181; Kirkland, WA, USA) [25]. The HRMECs were cultured in M199 medium supplemented with 20% fetal bovine serum, 3 ng/ml FGF-basic, 10 units/ml heparin, and 1% streptomycin/penicillin at 37°C. To establish the hypoxia model, cell culture was performed in a hypoxia chamber filled with an anaerobic gas mixture of 94% N2, 5% CO2, and 1% O2 for 12 h. Control groups were cultured in 95% air and 5% CO2.

Quantitative real-time PCR

Total RNA was extracted from HRMECs by using TRIzol reagent. RNA was reverse-transcribed into complementary DNA using PrimeScript RT Master Mix or miRNA First-Strand Synthesis kits (TaKaRa, Kumamoto, Japan), following the manufacturer’s protocols. qRT-PCR was performed by TB Green Premix ExTaq (TaKaRa, Kumamoto, Japan). The oligonucleotide primers used for PCR amplification were purchased from BioSune Biotechnology (Shanghai, China) and are listed in Supplementary Table S1. ACTB and U6 served as internal control for gene and miRNA expression in analysis, respectively. The relative expression of the DE-miRNAs and the DE-mRNAs was calculated using the 2−ΔΔCt method.

Statistical analysis

The data in the present study are presented as the mean ± SD. The results were analyzed by the two-tailed Student’s t test, and P<0.05 was considered to show a statistical difference.


Identification of the DE-mRNA and the DE-miRNA

According to the analysis of Liu et al. for GSE75665 [21], a total of 807 genes were differentially expressed between the non-AMS and AMS groups, with 271 up-regulated and 535 down-regulated DE-mRNAs (Supplementary Table S2). During an analysis of the GSE90500 dataset, we obtained 16 DE-miRNAs (8 up-regulated and 8 down-regulated) by using the GEO2R online tool (P <0.01 and |logFC| > 2) (Table 1). These DE-miRNAs are indicated by a volcano plot (Figure 2). MiR-495 includes miR-495-3p and miR-495-5p, and we discussed them separately in the study, so there were 17 DE-miRNAs (9 up-regulated and 8 down-regulated).

Volcano plot of DE-miRNAs
Figure 2
Red dots indicate up-regulated miRNAs and green dots indicate down-regulated miRNAs with P<0.01 and |logFC| > 2. Black dots indicate the genes which were not differentially expressed in AMS vs. non-AMS.Volcano plot of DE-miRNAs
Table 1
Up- and down-regulated DE-miRNA

MiRNA–mRNA network

The target mRNAs of the DE-miRNAs were predicted by the three programs. Besides, by comparing the targets with DE-mRNAs, we only selected the overlapping genes as target DEGs. As shown in Figure 3A and B, 122 out of 5652 predicted targets of up-regulated DE-miRNAs overlapped down-regulated DE-mRNAs, and 37 out of 3857 predicted targets of down-regulated DE-miRNAs overlapped up-regulated DE-mRNAs. After removing 1 duplicate, we got 157 target DEGs, and the expression of these DEGs was visualized in Figure 3C. Then, we constructed the miRNA–mRNA regulatory network (Figure 4), which consisted of 240 miRNA–mRNA pairs in total. Among them, there were 196 up-regulated miRNA–mRNA pairs and 44 down-regulated miRNA–mRNA pairs. Among the target DEGs, there were 29 transcription factors and 128 non-transcription factors. Table 2 listed the target DEGs of DE-miRNAs.

The mRNAs related to AMS
Figure 3
(A) Overlapped genes between targets of up-regulated DE-miRNAs and down-regulated DE-mRNAs. (B) Overlapped genes between targets of down-regulated DE-miRNAs and up-regulated DE-mRNAs. (C) Heatmap plot of 157 target DEGs. Green means down-regulation, while red means up-regulation.The mRNAs related to AMS
The miRNA–mRNA network in AMS
Figure 4
The diamonds stand for DE-miRNAs, and the ellipses stand for target DEGs. The red and blue colors represent up-regulation and down-regulation, respectively. The yellow ellipses represent transcription factors. The lines indicate the regulation relationship between DE-miRNAs and target DEGs.The miRNA–mRNA network in AMS
Table 2
The target DEGs of DE-miRNA
MiRNAUp/DownCountTarget DEGs
hsa-miR-369-3pUp18IL10, IL7R, ANLN, CCNE2, CDH20, CELF2, CXADR, KCMF1, KL, MDFIC, MSN, NFIB, PDE7B, PTAR1, STEAP2, TLN1, TP63, WNT16
hsa-miR-495-5pUp10DDX3X, DGKH, GOLGA8A, IL6ST, KL, SOX5, TPBG, WWTR1, ZC3H12C, ZFP36L1
hsa-miR-4791Up7FOS, CXADR, FZD8, NUCKS1, RASSF6, SLC38A1, ZNF460
hsa-miR-23b-5pDown3SCUBE1, SPRED3, RFX2
hsa-miR-1304-3pDown10C2ORF48, CDC42EP4, GALNT5, LRP10, NDUFV3, PRICKLE2, RNF152, SATB2, CEACAM8, SEC62
hsa-miR-1183Down8CADM2, DNAJC6, GALNT5, GLUL, ZPBP2, CYP26B1, HDLB, SEC62
hsa-miR-15b-5pDown8BACE2, GDPD5, SPSB4, SYT3, ZBTB44, CARM1, CYP26B1, LYNX1
hsa-miR-1258Down4SEC14L1, TMBIM1, ZBTB44, TMPRSS4
hsa-miR-3144-3pDown3GALNT5, SATB2, FAM166B
hsa-miR-155-5pDown7CDC42EP4, GLUL, MARVELD3, SATB2, ZDHHC2, MYO1D, TRIP13

Functional annotation analysis of the target DEGs

To further explore the biological function of the target DEGs, GO enrichment and KEGG pathway analysis were performed using DAVID. GO analysis results showed that acrosome assembly was the only significantly enriched BP for the up-regulated target DEGs. As shown in Figure 5A and Table 3, the down-regulated target DEGs were mainly enriched in the process of RNA expression, such as positive regulation of transcription from RNA polymerase II promoter, positive regulation of gene expression at BP level. On the CC level, the down-regulated DEGs were mainly enriched in the nucleus, cytoplasm, and plasma membrane. For the MF, the down-regulated DEGs were mainly enriched in sequence-specific DNA binding, transcriptional activator activity, and RNA polymerase II core promoter proximal region sequence-specific binding. In addition, KEGG pathway analyses indicated that down-regulated DEGs were significantly enriched in the Hippo signaling pathway (Figure 5B and Table 4), which is associated with the proliferation and migration of vascular endothelial cells [26,27] and can be deactivated by hypoxia [28]. In addition, as shown in Figure 6A–D, heatmaps illustrated the allocation of down-regulated DEGs significantly enriched four KEGG pathways.

GO enrichment and KEGG pathway analysis of the down-regulated target DEGs
Figure 5
(A) Significant enriched GO terms of down-regulated target DEGs. Green bars stand for the BP, blue bars stand for the CC, and red bars stand for the MF. (B) KEGG pathway analysis of the down-regulated target DEGs.GO enrichment and KEGG pathway analysis of the down-regulated target DEGs
Heatmaps of down-regulated DEGs significantly enriched four KEGG pathways
Figure 6
Viral carcinogenesis (A), pathways in cancer (B), Hepatitis B (C), and Hippo signaling pathway (D). Green means down-regulation, while red means up-regulation.Heatmaps of down-regulated DEGs significantly enriched four KEGG pathways
Table 3
GO enrichment terms of down-regulated target DEGs in AMS
GO IDGO TermCategoryCountP-valueFDR
GO:0045944Positive regulation of transcription from RNA polymerase II promoterBP155.22E-037.91E+00
GO:0006366Transcription from RNA polymerase II promoterBP92.05E-022.78E+01
GO:0010628Positive regulation of gene expressionBP81.78E-032.77E+00
GO:0033077T-cell differentiation in thymusBP52.82E-054.45E-02
GO:0034097Response to cytokineBP44.91E-037.47E+00
GO:0030182Neuron differentiationBP42.51E-023.30E+01
GO:0010468Regulation of gene expressionBP42.86E-023.67E+01
GO:0035264Multicellular organism growthBP41.60E-022.24E+01
GO:0048041Focal adhesion assemblyBP31.09E-021.58E+01
GO:0048538Thymus developmentBP33.28E-024.09E+01
GO:0042100B-cell proliferationBP31.89E-022.60E+01
GO:1900740Positive regulation of protein insertion into mitochondrial membrane involved in apoptotic signaling pathwayBP31.67E-022.33E+01
GO:0008625Extrinsic apoptotic signaling pathway via death domain receptorsBP32.61E-023.41E+01
GO:0045727Positive regulation of translationBP34.80E-025.39E+01
GO:0031532Actin cytoskeleton reorganizationBP33.86E-024.62E+01
GO:0034446Substrate adhesion-dependent cell spreadingBP32.61E-023.41E+01
GO:0032835Glomerulus developmentBP31.51E-032.36E+00
GO:0072577Endothelial cell apoptotic processBP24.54E-025.19E+01
GO:0046632α-β T-cell differentiationBP23.90E-024.66E+01
GO:0005886Plasma membraneCC361.88E-022.07E+01
GO:0005925Focal adhesionCC81.00E-021.15E+01
GO:0005913Cell–cell adherens junctionCC71.42E-021.60E+01
GO:0015629Actin cytoskeletonCC54.44E-024.24E+01
GO:0030018Z discCC43.54E-023.56E+01
GO:0044300Cerebellar mossy fiberCC23.01E-023.11E+01
GO:0043565Sequence-specific DNA bindingMF91.81E-022.14E+01
GO:0001077Transcriptional activator activity, RNA polymerase II core promoter proximal region sequence-specific bindingMF88.05E-041.05E+00
GO:0098641Cadherin binding involved in cell-cell adhesionMF71.08E-021.34E+01
GO:0000166Nucleotide bindingMF72.44E-022.77E+01
GO:0003779Actin bindingMF63.31E-023.58E+01
GO:0003730mRNA 3′-UTR bindingMF43.79E-034.87E+00
Table 4
Enriched KEGG pathway of down-regulated target DEGs associated with AMS
hsa04390Hippo signaling pathway62.65E+00AJUBA, FZD8, WNT16, RASSF6, YWHAB, WWTR1
hsa05161Hepatitis B62.23E+00CCNE2, FOS, DDX3X, BCL2, MAP3K1, YWHAB
hsa05200Pathways in cancer73.27E+01CCNE2, FZD8, FOS, WNT16, GNAI1, BCL2, MMP1
hsa05203Viral carcinogenesis53.60E+01CCNE2, CDK1, DDX3X, IL6ST, YWHAB

PPI network of the target DEGs

Based on information from the STRING database, the PPI network was constructed and contained 59 nodes and 99 interactions (Figure 7). The details of the PPI network are described in Supplementary Table S3. In this network, with the criteria of filtering degree > 5, eight genes were defined as hub genes, including IL10, CDK1, FOS, HNRNPH1, IL7R, BCL11B, FCGR2A, and DDX3X (Table 5).

PPI network construction and module analysis
Figure 7
The red nodes represent the up-regulated DEGs. The blue nodes represent the down-regulated DEGs. The size of the nodes indicates the degree of the DEGs, and the lines represent the interaction between gene-encoded proteins.PPI network construction and module analysis
Table 5
Top eight highest degree nodes

Validation of hub genes in HAR cell model

To validate the results of bioinformatics analysis, we chose the nine up-regulated DE-miRNAs and eight hub genes for qRT-PCR under normoxia and hypoxia in HRMECs. As shown in Figure 8A and B, miR-3177-3p, miR-369-3p, miR-603, miR-495-3p, miR-495-5p, miR-4791, and miR-424-5p were significantly up-regulated in hypoxia-cultured HRMECs compared with normoxia (P<0.05). FOS, IL10, and IL7R were strongly down-regulated under hypoxia conditions (P<0.05), which was also consistent with our bioinformatics results. However, the expression of DDX3X, BCL1B, HNRNPH1, and CDK1 was higher than that in normoxia-cultured HRMECs.

Relative expression levels of DE-miRNAs and hub genes in HRMECs under normoxia and hypoxia conditions by qRT-PCR
Figure 8
(A) The expression of nine up-regulated miRNAs. (B) The expression of eight DEGs. *P<0.05; **P<0.01; ***P<0.001.Relative expression levels of DE-miRNAs and hub genes in HRMECs under normoxia and hypoxia conditions by qRT-PCR


HAR is an ocular manifestation of AAI [1], and studies have found a statistically significant correlation between HAR and HACE [3]. Some researchers have proposed that AMS is a mild form of HACE [29,30]. Since the retina and optic nerve act as the directly visible part of the brain, belonging to the central nervous system due to its embryonic origin [31,32], we obtained the DE-miRNAs and DE-mRNAs of HAR by analyzing the AMS datasets and then explored the possible potential target genes of HAR.

MiRNAs can regulate many BPs by negatively regulating the expression of target genes and are thus important targets for studying the genetic susceptibility of diseases [19,33]. Through a series of bioinformatics analyses and experiments, we selected 16 DE-miRNAs and 157 target DEGs related to AMS from GSE90500 and GSE75665 and constructed an miRNA–mRNA network containing 240 relationship pairs. The top eight hub genes—IL7R, FOS, IL10, FCGR2A, DDX3X, CDK1, BCL11B and HNRNPH1—filtered from the PPI network may serve as candidate biomarkers of HAR. The expression of miR-3177-3p, miR-369-3p, miR-603, miR-495-3p, miR-495-5p, miR-4791, miR-424-5p, FOS, IL10, and IL7R were confirmed by qRT-PCR in our HAR cell model to have the same trend as that found in our bioinformatics results. These miRNAs and mRNAs may play an important role in the development of HAR.

FOS, consisting of c-Fos, FosB, Fra1, and Fra2 [34], was significantly down-regulated in our hypoxia-cultured HRMECs. In the lung tissues of hypoxic rats, the expression of c-Fos and FosB was significantly decreased [35], which is consistent with our results. c-Fos promotes the development of inflammatory diseases such as arthritis [36], but it also inhibits inflammation in myeloid and lymphoid cell lineages [37]. In studies on retinal development, it has been found that c-Fos induces inflammatory signals in photoreceptor cells, leading to an increase in VEGF, thus promoting the formation of new blood vessels [38]. Therefore, FOS may be associated with the leakage of the peripheral retinal vessels and retinal hemorrhages in HAR through the regulation of inflammatory constituents. The bioinformatics analysis showed that miR-603 and miR-4791 may down-regulate FOS, and miR-603 can promote the growth of glioma cells via the Wnt/β-catenin pathway [39]. The Wnt pathway mediates many complex biological and pathological processes of retinopathy [40], but there are no reports about miR-603 and this pathway in the development of HAR. Further research about the role of FOS and miR-603 in HAR is needed.

Another inflammation-related gene, IL10, could reduce tissue damage and promote immune tolerance [41]. It was down-regulated in our study, and maybe serve as a biomarker in HAR. Under hypoxic conditions, the decrease in IL10 indicates hypoxia reduces the ability to inhibit immunity and inflammation [42]. However, in the study of Kang et al., mice adapted to the low oxygen environment (10% O2 ) had higher IL10 levels [43]. Therefore, the role of anti-inflammatory markers IL10 in HAR needs further investigation.

IL7R mediates the signaling of IL7, affecting the role of IL7 in regulating the development of T and B cells and the balance of the internal environment of T cells [44], and IL7R is involved in the pathogenesis of several autoimmune diseases, such as rheumatoid arthritis, multiple sclerosis, and systemic lupus erythematosus etc. [45] High altitude also could cause changes in different immune cells [46]. According to our bioinformatics analysis, we found that IL7R was regulated by miR-369-3p and miR-424-5p. It is noted that the induction of miR-369-3p expression inhibits the release of inflammatory cytokines [47]. Furthermore, up-regulated miR-424-5p in hypoxic vascular endothelial cells was shown to stabilize hypoxia-inducible factor-1α (HIF-1α) and fine-tune VEGF, promoting angiogenesis and preserving the integrity of the endothelial barrier [48]. Although the down-regulated IL7R and the up-regulated miR-369-3p and miR-424-5p may play an important role in the progression of HAR, the specific mechanism is not clear and remains to be clarified.

High-altitude retinal hemorrhages and vascular leakage may be caused by the retinal microcirculation disorder under hypoxia conditions at high altitude, but the specific mechanism is not clear. The expression of DDX3X was found to be correlated with overexpression of HIF-1α in breast cancer, indicating the oncogenic role of DDX3X [49]. CDK1 has been shown to promote tumor angiogenesis by stabilizing HIF-1α, and its disruption could inhibit retinal angiogenesis by inducing cell cycle arrest and apoptosis [50]. However, inhibition of miR-495 improves angiogenesis and blood flow recovery after ischemia [51]. Although the expression of DDX3X and CDK1 did not down-regulate in HRMECs as predicted, these genes had high degrees in the PPI network, indicating their potential roles in the development of HAR. It worth mentioning that very little was reported on miR-3177-3p, so the specific function of miR-3177-3p needs further research.

To our knowledge, our study is the first to explore the DE-miRNAs and DEGs of HAR and establish the miRNA–mRNA network through comprehensive bioinformatics analyses. However, there are still some limitations in the present study. First, the data of the GSE90500 and GSE75665 are collected from blood samples, which may be the reason why the expression of several hub genes in HRMECs was inconsistent with the bioinformatics analyses. Then, there are some differences between the analytical methods of GSE90500 and GSE75665, which reduces the reliability of this study to a certain extent. Third, the miRNA–mRNA relationships in HAR were based on target prediction and statistical evidence, which need experimental investigations. Thus, Further studies are required to verify the specific functions of the identified DE-miRNAs and hub genes in HAR.

In conclusion, the present study identified FOS, IL10, IL7R and seven miRNAs as candidate biomarkers of HAR, which may participate in the development of HAR through possible pathways. It has been proposed that HAR is related to individual sensitivity, promoting the effective prevention, and treatment of HAR in the future. Nevertheless, we carried out the systematic and comprehensive bioinformatics analyses to identify the new miRNA and genes of HAR and performed experimental validation in the HAR cell model, which is helpful for understanding the gene changes in HAR and provides novel candidate biomarkers for the diagnosis of HAR.

Data Availability

All data associated with the present paper are already included in the manuscript.

Competing Interests

The authors declare that there are no competing interests associated with the manuscript.


This work was supported by the National Natural Science Foundation of China [grant numbers 81970811, 81770947]; the National Science and Technology Major Project of China [grant number 2017ZX09304010]; the National Key R&D Program of China [grant numbers 2016YFC0904800, 2019YFC0840607]; and the Tibet Natural Science Foundation of China [grant numbers XZ2017ZR-ZYZ09, XZ2018ZRG-95, XZ2018ZRG-100(Z)].

Author Contribution

All authors contributed significantly to this work. Tong Su, Chufeng Cu, and Deji Draga were involved in all aspects of research, analyzed data and wrote the paper. Chuandi Zhou performed the experiments. Thashi Lhamo was involved in research design. Qinghua Qiu and Zhi Zheng contributed to research design, supervision of research, data analysis and revised the manuscript. All authors have read and approved the final manuscript.


AAIacute high-altitude illness
AMSacute mountain sickness
BPbiological process
CCcellular component
DAVIDDatabase for Annotation, Visualization, and Integrated Discovery
DE-miRNAdifferentially expressed miRNA
DE-mRNAdifferentially expressed mRNA
DEGdifferentially expressed gene
GEOGene Expression Omnibus
GOGene Ontology
HACEhigh-altitude cerebral edema
HAPEhigh-altitude pulmonary edema
HARhigh-altitude retinopathy
HIF-1αhypoxia-inducible factor-1α
HRMEChuman retinal microvascular endothelial cell
KEGGKyoto Encyclopedia of Genes and Genomes
MFmolecular function
qRT-PCRuantitative RT-PCR
STRINGSearch Tool for the Retrieval of Interacting Genes



Houston C.S. (1976) High altitude illness. Disease with protean manifestations. JAMA 236, 21932195 doi: 10.1001/jama.1976.03270200031025


Seth R.K., and Adelman R.A. (2010) High-altitude retinopathy and optical coherence tomography findings. Semin. Ophthalmol. 25, 1315 doi: 10.3109/08820538.2010.481560


Wiedman M., and Tabin G.C. (1999) High-altitude retinopathy and altitude illness. Ophthalmology 106, 19241926, discussion 1927 doi: 10.1016/S0161-6420(99)90402-5


Tian X., et al. (2018) Retinal changes following rapid ascent to a high-altitude environment. Eye (Lond.) 32, 370374 doi: 10.1038/eye.2017.195


Honigman B., et al. (2001) High altitude retinal hemorrhages in a Colorado skier. High Alt. Med. Biol. 2, 539544 doi: 10.1089/152702901753397117


McFadden D.M., et al. (1981) High-altitude retinopathy. JAMA 245, 581586 doi: 10.1001/jama.1981.03310310023016


Li K., et al. (2017) Transcriptome reveals the overexpression of a kallikrein gene cluster (KLK1/3/7/8/12) in the Tibetans with high altitude-associated polycythemia. Int. J. Mol. Med. 39, 287296 doi: 10.3892/ijmm.2016.2830


MacInnis M.J., and Koehle M.S. (2016) Evidence for and against genetic predispositions to acute and chronic altitude illnesses. High Alt. Med. Biol. 17, 281293 doi: 10.1089/ham.2016.0024


Jin T., et al. (2019) Association between the IL1R2 rs2072472 polymorphism and high-altitude pulmonary edema risk. Mol. Genet. Genom. Med. 7, e542 doi: 10.1002/mgg3.542


Hanaoka M., et al. (2009) Polymorphisms of human vascular endothelial growth factor gene in high-altitude pulmonary oedema susceptible subjects. Respirology 14, 4652 doi: 10.1111/j.1440-1843.2008.01420.x


Rong H., et al. (2017) Association between regulator of telomere elongation helicase1 (RTEL1) gene and HAPE risk: a case-control study. Medicine (Baltimore) 96, e8222 doi: 10.1097/MD.0000000000008222


Pan H., et al. (2018) Mutations in EPAS1 in congenital heart disease in Tibetans. Biosci. Rep. 38, 6 doi: 10.1042/BSR20181389.


He Y., et al. (2016) Telomere length-related gene ACYP2 polymorphism is associated with the risk of HAPE in Chinese Han population. J. Gene Med. 18, 244249 doi: 10.1002/jgm.2896


Mishra A., et al. (2012) CYBA and GSTP1 variants associate with oxidative stress under hypobaric hypoxia as observed in high-altitude pulmonary oedema. Clin. Sci. (Lond.) 122, 299309 doi: 10.1042/CS20110205


Richalet J.P., et al. (2012) Physiological risk factors for severe high-altitude illness: a prospective cohort study. Am. J. Respir. Crit. Care Med. 185, 192198 doi: 10.1164/rccm.201108-1396OC


Schneider M., et al. (2002) Acute mountain sickness: influence of susceptibility, preexposure, and ascent rate. Med. Sci. Sports Exerc. 34, 18861891 doi: 10.1097/00005768-200212000-00005


Bloch K.E., et al. (2009) Effect of ascent protocol on acute mountain sickness and success at Muztagh Ata, 7546 m. High Alt. Med. Biol. 10, 2532 doi: 10.1089/ham.2008.1043


Lou W., et al. (2018) MicroRNA regulation of liver cancer stem cells{Lou, 2018 #42}. Am J. Cancer Res. 8, 11261141


Nunez-Iglesias J., et al. (2010) Joint genome-wide profiling of miRNA and mRNA expression in Alzheimer’s disease cortex reveals altered miRNA regulation. PLoS ONE 5, e8898 doi: 10.1371/journal.pone.0008898


Alam P., , Saini N., and Pasha M.A. (2015) MicroRNAs: an apparent switch for high-altitude pulmonary edema. MicroRNA 4, 158167 doi: 10.2174/2211536604666151103121633


Liu B., et al. (2017) IL-10 dysregulation in acute mountain sickness revealed by transcriptome analysis. Front. Immunol. 8, 628 doi: 10.3389/fimmu.2017.00628


Smoot M.E., et al. (2011) Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics 27, 431432 doi: 10.1093/bioinformatics/btq675


Huang da W., , Sherman B.T., and Lempicki R.A. (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 4457 doi: 10.1038/nprot.2008.211


Szklarczyk D., et al. (2011) The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 39, D561D568 doi: 10.1093/nar/gkq973


Gu C., et al. (2019) miR-590-3p inhibits pyroptosis in diabetic retinopathy by targeting NLRP1 and inactivating the NOX4 signaling pathway. Invest. Ophthalmol. Vis. Sci. 60, 42154223 doi: 10.1167/iovs.19-27825


Moon K.H., and Kim J.W. (2018) Hippo signaling circuit and divergent tissue growth in mammalian eye. Mol. Cells 41, 257263


Kim J., et al. (2017) YAP/TAZ regulates sprouting angiogenesis and vascular barrier maturation. J. Clin. Invest. 127, 34413461 doi: 10.1172/JCI93825


Ma B., et al. (2015) Hypoxia regulates Hippo signalling through the SIAH2 ubiquitin E3 ligase. Nat. Cell Biol. 17, 95103 doi: 10.1038/ncb3073


Bartsch P., et al. (2004) Acute mountain sickness: controversies and advances. High Alt. Med. Biol. 5, 110124 doi: 10.1089/1527029041352108


Hackett P.H., and Roach R.C. (2001) High-altitude illness. N. Engl. J. Med. 345, 107114 doi: 10.1056/NEJM200107123450206


Hughes S., , Yang H., and Chan-Ling T. (2000) Vascularization of the human fetal retina: roles of vasculogenesis and angiogenesis. Invest. Ophthalmol. Vis. Sci. 41, 12171228


Patton N., et al. (2005) Retinal vascular image analysis as a potential screening tool for cerebrovascular disease: a rationale based on homology between cerebral and retinal microvasculatures. J. Anat. 206, 319348 doi: 10.1111/j.1469-7580.2005.00395.x


Beermann J., et al. (2016) Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches. Physiol. Rev. 96, 12971325 doi: 10.1152/physrev.00041.2015


Tkach V., et al. (2003) Role of the Fos family members, c-Fos, Fra-1 and Fra-2, in the regulation of cell motility. Oncogene 22, 50455054 doi: 10.1038/sj.onc.1206570


Singh M., et al. (2018) The MAPK-activator protein-1 signaling regulates changes in lung tissue of rat exposed to hypobaric hypoxia. J. Cell. Physiol. 233, 68516865 doi: 10.1002/jcp.26556


Shiozawa S., and Tsumiyama K. (2009) Pathogenesis of rheumatoid arthritis and c-Fos/AP-1. Cell Cycle 8, 15391543 doi: 10.4161/cc.8.10.8411


Ray N., et al. (2006) c-Fos suppresses systemic inflammatory response to endotoxin. Int. Immunol. 18, 671677 doi: 10.1093/intimm/dxl004


Sun Y., et al. (2017) Inflammatory signals from photoreceptor modulate pathological retinal angiogenesis via c-Fos. J. Exp. Med. 214, 17531767 doi: 10.1084/jem.20161645


Guo M., et al. (2015) miR-603 promotes glioma cell growth via Wnt/beta-catenin pathway by inhibiting WIF1 and CTNNBIP1. Cancer Lett. 360, 7686 doi: 10.1016/j.canlet.2015.02.003


Huang H.M., et al. (2015) Activating the Wnt/beta-catenin pathway did not protect immature retina from hypoxic-ischemic injury. Invest. Ophthalmol. Vis. Sci. 56, 43004308 doi: 10.1167/iovs.14-16176


Hilgenberg E., et al. (2014) Interleukin-10-producing B cells and the regulation of immunity. Curr. Top. Microbiol. Immunol. 380, 6992 doi: 10.1007/978-3-662-43492-5_4


Yuhai G.U., and Zhen Z. (2015) Significance of the changes occurring in the levels of interleukins, SOD and MDA in rat pulmonary tissue following exposure to different altitudes and exposure times. Exp. Ther. Med. 10, 915920 doi: 10.3892/etm.2015.2604


Kang J.G., et al. (2016) Low ambient oxygen prevents atherosclerosis. J. Mol. Med. (Berl.) 94, 277286 doi: 10.1007/s00109-016-1386-3


Mai H.L., et al. (2014) IL-7 receptor blockade following T cell depletion promotes long-term allograft survival. J. Clin. Invest. 124, 17231733 doi: 10.1172/JCI66287


Wang X.S., et al. (2013) Perspectives of the relationship between IL-7 and autoimmune diseases. Clin. Rheumatol. 32, 17031709 doi: 10.1007/s10067-013-2360-x


Khanna K., et al. (2018) High-altitude-induced alterations in gut-immune axis: a review. Int. Rev. Immunol. 37, 119126 doi: 10.1080/08830185.2017.1407763


Galleggiante V., et al. (2019) Quercetin-induced miR-369-3p suppresses chronic inflammatory response targeting C/EBP-beta. Mol. Nutr. Food Res. 63, e1801390 doi: 10.1002/mnfr.201801390


Ghosh G., et al. (2010) Hypoxia-induced microRNA-424 expression in human endothelial cells regulates HIF-alpha isoforms and promotes angiogenesis. J. Clin. Invest. 120, 41414154 doi: 10.1172/JCI42980


Bol G.M., et al. (2013) Expression of the RNA helicase DDX3 and the hypoxia response in breast cancer. PLoS ONE 8, e63548 doi: 10.1371/journal.pone.0063548


Gao X., et al. (2019) Cyclin-dependent kinase 1 disruption inhibits angiogenesis by inducing cell cycle arrest and apoptosis. Exp. Ther. Med. 18, 30623070 doi: 10.3892/etm.2019.7883


Welten S.M., et al. (2014) Inhibition of 14q32 microRNAs miR-329, miR-487b, miR-494, and miR-495 increases neovascularization and blood flow recovery after ischemia. Circ. Res. 115, 696708 doi: 10.1161/CIRCRESAHA.114.304747
This is an open access article published by Portland Press Limited on behalf of the Biochemical Society and distributed under the Creative Commons Attribution License 4.0 (CC BY). analysis of miRNA–mRNA network in high altitude retinopathy by bioinformatics analysis&author=Tong Su,Chufeng Gu,Deji Draga,Chuandi Zhou,Thashi Lhamo,Zhi Zheng,Qinghua Qiu,&keyword=Bioinformatics analysis,Differentially expressed genes,High-altitude retinopathy,MiRNA-mRNA network,&subject=Bioinformatics,Gene Expression & Regulation,Molecular Interactions,Research Articles,