Skip to main content


DNA methylation exploration for ARDS: a multi-omics and multi-microarray interrelated analysis



Despite advances in clinical management, there are currently no novel therapeutic targets for acute respiratory distress syndrome (ARDS). DNA methylation, as a reversible process involved in the development and progression of many diseases, would be used as potential therapeutic targets to improve the treatment strategies of ARDS. However, the meaningful DNA methylation sites associated with ARDS still remain largely unknown. We sought to determine the difference in DNA methylation between ARDS patients and healthy participants, and simultaneously, the feasible DNA methylation markers for potential therapeutic targets were also explored.


Microarray data of human blood samples for ARDS and healthy participants up to June 2019 was searched in GEO database. The difference analyses between ARDS and healthy population were performed through limma R package, and furthermore, interrelated analyses of DNA methylation and transcript were accomplished by VennDiagram R package. Perl and sva R package were used to merge microarray data and decrease heterogeneities among different studies. The biological function of screened methylation sites and their regulating genes were annotated according to UniProt database and Pubmed database. GO term and KEGG pathway enrichment analyses were conducted using DAVID 6.8 and KOBAS 3.0. The meaningful DNA methylation markers to distinguish ARDS from healthy controls were explored through ROC (receiver operating characteristic curves) analyses.


Five datasets in GEO databases (one DNA methylation dataset, three mRNA datasets, and one mRNA dataset of healthy people) were enrolled in present analyses finally, and the series were GSE32707, GSE66890, GSE10474, GSE61672, and GSE67530. These databases included 99 patients with ARDS (within 48 h of onset) and 136 healthy participants. Difference analyses indicated 44,439 DNA methylation alterations and 29 difference mRNAs between ARDS and healthy controls. 40 methylation variations regulated transcription of 16 genes was explored via interrelated analysis. According to the functional annotations, 30 DNA methylation sites were related to the imbalance of inflammation or immunity, endothelial function, epithelial function and/or coagulation function. cg03341377, cg24310395, cg07830557 and cg08418670, with AUC up to 0.99, might be the meaningful characteristics with the highest performance to distinguish ARDS from healthy controls.


44,439 DNA methylation alterations and 29 difference mRNAs exist between ARDS and healthy controls. 30 DNA methylation sites may regulate transcription of 10 genes, which take part in pathogenesis of ARDS. These findings could be intervention targets, with validation experiments to be warranted to assess these further.


Acute respiratory distress syndrome (ARDS) is a life-threatening form of respiratory failure that accounts for 10% of intensive care unit admissions. In spite of improvements in basic science and clinical research, treatment options for ARDS are still limited and the mortality of severe ARDS is 40–46% [1]. Most innovative therapies for ARDS have failed in the last decade, and clearly, there is a robust need for better insight in disease pathogenesis and subsequent emerging treatment strategies [1, 2].

The majority of studies focus on genomic or transcriptomic responses in ARDS [3,4,5]. Accumulating evidences demonstrate that the epigenetic alterations especial DNA methylation involved in the development and progression of many diseases, including various cancers, lupus, diabetes, asthma, and a variety of neurological disorders [6, 7]. It was reported that DNA methylation can directly block transcription by inhibiting the binding of specific transcription factors to their target sequences on the candidate gene, which could result in an obvious variation in transcriptome and lead to the occurrence and development of diseases eventually [8,9,10]. In this scenario, valuable DNA methylation variations could be used as biomarkers for molecular diagnosis of disease. Moreover, DNA methylation could also be potential target to explore new treatment strategy [6,7,8,9,10]. For instance, Dhas et al. showed that specific DNA methylation might be potential marker for prognosis of neonatal sepsis which may improve the treatment strategies [11].

To date, studies on DNA methylation about ARDS is insufficient. Szilagyi et al. showed that myosin light chain kinase (MYLK) epigenetic variations were implicated in ARDS pathogenesis and might influence ARDS mortality [12]. In this study, researchers only focus on single genetic variations and the other underlying methylated CpGs related genes still unclear. Fortunately, the author had uploaded the DNA methylation microarray data of patients with ARDS to the Gene Expression Omnibus (GEO) databases, which made it possible for us to explore DNA methylation in ARDS comprehensively. Beside these, there are three mRNA microarray data of patients with ARDS in GEO databases which provided by Kangelaris et al. [3], Dolinay et al. [4] and Howrylak et al. [5]. All of these public data will be helpful to us perform integrative analyses of DNA methylation and mRNA in ARDS.

In the current study, in order to find DNA methylation alterations and mRNA expression differences between patients with ARDS and health controls, we utilized difference analysis on DNA methylation microarray and mRNA microarray data. Furthermore, for exploration of the potentially meaningful DNA methylation alterations, we performed integrative analyses of DNA methylation variations and mRNA expression differences. In addition, to evaluate the value of these methylation alterations to distinguish ARDS, receiver operating characteristic curves (ROC) would be calculated.


Systematic search and data selection

The GEO was searched for all expression microarray that matched terms of ARDS. Clinical studies of ARDS using peripheral blood of adult were retained. Datasets that utilized endotoxin or lipopolysaccharide infusion as vitro or animal models for ARDS were excluded. Clinical-gene expression microarray derived from sorted cells was also excluded (Fig. 1).

Fig. 1

The flow-process diagram. The flow-process diagram on the left was how to screen microarray datasets

Expression microarray datasets of healthy participants were searched and set as the control group that matched terms of health and GPL10558 (the same platform number as included gene expression microarray of ARDS). Moreover, where longitudinal data were available for patients admitted with ARDS, we only included data derived from the within 48 h of onset.

Gene expression normalization

ComBat normalization in the sva R package and Perl was used to co-normalize these cohorts into a single cohort, after re-normalizing from raw data. All datasets were downloaded as txt files, and DNA methylation microarray were re-normalized and batch corrected through minfi, impute, and wateRmelon R packages. Outputs from mRNA array were normal-exponential background corrected and then between-arrays quantile normalized using limma R package. For compatibility with microarray studies, expression was normalized using a weighted linear regression, and the estimated precision weights of each observation were then multiplied with the corresponding log2 to yield final gene expression values.

Integrative analysis

There are two parts to the integrative analysis shown in Additional file 1: Figure S1: on the one hand, the differentially expressed DNA methylation sites and mRNAs were determined using the R package limma, which implements an empirical Bayesian approach to estimate gene expression changes using moderated t-tests. A log fold change > 0 was defined as hypermethylation or up-regulated genes, and a log fold change < 0 was defined as demethylation or down-regulated genes. For elimination of the influence of human species on DNA methylation, the intersection analyses of the DNA methylation alterations of the melanoderm and Caucasian groups were taken, Additional file 2: Figure S2.

On the other hand, the methylation sites with biological function and methylated genes were defined as hypermethylation accompanying with down-regulated genes and demethylation accompanying with up-regulated genes, simultaneously. In this scenario, we respectively took the intersection of hypermethylation sites and down-regulated genes, and the intersection of demethylation sites and up-regulated genes, Additional file 3: Figure S3.

The differentially expressed DNA methylation sites and mRNA were identified by significance criteria (adjusted P value < 0.05) as implemented in the R package limma. The intersection analyses were performed through VennDiagram R package.

Functional annotation and enrichment of screened methylation alterations

Functional annotations and gene-annotation enrichment analyses of genes regulated by screened methylation variations were referenced the UniProt database, KOBAS 3.0 online database, DAVID 6.8 online database and PubMed database.

For the exploration of the interaction of screened genes, STRING online database was used to plot the protein–protein interaction network (PPI).

Evaluation the value of methylation alterations to distinguish ARDS from healthy volunteers

To identify novel methylation alterations which would to distinguish ARDS from healthy volunteers, the receiver operating characteristic curves (ROCs) were performed to calculate the area under the curve (AUC) on screened DNA methylation variations using ROCR R package.

Software and versions

Perl 64 was used to merge data; R x64 3.4.4 was conducted to process data, analyse data and plot diagrams; Cytoscape 3.6.1 was performed to plot network diagrams.


Characteristics of the datasets

After search strategy and inclusive criteria, 5 studies containing one DNA methylation dataset (GSE67530) [12], 3 mRNA datasets of patients with ARDS (GSE32707, GSE66890 and GSE10474) [3,4,5], and one mRNA dataset of healthy people (GSE61672) [13] were used to build the methylation and mRNA expression profiling datasets, Table 1.

Table 1 Microarray data information

DNA methylation microarray was the M-values (log2 ratio of the intensities of modified probe vs unmodified probe) of CpG probes (485,577 probes) that passed quality control and batch corrected. Gene microarray datasets were mRNA expression profiling after quality control and batch correction, with the median number of mRNA probes assayed as 25,128 (ranging from 22,277 to 47,220). Our integrated dataset included a total of 99 patients with ARDS and 136 health participants. The number of samples investigated ranged from 13 to 106 cases (median 30) across the studies. All datasets were from United States, and there was no statistical difference among datasets in age, and male.

Screening of differentially expressed genes and methylation sites

In total, 22,654 hypermethylation sites and 21,785 demethylation sites were identified in Manhattan chart (Fig. 2). Sorted by adjusted P value from small to large, the top 10 methylation variations between patients with ARDS and healthy volunteers were cg17078393, cg04794690, cg01564818, cg07748255, cg07369374, cg23856138, cg21726551, cg25032321, cg10431989 and cg26852712. Obviously, 9 up-regulated mRNA and 20 down-regulated mRNA were identified in volcano plot (Fig. 2 and Additional file 4: Table S1).

Fig. 2

Differential methylation sites and mRNA expression. The Manhattan Figure at the top shows the differential methylation sites between the patients with ARDS and the healthy controls. The red line shows the 10 sites with minimum P values. The volcano figure at the bottom shows the mRNA expression. The red points indicate high-expression mRNA in ARDS, and the green points indicate low-expression mRNA according to the threshold of the FDR

Exploration of methylation sites with biological function and methylated genes

16 genes were differentially methylated by interrelated analyses criteria, including 32 hypermethylation sites and 8 demethylation sites, demonstrated in Table 2 and Fig. 3. Overall, according to the function annotation, 30 methylation alterations were finally screened to be involved in the pathogenesis of ARDS through their potentially regulated 10 genes (Fig. 4), with the mRNA expression of 10 potentially regulated genes demonstrated in Additional file 5: Figure S4.

Table 2 Relationship between methylation sites and mRNA, with AUC of methylation sites
Fig. 3

Circle diagram for the interrelated analysis and PPI analysis results. This diagram shows the positions of methylation sites and their target genes on 24 chromosomes after screening by preliminary interrelated analysis between DNA methylation and mRNA microarray datasets. The red lines represent the relationship between DNA methylation and their target genes. The blue lines represent interactions between proteins coded by these genes according to Protein Interaction Analyses

Fig. 4

Network diagram for the relationship between the screened genes and the 4 mechanisms of ARDS. Colour depth represents the log-fold change of genes. Function represents the 4 mechanisms of ARDS. Genes in cluster A was confirmed to be related to the pathogenesis of ARDS in both ARDS animal models and in vitro experiments. Genes in cluster B was confirmed to be related to the pathogenesis of inflammation or immunity, endothelial function, epithelial function, and coagulation function in in vitro experiments. Genes in cluster C might be related to the above mechanism because these genes are related to immune regulation, angiogenesis and the cell cycle

Peptidase inhibitor 3 (PI3), C-X3-C motif chemokine receptor 1 (CX3CR1) and FYN proto-oncogene, Src family tyrosine kinase (FYN) have been confirmed to be related to the pathogenesis of inflammation or immunity, endothelial function, epithelial function, and coagulation function in both ARDS animal models and in vitro experiments. It was noteworthy that the mRNA level of PI3 in the plasma of patients with ARDS was significantly decreased.

The dual specificity phosphatase 6 (DUSP6), DNA polymerase beta (POLB); and solute carrier family 3 member 2 (SLC3A2) genes only had in vitro experimental evidence, and paired immunoglobin like type 2 receptor alpha (PILRA), SRSF protein kinase 2 (SRPK2), ring finger protein 1 (RNF1) and tripartite motif containing 33 (TRIM33) might be related to the pathogenesis of ARDS, as shown in Additional file 6: Table S2.

Gene Ontology (GO) and pathway enrichment analyses indicated that the 10 screened genes were enriched in 9 functions and 25 pathways, as shown in Additional file 7: Figure S5, Additional file 8: Figure S6. These functions and pathways included 2 types. The first type was related to endothelial and epithelial apoptosis and repair, including mTOR and MAPK pathways. The second type was related to inflammation or immunity, including adaptive immune response and leukocyte migration.

Assessment of diagnostic efficacy on screened methylation sites

AUC was calculated on screened DNA methylation alterations based on ROC analyses, Table 2. The AUC on cg03341377, cg24310395, cg07830557 and cg08418670 were up to 0.99, which meant that these DNA methylation alterations could be most significant characteristics to distinguish ARDS from control subjects.


DNA methylation has been most studied and is classically associated with gene silencing via hypermethylation of CpG islands located in promoter regions of various diseases suppressor genes. In addition, DNA methylation is an attractive investigative tool for the study given that methylation is a reversible process. For instance, in myelodysplastic syndrome, demethylating agents have been tested in phase I clinical trials [14]. However, DNA methylation research on ARDS was inadequate and recent studies only focused on single genetic variations [15, 16]. Detailed studies of other novel candidates might lead to identification of unsuspected evolutionarily conserved mechanisms triggered by ARDS.

To our knowledge, this was the first study of all-encompassing analysis for DNA methylation alterations in ARDS, indicated in Additional file 9: Table S3. As furthermore research targets for improvement of therapies in ARDS, the present study uncovered 30 methylation alterations and their regulating DNA in ARDS regions. Meanwhile, as diagnostic molecules for ARDS, the clinicians might be interested in the top 10 difference methylation sites and 4 high diagnostic efficacy CpG islands variations.

Due to microarray experiments from human blood samples as the object of this study, 7 animal model studies and 10 studies not on blood samples were eliminated, leaving only 4 studies of microarray experiments and 1 study of transcriptome sequencing. In the consideration of comorbidity of hematologic malignancy existing in patients of transcriptome sequencing study, the sequencing data were eventually removed. Previous studies utilized only one single database, with a relatively small sample size. Transcriptome microarray datasets were not found healthy participants, which means in the need of searching healthy volunteers data. The data of healthy participants (generalized anxiety disorder score = 0) from one dataset related to anxiety disorder was extracted, according to matched terms of health and GPL10558 (the same platform number as included gene expression microarray of ARDS) in consideration of compatibility with microarray studies.

After merging these cohorts into a single cohort by Perl, batch correction and co-normalization were performed through ComBat normalization in the sva R package, which was a frequently-used method to decease heterogeneity among microarray studies [17,18,19,20]. In addition, the threshold values of screening differential mRNAs were P value < 0.05 in previous studies, which might lead to partially false positives. The significance criteria of adjusted P value < 0.05 aimed for reduction of false positives [21, 22]. Strikingly, it has been reported that the most epigenetic mutations appear unrelated to biological function [23].

Since DNA methylation could directly block transcription by inhibiting the binding of specific transcription factors to their target sequences on the candidate gene in the result of an obvious variation in transcriptome, integration with methylation data and matched transcriptome data was conducted to screen meaningful methylation mutations and methylated DNA regions [24].

Until now, the most significant 10 methylation variations found through difference analyses were lack of studies, with the function of their methylated gene related to immunity and cell proliferation, which could be appealing to researchers for further studies [25,26,27,28,29].

Tejera et al. [30] and Wang et al. [31] had confirmed that the mRNA level of PI3 in patients with ARDS was significantly decreased, which in turn might increase the activity of neutrophil elastase indirectly for PI3 plays a central role in controlling the excessive activity of neutrophil elastase. Liu and his colleagues [32] had showed that CX3CR1 was the crucial molecule that regulates the EGFR, Src and FAK pathways, which are crucial for epithelial and endothelial cell growth. The reduction of CX3CR1 in patients with ARDS could lead to damage to epithelial and endothelial cell regeneration. However, the mechanism of down-regulated PI3 and CX3CR1 in ARDS is unclear. Our study showed that the hypermethylation on cg24476939 of PI3 and 4 hypermethylation sites of CX3CR1 might explain down-regulation of PI3 and CX3CR1 in ARDS. Similar to PI3 and CX3CR1, FYN, DUSP6, POLB and SLC3A2, PILRA, SRPK2 and RNF1 were confirmed to be related to the pathogenesis of ARDS, including inflammation or immunity, endothelial–epithelial barrier, and coagulation function [33,34,35,36]. The present study verified above variations on mRNA levels in ARDS. Moreover, methylation alterations in these DNA regions of molecules might be responsible for the mRNA level variation in ARDS, which might be potential therapeutic targets for ARDS.

In addition, the ROC analyses on screened methylation alterations indicated that the AUC on cg03341377, cg24310395, cg07830557 and cg08418670 was up to 0.99, which meant that these DNA methylation alterations were valuable targets to distinguish ARDS from healthy controls, Table 2.

There are several limitations to the present study. First, as a retrospective study of primarily publically available data, we are not able to control for demographics, infection, patient severity, or individual treatment, and there was no analysis combined with clinical data because the microarray data did not provide clinical information. Second, although 4 microarray datasets were merged, a larger sample size may have been better. The ideal method of transcriptome and genome interrelated analysis should measure DNA methylation and mRNA in the same samples and then perform a co-expression analysis. However, in the public database, there was no such completed microarray data regarding ARDS. Besides, for the limitation of bioinformatics, the further studies and experiments should be conducted to sequentially verify and research these methylation alterations screened by our analyses.

Availability of data and materials

All data generated and/or analysed during this study can be found in GEO database. The series of these studies were GSE32707, GSE66890, GSE10474, GSE61672 and GSE67530.



acute respiratory distress syndrome


gene expression omnibus


 myosin light chain kinase


supplemental material


differentially expressed genes


protein–protein interaction network


gene ontology


Kyoto Encyclopaedia of Genes and Genomes


receiver operating characteristic curves


area under the curve


peptidase inhibitor 3


C-X3-C motif chemokine receptor 1


FYN proto-oncogene, Src family tyrosine kinase


dual specificity phosphatase 6


DNA polymerase beta


solute carrier family 3 member 2


paired immunoglobin like type 2 receptor alpha


SRSF protein kinase 2


ring finger protein 1


tripartite motif containing 33


  1. 1.

    Fan E, Brodie D, Slutsky AS. Acute respiratory distress syndrome: advances in diagnosis and treatment. JAMA. 2018;319(7):698–710.

  2. 2.

    Thompson BT, Chambers RC, Liu KD. Acute respiratory distress syndrome. N Engl J Med. 2017;377(19):1904–5.

  3. 3.

    Kangelaris KN, Prakash A, Liu KD, et al. Increased expression of neutrophil-related genes in patients with early sepsis-induced ARDS. Am J Physiol Lung Cell Mol Physiol. 2015;308(11):L1102–13.

  4. 4.

    Dolinay T, Kim YS, Howrylak J, et al. Inflammasome-regulated cytokines are critical mediators of acute lung injury. Am J Respir Crit Care Med. 2012;185(11):1225–34.

  5. 5.

    Howrylak JA, Dolinay T, Lucht L, et al. Discovery of the gene signature for acute lung injury in patients with sepsis. Physiol Genomics. 2009;37(2):133–9.

  6. 6.

    Koch A, Joosten SC, Feng Z, et al. Author correction: analysis of DNA methylation in cancer: location revisited. Nat Rev Clin Oncol. 2018;15(7):467.

  7. 7.

    Vogel CA, LaSalle J. The landscape of DNA methylation amid a perfect storm of autism aetiologies. Nat Rev Neurosci. 2016;17(7):411–23.

  8. 8.

    Bansal A, Pinney SE. DNA methylation and its role in the pathogenesis of diabetes. Pediatr Diabetes. 2017;18(3):167–77.

  9. 9.

    Yang IV, Lozupone CA, Schwartz DA. The environment, epigenome, and asthma. J Allergy Clin Immunol. 2017;140(1):14–23.

  10. 10.

    Long H, Yin H, Wang L, et al. The critical role of epigenetics in systemic lupus erythematosus and autoimmunity. J Autoimmun. 2016;74:118–38.

  11. 11.

    Dhas BB, Antony HA, Bhat V, et al. Global DNA methylation in neonatal sepsis. Indian J Pediatr. 2015;82(4):340–4.

  12. 12.

    Szilagyi KL, Liu C, Zhang X, et al. Epigenetic contribution of the myosin light chain kinase gene to the risk for acute respiratory distress syndrome. Transl Res. 2017;180:12–21.

  13. 13.

    Wingo AP, Gibson G. Blood gene expression profiles suggest altered immune function associated with symptoms of generalized anxiety disorder. Brain Behav Immun. 2015;43:184–91.

  14. 14.

    Issa JJ, Roboz G, Rizzieri D, et al. Safety and tolerability of guadecitabine (SGI-110) in patients with myelodysplastic syndrome and acute myeloid leukaemia: a multicentre, randomised, dose-escalation phase 1 study. Lancet Oncol. 2015;16(9):1099–110.

  15. 15.

    Huang X, Kong G, Li Y, et al. Decitabine and 5-azacitidine both alleviate LPS induced ARDS through anti-inflammatory/antioxidant activity and protection of glycocalyx and inhibition of MAPK pathways in mice. Biomed Pharmacother. 2016;84:447–53.

  16. 16.

    Lei C, Jiao Y, He B, et al. RIP140 down-regulation alleviates acute lung injury via the inhibition of LPS-induced PPARgamma promoter methylation. Pulm Pharmacol Ther. 2016;37:57–64.

  17. 17.

    Xia L, Su X, Shen J, et al. ANLN functions as a key candidate gene in cervical cancer as determined by integrated bioinformatic analysis. Cancer Manag Res. 2018;10:663–70.

  18. 18.

    Leek JT, Johnson WE, Parker HS, et al. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28(6):882–3.

  19. 19.

    Leek JT. svaseq: removing batch effects and other unwanted noise from sequencing data. Nucleic Acids Res. 2014.

  20. 20.

    Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.

  21. 21.

    Colquhoun D. The reproducibility of research and the misinterpretation of p-values. R Soc Open Sci. 2017;4(12):171085.

  22. 22.

    Aggarwal S, Yadav AK. False discovery rate estimation in proteomics. Methods Mol Biol. 2016;1362:119–28.

  23. 23.

    Shah SP, Roth A, Goya R, et al. The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012;486(7403):395–9.

  24. 24.

    Sanford T, Meng MV, Railkar R, et al. Integrative analysis of the epigenetic basis of muscle-invasive urothelial carcinoma. Clin Epigenetics. 2018;10:19.

  25. 25.

    Gelkop S, Gish GD, Babichev Y, et al. T cell activation-induced CrkII binding to the Zap70 protein tyrosine kinase is mediated by Lck-dependent phosphorylation of Zap70 tyrosine 315. J Immunol. 2005;175(12):8123–32.

  26. 26.

    Bala S, Kumar A, Soni S, et al. Emp is a component of the nuclear matrix of mammalian cells and undergoes dynamic rearrangements during cell division. Biochem Biophys Res Commun. 2006;342(4):1040–8.

  27. 27.

    Jiang M, Gao Y, Yang T, et al. Cyclin Y, a novel membrane-associated cyclin, interacts with PFTK1. FEBS Lett. 2009;583(13):2171–8.

  28. 28.

    Bohnekamp J, Schoneberg T. Cell adhesion receptor GPR133 couples to Gs protein. J Biol Chem. 2011;286(49):41912–6.

  29. 29.

    Kuo CJ, LaMontagne KJ, Garcia-Cardena G, et al. Oligomerization-dependent regulation of motility and morphogenesis by the collagen XVIII NC1/endostatin domain. J Cell Biol. 2001;152(6):1233–46.

  30. 30.

    Tejera P, O’Mahony DS, Owen CA, et al. Functional characterization of polymorphisms in the peptidase inhibitor 3 (elafin) gene and validation of their contribution to risk of acute respiratory distress syndrome. Am J Respir Cell Mol Biol. 2014;51(2):262–72.

  31. 31.

    Wang Z, Beach D, Su L, et al. A genome-wide expression analysis in blood identifies pre-elafin as a biomarker in ARDS. Am J Respir Cell Mol Biol. 2008;38(6):724–32.

  32. 32.

    Liu P, Liang Y, Jiang L, et al. CX3CL1/fractalkine enhances prostate cancer spinal metastasis by activating the Src/FAK pathway. Int J Oncol. 2018;53(4):1544–56.

  33. 33.

    Suresh K, Servinsky L, Reyes J, et al. CD36 mediates H2O2-induced calcium influx in lung microvascular endothelial cells. Am J Physiol Lung Cell Mol Physiol. 2017;312(1):L143–53.

  34. 34.

    Ikeda K, Kinoshita M, Kayama H, et al. Slc3a2 mediates branched-chain amino-acid-dependent maintenance of regulatory T cells. Cell Rep. 2017;21(7):1824–38.

  35. 35.

    Kobayashi H, Butler JM, O’Donnell R, et al. Angiocrine factors from Akt-activated endothelial cells balance self-renewal and differentiation of haematopoietic stem cells. Nat Cell Biol. 2010;12(11):1046–56.

  36. 36.

    Vaz M, Hwang SY, Kagiampakis I, et al. Chronic cigarette smoke-induced epigenomic changes precede sensitization of bronchial epithelial cells to single-step transformation by KRAS mutations. Cancer Cell. 2017;32(3):360–76.

Download references


We would like to thank all of the microarray data contributors of this study and all of the patients and volunteers who participated in this study.


Supported in part by grants from the National Natural Science Foundation of China (Grant Number: 81571847) and the projects of Jiangsu Province’s Medical Key Discipline (ZDXKA2016025). The funding sources had no role in the design and implementation of the study; the collection, management, analysis, and interpretation of the data; or the preparation, review, or approval of the article.

Author information

SZ and HQ had full access to all of the data in the study and take responsibility for the integrity and the accuracy of the data analysis. SZ and HQ performed the data download, bioinformatic analysis, and preparation of the article for publication. Bioinformatics engineer LW checked the methodology, the bioinformatic analysis process and the results. All authors participated in writing the article and preparing the figures. All authors read and approved the final manuscript.

Correspondence to Haibo Qiu.

Ethics declarations

Ethics approval and consent to participate

This study was reviewed and approved by the Institutional Ethics Committee of Zhongda Hospital. Institutional Ethics Committee of Zhongda Hospital classified it as exemption from Review. All the participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1. The flow-process diagram of the integrative analysis for exploration of the potentially meaningful DNA methylation alterations.

Additional file 2: Figure S2. The Venn diagram for the intersection of methylation sites between Melanoderm and Caucasian. The left Venn diagram was the intersection of hypermethylation sites between Melanoderm and Caucasian. The right Venn diagram was the intersection of demethylation sites between Melanoderm and Caucasian.

Additional file 3: Figure S3. The Venn Diagram for the intersection of methylation alterations and mRNA variations. The left Venn Diagram was the intersection of the genes between demethylation sites and high expression mRNA. The right Venn Diagram for the intersection of the genes between hypermethylation sites and low expression mRNA.

Additional file 4: Table S1. Log fold change, P value and adjusted P value of differential mRNA.

Additional file 5: Figure S4. The Violin Diagram for mRNA expression of 10 screened genes.

Additional file 6: Table S2. Annotations link in UniProt database and Search strategy in Pubmed database for screened genes.

Additional file 7: Figure S5. The Bar Diagram for GO term enrichment analyses of 10 screened genes.

Additional file 8: Figure S6. The Bubble Diagram for KEGG pathway enrichment analyses of 10 screened genes.

Additional file 9: Table S3. Our study provides ARDS researchers with five valuable research points.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, S., Wu, Z., Xie, J. et al. DNA methylation exploration for ARDS: a multi-omics and multi-microarray interrelated analysis. J Transl Med 17, 345 (2019).

Download citation


  • ARDS
  • DNA methylation
  • mRNA
  • Multi-omics
  • Interrelated analysis