Imaging Genomics of Glioma Revisited: Analytic Methods to Understand Spatial and Temporal Heterogeneity

An improved understanding of the cellular and molecular biologic processes responsible for brain tumor development, growth, and resistance to therapy is fundamental to improving clinical outcomes. Imaging genomics is the study of the relationships between microscopic, genetic, and molecular biologic features and macroscopic imaging features. Imaging genomics is beginning to shift clinical paradigms for diagnosing and treating brain tumors. This article provides an overview of imaging genomics in gliomas, in which imaging data including hallmarks such as IDH-mutation, MGMT methylation, and EGFR-mutation status can provide critical insights into the pretreatment and posttreatment stages. This article will accomplish the following: 1) review the methods used in imaging genomics, including visual analysis, quantitative analysis, and radiomics analysis; 2) recommend suitable analytic methods for imaging genomics according to biologic characteristics; 3) discuss the clinical applicability of imaging genomics; and 4) introduce subregional tumor habitat analysis with the goal of guiding future radiogenetics research endeavors toward translation into critically needed clinical applications.

Visual Analytics for Hypothesis-Driven Exploration in Computational Pathology

Recent advances in computational and algorithmic power are evolving the field of medical imaging rapidly. In cancer research, many new directions are sought to characterize patients with additional imaging features derived from radiology and pathology images. The emerging field of Computational Pathology targets the high-throughput extraction and analysis of the spatial distribution of cells from digital histopathology images. The associated morphological and architectural features allow researchers to quantify and characterize new imaging biomarkers for cancer diagnosis, prognosis, and treatment decisions. However, while the image feature space grows, exploration and analysis become more difficult and ineffective. There is a need for dedicated interfaces for interactive data manipulation and visual analysis of computational pathology and clinical data. For this purpose, we present IIComPath, a visual analytics approach that enables clinical researchers to formulate hypotheses and create computational pathology pipelines involving cohort construction, spatial analysis of image-derived features, and cohort analysis. We demonstrate our approach through use cases that investigate the prognostic value of current diagnostic features and new computational pathology biomarkers.

A Latent Gaussian Copula Model for Mixed Data Analysis in Brain Imaging Genetics

Recent advances in imaging genetics make it possible to combine different types of data including medical images like functional magnetic resonance imaging (fMRI) and genetic data like single nucleotide polymorphisms (SNPs) for comprehensive diagnosis of mental disorders. Understanding complex interactions among these heterogeneous data may give rise to a new perspective, while at the same time demand statistical models for their integration. Various graphical models have been proposed for the study of interaction or association networks with continuous, binary, and count data as well as the mixture of them. However, limited efforts have been made for the multinomial case, for instance, SNP data. Our goal is therefore to fill the void by developing a graphical model for the integration of fMRI image and SNP data, which can provide deeper understanding of the unknown neurogenetic mechanism. In this article, we propose a latent Gaussian copula model for mixed data containing multinomial components. We assume that the discrete variable is obtained by discretizing a latent (unobserved) continuous variable and then create a semi-rank based estimator of the graph structure. The simulation results demonstrate that the proposed latent correlation has more steady and accurate performance than several existing methods in detecting graph structure. When applying to a real schizophrenia data consisting of SNP array and fMRI image collected by the Mind Clinical Imaging Consortium (MCIC), the proposed method reveals a set of distinct SNP-brain associations, which are verified to be biologically significant. The proposed model is statistically promising in handling mixed types of data including multinomial components, which can find widespread applications. To promote reproducible research, the R code is available at https://github.com/Aiying0512/LGCM.

Identify Consistent Cross-Modality Imaging Genetic Patterns via Discriminant Sparse Canonical Correlation Analysis

Sparse canonical correlation analysis (SCCA) is a bi-multivariate technique used in imaging genetics to identify complex multi-SNP-multi-QT associations. However, the traditional SCCA algorithm has been designed to seek a linear correlation between the SNP genotype and brain imaging phenotype, ignoring the discriminant similarity information between within-class subjects in brain imaging genetics association analysis. In addition, multi-modality brain imaging phenotypes are extracted from different perspectives and imaging markers from the same region consistently showing up in multimodalities may provide more insights for the mechanistic understanding of diseases. In this paper, a novel multi-modality discriminant SCCA algorithm (MD-SCCA) is proposed to overcome these limitations as well as to improve learning results by incorporating valuable discriminant similarity information into the SCCA algorithm. Specifically, we first extract the discriminant similarity information between within-class subjects by the sparse representation. Second, the discriminant similarity information is enforced within SCCA to construct a discriminant SCCA algorithm (D-SCCA). At last, the MD-SCCA algorithm is adopted to fully explore the relationships among different modalities of different subjects. In experiments, both synthetic dataset and real data from the Alzheimer's Disease Neuroimaging Initiative database are used to test the performance of our algorithm. The empirical results have demonstrated that the proposed algorithm not only produces improved cross-validation performances but also identifies consistent cross-modality imaging genetic biomarkers.

Tissue clearing techniques for three-dimensional optical imaging of intact human prostate and correlations with multi-parametric MRI

BACKGROUND

Tissue clearing technologies have enabled remarkable advancements for in situ characterization of tissues and exploration of the three-dimensional (3D) relationships between cells, however, these studies have predominantly been performed in non-human tissues and correlative assessment with clinical imaging has yet to be explored. We sought to evaluate the feasibility of tissue clearing technologies for 3D imaging of intact human prostate and the mapping of structurally and molecularly preserved pathology data with multi-parametric volumetric MR imaging (mpMRI).

METHODS

Whole-mount prostates were processed with either hydrogel-based CLARITY or solvent-based iDISCO. The samples were stained with a nuclear dye or fluorescently labeled with antibodies against androgen receptor, alpha-methylacyl coenzyme-A racemase, or p63, and then imaged with 3D confocal microscopy. The apparent diffusion coefficient and K^trans maps were computed from preoperative mpMRI.

RESULTS

Quantitative analysis of cleared normal and tumor prostate tissue volumes displayed differences in 3D tissue architecture, marker-specific cell staining, and cell densities that were significantly correlated with mpMRI measurements in this initial, pilot cohort.

CONCLUSIONS

3D imaging of human prostate volumes following tissue clearing is a feasible technique for quantitative radiology-pathology correlation analysis with mpMRI and provides an opportunity to explore functional relationships between cellular structures and cross-sectional clinical imaging.

Multi-Task Sparse Canonical Correlation Analysis with Application to Multi-Modal Brain Imaging Genetics

Brain imaging genetics studies the genetic basis of brain structures and functionalities via integrating genotypic data such as single nucleotide polymorphisms (SNPs) and imaging quantitative traits (QTs). In this area, both multi-task learning (MTL) and sparse canonical correlation analysis (SCCA) methods are widely used since they are superior to those independent and pairwise univariate analysis. MTL methods generally incorporate a few of QTs and could not select features from multiple QTs; while SCCA methods typically employ one modality of QTs to study its association with SNPs. Both MTL and SCCA are computational expensive as the number of SNPs increases. In this paper, we propose a novel multi-task SCCA (MTSCCA) method to identify bi-multivariate associations between SNPs and multi-modal imaging QTs. MTSCCA could make use of the complementary information carried by different imaging modalities. MTSCCA enforces sparsity at the group level via the G2,1-norm, and jointly selects features across multiple tasks for SNPs and QTs via the l2,1-norm. A fast optimization algorithm is proposed using the grouping information of SNPs. Compared with conventional SCCA methods, MTSCCA obtains better correlation coefficients and canonical weights patterns. In addition, MTSCCA runs very fast and easy-to-implement, indicating its potential power in genome-wide brain-wide imaging genetics.

Integration of Imaging (epi)Genomics Data for the Study of Schizophrenia Using Group Sparse Joint Nonnegative Matrix Factorization

Schizophrenia (SZ) is a complex disease. Single nucleotide polymorphism (SNP), brain activity measured by functional magnetic resonance imaging (fMRI) and DNA methylation are all important biomarkers that can be used for the study of SZ. To our knowledge, there has been little effort to combine these three datasets together. In this study, we propose a group sparse joint nonnegative matrix factorization (GSJNMF) model to integrate SNP, fMRI, and DNA methylation for the identification of multi-dimensional modules associated with SZ, which can be used to study regulatory mechanisms underlying SZ at multiple levels. The proposed GSJNMF model projects multiple types of data onto a common feature space, in which heterogeneous variables with large coefficients on the same projected bases are used to identify multi-dimensional modules. We also incorporate group structure information available from each dataset. The genomic factors in such modules have significant correlations or functional associations with several brain activities. At the end, we have applied the method to the analysis of real data collected from the Mind Clinical Imaging Consortium (MCIC) for the study of SZ and identified significant biomarkers. These biomarkers were further used to discover genes and corresponding brain regions, which were confirmed to be significantly associated with SZ.

MRI-based radiogenomics analysis for predicting genetic alterations in oncogenic signalling pathways in invasive breast carcinoma

AIM

To investigate the effect of radiomics in the assessment of alterations in canonical cancer pathways in breast cancer.

MATERIALS AND METHODS

Eighty-eight biopsy-proven breast cancer cases were included in the present study. Radiomics features were extracted from T1-weighted sagittal dynamic contrast-enhanced magnetic resonance imaging (MRI) images. Radiomics signatures were developed to predict genetic alterations in the cell cycle, Myc, PI3K, RTK/RAS, and p53 signalling pathways by using hypothesis testing combined with least absolute shrinkage and selection operator (LASSO) regression analysis. The predictive powers of the models were examined by the area under the curve (AUC) of the receiver operating characteristic curve.

RESULTS

A total of 5,234 radiomics features were obtained from MRI images based on the tumour region of interest. Hypothesis tests screened 250, 229, 156, 785, and 319 radiomics features that were differentially displayed between cell cycle, Myc, PI3K, RTK/RAS, and p53 alterations and no alteration status. According to the LASSO algorithm, 11, 12, 12, 15, and 13 features were identified for the construction of the radiomics signatures to predict cell cycle, Myc, PI3K, RTK/RAS, and p53 alterations, with AUC values of 0.933, 0.926, 0.956, 0.940, and 0.886, respectively. The cell cycle radiomics score correlated closely with the RTK/RAS and p53 radiomics scores. These signatures were also dysregulated in patients with different oestrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2 statuses.

CONCLUSION

MRI-based radiogenomics analysis exhibits excellent performance in predicting genetic pathways alterations, thus providing a novel approach for non-invasively obtaining genetic-level molecular characteristics of tumours.