Cell subpopulation deconvolution reveals breast cancer heterogeneity based on DNA methylation signature

Experimental models of endocrine responsive breast cancer: strengths, limitations, and use

Breast cancers characterized by expression of estrogen receptor-alpha (ER; ESR1) represent approximately 70% of all new cases and comprise the largest molecular subtype of this disease. Despite this high prevalence, the number of adequate experimental models of ER+ breast cancer is relatively limited. Nonetheless, these models have proved very useful in advancing understanding of how cells respond to and resist endocrine therapies, and how the ER acts as a transcription factor to regulate cell fate signaling. We discuss the primary experimental models of ER+ breast cancer including 2D and 3D cultures of established cell lines, cell line- and patient-derived xenografts, and chemically induced rodent models, with a consideration of their respective general strengths and limitations. What can and cannot be learned easily from these models is also discussed, and some observations on how these models may be used more effectively are provided. Overall, despite their limitations, the panel of models currently available has enabled major advances in the field, and these models remain central to the ability to study mechanisms of therapy action and resistance and for hypothesis testing that would otherwise be intractable or unethical in human subjects.

Landscape of Immune Microenvironment Under Immune Cell Infiltration Pattern in Breast Cancer

Background

Increasing evdence supports the suggestion that the immune cell infiltration (ICI) patterns play a pivotal role in tumor progression in breast cancer (BRCA). Nonetheless, there has been no comprehensive analysis of the ICI patterns effects on the clinical outcomes and immunotherapy.

Methods

Multiomic data for BRCA samples were downloaded from TCGA. ESTIMATE algorithm, ssGSEA method, and CIBERSORT analysis were used to uncover the landscape of the tumor immune microenvironment (TIME). BRCA subtypes based on the ICI pattern were identified by consensus clustering and principal-component analysis was performed to obtain the ICI scores to quantify the ICI patterns in individual tumors. Their prognostic value was validated by the Kaplan-Meier survival curves. Gene set enrichment analysis (GSEA) was applied for functional annotation. Immunophenoscore (IPS) was employed to explore the immunotherapeutic role of the ICI scores. Finally, the mutation data was analyzed by using the “maftools” R package.

Results

Three different immune infiltration patterns with a distinct prognosis and biological signature were recognized among 1,198 BRCA samples. The characteristics of TIME under these three patterns were highly consistent with three known immune profiles: immune- excluded, immune-desert, and immune-inflamed phenotypes, respectively. The identification of the ICI patterns within individual tumors based on the ICI score, developed under the ICI-related signature genes, contributed into dissecting biological processes, clinical outcome, immune cells infiltration, immunotherapeutic effect, and genetic variation. High ICI score subtype, characterized with a suppression of immunity, suggested an immune-exhausted phenotype. Abundant effective immune cells were discovered in the low ICI score patients, which corresponded to an immune-activated phenotype and might present an immunotherapeutic advantage. Immunophenoscore was implemented as a surrogate of immunotherapeutic outcome, low-ICI scores samples obtained a significantly higher immunophenoscore. Enrichment of the JAK/STAT and VEGF signal pathways were activated in the ICI low-score subgroup. Finally, the synergistic effect between the ICI score and the tumor mutation burden (TMB) was confirmed.

Conclusion

This work comprehensively elucidated that the ICI patterns served as an indispensable player in complexity and diversity of TIME. Quantitative identification of the ICI patterns in individual tumor will contribute into mapping the landscape of TIME further optimizing precision immunotherapy.

DeCompress: tissue compartment deconvolution of targeted mRNA expression panels using compressed sensing

Targeted mRNA expression panels, measuring up to 800 genes, are used in academic and clinical settings due to low cost and high sensitivity for archived samples. Most samples assayed on targeted panels originate from bulk tissue comprised of many cell types, and cell-type heterogeneity confounds biological signals. Reference-free methods are used when cell-type-specific expression references are unavailable, but limited feature spaces render implementation challenging in targeted panels. Here, we present DeCompress, a semi-reference-free deconvolution method for targeted panels. DeCompress leverages a reference RNA-seq or microarray dataset from similar tissue to expand the feature space of targeted panels using compressed sensing. Ensemble reference-free deconvolution is performed on this artificially expanded dataset to estimate cell-type proportions and gene signatures. In simulated mixtures, four public cell line mixtures, and a targeted panel (1199 samples; 406 genes) from the Carolina Breast Cancer Study, DeCompress recapitulates cell-type proportions with less error than reference-free methods and finds biologically relevant compartments. We integrate compartment estimates into cis-eQTL mapping in breast cancer, identifying a tumor-specific cis-eQTL for CCR3 (C-C Motif Chemokine Receptor 3) at a risk locus. DeCompress improves upon reference-free methods without requiring expression profiles from pure cell populations, with applications in genomic analyses and clinical settings.

Understanding breast cancer heterogeneity through non-genetic heterogeneity

Intricacy in treatment and diagnosis of breast cancer has been an obstacle due to genotype and phenotype heterogeneity. Understanding of non-genetic heterogeneity mechanisms along with considering role of genetic heterogeneity may fill the gaps in landscape painting of heterogeneity. The main factors contribute to non-genetic heterogeneity including: transcriptional pulsing/bursting or discontinuous transcriptions, stochastic partitioning of components at cell division and various signal transduction from tumor ecosystem. Throughout this review, we desired to provide a conceptual framework focused on non-genetic heterogeneity, which has been intended to offer insight into prediction, diagnosis and treatment of breast cancer.

A framework for transcriptome-wide association studies in breast cancer in diverse study populations

BACKGROUND

The relationship between germline genetic variation and breast cancer survival is largely unknown, especially in understudied minority populations who often have poorer survival. Genome-wide association studies (GWAS) have interrogated breast cancer survival but often are underpowered due to subtype heterogeneity and clinical covariates and detect loci in non-coding regions that are difficult to interpret. Transcriptome-wide association studies (TWAS) show increased power in detecting functionally relevant loci by leveraging expression quantitative trait loci (eQTLs) from external reference panels in relevant tissues. However, ancestry- or race-specific reference panels may be needed to draw correct inference in ancestrally diverse cohorts. Such panels for breast cancer are lacking.

RESULTS

We provide a framework for TWAS for breast cancer in diverse populations, using data from the Carolina Breast Cancer Study (CBCS), a population-based cohort that oversampled black women. We perform eQTL analysis for 406 breast cancer-related genes to train race-stratified predictive models of tumor expression from germline genotypes. Using these models, we impute expression in independent data from CBCS and TCGA, accounting for sampling variability in assessing performance. These models are not applicable across race, and their predictive performance varies across tumor subtype. Within CBCS (N = 3,828), at a false discovery-adjusted significance of 0.10 and stratifying for race, we identify associations in black women near AURKA, CAPN13, PIK3CA, and SERPINB5 via TWAS that are underpowered in GWAS.

CONCLUSIONS

We show that carefully implemented and thoroughly validated TWAS is an efficient approach for understanding the genetics underpinning breast cancer outcomes in diverse populations.

A Novel Prostate Cell Type-Specific Gene Signature to Interrogate Prostate Tumor Differentiation Status and Monitor Therapeutic Response (Running Title: Phenotypic Classification of Prostate Tumors)

In this study, we extracted prostate cell-specific gene sets (metagenes) to define the epithelial differentiation status of prostate cancers and, using a deconvolution-based strategy, interrogated thousands of primary and metastatic tumors in public gene profiling datasets. We identified a subgroup of primary prostate tumors with low luminal epithelial enrichment (LumE^low). LumE^low tumors were associated with higher Gleason score and mutational burden, reduced relapse-free and overall survival, and were more likely to progress to castration-resistant prostate cancer (CRPC). Using discriminant function analysis, we generate a predictive 10-gene classifier for clinical implementation. This mini-classifier predicted with high accuracy the luminal status in both primary tumors and CRPCs. Immunohistochemistry for COL4A1, a low-luminal marker, sustained the association of attenuated luminal phenotype with metastatic disease. We found also an association of LumE score with tumor phenotype in genetically engineered mouse models (GEMMs) of prostate cancer. Notably, the metagene approach led to the discovery of drugs that could revert the low luminal status in prostate cell lines and mouse models. This study describes a novel tool to dissect the intrinsic heterogeneity of prostate tumors and provide predictive information on clinical outcome and treatment response in experimental and clinical samples.

Towards precision medicine: advances in 5-hydroxymethylcytosine cancer biomarker discovery in liquid biopsy

Robust and clinically convenient biomarkers for cancer diagnosis, early detection, and prognosis have great potential to improve patient survival and are the key to precision medicine. The advent of next-generation sequencing technologies enables a more sensitive and comprehensive profiling of genetic and epigenetic information in tumor-derived materials. Researchers are now able to monitor the dynamics of tumorigenesis in new dimensions, such as using circulating cell-free DNA (cfDNA) and tumor DNA (ctDNA). Mutation-based assays in liquid biopsy cannot always provide consistent results across studies due partly to intra- and inter-tumoral heterogeneity as well as technical limitations. In contrast, epigenetic analysis of patient-derived cfDNA is a promising alternative, especially for early detection and disease surveillance, because epigenetic modifications are tissue-specific and reflect the dynamic process of cancer progression. Therefore, cfDNA-based epigenetic assays are emerging to be a highly sensitive, minimally invasive tool for cancer diagnosis and prognosis with great potential in future precise care of cancer patients. The major obstacle for applying epigenetic analysis of cfDNA, however, has been the lack of enabling techniques with high sensitivity and technical robustness. In this review, we summarized the advances in epigenome-wide profiling of 5-hydroxymethylcytosine (5hmC) in cfDNA, focusing on the detection approaches and potential role as biomarkers in different cancer types.

A Hybrid Clustering Algorithm for Identifying Cell Types from Single-Cell RNA-Seq Data

Single-cell RNA sequencing (scRNA-seq) has recently brought new insight into cell differentiation processes and functional variation in cell subtypes from homogeneous cell populations. A lack of prior knowledge makes unsupervised machine learning methods, such as clustering, suitable for analyzing scRNA-seq . However, there are several limitations to overcome, including high dimensionality, clustering result instability, and parameter adjustment complexity. In this study, we propose a method by combining structure entropy and k nearest neighbor to identify cell subpopulations in scRNA-seq data. In contrast to existing clustering methods for identifying cell subtypes, minimized structure entropy results in natural communities without specifying the number of clusters. To investigate the performance of our model, we applied it to eight scRNA-seq datasets and compared our method with three existing methods (nonnegative matrix factorization, single-cell interpretation via multikernel learning, and structural entropy minimization principle). The experimental results showed that our approach achieves, on average, better performance in these datasets compared to the benchmark methods.

Integrative Analysis of Omics Big Data

The diversity and huge omics data take biology and biomedicine research and application into a big data era, just like that popular in human society a decade ago. They are opening a new challenge from horizontal data ensemble (e.g., the similar types of data collected from different labs or companies) to vertical data ensemble (e.g., the different types of data collected for a group of person with match information), which requires the integrative analysis in biology and biomedicine and also asks for emergent development of data integration to address the great changes from previous population-guided to newly individual-guided investigations.Data integration is an effective concept to solve the complex problem or understand the complicate system. Several benchmark studies have revealed the heterogeneity and trade-off that existed in the analysis of omics data. Integrative analysis can combine and investigate many datasets in a cost-effective reproducible way. Current integration approaches on biological data have two modes: one is "bottom-up integration" mode with follow-up manual integration, and the other one is "top-down integration" mode with follow-up in silico integration.This paper will firstly summarize the combinatory analysis approaches to give candidate protocol on biological experiment design for effectively integrative study on genomics and then survey the data fusion approaches to give helpful instruction on computational model development for biological significance detection, which have also provided newly data resources and analysis tools to support the precision medicine dependent on the big biomedical data. Finally, the problems and future directions are highlighted for integrative analysis of omics big data.

Whole‑genome scale identification of methylation markers specific for cerebral palsy in monozygotic discordant twins

Cerebral palsy (CP) is a severe type of brain disease affecting movement and posture. Although CP has strong genetic and environmental components, considerable differences in the methylome between monozygotic (MZ) twins discordant for CP implicates epigenetic contributors as well. In order to determine the differences in methylation in patients with CP without interference of the interindividual genomic variation, four pairs of MZ twins discordant for CP were profiled for DNA methylation changes using reduced representation bisulfite sequencing on the genomic-scale. Similar DNA methylation patterns were observed in all samples. However, MZ twins demonstrated higher correlations and closer evolutionary associations compared with the other samples, indicating a stable methylome of MZ twins. A total of 190 differentially methylated genes (DMGs) were identified using Student's t-test, of which 37 genes were hypermethylated in the CP group while the remainders were hypomethylated compared with control group. The identified DMGs were enriched in several cerebral abnormalities, including cerebral cortical atrophy and cerebral atrophy, suggesting that the occurrence of CP may be associated with the methylation alterations. The neighboring genes of DMGs in the protein-protein interaction network were enriched in numerous important functions in essential processes. The results of the present study identified important genes that may epigenetically contribute to the occurrence and development of CP in MZ twins, suggesting that the different prevalence of CP in identical twins may be associated with DNA methylation alterations.

DNA methylation dynamics: identification and functional annotation

DNA methylation is an epigenetic modification of cytosines that undergoes dynamic changes in a temporal, spatial and cell-type-specific manner. Recent advances in technology have permitted the profiling of high-throughput methylomes in large numbers of biological samples. Various computational tools have been developed to identify and analyze DNA methylation dynamics in a variety of critical biological processes. As DNA methylation is becoming increasingly viewed as a dynamic process, the mechanisms governing DNA methylation dynamics and its roles in the transcriptional regulatory network are of great interest. It has been reported that DNA methylation dynamics plays essential roles in multiple biological processes, including development and cancer. As a functional event, the dynamics of DNA methylation have become increasingly relevant to many researchers. Here, we review state-of-the-art advances at three levels (genome-wide identification, regulatory mechanism investigation and the functional annotation) in the field of DNA methylation dynamics, as well as the future perspective of DNA methylation dynamics.