Nanopore sequencing for N1-methylpseudouridine in RNA reveals sequence-dependent discrimination of the modified nucleotide triphosphate during transcription.

Direct RNA sequencing with a commercial nanopore platform was used to sequence RNA containing uridine (U), pseudouridine (Ψ) or N1-methylpseudouridine (m1Ψ) in >100 different 5-nucleotide contexts. The base calling data for Ψ or m1Ψ were similar but different from U allowing their detection. Understanding the nanopore signatures for Ψ and m1Ψ enabled a running start T7 RNA polymerase assay to study the selection of UTP versus ΨTP or m1ΨTP competing mixtures in all possible adjacent sequence contexts. A significant sequence context dependency was observed for T7 RNA polymerase with insertion yields for ΨTP versus UTP spanning a range of 20-65%, and m1ΨTP versus UTP producing variable yields that differ by 15-70%. Experiments with SP6 RNA polymerase, as well as chemically-modified triphosphates and DNA templates provide insight to explain the observations. The SP6 polymerase introduced m1ΨTP when competed with UTP with a smaller window of yields (15-30%) across all sequence contexts studied. These results may aid in future efforts that employ RNA polymerases to make therapeutic mRNAs with sub-stoichiometric amounts of m1Ψ.

DNA Damage Accelerates G-Quadruplex Folding in a Duplex-G-Quadruplex-Duplex Context.

Molecular details for the impact of DNA damage on folding of potential G-quadruplex sequences (PQSs) to noncanonical DNA structures involved in gene regulation are poorly understood. Here, the effects of DNA base damage and strand breaks on PQS folding kinetics were studied in the context of the VEGF promoter sequence embedded between two DNA duplex anchors, termed a duplex-G-quadruplex-duplex (DGD) motif. This DGD scaffold imposes constraints on the PQS folding process that more closely mimic those found in genomic DNA. Folding kinetics were monitored by circular dichroism (CD) to find folding half-lives ranging from 2 s to 12 min depending on the DNA damage type and sequence position. The presence of Mg2+ ions and G-quadruplex (G4)-binding protein APE1 facilitated the folding reactions. A strand break placing all four G runs required for G4 formation on one side of the break accelerated the folding rate by >150-fold compared to the undamaged sequence. Combined 1D 1H NMR and CD analyses confirmed that isothermal folding of the VEGF-DGD constructs yielded spectral signatures that suggest the formation of G4 motifs and demonstrated a folding dependency on the nature and location of DNA damage. Importantly, the PQS folding half-lives measured are relevant to replication, transcription, and DNA repair time frames.

Bisulfite and Nanopore Sequencing for Pseudouridine in RNA.

Nucleophilic addition of bisulfite to pyrimidine bases has been known for a half century, and the reaction has been in use for at least a quarter of a century for identifying 5-methylcytidine in DNA. This account focuses on the chemistry of bisulfite with pseudouridine, an isomer of the RNA nucleoside uridine in which the uracil base is connected to C1' of ribose via C5 instead of N1. Pseudouridine, Ψ, is the most common nucleotide modification found in cellular RNA overall, in part due to its abundance in rRNAs and tRNAs. It has a stabilizing influence on RNA structure because N1 is now available for additional hydrogen bonding and because the heterocycle is slightly better at π stacking. The isomerization of U to Ψ in RNA strands is catalyzed by 13 different enzymes in humans and 11 in E. coli; some of these enzymes are implicated in disease states which is testament to the biological importance of pseudouridine in cells. Recently, pseudouridine came into the limelight as the key modification that, after N1 methylation, enables mRNA vaccines to be delivered efficiently into human tissue with minimal generation of a deleterious immunogenic response. Here we describe the bisulfite reaction with pseudouridine which gives rise to a chemical sequencing method to map the modified base in the epitranscriptome. Unlike the reaction with cytidine, the addition of bisulfite to Ψ leads irreversibly to form an adduct that is bypassed during cDNA synthesis by reverse transcriptases yielding a characteristic deletion signature. Although there were hints to the structure of the bisulfite adduct(s) 30 to 50 years ago, it took modern spectroscopic and computational methods to solve the mystery. Raman spectroscopy along with extensive NMR, ECD, and computational work led to the assignment of the major product as the (R) diastereomer of an oxygen adduct at C1' of a ring-opened pseudouridine. Mechanistically, this arose from a succession of conjugate addition, E2 elimination, and a [2,3] sigmatropic rearrangement, all of which are stereodefined reactions. A minor reaction with excess bisulfite led to the (S) isomer of a S-adducted SO3- group. Understanding structure and mechanism aided the design of a Ψ-specific sequencing reaction and guided attempts to improve the utility and specificity of the method. Separately, we have been investigating the use of nanopore direct RNA sequencing, a single-molecule method that directly analyzes RNA strands isolated from cells after end-ligation of adaptor sequences. By combining the electrical current and base-calling data from the nanopore with dwell-time analysis from the helicase employed to deliver RNA to the nanopore, we were able to map Ψ sites in nearly all sequence contexts. This analysis was employed to find Ψ residues in the SARS-CoV-2 vRNA, to analyze the sequence context effects of mRNA vaccine synthesis via in vitro transcription, and to evaluate the impact of stress on chemical modifications in the E. coli ribosome. Most recently, we found that bisulfite treatment of RNA leading to Ψ adducts could modulate the nanopore signal to help in mapping modifications of low occupancy.

RNA polymerase II stalls on oxidative DNA damage via a torsion-latch mechanism involving lone pair-π and CH-π interactions.

Oxidation of guanine generates several types of DNA lesions, such as 8-oxoguanine (8OG), 5-guanidinohydantoin (Gh), and spiroiminodihydantoin (Sp). These guanine-derived oxidative DNA lesions interfere with both replication and transcription. However, the molecular mechanism of transcription processing of Gh and Sp remains unknown. In this study, by combining biochemical and structural analysis, we revealed distinct transcriptional processing of these chemically related oxidized lesions: 8OG allows both error-free and error-prone bypass, whereas Gh or Sp causes strong stalling and only allows slow error-prone incorporation of purines. Our structural studies provide snapshots of how polymerase II (Pol II) is stalled by a nonbulky Gh lesion in a stepwise manner, including the initial lesion encounter, ATP binding, ATP incorporation, jammed translocation, and arrested states. We show that while Gh can form hydrogen bonds with adenosine monophosphate (AMP) during incorporation, this base pair hydrogen bonding is not sufficient to hold an ATP substrate in the addition site and is not stable during Pol II translocation after the chemistry step. Intriguingly, we reveal a unique structural reconfiguration of the Gh lesion in which the hydantoin ring rotates ∼90° and is perpendicular to the upstream base pair planes. The perpendicular hydantoin ring of Gh is stabilized by noncanonical lone pair-π and CH-π interactions, as well as hydrogen bonds. As a result, the Gh lesion, as a functional mimic of a 1,2-intrastrand crosslink, occupies canonical -1 and +1 template positions and compromises the loading of the downstream template base. Furthermore, we suggest Gh and Sp lesions are potential targets of transcription-coupled repair.

Riboflavin Stabilizes Abasic, Oxidized G-Quadruplex Structures.

The G-quadruplex is a noncanonical fold of DNA commonly found at telomeres and within gene promoter regions of the genome. These guanine-rich sequences are highly susceptible to damages such as base oxidation and depurination, leading to abasic sites. In the present work, we address whether a vacancy, such as an abasic site, in a G-quadruplex serves as a specific ligand recognition site. When the G-tetrad is all guanines, the vacant (abasic) site is recognized and bound by free guanine nucleobase. However, we aim to understand whether the preference for a specific ligand recognition changes with the presence of a guanine oxidation product 8-oxo-7,8-dihydroguanine (OG) adjacent to the vacancy in the tetrad. Using molecular dynamics simulation, circular dichroism, and nuclear magnetic resonance, we examined the ability for riboflavin to stabilize abasic site-containing G-quadruplex structures. Through structural and free energy binding analysis, we observe riboflavin's ability to stabilize an abasic site-containing G-quadruplex only in the presence of an adjacent OG-modified base. Further, when compared to simulation with the vacancy filled by free guanine, we observe that the free guanine nucleobase is pushed outside of the tetrad by OG to interact with other parts of the structure, including loop residues. These results support the preference of riboflavin over free guanine to fill an OG-adjacent G-quadruplex abasic vacancy.

Second Harmonic Generation Interrogation of the Endonuclease APE1 Binding Interaction with G-Quadruplex DNA.

The binding interaction between the DNA repair enzyme apurinic/apyrimidinic endonuclease-1 (APE1) with promoter G-quadruplex (G4) folds bearing an abasic site (AP) can serve as a gene regulatory switch during oxidative stress. Prior fluorescence-based analysis in solution suggested APE1 binds the VEGF promoter G4 but whether this interaction was specific or not remained an open question. Second harmonic generation (SHG) was used in this work to measure the noncanonical DNA-protein binding interaction in a label-free assay with high sensitivity to demonstrate the interaction is ordered and specific. The binding of APE1 to the VEGF promoter G4 with AP sites modeled by a tetrahydrofuran analogue produced dissociation constants of ∼100 nM that differed from duplex and single-stranded DNA control studies. The SHG measurements confirmed APE1 binds the VEGF G4 folds in a specific manner resolving a remaining question regarding how this endonuclease with gene regulatory features engages G4 folds. The studies demonstrate the power of SHG to interrogate noncanonical DNA-protein interactions providing a foundational example for the use of this analytical method in future biochemical analyses.

Identification of the Major Product of Guanine Oxidation in DNA by Ozone.

Ozonolysis of guanosine formed the 5-carboxamido-5-formamido-2-iminohydantoin (2Ih) nucleoside along with trace spiroiminodihydantoin (Sp). On the basis of literature precedent, we propose an unconventional ozone mechanism involving incorporation of only one oxygen atom of O3 to form 2Ih with evolution of singlet oxygen responsible for Sp formation. The increased yield of Sp in the buffered 1O2-stabilizing solvent D2O, formation of 2Ih in a short oligodeoxynucleotide, and 18O-isotope labeling provided evidence to support this mechanism. The elusiveness and challenges of working with 2Ih are described in a series of studies on the significant context effects on the half-life of the 2Ih glycosidic bond.

Insights into the 5-Carboxamido-5-Formamido-2-Iminohydantoin Structural Isomerization Equilibria.

Exposure of DNA to oxidants results in modification of the electron-rich guanine heterocycle including formation of the mutagenic 5-carboxamido-5-formamido-2-iminohydantoin (2Ih) lesion. Previously thought to exist solely as a pair of diastereomers, we found under biologically relevant conditions that 2Ih reversibly closes to a formerly hypothetical intermediate and opens into a newly discovered regioisomer. In a nucleoside model, only ∼80% of 2Ih existed as the structure studied over the last 20 years with significant isomeric products persisting in buffered aqueous solution.

Transcriptome-wide profiling of multiple RNA modifications simultaneously at single-base resolution.

The breadth and importance of RNA modifications are growing rapidly as modified ribonucleotides can impact the sequence, structure, function, stability, and fate of RNAs and their interactions with other molecules. Therefore, knowing cellular RNA modifications at single-base resolution could provide important information regarding cell status and fate. A current major limitation is the lack of methods that allow the reproducible profiling of multiple modifications simultaneously, transcriptome-wide and at single-base resolution. Here we developed RBS-Seq, a modification of RNA bisulfite sequencing that enables the sensitive and simultaneous detection of m5C, Ψ, and m1A at single-base resolution transcriptome-wide. With RBS-Seq, m5C and m1A are accurately detected based on known signature base mismatches and are detected here simultaneously along with Ψ sites that show a 1-2 base deletion. Structural analyses revealed the mechanism underlying the deletion signature, which involves Ψ-monobisulfite adduction, heat-induced ribose ring opening, and Mg2+-assisted reorientation, causing base-skipping during cDNA synthesis. Detection of each of these modifications through a unique chemistry allows high-precision mapping of all three modifications within the same RNA molecule, enabling covariation studies. Application of RBS-Seq on HeLa RNA revealed almost all known m5C, m1A, and ψ sites in tRNAs and rRNAs and provided hundreds of new m5C and Ψ sites in noncoding RNAs and mRNAs. However, our results diverge greatly from earlier work, suggesting ∼10-fold fewer m5C sites in noncoding and coding RNAs and the absence of substantial m1A in mRNAs. Taken together, the approaches and refined datasets in this work will greatly enable future epitranscriptome studies.

Collateral Damage Occurs When Using Photosensitizer Probes to Detect or Modulate Nucleic Acid Modifications.

Nucleic acids are chemically modified to fine-tune their properties for biological function. Chemical tools for selective tagging of base modifications enables new approaches; the photosensitizers riboflavin and anthraquinone were previously proposed to oxidize N6 -methyladenine (m6 A) or 5-methylcytosine (5mdC) selectively. Herein, riboflavin, anthraquinone, or Rose Bengal were allowed to react with the canonical nucleosides dA, dC, dG, and dT, and the modified bases 5mdC, m6 A, 8-oxoguanine (dOG), and 8-oxoadenine (dOA) to rank their reactivities. The nucleoside studies reveal that dOG is the most reactive and that the native nucleoside dG is higher or similar in reactivity to 5mdC or m6 A; competition in both single- and double-stranded DNA of dG vs. 5mdC or 6mdA for oxidant confirmed that dG is favorably oxidized. Thus, photosensitizers are promiscuous nucleic acid oxidants with poor chemoselectivity that will negatively impact attempts at targeted oxidation of modified nucleotides in cells.

Hysteresis in poly-2'-deoxycytidine i-motif folding is impacted by the method of analysis as well as loop and stem lengths.

In DNA, i-motif (iM) folds occur under slightly acidic conditions when sequences rich in 2'-deoxycytidine (dC) nucleotides adopt consecutive dC self base pairs. The pH stability of an iM is defined by the midpoint in the pH transition (pHT ) between the folded and unfolded states. Two different experiments to determine pHT values via circular dichroism (CD) spectroscopy were performed on poly-dC iMs of length 15, 19, or 23 nucleotides. These experiments demonstrate two points: (1) pHT values were dependent on the titration experiment performed, and (2) pH-induced denaturing or annealing processes produced isothermal hysteresis in the pHT values. These results in tandem with model iMs with judicious mutations of dC to thymidine to favor particular folds found the hysteresis was maximal for the shorter poly-dC iMs and those with an even number of base pairs, while the hysteresis was minimal for longer poly-dC iMs and those with an odd number of base pairs. Experiments to follow the iM folding via thermal changes identified thermal hysteresis between the denaturing and annealing cycles. Similar trends were found to those observed in the CD experiments. The results demonstrate that the method of iM analysis can impact the pHT parameter measured, and hysteresis was observed in the pHT and Tm values.

Location dependence of the transcriptional response of a potential G-quadruplex in gene promoters under oxidative stress.

Oxidation of the guanine (G) heterocycle to 8-oxo-7,8-dihydroguanine (OG) in mammalian gene promoters was demonstrated to induce transcription. Potential G-quadruplex forming sequences (PQSs) in promoters have a high density of G nucleotides rendering them highly susceptible to oxidation and possible gene activation. The VEGF PQS with OG or an abasic site were synthesized at key locations in the SV40 or HSV-TK model promoters to determine the location dependency in the gene expression profile in human cells. The PQS location with respect to the transcription start site (TSS) and strand of occupancy (coding versus non-coding strand) are key parameters that determine the magnitude and direction in which gene expression changes with the chemically modified VEGF PQS. The greatest impact observed for OG or F in the PQS context in these promoters was within ∼200 bp of the TSS. Established PQSs found to occur naturally in a similar location relative to the TSS for possible oxidation-induced gene activation include c-MYC, KRAS, c-KIT, HIF-1α, PDGF-A and hTERT. The studies provide experimental constraints that were used to probe bioinformatic data regarding PQSs in the human genome for those that have the possibility to be redox switches for gene regulation.

On the irrelevancy of hydroxyl radical to DNA damage from oxidative stress and implications for epigenetics.

Contrary to frequent reports in the literature, hydroxyl radical is not a key species participating in endogenous oxidative DNA damage. Instead, carbonate radical anion is formed from the Fenton reaction under cellular conditions and from decomposition of nitrosoperoxycarbonate generated during inflammation. Carbonate radical anion is a potent one-electron oxidant capable of generating base radical cations that can migrate over long distances in duplex DNA, ultimately generating 8-oxo-7,8-dihydroguanine at a redox-sensitive sequence such as GGG. Such a mechanism enables G-quadruplex-forming sequences to act as long-range sensors of oxidative stress, impacting gene expression via the DNA repair mechanism that reads and ultimately erases the oxidized base. With a writing, reading and erasing mechanism in place, oxidative 'damage' to DNA might be relabeled as 'epigenetic' modifications.

Cruciform DNA Sequences in Gene Promoters Can Impact Transcription upon Oxidative Modification of 2'-Deoxyguanosine.

Sequences of DNA typically adopt B-form duplexes in genomes, although noncanonical structures such as G-quadruplexes, i-motifs, Z-DNA, and cruciform structures can occur. A challenge is to determine the functions of these various structures in cellular processes. We and others have hypothesized that G-rich G-quadruplex-forming sequences in human genome promoters serve to sense oxidative damage generated during oxidative stress impacting gene regulation. Herein, chemical tools and a cell-based assay were used to study the oxidation of guanine to 8-oxo-7,8-dihydroguanine (OG) in the context of a cruciform-forming sequence in a gene promoter to determine the impact on transcription. We found that OG in the nontemplate strand in the loop of a cruciform-forming sequence could induce gene expression; conversely when OG was in the same sequence on the template strand, gene expression was inhibited. A model for the transcriptional changes observed is proposed in which OG focuses the DNA repair process on the promoter to impact expression. Our cellular and biophysical studies and literature sources support the idea that removal of OG from duplex DNA by OGG1 yields an abasic site (AP) that triggers a structural shift to the cruciform fold. The AP-bearing cruciform structure is presented to APE1, which functions as a conduit between DNA repair and gene regulation. The significance is enhanced by a bioinformatic study of all human gene promoters and transcription termination sites for inverted repeats (IRs). Comparison of the two regions showed that promoters have stable and G-rich IRs at a low frequency and termination sites have many AT-rich IRs with low stability.

Interplay of Guanine Oxidation and G-Quadruplex Folding in Gene Promoters.

Living in an oxygen atmosphere demands an ability to thrive in the presence of reactive oxygen species (ROS). Aerobic organisms have successfully found solutions to the oxidative threats imposed by ROS by evolving an elaborate detoxification system, upregulating ROS during inflammation, and utilizing ROS as messenger molecules. In this Perspective, recent studies are discussed that demonstrate ROS as signaling molecules for gene regulation by combining two emergent properties of the guanine (G) heterocycle in DNA, namely, oxidation sensitivity and a propensity for G-quadruplex (G4) folding, both of which depend upon sequence context. In human gene promoters, this results from an elevated 5'-GG-3' dinucleotide frequency and GC enrichment near transcription start sites. Oxidation of DNA by ROS drives conversion of G to 8-oxo-7,8-dihydroguanine (OG) to mark target promoters for base excision repair initiated by OG-glycosylase I (OGG1). Sequence-dependent mechanisms for gene activation are available to OGG1 to induce transcription. Either OGG1 releases OG to yield an abasic site driving formation of a non-canonical fold, such as a G4, to be displayed to apurinic/apyrimidinic 1 (APE1) and stalling on the fold to recruit activating factors, or OGG1 binds OG and facilitates activator protein recruitment. The mechanisms described drive induction of stress response, DNA repair, or estrogen-induced genes, and these pathways are novel potential anticancer targets for therapeutic intervention. Chemical concepts provide a framework to discuss the regulatory or possible epigenetic potential of the OG modification in DNA, in which DNA "damage" and non-canonical folds collaborate to turn on or off gene expression. The next steps for scientific discovery in this growing field are discussed.

Oxidative DNA damage is epigenetic by regulating gene transcription via base excision repair.

Reactive oxygen species (ROS) have emerged as important cellular-signaling agents for cellular survival. Herein, we demonstrate that ROS-mediated oxidation of DNA to yield 8-oxo-7,8-dihydroguanine (OG) in gene promoters is a signaling agent for gene activation. Enhanced gene expression occurs when OG is formed in guanine-rich, potential G-quadruplex-forming sequences (PQS) in promoter-coding strands, initiating base excision repair (BER) by 8-oxoguanine DNA glycosylase (OGG1), yielding an abasic site (AP). The AP enables melting of the duplex to unmask the PQS, adopting a G-quadruplex fold in which apurinic/apyrimidinic endonuclease 1 (APE1) binds, but inefficiently cleaves, the AP for activation of vascular endothelial growth factor (VEGF) or endonuclease III-like protein 1 (NTHL1) genes. These details were mapped via synthesis of OG and AP analogs at single-nucleotide precision within the promoter of a luciferase reporter system. The reporters were analyzed in human and mouse cells while selectively knocking out or down critical BER proteins to identify the impact on luciferase expression. Identification of the oxidatively modified DNA base OG to guide BER activity in a gene promoter and impact cellular phenotype ascribes an epigenetic role to OG.

Human NEIL3 Gene Expression Regulated by Epigenetic-Like Oxidative DNA Modification.

The NEIL3 DNA repair gene is induced in cells or animal models experiencing oxidative or inflammatory stress along with oxidation of guanine (G) to 8-oxo-7,8-dihydroguanine (OG) in the genome. We hypothesize that a G-rich promoter element that is a potential G-quadruplex-forming sequence (PQS) in NEIL3 is a site for introduction of OG with epigenetic-like potential for gene regulation. Activation occurs when OG is formed in the NEIL3 PQS located near the transcription start site. Oxidative stress either introduced by TNFα or synthetically incorporated into precise locations focuses the base excision repair process to read and catalyze removal of OG via OG-glycosylase I (OGG1), yielding an abasic site (AP). Thermodynamic studies showed that AP destabilizes the duplex, enabling a structural transition of the sequence to a G-quadruplex (G4) fold that positions the AP in a loop facilitated by the NEIL3 PQS having five G runs in which the four unmodified runs adopt a stable G4. This presents AP to apurinic/apyrimidinic endonuclease 1 (APE1) that poorly cleaves the AP backbone in this context according to in vitro studies, allowing the protein to function as a trans activator of transcription. The proposal is supported by chemical studies in cellulo and in vitro. Activation of NEIL3 expression via the proposed mechanism allows cells to respond to mutagenic DNA damage removed by NEIL3 associated with oxidative or inflammatory stress. Lastly, inspection of many mammalian genomes identified conservation of the NEIL3 PQS, suggesting this sequence was favorably selected to function as a redox switch with OG as the epigenetic-like regulatory modification.

Effect of Oxidative Damage on Charge and Spin Transport in DNA.

A Hall device was used for measuring spin polarization on electrons that are either reorganized within the molecules or transmitted through the self-assembled monolayers of DNA adsorbed on the device surface. We were able to observe spin-dependent charge polarization and charge transport through double-stranded DNA of various lengths and through double-stranded DNA containing oxidative damage. We found enhancement in the spin-dependent transport through oxidatively damaged DNA. This phenomenon can be rationalized either by assuming that the damaged DNA is characterized by a higher barrier for conduction or by charge transfer through the DNA being conducted through at least two channels, one involves the bases and is highly conductive but less spin selective, while the other pathway is mainly through the ribophosphate backbone and it is the minor one in terms of charge transmission efficiency, but it is highly spin selective.

Structural Elucidation of Bisulfite Adducts to Pseudouridine That Result in Deletion Signatures during Reverse Transcription of RNA.

The recent report of RBS-Seq to map simultaneously the epitranscriptomic modifications N1-methyladenosine, 5-methylcytosine, and pseudouridine (Ψ) via bisulfite treatment of RNA provides a key advance to locate these important modifications. The locations of Ψ were found by a deletion signature generated during cDNA synthesis after bisulfite treatment for which the chemical details of the reaction are poorly understood. In the present work, the bisulfite reaction with Ψ was explored to identify six isomers of bisulfite adducted to Ψ. We found four of these adducts involved the heterocyclic ring, similar to the reaction with other pyrimidines. The remaining two adducts were bonded to the 1' carbon, which resulted in opening of the ribose ring. The utilization of complementary 1D- and 2D-NMR, Raman, and electronic circular dichroism spectroscopies led to the assignment of the two ribose adducts being the constitutional isomers of an S- and an O-adduct of bisulfite to the ribose, and these are the final products after heating. A mechanistic proposal is provided to rationalize chemically the formation and stereochemistries of all six isomeric bisulfite adducts to Ψ; conversion of intermediate adducts to the two final products is proposed to involve E2, SN2', and [2,3]-sigmatropic shift reactions. Lastly, a synthetic RNA template with Ψ at a known location was treated with bisulfite, leading to a deletion signature after reverse transcription, supporting the RBS-Seq report. This classical bisulfite reaction used for epigenomic and epitranscriptomic sequencing diverges from the C nucleoside Ψ to form stable bisulfite end products that yield signatures for next-generation sequencing.