Magnetic resonance image analysis by information theoretic criteria and stochastic site models

BACOM2.0 facilitates absolute normalization and quantification of somatic copy number alterations in heterogeneous tumor

Most published copy number datasets on solid tumors were obtained from specimens comprised of mixed cell populations, for which the varying tumor-stroma proportions are unknown or unreported. The inability to correct for signal mixing represents a major limitation on the use of these datasets for subsequent analyses, such as discerning deletion types or detecting driver aberrations. We describe the BACOM2.0 method with enhanced accuracy and functionality to normalize copy number signals, detect deletion types, estimate tumor purity, quantify true copy numbers, and calculate average-ploidy value. While BACOM has been validated and used with promising results, subsequent BACOM analysis of the TCGA ovarian cancer dataset found that the estimated average tumor purity was lower than expected. In this report, we first show that this lowered estimate of tumor purity is the combined result of imprecise signal normalization and parameter estimation. Then, we describe effective allele-specific absolute normalization and quantification methods that can enhance BACOM applications in many biological contexts while in the presence of various confounders. Finally, we discuss the advantages of BACOM in relation to alternative approaches. Here we detail this revised computational approach, BACOM2.0, and validate its performance in real and simulated datasets.

Unsupervised deconvolution of dynamic imaging reveals intratumor vascular heterogeneity and repopulation dynamics

With the existence of biologically distinctive malignant cells originated within the same tumor, intratumor functional heterogeneity is present in many cancers and is often manifested by the intermingled vascular compartments with distinct pharmacokinetics. However, intratumor vascular heterogeneity cannot be resolved directly by most in vivo dynamic imaging. We developed multi-tissue compartment modeling (MTCM), a completely unsupervised method of deconvoluting dynamic imaging series from heterogeneous tumors that can improve vascular characterization in many biological contexts. Applying MTCM to dynamic contrast-enhanced magnetic resonance imaging of breast cancers revealed characteristic intratumor vascular heterogeneity and therapeutic responses that were otherwise undetectable. MTCM is readily applicable to other dynamic imaging modalities for studying intratumor functional and phenotypic heterogeneity, together with a variety of foreseeable applications in the clinic.

Integration of Network Biology and Imaging to Study Cancer Phenotypes and Responses

Ever growing "omics" data and continuously accumulated biological knowledge provide an unprecedented opportunity to identify molecular biomarkers and their interactions that are responsible for cancer phenotypes that can be accurately defined by clinical measurements such as in vivo imaging. Since signaling or regulatory networks are dynamic and context-specific, systematic efforts to characterize such structural alterations must effectively distinguish significant network rewiring from random background fluctuations. Here we introduced a novel integration of network biology and imaging to study cancer phenotypes and responses to treatments at the molecular systems level. Specifically, Differential Dependence Network (DDN) analysis was used to detect statistically significant topological rewiring in molecular networks between two phenotypic conditions, and in vivo Magnetic Resonance Imaging (MRI) was used to more accurately define phenotypic sample groups for such differential analysis. We applied DDN to analyze two distinct phenotypic groups of breast cancer and study how genomic instability affects the molecular network topologies in high-grade ovarian cancer. Further, FDA-approved arsenic trioxide (ATO) and the ND2-SmoA1 mouse model of Medulloblastoma (MB) were used to extend our analyses of combined MRI and Reverse Phase Protein Microarray (RPMA) data to assess tumor responses to ATO and to uncover the complexity of therapeutic molecular biology.

Tissue-specific compartmental analysis for dynamic contrast-enhanced MR imaging of complex tumors

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides a noninvasive method for evaluating tumor vasculature patterns based on contrast accumulation and washout. However, due to limited imaging resolution and tumor tissue heterogeneity, tracer concentrations at many pixels often represent a mixture of more than one distinct compartment. This pixel-wise partial volume effect (PVE) would have profound impact on the accuracy of pharmacokinetics studies using existing compartmental modeling (CM) methods. We, therefore, propose a convex analysis of mixtures (CAM) algorithm to explicitly mitigate PVE by expressing the kinetics in each pixel as a nonnegative combination of underlying compartments and subsequently identifying pure volume pixels at the corners of the clustered pixel time series scatter plot simplex. The algorithm is supported theoretically by a well-grounded mathematical framework and practically by plug-in noise filtering and normalization preprocessing. We demonstrate the principle and feasibility of the CAM-CM approach on realistic synthetic data involving two functional tissue compartments, and compare the accuracy of parameter estimates obtained with and without PVE elimination using CAM or other relevant techniques. Experimental results show that CAM-CM achieves a significant improvement in the accuracy of kinetic parameter estimation. We apply the algorithm to real DCE-MRI breast cancer data and observe improved pharmacokinetic parameter estimation, separating tumor tissue into regions with differential tracer kinetics on a pixel-by-pixel basis and revealing biologically plausible tumor tissue heterogeneity patterns. This method combines the advantages of multivariate clustering, convex geometry analysis, and compartmental modeling approaches. The open-source MATLAB software of CAM-CM is publicly available from the Web.

Genetic algorithms for finite mixture model based voxel classification in neuroimaging

Finite mixture models (FMMs) are an indispensable tool for unsupervised classification in brain imaging. Fitting an FMM to the data leads to a complex optimization problem. This optimization problem is difficult to solve by standard local optimization methods, such as the expectation-maximization (EM) algorithm, if a principled initialization is not available. In this paper, we propose a new global optimization algorithm for the FMM parameter estimation problem, which is based on real coded genetic algorithms. Our specific contributions are two-fold: 1) we propose to use blended crossover in order to reduce the premature convergence problem to its minimum and 2) we introduce a completely new permutation operator specifically meant for the FMM parameter estimation. In addition to improving the optimization results, the permutation operator allows for imposing biologically meaningful constraints to the FMM parameter values. We also introduce a hybrid of the genetic algorithm and the EM algorithm for efficient solution of multidimensional FMM fitting problems. We compare our algorithm to the self-annealing EM-algorithm and a standard real coded genetic algorithm with the voxel classification tasks within the brain imaging. The algorithms are tested on synthetic data as well as real three-dimensional image data from human magnetic resonance imaging, positron emission tomography, and mouse brain MRI. The tissue classification results by our method are shown to be consistently more reliable and accurate than with the competing parameter estimation methods.

An adaptive markov model-based method to cluster validation in image segmentation

The number of class should be detected as part of the parameter estimation procedure prior to image segmentation for segmentation algorithms. It is very important in theory and application for estimating the class number correctly. In this paper, an adaptive total energy criterion (ATEC) to cluster validation is proposed based on the Markov random field (MRF) in the image segmentation. The criterion is composed of two parts: one part is inner-energy, which describes the difference of data in the same class; another is inter-class energy, which describes the edge information. The correct class number can be obtained by minimizing the ATEC. The parameters are estimated by expectation maximum (EM) algorithm and maximum psedu-likelihood (MPL) algorithm. The complex computation is optimized by the mixture of simulated algorithm (SA) and iterated conditional mode (ICM). The experiments show that the class number can be automatically detected by adjusting the hyper-parameter in MRF. As a by-product, the segmentation can be obtained by the maximum a posteriori (MAP).

Bayesian image segmentation using local iso-intensity structural orientation

Image segmentation is a fundamental problem in early computer vision. In segmentation of flat shaded, nontextured objects in real-world images, objects are usually assumed to be piecewise homogeneous. This assumption, however, is not always valid with images such as medical images. As a result, any techniques based on this assumption may produce less-than-satisfactory image segmentation. In this work, we relax the piecewise homogeneous assumption. By assuming that the intensity nonuniformity is smooth in the imaged objects, a novel algorithm that exploits the coherence in the intensity profile to segment objects is proposed. The algorithm uses a novel smoothness prior to improve the quality of image segmentation. The formulation of the prior is based on the coherence of the local structural orientation in the image. The segmentation process is performed in a Bayesian framework. Local structural orientation estimation is obtained with an orientation tensor. Comparisons between the conventional Hessian matrix and the orientation tensor have been conducted. The experimental results on the synthetic images and the real-world images have indicated that our novel segmentation algorithm produces better segmentations than both the global thresholding with the maximum likelihood estimation and the algorithm with the multilevel logistic MRF model.

Fast and robust parameter estimation for statistical partial volume models in brain MRI

Due to the finite spatial resolution of imaging devices, a single voxel in a medical image may be composed of mixture of tissue types, an effect known as partial volume effect (PVE). Partial volume estimation, that is, the estimation of the amount of each tissue type within each voxel, has received considerable interest in recent years. Much of this work has been focused on the mixel model, a statistical model of PVE. We propose a novel trimmed minimum covariance determinant (TMCD) method for the estimation of the parameters of the mixel PVE model. In this method, each voxel is first labeled according to the most dominant tissue type. Voxels that are prone to PVE are removed from this labeled set, following which robust location estimators with high breakdown points are used to estimate the mean and the covariance of each tissue class. Comparisons between different methods for parameter estimation based on classified images as well as expectation--maximization-like (EM-like) procedure for simultaneous parameter and partial volume estimation are reported. The robust estimators based on a pruned classification as presented here are shown to perform well even if the initial classification is of poor quality. The results obtained are comparable to those obtained using the EM-like procedure, but require considerably less computation time. Segmentation results of real data based on partial volume estimation are also reported. In addition to considering the parameter estimation problem, we discuss differences between different approximations to the complete mixel model. In summary, the proposed TMCD method allows for the accurate, robust, and efficient estimation of partial volume model parameters, which is crucial to a variety of brain MRI data analysis procedures such as the accurate estimation of tissue volumes and the accurate delineation of the cortical surface.

Electrotonically mediated oscillatory patterns in neuronal ensembles: an in vitro voltage-dependent dye-imaging study in the inferior olive

Spatiotemporal profiles of ensemble subthreshold neuronal oscillation were studied in brainstem slices using high-speed voltage-sensitive dye imaging. After local electrical stimuli, the overall voltage profile demonstrated coherent oscillatory waves that spread over the inferior olive (IO). These oscillations were also observed in concurrently obtained intracellular recordings from IO neurons. Over the first few seconds after the stimuli, the optically recorded oscillations clustered into coherent groups comprising hundreds of neurons. Statistical analysis of the spatial profiles of these clusters revealed size fluctuation around stable core regions that were surrounded by a rim the diameter of which varied in time during the oscillation period. The neuronal ensemble oscillations were calcium derived and had an average frequency range of 1-7 Hz. This rhythmic response demonstrated a different spatiotemporal distribution in the presence of picrotoxin, which induced the merging of neuronal clusters into larger areas of coherent activity. The possibility that such clustering is a consequence of intrinsic oscillations in ensembles of coupled neurons was tested using mathematical modeling.