Convex Analysis of Mixtures for Separating Non-negative Well-grounded Sources

Comprehensive evaluation of deconvolution methods for human brain gene expression

Transcriptome deconvolution aims to estimate the cellular composition of an RNA sample from its gene expression data, which in turn can be used to correct for composition differences across samples. The human brain is unique in its transcriptomic diversity, and comprises a complex mixture of cell-types, including transcriptionally similar subtypes of neurons. Here, we carry out a comprehensive evaluation of deconvolution methods for human brain transcriptome data, and assess the tissue-specificity of our key observations by comparison with human pancreas and heart. We evaluate eight transcriptome deconvolution approaches and nine cell-type signatures, testing the accuracy of deconvolution using in silico mixtures of single-cell RNA-seq data, RNA mixtures, as well as nearly 2000 human brain samples. Our results identify the main factors that drive deconvolution accuracy for brain data, and highlight the importance of biological factors influencing cell-type signatures, such as brain region and in vitro cell culturing.

Transcriptome deconvolution aims to estimate cellular composition based on gene expression data. Here the authors evaluate deconvolution methods for human brain transcriptome and conclude that partial deconvolution algorithms work best, but that appropriate cell-type signatures are also important.

Radiomic Analysis of Pharmacokinetic Heterogeneity Within Tumor Based on the Unsupervised Decomposition of Dynamic Contrast-Enhanced MRI for Predicting Histological Characteristics of Breast Cancer

BACKGROUND

Breast tumor heterogeneity is associated with histological characteristics. However, pharmacokinetic (PK) heterogeneity within tumor might merit further exploration.

PURPOSE

To enhance the predictive power of molecular subtypes, Ki-67, and tumor grade by analyzing PK heterogeneity within tumor based on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

STUDY TYPE

Retrospective.

POPULATION

Two hundred and eight biopsy-proven breast cancer patients, randomly divided into a training cohort (N = 144) and a testing cohort (N = 64).

FIELD STRENGTH/SEQUENCE

T₁ -weighted DCE-MRI at 3.0 T.

ASSESSMENT

A convex analysis of mixtures-compartmental modeling decomposition method was used to estimate the PK parameter (i.e., the volume transfer constant K^trans ) in tumor subregions with distinct physiological kinetic patterns, including fast-flow kinetics, slow-flow kinetics, and plasma input. Radiomic features based on the PK parameter were calculated from each tumor subregion.

STATISTICAL TESTS

The training cohort was used to build random forest classifiers based on the optimal features determined by the 5-fold cross-validation method. The performance was assessed on the testing cohort using the area under the receiver operating characteristic curve (AUC). The AUCs derived from the tumor subregion-based PK parameter were compared with those of the original images of the entire tumor using the DeLong test. A P-value of <0.05 was considered statistically significant.

RESULTS

The tumor subregion-based PK parameter, which yielded the highest AUCs of 0.8782, 0.7568, 0.7019, 0.7963, 0.8080, and 0.7375 for luminal A, luminal B, basal-like, human epidermal growth factor receptor 2, Ki-67, and tumor grade, respectively, obtained better diagnostic performance than the original images in the entire tumor (highest AUCs = 0.8612, 0.6191, 0.5593, 0.7704, 0.7494, and 0.6261, respectively). In particular, statistically significant improvement in the diagnostic performance was obtained for luminal B.

DATA CONCLUSION

Radiomic analysis of PK heterogeneity within tumor can enhance the predictive performance of radiomic models compared with that of the entire tumor.

LEVEL OF EVIDENCE

4 TECHNICAL EFFICACY STAGE: 3.

Molecular characterization of projection neuron subtypes in the mouse olfactory bulb

Projection neurons (PNs) in the mammalian olfactory bulb (OB) receive input from the nose and project to diverse cortical and subcortical areas. Morphological and physiological studies have highlighted functional heterogeneity, yet no molecular markers have been described that delineate PN subtypes. Here, we used viral injections into olfactory cortex and fluorescent nucleus sorting to enrich PNs for high-throughput single nucleus and bulk RNA deep sequencing. Transcriptome analysis and RNA in situ hybridization identified distinct mitral and tufted cell populations with characteristic transcription factor network topology, cell adhesion, and excitability-related gene expression. Finally, we describe a new computational approach for integrating bulk and snRNA-seq data and provide evidence that different mitral cell populations preferentially project to different target regions. Together, we have identified potential molecular and gene regulatory mechanisms underlying PN diversity and provide new molecular entry points into studying the diverse functional roles of mitral and tufted cell subtypes.

Identification of Putative Early Atherosclerosis Biomarkers by Unsupervised Deconvolution of Heterogeneous Vascular Proteomes

Coronary artery disease remains a leading cause of death in industrialized nations, and early detection of disease is a critical intervention target to effectively treat patients and manage risk. Proteomic analysis of mixed tissue homogenates may obscure subtle protein changes that occur uniquely in underlying tissue subtypes. The unsupervised 'convex analysis of mixtures' (CAM) tool has previously been shown to effectively segregate cellular subtypes from mixed expression data. In this study, we hypothesized that CAM would identify proteomic information specifically informative to early atherosclerosis lesion involvement that could lead to potential markers of early disease detection. We quantified the proteome of 99 paired abdominal aorta (AA) and left anterior descending coronary artery (LAD) specimens (N = 198 specimens total) acquired during autopsy of young adults free of diagnosed cardiac disease. The CAM tool was then used to segregate protein subsets uniquely associated with different underlying tissue types, yielding markers of normal and fibrous plaque (FP) tissues in LAD and AA (N = 62 lesions markers). CAM-derived FP marker expression was validated against pathologist estimated luminal surface involvement of FP, as well as in an orthogonal cohort of "pure" fibrous plaque, fatty streak, and normal vascular specimens. A targeted mass spectrometry (MS) assay quantified 39 of 62 CAM-FP markers in plasma from women with angiographically verified coronary artery disease (CAD, N = 46) or free from apparent CAD (control, N = 40). Elastic net variable selection with logistic regression reduced this list to 10 proteins capable of classifying CAD status in this cohort with <6% misclassification error, and a mean area under the receiver operating characteristic curve of 0.992 (confidence interval 0.968-0.998) after cross validation. The proteomics-CAM workflow identified lesion-specific molecular biomarker candidates by distilling the most representative molecules from heterogeneous tissue types.

Mathematical modelling of transcriptional heterogeneity identifies novel markers and subpopulations in complex tissues

Tissue heterogeneity is both a major confounding factor and an underexploited information source. While a handful of reports have demonstrated the potential of supervised computational methods to deconvolute tissue heterogeneity, these approaches require a priori information on the marker genes or composition of known subpopulations. To address the critical problem of the absence of validated marker genes for many (including novel) subpopulations, we describe convex analysis of mixtures (CAM), a fully unsupervised in silico method, for identifying subpopulation marker genes directly from the original mixed gene expressions in scatter space that can improve molecular analyses in many biological contexts. Validated with predesigned mixtures, CAM on the gene expression data from peripheral leukocytes, brain tissue, and yeast cell cycle, revealed novel marker genes that were otherwise undetectable using existing methods. Importantly, CAM requires no a priori information on the number, identity, or composition of the subpopulations present in mixed samples, and does not require the presence of pure subpopulations in sample space. This advantage is significant in that CAM can achieve all of its goals using only a small number of heterogeneous samples, and is more powerful to distinguish between phenotypically similar subpopulations.