Search

In Alzheimer's disease (AD), the microtubule-binding protein tau becomes abnormally hyperphosphorylated and aggregated in selective brain regions such as the cortex and hippocampus 1-3 . However, other brain regions like the cerebellum and brain stem remain largely intact despite the universal expression of tau throughout the brain. Here, we found that an understudied splice isoform of tau termed "big tau" is significantly more abundant in the brain regions less vulnerable to tau pathology compared to tau pathology-vulnerable regions. We used various cellular and animal models to demonstrate that big tau possesses multiple properties that can resist AD-related pathological changes. Importantly, human AD patients show a higher expression level of pathology-resisting big tau in the cerebellum, the brain region spared from tau pathology. Our study examines the unique properties of big tau, expanding our current understanding of tau pathophysiology. Altogether, our data suggest that alternative splicing to favor big tau is a viable strategy to modulate tau pathology.

Background: Regulation of the thermogenic response by brown adipose tissue (BAT) is an important component of energy homeostasis with implications for the treatment of obesity and diabetes. Our preliminary analyses of RNA-Seq data uncovered many nodes representing epigenetic modifiers that are altered in BAT in response to chronic thermogenic activation. Thus, we hypothesized that chronic thermogenic activation broadly alters epigenetic modifications of DNA and histones in BAT., Results: Motivated to understand how BAT function is regulated epigenetically, we developed a novel method for the first-ever unbiased top-down proteomic quantitation of histone modifications in BAT and validated our results with a multi-omic approach. To test our hypothesis, wildtype male C57BL/6J mice were housed under chronic conditions of thermoneutral temperature (TN, 28°C), mild cold/room temperature (RT, 22°C), or severe cold (SC, 8°C) and BAT was analyzed for DNA methylation and histone modifications. Methylation of promoters and intragenic regions in genomic DNA decrease in response to chronic cold exposure. Integration of DNA methylation and RNA expression datasets suggest a role for epigenetic modification of DNA in regulation of gene expression in response to cold. In response to cold housing, we observe increased bulk acetylation of histones H3.2 and H4, increased histone H3.2 proteoforms with di- and trimethylation of lysine 9 (K9me2 and K9me3), and increased histone H4 proteoforms with acetylation of lysine 16 (K16ac) in BAT., Conclusions: Our results reveal global epigenetically-regulated transcriptional "on" and "off" signals in murine BAT in response to varying degrees of chronic cold stimuli and establish a novel methodology to quantitatively study histones in BAT, allowing for direct comparisons to decipher mechanistic changes during the thermogenic response. Additionally, we make histone PTM and proteoform quantitation, RNA splicing, RRBS, and transcriptional footprint datasets available as a resource for future research., (© 2024. The Author(s).)

Cellular responses to the steroid hormones, estrogen (E2), and progesterone (P4) are governed by their cognate receptor's transcriptional output. However, the feed-forward mechanisms that shape cell-type-specific transcriptional fulcrums for steroid receptors are unidentified. Herein, we found that a common feed-forward mechanism between GREB1 and steroid receptors regulates the differential effect of GREB1 on steroid hormones in a physiological or pathological context. In physiological (receptive) endometrium, GREB1 controls P4-responses in uterine stroma, affecting endometrial receptivity and decidualization, while not affecting E2-mediated epithelial proliferation. Of mechanism, progesterone-induced GREB1 physically interacts with the progesterone receptor, acting as a cofactor in a positive feedback mechanism to regulate P4-responsive genes. Conversely, in endometrial pathology (endometriosis), E2-induced GREB1 modulates E2-dependent gene expression to promote the growth of endometriotic lesions in mice. This differential action of GREB1 exerted by a common feed-forward mechanism with steroid receptors advances our understanding of mechanisms that underlie cell- and tissue-specific steroid hormone actions., (© 2024. The Author(s).)

Alternative polyadenylation (APA) is a key post-transcriptional regulatory mechanism; yet, its regulation and impact on human diseases remain understudied. Existing bulk RNA sequencing (RNA-seq)-based APA methods predominantly rely on predefined annotations, severely impacting their ability to decode novel tissue- and disease-specific APA changes. Furthermore, they only account for the most proximal and distal cleavage and polyadenylation sites (C/PASs). Deconvoluting overlapping C/PASs and the inherent noisy 3' UTR coverage in bulk RNA-seq data pose additional challenges. To overcome these limitations, we introduce PolyAMiner-Bulk, an attention-based deep learning algorithm that accurately recapitulates C/PAS sequence grammar, resolves overlapping C/PASs, captures non-proximal-to-distal APA changes, and generates visualizations to illustrate APA dynamics. Evaluation on multiple datasets strongly evinces the performance merit of PolyAMiner-Bulk, accurately identifying more APA changes compared with other methods. With the growing importance of APA and the abundance of bulk RNA-seq data, PolyAMiner-Bulk establishes a robust paradigm of APA analysis., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease with limited effective treatment options, potentiating the importance of uncovering novel drug targets. Here, we target cleavage and polyadenylation specificity factor 3 (CPSF3), the 3' endonuclease that catalyzes mRNA cleavage during polyadenylation and histone mRNA processing. We find that CPSF3 is highly expressed in PDAC and is associated with poor prognosis. CPSF3 knockdown blocks PDAC cell proliferation and colony formation in vitro and tumor growth in vivo. Chemical inhibition of CPSF3 by the small molecule JTE-607 also attenuates PDAC cell proliferation and colony formation, while it has no effect on cell proliferation of nontransformed immortalized control pancreatic cells. Mechanistically, JTE-607 induces transcriptional readthrough in replication-dependent histones, reduces core histone expression, destabilizes chromatin structure, and arrests cells in the S-phase of the cell cycle. Therefore, CPSF3 represents a potential therapeutic target for the treatment of PDAC., (© 2024 Alahmari et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Matrin-3 (MATR3) is an RNA-binding protein implicated in neurodegenerative and neurodevelopmental diseases. However, little is known regarding the role of MATR3 in cryptic splicing within the context of functional genes and how disease-associated variants impact this function. We show that loss of MATR3 leads to cryptic exon inclusion in many transcripts. We reveal that ALS-linked S85C pathogenic variant reduces MATR3 solubility but does not impair RNA binding. In parallel, we report a novel neurodevelopmental disease-associated M548T variant, located in the RRM2 domain, which reduces protein solubility and impairs RNA binding and cryptic splicing repression functions of MATR3. Altogether, our research identifies cryptic events within functional genes and demonstrates how disease-associated variants impact MATR3 cryptic splicing repression function., (© 2024 The Authors. FEBS Letters published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Purpose: Gliomas and their surrounding microenvironment constantly interact to promote tumorigenicity, yet the underlying posttranscriptional regulatory mechanisms that govern this interplay are poorly understood., Methods: Utilizing our established PAC-seq approach and PolyAMiner bioinformatic analysis pipeline, we deciphered the NUDT21-mediated differential APA dynamics in glioma cells., Results: We identified LAMC1 as a critical NUDT21 alternative polyadenylation (APA) target, common in several core glioma-driving signaling pathways. qRT-PCR analysis confirmed that NUDT21-knockdown in glioma cells results in the preferred usage of the proximal polyA signal (PAS) of LAMC1. Functional studies revealed that NUDT21-knockdown-induced 3'UTR shortening of LAMC1 is sufficient to cause translational gain, as LAMC1 protein is upregulated in these cells compared to their respective controls. We demonstrate that 3'UTR shortening of LAMC1 after NUDT21 knockdown removes binding sites for miR-124/506, thereby relieving potent miRNA-based repression of LAMC1 expression. Remarkably, we report that the knockdown of NUDT21 significantly promoted glioma cell migration and that co-depletion of LAMC1 with NUDT21 abolished this effect. Lastly, we observed that LAMC1 3'UTR shortening predicts poor prognosis of low-grade glioma patients from The Cancer Genome Atlas., Conclusion: This study identifies NUDT21 as a core alternative polyadenylation factor that regulates the tumor microenvironment through differential APA and loss of miR-124/506 inhibition of LAMC1. Knockdown of NUDT21 in GBM cells mediates 3'UTR shortening of LAMC1, contributing to an increase in LAMC1, increased glioma cell migration/invasion, and a poor prognosis., (© 2023. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Mutations in MeCP2 result in a crippling neurological disease, but we lack a lucid picture of MeCP2's molecular role. Individual transcriptomic studies yield inconsistent differentially expressed genes. To overcome these issues, we demonstrate a methodology to analyze all modern public data. We obtained relevant raw public transcriptomic data from GEO and ENA, then homogeneously processed it (QC, alignment to reference, differential expression analysis). We present a web portal to interactively access the mouse data, and we discovered a commonly perturbed core set of genes that transcends the limitations of any individual study. We then found functionally distinct, consistently up- and downregulated subsets within these genes and some bias to their location. We present this common core of genes as well as focused cores for up, down, cell fraction models, and some tissues. We observed enrichment for this mouse core in other species MeCP2 models and observed overlap with ASD models. By integrating and examining transcriptomic data at scale, we have uncovered the true picture of this dysregulation. The vast scale of these data enables us to analyze signal-to-noise, evaluate a molecular signature in an unbiased manner, and demonstrate a framework for future disease focused informatics work.

Alternative polyadenylation (APA) creates distinct transcripts from the same gene by cleaving the pre-mRNA at poly(A) sites that can lie within the 3' untranslated region (3'UTR), introns, or exons. Most studies focus on APA within the 3'UTR; however, here, we show that CPSF6 insufficiency alters protein levels and causes a developmental syndrome by deregulating APA throughout the transcript. In neonatal humans and zebrafish larvae, CPSF6 insufficiency shifts poly(A) site usage between the 3'UTR and internal sites in a pathway-specific manner. Genes associated with neuronal function undergo mostly intronic APA, reducing their expression, while genes associated with heart and skeletal function mostly undergo 3'UTR APA and are up-regulated. This suggests that, under healthy conditions, cells toggle between internal and 3'UTR APA to modulate protein expression.

PGC-1α is well established as a metazoan transcriptional coactivator of cellular adaptation in response to stress. However, the mechanisms by which PGC-1α activates gene transcription are incompletely understood. Here, we report that PGC-1α serves as a scaffold protein that physically and functionally connects the DNA-binding protein estrogen-related receptor α (ERRα), cap-binding protein 80 (CBP80), and Mediator to overcome promoter-proximal pausing of RNAPII and transcriptionally activate stress-response genes. We show that PGC-1α promotes pausing release in a two-arm mechanism (1) by recruiting the positive transcription elongation factor b (P-TEFb) and (2) by outcompeting the premature transcription termination complex Integrator. Using mice homozygous for five amino acid changes in the CBP80-binding motif (CBM) of PGC-1α that destroy CBM function, we show that efficient differentiation of primary myoblasts to myofibers and timely skeletal muscle regeneration after injury require PGC-1α binding to CBP80. Our findings reveal how PGC-1α activates stress-response gene transcription in a previously unanticipated pre-mRNA quality-control pathway., Competing Interests: Declaration of interests K.A. is a member of the Advisory Board of Molecular Cell., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Recent genome-wide association studies corroborate classical research on developmental programming indicating that obesity is primarily a neurodevelopmental disease strongly influenced by nutrition during critical ontogenic windows. Epigenetic mechanisms regulate neurodevelopment; however, little is known about their role in establishing and maintaining the brain's energy balance circuitry. We generated neuron and glia methylomes and transcriptomes from male and female mouse hypothalamic arcuate nucleus, a key site for energy balance regulation, at time points spanning the closure of an established critical window for developmental programming of obesity risk. We find that postnatal epigenetic maturation is markedly cell type and sex specific and occurs in genomic regions enriched for heritability of body mass index in humans. Our results offer a potential explanation for both the limited ontogenic windows for and sex differences in sensitivity to developmental programming of obesity and provide a rich resource for epigenetic analyses of developmental programming of energy balance.

Alternative polyadenylation (APA) regulates gene expression by cleavage and addition of poly(A) sequence at different polyadenylation sites (PAS) in 3'UTR, thus, generating transcript isoforms with different lengths. Cleavage stimulating factor 64 (CstF64) is an APA regulator which plays a role in PAS selection and determines the length of 3'UTR. CstF64 favors the use of proximal PAS, resulting in 3'UTR shortening, which enhances the protein expression by increasing the stability of the target genes. The aim of this study is to investigate the role of CstF64 in cardiac fibrosis, a key event leading to heart failure (HF). We determined the expression of CstF64, key profibrotic genes, and their 3'UTR changes by calculating distal PAS (dPAS) usage in left ventricular (LV) tissues and cardiac fibroblasts from HF patients. CstF64 was upregulated in HF LV tissues and cardiac fibroblasts along with increased deposition of fibrosis genes such as COL1A and FN1 and significant shortening in their 3'UTR. In addition, HF cardiac fibroblasts showed increased transforming growth factor receptor β1 (TGFβR1) expression consistent with significant shortening in 3'UTR of TGFβR1. Upon knockdown of CstF64 from HF fibroblasts, downregulation in pro-fibrotic genes corresponding to lengthening in their 3'UTR was observed. Our finding suggests an important role of CstF64 in myofibroblast activation and promotion of cardiac fibrosis during HF through APA. Therefore, targeting CstF64 mediated RNA processing approach in human HF could provide a new therapeutic treatment strategy for limiting fibrotic remodeling., Competing Interests: Declaration of competing interest There are no conflicts of interest., (Copyright © 2022 Elsevier Inc. All rights reserved.)

MeCP2 is associated with Rett syndrome (RTT), MECP2 duplication syndrome, and a number of conditions with isolated features of these diseases, including autism, intellectual disability, and motor dysfunction. MeCP2 is known to broadly bind methylated DNA, but the precise molecular mechanism driving disease pathogenesis remains to be determined. Using proximity-dependent biotinylation (BioID), we identified a transcription factor 20 (TCF20) complex that interacts with MeCP2 at the chromatin interface. Importantly, RTT-causing mutations in MECP2 disrupt this interaction. TCF20 and MeCP2 are highly coexpressed in neurons and coregulate the expression of key neuronal genes. Reducing Tcf20 partially rescued the behavioral deficits caused by MECP2 overexpression, demonstrating a functional relationship between MeCP2 and TCF20 in MECP2 duplication syndrome pathogenesis. We identified a patient exhibiting RTT-like neurological features with a missense mutation in the PHF14 subunit of the TCF20 complex that abolishes the MeCP2-PHF14-TCF20 interaction. Our data demonstrate the critical role of the MeCP2-TCF20 complex for brain function., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

In the messenger RNA (mRNA) maturation process, the 3'-end of pre-mRNA is cleaved and a poly(A) sequence is added, this is an important determinant of mRNA stability and its cellular functions. More than 60%-70% of human genes have three or more polyadenylation (APA) sites and can be cleaved at different sites, generating mRNA transcripts of varying lengths. This phenomenon is termed as alternative cleavage and polyadenylation (APA) and it plays role in key biological processes like gene regulation, cell proliferation, senescence, and also in various human diseases. Loss of regulatory microRNA binding sites and interactions with RNA-binding proteins leading to APA are largely investigated in human diseases. However, the functions of the core APA machinery and related factors during disease conditions remain largely unknown. In this review, we discuss the roles of polyadenylation machinery in relation to brain disease, cardiac failure, pulmonary fibrosis, cancer, infectious conditions, and other human diseases. Collectively, we believe this review will be a useful avenue for understanding the emerging role of APA in the pathobiology of various human diseases., (© 2021 Wiley Periodicals LLC.)

A critical question in neurodegeneration is why the accumulation of disease-driving proteins causes selective neuronal loss despite their brain-wide expression. In Spinocerebellar ataxia type 1 (SCA1), accumulation of polyglutamine-expanded Ataxin-1 (ATXN1) causes selective degeneration of cerebellar and brainstem neurons. Previous studies revealed that inhibiting Msk1 reduces phosphorylation of ATXN1 at S776 as well as its levels leading to improved cerebellar function. However, there are no regulators that modulate ATXN1 in the brainstem-the brain region whose pathology is most closely linked to premature death. To identify new regulators of ATXN1, we performed genetic screens and identified a transcription factor-kinase axis (ZBTB7B-RSK3) that regulates ATXN1 levels. Unlike MSK1, RSK3 is highly expressed in the human and mouse brainstems where it regulates Atxn1 by phosphorylating S776. Reducing Rsk3 rescues brainstem-associated pathologies and deficits, and lowering Rsk3 and Msk1 together improves cerebellar and brainstem function in an SCA1 mouse model. Our results demonstrate that selective vulnerability of brain regions in SCA1 is governed by region-specific regulators of ATXN1, and targeting multiple regulators could rescue multiple degenerating brain areas., (© 2021 The Authors.)

An increasing number of investigations have established alternative polyadenylation (APA) as a key mechanism of gene regulation through altering the length of 3' untranslated region (UTR) and generating distinct mRNA termini. Further, appreciation for the significance of APA in disease contexts propelled the development of several 3' sequencing techniques. While these RNA sequencing technologies have advanced APA analysis, the intrinsic limitation of 3' read coverage and lack of appropriate computational tools constrain precise mapping and quantification of polyadenylation sites. Notably, Poly(A)-ClickSeq (PAC-seq) overcomes limiting factors such as poly(A) enrichment and 3' linker ligation steps using click-chemistry. Here we provide an updated PolyA-miner protocol, a computational approach to analyze PAC-seq or other 3'-Seq datasets. As a key practical constraint, we also provide a detailed account on the impact of sequencing depth on the number of detected polyadenylation sites and APA changes. This protocol is also updated to handle unique molecular identifiers used to address PCR duplication potentially observed in PAC-seq., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Over the past 15 years, investigations into alternative polyadenylation (APA) and its function in cellular physiology and pathology have greatly expanded due to the emergent appreciation of its key role in driving transcriptomic diversity. This growth has necessitated the development of new technologies capable of monitoring cleavage and polyadenylation events genome-wide. Advancements in approaches include both the creation of computational tools to re-analyze RNA-seq to identify APA events as well as targeted sequencing approaches customized to focus on the 3'-end of mRNA. Here we describe a streamlined protocol for polyA-Click-seq (PAC-seq), which utilizes click-chemistry to create mRNA 3'-ends sequencing libraries. Importantly, we offer additional considerations not present in our previous protocols including the use of spike-ins, unique molecular identifier primers, and guidance for appropriate depth of PAC-seq. In conjunction with the companion chapter on PolyA-miner (Yalamanchili et al., 2021) to computationally analyze PAC-seq data, we provide a complete experimental pipeline to analyze mRNA 3'-end usage in eukaryotic cells., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Background: Tau neurofibrillary tangle pathology characterizes Alzheimer's disease and other neurodegenerative tauopathies. Brain gene expression profiles can reveal mechanisms; however, few studies have systematically examined both the transcriptome and proteome or differentiated Tau- versus age-dependent changes., Methods: Paired, longitudinal RNA-sequencing and mass-spectrometry were performed in a Drosophila model of tauopathy, based on pan-neuronal expression of human wildtype Tau (Tau WT ) or a mutant form causing frontotemporal dementia (Tau R406W ). Tau-induced, differentially expressed transcripts and proteins were examined cross-sectionally or using linear regression and adjusting for age. Hierarchical clustering was performed to highlight network perturbations, and we examined overlaps with human brain gene expression profiles in tauopathy., Results: Tau WT induced 1514 and 213 differentially expressed transcripts and proteins, respectively. Tau R406W had a substantially greater impact, causing changes in 5494 transcripts and 697 proteins. There was a ~ 70% overlap between age- and Tau-induced changes and our analyses reveal pervasive bi-directional interactions. Strikingly, 42% of Tau-induced transcripts were discordant in the proteome, showing opposite direction of change. Tau-responsive gene expression networks strongly implicate innate immune activation. Cross-species analyses pinpoint human brain gene perturbations specifically triggered by Tau pathology and/or aging, and further differentiate between disease amplifying and protective changes., Conclusions: Our results comprise a powerful, cross-species functional genomics resource for tauopathy, revealing Tau-mediated disruption of gene expression, including dynamic, age-dependent interactions between the brain transcriptome and proteome.

Almost 70% of human genes undergo alternative polyadenylation (APA) and generate mRNA transcripts with varying lengths, typically of the 3' untranslated regions (UTR). APA plays an important role in development and cellular differentiation, and its dysregulation can cause neuropsychiatric diseases and increase cancer severity. Increasing awareness of APA's role in human health and disease has propelled the development of several 3' sequencing (3'Seq) techniques that allow for precise identification of APA sites. However, despite the recent data explosion, there are no robust computational tools that are precisely designed to analyze 3'Seq data. Analytical approaches that have been used to analyze these data predominantly use proximal to distal usage. With about 50% of human genes having more than two APA isoforms, current methods fail to capture the entirety of APA changes and do not account for non-proximal to non-distal changes. Addressing these key challenges, this study demonstrates PolyA-miner, an algorithm to accurately detect and assess differential alternative polyadenylation specifically from 3'Seq data. Genes are abstracted as APA matrices, and differential APA usage is inferred using iterative consensus non-negative matrix factorization (NMF) based clustering. PolyA-miner accounts for all non-proximal to non-distal APA switches using vector projections and reflects precise gene-level 3'UTR changes. It can also effectively identify novel APA sites that are otherwise undetected when using reference-based approaches. Evaluation on multiple datasets-first-generation MicroArray Quality Control (MAQC) brain and Universal Human Reference (UHR) PolyA-seq data, recent glioblastoma cell line NUDT21 knockdown Poly(A)-ClickSeq (PAC-seq) data, and our own mouse hippocampal and human stem cell-derived neuron PAC-seq data-strongly supports the value and protocol-independent applicability of PolyA-miner. Strikingly, in the glioblastoma cell line data, PolyA-miner identified more than twice the number of genes with APA changes than initially reported. With the emerging importance of APA in human development and disease, PolyA-miner can significantly improve data analysis and help decode the underlying APA dynamics., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We previously showed that NUDT21 -spanning copy-number variations (CNVs) are associated with intellectual disability (Gennarino et al., 2015). However, the patients' CNVs also included other genes. To determine if reduced NUDT21 function alone can cause disease, we generated Nudt21 +/- mice to mimic NUDT21 -deletion patients. We found that although these mice have 50% reduced Nudt21 mRNA, they only have 30% less of its cognate protein, CFIm25. Despite this partial protein-level compensation, the Nudt21 +/- mice have learning deficits, cortical hyperexcitability, and misregulated alternative polyadenylation (APA) in their hippocampi. Further, to determine the mediators driving neural dysfunction in humans, we partially inhibited NUDT21 in human stem cell-derived neurons to reduce CFIm25 by 30%. This induced APA and protein level misregulation in hundreds of genes, a number of which cause intellectual disability when mutated. Altogether, these results show that disruption of NUDT21 -regulated APA events in the brain can cause intellectual disability., Competing Interests: CA, HY, PJ, Mv, AS, NE, AL, ML, BB, SH, QW, AS, JT, AM, EW, ZL No competing interests declared, HZ Senior editor, eLife, (© 2020, Alcott et al.)

In Alzheimer's disease (AD), spliceosomal proteins with critical roles in RNA processing aberrantly aggregate and mislocalize to Tau neurofibrillary tangles. We test the hypothesis that Tau-spliceosome interactions disrupt pre-mRNA splicing in AD. In human postmortem brain with AD pathology, Tau coimmunoprecipitates with spliceosomal components. In Drosophila, pan-neuronal Tau expression triggers reductions in multiple core and U1-specific spliceosomal proteins, and genetic disruption of these factors, including SmB, U1-70K, and U1A, enhances Tau-mediated neurodegeneration. We further show that loss of function in SmB, encoding a core spliceosomal protein, causes decreased survival, progressive locomotor impairment, and neuronal loss, independent of Tau toxicity. Lastly, RNA sequencing reveals a similar profile of mRNA splicing errors in SmB mutant and Tau transgenic flies, including intron retention and non-annotated cryptic splice junctions. In human brains, we confirm cryptic splicing errors in association with neurofibrillary tangle burden. Our results implicate spliceosome disruption and the resulting transcriptome perturbation in Tau-mediated neurodegeneration in AD., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Nuclear receptor corepressor 1 (NCOR1) and NCOR2 (also known as SMRT) regulate gene expression by activating histone deacetylase 3 through their deacetylase activation domain (DAD). We show that mice with DAD knock-in mutations have memory deficits, reduced anxiety levels, and reduced social interactions. Mice with NCOR1 and NORC2 depletion specifically in GABAergic neurons (NS-V mice) recapitulated the memory deficits and had reduced GABA A receptor subunit α2 (GABRA2) expression in lateral hypothalamus GABAergic (LH GABA ) neurons. This was associated with LH GABA neuron hyperexcitability and impaired hippocampal long-term potentiation, through a monosynaptic LH GABA to CA3 GABA projection. Optogenetic activation of this projection caused memory deficits, whereas targeted manipulation of LH GABA or CA3 GABA neuron activity reversed memory deficits in NS-V mice. We describe de novo variants in NCOR1, NCOR2 or HDAC3 in patients with intellectual disability or neurodevelopmental disorders. These findings identify a hypothalamus-hippocampus projection that may link endocrine signals with synaptic plasticity through NCOR-mediated regulation of GABA signaling.

RNA splicing entails the coordinated interaction of more than 150 proteins in the spliceosome, one of the most complex of the cell's molecular machines. We previously discovered that the RNA-binding motif protein 17 (RBM17), a component of the spliceosome, is essential for survival and cell maintenance. Here, we find that it interacts with the spliceosomal factors U2SURP and CHERP and that they reciprocally regulate each other's stability, both in mouse and in human cells. Individual knockdown of each of the three proteins induces overlapping changes in splicing and gene expression of transcripts enriched for RNA-processing factors. Our results elucidate the function of RBM17, U2SURP, and CHERP and link the activity of the spliceosome to the regulation of downstream RNA-binding proteins. These data support the hypothesis that, beyond driving constitutive splicing, spliceosomal factors can regulate alternative splicing of specific targets., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Recent studies have suggested that genes longer than 100 kb are more likely to be misregulated in neurological diseases associated with synaptic dysfunction, such as autism and Rett syndrome. These length-dependent transcriptional changes are modest in MeCP2-mutant samples, but, given the low sensitivity of high-throughput transcriptome profiling technology, here we re-evaluate the statistical significance of these results. We find that the apparent length-dependent trends previously observed in MeCP2 microarray and RNA-sequencing datasets disappear after estimating baseline variability from randomized control samples. This is particularly true for genes with low fold changes. We find no bias with NanoString technology, so this long gene bias seems to be particular to polymerase chain reaction amplification-based platforms. In contrast, authentic long gene effects, such as those caused by topoisomerase inhibition, can be detected even after adjustment for baseline variability. We conclude that accurate characterization of length-dependent (or other) trends requires establishing a baseline from randomized control samples.

Transposable elements (TEs) are DNA sequences which are capable of moving from one location to another and represent a large proportion (45%) of the human genome. TEs have functional roles in a variety of biological phenomena such as cancer, neurodegenerative disease, and aging. Rapid development in RNA-sequencing technology has enabled us, for the first time, to study the activity of TE at the systems level.However, efficient TE analysis tools are not yet developed. In this work, we developed SalmonTE, a fast and reliable pipeline for the quantification of TEs from RNA-seq data. We benchmarked our tool against TEtranscripts, a widely used TE quantification method, and three other quantification methods using several RNA-seq datasets from Drosophila melanogaster and human cell-line. We achieved 20 times faster execution speed without compromising the accuracy. This pipeline will enable the biomedical research community to quantify and analyze TEs from large amounts of data and lead to novel TE centric discoveries.

RNA sequencing (RNA-seq) is a high-throughput technology that provides unique insights into the transcriptome. It has a wide variety of applications in quantifying genes/isoforms and in detecting non-coding RNA, alternative splicing, and splice junctions. It is extremely important to comprehend the entire transcriptome for a thorough understanding of the cellular system. Several RNA-seq analysis pipelines have been proposed to date. However, no single analysis pipeline can capture dynamics of the entire transcriptome. Here, we compile and present a robust and commonly used analytical pipeline covering the entire spectrum of transcriptome analysis, including quality checks, alignment of reads, differential gene/transcript expression analysis, discovery of cryptic splicing events, and visualization. Challenges, critical parameters, and possible downstream functional analysis pipelines associated with each step are highlighted and discussed. This unit provides a comprehensive understanding of state-of-the-art RNA-seq analysis pipeline and a greater understanding of the transcriptome. © 2017 by John Wiley & Sons, Inc., (Copyright © 2017 John Wiley & Sons, Inc.)

Splicing regulation is an important step of post-transcriptional gene regulation. It is a highly dynamic process orchestrated by RNA-binding proteins (RBPs). RBP dysfunction and global splicing dysregulation have been implicated in many human diseases, but the in vivo functions of most RBPs and the splicing outcome upon their loss remain largely unexplored. Here we report that constitutive deletion of Rbm17, which encodes an RBP with a putative role in splicing, causes early embryonic lethality in mice and that its loss in Purkinje neurons leads to rapid degeneration. Transcriptome profiling of Rbm17-deficient and control neurons and subsequent splicing analyses using CrypSplice, a new computational method that we developed, revealed that more than half of RBM17-dependent splicing changes are cryptic. Importantly, RBM17 represses cryptic splicing of genes that likely contribute to motor coordination and cell survival. This finding prompted us to re-analyze published datasets from a recent report on TDP-43, an RBP implicated in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), as it was demonstrated that TDP-43 represses cryptic exon splicing to promote cell survival. We uncovered a large number of TDP-43-dependent splicing defects that were not previously discovered, revealing that TDP-43 extensively regulates cryptic splicing. Moreover, we found a significant overlap in genes that undergo both RBM17- and TDP-43-dependent cryptic splicing repression, many of which are associated with survival. We propose that repression of cryptic splicing by RBPs is critical for neuronal health and survival. CrypSplice is available at www.liuzlab.org/CrypSplice., (© The Author 2016. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Background & Aims: The transcription factor atonal homolog 1 (ATOH1) controls the fate of intestinal progenitors downstream of the Notch signaling pathway. Intestinal progenitors that escape Notch activation express high levels of ATOH1 and commit to a secretory lineage fate, implicating ATOH1 as a gatekeeper for differentiation of intestinal epithelial cells. Although some transcription factors downstream of ATOH1, such as SPDEF, have been identified to specify differentiation and maturation of specific cell types, the bona fide transcriptional targets of ATOH1 still largely are unknown. Here, we aimed to identify ATOH1 targets and to identify transcription factors that are likely to co-regulate gene expression with ATOH1., Methods: We used a combination of chromatin immunoprecipitation and messenger RNA-based high-throughput sequencing (ChIP-seq and RNA-seq), together with cell sorting and transgenic mice, to identify direct targets of ATOH1, and establish the epistatic relationship between ATOH1 and SPDEF., Results: By using unbiased genome-wide approaches, we identified more than 700 genes as ATOH1 transcriptional targets in adult small intestine and colon. Ontology analysis indicated that ATOH1 directly regulates genes involved in specification and function of secretory cells. De novo motif analysis of ATOH1 targets identified SPDEF as a putative transcriptional co-regulator of ATOH1. Functional epistasis experiments in transgenic mice show that SPDEF amplifies ATOH1-dependent transcription but cannot independently initiate transcription of ATOH1 target genes., Conclusions: This study unveils the direct targets of ATOH1 in the adult intestines and illuminates the transcriptional events that initiate the specification and function of intestinal secretory lineages.

Decreasing pH due to anthropogenic CO2 inputs, called ocean acidification (OA), can make coastal environments unfavorable for oysters. This is a serious socioeconomical issue for China which supplies >70% of the world's edible oysters. Here, we present an iTRAQ-based protein profiling approach for the detection and quantification of proteome changes under OA in the early life stage of a commercially important oyster, Crassostrea hongkongensis. Availability of complete genome sequence for the pacific oyster (Crassostrea gigas) enabled us to confidently quantify over 1500 proteins in larval oysters. Over 7% of the proteome was altered in response to OA at pHNBS 7.6. Analysis of differentially expressed proteins and their associated functional pathways showed an upregulation of proteins involved in calcification, metabolic processes, and oxidative stress, each of which may be important in physiological adaptation of this species to OA. The downregulation of cytoskeletal and signal transduction proteins, on the other hand, might have impaired cellular dynamics and organelle development under OA. However, there were no significant detrimental effects in developmental processes such as metamorphic success. Implications of the differentially expressed proteins and metabolic pathways in the development of OA resistance in oyster larvae are discussed. The MS proteomics data have been deposited to the ProteomeXchange with identifiers PXD002138 (http://proteomecentral.proteomexchange.org/dataset/PXD002138)., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Conventionally, overall gene expressions from microarrays are used to infer gene networks, but it is challenging to account splicing isoforms. High-throughput RNA Sequencing has made splice variant profiling practical. However, its true merit in quantifying splicing isoforms and isoform-specific exon expressions is not well explored in inferring gene networks. This study demonstrates SpliceNet, a method to infer isoform-specific co-expression networks from exon-level RNA-Seq data, using large dimensional trace. It goes beyond differentially expressed genes and infers splicing isoform network changes between normal and diseased samples. It eases the sample size bottleneck; evaluations on simulated data and lung cancer-specific ERBB2 and MAPK signaling pathways, with varying number of samples, evince the merit in handling high exon to sample size ratio datasets. Inferred network rewiring of well established Bcl-x and EGFR centered networks from lung adenocarcinoma expression data is in good agreement with literature. Gene level evaluations demonstrate a substantial performance of SpliceNet over canonical correlation analysis, a method that is currently applied to exon level RNA-Seq data. SpliceNet can also be applied to exon array data. SpliceNet is distributed as an R package available at http://www.jjwanglab.org/SpliceNet., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Inferring gene regulatory networks from gene expression data at whole genome level is still an arduous challenge, especially in higher organisms where the number of genes is large but the number of experimental samples is small. It is reported that the accuracy of current methods at genome scale significantly drops from Escherichia coli to Saccharomyces cerevisiae due to the increase in number of genes. This limits the applicability of current methods to more complex genomes, like human and mouse. Least absolute shrinkage and selection operator (LASSO) is widely used for gene regulatory network inference from gene expression profiles. However, the accuracy of LASSO on large genomes is not satisfactory. In this study, we apply two extended models of LASSO, L0 and L1/2 regularization models to infer gene regulatory network from both high-throughput gene expression data and transcription factor binding data in mouse embryonic stem cells (mESCs). We find that both the L0 and L1/2 regularization models significantly outperform LASSO in network inference. Incorporating interactions between transcription factors and their targets remarkably improved the prediction accuracy. Current study demonstrates the efficiency and applicability of these two models for gene regulatory network inference from integrative omics data in large genomes. The applications of the two models will facilitate biologists to study the gene regulation of higher model organisms in a genome-wide scale., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Motivation: Inferring gene-regulatory networks is very crucial in decoding various complex mechanisms in biological systems. Synthesis of a fully functional transcriptional factor/protein from DNA involves series of reactions, leading to a delay in gene regulation. The complexity increases with the dynamic delay induced by other small molecules involved in gene regulation, and noisy cellular environment. The dynamic delay in gene regulation is quite evident in high-temporal live cell lineage-imaging data. Although a number of gene-network-inference methods are proposed, most of them ignore the associated dynamic time delay., Results: Here, we propose DDGni (dynamic delay gene-network inference), a novel gene-network-inference algorithm based on the gapped local alignment of gene-expression profiles. The local alignment can detect short-term gene regulations, that are usually overlooked by traditional correlation and mutual Information based methods. DDGni uses 'gaps' to handle the dynamic delay and non-uniform sampling frequency in high-temporal data, like live cell imaging data. Our algorithm is evaluated on synthetic and yeast cell cycle data, and Caenorhabditis elegans live cell imaging data against other prominent methods. The area under the curve of our method is significantly higher when compared to other methods on all three datasets., Availability: The program, datasets and supplementary files are available at http://www.jjwanglab.org/DDGni/.

One of the fundamental goals of genetics is to understand gene functions and their associated phenotypes. To achieve this goal, in this study we developed a computational algorithm that uses orthology and protein-protein interaction information to infer gene-phenotype associations for multiple species. Furthermore, we developed a web server that provides genome-wide phenotype inference for six species: fly, human, mouse, worm, yeast, and zebrafish. We evaluated our inference method by comparing the inferred results with known gene-phenotype associations. The high Area Under the Curve values suggest a significant performance of our method. By applying our method to two human representative diseases, Type 2 Diabetes and Breast Cancer, we demonstrated that our method is able to identify related Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways. The web server can be used to infer functions and putative phenotypes of a gene along with the candidate genes of a phenotype, and thus aids in disease candidate gene discovery. Our web server is available at http://jjwanglab.org/PhenoPPIOrth.

Background: Large amounts of data are being generated by high-throughput genome sequencing methods. But the rate of the experimental functional characterization falls far behind. To fill the gap between the number of sequences and their annotations, fast and accurate automated annotation methods are required. Many methods, such as GOblet, GOFigure, and Gotcha, are designed based on the BLAST search. Unfortunately, the sequence coverage of these methods is low as they cannot detect the remote homologues. Adding to this, the lack of annotation specificity advocates the need to improve automated protein function prediction., Results: We designed a novel automated protein functional assignment method based on the neural response algorithm, which simulates the neuronal behavior of the visual cortex in the human brain. Firstly, we predict the most similar target protein for a given query protein and thereby assign its GO term to the query sequence. When assessed on test set, our method ranked the actual leaf GO term among the top 5 probable GO terms with accuracy of 86.93%., Conclusions: The proposed algorithm is the first instance of neural response algorithm being used in the biological domain. The use of HMM profiles along with the secondary structure information to define the neural response gives our method an edge over other available methods on annotation accuracy. Results of the 5-fold cross validation and the comparison with PFP and FFPred servers indicate the prominent performance by our method. The program, the dataset, and help files are available at http://www.jjwanglab.org/NRProF/.

We have demonstrated that the natural flavone acacetin selectively inhibits ultra-rapid delayed rectifier potassium current (I(Kur)) in human atria. However, molecular determinants of this ion channel blocker are unknown. The present study was designed to investigate the molecular determinants underlying the ability of acacetin to block hKv1.5 channels (coding I(Kur)) in human atrial myocytes using the whole-cell patch voltage-clamp technique to record membrane current in HEK 293 cells stably expressing the hKv1.5 gene or transiently expressing mutant hKv1.5 genes generated by site-directed mutagenesis. It was found that acacetin blocked hKv1.5 channels by binding to both closed and open channels. The blockade of hKv1.5 channels by acacetin was use- and frequency-dependent, and the IC(50) of acacetin for inhibiting hKv1.5 was 3.5, 3.1, 2.9, 2.1, and 1.7μM, respectively, at 0.2, 0.5, 1, 3, and 4Hz. The mutagenesis study showed that the hKv1.5 mutants V505A, I508A, and V512A in the S6-segment remarkably reduced the channel blocking properties by acacetin (IC(50), 29.5μM for V505A, 19.1μM for I508A, and 6.9μM for V512A). These results demonstrate the novel information that acacetin mainly blocks open hKv1.5 channels by binding to their S6 domain. The use- and rate-dependent blocking of hKv1.5 by acacetin is beneficial for anti-atrial fibrillation., (2011 Elsevier Ltd. All rights reserved.)

Motivation: The interaction between transcription factor (TF) and transcription factor binding site (TFBS) is essential for gene regulation. Mutation in either the TF or the TFBS may weaken their interaction and thus result in abnormalities. To maintain such vital interaction, a mutation in one of the interacting partners might be compensated by a corresponding mutation in its binding partner during the course of evolution. Confirming this co-evolutionary relationship will guide us in designing protein sequences to target a specific DNA sequence or in predicting TFBS for poorly studied proteins, or even correcting and rescuing disease mutations in clinical applications., Results: Based on six, publicly available, experimentally validated TF-TFBS binding datasets for the basic Helix-Loop-Helix (bHLH) family, Homeo family, High-Mobility Group (HMG) family and Transient Receptor Potential channels (TRP) family, we showed that the evolutions of the TFs and their TFBSs are significantly correlated across eukaryotes. We further developed a mutual information-based method to identify co-evolved protein residues and DNA bases. This research sheds light on the dynamic relationship between TF and TFBS during their evolution. The same principle and strategy can be applied to co-evolutionary studies on protein-DNA interactions in other protein families., Availability: All the datasets, scripts and other related files have been made freely available at: http://jjwanglab.org/co-evo., Contact: junwen@uw.edu., Supplementary Information: Supplementary data are available at Bioinformatics online.