Do RNA Chemical Modifications Encode a Second Layer of Gene Regulation?
Do RNA chemical modifications such as m6A, m5C, pseudouridine, and A-to-I editing form a higher-order regulatory system that encodes transcript fate, ribosome specialization, and chromatin–RNA feedback loops? Explore three testable hypotheses on epitranscriptomic combinatorial codes, m6A–chromatin crosstalk in cancer apoptosis resistance, and rRNA modification-driven ribosome heterogeneity shaping stress-responsive translation.
Scientific Hypothesis Generation
Do RNA Chemical Modifications Encode a Second Layer of Gene Regulation?
Over 170 chemical modifications decorate RNA molecules in human cells, yet most remain poorly characterised in terms of function, coordination, and disease relevance. The epitranscriptome sits at the intersection of chromatin biology, translational control, and ribosome biology, raising questions that existing frameworks cannot answer.
Hypothesis 1
m6A-chromatin crosstalk operates as a bidirectional feedback loop that, when decoupled, confers apoptosis resistance in solid tumours
The Gap
Histone modifications such as H3K36me3 are known to recruit the METTL3/METTL14/WTAP methyltransferase complex to chromatin, guiding co-transcriptional m6A deposition on nascent pre-mRNA. Separately, m6A reader proteins like YTHDC1 have been shown to recruit chromatin remodellers. However, no study has mapped the complete bidirectional circuit connecting specific histone marks, m6A site selection, reader-mediated chromatin feedback, and downstream phenotypic consequences in a single disease model.
The apoptosis context is particularly underexplored: apoptotic stimulus-specific m6A sites on pro- and antiapoptotic transcripts remain largely uncharacterised, and it is unclear whether disruption of the chromatin-to-m6A-to-chromatin loop is causally linked to apoptosis evasion.
The Claim
We hypothesise that H3K36me3-guided m6A deposition on pro-apoptotic transcripts (such as BAX and BIM mRNAs) creates a positive feedback loop in which YTHDC1-mediated recruitment of the SWI/SNF chromatin remodelling complex to the corresponding gene loci sustains an open chromatin state, thereby maintaining transcription of these same pro-apoptotic genes.
In tumour cells that acquire resistance to apoptosis, this loop is broken at the METTL3 recruitment step. Loss of H3K36me3 at pro-apoptotic loci (through SETD2 mutation or EZH2 overexpression) reduces m6A deposition, silences YTHDC1-mediated chromatin feedback, and shifts the chromatin to a repressive state at these genes. The net result is selective transcriptional silencing of pro-apoptotic mRNAs while antiapoptotic transcripts (BCL2, MCL1) remain m6A-marked and translationally active.
This model predicts that SETD2-mutant tumours should show a coordinated loss of both H3K36me3 and m6A at pro-apoptotic gene bodies, measurable at single-locus resolution.
Why It's Testable Now
CUT&Tag for histone marks and m6A-seq (or miCLIP) can now be performed from the same cell population, enabling locus-resolved comparison of H3K36me3 and m6A co-occupancy. CRISPR-mediated SETD2 knockout panels in isogenic cell lines provide the perturbation model, while nanopore direct RNA sequencing detects m6A changes on native transcripts without antibody bias.
The Intriguing Outcome
If confirmed, this would establish a mechanistic link between a chromatin lesion (SETD2 loss) and epitranscriptomic silencing of apoptotic programmes, suggesting that m6A-restoring strategies (such as METTL3-tethering via dCas13) could re-sensitise resistant tumours to apoptosis-inducing therapies.
It would also reframe SETD2 mutations in clear cell renal carcinoma, lung adenocarcinoma, and other cancers as epitranscriptomic driver events rather than purely epigenomic ones.
Thesis Entry Points
- Perform paired CUT&Tag (H3K36me3) and miCLIP-seq in SETD2-wildtype vs. SETD2-knockout isogenic renal carcinoma cells; quantify locus-specific m6A loss at BAX, BIM, PUMA gene bodies using a differential methylation pipeline, with RT-qPCR and Western blot validation of transcript and protein output.
- Use YTHDC1-CLIP in the same cell pairs to map changes in reader occupancy at pro-apoptotic loci; assess SWI/SNF subunit (SMARCA2/SMARCA4) recruitment by ChIP-qPCR at those loci before and after SETD2 ablation, measuring ATAC-seq accessibility as a readout of chromatin state.
- Treat SETD2-null cells with a dCas13-METTL3 fusion targeted to BAX/BIM 3'UTRs; measure restored m6A by SELECT assay, monitor apoptosis re-sensitisation by Annexin V/PI flow cytometry after cisplatin or TRAIL challenge, and quantify chromatin reopening by ATAC-seq at the target loci.
Novelty Signal
Emerging: m6A-chromatin crosstalk is actively studied, but the bidirectional feedback model linking a specific histone writer (SETD2) to apoptosis resistance through an epitranscriptomic relay has not been tested in an integrated, locus-resolved framework.
Hypothesis 2
Co-occurring RNA modifications on single mRNA molecules form a combinatorial code that determines transcript fate
The Gap
Individual RNA modifications (m6A, m5C, pseudouridine, A-to-I editing) have each been linked to changes in RNA stability, splicing, translation efficiency, and protein binding. Yet nearly all functional studies examine one modification type in isolation, using bulk sequencing that averages signals across millions of transcript copies.
Whether two or more distinct modifications co-occur on the same physical molecule, and if so, whether their combinations produce non-additive functional effects, remains unknown. The concept of a combinatorial epitranscriptomic code, analogous to the histone code, has been proposed but lacks direct experimental evidence at the single-molecule level.
The Claim
We hypothesise that individual mRNA molecules carry stereotyped combinations of m6A and pseudouridine marks, and that these combinations, rather than individual marks alone, determine whether a given transcript is routed to translation, storage in P-bodies, or degradation via the CCR4-NOT complex.
Specifically, we predict that mRNAs bearing m6A in the 3'UTR together with pseudouridine in the coding sequence are preferentially stabilised and translated, whereas transcripts carrying m6A without co-occurring pseudouridine are preferentially targeted for YTHDF2-mediated decay. This asymmetry arises because pseudouridine-dependent structural stabilisation of the coding region protects the transcript from endonucleolytic cleavage, allowing YTHDF1-mediated translational enhancement at the m6A site to dominate over the competing decay pathway.
The prediction is that disrupting pseudouridine synthase PUS7, which modifies coding-region uridines in a subset of transcripts, will shift the balance from translation to decay for m6A-positive transcripts, measurable as a decrease in polysome association and an increase in CCR4-NOT co-immunoprecipitation.
Why It's Testable Now
Oxford Nanopore direct RNA sequencing can simultaneously detect m6A and pseudouridine on individual native RNA molecules without antibody enrichment. Machine learning classifiers (m6Anet, Nanocompore) trained on known modification signatures enable multi-mark co-detection from a single sequencing run, making single-molecule combinatorial profiling feasible for the first time.
The Intriguing Outcome
Demonstration of a combinatorial code would fundamentally expand the regulatory vocabulary of the epitranscriptome from individual marks to mark combinations, each with distinct downstream consequences. It would also explain contradictory findings in the m6A field, where the same modification has been associated with both transcript stabilisation and decay depending on the experimental system.
Therapeutically, it would suggest that modulating one modification type (e.g., inhibiting PUS7) could redirect the fate of a specific class of m6A-bearing transcripts, opening a new axis for epitranscriptomic drug design.
Thesis Entry Points
- Perform nanopore direct RNA sequencing on poly(A)-selected RNA from HEK293T cells; apply m6Anet and Nanocompore to co-detect m6A and pseudouridine on individual reads; cluster transcripts by modification pattern and correlate each cluster with ribosome profiling (Ribo-seq) translation efficiency and RNA half-life (SLAM-seq) data from the same cell line.
- Generate PUS7-knockout HEK293T cells by CRISPR; repeat nanopore sequencing to confirm loss of coding-region pseudouridine on target transcripts; measure changes in polysome distribution of m6A-positive transcripts by sucrose gradient fractionation followed by RT-qPCR, and assess CCR4-NOT recruitment by CNOT1 co-immunoprecipitation and RT-qPCR of bound mRNAs.
- Reconstitute the combinatorial effect in vitro using synthetic reporter mRNAs bearing defined m6A and pseudouridine patterns; transfect into cells and measure translation output (NanoLuc luciferase) and mRNA decay kinetics (actinomycin D chase, RT-qPCR) as a function of modification combination, with YTHDF1/YTHDF2 knockdown controls to dissect reader specificity.
Novelty Signal
Frontier: No published study has simultaneously mapped two or more distinct RNA modification types on single molecules and correlated the combinatorial patterns with functional transcript fate outcomes.
Hypothesis 3
Differential rRNA modification patterns generate functionally specialised ribosome subpopulations with distinct translational outputs
The Gap
Eukaryotic ribosomes carry over 100 rRNA modifications, predominantly 2'-O-methylation (Nm) and pseudouridylation, installed by snoRNA-guided complexes during ribosome biogenesis. These modifications cluster at functionally critical sites: the peptidyl transferase centre, the decoding centre, and inter-subunit bridges. Recent nanopore sequencing of yeast rRNA revealed that not all ribosomes are identically modified, with concerted modification patterns at functional centres.
However, whether these differentially modified ribosome populations translate distinct mRNA pools, and whether the modification pattern is actively regulated in response to cellular signals, remains unanswered.
The Claim
We hypothesise that snoRNA expression changes during cellular stress (e.g., hypoxia or nutrient deprivation) alter the rRNA modification landscape at the decoding centre (specifically at positions A1832 and C1843 in human 18S rRNA), generating a ribosome subpopulation with altered codon-reading fidelity. These "stress ribosomes" preferentially translate mRNAs enriched in non-optimal codons, which include many stress-response and autophagy-related transcripts (ATG genes, HIF1A, BNIP3).
The mechanism involves reduced 2'-O-methylation at A1832 (guided by SNORD42A), which loosens the geometry of the decoding centre and increases tolerance for wobble base pairing. This shifts the translational output toward transcripts that rely on wobble-decoded codons, effectively creating a stress-adaptive translational programme without transcriptional changes.
The prediction is that SNORD42A knockdown under normal growth conditions should phenocopy the stress-ribosome translational signature, measurable as increased translation efficiency of ATG/HIF1A transcripts by Ribo-seq, even in the absence of hypoxic stimulation.
Why It's Testable Now
Nanopore direct rRNA sequencing, as demonstrated in yeast by Bailey et al. (2022), can profile modification status at dozens of sites on individual rRNA molecules. Combined with Ribo-seq and polysome profiling, the translational consequences of modification heterogeneity can be linked to specific ribosome subpopulations in human cells.
The Intriguing Outcome
If validated, this would establish rRNA modification as a tuneable switch for translational reprogramming during stress, operating independently of mRNA transcription. It would add a new regulatory layer to the cellular stress response and place ribosome heterogeneity at the centre of translational control.
Clinically, it could explain how tumour cells in hypoxic microenvironments sustain translation of survival genes despite global translational suppression, and suggest snoRNA-targeted antisense oligonucleotides as a strategy to disrupt this adaptive programme.
Thesis Entry Points
- Subject HeLa or A549 cells to hypoxia (1% O2, 24 h); perform nanopore direct rRNA sequencing on purified 40S subunits; compare 2'-O-methylation levels at A1832 and C1843 between normoxic and hypoxic conditions using Nanocompore or ELIGOS2, with RiboMeth-seq as orthogonal validation.
- Generate stable SNORD42A-knockdown cell lines using CRISPRi targeting the snoRNA host gene promoter; confirm reduced 2'-O-methylation at A1832 by RiboMeth-seq; perform Ribo-seq under normoxic conditions and compare translational efficiency of stress-response transcripts (ATG5, ATG7, HIF1A, BNIP3) against a scrambled-guide control, using differential translation efficiency analysis (deltaTE).
- Fractionate polysome gradients from hypoxic cells; isolate monosome vs. heavy polysome fractions; perform nanopore rRNA sequencing on each fraction to determine whether hypo-modified ribosomes are enriched in polysomes translating stress-response mRNAs (identified by polysome-fraction RT-qPCR), testing the prediction that modification state and mRNA identity co-segregate.
Novelty Signal
Open field: Ribosome heterogeneity driven by rRNA modification has been documented in yeast, but the causal link between stress-induced snoRNA changes, site-specific rRNA hypo-modification, and selective translation of stress-response mRNAs has not been tested in mammalian cells.
Frequently asked questions
What are RNA chemical modifications?
RNA chemical modifications are covalent additions to RNA nucleotides that alter their structure and function without changing the underlying sequence. Over 170 types have been identified, including m6A, m5C, pseudouridine, and A-to-I editing. They regulate RNA stability, splicing, translation, and protein interactions.
What is the epitranscriptome?
The epitranscriptome refers to the complete set of biochemical modifications present on RNA molecules in a cell. It represents a post-transcriptional regulatory layer analogous to the epigenome on DNA. Epitranscriptomic marks are dynamic and reversible, installed and removed by dedicated writer and eraser enzymes.
How does m6A modification crosstalk with chromatin?
The m6A methyltransferase complex (METTL3/METTL14/WTAP) is recruited to chromatin co-transcriptionally, where specific histone marks such as H3K36me3 guide m6A deposition on nascent RNA. Conversely, m6A-marked transcripts can recruit chromatin remodellers back to DNA, creating a bidirectional feedback loop between histone and RNA modification states.
Can multiple RNA modifications co-occur on the same transcript?
Yes. Recent nanopore direct RNA sequencing studies show that individual transcript molecules can carry multiple distinct modifications simultaneously. Whether these co-occurring marks form a combinatorial code with predictable functional outputs is an active area of investigation.
What is ribosome heterogeneity?
Ribosome heterogeneity refers to the existence of structurally and functionally distinct populations of ribosomes within a single cell. Differences in rRNA modification patterns, ribosomal protein composition, and associated factors may produce ribosomes with different translational preferences or efficiencies.
How does BioSkepsis generate scientific hypotheses?
BioSkepsis synthesises peer-reviewed literature into structured research threads, identifies knowledge gaps and mechanistic unknowns, and formulates falsifiable hypotheses grounded in the existing evidence base. Each hypothesis specifies molecular entities, testable predictions, and experimental entry points.
What technologies enable single-molecule epitranscriptomics?
Oxford Nanopore direct RNA sequencing reads native RNA molecules without reverse transcription or amplification, preserving chemical modifications as electrical signal signatures. Combined with machine learning tools such as m6Anet, it enables detection of modification status at individual sites across single molecules.
Generate Hypotheses from Peer-Reviewed Literature
BioSkepsis synthesises published research into structured, testable hypotheses. Identify knowledge gaps, plan experiments, and accelerate your next grant application.
Start freeSources and further reading
- Lee H, Park YJ, Seo PJ. A new epigenetic crosstalk: chemical modification information flow. Genomics Genet Genom. 2023;1(1):e202200033. doi:10.1002/ggn2.202200033. PMID: 37829606
- Wang Y, Huang H, Chen J, Weng H. Crosstalk between histone/DNA modifications and RNA N6-methyladenosine modification. Curr Opin Genet Dev. 2024;86:102205. doi:10.1016/j.gde.2024.102205. PMID: 38759337
- Bailey AD, Talkish J, Ding H, et al. Concerted modification of nucleotides at functional centers of the ribosome revealed by single-molecule RNA modification profiling. eLife. 2022;11:e76562. doi:10.7554/eLife.76562. PMID: 35175196
- Xu Z, Xie T, Sui X, et al. Crosstalk between histone and m6A modifications and emerging roles of m6A RNA methylation. Front Genet. 2022;13:908289. doi:10.3389/fgene.2022.908289. PMID: 35783291
- Haussmann IU, Bodi Z, Sanchez-Moran E, et al. m6A potentiates Sxl alternative pre-mRNA splicing for robust Drosophila sex determination. Nature. 2016;540:301-304. PMID: 27919081
- Schwartz S, Bernstein DA, Mumbach MR, et al. Transcriptome-wide mapping reveals widespread dynamic-regulated pseudouridylation of ncRNA and mRNA. Cell. 2014;159(1):148-162. PMID: 25219674
- Huang H, Weng H, Zhou K, et al. Histone H3 trimethylation at lysine 36 guides m6A RNA modification co-transcriptionally. Nature. 2019;567:414-419. PMID: 30867593
- Hendra C, Pratanwanich PN, Poh YS, et al. Detection of m6A from direct RNA sequencing using a multiple instance learning framework. Nat Methods. 2022;19:1590-1598. PMID: 36357692