Multifactorial screening design and analysis of SELDI-TOF ProteinChip array optimization experiments

Surface-enhanced laser desorption/ionization time-of-flight mass spectrometry is a powerful tool for rapidly generating protein expression data (peptide and protein profiles) from a large number of samples. However, as with any technology, it must be optimized and reproducible for one to have confidence in the results. Using a classical statistical method called the fractional factorial design of experiments, we assessed the effects of 11 different experimental factors. We also developed several metrics that reflect trace quality and reproducibility. These were used to measure the effect of each individual factor, and the interactions between factors, to determine optimal factor settings and thus ultimately produce the best possible traces. Significant improvements to output traces were seen by simultaneously altering several parameters, either in the sample preparation procedure or during the matrix preparation and application procedure. This has led to the implementation of an improved method that gives a better quality, reproducible, and robust output.     

A logic-reasoning based system to harness bioprocess experimental data and knowledge for design

Bioprocess design requires substantial resources during the experimental investigation of the options for each bioprocess step. This is both time-consuming and expensive. The amount of data available has increased exponentially since the expansion of new biological drug development. Data are heterogeneous, sometimes inconsistent and incomplete, making them difficult to be systematically utilised for analysis for any new bioprocess design. In this paper, we report a novel computational method that harnesses the bioprocess experimental data to assist design decision making, and perhaps identify further needed experiments. First, we develop a new data representation structure to capture the experimental data systematically. Then the ontology for modelling the relationship of data properties is created. A computational system has been developed to search relevant data, or to predict required process conditions, or to suggest a new set of experiments for process development. A prototype for harnessing centrifugation experimental data has been built, and is then used to illustrate the method and demonstrate the type of results that can be obtained. Evaluations show that such a system has significant potential to mine the relevant experimental data to assist new drug bioprocess development, which should reduce process development time and cost.     

Evaluation of experimental designs for two-color cDNA microarrays.

The major goal of two-color cDNA microarray experiments is to measure the relative gene expression level (i.e., relative amount of mRNA) of each gene between samples in studies of gene expression. More specifically, given an N-sample experiment, we need all N(N - 1)/2 relative expression levels of all sample pairs of each gene for identification of the differentially expressed genes and for clustering of gene expression patterns. However, the intensities observed from two-color cDNA microarray experiments do not simply represent the relative gene expression level. They are composed of signal (gene expression level), noise, and other factors. In discussions on the experimental design of two-color cDNA microarray experiments, little attention has been given to the fact that different combinations of test and control samples will produce microarray intensities data with varying intrinsic composition of factors. As a consequence, not all experimental designs for two-color cDNA microarray experiments are able to provide all possible relative gene expression levels. This phenomenon has never been addressed. To obtain all possible relative gene expression levels, a novel method for two-color cDNA microarray experimental design evaluation is necessary that will allow the making of an accurate choice. In this study, we propose a model-based approach to illustrate how the factor composition of microarray intensities changed with different experimental designs in two-color cDNA microarray experiments. By analyzing 12 experimental designs (including 5 general forms), we demonstrate that not all experimental designs are able to provide all possible relative gene expression levels due to the differences in factor composition. Our results indicate that whether an experimental design can provide all possible relative expression levels of all sample pairs for each gene should be the first criterion to be considered in an evaluation of experimental designs for two-color cDNA microarray experiments.     

Towards GAG glycomics: analysis of highly sulfated heparins by MALDI-TOF mass spectrometry

Glycomics is a developing field that provides structural information on complex populations of glycans isolated from tissues, cells and organs. Strategies employing matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) are central to glycomic analysis. Current MALDI-based glycomic strategies are capable of efficiently analyzing glycoprotein and glycosphingolipid glycomes but little attention has been paid to devising glycomic methodologies suited to the analysis of glycosaminoglycan (GAG) polysaccharides which pose special problems for MALDI analysis because of their high level of sulfation and large size. In this paper, we describe MALDI strategies that have been optimized for the analysis of highly sulfated GAG-derived oligosaccharides. A crystalline matrix norharmane, as well as an ionic liquid 1-methylimidazolium alpha-cyano-4-hydroxycinnamate (ImCHCA), have been used for the analysis of heparin di-, tetra-, hexa- and decasaccharides carrying from 2 to 13 sulfate groups. Information about the maximum number of sulfate groups is obtained using the ionic liquid whereas MALDI-TOF/TOF MS/MS experiments using norharmane allowed the determination of the nature of the glycosidic backbone, and more precise information about the presence and the position in the sequence of N-acetylated residues.     

Topology of gene expression networks as revealed by data mining and modeling

MOTIVATION: Interpretation of high-throughput gene expression profiling requires a knowledge of the design principles underlying the networks that sustain cellular machinery. Recently a novel approach based on the study of network topologies has been proposed. This methodology has proven to be useful for the analysis of a variety of biological systems, including metabolic networks, networks of protein-protein interactions, and gene networks that can be derived from gene expression data. In the present paper, we focus on several important issues related to the topology of gene expression networks that have not yet been fully studied. RESULTS: The networks derived from gene expression profiles for both time series experiments in yeast and perturbation experiments in cell lines are studied. We demonstrate that independent from the experimental organism (yeast versus cell lines) and the type of experiment (time courses versus perturbations) the extracted networks have similar topological characteristics suggesting together with the results of other common principles of the structural organization of biological networks. A novel computational model of network growth that reproduces the basic design principles of the observed networks is presented. Advantage of the model is that it provides a general mechanism to generate networks with different types of topology by a variation of a few parameters. We investigate the robustness of the network structure to random damages and to deliberate removal of the most important parts of the system and show a surprising tolerance of gene expression networks to both kinds of disturbance.     

Assessment method for a power analysis to identify differentially expressed pathways

Gene expression data can provide a very rich source of information for elucidating the biological function on the pathway level if the experimental design considers the needs of the statistical analysis methods. The purpose of this paper is to provide a comparative analysis of statistical methods for detecting the differentially expression of pathways (DEP). In contrast to many other studies conducted so far, we use three novel simulation types, producing a more realistic correlation structure than previous simulation methods. This includes also the generation of surrogate data from two large-scale microarray experiments from prostate cancer and ALL. As a result from our comprehensive analysis of 41,004 parameter configurations, we find that each method should only be applied if certain conditions of the data from a pathway are met. Further, we provide method-specific estimates for the optimal sample size for microarray experiments aiming to identify DEP in order to avoid an underpowered design. Our study highlights the sensitivity of the studied methods on the parameters of the system.     

Result-driven strategies for protein identification and quantitation--a way to optimize experimental design and derive reliable results

Uni- or multidimensional microcapillary liquid chromatography (microLC) matrix-assisted laser desorption/ionization (MALDI) tandem mass spectrometry (MS/MS) approaches have gained significant attention for quantifying and identifying proteins in complex biological samples. The off-line coupling of microLC with MS quantitation and MS/MS identification methods makes new result-dependent workflows possible. A relational database is used to store the results from multiple high performance liquid chromatography runs, including information about MALDI plate positions, and both peptide and protein quantitations, and identifications. Unlike electrospray methodology, where all the decisions about which peptide to fragment, must be made during peptide fractionations, in the MALDI experiments the samples are effectively "frozen in time". Therefore, additional MS and MS/MS spectra can be acquired, to promote more accurate quantitation or additional identifications until reliable results are derived that meet experimental design criteria. In the case of what can be designated the expression-dependent workflow, quantitation can be detached from identification and only peak pairs with biological relevant expression changes can be selected for further MS/MS analyses. Alternatively, additional MS/MS data can be acquired to confirm tentative peptide mass fingerprint hits in what is designated a search result-dependent workflow. In the MS data-dependent workflow, the goal is to collect as many meaningful spectra as possible by judiciously adjusting the acquisition parameters based on characteristics of the parent masses. This level of sophistication requires the development of innovative algorithms for these three result-dependent workflows that make MS and MS/MS analysis more efficient and also add confidence to experimental results.     

Computational analysis of microarray data

Microarray experiments are providing unprecedented quantities of genome-wide data on gene-expression patterns. Although this technique has been enthusiastically developed and applied in many biological contexts, the management and analysis of the millions of data points that result from these experiments has received less attention. Sophisticated computational tools are available, but the methods that are used to analyse the data can have a profound influence on the interpretation of the results. A basic understanding of these computational tools is therefore required for optimal experimental design and meaningful data analysis.     

Comparison of Different Three-site Interaction Potentials for Liquid Acetonitrile

Abstract

Computer simulations of liquid acetonitrile at normal room conditions are reported. Both static and dynamic properties are analysed. Special attention is paid to the dielectric properties. A three-site interaction potential has been derived from ab initio calculations on the gas phase dimer and a comparison with different three-site interaction potentials available in the literature is presented. The suitability of three-site models to reproduce the properties of the real liquid is discussed by comparing computer simulation results with experimental data.     

Proteomic biomarker identification for diagnosis of early relapse in ovarian cancer

Ovarian cancer recurs at the rate of 75% within a few months or several years later after therapy. Early recurrence, though responding better to treatment, is difficult to detect. Surface-enhanced laser desorption/ionization time-of-flight (SELDI-TOF) mass spectrometry has showed the potential to accurately identify disease biomarkers to help early diagnosis. A major challenge in the interpretation of SELDI-TOF data is the high dimensionality of the feature space. To tackle this problem, we have developed a multi-step data processing method composed of t-test, binning and backward feature selection. A new algorithm, support vector machine-Markov blanket/recursive feature elimination (SVM-MB/RFE) is presented for the backward feature selection. This method is an integration of minimum weight feature elimination by SVM-RFE and information theory based redundant/irrelevant feature removal by Markov Blanket. Subsequently, SVM was used for classification. We conducted the biomarker selection algorithm on 113 serum samples to identify early relapse from ovarian cancer patients after primary therapy. To validate the performance of the proposed algorithm, experiments were carried out in comparison with several other feature selection and classification algorithms.     

Multiplexed microsphere diagnostic tools in gene expression applications: factors and futures

Microarrays have received significant attention in recent years as scientists have firstly identified factors that can produce reduced confidence in gene expression data obtained on these platforms, and secondly sought to establish laboratory practices and a set of standards by which data are reported with integrity. Microsphere-based assays represent a new generation of diagnostics in this field capable of providing substantial quantitative and qualitative information from gene expression profiling. However, for gene expression profiling, this type of platform is still in the demonstration phase, with issues arising from comparative studies in the literature not yet identified. It is desirable to identify potential parameters that are established as important in controlling the information derived from microsphere-based hybridizations to quantify gene expression. As these evolve, a standard set of parameters will be established that are required to be provided when data are submitted for publication. Here we initiate this process by identifying a number of parameters we have found to be important in microsphere-based assays designed for the quantification of low abundant genes which are variable between studies.     

Analyzing time series gene expression data

MOTIVATION: Time series expression experiments are an increasingly popular method for studying a wide range of biological systems. However, when analyzing these experiments researchers face many new computational challenges. Algorithms that are specifically designed for time series experiments are required so that we can take advantage of their unique features (such as the ability to infer causality from the temporal response pattern) and address the unique problems they raise (e.g. handling the different non-uniform sampling rates). RESULTS: We present a comprehensive review of the current research in time series expression data analysis. We divide the computational challenges into four analysis levels: experimental design, data analysis, pattern recognition and networks. For each of these levels, we discuss computational and biological problems at that level and point out some of the methods that have been proposed to deal with these issues. Many open problems in all these levels are discussed. This review is intended to serve as both, a point of reference for experimental biologists looking for practical solutions for analyzing their data, and a starting point for computer scientists interested in working on the computational problems related to time series expression analysis.     

Bioavailability of phenolic acids

Two large classes of phenolic acids were comprised in this review: benzoic acid derivatives and cinnamic acid derivatives. They have been found to be very extended in fruits and vegetables at different concentrations. For example, hydroxycinnamic acids concentration was higher than that found for hydroxybenzoic acids. Concerning their consumption, hydroxycinnamic acids provide larger contributions to the total polyphenol intake than benzoic acid derivatives or flavonoids. This phenolic acid intake is led by the coffee intake since it has very rich concentrations in hydroxycinnamic acids. Moreover, several experimental and epidemiological studies report the protection of phenolic acids against various degenerative diseases. However, despite all these interesting attributions and even if phenolic acids are the main polyphenols consumed, their bioavailability has not received as attention as that flavonoids. This concept is an essential step to understand the health-promoting properties of phenolic acids and to serve as tool to design in vivo and in vitro experiments to know their biological properties. Therefore, a compilation of bioavailability data of phenolic acids have been presented here paying attention to the two types of phenolic acid bioavailability, direct and indirect derived from the direct phenolic acid and flavonoid consumption, respectively. Then, a new relevant concept which may be named as total bioavailability of phenolic acids includes the direct absorption and metabolism of phenolic acids from food consumption and phenolic acids bioavailability as a result of the cleavage on the main skeleton ring of flavonoids by the gut microflora.     

遗传与基因表达数据的整合—— eQTL的方法及应用

高通量的基因型分析和芯片技术的发展使人们能够进一步研究哪些遗传差异最终影响基因的表达。通过表达数量性状座位(eQTL)作图方法可对基因表达水平的遗传基础进行解析。与传统的QTL分析方法一样, eQTL的主要目标是鉴别表达性状座位所在的染色体区域。但由于表达谱数据成千上万, 而传统的QTL分析方法最多分析几十个性状, 因此需要考虑这类实验设计的特点以及统计分析方法。本文详细介绍了eQTL定位过程及其研究方法, 重点从个体选择、基因芯片实验设计、基因表达数据的获得与标准化、作图方法及结果分析等方面进行了综述, 指出了当前eQTL研究存在的问题和局限性。最后介绍了eQTL研究在估计基因表达遗传率、挖掘候选基因、构建基因调控网络、理解基因间及基因与环境的互作的应用进展。     

Modern Tumor Marker Discovery in Urology: Surface Enhanced Laser Desorption and Ionization (SELDI)

During the last two decades, biomarker research has benefited from the introduction of new proteomic analytical techniques. In this article, we review the application of surface enhanced laser desorption/ionization time-of-flight (SELDI-TOF) mass spectroscopy in urologic cancer research. After reviewing the literature from MEDLINE on proteomics and urologic oncology, we found that SELDI-TOF is an emerging proteomic technology in biomarker discovery that allows for rapid and sensitive analysis of complex protein mixtures. SELDI-TOF is a novel proteomic technology that has the potential to contribute further to the understanding and clinical exploitation of new, clinically relevant biomarkers.     

Phylogenetic information and experimental design in molecular systematics.

Despite the widespread perception that evolutionary inference from molecular sequences is a statistical problem, there has been very little attention paid to questions of experimental design. Previous consideration of this topic has led to little more than an empirical folklore regarding the choice of suitable genes for analysis, and to dispute over the best choice of taxa for inclusion in data sets. I introduce what I believe are new methods that permit the quantification of phylogenetic information in a sequence alignment. The methods use likelihood calculations based on Markov-process models of nucleotide substitution allied with phylogenetic trees, and allow a general approach to optimal experimental design. Two examples are given, illustrating realistic problems in experimental design in molecular phylogenetics and suggesting more general conclusions about the choice of genomic regions, sequence lengths and taxa for evolutionary studies.     

Optimal sampling time selection for parameter estimation in dynamic pathway modeling

Systems Biology is an emerging research area, which considers mathematical representations of inter- and intra-cellular dynamics. Among the many research problems that have been addressed, dynamic modeling of signal transduction pathways has received increasing attention. The usual approach to represent intra-cellular dynamics are nonlinear, usually ordinary, differential equations. The purpose of the models is to test and generate hypothesis of specific pathways and it is therefore required to estimate model parameters from experimental data. The experiments to generate data are complex and expensive, as a consequence of which the time series available are usually rather short, with few if any replicates. Almost certainly, not all variables one would like to include in a model can be measured. Parameter estimation is therefore an important research problem in Systems Biology and the focus of this paper. In particular, we are interested in optimizing the sampling time selection in order to minimize the variance of the parameter estimation error. With few sampling time points feasible, their selection is of practical importance in experimental design. Finally, the theoretical results are supported with an application.     

Microarray data quality analysis: lessons from the AFGC project

Genome-wide expression profiling with DNA microarrays has and will provide a great deal of data to the plant scientific community. However, reliability concerns have required the development data quality tests for common systematic biases. Fortunately, most large-scale systematic biases are detectable and some are correctable by normalization. Technical replication experiments and statistical surveys indicate that these biases vary widely in severity and appearance. As a result, no single normalization or correction method currently available is able to address all the issues. However, careful sequence selection, array design, experimental design and experimental annotation can substantially improve the quality and biological of microarray data. In this review, we discuss these issues with reference to examples from the Arabidopsis Functional Genomics Consortium (AFGC) microarray project.