GoArrays: highly dynamic and efficient microarray probe design

MOTIVATION: The use of oligonucleotide microarray technology requires a very detailed attention to the design of specific probes spotted on the solid phase. These problems are far from being commonplace since they refer to complex physicochemical constraints. Whereas there are more and more publicly available programs for microarray oligonucleotide design, most of them use the same algorithm or criteria to design oligos, with only little variation. RESULTS: We show that classical approaches used in oligo design software may be inefficient under certain experimental conditions, especially when dealing with complex target mixtures. Indeed, our biological model is a human obligate parasite, the microsporidia Encephalitozoon cuniculi. Targets that are extracted from biological samples are composed of a mixture of pathogen transcripts and host cell transcripts. We propose a new approach to design oligonucleotides which combines good specificity with a potentially high sensitivity. This approach is original in the biological point of view as well as in the algorithmic point of view. We also present an experimental validation of this new strategy by comparing results obtained with standard oligos and with our composite oligos. A specific E.cuniculi microarray will overcome the difficulty to discriminate the parasite mRNAs from the host cell mRNAs demonstrating the power of the microarray approach to elucidate the lifestyle of an intracellular pathogen using mix mRNAs

Physico-chemical foundations underpinning microarray and next-generation sequencing experiments

Hybridization of nucleic acids on solid surfaces is a key process involved in high-throughput technologies such as microarrays and, in some cases, next-generation sequencing (NGS). A physical understanding of the hybridization process helps to determine the accuracy of these technologies. The goal of a widespread research program is to develop reliable transformations between the raw signals reported by the technologies and individual molecular concentrations from an ensemble of nucleic acids. This research has inputs from many areas, from bioinformatics and biostatistics, to theoretical and experimental biochemistry and biophysics, to computer simulations. A group of leading researchers met in Ploen Germany in 2011 to discuss present knowledge and limitations of our physico-chemical understanding of high-throughput nucleic acid technologies. This meeting inspired us to write this summary, which provides an overview of the state-of-the-art approaches based on physico-chemical foundation to modeling of the nucleic acids hybridization process on solid surfaces. In addition, practical application of current knowledge is emphasized

Development of a DNA-based microarray for the detection of zoonotic pathogens in rodent species

The demand for diagnostic tools that allow simultaneous screening of samples for multiple pathogens is increasing because they overcome the limitations of other methods, which can only screen for a single or a few pathogens at a time. Microarrays offer the advantages of being capable to test a large number of samples simultaneously, screening for multiple pathogen types per sample and having comparable sensitivity to existing methods such as PCR. Array design is often considered the most important process in any microarray experiment and can be the deciding factor in the success of a study. There are currently no microarrays for simultaneous detection of rodent-borne pathogens. The aim of this report is to explicate the design, development and evaluation of a microarray platform for use as a screening tool that combines ease of use and rapid identification of a number of rodent-borne pathogens of zoonotic importance. Nucleic acid was amplified by multiplex biotinylation PCR prior to hybridisation onto microarrays. The array sensitivity was comparable to standard PCR, though less sensitive than real-time PCR. The array presented here is a prototype microarray identification system for zoonotic pathogens that can infect rodent species

Alternative splicing and bioinformatic analysis of human U12-type introns

U12-type introns exist, albeit rarely, in a variety of multicellular organisms. Splicing of U12 intron-containing precursor mRNAs takes place in the U12-type spliceosome that is distinct from the major U2-type spliceosome. Due to incompatibility of these two spliceosomes, alternative splicing involving a U12-type intron may give rise to a relatively complicated impact on gene expression. We studied alternative U12-type intron splicing in an attempt to gain more mechanistic insights. First, we characterized mutually exclusive exon selection of the human JNK2 gene, which involves an unusual intron possessing the U12-type 5′ splice site and the U2-type 3′ splice site. We demonstrated that the long and evolutionary conserved polypyrimidine tract of this hybrid intron provides important signals for inclusion of its downstream alternative exon. In addition, we examined the effects of single nucleotide polymorphisms in the human WDFY1 U12-type intron on pre-mRNA splicing. These results provide mechanistic implications on splice-site selection of U12-type intron splicing. We finally discuss the potential effects of splicing of a U12-type intron with genetic defects or within a set of genes encoding RNA processing factors on global gene expression

Application of Equilibrium Models of Solution Hybridization to Microarray Design and Analysis

Background: The probe percent bound value, calculated using multi-state equilibrium models of solution hybridization, is shown to be useful in understanding the hybridization behavior of microarray probes having 50 nucleotides, with and without mismatches. These longer oligonucleotides are in widespread use on microarrays, but there are few controlled studies of their interactions with mismatched targets compared to 25-mer based platforms. Principal Findings: 50-mer oligonucleotides with centrally placed single, double and triple mismatches were spotted on an array. Over a range of target concentrations it was possible to discriminate binding to perfect matches and mismatches, and the type of mismatch could be predicted accurately in the concentration midrange (100 pM to 200 pM) using solution hybridization modeling methods. These results have implications for microarray design, optimization and analysis methods. Conclusions: Our results highlight the importance of incorporating biophysical factors in both the design and the analysis of microarrays. Use of the probe ‘‘percent bound’ ’ value predicted by equilibrium models of hybridization is confirmed to be important for predicting and interpreting the behavior of long oligonucleotide arrays, as has been shown for shor

Model-based probe set optimization for high-performance microarrays

A major challenge in microarray design is the selection of highly specific oligonucleotide probes for all targeted genes of interest, while maintaining thermodynamic uniformity at the hybridization temperature. We introduce a novel microarray design framework (Thermodynamic Model-based Oligo Design Optimizer, TherMODO) that for the first time incorporates a number of advanced modelling features: (i) A model of position-dependent labelling effects that is quantitatively derived from experiment. (ii) Multi-state thermodynamic hybridization models of probe binding behaviour, including potential cross-hybridization reactions. (iii) A fast calibrated sequence-similarity-based heuristic for cross-hybridization prediction supporting large-scale designs. (iv) A novel compound score formulation for the integrated assessment of multiple probe design objectives. In contrast to a greedy search for probes meeting parameter thresholds, this approach permits an optimization at the probe set level and facilitates the selection of highly specific probe candidates while maintaining probe set uniformity. (v) Lastly, a flexible target grouping structure allows easy adaptation of the pipeline to a variety of microarray application scenarios. The algorithm and features are discussed and demonstrated on actual design runs. Source code is available on request

Multi-stringency wash of partially hybridized 60-mer probes reveals that the stringency along the probe decreases with distance from the microarray surface

Here, we describe a multi-parametric study of DNA hybridization to probes with 20–70% G + C content. Probes were designed towards 71 different sites/mutations in the phenylalanine hydroxylase gene. Seven probe lengths, three spacer lengths and six stringencies were systematically varied. The three spacer lengths were obtained by placing the gene-specific sequence in discrete steps along the 60-mer probes. The study was performed using Agilent 8 × 15 000 probes custom-made arrays and a home-built array washer providing different stringencies to each of the eight sub-arrays on the slides. Investigation of hybridization signals, specificity and dissociation curves indicated that probes close to the surface were influenced by an additional stringency provided by the microarray surface. Consistent with this, probes close to the surface required 4 × SSC, while probes placed away from the surface required 0.35 × SSC wash buffers in order to give accurate genotyping results. Multiple step dissociation was frequently observed for probes placed furthest away from surface, but not for probes placed proximal to the surface, which is consistent with the hypothesis that there is different stringency along the 60-mer. The results have impact on design of probes for genotyping, gene expression and comparative genome hybridization analysis

Investigating the global genomic diversity of Escherichia coli using a multi-genome DNA microarray platform with novel gene prediction strategies

<p>Abstract</p> <p>Background</p> <p>The gene content of a diverse group of 183 unique <it>Escherichia coli </it>and <it>Shigella </it>isolates was determined using the Affymetrix GeneChip^®<it>E. coli </it>Genome 2.0 Array, originally designed for transcriptome analysis, as a genotyping tool. The probe set design utilized by this array provided the opportunity to determine the gene content of each strain very accurately and reliably. This array constitutes 10,112 independent genes representing four individual <it>E. coli </it>genomes, therefore providing the ability to survey genes of several different pathogen types. The entire ECOR collection, 80 EHEC-like isolates, and a diverse set of isolates from our FDA strain repository were included in our analysis.</p> <p>Results</p> <p>From this study we were able to define sets of genes that correspond to, and therefore define, the EHEC pathogen type. Furthermore, our sampling of 63 unique strains of O157:H7 showed the ability of this array to discriminate between closely related strains. We found that individual strains of O157:H7 differed, on average, by 197 probe sets. Finally, we describe an analysis method that utilizes the power of the probe sets to determine accurately the presence/absence of each gene represented on this array.</p> <p>Conclusions</p> <p>These elements provide insights into understanding the microbial diversity that exists within extant <it>E. coli </it>populations. Moreover, these data demonstrate that this novel microarray-based analysis is a powerful tool in the field of molecular epidemiology and the newly emerging field of microbial forensics.</p

An evaluation of custom microarray applications: the oligonucleotide design challenge

The increase in feature resolution and the availability of multipack formats from microarray providers has opened the way to various custom genomic applications. However, oligonucleotide design and selection remains a bottleneck of the microarray workflow. Several tools are available to perform this work, and choosing the best one is not an easy task, nor are the choices obvious. Here we review the oligonucleotide design field to help users make their choice. We have first performed a comparative evaluation of the available solutions based on a set of criteria including: ease of installation, user-friendly access, the number of parameters and settings available. In a second step, we chose to submit two real cases to a selection of programs. Finally, we used a set of tests for the in silico benchmark of the oligo sets obtained from each type of software. We show that the design software must be selected according to the goal of the scientist, depending on factors such as the organism used, the number of probes required and their localization on the target sequence. The present work provides keys to the choice of the most relevant software, according to the various parameters we tested

Intraspecific variation in herbivore community composition and transcriptional profiles in field-grown Brassica oleracea cultivars

Intraspecific differences in plant defence traits are often correlated with variation in transcriptional profiles and can affect the composition of herbivore communities on field-grown plants. However, most studies on transcriptional profiling of plant–herbivore interactions have been carried out under controlled conditions in the laboratory or greenhouse and only a few examine intraspecific transcriptional variation. Here, intraspecific variation in herbivore community composition and transcriptional profiles between two Brassica oleracea cultivars grown in the field is addressed. Early in the season, no differences in community composition were found for naturally occurring herbivores, whereas cultivars differed greatly in abundance, species richness, and herbivore community later in the season. Genome-wide transcriptomic analysis using an Arabidopsis thaliana oligonucleotide microarray showed clear differences for the expression levels of 26 genes between the two cultivars later in the season. Several defence-related genes showed higher levels of expression in the cultivar that harboured the lowest numbers of herbivores. Our study shows that herbivore community composition develops differentially throughout the season on the two B. oleracea cultivars grown in the field. The correlation between the differences in herbivore communities and differential expression of particular defence-related genes is discussed