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

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Loop-extruding Smc5/6 organizes transcription-induced positive DNA supercoils.

The structural maintenance of chromosomes (SMC) protein complexes-cohesin, condensin, and the Smc5/6 complex (Smc5/6)-are essential for chromosome function. At the molecular level, these complexes fold DNA by loop extrusion. Accordingly, cohesin creates chromosome loops in interphase, and condensin compacts mitotic chromosomes. However, the role of Smc5/6's recently discovered DNA loop extrusion activity is unknown. Here, we uncover that Smc5/6 associates with transcription-induced positively supercoiled DNA at cohesin-dependent loop boundaries on budding yeast (Saccharomyces cerevisiae) chromosomes. Mechanistically, single-molecule imaging reveals that dimers of Smc5/6 specifically recognize the tip of positively supercoiled DNA plectonemes and efficiently initiate loop extrusion to gather the supercoiled DNA into a large plectonemic loop. Finally, Hi-C analysis shows that Smc5/6 links chromosomal regions containing transcription-induced positive supercoiling in cis. Altogether, our findings indicate that Smc5/6 controls the three-dimensional organization of chromosomes by recognizing and initiating loop extrusion on positively supercoiled DNA.

Highly rigid H3.1/H3.2-H3K9me3 domains set a barrier for cell fate reprogramming in trophoblast stem cells.

The placenta is a highly evolved, specialized organ in mammals. It differs from other organs in that it functions only for fetal maintenance during gestation. Therefore, there must be intrinsic mechanisms that guarantee its unique functions. To address this question, we comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). Our genome-wide, high-throughput analyses revealed that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures with a high degree of histone H3.1/3.2-H3K9me3 accumulation, which we termed TSC-defined highly heterochromatinized domains (THDs). Importantly, depletion of THDs by knockdown of CAF1, an H3.1/3.2 chaperone, resulted in down-regulation of TSC markers, such as Cdx2 and Elf5, and up-regulation of the pluripotent marker Oct3/4, indicating that THDs maintain the trophoblastic nature of TSCs. Furthermore, our nuclear transfer technique revealed that THDs are highly resistant to genomic reprogramming. However, when H3K9me3 was removed, the TSC genome was fully reprogrammed, giving rise to the first TSC cloned offspring. Interestingly, THD-like domains are also present in mouse and human placental cells in vivo, but not in other cell types. Thus, THDs are genomic architectures uniquely developed in placental lineage cells, which serve to protect them from fate reprogramming to stably maintain placental function.

Morc1 reestablishes H3K9me3 heterochromatin on piRNA-targeted transposons in gonocytes.

To maintain fertility, male mice re-repress transposable elements (TEs) that were de-silenced in the early gonocytes before their differentiation into spermatogonia. However, the mechanism of TE silencing re-establishment remains unknown. Here, we found that the DNA-binding protein Morc1, in cooperation with the methyltransferase SetDB1, deposits the repressive histone mark H3K9me3 on a large fraction of activated TEs, leading to heterochromatin. Morc1 also triggers DNA methylation, but TEs targeted by Morc1-driven DNA methylation only slightly overlapped with those repressed by Morc1/SetDB1-dependent heterochromatin formation, suggesting that Morc1 silences TEs in two different manners. In contrast, TEs regulated by Morc1 and Miwi2, the nuclear PIWI-family protein, almost overlapped. Miwi2 binds to PIWI-interacting RNAs (piRNAs) that base-pair with TE mRNAs via sequence complementarity, while Morc1 DNA binding is not sequence specific, suggesting that Miwi2 selects its targets, and then, Morc1 acts to repress them with cofactors. A high-ordered mechanism of TE repression in gonocytes has been identified.

Clover: An unbiased method for prioritizing differentially expressed genes using a data-driven approach.

Identifying key genes from a list of differentially expressed genes (DEGs) is a critical step in transcriptome analysis. However, current methods, including Gene Ontology analysis and manual annotation, essentially rely on existing knowledge, which is highly biased depending on the extent of the literature. As a result, understudied genes, some of which may be associated with important molecular mechanisms, are often ignored or remain obscure. To address this problem, we propose Clover, a data-driven scoring method to specifically highlight understudied genes. Clover aims to prioritize genes associated with important molecular mechanisms by integrating three metrics: the likelihood of appearing in the DEG list, tissue specificity, and number of publications. We applied Clover to Alzheimer's disease data and confirmed that it successfully detected known associated genes. Moreover, Clover effectively prioritized understudied but potentially druggable genes. Overall, our method offers a novel approach to gene characterization and has the potential to expand our understanding of gene functions. Clover is an open-source software written in Python3 and available on GitHub at https://github.com/G708/Clover.

Histone H3.3 chaperone HIRA renders stress-responsive genes poised for prospective lethal stresses in acquired tolerance.

Appropriate responses to environmental challenges are imperative for the survival of all living organisms. Exposure to low-dose stresses is recognized to yield increased cellular fitness, a phenomenon termed hormesis. However, our molecular understanding of how cells respond to low-dose stress remains profoundly limited. Here we report that histone variant H3.3-specific chaperone, HIRA, is required for acquired tolerance, where low-dose heat stress exposure confers resistance to subsequent lethal heat stress. We found that human HIRA activates stress-responsive genes, including HSP70, by depositing histone H3.3 following low-dose stresses. These genes are also marked with histone H3 Lys-4 trimethylation and H3 Lys-9 acetylation, both active chromatin markers. Moreover, depletion of HIRA greatly diminished acquired tolerance, both in normal diploid fibroblasts and in HeLa cells. Collectively, our study revealed that HIRA is required for eliciting adaptive stress responses under environmental fluctuations and is a master regulator of stress tolerance.

CohesinDB: a comprehensive database for decoding cohesin-related epigenomes, 3D genomes and transcriptomes in human cells.

Cohesin is a multifunctional protein responsible for transcriptional regulation and chromatin organization. Cohesin binds to chromatin at tens of thousands of distinct sites in a conserved or tissue-specific manner, whereas the function of cohesin varies greatly depending on the epigenetic properties of specific chromatin loci. Cohesin also extensively mediates cis-regulatory modules (CRMs) and chromatin loops. Even though next-generation sequencing technologies have provided a wealth of information on different aspects of cohesin, the integration and exploration of the resultant massive cohesin datasets are not straightforward. Here, we present CohesinDB (https://cohesindb.iqb.u-tokyo.ac.jp), a comprehensive multiomics cohesin database in human cells. CohesinDB includes 2043 epigenomics, transcriptomics and 3D genomics datasets from 530 studies involving 176 cell types. By integrating these large-scale data, CohesinDB summarizes three types of 'cohesin objects': 751 590 cohesin binding sites, 957 868 cohesin-related chromatin loops and 2 229 500 cohesin-related CRMs. Each cohesin object is annotated with locus, cell type, classification, function, 3D genomics and cis-regulatory information. CohesinDB features a user-friendly interface for browsing, searching, analyzing, visualizing and downloading the desired information. CohesinDB contributes a valuable resource for all researchers studying cohesin, epigenomics, transcriptional regulation and chromatin organization.

HiC1Dmetrics: framework to extract various one-dimensional features from chromosome structure data.

Eukaryotic genomes are organized in a three-dimensional spatial structure. In this regard, the development of chromosome conformation capture methods has enabled studies of chromosome organization on a genomic scale. Hi-C, the high-throughput chromosome conformation capture method, can reveal a population-averaged, hierarchical chromatin structure. The typical Hi-C analysis uses a two-dimensional (2D) contact matrix that indicates contact frequencies between all possible genomic position pairs. Oftentimes, however, such a 2D matrix is not amenable to handling quantitative comparisons, visualizations and integrations across multiple datasets. Although several one-dimensional (1D) metrics have been proposed to depict structural information in Hi-C data, their effectiveness is still underappreciated. Here, we first review the currently available 1D metrics for individual Hi-C samples or two-sample comparisons and then discuss their validity and suitable analysis scenarios. We also propose several new 1D metrics to identify additional unique features of chromosome structures. We highlight that the 1D metrics are reproducible and robust for comparing and visualizing multiple Hi-C samples. Moreover, we show that 1D metrics can be easily combined with epigenome tracks to annotate chromatin states in greater details. We develop a new framework, called HiC1Dmetrics, to summarize all 1D metrics discussed in this study. HiC1Dmetrics is open-source (github.com/wangjk321/HiC1Dmetrics) and can be accessed from both command-line and web-based interfaces. Our tool constitutes a useful resource for the community of chromosome-organization researchers.

Transient Methionine Deprivation Triggers Histone Modification and Potentiates Differentiation of Induced Pluripotent Stem Cells.

Human induced pluripotent stem cells (iPSCs) require high levels of methionine (Met). Met deprivation results in a rapid decrease in intracellular S-adenosyl-methionine (SAM), poising human iPSCs for differentiation and leading to the apoptosis of undifferentiated cells. Met deprivation triggers rapid metabolic changes, including SAM, followed by reversible epigenetic modifications. Here, we show that short-term Met deprivation impairs the pluripotency network through epigenetic modification in a 3D suspension culture. The trimethylation of lysine 4 on histone H3 (H3K4me3) was drastically affected compared with other histone modifications. Short-term Met deprivation specifically affects the transcription start site (TSS) region of genes, such as those involved in the transforming growth factor ß pathway and cholesterol biosynthetic process, besides key pluripotent genes such as NANOG and POU5F1. The expression levels of these genes decreased, correlating with the loss of H3K4me3 marks. Upon differentiation, Met deprivation triggers the upregulation of various lineage-specific genes, including key definitive endoderm genes, such as GATA6. Upon differentiation, loss of H3K27me3 occurs in many endodermal genes, switching from a bivalent to a monovalent (H3K4me3) state. In conclusion, Met metabolism maintains the pluripotent network with histone marks, and their loss potentiates differentiation.

Codependency and mutual exclusivity for gene community detection from sparse single-cell transcriptome data.

Single-cell RNA-seq (scRNA-seq) can be used to characterize cellular heterogeneity in thousands of cells. The reconstruction of a gene network based on coexpression patterns is a fundamental task in scRNA-seq analyses, and the mutual exclusivity of gene expression can be critical for understanding such heterogeneity. Here, we propose an approach for detecting communities from a genetic network constructed on the basis of coexpression properties. The community-based comparison of multiple coexpression networks enables the identification of functionally related gene clusters that cannot be fully captured through differential gene expression-based analysis. We also developed a novel metric referred to as the exclusively expressed index (EEI) that identifies mutually exclusive gene pairs from sparse scRNA-seq data. EEI quantifies and ranks the exclusive expression levels of all gene pairs from binary expression patterns while maintaining robustness against a low sequencing depth. We applied our methods to glioblastoma scRNA-seq data and found that gene communities were partially conserved after serum stimulation despite a considerable number of differentially expressed genes. We also demonstrate that the identification of mutually exclusive gene sets with EEI can improve the sensitivity of capturing cellular heterogeneity. Our methods complement existing approaches and provide new biological insights, even for a large, sparse dataset, in the single-cell analysis field.

CRISPR/Cas9 Screening for Identification of Genes Required for the Growth of Ovarian Clear Cell Carcinoma Cells.

Epithelial ovarian cancer is classified into four major histological subtypes: serous, clear cell, endometrioid and mucinous. Ovarian clear cell carcinoma (OCCC) responds poorly to conventional chemotherapies and shows poor prognosis. Thus, there is a need to develop new drugs for the treatment of OCCC. In this study, we performed CRISPR/Cas9 screens against OCCC cell lines and identified candidate genes important for their proliferation. We found that quite different genes are required for the growth of ARID1A and PIK3CA mutant and wild-type OCCC cell lines, respectively. Furthermore, we found that the epigenetic regulator KDM2A and the translation regulator PAIP1 may play important roles in the growth of ARID1A and PIK3CA mutant, but not wild-type, OCCC cells. The results of our CRISPR/Cas9 screening may be useful in elucidating the molecular mechanism of OCCC tumorigenesis and in developing OCCC-targeted drugs.

PHOSPHATE exporter XPR1/SLC53A1 is required for the tumorigenicity of epithelial ovarian cancer.

Ovarian cancer is the fifth most common cause of cancer-related death in women. Ovarian clear cell carcinoma (OCCC) is a chemotherapy-resistant epithelial ovarian cancer with poor prognosis. As a basis for the development of therapeutic agents that could improve the prognosis of OCCC, we performed a screen for proteins critical for the tumorigenicity of OCCC using the CRISPR/Cas9 system. Here we show that knockdown of the phosphate exporter XPR1/SLC53A1 induces the growth arrest and apoptosis of OCCC cells in vitro. Moreover, we show that knockdown of XPR1/SLC53A1 inhibits the proliferation of OCCC cells xenografted into immunocompromised mice. These results suggest that XPR1/SLC53A1 plays a critical role in the tumorigenesis of OCCC cells. We speculate that XPR1/SLC53A1 might be a promising molecular target for the therapeutic treatment of OCCC.

Missense Mutations in NKAP Cause a Disorder of Transcriptional Regulation Characterized by Marfanoid Habitus and Cognitive Impairment.

NKAP is a ubiquitously expressed nucleoplasmic protein that is currently known as a transcriptional regulatory molecule via its interaction with HDAC3 and spliceosomal proteins. Here, we report a disorder of transcriptional regulation due to missense mutations in the X chromosome gene, NKAP. These mutations are clustered in the C-terminal region of NKAP where NKAP interacts with HDAC3 and post-catalytic spliceosomal complex proteins. Consistent with a role for the C-terminal region of NKAP in embryogenesis, nkap mutant zebrafish with a C-terminally truncated NKAP demonstrate severe developmental defects. The clinical features of affected individuals are highly conserved and include developmental delay, hypotonia, joint contractures, behavioral abnormalities, Marfanoid habitus, and scoliosis. In affected cases, transcriptome analysis revealed the presence of a unique transcriptome signature, which is characterized by the downregulation of long genes with higher exon numbers. These observations indicate the critical role of NKAP in transcriptional regulation and demonstrate that perturbations of the C-terminal region lead to developmental defects in both humans and zebrafish.

Methods for ChIP-seq analysis: A practical workflow and advanced applications.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a central method in epigenomic research. Genome-wide analysis of histone modifications, such as enhancer analysis and genome-wide chromatin state annotation, enables systematic analysis of how the epigenomic landscape contributes to cell identity, development, lineage specification, and disease. In this review, we first present a typical ChIP-seq analysis workflow, from quality assessment to chromatin-state annotation. We focus on practical, rather than theoretical, approaches for biological studies. Next, we outline various advanced ChIP-seq applications and introduce several state-of-the-art methods, including prediction of gene expression level and chromatin loops from epigenome data and data imputation. Finally, we discuss recently developed single-cell ChIP-seq analysis methodologies that elucidate the cellular diversity within complex tissues and cancers.

The novel lncRNA CALIC upregulates AXL to promote colon cancer metastasis.

Long non-coding RNAs (lncRNAs) are aberrantly expressed in many disease conditions, including cancer. Accumulating evidence indicates that some lncRNAs may play critical roles in cancer progression and metastasis. Here, we identify a set of lncRNAs that are upregulated in metastatic subpopulations isolated from colon cancer HCT116 cells in vivo and show that one of these lncRNAs, which we name CALIC, is required for the metastatic activity of colon cancer cells. We show that CALIC associates with the RNA-binding protein hnRNP-L and imparts specificity to hnRNP-L-mediated gene expression. Furthermore, we demonstrate that the CALIC/hnRNP-L complex upregulates the tyrosine kinase receptor AXL and that knockdown of CALIC or AXL using shRNA in colon cancer cells attenuates their ability to form metastases in mice. These results suggest that the CALIC/hnRNP-L complex enhances the metastatic potential of colon cancer cells.

Replisome stability at defective DNA replication forks is independent of S phase checkpoint kinases.

The S phase checkpoint pathway preserves genome stability by protecting defective DNA replication forks, but the underlying mechanisms are still understood poorly. Previous work with budding yeast suggested that the checkpoint kinases Mec1 and Rad53 might prevent collapse of the replisome when nucleotide concentrations are limiting, thereby allowing the subsequent resumption of DNA synthesis. Here we describe a direct analysis of replisome stability in budding yeast cells lacking checkpoint kinases, together with a high-resolution view of replisome progression across the genome. Surprisingly, we find that the replisome is stably associated with DNA replication forks following replication stress in the absence of Mec1 or Rad53. A component of the replicative DNA helicase is phosphorylated within the replisome in a Mec1-dependent manner upon replication stress, and our data indicate that checkpoint kinases control replisome function rather than stability, as part of a multifaceted response that allows cells to survive defects in chromosome replication.

MRG15-mediated tethering of PALB2 to unperturbed chromatin protects active genes from genotoxic stress.

The partner and localiser of BRCA2 (PALB2) plays important roles in the maintenance of genome integrity and protection against cancer. Although PALB2 is commonly described as a repair factor recruited to sites of DNA breaks, recent studies provide evidence that PALB2 also associates with unperturbed chromatin. Here, we investigated the previously poorly described role of chromatin-associated PALB2 in undamaged cells. We found that PALB2 associates with active genes through its major binding partner, MRG15, which recognizes histone H3 trimethylated at lysine 36 (H3K36me3) by the SETD2 methyltransferase. Missense mutations that ablate PALB2 binding to MRG15 confer elevated sensitivity to the topoisomerase inhibitor camptothecin (CPT) and increased levels of aberrant metaphase chromosomes and DNA stress in gene bodies, which were suppressed by preventing DNA replication. Remarkably, the level of PALB2 at genic regions was frequently decreased, rather than increased, upon CPT treatment. We propose that the steady-state presence of PALB2 at active genes, mediated through the SETD2/H3K36me3/MRG15 axis, ensures an immediate response to DNA stress and therefore effective protection of these regions during DNA replication. This study provides a conceptual advance in demonstrating that the constitutive chromatin association of repair factors plays a key role in the maintenance of genome stability and furthers our understanding of why PALB2 defects lead to human genome instability syndromes.

Telomere-binding protein Taz1 controls global replication timing through its localization near late replication origins in fission yeast.

In eukaryotes, the replication of chromosome DNA is coordinated by a replication timing program that temporally regulates the firing of individual replication origins. However, the molecular mechanism underlying the program remains elusive. Here, we report that the telomere-binding protein Taz1 plays a crucial role in the control of replication timing in fission yeast. A DNA element located proximal to a late origin in the chromosome arm represses initiation from the origin in early S phase. Systematic deletion and substitution experiments demonstrated that two tandem telomeric repeats are essential for this repression. The telomeric repeats recruit Taz1, a counterpart of human TRF1 and TRF2, to the locus. Genome-wide analysis revealed that Taz1 regulates about half of chromosomal late origins, including those in subtelomeres. The Taz1-mediated mechanism prevents Dbf4-dependent kinase (DDK)-dependent Sld3 loading onto the origins. Our results demonstrate that the replication timing program in fission yeast uses the internal telomeric repeats and binding of Taz1.

Assembly of Slx4 signaling complexes behind DNA replication forks.

Obstructions to replication fork progression, referred to collectively as DNA replication stress, challenge genome stability. In Saccharomyces cerevisiae, cells lacking RTT107 or SLX4 show genome instability and sensitivity to DNA replication stress and are defective in the completion of DNA replication during recovery from replication stress. We demonstrate that Slx4 is recruited to chromatin behind stressed replication forks, in a region that is spatially distinct from that occupied by the replication machinery. Slx4 complex formation is nucleated by Mec1 phosphorylation of histone H2A, which is recognized by the constitutive Slx4 binding partner Rtt107. Slx4 is essential for recruiting the Mec1 activator Dpb11 behind stressed replication forks, and Slx4 complexes are important for full activity of Mec1. We propose that Slx4 complexes promote robust checkpoint signaling by Mec1 by stably recruiting Dpb11 within a discrete domain behind the replication fork, during DNA replication stress.