Search

Background: ADAR1 (adenosine deaminase acting on RNA-1)-mediated adenosine to inosine (A-to-I) RNA editing plays an essential role for distinguishing endogenous from exogenous RNAs, preventing autoinflammatory ADAR1 also regulates cellular processes by recoding specific mRNAs, thereby altering protein functions, but may also act in an editing-independent manner. The specific role of ADAR1 in cardiomyocytes and its mode of action in the heart is not fully understood. To determine the role of ADAR1 in the heart, we used different mutant mouse strains, which allows to distinguish immunogenic, editing-dependent, and editing-independent functions of ADAR1. Methods: Different Adar1 -mutant mouse strains were employed for gene deletion or specific inactivation of ADAR1 enzymatic activity in cardiomyocytes, either alone or in combination with Ifih1 (interferon induced with helicase C domain 1) or Irf7 (interferon regulatory factor 7) gene inactivation. Mutant mice were investigated by immunofluorescence, Western blot, RNAseq, proteomics, and functional MRI analysis. Results: Inactivation of Adar1 in cardiomyocytes resulted in late-onset autoinflammatory myocarditis progressing into dilated cardiomyopathy and heart failure at 6 months of age. Adar1 depletion activated interferon signaling genes but not NFκB (nuclear factor kappa B) signaling or apoptosis and reduced cardiac hypertrophy during pressure overload via induction of Irf7 . Additional inactivation of the cytosolic RNA sensor MDA5 (melanoma differentiation-associated gene 5; encoded by the Ifih1 gene) in Adar1 mutant mice prevented activation of interferon signaling gene and delayed heart failure but did not prevent lethality after 8.5 months. In contrast, compound mutants only expressing catalytically inactive ADAR1 in an Ifih1 -mutant background were completely normal. Inactivation of Irf7 attenuated the phenotype of Adar1 -deficient cardiomyocytes to a similar extent as Ifih1 depletion, identifying IRF7 as the main mediator of autoinflammatory responses caused by the absence of ADAR1 in cardiomyocytes. Conclusions: Enzymatically active ADAR1 prevents IRF7-mediated autoinflammatory reactions in the heart triggered by endogenous nonedited RNAs. In addition to RNA editing, ADAR1 also serves editing-independent roles in the heart required for long-term cardiac function and survival.

Circular RNAs are endogenous long-lived and abundant non-coding species. Despite their prevalence, only few circRNAs have been dissected mechanistically to date. Here, we catalogued nascent RNA-enriched circRNAs from primary human cells and functionally assign a role to circRAB3IP in sustaining cellular homeostasis. We combined “omics” and functional experiments to show how circRAB3IP depletion deregulates hundreds of genes, suppresses cell cycle progression, and induces senescence-associated gene expression changes. Conversely, excess circRAB3IP delivered to endothelial cells via extracellular vesicles suffices for accelerating their division. We attribute these effects to the interplay between circRAB3IP and the general splicing factor SF3B1, which impacts transcript usage of cell cycle-related genes. Together, our findings link the maintenance of cell homeostasis to the tightly regulated titers of a single circRNA.

Background: The human heart has limited capacity to generate new cardiomyocytes and this capacity declines with age. Because loss of cardiomyocytes may contribute to heart failure, it is crucial to explore stimuli of endogenous cardiac regeneration to favorably shift the balance between loss of cardiomyocytes and the birth of new cardiomyocytes in the aged heart. We have previously shown that cardiomyogenesis can be activated by exercise in the young adult mouse heart. Whether exercise also induces cardiomyogenesis in aged hearts, however, is still unknown. Here, we aim to investigate the effect of exercise on the generation of new cardiomyocytes in the aged heart. Methods: Aged (20-month-old) mice were subjected to an 8-week voluntary running protocol, and age-matched sedentary animals served as controls. Cardiomyogenesis in aged hearts was assessed on the basis of 15 N-thymidine incorporation and multi-isotope imaging mass spectrometry. We analyzed 1793 cardiomyocytes from 5 aged sedentary mice and compared these with 2002 cardiomyocytes from 5 aged exercised mice, followed by advanced histology and imaging to account for ploidy and nucleation status of the cell. RNA sequencing and subsequent bioinformatic analyses were performed to investigate transcriptional changes induced by exercise specifically in aged hearts in comparison with young hearts. Results: Cardiomyogenesis was observed at a significantly higher frequency in exercised compared with sedentary aged hearts on the basis of the detection of mononucleated/diploid 15 N-thymidine–labeled cardiomyocytes. No mononucleated/diploid 15 N-thymidine–labeled cardiomyocyte was detected in sedentary aged mice. The annual rate of mononucleated/diploid 15 N-thymidine–labeled cardiomyocytes in aged exercised mice was 2.3% per year. This compares with our previously reported annual rate of 7.5% in young exercised mice and 1.63% in young sedentary mice. Transcriptional profiling of young and aged exercised murine hearts and their sedentary controls revealed that exercise induces pathways related to circadian rhythm, irrespective of age. One known oscillating transcript, however, that was exclusively upregulated in aged exercised hearts, was isoform 1.4 of regulator of calcineurin, whose regulation and functional role were explored further. Conclusions: Our data demonstrate that voluntary running in part restores cardiomyogenesis in aged mice and suggest that pathways associated with circadian rhythm may play a role in physiologically stimulated cardiomyogenesis.

Tumor purity of matched tumor/cell line pairs: Tumor purity of GBM bulk tumor samples and matched gliomasphere cell lines was determined either using ESTIMATE for deconvolution of gene expression profiles (A) or pTERT MAF (B, C). Connected dots indicate matched pairs.

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.

Primers and annealing temperatures used in quantitative allele-specific real-time PCR.

Telomerase promoter genotype and MGMT status in matched glioblastoma tumor vs. cell line pairs.

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 detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed by computational detection tools. During the last decade, a plethora of such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools were used and detected over 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were empirically validated using three orthogonal methods. Generally, tool-specific precision values are high and similar (median of 98.8%, 96.3%, and 95.5% for qPCR, RNase R, and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant tool differentiators. Furthermore, we demonstrate the complementarity of tools through the increase in detection sensitivity by considering the union of highly-precise tool combinations while keeping the number of false discoveries low. Finally, based on the benchmarking results, recommendations are put forward for circRNA detection and validation.

Background and Hypothesis: Expression of the cardiomyocyte (CM) protein S100A1, which improves contractile performance of the heart, sharply increases during postnatal myocardial maturation but declines rapidly in failing hearts. We therefore hypothesized that CMs are wired with transcriptional factors (TFs) that positively and negatively regulate S100A1’s gene locus activity. Understanding these reciprocal circuits may be relevant for advanced therapeutic modulation of S100A1’s abundance in diseased hearts. Methods and Results: H9C2 rat cardiomyoblasts, an animal-free in vitro tool, displayed a strong concordant rise in S100A1 mRNA and protein levels (8.1+/-1.1 vs. cont.; n=9, p2-fold enrichment over control were selected delivering EGR1 and SP1 as top hits. Subsequent siRNA-mediated silencing of both TFs yielded an EGR1 and SP1 knock-down dose-dependent inhibition (80%; p Conclusion: Our study identified the known cardioprotective EGR1 and maladaptive SP1 as novel positive and negative TF regulators of CM S100A1 expression and targets for therapeutic S100A1 gene locus modulation i.e., by SP1 TFBS gene-editing.

Post-transcriptional gene regulation is accomplished by the interplay of the transcriptome with RNA-binding proteins, which occurs in a dynamic manner in response to altered cellular conditions. Recording the combined occupancy of all proteins binding to the transcriptome offers the opportunity to interrogate if a particular treatment leads to any interaction changes, pointing to sites in RNA that undergo post-transcriptional regulation. Here, we establish a method to monitor protein occupancy in a transcriptome-wide fashion by RNA sequencing. To this end, peptide-enhanced pull-down for RNA sequencing (or PEPseq) uses metabolic RNA labelling with 4-thiouridine (4SU) for light-induced protein-RNA crosslinking, and N-hydroxysuccinimide (NHS) chemistry to isolate protein-crosslinked RNA fragments across all long RNA biotypes. We use PEPseq to investigate changes in protein occupancy during the onset of arsenite-induced translational stress in human cells and reveal evidence for ribosome stalling and depletion from stress granules for a distinct set of mRNAs, many coding for ribosomal proteins. We use quantitative proteomics to demonstrate that translation of these mRNAs remains repressed during the initial hours of recovery after arsenite stress. Thus, we present PEPseq as a discovery platform for the unbiased investigation of post-transcriptional regulation.

We introduce Single-cell Nanopore Spatial Transcriptomics (scNaST), a software suite to facilitate the analysis of spatial gene expression from second- and third-generation sequencing, allowing to generate a full-length near-single-cell transcriptional landscape of the tissue microenvironment. Taking advantage of the Visium Spatial platform, we adapted a strategy recently developed to assign barcodes to long-read single-cell sequencing data for spatial capture technology. Here, we demonstrate our workflow using four short axis sections of the mouse heart following myocardial infarction. We constructed a de novo transcriptome using long-read data, and successfully assigned 19,794 transcript isoforms in total, including clinically-relevant, but yet uncharacterized modes of transcription, such as intron retention or antisense overlapping transcription. We showed a higher transcriptome complexity in the healthy regions, and identified intron retention as a mode of transcription associated with the infarct area. Our data revealed a clear regional isoform switching among differentially used transcripts for genes involved in cardiac muscle contraction and tissue morphogenesis. Molecular signatures involved in cardiac remodeling integrated with morphological context may support the development of new therapeutics towards the treatment of heart failure and the reduction of cardiac complications.

Coronavirus disease 2019 (COVID-19) has a highly variable disease severity. Possible associations between peripheral blood signatures and disease severity have been investigated since the emergence of the pandemic. Although there appear to be several signatures identified based on exploratory analyses of single-cell omics data, there are no state-of-the-art validated models to predict COVID-19 severity from comprehensive transcriptome profiling of Peripheral Blood Mononuclear Cells (PBMCs) across multiple sites. In this manuscript, we present a computational workflow based on a Multilayer perceptron (MLP) network that predicts disease severity from PBMCs single-cell RNA-seq data of COVID-19 patients. The study includes patient cohorts from different sites: the University Hospital in Bonn, the Charité in Berlin, the Stanford University COVID-19 Biobanking studies, and three Korean medical centers; all accompanied by their severity status in terms of mechanical ventilation necessity. Training and model validation are performed on randomly selected samples solely from the Berlin and Bonn cohorts, while testing is performed on completely unseen samples from the Stanford and Korean datasets. Our proposed predictive model shows a high area under the receiver operating characteristic (AUROC) curve on the testing datasets (1 (CI:1-1) Korea, 0.86 (CI:0.81-0.9) Stanford), proving our model's robustness. Moreover, we identified a number of features that contributed strongly to the prediction, some of which have been reported previously to have a strong relation with the severity state. In summary, we could show that the expression of 15 genes and the cell proportion profile of 29 PBMC cell types alone are sufficient to distinguish between severe and mild COVID-19 disease states. Our model is publically available at https://github.com/dieterich-lab/ImmunOMICS.

The catalytically active component of ribosomes, rRNA, is long studied and heavily modified. However, little is known about functional and pathological consequences of changes in human rRNA modification status. Direct RNA-sequencing on the Nanopore platform enables the direct assessment of rRNA modifications. We established a targeted Nanopore direct rRNA sequencing approach and applied it to CRISPR-Cas9 engineered HCT116 cells, lacking specific enzymatic activities required to establish defined rRNA base modifications. We analyzed these sequencing data along with wild type samples and in vitro transcribed reference sequences to specifically detect changes in modification status. We show for the first time that direct RNA-sequencing is feasible on smaller, i.e. Flongle, flow cells. Our targeted approach reduces RNA input requirements, making it accessible to the analysis of limited samples such as patient derived material. The analysis of rRNA modifications during cardiomyocyte differentiation of human induced pluripotent stem cells, and of heart biopsies from cardiomyopathy patients revealed altered modifications of specific sites, among them pseudouridines, 2’-O-methylation of riboses and an acetylation of cytidine that was validated by orthogonal methods and the use of a genetic model. Targeted direct rRNA-sequencing analysis with JACUSA2 opens up the possibility to analyze dynamic changes in rRNA modifications in a wide range of biological and clinical samples.

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.

ONT long-read sequencing provides real-time monitoring and controlling of individual nanopores. Adaptive sampling enriches or depletes specific sequences in Nanopore DNA sequencing, but was not applicable to direct sequencing of RNA so far.Here, we identify essential parameter settings for direct RNA sequencing (DRS). We demonstrate the superior performance of depletion over enrichment and show that adaptive sampling efficiently depletes specific transcripts in transcriptome-wide sequencing applications. Specifically, we applied our adaptive sampling approach to polyA+ RNA samples from human cardiomyocytes and mouse whole heart tissue. Herein, we show more than 2.5-fold depletion of highly abundant mitochondrial-encoded transcripts that in normal sequencing account for up to 40% of sequenced bases in heart tissue samples.

A whole series of high-throughput antibody-free methods for RNA modification detection from sequencing data emerged lately. We present JACUSA2 as a versatile software solution and comprehensive analysis framework for RNA modification detection assays that are based on either the Illumina or Nanopore platform. Importantly, JACUSA2 can integrate information from multiple experiments (e.g. replicates and different conditions) and different library types (e.g. first- or secondstrand libraries). We demonstrate its utility by example, showing three analysis workflows for m6A detection on published data sets: 1) MazF m6a-sensitive RNA digestion (FTO+ vs FTO-), 2) DART-seq (YTHwt vs YTHmut) and 3) Nanopore profiling (METTL3 +/+ vs -/-). All assays have been conducted in HEK293 cells and complement one another.

Aims: Takotsubo cardiomyopathy (TTC) presents an acute heart failure (AHF) syndrome that mimics the symptoms of acute myocardial infarction and is often preceded by an episode of emotional and/or physical stress. There is no specific treatment available, and evidence is based mostly on retrospective data. Here, we show that a high-dose of epinephrine (EPI) recapitulates numerous features of human TTC in mice, enabling genetic KO studies. By gene set enrichment network analysis and studies involving myocardial KO of the endogenous calcineurin inhibitor calcium-calmodulin-dependent kinase II (CaMKII), we identify cyclosporine A (CSA) as a novel potential therapy for TTC.Methods and results: high-dose EPI caused reversible AHF with significant ST-segment changes and high-sensitive Troponin T (hs-TnT) elevation. Male mice displayed markedly elevated hs-TnT and mortality when compared to females, mimicking the clinical syndrome. Left ventricular (LV) catecholamines revealed a distinct pattern of blunted norepinephrine as opposed to elevated EPI in male vs. female mice. Gene set enrichment analysis of LV RNA-sequencing after EPI revealed calcineurin-dependent pro-inflammatory gene networks to be significantly enriched, particularly in male hearts. Cardiac-specific deletion of CaMKII caused exacerbated AHF and myocardial damage, that was attenuated by CSA. In human peripheral blood mononuclear cells from TTC patients, expression of the calcineurin reporter regulator of calcineurin 1 isoform 4 (rcan1-4) was significantly increased compared to age- and gender matched controls, indicating human relevance of this pathway.Conclusion: Combining TTC disease modeling and mouse genetics identified calcineurin inhibition as a novel specific therapeutic concept to treat human TTC.

The paralogous human proteins UPF3A and UPF3B are involved in recognizing mRNAs targeted by nonsense-mediated mRNA decay (NMD). UPF3B has been demonstrated to support NMD, presumably by bridging an exon junction complex (EJC) to the NMD factor UPF2. The role of UPF3A has been described either as a weak NMD activator or an NMD inhibitor. Here, we present a comprehensive functional analysis of UPF3A and UPF3B in human cells using combinatory experimental approaches. Overexpression or knockout of UPF3A as well as knockout of UPF3B did not substantially change global NMD activity. In contrast, the co-depletion of UPF3A and UPF3B resulted in a marked NMD inhibition and a transcriptome-wide upregulation of NMD substrates, demonstrating a functional redundancy between both NMD factors. In rescue experiments, UPF2 or EJC binding-deficient UPF3B largely retained NMD activity. However, combinations of different mutants, including deletion of the middle domain, showed additive or synergistic effects and therefore failed to maintain NMD. Collectively, UPF3A and UPF3B emerge as fault-tolerant, functionally redundant NMD activators in human cells.

Internal modifications of mRNA have emerged as widespread and versatile regulatory mechanism to control gene expression at the post-transcriptional level. Current insights rely on the ability to make a modified nucleoside amenable to sequencing. Most of the modifications are methylations involving the co-factor S-adenosyl-L-methionine (SAM), however, simultaneous detection of different methylation sites in the same sample has remained elusive. We present metabolic labeling with propargyl-selenohomocysteine (PSH) in combination with click chemistry to detect N6- methyladenosine (m6A) and 5-methylcytidine (m5C) sites in mRNA with single nucleotide precision in the same sequencing run (MePMe-seq). Intracellular formation of the corresponding SAM analogue leads to detectable levels of N6-propargyl-A (prop6A) and 5-propargyl-C (prop5C). MePMe-seq overcomes the problems of antibodies for enrichment and sequence-motifs for evaluation, limiting previous methodologies. The joint evaluation of m6A and m5C sites opens the door to study their interconnectivity and improve our understanding of mechanisms and functions of the RNA methylome.

Ageing is a complex process with common and distinct features across tissues. Unveiling the underlying processes driving ageing in individual tissues is indispensable to decipher the mechanisms of organismal longevity. Caenorhabditis elegans is a well-established model organism that has spearheaded ageing research with the discovery of numerous genetic pathways controlling its lifespan. However, it remains challenging to dissect the ageing of worm tissues due to the limited description of tissue pathology and access to tissue-specific molecular changes during ageing. In this study, we isolated cells from five major tissues in young and old worms and profiled the age-induced transcriptomic changes within these tissues. We observed a striking diversity of ageing across tissues and identified different sets of longevity regulators therein. In addition, we found novel tissue-specific factors, including irx-1 and myrf-2, which control the integrity of the intestinal barrier and sarcomere structure during ageing respectively. This study demonstrates the complexity of ageing across worm tissues and highlights the power of tissue-specific transcriptomic profiling during ageing, which can serve as a resource to the field.

RNA modifications exist in all kingdom of life. Several different types of base or ribose modifications are now summarized under the term the ”epitranscriptome”. With the advent of high-throughput sequencing technologies, much progress has been made in understanding RNA modification biology and how these modifications can influence many aspects of RNA life. The most widespread internal modification on mRNA is m6A, which has been implicated in physiological processes as well as disease pathogenesis. Here, we provide a workflow for the mapping of m6A sites using Nanopore direct RNA sequencing data. Our strategy employs pairwise comparison of base calling error profiles with JACUSA2. We outline a general strategy for RNA modification detection on mRNA and describe two specific use cases on m6A detection in detail. Use case 1: a sample of interest with modifications (e.g. ”wild type” sample) is compared to a sample lacking a specific modification type (e.g. ”knock out” sample, here METTL3 -KO) or Use case 2: a sample of interest with modifications is compared to a sample lacking all modifications (e.g. in vitro transcribed cDNA). We provide a detailed protocol on experimental and computational aspects. Extensive online material provides a snakemake pipeline to identify m6A positions in mRNA and to validate the results against a miCLIP-derived m6A reference set. The general strategy is flexible and can be easily adapted by users in different application scenarios.