The genetic landscape of origins of replication in P. falciparum.

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

TAp73 is a central transcriptional regulator of airway multiciliogenesis.

Motile multiciliated cells (MCCs) have critical roles in respiratory health and disease and are essential for cleaning inhaled pollutants and pathogens from airways. Despite their significance for human disease, the transcriptional control that governs multiciliogenesis remains poorly understood. Here we identify TP73, a p53 homolog, as governing the program for airway multiciliogenesis. Mice with TP73 deficiency suffer from chronic respiratory tract infections due to profound defects in ciliogenesis and complete loss of mucociliary clearance. Organotypic airway cultures pinpoint TAp73 as necessary and sufficient for basal body docking, axonemal extension, and motility during the differentiation of MCC progenitors. Mechanistically, cross-species genomic analyses and complete ciliary rescue of knockout MCCs identify TAp73 as the conserved central transcriptional integrator of multiciliogenesis. TAp73 directly activates the key regulators FoxJ1, Rfx2, Rfx3, and miR34bc plus nearly 50 structural and functional ciliary genes, some of which are associated with human ciliopathies. Our results position TAp73 as a novel central regulator of MCC differentiation.

Chromatin immunoprecipitation indirect peaks highlight long-range interactions of insulator proteins and Pol II pausing.

Eukaryotic chromosomes are partitioned into topologically associating domains (TADs) that are demarcated by distinct insulator-binding proteins (IBPs) in Drosophila. Whether IBPs regulate specific long-range contacts and how this may impact gene expression remains unclear. Here we identify "indirect peaks" of multiple IBPs that represent their distant sites of interactions through long-range contacts. Indirect peaks depend on protein-protein interactions among multiple IBPs and their common cofactors, including CP190, as confirmed by high-resolution analyses of long-range contacts. Mutant IBPs unable to interact with CP190 impair long-range contacts as well as the expression of hundreds of distant genes that are specifically flanked by indirect peaks. Regulation of distant genes strongly correlates with RNAPII pausing, highlighting how this key transcriptional stage may trap insulator-based long-range interactions. Our data illustrate how indirect peaks may decipher gene regulatory networks through specific long-range interactions.

Structural and functional insights into CWC27/CWC22 heterodimer linking the exon junction complex to spliceosomes.

Human CWC27 is an uncharacterized splicing factor and mutations in its gene are linked to retinal degeneration and other developmental defects. We identify the splicing factor CWC22 as the major CWC27 partner. Both CWC27 and CWC22 are present in published Bact spliceosome structures, but no interacting domains are visible. Here, the structure of a CWC27/CWC22 heterodimer bound to the exon junction complex (EJC) core component eIF4A3 is solved at 3Å-resolution. According to spliceosomal structures, the EJC is recruited in the C complex, once CWC27 has left. Our 3D structure of the eIF4A3/CWC22/CWC27 complex is compatible with the Bact spliceosome structure but not with that of the C complex, where a CWC27 loop would clash with the EJC core subunit Y14. A CWC27/CWC22 building block might thus form an intermediate landing platform for eIF4A3 onto the Bact complex prior to its conversion into C complex. Knock-down of either CWC27 or CWC22 in immortalized retinal pigment epithelial cells affects numerous common genes, indicating that these proteins cooperate, targeting the same pathways. As the most up-regulated genes encode factors involved in inflammation, our findings suggest a possible link to the retinal degeneration associated with CWC27 deficiencies.

Epigenetic alterations in longevity regulators, reduced life span, and exacerbated aging-related pathology in old father offspring mice.

Advanced age is not only a major risk factor for a range of disorders within an aging individual but may also enhance susceptibility for disease in the next generation. In humans, advanced paternal age has been associated with increased risk for a number of diseases. Experiments in rodent models have provided initial evidence that paternal age can influence behavioral traits in offspring animals, but the overall scope and extent of paternal age effects on health and disease across the life span remain underexplored. Here, we report that old father offspring mice showed a reduced life span and an exacerbated development of aging traits compared with young father offspring mice. Genome-wide epigenetic analyses of sperm from aging males and old father offspring tissue identified differentially methylated promoters, enriched for genes involved in the regulation of evolutionarily conserved longevity pathways. Gene expression analyses, biochemical experiments, and functional studies revealed evidence for an overactive mTORC1 signaling pathway in old father offspring mice. Pharmacological mTOR inhibition during the course of normal aging ameliorated many of the aging traits that were exacerbated in old father offspring mice. These findings raise the possibility that inherited alterations in longevity pathways contribute to intergenerational effects of aging in old father offspring mice.

The histone H2B monoubiquitination regulatory pathway is required for differentiation of multipotent stem cells.

Extensive changes in posttranslational histone modifications accompany the rewiring of the transcriptional program during stem cell differentiation. However, the mechanisms controlling the changes in specific chromatin modifications and their function during differentiation remain only poorly understood. We show that histone H2B monoubiquitination (H2Bub1) significantly increases during differentiation of human mesenchymal stem cells (hMSCs) and various lineage-committed precursor cells and in diverse organisms. Furthermore, the H2B ubiquitin ligase RNF40 is required for the induction of differentiation markers and transcriptional reprogramming of hMSCs. This function is dependent upon CDK9 and the WAC adaptor protein, which are required for H2B monoubiquitination. Finally, we show that RNF40 is required for the resolution of the H3K4me3/H3K27me3 bivalent poised state on lineage-specific genes during the transition from an inactive to an active chromatin conformation. Thus, these data indicate that H2Bub1 is required for maintaining multipotency of hMSCs and plays a central role in controlling stem cell differentiation.

Sodium butyrate rescues dopaminergic cells from alpha-synuclein-induced transcriptional deregulation and DNA damage.

Alpha-synuclein (aSyn) is considered a major culprit in Parkinson's disease (PD) pathophysiology. However, the precise molecular function of the protein remains elusive. Recent evidence suggests that aSyn may play a role on transcription regulation, possibly by modulating the acetylation status of histones. Our study aimed at evaluating the impact of wild-type (WT) and mutant A30P aSyn on gene expression, in a dopaminergic neuronal cell model, and decipher potential mechanisms underlying aSyn-mediated transcriptional deregulation. We performed gene expression analysis using RNA-sequencing in Lund Human Mesencephalic (LUHMES) cells expressing endogenous (control) or increased levels of WT or A30P aSyn. Compared to control cells, cells expressing both aSyn variants exhibited robust changes in the expression of several genes, including downregulation of major genes involved in DNA repair. WT aSyn, unlike A30P aSyn, promoted DNA damage and increased levels of phosphorylated p53. In dopaminergic neuronal cells, increased aSyn expression led to reduced levels of acetylated histone 3. Importantly, treatment with sodium butyrate, a histone deacetylase inhibitor (HDACi), rescued WT aSyn-induced DNA damage, possibly via upregulation of genes involved in DNA repair. Overall, our findings provide novel and compelling insight into the mechanisms associated with aSyn neurotoxicity in dopaminergic cells, which could be ameliorated with an HDACi. Future studies will be crucial to further validate these findings and to define novel possible targets for intervention in PD.

Insulators recruit histone methyltransferase dMes4 to regulate chromatin of flanking genes.

Chromosomal domains in Drosophila are marked by the insulator-binding proteins (IBPs) dCTCF/Beaf32 and cofactors that participate in regulating long-range interactions. Chromosomal borders are further enriched in specific histone modifications, yet the role of histone modifiers and nucleosome dynamics in this context remains largely unknown. Here, we show that IBP depletion impairs nucleosome dynamics specifically at the promoters and coding sequence of genes flanked by IBP binding sites. Biochemical purification identifies the H3K36 histone methyltransferase NSD/dMes-4 as a novel IBP cofactor, which specifically co-regulates the chromatin accessibility of hundreds of genes flanked by dCTCF/Beaf32. NSD/dMes-4 presets chromatin before the recruitment of transcriptional activators including DREF that triggers Set2/Hypb-dependent H3K36 trimethylation, nucleosome positioning, and RNA splicing. Our results unveil a model for how IBPs regulate nucleosome dynamics and gene expression through NSD/dMes-4, which may regulate H3K27me3 spreading. Our data uncover how IBPs dynamically regulate chromatin organization depending on distinct cofactors.

FORK-seq: Single-Molecule Profiling of DNA Replication.

Most genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200 nucleotide resolution using a nanopore current interpretation tool allowing the quantification of BrdU incorporation. Along pulse-chased replication intermediates from Saccharomyces cerevisiae, we can orient replication tracks and reproduce population-based replication directionality profiles. Additionally, we can map individual initiation and termination events. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

Genome-wide mapping of individual replication fork velocities using nanopore sequencing.

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

BEAF regulates cell-cycle genes through the controlled deposition of H3K9 methylation marks into its conserved dual-core binding sites.

Chromatin insulators/boundary elements share the ability to insulate a transgene from its chromosomal context by blocking promiscuous enhancer-promoter interactions and heterochromatin spreading. Several insulating factors target different DNA consensus sequences, defining distinct subfamilies of insulators. Whether each of these families and factors might possess unique cellular functions is of particular interest. Here, we combined chromatin immunoprecipitations and computational approaches to break down the binding signature of the Drosophila boundary element-associated factor (BEAF) subfamily. We identify a dual-core BEAF binding signature at 1,720 sites genome-wide, defined by five to six BEAF binding motifs bracketing 200 bp AT-rich nuclease-resistant spacers. Dual-cores are tightly linked to hundreds of genes highly enriched in cell-cycle and chromosome organization/segregation annotations. siRNA depletion of BEAF from cells leads to cell-cycle and chromosome segregation defects. Quantitative RT-PCR analyses in BEAF-depleted cells show that BEAF controls the expression of dual core-associated genes, including key cell-cycle and chromosome segregation regulators. beaf mutants that impair its insulating function by preventing proper interactions of BEAF complexes with the dual-cores produce similar effects in embryos. Chromatin immunoprecipitations show that BEAF regulates transcriptional activity by restricting the deposition of methylated histone H3K9 marks in dual-cores. Our results reveal a novel role for BEAF chromatin dual-cores in regulating a distinct set of genes involved in chromosome organization/segregation and the cell cycle.

FORK-seq: replication landscape of the Saccharomyces cerevisiae genome by nanopore sequencing.

Genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200-nucleotide resolution. By quantifying BrdU incorporation along pulse-chased replication intermediates from Saccharomyces cerevisiae, we orient 58,651 replication tracks reproducing population-based replication directionality profiles and map 4964 and 4485 individual initiation and termination events, respectively. Although most events cluster at known origins and fork merging zones, 9% and 18% of initiation and termination events, respectively, occur at many locations previously missed. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

RNF40 regulates gene expression in an epigenetic context-dependent manner.

BACKGROUND: Monoubiquitination of H2B (H2Bub1) is a largely enigmatic histone modification that has been linked to transcriptional elongation. Because of this association, it has been commonly assumed that H2Bub1 is an exclusively positively acting histone modification and that increased H2Bub1 occupancy correlates with increased gene expression. In contrast, depletion of the H2B ubiquitin ligases RNF20 or RNF40 alters the expression of only a subset of genes. RESULTS: Using conditional Rnf40 knockout mouse embryo fibroblasts, we show that genes occupied by low to moderate amounts of H2Bub1 are selectively regulated in response to Rnf40 deletion, whereas genes marked by high levels of H2Bub1 are mostly unaffected by Rnf40 loss. Furthermore, we find that decreased expression of RNF40-dependent genes is highly associated with widespread narrowing of H3K4me3 peaks. H2Bub1 promotes the broadening of H3K4me3 to increase transcriptional elongation, which together lead to increased tissue-specific gene transcription. Notably, genes upregulated following Rnf40 deletion, including Foxl2, are enriched for H3K27me3, which is decreased following Rnf40 deletion due to decreased expression of the Ezh2 gene. As a consequence, increased expression of some RNF40-"suppressed" genes is associated with enhancer activation via FOXL2. CONCLUSION: Together these findings reveal the complexity and context-dependency whereby one histone modification can have divergent effects on gene transcription. Furthermore, we show that these effects are dependent upon the activity of other epigenetic regulatory proteins and histone modifications.

Genome-wide chromatin and gene expression profiling during memory formation and maintenance in adult mice.

Recent evidence suggests that the formation and maintenance of memory requires epigenetic changes. In an effort to understand the spatio-temporal extent of learning and memory-related epigenetic changes we have charted genome-wide histone and DNA methylation profiles, in two different brain regions, two cell types, and three time-points, before and after learning. In this data descriptor we provide detailed information on data generation, give insights into the rationale of experiments, highlight necessary steps to assess data quality, offer guidelines for future use of the data and supply ready-to-use code to replicate the analysis results. The data provides a blueprint of the gene regulatory network underlying short- and long-term memory formation and maintenance. This 'healthy' gene regulatory network of learning can now be compared to changes in neurological or psychiatric diseases, providing mechanistic insights into brain disorders and highlighting potential therapeutic avenues.

DNA methylation changes in plasticity genes accompany the formation and maintenance of memory.

The ability to form memories is a prerequisite for an organism's behavioral adaptation to environmental changes. At the molecular level, the acquisition and maintenance of memory requires changes in chromatin modifications. In an effort to unravel the epigenetic network underlying both short- and long-term memory, we examined chromatin modification changes in two distinct mouse brain regions, two cell types and three time points before and after contextual learning. We found that histone modifications predominantly changed during memory acquisition and correlated surprisingly little with changes in gene expression. Although long-lasting changes were almost exclusive to neurons, learning-related histone modification and DNA methylation changes also occurred in non-neuronal cell types, suggesting a functional role for non-neuronal cells in epigenetic learning. Finally, our data provide evidence for a molecular framework of memory acquisition and maintenance, wherein DNA methylation could alter the expression and splicing of genes involved in functional plasticity and synaptic wiring.

Replication-fork stalling and processing at a single psoralen interstrand crosslink in Xenopus egg extracts.

Interstrand crosslink (ICL)-inducing agents block the separation of the two DNA strands. They prevent transcription and replication and are used in clinics for the treatment of cancer and skin diseases. Here, we have introduced a single psoralen ICL at a specific site in plasmid DNA using a triplex-forming-oligonucleotide (TFO)-psoralen conjugate and studied its repair in Xenopus egg extracts that support nuclear assembly and replication of plasmid DNA. Replication forks arriving from either side stalled at the psoralen ICL. In contrast to previous observations with other ICL-inducing agents, the leading strands advanced up to the lesion without any prior pausing. Subsequently, incisions were introduced on one parental strand on both sides of the ICL. These incisions could be detected whether one or both forks reached the ICL. Using small molecule inhibitors, we found that the ATR-Chk1 pathway, but not the ATM-Chk2 pathway, stimulated both the incision step and the subsequent processing of the broken replication intermediates. Our results highlight both similarities and differences in fork stalling and repair induced by psoralen and by other ICL-forming agents.

Estrogen-dependent gene transcription in human breast cancer cells relies upon proteasome-dependent monoubiquitination of histone H2B.

The estrogen receptor-α (ERα) determines the phenotype of breast cancers where it serves as a positive prognostic indicator. ERα is a well-established target for breast cancer therapy, but strategies to target its function remain of interest to address therapeutic resistance and further improve treatment. Recent findings indicate that proteasome inhibition can regulate estrogen-induced transcription, but how ERα function might be regulated was uncertain. In this study, we investigated the transcriptome-wide effects of the proteasome inhibitor bortezomib on estrogen-regulated transcription in MCF7 human breast cancer cells and showed that bortezomib caused a specific global decrease in estrogen-induced gene expression. This effect was specific because gene expression induced by the glucocorticoid receptor was unaffected by bortezomib. Surprisingly, we observed no changes in ERα recruitment or assembly of its transcriptional activation complex on ERα target genes. Instead, we found that proteasome inhibition caused a global decrease in histone H2B monoubiquitination (H2Bub1), leading to transcriptional elongation defects on estrogen target genes and to decreased chromatin dynamics overall. In confirming the functional significance of this link, we showed that RNA interference-mediated knockdown of the H2B ubiquitin ligase RNF40 decreased ERα-induced gene transcription. Surprisingly, RNF40 knockdown also supported estrogen-independent cell proliferation and activation of cell survival signaling pathways. Most importantly, we found that H2Bub1 levels decrease during tumor progression. H2Bub1 was abundant in normal mammary epithelium and benign breast tumors but absent in most malignant and metastatic breast cancers. Taken together, our findings show how ERα activity is blunted by bortezomib treatment as a result of reducing the downstream ubiquitin-dependent function of H2Bub1. In supporting a tumor suppressor role for H2Bub1 in breast cancer, our findings offer a rational basis to pursue H2Bub1-based therapies for future management of breast cancer.