Non-canonical functions of UHRF1 maintain DNA methylation homeostasis in cancer cells.

DNA methylation is an essential epigenetic chromatin modification, and its maintenance in mammals requires the protein UHRF1. It is yet unclear if UHRF1 functions solely by stimulating DNA methylation maintenance by DNMT1, or if it has important additional functions. Using degron alleles, we show that UHRF1 depletion causes a much greater loss of DNA methylation than DNMT1 depletion. This is not caused by passive demethylation as UHRF1-depleted cells proliferate more slowly than DNMT1-depleted cells. Instead, bioinformatics, proteomics and genetics experiments establish that UHRF1, besides activating DNMT1, interacts with DNMT3A and DNMT3B and promotes their activity. In addition, we show that UHRF1 antagonizes active DNA demethylation by TET2. Therefore, UHRF1 has non-canonical roles that contribute importantly to DNA methylation homeostasis; these findings have practical implications for epigenetics in health and disease.

A highly efficient scheme for library preparation from single-stranded DNA.

Although methods for sequencing library preparation from double-stranded DNA are well established, those from single-stranded DNA (ssDNA) have not been well studied. Further, the existing methods have limitations in efficiency and yield. Therefore, we developed a highly efficient procedure for sequencing library preparation from ssDNA. In this method, the first adaptor tagging of ssDNA is performed using terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation, which we reported recently. After complementary strand synthesis using the adaptor-tagged ssDNA, second adaptor tagging via Vaccinia virus topoisomerase I (VTopoI or TOPO)-based adaptor ligation is performed. With additional steps for degradation, repression, and removal of the adaptor dimer, the proposed TACS-TOPO scheme realizes adaptor dimer-free sequencing library preparation from ssDNA samples of 24 pg. The TACS-TOPO scheme was successfully applied to cell-free DNA analysis with amplification-free library preparation from 50 µL of human serum. A modified TACS-TOPO scheme was also applied to DNA extracted from ancient human bones, bringing two to eight times more library yields than those using a conventional library preparation protocol. The procedures for preparing VTopoI and its complex with a double-stranded oligonucleotide adaptor are also described. Overall, the proposed TACS-TOPO scheme can facilitate practical and sensitive sequencing analysis of ssDNA.

A genetic screen identifies BEND3 as a regulator of bivalent gene expression and global DNA methylation.

Epigenetic mechanisms are essential to establish and safeguard cellular identities in mammals. They dynamically regulate the expression of genes, transposable elements and higher-order chromatin structures. Consequently, these chromatin marks are indispensable for mammalian development and alterations often lead to disease, such as cancer. Bivalent promoters are especially important during differentiation and development. Here we used a genetic screen to identify new regulators of a bivalent repressed gene. We identify BEND3 as a regulator of hundreds of bivalent promoters, some of which it represses, and some of which it activates. We show that BEND3 is recruited to a CpG-containg consensus site that is present in multiple copies in many bivalent promoters. Besides having direct effect on the promoters it binds, the loss of BEND3 leads to genome-wide gains of DNA methylation, which are especially marked at regions normally protected by the TET enzymes. DNA hydroxymethylation is reduced in Bend3 mutant cells, possibly as consequence of altered gene expression leading to diminished alpha-ketoglutarate production, thus lowering TET activity. Our results clarify the direct and indirect roles of an important chromatin regulator, BEND3, and, more broadly, they shed light on the regulation of bivalent promoters.

Crucial role of iron in epigenetic rewriting during adipocyte differentiation mediated by JMJD1A and TET2 activity.

Iron metabolism is closely associated with the pathogenesis of obesity. However, the mechanism of the iron-dependent regulation of adipocyte differentiation remains unclear. Here, we show that iron is essential for rewriting of epigenetic marks during adipocyte differentiation. Iron supply through lysosome-mediated ferritinophagy was found to be crucial during the early stage of adipocyte differentiation, and iron deficiency during this period suppressed subsequent terminal differentiation. This was associated with demethylation of both repressive histone marks and DNA in the genomic regions of adipocyte differentiation-associated genes, 　including Pparg, which encodes PPARγ, the master regulator of adipocyte differentiation. In addition, we identified several epigenetic demethylases to be responsible for iron-dependent adipocyte differentiation, with the histone demethylase jumonji domain-containing 1A and the DNA demethylase ten-eleven translocation 2 as the major enzymes. The interrelationship between repressive histone marks and DNA methylation was indicated by an integrated genome-wide association analysis, and was also supported by the findings that both histone and DNA demethylation were suppressed by either the inhibition of lysosomal ferritin flux or the knockdown of iron chaperone poly(rC)-binding protein 2. In summary, epigenetic regulations through iron-dependent control of epigenetic enzyme activities play an important role in the organized gene expression mechanisms of adipogenesis.

Persistent epigenetic alterations in transcription factors after a sustained virological response in hepatocellular carcinoma.

Background and Aim: The risk of hepatocellular carcinoma (HCC) persists in a condition of sustained virologic response (SVR) after hepatitis C virus (HCV) eradication. Comprehensive molecular analyses were performed to test the hypothesis that epigenetic abnormalities present after an SVR play a role in hepatocarcinogenesis. Methods: Whole-genome methylome and RNA sequencing were performed on HCV, SVR, and healthy liver tissue. Integrated analysis of the sequencing data focused on expression changes in transcription factors and their target genes, commonly found in HCV and SVR. Identified expression changes were validated in demethylated cultured HCC cell lines and an independent validation cohort. Results: The coincidence rates of the differentially methylated regions between the HCV and SVR groups were 91% in the hypomethylated and 71% in the hypermethylated regions in tumorous tissues, and 37% in the hypomethylated and 36% in the hypermethylated regions in non-tumorous tissues. These results indicate that many epigenomic abnormalities persist even after an SVR was achieved. Integrated analysis identified 61 transcription factors and 379 other genes that had methylation abnormalities and gene expression changes in both groups. Validation cohort specified gene expression changes for 14 genes, and gene ontology pathway analysis revealed apoptotic signaling and inflammatory response were associated with these genes. Conclusion: This study demonstrates that DNA methylation abnormalities, retained after HCV eradication, affect the expression of transcription factors and their target genes. These findings suggest that DNA methylation in SVR patients may be functionally important in carcinogenesis, and could serve as biomarkers to predict HCC occurrence.

Clonal germinal center B cells function as a niche for T-cell lymphoma.

Angioimmunoblastic T-cell lymphoma (AITL) is proposed to be initiated by age-related clonal hematopoiesis (ACH) with TET2 mutations, whereas the G17V RHOA mutation in immature cells with TET2 mutations promotes the development of T follicular helper (TFH)-like tumor cells. Here, we investigated the mechanism by which TET2-mutant immune cells enable AITL development using mouse models and human samples. Among the 2 mouse models, mice lacking Tet2 in all the blood cells (Mx-Cre × Tet2flox/flox × G17V RHOA transgenic mice) spontaneously developed AITL for approximately up to a year, while mice lacking Tet2 only in the T cells (Cd4-Cre × Tet2flox/flox × G17V RHOA transgenic mice) did not. Therefore, Tet2-deficient immune cells function as a niche for AITL development. Single-cell RNA-sequencing (scRNA-seq) of >50 000 cells from mouse and human AITL samples revealed significant expansion of aberrant B cells, exhibiting properties of activating light zone (LZ)-like and proliferative dark zone (DZ)-like germinal center B (GCB) cells. The GCB cells in AITL clonally evolved with recurrent mutations in genes related to core histones. In silico network analysis using scRNA-seq data identified Cd40-Cd40lg as a possible mediator of GCB and tumor cell cluster interactions. Treatment of AITL model mice with anti-Cd40lg inhibitory antibody prolonged survival. The genes expressed in aberrantly expanded GCB cells in murine tumors were also broadly expressed in the B-lineage cells of TET2-mutant human AITL. Therefore, ACH-derived GCB cells could undergo independent clonal evolution and support the tumorigenesis in AITL via the CD40-CD40LG axis.

Pan-cancer methylome analysis for cancer diagnosis and classification of cancer cell of origin.

The accurate and early diagnosis and classification of cancer origin from either tissue or liquid biopsy is crucial for selecting the appropriate treatment and reducing cancer-related mortality. Here, we established the CAncer Cell-of-Origin (CACO) methylation panel using the methylation data of the 28 types of cancer in The Cancer Genome Atlas (7950 patients and 707 normal controls) as well as healthy whole blood samples (95 subjects). We showed that the CACO methylation panel had high diagnostic potential with high sensitivity and specificity in the discovery (maximum AUC = 0.998) and validation (maximum AUC = 1.000) cohorts. Moreover, we confirmed that the CACO methylation panel could identify the cancer cell type of origin using the methylation profile from liquid as well as tissue biopsy, including primary, metastatic, and multiregional cancer samples and cancer of unknown primary, independent of the methylation analysis platform and specimen preparation method. Together, the CACO methylation panel can be a powerful tool for the classification and diagnosis of cancer.

The Evolving Genomic Landscape of Esophageal Squamous Cell Carcinoma Under Chemoradiotherapy.

Esophageal squamous cell carcinoma (ESCC) often recurs after chemoradiotherapy, and the prognosis of ESCC after chemoradiotherapy has not improved over the past few decades. The mutation process in chemoradiotherapy-resistant clones and the functional relevance of genetic alterations remain unclear. To address these problems, we performed whole-exome sequencing of 52 tumor samples from 33 patients with ESCC who received radiotherapy combined with 5-fluorouracil/platinum. In multiregion analyses of pretreatment and locally recurrent lesions from five cases, most driver gene-altered clones remained under chemoradiotherapy selection pressure, while few driver gene alterations were acquired at recurrence. The mutation signatures of recurrent ESCC, including increased deletion frequency and platinum dose-dependent base substitution signatures, were substantially different from those of primary ESCC and reflected the iatrogenic impacts of chemoradiotherapy. Single-region analysis of 28 pretreatment tumors indicated that focal copy-number gain at the MYC locus was significantly associated with poor progression-free survival and overall survival after chemoradiotherapy. MYC gain remained throughout the chemoradiotherapy course and potentially contributes to intrinsic resistance to chemoradiotherapy. Consistent with these findings, MYC copy number and mRNA and protein levels in ESCC cell lines correlated positively with resistance to radiotherapy, and MYC knockdown improved sensitivity to radiotherapy. Overall, these data characterize the clonal evolution process induced by chemoradiotherapy and clinically relevant associations for genetic alterations in ESCC. These findings increase our understanding of therapeutic resistance and support the rationale for precision chemoradiotherapy. SIGNIFICANCE: Whole-exome sequencing reveals the genetic evolution of ESCC during chemoradiotherapy, highlighting MYC gain in pretreatment tumors as a potential marker of therapy resistance.

Limited rejuvenation of aged hematopoietic stem cells in young bone marrow niche.

Hematopoietic stem cells (HSCs) exhibit functional alterations, such as reduced regenerative capacity and myeloid-biased differentiation, with age. The HSC niche, which is essential for the maintenance of HSCs, also undergoes marked changes with aging. However, it has been technically challenging to directly evaluate the contribution of niche aging to age-associated HSC alterations without niche-damaging myeloablation in HSC transplantation assays. We herein transplanted an excess of aged HSCs into young mice without preconditioning. Although aged HSCs successfully engrafted in the intact young bone marrow niche, they poorly regenerated downstream progenitors and exhibited persistent myeloid-biased differentiation, resulting in no significant functional rejuvenation. Transcriptome and methylome analyses revealed that the young niche largely restored the transcriptional profile of aged HSCs, but not their DNA methylation profiles. Therefore, the restoration of the young niche is insufficient for rejuvenating HSC functions, highlighting a key role for age-associated cell-intrinsic defects in HSC aging.

Base-resolution methylomes of gliomas bearing histone H3.3 mutations reveal a G34 mutant-specific signature shared with bone tumors.

Two recurrent mutations, K27M and G34R/V, in H3F3A, encoding non-canonical histone H3.3, are reported in pediatric and young adult gliomas, whereas G34W mutation is prevalent in bone tumors. In contrast to K27M mutation, it remains elusive how G34 mutations affect the epigenome. Here we performed whole-genome bisulfite sequencing of four G34R-mutated gliomas and the G34V-mutated glioma cell line KNS-42 for comparison with gliomas harboring K27M and no mutations in H3F3A and with G34W-mutated bone tumors. G34R-mutated gliomas exhibited lower global methylation levels, similar CpG island (CGI) methylation levels, and compromised hypermethylation of telomere-proximal CGIs, compared to the other two glioma subgroups. Hypermethylated regions specific to G34R-mutated gliomas were enriched for CGIs, including those of OLIG1, OLIG2, and canonical histone genes in the HIST1 cluster. They were notably hypermethylated in osteosarcomas with, but not without, G34W mutation. Independent component analysis revealed that G34 mutation-specific components shared a significant similarity between glioma and osteosarcoma, suggesting that G34 mutations exert characteristic methylomic effects regardless of the tumor tissue-of-origin. CRISPR/Cas9-mediated disruption of G34V-allele in KNS-42 cells led to demethylation of a subset of CGIs hypermethylated in G34R-mutated gliomas. These findings will provide a basis for elucidating epigenomic roles of G34 oncohistone in tumorigenesis.

Tet2 and Tet3 in B cells are required to repress CD86 and prevent autoimmunity.

A contribution of epigenetic modifications to B cell tolerance has been proposed but not directly tested. Here we report that deficiency of ten-eleven translocation (Tet) DNA demethylase family members Tet2 and Tet3 in B cells led to hyperactivation of B and T cells, autoantibody production and lupus-like disease in mice. Mechanistically, in the absence of Tet2 and Tet3, downregulation of CD86, which normally occurs following chronic exposure of self-reactive B cells to self-antigen, did not take place. The importance of dysregulated CD86 expression in Tet2- and Tet3-deficient B cells was further demonstrated by the restriction, albeit not complete, on aberrant T and B cell activation following anti-CD86 blockade. Tet2- and Tet3-deficient B cells had decreased accumulation of histone deacetylase 1 (HDAC1) and HDAC2 at the Cd86 locus. Thus, our findings suggest that Tet2- and Tet3-mediated chromatin modification participates in repression of CD86 on chronically stimulated self-reactive B cells, which contributes, at least in part, to preventing autoimmunity.

Base-Resolution Methylome of Retinal Pigment Epithelial Cells Used in the First Trial of Human Induced Pluripotent Stem Cell-Based Autologous Transplantation.

The first-in-human trial of induced pluripotent stem cell (iPSC)-based autologous transplantation was successfully performed on a female patient with age-related macular degeneration. Here we delineated the base-resolution methylome of the iPSC-derived retinal pigment epithelium (iRPE) used in this trial. The methylome of iRPE closely resembled that of native RPE (nRPE), although partially methylated domains (PMDs) emerged in iRPE but not nRPE. Most differentially methylated regions between iRPE and nRPE appeared to originate from (de)methylation errors during differentiation, whereas errors at reprogramming resulted in aberrant genomic imprinting and X chromosome reactivation. Moreover, non-CpG methylation was prominent in nRPE but not iRPE. Intriguingly, xenotransplantation to mouse remodeled the iRPE methylome to demethylate a subset of suppressed genes and accumulate non-CpG methylation, but failed to resolve PMDs and hypermethylated CpG islands. Although the impacts of these alterations remain elusive, our findings should provide a useful guide for methylome analyses of other iPSC-derived cells.

Highly efficient single-stranded DNA ligation technique improves low-input whole-genome bisulfite sequencing by post-bisulfite adaptor tagging.

Whole-genome bisulfite sequencing (WGBS) is the current gold standard of methylome analysis. Post-bisulfite adaptor tagging (PBAT) is an increasingly popular WGBS protocol because of high sensitivity and low bias. PBAT originally relied on two rounds of random priming for adaptor-tagging of single-stranded DNA (ssDNA) to attain high efficiency but at a cost of library insert length. To overcome this limitation, we developed terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation as an alternative to random priming. In this method, TdT attaches adenylates to the 3'-end of input ssDNA, which are then utilized by RNA ligase as an efficient connector to the ssDNA adaptor. A protocol that uses TACS ligation instead of the second random priming step substantially increased the lengths of PBAT library fragments. Moreover, we devised a dual-library strategy that splits the input DNA to prepare two libraries with reciprocal adaptor polarity, combining them prior to sequencing. This strategy ensured an ideal base-color balance to eliminate the need for DNA spike-in for color compensation, further improving the throughput and quality of WGBS. Adopting the above strategies to the HiSeq X Ten and NovaSeq 6000 platforms, we established a cost-effective, high-quality WGBS, which should accelerate various methylome analyses.

Epigenetic landscape influences the liver cancer genome architecture.

The accumulations of different types of genetic alterations such as nucleotide substitutions, structural rearrangements and viral genome integrations and epigenetic alterations contribute to carcinogenesis. Here, we report correlation between the occurrence of epigenetic features and genetic aberrations by whole-genome bisulfite, whole-genome shotgun, long-read, and virus capture sequencing of 373 liver cancers. Somatic substitutions and rearrangement breakpoints are enriched in tumor-specific hypo-methylated regions with inactive chromatin marks and actively transcribed highly methylated regions in the cancer genome. Individual mutation signatures depend on chromatin status, especially, signatures with a higher transcriptional strand bias occur within active chromatic areas. Hepatitis B virus (HBV) integration sites are frequently detected within inactive chromatin regions in cancer cells, as a consequence of negative selection for integrations in active chromatin regions. Ultra-high structural instability and preserved unmethylation of integrated HBV genomes are observed. We conclude that both precancerous and somatic epigenetic features contribute to the cancer genome architecture.

The AP-1 transcription factor JunB is required for Th17 cell differentiation.

Interleukin (IL)-17-producing T helper (Th17) cells are crucial for host defense against extracellular microbes and pathogenesis of autoimmune diseases. Here we show that the AP-1 transcription factor JunB is required for Th17 cell development. Junb-deficient CD4+ T cells are able to develop in vitro into various helper T subsets except Th17. The RNA-seq transcriptome analysis reveals that JunB is crucial for the Th17-specific gene expression program. Junb-deficient mice are completely resistant to experimental autoimmune encephalomyelitis, a Th17-mediated inflammatory disease, and naive T helper cells from such mice fail to differentiate into Th17 cells. JunB appears to activate Th17 signature genes by forming a heterodimer with BATF, another AP-1 factor essential for Th17 differentiation. The mechanism whereby JunB controls Th17 cell development likely involves activation of the genes for the Th17 lineage-specifying orphan receptors RORγt and RORα and reduced expression of Foxp3, a transcription factor known to antagonize RORγt function.

Changepoint detection in base-resolution methylome data reveals a robust signature of methylated domain landscape.

BACKGROUND: Base-resolution methylome data generated by whole-genome bisulfite sequencing (WGBS) is often used to segment the genome into domains with distinct methylation levels. However, most segmentation methods include many parameters to be carefully tuned and/or fail to exploit the unsurpassed resolution of the data. Furthermore, there is no simple method that displays the composition of the domains to grasp global trends in each methylome. RESULTS: We propose to use changepoint detection for domain demarcation based on base-resolution methylome data. While the proposed method segments the methylome in a largely comparable manner to conventional approaches, it has only a single parameter to be tuned. Furthermore, it fully exploits the base-resolution of the data to enable simultaneous detection of methylation changes in even contrasting size ranges, such as focal hypermethylation and global hypomethylation in cancer methylomes. We also propose a simple plot termed methylated domain landscape (MDL) that globally displays the size, the methylation level and the number of the domains thus defined, thereby enabling one to intuitively grasp trends in each methylome. Since the pattern of MDL often reflects cell lineages and is largely unaffected by data size, it can serve as a novel signature of methylome. CONCLUSIONS: Changepoint detection in base-resolution methylome data followed by MDL plotting provides a novel method for methylome characterization and will facilitate global comparison among various WGBS data differing in size and even species origin.

A comprehensive analysis of allelic methylation status of CpG islands on human chromosome 21q.

Approximately half of all human genes have CpG islands (CGIs)around their promoter regions. Although CGIs usually escape methylation, those on Chromosome X in females and those in the vicinity of imprinted genes are exceptions: They have both methylated and unmethylated alleles to display a "composite" pattern in methylation analysis. In addition, aberrant methylation of CGIs is known to often occur in cancer cells. Here we developed a simple HpaII-McrBC PCR method for discrimination of full, null, incomplete, and composite methylation patterns, and applied it to all computationally identified CGIs on human Chromosome 21q. This comprehensive analysis revealed that, although most CGIs (103 out of 149)escape methylation, a sizable fraction (31 out of 149)are fully methylated even in normal peripheral blood cells. Furthermore, we identified seven CGIs showing the composite methylation, and demonstrated that three of them are indeed methylated monoallelically. Further analyses using informative pedigrees revealed that two of the three are subject to maternal allele-specific methylation. Intriguingly, the other CGI is methylated in an allele-specific but parental-origin-independent manner. Thus, the cell seems to have a broader repertoire of methylating CGIs than previously thought, and our approach may contribute to uncover novel modes of allelic methylation.