Biomedical subjects

Mark Girolami

Publications and source records attributed to Mark Girolami.

5 recordsLinked to original sources

PubMed2006

Accurate prediction of BRCA1 and BRCA2 heterozygous genotype using expression profiling after induced DNA damage.

PURPOSE: In this study, the differential gene expression changes following radiation-induced DNA damage in healthy cells from BRCA1/BRCA1 mutation carriers have been compared with controls using high-density microarray technology. We aimed to establish if BRCA1/BRCA2 mutation carriers could be distinguished from noncarriers based on expression profiling of normal cells. EXPERIMENTAL DESIGN: Short-term primary fibroblast cultures were established from skin biopsies from 10 BRCA1 and 10 BRCA2 mutation carriers and 10 controls, all of whom had previously had breast cancer. The cells were subjected to 15 Gy ionizing irradiation to induce DNA damage. RNA was extracted from all cell cultures, preirradiation and at 1 hour postirradiation. For expression profiling, 15 K spotted cDNA microarrays manufactured by the Cancer Research UK DNA Microarray Facility were used. Statistical feature selection was used with a support vector machine (SVM) classifier to determine the best feature set for predicting BRCA1 or BRCA2 heterozygous genotype. To investigate prediction accuracy, a nonprobabilistic classifier (SVM) and a probabilistic Gaussian process classifier were used. RESULTS: In the task of distinguishing BRCA1 and BRCA2 mutation carriers from noncarriers and from each other following radiation-induced DNA damage, the SVM achieved 90%, and the Gaussian process classifier achieved 100% accuracy. This effect could not be achieved without irradiation. In addition, the SVM identified a set of BRCA genotype predictor genes. CONCLUSIONS: We conclude that after irradiation-induced DNA damage, BRCA1 and BRCA2 mutation carrier cells have a distinctive expression phenotype, and this may have a future role in predicting genotypes, with application to clinical detection and classification of mutations.

Adult↗

PubMed2005

A Bayesian regression approach to the inference of regulatory networks from gene expression data.

MOTIVATION: There is currently much interest in reverse-engineering regulatory relationships between genes from microarray expression data. We propose a new algorithmic method for inferring such interactions between genes using data from gene knockout experiments. The algorithm we use is the Sparse Bayesian regression algorithm of Tipping and Faul. This method is highly suited to this problem as it does not require the data to be discretized, overcomes the need for an explicit topology search and, most importantly, requires no heuristic thresholding of the discovered connections. RESULTS: Using simulated expression data, we are able to show that this algorithm outperforms a recently published correlation-based approach. Crucially, it does this without the need to set any ad hoc threshold on possible connections.

Algorithms↗

PubMed2004

Biologically valid linear factor models of gene expression.

MOTIVATION: The identification of physiological processes underlying and generating the expression pattern observed in microarray experiments is a major challenge. Principal component analysis (PCA) is a linear multivariate statistical method that is regularly employed for that purpose as it provides a reduced-dimensional representation for subsequent study of possible biological processes responding to the particular experimental conditions. Making explicit the data assumptions underlying PCA highlights their lack of biological validity thus making biological interpretation of the principal components problematic. A microarray data representation which enables clear biological interpretation is a desirable analysis tool. RESULTS: We address this issue by employing the probabilistic interpretation of PCA and proposing alternative linear factor models which are based on refined biological assumptions. A practical study on two well-understood microarray datasets highlights the weakness of PCA and the greater biological interpretability of the linear models we have developed.

Algorithms↗

PubMed2002

Orthogonal series density estimation and the kernel eigenvalue problem.

Kernel principal component analysis has been introduced as a method of extracting a set of orthonormal nonlinear features from multivariate data, and many impressive applications are being reported within the literature. This article presents the view that the eigenvalue decomposition of a kernel matrix can also provide the discrete expansion coefficients required for a nonparametric orthogonal series density estimator. In addition to providing novel insights into nonparametric density estimation, this article provides an intuitively appealing interpretation for the nonlinear features extracted from data using kernel principal component analysis.

Journal Article↗

PubMedDate not supplied

The latent process decomposition of cDNA microarray data sets.

We present a new computational technique (a software implementation, data sets, and supplementary information are available at http://www.enm.bris.ac.uk/lpd/) which enables the probabilistic analysis of cDNA microarray data and we demonstrate its effectiveness in identifying features of biomedical importance. A hierarchical Bayesian model, called Latent Process Decomposition (LPD), is introduced in which each sample in the data set is represented as a combinatorial mixture over a finite set of latent processes, which are expected to correspond to biological processes. Parameters in the model are estimated using efficient variational methods. This type of probabilistic model is most appropriate for the interpretation of measurement data generated by cDNA microarray technology. For determining informative substructure in such data sets, the proposed model has several important advantages over the standard use of dendrograms. First, the ability to objectively assess the optimal number of sample clusters. Second, the ability to represent samples and gene expression levels using a common set of latent variables (dendrograms cluster samples and gene expression values separately which amounts to two distinct reduced space representations). Third, in constrast to standard cluster models, observations are not assigned to a single cluster and, thus, for example, gene expression levels are modeled via combinations of the latent processes identified by the algorithm. We show this new method compares favorably with alternative cluster analysis methods. To illustrate its potential, we apply the proposed technique to several microarray data sets for cancer. For these data sets it successfully decomposes the data into known subtypes and indicates possible further taxonomic subdivision in addition to highlighting, in a wholly unsupervised manner, the importance of certain genes which are known to be medically significant. To illustrate its wider applicability, we also illustrate its performance on a microarray data set for yeast.

Algorithms↗