Inherited DNA methylation primes the establishment of accessible chromatin during genome activation

Binding of PtoRAP2.12 to demethylated and accessible chromatin regions in the PtoGntK promoter stimulates growth of poplar

DNA methylation is an essential epigenetic modification for gene regulation in plant growth and development. However, the precise mechanisms of DNA methylation remain poorly understood, especially in woody plants. We employed whole-genome bisulfite sequencing (WGBS), assays for transposase-accessible chromatin using sequencing (ATAC-seq), and RNA-Seq to investigate epigenetic regulatory relationships in Populus tomentosa treated with DNA methylation inhibitor 5-azacitidine. Expression-quantitative trait methylation analysis (eQTM), epigenome-wide association study (EWAS), and joint linkage-linkage disequilibrium mapping were used to explore the epigenetic regulatory genes, and using CRISPR/Cas9 to identify the role of candidate genes. Plant developmental abnormalities occurred when DNA methylation levels were substantially reduced. DNA methylation regulated 112 expressed genes via chromatin accessibility, of which 61 genes were significantly influenced by DNA methylation variation at the population level. One DNA methylation-regulated gene, PtoGntK, was located in a major quantitative trait locus (QTL) for poplar growth. Overexpression and CRISPR/Cas9 of PtoGntK revealed it affected poplar height and stem diameter. The PtoRAP2.12 was found to bind to the demethylated accessible region in the PtoGntK promoter, thereby promoting growth in poplar. This study identified key genes with epigenetic regulation for plant growth and provides insights into epigenetic regulation mechanisms in woody plants.

Exploring the versatility of zygotic genome regulators: A comparative and functional analysis

The activation of the zygotic genome constitutes an essential process during early embryogenesis that determines the overall progression of embryonic development. Zygotic genome activation (ZGA) is tightly regulated, involving a delicate interplay of activators and repressors, to precisely control the timing and spatial pattern of gene expression. While regulators of ZGA vary across species, they accomplish comparable outcomes. Recent studies have shed light on the unanticipated roles of ZGA components both during and after ZGA. Moreover, different ZGA regulators seem to have acquired unique functional modalities to manifest their regulatory potential. In this review, we explore these observations to assess whether these are simply anecdotal or contribute to a broader regulatory framework that employs a versatile means to arrive at the conserved outcome.

DNA methylation reprogramming in teleosts

Early embryonic development is crucially important but also remarkably diverse among animal taxa. Axis formation and cell lineage specification occur due to both spatial and temporal control of gene expression. This complex system involves various signaling pathways and developmental genes such as transcription factors as well as other molecular interactants that maintain cellular states, including several types of epigenetic marks. 5mC DNA methylation, the chemical modification of cytosines in eukaryotes, represents one such mark. By influencing the compaction of chromatin (a high-order DNA structure), DNA methylation can either repress or induce transcriptional activity. Mammals exhibit a reprogramming of DNA methylation from the parental genomes in the zygote following fertilization, and later in primordial germ cells (PGCs). Whether these periods of methylation reprogramming are evolutionarily conserved, or an innovation in mammals, is an emerging question. Looking into these processes in other vertebrate lineages is thus important, and teleost fish, with their extensive species richness, phenotypic diversity, and multiple rounds of whole genome duplication, provide the perfect research playground for answering such a question. This review aims to present a concise state of the art of DNA methylation reprogramming in early development in fish by summarizing findings from different research groups investigating methylation reprogramming patterns in teleosts, while keeping in mind the ramifications of the methodology used, then comparing those patterns to reprogramming patterns in mammals.

H3K4me2 distinguishes a distinct class of enhancers during the maternal-to-zygotic transition

After egg fertilization, an initially silent embryonic genome is transcriptionally activated during the maternal-to-zygotic transition. In zebrafish, maternal vertebrate pluripotency factors Nanog, Pou5f3 (OCT4 homolog), and Sox19b (SOX2 homolog) (NPS) play essential roles in orchestrating embryonic genome activation, acting as "pioneers" that open condensed chromatin and mediate acquisition of activating histone modifications. However, some embryonic gene transcription still occurs in the absence of these factors, suggesting the existence of other mechanisms regulating genome activation. To identify chromatin signatures of these unknown pathways, we profiled the histone modification landscape of zebrafish embryos using CUT&RUN. Our regulatory map revealed two subclasses of enhancers distinguished by presence or absence of H3K4me2. Enhancers lacking H3K4me2 tend to require NPS factors for de novo activation, while enhancers bearing H3K4me2 are epigenetically bookmarked by DNA hypomethylation to recapitulate gamete activity in the embryo, independent of NPS pioneering. Thus, parallel enhancer activation pathways combine to induce transcriptional reprogramming to pluripotency in the early embryo.

ULI-ssDRIP-seq revealed R-loop dynamics during vertebrate early embryogenesis

R-loop, a chromatin structure containing one RNA:DNA hybrid and one unpaired single-stranded DNA, plays multiple biological roles. However, due to technical limitations, the landscapes and potential functions of R-loops during embryogenesis remain elusive. Here, we developed a quantitative and high-resolution ultra-low input R-loop profiling method, named ULI-ssDRIP-seq, which can map global R-loops with as few as 1000 cells. By using ULI-ssDRIP-seq, we reveal the R-loop dynamics in the zebrafish from gametes to early embryos. In oocytes, the R-loop level is relatively low in most regions of the nuclear genome, except maternal-inherited rDNA and mitochondrial genome. The correlation between R-loop and CG methylation dynamics during early development is relatively weak. Furthermore, either up- or down-regulation of global R-loops by knockdown or overexpression of RNase H1 causes a delay of embryonic development with dramatic expression changes in zygotic and maternal genes. This study provides comprehensive R-loop landscapes during early vertebrate embryogenesis and demonstrates the implication of R-loops in embryonic development.

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

Rhabdomyosarcoma fusion oncoprotein initially pioneers a neural signature in vivo

Fusion-positive rhabdomyosarcoma is an aggressive pediatric cancer molecularly characterized by arrested myogenesis. The defining genetic driver, PAX3::FOXO1, functions as a chimeric gain-of-function transcription factor. An incomplete understanding of PAX3::FOXO1's in vivo epigenetic mechanisms has hindered therapeutic development. Here, we establish a PAX3::FOXO1 zebrafish injection model and semi-automated ChIP-seq normalization strategy to evaluate how PAX3::FOXO1 initially interfaces with chromatin in a developmental context. We investigated PAX3::FOXO1's recognition of chromatin and subsequent transcriptional consequences. We find that PAX3::FOXO1 interacts with inaccessible chromatin through partial/homeobox motif recognition consistent with pioneering activity. However, PAX3::FOXO1-genome binding through a composite paired-box/homeobox motif alters chromatin accessibility and redistributes H3K27ac to activate neural transcriptional programs. We uncover neural signatures that are highly representative of clinical rhabdomyosarcoma gene expression programs that are enriched following chemotherapy. Overall, we identify partial/homeobox motif recognition as a new mode for PAX3::FOXO1 pioneer function and identify neural signatures as a potentially critical PAX3::FOXO1 tumor initiation event.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

ADNP modulates SINE B2-derived CTCF-binding sites during blastocyst formation in mice

CTCF is crucial for chromatin structure and transcription regulation in early embryonic development. However, the kinetics of CTCF chromatin occupation in preimplantation embryos have remained unclear. In this study, we used CUT&RUN technology to investigate CTCF occupancy in mouse preimplantation development. Our findings revealed that CTCF begins binding to the genome prior to zygotic genome activation (ZGA), with a preference for CTCF-anchored chromatin loops. Although the majority of CTCF occupancy is consistently maintained, we identified a specific set of binding sites enriched in the mouse-specific short interspersed element (SINE) family B2 that are restricted to the cleavage stages. Notably, we discovered that the neuroprotective protein ADNP counteracts the stable association of CTCF at SINE B2-derived CTCF-binding sites. Knockout of Adnp in the zygote led to impaired CTCF binding signal recovery, failed deposition of H3K9me3, and transcriptional derepression of SINE B2 during the morula-to-blastocyst transition, which further led to unfaithful cell differentiation in embryos around implantation. Our analysis highlights an ADNP-dependent restriction of CTCF binding during cell differentiation in preimplantation embryos. Furthermore, our findings shed light on the functional importance of transposable elements (TEs) in promoting genetic innovation and actively shaping the early embryo developmental process specific to mammals.

miR-430 regulates zygotic mRNA during zebrafish embryogenesis

BACKGROUND

Early embryonic developmental programs are guided by the coordinated interplay between maternally inherited and zygotically manufactured RNAs and proteins. Although these processes happen concomitantly and affecting gene function during this period is bound to affect both pools of mRNAs, it has been challenging to study their expression dynamics separately.

RESULTS

By employing SLAM-seq, a nascent mRNA labeling transcriptomic approach, in a developmental time series we observe that over half of the early zebrafish embryo transcriptome consists of maternal-zygotic genes, emphasizing their pivotal role in early embryogenesis. We provide an hourly resolution of de novo transcriptional activation events and follow nascent mRNA trajectories, finding that most de novo transcriptional events are stable throughout this period. Additionally, by blocking microRNA-430 function, a key post transcriptional regulator during zebrafish embryogenesis, we directly show that it destabilizes hundreds of de novo transcribed mRNAs from pure zygotic as well as maternal-zygotic genes. This unveils a novel miR-430 function during embryogenesis, fine-tuning zygotic gene expression.

CONCLUSION

These insights into zebrafish early embryo transcriptome dynamics emphasize the significance of post-transcriptional regulators in zygotic genome activation. The findings pave the way for future investigations into the coordinated interplay between transcriptional and post-transcriptional landscapes required for the establishment of animal cell identities and functions.

Analysis of Chromatin Accessibility and DNA Methylation to Reveal the Functions of Epigenetic Modifications in Cyprinus carpio Gonads

Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.

The Novel Role of Zfp296 in Mammalian Embryonic Genome Activation as an H3K9me3 Modulator

The changes in epigenetic modifications during early embryonic development significantly impact mammalian embryonic genome activation (EGA) and are species-conserved to some degree. Here, we reanalyzed the published RNA-Seq of human, mouse, and goat early embryos and found that Zfp296 (zinc finger protein 296) expression was higher at the EGA stage than at the oocyte stage in all three species (adjusted p-value < 0.05 |log2(foldchange)| ≥ 1). Subsequently, we found that Zfp296 was conserved across human, mouse, goat, sheep, pig, and bovine embryos. In addition, we identified that ZFP296 interacts with the epigenetic regulators KDM5B, SMARCA4, DNMT1, DNMT3B, HP1β, and UHRF1. The Cys₂-His₂(C2H2) zinc finger domain TYPE2 TYPE3 domains of ZFP296 co-regulated the modification level of the trimethylation of lysine 9 on the histone H3 protein subunit (H3K9me3). According to ChIP-seq analysis, ZFP296 was also enriched in Trim28, Suv39h1, Setdb1, Kdm4a, and Ehmt2 in the mESC genome. Then, knockdown of the expression of Zfp296 at the late zygote of the mouse led to the early developmental arrest of the mouse embryos and failure resulting from a decrease in H3K9me3. Together, our results reveal that Zfp296 is an H3K9me3 modulator which is essential to the embryonic genome activation of mouse embryos.

The miR-430 locus with extreme promoter density forms a transcription body during the minor wave of zygotic genome activation

In anamniote embryos, the major wave of zygotic genome activation starts during the mid-blastula transition. However, some genes escape global genome repression, are activated substantially earlier, and contribute to the minor wave of genome activation. The mechanisms underlying the minor wave of genome activation are little understood. We explored the genomic organization and cis-regulatory mechanisms of a transcription body, in which the minor wave of genome activation is first detected in zebrafish. We identified the miR-430 cluster as having excessive copy number and the highest density of Pol-II-transcribed promoters in the genome, and this is required for forming the transcription body. However, this transcription body is not essential for, nor does it encompasse, minor wave transcription globally. Instead, distinct minor-wave-specific promoter architecture suggests that promoter-autonomous mechanisms regulate the minor wave of genome activation. The minor-wave-specific features also suggest distinct transcription initiation mechanisms between the minor and major waves of genome activation.

Epigenomic dynamics of early Xenopus Embryos

How the embryonic genome regulates accessibility to transcription factors is one of the major questions in understanding the spatial and temporal dynamics of gene expression during embryogenesis. Epigenomic analyses of embryonic chromatin provide molecular insights into cell-specific gene activities and genomic architectures. In recent years, significant advances have been made to elucidate the dynamic changes behind the activation of the zygotic genome in various model organisms. Here we provide an overview of the recent epigenomic studies pertaining to early Xenopus development.

MiR-202-3p determines embryo viability during mid-blastula transition

Developmental growth is an intricate process involving the coordinated regulation of the expression of various genes, and microRNAs (miRNAs) play crucial roles in diverse processes throughout animal development. The mid-blastula transition (MBT) is a developmental milestone when maternal RNAs are cleared and the zygotic genome programmed asynchronous cell division begins to drive embryogenesis. While mechanisms underlying MBT have been intensively revealed, factors regulating cell proliferation at the transition remain largely unknown. We report here a microRNA, miR-202-3p to be a key factor that determines embryonic fate during MBT in zebrafish. A miR-202-3p antagomir specifically terminated embryo development at the mid-blastula stage. In vivo deletion of the miR-202 locus recapitulated the fatal phenotypes, which were rescued only by miR-202-3p or its precursor. Transcriptome comparison revealed >250 RNAs including both maternal and zygotic origins were dysregulated at MBT in the miR-202^−/− embryos, corresponding with arrays of homeostatic disorders leading to massive apoptosis. A trio of genes: nfkbiaa, perp and mgll, known to be intimately involved with cell proliferation and survival, were identified as direct targets of miR-202-3p. Importantly, over- or under-expression of any of the trio led to developmental delay or termination at the blastula or gastrula stages. Furthermore, nfkbiaa and perp were shown to inter-regulate each other. Thus, miR-202-3p mediates a regulatory network whose components interact closely during MBT to determine embryonic viability and development.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

Comparative methylomics and chromatin accessibility analysis in uncovers regulation of genic transcription and mechanisms of key floral scent production

Linalool and ionone are two important aromatic components in sweet osmanthus petals, and the regulatory mechanisms that produce these two components remain unclear. In this study, we employed whole-genome methylation sequencing and ATAC-seq technology to analyze the genomic DNA methylation status and chromatin accessibility of the sweet osmanthus cultivars 'Zaohuang' and 'Chenghong Dangui'. Results showed that the promoter region of TPS2, a key gene in the linalool synthesis pathway, was less methylated in 'Chenghong Dangui' than in 'Zaohuang'. The chromatin was more accessible in 'Chenghong Dangui' than in 'Zaohuang', which resulted in a much stronger expression of this gene in 'Chenghong Dangui' than in 'Zaohuang'. This eventually led to a high quantity of linalool and its oxides in the petals of 'Chenghong Dangui', but there were lower levels present in the petals of 'Zaohuang'. These results suggest that DNA methylation and chromatin accessibility play major roles in linalool synthesis in sweet osmanthus. The methylation level of the promoter region of CCD4, a key gene for ionone synthesis, was higher in 'Zaohuang' than in 'Chenghong Dangui'. The chromatin accessibility was lower in 'Zaohuang' than in 'Chenghong Dangui', although the expression of this gene was significantly higher in 'Zaohuang' than in 'Chenghong Dangui'. ChIP-seq analysis and a series of experiments showed that the differential expression of CCD4 and CCD1 in the two cultivars may predominantly be the result of regulation by ERF2 and other transcription factors. However, a 183-bp deletion involving the CCD4 promoter region in 'Chenghong Dangui' may be the main reason for the low expression of this gene in its petals. This study provides an important theoretical basis for improving selective breeding of key floral fragrance components in sweet osmanthus.

The landscape of pioneer factor activity reveals the mechanisms of chromatin reprogramming and genome activation

Upon fertilization, embryos undergo chromatin reprogramming and genome activation; however, the mechanisms that regulate these processes are poorly understood. Here, we generated a triple mutant for Nanog, Pou5f3, and Sox19b (NPS) in zebrafish and found that NPS pioneer chromatin opening at >50% of active enhancers. NPS regulate acetylation across core histones at enhancers and promoters, and their function in gene activation can be bypassed by recruiting histone acetyltransferase to individual genes. NPS pioneer chromatin opening individually, redundantly, or additively depending on sequence context, and we show that high nucleosome occupancy facilitates NPS pioneering activity. Nucleosome position varies based on the input of different transcription factors (TFs), providing a flexible platform to modulate pioneering activity. Altogether, our results illuminate the sequence of events during genome activation and offer a conceptual framework to understand how pioneer factors interpret the genome and integrate different TF inputs across cell types and developmental transitions.

Primordial Germ Cell Development in the Poeciliid, Gambusia holbrooki, Reveals Shared Features Between Lecithotrophs and Matrotrophs

Metazoans exhibit two modes of primordial germ cell (PGC) specification that are interspersed across taxa. However, the evolutionary link between the two modes and the reproductive strategies of lecithotrophy and matrotrophy is poorly understood. As a first step to understand this, the spatio-temporal expression of teleostean germ plasm markers was investigated in Gambusia holbrooki, a poecilid with shared lecitho- and matrotrophy. A group of germ plasm components was detected in the ovum suggesting maternal inheritance mode of PGC specification. However, the strictly zygotic activation of dnd-β and nanos1 occurred relatively early, reminiscent of models with induction mode (e.g., mice). The PGC clustering, migration and colonisation patterns of G. holbrooki resembled those of zebrafish, medaka and mice at blastula, gastrula and somitogenesis, respectively-recapitulating features of advancing evolutionary nodes with progressive developmental stages. Moreover, the expression domains of PGC markers in G. holbrooki were either specific to teleost (vasa expression in developing PGCs), murine models (dnd spliced variants) or shared between the two taxa (germline and somatic expression of piwi and nanos1). Collectively, the results suggest that the reproductive developmental adaptations may reflect a transition from lecithotrophy to matrotrophy.

CpG island reconfiguration for the establishment and synchronization of polycomb functions upon exit from naive pluripotency

Polycomb group (PcG) proteins are essential for post-implantation development by depositing repressive histone modifications at promoters, mainly CpG islands (CGIs), of developmental regulator genes. However, promoter PcG marks are erased after fertilization and de novo established in peri-implantation embryos, coinciding with the transition from naive to primed pluripotency. Nevertheless, the molecular basis for this establishment remains unknown. In this study, we show that the expression of the long KDM2B isoform (KDM2BLF), which contains the demethylase domain, is specifically induced at peri-implantation and that its H3K36me2 demethylase activity is required for PcG enrichment at CGIs. Moreover, KDM2BLF interacts with BRG1/BRM-associated factor (BAF) and stabilizes BAF occupancy at CGIs for subsequent gain of accessibility, which precedes PcG enrichment. Consistently, KDM2BLF inactivation results in significantly delayed post-implantation development. In summary, our data unveil dynamic chromatin configuration of CGIs during exit from naive pluripotency and provide a conceptual framework for the spatiotemporal establishment of PcG functions.

Antibody-free profiling of transcription factor occupancy during early embryogenesis by FitCUT&RUN

Key transcription factors (TFs) play critical roles in zygotic genome activation (ZGA) during early embryogenesis, while genome-wide occupancies of only a few factors have been profiled during ZGA due to the limitation of cell numbers or the lack of high-quality antibodies. Here, we present FitCUT&RUN, a modified CUT&RUN method, in which an Fc fragment of immunoglobulin G is used for tagging, to profile TF occupancy in an antibody-free manner and demonstrate its reliability and robustness using as few as five thousand K562 cells. We applied FitCUT&RUN to zebrafish undergoing embryogenesis to generate reliable occupancy profiles of three known activators of zebrafish ZGA: Nanog, Pou5f3 and Sox19b. By profiling the time-series occupancy of Nanog during zebrafish ZGA, we observed a clear trend toward a gradual increase in Nanog occupancy and found that Nanog occupancy prior to the major phase of ZGA is critical for the activation of a significant proportion of early transcribed genes. Our results further suggested that the sequential binding of Nanog may be controlled by replication timing and the presence of Nanog motifs.

Promoter methylation-mediated repression of UNC5 receptors and the associated clinical significance in human colorectal cancer

Background

Deregulated methylation of tumor suppressor genes is a hallmark event in colorectal cancer (CRC) carcinogenesis. UNC5 receptors, down-regulated in various human malignancies due to epigenetic alterations, have been proposed as putative tumor suppressor genes. In this study, we focused on the methylation-mediated inhibition of UNC5 receptors and the associated clinical significance in CRC.

Methods

Methylation and expression analysis was performed in TCGA datasets. And the results were confirmed in vitro in CRC cell lines treated with 5-aza-deoxycytidine. Then, the expression and epigenetic alterations of UNC5 receptors were evaluated in clinical specimens. Moreover, the diagnostic and prognostic values of the methylation alterations were also analyzed.

Results

Methylation-mediated repression was observed in UNC5C and UNC5D, but not in UNC5A and UNC5B, which was confirmed in CRC cell lines. Except for UNC5B, significantly elevated methylation was observed in UNC5A, UNC5C, and UNC5D in CRC. The discrimination efficiency of the three receptors was comparable with that of SEPT9. Kaplan–Meier curve survival analysis showed that hypermethylation of UNC5A, UNC5C and UNC5D was associated with poor progression-free and overall survival. Moreover, methylation levels of UNC5C and UNC5D were independent predictors of CRC progression-free (P = 0.001, P = 0.003, respectively) and overall survival (P = 0.008, P = 0.004, respectively).

Conclusions

Hypermethylation of UNC5C and UNC5D mediates the repression and has promising diagnostic and prognostic values in CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01211-5.

Key features of the POU transcription factor Oct4 from an evolutionary perspective

Oct4, a class V POU-domain protein that is encoded by the Pou5f1 gene, is thought to be a key transcription factor in the early development of mammals. This transcription factor plays indispensable roles in pluripotent stem cells as well as in the acquisition of pluripotency during somatic cell reprogramming. Oct4 has also been shown to play a role as a pioneer transcription factor during zygotic genome activation (ZGA) from zebrafish to human. However, during the past decade, several studies have brought these conclusions into question. It was clearly shown that the first steps in mouse development are not affected by the loss of Oct4. Subsequently, the role of Oct4 as a genome activator was brought into doubt. It was also found that the reprogramming of somatic cells into induced pluripotent stem cells (iPSCs) could proceed without Oct4. In this review, we summarize recent findings, reassess the role of Oct4 in reprogramming and ZGA, and point to structural features that may underlie this role. We speculate that pluripotent stem cells resemble neural stem cells more closely than previously thought. Oct4 orthologs within the POUV class hold key roles in genome activation during early development of species with late ZGA. However, in Placentalia, eutherian-specific proteins such as Dux overtake Oct4 in ZGA and endow them with the formation of an evolutionary new tissue-the placenta.

Control of zygotic genome activation in Xenopus

The fertilized frog egg contains all the materials needed to initiate development of a new organism, including stored RNAs and proteins deposited during oogenesis, thus the earliest stages of development do not require transcription. The onset of transcription from the zygotic genome marks the first genetic switch activating the gene regulatory network that programs embryonic development. Zygotic genome activation occurs after an initial phase of transcriptional quiescence that continues until the midblastula stage, a period called the midblastula transition, which was first identified in Xenopus. Activation of transcription is programmed by maternally supplied factors and is regulated at multiple levels. A similar switch exists in most animals and is of great interest both to developmental biologists and to those interested in understanding nuclear reprogramming. Here we review in detail our knowledge on this major switch in transcription in Xenopus and place recent discoveries in the context of a decades old problem.

Pan-cancer chromatin analysis of the human vtRNA genes uncovers their association with cancer biology

Background: The vault RNAs (vtRNAs) are a class of 84-141-nt eukaryotic non-coding RNAs transcribed by RNA polymerase III, associated to the ribonucleoprotein complex known as vault particle. Of the four human vtRNA genes, vtRNA1-1, vtRNA1-2 and vtRNA1-3, clustered at locus 1, are integral components of the vault particle, while vtRNA2-1 is a more divergent homologue located in a second locus. Gene expression studies of vtRNAs in large cohorts have been hindered by their unsuccessful sequencing using conventional transcriptomic approaches.

Methods: VtRNA expression in The Cancer Genome Atlas (TCGA) Pan-Cancer cohort was estimated using the genome-wide DNA methylation and chromatin accessibility data (ATAC-seq) of their genes as surrogate variables. The association between vtRNA expression and patient clinical outcome, immune subtypes and transcriptionally co-regulated gene programs was analyzed in the dataset.

Results: VtRNAs promoters are enriched in transcription factors related to viral infection. VtRNA2-1 is likely the most independently regulated homologue. VtRNA1-1 has the most accessible chromatin, followed by vtRNA1-2, vtRNA2-1 and vtRNA1-3. VtRNA1-1 and vtRNA1-3 chromatin status does not significantly change in cancer tissues. Meanwhile, vtRNA2-1 and vtRNA1-2 expression is widely deregulated in neoplastic tissues and its alteration is compatible with a broad oncogenic role for vtRNA1-2, and both tumor suppressor and oncogenic functions for vtRNA2-1. Yet, vtRNA1-1, vtRNA1-2 and vtRNA2-1 promoter DNA methylation predicts a shorter patient overall survival cancer-wide. In addition, gene ontology analyses of vtRNAs co-regulated genes identify a chromosome regulatory domain, epithelial differentiation, immune and thyroid cancer gene sets for specific vtRNAs. Furthermore, vtRNA expression patterns are associated with cancer immune subtypes and vtRNA1-2 expression is positively associated with cell proliferation and wound healing.

Conclusions: Our study presents the landscape of vtRNA chromatin status cancer-wide, identifying co-regulated gene networks and ontological pathways associated with the different vtRNA genes that may account for their diverse roles in cancer.

Genome-wide variation in DNA methylation linked to developmental stage and chromosomal suppression of recombination in white-throated sparrows

Much of our knowledge on regulatory impacts of DNA methylation has come from laboratory‐bred model organisms, which may not exhibit the full extent of variation found in wild populations. Here, we investigated naturally‐occurring variation in DNA methylation in a wild avian species, the white‐throated sparrow (Zonotrichia albicollis). This species offers exceptional opportunities for studying the link between genetic differentiation and phenotypic traits because of a nonrecombining chromosome pair linked to both plumage and behavioural phenotypes. Using novel single‐nucleotide resolution methylation maps and gene expression data, we show that DNA methylation and the expression of DNA methyltransferases are significantly higher in adults than in nestlings. Genes for which DNA methylation varied between nestlings and adults were implicated in development and cell differentiation and were located throughout the genome. In contrast, differential methylation between plumage morphs was concentrated in the nonrecombining chromosome pair. Interestingly, a large number of CpGs on the nonrecombining chromosome, localized to transposable elements, have undergone dramatic loss of DNA methylation since the split of the ZAL2 and ZAL2^m chromosomes. Changes in methylation predicted changes in gene expression for both chromosomes. In summary, we demonstrate changes in genome‐wide DNA methylation that are associated with development and with specific functional categories of genes in white‐throated sparrows. Moreover, we observe substantial DNA methylation reprogramming associated with the suppression of recombination, with implications for genome integrity and gene expression divergence. These results offer an unprecedented view of ongoing epigenetic reprogramming in a wild population.

see also the Perspective by Jordan A. Anderson and Jenny Tung.

A map of cis-regulatory elements and 3D genome structures in zebrafish

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species¹. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

EGCG Attenuates Renal Damage via Reversing Klotho Hypermethylation in Diabetic db/db Mice and HK-2 Cells

To explore whether epigallocatechin-3-gallate (EGCG) improves renal damage in diabetic db/db mice and high-glucose- (HG-) induced injury in HK-2 cells by regulating the level of Klotho gene promoter methylation. Western blotting was used to detect the protein expression levels of DNA methyltransferase 1 (DNMT1), DNMT3a, DNMT3b, transforming growth factor-β1 (TGF-β1), α-smooth muscle actin (α-SMA), and Klotho. The methylation level of the Klotho gene promoter was detected by pyrosequencing. Chromatin immunoprecipitation was used to detect the binding of the Klotho gene promoter to DNMT1 and DNMT3a. The expression of oxidative stress markers (reactive oxygen species (ROS), superoxide dismutase (SOD), malondialdehyde (MDA), catalase (CAT), and 8-hydroxy-2′-deoxyguanosine (8-OHdG)) and inflammatory cytokines (interleukin-1β (IL-1β), IL-6, and tumor necrosis factor-α (TNF-α)) in kidney homogenates was also measured using ELISA. Klotho and DNMT3b protein expression was upregulated, while DNMT1, DNMT3a, TGF-β1, and α-SMA protein expression was downregulated after EGCG treatment. EGCG treatment also reduced the methylation level of the Klotho gene promoter as well as the binding of DNMT3a to the Klotho gene promoter. In addition, EGCG treatment significantly decreased the levels of ROS, MDA, 8-OHdG, IL-1β, IL-6, and TNF-α and increased the levels of CAT and SOD. Under HG conditions, EGCG regulated Klotho gene promoter methylation via DNMT3a and decreased the methylation level of the Klotho gene promoter, thereby upregulating the expression of the Klotho protein to exert its protective effect.

Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis

Medaka (Oryzias latipes) has become an important vertebrate model widely used in genetics, developmental biology, environmental sciences, and many other fields. A high-quality genome sequence and a variety of genetic tools are available for this model organism. However, existing genome annotation is still rudimentary, as it was mainly based on computational prediction and short-read RNA-seq data. Here we report a dynamic transcriptome landscape of medaka embryogenesis profiled by long-read RNA-seq, short-read RNA-seq, and ATAC-seq. By integrating these data sets, we constructed a much-improved gene model set including about 17,000 novel isoforms and identified 1600 transcription factors, 1100 long noncoding RNAs, and 150,000 potential cis-regulatory elements as well. Time-series data sets provided another dimension of information. With the expression dynamics of genes and accessibility dynamics of cis-regulatory elements, we investigated isoform switching, as well as regulatory logic between accessible elements and genes, during embryogenesis. We built a user-friendly medaka omics data portal to present these data sets. This resource provides the first comprehensive omics data sets of medaka embryogenesis. Ultimately, we term these three assays as the minimum ENCODE toolbox and propose the use of it as the initial and essential profiling genomic assays for model organisms that have limited data available. This work will be of great value for the research community using medaka as the model organism and many others as well.

Histones, Their Variants and Post-translational Modifications in Zebrafish Development

Complex multi-cellular organisms are shaped starting from a single-celled zygote, owing to elaborate developmental programs. These programs involve several layers of regulation to orchestrate the establishment of progressively diverging cell type-specific gene expression patterns. In this scenario, epigenetic modifications of chromatin are central in influencing spatiotemporal patterns of gene transcription. In fact, it is generally recognized that epigenetic changes of chromatin states impact on the accessibility of genomic DNA to regulatory proteins. Several lines of evidence highlighted that zebrafish is an excellent vertebrate model for research purposes in the field of developmental epigenetics. In this review, I focus on the dynamic roles recently emerged for histone post-translational modifications (PTMs), histone modifying enzymes, histone variants and histone themselves in the coordination between the precise execution of transcriptional programs and developmental progression in zebrafish. In particular, I first outline a synopsis of the current state of knowledge in this field during early embryogenesis. Then, I present a survey of histone-based epigenetic mechanisms occurring throughout morphogenesis, with a stronger emphasis on cardiac formation. Undoubtedly, the issues addressed in this review take on particular importance in the emerging field of comparative biology of epigenetics, as well as in translational research.

Chromatin dynamics at the maternal to zygotic transition: recent advances from the zebrafish model

Early animal development is characterized by intense reorganization of the embryonic genome, including large-scale changes in chromatin structure and in the DNA and histone modifications that help shape this structure. Particularly profound shifts in the chromatin landscape are associated with the maternal-to-zygotic transition, when the zygotic genome is first transcribed and maternally loaded transcripts are degraded. The accessibility of the early zebrafish embryo facilitates the interrogation of chromatin during this critical window of development, making it an important model for early chromatin regulation. Here, we review our current understanding of chromatin dynamics during early zebrafish development, highlighting new advances as well as similarities and differences between early chromatin regulation in zebrafish and other species.

From mother to embryo: A molecular perspective on zygotic genome activation

In animals, the early embryo is mostly transcriptionally silent and development is fueled by maternally supplied mRNAs and proteins. These maternal products are important not only for survival, but also to gear up the zygote's genome for activation. Over the last three decades, research with different model organisms and experimental approaches has identified molecular factors and proposed mechanisms for how the embryo transitions from being transcriptionally silent to transcriptionally competent. In this chapter, we discuss the molecular players that shape the molecular landscape of ZGA and provide insights into their mode of action in activating the transcription program in the developing embryo.

Chromatin accessibility established by Pou5f3, Sox19b and Nanog primes genes for activity during zebrafish genome activation

In many organisms, early embryonic development is driven by maternally provided factors until the controlled onset of transcription during zygotic genome activation. The regulation of chromatin accessibility and its relationship to gene activity during this transition remain poorly understood. Here, we generated chromatin accessibility maps with ATAC-seq from genome activation until the onset of lineage specification. During this period, chromatin accessibility increases at regulatory elements. This increase is independent of RNA polymerase II-mediated transcription, with the exception of the hypertranscribed miR-430 locus. Instead, accessibility often precedes the transcription of associated genes. Loss of the maternal transcription factors Pou5f3, Sox19b, and Nanog, which are known to be required for zebrafish genome activation, results in decreased accessibility at regulatory elements. Importantly, the accessibility of regulatory regions, especially when established by Pou5f3, Sox19b and Nanog, is predictive for future transcription. Our results show that the maternally provided transcription factors Pou5f3, Sox19b, and Nanog open up chromatin and prime genes for activity during zygotic genome activation in zebrafish.

DNA methylation stability in fish spermatozoa upon external constraint: Impact of fish hormonal stimulation and sperm cryopreservation

Highly differentiated mature spermatozoa carry not only genetic but also epigenetic information that is to be transmitted to the embryo. DNA methylation is one epigenetic actor associated with sperm nucleus compaction, gene silencing, and prepatterning of embryonic gene expression. Therefore, the stability of this mark toward reproductive biotechnologies is a major issue in animal production. The present work explored the impact of hormonal induction of spermiation and sperm cryopreservation in two cyprinids, the goldfish (Carassius auratus) and the zebrafish (Danio rerio), using LUminometric Methylation Assay (LUMA). We showed that while goldfish hormonal treatment did increase sperm production, it did not alter global DNA methylation of spermatozoa. Different sperm samples repeatedly collected from the same males for 2 months also showed the same global DNA methylation level. Similarly, global DNA methylation was not affected after cryopreservation of goldfish spermatozoa with methanol, whereas less efficient cryoprotectants (dimethylsulfoxide and 1,2-propanediol) decreased DNA methylation. In contrast, cryopreservation of zebrafish spermatozoa with methanol induced a slight, but significant, increase in global DNA methylation. In the less compact nuclei, that is, goldfish fin somatic cells, cryopreservation did not change global DNA methylation regardless of the choice of cryoprotectant. To conclude, global DNA methylation is a robust parameter with respect to biotechnologies such as hormonal induction of spermiation and sperm cryopreservation, but it can be altered when the best sperm manipulation conditions are not met.

Mechanisms regulating zygotic genome activation

Following fertilization, the two specified gametes must unite to create an entirely new organism. The genome is initially transcriptionally quiescent, allowing the zygote to be reprogrammed into a totipotent state. Gradually, the genome is activated through a process known as the maternal-to-zygotic transition, which enables zygotic gene products to replace the maternal supply that initiated development. This essential transition has been broadly characterized through decades of research in several model organisms. However, we still lack a full mechanistic understanding of how genome activation is executed and how this activation relates to the reprogramming of the zygotic chromatin architecture. Recent work highlights the central role of transcriptional activators and suggests that these factors may coordinate transcriptional activation with other developmental changes.

The Spatiotemporal Control of Zygotic Genome Activation

One of the earliest and most significant events in embryonic development is zygotic genome activation (ZGA). In several species, bulk transcription begins at the midblastula transition (MBT) when, after a certain number of cleavages, the embryo attains a particular nuclear-to-cytoplasmic (N/C) ratio, maternal repressors become sufficiently diluted, and the cell cycle slows down. Here we resolve the frog ZGA in time and space by profiling RNA polymerase II (RNAPII) engagement and its transcriptional readout. We detect a gradual increase in both the quantity and the length of RNAPII elongation before the MBT, revealing that >1,000 zygotic genes disregard the N/C timer for their activation and that the sizes of newly transcribed genes are not necessarily constrained by cell cycle duration. We also find that Wnt, Nodal, and BMP signaling together generate most of the spatiotemporal dynamics of regional ZGA, directing the formation of orthogonal body axes and proportionate germ layers.

The maternal-to-zygotic transition revisited

The development of animal embryos is initially directed by maternal gene products. Then, during the maternal-to-zygotic transition (MZT), developmental control is handed to the zygotic genome. Extensive research in both vertebrate and invertebrate model organisms has revealed that the MZT can be subdivided into two phases, during which very different modes of gene regulation are implemented: initially, regulation is exclusively post-transcriptional and post-translational, following which gradual activation of the zygotic genome leads to predominance of transcriptional regulation. These changes in the gene expression program of embryos are precisely controlled and highly interconnected. Here, we review current understanding of the mechanisms that underlie handover of developmental control during the MZT.

Widespread Enhancer Dememorization and Promoter Priming during Parental-to-Zygotic Transition

The epigenome plays critical roles in controlling gene expression and development. However, how the parental epigenomes transit to the zygotic epigenome in early development remains elusive. Here we show that parental-to-zygotic transition in zebrafish involves extensive erasure of parental epigenetic memory, starting with methylating gametic enhancers. Surprisingly, this occurs even prior to fertilization for sperm. Both parental enhancers lose histone marks by the 4-cell stage, and zygotic enhancers are not activated until around zygotic genome activation (ZGA). By contrast, many promoters remain hypomethylated and, unexpectedly, acquire histone acetylation before ZGA at as early as the 4-cell stage. They then resolve into either activated or repressed promoters upon ZGA. Maternal depletion of histone acetyltransferases results in aberrant ZGA and early embryonic lethality. Finally, such reprogramming is largely driven by maternal factors, with zygotic products mainly contributing to embryonic enhancer activation. These data reveal widespread enhancer dememorization and promoter priming during parental-to-zygotic transition.

Causal Gene Regulatory Network Modeling and Genomics: Second-Generation Challenges

Gene regulatory network modeling has played a major role in advancing the understanding of developmental systems, by crystallizing structures of relevant extant information, by formally posing hypothetical functional relationships between network elements, and by offering clear predictive tests to improve understanding of the mechanisms driving developmental progression. Both ordinary differential equation (ODE)-based and Boolean models have also been highly successful in explaining dynamics within subcircuits of more complex processes. In a very small number of cases, gene regulatory network models of much more global scope have been proposed that successfully predict the dynamics of the processes establishing most of an embryonic body plan. Can such successes be expanded to very different developmental systems, including post-embryonic mammalian systems? This perspective discusses several problems that must be solved in more quantitative and predictive theoretical terms, to make this possible. These problems include: the effects of cellular history on chromatin state and how these affect gene accessibility; the dose dependence of activities of many transcription factors (a problem for Boolean models); stochasticity of some transcriptional outputs (a problem for simple ODE models); response timing delays due to epigenetic remodeling requirements; functionally different kinds of repression; and the regulatory syntax that governs responses of genes with multiple enhancers.

Packaging development: how chromatin controls transcription in zebrafish embryogenesis

How developmental gene expression is activated, co-ordinated and maintained is one of the biggest questions in developmental biology. While transcription factors lead the way in directing developmental gene expression, their accessibility to the correct repertoire of genes can depend on other factors such as DNA methylation, the presence of particular histone variants and post-translational modifications of histones. Collectively, factors that modify DNA or affect its packaging and accessibility contribute to a chromatin landscape that helps to control the timely expression of developmental genes. Zebrafish, perhaps better known for their strength as a model of embryology and organogenesis during development, are coming to the fore as a powerful model for interpreting the role played by chromatin in gene expression. Several recent advances have shown that zebrafish exhibit both similarities and differences to other models (and humans) in the way that they employ chromatin mechanisms of gene regulation. Here, I review how chromatin influences developmental transcriptional programmes during early zebrafish development, patterning and organogenesis. Lastly, I briefly highlight the importance of zebrafish chromatin research towards the understanding of human disease and transgenerational inheritance.

Continued Activity of the Pioneer Factor Zelda Is Required to Drive Zygotic Genome Activation

Reprogramming cell fate during the first stages of embryogenesis requires that transcriptional activators gain access to the genome and remodel the zygotic transcriptome. Nonetheless, it is not clear whether the continued activity of these pioneering factors is required throughout zygotic genome activation or whether they are only required early to establish cis-regulatory regions. To address this question, we developed an optogenetic strategy to rapidly and reversibly inactivate the master regulator of genome activation in Drosophila, Zelda. Using this strategy, we demonstrate that continued Zelda activity is required throughout genome activation. We show that Zelda binds DNA in the context of nucleosomes and suggest that this allows Zelda to occupy the genome despite the rapid division cycles in the early embryo. These data identify a powerful strategy to inactivate transcription factor function during development and suggest that reprogramming in the embryo may require specific, continuous pioneering functions to activate the genome.

Pou5f3, SoxB1, and Nanog remodel chromatin on high nucleosome affinity regions at zygotic genome activation

The zebrafish embryo is transcriptionally mostly quiescent during the first 10 cell cycles, until the main wave of zygotic genome activation (ZGA) occurs, accompanied by fast chromatin remodeling. At ZGA, homologs of the mammalian stem cell transcription factors (TFs) Pou5f3, Nanog, and Sox19b bind to thousands of developmental enhancers to initiate transcription. So far, how these TFs influence chromatin dynamics at ZGA has remained unresolved. To address this question, we analyzed nucleosome positions in wild-type and maternal-zygotic (MZ) mutants for pou5f3 and nanog by MNase-seq. We show that Nanog, Sox19b, and Pou5f3 bind to the high nucleosome affinity regions (HNARs). HNARs are spanning over 600 bp, featuring high in vivo and predicted in vitro nucleosome occupancy and high predicted propeller twist DNA shape value. We suggest a two-step nucleosome destabilization-depletion model, in which the same intrinsic DNA properties of HNAR promote both high nucleosome occupancy and differential binding of TFs. In the first step, already before ZGA, Pou5f3 and Nanog destabilize nucleosomes at HNAR centers genome-wide. In the second step, post-ZGA, Nanog, Pou5f3, and SoxB1 maintain open chromatin state on the subset of HNARs, acting synergistically. Nanog binds to the HNAR center, whereas the Pou5f3 stabilizes the flanks. The HNAR model will provide a useful tool for genome regulatory studies in a variety of biological systems.

Evolution of embryonic cis-regulatory landscapes between divergent Phallusia and Ciona ascidians

Ascidian species of the Phallusia and Ciona genera are distantly related, their last common ancestor dating several hundred million years ago. Although their genome sequences have extensively diverged since this radiation, Phallusia and Ciona species share almost identical early morphogenesis and stereotyped cell lineages. Here, we explored the evolution of transcriptional control between P. mammillata and C. robusta. We combined genome-wide mapping of open chromatin regions in both species with a comparative analysis of the regulatory sequences of a test set of 10 pairs of orthologous early regulatory genes with conserved expression patterns. We find that ascidian chromatin accessibility landscapes obey similar rules as in other metazoa. Open-chromatin regions are short, highly conserved within each genus and cluster around regulatory genes. The dynamics of chromatin accessibility and closest-gene expression are strongly correlated during early embryogenesis. Open-chromatin regions are highly enriched in cis-regulatory elements: 73% of 49 open chromatin regions around our test genes behaved as either distal enhancers or proximal enhancer/promoters following electroporation in Phallusia eggs. Analysis of this datasets suggests a pervasive use in ascidians of "shadow" enhancers with partially overlapping activities. Cross-species electroporations point to a deep conservation of both the trans-regulatory logic between these distantly-related ascidians and the cis-regulatory activities of individual enhancers. Finally, we found that the relative order and approximate distance to the transcription start site of open chromatin regions can be conserved between Ciona and Phallusia species despite extensive sequence divergence, a property that can be used to identify orthologous enhancers, whose regulatory activity can partially diverge.