Selection of oligonucleotide probes for protein coding sequences

MOTIVATION: Large arrays of oligonucleotide probes have become popular tools for analyzing RNA expression. However to date most oligo collections contain poorly validated sequences or are biased toward untranslated regions (UTRs). Here we present a strategy for picking oligos for microarrays that focus on a design universe consisting exclusively of protein coding regions. We describe the constraints in oligo design that are imposed by this strategy, as well as a software tool that allows the strategy to be applied broadly. RESULT: In this work we sequentially apply a variety of simple filters to candidate sequences for oligo probes. The primary filter is a rejection of probes that contain contiguous identity with any other sequence in the sample universe that exceeds a pre-established threshold length. We find that rejection of oligos that contain 15 bases of perfect match with other sequences in the design universe is a feasible strategy for oligo selection for probe arrays designed to interrogate mammalian RNA populations. Filters to remove sequences with low complexity and predicted poor probe accessibility narrow the candidate probe space only slightly. Rejection based on global sequence alignment is performed as a secondary, rather than primary, test, leading to an algorithm that is computationally efficient. Splice isoforms pose unique challenges and we find that isoform prevalence will for the most part have to be determined by analysis of the patterns of hybridization of partially redundant oligonucleotides. AVAILABILITY: The oligo design program OligoPicker and its source code are freely available at our website.     

Thermodynamic properties of DNA sequences: characteristic values for the human genome

MOTIVATION: Central to many molecular biology techniques as ubiquitous as PCR and Southern blotting is the design of oligonucleotide (oligo) probes and/or primers possessing specific thermodynamic properties. Here, we use validated theoretical methods to generate distributions of predicted thermodynamic properties for DNA oligos of various lengths. These distributions facilitate immediate appreciation of typical thermodynamic values for oligos of various lengths. RESULTS: Distributions of melting temperature (Tm), free energy (DeltaG(T)o), and fraction hybridized or fraction bound (Fb), are presented for oligos of length 10-50 bases sampled from the human genome. The effects of changing temperature, oligo and salt concentrations, constraining G+C content, and introducing mismatches are exemplified. Our results provide the first survey of typical and limiting thermodynamic values evaluated on a genomic scale. Described numbers comprise useful 'rules of thumb' that are applicable to most technologies dependent upon DNA oligo design.     

Oligo replication advantage driven by GC content and Gibbs free energy

Large scale DNA oligo pools are emerging as a novel material in a variety of advanced applications. However, GC content and length cause significant bias in amplification of oligos. We systematically explored the amplification of one oligo pool comprising of over ten thousand distinct strands with moderate GC content in the range of 35–65%. Uniqual amplification of oligos result to the increased Gini index of the oligo distribution while a few oligos greatly increased their proportion after 60 cycles of PCR. However, the significantly enriched oligos all have relatively high GC content. Further thermodynamic analysis demonstrated that a high value of both GC content and Gibbs free energy could improve the replication of specific oligos during biased amplification. Therefore, this double-G (GC content and Gibbs free energy) driven replication advantage can be used as a guiding principle for the sequence design for a variety of applications, particularly for data storage.

     

Efficient selection of unique and popular oligos for large EST databases

MOTIVATION: Expressed sequence tag (EST) databases have grown exponentially in recent years and now represent the largest collection of genetic sequences. An important application of these databases is that they contain information useful for the design of gene-specific oligonucleotides (or simply, oligos) that can be used in PCR primer design, microarray experiments and genomic library screening. RESULTS: In this paper, we study two complementary problems concerning the selection of short oligos, e.g. 20-50 bases, from a large database of tens of thousands of ESTs: (i) selection of oligos each of which appears (exactly) in one unigene but does not appear (exactly or approximately) in any other unigene and (ii) selection of oligos that appear (exactly or approximately) in many unigenes. The first problem is called the unique oligo problem and has applications in PCR primer and microarray probe designs, and library screening for gene-rich clones. The second is called the popular oligo problem and is also useful in screening genomic libraries. We present an efficient algorithm to identify all unique oligos in the unigenes and an efficient heuristic algorithm to enumerate the most popular oligos. By taking into account the distribution of the frequencies of the words in the unigene database, the algorithms have been engineered carefully to achieve remarkable running times on regular PCs. Each of the algorithms takes only a couple of hours (on a 1.2 GHz CPU, 1 GB RAM machine) to run on a dataset 28 Mb of barley unigenes from the HarvEST database. We present simulation results on the synthetic data and a preliminary analysis of the barley unigene database. AVAILABILITY: Available on request from the authors.     

Rapid nuclear accumulation of injected oligodeoxyribonucleotides.

The intracellular transport and fate of nucleic acids is poorly understood. To study this process, we injected fluorescent oligodeoxyribonucleotides (oligos) into the cytoplasm of CV-1 epithelial cells and primary human fibroblasts. Rapid nuclear accumulation was found with the phosphodiester (PD), phosphorothioate (PT), and methylphosphonate (MP) forms of a 28-mer oligo complimentary to the rev mRNA of the human immunodeficiency virus type 1. Migration of the oligos in the cytoplasm was slower than diffusion of a coinjected dextran, but the oligos freely diffused into the nucleus. Nuclear incorporation was temperature but not energy dependent. The intranuclear distribution of the oligos was influenced by the chemistry of internucleoside linkages. The PD oligos and, to a lesser extent, the PT oligos colocalized with small nuclear ribonucleoproteins (snRNPs), whereas the MP oligos colocalized with concentrated regions of genomic DNA. These data have important implications for our understanding of the transport and accumulation of exogenous nucleic acids in mammalian nuclei, and the assay described could potentially be used for testing the efficacy of oligos designed as therapeutic agents.     

Oligonucleotide recombination in Gram-negative bacteria

This report describes several key aspects of a novel form of RecA-independent homologous recombination. We found that synthetic single-stranded DNA oligonucleotides (oligos) introduced into bacteria by transformation can site-specifically recombine with bacterial chromosomes in the absence of any additional phage-encoded functions. Oligo recombination was tested in four genera of Gram-negative bacteria and in all cases evidence for recombination was apparent. The experiments presented here were designed with an eye towards learning to use oligo recombination in order to bootstrap identification and development of phage-encoded recombination systems for recombineering in a wide range of bacteria. The results show that oligo concentration and sequence have the greatest influence on recombination frequency, while oligo length was less important. Apart from the utility of oligo recombination, these findings also provide insights regarding the details of recombination mediated by phage-encoded functions. Establishing that oligos can recombine with bacterial genomes provides a link to similar observations of oligo recombination in archaea and eukaryotes suggesting the possibility that this process is evolutionary conserved.     

Bacterial DNA polymerases participate in oligonucleotide recombination

Synthetic single‐strand oligonucleotides (oligos) with homology to genomic DNA have proved to be highly effective for constructing designed mutations in targeted genomes, a process referred to as recombineering. The cellular functions important for this type of homologous recombination have yet to be determined. Towards this end, we have identified Escherichia coli functions that process the recombining oligo and affect bacteriophage λ Red‐mediated oligo recombination. To determine the nature of oligo processing during recombination, each oligo contained multiple nucleotide changes: a single base change allowing recombinant selection, and silent changes serving as genetic markers to determine the extent of oligo processing during the recombination. Such oligos were often not incorporated into the host chromosome intact; many were partially degraded in the process of recombination. The position and number of these silent nucleotide changes within the oligo strongly affect both oligo processing and recombination frequency. Exonucleases, especially those associated with DNA Polymerases I and III, affect inheritance of the silent nucleotide changes in the oligos. We demonstrate for the first time that the major DNA polymerases (Pol I and Pol III) and DNA ligase are directly involved with oligo recombination.     

An experimental and theoretical study of the inhibition of Escherichia coli lac operon gene expression by antigene oligonucleotides

Previously, we have developed a genetically structured mathematical model to describe the inhibition of Escherichia coli lac operon gene expression by antigene oligos. Our model predicted that antigene oligos targeted to the operator region of the lac operon would have a significant inhibitory effect on beta-galactosidase production. In this investigation, the E. coli lac operon gene expression in the presence of antigene oligos was studied experimentally. A 21-mer oligo, which was designed to form a triplex with the operator, was found to be able to specifically inhibit beta-galactosidase production in a dose-dependent manner. In contrast to the 21-mer triplex-forming oligonucleotide (TFO), several control oligos showed no inhibitory effect. The ineffectiveness of the various control oligos, along with the fact that the 21-mer oligo has no homology sequence with lacZYA, and no mRNA is transcribed from the operator, suggests that the 21-mer oligo inhibits target gene expression by an antigene mechanism. To simulate the kinetics of lac operon gene expression in the presence of antigene oligos, a genetically structured kinetic model, which includes transport of oligo into the cell, growth of bacteria cells, and lac operon gene expression, was developed. Predictions of the kinetic model fit the experimental data quite well after adjustment of the value of the oligonucleotide transport rate constant (9.0 x 10(-)(3) min(-)(1)) and oligo binding affinity constant (1.05 x 10(6) M(-)(1)). Our values for these two adjusted parameters are in the range of reported literature values.     

Fast large scale oligonucleotide selection using the longest common factor approach

We present a fast method that selects oligonucleotide probes (such as DNA 25-mers) for microarray experiments on a truly large scale. For example, reliable oligos for human genes can be found within four days, a speedup of one to two orders of magnitude compared to previous approaches. This speed is attained by using the longest common substring as a specificity measure for candidate oligos. We present a space- and time-efficient algorithm, based on a suffix array with additional information, to compute matching statistics (lengths of longest matches) between all candidate oligos and all remaining sequences. With the matching statistics available, we show how to incorporate constraints such as oligo length, melting temperature, and self-complementarity into the selection process at a postprocessing stage. As a result, we can now design custom oligos for any sequenced genome, just as the technology for on-site chip synthesis is becoming increasingly mature.     

Patterns of intracellular compartmentalization, trafficking and acidification of 5'-fluorescein labeled phosphodiester and phosphorothioate oligodeoxynucleotides in HL60 cells.

We have examined the intracellular compartmentalization and trafficking of fluorescein labeled (F) phosphodiester (PO) and phosphorothioate (PS) oligodeoxynucleotides (oligos) in HL60 cells. A series of F-oligos (PO and PS) were incubated for 6 hrs. with HL60 cells and the mean intracellular fluorescence determined by flow cytometry. The F signal was normalized by the addition of the ionophore monensin. An increase in signal intensity following addition of monensin indicated that the oligo was resident in an acidic intracellular environment. F-PS, but not F-PO oligos were found to reside in an acidic environment. An exception was a PO homopolymer of 15 cytidine bases (FOdC15) which was acidified. Using two different methods, the average resident intracellular pH of F-PS oligos and F-OdC15 was shown to be approximately 1 pH unit lower than that of F-PO oligos. Acidification of F-PS oligos could be blocked by the antibiotic bafilomycin, indicating that acidification was occurring in endosomes or vacuoles. F-PO and F-PS oligos were effluxed from HL60 cells from two intracellular compartments. However, approximately 60% of internalized F-PO oligo resided in a 'shallow' compartment that was turned over rapidly (t1/2 = 5-10 min.) whereas only 20% of F-PS oligo resided in this compartment. Conversely, approximately 80% of the internalized F-PS oligo but only 40% of F-PO oligo resided in a 'deep' compartment that turned over with t1/2 = 2-5 hrs. This report is the first quantitative demonstration that PO and PS oligos, and PO oligos of different sequences are trafficked differently by HL60 cells.     

Specificity of the nick-closing activity of bacteriophage T4 DNA ligase

Bacteriophage T4 DNA ligase effectively joins two adjacent, short synthetic oligodeoxyribonucleotides (oligos), as guided by complementary oligo, plasmid and genomic DNA templates. When a single bp mismatch exists at either side of the ligation junction, the efficiency of the enzyme to ligate the two oligos decreases. Mismatch ligation is approximately five-fold greater if the mismatch occurs at the 3' side rather than at the 5' side of the junction. During mismatch ligation the 5' adenylate of the 3' oligo accumulates in the reaction. The level of the adenylate formation correlates closely with the level of the mismatch ligation. Both mismatch ligation and adenylate formation are suppressed at elevated temperatures and in the presence of 200 mM NaCl or 2-5 mM spermidine. The apparent Km for the oligo template in the absence of salt is 0.05 microM, whereas the Km increases to 0.2 microM in the presence of 200 mM of NaCl. In this report, we demonstrate these properties of T4 DNA ligase for oligo pairs complementary to the beta-globin gene at the sequence surrounding the single bp mutation responsible for sickle-cell anemia. Because of the highly specific nature of the nick-closing reaction, ligation of short oligos with DNA ligase can be used to distinguish two DNA templates differing by a single nucleotide.     

Transformation with oligonucleotides creating clustered changes in the yeast genome

We have studied single-strand oligonucleotide (oligo) transformation of yeast by using 40-nt long oligos that create multiple base changes to the yeast genome spread throughout the length of the oligos, making it possible to measure the portions of an oligo that are incorporated during transformation. Although the transformation process is greatly inhibited by DNA mismatch repair (MMR), the pattern of incorporation is essentially the same in the presence or absence of MMR, whether the oligo anneals to the leading or lagging strand of DNA replication, or whether phosphorothioate linkages are used at either end. A central core of approximately 15 nt is incorporated with a frequency of >90%; the ends are incorporated with a lower frequency, and loss of the two ends appears to be by different mechanisms. Bases that are 5-10 nt from the 5' end are generally lost with a frequency of >95%, likely through a process involving flap excision. On the 3' end, bases 5-10 nt from the 3' end are lost about 1/3 of the time. These results indicate that oligos can be used to create multiple simultaneous changes to the yeast genome, even in the presence of MMR.     

Sensitivity of 70-mer oligonucleotides and cDNAs for microarray analysis of gene expression in Arabidopsis and its related species

Synthetic oligonucleotides (oligos) represent an attractive alternative to cDNA amplicons for spotted microarray analysis in a number of model organisms, including Arabidopsis, C. elegans, Drosophila, human, mouse and yeast. However, little is known about the relative effectiveness of 60-70-mer oligos and cDNAs for detecting gene expression changes. Using 192 pairs of Arabidopsis thaliana cDNAs and corresponding 70-mer oligos, we performed three sets of dye-swap experiments and used analysis of variance (anova) to compare sources of variation and sensitivities for detecting gene expression changes in A. thaliana, A. arenosa and Brassica oleracea. Our major findings were: (1) variation among different RNA preparations from the same tissue was small, but large variation among dye-labellings and slides indicates the need to replicate these factors; (2) sources of variation were similar for experiments with all three species, suggesting these feature types are effective for analysing gene expression in related species; (3) oligo and cDNA features had similar sensitivities for detecting expression changes and they identified a common subset of significant genes, but results from quantitative RT-PCR did not support the use of one over the other. These findings indicate that spotted oligos are at least as effective as cDNAs for microarray analyses of gene expression. We are using oligos designed from approximately 26,000 annotated genes of A. thaliana to study gene expression changes in Arabidopsis and Brassica polyploids.     

Pseudo-cyclic oligonucleotides: in vitro and in vivo properties

We have designed and studied antisense oligodeoxynucleotides (oligonucleotides; oligos) which we call ‘pseudo-cyclic oligonucleotides’ (PCOs). PCOs contain two oligonucleotide segments attached through their 3′-3′- or 5′-5′-ends. One of the segments of the PCO is an antisense oligo complementary to a target mRNA, and the other is a short protective oligo that is 5–8 nucleotides long and complementary to the 3′- or 5′-end of the antisense oligo. As a result of complementarity between the antisense and protective oligo segments, PCOs form intramolecular pseudo-cyclic structures in the absence of the target RNA. The antisense oligo segment of PCOs used for the studies described here is complementary to an 18-nucleotide-long site on the mRNA of the protein kinase A regulatory subunit RI (PKA-RI). Thermal melting studies of PCOs in the absence and presence of the complementary RNA suggest that the pseudo-cyclic structures formed in the absence of the target RNA dissociate, bind to the target RNA, and form heteroduplexes. The results of RNase H cleavage assays suggest that PCOs bind to complementary RNA and activate RNase H in a manner similar to that of an 18-mer conventional antisense PS-oligo. In snake venom (a 3′-exonuclease) or spleen (a 5′-exonuclease) phosphodiesterase digestion studies, PCOs are more stable than conventional antisense oligos because of the presence of 3′-3′- or 5′-5′-linkages and the formation of intramolecular pseudo-cyclic structures. PCOs with a phosphorothioate antisense oligo segment inhibited cell growth of MDA-MB-468 and GEO cancer cell lines similar to that of the conventional antisense PS-oligo, suggesting efficient cellular uptake and target binding. The nuclease stability studies in mice suggest that PCOs have higher in vivo stability than antisense PS-oligos. The studies in mice showed similar pharmacokinetic and tissue distribution profiles for PCOs to those of antisense PS-oligos in general, but rapid elimination from selected tissues.     

A more efficient and specific strategy in the ablation of mRNA in Xenopus laevis using mixtures of antisense oligos.

Previously, antisense oligodeoxyribonucleotides (oligos) have been used to ablate specific mRNAs from the maternal RNA pool of Xenopus laevis oocytes. However, this strategy is limited by the dose of oligo which can be used and the fact that 100% cleavage of the target RNA is rare. Further, non-specific cleavage of other RNAs can also occur. We demonstrate that the use of several oligos against the histone H4 RNA results in a marked improvement in the efficiency of target degradation, due to synergistic action between oligos and the existence of RNA in at least two different secondary structures. We show, by using a set of overlapping oligos complementary to the entire H4 RNA, that the amount of oligo required for efficient target ablation is greatly lowered and non-specific effects are reduced.     

Flexible chromosome painting based on multiplex PCR of oligonucleotides and its application for comparative chromosome analyses in Cucumis

Chromosome painting is a powerful technique for chromosome and genome studies. We developed a flexible chromosome painting technique based on multiplex PCR of a synthetic oligonucleotide (oligo) library in cucumber (Cucumis sativus L., 2n = 14). Each oligo in the library was associated with a universal as well as nested specific primers for amplification, which allow the generation of different probes from the same oligo library. We were also able to generate double‐stranded labelled oligos, which produced much stronger signals than single‐stranded labelled oligos, by amplification using fluorophore‐conjugated primer pairs. Oligos covering cucumber chromosome 1 (Chr1) and chromosome 4 (Chr4) consisting of eight segments were synthesized in one library. Different oligo probes generated from the library painted the corresponding chromosomes/segments unambiguously, especially on pachytene chromosomes. This technique was then applied to study the homoeologous relationships among cucumber, C. hystrix and C. melo chromosomes based on cross‐species chromosome painting using Chr4 probes. We demonstrated that the probe was feasible to detect interspecies chromosome homoeologous relationships and chromosomal rearrangement events. Based on its advantages and great convenience, we anticipate that this flexible oligo‐painting technique has great potential for the studies of the structure, organization, and evolution of chromosomes in any species with a sequenced genome.     

Supersequences of masks for oligo-chips

On a very small surface, a chip, several thousands of oligonucleotides, can be synthesized using a mask technology. Given a set of oligos, the problems tackled here are: What is the minimum number of masks necessary to synthesize one copy of each oligo? How long will be the series of masks if each oligo is synthesized twice, such that the two copies are realized with two completely different series of masks? We establish that, for 20,000 oligos of 20 bases a single copy can be synthesized with around 67 masks and for two copies with less than the double.