Model-based analysis of non-specific binding for background correction of high-density oligonucleotide microarrays

Motivation: High-density DNA microarrays provide us with usefultools for analyzing DNA and RNA comprehensively. However, thebackground signal caused by the non-specific binding (NSB) betweenprobe and target makes it difficult to obtain accurate measurements.To remove the background signal, there is a set of backgroundprobes on Affymetrix Exon arrays to represent the amount ofnon-specific signals, and an accurate estimation of non-specificsignals using these background probes is desirable for improvementof microarray analyses. Results: We developed a thermodynamic model of NSB on shortnucleotide microarrays in which the NSBs are modeled by duplexformation of probes and multiple hypothetical targets. We fittedthe observed signal intensities of the background probes withthose expected by the model to obtain the model parameters.As a result, we found that the presented model can improve theaccuracy of prediction of non-specific signals in comparisonwith previously proposed methods. This result will provide auseful method to correct for the background signal in oligonucleotidemicroarray analysis. Availability: The software is implemented in the R languageand can be downloaded from our website (http://www-shimizu.ist.osaka-u.ac.jp/shimizu_lab/MSNS/). Contact: furusawa{at}ist.osaka-u.ac.jp Supplementary information: Supplementary data are availableat Bioinformatics online. The authors wish it to be known that, in their opinion, thefirst two authors should be regarded as joint First Authors. Associate Editor: Trey Ideker     

From pull-down data to protein interaction networks and complexes with biological relevance

Motivation: Recent improvements in high-throughput Mass Spectrometry(MS) technology have expedited genome-wide discovery of protein–proteininteractions by providing a capability of detecting proteincomplexes in a physiological setting. Computational inferenceof protein interaction networks and protein complexes from MSdata are challenging. Advances are required in developing robustand seamlessly integrated procedures for assessment of protein–proteininteraction affinities, mathematical representation of proteininteraction networks, discovery of protein complexes and evaluationof their biological relevance. Results: A multi-step but easy-to-follow framework for identifyingprotein complexes from MS pull-down data is introduced. It assessesinteraction affinity between two proteins based on similarityof their co-purification patterns derived from MS data. It constructsa protein interaction network by adopting a knowledge-guidedthreshold selection method. Based on the network, it identifiesprotein complexes and infers their core components using a graph-theoreticalapproach. It deploys a statistical evaluation procedure to assessbiological relevance of each found complex. On Saccharomycescerevisiae pull-down data, the framework outperformed othermore complicated schemes by at least 10% in F₁-measure and identified610 protein complexes with high-functional homogeneity basedon the enrichment in Gene Ontology (GO) annotation. Manual examinationof the complexes brought forward the hypotheses on cause offalse identifications. Namely, co-purification of differentprotein complexes as mediated by a common non-protein molecule,such as DNA, might be a source of false positives. Protein identificationbias in pull-down technology, such as the hydrophilic bias couldresult in false negatives. Contact: samatovan{at}ornl.gov Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: Jonathan Wren Present address: Department of Biomedical Informatics, VanderbiltUniversity, Nashville, TN 37232. The authors wish it to be known that, in their opinion, thefirst two authors should be regarded as joint First Authors.     

Power enhancement via multivariate outlier testing with gene expression arrays

Motivation: As the use of microarrays in human studies continuesto increase, stringent quality assurance is necessary to ensureaccurate experimental interpretation. We present a formal approachfor microarray quality assessment that is based on dimensionreduction of established measures of signal and noise componentsof expression followed by parametric multivariate outlier testing. Results: We applied our approach to several data resources.First, as a negative control, we found that the Affymetrix andIllumina contributions to MAQC data were free from outliersat a nominal outlier flagging rate of =0.01. Second, we createda tunable framework for artificially corrupting intensity datafrom the Affymetrix Latin Square spike-in experiment to allowinvestigation of sensitivity and specificity of quality assurance(QA) criteria. Third, we applied the procedure to 507 Affymetrixmicroarray GeneChips processed with RNA from human peripheralblood samples. We show that exclusion of arrays by this approachsubstantially increases inferential power, or the ability todetect differential expression, in large clinical studies. Availability: http://bioconductor.org/packages/2.3/bioc/html/arrayMvout.htmland http://bioconductor.org/packages/2.3/bioc/html/affyContam.htmlaffyContam (credentials: readonly/readonly) Contact: aasare{at}immunetolerance.org; stvjc{at}channing.harvard.edu The authors wish it to be known that, in their opinion, thefirst two authors should be regarded as joint First Authors. Associate Editor: Trey Ideker     

PlasmoGF: an integrated system for comparative genomics and phylogenetic analysis of Plasmodium gene families

Summary: Malaria, one of the world's most common diseases, iscaused by the intracellular protozoan parasite known as Plasmodium.Recently, with the arrival of several malaria parasite genomes,we established an integrated system named PlasmoGF for comparativegenomics and phylogenetic analysis of Plasmodium gene families.Gene families were clustered using the Markov Cluster algorithmimplemented in TribeMCL program and could be searched usingkeywords, gene-family information, domain composition, GeneOntology and BLAST. Moreover, a number of useful bioinformaticstools were implemented to facilitate the analysis of these putativePlasmodium gene families, including gene retrieval, annotation,sequence alignment, phylogeny construction and visualization.In the current version, PlasmoGF contained 8980 sets of genefamilies derived from six malaria parasite genomes: Plasmodium.falciparum, P. berghei, P. knowlesi, P. chabaudi, P. vivax andP. yoelii. The availability of such a highly integrated systemwould be of great interest for the community of researchersworking on malaria parasite phylogenomics. Availability: PlasmoGF is freely available at http://bioinformatics.zj.cn/pgf/ Contact: xiaokunli{at}163.net; baoqy{at}genomics.org.cn; fuz3{at}psu.edu Associate Editor: Jonathan Wren The authors wish it to be known that, in their opinion, thefirst two authors should be regarded as joint First Authors.     

BicOverlapper: a tool for bicluster visualization

Summary: BicOverlapper is a tool to visualize biclusters fromgene-expression matrices in a way that helps to compare biclusteringmethods, to unravel trends and to highlight relevant genes andconditions. A visual approach can complement biological andstatistical analysis and reduce the time spent by specialistsinterpreting the results of biclustering algorithms. The techniqueis based on a force-directed graph where biclusters are representedas flexible overlapped groups of genes and conditions. Availability: The BicOverlapper software and supplementary materialare available at http://vis.usal.es/bicoverlapper Contact: rodri{at}usal.es Associate Editor: John Quackenbush The first two authors should be reported as joint first authors.     

A servo-controlled respiration system for inhalation studies in anesthetized animals

Rosenthal, Frank S., and Changhong Li. Aservo-controlled respiration system for inhalation studies inanesthetized animals. J. Appl.Physiol. 83(5): 1768-1774, 1997.To facilitate aerosol deposition experiments and aerosol exposures in anesthetized animals, a servo-controlled respiration system was developed and tested. The system induces ventilation by varying extrathoracic pressure in a whole body respirator in which an intubated animal isplaced. The pressure inside the whole body respirator is varied with athree-way servo-controlled spool valve connected to sources of positiveand negative pressure. A computer-based system detects respiratory flowand computes the controlling signal for the valve by using aproportional-integral-derivative algorithm, to achieve desired patternsof flow and volume vs. time. The system was used with dogs and found toaccurately induce various single-breath breathing patterns involvingconstant-flow inspirations and expirations as well as breath-holdperiods. A similar system was used to induced repeated breaths withdesired parameters for continuous exposure to particles and forventilation of animals between experiments.

     

A comparison of meta-analysis methods for detecting differentially expressed genes in microarray experiments

Motivation: The proliferation of public data repositories createsa need for meta-analysis methods to efficiently evaluate, integrateand validate related datasets produced by independent groups.A t-based approach has been proposed to integrate effect sizefrom multiple studies by modeling both intra- and between-studyvariation. Recently, a non-parametric ‘rank product’method, which is derived based on biological reasoning of fold-changecriteria, has been applied to directly combine multiple datasetsinto one meta study. Fisher's Inverse ² method, which only dependson P-values from individual analyses of each dataset, has beenused in a couple of medical studies. While these methods addressthe question from different angles, it is not clear how theycompare with each other. Results: We comparatively evaluate the three methods; t-basedhierarchical modeling, rank products and Fisher's Inverse ²test with P-values from either the t-based or the rank productmethod. A simulation study shows that the rank product method,in general, has higher sensitivity and selectivity than thet-based method in both individual and meta-analysis, especiallyin the setting of small sample size and/or large between-studyvariation. Not surprisingly, Fisher's ² method highly dependson the method used in the individual analysis. Application toreal datasets demonstrates that meta-analysis achieves morereliable identification than an individual analysis, and rankproducts are more robust in gene ranking, which leads to a muchhigher reproducibility among independent studies. Though t-basedmeta-analysis greatly improves over the individual analysis,it suffers from a potentially large amount of false positiveswhen P-values serve as threshold. We conclude that careful meta-analysisis a powerful tool for integrating multiple array studies. Contact: fxhong{at}jimmy.harvard.edu Supplementary information: Supplementary data are availableat Bioinformatics online. Associate Editor: David Rocke Present address: Department of Biostatistics and ComputationalBiology, Dana-Farber Cancer Institute, Harvard School of PublicHealth, 44 Binney Street, Boston, MA 02115, USA.     

ISD: a software package for Bayesian NMR structure calculation

Summary: The conventional approach to calculating biomolecularstructures from nuclear magnetic resonance (NMR) data is oftenviewed as subjective due to its dependence on rules of thumbfor deriving geometric constraints and suitable values for theoryparameters from noisy experimental data. As a result, it canbe difficult to judge the precision of an NMR structure in anobjective manner. The inferential Structure determination (ISD)framework, which has been introduced recently, addresses thisproblem by using Bayesian inference to derive a probabilitydistribution that represents both the unknown structure andits uncertainty. It also determines additional unknowns, suchas theory parameters, that normally need to be chosen empirically.Here we give an overview of the ISD software package, whichimplements this methodology. Availability: http://www.bioc.cam.ac.uk/isd Contact: wolfgang.rieping{at}bioc.cam.ac.uk, michael.habeck{at}tuebingen.mpg.de Associate Editor: Alfonso Valencia     

Co-operativity in non-covalent interactions in ternary complexes: a comprehensive electronic structure theory based investigation

The structure and stability of various ternary complexes in which an extended aromatic system such as coronene interacts with ions/atoms/molecules on opposite faces of the π-electron cloud were investigated using ab initio calculations. By characterizing the nature of the intermolecular interactions using an energy decomposition analysis, it was shown that there is an interplay between various types of interactions and that there are co-operativity effects, particularly when different types of interactions coexist in the same system.

Open image in new window

Graphical abstract Weak OH-π, π-π and van der Waals-π ternary systems are stabilized through dispersion interactions. Cation-π ternary systems are stabilized by through-space electrostatic interactions.

     

Effects of blue light deficiency on acclimation of light energy partitioning in PSII and CO2 assimilation capacity to high irradiance in spinach leaves

Blue light effects on the acclimation of energy partitioningcharacteristics in PSII and CO₂ assimilation capacity in spinachto high growth irradiance were investigated. Plants were grownhydroponically in different light treatments that were a combinationof two light qualities and two irradiances, i.e. white lightand blue-deficient light at photosynthetic photon flux densities(PPFDs) of 100 and 500 µmol m^–2 s^–1. The CO₂assimilation rate, the quantum efficiency of PSII (PSII) andthermal dissipation activity / in young, fully expanded leaves were measured under 1,600 µmol m^–2 s^–1white light. The CO₂ assimilation rate and PSII were higher,while / was lower in plants grown under high irradiancethan in plants grown under low irradiance. These responses wereobserved irrespective of the presence or absence of blue lightduring growth. The extent of the increase in the CO₂ assimilationrate and PSII and the decrease in / by high growth irradiance was smaller under blue light-deficient conditions. These resultsindicate that blue light helps to boost the acclimation responsesof energy partitioning in PSII and CO₂ assimilation to highirradiance. Similarly, leaf N, Cyt f and Chl contents per unitleaf area increased by high growth irradiance, and the extentof the increment in leaf N, Cyt f and Chl was smaller underblue light-deficient conditions. Regression analysis showedthat the differences in energy partitioning in PSII and CO₂assimilation between plants grown under high white light andhigh blue-deficient light were closely related to the differencein leaf N.     

Comparative molecular field analysis and molecular dynamics studies of α/β hydrolase domain containing 6 (ABHD6) inhibitors

The endocannabinoid system remains an attractive molecular target for pharmacological intervention due to its roles in the central nervous system in learning, thinking, emotional function, regulation of food intake or pain sensation, as well as in the peripheral nervous system, where it modulates the action of cardiovascular, immune, metabolic or reproductive function. α/β hydrolase domain containing 6 (ABHD6)—an enzyme forming part of the endocannabinoid system—is a newly discovered post-genomic protein acting as a 2-AG (2-arachidonoylglycerol) serine hydrolase. We have recently reported a series of 1,2,5-thiadiazole carbamates as potent and selective ABHD6 inhibitors. Here, we present comparative molecular field analysis (CoMFA) and molecular dynamics studies of these compounds. First, we performed a homology modeling study of ABHD6 based on the assumption that the catalytic triad of ABHD6 comprises Ser148–His306–Asp 278 and the oxyanion hole is formed by Met149 and Phe80. A total of 42 compounds was docked to the homology model using the Glide module from the Schrödinger suite of software and the selected docking poses were used for CoMFA alignment. A model with the following statistics was obtained: R ² = 0.98, Q ² = 0.55. In order to study the molecular interactions of the inhibitors with ABHD6 in detail, molecular dynamics was performed with the Desmond program. It was found that, during the simulations, the hydrogen bond between the inhibitor carbonyl group and the main chain of Phe80 is weakened, whereas a new hydrogen bond with the side chain of Ser148 is formed, facilitating the possible formation of a covalent bond.

Open image in new window

Graphical Abstract Left–right: Docking pose of 1 in the binding pocket of α/β hydrolase domain containing 6 (ABHD6) selected for molecular alignment; CoMFA steric and electrostatic contour fields; changes in potential energy of the complex during simulations for the complex of 6 and ABHD6

     

Phosphorus: a limiting nutrient for humanity?

1. Download : Download high-res image (82KB)
2. Download : Download full-size image

Highlights? P is essential to build key molecules, such as RNA needed for rapid growth. ? P is often limiting to productivity in ecosystems, including to crops. ? Large amounts of P are mined annually to make fertilizer to produce food. ? Current P use is unsustainable owing to ecosystem impacts and uncertainties about fertilizer production. ? A food system redesign to combine P-efficient crops with P recycling is needed.     

Slider--maximum use of probability information for alignment of short sequence reads and SNP detection

Motivation: A plethora of alignment tools have been createdthat are designed to best fit different types of alignment conditions.While some of these are made for aligning Illumina SequenceAnalyzer reads, none of these are fully utilizing its probability(prb) output. In this article, we will introduce a new alignmentapproach (Slider) that reduces the alignment problem space byutilizing each read base's probabilities given in the prb files. Results: Compared with other aligners, Slider has higher alignmentaccuracy and efficiency. In addition, given that Slider matchesbases with probabilities other than the most probable, it significantlyreduces the percentage of base mismatches. The result is thatits SNP predictions are more accurate than other SNP predictionapproaches used today that start from the most probable sequence,including those using base quality. Contact: nmalhis{at}bcgsc.ca Supplementary information and availability: http://www.bcgsc.ca/platform/bioinfo/software/slider Associate Editor: Dmitrij Frishman     

Ghrelin regulates phasic dopamine and nucleus accumbens signaling evoked by food‐predictive stimuli

Environmental stimuli that signal food availability hold powerful sway over motivated behavior and promote feeding, in part, by activating the mesolimbic system. These food‐predictive cues evoke brief (phasic) changes in nucleus accumbens (NAc) dopamine concentration and in the activity of individual NAc neurons. Phasic fluctuations in mesolimbic signaling have been directly linked to goal‐directed behaviors, including behaviors elicited by food‐predictive cues. Food‐seeking behavior is also strongly influenced by physiological state (i.e., hunger vs. satiety). Ghrelin, a stomach hormone that crosses the blood‐brain barrier, is linked to the perception of hunger and drives food intake, including intake potentiated by environmental cues. Notwithstanding, whether ghrelin regulates phasic mesolimbic signaling evoked by food‐predictive stimuli is unknown. Here, rats underwent Pavlovian conditioning in which one cue predicted the delivery of rewarding food (CS+) and a second cue predicted nothing (CS?). After training, we measured the effect of ghrelin infused into the lateral ventricle (LV) on sub‐second fluctuations in NAc dopamine using fast‐scan cyclic voltammetry and individual NAc neuron activity using in vivo electrophysiology in separate groups of rats. LV ghrelin augmented both phasic dopamine and phasic increases in the activity of NAc neurons evoked by the CS+. Importantly, ghrelin did not affect the dopamine nor NAc neuron response to the CS?, suggesting that ghrelin selectively modulated mesolimbic signaling evoked by motivationally significant stimuli. These data demonstrate that ghrelin, a hunger signal linked to physiological state, can regulate cue‐evoked mesolimbic signals that underlie food‐directed behaviors.

     

Identification of genes expressed in vascular tissues using NPA-induced vascular overgrowth in Arabidopsis

The genetic basis of vascular differentiation and function isrelatively poorly understood, partly due to the difficulty ofscreening for mutants defective in internal vascular tissues.Here we present an approach based on a predicted increase invascular-related gene expression in response to an auxin transportinhibitor-induced vascular overgrowth. We used microarray analysesto identify 336 genes that were up-regulated 2-fold in shoottissues of Arabidopsis thaliana showing vascular overgrowth.Promoter–marker gene fusions revealed that 38 out of 40genes with 4-fold up-regulation in vascular overgrowth tissueshad vascular-related expression in transgenic Arabidopsis plants.Obtained expression patterns included cambial tissues and differentiatingxylem, phloem and fibers. A total of 15 genes were found tohave vascular-specific expression patterns in the leaves and/orinflorescence stems. This study provides empirical evidenceof the efficiency of the approach and describes for the firsttime the in situ expression patterns of the majority of theassessed genes.     

The encapsulation of the gemcitabine anticancer drug into grapheme nest: a theoretical study

The efficient transport of a drug molecule until its target cell constitutes a significant challenge for delivery processes. To achieve such objectives, solid nanocapsules that protect the immune system during the transport should be developed and controlled at the nanoscale level. From this point of view, nanostructures based on graphene sheets could present some promising properties due to their ultimate size and dimension. In this work, we present theoretical results using DFT calculations, dealing with a graphene-based delivery system. Indeed, we demonstrate the stability of the gemcitabine anticancer molecule when it is encapsulated into two concave graphene sheets organized as a nest. Quantum calculations showed that the most stable state is located inside the nest, which is then formed by two layers distanced 6 Å from each other. For all the optimized systems, we focused on the dependence of the interaction energy on the molecule displacements during its entrance in the graphene nest and its exit from it. We also analyzed their consequence on the local morphological and electronic charge properties.

Open image in new window

Graphical Abstract Adsorption energy (in eV) of gemcitabine drug during its encapsulation inside the nest of grapheneand its release from it