Stem-cell states converge in multistage cutaneous squamous cell carcinoma development.

Stem cells play a critical role in cancer development by contributing to cell heterogeneity, lineage plasticity, and drug resistance. We created gene expression networks from hundreds of mouse tissue samples (both normal and tumor) and integrated these with lineage tracing and single-cell RNA-seq, to identify convergence of cell states in premalignant tumor cells expressing markers of lineage plasticity and drug resistance. Two of these cell states representing multilineage plasticity or proliferation were inversely correlated, suggesting a mutually exclusive relationship. Treatment of carcinomas in vivo with chemotherapy repressed the proliferative state and activated multilineage plasticity whereas inhibition of differentiation repressed plasticity and potentiated responses to cell cycle inhibitors. Manipulation of this cell state transition point may provide a source of potential combinatorial targets for cancer therapy.

Embracing cancer complexity: Hallmarks of systemic disease.

The last 50 years have witnessed extraordinary developments in understanding mechanisms of carcinogenesis, synthesized as the hallmarks of cancer. Despite this logical framework, our understanding of the molecular basis of systemic manifestations and the underlying causes of cancer-related death remains incomplete. Looking forward, elucidating how tumors interact with distant organs and how multifaceted environmental and physiological parameters impinge on tumors and their hosts will be crucial for advances in preventing and more effectively treating human cancers. In this perspective, we discuss complexities of cancer as a systemic disease, including tumor initiation and promotion, tumor micro- and immune macro-environments, aging, metabolism and obesity, cancer cachexia, circadian rhythms, nervous system interactions, tumor-related thrombosis, and the microbiome. Model systems incorporating human genetic variation will be essential to decipher the mechanistic basis of these phenomena and unravel gene-environment interactions, providing a modern synthesis of molecular oncology that is primed to prevent cancers and improve patient quality of life and cancer outcomes.

Mutations, Bottlenecks, and Clonal Sweeps: How Environmental Carcinogens and Genomic Changes Shape Clonal Evolution during Tumor Progression.

The transition from a single, initiated cell to a full-blown malignant tumor involves significant genomic evolution. Exposure to carcinogens-whether directly mutagenic or not-can drive progression toward malignancy, as can stochastic acquisition of cancer-promoting genetic events. Mouse models using both carcinogens and germline genetic manipulations have enabled precise inquiry into the evolutionary dynamics that take place as a tumor progresses from benign to malignant to metastatic stages. Tumor progression is characterized by changes in somatic point mutations and copy-number alterations, even though any single tumor can itself have a high or low burden of genomic alterations. Further, lineage-tracing, single-cell analyses and CRISPR barcoding have revealed the distinct clonal dynamics within benign and malignant tumors. Application of these tools in a range of mouse models can shed unique light on the patterns of clonal evolution that take place in both mouse and human tumors.

Gene networks reveal stem-cell state convergence during preneoplasia and progression to malignancy in multistage skin carcinogenesis.

Adult mammalian stem cells play critical roles in normal tissue homeostasis, as well as in tumor development, by contributing to cell heterogeneity, plasticity, and development of drug resistance. The relationship between different types of normal and cancer stem cells is highly controversial and poorly understood. Here, we carried out gene expression network analysis of normal and tumor samples from genetically heterogeneous mice to create network metagenes for visualization of stem-cell networks, rather than individual stem-cell markers, at the single-cell level during multistage carcinogenesis. We combined this approach with lineage tracing and single-cell RNASeq of stem cells and their progeny, identifying a previously unrecognized hierarchy in which Lgr6+ stem cells from tumors generate progeny that express a range of other stem-cell markers including Sox2, Pitx1, Foxa1, Klf5, and Cd44. Our data identify a convergence of multiple stem-cell and tumor-suppressor pathways in benign tumor cells expressing markers of lineage plasticity and oxidative stress. This same single-cell population expresses network metagenes corresponding to markers of cancer drug resistance in human tumors of the skin, lung and prostate. Treatment of mouse squamous carcinomas in vivo with the chemotherapeutic cis-platin resulted in elevated expression of the genes that mark this cell population. Our data have allowed us to create a simplified model of multistage carcinogenesis that identifies distinct stem-cell states at different stages of tumor progression, thereby identifying networks involved in lineage plasticity, drug resistance, and immune surveillance, providing a rich source of potential targets for cancer therapy.

Circadian PERformance in breast cancer: a germline and somatic genetic study of PER3VNTR polymorphisms and gene co-expression.

Polymorphisms in the PER3 gene have been associated with several human disease phenotypes, including sleep disorders and cancer. In particular, the long allele of a variable number of tandem repeat (VNTR) polymorphism has been previously linked to an increased risk of breast cancer. Here we carried out a combined germline and somatic genetic analysis of the role of the PER3VNRT polymorphism in breast cancer. The combined data from 8284 individuals showed a non-significant trend towards increased breast cancer risk in the 5-repeat allele homozygous carriers (OR = 1.17, 95% CI: 0.97-1.42). We observed allelic imbalance at the PER3 locus in matched blood and tumor DNA samples, showing a significant retention of the long variant (risk) allele in tumor samples, and a preferential loss of the short repetition allele (p = 0.0005). Gene co-expression analysis in healthy and tumoral breast tissue samples uncovered significant associations between PER3 expression levels with those from genes which belong to several cancer-associated pathways. Finally, relapse-free survival (RFS) analysis showed that low expression levels of PER3 were linked to a significant lower RSF in luminal A (p = 3 × 10-12) but not in the rest of breast cancer subtypes.

Targeting KRAS4A splicing through the RBM39/DCAF15 pathway inhibits cancer stem cells.

The commonly mutated human KRAS oncogene encodes two distinct KRAS4A and KRAS4B proteins generated by differential splicing. We demonstrate here that coordinated regulation of both isoforms through control of splicing is essential for development of Kras mutant tumors. The minor KRAS4A isoform is enriched in cancer stem-like cells, where it responds to hypoxia, while the major KRAS4B is induced by ER stress. KRAS4A splicing is controlled by the DCAF15/RBM39 pathway, and deletion of KRAS4A or pharmacological inhibition of RBM39 using Indisulam leads to inhibition of cancer stem cells. Our data identify existing clinical drugs that target KRAS4A splicing, and suggest that levels of the minor KRAS4A isoform in human tumors can be a biomarker of sensitivity to some existing cancer therapeutics.

The critical roles of somatic mutations and environmental tumor-promoting agents in cancer risk.

Cancer is driven by genomic mutations in 'cancer driver' genes, which have essential roles in tumor development. These mutations may be caused by exposure to mutagens in the environment or by endogenous DNA-replication errors in tissue stem cells. Recent observations of abundant mutations, including cancer driver mutations, in histologically normal human tissues suggest that mutations alone are not sufficient for tumor development, thus prompting the question of how single mutant cells give rise to neoplasia. In a concept supported by decades-old data from mouse tumor models, non-mutagenic tumor-promoting agents have been posited to activate the proliferation of dormant mutated cells, thus generating actively growing lesions, with the promotion stage as the rate-limiting step in tumor formation. Non-mutagenic promoting agents, either endogenous or environmental, may therefore have a more important role in human cancer etiology than previously thought.

The mutational signature profile of known and suspected human carcinogens in mice.

Epidemiological studies have identified many environmental agents that appear to significantly increase cancer risk in human populations. By analyzing tumor genomes from mice chronically exposed to 1 of 20 known or suspected human carcinogens, we reveal that most agents do not generate distinct mutational signatures or increase mutation burden, with most mutations, including driver mutations, resulting from tissue-specific endogenous processes. We identify signatures resulting from exposure to cobalt and vinylidene chloride and link distinct human signatures (SBS19 and SBS42) with 1,2,3-trichloropropane, a haloalkane and pollutant of drinking water, and find these and other signatures in human tumor genomes. We define the cross-species genomic landscape of tumors induced by an important compendium of agents with relevance to human health.

Integration of multiple biological contexts reveals principles of synthetic lethality that affect reproducibility.

Synthetic lethal screens have the potential to identify new vulnerabilities incurred by specific cancer mutations but have been hindered by lack of agreement between studies. In the case of KRAS, we identify that published synthetic lethal screen hits significantly overlap at the pathway rather than gene level. Analysis of pathways encoded as protein networks could identify synthetic lethal candidates that are more reproducible than those previously reported. Lack of overlap likely stems from biological rather than technical limitations as most synthetic lethal phenotypes are strongly modulated by changes in cellular conditions or genetic context, the latter determined using a pairwise genetic interaction map that identifies numerous interactions that suppress synthetic lethal effects. Accounting for pathway, cellular and genetic context nominates a DNA repair dependency in KRAS-mutant cells, mediated by a network containing BRCA1. We provide evidence for why most reported synthetic lethals are not reproducible which is addressable using a multi-faceted testing framework.

Mutational signatures in tumours induced by high and low energy radiation in Trp53 deficient mice.

Ionising radiation (IR) is a recognised carcinogen responsible for cancer development in patients previously treated using radiotherapy, and in individuals exposed as a result of accidents at nuclear energy plants. However, the mutational signatures induced by distinct types and doses of radiation are unknown. Here, we analyse the genetic architecture of mammary tumours, lymphomas and sarcomas induced by high (56Fe-ions) or low (gamma) energy radiation in mice carrying Trp53 loss of function alleles. In mammary tumours, high-energy radiation is associated with induction of focal structural variants, leading to genomic instability and Met amplification. Gamma-radiation is linked to large-scale structural variants and a point mutation signature associated with oxidative stress. The genomic architecture of carcinomas, sarcomas and lymphomas arising in the same animals are significantly different. Our study illustrates the complex interactions between radiation quality, germline Trp53 deficiency and tissue/cell of origin in shaping the genomic landscape of IR-induced tumours.

KRAS4A directly regulates hexokinase 1.

The most frequently mutated oncogene in cancer is KRAS, which uses alternative fourth exons to generate two gene products (KRAS4A and KRAS4B) that differ only in their C-terminal membrane-targeting region1. Because oncogenic mutations occur in exons 2 or 3, two constitutively active KRAS proteins-each capable of transforming cells-are encoded when KRAS is activated by mutation2. No functional distinctions among the splice variants have so far been established. Oncogenic KRAS alters the metabolism of tumour cells3 in several ways, including increased glucose uptake and glycolysis even in the presence of abundant oxygen4 (the Warburg effect). Whereas these metabolic effects of oncogenic KRAS have been explained by transcriptional upregulation of glucose transporters and glycolytic enzymes3-5, it is not known whether there is direct regulation of metabolic enzymes. Here we report a direct, GTP-dependent interaction between KRAS4A and hexokinase 1 (HK1) that alters the activity of the kinase, and thereby establish that HK1 is an effector of KRAS4A. This interaction is unique to KRAS4A because the palmitoylation-depalmitoylation cycle of this RAS isoform enables colocalization with HK1 on the outer mitochondrial membrane. The expression of KRAS4A in cancer may drive unique metabolic vulnerabilities that can be exploited therapeutically.

Niche-Specific Factors Dynamically Regulate Sebaceous Gland Stem Cells in the Skin.

Oil-secreting sebaceous glands (SGs) are critical for proper skin function; however, it remains unclear how different factors act together to modulate SG stem cells. Here, we provide functional evidence that each SG lobe is serviced by its own dedicated stem cell population. Upon ablating Notch signaling in different skin subcompartments, we find that this pathway exerts dual counteracting effects on SGs. Suppressing Notch in SG progenitors traps them in a hybrid state where stem and differentiation features become intermingled. In contrast, ablating Notch outside of the SG stem cell compartment indirectly drives SG expansion. Finally, we report that a K14:K5→K14:K79 keratin shift occurs during SG differentiation. Deleting K79 destabilizes K14 in sebocytes, and attenuates SGs and eyelid meibomian glands, leading to corneal ulceration. Altogether, our findings demonstrate that SGs integrate diverse signals from different niches and suggest that mutations incurred within one stem cell compartment can indirectly influence another.

A model for RAS mutation patterns in cancers: finding the sweet spot.

The three RAS genes - HRAS, NRAS and KRAS - are collectively mutated in one-third of human cancers, where they act as prototypic oncogenes. Interestingly, there are rather distinct patterns to RAS mutations; the isoform mutated as well as the position and type of substitution vary between different cancers. As RAS genes are among the earliest, if not the first, genes mutated in a variety of cancers, understanding how these mutation patterns arise could inform on not only how cancer begins but also the factors influencing this event, which has implications for cancer prevention. To this end, we suggest that there is a narrow window or 'sweet spot' by which oncogenic RAS signalling can promote tumour initiation in normal cells. As a consequence, RAS mutation patterns in each normal cell are a product of the specific RAS isoform mutated, as well as the position of the mutation and type of substitution to achieve an ideal level of signalling.

Multicolour lineage tracing reveals clonal dynamics of squamous carcinoma evolution from initiation to metastasis.

Tumour cells are subjected to evolutionary selection pressures during progression from initiation to metastasis. We analysed the clonal evolution of squamous skin carcinomas induced by DMBA/TPA treatment using the K5CreER-Confetti mouse and stage-specific lineage tracing. We show that benign tumours are polyclonal, but only one population contains the Hras driver mutation. Thus, benign papillomas are monoclonal in origin but recruit neighbouring epithelial cells during growth. Papillomas that never progress to malignancy retain several distinct clones, whereas progression to carcinoma is associated with a clonal sweep. Newly generated clones within carcinomas demonstrate intratumoural invasion and clonal intermixing, often giving rise to metastases containing two or more distinct clones derived from the matched primary tumour. These data demonstrate that late-stage tumour progression and dissemination are governed by evolutionary selection pressures that operate at a multicellular level and, therefore, differ from the clonal events that drive initiation and the benign-malignant transition.

Age, estrogen, and immune response in breast adenocarcinoma and adjacent normal tissue.

Chronic inflammation promotes breast tumor growth and invasion by accelerating angiogenesis and tissue remodeling in the tumor microenvironment. There is a complex relationship between inflammation and estrogen, which drives the growth of 70 percent of breast tumors. While low levels of estrogen exposure stimulate macrophages and other inflammatory cell populations, very high levels are immune suppressive. Breast tumor incidence is increased by obesity and age, which interact to influence inflammatory cell populations in normal breast tissue. To characterize the impact of these factors on tumors and the tumor microenvironment, we measured gene expression in 195 breast adenocarcinomas and matched adjacent normal breast tissue samples collected at Akershus University Hospital (AHUS). Age and Body Mass Index (BMI) were independently associated with inflammation in adjacent normal tissue but not tumors. Estrogen Receptor (ER)-negative tumors had elevated macrophage expression compared with matched normal tissue, but ER-positive tumors showed an unexpected decrease in macrophage expression. We found an inverse relationship between the increase in tumor estrogen pathway expression compared with adjacent normal tissue and tumor macrophage score. We validated this finding in 126 breast tumor-normal pairs from the previously published METABRIC cohort. We developed a novel statistic, the Rewiring Coefficient, to quantify the rewiring of gene co-expression networks at the level of individual genes. Differential correlation analysis demonstrated distinct pathways were rewired during tumorigenesis. Our data support an immune suppressive effect of high doses of estrogen signaling in breast tumor microenvironment, suggesting that this effect contributes to the greater presence of prognostic and therapeutically relevant immune cells in ER-negative tumors.

Lgr6 is a stem cell marker in mouse skin squamous cell carcinoma.

The G-protein-coupled receptors LGR4, LGR5 and LGR6 are Wnt signaling mediators, but their functions in squamous cell carcinoma (SCC) are unclear. Using lineage tracing in Lgr5-EGFP-CreERT2/Rosa26-Tomato and Lgr6-EGFP-CreERT2/Rosa26-Tomato reporter mice, we demonstrate that Lgr6, but not Lgr5, acts as an epithelial stem cell marker in SCCs in vivo. We identify, by single-molecule in situ hybridization and cell sorting, rare cells positive for Lgr6 expression in immortalized keratinocytes and show that their frequency increases in advanced SCCs. Lgr6 expression is enriched in cells with stem cell characteristics, and Lgr6 downregulation in vivo causes increased epidermal proliferation with expanded lineage tracing from epidermal stem cells positive for Lgr6 expression. Surprisingly, mice with germline knockout of Lgr6 are predisposed to SCC development, through a mechanism that includes compensatory upregulation of Lgr5. These data provide a model for human patients with germline loss-of-function mutations in Wnt pathway genes, including RSPO1 or LGR4, who show increased susceptibility to squamous tumor development.

Panx3 links body mass index and tumorigenesis in a genetically heterogeneous mouse model of carcinogen-induced cancer.

BACKGROUND: Body mass index (BMI) has been implicated as a primary factor influencing cancer development. However, understanding the relationship between these two complex traits has been confounded by both environmental and genetic heterogeneity. METHODS: In order to gain insight into the genetic factors linking BMI and cancer, we performed chemical carcinogenesis on a genetically heterogeneous cohort of interspecific backcross mice ((Mus Spretus × FVB/N) F1 × FVB/N). Using this cohort, we performed quantitative trait loci (QTL) analysis to identify regions linked to BMI. We then performed an integrated analysis incorporating gene expression, sequence comparison between strains, and gene expression network analysis to identify candidate genes influencing both tumor development and BMI. RESULTS: Analysis of QTL linked to tumorigenesis and BMI identified several loci associated with both phenotypes. Exploring these loci in greater detail revealed a novel relationship between the Pannexin 3 gene (Panx3) and both BMI and tumorigenesis. Panx3 is positively associated with BMI and is strongly tied to a lipid metabolism gene expression network. Pre-treatment Panx3 gene expression levels in normal skin are associated with tumor susceptibility and inhibition of Panx function strongly influences inflammation. CONCLUSIONS: These studies have identified several genetic loci that influence both BMI and carcinogenesis and implicate Panx3 as a candidate gene that links these phenotypes through its effects on inflammation and lipid metabolism.

Gene Expression Architecture of Mouse Dorsal and Tail Skin Reveals Functional Differences in Inflammation and Cancer.

Inherited germline polymorphisms can cause gene expression levels in normal tissues to differ substantially between individuals. We present an analysis of the genetic architecture of normal adult skin from 470 genetically unique mice, demonstrating the effect of germline variants, skin tissue location, and perturbation by exogenous inflammation or tumorigenesis on gene signaling pathways. Gene networks related to specific cell types and signaling pathways, including sonic hedgehog (Shh), Wnt, Lgr family stem cell markers, and keratins, differed at these tissue sites, suggesting mechanisms for the differential susceptibility of dorsal and tail skin to development of skin diseases and tumorigenesis. The Pten tumor suppressor gene network is rewired in premalignant tumors compared to normal tissue, but this response to perturbation is lost during malignant progression. We present a software package for expression quantitative trait loci (eQTL) network analysis and demonstrate how network analysis of whole tissues provides insights into interactions between cell compartments and signaling molecules.