Dihydrouridine (D) is one of the most abundant tRNA modifications yet one of the most difficult to map quantitatively. Prior methods (D-seq, Rho-seq, AlkAniline-Seq) relied on RT truncation or chemical cleavage signals that cannot provide accurate stoichiometry and fail in densely modified D-loop regions. The authors developed CRACI-Seq to overcome these limitations, targeting quantitative single-base D mapping across the entire transcriptome.

Figure 1. CRACI-Seq method principle and validation.
KBH₄ reduction of dihydrouridine (D) to tetrahydrouridine enables HIV RT with high dGTP/dNTP ratio to incorporate C opposite the reduced site, generating a T→C misincorporation signature. Validation in synthetic RNA controls demonstrates ~96% conversion efficiency and single-base quantitative resolution. Adapted from Ju et al. 2025 (CC BY 4.0).
Methods: The study used total RNA from HEK293T cells (human), mouse tissues, and Arabidopsis seedlings. CRACI-Seq was performed with KBH₄ reduction followed by HIV RT under high dGTP conditions, Illumina sequencing, and custom bioinformatic analysis for T→C misincorporation calling. siRNA KD of each DUS writer (DUS1L–DUS4L) was used for enzyme assignment. LC-MS/MS provided orthogonal validation of D site stoichiometry.

Figure 2. DUS writer enzyme assignment by CRACI-Seq.
siRNA knockdown of individual DUS enzymes (DUS1L–DUS4L) followed by CRACI-Seq reveals writer-specific D site assignments. DUS1L KD → reduced D16/D17, DUS2L KD → reduced D20 (including mitochondrial D sites), DUS3L KD → reduced D47, DUS4L KD → reduced D20a. Heatmap shows site-specific stoichiometry changes upon each DUS KD. Adapted from Ju et al. 2025 (CC BY 4.0).
Results: CRACI-Seq identified D at positions D16, D17, D20, D20a, D20b, and D47 across cytoplasmic tRNAs with quantitative stoichiometry (D16/D17 >70%; D20a ranging 20–100%). Novel D sites were discovered in mitochondrial tRNAs (mt-tRNAAsn, mt-tRNAGln, mt-tRNALeu(UUR)) at moderate stoichiometries (<40%). DUS writer assignment revealed DUS1L→D16/17, DUS2L→D20, DUS3L→D47, DUS4L→D20a. Notably, D20a was found to cis-regulate D20 installation on the same tRNA. Only 8 high-confidence mRNA D sites were identified (10–40% stoichiometry), confirming that D is predominantly a tRNA modification.
D site detection and stoichiometry quantification — T→C misincorporation signals at D16, D17, D20, D47 in human tRNAs with background-corrected stoichiometry values.
tRNA D-loop heatmap — per-isodecoder D stoichiometry at D16, D17, D20, D20a, D20b, and D47 positions, with hierarchical clustering and isotype family annotation.
Differential D analysis — volcano plot and sample-level comparison showing D20 stoichiometry reduction upon DUS2L KD, with DUS writer assignment matrix.