Abstract
Whether individual transcripts carry intrinsic features that predetermine their response to perturbations is unknown. Here we used nanopore direct RNA sequencing of male mouse dorsal root ganglia (DRG) to simultaneously profile N6-methyladenosine (m6A) modifications, poly(A) tail dynamics, and full-length isoform identity from mice treated with bortezomib, a proteasome inhibitor that causes painful peripheral neuropathy. Machine learning revealed that transcript-intrinsic features predetermine the magnitude of perturbation-induced m6A loss (R² = 0.983). Expression level contributed just 2.6% of predictive importance. Bortezomib removed a fixed ∼73.5% fraction of m6A marks, meaning absolute loss scaled linearly with baseline density and a transcript’s epitranscriptomic fate was encoded in its architecture before drug exposure. Unsupervised clustering identified four response programs where the dominant m6A erosion cluster enriched for oxidative phosphorylation (OXPHOS, p = 1.0 × 10⁻¹⁷) and proteasome (p = 2.8 × 10⁻¹⁰) genes, recapitulating bortezomib’s established mechanisms without prior biological knowledge. Isoform-resolved analysis uncovered m6A remodeling patterns suggesting post-transcriptional regulation of glycolytic and OXPHOS genes, and Western blot confirmed protein-level suppression of OXPHOS components. Integration with single-nuclei sequencing showed sensory neurons carried 2.2-fold greater m6A loss burden than non-neuronal cells, a direct consequence of architectural determinism applied to cell-type-specific transcriptomes. These findings establish that epitranscriptomic bortezomib response is predetermined by transcript architecture, with pathway specificity and cell-type vulnerability emerging as downstream consequences of intrinsic RNA structure.
Full text
1,880 characters
· extracted from
oa-doi-fallback
· click to expand
Abstract
Whether individual transcripts carry intrinsic features that predetermine their response to perturbations is unknown. Here we used nanopore direct RNA sequencing of male mouse dorsal root ganglia (DRG) to simultaneously profile N6-methyladenosine (m6A) modifications, poly(A) tail dynamics, and full-length isoform identity from mice treated with bortezomib, a proteasome inhibitor that causes painful peripheral neuropathy. Machine learning revealed that transcript-intrinsic features predetermine the magnitude of perturbation-induced m6A loss (R² = 0.983). Expression level contributed just 2.6% of predictive importance. Bortezomib removed a fixed ∼73.5% fraction of m6A marks, meaning absolute loss scaled linearly with baseline density and a transcript’s epitranscriptomic fate was encoded in its architecture before drug exposure. Unsupervised clustering identified four response programs where the dominant m6A erosion cluster enriched for oxidative phosphorylation (OXPHOS, p = 1.0 × 10⁻¹⁷) and proteasome (p = 2.8 × 10⁻¹⁰) genes, recapitulating bortezomib’s established mechanisms without prior biological knowledge. Isoform-resolved analysis uncovered m6A remodeling patterns suggesting post-transcriptional regulation of glycolytic and OXPHOS genes, and Western blot confirmed protein-level suppression of OXPHOS components. Integration with single-nuclei sequencing showed sensory neurons carried 2.2-fold greater m6A loss burden than non-neuronal cells, a direct consequence of architectural determinism applied to cell-type-specific transcriptomes. These findings establish that epitranscriptomic bortezomib response is predetermined by transcript architecture, with pathway specificity and cell-type vulnerability emerging as downstream consequences of intrinsic RNA structure.
Competing Interest Statement
The authors have declared no competing interest.
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.