Retinopathy of prematurity (ROP) is a leading cause of childhood blindness. At present, the molecular mechanisms underlying ROP are still far from being clearly understood. Circular RNAs (circRNAs), a novel class of noncoding RNAs, have been reported to serve vital regulatory roles in several human diseases. However, it is still unclear how circRNAs are involved in ROP. In the present study, oxygen-induced retinopathy (OIR) murine retinal samples and paired normal tissues were chosen for high-throughput transcriptome RNA sequencing and bioinformatic analyses. As a result, a total of 236 differentially expressed circRNAs, 14 differentially expressed miRNAs, and 9,756 differentially expressed mRNAs were identified in the OIR samples. Gene ontology analysis showed that angiogenesis ranked in the top five upregulated biological processes associated with differential mRNA expression. Then, 66 co-expression pairs of circRNA-mRNA were predicted according to the mRNAs that were enriched in angiogenesis. Furthermore, coregulation prediction was separately performed to identify the differentially expressed miRNAs that targeted angiogenesis-associated circRNAs or mRNAs. Finally, nine differentially expressed circRNAs were predicted to be competing endogenous RNAs by constructing a circRNA-miRNA-mRNA network followed by reverse transcription-quantitative PCR validation. The results of the present study suggest that the identified set of circRNA transcripts and the potential regulatory mechanisms for the development of ROP are worthy of functional studies.
Keywords: RNA sequencing; angiogenesis; circular RNA; oxygen-induced retinopathy; retinopathy of prematurity.