Many forms of cancer have multiple subtypes with different causes and clinical outcomes. Somatic tumor genome sequences provide a rich new source of data for uncovering these subtypes but have proven ...difficult to compare, as two tumors rarely share the same mutations. Here we introduce network-based stratification (NBS), a method to integrate somatic tumor genomes with gene networks. This approach allows for stratification of cancer into informative subtypes by clustering together patients with mutations in similar network regions. We demonstrate NBS in ovarian, uterine and lung cancer cohorts from The Cancer Genome Atlas. For each tissue, NBS identifies subtypes that are predictive of clinical outcomes such as patient survival, response to therapy or tumor histology. We identify network regions characteristic of each subtype and show how mutation-derived subtypes can be used to train an mRNA expression signature, which provides similar information in the absence of DNA sequence.
Stroma is a poorly defined non-parenchymal component of virtually every organ with key roles in organ development, homeostasis, and repair. Studies of the bone marrow stroma have defined individual ...populations in the stem cell niche regulating hematopoietic regeneration and capable of initiating leukemia. Here, we use single-cell RNA sequencing (scRNA-seq) to define a cellular taxonomy of the mouse bone marrow stroma and its perturbation by malignancy. We identified seventeen stromal subsets expressing distinct hematopoietic regulatory genes spanning new fibroblastic and osteoblastic subpopulations including distinct osteoblast differentiation trajectories. Emerging acute myeloid leukemia impaired mesenchymal osteogenic differentiation and reduced regulatory molecules necessary for normal hematopoiesis. These data suggest that tissue stroma responds to malignant cells by disadvantaging normal parenchymal cells. Our taxonomy of the stromal compartment provides a comprehensive bone marrow cell census and experimental support for cancer cell crosstalk with specific stromal elements to impair normal tissue function and thereby enable emergent cancer.
Display omitted
•A single-cell bone marrow stromal cell atlas at steady-state and emergent leukemia•A portal for comparative cell and molecular analyses of bone marrow stromal cells•Distinction among putative niche cells and types of osteolineage differentiation•Leukemia remodeling of stroma to the disadvantage of normal hematopoietic cells
Molecular definition of the cell populations comprising bone marrow stroma is provided with single-cell resolution. Seventeen distinct cell subsets with new mesenchymal, pericyte, fibroblast, and endothelial subpopulations, new inferred osteolineage differentiation trajectories, and distinctions among Lepr, Nestin, and NG2-expressing HSC niche populations are defined. The changes imposed by an emergent AML indicate impaired stromal differentiation among mesenchymal cells and decreased expression of hematopoietic regulatory genes consistent with cancer cells indirectly impairing the normal hematopoietic cells with whom they must compete.
The diversity of mesenchymal cell types in the lung that influence epithelial homeostasis and regeneration is poorly defined. We used genetic lineage tracing, single-cell RNA sequencing, and organoid ...culture approaches to show that Lgr5 and Lgr6, well-known markers of stem cells in epithelial tissues, are markers of mesenchymal cells in the adult lung. Lgr6+ cells comprise a subpopulation of smooth muscle cells surrounding airway epithelia and promote airway differentiation of epithelial progenitors via Wnt-Fgf10 cooperation. Genetic ablation of Lgr6+ cells impairs airway injury repair in vivo. Distinct Lgr5+ cells are located in alveolar compartments and are sufficient to promote alveolar differentiation of epithelial progenitors through Wnt activation. Modulating Wnt activity altered differentiation outcomes specified by mesenchymal cells. This identification of region- and lineage-specific crosstalk between epithelium and their neighboring mesenchymal partners provides new understanding of how different cell types are maintained in the adult lung.
Display omitted
•Lgr5 and Lgr6 mark mesenchymal cells in adult lungs•Single-cell transcriptome analysis defines mesenchymal heterogeneity•Distinct mesenchymal niches drive airway and alveolar differentiation•Wnt activity affects epithelial differentiation specified by mesenchymal cells
Heterogeneous mesenchymal cell populations in the lung play a central role in epithelial maintenance and alveolar differentiation.
Genome-wide association studies (GWAS) have revealed risk alleles for ulcerative colitis (UC). To understand their cell type specificities and pathways of action, we generate an atlas of 366,650 ...cells from the colon mucosa of 18 UC patients and 12 healthy individuals, revealing 51 epithelial, stromal, and immune cell subsets, including BEST4+ enterocytes, microfold-like cells, and IL13RA2+IL11+ inflammatory fibroblasts, which we associate with resistance to anti-TNF treatment. Inflammatory fibroblasts, inflammatory monocytes, microfold-like cells, and T cells that co-express CD8 and IL-17 expand with disease, forming intercellular interaction hubs. Many UC risk genes are cell type specific and co-regulated within relatively few gene modules, suggesting convergence onto limited sets of cell types and pathways. Using this observation, we nominate and infer functions for specific risk genes across GWAS loci. Our work provides a framework for interrogating complex human diseases and mapping risk variants to cell types and pathways.
Display omitted
•51 cell subsets in colon mucosa of 18 ulcerative colitis and 12 healthy individuals•M-like cells, inflammatory monocytes and fibroblasts, and CD8+IL-17+ T cells expand in disease•Oncostatin M circuit in inflammatory monocytes and fibroblasts may affect drug response•Co-expression of genes within cells allows inference of causal genes across risk loci
Single-cell analyses of colon biopsy specimens from patients with ulcerative colitis delineate how expression patterns and shifting cell populations may shape disease and drug resistance, and provide a framework for linking GWAS risk loci with specific cell types and functional pathways.
Epithelial-mesenchymal transition (EMT) programs operate within carcinoma cells, where they generate phenotypes associated with malignant progression. In their various manifestations, EMT programs ...enable epithelial cells to enter into a series of intermediate states arrayed along the E-M phenotypic spectrum. At present, we lack a coherent understanding of how carcinoma cells control their entrance into and continued residence in these various states, and which of these states favour the process of metastasis. Here we characterize a layer of EMT-regulating machinery that governs E-M plasticity (EMP). This machinery consists of two chromatin-modifying complexes, PRC2 and KMT2D-COMPASS, which operate as critical regulators to maintain a stable epithelial state. Interestingly, loss of these two complexes unlocks two distinct EMT trajectories. Dysfunction of PRC2, but not KMT2D-COMPASS, yields a quasi-mesenchymal state that is associated with highly metastatic capabilities and poor survival of patients with breast cancer, suggesting that great caution should be applied when PRC2 inhibitors are evaluated clinically in certain patient cohorts. These observations identify epigenetic factors that regulate EMP, determine specific intermediate EMT states and, as a direct consequence, govern the metastatic ability of carcinoma cells.
Metastasis is the leading cause of death in people with lung cancer, yet the molecular effectors underlying tumor dissemination remain poorly defined. Through the development of an in vivo ...spontaneous lung cancer metastasis model, we show that the developmentally regulated transcriptional repressor Capicua (CIC) suppresses invasion and metastasis. Inactivation of CIC relieves repression of its effector ETV4, driving ETV4-mediated upregulation of MMP24, which is necessary and sufficient for metastasis. Loss of CIC, or an increase in levels of its effectors ETV4 and MMP24, is a biomarker of tumor progression and worse outcomes in people with lung and/or gastric cancer. Our findings reveal CIC as a conserved metastasis suppressor, highlighting new anti-metastatic strategies that could potentially improve patient outcomes.
The implications of intratumor heterogeneity on the four consensus molecular subtypes (CMS) of colorectal cancer (CRC) are not well known. Here, we use single-cell RNA sequencing (scRNASeq) to build ...an algorithm to assign CMS classification to individual cells, which we use to explore the distributions of CMSs in tumor and non-tumor cells. A dataset of colorectal tumors with bulk RNAseq (
= 3232) was used to identify CMS specific-marker gene sets. These gene sets were then applied to a discovery dataset of scRNASeq profiles (
= 10) to develop an algorithm for single-cell CMS (scCMS) assignment, which recapitulated the intrinsic biology of all four CMSs. The single-cell CMS assignment algorithm was used to explore the scRNASeq profiles of two prospective CRC tumors with mixed CMS via bulk sequencing. We find that every CRC tumor contains individual cells of each scCMS, as well as many individual cells that have enrichment for features of more than one scCMS (called mixed cells). scCMS4 and scCMS1 cells dominate stroma and immune cell clusters, respectively, but account for less than 3% epithelial cells. These data imply that CMS1 and CMS4 are driven by the transcriptomic contribution of immune and stromal cells, respectively, not tumor cells.
Head and neck squamous cell carcinoma (HNSCC) is characterized by aggressive behavior with a propensity for metastasis and recurrence. Here we report a comprehensive analysis of the molecular and ...clinical features of HNSCC that govern patient survival. We find that TP53 mutation is frequently accompanied by loss of chromosome 3p and that the combination of these events is associated with a surprising decrease in survival time (1.9 years versus >5 years for TP53 mutation alone). The TP53-3p interaction is specific to chromosome 3p and validates in HNSCC and pan-cancer cohorts. In human papillomavirus (HPV)-positive tumors, in which HPV inactivates TP53, 3p deletion is also common and is associated with poor outcomes. The TP53-3p event is modified by mir-548k expression, which decreases survival further, and is mutually exclusive with mutations affecting RAS signaling. Together, the identified markers underscore the molecular heterogeneity of HNSCC and enable a new multi-tiered classification of this disease.
African-American men have the highest incidence of and mortality from prostate cancer. Whether a biological basis exists for this disparity remains unclear. Exome sequencing (
= 102) and targeted ...validation (
= 90) of localized primary hormone-naïve prostate cancer in African-American men identified several gene mutations not previously observed in this context, including recurrent loss-of-function mutations in
, an ETS transcriptional repressor, in 5% of cases. Analysis of existing prostate cancer cohorts revealed
deletions in 3% of primary prostate cancers and mutations or deletions in
in 3% to 5% of lethal castration-resistant prostate cancers. Knockdown of
confers increased anchorage-independent growth and generates a gene expression signature associated with oncogenic ETS activation and androgen signaling. Together, these results suggest that
is a prostate cancer tumor-suppressor gene. More generally, our findings support the application of systematic cancer genomic characterization in settings of broader ancestral diversity to enhance discovery and, eventually, therapeutic applications.
Systematic genomic sequencing of prostate cancer in African-American men revealed new insights into prostate cancer, including the identification of
as a prostate cancer gene; somatic copy-number alteration differences; and uncommon
and
alterations. This study highlights the importance of inclusion of underrepresented minorities in cancer sequencing studies.
.
Massively parallel sequencing has permitted an unprecedented examination of the cancer exome, leading to predictions that all genes important to cancer will soon be identified by genetic analysis of ...tumours. To examine this potential, here we evaluate the ability of state-of-the-art sequence analysis methods to specifically recover known cancer genes. While some cancer genes are identified by analysis of recurrence, spatial clustering or predicted impact of somatic mutations, many remain undetected due to lack of power to discriminate driver mutations from the background mutational load (13-60% recall of cancer genes impacted by somatic single-nucleotide variants, depending on the method). Cancer genes not detected by mutation recurrence also tend to be missed by all types of exome analysis. Nonetheless, these genes are implicated by other experiments such as functional genetic screens and expression profiling. These challenges are only partially addressed by increasing sample size and will likely hold even as greater numbers of tumours are analysed.