Abundant Neural circRNA Cdr1as Is Not Indispensable for Retina Maintenance

Introduction

Circular RNAs (circRNAs) are a class of stable, covalently closed RNA molecules produced from precursor mRNAs (pre-mRNAs) through back-splicing reactions (Jeck et al., 2013; Lasda and Parker, 2014; Ebbesen et al., 2017). Recently, a large number of circRNAs with complex tissue- and stage-specific expression patterns have been identified in eukaryotes (Salzman et al., 2013; Guo et al., 2014; Veno et al., 2015). Recent studies have demonstrated the biological functions of circRNAs at the molecular level (Chen, 2016; Li et al., 2018), including sequestration of microRNAs or associated proteins (Hansen et al., 2013; Li et al., 2017), splicing interference of their linear cognates (Ashwal-Fluss et al., 2014), modulation of transcription of parent genes (Li et al., 2015), and translation to produce peptides (Wang and Wang, 2015). However, considering the low efficiency of back-splicing and special structure of circRNAs, the function of most individual circRNAs remains elusive.

Circular RNAs are reported to be highly enriched in the central nervous system (CNS) and regulate synaptic function (Rybak-Wolf et al., 2015; You et al., 2015). As the most extensively characterized circRNA, Cdr1as is highly abundant in neurons and acts as a post-transcriptional regulator with many conserved binding sites for miR-7 and miR-671 (Hansen et al., 2013; Memczak et al., 2013). Expression of human Cdr1as in zebrafish caused midbrain defects, similar to miR-7 knockdown, and Cdr1as was once considered to possess important regulatory function (Memczak et al., 2013). However, Cdr1as knockout (KO) mouse, the first and sole circRNA KO mouse model, displayed a mild neuropsychiatric phenotype reflected in sensorimotor gating deficit and dysfunctional synaptic transmission (Piwecka et al., 2017).

The retina is a key part of the CNS and is responsible for vision production. As the abundant circRNA in the retina, Cdr1as is upregulated during retinal development (Chen et al., 2020), indicating the regulatory potential of Cdr1as in the retina. Whether and how Cdr1as affects retinal function remains unknown. To address this question, we generated a Cdr1as KO model using the CRISPR/Cas9 strategy. We found that Cdr1as deletion in mice caused mild alterations both in retinal phenotypic and vision functions as late as P300. In addition, miR-7, its target genes, and immediate early genes (IEGs) were deregulated in the Cdr1as KO retina. Our results demonstrated for the first time that Cdr1as ablation led to a mild retinal consequence in mice, indicating that Cdr1as abundance is not indispensable for retinal development and maintenance.

Results

Spatial Expression Pattern of Cdr1as in Retina and Other Tissues

It has been reported that Cdr1as is highly abundant in neurons, but is scarce in other tissues (Piwecka et al., 2017). To validate its spatial expression pattern, we investigated the expression of Cdr1as in different tissues by quantitative RT-PCR (qRT-PCR). As expected, Cdr1as was highly expressed in the brain and retina, and was expressed at low levels in other tissues, such as the lung, heart, liver, spleen, kidney, and muscle (Supplementary Figure 1A). To evaluate the localization of Cdr1as, we performed in situ hybridization (ISH) in the mouse retina at P120 using the BaseScope assay (Supplementary Figure 1B). We found Cdr1as located in the inner retina, predominantly in the inner nuclear layer (INL) near the inner plexiform layer (IPL). ISH and immunofluorescence (IF) staining demonstrated a portion of Cdr1as expressed in TH+ amacrine cells, which could release dopamine neurotransmitter (Supplementary Figure 1C). Together with our previous findings that Cdr1as increased sharply from P1 to P14 (Supplementary Figure 1D), these results indicated the potential role of Cdr1as during retinal development. However, we could not find any obvious abnormal phenotype during the retinal development, these may due to the compensation mechanisms in organisms.

Successful Knockout of Cdr1as in Mice

To investigate the role of Cdr1as in the retina, we attempted to generate KO mice by using CRISPR/Cas9. Because the linear transcript of Cdr1as cannot be detected, it is feasible to explore the function of Cdr1as using a KO mouse model without concerning the phenotype caused by linear transcript interruption. In this study, we used the CRISPR/Cas9 strategy to generate the Cdr1as KO mouse model. Two sgRNAs were designed to bind upstream and downstream of Cdr1as splice sites (Figure 1A and Supplementary Table 1). The F1/R1 and F1/R5 primer pairs capable of distinguishing homozygotes from the heterozygotes and wild-type (WT) were designed for genotyping (Figure 1B and Supplementary Table 2). The head-to-tail junction sequence was confirmed by Sanger sequencing (Supplementary Figure 2). qRT-PCR experiments showed that Cdr1as disappeared in retinal and brain tissues of KO mice (Figure 1C). Additionally, ISH further validated the successful deletion of Cdr1as in the retina of KO mice (Figure 1D). The overall survival and life span of mice revealed no difference between Cdr1as KO and WT mice.

FIGURE 1

Figure 1. Generation and identification of Cdr1as KO mice. (A) Cdr1as was removing using CRISPR/Cas9 strategy, and two sgRNAs targeting different regions of Cdr1as were designed. PAM sequences were labeled in red. (B) Genotyping strategy for the WT, Hete, and Homo strains using F1/R1 and F1/R5, which are shown in (A). (C) QPCR analysis of Cdr1as expression in the retinas and brains of WT and Cdr1as KO mice at P30. (D) BaseScope assay showing the expression of Cdr1as in WT and Cdr1as KO retinas at P30; red dots show the sense of ISH image, which has been labeled by black arrowhead. ONL, outer nuclear layer; OPL, outer plexiform layer, INL, inner nuclear layer; IPL, inner plexiform layer; GCL, ganglion cell layer.

Cdr1as Ablation Does Not Alter Retinal Structure

To assess the retinal structure in Cdr1as KO mice, fundus photography and high-resolution spectral-domain optical coherence tomography (SD-OCT) were performed to detect the retinal morphologies and organization in adult and aged mice. Fundus photographs showed normal fundus appearances in both WT and KO mice at P70, P360, and P500 (Figure 2A and Supplementary Figure 3A). In vivo OCT imaging displayed that the inner retinal layer was unaltered obviously, even though a few changes happened at certain locations in Cdr1as KO mice (Figure 2B and Supplementary Figures 3B–D). We further examined the photoreceptor, retinal neurons, and synaptic structure by IF staining in adult mice. Photoreceptor cells, immunolabeled for recoverin and cone-arrestin, showed a preserved pattern in Cdr1as KO as in WT retinas. Consistent with this, immunostaining for calbindin, pkcα, pax6, and rbpms to label horizontal cells, bipolar cells, amacrine cells, and ganglion cells were also unchanged with the Cdr1as deletion. Additionally, immunostaining patterns for synaptic vGlut1, α-synuclein, and ctbp2 demonstrated similar patterns in the plexiform layer (Figure 2C and Supplementary Figures 4A–C).

FIGURE 2

Figure 2. Structural phenotypes of retina in the adult Cdr1as KO mice. (A) Fundus photographs of WT and Cdr1as KO retinas in adult (n = 4) (left), OCT examination of WT and Cdr1as KO retinas at P70 (right). (B) Quantification of the thickness of the INL and IPL+RGC from the SD-OCT image in (A) (n = 5 for WT, n = 4 for KO). (C) Immunostaining of retinal neurons and synapses in WT and Cdr1as retinas (n = 3, scale bar: 25 μm).

Cdr1as Knockout Resulted in Slight Changes of Retinal Function

To investigate the impact of Cdr1as deletion on visional function, we first monitored electroretinography (ERG) responses in adult and aged mice. The b wave amplitudes of both scotopic and photopic ERG response were unchanged at P150 and P300 (Figures 3A,B). However, ERG results of elder mice (P500) showed that the photopic b wave amplitude increased significantly (Figure 3C). This phenotype is similar to the effects mediated by dopamine D2 receptor (D2R) KO in mice (Lavoie et al., 2014; Tian et al., 2015) and D2R antagonist in goldfish (Kim and Jung, 2012) and cats (Schneider and Zrenner, 1991). Indicating the function of dopamine impaired in Cdr1as KO mice. Due to partial Cdr1as expressed in certain amacrine cells, which could release dopamine. Cdr1as deletion results in abnormal function of amacrine and then further affects dopamine function. The major ERG wave components, oscillatory potentials (OPS), which are displayed by certain amacrine cells in the retina were also analyzed at P240. As a result, there was no significantly difference between WT and KO mice, even though the slight lower amplitude of OPS was observed in Cdr1as KO mice (Supplementary Figure 3E). Considering that Cdr1as is expressed in the inner retina, the visual evoked potential (VEP) was applied to assess the electric conduction from the retina to visual cortices at P240 and P500. The amplitudes from the N1 to P1 peak in Cdr1as KO mice were slightly decreased, but no difference was observed in two groups (Supplementary Figure 3F).

FIGURE 3

Figure 3. Assessment of retinal function in the adult and aged Cdr1as KO mice. (A) b-Wave amplitudes of scotopic and photopic for the Cdr1as KO and WT retinas at P150 (n = 5 for WT, n = 6 for KO). (B) b-Wave amplitudes of scotopic and photopic for the Cdr1as KO and WT retinas at P300 (n = 13 for WT, n = 13 for KO). (C) b-Wave amplitudes of scotopic and photopic for the Cdr1as KO and WT retinas at P500 (n = 10 for WT, n = 10 for KO). (D) OMR examination of Cdr1as KO and WT mice at P360 (n = 3 for WT, n = 4 for KO) and P500 (n = 7 for WT, n = 9 for KO). The normalized values represent means ± SEM. *P < 0.05, **P < 0.01, Student’s t-test.

Additionally, an optomotor response (OMR) assay was performed to assess the animal photopic response to a visual stimulus at P360 and P500. We found that both Cdr1as KO and WT mice responded to the rotating grating through head tracking, the vision contrasts sensitivity showed no significant alteration at P360 but decreased significantly at P500 (Figure 3D). The reduced vision acuity may be due to the reduction of dopamine in the retina, which has been reported previously (Witkovsky, 2004). Taken together, these results indicated that the deletion of Cdr1as has a very mild effect on visual function.

Retinal miRNAs and Their Target Genes Are Deregulated in Cdr1as KO Retinas

Following the assessment of the retinal structure and function in Cdr1as KO mice, molecular phenotypes were further investigated. Considering the strong sponge efficiency of Cdr1as for miRNAs, we evaluated the expression pattern of retinal miRNAs in Cdr1as KO mice. Whole-retinal small RNA-seq from adult mice showed that the expression of 25 miRNAs, such as miR-7, miR-344c, miR-322, miR-326, miR-122, and miR-7048 was significantly changed (fold change > 1.5, p < 0.05) between the Cdr1as KO and WT groups (Figure 4A), miR-7, the highly expressed miRNA, was confirmed by qRT-PCR and was downregulated in Cdr1as KO retinas (Figure 4B). The miR-7 target genes we tested, including Nr4a3, Klf4, Irs2, α-synuclein, and c-fos, displayed an increasing expression tendency (Figure 4C). Among the target genes, Nr4a3 and c-fos are IEGs, which have been linked to neural activity and could respond to different stimuli (Caputto and Guido, 2000; Piwecka et al., 2017). Other IEGs, such as Egr1, Egr4, and Arc, were also upregulated in KO retinas (Supplementary Figure 5). These results suggested that degeneration in Cdr1as KO retinas was more likely to occur under certain stimuli.

FIGURE 4

Figure 4. miRNAs and genes expression changes in Cdr1as KO retina. (A) Differentially expressed genes between Cdr1as and WT retinas in adult mouse (n = 3). The red dots represent the upregulated genes, and green dots represent the downregulated genes. (B) miR-7 expression changes in Cdr1as and WT retinas (n = 3). (C) Changes in the expression of target genes of miR-7 in Cdr1as and WT retinas (n = 3). The normalized values represent means ± SEM. *P < 0.05, **P < 0.01, ***P < 0.001, Student’s t-test.

In addition, we performed the GO and KEGG to analyze the 1686 target genes of 25 miRNAs that changed significantly. As was shown in Figure 5A, several GO terms were found to be significantly enriched, including nervous system development, regulation of axon regeneration, and neuron projection development, suggesting that some biological processes of the neuron system may be affected in the KO mice. The subsequent KEGG pathway analysis showed that PI3K-Akt, MAPK, Hif-1, Wnt, mTOR, and other important signaling pathways were significantly enriched (Figure 5B).

FIGURE 5

Figure 5. Bioinformatics analysis of the target genes of 25 miRNAs that changed significantly in Cdr1as KO retina using GO (A) and KEGG pathway analysis (B).

Discussion

In this study, we successfully generated a Cdr1as KO mouse model and first explored the function of Cdr1as in the retina. Despite of the seemingly normal appearance, Cdr1as KO mice had mild retinal phenotypes. We found that miR-7, which is highly sponged by Cdr1as, was downregulated and IEGs were upregulated (Supplementary Figure 5). These data suggested that Cdr1as KO retinas might be more vulnerable to degenerative alterations. The function of Cdr1as retina could be studied under different stimuli, such as light and circadian regulation.

The localization of Cdr1as was partially similar to that of α-synuclein (Surguchov et al., 2001; Chandra et al., 2005), which is aggregated in Parkinson’s disease retinas. miR-7 was reported to repress α-synuclein expression and protect cells against oxidative stress (Junn et al., 2009). We also found phospho-α-synuclein increased in Cdr1as KO retina (Supplementary Figure 6). Whether α-synuclein aggregation induced by miR-7 downregulation could cause the retina to be more vulnerable to stress needs further study.

Among the downregulated miRNAs in Cdr1as KO retinas, miR-326, miR-322, miR-344, miR-7048, and miR-122 should be mentioned here. CDC42 which is a member of the Rho GTPase family is a target of miR-326, and CDC42 dysregulation is linked to neuronal diseases such as Alzheimer’s and Parkinson’s disease (Zhu et al., 2000; Mitchell et al., 2007). CDC42 was upregulated in hippocampal neurons of Alzheimer’s patients (Zhu et al., 2000). Abnormal CDC42 may lead to altered retinal function. RNA helicase Ddx3x was found to be the direct target of miR-322 (Che et al., 2019). Ddx3x was reported to be essential for the even distribution of cells across layers, and Ddx3x was shown to be involved in the network that could control miR-183C accumulation in the photoreceptor layer architecture (Krol et al., 2015). Dysregulation of Ddx3x may contribute to the abnormalities in retinal function. miR-344 directly targets glycogen synthase kinase 3 beta (GSK3β) (Chen et al., 2014), which was reported to regulate Wnt signaling via phosphorylation of β-catenin (Cohen and Frame, 2001). Additionally, inactivation of GSK3β led to β-catenin stabilization and MG proliferation without retinal injury (Yao et al., 2016). miR-7048 and its target Ascl1 were involved in enhanced neuronal regeneration after injury (Lisi et al., 2017). miR-122 was downregulated in canine models of retinal degeneration (Kok et al., 2015). miR-122 was also reported to be downregulated in extracellular vesicles from ARPE19 cells under oxidative conditions (Genini et al., 2014). We speculate that dysregulated miRNA-induced target genes alteration could have some effects on the retina. We also predicted the targets of 25 miRNA that changed significantly using seven tools and further performed bioinformatics analysis based on these 1686 target genes. GO and KEGG pathway analysis showed that these genes were related to neuron system and several important signaling pathway. These results broadly reflected the biological effect of Cdr1as in the retina.

There are several unknowns in the present study. First, since the location of Cdr1as is mainly in the inner retina while the ERG response is derived from photoreceptors to bipolar and amacrine cells, it seems that ERG is not a reasonable method to assess the retinal function of Cdr1as KO retinas. In contrast, the VEP results in our study seemed to be more believable compared to the ERG data. Second, we do not know if any ultrastructural changes in the retina, which could be done by electron microscopy analysis in the future. Third, the abundant expression of Cdr1as in retina seems to be mismatch with the mild retinal phenotype in the Cdr1as KO mice.

As mentioned above, Cdr1as is the most abundant in the mammalian brain, followed by the retina and is expressed at least level or absent in other tissues. It has been reported that loss-of-function Cdr1as in vivo causes miRNA deregulation, mild physiological consequences, and impaired sensorimotor gating (Piwecka et al., 2017). In contrast, the investigation of Cdr1as in zebrafish experiments showed the potential regulatory role of Cdr1as (Memczak et al., 2013); these differences might be due to the different experimental approaches (Oltra et al., 2019). Additionally, the slight effect of Cdr1as on the brain and retina further confirmed that as epigenetic factors, the regulation of some circRNAs or miRNAs is typically relatively minor, even though their high expression in eukaryotes (Jin et al., 2009; Piwecka et al., 2017; Kleaveland et al., 2018).

Taken together, we demonstrated for the first time that Cdr1as is abundantly expressed in the inner retina and that its KO altered retinal miRNA expression patterns as well as the expression of their target genes but had a slight influence on retinal morphogenesis and function.

Materials and Methods

Animals

The Cdr1as KO mice and C57BL/6J mice were bred and maintained in the animal facility of Wenzhou Medical University with a 12-h light/12-h dark cycle and had free access to food and water. All experiments and procedures about mice were approved by the Institutional Animal Care and Use Committee.

Cdr1as KO Mice

The Cdr1as KO mice were generated and maintained on the C57BL/6J background with CRISPR/Cas9-mediated genome editing technology. The sgRNAs (sgRNA1 and sgRNA4) to mouse Cdr1as and Cas9 mRNA were co-injected into fertilized mouse eggs to generate targeted KO offspring. F0 founders were identified by PCR followed by sequence analysis, which were bred to WT mice to test germline transmission and F1 animal generation.

Fundus Photography and High-Resolution Spectral-Domain Optical Coherence Tomography (SD-OCT)

Cdr1as KO and C57BL/6J mice were anesthetized intraperitoneally with pentobarbital sodium. Before the examination, 2.5% hydroxypropyl methylcellulose was dropped into the eyes to improve the connection with the machine (Micron IV, Phoenix Research Labs), and then fundus photography was performed. For SD-OCT measurements, images crossing through the optic nerve were obtained and collected for each eye. The thicknesses of the different retinal layers were measured using Insight software (Pleasanton, CA, United States).

Electroretinography (ERG)

Electroretinography responses in both eyes of mice were carried out as described in the instrument manual (Phoenix Research Laboratories) and as previously described (Jin et al., 2014; Xiang et al., 2017). In brief, mice were dark-adapted overnight and then anesthetized intraperitoneally with pentobarbital sodium. Pupils were dilated with 0.5% tropicamide. A drop of 1% methylcellulose was applied on the cornea to improve the conjunction with the gold wire loop electrode. Ground electrodes and referential needles were punctured into the tail and scalp, respectively. Scotopic ERG was recorded at −2.2 and 0.3 log cd-s/m² stimulus intensity with a 30-s interstimulus interval. Photopic ERG was measured at 0.65 log cd-s/m² with a 0.4-s interstimulus interval after 10 min of light adaptation with a background illumination of 30 cd/m².

Immunofluorescence Staining

Whole eyeballs of mice were extracted immediately after euthanasia. After removing the cornea and lens, the eyecups were fixed in 4% paraformaldehyde for 2 h. Then, retinas were dehydrated in 30% (wt/vol) sucrose and then embedded in embedding medium (Neg-50, Thermo). Sections with 12-μm-thick cryosection slides were cut and washed with PBS, blocked in blocking buffer [4% bovine serum albumin (BSA), 0.5% Triton X-100 in PBS] for 1 h, treated with primary antibody at 4°C overnight, and then incubated with secondary antibody at room temperature for 1 h. The following primary and secondary antibodies were used: mouse anti-Rhodopsin (1:500, Sigma), rabbit anti-Cone-arrestin (1:50, Millipore), mouse anti-Calbindin (1:200, BD), rabbit anti-Recoverin (1:500, Millipore), mouse anti-α-synuclein (1:200), guinea pig anti-Vglut1 (1:200, Millipore), rabbit anti-Pkcα (1:100, Sigma), mouse anti-Ctbp2 (1:100, BD Biosciences), mouse anti-Ctbp2 (1:200, BD), rabbit anti-Pax6 (1:200, Sigma), mouse anti-Rbpms(1:50, Santa Cruz) and donkey anti-rabbit IgG conjugated to Alexa Fluor 488 (1:200, Life Technologies), donkey anti-rabbit IgG conjugated to Alexa Fluor 594 (1:200, Life Technologies), donkey anti-mouse IgG conjugated to Alexa Fluor 594 (1:200, Li-cor Biosciences), and goat anti-guinea pig IgG conjugated to Alexa Fluor 568 (1:200, Abcam). Nuclei were stained with 4, 6-diamidino-2-phenylindole (DAPI, 1:3000, Invitrogen, Carlsbad, CA, United States), Morphologies of the stained retina were imaged using a Leica SP8 laser scanning confocal microscope (Leica, Wetzlar, Germany).

In situ Hybridization

BaseScope is an efficient method of ISH. ACdr1as probe targeting the junction site was designed and BaseScope^TM Reagent Kit—RED supplied by Advanced Cell Diagnostics (ACD) was used. Slides of mouse retinas were dried completely at RT and treated with hydrogen peroxide for 10 min and then protease plus for another 30 min at room temperature; then, probe hybridization was performed strictly according to the manufacturer’s protocol. Images were acquired using a microscope (Nikon Eclipse).

Visual Evoked Potential (VEP)

Visual evoked potential was recorded using the same equipment as ERG. After 5 min of dark adaptation in cages, the mice were anesthetized as previously described. Pupils were dilated with 0.5% tropicamide. A drop of 1% methylcellulose was applied on the cornea to avoid eye drying. During each VEP session, body temperature was maintained at 37°C using a heating pad. The active electrode was positioned on the back of the head at the location of the visual cortex; the reference electrode was put into the cheek, and the ground electrode was placed into the tail. VEPs were evoked by continuousflash of 1.4 and 2.0 Hz, 5 cds/m² white light. VEP signals were recorded by a commercial system (RETIport, Roland Consult GmbH, Germany).

Optomotor Response (OMR)

Mice were placed on a raised platform surrounded by a motorized drum with vertical black and white stripes; the drum could rotate clockwise or anticlockwise. The stripe pattern slowly rotated around the animal at a speed of 12°/s, and triggered the optomotor reflex. Animal behavior was monitored by camera, and the behavior was automatically detected and then analyzed by OptoDrum software (Striatech, Germany). The stimulus pattern was continuously and automatically adjusted during the experiment to find the Cdr1as KO and WT mouse visual thresholds (visual acuity or contrast sensitivity).

RNA Isolation, qRT-PCR, and Small RNA Sequencing

The retinas were isolated from mice and collected in TRIzol (Invitrogen, United States). RNAs were extracted using the RNeasy Kit (Qiagen). For miRNA, qRT-PCR was carried out using the Bulge-Loop miRNA qRT-PCR Starter Kit according to the manufacturer’s protocol. For mRNAs, complementary DNA (cDNA) was synthesized using random primers (Promega) and quantified by FastStart Universal SYBR Green Master Mix (Roche). GAPDH was used as the reference gene. Primers for genes were listed in Supplementary Table 3. For small RNA sequencing, retinas were isolated from 2M WT and Cdr1as KO mice (n = 3); the integrity and quantity of RNA was assessed using Agilent 2200 TapeStation and Qubit2.0, respectively. Small RNA libraries were constructed and sequenced by HiSeq 2500 (Illumina, United States) at Ribobio Co. Ltd. (Ribobio, China). miRDeep2 was used to identify known mature miRNA based on miRBase21 and predict novel miRNAs. The expression levels of miRNAs were normalized by RPM, RPM = (number of reads mapping to miRNA/number of reads in clean data) × 10⁶. Differential expression between WT and Cdr1as KO retina was calculated by edge R algorithm according to the criteria of | Fold Change| ≥ 1.5 and P-value < 0.05.

GO and KEGG Pathway Analysis

GO and KEGG pathway analysis was based on the targets of different expressed 25 miRNAs in Cdr1as KO retina. The target genes of miRNAs were predicted using seven tools, including PITA¹, RNA22², miRNAmap³, microT⁴, miRanda⁵, PicTar⁶, and TargetScan⁷. Target genes predicted by at least three tools and verified by CLIP assay (Zhou et al.; Li et al., 2014) were retained for further GO and KEGG pathway analysis⁸.

Western Blots

Retinas (P300) were isolated, collected, and lyzed in lysis buffer containing 1× PMSF. Protein was then extracted and quantified using a BCA protein assay kit (Invitrogen). Proteins were separated using SDS-PAGE and then analyzed by anti-α-synuclein (1:2000, Abcam), anti-phospho-α-synuclein (1:2000, Wako), and anti-GAPDH (1:1000, KangChen Biotech).

Statistical Analysis

The values shown in the graphs represent averages of several independent experiments and the actual number of samples for each experiment stated in the figure legends. The results are represented as mean ± SEM. The statistical significance was assessed by a two-tailed student’s t-test. ^∗P < 0.05; ^∗∗P < 0.005; ^∗∗∗P < 0.001.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics Statement

The animal study was reviewed and approved by the Institutional Animal Care and Use Committee WMU.

Author Contributions

Z-BJ conceived and supervised the whole study. X-JC, M-LL, Y-HW, HM, ZW, X-YW, and C-JZ performed the experiments. XX participated in the data interpretation. SB, Z-HX, and HZ performed bioinformatics statistical analysis. X-JC wrote the manuscript. Z-BJ revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Key R&D Program of China (2017YFA0105300) and Natural Science Foundation of China (81870690 and 81970838).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2020.565543/full#supplementary-material

Footnotes

1. ^ http://genie.weizmann.ac.il/pubs/mir07/mir07data.html
2. ^ http://cbcsrv.watson.ibm.com/rna22.html
3. ^ http://mirnamap.mbc.nctu.edu.tw/
4. ^ http://www.microrna.gr/microT
5. ^ http://www.microrna.org/microrna/home.do
6. ^ http://www.pictar.org/
7. ^ http://www.targetscan.org/
8. ^ https://david.ncifcrf.gov/summary.jsp

References

Ashwal-Fluss, R., Meyer, M., Pamudurti, N. R., Ivanov, A., Bartok, O., Hanan, M., et al. (2014). circRNA biogenesis competes with pre-mRNA splicing. Mol. Cell. 56, 55–66. doi: 10.1016/j.molcel.2014.08.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Caputto, B. L., and Guido, M. E. (2000). Immediate early gene expression within the visual system: light and circadian regulation in the retina and the suprachiasmatic nucleus. Neurochem. Res. 25, 153–162.

Google Scholar

Chandra, S., Gallardo, G., Fernandez-Chacon, R., Schluter, O. M., and Thomas, C. S. (2005). Sudhof, Alpha-synuclein cooperates with CSPalpha in preventing neurodegeneration. Cell 123, 383–396. doi: 10.1016/j.cell.2005.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Che, Q., Wang, W., Duan, P., Fang, F., Liu, C., Zhou, T., et al. (2019). Downregulation of miR-322 promotes apoptosis of GC-2 cell by targeting Ddx3x. Reprod. Biol. Endocrinol. 17:63.

Google Scholar

Chen, H., Wang, S., Chen, L., Chen, Y., Wu, M., Zhang, Y., et al. (2014). MicroRNA-344 inhibits 3T3-L1 cell differentiation via targeting GSK3beta of Wnt/beta-catenin signaling pathway. FEBS Lett. 588, 429–435. doi: 10.1016/j.febslet.2013.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L. L. (2016). The biogenesis and emerging roles of circular RNAs. Nat. Rev. Mol. Cell. Biol. 17, 205–211. doi: 10.1038/nrm.2015.32

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X. J., Zhang, Z. C., Wang, X. Y., Zhao, H. Q., Li, M. L., and Ma, Y. (2020). The circular RNome of developmental retina in mice. Mol. Ther. Nucleic Acids 19, 339–349. doi: 10.1016/j.omtn.2019.11.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, P., and Frame, S. (2001). The renaissance of GSK3. Nat. Rev. Mol. Cell. Biol. 2, 769–776. doi: 10.1038/35096075

PubMed Abstract | CrossRef Full Text | Google Scholar

Ebbesen, K. K., Hansen, T. B., and Kjems, J. (2017). Insights into circular RNA biology. RNA Biol. 14, 1035–1045. doi: 10.1080/15476286.2016.1271524

PubMed Abstract | CrossRef Full Text | Google Scholar

Genini, S., Guziewicz, K. E., Beltran, W. A., and Aguirre, G. D. (2014). Altered miRNA expression in canine retinas during normal development and in models of retinal degeneration. BMC Genomics 15:172. doi: 10.1186/1471-2164-15-172

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, J. U., Agarwal, V., Guo, H., and Bartel, D. P. (2014). Expanded identification and characterization of mammalian circular RNAs. Genome Biol. 15:409.

Google Scholar

Hansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., et al. (2013). Natural RNA circles function as efficient microRNA sponges. Nature 495, 384–388. doi: 10.1038/nature11993

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeck, W. R., Sorrentino, J. A., Wang, K., Slevin, M. K., Burd, C. E., and Liu, J. (2013). Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 19, 141–157. doi: 10.1261/rna.035667.112

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Z. B., Hirokawa, G., Gui, L., Takahashi, R., Osakada, F., Hiura, Y., et al. (2009). Targeted deletion of miR-182, an abundant retinal microRNA. Mol. Vis. 15, 523–533.

Google Scholar

Jin, Z. B., Huang, X. F., Lv, J. N., Xiang, L., Li, D. Q., Chen, J., et al. (2014). SLC7A14 linked to autosomal recessive retinitis pigmentosa. Nat. Commun. 5:3517.

Google Scholar

Junn, E., Lee, K. W., Jeong, B. S., Chan, T. W., Im, J. Y., and Mouradian, M. M. (2009). Repression of alpha-synuclein expression and toxicity by microRNA-7. Proc. Natl. Acad. Sci. U.S.A. 106, 13052–13057. doi: 10.1073/pnas.0906277106

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D. Y., and Jung, C. S. (2012). Gap junction contributions to the goldfish lectroretinogram at the photopic illumination level. Korean J. Physiol. Pharmacol. 16, 219–224. doi: 10.4196/kjpp.2012.16.3.219

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleaveland, B., Shi, C. Y., Stefano, J., and Bartel, D. P. (2018). A network of noncoding regulatory RNAs Acts in the mammalian brain. Cell 174, 350.e17–362.e17.

Google Scholar

Kok, F. O., Shin, M., Ni, C. W., Gupta, A., Grosse, A. S., and Kirchmaier, B. C. (2015). Reverse genetic screening reveals poor correlation between morpholino-induced and mutant phenotypes in zebrafish. Dev. Cell. 32, 97–108. doi: 10.1016/j.devcel.2014.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Krol, J., Krol, I., Alvarez, C. P., Fiscella, M., Hierlemann, A., Roska, B., et al. (2015). A network comprising short and long noncoding RNAs and RNA helicase controls mouse retina architecture. Nat. Commun. 6:7305.

Google Scholar

Lasda, E., and Parker, R. (2014). Circular RNAs: diversity of form and function. RNA 20, 1829–1842. doi: 10.1261/rna.047126.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Lavoie, J., Illiano, P., Sotnikova, T. D., Gainetdinov, R. R., Beaulieu, J. M., and Hebert, M. (2014). The electroretinogram as a biomarker of central dopamine and serotonin: potential relevance to psychiatric disorders. Biol. Psychiatry 75, 479–486. doi: 10.1016/j.biopsych.2012.11.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 42, 92–97.

Google Scholar

Li, X., Liu, C. X., Xue, W., Zhang, Y., Jiang, S., Yin, Q. F., et al. (2017). Coordinated circRNA Biogenesis and Function with NF90/NF110 in Viral Infection. Mol. Cell. 67, 214–227. doi: 10.1016/j.molcel.2017.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Yang, L., and Chen, L. L. (2018). The Biogenesis, Functions, and Challenges of Circular RNAs. Mol. Cell. 71, 428–442. doi: 10.1016/j.molcel.2018.06.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Huang, C., Bao, C., Chen, L., Lin, M., Wang, X., et al. (2015). Exon-intron circular RNAs regulate transcription in the nucleus. Nat. Struct. Mol. Biol. 22, 256–264. doi: 10.1038/nsmb.2959

PubMed Abstract | CrossRef Full Text | Google Scholar

Lisi, V., Singh, B., Giroux, M., Guzman, E., Painter, M. W., and Cheng, Y. C. (2017). Enhanced neuronal regeneration in the CAST/Ei mouse strain is linked to expression of differentiation markers after injury. Cell. Rep. 20, 1136–1147. doi: 10.1016/j.celrep.2017.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., et al. (2013). Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 495, 333–338. doi: 10.1038/nature11928

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, D. C., Bryan, B. A., Liu, J. P., Liu, W. B., Zhang, L., Qu, J., et al. (2007). Developmental expression of three small GTPases in the mouse eye. Mol. Vis. 13, 1144–1153.

Google Scholar

Oltra, M., Vidal-Gil, L., Maisto, R., Oltra, S. S., Romero, F. J., Sancho-Pelluz, J., et al. (2019). miR302a and 122 are deregulated in small extracellular vesicles from ARPE-19 cells cultured with H2O2. Sci. Rep. 9:17954.

Google Scholar

Piwecka, M., Glazar, P., Memczak, S., Wolf, S. A., Filipchyk, A., and Klironomos, F. (2017). Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function. Science 357:eaam8526. doi: 10.1126/science.aam8526

PubMed Abstract | CrossRef Full Text | Google Scholar

Rybak-Wolf, A., Stottmeister, C., Glazar, P., Jens, M., Pino, N., Giusti, S., et al. (2015). Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol. Cell. 58, 870–885. doi: 10.1016/j.molcel.2015.03.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Salzman, J., Chen, R. E., Olsen, M. N., Wang, P. L., and Brown, P. O. (2013). Cell-type specific features of circular RNA expression. PLoS Genet. 9:e1003777. doi: 10.1371/journal.pgen.1003777

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, T., and Zrenner, E. (1991). Effects of D-1 and D-2 dopamine antagonists on ERG and optic nerve response of the cat. Exp. Eye Res. 52, 425–430. doi: 10.1016/0014-4835(91)90038-g

CrossRef Full Text | Google Scholar

Surguchov, A., McMahan, B., Masliah, E., and Surgucheva, I. (2001). Synucleins in ocular tissues. J. Neurosci. Res. 65, 68–77. doi: 10.1002/jnr.1129

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, N., Xu, H. P., and Wang, P. (2015). Dopamine D2 receptors preferentially regulate the development of light responses of the inner retina. Eur. J. Neurosci. 41, 17–30. doi: 10.1111/ejn.12783

PubMed Abstract | CrossRef Full Text | Google Scholar

Veno, M. T., Hansen, T. B., Veno, S. T., Clausen, B. H., Grebing, M., Finsen, B., et al. (2015). Spatio-temporal regulation of circular RNA expression during porcine embryonic brain development. Genome Biol. 16:245.

Google Scholar

Wang, Y., and Wang, Z. (2015). Efficient backsplicing produces translatable circular mRNAs. RNA 21, 172–179. doi: 10.1261/rna.048272.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Witkovsky, P. (2004). Dopamine and retinal function. Doc. Ophthalmol. 108, 17–40. doi: 10.1023/b:doop.0000019487.88486.0a

CrossRef Full Text | Google Scholar

Xiang, L., Chen, X. J., Wu, K. C., Zhang, C. J., Zhou, G. H., Lv, J. N., et al. (2017). miR-183/96 plays a pivotal regulatory role in mouse photoreceptor maturation and maintenance. Proc. Natl. Acad. Sci. U.S.A. 114, 6376–6381. doi: 10.1073/pnas.1618757114

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, K., Qiu, S., Tian, L., Snider, W. D., Flannery, J. G., Schaffer, D. V., et al. (2016). Wnt regulates proliferation and neurogenic potential of muller glial cells via a Lin28/let-7 miRNA-dependent pathway in adult mammalian retinas. Cell. Rep. 17, 165–178. doi: 10.1016/j.celrep.2016.08.078

PubMed Abstract | CrossRef Full Text | Google Scholar

You, X., Vlatkovic, I., Babic, A., Will, T., Epstein, I., Tushev, G., et al. (2015). Neural circular RNAs are derived from synaptic genes and regulated by development and plasticity. Nat. Neurosci. 18, 603–610. doi: 10.1038/nn.3975

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, K. R., Liu, S., Cai, L., Bin, L., et al. ENCORI: The Encyclopedia of RNA Interactomes.

Google Scholar

Zhu, X., Raina, A. K., Boux, H., Simmons, Z. L., Takeda, A., and Smith, M. A. (2000). Activation of oncogenic pathways in degenerating neurons in Alzheimer disease. Int. J. Dev. Neurosci. 18, 433–437. doi: 10.1016/s0736-5748(00)00010-1

CrossRef Full Text | Google Scholar