Time-resolved, integrated analysis of clonally evolving genomes.

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

Differentiation-related epigenomic changes define clinically distinct keratinocyte cancer subclasses.

Keratinocyte cancers (KC) are the most prevalent malignancies in fair-skinned populations, posing a significant medical and economic burden to health systems. KC originate in the epidermis and mainly comprise basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC). Here, we combined single-cell multi-omics, transcriptomics, and methylomics to investigate the epigenomic dynamics during epidermal differentiation. We identified ~3,800 differentially accessible regions between undifferentiated and differentiated keratinocytes, corresponding to regulatory regions associated with key transcription factors. DNA methylation at these regions defined AK/cSCC subtypes with epidermal stem cell- or keratinocyte-like features. Using cell-type deconvolution tools and integration of bulk and single-cell methylomes, we demonstrate that these subclasses are consistent with distinct cells-of-origin. Further characterization of the phenotypic traits of the subclasses and the study of additional unstratified KC entities uncovered distinct clinical features for the subclasses, linking invasive and metastatic KC cases with undifferentiated cells-of-origin. Our study provides a thorough characterization of the epigenomic dynamics underlying human keratinocyte differentiation and uncovers novel links between KC cells-of-origin and their prognosis.

DAZAP2 acts as specifier of the p53 response to DNA damage.

The DNA damage-responsive tumor suppressors p53 and HIPK2 are well established regulators of cell fate decision-making and regulate the cellular sensitivity to DNA-damaging drugs. Here, we identify Deleted in Azoospermia-associated protein 2 (DAZAP2), a small adaptor protein, as a novel regulator of HIPK2 and specifier of the DNA damage-induced p53 response. Knock-down or genetic deletion of DAZAP2 strongly potentiates cancer cell chemosensitivity both in cells and in vivo using a mouse tumour xenograft model. In unstressed cells, DAZAP2 stimulates HIPK2 polyubiquitination and degradation through interplay with the ubiquitin ligase SIAH1. Upon DNA damage, HIPK2 site-specifically phosphorylates DAZAP2, which terminates its HIPK2-degrading function and triggers its re-localization to the cell nucleus. Interestingly, nuclear DAZAP2 interacts with p53 and specifies target gene expression through modulating a defined subset of p53 target genes. Furthermore, our results suggest that DAZAP2 co-occupies p53 response elements to specify target gene expression. Collectively, our findings propose DAZAP2 as novel regulator of the DNA damage-induced p53 response that controls cancer cell chemosensitivity.

The microbiota programs DNA methylation to control intestinal homeostasis and inflammation.

Although much research has been done on the diversity of the gut microbiome, little is known about how it influences intestinal homeostasis under normal and pathogenic conditions. Epigenetic mechanisms have recently been suggested to operate at the interface between the microbiota and the intestinal epithelium. We performed whole-genome bisulfite sequencing on conventionally raised and germ-free mice, and discovered that exposure to commensal microbiota induced localized DNA methylation changes at regulatory elements, which are TET2/3-dependent. This culminated in the activation of a set of 'early sentinel' response genes to maintain intestinal homeostasis. Furthermore, we demonstrated that exposure to the microbiota in dextran sodium sulfate-induced acute inflammation results in profound DNA methylation and chromatin accessibility changes at regulatory elements, leading to alterations in gene expression programs enriched in colitis- and colon-cancer-associated functions. Finally, by employing genetic interventions, we show that microbiota-induced epigenetic programming is necessary for proper intestinal homeostasis in vivo.

DNA (de)methylation in embryonic stem cells controls CTCF-dependent chromatin boundaries.

Coordinated changes of DNA (de)methylation, nucleosome positioning, and chromatin binding of the architectural protein CTCF play an important role for establishing cell-type-specific chromatin states during differentiation. To elucidate molecular mechanisms that link these processes, we studied the perturbed DNA modification landscape in mouse embryonic stem cells (ESCs) carrying a double knockout (DKO) of the Tet1 and Tet2 dioxygenases. These enzymes are responsible for the conversion of 5-methylcytosine (5mC) into its hydroxymethylated (5hmC), formylated (5fC), or carboxylated (5caC) forms. We determined changes in nucleosome positioning, CTCF binding, DNA methylation, and gene expression in DKO ESCs and developed biophysical models to predict differential CTCF binding. Methylation-sensitive nucleosome repositioning accounted for a significant portion of CTCF binding loss in DKO ESCs, whereas unmethylated and nucleosome-depleted CpG islands were enriched for CTCF sites that remained occupied. A number of CTCF sites also displayed direct correlations with the CpG modification state: CTCF was preferentially lost from sites that were marked with 5hmC in wild-type (WT) cells but not from 5fC-enriched sites. In addition, we found that some CTCF sites can act as bifurcation points defining the differential methylation landscape. CTCF loss from such sites, for example, at promoters, boundaries of chromatin loops, and topologically associated domains (TADs), was correlated with DNA methylation/demethylation spreading and can be linked to down-regulation of neighboring genes. Our results reveal a hierarchical interplay between cytosine modifications, nucleosome positions, and DNA sequence that determines differential CTCF binding and regulates gene expression.

Methylation profiling identifies two subclasses of squamous cell carcinoma related to distinct cells of origin.

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer and usually progresses from a UV-induced precancerous lesion termed actinic keratosis (AK). Despite various efforts to characterize these lesions molecularly, the etiology of AK and its progression to cSCC remain partially understood. Here, we use Infinium MethylationEPIC BeadChips to interrogate the DNA methylation status in healthy, AK and cSCC epidermis samples. Importantly, we show that AK methylation patterns already display classical features of cancer methylomes and are highly similar to cSCC profiles. Further analysis identifies typical features of stem cell methylomes, such as reduced DNA methylation age, non-CpG methylation, and stem cell-related keratin and enhancer methylation patterns. Interestingly, this signature is detected only in half of the samples, while the other half shows patterns more closely related to healthy epidermis. These findings suggest the existence of two subclasses of AK and cSCC emerging from distinct keratinocyte differentiation stages.

mIDH-associated DNA hypermethylation in acute myeloid leukemia reflects differentiation blockage rather than inhibition of TET-mediated demethylation.

Isocitrate dehydrogenases 1 and 2 (IDH1/2) are recurrently mutated in acute myeloid leukemia (AML), but their mechanistic role in leukemogenesis is poorly understood. The inhibition of TET enzymes by D-2-hydroxyglutarate (D-2-HG), which is produced by mutant IDH1/2 (mIDH1/2), has been suggested to promote epigenetic deregulation during tumorigenesis. In addition, mIDH also induces a differentiation block in various cell culture and mouse models. Here we analyze the genomic methylation patterns of AML patients with mIDH using Infinium 450K data from a large AML cohort and found that mIDH is associated with pronounced DNA hypermethylation at tens of thousands of CpGs. Interestingly, however, myeloid leukemia cells overexpressing mIDH, cells that were cultured in the presence of D-2-HG or TET2 mutant AML patients did not show similar methylation changes. In further analyses, we also characterized the methylation landscapes of myeloid progenitor cells and analyzed their relationship to mIDH-associated hypermethylation. Our findings identify the differentiation state of myeloid cells, rather than inhibition of TET-mediated DNA demethylation, as a major factor of mIDH-associated hypermethylation in AML. Furthermore, our results are also important for understanding the mode of action of currently developed mIDH inhibitors.

Tet1 and Tet2 Protect DNA Methylation Canyons against Hypermethylation.

DNA methylation is a dynamic epigenetic modification with an important role in cell fate specification and reprogramming. The Ten eleven translocation (Tet) family of enzymes converts 5-methylcytosine to 5-hydroxymethylcytosine, which promotes passive DNA demethylation and functions as an intermediate in an active DNA demethylation process. Tet1/Tet2 double-knockout mice are characterized by developmental defects and epigenetic instability, suggesting a requirement for Tet-mediated DNA demethylation for the proper regulation of gene expression during differentiation. Here, we used whole-genome bisulfite and transcriptome sequencing to characterize the underlying mechanisms. Our results uncover the hypermethylation of DNA methylation canyons as the genomic key feature of Tet1/Tet2 double-knockout mouse embryonic fibroblasts. Canyon hypermethylation coincided with disturbed regulation of associated genes, suggesting a mechanistic explanation for the observed Tet-dependent differentiation defects. Based on these results, we propose an important regulatory role of Tet-dependent DNA demethylation for the maintenance of DNA methylation canyons, which prevents invasive DNA methylation and allows functional regulation of canyon-associated genes.

Chronic inflammation induces a novel epigenetic program that is conserved in intestinal adenomas and in colorectal cancer.

Chronic inflammation represents a major risk factor for tumor formation, but the underlying mechanisms have remained largely unknown. Epigenetic mechanisms can record the effects of environmental challenges on the genome level and could therefore play an important role in the pathogenesis of inflammation-associated tumors. Using single-base methylation maps and transcriptome analyses of a colitis-induced mouse colon cancer model, we identified a novel epigenetic program that is characterized by hypermethylation of DNA methylation valleys that are characterized by low CpG density and active chromatin marks. This program is conserved and functional in mouse intestinal adenomas and results in silencing of active intestinal genes that are involved in gastrointestinal homeostasis and injury response. Further analyses reveal that the program represents a prominent feature of human colorectal cancer and can be used to correctly classify colorectal cancer samples with high accuracy. Together, our results show that inflammatory signals establish a novel epigenetic program that silences a specific set of genes that contribute to inflammation-induced cellular transformation.

Quantitative determination of decitabine incorporation into DNA and its effect on mutation rates in human cancer cells.

Decitabine (5-aza-2'-deoxycytidine) is a DNA methyltransferase inhibitor and an archetypal epigenetic drug for the therapy of myeloid leukemias. The mode of action of decitabine strictly depends on the incorporation of the drug into DNA. However, DNA incorporation and ensuing genotoxic effects of decitabine have not yet been investigated in human cancer cell lines or in models related to the approved indication of the drug. Here we describe a robust assay for the quantitative determination of decitabine incorporation rates into DNA from human cancer cells. Using a panel of human myeloid leukemia cell lines we show appreciable amounts of decitabine incorporation that closely correlated with cellular drug uptake. Decitabine incorporation was also detectable in primary cells from myeloid leukemia patients, indicating that the assay is suitable for biomarker analyses to predict drug responses in patients. Finally, we also used next-generation sequencing to comprehensively analyze the effects of decitabine incorporation on the DNA sequence level. Interestingly, this approach failed to reveal significant changes in the rates of point mutations and genome rearrangements in myeloid leukemia cell lines. These results indicate that standard rates of decitabine incorporation are not genotoxic in myeloid leukemia cells.

Loss of Tet enzymes compromises proper differentiation of embryonic stem cells.

Tet enzymes (Tet1/2/3) convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) and are dynamically expressed during development. Whereas loss of individual Tet enzymes or combined deficiency of Tet1/2 allows for embryogenesis, the effect of complete loss of Tet activity and 5hmC marks in development is not established. We have generated Tet1/2/3 triple-knockout (TKO) mouse embryonic stem cells (ESCs) and examined their developmental potential. Combined deficiency of all three Tets depleted 5hmC and impaired ESC differentiation, as seen in poorly differentiated TKO embryoid bodies (EBs) and teratomas. Consistent with impaired differentiation, TKO ESCs contributed poorly to chimeric embryos, a defect rescued by Tet1 reexpression, and could not support embryonic development. Global gene-expression and methylome analyses of TKO EBs revealed promoter hypermethylation and deregulation of genes implicated in embryonic development and differentiation. These findings suggest a requirement for Tet- and 5hmC-mediated DNA demethylation in proper regulation of gene expression during ESC differentiation and development.

Combined deficiency of Tet1 and Tet2 causes epigenetic abnormalities but is compatible with postnatal development.

Tet enzymes (Tet1/2/3) convert 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) in various embryonic and adult tissues. Mice mutant for either Tet1 or Tet2 are viable, raising the question of whether these enzymes have overlapping roles in development. Here we have generated Tet1 and Tet2 double-knockout (DKO) embryonic stem cells (ESCs) and mice. DKO ESCs remained pluripotent but were depleted of 5hmC and caused developmental defects in chimeric embryos. While a fraction of double-mutant embryos exhibited midgestation abnormalities with perinatal lethality, viable and overtly normal Tet1/Tet2-deficient mice were also obtained. DKO mice had reduced 5hmC and increased 5mC levels and abnormal methylation at various imprinted loci. Nevertheless, animals of both sexes were fertile, with females having smaller ovaries and reduced fertility. Our data show that loss of both enzymes is compatible with development but promotes hypermethylation and compromises imprinting. The data also suggest a significant contribution of Tet3 to hydroxylation of 5mC during development.

Characterization of genome methylation patterns in the desert locust Schistocerca gregaria.

DNA methylation is a widely conserved epigenetic modification. The analysis of genome-scale DNA methylation patterns in various organisms suggests that major features of animal methylomes are widely conserved. However, based on the variation of DNA methyltransferase genes in invertebrates, it has also been proposed that DNA methylation could provide a molecular mechanism for ecological adaptation. We have now analyzed the methylome of the desert locust, Schistocerca gregaria, which represents an organism with a high degree of phenotypic plasticity. Using genome-scale bisulfite sequencing, we show here that the S. gregaria methylome is characterized by CpG- and exon-specific methylation and thus shares two major features with other animal methylomes. In contrast to other invertebrates, however, overall methylation levels were substantially higher and a significant fraction of transposons was methylated. Additionally, genic sequences were densely methylated in a pronounced bimodal pattern, suggesting a role for DNA methylation in the regulation of locust gene expression. Our results thus uncover a unique pattern of genome methylation in locusts and provide an important foundation for investigating the role of DNA methylation in locust phase polyphenism.

Hydroxylation of 5-methylcytosine by TET2 maintains the active state of the mammalian HOXA cluster.

Differentiation is accompanied by extensive epigenomic reprogramming, leading to the repression of stemness factors and the transcriptional maintenance of activated lineage-specific genes. Here we use the mammalian Hoxa cluster of developmental genes as a model system to follow changes in DNA modification patterns during retinoic acid-induced differentiation. We find the inactive cluster to be marked by defined patterns of 5-methylcytosine (5mC). Upon the induction of differentiation, the active anterior part of the cluster becomes increasingly enriched in 5-hydroxymethylcytosine (5hmC), following closely the colinear activation pattern of the gene array, which is paralleled by the reduction of 5mC. Depletion of the 5hmC generating dioxygenase Tet2 impairs the maintenance of Hoxa activity and partially restores 5mC levels. Our results indicate that gene-specific 5mC-5hmC conversion by Tet2 is crucial for the maintenance of active chromatin states at lineage-specific loci.

Dnmt3a protects active chromosome domains against cancer-associated hypomethylation.

Changes in genomic DNA methylation patterns are generally assumed to play an important role in the etiology of human cancers. The Dnmt3a enzyme is required for the establishment of normal methylation patterns, and mutations in Dnmt3a have been described in leukemias. Deletion of Dnmt3a in a K-ras-dependent mouse lung cancer model has been shown to promote tumor progression, which suggested that the enzyme might suppress tumor development by stabilizing DNA methylation patterns. We have used whole-genome bisulfite sequencing to comprehensively characterize the methylomes from Dnmt3a wildtype and Dnmt3a-deficient mouse lung tumors. Our results show that profound global methylation changes can occur in K-ras-induced lung cancer. Dnmt3a wild-type tumors were characterized by large hypomethylated domains that correspond to nuclear lamina-associated domains. In contrast, Dnmt3a-deficient tumors showed a uniformly hypomethylated genome. Further data analysis revealed that Dnmt3a is required for efficient maintenance methylation of active chromosome domains and that Dnmt3a-deficient tumors show moderate levels of gene deregulation in these domains. In summary, our results uncover conserved features of cancer methylomes and define the role of Dnmt3a in maintaining DNA methylation patterns in cancer.

The head of Bartonella adhesin A is crucial for host cell interaction of Bartonella henselae.

Human pathogenic Bartonella henselae cause cat scratch disease and vasculoproliferative disorders (e.g. bacillary angiomatosis). Expression of Bartonella adhesin A (BadA) is crucial for bacterial autoagglutination, adhesion to host cells, binding to extracellular matrix proteins and proangiogenic reprogramming via activation of hypoxia inducible factor (HIF)-1. Like the prototypic Yersinia adhesin A, BadA belongs to the class of trimeric autotransporter adhesins and is constructed modularly consisting of a head, a long and repetitive neck-stalk module and a membrane anchor. Until now, the exact biological role of these domains is not known. Here, we analysed the function of the BadA head by truncating the repetitive neck-stalk module of BadA (B. henselae badA(-)/pHN23). Like B. henselae Marseille wild type, B. henselae badA(-)/pHN23 showed autoagglutination, adhesion to collagen and endothelial cells and activation of HIF-1 in host cells. Remarkably, B. henselae badA(-)/pHN23 did not bind to fibronectin (Fn) suggesting a crucial role of the deleted stalk domain in Fn binding. Additionally, the recombinantly expressed BadA head adhered to human umbilical vein endothelial cells and to a lesser degree to epithelial (HeLa 229) cells. Our data suggest that the head represents the major functional domain of BadA responsible for host adhesion and angiogenic reprogramming.

Complete genome sequence of the myxobacterium Sorangium cellulosum.

The genus Sorangium synthesizes approximately half of the secondary metabolites isolated from myxobacteria, including the anti-cancer metabolite epothilone. We report the complete genome sequence of the model Sorangium strain S. cellulosum So ce56, which produces several natural products and has morphological and physiological properties typical of the genus. The circular genome, comprising 13,033,779 base pairs, is the largest bacterial genome sequenced to date. No global synteny with the genome of Myxococcus xanthus is apparent, revealing an unanticipated level of divergence between these myxobacteria. A large percentage of the genome is devoted to regulation, particularly post-translational phosphorylation, which probably supports the strain's complex, social lifestyle. This regulatory network includes the highest number of eukaryotic protein kinase-like kinases discovered in any organism. Seventeen secondary metabolite loci are encoded in the genome, as well as many enzymes with potential utility in industry.

Analysis of Bartonella adhesin A expression reveals differences between various B. henselae strains.

Bartonella henselae causes cat scratch disease and the vasculoproliferative disorders bacillary angiomatosis and peliosis hepatis in humans. One of the best known pathogenicity factors of B. henselae is Bartonella adhesin A (BadA), which is modularly constructed, consisting of head, neck/stalk, and membrane anchor domains. BadA is important for the adhesion of B. henselae to extracellular-matrix proteins and endothelial cells (ECs). In this study, we analyzed different B. henselae strains for BadA expression, autoagglutination, fibronectin (Fn) binding, and adhesion to ECs. We found that the B. henselae strains Marseille, ATCC 49882, Freiburg 96BK3 (FR96BK3), FR96BK38, and G-5436 express BadA. Remarkably, BadA expression was lacking in a B. henselae ATCC 49882 variant, in strains ATCC 49793 and Berlin-1, and in the majority of bacteria of strain Berlin-2. Adherence of B. henselae to ECs and Fn reliably correlated with BadA expression. badA was present in all tested strains, although the length of the gene varied significantly due to length variations of the stalk region. Sequencing of the promoter, head, and membrane anchor regions revealed only minor differences that did not correlate with BadA expression, apart from strain Berlin-1, in which a 1-bp deletion led to a frameshift in the head region of BadA. Our data suggest that, apart from the identified genetic modifications (frameshift deletion and recombination), other so-far-unknown regulatory mechanisms influence BadA expression. Because of variations between and within different B. henselae isolates, BadA expression should be analyzed before performing infection experiments with B. henselae.