Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation.

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

Reconstructing DNA methylation maps of ancient populations.

Studying premortem DNA methylation from ancient DNA (aDNA) provides a proxy for ancient gene activity patterns, and hence valuable information on evolutionary changes in gene regulation. Due to statistical limitations, current methods to reconstruct aDNA methylation maps are constrained to high-coverage shotgun samples, which comprise a small minority of available ancient samples. Most samples are sequenced using in-situ hybridization capture sequencing which targets a predefined set of genomic positions. Here, we develop methods to reconstruct aDNA methylation maps of samples that were not sequenced using high-coverage shotgun sequencing, by way of pooling together individuals to obtain a DNA methylation map that is characteristic of a population. We show that the resulting DNA methylation maps capture meaningful biological information and allow for the detection of differential methylation across populations. We offer guidelines on how to carry out comparative studies involving ancient populations, and how to control the rate of falsely discovered differentially methylated regions. The ability to reconstruct DNA methylation maps of past populations allows for the development of a whole new frontier in paleoepigenetic research, tracing DNA methylation changes throughout human history, using data from thousands of ancient samples.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness.

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality. Senescent human beta-like cells in culture undergo chromatin reorganization that leads to activation of enhancers regulating functional maturation genes and acquisition of glucose-stimulated insulin secretion capacity. Strikingly, Interferon-stimulated genes are elevated in senescent human beta cells, but genes encoding senescence-associated secretory phenotype (SASP) cytokines are not. Senescent beta cells in culture and in human tissue show elevated levels of cytoplasmic DNA, contributing to their increased interferon responsiveness. Human beta-cell senescence thus involves chromatin-driven upregulation of a functional-maturation program, and increased responsiveness of interferon-stimulated genes, changes that could increase both insulin secretion and immune reactivity.

PRC2-independent actions of H3.3K27M in embryonic stem cell differentiation.

The histone H3 variant, H3.3, is localized at specific regions in the genome, especially promoters and active enhancers, and has been shown to play important roles in development. A lysine to methionine substitution in position 27 (H3.3K27M) is a main cause of Diffuse Intrinsic Pontine Glioma (specifically Diffuse Midline Glioma, K27M-mutant), a lethal type of pediatric cancer. H3.3K27M has a dominant-negative effect by inhibiting the Polycomb Repressor Complex 2 (PRC2) activity. Here, we studied the immediate, genome-wide, consequences of the H3.3K27M mutation independent of PRC2 activity. We developed Doxycycline (Dox)-inducible mouse embryonic stem cells (ESCs) carrying a single extra copy of WT-H3.3, H3.3K27M and H3.3K27L, all fused to HA. We performed RNA-Seq and ChIP-Seq at different times following Dox induction in undifferentiated and differentiated ESCs. We find increased binding of H3.3 around transcription start sites in cells expressing both H3.3K27M and H3.3K27L compared with WT, but not in cells treated with PRC2 inhibitors. Differentiated cells carrying either H3.3K27M or H3.3K27L retain expression of ESC-active genes, in expense of expression of genes related to neuronal differentiation. Taken together, our data suggest that a modifiable H3.3K27 is required for proper histone incorporation and cellular maturation, independent of PRC2 activity.

Impaired activation of transposable elements in SARS-CoV-2 infection.

Emerging evidence shows that transposable elements (TEs) are induced in response to viral infections. This TE induction is suggested to trigger a robust and durable interferon response, providing a host defense mechanism. Here, we analyze TE expression changes in response to SARS-CoV-2 infection in different human cellular models. Unlike other viruses, SARS-CoV-2 infection does not lead to global upregulation of TEs in primary cells. We report a correlation between TEs activation and induction of interferon-related genes, suggesting that failure to activate TEs may account for the weak interferon response. Moreover, we identify two variables that explain most of the observed diverseness in immune responses: basal expression levels of TEs in the pre-infected cells and the viral load. Finally, analyzing the SARS-CoV-2 interactome and the epigenetic landscape around the TEs activated following infection, we identify SARS-CoV-2 interacting proteins, which may regulate chromatin structure and TE transcription. This work provides a possible functional explanation for SARS-CoV-2 success in its fight against the host immune system and suggests that TEs could serve as potential drug targets for COVID-19.

Open chromatin in pluripotency and reprogramming.

Pluripotent stem cells can be derived from embryos or induced from adult cells by reprogramming. They are unique among stem cells in that they can give rise to all cell types of the body. Recent findings indicate that a particularly 'open' chromatin state contributes to maintenance of pluripotency. Two principles are emerging: specific factors maintain a globally open chromatin state that is accessible for transcriptional activation; and other chromatin regulators contribute locally to the silencing of lineage-specific genes until differentiation is triggered. These same principles may apply during reacquisition of an open chromatin state upon reprogramming to pluripotency, and during de-differentiation in cancer.

Histone H1 eviction by the histone chaperone SET reduces cell survival following DNA damage.

Many chromatin remodeling and modifying proteins are involved in the DNA damage response, where they stimulate repair or induce DNA damage signaling. Interestingly, we identified that downregulation of the histone H1 (H1)-interacting protein SET results in increased resistance to a wide variety of DNA damaging agents. We found that this increased resistance does not result from alleviation of an inhibitory effect of SET on DNA repair but, rather, is the consequence of a suppressed apoptotic response to DNA damage. Furthermore, we provide evidence that the histone chaperone SET is responsible for the eviction of H1 from chromatin. Knockdown of H1 in SET-depleted cells resulted in re-sensitization of cells to DNA damage, suggesting that the increased DNA damage resistance in SET-depleted cells is the result of enhanced retention of H1 on chromatin. Finally, clonogenic survival assays showed that SET and p53 act epistatically in the attenuation of DNA damage-induced cell death. Taken together, our data indicate a role for SET in the DNA damage response as a regulator of cell survival following genotoxic stress.This article has an associated First Person interview with the first author of the paper.

AUTS2 isoforms control neuronal differentiation.

Mutations in AUTS2 are associated with autism, intellectual disability, and microcephaly. AUTS2 is expressed in the brain and interacts with polycomb proteins, yet it is still unclear how mutations in AUTS2 lead to neurodevelopmental phenotypes. Here we report that when neuronal differentiation is initiated, there is a shift in expression from a long isoform to a short AUTS2 isoform. Yeast two-hybrid screen identified the splicing factor SF3B1 as an interactor of both isoforms, whereas the polycomb group proteins, PCGF3 and PCGF5, were found to interact exclusively with the long AUTS2 isoform. Reporter assays showed that the first exons of the long AUTS2 isoform function as a transcription repressor, but the part that consist of the short isoform acts as a transcriptional activator, both influenced by the cellular context. The expression levels of PCGF3 influenced the ability of the long AUTS2 isoform to activate or repress transcription. Mouse embryonic stem cells (mESCs) with heterozygote mutations in Auts2 had an increase in cell death during in vitro corticogenesis, which was significantly rescued by overexpressing the human AUTS2 transcripts. mESCs with a truncated AUTS2 protein (missing exons 12-20) showed premature neuronal differentiation, whereas cells overexpressing AUTS2, especially the long transcript, showed increase in expression of pluripotency markers and delayed differentiation. Taken together, our data suggest that the precise expression of AUTS2 isoforms is essential for regulating transcription and the timing of neuronal differentiation.

Higher chromatin mobility supports totipotency and precedes pluripotency in vivo.

The fusion of the gametes upon fertilization results in the formation of a totipotent cell. Embryonic chromatin is expected to be able to support a large degree of plasticity. However, whether this plasticity relies on a particular conformation of the embryonic chromatin is unknown. Moreover, whether chromatin plasticity is functionally linked to cellular potency has not been addressed. Here, we adapted fluorescence recovery after photobleaching (FRAP) in the developing mouse embryo and show that mobility of the core histones H2A, H3.1, and H3.2 is unusually high in two-cell stage embryos and decreases as development proceeds. The transition toward pluripotency is accompanied by a decrease in histone mobility, and, upon lineage allocation, pluripotent cells retain higher mobility than the differentiated trophectoderm. Importantly, totipotent two-cell-like embryonic stem cells also display high core histone mobility, implying that reprogramming toward totipotency entails changes in chromatin mobility. Our data suggest that changes in chromatin dynamics underlie the transitions in cellular plasticity and that higher chromatin mobility is at the nuclear foundations of totipotency.

Transcription Factor Binding in Embryonic Stem Cells Is Constrained by DNA Sequence Repeat Symmetry.

Transcription factor (TF) recognition is dictated by the underlying DNA motif sequence specific for each TF. Here, we reveal that DNA sequence repeat symmetry plays a central role in defining TF-DNA-binding preferences. In particular, we find that different TFs bind similar symmetry patterns in the context of different developmental layers. Most TFs possess dominant preferences for similar DNA repeat symmetry types. However, in some cases, preferences of specific TFs are changed during differentiation, suggesting the importance of information encoded outside of known motif regions. Histone modifications also exhibit strong preferences for similar DNA repeat symmetry patterns unique to each type of modification. Next, using an in vivo reporter assay, we show that gene expression in embryonic stem cells can be positively modulated by the presence of genomic and computationally designed DNA oligonucleotides containing identified nonconsensus-repetitive sequence elements. This supports the hypothesis that certain nonconsensus-repetitive patterns possess a functional ability to regulate gene expression. We also performed a solution NMR experiment to probe the stability of double-stranded DNA via imino proton resonances for several double-stranded DNA sequences characterized by different repetitive patterns. We suggest that such local stability might play a key role in determining TF-DNA binding preferences. Overall, our findings show that despite the enormous sequence complexity of the TF-DNA binding landscape in differentiating embryonic stem cells, this landscape can be quantitatively characterized in simple terms using the notion of DNA sequence repeat symmetry.

Glioblastoma initiating cells are sensitive to histone demethylase inhibition due to epigenetic deregulation.

Tumor-initiating cells are a subpopulation of cells that have self-renewal capacity to regenerate a tumor. Here, we identify stem cell-like chromatin features in human glioblastoma initiating cells (GICs) and link them to a loss of the repressive histone H3 lysine 9 trimethylation (H3K9me3) mark. Increasing H3K9me3 levels by histone demethylase inhibition led to cell death in GICs but not in their differentiated counterparts. The induction of apoptosis was accompanied by a loss of the activating H3 lysine 9 acetylation (H3K9ac) modification and accumulation of DNA damage and downregulation of DNA damage response genes. Upon knockdown of histone demethylases, KDM4C and KDM7A both differentiation and DNA damage were induced. Thus, the H3K9me3-H3K9ac equilibrium is crucial for GIC viability and represents a chromatin feature that can be exploited to specifically target this tumor subpopulation.

A hyperdynamic H3.3 nucleosome marks promoter regions in pluripotent embryonic stem cells.

Histone variants and their chaperones are key regulators of eukaryotic transcription, and are critical for normal development. The histone variant H3.3 has been shown to play important roles in pluripotency and differentiation, and although its genome-wide patterns have been investigated, little is known about the role of its dynamic turnover in transcriptional regulation. To elucidate the role of H3.3 dynamics in embryonic stem cell (ESC) biology, we generated mouse ESC lines carrying a single copy of a doxycycline (Dox)-inducible HA-tagged version of H3.3 and monitored the rate of H3.3 incorporation by ChIP-seq at varying time points following Dox induction, before and after RA-induced differentiation. Comparing H3.3 turnover profiles in ESCs and RA-treated cells, we identified a hyperdynamic H3.3-containing nucleosome at the -1 position in promoters of genes expressed in ESCs. This dynamic nucleosome is restricted and shifted downstream into the +1 position following differentiation. We suggest that histone turnover dynamics provides an additional mechanism involved in expression regulation, and that a hyperdynamic -1 nucleosome marks promoters in ESCs. Our data provide evidence for regional regulation of H3.3 turnover in ESC promoters, and calls for testing, in high resolution, the dynamic behavior of additional histone variants and other structural chromatin proteins.

Systematic identification of gene family regulators in mouse and human embryonic stem cells.

Pluripotent self-renewing embryonic stem cells (ESCs) have been the focus of a growing number of high-throughput experiments, revealing the genome-wide locations of hundreds of transcription factors and histone modifications. While most of these datasets were used in a specific context, all datasets combined offer a comprehensive view of chromatin characteristics and regulatory elements that govern cell states. Here, using hundreds of datasets in ESCs, we generated colocalization maps of chromatin proteins and modifications, and built a discovery pipeline for regulatory proteins of gene families. By comparing genome-wide binding data with over-expression and knockdown analysis of hundreds of genes, we discovered that the pluripotency-related factor NR5A2 separates mitochondrial from cytosolic ribosomal genes, regulating their expression. We further show that genes with a common chromatin profile are enriched for distinct Gene Ontology (GO) categories. Our approach can be generalized to reveal common regulators of any gene group; discover novel gene families, and identify common genomic elements based on shared chromatin features.

Chromatin remodeling and bivalent histone modifications in embryonic stem cells.

Pluripotent embryonic stem cells (ESCs) are characterized by distinct epigenetic features including a relative enrichment of histone modifications related to active chromatin. Among these is tri-methylation of lysine 4 on histone H3 (H3K4me3). Several thousands of the H3K4me3-enriched promoters in pluripotent cells also contain a repressive histone mark, namely H3K27me3, a situation referred to as "bivalency". While bivalent promoters are not unique to pluripotent cells, they are relatively enriched in these cell types, largely marking developmental and lineage-specific genes which are silent but poised for immediate action. The H3K4me3 and H3K27me3 modifications are catalyzed by lysine methyltransferases which are usually found within, although not entirely limited to, the Trithorax group (TrxG) and Polycomb group (PcG) protein complexes, respectively, but these do not provide selective bivalent specificity. Recent studies highlight the family of ATP-dependent chromatin remodeling proteins as regulators of bivalent domains. Here, we discuss bivalency in general, describe the machineries that catalyze bivalent chromatin domains, and portray the emerging connection between bivalency and the action of different families of chromatin remodelers, namely INO80, esBAF, and NuRD, in pluripotent cells. We posit that chromatin remodeling proteins may enable "bivalent specificity", often selectively acting on, or selectively depleted from, bivalent domains.

PAX6 regulates melanogenesis in the retinal pigmented epithelium through feed-forward regulatory interactions with MITF.

During organogenesis, PAX6 is required for establishment of various progenitor subtypes within the central nervous system, eye and pancreas. PAX6 expression is maintained in a variety of cell types within each organ, although its role in each lineage and how it acquires cell-specific activity remain elusive. Herein, we aimed to determine the roles and the hierarchical organization of the PAX6-dependent gene regulatory network during the differentiation of the retinal pigmented epithelium (RPE). Somatic mutagenesis of Pax6 in the differentiating RPE revealed that PAX6 functions in a feed-forward regulatory loop with MITF during onset of melanogenesis. PAX6 both controls the expression of an RPE isoform of Mitf and synergizes with MITF to activate expression of genes involved in pigment biogenesis. This study exemplifies how one kernel gene pivotal in organ formation accomplishes a lineage-specific role during terminal differentiation of a single lineage.

Transcriptional competence in pluripotency.

Embryonic stem (ES) cells possess a globally open, decondensed chromatin structure that, together with trans-acting factors, supports transcriptional competence of developmentally regulated genes. However, our understanding of the mechanisms that establish transcriptional competence of specific genes is limited. In this issue of Genes & Development, Xu and colleagues (pp. 2824-2838) show that tissue-specific enhancers are actively marked by an unmethylated window in ES cells and induced pluripotent stem (iPS) cells. They propose a model and present supporting evidence to demonstrate the active involvement of pioneer transcription factors in this process. This work marks an important step toward the understanding of the mechanisms that define and maintain pluripotency, and calls for the identification of the factors that participate in the establishment of transcriptional competence in pluripotent cells.

A high-throughput chemical screen with FDA approved drugs reveals that the antihypertensive drug Spironolactone impairs cancer cell survival by inhibiting homology directed repair.

DNA double-strand breaks (DSBs) are the most severe type of DNA damage. DSBs are repaired by non-homologous end-joining or homology directed repair (HDR). Identifying novel small molecules that affect HDR is of great importance both for research use and therapy. Molecules that elevate HDR may improve gene targeting whereas inhibiting molecules can be used for chemotherapy, since some of the cancers are more sensitive to repair impairment. Here, we performed a high-throughput chemical screen for FDA approved drugs, which affect HDR in cancer cells. We found that HDR frequencies are increased by retinoic acid and Idoxuridine and reduced by the antihypertensive drug Spironolactone. We further revealed that Spironolactone impairs Rad51 foci formation, sensitizes cancer cells to DNA damaging agents, to Poly (ADP-ribose) polymerase (PARP) inhibitors and cross-linking agents and inhibits tumor growth in xenografts, in mice. This study suggests Spironolactone as a new candidate for chemotherapy.

Genes related to differentiation are correlated with the gene regulatory network structure.

MOTIVATION: Many secondary messengers, receptors and transcription factors are related to cell differentiation. Their role in cell differentiation can be affected by their position in the gene regulatory network. Here, we test whether the properties of the gene regulatory network can highlight which genes and proteins are associated with cell differentiation. We use a previously developed purely theoretical algorithm built to detect nodes that can induce a state change in Boolean gene regulatory networks, and show that most genes predicted to participate in differentiation in the theoretical framework are also experimentally known to be associated with such differentiation. These results show that genes related to differentiation are associated with specific features of the genetic regulatory network. The proposed algorithm produces a better classification than simple network measures such as the nodes degree or centrality. Boolean networks were used in many previous theoretical models. Here, we show a direct application of such networks to the detection of genes and subnetworks related to differentiation. The subnetwork emerging from the genes and edges that are predicted to be associated with differentiation are the most active molecular pathways experimentally described to be involved in cell differentiation. AVAILABILITY AND IMPLEMENTATION: http://peptibase.cs.biu.ac.il/homepage/Boolean_network_conversion_code.zip.