Nonequilibrium effects in DNA microarrays: a multiplatform study

It has recently been shown that in some DNA microarrays the time needed to reach thermal equilibrium may largely exceed the typical experimental time, which is about 15h in standard protocols (Hooyberghs et al. Phys. Rev. E 81, 012901 (2010)). In this paper we discuss how this breakdown of thermodynamic equilibrium could be detected in microarray experiments without resorting to real time hybridization data, which are difficult to implement in standard experimental conditions. The method is based on the analysis of the distribution of fluorescence intensities I from different spots for probes carrying base mismatches. In thermal equilibrium and at sufficiently low concentrations, log I is expected to be linearly related to the hybridization free energy $\Delta G$ with a slope equal to $1/RT_{exp}$, where $T_{exp}$ is the experimental temperature and R is the gas constant. The breakdown of equilibrium results in the deviation from this law. A model for hybridization kinetics explaining the observed experimental behavior is discussed, the so-called 3-state model. It predicts that deviations from equilibrium yield a proportionality of $\log I$ to $\Delta G/RT_{eff}$. Here, $T_{eff}$ is an effective temperature, higher than the experimental one. This behavior is indeed observed in some experiments on Agilent arrays. We analyze experimental data from two other microarray platforms and discuss, on the basis of the results, the attainment of equilibrium in these cases. Interestingly, the same 3-state model predicts a (dynamical) saturation of the signal at values below the expected one at equilibrium.Comment: 27 pages, 9 figures, 1 tabl

Decoding of Superimposed Traces Produced by Direct Sequencing of Heterozygous Indels

Direct Sanger sequencing of a diploid template containing a heterozygous insertion or deletion results in a difficult-to-interpret mixed trace formed by two allelic traces superimposed onto each other. Existing computational methods for deconvolution of such traces require knowledge of a reference sequence or the availability of both direct and reverse mixed sequences of the same template. We describe a simple yet accurate method, which uses dynamic programming optimization to predict superimposed allelic sequences solely from a string of letters representing peaks within an individual mixed trace. We used the method to decode 104 human traces (mean length 294 bp) containing heterozygous indels 5 to 30 bp with a mean of 99.1% bases per allelic sequence reconstructed correctly and unambiguously. Simulations with artificial sequences have demonstrated that the method yields accurate reconstructions when (1) the allelic sequences forming the mixed trace are sufficiently similar, (2) the analyzed fragment is significantly longer than the indel, and (3) multiple indels, if present, are well-spaced. Because these conditions occur in most encountered DNA sequences, the method is widely applicable. It is available as a free Web application Indelligent at http://ctap.inhs.uiuc.edu/dmitriev/indel.asp

Novel computational methods for increasing PCR primer design effectiveness in directed sequencing

<p>Abstract</p> <p>Background</p> <p>Polymerase chain reaction (PCR) is used in directed sequencing for the discovery of novel polymorphisms. As the first step in PCR directed sequencing, effective PCR primer design is crucial for obtaining high-quality sequence data for target regions. Since current computational primer design tools are not fully tuned with stable underlying laboratory protocols, researchers may still be forced to iteratively optimize protocols for failed amplifications after the primers have been ordered. Furthermore, potentially identifiable factors which contribute to PCR failures have yet to be elucidated. This inefficient approach to primer design is further intensified in a high-throughput laboratory, where hundreds of genes may be targeted in one experiment.</p> <p>Results</p> <p>We have developed a fully integrated computational PCR primer design pipeline that plays a key role in our high-throughput directed sequencing pipeline. Investigators may specify target regions defined through a rich set of descriptors, such as Ensembl accessions and arbitrary genomic coordinates. Primer pairs are then selected computationally to produce a minimal amplicon set capable of tiling across the specified target regions. As part of the tiling process, primer pairs are computationally screened to meet the criteria for success with one of two PCR amplification protocols. In the process of improving our sequencing success rate, which currently exceeds 95% for exons, we have discovered novel and accurate computational methods capable of identifying primers that may lead to PCR failures. We reveal the laboratory protocols and their associated, empirically determined computational parameters, as well as describe the novel computational methods which may benefit others in future primer design research.</p> <p>Conclusion</p> <p>The high-throughput PCR primer design pipeline has been very successful in providing the basis for high-quality directed sequencing results and for minimizing costs associated with labor and reprocessing. The modular architecture of the primer design software has made it possible to readily integrate additional primer critique tests based on iterative feedback from the laboratory. As a result, the primer design software, coupled with the laboratory protocols, serves as a powerful tool for low and high-throughput primer design to enable successful directed sequencing.</p

Paroxysmal exercise-induced dyskinesia and epilepsy is due to mutations in SLC2A1, encoding the glucose transporter GLUT1

Paroxysmal exercise-induced dyskinesia (PED) can occur in isolation or in association with epilepsy, but the genetic causes and pathophysiological mechanisms are still poorly understood. We performed a clinical evaluation and genetic analysis in a five-generation family with co-occurrence of PED and epilepsy (n = 39), suggesting that this combination represents a clinical entity. Based on a whole genome linkage analysis we screened SLC2A1, encoding the glucose transporter of the blood-brain-barrier, GLUT1 and identified heterozygous missense and frameshift mutations segregating in this and three other nuclear families with a similar phenotype. PED was characterized by choreoathetosis, dystonia or both, affecting mainly the legs. Predominant epileptic seizure types were primary generalized. A median CSF/blood glucose ratio of 0.52 (normal >0.60) in the patients and a reduced glucose uptake by mutated transporters compared with the wild-type as determined in Xenopus oocytes confirmed a pathogenic role of these mutations. Functional imaging studies implicated alterations in glucose metabolism in the corticostriate pathways in the pathophysiology of PED and in the frontal lobe cortex in the pathophysiology of epileptic seizures. Three patients were successfully treated with a ketogenic diet. In conclusion, co-occurring PED and epilepsy can be due to autosomal dominant heterozygous SLC2A1 mutations, expanding the phenotypic spectrum associated with GLUT1 deficiency and providing a potential new treatment option for this clinical syndrome

A Deletion in Exon 9 of the LIPH Gene Is Responsible for the Rex Hair Coat Phenotype in Rabbits (Oryctolagus cuniculus)

The fur of common rabbits is constituted of 3 types of hair differing in length and diameter while that of rex animals is essentially made up of amazingly soft down-hair. Rex short hair coat phenotypes in rabbits were shown to be controlled by three distinct loci. We focused on the “r1” mutation which segregates at a simple autosomal-recessive locus in our rabbit strains. A positional candidate gene approach was used to identify the rex gene and the corresponding mutation. The gene was primo-localized within a 40 cM region on rabbit chromosome 14 by genome scanning families of 187 rabbits in an experimental mating scheme. Then, fine mapping refined the region to 0.5 cM (Z = 78) by genotyping an additional 359 offspring for 94 microsatellites present or newly generated within the first defined interval. Comparative mapping pointed out a candidate gene in this 700 kb region, namely LIPH (Lipase Member H). In humans, several mutations in this major gene cause alopecia, hair loss phenotypes. The rabbit gene structure was established and a deletion of a single nucleotide was found in LIPH exon 9 of rex rabbits (1362delA). This mutation results in a frameshift and introduces a premature stop codon potentially shortening the protein by 19 amino acids. The association between this deletion and the rex phenotype was complete, as determined by its presence in our rabbit families and among a panel of 60 rex and its absence in all 60 non-rex rabbits. This strongly suggests that this deletion, in a homozygous state, is responsible for the rex phenotype in rabbits

A Deletion in Exon 9 of the LIPH Gene Is Responsible for the Rex Hair Coat Phenotype in Rabbits (Oryctolagus cuniculus)

The fur of common rabbits is constituted of 3 types of hair differing in length and diameter while that of rex animals is essentially made up of amazingly soft down-hair. Rex short hair coat phenotypes in rabbits were shown to be controlled by three distinct loci. We focused on the “r1” mutation which segregates at a simple autosomal-recessive locus in our rabbit strains. A positional candidate gene approach was used to identify the rex gene and the corresponding mutation. The gene was primo-localized within a 40 cM region on rabbit chromosome 14 by genome scanning families of 187 rabbits in an experimental mating scheme. Then, fine mapping refined the region to 0.5 cM (Z = 78) by genotyping an additional 359 offspring for 94 microsatellites present or newly generated within the first defined interval. Comparative mapping pointed out a candidate gene in this 700 kb region, namely LIPH (Lipase Member H). In humans, several mutations in this major gene cause alopecia, hair loss phenotypes. The rabbit gene structure was established and a deletion of a single nucleotide was found in LIPH exon 9 of rex rabbits (1362delA). This mutation results in a frameshift and introduces a premature stop codon potentially shortening the protein by 19 amino acids. The association between this deletion and the rex phenotype was complete, as determined by its presence in our rabbit families and among a panel of 60 rex and its absence in all 60 non-rex rabbits. This strongly suggests that this deletion, in a homozygous state, is responsible for the rex phenotype in rabbits

Position dependent mismatch discrimination on DNA microarrays – experiments and model

<p>Abstract</p> <p>Background</p> <p>The propensity of oligonucleotide strands to form stable duplexes with complementary sequences is fundamental to a variety of biological and biotechnological processes as various as microRNA signalling, microarray hybridization and PCR. Yet our understanding of oligonucleotide hybridization, in particular in presence of surfaces, is rather limited. Here we use oligonucleotide microarrays made in-house by optically controlled DNA synthesis to produce probe sets comprising all possible single base mismatches and base bulges for each of 20 sequence motifs under study.</p> <p>Results</p> <p>We observe that mismatch discrimination is mostly determined by the defect position (relative to the duplex ends) as well as by the sequence context. We investigate the thermodynamics of the oligonucleotide duplexes on the basis of double-ended molecular zipper. Theoretical predictions of defect positional influence as well as long range sequence influence agree well with the experimental results.</p> <p>Conclusion</p> <p>Molecular zipping at thermodynamic equilibrium explains the binding affinity of mismatched DNA duplexes on microarrays well. The position dependent nearest neighbor model (PDNN) can be inferred from it. Quantitative understanding of microarray experiments from first principles is in reach.</p

An efficient algorithm for the stochastic simulation of the hybridization of DNA to microarrays

<p>Abstract</p> <p>Background</p> <p>Although oligonucleotide microarray technology is ubiquitous in genomic research, reproducibility and standardization of expression measurements still concern many researchers. Cross-hybridization between microarray probes and non-target ssDNA has been implicated as a primary factor in sensitivity and selectivity loss. Since hybridization is a chemical process, it may be modeled at a population-level using a combination of material balance equations and thermodynamics. However, the hybridization reaction network may be exceptionally large for commercial arrays, which often possess at least one reporter per transcript. Quantification of the kinetics and equilibrium of exceptionally large chemical systems of this type is numerically infeasible with customary approaches.</p> <p>Results</p> <p>In this paper, we present a robust and computationally efficient algorithm for the simulation of hybridization processes underlying microarray assays. Our method may be utilized to identify the extent to which nucleic acid targets (e.g. cDNA) will cross-hybridize with probes, and by extension, characterize probe robustnessusing the information specified by MAGE-TAB. Using this algorithm, we characterize cross-hybridization in a modified commercial microarray assay.</p> <p>Conclusions</p> <p>By integrating stochastic simulation with thermodynamic prediction tools for DNA hybridization, one may robustly and rapidly characterize of the selectivity of a proposed microarray design at the probe and "system" levels. Our code is available at <url>http://www.laurenzi.net</url>.</p

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