The earliest transcribed zygotic genes are short, newly evolved, and different across species

Parent-of-origin regulation by maternal auts2 shapes neurodevelopment and behavior in fish

BACKGROUND

Parental experience can influence progeny behavior through gamete-mediated non-genetic inheritance, that is, mechanisms that do not involve changes in inherited DNA sequence. However, underlying mechanisms remain poorly understood in vertebrates, especially for maternal effects. Here, we use the medaka, a model fish species, to investigate the role of auts2a, the ortholog of human AUTS2, a gene repressed in the fish oocyte following maternal stress and associated with neurodevelopmental disorders.

RESULTS

We show that auts2a expression in the oocyte influences long-term progeny behavior, including anxiety-like behavior and environment recognition capabilities. Using single-nuclei RNA-sequencing, we reveal that maternal auts2a influences gene expression in neural cell populations during neurodevelopment. We also show that maternal auts2a knock-out triggers differences in maternally inherited factors, including early embryonic transcriptional and post-transcriptional regulators.

CONCLUSIONS

Together, our results reveal the unsuspected role of an autism-related gene expressed in the mother's oocyte in shaping progeny neurodevelopment and behavior. Finally, we report that auts2a/AUTS2 is part of a group of evolutionarily conserved genes associated with human neurodevelopmental disorders and expressed in oocytes across species, from fish to mammals. These findings raise important questions about their potential role in the non-genetic regulation of progeny neurodevelopment and behavior in vertebrates.

Zygotic activation of transposable elements during zebrafish early embryogenesis

Although previous studies have shown that transposable elements (TEs) are conservatively activated to play key roles during early embryonic development, the details of zygotic TE activation (ZTA) remain poorly understood. Here, we employ long-read sequencing to precisely identify that only a small subset of TE loci are activated among numerous copies, allowing us to map their hierarchical transcriptional cascades at the single-locus and single-transcript level. Despite the heterogeneity of ZTA across family, subfamily, locus, and transcript levels, our findings reveal that ZTA follows a markedly different pattern from conventional zygotic gene activation (ZGA): ZTA occurs significantly later than ZGA and shows a pronounced bias for nuclear localization of TE transcripts. This study advances our understanding of TE activation by providing a high-resolution view of TE copies and creating a comprehensive catalog of thousands of previously unannotated transcripts and genes that are activated during early zebrafish embryogenesis. Among these genes, we highlight two that are essential for zebrafish development.

The maternal-to-zygotic transition: reprogramming of the cytoplasm and nucleus

A fertilized egg is initially transcriptionally silent and relies on maternally provided factors to initiate development. For embryonic development to proceed, the oocyte-inherited cytoplasm and the nuclear chromatin need to be reprogrammed to create a permissive environment for zygotic genome activation (ZGA). During this maternal-to-zygotic transition (MZT), which is conserved in metazoans, transient totipotency is induced and zygotic transcription is initiated to form the blueprint for future development. Recent technological advances have enhanced our understanding of MZT regulation, revealing common themes across species and leading to new fundamental insights about transcription, mRNA decay and translation.

Muskelin is a substrate adaptor of the highly regulated Drosophila embryonic CTLH E3 ligase

The maternal-to-zygotic transition (MZT) is a conserved developmental process where the maternally-derived protein and mRNA cache is replaced with newly made zygotic gene products. We have previously shown that in Drosophila the deposited RNA-binding proteins ME31B, Cup, and Trailer Hitch are ubiquitylated by the CTLH E3 ligase and cleared. However, the organization and regulation of the CTLH complex remain poorly understood in flies because Drosophila lacks an identifiable substrate adaptor, and the mechanisms restricting the degradation of ME31B and its cofactors to the MZT are unknown. Here, we show that the developmental regulation of the CTLH complex is multi-pronged, including transcriptional control by OVO and autoinhibition of the E3 ligase. One major regulatory target is the subunit Muskelin, which we demonstrate is a substrate adaptor for the Drosophila CTLH complex. Finally, we find that Muskelin has few targets beyond the three known RNA-binding proteins, showing exquisite target specificity. Thus, multiple levels of integrated regulation restrict the activity of the embryonic CTLH complex to early embryogenesis, during which time it regulates three important RNA-binding proteins.

Quantifying Nascent Transcription in Early Embryogenesis

The early embryonic genome exists in a dormant state following fertilization, and it then subsequently undergoes broad activation of zygotic transcription at the early stages of development. A major challenge is the detection of newly made zygotic transcripts and the determination of their activation onset time due to the presence of large and predominantly maternal pool of RNAs. Here we describe a detailed method to measure the zygotic transcription during zygotic genome activation (ZGA) of Xenopus early embryos using metabolic labeling of nascent transcripts with 5-ethynyl uridine (5-EU) followed by purifying and sequencing the nascent EU-RNAs (EU-RNA-seq). This method is highly sensitive in detecting early zygotic transcripts that are not detected by total RNA-seq and determines the actual onset time of transcriptional activation for zygotic genes. The method is applicable to a wide variety of embryonic model systems and has already afforded novel insights into gene regulation in early embryogenesis.

Rise and SINE: roles of transcription factors and retrotransposons in zygotic genome activation

In sexually reproducing organisms, life begins with the fusion of transcriptionally silent gametes, the oocyte and sperm. Although initiation of transcription in the embryo, known as zygotic genome activation (ZGA), is universally required for development, the transcription factors regulating this process are poorly conserved. In this Perspective, we discuss recent insights into the mechanisms of ZGA in totipotent mammalian embryos, namely ZGA regulation by several transcription factors, including by orphan nuclear receptors (OrphNRs) such as the pioneer transcription factor NR5A2, and by factors of the DUX, TPRX and OBOX families. We performed a meta-analysis and compiled a list of pan-ZGA genes, and found that most of these genes are indeed targets of the above transcription factors. Remarkably, more than a third of these ZGA genes appear to be regulated both by OrphNRs such as NR5A2 and by OBOX proteins, whose motifs co-occur in SINE B1 retrotransposable elements, which are enriched near ZGA genes. We propose that ZGA in mice is activated by recruitment of multiple transcription factors to SINE B1 elements that function as enhancers, and discuss a potential relevance of this mechanism to Alu retrotransposable elements in human ZGA.

Unveiling axolotl transcriptome for tissue regeneration with high-resolution annotation via long-read sequencing

Axolotls are known for their remarkable regeneration ability. Exploring their transcriptome provides insight into regenerative mechanisms. However, the current annotation of the axolotl transcriptome is limited, leaving the role of unannotated transcripts in regeneration unknown. To discourse this challenge, we exploited long-read sequencing technology, which enables direct observation of full-length RNA transcripts, greatly enhancing the coverage and accuracy of axolotl transcriptome annotation. By utilizing this method, we identified 222 novel gene loci and 4775 novel transcripts, which were quantified using short-read sequencing data. Through the inclusive analysis, we discovered novel homologs, potential functional proteins, noncoding RNAs, and alternative splicing events in key regeneration pathways. In particular, we identified novel transcripts with high protein-coding potential implicated in cell cycle regulation and musculoskeletal development, and regeneration were identified. Interestingly, alternative splice variants were also detected across diverse pathways critical to regeneration. This specifies that these novel transcripts potentially play vital roles underpinning the robust regenerative capacities of axolotls. Single-cell transcriptomic analysis further revealed these isoforms to predominantly exist in axolotl limb chondrocytes and mature tissue cell populations. Overall, the findings significantly advanced consideration of the axolotl transcriptome and provided a new perspective for understanding the mechanisms of regenerative abilities of axolotls.

The dynamics and functional impact of tRNA repertoires during early embryogenesis in zebrafish

Embryogenesis entails dramatic shifts in mRNA translation and turnover that reprogram gene expression during cellular proliferation and differentiation. Codon identity modulates mRNA stability during early vertebrate embryogenesis, but how the composition of tRNA pools is matched to translational demand is unknown. By quantitative profiling of tRNA repertoires in zebrafish embryos during the maternal-to-zygotic transition, we show that zygotic tRNA repertoires are established after the onset of gastrulation, succeeding the major wave of zygotic mRNA transcription. Maternal and zygotic tRNA pools are distinct, but their reprogramming does not result in a better match to the codon content of the zygotic transcriptome. Instead, we find that an increase in global translation at gastrulation sensitizes decoding rates to tRNA supply, thus destabilizing maternal mRNAs enriched in slowly translated codons. Translational activation and zygotic tRNA expression temporally coincide with an increase of TORC1 activity at gastrulation, which phosphorylates and inactivates the RNA polymerase III repressor Maf1a/b. Our data indicate that a switch in global translation, rather than tRNA reprogramming, determines the onset of codon-dependent maternal mRNA decay during zebrafish embryogenesis.

Protein profiling of zebrafish embryos unmasks regulatory layers during early embryogenesis

The maternal-to-zygotic transition is crucial in embryonic development, marked by the degradation of maternally provided mRNAs and initiation of zygotic gene expression. However, the changes occurring at the protein level during this transition remain unclear. Here, we conducted protein profiling throughout zebrafish embryogenesis using quantitative mass spectrometry, integrating transcriptomics and translatomics datasets. Our data show that, unlike RNA changes, protein changes are less dynamic. Further, increases in protein levels correlate with mRNA translation, whereas declines in protein levels do not, suggesting active protein degradation processes. Interestingly, proteins from pure zygotic genes are present at fertilization, challenging existing mRNA-based gene classifications. As a proof of concept, we utilized CRISPR-Cas13d to target znf281b mRNA, a gene whose protein significantly accumulates within the first 2 h post-fertilization, demonstrating its crucial role in development. Consequently, our protein profiling, coupled with CRISPR-Cas13d, offers a complementary approach to unraveling maternal factor function during embryonic development.

Epigenetic inheritance and gene expression regulation in early Drosophila embryos

Precise spatiotemporal regulation of gene expression is of paramount importance for eukaryotic development. The maternal-to-zygotic transition (MZT) during early embryogenesis in Drosophila involves the gradual replacement of maternally contributed mRNAs and proteins by zygotic gene products. The zygotic genome is transcriptionally activated during the first 3 hours of development, in a process known as "zygotic genome activation" (ZGA), by the orchestrated activities of a few pioneer factors. Their decisive role during ZGA has been characterized in detail, whereas the contribution of chromatin factors to this process has been historically overlooked. In this review, we aim to summarize the current knowledge of how chromatin regulation impacts the first stages of Drosophila embryonic development. In particular, we will address the following questions: how chromatin factors affect ZGA and transcriptional silencing, and how genome architecture promotes the integration of these processes early during development. Remarkably, certain chromatin marks can be intergenerationally inherited, and their presence in the early embryo becomes critical for the regulation of gene expression at later stages. Finally, we speculate on the possible roles of these chromatin marks as carriers of epialleles during transgenerational epigenetic inheritance (TEI).

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.

Muskelin acts as a substrate receptor of the highly regulated Drosophila CTLH E3 ligase during the maternal-to-zygotic transition

The maternal-to-zygotic transition (MZT) is a conserved developmental process where the maternally-derived protein and mRNA cache is replaced with newly made zygotic gene products. We have previously shown that in Drosophila the deposited RNA-binding proteins ME31B, Cup, and Trailer Hitch (TRAL) are ubiquitylated by the CTLH E3 ligase and cleared. However, the organization and regulation of the CTLH complex remain poorly understood in flies. In particular, Drosophila lacks an identifiable substrate adaptor, and the mechanisms restricting degradation of ME31B and its cofactors to the MZT are unknown. Here, we show that the developmental specificity of the CTLH complex is mediated by multipronged regulation, including transcriptional control by the transcription factor OVO and autoinhibition of the E3 ligase. One major regulatory target is the subunit Muskelin, which we demonstrate acts as a substrate adaptor for the Drosophila CTLH complex. Although conserved, Muskelin has structural roles in other species, suggesting a surprising functional plasticity. Finally, we find that Muskelin has few targets beyond the three known RNA binding proteins, showing exquisite target specificity. Thus, multiple levels of integrated regulation restrict the activity of the embryonic CTLH complex to early embryogenesis, seemingly with the goal of regulating three important RNA binding proteins.

Quantitative proteomics reveals the dynamic proteome landscape of zebrafish embryos during the maternal-to-zygotic transition

Maternal-to-zygotic transition (MZT) is central to early embryogenesis. However, its underlying molecular mechanisms are still not well described. Here, we revealed the expression dynamics of 5,000 proteins across four stages of zebrafish embryos during MZT, representing one of the most systematic surveys of proteome landscape of the zebrafish embryos during MZT. Nearly 700 proteins were differentially expressed and were divided into six clusters according to their expression patterns. The proteome expression profiles accurately reflect the main events that happen during the MZT, i.e., zygotic genome activation (ZGA), clearance of maternal mRNAs, and initiation of cellular differentiation and organogenesis. MZT is modulated by many proteins at multiple levels in a collaborative fashion, i.e., transcription factors, histones, histone-modifying enzymes, RNA helicases, and P-body proteins. Significant discrepancies were discovered between zebrafish proteome and transcriptome profiles during the MZT. The proteome dynamics database will be a valuable resource for bettering our understanding of MZT.

CRISPR-RfxCas13d screening uncovers Bckdk as a post-translational regulator of the maternal-to-zygotic transition in teleosts

The Maternal-to-Zygotic transition (MZT) is a reprograming process encompassing zygotic genome activation (ZGA) and the clearance of maternally-provided mRNAs. While some factors regulating MZT have been identified, there are thousands of maternal RNAs whose function has not been ascribed yet. Here, we have performed a proof-of-principle CRISPR-RfxCas13d maternal screening targeting mRNAs encoding protein kinases and phosphatases in zebrafish and identified Bckdk as a novel post-translational regulator of MZT. Bckdk mRNA knockdown caused epiboly defects, ZGA deregulation, H3K27ac reduction and a partial impairment of miR-430 processing. Phospho-proteomic analysis revealed that Phf10/Baf45a, a chromatin remodeling factor, is less phosphorylated upon Bckdk depletion. Further, phf10 mRNA knockdown also altered ZGA and Phf10 constitutively phosphorylated rescued the developmental defects observed after bckdk mRNA depletion. Altogether, our results demonstrate the competence of CRISPR-RfxCas13d screenings to uncover new regulators of early vertebrate development and shed light on the post-translational control of MZT mediated by protein phosphorylation.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

Comparative Analysis of Maternal Gene Expression Patterns Unravels Evolutionary Signatures Across Reproductive Modes

Maternal genes have a pivotal role in regulating metazoan early development. As such their functions have been extensively studied since the dawn of developmental biology. The temporal and spatial dynamics of their transcripts have been thoroughly described in model organisms and their functions have been undergoing heavy investigations. Yet, less is known about the evolutionary changes shaping their presence within diverse oocytes. Due to their unique maternal inheritance pattern, a high degree is predicted to be present when it comes to their expression. Insofar only limited and conflicting results have emerged around it. Here, we set out to elucidate which evolutionary changes could be detected in the maternal gene expression patterns using phylogenetic comparative methods on RNAseq data from 43 species. Using normalized gene expression values and fold change information throughout early development we set out to find the best-fitting evolutionary model. Through modeling, we find evidence supporting both the high degree of divergence and constraint on gene expression values, together with their temporal dynamics. Furthermore, we find that maternal gene expression alone can be used to explain the reproductive modes of different species. Together, these results suggest a highly dynamic evolutionary landscape of maternal gene expression. We also propose a possible functional dichotomy of maternal genes which is influenced by the reproductive strategy undertaken by examined species.

Detection of newly synthesized RNA reveals transcriptional reprogramming during ZGA and a role of Obox3 in totipotency acquisition

Zygotic genome activation (ZGA) after fertilization enables the maternal-to-zygotic transition. However, the global view of ZGA, particularly at initiation, is incompletely understood. Here, we develop a method to capture and sequence newly synthesized RNA in early mouse embryos, providing a view of transcriptional reprogramming during ZGA. Our data demonstrate that major ZGA gene activation begins earlier than previously thought. Furthermore, we identify a set of genes activated during minor ZGA, the promoters of which show enrichment of the Obox factor motif, and find that Obox3 or Obox5 overexpression in mouse embryonic stem cells activates ZGA genes. Notably, the expression of Obox factors is severely impaired in somatic cell nuclear transfer (SCNT) embryos, and restoration of Obox3 expression corrects the ZGA profile and greatly improves SCNT embryo development. Hence, our study reveals dynamic transcriptional reprogramming during ZGA and underscores the crucial role of Obox3 in facilitating totipotency acquisition.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Transcription bodies regulate gene expression by sequestering CDK9

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells. It remains unclear, however, if and how transcription bodies affect gene expression. Here we disrupted the formation of two prominent endogenous transcription bodies that mark the onset of zygotic transcription in zebrafish embryos and analysed the effect on gene expression using enriched SLAM-seq and live-cell imaging. We find that the disruption of transcription bodies results in the misregulation of hundreds of genes. Here we focus on genes that are upregulated. These genes have accessible chromatin and are poised to be transcribed in the presence of the two transcription bodies, but they do not go into elongation. Live-cell imaging shows that disruption of the two large transcription bodies enables these poised genes to be transcribed in ectopic transcription bodies, suggesting that the large transcription bodies sequester a pause release factor. Supporting this hypothesis, we find that CDK9-the kinase that releases paused polymerase II-is highly enriched in the two large transcription bodies. Overexpression of CDK9 in wild-type embryos results in the formation of ectopic transcription bodies and thus phenocopies the removal of the two large transcription bodies. Taken together, our results show that transcription bodies regulate transcription by sequestering machinery, thereby preventing genes elsewhere in the nucleus from being transcribed.

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.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting the hypothesis that the Chr4R transcriptome might be different from the rest of the genome. We conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes and identified four regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brain and liver validated few transcripts from Region-2 in somatic cells, but without sex-specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. In Region-2, protein-coding genes lack human orthologs; it has zinc finger genes expressed early in zygotic genome activation; it has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and an distinct set of repetitive elements. The colocalization of 1) genes silenced in ovaries but not in testes that are 2) expressed in embryos briefly at the onset of zygotic genome activation; 3) maternal-specific genes for translation machinery; 4) maternal-specific spliceosome components; and 4) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a Maternal-to-Zygotic-Transition Gene Regulatory Block.

Effects of Insulin on Proliferation, Apoptosis, and Ferroptosis in Primordial Germ Cells via PI3K-AKT-mTOR Signaling Pathway

Primordial germ cells (PGCs) are essential for the genetic modification, resource conservation, and recovery of endangered breeds in chickens and need to remain viable and proliferative in vitro. Therefore, there is an urgent need to elucidate the functions of the influencing factors and their regulatory mechanisms. In this study, PGCs collected from Rugao yellow chicken embryonic eggs at Day 5.5 were cultured in media containing 0, 5, 10, 20, 50, and 100 μg/mL insulin. The results showed that insulin regulates cell proliferation in PGCs in a dose-dependent way, with an optimal dose of 10 μg/mL. Insulin mediates the mRNA expression of cell cycle-, apoptosis-, and ferroptosis-related genes. Insulin at 50 μg/mL and 100 μg/mL slowed down the proliferation with elevated ion content and GSH/oxidized glutathione (GSSG) in PGCs compared to 10 μg/mL. In addition, insulin activates the PI3K/AKT/mTOR pathway dose dependently. Collectively, this study demonstrates that insulin reduces apoptosis and ferroptosis and enhances cell proliferation in a dose-dependent manner via the PI3K-AKT-mTOR signaling pathway in PGCs, providing a new addition to the theory of the regulatory role of the growth and proliferation of PGC in vitro cultures.

Activator-blocker model of transcriptional regulation by pioneer-like factors

Zygotic genome activation (ZGA) in the development of flies, fish, frogs and mammals depends on pioneer-like transcription factors (TFs). Those TFs create open chromatin regions, promote histone acetylation on enhancers, and activate transcription. Here, we use the panel of single, double and triple mutants for zebrafish genome activators Pou5f3, Sox19b and Nanog, multi-omics and mathematical modeling to investigate the combinatorial mechanisms of genome activation. We show that Pou5f3 and Nanog act differently on synergistic and antagonistic enhancer types. Pou5f3 and Nanog both bind as pioneer-like TFs on synergistic enhancers, promote histone acetylation and activate transcription. Antagonistic enhancers are activated by binding of one of these factors. The other TF binds as non-pioneer-like TF, competes with the activator and blocks all its effects, partially or completely. This activator-blocker mechanism mutually restricts widespread transcriptional activation by Pou5f3 and Nanog and prevents premature expression of late developmental regulators in the early embryo.

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

miR-430 microRNA Family in Fishes: Molecular Characterization and Evolution

The miR-430 microRNA family has been described in multiple fish species as one of the first microRNAs expressed by the zygote. It has been suggested that this family is implicated in maternal mRNA elimination, but may also play a role in steroidogenesis, sexual differentiation, and flatfish metamorphosis. The miR-430 sequences have been found in multiple-copy tandem clusters but evidence of their conservation outside of teleost fishes is scarce. In the present study, we have characterized the tandem repeats organization of these microRNAs in different fish species, both model and of interest in aquaculture. A phylogenetic analysis of this family has allowed us to identify that the miR-430 duplication, which took place before the Chondrostei and Neopterygii groups' divergence, has resulted in three variants ("a", "b", and "c"). According to our data, variant "b" is the most closely related to the ancestral sequence. Furthermore, we have detected isolated instances of the miR-430 repeat subunit in some species, which suggests that this microRNA family may be affected by DNA rearrangements. This study provides new data about the abundance, variability, and organization of the miR-430 family in fishes.

Chromatin expansion microscopy reveals nanoscale organization of transcription and chromatin

Nanoscale chromatin organization regulates gene expression. Although chromatin is notably reprogrammed during zygotic genome activation (ZGA), the organization of chromatin regulatory factors during this universal process remains unclear. In this work, we developed chromatin expansion microscopy (ChromExM) to visualize chromatin, transcription, and transcription factors in vivo. ChromExM of embryos during ZGA revealed how the pioneer factor Nanog interacts with nucleosomes and RNA polymerase II (Pol II), providing direct visualization of transcriptional elongation as string-like nanostructures. Blocking elongation led to more Pol II particles clustered around Nanog, with Pol II stalled at promoters and Nanog-bound enhancers. This led to a new model termed "kiss and kick", in which enhancer-promoter contacts are transient and released by transcriptional elongation. Our results demonstrate that ChromExM is broadly applicable to study nanoscale nuclear organization.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus. Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities similar to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far particular to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Actin filaments accumulated in the nucleus remain in the vicinity of condensing chromosomes in the zebrafish early embryo

In the cytoplasm, filamentous actin (F-actin) plays a critical role in cell regulation, including cell migration, stress fiber formation, and cytokinesis. Recent studies have shown that actin filaments that form in the nucleus are associated with diverse functions. Here, using live imaging of an F-actin-specific probe, superfolder GFP-tagged utrophin (UtrCH-sfGFP), we demonstrated the dynamics of nuclear actin in zebrafish (Danio rerio) embryos. In early zebrafish embryos up to around the high stage, UtrCH-sfGFP increasingly accumulated in nuclei during the interphase and reached a peak during the prophase. After nuclear envelope breakdown (NEBD), patches of UtrCH-sfGFP remained in the vicinity of condensing chromosomes during the prometaphase to metaphase. When zygotic transcription was inhibited by injecting α-amanitin, the nuclear accumulation of UtrCH-sfGFP was still observed at the sphere and dome stages, suggesting that zygotic transcription may induce a decrease in nuclear F-actin. The accumulation of F-actin in nuclei may contribute to proper mitotic progression of large cells with rapid cell cycles in zebrafish early embryos, by assisting in NEBD, chromosome congression, and/or spindle assembly.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus . Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities analogous to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far unique to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Intron dynamics reveal principles of gene regulation during the maternal-to-zygotic transition

The maternal-to-zygotic transition (MZT) is a conserved embryonic process in animals where developmental control shifts from the maternal to zygotic genome. A key step in this transition is zygotic transcription, and deciphering the MZT requires classifying newly transcribed genes. However, due to current technological limitations, this starting point remains a challenge for studying many species. Here, we present an alternative approach that characterizes transcriptome changes based solely on RNA-seq data. By combining intron-mapping reads and transcript-level quantification, we characterized transcriptome dynamics during the Drosophila melanogaster MZT. Our approach provides an accessible platform to investigate transcriptome dynamics that can be applied to the MZT in nonmodel organisms. In addition to classifying zygotically transcribed genes, our analysis revealed that over 300 genes express different maternal and zygotic transcript isoforms due to alternative splicing, polyadenylation, and promoter usage. The vast majority of these zygotic isoforms have the potential to be subject to different regulatory control, and over two-thirds encode different proteins. Thus, our analysis reveals an additional layer of regulation during the MZT, where new zygotic transcripts can generate additional proteome diversity.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Stress-sensitive dynamics of miRNAs and Elba1 in Drosophila embryogenesis

Early-life stress can result in life-long effects that impact adult health and disease risk, but little is known about how such programming is established and maintained. Here, we show that such epigenetic memories can be initiated in the Drosophila embryo before the major wave of zygotic transcription, and higher-order chromatin structures are established. An early short heat shock results in elevated levels of maternal miRNA and reduced levels of a subgroup of zygotic genes in stage 5 embryos. Using a Dicer-1 mutant, we show that the stress-induced decrease in one of these genes, the insulator-binding factor Elba1, is dependent on functional miRNA biogenesis. Reduction in Elba1 correlates with the upregulation of early developmental genes and promotes a sustained weakening of heterochromatin in the adult fly as indicated by an increased expression of the PEV w^m4h reporter. We propose that maternal miRNAs, retained in response to an early embryonic heat shock, shape the subsequent de novo heterochromatin establishment that occurs during early development via direct or indirect regulation of some of the earliest expressed genes, including Elba1.

SLAMseq resolves the kinetics of maternal and zygotic gene expression during early zebrafish embryogenesis

The maternal-to-zygotic transition (MZT) is a key developmental process in metazoan embryos that involves the activation of zygotic transcription (ZGA) and degradation of maternal transcripts. We employed metabolic mRNA sequencing (SLAMseq) to deconvolute the compound embryonic transcriptome in zebrafish. While mitochondrial zygotic transcripts prevail prior to MZT, we uncover the spurious transcription of hundreds of short and intron-poor genes as early as the 2-cell stage. Upon ZGA, most zygotic transcripts originate from thousands of maternal-zygotic (MZ) genes that are transcribed at rates comparable to those of hundreds of purely zygotic genes and replenish maternal mRNAs at distinct timescales. Rapid replacement of MZ transcripts involves transcript decay features unrelated to major maternal degradation pathways and promotes de novo synthesis of the core gene expression machinery by increasing poly(A)-tail length and translation efficiency. SLAMseq hence provides insights into the timescales, molecular features, and regulation of MZT during zebrafish embryogenesis.

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.

Combined modelling of mRNA decay dynamics and single-molecule imaging in the Drosophila embryo uncovers a role for P-bodies in 5' to 3' degradation

Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. Many methods for determining mRNA half-lives rely on transcriptional inhibition or metabolic labelling. Here, we use a non-invasive method for estimating half-lives for hundreds of mRNAs in the early Drosophila embryo. This approach uses the intronic and exonic reads from a total RNA-seq time series and Gaussian process regression to model the dynamics of premature and mature mRNAs. We show how regulation of mRNA stability is used to establish a range of mature mRNA dynamics during embryogenesis, despite shared transcription profiles. Using single-molecule imaging, we provide evidence that, for the mRNAs tested, there is a correlation between short half-life and mRNA association with P-bodies. Moreover, we detect an enrichment of mRNA 3' ends in P-bodies in the early embryo, consistent with 5' to 3' degradation occurring in P-bodies for at least a subset of mRNAs. We discuss our findings in relation to recently published data suggesting that the primary function of P-bodies in other biological contexts is mRNA storage.

Nanog organizes transcription bodies

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells.^1-3 How proteins of the transcriptional machinery come together to form such bodies, however, is unclear. Here, we take advantage of two large, isolated, and long-lived transcription bodies that reproducibly form during early zebrafish embryogenesis to characterize the dynamics of transcription body formation. Once formed, these transcription bodies are enriched for initiating and elongating RNA polymerase II, as well as the transcription factors Nanog and Sox19b. Analyzing the events leading up to transcription, we find that Nanog and Sox19b cluster prior to transcription. The clustering of transcription factors is sequential; Nanog clusters first, and this is required for the clustering of Sox19b and the initiation of transcription. Mutant analysis revealed that both the DNA-binding domain as well as one of the two intrinsically disordered regions of Nanog are required to organize the two bodies of transcriptional activity. Taken together, our data suggest that the clustering of transcription factors dictates the formation of transcription bodies.

Comprehensive maturity of nuclear pore complexes regulates zygotic genome activation

Nuclear pore complexes (NPCs) are channels for nucleocytoplasmic transport of proteins and RNAs. However, it remains unclear whether composition, structure, and permeability of NPCs dynamically change during the cleavage period of vertebrate embryos and affect embryonic development. Here, we report that the comprehensive NPC maturity (CNM) controls the onset of zygotic genome activation (ZGA) during zebrafish early embryogenesis. We show that more nucleoporin proteins are recruited to and assembled into NPCs with development, resulting in progressive increase of NPCs in size and complexity. Maternal transcription factors (TFs) transport into nuclei more efficiently with increasing CNM. Deficiency or dysfunction of Nup133 or Ahctf1/Elys impairs NPC assembly, maternal TFs nuclear transport, and ZGA onset, while nup133 overexpression promotes these processes. Therefore, CNM may act as a molecular timer for ZGA by controlling nuclear transport of maternal TFs that reach nuclear concentration thresholds at a given time to initiate ZGA.

Nascent transcriptome reveals orchestration of zygotic genome activation in early embryogenesis

Early embryo development requires maternal-to-zygotic transition, during which transcriptionally silent nuclei begin widespread gene expression during zygotic genome activation (ZGA).1-3 ZGA is vital for early cell fating and germ-layer specification,3,4 and ZGA timing is regulated by multiple mechanisms.1-5 However, controversies remain about whether these mechanisms are interrelated and vary among species6-10 and whether the timing of germ-layer-specific gene activation is temporally ordered.11,12 In some embryonic models, widespread ZGA onset is spatiotemporally graded,13,14 yet it is unclear whether the transcriptome follows this pattern. A major challenge in addressing these questions is to accurately measure the timing of each gene activation. Here, we metabolically label and identify the nascent transcriptome using 5-ethynyl uridine (5-EU) in Xenopus blastula embryos. We find that EU-RNA-seq outperforms total RNA-seq in detecting the ZGA transcriptome, which is dominated by transcription from maternal-zygotic genes, enabling improved ZGA timing determination. We uncover discrete spatiotemporal patterns for individual gene activation, a majority following a spatial pattern of ZGA that is correlated with a cell size gradient.14 We further reveal that transcription necessitates a period of developmental progression and that ZGA can be precociously induced by cycloheximide, potentially through elongation of interphase. Finally, most ectodermal genes are activated earlier than endodermal genes, suggesting a temporal orchestration of germ-layer-specific genes, potentially linked to the spatially graded pattern of ZGA. Together, our study provides fundamental new insights into the composition and dynamics of the ZGA transcriptome, mechanisms regulating ZGA timing, and its role in the onset of early cell fating.

Natural variation in the maternal and zygotic mRNA complements of the early embryo in Drosophila melanogaster

Background

Maternal gene products supplied to the egg during oogenesis drive the earliest events of development in all metazoans. After the initial stages of embryogenesis, maternal transcripts are degraded as zygotic transcription is activated; this is known as the maternal to zygotic transition (MZT). Recently, it has been shown that the expression of maternal and zygotic transcripts have evolved in the Drosophila genus over the course of 50 million years. However, the extent of natural variation of maternal and zygotic transcripts within a species has yet to be determined. We asked how the maternal and zygotic pools of mRNA vary within and between populations of D. melanogaster. In order to maximize sampling of genetic diversity, African lines of D. melanogaster originating from Zambia as well as DGRP lines originating from North America were chosen for transcriptomic analysis.

Results

Generally, we find that maternal transcripts are more highly conserved, and zygotic transcripts evolve at a higher rate. We find that there is more within-population variation in transcript abundance than between populations and that expression variation is highest post- MZT between African lines.

Conclusions

Determining the natural variation of gene expression surrounding the MZT in natural populations of D. melanogaster gives insight into the extent of how a tightly regulated process may vary within a species, the extent of developmental constraint at both stages and on both the maternal and zygotic genomes, and reveals expression changes allowing this species to adapt as it spread across the world.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08839-4.

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.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

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.

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.

Volume segregation programming in a nematode's early embryogenesis

Nematode species are well-known for their invariant cell lineage pattern during development. Combining knowledge about the fate specification induced by asymmetric division and the anti-correlation between cell cycle length and cell volume in Caenorhabditis elegans, we propose a minimal model to simulate lineage initiation by altering cell volume segregation ratio in each division, and quantify the derived pattern's performance in proliferation speed, fate diversity, and space robustness. The stereotypic pattern in C. elegans embryo is found to be one of the most optimal solutions taking minimum time to achieve the cell number before gastrulation, by programming asymmetric divisions as a strategy.

Chromatin remodelling and retrotransposons activities during regeneration in Drosophila

Regeneration is a response mechanism aiming to reconstruct lost or damaged structures. To achieve this, the cells repopulating the lost tissue often have to change their original identity, a process that involves chromatin remodelling.We have analysed the issue of chromatin remodelling during regeneration in the wing disc of Drosophila . In this disc the ablation of the central region (the pouch) induces the regenerative response of the cells from the lateral region (the hinge), which reconstitute the wing pouch. We have examined euchromatin and heterochromatin histone marks during the process and find that heterochromatin marks disappear but are recovered when regeneration is complete. Euchromatin marks are not modified. We also describe the transcription of two retrotransposons, Roo and F-element in the regenerating cells. We have established a temporal correlation between the alterations of heterochromatin marks and the levels of transcription of two retrotransposons, Roo and F-element, both during embryonic development and in the regeneration process.

Conserved germ plasm characteristics across the Danio and Devario lineages

In many animal species, germ cell specification requires the inheritance of germ plasm, a biomolecular condensate containing maternally derived RNAs and proteins. Most studies of germ plasm composition and function have been performed in widely evolutionarily divergent model organisms, such as Caenorhabditis elegans, Drosophila, Xenopus laevis, and Danio rerio (zebrafish). In zebrafish, 12 RNAs localize to germ plasm at the furrows of the early embryo. Here, we tested for the presence of these RNAs in three additional species within the Danionin clade: Danio kyathit, Danio albolineatus, and Devario aequipinnatus. By visualizing nanos RNA, we find that germ plasm segregation patterns during early embryogenesis are conserved across these species. Ten additional germ plasm RNAs exhibit localization at the furrows of early embryos in all three non-zebrafish Danionin species, consistent with germ plasm localization. One component of zebrafish germ plasm, ca15b, lacked specific localization in embryos of the more distantly related D. aequipinnatus. Our findings show that within a subset of closely related Danionin species, the vast majority of germ plasm RNA components are conserved. At the same time, the lack of ca15b localization in D. aequipinnatus germ plasm highlights the potential for the divergence of germ plasm composition across a restricted phylogenetic space.

Temporal Transcriptome Analysis Reveals Dynamic Expression Profiles of Gametes and Embryonic Development in Japanese Flounder (Paralichthys olivaceus)

The maternal-to-zygotic transition (MZT) is a crucial event in embryo development. While the features of the MZT across species are shared, the stage of this transition is different among species. We characterized MZT in a flatfish species, Japanese flounder (Paralichthys olivaceus). In this study, we analyzed the 551.57 GB transcriptome data of two types of gametes (sperms and eggs) and 10 embryo developmental stages in Japanese flounder. We identified 2512 maternal factor-related genes and found that most of those maternal factor-related genes expression decreased at the low blastula (LB) stage and remained silent in the subsequent embryonic development period. Meanwhile, we verified that the zygotic genome transcription might occur at the 128-cell stage and large-scale transcription began at the LB stage, which indicates the LB stage is the major wave zygotic genome activation (ZGA) occurs. In addition, we indicated that the Wnt signaling pathway, playing a diverse role in embryonic development, was involved in the ZGA and the axis formation. The results reported the list of the maternal genes in Japanese flounder and defined the stage of MZT, contributing to the understanding of the details of MZT during Japanese flounder embryonic development.