Dysregulation of Alternative Poly-adenylation as a Potential Player in Autism Spectrum Disorder

Autism-related KLHL17 and SYNPO act in concert to control activity-dependent dendritic spine enlargement and the spine apparatus

Dendritic spines, the tiny and actin-rich protrusions emerging from dendrites, are the subcellular locations of excitatory synapses in the mammalian brain that control synaptic activity and plasticity. Dendritic spines contain a specialized form of endoplasmic reticulum (ER), i.e., the spine apparatus, required for local calcium signaling and that is involved in regulating dendritic spine enlargement and synaptic plasticity. Many autism-linked genes have been shown to play critical roles in synaptic formation and plasticity. Among them, KLHL17 is known to control dendritic spine enlargement during development. As a brain-specific disease-associated gene, KLHL17 is expected to play a critical role in the brain, but it has not yet been well characterized. In this study, we report that KLHL17 expression in mice is strongly regulated by neuronal activity and KLHL17 modulates the synaptic distribution of synaptopodin (SYNPO), a marker of the spine apparatus. Both KLHL17 and SYNPO are F-actin-binding proteins linked to autism. SYNPO is known to maintain the structure of the spine apparatus in mature spines and contributes to synaptic plasticity. Our super-resolution imaging using expansion microscopy demonstrates that SYNPO is indeed embedded into the ER network of dendritic spines and that KLHL17 is closely adjacent to the ER/SYNPO complex. Using mouse genetic models, we further show that Klhl17 haploinsufficiency and knockout result in fewer dendritic spines containing ER clusters and an alteration of calcium events at dendritic spines. Accordingly, activity-dependent dendritic spine enlargement and neuronal activation (reflected by extracellular signal-regulated kinase (ERK) phosphorylation and C-FOS expression) are impaired. In addition, we show that the effect of disrupting the KLHL17 and SYNPO association is similar to the results of Klhl17 haploinsufficiency and knockout, further strengthening the evidence that KLHL17 and SYNPO act together to regulate synaptic plasticity. In conclusion, our findings unravel a role for KLHL17 in controlling synaptic plasticity via its regulation of SYNPO and synaptic ER clustering and imply that impaired synaptic plasticity contributes to the etiology of KLHL17-related disorders.

Deciphering the impact of genetic variation on human polyadenylation using APARENT2

BACKGROUND

3'-end processing by cleavage and polyadenylation is an important and finely tuned regulatory process during mRNA maturation. Numerous genetic variants are known to cause or contribute to human disorders by disrupting the cis-regulatory code of polyadenylation signals. Yet, due to the complexity of this code, variant interpretation remains challenging.

RESULTS

We introduce a residual neural network model, APARENT2, that can infer 3'-cleavage and polyadenylation from DNA sequence more accurately than any previous model. This model generalizes to the case of alternative polyadenylation (APA) for a variable number of polyadenylation signals. We demonstrate APARENT2's performance on several variant datasets, including functional reporter data and human 3' aQTLs from GTEx. We apply neural network interpretation methods to gain insights into disrupted or protective higher-order features of polyadenylation. We fine-tune APARENT2 on human tissue-resolved transcriptomic data to elucidate tissue-specific variant effects. By combining APARENT2 with models of mRNA stability, we extend aQTL effect size predictions to the entire 3' untranslated region. Finally, we perform in silico saturation mutagenesis of all human polyadenylation signals and compare the predicted effects of [Formula: see text] million variants against gnomAD. While loss-of-function variants were generally selected against, we also find specific clinical conditions linked to gain-of-function mutations. For example, we detect an association between gain-of-function mutations in the 3'-end and autism spectrum disorder. To experimentally validate APARENT2's predictions, we assayed clinically relevant variants in multiple cell lines, including microglia-derived cells.

CONCLUSIONS

A sequence-to-function model based on deep residual learning enables accurate functional interpretation of genetic variants in polyadenylation signals and, when coupled with large human variation databases, elucidates the link between functional 3'-end mutations and human health.

Selective poly adenylation predicts the efficacy of immunotherapy in patients with lung adenocarcinoma by multiple omics research

The aim of this study was to find the application value of selective polyadenylation in immune cell infiltration, biological transcription function and risk assessment of survival and prognosis in lung adenocarcinoma (LUAD). The processed original mRNA expression data of LUAD were downloaded, and the expression profiles of 594 patient samples were collected. The (APA) events in TCGA-NA-SEQ data were evaluated by polyadenylation site use Index (PDUI) values, and the invasion of stromal cells and immune cells and tumor purity were calculated to group and select the differential genes. Lasso regression and stratified analysis were used to examine the role of risk scores in predicting patient outcomes. The study also used the GDSC database to predict the chemotherapeutic sensitivity of each tumor sample and used a regression method to obtain an IC50 estimate for each specific chemotherapeutic drug treatment. Then CIBERSORT algorithm was used to conduct Spearman correlation analysis, immune regulatory factor analysis and TIDE immune system function analysis for gene expression level and immune cell content. Finally, the Kaplan-Meier curve was used to analyze the correlation between stromal score and the immune score of LUAD. In this study, APA's LUAD risk score prognostic model was constructed. KM survival analysis showed that immune score affected the prognosis of LUAD patients ( P  = 0.027) but the matrix score was not statistically significant ( P  = 0.1). We extracted 108 genes with APA events from 827 different genes and based on PUDI clustering and heat map, the survival rate of patients in the four groups was significantly different ( P  = 0.05). Multiple omics studies showed that risk score was significantly positively correlated with Macrophages M0, T cells Follicular helper, B cells naive and NK cells resting. It is significantly negatively correlated with dendritic cells resting, mast cells resting, monocyte, T cells CD4 memory resting and B cells memory. We further explored the relationship between the expression of immunosuppressor genes and risk score and found that ADORA2A, BTLA, CD160, CD244, CD274, CD96, CSF1R and CTLA4 genes were highly correlated with the risk score. Selective poly adenylation plays an important role in the development and progression of LUAD, immune invasion, tumor cell invasion and metastasis and biological transcription, and affects the survival and prognosis of LUAD patients.

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

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

RNA Modifications and RNA Metabolism in Neurological Disease Pathogenesis

The intrinsic cellular heterogeneity and molecular complexity of the mammalian nervous system relies substantially on the dynamic nature and spatiotemporal patterning of gene expression. These features of gene expression are achieved in part through mechanisms involving various epigenetic processes such as DNA methylation, post-translational histone modifications, and non-coding RNA activity, amongst others. In concert, another regulatory layer by which RNA bases and sugar residues are chemically modified enhances neuronal transcriptome complexity. Similar RNA modifications in other systems collectively constitute the cellular epitranscriptome that integrates and impacts various physiological processes. The epitranscriptome is dynamic and is reshaped constantly to regulate vital processes such as development, differentiation and stress responses. Perturbations of the epitranscriptome can lead to various pathogenic conditions, including cancer, cardiovascular abnormalities and neurological diseases. Recent advances in next-generation sequencing technologies have enabled us to identify and locate modified bases/sugars on different RNA species. These RNA modifications modulate the stability, transport and, most importantly, translation of RNA. In this review, we discuss the formation and functions of some frequently observed RNA modifications—including methylations of adenine and cytosine bases, and isomerization of uridine to pseudouridine—at various layers of RNA metabolism, together with their contributions to abnormal physiological conditions that can lead to various neurodevelopmental and neurological disorders.

Alternative polyadenylation: An enigma of transcript length variation in health and disease

Alternative polyadenylation (APA) is a molecular mechanism during a pre-mRNA processing that involves usage of more than one polyadenylation site (PA-site) generating transcripts of varying length from a single gene. The location of a PA-site affects transcript length and coding potential of an mRNA contributing to both mRNA and protein diversification. This variation in the transcript length affects mRNA stability and translation, mRNA subcellular and tissue localization, and protein function. APA is now considered as an important regulatory mechanism in the pathophysiology of human diseases. An important consequence of the changes in the length of 3'-untranslated region (UTR) from disease-induced APA is altered protein expression. Yet, the relationship between 3'-UTR length and protein expression remains a paradox in a majority of diseases. Here, we review occurrence of APA, mechanism of PA-site selection, and consequences of transcript length variation in different diseases. Emerging evidence reveals coordinated involvement of core RNA processing factors including poly(A) polymerases in the PA-site selection in diseases-associated APAs. Targeting such APA regulators will be therapeutically significant in combating drug resistance in cancer and other complex diseases. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease Translation > Regulation.

flexiMAP: a regression-based method for discovering differential alternative polyadenylation events in standard RNA-seq data

Motivation

We present flexible Modeling of Alternative PolyAdenylation (flexiMAP), a new beta-regression-based method implemented in R, for discovering differential alternative polyadenylation events in standard RNA-seq data.

Results

We show, using both simulated and real data, that flexiMAP exhibits a good balance between specificity and sensitivity and compares favourably to existing methods, especially at low fold changes. In addition, the tests on simulated data reveal some hitherto unrecognized caveats of existing methods. Importantly, flexiMAP allows modeling of multiple known covariates that often confound the results of RNA-seq data analysis.

Availability and implementation

The flexiMAP R package is available at: https://github.com/kszkop/flexiMAP. Scripts and data to reproduce the analysis in this paper are available at: https://doi.org/10.5281/zenodo.3689788.

Supplementary information

Supplementary data are available at Bioinformatics online.

Brain-immune interactions in neuropsychiatric disorders: Lessons from transcriptome studies for molecular targeting

Understanding the pathophysiological mechanisms of neuropsychiatric disorders has been a challenging quest for neurobiologists. Recent years have witnessed enormous technological advances in the field of neuroimmunology, blurring boundaries between the central nervous system and the periphery. Consequently, the discipline has expanded to cover interactions between the nervous and immune systems in health and diseases. The complex interplay between the peripheral and central immune pathways in neuropsychiatric disorders has recently been documented in various studies, but the genetic determinants remain elusive. Recent transcriptome studies have identified dysregulated genes involved in peripheral immune cell activation, blood-brain barrier integrity, glial cell activation, and synaptic plasticity in major depressive disorder, bipolar disorder, autism spectrum disorder, and schizophrenia. Herein, the key transcriptomic techniques applied in investigating differentially expressed genes and pathways responsible for altered brain-immune interactions in neuropsychiatric disorders are discussed. The application of transcriptomics that can aid in identifying molecular targets in various neuropsychiatric disorders is highlighted.

A missense mutation in the CSTF2 gene that impairs the function of the RNA recognition motif and causes defects in 3' end processing is associated with intellectual disability in humans

CSTF2 encodes an RNA-binding protein that is essential for mRNA cleavage and polyadenylation (C/P). No disease-associated mutations have been described for this gene. Here, we report a mutation in the RNA recognition motif (RRM) of CSTF2 that changes an aspartic acid at position 50 to alanine (p.D50A), resulting in intellectual disability in male patients. In mice, this mutation was sufficient to alter polyadenylation sites in over 1300 genes critical for brain development. Using a reporter gene assay, we demonstrated that C/P efficiency of CSTF2^D50A was lower than wild type. To account for this, we determined that p.D50A changed locations of amino acid side chains altering RNA binding sites in the RRM. The changes modified the electrostatic potential of the RRM leading to a greater affinity for RNA. These results highlight the significance of 3′ end mRNA processing in expression of genes important for brain plasticity and neuronal development.