Systematic evaluation of parameters in RNA bisulfite sequencing data generation and analysis

Splice_sim: a nucleotide conversion-enabled RNA-seq simulation and evaluation framework

Nucleotide conversion RNA sequencing techniques interrogate chemical RNA modifications in cellular transcripts, resulting in mismatch-containing reads. Biases in mapping the resulting reads to reference genomes remain poorly understood. We present splice_sim, a splice-aware RNA-seq simulation and evaluation pipeline that introduces user-defined nucleotide conversions at set frequencies, creates mixture models of converted and unconverted reads, and calculates mapping accuracies per genomic annotation. By simulating nucleotide conversion RNA-seq datasets under realistic experimental conditions, including metabolic RNA labeling and RNA bisulfite sequencing, we measure mapping accuracies of state-of-the-art spliced-read mappers for mouse and human transcripts and derive strategies to prevent biases in the data interpretation.

Prediction of m6A and m5C at single-molecule resolution reveals a transcriptome-wide co-occurrence of RNA modifications

The epitranscriptome embodies many new and largely unexplored functions of RNA. A significant roadblock hindering progress in epitranscriptomics is the identification of more than one modification in individual transcript molecules. We address this with CHEUI (CH3 (methylation) Estimation Using Ionic current). CHEUI predicts N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in individual molecules from the same sample, the stoichiometry at transcript reference sites, and differential methylation between any two conditions. CHEUI processes observed and expected nanopore direct RNA sequencing signals to achieve high single-molecule, transcript-site, and stoichiometry accuracies in multiple tests using synthetic RNA standards and cell line data. CHEUI's capability to identify two modification types in the same sample reveals a co-occurrence of m6A and m5C in individual mRNAs in cell line and tissue transcriptomes. CHEUI provides new avenues to discover and study the function of the epitranscriptome.

Dynamics of RNA m5C modification during brain development

Post-transcriptional RNA modifications have been recognized as key regulators of neuronal differentiation and synapse development in the mammalian brain. While distinct sets of 5-methylcytosine (m⁵C) modified mRNAs have been detected in neuronal cells and brain tissues, no study has been performed to characterize methylated mRNA profiles in the developing brain. Here, together with regular RNA-seq, we performed transcriptome-wide bisulfite sequencing to compare RNA cytosine methylation patterns in neural stem cells (NSCs), cortical neuronal cultures, and brain tissues at three postnatal stages. Among 501 m⁵C sites identified, approximately 6% are consistently methylated across all five conditions. Compared to m⁵C sites identified in NSCs, 96% of them were hypermethylated in neurons and enriched for genes involved in positive transcriptional regulation and axon extension. In addition, brains at the early postnatal stage demonstrated substantial changes in both RNA cytosine methylation and gene expression of RNA cytosine methylation readers, writers, and erasers. Furthermore, differentially methylated transcripts were significantly enriched for genes regulating synaptic plasticity. Altogether, this study provides a brain epitranscriptomic dataset as a new resource and lays the foundation for further investigations into the role of RNA cytosine methylation during brain development.

Quantification of substoichiometric modification reveals global tsRNA hypomodification, preferences for angiogenin-mediated tRNA cleavage, and idiosyncratic epitranscriptomes of human neuronal cell-lines

Modification of tRNA is an integral part of the epitranscriptome with a particularly pronounced potential to generate diversity in RNA expression. Eukaryotic tRNA contains modifications in up to 20% of their nucleotides, but not all sites are always fully modified. Combinations and permutations of partially modified sites in tRNAs can generate a plethora of tRNA isoforms, termed modivariants. Here, we investigate the stoichiometry of incompletely modified sites in tRNAs from human cell lines for their information content. Using a panel of RNA modification mapping methods, we assess the stoichiometry of sites that contain the modifications 5-methylcytidine (m⁵C), 2'-O-ribose methylation (Nm), 3-methylcytidine (m³C), 7-methylguanosine (m⁷G), and Dihydrouridine (D). We discovered that up to 75% of sites can be incompletely modified and that the differential modification status of a cellular tRNA population holds information that allows to discriminate e.g. different cell lines. As a further aspect, we investigated potential causal connectivity between tRNA modification and its processing into tRNA fragments (tiRNAs and tRFs). Upon exposure of cultured living cells to cell-penetrating angiogenin, the modification patterns of the corresponding RNA populations was changed. Importantly, we also found that tsRNAs were significantly less modified than their parent tRNAs at numerous sites, suggesting that tsRNAs might derive chiefly from hypomodified tRNAs.

Folate regulates RNA m5C modification and translation in neural stem cells

BACKGROUND

Folate is an essential B-group vitamin and a key methyl donor with important biological functions including DNA methylation regulation. Normal neurodevelopment and physiology are sensitive to the cellular folate levels. Either deficiency or excess of folate may lead to neurological disorders. Recently, folate has been linked to tRNA cytosine-5 methylation (m⁵C) and translation in mammalian mitochondria. However, the influence of folate intake on neuronal mRNA m⁵C modification and translation remains largely unknown. Here, we provide transcriptome-wide landscapes of m⁵C modification in poly(A)-enriched RNAs together with mRNA transcription and translation profiles for mouse neural stem cells (NSCs) cultured in three different concentrations of folate.

RESULTS

NSCs cultured in three different concentrations of folate showed distinct mRNA methylation profiles. Despite uncovering only a few differentially expressed genes, hundreds of differentially translated genes were identified in NSCs with folate deficiency or supplementation. The differentially translated genes induced by low folate are associated with cytoplasmic translation and mitochondrial function, while the differentially translated genes induced by high folate are associated with increased neural stem cell proliferation. Interestingly, compared to total mRNAs, polysome mRNAs contained high levels of m⁵C. Furthermore, an integrative analysis indicated a transcript-specific relationship between RNA m⁵C methylation and mRNA translation efficiency.

CONCLUSIONS

Altogether, our study reports a transcriptome-wide influence of folate on mRNA m⁵C methylation and translation in NSCs and reveals a potential link between mRNA m⁵C methylation and mRNA translation.

Neuronal Depolarization Induced RNA m5C Methylation Changes in Mouse Cortical Neurons

Neuronal activity is accomplished via substantial changes in gene expression, which may be accompanied by post-transcriptional modifications including RNA cytosine-5 methylation (m5C). Despite several reports on the transcriptome profiling of activated neurons, the dynamics of neuronal mRNA m5C modification in response to environmental stimuli has not been explored. Here, we provide transcriptome-wide maps of m5C modification, together with gene expression profiles, for mouse cortical neurons at 0 h, 2 h, and 6 h upon membrane depolarization. Thousands of differentially expressed genes (DEGs) were identified during the neuronal depolarization process. In stimulated neurons, the majority of early response genes were found to serve as expression regulators of late response genes, which are involved in signaling pathways and diverse synaptic functions. With RNA bisulfite sequencing data, a union set of 439 m5C sites was identified with high confidence, and approximately 30% of them were shared by neurons at all three time points. Interestingly, over 41% of the m5C sites showed increased methylation upon neuronal activation and were enriched in transcripts coding for proteins with synaptic functions. In addition, a modest negative correlation was observed between RNA expression and methylation. In summary, our study provided dynamic transcriptome-wide landscapes of RNA m5C methylation in neurons, and revealed that mRNA m5C methylation is associated with the regulation of gene expression.