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.     

PROBEmer: A web-based software tool for selecting optimal DNA oligos

PROBEmer (http://probemer.cs.loyola.edu) is a web-based software tool that enables a researcher to select optimal oligos for PCR applications and multiplex detection platforms including oligonucleotide microarrays and bead-based arrays. Given two groups of nucleic-acid sequences, a target group and a non-target group, the software identifies oligo sequences that occur in members of the target group, but not in the non-target group. To help predict potential cross hybridization, PROBEmer computes all near neighbors in the non-target group and displays their alignments. The software has been used to obtain genus-specific prokaryotic probes based on the 16S rRNA gene, gene-specific probes for expression analyses and PCR primers. In this paper, we describe how to use PROBEmer, the computational methods it employs, and experimental results for oligos identified by this software tool.     

Mprobe 2.0: computer-aided probe design for oligonucleotide microarray

DNA chips have proven to be effective tools in detecting gene expression levels. Compared with DNA chips using complementary DNA as probes, oligonucleotide microarrays using oligonucleotides as probes have attracted great attention because of their well known advantages. The design of gene-specific probes for each target is essential to the development of oligonucleotide microarrays. We have previously reported the development of a probe design software termed Mprobe 1.0. Here, we present a new version of this software, termed Mprobe 2.0. Several new features are included in Mprobe 2.0. Firstly, a paradox-based sequence database management system has been developed and integrated into the software, which consequently allows interoperability with sequences in GenBank, EMBL, and FASTA formats. Secondly, in contrast to setting a fixed threshold for the secondary structure of probes in Mprobe 1.0 and other related software, Mprobe 2.0 employs a different method. After parameters such as GC type, probe melting temperature and GC contents have been evaluated, candidate probes are sorted by the free energy from high to low value, followed by specificity analysis. Thirdly, Mprobe 2.0 provides users with substantial parameter options in the visual mode. Mprobe 2.0 possesses an easier interface for users to manage sequences annotated in different formats and design the optimal probes for oligonucleotide microarrays and other applications. AVAILABILITY: The program is free for non-commercial users and can be downloaded from the web page http://www.biosun.org.cn/mprobe/ CONTACT: Wuju Li (wujuli@yahoo.com or liwj@nic.bmi.ac.cn).     

Mathematical tools to optimize the design of oligonucleotide probes and primers

The identification and quantification of specific organisms in mixed microbial communities often relies on the ability to design oligonucleotide probes and primers with high specificity and sensitivity. The design of these oligonucleotides (or “oligos” for short) shares many of the same principles in spite of their widely divergent applications. Three common molecular biology technologies that require oligonucleotide design are polymerase chain reaction (PCR), fluorescence in situ hybridization (FISH), and DNA microarrays. This article reviews techniques and software available for the design and optimization of oligos with the goal of targeting a specific group of organisms within mixed microbial communities. Strategies for enhancing specificity without compromising sensitivity are described, as well as design tools well suited for this purpose.     

Using DNA microarrays to study gene expression in closely related species

MOTIVATION: Comparisons of gene expression levels within and between species have become a central tool in the study of the genetic basis for phenotypic variation, as well as in the study of the evolution of gene regulation. DNA microarrays are a key technology that enables these studies. Currently, however, microarrays are only available for a small number of species. Thus, in order to study gene expression levels in species for which microarrays are not available, researchers face three sets of choices: (i) use a microarray designed for another species, but only compare gene expression levels within species, (ii) construct a new microarray for every species whose gene expression profiles will be compared or (iii) build a multi-species microarray with probes from each species of interest. Here, we use data collected using a multi-primate cDNA array to evaluate the reliability of each approach. RESULTS: We find that, for inter-species comparisons, estimates of expression differences based on multi-species microarrays are more accurate than those based on multiple species-specific arrays. We also demonstrate that within-species expression differences can be estimated using a microarray for a closely related species, without discernible loss of information. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Genome-directed primers for selective labeling of bacterial transcripts for DNA microarray analysis

DNA microarrays have the ability to analyze the expression of thousands of the same set of genes under at least two different experimental conditions. However, DNA microarrays require substantial amounts of RNA to generate the probes, especially when bacterial RNA is used for hybridization (50 microg of bacterial total RNA contains approximately 2 microg of mRNA). We have developed a computer-based algorithm for prediction of the minimal number of primers to specifically anneal to all genes in a given genome. The algorithm predicts, for example, that 37 oligonucleotides should prime all genes in the Mycobacterium tuberculosis genome. We tested the usefulness of the genome-directed primers (GDPs) in comparison to random primers for gene expression profiling using DNA microarrays. Both types of primers were used to generate fluorescent-labeled probes and to hybridize to an array of 960 mycobacterial genes. Compared to random-primer probes, the GDP probes were more sensitive and more specific, especially when mammalian RNA samples were spiked with mycobacterial RNA. The GDPs were used for gene expression profiling of mycobacterial cultures grown to early log or stationary growth phases. This approach could be useful for accurate genome-wide expression analysis, especially for in vivo gene expression profiling, as well as directed amplification of sequenced genomes.     

On-chip hybridization kinetics for optimization of gene expression experiments

DNA microarray technology is a powerful tool for getting an overview of gene expression in biological samples. Although the successful use of microarray-based expression analysis was demonstrated in a number of applications, the main problem with this approach is the fact that expression levels deduced from hybridization experiments do not necessarily correlate with RNA concentrations. Moreover oligonucleotide probes corresponding to the same gene can give different hybridization signals. Apart from cross-hybridizations and differential splicing, this could be due to secondary structures of probes or targets. In addition, for low-copy genes, hybridization equilibrium may be reached after hybridization times much longer than the one commonly used (overnight, i.e., 15 h). Thus, hybridization signals could depend on kinetic properties of the probe, which may vary between different oligonucleotide probes immobilized on the same microarray. To validate this hypothesis, on-chip hybridization kinetics and duplex thermostability analysis were performed using oligonucleotide microarrays containing 50-mer probes corresponding to 10 mouse genes. We demonstrate that differences in hybridization kinetics between the probes exist and can influence the interpretation of expression data. In addition, we show that using on-chip hybridization kinetics, quantification of targets is feasible using calibration curves.     

A general framework for designing and validating oligomer-based DNA microarrays and its application to Clostridium acetobutylicum

While DNA microarray analysis is widely accepted as an essential tool for modern biology, its use still eludes many researchers for several reasons, especially when microarrays are not commercially available. In that case, the design, construction, and use of microarrays for a sequenced organism constitute substantial, time-consuming, and expensive tasks. Recently, it has become possible to construct custom microarrays using industrial manufacturing processes, which offer several advantages, including speed of manufacturing, quality control, no up-front setup costs, and need-based microarray ordering. Here, we describe a strategy for designing and validating DNA microarrays manufactured using a commercial process. The 22K microarrays for the solvent producer Clostridium acetobutylicum ATCC 824 are based on in situ-synthesized 60-mers employing the Agilent technology. The strategy involves designing a large library of possible oligomer probes for each target (i.e., gene or DNA sequence) and experimentally testing and selecting the best probes for each target. The degenerate C. acetobutylicum strain M5 lacking the pSOL1 megaplasmid (with 178 annotated open reading frames [genes]) was used to estimate the level of probe cross-hybridization in the new microarrays and to establish the minimum intensity for a gene to be considered expressed. Results obtained using this microarray design were consistent with previously reported results from spotted cDNA-based microarrays. The proposed strategy is applicable to any sequenced organism.     

A method for cross-species gene expression analysis with high-density oligonucleotide arrays

DNA microarrays have been widely used in gene expression analysis of biological processes. Due to a lack of sequence information, the applications have been largely restricted to humans and a few model organisms. Presented within this study are results of the cross-species hybridization with Affymetrix human high-density oligonucleotide arrays or GeneChip® using distantly related mammalian species; cattle, pig and dog. Based on the unique feature of the Affymetrix GeneChip® where every gene is represented by multiple probes, we hypothesized that sequence conservation within mammals is high enough to generate sufficient signals from some of the probes for expression analysis. We demonstrated that while overall hybridization signals are low for cross-species hybridization, a few probes of most genes still generated signals equivalent to the same-species hybridization. By masking the poorly hybridized probes electronically, the remaining probes provided reliable data for gene expression analysis. We developed an algorithm to select the reliable probes for analysis utilizing the match/mismatch feature of GeneChip®. When comparing gene expression between two tissues using the selected probes, we found a linear correlation between the cross-species and same-species hybridization. In addition, we validated cross-species hybridization results by quantitative PCR using randomly selected genes. The method shown herein could be applied to both plant and animal research.     

Comparative genomics and experimental promoter analysis reveal functional liver-specific elements in mammalian hepatic lipase genes

Background

Array-based comparative genomic hybridization (aCGH) is a high-throughput method for measuring genome-wide DNA copy number changes. Current aCGH methods have limited resolution, sensitivity and reproducibility. Microarrays for aCGH are available only for a few organisms and combination of aCGH data with expression data is cumbersome.

Results

We present a novel method of using commercial oligonucleotide expression microarrays for aCGH, enabling DNA copy number measurements and expression profiles to be combined using the same platform. This method yields aCGH data from genomic DNA without complexity reduction at a median resolution of approximately 17,500 base pairs. Due to the well-defined nature of oligonucleotide probes, DNA amplification and deletion can be defined at the level of individual genes and can easily be combined with gene expression data.

Conclusion

A novel method of gene resolution analysis of copy number variation (graCNV) yields high-resolution maps of DNA copy number changes and is applicable to a broad range of organisms for which commercial oligonucleotide expression microarrays are available. Due to the standardization of oligonucleotide microarrays, graCNV results can reliably be compared between laboratories and can easily be combined with gene expression data using the same platform.     

Assessment of cry1 gene contents of Bacillus thuringiensis strains by use of DNA microarrays

A single Bacillus thuringiensis strain can harbor numerous different insecticidal crystal protein (cry) genes from 46 known classes or primary ranks. The cry1 primary rank is the best known and contains the highest number of cry genes which currently totals over 130. We have designed an oligonucleotide-based DNA microarray (cryArray) to test the feasibility of using microarrays to identify the cry gene content of B. thuringiensis strains. Specific 50-mer oligonucleotide probes representing the cry1 primary and tertiary ranks were designed based on multiple cry gene sequence alignments. To minimize false-positive results, a consentaneous approach was adopted in which multiple probes against a specific gene must unanimously produce positive hybridization signals to confirm the presence of a particular gene. In order to validate the cryArray, several well-characterized B. thuringiensis strains including isolates from a Mexican strain collection were tested. With few exceptions, our probes performed in agreement with known or PCR-validated results. In one case, hybridization of primary- but not tertiary-ranked cry1I probes indicated the presence of a novel cry1I gene. Amplification and partial sequencing of the cry1I gene in strains IB360 and IB429 revealed the presence of a cry1Ia gene variant. Since a single microarray hybridization can replace hundreds of individual PCRs, DNA microarrays should become an excellent tool for the fast screening of new B. thuringiensis isolates presenting interesting insecticidal activity.     

Cutting-edge technology. I. Global gene expression profiling using DNA microarrays

Having the complete human genomic sequence poses a new challenge: to use genomic structural information to display and analyze biological processes on a genome-wide scale to assign gene function. DNA microarrays are a miniaturized, ordered arrangement of nucleic acid fragments from individual genes located at defined positions on a solid support, enabling the analysis of thousands of genes in parallel by specific hybridization. This review describes technical aspects, discusses relevant applications, and suggests factors affecting the use of this technology and how it fits in the grand scheme of meeting the needs of the postgenomic era.     

A General Framework for Designing and Validating Oligomer-Based DNA Microarrays and Its Application to Clostridium acetobutylicum

While DNA microarray analysis is widely accepted as an essential tool for modern biology, its use still eludes many researchers for several reasons, especially when microarrays are not commercially available. In that case, the design, construction, and use of microarrays for a sequenced organism constitute substantial, time-consuming, and expensive tasks. Recently, it has become possible to construct custom microarrays using industrial manufacturing processes, which offer several advantages, including speed of manufacturing, quality control, no up-front setup costs, and need-based microarray ordering. Here, we describe a strategy for designing and validating DNA microarrays manufactured using a commercial process. The 22K microarrays for the solvent producer Clostridium acetobutylicum ATCC 824 are based on in situ-synthesized 60-mers employing the Agilent technology. The strategy involves designing a large library of possible oligomer probes for each target (i.e., gene or DNA sequence) and experimentally testing and selecting the best probes for each target. The degenerate C. acetobutylicum strain M5 lacking the pSOL1 megaplasmid (with 178 annotated open reading frames [genes]) was used to estimate the level of probe cross-hybridization in the new microarrays and to establish the minimum intensity for a gene to be considered expressed. Results obtained using this microarray design were consistent with previously reported results from spotted cDNA-based microarrays. The proposed strategy is applicable to any sequenced organism.