Sensitive and accurate analysis of gene expression signatures enabled by oligonucleotide-labelled cDNA.

High-throughput RNA sequencing offers a comprehensive analysis of transcriptome complexity originated from regulatory events, such as differential gene expression, alternative polyadenylation and others, and allows the increase in diagnostic capacity and precision. For gene expression profiling applications that do not specifically require information on alternative splicing events, the mRNA 3' termini counting approach is a cost-effective alternative to whole transcriptome sequencing. Here, we report MTAS-seq (mRNA sequencing via terminator-assisted synthesis) - a novel RNA-seq library preparation method directed towards mRNA 3' termini. We demonstrate the specific enrichment for 3'-terminal regions by simple and quick single-tube protocol with built-in molecular barcoding to enable accurate estimation of transcript abundance. To achieve that, we synthesized oligonucleotide-modified dideoxynucleotides which enable the generation of cDNA libraries at the reverse transcription step. We validated the performance of MTAS-seq on well-characterized reference bulk RNA and further tested it with eukaryotic cell lysates.

UbaLAI is a monomeric Type IIE restriction enzyme.

Type II restriction endonucleases (REases) form a large and highly diverse group of enzymes. Even REases specific for a common recognition site often vary in their oligomeric structure, domain organization and DNA cleavage mechanisms. Here we report biochemical and structural characterization of the monomeric restriction endonuclease UbaLAI, specific for the pseudosymmetric DNA sequence 5'-CC/WGG-3' (where W = A/T, and '/' marks the cleavage position). We present a 1.6 Å co-crystal structure of UbaLAI N-terminal domain (UbaLAI-N) and show that it resembles the B3-family domain of EcoRII specific for the 5'-CCWGG-3' sequence. We also find that UbaLAI C-terminal domain (UbaLAI-C) is closely related to the monomeric REase MvaI, another enzyme specific for the 5'-CCWGG-3' sequence. Kinetic studies of UbaLAI revealed that it requires two recognition sites for optimal activity, and, like other type IIE enzymes, uses one copy of a recognition site to stimulate cleavage of a second copy. We propose that during the reaction UbaLAI-N acts as a handle that tethers the monomeric UbaLAI-C domain to the DNA, thereby helping UbaLAI-C to perform two sequential DNA nicking reactions on the second recognition site during a single DNA-binding event. A similar reaction mechanism may be characteristic to other monomeric two-domain REases.

Screening for catalytically active Type II restriction endonucleases using segregation-induced methylation deficiency.

Type II restriction endonucleases (REases) are one of the basic tools of recombinant DNA technology. They also serve as models for elucidation of mechanisms for both site-specific DNA recognition and cleavage by proteins. However, isolation of catalytically active mutants from their libraries is challenging due to the toxicity of REases in the absence of protecting methylation, and techniques explored so far had limited success. Here, we present an improved SOS induction-based approach for in vivo screening of active REases, which we used to isolate a set of active variants of the catalytic mutant, Cfr10I(E204Q). Detailed characterization of plasmids from 64 colonies screened from the library of ∼200,000 transformants revealed 29 variants of cfr10IR gene at the level of nucleotide sequence and 15 variants at the level of amino acid sequence, all of which were able to induce SOS response. Specific activity measurements of affinity-purified mutants revealed >200-fold variance among them, ranging from 100% (wild-type isolates) to 0.5% (S188C mutant), suggesting that the technique is equally suited for screening of mutants possessing high or low activity and confirming that it may be applied for identification of residues playing a role in catalysis.

Three-stage biochemical selection: cloning of prototype class IIS/IIC/IIG restriction endonuclease-methyltransferase TsoI from the thermophile Thermus scotoductus.

BACKGROUND: In continuing our research into the new family of bifunctional restriction endonucleases (REases), we describe the cloning of the tsoIRM gene. Currently, the family includes six thermostable enzymes: TaqII, Tth111II, TthHB27I, TspGWI, TspDTI, TsoI, isolated from various Thermus sp. and two thermolabile enzymes: RpaI and CchII, isolated from mesophilic bacteria Rhodopseudomonas palustris and Chlorobium chlorochromatii, respectively. The enzymes have several properties in common. They are large proteins (molecular size app. 120 kDa), coded by fused genes, with the REase and methyltransferase (MTase) in a single polypeptide, where both activities are affected by S-adenosylmethionine (SAM). They recognize similar asymmetric cognate sites and cleave at a distance of 11/9 nt from the recognition site. Thus far, we have cloned and characterised TaqII, Tth111II, TthHB27I, TspGWI and TspDTI. RESULTS: TsoI REase, which originate from thermophilic Thermus scotoductus RFL4 (T. scotoductus), was cloned in Escherichia coli (E. coli) using two rounds of biochemical selection of the T. scotoductus genomic library for the TsoI methylation phenotype. DNA sequencing of restriction-resistant clones revealed the common open reading frame (ORF) of 3348 bp, coding for a large polypeptide of 1116 aminoacid (aa) residues, which exhibited a high level of similarity to Tth111II (50% identity, 60% similarity). The ORF was PCR-amplified, subcloned into a pET21 derivative under the control of a T7 promoter and was subjected to the third round of biochemical selection in order to isolate error-free clones. Induction experiments resulted in synthesis of an app. 125 kDa protein, exhibiting TsoI-specific DNA cleavage. Also, the wild-type (wt) protein was purified and reaction optima were determined. CONCLUSIONS: Previously we identified and cloned the Thermus family RM genes using a specially developed method based on partial proteolysis of thermostable REases. In the case of TsoI the classic biochemical selection method was successful, probably because of the substantially lower optimal reaction temperature of TsoI (app. 10-15°C). That allowed for sufficient MTase activity in vivo in recombinant E. coli. Interestingly, TsoI originates from bacteria with a high optimum growth temperature of 67°C, which indicates that not all bacterial enzymes match an organism's thermophilic nature, and yet remain functional cell components. Besides basic research advances, the cloning and characterisation of the new prototype REase from the Thermus sp. family enzymes is also of practical importance in gene manipulation technology, as it extends the range of available DNA cleavage specificities.

Fusion sequencing via terminator-assisted synthesis (FTAS-seq) identifies TMPRSS2 fusion partners in prostate cancer.

Genetic rearrangements that fuse an androgen-regulated promoter area with a protein-coding portion of an originally androgen-unaffected gene are frequent in prostate cancer, with the fusion between transmembrane serine protease 2 (TMPRSS2) and ETS transcription factor ERG (ERG) (TMPRSS2-ERG fusion) being the most prevalent. Conventional hybridization- or amplification-based methods can test for the presence of expected gene fusions, but the exploratory analysis of currently unknown fusion partners is often cost-prohibitive. Here, we developed an innovative next-generation sequencing (NGS)-based approach for gene fusion analysis termed fusion sequencing via terminator-assisted synthesis (FTAS-seq). FTAS-seq can be used to enrich the gene of interest while simultaneously profiling the whole spectrum of its 3'-terminal fusion partners. Using this novel semi-targeted RNA-sequencing technique, we were able to identify 11 previously uncharacterized TMPRSS2 fusion partners and capture a range of TMPRSS2-ERG isoforms. We tested the performance of FTAS-seq with well-characterized prostate cancer cell lines and utilized the technique for the analysis of patient RNA samples. FTAS-seq chemistry combined with appropriate primer panels holds great potential as a tool for biomarker discovery that can support the development of personalized cancer therapies.

Related bifunctional restriction endonuclease-methyltransferase triplets: TspDTI, Tth111II/TthHB27I and TsoI with distinct specificities.

BACKGROUND: We previously defined a family of restriction endonucleases (REases) from Thermus sp., which share common biochemical and biophysical features, such as the fusion of both the nuclease and methyltransferase (MTase) activities in a single polypeptide, cleavage at a distance from the recognition site, large molecular size, modulation of activity by S-adenosylmethionine (SAM), and incomplete cleavage of the substrate DNA. Members include related thermophilic REases with five distinct specificities: TspGWI, TaqII, Tth111II/TthHB27I, TspDTI and TsoI. RESULTS: TspDTI, TsoI and isoschizomers Tth111II/TthHB27I recognize different, but related sequences: 5'-ATGAA-3', 5'-TARCCA-3' and 5'-CAARCA-3' respectively. Their amino acid sequences are similar, which is unusual among REases of different specificity. To gain insight into this group of REases, TspDTI, the prototype member of the Thermus sp. enzyme family, was cloned and characterized using a recently developed method for partially cleaving REases. CONCLUSIONS: TspDTI, TsoI and isoschizomers Tth111II/TthHB27I are closely related bifunctional enzymes. They comprise a tandem arrangement of Type I-like domains, like other Type IIC enzymes (those with a fusion of a REase and MTase domains), e.g. TspGWI, TaqII and MmeI, but their sequences are only remotely similar to these previously characterized enzymes. The characterization of TspDTI, a prototype member of this group, extends our understanding of sequence-function relationships among multifunctional restriction-modification enzymes.

Advanced preparation of fragment libraries enabled by oligonucleotide-modified 2',3'-dideoxynucleotides.

The ever-growing demand for inexpensive, rapid, and accurate exploration of genomes calls for refinement of existing sequencing techniques. The development of next-generation sequencing (NGS) was a revolutionary milestone in genome analysis. While modified nucleotides already were inherent tools in sequencing and imaging, further modification of nucleotides enabled the expansion into even more diverse applications. Herein we describe the design and synthesis of oligonucleotide-tethered 2',3'-dideoxynucleotide (ddONNTP) terminators bearing universal priming sites attached to the nucleobase, as well as their enzymatic incorporation and performance in read-through assays. In the context of NGS library preparation, the incorporation of ddONNTP fulfills two requirements at once: the fragmentation step is integrated into the workflow and the obtained fragments are readily labeled by platform-specific adapters. DNA polymerases can incorporate ddONNTP nucleotides, as shown by primer extension assays. More importantly, reading through the unnatural linkage during DNA synthesis was demonstrated, with 25-30% efficiency in single-cycle extension.

Fluoride-cleavable, fluorescently labelled reversible terminators: synthesis and use in primer extension.

Fluorescent 2'-deoxynucleotides containing a protecting group at the 3'-O-position are reversible terminators that enable array-based DNA sequencing-by-synthesis (SBS) approaches. Herein, we describe the synthesis and full characterisation of four reversible terminators bearing a 3'-blocking moiety and a linker-dye system that is removable under the same fluoride-based treatment. Each nucleotide analogue has a different fluorophore attached to the base through a fluoride-cleavable linker and a 2-cyanoethyl moiety as the 3'-blocking group, which can be removed by using a fluoride treatment as well. Furthermore, we identified a DNA polymerase, namely, RevertAid M-MuLV reverse transcriptase, which can incorporate the four modified reversible terminators. The synthesised nucleotides and the optimised DNA polymerase were used on CodeLink slides spotted with hairpin oligonucleotides to demonstrate their potential in a cyclic reversible terminating approach.

Enzymatic Synthesis of Chimeric DNA Oligonucleotides by in Vitro Transcription with dTTP, dCTP, dATP, and 2'-Fluoro Modified dGTP.

Efficient ways to produce single-stranded DNA are of great interest for diverse applications in molecular biology and nanotechnology. In the present study, we selected T7 RNA polymerase mutants with reduced substrate specificity to employ an in vitro transcription reaction for the synthesis of chimeric DNA oligonucleotides, either individually or in pools. We performed in vitro evolution based on fluorescence-activated droplet sorting and identified mutations V783M, V783L, V689Q, and G555L as novel variants leading to relaxed substrate discrimination. Transcribed chimeric oligonucleotides were tested in PCR, and the quality of amplification products as well as fidelity of oligonucleotide synthesis were assessed by NGS. We concluded that enzymatically produced chimeric DNA transcripts contain significantly fewer deletions and insertions compared to chemically synthesized counterparts and can successfully serve as PCR primers, making the evolved enzymes superior for simple and cheap one-pot synthesis of multiple chimeric DNA oligonucleotides in parallel using a plethora of premixed templates.

MuA-based Molecular Indexing for Rare Mutation Detection by Next-Generation Sequencing.

Detection of low-frequency mutations in cancer genomes or other heterogeneous cell populations requires high-fidelity sequencing. Molecular barcoding is one of the key technologies that enables the differentiation of true mutations from errors, which can be caused by sequencing or library preparation processes. However, current approaches where barcodes are introduced via primer extension or adaptor ligation do not utilize the full power of barcoding, due to complicated library preparation workflows and biases. Here we demonstrate the remarkable tolerance of MuA transposase to the presence of multiple replacements in transposon sequence, and explore this unique feature to engineer the MuA transposome complex with randomised nucleotides in 12 transposon positions, which can be introduced as a barcode into the target molecule after transposition event. We applied the approach of Unique MuA-based Molecular Indexing (UMAMI) to assess the power of rare mutation detection by shortgun sequencing on the Illumina platform. Our results show that UMAMI allows detection of rare mutations readily and reliably, and in this paper we report error rate values for the number of thermophilic DNA polymerases measured by using UMAMI.

High-resolution microbiome analysis enabled by linking of 16S rRNA gene sequences with adjacent genomic contexts.

Sequence-based characterization of bacterial communities has long been a hostage of limitations of both 16S rRNA gene and whole metagenome sequencing. Neither approach is universally applicable, and the main efforts to resolve constraints have been devoted to improvement of computational prediction tools. Here, we present semi-targeted 16S rRNA sequencing (st16S-seq), a method designed for sequencing V1-V2 regions of the 16S rRNA gene along with the genomic locus upstream of the gene. By in silico analysis of 13 570 bacterial genome assemblies, we show that genome-linked 16S rRNA sequencing is superior to individual hypervariable regions or full-length gene sequences in terms of classification accuracy and identification of gene copy numbers. Using mock communities and soil samples we experimentally validate st16S-seq and benchmark it against the established microbial classification techniques. We show that st16S-seq delivers accurate estimation of 16S rRNA gene copy numbers, enables taxonomic resolution at the species level and closely approximates community structures obtainable by whole metagenome sequencing.

Bioinformatic and partial functional analysis of pEspA and pEspB, two plasmids from Exiguobacterium arabatum sp. nov. RFL1109.

The complete nucleotide sequences of two plasmids from Exiguobacterium arabatum sp. nov. RFL1109, pEspA (4563bp) and pEspB (38,945bp), have been determined. Five ORFs were identified in the pEspA plasmid, and putative functions were assigned to two of them. Using deletion mapping approach, the Rep-independent replication region of pEspA, which functions in Bacillus subtilis, was localized within a 0.6kb DNA region. Analysis of the pEspB sequence revealed 42 ORFs. From these, function of two genes encoding enzymes of the Lsp1109I restriction-modification system was confirmed experimentally, while putative functions of another 18 ORFs were suggested based on comparative analysis. Three functional regions have been proposed for the pEspB plasmid: the putative conjugative transfer region, the region involved in plasmid replication and maintenance, and the region responsible for transposition of the IS21 family-like transposable elements.

Author Correction: Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles.

The original version of this Article contained errors in the units of concentration of three reagents listed in the Methods. These errors have all been corrected in both the PDF and HTML versions of the Article.

Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles.

Microbial single-cell genomics can be used to provide insights into the metabolic potential, interactions, and evolution of uncultured microorganisms. Here we present WGA-X, a method based on multiple displacement amplification of DNA that utilizes a thermostable mutant of the phi29 polymerase. WGA-X enhances genome recovery from individual microbial cells and viral particles while maintaining ease of use and scalability. The greatest improvements are observed when amplifying high G+C content templates, such as those belonging to the predominant bacteria in agricultural soils. By integrating WGA-X with calibrated index-cell sorting and high-throughput genomic sequencing, we are able to analyze genomic sequences and cell sizes of hundreds of individual, uncultured bacteria, archaea, protists, and viral particles, obtained directly from marine and soil samples, in a single experiment. This approach may find diverse applications in microbiology and in biomedical and forensic studies of humans and other multicellular organisms.Single-cell genomics can be used to study uncultured microorganisms. Here, Stepanauskas et al. present a method combining improved multiple displacement amplification and FACS, to obtain genomic sequences and cell size information from uncultivated microbial cells and viral particles in environmental samples.

MnlI--The member of H-N-H subtype of Type IIS restriction endonucleases.

The Type IIS restriction endonuclease MnlI recognizes the non-palindromic nucleotide sequence 5'-CCTC(N)7/6 downward arrow and cleaves DNA strands as indicated by the arrow. The genes encoding MnlI restriction-modification system were cloned and sequenced. It comprises N6-methyladenine and C5-methylcytosine methyltransferases and the restriction endonuclease. Biochemical studies revealed that MnlI restriction endonuclease cleaves double- and single-stranded DNA, and that it prefers different metal ions for hydrolysis of these substrates. Mg2+ ions were shown to be required for the specific cleavage of double-stranded DNA, whereas Ni2+ and some other transition metal ions were preferred for nonspecific cleavage of single-stranded DNA. The C-terminal part of MnlI restriction endonuclease revealed an intriguing similarity with the H-N-H type nucleolytic domain of bacterial toxins, Colicin E7 and Colicin E9. Alanine replacements in the conserved sequence motif 306Rx3ExHHx14Nx8H greatly reduced specific activity of MnlI, and some mutations even completely inactivated the enzyme. However, none of these mutations had effect on MnlI binding to the specific DNA, and on its oligomerisation state as well. We interpret the presented experimental evidence as a suggestion that the motif 306Rx3ExHHx14Nx8H represents the active site of MnlI. Consequentially, MnlI seems to be the member of Type IIS with the active site of the H-N-H type.

Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease.

The restriction endonuclease (REase) R. HphI is a Type IIS enzyme that recognizes the asymmetric target DNA sequence 5'-GGTGA-3' and in the presence of Mg(2+) hydrolyzes phosphodiester bonds in both strands of the DNA at a distance of 8 nucleotides towards the 3' side of the target, producing a 1 nucleotide 3'-staggered cut in an unspecified sequence at this position. REases are typically ORFans that exhibit little similarity to each other and to any proteins in the database. However, bioinformatics analyses revealed that R.HphI is a member of a relatively big sequence family with a conserved C-terminal domain and a variable N-terminal domain. We predict that the C-terminal domains of proteins from this family correspond to the nuclease domain of the HNH superfamily rather than to the most common PD-(D/E)XK superfamily of nucleases. We constructed a three-dimensional model of the R.HphI catalytic domain and validated our predictions by site-directed mutagenesis and studies of DNA-binding and catalytic activities of the mutant proteins. We also analyzed the genomic neighborhood of R.HphI homologs and found that putative nucleases accompanied by a DNA methyltransferase (i.e. predicted REases) do not form a single group on a phylogenetic tree, but are dispersed among free-standing putative nucleases. This suggests that nucleases from the HNH superfamily were independently recruited to become REases in the context of RM systems multiple times in the evolution and that members of the HNH superfamily may be much more frequent among the so far unassigned REase sequences than previously thought.

Esp1396I restriction-modification system: structural organization and mode of regulation.

Esp1396I restriction-modification (RM) system recognizes an interrupted palindromic DNA sequence 5'-CCA(N)(5)TGG-3'. The Esp1396I RM system was found to reside on pEsp1396, a 5.6 kb plasmid naturally occurring in Enterobacter sp. strain RFL1396. The nucleotide sequence of the entire 5622 bp pEsp1396 plasmid was determined on both strands. Identified genes for DNA methyltransferase (esp1396IM) and restriction endonuclease (esp1396IR) are transcribed convergently. The restriction endonuclease gene is preceded by the small ORF (esp1396IC) that possesses a strong helix-turn-helix motif and resembles regulatory proteins found in PvuII, BamHI and few other RM systems. Gene regulation studies revealed that C.Esp1396I acts as both a repressor of methylase expression and an activator of regulatory protein and restriction endonuclease expression. Our data indicate that C protein from Esp1396I RM system activates the expression of the Enase gene, which is co-transcribed from the promoter of regulatory gene, by the mechanism of coupled translation.

Circular permutation of DNA cytosine-N4 methyltransferases: in vivo coexistence in the BcnI system and in vitro probing by hybrid formation.

Sequence analysis of the BcnI restriction-modification system from Bacillus centrosporus revealed four open reading frames (bcnIC, bcnIR, bcnIB and bcnIA) that are arranged as two converging collinear pairs. One pair encodes a putative small regulatory protein, C.BcnI, and the restriction endonuclease R.BcnI. The other two gene products are the DNA cytosine-N4 methyltransferases M.BcnIA and M.BcnIB, which differ by circular permutation of conserved sequence motifs. The BcnI methyltransferases are isospecific on double-stranded DNA [methylation specificity CC(C/G)GG], but M.BcnIA can also methylate the target sites in single-stranded DNA. Functional analysis shows that bcnIA is dispensable (bcnIB is capable of protecting the DNA against the in vivo activity of bcnIR); in contrast, no stable clones were obtained if bcnIB alone was deleted from the system. By analogy with the DpnII system, the second methylase M.BcnIA may play a role in the transformation proficiency of its gram-positive host. The interchangeability of homologous elements in the beta class of cytosine-N4 methylases was probed by hybrid formation between M.BcnIB and its closest homolog M.Cfr9I (CCCGGG) employing a novel semi-random strategy combined with selection for catalytic activity. The fusion points in the active hybrids mapped in a narrow region located between sequence motifs X and I. Our data illustrate that recombination of two related sequences by circular permutation may serve as an evolutionary mechanism for creating new specificities of amino MTases.

Engineering of restriction endonucleases: using methylation activity of the bifunctional endonuclease Eco57I to select the mutant with a novel sequence specificity.

Type II restriction endonucleases (REs) are widely used tools in molecular biology, biotechnology and diagnostics. Efforts to generate new specificities by structure-guided design and random mutagenesis have been unsuccessful so far. We have developed a new procedure called the methylation activity-based selection (MABS) for generating REs with a new specificity. MABS uses a unique property of bifunctional type II REs to methylate DNA targets they recognize. The procedure includes three steps: (1) conversion of a bifunctional RE into a monofunctional DNA-modifying enzyme by cleavage center disruption; (2) mutagenesis and selection of mutants with altered DNA modification specificity based on their ability to protect predetermined DNA targets; (3) reconstitution of the cleavage center's wild-type structure. The efficiency of the MABS technique was demonstrated by altering the sequence specificity of the bifunctional RE Eco57I from 5'-CTGAAG to 5'-CTGRAG, and thus generating the mutant restriction endonuclease (and DNA methyltransferase) of a specificity not known before. This study provides evidence that MABS is a promising technique for generation of REs with new specificities.