Search

RNA sequencing is a powerful tool to analyze gene expression transcriptome wide. However, RNA sequencing in general and especially the recently developed methods of long read RNA sequencing are still low-throughput and cost-intensive. Here, one important design choice is to concentrate the sequencing capacity on specific parts of the transcriptome. Especially, abundant transcripts as ribosomal RNAs may dominate the available sequencing space, if not removed prior to sequencing. Several methods exist to reduce ribosomal RNA read numbers: either based on enrichment of the relevant fraction (polyA+ RNA) or depletion, respectively degradation of ribosomal RNAs. Furthermore, commercial kits are available to deplete globin transcripts from blood samples. However, so far, no solution exists to deal with other tissue-specific highly abundant transcripts. This is especially of interest in the heart and other muscle derived samples, where reads originating from mitochondrial RNAs make up to 30% of reads in polyA+ selected libraries and around 70% in single cell sequencing experiments. We present a simple method to diminish sequencing of mitochondrial RNAs in Oxford Nanopore direct RNA sequencing libraries by RNase H based clipping of the polyA tail. We show that mt-clipping enables enhanced detection of cytoplasmic mRNAs, among them genes involved in heart development and pathogenesis. Mt-clipping may be applied as well to other sequencing protocols that are based on oligo(dT) priming and can be easily adapted to other tissue-specific high-abundant transcripts.

The exon junction complex (EJC) is an RNA-binding multi-protein complex with critical functions in post-transcriptional gene regulation. It is deposited on the mRNA during splicing and regulates diverse processes including pre-mRNA splicing and nonsense-mediated mRNA decay (NMD) via various interacting proteins. The peripheral EJC-binding protein RNPS1 was reported to serve two insufficiently characterized functions: suppressing mis-splicing of cryptic splice sites and activating NMD in the cytoplasm. The analysis of transcriptome-wide effects of EJC and RNPS1 knockdowns in different human cell lines supports the conclusion that RNPS1 can moderately influence NMD activity, but is not a globally essential NMD factor. However, numerous aberrant splicing events strongly suggest that the main function of RNPS1 is splicing regulation. Rescue analyses revealed that the RRM and C-terminal domain of RNPS1 both contribute partially to regulate RNPS1-dependent splicing events. We defined the RNPS1 core interactome using complementary immunoprecipitations and proximity labeling, which identified interactions with splicing-regulatory factors that are dependent on the C-terminus or the RRM domain of RNPS1. Thus, RNPS1 emerges as a multifunctional splicing regulator that promotes correct and efficient splicing of different vulnerable splicing events via the formation of diverse splicing-promoting complexes.

Clinically translatable large animal models have become indispensable for cardiovascular research, clinically relevant proof of concept studies and for novel therapeutic interventions. In particular, the pig has emerged as an essential cardiovascular disease model, because its heart, circulatory system, and blood supply are anatomically and functionally similar to that of humans. Currently, molecular and omics-based studies in the pig are hampered by the incompleteness of the genome and the lack of diversity of the corresponding transcriptome annotation. Here, we employed Nanopore long-read sequencing and in-depth proteomics on top of Illumina RNA-seq to enhance the pig cardiac transcriptome annotation. We assembled 15,926 transcripts, stratified into coding and non-coding, and validated our results by complementary mass spectrometry. A manual review of several gene loci, which are associated with cardiac function, corroborated the utility of our enhanced annotation. All our data are available for download and are provided as tracks for integration in genome browsers. We deem this resource as highly valuable for molecular research in an increasingly relevant large animal model.

Chemical modifications are found on almost all RNAs and affect their coding and noncoding functions. The identification of m6A on mRNA and its important role in gene regulation stimulated the field to investigate whether additional modifications are present on mRNAs. Indeed, modifications including m1A, m5C, m7G, 2′-OMe, and Ψ were detected. However, since their abundances are low and tools used for their corroboration are often not well characterized, their physiological relevance remains largely elusive. Antibodies targeting modified nucleotides are often used but have limitations such as low affinity or specificity. Moreover, they are not always well characterized and due to the low abundance of the modification, particularly on mRNAs, generated data sets might resemble noise rather than specific modification patterns. Therefore, it is critical that the affinity and specificity is rigorously tested using complementary approaches. Here, we provide an experimental toolbox that allows for testing antibody performance prior to their use.

The exon junction complex (EJC) is an essential constituent and regulator of spliced messenger ribonucleoprotein particles (mRNPs) in metazoans. As a core component of the EJC, CASC3 was described to be pivotal for EJC-dependent nuclear and cytoplasmic processes. However, recent evidence suggests that CASC3 functions differently from other EJC core proteins. Here, we have established human CASC3 knockout cell lines to elucidate the cellular role of CASC3. In the knockout cells, overall EJC composition and EJC-dependent splicing are unchanged. A transcriptome-wide analysis reveals that hundreds of mRNA isoforms targeted by nonsense-mediated decay (NMD) are upregulated. Mechanistically, recruiting CASC3 to reporter mRNAs by direct tethering or via binding to the EJC stimulates mRNA decay and endonucleolytic cleavage at the termination codon. Building on existing EJC-NMD models, we propose that CASC3 equips the EJC with the persisting ability to communicate with the NMD machinery in the cytoplasm. Collectively, our results characterize CASC3 as a peripheral EJC protein that tailors the transcriptome by promoting the degradation of EJC-dependent NMD substrates.

Funder: Wellcome Trust, Funder: Medical Research Council, Pausing of RNA polymerase II (Pol II) close to promoters is a common regulatory step in RNA synthesis, and is coordinated by a ribonucleoprotein complex scaffolded by the noncoding RNA RN7SK. The function of RN7SK-regulated gene transcription in adult tissue homoeostasis is currently unknown. Here, we deplete RN7SK during mouse and human epidermal stem cell differentiation. Unexpectedly, loss of this small nuclear RNA specifically reduces transcription of numerous cell cycle regulators leading to cell cycle exit and differentiation. Mechanistically, we show that RN7SK is required for efficient transcription of highly expressed gene pairs with bidirectional promoters, which in the epidermis co-regulated cell cycle and chromosome organization. The reduction in transcription involves impaired splicing and RNA decay, but occurs in the absence of chromatin remodelling at promoters and putative enhancers. Thus, RN7SK is directly required for efficient Pol II transcription of highly transcribed bidirectional gene pairs, and thereby exerts tissue-specific functions, such as maintaining a cycling cell population in the epidermis.

The mechanistic target of rapamycin (mTOR) promotes pathological remodeling in the heart by activating ribosomal biogenesis and mRNA translation. Inhibition of mTOR in cardiomyocytes is protective; however, a detailed role of mTOR in translational regulation of specific mRNA networks in the diseased heart is unknown. We performed cardiomyocyte genome‐wide sequencing to define mTOR‐dependent gene expression control at the level of mRNA translation. We identify the muscle‐specific protein Cullin‐associated NEDD8‐dissociated protein 2 (Cand2) as a translationally upregulated gene, dependent on the activity of mTOR. Deletion of Cand2 protects the myocardium against pathological remodeling. Mechanistically, we show that Cand2 links mTOR signaling to pathological cell growth by increasing Grk5 protein expression. Our data suggest that cell‐type‐specific targeting of mTOR might have therapeutic value against pathological cardiac remodeling., Genome‐wide translational profiling identifies mTORC1‐dependent genes in cardiomyocytes in response to neurohumoral stimulation. Expression of the muscle‐specific gene Cand2 is controlled by mTORC1 and Cand2 regulates cardiac function and pathological hypertrophy.

Background Circular RNAs (circRNAs) are a special group of non-coding RNAs formed by back-splicing. Mostly cytosolic in eukaryotic cells, circRNAs are particularly enriched and conserved in neurons and originate from protein-coding genes to function as global regulators of gene expression. Recent studies showed that circRNAs partake in brain physiology and pathology, contributing to neurological disorders, such as myotonic dystrophy type 1, Parkinson’s and Alzheimer’s Diseases. Here, we identified the first ever known brain enriched RNA circle originating from the Huntington’s Disease gene HTT (CircHTT: 484 nt, Ex 2-6, chr4:3088665-3109150) which is conserved in mouse and minipig. Aims and Methods To functionally characterize circHTT and its implication for HD pathogenesis, we first studied its expression pattern in human and mouse body districts and in iPS-derived neuronal cell lines with different CAG repeats. Then, overexpression and down-regulation of the circle were used to determine the effects on transcription of the HD gene and translation of wild-type and mutant protein. Results CircHTT expression increases significantly with increasing CAG repeats in terminally differentiated cortical neurons. Furthermore, circHtt is significantly more expressed in brain districts of Q111 and zQ175 knock-in mice. These findings indicate that the expression of circHTT/Htt occurs in a CAG repeat dependent manner in neuronal cells, typical hallmark of HD pathology. Strikingly, overexpression of the circular RNA in human HEK293T, PC3, mouse STHdh Q7/7, Q7/111 and Q111/111 cell lines consistently increase wild-type huntingtin, while decreasing mutant huntingtin protein, with no alteration of the HTT/Htt transcript level. Conclusions In conclusion, we identified a brain enriched RNA circle originating from the HD gene locus that is sensitive to the HD mutation and may modulate huntingtin expression. Our observations might pave the way to new trials of therapeutic intervention.

Synaptic scaling is a form of homeostatic plasticity which allows neurons to adjust their action potential firing rate in response to chronic alterations in neural activity. Synaptic scaling requires profound changes in gene expression, but the relative contribution of local and cell‐wide mechanisms is controversial. Here we perform a comprehensive multi‐omics characterization of the somatic and process compartments of primary rat hippocampal neurons during synaptic scaling. We uncover both highly compartment‐specific and correlating changes in the neuronal transcriptome and proteome. Whereas downregulation of crucial regulators of neuronal excitability occurs primarily in the somatic compartment, structural components of excitatory postsynapses are mostly downregulated in processes. Local inhibition of protein synthesis in processes during scaling is confirmed for candidate synaptic proteins. Motif analysis further suggests an important role for trans‐acting post‐transcriptional regulators, including RNA‐binding proteins and microRNAs, in the local regulation of the corresponding mRNAs. Altogether, our study indicates that, during synaptic scaling, compartmentalized gene expression changes might co‐exist with neuron‐wide mechanisms to allow synaptic computation and homeostasis., Downscaling of excitatory synapses in response to high activity involves compartment‐specific changes in the neuron's transcriptome and proteome. This study supports the view that local gene regulation is important in both Hebbian and homeostatic forms of synaptic plasticity.

Coronavirus disease 2019 (COVID-19) spawned a global health crisis in late 2019 and is caused by the novel coronavirus SARS-CoV-2. SARS-CoV-2 infection can lead to elevated markers of endothelial dysfunction associated with higher risk of mortality. It is unclear whether endothelial dysfunction is caused by direct infection of endothelial cells or is mainly secondary to inflammation. Here, we investigate whether different types of endothelial cells are susceptible to SARS-CoV-2. Human endothelial cells from different vascular beds including umbilical vein endothelial cells, coronary artery endothelial cells (HCAEC), cardiac and lung microvascular endothelial cells, or pulmonary arterial cells were inoculated in vitro with SARS-CoV-2. Viral spike protein was only detected in HCAECs after SARS-CoV-2 infection but not in the other endothelial cells tested. Consistently, only HCAEC expressed the SARS-CoV-2 receptor angiotensin-converting enzyme 2 (ACE2), required for virus infection. Infection with the SARS-CoV-2 variants B.1.1.7, B.1.351, and P.2 resulted in significantly higher levels of viral spike protein. Despite this, no intracellular double-stranded viral RNA was detected and the supernatant did not contain infectious virus. Analysis of the cellular distribution of the spike protein revealed that it co-localized with endosomal calnexin. SARS-CoV-2 infection did induce the ER stress gene EDEM1, which is responsible for clearance of misfolded proteins from the ER. Whereas the wild type of SARS-CoV-2 did not induce cytotoxic or pro-inflammatory effects, the variant B.1.1.7 reduced the HCAEC cell number. Of the different tested endothelial cells, HCAECs showed highest viral uptake but did not promote virus replication. Effects on cell number were only observed after infection with the variant B.1.1.7, suggesting that endothelial protection may be particularly important in patients infected with this variant. Supplementary Information The online version contains supplementary material available at 10.1007/s00395-021-00882-8.

The current ecosystem of single-cell RNA-seq platforms is rapidly expanding, but robust solutions for single-cell and single-molecule full-length RNA sequencing are virtually absent. A high-throughput solution that covers all aspects is necessary to study the complex life of mRNA on the single-cell level. The Nanopore platform offers long read sequencing and can be integrated with the popular single-cell sequencing method on the 10× Chromium platform. However, the high error-rate of Nanopore reads poses a challenge in downstream processing (e.g., for cell barcode assignment). We propose a solution to this particular problem by using a hybrid sequencing approach on Nanopore and Illumina platforms. Our software ScNapBar enables cell barcode assignment with high accuracy, especially if sequencing saturation is low. ScNapBar uses unique molecular identifier (UMI) or Naïve Bayes probabilistic approaches in the barcode assignment, depending on the available Illumina sequencing depth. We have benchmarked the two approaches on simulated and real Nanopore data sets. We further applied ScNapBar to pools of cells with an active or a silenced nonsense-mediated RNA decay pathway. Our Nanopore read assignment distinguishes the respective cell populations and reveals characteristic nonsense-mediated mRNA decay events depending on cell status.

Nucleic acids research : NAR 49(16), e92 (2021). doi:10.1093/nar/gkab485, Published by Oxford Univ. Press, Oxford

In the adult heart, the epicardium becomes activated after injury, contributing to cardiac healing by secretion of paracrine factors. Here, we analyzed by single-cell RNA sequencing combined with RNA in situ hybridization and lineage tracing of Wilms tumor protein 1-positive (WT1+) cells, the cellular composition, location, and hierarchy of epicardial stromal cells (EpiSC) in comparison to activated myocardial fibroblasts/stromal cells in infarcted mouse hearts. We identified 11 transcriptionally distinct EpiSC populations, which can be classified into three groups, each containing a cluster of proliferating cells. Two groups expressed cardiac specification markers and sarcomeric proteins suggestive of cardiomyogenic potential. Transcripts of hypoxia-inducible factor (HIF)-1α and HIF-responsive genes were enriched in EpiSC consistent with an epicardial hypoxic niche. Expression of paracrine factors was not limited to WT1+ cells but was a general feature of activated cardiac stromal cells. Our findings provide the cellular framework by which myocardial ischemia may trigger in EpiSC the formation of cardioprotective/regenerative responses.

Catecholamines stimulate the first step of lipolysis by PKA-dependent release of the lipid droplet-associated protein ABHD5 from perilipin to co-activate the lipase ATGL. Here, we unmask a yet unrecognized proteolytic and cardioprotective function of ABHD5. ABHD5 acts in vivo and in vitro as a serine protease cleaving HDAC4. Through the production of an N-terminal polypeptide of HDAC4 (HDAC4-NT), ABHD5 inhibits MEF2-dependent gene expression and thereby controls glucose handling. ABHD5-deficiency leads to neutral lipid storage disease in mice. Cardiac-specific gene therapy of HDAC4-NT does not protect from intra-cardiomyocyte lipid accumulation but strikingly from heart failure, thereby challenging the concept of lipotoxicity-induced heart failure. ABHD5 levels are reduced in failing human hearts and murine transgenic ABHD5 expression protects from pressure-overload induced heart failure. These findings represent a conceptual advance by connecting lipid with glucose metabolism through HDAC4 proteolysis and enable new translational approaches to treat cardiometabolic disease.

Diseases of the renal filtration unit—the glomerulus—are the most common cause of chronic kidney disease. Podocytes are the pivotal cell type for the function of this filter and focal-segmental glomerulosclerosis (FSGS) is a classic example of a podocytopathy leading to proteinuria and glomerular scarring. Currently, no targeted treatment of FSGS is available. This lack of therapeutic strategies is explained by a limited understanding of the defects in podocyte cell biology leading to FSGS. To date, most studies in the field have focused on protein-coding genes and their gene products. However, more than 80% of all transcripts produced by mammalian cells are actually non-coding. Here, long non-coding RNAs (lncRNAs) are a relatively novel class of transcripts and have not been systematically studied in FSGS to date. The appropriate tools to facilitate lncRNA research for the renal scientific community are urgently required due to a row of challenges compared to classical analysis pipelines optimized for coding RNA expression analysis. Here, we present the bioinformatic pipeline CALINCA as a solution for this problem. CALINCA automatically analyzes datasets from murine FSGS models and quantifies both annotated and de novo assembled lncRNAs. In addition, the tool provides in-depth information on podocyte specificity of these lncRNAs, as well as evolutionary conservation and expression in human datasets making this pipeline a crucial basis to lncRNA studies in FSGS.

In the adult heart, the epicardium becomes activated after injury, contributing to cardiac healing by secretion of paracrine factors. Here we analyzed by single-cell RNA sequencing combined with RNA in situ hybridization and lineage tracing of WT1+ cells the cellular composition, location, and hierarchy of epicardial stromal cells (EpiSC) in comparison to activated myocardial fibroblasts/stromal cells in infarcted mouse hearts. We identified 11 transcriptionally distinct EpiSC populations, that can be classified in three groups each containing a cluster of proliferating cells. Two groups expressed cardiac specification makers and sarcomeric proteins suggestive of cardiomyogenic potential. Transcripts of HIF-1α and HIF-responsive genes were enriched in EpiSC consistent with an epicardial hypoxic niche. Expression of paracrine factors was not limited to WT1+ cells but was a general feature of activated cardiac stromal cells. Our findings provide the cellular framework by which myocardial ischemia may trigger in EpiSC the formation of cardioprotective/regenerative responses.

N6‐methyladenosine (m6A) regulates a variety of physiological processes through modulation of RNA metabolism. This modification is particularly enriched in the nervous system of several species, and its dysregulation has been associated with neurodevelopmental defects and neural dysfunctions. In Drosophila, loss of m6A alters fly behavior, albeit the underlying molecular mechanism and the role of m6A during nervous system development have remained elusive. Here we find that impairment of the m6A pathway leads to axonal overgrowth and misguidance at larval neuromuscular junctions as well as in the adult mushroom bodies. We identify Ythdf as the main m6A reader in the nervous system, being required to limit axonal growth. Mechanistically, we show that the m6A reader Ythdf directly interacts with Fmr1, the fly homolog of Fragile X mental retardation RNA binding protein (FMRP), to inhibit the translation of key transcripts involved in axonal growth regulation. Altogether, this study demonstrates that the m6A pathway controls development of the nervous system and modulates Fmr1 target transcript selection., Proper neuromuscular junction formation in flies depends on a key mRNA modification guiding translational repression by the Fragile X mental retardation (FMRP) RNA‐binding protein.

Similar to DNA modifications, an epitranscriptome of mRNA modifications emerged over the last decade. Specifically, the discovery of the N6-methyladenosine (m6A) as reversible posttranscriptional modifications in mRNAs renewed interest in mRNA modifications and their role in health and disease. First m6A mRNA methylomes of human cardiac tissue from healthy patients and heart failure patients and of animal model heart samples were reported and shed point toward an important role of epitranscriptomic control on cell growth. More studies are clearly needed to fully understand this novel stress–response mechanism in the heart for maintaining normal cardiac function. Moreover, the role of mRNA modifications on physiological growth in the myocardium is unknown. In this chapter, we review the current state of the art of epitranscriptomic research in the heart and highlight innovative experimental and computational approaches, which will definitely drive this field forward.

Splicing varies across brain regions, but the single-cell resolution of regional variation is unclear. We present a single-cell investigation of differential isoform expression (DIE) between brain regions using single-cell long-read sequencing in mouse hippocampus and prefrontal cortex in 45 cell types at postnatal day 7 (www.isoformAtlas.com). Isoform tests for DIE show better performance than exon tests. We detect hundreds of DIE events traceable to cell types, often corresponding to functionally distinct protein isoforms. Mostly, one cell type is responsible for brain-region specific DIE. However, for fewer genes, multiple cell types influence DIE. Thus, regional identity can, although rarely, override cell-type specificity. Cell types indigenous to one anatomic structure display distinctive DIE, e.g. the choroid plexus epithelium manifests distinct transcription-start-site usage. Spatial transcriptomics and long-read sequencing yield a spatially resolved splicing map. Our methods quantify isoform expression with cell-type and spatial resolution and it contributes to further our understanding of how the brain integrates molecular and cellular complexity., Alternative RNA splicing varies across the brain. Its mapping at single cell resolution is unclear. Here, the authors provide a spatial and single-cell splicing atlas reporting brain region- and cell type-specific expression of different isoforms in the postnatal mouse brain.

The mechanistic target of rapamycin (mTOR) is a key regulator of pathological remodeling in the heart by activating ribosomal biogenesis and mRNA translation. Inhibition of mTOR in cardiomyocytes is protective, however, a detailed role of mTOR in translational regulation of specific mRNA networks in the diseased heart is largely unknown. A cardiomyocyte genome-wide sequencing approach was used to define mTOR-dependent post-transcriptional gene expression control at the level of mRNA translation. This approach identified the muscle-specific protein Cullin-associated NEDD8-dissociated protein 2 (Cand2) as a translationally upregulated gene, dependent on the activity of mTOR. Deletion of Cand2 protects the myocardium against pathological remodeling. Mechanistically, we found that Cand2 links mTOR signaling to pathological cell growth by increasing Grk5 protein expression. Our data suggest that cell-type-specific targeting of mTOR might have therapeutic value for adverse pathological cardiac remodeling.

Synaptic scaling is a form of homeostatic plasticity which allows neurons to adjust their action potential firing rate in response to chronic alterations in neural activity. Synaptic scaling requires profound changes in gene expression, but the relative contribution of local and cell-wide mechanisms is controversial. Here we performed a comprehensive multi-omics characterization of the somatic and process compartments of primary rat hippocampal neurons during synaptic scaling. Thereby, we uncovered both highly compartment-specific and correlated changes in the neuronal transcriptome and proteome. Whereas downregulation of crucial regulators of neuronal excitability occurred primarily in the somatic compartment, structural components of excitatory postsynapses were mostly downregulated in processes. Motif analysis further suggests an important role for trans-acting post-transcriptional regulators, including RNA-binding proteins and microRNAs, in the local regulation of the corresponding mRNAs. Altogether, our study indicates that compartmentalized gene expression changes are widespread in synaptic scaling and might co-exist with neuron-wide mechanisms to allow synaptic computation and homeostasis.

The activation induced cytidine deaminase/apolipoprotein B editing complex (AID/APOBEC) family comprises several nucleic acid editors with roles ranging from antibody diversification to mRNA editing. APOBEC2, an evolutionarily conserved member of this family, has neither an established substrate nor a mechanism of action, however genetic evidence suggests functional relevance in tissues such as muscle. Here, we demonstrate that in muscle, APOBEC2 does not have any of the attributed molecular functions of the AID/APOBEC family, such as RNA editing, DNA demethylation, or DNA mutation. Instead, we show that APOBEC2 occupies chromatin at promoter regions of certain genes, whose expression is repressed during muscle cell differentiation. We further demonstrate that APOBEC2 on one hand binds promoter region DNA directly and in a sequence specific fashion, while on the other it interacts with HDAC transcriptional corepressor complexes. Therefore, APOBEC2, by actively repressing the expression of non-myogenesis pathway genes, plays a key role in enforcing the proper establishment of muscle cell fate.

Our understanding of the transition from physiological to pathological cardiac hypertrophy remains elusive and largely based on reductionist hypotheses. Here, we profiled the translatomes of 15 mouse hearts to provide a molecular blueprint of altered gene networks in early cardiac remodeling. Using co-expression analysis, we reveal how sub-networks are orchestrated into functional modules associated with pathological phenotypes. We show how transcriptome networks are only partially reproducible at the translatome level. We find unappreciated hub genes and genes in the transcriptional network that were rewired in the translational network, and associated with semantically different subsets of enriched functional terms, providing novel insights into the complexity of the organization of in vivo cardiac regulatory networks.

N6-methyladenosine (m6A) regulates a variety of physiological processes through modulation of RNA metabolism. The modification is particularly enriched in the nervous system of several species, and its dysregulation has been associated with neurodevelopmental defects and neural dysfunctions. InDrosophila, loss of m6A alters fly behavior albeit the underlying mechanism and the role of m6A during nervous system development have remained elusive. Here we find that impairment of the m6A pathway leads to axonal overgrowth and misguidance at larval neuromuscular junctions as well as in the adult mushroom bodies. We identify Ythdf as the main m6A reader in the nervous system being required for limiting axonal growth. Mechanistically, we show that Ythdf directly interacts with Fragile X mental retardation protein to inhibit the translation of key transcripts involved in axonal growth regulation. Altogether, this study demonstrates that the m6A pathway controls development of the nervous system by modulating Fmr1 target selection.

Obesity and type 2 diabetes mellitus are global emergencies and long noncoding RNAs (lncRNAs) are regulatory transcripts with elusive functions in metabolism. Here we show that a high fraction of lncRNAs, but not protein-coding mRNAs, are repressed during diet-induced obesity (DIO) and refeeding, whilst nutrient deprivation induced lncRNAs in mouse liver. Similarly, lncRNAs are lost in diabetic humans. LncRNA promoter analyses, global cistrome and gain-of-function analyses confirm that increased MAFG signaling during DIO curbs lncRNA expression. Silencing Mafg in mouse hepatocytes and obese mice elicits a fasting-like gene expression profile, improves glucose metabolism, de-represses lncRNAs and impairs mammalian target of rapamycin (mTOR) activation. We find that obesity-repressed LincIRS2 is controlled by MAFG and observe that genetic and RNAi-mediated LincIRS2 loss causes elevated blood glucose, insulin resistance and aberrant glucose output in lean mice. Taken together, we identify a MAFG-lncRNA axis controlling hepatic glucose metabolism in health and metabolic disease., Despite widespread transcription of LncRNA in mammalian systems, their contribution to metabolic homeostasis at the cellular and tissue level remains elusive. Here Pradas-Juni et al. describe a transcription factor–LncRNA pathway that couples hepatocyte nutrient sensing to regulation of glucose metabolism in mice.

For decades, cardiovascular disease (CVD) has been the leading cause of death throughout most developed countries. Several studies relate RNA splicing, and more recently also circular RNAs (circRNAs), to CVD. CircRNAs originate from linear transcripts and have been shown to exhibit tissue-specific expression profiles. Here, we present an in-depth analysis of sequence, structure, modification, and cardiac circRNA interactions. We used human induced pluripotent stem cell-derived cardiac myocytes (hiPSC-CMs), human healthy and diseased (ischemic cardiomyopathy, dilated cardiomyopathy) cardiac tissue, and human umbilical vein endothelial cells (HUVECs) to profile circRNAs. We identified shared circRNAs across all samples, as well as model-specific circRNA signatures. Based on these circRNAs, we identified 63 positionally conserved and expressed circRNAs in human, pig, and mouse hearts. Furthermore, we found that the sequence of circRNAs can deviate from the sequence derived from the genome sequence, an important factor in assessing potential functions. Integration of additional data yielded evidence for m6A-methylation of circRNAs, potentially linked to translation, as well as, circRNAs overlapping with potential Argonaute 2 binding sites, indicating potential association with the RISC complex. Moreover, we describe, for the first time in cardiac model systems, a sub class of circRNAs containing the start codon of their primary transcript (AUG circRNAs) and observe an enrichment for m6A-methylation for AUG circRNAs.

Until recently considered as rare, circular RNAs (circRNAs) are emerging as important regulators of gene expression. They are ubiquitously expressed and represent a novel branch of the family of non-coding RNAs. Recent investigations showed that circRNAs are regulated in the cardiovascular system and participate in its physiological and pathological development. In this review article, we will provide an overview of the role of circRNAs in cardiovascular health and disease. After a description of the biogenesis of circRNAs, we will summarize what is known of the expression, regulation and function of circRNAs in the cardiovascular system. We will then address some technical aspects of circRNAs research, discussing how artificial intelligence may aid in circRNAs research. Finally, the potential of circRNAs as biomarkers of cardiovascular disease will be addressed and directions for future research will be proposed. Keywords: Non-coding RNAs, Circular RNAs, Cardiovascular system, Cardiovascular disease, Artificial intelligence, Biomarker

Synaptic downscaling is a homeostatic mechanism that allows neurons to reduce firing rates during chronically elevated network activity. Although synaptic downscaling is important in neural circuit development and epilepsy, the underlying mechanisms are poorly described. We performed small RNA profiling in picrotoxin (PTX)‐treated hippocampal neurons, a model of synaptic downscaling. Thereby, we identified eight microRNAs (miRNAs) that were increased in response to PTX, including miR‐129‐5p, whose inhibition blocked synaptic downscaling in vitro and reduced epileptic seizure severity in vivo . Using transcriptome, proteome, and bioinformatic analysis, we identified the calcium pump Atp2b4 and doublecortin (Dcx) as miR‐129‐5p targets. Restoring Atp2b4 and Dcx expression was sufficient to prevent synaptic downscaling in PTX‐treated neurons. Furthermore, we characterized a functional crosstalk between miR‐129‐5p and the RNA‐binding protein (RBP) Rbfox1. In the absence of PTX, Rbfox1 promoted the expression of Atp2b4 and Dcx. Upon PTX treatment, Rbfox1 expression was downregulated by miR‐129‐5p, thereby allowing the repression of Atp2b4 and Dcx. We therefore identified a novel activity‐dependent miRNA/RBP crosstalk during synaptic scaling, with potential implications for neural network homeostasis and epileptogenesis.

Promoter methylation status of O-6-methylguanine-DNA methyltransferase (MGMT), a DNA repair enzyme, is a critical biomarker in glioblastoma (GBM), as treatment decisions and clinical trial inclusion rely on its accurate assessment. However, interpretation of results is complicated by poor interassay reproducibility as well as a weak correlation between methylation status and expression levels of MGMT. This study systematically investigates the influence of tumor purity on tissue subjected to MGMT analysis. A quantitative, allele-specific real-time PCR (qAS-PCR) assay was developed to determine genotype and mutant allele frequency of telomerase promoter (pTERT) mutations as a direct measure of tumor purity. We studied tumor purity, pTERT mutation by Sanger sequencing, MGMT methylation by pyrosequencing, IDH1 mutation status, and clinical parameters in a cohort of high-grade gliomas (n = 97). The qAS-PCR reliably predicted pTERT genotype and tumor purity compared with independent methods. Tumor purity positively and significantly correlated with the extent of methylation in MGMT methylated GBMs. Extent of MGMT methylation differed significantly with respect to pTERT mutation hotspot (C228T vs. C250T). Interestingly, frontal lobe tumors showed greater tumor purity than those in other locations. Above all, tumor purity was identified as an independent prognostic factor in GBM. In conclusion, we determined mutual associations of tumor purity with MGMT methylation and pTERT mutations and found that the extent of MGMT methylation reflects tumor purity. In turn, tumor purity is prognostic in IDH1 wild-type GBM. Implications: Tumor purity is an independent prognostic marker in glioblastoma and is associated with the extent of MGMT methylation. Mol Cancer Res; 15(5); 532–40. ©2017 AACR.

Ribosome profiling via high-throughput sequencing (ribo-seq) is a promising new technique for characterizing the occupancy of ribosomes on messenger RNA (mRNA) at base-pair resolution. The ribosome is responsible for translating mRNA into proteins, so information about its occupancy offers a detailed view of ribosome density and position which could be used to discover new translated open reading frames (ORFs), among other things. In this work, we propose Rp-Bp, an unsupervised Bayesian approach to predict translated ORFs from ribosome profiles. We use state-of-the-art Markov chain Monte Carlo techniques to estimate posterior distributions of the likelihood of translation of each ORF. Hence, an important feature of Rp-Bp is its ability to incorporate and propagate uncertainty in the prediction process. A second novel contribution is automatic Bayesian selection of read lengths and ribosome P-site offsets (BPPS). We empirically demonstrate that our read length selection technique modestly improves sensitivity by identifying more canonical and non-canonical ORFs. Proteomics- and quantitative translation initiation sequencing-based validation verifies the high quality of all of the predictions. Experimental comparison shows that Rp-Bp results in more peptide identifications and proteomics-validated ORF predictions compared to another recent tool for translation prediction.

Macroautophagy, a key player in protein quality control, is proposed to be systematically impaired in distinct tissues and causes coordinated disruption of protein homeostasis and ageing throughout the body. Although tissue-specific changes in autophagy and ageing have been extensively explored, the mechanism underlying the inter-tissue regulation of autophagy with ageing is poorly understood. Here, we show that a secreted microRNA, mir-83/miR-29, controls the age-related decrease in macroautophagy across tissues in Caenorhabditis elegans. Upregulated in the intestine by hsf-1/HSF1 with age, mir-83 is transported across tissues potentially via extracellular vesicles and disrupts macroautophagy by suppressing CUP-5/MCOLN, a vital autophagy regulator, autonomously in the intestine as well as non-autonomously in body wall muscle. Mutating mir-83 thereby enhances macroautophagy in different tissues, promoting protein homeostasis and longevity. These findings thus identify a microRNA-based mechanism to coordinate the decreasing macroautophagy in various tissues with age., Decreased autophagy is a hallmark of ageing, but its inter-tissue regulation is poorly understood. Here, Zhou et al. identify mir-83 in C. elegans, which is transported across tissues and suppresses autophagy, contributing to age-related decline.

Adenosine-to-inosine RNA editing and pre-mRNA splicing largely occur cotranscriptionally and influence each other. Here, we use mice deficient in either one of the two editing enzymes ADAR (ADAR1) or ADARB1 (ADAR2) to determine the transcriptome-wide impact of RNA editing on splicing across different tissues. We find that ADAR has a 100× higher impact on splicing than ADARB1, although both enzymes target a similar number of substrates with a large common overlap. Consistently, differentially spliced regions frequently harbor ADAR editing sites. Moreover, catalytically dead ADAR also impacts splicing, demonstrating that RNA binding of ADAR affects splicing. In contrast, ADARB1 editing sites are found enriched 5′ of differentially spliced regions. Several of these ADARB1-mediated editing events change splice consensus sequences, therefore strongly influencing splicing of some mRNAs. A significant overlap between differentially edited and differentially spliced sites suggests evolutionary selection toward splicing being regulated by editing in a tissue-specific manner.

UBR5 is an E3 ubiquitin ligase involved in distinct processes such as transcriptional regulation and development. UBR5 is highly upregulated in embryonic stem cells (ESCs), whereas its expression decreases with differentiation, suggesting a role for UBR5 in ESC function. However, little is known about how UBR5 regulates ESC identity. Here, we define the protein interactome of UBR5 in ESCs and find interactions with distinct components of the H/ACA ribonucleoprotein complex, which is required for proper maturation of ribosomal RNA (rRNA). Notably, loss of UBR5 induces an abnormal accumulation of rRNA processing intermediates, resulting in diminished ribosomal levels. Consequently, lack of UBR5 triggers an increase in p53 levels and a concomitant decrease in cellular proliferation rates. Thus, our results indicate a link between UBR5 and rRNA maturation.

Invasive mycoses remain underdiagnosed and difficult to treat. Hospitalized individuals with compromised immunity increase in number and constitute the main risk group for severe fungal infections. Current antifungal therapy is hampered by slow and insensitive diagnostics and frequent toxic side effects of standard antifungal drugs. Identification of new antifungal compounds with high efficacy and low toxicity is therefore urgently required. We investigated the antifungal activity of tempol, a cell-permeable nitroxide. To narrow down possible mode of action we used RNA-seq technology and metabolomics to probe for pathways specifically disrupted in the human fungal pathogen Candida albicans due to tempol administration. We found genes upregulated which are involved in iron homeostasis, mitochondrial stress, steroid synthesis, and amino acid metabolism. In an ex vivo whole blood infection, tempol treatment reduced C. albicans colony forming units and at the same time increased the release of pro-inflammatory cytokines, such as interleukin 8 (IL-8, monocyte chemoattractant protein-1, and macrophage migration inhibitory factor). In a systemic mouse model, tempol was partially protective with a significant reduction of fungal burden in the kidneys of infected animals during infection onset. The results obtained propose tempol as a promising new antifungal compound and open new opportunities for the future development of novel therapies. Originally included in thesis in manuscript form

Supplemental Digital Content is available in the text., Rationale: Gene expression profiles have been mainly determined by analysis of transcript abundance. However, these analyses cannot capture posttranscriptional gene expression control at the level of translation, which is a key step in the regulation of gene expression, as evidenced by the fact that transcript levels often poorly correlate with protein levels. Furthermore, genome-wide transcript profiling of distinct cell types is challenging due to the fact that lysates from tissues always represent a mixture of cells. Objectives: This study aimed to develop a new experimental method that overcomes both limitations and to apply this method to perform a genome-wide analysis of gene expression on the translational level in response to pressure overload. Methods and Results: By combining ribosome profiling (Ribo-seq) with a ribosome-tagging approach (Ribo-tag), it was possible to determine the translated transcriptome in specific cell types from the heart. After pressure overload, we monitored the cardiac myocyte translatome by purifying tagged cardiac myocyte ribosomes from cardiac lysates and subjecting the ribosome-protected mRNA fragments to deep sequencing. We identified subsets of mRNAs that are regulated at the translational level and found that translational control determines early changes in gene expression in response to cardiac stress in cardiac myocytes. Translationally controlled transcripts are associated with specific biological processes related to translation, protein quality control, and metabolism. Mechanistically, Ribo-seq allowed for the identification of upstream open reading frames in transcripts, which we predict to be important regulators of translation. Conclusions: This method has the potential to (1) provide a new tool for studying cell-specific gene expression at the level of translation in tissues, (2) reveal new therapeutic targets to prevent cellular remodeling, and (3) trigger follow-up studies that address both, the molecular mechanisms involved in the posttranscriptional control of gene expression in cardiac cells, and the protective functions of proteins expressed in response to cellular stress.

Circular RNAs (circRNAs) are abundant and accumulate with age in neurons of diverse species. However, only few circRNAs have been functionally characterized, and their role during aging has not been addressed. Here, we use transcriptome profiling during aging and find that accumulation of circRNAs is slowed down in long-lived insulin mutant flies. Next, we characterize the in vivo function of a circRNA generated by the sulfateless gene (circSfl), which is consistently upregulated, particularly in the brain and muscle, of diverse long-lived insulin mutants. Strikingly, lifespan extension of insulin mutants is dependent on circSfl, and overexpression of circSfl alone is sufficient to extend the lifespan. Moreover, circSfl is translated into a protein that shares the N terminus and potentially some functions with the full-length Sfl protein encoded by the host gene. Our study demonstrates that insulin signaling affects global circRNA accumulation and reveals an important role of circSfl during aging in vivo.

Background RNA-binding proteins (RBPs) are fundamental regulators of cellular biology that affect all steps in the generation and processing of RNA molecules. Recent evidence suggests that regulation of RBPs that modulate both RNA stability and translation may have a profound effect on the proteome. However, regulation of RBPs in clinically relevant experimental conditions has not been studied systematically. Methods We used RNA interactome capture, a method for the global identification of RBPs to characterize the global RNA‐binding proteome (RBPome) associated with polyA-tailed RNA species in murine ciliated epithelial cells of the inner medullary collecting duct. To study regulation of RBPs in a clinically relevant condition, we analyzed hypoxia-associated changes of the RBPome. Results We identified >1000 RBPs that had been previously found using other systems. In addition, we found a number of novel RBPs not identified by previous screens using mouse or human cells, suggesting that these proteins may be specific RBPs in differentiated kidney epithelial cells. We also found quantitative differences in RBP-binding to mRNA that were associated with hypoxia versus normoxia. Conclusions These findings demonstrate the regulation of RBPs through environmental stimuli and provide insight into the biology of hypoxia-response signaling in epithelial cells in the kidney. A repository of the RBPome and proteome in kidney tubular epithelial cells, derived from our findings, is freely accessible online, and may contribute to a better understanding of the role of RNA-protein interactions in kidney tubular epithelial cells, including the response of these cells to hypoxia.

CircRNAs are a subgroup of RNAs which form a circular molecule. During the splicing process the 3’ splice donor loops back to form a covalent bond with an upstream 5’ splice acceptor instead of the downstream 5’ acceptor. The majority of circRNAs are non-coding splice isoforms of protein coding genes showing tissue and time specific expression. Because circRNAs have no 5’ cap nor 3’ poly-A tail, they degrade slower than their linear host genes. Only few studies were able to show a specific function for circRNAs.In this work we sequenced 23 samples from three tissues (brain, muscle, and gut) at three (two for gut) time points throughout the life of the naturally short-lived African turquoise killifish (Nothobranchius furzeri). We identified 1810 unique circRNAs, half of which are conserved with humans and mice. With this study we provide a comprehensive atlas to the circRNA landscape in the ageing African turquoise killifish, which serves as a novel resource to the circRNA as well as the ageing community.

RHOA is a member of the Rho family of GTPases that are involved in fundamental cellular processes including cell adhesion, migration, and proliferation. RHOA can stimulate the formation of stress fibers and focal adhesions and is a key regulator of actomyosin dynamics in various tissues. In a Genematcher-facilitated collaboration, we were able to identify four unrelated individuals with a specific phenotype characterized by hypopigmented areas of the skin, dental anomalies, body asymmetry, and limb length discrepancy due to hemihypotrophy of one half of the body, as well as brain magnetic resonance imaging (MRI) anomalies. Using whole-exome and ultra-deep amplicon sequencing and comparing genomic data of affected and unaffected areas of the skin, we discovered that all four individuals carried the identical RHOA missense variant, c.139G>A; p.Glu47Lys, in a postzygotic state. Molecular modeling and in silico analysis of the affected p.Glu47Lys residue in RHOA indicated that this exchange is predicted to specifically alter the interaction of RHOA with its downstream effectors containing a PKN-type binding domain and thereby disrupts its ability to activate signaling. Our findings indicate that the recurrent postzygotic RHOA missense variant p.Glu47Lys causes a specific mosaic disorder in humans. peerReviewed

Rationale: Endoplasmic reticulum (ER) stress dysregulates ER proteostasis, which activates the transcription factor, ATF6 (activating transcription factor 6α), an inducer of genes that enhance protein folding and restore ER proteostasis. Because of increased protein synthesis, it is possible that protein folding and ER proteostasis are challenged during cardiac myocyte growth. However, it is not known whether ATF6 is activated, and if so, what its function is during hypertrophic growth of cardiac myocytes. Objective: To examine the activity and function of ATF6 during cardiac hypertrophy. Methods and Results: We found that ER stress and ATF6 were activated and ATF6 target genes were induced in mice subjected to an acute model of transverse aortic constriction, or to free-wheel exercise, both of which promote adaptive cardiac myocyte hypertrophy with preserved cardiac function. Cardiac myocyte-specific deletion of Atf6 (ATF6 cKO [conditional knockout]) blunted transverse aortic constriction and exercise-induced cardiac myocyte hypertrophy and impaired cardiac function, demonstrating a role for ATF6 in compensatory myocyte growth. Transcript profiling and chromatin immunoprecipitation identified RHEB (Ras homologue enriched in brain) as an ATF6 target gene in the heart. RHEB is an activator of mTORC1 (mammalian/mechanistic target of rapamycin complex 1), a major inducer of protein synthesis and subsequent cell growth. Both transverse aortic constriction and exercise upregulated RHEB , activated mTORC1, and induced cardiac hypertrophy in wild type mouse hearts but not in ATF6 cKO hearts. Mechanistically, knockdown of ATF6 in neonatal rat ventricular myocytes blocked phenylephrine- and IGF1 (insulin-like growth factor 1)-mediated RHEB induction, mTORC1 activation, and myocyte growth, all of which were restored by ectopic RHEB expression. Moreover, adeno-associated virus 9- RHEB restored cardiac growth to ATF6 cKO mice subjected to transverse aortic constriction. Finally, ATF6 induced RHEB in response to growth factors, but not in response to other activators of ATF6 that do not induce growth, indicating that ATF6 target gene induction is stress specific. Conclusions: Compensatory cardiac hypertrophy activates ER stress and ATF6, which induces RHEB and activates mTORC1. Thus, ATF6 is a previously unrecognized link between growth stimuli and mTORC1-mediated cardiac growth.

m 6 A is the most abundant internal modification on mRNA. Recent improvements of high-throughput sequencing techniques enables its detection at the transcriptome level, even at the nucleotide resolution. However most current techniques require large amounts of starting material to detect the modification. Here, we describe a complementary technique of standard meRIP-seq/miCLIP-seq approaches to identify methylated RNA using a low amount of material. We believe this approach can be applied in vivo to identify methylated targets in specific tissues or subpopulations of cells.

Circular RNAs (circRNAs) are a recent addition to the expanding universe of RNA species and originate through back-splicing events from linear primary transcripts. CircRNAs show specific expression profiles with regards to cell type and developmental stage. Importantly, only few circRNAs have been functionally characterized to date. The detection of circRNAs from RNA sequencing data is a complex computational workflow that, depending on tissue and condition typically yields candidate sets of hundreds or thousands of circRNA candidates. Here, we provide an overview on different computational analysis tools and pipelines that became available throughout the last years. We outline technical and experimental requirements that are common to all approaches and point out potential pitfalls during the computational analysis. Although computational prediction of circRNAs has become quite mature in recent years, we provide a set of valuable validation strategies, in silico as well as in vitro-based approaches. In addition to circRNA detection via back-splicing junction, we present available analysis pipelines for delineating the primary sequence and for predicting possible functions of circRNAs. Finally, we outline the most important web resources for circRNA research. This article is categorized under: RNA Methods > RNA Analyses in vitro and In Silico RNA Evolution and Genomics > Computational Analyses of RNA.

Background Huntington’s Disease (HD), is a hereditary, fatal neurodegenerative disorder due to CAG trinucleotide expansion in exon 1 of the HD gene (HTT). Recent findings revealed that the process of alternative splicing (AS) might be compromised in HD. AS regulation is crucial not only to the establishment of a repertoire of protein coding isoforms extremely relevant for the proper physiological characteristics of the nervous system, but also to the biogenesis of circular RNAs (circRNAs), unusually stable, highly expressed non-coding RNAs produced by the circularization of exons which, seems to have important roles during neuronal development and functioning. Aims Therefore, the overarching goal of our research was to take advantage of genetically engineered Htt knock-in (KI) mouse cellular and in vivo model systems bearing normal or pathological Htt CAG repeat lengths to discover alterations in linear and back-splicing events that might reflect on the levels of expression as well as the composition of a repertoire of proteins, thus contributing to HD striatal vulnerability and pathogenesis. Methods/techniques To discover alterations in linear splicing, we resourced to a publicly available dataset (Langfelder P. et al., 2016, Nature Neuroscience) analyzing full transcriptomic profiling dataset for striatum, cortex and liver at different developmental time points, identifying and quantifying the total numbers of differential AS events for each genotype/sample type and time point. Results/outcome Our analysis demonstrated that specifically for the striatum – most vulnerable to HD degeneration – the total number of detected, differential AS events increased significantly at 6 and 10 month with increasing CAG expansion. Interestingly, not all the AS events increased with CAG-expansion and age, but only a specific subtype of AS event – EXON SKIPPING – revealed to be strongly and specifically affected. On the other hand, we identified and characterized the first circular RNA originating locally from the HTT locus, expressed in whole body, but predominantly in the brain and spinal cord and presenting augmented expression level with increasing HTT CAG size. Moreover, at genome-wide level, data analysis of circRNA-seq showed a decreased circRNAs production when mutant huntingtin is expressed. Specifically, 12 circRNAs, identified by stringent cut-off criteria, showed continuous decreased expression following CAG expansion in neuronal progenitor cells. Since, most of the CAG-sensitive circRNAs have annotated human orthologues, their expression can be further characterized and functionally studied in human cell lines. Conclusions In conclusion, our results support the idea that AS machinery is responding to HD mutational process altering both linear and back-splicing events locally at the HTT locus, but also at genome-wide level. This knowledge will pave the way to new trials of therapeutic intervention aimed to possibly target spliceosomal-circRNAs alterations through specific drugs or genetic manipulation.

Down-regulation of cardiac inotropic EF-hand calcium sensor protein S100A1 expression is a hallmark of human heart failure and drives progression and mortality. An AAV-S100A1 gene replacement therapy, which has been shown in preclinical models to rescue chronic heart failure, is in preparation for first clinical trials. However, the nuclear mechanisms that convey S100A1 gene locus silencing in the failing heart have not been determined yet. Given the steady increase of S100A1 expression in the developing heart until young adulthood and its decrease in the aged and failing heart, we developed a novel approach to decipher the mechanisms of S100a1 transcriptional regulation. We performed serial RNAseq analyses in an in vitro cardiac differentiation model of H9c2 cells that displays robust concurrent S100A1 and e.g. sarcoplasmic reticulum calcium ATPase 2 (SERCA2a) mRNA and protein during differentiation. Potential transcription factors (TF) were derived by computational motif analysis from S100A1’s gene locus regulatory elements that delivered more than 200 candidates. Subsequent computational TF footprint binding motive analysis from all differentially expressed transcripts of the serial RNAseq data provided a smaller group of less than 30 putatively TF families in the cH9C2 model that is actively during the in vitro biological process. Additional comparative TF expression analyses with developing rat hearts then enabled subtraction of a final TF group of 10 members for RNA interference screens in the cH9C2 model using S100A1 RNA expression as a read-out. Using a siRNA-mediated TF knock-down screen, we identified a small group of TFs such as zink finger and coiled-coil TFs (e.g. KLFs) and signal transducer and activator of transcription (STATs) TFs with an hitherto unknown ability to independently regulate S100A1 expression. Of note, other inotropic factors such as SERCA2a were found to be positively co-regulated by some of the TFs, such as KLFs, that enhance S100A1 expression. Further analyses were performed to confirm these results. In summary, we present a feasible computational filter and experimental strategy to identify relevant TF networks which we have successfully applied to identify cardiac regulatory networks.

Aberrant synaptic function is thought to underlie social deficits in neurodevelopmental disorders such as autism and schizophrenia. Although microRNAs have been shown to regulate synapse development and plasticity, their potential involvement in the control of social behaviour in mammals remains unexplored. Here, we show that deletion of the placental mammal‐specific miR379‐410 cluster in mice leads to hypersocial behaviour, which is accompanied by increased excitatory synaptic transmission, and exaggerated expression of ionotropic glutamate receptor complexes in the hippocampus. Bioinformatic analyses further allowed us to identify five “hub” microRNAs whose deletion accounts largely for the upregulation of excitatory synaptic genes observed, including Cnih2, Dlgap3, Prr7 and Src. Thus, the miR379‐410 cluster acts a natural brake for sociability, and interfering with specific members of this cluster could represent a therapeutic strategy for the treatment of social deficits in neurodevelopmental disorders.

Circular RNAs (circRNAs) have been first described as "scrambled exons" in the 1990s. CircRNAs originate from back splicing or exon skipping of linear RNA templates and have continuously gained attention in recent years due to the availability of high-throughput whole-transcriptome sequencing methods. Numerous manuscripts describe thousands of circRNAs throughout uni- and multicellular eukaryote species and demonstrated that they are conserved, stable, and abundant in specific tissues or conditions. This manuscript provides a walk-through of our bioinformatics toolbox, which covers all aspects of in silico circRNA analysis, starting from raw sequencing data and back-splicing junction discovery to circRNA quantitation and reconstruction of internal the circRNA structure.

Gene expression profiling provides a tool to analyze the internal states of cells or organisms, and their responses to perturbations. While global measurements of mRNA levels have thus been widely used for many years, it is only through the recent development of the ribosome profiling technique that an analogous examination of global mRNA translation programs has become possible. Ribosome profiling reveals which RNAs are being translated to what extent and where the translated open reading frames are located. In addition, different modes of translation regulation can be distinguished and characterized. Here, we present an optimized, step-by-step protocol for ribosome profiling. Although established in Caenorhabditis elegans, our protocol and optimization approaches should be equally usable for other model organisms or cell culture with little adaptation. Next to providing a protocol, we compare two different methods for isolation of single ribosomes and two different library preparations, and describe strategies to optimize the RNase digest and to reduce ribosomal RNA contamination in the libraries. Moreover, we discuss bioinformatic strategies to evaluate the quality of the data and explain how the data can be analyzed for different applications. In sum, this article seeks to facilitate the understanding, execution, and optimization of ribosome profiling experiments.

In metazoan cells, spliced mRNAs are marked by the exon junction complex (EJC), a multi-protein complex that serves as a key regulator of post-transcriptional mRNA metabolism. Deposition of EJCs on mRNA is intimately linked to the splicing process. The spliceosomal protein CWC22 directly binds the core EJC-protein eIF4A3, guides it to the spliceosome and initiates EJC assembly. In addition, CWC22 is involved in the splicing process itself, but the molecular details of its dual function remain elusive. Here we analyze the mechanisms, by which CWC22 co-regulates pre-mRNA splicing and EJC assembly. We show that the core of CWC22 is sufficient to mediate both pre-mRNA splicing and EJC assembly. Nonetheless, both processes can be functionally uncoupled with an eIF4A3-binding deficient mutant of CWC22, which impedes EJC assembly. A C-terminal domain of CWC22 strongly enhances its spliceosomal interaction and likely regulates its function. High-throughput RNA-sequencing identifies global defects of pre-mRNA splicing and downregulation of diverse gene expression pathways in CWC22-depleted cells. We propose a model, in which CWC22 represents an integral component of the spliceosome and orchestrates pre-mRNA splicing and eIF4A3 binding to achieve global assembly of exon junction complexes.

Mitofusin 2 plays an unexpected role in maintaining the terpenoid biosynthesis pathway and is necessary for mitochondrial coenzyme Q biosynthesis., Mitochondria form a dynamic network within the cell as a result of balanced fusion and fission. Despite the established role of mitofusins (MFN1 and MFN2) in mitochondrial fusion, only MFN2 has been associated with metabolic and neurodegenerative diseases, which suggests that MFN2 is needed to maintain mitochondrial energy metabolism. The molecular basis for the mitochondrial dysfunction encountered in the absence of MFN2 is not understood. Here we show that loss of MFN2 leads to impaired mitochondrial respiration and reduced ATP production, and that this defective oxidative phosphorylation process unexpectedly originates from a depletion of the mitochondrial coenzyme Q pool. Our study unravels an unexpected and novel role for MFN2 in maintenance of the terpenoid biosynthesis pathway, which is necessary for mitochondrial coenzyme Q biosynthesis. The reduced respiratory chain function in cells lacking MFN2 can be partially rescued by coenzyme Q10 supplementation, which suggests a possible therapeutic strategy for patients with diseases caused by mutations in the Mfn2 gene.