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

Chengqiang Wang; Sheng Zhang; Lina Zhu; Jun Duan; Bo Huang; Xiaoying Zhang

doi:10.1055/s-0041-1726124

American Journal of Perinatology, Table of Contents

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

¹Public Health, Guilin Medical University, Lingui, Guilin, People's Republic of China

,

Sheng Zhang

²Affiliated BaYi Children's Hospital, Seventh Medical Center of People's Liberation Army General Hospital, Dongcheng, Beijing, People's Republic of China

³Beijing Key Laboratory of Pediatric Organ Failure, Dongcheng, Beijing, People's Republic of China

,

Lina Zhu

²Affiliated BaYi Children's Hospital, Seventh Medical Center of People's Liberation Army General Hospital, Dongcheng, Beijing, People's Republic of China

,

Jun Duan

⁴Department of Pediatrics, the First Affiliated Hospital of Anhui Medical University, Shushan, Hefei, People's Republic of China

,

Bo Huang

¹Public Health, Guilin Medical University, Lingui, Guilin, People's Republic of China

,

Xiaoying Zhang

¹Public Health, Guilin Medical University, Lingui, Guilin, People's Republic of China

²Affiliated BaYi Children's Hospital, Seventh Medical Center of People's Liberation Army General Hospital, Dongcheng, Beijing, People's Republic of China

› Author Affiliations

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.

Keywords

bronchopulmonary dysplasia - preterm infant - hyperoxia - microRNA - bioinformatics - immune regulation

Full Text

References

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

Supplementary Material

Supplementary Material