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

Modulation of netrin-1/DCC signaling pathway by Jiawei Kongsheng Zhenzhong Pill improves synaptic structural plasticity in PSD rats

Jiawei Kongsheng Zhenzhong Pill(JKZP) is based on Kongsheng Zhenzhong Pill contained in the Tang Dynasty's "Thousand Golden Prescriptions," which exhibited good anti-ischemic and antidepressant effects in the previous study. However, its specific effects on post-stroke depression (PSD) and the mechanism are not clear. This study aimed to investigate the effects of JKZP in the treatment of PSD and related mechanisms. The decoction of JKZP was first analyzed for its medicinal chemical composition and screened for representative components of JKZP. The Middle cerebral artery occlusion (MCAO) method combined with solitary rearing and chronic unpredictable mild stress (CUMS) was used to establish a rat model of PSD, and to observe the effects of JKZP on the behavior and synaptic plasticity of PSD rats, and to investigate the mechanism of JKZP in the treatment of PSD by detecting the mRNA level, protein expression and activity of Netrin-1/DCC signaling pathway-related proteins. The results showed that the JKZP decoction contained loganin, β-asarone and other pharmaceutical ingredients, which have been reported to protect against cerebral ischemic injury and antidepressant effects. JKZP significantly improved the depression-like behavior of PSD rats and improved the damage to pyramidal neurons in the medial prefrontal cortex (mPFC) of PSD rats. Moreover, JKZP increased the density of dendritic spines in the mPFC of PSD rats, improved synaptic gap width and thickness of the post-synaptic density, and increased the number of synaptic vesicles. The results of Real-Time quantitative reverse transcription PCR (qRT-PCR), Western blotting, and pull-down assays revealed that JKZP increased netrin-1, deleted in colorectal cancer (DCC), and focal adhesion kinase (FAK) mRNA and protein expression, elevated the p-FAK/FAK ratio, and decreased myosin II protein expression and Ras homolog gene family member A (RhoA-GTP) activity in the mPFC of PSD rats. Taken together, JKZP can affect synaptic structural remodeling and improve depressive manifestations and neuronal damage in PSD rats by regulating the expression and activity of signaling molecules related to the netrin-1/DCC signaling pathway.

Heterozygous CELF4 variants in the N-term region crucial for the RNA-binding activity lead to neurodevelopmental disorder and obesity

RNA-binding proteins play a key role in post-transcriptional events, such as mRNA splicing, transport, stability, translation and decay. Dysregulation of RNA life can have dramatic consequences. CELF RNA-binding proteins appear to be essential during embryo development. In this study, we identified 15 patients with heterozygous missense or loss-of-function variants in the CELF4 gene by exome or genome sequencing. All variants affecting the N-terminus of the protein are essential and sufficient for the RNA-binding and splicing activity or RRM domains. Most patients presented with neurodevelopmental disorders including global developmental delay/intellectual disability (11/14), seizures (9/15) and overweight/obesity (10/14) that began in childhood. Clinical features are similar to the reported celf4-mouse mutant phenotype. This study highlights the essential role of CELF4 in development and its involvement as a novel etiology of neurodevelopmental disorders with obesity.

Neonatal Encephalopathy: Novel Phenotypes and Genotypes Identified by Genome Sequencing

Background and Objectives

Neonatal encephalopathy (NE) is characterized by an abnormal level of consciousness with or without seizures in the neonatal period. It affects 1-6/1,000 live term newborns. We applied genome sequencing (GS) in term newborns with NE to investigate the underlying genetic causes.

Methods

We enrolled term newborns according to inclusion/exclusion criteria during their Neonatal Intensive Care admission. We performed GS trio and applied bioinformatic tools. We developed pipelines for manual filters. We applied in silico prediction tools, protein 3D modeling, and functional characterization to assess the pathogenicity of variants.

Results

Seventeen newborns fulfilled inclusion criteria. We identified 12 variants in 10 genes. We classified 4 variants in PPP2R5D, BCOR, CFL2, and SCN2A (previously established disease genes) as pathogenic/likely pathogenic; 7 variants in DST (previously established disease gene), STAB2, CELF4, SORCS2, CTNND2, and ASTN1 (5 candidate genes) as variants of uncertain significance (VUS); and one variant in STAB2 as likely benign. The CELF4 and ASTN1 copy number variants (CNVs) resulted in structural changes in protein 3D models. The functional characterization of SORCS2 VUS revealed disruption of SorCS2 dimer formation and confirmed its pathogenicity. The functional characterization of STAB2 variants updated their characterization from VUS/likely benign to benign. The CTNND2 VUS resulted in a shift in 3D protein structure. We were not able to perform protein 3D modeling and functional characterization of two DST VUS. We are not certain whether CTNND2 and DST variants may be causative of NE in our study.

Discussion

The diagnostic rate of research GS was 41% in our prospective study. We broaden the phenotypic spectrum of PPP2R5D-associated Hogue-Janssens syndrome 1, CFL2-associated nemaline myopathy 7, and BCOR-associated oculo-facio-cardio-dental syndrome to include NE and/or neonatal seizures. We identified 3 candidate genes (SORCS2, CELF4, ASTN1) that may cause NE. We believe that protein 3D modeling is an important tool to assess the pathogenicity of CNVs and may advance the discoveries of novel genetic diseases. However, functional characterization of missense variants is essential for discoveries of novel genetic diseases. It seems that GS can help identify more candidate genes compared with ES.

GAN-WGCNA: Calculating gene modules to identify key intermediate regulators in cocaine addiction

Understanding time-series interplay of genes is essential for diagnosis and treatment of disease. Spatio-temporally enriched NGS data contain important underlying regulatory mechanisms of biological processes. Generative adversarial networks (GANs) have been used to augment biological data to describe hidden intermediate time-series gene expression profiles during specific biological processes. Developing a pipeline that uses augmented time-series gene expression profiles is needed to provide an unbiased systemic-level map of biological processes and test for the statistical significance of the generated dataset, leading to the discovery of hidden intermediate regulators. Two analytical methods, GAN-WGCNA (weighted gene co-expression network analysis) and rDEG (rescued differentially expressed gene), interpreted spatiotemporal information and screened intermediate genes during cocaine addiction. GAN-WGCNA enables correlation calculations between phenotype and gene expression profiles and visualizes time-series gene module interplay. We analyzed a transcriptome dataset of two weeks of cocaine self-administration in C57BL/6J mice. Utilizing GAN-WGCNA, two genes (Alcam and Celf4) were selected as missed intermediate significant genes that showed high correlation with addiction behavior. Their correlation with addictive behavior was observed to be notably significant in aspect of statistics, and their expression and co-regulation were comprehensively mapped in terms of time, brain region, and biological process.

A genome-wide association study of occupational creativity and its relations with well-being and career success

Creativity is one defining characteristic of human species. There have been mixed findings on how creativity relates to well-being, and little is known about its relationship with career success. We conduct a large-scale genome-wide association study to examine the genetic architecture of occupational creativity, and its genetic correlations with well-being and career success. The SNP-h2 estimates range from 0.08 (for managerial creativity) to 0.22 (for artistic creativity). We record positive genetic correlations between occupational creativity with autism, and positive traits and well-being variables (e.g., physical height, and low levels of neuroticism, BMI, and non-cancer illness). While creativity share positive genetic overlaps with indicators of high career success (i.e., income, occupational status, and job satisfaction), it also has a positive genetic correlation with age at first birth and a negative genetic correlation with number of children, indicating creativity-related genes may reduce reproductive success.

Identifying regulatory elements and their RNA-binding proteins in the 3' untranslated regions of papillomavirus late mRNAs

Human papillomaviruses (HPVs) infect cutaneous and mucosal epithelia to cause benign (warts) and malignant lesions (e.g. cervical cancer). Bovine papillomaviruses (BPVs) infect fibroblasts to cause fibropapillomas but can also infect cutaneous epithelial cells. For HPV-1, -16, -31 and BPV-1, cis-acting RNA elements in the late 3' untranslated region (3'UTR) control expression of virus proteins by binding host cell proteins. The present study compared the effects on gene expression of the cis-acting elements of seven PV late 3'UTRs (HPV-6b, -11, -16, -31 and BPV-1, -3 and -4) representing a range of different genera and species and pathological properties. pSV-beta-galactosidase reporter plasmids containing the late 3'UTRs from seven PVs were transiently transfected into cervical adenocarcinoma HeLa cells, and reporter gene expression quantified by reverse transcription-quantitative PCR and a beta-galactosidase assay. All elements inhibited gene expression in keratinocytes. Cancer-related types HPV-16 and -31, had the greatest inhibitory activity whereas the lowest inhibition was found in the non-cancer related types, BPV-3 and HPV-11. Using RBPmap version 1.1, bioinformatics predictions of factors binding the elements identified proteins which function mainly in mRNA splicing. Markedly, in terms of protein binding motifs, BPV late 3'UTR elements were similar to those of HPV-1a but not to other HPVs. Using HPV-1a as a model and siRNA depletion, the bioinformatics predictions were tested and it was found that PABPC4 was responsible for some of the 3'UTR repressive activity. The data revealed candidate proteins that could control PV late gene expression.

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.

Social Determinants of Health in Cardio-Oncology: Multi-Level Strategies to Overcome Disparities in Care: JACC: CardioOncology State-of-the-Art Review

Addressing the need for more equitable cardio-oncology care requires attention to existing disparities in cardio-oncologic disease prevention and outcomes. This is particularly important among those affected by adverse social determinants of health (SDOH). The intricate relationship of SDOH, cancer diagnosis, and outcomes from cardiotoxicities associated with oncologic therapies is influenced by sociopolitical, economic, and cultural factors. Furthermore, mechanisms in cell signaling and epigenetic effects on gene expression link adverse SDOH to cancer and the CVD-related complications of oncologic therapies. To mitigate these disparities, a multifaceted strategy is needed that includes attention to health care access, policy, and community engagement for improved disease screening and management. Interdisciplinary teams must also promote cultural humility and competency and leverage new health technology to foster collaboration in addressing the impact of adverse SDOH in cardio-oncologic outcomes.

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.

RNA binding protein CUGBP1 mediates the liver metastasis of colorectal cancer by regulating the ErbB signal pathway

BACKGROUND

The CUGBP1 (CELF1) is differentially expressed in liver metastasis and no liver metastasis colorectal cancers (CRC) tissues and the function of CUGBP1 in CRC is still unclear.

METHODS

Five cases of colorectal adenocarcinoma and 6 cases of liver metastatic CRC lesions were collected and subjected to cDNA microarray and bioinformatical analyses. The quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to confirm the result. Cell function assays were used to study the function of CUGBP1, and the western blot was used to discover the change of the downstream molecules.

RESULTS

CUGBP1 was significantly elevated in liver metastatic CRC lesions. Besides, the CUGBP1 can promote proliferation, colony formation, invasion, metastasis abilities as well as increase the apoptosis rates of CRC cells. ERBB2 was positively related to the CUGBP1. Western blot results found that silence of CUGBP1 decreased the protein level of p-AKT and p-ERK without influence the expression level of total protein of AKT and ERK.

CONCLUSIONS

CUGBP1 can promote liver metastasis of CRC by promoting the phosphorylation of AKT and ERK through the ErbB signaling pathway. CUGBP1 is a potential biomarker for early detection of CRC and maybe a novel therapeutic target of CRC treatment, especially in liver metastasis.

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.

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/.