Discovery of cytosine deaminases enables base-resolution methylome mapping using a single enzyme.

Deaminases have important uses in modification detection and genome editing. However, the range of applications is limited by the small number of characterized enzymes. To expand the toolkit of deaminases, we developed an in vitro approach that bypasses a major hurdle with their toxicity in cells. We assayed 175 putative cytosine deaminases on a variety of substrates and found a broad range of activity on double- and single-stranded DNA in various sequence contexts, including CpG-specific deaminases and enzymes without sequence preference. We also characterized enzyme selectivity across six DNA modifications and reported enzymes that do not deaminate modified cytosines. The detailed analysis of diverse deaminases opens new avenues for biotechnological and medical applications. As a demonstration, we developed SEM-seq, a non-destructive single-enzyme methylation sequencing method using a modification-sensitive double-stranded DNA deaminase. The streamlined protocol enables accurate, base-resolution methylome mapping of scarce biological material, including cell-free DNA and 10 pg input DNA.

Nondestructive enzymatic deamination enables single-molecule long-read amplicon sequencing for the determination of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution.

The predominant methodology for DNA methylation analysis relies on the chemical deamination by sodium bisulfite of unmodified cytosine to uracil to permit the differential readout of methylated cytosines. Bisulfite treatment damages the DNA, leading to fragmentation and loss of long-range methylation information. To overcome this limitation of bisulfite-treated DNA, we applied a new enzymatic deamination approach, termed enzymatic methyl-seq (EM-seq), to long-range sequencing technologies. Our methodology, named long-read enzymatic modification sequencing (LR-EM-seq), preserves the integrity of DNA, allowing long-range methylation profiling of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) over multikilobase length of genomic DNA. When applied to known differentially methylated regions (DMRs), LR-EM-seq achieves phasing of >5 kb, resulting in broader and better defined DMRs compared with that previously reported. This result showed the importance of phasing methylation for biologically relevant questions and the applicability of LR-EM-seq for long-range epigenetic analysis at single-molecule and single-nucleotide resolution.

Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA.

Bisulfite sequencing detects 5mC and 5hmC at single-base resolution. However, bisulfite treatment damages DNA, which results in fragmentation, DNA loss, and biased sequencing data. To overcome these problems, enzymatic methyl-seq (EM-seq) was developed. This method detects 5mC and 5hmC using two sets of enzymatic reactions. In the first reaction, TET2 and T4-BGT convert 5mC and 5hmC into products that cannot be deaminated by APOBEC3A. In the second reaction, APOBEC3A deaminates unmodified cytosines by converting them to uracils. Therefore, these three enzymes enable the identification of 5mC and 5hmC. EM-seq libraries were compared with bisulfite-converted DNA, and each library type was ligated to Illumina adaptors before conversion. Libraries were made using NA12878 genomic DNA, cell-free DNA, and FFPE DNA over a range of DNA inputs. The 5mC and 5hmC detected in EM-seq libraries were similar to those of bisulfite libraries. However, libraries made using EM-seq outperformed bisulfite-converted libraries in all specific measures examined (coverage, duplication, sensitivity, etc.). EM-seq libraries displayed even GC distribution, better correlations across DNA inputs, increased numbers of CpGs within genomic features, and accuracy of cytosine methylation calls. EM-seq was effective using as little as 100 pg of DNA, and these libraries maintained the described advantages over bisulfite sequencing. EM-seq library construction, using challenging samples and lower DNA inputs, opens new avenues for research and clinical applications.

RNase H-based analysis of synthetic mRNA 5' cap incorporation.

Advances in mRNA synthesis and lipid nanoparticles technologies have helped make mRNA therapeutics and vaccines a reality. The 5' cap structure is a crucial modification required to functionalize synthetic mRNA for efficient protein translation in vivo and evasion of cellular innate immune responses. The extent of 5' cap incorporation is one of the critical quality attributes in mRNA manufacturing. RNA cap analysis involves multiple steps: generation of predefined short fragments from the 5' end of the kilobase-long synthetic mRNA molecules using RNase H, a ribozyme or a DNAzyme, enrichment of the 5' cleavage products, and LC-MS intact mass analysis. In this paper, we describe (1) a framework to design site-specific RNA cleavage using RNase H; (2) a method to fluorescently label the RNase H cleavage fragments for more accessible readout methods such as gel electrophoresis or high-throughput capillary electrophoresis; (3) a simplified method for post-RNase H purification using desthiobiotinylated oligonucleotides and streptavidin magnetic beads followed by elution using water. By providing a design framework for RNase H-based RNA 5' cap analysis using less resource-intensive analytical methods, we hope to make RNA cap analysis more accessible to the scientific community.

Identification of high-confidence human poly(A) RNA isoform scaffolds using nanopore sequencing.

Nanopore sequencing devices read individual RNA strands directly. This facilitates identification of exon linkages and nucleotide modifications; however, using conventional direct RNA nanopore sequencing, the 5' and 3' ends of poly(A) RNA cannot be identified unambiguously. This is due in part to RNA degradation in vivo and in vitro that can obscure transcription start and end sites. In this study, we aimed to identify individual full-length human RNA isoforms among ∼4 million nanopore poly(A)-selected RNA reads. First, to identify RNA strands bearing 5' m7G caps, we exchanged the biological cap for a modified cap attached to a 45-nt oligomer. This oligomer adaptation method improved 5' end sequencing and ensured correct identification of the 5' m7G capped ends. Second, among these 5'-capped nanopore reads, we screened for features consistent with a 3' polyadenylation site. Combining these two steps, we identified 294,107 individual high-confidence full-length RNA scaffolds from human GM12878 cells, most of which (257,721) aligned to protein-coding genes. Of these, 4876 scaffolds indicated unannotated isoforms that were often internal to longer, previously identified RNA isoforms. Orthogonal data for m7G caps and open chromatin, such as CAGE and DNase-HS seq, confirmed the validity of these high-confidence RNA scaffolds.

Human RNase 4 improves mRNA sequence characterization by LC-MS/MS.

With the rapid growth of synthetic messenger RNA (mRNA)-based therapeutics and vaccines, the development of analytical tools for characterization of long, complex RNAs has become essential. Tandem liquid chromatography-mass spectrometry (LC-MS/MS) permits direct assessment of the mRNA primary sequence and modifications thereof without conversion to cDNA or amplification. It relies upon digestion of mRNA with site-specific endoribonucleases to generate pools of short oligonucleotides that are then amenable to MS-based sequence analysis. Here, we showed that the uridine-specific human endoribonuclease hRNase 4 improves mRNA sequence coverage, in comparison with the benchmark enzyme RNase T1, by producing a larger population of uniquely mappable cleavage products. We deployed hRNase 4 to characterize mRNAs fully substituted with 1-methylpseudouridine (m1Ψ) or 5-methoxyuridine (mo5U), as well as mRNAs selectively depleted of uridine-two key strategies to reduce synthetic mRNA immunogenicity. Lastly, we demonstrated that hRNase 4 enables direct assessment of the 5' cap incorporation into in vitro transcribed mRNA. Collectively, this study highlights the power of hRNase 4 to interrogate mRNA sequence, identity, and modifications by LC-MS/MS.

Chemical capping improves template switching and enhances sequencing of small RNAs.

Template-switching reverse transcription is widely used in RNA sequencing for low-input and low-quality samples, including RNA from single cells or formalin-fixed paraffin-embedded (FFPE) tissues. Previously, we identified the native eukaryotic mRNA 5' cap as a key structural element for enhancing template switching efficiency. Here, we introduce CapTS-seq, a new strategy for sequencing small RNAs that combines chemical capping and template switching. We probed a variety of non-native synthetic cap structures and found that an unmethylated guanosine triphosphate cap led to the lowest bias and highest efficiency for template switching. Through cross-examination of different nucleotides at the cap position, our data provided unequivocal evidence that the 5' cap acts as a template for the first nucleotide in reverse transcriptase-mediated post-templated addition to the emerging cDNA-a key feature to propel template switching. We deployed CapTS-seq for sequencing synthetic miRNAs, human total brain and liver FFPE RNA, and demonstrated that it consistently improves library quality for miRNAs in comparison with a gold standard template switching-based small RNA-seq kit.

Nanopore ReCappable sequencing maps SARS-CoV-2 5' capping sites and provides new insights into the structure of sgRNAs.

The SARS-CoV-2 virus has a complex transcriptome characterised by multiple, nested subgenomic RNAsused to express structural and accessory proteins. Long-read sequencing technologies such as nanopore direct RNA sequencing can recover full-length transcripts, greatly simplifying the assembly of structurally complex RNAs. However, these techniques do not detect the 5' cap, thus preventing reliable identification and quantification of full-length, coding transcript models. Here we used Nanopore ReCappable Sequencing (NRCeq), a new technique that can identify capped full-length RNAs, to assemble a complete annotation of SARS-CoV-2 sgRNAs and annotate the location of capping sites across the viral genome. We obtained robust estimates of sgRNA expression across cell lines and viral isolates and identified novel canonical and non-canonical sgRNAs, including one that uses a previously un-annotated leader-to-body junction site. The data generated in this work constitute a useful resource for the scientific community and provide important insights into the mechanisms that regulate the transcription of SARS-CoV-2 sgRNAs.

Pathways of thymidine hypermodification.

The DNAs of bacterial viruses are known to contain diverse, chemically complex modifications to thymidine that protect them from the endonuclease-based defenses of their cellular hosts, but whose biosynthetic origins are enigmatic. Up to half of thymidines in the Pseudomonas phage M6, the Salmonella phage ViI, and others, contain exotic chemical moieties synthesized through the post-replicative modification of 5-hydroxymethyluridine (5-hmdU). We have determined that these thymidine hypermodifications are derived from free amino acids enzymatically installed on 5-hmdU. These appended amino acids are further sculpted by various enzyme classes such as radical SAM isomerases, PLP-dependent decarboxylases, flavin-dependent lyases and acetyltransferases. The combinatorial permutations of thymidine hypermodification genes found in viral metagenomes from geographically widespread sources suggests an untapped reservoir of chemical diversity in DNA hypermodifications.

Phage-encoded ten-eleven translocation dioxygenase (TET) is active in C5-cytosine hypermodification in DNA.

TET/JBP (ten-eleven translocation/base J binding protein) enzymes are iron(II)- and 2-oxo-glutarate-dependent dioxygenases that are found in all kingdoms of life and oxidize 5-methylpyrimidines on the polynucleotide level. Despite their prevalence, few examples have been biochemically characterized. Among those studied are the metazoan TET enzymes that oxidize 5-methylcytosine in DNA to hydroxy, formyl, and carboxy forms and the euglenozoa JBP dioxygenases that oxidize thymine in the first step of base J biosynthesis. Both enzymes have roles in epigenetic regulation. It has been hypothesized that all TET/JBPs have their ancestral origins in bacteriophages, but only eukaryotic orthologs have been described. Here we demonstrate the 5mC-dioxygenase activity of several phage TETs encoded within viral metagenomes. The clustering of these TETs in a phylogenetic tree correlates with the sequence specificity of their genomically cooccurring cytosine C5-methyltransferases, which install the methyl groups upon which TETs operate. The phage TETs favor Gp5mC dinucleotides over the 5mCpG sites targeted by the eukaryotic TETs and are found within gene clusters specifying complex cytosine modifications that may be important for DNA packaging and evasion of host restriction.

Enzymatic characterization of mRNA cap adenosine-N6 methyltransferase PCIF1 activity on uncapped RNAs.

The phosphorylated RNA polymerase II CTD interacting factor 1 (PCIF1) is a methyltransferase that adds a methyl group to the N6-position of 2'O-methyladenosine (Am), generating N6, 2'O-dimethyladenosine (m6Am) when Am is the cap-proximal nucleotide. In addition, PCIF1 has ancillary methylation activities on internal adenosines (both A and Am), although with much lower catalytic efficiency relative to that of its preferred cap substrate. The PCIF1 preference for 2'O-methylated Am over unmodified A nucleosides is due mainly to increased binding affinity for Am. Importantly, it was recently reported that PCIF1 can methylate viral RNA. Although some viral RNA can be translated in the absence of a cap, it is unclear what roles PCIF1 modifications may play in the functionality of viral RNAs. Here we show, using in vitro assays of binding and methyltransfer, that PCIF1 binds an uncapped 5'-Am oligonucleotide with approximately the same affinity as that of a cap analog (KM = 0.4 versus 0.3 µM). In addition, PCIF1 methylates the uncapped 5'-Am with activity decreased by only fivefold to sixfold compared with its preferred capped substrate. We finally discuss the relationship between PCIF1-catalyzed RNA methylation, shown here to have broader substrate specificity than previously appreciated, and that of the RNA demethylase fat mass and obesity-associated protein (FTO), which demonstrates PCIF1-opposing activities on capped RNAs.

A sensitive approach to map genome-wide 5-hydroxymethylcytosine and 5-formylcytosine at single-base resolution.

Mapping genome-wide 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC) at single-base resolution is important to understand their biological functions. We present a cost-efficient mapping method that combines 5hmC-specific restriction enzyme PvuRts1I with a 5hmC chemical labeling enrichment method. The sensitive method enables detection of low-abundance 5hmC sites, providing a more complete 5hmC landscape than available bisulfite-based methods. This method generated a genome-wide 5fC map at single-base resolution. Parallel analyses revealed that 5hmC and 5fC in non-CpG context exhibit lower abundance, more dynamically, than those in CpG context. In the genic region, distribution of 5hmCpG and 5fCpG differed from 5hmCH and 5fCH (H = A, T, C). 5hmC and 5fC were distributed distinctly at regulatory protein-DNA binding sites, depleted in permissive transcription factor binding sites, and enriched at active and poised enhancers. This sensitive bisulfite conversion-free method can be applied to biological samples with limited starting material or low-abundance cytosine modifications.

E. coli RNase I exhibits a strong Ca2+-dependent inherent double-stranded RNase activity.

Since its initial characterization, Escherichia coli RNase I has been described as a single-strand specific RNA endonuclease that cleaves its substrate in a largely sequence independent manner. Here, we describe a strong calcium (Ca2+)-dependent activity of RNase I on double-stranded RNA (dsRNA), and a Ca2+-dependent novel hybridase activity, digesting the RNA strand in a DNA:RNA hybrid. Surprisingly, Ca2+ does not affect the activity of RNase I on single stranded RNA (ssRNA), suggesting a specific role for Ca2+ in the modulation of RNase I activity. Mutation of a previously overlooked Ca2+ binding site on RNase I resulted in a gain-of-function enzyme that is highly active on dsRNA and could no longer be stimulated by the metal. In summary, our data imply that native RNase I contains a bound Ca2+, allowing it to target both single- and double-stranded RNAs, thus having a broader substrate specificity than originally proposed for this traditional enzyme. In addition, the finding that the dsRNase activity, and not the ssRNase activity, is associated with the Ca2+-dependency of RNase I may be useful as a tool in applied molecular biology.

Genetic and biased agonist-mediated reductions in ß-arrestin recruitment prolong cAMP signaling at glucagon family receptors.

Receptors for the peptide hormones glucagon-like peptide-1 (GLP-1R), glucose-dependent insulinotropic polypeptide (GIPR), and glucagon (GCGR) are important regulators of insulin secretion and energy metabolism. GLP-1R agonists have been successfully deployed for the treatment of type 2 diabetes, but it has been suggested that their efficacy is limited by target receptor desensitization and downregulation due to recruitment of ß-arrestins. Indeed, recently described GLP-1R agonists with reduced ß-arrestin-2 recruitment have delivered promising results in preclinical and clinical studies. We therefore aimed to determine if the same phenomenon could apply to the closely related GIPR and GCGR. In HEK293 cells depleted of both ß-arrestin isoforms the duration of G protein-dependent cAMP/PKA signaling was increased in response to the endogenous ligand for each receptor. Moreover, in wildtype cells, "biased" GLP-1, GCG, and GIP analogs with selective reductions in ß-arrestin-2 recruitment led to reduced receptor endocytosis and increased insulin secretion over a prolonged stimulation period, although the latter effect was only seen at high agonist concentrations. Biased GCG analogs increased the duration of cAMP signaling, but this did not lead to increased glucose output from hepatocytes. Our study provides a rationale for the development of GLP-1R, GIPR, and GCGR agonists with reduced ß-arrestin recruitment, but further work is needed to maximally exploit this strategy for therapeutic purposes.

Exploiting Endogenous Enzymes for Cancer-Cell Selective Metabolic Labeling of RNA in Vivo.

Tissues and organs are composed of many diverse cell types, making cell-specific gene expression profiling a major challenge. Herein we report that endogenous enzymes, unique to a cell of interest, can be utilized to enable cell-specific metabolic labeling of RNA. We demonstrate that appropriately designed "caged" nucleosides can be rendered active by serving as a substrate for cancer-cell specific enzymes to enable RNA metabolic labeling, only in cancer cells. We envision that the ease and high stringency of our approach will enable expression analysis of tumor cells in complex environments.

Agonist-induced membrane nanodomain clustering drives GLP-1 receptor responses in pancreatic beta cells.

The glucagon-like peptide-1 receptor (GLP-1R), a key pharmacological target in type 2 diabetes (T2D) and obesity, undergoes rapid endocytosis after stimulation by endogenous and therapeutic agonists. We have previously highlighted the relevance of this process in fine-tuning GLP-1R responses in pancreatic beta cells to control insulin secretion. In the present study, we demonstrate an important role for the translocation of active GLP-1Rs into liquid-ordered plasma membrane nanodomains, which act as hotspots for optimal coordination of intracellular signaling and clathrin-mediated endocytosis. This process is dynamically regulated by agonist binding through palmitoylation of the GLP-1R at its carboxyl-terminal tail. Biased GLP-1R agonists and small molecule allosteric modulation both influence GLP-1R palmitoylation, clustering, nanodomain signaling, and internalization. Downstream effects on insulin secretion from pancreatic beta cells indicate that these processes are relevant to GLP-1R physiological actions and might be therapeutically targetable.

Enhanced expression and purification of nucleotide-specific ribonucleases MC1 and Cusativin.

Combinations of ribonucleases (RNases) are commonly used to digest RNA into oligoribonucleotide fragments prior to liquid chromatography-mass spectrometry (LC-MS) analysis. The distribution of the RNase target sequences or nucleobase sites within an RNA molecule is critical for achieving a high mapping coverage. Cusativin and MC1 are nucleotide-specific endoribonucleases encoded in the cucumber and bitter melon genomes, respectively. Their high specificity for cytidine (Cusativin) and uridine (MC1) make them ideal molecular biology tools for RNA modification mapping. However, heterogenous recombinant expression of either enzyme has been challenging because of their high toxicity to expression hosts and the requirement of posttranslational modifications. Here, we present two highly efficient and time-saving protocols that overcome these hurdles and enhance the expression and purification of these RNases. We first purified MC1 and Cusativin from bacteria by expressing and shuttling both enzymes to the periplasm as MBP-fusion proteins in T7 Express lysY/IqE. coli strain at low temperature. The RNases were enriched using amylose affinity chromatography, followed by a subsequent purification via a C-terminal 6xHIS tag. This fast, two-step purification allows for the purification of highly active recombinant RNases significantly surpassing yields reported in previous studies. In addition, we expressed and purified a Cusativin-CBD fusion enzyme in P. pastoris using chitin magnetic beads. Both Cusativin variants exhibited a similar sequence preference, suggesting that neither posttranslational modifications nor the epitope-tags have a substantial effect on the sequence specificity of the enzyme.

Identification and biosynthesis of thymidine hypermodifications in the genomic DNA of widespread bacterial viruses.

Certain viruses of bacteria (bacteriophages) enzymatically hypermodify their DNA to protect their genetic material from host restriction endonuclease-mediated cleavage. Historically, it has been known that virion DNAs from the Delftia phage ΦW-14 and the Bacillus phage SP10 contain the hypermodified pyrimidines α-putrescinylthymidine and α-glutamylthymidine, respectively. These bases derive from the modification of 5-hydroxymethyl-2'-deoxyuridine (5-hmdU) in newly replicated phage DNA via a pyrophosphorylated intermediate. Like ΦW-14 and SP10, the Pseudomonas phage M6 and the Salmonella phage ViI encode kinase homologs predicted to phosphorylate 5-hmdU DNA but have uncharacterized nucleotide content [Iyer et al. (2013) Nucleic Acids Res 41:7635-7655]. We report here the discovery and characterization of two bases, 5-(2-aminoethoxy)methyluridine (5-NeOmdU) and 5-(2-aminoethyl)uridine (5-NedU), in the virion DNA of ViI and M6 phages, respectively. Furthermore, we show that recombinant expression of five gene products encoded by phage ViI is sufficient to reconstitute the formation of 5-NeOmdU in vitro. These findings point to an unexplored diversity of DNA modifications and the underlying biochemistry of their formation.

Non-templated addition and template switching by Moloney murine leukemia virus (MMLV)-based reverse transcriptases co-occur and compete with each other.

Single-cell RNA-Seq (scRNA-Seq) has led to an unprecedented understanding of gene expression and regulation in individual cells. Many scRNA-Seq approaches rely upon the template switching property of Moloney murine leukemia virus (MMLV)-type reverse transcriptases. Template switching is believed to happen in a sequential process involving nontemplated addition of three protruding nucleotides (+CCC) to the 3'-end of the nascent cDNA, which can then anneal to the matching rGrGrG 3'-end of the template-switching oligo (TSO), allowing the reverse transcriptase (RT) to switch templates and continue copying the TSO sequence. In this study, we present a detailed analysis of template switching biases with respect to the RNA template, specifically of the role of the sequence and nature of its 5'-end (capped versus noncapped) in these biases. Our findings confirmed that the presence of a 5'-m7G cap enhances template switching efficiency. We also profiled the composition of the nontemplated addition in the absence of TSO and observed that the 5'-end of RNA template influences the terminal transferase activity of the RT. Furthermore, we found that designing new TSOs that pair with the most common nontemplated additions did little to improve template switching efficiency. Our results provide evidence suggesting that, in contrast to the current understanding of the template switching process, nontemplated addition and template switching are concurrent and competing processes.

Light-activated chemical probing of nucleobase solvent accessibility inside cells.

The discovery of functional RNAs that are critical for normal and disease physiology continues to expand at a breakneck pace. Many RNA functions are controlled by the formation of specific structures, and an understanding of each structural component is necessary to elucidate its function. Measuring solvent accessibility intracellularly with experimental ease is an unmet need in the field. Here, we present a novel method for probing nucleobase solvent accessibility, Light Activated Structural Examination of RNA (LASER). LASER depends on light activation of a small molecule, nicotinoyl azide (NAz), to measure solvent accessibility of purine nucleobases. In vitro, this technique accurately monitors solvent accessibility and identifies rapid structural changes resulting from ligand binding in a metabolite-responsive RNA. LASER probing can further identify cellular RNA-protein interactions and unique intracellular RNA structures. Our photoactivation technique provides an adaptable framework to structurally characterize solvent accessibility of RNA in many environments.