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.

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.

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.

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.

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.

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.

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.

A novel role for the histone acetyltransferase Hat1 in the CENP-A/CID assembly pathway in Drosophila melanogaster.

The incorporation of CENP-A into centromeric chromatin is an essential prerequisite for kinetochore formation. Yet, the molecular mechanisms governing this process are surprisingly divergent in different organisms. While CENP-A loading mechanisms have been studied in some detail in mammals, there are still large gaps to our understanding of CENP-A/Cid loading pathways in Drosophila. Here, we report on the characterization and delineation of at least three different CENP-A preloading complexes in Drosophila. Two complexes contain the CENP-A chaperones CAL1, FACT and/or Caf1/Rbap48. Notably, we identified a novel complex consisting of the histone acetyltransferase Hat1, Caf1 and CENP-A/H4. We show that Hat1 is required for proper CENP-A loading into chromatin, since knock-down in S2 cells leads to reduced incorporation of newly synthesized CENP-A. In addition, we demonstrate that CENP-A/Cid interacts with the HAT1 complex via an N-terminal region, which is acetylated in cytoplasmic but not in nuclear CENP-A. Since Hat1 is not responsible for acetylation of CENP-A/Cid, these results suggest a histone acetyltransferase activity-independent escort function for Hat1. Thus, our results point toward intriguing analogies between the complex processing pathways of newly synthesized CENP-A and canonical histones.

meRanTK: methylated RNA analysis ToolKit.

UNLABELLED: The significance and function of posttranscriptional cytosine methylation in poly(A)RNA attracts great interest but is still poorly understood. High-throughput sequencing of RNA treated with bisulfite (RNA-BSseq) or subjected to enrichment techniques like Aza-IP or miCLIP enables transcriptome wide studies of this particular modification at single base pair resolution. However, to date, there are no specialized software tools available for the analysis of RNA-BSseq or Aza-IP data. Therefore, we developed meRanTK, the first publicly available tool kit which addresses the special demands of high-throughput RNA cytosine methylation data analysis. It provides fast and easy to use splice-aware bisulfite sequencing read mapping, comprehensive methylation calling and identification of differentially methylated cytosines by statistical analysis of single- and multi-replicate experiments. Application of meRanTK to RNA-BSseq or Aza-IP data produces accurate results in standard compliant formats. AVAILABILITY AND IMPLEMENTATION: meRanTK, source code and test data are released under the GNU GPLv3+ license and are available at http://icbi.at/software/meRanTK/ CONTACT: dietmar.rieder@i-med.ac.at.

Synthesis, Thermodynamic Properties, and Crystal Structure of RNA Oligonucleotides Containing 5-Hydroxymethylcytosine.

5-Hydroxymethylcytosine (hm5C) is an RNA modification that has attracted significant interest because of the finding that RNA hydroxymethylation can favor mRNA translation. For insight into the mechanistic details of hm5C function to be obtained, the availability of RNAs containing this modification at defined positions that can be used for in vitro studies is highly desirable. In this work, we present an eight-step route to 5-hydroxymethylcytidine (hm5rC) phosphoramidite for solid-phase synthesis of modified RNA oligonucleotides. Furthermore, we examined the effects of hm5rC on RNA duplex stability and its impact on structure formation using the sarcin-ricin loop (SRL) motif. Thermal denaturation experiments revealed that hm5rC increases RNA duplex stability. By contrast, when cytosine within an UNCG tetraloop motif was replaced by hm5rC, the thermodynamic stability of the corresponding hairpin fold was attenuated. Importantly, incorporation of hm5rC into the SRL motif resulted in an RNA crystal structure at 0.85 Å resolution. Besides changes in the hydration pattern at the site of modification, a slight opening of the hm5rC-G pair compared to the unmodified C-G in the native structure was revealed.

Osmium-Mediated Transformation of 4-Thiouridine to Cytidine as Key To Study RNA Dynamics by Sequencing.

To understand the functional roles of RNA in the cell, it is essential to elucidate the dynamics of their production, processing and decay. A recent method for assessing mRNA dynamics is metabolic labeling with 4-thiouridine (4sU), followed by thio-selective attachment of affinity tags. Detection of labeled transcripts by affinity purification and hybridization to microarrays or by deep sequencing then reveals RNA expression levels. Here, we present a novel sequencing method (TUC-seq) that eliminates affinity purification and allows for direct assessment of 4sU-labeled RNA. It employs an OsO4 -mediated transformation to convert 4sU into cytosine. We exemplify the utility of the new method for verification of endogenous 4sU in tRNAs and for the detection of pulse-labeled mRNA of seven selected genes in mammalian cells to determine the relative abundance of the new transcripts. The results prove TUC-seq as a straight-forward and highly versatile method for studies of cellular RNA dynamics.

Differential neuronal targeting of a new and two known calcium channel ß4 subunit splice variants correlates with their regulation of gene expression.

The ß subunits of voltage-gated calcium channels regulate surface expression and gating of CaV1 and CaV2 α1 subunits and thus contribute to neuronal excitability, neurotransmitter release, and calcium-induced gene regulation. In addition, certain ß subunits are targeted into the nucleus, where they interact directly with the epigenetic machinery. Whereas their involvement in this multitude of functions is reflected by a great molecular heterogeneity of ß isoforms derived from four genes and abundant alternative splicing, little is known about the roles of individual ß variants in specific neuronal functions. In the present study, an alternatively spliced ß4 subunit lacking the variable N terminus (ß4e) is identified. It is highly expressed in mouse cerebellum and cultured cerebellar granule cells (CGCs) and modulates P/Q-type calcium currents in tsA201 cells and CaV2.1 surface expression in neurons. Compared with the other two known full-length ß4 variants (ß4a and ß4b), ß4e is most abundantly expressed in the distal axon, but lacks nuclear-targeting properties. To determine the importance of nuclear targeting of ß4 subunits for transcriptional regulation, we performed whole-genome expression profiling of CGCs from lethargic (ß4-null) mice individually reconstituted with ß4a, ß4b, and ß4e. Notably, the number of genes regulated by each ß4 splice variant correlated with the rank order of their nuclear-targeting properties (ß4b > ß4a > ß4e). Together, these findings support isoform-specific functions of ß4 splice variants in neurons, with ß4b playing a dual role in channel modulation and gene regulation, whereas the newly detected ß4e variant serves exclusively in calcium-channel-dependent functions.

Dynamics of nucleosomal structures measured by high-speed atomic force microscopy.

The accessibility of DNA is determined by the number, position, and stability of nucleosomes, complexes consisting of a core of 8 histone proteins with DNA wrapped around it. Since the structure and dynamics of nucleosomes affects essential cellular processes, they are the subject of many current studies. Here, high-speed atomic force microscopy is used to visualize dynamic processes in nucleosomes and tetrasomes (subnucleosomal structures that contain 4 rather than 8 histones in the protein core). Nucleosomes can spontaneously disassemble in a process (at a 1 second timescale). For tetrasomes, multiple dynamic phenomena are observed. For example, during disassembly the formation of a DNA loop (∼25 nm in length) is seen, which remains stable for several minutes. For intact tetrasomes, dynamics in the form of sliding and reversible hopping between stable positions along the DNA are observed. The data emphasize that tetrasomes are not merely static objects but highly dynamic. Since tetrasomes (in contrast to nucleosomes) can stay on the DNA during transcription, the observed tetrasome dynamics is relevant for an understanding of the nucleosomal dynamics during transcription. These results illustrate the diversity of nucleosome dynamics and demonstrate the ability of high speed AFM to characterize protein-DNA interactions.

Long non-coding RNAs as targets for cytosine methylation.

Post-synthetic modifications of nucleic acids have long been known to affect their functional and structural properties. For instance, numerous different chemical modifications modulate the structural organization, stability or translation efficiency of tRNAs and rRNAs. In contrast, little is known about modifications of poly(A)RNAs. Here, we demonstrate for the first time that the two well-studied regulatory long non-coding RNAs HOTAIR and XIST are targets of site-specific cytosine methylation. In both XIST and HOTAIR, we found methylated cytosines located within or near functionally important regions that are known to mediate interaction with chromatin-associated protein complexes. We show that cytosine methylation in the XIST A structure strongly affects binding to the chromatin-modifying complex PRC2 in vitro. These results suggest that cytosine methylation may serve as a general strategy to regulate the function of long non-coding RNAs.

The chromodomains of CHD1 are critical for enzymatic activity but less important for chromatin localization.

The molecular motor protein CHD1 has been implicated in the regulation of transcription and in the transcription-independent genome-wide incorporation of H3.3 into paternal chromatin in Drosophila melanogaster. A key feature of CHD1 is the presence of two chromodomains, which can bind to histone H3 methylated at lysine 4 and thus might serve to recruit and/or maintain CHD1 at the chromatin. Here, we describe genetic and biochemical approaches to the study of the Drosophila CHD1 chromodomains. We found that overall localization of CHD1 on polytene chromosomes does not appreciably change in chromodomain-mutant flies. In contrast, the chromodomains are important for transcription-independent activities of CHD1 during early embryonic development as well as for transcriptional regulation of several heat shock genes. However, neither CHD1 nor its chromodomains are needed for RNA polymerase II localization and H3K4 methylation but loss of CHD1 decreases transcription-induced histone eviction at the Hsp70 gene in vivo. Chromodomain mutations negatively affect the chromatin assembly activities of CHD1 in vitro, and they appear to be involved in linking the ATP-dependent motor to the chromatin assembly function of CHD1.

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.

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.

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.