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.

Hypoxia-induced immortalization of primary cells depends on Tfcp2L1 expression

Cellular senescence is a stress response mechanism that induces proliferative arrest. Hypoxia can bypass senescence and extend the lifespan of primary cells, mainly by decreasing oxidative damage. However, how hypoxia promotes these effects prior to malignant transformation is unknown. Here we observed that the lifespan of mouse embryonic fibroblasts (MEFs) is increased when they are cultured in hypoxia by reducing the expression of p16INK4a, p15INK4b and p21Cip1. We found that proliferating MEFs in hypoxia overexpress Tfcp2l1, which is a main regulator of pluripotency and self-renewal in embryonic stem cells, as well as stemness genes including Oct3/4, Sox2 and Nanog. Tfcp2l1 expression is lost during culture in normoxia, and its expression in hypoxia is regulated by Hif1α. Consistently, its overexpression in hypoxic levels increases the lifespan of MEFs and promotes the overexpression of stemness genes. ATAC-seq and Chip-seq experiments showed that Tfcp2l1 regulates genes that control proliferation and stemness such as Sox2, Sox9, Jarid2 and Ezh2. Additionally, Tfcp2l1 can replicate the hypoxic effect of increasing cellular reprogramming. Altogether, our data suggest that the activation of Tfcp2l1 by hypoxia contributes to immortalization prior to malignant transformation, facilitating tumorigenesis and dedifferentiation by regulating Sox2, Sox9, and Jarid2.

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.

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.

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.

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.

Yra1 is an essential nuclear RNA-binding protein that plays a role in mRNA export in Saccharomyces cerevisiae. The cellular levels of Yra1 are tightly auto-regulated by splicing of an unusual intron in its pre-mRNA, removal of which causes Yra1 overexpression that results in a dominant-negative growth defect and mRNA export defect. We wondered whether or not YRA1 overexpression has an effect on genome integrity that could explain the loss of cell viability. Our analyses reveal that YRA1 overexpression causes DNA damage, confers a hyperrecombination phenotype that depends on transcription and that is mediated by RNA:DNA hybrids. YRA1 overexpression also leads to a cell senescence-like phenotype and telomere shortening. We show by ChIP-chip analysis that Yra1 binds to active chromatin and Y’ telomeric regions when it is overexpressed, in agreement with a possible role of this mRNP factor in the maintenance of telomere integrity. Our data indicate that YRA1 overexpression correlates with replication impairment as inferred by the reduction of BrdU incorporation and the increase of Rrm3 recruitment, a helicase involved in replication fork progression, at transcribed genes and Y’ regions. We conclude that the stoichiometry of specific RNA-binding factors such as Yra1 at telomeres is critical for genome integrity and for preventing transcription-replication conflicts.

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.