SPAAC-NAD-seq, a sensitive and accurate method to profile NAD+-capped transcripts

Abstract

Significance The m7G cap is the canonical RNA cap in eukaryotes, but other noncanonical RNA caps exist, including the NAD+ cap. NAD captureSeq has been widely used to profile NAD+-capped RNAs (NAD-RNAs) in prokaryotes and eukaryotes. However, NAD captureSeq reacts at a low level with m7G-RNAs and introduces copper ions that cause RNA fragmentation, resulting in reduced sensitivity and loss of full-length sequence information. To address these issues, we developed the copper-free SPAAC-NAD-seq, which utilizes the strain-promoted azide–alkyne cycloaddition reaction to capture NAD-RNAs followed by high-throughput sequencing. Compared to NAD captureSeq, SPAAC-NAD-seq is more sensitive and retains full-length sequence information. This optimized technique provides a useful tool to profile NAD-RNAs in prokaryotes and, when combined with m7G-RNA depletion, in eukaryotes. Nicotinamide adenine diphosphate (NAD+) is a novel messenger RNA 5′ cap in Escherichia coli, yeast, mammals, and Arabidopsis. Transcriptome-wide identification of NAD+-capped RNAs (NAD-RNAs) was accomplished through NAD captureSeq, which combines chemoenzymatic RNA enrichment with high-throughput sequencing. NAD-RNAs are enzymatically converted to alkyne-RNAs that are then biotinylated using a copper-catalyzed azide–alkyne cycloaddition (CuAAC) reaction. Originally applied to E. coli RNA, which lacks the m7G cap, NAD captureSeq was then applied to eukaryotes without extensive verification of its specificity for NAD-RNAs vs. m7G-capped RNAs (m7G-RNAs). In addition, the Cu2+ ion in the CuAAC reaction causes RNA fragmentation, leading to greatly reduced yield and loss of full-length sequence information. We developed an NAD-RNA capture scheme utilizing the copper-free, strain-promoted azide–alkyne cycloaddition reaction (SPAAC). We examined the specificity of CuAAC and SPAAC reactions toward NAD-RNAs and m7G-RNAs and found that both prefer the former, but also act on the latter. We demonstrated that SPAAC-NAD sequencing (SPAAC-NAD-seq), when combined with immunodepletion of m7G-RNAs, enables NAD-RNA identification with accuracy and sensitivity, leading to the discovery of new NAD-RNA profiles in Arabidopsis. Furthermore, SPAAC-NAD-seq retained full-length sequence information. Therefore, SPAAC-NAD-seq would enable specific and efficient discovery of NAD-RNAs in prokaryotes and, when combined with m7G-RNA depletion, in eukaryotes.