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.

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.

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.

Dnmt2-dependent methylomes lack defined DNA methylation patterns.

Several organisms have retained methyltransferase 2 (Dnmt2) as their only candidate DNA methyltransferase gene. However, information about Dnmt2-dependent methylation patterns has been limited to a few isolated loci and the results have been discussed controversially. In addition, recent studies have shown that Dnmt2 functions as a tRNA methyltransferase, which raised the possibility that Dnmt2-only genomes might be unmethylated. We have now used whole-genome bisulfite sequencing to analyze the methylomes of Dnmt2-only organisms at single-base resolution. Our results show that the genomes of Schistosoma mansoni and Drosophila melanogaster lack detectable DNA methylation patterns. Residual unconverted cytosine residues shared many attributes with bisulfite deamination artifacts and were observed at comparable levels in Dnmt2-deficient flies. Furthermore, genetically modified Dnmt2-only mouse embryonic stem cells lost the DNA methylation patterns found in wild-type cells. Our results thus uncover fundamental differences among animal methylomes and suggest that DNA methylation is dispensable for a considerable number of eukaryotic organisms.

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.

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.

Development of an epigenetic clock to predict visual age progression of human skin.

Aging is a complex process characterized by the gradual decline of physiological functions, leading to increased vulnerability to age-related diseases and reduced quality of life. Alterations in DNA methylation (DNAm) patterns have emerged as a fundamental characteristic of aged human skin, closely linked to the development of the well-known skin aging phenotype. These changes have been correlated with dysregulated gene expression and impaired tissue functionality. In particular, the skin, with its visible manifestations of aging, provides a unique model to study the aging process. Despite the importance of epigenetic age clocks in estimating biological age based on the correlation between methylation patterns and chronological age, a second-generation epigenetic age clock, which correlates DNAm patterns with a particular phenotype, specifically tailored to skin tissue is still lacking. In light of this gap, we aimed to develop a novel second-generation epigenetic age clock explicitly designed for skin tissue to facilitate a deeper understanding of the factors contributing to individual variations in age progression. To achieve this, we used methylation patterns from more than 370 female volunteers and developed the first skin-specific second-generation epigenetic age clock that accurately predicts the skin aging phenotype represented by wrinkle grade, visual facial age, and visual age progression, respectively. We then validated the performance of our clocks on independent datasets and demonstrated their broad applicability. In addition, we integrated gene expression and methylation data from independent studies to identify potential pathways contributing to skin age progression. Our results demonstrate that our epigenetic age clock, VisAgeX, specifically predicting visual age progression, not only captures known biological pathways associated with skin aging, but also adds novel pathways associated with skin aging.

Context-dependent DNA methylation signatures in animal livestock.

DNA methylation is an important epigenetic modification that is widely conserved across animal genomes. It is widely accepted that DNA methylation patterns can change in a context-dependent manner, including in response to changing environmental parameters. However, this phenomenon has not been analyzed in animal livestock yet, where it holds major potential for biomarker development. Building on the previous identification of population-specific DNA methylation in clonal marbled crayfish, we have now generated numerous base-resolution methylomes to analyze location-specific DNA methylation patterns. We also describe the time-dependent conversion of epigenetic signatures upon transfer from one environment to another. We further demonstrate production system-specific methylation signatures in shrimp, river-specific signatures in salmon and farm-specific signatures in chicken. Together, our findings provide a detailed resource for epigenetic variation in animal livestock and suggest the possibility for origin tracing of animal products by epigenetic fingerprinting.

The sequence of the pYV virulence plasmid from Yersinia enterocolitica strain WA-314 biogroup 1B serotype O:8.

Our laboratory strain Yersinia enterocolitica strain WA-314 biogroup 1B serotype O:8 displayed a different adhesion behavior to host cells compared to other Yersinia strains. To investigate whether this is based on differences in the gene content of the large pYV virulence plasmid which contains the major Yersinia adhesin YadA, we set out to sequence pYV(WA-314). pYV(WA-314) is very similar to pYV127/90, with a notable difference in the length of the Type III secretion system component YscP, which determines the needle length of the system. We found that we could improve the annotation of proteins previously described as "hypothetical" in pYV127/90 and other pYV plasmids, and show that pYV plasmids contain several and seemingly redundant plasmid partitioning and stabilization systems, explaining why these plasmids are not easily lost in laboratory cultures of Yersinia strains.

A chicken DNA methylation clock for the prediction of broiler health.

The domestic chicken (Gallus gallus domesticus) is the globally most important source of commercially produced meat. While genetic approaches have played an important role in the development of chicken stocks, little is known about chicken epigenetics. We have systematically analyzed the chicken DNA methylation machinery and DNA methylation landscape. While overall DNA methylation distribution was similar to mammals, sperm DNA appeared hypomethylated, which correlates with the absence of the DNMT3L cofactor in the chicken genome. Additional analysis revealed the presence of low-methylated regions, which are conserved gene regulatory elements that show tissue-specific methylation patterns. We also used whole-genome bisulfite sequencing to generate 56 single-base resolution methylomes from various tissues and developmental time points to establish an LMR-based DNA methylation clock for broiler chicken age prediction. This clock was used to demonstrate epigenetic age acceleration in animals with experimentally induced inflammation. Our study provides detailed insights into the chicken methylome and suggests a novel application of the DNA methylation clock as a marker for livestock health.

Location-Dependent DNA Methylation Signatures in a Clonal Invasive Crayfish.

DNA methylation is an important epigenetic modification that has been repeatedly implied in organismal adaptation. However, many previous studies that have linked DNA methylation patterns to environmental parameters have been limited by confounding factors, such as cell-type heterogeneity and genetic variation. In this study, we analyzed DNA methylation variation in marbled crayfish, a clonal and invasive freshwater crayfish that is characterized by a largely tissue-invariant methylome and negligible genetic variation. Using a capture-based subgenome bisulfite sequencing approach that covers a small, variably methylated portion of the marbled crayfish genome, we identified specific and highly localized DNA methylation signatures for specimens from geographically and ecologically distinct wild populations. These results were replicated both biologically and technically by re-sampling at different time points and by using independent methodology. Finally, we show specific methylation signatures for laboratory animals and for laboratory animals that were reared at a lower temperature. Our results thus demonstrate the existence of context-dependent DNA methylation signatures in a clonal animal.

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.

Single-cell transcriptomes of the human skin reveal age-related loss of fibroblast priming.

Fibroblasts are an essential cell population for human skin architecture and function. While fibroblast heterogeneity is well established, this phenomenon has not been analyzed systematically yet. We have used single-cell RNA sequencing to analyze the transcriptomes of more than 5,000 fibroblasts from a sun-protected area in healthy human donors. Our results define four main subpopulations that can be spatially localized and show differential secretory, mesenchymal and pro-inflammatory functional annotations. Importantly, we found that this fibroblast 'priming' becomes reduced with age. We also show that aging causes a substantial reduction in the predicted interactions between dermal fibroblasts and other skin cells, including undifferentiated keratinocytes at the dermal-epidermal junction. Our work thus provides evidence for a functional specialization of human dermal fibroblasts and identifies the partial loss of cellular identity as an important age-related change in the human dermis. These findings have important implications for understanding human skin aging and its associated phenotypes.

Epigenetic deregulation of lamina-associated domains in Hutchinson-Gilford progeria syndrome.

BACKGROUND: Hutchinson-Gilford progeria syndrome (HGPS) is a progeroid disease characterized by the early onset of age-related phenotypes including arthritis, loss of body fat and hair, and atherosclerosis. Cells from affected individuals express a mutant version of the nuclear envelope protein lamin A (termed progerin) and have previously been shown to exhibit prominent histone modification changes. METHODS: Here, we analyze the possibility that epigenetic deregulation of lamina-associated domains (LADs) is involved in the molecular pathology of HGPS. To do so, we studied chromatin accessibility (Assay for Transposase-accessible Chromatin (ATAC)-see/-seq), DNA methylation profiles (Infinium MethylationEPIC BeadChips), and transcriptomes (RNA-seq) of nine primary HGPS fibroblast cell lines and six additional controls, two parental and four age-matched healthy fibroblast cell lines. RESULTS: Our ATAC-see/-seq data demonstrate that primary dermal fibroblasts from HGPS patients exhibit chromatin accessibility changes that are enriched in LADs. Infinium MethylationEPIC BeadChip profiling further reveals that DNA methylation alterations observed in HGPS fibroblasts are similarly enriched in LADs and different from those occurring during healthy aging and Werner syndrome (WS), another premature aging disease. Moreover, HGPS patients can be stratified into two different subgroups according to their DNA methylation profiles. Finally, we show that the epigenetic deregulation of LADs is associated with HGPS-specific gene expression changes. CONCLUSIONS: Taken together, our results strongly implicate epigenetic deregulation of LADs as an important and previously unrecognized feature of HGPS, which contributes to disease-specific gene expression. Therefore, they not only add a new layer to the study of epigenetic changes in the progeroid syndrome, but also advance our understanding of the disease's pathology at the cellular level.

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.

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.

Who ate whom? Adaptive Helicobacter genomic changes that accompanied a host jump from early humans to large felines.

Helicobacter pylori infection of humans is so old that its population genetic structure reflects that of ancient human migrations. A closely related species, Helicobacter acinonychis, is specific for large felines, including cheetahs, lions, and tigers, whereas hosts more closely related to humans harbor more distantly related Helicobacter species. This observation suggests a jump between host species. But who ate whom and when did it happen? In order to resolve this question, we determined the genomic sequence of H. acinonychis strain Sheeba and compared it to genomes from H. pylori. The conserved core genes between the genomes are so similar that the host jump probably occurred within the last 200,000 (range 50,000-400,000) years. However, the Sheeba genome also possesses unique features that indicate the direction of the host jump, namely from early humans to cats. Sheeba possesses an unusually large number of highly fragmented genes, many encoding outer membrane proteins, which may have been destroyed in order to bypass deleterious responses from the feline host immune system. In addition, the few Sheeba-specific genes that were found include a cluster of genes encoding sialylation of the bacterial cell surface carbohydrates, which were imported by horizontal genetic exchange and might also help to evade host immune defenses. These results provide a genomic basis for elucidating molecular events that allow bacteria to adapt to novel animal hosts.

The methylome of the marbled crayfish links gene body methylation to stable expression of poorly accessible genes.

BACKGROUND: The parthenogenetic marbled crayfish (Procambarus virginalis) is a novel species that has rapidly invaded and colonized various different habitats. Adaptation to different environments appears to be independent of the selection of genetic variants, but epigenetic programming of the marbled crayfish genome remains to be understood. RESULTS: Here, we provide a comprehensive analysis of DNA methylation in marbled crayfish. Whole-genome bisulfite sequencing of multiple replicates and different tissues revealed a methylation pattern that is characterized by gene body methylation of housekeeping genes. Interestingly, this pattern was largely tissue invariant, suggesting a function that is unrelated to cell fate specification. Indeed, integrative analysis of DNA methylation, chromatin accessibility and mRNA expression patterns revealed that gene body methylation correlated with limited chromatin accessibility and stable gene expression, while low-methylated genes often resided in chromatin with higher accessibility and showed increased expression variation. Interestingly, marbled crayfish also showed reduced gene body methylation and higher gene expression variability when compared with their noninvasive mother species, Procambarus fallax. CONCLUSIONS: Our results provide novel insights into invertebrate gene body methylation and its potential role in adaptive gene regulation.