HITS-CLIP and integrative modeling define the Rbfox splicing-regulatory network linked to brain development and autism

Post-transcriptional mechanisms controlling neurogenesis and direct neuronal reprogramming

Neurogenesis is a tightly regulated process in time and space both in the developing embryo and in adult neurogenic niches. A drastic change in the transcriptome and proteome of radial glial cells or neural stem cells towards the neuronal state is achieved due to sophisticated mechanisms of epigenetic, transcriptional, and post-transcriptional regulation. Understanding these neurogenic mechanisms is of major importance, not only for shedding light on very complex and crucial developmental processes, but also for the identification of putative reprogramming factors, that harbor hierarchically central regulatory roles in the course of neurogenesis and bare thus the capacity to drive direct reprogramming towards the neuronal fate. The major transcriptional programs that orchestrate the neurogenic process have been the focus of research for many years and key neurogenic transcription factors, as well as repressor complexes, have been identified and employed in direct reprogramming protocols to convert non-neuronal cells, into functional neurons. The post-transcriptional regulation of gene expression during nervous system development has emerged as another important and intricate regulatory layer, strongly contributing to the complexity of the mechanisms controlling neurogenesis and neuronal function. In particular, recent advances are highlighting the importance of specific RNA binding proteins that control major steps of mRNA life cycle during neurogenesis, such as alternative splicing, polyadenylation, stability, and translation. Apart from the RNA binding proteins, microRNAs, a class of small non-coding RNAs that block the translation of their target mRNAs, have also been shown to play crucial roles in all the stages of the neurogenic process, from neural stem/progenitor cell proliferation, neuronal differentiation and migration, to functional maturation. Here, we provide an overview of the most prominent post-transcriptional mechanisms mediated by RNA binding proteins and microRNAs during the neurogenic process, giving particular emphasis on the interplay of specific RNA binding proteins with neurogenic microRNAs. Taking under consideration that the molecular mechanisms of neurogenesis exert high similarity to the ones driving direct neuronal reprogramming, we also discuss the current advances in in vitro and in vivo direct neuronal reprogramming approaches that have employed microRNAs or RNA binding proteins as reprogramming factors, highlighting the so far known mechanisms of their reprogramming action.

The Rbfox1/LASR complex controls alternative pre-mRNA splicing by recognition of multi-part RNA regulatory modules

The Rbfox proteins regulate alternative pre-mRNA splicing by binding to the RNA element GCAUG. In the nucleus, most of Rbfox is bound to LASR, a complex of RNA-binding proteins that recognize additional RNA motifs. However, it remains unclear how the different subunits of the Rbfox/LASR complex act together to bind RNA and regulate splicing. We used a nuclease-protection assay to map the transcriptome-wide footprints of Rbfox1/LASR on nascent cellular RNA. In addition to GCAUG, Rbfox1/LASR binds RNA containing motifs for LASR subunits hnRNPs M, H/F, C, and Matrin3. These elements are often arranged in tandem, forming multi-part modules of RNA motifs. To distinguish contact sites of Rbfox1 from the LASR subunits, we analyzed a mutant Rbfox1(F125A) that has lost RNA binding but remains associated with LASR. Rbfox1(F125A)/LASR complexes no longer interact with GCAUG but retain binding to RNA elements for LASR. Splicing analyses reveal that in addition to activating exons through adjacent GCAUG elements, Rbfox can also stimulate exons near binding sites for LASR subunits. Mini-gene experiments demonstrate that these diverse elements produce a combined regulatory effect on a target exon. These findings illuminate how a complex of RNA-binding proteins can decode combinatorial splicing regulatory signals by recognizing groups of tandem RNA elements.

Accelerated evolution in the human lineage led to gain and loss of transcriptional enhancers in the RBFOX1 locus

A long-standing goal of evolutionary biology is to decode how changes in gene regulatory networks contribute to human-specific traits. Human accelerated regions (HARs) are prime candidates for driving gene regulatory modifications in human development. The RBFOX1 locus is densely populated with HARs, providing a set of potential regulatory elements that could have changed its expression in the human lineage. Here, we examined the role of RBFOX1-HARs using transgenic zebrafish reporter assays and identified 15 transcriptional enhancers that are active in the developing nervous system, 9 of which displayed differential activity between the human and chimpanzee sequences. The engineered loss of two selected RBFOX1-HARs in knockout mouse models modified Rbfox1 expression at specific developmental stages and tissues in the brain, influencing the expression and splicing of a high number of Rbfox1 target genes. Our results provided insight into the spatial and temporal changes in gene expression driven by RBFOX1-HARs.

Mapping RNA-protein interactions with subcellular resolution using colocalization CLIP

RNA-binding proteins (RBPs) are essential for RNA metabolism and profoundly impact health and disease. The subcellular organization of RBP interaction networks with target RNAs remains largely unexplored. Here, we develop colocalization CLIP (coCLIP), a method that combines cross-linking and immunoprecipitation (CLIP) with proximity labeling, to explore in-depth the subcellular RNA interactions of the RBP human antigen R (HuR). Using this method, we uncover HuR's dynamic and location-specific interactions with RNA, revealing alterations in sequence preferences and interactions in the nucleus, cytosol, or stress granule (SG) compartments. We uncover HuR's unique binding preferences within SGs during arsenite stress, illuminating intricate interactions that conventional methodologies cannot capture. Overall, coCLIP provides a powerful method for revealing RBP-RNA interactions based on localization and lays the foundation for an advanced understanding of RBP models that incorporate subcellular location as a critical determinant of their functions.

Reverse engineering neuron type-specific and type-orthogonal splicing-regulatory networks using single-cell transcriptomes

Cell type-specific alternative splicing (AS) enables differential gene isoform expression between diverse neuron types with distinct identities and functions. Current studies linking individual RNA-binding proteins (RBPs) to AS in a few neuron types underscore the need for holistic modeling. Here, we use network reverse engineering to derive a map of the neuron type-specific AS regulatory landscape from 133 mouse neocortical cell types defined by single-cell transcriptomes. This approach reliably inferred the regulons of 350 RBPs and their cell type-specific activities. Our analysis revealed driving factors delineating neuronal identities, among which we validated Elavl2 as a key RBP for MGE-specific splicing in GABAergic interneurons using an in vitro ESC differentiation system. We also identified a module of exons and candidate regulators specific for long- and short-projection neurons across multiple neuronal classes. This study provides a resource for elucidating splicing regulatory programs that drive neuronal molecular diversity, including those that do not align with gene expression-based classifications.

m6a methylation orchestrates IMP1 regulation of microtubules during human neuronal differentiation

Neuronal differentiation requires building a complex intracellular architecture, and therefore the coordinated regulation of defined sets of genes. RNA-binding proteins (RBPs) play a key role in this regulation. However, while their action on individual mRNAs has been explored in depth, the mechanisms used to coordinate gene expression programs shaping neuronal morphology are poorly understood. To address this, we studied how the paradigmatic RBP IMP1 (IGF2BP1), an essential developmental factor, selects and regulates its RNA targets during the human neuronal differentiation. We perform a combination of system-wide and molecular analyses, revealing that IMP1 developmentally transitions to and directly regulates the expression of mRNAs encoding essential regulators of the microtubule network, a key component of neuronal morphology. Furthermore, we show that m6A methylation drives the selection of specific IMP1 mRNA targets and their protein expression during the developmental transition from neural precursors to neurons, providing a molecular principle for the onset of target selectivity.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

Structure-based prediction and characterization of photo-crosslinking in native protein-RNA complexes

UV-crosslinking of protein and RNA in direct contacts has been widely used to study protein-RNA complexes while our understanding of the photo-crosslinking mechanisms remains poor. This knowledge gap is due to the challenge of precisely mapping the crosslink sites in protein and RNA simultaneously in their native sequence and structural contexts. Here we systematically analyze protein-RNA interactions and photo-crosslinking by bridging crosslinked nucleotides and amino acids mapped using different assays with protein-RNA complex structures. We developed a computational method PxR3D-map which reliably predicts crosslink sites using structural information characterizing protein-RNA interaction interfaces. Analysis of the informative features revealed that photo-crosslinking is facilitated by base stacking with not only aromatic residues, but also dipeptide bonds that involve glycine, and distinct mechanisms are utilized by different RNA-binding domains. Our work suggests protein-RNA photo-crosslinking is highly selective in the cellular environment, which can guide data interpretation and further technology development for UV-crosslinking-based assays.

The spectrum of pre-mRNA splicing in autism

Disruptions in spatiotemporal gene expression can result in atypical brain function. Specifically, autism spectrum disorder (ASD) is characterized by abnormalities in pre-mRNA splicing. Abnormal splicing patterns have been identified in the brains of individuals with ASD, and mutations in splicing factors have been found to contribute to neurodevelopmental delays associated with ASD. Here we review studies that shed light on the importance of splicing observed in ASD and that explored the intricate relationship between splicing factors and ASD, revealing how disruptions in pre-mRNA splicing may underlie ASD pathogenesis. We provide an overview of the research regarding all splicing factors associated with ASD and place a special emphasis on five specific splicing factors-HNRNPH2, NOVA2, WBP4, SRRM2, and RBFOX1-known to impact the splicing of ASD-related genes. In the discussion of the molecular mechanisms influenced by these splicing factors, we lay the groundwork for a deeper understanding of ASD's complex etiology. Finally, we discuss the potential benefit of unraveling the connection between splicing and ASD for the development of more precise diagnostic tools and targeted therapeutic interventions. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Evolution and Genomics > Computational Analyses of RNA RNA-Based Catalysis > RNA Catalysis in Splicing and Translation.

TCF4 and RBFOX1 as peripheral biomarkers for the differential diagnosis and treatment of major depressive disorder

BACKGROUND

Recent genome-wide association studies on major depressive disorder (MDD) have indicated the involvement of LRFN5 and OLFM4; however, the expression levels and roles of these molecules in MDD remain unclear. The present study aimed to determine the serum levels of TCF4 and RBFOX1 in patients with MDD and to investigate whether these molecules could be used as biomarkers for MDD diagnosis.

METHODS

The study included 99 drug-naïve MDD patients, 90 drug-treated MDD patients, and 81 healthy controls (HCs). Serum TCF4 and RBFOX1 levels were measured by ELISA. Pearson's correlation analysis was conducted to determine the association between TCF4/RBFOX1 and clinical variables. Linear support vector machine classifier was used to evaluate the diagnostic capabilities of TCF4 and RBFOX1.

RESULTS

Serum TCF4 and RBFOX1 levels were substantially higher in MDD patients than in HCs and significantly lower in drug-treated MDD patients than in drug-naïve MDD patients. Moreover, serum TCF4 and RBFOX1 levels were associated with the Hamilton Depression Scale score, duration of illness, serum lipids levels, and hepatic function. Thus, both these molecules showed potential as biomarkers for MDD. TCF4 and RBFOX1 combination exhibited a higher diagnostic performance, with the mean area under the curve values of 0.9861 and 0.9936 in the training and testing sets, respectively.

LIMITATIONS

Small sample size and investigation of only the peripheral nervous system.

CONCLUSIONS

TCF4 and RBFOX1 may be involved in the pathogenesis of MDD, and their combination may serve as a diagnostic biomarker panel for MDD.

Single-cell long-read sequencing in human cerebral organoids uncovers cell-type-specific and autism-associated exons

Dysregulation of alternative splicing has been repeatedly associated with neurodevelopmental disorders, but the extent of cell-type-specific splicing in human neural development remains largely uncharted. Here, single-cell long-read sequencing in induced pluripotent stem cell (iPSC)-derived cerebral organoids identifies over 31,000 uncatalogued isoforms and 4,531 cell-type-specific splicing events. Long reads uncover coordinated splicing and cell-type-specific intron retention events, which are challenging to study with short reads. Retained neuronal introns are enriched in RNA splicing regulators, showing shorter lengths, higher GC contents, and weaker 5' splice sites. We use this dataset to explore the biological processes underlying neurological disorders, focusing on autism. In comparison with prior transcriptomic data, we find that the splicing program in autistic brains is closer to the progenitor state than differentiated neurons. Furthermore, cell-type-specific exons harbor significantly more de novo mutations in autism probands than in siblings. Overall, these results highlight the importance of cell-type-specific splicing in autism and neuronal gene regulation.

Modulation of insulin secretion by RBFOX2-mediated alternative splicing

Insulin secretion is a tightly regulated process that is vital for maintaining blood glucose homeostasis. Although the molecular components of insulin granule trafficking and secretion are well established, how they are regulated to rapidly fine-tune secretion in response to changing environmental conditions is not well characterized. Recent studies have determined that dysregulation of RNA-binding proteins (RBPs) and aberrant mRNA splicing occurs at the onset of diabetes. We demonstrate that the RBP, RBFOX2, is a critical regulator of insulin secretion through the alternative splicing of genes required for insulin granule docking and exocytosis. Conditional mutation of Rbfox2 in the mouse pancreas results in decreased insulin secretion and impaired blood glucose homeostasis. Consistent with defects in secretion, we observe reduced insulin granule docking and corresponding splicing defects in the SNARE complex components. These findings identify an additional mechanism for modulating insulin secretion in both healthy and dysfunctional pancreatic β cells.

Mapping RNA-Protein Interactions with Subcellular Resolution Using Colocalization CLIP

RNA binding proteins (RBPs) are essential for RNA metabolism and profoundly impact health and disease. The subcellular organization of RBP interaction networks with target RNAs remains largely unexplored. Here, we develop colocalization CLIP, a method that combines CrossLinking and ImmunoPrecipitation (CLIP) with proximity labeling, to explore in-depth the subcellular RNA interactions of the well-studied RNA-binding protein HuR. Using this method, we uncover HuR's dynamic and location-specific interactions with RNA, revealing alterations in sequence preferences and interactions in the nucleus, cytosol, or stress granule compartments. We uncover HuR's unique binding preferences within stress granules during arsenite stress, illuminating intricate interactions that conventional methodologies cannot capture. Overall, coCLIP provides a powerful method for revealing RBP:RNA interactions based on localization and lays the foundation for an advanced understanding of RBP models that incorporate subcellular location as a critical determinant of their functions.

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

Schizophrenia-associated NRXN1 deletions induce developmental-timing- and cell-type-specific vulnerabilities in human brain organoids

De novo mutations and copy number deletions in NRXN1 (2p16.3) pose a significant risk for schizophrenia (SCZ). It is unclear how NRXN1 deletions impact cortical development in a cell type-specific manner and disease background modulates these phenotypes. Here, we leveraged human pluripotent stem cell-derived forebrain organoid models carrying NRXN1 heterozygous deletions in isogenic and SCZ patient genetic backgrounds and conducted single-cell transcriptomic analysis over the course of brain organoid development from 3 weeks to 3.5 months. Intriguingly, while both deletions similarly impacted molecular pathways associated with ubiquitin-proteasome system, alternative splicing, and synaptic signaling in maturing glutamatergic and GABAergic neurons, SCZ-NRXN1 deletions specifically perturbed developmental trajectories of early neural progenitors and accumulated disease-specific transcriptomic signatures. Using calcium imaging, we found that both deletions led to long-lasting changes in spontaneous and synchronous neuronal networks, implicating synaptic dysfunction. Our study reveals developmental-timing- and cell-type-dependent actions of NRXN1 deletions in unique genetic contexts.

PIE-seq: identifying RNA-binding protein targets by dual RNA-deaminase editing and sequencing

RNA-binding proteins (RBPs) are essential for gene regulation, but it remains a challenge to identify their RNA targets across cell types. Here we present PIE-Seq to investigate Protein-RNA Interaction with dual-deaminase Editing and Sequencing by conjugating C-to-U and A-to-I base editors to RBPs. We benchmark PIE-Seq and demonstrate its sensitivity in single cells, its application in the developing brain, and its scalability with 25 human RBPs. Bulk PIE-Seq identifies canonical binding features for RBPs such as PUM2 and NOVA1, and nominates additional target genes for most tested RBPs such as SRSF1 and TDP-43/TARDBP. Homologous RBPs frequently edit similar sequences and gene sets in PIE-Seq while different RBP families show distinct targets. Single-cell PIE-PUM2 uncovers comparable targets to bulk samples and applying PIE-PUM2 to the developing mouse neocortex identifies neural-progenitor- and neuron-specific target genes such as App. In summary, PIE-Seq provides an orthogonal approach and resource to uncover RBP targets in mice and human cells.

Dosage sensitivity to Pumilio1 variants in the mouse brain reflects distinct molecular mechanisms

Different mutations in the RNA-binding protein Pumilio1 (PUM1) cause divergent phenotypes whose severity tracks with dosage: a mutation that reduces PUM1 levels by 25% causes late-onset ataxia, whereas haploinsufficiency causes developmental delay and seizures. Yet PUM1 targets are derepressed to equal degrees in both cases, and the more severe mutation does not hinder PUM1's RNA-binding ability. We therefore considered the possibility that the severe mutation might disrupt PUM1 interactions, and identified PUM1 interactors in the murine brain. We find that mild PUM1 loss derepresses PUM1-specific targets, but the severe mutation disrupts interactions with several RNA-binding proteins and the regulation of their targets. In patient-derived cell lines, restoring PUM1 levels restores these interactors and their targets to normal levels. Our results demonstrate that dosage sensitivity does not always signify a linear relationship with protein abundance but can involve distinct mechanisms. We propose that to understand the functions of RNA-binding proteins in a physiological context will require studying their interactions as well as their targets.

Upregulation of SYNGAP1 expression in mice and human neurons by redirecting alternative splicing

The Ras GTPase-activating protein SYNGAP1 plays a central role in synaptic plasticity, and de novo SYNGAP1 mutations are among the most frequent causes of autism and intellectual disability. How SYNGAP1 is regulated during development and how to treat SYNGAP1-associated haploinsufficiency remain challenging questions. Here, we characterize an alternative 3' splice site (A3SS) of SYNGAP1 that induces nonsense-mediated mRNA decay (A3SS-NMD) in mouse and human neural development. We demonstrate that PTBP1/2 directly bind to and promote SYNGAP1 A3SS inclusion. Genetic deletion of the Syngap1 A3SS in mice upregulates Syngap1 protein and alleviates the long-term potentiation and membrane excitability deficits caused by a Syngap1 knockout allele. We further report a splice-switching oligonucleotide (SSO) that converts SYNGAP1 unproductive isoform to the functional form in human iPSC-derived neurons. This study describes the regulation and function of SYNGAP1 A3SS-NMD, the genetic rescue of heterozygous Syngap1 knockout mice, and the development of an SSO to potentially alleviate SYNGAP1-associated haploinsufficiency.

The mRNA-binding protein DDX3 mediates TGF-β1 upregulation of translation and promotes pulmonary fibrosis

Pulmonary fibrosis is potentiated by a positive feedback loop involving the extracellular sialidase enzyme neuraminidase 3 (NEU3) causing release of active TGF-β1 and TGF-β1 upregulating NEU3 by increasing translation without affecting mRNA levels. In this report, we elucidate the TGF-β1 upregulation of the translation mechanism. In human lung fibroblasts, TGF-β1 increased levels of proteins, including NEU3, by increasing translation of the encoding mRNAs without significantly affecting levels of these mRNAs. A total of 180 of these mRNAs shared a common 20-nucleotide motif. Deletion of this motif from NEU3 mRNA eliminated the TGF-β1 upregulation of NEU3 translation, while insertion of this motif in 2 mRNAs insensitive to TGF-β1 caused TGF-β1 to upregulate their translation. RNA-binding proteins including DEAD box helicase 3, X-linked (DDX3), bind the RNA motif, and TGF-β1 regulates their protein levels and/or binding to the motif. We found that DDX3 was upregulated in the fibrotic lesions in patients with pulmonary fibrosis, and inhibiting DDX3 in fibroblasts reduced TGF-β1 upregulation of NEU3 levels. In the mouse bleomycin model of pulmonary fibrosis, injections of the DDX3 inhibitor RK-33 potentiated survival and reduced lung inflammation, fibrosis, and tissue levels of DDX3, TGF-β1, and NEU3. These results suggest that inhibiting an mRNA-binding protein that mediates TGF-β1 upregulation of translation can reduce pulmonary fibrosis.

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.

Coordination of Alternative Splicing and Alternative Polyadenylation revealed by Targeted Long-Read Sequencing

Nervous system development is associated with extensive regulation of alternative splicing (AS) and alternative polyadenylation (APA). AS and APA have been extensively studied in isolation, but little is known about how these processes are coordinated. Here, the coordination of cassette exon (CE) splicing and APA in Drosophila was investigated using a targeted long-read sequencing approach we call Pull-a-Long-Seq (PL-Seq). This cost-effective method uses cDNA pulldown and Nanopore sequencing combined with an analysis pipeline to resolve the connectivity of alternative exons to alternative 3' ends. Using PL-Seq, we identified genes that exhibit significant differences in CE splicing depending on connectivity to short versus long 3'UTRs. Genomic long 3'UTR deletion was found to alter upstream CE splicing in short 3'UTR isoforms and ELAV loss differentially affected CE splicing depending on connectivity to alternative 3'UTRs. This work highlights the importance of considering connectivity to alternative 3'UTRs when monitoring AS events.

Approaches for Modes of Action Study of Long Non-Coding RNAs: From Single Verification to Genome-Wide Determination

Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides (nt) that are not translated into known functional proteins. This broad definition covers a large collection of transcripts with diverse genomic origins, biogenesis, and modes of action. Thus, it is very important to choose appropriate research methodologies when investigating lncRNAs with biological significance. Multiple reviews to date have summarized the mechanisms of lncRNA biogenesis, their localization, their functions in gene regulation at multiple levels, and also their potential applications. However, little has been reviewed on the leading strategies for lncRNA research. Here, we generalize a basic and systemic mind map for lncRNA research and discuss the mechanisms and the application scenarios of ‘up-to-date’ techniques as applied to molecular function studies of lncRNAs. Taking advantage of documented lncRNA research paradigms as examples, we aim to provide an overview of the developing techniques for elucidating lncRNA interactions with genomic DNA, proteins, and other RNAs. In the end, we propose the future direction and potential technological challenges of lncRNA studies, focusing on techniques and applications.

A miR-124-mediated post-transcriptional mechanism controlling the cell fate switch of astrocytes to induced neurons

The microRNA (miRNA) miR-124 has been employed supplementary to neurogenic transcription factors (TFs) and other miRNAs to enhance direct neurogenic conversion. The aim of this study was to investigate whether miR-124 is sufficient to drive direct reprogramming of astrocytes to induced neurons (iNs) on its own and elucidate its independent mechanism of reprogramming action. Our data show that miR-124 is a potent driver of the reprogramming switch of astrocytes toward an immature neuronal fate by directly targeting the RNA-binding protein Zfp36L1 implicated in ARE-mediated mRNA decay and subsequently derepressing Zfp36L1 neurogenic interactome. To this end, miR-124 contribution in iNs' production largely recapitulates endogenous neurogenesis pathways, being further enhanced upon addition of the neurogenic compound ISX9, which greatly improves iNs' differentiation and functional maturation. Importantly, miR-124 is potent in guiding direct conversion of reactive astrocytes to immature iNs in vivo following cortical trauma, while ISX9 supplementation confers a survival advantage to newly produced iNs.

Independent Associated SNPs at SORCS3 and Its Protein Interactors for Multiple Brain-Related Disorders and Traits

Sortilin-related vacuolar protein sorting 10 (VPS10) domain containing receptor 3 (SORCS3) is a neuron-specific transmembrane protein involved in the trafficking of proteins between intracellular vesicles and the plasma membrane. Genetic variation at SORCS3 is associated with multiple neuropsychiatric disorders and behavioural phenotypes. Here, we undertake a systematic search of published genome-wide association studies to identify and catalogue associations between SORCS3 and brain-related disorders and traits. We also generate a SORCS3 gene-set based on protein-protein interactions and investigate the contribution of this gene-set to the heritability of these phenotypes and its overlap with synaptic biology. Analysis of association signals at SORSC3 showed individual SNPs to be associated with multiple neuropsychiatric and neurodevelopmental brain-related disorders and traits that have an impact on the experience of feeling, emotion or mood or cognitive function, while multiple LD-independent SNPs were associated with the same phenotypes. Across these SNPs, alleles associated with the more favourable outcomes for each phenotype (e.g., decreased risk of neuropsychiatric illness) were associated with increased expression of the SORCS3 gene. The SORCS3 gene-set was enriched for heritability contributing to schizophrenia (SCZ), bipolar disorder (BPD), intelligence (IQ) and education attainment (EA). Eleven genes from the SORCS3 gene-set were associated with more than one of these phenotypes at the genome-wide level, with RBFOX1 associated with SCZ, IQ and EA. Functional annotation revealed that the SORCS3 gene-set is enriched for multiple ontologies related to the structure and function of synapses. Overall, we find many independent association signals at SORCS3 with brain-related disorders and traits, with the effect possibly mediated by reduced gene expression, resulting in a negative impact on synaptic function.

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.

Self-supervised graph representation learning integrates multiple molecular networks and decodes gene-disease relationships

Leveraging molecular networks to discover disease-relevant modules is a long-standing challenge. With the accumulation of interactomes, there is a pressing need for powerful computational approaches to handle the inevitable noise and context-specific nature of biological networks. Here, we introduce Graphene, a two-step self-supervised representation learning framework tailored to concisely integrate multiple molecular networks and adapted to gene functional analysis via downstream re-training. In practice, we first leverage GNN (graph neural network) pre-training techniques to obtain initial node embeddings followed by re-training Graphene using a graph attention architecture, achieving superior performance over competing methods for pathway gene recovery, disease gene reprioritization, and comorbidity prediction. Graphene successfully recapitulates tissue-specific gene expression across disease spectrum and demonstrates shared heritability of common mental disorders. Graphene can be updated with new interactomes or other omics features. Graphene holds promise to decipher gene function under network context and refine GWAS (genome-wide association study) hits and offers mechanistic insights via decoding diseases from genome to networks to phenotypes.

Integrating de novo and inherited variants in 42,607 autism cases identifies mutations in new moderate-risk genes

To capture the full spectrum of genetic risk for autism, we performed a two-stage analysis of rare de novo and inherited coding variants in 42,607 autism cases, including 35,130 new cases recruited online by SPARK. We identified 60 genes with exome-wide significance (P < 2.5 × 10^-6), including five new risk genes (NAV3, ITSN1, MARK2, SCAF1 and HNRNPUL2). The association of NAV3 with autism risk is primarily driven by rare inherited loss-of-function (LoF) variants, with an estimated relative risk of 4, consistent with moderate effect. Autistic individuals with LoF variants in the four moderate-risk genes (NAV3, ITSN1, SCAF1 and HNRNPUL2; n = 95) have less cognitive impairment than 129 autistic individuals with LoF variants in highly penetrant genes (CHD8, SCN2A, ADNP, FOXP1 and SHANK3) (59% vs 88%, P = 1.9 × 10^-6). Power calculations suggest that much larger numbers of autism cases are needed to identify additional moderate-risk genes.

The germ cell-specific RNA binding protein RBM46 is essential for spermatogonial differentiation in mice

Control over gene expression is exerted, in multiple stages of spermatogenesis, at the post-transcriptional level by RNA binding proteins (RBPs). We identify here an essential role in mammalian spermatogenesis and male fertility for 'RNA binding protein 46' (RBM46). A highly evolutionarily conserved gene, Rbm46 is also essential for fertility in both flies and fish. We found Rbm46 expression was restricted to the mouse germline, detectable in males in the cytoplasm of premeiotic spermatogonia and meiotic spermatocytes. To define its requirement for spermatogenesis, we generated Rbm46 knockout (KO, Rbm46-/-) mice; although male Rbm46-/- mice were viable and appeared grossly normal, they were infertile. Testes from adult Rbm46-/- mice were small, with seminiferous tubules containing only Sertoli cells and few undifferentiated spermatogonia. Using genome-wide unbiased high throughput assays RNA-seq and 'enhanced crosslinking immunoprecipitation' coupled with RNA-seq (eCLIP-seq), we discovered RBM46 could bind, via a U-rich conserved consensus sequence, to a cohort of mRNAs encoding proteins required for completion of differentiation and subsequent meiotic initiation. In summary, our studies support an essential role for RBM46 in regulating target mRNAs during spermatogonia differentiation prior to the commitment to meiosis in mice.

Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs

Measuring mRNA decay in tumours is a prohibitive challenge, limiting our ability to map the post-transcriptional programs of cancer. Here, using a statistical framework to decouple transcriptional and post-transcriptional effects in RNA-seq data, we uncover the mRNA stability changes that accompany tumour development and progression. Analysis of 7760 samples across 18 cancer types suggests that mRNA stability changes are ~30% as frequent as transcriptional events, highlighting their widespread role in shaping the tumour transcriptome. Dysregulation of programs associated with >80 RNA-binding proteins (RBPs) and microRNAs (miRNAs) drive these changes, including multi-cancer inactivation of RBFOX and miR-29 families. Phenotypic activation or inhibition of RBFOX1 highlights its role in calcium signaling dysregulation, while modulation of miR-29 shows its impact on extracellular matrix organization and stemness genes. Overall, our study underlines the integral role of mRNA stability in shaping the cancer transcriptome, and provides a resource for systematic interrogation of cancer-associated stability pathways.

The role of mRNA stability in shaping the cancer transcriptome is revealed using a statistical analysis of transcriptomic data.

De novo variants in genes regulating stress granule assembly associate with neurodevelopmental disorders

Stress granules (SGs) are cytoplasmic assemblies in response to a variety of stressors. We report a new neurodevelopmental disorder (NDD) with common features of language problems, intellectual disability, and behavioral issues caused by de novo likely gene-disruptive variants in UBAP2L, which encodes an essential regulator of SG assembly. Ubap2l haploinsufficiency in mouse led to social and cognitive impairments accompanied by disrupted neurogenesis and reduced SG formation during early brain development. On the basis of data from 40,853 individuals with NDDs, we report a nominally significant excess of de novo variants within 29 genes that are not implicated in NDDs, including 3 essential genes (G3BP1, G3BP2, and UBAP2L) in the core SG interaction network. We validated that NDD-related de novo variants in newly implicated and known NDD genes, such as CAPRIN1, disrupt the interaction of the core SG network and interfere with SG formation. Together, our findings suggest the common SG pathology in NDDs.

Neuronal splicing regulator RBFOX3 mediates seizures via regulating Vamp1 expression preferentially in NPY-expressing GABAergic neurons

Epilepsy is a common neurological disorder, which has been linked to mutations or deletions of RNA binding protein, fox-1 homolog (Caenorhabditis elegans) 3 (RBFOX3)/NeuN, a neuronal splicing regulator. However, the mechanism of seizure mediation by RBFOX3 remains unknown. Here, we show that mice with deletion of Rbfox3 in gamma-aminobutyric acid (GABA) ergic neurons exhibit spontaneous seizures and high premature mortality due to increased presynaptic release, postsynaptic potential, neuronal excitability, and synaptic transmission in hippocampal dentate gyrus granule cells (DGGCs). Attenuating early excitatory gamma-aminobutyric acid (GABA) action by administering bumetanide, an inhibitor of early GABA depolarization, rescued premature mortality. Rbfox3 deletion reduced hippocampal expression of vesicle-associated membrane protein 1 (VAMP1), a GABAergic neuron-specific presynaptic protein. Postnatal restoration of VAMP1 rescued premature mortality and neuronal excitability in DGGCs. Furthermore, Rbfox3 deletion in GABAergic neurons showed fewer neuropeptide Y (NPY)-expressing GABAergic neurons. In addition, deletion of Rbfox3 in NPY-expressing GABAergic neurons lowered intrinsic excitability and increased seizure susceptibility. Our results establish RBFOX3 as a critical regulator and possible treatment path for epilepsy.

Identifying cell state-associated alternative splicing events and their coregulation

Alternative splicing shapes the transcriptome and contributes to each cell's unique identity, but single-cell RNA sequencing (scRNA-seq) has struggled to capture the impact of alternative splicing. We previously showed that low recovery of mRNAs from single cells led to erroneous conclusions about the cell-to-cell variability of alternative splicing. Here, we present a method, Psix, to confidently identify splicing that changes across a landscape of single cells, using a probabilistic model that is robust against the data limitations of scRNA-seq. Its autocorrelation-inspired approach finds patterns of alternative splicing that correspond to patterns of cell identity, such as cell type or developmental stage, without the need for explicit cell clustering, labeling, or trajectory inference. Applying Psix to data that follow the trajectory of mouse brain development, we identify exons whose alternative splicing patterns cluster into modules of coregulation. We show that the exons in these modules are enriched for binding by distinct neuronal splicing factors and that their changes in splicing correspond to changes in expression of these splicing factors. Thus, Psix reveals cell type-dependent splicing patterns and the wiring of the splicing regulatory networks that control them. Our new method will enable scRNA-seq analysis to go beyond transcription to understand the roles of post-transcriptional regulation in determining cell identity.

RNA and neuronal function: the importance of post-transcriptional regulation

The brain represents an organ with a particularly high diversity of genes that undergo post-transcriptional gene regulation through multiple mechanisms that affect RNA metabolism and, consequently, brain function. This vast regulatory process in the brain allows for a tight spatiotemporal control over protein expression, a necessary factor due to the unique morphologies of neurons. The numerous mechanisms of post-transcriptional regulation or translational control of gene expression in the brain include alternative splicing, RNA editing, mRNA stability and transport. A large number of trans-elements such as RNA-binding proteins and micro RNAs bind to specific cis-elements on transcripts to dictate the fate of mRNAs including its stability, localization, activation and degradation. Several trans-elements are exemplary regulators of translation, employing multiple cofactors and regulatory machinery so as to influence mRNA fate. Networks of regulatory trans-elements exert control over key neuronal processes such as neurogenesis, synaptic transmission and plasticity. Perturbations in these networks may directly or indirectly cause neuropsychiatric and neurodegenerative disorders. We will be reviewing multiple mechanisms of gene regulation by trans-elements occurring specifically in neurons.

DeepND: Deep multitask learning of gene risk for comorbid neurodevelopmental disorders

Autism spectrum disorder and intellectual disability are comorbid neurodevelopmental disorders with complex genetic architectures. Despite large-scale sequencing studies, only a fraction of the risk genes was identified for both. We present a network-based gene risk prioritization algorithm, DeepND, that performs cross-disorder analysis to improve prediction by exploiting the comorbidity of autism spectrum disorder (ASD) and intellectual disability (ID) via multitask learning. Our model leverages information from human brain gene co-expression networks using graph convolutional networks, learning which spatiotemporal neurodevelopmental windows are important for disorder etiologies and improving the state-of-the-art prediction in single- and cross-disorder settings. DeepND identifies the prefrontal and motor-somatosensory cortex (PFC-MFC) brain region and periods from early- to mid-fetal and from early childhood to young adulthood as the highest neurodevelopmental risk windows for ASD and ID. We investigate ASD- and ID-associated copy-number variation (CNV) regions and report our findings for several susceptibility gene candidates. DeepND can be generalized to analyze any combinations of comorbid disorders.

•

DeepND can co-analyze comorbid neurodevelopmental disorders to discover risk genes

•

The approach employs multitask learning to learn shared and disorder-specific weights

•

DeepND uses graph convolution to process gene interactions in multiple networks

•

The model includes a mixture-of-experts model to detect informative networks

While risk-gene-discovery algorithms have complemented exome/genome-sequencing studies of neurodevelopmental disorders, they are not capable of co-analyzing multiple comorbid conditions like autism and intellectual disability. A common approach is analyzing disorders one by one and comparing the outcomes. With this approach, the method does not utilize cross-disorder interactions and is bound by limited evidence per disorder. We address this gap with a technique, Deep Neurodevelopmental Disorders (DeepND), that uses multitask learning to co-analyze data from multiple disorders to learn shared and disorder-specific patterns. DeepND includes graph convolutional neural networks that process gene-interaction information from multiple networks. DeepND also learns which networks are important for disorder etiologies. Based on this, we propose an interpretable risk-gene-discovery algorithm for neuropsychiatric disorders.

MECP2-related pathways are dysregulated in a cortical organoid model of myotonic dystrophy

Myotonic dystrophy type 1 (DM1) is a multisystem, autosomal-dominant inherited disorder caused by CTG microsatellite repeat expansions (MREs) in the 3' untranslated region of the dystrophia myotonica-protein kinase (DMPK) gene. Despite its prominence as the most common adult-onset muscular dystrophy, patients with congenital to juvenile-onset forms of DM1 can present with debilitating neurocognitive symptoms along the autism spectrum, characteristic of possible in utero cortical defects. However, the molecular mechanism by which CTG MREs lead to these developmental central nervous system (CNS) manifestations is unknown. Here, we showed that CUG foci found early in the maturation of three-dimensional (3D) cortical organoids from DM1 patient-derived induced pluripotent stem cells (iPSCs) cause hyperphosphorylation of CUGBP Elav-like family member 2 (CELF2) protein. Integrative single-cell RNA sequencing and enhanced cross-linking and immunoprecipitation (eCLIP) analysis revealed that reduced CELF2 protein-RNA substrate interactions results in misregulation of genes critical for excitatory synaptic signaling in glutamatergic neurons, including key components of the methyl-CpG binding protein 2 (MECP2) pathway. Comparisons to MECP2(y/-) cortical organoids revealed convergent molecular and cellular defects such as glutamate toxicity and neuronal loss. Our findings provide evidence suggesting that early-onset DM1 might involve neurodevelopmental disorder-associated pathways and identify N-methyl-d-aspartic acid (NMDA) antagonists as potential treatment avenues for neuronal defects in DM1.

Complex regulation of Gephyrin splicing is a determinant of inhibitory postsynaptic diversity

Gephyrin (GPHN) regulates the clustering of postsynaptic components at inhibitory synapses and is involved in pathophysiology of neuropsychiatric disorders. Here, we uncover an extensive diversity of GPHN transcripts that are tightly controlled by splicing during mouse and human brain development. Proteomic analysis reveals at least a hundred isoforms of GPHN incorporated at inhibitory Glycine and gamma-aminobutyric acid A receptors containing synapses. They exhibit different localization and postsynaptic clustering properties, and altering the expression level of one isoform is sufficient to affect the number, size, and density of inhibitory synapses in cerebellar Purkinje cells. Furthermore, we discovered that splicing defects reported in neuropsychiatric disorders are carried by multiple alternative GPHN transcripts, demonstrating the need for a thorough analysis of the GPHN transcriptome in patients. Overall, we show that alternative splicing of GPHN is an important genetic variation to consider in neurological diseases and a determinant of the diversity of postsynaptic inhibitory synapses.

The protein gephyrin is involved in organizing synapses. Here, the authors show how different transcripts of gephyrin form and regulate inhibitory synapses.

Ribonomics Approaches to Identify RBPome in Plants and Other Eukaryotes: Current Progress and Future Prospects

RNA-binding proteins (RBPs) form complex interactions with RNA to regulate the cell’s activities including cell development and disease resistance. RNA-binding proteome (RBPome) aims to profile and characterize the RNAs and proteins that interact with each other to carry out biological functions. Generally, RNA-centric and protein-centric ribonomic approaches have been successfully developed to profile RBPome in different organisms including plants and animals. Further, more and more novel methods that were firstly devised and applied in mammalians have shown great potential to unravel RBPome in plants such as RNA-interactome capture (RIC) and orthogonal organic phase separation (OOPS). Despise the development of various robust and state-of-the-art ribonomics techniques, genome-wide RBP identifications and characterizations in plants are relatively fewer than those in other eukaryotes, indicating that ribonomics techniques have great opportunities in unraveling and characterizing the RNA–protein interactions in plant species. Here, we review all the available approaches for analyzing RBPs in living organisms. Additionally, we summarize the transcriptome-wide approaches to characterize both the coding and non-coding RBPs in plants and the promising use of RBPome for booming agriculture.

Nucleotide-amino acid π-stacking interactions initiate photo cross-linking in RNA-protein complexes

Photo-induced cross-linking is a mainstay technique to characterize RNA-protein interactions. However, UV-induced cross-linking between RNA and proteins at “zero-distance” is poorly understood. Here, we investigate cross-linking of the RBFOX alternative splicing factor with its hepta-ribonucleotide binding element as a model system. We examine the influence of nucleobase, nucleotide position and amino acid composition using CLIR-MS technology (crosslinking-of-isotope-labelled-RNA-and-tandem-mass-spectrometry), that locates cross-links on RNA and protein with site-specific resolution. Surprisingly, cross-linking occurs only at nucleotides that are π-stacked to phenylalanines. Notably, this π-stacking interaction is also necessary for the amino-acids flanking phenylalanines to partake in UV-cross-linking. We confirmed these observations in several published datasets where cross-linking sites could be mapped to a high resolution structure. We hypothesize that π-stacking to aromatic amino acids activates cross-linking in RNA-protein complexes, whereafter nucleotide and peptide radicals recombine. These findings will facilitate interpretation of cross-linking data from structural studies and from genome-wide datasets generated using CLIP (cross-linking-and-immunoprecipitation) methods.

Although UV-induced cross-linking is a widely used method to study RNA-protein complexes, the cross-linking reactions are poorly understood. Here, the authors show that π-stacking interactions between nucleobases and aromatic amino acids play a key role in the cross-linking process.

CircMiMi: a stand-alone software for constructing circular RNA-microRNA-mRNA interactions across species

Background

Circular RNAs (circRNAs) are a class of non-coding RNAs formed by pre-mRNA back-splicing, which are widely expressed in animal/plant cells and often play an important role in regulating microRNA (miRNA) activities. While numerous databases have collected a large amount of predicted circRNA candidates and provided the corresponding circRNA-regulated interactions, a stand-alone package for constructing circRNA-miRNA-mRNA interactions based on user-identified circRNAs across species is lacking.

Results

We present CircMiMi (circRNA-miRNA-mRNA interactions), a modular, Python-based software to identify circRNA-miRNA-mRNA interactions across 18 species (including 16 animals and 2 plants) with the given coordinates of circRNA junctions. The CircMiMi-constructed circRNA-miRNA-mRNA interactions are derived from circRNA-miRNA and miRNA-mRNA axes with the support of computational predictions and/or experimental data. CircMiMi also allows users to examine alignment ambiguity of back-splice junctions for checking circRNA reliability and examine reverse complementary sequences residing in the sequences flanking the circularized exons for investigating circRNA formation. We further employ CircMiMi to identify circRNA-miRNA-mRNA interactions based on the circRNAs collected in NeuroCirc, a large-scale database of circRNAs in the human brain. We construct circRNA-miRNA-mRNA interactions comprising differentially expressed circRNAs, and miRNAs in autism spectrum disorder (ASD) and cross-species analyze the relevance of the targets to ASD. We thus provide a rich set of ASD-associated circRNA-miRNA-mRNA axes and a useful starting point for investigation of regulatory mechanisms in ASD pathophysiology.

Conclusions

CircMiMi allows users to identify circRNA-mediated interactions in multiple species, shedding light on regulatory roles of circRNAs. The software package and web interface are freely available at https://github.com/TreesLab/CircMiMi and http://circmimi.genomics.sinica.edu.tw/, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04692-0.

Leveraging Gene-Level Prediction as Informative Covariate in Hypothesis Weighting Improves Power for Rare Variant Association Studies

Gene-based rare variant association studies (RVASs) have low power due to the infrequency of rare variants and the large multiple testing burden. To correct for multiple testing, traditional false discovery rate (FDR) procedures which depend solely on P-values are often used. Recently, Independent Hypothesis Weighting (IHW) was developed to improve the detection power while maintaining FDR control by leveraging prior information for each hypothesis. Here, we present a framework to increase power of gene-based RVASs by incorporating prior information using IHW. We first build supervised machine learning models to assign each gene a prediction score that measures its disease risk, using the input of multiple biological features, fed with high-confidence risk genes and local background genes selected near GWAS significant loci as the training set. Then we use the prediction scores as covariates to prioritize RVAS results via IHW. We demonstrate the effectiveness of this framework through applications to RVASs in schizophrenia and autism spectrum disorder. We found sizeable improvements in the number of significant associations compared to traditional FDR approaches, and independent evidence supporting the relevance of the genes identified by our framework but not traditional FDR, demonstrating the potential of our framework to improve power of gene-based RVASs.

The m6A reading protein YTHDF3 potentiates tumorigenicity of cancer stem-like cells in ocular melanoma through facilitating CTNNB1 translation

N6-methyladenosine (m6A) is the most universal internal RNA modification on messenger RNAs and regulates the fate and functions of m6A-modified transcripts through m6A-specific binding proteins. Nevertheless, the functional role and potential mechanism of the m6A reading proteins in ocular melanoma tumorigenicity, especially cancer stem-like cell (CSC) properties, remain to be elucidated. Herein, we demonstrated that the m6A reading protein YTHDF3 promotes the translation of the target transcript CTNNB1, contributing to ocular melanoma propagation and migration through m6A methylation. YTHDF3 is highly expressed in ocular melanoma stem-like cells and abundantly enriched in ocular melanoma tissues, which is related to poor clinical prognosis. Moreover, YTHDF3 is required for the maintenance of CSC properties and tumor initiation capacity in ocular melanoma both in vitro and in vivo. Ocular melanoma cells with targeted YTHDF3 knockdown exhibited inhibitory tumor proliferation and migration abilities. Transcriptome-wide mapping of m6A peaks and YTHDF3 binding peaks on mRNAs revealed a key target gene candidate, CTNNB1. Mechanistically, YTHDF3 enhances CTNNB1 translation through recognizing and binding the m6A peaks on CTNNB1 mRNA.

YTHDC2 control of gametogenesis requires helicase activity but not m6A binding

Mechanisms regulating meiotic progression in mammals are poorly understood. The N⁶-methyladenosine (m⁶A) reader and 3' → 5' RNA helicase YTHDC2 switches cells from mitotic to meiotic gene expression programs and is essential for meiotic entry, but how this critical cell fate change is accomplished is unknown. Here, we provide insight into its mechanism and implicate YTHDC2 in having a broad role in gene regulation during multiple meiotic stages. Unexpectedly, mutation of the m⁶A-binding pocket of YTHDC2 had no detectable effect on gametogenesis and mouse fertility, suggesting that YTHDC2 function is m⁶A-independent. Supporting this conclusion, CLIP data defined YTHDC2-binding sites on mRNA as U-rich and UG-rich motif-containing regions within 3' UTRs and coding sequences, distinct from the sites that contain m⁶A during spermatogenesis. Complete loss of YTHDC2 during meiotic entry did not substantially alter translation of its mRNA binding targets in whole-testis ribosome profiling assays but did modestly affect their steady-state levels. Mutation of the ATPase motif in the helicase domain of YTHDC2 did not affect meiotic entry, but it blocked meiotic prophase I progression, causing sterility. Our findings inform a model in which YTHDC2 binds transcripts independent of m⁶A status and regulates gene expression during multiple stages of meiosis by distinct mechanisms.

A matter of space and time: Emerging roles of disease-associated proteins in neural development

Recent genetic studies of neurodevelopmental disorders point to synaptic proteins and ion channels as key contributors to disease pathogenesis. Although many of these proteins, such as the L-type calcium channel Cav1.2 or the postsynaptic scaffolding protein SHANK3, have well-studied functions in mature neurons, new evidence indicates that they may subserve novel, distinct roles in immature cells as the nervous system is assembled in prenatal development. Emerging tools and technologies, including single-cell sequencing and human cellular models of disease, are illuminating differential isoform utilization, spatiotemporal expression, and subcellular localization of ion channels and synaptic proteins in the developing brain compared with the adult, providing new insights into the regulation of developmental processes. We propose that it is essential to consider the temporally distinct and cell-specific roles of these proteins during development and maturity in our framework for understanding neuropsychiatric disorders.

De novo mutations identified by whole-genome sequencing implicate chromatin modifications in obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a chronic anxiety disorder with a substantial genetic basis and a broadly undiscovered etiology. Recent studies of de novo mutation (DNM) exome-sequencing studies for OCD have reinforced the hypothesis that rare variation contributes to the risk. We performed, to our knowledge, the first whole-genome sequencing on 53 parent-offspring families with offspring affected with OCD to investigate all rare de novo variants and insertions/deletions. We observed higher mutation rates in promoter-anchored chromatin loops (empirical P = 0.0015) and regions with high frequencies of histone marks (empirical P = 0.0001). Mutations affecting coding regions were significantly enriched within coexpression modules of genes involved in chromatin modification during human brain development. Four genes—SETD5, KDM3B, ASXL3, and FBL—had strong aggregated evidence and functionally converged on transcription’s epigenetic regulation, suggesting an important OCD risk mechanism. Our data characterized different genome-wide DNMs and highlighted the contribution of chromatin modification in the etiology of OCD.

Trans-genetic effects of circular RNA expression quantitative trait loci and potential causal mechanisms in autism

Genetic risk variants and transcriptional expression changes in autism spectrum disorder (ASD) were widely investigated, but their causal relationship remains largely unknown. Circular RNAs (circRNAs) are abundant in brain and often serve as upstream regulators of mRNAs. By integrating RNA-sequencing with genotype data from autistic brains, we assessed expression quantitative trait loci of circRNAs (circQTLs) that cis-regulated expression of nearby circRNAs and trans-regulated expression of distant genes (trans-eGenes) simultaneously. We thus identified 3619 circQTLs that were also trans-eQTLs and constructed 19,804 circQTL-circRNA-trans-eGene regulatory axes. We conducted two different types of approaches, mediation and partial correlation tests (MPT), to determine the axes with mediation effects of circQTLs on trans-eGene expression through circRNA expression. We showed that the mediation effects of the circQTLs (trans-eQTLs) on circRNA expression were positively correlated with the magnitude of circRNA-trans-eGene correlation of expression profile. The positive correlation became more significant after adjustment for the circQTLs. Of the 19,804 axes, 8103 passed MPT. Meanwhile, we performed causal inference test (CIT) and identified 2070 circQTL-trans-eGene-ASD diagnosis propagation paths. We showed that the CIT-passing genes were significantly enriched for ASD risk genes, genes encoding postsynaptic density proteins, and other ASD-relevant genes, supporting the relevance of the CIT-passing genes to ASD pathophysiology. Integration of MPT- and CIT-passing axes further constructed 352 circQTL-circRNA-trans-eGene-ASD diagnosis propagation paths, wherein the circRNA-trans-eGene axes may act as causal mediators for the circQTL-ASD diagnosis associations. These analyses were also successfully applied to an independent dataset from schizophrenia brains. Collectively, this study provided the first framework for systematically investigating trans-genetic effects of circQTLs and inferring the corresponding causal relations in diseases. The identified circQTL-circRNA-trans-eGene regulatory interactions, particularly the internal modules that were previously implicated in the examined disorders, also provided a helpful dataset for further investigating causative biology and cryptic regulatory mechanisms underlying the neuropsychiatric diseases.

A survey of RNA editing at single-cell resolution links interneurons to schizophrenia and autism

Conversion of adenosine to inosine in RNA by ADAR enzymes, termed "RNA editing," is essential for healthy brain development. Editing is dysregulated in neuropsychiatric diseases, but has not yet been investigated at scale at the level of individual neurons. We quantified RNA editing sites in nuclear transcriptomes of 3055 neurons from six cortical regions of a neurotypical female donor, and found 41,930 sites present in at least ten nuclei. Most sites were located within Alu repeats in introns or 3' UTRs, and approximately 80% were cataloged in public RNA editing databases. We identified 9285 putative novel editing sites, 29% of which were also detectable in unrelated donors. Intersection with results from bulk RNA-seq studies provided cell-type and spatial context for 1730 sites that are differentially edited in schizophrenic brain donors, and 910 such sites in autistic donors. Autism-related genes were also enriched with editing sites predicted to modify RNA structure. Inhibitory neurons showed higher overall transcriptome editing than excitatory neurons, and the highest editing rates were observed in the frontal cortex. We used generalized linear models to identify differentially edited sites and genes between cell types. Twenty nine genes were preferentially edited in excitatory neurons, and 43 genes were edited more heavily in inhibitory neurons, including RBFOX1, its target genes, and genes in the autism-associated Prader-Willi locus (15q11). The abundance of SNORD115/116 genes from locus 15q11 was positively associated with editing activity across the transcriptome. We contend that insufficient editing of autism-related genes in inhibitory neurons may contribute to the specific perturbation of those cells in autism.

Postsynaptic autism spectrum disorder genes and synaptic dysfunction

This review provides an overview of the synaptic dysfunction of neuronal circuits and the ensuing behavioral alterations caused by mutations in autism spectrum disorder (ASD)-linked genes directly or indirectly affecting the postsynaptic neuronal compartment. There are plenty of ASD risk genes, that may be broadly grouped into those involved in gene expression regulation (epigenetic regulation and transcription) and genes regulating synaptic activity (neural communication and neurotransmission). Notably, the effects mediated by ASD-associated genes can vary extensively depending on the developmental time and/or subcellular site of expression. Therefore, in order to gain a better understanding of the mechanisms of disruptions in postsynaptic function, an effort to better model ASD in experimental animals is required to improve standardization and increase reproducibility within and among studies. Such an effort holds promise to provide deeper insight into the development of these disorders and to improve the translational value of preclinical studies.

Pleiotropy between language impairment and broader behavioral disorders-an investigation of both common and rare genetic variants

BACKGROUND

Language plays a major role in human behavior. For this reason, neurodevelopmental and psychiatric disorders in which linguistic ability is impaired could have a big impact on the individual's social interaction and general wellbeing. Such disorders tend to have a strong genetic component, but most past studies examined mostly the linguistic overlaps across these disorders; investigations into their genetic overlaps are limited. The aim of this study was to assess the potential genetic overlap between language impairment and broader behavioral disorders employing methods capturing both common and rare genetic variants.

METHODS

We employ polygenic risk scores (PRS) trained on specific language impairment (SLI) to evaluate genetic overlap across several disorders in a large case-cohort sample comprising ~13,000 autism spectrum disorder (ASD) cases, including cases of childhood autism and Asperger's syndrome, ~15,000 attention deficit/hyperactivity disorder (ADHD) cases, ~3000 schizophrenia cases, and ~21,000 population controls. We also examine rare variants in SLI/language-related genes in a subset of the sample that was exome-sequenced using the SKAT-O method.

RESULTS

We find that there is little evidence for genetic overlap between SLI and ADHD, schizophrenia, and ASD, the latter being in line with results of linguistic analyses in past studies. However, we observe a small, significant genetic overlap between SLI and childhood autism specifically, which we do not observe for SLI and Asperger's syndrome. Moreover, we observe that childhood autism cases have significantly higher SLI-trained PRS compared to Asperger's syndrome cases; these results correspond well to the linguistic profiles of both disorders. Our rare variant analyses provide suggestive evidence of association for specific genes with ASD, childhood autism, and schizophrenia.

CONCLUSIONS

Our study provides, for the first time, to our knowledge, genetic evidence for ASD subtypes based on risk variants for language impairment.

RBFOX2 is critical for maintaining alternative polyadenylation patterns and mitochondrial health in rat myoblasts

RBFOX2, which has a well-established role in alternative splicing, is linked to heart diseases. However, it is unclear whether RBFOX2 has other roles in RNA processing that can influence gene expression in muscle cells, contributing to heart disease. Here, we employ both 3'-end and nanopore cDNA sequencing to reveal a previously unrecognized role for RBFOX2 in maintaining alternative polyadenylation (APA) signatures in myoblasts. RBFOX2-mediated APA modulates mRNA levels and/or isoform expression of a collection of genes, including contractile and mitochondrial genes. Depletion of RBFOX2 adversely affects mitochondrial health in myoblasts, correlating with disrupted APA of mitochondrial gene Slc25a4. Mechanistically, RBFOX2 regulation of Slc25a4 APA is mediated through consensus RBFOX2 binding motifs near the distal polyadenylation site, enforcing the use of the proximal polyadenylation site. In sum, our results unveil a role for RBFOX2 in fine-tuning expression of mitochondrial and contractile genes via APA in myoblasts relevant to heart diseases.