Mutational signature in colorectal cancer caused by genotoxic pks+ E. coli.

Various species of the intestinal microbiota have been associated with the development of colorectal cancer1,2, but it has not been demonstrated that bacteria have a direct role in the occurrence of oncogenic mutations. Escherichia coli can carry the pathogenicity island pks, which encodes a set of enzymes that synthesize colibactin3. This compound is believed to alkylate DNA on adenine residues4,5 and induces double-strand breaks in cultured cells3. Here we expose human intestinal organoids to genotoxic pks+ E. coli by repeated luminal injection over five months. Whole-genome sequencing of clonal organoids before and after this exposure revealed a distinct mutational signature that was absent from organoids injected with isogenic pks-mutant bacteria. The same mutational signature was detected in a subset of 5,876 human cancer genomes from two independent cohorts, predominantly in colorectal cancer. Our study describes a distinct mutational signature in colorectal cancer and implies that the underlying mutational process results directly from past exposure to bacteria carrying the colibactin-producing pks pathogenicity island.

A recurrent somatic missense mutation in GNAS gene identified in familial thyroid follicular cell carcinomas in German longhaired pointer dogs.

BACKGROUND: We previously reported a familial thyroid follicular cell carcinoma (FCC) in a large number of Dutch German longhaired pointers and identified two deleterious germline mutations in the TPO gene associated with disease predisposition. However, the somatic mutation profile of the FCC in dogs has not been investigated at a genome-wide scale. RESULTS: Herein, we comprehensively investigated the somatic mutations that potentially contribute to the inherited tumor formation and progression using high depth whole-genome sequencing. A GNAS p.A204D missense mutation was identified in 4 out of 7 FCC tumors by whole-genome sequencing and in 20 out of 32 dogs' tumors by targeted sequencing. In contrast to this, in the human TC, mutations in GNAS gene have lower prevalence. Meanwhile, the homologous somatic mutation in humans has not been reported. These findings suggest a difference in the somatic mutation landscape between TC in these dogs and human TC. Moreover, tumors with the GNAS p.A204D mutation had a significantly lower somatic mutation burden in these dogs. Somatic structural variant and copy number alterations were also investigated, but no potential driver event was identified. CONCLUSION: This study provides novel insight in the molecular mechanism of thyroid carcinoma development in dogs. German longhaired pointers carrying GNAS mutations in the tumor may be used as a disease model for the development and testing of novel therapies to kill the tumor with somatic mutations in the GNAS gene.

MutationalPatterns: the one stop shop for the analysis of mutational processes.

BACKGROUND: The collective of somatic mutations in a genome represents a record of mutational processes that have been operative in a cell. These processes can be investigated by extracting relevant mutational patterns from sequencing data. RESULTS: Here, we present the next version of MutationalPatterns, an R/Bioconductor package, which allows in-depth mutational analysis of catalogues of single and double base substitutions as well as small insertions and deletions. Major features of the package include the possibility to perform regional mutation spectra analyses and the possibility to detect strand asymmetry phenomena, such as lesion segregation. On top of this, the package also contains functions to determine how likely it is that a signature can cause damaging mutations (i.e., mutations that affect protein function). This updated package supports stricter signature refitting on known signatures in order to prevent overfitting. Using simulated mutation matrices containing varied signature contributions, we showed that reliable refitting can be achieved even when only 50 mutations are present per signature. Additionally, we incorporated bootstrapped signature refitting to assess the robustness of the signature analyses. Finally, we applied the package on genome mutation data of cell lines in which we deleted specific DNA repair processes and on large cancer datasets, to show how the package can be used to generate novel biological insights. CONCLUSIONS: This novel version of MutationalPatterns allows for more comprehensive analyses and visualization of mutational patterns in order to study the underlying processes. Ultimately, in-depth mutational analyses may contribute to improved biological insights in mechanisms of mutation accumulation as well as aid cancer diagnostics. MutationalPatterns is freely available at http://bioconductor.org/packages/MutationalPatterns .

DNA mismatch repair controls the mutagenicity of Polymerase ζ-dependent translesion synthesis at methylated guanines.

By replicating damaged nucleotides, error-prone DNA translesion synthesis (TLS) enables the completion of replication, albeit at the expense of fidelity. TLS of helix-distorting DNA lesions, that usually have reduced capacity of basepairing, comprises insertion opposite the lesion followed by extension, the latter in particular by polymerase ζ (Pol ζ). However, little is known about involvement of Pol ζ in TLS of non- or poorly-distorting, but miscoding, lesions such as O6-methyldeoxyguanosine (O6-medG). Using purified Pol ζ we describe that the enzyme can misincorporate thymidine opposite O6-medG and efficiently extend from terminal mismatches, suggesting its involvement in the mutagenicity of O6-medG. Surprisingly, O6-medG lesions induced by the methylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) appeared more, rather than less, mutagenic in Pol ζ-deficient mouse embryonic fibroblasts (MEFs) than in wild type MEFs. This suggested that in vivo Pol ζ participates in non-mutagenic TLS of O6-medG. However, we found that the Pol ζ-dependent misinsertions at O6-medG lesions are efficiently corrected by DNA mismatch repair (MMR), which masks the error-proneness of Pol ζ. We also found that the MNNG-induced mutational signature is determined by the adduct spectrum, and modulated by MMR. The signature mimicked single base substitution signature 11 in the catalogue of somatic mutations in cancer, associated with treatment with the methylating drug temozolomide. Our results unravel the individual roles of the major contributors to methylating drug-induced mutagenesis. Moreover, these results warrant caution as to the classification of TLS as mutagenic or error-free based on in vitro data or on the analysis of mutations induced in MMR-proficient cells.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene.

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox.

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole-genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole-genome sequencing (WGS) but still generate hundreds of artifacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolbox (PTATO), to accurately detect single base substitutions, insertions-deletions (indels), and structural variants in PTA-based WGS data. PTATO includes a machine learning approach and filtering based on recurrence to distinguish PTA artifacts from true mutations with high sensitivity (up to 90%), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that hematopoietic stem cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS, have normal somatic single base substitution burdens but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.

Mutation accumulation in mtDNA of cancers resembles mutagenesis in normal stem cells.

Mitochondria are small organelles that play an essential role in the energy production of eukaryotic cells. Defects in their genomes are associated with diseases, such as aging and cancer. Here, we analyzed the mitochondrial genomes of 532 whole-genome sequencing samples from cancers and normal clonally expanded single cells. We show that the mitochondria of normal cells accumulate mutations with age and that most of the mitochondrial mutations found in cancer are the result of healthy mutation accumulation. We also show that the normal HSPCs of patients with leukemia have an increased mitochondrial mutation load. Finally, we show that secondary pediatric cancers and chemotherapy treatments do not impact the mitochondrial mutation load and mtDNA copy numbers of most cells, suggesting that damage to the mitochondrial genome is not a major driver for carcinogenesis. Overall, these findings may contribute to our understanding of mitochondrial genomes and their role in cancer.

Human induced pluripotent stem cells display a similar mutation burden as embryonic pluripotent cells in vivo.

Induced pluripotent stem cells (iPSCs) hold great promise for regenerative medicine, but genetic instability is a major concern. Embryonic pluripotent cells also accumulate mutations during early development, but how this relates to the mutation burden in iPSCs remains unknown. Here, we directly compared the mutation burden of cultured iPSCs with their isogenic embryonic cells during human embryogenesis. We generated developmental lineage trees of human fetuses by phylogenetic inference from somatic mutations in the genomes of multiple stem cells, which were derived from different germ layers. Using this approach, we characterized the mutations acquired pre-gastrulation and found a rate of 1.65 mutations per cell division. When cultured in hypoxic conditions, iPSCs generated from fetal stem cells of the assessed fetuses displayed a similar mutation rate and spectrum. Our results show that iPSCs maintain a genomic integrity during culture at a similar degree as their pluripotent counterparts do in vivo.

Mesenchymal tumor organoid models recapitulate rhabdomyosarcoma subtypes.

Rhabdomyosarcomas (RMS) are mesenchyme-derived tumors and the most common childhood soft tissue sarcomas. Treatment is intense, with a nevertheless poor prognosis for high-risk patients. Discovery of new therapies would benefit from additional preclinical models. Here, we describe the generation of a collection of 19 pediatric RMS tumor organoid (tumoroid) models (success rate of 41%) comprising all major subtypes. For aggressive tumors, tumoroid models can often be established within 4-8 weeks, indicating the feasibility of personalized drug screening. Molecular, genetic, and histological characterization show that the models closely resemble the original tumors, with genetic stability over extended culture periods of up to 6 months. Importantly, drug screening reflects established sensitivities and the models can be modified by CRISPR/Cas9 with TP53 knockout in an embryonal RMS model resulting in replicative stress drug sensitivity. Tumors of mesenchymal origin can therefore be used to generate organoid models, relevant for a variety of preclinical and clinical research questions.

The Dynamics of Somatic Mutagenesis During Life in Humans.

From conception to death, human cells accumulate somatic mutations in their genomes. These mutations can contribute to the development of cancer and non-malignant diseases and have also been associated with aging. Rapid technological developments in sequencing approaches in the last few years and their application to normal tissues have greatly advanced our knowledge about the accumulation of these mutations during healthy aging. Whole genome sequencing studies have revealed that there are significant differences in mutation burden and patterns across tissues, but also that the mutation rates within tissues are surprisingly constant during adult life. In contrast, recent lineage-tracing studies based on whole-genome sequencing have shown that the rate of mutation accumulation is strongly increased early in life before birth. These early mutations, which can be shared by many cells in the body, may have a large impact on development and the origin of somatic diseases. For example, cancer driver mutations can arise early in life, decades before the detection of the malignancy. Here, we review the recent insights in mutation accumulation and mutagenic processes in normal tissues. We compare mutagenesis early and later in life and discuss how mutation rates and patterns evolve during aging. Additionally, we outline the potential impact of these mutations on development, aging and disease.

Mutation signatures of pediatric acute myeloid leukemia and normal blood progenitors associated with differential patient outcomes.

Acquisition of oncogenic mutations with age is believed to be rate limiting for carcinogenesis. However, the incidence of leukemia in children is higher than in young adults. Here we compare somatic mutations across pediatric acute myeloid leukemia (pAML) patient-matched leukemic blasts and hematopoietic stem and progenitor cells (HSPCs), as well as HSPCs from age-matched healthy donors. HSPCs in the leukemic bone marrow have limited genetic relatedness and share few somatic mutations with the cell-of-origin of the malignant blasts, suggesting polyclonal hematopoiesis in pAML patients. Compared to normal HSPCs, a subset of pAML cases harbored more somatic mutations and a distinct composition of mutational process signatures. We hypothesize these cases might have arisen from a more committed progenitor. This subset had better outcomes than pAML cases with mutation burden comparable to age-matched healthy HSPCs. Our study provides insights into the etiology and patient stratification of pAML.

Antiviral treatment causes a unique mutational signature in cancers of transplantation recipients.

Genetic instability is a major concern for successful application of stem cells in regenerative medicine. However, the mutational consequences of the most applied stem cell therapy in humans, hematopoietic stem cell transplantation (HSCT), remain unknown. Here we characterized the mutation burden of hematopoietic stem and progenitor cells (HSPCs) of human HSCT recipients and their donors using whole-genome sequencing. We demonstrate that the majority of transplanted HSPCs did not display altered mutation accumulation. However, in some HSCT recipients, we identified multiple HSPCs with an increased mutation burden after transplantation. This increase could be attributed to a unique mutational signature caused by the antiviral drug ganciclovir. Using a machine learning approach, we detected this signature in cancer genomes of individuals who received HSCT or solid organ transplantation earlier in life. Antiviral treatment with nucleoside analogs can cause enhanced mutagenicity in transplant recipients, which may ultimately contribute to therapy-related carcinogenesis.

Mutation accumulation and developmental lineages in normal and Down syndrome human fetal haematopoiesis.

Children show a higher incidence of leukemia compared to young adolescents, yet their cells have less age-related (oncogenic) somatic mutations. Newborns with Down syndrome have an even higher risk of developing leukemia, which is thought to be driven by mutations that accumulate during fetal development. To characterize mutation accumulation in individual stem and progenitor cells of Down syndrome and karyotypically normal fetuses, we clonally expanded single cells and performed whole-genome sequencing. We found a higher mutation rate in haematopoietic stem and progenitor cells during fetal development compared to the post-infant rate. In fetal trisomy 21 cells the number of somatic mutations is even further increased, which was already apparent during the first cell divisions of embryogenesis before gastrulation. The number and types of mutations in fetal trisomy 21 haematopoietic stem and progenitor cells were similar to those in Down syndrome-associated myeloid preleukemia and could be attributed to mutational processes that were active during normal fetal haematopoiesis. Finally, we found that the contribution of early embryonic cells to human fetal tissues can vary considerably between individuals. The increased mutation rates found in this study, may contribute to the increased risk of leukemia early during life and the higher incidence of leukemia in Down syndrome.

Characterizing Mutational Load and Clonal Composition of Human Blood.

Hematopoietic stem and progenitor cells (HSPCs) gradually accumulate DNA mutations during a lifespan, which can contribute to age-associated diseases such as leukemia. Characterizing mutation accumulation can improve understanding of the etiology of age-associated diseases. Presented here is a method to catalogue somatic mutations in individual HSPCs, which is based on whole-genome sequencing (WGS) of clonal primary cell cultures. Mutations that are present in the original cell are shared by all cells in the clonal culture, whereas mutations acquired in vitro after cell sorting are present in a subset of cells. Therefore, this method allows for accurate detection of somatic mutations present in the genomes of individual HSPCs, which accumulate during life. These catalogues of somatic mutations can provide valuable insights into mutational processes active in the hematopoietic tissue and how these processes contribute to leukemogenesis. In addition, by assessing somatic mutations that are shared between multiple HSPCs of the same individual, clonal lineage relationships and population dynamics of blood populations can be determined. As this approach relies on in vitro expansion of single cells, the method is limited to hematopoietic cells with sufficient replicative potential.