Structural bias in T4 RNA ligase-mediated 3'-adapter ligation

T4 RNA ligases are commonly used to attach adapters to RNAs, but large differences in ligation efficiency make detection and quantitation problematic. We developed a ligation selection strategy using random RNAs in combination with high-throughput sequencing to gain insight into the differences in efficiency of ligating pre-adenylated DNA adapters to RNA 3'-ends. After analyzing biases in RNA sequence, secondary structure and RNA-adapter cofold structure, we conclude that T4 RNA ligases do not show significant primary sequence preference in RNA substrates, but are biased against structural features within RNAs and adapters. Specifically, RNAs with less than three unstructured nucleotides at the 3'-end and RNAs that are predicted to cofold with an adapter in unfavorable structures are likely to be poorly ligated. The effect of RNA-adapter cofold structures on ligation is supported by experiments where the ligation efficiency of specific miRNAs was changed by designing adapters to alter cofold structure. In addition, we show that using adapters with randomized regions results in higher ligation efficiency and reduced ligation bias. We propose that using randomized adapters may improve RNA representation in experiments that include a 3'-adapter ligation step.     

Quantitative bias in Illumina TruSeq and a novel post amplification barcoding strategy for multiplexed DNA and small RNA deep sequencing

Here we demonstrate a method for unbiased multiplexed deep sequencing of RNA and DNA libraries using a novel, efficient and adaptable barcoding strategy called Post Amplification Ligation-Mediated (PALM). PALM barcoding is performed as the very last step of library preparation, eliminating a potential barcode-induced bias and allowing the flexibility to synthesize as many barcodes as needed. We sequenced PALM barcoded micro RNA (miRNA) and DNA reference samples and evaluated the quantitative barcode-induced bias in comparison to the same reference samples prepared using the Illumina TruSeq barcoding strategy. The Illumina TruSeq small RNA strategy introduces the barcode during the PCR step using differentially barcoded primers, while the TruSeq DNA strategy introduces the barcode before the PCR step by ligation of differentially barcoded adaptors. Results show virtually no bias between the differentially barcoded miRNA and DNA samples, both for the PALM and the TruSeq sample preparation methods. We also multiplexed miRNA reference samples using a pre-PCR barcode ligation. This barcoding strategy results in significant bias.     

Mismatch discrimination and sequence bias during end-joining by DNA ligases

DNA ligases, critical enzymes for in vivo genome maintenance and modern molecular biology, catalyze the joining of adjacent 3′-OH and 5′-phosphorylated ends in DNA. To determine whether DNA annealing equilibria or properties intrinsic to the DNA ligase enzyme impact end-joining ligation outcomes, we used a highly multiplexed, sequencing-based assay to profile mismatch discrimination and sequence bias for several ligases capable of efficient end-joining. Our data reveal a spectrum of fidelity and bias, influenced by both the strength of overhang annealing as well as sequence preferences and mismatch tolerances that vary both in degree and kind between ligases. For example, while T7 DNA ligase shows a strong preference for ligating high GC sequences, other ligases show little GC-dependent bias, with human DNA Ligase 3 showing almost none. Similarly, mismatch tolerance varies widely among ligases, and while all ligases tested were most permissive of G:T mismatches, some ligases also tolerated bulkier purine:purine mismatches. These comprehensive fidelity and bias profiles provide insight into the biology of end-joining reactions and highlight the importance of ligase choice in application design.     

Bias in Ligation-Based Small RNA Sequencing Library Construction Is Determined by Adaptor and RNA Structure

High-throughput sequencing (HTS) has become a powerful tool for the detection of and sequence characterization of microRNAs (miRNA) and other small RNAs (sRNA). Unfortunately, the use of HTS data to determine the relative quantity of different miRNAs in a sample has been shown to be inconsistent with quantitative PCR and Northern Blot results. Several recent studies have concluded that the major contributor to this inconsistency is bias introduced during the construction of sRNA libraries for HTS and that the bias is primarily derived from the adaptor ligation steps, specifically where single stranded adaptors are sequentially ligated to the 3’ and 5’-end of sRNAs using T4 RNA ligases. In this study we investigated the effects of ligation bias by using a pool of randomized ligation substrates, defined mixtures of miRNA sequences and several combinations of adaptors in HTS library construction. We show that like the 3’ adaptor ligation step, the 5’ adaptor ligation is also biased, not because of primary sequence, but instead due to secondary structures of the two ligation substrates. We find that multiple secondary structural factors influence final representation in HTS results. Our results provide insight about the nature of ligation bias and allowed us to design adaptors that reduce ligation bias and produce HTS results that more accurately reflect the actual concentrations of miRNAs in the defined starting material.     

Conserved residues in domain Ia are required for the reaction of Escherichia coli DNA ligase with NAD+.

NAD(+)-dependent DNA ligases are present in all bacteria and are essential for growth. Their unique substrate specificity compared with ATP-dependent human DNA ligases recommends the NAD(+) ligases as targets for the development of new broad-spectrum antibiotics. A plausible strategy for drug discovery is to identify the structural components of bacterial DNA ligase that interact with NAD(+) and then to isolate small molecules that recognize these components and thereby block the binding of NAD(+) to the ligase. The limitation to this strategy is that the structural determinants of NAD(+) specificity are not known. Here we show that reactivity of Escherichia coli DNA ligase (LigA) with NAD(+) requires N-terminal domain Ia, which is unique to, and conserved among, NAD(+) ligases but absent from ATP-dependent ligases. Deletion of domain Ia abolished the sealing of 3'-OH/5'-PO(4) nicks and the reaction with NAD(+) to form ligase-adenylate but had no effect on phosphodiester formation at a preadenylated nick. Alanine substitutions at conserved residues within domain Ia either reduced (His-23, Tyr-35) or abolished (Tyr-22, Asp-32, Asp-36) sealing of a 5'-PO(4) nick and adenylyl transfer from NAD(+) without affecting ligation of pre-formed DNA-adenylate. We suggest that these five side chains comprise a binding site for the nicotinamide mononucleotide moiety of NAD(+). Structure-activity relationships were clarified by conservative substitutions.     

RNA沉默在植物生物逆境反应中的作用

RNA沉默是真核生物共有的基因表达调节机制和防御机制。在植物RNA沉默中, 一些小RNAs, 如微小 RNAs和小干扰RNAs, 在植物防御病毒、细菌或食草动物的反应中具有重要作用。为了抑制宿主的RNA沉默系统, 植物病毒或细菌进化出了在RNA沉默不同阶段起作用的病毒沉默抑制子或细菌沉默抑制子, 来克服寄主的RNA沉默反应。文章就植物RNA沉默、病毒沉默抑制子、细菌沉默抑制子及其相关防御反应的一些新进展做一概述。     

Elimination of Ligation Dependent Artifacts in T4 RNA Ligase to Achieve High Efficiency and Low Bias MicroRNA Capture

Adapter ligation is a critical first step in many microRNA analysis methods including microarray, qPCR, and sequencing. Previous studies have shown that ligation bias can have dramatic effects on both the fidelity of expression profiles and reproducibility across samples. We have developed a method for high efficiency and low bias microRNA capture by 3′ adapter ligation using T4 RNA ligase that does not require pooled adapters. Using a panel of 20 microRNA, we investigated the effects of ligase type, PEG concentration, ligase amount, adapter concentration, incubation time, incubation temperature, and adapter design on capture efficiency and bias. Of these factors, high PEG% was found to be critical in suppressing ligation bias. We obtained high average capture efficiency and low CV across the 20 microRNA panel, both in idealized buffer conditions (86%±10%) and total RNA spiking conditions (64%±17%). We demonstrate that this method is reliable across microRNA species that previous studies have had difficulty capturing and that our adapter design performs significantly better than the common adapter designs. Further, we demonstrate that the optimization methodology must be specifically designed for minimizing bias in order to obtain the ideal reaction parameters.     

Assessing an Improved Protocol for Plasma microRNA Extraction

The first step in biomarkers discovery is to identify the best protocols for their purification and analysis. This issue is critical when considering peripheral blood samples (plasma and serum) that are clinically interesting but meet several methodological problems, mainly complexity and low biomarker concentration. Analysis of small molecules, such as circulating microRNAs, should overcome these disadvantages. The present study describes an optimal RNA extraction method of microRNAs from human plasma samples. Different reagents and commercially available kits have been analyzed, identifying also the best pre-analytical conditions for plasma isolation. Between all of them, the column-based approaches were shown to be the most effective. In this context, miRNeasy Serum/Plasma Kit (from Qiagen) rendered more concentrated RNA, that was better suited for microarrays studies and did not require extra purification steps for sample concentration and purification than phenol based extraction methods. We also present evidences that the addition of low doses of an RNA carrier before starting the extraction process improves microRNA purification while an already published carrier dose can result in significant bias over microRNA profiles. Quality controls for best protocol selection were developed by spectrophotometry measurement of contaminants and microfluidics electrophoresis (Agilent 2100 Bioanalyzer) for RNA integrity. Selected donor and patient plasma samples and matched biopsies were tested by Affymetrix microarray technology to compare differentially expressed microRNAs. In summary, this study defines an optimized protocol for microRNA purification from human blood samples, increasing the performance of assays and shedding light over the best way to discover and use these biomarkers in clinical practice.     

RNA capping enzyme and DNA ligase: a superfamily of covalent nucleotidyl transferases

mRNA capping entails GMP transfer from GTP to a 5' diphosphate RNA end to form the structure G(5')ppp(5')N. A similar reaction involving AMP transfer to the 5' monophosphate end of DNA or RNA occurs during strand joining by polynucleotide ligases. In both cases, nucleotidyl transfer occurs through a covalent lysyl-NMP intermediate. Sequence conservation among capping enzymes and ATP-dependent ligases in the vicinity of the active site lysine (KxDG) and at five other co-linear motifs suggests a common structural basis for covalent catalysis. Mutational studies support this view. We propose that the cellular and DNA virus capping enzymes and ATP-dependent ligases constitute a protein superfamily evolved from a common ancestral enzyme. Within this superfamily, the cellular capping enzymes display more extensive similarity to the ligases than they do to the poxvirus capping enzymes. Recent studies suggest that eukaryotic RNA viruses have evolved alternative pathways of cap metabolism catalysed by structurally unrelated enzymes that nonetheless employ a phosphoramidate intermediate. Comparative analysis of these enzymes, particularly at the structural level, should illuminate the shared reaction mechanism while clarifying the basis for nucleotide specificity and end recognition. The capping enzymes merit close attention as potential targets for antiviral therapy.     

Assessment of Whole Genome Amplification for Sequence Capture and Massively Parallel Sequencing

Exome sequence capture and massively parallel sequencing can be combined to achieve inexpensive and rapid global analyses of the functional sections of the genome. The difficulties of working with relatively small quantities of genetic material, as may be necessary when sharing tumor biopsies between collaborators for instance, can be overcome using whole genome amplification. However, the potential drawbacks of using a whole genome amplification technology based on random primers in combination with sequence capture followed by massively parallel sequencing have not yet been examined in detail, especially in the context of mutation discovery in tumor material. In this work, we compare mutations detected in sequence data for unamplified DNA, whole genome amplified DNA, and RNA originating from the same tumor tissue samples from 16 patients diagnosed with non-small cell lung cancer. The results obtained provide a comprehensive overview of the merits of these techniques for mutation analysis. We evaluated the identified genetic variants, and found that most (74%) of them were observed in both the amplified and the unamplified sequence data. Eighty-nine percent of the variations found by WGA were shared with unamplified DNA. We demonstrate a strategy for avoiding allelic bias by including RNA-sequencing information.     

WWP1 E3 Ligase Targets LATS1 for Ubiquitin-Mediated Degradation in Breast Cancer Cells

The Large Tumor Suppressor 1 (LATS1) is a serine/threonine kinase and tumor suppressor found down-regulated in various human cancers. LATS1 has recently been identified as a central player of the emerging Hippo signaling pathway, which plays important roles in organ size control, tumorigenesis, and stem cell differentiation and renewal, etc. Although mounting evidence supports a role of LATS1 in tumor suppression and tumorigenesis, how LATS1 is regulated at the molecular level is not fully understood. Recently several positive regulators of LATS1 (Mst1/2, MOB1, Kibra, etc) have been identified but how LATS1 is negatively regulated is still largely unknown. We have recently identified Itch, a member of the NEDD4-like family E3 ubiquitin ligases, as a novel negative regulator of LATS1. However, whether other ubiquitin ligases modulate LATS1 stability and function is unclear. By screening many E3 ligases of the NEDD4-like family using over-expression and short-interference RNA knockdown approaches, we have identified WWP1 E3 ligase as another novel negative regulator of LATS1. We have provided in vitro and in vivo evidence that WWP1 is essential for LATS1 stability and negatively regulate LATS1 by promoting LATS1 degradation through polyubiquitination and the 26S proteasome pathway. Importantly, we also showed that degradation of LATS1 is critical in mediating WWP1-induced increased cell proliferation in breast cancer cells. Since WWP1 is an oncogene and LATS1 is a tumor suppressor gene in breast cancer, our studies provide a promising therapeutic strategy in which developed drugs targeting WWP1 cause activation of LATS1 in suppressing breast cancer cell growth.     

外泌体在小RNA基因治疗中的研究进展

小RNA,包括小干扰RNA以及微小RNA,已成为多种疾病的潜在治疗药物。目前,小RNA的运输载体主要是病毒或者合成试剂。然而这类载体往往毒性高且特异性低。外泌体是由内源细胞分泌出来的天然纳米材料,本身能够穿越生物膜并在细胞间传递小RNA。以外泌体为基础的小RNA递送作为一种新的转运方式,能够克服低效率,低特异性以及免疫反应等缺陷,有望成为新型载体。本文简要论述了以外泌体为载体的小RNA递送系统在临床治疗研究中的前沿进展。     

Structure and mechanism of RNA ligase

T4 RNA ligase 2 (Rnl2) exemplifies an RNA ligase family that includes the RNA editing ligases (RELs) of Trypanosoma and Leishmania. The Rnl2/REL enzymes are defined by essential signature residues and a unique C-terminal domain, which we show is essential for sealing of 3'-OH and 5'-PO4 RNA ends by Rnl2, but not for ligase adenylation or phosphodiester bond formation at a preadenylated AppRNA end. The N-terminal segment Rnl2(1-249) of the 334 aa Rnl2 protein comprises an autonomous adenylyltransferase/AppRNA ligase domain. We report the 1.9 A crystal structure of the ligase domain with AMP bound at the active site, which reveals a shared fold, catalytic mechanism, and evolutionary history for RNA ligases, DNA ligases, and mRNA capping enzymes.     

Transcriptional activation of microRNA-34a by NF-kappa B in human esophageal cancer cells

Background

RNA ligases are essential reagents for many methods in molecular biology including NextGen RNA sequencing. To prevent ligation of RNA to itself, ATP independent mutant ligases, defective in self-adenylation, are often used in combination with activated pre-adenylated linkers. It is important that these ligases not have de-adenylation activity, which can result in activation of RNA and formation of background ligation products. An additional useful feature is for the ligase to be active at elevated temperatures. This has the advantage or reducing preferences caused by structures of single-stranded substrates and linkers.

Results

To create an RNA ligase with these desirable properties we performed mutational analysis of the archaeal thermophilic RNA ligase from Methanobacterium thermoautotrophicum. We identified amino acids essential for ATP binding and reactivity but dispensable for phosphodiester bond formation with 5’ pre-adenylated donor substrate. The motif V lysine mutant (K246A) showed reduced activity in the first two steps of ligation reaction. The mutant has full ligation activity with pre-adenylated substrates but retained the undesirable activity of deadenylation, which is the reverse of step 2 adenylation. A second mutant, an alanine substitution for the catalytic lysine in motif I (K97A) abolished activity in the first two steps of the ligation reaction, but preserved wild type ligation activity in step 3. The activity of the K97A mutant is similar with either pre-adenylated RNA or single-stranded DNA (ssDNA) as donor substrates but we observed two-fold preference for RNA as an acceptor substrate compared to ssDNA with an identical sequence. In contrast, truncated T4 RNA ligase 2, the commercial enzyme used in these applications, is significantly more active using pre-adenylated RNA as a donor compared to pre-adenylated ssDNA. However, the T4 RNA ligases are ineffective in ligating ssDNA acceptors.

Conclusions

Mutational analysis of the heat stable RNA ligase from Methanobacterium thermoautotrophicum resulted in the creation of an ATP independent ligase. The K97A mutant is defective in the first two steps of ligation but retains full activity in ligation of either RNA or ssDNA to a pre-adenylated linker. The ability of the ligase to function at 65°C should reduce the constraints of RNA secondary structure in RNA ligation experiments.     

Structural and Functional Analysis of Viral siRNAs

A large amount of short interfering RNA (vsiRNA) is generated from plant viruses during infection, but the function, structure and biogenesis of these is not understood. We profiled vsiRNAs using two different high-throughput sequencing platforms and also developed a hybridisation based array approach. The profiles obtained through the Solexa platform and by hybridisation were very similar to each other but different from the 454 profile. Both deep sequencing techniques revealed a strong bias in vsiRNAs for the positive strand of the virus and identified regions on the viral genome that produced vsiRNA in much higher abundance than other regions. The hybridisation approach also showed that the position of highly abundant vsiRNAs was the same in different plant species and in the absence of RDR6. We used the Terminator 5′-Phosphate-Dependent Exonuclease to study the 5′ end of vsiRNAs and showed that a perfect control duplex was not digested by the enzyme without denaturation and that the efficiency of the Terminator was strongly affected by the concentration of the substrate. We found that most vsiRNAs have 5′ monophosphates, which was also confirmed by profiling short RNA libraries following either direct ligation of adapters to the 5′ end of short RNAs or after replacing any potential 5′ ends with monophosphates. The Terminator experiments also showed that vsiRNAs were not perfect duplexes. Using a sensor construct we also found that regions from the viral genome that were complementary to non-abundant vsiRNAs were targeted in planta just as efficiently as regions recognised by abundant vsiRNAs. Different high-throughput sequencing techniques have different reproducible sequence bias and generate different profiles of short RNAs. The Terminator exonuclease does not process double stranded RNA, and because short RNAs can quickly re-anneal at high concentration, this assay can be misleading if the substrate is not denatured and not analysed in a dilution series. The sequence profiles and Terminator digests suggest that CymRSV siRNAs are produced from the structured positive strand rather than from perfect double stranded RNA or by RNA dependent RNA polymerase.