Whole transcriptome sequencing identifies key circRNAs, lncRNAs, and miRNAs regulating neurogenesis in developing mouse retina

Circular RNA expression profile and functional analysis of circUvrag in light-induced photoreceptor degeneration

BACKGROUND

Circular RNAs (circRNAs) are implicated in retinal pathophysiology; however, their expression profiles and functions in photoreceptor apoptosis are largely unknown. We explored circRNA-expression profiles and circUvrag (host gene: Uvrag, ultraviolet radiation resistance associated gene) function in light-induced photoreceptor apoptosis.

METHODS

Sprague-Dawley rats and 661 W photoreceptor cells were exposed to blue light to establish light-induced photoreceptor degeneration. Differentially expressed circRNAs were identified using microarrays. Potential functions of dysregulated circRNAs were analysed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. CircUvrag expression and localization were evaluated using quantitative RT-PCR and fluorescence in situ hybridization, respectively. CircUvrag overexpression and knockdown were induced using a plasmid and a small interfering RNA, respectively, and retinal function and structure were assessed using scotopic electroretinography, haematoxylin-eosin staining, and TUNEL staining. Microglial migration was assessed using IBA1 immunostaining. The apoptosis ratio of photoreceptor cells in vitro was detected using flow cytometry.

RESULTS

We identified 764 differentially expressed circRNAs, which were potentially related with the development of retinal structures, including neurons, dendrites, and synapses, and might participate in nervous-system pathophysiology. Light exposure enriched circUvrag in the cytoplasm of photoreceptors in the outer nuclear layer (ONL). CircUvrag knockdown decreased photoreceptor apoptosis and microglial migration to the ONL after light exposure, preserving ONL thickness and a-wave amplitude. In vitro, circUvrag knockdown inhibited photoreceptor apoptosis, although circUvrag overexpression slightly promoted photoreceptor apoptosis.

CONCLUSIONS

CircUvrag knockdown attenuated light-induced photoreceptor apoptosis, and might be a potential target in retinal degeneration.

Systemic whole transcriptome analysis identified underlying molecular characteristics and regulatory networks implicated in the retina following optic nerve injury

Optic nerve injuries are severely disrupt the structural and functional integrity of the retina, often leading to visual impairment or blindness. Despite the profound impact of these injuries, the molecular mechanisms involved remain poorly understood. In this study, we performed a comprehensive whole-transcriptome analysis of mouse retina samples after optic nerve crush (ONC) to elucidate changes in gene expression and regulatory networks. Transcriptome analysis revealed a variety of molecular alterations, including 256 mRNAs, 530 lncRNAs, and 37 miRNAs, associated with metabolic, inflammatory, signaling, and biosynthetic pathways in the injured retina. The integrated analysis of co-expression and protein-protein interactions identified an active interconnected module comprising 5 co-expressed proteins (Fga, Serpina1a, Hpd, Slc38a4, and Ahsg) associated with the complement and coagulation cascades. Finally, 5 mRNAs (Fga, Serpinala, Hpd, Slc38a4, and Ahsg), 2 miRNAs (miR-671-5p and miR-3057-5p), and 6 lncRNAs (MSTRG. 1830.1, Gm10814, A530013C23Rik, Gm40634, MSTRG.9514.1, A330023F24Rik) were identified by qPCR in the injured retina, and some of them were validated as critical components of a ceRNA network active in 661W and HEK293T cells through dual-luciferase reporter assays. In conclusion, our study provides comprehensive insight into the complex and dynamic biological mechanisms involved in retinal injury responses and highlights promising potential targets to enhance neuroprotection and restore vision.

Circulating Small Extracellular Vesicles Involved in Systemic Regulation Respond to RGC Degeneration in Glaucoma

Glaucoma is a leading cause of irreversible blindness worldwide and is characterized by progressive retinal ganglion cell (RGC) degeneration and vision loss. Since irreversible neurodegeneration occurs before diagnosable, early diagnosis and effective neuroprotection are critical for glaucoma management. Small extracellular vesicles (sEVs) are demonstrated to be potential novel biomarkers and therapeutics for a variety of diseases. In this study, it is found that intravitreal injection of circulating plasma-derived sEVs (PDEV) from glaucoma patients ameliorated retinal degeneration in chronic ocular hypertension (COH) mice. Moreover, it is found that PDEV-miR-29s are significantly upregulated in glaucoma patients and are associated with visual field defects in progressed glaucoma. Subsequently, in vivo and in vitro experiments are conducted to investigate the possible function of miR-29s in RGC pathophysiology. It is showed that the overexpression of miR-29b-3p effectively prevents RGC degeneration in COH mice and promotes the neuronal differentiation of human induced pluripotent stem cells (hiPSCs). Interestingly, engineered sEVs with sufficient miR-29b-3p delivery exhibit more effective RGC protection and neuronal differentiation efficiency. Thus, elevated PDEV-miR-29s may imply systemic regulation to prevent RGC degeneration in glaucoma patients. This study provides new insights into PDEV-based glaucoma diagnosis and therapeutic strategies for neurodegenerative diseases.

Integrative transcriptomic profiling of ncRNAs and mRNAs in developing mouse lens

In recent years, burgeoning research has underscored the pivotal role of non-coding RNA in orchestrating the growth, development, and pathogenesis of various diseases across organisms. However, despite these advances, our understanding of the specific contributions of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) to lens development remains notably limited. Clarifying the intricate gene regulatory networks is imperative for unraveling the molecular underpinnings of lens-related disorders. In this study, we aimed to address this gap by conducting a comprehensive analysis of the expression profiles of messenger RNAs (mRNAs), lncRNAs, and circRNAs at critical developmental time points of the mouse lens, encompassing both embryonic (E10.5, E12.5, and E16.5) and postnatal stages (P0.5, P10.5, and P60). Leveraging RNA-sequencing technology, we identified key transcripts pivotal to lens development. Our analysis revealed differentially expressed (DE) mRNAs, lncRNAs, and circRNAs across various developmental stages. Particularly noteworthy, there were 1831 co-differentially expressed (CO-DE) mRNAs, 150 CO-DE lncRNAs, and 13 CO-DE circRNAs identified during embryonic stages. Gene Ontology (GO) enrichment analysis unveiled associations primarily related to lens development, DNA conformational changes, and angiogenesis among DE mRNAs and lncRNAs. Furthermore, employing protein-protein interaction networks, mRNA-lncRNA co-expression networks, and circRNA-microRNA-mRNA networks, we predicted candidate key molecules implicated in lens development. Our findings underscore the pivotal roles of lncRNAs and circRNAs in this process, offering fresh insights into the pathogenesis of lens-related disorders and paving the way for future exploration in this field.

Multiple transcriptome analyses reveal mouse testis developmental dynamics

The testes are the organs of gamete production and testosterone synthesis. Up to date, no model system is available for mammalian testicular development, and only few studies have characterized the mouse testis transcriptome from no more than three postnatal ages. To describe the transcriptome landscape of the developing mouse testis and identify the potential molecular mechanisms underlying testis maturation, we examined multiple RNA-seq data of mouse testes from 3-week-old (puberty) to 11-week-old (adult). Sperm cells appeared as expected in 5-week-old mouse testis, suggesting the proper sample collection. The principal components analysis revealed the genes from 3w to 4w clustered away from other timepoints, indicating they may be the important nodes for testicular development. The pairwise comparisons at two adjacent timepoints identified 7,612 differentially expressed genes (DEGs), resulting in 58 unique mRNA expression patterns. Enrichment analysis identified functions in tissue morphogenesis (3-4w), regulation of peptidase activity (4-5w), spermatogenesis (7-8w), and antigen processing (10-11w), suggesting distinct functions in different developmental periods. 50 hub genes and 10 gene cluster modules were identified in the testis maturation process by protein-protein interaction (PPI) network analysis, and the miRNA-lncRNA-mRNA, miRNA-circRNA-mRNA and miRNA-circRNA-lncRNA-mRNA competing endogenous RNA (ceRNA) networks were constructed. The results suggest that testis maturation is a complex developmental process modulated by various molecules, and that some potential RNA-RNA interactions may be involved in specific developmental stages. In summary, this study provides an update on the molecular basis of testis development, which may help to understand the molecular mechanisms of mouse testis development and provide guidance for mouse reproduction.

Comprehensive transcriptional analysis of early dorsal skin development in pigs

Porcine skin is similar to human skin in physiology, anatomy and histology and is often used as a model animal for human skin research. There are few studies on the transcriptome aspects of pig skin during the embryonic period. In this study, RNA sequencing was performed on the dorsal skin of Chenghua sows at embryonic day 56 (E56), embryonic day 76 (E76), embryonic day 105 (E105), and 3 days after birth (D3) to explore RNA changes in pig dorsal skin at four ages. A number of skin-related differential genes were identified by intercomparison between RNAs at four time points, and KEGG functional analysis showed that these differential genes were mainly enriched in metabolic and developmental, immune, and disease pathways, and the pathways enriched in GO analysis were highly overlapping. Collagen is an important part of the skin, with type I collagen making up the largest portion. In this study, collagen type I alpha 1 (COL1A1) and collagen type I alpha 2 (COL1A2) were significantly upregulated at four time points. In addition, lncRNA-miRNA-mRNA and miRNA-circRNA coexpression networks were constructed. The data obtained may help to explain age-related changes in transcriptional patterns during skin development and provide further references for understanding human skin development at the molecular level.

Sstr2 Defines the Cone Differentiation-Competent Late-Stage Retinal Progenitor Cells in the Developing Mouse Retina

Cone cell death is a characteristic shared by various retinal degenerative disorders, such as cone-rod dystrophy, Stargardt disease, achromatopsia, and retinitis pigmentosa. This leads to conditions like color blindness and permanently impaired visual acuity. Stem cell therapy focused on photoreceptor replacement holds promise for addressing these conditions. However, identifying surface markers that aid in enriching retinal progenitor cells (RPCs) capable of differentiating into cones remains a complex task. In this study, we employed single-cell RNA sequencing to scrutinize the transcriptome of developing retinas in C57BL/6J mice. This revealed the distinctive expression of somatostatin receptor 2 (Sstr2), a surface protein, in late-stage RPCs exhibiting the potential for photoreceptor differentiation. In vivo lineage tracing experiments verified that Sstr2+ cells within the late embryonic retina gave rise to cones, amacrine and horizontal cells during the developmental process. Furthermore, Sstr2+ cells that were isolated from the late embryonic mouse retina displayed RPC markers and exhibited the capability to differentiate into cones in vitro. Upon subretinal transplantation into both wild-type and retinal degeneration 10 (rd10) mice, Sstr2+ cells survived and expressed cone-specific markers. This study underscores the ability of Sstr2 to enrich late-stage RPCs primed for cone differentiation to a large extent. It proposes the utility of Sstr2 as a biomarker for RPCs capable of generating cones for transplantation purposes.

Differentially expressed genes in orbital adipose/connective tissue of thyroid-associated orbitopathy

Background

Thyroid-associated orbitopathy (TAO) is a disease associated with autoimmune thyroid disorders and it can lead to proptosis, diplopia, and vision-threatening compressive optic neuropathy. To comprehensively understand the molecular mechanisms underlying orbital adipogenesis in TAO, we characterize the intrinsic molecular properties of orbital adipose/connective tissue from patients with TAO and control individuals.

Methods

RNA sequencing analysis (RNA-seq) was performed to measure the gene expression of orbital adipose/connective tissues of TAO patients. Differentially expressed genes (DEGs) were detected and analyzed through Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and Gene Set Enrichment Analysis (GSEA). The protein-protein interaction (PPI) network was constructed using the STRING database, and hub genes were identified by the Cytoscape plug-in, cytoHubba. We validated several top DEGs through quantitative real-time polymerase chain reaction (qRT-PCR).

Results

We identified 183 DEGs in adipose tissue between TAO patients (n = 3) and control patients (n = 3) through RNA sequencing, including 114 upregulated genes and 69 downregulated genes. The PPI network of these DEGs had 202 nodes and 743 edges. PCR-based validation results of orbital adipose tissue showed multiple top-ranked genes in TAO patients (n = 4) are immune and inflammatory response genes compared with the control individual (n = 4). They include ceruloplasmin isoform x3 (CP), alkaline tissue-nonspecific isozyme isoform x1 (ALPL), and angiotensinogen (AGT), which were overrepresented by 2.27- to 6.40-fold. Meanwhile, protein mab-21-like 1 (MAB21L1), phosphoinositide 3-kinase gamma-subunit (PIK3C2G), and clavesin-2 (CLVS2) decreased by 2.6% to 32.8%. R-spondin 1 (RSPO1), which is related to oogonia differentiation and developmental angiogenesis, was significantly downregulated in the orbital muscle tissues of patients with TAO compared with the control groups (P = 0.024).

Conclusions

Our results suggest that there are genetic differences in orbital adipose-connective tissues derived from TAO patients. The upregulation of the inflammatory response in orbital fat of TAO may be consistent with the clinical phenotype like eyelid edema, exophthalmos, and excess tearing. Downregulation of MAB21L1, PIK3C2G, and CLVS2 in TAO tissue demonstrates dysregulation of differentiation, oxidative stress, and developmental pathways.

Long Noncoding RNA MIAT Modulates Chronic Retinal Ischemia-Reperfusion Injury in Mice via the microRNA-203-3p/SNAI2 Axis

Retinal ischemia-reperfusion injury (RIRI) is a vital pathological process of multiple ocular diseases. This study aimed at investigating the effects of the MIAT/miR-203-3p/SNAI2 axis on RIRI. RIRI was produced by inducing an exceedingly high intraocular pressure (IOP) in mice. Mouse retinal ganglion cells (RGCs) were subjected to oxygen-glucose deprivation/reoxygenation (OGD/R) to mimic in vitro models. Relevant oligonucleotides or plasmids were transfected into OGD/R-induced RGCs in vitro or injected into RIRI mice models in vivo via a vitreous cavity. The findings of our paper indicated that MIAT and SNAI2 were highly expressed and miR-203-3p was lowly expressed in mouse RIRI tissues and OGD/R-induced RGCs. Interfering MIAT promoted the viability of OGD/R-induced RGCs, decreased apoptosis, and reduced oxidative stress in vitro. Silencing MIAT increased retinal neuronal cell numbers and decreased retinal neuronal cell apoptosis in mouse RIRI tissues in vivo. MIAT sponged miR-203-3p, and miR-203-3p targeted and inhibited SNAI2 expression. SNAI2 up-regulation or miR-203-3p down-regulation reversed the protective effects of MIAT down-regulation on RIRI in mice and OGD/R-induced RGCs. MIAT sponges miR-203-3p upregulated the expression of SNAI2, thereby promoting RIRI in mice. In summary, MIAT may be a therapeutic target for the treatment of chronic RIRI.

Integrated Transcriptome Analysis of Long Noncoding RNA and mRNA in Developing and Aging Mouse Retina

Mice have emerged as a widely employed model for investigating various retinal diseases. However, the availability of comprehensive datasets capturing the entire developmental and aging stages of the mouse retina, particularly during the elderly period, encompassing integrated lncRNA and mRNA expression profiles, is limited. In this study, we assembled a total of 18 retina samples from mice across 6 distinct stages of development and aging (5 days, 3 weeks, 6 weeks, 10 weeks, 6 months, and 15 months) to conduct integrated lncRNA and mRNA sequencing analysis. This invaluable dataset offers a comprehensive transcriptomic resource of mRNA and lncRNA expression profiles during the natural progression of retinal development and aging. The discoveries stemming from this investigation will significantly contribute to the elucidation of the underlying molecular mechanisms associated with various retinal diseases, such as congenital retinal dysplasia and retinal degenerative diseases.

Analysis of shared ceRNA networks and related-hub genes in rats with primary and secondary photoreceptor degeneration

Introduction

Photoreceptor degenerative diseases are characterized by the progressive death of photoreceptor cells, resulting in irreversible visual impairment. However, the role of competing endogenous RNA (ceRNA) in photoreceptor degeneration is unclear. We aimed to explore the shared ceRNA regulation network and potential molecular mechanisms between primary and secondary photoreceptor degenerations.

Methods

We established animal models for both types of photoreceptor degenerations and conducted retina RNA sequencing to identify shared differentially expressed long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs). Using ceRNA regulatory principles, we constructed a shared ceRNA network and performed function enrichment and protein-protein interaction (PPI) analyses to identify hub genes and key pathways. Immune cell infiltration and drug-gene interaction analyses were conducted, and hub gene expression was validated by quantitative real-time polymerase chain reaction (qRT-PCR).

Results

We identified 37 shared differentially expressed lncRNAs, 34 miRNAs, and 247 mRNAs and constructed a ceRNA network consisting of 3 lncRNAs, 5 miRNAs, and 109 mRNAs. Furthermore, we examined 109 common differentially expressed genes (DEGs) through functional annotation, PPI analysis, and regulatory network analysis. We discovered that these diseases shared the complement and coagulation cascades pathway. Eight hub genes were identified and enriched in the immune system process. Immune infiltration analysis revealed increased T cells and decreased B cells in both photoreceptor degenerations. The expression of hub genes was closely associated with the quantities of immune cell types. Additionally, we identified 7 immune therapeutical drugs that target the hub genes.

Discussion

Our findings provide new insights and directions for understanding the common mechanisms underlying the development of photoreceptor degeneration. The hub genes and related ceRNA networks we identified may offer new perspectives for elucidating the mechanisms and hold promise for the development of innovative treatment strategies.

LncRNA MALAT 1/miR-625-3p/HIF-1α axis regulates the EMT of hypoxia-induced RPE cells by activating NF-κB/snail signaling

The retina may undergo structural and functional damage as a result of hypoxia, which could lead to permanent blindness. As competing endogenous RNAs (ceRNAs), lncRNAs are essential in eye disorders. The biological function of lncRNA MALAT 1 and its potential mechanisms in hypoxic-ischemic retinal diseases are still unknown. MALAT 1 and miR-625-3p expression alterations in hypoxia-treated RPE cells were examined using qRT-PCR. The target binding relationships between MALAT 1 and miR-625-3p, as well as between miR-625-3p and HIF-1α, were identified utilizing bioinformatics analysis and dual luciferase reporter assay. We observed that si-MALAT 1 and miR-625-3p mimic both reduced apoptosis and epithelial-mesenchymal transition (EMT) in hypoxic RPE cells, whereas si-MALAT 1 was reversed by miR-625-3p inhibitor. Further, we carried out a mechanistic investigation, and rescue assays demonstrated that MALAT 1 sponging miR-625-3p influenced HIF-1α expression and consequently took part in the NF-κB/Snail signaling pathway, which regulated apoptosis and EMT. In conclusion, our research had shown that the MALAT 1/miR-625-3p/HIF-1α axis drove the progression of hypoxic-ischemic retinal disorders and may serve as a promising predictive biomarker for their therapeutic and diagnostic targets.

The landscape of the long non-coding RNAs in developing mouse retinas

BACKGROUND

The long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes. Nevertheless, a global view of its expression and function in the mouse retina, a crucial model for neurogenesis study, still needs to be made available.

RESULTS

Herein, by integrating the established gene models and the result from ab initio prediction using short- and long-read sequencing, we characterized 4,523 lncRNA genes (MRLGs) in developing mouse retinas (from the embryonic day of 12.5 to the neonatal day of P28), which was so far the most comprehensive collection of retinal lncRNAs. Next, derived from transcriptomics analyses of different tissues and developing retinas, we found that the MRLGs were highly spatiotemporal specific in expression and played essential roles in regulating the genesis and function of mouse retinas. In addition, we investigated the expression of MRLGs in some mouse mutants and revealed that 97 intergenic MRLGs might be involved in regulating differentiation and development of retinal neurons through Math5, Isl1, Brn3b, NRL, Onecut1, or Onecut2 mediated pathways.

CONCLUSIONS

In summary, this work significantly enhanced our knowledge of lncRNA genes in mouse retina development and provided valuable clues for future exploration of their biological roles.

MicroRNAs in the Mouse Developing Retina

The retina is among the highest organized tissues of the central nervous system. To achieve such organization, a finely tuned regulation of developmental processes is required to form the retinal layers that contain the specialized neurons and supporting glial cells to allow precise phototransduction. MicroRNAs are a class of small RNAs with undoubtful roles in fundamental biological processes, including neurodevelopment of the brain and the retina. This review provides a short overview of the most important findings regarding microRNAs in the regulation of retinal development, from the developmental-dependent rearrangement of the microRNA expression program to the key roles of particular microRNAs in the differentiation and maintenance of retinal cell subtypes.

Involvement of CircRNA Expression Profile in Diabetic Retinopathy and Its Potential Diagnostic Value

Background: Circular RNAs (circRNAs), a class of non-coding and undegradable RNAs, play many pathological functions by acting as miRNA sponges, interacting with RNA-binding proteins, and others. The recent literature indicates that circRNAs possess the advanced superiority for the early screening of diabetic retinopathy (DR).

Methods: CircRNA sources of peripheral blood mononuclear cells (PBMCs) from healthy controls (n = 4), diabetes mellitus patients (DM) (n = 4), and DR patients (n = 4) were extracted for circular RNA microarray analysis. Enriched biological modules and signaling pathways were analyzed by Gene Ontology Enrichment and Kyoto Encyclopedia of Genes and Genomes analysis, respectively. Real-time quantitative reverse transcription PCR (RT-qPCR) was performed to validate differentiated levels of several circRNAs (fold change ≥2, p < .05) in different groups of healthy control subjects (n = 20), DM patients (n = 60), and DR patients (n = 42). Based on our clinical data from DR, the diagnostic performance of candidate circRNAs was measured by operating characteristic curves (ROCs). Subsequently, their circRNA–miRNA networks were constructed by bioinformatics analysis.

Results: Circular RNA microarray analysis was performed, and 2,452 and 289 circRNAs were screened with differential expression in DR patients compared to healthy controls and DM patients, respectively. Enrichment analyses showed that circRNAs in DR patients were enriched in extracellular matrix (ECM)–receptor interaction and focal adhesion pathways. The top 5 differential circRNAs in circRNA microarray analysis were subsequently quantified and verified by RT-qPCR. Consistently, a significant 2.2-fold reduction of hsa_circ_0095008 and 1.7-fold increase in hsa_circ_0001883 were identified in DR patients compared to DM patients. Meanwhile, the area under curves of hsa_circ_0095008 and hsa_circ_0001883 were 0.6710 (95% CI, 0.5646–0.7775) (p = 0.003399) and 0.6071 (95% CI, 0.4953–0.7189) (p = 0.06644), respectively, indicating a good diagnostic value.

Conclusion: Our study provided a new sight for the pathological mechanism of DR and revealed the potential value of hsa_circ_0095008 and hsa_circ_0001883 as diagnostic biomarkers for the early diagnosis of DR patients.