Correcting 4sU induced quantification bias in nucleotide conversion RNA-seq data.

Nucleoside analogues like 4-thiouridine (4sU) are used to metabolically label newly synthesized RNA. Chemical conversion of 4sU before sequencing induces T-to-C mismatches in reads sequenced from labelled RNA, allowing to obtain total and labelled RNA expression profiles from a single sequencing library. Cytotoxicity due to extended periods of labelling or high 4sU concentrations has been described, but the effects of extensive 4sU labelling on expression estimates from nucleotide conversion RNA-seq have not been studied. Here, we performed nucleotide conversion RNA-seq with escalating doses of 4sU with short-term labelling (1h) and over a progressive time course (up to 2h) in different cell lines. With high concentrations or at later time points, expression estimates were biased in an RNA half-life dependent manner. We show that bias arose by a combination of reduced mappability of reads carrying multiple conversions, and a global, unspecific underrepresentation of labelled RNA emerging during library preparation and potentially global reduction of RNA synthesis. We developed a computational tool to rescue unmappable reads, which performed favourably compared to previous read mappers, and a statistical method, which could fully remove remaining bias. All methods developed here are freely available as part of our GRAND-SLAM pipeline and grandR package.

No evidence for epitranscriptomic m5C modification of SARS-CoV-2, HIV and MLV viral RNA.

The addition of chemical groups to cellular RNA to modulate RNA fate and/or function is summarized under the term epitranscriptomic modification. More than 170 different modifications have been identified on cellular RNA, such as tRNA, rRNA and, to a lesser extent, on other RNA types. Recently, epitranscriptomic modification of viral RNA has received considerable attention as a possible additional mechanism regulating virus infection and replication. N6-methyladenosine (m6A) and C5-methylcytosine (m5C) have been most broadly studied in different RNA viruses. Various studies, however, reported varying results with regard to number and extent of the modification. Here we investigated the m5C methylome of SARS-CoV-2, and we reexamined reported m5C sites in HIV and MLV. Using a rigorous bisulfite-sequencing protocol and stringent data analysis, we found no evidence for the presence of m5C in these viruses. The data emphasize the necessity for optimizing experimental conditions and bioinformatic data analysis.

6-Thioguanosine Monophosphate Prodrugs Display Enhanced Performance against Thiopurine-Resistant Leukemia and Breast Cancer Cells.

Thiopurines are in widespread clinical use for the treatment of immunological disorders and certain cancers. However, treatment failure due to resistance or adverse drug reactions are common, asking for new therapeutic strategies. We investigated the potential of 6-thioguanosine monophosphate (6sGMP) prodrugs to overcome resistance to 6-thioguanine. We successfully developed synthetic routes toward diverse 6sGMP prodrugs, tested their proliferation inhibitory potential in different cell lines, and examined their mode of action. Our results show that 4-acetyloxybenzyl- and cycloSaligenyl-derivatized 6sGMP prodrugs are effective antiproliferative compounds in cells that are resistant to thiopurines. We find that resistance is related to the expression of salvage pathway enzyme HGPRT. Using TUC-seq DUAL, we demonstrate the intracellular conversion of 6sGMP prodrugs into bioactive 6sGTPs. Thus, our study offers a promising strategy for thiopurine therapy by using 6sGMP prodrugs, and it suggests TUC-seq DUAL as a simple and fast method to measure the success of thiopurine therapy.

Regulation of sensory perception and motor abilities by brain-specific action of chromatin remodeling factor CHD1.

The ATP-dependent chromatin remodeling factor CHD1 (chromodomain-helicase-DNA binding protein 1) is involved in both the de novo assembly and the remodeling of chromatin. Recently, we discovered a crucial role of CHD1 in the incorporation of the histone variant H3.3 in the fly brain illustrated by widespread transcriptional upregulation and shortened lifespan in Chd1-mutant animals. Because many genes linked to sensory perception were dysregulated in Chd1-mutant heads, we studied the role of CHD1 in these processes. Here we show that Chd1-mutant flies have severe defects in their response behavior to olfactory and gustatory but not visual stimuli. Further analyses suggested that poor performance in gustatory response assays was caused by reduced motivation for foraging and feeding rather than defects in taste perception. Moreover, we show that shortened lifespan of Chd1-mutant flies is accompanied by indications of premature functional aging as suggested by defects in negative geotaxis and exploratory walking assays. The latter phenotype was rescued by neuronal re-expression of Chd1, while the olfactory defects were not. Interestingly, we found evidence for indirect regulation of the non-neuronal expression of odorant binding proteins (Obp) by neuronal expression of Chd1. Together, these results emphasize the crucial role of CHD1 activity controlling diverse neuronal processes thereby affecting healthy lifespan.

Towards a comprehensive understanding of RNA deamination: synthesis and properties of xanthosine-modified RNA.

Nucleobase deamination, such as A-to-I editing, represents an important posttranscriptional modification of RNA. When deamination affects guanosines, a xanthosine (X) containing RNA is generated. However, the biological significance and chemical consequences on RNA are poorly understood. We present a comprehensive study on the preparation and biophysical properties of X-modified RNA. Thermodynamic analyses revealed that base pairing strength is reduced to a level similar to that observed for a G•U replacement. Applying NMR spectroscopy and X-ray crystallography, we demonstrate that X can form distinct wobble geometries with uridine depending on the sequence context. In contrast, X pairing with cytidine occurs either through wobble geometry involving protonated C or in Watson-Crick-like arrangement. This indicates that the different pairing modes are of comparable stability separated by low energetic barriers for switching. Furthermore, we demonstrate that the flexible pairing properties directly affect the recognition of X-modified RNA by reverse transcription enzymes. Primer extension assays and PCR-based sequencing analysis reveal that X is preferentially read as G or A and that the ratio depends on the type of reverse transcriptase. Taken together, our results elucidate important properties of X-modified RNA paving the way for future studies on its biological significance.

Synthesis of 4-thiouridines with prodrug functionalization for RNA metabolic labeling.

Metabolic labeling has emerged as a powerful tool to endow RNA with reactive handles allowing for subsequent chemical derivatization and processing. Recently, thiolated nucleosides, such as 4-thiouridine (4sU), have attracted great interest in metabolic labeling-based RNA sequencing approaches (TUC-seq, SLAM-seq, TimeLapse-seq) to study cellular RNA expression and decay dynamics. For these and other applications (e.g. PAR-CLIP), thus far only the naked nucleoside 4sU has been applied. Here we examined the concept of derivatizing 4sU into a 5'-monophosphate prodrug that would allow for cell permeation and potentially improve labeling efficiency by bypassing the rate-limiting first step of 5' phosphorylation of the nucleoside into the ultimately bioactive 4sU triphosphate (4sUTP). To this end, we developed robust synthetic routes towards diverse 4sU monophosphate prodrugs. Using metabolic labeling assays, we found that most of the newly introduced 4sU prodrugs were well tolerated by the cells. One derivative, the bis(4-acetyloxybenzyl) 5'-monophosphate of 4sU, was also efficiently incorporated into nascent RNA.

CHD1 controls H3.3 incorporation in adult brain chromatin to maintain metabolic homeostasis and normal lifespan.

The ATP-dependent chromatin remodeling factor CHD1 is essential for the assembly of variant histone H3.3 into paternal chromatin during sperm chromatin remodeling in fertilized eggs. It remains unclear, however, if CHD1 has a similar role in normal diploid cells. Using a specifically tailored quantitative mass spectrometry approach, we show that Chd1 disruption results in reduced H3.3 levels in heads of Chd1 mutant flies. Chd1 deletion perturbs brain chromatin structure in a similar way as H3.3 deletion and leads to global de-repression of transcription. The physiological consequences are reduced food intake, metabolic alterations, and shortened lifespan. Notably, brain-specific CHD1 expression rescues these phenotypes. We further demonstrate a strong genetic interaction between Chd1 and H3.3 chaperone Hira. Thus, our findings establish CHD1 as a factor required for the assembly of H3.3-containing chromatin in adult cells and suggest a crucial role for CHD1 in the brain as a regulator of organismal health and longevity.

Spt6 is a maintenance factor for centromeric CENP-A.

Replication and transcription of genomic DNA requires partial disassembly of nucleosomes to allow progression of polymerases. This presents both an opportunity to remodel the underlying chromatin and a danger of losing epigenetic information. Centromeric transcription is required for stable incorporation of the centromere-specific histone dCENP-A in M/G1 phase, which depends on the eviction of previously deposited H3/H3.3-placeholder nucleosomes. Here we demonstrate that the histone chaperone and transcription elongation factor Spt6 spatially and temporarily coincides with centromeric transcription and prevents the loss of old CENP-A nucleosomes in both Drosophila and human cells. Spt6 binds directly to dCENP-A and dCENP-A mutants carrying phosphomimetic residues alleviate this association. Retention of phosphomimetic dCENP-A mutants is reduced relative to wildtype, while non-phosphorylatable dCENP-A retention is increased and accumulates at the centromere. We conclude that Spt6 acts as a conserved CENP-A maintenance factor that ensures long-term stability of epigenetic centromere identity during transcription-mediated chromatin remodeling.

Thioguanosine Conversion Enables mRNA-Lifetime Evaluation by RNA Sequencing Using Double Metabolic Labeling (TUC-seq DUAL).

Temporal information about cellular RNA populations is essential to understand the functional roles of RNA. We have developed the hydrazine/NH4 Cl/OsO4 -based conversion of 6-thioguanosine (6sG) into A', where A' constitutes a 6-hydrazino purine derivative. A' retains the Watson-Crick base-pair mode and is efficiently decoded as adenosine in primer extension assays and in RNA sequencing. Because 6sG is applicable to metabolic labeling of freshly synthesized RNA and because the conversion chemistry is fully compatible with the conversion of the frequently used metabolic label 4-thiouridine (4sU) into C, the combination of both modified nucleosides in dual-labeling setups enables high accuracy measurements of RNA decay. This approach, termed TUC-seq DUAL, uses the two modified nucleosides in subsequent pulses and their simultaneous detection, enabling mRNA-lifetime evaluation with unprecedented precision.

Thiouridine-to-Cytidine Conversion Sequencing (TUC-Seq) to Measure mRNA Transcription and Degradation Rates.

The study of RNA dynamics, specifically RNA transcription and decay rates, has gained increasing attention in recent years because various mechanisms have been discovered that affect mRNA half-life, thereby ultimately controlling protein output. Therefore, there is a need for methods enabling minimally invasive, simple and high-throughput determination of RNA stability that can be applied to determine RNA transcription and decay rates in cells and organisms. We have recently developed a protocol which we named TUC-seq to directly distinguish newly synthesized transcripts from the preexisting pool of transcripts by metabolic labeling of nascent RNAs with 4-thiouridine (4sU) followed by osmium tetroxide-mediated conversion of 4sU to cytidine (C) and direct sequencing. In contrast to other related methods (SLAM-seq, TimeLapse-seq), TUC-seq converts 4sU to a native C instead of an alkylated or otherwise modified nucleoside derivative. TUC-seq can be applied to any cell type that is amenable to 4sU labeling. By employing different labeling strategies (pulse or pulse-chase labeling), it is suitable for a broad field of applications and provides a fast and highly efficient means to determine mRNA transcription and decay rates.

DNA Sequence Is a Major Determinant of Tetrasome Dynamics.

Eukaryotic genomes are hierarchically organized into protein-DNA assemblies for compaction into the nucleus. Nucleosomes, with the (H3-H4)2 tetrasome as a likely intermediate, are highly dynamic in nature by way of several different mechanisms. We have recently shown that tetrasomes spontaneously change the direction of their DNA wrapping between left- and right-handed conformations, which may prevent torque buildup in chromatin during active transcription or replication. DNA sequence has been shown to strongly affect nucleosome positioning throughout chromatin. It is not known, however, whether DNA sequence also impacts the dynamic properties of tetrasomes. To address this question, we examined tetrasomes assembled on a high-affinity DNA sequence using freely orbiting magnetic tweezers. In this context, we also studied the effects of mono- and divalent salts on the flipping dynamics. We found that neither DNA sequence nor altered buffer conditions affect overall tetrasome structure. In contrast, tetrasomes bound to high-affinity DNA sequences showed significantly altered flipping kinetics, predominantly via a reduction in the lifetime of the canonical state of left-handed wrapping. Increased mono- and divalent salt concentrations counteracted this behavior. Thus, our study indicates that high-affinity DNA sequences impact not only the positioning of the nucleosome but that they also endow the subnucleosomal tetrasome with enhanced conformational plasticity. This may provide a means to prevent histone loss upon exposure to torsional stress, thereby contributing to the integrity of chromatin at high-affinity sites.

Phosphorylation of Drosophila CENP-A on serine 20 regulates protein turn-over and centromere-specific loading.

Centromeres are specialized chromosomal regions epigenetically defined by the presence of the histone H3 variant CENP-A. CENP-A is required for kinetochore formation which is essential for chromosome segregation during mitosis. Spatial restriction of CENP-A to the centromere is tightly controlled. Its overexpression results in ectopic incorporation and the formation of potentially deleterious neocentromeres in yeast, flies and in various human cancers. While the contribution of posttranslational modifications of CENP-A to these processes has been studied in yeast and mammals to some extent, very little is known about Drosophila melanogaster. Here, we show that CENP-A is phosphorylated at serine 20 (S20) by casein kinase II and that in mitotic cells, the phosphorylated form is enriched on chromatin. Importantly, our results reveal that S20 phosphorylation regulates the turn-over of prenucleosomal CENP-A by the SCFPpa-proteasome pathway and that phosphorylation promotes removal of CENP-A from ectopic but not from centromeric sites in chromatin. We provide multiple lines of evidence for a crucial role of S20 phosphorylation in controlling restricted incorporation of CENP-A into centromeric chromatin in flies. Modulation of the phosphorylation state of S20 may provide the cells with a means to fine-tune CENP-A levels in order to prevent deleterious loading to extra-centromeric sites.

Publisher Correction: RNA modification: Getting a hold on cytosine methylation in mRNA.

In the version of this article initially published, the author order (Liu, J., Huang, T., Chen, W., Gu, N. & Zhang, R.) for the linked Resource (ref. 7) was incorrect, and all named citations of that reference in the text ("Liu et al.") were incorrect. The correct author order is "Huang, T., Chen, W., Liu, J., Gu, N. & Zhang, R." and all named citations should reflect that ("Huang et al."). The error has been corrected in the HTML and PDF versions of the article.

Role for Chromatin Remodeling Factor Chd1 in Learning and Memory.

Precise temporal and spatial regulation of gene expression in the brain is a prerequisite for cognitive processes such as learning and memory. Epigenetic mechanisms that modulate the chromatin structure have emerged as important regulators in this context. While posttranslational modification of histones or the modification of DNA bases have been examined in detail in many studies, the role of ATP-dependent chromatin remodeling factors (ChRFs) in learning- and memory-associated gene regulation has largely remained obscure. Here we present data that implicate the highly conserved chromatin assembly and remodeling factor Chd1 in memory formation and the control of immediate early gene (IEG) response in the hippocampus. We used various paradigms to assess short-and long-term memory in mice bearing a mutated Chd1 gene that gives rise to an N-terminally truncated protein. Our data demonstrate that the Chd1 mutation negatively affects long-term object recognition and short- and long-term spatial memory. We found that Chd1 regulates hippocampal expression of the IEG early growth response 1 (Egr1) and activity-regulated cytoskeleton-associated (Arc) but not cFos and brain derived neurotrophic factor (Bdnf), because the Chd1-mutation led to dysregulation of Egr1 and Arc expression in naive mice and in mice analyzed at different stages of object location memory (OLM) testing. Of note, Chd1 likely regulates Egr1 in a direct manner, because chromatin immunoprecipitation (ChIP) assays revealed enrichment of Chd1 upon stimulation at the Egr1 genomic locus in the hippocampus and in cultured cells. Together these data support a role for Chd1 as a critical regulator of molecular mechanisms governing memory-related processes, and they show that this function involves the N-terminal serine-rich region of the protein.

The dynamic RNA modification 5-methylcytosine and its emerging role as an epitranscriptomic mark.

It is a well-known fact that RNA is the target of a plethora of modifications which currently amount to over a hundred. The vast majority of these modifications was observed in the two most abundant classes of RNA, rRNA and tRNA. With the recent advance in mapping technologies, modifications have been discovered also in mRNA and in less abundant non-coding RNA species. These developments have sparked renewed interest in elucidating the nature and functions of those "epitransciptomic" modifications in RNA. N6-methyladenosine (m6 A) is the best understood and most frequent mark of mRNA with demonstrated functions ranging from pre-mRNA processing, translation, miRNA biogenesis to mRNA decay. By contrast, much less research has been conducted on 5-methylcytosine (m5C), which was detected in tRNAs and rRNAs and more recently in poly(A)RNAs. In this review, we discuss recent developments in the discovery of m5C RNA methylomes, the functions of m5C as well as the proteins installing, translating and manipulating this modification. Although our knowledge about m5C in RNA transcripts is just beginning to consolidate, it has become clear that cytosine methylation represents a powerful mechanistic strategy to regulate cellular processes on an epitranscriptomic level. This article is categorized under: RNA Processing > RNA Editing and Modification RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > tRNA Processing RNA Turnover and Surveillance > Regulation of RNA Stability.

Bisulfite Sequencing of RNA for Transcriptome-Wide Detection of 5-Methylcytosine.

A powerful method to determine the methylation status of specific cytosine residues within RNA is bisulfite sequencing. In combination with high-throughput sequencing methods cytosine methylation can be determined at nucleotide resolution on a transcriptome-wide level. Nevertheless, several critical aspects need to be considered before starting such a project. Below we describe a detailed step-by-step protocol for planning and performing a transcriptome-wide bisulfite sequencing experiment and subsequent data analysis to determine methyl-cytosine in poly(A)RNA from cells and tissues.

Superior cellular activities of azido- over amino-functionalized ligands for engineered preQ1 riboswitches in E.coli.

For this study, we utilized class-I and class-II preQ1-sensing riboswitches as model systems to decipher the structure-activity relationship of rationally designed ligand derivatives in vitro and in vivo. We found that synthetic preQ1 ligands with amino-modified side chains that protrude from the ligand-encapsulating binding pocket, and thereby potentially interact with the phosphate backbone in their protonated form, retain or even increase binding affinity for the riboswitches in vitro. They, however, led to significantly lower riboswitch activities in a reporter system in vivo in E. coli. Importantly, when we substituted the amino- by azido-modified side chains, the cellular activities of the ligands were restored for the class-I conditional gene expression system and even improved for the class-II counterpart. Kinetic analysis of ligand binding in vitro revealed enhanced on-rates for amino-modified derivatives while they were attenuated for azido-modified variants. This shows that neither high affinities nor fast on-rates are necessarily translated into efficient cellular activities. Taken together, our comprehensive study interconnects in vitro kinetics and in vitro thermodynamics of RNA-ligand binding with the ligands' in vivo performance and thereby encourages azido- rather than amino-functionalized design for enhanced cellular activity.

Modification of the histone tetramer at the H3-H3 interface impacts tetrasome conformations and dynamics.

Nucleosomes consisting of a short piece of deoxyribonucleic acid (DNA) wrapped around an octamer of histone proteins form the fundamental unit of chromatin in eukaryotes. Their role in DNA compaction comes with regulatory functions that impact essential genomic processes such as replication, transcription, and repair. The assembly of nucleosomes obeys a precise pathway in which tetramers of histones H3 and H4 bind to the DNA first to form tetrasomes, and two dimers of histones H2A and H2B are subsequently incorporated to complete the complex. As viable intermediates, we previously showed that tetrasomes can spontaneously flip between a left-handed and right-handed conformation of DNA-wrapping. To pinpoint the underlying mechanism, here we investigated the role of the H3-H3 interface for tetramer flexibility in the flipping process at the single-molecule level. Using freely orbiting magnetic tweezers, we studied the assembly and structural dynamics of individual tetrasomes modified at the cysteines close to this interaction interface by iodoacetamide (IA) in real time. While such modification did not affect the structural properties of the tetrasomes, it caused a 3-fold change in their flipping kinetics. The results indicate that the IA-modification enhances the conformational plasticity of tetrasomes. Our findings suggest that subnucleosomal dynamics may be employed by chromatin as an intrinsic and adjustable mechanism to regulate DNA supercoiling.

RNA cytosine methyltransferase Nsun3 regulates embryonic stem cell differentiation by promoting mitochondrial activity.

Chemical modifications of RNA have been attracting increasing interest because of their impact on RNA fate and function. Therefore, the characterization of enzymes catalyzing such modifications is of great importance. The RNA cytosine methyltransferase NSUN3 was recently shown to generate 5-methylcytosine in the anticodon loop of mitochondrial tRNAMet. Further oxidation of this position is required for normal mitochondrial translation and function in human somatic cells. Because embryonic stem cells (ESCs) are less dependent on oxidative phosphorylation than somatic cells, we examined the effects of catalytic inactivation of Nsun3 on self-renewal and differentiation potential of murine ESCs. We demonstrate that Nsun3-mutant cells show strongly reduced mt-tRNAMet methylation and formylation as well as reduced mitochondrial translation and respiration. Despite the lower dependence of ESCs on mitochondrial activity, proliferation of mutant cells was reduced, while pluripotency marker gene expression was not affected. By contrast, ESC differentiation was skewed towards the meso- and endoderm lineages at the expense of neuroectoderm. Wnt3 was overexpressed in early differentiating mutant embryoid bodies and in ESCs, suggesting that impaired mitochondrial function disturbs normal differentiation programs by interfering with cellular signalling pathways. Interestingly, basal levels of reactive oxygen species (ROS) were not altered in ESCs, but Nsun3 inactivation attenuated induction of mitochondrial ROS upon stress, which may affect gene expression programs upon differentiation. Our findings not only characterize Nsun3 as an important regulator of stem cell fate but also provide a model system to study the still incompletely understood interplay of mitochondrial function with stem cell pluripotency and differentiation.