Transite: A Computational Motif-Based Analysis Platform That Identifies RNA-Binding Proteins Modulating Changes in Gene Expression

ePRINT: exonuclease assisted mapping of protein-RNA interactions

RNA-binding proteins (RBPs) regulate key aspects of RNA processing including alternative splicing, mRNA degradation and localization by physically binding RNA molecules. Current methods to map these interactions, such as CLIP, rely on purifying single proteins at a time. Our new method, ePRINT, maps RBP-RNA interaction networks on a global scale without purifying individual RBPs. ePRINT uses exoribonuclease XRN1 to precisely map the 5' end of the RBP binding site and uncovers direct and indirect targets of an RBP of interest. Importantly, ePRINT can also uncover RBPs that are differentially activated between cell fate transitions, including neural progenitor differentiation into neurons.

Transcriptome meta-analysis-based identification of hub transcription factors and RNA-binding proteins potentially orchestrating gene regulatory cascades and crosstalk in response to abiotic stresses in Arabidopsis thaliana

Deteriorating climatic conditions and increasing human population necessitate the development of robust plant varieties resistant to harsh environments. Manipulation of regulatory proteins such as transcription factors (TFs) and RNA-binding proteins (RBPs) would be a beneficial strategy in this regard. Further, understanding the complex interconnections between different classes of regulatory molecules would be essential for the identification of candidate genes/proteins for trait improvement. Most studies to date have analysed the roles of TFs or RBPs individually, in conferring stress resilience. However, it would be important to identify dominant/upstream TFs and RBPs inducing widespread transcriptomic alterations through other regulators (i.e., other TFs/RBPs targeted by the upstream regulators). To this end, the present study employed a transcriptome meta-analysis and computational approaches to obtain a comprehensive overview of regulatory interactions. This work identified dominant TFs and RBPs potentially influencing stress-mediated differential expression of other regulators, which could in turn influence gene expression, and consequently, physiological responses. Twenty transcriptomic studies [related to (i) UV radiation, (ii) wounding, (iii) salinity, (iv) cold, and (v) drought stresses in Arabidopsis thaliana] were analysed for differential gene expression, followed by the identification of differentially expressed TFs and RBPs. Subsequently, other TFs and RBPs which could be influencing these regulators were identified, and their interaction networks and hub nodes were analysed. As a result, an interacting module of Basic Leucine Zipper (bZIP) family TFs as well as Heterogeneous nuclear ribonucleoproteins (hnRNP) and Glycine-rich protein (GRP) family RBPs (among other TFs and RBPs) were shown to potentially influence the stress-induced differential expression of other TFs and RBPs under all the considered stress conditions. Some of the identified hub TFs and RBPs are known to be of major importance in orchestrating stress-induced transcriptomic changes influencing a variety of physiological processes from seed germination to senescence. This study highlighted the gene/protein candidates that could be considered for multiplexed genetic manipulation - a promising approach to develop robust, multi-stress-resilient plant varieties.

Dynamic subtype- and context-specific subcellular RNA regulation in growth cones of developing neurons of the cerebral cortex

Molecular mechanisms that cells employ to compartmentalize function via localization of function-specific RNA and translation are only partially elucidated. We investigate long-range projection neurons of the cerebral cortex as highly polarized exemplars to elucidate dynamic regulation of RNA localization, stability, and translation within growth cones (GCs), leading tips of growing axons. Comparison of GC-localized transcriptomes between two distinct subtypes of projection neurons- interhemispheric-callosal and corticothalamic- across developmental stages identifies both distinct and shared subcellular machinery, and intriguingly highlights enrichment of genes associated with neurodevelopmental and neuropsychiatric disorders. Developmental context-specific components of GC-localized transcriptomes identify known and novel potential regulators of distinct phases of circuit formation: long-distance growth, target area innervation, and synapse formation. Further, we investigate mechanisms by which transcripts are enriched and dynamically regulated in GCs, and identify GC-enriched motifs in 3' untranslated regions. As one example, we identify cytoplasmic adenylation element binding protein 4 (CPEB4), an RNA binding protein regulating localization and translation of mRNAs encoding molecular machinery important for axonal branching and complexity. We also identify RNA binding motif single stranded interacting protein 1 (RBMS1) as a dynamically expressed regulator of RNA stabilization that enables successful callosal circuit formation. Subtly aberrant associative and integrative cortical circuitry can profoundly affect cortical function, often causing neurodevelopmental and neuropsychiatric disorders. Elucidation of context-specific subcellular RNA regulation for GC- and soma-localized molecular controls over precise circuit development, maintenance, and function offers generalizable insights for other polarized cells, and might contribute substantially to understanding neurodevelopmental and behavioral-cognitive disorders and toward targeted therapeutics.

FAK and p130Cas modulate stiffness-mediated early transcription and cellular metabolism

Cellular metabolism is influenced by the stiffness of the extracellular matrix. Focal adhesion kinase (FAK) and its binding partner, p130Cas, transmit biomechanical signals about substrate stiffness to the cell to regulate a variety of cellular responses, but their roles in early transcriptional and metabolic responses remain largely unexplored. We cultured mouse embryonic fibroblasts with or without siRNA-mediated FAK or p130Cas knockdown and assessed the early transcriptional responses of these cells to placement on soft and stiff substrates by RNA sequencing and bioinformatics analyses. Exposure to the stiff ECM altered the expression of genes important for metabolic and biosynthetic processes, and these responses were influenced by knockdown of FAK and p130Cas. Our findings reveal that FAK-p130Cas signaling mechanotransduces ECM stiffness to early transcriptional changes that alter cellular metabolism and biosynthesis.

Polygenic risk for schizophrenia converges on alternative polyadenylation as molecular mechanism underlying synaptic impairment

Schizophrenia (SCZ) is a genetically heterogenous psychiatric disorder of highly polygenic nature. Correlative evidence from genetic studies indicate that the aggregated effects of distinct genetic risk factor combinations found in each patient converge onto common molecular mechanisms. To prove this on a functional level, we employed a reductionistic cellular model system for polygenic risk by differentiating induced pluripotent stem cells (iPSCs) from 104 individuals with high polygenic risk load and controls into cortical glutamatergic neurons (iNs). Multi-omics profiling identified widespread differences in alternative polyadenylation (APA) in the 3' untranslated region of many synaptic transcripts between iNs from SCZ patients and healthy donors. On the cellular level, 3'APA was associated with a reduction in synaptic density of iNs. Importantly, differential APA was largely conserved between postmortem human prefrontal cortex from SCZ patients and healthy donors, and strongly enriched for transcripts related to synapse biology. 3'APA was highly correlated with SCZ polygenic risk and affected genes were significantly enriched for SCZ associated common genetic variation. Integrative functional genomic analysis identified the RNA binding protein and SCZ GWAS risk gene PTBP2 as a critical trans-acting factor mediating 3'APA of synaptic genes in SCZ subjects. Functional characterization of PTBP2 in iNs confirmed its key role in 3'APA of synaptic transcripts and regulation of synapse density. Jointly, our findings show that the aggregated effects of polygenic risk converge on 3'APA as one common molecular mechanism that underlies synaptic impairments in SCZ.

Nucleocytoplasmic mRNA redistribution accompanies RNA binding protein mislocalization in ALS motor neurons and is restored by VCP ATPase inhibition

Amyotrophic lateral sclerosis (ALS) is characterized by nucleocytoplasmic mislocalization of the RNA-binding protein (RBP) TDP-43. However, emerging evidence suggests more widespread mRNA and protein mislocalization. Here, we employed nucleocytoplasmic fractionation, RNA sequencing, and mass spectrometry to investigate the localization of mRNA and protein in induced pluripotent stem cell-derived motor neurons (iPSMNs) from ALS patients with TARDBP and VCP mutations. ALS mutant iPSMNs exhibited extensive nucleocytoplasmic mRNA redistribution, RBP mislocalization, and splicing alterations. Mislocalized proteins exhibited a greater affinity for redistributed transcripts, suggesting a link between RBP mislocalization and mRNA redistribution. Notably, treatment with ML240, a VCP ATPase inhibitor, partially restored mRNA and protein localization in ALS mutant iPSMNs. ML240 induced changes in the VCP interactome and lysosomal localization and reduced oxidative stress and DNA damage. These findings emphasize the link between RBP mislocalization and mRNA redistribution in ALS motor neurons and highlight the therapeutic potential of VCP inhibition.

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.

Dynamically regulated transcription factors are encoded by highly unstable mRNAs in the Drosophila larval brain

The level of each RNA species depends on the balance between its rates of production and decay. Although previous studies have measured RNA decay across the genome in tissue culture and single-celled organisms, few experiments have been performed in intact complex tissues and organs. It is therefore unclear whether the determinants of RNA decay found in cultured cells are preserved in an intact tissue, and whether they differ between neighboring cell types and are regulated during development. To address these questions, we measured RNA synthesis and decay rates genome wide via metabolic labeling of whole cultured Drosophila larval brains using 4-thiouridine. Our analysis revealed that decay rates span a range of more than 100-fold, and that RNA stability is linked to gene function, with mRNAs encoding transcription factors being much less stable than mRNAs involved in core metabolic functions. Surprisingly, among transcription factor mRNAs there was a clear demarcation between more widely used transcription factors and those that are expressed only transiently during development. mRNAs encoding transient transcription factors are among the least stable in the brain. These mRNAs are characterized by epigenetic silencing in most cell types, as shown by their enrichment with the histone modification H3K27me3. Our data suggest the presence of an mRNA destabilizing mechanism targeted to these transiently expressed transcription factors to allow their levels to be regulated rapidly with high precision. Our study also demonstrates a general method for measuring mRNA transcription and decay rates in intact organs or tissues, offering insights into the role of mRNA stability in the regulation of complex developmental programs.

Local translation in microglial processes is required for efficient phagocytosis

Neurons, astrocytes and oligodendrocytes locally regulate protein translation within distal processes. Here, we tested whether there is regulated local translation within peripheral microglial processes (PeMPs) from mouse brain. We show that PeMPs contain ribosomes that engage in de novo protein synthesis, and these are associated with transcripts involved in pathogen defense, motility and phagocytosis. Using a live slice preparation, we further show that acute translation blockade impairs the formation of PeMP phagocytic cups, the localization of lysosomal proteins within them, and phagocytosis of apoptotic cells and pathogen-like particles. Finally, PeMPs severed from their somata exhibit and require de novo local protein synthesis to effectively surround pathogen-like particles. Collectively, these data argue for regulated local translation in PeMPs and indicate a need for new translation to support dynamic microglial functions.

Single-cell quantification of ribosome occupancy in early mouse development

Translation regulation is critical for early mammalian embryonic development¹. However, previous studies had been restricted to bulk measurements², precluding precise determination of translation regulation including allele-specific analyses. Here, to address this challenge, we developed a novel microfluidic isotachophoresis (ITP) approach, named RIBOsome profiling via ITP (Ribo-ITP), and characterized translation in single oocytes and embryos during early mouse development. We identified differential translation efficiency as a key mechanism regulating genes involved in centrosome organization and N⁶-methyladenosine modification of RNAs. Our high-coverage measurements enabled, to our knowledge, the first analysis of allele-specific ribosome engagement in early development. These led to the discovery of stage-specific differential engagement of zygotic RNAs with ribosomes and reduced translation efficiency of transcripts exhibiting allele-biased expression. By integrating our measurements with proteomics data, we discovered that ribosome occupancy in germinal vesicle-stage oocytes is the predominant determinant of protein abundance in the zygote. The Ribo-ITP approach will enable numerous applications by providing high-coverage and high-resolution ribosome occupancy measurements from ultra-low input samples including single cells.

New Functions of Intracellular LOXL2: Modulation of RNA-Binding Proteins

Lysyl oxidase-like 2 (LOXL2) was initially described as an extracellular enzyme involved in extracellular matrix remodeling. Nevertheless, numerous recent reports have implicated intracellular LOXL2 in a wide variety of processes that impact on gene transcription, development, differentiation, proliferation, migration, cell adhesion, and angiogenesis, suggesting multiple different functions for this protein. In addition, increasing knowledge about LOXL2 points to a role in several types of human cancer. Moreover, LOXL2 is able to induce the epithelial-to-mesenchymal transition (EMT) process-the first step in the metastatic cascade. To uncover the underlying mechanisms of the great variety of functions of intracellular LOXL2, we carried out an analysis of LOXL2's nuclear interactome. This study reveals the interaction of LOXL2 with numerous RNA-binding proteins (RBPs) involved in several aspects of RNA metabolism. Gene expression profile analysis of cells silenced for LOXL2, combined with in silico identification of RBPs' targets, points to six RBPs as candidates to be substrates of LOXL2's action, and that deserve a more mechanistic analysis in the future. The results presented here allow us to hypothesize novel LOXL2 functions that might help to comprehend its multifaceted role in the tumorigenic process.

Buffering of transcription rate by mRNA half-life is a conserved feature of Rett syndrome models

Transcriptional changes in Rett syndrome (RTT) are assumed to directly correlate with steady-state mRNA levels, but limited evidence in mice suggests that changes in transcription can be compensated by post-transcriptional regulation. We measure transcription rate and mRNA half-life changes in RTT patient neurons using RATEseq, and re-interpret nuclear and whole-cell RNAseq from Mecp2 mice. Genes are dysregulated by changing transcription rate or half-life and are buffered when both change. We utilized classifier models to predict the direction of transcription rate changes and find that combined frequencies of three dinucleotides are better predictors than CA and CG. MicroRNA and RNA-binding Protein (RBP) motifs are enriched in 3'UTRs of genes with half-life changes. Nuclear RBP motifs are enriched on buffered genes with increased transcription rate. We identify post-transcriptional mechanisms in humans and mice that alter half-life or buffer transcription rate changes when a transcriptional modulator gene is mutated in a neurodevelopmental disorder.

Integrated transcriptome and proteome analysis reveals posttranscriptional regulation of ribosomal genes in human brain organoids

During development of the human cerebral cortex, multipotent neural progenitors generate excitatory neurons and glial cells. Investigations of the transcriptome and epigenome have revealed important gene regulatory networks underlying this crucial developmental event. However, the posttranscriptional control of gene expression and protein abundance during human corticogenesis remains poorly understood. We addressed this issue by using human telencephalic brain organoids grown using a dual reporter cell line to isolate neural progenitors and neurons and performed cell class and developmental stage-specific transcriptome and proteome analysis. Integrating the two datasets revealed modules of gene expression during human corticogenesis. Investigation of one such module uncovered mTOR-mediated regulation of translation of the 5'TOP element-enriched translation machinery in early progenitor cells. We show that in early progenitors partial inhibition of the translation of ribosomal genes prevents precocious translation of differentiation markers. Overall, our multiomics approach proposes novel posttranscriptional regulatory mechanisms crucial for the fidelity of cortical development.

Computational analysis of crosstalk between transcriptional regulators and RNA-binding proteins suggests mutual regulation of polycomb proteins and SRSF1 influencing adult hippocampal neurogenesis

BACKGROUND

Adult hippocampal neurogenesis (AHN) is a clinically significant neural phenomenon. Understanding its molecular regulation would be important. In this regard, most studies have focused on transcriptional regulators (TRs), epigenetic modifiers, or non-coding RNAs. RNA-binding proteins (RBPs) have emerged as dominant molecular regulators. It would be significant to understand the potential cross-talk between RBPs and TRs, which could influence AHN.

METHODS

The present study employed computational analyses to identify RBPs and TRs regulating AHN, followed by the analysis of their interaction networks and detection of hub proteins. Next, the potential mutual regulation of hub TRs and RBPs was analyzed. Additionally, hippocampal genes differentially expressed upon exercise were analyzed for potential regulation by the identified TRs and RBPs.

RESULTS

105 TRs and 26 RBPs were found to influence AHN, which could also form interactive networks. Polycomb complex proteins were among the TR network hubs, while HNRNP and SRSF family members were among the hub RBPs. Further, the polycomb complex proteins and SRSF1 could have a mutual regulatory relationship, suggesting a cross-talk between epigenetic/transcriptional and post-transcriptional regulatory pathways. A number of exercise-induced hippocampal genes were also found to be potential targets of the identified TRs and RBPs.

CONCLUSION

SRSF1 may influence post-transcriptional stability, localization, and alternative splicing patterns of polycomb complex transcripts, and the polycomb proteins may in turn epigenetically influence the SRSF1. Further experimental validation of these regulatory loops/networks could provide novel insights into the molecular regulation of AHN, and unravel new targets for disease-treatment.

MAPKAPK2-centric transcriptome profiling reveals its major role in governing molecular crosstalk of IGFBP2, MUC4, and PRKAR2B during HNSCC pathogenesis

Transcriptome analysis of head and neck squamous cell carcinoma (HNSCC) has been pivotal to comprehending the convoluted biology of HNSCC tumors. MAPKAPK2 or MK2 is a critical modulator of the mRNA turnover of crucial genes involved in HNSCC progression. However, MK2-centric transcriptome profiles of tumors are not well known. This study delves into HNSCC progression with MK2 at the nexus to delineate the biological relevance and intricate crosstalk of MK2 in the tumor milieu. We performed next-generation sequencing-based transcriptome profiling of HNSCC cells and xenograft tumors to ascertain mRNA expression profiles in MK2-wild type and MK2-knockdown conditions. The findings were validated using gene expression assays, immunohistochemistry, and transcript turnover studies. Here, we identified a pool of crucial MK2-regulated candidate genes by annotation and differential gene expression analyses. Regulatory network and pathway enrichment revealed their significance and involvement in the HNSCC pathogenesis. Additionally, 3'-UTR-based filtering recognized important MK2-regulated downstream target genes and validated them by nCounter gene expression assays. Finally, immunohistochemistry and transcript stability studies revealed the putative role of MK2 in regulating the transcript turnover of IGFBP2, MUC4, and PRKAR2B in HNSCC. Conclusively, MK2-regulated candidate genes were identified in this study, and their plausible involvement in HNSCC pathogenesis was elucidated. These genes possess investigative values as targets for diagnosis and therapeutic interventions for HNSCC.

Genome-wide landscape of RNA-binding protein (RBP) networks as potential molecular regulators of psychiatric co-morbidities: a computational analysis

Background

Psychiatric disorders are a major burden on global health. These illnesses manifest as co-morbid conditions, further complicating the treatment. There is a limited understanding of the molecular and regulatory basis of psychiatric co-morbidities. The existing research in this regard has largely focused on epigenetic modulators, non-coding RNAs, and transcription factors. RNA-binding proteins (RBPs) functioning as multi-protein complexes are now known to be predominant controllers of multiple gene regulatory processes. However, their involvement in gene expression dysregulation in psychiatric co-morbidities is yet to be understood.

Results

Ten RBPs (QKI, ELAVL2, EIF2S1, SRSF3, IGF2BP2, EIF4B, SNRNP70, FMR1, DAZAP1, and MBNL1) were identified to be associated with psychiatric disorders such as schizophrenia, major depression, and bipolar disorders. Analysis of transcriptomic changes in response to individual depletion of these RBPs showed the potential influence of a large number of RBPs driving differential gene expression, suggesting functional cross-talk giving rise to multi-protein networks. Subsequent transcriptome analysis of post-mortem human brain samples from diseased and control individuals also suggested the involvement of ~ 100 RBPs influencing gene expression changes. These RBPs were found to regulate various processes including transcript splicing, mRNA transport, localization, stability, and translation. They were also found to form an extensive interactive network. Further, hnRNP, SRSF, and PCBP family RBPs, Matrin3, U2AF2, KHDRBS1, PTBP1, and also PABPN1 were found to be the hub proteins of the RBP network.

Conclusions

Extensive RBP networks involving a few hub proteins could result in transcriptome-wide dysregulation of post-transcriptional modifications, potentially driving multiple psychiatric disorders. Understanding the functional involvement of RBP networks in psychiatric disorders would provide insights into the molecular basis of psychiatric co-morbidities.

Roles of RNA-binding proteins in neurological disorders, COVID-19, and cancer

RNA-binding proteins (RBPs) have emerged as important players in multiple biological processes including transcription regulation, splicing, R-loop homeostasis, DNA rearrangement, miRNA function, biogenesis, and ribosome biogenesis. A large number of RBPs had already been identified by different approaches in various organisms and exhibited regulatory functions on RNAs' fate. RBPs can either directly or indirectly interact with their target RNAs or mRNAs to assume a key biological function whose outcome may trigger disease or normal biological events. They also exert distinct functions related to their canonical and non-canonical forms. This review summarizes the current understanding of a wide range of RBPs' functions and highlights their emerging roles in the regulation of diverse pathways, different physiological processes, and their molecular links with diseases. Various types of diseases, encompassing colorectal carcinoma, non-small cell lung carcinoma, amyotrophic lateral sclerosis, and Severe acute respiratory syndrome coronavirus 2, aberrantly express RBPs. We also highlight some recent advances in the field that could prompt the development of RBPs-based therapeutic interventions.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

Global Exploration of RNA-Binding Proteins in Exercise-Induced Adult Hippocampal Neurogenesis: A Transcriptome Meta-analysis and Computational Study

Voluntary physical exercise is a robust enhancer of adult hippocampal neurogenesis (AHN). A complete understanding of the molecular regulation of AHN is important in order to exploit the benefits of the process toward therapeutic approaches. Several factors such as epigenetic modifiers, non-coding RNAs, and transcription factors have been reported to regulate AHN. However, there is a limited understanding of the impact of RNA-binding proteins (RBPs) on exercise-mediated AHN, in spite of their well-documented significance in embryonic neurogenesis. The present study is the first global analysis to catalog the potential RBPs influencing exercise-mediated AHN. Here, a transcriptome meta-analysis was conducted to study exercise-mediated gene expression modulation in hippocampi of adult mice. Next, potential RBPs influencing transcriptome-wide expression changes via untranslated regions (UTRs) were identified. Among other RBPs, MATR3, Musashi, TIA1, and FXR2 (known critical modulators of neurogenesis) were found to potentially regulate gene expression patterns. Subsequently, binding sites of known neurogenesis-regulating RBPs were identified in the UTRs of AHN-associated genes modulated by exercise. Finally, a number of RBPs including RBFOX1, RBFOX3, and QKI (known regulators of neurogenesis) were found to be highly expressed in mouse hippocampal formation and also potentially interact with other RBPs, suggesting their combinatorial functioning in exercise-induced AHN. Thus, the present meta-analysis-based computational study identified several RBPs potentially important in exercise-induced AHN, which could form a foundation for further experiments to unravel RBP-mediated regulation of AHN.

Liver Transcriptome Dynamics During Hibernation Are Shaped by a Shifting Balance Between Transcription and RNA Stability

Hibernators dramatically lower metabolism to save energy while fasting for months. Prolonged fasting challenges metabolic homeostasis, yet small-bodied hibernators emerge each spring ready to resume all aspects of active life, including immediate reproduction. The liver is the body’s metabolic hub, processing and detoxifying macromolecules to provide essential fuels to brain, muscle and other organs throughout the body. Here we quantify changes in liver gene expression across several distinct physiological states of hibernation in 13-lined ground squirrels, using RNA-seq to measure the steady-state transcriptome and GRO-seq to measure transcription for the first time in a hibernator. Our data capture key timepoints in both the seasonal and torpor-arousal cycles of hibernation. Strong positive correlation between transcription and the transcriptome indicates that transcriptional control dominates the known seasonal reprogramming of metabolic gene expression in liver for hibernation. During the torpor-arousal cycle, however, discordance develops between transcription and the steady-state transcriptome by at least two mechanisms: 1) although not transcribed during torpor, some transcripts are unusually stable across the torpor bout; and 2) unexpectedly, on some genes, our data suggest continuing, slow elongation with a failure to terminate transcription across the torpor bout. While the steady-state RNAs corresponding to these read through transcripts did not increase during torpor, they did increase shortly after rewarming despite their simultaneously low transcription. Both of these mechanisms would assure the immediate availability of functional transcripts upon rewarming. Integration of transcriptional, post-transcriptional and RNA stability control mechanisms, all demonstrated in these data, likely initiate a serial gene expression program across the short euthermic period that restores the tissue and prepares the animal for the next bout of torpor.

Dynamic RNA Regulation in the Brain Underlies Physiological Plasticity in a Hibernating Mammal

Hibernation is a physiological and behavioral phenotype that minimizes energy expenditure. Hibernators cycle between profound depression and rapid hyperactivation of multiple physiological processes, challenging our concept of mammalian homeostasis. How the hibernator orchestrates and survives these extremes while maintaining cell to organismal viability is unknown. Here, we enhance the genome integrity and annotation of a model hibernator, the 13-lined ground squirrel. Our new assembly brings this genome to near chromosome-level contiguity and adds thousands of previously unannotated genes. These new genomic resources were used to identify 6,505 hibernation-related, differentially-expressed and processed transcripts using RNA-seq data from three brain regions in animals whose physiological status was precisely defined using body temperature telemetry. A software tool, squirrelBox, was developed to foster further data analyses and visualization. SquirrelBox includes a comprehensive toolset for rapid visualization of gene level and cluster group dynamics, sequence scanning of k-mer and domains, and interactive exploration of gene lists. Using these new tools and data, we deconvolute seasonal from temperature-dependent effects on the brain transcriptome during hibernation for the first time, highlighting the importance of carefully timed samples for studies of differential gene expression in hibernation. The identified genes include a regulatory network of RNA binding proteins that are dynamic in hibernation along with the composition of the RNA pool. In addition to passive effects of temperature, we provide evidence for regulated transcription and RNA turnover during hibernation. Significant alternative splicing, largely temperature dependent, also occurs during hibernation. These findings form a crucial first step and provide a roadmap for future work toward defining novel mechanisms of tissue protection and metabolic depression that may 1 day be applied toward improving human health.