An automated algorithm for sequence confirmation of chemically modified oligonucleotides by tandem mass spectrometry

We have developed a tandem mass spectrometry (MS/MS) data analysis program for confirmation of sequence of chemically modified oligonucleotides. The method is based on the analysis of deconvoluted MS/MS data for fragment ions from three charge states and comparison of these data against a set of computer-generated masses from expected fragmentation patterns. The algorithm compares the experimental masses not only against the fragment set predicted for the expected sequence but also against a wider test set covering all next-neighbor position switches of the original sequence and all pairwise swaps of nucleosides, which in synthesis would result in molecules with masses within a preset mass tolerance. The algorithm is capable of identifying incorrect sequences that would not be distinguished by identity testing with electrospray ionization mass spectrometry. The method has been tested with permutations of the two 21-mer single strands of a chemically modified short interfering RNA containing 2′-O-methyl and phosphorothioate linkages. For both strands, challenge sequences were synthesized and tested with the premise that they were the original sequences. The algorithm correctly reported the locations of next-neighbor position switches and nucleoside swaps. The results confirm the approach as useful for MS/MS-based identity test methods for synthetic oligonucleotides.     

反义寡核苷酸药物癌泰得的定性分析方法研究

 癌泰得 (ACTCACTCAGGCCTCAGACT)为端粒酶表达抑制活性反义寡核苷酸 .为了探讨其定性检测手段 ,通过阴离子交换高效液相色谱和毛细管凝胶电泳分析方法 ,确定了该硫代寡核苷酸以及与其有关的短序列和部分未被硫代类似物的保留时间 ,并分析了不同混合物样品 .结果表明 ,阴离子交换高效液相色谱对硫代寡核苷酸骨架上的差异非常敏感 ,可很好地分离长度相同的硫代和未完全硫代类似物 ,并且随未被硫代磷酸基数目增加 ,保留时间依次缩短 .阴离子交换高效液相色谱对硫代寡核苷酸的长度不敏感 ,不能分离相差一个碱基的硫代寡核苷酸 .毛细管凝胶电泳可很好地分离长度相差一个碱基的硫代寡核苷酸 ,不能分离同长硫代和部分硫代寡核苷酸 .高效液相色谱结合毛细管凝胶电泳可有效地确定癌泰得的纯度和修饰程度 .     

Microarrays for identifying binding sites and probing structure of RNAs

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure–function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.     

Oligonucleotide mapping: evaluation of its sensitivity by computer-simulation.

A frequently used method of comparing large RNA molecules employs the two-dimensional display of oligonucleotides generated through the action of specific RNases (oligonucleotide mapping, fingerprinting). Using computer simulations and simple analytic expressions the number of large RNase T1-resistant oligonucleotides obtained from random RNA sequences can be estimated. The computer simulations also permit estimation of the number of large oligonucleotides which remain unchanged as random variations are introduced into a random RNA sequence. In addition, computer analysis also provides a means of estimating statistical confidence limits to be used in a quantitative comparison of fingerprints of different RNA molecules. The model shows that two RNA sequences which differ overall by 1%, 5% or 10% share, on average, only 85%, 50% or 25%, respectively, of their large oligonucleotides. Thus, the use of fingerprint analysis is recommended only when closely related RNAs or regions of RNAs are compared (sequence homology greater than 90%).     

Large multiple organism gene finding by collapsed Gibbs sampling.

The Gibbs sampling method has been widely used for sequence analysis after it was successfully applied to the problem of identifying regulatory motif sequences upstream of genes. Since then, numerous variants of the original idea have emerged: however, in all cases the application has been to finding short motifs in collections of short sequences (typically less than 100 nucleotides long). In this paper, we introduce a Gibbs sampling approach for identifying genes in multiple large genomic sequences up to hundreds of kilobases long. This approach leverages the evolutionary relationships between the sequences to improve the gene predictions, without explicitly aligning the sequences. We have applied our method to the analysis of genomic sequence from 14 genomic regions, totaling roughly 1.8 Mb of sequence in each organism. We show that our approach compares favorably with existing ab initio approaches to gene finding, including pairwise comparison based gene prediction methods which make explicit use of alignments. Furthermore, excellent performance can be obtained with as little as four organisms, and the method overcomes a number of difficulties of previous comparison based gene finding approaches: it is robust with respect to genomic rearrangements, can work with draft sequence, and is fast (linear in the number and length of the sequences). It can also be seamlessly integrated with Gibbs sampling motif detection methods.     

Reducing the number of microlocations in oligonucleotide microchip matrices by the application of degenerate oligonucleotides

The application of degenerate oligonucleotides to DNA Sequencing by Hybridisation with Oligonucleotide Matrix (SHOM) is proposed. The use of degenerate oligonucleotides is regarded as an example of pooling methods that are suitable for various laboratory procedures requiring numerous samples to be assayed. As each DNA sequence coded by four letters (A, G, C, T) may be defined by two sequences: a sequence coded by W and S (W-weak-A or T, S-strong-G or C) and a sequence coded by R and Y (R-purine-A or G, Y-pirymidine-T or C), n4n -nucleotide sequences may be defined with the help of 2xn2sequences. In the place of the originally described microchip matrix composed of all possible unambiguous octanucleotides (4(8)=65 536) attached to the equal number of 65 536 microlocations a matrix composed of 512 microlocations containing 256 2(8)-degenerate octanucleotides is proposed. The matrix contains all 256 possible octanucleotides coded by W and S variations and all 256 possible octanucleotides coded by R and Y variations. The 512 256-degenerate octanucleotides allows to retrieve the same information as 65 536 unambiguous octanucleotides. A variant of the DNA sequence reconstruction method applicable to this system is presented. The use of degenerate oligonucleotides also gives the possibility to apply matrices composed of longer oligonucleotides without increasing the number of microlocations in matrices, which would enable increasing the length of unambiguously reconstructed sequence, e.g. a matrix comprising 131 072 16-mer oligonucleotides i.e. 65 536 65 536-fold degenerate oligonucleotide coded by W and S variations and 65 536 65 536-fold degenerate oligonucleotide coded by R and Y variations could replace one matrix comprising all possible unambiguous 16-mer oligonucleotides (ca. 4.3x10(9)).     

Design, construction, and validation of a modular library of sequence diversity standards for polymerase chain reaction

Methods to measure the sequence diversity of polymerase chain reaction (PCR)-amplified DNA lack standards for use as assay calibrators and controls. Here we present a general and economical method for developing customizable DNA standards of known sequence diversity. Standards ranging from 1 to 25,000 sequences were generated by directional ligation of oligonucleotide “words” of standard length and GC content and then amplified by PCR. The sequence accuracy and diversity of the library were validated using AmpliCot analysis (DNA hybridization kinetics) and Illumina sequencing. The library has the following features: (i) pools containing tens of thousands of sequences can be generated from the ligation of relatively few commercially synthesized short oligonucleotides; (ii) each sequence differs from all others in the library at a minimum of three nucleotide positions, permitting discrimination between different sequences by either sequencing or hybridization; (iii) all sequences have identical length, GC content, and melting temperature; (iv) the identity of each standard can be verified by restriction digestion; and (v) once made, the ends of the library may be cleaved and replaced with sequences to match any PCR primer pair. These standards should greatly improve the accuracy and reproducibility of sequence diversity measurements.     

Recent developments in methods for RNA sequencing using in vitro 32P-labeling

A variety of approaches that utilize in vitro 32P-labeling of RNA and of oligonucleotides in the sequence analysis of RNAs are described. These include 1) methods for 5'- and 3'- end labeling of RNAs; 2) end labeling and sequencing of oligonucleotides present in complete T1 RNase or pancreatic RNase digests of RNA; 3) use of random endonucleases, such as nuclease P1, for terminal sequence analysis of end labeled RNAs; and 4) use of base specific enzymes or chemical reagents in the sequence analysis of end-labeled RNAs. Also described is an approach to RNA sequencing, applied so far to tRNAs, which is based on partial and random alkaline cleavage of an RNA to generate a series of overlapping oligonucleotide fragments, all containing the original 3'-end of the RNA. Analysis of the 5'- end group of each of these oligonucleotides (following 5'-end labeling with 32P) provides the sequence of most of the tRNA. The above methods have been used to derive the sequences of several tRNAs, the ribosomal 5S and 5 x 8S RNAs, a viroid RNA, and large segments of both prokaryotic and eukaryotic ribosomal and messenger RNAs.     

Linguistic Measure of Taxonomic and Functional Relatedness of Nucleotide Sequences

Abstract

The frequencies of “words”, oligonucleotides within nucleotide sequences, reflect the genetic information contained in the sequence “texts”. Nucleotide sequences are characteristically represented by their contrast word vocabularies. Comparison of the sequences by correlating their contrast vocabularies is shown to reflect well the relatedness (unrelatedness) between the sequences. A single value, the linguistic similarity between the sequences, is suggested asa measure of sequence relatedness. Sequences as short as 1000 bases can be characterized and quantitatively related to other sequences by this technique. The linguistic sequence similarity value is used for analysis of taxonomically and functionally diverse nucleotide sequences. The similarity value is shown to be very sensitive to the relatedness of the source species, thus providing a convenient tool for taxonomic classification of species by their sequence vocabularies. Functionally diverse sequences appear distinct by their linguistic similarity values. This can be a basis for a quick screening technique for functional characterization of the sequences and for mapping functionally distinct regions in long sequences.     

Sensitive and selective amplification of methylated DNA sequences using helper-dependent chain reaction in combination with a methylation-dependent restriction enzymes

We have developed a novel technique for specific amplification of rare methylated DNA fragments in a high background of unmethylated sequences that avoids the need of bisulphite conversion. The methylation-dependent restriction enzyme GlaI is used to selectively cut methylated DNA. Then targeted fragments are tagged using specially designed ‘helper’ oligonucleotides that are also used to maintain selection in subsequent amplification cycles in a process called ‘helper-dependent chain reaction’. The process uses disabled primers called ‘drivers’ that can only prime on each cycle if the helpers recognize specific sequences within the target amplicon. In this way, selection for the sequence of interest is maintained throughout the amplification, preventing amplification of unwanted sequences. Here we show how the method can be applied to methylated Septin 9, a promising biomarker for early diagnosis of colorectal cancer. The GlaI digestion and subsequent amplification can all be done in a single tube. A detection sensitivity of 0.1% methylated DNA in a background of unmethylated DNA was achieved, which was similar to the well-established Heavy Methyl method that requires bisulphite-treated DNA.     

Nucleotide sequence analysis with polynucleotide kinase and nucleotide `mapping'' methods. 5′-Terminal sequence of deoxyribonucleic acid from bacteriophages λ and 424

The polynucleotide kinase reaction was used in analyses of complex mixtures of oligodeoxynucleotides which were fractionated by various two-dimensional nucleotide ;mapping' procedures. Parallel ionophoretic analyses on DEAE-cellulose paper, pH2, and AE-cellulose paper, pH3.5, of venom phosphodiesterase partial digests of 5'-terminally labelled oligonucleotides enabled the sequence of the nucleotides to be deduced uniquely. A ;diagonal ionophoresis' method has been used with mixtures of nucleotides. Application of these methods to 5'-terminally labelled DNA from bacteriophage lambda gave the terminal sequences pA-G-G-T-C-G and pG-G-G-C-G. Identical 5'-terminal sequences were found with DNA from bacteriophage 424.     

Prediction of the frequencies of restriction endonuclease recognition sequences using di- and mononucleotide frequencies

The calculation of probabilities of nucleotide sequences from the frequencies of dinucleotides is described. The dinucleotide and mononucleotide frequencies used can be obtained from nearest neighbor analysis or from databank sequences. If dinucleotide and mononucleotide frequencies from nearest neighbor analysis are used, probabilities for oligonucleotides can be calculated for genomes in which there is little or no sequence data. Within a given genome, a broad range of probabilities for hexanucleotide palindromes with the same base composition is predicted and shown (14).     

SNP-specific extraction of haplotype-resolved targeted genomic regions

The availability of genotyping platforms for comprehensive genetic analysis of complex traits has resulted in a plethora of studies reporting the association of specific single-nucleotide polymorphisms (SNPs) with common diseases or drug responses. However, detailed genetic analysis of these associated regions that would correlate particular polymorphisms to phenotypes has lagged. This is primarily due to the lack of technologies that provide additional sequence information about genomic regions surrounding specific SNPs, preferably in haploid form. Enrichment methods for resequencing should have the specificity to provide DNA linked to SNPs of interest with sufficient quality to be used in a cost-effective and high-throughput manner. We describe a simple, automated method of targeting specific sequences of genomic DNA that can directly be used in downstream applications. The method isolates haploid chromosomal regions flanking targeted SNPs by hybridizing and enzymatically elongating oligonucleotides with biotinylated nucleotides based on their selective binding to unique sequence elements that differentiate one allele from any other differing sequence. The targeted genomic region is captured by streptavidin-coated magnetic particles and analyzed by standard genotyping, sequencing or microarray analysis. We applied this technology to determine contiguous molecular haplotypes across a ~150 kb genomic region of the major histocompatibility complex.     

PCR fragmentation of DNA

A method has been developed to prepare random DNA fragments using PCR. First, two cycles are carried out at 16 degrees C with the Klenow's fragment and oligonucleotides (random primers) with random 3'-sequences and the 5'-constant part containing the site for cloning with the site-specific endonuclease. The random primers can link to any DNA site, and random DNA fragments are formed during DNA synthesis. During the second cycle, after denaturation of the DNA and addition of the Klenow's fragment, the random primers can link to newly synthesized DNA strands, and after DNA synthesis single-stranded DNA fragments are produced which have a constant primer sequence at the 5'-end and a complementary to it sequence at the 3'-end. During the third cycle, the constant primer is added and double-stranded fragments with the constant primer sequences at both ends are formed during DNA synthesis. Incubation for 1 h at 37 degrees C degrades the oligonucleotides used at the first stage due to endonuclease activity of the Klenow's fragment. Then routine PCR amplification is carried out using the constant primer. This method is more advantageous than hydrodynamic methods of DNA fragmentation widely used for "shotgun" cloning.     

Enzymatic assignment of diastereomeric purity of stereodefined phosphorothioate oligonucleotides.

Enzymatic hydrolysis of stereoregular oligodeoxyribonucleoside phosphorothioates (PS-oligos) synthesized via the oxathiaphospholane method has been used for assignment of their diastereomeric purity. For this purpose, two well-known enzymes of established diastereoselectivity, nuclease P1 and snake venom phosphodiesterase (svPDE) have been used. However, because of some disadvantageous properties of svPDE, a search for other [Rp]-specific endonucleases was undertaken. Extracellular bacterial endonuclease isolated from Serratia marcescens accepts PS-oligos as substrates and hydrolyzes phosphorothioate bonds of the [Rp] configuration, whereas internucleotide [Sp]-phosphorothioates are resistant to its action. Cleavage experiments carried out with the use of unmodified and phosphorothioate oligonucleotides of different sequences demonstrate that the Serratia nuclease is more selective in recognition and hydrolysis of oligodeoxyribonucleotides than previously reported. The substrate specificity exhibited by the enzyme is influenced not only by the nucleotide sequence at the cleavage site but also by the length and base sequence of flanking sequences. The Serratia nuclease can be useful for analysis of diastereomeric purity of stereodefined phosphorothioate oligonucleotides, but because of its sequence preferences, the use of this enzyme in conjunction with svPDE is more reliable.     

Linguistic measure of taxonomic and functional relatedness of nucleotide sequences

The frequencies of "words", oligonucleotides within nucleotide sequences, reflect the genetic information contained in the sequence "texts". Nucleotide sequences are characteristically represented by their contrast word vocabularies. Comparison of the sequences by correlating their contrast vocabularies is shown to reflect well the relatedness (unrelatedness) between the sequences. A single value, the linguistic similarity between the sequences, is suggested as a measure of sequence relatedness. Sequences as short as 1000 bases can be characterized and quantitatively related to other sequences by this technique. The linguistic sequence similarity value is used for analysis of taxonomically and functionally diverse nucleotide sequences. The similarity value is shown to be very sensitive to the relatedness of the source species, thus providing a convenient tool for taxonomic classification of species by their sequence vocabularies. Functionally diverse sequences appear distinct by their linguistic similarity values. This can be a basis for a quick screening technique for functional characterization of the sequences and for mapping functionally distinct regions in long sequences.     

Interactions of Escherichia coli RNA with bacteriophage MS2 coat protein: genomic SELEX

Genomic SELEX is a method for studying the network of nucleic acid–protein interactions within any organism. Here we report the discovery of several interesting and potentially biologically important interactions using genomic SELEX. We have found that bacteriophage MS2 coat protein binds several Escherichia coli mRNA fragments more tightly than it binds the natural, well-studied, phage mRNA site. MS2 coat protein binds mRNA fragments from rffG (involved in formation of lipopolysaccharide in the bacterial outer membrane), ebgR (lactose utilization repressor), as well as from several other genes. Genomic SELEX may yield experimentally induced artifacts, such as molecules in which the fixed sequences participate in binding. We describe several methods (annealing of oligonucleotides complementary to fixed sequences or switching fixed sequences) to eliminate some, or almost all, of these artifacts. Such methods may be useful tools for both randomized sequence SELEX and genomic SELEX.     

Physical map of the Kirsten sarcoma virus genome as determined by fingerprinting RNase T1-resistant oligonucleotides.

From analysis of the large RNase T1-resistant oligonucleotides of Kirsten sarcoma virus (Ki-SV), a physical map of the virus genome was deduced. Kirsten murine leukemia virus (Ki-MuLV) sequences were detected in T1 oligonucleotides closest to the 3' end of the viral RNA and extended approximately 1,000 nucleotides into the genome. The rat genetic sequences started at this point and extended all the way to the very 5' end of the RNA molecules, where a small stretch of Ki-MuLV sequence was detected. By comparison of the fingerprints of Ki-SV RNA and the RNA of the endogenous rat src genetic sequences, it was found that more than 50% of the T1 oligonucleotides were similar between Ki-SV and the endogenous rat src RNA, suggesting an identical primary nucleotide sequence in over 50% of the viral genomes. The results indicate that Ki-SV arose by recombination between the 5' and 3' ends of Ki-MuLV and a large portion of the homologous sequences of the endogenous rat src RNA.