A comparative analysis of rod bipolar cell transcriptomes identifies novel genes implicated in night vision

Temporal Regulation of Myopia and Inflammation-Associated Pathways in the Interphotoreceptor Retinoid-Binding Protein Knockout Mouse Model

PURPOSE

Myopia is a complex disorder with etiology involving an interplay between several genetic and environmental factors. Interphotoreceptor retinoid-binding protein (IRBP) is found in the subretinal space and is crucial in the visual cycle. The interphotoreceptor retinoid-binding protein knockout mouse (IRBP KO) was established as a model system to understand myopia and retinal degeneration. The current study investigated genes associated with myopia, retinal homeostasis, and inflammation in IRBP KO.

METHODS

RNA from retinas of congenic IRBP KO and wild-type C57BL/6J (WT) mice at postnatal day 5 (P5), P40, and P213 were subjected to digital droplet PCR (ddPCR) using a Bio-Rad automated droplet generator and QX200 reader. Target genes were selected based on genome-wide association studies, animal models, myopia studies, and other genes associated with retinal homeostasis and inflammation. HPRT, a housekeeping gene, was used for normalization. An average expression ratio (target/HPRT) and standard deviation (SD) were calculated. ANOVA assessed statistical significance, and a p < 0.05 was considered significant.

RESULTS

The ddPCR data analysis indicated that numerous myopia and inflammation-associated genes were differentially regulated in IRBP KO retinas with distinct temporal variation (upregulated at P5, decreased at P40, and no change at P213 relative to WT). C1qa, Gjd2, Sntb1, and Vsx2 emerged as top genetic candidate pathways. Compared with WT, immunoblotting analysis of C1qa showed no significant differences at P5 but significantly increased protein levels at P7 in IRBP KOs. Vsx2 remained unaltered at P5 and P7 in KO when compared with WT.

CONCLUSIONS

Data analysis indicated significant contributions from C1q, Gjd2, Sntb1, and Vsx2 genes in IRBP deficiency.

CNS/PNS proteoglycans functionalize neuronal and astrocyte niche microenvironments optimizing cellular activity by preserving membrane polarization dynamics, ionic microenvironments, ion fluxes, neuronal activation, and network neurotransductive capacity

Central and peripheral nervous system (CNS/PNS) proteoglycans (PGs) have diverse functional roles, this study examined how these control cellular behavior and tissue function. The CNS/PNS extracellular matrix (ECM) is a dynamic, responsive, highly interactive, space-filling, cell supportive, stabilizing structure maintaining tissue compartments, ionic microenvironments, and microgradients that regulate neuronal activity and maintain the neuron in an optimal ionic microenvironment. The CNS/PNS contains a high glycosaminoglycan content (60% hyaluronan, HA) and a diverse range of stabilizing PGs. Immobilization of HA in brain tissues by HA interactive hyalectan PGs preserves tissue hydration and neuronal activity, a paucity of HA in brain tissues results in a pro-convulsant epileptic phenotype. Diverse CS, KS, and HSPGs stabilize the blood-brain barrier and neurovascular unit, provide smart gel neurotransmitter neuron vesicle storage and delivery, organize the neuromuscular junction basement membrane, and provide motor neuron synaptic plasticity, and photoreceptor and neuron synaptic functions. PG-HA networks maintain ionic fluxes and microgradients and tissue compartments that contribute to membrane polarization dynamics essential to neuronal activation and neurotransduction. Hyalectans form neuroprotective perineuronal nets contributing to synaptic plasticity, memory, and cognitive learning. Sialoglycoprotein associated with cones and rods (SPACRCAN), an HA binding CSPG, stabilizes the inter-photoreceptor ECM. HSPGs pikachurin and eyes shut stabilize the photoreceptor synapse aiding in phototransduction and neurotransduction with retinal bipolar neurons crucial to visual acuity. This is achieved through Laminin G motifs in pikachurin, eyes shut, and neurexins that interact with the dystroglycan-cytoskeleton-ECM-stabilizing synaptic interconnections, neuronal interactive specificity, and co-ordination of regulatory action potentials in neural networks.

Specific photoreceptor cell fate pathways are differentially altered in NR2E3-associated diseases

Mutations in NR2E3, a gene encoding an orphan nuclear transcription factor, cause two retinal dystrophies with a distinct phenotype, but the precise role of NR2E3 in rod and cone transcriptional networks remains unclear. To dissect NR2E3 function, we performed scRNA-seq in the retinas of wildtype and two different Nr2e3 mouse models that show phenotypes similar to patients carrying NR2E3 mutations. Our results reveal that rod and cone populations are not homogeneous and can be separated into different sub-classes. We identify a previously unreported cone pathway that generates hybrid cones co-expressing both cone- and rod-related genes. In mutant retinas, this hybrid cone subpopulation is more abundant and includes a subpopulation of rods transitioning towards a cone cell fate. Hybrid photoreceptors with high misexpression of cone- and rod-related genes are prone to regulated necrosis. Overall, our results shed light on the role of NR2E3 in modulating photoreceptor differentiation towards cone and rod fates and explain how different mutations in NR2E3 lead to distinct visual disorders in humans.

Cross-species scRNA-seq reveals the cellular landscape of retina and early alterations in type 2 diabetes mice

Single-cell RNA sequencing (scRNA-seq) analysis have provided an unprecedented resolution for the studies on diabetic retinopathy (DR). However, the early changes in the retina in diabetes remain unclear. A total of 8 human and mouse scRNA-seq datasets, containing 276,402 cells were analyzed individually to comprehensively delineate the retinal cell atlas. The neural retinas were isolated from the type 2 diabetes (T2D) and control mice, and scRNA-seq analysis was conducted to evaluate the early effects of diabetes on the retina. Bipolar cell (BC) heterogeneity were identified. We found some stable BCs across multiple datasets, and explored their biological functions. A new RBC subtype (Car8_RBC) in the mouse retina was validated using the multi-color immunohistochemistry. AC149090.1 was significantly upregulated in the rod cells, ON cone BCs (CBCs), OFF CBCs, and RBCs in T2D mice. Additionally, the interneurons, especially BCs, were the most vulnerable cells to diabetes by integrating scRNA-seq and genome-wide association studies (GWAS) analyses. In conclusion, this study delineated a cross-species retinal cell atlas and uncovered the early pathological alterations in the retina of T2D mice.

Shedding light on myopia by studying complete congenital stationary night blindness

Myopia is the most common eye disorder, caused by heterogeneous genetic and environmental factors. Rare progressive and stationary inherited retinal disorders are often associated with high myopia. Genes implicated in myopia encode proteins involved in a variety of biological processes including eye morphogenesis, extracellular matrix organization, visual perception, circadian rhythms, and retinal signaling. Differentially expressed genes (DEGs) identified in animal models mimicking myopia are helpful in suggesting candidate genes implicated in human myopia. Complete congenital stationary night blindness (cCSNB) in humans and animal models represents an ON-bipolar cell signal transmission defect and is also associated with high myopia. Thus, it represents also an interesting model to identify myopia-related genes, as well as disease mechanisms. While the origin of night blindness is molecularly well established, further research is needed to elucidate the mechanisms of myopia development in subjects with cCSNB. Using whole transcriptome analysis on three different mouse models of cCSNB (in Gpr179^-/-, Lrit3^-/- and Grm6^-/-), we identified novel actors of the retinal signaling cascade, which are also novel candidate genes for myopia. Meta-analysis of our transcriptomic data with published transcriptomic databases and genome-wide association studies from myopia cases led us to propose new biological/cellular processes/mechanisms potentially at the origin of myopia in cCSNB subjects. The results provide a foundation to guide the development of pharmacological myopia therapies.

Transient expression of a GABA receptor subunit during early development is critical for inhibitory synapse maturation and function

Developing neural circuits, including GABAergic circuits, switch receptor types. But the role of early GABA receptor expression for establishment of functional inhibitory circuits remains unclear. Tracking the development of GABAergic synapses across axon terminals of retinal bipolar cells (BCs), we uncovered a crucial role of early GABA_A receptor expression for the formation and function of presynaptic inhibitory synapses. Specifically, early α3-subunit-containing GABA_A (GABA_Aα3) receptors are a key developmental organizer. Before eye opening, GABA_Aα3 gives way to GABA_Aα1 at individual BC presynaptic inhibitory synapses. The developmental downregulation of GABA_Aα3 is independent of GABA_Aα1 expression. Importantly, lack of early GABA_Aα3 impairs clustering of GABA_Aα1 and formation of functional GABA_A synapses across mature BC terminals. This impacts the sensitivity of visual responses transmitted through the circuit. Lack of early GABA_Aα3 also perturbs aggregation of LRRTM4, the organizing protein at GABAergic synapses of rod BC terminals, and their arrangement of output ribbon synapses.

Evolutionary Perspective and Expression Analysis of Intronless Genes Highlight the Conservation of Their Regulatory Role

The structure of eukaryotic genes is generally a combination of exons interrupted by intragenic non-coding DNA regions (introns) removed by RNA splicing to generate the mature mRNA. A fraction of genes, however, comprise a single coding exon with introns in their untranslated regions or are intronless genes (IGs), lacking introns entirely. The latter code for essential proteins involved in development, growth, and cell proliferation and their expression has been proposed to be highly specialized for neuro-specific functions and linked to cancer, neuropathies, and developmental disorders. The abundant presence of introns in eukaryotic genomes is pivotal for the precise control of gene expression. Notwithstanding, IGs exempting splicing events entail a higher transcriptional fidelity, making them even more valuable for regulatory roles. This work aimed to infer the functional role and evolutionary history of IGs centered on the mouse genome. IGs consist of a subgroup of genes with one exon including coding genes, non-coding genes, and pseudogenes, which conform approximately 6% of a total of 21,527 genes. To understand their prevalence, biological relevance, and evolution, we identified and studied 1,116 IG functional proteins validating their differential expression in transcriptomic data of embryonic mouse telencephalon. Our results showed that overall expression levels of IGs are lower than those of MEGs. However, strongly up-regulated IGs include transcription factors (TFs) such as the class 3 of POU (HMG Box), Neurog1, Olig1, and BHLHe22, BHLHe23, among other essential genes including the β-cluster of protocadherins. Most striking was the finding that IG-encoded BHLH TFs fit the criteria to be classified as microproteins. Finally, predicted protein orthologs in other six genomes confirmed high conservation of IGs associated with regulating neural processes and with chromatin organization and epigenetic regulation in Vertebrata. Moreover, this study highlights that IGs are essential modulators of regulatory processes, such as the Wnt signaling pathway and biological processes as pivotal as sensory organ developing at a transcriptional and post-translational level. Overall, our results suggest that IG proteins have specialized, prevalent, and unique biological roles and that functional divergence between IGs and MEGs is likely to be the result of specific evolutionary constraints.

LRRTM4 is a member of the transsynaptic complex between rod photoreceptors and bipolar cells

Leucine rich repeat transmembrane (LRRTM) proteins are synaptic adhesion molecules with roles in synapse formation and signaling. LRRTM4 transcripts were previously shown to be enriched in rod bipolar cells (BCs), secondary neurons of the retina that form synapses with rod photoreceptors. Using two different antibodies, LRRTM4 was found to reside primarily at rod BC dendritic tips, where it colocalized with the transduction channel protein, TRPM1. LRRTM4 was not detected at dendritic tips of ON-cone BCs. Following somatic knockout of LRRTM4 in BCs by subretinal injection and electroporation of CRISPR/Cas9, LRRTM4 was abolished or reduced in the dendritic tips of transfected cells. Knockout cells had a normal complement of TRPM1 at their dendritic tips, while GPR179 accumulation was partially reduced. In experiments with heterologously expressed protein, the extracellular domain of LRRTM4 was found to engage in heparan-sulfate dependent binding with pikachurin. These results implicate LRRTM4 in the GPR179-pikachurin-dystroglycan transsynaptic complex at rod synapses.

Expression and Localization of Kcne2 in the Vertebrate Retina

Purpose

To characterize the retinal expression and localization of Kcne2, an ancillary (β) ion-channel subunit with an important role in fine-tuning cellular excitability.

Methods

We analyzed available single-cell transcriptome data from tens of thousands of murine retinal cells for cell-type-specific expression of Kcne2 using state-of-the-art bioinformatics techniques. This evidence at the transcriptome level was complemented with a comprehensive immunohistochemical characterization of mouse retina (C57BL/6, ages 8-12 weeks) employing co-labeling techniques and cell-type-specific antibody markers. We furthermore examined how conserved the Kcne2 localization pattern in the retina was across species by performing immunostaining on zebrafish, cowbird, sheep, mice, and macaque.

Results

Kcne2 is distinctly expressed in cone photoreceptors and rod bipolar cells. At a subcellular level, the bulk of Kcne2 immunoreactivity can be observed in the outer plexiform layer. Here, it localizes into cone pedicles and likely the postsynaptic membrane of the rod bipolar cells. Thus, the vast majority of Kcne2 immunoreactivity is observed in a thin band in the outer plexiform layer. In addition to this, faint Kcne2 immunoreactivity can also be observed in cone inner segments and the somata of a small subset of cone ON bipolar cells. Strikingly, the localization of Kcne2 in the outer plexiform layer was preserved among all of the species studied, spanning at least 300 million years of evolution of the vertebrate kingdom.

Conclusions

The data we present here suggest an important and specific role for Kcne2 in the highly specialized photoreceptor-bipolar cell synapse.

LRRTM4: A Novel Regulator of Presynaptic Inhibition and Ribbon Synapse Arrangements of Retinal Bipolar Cells

LRRTM4 is a transsynaptic adhesion protein regulating glutamatergic synapse assembly on dendrites of central neurons. In the mouse retina, we find that LRRTM4 is enriched at GABAergic synapses on axon terminals of rod bipolar cells (RBCs). Knockout of LRRTM4 reduces RBC axonal GABA_A and GABA_C receptor clustering and disrupts presynaptic inhibition onto RBC terminals. LRRTM4 removal also perturbs the stereotyped output synapse arrangement at RBC terminals. Synaptic ribbons are normally apposed to two distinct postsynaptic "dyad" partners, but in the absence of LRRTM4, "monad" and "triad" arrangements are also formed. RBCs from retinas deficient in GABA release also demonstrate dyad mis-arrangements but maintain LRRTM4 expression, suggesting that defects in dyad organization in the LRRTM4 knockout could originate from reduced GABA receptor function. LRRTM4 is thus a key synapse organizing molecule at RBC terminals, where it regulates function of GABAergic synapses and assembly of RBC synaptic dyads.

Expression and distribution of trophoblast glycoprotein in the mouse retina

We recently identified the leucine-rich repeat (LRR) adhesion protein, trophoblast glycoprotein (TPBG), as a novel PKCα-dependent phosphoprotein in retinal rod bipolar cells (RBCs). Since TPBG has not been thoroughly examined in the retina, this study characterizes the localization and expression patterns of TPBG in the developing and adult mouse retina using two antibodies, one against the N-terminal LRR domain and the other against the C-terminal PDZ-interacting motif. Both antibodies labeled RBC dendrites in the outer plexiform layer and axon terminals in the IPL, as well as a putative amacrine cell with their cell bodies in the inner nuclear layer (INL) and a dense layer in the middle of the inner plexiform layer (IPL). In live transfected HEK293 cells, TPBG was localized to the plasma membrane with the N-terminal LRR domain facing the extracellular space. TPBG immunofluorescence in RBCs was strongly altered by the loss of TRPM1 in the adult retina, with significantly less dendritic and axon terminal labeling in TRPM1 knockout compared to wild type, despite no change in total TPBG detected by immunoblotting. During retinal development, TPBG expression increases dramatically just prior to eye opening with a time course closely correlated with that of TRPM1 expression. In the retina, LRR proteins have been implicated in the development and maintenance of functional bipolar cell synapses, and TPBG may play a similar role in RBCs.

An integrated transcriptional analysis of the developing human retina

The scarcity of embryonic/foetal material as a resource for direct study means that there is still limited understanding of human retina development. Here, we present an integrated transcriptome analysis combined with immunohistochemistry in human eye and retinal samples from 4 to 19 post-conception weeks. This analysis reveals three developmental windows with specific gene expression patterns that informed the sequential emergence of retinal cell types and enabled identification of stage-specific cellular and biological processes, and transcriptional regulators. Each stage is characterised by a specific set of alternatively spliced transcripts that code for proteins involved in the formation of the photoreceptor connecting cilium, pre-mRNA splicing and epigenetic modifiers. Importantly, our data show that the transition from foetal to adult retina is characterised by a large increase in the percentage of mutually exclusive exons that code for proteins involved in photoreceptor maintenance. The circular RNA population is also defined and shown to increase during retinal development. Collectively, these data increase our understanding of human retinal development and the pre-mRNA splicing process, and help to identify new candidate disease genes.