Differentially expressed miRNAs in oxygen‑induced retinopathy newborn mouse models

Transcriptome analysis of AAV-induced retinopathy models expressing human VEGF, TNF-α, and IL-6 in murine eyes

Retinopathies are multifactorial diseases with complex pathologies that eventually lead to vision loss. Animal models facilitate the understanding of the pathophysiology and identification of novel treatment options. However, each animal model reflects only specific disease aspects and understanding of the specific molecular changes in most disease models is limited. Here, we conducted transcriptome analysis of murine ocular tissue transduced with recombinant Adeno-associated viruses (AAVs) expressing either human VEGF-A, TNF-α, or IL-6. VEGF expression led to a distinct regulation of extracellular matrix (ECM)-associated genes. In contrast, both TNF-α and IL-6 led to more comparable gene expression changes in interleukin signaling, and the complement cascade, with TNF-α-induced changes being more pronounced. Furthermore, integration of single cell RNA-Sequencing data suggested an increase of endothelial cell-specific marker genes by VEGF, while TNF-α expression increased the expression T-cell markers. Both TNF-α and IL-6 expression led to an increase in macrophage markers. Finally, transcriptomic changes in AAV-VEGF treated mice largely overlapped with gene expression changes observed in the oxygen-induced retinopathy model, especially regarding ECM components and endothelial cell-specific gene expression. Altogether, our study represents a valuable investigation of gene expression changes induced by VEGF, TNF-α, and IL-6 and will aid researchers in selecting appropriate animal models for retinopathies based on their agreement with the human pathophysiology.

MicroRNA Expression Analysis of Mice Retinas with Oxygen-Induced Retinopathy by RNA Sequencing

Purpose. To characterize the microRNA (miRNA) expression profiles in the retinas of mice with oxygen-induced retinopathy by RNA sequencing and to ascertain miRNAs associated with retinal neovascularization. Methods. Retina samples were obtained from 3 groups (6 retinas/group) of OIR mice and normal mice at P17. RNA was isolated from 24 retina samples and then detected on an Illumina HiSeq. Twelve retina samples were used for quantitative polymerase chain reaction to validate the RNA sequencing. Bioinformatics analyses were performed. Result. The RNA sequence showed that 565 miRNAs were detected in the retina of OIR mice and 583 miRNAs in the retina of normal control mice. A total of 553 miRNAs were expressed in both groups. Thirty-eight miRNAs showed altered expression in both groups ( $p\le 0.05$ ). Compared with the control group, 2 miRNAs were significantly upregulated in the OIR group, while 36 miRNAs were significantly downregulated. Meanwhile, 2 candidate miRNAs (miR-181a-5p and miR-21a-5p) with significant differences in miRNA expression ( $p<0.01$ ) were selected for validation. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to confirm the relative expression of the two miRNAs. Bioinformatics analyses showed that pathways involved in ischemic retinopathy (such as TGF-β, Ras, Hippo, PI3K-Akt, VEGF, and HIF-1 signaling pathways) were enriched. Conclusions. Our study provided an overall view of miRNA profiling in the OIR retina. These miRNA profiles provide a valuable framework for the potential therapy of retinal angiogenesis.

MicroRNAs and the HIF/VEGF axis in ocular neovascular diseases

Ocular neovascular diseases, such as proliferative diabetic retinopathy, retinopathy of prematurity and neovascular age-related macular degeneration, are the leading causes of visual impairment worldwide. The hypoxia-inducible factors and vascular endothelial growth factors are key molecular promoters of ocular neovascularization. Moreover, the role of microRNAs as regulators of angiogenesis has been expanding, particularly hypoxia-associated microRNA; hypoxamiRs. This review provides a summary of hypoxamiRs that directly and specifically target HIF1A and VEGF mRNAs, thus critically involved in the regulation of ocular neovascular pathologies. The discussed microRNAs highlight putative diagnostic markers and therapeutic agents in choroidal and retinal angiogenic diseases, including proliferative diabetic retinopathy, retinopathy of prematurity and neovascular age-related macular degeneration.

Maternally expressed gene 3 regulates retinal neovascularization in retinopathy of prematurity

The mouse model of oxygen induced retinopathy is suitable for the study of various retinal neovascularization diseases, including retinopathy of prematurity. The maternally expressed gene 3 (MEG3) has been demonstrated to have an inhibitory effect on diabetic retinopathy. In this study, we investigated the role of MEG3 overexpression in oxygen-induced retinopathy in mice. The results showed that MEG3 overexpression effectively inhibited the production of retinal neovascularization in oxygen-induced retinopathy mice. It acts by down-regulating the expression of phosphoinositide 3-kinase, serine/threonine kinase, and vascular endothelial growth factor and pro-inflammatory factors. MEG3 overexpression lentivirus has a future as a new method for the clinical treatment of retinopathy of prematurity. The animal experiments were approved by the Animal Ethics Committee of Shengjing Hospital of China Medical University, China (approval No. 2016PS074K) on February 25, 2016.

Competing endogenous RNA network associated with oxygen-induced retinopathy: Expression of the network and identification of the MALAT1/miR-124-3p/EGR1 regulatory axis

Retinopathy of prematurity (ROP) is a severe retinal dysfunction in prematurely born babies. The relationship between non-coding RNAs and retinopathy of prematurity (ROP) remain unclear. Microarray analysis of lncRNAs, miRNAs, and mRNAs was conducted in a mouse model of ROP. A competing endogenous RNA (ceRNA) network was constructed. The relationship among MALAT1, miR-124-3p, and Early growth response protein 1 (EGR1) was assessed in hypoxia-induced primary human umbilical vein endothelial cells (HUVECs) and ROP mouse model. In the study, we found 2252 lncRNAs, 1239 mRNAs, and 36 miRNAs were differentially regulated. ceRNA network consisting of 21 lncRNAs, 10 miRNAs, and 19 mRNAs was established. Of the most down-regulated miRNAs, miR-124-3p was selected for additional study. miR-124-3p ceased the migration and proliferation of primary HUVECs in hypoxic conditions, and directly suppressed EGR1. Additionally, MALAT1 directly sponged miR-124-3p. Knockdown of MALAT1 decreased EGR1 expression and inhibited the migration and proliferation of primary HUVECs in hypoxia. Furthermore, these changes were rescued by depletion of miR-124-3p. In vivo, intravitreal injection of miR-124-3p, shMALAT1 decreased EGR1 expression and markedly suppressed retinal neovascularization in OIR models. Intravitreal injection of shMALAT1 and miR-124-3p antagomir at the same time can promote retinal neovascularization, which reversed the suppression of retinal neovascularization functioned by shMALAT1. In conclusion, the expression profiles of lncRNAs and miRNAs and the ceRNA network in a mouse model of ROP may be indicative of the underlying mechanisms of retinal angiogenesis and neural activity. The MALAT1/miR-124-3p/EGR1 regulatory axis is partly responsible for retinal neovascularization, which may provide a novel theoretical basis for the pathogenesis of ROP.

Cluster Analysis of Early Postnatal Biochemical Markers May Predict Development of Retinopathy of Prematurity

Purpose

Growth factors and inflammatory and angiogenetic proteins are involved in the development of retinopathy of prematurity (ROP). However, no early biochemical markers are in clinical use to predict ROP. By performing cluster analysis of multiple biomarkers, we aimed to determine patient groups with high and low risk for developing ROP.

Methods

In total, 202 protein markers in plasma were quantified by proximity extension assay from 35 extremely preterm infants on day 2 of life. Infants were sorted in groups by automated two-dimensional hierarchical clustering of all biomarkers. ROP was classified as stages I to III with or without surgical treatment. Predictive biomarkers were evaluated by analysis of variance and detected differences by two-sided paired t-test with Bonferroni corrections for multiple comparisons.

Results

Differences in 39 biochemical markers divided infants without ROP into two control groups (control 1, n = 7; control 2, n = 5; P < 0.05). Sixty-six biochemical markers defined differences between the control groups (n = 13) and all ROP infants (n = 23; P < 0.05). PARK7, VIM, MPO, CD69, and NEMO were markedly increased in control 1 compared to all ROP infants (P < 0.001). Lower TNFRSF4 and higher HER2 and GAL appeared in infants with ROP as compared to control 1 and/or 2 (P < 0.05, respectively).

Conclusions

Our data suggest that early elevated levels of PARK7, VIM, MPO, CD69, and NEMO may be associated with lower risk of developing ROP. Lower levels of TNFRSF4 with higher levels of HER2 and GAL may predict ROP development.

Translational Relevance

Cluster analysis of early postnatal biomarkers may help to identify infants with low or high risk of developing ROP.

MicroRNA-96 Promotes Vascular Repair in Oxygen-Induced Retinopathy-A Novel Uncovered Vasoprotective Function

Background and Aims

Vascular degeneration is a hallmark in the pathogenesis of oxygen-induced retinopathy (OIR). Dysregulation of microRNAs (miRNAs), key regulators of genes expressions, has been implicated in the regulation of ocular angiogenesis. However, miRNAs specific functions in impaired vascular development during OIR are poorly understood. Herein, we identified miR-96 as one of the most highly expressed miRNAs in the retina and choroid during vascular development and investigated the potential role of miR-96 on microvascular degeneration in a rat OIR model.

Methods and Results

Next generation sequencing (NGS) and qRT-PCR analysis showed that miR-96 maintain high levels of expression during ocular vascular development. Nevertheless, miR-96 was significantly downregulated in the retina and choroid of OIR rats (80% O₂ from P5 to P10) during the phase of microvascular degeneration. Similarly, human retinal microvascular endothelial cells (HRMEC) subjected to hyperoxia (80% O₂) showed a significant downregulation of miR-96 evaluated by qPCR. Interestingly, HRMEC supplemented with miR-96 regulated positively the expression of several key angiogenic factors including VEGF and ANG-2. To explore the angiogenic activity of miR-96 on HRMEC, we performed a gain/loss of function study. In a similar way to hyperoxia exposure, we observed a robust angiogenic impairment (tubulogenesis and migration) on HRMEC transfected with an antagomiR-96. Conversely, overexpression of miR-96 stimulated the angiogenic activity of HRMEC and protected against hyperoxia-induced endothelial dysfunction. Finally, we evaluated the potential vasoprotective function of miR-96 in OIR animals. Rat pups intravitreally supplemented with miR-96 mimic (1 mg/kg) displayed a significant preservation of retinal/choroidal microvessels at P10 compared to controls. This result was consistent with the maintenance of physiologic levels of VEGF and ANG-2 in the OIR retina.

Conclusion

This study demonstrates that miR-96 regulates the expression of angiogenic factors (VEGF/ANG-2) associated to the maintenance of retinal and choroidal microvasculature during physiological and pathological conditions. Intravitreal supplementation of miR-96 mimic could constitute a novel therapeutic strategy to improve vascular repair in OIR and other ischemic retinopathies.

Oxygen-induced circRNA profiles and coregulatory networks in a retinopathy of prematurity mouse model

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.

MicroRNA expression profile in retina and choroid in oxygen-induced retinopathy model

Background

Ischemic retinopathies (IRs) are leading causes of visual impairment. They are characterized by an initial phase of microvascular degeneration and a second phase of aberrant pre-retinal neovascularization (NV). microRNAs (miRNAs) regulate gene expression, and a number play a role in normal and pathological NV. But, post-transcriptional modulation of miRNAs in the eye during the development of IRs has not been systematically evaluated.

Aims & methods

Using Next Generation Sequencing (NGS) we profiled miRNA expression in the retina and choroid during vasodegenerative and NV phases of oxygen-induced retinopathy (OIR).

Results

Approximately 20% of total miRNAs exhibited altered expression (up- or down-regulation); 6% of miRNA were found highly expressed in retina and choroid of rats subjected to OIR. During OIR-induced vessel degeneration phase, miR-199a-3p, -199a-5p, -1b, -126a-3p displayed a robust decreased expression (> 85%) in the retina. While in the choroid, miR-152-3p, -142-3p, -148a-3p, -532-3p were upregulated (>200%) and miR-96-5p, -124-3p, -9a-3p, -190b-5p, -181a-1-3p, -9a-5p, -183-5p were downregulated (>70%) compared to controls. During peak pathological NV, miR-30a-5p, -30e-5p and 190b-5p were markedly reduced (>70%), and miR-30e-3p, miR-335, -30b-5p strongly augmented (by up to 300%) in the retina. Whereas in choroid, miR-let-7f-5p, miR-126a-5p and miR-101a-3p were downregulated by (>81%), and miR-125a-5p, let-7e-5p and let-7g-5p were upregulated by (>570%) during NV. Changes in miRNA observed using NGS were validated using qRT-PCR for the 24 most modulated miRNAs. In silico approach to predict miRNA target genes (using algorithms of miRSystem database) identified potential new target genes with pro-inflammatory, apoptotic and angiogenic properties.

Conclusion

The present study is the first comprehensive description of retinal/choroidal miRNAs profiling in OIR (using NGS technology). Our results provide a valuable framework for the characterization and possible therapeutic potential of specific miRNAs involved in ocular IR-triggered inflammation, angiogenesis and degeneration.

Identification of altered microRNAs in retinas of mice with oxygen-induced retinopathy

AIM

To identify disease-related miRNAs in retinas of mice with oxygen-induced retinopathy (OIR), and to explore their potential roles in retinal pathological neovascularization.

METHODS

The retinal miRNA expression profile in mice with OIR and room air controls at postnatal day 17 (P17) were determined through miRNA microarray analysis. Several miRNAs were significantly up- and down-regulated in retinas of mice with OIR compared to controls by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Two databases including Targetscan7.1 and MirdbV5 were used to predict target genes that associated with those significantly altered miRNAs in retinas of mice with OIR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were also conducted to identify possible biological functions of the target genes.

RESULTS

In comparison with room air controls, 3 and 8 miRNAs were significantly up- and down-regulated, respectively, in retinas of mice with OIR. The qRT-PCR data confirmed that mmu-miR-350-3p and mmu-miR-202-3p were significantly up-regulated, while mmu-miR-711 and mmu-miR-30c-1-3p were significantly down-regulated in mice with OIR compared to controls. GO analysis demonstrated that the identified target genes were related to functions such as cellular macromolecule metabolic process. KEGG pathway analysis showed a group of pathways, such as Wnt signaling pathway, transcriptional misregulation in cancer, Mucin type O-glycan biosynthesis, and mitogen-activated protein kinase (MAPK) signaling pathway might be involved in pathological process of retinal neovascularization.

CONCLUSION

Our findings suggest that the differentially expressed miRNAs in retinas of mice with OIR might provide potential therapeutic targets for treating retinal neovascularization.

Differentially expressed miRNAs in premature infants with retinopathy-a bioinformatics analysis

AIM

To reveal the role of miRNAs in retinopathy of prematurity (ROP) by bioinformatics analysis.

METHODS

The raw data of this study came from the researches of Wang et al and Zhao et al who analyzed the microRNA (miRNA) expression profile between ROP and controls. Based on the identified differentially expressed miRNAs, the related target genes, lncRNA and circRNA were predicted. Then we performed functional enrichment analysis to further analyze the functions of target genes.

RESULTS

Hsa-miRNA-128-3p and hsa-miRNA-9-5p showed significantly different expression in both studies. LncRNA of POLDIP2, GAS5, NEFL and UHRF1, circRNA of ZNF280C_hsa_circ_001211 and SIAE_hsa_circ_002083, tar-get gene of QKI showed meaningful differential expression in ROP. Enrichment analysis showed that TGF-β signaling pathway, PI3K-Akt signaling pathway and MAPK signaling pathway might play important roles in the prog-ress of ROP.

CONCLUSION

This research may provide a comprehensive bioinformatics analysis of differentially expressed miRNAs which are possibly involved in ROP.