Methods and applications of genome-wide profiling of DNA damage and rare mutations.

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.

Deficiency of the Polycomb Protein RYBP and TET Methylcytosine Oxidases Promotes Extensive CpG Island Hypermethylation and Malignant Transformation.

Hypermethylation of CpG islands (CGI) is a common feature of cancer cells and predominantly affects Polycomb-associated genomic regions. Elucidating the underlying mechanisms leading to DNA hypermethylation in human cancer could help identify chemoprevention strategies. Here, we evaluated the role of Polycomb complexes and 5-methylcytosine (5mC) oxidases in protecting CGIs from DNA methylation and observed that four genes coding for components of Polycomb repressive complex 1 (PRC1) are downregulated in tumors. Inactivation of RYBP, a key activator of variant PRC1 complexes, in combination with all three 5mC oxidases (TET proteins) in nontumorigenic bronchial epithelial cells led to widespread hypermethylation of Polycomb-marked CGIs affecting almost 4,000 target genes, which closely resembled the DNA hypermethylation landscape observed in human squamous cell lung tumors. The RYBP- and TET-deficient cells showed methylation-associated aberrant regulation of cancer-relevant pathways, including defects in the Hippo tumor suppressor network. Notably, the quadruple knockout cells acquired a transformed phenotype, including anchorage-independent growth and formation of squamous cell carcinomas in mice. This work provides a mechanism promoting hypermethylation of CGIs and shows that such hypermethylation can lead to cell transformation. The breakdown of a two-pronged protection mechanism can be a route towards genome-wide hypermethylation of CGIs in tumors. SIGNIFICANCE: Dysfunction of the Polycomb component RYBP in combination with loss of 5-methylcytosine oxidases promotes widespread hypermethylation of CpG islands in bronchial cells and induces tumorigenesis, resembling changes seen in human lung tumors.

Circle Damage Sequencing for Whole-Genome Analysis of DNA Damage.

There are many sources of endogenous and exogenous DNA damage. Damaged bases represent a threat to genome integrity and may interfere with normal cellular processes such as replication and transcription. To understand the specificity and biological consequences of DNA damage, it is essential to employ methods that are sensitive enough to detect damaged DNA bases at the level of single nucleotide resolution and genome-wide. Here we describe in detail a method we developed for this purpose, circle damage sequencing (CD-seq). This method is based on the circularization of genomic DNA that contains damaged bases and conversion of the damaged sites into double-strand breaks using specific DNA repair enzymes. Library sequencing of the opened circles yields the precise positions of the DNA lesions that are present. CD-seq can be adopted to various types of DNA damage as long as a specific cleavage scheme can be designed.

TRIM28 secures skeletal stem cell fate during skeletogenesis by silencing neural gene expression and repressing GREM1/AKT/mTOR signaling axis.

Long bones are generated by mesoderm-derived skeletal progenitor/stem cells (SSCs) through endochondral ossification, a process of sequential chondrogenic and osteogenic differentiation tightly controlled by the synergy between intrinsic and microenvironment cues. Here, we report that loss of TRIM28, a transcriptional corepressor, in mesoderm-derived cells expands the SSC pool, weakens SSC osteochondrogenic potential, and endows SSCs with properties of ectoderm-derived neural crest cells (NCCs), leading to severe defects of skeletogenesis. TRIM28 preferentially enhances H3K9 trimethylation and DNA methylation on chromatin regions more accessible in NCCs; loss of this silencing upregulates neural gene expression and enhances neurogenic potential. Moreover, TRIM28 loss causes hyperexpression of GREM1, which is an extracellular signaling factor promoting SSC self-renewal and SSC neurogenic potential by activating AKT/mTORC1 signaling. Our results suggest that TRIM28-mediated chromatin silencing establishes a barrier for maintaining the SSC lineage trajectory and preventing a transition to ectodermal fate by regulating both intrinsic and microenvironment cues.

UVA Radiation, DNA Damage, and Melanoma.

Melanoma is a lethal type of skin tumor that has been linked with sunlight exposure chiefly in fair-skinned human populations. Wavelengths from the sun that can reach the earth's surface include UVA radiation (320-400 nm) and UVB radiation (280-320 nm). UVB effectively induces the formation of dimeric DNA photoproducts, preferentially the cyclobutane pyrimidine dimers (CPDs). The characteristic UVB signature mutations in the form of C to T mutations at dipyrimidine sequences are prevalent in melanoma tumor genomes and have been ascribed to deamination of cytosines within CPDs before DNA polymerase bypass. However, evidence from epidemiological, animal, and other experimental studies also suggest that UVA radiation may participate in melanoma formation. The DNA damage relevant for UVA includes specific types of CPDs at TT sequences and perhaps oxidative DNA damage to guanine, both induced by direct or indirect, photosensitization-mediated chemical and biophysical processes. We summarize the evidence for a potential role of UVA in melanoma and discuss some of the mechanistic pathways of how UVA may induce mutagenesis in melanocytes.

Z-DNA is remodelled by ZBTB43 in prospermatogonia to safeguard the germline genome and epigenome.

Mutagenic purine-pyrimidine repeats can adopt the left-handed Z-DNA conformation. DNA breaks at potential Z-DNA sites can lead to somatic mutations in cancer or to germline mutations that are transmitted to the next generation. It is not known whether any mechanism exists in the germ line to control Z-DNA structure and DNA breaks at purine-pyrimidine repeats. Here we provide genetic, epigenomic and biochemical evidence for the existence of a biological process that erases Z-DNA specifically in germ cells of the mouse male foetus. We show that a previously uncharacterized zinc finger protein, ZBTB43, binds to and removes Z-DNA, preventing the formation of DNA double-strand breaks. By removing Z-DNA, ZBTB43 also promotes de novo DNA methylation at CG-containing purine-pyrimidine repeats in prospermatogonia. Therefore, the genomic and epigenomic integrity of the species is safeguarded by remodelling DNA structure in the mammalian germ line during a critical window of germline epigenome reprogramming.

Concordance of hydrogen peroxide-induced 8-oxo-guanine patterns with two cancer mutation signatures of upper GI tract tumors.

Oxidative DNA damage has been linked to inflammation, cancer, and aging. Here, we have mapped two types of oxidative DNA damage, oxidized guanines produced by hydrogen peroxide and oxidized thymines created by potassium permanganate, at a single-base resolution. 8-Oxo-guanine occurs strictly dependent on the G/C sequence context and shows a pronounced peak at transcription start sites (TSSs). We determined the trinucleotide sequence pattern of guanine oxidation. This pattern shows high similarity to the cancer-associated single-base substitution signatures SBS18 and SBS36. SBS36 is found in colorectal cancers that carry mutations in MUTYH, encoding a repair enzyme that operates on 8-oxo-guanine mispairs. SBS18 is common in inflammation-associated upper gastrointestinal tract tumors including esophageal and gastric adenocarcinomas. Oxidized thymines induced by permanganate occur with a distinct dinucleotide specificity, 5'T-A/C, and are depleted at the TSS. Our data suggest that two cancer mutational signatures, SBS18 and SBS36, are caused by reactive oxygen species.

The major mechanism of melanoma mutations is based on deamination of cytosine in pyrimidine dimers as determined by circle damage sequencing.

Sunlight-associated melanomas carry a unique C-to-T mutation signature. UVB radiation induces cyclobutane pyrimidine dimers (CPDs) as the major form of DNA damage, but the mechanism of how CPDs cause mutations is unclear. To map CPDs at single-base resolution genome wide, we developed the circle damage sequencing (circle-damage-seq) method. In human cells, CPDs form preferentially in a tetranucleotide sequence context (5'-Py-T<>Py-T/A), but this alone does not explain the tumor mutation patterns. To test whether mutations arise at CPDs by cytosine deamination, we specifically mapped UVB-induced cytosine-deaminated CPDs. Transcription start sites (TSSs) were protected from CPDs and deaminated CPDs, but both lesions were enriched immediately upstream of the TSS, suggesting a mutation-promoting role of bound transcription factors. Most importantly, the genomic dinucleotide and trinucleotide sequence specificity of deaminated CPDs matched the prominent mutation signature of melanomas. Our data identify the cytosine-deaminated CPD as the leading premutagenic lesion responsible for mutations in melanomas.

Purification of TET Proteins.

The 5-methylcytosine (5mC) oxidation pathway mediated by TET proteins involves step-wise oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). 5fC and 5caC can be removed from DNA by base excision repair and the completion of this pathway results in "demethylation" of 5mC by converting the modified base back into cytosine. In vitro studies with TET proteins aimed at analyzing their DNA substrate specificities and their activity within defined chromatin templates are relatively limited. Here we describe purification methods for mammalian TET proteins based on expression in insect cells or in 293T cells. We also briefly summarize a method that can be used to monitor 5-methylcytosine oxidase activity of the purified TET proteins in vitro.

Reprogramming of DNA methylation at NEUROD2-bound sequences during cortical neuron differentiation.

The characteristics of DNA methylation changes that occur during neurogenesis in vivo remain unknown. We used whole-genome bisulfite sequencing to quantitate DNA cytosine modifications in differentiating neurons and their progenitors isolated from mouse brain at the peak of embryonic neurogenesis. Localized DNA hypomethylation was much more common than hypermethylation and often occurred at putative enhancers within genes that were upregulated in neurons and encoded proteins crucial for neuronal differentiation. The hypomethylated regions strongly overlapped with mapped binding sites of the key neuronal transcription factor NEUROD2. The 5-methylcytosine oxidase ten-eleven translocation 2 (TET2) interacted with NEUROD2, and its reaction product 5-hydroxymethylcytosine accumulated at the demethylated regions. NEUROD2-targeted differentially methylated regions retained higher methylation levels in Neurod2 knockout mice, and inducible expression of NEUROD2 caused TET2-associated demethylation at its in vivo binding sites. The data suggest that the reorganization of DNA methylation in developing neurons involves NEUROD2 and TET2-mediated DNA demethylation.

An Intrinsic Epigenetic Barrier for Functional Axon Regeneration.

Mature neurons in the adult peripheral nervous system can effectively switch from a dormant state with little axonal growth to robust axon regeneration upon injury. The mechanisms by which injury unlocks mature neurons' intrinsic axonal growth competence are not well understood. Here, we show that peripheral sciatic nerve lesion in adult mice leads to elevated levels of Tet3 and 5-hydroxylmethylcytosine in dorsal root ganglion (DRG) neurons. Functionally, Tet3 is required for robust axon regeneration of DRG neurons and behavioral recovery. Mechanistically, peripheral nerve injury induces DNA demethylation and upregulation of multiple regeneration-associated genes in a Tet3- and thymine DNA glycosylase-dependent fashion in DRG neurons. In addition, Pten deletion-induced axon regeneration of retinal ganglion neurons in the adult CNS is attenuated upon Tet1 knockdown. Together, our study suggests an epigenetic barrier that can be removed by active DNA demethylation to permit axon regeneration in the adult mammalian nervous system.

Tet3 Reads 5-Carboxylcytosine through Its CXXC Domain and Is a Potential Guardian against Neurodegeneration.

We report that the mammalian 5-methylcytosine (5mC) oxidase Tet3 exists as three major isoforms and characterized the full-length isoform containing an N-terminal CXXC domain (Tet3FL). This CXXC domain binds to unmethylated CpGs, but, unexpectedly, its highest affinity is toward 5-carboxylcytosine (5caC). We determined the crystal structure of the CXXC domain-5caC-DNA complex, revealing the structural basis of the binding specificity of this domain as a reader of CcaCG sequences. Mapping of Tet3FL in neuronal cells shows that Tet3FL is localized precisely at the transcription start sites (TSSs) of genes involved in lysosome function, mRNA processing, and key genes of the base excision repair pathway. Therefore, Tet3FL may function as a regulator of 5caC removal by base excision repair. Active removal of accumulating 5mC from the TSSs of genes coding for lysosomal proteins by Tet3FL in postmitotic neurons of the brain may be important for preventing neurodegenerative diseases.

The DNA methylation landscape of human melanoma.

Using MIRA-seq, we have characterized the DNA methylome of metastatic melanoma and normal melanocytes. Individual tumors contained several thousand hypermethylated regions. We discovered 179 tumor-specific methylation peaks present in all (27/27) melanomas that may be effective disease biomarkers, and 3113 methylation peaks were seen in >40% of the tumors. We found that 150 of the approximately 1200 tumor-associated methylation peaks near transcription start sites (TSSs) were marked by H3K27me3 in melanocytes. DNA methylation in melanoma was specific for distinct H3K27me3 peaks rather than for broadly covered regions. However, numerous H3K27me3 peak-associated TSS regions remained devoid of DNA methylation in tumors. There was no relationship between BRAF mutations and the number of methylation peaks. Gene expression analysis showed upregulated immune response genes in melanomas presumably as a result of lymphocyte infiltration. Down-regulated genes were enriched for melanocyte differentiation factors; e.g., KIT, PAX3 and SOX10 became methylated and downregulated in melanoma.

Longitudinal epigenetic and gene expression profiles analyzed by three-component analysis reveal down-regulation of genes involved in protein translation in human aging.

Data on biological mechanisms of aging are mostly obtained from cross-sectional study designs. An inherent disadvantage of this design is that inter-individual differences can mask small but biologically significant age-dependent changes. A serially sampled design (same individual at different time points) would overcome this problem but is often limited by the relatively small numbers of available paired samples and the statistics being used. To overcome these limitations, we have developed a new vector-based approach, termed three-component analysis, which incorporates temporal distance, signal intensity and variance into one single score for gene ranking and is combined with gene set enrichment analysis. We tested our method on a unique age-based sample set of human skin fibroblasts and combined genome-wide transcription, DNA methylation and histone methylation (H3K4me3 and H3K27me3) data. Importantly, our method can now for the first time demonstrate a clear age-dependent decrease in expression of genes coding for proteins involved in translation and ribosome function. Using analogies with data from lower organisms, we propose a model where age-dependent down-regulation of protein translation-related components contributes to extend human lifespan.

MIRA-seq for DNA methylation analysis of CpG islands.

AIM: To develop a reliable method for whole genome analysis of DNA methylation. MATERIALS & METHODS: Genome-scale analysis of DNA methylation includes affinity-based approaches such as enrichment using methyl-CpG-binding proteins. One of these methods, the methylated-CpG island recovery assay (MIRA), is based on the high affinity of the MBD2b-MBD3L1 complex for CpG-methylated DNA. Here we provide a detailed description of MIRA and combine it with next generation sequencing platforms (MIRA-seq). RESULTS: We assessed the performance of MIRA-seq and compared the data with whole genome bisulfite sequencing. CONCLUSION: MIRA-seq is a reliable, genome-scale DNA methylation analysis platform for scoring DNA methylation differences at CpG-rich genomic regions. The method is not limited by primer or probe design and is cost effective.

Tet-mediated formation of 5-hydroxymethylcytosine in RNA.

Oxidation of 5-methylcytosine in DNA by ten-eleven translocation (Tet) family of enzymes has been demonstrated to play a significant role in epigenetic regulation in mammals. We found that Tet enzymes also possess the activity of catalyzing the formation of 5-hydroxymethylcytidine (5-hmrC) in RNA in vitro. In addition, the catalytic domains of all three Tet enzymes as well as full-length Tet3 could induce the formation of 5-hmrC in human cells. Moreover, 5-hmrC was present at appreciable levels (∼1 per 5000 5-methylcytidine) in RNA of mammalian cells and tissues. Our results suggest the involvement of this oxidation in RNA biology.

The role of 5-hydroxymethylcytosine in human cancer.

The patterns of DNA methylation in human cancer cells are highly abnormal and often involve the acquisition of DNA hypermethylation at hundreds or thousands of CpG islands that are usually unmethylated in normal tissues. The recent discovery of 5-hydroxymethylcytosine (5hmC) as an enzymatic oxidation product of 5-methylcytosine (5mC) has led to models and experimental data in which the hypermethylation and 5mC oxidation pathways seem to be connected. Key discoveries in this setting include the findings that several genes coding for proteins involved in the 5mC oxidation reaction are mutated in human tumors, and that a broad loss of 5hmC occurs across many types of cancer. In this review, we will summarize current knowledge and discuss models of the potential roles of 5hmC in human cancer biology.

Formation of cyclobutane pyrimidine dimers at dipyrimidines containing 5-hydroxymethylcytosine.

Much of the cancer-causing effects of ultraviolet radiation from the sun have been linked to the formation of dimerized DNA bases. These dimeric DNA photoproducts include the cyclobutane pyrimidine dimers (CPDs) and the pyrimidine(6-4)pyrimidone photoproducts [(6-4)PPs]. CPDs are highly mutagenic and are produced in substantial quantities by UVB radiation. These dimers can form between any two adjacent pyrimidines and can involve thymine, cytosine, or 5-methylcytosine. Very recently, a sixth DNA base, 5-hydroxymethylcytosine (5hmC) has been identified and characterized as a normal component of mammalian DNA. Here, we investigated the formation of CPDs at different DNA sequences containing 5hmC following irradiation with UVA, UVB, or UVC light sources. We show that the formation of CPDs at dipyrimidines containing 5hmC occurs at different DNA sequences but is not enhanced relative to cytosine or 5-methylcytosines at the same sequence positions. In fact, in some sequence contexts, CPDs containing 5hmC are formed at very low levels. Nonetheless, CPD formation at 5hmC pyrimidines is expected to be biologically relevant since three types of human skin-derived cells, fibroblasts, keratinocytes and melanocytes, all contain detectable levels of this modified base.

5-hydroxymethylcytosine and its potential roles in development and cancer.

Only a few years ago it was demonstrated that mammalian DNA contains oxidized forms of 5-methylcytosine (5mC). The base 5-hydroxymethylcytosine (5hmC) is the most abundant of these oxidation products and is referred to as the sixth DNA base. 5hmC is produced from 5mC in an enzymatic pathway involving three 5mC oxidases, Ten-eleven translocation (TET)1, TET2, and TET3. The biological role of 5hmC is still unclear. Current models propose that 5hmC is an intermediate base in an active or passive DNA demethylation process that operates during important reprogramming phases of mammalian development. Tumors originating in various human tissues have strongly depleted levels of 5hmC. Apparently, 5hmC cannot be maintained in proliferating cells. Furthermore, mutations in the TET2 gene are commonly observed in human myeloid malignancies. Since TET proteins and many lysine demethylases require 2-oxoglutarate as a cofactor, aberrations in cofactor biochemical pathways, including mutations in isocitrate dehydrogenase (IDH), may affect levels of 5hmC and 5mC in certain types of tumors, either directly or indirectly. We discuss current data and models of the function of 5hmC in general, with special emphasis on its role in mechanisms of development and cancer.

Dynamics of 5-hydroxymethylcytosine and chromatin marks in Mammalian neurogenesis.

DNA methylation in mammals is highly dynamic during germ cell and preimplantation development but is relatively static during the development of somatic tissues. 5-hydroxymethylcytosine (5hmC), created by oxidation of 5-methylcytosine (5mC) by Tet proteins and most abundant in the brain, is thought to be an intermediary toward 5mC demethylation. We investigated patterns of 5mC and 5hmC during neurogenesis in the embryonic mouse brain. 5hmC levels increase during neuronal differentiation. In neuronal cells, 5hmC is not enriched at enhancers but associates preferentially with gene bodies of activated neuronal function-related genes. Within these genes, gain of 5hmC is often accompanied by loss of H3K27me3. Enrichment of 5hmC is not associated with substantial DNA demethylation, suggesting that 5hmC is a stable epigenetic mark. Functional perturbation of the H3K27 methyltransferase Ezh2 or of Tet2 and Tet3 leads to defects in neuronal differentiation, suggesting that formation of 5hmC and loss of H3K27me3 cooperate to promote brain development.