CELF4 regulates translation and local abundance of a vast set of mRNAs, including genes associated with regulation of synaptic function

The role of CELF family in neurodevelopment and neurodevelopmental disorders

RNA-binding proteins (RBPs) bind to RNAs and are crucial for regulating RNA splicing, stability, translation, and transport. Among these proteins, the CUGBP Elav-like family (CELF) is a highly conserved group crucial for posttranscriptional regulation by binding to CUG repeats. Comprising CELF1-6, this family exhibits diverse expression patterns and functions. Dysregulation of CELF has been implicated in various neural disorders, encompassing both neurodegenerative and neurodevelopmental conditions, such as Alzheimer's disease and autism. This article aims to provide a comprehensive summary of the CELF family's role in neurodevelopment and neurodevelopmental disorders. Understanding CELF's mechanisms may offer clues for potential therapeutic strategies by regulating their targets in neurodevelopmental disorders.

Altered CELF4 splicing factor enhances pancreatic neuroendocrine tumors aggressiveness influencing mTOR and everolimus response

Pancreatic neuroendocrine tumors (PanNETs) comprise a heterogeneous group of tumors with growing incidence. Recent molecular analyses provided a precise picture of their genomic and epigenomic landscape. Splicing dysregulation is increasingly regarded as a novel cancer hallmark influencing key tumor features. We have previously demonstrated that splicing machinery is markedly dysregulated in PanNETs. Here, we aimed to elucidate the molecular and functional implications of CUGBP ELAV-like family member 4 (CELF4), one of the most altered splicing factors in PanNETs. CELF4 expression was determined in 20 PanNETs, comparing tumor and non-tumoral adjacent tissue. An RNA sequencing (RNA-seq) dataset was analyzed to explore CELF4-linked interrelations among clinical features, gene expression, and splicing events. Two PanNET cell lines were employed to assess CELF4 function in vitro and in vivo. PanNETs display markedly upregulated CELF4 expression, which is closely associated with malignancy features, altered expression of key tumor players, and distinct splicing event profiles. Modulation of CELF4 influenced proliferation in vitro and reduced in vivo xenograft tumor growth. Interestingly, functional assays and RNA-seq analysis revealed that CELF4 silencing altered mTOR signaling pathway, enhancing the effect of everolimus. We demonstrate that CELF4 is dysregulated in PanNETs, where it influences tumor development and aggressiveness, likely by modulating the mTOR pathway, suggesting its potential as therapeutic target.

Genetic architecture of the structural connectome

Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.

RNA-binding proteins in pain

RNAs are meticulously controlled by proteins. Through direct and indirect associations, every facet in the brief life of an mRNA is subject to regulation. RNA-binding proteins (RBPs) permeate biology. Here, we focus on their roles in pain. Chronic pain is among the largest challenges facing medicine and requires new strategies. Mounting pharmacologic and genetic evidence obtained in pre-clinical models suggests fundamental roles for a broad array of RBPs. We describe their diverse roles that span RNA modification, splicing, stability, translation, and decay. Finally, we highlight opportunities to expand our understanding of regulatory interactions that contribute to pain signaling. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Regulation RNA in Disease and Development > RNA in Disease.

The putative RNA helicase DDX1 associates with the nuclear RNA exosome and modulates RNA/DNA hybrids (R-loops)

The RNA exosome is a ribonuclease complex that mediates both RNA processing and degradation. This complex is evolutionarily conserved, ubiquitously expressed, and required for fundamental cellular functions, including rRNA processing. The RNA exosome plays roles in regulating gene expression and protecting the genome, including modulating the accumulation of RNA-DNA hybrids (R-loops). The function of the RNA exosome is facilitated by cofactors, such as the RNA helicase MTR4, which binds/remodels RNAs. Recently, missense mutations in RNA exosome subunit genes have been linked to neurological diseases. One possibility to explain why missense mutations in genes encoding RNA exosome subunits lead to neurological diseases is that the complex may interact with cell- or tissue-specific cofactors that are impacted by these changes. To begin addressing this question, we performed immunoprecipitation of the RNA exosome subunit, EXOSC3, in a neuronal cell line (N2A), followed by proteomic analyses to identify novel interactors. We identified the putative RNA helicase, DDX1, as an interactor. DDX1 plays roles in double-strand break repair, rRNA processing, and R-loop modulation. To explore the functional connections between EXOSC3 and DDX1, we examined the interaction following double-strand breaks and analyzed changes in R-loops in N2A cells depleted for EXOSC3 or DDX1 by DNA/RNA immunoprecipitation followed by sequencing. We find that EXOSC3 interaction with DDX1 is decreased in the presence of DNA damage and that loss of EXOSC3 or DDX1 alters R-loops. These results suggest EXOSC3 and DDX1 interact during events of cellular homeostasis and potentially suppress unscrupulous expression of genes promoting neuronal projection.

Hypermethylated CDO1 and CELF4 in cytological specimens as triage strategy biomarkers in endometrial malignant lesions

Objective

Developing a non-invasive and reliable triage test for endometrial malignant lesions is an important goal, as it could help to reduce the number of invasive diagnostic procedures required and improve patient survival. We aimed to estimate the diagnostic value of DNA methylation levels in cervical cytological samples of endometrial cancer (EC) and endometrial atypical hyperplasia (AH).

Methods

A total of 607 women who had indications for endometrial biopsy in the Department of Obstetrics and Gynecology of Cangzhou Central Hospital from October 2022 to April 2023 were enrolled in this study. The cervical exfoliated cells were collected for gene methylation before endometrial biopsy. Clinical information, tumor biomarkers, and endometrial thickness (ET) of transvaginal ultrasonography (TVS) were also collected. With endometrial histopathology as the gold standard, multivariate unconditional logistic regression was applied to analyze the risk factors of endometrial malignant lesions. The role of cysteine dioxygenase type 1 (CDO1) and CUGBP Elav-like family member 4 (CELF4) gene methylation as a triage strategy biomarker in endometrial malignant lesions was specifically explored.

Results

Multivariate logistic regression analysis showed that premenopausal ET ≥ 11 mm or postmenopausal ET ≥ 5 mm, CDO1 ΔCt ≤ 8.4, or CELF4 ΔCt ≤ 8.8 were the risk factors for AH and EC, with odds ratios (ORs) (95%CI) of 5.03 (1.83-13.82) and 6.92 (1.10-43.44), respectively (p-values < 0.05). The sensitivity and specificity of CDO1/CELF4 dual-gene methylation assay for AH and EC reached 84.9% (95%CI: 75.3%-94.5%) and 86.6% (95%CI: 83.8%-89.5%), respectively. ET combined with DNA methylation detection further improved the specificity to (94.9%, 95%CI: 93.1%-96.8%).

Conclusion

The accuracy of cervical cytology DNA methylation is superior to that of other clinical indicators in the non-invasive examination of endometrial malignant lesions. DNA methylation combined with TVS can further improve the specificity and is a promising biomarker triage strategy in women with suspected endometrial lesions.

Members of the CUGBP Elav-like family of RNA-binding proteins are expressed in distinct populations of primary sensory neurons

Primary sensory dorsal root ganglia (DRG) neurons are diverse, with distinct populations that respond to specific stimuli. Previously, we observed that functionally distinct populations of DRG neurons express mRNA transcript variants with different 3' untranslated regions (3'UTRs). 3'UTRs harbor binding sites for interaction with RNA-binding proteins (RBPs) for transporting mRNAs to subcellular domains, modulating transcript stability, and regulating the rate of translation. In the current study, analysis of publicly available single-cell RNA-sequencing data generated from adult mice revealed that 17 3'UTR-binding RBPs were enriched in specific populations of DRG neurons. This included four members of the CUG triplet repeat (CUGBP) Elav-like family (CELF): CELF2 and CELF4 were enriched in peptidergic, CELF6 in both peptidergic and nonpeptidergic, and CELF3 in tyrosine hydroxylase-expressing neurons. Immunofluorescence studies confirmed that 60% of CELF4+ neurons are small-diameter C fibers and 33% medium-diameter myelinated (likely Aδ) fibers and showed that CELF4 is distributed to peripheral termini. Coexpression analyses using transcriptomic data and immunofluorescence revealed that CELF4 is enriched in nociceptive neurons that express GFRA3, CGRP, and the capsaicin receptor TRPV1. Reanalysis of published transcriptomic data from macaque DRG revealed a highly similar distribution of CELF members, and reanalysis of single-nucleus RNA-sequencing data derived from mouse and rat DRG after sciatic injury revealed differential expression of CELFs in specific populations of sensory neurons. We propose that CELF RBPs may regulate the fate of mRNAs in populations of nociceptors, and may play a role in pain and/or neuronal regeneration following nerve injury.

Celf4 controls mRNA translation underlying synaptic development in the prenatal mammalian neocortex

Abnormalities in neocortical and synaptic development are linked to neurodevelopmental disorders. However, the molecular and cellular mechanisms governing initial synapse formation in the prenatal neocortex remain poorly understood. Using polysome profiling coupled with snRNAseq on human cortical samples at various fetal phases, we identify human mRNAs, including those encoding synaptic proteins, with finely controlled translation in distinct cell populations of developing frontal neocortices. Examination of murine and human neocortex reveals that the RNA binding protein and translational regulator, CELF4, is expressed in compartments enriched in initial synaptogenesis: the marginal zone and the subplate. We also find that Celf4/CELF4-target mRNAs are encoded by risk genes for adverse neurodevelopmental outcomes translating into synaptic proteins. Surprisingly, deleting Celf4 in the forebrain disrupts the balance of subplate synapses in a sex-specific fashion. This highlights the significance of RNA binding proteins and mRNA translation in evolutionarily advanced synaptic development, potentially contributing to sex differences.

The predicted RNA-binding protein regulome of axonal mRNAs

Neurons are morphologically complex cells that rely on the compartmentalization of protein expression to develop and maintain their cytoarchitecture. The targeting of RNA transcripts to axons is one of the mechanisms that allows rapid local translation of proteins in response to extracellular signals. 3' Untranslated regions (UTRs) of mRNA are noncoding sequences that play a critical role in determining transcript localization and translation by interacting with specific RNA-binding proteins (RBPs). However, how 3' UTRs contribute to mRNA metabolism and the nature of RBP complexes responsible for these functions remains elusive. We performed 3' end sequencing of RNA isolated from cell bodies and axons of sympathetic neurons exposed to either nerve growth factor (NGF) or neurotrophin 3 (NTF3, also known as NT-3). NGF and NTF3 are growth factors essential for sympathetic neuron development through distinct signaling mechanisms. Whereas NTF3 acts mostly locally, NGF signal is retrogradely transported from axons to cell bodies. We discovered that both NGF and NTF3 affect transcription and alternative polyadenylation in the nucleus and induce the localization of specific 3' UTR isoforms to axons, including short 3' UTR isoforms found exclusively in axons. The integration of our data with CLIP sequencing data supports a model whereby long 3' UTR isoforms associate with RBP complexes in the nucleus and, upon reaching the axons, are remodeled locally into shorter isoforms. Our findings shed new light into the complex relationship between nuclear polyadenylation, mRNA localization, and local 3' UTR remodeling in developing neurons.

Molecular features driving cellular complexity of human brain evolution

Human-specific genomic changes contribute to the unique functionalities of the human brain1-5. The cellular heterogeneity of the human brain6,7 and the complex regulation of gene expression highlight the need to characterize human-specific molecular features at cellular resolution. Here we analysed single-nucleus RNA-sequencing and single-nucleus assay for transposase-accessible chromatin with sequencing datasets for human, chimpanzee and rhesus macaque brain tissue from posterior cingulate cortex. We show a human-specific increase of oligodendrocyte progenitor cells and a decrease of mature oligodendrocytes across cortical tissues. Human-specific regulatory changes were accelerated in oligodendrocyte progenitor cells, and we highlight key biological pathways that may be associated with the proportional changes. We also identify human-specific regulatory changes in neuronal subtypes, which reveal human-specific upregulation of FOXP2 in only two of the neuronal subtypes. We additionally identify hundreds of new human accelerated genomic regions associated with human-specific chromatin accessibility changes. Our data also reveal that FOS::JUN and FOX motifs are enriched in the human-specifically accessible chromatin regions of excitatory neuronal subtypes. Together, our results reveal several new mechanisms underlying the evolutionary innovation of human brain at cell-type resolution.

The RNA helicase DDX1 associates with the nuclear RNA exosome and modulates R-loops

The RNA exosome is a ribonuclease complex that mediates both RNA processing and degradation. This complex is evolutionarily conserved, ubiquitously expressed, and required for fundamental cellular functions, including rRNA processing. The RNA exosome plays roles in regulating gene expression and protecting the genome, including modulating the accumulation of RNA-DNA hybrids (R-loops). The function of the RNA exosome is facilitated by cofactors, such as the RNA helicase MTR4, which binds/remodels RNAs. Recently, missense mutations in RNA exosome subunit genes have been linked to neurological diseases. One possibility to explain why missense mutations in genes encoding RNA exosome subunits lead to neurological diseases is that the complex may interact with cell- or tissue-specific cofactors that are impacted by these changes. To begin addressing this question, we performed immunoprecipitation of the RNA exosome subunit, EXOSC3, in a neuronal cell line (N2A) followed by proteomic analyses to identify novel interactors. We identified the putative RNA helicase, DDX1, as an interactor. DDX1 plays roles in double-strand break repair, rRNA processing, and R-loop modulation. To explore the functional connections between EXOSC3 and DDX1, we examined the interaction following double-strand breaks, and analyzed changes in R-loops in N2A cells depleted for EXOSC3 or DDX1 by DNA/RNA immunoprecipitation followed by sequencing (DRIP-Seq). We find that EXOSC3 interaction with DDX1 is decreased in the presence of DNA damage and that loss of EXOSC3 or DDX1 alters R-loops. These results suggest EXOSC3 and DDX1 interact during events of cellular homeostasis and potentially suppress unscrupulous expression of genes promoting neuronal projection.

Adaptation of the Porcine Pituitary Transcriptome, Spliceosome and Editome during Early Pregnancy

The physiological mechanisms of the porcine reproduction are relatively well-known. However, transcriptomic changes and the mechanisms accompanying transcription and translation processes in various reproductive organs, as well as their dependence on hormonal status, are still poorly understood. The aim of this study was to gain a principal understanding of alterations within the transcriptome, spliceosome and editome occurring in the pituitary of the domestic pig (Sus scrofa domestica L.), which controls basic physiological processes in the reproductive system. In this investigation, we performed extensive analyses of data obtained by high-throughput sequencing of RNA from the gilts' pituitary anterior lobes during embryo implantation and the mid-luteal phase of the estrous cycle. During analyses, we obtained detailed information on expression changes of 147 genes and 43 long noncoding RNAs, observed 784 alternative splicing events and also found the occurrence of 8729 allele-specific expression sites and 122 RNA editing events. The expression profiles of the selected 16 phenomena were confirmed by PCR or qPCR techniques. As a final result of functional meta-analysis, we acquired knowledge regarding intracellular pathways that induce changes in the processes accompanying transcription and translation regulation, which may induce modifications in the secretory activity of the porcine adenohypophyseal cells.

Molecular pathways of major depressive disorder converge on the synapse

Major depressive disorder (MDD) is a psychiatric disease of still poorly understood molecular etiology. Extensive studies at different molecular levels point to a high complexity of numerous interrelated pathways as the underpinnings of depression. Major systems under consideration include monoamines, stress, neurotrophins and neurogenesis, excitatory and inhibitory neurotransmission, mitochondrial dysfunction, (epi)genetics, inflammation, the opioid system, myelination, and the gut-brain axis, among others. This review aims at illustrating how these multiple signaling pathways and systems may interact to provide a more comprehensive view of MDD's neurobiology. In particular, considering the pattern of synaptic activity as the closest physical representation of mood, emotion, and conscience we can conceptualize, each pathway or molecular system will be scrutinized for links to synaptic neurotransmission. Models of the neurobiology of MDD will be discussed as well as future actions to improve the understanding of the disease and treatment options.

RNA splicing regulators play critical roles in neurogenesis

Alternative RNA splicing increases transcript diversity in different cell types and under varying conditions. It is executed with the help of RNA splicing regulators (RSRs), which are operationally defined as RNA-binding proteins (RBPs) that regulate alternative splicing, but not directly catalyzing the chemical reactions of splicing. By systematically searching for RBPs and manually identifying those that regulate splicing, we curated 305 RSRs in the human genome. Surprisingly, most of the RSRs are involved in neurogenesis. Among these RSRs, we focus on nine families (PTBP, NOVA, RBFOX, ELAVL, CELF, DBHS, MSI, PCBP, and MBNL) that play essential roles in the neurogenic pathway. A better understanding of their functions will provide novel insights into the role of splicing in brain development, health, and disease. This comprehensive review serves as a stepping-stone to explore the diverse and complex set of RSRs as fundamental regulators of neural development. This article is categorized under: RNA-Based Catalysis &gt; RNA Catalysis in Splicing and Translation RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Interactions: Functional Implications RNA Processing &gt; Splicing Regulation/Alternative Splicing.

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.

Krüppel-like Transcription Factor 7 Is a Causal Gene in Autism Development

Background: Autism spectrum disorder (ASD) is a complex neurodevelopmental disease. To date, more than 1000 genes have been shown to be associated with ASD, and only a few of these genes account for more than 1% of autism cases. Klf7 is an important transcription factor of cell proliferation and differentiation in the nervous system, but whether klf7 is involved in autism is unclear. Methods: We first performed ChIP-seq analysis of klf7 in N2A cells, then performed behavioral tests and RNA-seq in klf7^+/− mice, and finally restored mice with adeno-associated virus (AAV)-mediated overexpression of klf7 in klf7^+/− mice. Results: Klf7 targeted genes are enriched with ASD genes, and 631 ASD risk genes are also differentially expressed in klf7^+/− mice which exhibited the core symptoms of ASD. When klf7 levels were increased in the central nervous system (CNS) in klf7^+/− adult mice, deficits in social interaction, repetitive behavior and majority of dysregulated ASD genes were rescued in the adults, suggesting transcriptional regulation. Moreover, knockdown of klf7 in human brain organoids caused dysregulation of 517 ASD risk genes, 344 of which were shared with klf7^+/− mice, including some high-confidence ASD genes. Conclusions: Our findings highlight a klf7 regulation of ASD genes and provide new insights into the pathogenesis of ASD and promising targets for further research on mechanisms and treatments.

Shared Genetics and Causality Between Decaffeinated Coffee Consumption and Neuropsychiatric Diseases: A Large-Scale Genome-Wide Cross-Trait Analysis and Mendelian Randomization Analysis

Coffee or caffeine consumption has been associated with neuropsychiatric disorders, implying a shared etiology. However, whether these associations reflect causality remains largely unknown. To understand the genetic structure of the association between decaffeinated coffee consumption (DCC) and neuropsychiatric traits, we examined the genetic correlation, causality, and shared genetic structure between DCC and neuropsychiatric traits using linkage disequilibrium score regression, bidirectional Mendelian randomization (MR), and genome-wide cross-trait meta-analysis in large GWAS Consortia for coffee consumption (N = 329,671) and 13 neuropsychiatric traits (sample size ranges from 36,052 to 500,199). We found strong positive genetic correlations between DCC and lifetime cannabis use (LCU; Rg = 0.48, P = 8.40 × 10^-19), alcohol use disorder identification test (AUDIT) total score (AUDIT_T; Rg = 0.40, P = 4.63 × 10^-13), AUDIT_C score (alcohol consumption component of the AUDIT; Rg = 0.40, P = 5.26 × 10^-11), AUDIT_P score (dependence and hazardous-use component of the AUDIT; Rg = 0.28, P = 1.36 × 10^-05), and strong negative genetic correlations between DCC and neuroticism (Rg = -0.15, P = 7.27 × 10^-05), major depressed diseases (MDD; Rg = -0.15, P = 0.0010), and insomnia (Rg= -0.15, P = 0.0007). In the cross-trait meta-analysis, we identified 6, 5, 1, 1, 2, 31, and 27 shared loci between DCC and Insomnia, LCU, AUDIT_T, AUDIT_C, AUDIT_P, neuroticism, and MDD, respectively, which were mainly enriched in bone marrow, lymph node, cervix, uterine, lung, and thyroid gland tissues, T cell receptor signaling pathway, antigen receptor-mediated signaling pathway, and epigenetic pathways. A large of TWAS-significant associations were identified in tissues that are part of the nervous system, digestive system, and exo-/endocrine system. Our findings further indicated a causal influence of liability to DCC on LCU and low risk of MDD (odds ratio: 0.90, P = 9.06 × 10^-5 and 1.27, P = 7.63 × 10^-4 respectively). We also observed that AUDIT_T and AUDIT_C were causally related to DCC (odds ratio: 1.83 per 1-SD increase in AUDIT_T, P = 1.67 × 10^-05, 1.80 per 1-SD increase in AUDIT_C, P = 5.09 × 10^-04). Meanwhile, insomnia and MDD had a causal negative influence on DCC (OR: 0.91, 95% CI: 0.86-0.95, P = 1.51 × 10^-04 for Insomnia; OR: 0.93, 95% CI: 0.89-0.99, P = 6.02 × 10^-04 for MDD). These findings provided evidence for the shared genetic basis and causality between DCC and neuropsychiatric diseases, and advance our understanding of the shared genetic mechanisms underlying their associations, as well as assisting with making recommendations for clinical works or health education.

CELF Family Proteins in Cancer: Highlights on the RNA-Binding Protein/Noncoding RNA Regulatory Axis

Post-transcriptional modifications to coding and non-coding RNAs are unquestionably a pivotal way in which human mRNA and protein diversity can influence the different phases of a transcript's life cycle. CELF (CUGBP Elav-like family) proteins are RBPs (RNA-binding proteins) with pleiotropic capabilities in RNA processing. Their responsibilities extend from alternative splicing and transcript editing in the nucleus to mRNA stability, and translation into the cytoplasm. In this way, CELF family members have been connected to global alterations in cancer proliferation and invasion, leading to their identification as potential tumor suppressors or even oncogenes. Notably, genetic variants, alternative splicing, phosphorylation, acetylation, subcellular distribution, competition with other RBPs, and ultimately lncRNAs, miRNAs, and circRNAs all impact CELF regulation. Discoveries have emerged about the control of CELF functions, particularly via noncoding RNAs, and CELF proteins have been identified as competing, antagonizing, and regulating agents of noncoding RNA biogenesis. On the other hand, CELFs are an intriguing example through which to broaden our understanding of the RBP/noncoding RNA regulatory axis. Balancing these complex pathways in cancer is undeniably pivotal and deserves further research. This review outlines some mechanisms of CELF protein regulation and their functional consequences in cancer physiology.

Expression Biomarkers of Pharmacological Treatment Outcomes in Women with Unipolar and Bipolar Depression

INTRODUCTION

This study aimed to find the expression biomarkers of pharmacological treatment response in a naturalistic hospital setting. Through gene expression profiling, we were able to find differentially-expressed genes (DEGs) in unipolar (UD) and bipolar (BD) depressed women.

METHODS

We performed gene expression profiling in hospitalized women with unipolar (n=24) and bipolar depression (n=32) who achieved clinical improvement after pharmacological treatment (without any restriction). To identify DEGs in peripheral blood mononuclear cells (PBMCs), we used the SurePrint G3 Microarray and GeneSpring software.

RESULTS

After pharmacological treatment, UD and BD varied in the number of regulated genes and ontological pathways. Also, the pathways of neurogenesis and synaptic transmission were significantly up-regulated. Our research focused on DEGs with a minimum fold change (FC) of more than 2. For both types of depression, 2 up-regulated genes, OPRM1 and CELF4 (p=0.013), were significantly associated with treatment response (defined as a 50% reduction on the Hamilton Depression Rating Scale [HDRS]). We also uncovered the SHANK3 (p=0.001) gene that is unique for UD and found that the RASGRF1 (p=0.010) gene may be a potential specific biomarker of treatment response for BD.

CONCLUSION

Based on transcriptomic profiling, we identified potential expression biomarkers of treatment outcomes for UD and BD. We also proved that the Ras-GEF pathway associated with long-term memory, female stress response, and treatment response modulation in animal studies impacts treatment efficacy in patients with BD. Further studies focused on the outlined genes may help provide predictive markers of treatment outcomes in UD and BD.

The Genetic Control of Stoichiometry Underlying Autism

Autism is a common and complex neurologic disorder whose scientific underpinnings have begun to be established in the past decade. The essence of this breakthrough has been a focus on families, where genetic analyses are strongest, versus large-scale, case-control studies. Autism genetics has progressed in parallel with technology, from analyses of copy number variation to whole-exome sequencing (WES) and whole-genome sequencing (WGS). Gene mutations causing complete loss of function account for perhaps one-third of cases, largely detected through WES. This limitation has increased interest in understanding the regulatory variants of genes that contribute in more subtle ways to the disorder. Strategies combining biochemical analysis of gene regulation, WGS analysis of the noncoding genome, and machine learning have begun to succeed. The emerging picture is that careful control of the amounts of transcription, mRNA, and proteins made by key brain genes-stoichiometry-plays a critical role in defining the clinical features of autism.

Novel ultra-rare exonic variants identified in a founder population implicate cadherins in schizophrenia

The identification of rare variants associated with schizophrenia has proven challenging due to genetic heterogeneity, which is reduced in founder populations. In samples from the Ashkenazi Jewish population, we report that schizophrenia cases had a greater frequency of novel missense or loss of function (MisLoF) ultra-rare variants (URVs) compared to controls, and the MisLoF URV burden was inversely correlated with polygenic risk scores in cases. Characterizing 141 "case-only" genes (MisLoF URVs in ≥3 cases with none in controls), the cadherin gene set was associated with schizophrenia. We report a recurrent case mutation in PCDHA3 that results in the formation of cytoplasmic aggregates and failure to engage in homophilic interactions on the plasma membrane in cultured cells. Modeling purifying selection, we demonstrate that deleterious URVs are greatly overrepresented in the Ashkenazi population, yielding enhanced power for association studies. Identification of the cadherin/protocadherin family as risk genes helps specify the synaptic abnormalities central to schizophrenia.

Exploring the biological role of postzygotic and germinal de novo mutations in ASD

De novo mutations (DNMs), including germinal and postzygotic mutations (PZMs), are a strong source of causality for Autism Spectrum Disorder (ASD). However, the biological processes involved behind them remain unexplored. Our aim was to detect DNMs (germinal and PZMs) in a Spanish ASD cohort (360 trios) and to explore their role across different biological hierarchies (gene, biological pathway, cell and brain areas) using bioinformatic approaches. For the majority of the analysis, a combined ASD cohort (N = 2171 trios) was created using previously published data by the Autism Sequencing Consortium (ASC). New plausible candidate genes for ASD such as FMR1 and NFIA were found. In addition, genes harboring PZMs were significantly enriched for miR-137 targets in comparison with germinal DNMs that were enriched in GO terms related to synaptic transmission. The expression pattern of genes with PZMs was restricted to early mid-fetal cortex. In contrast, the analysis of genes with germinal DNMs revealed a spatio-temporal window from early to mid-fetal development stages, with expression in the amygdala, cerebellum, cortex and striatum. These results provide evidence of the pathogenic role of PZMs and suggest the existence of distinct mechanisms between PZMs and germinal DNMs that are influencing ASD risk.

Massively Parallel Reporter Assays: Defining Functional Psychiatric Genetic Variants Across Biological Contexts

Neuropsychiatric phenotypes have long been known to be influenced by heritable risk factors, directly confirmed by the past decade of genetic studies that have revealed specific genetic variants enriched in disease cohorts. However, the initial hope that a small set of genes would be responsible for a given disorder proved false. The more complex reality is that a given disorder may be influenced by myriad small-effect noncoding variants and/or by rare but severe coding variants, many de novo. Noncoding genomic sequences-for which molecular functions cannot usually be inferred-harbor a large portion of these variants, creating a substantial barrier to understanding higher-order molecular and biological systems of disease. Fortunately, novel genetic technologies-scalable oligonucleotide synthesis, RNA sequencing, and CRISPR (clustered regularly interspaced short palindromic repeats)-have opened novel avenues to experimentally identify biologically significant variants en masse. Massively parallel reporter assays (MPRAs) are an especially versatile technique resulting from such innovations. MPRAs are powerful molecular genetics tools that can be used to screen thousands of untranscribed or untranslated sequences and their variants for functional effects in a single experiment. This approach, though underutilized in psychiatric genetics, has several useful features for the field. We review methods for assaying putatively functional genetic variants and regions, emphasizing MPRAs and the opportunities they hold for dissection of psychiatric polygenicity. We discuss literature applying functional assays in neurogenetics, highlighting strengths, caveats, and design considerations-especially regarding disease-relevant variables (cell type, neurodevelopment, and sex), and we ultimately propose applications of MPRA to both computational and experimental neurogenetics of polygenic disease risk.

An endogenous PI3K interactome promoting astrocyte-mediated neuroprotection identifies a novel association with RNA-binding protein ZC3H14

Astrocytes can support neuronal survival through a range of secreted signals that protect against neurotoxicity, oxidative stress, and apoptotic cascades. Thus, analyzing the effects of the astrocyte secretome may provide valuable insight into these neuroprotective mechanisms. Previously, we characterized a potent neuroprotective activity mediated by retinal astrocyte conditioned media (ACM) on retinal and cortical neurons in metabolic stress models. However, the molecular mechanism underlying this complex activity in neuronal cells has remained unclear. Here, a chemical genetics screen of kinase inhibitors revealed phosphoinositide 3-kinase (PI3K) as a central player transducing ACM-mediated neuroprotection. To identify additional proteins contributing to the protective cascade, endogenous PI3K was immunoprecipitated from neuronal cells exposed to ACM or control media, followed by MS/MS proteomic analyses. These data pointed toward a relatively small number of proteins that coimmunoprecipitated with PI3K, and surprisingly only five were regulated by the ACM signal. These hits included expected PI3K interactors, such as the platelet-derived growth factor receptor A (PDGFRA), as well as novel RNA-binding protein interactors ZC3H14 (zinc finger CCCH-type containing 14) and THOC1 (THO complex protein 1). In particular, ZC3H14 has recently emerged as an important RNA-binding protein with multiple roles in posttranscriptional regulation. In validation studies, we show that PI3K recruitment of ZC3H14 is necessary for PDGF-induced neuroprotection and that this interaction is present in primary retinal ganglion cells. Thus, we identified a novel non-cell autonomous neuroprotective signaling cascade mediated through PI3K that requires recruitment of ZC3H14 and may present a promising strategy to promote astrocyte-secreted prosurvival signals.

CLIP and Massively Parallel Functional Analysis of CELF6 Reveal a Role in Destabilizing Synaptic Gene mRNAs through Interaction with 3' UTR Elements

CELF6 is a CELF-RNA-binding protein, and thus part of a protein family with roles in human disease; however, its mRNA targets in the brain are largely unknown. Using cross-linking immunoprecipitation and sequencing (CLIP-seq), we define its CNS targets, which are enriched for 3' UTRs in synaptic protein-coding genes. Using a massively parallel reporter assay framework, we test the consequence of CELF6 expression on target sequences, with and without mutating putative binding motifs. Where CELF6 exerts an effect on sequences, it is largely to decrease RNA abundance, which is reversed by mutating UGU-rich motifs. This is also the case for CELF3-5, with a protein-dependent effect on magnitude. Finally, we demonstrate that targets are derepressed in CELF6-mutant mice, and at least two key CNS proteins, FOS and FGF13, show altered protein expression levels and localization. Our works find, in addition to previously identified roles in splicing, that CELF6 is associated with repression of its CNS targets via the 3' UTR in vivo.

CLIP and Massively Parallel Functional Analysis of CELF6 Reveal a Role in Destabilizing Synaptic Gene mRNAs through Interaction with 3' UTR Elements

CELF6 is a CELF-RNA-binding protein, and thus part of a protein family with roles in human disease; however, its mRNA targets in the brain are largely unknown. Using cross-linking immunoprecipitation and sequencing (CLIP-seq), we define its CNS targets, which are enriched for 3′ UTRs in synaptic protein-coding genes. Using a massively parallel reporter assay framework, we test the consequence of CELF6 expression on target sequences, with and without mutating putative binding motifs. Where CELF6 exerts an effect on sequences, it is largely to decrease RNA abundance, which is reversed by mutating UGU-rich motifs. This is also the case for CELF3–5, with a protein-dependent effect on magnitude. Finally, we demonstrate that targets are derepressed in CELF6-mutant mice, and at least two key CNS proteins, FOS and FGF13, show altered protein expression levels and localization. Our works find, in addition to previously identified roles in splicing, that CELF6 is associated with repression of its CNS targets via the 3′ UTR in vivo.

Rieger et al. assay the function of the RNA-binding protein CELF6 by defining its targets in the brain. They show that CELF6 largely binds 3′ UTRs of synaptic mRNAs. Using a massively parallel reporter assay, they further show that CELF6 and other CELFs are associated with lower mRNA abundance and that targets are derepressed in Celf6-knockout mice in vivo.

Genome-wide scan for selection signatures reveals novel insights into the adaptive capacity in local North African cattle

Natural-driven selection is supposed to have left detectable signatures on the genome of North African cattle which are often characterized by the fixation of genetic variants associated with traits under selection pressure and/or an outstanding genetic differentiation with other populations at particular loci. Here, we investigate the population genetic structure and we provide a first outline of potential selection signatures in North African cattle using single nucleotide polymorphism genotyping data. After comparing our data to African, European and indicine cattle populations, we identified 36 genomic regions using three extended haplotype homozygosity statistics and 92 outlier markers based on Bayescan test. The 13 outlier windows detected by at least two approaches, harboured genes (e.g. GH1, ACE, ASIC3, HSPH1, MVD, BCL2, HIGD2A, CBFA2T3) that may be involved in physiological adaptations required to cope with environmental stressors that are typical of the North African area such as infectious diseases, extended drought periods, scarce food supply, oxygen scarcity in the mountainous areas and high-intensity solar radiation. Our data also point to candidate genes involved in transcriptional regulation suggesting that regulatory elements had also a prominent role in North African cattle response to environmental constraints. Our study yields novel insights into the unique adaptive capacity in these endangered populations emphasizing the need for the use of whole genome sequence data to gain a better understanding of the underlying molecular mechanisms.

Translational Control during Mammalian Neocortex Development and Postembryonic Neuronal Function

The control of mRNA translation has key roles in the regulation of gene expression and biological processes such as mammalian cellular differentiation and identity. Methodological advances in the last decade have resulted in considerable progress towards understanding how translational control contributes to the regulation of diverse biological phenomena. In this review, we discuss recent findings in the involvement of translational control in the mammalian neocortex development and neuronal biology. We focus on regulatory mechanisms that modulate translational efficiency during neural stem cells self-renewal and differentiation, as well as in neuronal-related processes such as synapse, plasticity, and memory.

RNA-binding proteins balance brain function in health and disease

Posttranscriptional gene expression including splicing, RNA transport, translation, and RNA decay provides an important regulatory layer in many if not all molecular pathways. Research in the last decades has positioned RNA-binding proteins (RBPs) right in the center of posttranscriptional gene regulation. Here, we propose interdependent networks of RBPs to regulate complex pathways within the central nervous system (CNS). These are involved in multiple aspects of neuronal development and functioning, including higher cognition. Therefore, it is not sufficient to unravel the individual contribution of a single RBP and its consequences but rather to study and understand the tight interplay between different RBPs. In this review, we summarize recent findings in the field of RBP biology and discuss the complex interplay between different RBPs. Second, we emphasize the underlying dynamics within an RBP network and how this might regulate key processes such as neurogenesis, synaptic transmission, and synaptic plasticity. Importantly, we envision that dysfunction of specific RBPs could lead to perturbation within the RBP network. This would have direct and indirect (compensatory) effects in mRNA binding and translational control leading to global changes in cellular expression programs in general and in synaptic plasticity in particular. Therefore, we focus on RBP dysfunction and how this might cause neuropsychiatric and neurodegenerative disorders. Based on recent findings, we propose that alterations in the entire regulatory RBP network might account for phenotypic dysfunctions observed in complex diseases including neurodegeneration, epilepsy, and autism spectrum disorders.

Systemic neuro-dysregulation in depression: Evidence from genome-wide association

Depression is the world's leading cause of disability. Greater understanding of the neurobiological basis of depression is necessary for developing novel treatments with improved efficacy and acceptance. Recently, major advances have been made in the search for genetic variants associated with depression which may help to elucidate etiological mechanisms. The present review has two major objectives. First, we offer a brief review of two major biological systems with strong evidence for involvement in depression pathology: neurotransmitter systems and the stress response. Secondly, we provide a synthesis of the functions of the 269 genes implicated by the most recent genome-wide meta-analysis, supporting the importance of these systems in depression and providing insights into other possible mechanisms involving neurodevelopment, neurogenesis, and neurodegeneration. Our goal is to undertake a broad, preliminary stock-taking of the most recent hypothesis-free findings and examine the weight of the evidence supporting these existing theories and highlighting novel directions. This qualitative review and accompanying gene function table provides a valuable resource and guide for basic and translational researchers, with suggestions for future mechanistic research, leveraging genetics to prioritize studies on the neurobiological processes involved in depression etiology and treatment.

Minor DNA methylation changes are observed in spermatozoa prepared using different protocols

BACKGROUND

DNA methylation patterns can show transgenerational inheritance and are influenced by lifestyle and environmental factors. It is suggested that these patterns can be changed by assisted reproductive technology.

OBJECTIVES

To evaluate the impact of two different sperm preparation methods, conventional density gradient centrifugation (DGC) vs. density gradient centrifugation followed by magnetic-activated cell sorting (MACS) of non-apoptotic spermatozoa, on sperm DNA methylation profile.

MATERIALS AND METHODS

We analyzed semen of patients included in our IVF treatment program. Half of the semen from each included patient was prepared for ICSI using the DGC method and the other half with DGC followed by MACS. The remaining samples were processed for DNA methylation analysis with reduced representation bisulfite sequencing (RRBS). In addition to the DNA methylation profile, we assessed the morphology and DNA fragmentation of spermatozoa.

RESULTS

RRBS analysis revealed that the average genome-wide methylation level was similar between both groups (DGC vs. MACS group) and ranged from 0.53 to 0.56. Furthermore, RRBS analysis identified 99 differentially methylated regions (DMRs) and 800 differentially methylated positions (DMPs). In the DGC group, 43 DMRs and 392 DMPs were hypermethylated whereas 56 DMRs and 408 DMPs were hypomethylated compared with those in the MACS group. When DMRs and DMPs were annotated to genes, 3 genes associated with imprinting were found: IGF2, PRDM16, and CLF4/BRUNOL4. The percentage of morphologically normal spermatozoa (MACS vs. DGC; 14.0 ± 10.8 vs. 13.2 ± 10.0; P = .335) and of spermatozoa with fragmented DNA of patients with RRBS analysis (22.9 ± 21.1% vs. 34.4 ± 21.2; P = .529) were also similar between groups.

DISCUSSION AND CONCLUSION

Although the average genome-wide level of sperm DNA methylation was similar in both sample groups, a distinctive number of methylation changes were observed in DMR and DMP levels. A larger number of samples should be analyzed and additional sperm preparation methods should be tested to confirm our findings.

Large-scale analysis of the position-dependent binding and regulation of human RNA binding proteins

Background

RNA binding proteins (RBPs) play essential roles in the regulation of RNA metabolism. Recent studies have disclosed that RBPs achieve their functions via binding to their targets in a position-dependent pattern on RNAs. However, few studies have systematically addressed the associations between the RBP's functions and their positional binding preferences.

Methods

Here, we present large-scale analyses on the functional targets of human RBPs by integrating the enhanced cross-linking and immunoprecipitation followed by sequencing (eCLIP-seq) datasets and the shRNA knockdown followed by RNA-seq datasets that are deposited in the integrated ENCyclopedia of DNA Elements in the human genome (ENCODE) data portal.

Results

We found that (1) binding to the translation termination site and the 3'untranslated region is important to most human RBPs in the RNA decay regulation; (2) RBPs' binding and regulation follow a cell-type specific pattern.

Conclusions

These analysis results show the strong relationship between the binding position and the functions of RBPs, which provides novel insights into the RBPs' regulation mechanisms.

Single-cell DNA methylation sequencing reveals epigenetic alterations in mouse oocytes superovulated with different dosages of gonadotropins

BACKGROUND

Epigenetic abnormalities caused by superovulation have recently attracted increasing attention. Superovulation with exogenous hormones may prevent oocytes from establishing an appropriate epigenetic state, and this effect may extend to the methylation programming in preimplantation embryos, as de novo DNA methylation is a function of developmental stage of follicles and oocyte size. Follicle-stimulating hormone (FSH) and human menopausal gonadotropin (hMG) are common gonadotropins used for superovulation, and appropriate concentrations of these gonadotropins might be necessary. However, no systematic study on the effects of DNA methylation alterations in oocytes associated with superovulation with different dosages of FSH/hMG at the single-cell level has yet been reported. In the current study, different dosages of FSH/hMG combined with human chorionic gonadotropin (hCG) were used in female mice to generate experimental groups, while naturally matured oocytes and oocytes superovulated with only hCG were respectively used as controls. Single-cell level DNA methylation sequencing was carried out on all these matured oocytes.

RESULTS

In this study, we revealed that the genome-wide methylation pattern and CG methylation level of the maternal imprinting control regions of all mature oocytes were globally conserved and stable. However, methylation alterations associated with superovulation were found at a specific set of loci, and the differentially methylated regions (DMRs) mainly occurred in regions other than promoters. Furthermore, some of the annotated genes in the DMRs were involved in biological processes such as glucose metabolism, nervous system development, cell cycle, cell proliferation, and embryo implantation and were altered in all dosages of FSH/hMG group (for example, Gfod2 and SYF2). Other genes were impaired only after high gonadotropin dosages (for instance, Sox17 and Phactr4).

CONCLUSIONS

In conclusion, the current study addressed the effects of superovulation on DNA methylation from the perspective of different dosages of gonadotropins at the single-cell level. We found that the genome-wide DNA methylation landscape was globally preserved irrespective of superovulation or of the kind and dosage of gonadotropins used, whereas the methylation alterations associated with superovulation occurred at a specific set of loci. These observed effects reflect that superovulation recruits oocytes that would not normally be ovulated or that have not undergone complete epigenetic maturation. Our results provide an important reference for the safety assessment of superovulation with different dosages of gonadotropins. However, it should be noted that this study has some limitations, as the sample number and library coverage of analyzed oocytes were relatively low. Future studies with larger sample sizes and high-coverage libraries that examine the effects of superovulation on embryo development and offspring health as well as the underlying mechanisms are still needed.

Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Common variants identified by large-scale genomewide association studies cannot account fully account for the heritability of schizophrenia (SCZ). Here, the authors report high-coverage whole-genome sequencing of 1162 SCZ cases and 936 controls and explore the contribution of different types of variants to SCZ.

Translational derepression of Elavl4 isoforms at their alternative 5' UTRs determines neuronal development

Neurodevelopment requires precise regulation of gene expression, including post-transcriptional regulatory events such as alternative splicing and mRNA translation. However, translational regulation of specific isoforms during neurodevelopment and the mechanisms behind it remain unknown. Using RNA-seq analysis of mouse neocortical polysomes, here we report translationally repressed and derepressed mRNA isoforms during neocortical neurogenesis whose orthologs include risk genes for neurodevelopmental disorders. We demonstrate that the translation of distinct mRNA isoforms of the RNA binding protein (RBP), Elavl4, in radial glia progenitors and early neurons depends on its alternative 5' UTRs. Furthermore, 5' UTR-driven Elavl4 isoform-specific translation depends on upstream control by another RBP, Celf1. Celf1 regulation of Elavl4 translation dictates development of glutamatergic neurons. Our findings reveal a dynamic interplay between distinct RBPs and alternative 5' UTRs in neuronal development and underscore the risk of post-transcriptional dysregulation in co-occurring neurodevelopmental disorders.

dSreg: a Bayesian model to integrate changes in splicing and RNA-binding protein activity

MOTIVATION

Alternative splicing (AS) is an important mechanism in the generation of transcript diversity across mammals. AS patterns are dynamically regulated during development and in response to environmental changes. Defects or perturbations in its regulation may lead to cancer or neurological disorders, among other pathological conditions. The regulatory mechanisms controlling AS in a given biological context are typically inferred using a two-step framework: differential AS analysis followed by enrichment methods. These strategies require setting rather arbitrary thresholds and are prone to error propagation along the analysis.

RESULTS

To overcome these limitations, we propose dSreg, a Bayesian model that integrates RNA-seq with data from regulatory features, e.g. binding sites of RNA-binding proteins. dSreg identifies the key underlying regulators controlling AS changes and quantifies their activity while simultaneously estimating the changes in exon inclusion rates. dSreg increased both the sensitivity and the specificity of the identified AS changes in simulated data, even at low read coverage. dSreg also showed improved performance when analyzing a collection of knock-down RNA-binding proteins' experiments from ENCODE, as opposed to traditional enrichment methods, such as over-representation analysis and gene set enrichment analysis. dSreg opens the possibility to integrate a large amount of readily available RNA-seq datasets at low coverage for AS analysis and allows more cost-effective RNA-seq experiments.

AVAILABILITY AND IMPLEMENTATION

dSreg was implemented in python using stan and is freely available to the community at https://bitbucket.org/cmartiga/dsreg.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A novel pathway regulates social hierarchy via lncRNA AtLAS and postsynaptic synapsin IIb

Dominance hierarchy is a fundamental phenomenon in grouped animals and human beings, however, the underlying regulatory mechanisms remain elusive. Here, we report that an antisense long non-coding RNA (lncRNA) of synapsin II, named as AtLAS, plays a crucial role in the regulation of social hierarchy. AtLAS is decreased in the prefrontal cortical excitatory pyramidal neurons of dominant mice; consistently, silencing or overexpression of AtLAS increases or decreases the social rank, respectively. Mechanistically, we show that AtLAS regulates alternative polyadenylation of synapsin II gene and increases synapsin 2b (syn2b) expression. Syn2b reduces AMPA receptor (AMPAR)-mediated excitatory synaptic transmission through a direct binding with AMPAR at the postsynaptic site via its unique C-terminal sequence. Moreover, a peptide disrupting the binding of syn2b with AMPARs enhances the synaptic strength and social ranks. These findings reveal a novel role for lncRNA AtLAS and its target syn2b in the regulation of social behaviors by controlling postsynaptic AMPAR trafficking.

Comprehensive transcriptome analysis of cochlear spiral ganglion neurons at multiple ages

Inner ear cochlear spiral ganglion neurons (SGNs) transmit sound information to the brainstem. Recent single cell RNA-Seq studies have revealed heterogeneities within SGNs. Nonetheless, much remains unknown about the transcriptome of SGNs, especially which genes are specifically expressed in SGNs. To address these questions, we needed a deeper and broader gene coverage than that in previous studies. We performed bulk RNA-Seq on mouse SGNs at five ages, and on two reference cell types (hair cells and glia). Their transcriptome comparison identified genes previously unknown to be specifically expressed in SGNs. To validate our dataset and provide useful genetic tools for this research field, we generated two knockin mouse strains: Scrt2-P2A-tdTomato and Celf4-3xHA-P2A-iCreER-T2A-EGFP. Our comprehensive analysis confirmed the SGN-selective expression of the candidate genes, testifying to the quality of our transcriptome data. These two mouse strains can be used to temporally label SGNs or to sort them.

Convergent Evidence That ZNF804A Is a Regulator of Pre-messenger RNA Processing and Gene Expression

Genome-wide association studies have linked common variation in ZNF804A with an increased risk of schizophrenia. However, little is known about the biology of ZNF804A and its role in schizophrenia. Here, we investigate the function of ZNF804A using a variety of complementary molecular techniques. We show that ZNF804A is a nuclear protein that interacts with neuronal RNA splicing factors and RNA-binding proteins including RBFOX1, which is also associated with schizophrenia, CELF3/4, components of the ubiquitin-proteasome system and the ZNF804A paralog, GPATCH8. GPATCH8 also interacts with splicing factors and is localized to nuclear speckles indicative of a role in pre-messenger RNA (mRNA) processing. Sequence analysis showed that GPATCH8 contains ultraconserved, alternatively spliced poison exons that are also regulated by RBFOX proteins. ZNF804A knockdown in SH-SY5Y cells resulted in robust changes in gene expression and pre-mRNA splicing converging on pathways associated with nervous system development, synaptic contact, and cell adhesion. We observed enrichment (P = 1.66 × 10-9) for differentially spliced genes in ZNF804A-depleted cells among genes that contain RBFOX-dependent alternatively spliced exons. Differentially spliced genes in ZNF804A-depleted cells were also enriched for genes harboring de novo loss of function mutations in autism spectrum disorder (P = 6.25 × 10-7, enrichment 2.16) and common variant alleles associated with schizophrenia (P = .014), bipolar disorder and schizophrenia (P = .003), and autism spectrum disorder (P = .005). These data suggest that ZNF804A and its paralogs may interact with neuronal-splicing factors and RNA-binding proteins to regulate the expression of a subset of synaptic and neurodevelopmental genes.

Getting to the Cores of Autism

The genetic architecture of autism spectrum disorder (ASD) is itself a diverse allelic spectrum that consists of rare de novo or inherited variants in hundreds of genes and common polygenic risk at thousands of loci. ASD susceptibility genes are interconnected at the level of transcriptional and protein networks, and many function as genetic regulators of neurodevelopment or synaptic proteins that regulate neural activity. So that the core underlying neuropathologies can be further elucidated, we emphasize the importance of first defining subtypes of ASD on the basis of the phenotypic signatures of genes in model systems and humans.

PRAS: Predicting functional targets of RNA binding proteins based on CLIP-seq peaks

RNA-protein interaction plays important roles in post-transcriptional regulation. Recent advancements in cross-linking and immunoprecipitation followed by sequencing (CLIP-seq) technologies make it possible to detect the binding peaks of a given RNA binding protein (RBP) at transcriptome scale. However, it is still challenging to predict the functional consequences of RBP binding peaks. In this study, we propose the Protein-RNA Association Strength (PRAS), which integrates the intensities and positions of the binding peaks of RBPs for functional mRNA targets prediction. We illustrate the superiority of PRAS over existing approaches on predicting the functional targets of two related but divergent CELF (CUGBP, ELAV-like factor) RBPs in mouse brain and muscle. We also demonstrate the potential of PRAS for wide adoption by applying it to the enhanced CLIP-seq (eCLIP) datasets of 37 RNA decay related RBPs in two human cell lines. PRAS can be utilized to investigate any RBPs with available CLIP-seq peaks. PRAS is freely available at http://ouyanglab.jax.org/pras/.

Posttranscriptional Gene Regulation of the GABA Receptor to Control Neuronal Inhibition

Behavior and higher cognition rely on the transfer of information between neurons through specialized contact sites termed synapses. Plasticity of neuronal circuits, a prerequisite to respond to environmental changes, is intrinsically coupled with the nerve cell’s ability to form, structurally modulate or remove synapses. Consequently, the synaptic proteome undergoes dynamic alteration on demand in a spatiotemporally restricted manner. Therefore, proper protein localization at synapses is essential for synaptic function. This process is regulated by: (i) protein transport and recruitment; (ii) local protein synthesis; and (iii) synaptic protein degradation. These processes shape the transmission efficiency of excitatory synapses. Whether and how these processes influence synaptic inhibition is, however, widely unknown. Here, we summarize findings on fundamental regulatory processes that can be extrapolated to inhibitory synapses. In particular, we focus on known aspects of posttranscriptional regulation and protein dynamics of the GABA receptor (GABAR). Finally, we propose that local (co)-translational control mechanism might control transmission of inhibitory synapses.

Human genes influence the interaction between Streptococcus mutans and host caries susceptibility: a genome-wide association study in children with primary dentition

Streptococcus mutans is a well-known cause of dental caries, due to its acidogenicity, aciduricity, and ability to synthesize exopolysaccharides in dental plaques. Intriguingly, not all children who carry S. mutans manifest caries, even with similar characteristics in oral hygiene, diet, and other environmental factors. This phenomenon suggests that host susceptibility potentially plays a role in the development of dental caries; however, the association between host genetics, S. mutans, and dental caries remains unclear. Therefore, this study examined the influence of host gene-by-S. mutans interaction on dental caries. Genome-wide association analyses were conducted in 709 US children (<13 years old), using the dbGap database acquired from the center for oral health research in appalachia (COHRA) and the Iowa Head Start programmes (GEIRS). A generalized estimating equation was used to examine the gene-by-S. mutans interaction effects on the outcomes (decayed and missing/filled primary teeth due to caries). Sequentially, the COHRA and GEIRS data were used to identify potential interactions and replicate the findings. Three loci at the genes interleukin 32 (IL32), galactokinase 2 (GALK2), and CUGBP, Elav-like family member 4 (CELF4) were linked to S. mutans carriage, and there was a severity of caries at a suggestive significance level among COHRA children (P < 9 × 10^-5), and at a nominal significance level among GEIRS children (P = 0.047-0.001). The genetic risk score that combined the three loci also significantly interacted with S. mutans (P < 0.000 1). Functional analyses indicated that the identified genes are involved in the host immune response, galactose carbohydrate metabolism, and food-rewarding system, which could potentially be used to identify children at high risk for caries and to develop personalized caries prevention strategies.

Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions

Major depression is a debilitating psychiatric illness that is typically associated with low mood and anhedonia. Depression has a heritable component that has remained difficult to elucidate with current sample sizes due to the polygenic nature of the disorder. To maximize sample size, we meta-analyzed data on 807,553 individuals (246,363 cases and 561,190 controls) from the three largest genome-wide association studies of depression. We identified 102 independent variants, 269 genes, and 15 genesets associated with depression, including both genes and gene pathways associated with synaptic structure and neurotransmission. An enrichment analysis provided further evidence of the importance of prefrontal brain regions. In an independent replication sample of 1,306,354 individuals (414,055 cases and 892,299 controls), 87 of the 102 associated variants were significant after multiple testing correction. These findings advance our understanding of the complex genetic architecture of depression and provide several future avenues for understanding etiology and developing new treatment approaches.

Biological annotation of genetic loci associated with intelligence in a meta-analysis of 87,740 individuals

Variance in IQ is associated with a wide range of health outcomes, and 1% of the population are affected by intellectual disability. Despite a century of research, the fundamental neural underpinnings of intelligence remain unclear. We integrate results from genome-wide association studies (GWAS) of intelligence with brain tissue and single cell gene expression data to identify tissues and cell types associated with intelligence. GWAS data for IQ (N = 78,308) were meta-analyzed with a study comparing 1247 individuals with mean IQ ~170 to 8185 controls. Genes associated with intelligence implicate pyramidal neurons of the somatosensory cortex and CA1 region of the hippocampus, and midbrain embryonic GABAergic neurons. Tissue-specific analyses find the most significant enrichment for frontal cortex brain expressed genes. These results suggest specific neuronal cell types and genes may be involved in intelligence and provide new hypotheses for neuroscience experiments using model systems.

The TIA1 RNA-Binding Protein Family Regulates EIF2AK2-Mediated Stress Response and Cell Cycle Progression

TIA1 and TIAL1 encode a family of U-rich element mRNA-binding proteins ubiquitously expressed and conserved in metazoans. Using PAR-CLIP, we determined that both proteins bind target sites with identical specificity in 3' UTRs and introns proximal to 5' as well as 3' splice sites. Double knockout (DKO) of TIA1 and TIAL1 increased target mRNA abundance proportional to the number of binding sites and also caused accumulation of aberrantly spliced mRNAs, most of which are subject to nonsense-mediated decay. Loss of PRKRA by mis-splicing triggered the activation of the double-stranded RNA (dsRNA)-activated protein kinase EIF2AK2/PKR and stress granule formation. Ectopic expression of PRKRA cDNA or knockout of EIF2AK2 in DKO cells rescued this phenotype. Perturbation of maturation and/or stability of additional targets further compromised cell cycle progression. Our study reveals the essential contributions of the TIA1 protein family to the fidelity of mRNA maturation, translation, and RNA-stress-sensing pathways in human cells.

Polygenic risk score for schizophrenia is not strongly associated with the expression of specific genes or gene sets

BACKGROUND

The polygenic risk score (PRS) is derived from single nucleotide polymorphisms (SNPs) including those that are genome-wide significant and also including a large number of others more weakly associated with schizophrenia. Such variants are widely dispersed, though concentrated near genes expressed in the brain, and it has been proposed that these SNP associations result from impacts on cell regulatory networks that ultimately affect the expression or function of a modest number of 'core' genes. A previous study showed association of some genome-wide association study-significant variants with expression of a number of genes, by examining pairwise correlations of gene expression with SNP genotypes.

METHODS

The present study used data downloaded from the CommonMind Consortium site, consisting of SNP genotypes and RNAseq expression data from the dorsolateral prefrontal cortex, to examine whether the expression of individual genes or sets of genes correlated with PRS in 207 controls and 209 schizophrenia cases.

RESULTS

Although the PRS was significantly associated with phenotype, the correlations with genes and gene sets followed distributions expected by chance. Thus, this analysis failed to show that the PRS captures a cumulative effect of multiple variants impacting the expression of a small number of genes and it failed to focus attention on a small number of genes of biological relevance.

CONCLUSION

The multiple SNP associations observed in schizophrenia may result from other mechanisms, including effects mediated indirectly through environmental risk factors.

Dysregulated Translation in Neurodevelopmental Disorders: An Overview of Autism-Risk Genes Involved in Translation

Regulated local translation-whereby specific mRNAs are transported and localized in subcellular domains where they are translated in response to regional signals-allows for remote control of gene expression to concentrate proteins in subcellular compartments. Neurons are highly polarized cells with unique features favoring local control for axonal pathfinding and synaptic plasticity, which are key processes involved in constructing functional circuits in the developing brain. Neurodevelopmental disorders are caused by genetic or environmental factors that disturb the nervous system's development during prenatal and early childhood periods. The growing list of genetic mutations that affect mRNA translation raises the question of whether aberrant translatomes in individuals with neurodevelopmental disorders share common molecular features underlying their stereotypical phenotypes and, vice versa, cause a certain degree of phenotypic heterogeneity. Here, we briefly give an overview of the role of local translation during neuronal development. We take the autism-risk gene list and discuss the molecules that (perhaps) are involved in mRNA transport and translation. Both exaggerated and suppressed translation caused by mutations in those genes have been identified or suggested. Finally, we discuss some proof-of-principle regimens for use in autism mouse models to correct dysregulated translation.

A weighted burden test using logistic regression for integrated analysis of sequence variants, copy number variants and polygenic risk score

Previously described methods of analysis allow variants in a gene to be weighted more highly according to rarity and/or predicted function and then for the variant contributions to be summed into a gene-wise risk score, which can be compared between cases and controls using a t-test. However, this does not allow incorporating covariates into the analysis. Schizophrenia is an example of an illness where there is evidence that different kinds of genetic variation can contribute to risk, including common variants contributing to a polygenic risk score (PRS), very rare copy number variants (CNVs) and sequence variants. A logistic regression approach has been implemented to compare the gene-wise risk scores between cases and controls, while incorporating as covariates population principal components, the PRS and the presence of pathogenic CNVs and sequence variants. A likelihood ratio test is performed, comparing the likelihoods of logistic regression models with and without this score. The method was applied to an ethnically heterogeneous exome-sequenced sample of 6000 controls and 5000 schizophrenia cases. In the raw analysis, the test statistic is inflated but inclusion of principal components satisfactorily controls for this. In this dataset, the inclusion of the PRS and effect from CNVs and sequence variants had only small effects. The set of genes which are FMRP targets showed some evidence for enrichment of rare, functional variants among cases (p = 0.0005). This approach can be applied to any disease in which different kinds of genetic and non-genetic risk factors make contributions to risk.