miRNA Profiling of Developing Rat Retina in the First Three Postnatal Weeks

The morphogenesis of the mammalian retina depends on the precise control of gene expression during development. Small non-coding RNAs, including microRNAs play profound roles in various physiological and pathological processes via gene expression regulation. A systematic analysis of the expression profile of small non-coding RNAs in developing Wistar rat retinas (postnatally day 5 (P5), P7, P10, P15 and P21) was executed using IonTorrent PGM next-generation sequencing technique to reveal the crucial players in the early postnatal retinogenesis. Our analysis reveals extensive regulatory potential of microRNAs during retinal development. We found a group of microRNAs that show constant high abundance (miR-19, miR-101; miR-181, miR-183, miR-124 and let-7) during the development process. Others are present only in the early stages (miR-20a, miR-206, miR-133, miR-466, miR-1247, miR-3582), or at later stages (miR-29, miR-96, miR-125, miR-344 or miR-664). Further miRNAs were detected which are differentially expressed in time. Finally, pathway enrichment analysis has revealed 850 predicted target genes that mainly participate in lipid-, amino acid- and glycan metabolisms in the examined time-period (P5–P21). P5–P7 transition revealed the importance of miRNAs in glutamatergic synapse and gap junction pathways. Significantly downregulated miRNAs rno-miR-30c1 and 2, rno-miR-205 and rno-miR-503 were detected to target Prkx (ENSRNOG00000003696), Adcy6 (ENSRNOG00000011587), Gnai3 (ENSRNOG00000019465) and Gja1 (ENSRNOG00000000805) genes. The dataset described here will be a valuable resource for clarifying new regulatory mechanisms for retinal development and will greatly contribute to our understanding of the divergence and function of microRNAs. Supplementary Information The online version contains supplementary material available at 10.1007/s10571-023-01347-3.


Introduction
Vertebrates retina is a multi-layered tissue that is formed as a consequence of a series of sequential cellular events during development: neural progenitor proliferation, cell fate specification, neuronal migration, neurite outgrowth and pathfinding, then eventually the formation of synaptic connections (Amini et al. 2018). Seven principal cell types are generated from a pool of multipotent retinal progenitor cells (RPCs) in a conserved sequence: retinal ganglion cells (RGCs) (embryonic day (E) 9-postnatal day (P) 2), then horizontal cells (E10-E15) and cones (E10-E20) followed by the majority of the amacrine cells (E10-P1) at embryonic state. The remaining amacrine cells and most of the rod photoreceptors (E16-P12) as well bipolar cells (E20-P12) and Müller glial cells (E18-P12) are generated postnatally (Rapaport et al. 2004;Reese 2011;Bassett and Wallace 2012;Hoon et al. 2014). Based on morphologies, physiological properties, and/or sublaminar positions of all major types of retinal neurons can be divided into two or more subtypes, except for rods. Amacrine cells and RGCs are known as the most diversified cell types (Masland 2004;Amini et al. 2018;Sharma et al. 2014).
In order to reach the cell number characteristic of the mature retina, each neuronal cell types have a well-characterized, unique proliferative peak. While postnatally generated cells have a normal or "bell-shaped" curve kinetics of cell production, the early generated cells can be considered as "plateau" shape (Rapaport et al. 2004). During postnatal development, amacrine cells show proliferation peaks on day P1, rod photoreceptors from birth to day P6, while Müller and bipolar cells show proliferation peaks on day P6 followed by a rapid, steady decrease until the opening of the eye, on days P12-P13 (Bagnoli et al. 2003;Bassett and Wallace 2012;Reese and Colello 1992;Rapaport et al. 2004). The development of the morphological characteristics of each cell type in the rat retina lasts from day P10 until day P21 (Johansson et al. 2000).
The migration and differentiation of retinal neurons into the appropriate layer is key for assembling the laminar architecture. Subsequently synaptic connections and functional neuronal circuits are formed. The cells of the outer neuroblast layer are only arranged in the outer and inner nuclear layer separated by the outer plexiform layer (OPL) around P6-P8 (Bagnoli et al. 2003;Amini et al. 2018). In the outer retina, photoreceptors contact horizontal cells and bipolar cells, while retinal ganglion cells synapse their presynaptic partners, the amacrine and bipolar cells in the inner retina forming the inner plexiform layer (IPL). Although synaptogenesis in the IPL begins on P3, while in the OPL on P5, the characteristics of a mature cell only occur when the surviving cells enter the appropriate layer. After migration from P7 on the first synaptic vesicles also appear. The expression of receptors, transporters and synaptic proteins are crucial in synaptic transmission (Bagnoli et al. 2003).
Recently, non-coding RNAs especially microRNAs are deeply researched as crucial factors that are members of the complex gene regulation network involved in generating retinal cell diversity (Xu et al. 2007;Sundermeier and Palczewski 2012;Xiang 2013). Now it is well accepted that retinal disorders at least partially have been connected to abnormal miRNAs expression such as diabetic retinopathy, retinitis pigmentosa, retinoblastoma and age-related macular degeneration (Maiorano and Hindges 2012;Andreeva and Cooper 2014). It is also important to know that studies on Dicer ablated animals demonstrated selective loss of miRNAs which caused developmental changes. In conditional knockout Dicer retinas the early born cell types such as RPCs were increased and prolonged, while if late RPC markers including Sox9 and Ascl1 were not expressed, late-born cell types including rods and Müller cells were not generated (Gao 2010;Georgi and Reh 2010;Hackler et al. 2010).
Approximately 1% of the human genome consists of miRNA genes (Andreeva and Cooper 2014). The detection of miRNA molecules nowadays is still very challenging because of their low abundance, small size and the high sequence similarities. However, the improvement in technology over the last decade offers new, highly sensitive methodologies to reveal the existing fine-tuning machinery. Several methods are available to qualify or quantify miRNA expression such as in situ hybridisation, Northern blotting, quantitative reverse transcription polymerase chain reaction (qRT-PCR), microarrays, bead-based flow-cytometry and small-RNA sequencing. While most of them support a miRNA-specific detection, the small RNA-sequencing methodology offers the only way to detect all miRNAs present via a high throughput technics (Pritchard et al. 2012;Kolanowska et al. 2018).
Although the structural development of the vertebrate retina has long been researched, knowledge on miRNAs is still very limited. Here, we describe the small RNA profile of postnatal retina in Wistar WU albino rat by Ion-Torrent PGM sequencing technology. These data provide new insight into the postnatal retinogenesis. Moreover, due to the possible similarity, our results may become a useful resource for developmental biology especially in small RNA studies of other sensory organs.

Animals
Animal handling, housing and experimental procedures were reviewed and approved by the ethics committee of University of Pécs (BAI/35/51-58/2016;PTE/43902/2016). To determine the experimental sample sizes a G-power analysis was conducted (Additional file 1). All efforts were made to minimize pain. Animals were housed as families (female + pups)/ cage distribution, not separated by gender and without the presence of an adult male. We kept max. two adult rats in a box (T3-H type, 378 × 217 × 180 mm polycarbonate box (Acéllabor Kft., Budapest, Hungary). Animals were maintained on a 12:12 h light/dark cycle and were provided with food and water ad libitum. The day of birth is recorded as postnatal day 0 (P0). The animals were randomly allocated to the groups (P5, P7, P10, P15, P21). Wistar WU rats were anesthetized by inhalation using Forane (Abbott Laboratories, Budapest Hungary) prior to sacrifice at the same hour of the day to avoid circadian variation. We had five groups (P5, P7, P10, P15, P21) and carried out two methodologies (sequencing and qPCR, see details later) to examine the retinal development changes of these animals. For sequencing experiment we used the following number of animals: P5 n = 3; P7 n = 6; P10 n = 3; P15 n = 3; P21 = 9. By the end of sequencing, biological sample library originated from RNAs of three individuals. Biological replicates mean at least two different biological samples, all of them originated from RNAs of three individuals. In that case the sequencing outcomes represents finally the miRNA profiles of at least six animals, while technical replicate came from the same barcoded library that was sequenced on another chip. In that case the RNA profile originated from the same three animals: their RNAs were pooled for library preparation and this library was barcoded, then used for template preparation, and finally this template was sequenced on the appropriate chip of the ION Torrent PGM Machine. Similar pooling method was applied by Hackler et al. For qPCR validation the number of the animals was 6, that contained in all cases the samples that were used for sequencing methodology. Each animal was measured in three replicates.
The purpose of this exploratory study was mainly to find over/under expressed small non-coding RNAs in consecutive time points in the developing retinas. As we have relatively few time points (3-5), we handled this simple time course experiment similar to an experiment with several groups (see limma, 9.6.1 guide in R statistic environment).

miRNA Isolation
Total RNA was extracted from the retinal tissues of various developmental stages using NucleoSpin miRNA kit (Macherey-Nagel, Düren, Germany, Cat: 740971.50) following manufacturer's instructions. The miRNA concentration was assessed using Qubit microRNA Assay Kit (ThermoFisher Scientific, MA, USA, Cat: Q32880). The RNA integrity number (RIN), an algorithm for judging the integrity of RNA samples, were evaluated using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California, USA, Cat: G2939BA) following the manufacturing instruction of the RNA 6000 Nano kit (Agilent Technologies, Santa Clara, California, USA, Cat: 5067-1511) and RIN > 7 was considered acceptable. MicroRNAs was also determined using the Agilent Small RNA Kit (Agilent Technologies, Santa Clara, California, USA, Cat: 5067-1548) to have deeper view in the 10-40-nucleotide size range.

Small RNA Sequencing
Small RNA library construction from pooled retinal samples (N = 3, in each age group) were carried out according to the Ion Total RNA-Seq Kit v2 (ThermoFisher Scientific, MA, USA, Cat: 4475936) protocol (Revision E) with slight modification and the experimenters were blinded during the experimental procedures. Enzymatic reaction was performed in half reaction volume, while cleaning procedures were executed with the accurate final volume according to the protocol.
3 μg RNA per three biological replicates in each age group was pooled and enriched for small RNA applying the Magnetic Bead Cleanup Module. Eluted small RNA from the beads was concentrated and used for the hybridisation step at 65 °C for 10 min and 16 °C for 5 min. Ligation reagents (6 μL) were then added and the samples were incubated at 16 °C for 24 h.
After ligation, reverse transcription (RT) was performed on a Veriti™ 96-Well Thermal Cycler (ThermoFisher Scientific, MA, USA, Cat: 4452300): 8 μL RT master mix was added to the ligated RNA samples (10 μL), tubes were incubated at 70 °C for 10 min, and then snap-cooled on ice, the RT enzyme (2 μL of 10× SuperScript® III Enzyme) was added, and reactions (20 μL) were incubated at 42 °C for 30 min. cDNA from the RT reaction was purified and size-selected according to the protocol of Magnetic Bead Cleanup Module and eluted in 12 μL of nuclease-free water.
For amplifying the cDNA, 3 μL of the purified cDNA was combined with 0.5 μL of Ion Xpress™ RNA 3′ Barcode Primer, 0.5 μL of Ion Xpress™ RNA-Seq Barcode BC primer (choose from BC01 to BC16 for different samples, ThermoFisher Scientific, MA, USA, Cat: 4475485), 22.5 μL of Platinum PCR SuperMix High Fidelity, reaction mix (26.5 μL) was then amplified using the following protocol: 94 °C for 2 min; (94 °C for 30 s, 50 °C for 30 s, and 68 °C for 30 s) for 2 cycles; (94 °C for 30 s, 62 °C for 30 s, and 68 °C for 30 s) for 16 cycles; 68 °C for 5 min. The amplified DNA for each sample was completed to 53 μL with nuclease-free water and samples were purified using the Magnetic Bead Cleanup Module and eluted in 15 μL of nuclease-free water.
The quality and quantity of the library was detected following the manufacturer's recommendations of High Sensitivity DNA Chip of Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, California, USA, Cat: 5067-4626). The ratio of miRNA ligation products in total ligation products was calculated using the formula for Barcoded libraries: [Area (94-114 bp)] ÷ [Area (50-300 bp)]. The molar concentration of the barcoded library was determined using size range 50-300 bp, then diluted to a final concentration of ~ 20 pM following the manufacturer's protocol. Equal volumes of three or six diluted barcoded libraries were combined for the next steps depending on the chip type.
For template preparation, the diluted library (15 μL) was used to generate template positive Ion Sphere™ Particles (ISPs) containing clonally amplified DNA. EmulsionPCR (emPCR) was carried out on the Ion OneTouch™ 2 System (ThermoFisher Scientific, MA, USA, Cat: 4474779) with the Ion PGM™ Hi-Q View OT2 200 Kit (ThermoFisher Scientific, MA, USA, Cat: A29900) following the recommended protocol. Template-positive ISPs were enriched with the Ion OneTouch™ ES (ThermoFisher Scientific, MA, USA, Cat: 4469495) following the manufacturer's recommendations.
The Ion 316 or 318 ™ Chip v2 (ThermoFisher Scientific, MA, USA, Cat: 4488149, 4488146) was used for sequencing on the Ion Torrent PGM™ instrument (ThermoFisher Scientific, MA, USA, Cat: 4462921) according to the protocol of Ion PGM™ Hi-Q View Sequencing Kit (ThermoFisher Scientific, MA, USA, Cat: A30044). Sequencing primer and Control Ion Spheres™ Particles of the Ion PGM sequencing kit were added to the enriched, template-positive ISPs. After annealing the sequencing primer, sequencing polymerase was added and a final volume of 30 μL was loaded onto the chips.
To assess reliability of miRNA workflow solution in 316 v2 chip (as Yao et al. 2012), two independent P7 samples were assayed (biological replicates) and P21 sample were sequenced as technical replicates in every sequencing procedure. Some samples were also sequenced in 318 v2 chip (P5, P10, P15 and P21) as biological replicate.

Computational Analysis
Ion Torrent Suite Platform was used to trim the raw sequence data and remove any residual sequencing adapter fragments that remained on the 5′ or 3′ ends. Reads were mapped to the non-coding RNAs from ENSEMBL [Rnor_6.0 (GCA_000001895.4)] using manufacturers' alignment tools for PGM, TMAP algorithm. These aligned BAM (Binary Alignment Map) files were further processed in Galaxy Web-based platform (Afgan et al. 2018) via Cufflinks, Cuffmerge and Cuffdiff (Version 2.2.1.3) application (Trapnell et al. 2010). The expression level of individual transcripts was calculated fragments per kilobase million (FPKM) method in Cufflinks, while geometric library method with 25 of average fragment length was applied in Cuffdiff. Pearson's correlation was used to assess the correlation between read counts per miRNA of biological replicates. Further ascertainment of the most important miRNAs were based on statistical analysis in R Studio Software environment (RStudio Team 2015).
DIANA miRPath 3.0 was used for miRNA target and pathway analysis (Vlachos et al. 2015). Experimentally validated miRNA targets were identified using DIANA-TargetScan and the Genes Union option (an a posteriori method) was used for pathway analysis. The integrated Fisher's Exact test followed by FDR adjustment were used for statistical analysis. Enriched pathways from the KEGG database were exported from the tool along with the corresponding FDRcorrected P values. A significance threshold of FDR = 0.05 was applied to the corrected P values.

Validation of the Sequencing Results
TaqMan® MicroRNA Assays (ThermoFisher Scientific, MA, USA, Cat: 4427975) were applied to validate the miRNA expression profiles. Isolated RNA from six retina in each age group were reverse-transcribed with TaqMan MicroRNA Reverse Transcription Kit (ThermoFisher Scientific, MA, USA, Cat: 4366596) according to the manufacturer's instructions. Briefly, 15 μL of final volume contained 5 μg purified total RNA (10 ng), 10× RT buffer, 100 mM dNTPs, 50 U/μL MultiScribe Reverse Transcriptase (MuLV), 20 U RNase inhibitor and 3 μL 5× stemloop RT-primers. The reactions were incubated in Veriti thermal cycler (ThermoFisher Scientific, MA, USA, Cat: 4452300) at 16 °C for 30 min, 42 °C for 30 min, 85 °C for 5 min and then held in 4 °C. RT-PCR was then preformed using the TaqMan Universal PCR Master Mix and the specific primers from TaqMan MicroRNA Assay specific for each miRNA. U6 small nuclear RNA (U6 snRNA) (Ther-moFisher Scientific, MA, USA, assay ID 001973) was used as an endogenous control. Studied miRNA were miRNA-9 (ThermoFisher Scientific, MA, USA, assay ID 000583) and miRNA-23 (ThermoFisher Scientific, MA, USA, assay ID000399). In a 20-μL PCR reaction, 1.33 µL of cDNA was added to 10 μL of TaqMan Universal PCR master mix II (ThermoFisher Scientific, MA, USA, Cat: 4440043) and 1.0 μL of the 20 × TaqMan miR-specific primers and probe 1 3 mix. The reactions were incubated in a 96-well optical plate at 95 °C for 10 min followed by 40 cycles of 95 °C for 15 s and 60 °C for 60 s, using a Step One Plus Detection system (ThermoFisher Scientific, MA, USA, Cat: 4376600). Negative controls (without template) were also established and all experiments were run in triplicate. Relative quantification was performed by the 2 −ΔΔCT method (Livak and Schmittgen 2001) and P21 was used as reference.

Availability of Data
The gene expression profiling data discussed in this publication have been deposited in the NCBI Gene Expression Omnibus and are accessible through GEO Series accession number GSE159168.

Results
To reveal the main miRNAs quantitative alterations during the early postnatal retinal development total RNA samples were extracted from Wistar WU rat retinal tissues of five developmental stages (P5, P7, P10, P15, P21). In order to monitor the general quality of extracted RNA samples RIN was calculated by Agilent Bioanalyzer Software. High quality and low degradation samples (RIN ≥ 7) were only processed further. To evaluate the reproducibility of the experimental procedures in 316 v2 chip two independent P7 samples were assayed (biological replicates) and P21 samples were sequenced as technical replicates in every sequencing procedure. The miRNA levels between the two independent P7 tissue samples and every P21 samples were characterized by a high Pearson's correlation coefficient (Pearson r = 0.92 and 1.00, respectively), that means high reproducibility of the sequencing (Fig. 1A). Some samples were also sequenced in 318 v2 chip (P5, P10, P15 and P21; as biological replicate) also with a high Pearson's correlation coefficient (Fig. 1B). Results were averaged to reduce experimental errors. In 316 v2 chip we obtained in total 143,796,965 bases and in average 12,182,661 reads for each sample after removal of low-quality reads and contaminants, with the peak length of each sample at about 24 nt (Additional file 2).
RNA reads were annotated based on their sequences mapped onto Rattus norvegicus reference from ENSEMBL (Rnor_6.0 (GCA_000001895.4)) (in average 420 178 reads/ sample), and their relative abundances were determined by their counts, normalized to the number of RPKM methods. To minimize the false positive signal, only reads that were detected in both 316 and 318 v2 sequencings (biological and technical replicates) were used for further bioinformatics analysis. Principal component analysis (PCA) revealed partial separation of samples based on postnatal day of development as shown in Additional file 3.
Based on the non-coding RNA classification methods of ENSEMBL our data generated by sequencings were classified into five functionally important categories including ribosomal RNAs (rRNAs), as well as small RNAs such as microRNAs, snoRNAs and snRNAs (Table 1).
It is well known the rRNAs playing important roles in the protein synthesis machinery, the highest rRNA number was 97 and the lowest value was 91 at the comparison of P10 and P15. The total expression levels for small RNAs that contribute to the biogenesis of rRNAs or to the protein synthesis, mainly characterized by the same tendency during the development period. Slight decrease of the number of non-coding RNAs at P10 suggests an important role of regulation of protein synthesis for the retinal synaptogenesis.
Focusing on the known miRNAs in non-coding RNAs from ENSEMBL, we identified approximately 272 miRNAs in at least one of the five developmental stages. The expression profile of miRNAs could provide indication of their potential functions during development. We observed that although there was no obvious difference in the total number of miRNAs detected in retina during development process, the expression level of different miRNAs in retina was very dynamic over stages (Fig. 2). The most abundant miRNAs at each developmental stages are collected in Table 2.
The differentially expressed miRNAs were collected on consecutive time points with P < 0.05 (Additional file 4). Interestingly, only downregulated miRNAs reached the cutoff at P5 vs. P7 and P10 vs. P15, while mainly upregulated miRNAs were detected at P7 vs. P10. There are some miR-NAs that show constant high abundance during the development process such as miR-19, miR-101; miR-181, miR-183, miR-124 and let-7. While there are miRNAs that appear as most abundant only in the early stages such as miR-20a, miR-206, miR-133b, miR-466, miR-1247 or miR-3582, others are characteristic with high abundance at later stages or increasing with development for example miR-29b, miR-96, mir-125, miR-344 or miR-664.
To further validate the sequencing results by quantitative polymerase chain reaction (qPCR) we choose miR-9 and miR-23 which based on our previous investigations and literature search seemed to have major importance during postnatal retinal development (Leucht et al. 2008;Arora et al. 2010;Gao 2010;La Torre et al. 2013;Qi 2016;Pöstyéni et al. 2021a, b).The relative expression level of miR-9 exhibited the upregulated pattern while that of mir-23 exhibited the downregulated pattern compared to P21. In addition, the expression patterns of miRNAs revealed by qPCR were mostly consistent with sequencing results (Fig. 3) with some exceptions. For example, miRNA-9 at P5 and P7 or mir-23 at P15, when the sequencing results suggest downregulation compared to P21, while qPCR data propose minor 1 3 upregulations. Despite of these minor discrepancies between qPCR and deep-sequencing results, our data are highly correlated and accurate in detection and quantification of the relative expression levels of miRNAs.

Prediction and Functional Analysis of Differentially Expressed miRNA Target Genes
To gain insight into the roles of the differentially expressed miRNAs during postnatal-development of the retina, DIANA web-tool was used to predict potential target genes and to apply pathway enrichment analysis. This enrichment analysis identifies pathways significantly enriched with genes belonging to all genes targeted by at least one selected miRNA. For differentially expressed miRNAs, 850 predicted target genes were annotated with function. It is important to know that each miRNA can have different targets, and different miRNAs can have the same target gene. Between enriched pathways targeted by miR-NAs cut off for the value of log2 fold change was more  1 3 or less than ± 1.5 (P < 0.05) there were as follows: lipid metabolisms (such as unsaturated fatty acid, arachidonic acid or glycerophospholipid), amino acid metabolisms (such as valine, leucine, isoleucine or lysin degradation, tyrosine metabolism) and glycan metabolisms in large numbers. The dataset of P5-P7 transition has shown the crucial role in glutamatergic synapse formation (Additional file 5). Among significantly downregulated miR-NAs rno-miR-30c1 and 2, rno-miR-205 and rno-miR-503 were detected to target Prkx (ENSRNOG00000003696), A d c y 6 ( E N S R N O G 0 0 0 0 0 0 1 1 5 8 7 ) , G n a i 3 ( E N S R N O G 0 0 0 0 0 0 1 9 4 6 5 ) a n d H o m e r 2 (ENSRNOG00000019297) genes as shown in Fig. 4 highlighted. Furthermore, the data analysis of P5-P7 transition has also revealed the importance in gap junction (rn04540) KEGG pathways (Additional file 5).

Discussion
Nowadays increasing number of data are generated by sequencing, especially to understand gene expression regulation at miRNA level in retina. However, there is a limited number of studies that focuses on normal developmental events only with the aim to reveal the fine-tuning roles of miRNAs in the gene regulatory network (Xiang 2013;Xu et al. 2007;Georgi and Reh 2010;Arora et al. 2010;Torre et al. 2013;Coolen et al. 2013;Kapsimali et al. 2007). It is already well-known that miRNAs are crucial for physiological and developmental processes in the central nervous system including in the retina (Sundermeier and Palczewski 2012;Maiorano and Hindges 2012;Gao 2010;Hackler et al. 2010;Linsen et al. 2010). Here, we analyzed expression of miRNAs in Wistar albino rat retina tissue during postnatal development and we found that miRNAs showed a wide diversity of expression pattern. Some of them seemed to be  -mir-20a, rno-mir-133, rno-mir-206, rno-mir-466c, while others exhibited a higher abundance at later days, indicating distinct roles played in controlling the development of late-born retinal cells (rod photoreceptors, bipolar cells and Muller glia) and in maintaining the retinal neurons and/or Muller glia (for example rno-mir-29, rno-mir-96, rno-mir-125, rno-mir-344 (Xu et al. 2007)). Again other miRNAs have shown constant high abundance such as let-7 as the most highly expressed rno-miRNA related to retinal-targets (CLASP2, DMD, DUSP1, ELOVL4, NLK, FZD4, RB1, RDH10, RGS16, SLC25A18; Arora et al. 2007).
Furthermore, it is also well-established that let-7, rnomir-125 and rno-mir-9 play an essential role in the developing mouse retina to regulate the transition from early to late stage retinal progenitors (Andreeva and Cooper 2014).
Other studies emphasize the importance of the rnomir-183/96/182 cluster as a sensory organ specific rno-miRNA. All three of these miRNAs have a similar developmental expression pattern and transcribed as a single polycistronic transcript. The sensory organ specificity of these miRNAs is also supported by sequence analysis of the upstream CpG island containing multiple predicted transcription factor binding sites characteristic of genes expressed in neurosensory cells, such as OTX1, Pou3F2 or Pax5 (Xu et al. 2007;Coolen et al. 2013). In a good agreement with these studies we have found a dramatic increase of rno-mir-183/96/182 by the 10th postnatal day. We have already mentioned that we were able to detect rno-mir-183 and rno-mir-96, but interestingly not rno-mir-182.
Our findings are also broadly in line with the microarray based study of microRNA expression in C57BL/6 mouse eye in the case of rno-miRNAs with increasing or decreasing expression tendency during development in wild-type animals [e.g. declining expression of rnomir-20; increased of rno-mir-29 or rno-mir-125 (Hackler et al. 2010)]. However, there are some differences comparing the results to the study of Hackler et al. (2010). While in our study rno-mir-16 showed constant high expression with gentle alterations (namely there is a slight decrease from P5 to P7, then a peak at P10 after that again a slight decrease to P21), in the developing C57BL/6 mouse retina this rno-miRNA was revealed with a decreasing tendency. At the same time, in the case of rno-mir-30 an increasing expression pattern was detected in wild-type mouse, our detailed time-scaling of postnatal days was able to detect a fluctuating profile: decreasing tendency to a minimum at P7, after that a sharp increment to P15 and finally a decline at P21 again. Finally, rno-mir-103 was assessed in our study to have a gradually increasing profile, while it was not even mentioned in the developing C57BL/6 mouse retina.
In order to better understand the role of miRNAs miRNAtarget interactions have been revealed. DIANA web-tool found their roles in various biological processes. All of our investigations the dataset of P5-P7 transition seems the most exciting from the aspect of retinogenesis. Although the development of the neural retina and visual system, especially synaptogenesis, is complex, studies have shown that the major events correspond to occur at a special period (Rapaport et al. 2004;Reese 2011;Bassett and Wallace 2012;Hoon et al. 2014;Kapsimali et al. 2007).
Our P5-P7 transition dataset seems to underline the possible roles of miRNAs in glutamatergic synapse and gap junction formation. There are remarkable changes at this time in the retina as we mentioned earlier. Formation of the OPL through the rearrangement of the neuroblast layers cells (Bagnoli et al. 2003;Amini et al. 2018), and the proliferation peaks of Müller and bipolar cells have also been observed at that time (Bagnoli et al. 2003;Bassett and Wallace 2012;Reese and Colello 1992;Rapaport et al. 2004) Furthermore the synaptogenesis mainly occurs in the OPL at this time and the first synaptic vesicles appear (Bagnoli et al. 2003). In the neuronal circuitry of the retina chemical transmission is dominated by glutamate and γ-aminobutyric acid (GABA), being the major excitatory and inhibitory neurotransmitters, respectively (Yang 2004). Glutamate is the neurotransmitter of the neurons (photoreceptors, bipolar cells, ganglion cells) of the vertical pathways through the entire retina (Massey 1990;Crooks and Kolb 1992;Thoreson and Witkovsky 1999;Haverkamp and Wässle 2004). Specific vesicular glutamate transporters (vGLUTs) load glutamate into synaptic vesicles. Studies from rodent retina have also revealed that remarkable maturational processes of bipolar cell happen at P6-P8: namely OFF and ON bipolar terminals first invade IPL; OFF-bipolar cell terminals express vGluT1 (P6), before ON-bipolar cell terminals (around P10); and rod bipolar cells also appear early (P6) (Sherry et al. 2003). These events are important because the inception of the ON and OFF pathways segregate at bipolar cell level (Yang 2004).
Although there are no data in the literature about the most of the potential target genes (Prkx, Adcy6, Gnai3, Homer2) predicted by our studies, we assume their essential role in retinogenesis based on the limited number of investigation demonstrate their importance in developmental processes in various tissues especially in the olfactory system (Li et al. 2005;Azaiez et al. 2015;Huang et al. 2016;Guo et al. 2017;Beer-Hammer et al. 2018). Finally, Gja1 [also known as con-nexin43 (Cx43)] is the most abundant gap junction protein in the central nervous system and is expressed primarily on astrocytes and Müller cells, thus responsible for maintaining the local homeostasis (Kerr et al. 2010;Völgyi et al. 2013).

Conclusions
The dataset described here will be a valuable resource for getting deeper insight into the gene regulatory network during retina development. Numerous genes are involved in the formation and differentiation of the retina are targeted by various miRNAs thus creating a complex regulatory network to usher the development of the visual system in stage specific manner. So, these data provide useful framework for studying the expression across development or the impact of pharmacological manipulations at a particular developmental stage, especially in P5-P7 transition.