Secondary structure in the target as a confounding factor in synthetic oligomer microarray design

BACKGROUND: Secondary structure in the target is a property not usually considered in software applications for design of optimal custom oligonucleotide probes. It is frequently assumed that eliminating self-complementarity, or screening for secondary structure in the probe, is sufficient to avoid interference with hybridization by stable secondary structures in the probe binding site. Prediction and thermodynamic analysis of secondary structure formation in a genome-wide set of transcripts from Brucella suis 1330 demonstrates that the properties of the target molecule have the potential to strongly influence the rate and extent of hybridization between transcript and tethered oligonucleotide probe in a microarray experiment. RESULTS: Despite the relatively high hybridization temperatures and 1M monovalent salt imposed in the modeling process to approximate hybridization conditions used in the laboratory, we find that parts of the target molecules are likely to be inaccessible to intermolecular hybridization due to the formation of stable intramolecular secondary structure. For example, at 65°C, 28 ± 7% of the average cDNA target sequence is predicted to be inaccessible to hybridization. We also analyzed the specific binding sites of a set of 70mer probes previously designed for Brucella using a freely available oligo design software package. 21 ± 13% of the nucleotides in each probe binding site are within a double-stranded structure in over half of the folds predicted for the cDNA target at 65°C. The intramolecular structures formed are more stable and extensive when an RNA target is modeled rather than cDNA. When random shearing of the target is modeled for fragments of 200, 100 and 50 nt, an overall destabilization of secondary structure is predicted, but shearing does not eliminate secondary structure. CONCLUSION: Secondary structure in the target is pervasive, and a significant fraction of the target is found in double stranded conformations even at high temperature. Stable structure in the target has the potential to interfere with hybridization and should be a factor in interpretation of microarray results, as well as an explicit criterion in array design. Inclusion of this property in an oligonucleotide design procedure would change the definition of an optimal oligonucleotide significantly

Comparative Genomic Analysis of Two Vibrio toranzoniae Strains with Different Virulence Capacity Reveals Clues on Its Pathogenicity for Fish

Vibrio toranzoniae is a Gram-negative bacterium of the Splendidus clade within the Vibrio genus. V. toranzoniae was first isolated from healthy clams in Galicia (Spain) but recently was also identified associated to disease outbreaks of red conger eel in Chile. Experimental challenges showed that the Chilean isolates were able to produce fish mortalities but not the strains isolated from clams. The aim of the present study was to determine the differences at the genomic level between the type strain of the species (CECT 7225T) and the strain R17, isolated from red conger eel in Chile, which could explain their different virulent capacity. The genome-based comparison showed high homology between both strains but differences were observed in certain gene clusters that include some virulence factors. Among these, we found that iron acquisition systems and capsule synthesis genes were the main differential features between both genomes that could explain the differences in the pathogenicity of the strains. Besides, the studied genomes presented genomic islands and toxins, and the R17 strain presented CRISPR sequences that are absent on the type strain. Taken together, this analysis provided important insights into virulence factors of V. toranzoniae that will lead to a better understanding of the pathogenic processThis research was supported in part by grants AGL2013-42628R and AGL2016-77539R from Ministerio de Economia y Competitividad (Spain)S

Background correction using dinucleotide affinities improves the performance of GCRMA

BACKGROUND: High-density short oligonucleotide microarrays are a primary research tool for assessing global gene expression. Background noise on microarrays comprises a significant portion of the measured raw data, which can have serious implications for the interpretation of the generated data if not estimated correctly. RESULTS: We introduce an approach to calculate probe affinity based on sequence composition, incorporating nearest-neighbor (NN) information. Our model uses position-specific dinucleotide information, instead of the original single nucleotide approach, and adds up to 10% to the total variance explained (R(2)) when compared to the previously published model. We demonstrate that correcting for background noise using this approach enhances the performance of the GCRMA preprocessing algorithm when applied to control datasets, especially for detecting low intensity targets. CONCLUSION: Modifying the previously published position-dependent affinity model to incorporate dinucleotide information significantly improves the performance of the model. The dinucleotide affinity model enhances the detection of differentially expressed genes when implemented as a background correction procedure in GeneChip preprocessing algorithms. This is conceptually consistent with physical models of binding affinity, which depend on the nearest-neighbor stacking interactions in addition to base-pairing

Molecular targets for rapid identification of Brucella spp

BACKGROUND: Brucella is an intracellular pathogen capable of infecting animals and humans. There are six recognized species of Brucella that differ in their host preference. The genomes of the three Brucella species have been recently sequenced. Comparison of the three revealed over 98% sequence similarity at the protein level and enabled computational identification of common and differentiating genes. We validated these computational predictions and examined the expression patterns of the putative unique and differentiating genes, using genomic and reverse transcription PCR. We then screened a set of differentiating genes against classical Brucella biovars and showed the applicability of these regions in the design of diagnostic tests. RESULTS: We have identified and tested set of molecular targets that are associated in unique patterns with each of the sequenced Brucella spp. A comprehensive comparison was made among the published genome sequences of B. abortus, B. melitensis and B. suis. The comparison confirmed published differences between the three Brucella genomes, and identified subsets of features that were predicted to be of interest in a functional comparison of B. melitensis and B. suis to B. abortus. Differentiating sequence regions from B. abortus, B. melitensis and B. suis were used to develop PCR primers to test for the existence and in vitro transcription of these genes in these species. Only B. suis is found to have a significant number of unique genes, but combinations of genes and regions that exist in only two out of three genomes and are therefore useful for diagnostics were identified and confirmed. CONCLUSION: Although not all of the differentiating genes identified were transcribed under steady state conditions, a group of genes sufficient to discriminate unambiguously between B. suis, B. melitensis, and B. abortus was identified. We present an overview of these genomic differences and the use of these features to discriminate among a number of Brucella biovars

Equifinality and preservation potential of complex eskers

Eskers are useful for reconstructing meltwater drainage systems of glaciers and ice sheets. However, our process understanding of eskers suffers from a disconnect between sporadic detailed morpho‐sedimentary investigations of abundant large‐scale ancient esker systems, and a small number of modern analogues where esker formation has been observed. This paper presents the results of detailed field and high‐resolution remote sensing studies into two esker systems that have recently emerged at Hørbyebreen, Svalbard, and one at Breiðamerkurjökull, Iceland. Despite the different glaciological settings (polythermal valley glacier vs. active temperate piedmont lobe), in all cases a distinctive planform morphology has developed, where ridges are orientated in two dominant directions corresponding to the direction of ice flow and the shape of the ice margin. These two orientations in combination form a cross‐cutting and locally rectilinear pattern. One set of ridges at Hørbyebreen is a hybrid of eskers and geometric ridges formed during a surge and/or jökulhlaup event. The other sets of ridges are eskers formed time‐transgressively at a retreating ice margin. The similar morphology of esker complexes formed in different ways on both glacier forelands implies equifinality, meaning that care should be taken when interpreting Quaternary esker patterns. The eskers at Hørbyebreen contain substantial ice‐cores with a high ice:sediment ratio, suggesting that they would be unlikely to survive after ice melt. The Breiðamerkurjökull eskers emerged from terrain characterized by buried ice that has melted out. Our observations lead us to conclude that eskers may reflect a wide range of processes at dynamic ice margins, including significant paraglacial adjustments. This work, as well as previous studies, confirms that constraints on esker morphology include: topographic setting (e.g. confined valley or broad plain); sediment and meltwater availability (including surges and jökulhlaups); position of formation (supraglacial, englacial or subglacial); and ice‐marginal dynamics such as channel abandonment, the formation of outwash heads or the burial and/or exhumation of dead ice

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

Sequencing by Hybridization of Long Targets

Sequencing by Hybridization (SBH) reconstructs an n-long target DNA sequence from its biochemically determined l-long subsequences. In the standard approach, the length of a uniformly random sequence that can be unambiguously reconstructed is limited to due to repetitive subsequences causing reconstruction degeneracies. We present a modified sequencing method that overcomes this limitation without the need for different types of biochemical assays and is robust to error

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038−2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pKa values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pKa values in closer agreement to those obtained experimentally for the more slowly converging residues

Accurate Estimates of Microarray Target Concentration from a Simple Sequence-Independent Langmuir Model

Background: Microarray technology is a commonly used tool for assessing global gene expression. Many models for estimation of target concentration based on observed microarray signal have been proposed, but, in general, these models have been complex and platform-dependent. Principal Findings: We introduce a universal Langmuir model for estimation of absolute target concentration from microarray experiments. We find that this sequence-independent model, characterized by only three free parameters, yields excellent predictions for four microarray platforms, including Affymetrix, Agilent, Illumina and a custom-printed microarray. The model also accurately predicts concentration for the MAQC data sets. This approach significantly reduces the computational complexity of quantitative target concentration estimates. Conclusions: Using a simple form of the Langmuir isotherm model, with a minimum of parameters and assumptions, and without explicit modeling of individual probe properties, we were able to recover absolute transcript concentrations with high R 2 on four different array platforms. The results obtained here suggest that with a ‘‘spiked-in’ ’ concentration serie

Spark plasma sintering of TiyNb1-yCxN1-x monolithic ceramics obtained by mechanically induced self-sustaining reaction

Nanometer-sized titanium-niobium carbonitride powders (TiyNb1-yCxN1-x) with different Ti/Nb atomic ratios were obtained by a mechanically induced self-sustaining reaction, and sintered by spark plasma sintering technique at 1500 degrees C for 1 min in a vacuum atmosphere. Mechanical properties such as hardness and Young's modulus were determined by nanoindentation technique and friction and wear coefficients assessed by ball-on-disk testing using alumina ball in dry sliding conditions. The fracture surface and wear tracks of samples were examined by scanning electron microscopy. Results showed that it is possible to obtain dense monolithic ceramics from the solid solution (TiyNb1-yCxN1-x) with good mechanical properties and excellent wear resistance. The optimum values of nanomechanical properties were found for the Ti0.3Nb0.7C0.5N0.5 ceramic composition, which exhibited a high hardness over 26.0 GPa and Young's modulus around 400 GPa. Crown Copyright (C) 2012 Published by Elsevier B.V. All rights reserved.The authors would like to thank Dr. Emilio Rayon for performing the nanoindentation analysis in the Materials Technology institute (ITM) of the Polytechnic University of Valencia and to acknowledge the financial support received from Spanish Ministry of Science and Innovation for FPI grant (MAT2006-01783). This work was supported by the Spanish government under grant (MAT2010-17046), which is financed in part by the European Regional Development Fund of 2007-2013. E. Chicardi was supported by CSIC through a JAE-Pre grant, which is financed in part by the European Social Fund (ESF).Borrell Tomás, MA.; Salvador Moya, MD.; Garcia-Rocha, V.; Fernandez, A.; Chicardi, E.; Gotor, FJ. (2012). Spark plasma sintering of TiyNb1-yCxN1-x monolithic ceramics obtained by mechanically induced self-sustaining reaction. Materials Science and Engineering: A. 543:173-179. https://doi.org/10.1016/j.msea.2012.02.071S17317954