Replication stress impairs chromosome segregation and preimplantation development in human embryos.

Human cleavage-stage embryos frequently acquire chromosomal aneuploidies during mitosis due to unknown mechanisms. Here, we show that S phase at the 1-cell stage shows replication fork stalling, low fork speed, and DNA synthesis extending into G2 phase. DNA damage foci consistent with collapsed replication forks, DSBs, and incomplete replication form in G2 in an ATR- and MRE11-dependent manner, followed by spontaneous chromosome breakage and segmental aneuploidies. Entry into mitosis with incomplete replication results in chromosome breakage, whole and segmental chromosome errors, micronucleation, chromosome fragmentation, and poor embryo quality. Sites of spontaneous chromosome breakage are concordant with sites of DNA synthesis in G2 phase, locating to gene-poor regions with long neural genes, which are transcriptionally silent at this stage of development. Thus, DNA replication stress in mammalian preimplantation embryos predisposes gene-poor regions to fragility, and in particular in the human embryo, to the formation of aneuploidies, impairing developmental potential.

p53 Represses the Mevalonate Pathway to Mediate Tumor Suppression.

There are still gaps in our understanding of the complex processes by which p53 suppresses tumorigenesis. Here we describe a novel role for p53 in suppressing the mevalonate pathway, which is responsible for biosynthesis of cholesterol and nonsterol isoprenoids. p53 blocks activation of SREBP-2, the master transcriptional regulator of this pathway, by transcriptionally inducing the ABCA1 cholesterol transporter gene. A mouse model of liver cancer reveals that downregulation of mevalonate pathway gene expression by p53 occurs in premalignant hepatocytes, when p53 is needed to actively suppress tumorigenesis. Furthermore, pharmacological or RNAi inhibition of the mevalonate pathway restricts the development of murine hepatocellular carcinomas driven by p53 loss. Like p53 loss, ablation of ABCA1 promotes murine liver tumorigenesis and is associated with increased SREBP-2 maturation. Our findings demonstrate that repression of the mevalonate pathway is a crucial component of p53-mediated liver tumor suppression and outline the mechanism by which this occurs.

Histone demethylase LSD1 is required for germinal center formation and BCL6-driven lymphomagenesis.

Germinal center (GC) B cells feature repression of many gene enhancers to establish their characteristic transcriptome. Here we show that conditional deletion of Lsd1 in GCs significantly impaired GC formation, associated with failure to repress immune synapse genes linked to GC exit, which are also direct targets of the transcriptional repressor BCL6. We found that BCL6 directly binds LSD1 and recruits it primarily to intergenic and intronic enhancers. Conditional deletion of Lsd1 suppressed GC hyperplasia caused by constitutive expression of BCL6 and significantly delayed BCL6-driven lymphomagenesis. Administration of catalytic inhibitors of LSD1 had little effect on GC formation or GC-derived lymphoma cells. Using a CRISPR-Cas9 domain screen, we found instead that the LSD1 Tower domain was critical for dependence on LSD1 in GC-derived B cells. These results indicate an essential role for LSD1 in the humoral immune response, where it modulates enhancer function by forming repression complexes with BCL6.

Ordered and deterministic cancer genome evolution after p53 loss.

Although p53 inactivation promotes genomic instability1 and presents a route to malignancy for more than half of all human cancers2,3, the patterns through which heterogenous TP53 (encoding human p53) mutant genomes emerge and influence tumorigenesis remain poorly understood. Here, in a mouse model of pancreatic ductal adenocarcinoma that reports sporadic p53 loss of heterozygosity before cancer onset, we find that malignant properties enabled by p53 inactivation are acquired through a predictable pattern of genome evolution. Single-cell sequencing and in situ genotyping of cells from the point of p53 inactivation through progression to frank cancer reveal that this deterministic behaviour involves four sequential phases-Trp53 (encoding mouse p53) loss of heterozygosity, accumulation of deletions, genome doubling, and the emergence of gains and amplifications-each associated with specific histological stages across the premalignant and malignant spectrum. Despite rampant heterogeneity, the deletion events that follow p53 inactivation target functionally relevant pathways that can shape genomic evolution and remain fixed as homogenous events in diverse malignant populations. Thus, loss of p53-the 'guardian of the genome'-is not merely a gateway to genetic chaos but, rather, can enable deterministic patterns of genome evolution that may point to new strategies for the treatment of TP53-mutant tumours.

Tumor Cell-Intrinsic Factors Underlie Heterogeneity of Immune Cell Infiltration and Response to Immunotherapy.

The biological and functional heterogeneity between tumors-both across and within cancer types-poses a challenge for immunotherapy. To understand the factors underlying tumor immune heterogeneity and immunotherapy sensitivity, we established a library of congenic tumor cell clones from an autochthonous mouse model of pancreatic adenocarcinoma. These clones generated tumors that recapitulated T cell-inflamed and non-T-cell-inflamed tumor microenvironments upon implantation in immunocompetent mice, with distinct patterns of infiltration by immune cell subsets. Co-injecting tumor cell clones revealed the non-T-cell-inflamed phenotype is dominant and that both quantitative and qualitative features of intratumoral CD8+ T cells determine response to therapy. Transcriptomic and epigenetic analyses revealed tumor-cell-intrinsic production of the chemokine CXCL1 as a determinant of the non-T-cell-inflamed microenvironment, and ablation of CXCL1 promoted T cell infiltration and sensitivity to a combination immunotherapy regimen. Thus, tumor cell-intrinsic factors shape the tumor immune microenvironment and influence the outcome of immunotherapy.

Phase and context shape the function of composite oncogenic mutations.

Cancers develop as a result of driver mutations1,2 that lead to clonal outgrowth and the evolution of disease3,4. The discovery and functional characterization of individual driver mutations are central aims of cancer research, and have elucidated myriad phenotypes5 and therapeutic vulnerabilities6. However, the serial genetic evolution of mutant cancer genes7,8 and the allelic context in which they arise is poorly understood in both common and rare cancer genes and tumour types. Here we find that nearly one in four human tumours contains a composite mutation of a cancer-associated gene, defined as two or more nonsynonymous somatic mutations in the same gene and tumour. Composite mutations are enriched in specific genes, have an elevated rate of use of less-common hotspot mutations acquired in a chronology driven in part by oncogenic fitness, and arise in an allelic configuration that reflects context-specific selective pressures. cis-acting composite mutations are hypermorphic in some genes in which dosage effects predominate (such as TERT), whereas they lead to selection of function in other genes (such as TP53). Collectively, composite mutations are driver alterations that arise from context- and allele-specific selective pressures that are dependent in part on gene and mutation function, and which lead to complex-often neomorphic-functions of biological and therapeutic importance.

α-Ketoglutarate links p53 to cell fate during tumour suppression.

The tumour suppressor TP53 is mutated in the majority of human cancers, and in over 70% of pancreatic ductal adenocarcinoma (PDAC)1,2. Wild-type p53 accumulates in response to cellular stress, and regulates gene expression to alter cell fate and prevent tumour development2. Wild-type p53 is also known to modulate cellular metabolic pathways3, although p53-dependent metabolic alterations that constrain cancer progression remain poorly understood. Here we find that p53 remodels cancer-cell metabolism to enforce changes in chromatin and gene expression that favour a premalignant cell fate. Restoring p53 function in cancer cells derived from KRAS-mutant mouse models of PDAC leads to the accumulation of α-ketoglutarate (αKG, also known as 2-oxoglutarate), a metabolite that also serves as an obligate substrate for a subset of chromatin-modifying enzymes. p53 induces transcriptional programs that are characteristic of premalignant differentiation, and this effect can be partially recapitulated by the addition of cell-permeable αKG. Increased levels of the αKG-dependent chromatin modification 5-hydroxymethylcytosine (5hmC) accompany the tumour-cell differentiation that is triggered by p53, whereas decreased 5hmC characterizes the transition from premalignant to de-differentiated malignant lesions that is associated with mutations in Trp53. Enforcing the accumulation of αKG in p53-deficient PDAC cells through the inhibition of oxoglutarate dehydrogenase-an enzyme of the tricarboxylic acid cycle-specifically results in increased 5hmC, tumour-cell differentiation and decreased tumour-cell fitness. Conversely, increasing the intracellular levels of succinate (a competitive inhibitor of αKG-dependent dioxygenases) blunts p53-driven tumour suppression. These data suggest that αKG is an effector of p53-mediated tumour suppression, and that the accumulation of αKG in p53-deficient tumours can drive tumour-cell differentiation and antagonize malignant progression.

Senescence induction dictates response to chemo- and immunotherapy in preclinical models of ovarian cancer.

High-grade serous ovarian carcinoma (HGSOC) is a cancer with dismal prognosis due to the limited effectiveness of existing chemo- and immunotherapies. To elucidate mechanisms mediating sensitivity or resistance to these therapies, we developed a fast and flexible autochthonous mouse model based on somatic introduction of HGSOC-associated genetic alterations into the ovary of immunocompetent mice using tissue electroporation. Tumors arising in these mice recapitulate the metastatic patterns and histological, molecular, and treatment response features of the human disease. By leveraging these models, we show that the ability to undergo senescence underlies the clinically observed increase in sensitivity of homologous recombination (HR)-deficient HGSOC tumors to platinum-based chemotherapy. Further, cGas/STING-mediated activation of a restricted senescence-associated secretory phenotype (SASP) was sufficient to induce immune infiltration and sensitize HR-deficient tumors to immune checkpoint blockade. In sum, our study identifies senescence propensity as a predictor of therapy response and defines a limited SASP profile that appears sufficient to confer added vulnerability to concurrent immunotherapy and, more broadly, provides a blueprint for the implementation of electroporation-based mouse models to reveal mechanisms of oncogenesis and therapy response in HGSOC.

BMP2/SMAD pathway activation in JAK2/p53-mutant megakaryocyte/erythroid progenitors promotes leukemic transformation.

Leukemic transformation (LT) of myeloproliferative neoplasm (MPN) has a dismal prognosis and is largely fatal. Mutational inactivation of TP53 is the most common somatic event in LT; however, the mechanisms by which TP53 mutations promote LT remain unresolved. Using an allelic series of mouse models of Jak2/Trp53 mutant MPN, we identify that only biallelic inactivation of Trp53 results in LT (to a pure erythroleukemia [PEL]). This PEL arises from the megakaryocyte-erythroid progenitor population. Importantly, the bone morphogenetic protein 2/SMAD pathway is aberrantly activated during LT and results in abnormal self-renewal of megakaryocyte-erythroid progenitors. Finally, we identify that Jak2/Trp53 mutant PEL is characterized by recurrent copy number alterations and DNA damage. Using a synthetic lethality strategy, by targeting active DNA repair pathways, we show that this PEL is highly sensitive to combination WEE1 and poly(ADP-ribose) polymerase inhibition. These observations yield new mechanistic insights into the process of p53 mutant LT and offer new, clinically translatable therapeutic approaches.

Mcm2 hypomorph leads to acute leukemia or hematopoietic stem cell failure, dependent on genetic context.

Minichromosome maintenance proteins (Mcm2-7) form a hexameric complex that unwinds DNA ahead of a replicative fork. The deficiency of Mcm proteins leads to replicative stress and consequent genomic instability. Mice with a germline insertion of a Cre cassette into the 3'UTR of the Mcm2 gene (designated Mcm2Cre ) have decreased Mcm2 expression and invariably develop precursor T-cell lymphoblastic leukemia/lymphoma (pre-T LBL), due to 100-1000 kb deletions involving important tumor suppressor genes. To determine whether mice that were protected from pre-T LBL would develop non-T-cell malignancies, we used two approaches. Mice engrafted with Mcm2Cre/Cre Lin- Sca-1+ Kit+ hematopoietic stem/progenitor cells did not develop hematologic malignancy; however, these mice died of hematopoietic stem cell failure by 6 months of age. Placing the Mcm2Cre allele onto an athymic nu/nu background completely prevented pre-T LBL and extended survival of these mice three-fold (median 296.5 vs. 80.5 days). Ultimately, most Mcm2Cre/Cre ;nu/nu mice developed B-cell precursor acute lymphoblastic leukemia (BCP-ALL). We identified recurrent deletions of 100-1000 kb that involved genes known or suspected to be involved in BCP-ALL, including Pax5, Nf1, Ikzf3, and Bcor. Moreover, whole-exome sequencing identified recurrent mutations of genes known to be involved in BCP-ALL progression, such as Jak1/Jak3, Ptpn11, and Kras. These findings demonstrate that an Mcm2Cre/Cre hypomorph can induce hematopoietic dysfunction via hematopoietic stem cell failure as well as a "deletor" phenotype affecting known or suspected tumor suppressor genes.

High resolution copy number inference in cancer using short-molecule nanopore sequencing.

Genome copy number is an important source of genetic variation in health and disease. In cancer, Copy Number Alterations (CNAs) can be inferred from short-read sequencing data, enabling genomics-based precision oncology. Emerging Nanopore sequencing technologies offer the potential for broader clinical utility, for example in smaller hospitals, due to lower instrument cost, higher portability, and ease of use. Nonetheless, Nanopore sequencing devices are limited in the number of retrievable sequencing reads/molecules compared to short-read sequencing platforms, limiting CNA inference accuracy. To address this limitation, we targeted the sequencing of short-length DNA molecules loaded at optimized concentration in an effort to increase sequence read/molecule yield from a single nanopore run. We show that short-molecule nanopore sequencing reproducibly returns high read counts and allows high quality CNA inference. We demonstrate the clinical relevance of this approach by accurately inferring CNAs in acute myeloid leukemia samples. The data shows that, compared to traditional approaches such as chromosome analysis/cytogenetics, short molecule nanopore sequencing returns more sensitive, accurate copy number information in a cost effective and expeditious manner, including for multiplex samples. Our results provide a framework for short-molecule nanopore sequencing with applications in research and medicine, which includes but is not limited to, CNAs.

Deletions linked to TP53 loss drive cancer through p53-independent mechanisms.

Mutations disabling the TP53 tumour suppressor gene represent the most frequent events in human cancer and typically occur through a two-hit mechanism involving a missense mutation in one allele and a 'loss of heterozygosity' deletion encompassing the other. While TP53 missense mutations can also contribute gain-of-function activities that impact tumour progression, it remains unclear whether the deletion event, which frequently includes many genes, impacts tumorigenesis beyond TP53 loss alone. Here we show that somatic heterozygous deletion of mouse chromosome 11B3, a 4-megabase region syntenic to human 17p13.1, produces a greater effect on lymphoma and leukaemia development than Trp53 deletion. Mechanistically, the effect of 11B3 loss on tumorigenesis involves co-deleted genes such as Eif5a and Alox15b (also known as Alox8), the suppression of which cooperates with Trp53 loss to produce more aggressive disease. Our results imply that the selective advantage produced by human chromosome 17p deletion reflects the combined impact of TP53 loss and the reduced dosage of linked tumour suppressor genes.

Complex rearrangements and oncogene amplifications revealed by long-read DNA and RNA sequencing of a breast cancer cell line.

The SK-BR-3 cell line is one of the most important models for HER2+ breast cancers, which affect one in five breast cancer patients. SK-BR-3 is known to be highly rearranged, although much of the variation is in complex and repetitive regions that may be underreported. Addressing this, we sequenced SK-BR-3 using long-read single molecule sequencing from Pacific Biosciences and develop one of the most detailed maps of structural variations (SVs) in a cancer genome available, with nearly 20,000 variants present, most of which were missed by short-read sequencing. Surrounding the important ERBB2 oncogene (also known as HER2), we discover a complex sequence of nested duplications and translocations, suggesting a punctuated progression. Full-length transcriptome sequencing further revealed several novel gene fusions within the nested genomic variants. Combining long-read genome and transcriptome sequencing enables an in-depth analysis of how SVs disrupt the genome and sheds new light on the complex mechanisms involved in cancer genome evolution.

DNAJB1-PRKACA fusion kinase interacts with ß-catenin and the liver regenerative response to drive fibrolamellar hepatocellular carcinoma.

A segmental deletion resulting in DNAJB1-PRKACA gene fusion is now recognized as the signature genetic event of fibrolamellar hepatocellular carcinoma (FL-HCC), a rare but lethal liver cancer that primarily affects adolescents and young adults. Here we implement CRISPR-Cas9 genome editing and transposon-mediated somatic gene transfer to demonstrate that expression of either the endogenous fusion protein or a chimeric cDNA leads to the formation of indolent liver tumors in mice that closely resemble human FL-HCC. Notably, overexpression of the wild-type PRKACA was unable to fully recapitulate the oncogenic activity of DNAJB1-PRKACA, implying that FL-HCC does not simply result from enhanced PRKACA expression. Tumorigenesis was significantly enhanced by genetic activation of ß-catenin, an observation supported by evidence of recurrent Wnt pathway mutations in human FL-HCC, as well as treatment with the hepatotoxin 3,5-diethoxycarbonyl-1,4-dihydrocollidine, which causes tissue injury, inflammation, and fibrosis. Our study validates the DNAJB1-PRKACA fusion kinase as an oncogenic driver and candidate drug target for FL-HCC, and establishes a practical model for preclinical studies to identify strategies to treat this disease.

Optimizing sparse sequencing of single cells for highly multiplex copy number profiling.

Genome-wide analysis at the level of single cells has recently emerged as a powerful tool to dissect genome heterogeneity in cancer, neurobiology, and development. To be truly transformative, single-cell approaches must affordably accommodate large numbers of single cells. This is feasible in the case of copy number variation (CNV), because CNV determination requires only sparse sequence coverage. We have used a combination of bioinformatic and molecular approaches to optimize single-cell DNA amplification and library preparation for highly multiplexed sequencing, yielding a method that can produce genome-wide CNV profiles of up to a hundred individual cells on a single lane of an Illumina HiSeq instrument. We apply the method to human cancer cell lines and biopsied cancer tissue, thereby illustrating its efficiency, reproducibility, and power to reveal underlying genetic heterogeneity and clonal phylogeny. The capacity of the method to facilitate the rapid profiling of hundreds to thousands of single-cell genomes represents a key step in making single-cell profiling an easily accessible tool for studying cell lineage.

Interactive analysis and assessment of single-cell copy-number variations.

We present Ginkgo (http://qb.cshl.edu/ginkgo), a user-friendly, open-source web platform for the analysis of single-cell copy-number variations (CNVs). Ginkgo automatically constructs copy-number profiles of cells from mapped reads and constructs phylogenetic trees of related cells. We validated Ginkgo by reproducing the results of five major studies. After comparing three commonly used single-cell amplification techniques, we concluded that degenerate oligonucleotide-primed PCR is the most consistent for CNV analysis.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility.

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that late replicating regions are established in association with the B compartment and the nuclear lamina from the first cell cycle after fertilization on both maternal and paternal genomes. Late replicating regions contain a relative paucity of active origins and few but long genes and low G/C content. In both bovine and mouse embryos, replication timing patterns are established prior to embryonic genome activation. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to separation of soma and germ cell lineages. Our studies show that the segregation of early and late replicating regions is among the first layers of genome organization established after fertilization.

Somatic mouse models of gastric cancer reveal genotype-specific features of metastatic disease.

Metastatic gastric carcinoma is a highly lethal cancer that responds poorly to conventional and molecularly targeted therapies. Despite its clinical relevance, the mechanisms underlying the behavior and therapeutic response of this disease are poorly understood owing, in part, to a paucity of tractable models. Here we developed methods to somatically introduce different oncogenic lesions directly into the murine gastric epithelium. Genotypic configurations observed in patients produced metastatic gastric cancers that recapitulated the histological, molecular and clinical features of all nonviral molecular subtypes of the human disease. Applying this platform to both wild-type and immunodeficient mice revealed previously unappreciated links between the genotype, organotropism and immune surveillance of metastatic cells, which produced distinct patterns of metastasis that were mirrored in patients. Our results establish a highly portable platform for generating autochthonous cancer models with flexible genotypes and host backgrounds, which can unravel mechanisms of gastric tumorigenesis or test new therapeutic concepts.

Plasticity-induced repression of Irf6 underlies acquired resistance to cancer immunotherapy in pancreatic ductal adenocarcinoma.

Acquired resistance to immunotherapy remains a critical yet incompletely understood biological mechanism. Here, using a mouse model of pancreatic ductal adenocarcinoma (PDAC) to study tumor relapse following immunotherapy-induced responses, we find that resistance is reproducibly associated with an epithelial-to-mesenchymal transition (EMT), with EMT-transcription factors ZEB1 and SNAIL functioning as master genetic and epigenetic regulators of this effect. Acquired resistance in this model is not due to immunosuppression in the tumor immune microenvironment, disruptions in the antigen presentation machinery, or altered expression of immune checkpoints. Rather, resistance is due to a tumor cell-intrinsic defect in T-cell killing. Molecularly, EMT leads to the epigenetic and transcriptional silencing of interferon regulatory factor 6 (Irf6), rendering tumor cells less sensitive to the pro-apoptotic effects of TNF-α. These findings indicate that acquired resistance to immunotherapy may be mediated by programs distinct from those governing primary resistance, including plasticity programs that render tumor cells impervious to T-cell killing.

Metastatic site influences driver gene function in pancreatic cancer.

Driver gene mutations can increase the metastatic potential of the primary tumor1-3, but their role in sustaining tumor growth at metastatic sites is poorly understood. A paradigm of such mutations is inactivation of SMAD4 - a transcriptional effector of TGFß signaling - which is a hallmark of multiple gastrointestinal malignancies4,5. SMAD4 inactivation mediates TGFß's remarkable anti- to pro-tumorigenic switch during cancer progression and can thus influence both tumor initiation and metastasis6-14. To determine whether metastatic tumors remain dependent on SMAD4 inactivation, we developed a mouse model of pancreatic ductal adenocarcinoma (PDAC) that enables Smad4 depletion in the pre-malignant pancreas and subsequent Smad4 reactivation in established metastases. As expected, Smad4 inactivation facilitated the formation of primary tumors that eventually colonized the liver and lungs. By contrast, Smad4 reactivation in metastatic disease had strikingly opposite effects depending on the tumor's organ of residence: suppression of liver metastases and promotion of lung metastases. Integrative multiomic analysis revealed organ-specific differences in the tumor cells' epigenomic state, whereby the liver and lungs harbored chromatin programs respectively dominated by the KLF and RUNX developmental transcription factors, with Klf4 depletion being sufficient to reverse Smad4's tumor-suppressive activity in liver metastases. Our results show how epigenetic states favored by the organ of residence can influence the function of driver genes in metastatic tumors. This organ-specific gene-chromatin interplay invites consideration of anatomical site in the interpretation of tumor genetics, with implications for the therapeutic targeting of metastatic disease.