Int J Med Sci 2019; 16(10):1356-1365. doi:10.7150/ijms.35149

Research Paper

Identifying circRNA-associated-ceRNA networks in retinal neovascularization in mice

Manjing Cao1,2, Lusi Zhang1,2, Jiang-Hui Wang3,4, Huilan Zeng1,2, Yingqian Peng1,2, Jingling Zou1,2, Jingming Shi1,2, Liwei Zhang1,2, Yun Li1,2, Shigeo Yoshida5, Luosheng Tang1,2, Yedi Zhou1,2 Corresponding address

1. Department of Ophthalmology, The Second Xiangya Hospital, Central South University, Changsha, Hunan 410011, China
2. Hunan Clinical Research Center of Ophthalmic Disease, Changsha, Hunan 410011, China
3. Centre for Eye Research Australia, Royal Victorian Eye and Ear Hospital, East Melbourne, Victoria, Australia
4. Ophthalmology, Department of Surgery, University of Melbourne, East Melbourne, Victoria, Australia
5. Department of Ophthalmology, Kurume University School of Medicine, Kurume, Fukuoka 830-0011, Japan

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
How to cite this article:
Cao M, Zhang L, Wang JH, Zeng H, Peng Y, Zou J, Shi J, Zhang L, Li Y, Yoshida S, Tang L, Zhou Y. Identifying circRNA-associated-ceRNA networks in retinal neovascularization in mice. Int J Med Sci 2019; 16(10):1356-1365. doi:10.7150/ijms.35149. Available from http://www.medsci.org/v16p1356.htm

Abstract

Retinal neovascularization is a complication which caused human vision loss severely. It has been shown that circular RNAs (circRNAs) play essential roles in gene regulation. However, circRNA expression profile and the underlying mechanisms in retinal neovascular diseases remain unclear. In the present study, we identified altered circRNAs in the retinas of oxygen-induced retinopathy (OIR) mouse model by microarray profiling. Microarray analysis revealed that 539 circRNAs were significantly altered in OIR retinas compared with controls. Among them, 185 up-regulated and 354 down-regulated circRNAs were identified. The expression levels of 4 altered circRNAs including mmu_circRNA_002573, mmu_circRNA_011180, mmu_circRNA_016108 and mmu_circRNA_22546 were validated by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Bioinformatic analysis with validated circRNAs such as competing endogenous RNA (ceRNA) regulatory networks with Gene Ontology (GO) enrichment analysis demonstrated that qRT-PCR validated circRNAs were associated with cellular process, cell part and phosphoric ester hydrolase activity. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that MAPK signaling pathway and renin-angiotensin system were related to validated circRNAs, suggesting these pathways may participate in pathological angiogenesis. The results together suggested that circRNAs were aberrantly expressed in OIR retinas and may play potential roles in retinal neovascular diseases.

Keywords: circular RNA (circRNA), expression profile, microarray, oxygen-induced retinopathy, retinal neovascularization

Introduction

Ischemia-induced retinal neovascularization is a key pathological process of many retinal vascular diseases like diabetic retinopathy and retinopathy of prematurity (ROP). It is also a leading cause of vision loss as well as blindness in developed countries and regions [1]. Studies have shown that a variety of altered genes are involved in the pathogenesis of retinal neovascularization [2, 3]. Several molecules, such as vascular endothelial growth factor (VEGF), basic fibroblast growth factor as well as periostin, have been identified as key mediators in pathological angiogenesis [4, 5]. In particular, VEGF has been widely studied and targeted in clinical applications where reports showed anti-VEGF treatment was remarkably efficient in patients [6]. However, some patients are not responsive to anti-VEGF therapy [7], suggesting further studies should aim to look for other therapeutic targets and biomarkers.

Circular RNA (circRNA), as a type of endogenous non-coding RNA, is characterized by covalently closed continuous loop structure lacking 5' cap and poly(A) tail at 3' ends [8]. CircRNA is widely expressed in serum and tumor tissues in mammals [9, 10]. CircRNA functions as a molecular sponge efficiently targeting miRNA and inhibiting miRNA transcription [11, 12]. Through targeting miRNA, circRNA regulates downstream gene expression and may play a crucial role in disease mechanisms [13]. A few studies have demonstrated that several circRNAs participate in the pathogenesis of diabetic retinopathy [14, 15], vascular endothelial cells proliferation and angiogenesis [16], suggesting circRNAs could serve as potential biomarkers for diagnosis of diabetic retinopathy and provide novel therapeutic targets to treat diabetic retinopathy. A study reported that circRNA-MYLK served as competing endogenous RNA (ceRNA) for miR-29a, contributing to epithelial-mesenchymal transition and the bladder cancer progression by activating VEGFA/VEGFR2 pathway [17]. Besides, circRNAs ZNF280C_hsa_circ_001211 and SIAE_ hsa_circ_002083 may participate in the key pathways of ROP pathogenesis [18]. However, the involvement of other potential circRNAs and the underlying mechanisms of ROP remain unclear.

We previously demonstrated differentially expressed mRNAs and long non-coding RNAs in the oxygen-induced retinopathy (OIR) mouse model [19], suggesting that non-coding RNAs played crucial and different roles in retinal neovascularization. In this study, we performed microarray analysis aiming to profile expressions of another kind of emerging non-coding RNA, circRNA in retinas of OIR mice. Furthermore, we interrogated the putative functions of the altered circRNAs by bioinformatic analysis.

Materials and Methods

Animals and statement of ethics

C57BL/6J mice (Hunan SJA Laboratory Animal, Changsha, China) were used in the present study. The animal experiments were conducted based on ARVO Statement for the Use of Animals in Ophthalmic and Vision Research, and were subjected to approval by the Institutional Animal Care and Use Committee of Central South University.

Oxygen-induced retinopathy mouse model

The OIR mouse model was established as described [20]. Pups were exposed to hyperoxia environment (75% oxygen) for 5 days at postnatal day 7 (P7), followed by return to the environment of room air at P12. Newborn pups kept continuously in the room air environment were used as the controls. The retina samples were collected from mice of both groups at P17 for analysis.

RNA isolation

Total RNAs were isolated from retinas (retinas from both eyes of the same mice were pooled as one sample) by using TRIzol reagent (Invitrogen, Carlsbad, USA). The concentrations were assessed by the NanoDrop ND-1000 (Thermo Scientific, Wilmington, DE, USA). The RNA integrity was examined by electrophoresis on a denaturing agarose gel prior to further investigation.

Microarray analysis of circRNAs

Microarray analysis is widely used to identify altered circRNAs in in vivo studies [21]. Three pairs of retinas in each group were used for the analysis. Sample labeling and microarray hybridization were conducted according to the standard protocols Arraystar (Rockville, MD, USA). Briefly, total RNAs were digested with Rnase R (Epicentre, Madison, WI, USA). A random priming method with the Arraystar Super RNA Labeling Kit was utilized to amplify the enriched circRNAs and to transcribe the circRNAs into fluorescent cRNA. Labeled cRNAs were hybridized onto Mouse circRNA Array V2 (8x15K, Arraystar). Then the slides were washed, and the arrays were then scanned by Agilent Scanner G2505C, and the images were analyzed by an Agilent Feature Extraction software (version 11.0.1.1). Detected circRNAs were regarded as significantly differentially expressed by the value of fold change≥1.5 and P<0.05.

Validation of quantitative real-time reverse transcription-polymerase chain reaction

To validate the data of microarray analysis, quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) was performed as previously described [22]. Briefly, total RNAs were transcribed into cDNAs by the SuperScript III Reverse Transcriptase kit (Invitrogen, Carlsbad, CA, USA). The ViiA 7 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) as well as 2× PCR Master Mix were used to perform the qRT-PCR. The primers used for qRT-PCR were listed in Table 1. The relative expression levels of circRNAs were normalized to GAPDH.

Bioinformatics analysis

According to TargetScan (http://www.targetscan.org) and miRanda (http://www.microrna.org), the miRNA target prediction software (Arraystar) was used to predict circRNA-miRNA interactions. Cytoscape was used to predict a circRNA-miRNA-mRNA network to further visualize the interactions. To forecast the functional annotation of target genes, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used for further bioinformatics analysis.

Statistical Analyses

In the present study, statistical differences were assessed by Student t-test, and P<0.05 was considered as statistically significant.

 Table 1 

Primers designed for qRT-PCR validation of selected circRNAs. Tm: temperature. bp: base pair

Gene nameForward and reverse primerTm (℃)Product length (bp)
GAPDHF:5' CACTGAGCAAGAGAGGCCCTAT 3'
R:5' GCAGCGAACTTTATTGATGGTATT 3'
60144
mmu_circRNA_011180F:5' CGTGACCACCCAGGAGACT 3'
R :5' GCCATGTTGTCCACTTTCTC 3'
6065
mmu_circRNA_016108F:5' AGAGGAGAACGTGCAGATG 3'
R :5' TCAGTGGGCAATGTTTCT 3'
6087
mmu_circRNA_22546F:5' CCAAGACCGATCACATCCC 3'
R:5'GAAAAAGAGACAAGTTCTTCCTTGT 3'
6078
mmu_circRNA_002573F:5' GTGTTTACTGCCTTTGATGTGGTTT 3'
R :5' CAAAATGCGTTTCCTGTGGC 3'
6063

Results

Expression profiling of altered circRNAs in OIR retinas

CircRNA expressions between OIR retinas and room air control retinas were measured by microarray. The analysis of hierarchical cluster identified different circRNA expression levels in OIR group and control group. A box plot showing the circRNA profiles (Fig. 1A) indicated similar distributions of all included samples in both groups. The variation of these detected circRNAs between each group was assessed by the scatter plot (Fig. 1B) and the volcano plot (Fig. 1C). In addition, the hierarchical cluster analysis showed the altered circRNA expression levels as well as the classifications in OIR group and control group. Totally 539 circRNAs were significantly altered in OIR retinas compared with control retinas (fold change≥1.5, P<0.05). Of which, 185 circRNAs were significantly up-regulated, and 354 circRNAs were significantly down-regulated in the OIR retinas (Fig. 1D). The top 20 significantly up- and down-regulated circRNAs were listed in Table 2-3. In particular, mmu_circRNA_007438 and mmu_circRNA_012434 were up/down-regulated circRNAs with the most significant alterations in OIR retinas.

Validation of altered circRNAs by qRT-PCR

Four of the significantly altered circRNAs including mmu_circRNA_002573, mmu_circRNA_011180, mmu_circRNA_016108 and mmu_circRNA_22546 were validated by qRT-PCR (Fig. 2). The results showed that the expressions of mmu_circRNA_011180 and mmu_circRNA_016108 were significantly increased in OIR retinas (3.68-fold, P<0.01 and 1.69-fold, P<0.05, respectively). Likewise, mmu_circRNA_002573 and mmu_circRNA_22546 were significantly decreased in OIR retinas (1.74-fold, P<0.05 and 2.10-fold, P<0.05, respectively). The qRT-PCR showed similar trend with microarray analysis, suggesting the reliability of our circRNAs expression profile by microarray.

Prediction of circRNA-miRNA interactions

CircRNA acted as microRNA sponge associating with related miRNAs, and together they made up circRNA-miRNA axis involving in disease pathogenesis. To determine the function of the validated circRNAs (mmu_circRNA_002573, mmu_circRNA_011180, mmu_circRNA_016108 and mmu_circRNA_22546), their target miRNAs were predicted by informatics analysis based on TargetScan as well as miRanda. The top 5 miRNAs related to each validated circRNA were shown in Fig. 3A, and the predicted interaction sites of mmu_circRNA_002573 and mmu_circRNA_016108 were also displayed in Fig. 3B-C.

Competing endogenous RNA (ceRNA) regulatory networks with GO enrichment and KEGG pathway analyses

To further elucidate the underlying mechanism of altered circRNAs we identified, construction of a circRNA-miRNA-mRNA regulatory network was generated by Cytoscape (Fig. 4). The network was constructed with 236 mRNAs, 4 circRNAs as decoys, and 42 predicted miRNAs in total. The network suggested that circRNAs could indirectly regulate miRNA target genes by competitively binding to miRNA through binding sites.

The GO analysis and KEGG pathway analysis were conducted to explore the biological function of the parental genes. The GO analysis revealed that the term with most genes and the term with the highest enrichment score were both cellular process (GO:0009987) for biological process (Fig.5A). For cellular component, the term “cell part” (GO: 0044464) was with most genes, and was the most significant enriched term as well (Fig.5B). Furthermore, binding (GO:0005488) was related with most genes and the term enriched most was phosphoric ester hydrolase activity (GO:0042578) (Fig.5C) for molecular function. Moreover, the enriched pathways analyzed by KEGG were MAPK signaling pathway, long-term potentiation, renin secretion, spliceosome, pentose phosphate, proteoglycans in cancer, fluid shear stress and atherosclerosis and steroid hormone biosynthesis (Fig.5D). These pathways were related to angiogenesis and endothelial cell motility.

 Table 2 

Top 20 significantly up-regulated circRNAs by microarray analysis.

circRNAP-valueFDRFC (abs)RegulationchromstrandcircRNA_typeGeneSymbol
mmu_circRNA_0074380.0147950.2774472.446349upchr4+exonicGm20459
mmu_circRNA_319920.0415230.3163412.365270upchr18+sense overlappingZbtb7c
mmu_circRNA_293190.0043260.2546562.306132upchr16-exonicSpidr
mmu_circRNA_360650.0304440.3051432.220547upchr3-exonic6530403H02Rik
mmu_circRNA_0043550.0086330.2674042.183836upchr7-sense overlappingSerpinh1
mmu_circRNA_203320.0356040.3099202.139975upchr1-sense overlappingFn1
mmu_circRNA_326280.0259830.3006992.079943upchr19-exonicHpse2
mmu_circRNA_252770.0347380.3090172.078075upchr12+intronicPsma3
mmu_circRNA_0161080.0269690.3018022.039924upchr8-exonicKlhl2
mmu_circRNA_0111800.0263660.3006991.964496upchr8-exonicDlc1
mmu_circRNA_195340.0238470.2976761.962567upchr9+sense overlappingClasp2
mmu_circRNA_0120820.0443360.3207101.919863upchr5+sense overlappingLrrc8d
mmu_circRNA_242340.0334580.3071101.917526upchr11-exonicKansl1
mmu_circRNA_451550.0406140.3150501.917268upchr9+exonicLrrc2
mmu_circRNA_191590.0166400.2818761.909133upchr16-sense overlappingSpidr
mmu_circRNA_402990.0182010.2841261.891718upchr6+sense overlappingZnrf2
mmu_circRNA_0077840.0208400.2924101.886580upchr9+exonicElovl5
mmu_circRNA_260960.0234280.2976761.876620upchr13-exonicAK045681
mmu_circRNA_439130.0442690.3207101.876312upchr9+exonicOpcml
mmu_circRNA_368320.0318960.3053131.862768upchr4-exonicPtpn3
 Table 3 

Top 20 significantly down-regulated circRNAs by microarray analysis.

circRNAP-valueFDRFC (abs)RegulationchromstrandcircRNA_typeGeneSymbol
mmu_circRNA_0124340.0113840.2746633.135515downchr9-sense overlappingFat3
mmu_circRNA_373280.0193370.2872813.042786downchr4+exonicDhcr24
mmu_circRNA_191750.0260230.3006992.988649downchr17+intronicTulp4
mmu_circRNA_0080090.0234700.2976762.958134downchr17+intronicTulp4
mmu_circRNA_419900.0063210.2602452.945088downchr7+exonicPgm2l1
mmu_circRNA_0143930.0149780.2774472.911356downchr5-exonicZfp644
mmu_circRNA_390990.0112330.2746632.817654downchr5-intronicZfp644
mmu_circRNA_0118650.0285590.3033672.801541downchr5-exonicZfp644
mmu_circRNA_391000.0248880.2981102.736680downchr5-intronicZfp644
mmu_circRNA_259300.0432730.3192722.616756downchr12+exonicRapgef5
mmu_circRNA_222650.0130770.2758942.613595downchr10-exonicTrappc10
mmu_circRNA_301960.0238380.2976762.583449downchr17+intronicTulp4
mmu_circRNA_333630.0061700.2598462.549741downchr2-exonicStrbp
mmu_circRNA_454530.0186460.2858782.541237downchrX+exonicOcrl
mmu_circRNA_192410.0118660.2746632.511977downchr2-sense overlappingStrbp
mmu_circRNA_419910.0187480.2858782.501057downchr7+exonicPgm2l1
mmu_circRNA_0148150.0230800.2970042.482177downchr4-exonicZfyve9
mmu_circRNA_0068600.0124240.2746632.431044downchr5-exonicZfp644
mmu_circRNA_189710.0046880.2546562.424547downchr10-sense overlappingL3mbtl3
mmu_circRNA_391010.0240250.2981102.383603downchr5-intronicZfp644
 Figure 1 

Altered circRNAs between OIR retinas and room air control retinas by microarray analysis. A. Box plot. The expression profiles of circRNAs were shown as a box plot after normalization. B. Scatter plot. A scatter plot is used to show the raw variations of the expression profile between the two groups. C. Volcano plot. The horizontal line represents P=0.05, and the red points represent the statistical altered circRNAs. The green lines in B and C represent the default significant fold change (1.5). D. Hierarchical cluster analysis (heat map) represents all of the altered circRNAs between the two groups. Red and green denoted high and low expression, respectively.

Int J Med Sci Image (Click on the image to enlarge.)
 Figure 2 

The relative expressions of candidate circRNAs in OIR group and control group for validation assessed by qRT-PCR.

Int J Med Sci Image (Click on the image to enlarge.)
 Figure 3 

Detailed annotation for circRNA-miRNA interaction. A. Top 5 predicted targets of mmu_circRNA_011180, mmu_circRNA_016108, mmu_circRNA_22546 and mmu_circRNA_002573. B. Predicted interaction sites of mmu_circRNA_002573. C. Predicted interaction sites of mmu_circRNA_016108. M: be predicted by miRanda; T: be predicted by TargetScan.

Int J Med Sci Image (Click on the image to enlarge.)

Discussion

Several studies have revealed roles of non-coding RNAs, particularly miRNA and lncRNA [13, 23]. Some miRNAs, such as miR-23 and miR-27 served as promotors for neovascularization, while other miRNAs like miR-24 and miR-31 acted as inhibitors [23]. Many miRNAs have been revealed to regulate crucial genes involved in the pathogenesis of ischemic retinopathy. For example, miRNA-126 maintained the integrity of the blood-retina barrier [24] and miR-29a inhibited retinal neovascularization to prevent ROP by down-regulating AGT [25]. However, another spectrum of non-coding RNA, circRNA, was considered to be a by-product of error splicing and was overlooked for a long time. With closed-loop structures, circRNAs are more stable than linear RNAs and play essential roles in the pathogenesis of human diseases [26]. A study reported that 529 circRNAs were abnormally expressed in diabetic retinas in human [15]. Another study identified altered circRNAs in serum samples from patients with or without type 2 diabetes mellitus [14]. Systemic dissection of circRNA profile in OIR model may provide clues to roles of circRNA in pathological neovascularization.

 Figure 4 

The predicted circRNA‐miRNA‐mRNA networks. The red color and light-blue color represent miRNA and mRNA, respectively. Yellow color and green color represent up-regulated and down-regulated circRNAs, respectively. Edges with T-shape arrow represent directed relationships, while edges without arrow represent undirected relationships.

Int J Med Sci Image (Click on the image to enlarge.)

A few studies showed an intricate interplay between circRNAs and miRNAs. CircRNAs could competitively bind to miRNAs by miRNA-binding sites and the target genes were regulated as a result [27, 28]. We examined circRNA-miRNA interactions and discovered that each selected circRNA, containing at least one miRNA binding sites, was able to interact with several miRNAs (Fig. 3). Our results suggested that circRNAs participated in OIR putatively through targeted miRNA and indirectly regulated gene expression. In addition, there is a subpopulation of circRNAs named exon-intron circRNAs (EIciRNAs), crosslink to RNA polymerase II, that remain in the nucleus and can regulate gene transcription directly in cis [29]. However, the role of EIciRNAs in OIR models largely remains unclear.

The four validated circRNAs by qRT-PCR were selected to draw a whole picture of ceRNA regulatory networks. This circRNA-miRNA-mRNA network may provide clues to the regulatory pathways in OIR for the cascade-amplifying synergistic effects of circRNA-miRNA and miRNA-mRNA. Furthermore, GO enrichment analysis and KEGG pathway analysis were also performed to functionally annotate the predicted target genes. The ceRNA network revealed a novel interaction between the altered circRNAs and 236 mRNAs, and also showed that circRNA can regulate mRNAs through circRNA-miRNA binding. This network served as a shred of credible evidence that circRNA played a role in the pathogenesis of OIR by indirectly targeting certain mRNAs. GO analysis revealed that the most significant enriched term was cellular process in biological process, while cell part was the first highly ranked enriched term in cellular component and phosphoric ester hydrolase activity in molecular function.

 Figure 5 

GO and KEGG pathway analyses of validated circRNAs. For A-C, parental genes according to the values in the enrichment score (left) and gene count (right) under the themes of biological processes, cellular components and molecular functions. The x‐ and y‐axis represent the top 10 significantly enriched terms. A represents biological processes. B represents cellular components. C represents molecular functions. D. The top enriched KEGG pathways of the significantly altered circRNA parental genes.

Int J Med Sci Image (Click on the image to enlarge.)

The KEGG pathway analysis revealed several important pathways related to angiogenesis, including MAPK signaling pathway and renin-angiotensin system (RAS). A number of studies reported that angiogenesis can be restrained via downregulation of MAPK signaling pathway [30-32]. Moreover, RAS was an essential growth factor which can stimulate angiogenesis in ischemia [33] and blockage of RAS could inhibit cancer angiogenesis [34]. Altogether, GO and KEGG analysis suggested that mmu_circRNA_002573, mmu_circRNA_011180, mmu_circRNA_016108 and mmu_circRNA_22546, may take part in the process of neovascularization via different pathways.

Although we systemically profiled circRNA expression in retinas of OIR mice, limitations were existed in the study. For example, the circRNA expression profile should be interpreted with caution due to a limited number of samples. Secondly, false-negative data could be presented because the small quantity of circRNAs may below the detection threshold [35]. Moreover, it is possible that the microarray-based screening might miss out some key circRNAs involved in retinal neovascularization, as microarray assays lack the sensitivity of advanced approaches such as next-generation sequencing. Apart from the limitations, there is a promising prospect for circRNAs in diagnosis and therapies since circRNAs are expressed widely and stably. For instance, circ_0005015 promoted retinal angiogenesis by regulating migration, proliferation as well as tube formation of the endothelial cells, and was able to act as potential biomarkers for diagnosis of diabetic retinopathy [15]. A study suggested that exogenous introduction of circRNAs may stimulate the immune system for therapeutic purposes [36]. In addition, circRNA HIPK3 promoted endothelial proliferation and vascular dysfunction by inhibiting miR-30a function [37]. After intravitreal injection of AAV vectors encoding circRNA in mice, multiple retinal cell layers were observed with robust transgene expression, potentially providing novel therapeutic platforms of circRNAs [38].

In conclusion, our study demonstrated that the interactions of circRNAs and miRNAs may indirectly regulate gene expression and altered circRNA may play a role in retinal neovascularization. Further research on downstream molecular functions and their mechanisms may expose potential therapeutic targets for ischemia-induced retinal neovascularization.

Abbreviations

ceRNA: competing endogenous RNA; circRNA: circular RNA; EIciRNAs: exon-intron circRNAs; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; lncRNA: long non-coding RNA; miRNA: microRNA; OIR: oxygen-induced retinopathy; RAS: renin-angiotensin system; ROP: retinopathy of prematurity; qRT-PCR: quantitative real-time reverse transcription-polymerase chain reaction; VEGF: vascular endothelial growth factor.

Acknowledgements

This work was supported by National Natural Science Foundation of China (No. 81800855, 81800856, 81700837 and 81500746), Natural Science Foundation of Hunan Province (No. 2018JJ3765 and 2018JJ3737), and Department of Science and Technology, Hunan (No.2015TP2007).

Competing Interests

The authors have declared that no competing interest exists.

References

1. Yoshida A, Yoshida S, Ishibashi T, Inomata H. Intraocular neovascularization. Histol Histopathol. 1999;14:1287-94

2. Ishikawa K, Yoshida S, Kadota K, Nakamura T, Niiro H, Arakawa S. et al. Gene expression profile of hyperoxic and hypoxic retinas in a mouse model of oxygen-induced retinopathy. Invest Ophthalmol Vis Sci. 2010;51:4307-19

3. Ishikawa K, Yoshida S, Kobayashi Y, Zhou Y, Nakama T, Nakao S. et al. Microarray analysis of gene expression in fibrovascular membranes excised from patients with proliferative diabetic retinopathy. Invest Ophthalmol Vis Sci. 2015;56:932-46

4. Li D, Xie K, Zhang L, Yao X, Li H, Xu Q. et al. Dual blockade of vascular endothelial growth factor (VEGF) and basic fibroblast growth factor (FGF-2) exhibits potent anti-angiogenic effects. Cancer Lett. 2016;377:164-73

5. Nakama T, Yoshida S, Ishikawa K, Kubo Y, Kobayashi Y, Zhou YD. et al. Therapeutic Effect of Novel Single-Stranded RNAi Agent Targeting Periostin in Eyes with Retinal Neovascularization. Mol Ther-Nucl Acids. 2017;6:279-89

6. Osaadon P, Fagan XJ, Lifshitz T, Levy J. A review of anti-VEGF agents for proliferative diabetic retinopathy. Eye. 2014;28:510-20

7. Stitt AW, Curtis TM, Chen M, Medina RJ, McKay GJ, Jenkins A. et al. The progress in understanding and treatment of diabetic retinopathy. Prog Retin Eye Res. 2016;51:156-86

8. Chen LL, Yang L. Regulation of circRNA biogenesis. RNA Biol. 2015;12:381-8

9. Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J. et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19:141-57

10. Guo JU, Agarwal V, Guo H, Bartel DP. Expanded identification and characterization of mammalian circular RNAs. Genome Biol. 2014;15:409

11. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK. et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495:384-8

12. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A. et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333-8

13. Beermann J, Piccoli MT, Viereck J, Thum T. Non-coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol Rev. 2016;96:1297-325

14. Gu Y, Ke G, Wang L, Zhou E, Zhu K, Wei Y. Altered Expression Profile of Circular RNAs in the Serum of Patients with Diabetic Retinopathy Revealed by Microarray. Ophthalmic Res. 2017;58:176-84

15. Zhang SJ, Chen X, Li CP, Li XM, Liu C, Liu BH. et al. Identification and Characterization of Circular RNAs as a New Class of Putative Biomarkers in Diabetes Retinopathy. Investigative Ophthalmology & Visual Science. 2017;58:6500

16. Li CY, Ma L, Yu B. Circular RNA hsa_circ_0003575 regulates oxLDL induced vascular endothelial cells proliferation and angiogenesis. Biomed Pharmacother. 2017;95:1514-9

17. Zhong Z, Huang M, Lv M, He Y, Duan C, Zhang L. et al. Circular RNA MYLK as a competing endogenous RNA promotes bladder cancer progression through modulating VEGFA/VEGFR2 signaling pathway. Cancer Lett. 2017;403:305-17

18. Yang Y, Pan JJ, Zhou XG, Zhou XY, Cheng R. Differentially expressed miRNAs in premature infants with retinopathy-a bioinformatics analysis. Int J Ophthalmol. 2018;11:773-9

19. Zhang L, Fu X, Zeng H, Wang JH, Peng Y, Zhao H. et al. Microarray Analysis of Long Non-Coding RNAs and Messenger RNAs in a Mouse Model of Oxygen-Induced Retinopathy. Int J Med Sci. 2019;16:537-47

20. Connor KM, Krah NM, Dennison RJ, Aderman CM, Chen J, Guerin KI. et al. Quantification of oxygen-induced retinopathy in the mouse: a model of vessel loss, vessel regrowth and pathological angiogenesis. Nat Protoc. 2009;4:1565-73

21. Yu N, Tong Y, Zhang D, Zhao S, Fan X, Wu L. et al. Circular RNA expression profiles in hippocampus from mice with perinatal glyphosate exposure. Biochem Biophys Res Commun. 2018;501:838-45

22. Wang Z, Xu P, Chen B, Zhang Z, Zhang C, Zhan Q. et al. Identifying circRNA-associated-ceRNA networks in the hippocampus of Abeta1-42-induced Alzheimer's disease-like rats using microarray analysis. Aging (Albany NY). 2018;10:775-88

23. Zhang Y, Cai S, Jia Y, Qi C, Sun J, Zhang H. et al. Decoding Noncoding RNAs: Role of MicroRNAs and Long Noncoding RNAs in Ocular Neovascularization. Theranostics. 2017;7:3155-67

24. Bai X, Luo J, Zhang X, Han J, Wang Z, Miao J. et al. MicroRNA-126 Reduces Blood-Retina Barrier Breakdown via the Regulation of VCAM-1 and BCL2L11 in Ischemic Retinopathy. Ophthalmic Res. 2017;57:173-85

25. Chen XK, Ouyang LJ, Yin ZQ, Xia YY, Chen XR, Shi H. et al. Effects of microRNA-29a on retinopathy of prematurity by targeting AGT in a mouse model. Am J Transl Res. 2017;9:791-801

26. Verduci L, Strano S, Yarden Y, Blandino G. The circRNA-microRNA code: emerging implications for cancer diagnosis and treatment. Mol Oncol. 2019;13:669-80

27. Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505:344-52

28. Xie L, Mao M, Xiong K, Jiang B. Circular RNAs: A Novel Player in Development and Disease of the Central Nervous System. Front Cell Neurosci. 2017;11:354

29. Li Z, Huang C, Bao C, Chen L, Lin M, Wang X. et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22:256-64

30. Duan J, Hu H, Feng L, Yang X, Sun Z. Silica nanoparticles inhibit macrophage activity and angiogenesis via VEGFR2-mediated MAPK signaling pathway in zebrafish embryos. Chemosphere. 2017;183:483-90

31. Chen HX, Xu XX, Tan BZ, Zhang Z, Zhou XD. MicroRNA-29b Inhibits Angiogenesis by Targeting VEGFA through the MAPK/ERK and PI3K/Akt Signaling Pathways in Endometrial Carcinoma. Cell Physiol Biochem. 2017;41:933-46

32. Zhai KF, Duan H, Cui CY, Cao YY, Si JL, Yang HJ. et al. Liquiritin from Glycyrrhiza uralensis Attenuating Rheumatoid Arthritis via Reducing Inflammation, Suppressing Angiogenesis, and Inhibiting MAPK Signaling Pathway. J Agric Food Chem. 2019;67:2856-64

33. Sarkar C, Ganju RK, Pompili VJ, Chakroborty D. Enhanced peripheral dopamine impairs post-ischemic healing by suppressing angiotensin receptor type 1 expression in endothelial cells and inhibiting angiogenesis. Angiogenesis. 2017;20:97-107

34. Miyajima A, Kosaka T, Kikuchi E, Oya M. Renin-angiotensin system blockade: Its contribution and controversy. Int J Urol. 2015;22:721-30

35. Huang JL, Qin MC, Zhou Y, Xu ZH, Yang SM, Zhang F. et al. Comprehensive analysis of differentially expressed profiles of Alzheimer's disease associated circular RNAs in an Alzheimer's disease mouse model. Aging (Albany NY). 2018;10:253-65

36. Chen YG, Kim MV, Chen X, Batista PJ, Aoyama S, Wilusz JE. et al. Sensing Self and Foreign Circular RNAs by Intron Identity. Mol Cell. 2017;67:228-38 e5

37. Shan K, Liu C, Liu BH, Chen X, Dong R, Liu X. et al. Circular Noncoding RNA HIPK3 Mediates Retinal Vascular Dysfunction in Diabetes Mellitus. Circulation. 2017;136:1629-42

38. Meganck RM, Borchardt EK, Castellanos Rivera RM, Scalabrino ML, Wilusz JE, Marzluff WF. et al. Tissue-Dependent Expression and Translation of Circular RNAs with Recombinant AAV Vectors In Vivo. Molecular therapy Nucleic acids. 2018;13:89-98

Author contact

Corresponding address Corresponding author: Yedi Zhou, MD, PhD; Department of Ophthalmology, The Second Xiangya Hospital, Central South University, Changsha, Hunan 410011, China; Telephone: +86-731-85292175; E-mail: zhouyediedu.cn


Received 2019-3-23
Accepted 2019-7-9
Published 2019-9-20