Search

Tasmanian devils have spawned two transmissible cancer lineages, named devil facial tumor 1 (DFT1) and devil facial tumor 2 (DFT2). We investigated the genetic diversity and evolution of these clones by analyzing 78 DFT1 and 41 DFT2 genomes relative to a newly assembled, chromosome-level reference. Time-resolved phylogenetic trees reveal that DFT1 first emerged in 1986 (1982 to 1989) and DFT2 in 2011 (2009 to 2012). Subclone analysis documents transmission of heterogeneous cell populations. DFT2 has faster mutation rates than DFT1 across all variant classes, including substitutions, indels, rearrangements, transposable element insertions, and copy number alterations, and we identify a hypermutated DFT1 lineage with defective DNA mismatch repair. Several loci show plausible evidence of positive selection in DFT1 or DFT2, including loss of chromosome Y and inactivation of MGA , but none are common to both cancers. This study reveals the parallel long-term evolution of two transmissible cancers inhabiting a common niche in Tasmanian devils.

Purpose: Glioblastomas are treated by surgical resection followed by radiotherapy [X-ray therapy (XRT)] and the alkylating chemotherapeutic agent temozolomide. Recently, inactivating mutations in the mismatch repair gene MSH6 were identified in two glioblastomas recurrent post-temozolomide. Because mismatch repair pathway inactivation is a known mediator of alkylator resistance in vitro, these findings suggested that MSH6 inactivation was causally linked to these two recurrences. However, the extent of involvement of MSH6 in glioblastoma is unknown. We sought to determine the overall frequency and clinical relevance of MSH6 alterations in glioblastomas.Experimental Design: The MSH6 gene was sequenced in 54 glioblastomas. MSH6 and O6-methylguanine methyltransferase (MGMT) immunohistochemistry was systematically scored in a panel of 46 clinically well-characterized glioblastomas, and the corresponding patient response to treatment evaluated.Results:MSH6 mutation was not observed in any pretreatment glioblastoma (0 of 40), whereas 3 of 14 recurrent cases had somatic mutations (P = 0.015). MSH6 protein expression was detected in all pretreatment (17 of 17) cases examined but, notably, expression was lost in 7 of 17 (41%) recurrences from matched post–XRT + temozolomide cases (P = 0.016). Loss of MSH6 was not associated with O6-methylguanine methyltransferase status. Measurements of in vivo tumor growth using three-dimensional reconstructed magnetic resonance imaging showed that MSH6-negative glioblastomas had a markedly increased rate of growth while under temozolomide treatment (3.17 versus 0.04 cc/mo for MSH6-positive tumors; P = 0.020).Conclusions: Loss of MSH6 occurs in a subset of post–XRT + temozolomide glioblastoma recurrences and is associated with tumor progression during temozolomide treatment, mirroring the alkylator resistance conferred by MSH6 inactivation in vitro. MSH6 deficiency may therefore contribute to the emergence of recurrent glioblastomas during temozolomide treatment.

Substitution signature exposures estimated in 640 breast cancer WGS and the ovarian sample (AOCS-166).

Mismatch repair (MMR)–deficient cancers have been discovered to be highly responsive to immune therapies such as PD-1 checkpoint blockade, making their definition in patients, where they may be relatively rare, paramount for treatment decisions. In this study, we utilized patterns of mutagenesis known as mutational signatures, which are imprints of the mutagenic processes associated with MMR deficiency, to identify MMR-deficient breast tumors from a whole-genome sequencing dataset comprising a cohort of 640 patients. We identified 11 of 640 tumors as MMR deficient, but only 2 of 11 exhibited germline mutations in MMR genes or Lynch Syndrome. Two additional tumors had a substantially reduced proportion of mutations attributed to MMR deficiency, where the predominant mutational signatures were related to APOBEC enzymatic activity. Overall, 6 of 11 of the MMR-deficient cases in this cohort were confirmed genetically or epigenetically as having abrogation of MMR genes. However, IHC analysis of MMR-related proteins revealed all but one of 10 samples available for testing as MMR deficient. Thus, the mutational signatures more faithfully reported MMR deficiency than sequencing of MMR genes, because they represent a direct pathophysiologic readout of repair pathway abnormalities. As whole-genome sequencing continues to become more affordable, it could be used to expose individually abnormal tumors in tissue types where MMR deficiency has been rarely detected, but also rarely sought. Cancer Res; 77(18); 4755–62. ©2017 AACR.

Mutations in MMR genes

Detection of MMR deficient breast cancers in whole exome sequencing data (WES)

Supplementary Table S1 from Estrogen Receptor Status Could Modulate the Genomic Pattern in Familial and Sporadic Breast Cancer

Supplementary Figures S1-S2 from Estrogen Receptor Status Could Modulate the Genomic Pattern in Familial and Sporadic Breast Cancer

Details of 11 MMR deficient breast cancers.

Supplementary Table S1 from Loss of the Mismatch Repair Protein MSH6 in Human Glioblastomas Is Associated with Tumor Progression during Temozolomide Treatment

This file contains the following Supplementary figures and Supplementary Material and Methods: Supplementary Figure 1: Contribution of base substitution signatures and indels in 640 breast cancer WGS. Supplementary Figure 2: Whole genome plots for 11 MMR deficient breast cancers Supplementary Figure 3: Genome plots for the two samples with germline PMS2 mutations. Supplementary Figure 4: Scatter plot demonstrating the relationship between total number of substitutions associated with the three MMR signatures 6, 20, and 26 and number of indels. Supplementary Materials and Methods 1. Mutational signatures background 2. Detection of MMR deficient cancers in other data sets 2.1 Detection of a genomic profile characteristic of PMS2 deficient tumours in Ovarian cancer. 2.2 Detection of MMR deficient breast cancers in WES

Supplementary Material from In vitro differential sensitivity of melanomas to phenothiazines is based on the presence of codon 600 BRAF mutation

Supplementary Table 6 from Novel Candidate Cancer Genes Identified by a Large-Scale Cross-Species Comparative Oncogenomics Approach

Supplementary Information from A Hypermutation Phenotype and Somatic MSH6 Mutations in Recurrent Human Malignant Gliomas after Alkylator Chemotherapy

Supplementary Glioma Genes from A Hypermutation Phenotype and Somatic MSH6 Mutations in Recurrent Human Malignant Gliomas after Alkylator Chemotherapy

Comparative genomic hybridization (CGH) can reveal important disease genes but the large regions identified could sometimes contain hundreds of genes. Here we combine high-resolution CGH analysis of 598 human cancer cell lines with insertion sites isolated from 1,005 mouse tumors induced with the murine leukemia virus (MuLV). This cross-species oncogenomic analysis revealed candidate tumor suppressor genes and oncogenes mutated in both human and mouse tumors, making them strong candidates for novel cancer genes. A significant number of these genes contained binding sites for the stem cell transcription factors Oct4 and Nanog. Notably, mice carrying tumors with insertions in or near stem cell module genes, which are thought to participate in cell self-renewal, died significantly faster than mice without these insertions. A comparison of the profile we identified to that induced with the Sleeping Beauty (SB) transposon system revealed significant differences in the profile of recurrently mutated genes. Collectively, this work provides a rich catalogue of new candidate cancer genes for functional analysis. Cancer Res; 70(3); 883–95

Supplementary Methods from Novel Candidate Cancer Genes Identified by a Large-Scale Cross-Species Comparative Oncogenomics Approach

Supplementary Glioma Gene Mutations from A Hypermutation Phenotype and Somatic MSH6 Mutations in Recurrent Human Malignant Gliomas after Alkylator Chemotherapy

Supplementary Tables 1-5 from Novel Candidate Cancer Genes Identified by a Large-Scale Cross-Species Comparative Oncogenomics Approach

Malignant gliomas have a very poor prognosis. The current standard of care for these cancers consists of extended adjuvant treatment with the alkylating agent temozolomide after surgical resection and radiotherapy. Although a statistically significant increase in survival has been reported with this regimen, nearly all gliomas recur and become insensitive to further treatment with this class of agents. We sequenced 500 kb of genomic DNA corresponding to the kinase domains of 518 protein kinases in each of nine gliomas. Large numbers of somatic mutations were observed in two gliomas recurrent after alkylating agent treatment. The pattern of mutations in these cases showed strong similarity to that induced by alkylating agents in experimental systems. Further investigation revealed inactivating somatic mutations of the mismatch repair gene MSH6 in each case. We propose that inactivating somatic mutations of MSH6 confer resistance to alkylating agents in gliomas in vivo and concurrently unleash accelerated mutagenesis in resistant clones as a consequence of continued exposure to alkylating agents in the presence of defective mismatch repair. The evidence therefore suggests that when MSH6 is inactivated in gliomas, alkylating agents convert from induction of tumor cell death to promotion of neoplastic progression. These observations highlight the potential of large scale sequencing for revealing and elucidating mutagenic processes operative in individual human cancers. (Cancer Res 2006; 66(8): 3987-91)

Supplementary Table 7 from Novel Candidate Cancer Genes Identified by a Large-Scale Cross-Species Comparative Oncogenomics Approach

The endometrium, the mucosal lining of the uterus, undergoes dynamic changes throughout the menstrual cycle in response to ovarian hormones. We have generated dense single-cell and spatial reference maps of the human uterus and three-dimensional endometrial organoid cultures. We dissect the signaling pathways that determine cell fate of the epithelial lineages in the lumenal and glandular microenvironments. Our benchmark of the endometrial organoids reveals the pathways and cell states regulating differentiation of the secretory and ciliated lineages both in vivo and in vitro. In vitro downregulation of WNT or NOTCH pathways increases the differentiation efficiency along the secretory and ciliated lineages, respectively. We utilize our cellular maps to deconvolute bulk data from endometrial cancers and endometriotic lesions, illuminating the cell types dominating in each of these disorders. These mechanistic insights provide a platform for future development of treatments for common conditions including endometriosis and endometrial carcinoma. Single-cell and spatial transcriptomic profiling of the human endometrium highlights pathways governing the proliferative and secretory phases of the menstrual cycle. Analyses of endometrial organoids show that WNT and NOTCH signaling modulate differentiation into the secretory and ciliated epithelial lineages, respectively. Correction in: Nature Genetics, 55, page 165 (2023)DOI: 10.1038/s41588-022-01287-6

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues. published

Esophageal squamous cell carcinoma (ESCC) shows remarkable variation in incidence that is not fully explained by known lifestyle and environmental risk factors. It has been speculated that an unknown exogenous exposure(s) could be responsible. Here we combine the fields of mutational signature analysis with cancer epidemiology to study 552 ESCC genomes from eight countries with varying incidence rates. Mutational profiles were similar across all countries studied. Associations between specific mutational signatures and ESCC risk factors were identified for tobacco, alcohol, opium and germline variants, with modest impacts on mutation burden. We find no evidence of a mutational signature indicative of an exogenous exposure capable of explaining differences in ESCC incidence. Apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like (APOBEC)-associated mutational signatures single-base substitution (SBS)2 and SBS13 were present in 88% and 91% of cases, respectively, and accounted for 25% of the mutation burden on average, indicating that APOBEC activation is a crucial step in ESCC tumor development.

The progression of chronic liver disease to hepatocellular carcinoma is caused by the acquisition of somatic mutations that affect 20–30 cancer genes1–8. Burdens of somatic mutations are higher and clonal expansions larger in chronic liver disease9–13 than in normal liver13–16, which enables positive selection to shape the genomic landscape9–13. Here we analysed somatic mutations from 1,590 genomes across 34 liver samples, including healthy controls, alcohol-related liver disease and non-alcoholic fatty liver disease. Seven of the 29 patients with liver disease had mutations in FOXO1, the major transcription factor in insulin signalling. These mutations affected a single hotspot within the gene, impairing the insulin-mediated nuclear export of FOXO1. Notably, six of the seven patients with FOXO1S22W hotspot mutations showed convergent evolution, with variants acquired independently by up to nine distinct hepatocyte clones per patient. CIDEB, which regulates lipid droplet metabolism in hepatocytes17–19, and GPAM, which produces storage triacylglycerol from free fatty acids20,21, also had a significant excess of mutations. We again observed frequent convergent evolution: up to fourteen independent clones per patient with CIDEB mutations and up to seven clones per patient with GPAM mutations. Mutations in metabolism genes were distributed across multiple anatomical segments of the liver, increased clone size and were seen in both alcohol-related liver disease and non-alcoholic fatty liver disease, but rarely in hepatocellular carcinoma. Master regulators of metabolic pathways are a frequent target of convergent somatic mutation in alcohol-related and non-alcoholic fatty liver disease. Whole-genome sequencing analysis of somatic mutations in liver samples from patients with chronic liver disease identifies driver mutations in metabolism-related genes such as FOXO1, and shows that these variants frequently exhibit convergent evolution.

Funder: Wellcome Clinical PhD fellowship, Funder: Wellcome PhD Studentship, Funder: Jean Shank/Pathological Society Intermediate Fellowship, Mutation accumulation in somatic cells contributes to cancer development and is proposed as a cause of aging. DNA polymerases Pol ε and Pol δ replicate DNA during cell division. However, in some cancers, defective proofreading due to acquired POLE/POLD1 exonuclease domain mutations causes markedly elevated somatic mutation burdens with distinctive mutational signatures. Germline POLE/POLD1 mutations cause familial cancer predisposition. Here, we sequenced normal tissue and tumor DNA from individuals with germline POLE/POLD1 mutations. Increased mutation burdens with characteristic mutational signatures were found in normal adult somatic cell types, during early embryogenesis and in sperm. Thus human physiology can tolerate ubiquitously elevated mutation burdens. Except for increased cancer risk, individuals with germline POLE/POLD1 mutations do not exhibit overt features of premature aging. These results do not support a model in which all features of aging are attributable to widespread cell malfunction directly resulting from somatic mutation burdens accrued during life.

Starting from the zygote, all cells in the human body continuously acquire mutations. Mutations shared between different cells imply a common progenitor and are thus naturally occurring markers for lineage tracing1,2. Here we reconstruct extensive phylogenies of normal tissues from three adult individuals using whole-genome sequencing of 511 laser capture microdissections. Reconstructed embryonic progenitors in the same generation of a phylogeny often contribute to different extents to the adult body. The degree of this asymmetry varies between individuals, with ratios between the two reconstructed daughter cells of the zygote ranging from 60:40 to 93:7. Asymmetries pervade subsequent generations and can differ between tissues in the same individual. The phylogenies resolve the spatial embryonic patterning of tissues, revealing contiguous patches of, on average, 301 crypts in the adult colonic epithelium derived from a most recent embryonic cell and also a spatial effect in brain development. Using data from ten additional men, we investigated the developmental split between soma and germline, with results suggesting an extraembryonic contribution to primordial germ cells. This research demonstrates that, despite reaching the same ultimate tissue patterns, early bottlenecks and lineage commitments lead to substantial variation in embryonic patterns both within and between individuals. Somatic mutations obtained from laser microdissected biopsies of human tissues are used to reconstruct the developmental phylogenies of these tissues back to the zygote.

Somatic mutations accumulate in healthy tissues as we age, giving rise to cancer and potentially contributing to ageing. To study somatic mutations in non-neoplastic tissues, we developed a series of protocols to sequence the genomes of small populations of cells isolated from histological sections. Here, we describe a complete workflow that combines laser-capture microdissection (LCM) with low-input genome sequencing, while circumventing the use of whole-genome amplification (WGA). The protocol is subdivided broadly into four steps: tissue processing, LCM, low-input library generation and mutation calling and filtering. The tissue processing and LCM steps are provided as general guidelines that might require tailoring based on the specific requirements of the study at hand. Our protocol for low-input library generation uses enzymatic rather than acoustic fragmentation to generate WGA-free whole-genome libraries. Finally, the mutation calling and filtering strategy has been adapted from previously published protocols to account for artifacts introduced via library creation. To date, we have used this workflow to perform targeted and whole-genome sequencing of small populations of cells (typically 100–1,000 cells) in thousands of microbiopsies from a wide range of human tissues. The low-input DNA protocol is designed to be compatible with liquid handling platforms and make use of equipment and expertise standard to any core sequencing facility. However, obtaining low-input DNA material via LCM requires specialized equipment and expertise. The entire protocol from tissue reception through whole-genome library generation can be accomplished in as little as 1 week, although 2–3 weeks would be a more typical turnaround time.