Insights into the expression profiles and functions of circRNAs in a newborn hyperoxia-induced rat bronchopulmonary dysplasia model

Background Bronchopulmonary dysplasia (BPD) is a severe chronic lung disease in preterm infants. Circular RNAs (circRNAs) are key regulators of various biological processes. The present study aimed to explore the biological roles of circRNAs in BPD pathogenesis. Methods A newborn BPD rat model was developed to construct a circRNA library; Illumina deep sequencing (Illumina, San Diego, CA, USA) was used to reveal differential expression of circRNAs in the hyperoxia-induced BPD rat models. Sanger sequencing and a reverse transcription-polymerase chain reaction were performed to confirm circRNAs that may be related to BPD. After miRNA binding-site prediction, we constructed a network diagram of circRNA-competing endogenous RNAs (ceRNAs) related to transforming growth factor (TGF)-beta and p53 pathways using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Results In total, 256 differentially expressed circRNAs were detected between the hyperoxia group and the normoxia group. Of these circRNAs, 195 were up-regulated and 61 were down-regulated. The differences of circRNA distribution between the two groups were analyzed and six circRNAs were validated in the tissue samples. GO analysis indicated that 6519 target genes were enriched in cell location and biological processes. KEGG pathway enrichment analysis showed that circRNAs involved in 242 KEGG pathways. A network diagram of circRNA-ceRNA related to TGF-beta and p53 pathways was constructed. Conclusions CircRNAs are differentially expressed between the BPD model and control group. Many target genes of circRNAs are involved in the developmental process, which suggests that BPD may be associated with pathways including extracellular matrix-receptor interaction, vascular endothelial growth factor signaling and vascular smooth muscle contraction.