Alternative polyadenylation transcriptome-wide association study identifies APA-linked susceptibility genes in brain disorders

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Leveraging Random Effects in Cistrome-Wide Association Studies for Decoding the Genetic Determinants of Prostate Cancer

Cistrome-wide association studies (CWAS) are pivotal for identifying genetic determinants of diseases by correlating genetically regulated cistrome states with phenotypes. Traditional CWAS typically develops a model based on cistrome and genotype data to associate predicted cistrome states with phenotypes. The random effect cistrome-wide association study (RECWAS), reevaluates the necessity of cistrome state prediction in CWAS. RECWAS utilizes either a linear model or marginal effect for initial feature selection, followed by kernel-based feature aggregation for association testing is introduced. Through simulations and analysis of prostate cancer data, a thorough evaluation of CWAS and RECWAS is conducted. The results suggest that RECWAS offers improved power compared to traditional CWAS, identifying additional genomic regions associated with prostate cancer. CWAS identified 102 significant regions, while RECWAS found 50 additional significant regions compared to CWAS, many of which are validated. Validation encompassed a range of biological evidence, including risk signals from the GWAS catalog, susceptibility genes from the DisGeNET database, and enhancer-domain scores. RECWAS consistently demonstrated improved performance over traditional CWAS in identifying genomic regions associated with prostate cancer. These findings demonstrate the benefits of incorporating kernel methods into CWAS and provide new insights for genetic discovery in complex diseases.

Large-Scale Alternative Polyadenylation-Wide Association Studies to Identify Putative Cancer Susceptibility Genes

Alternative polyadenylation (APA) modulates mRNA processing in the 3'-untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. In this study, we conducted large APA-wide association studies to investigate associations between APA levels and cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA sequencing data from 1,337 samples from the Genotype-Tissue Expression project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European ancestry populations: breast, ovarian, prostate, colorectal, lung, and pancreatic cancers. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3'-UTR APA quantitative trait loci and colocalization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3'-UTR variants demonstrated that the risk alleles of 3'-UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the posttranscriptional activities of their target genes compared with reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers. Significance: Systematic evaluation of associations of alternative polyadenylation with cancer risk reveals 58 putative susceptibility genes, highlighting the contribution of genetically regulated alternative polyadenylation of 3'UTRs to genetic susceptibility to cancer.

Sleep and diurnal alternative polyadenylation sites associated with human APA-linked brain disorders

Disruption of sleep and circadian rhythms are a comorbid feature of many pathologies, and can negatively influence many health conditions, including neurodegenerative disease, metabolic illness, cancer, and various neurological disorders. Genetic association studies linking sleep and circadian disturbances with disease susceptibility have mainly focused on changes in gene expression due to mutations, such as single-nucleotide polymorphisms. The interaction between sleep and/or circadian rhythms with the use of Alternative Polyadenylation (APA) has been largely undescribed, particularly in the context of other disorders. APA is a process that generates various transcript isoforms of the same gene affecting its mRNA translation, stability, localization, and subsequent function. Here we identified unique APAs expressed in rat brain over time-of-day, immediately following sleep deprivation, and the subsequent recovery period. From these data, we performed a secondary analysis of these sleep- or time-of-day associated PASs with recently described APA-linked human brain disorder susceptibility genes.

Identification, Design, and Application of Noncoding Cis-Regulatory Elements

Cis-regulatory elements (CREs) play a pivotal role in orchestrating interactions with trans-regulatory factors such as transcription factors, RNA-binding proteins, and noncoding RNAs. These interactions are fundamental to the molecular architecture underpinning complex and diverse biological functions in living organisms, facilitating a myriad of sophisticated and dynamic processes. The rapid advancement in the identification and characterization of these regulatory elements has been marked by initiatives such as the Encyclopedia of DNA Elements (ENCODE) project, which represents a significant milestone in the field. Concurrently, the development of CRE detection technologies, exemplified by massively parallel reporter assays, has progressed at an impressive pace, providing powerful tools for CRE discovery. The exponential growth of multimodal functional genomic data has necessitated the application of advanced analytical methods. Deep learning algorithms, particularly large language models, have emerged as invaluable tools for deconstructing the intricate nucleotide sequences governing CRE function. These advancements facilitate precise predictions of CRE activity and enable the de novo design of CREs. A deeper understanding of CRE operational dynamics is crucial for harnessing their versatile regulatory properties. Such insights are instrumental in refining gene therapy techniques, enhancing the efficacy of selective breeding programs, pushing the boundaries of genetic innovation, and opening new possibilities in microbial synthetic biology.

Susceptibility gene identification and risk evaluation model construction by transcriptome-wide association analysis for salt sensitivity of blood pressure

BACKGROUND

Salt sensitivity of blood pressure (SSBP) is an intermediate phenotype of hypertension and is a predictor of long-term cardiovascular events and death. However, the genetic structures of SSBP are uncertain, and it is difficult to precisely diagnose SSBP in population. So, we aimed to identify genes related to susceptibility to the SSBP, construct a risk evaluation model, and explore the potential functions of these genes.

METHODS AND RESULTS

A genome-wide association study of the systemic epidemiology of salt sensitivity (EpiSS) cohort was performed to obtain summary statistics for SSBP. Then, we conducted a transcriptome-wide association study (TWAS) of 12 tissues using FUSION software to predict the genes associated with SSBP and verified the genes with an mRNA microarray. The potential roles of the genes were explored. Risk evaluation models of SSBP were constructed based on the serial P value thresholds of polygenetic risk scores (PRSs), polygenic transcriptome risk scores (PTRSs) and their combinations of the identified genes and genetic variants from the TWAS. The TWAS revealed that 2605 genes were significantly associated with SSBP. Among these genes, 69 were differentially expressed according to the microarray analysis. The functional analysis showed that the genes identified in the TWAS were enriched in metabolic process pathways. The PRSs were correlated with PTRSs in the heart atrial appendage, adrenal gland, EBV-transformed lymphocytes, pituitary, artery coronary, artery tibial and whole blood. Multiple logistic regression models revealed that a PRS of P < 0.05 had the best predictive ability compared with other PRSs and PTRSs. The combinations of PRSs and PTRSs did not significantly increase the prediction accuracy of SSBP in the training and validation datasets.

CONCLUSIONS

Several known and novel susceptibility genes for SSBP were identified via multitissue TWAS analysis. The risk evaluation model constructed with the PRS of susceptibility genes showed better diagnostic performance than the transcript levels, which could be applied to screen for SSBP high-risk individuals.

PABPN1 loss-of-function causes APA-shift in oculopharyngeal muscular dystrophy

Alternative polyadenylation (APA) at the 3' UTR of transcripts contributes to the cell transcriptome. APA is suppressed by the nuclear RNA-binding protein PABPN1. Aging-associated reduced PABPN1 levels in skeletal muscles lead to muscle wasting. Muscle weakness in oculopharyngeal muscular dystrophy (OPMD) is caused by short alanine expansion in PABPN1 exon1. The expanded PABPN1 forms nuclear aggregates, an OPMD hallmark. Whether the expanded PABPN1 affects APA and how it contributes to muscle pathology is unresolved. To investigate these questions, we developed a procedure including RNA library preparation and a simple pipeline calculating the APA-shift ratio as a readout for PABPN1 activity. Comparing APA-shift results to previously published PAS utilization and APA-shift results, we validated this procedure. The procedure was then applied on the OPMD cell model and on RNA from OPMD muscles. APA-shift was genome-wide in the mouse OPMD model, primarily affecting muscle transcripts. In OPMD individuals, APA-shift was enriched with muscle transcripts. In an OPMD cell model APA-shift was not significant. APA-shift correlated with reduced expression levels of a subset of PABPN1 isoforms, whereas the expression of the expanded PABPN1 did not correlate with APA-shift. PABPN1 activity is not affected by the expression of expanded PABPN1, but rather by reduced PABPN1 expression levels. In muscles, PABPN1 activity initially affects muscle transcripts. We suggest that muscle weakness in OPMD is caused by PABPN1 loss-of-function leading to APA-shift that primarily affects in muscle transcripts.

The deubiquitinase function of ataxin-3 and its role in the pathogenesis of Machado-Joseph disease and other diseases

Machado-Joseph disease (MJD) is a devastating and incurable neurodegenerative disease characterised by progressive ataxia, difficulty speaking and swallowing. Consequently, affected individuals ultimately become wheelchair dependent, require constant care, and face a shortened life expectancy. The monogenic cause of MJD is expansion of a trinucleotide (CAG) repeat region within the ATXN3 gene, which results in polyglutamine (polyQ) expansion within the resultant ataxin-3 protein. While it is well established that the ataxin-3 protein functions as a deubiquitinating (DUB) enzyme and is therefore critically involved in proteostasis, several unanswered questions remain regarding the impact of polyQ expansion in ataxin-3 on its DUB function. Here we review the current literature surrounding ataxin-3's DUB function, its DUB targets, and what is known regarding the impact of polyQ expansion on ataxin-3's DUB function. We also consider the potential neuroprotective effects of ataxin-3's DUB function, and the intersection of ataxin-3's role as a DUB enzyme and regulator of gene transcription. Ataxin-3 is the principal pathogenic protein in MJD and also appears to be involved in cancer. As aberrant deubiquitination has been linked to both neurodegeneration and cancer, a comprehensive understanding of ataxin-3's DUB function is important for elucidating potential therapeutic targets in these complex conditions. In this review, we aim to consolidate knowledge of ataxin-3 as a DUB and unveil areas for future research to aid therapeutic targeting of ataxin-3's DUB function for the treatment of MJD and other diseases.

FMRP-mediated spatial regulation of physiologic NMD targets in neuronal cells

In non-polarized cells, nonsense-mediated mRNA decay (NMD) generally begins during the translation of newly synthesized mRNAs after the mRNAs are exported to the cytoplasm. Binding of the FMRP translational repressor to UPF1 on NMD targets mainly inhibits NMD. However, in polarized cells like neurons, FMRP additionally localizes mRNAs to cellular projections. Here, we review the literature and evaluate available transcriptomic data to conclude that, in neurons, the translation of physiologic NMD targets bound by FMRP is partially inhibited until the mRNAs localize to projections. There, FMRP displacement in response to signaling induces a burst in protein synthesis followed by rapid mRNA decay.

Identification of sleep and circadian alternative polyadenylation sites associated with APA-linked human brain disorders

Sleep and circadian rhythm disruptions are comorbid features of many pathologies and can negatively influence numerous health conditions, including degenerative diseases, metabolic illnesses, cancer, and various neurological disorders. Genetic association studies linking sleep and circadian disturbances with disease susceptibility have mainly focused on changes in gene expression due to mutations, such as single-nucleotide polymorphisms. Thus, associations between sleep and/or circadian rhythm and alternative polyadenylation (APA), particularly in the context of other health challenges, are largely undescribed. APA is a process that generates various transcript isoforms from the same gene, resulting in effects on mRNA translation, stability, localization, and subsequent function. Here, we have identified unique APAs in rat brain that exhibit time-of-day-dependent oscillations in expression as well as APAs that are altered by sleep deprivation and the subsequent recovery period. Genes affected by APA usage include Mapt/Tau, Ntrk2, Homer1A, Sin3band Sorl. Sorl1 has two APAs which cycle with a 24 h period, one additional APA cycles with a 12 h period and one more that is reduced during recovery sleep. Finally, we compared sleep- or circadian-associated APAs with recently described APA-linked brain disorder susceptibility genes and found 46 genes in common.

Neuroepigenetics of ageing and neurodegeneration-associated dementia: An updated review

Gene expression is tremendously altered in the brain during memory acquisition, recall, and forgetfulness. However, non-genetic factors, including environmental elements, epigenetic changes, and lifestyle, have grabbed significant attention in recent years regarding the etiology of neurodegenerative diseases (NDD) and age-associated dementia. Epigenetic modifications are essential in regulating gene expression in all living organisms in a DNA sequence-independent manner. The genes implicated in ageing and NDD-related memory disorders are epigenetically regulated by processes such as DNA methylation, histone acetylation as well as messenger RNA editing machinery. The physiological and optimal state of the epigenome, especially within the CNS of humans, plays an intricate role in helping us adjust to the changing environment, and alterations in it cause many brain disorders, but the mechanisms behind it still need to be well understood. When fully understood, these epigenetic landscapes could act as vital targets for pharmacogenetic rescue strategies for treating several diseases, including neurodegeneration- and age-induced dementia. Keeping this objective in mind, this updated review summarises the epigenetic changes associated with age and neurodegeneration-associated dementia.

Coordination of alternative splicing and alternative polyadenylation revealed by targeted long read sequencing

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

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

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

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