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.

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.

Protocol for identifying differentially expressed genes using the RumBall RNA-seq analysis platform

Here, we present a protocol for the identification of differentially expressed genes through RNA sequencing analysis. Starting with FASTQ files from public datasets, this protocol leverages RumBall within a self-contained Docker system. We describe the steps for software setup, obtaining data, read mapping, sample normalization, statistical modeling, and gene ontology enrichment. We then detail procedures for interpreting results with plots and tables. RumBall internally utilizes popular tools, ensuring a comprehensive understanding of the analysis process.

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.

Modeling embryo-endometrial interface recapitulating human embryo implantation

The initiation of human pregnancy is marked by the implantation of an embryo into the uterine environment; however, the underlying mechanisms remain largely elusive. To address this knowledge gap, we developed hormone-responsive endometrial organoids (EMO), termed apical-out (AO)-EMO, which emulate the in vivo architecture of endometrial tissue. The AO-EMO comprise an exposed apical epithelium surface, dense stromal cells, and a self-formed endothelial network. When cocultured with human embryonic stem cell-derived blastoids, the three-dimensional feto-maternal assembloid system recapitulates critical implantation stages, including apposition, adhesion, and invasion. Endometrial epithelial cells were subsequently disrupted by syncytial cells, which invade and fuse with endometrial stromal cells. We validated this fusion of syncytiotrophoblasts and stromal cells using human blastocysts. Our model provides a foundation for investigating embryo implantation and feto-maternal interactions, offering valuable insights for advancing reproductive medicine.

Churros: a Docker-based pipeline for large-scale epigenomic analysis

The epigenome, which reflects the modifications on chromatin or DNA sequences, provides crucial insight into gene expression regulation and cellular activity. With the continuous accumulation of epigenomic datasets such as chromatin immunoprecipitation followed by sequencing (ChIP-seq) data, there is a great demand for a streamlined pipeline to consistently process them, especially for large-dataset comparisons involving hundreds of samples. Here, we present Churros, an end-to-end epigenomic analysis pipeline that is environmentally independent and optimized for handling large-scale data. We successfully demonstrated the effectiveness of Churros by analyzing large-scale ChIP-seq datasets with the hg38 or Telomere-to-Telomere (T2T) human reference genome. We found that applying T2T to the typical analysis workflow has important impacts on read mapping, quality checks, and peak calling. We also introduced a useful feature to study context-specific epigenomic landscapes. Churros will contribute a comprehensive and unified resource for analyzing large-scale epigenomic data.

Cell cycle-dependent gene networks for cell proliferation activated by nuclear CK2α complexes

Nuclear expression of protein kinase CK2α is reportedly elevated in human carcinomas, but mechanisms underlying its variable localization in cells are poorly understood. This study demonstrates a functional connection between nuclear CK2 and gene expression in relation to cell proliferation. Growth stimulation of quiescent human normal fibroblasts and phospho-proteomic analysis identified a pool of CK2α that is highly phosphorylated at serine 7. Phosphorylated CK2α translocates into the nucleus, and this phosphorylation appears essential for nuclear localization and catalytic activity. Protein signatures associated with nuclear CK2 complexes reveal enrichment of apparently unique transcription factors and chromatin remodelers during progression through the G1 phase of the cell cycle. Chromatin immunoprecipitation-sequencing profiling demonstrated recruitment of CK2α to active gene loci, more abundantly in late G1 phase than in early G1, notably at transcriptional start sites of core histone genes, growth stimulus-associated genes, and ribosomal RNAs. Our findings reveal that nuclear CK2α complexes may be essential to facilitate progression of the cell cycle, by activating histone genes and triggering ribosomal biogenesis, specified in association with nuclear and nucleolar transcriptional regulators.

FXYD3 functionally demarcates an ancestral breast cancer stem cell subpopulation with features of drug-tolerant persisters

The heterogeneity of cancer stem cells (CSCs) within tumors presents a challenge in therapeutic targeting. To decipher the cellular plasticity that fuels phenotypic heterogeneity, we undertook single-cell transcriptomics analysis in triple-negative breast cancer (TNBC) to identify subpopulations in CSCs. We found a subpopulation of CSCs with ancestral features that is marked by FXYD domain-containing ion transport regulator 3 (FXYD3), a component of the Na+/K+ pump. Accordingly, FXYD3+ CSCs evolve and proliferate, while displaying traits of alveolar progenitors that are normally induced during pregnancy. Clinically, FXYD3+ CSCs were persistent during neoadjuvant chemotherapy, hence linking them to drug-tolerant persisters (DTPs) and identifying them as crucial therapeutic targets. Importantly, FXYD3+ CSCs were sensitive to senolytic Na+/K+ pump inhibitors, such as cardiac glycosides. Together, our data indicate that FXYD3+ CSCs with ancestral features are drivers of plasticity and chemoresistance in TNBC. Targeting the Na+/K+ pump could be an effective strategy to eliminate CSCs with ancestral and DTP features that could improve TNBC prognosis.

Nutrient-driven dedifferentiation of enteroendocrine cells promotes adaptive intestinal growth in Drosophila

Post-developmental organ resizing improves organismal fitness under constantly changing nutrient environments. Although stem cell abundance is a fundamental determinant of adaptive resizing, our understanding of its underlying mechanisms remains primarily limited to the regulation of stem cell division. Here, we demonstrate that nutrient fluctuation induces dedifferentiation in the Drosophila adult midgut to drive adaptive intestinal growth. From lineage tracing and single-cell RNA sequencing, we identify a subpopulation of enteroendocrine (EE) cells that convert into functional intestinal stem cells (ISCs) in response to dietary glucose and amino acids by activating the JAK-STAT pathway. Genetic ablation of EE-derived ISCs severely impairs ISC expansion and midgut growth despite the retention of resident ISCs, and in silico modeling further indicates that EE dedifferentiation enables an efficient increase in the midgut cell number while maintaining epithelial cell composition. Our findings identify a physiologically induced dedifferentiation that ensures ISC expansion during adaptive organ growth in concert with nutrient conditions.

Context-dependent perturbations in chromatin folding and the transcriptome by cohesin and related factors

Cohesin regulates gene expression through context-specific chromatin folding mechanisms such as enhancer-promoter looping and topologically associating domain (TAD) formation by cooperating with factors such as cohesin loaders and the insulation factor CTCF. We developed a computational workflow to explore how three-dimensional (3D) structure and gene expression are regulated collectively or individually by cohesin and related factors. The main component is CustardPy, by which multi-omics datasets are compared systematically. To validate our methodology, we generated 3D genome, transcriptome, and epigenome data before and after depletion of cohesin and related factors and compared the effects of depletion. We observed diverse effects on the 3D genome and transcriptome, and gene expression changes were correlated with the splitting of TADs caused by cohesin loss. We also observed variations in long-range interactions across TADs, which correlated with their epigenomic states. These computational tools and datasets will be valuable for 3D genome and epigenome studies.

Phase-separated nuclear bodies of nucleoporin fusions promote condensation of MLL1/CRM1 and rearrangement of 3D genome structure

NUP98 and NUP214 form chimeric fusion proteins that assemble into phase-separated nuclear bodies containing CRM1, a nuclear export receptor. However, these nuclear bodies' function in controlling gene expression remains elusive. Here, we demonstrate that the nuclear bodies of NUP98::HOXA9 and SET::NUP214 promote the condensation of mixed lineage leukemia 1 (MLL1), a histone methyltransferase essential for the maintenance of HOX gene expression. These nuclear bodies are robustly associated with MLL1/CRM1 and co-localized on chromatin. Furthermore, whole-genome chromatin-conformation capture analysis reveals that NUP98::HOXA9 induces a drastic alteration in high-order genome structure at target regions concomitant with the generation of chromatin loops and/or rearrangement of topologically associating domains in a phase-separation-dependent manner. Collectively, these results show that the phase-separated nuclear bodies of nucleoporin fusion proteins can enhance the activation of target genes by promoting the condensation of MLL1/CRM1 and rearrangement of the 3D genome structure.

Comprehensive multiomics analyses reveal pervasive involvement of aberrant cohesin binding in transcriptional and chromosomal disorder of cancer cells

Chromatin organization, whose malfunction causes various diseases including cancer, is fundamentally controlled by cohesin. While cancer cells have been found with mutated or misexpressed cohesin genes, there is no comprehensive survey about the presence and role of abnormal cohesin binding in cancer cells. Here, we systematically identified ∼1% of cohesin-binding sites (701-2,633) as cancer-aberrant binding sites of cohesin (CASs). We integrated CASs with large-scale transcriptomics, epigenomics, 3D genomics, and clinical information. CASs represent tissue-specific epigenomic signatures enriched for cancer-dysregulated genes with functional and clinical significance. CASs exhibited alterations in chromatin compartments, loops within topologically associated domains, and cis-regulatory elements, indicating that CASs induce dysregulated genes through misguided chromatin structure. Cohesin depletion data suggested that cohesin binding at CASs actively regulates cancer-dysregulated genes. Overall, our comprehensive investigation suggests that aberrant cohesin binding is an essential epigenomic signature responsible for dysregulated chromatin structure and transcription in cancer cells.

Mesenchymal loss of p53 alters stem cell capacity and models human soft tissue sarcoma traits

Soft tissue sarcomas (STSs) are a heterogeneous group of tumors that originate from mesenchymal cells. p53 is frequently mutated in human STS. In this study, we found that the loss of p53 in mesenchymal stem cells (MSCs) mainly causes adult undifferentiated soft tissue sarcoma (USTS). MSCs lacking p53 show changes in stem cell properties, including differentiation, cell cycle progression, and metabolism. The transcriptomic changes and genetic mutations in murine p53-deficient USTS mimic those seen in human STS. Furthermore, single-cell RNA sequencing revealed that MSCs undergo transcriptomic alterations with aging-a risk factor for certain types of USTS-and that p53 signaling decreases simultaneously. Moreover, we found that human STS can be transcriptomically classified into six clusters with different prognoses, different from the current histopathological classification. This study paves the way for understanding MSC-mediated tumorigenesis and provides an efficient mouse model for sarcoma studies.

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.

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.

Cohesin-dependent chromosome loop extrusion is limited by transcription and stalled replication forks

Genome function depends on regulated chromosome folding, and loop extrusion by the protein complex cohesin is essential for this multilayered organization. The chromosomal positioning of cohesin is controlled by transcription, and the complex also localizes to stalled replication forks. However, the role of transcription and replication in chromosome looping remains unclear. Here, we show that reduction of chromosome-bound RNA polymerase weakens normal cohesin loop extrusion boundaries, allowing cohesin to form new long-range chromosome cis interactions. Stress response genes induced by transcription inhibition are also shown to act as new loop extrusion boundaries. Furthermore, cohesin loop extrusion during early S phase is jointly controlled by transcription and replication units. Together, the results reveal that replication and transcription machineries are chromosome-folding regulators that block the progression of loop-extruding cohesin, opening for new perspectives on cohesin's roles in genome function and stability.

Large-scale multi-omics analysis suggests specific roles for intragenic cohesin in transcriptional regulation

Cohesin, an essential protein complex for chromosome segregation, regulates transcription through a variety of mechanisms. It is not a trivial task to assign diverse cohesin functions. Moreover, the context-specific roles of cohesin-mediated interactions, especially on intragenic regions, have not been thoroughly investigated. Here we perform a comprehensive characterization of cohesin binding sites in several human cell types. We integrate epigenomic, transcriptomic and chromatin interaction data to explore the context-specific functions of intragenic cohesin related to gene activation. We identify a specific subset of cohesin binding sites, decreased intragenic cohesin sites (DICs), which are negatively correlated with transcriptional regulation. A subgroup of DICs is enriched with enhancer markers and RNA polymerase II, while the others are more correlated to chromatin architecture. DICs are observed in various cell types, including cells from patients with cohesinopathy. We also implement machine learning to our data and identified genomic features for isolating DICs from all cohesin sites. These results suggest a previously unidentified function of cohesin on intragenic regions for transcriptional regulation.

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.

The 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.

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

Here, Hada et al. comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). They used genome-wide, high-throughput analyses to show that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures that have a high degree of histone H3.1/3.2–H3K9me3 accumulation, termed TSC-defined highly heterochromatinized domains (THDs), and are uniquely developed in placental lineage cells that serve to protect them from fate reprogramming to stably maintain placental function.

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.

The splicing factor DHX38/PRP16 is required for ovarian clear cell carcinoma tumorigenesis, as revealed by a CRISPR-Cas9 screen

Certain cancers, such as ovarian clear cell carcinoma (OCCC), display high levels of genetic variation between patients, making it difficult to develop effective therapies. In order to identify novel genes critical to OCCC growth, we carried out a comprehensive CRISPR‐Cas9 knockout screen against cell growth using an OCCC cell line and a normal ovarian surface epithelium cell line. We identified the gene encoding DHX38/PRP16, an ATP‐dependent RNA helicase involved in splicing, as critical for the growth and tumorigenesis of OCCC. DHX38/PRP16 knockdown in OCCC cells, but not normal cells, induces apoptosis and impairs OCCC tumorigenesis in a mouse model. Our results suggest that DHX38/PRP16 may play a role in OCCC tumorigenesis and could potentially be a promising therapeutic target.

Development of human iPSC-derived quiescent hepatic stellate cell-like cells for drug discovery and in vitro disease modeling

Hepatic stellate cells (HSCs) play a central role in the progression of liver fibrosis by producing extracellular matrices. The development of drugs to suppress liver fibrosis has been hampered by the lack of human quiescent HSCs (qHSCs) and an appropriate in vitro model that faithfully recapitulates HSC activation. In the present study, we developed a culture system to generate qHSC-like cells from human-induced pluripotent stem cells (hiPSCs) that can be converted into activated HSCs in culture. To monitor the activation process, a red fluorescent protein (RFP) gene was inserted in hiPSCs downstream of the activation marker gene actin alpha 2 smooth muscle (ACTA2). Using qHSC-like cells derived from RFP reporter iPSCs, we screened a repurposing chemical library and identified therapeutic candidates that prevent liver fibrosis. Hence, hiPSC-derived qHSC-like cells will be a useful tool to study the mechanism of HSC activation and to identify therapeutic agents.

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.

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.

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.

Toward understanding the dynamic state of 3D genome

The three-dimensional (3D) genome organization and its role in biological activities have been investigated for over a decade in the field of cell biology. Recent studies using live-imaging and polymer simulation have suggested that the higher-order chromatin structures are dynamic; the stochastic fluctuations of nucleosomes and genomic loci cannot be captured by bulk-based chromosome conformation capture techniques (Hi-C). In this review, we focus on the physical nature of the 3D genome architecture. We first describe how to decode bulk Hi-C data with polymer modeling. We then introduce our recently developed PHi-C method, a computational tool for modeling the fluctuations of the 3D genome organization in the presence of stochastic thermal noise. We also present another new method that analyzes the dynamic rheology property (represented as microrheology spectra) as a measure of the flexibility and rigidity of genomic regions over time. By applying these methods to real Hi-C data, we highlighted a temporal hierarchy embedded in the 3D genome organization; chromatin interaction boundaries are more rigid than the boundary interior, while functional domains emerge as dynamic fluctuations within a particular time interval. Our methods may bridge the gap between live-cell imaging and Hi-C data and elucidate the nature of the dynamic 3D genome organization.

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.

Publisher Correction: Paternal restraint stress affects offspring metabolism via ATF-2 dependent mechanisms in Drosophila melanogaster germ cells

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

PHi-C: deciphering Hi-C data into polymer dynamics

Genomes are spatiotemporally organized within the cell nucleus. Genome-wide chromosome conformation capture (Hi-C) technologies have uncovered the 3D genome organization. Furthermore, live-cell imaging experiments have revealed that genomes are functional in 4D. Although computational modeling methods can convert 2D Hi-C data into population-averaged static 3D genome models, exploring 4D genome nature based on 2D Hi-C data remains lacking. Here, we describe a 4D simulation method, PHi-C (polymer dynamics deciphered from Hi-C data), that depicts 4D genome features from 2D Hi-C data by polymer modeling. PHi-C allows users to interpret 2D Hi-C data as physical interaction parameters within single chromosomes. The physical interaction parameters can then be used in the simulations and analyses to demonstrate dynamic characteristics of genomic loci and chromosomes as observed in live-cell imaging experiments. PHi-C is available at https://github.com/soyashinkai/PHi-C.

Comprehensive epigenome characterization reveals diverse transcriptional regulation across human vascular endothelial cells

Background

Endothelial cells (ECs) make up the innermost layer throughout the entire vasculature. Their phenotypes and physiological functions are initially regulated by developmental signals and extracellular stimuli. The underlying molecular mechanisms responsible for the diverse phenotypes of ECs from different organs are not well understood.

Results

To characterize the transcriptomic and epigenomic landscape in the vascular system, we cataloged gene expression and active histone marks in nine types of human ECs (generating 148 genome-wide datasets) and carried out a comprehensive analysis with chromatin interaction data. We developed a robust procedure for comparative epigenome analysis that circumvents variations at the level of the individual and technical noise derived from sample preparation under various conditions. Through this approach, we identified 3765 EC-specific enhancers, some of which were associated with disease-associated genetic variations. We also identified various candidate marker genes for each EC type. We found that the nine EC types can be divided into two subgroups, corresponding to those with upper-body origins and lower-body origins, based on their epigenomic landscape. Epigenomic variations were highly correlated with gene expression patterns, but also provided unique information. Most of the deferentially expressed genes and enhancers were cooperatively enriched in more than one EC type, suggesting that the distinct combinations of multiple genes play key roles in the diverse phenotypes across EC types. Notably, many homeobox genes were differentially expressed across EC types, and their expression was correlated with the relative position of each organ in the body. This reflects the developmental origins of ECs and their roles in angiogenesis, vasculogenesis and wound healing.

Conclusions

This comprehensive analysis of epigenome characterization of EC types reveals diverse transcriptional regulation across human vascular systems. These datasets provide a valuable resource for understanding the vascular system and associated diseases.

Sensitive and robust assessment of ChIP-seq read distribution using a strand-shift profile

Motivation

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) can detect read-enriched DNA loci for point-source (e.g. transcription factor binding) and broad-source factors (e.g. various histone modifications). Although numerous quality metrics for ChIP-seq data have been developed, the 'peaks' thus obtained are still difficult to assess with respect to signal-to-noise ratio (S/N) and the percentage of false positives.

Results

We developed a quality-assessment tool for ChIP-seq data, strand-shift profile (SSP), which quantifies S/N and peak reliability without peak calling. We validated SSP in-depth using  ≥ 1000 publicly available ChIP-seq datasets along with virtual data to demonstrate that SSP provides a quantifiable and sensitive score to different S/Ns for both point- and broad-source factors, which can be standardized across diverse cell types and read depths. SSP also provides an effective criterion to judge whether a specific normalization or a rejection is required for each sample, which cannot be estimated by quality metrics currently available. Finally, we show that 'hidden-duplicate reads' cause aberrantly high S/Ns, and SSP provides an additional metric to avoid them, which can also contribute to estimation of peak mode (point- or broad-source) of samples.

Availability and implementation

SSP is open source software written in C++ and can be downloaded at https://github.com/rnakato/SSP.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Combined Loss of JMJD1A and JMJD1B Reveals Critical Roles for H3K9 Demethylation in the Maintenance of Embryonic Stem Cells and Early Embryogenesis

Histone H3 lysine 9 (H3K9) methylation is unevenly distributed in mammalian chromosomes. However, the molecular mechanism controlling the uneven distribution and its biological significance remain to be elucidated. Here, we show that JMJD1A and JMJD1B preferentially target H3K9 demethylation of gene-dense regions of chromosomes, thereby establishing an H3K9 hypomethylation state in euchromatin. JMJD1A/JMJD1B-deficient embryos died soon after implantation accompanying epiblast cell death. Furthermore, combined loss of JMJD1A and JMJD1B caused perturbed expression of metabolic genes and rapid cell death in embryonic stem cells (ESCs). These results indicate that JMJD1A/JMJD1B-meditated H3K9 demethylation has critical roles for early embryogenesis and ESC maintenance. Finally, genetic rescue experiments clarified that H3K9 overmethylation by G9A was the cause of the cell death and perturbed gene expression of JMJD1A/JMJD1B-depleted ESCs. We summarized that JMJD1A and JMJD1B, in combination, ensure early embryogenesis and ESC viability by establishing the correct H3K9 methylated epigenome.

Statistical Analysis and Quality Assessment of ChIP-seq Data with DROMPA

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis can detect protein/DNA-binding and histone-modification sites across an entire genome. As there are various factors during sample preparation that affect the obtained results, multilateral quality assessments are essential. Here, we describe a step-by-step protocol using DROMPA, a program for user-friendly ChIP-seq pipelining. DROMPA can be used for quality assessment, data normalization, visualization, peak calling, and multiple statistical analyses.

Recent advances in ChIP-seq analysis: from quality management to whole-genome annotation

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis can detect protein/DNA-binding and histone-modification sites across an entire genome. Recent advances in sequencing technologies and analyses enable us to compare hundreds of samples simultaneously; such large-scale analysis has potential to reveal the high-dimensional interrelationship level for regulatory elements and annotate novel functional genomic regions de novo. Because many experimental considerations are relevant to the choice of a method in a ChIP-seq analysis, the overall design and quality management of the experiment are of critical importance. This review offers guiding principles of computation and sample preparation for ChIP-seq analyses, highlighting the validity and limitations of the state-of-the-art procedures at each step. We also discuss the latest challenges of single-cell analysis that will encourage a new era in this field.

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.

Genome-Wide Target Analyses of Otx2 Homeoprotein in Postnatal Cortex

Juvenile brain has a unique time window, or critical period, in which neuronal circuits are remodeled by experience. Mounting evidence indicates the importance of neuronal circuit rewiring in various neurodevelopmental disorders of human cognition. We previously showed that Otx2 homeoprotein, essential for brain formation, is recaptured during postnatal maturation of parvalbumin-positive interneurons (PV cells) to activate the critical period in mouse visual cortex. Cortical Otx2 is the only interneuron-enriched transcription factor known to regulate the critical period, but its downstream targets remain unknown. Here, we used ChIP-seq (chromatin immunoprecipitation sequencing) to identify genome-wide binding sites of Otx2 in juvenile mouse cortex, and interneuron-specific RNA-seq to explore the Otx2-dependent transcriptome. Otx2-bound genes were associated with human diseases such as schizophrenia as well as critical periods. Of these genes, expression of neuronal factors involved in transcription, signal transduction and mitochondrial function was moderately and broadly affected in Otx2-deficient interneurons. In contrast to reported binding sites in the embryo, genes encoding potassium ion transporters such as KV3.1 had juvenile cortex-specific binding sites, suggesting that Otx2 is involved in regulating fast-spiking properties during PV cell maturation. Moreover, transcripts of oxidative resistance-1 (Oxr1), whose promoter has Otx2 binding sites, were markedly downregulated in Otx2-deficient interneurons. Therefore, an important role of Otx2 may be to protect the cells from the increased oxidative stress in fast-spiking PV cells. Our results suggest that coordinated expression of Otx2 targets promotes PV cell maturation and maintains its function in neuronal plasticity and disease.

Decreased cohesin in the brain leads to defective synapse development and anxiety-related behavior

Abnormal epigenetic regulation can cause the nervous system to develop abnormally. Here, we sought to understand the mechanism by which this occurs by investigating the protein complex cohesin, which is considered to regulate gene expression and, when defective, is associated with higher-level brain dysfunction and the developmental disorder Cornelia de Lange syndrome (CdLS). We generated conditional Smc3-knockout mice and observed greater dendritic complexity and larger numbers of immature synapses in the cerebral cortex of Smc3+/- mice. Smc3+/- mice also exhibited more anxiety-related behavior, which is a symptom of CdLS. Further, a gene ontology analysis after RNA-sequencing suggested the enrichment of immune processes, particularly the response to interferons, in the Smc3+/- mice. Indeed, fewer synapses formed in their cortical neurons, and this phenotype was rescued by STAT1 knockdown. Thus, low levels of cohesin expression in the developing brain lead to changes in gene expression that in turn lead to a specific and abnormal neuronal and behavioral phenotype.

Replication-Coupled PCNA Unloading by the Elg1 Complex Occurs Genome-wide and Requires Okazaki Fragment Ligation

The sliding clamp PCNA is a crucial component of the DNA replication machinery. Timely PCNA loading and unloading are central for genome integrity and must be strictly coordinated with other DNA processing steps during replication. Here, we show that the S. cerevisiae Elg1 replication factor C-like complex (Elg1-RLC) unloads PCNA genome-wide following Okazaki fragment ligation. In the absence of Elg1, PCNA is retained on chromosomes in the wake of replication forks, rather than at specific sites. Degradation of the Okazaki fragment ligase Cdc9 leads to PCNA accumulation on chromatin, similar to the accumulation caused by lack of Elg1. We demonstrate that Okazaki fragment ligation is the critical prerequisite for PCNA unloading, since Chlorella virus DNA ligase can substitute for Cdc9 in yeast and simultaneously promotes PCNA unloading. Our results suggest that Elg1-RLC acts as a general PCNA unloader and is dependent upon DNA ligation during chromosome replication.

The Deubiquitinating Enzyme USP7 Regulates Androgen Receptor Activity by Modulating Its Binding to Chromatin

The androgen receptor (AR), a nuclear receptor superfamily transcription factor, plays a key role in prostate cancer. AR signaling is the principal target for prostate cancer treatment, but current androgen-deprivation therapies cannot completely abolish AR signaling because of the heterogeneity of prostate cancers. Therefore, unraveling the mechanism of AR reactivation in androgen-depleted conditions can identify effective prostate cancer therapeutic targets. Increasing evidence indicates that AR activity is mediated by the interplay of modifying/demodifying enzymatic co-regulators. To better understand the mechanism of AR transcriptional activity regulation, we used antibodies against AR for affinity purification and identified the deubiquitinating enzyme ubiquitin-specific protease 7, USP7 as a novel AR co-regulator in prostate cancer cells. We showed that USP7 associates with AR in an androgen-dependent manner and mediates AR deubiquitination. Sequential ChIP assays indicated that USP7 forms a complex with AR on androgen-responsive elements of target genes upon stimulation with the androgen 5α-dihydrotestosterone. Further investigation indicated that USP7 is necessary to facilitate androgen-activated AR binding to chromatin. Transcriptome profile analysis of USP7-knockdown LNCaP cells also revealed the essential role of USP7 in the expression of a subset of androgen-responsive genes. Hence, inhibition of USP7 represents a compelling therapeutic strategy for the treatment of prostate cancer.

Esco1 Acetylates Cohesin via a Mechanism Different from That of Esco2

Sister chromatid cohesion is mediated by cohesin and is essential for accurate chromosome segregation. The cohesin subunits SMC1, SMC3, and Rad21 form a tripartite ring within which sister chromatids are thought to be entrapped. This event requires the acetylation of SMC3 and the association of sororin with cohesin by the acetyltransferases Esco1 and Esco2 in humans, but the functional mechanisms of these acetyltransferases remain elusive. Here, we showed that Esco1 requires Pds5, a cohesin regulatory subunit bound to Rad21, to form cohesion via SMC3 acetylation and the stabilization of the chromatin association of sororin, whereas Esco2 function was not affected by Pds5 depletion. Consistent with the functional link between Esco1 and Pds5, Pds5 interacted exclusively with Esco1, and this interaction was dependent on a unique and conserved Esco1 domain. Crucially, this interaction was essential for SMC3 acetylation and sister chromatid cohesion. Esco1 localized to cohesin localization sites on chromosomes throughout interphase in a manner that required the Esco1-Pds5 interaction, and it could acetylate SMC3 before and after DNA replication. These results indicate that Esco1 acetylates SMC3 via a mechanism different from that of Esco2. We propose that, by interacting with a unique domain of Esco1, Pds5 recruits Esco1 to chromatin-bound cohesin complexes to form cohesion. Furthermore, Esco1 acetylates SMC3 independently of DNA replication.

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.

Germline gain-of-function mutations in AFF4 cause a developmental syndrome functionally linking the super elongation complex and cohesin

Transcriptional elongation is critical for gene expression regulation during embryogenesis. The super elongation complex (SEC) governs this process by mobilizing paused RNA polymerase II (RNAP2). Using exome sequencing, we discovered missense mutations in AFF4, a core component of the SEC, in three unrelated probands with a new syndrome that phenotypically overlaps Cornelia de Lange syndrome (CdLS) that we have named CHOPS syndrome (C for cognitive impairment and coarse facies, H for heart defects, O for obesity, P for pulmonary involvement and S for short stature and skeletal dysplasia). Transcriptome and chromatin immunoprecipitation sequencing (ChIP-seq) analyses demonstrated similar alterations of genome-wide binding of AFF4, cohesin and RNAP2 in CdLS and CHOPS syndrome. Direct molecular interaction of the SEC, cohesin and RNAP2 was demonstrated. These data support a common molecular pathogenesis for CHOPS syndrome and CdLS caused by disturbance of transcriptional elongation due to alterations in genome-wide binding of AFF4 and cohesin.

Histone variants enriched in oocytes enhance reprogramming to induced pluripotent stem cells

Expression of Oct3/4, Sox2, Klf4, and c-Myc (OSKM) can reprogram somatic cells into induced pluripotent stem cells (iPSCs). Somatic cell nuclear transfer (SCNT) can also be used for reprogramming, suggesting that factors present in oocytes could potentially augment OSKM-mediated induction of pluripotency. Here, we report that two histone variants, TH2A and TH2B, which are highly expressed in oocytes and contribute to activation of the paternal genome after fertilization, enhance OSKM-dependent generation of iPSCs and can induce reprogramming with Klf4 and Oct3/4 alone. TH2A and TH2B are enriched on the X chromosome during the reprogramming process, and their expression in somatic cells increases the DNase I sensitivity of chromatin. In addition, Xist deficiency, which was reported to enhance SCNT reprogramming efficiency, stimulates iPSC generation using TH2A/TH2B in conjunction with OSKM, but not OSKM alone. Thus, TH2A/TH2B may enhance reprogramming by introducing processes that normally operate in zygotes and during SCNT.

The chromosomal association of the Smc5/6 complex depends on cohesion and predicts the level of sister chromatid entanglement

The cohesin complex, which is essential for sister chromatid cohesion and chromosome segregation, also inhibits resolution of sister chromatid intertwinings (SCIs) by the topoisomerase Top2. The cohesin-related Smc5/6 complex (Smc5/6) instead accumulates on chromosomes after Top2 inactivation, known to lead to a buildup of unresolved SCIs. This suggests that cohesin can influence the chromosomal association of Smc5/6 via its role in SCI protection. Using high-resolution ChIP-sequencing, we show that the localization of budding yeast Smc5/6 to duplicated chromosomes indeed depends on sister chromatid cohesion in wild-type and top2-4 cells. Smc5/6 is found to be enriched at cohesin binding sites in the centromere-proximal regions in both cell types, but also along chromosome arms when replication has occurred under Top2-inhibiting conditions. Reactivation of Top2 after replication causes Smc5/6 to dissociate from chromosome arms, supporting the assumption that Smc5/6 associates with a Top2 substrate. It is also demonstrated that the amount of Smc5/6 on chromosomes positively correlates with the level of missegregation in top2-4, and that Smc5/6 promotes segregation of short chromosomes in the mutant. Altogether, this shows that the chromosomal localization of Smc5/6 predicts the presence of the chromatid segregation-inhibiting entities which accumulate in top2-4 mutated cells. These are most likely SCIs, and our results thus indicate that, at least when Top2 is inhibited, Smc5/6 facilitates their resolution.

Dimeric combinations of MafB, cFos and cJun control the apoptosis-survival balance in limb morphogenesis

Apoptosis is an important mechanism for sculpting morphology. However, the molecular cascades that control apoptosis in developing limb buds remain largely unclear. Here, we show that MafB was specifically expressed in apoptotic regions of chick limb buds, and MafB/cFos heterodimers repressed apoptosis, whereas MafB/cJun heterodimers promoted apoptosis for sculpting the shape of the limbs. Chromatin immunoprecipitation sequencing in chick limb buds identified potential target genes and regulatory elements controlled by Maf and Jun. Functional analyses revealed that expression of p63 and p73, key components known to arrest the cell cycle, was directly activated by MafB and cJun. Our data suggest that dimeric combinations of MafB, cFos and cJun in developing chick limb buds control the number of apoptotic cells, and that MafB/cJun heterodimers lead to apoptosis via activation of p63 and p73.

BRCA2 coordinates the activities of cell-cycle kinases to promote genome stability

Numerous human genome instability syndromes, including cancer, are closely associated with events arising from malfunction of the essential recombinase Rad51. However, little is known about how Rad51 is dynamically regulated in human cells. Here, we show that the breast cancer susceptibility protein BRCA2, a key Rad51 binding partner, coordinates the activity of the central cell-cycle drivers CDKs and Plk1 to promote Rad51-mediated genome stability control. The soluble nuclear fraction of BRCA2 binds Plk1 directly in a cell-cycle- and CDK-dependent manner and acts as a molecular platform to facilitate Plk1-mediated Rad51 phosphorylation. This phosphorylation is important for enhancing the association of Rad51 with stressed replication forks, which in turn protects the genomic integrity of proliferating human cells. This study reveals an elaborate but highly organized molecular interplay between Rad51 regulators and has significant implications for understanding tumorigenesis and therapeutic resistance in patients with BRCA2 deficiency.

A cohesin-independent role for NIPBL at promoters provides insights in CdLS

The cohesin complex is crucial for chromosome segregation during mitosis and has recently also been implicated in transcriptional regulation and chromatin architecture. The NIPBL protein is required for the loading of cohesin onto chromatin, but how and where cohesin is loaded in vertebrate cells is unclear. Heterozygous mutations of NIPBL were found in 50% of the cases of Cornelia de Lange Syndrome (CdLS), a human developmental syndrome with a complex phenotype. However, no defects in the mitotic function of cohesin have been observed so far and the links between NIPBL mutations and the observed developmental defects are unclear. We show that NIPBL binds to chromatin in somatic cells with a different timing than cohesin. Further, we observe that high-affinity NIPBL binding sites localize to different regions than cohesin and almost exclusively to the promoters of active genes. NIPBL or cohesin knockdown reduce transcription of these genes differently, suggesting a cohesin-independent role of NIPBL for transcription. Motif analysis and comparison to published data show that NIPBL co-localizes with a specific set of other transcription factors. In cells derived from CdLS patients NIPBL binding levels are reduced and several of the NIPBL-bound genes have previously been observed to be mis-expressed in CdLS. In summary, our observations indicate that NIPBL mutations might cause developmental defects in different ways. First, defects of NIPBL might lead to cohesin-loading defects and thereby alter gene expression and second, NIPBL deficiency might affect genes directly via its role at the respective promoters.

The dynamics of genome replication using deep sequencing

Eukaryotic genomes are replicated from multiple DNA replication origins. We present complementary deep sequencing approaches to measure origin location and activity in Saccharomyces cerevisiae. Measuring the increase in DNA copy number during a synchronous S-phase allowed the precise determination of genome replication. To map origin locations, replication forks were stalled close to their initiation sites; therefore, copy number enrichment was limited to origins. Replication timing profiles were generated from asynchronous cultures using fluorescence-activated cell sorting. Applying this technique we show that the replication profiles of haploid and diploid cells are indistinguishable, indicating that both cell types use the same cohort of origins with the same activities. Finally, increasing sequencing depth allowed the direct measure of replication dynamics from an exponentially growing culture. This is the first time this approach, called marker frequency analysis, has been successfully applied to a eukaryote. These data provide a high-resolution resource and methodological framework for studying genome biology.

Chromatin immunoprecipitation protocol for mammalian cells

The genomic approach (ChIP-seq) we introduce here is now a widely used powerful tool to explore protein-DNA interaction at genome-wide level in high resolution. This technology opens up the way to understand how local event mediated by protein-protein or protein-DNA interactions lead to the dynamic changes of overall chromosome structure and how variety of proteins make a regulatory network for the faithful execution of various chromosomal functions (i.e., transcription, replication, recombination, repair, and partition).