Identification of antisense nucleic acid hybridization sites in mRNA molecules with self-quenching fluorescent reporter molecules

We describe a physical mRNA mapping strategy employing fluorescent self-quenching reporter molecules (SQRMs) that facilitates the identification of mRNA sequence accessible for hybridization with antisense nucleic acids in vitro and in vivo, real time. SQRMs are 20–30 base oligodeoxynucleotides with 5–6 bp complementary ends to which a 5′ fluorophore and 3′ quenching group are attached. Alone, the SQRM complementary ends form a stem that holds the fluorophore and quencher in contact. When the SQRM forms base pairs with its target, the structure separates the fluorophore from the quencher. This event can be reported by fluorescence emission when the fluorophore is excited. The stem–loop of the SQRM suggests that SQRM be made to target natural stem–loop structures formed during mRNA synthesis. The general utility of this method is demonstrated by SQRM identification of targetable sequence within c-myb and bcl-6 mRNA. Corresponding antisense oligonucleotides reduce these gene products in cells.     

Replication protein A modulates its interface with the primed DNA template during RNA-DNA primer elongation in replicating SV40 chromosomes

The eukaryal single-stranded DNA binding protein replication protein A (RPA) binds short oligonucleotides with high affinity but exhibits low cooperativity in binding longer templates, opposite to prokaryal counterparts. This discrepancy could reflect the smaller size of the replicative template portion availed to RPA. According to current models, this portion accommodates an RNA–DNA primer (RDP) of <40 nt (nested discontinuity) or a several-fold longer Okazaki fragment (initiation zone). Previous in situ UV-crosslinking revealed that RPA also interacts with nascent DNA, especially growing RDPs. Here we compare nascent SV40 DNA chains UV-crosslinked to the middle and large RPA subunits and use the data to re-examine the two models. The middle subunit interacted with the nascent chains after a few DNA residues were added to the RNA primer while the large subunit became accessible after extension by several more. Upon RDP maturation, the middle subunit disengaged while the large subunit remained accessible during further limited extension. A corresponding shift in preference in favor of the large subunit has been reported for purified RPA and synthetic gapped duplexes upon reduction of the gap from 19 to 9 nt. Combined, these facts support the proposal that the mature RDP faces downstream a correspondingly small gap, possibly created by removal of the RNA primer moiety from an adjacent, previously synthesized RDP (nested discontinuity) but insufficient for continuous elongation of the RDP into an Okazaki fragment (initiation zone).     

Thermodynamic criteria for high hit rate antisense oligonucleotide design

Antisense oligonucleotides are used for therapeutic applications and in functional genomic studies. In practice, however, many of the oligonucleotides complementary to an mRNA have little or no antisense activity. Theoretical strategies to improve the ‘hit rate’ in antisense screens will reduce the cost of discovery and may lead to identification of antisense oligonucleotides with increased potency. Statistical analysis performed on data collected from more than 1000 experiments with phosphorothioate-modified oligonucleotides revealed that the oligo-probes, which form stable duplexes with RNA (ΔG^o₃₇ ≤ –30 kcal/mol) and have small self-interaction potential, are more frequently efficient than molecules that form less stable oligonucleotide–RNA hybrids or more stable self-structures. To achieve optimal statistical preference, the values for self-interaction should be (ΔG^o₃₇) ≥ –8 kcal/mol for inter-oligonucleotide pairing and (ΔG^o₃₇) ≥ –1.1 kcal/mol for intra-molecular pairing. Selection of oligonucleotides with these thermodynamic values in the analyzed experiments would have increased the ‘hit rate’ by as much as 6-fold.     

Enhanced oligonucleotide-nanoparticle conjugate stability using thioctic acid modified oligonucleotides

Metallic nanoparticles of gold functionalized with oligonucleotides conventionally use a terminal thiol modification and have been used in a wide range of applications. Although readily available, the oligonucleotide–nanoparticle conjugates prepared in this way suffer from a lack of stability when exposed to a variety of small molecules or elevated temperatures. If silver is used in place of gold then this lack of stability is even more pronounced. In this study we report the synthesis of highly stabilized oligonucleotide–nanoparticle conjugates using a simple oligonucleotide modification. A modified solid support was used to generate 3′-thioctic acid modified oligonucleotides by treatment with an N-hydroxysuccimidyl ester of thioctic acid. Unusually, both gold and silver nanoparticles have been investigated in this study and show that these disulphide-modified oligonucleotide probes offer significant improvements in nanoparticle stability when treated with dithiothreitol (DTT) compared with monothiol analogues. This is a significant advance in oligonucleotide–nanoparticle conjugate stability and for the first time allows silver nanoparticles to be prepared that are more stable than standard gold-thiol functionalized nanoparticles. This opens up the possibility of using silver nanoparticles functionalized with oligonucleotides as an alternative to gold.     

In vitro selection of oligonucleotides that bind double-stranded DNA in the presence of triplex-stabilizing agents

A SELEX approach has been developed in order to select oligonucleotides that bind double-stranded DNA in the presence of a triplex-stabilizing agent, and was applied to a target sequence containing an oligopurine–oligopyrimidine stretch. After only seven rounds of selection, the process led to the identification of oligonucleotides that were able to form triple helices within the antiparallel motif. Inspection of the selected sequences revealed that, contrary to GC base pair which were always recognized by guanines, recognition of AT base pair could be achieved by either adenine or thymine, depending on the sequence context. While thymines are strongly preferred for several positions, some others can accommodate the presence of adenines. These results contribute to set the rules for designing oligonucleotides that form stable triple helices in the presence of triplex-stabilizing agents at physiological pH. They set the basis for further experiments regarding extension of potential target sequences for triple-helix formation or recognition of ligand–DNA complexes.     

Determination of optimal sites of antisense oligonucleotide cleavage within TNFalpha mRNA

Antisense oligonucleotides provide a powerful tool in order to determine the consequences of the reduced expression of a selected target gene and may include target validation and therapeutic applications. Methods of predicting optimum antisense sites are not always effective. We have compared the efficacy of antisense oligonucleotides, which were selected in vitro using random combinatorial oligonucleotide libraries of differing length and complexity, upon putative target sites within TNFα mRNA. The relationship of specific target site accessibility and oligonucleotide efficacy with respect to these parameters proved to be complex. Modification of the length of the recognition sequence of the oligonucleotide library illustrated that independent target sites demonstrated a preference for antisense oligonucleotides of a defined and independent optimal length. The efficacy of antisense oligonucleotide sequences selected in vitro paralleled that observed in phorbol 12-myristate 13-acetate (PMA)-activated U937 cells. The application of methylphosphonate:phosphodiester chimaeric oligonucleotides to U937 cells reduced mRNA levels to up to 19.8% that of the untreated cell population. This approach provides a predictive means to profile any mRNA of known sequence with respect to the identification and optimisation of sites accessible to antisense oligonucleotide activity.     

Osprey: a comprehensive tool employing novel methods for the design of oligonucleotides for DNA sequencing and microarrays

We have developed a software package called Osprey for the calculation of optimal oligonucleotides for DNA sequencing and the creation of microarrays based on either PCR-products or directly spotted oligomers. It incorporates a novel use of position-specific scoring matrices, for the sensitive and specific identification of secondary binding sites anywhere in the target sequence. Using accelerated hardware is faster and more efficient than the traditional pairwise alignments used in most oligo-design software. Osprey consists of a module for target site selection based on user input, novel utilities for dealing with problematic sequences such as repeats, and a common code base for the identification of optimal oligonucleotides from the target list. Overall, these improvements provide a program that, without major increases in run time, reflects current DNA thermodynamics models, improves specificity and reduces the user's data preprocessing and parameterization requirements. Using a TimeLogic™ hardware accelerator, we report up to 50-fold reduction in search time versus a linear search strategy. Target sites may be derived from computer analysis of DNA sequence assemblies in the case of sequencing efforts, or genome or EST analysis in the case of microarray development in both prokaryotes and eukaryotes.     

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.     

Development and comparison of procedures for the selection of delta ribozyme cleavage sites within the hepatitis B virus

Delta ribozyme possesses several unique features related to the fact that it is the only catalytic RNA known to be naturally active in human cells. This makes it attractive as a therapeutic tool for the inactivation of clinically relevant RNAs. However, several hurdles must be overcome prior to the development of useful gene-inactivation systems based on delta ribozyme. We have developed three procedures for the selection of potential delta ribozyme target sites within the hepatitis B virus (HBV) pregenome: (i) the use of bioinformatic tools coupled to biochemical assays; (ii) RNase H hydrolysis with a pool of oligonucleotides; and (iii) cleavage assays with a pool of ribozymes. The results obtained with delta ribozyme show that these procedures are governed by several rules, some of which are different from those both for other catalytic RNAs and antisense oligonucleotides. Together, these procedures identified 12 sites in the HBV pregenome that can be cleaved by delta ribozymes, although with different efficiencies. Clearly, both target site accessibility and the ability to form an active ribozyme–substrate complex constitute interdependent factors that can best be addressed using a combinatorial library of either oligonucleotides or ribozymes.     

Analysis of a nuclease activity of catalytic domain of Thermus thermophilus MutS2 by high-accuracy mass spectrometry

Electrospray ionization with Fourier-transform ion cyclotron resonance mass spectrometry (ESI–FT ICR MS) is a powerful tool for analyzing the precise structural features of biopolymers, including oligonucleotides. Here, we described the detailed characterization of a newly discovered nuclease activity of the C-terminal domain of Thermus thermophilus MutS2 (ttMutS2). Using this method, the length, nucleotide content and nature of the 5′- and 3′-termini of the product oligonucleotides were accurately identified. It is revealed that the C-terminal domain of ttMutS2 incised the phosphate backbone of oligodeoxynucleotides non-sequence-specifically at the 3′ side of the phosphates. The simultaneous identification of the innumerable fragments was achieved by the extremely high-accuracy of ESI–FT ICR MS.     

Clamping down on weak terminal base pairs: oligonucleotides with molecular caps as fidelity-enhancing elements at the 5'- and 3'-terminal residues

The base-pairing fidelity of oligonucleotides depends on the identity of the nucleobases involved and the position of matched or mismatched base pairs in the duplex. Nucleobases forming weak base pairs, as well as a terminal position favor mispairing. We have searched for 5′-appended acylamido caps that enhance the stability and base-pairing fidelity of oligonucleotides with a 5′-terminal 2′-deoxyadenosine residue using combinatorial synthesis and MALDI-monitored nuclease selections. This provided the residue of 4-(pyren-1-yl)butyric acid as a lead. Lead optimization gave (S)-N-(pyren-1-ylmethyl)pyrrolidine-3-phosphate as a cap that increases duplex stability and base-pairing fidelity. For the duplex of 5′-AGGTTGAC-3′ with its fully complementary target, this cap gives an increase in the UV melting point T_m of +10.9°C. The T_m is 6.3–8.3°C lower when a mismatched nucleobase faces the 5′-terminal dA residue. The optimized cap can be introduced via automated DNA synthesis. It was combined with an anthraquinone carboxylic acid residue as a cap for the 3′-terminal residue. A doubly capped dodecamer thus prepared gives a melting point decrease for double-terminal mismatches that is 5.7–5.9°C greater than that for the unmodified control duplex.     

A gene-specific DNA sequencing chip for exploring molecular evolutionary change

Sequencing by hybridization (SBH) approaches to DNA sequencing face two conflicting constraints. First, in order to ensure that the target DNA binds reliably, the oligonucleotide probes that are attached to the chip array must be >15 bp in length. Secondly, the total number of possible 15 bp oligonucleotides is too large (>4¹⁵) to fit on a chip with current technology. To circumvent the conflict between these two opposing constraints, we present a novel gene-specific DNA chip design. Our design is based on the idea that not all conceivable oligonucleotides need to be placed on a chip— only those that capture sequence combinations occurring in nature. Our approach uses a training set of aligned sequences that code for the gene in question. We compute the minimum number of oligonucleotides (generally 15–30 bp in length) that need to be placed on a DNA chip to capture the variation implied by the training set using a graph search algorithm. We tested the approach in silico using cytochrome-b sequences. Results indicate that on average, 98% of the sequence of an unknown target can be determined using the approach.     

Identification of the First ATRIP–Deficient Patient and Novel Mutations in ATR Define a Clinical Spectrum for ATR–ATRIP Seckel Syndrome

A homozygous mutational change in the Ataxia-Telangiectasia and RAD3 related (ATR) gene was previously reported in two related families displaying Seckel Syndrome (SS). Here, we provide the first identification of a Seckel Syndrome patient with mutations in ATRIP, the gene encoding ATR–Interacting Protein (ATRIP), the partner protein of ATR required for ATR stability and recruitment to the site of DNA damage. The patient has compound heterozygous mutations in ATRIP resulting in reduced ATRIP and ATR expression. A nonsense mutational change in one ATRIP allele results in a C-terminal truncated protein, which impairs ATR–ATRIP interaction; the other allele is abnormally spliced. We additionally describe two further unrelated patients native to the UK with the same novel, heterozygous mutations in ATR, which cause dramatically reduced ATR expression. All patient-derived cells showed defective DNA damage responses that can be attributed to impaired ATR–ATRIP function. Seckel Syndrome is characterised by microcephaly and growth delay, features also displayed by several related disorders including Majewski (microcephalic) osteodysplastic primordial dwarfism (MOPD) type II and Meier-Gorlin Syndrome (MGS). The identification of an ATRIP–deficient patient provides a novel genetic defect for Seckel Syndrome. Coupled with the identification of further ATR–deficient patients, our findings allow a spectrum of clinical features that can be ascribed to the ATR–ATRIP deficient sub-class of Seckel Syndrome. ATR–ATRIP patients are characterised by extremely severe microcephaly and growth delay, microtia (small ears), micrognathia (small and receding chin), and dental crowding. While aberrant bone development was mild in the original ATR–SS patient, some of the patients described here display skeletal abnormalities including, in one patient, small patellae, a feature characteristically observed in Meier-Gorlin Syndrome. Collectively, our analysis exposes an overlapping clinical manifestation between the disorders but allows an expanded spectrum of clinical features for ATR–ATRIP Seckel Syndrome to be defined.     

Improved nearest-neighbor parameters for the stability of RNA/DNA hybrids under a physiological condition

The stability of Watson–Crick paired RNA/DNA hybrids is important for designing optimal oligonucleotides for ASO (Antisense Oligonucleotide) and CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)–Cas9 techniques. Previous nearest-neighbour (NN) parameters for predicting hybrid stability in a 1 M NaCl solution, however, may not be applicable for predicting stability at salt concentrations closer to physiological condition (e.g. ∼100 mM Na⁺ or K⁺ in the presence or absence of Mg²⁺). Herein, we report measured thermodynamic parameters of 38 RNA/DNA hybrids at 100 mM NaCl and derive new NN parameters to predict duplex stability. Predicted ΔG°₃₇ and T_m values based on the established NN parameters agreed well with the measured values with 2.9% and 1.1°C deviations, respectively. The new results can also be used to make precise predictions for duplexes formed in 100 mM KCl or 100 mM NaCl in the presence of 1 mM Mg²⁺, which can mimic an intracellular and extracellular salt condition, respectively. Comparisons of the predicted thermodynamic parameters with published data using ASO and CRISPR–Cas9 may allow designing shorter oligonucleotides for these techniques that will diminish the probability of non-specific binding and also improve the efficiency of target gene regulation.     

Predicting future community-level ocular Chlamydia trachomatis infection prevalence using serological,clinical, molecular,and geospatial data

Trachoma is an infectious disease characterized by repeated exposures to Chlamydia trachomatis (Ct) that may ultimately lead to blindness. Efficient identification of communities with high infection burden could help target more intensive control efforts. We hypothesized that IgG seroprevalence in combination with geospatial layers, machine learning, and model-based geostatistics would be able to accurately predict future community-level ocular Ct infections detected by PCR. We used measurements from 40 communities in the hyperendemic Amhara region of Ethiopia to assess this hypothesis. Median Ct infection prevalence among children 0–5 years old increased from 6% at enrollment, in the context of recent mass drug administration (MDA), to 29% by month 36, following three years without MDA. At baseline, correlation between seroprevalence and Ct infection was stronger among children 0–5 years old (ρ = 0.77) than children 6–9 years old (ρ = 0.48), and stronger than the correlation between active trachoma and Ct infection (0-5y ρ = 0.56; 6-9y ρ = 0.40). Seroprevalence was the strongest concurrent predictor of infection prevalence at month 36 among children 0–5 years old (cross-validated R² = 0.75, 95% CI: 0.58–0.85), though predictive performance declined substantially with increasing temporal lag between predictor and outcome measurements. Geospatial variables, a spatial Gaussian process, and stacked ensemble machine learning did not meaningfully improve predictions. Serological markers among children 0–5 years old may be an objective tool for identifying communities with high levels of ocular Ct infections, but accurate, future prediction in the context of changing transmission remains an open challenge.     

A Mixed-Method Approach for Quantifying Illegal Fishing and Its Impact on an Endangered Fish Species

Illegal harvest is recognized as a widespread problem in natural resource management. The use of multiple methods for quantifying illegal harvest has been widely recommended yet infrequently applied. We used a mixed-method approach to evaluate the extent, character, and motivations of illegal gillnet fishing in Lake Hovsgol National Park, Mongolia and its impact on the lake’s fish populations, especially that of the endangered endemic Hovsgol grayling (Thymallus nigrescens). Surveys for derelict fishing gear indicate that gillnet fishing is widespread and increasing and that fishers generally use 3–4 cm mesh gillnet. Interviews with resident herders and park rangers suggest that many residents fish for subsistence during the spring grayling spawning migration and that some residents fish commercially year-round. Interviewed herders and rangers generally agree that fish population sizes are decreasing but are divided on the causes and solutions. Biological monitoring indicates that the gillnet mesh sizes used by fishers efficiently target Hovsgol grayling. Of the five species sampled in the monitoring program, only burbot (Lota lota) showed a significant decrease in population abundance from 2009–2013. However, grayling, burbot, and roach (Rutilus rutilus) all showed significant declines in average body size, suggesting a negative fishing impact. Data-poor stock assessment methods suggest that the fishing effort equivalent to each resident family fishing 50-m of gillnet 11–15 nights per year would be sufficient to overexploit the grayling population. Results from the derelict fishing gear survey and interviews suggest that this level of effort is not implausible. Overall, we demonstrate the ability for a mixed-method approach to effectively describe an illegal fishery and suggest that these methods be used to assess illegal fishing and its impacts in other protected areas.     

Long insert whole genome sequencing for copy number variant and translocation detection

As next-generation sequencing continues to have an expanding presence in the clinic, the identification of the most cost-effective and robust strategy for identifying copy number changes and translocations in tumor genomes is needed. We hypothesized that performing shallow whole genome sequencing (WGS) of 900–1000-bp inserts (long insert WGS, LI-WGS) improves our ability to detect these events, compared with shallow WGS of 300–400-bp inserts. A priori analyses show that LI-WGS requires less sequencing compared with short insert WGS to achieve a target physical coverage, and that LI-WGS requires less sequence coverage to detect a heterozygous event with a power of 0.99. We thus developed an LI-WGS library preparation protocol based off of Illumina’s WGS library preparation protocol and illustrate the feasibility of performing LI-WGS. We additionally applied LI-WGS to three separate tumor/normal DNA pairs collected from patients diagnosed with different cancers to demonstrate our application of LI-WGS on actual patient samples for identification of somatic copy number alterations and translocations. With the evolution of sequencing technologies and bioinformatics analyses, we show that modifications to current approaches may improve our ability to interrogate cancer genomes.     

Plant microRNA-Target Interaction Identification Model Based on the Integration of Prediction Tools and Support Vector Machine

Background

Confident identification of microRNA-target interactions is significant for studying the function of microRNA (miRNA). Although some computational miRNA target prediction methods have been proposed for plants, results of various methods tend to be inconsistent and usually lead to more false positive. To address these issues, we developed an integrated model for identifying plant miRNA–target interactions.

Results

Three online miRNA target prediction toolkits and machine learning algorithms were integrated to identify and analyze Arabidopsis thaliana miRNA-target interactions. Principle component analysis (PCA) feature extraction and self-training technology were introduced to improve the performance. Results showed that the proposed model outperformed the previously existing methods. The results were validated by using degradome sequencing supported Arabidopsis thaliana miRNA-target interactions. The proposed model constructed on Arabidopsis thaliana was run over Oryza sativa and Vitis vinifera to demonstrate that our model is effective for other plant species.

Conclusions

The integrated model of online predictors and local PCA-SVM classifier gained credible and high quality miRNA-target interactions. The supervised learning algorithm of PCA-SVM classifier was employed in plant miRNA target identification for the first time. Its performance can be substantially improved if more experimentally proved training samples are provided.     

Comparison of different antisense strategies in mammalian cells using locked nucleic acids, 2'-O-methyl RNA,phosphorothioates and small interfering RNA

Locked nucleic acids (LNAs) and double-stranded small interfering RNAs (siRNAs) are rather new promising antisense molecules for cell culture and in vivo applications. Here, we compare LNA–DNA–LNA gapmer oligonucleotides and siRNAs with a phosphorothioate and a chimeric 2′-O-methyl RNA–DNA gapmer with respect to their capacities to knock down the expression of the vanilloid receptor subtype 1 (VR1). LNA–DNA–LNA gapmers with four or five LNAs on either side and a central stretch of 10 or 8 DNA monomers in the center were found to be active gapmers that inhibit gene expression. A comparative co-transfection study showed that siRNA is the most potent inhibitor of VR1–green fluorescent protein (GFP) expression. A specific inhibition was observed with an estimated IC₅₀ of 0.06 nM. An LNA gapmer was found to be the most efficient single-stranded antisense oligonucleotide, with an IC₅₀ of 0.4 nM being 175-fold lower than that of commonly used phosphorothioates (IC₅₀ ~70 nM). In contrast, the efficiency of a 2′-O-methyl-modified oligonucleotide (IC₅₀ ~220 nM) was 3-fold lower compared with the phosphorothioate. The high potency of siRNAs and chimeric LNA–DNA oligonucleotides make them valuable candidates for cell culture and in vivo applications targeting the VR1 mRNA.     

Computational generation and screening of RNA motifs in large nucleotide sequence pools

Although identification of active motifs in large random sequence pools is central to RNA in vitro selection, no systematic computational equivalent of this process has yet been developed. We develop a computational approach that combines target pool generation, motif scanning and motif screening using secondary structure analysis for applications to 10¹²–10¹⁴-sequence pools; large pool sizes are made possible using program redesign and supercomputing resources. We use the new protocol to search for aptamer and ribozyme motifs in pools up to experimental pool size (10¹⁴ sequences). We show that motif scanning, structure matching and flanking sequence analysis, respectively, reduce the initial sequence pool by 6–8, 1–2 and 1 orders of magnitude, consistent with the rare occurrence of active motifs in random pools. The final yields match the theoretical yields from probability theory for simple motifs and overestimate experimental yields, which constitute lower bounds, for aptamers because screening analyses beyond secondary structure information are not considered systematically. We also show that designed pools using our nucleotide transition probability matrices can produce higher yields for RNA ligase motifs than random pools. Our methods for generating, analyzing and designing large pools can help improve RNA design via simulation of aspects of in vitro selection.