Am J Perinatol 2022; 39(15): 1702-1710
DOI: 10.1055/s-0041-1726124
Original Article

Integrated MicroRNA–mRNA Analyses of Distinct Expression Profiles in Hyperoxia-Induced Bronchopulmonary Dysplasia in Neonatal Mice

Chengqiang Wang
1   Public Health, Guilin Medical University, Lingui, Guilin, People's Republic of China
,
Sheng Zhang
2   Affiliated BaYi Children's Hospital, Seventh Medical Center of People's Liberation Army General Hospital, Dongcheng, Beijing, People's Republic of China
3   Beijing Key Laboratory of Pediatric Organ Failure, Dongcheng, Beijing, People's Republic of China
,
Lina Zhu
2   Affiliated BaYi Children's Hospital, Seventh Medical Center of People's Liberation Army General Hospital, Dongcheng, Beijing, People's Republic of China
,
Jun Duan
4   Department of Pediatrics, the First Affiliated Hospital of Anhui Medical University, Shushan, Hefei, People's Republic of China
,
Bo Huang
1   Public Health, Guilin Medical University, Lingui, Guilin, People's Republic of China
,
1   Public Health, Guilin Medical University, Lingui, Guilin, People's Republic of China
2   Affiliated BaYi Children's Hospital, Seventh Medical Center of People's Liberation Army General Hospital, Dongcheng, Beijing, People's Republic of China
› Author Affiliations
Funding This work was supported by the project supported by the National Natural Science Foundation of China (81270059); the Science and Technology Base and Talent Project of Guangxi Province (AD19245006).

Abstract

Objective Bronchopulmonary dysplasia (BPD) is a common chronic lung disease of preterm neonates; the underlying pathogenesis is not fully understood. Recent studies suggested microRNAs (miRNAs) may be involved in BPD.

Study Design miRNA and mRNA microarrays were performed to analyze the expression profiles of miRNA and mRNA in BPD and control lung tissues after oxygen and air exposure on day 21. Bioinformatics methods, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), were performed to predict the potential functions of differentially expressed genes. Then, miRNA–mRNA regulatory network was constructed by protein–protein interaction (PPI) data and TarBase data.

Results Our results showed that a total of 192 differentially expressed miRNAs (74 downregulated and 118 upregulated) and 1,225 differentially expressed mRNAs (479 downregulated and 746 upregulated) were identified between BPD mice and normoxia-control mice. GO and KEGG analysis showed that for downregulated genes, the top significant enriched GO terms and KEGG pathways were both mainly related to immune and inflammation processes; for upregulated genes, the top significant enriched GO terms and KEGG pathways were both mainly related to extracellular matrix (ECM) remodeling. PPI network and miRNA–mRNA regulatory network construction revealed that the key genes and pathways associated with inflammation and immune regulation.

Conclusion Our findings revealed the integrated miRNA–mRNA data of distinct expression profiles in hyperoxia-induced BPD mice, and may provide some clues of the potential biomarkers for BPD, and provide novel insights into the development of new promising biomarkers for the treatment of BPD.

Key Points

  • Integrated advanced bioinformatics methods may offer a better way to understand the molecular expression profiles involved in BPD.

  • ECM remodeling, inflammation, and immune regulation may be essential to BPD.

  • The miRNA–mRNA regulatory network construction may contribute to develop new biomarkers for the treatment of BPD.

Authors' Contributions

X.Z. designed and supervised the study; C.W., S.Z., L.Z., J.D., and B.H. performed the experiments and analysis work; and X.Z. prepared the manuscript. All authors reviewed the final paper.


Supplementary Material



Publication History

Received: 10 September 2020

Accepted: 02 January 2021

Article published online:
23 March 2021

© 2021. Thieme. All rights reserved.

Thieme Medical Publishers, Inc.
333 Seventh Avenue, 18th Floor, New York, NY 10001, USA

 
  • References

  • 1 Zhang X, Peng W, Zhang S. et al. MicroRNA expression profile in hyperoxia-exposed newborn mice during the development of bronchopulmonary dysplasia. Respir Care 2011; 56 (07) 1009-1015
  • 2 Villamor-Martinez E, Álvarez-Fuente M, Ghazi AMT. et al. Association of chorioamnionitis with bronchopulmonary dysplasia among preterm infants: a systematic review, meta-analysis, and metaregression. JAMA Netw Open 2019; 2 (11) e1914611
  • 3 Nold MF, Mangan NE, Rudloff I. et al. Interleukin-1 receptor antagonist prevents murine bronchopulmonary dysplasia induced by perinatal inflammation and hyperoxia. Proc Natl Acad Sci USA 2013; 110 (35) 14384-14389
  • 4 Gronbach J, Shahzad T, Radajewski S. et al. The potentials and caveats of mesenchymal stromal cell-based therapies in the preterm infant. Stem Cells Int 2018; 2018: 9652897
  • 5 Ekiz HA, Ramstead AG, Lee SH. et al. T cell-expressed microRNA-155 reduces lifespan in a mouse model of age-related chronic inflammation. J Immunol 2020; 204 (08) 2064-2075
  • 6 Pham TT, Ban J, Lee K. et al. MicroRNA gga-miR-10a-mediated transcriptional regulation of the immune genes in necrotic enteritis afflicted chickens. Dev Comp Immunol 2020; 102: 103472
  • 7 Coarfa C, Zhang Y, Maity S. et al. Sexual dimorphism of the pulmonary transcriptome in neonatal hyperoxic lung injury: identification of angiogenesis as a key pathway. Am J Physiol Lung Cell Mol Physiol 2017; 313 (06) L991-L1005
  • 8 Butler B, De Dios R, Nguyen L, McKenna S, Ghosh S, Wright CJ. Developmentally regulated innate immune NFκB signaling mediates IL-1α expression in the perinatal murine lung. Front Immunol 2019; 10: 1555
  • 9 Shivanna B, Maity S, Zhang S. et al. Gene expression profiling identifies cell proliferation and inflammation as the predominant pathways regulated by aryl hydrocarbon receptor in primary human fetal lung cells exposed to hyperoxia. Toxicol Sci 2016; 152 (01) 155-168
  • 10 Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP. Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res 2003; 31 (04) e15
  • 11 Irizarry RA, Hobbs B, Collin F. et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003; 4 (02) 249-264
  • 12 Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003; 19 (02) 185-193
  • 13 Harris MA, Clark J, Ireland A. et al; Gene Ontology Consortium. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res 2004; 32 (Database issue): D258-D261
  • 14 Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000; 28 (01) 27-30
  • 15 Li Y, Li Q, Wang C, Li S, Yu L. Long noncoding RNA expression profile in BV2 microglial cells exposed to lipopolysaccharide. BioMed Res Int 2019; 2019: 5387407
  • 16 Wang W, Wang T, Wang Y. et al. Integration of gene expression profile data to verify hub genes of patients with Stanford A aortic dissection. BioMed Res Int 2019; 2019: 3629751
  • 17 Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNA.org resource: targets and expression. Nucleic Acids Res 2008; 36 (Database issue): D149-D153
  • 18 Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. eLife 2015; 4: 4
  • 19 Krek A, Grün D, Poy MN. et al. Combinatorial microRNA target predictions. Nat Genet 2005; 37 (05) 495-500
  • 20 Jiang S, Fang X, Liu M, Ni Y, Ma W, Zhao R. MiR-20b down-regulates intestinal ferroportin expression in vitro and in vivo. Cells 2019; 8 (10) E1135
  • 21 Kerrien S, Aranda B, Breuza L. et al. The IntAct molecular interaction database in 2012. Nucleic Acids Res 2012; 40 (Database issue): D841-D846
  • 22 Pagel P, Kovac S, Oesterheld M. et al. The MIPS mammalian protein-protein interaction database. Bioinformatics 2005; 21 (06) 832-834
  • 23 Chou CH, Shrestha S, Yang CD. et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res 2018; 46 (D1): D296-D302
  • 24 Sethupathy P, Corda B, Hatzigeorgiou AG. TarBase: a comprehensive database of experimentally supported animal microRNA targets. RNA 2006; 12 (02) 192-197
  • 25 Nepusz G.C.a.T. The igraph software package for complex network research. InterJournal Complex Systems 2006; 1695 (05) 1-9
  • 26 Saçar Demirci MD, Yousef M, Allmer J. Computational prediction of functional microRNA-mRNA interactions. Methods Mol Biol 2019; 1912: 175-196
  • 27 Shannon P, Markiel A, Ozier O. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003; 13 (11) 2498-2504
  • 28 Zhang X, Liu S, Hu T, Liu S, He Y, Sun S. Up-regulated microRNA-143 transcribed by nuclear factor kappa B enhances hepatocarcinoma metastasis by repressing fibronectin expression. Hepatology 2009; 50 (02) 490-499
  • 29 Kumar VHS, Wang H, Nielsen L. Adaptive immune responses are altered in adult mice following neonatal hyperoxia. Physiol Rep 2018;6(02):
  • 30 Balany J, Bhandari V. Understanding the impact of infection, inflammation, and their persistence in the pathogenesis of bronchopulmonary dysplasia. Front Med (Lausanne) 2015; 2: 90
  • 31 Kasat K, Patel H, Predtechenska O, Vancurova I, Davidson D. Anti-inflammatory actions of endogenous and exogenous interleukin-10 versus glucocorticoids on macrophage functions of the newly born. J Perinatol 2014; 34 (05) 380-385
  • 32 Kramer BW, Kallapur S, Newnham J, Jobe AH. Prenatal inflammation and lung development. Semin Fetal Neonatal Med 2009; 14 (01) 2-7
  • 33 Zhang X, Xu J, Wang J. et al. Reduction of microRNA-206 contributes to the development of bronchopulmonary dysplasia through up-regulation of fibronectin 1. PLoS One 2013; 8 (09) e74750
  • 34 Zhang X, Wang H, Shi Y. et al. Role of bone marrow-derived mesenchymal stem cells in the prevention of hyperoxia-induced lung injury in newborn mice. Cell Biol Int 2012; 36 (06) 589-594
  • 35 Bhaskaran M, Xi D, Wang Y. et al. Identification of microRNAs changed in the neonatal lungs in response to hyperoxia exposure. Physiol Genomics 2012; 44 (20) 970-980
  • 36 Olave N, Lal CV, Halloran B. et al. Regulation of alveolar septation by microRNA-489. Am J Physiol Lung Cell Mol Physiol 2016; 310 (05) L476-L487
  • 37 Wu YT, Chen WJ, Hsieh WS. et al. MicroRNA expression aberration associated with bronchopulmonary dysplasia in preterm infants: a preliminary study. Respir Care 2013; 58 (09) 1527-1535
  • 38 Dong J, Carey WA, Abel S. et al. MicroRNA-mRNA interactions in a murine model of hyperoxia-induced bronchopulmonary dysplasia. BMC Genomics 2012; 13: 204
  • 39 Carrera P, Di Resta C, Volonteri C. et al. Exome sequencing and pathway analysis for identification of genetic variability relevant for bronchopulmonary dysplasia (BPD) in preterm newborns: a pilot study. Clin Chim Acta 2015; 451 (pt A): 39-45
  • 40 Mao Q, Gundavarapu S, Patel C, Tsai A, Luks FI, De Paepe ME. The Fas system confers protection against alveolar disruption in hyperoxia-exposed newborn mice. Am J Respir Cell Mol Biol 2008; 39 (06) 717-729
  • 41 Ma CJ, Liu X, Che L, Liu ZH, Samartzis D, Wang HQ. Stem cell therapies for intervertebral disc degeneration: immune privilege reinforcement by Fas/FasL regulating machinery. Curr Stem Cell Res Ther 2015; 10 (04) 285-295