Retinal transcriptome profiling at transcription start sites: a cap analysis of gene expression early after axonal injury

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.

Non-invasive harvesting of conjunctival cells for whole transcriptome sequencing

The purpose of this study was to determine the feasibility of using a non-invasive technique, the EYEPRIM™ conjunctival cell impression device, to harvest sufficient RNA from conjunctival cells for the whole-transcriptome sequencing. Conjunctival cells from 40 participants were collected using an EYEPRIM™ conjunctival cell impression device. RNA was extracted from the samples, followed by library construction and transcriptome sequencing. Quality checks were performed for each technical step of the experiment, and the feasibility of this procedure was examined. RNA of sufficient yield and quality was successfully extracted following additional disruption and homogenization of the conjunctival cells and collection of two impression samples per eye. Successful library preparation and RNA sequencing were performed, with all 40 samples passing the various quality checks used for each step. In conclusion, harvesting cells from the ocular surface using an impression cytology device yields good quality and sufficient mRNA for whole transcriptome sequencing to study diseases of the eye. This technique provides a convenient alternative to using post-mortem tissues or surgical excisions.

Whole Exome Sequencing Reveals Novel Candidate Genes in Familial Forms of Glaucomatous Neurodegeneration

Glaucoma is the largest cause of irreversible blindness with a multifactorial genetic etiology. This study explores novel genes and gene networks in familial forms of primary open angle glaucoma (POAG) and primary angle closure glaucoma (PACG) to identify rare mutations with high penetrance. Thirty-one samples from nine MYOC-negative families (five POAG and four PACG) underwent whole-exome sequencing and analysis. A set of prioritized genes and variations were screened in an independent validation cohort of 1536 samples and the whole-exome data from 20 sporadic patients. The expression profiles of the candidate genes were analyzed in 17 publicly available expression datasets from ocular tissues and single cells. Rare, deleterious SNVs in AQP5, SRFBP1, CDH6 and FOXM1 from POAG families and in ACACB, RGL3 and LAMA2 from PACG families were found exclusively in glaucoma cases. AQP5, SRFBP1 and CDH6 also revealed significant altered expression in glaucoma in expression datasets. Single-cell expression analysis revealed enrichment of identified candidate genes in retinal ganglion cells and corneal epithelial cells in POAG; whereas for PACG families, retinal ganglion cells and Schwalbe's Line showed enriched expression. Through an unbiased exome-wide search followed by validation, we identified novel candidate genes for familial cases of POAG and PACG. The SRFBP1 gene found in a POAG family is located within the GLC1M locus on Chr5q. Pathway analysis of candidate genes revealed enrichment of extracellular matrix organization in both POAG and PACG.

Differential Retinal Ganglion Cell Vulnerability, A Critical Clue for the Identification of Neuroprotective Genes in Glaucoma

Retinal ganglion cells (RGCs) are the neurons in the retina which directly project to the brain and transmit visual information along the optic nerve. Glaucoma, one of the leading causes of blindness, is characterized by elevated intraocular pressure (IOP) and degeneration of the optic nerve, which is followed by RGC death. Currently, there are no clinical therapeutic drugs or molecular interventions that prevent RGC death outside of IOP reduction. In order to overcome these major barriers, an increased number of studies have utilized the following combined analytical methods: well-established rodent models of glaucoma including optic nerve injury models and transcriptomic gene expression profiling, resulting in the successful identification of molecules and signaling pathways relevant to RGC protection. In this review, we present a comprehensive overview of pathological features in a variety of animal models of glaucoma and top differentially expressed genes (DEGs) depending on disease progression, RGC subtypes, retinal regions or animal species. By comparing top DEGs among those different transcriptome profiles, we discuss whether commonly listed DEGs could be defined as potential novel therapeutic targets in glaucoma, which will facilitate development of future therapeutic neuroprotective strategies for treatments of human patients in glaucoma.

Transcriptional profiling of the chick retina identifies down-regulation of VIP and UTS2B genes during early lens-induced myopia

Gene expression of the chick retina was examined during the early development of lens-induced myopia (LIM) using whole transcriptome sequencing. Monocular treatment of the right eyes with −10 diopter (D)...

Construction and analysis of mRNA, lncRNA, and transcription factor regulatory networks after retinal ganglion cell injury

Retinal ganglion cell (RGC) injury is a critical pathological feature of several optic neurodegenerative diseases. The regulatory mechanisms underlying RGC injury remain poorly understood. Recent evidence has highlighted the important roles of long noncoding RNAs (lncRNAs) in degenerative neuropathy but few studies have focused on lncRNAs associated with RGC injury. In this study, we analyzed dysregulated lncRNAs associated with RGC injury, their potential regulatory functions, and the molecular mechanisms underlying the regulation of lncRNAs and transcription factors (TFs). We analyzed lncRNA and mRNA profiles in the GSE142881 dataset associated with RGC injury and identified 1049 differentially expressed genes (DEGs), with 18 differentially expressed (DE) TFs among 883 DE mRNAs and 312 DE lncRNAs. The predicted DE lncRNAs and DE mRNAs were used to construct a lncRNA-mRNA co-expression network. Functional enrichment analysis was performed to explore the functions of the lncRNAs and mRNAs. The co-expression network between DE lncRNAs and DE mRNAs was highly enriched in inflammatory and immune-related pathways, indicating that they play role in the process of RGC injury. Among the DE mRNAs, we screened 18 DE TFs, including activating transcription factor 3 (ATF3), associated with RGC injury. Co-expression analysis predicted that 13 lncRNAs were potential binding targets of ATF3. The screening of the potential targets of these 13 lncRNAs showed that they were also significantly enriched in functional pathways associated with inflammation and apoptosis. After analysis, we constructed the mRNA-ATF3-lncRNA regulatory network after RGCs injury. In summary, we identified the gene module associated with immune and inflammatory responses after optic nerve injury and constructed a regulatory network of lncRNA-TF-mRNA. The results indicate that lncRNAs, by binding to TFs, can regulate downstream genes and function during RGC injury. The results provide a foundation for further studies of the mechanism of RGC injury and provide insight into the clinical diagnosis and investigation direction of neurodegenerative diseases such as traumatic optic neuropathy and glaucoma.

Commonalities of optic nerve injury and glaucoma-induced neurodegeneration: Insights from transcriptome-wide studies

Glaucoma is a collection of diseases that lead to an irreversible vision loss due to damage of retinal ganglion cells (RGCs). Although the underlying events leading to RGC death are not fully understood, recent research efforts are beginning to define the genetic changes that play a critical role in the initiation and progression of glaucomatous injury and RGC death. Several genetic and experimental animal models have been developed to mimic glaucomatous neurodegeneration. These models differ in many respects but all result in the loss of RGCs. Assessing transcriptional changes across different models could provide a more complete perspective on the molecular drivers of RGC degeneration. For the past several decades, changes in the retinal transcriptome during neurodegeneration process were defined using microarray methods, RNA sequencing and now single cell RNA sequencing. It is understood that these methods have strengths and weaknesses due to technical differences and variations in the analytical tools used. In this review, we focus on the use of transcriptome-wide expression profiling of the changes occurring as RGCs are lost across different glaucoma models. Commonalities of optic nerve crush and glaucoma-induced neurodegeneration are identified and discussed.

Insights into the regulatory molecules involved in glaucoma pathogenesis

Glaucoma is an important cause of irreversible blindness, characterized by optic nerve anomalies. Increased intraocular pressure (IOP) and aging are major risk factors. Retinal ganglion cells and trabecular meshwork cells are certainly involved in the etiology of glaucoma. Glaucoma is usually a complex disease, and various genes and functions may contribute to its etiology. Among these may be genes that encode regulatory molecules. In this review, regulatory molecules including 18 transcription factors (TFs), 195 microRNAs (miRNAs), 106 long noncoding RNAs (lncRNAs), and two circular RNAs (circRNAs) that are reasonable candidates for having roles in glaucoma pathogenesis are described. The targets of the regulators are reported. Glaucoma-related features including apoptosis, stress responses, immune functions, ECM properties, IOP, and eye development are affected by the targeted genes. The targeted genes that are frequently targeted by multiple regulators most often affect apoptosis and the related features of cell death and cell survival. BCL2, CDKN1A, and TP53 are among the frequent targets of three types of glaucoma-relevant regulators, TFs, miRNAs, and lncRNAs. TP53 was itself identified as a glaucoma-relevant TF. Several of the glaucoma-relevant TFs are themselves among frequent targets of regulatory molecules, which is consistent with existence of a complex network involved in glaucoma pathogenesis.

Comprehensive identification of mRNA isoforms reveals the diversity of neural cell-surface molecules with roles in retinal development and disease

Genes encoding cell-surface proteins control nervous system development and are implicated in neurological disorders. These genes produce alternative mRNA isoforms which remain poorly characterized, impeding understanding of how disease-associated mutations cause pathology. Here we introduce a strategy to define complete portfolios of full-length isoforms encoded by individual genes. Applying this approach to neural cell-surface molecules, we identify thousands of unannotated isoforms expressed in retina and brain. By mass spectrometry we confirm expression of newly-discovered proteins on the cell surface in vivo. Remarkably, we discover that the major isoform of a retinal degeneration gene, CRB1, was previously overlooked. This CRB1 isoform is the only one expressed by photoreceptors, the affected cells in CRB1 disease. Using mouse mutants, we identify a function for this isoform at photoreceptor-glial junctions and demonstrate that loss of this isoform accelerates photoreceptor death. Therefore, our isoform identification strategy enables discovery of new gene functions relevant to disease.

Ingenuity pathway analysis of differentially expressed genes involved in signaling pathways and molecular networks in RhoE gene‑edited cardiomyocytes

RhoE/Rnd3 is an atypical member of the Rho superfamily of proteins, However, the global biological function profile of this protein remains unsolved. In the present study, a RhoE‑knockout H9C2 cardiomyocyte cell line was established using CRISPR/Cas9 technology, following which differentially expressed genes (DEGs) between the knockout and wild‑type cell lines were screened using whole genome expression gene chips. A total of 829 DEGs, including 417 upregulated and 412 downregulated, were identified using the threshold of fold changes ≥1.2 and P<0.05. Using the ingenuity pathways analysis system with a threshold of ‑Log (P‑value)>2, 67 canonical pathways were found to be enriched. Many of the detected signaling pathways, including that of oncostatin M signaling, were found to be associated with the inflammatory response. Subsequent disease and function analysis indicated that apart from cardiovascular disease and development function, RhoE may also be involved in other diseases and function, including organismal survival, cancer, organismal injury and abnormalities, cell‑to‑cell signaling and interaction, and molecular transport. In addition, 885 upstream regulators were enriched, including 59 molecules that were predicated to be strongly activated (Z‑score >2) and 60 molecules that were predicated to be significantly inhibited (Z‑scores <‑2). In particular, 33 regulatory effects and 25 networks were revealed to be associated with the DEGs. Among them, the most significant regulatory effects were 'adhesion of endothelial cells' and 'recruitment of myeloid cells' and the top network was 'neurological disease', 'hereditary disorder, organismal injury and abnormalities'. In conclusion, the present study successfully edited the RhoE gene in H9C2 cells using CRISPR/Cas9 technology and subsequently analyzed the enriched DEGs along with their associated canonical signaling pathways, diseases and functions classification, upstream regulatory molecules, regulatory effects and interaction networks. The results of the present study should facilitate the discovery of the global biological and functional properties of RhoE and provide new insights into role of RhoE in human diseases, especially those in the cardiovascular system.

Bioenergetic shift and actin cytoskeleton remodelling as acute vascular adaptive mechanisms to angiotensin II in murine retina and ophthalmic artery

Ocular vascular dysfunction is a major contributing factor to the pathogenesis of glaucoma. In recent years, there has been a renewed interest in the role of angiotensin II (Ang II) in mediating the disease progression. Despite its (patho)physiological importance, the molecular mechanisms underlying Ang II-mediated oxidative stress remain largely unexplored in the ocular vasculature. Here, we provide the first direct evidence of the alterations of proteome and signalling pathways underlying Ang II-elicited oxidative insult independent of arterial pressure changes in the ophthalmic artery (OA) and retina (R) employing an in vitro experimental model. Both R and OA were isolated from male C57Bl/6J mice (n = 15/group; n = 5/biological replicate) and incubated overnight in medium containing either vehicle or Ang II (0.1 μM) at physiological conditions. Label-free quantitative mass spectrometry (MS)-based proteomics analysis identified a differential expression of 107 and 34 proteins in the R and OA, respectively. Statistical and bioinformatics analyses revealed that protein clusters involved in actin cytoskeleton and integrin-linked kinase signalling were significantly activated in the OA. Conversely, a large majority of differentially expressed retinal proteins were involved in dysregulation of numerous energy-producing and metabolic signalling pathways, hinting to a possible shift in retinal cell bioenergetics. Particularly, Ang II-mediated downregulation of septin-7 (Sept7; p < 0.01) and superoxide dismutase [Cu-Zn] (Sod1; p < 0.05), and upregulation of troponin T, fast skeletal muscle (Tnnt3; p < 0.05) and tropomyosin alpha-3 chain (Tpm3; p < 0.01) in the OA, and significant decreased expressions of two crystallin proteins (Cryab; p < 0.05 and Crybb2; p < 0.0001) in the R were verified at the mRNA level, corroborating our proteomics findings. In summary, these results demonstrated that exogenous application of Ang II over an acute time period caused impairment of retinal bioenergetics and cellular demise, and actin cytoskeleton-mediated vascular remodelling in the OA.

•

Acute Ang II stimulation elicits oxidative stress in ocular vasculature without pressor effect. .

•

Dysregulation of energy-producing and metabolic pathways are implicated in the retina. .

•

Actin cytoskeleton remodelling are vascular adaptation processes in the ophthalmic artery. .

Systems Genetics of Optic Nerve Axon Necrosis During Glaucoma

In this study, we identify genomic regions that modulate the number of necrotic axons in optic nerves of a family of mice, some of which have severe glaucoma, and define a set of high priority positional candidate genes that modulate retinal ganglion cell (RGC) axonal degeneration. A large cohort of the BXD family were aged to greater than 13 months of age. Optic nerves from 74 strains and the DBA/2J (D2) parent were harvested, sectioned, and stained with p-phenylenediamine. Numbers of necrotic axons per optic nerve cross-section were counted from 1 to 10 replicates per genotype. Strain means and standard errors were uploaded into GeneNetwork 2 for mapping and systems genetics analyses (Trait 18614). The number of necrotic axons per nerve ranged from only a few hundred to more than 4,000. Using conventional interval mapping as well as linear mixed model mapping, we identified a single locus on chromosome 12 between 109 and 112.5 Mb with a likelihood ratio statistic (LRS) of ~18.5 (p genome-wide ~0.1). Axon necrosis is not linked to locations of major known glaucoma genes in this family, including Gpnmb, Tyrp1, Cdh11, Pou6f2, and Cacna2d1. This indicates that although these genes contribute to pigmentary dispersion or elevated IOP, none directly modulates axon necrosis. Of 156 positional candidates, eight genes—CDC42 binding protein kinase beta (Cdc42bpb); eukaryotic translation initiation factor 5 (Eif5); BCL2-associated athanogene 5 (Bag5); apoptogenic 1, mitochondrial (Apopt1); kinesin light chain 1 (Klc1); X-ray repair cross complementing 3 (Xrcc3); protein phosphatase 1, regulatory subunit 13B (Ppp1r13b); and transmembrane protein 179 (Tmem179)—passed stringent criteria and are high priority candidates. Several candidates are linked to mitochondria and/or axons, strengthening their plausible role as modulators of ON necrosis. Additional studies are required to validate and/or eliminate plausible candidates. Surprisingly, IOP and ON necrosis are inversely correlated across the BXD family in mice >13 months of age and these two traits share few genes among their top ocular and retinal correlates. These data suggest that the two traits are independently modulated or that a more complex and multifaceted approach is required to reveal their association.

Maser: one-stop platform for NGS big data from analysis to visualization

A major challenge in analyzing the data from high-throughput next-generation sequencing (NGS) is how to handle the huge amounts of data and variety of NGS tools and visualize the resultant outputs. To address these issues, we developed a cloud-based data analysis platform, Maser (Management and Analysis System for Enormous Reads), and an original genome browser, Genome Explorer (GE). Maser enables users to manage up to 2 terabytes of data to conduct analyses with easy graphical user interface operations and offers analysis pipelines in which several individual tools are combined as a single pipeline for very common and standard analyses. GE automatically visualizes genome assembly and mapping results output from Maser pipelines, without requiring additional data upload. With this function, the Maser pipelines can graphically display the results output from all the embedded tools and mapping results in a web browser. Therefore Maser realized a more user-friendly analysis platform especially for beginners by improving graphical display and providing the selected standard pipelines that work with built-in genome browser. In addition, all the analyses executed on Maser are recorded in the analysis history, helping users to trace and repeat the analyses. The entire process of analysis and its histories can be shared with collaborators or opened to the public. In conclusion, our system is useful for managing, analyzing, and visualizing NGS data and achieves traceability, reproducibility, and transparency of NGS analysis.

Database URL: http://cell-innovation.nig.ac.jp/maser/

Identification and prediction of alternative transcription start sites that generate rod photoreceptor-specific transcripts from ubiquitously expressed genes

Transcriptome complexity is substantially increased by the use of multiple transcription start sites for a given gene. By utilizing a rod photoreceptor-specific chromatin signature, and the RefSeq database of established transcription start sites, we have identified essentially all known rod photoreceptor genes as well as a group of novel genes that have a high probability of being expressed in rod photoreceptors. Approximately half of these novel rod genes are transcribed into multiple mRNA and/or protein isoforms through alternative transcriptional start sites (ATSS), only one of which has a rod-specific epigenetic signature and gives rise to a rod transcript. This suggests that, during retina development, some genes use ATSS to regulate cell type and temporal specificity, effectively generating a rod transcript from otherwise ubiquitously expressed genes. Biological confirmation of the relationship between epigenetic signatures and gene expression, as well as comparison of our genome-wide chromatin signature maps with available data sets for retina, namely a ChIP-on-Chip study of Polymerase-II (Pol-II) binding sites, ChIP-Seq studies for NRL- and CRX- binding sites and DHS (University of Washington data, available on UCSC mouse Genome Browser as a part of ENCODE project) fully support our hypothesis and together accurately identify and predict an array of new rod transcripts. The same approach was used to identify a number of TSS that are not currently in RefSeq. Biological conformation of the use of some of these TSS suggests that this method will be valuable for exploring the range of transcriptional complexity in many tissues. Comparison of mouse and human genome-wide data indicates that most of these alternate TSS appear to be present in both species, indicating that our approach can be useful for identification of regulatory regions that might play a role in human retinal disease.

Transcriptome profiling of the rat retina after optic nerve transection

Glaucoma is a group of eye diseases characterized by alterations in the contour of the optic nerve head (ONH), with corresponding visual field defects and progressive loss of retinal ganglion cells (RGCs). This progressive RGC death is considered to originate in axonal injury caused by compression of the axon bundles in the ONH. However, the molecular pathomechanisms of axonal injury-induced RGC death are not yet well understood. Here, we used RNA sequencing (RNA-seq) to examine transcriptome changes in rat retinas 2 days after optic nerve transection (ONT), and then used computational techniques to predict the resulting alterations in the transcriptional regulatory network. RNA-seq revealed 267 differentially expressed genes after ONT, 218 of which were annotated and 49 unannotated. We also identified differentially expressed transcripts, including potentially novel isoforms. An in silico pathway analysis predicted that CREB1 was the most significant upstream regulator. Thus, this study identified genes and pathways that may be involved in the pathomechanisms of axonal injury. We believe that our data should serve as a valuable resource to understand the molecular processes that define axonal injury-driven RGC death and to discover novel therapeutic targets for glaucoma.

Lens ER-stress response during cataract development in Mip-mutant mice

Major intrinsic protein (MIP) is a functional water-channel (AQP0) that also plays a key role in establishing lens fiber cell architecture. Genetic variants of MIP have been associated with inherited and age-related forms of cataract; however, the underlying pathogenic mechanisms are unclear. Here we have used lens transcriptome profiling by microarray-hybridization and qPCR to identify pathogenic changes during cataract development in Mip-mutant (Lop/+) mice. In postnatal Lop/+ lenses (P7) 99 genes were up-regulated and 75 were down-regulated (>2-fold, p=<0.05) when compared with wild-type. A pathway analysis of up-regulated genes in the Lop/+ lens (P7) was consistent with endoplasmic reticulum (ER)-stress and activation of the unfolded protein response (UPR). The most up-regulated UPR genes (>4-fold) in the Lop/+ lens included Chac1>Ddit3>Atf3>Trib3>Xbp1 and the most down-regulated genes (>5-fold) included two anti-oxidant genes, Hspb1 and Hmox1. Lop/+ lenses were further characterized by abundant TUNEL-positive nuclei within central degenerating fiber cells, glutathione depletion, free-radical overproduction, and calpain hyper-activation. These data suggest that Lop/+ lenses undergo proteotoxic ER-stress induced cell-death resulting from prolonged activation of the Eif2ak3/Perk-Atf4-Ddit3-Chac1 branch of the UPR coupled with severe oxidative-stress.

In vivo cellular imaging of various stress/response pathways using AAV following axonal injury in mice

Glaucoma, a leading cause of blindness worldwide, is instigated by various factors, including axonal injury, which eventually leads to a progressive loss of retinal ganglion cells (RGCs). To study various pathways reportedly involved in the pathogenesis of RGC death caused by axonal injury, seven pathways were investigated. Pathway-specific fluorescent protein-coded reporters were each packaged into an adeno-associated virus (AAV). After producing axonal injury in the eye, injected with AAV to induce RGC death, the temporal activity of each stress-related pathway was monitored in vivo through the detection of fluorescent RGCs using confocal ophthalmoscopy. We identified the activation of ATF6 and MCP-1 pathways involved in endoplasmic reticulum stress and macrophage recruitment, respectively, as early markers of RGC stress that precede neuronal death. Conversely, inflammatory responses probed by NF-κB and cell-death-related pathway p53 were most prominent in the later phases, when RGC death was already ongoing. AAV-mediated delivery of stress/response reporters followed by in vivo cellular imaging is a powerful strategy to characterize the temporal aspects of complex molecular pathways involved in retinal diseases. The identification of promoter elements that are activated before the death of RGCs enables the development of pre-emptive gene therapy, exclusively targeting the early phases of diseased cells.