Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development.

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

Essential role of PLD2 in hypoxia-induced stemness and therapy resistance in ovarian tumors.

BACKGROUND: Hypoxia in solid tumors is an important source of chemoresistance that can determine poor patient prognosis. Such chemoresistance relies on the presence of cancer stem cells (CSCs), and hypoxia promotes their generation through transcriptional activation by HIF transcription factors. METHODS: We used ovarian cancer (OC) cell lines, xenograft models, OC patient samples, transcriptional databases, induced pluripotent stem cells (iPSCs) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). RESULTS: Here, we show that hypoxia induces CSC formation and chemoresistance in ovarian cancer through transcriptional activation of the PLD2 gene. Mechanistically, HIF-1α activates PLD2 transcription through hypoxia response elements, and both hypoxia and PLD2 overexpression lead to increased accessibility around stemness genes, detected by ATAC-seq, at sites bound by AP-1 transcription factors. This in turn provokes a rewiring of stemness genes, including the overexpression of SOX2, SOX9 or NOTCH1. PLD2 overexpression also leads to decreased patient survival, enhanced tumor growth and CSC formation, and increased iPSCs reprograming, confirming its role in dedifferentiation to a stem-like phenotype. Importantly, hypoxia-induced stemness is dependent on PLD2 expression, demonstrating that PLD2 is a major determinant of de-differentiation of ovarian cancer cells to stem-like cells in hypoxic conditions. Finally, we demonstrate that high PLD2 expression increases chemoresistance to cisplatin and carboplatin treatments, both in vitro and in vivo, while its pharmacological inhibition restores sensitivity. CONCLUSIONS: Altogether, our work highlights the importance of the HIF-1α-PLD2 axis for CSC generation and chemoresistance in OC and proposes an alternative treatment for patients with high PLD2 expression.

MacroH2As regulate enhancer-promoter contacts affecting enhancer activity and sensitivity to inflammatory cytokines.

MacroH2A histone variants have a function in gene regulation that is poorly understood at the molecular level. We report that macroH2A1.2 and macroH2A2 modulate the transcriptional ground state of cancer cells and how they respond to inflammatory cytokines. Removal of macroH2A1.2 and macroH2A2 in hepatoblastoma cells affects the contact frequency of promoters and distal enhancers coinciding with changes in enhancer activity or preceding them in response to the cytokine tumor necrosis factor alpha. Although macroH2As regulate genes in both directions, they globally facilitate the nuclear factor κB (NF-κB)-mediated response. In contrast, macroH2As suppress the response to the pro-inflammatory cytokine interferon gamma. MacroH2A2 has a stronger contribution to gene repression than macroH2A1.2. Taken together, our results suggest that macroH2As have a role in regulating the response of cancer cells to inflammatory signals on the level of chromatin structure. This is likely relevant for the interaction of cancer cells with immune cells of their microenvironment.

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression.

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease.

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Genome-Wide Mapping of Histone H3 Lysine 4 Trimethylation (H3K4me3) and Its Involvement in Fatty Acid Biosynthesis in Sunflower Developing Seeds.

Histone modifications are of paramount importance during plant development. Investigating chromatin remodeling in developing oilseeds sheds light on the molecular mechanisms controlling fatty acid metabolism and facilitates the identification of new functional regions in oil crop genomes. The present study characterizes the epigenetic modifications H3K4me3 in relationship with the expression of fatty acid-related genes and transcription factors in developing sunflower seeds. Two master transcriptional regulators identified in this analysis, VIV1 (homologous to Arabidopsis ABI3) and FUS3, cooperate in the regulation of WRINKLED 1, a transcriptional factor regulating glycolysis, and fatty acid synthesis in developing oilseeds.

Functional Conservation of Divergent p63-Bound cis-Regulatory Elements.

The transcription factor p63 is an essential regulator of vertebrate ectoderm development, including epidermis, limbs, and craniofacial tissues. Here, we have investigated the evolutionary conservation of p63 binding sites (BSs) between zebrafish and human. First, we have analyzed sequence conservation of p63 BSs by comparing ChIP-seq data from human keratinocytes and zebrafish embryos, observing a very poor conservation. Next, we compared the gene regulatory network orchestrated by p63 in both species and found a high overlap between them, suggesting a high degree of functional conservation during evolution despite sequence divergence and the large evolutionary distance. Finally, we used transgenic reporter assays in zebrafish embryos to functionally validate a set of equivalent p63 BSs from zebrafish and human located close to genes involved in epidermal development. Reporter expression was driven by human and zebrafish BSs to many common tissues related to p63 expression domains. Therefore, we conclude that the gene regulatory network controlled by p63 is highly conserved across vertebrates despite the fact that p63-bound regulatory elements show high divergence.

Pioneer and repressive functions of p63 during zebrafish embryonic ectoderm specification.

The transcription factor p63 is a master regulator of ectoderm development. Although previous studies show that p63 triggers epidermal differentiation in vitro, the roles of p63 in developing embryos remain poorly understood. Here, we use zebrafish embryos to analyze in vivo how p63 regulates gene expression during development. We generate tp63-knock-out mutants that recapitulate human phenotypes and show down-regulated epidermal gene expression. Following p63-binding dynamics, we find two distinct functions clearly separated in space and time. During early development, p63 binds enhancers associated to neural genes, limiting Sox3 binding and reducing neural gene expression. Indeed, we show that p63 and Sox3 are co-expressed in the neural plate border. On the other hand, p63 acts as a pioneer factor by binding non-accessible chromatin at epidermal enhancers, promoting their opening and epidermal gene expression in later developmental stages. Therefore, our results suggest that p63 regulates cell fate decisions during vertebrate ectoderm specification.

The Nup84 complex coordinates the DNA damage response to warrant genome integrity.

DNA lesions interfere with cellular processes such as transcription and replication and need to be adequately resolved to warrant genome integrity. Beyond their primary role in molecule transport, nuclear pore complexes (NPCs) function in other processes such as transcription, nuclear organization and DNA double strand break (DSB) repair. Here we found that the removal of UV-induced DNA lesions by nucleotide excision repair (NER) is compromised in the absence of the Nup84 nuclear pore component. Importantly, nup84Δ cells show an exacerbated sensitivity to UV in early S phase and delayed replication fork progression, suggesting that unrepaired spontaneous DNA lesions persist during S phase. In addition, nup84Δ cells are defective in the repair of replication-born DSBs by sister chromatid recombination (SCR) and rely on post-replicative repair functions for normal proliferation, indicating dysfunctions in the cellular pathways that enable replication on damaged DNA templates. Altogether, our data reveal a central role of the NPC in the DNA damage response to facilitate replication progression through damaged DNA templates by promoting efficient NER and SCR and preventing chromosomal rearrangements.

Excess of Yra1 RNA-Binding Factor Causes Transcription-Dependent Genome Instability, Replication Impairment and Telomere Shortening.

Yra1 is an essential nuclear factor of the evolutionarily conserved family of hnRNP-like export factors that when overexpressed impairs mRNA export and cell growth. To investigate further the relevance of proper Yra1 stoichiometry in the cell, we overexpressed Yra1 by transforming yeast cells with YRA1 intron-less constructs and analyzed its effect on gene expression and genome integrity. We found that YRA1 overexpression induces DNA damage and leads to a transcription-associated hyperrecombination phenotype that is mediated by RNA:DNA hybrids. In addition, it confers a genome-wide replication retardation as seen by reduced BrdU incorporation and accumulation of the Rrm3 helicase. In addition, YRA1 overexpression causes a cell senescence-like phenotype and telomere shortening. ChIP-chip analysis shows that overexpressed Yra1 is loaded to transcribed chromatin along the genome and to Y' telomeric regions, where Rrm3 is also accumulated, suggesting an impairment of telomere replication. Our work not only demonstrates that a proper stoichiometry of the Yra1 mRNA binding and export factor is required to maintain genome integrity and telomere homeostasis, but suggests that the cellular imbalance between transcribed RNA and specific RNA-binding factors may become a major cause of genome instability mediated by co-transcriptional replication impairment.

R loops: new modulators of genome dynamics and function.

R loops are nucleic acid structures composed of an RNA-DNA hybrid and a displaced single-stranded DNA. Recently, evidence has emerged that R loops occur more often in the genome and have greater physiological relevance, including roles in transcription and chromatin structure, than was previously predicted. Importantly, however, R loops are also a major threat to genome stability. For this reason, several DNA and RNA metabolism factors prevent R-loop formation in cells. Dysfunction of these factors causes R-loop accumulation, which leads to replication stress, genome instability, chromatin alterations or gene silencing, phenomena that are frequently associated with cancer and a number of genetic diseases. We review the current knowledge of the mechanisms controlling R loops and their putative relationship with disease.

A genome-wide function of THSC/TREX-2 at active genes prevents transcription-replication collisions.

The THSC/TREX-2 complex of Saccharomyces cerevisiae mediates the anchoring of transcribed genes to the nuclear pore, linking transcription elongation with mRNA export and genome stability, as shown for specific reporters. However, it is still unknown whether the function of TREX-2 is global and the reason for its relevant role in genome integrity. Here, by studying two TREX-2 representative subunits, Thp1 and Sac3, we show that TREX-2 has a genome-wide role in gene expression. Both proteins show similar distributions along the genome, with a gradient disposition at active genes that increases towards the 3' end. Thp1 and Sac3 have a relevant impact on the expression of long, G+C-rich and highly transcribed genes. Interestingly, replication impairment detected by the genome-wide accumulation of the replicative Rrm3 helicase is increased preferentially at highly expressed genes in the thp1Δ and sac3Δ mutants analyzed. Therefore, our work provides evidence of a function of TREX-2 at the genome-wide level and suggests a role for TREX-2 in preventing transcription-replication conflicts, as a source of genome instability derived from a defective messenger ribonucleoprotein particle (mRNP) biogenesis.

Npl3, a new link between RNA-binding proteins and the maintenance of genome integrity.

The mRNA is co-transcriptionally bound by a number of RNA-binding proteins (RBPs) that contribute to its processing and formation of an export-competent messenger ribonucleoprotein particle (mRNP). In the last few years, increasing evidence suggests that RBPs play a key role in preventing transcription-associated genome instability. Part of this instability is mediated by the accumulation of co-transcriptional R loops, which may impair replication fork (RF) progression due to collisions between transcription and replication machineries. In addition, some RBPs have been implicated in DNA repair and/or the DNA damage response (DDR). Recently, the Npl3 protein, one of the most abundant heterogeneous nuclear ribonucleoproteins (hnRNPs) in yeast, has been shown to prevent transcription-associated genome instability and accumulation of RF obstacles, partially associated with R-loop formation. Interestingly, Npl3 seems to have additional functions in DNA repair, and npl3∆ mutants are highly sensitive to genotoxic agents, such as the antitumor drug trabectedin. Here we discuss the role of Npl3 in particular, and RBPs in general, in the connection of transcription with replication and genome instability, and its effect on the DDR.

The Npl3 hnRNP prevents R-loop-mediated transcription-replication conflicts and genome instability.

Transcription is a major obstacle for replication fork (RF) progression and a cause of genome instability. Part of this instability is mediated by cotranscriptional R loops, which are believed to increase by suboptimal assembly of the nascent messenger ribonucleoprotein particle (mRNP). However, no clear evidence exists that heterogeneous nuclear RNPs (hnRNPs), the basic mRNP components, prevent R-loop stabilization. Here we show that yeast Npl3, the most abundant RNA-binding hnRNP, prevents R-loop-mediated genome instability. npl3Δ cells show transcription-dependent and R-loop-dependent hyperrecombination and genome-wide replication obstacles as determined by accumulation of the Rrm3 helicase. Such obstacles preferentially occur at long and highly expressed genes, to which Npl3 is preferentially bound in wild-type cells, and are reduced by RNase H1 overexpression. The resulting replication stress confers hypersensitivity to double-strand break-inducing agents. Therefore, our work demonstrates that mRNP factors are critical for genome integrity and opens the option of using them as therapeutic targets in anti-cancer treatment.

R loops are linked to histone H3 S10 phosphorylation and chromatin condensation.

R loops are transcription byproducts that constitute a threat to genome integrity. Here we show that R loops are tightly linked to histone H3 S10 phosphorylation (H3S10P), a mark of chromatin condensation. Chromatin immunoprecipitation (ChIP)-on-chip (ChIP-chip) analyses reveal H3S10P accumulation at centromeres, pericentromeric chromatin, and a large number of active open reading frames (ORFs) in R-loop-accumulating yeast cells, better observed in G1. Histone H3S10 plays a key role in maintaining genome stability, as scored by ectopic recombination and plasmid loss, Rad52 foci, and Rad53 checkpoint activation. H3S10P coincides with the presence of DNA-RNA hybrids, is suppressed by ribonuclease H overexpression, and causes reduced accessibility of restriction endonucleases, implying a tight connection between R loops, H3S10P, and chromatin compaction. Such histone modifications were also observed in R-loop-accumulating Caenorhabditis elegans and HeLa cells. We therefore provide a role of RNA in chromatin structure essential to understand how R loops modulate genome dynamics.