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.

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.

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.

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.

A novel enrichment strategy reveals unprecedented number of novel transcription start sites at single base resolution in a model prokaryote and the gut microbiome.

BACKGROUND: The initiating nucleotide found at the 5' end of primary transcripts has a distinctive triphosphorylated end that distinguishes these transcripts from all other RNA species. Recognizing this distinction is key to deconvoluting the primary transcriptome from the plethora of processed transcripts that confound analysis of the transcriptome. The currently available methods do not use targeted enrichment for the 5'end of primary transcripts, but rather attempt to deplete non-targeted RNA. RESULTS: We developed a method, Cappable-seq, for directly enriching for the 5' end of primary transcripts and enabling determination of transcription start sites at single base resolution. This is achieved by enzymatically modifying the 5' triphosphorylated end of RNA with a selectable tag. We first applied Cappable-seq to E. coli, achieving up to 50 fold enrichment of primary transcripts and identifying an unprecedented 16539 transcription start sites (TSS) genome-wide at single base resolution. We also applied Cappable-seq to a mouse cecum sample and identified TSS in a microbiome. CONCLUSIONS: Cappable-seq allows for the first time the capture of the 5' end of primary transcripts. This enables a unique robust TSS determination in bacteria and microbiomes.  In addition to and beyond TSS determination, Cappable-seq depletes ribosomal RNA and reduces the complexity of the transcriptome to a single quantifiable tag per transcript enabling digital profiling of gene expression in any microbiome.

High-resolution nucleosome mapping of targeted regions using BAC-based enrichment.

We report a target enrichment method to map nucleosomes of large genomes at unprecedented coverage and resolution by deeply sequencing locus-specific mononucleosomal DNA enriched via hybridization with bacterial artificial chromosomes. We achieved ≈ 10 000-fold enrichment of specific loci, which enabled sequencing nucleosomes at up to ≈ 500-fold higher coverage than has been reported in a mammalian genome. We demonstrate the advantages of generating high-sequencing coverage for mapping the center of discrete nucleosomes, and we show the use of the method by mapping nucleosomes during T cell differentiation using nuclei from effector T-cells differentiated from clonal, isogenic, naïve, primary murine CD4 and CD8 T lymphocytes. The analysis reveals that discrete nucleosomes exhibit cell type-specific occupancy and positioning depending on differentiation status and transcription. This method is widely applicable to mapping many features of chromatin and discerning its landscape in large genomes at unprecedented resolution.

Nucleosome mapping across the CFTR locus identifies novel regulatory factors.

Nucleosome positioning on the chromatin strand plays a critical role in regulating accessibility of DNA to transcription factors and chromatin modifying enzymes. Hence, detailed information on nucleosome depletion or movement at cis-acting regulatory elements has the potential to identify predicted binding sites for trans-acting factors. Using a novel method based on enrichment of mononucleosomal DNA by bacterial artificial chromosome hybridization, we mapped nucleosome positions by deep sequencing across 250 kb, encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. CFTR shows tight tissue-specific regulation of expression, which is largely determined by cis-regulatory elements that lie outside the gene promoter. Although multiple elements are known, the repertoire of transcription factors that interact with these sites to activate or repress CFTR expression remains incomplete. Here, we show that specific nucleosome depletion corresponds to well-characterized binding sites for known trans-acting factors, including hepatocyte nuclear factor 1, Forkhead box A1 and CCCTC-binding factor. Moreover, the cell-type selective nucleosome positioning is effective in predicting binding sites for novel interacting factors, such as BAF155. Finally, we identify transcription factor binding sites that are overrepresented in regions where nucleosomes are depleted in a cell-specific manner. This approach recognizes the glucocorticoid receptor as a novel trans-acting factor that regulates CFTR expression in vivo.

The structural basis of mRNA recognition and binding by yeast pseudouridine synthase PUS1.

The chemical modification of RNA bases represents a ubiquitous activity that spans all domains of life. Pseudouridylation is the most common RNA modification and is observed within tRNA, rRNA, ncRNA and mRNAs. Pseudouridine synthase or 'PUS' enzymes include those that rely on guide RNA molecules and others that function as 'stand-alone' enzymes. Among the latter, several have been shown to modify mRNA transcripts. Although recent studies have defined the structural requirements for RNA to act as a PUS target, the mechanisms by which PUS1 recognizes these target sequences in mRNA are not well understood. Here we describe the crystal structure of yeast PUS1 bound to an RNA target that we identified as being a hot spot for PUS1-interaction within a model mRNA at 2.4 Å resolution. The enzyme recognizes and binds both strands in a helical RNA duplex, and thus guides the RNA containing the target uridine to the active site for subsequent modification of the transcript. The study also allows us to show the divergence of related PUS1 enzymes and their corresponding RNA target specificities, and to speculate on the basis by which PUS1 binds and modifies mRNA or tRNA substrates.

Genome-Wide Profiling of tRNA Using an Unexplored Reverse Transcriptase with High Processivity.

Monitoring the dynamic changes of cellular tRNA pools is challenging, due to the extensive post-transcriptional modifications of individual species. The most critical component in tRNAseq is a processive reverse transcriptase (RT) that can read through each modification with high efficiency. Here we show that the recently developed group-II intron RT Induro has the processivity and efficiency necessary to profile tRNA dynamics. Using our Induro-tRNAseq, simpler and more comprehensive than the best methods to date, we show that Induro progressively increases readthrough of tRNA over time and that the mechanism of increase is selective removal of RT stops, without altering the misincorporation frequency. We provide a parallel dataset of the misincorporation profile of Induro relative to the related TGIRT RT to facilitate the prediction of non-annotated modifications. We report an unexpected modification profile among human proline isoacceptors, absent from mouse and lower eukaryotes, that indicates new biology of decoding proline codons.

An RIG-I-Like RNA helicase mediates antiviral RNAi downstream of viral siRNA biogenesis in Caenorhabditis elegans.

Dicer ribonucleases of plants and invertebrate animals including Caenorhabditis elegans recognize and process a viral RNA trigger into virus-derived small interfering RNAs (siRNAs) to guide specific viral immunity by Argonaute-dependent RNA interference (RNAi). C. elegans also encodes three Dicer-related helicase (drh) genes closely related to the RIG-I-like RNA helicase receptors which initiate broad-spectrum innate immunity against RNA viruses in mammals. Here we developed a transgenic C. elegans strain that expressed intense green fluorescence from a chromosomally integrated flock house virus replicon only after knockdown or knockout of a gene required for antiviral RNAi. Use of the reporter nematode strain in a feeding RNAi screen identified drh-1 as an essential component of the antiviral RNAi pathway. However, RNAi induced by either exogenous dsRNA or the viral replicon was enhanced in drh-2 mutant nematodes, whereas exogenous RNAi was essentially unaltered in drh-1 mutant nematodes, indicating that exogenous and antiviral RNAi pathways are genetically distinct. Genetic epistatic analysis shows that drh-1 acts downstream of virus sensing and viral siRNA biogenesis to mediate specific antiviral RNAi. Notably, we found that two members of the substantially expanded subfamily of Argonautes specific to C. elegans control parallel antiviral RNAi pathways. These findings demonstrate both conserved and unique strategies of C. elegans in antiviral defense.

Expression of Human ACE2 N-terminal Domain, Part of the Receptor for SARS-CoV-2, in Fusion With Maltose-Binding Protein, E. coli Ribonuclease I and Human RNase A.

The SARS-CoV-2 viral genome contains a positive-strand single-stranded RNA of ∼30 kb. Human ACE2 protein is the receptor for SARS-CoV-2 virus attachment and infection. We propose to use ribonucleases (RNases) as antiviral agents to destroy the viral genome in vitro. In the virions, the RNA is protected by viral capsid proteins, membrane proteins, and nucleocapsid proteins. To utilize RNases as antiviral strategy, we set out to construct RNase fusion with human ACE2 receptor N-terminal domain (ACE2NTD). We expressed six proteins in E. coli cells: (1) MBP-ACE2NTD, (2) ACE2NTD-GFP, (3) RNase I (6×His), (4) RNase III (6×His), (5) RNase I-ACE2NTD (6×His), and (6) human RNase A-ACE2NTD (6×His). We evaluated fusion expression in different E. coli strains, partially purified MBP-ACE2NTD protein from the soluble fraction of bacterial cell lysate, and refolded MBP-ACE2NTD protein from inclusion body. The engineered RNase I-ACE2NTD (6×His) and hRNase A-ACE2NTD (6×His) fusions are active in cleaving SARS-CoV-2 RNA fragment in vitro. The recombinant RNase I (6×His) and RNase III (6×His) are active in cleaving RNA and dsRNA in test tube. This study provides a proof-of-concept for construction of fusion protein between human receptor and nuclease that may be used to degrade viral nucleic acids.

A Microbiome DNA Enrichment Method for Next-Generation Sequencing Sample Preparation.

"Microbiome" is used to describe the communities of microorganisms and their genes in a particular environment, including communities in association with a eukaryotic host or part of a host. One challenge in microbiome analysis concerns the presence of host DNA in samples. Removal of host DNA before sequencing results in greater sequence depth of the intended microbiome target population. This unit describes a novel method of microbial DNA enrichment in which methylated host DNA such as human genomic DNA is selectively bound and separated from microbial DNA before next-generation sequencing (NGS) library construction. This microbiome enrichment technique yields a higher fraction of microbial sequencing reads and improved read quality resulting in a reduced cost of downstream data generation and analysis. © 2016 by John Wiley & Sons, Inc.

Genome and metagenome sequencing: Using the human methyl-binding domain to partition genomic DNA derived from plant tissues.

PREMISE OF THE STUDY: Variation in the distribution of methylated CpG (methyl-CpG) in genomic DNA (gDNA) across the tree of life is biologically interesting and useful in genomic studies. We illustrate the use of human methyl-CpG-binding domain (MBD2) to fractionate angiosperm DNA into eukaryotic nuclear (methyl-CpG-rich) vs. organellar and prokaryotic (methyl-CpG-poor) elements for genomic and metagenomic sequencing projects. • METHODS: MBD2 has been used to enrich prokaryotic DNA in animal systems. Using gDNA from five model angiosperm species, we apply a similar approach to identify whether MBD2 can fractionate plant gDNA into methyl-CpG-depleted vs. enriched methyl-CpG elements. For each sample, three gDNA libraries were sequenced: (1) untreated gDNA, (2) a methyl-CpG-depleted fraction, and (3) a methyl-CpG-enriched fraction. • RESULTS: Relative to untreated gDNA, the methyl-depleted libraries showed a 3.2-11.2-fold and 3.4-11.3-fold increase in chloroplast DNA (cpDNA) and mitochondrial DNA (mtDNA), respectively. Methyl-enriched fractions showed a 1.8-31.3-fold and 1.3-29.0-fold decrease in cpDNA and mtDNA, respectively. • DISCUSSION: The application of MBD2 enabled fractionation of plant gDNA. The effectiveness was particularly striking for monocot gDNA (Poaceae). When sufficiently effective on a sample, this approach can increase the cost efficiency of sequencing plant genomes as well as prokaryotes living in or on plant tissues.

High-density nucleosome occupancy map of human chromosome 9p21-22 reveals chromatin organization of the type I interferon gene cluster.

Genome-wide investigations have dramatically increased our understanding of nucleosome positioning and the role of chromatin in gene regulation, yet some genomic regions have been poorly represented in human nucleosome maps. One such region is represented by human chromosome 9p21-22, which contains the type I interferon gene cluster that includes 16 interferon alpha genes and the single interferon beta, interferon epsilon, and interferon omega genes. A high-density nucleosome mapping strategy was used to generate locus-wide maps of the nucleosome organization of this biomedically important locus at a steady state and during a time course of infection with Sendai virus, an inducer of interferon gene expression. Detailed statistical and computational analysis illustrates that nucleosomes in this locus exhibit preferences for particular dinucleotide and oligomer DNA sequence motifs in vivo, which are similar to those reported for lower eukaryotic nucleosome-DNA interactions. These data were used to visualize the region's chromatin architecture and reveal features that are common to the organization of all the type I interferon genes, indicating a common nucleosome-mediated gene regulatory paradigm. Additionally, this study clarifies aspects of the dynamic changes that occur with the nucleosome occupying the transcriptional start site of the interferon beta gene after virus infection.

A method for selectively enriching microbial DNA from contaminating vertebrate host DNA.

DNA samples derived from vertebrate skin, bodily cavities and body fluids contain both host and microbial DNA; the latter often present as a minor component. Consequently, DNA sequencing of a microbiome sample frequently yields reads originating from the microbe(s) of interest, but with a vast excess of host genome-derived reads. In this study, we used a methyl-CpG binding domain (MBD) to separate methylated host DNA from microbial DNA based on differences in CpG methylation density. MBD fused to the Fc region of a human antibody (MBD-Fc) binds strongly to protein A paramagnetic beads, forming an effective one-step enrichment complex that was used to remove human or fish host DNA from bacterial and protistan DNA for subsequent sequencing and analysis. We report enrichment of DNA samples from human saliva, human blood, a mock malaria-infected blood sample and a black molly fish. When reads were mapped to reference genomes, sequence reads aligning to host genomes decreased 50-fold, while bacterial and Plasmodium DNA sequences reads increased 8-11.5-fold. The Shannon-Wiener diversity index was calculated for 149 bacterial species in saliva before and after enrichment. Unenriched saliva had an index of 4.72, while the enriched sample had an index of 4.80. The similarity of these indices demonstrates that bacterial species diversity and relative phylotype abundance remain conserved in enriched samples. Enrichment using the MBD-Fc method holds promise for targeted microbiome sequence analysis across a broad range of sample types.

Analysis of the C. elegans Argonaute family reveals that distinct Argonautes act sequentially during RNAi.

Argonaute (AGO) proteins interact with small RNAs to mediate gene silencing. C. elegans contains 27 AGO genes, raising the question of what roles these genes play in RNAi and related gene-silencing pathways. Here we describe 31 deletion alleles representing all of the previously uncharacterized AGO genes. Analysis of single- and multiple-AGO mutant strains reveals functions in several pathways, including (1) chromosome segregation, (2) fertility, and (3) at least two separate steps in the RNAi pathway. We show that RDE-1 interacts with trigger-derived sense and antisense RNAs to initiate RNAi, while several other AGO proteins interact with amplified siRNAs to mediate downstream silencing. Overexpression of downstream AGOs enhances silencing, suggesting that these proteins are limiting for RNAi. Interestingly, these AGO proteins lack key residues required for mRNA cleavage. Our findings support a two-step model for RNAi, in which functionally and structurally distinct AGOs act sequentially to direct gene silencing.

The dsRNA binding protein RDE-4 interacts with RDE-1, DCR-1, and a DExH-box helicase to direct RNAi in C. elegans.

Double-stranded (ds) RNA induces potent gene silencing, termed RNA interference (RNAi). At an early step in RNAi, an RNaseIII-related enzyme, Dicer (DCR-1), processes long-trigger dsRNA into small interfering RNAs (siRNAs). DCR-1 is also required for processing endogenous regulatory RNAs called miRNAs, but how DCR-1 recognizes its endogenous and foreign substrates is not yet understood. Here we show that the C. elegans RNAi pathway gene, rde-4, encodes a dsRNA binding protein that interacts during RNAi with RNA identical to the trigger dsRNA. RDE-4 protein also interacts in vivo with DCR-1, RDE-1, and a conserved DExH-box helicase. Our findings suggest a model in which RDE-4 and RDE-1 function together to detect and retain foreign dsRNA and to present this dsRNA to DCR-1 for processing.