Zygotic genome activation during the maternal-to-zygotic transition.

Embryogenesis depends on a highly coordinated cascade of genetically encoded events. In animals, maternal factors contributed by the egg cytoplasm initially control development, whereas the zygotic nuclear genome is quiescent. Subsequently, the genome is activated, embryonic gene products are mobilized, and maternal factors are cleared. This transfer of developmental control is called the maternal-to-zygotic transition (MZT). In this review, we discuss recent advances toward understanding the scope, timing, and mechanisms that underlie zygotic genome activation at the MZT in animals. We describe high-throughput techniques to measure the embryonic transcriptome and explore how regulation of the cell cycle, chromatin, and transcription factors together elicits specific patterns of embryonic gene expression. Finally, we illustrate the interplay between zygotic transcription and maternal clearance and show how these two activities combine to reprogram two terminally differentiated gametes into a totipotent embryo.

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.

Corrigendum: Cloche is a bHLH-PAS transcription factor that drives haemato-vascular specification.

This corrects the article DOI: 10.1038/nature18614.

Zinc protection of fertilized eggs is an ancient feature of sexual reproduction in animals.

One of the earliest and most prevalent barriers to successful reproduction is polyspermy, or fertilization of an egg by multiple sperm. To prevent these supernumerary fertilizations, eggs have evolved multiple mechanisms. It has recently been proposed that zinc released by mammalian eggs at fertilization may block additional sperm from entering. Here, we demonstrate that eggs from amphibia and teleost fish also release zinc. Using Xenopus laevis as a model, we document that zinc reversibly blocks fertilization. Finally, we demonstrate that extracellular zinc similarly disrupts early embryonic development in eggs from diverse phyla, including Cnidaria, Echinodermata, and Chordata. Our study reveals that a fundamental strategy protecting human eggs from fertilization by multiple sperm may have evolved more than 650 million years ago.

Optimized design of antisense oligomers for targeted rRNA depletion.

RNA sequencing (RNA-seq) is extensively used to quantify gene expression transcriptome-wide. Although often paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA), many applications require alternate approaches to counteract the high proportion of ribosomal RNA (rRNA) in total RNA. Recently, digestion using RNaseH and antisense DNA oligomers tiling target rRNAs has emerged as an alternative to commercial rRNA depletion kits. Here, we present a streamlined, more economical RNaseH-mediated rRNA depletion with substantially lower up-front costs, using shorter antisense oligos only sparsely tiled along the target RNA in a 5-min digestion reaction. We introduce a novel Web tool, Oligo-ASST, that simplifies oligo design to target regions with optimal thermodynamic properties, and additionally can generate compact, common oligo pools that simultaneously target divergent RNAs, e.g. across different species. We demonstrate the efficacy of these strategies by generating rRNA-depletion oligos for Xenopus laevis and for zebrafish, which expresses two distinct versions of rRNAs during embryogenesis. The resulting RNA-seq libraries reduce rRNA to <5% of aligned reads, on par with poly(A) selection, and also reveal expression of many non-adenylated RNA species. Oligo-ASST is freely available at https://mtleelab.pitt.edu/oligo to design antisense oligos for any taxon or to target any abundant RNA for depletion.

Genome wide analysis of 3' UTR sequence elements and proteins regulating mRNA stability during maternal-to-zygotic transition in zebrafish.

Posttranscriptional regulation plays a crucial role in shaping gene expression. During the maternal-to-zygotic transition (MZT), thousands of maternal transcripts are regulated. However, how different cis-elements and trans-factors are integrated to determine mRNA stability remains poorly understood. Here, we show that most transcripts are under combinatorial regulation by multiple decay pathways during zebrafish MZT. By using a massively parallel reporter assay, we identified cis-regulatory sequences in the 3' UTR, including U-rich motifs that are associated with increased mRNA stability. In contrast, miR-430 target sequences, UAUUUAUU AU-rich elements (ARE), CCUC, and CUGC elements emerged as destabilizing motifs, with miR-430 and AREs causing mRNA deadenylation upon genome activation. We identified trans-factors by profiling RNA-protein interactions and found that poly(U)-binding proteins are preferentially associated with 3' UTR sequences and stabilizing motifs. We show that this activity is antagonized by C-rich motifs and correlated with protein binding. Finally, we integrated these regulatory motifs into a machine learning model that predicts reporter mRNA stability in vivo.

Cloche is a bHLH-PAS transcription factor that drives haemato-vascular specification.

Vascular and haematopoietic cells organize into specialized tissues during early embryogenesis to supply essential nutrients to all organs and thus play critical roles in development and disease. At the top of the haemato-vascular specification cascade lies cloche, a gene that when mutated in zebrafish leads to the striking phenotype of loss of most endothelial and haematopoietic cells and a significant increase in cardiomyocyte numbers. Although this mutant has been analysed extensively to investigate mesoderm diversification and differentiation and continues to be broadly used as a unique avascular model, the isolation of the cloche gene has been challenging due to its telomeric location. Here we used a deletion allele of cloche to identify several new cloche candidate genes within this genomic region, and systematically genome-edited each candidate. Through this comprehensive interrogation, we succeeded in isolating the cloche gene and discovered that it encodes a PAS-domain-containing bHLH transcription factor, and that it is expressed in a highly specific spatiotemporal pattern starting during late gastrulation. Gain-of-function experiments show that it can potently induce endothelial gene expression. Epistasis experiments reveal that it functions upstream of etv2 and tal1, the earliest expressed endothelial and haematopoietic transcription factor genes identified to date. A mammalian cloche orthologue can also rescue blood vessel formation in zebrafish cloche mutants, indicating a highly conserved role in vertebrate vasculogenesis and haematopoiesis. The identification of this master regulator of endothelial and haematopoietic fate enhances our understanding of early mesoderm diversification and may lead to improved protocols for the generation of endothelial and haematopoietic cells in vivo and in vitro.

RESA identifies mRNA-regulatory sequences at high resolution.

Gene expression is extensively regulated at the levels of mRNA stability, localization and translation. However, decoding functional RNA-regulatory features remains a limitation to understanding post-transcriptional regulation in vivo. Here, we developed RNA-element selection assay (RESA), a method that selects RNA elements on the basis of their activity in vivo and uses high-throughput sequencing to provide a quantitative measurement of their regulatory functions at near-nucleotide resolution. We implemented RESA to identify sequence elements modulating mRNA stability during zebrafish embryogenesis. RESA provides a sensitive and quantitative measure of microRNA activity in vivo and also identifies novel regulatory sequences. To uncover specific sequence requirements within regulatory elements, we developed a bisulfite-mediated nucleotide-conversion strategy for large-scale mutational analysis (RESA-bisulfite). Finally, we used the versatile RESA platform to map candidate protein-RNA interactions in vivo (RESA-CLIP).

Nanog, Pou5f1 and SoxB1 activate zygotic gene expression during the maternal-to-zygotic transition.

After fertilization, maternal factors direct development and trigger zygotic genome activation (ZGA) at the maternal-to-zygotic transition (MZT). In zebrafish, ZGA is required for gastrulation and clearance of maternal messenger RNAs, which is in part regulated by the conserved microRNA miR-430. However, the factors that activate the zygotic program in vertebrates are unknown. Here we show that Nanog, Pou5f1 (also called Oct4) and SoxB1 regulate zygotic gene activation in zebrafish. We identified several hundred genes directly activated by maternal factors, constituting the first wave of zygotic transcription. Ribosome profiling revealed that nanog, sox19b and pou5f1 are the most highly translated transcription factors pre-MZT. Combined loss of these factors resulted in developmental arrest before gastrulation and a failure to activate >75% of zygotic genes, including miR-430. Our results demonstrate that maternal Nanog, Pou5f1 and SoxB1 are required to initiate the zygotic developmental program and induce clearance of the maternal program by activating miR-430 expression.

Identification of small ORFs in vertebrates using ribosome footprinting and evolutionary conservation.

Identification of the coding elements in the genome is a fundamental step to understanding the building blocks of living systems. Short peptides (< 100 aa) have emerged as important regulators of development and physiology, but their identification has been limited by their size. We have leveraged the periodicity of ribosome movement on the mRNA to define actively translated ORFs by ribosome footprinting. This approach identifies several hundred translated small ORFs in zebrafish and human. Computational prediction of small ORFs from codon conservation patterns corroborates and extends these findings and identifies conserved sequences in zebrafish and human, suggesting functional peptide products (micropeptides). These results identify micropeptide-encoding genes in vertebrates, providing an entry point to define their function in vivo.

Divergence of RNA localization between rat and mouse neurons reveals the potential for rapid brain evolution.

BACKGROUND: Neurons display a highly polarized architecture. Their ability to modify their features under intracellular and extracellular stimuli, known as synaptic plasticity, is a key component of the neurochemical basis of learning and memory. A key feature of synaptic plasticity involves the delivery of mRNAs to distinct sub-cellular domains where they are locally translated. Regulatory coordination of these spatio-temporal events is critical for synaptogenesis and synaptic plasticity as defects in these processes can lead to neurological diseases. In this work, using microdissected dendrites from primary cultures of hippocampal neurons of two mouse strains (C57BL/6 and Balb/c) and one rat strain (Sprague-Dawley), we investigate via microarrays, subcellular localization of mRNAs in dendrites of neurons to assay the evolutionary differences in subcellular dendritic transcripts localization. RESULTS: Our microarray analysis highlighted significantly greater evolutionary diversification of RNA localization in the dendritic transcriptomes (81% gene identity difference among the top 5% highly expressed genes) compared to the transcriptomes of 11 different central nervous system (CNS) and non-CNS tissues (average of 44% gene identity difference among the top 5% highly expressed genes). Differentially localized genes include many genes involved in CNS function. CONCLUSIONS: Species differences in sub-cellular localization may reflect non-functional neutral drift. However, the functional categories of mRNA showing differential localization suggest that at least part of the divergence may reflect activity-dependent functional differences of neurons, mediated by species-specific RNA subcellular localization mechanisms.

Intron retention facilitates splice variant diversity in calcium-activated big potassium channel populations.

We report that the stress axis-regulated exon (STREX)-containing calcium-activated big potassium (BKCa) channel splice variant expression and physiology are regulated in part by cytoplasmic splicing and intron retention. NextGen sequencing of the mRNA complement of pooled hippocampal dendrite samples found intron 17a (i17a), the intron immediately preceding STREX, in the BKCa mRNA. Further molecular analyses of i17a revealed that the majority of i17a-containing BKCa channel mRNAs associate with STREX. i17a siRNA treatment followed by STREX protein immunocytochemistry demonstrated both reduced levels and altered subcellular distribution of STREX-containing BKCa channel protein. Selective reduction of i17a-BKCa or STREX-BKCa mRNAs induced similar changes in the burst firing properties of hippocampal neurons. Collectively, these data show that STREX splice variant regulation via cytoplasmic splicing and intron retention helps generate STREX-dependent BKCa current diversity in hippocampal neurons.

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.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis.

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Transcriptome transfer produces a predictable cellular phenotype.

Cellular phenotype is the conglomerate of multiple cellular processes involving gene and protein expression that result in the elaboration of a cell's particular morphology and function. It has been thought that differentiated postmitotic cells have their genomes hard wired, with little ability for phenotypic plasticity. Here we show that transfer of the transcriptome from differentiated rat astrocytes into a nondividing differentiated rat neuron resulted in the conversion of the neuron into a functional astrocyte-like cell in a time-dependent manner. This single-cell study permits high resolution of molecular and functional components that underlie phenotype identity. The RNA population from astrocytes contains RNAs in the appropriate relative abundances that give rise to regulatory RNAs and translated proteins that enable astrocyte identity. When transferred into the postmitotic neuron, the astrocyte RNA population converts 44% of the neuronal host cells into the destination astrocyte-like phenotype. In support of this observation, quantitative measures of cellular morphology, single-cell PCR, single-cell microarray, and single-cell functional analyses have been performed. The host-cell phenotypic changes develop over many weeks and are persistent. We call this process of RNA-induced phenotype changes, transcriptome-induced phenotype remodeling.

RNA degradation is required for the germ-cell to maternal transition in Drosophila.

In sexually reproducing animals, the oocyte contributes a large supply of RNAs that are essential to launch development upon fertilization. The mechanisms that regulate the composition of the maternal RNA contribution during oogenesis are unclear. Here, we show that a subset of RNAs expressed during the early stages of oogenesis is subjected to regulated degradation during oocyte specification. Failure to remove these RNAs results in oocyte dysfunction and death. We identify the RNA-degrading Super Killer complex and No-Go Decay factor Pelota as key regulators of oogenesis via targeted degradation of specific RNAs expressed in undifferentiated germ cells. These regulators target RNAs enriched for cytidine sequences that are bound by the polypyrimidine tract binding protein Half pint. Thus, RNA degradation helps orchestrate a germ cell-to-maternal transition that gives rise to the maternal contribution to the zygote.

Self containment, a property of modular RNA structures, distinguishes microRNAs.

RNA molecules will tend to adopt a folded conformation through the pairing of bases on a single strand; the resulting so-called secondary structure is critical to the function of many types of RNA. The secondary structure of a particular substring of functional RNA may depend on its surrounding sequence. Yet, some RNAs such as microRNAs retain their specific structures during biogenesis, which involves extraction of the substructure from a larger structural context, while other functional RNAs may be composed of a fusion of independent substructures. Such observations raise the question of whether particular functional RNA substructures may be selected for invariance of secondary structure to their surrounding nucleotide context. We define the property of self containment to be the tendency for an RNA sequence to robustly adopt the same optimal secondary structure regardless of whether it exists in isolation or is a substring of a longer sequence of arbitrary nucleotide content. We measured degree of self containment using a scoring method we call the self-containment index and found that miRNA stem loops exhibit high self containment, consistent with the requirement for structural invariance imposed by the miRNA biogenesis pathway, while most other structured RNAs do not. Further analysis revealed a trend toward higher self containment among clustered and conserved miRNAs, suggesting that high self containment may be a characteristic of novel miRNAs acquiring new genomic contexts. We found that miRNAs display significantly enhanced self containment compared to other functional RNAs, but we also found a trend toward natural selection for self containment in most functional RNA classes. We suggest that self containment arises out of selection for robustness against perturbations, invariance during biogenesis, and modular composition of structural function. Analysis of self containment will be important for both annotation and design of functional RNAs. A Python implementation and Web interface to calculate the self-containment index are available at http://kim.bio.upenn.edu/software/.

The Paf1 Complex Broadly Impacts the Transcriptome of Saccharomyces cerevisiae.

The Polymerase Associated Factor 1 complex (Paf1C) is a multifunctional regulator of eukaryotic gene expression important for the coordination of transcription with chromatin modification and post-transcriptional processes. In this study, we investigated the extent to which the functions of Paf1C combine to regulate the Saccharomyces cerevisiae transcriptome. While previous studies focused on the roles of Paf1C in controlling mRNA levels, here, we took advantage of a genetic background that enriches for unstable transcripts, and demonstrate that deletion of PAF1 affects all classes of Pol II transcripts including multiple classes of noncoding RNAs (ncRNAs). By conducting a de novo differential expression analysis independent of gene annotations, we found that Paf1 positively and negatively regulates antisense transcription at multiple loci. Comparisons with nascent transcript data revealed that many, but not all, changes in RNA levels detected by our analysis are due to changes in transcription instead of post-transcriptional events. To investigate the mechanisms by which Paf1 regulates protein-coding genes, we focused on genes involved in iron and phosphate homeostasis, which were differentially affected by PAF1 deletion. Our results indicate that Paf1 stimulates phosphate gene expression through a mechanism that is independent of any individual Paf1C-dependent histone modification. In contrast, the inhibition of iron gene expression by Paf1 correlates with a defect in H3 K36 trimethylation. Finally, we showed that one iron regulon gene, FET4, is coordinately controlled by Paf1 and transcription of upstream noncoding DNA. Together, these data identify roles for Paf1C in controlling both coding and noncoding regions of the yeast genome.

PLC and IP3-evoked Ca2+ release initiate the fast block to polyspermy in Xenopus laevis eggs.

The prevention of polyspermy is essential for the successful progression of normal embryonic development in most sexually reproducing species. In external fertilizers, the process of fertilization induces a depolarization of the egg's membrane within seconds, which inhibits supernumerary sperm from entering an already-fertilized egg. This fast block requires an increase of intracellular Ca2+ in the African clawed frog, Xenopus laevis, which in turn activates an efflux of Cl- that depolarizes the cell. Here we seek to identify the source of this intracellular Ca2+ Using electrophysiology, pharmacology, bioinformatics, and developmental biology, we explore the requirement for both Ca2+ entry into the egg from the extracellular milieu and Ca2+ release from an internal store, to mediate fertilization-induced depolarization. We report that although eggs express Ca2+-permeant ion channels, blockade of these channels does not alter the fast block. In contrast, insemination of eggs in the presence of Xestospongin C-a potent inhibitor of inositol 1,4,5-trisphosphate (IP3)-induced Ca2+ release from the endoplasmic reticulum (ER)-completely inhibits fertilization-evoked depolarization and increases the incidence of polyspermy. Inhibition of the IP3-generating enzyme phospholipase C (PLC) with U73122 similarly prevents fertilization-induced depolarization and increases polyspermy. Together, these results demonstrate that fast polyspermy block after fertilization in X. laevis eggs is mediated by activation of PLC, which increases IP3 and evokes Ca2+ release from the ER. This ER-derived Ca2+ then activates a Cl- channel to induce the fast polyspermy block. The PLC-induced cascade of events represents one of the earliest known signaling pathways initiated by fertilization.