A Bayesian two-way latent structure model for genomic data integration reveals few pan-genomic cluster subtypes in a breast cancer cohort

NetMIM: network-based multi-omics integration with block missingness for biomarker selection and disease outcome prediction

Compared with analyzing omics data from a single platform, an integrative analysis of multi-omics data provides a more comprehensive understanding of the regulatory relationships among biological features associated with complex diseases. However, most existing frameworks for integrative analysis overlook two crucial aspects of multi-omics data. Firstly, they neglect the known dependencies among biological features that exist in highly credible biological databases. Secondly, most existing integrative frameworks just simply remove the subjects without full omics data to handle block missingness, resulting in decreasing statistical power. To overcome these issues, we propose a network-based integrative Bayesian framework for biomarker selection and disease outcome prediction based on multi-omics data. Our framework utilizes Dirac spike-and-slab variable selection prior to identifying a small subset of biomarkers. The incorporation of gene pathway information improves the interpretability of feature selection. Furthermore, with the strategy in the FBM (stand for "full Bayesian model with missingness") model where missing omics data are augmented via a mechanistic model, our framework handles block missingness in multi-omics data via a data augmentation approach. The real application illustrates that our approach, which incorporates existing gene pathway information and includes subjects without DNA methylation data, results in more interpretable feature selection results and more accurate predictions.

Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

BCL2A1 and CCL18 Are Predictive Biomarkers of Cisplatin Chemotherapy and Immunotherapy in Colon Cancer Patients

Background: Cisplatin enhances the antitumor T cell response, and the combination of PD-L1 blockade produces a synergistic therapeutic effect. However, the clinical correlation between cisplatin and immunotherapy in colon cancer (CC) is unknown.

Methods: Using the “pRRophetic” package, we calculated the IC50 of cisplatin. The correlation between cisplatin IC50, cisplatin resistance–related genes (CCL18 and BCL2A1), and immunotherapy were preliminarily verified in TCGA and further validated in independent cohorts (GSE39582 and GSE17538), cisplatin-resistant CC cell line DLD1, and our own clinical specimens. Classification performance was evaluated using the AUC value of the ROC curve. Scores of immune signatures, autophagy, ferroptosis, and stemness were quantified using the ssGSEA algorithm.

Results: Based on respective medians of three CC cohorts, patients were divided into high- and low-IC50 groups. Compared with the high IC50 group, the low-IC50 group had significantly higher tumor microenvironment (TME) scores and lower tumor purity. Most co-signaling molecules were upregulated in low IC50 group. CC patients with good immunotherapy efficacy (MSI, dMMR, and more TMB) were more attributable to the low-IC50 group. Among seven shared differentially expressed cisplatin resistance–related genes, CCL18 and BCL2A1 had the best predictive efficacy of the above immunotherapy biomarkers. For wet experimental verification, compared with cisplatin-resistant DLD1, similar to PD-L1, CCL18 and BCL2A1 were significantly upregulated in wild-type DLD1. In our own CC tissues, the mRNA expression of CCL18, BCL2A1, and PD-L1 in dMMR were significantly increased. The high group of CCL18 or BCL2A1 had a higher proportion of MSI, dMMR, and more TMB. IC50, CCL18, BCL2A1, and PD-L1 were closely related to scores of immune-related pathways, immune signatures, autophagy, ferroptosis, and stemness. The microRNA shared by BCL2A1 and PD-L1, hsa-miR-137, were significantly associated with CCL18, BCL2A1, and PD-L1, and downregulated in low-IC50 group. The activity of the TOLL-like receptor signaling pathway affected the sensitivity of CC patients to cisplatin and immunotherapy. For subtype analysis, immune C2, immune C6, HM-indel, HM-SNV, C18, and C20 were equally sensitive to cisplatin chemotherapy and immunotherapy.

Conclusions: CC patients sensitive to cisplatin chemotherapy were also sensitive to immunotherapy. CCL18 and BCL2A1 were novel biomarkers for cisplatin and immunotherapy.

Vertical integration methods for gene expression data analysis

Gene expression data have played an essential role in many biomedical studies. When the number of genes is large and sample size is limited, there is a 'lack of information' problem, leading to low-quality findings. To tackle this problem, both horizontal and vertical data integrations have been developed, where vertical integration methods collectively analyze data on gene expressions as well as their regulators (such as mutations, DNA methylation and miRNAs). In this article, we conduct a selective review of vertical data integration methods for gene expression data. The reviewed methods cover both marginal and joint analysis and supervised and unsupervised analysis. The main goal is to provide a sketch of the vertical data integration paradigm without digging into too many technical details. We also briefly discuss potential pitfalls, directions for future developments and application notes.

Bipartite graph-based approach for clustering of cell lines by gene expression-drug response associations

MOTIVATION

In pharmacogenomic studies, the biological context of cell lines influences the predictive ability of drug-response models and the discovery of biomarkers. Thus, similar cell lines are often studied together based on prior knowledge of biological annotations. However, this selection approach is not scalable with the number of annotations, and the relationship between gene-drug association patterns and biological context may not be obvious.

RESULTS

We present a procedure to compare cell lines based on their gene-drug association patterns. Starting with a grouping of cell lines from biological annotation, we model gene-drug association patterns for each group as a bipartite graph between genes and drugs. This is accomplished by applying sparse canonical correlation analysis (SCCA) to extract the gene-drug associations, and using the canonical vectors to construct the edge weights. Then, we introduce a nuclear norm-based dissimilarity measure to compare the bipartite graphs. Accompanying our procedure is a permutation test to evaluate the significance of similarity of cell line groups in terms of gene-drug associations. In the pharmacogenomics datasets CTRP2, GDSC2, and CCLE, hierarchical clustering of carcinoma groups based on this dissimilarity measure uniquely reveals clustering patterns driven by carcinoma subtype rather than primary site. Next, we show that the top associated drugs or genes from SCCA can be used to characterize the clustering patterns of haematopoietic and lymphoid malignancies. Finally, we confirm by simulation that when drug responses are linearly-dependent on expression, our approach is the only one that can effectively infer the true hierarchy compared to existing approaches.

AVAILABILITY

Bipartite graph-based hierarchical clustering is implemented in R and can be obtained from CRAN: https://CRAN.R-project.org/package=hierBipartite. The source code is available at https://github.com/CalvinTChi/hierBipartite.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrating multi-OMICS data through sparse canonical correlation analysis for the prediction of complex traits: a comparison study

Motivation

Recent developments in technology have enabled researchers to collect multiple OMICS datasets for the same individuals. The conventional approach for understanding the relationships between the collected datasets and the complex trait of interest would be through the analysis of each OMIC dataset separately from the rest, or to test for associations between the OMICS datasets. In this work we show that integrating multiple OMICS datasets together, instead of analysing them separately, improves our understanding of their in-between relationships as well as the predictive accuracy for the tested trait. Several approaches have been proposed for the integration of heterogeneous and high-dimensional (p≫n) data, such as OMICS. The sparse variant of canonical correlation analysis (CCA) approach is a promising one that seeks to penalize the canonical variables for producing sparse latent variables while achieving maximal correlation between the datasets. Over the last years, a number of approaches for implementing sparse CCA (sCCA) have been proposed, where they differ on their objective functions, iterative algorithm for obtaining the sparse latent variables and make different assumptions about the original datasets.

Results

Through a comparative study we have explored the performance of the conventional CCA proposed by Parkhomenko et al., penalized matrix decomposition CCA proposed by Witten and Tibshirani and its extension proposed by Suo et al. The aforementioned methods were modified to allow for different penalty functions. Although sCCA is an unsupervised learning approach for understanding of the in-between relationships, we have twisted the problem as a supervised learning one and investigated how the computed latent variables can be used for predicting complex traits. The approaches were extended to allow for multiple (more than two) datasets where the trait was included as one of the input datasets. Both ways have shown improvement over conventional predictive models that include one or multiple datasets.

Availability and implementation

https://github.com/theorod93/sCCA.

Supplementary information

Supplementary data are available at Bioinformatics online.

MONET: Multi-omic module discovery by omic selection

Recent advances in experimental biology allow creation of datasets where several genome-wide data types (called omics) are measured per sample. Integrative analysis of multi-omic datasets in general, and clustering of samples in such datasets specifically, can improve our understanding of biological processes and discover different disease subtypes. In this work we present MONET (Multi Omic clustering by Non-Exhaustive Types), which presents a unique approach to multi-omic clustering. MONET discovers modules of similar samples, such that each module is allowed to have a clustering structure for only a subset of the omics. This approach differs from most existent multi-omic clustering algorithms, which assume a common structure across all omics, and from several recent algorithms that model distinct cluster structures. We tested MONET extensively on simulated data, on an image dataset, and on ten multi-omic cancer datasets from TCGA. Our analysis shows that MONET compares favorably with other multi-omic clustering methods. We demonstrate MONET's biological and clinical relevance by analyzing its results for Ovarian Serous Cystadenocarcinoma. We also show that MONET is robust to missing data, can cluster genes in multi-omic dataset, and reveal modules of cell types in single-cell multi-omic data. Our work shows that MONET is a valuable tool that can provide complementary results to those provided by existent algorithms for multi-omic analysis.

Contrasting DCIS and invasive breast cancer by subtype suggests basal-like DCIS as distinct lesions

Ductal carcinoma in situ (DCIS) is a non-invasive type of breast cancer with highly variable potential of becoming invasive and affecting mortality. Currently, many patients with DCIS are overtreated due to the lack of specific biomarkers that distinguish low risk lesions from those with a higher risk of progression. In this study, we analyzed 57 pure DCIS and 313 invasive breast cancers (IBC) from different patients. Three levels of genomic data were obtained; gene expression, DNA methylation, and DNA copy number. We performed subtype stratified analyses and identified key differences between DCIS and IBC that suggest subtype specific progression. Prominent differences were found in tumors of the basal-like subtype: Basal-like DCIS were less proliferative and showed a higher degree of differentiation than basal-like IBC. Also, core basal tumors (characterized by high correlation to the basal-like centroid) were not identified amongst DCIS as opposed to IBC. At the copy number level, basal-like DCIS exhibited fewer copy number aberrations compared with basal-like IBC. An intriguing finding through analysis of the methylome was hypermethylation of multiple protocadherin genes in basal-like IBC compared with basal-like DCIS and normal tissue, possibly caused by long range epigenetic silencing. This points to silencing of cell adhesion-related genes specifically in IBC of the basal-like subtype. Our work confirms that subtype stratification is essential when studying progression from DCIS to IBC, and we provide evidence that basal-like DCIS show less aggressive characteristics and question the assumption that basal-like DCIS is a direct precursor of basal-like invasive breast cancer.