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.

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.

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.

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.

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.

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.

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.

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.

The TMEM16A channel mediates the fast polyspermy block in Xenopus laevis.

In externally fertilizing animals, such as sea urchins and frogs, prolonged depolarization of the egg immediately after fertilization inhibits the entry of additional sperm-a phenomenon known as the fast block to polyspermy. In the African clawed frog Xenopus laevis, this depolarization is driven by Ca2+-activated Cl- efflux. Although the prominent Ca2+-activated Cl- currents generated in immature X. laevis oocytes are mediated by X. laevis transmembrane protein 16a (xTMEM16A) channels, little is known about the channels that contribute to the fast block in mature eggs. Moreover, the gamete undergoes a gross transformation as it develops from an immature oocyte into a fertilization-competent egg. Here, we report the results of our approach to identify the Ca2+-activated Cl- channel that triggers the fast block. By querying published proteomic and RNA-sequencing data, we identify two Ca2+-activated Cl- channels expressed in fertilization-competent X. laevis eggs: xTMEM16A and X. laevis bestrophin 2A (xBEST2A). By exogenously expressing xTMEM16A and xBEST2A in axolotl cells lacking endogenous Ca2+-activated currents, we characterize the effect of inhibitors on currents mediated by these channels. None of the inhibitors tested block xBEST2A currents specifically. However, 2-(4-chloro-2-methylphenoxy)-N-[(2-methoxyphenyl)methylideneamino]-acetamide (Ani9) and N-((4-methoxy)-2-naphthyl)-5-nitroanthranilic acid (MONNA) each reduce xTMEM16A currents by more than 70% while only nominally inhibiting those generated by xBEST2A. Using whole-cell recordings during fertilization, we find that Ani9 and MONNA effectively diminish fertilization-evoked depolarizations. Additionally, these inhibitors lead to increased polyspermy in X. laevis embryos. These results indicate that fertilization activates TMEM16A channels in X. laevis eggs and induces the earliest known event triggered by fertilization: the fast block to polyspermy.

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

This corrects the article DOI: 10.1038/nature18614.

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).

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.

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.

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.

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.

Subcellular RNA sequencing reveals broad presence of cytoplasmic intron-sequence retaining transcripts in mouse and rat neurons.

Recent findings have revealed the complexity of the transcriptional landscape in mammalian cells. One recently described class of novel transcripts are the Cytoplasmic Intron-sequence Retaining Transcripts (CIRTs), hypothesized to confer post-transcriptional regulatory function. For instance, the neuronal CIRT KCNMA1i16 contributes to the firing properties of hippocampal neurons. Intronic sub-sequence retention within IL1-ß mRNA in anucleate platelets has been implicated in activity-dependent splicing and translation. In a recent study, we showed CIRTs harbor functional SINE ID elements which are hypothesized to mediate dendritic localization in neurons. Based on these studies and others, we hypothesized that CIRTs may be present in a broad set of transcripts and comprise novel signals for post-transcriptional regulation. We carried out a transcriptome-wide survey of CIRTs by sequencing micro-dissected subcellular RNA fractions. We sequenced two batches of 150-300 individually dissected dendrites from primary cultures of hippocampal neurons in rat and three batches from mouse hippocampal neurons. After statistical processing to minimize artifacts, we found a broad prevalence of CIRTs in the neurons in both species (44-60% of the expressed transcripts). The sequence patterns, including stereotypical length, biased inclusion of specific introns, and intron-intron junctions, suggested CIRT-specific nuclear processing. Our analysis also suggested that these cytoplasmic intron-sequence retaining transcripts may serve as a primary transcript for ncRNAs. Our results show that retaining intronic sequences is not isolated to a few loci but may be a genome-wide phenomenon for embedding functional signals within certain mRNA. The results hypothesize a novel source of cis-sequences for post-transcriptional regulation. Our results hypothesize two potentially novel splicing pathways: one, within the nucleus for CIRT biogenesis; and another, within the cytoplasm for removing CIRT sequences before translation. We also speculate that release of CIRT sequences prior to translation may form RNA-based signals within the cell potentially comprising a novel class of signaling pathways.

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.

Ribosome profiling shows that miR-430 reduces translation before causing mRNA decay in zebrafish.

MicroRNAs regulate gene expression through deadenylation, repression, and messenger RNA (mRNA) decay. However, the contribution of each mechanism in non-steady-state situations remains unclear. We monitored the impact of miR-430 on ribosome occupancy of endogenous mRNAs in wild-type and dicer mutant zebrafish embryos and found that miR-430 reduces the number of ribosomes on target mRNAs before causing mRNA decay. Translational repression occurs before complete deadenylation, and disrupting deadenylation with use of an internal polyadenylate tail did not block target repression. Lastly, we observed that ribosome density along the length of the message remains constant, suggesting that translational repression occurs by reducing the rate of initiation rather than affecting elongation or causing ribosomal drop-off. These results show that miR-430 regulates translation initiation before inducing mRNA decay during zebrafish development.