Joint and individual analysis of breast cancer histologic images and genomic covariates

A key challenge in modern data analysis is understanding connections between complex and differing modalities of data. For example, two of the main approaches to the study of breast cancer are histopathology (analyzing visual characteristics of tumors) and genetics. While histopathology is the gold standard for diagnostics and there have been many recent breakthroughs in genetics, there is little overlap between these two fields. We aim to bridge this gap by developing methods based on Angle-based Joint and Individual Variation Explained (AJIVE) to directly explore similarities and differences between these two modalities. Our approach exploits Convolutional Neural Networks (CNNs) as a powerful, automatic method for image feature extraction to address some of the challenges presented by statistical analysis of histopathology image data. CNNs raise issues of interpretability that we address by developing novel methods to explore visual modes of variation captured by statistical algorithms (e.g. PCA or AJIVE) applied to CNN features. Our results provide many interpretable connections and contrasts between histopathology and genetics

Passing to the Limit in a Wasserstein Gradient Flow: From Diffusion to Reaction

We study a singular-limit problem arising in the modelling of chemical reactions. At finite {\epsilon} > 0, the system is described by a Fokker-Planck convection-diffusion equation with a double-well convection potential. This potential is scaled by 1/{\epsilon}, and in the limit {\epsilon} -> 0, the solution concentrates onto the two wells, resulting into a limiting system that is a pair of ordinary differential equations for the density at the two wells. This convergence has been proved in Peletier, Savar\'e, and Veneroni, SIAM Journal on Mathematical Analysis, 42(4):1805-1825, 2010, using the linear structure of the equation. In this paper we re-prove the result by using solely the Wasserstein gradient-flow structure of the system. In particular we make no use of the linearity, nor of the fact that it is a second-order system. The first key step in this approach is a reformulation of the equation as the minimization of an action functional that captures the property of being a curve of maximal slope in an integrated form. The second important step is a rescaling of space. Using only the Wasserstein gradient-flow structure, we prove that the sequence of rescaled solutions is pre-compact in an appropriate topology. We then prove a Gamma-convergence result for the functional in this topology, and we identify the limiting functional and the differential equation that it represents. A consequence of these results is that solutions of the {\epsilon}-problem converge to a solution of the limiting problem.Comment: Added two sections, corrected minor typos, updated reference

BrainPrint: Identifying Subjects by Their Brain

Introducing BrainPrint, a compact and discriminative representation of anatomical structures in the brain. BrainPrint captures shape information of an ensemble of cortical and subcortical structures by solving the 2D and 3D Laplace-Beltrami operator on triangular (boundary) and tetrahedral (volumetric) meshes. We derive a robust classifier for this representation that identifies the subject in a new scan, based on a database of brain scans. In an example dataset containing over 3000 MRI scans, we show that BrainPrint captures unique information about the subject’s anatomy and permits to correctly classify a scan with an accuracy of over 99.8%. All processing steps for obtaining the compact representation are fully automated making this processing framework particularly attractive for handling large datasets.Alexander von Humboldt-StiftungAthinoula A. Martinos Center for Biomedical Imaging (P41-RR014075)Athinoula A. Martinos Center for Biomedical Imaging (P41-EB015896)National Alliance for Medical Image Computing (U.S.) (U54-EB005149)Neuroimaging Analysis Center (U.S.) (P41-EB015902

The visual perception of natural motion: abnormal task-related neural activity in DYT1 dystonia

Although primary dystonia is defined by its characteristic motor manifestations, non-motor signs and symptoms have increasingly been recognized in this disorder. Recent neuroimaging studies have related the motor features of primary dystonia to connectivity changes in cerebello-thalamo-cortical pathways. It is not known, however, whether the non-motor manifestations of the disorder are associated with similar circuit abnormalities. To explore this possibility, we used functional magnetic resonance imaging to study primary dystonia and healthy volunteer subjects while they performed a motion perception task in which elliptical target trajectories were visually tracked on a computer screen. Prior functional magnetic resonance imaging studies of healthy subjects performing this task have revealed selective activation of motor regions during the perception of \u27natural\u27 versus \u27unnatural\u27 motion (defined respectively as trajectories with kinematic properties that either comply with or violate the two-thirds power law of motion). Several regions with significant connectivity changes in primary dystonia were situated in proximity to normal motion perception pathways, suggesting that abnormalities of these circuits may also be present in this disorder. To determine whether activation responses to natural versus unnatural motion in primary dystonia differ from normal, we used functional magnetic resonance imaging to study 10 DYT1 dystonia and 10 healthy control subjects at rest and during the perception of \u27natural\u27 and \u27unnatural\u27 motion. Both groups exhibited significant activation changes across perceptual conditions in the cerebellum, pons, and subthalamic nucleus. The two groups differed, however, in their responses to \u27natural\u27 versus \u27unnatural\u27 motion in these regions. In healthy subjects, regional activation was greater during the perception of natural (versus unnatural) motion (P \u3c 0.05). By contrast, in DYT1 dystonia subjects, activation was relatively greater during the perception of unnatural (versus natural) motion (P \u3c 0.01). To explore the microstructural basis for these functional changes, the regions with significant interaction effects (i.e. those with group differences in activation across perceptual conditions) were used as seeds for tractographic analysis of diffusion tensor imaging scans acquired in the same subjects. Fibre pathways specifically connecting each of the significant functional magnetic resonance imaging clusters to the cerebellum were reconstructed. Of the various reconstructed pathways that were analysed, the ponto-cerebellar projection alone differed between groups, with reduced fibre integrity in dystonia (P \u3c 0.001). In aggregate, the findings suggest that the normal pattern of brain activation in response to motion perception is disrupted in DYT1 dystonia. Thus, it is unlikely that the circuit changes that underlie this disorder are limited to primary sensorimotor pathways

Simulation-Based Joint Estimation of Body Deformation and Elasticity Parameters for Medical Image Analysis

Estimation of tissue stiffness is an important means of noninvasive cancer detection. Existing elasticity reconstruction methods usually depend on a dense displacement field (inferred from ultrasound or MR images) and known external forces. Many imaging modalities, however, cannot provide details within an organ and therefore cannot provide such a displacement field. Furthermore, force exertion and measurement can be difficult for some internal organs, making boundary forces another missing parameter. We propose a general method for estimating elasticity and boundary forces automatically using an iterative optimization framework, given the desired (target) output surface. During the optimization, the input model is deformed by the simulator, and an objective function based on the distance between the deformed surface and the target surface is minimized numerically. The optimization framework does not depend on a particular simulation method and is therefore suitable for different physical models. We show a positive correlation between clinical prostate cancer stage (a clinical measure of severity) and the recovered elasticity of the organ. Since the surface correspondence is established, our method also provides a non-rigid image registration, where the quality of the deformation fields is guaranteed, as they are computed using a physics-based simulation

Restoration of DWI Data Using a Rician LMMSE Estimator

This paper introduces and analyzes a linear minimum mean square error (LMMSE) estimator using a Rician noise model and its recursive version (RLMMSE) for the restoration of diffusion weighted images. A method to estimate the noise level based on local estimations of mean or variance is used to automatically parametrize the estimator. The restoration performance is evaluated using quality indexes and compared to alternative estimation schemes. The overall scheme is simple, robust, fast, and improves estimations. Filtering diffusion weighted magnetic resonance imaging (DW-MRI) with the proposed methodology leads to more accurate tensor estimations. Real and synthetic datasets are analyzed

Using the MitoB method to assess levels of reactive oxygen species in ecological studies of oxidative stress

In recent years evolutionary ecologists have become increasingly interested in the effects of reactive oxygen species (ROS) on the life-histories of animals. ROS levels have mostly been inferred indirectly due to the limitations of estimating ROS from in vitro methods. However, measuring ROS (hydrogen peroxide, H2O2) content in vivo is now possible using the MitoB probe. Here, we extend and refine the MitoB method to make it suitable for ecological studies of oxidative stress using the brown trout Salmo trutta as model. The MitoB method allows an evaluation of H2O2 levels in living organisms over a timescale from hours to days. The method is flexible with regard to the duration of exposure and initial concentration of the MitoB probe, and there is no transfer of the MitoB probe between fish. H2O2 levels were consistent across subsamples of the same liver but differed between muscle subsamples and between tissues of the same animal. The MitoB method provides a convenient method for measuring ROS levels in living animals over a significant period of time. Given its wide range of possible applications, it opens the opportunity to study the role of ROS in mediating life history trade-offs in ecological settings