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.

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' m⁷G 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' m⁷G 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 m⁷G caps and open chromatin, such as CAGE and DNase-HS seq, confirmed the validity of these high-confidence RNA scaffolds.

Expression and Purification of BsaXI Restriction Endonuclease and Engineering New Specificity From BsaXI Specificity Subunit

It is stated that BsaXI is a Type IIB restriction endonuclease (REase) that cleaves both sides of its recognition sequence 5′↓N9 AC N5 CTCC N10↓ 3′ (complement strand 5′ ↓N7 GGAG N5 GT N12↓ 3′), creating 3-base 3′ overhangs. Here we report the cloning and expression of bsaXIS and bsaXIRM genes in Escherichia coli. The BsaXI activity was successfully reconstituted by mixing the BsaXI RM fusion subunit with the BsaXI S subunit and the enzyme complex further purified by chromatography over 6 columns. As expected, the S subunit consisted of two subdomains encoding TRD1-CR1 [target recognition domain (TRD), conserved region (CR)] for 5′ AC 3′, and TRD2-CR2 presumably specifying 5′ CTCC 3′. TRD1-CR1 (TRD2-CR2 deletion) or duplication of TRD1 (TRD1-CR1-TRD1-CR2) both generated a new specificity 5′ AC N5 GT 3′ when the S variants were complexed with the RM subunits. The circular permutation of TRD1 and TRD2, i.e., the relocation of TRD2-CR2 to the N-terminus and TRD1-CR1 to the C-terminus generated the same specificity with the RM subunits, although some wobble cleavage was detected. The TRD2 domain in the BsaXI S subunit can be substituted by a close homolog (∼59% sequence identity) and generated the same specificity. However, TRD2-CR2 domain alone failed to express in E. coli, but CR1-TRD2-CR2 protein could be expressed and purified which showed partial nicking activity with the RM subunits. This work demonstrated that like Type I restriction systems, the S subunit of a Type IIB system could also be manipulated to create new specificities. The genome mining of BsaXI TRD2 homologs in GenBank found more than 36 orphan TRD2 homologs, implying that quite a few orphan TRD2s are present in microbial genomes that may be potentially paired with other TRDs to create new restriction specificities.