Bayesian modelling of shared gene function

Predictive and prognostic molecular markers for cancer medicine

Over the last 10 years there has been an explosion of information about the molecular biology of cancer. A challenge in oncology is to translate this information into advances in patient care. While there are well-formed routes for translating new molecular information into drug therapy, the routes for translating new information into sensitive and specific diagnostic, prognostic and predictive tests are still being developed. Similarly, the science of using tumor molecular profiles to select clinical trial participants or to optimize therapy for individual patients is still in its infancy. This review will summarize the current technologies for predicting treatment response and prognosis in cancer medicine, and outline what the future may hold. It will also highlight the potential importance of methods that can integrate molecular, histopathological and clinical information into a synergistic understanding of tumor progression. While these possibilities are without doubt exciting, significant challenges remain if we are to implement them with a strong evidence base in a widely available and cost-effective manner.

The impact of quantitative optimization of hybridization conditions on gene expression analysis

Background

With the growing availability of entire genome sequences, an increasing number of scientists can exploit oligonucleotide microarrays for genome-scale expression studies. While probe-design is a major research area, relatively little work has been reported on the optimization of microarray protocols.

Results

As shown in this study, suboptimal conditions can have considerable impact on biologically relevant observations. For example, deviation from the optimal temperature by one degree Celsius lead to a loss of up to 44% of differentially expressed genes identified. While genes from thousands of Gene Ontology categories were affected, transcription factors and other low-copy-number regulators were disproportionately lost. Calibrated protocols are thus required in order to take full advantage of the large dynamic range of microarrays.

For an objective optimization of protocols we introduce an approach that maximizes the amount of information obtained per experiment. A comparison of two typical samples is sufficient for this calibration. We can ensure, however, that optimization results are independent of the samples and the specific measures used for calibration. Both simulations and spike-in experiments confirmed an unbiased determination of generally optimal experimental conditions.

Conclusions

Well calibrated hybridization conditions are thus easily achieved and necessary for the efficient detection of differential expression. They are essential for the sensitive pro filing of low-copy-number molecules. This is particularly critical for studies of transcription factor expression, or the inference and study of regulatory networks.

Biological assessment of robust noise models in microarray data analysis

Motivation: Although several recently proposed analysis packages for microarray data can cope with heavy-tailed noise, many applications rely on Gaussian assumptions. Gaussian noise models foster computational efficiency. This comes, however, at the expense of increased sensitivity to outlying observations. Assessing potential insufficiencies of Gaussian noise in microarray data analysis is thus important and of general interest.

Results: We propose to this end assessing different noise models on a large number of microarray experiments. The goodness of fit of noise models is quantified by a hierarchical Bayesian analysis of variance model, which predicts normalized expression values as a mixture of a Gaussian density and t-distributions with adjustable degrees of freedom. Inference of differentially expressed genes is taken into consideration at a second mixing level. For attaining far reaching validity, our investigations cover a wide range of analysis platforms and experimental settings. As the most striking result, we find irrespective of the chosen preprocessing and normalization method in all experiments that a heavy-tailed noise model is a better fit than a simple Gaussian. Further investigations revealed that an appropriate choice of noise model has a considerable influence on biological interpretations drawn at the level of inferred genes and gene ontology terms. We conclude from our investigation that neglecting the over dispersed noise in microarray data can mislead scientific discovery and suggest that the convenience of Gaussian-based modelling should be replaced by non-parametric approaches or other methods that account for heavy-tailed noise.

Contact:peter.sykacek@boku.ac.at

Availability:http://bioinf.boku.ac.at/alexp/robmca.html.