Single-nucleus RNA sequencing of developing superior colliculus identifies neuronal diversity and candidate mediators of circuit assembly.

The superior colliculus (SC) is a sensorimotor structure in the midbrain that integrates input from multiple sensory modalities to initiate motor commands. It undergoes well-characterized steps of circuit assembly during development, rendering the mouse SC a popular model to study establishment of neural connectivity. Here we perform single-nucleus RNA-sequencing analysis of the mouse SC isolated at various developmental time points. Our study provides a transcriptomic landscape of the cell types that comprise the SC across murine development with particular emphasis on neuronal heterogeneity. We report a repertoire of genes differentially expressed across the different postnatal ages, many of which are known to regulate axon guidance and synapse formation. Using these data, we find that Pax7 expression is restricted to a subset of GABAergic neurons. Our data provide a valuable resource for interrogating the mechanisms of circuit development and identifying markers for manipulating specific SC neuronal populations and circuits.

Single-Nucleus RNA Sequencing of Developing and Mature Superior Colliculus Identifies Neuronal Diversity and Candidate Mediators of Circuit Assembly.

The superior colliculus (SC) is a sensorimotor structure in the midbrain that integrates input from multiple sensory modalities to initiate motor commands. It undergoes well-characterized steps of circuit assembly during development, rendering the mouse SC a popular model to study establishment and refinement of neural connectivity. Here we performed single nucleus RNA-sequencing analysis of the mouse SC isolated at various developmental time points. Our study provides a transcriptomic landscape of the cell types that comprise the SC across murine development with particular emphasis on neuronal heterogeneity. We used these data to identify Pax7 as a marker for an anatomically homogeneous population of GABAergic neurons. Lastly, we report a repertoire of genes differentially expressed across the different postnatal ages, many of which are known to regulate axon guidance and synapse formation. Our data provide a valuable resource for interrogating the mechanisms of circuit development, and identifying markers for manipulating specific SC neuronal populations and circuits.

Retinal ganglion cell expression of cytokine enhances occupancy of NG2 cell-derived astrocytes at the nerve injury site: Implication for axon regeneration.

Following injury in the central nervous system, a population of astrocytes occupy the lesion site, form glial bridges and facilitate axon regeneration. These astrocytes originate primarily from resident astrocytes or NG2+ oligodendrocyte progenitor cells. However, the extent to which these cell types give rise to the lesion-filling astrocytes, and whether the astrocytes derived from different cell types contribute similarly to optic nerve regeneration remain unclear. Here we examine the distribution of astrocytes and NG2+ cells in an optic nerve crush model. We show that optic nerve astrocytes partially fill the injury site over time after a crush injury. Viral mediated expression of a growth-promoting factor, ciliary neurotrophic factor (CNTF), in retinal ganglion cells (RGCs) promotes axon regeneration without altering the lesion size or the degree of lesion-filling GFAP+ cells. Strikingly, using inducible NG2CreER driver mice, we found that CNTF overexpression in RGCs increases the occupancy of NG2+ cell-derived astrocytes in the optic nerve lesion. An EdU pulse-chase experiment shows that the increase in NG2 cell-derived astrocytes is not due to an increase in cell proliferation. Lastly, we performed RNA-sequencing on the injured optic nerve and reveal that CNTF overexpression in RGCs results in significant changes in the expression of distinct genes, including those that encode chemokines, growth factor receptors, and immune cell modulators. Even though CNTF-induced axon regeneration has long been recognized, this is the first evidence of this procedure affecting glial cell fate at the optic nerve crush site. We discuss possible implication of these results for axon regeneration.

Identification of long noncoding RNAs in injury-resilient and injury-susceptible mouse retinal ganglion cells.

BACKGROUND: Emerging evidence indicates that long noncoding RNAs (lncRNAs) are important regulators of various biological processes, and their expression can be altered following certain pathological conditions, including central nervous system injury. Retinal ganglion cells (RGCs), whose axons form the optic nerve, are a heterogeneous population of neurons with more than 40 molecularly distinct subtypes in mouse. While most RGCs, including the ON-OFF direction-selective RGCs (ooDSGCs), are vulnerable to axonal injury, a small population of RGCs, including the intrinsically photosensitive RGCs (ipRGCs), are more resilient. RESULTS: By performing systematic analyses on RNA-sequencing data, here we identify lncRNAs that are expressed in ooDSGCs and ipRGCs with and without axonal injury. Our results reveal a repertoire of different classes of lncRNAs, including long intergenic noncoding RNAs and antisense ncRNAs that are differentially expressed between these RGC types. Strikingly, we also found dozens of lncRNAs whose expressions are altered markedly in response to axonal injury, some of which are expressed exclusively in either one of the types. Moreover, analyses into these lncRNAs unraveled their neighboring coding genes, many of which encode transcription factors and signaling molecules, suggesting that these lncRNAs may act in cis to regulate important biological processes in these neurons. Lastly, guilt-by-association analysis showed that lncRNAs are correlated with apoptosis associated genes, suggesting potential roles for these lncRNAs in RGC survival. CONCLUSIONS: Overall, the results of this study reveal RGC type-specific expression of lncRNAs and provide a foundation for future investigation of the function of lncRNAs in regulating neuronal type specification and survival.

Neural Cadherin Plays Distinct Roles for Neuronal Survival and Axon Growth under Different Regenerative Conditions.

Growing axons in the CNS often migrate along specific pathways to reach their targets. During embryonic development, this migration is guided by different types of cell adhesion molecules (CAMs) present on the surface of glial cells or other neurons, including the neural cadherin (NCAD). Axons in the adult CNS can be stimulated to regenerate, and travel long distances. Crucially, however, while a few axons are guided effectively through the injured nerve under certain conditions, most axons never migrate properly. The molecular underpinnings of the variable growth, and the glial CAMs that are responsible for CNS axon regeneration remain unclear. Here we used optic nerve crush to demonstrate that NCAD plays multifaceted functions in facilitating CNS axon regeneration. Astrocyte-specific deletion of NCAD dramatically decreases regeneration induced by phosphatase and tensin homolog (PTEN) ablation in retinal ganglion cells (RGCs). Consistent with NCAD's tendency to act as homodimers, deletion of NCAD in RGCs also reduces regeneration. Deletion of NCAD in astrocytes neither alters RGCs' mammalian target of rapamycin complex 1 (mTORC1) activity nor lesion size, two factors known to affect regeneration. Unexpectedly, however, we find that NCAD deletion in RGCs reduces PTEN-deletion-induced RGC survival. We further show that NCAD deletion, in either astrocytes or RGCs, has negligible effects on the regeneration induced by ciliary neurotrophic factor (CNTF), suggesting that other CAMs are critical under this regenerative condition. Consistent with this notion, CNTF induces expression various integrins known to mediate cell adhesion. Together, our study reveals multilayered functions of NCAD and a molecular basis of variability in guided axon growth.

Thrombospondin-1 Mediates Axon Regeneration in Retinal Ganglion Cells.

Neuronal subtypes show diverse injury responses, but the molecular underpinnings remain elusive. Using transgenic mice that allow reliable visualization of axonal fate, we demonstrate that intrinsically photosensitive retinal ganglion cells (ipRGCs) are both resilient to cell death and highly regenerative. Using RNA sequencing (RNA-seq), we show genes that are differentially expressed in ipRGCs and that associate with their survival and axon regeneration. Strikingly, thrombospondin-1 (Thbs1) ranked as the most differentially expressed gene, along with the well-documented injury-response genes Atf3 and Jun. THBS1 knockdown in RGCs eliminated axon regeneration. Conversely, RGC overexpression of THBS1 enhanced regeneration in both ipRGCs and non-ipRGCs, an effect that was dependent on syndecan-1, a known THBS1-binding protein. All structural domains of the THBS1 were not equally effective; the trimerization and C-terminal domains promoted regeneration, while the THBS type-1 repeats were dispensable. Our results identify cell-type-specific induction of Thbs1 as a novel gene conferring high regenerative capacity.

Chromatin Landscape Distinguishes the Genomic Loci of Hundreds of Androgen-Receptor-Associated LincRNAs From the Loci of Non-associated LincRNAs.

Cell signaling events triggered by androgen hormone in prostate cells is dependent on activation of the androgen receptor (AR) transcription factor. Androgen hormone binding to AR promotes its displacement from the cytoplasm to the nucleus and AR binding to DNA motifs, thus inducing activatory and inhibitory transcriptional programs through a complex regulatory mechanism not yet fully understood. In this work, we performed RNA-seq deep-sequencing of LNCaP prostate cancer cells and found over 7000 expressed long intergenic non-coding RNAs (lincRNAs), of which ∼4000 are novel lincRNAs, and 258 lincRNAs have their expression activated by androgen. Immunoprecipitation of AR, followed by large-scale sequencing of co-immunoprecipitated RNAs (RIP-Seq) has identified in the LNCaP cell line a total of 619 lincRNAs that were significantly enriched (FDR < 10%, DESeq2) in the anti-Androgen Receptor (antiAR) fraction in relation to the control fraction (non-specific IgG), and we named them Androgen-Receptor-Associated lincRNAs (ARA-lincRNAs). A genome-wide analysis showed that protein-coding gene neighbors to ARA-lincRNAs had a significantly higher androgen-induced change in expression than protein-coding genes neighboring lincRNAs not associated to AR. To find relevant epigenetic signatures enriched at the ARA-lincRNAs' transcription start sites (TSSs) we used a machine learning approach and identified that the ARA-lincRNA genomic loci in LNCaP cells are significantly enriched with epigenetic marks that are characteristic of in cis enhancer RNA regulators, and that the H3K27ac mark of active enhancers is conspicuously enriched at the TSS of ARA-lincRNAs adjacent to androgen-activated protein-coding genes. In addition, LNCaP topologically associating domains (TADs) that comprise chromatin regions with ARA-lincRNAs exhibit transcription factor contents, epigenetic marks and gene transcriptional activities that are significantly different from TADs not containing ARA-lincRNAs. This work highlights the possible involvement of hundreds of lincRNAs working in synergy with the AR on the genome-wide androgen-induced gene regulatory program in prostate cells.

Chromatin landscape distinguishes the genomic loci of hundreds of Androgen-Receptor-Associated LincRNAs from the loci of non-associated LincRNAs

Cell signaling events triggered by androgen hormone in prostate cells is dependent on activation of the androgen receptor (AR) transcription factor. Androgen hormone binding to AR promotes its displacement from the cytoplasm to the nucleus and AR binding to DNA motifs, thus inducing activatory and inhibitory transcriptional programs through a complex regulatory mechanism not yet fully understood. In this work, we performed RNA-seq deep-sequencing of LNCaP prostate cancer cells and found over 7000 expressed long intergenic non-coding RNAs (lincRNAs), of which similar to 4000 are novel lincRNAs, and 258 lincRNAs have their expression activated by androgen. Immunoprecipitation of AR, followed by large-scale sequencing of co-immunoprecipitated RNAs (RIP-Seq) has identified in the LNCaP cell line a total of 619 lincRNAs that were significantly enriched (FDR < 10%, DESeq2) in the anti-Androgen Receptor (antiAR) fraction in relation to the control fraction (non-specific IgG), and we named them Androgen-Receptor-Associated lincRNAs (ARA-lincRNAs). A genome-wide analysis showed that protein-coding gene neighbors to ARA-lincRNAs had a significantly higher androgen-induced change in expression than protein-coding genes neighboring lincRNAs not associated to AR. To find relevant epigenetic signatures enriched at the ARA-lincRNAs' transcription start sites (TSSs) we used a machine learning approach and identified that the ARA-lincRNA genomic loci in LNCaP cells are significantly enriched with epigenetic marks that are characteristic of in cis enhancer RNA regulators, and that the H3K27ac mark of active enhancers is conspicuously enriched at the TSS of ARA-lincRNAs adjacent to androgen-activated protein-coding genes. In addition, LNCaP topologically associating domains (TADs) that comprise chromatin regions with ARA-lincRNAs exhibit transcription factor contents, epigenetic marks and gene transcriptional activities that are significantly different from TADs not containing ARA-lincRNAs. This work highlights the possible involvement of hundreds of lincRNAs working in synergy with the AR on the genome-wide androgen-induced gene regulatory program in prostate cells.

Chromatin landscape distinguishes the genomic loci of hundreds of Androgen-Receptor-Associated LincRNAs from the loci of non-associated LincRNAs

Cell signaling events triggered by androgen hormone in prostate cells is dependent on activation of the androgen receptor (AR) transcription factor. Androgen hormone binding to AR promotes its displacement from the cytoplasm to the nucleus and AR binding to DNA motifs, thus inducing activatory and inhibitory transcriptional programs through a complex regulatory mechanism not yet fully understood. In this work, we performed RNA-seq deep-sequencing of LNCaP prostate cancer cells and found over 7000 expressed long intergenic non-coding RNAs (lincRNAs), of which similar to 4000 are novel lincRNAs, and 258 lincRNAs have their expression activated by androgen. Immunoprecipitation of AR, followed by large-scale sequencing of co-immunoprecipitated RNAs (RIP-Seq) has identified in the LNCaP cell line a total of 619 lincRNAs that were significantly enriched (FDR < 10%, DESeq2) in the anti-Androgen Receptor (antiAR) fraction in relation to the control fraction (non-specific IgG), and we named them Androgen-Receptor-Associated lincRNAs (ARA-lincRNAs). A genome-wide analysis showed that protein-coding gene neighbors to ARA-lincRNAs had a significantly higher androgen-induced change in expression than protein-coding genes neighboring lincRNAs not associated to AR. To find relevant epigenetic signatures enriched at the ARA-lincRNAs' transcription start sites (TSSs) we used a machine learning approach and identified that the ARA-lincRNA genomic loci in LNCaP cells are significantly enriched with epigenetic marks that are characteristic of in cis enhancer RNA regulators, and that the H3K27ac mark of active enhancers is conspicuously enriched at the TSS of ARA-lincRNAs adjacent to androgen-activated protein-coding genes. In addition, LNCaP topologically associating domains (TADs) that comprise chromatin regions with ARA-lincRNAs exhibit transcription factor contents, epigenetic marks and gene transcriptional activities that are significantly different from TADs not containing ARA-lincRNAs. This work highlights the possible involvement of hundreds of lincRNAs working in synergy with the AR on the genome-wide androgen-induced gene regulatory program in prostate cells.

Global analysis of biogenesis, stability and sub-cellular localization of lncRNAs mapping to intragenic regions of the human genome.

Long noncoding RNAs (lncRNAs) that map to intragenic regions of the human genome with the same (intronic lncRNAs) or opposite orientation (antisense lncRNAs) relative to protein-coding mRNAs have been largely dismissed from biochemical and functional characterization due to the belief that they are mRNA precursors, byproducts of RNA splicing or simply transcriptional noise. In this work, we used a custom microarray to investigate aspects of the biogenesis, processing, stability, evolutionary conservation, and cellular localization of ∼ 6,000 intronic lncRNAs and ∼ 10,000 antisense lncRNAs. Most intronic (2,903 of 3,427, 85%) and antisense lncRNAs (4,945 of 5,214, 95%) expressed in HeLa cells showed evidence of 5' cap modification, compatible with their transcription by RNAP II. Antisense lncRNAs (median t1/2 = 3.9 h) were significantly (p < 0.0001) more stable than mRNAs (median t1/2 = 3.2 h), whereas intronic lncRNAs (median t1/2 = 2.1 h) comprised a more heterogeneous class that included both stable (t1/2 > 3 h) and unstable (t1/2 < 1 h) transcripts. Intragenic lncRNAs display evidence of evolutionary conservation, have little/no coding potential and were ubiquitously detected in the cytoplasm. Notably, a fraction of the intronic and antisense lncRNAs (13 and 15%, respectively) were expressed from loci at which the corresponding host mRNA was not detected. The abundances of a subset of intronic/antisense lncRNAs were correlated (r ≥ |0.8|) with those of genes encoding proteins involved in cell division and DNA replication. Taken together, the findings of this study contribute novel biochemical and genomic information regarding intronic and antisense lncRNAs, supporting the notion that these classes include independently transcribed RNAs with potentials for exerting regulatory functions in the cell.

Long non-coding RNA INXS is a critical mediator of BCL-XS induced apoptosis.

BCL-X mRNA alternative splicing generates pro-apoptotic BCL-XS or anti-apoptotic BCL-XL gene products and the mechanism that regulates splice shifting is incompletely understood. We identified and characterized a long non-coding RNA (lncRNA) named INXS, transcribed from the opposite genomic strand of BCL-X, that was 5- to 9-fold less abundant in tumor cell lines from kidney, liver, breast and prostate and in kidney tumor tissues compared with non-tumors. INXS is an unspliced 1903 nt-long RNA, is transcribed by RNA polymerase II, 5'-capped, nuclear enriched and binds Sam68 splicing-modulator. Three apoptosis-inducing agents increased INXS lncRNA endogenous expression in the 786-O kidney tumor cell line, increased BCL-XS/BCL-XL mRNA ratio and activated caspases 3, 7 and 9. These effects were abrogated in the presence of INXS knockdown. Similarly, ectopic INXS overexpression caused a shift in splicing toward BCL-XS and activation of caspases, thus leading to apoptosis. BCL-XS protein accumulation was detected upon INXS overexpression. In a mouse xenograft model, intra-tumor injections of an INXS-expressing plasmid caused a marked reduction in tumor weight, and an increase in BCL-XS isoform, as determined in the excised tumors. We revealed an endogenous lncRNA that induces apoptosis, suggesting that INXS is a possible target to be explored in cancer therapies.

The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation.

The down-regulation of the tumor-suppressor gene RASSF1A has been shown to increase cell proliferation in several tumors. RASSF1A expression is regulated through epigenetic events involving the polycomb repressive complex 2 (PRC2); however, the molecular mechanisms modulating the recruitment of this epigenetic modifier to the RASSF1 locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand on the RASSF1 gene locus in several cell lines and tissues and binds PRC2. ANRASSF1 is transcribed through RNA polymerase II and is 5'-capped and polyadenylated; it exhibits nuclear localization and has a shorter half-life compared with other lncRNAs that bind PRC2. ANRASSF1 endogenous expression is higher in breast and prostate tumor cell lines compared with non-tumor, and an opposite pattern is observed for RASSF1A. ANRASSF1 ectopic overexpression reduces RASSF1A abundance and increases the proliferation of HeLa cells, whereas ANRASSF1 silencing causes the opposite effects. These changes in ANRASSF1 levels do not affect the RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase in both PRC2 occupancy and histone H3K27me3 repressive marks, specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression was detected on PRC2 occupancy and histone H3K27me3 at the promoter regions of RASSF1C and the four other neighboring genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrated that ANRASSF1 forms an RNA/DNA hybrid and recruits PRC2 to the RASSF1A promoter. Together, these results demonstrate a novel mechanism of epigenetic repression of the RASSF1A tumor suppressor gene involving antisense unspliced lncRNA, in which ANRASSF1 selectively represses the expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the human genome might contribute to a location-specific epigenetic modulation of genes.