mRNA expression, splicing and editing in the embryonic and adult mouse cerebral cortex

Alternative splicing categorizes organ development by stage and reveals unique human splicing variants linked to neuromuscular disorders

Alternative splicing (AS) diversifies protein expression and contributes to species-specific differences in organ development. Here, we focused on stage-specific splicing variants and their correlation with disease in human compared to mouse during brain and heart development. Temporal transcriptomic analysis revealed that splicing factors (SFs) can accurately classify organ developmental stages, and 5 SFs were identified specifically upregulated in human during organogenesis. Additionally, inter-stage splicing variations were identified across analogous human and mouse developmental stages. Developmentally dynamic alternative splicing genes (devASGs) were enriched in various neurodevelopmental disorders in both species, with the most significant changes observed in human newborn brain and 16 weeks post-conception heart. Intriguingly, diseases specifically enriched in humans were primarily associated with neuro-muscular dysfunction, and human-specific neuromuscular devASGs were linked to mannose glycosylation and ciliary motility. These findings highlight the significance of SFs and AS events in organogenesis, and inform the selection of appropriate models for translational research.

Transcriptome-wide alternative mRNA splicing analysis reveals post-transcriptional regulation of neuronal differentiation

Alternative splicing (AS) plays an important role in neuronal development, function, and disease. Efforts to analyze the transcriptome of AS in neurons on a wide scale are currently limited. We characterized the transcriptome-wide AS changes in SH-SY5Y neuronal differentiation model, which is widely used to study neuronal function and disorders. Our analysis revealed global changes in five AS programs that drive neuronal differentiation. Motif analysis revealed the contribution of RNA-binding proteins (RBPs) to the regulation of AS during neuronal development. We concentrated on the primary alternative splicing program that occurs during differentiation, specifically on events involving exon skipping (SE). Motif analysis revealed motifs for polypyrimidine tract-binding protein 1 (PTB) and ELAV-like RNA binding protein 1 (HuR/ELAVL1) to be the top enriched in SE events, and their protein levels were downregulated after differentiation. shRNA knockdown of either PTB and HuR was associated with enhanced neuronal differentiation and transcriptome-wide exon skipping events that drive the process of differentiation. At the level of gene expression, we observed only modest changes, indicating predominant post-transcriptional effects of PTB and HuR. We also observed that both RBPs altered cellular responses to oxidative stress, in line with the differentiated phenotype observed after either gene knockdown. Our work characterizes the AS changes in a widely used and important model of neuronal development and neuroscience research and reveals intricate post-transcriptional regulation of neuronal differentiation.

Computational Comparison of Differential Splicing Tools for Targeted RNA Long-Amplicon Sequencing (rLAS)

RNA sequencing (RNA-Seq) is a powerful technique for the quantification of transcripts and the analysis of alternative splicing. Previously, our laboratory developed the targeted RNA long-amplicon sequencing (rLAS) method, which has the advantage of allowing deep analysis of targeted specific transcripts. The computational tools for analyzing RNA-Seq data have boosted alternative splicing research by detecting and quantifying splicing events. However, the performance of these splicing tools has not yet been investigated for rLAS. Here, we evaluated the performance of four splicing tools (MAJIQ, rMATS, MISO, and SplAdder) using samples with different types of known splicing events (exon-skipping, multiple-exon-skipping, alternative 5' splicing, and alternative 3' splicing). MAJIQ was able to detect all of the types of events, but it was unable to detect one of the exon-skipping events. On the other hand, rMATS was able to detect all of the exon-skipping events. However, rMATS failed to detect other types of events besides exon-skipping events. Both MISO and SplAdder were unable to detect any of the events. These results indicate that MAJIQ presents better performance for the different types of splicing events in rLAS and that rMATS shows better performance for exon-skipping splicing events.

Next-Generation Mapping of the ACINUS-Mediated Alternative Splicing Machinery and Its Regulation by O-glycosylation in Arabidopsis

Alternative splicing (AS) is a key mechanism of gene regulation, but the full repertoire of proteins involved and the regulatory mechanisms governing this process remain poorly understood. Using TurboID-based proximity labeling coupled with mass spectrometry (PL-MS), we comprehensively mapped the Arabidopsis AS machinery, focusing on the evolutionarily conserved splicing factor ACINUS, its paralog PININ, and the stable interactor SR45. We identified 298 high-confidence components, including both established and novel interactors, providing strong evidence that alternative splicing is coupled to transcription and that multiple RNA processing steps occur simultaneously in plants. Bioinformatic analysis reveals high redundancy, conserved mechanisms, and unique plant-specific features. Selected known and novel interactors were validated by AS readouts and phenotypic analysis, which also revealed a coordinated influence on splicing. Furthermore, a systematic evaluation of O-glycosylation double mutants revealed that SECRET AGENT (O-GlcNAc transferase) and SPINDLY (O-fucose transferase) modulate AS through both ACINUS-dependent and -independent pathways. Our results reveal the conserved as well as plant-specific AS regulatory network and highlight the global role of sugar modification in RNA processing.

Cell-Type-Specific Splicing of Transcription Regulators and Ptbp1 by Rbfox1/2/3 in the Developing Neocortex

How master splicing regulators cross talk with each other and to what extent transcription regulators are differentially spliced remain unclear in the developing brain. Here, cell-type-specific RNA-Seq analyses of the developing neocortex uncover variable expression of the Rbfox1/2/3 genes and enriched alternative splicing events in transcription regulators, altering protein isoforms or inducing nonsense-mediated mRNA decay. Transient expression of Rbfox proteins in radial glial progenitors induces neuronal splicing events preferentially in transcription regulators such as Meis2 and Tead1 Surprisingly, Rbfox proteins promote the inclusion of a mammal-specific alternative exon and a previously undescribed poison exon in Ptbp1 Simultaneous ablation of Rbfox1/2/3 in the neocortex downregulates neuronal isoforms and disrupts radial neuronal migration. Furthermore, the progenitor isoform of Meis2 promotes Tgfb3 transcription, while the Meis2 neuron isoform promotes neuronal differentiation. These observations indicate that transcription regulators are differentially spliced between cell types in the developing neocortex. (The sex has not been reported to affect cortical neurogenesis in mice, and embryos of both sexes were studied without distinguishing one or the other.).

Alternative splicing of modulatory immune receptors in T lymphocytes: a newly identified and targetable mechanism for anticancer immunotherapy

Alternative splicing (AS) is a mechanism that generates translational diversity within a genome. Equally important is the dynamic adaptability of the splicing machinery, which can give preference to one isoform over others encoded by a single gene. These isoform preferences change in response to the cell's state and function. Particularly significant is the impact of physiological alternative splicing in T lymphocytes, where specific isoforms can enhance or reduce the cells' reactivity to stimuli. This process makes splicing isoforms defining features of cell states, exemplified by CD45 splice isoforms, which characterize the transition from naïve to memory states. Two developments have accelerated the use of AS dynamics for therapeutic interventions: advancements in long-read RNA sequencing and progress in nucleic acid chemical modifications. Improved oligonucleotide stability has enabled their use in directing splicing to specific sites or modifying sequences to enhance or silence particular splicing events. This review highlights immune regulatory splicing patterns with potential significance for enhancing anticancer immunotherapy.

Developmental correspondence of juvenile stages across the locust, harlequin ladybird, and diamondback moth

Insect metamorphosis is a captivating aspect of animal research. To address the controversy regarding the developmental correspondence between hemimetabolous and holometabolous insects, we utilized non-destructive micro-computed tomography (CT) imaging and RNA-seq to examine wing growth and transcriptome profiles across juvenile stages in the Locusta migratoria, Harmonia axyridis, and Plutella xylostella, with distinct metamorphic strategies. Micro-CT revealed that over 88% of wing volume increase in ladybirds and moths occurs during the prepupal-pupal transition, similar to locust nymphs. Developmental transcriptome clustering demonstrated that gene expression patterns more closely resembled those between ladybird/moth prepupae/pupae and locust nymphs, whereas holometabolous larvae exhibited distinct profiles. Notably, gene expression specificity increased across juvenile stages with more recent species divergence. Genes highly expressed around the prepupal/pupal stages accumulated higher evolutionary rates. These integrated findings unveil commonalities in juvenile stage development among the locust, ladybird, and moth, providing insights into the evolution of metamorphosis within neopteran insects.

Cell-type-specific splicing of transcription regulators and Ptbp1 by Rbfox1/2/3 in the developing neocortex

How master splicing regulators crosstalk with each other and to what extent transcription regulators are differentially spliced remain unclear in the developing brain. Here, cell-type-specific RNA-Seq of the developing neocortex uncover that transcription regulators are enriched for differential splicing, altering protein isoforms or inducing nonsense-mediated mRNA decay. Transient expression of Rbfox proteins in radial glia progenitors induces neuronal splicing events preferentially in transcription regulators such as Meis2 and Tead1. Surprisingly, Rbfox proteins promote the inclusion of a mammal-specific alternative exon and a previously undescribed poison exon in Ptbp1. Simultaneous ablation of Rbfox1/2/3 in the neocortex downregulates neuronal isoforms and disrupts radial neuronal migration. Furthermore, the progenitor isoform of Meis2 promotes Tgfb3 transcription, while the Meis2 neuron isoform promotes neuronal differentiation. These observations indicate that transcription regulators are differentially spliced between cell types in the developing neocortex.

Splicing is dynamically regulated during limb development

Modifications to highly conserved developmental gene regulatory networks are thought to underlie morphological diversification in evolution and contribute to human congenital malformations. Relationships between gene expression and morphology have been extensively investigated in the limb, where most of the evidence for alterations to gene regulation in development consists of pre-transcriptional mechanisms that affect expression levels, such as epigenetic alterations to regulatory sequences and changes to cis-regulatory elements. Here we report evidence that alternative splicing (AS), a post-transcriptional process that modifies and diversifies mRNA transcripts, is dynamic during limb development in two mammalian species. We evaluated AS patterns in mouse (Mus musculus) and opossum (Monodelphis domestica) across the three key limb developmental stages: the ridge, bud, and paddle. Our data show that splicing patterns are dynamic over developmental time and suggest differences between the two mammalian taxa. Additionally, multiple key limb development genes, including Fgf8, are differentially spliced across the three stages in both species, with expression levels of the conserved splice variants, Fgf8a and Fgf8b, changing across developmental time. Our data demonstrates that AS is a critical mediator of mRNA diversity in limb development and provides an additional mechanism for evolutionary tweaking of gene dosage.

SOX2 interacts with hnRNPK to modulate alternative splicing in mouse embryonic stem cells

BACKGROUND

SOX2 is a determinant transcription factor that governs the balance between stemness and differentiation by influencing transcription and splicing programs. The role of SOX2 is intricately shaped by its interactions with specific partners. In the interactome of SOX2 in mouse embryonic stem cells (mESCs), there is a cohort of heterogeneous nuclear ribonucleoproteins (hnRNPs) that contributes to multiple facets of gene expression regulation. However, the cross-talk between hnRNPs and SOX2 in gene expression regulation remains unclear.

RESULTS

Here we demonstrate the indispensable role of the co-existence of SOX2 and heterogeneous nuclear ribonucleoprotein K (hnRNPK) in the maintenance of pluripotency in mESCs. While hnRNPK directly interacts with the SOX2-HMG DNA-binding domain and induces the collapse of the transcriptional repressor 7SK small nuclear ribonucleoprotein (7SK snRNP), hnRNPK does not influence SOX2-mediated transcription, either by modulating the interaction between SOX2 and its target cis-regulatory elements or by facilitating transcription elongation as indicated by the RNA-seq analysis. Notably, hnRNPK enhances the interaction of SOX2 with target pre-mRNAs and collaborates with SOX2 in regulating the alternative splicing of a subset of pluripotency genes.

CONCLUSIONS

These data reveal that SOX2 and hnRNPK have a direct protein-protein interaction, and shed light on the molecular mechanisms by which hnRNPK collaborates with SOX2 in alternative splicing in mESCs.

Loss-of-function mutation in PRMT9 causes abnormal synapse development by dysregulation of RNA alternative splicing

Protein arginine methyltransferase 9 (PRMT9) is a recently identified member of the PRMT family, yet its biological function remains largely unknown. Here, by characterizing an intellectual disability associated PRMT9 mutation (G189R) and establishing a Prmt9 conditional knockout (cKO) mouse model, we uncover an important function of PRMT9 in neuronal development. The G189R mutation abolishes PRMT9 methyltransferase activity and reduces its protein stability. Knockout of Prmt9 in hippocampal neurons causes alternative splicing of ~1900 genes, which likely accounts for the aberrant synapse development and impaired learning and memory in the Prmt9 cKO mice. Mechanistically, we discover a methylation-sensitive protein-RNA interaction between the arginine 508 (R508) of the splicing factor 3B subunit 2 (SF3B2), the site that is exclusively methylated by PRMT9, and the pre-mRNA anchoring site, a cis-regulatory element that is critical for RNA splicing. Additionally, using human and mouse cell lines, as well as an SF3B2 arginine methylation-deficient mouse model, we provide strong evidence that SF3B2 is the primary methylation substrate of PRMT9, thus highlighting the conserved function of the PRMT9/SF3B2 axis in regulating pre-mRNA splicing.

Target-capture full-length double-stranded cDNA long-read sequencing through Nanopore revealed novel intron retention in patient with tuberous sclerosis complex

Tuberous sclerosis complex (TSC) is a relatively common autosomal dominant disorder characterized by multiple dysplastic organ lesions and neuropsychiatric symptoms caused by loss-of-function mutation of either TSC1 or TSC2. The genetic diagnosis of inherited diseases, including TSC, in the clinical field is widespread using next-generation sequencing. The mutations in protein-coding exon tend to be verified because mutations directly cause abnormal protein. However, it is relatively difficult to verify mutations in the intron region because it is required to investigate whether the intron mutations affect the abnormal splicing of transcripts. In this study, we developed a target-capture full-length double-stranded cDNA sequencing method using Nanopore long-read sequencer (Nanopore long-read target sequencing). This method revealed the occurrence of intron mutation in the TSC2 gene and found that the intron mutation produces novel intron retention splicing transcripts that generate truncated proteins. The protein-coding transcripts were decreased due to the expression of the novel intron retention transcripts, which caused TSC in patients with the intron mutation. Our results indicate that Nanopore long-read target sequencing is useful for the detection of mutations and confers information on the full-length alternative splicing of transcripts for genetic diagnosis.

Bi-allelic SNAPC4 variants dysregulate global alternative splicing and lead to neuroregression and progressive spastic paraparesis

The vast majority of human genes encode multiple isoforms through alternative splicing, and the temporal and spatial regulation of those isoforms is critical for organismal development and function. The spliceosome, which regulates and executes splicing reactions, is primarily composed of small nuclear ribonucleoproteins (snRNPs) that consist of small nuclear RNAs (snRNAs) and protein subunits. snRNA gene transcription is initiated by the snRNA-activating protein complex (SNAPc). Here, we report ten individuals, from eight families, with bi-allelic, deleterious SNAPC4 variants. SNAPC4 encoded one of the five SNAPc subunits that is critical for DNA binding. Most affected individuals presented with delayed motor development and developmental regression after the first year of life, followed by progressive spasticity that led to gait alterations, paraparesis, and oromotor dysfunction. Most individuals had cerebral, cerebellar, or basal ganglia volume loss by brain MRI. In the available cells from affected individuals, SNAPC4 abundance was decreased compared to unaffected controls, suggesting that the bi-allelic variants affect SNAPC4 accumulation. The depletion of SNAPC4 levels in HeLa cell lines via genomic editing led to decreased snRNA expression and global dysregulation of alternative splicing. Analysis of available fibroblasts from affected individuals showed decreased snRNA expression and global dysregulation of alternative splicing compared to unaffected cells. Altogether, these data suggest that these bi-allelic SNAPC4 variants result in loss of function and underlie the neuroregression and progressive spasticity in these affected individuals.

Programmable macromolecule-based RNA-targeting therapies to treat human neurological disorders

Disruptions in RNA processing play critical roles in the pathogenesis of neurological diseases. In this Perspective, we discuss recent progress in the development of RNA-targeting therapeutic modalities. We focus on progress, limitations, and opportunities in a new generation of therapies engineered from RNA binding proteins and other endogenous RNA regulatory macromolecules to treat human neurological disorders.

Post-Transcriptional Modification by Alternative Splicing and Pathogenic Splicing Variants in Cardiovascular Development and Congenital Heart Defects

Advancements in genomics, bioinformatics, and genome editing have uncovered new dimensions in gene regulation. Post-transcriptional modifications by the alternative splicing of mRNA transcripts are critical regulatory mechanisms of mammalian gene expression. In the heart, there is an expanding interest in elucidating the role of alternative splicing in transcriptome regulation. Substantial efforts were directed toward investigating this process in heart development and failure. However, few studies shed light on alternative splicing products and their dysregulation in congenital heart defects (CHDs). While elegant reports showed the crucial roles of RNA binding proteins (RBPs) in orchestrating splicing transitions during heart development and failure, the impact of RBPs dysregulation or genetic variation on CHDs has not been fully addressed. Herein, we review the current understanding of alternative splicing and RBPs' roles in heart development and CHDs. Wediscuss the impact of perinatal splicing transition and its dysregulation in CHDs. We further summarize the discoveries made of causal splicing variants in key transcription factors that are implicated in CHDs. An improved understanding of the roles of alternative splicing in heart development and CHDs may potentially inform novel preventive and therapeutic advancements for newborn infants with CHDs.

Transcriptomic analysis of the upper lip and primary palate development in mice

Background: Normal fusion of the upper lip and primary palate is a complex process involving a series of characteristic and orderly regulated cellular events. Cleft lip with or without palate (CL/P), one of the most common congenital malformations, may be induced by abnormalities in any of these events. However, less is known about the precise regulatory process in the fusion of the upper lip and primary palate. Methods: Lambdoidal junction tissues of mice from embryonic days 10.5, 11.5, and 12.5- three key fusion stages-were acquired for RNA sequencing. Results: Gene expression profiles in distinct fusion stages of mice were identified. Some of the differentially expressed genes (DEGs) have been reported to affect upper lip and primary palate development. However, other DEGs, such as Krt5, Pax1, Ambn, Hey2, and Tnmd, have not previously been investigated. Gene set enrichment analysis (GSEA) of these DEGs revealed the sequential intensification of Wnt, PI3K-Akt, MAPK, Hippo, and TGF-beta signaling pathways and identified relatively highly expressed genes including Tnn, Wnt3a, and Wnt16. We also observed substantial alternative splicing events during the fusion of the upper lip and primary palate and identified potentially important genes including Gtpbp8, Armcx1, Tle3, and Numa1. Protein-protein interaction (PPI) network analysis identified a series of hub genes, including Col1a2, Fos, Bmp2, Shh, Col1a1, Wnt3a, Anxa1, Gem, etc. Conclusion: Overall, the results of this study provided a comprehensive analysis of the development of the upper lip and primary palate. Our work provides insight into future studies of normal upper lip and primary palate development and the etiology of CL/P.

RNA epitranscriptomics dysregulation: A major determinant for significantly increased risk of ASD pathogenesis

Autism spectrum disorders (ASDs) are perhaps the most severe, intractable and challenging child psychiatric disorders. They are complex, pervasive and highly heterogeneous and depend on multifactorial neurodevelopmental conditions. Although the pathogenesis of autism remains unclear, it revolves around altered neurodevelopmental patterns and their implications for brain function, although these cannot be specifically linked to symptoms. While these affect neuronal migration and connectivity, little is known about the processes that lead to the disruption of specific laminar excitatory and inhibitory cortical circuits, a key feature of ASD. It is evident that ASD has multiple underlying causes and this multigenic condition has been considered to also dependent on epigenetic effects, although the exact nature of the factors that could be involved remains unclear. However, besides the possibility for differential epigenetic markings directly affecting the relative expression levels of individual genes or groups of genes, there are at least three mRNA epitranscriptomic mechanisms, which function cooperatively and could, in association with both genotypes and environmental conditions, alter spatiotemporal proteins expression patterns during brain development, at both quantitative and qualitative levels, in a tissue-specific, and context-dependent manner. As we have already postulated, sudden changes in environmental conditions, such as those conferred by maternal inflammation/immune activation, influence RNA epitranscriptomic mechanisms, with the combination of these processes altering fetal brain development. Herein, we explore the postulate whereby, in ASD pathogenesis, RNA epitranscriptomics might take precedence over epigenetic modifications. RNA epitranscriptomics affects real-time differential expression of receptor and channel proteins isoforms, playing a prominent role in central nervous system (CNS) development and functions, but also RNAi which, in turn, impact the spatiotemporal expression of receptors, channels and regulatory proteins irrespective of isoforms. Slight dysregulations in few early components of brain development, could, depending upon their extent, snowball into a huge variety of pathological cerebral alterations a few years after birth. This may very well explain the enormous genetic, neuropathological and symptomatic heterogeneities that are systematically associated with ASD and psychiatric disorders at large.

Derivation of splice junction-specific antibodies using a unique hapten targeting strategy and directed evolution

Alternative splicing of RNA occurs frequently in eukaryotic cells and can result in multiple protein isoforms that are nearly identical in amino acid sequence, but have unique biological roles. Moreover, the relative abundance of these unique isoforms can be correlative with diseased states and potentially used as biomarkers or therapeutic targets. However, due to high sequence similarities among isoforms, current proteomic methods are incapable of differentiating native protein isoforms derived from most alternative splicing events. Herein, a strategy employing a nonsynonymous, non-native amino acid (nnAA) pseudo-hapten (i.e. an amino acid or amino acid derivative that is different from the native amino acid at a particular position) as a targeting epitope in splice junction-spanning peptides was successful in directed antibody derivation. After isolating nnAA-specific antibodies, directed evolution reduced the antibody's binding dependence on the nnAA pseudo-hapten and improved binding to the native splice junction epitope. The resulting antibodies demonstrated codependent binding affinity to each exon of the splice junction and thus are splice junction- and isoform-specific. Furthermore, epitope scanning demonstrated that positioning of the nnAA pseudo-hapten within a peptide antigen can be exploited to predetermine the isolated antibody's specificity at, or near, amino acid resolution. Thus, this nnAA targeting strategy has the potential to robustly derive splice junction- and site-specific antibodies that can be used in a wide variety of research endeavors to unambiguously differentiate native protein isoforms.

Truncated jarid2 and kdm6b transcripts are associated with temperature-induced sex reversal during development in a dragon lizard

Sex determination and differentiation in reptiles are complex. In the model species, Pogona vitticeps, high incubation temperature can cause male to female sex reversal. To elucidate the epigenetic mechanisms of thermolabile sex, we used an unbiased genome-wide assessment of intron retention during sex reversal. The previously implicated chromatin modifiers (jarid2 and kdm6b) were two of three genes to display sex reversal-specific intron retention. In these species, embryonic intron retention resulting in C-terminally truncated jarid2 and kdm6b isoforms consistently occurs at low temperatures. High-temperature sex reversal is uniquely characterized by a high prevalence of N-terminally truncated isoforms of jarid2 and kdm6b, which are not present at low temperatures, or in two other reptiles with temperature-dependent sex determination. This work verifies that chromatin-modifying genes are involved in highly conserved temperature responses and can also be transcribed into isoforms with new sex-determining roles.

Nociception and hypersensitivity involve distinct neurons and molecular transducers in Drosophila

SignificanceFunctional plasticity of the nociceptive circuit is a remarkable feature and is of clinical relevance. As an example, nociceptors lower their threshold upon tissue injury, a process known as allodynia that would facilitate healing by guarding the injured areas. However, long-lasting hypersensitivity could lead to chronic pain, a debilitating disease not effectively treated. Therefore, it is crucial to dissect the mechanisms underlying basal nociception and nociceptive hypersensitivity. In both vertebrate and invertebrate species, conserved transient receptor potential (Trp) channels are the primary transducers of noxious stimuli. Here, we provide a precedent that in Drosophila larvae, heat sensing in the nociception and hypersensitivity states is mediated by distinct heat-sensitive neurons and TrpA1 alternative isoforms.

Landscape of adenosine-to-inosine RNA recoding across human tissues

RNA editing by adenosine deaminases changes the information encoded in the mRNA from its genomic blueprint. Editing of protein-coding sequences can introduce novel, functionally distinct, protein isoforms and diversify the proteome. The functional importance of a few recoding sites has been appreciated for decades. However, systematic methods to uncover these sites perform poorly, and the full repertoire of recoding in human and other mammals is unknown. Here we present a new detection approach, and analyze 9125 GTEx RNA-seq samples, to produce a highly-accurate atlas of 1517 editing sites within the coding region and their editing levels across human tissues. Single-cell RNA-seq data shows protein recoding contributes to the variability across cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation. This comprehensive set can facilitate understanding of the role of recoding in human physiology and diseases.

Gabay et al. provide a highly-accurate atlas of recoding by A-to-I RNA editing in human, profiled across tissues and cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation.

NOVA2 regulates neural circRNA biogenesis

Circular RNAs (circRNAs) are highly expressed in the brain and their expression increases during neuronal differentiation. The factors regulating circRNAs in the developing mouse brain are unknown. NOVA1 and NOVA2 are neural-enriched RNA-binding proteins with well-characterized roles in alternative splicing. Profiling of circRNAs from RNA-seq data revealed that global circRNA levels were reduced in embryonic cortex of Nova2 but not Nova1 knockout mice. Analysis of isolated inhibitory and excitatory cortical neurons lacking NOVA2 revealed an even more dramatic reduction of circRNAs and establishes a widespread role for NOVA2 in enhancing circRNA biogenesis. To investigate the cis-elements controlling NOVA2-regulation of circRNA biogenesis, we generated a backsplicing reporter based on the Efnb2 gene. We found that NOVA2-mediated backsplicing of circEfnb2 was impaired when YCAY clusters located in flanking introns were mutagenized. CLIP (cross-linking and immunoprecipitation) and additional reporter analyses demonstrated the importance of NOVA2 binding sites located in both flanking introns of circRNA loci. NOVA2 is the first RNA-binding protein identified to globally promote circRNA biogenesis in the developing brain.

Functional analysis of CASK transcript variants expressed in human brain

The calcium-/calmodulin dependent serine protein kinase (CASK) belongs to the membrane-associated guanylate kinases (MAGUK) family of proteins. It fulfils several different cellular functions, ranging from acting as a scaffold protein to transcription control, as well as regulation of receptor sorting. CASK functions depend on the interaction with a variety of partners, for example neurexin, liprin-α, Tbr1 and SAP97. So far, it is uncertain how these seemingly unrelated interactions and resulting functions of CASK are regulated. Here, we show that alternative splicing of CASK can guide the binding affinity of CASK isoforms to distinct interaction partners. We report seven different variants of CASK expressed in the fetal human brain. Four out of these variants are not present in the NCBI GenBank database as known human variants. Functional analyses showed that alternative splicing affected the affinities of CASK variants for several of the tested interaction partners. Thus, we observed a clear correlation of the presence of one splice insert with poor binding of CASK to SAP97, supported by molecular modelling. The alternative splicing and distinct properties of CASK variants in terms of protein-protein interaction should be taken into consideration for future studies.

A porcine brain-wide RNA editing landscape

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is an essential post-transcriptional modification. Although hundreds of thousands of RNA editing sites have been reported in mammals, brain-wide analysis of the RNA editing in the mammalian brain remains rare. Here, a genome-wide RNA-editing investigation is performed in 119 samples, representing 30 anatomically defined subregions in the pig brain. We identify a total of 682,037 A-to-I RNA editing sites of which 97% are not identified before. Within the pig brain, cerebellum and olfactory bulb are regions with most edited transcripts. The editing level of sites residing in protein-coding regions are similar across brain regions, whereas region-distinct editing is observed in repetitive sequences. Highly edited conserved recoding events in pig and human brain are found in neurotransmitter receptors, demonstrating the evolutionary importance of RNA editing in neurotransmission functions. Although potential data biases caused by age, sex or health status are not considered, this study provides a rich resource to better understand the evolutionary importance of post-transcriptional RNA editing.

Alternative splicing in Alzheimer's disease

Alzheimer's disease (AD) is the most frequent neurodegenerative disorder in the elderly, occurring in approximately 20% of people older than 80. The molecular causes of AD are still poorly understood. However, recent studies have shown that Alternative Splicing (AS) is involved in the gene expression reprogramming associated with the functional changes observed in AD patients. In particular, mutations in cis-acting regulatory sequences as well as alterations in the activity and sub-cellular localization of trans-acting splicing factors and components of the spliceosome machinery are associated with splicing abnormalities in AD tissues, which may influence the onset and progression of the disease. In this review, we discuss the current molecular understanding of how alterations in the AS process contribute to AD pathogenesis. Finally, recent therapeutic approaches targeting aberrant AS regulation in AD are also reviewed.

Alternative splicing of clathrin heavy chain contributes to the switch from coated pits to plaques

Clathrin function directly derives from its coat structure, and while endocytosis is mediated by clathrin-coated pits, large plaques contribute to cell adhesion. Here, we show that the alternative splicing of a single exon of the clathrin heavy chain gene (CLTC exon 31) helps determine the clathrin coat organization. Direct genetic control was demonstrated by forced CLTC exon 31 skipping in muscle cells that reverses the plasma membrane content from clathrin plaques to pits and by promoting exon inclusion that stimulated flat plaque assembly. Interestingly, mis-splicing of CLTC exon 31 found in the severe congenital form of myotonic dystrophy was associated with reduced plaques in patient myotubes. Moreover, forced exclusion of this exon in WT mice muscle induced structural disorganization and reduced force, highlighting the contribution of this splicing event for the maintenance of tissue homeostasis. This genetic control on clathrin assembly should influence the way we consider how plasticity in clathrin-coated structures is involved in muscle development and maintenance.

Rab27a-Dependent Paracrine Communication Controls Dendritic Spine Formation and Sensory Responses in the Barrel Cortex

Rab27a is an evolutionarily conserved small GTPase that regulates vesicle trafficking, and copy number variants of RAB27a are associated with increased risk of autism. However, the function of Rab27a on brain development is unknown. Here, we identified a form of paracrine communication that regulates spine development between distinct populations of developing cortical neurons. In the developing somatosensory cortex of mice, we show that decreasing Rab27a levels in late-born pyramidal neurons destined for layer (L) 2/3 had no cell-autonomous effect on their synaptic integration but increased excitatory synaptic transmission onto L4 neurons that receive somatosensory information. This effect resulted in an increased number of L4 neurons activated by whisker stimulation in juvenile mice. In addition, we found that Rab27a, the level of which decreases as neurons mature, regulates the release of small extracellular vesicles (sEVs) in developing neurons in vitro and decreasing Rab27a levels led to the accumulation of CD63-positive vesicular compartments in L2/3 neurons in vivo. Together, our study reveals that Rab27a-mediated paracrine communication regulates the development of synaptic connectivity, ultimately tuning responses to sensory stimulation, possibly via controlling the release of sEVs.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

Target-capture full-length double-strand cDNA sequencing for alternative splicing analysis

Alternative splicing is a regulated process by which eukaryotic genes may produce diverse biological products. Defects in the process typically affect cellular function and can lead to disease. Next-generation sequencing (NGS) technologies have been developed to detect alternative splicing events; however, the alternative splicing events detected by standard RNA-Seq may or may not be derived from full-length RNA. The SMARTer method provides full-length double-strand cDNA synthesis, and the resulting gene expression patterns correlate strongly with standard RNA-Seq. However, it also yields non-specific genomic DNA amplification. We improved the SMARTer method by employing a target-capture full-length double-strand cDNA sequencing method. High-fidelity, full-length cDNA is generated by the SMARTer method, followed by target-specific capture with exon probes. The expression pattern observed with this SMARTer Capture method was highly correlated with the results of the original SMARTer method. The number and accuracy of the detected splicing events were increased by eliminating non-specific genomic DNA amplification by the SMARTer Capture. Compared to the original SMARTer method, the SMARTer Capture provided 4-fold greater detection of alternative splicing events at the same read number, and it took less than 1/100 of read number to detect the same number of splicing events. The percent splicing in index (PSI) of the SMARTer Capture is highly correlated with the PSI of the SMARTer. These results indicate that the SMARTer Capture represents an improvement of the SMARTer method to accurately characterize alternative splicing repertories in targeted genes without biases.

Scaling of gene transcriptional gradients with brain size across mouse development

The structure of the adult brain is the result of complex physical mechanisms acting in three-dimensional space through development. Consequently, the brain's spatial embedding plays a key role in its organization, including the gradient-like patterning of gene expression that encodes the molecular underpinning of functional specialization. However, we do not yet understand how changes in brain shape and size that occur across development influence the brain's transcriptional architecture. Here we investigate the spatial embedding of transcriptional patterns of over 1800 genes across seven time points through mouse-brain development using data from the Allen Developing Mouse Brain Atlas. We find that transcriptional similarity decreases exponentially with separation distance across all developmental time points, with a correlation length scale that follows a power-law scaling relationship with a linear dimension of brain size. This scaling suggests that the mouse brain achieves a characteristic balance between local molecular similarity (homogeneous gene expression within a specialized brain area) and longer-range diversity (between functionally specialized brain areas) throughout its development. Extrapolating this mouse developmental scaling relationship to the human cortex yields a prediction consistent with the value measured from microarray data. We introduce a simple model of brain growth as spatially autocorrelated gene-expression gradients that expand through development, which captures key features of the mouse developmental data. Complementing the well-known exponential distance rule for structural connectivity, our findings characterize an analogous exponential distance rule for transcriptional gradients that scales across mouse brain development, providing new understanding of spatial constraints on the brain's molecular patterning.

Hypoxia-Ischemia Induced Age-Dependent Gene Transcription Effects at Two Development Stages in the Neonate Mouse Brain

Human brain lesions in the perinatal period result in life-long neuro-disabilities impairing sensory-motor, cognitive, and behavior functions for years. Topographical aspects of brain lesions depend on gestational age at the time of insult in preterm or term infants and impaired subsequent steps of brain development and maturation. In mice, the Rice-Vannucci procedure of neonate hypoxia-ischemia (HI) was used at 5 days (P5) or P10, mimicking the development of 30 week-gestation fetus/preterm newborn, or full-term infant, respectively. Transcription response to HI was assessed at 3, 6, 12, and 24 h after insult, using micro-array technology. Statistical Pathway and Gene Ontology terms enrichments were investigated using DAVID®, Revigo® and Ingenuity Pathway Analysis (IPA®) to identify a core of transcription response to HI, age-specific regulations, and interactions with spontaneous development. Investigations were based on direction, amplitude, and duration of responses, basal expression, and annotation. Five major points deserve attention; (i) inductions exceeded repressions (60/40%) at both ages, (ii) only 20.3% (393/1938 records) were common to P5 and P10 mice, (iii) at P5, HI effects occurred early and decreased 24 h after insult whereas they were delayed at P10 and increased 24 h after insult, (iv) common responses at P5 and P10 involved inflammation, immunity, apoptosis, and angiogenesis. (v) age-specific effects occurred with higher statistical significance at P5 than at P10. Transient repression of 12 genes encoding cholesterol biosynthesis enzymes was transiently observed 12 h after HI at P5. Synaptogenesis appeared inhibited at P5 while induced at P10, showing reciprocal effects on glutamate receptors. Specific involvement of Il-1 (interleukin-1) implicated in the firing of inflammation was observed at P10. This study pointed out age-differences in HI responses kinetics, e.g., a long-lasting inflammatory response at P10 compared to P5. Whether the specific strong depression of cholesterol biosynthesis genes that could account for white matter-specific vulnerability at P5 or prevent delayed inflammation needs further investigation. Determination of putative involvement of Il-1 and the identification of upstream regulators involved in the delayed inflammation firing at P10 appears promising routes of research in the understandings of age-dependent vulnerabilities in the neonatal brain.

Transcriptome Analysis Identifies SenZfp536, a Sense LncRNA that Suppresses Self-renewal of Cortical Neural Progenitors

Long non-coding RNAs (lncRNAs) regulate transcription to control development and homeostasis in a variety of tissues and organs. However, their roles in the development of the cerebral cortex have not been well elucidated. Here, a bioinformatics pipeline was applied to delineate the dynamic expression and potential cis-regulating effects of mouse lncRNAs using transcriptome data from 8 embryonic time points and sub-regions of the developing cerebral cortex. We further characterized a sense lncRNA, SenZfp536, which is transcribed downstream of and partially overlaps with the protein-coding gene Zfp536. Both SenZfp536 and Zfp536 were predominantly expressed in the proliferative zone of the developing cortex. Zfp536 was cis-regulated by SenZfp536, which facilitates looping between the promoter of Zfp536 and the genomic region that transcribes SenZfp536. Surprisingly, knocking down or activating the expression of SenZfp536 increased or compromised the proliferation of cortical neural progenitor cells (NPCs), respectively. Finally, overexpressing Zfp536 in cortical NPCs reversed the enhanced proliferation of cortical NPCs caused by SenZfp536 knockdown. The study deepens our understanding of how lncRNAs regulate the propagation of cortical NPCs through cis-regulatory mechanisms.

Statistical inference of differential RNA-editing sites from RNA-sequencing data by hierarchical modeling

MOTIVATION

RNA-sequencing (RNA-seq) enables global identification of RNA-editing sites in biological systems and disease. A salient step in many studies is to identify editing sites that statistically associate with treatment (e.g. case versus control) or covary with biological factors, such as age. However, RNA-seq has technical features that incumbent tests (e.g. t-test and linear regression) do not consider, which can lead to false positives and false negatives.

RESULTS

In this study, we demonstrate the limitations of currently used tests and introduce the method, RNA-editing tests (REDITs), a suite of tests that employ beta-binomial models to identify differential RNA editing. The tests in REDITs have higher sensitivity than other tests, while also maintaining the type I error (false positive) rate at the nominal level. Applied to the GTEx dataset, we unveil RNA-editing changes associated with age and gender, and differential recoding profiles between brain regions.

AVAILABILITY AND IMPLEMENTATION

REDITs are implemented as functions in R and freely available for download at https://github.com/gxiaolab/REDITs. The repository also provides a code example for leveraging parallelization using multiple cores.

Specific histone modifications associate with alternative exon selection during mammalian development

Alternative splicing (AS) is frequent during early mouse embryonic development. Specific histone post-translational modifications (hPTMs) have been shown to regulate exon splicing by either directly recruiting splice machinery or indirectly modulating transcriptional elongation. In this study, we hypothesized that hPTMs regulate expression of alternatively spliced genes for specific processes during differentiation. To address this notion, we applied an innovative machine learning approach to relate global hPTM enrichment to AS regulation during mammalian tissue development. We found that specific hPTMs, H3K36me3 and H3K4me1, play a role in skipped exon selection among all the tissues and developmental time points examined. In addition, we used iterative random forest model and found that interactions of multiple hPTMs most strongly predicted splicing when they included H3K36me3 and H3K4me1. Collectively, our data demonstrated a link between hPTMs and alternative splicing which will drive further experimental studies on the functional relevance of these modifications to alternative splicing.

Shiny-DEG: A Web Application to Analyze and Visualize Differentially Expressed Genes in RNA-seq

RNA-seq analysis has become one of the most widely used methods for biological and medical experiments, aiming to identify differentially expressed genes at a large scale. However, due to lack of programming skills and statistical background, it is difficult for biologists including faculty and students to fully understand what the RNA-seq results are and how to interpret them. In recent years, even though, there are several programs or websites that assist researchers to analyze and visualize NGS results, they have several limitations. Therefore, Shiny-DEG, a web application that facilitates the exploration and visualization of differentially expressed genes from RNA-seq, was developed. It integrates multi-factor design experiments, allows users to modify the parameters interactively according to experiments purpose and all analysis results can be downloaded directly, aiming to further assisting the interpretation and explanation of the biological questions. Therefore, it serves better for biologists without programming skills. Overall, this project is of great significance to reveal the mechanism of transcriptome differences.

Sierra: discovery of differential transcript usage from polyA-captured single-cell RNA-seq data

High-throughput single-cell RNA-seq (scRNA-seq) is a powerful tool for studying gene expression in single cells. Most current scRNA-seq bioinformatics tools focus on analysing overall expression levels, largely ignoring alternative mRNA isoform expression. We present a computational pipeline, Sierra, that readily detects differential transcript usage from data generated by commonly used polyA-captured scRNA-seq technology. We validate Sierra by comparing cardiac scRNA-seq cell types to bulk RNA-seq of matched populations, finding significant overlap in differential transcripts. Sierra detects differential transcript usage across human peripheral blood mononuclear cells and the Tabula Muris, and 3 'UTR shortening in cardiac fibroblasts. Sierra is available at https://github.com/VCCRI/Sierra .

Pig-specific RNA editing during early embryo development revealed by genome-wide comparisons

Posttranscriptional modification of mRNA sequences through RNA editing can increase transcriptome and proteome diversity in eukaryotes. Studies of fetal and adult tissues showed that adenosine‐to‐inosine RNA editing plays a crucial role in early human development, but there is a lack of global understanding of dynamic RNA editing during mammalian early embryonic development. Therefore, here we used RNA sequencing data from human, pig and mouse during early embryonic development to detect edited genes that may regulate stem cell pluripotency. We observed that although most of the RNA editing sites are located in intergenic, intron and UTR, a few editing sites are in coding regions and may result in nonsynonymous amino acid changes. Some editing sites are predicted to change the structure of a protein. We also report that HNF1A, TBX3, ACLY, ECI1 and ERDR1 are related to embryonic development and cell division.

Role of pericyte-derived SENP1 in neuronal injury after brain ischemia

Aims

SUMOylation is a posttranslational modification related to multiple human diseases. SUMOylation can be reversed by classes of proteases known as the sentrin/SUMO‐specific proteases (SENPs). In the present study, we investigate the potential role of SENP1 in pericytes in the brain ischemia.

Methods

Pericyte‐specific deletion of senp1 mice (Cspg4‐Cre; senp1^f/f) were used for brain function and neuronal damage evaluation following brain ischemia. The cerebral blood vessels of diameter, velocity, and flux were performed in living mice by two‐photon laser scanning microscopy (TPLSM). Biochemical analysis and immunohistochemistry methods were used to address the role and mechanism of pericyte‐specific SENP1 in the pathological process of brain ischemia. A coculture model of HBVPs and HBMECs mimicked the BBB in vitro and was used to evaluate BBB integrity after glucose deprivation.

Results

Our results showed that senp1‐specific deletion in pericytes did not affect the motor function and cognitive function of mice. However, the pericyte‐specific deletion of senp1 aggravated the infarct size and motor deficit following focal brain ischemia. Consistently, the TPLSM data demonstrated that SENP1 deletion in pericytes accelerated thrombosis formation in brain microvessels. We also found that pericyte‐specific deletion of senp1 exaggerated the neuronal damage significantly following brain ischemia in mice. Moreover, SENP1 knockdown in pericytes could activate the apoptosis signaling and disrupt the barrier integrity in vitro coculture model.

Conclusions

Our findings revealed that targeting SENP1 in pericytes may represent a novel therapeutic strategy for neurovascular protection in stroke.

Asian Zika Virus Isolate Significantly Changes the Transcriptional Profile and Alternative RNA Splicing Events in a Neuroblastoma Cell Line

The alternative splicing of pre-mRNAs expands a single genetic blueprint to encode multiple, functionally diverse protein isoforms. Viruses have previously been shown to interact with, depend on, and alter host splicing machinery. The consequences, however, incited by viral infection on the global alternative slicing (AS) landscape are under-appreciated. Here, we investigated the transcriptional and alternative splicing profile of neuronal cells infected with a contemporary Puerto Rican Zika virus (ZIKV^PR) isolate, an isolate of the prototypical Ugandan ZIKV (ZIKV^MR), and dengue virus 2 (DENV2). Our analyses revealed that ZIKV^PR induced significantly more differential changes in expressed genes compared to ZIKV^MR or DENV2, despite all three viruses showing equivalent infectivity and viral RNA levels. Consistent with the transcriptional profile, ZIKV^PR induced a higher number of alternative splicing events compared to ZIKV^MR or DENV2, and gene ontology analyses highlighted alternative splicing changes in genes associated with mRNA splicing. In summary, we show that ZIKV affects cellular RNA homeostasis not only at the transcriptional levels but also through the alternative splicing of cellular transcripts. These findings could provide new molecular insights into the neuropathologies associated with this virus.

Regulating Divergent Transcriptomes through mRNA Splicing and Its Modulation Using Various Small Compounds

Human transcriptomes are more divergent than genes and contribute to the sophistication of life. This divergence is derived from various isoforms arising from alternative splicing. In addition, alternative splicing regulated by spliceosomal factors and RNA structures, such as the RNA G-quadruplex, is important not only for isoform diversity but also for regulating gene expression. Therefore, abnormal splicing leads to serious diseases such as cancer and neurodegenerative disorders. In the first part of this review, we describe the regulation of divergent transcriptomes using alternative mRNA splicing. In the second part, we present the relationship between the disruption of splicing and diseases. Recently, various compounds with splicing inhibitor activity were established. These splicing inhibitors are recognized as a biological tool to investigate the molecular mechanism of splicing and as a potential therapeutic agent for cancer treatment. Food-derived compounds with similar functions were found and are expected to exhibit anticancer effects. In the final part, we describe the compounds that modulate the messenger RNA (mRNA) splicing process and their availability for basic research and future clinical potential.

Uncovering the mouse olfactory long non-coding transcriptome with a novel machine-learning model

Very little is known about long non-coding RNAs (lncRNAs) in the mammalian olfactory sensory epithelia. Deciphering the non-coding transcriptome in olfaction is relevant because these RNAs have been shown to play a role in chromatin modification and nuclear architecture reorganization, processes that accompany olfactory differentiation and olfactory receptor gene choice, one of the most poorly understood gene regulatory processes in mammals. In this study, we used a combination of in silico and ex vivo approaches to uncover a comprehensive catalogue of olfactory lncRNAs and to investigate their expression in the mouse olfactory organs. Initially, we used a novel machine-learning lncRNA classifier to discover hundreds of annotated and unannotated lncRNAs, some of which were predicted to be preferentially expressed in the main olfactory epithelium and the vomeronasal organ, the most important olfactory structures in the mouse. Moreover, we used whole-tissue and single-cell RNA sequencing data to discover lncRNAs expressed in mature sensory neurons of the main epithelium. Candidate lncRNAs were further validated by in situ hybridization and RT-PCR, leading to the identification of lncRNAs found throughout the olfactory epithelia, as well as others exquisitely expressed in subsets of mature olfactory neurons or progenitor cells.

Polymodal Nociception in Drosophila Requires Alternative Splicing of TrpA1

Transcripts of noxious stimulus-detecting TrpA1 channels are alternatively spliced. Despite the importance of nociception for survival, the in vivo significance of expressing different TrpA1 isoforms is largely unknown. Here, we develop a novel genetic approach to generate Drosophila knockin strains expressing single TrpA1 isoforms. Drosophila TrpA1 mediates heat and UVC-triggered nociception. We show that TrpA1-C and TrpA1-D, two alternative isoforms, are co-expressed in nociceptors. When examined in heterologous cells, both TrpA1-C and TrpA1-D are activated by heat and UVC. By contrast, analysis of knockin flies reveals the striking functional specificity; TrpA1-C mediates UVC-nociception, whereas TrpA1-D mediates heat-nociception. Therefore, in vivo functions of TrpA1-C and TrpA1-D are different from each other and are different from their in vitro properties. Our results indicate that a given sensory stimulus preferentially activates a single TrpA1 isoform in vivo and that polymodal nociception requires co-expression of TrpA1 isoforms, providing novel insights of how alternative splicing regulates nociception.

Alternative RNA Splicing Associated With Mammalian Neuronal Differentiation

Alternative pre-mRNA splicing (AS) produces multiple isoforms of mRNAs and proteins from a single gene. It is most prevalent in the mammalian brain and is thought to contribute to the formation and/or maintenance of functional complexity of the brain. Increasing evidence has documented the significant changes of AS between different regions or different developmental stages of the brain, however, the dynamics of AS and the possible function of it during neural progenitor cell (NPC) differentiation is less well known. Here, using purified NPCs and their progeny neurons isolated from the embryonic mouse cerebral cortex, we characterized the global differences of AS events between the 2 cell types by deep sequencing. The sequencing results revealed cell type-specific AS in NPCs and neurons that are important for distinct functions pertinent to the corresponding cell type. Our data may serve as a resource useful for further understanding how AS contributes to molecular regulations in NPCs and neurons during cortical development.

Translating neural stem cells to neurons in the mammalian brain

The mammalian neocortex underlies our perception of sensory information, performance of motor activities, and higher-order cognition. During mammalian embryogenesis, radial glial precursor cells sequentially give rise to diverse populations of excitatory cortical neurons, followed by astrocytes and oligodendrocytes. A subpopulation of these embryonic neural precursors persists into adulthood as neural stem cells, which give rise to inhibitory interneurons and glia. Although the intrinsic mechanisms instructing the genesis of these distinct progeny have been well-studied, most work to date has focused on transcriptional, epigenetic, and cell-cycle control. Recent studies, however, have shown that posttranscriptional mechanisms also regulate the cell fate choices of transcriptionally primed neural precursors during cortical development. These mechanisms are mediated primarily by RNA-binding proteins and microRNAs that coordinately regulate mRNA translation, stability, splicing, and localization. Together, these findings point to an extensive network of posttranscriptional control and provide insight into both normal cortical development and disease. They also add another layer of complexity to brain development and raise important biological questions for future investigation.

Exposure to Ionizing Radiation Triggers Prolonged Changes in Circular RNA Abundance in the Embryonic Mouse Brain and Primary Neurons

The exposure of mouse embryos in utero and primary cortical neurons to ionizing radiation results in the P53-dependent activation of a subset of genes that is highly induced during brain development and neuronal maturation, a feature that these genes reportedly share with circular RNAs (circRNAs). Interestingly, some of these genes are predicted to express circular transcripts. In this study, we validated the abundance of the circular transcript variants of four P53 target genes (Pvt1, Ano3, Sec14l5, and Rnf169). These circular variants were overall more stable than their linear counterparts. They were furthermore highly enriched in the brain and their transcript levels continuously increase during subsequent developmental stages (from embryonic day 12 until adulthood), while no further increase could be observed for linear mRNAs beyond post-natal day 30. Finally, whereas radiation-induced expression of P53 target mRNAs peaks early after exposure, several of the circRNAs showed prolonged induction in irradiated embryonic mouse brain, primary mouse cortical neurons, and mouse blood. Together, our results indicate that the circRNAs from these P53 target genes are induced in response to radiation and they corroborate the findings that circRNAs may represent biomarkers of brain age. We also propose that they may be superior to mRNA as long-term biomarkers for radiation exposure.

Long non-coding RNA LncKdm2b regulates cortical neuronal differentiation by cis-activating Kdm2b

The mechanisms underlying spatial and temporal control of cortical neurogenesis of the brain are largely elusive. Long non-coding RNAs (lncRNAs) have emerged as essential cell fate regulators. Here we found LncKdm2b (also known as Kancr), a lncRNA divergently transcribed from a bidirectional promoter of Kdm2b, is transiently expressed during early differentiation of cortical projection neurons. Interestingly, Kdm2b’s transcription is positively regulated in cis by LncKdm2b, which has intrinsic-activating function and facilitates a permissive chromatin environment at the Kdm2b’s promoter by associating with hnRNPAB. Lineage tracing experiments and phenotypic analyses indicated LncKdm2b and Kdm2b are crucial in proper differentiation and migration of cortical projection neurons. These observations unveiled a lncRNA-dependent machinery in regulating cortical neuronal differentiation.

Gene expression across mammalian organ development

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, rhesus macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified correspondences of developmental stages across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

Alcohol exposure alters pre-mRNA splicing of antiapoptotic Mcl-1L isoform and induces apoptosis in neural progenitors and immature neurons

Alternative splicing and expression of splice variants of genes in the brain may lead to the modulation of protein functions, which may ultimately influence behaviors associated with alcohol dependence and neurotoxicity. We recently showed that ethanol exposure can lead to pre-mRNA missplicing of Mcl-1, a pro-survival member of the Bcl-2 family, by downregulating the expression levels of serine/arginine rich splicing factor 1 (SRSF1). Little is known about the physiological expression of these isoforms in neuronal cells and their role in toxicity induced by alcohol exposure during the developmental period. In order to investigate the impact of alcohol exposure on alternative splicing of Mcl-1 pre-mRNA and its role in neurotoxicity, we developed a unique primary human neuronal culture model where neurospheres (hNSPs), neural progenitors (hNPCs), immature neurons, and mature neurons were cultured from the matching donor fetal brain tissues. Our data suggest that neural progenitors and immature neurons are highly sensitive to the toxic effects of ethanol, while mature neuron cultures showed resistance to ethanol exposure. Further analysis of Mcl-1 pre-mRNA alternative splicing by semi-quantitative and quantitative analysis revealed that ethanol exposure causes a significant decrease in Mcl-1L/Mcl-1S ratio in a dose and time dependent manner in neural progenitors. Interestingly, ectopic expression of Mcl-1L isoform in neural progenitors was able to recover the viability loss and apoptosis induced by alcohol exposure. Altogether, these observations suggest that alternative splicing of Mcl-1 may play a crucial role in neurotoxicity associated with alcohol exposure in the developing fetal brain.