Search

The Balbiani body (Bb) is the first marker of polarity in vertebrate oocytes. The Bb is a conserved structure found in diverse animals including insects, fish, amphibians, and mammals. During early zebrafish oogenesis, the Bb assembles as a transient aggregate of mRNA, proteins, and membrane-bound organelles at the presumptive vegetal side of the oocyte. As the early oocyte develops, the Bb appears to grow slowly, until at the end of stage I of oogenesis it disassembles and deposits its cargo of localized mRNAs and proteins. In fish and frogs, this cargo includes the germ plasm as well as gene products required to specify dorsal tissues of the future embryo. We demonstrate that the Bb is a stable, solid structure that forms a size exclusion barrier similar to other biological hydrogels. Despite its central role in oocyte polarity, little is known about the mechanism behind the Bb's action. Analysis of the few known protein components of the Bb is insufficient to explain how the Bb assembles, translocates, and disassembles. We isolated Bbs from zebrafish oocytes and performed mass spectrometry to define the Bb proteome. We successfully identified 77 proteins associated with the Bb sample, including known Bb proteins and novel RNA-binding proteins. In particular, we identified Cirbpa and Cirbpb, which have both an RNA-binding domain and a predicted self-aggregation domain. In stage I oocytes, Cirbpa and Cirbpb localize to the Bb rather than the nucleus (as in somatic cells), indicating that they may have a specialized function in the germ line. Both the RNA-binding domain and the self-aggregation domain are sufficient to localize to the Bb, suggesting that Cirbpa and Cirbpb interact with more than just their mRNA targets within the Bb. We propose that Cirbp proteins crosslink mRNA cargo and proteinaceous components of the Bb as it grows. Beyond Cirbpa and Cirbpb, our proteomics dataset presents many candidates for further study, making it a valuable resource for building a comprehensive mechanism for Bb function at a protein level.

Thousands of long intergenic non-coding RNAs (lincRNAs) are transcribed throughout the vertebrate genome. A subset of lincRNAs enriched in developing brains has recently been found to contain cryptic open reading frames and are speculated to encode micropeptides. However, systematic identification and functional assessment of these transcripts have been hindered by technical challenges caused by their small size. Here we show that two putative lincRNAs (linc-mipep and linc-wrb) encode micropeptides with homology to the vertebrate-specific chromatin architectural protein, Hmgn1, and demonstrate that they are required for development of vertebrate-specific brain cell types. Specifically, we show that NMDA receptor-mediated pathways are dysregulated in zebrafish lacking these micropeptides and that their loss preferentially alters the gene regulatory networks that establish cerebellar cells and oligodendrocytes – evolutionarily newer cell types that develop postnatally in humans. These findings highlight the power of screening for unexplored micropeptide functions by revealing a key missing link in the evolution of vertebrate brain cell development and illustrating a genetic basis for how some neural cell types are more susceptible to chromatin disruptions, with implications for neurodevelopmental disorders and disease.

Although gene expression is tightly regulated during embryonic development, the impact of translational control has received less experimental attention. Here, we find that eukaryotic translation initiation factor-3 (eIF3) is required for Shh-mediated tissue patterning. Analysis of loss-of-function eIF3 subunit c (Eif3c) mice reveal a unique sensitivity to the Shh receptor Patched 1 (Ptch1) dosage. Genome-wide in vivo enhanced cross-linking immunoprecipitation sequence (eCLIP-seq) shows unexpected specificity for eIF3 binding to a pyrimidine rich motif present in subsets of 5’-UTRs and a corresponding change in the translation of these transcripts by ribosome profiling in Eif3c loss-of-function embryos. We further find that while Eif3c loss-of-function embryos do not show a global decrease in protein synthesis, translation of Ptch1 through this pyrimidine rich motif is specifically sensitive to eIF3 amount. Altogether, this work uncovers hidden specificity of housekeeping translation initiation machinery for the translation of key developmental signaling transcripts.

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3' end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3' end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome. We thank all the members of the Novoa lab for their valuable insights and discussion. O.B. is supported by a UNSW International PhD fellowship and Australian Government Research Training Program Scholarship. G.D. is part of the ROPES ITN, which received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement number 956810. A.D.-T. is supported by an FPI-SO fellowship (PRE2019-088498). This work was supported by the Australian Research Council (DP180103571 to E.M.N.), the Spanish Ministry of Economy, Industry and Competitiveness (MEIC) (PGC2018-098152-A-100 to E.M.N.), the European Research Council (ERC-StG-2021 No 101042103 to E.M.N.) and the NIH (R35GM122580 and RO1HD00035 to A.J.G). We acknowledge the support of the MEIC to the EMBL partnership, Centro de Excelencia Severo Ochoa and CERCA Programme/Generalitat de Catalunya. We thank T. Mercer (Garvan Institute, Australia) for providing us with the sequins used as spike-ins in Nano3P-seq zebrafish runs. We thank the Sdelci lab (CRG, Spain) for providing us with HeLa cell lines pellets used in this work.

Oogenesis produces functional eggs and is essential for fertility, embryonic development, and reproduction. The zebrafish ovary is an excellent model to study oogenesis in vertebrates, and recent studies have identified multiple regulators in oocyte development through forward genetic screens, as well as reverse genetics by CRISPR mutagenesis. However, many developmental steps in oogenesis, in zebrafish and other species, remain poorly understood, and their underlying mechanisms are unknown. Here, we take a genomic approach to systematically uncover biological activities throughout oogenesis. We performed transcriptomic analysis on five stages of oogenesis, from the onset of oocyte differentiation through Stage III, which precedes oocyte maturation. These transcriptomes revealed thousands of differentially expressed genes across stages of oogenesis. We analyzed trends of gene expression dynamics along oogenesis, as well as their expression in pair-wise comparisons between stages. We determined their functionally enriched terms, identifying uniquely characteristic biological activities in each stage. These data identified two prominent developmental phases in oocyte differentiation and traced the accumulation of maternally deposited embryonic regulator transcripts in the developing oocyte. Our analysis provides the first molecular description for oogenesis in zebrafish, which we deposit online as a resource for the community. Further, the presence of multiple gene paralogs in zebrafish, and the exclusive curation by many bioinformatic tools of the single paralogs present in humans, challenge zebrafish genomic analyses. We offer an approach for converting zebrafish gene name nomenclature to the human nomenclature for supporting genomic analyses generally in zebrafish. Altogether, our work provides a valuable resource as a first step to uncover oogenesis mechanisms and candidate regulators and track accumulating transcripts of maternal regulators of embryonic development.

Messenger RNA (mRNA) compartmentalization within the cytosol is well-recognized as a key mechanism of local translation-mediated regulation of protein levels, but whether such localization could be a means of exercising non-coding mRNA function is unknown. Here, we explore non-coding functions for mRNAs associated with focal adhesions (FAs), cellular structures responsible for mediating cell adhesion and response to changes in the extracellular matrix (ECM). Using high-throughput single molecule imaging and genomic profiling approaches, we find that mRNAs with distinct sequence characteristics localize to FAs in different human cell types. Notably, ∼85% of FA-mRNAs are not translationally active at steady state or under conditions of FA dissolution or activation. Untranslated mRNA sequences are anchored to FA based on their functional states by the RNA binding protein, G3BP1, forming biomolecular granules. Removing RNA or G3BP1, but not blocking new polypeptide synthesis, dramatically changes FA protein composition and organization, resulting in loss of cell contractility and cellular ability to adapt to changing ECM. We have therefor uncovered a novel, non-coding role for mRNAs as scaffolds to maintain FA structure and function, broadening our understating of noncanonical mRNA functions.

Zygotic genome activation (ZGA) initiates regionalized transcription underlying distinct cellular identities. ZGA is dependent upon dynamic chromatin architecture sculpted by conserved DNA-binding proteins. However, the direct mechanistic link between the onset of ZGA and the tissue-specific transcription remains unclear. Here, we have addressed the involvement of chromatin organizer Satb2 in orchestrating both processes during zebrafish embryogenesis. Integrative analysis of transcriptome, genome-wide occupancy and chromatin accessibility reveals contrasting molecular activities of maternally deposited and zygotically synthesized Satb2. Maternal Satb2 prevents premature transcription of zygotic genes by influencing the interplay between the pluripotency factors. By contrast, zygotic Satb2 activates transcription of the same group of genes during neural crest development and organogenesis. Thus, our comparative analysis of maternal versus zygotic function of Satb2 underscores how these antithetical activities are temporally coordinated and functionally implemented highlighting the evolutionary implications of the biphasic and bimodal regulation of landmark developmental transitions by a single determinant., Activation of the zygotic genome is a critical transition during development, though the link to tissue-specific gene regulation remains unclear. Here the authors demonstrate distinct functions for Satb2 before and after zygotic genome activation, highlighting the temporal coordination of these roles.

Rotating cilia at the vertebrate left-right organizer (LRO) generate an asymmetric leftward flow, which is sensed by cells at the left LRO margin. Ciliary activity of the calcium channel Pkd2 is crucial for flow sensing. How this flow signal is further processed and relayed to the laterality-determining Nodal cascade in the left lateral plate mesoderm (LPM) is largely unknown. We previously showed that flow down-regulates mRNA expression of the Nodal inhibitor Dand5 in left sensory cells. De-repression of the co-expressed Nodal, complexed with the TGFß growth factor Gdf3, drives LPM Nodal cascade induction. Here, we show that post-transcriptional repression of dand5 is a central process in symmetry breaking of Xenopus, zebrafish and mouse. The RNA binding protein Bicc1 was identified as a post-transcriptional regulator of dand5 and gdf3 via their 3′-UTRs. Two distinct Bicc1 functions on dand5 mRNA were observed at pre- and post-flow stages, affecting mRNA stability or flow induced translational inhibition, respectively. To repress dand5, Bicc1 co-operates with Dicer1, placing both proteins in the process of flow sensing. Intriguingly, Bicc1 mediated translational repression of a dand5 3′-UTR mRNA reporter was responsive to pkd2, suggesting that a flow induced Pkd2 signal triggers Bicc1 mediated dand5 inhibition during symmetry breakage., The authors show that post-transcriptional regulation of the cilia-driven leftward flow target dand5 is central to symmetry breakage in frog, fish and mouse and is mediated by a 139 nt Bicc1 responsive element in the dand5 3′UTR, and they present evidence that Pkd2 regulates this Bicc1/dand5 module.

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.

Genetic robustness, or the ability of an organism to maintain fitness in the presence of harmful mutations, can be achieved via protein feedback loops. Previous work has suggested that organisms may also respond to mutations by transcriptional adaptation, a process by which related gene(s) are upregulated independently of protein feedback loops. However, the prevalence of transcriptional adaptation and its underlying molecular mechanisms are unknown. Here, by analysing several models of transcriptional adaptation in zebrafish and mouse, we uncover a requirement for mutant mRNA degradation. Alleles that fail to transcribe the mutated gene do not exhibit transcriptional adaptation, and these alleles give rise to more severe phenotypes than alleles displaying mutant mRNA decay. Transcriptome analysis in alleles displaying mutant mRNA decay reveals the upregulation of a substantial proportion of the genes that exhibit sequence similarity with the mutated gene's mRNA, suggesting a sequence-dependent mechanism. These findings have implications for our understanding of disease-causing mutations, and will help in the design of mutant alleles with minimal transcriptional adaptation-derived compensation.

The requirement for Cas nucleases to recognize a specific PAM is a major restriction for genome editing. SpCas9 variants SpG and SpRY, recognizing NGN and NRN PAM, respectively, have contributed to increase the number of editable genomic sites in cell cultures and plants. However, their use has not been demonstrated in animals.We have characterized and optimized the activity of SpG and SpRY in zebrafish and C. elegans. Delivered as mRNA-gRNA or ribonucleoprotein (RNP) complexes, SpG and SpRY were able to induce mutations in vivo, albeit at a lower rate than SpCas9 in equivalent formulations. This lower activity was overcome by optimizing mRNA-gRNA or RNP concentration, leading to efficient mutagenesis at regions inaccessible to SpCas9. We also found that the CRISPRscan algorithm can predict SpG and SpRY activity in vivo. Finally, we applied SpG and SpRY to generate knock-ins by homology-directed repair. Altogether, our results expand the CRISPR-Cas targeting genomic landscape in animals.

Improved molecular screening and diagnostic tools are needed to substantially increase SARS-CoV-2 testing capacity and throughput while reducing the time to receive test results. Here we developed multiplex reverse transcriptase polymerase chain reaction (m-RT-PCR) for detection of SARS-CoV-2 using rapid DNA electrophoresis and alternatively using multiplex viral sequencing (mVseq). For RNA specimens extracted from nasopharyngeal (NP) swabs in viral transport media (VTM), our assays achieved a sensitivity for SARS-CoV-2 detection corresponding to cycle threshold (Ct) of 37.2 based on testing of these specimens using quantitative reverse transcription PCR (RT-qPCR). For NP swab-VTM specimens without prior RNA extraction, sensitivity was reduced to Ct of 31.6, which was due to lower concentration of SARS-CoV-2 genome copies in VTM compared to RNA-extracted samples. Assay turnaround time was 60 minutes using rapid gel electrophoresis, 90 minutes using Agilent Bioanalyzer, and 24-48 hours using Illumina sequencing, the latter of which required a second PCR to produce a sequence-ready library using m-RT-PCR products as the template. Our assays can be employed for high-throughput sequencing-based detection of SARS-CoV-2 directly from a clinical specimen without RNA isolation, while ease-of-use and low cost of the electrophoresis-based readout enables screening, particularly in resource-constrained settings.

Genome-wide chromatin reprogramming is a fundamental requirement for establishing developmental competence in the newly-formed zygote. In zebrafish, Nanog, Pou5f3 and Sox19b play partially redundant roles in zygotic genome activation, however their interplay in establishing chromatin competency, the context in which they do so and their mechanism of action remain poorly defined. Here, we generated a triple maternal-zygoticnanog-/-;pou5f3-/-;sox19b-/-mutant and assessed the causal relationship between transcription factor (TF) occupancy, chromatin accessibility and genome activation. Analyses of this triple mutant and combinatorial rescues revealed highly synergistic and context-dependent activity of Nanog, Pou5f3, and Sox19b (NPS) in establishing chromatin competency at >50% of active enhancers. Motif analysis revealed a network of TFs that depend on NPS for establishing chromatin accessibility, including the endodermal determinant Eomesa, whose binding we show is regulated by NPS pioneer-like activity. Finally, we demonstrated that NPS play an essential role in establishing H3K27ac and H3K18ac at enhancers and promoters, and that their function in transcriptional activation can be bypassed by targeted recruitment of histone acetyltransferases to individual genes. Altogether, our findings reveal a large network of TFs that function to establish developmental competency, many of which depend on the synergistic and highly context-dependent role of NPS in establishing chromatin accessibility and regulating histone acetylation in order to activate the genome.

During the maternal-to-zygotic transition (MZT), multiple mechanisms precisely control massive decay of maternal mRNAs. N

The xbp-1 mRNA encodes the XBP-1 transcription factor, a critical part of the unfolded protein response. Here we report that an RNA fragment produced from xbp-1 mRNA cleavage is a biologically active non-coding RNA (ncRNA) in Caenorhabditis elegans neurons, providing the first example of ncRNA derived from mRNA cleavage. We show that the xbp-1 ncRNA is crucial for axon regeneration in vivo, and that it acts independently of the protein-coding function of the xbp-1 transcript. Structural analysis indicates that the function of the xbp-1 ncRNA depends on a single RNA stem; and this stem forms only in the cleaved xbp-1 ncRNA fragment. Disruption of this stem abolishes the non-coding but not coding function of the endogenous xbp-1 transcript. Thus, cleavage of the xbp-1 mRNA bifurcates it into a coding and a non-coding pathway; modulation of the two pathways may allow neurons to fine-tune their response to injury and other stresses.Graphic abstract

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.

RNA folding plays a crucial role in RNA function. However, knowledge of the global structure of the transcriptome is limited to cellular systems at steady state, thus hindering the understanding of RNA structure dynamics during biological transitions and how it influences gene function. Here, we characterized mRNA structure dynamics during zebrafish development. We observed that on a global level, translation guides structure rather than structure guiding translation. We detected a decrease in structure in translated regions and identified the ribosome as a major remodeler of RNA structure in vivo. In contrast, we found that 3′ untranslated regions (UTRs) form highly folded structures in vivo, which can affect gene expression by modulating microRNA activity. Furthermore, dynamic 3′-UTR structures contain RNA-decay elements, such as the regulatory elements in nanog and ccna1, two genes encoding key maternal factors orchestrating the maternal-to-zygotic transition. These results reveal a central role of RNA structure dynamics in gene regulatory programs. Characterization of mRNA structure during the zebrafish maternal-to-zygotic transition identifies the ribosome as a major RNA structure remodeler in vivo and reveals that structural dynamics can affect gene expression, partly by modulating miRNA activity.

Cpf1 is a novel class of CRISPR-Cas DNA endonucleases, with a wide range of activity across different eukaryotic systems. Yet, the underlying determinants of this variability are poorly understood. Here, we demonstrate that LbCpf1, but not AsCpf1, ribonucleoprotein complexes allow efficient mutagenesis in zebrafish and Xenopus. We show that temperature modulates Cpf1 activity by controlling its ability to access genomic DNA. This effect is stronger on AsCpf1, explaining its lower efficiency in ectothermic organisms. We capitalize on this property to show that temporal control of the temperature allows post-translational modulation of Cpf1-mediated genome editing. Finally, we determine that LbCpf1 significantly increases homology-directed repair in zebrafish, improving current approaches for targeted DNA integration in the genome. Together, we provide a molecular understanding of Cpf1 activity in vivo and establish Cpf1 as an efficient and inducible genome engineering tool across ectothermic species., Cpf1 is a promising addition to the CRISPR toolkit but displays wide variability of activity in different eurkaryotes. Here the authors identify temperature as a modulator of activity and use this to efficiently edit ectothermic vertebrate species.

Rotating cilia at the vertebrate left-right organizer (LRO) generate an asymmetric leftward flow, which is sensed by cells at the left LRO margin. How the flow signal is processed and relayed to the laterality-determining Nodal cascade in the left lateral plate mesoderm (LPM) is largely unknown. We previously showed that flow down-regulates mRNA expression of the Nodal inhibitor Dand5 in left sensory cells. De-repression of the co-expressed Nodal drives LPM Nodal cascade induction. Here, we identify the mechanism of dand5 downregulation, finding that its posttranscriptional repression is a central process in symmetry breaking. Specifically, the RNA binding protein Bicc1 interacts with a proximal element in the 3’-UTR of dand5 to repress translation in a dicer1-dependent manner. The bicc1/dicer1 module acts downstream of flow, as LRO ciliation was not affected upon its loss. Loss of bicc1 or dicer1 was rescued by parallel knockdown of dand5, placing both genes in the process of flow sensing.

Zebrafish are able to regenerate heart muscle after catastrophic injury. Fervent research is underway to discover critical molecules with therapeutic potential, yet many principal biological processes that underlie regeneration likely remain undiscovered. Changes in the abundance of dozens to thousands of mRNA have been reported using both targeted approaches like in-situ hybridization and unbiased approaches like RNA-seq. However, whether the abundance of mRNA represents the entirety of gene expression control remains unknown. Regulation of translation may alter the abundance of a critical protein product without requiring an intervening transcription event. Here, we have taken advantage of ribosome profiling to document the changes in translation efficiency that accompany heart regeneration. Using high-throughput sequencing of the ~28 nucleotide fragments of mRNA protected from enzymatic degradation, we find thousands of transcripts that change in their association with ribosomes. We also find that key components of the translational machinery are transcriptionally induced during regeneration, supporting a general model in which regulation of translation is a fundamental component of zebrafish heart regeneration.

The awakening of the genome after fertilization is a cornerstone of animal development. However, the mechanisms that activate the silent genome after fertilization are poorly understood. Here, we show that transcriptional competency is regulated by Brd4- and P300-dependent histone acetylation in zebrafish. Live imaging of transcription revealed that genome activation, beginning at the miR-430 locus, is gradual and stochastic. We show that genome activation does not require slowdown of the cell cycle and is regulated through the translation of maternally inherited mRNAs. Among these, the enhancer regulators P300 and Brd4 can prematurely activate transcription and restore transcriptional competency when maternal mRNA translation is blocked, whereas inhibition of histone acetylation blocks genome activation. We conclude that P300 and Brd4 are sufficient to trigger genome-wide transcriptional competency by regulating histone acetylation on the first zygotic genes in zebrafish. This mechanism is critical for initiating zygotic development and developmental reprogramming., This research was supported by the Fonds de Recherche du Québec Santé (postdoctoral fellowship to J.-D.B.), the Surdna Foundation and the Yale Genetics Venture Fund (postdoctoral fellowship to L.M.), and the NIH (grants R01 HD074078, GM103789, GM102251, GM101108, GM081602, R35 GM122580, and 4DNucleome program to A.J.G.). M.A.M.-M. was supported by Programa de Movilidad en Areas de Investigación priorizadas por la Consejería de Igualdad, Salud y Políticas Sociales de la Junta de Andalucía and is currently supported by Ramon y Cajal program from Ministerio de Ciencia, Innovación y Universidades of the Spanish government (RYC-2017-23041). The Giraldez lab is supported by the Howard Hughes Medical Institute Faculty Scholar program and was supported by the Pew Scholars Program in the Biomedical Sciences, March of Dimes 1-FY12-230, the Yale Scholars Program, and Whitman fellowship funds provided by E.E. Just, Lucy B. Lemann, Evelyn and Melvin Spiegel, The Haldan Keffer Hartline and Edward F. MacNichol, Jr., of the Marine Biological Laboratory in Woods Hole, Massachusetts to A.J.G.

Cellular transitions require dramatic changes in gene expression that are supported by regulated mRNA decay and new transcription. The maternal‐to‐zygotic transition is a conserved developmental progression during which thousands of maternal mRNAs are cleared by post‐transcriptional mechanisms. Although some maternal mRNAs are targeted for degradation by microRNAs, this pathway does not fully explain mRNA clearance. We investigated how codon identity and translation affect mRNA stability during development and homeostasis. We show that the codon triplet contains translation‐dependent regulatory information that influences transcript decay. Codon composition shapes maternal mRNA clearance during the maternal‐to‐zygotic transition in zebrafish, Xenopus , mouse, and Drosophila , and gene expression during homeostasis across human tissues. Some synonymous codons show consistent stabilizing or destabilizing effects, suggesting that amino acid composition influences mRNA stability. Codon composition affects both polyadenylation status and translation efficiency. Thus, the ribosome interprets two codes within the mRNA: the genetic code which specifies the amino acid sequence and a conserved “codon optimality code” that shapes mRNA stability and translation efficiency across vertebrates.

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.

During the maternal-to-zygotic transition (MZT), multiple mechanisms precisely control massive decay of maternal mRNAs. N6-methyladenosine (m6A) is known to regulate mRNA decay, yet how this modification promotes maternal transcript degradation remains unclear. Here, we find that m6A promotes maternal mRNA deadenylation. Yet, genetic loss of m6A readers Ythdf2 and Ythdf3 did not impact global maternal mRNA clearance, zygotic genome activation, or the onset of gastrulation, challenging the view that Ythdf2 alone is critical to developmental timing. We reveal that Ythdf proteins function redundantly during zebrafish oogenesis and development, as double Ythdf2 and Ythdf3 deletion prevented female gonad formation and triple Ythdf mutants were lethal. Finally, we show that the microRNA miR-430 functions additively with methylation to promote degradation of common transcript targets. Together these findings reveal that m6A facilitates maternal mRNA deadenylation and that multiple pathways and readers act in concert to mediate these effects of methylation on RNA stability.

Current methods for comparing scRNA-seq datasets collected in multiple conditions focus on discrete regions of the transcriptional state space, such as clusters of cells. Here, we quantify the effects of perturbations at the single-cell level using a continuous measure of the effect of a perturbation across the transcriptomic space. We describe this space as a manifold and develop a relative likelihood estimate of observing each cell in each of the experimental conditions using graph signal processing. This likelihood estimate can be used to identify cell populations specifically affected by a perturbation. We also develop vertex frequency clustering to extract populations of affected cells at the level of granularity that matches the perturbation response. The accuracy of our algorithm to identify clusters of cells that are enriched or depleted in each condition is on average 57% higher than the next best-performing algorithm tested. Gene signatures derived from these clusters are more accurate compared to six alternative algorithms in ground-truth comparisons.

The awakening of the zygote genome, signaling the transition from maternal transcriptional control to zygotic control, is a watershed in embryonic development, but the factors and mechanisms controlling this transition are still poorly understood. By combining CRISPR-Cas9-mediated live imaging of the first transcribed genes (miR-430), chromatin and transcription analysis during zebrafish embryogenesis, we observed that genome activation is gradual and stochastic, and the active state is inherited in daughter cells. We discovered that genome activation is regulated through both translation of maternal mRNAs and the effects of these factors on the chromatin. We show that chemical inhibition of H3K27Ac writer (P300) and reader (Brd4) block genome activation, while induction of a histone acetylation prematurely activates transcription, and restore genome activation in embryos where translation of maternal mRNAs is impaired, demonstrating that they are limiting factors for the activation of the genome. In contrast to current models, we do not observe triggering of genome activation by a reduction of the nuclear-cytoplasmic (N/C) ratio or slower cell division. We conclude that genome activation is controlled by a time-dependent mechanism involving the translation of maternal mRNAs and the regulation of histone acetylation through P300 and Brd4. This mechanism is critical to initiating zygotic development and developmental reprogramming.

Genetic compensation by transcriptional modulation of related gene(s) (also known as transcriptional adaptation) has been reported in numerous systems1–3; however, whether and how such a response can be activated in the absence of protein feedback loops is unknown. Here, we develop and analyze several models of transcriptional adaptation in zebrafish and mouse that we show are not caused by loss of protein function. We find that the increase in transcript levels is due to enhanced transcription, and observe a correlation between the levels of mutant mRNA decay and transcriptional upregulation of related genes. To assess the role of mutant mRNA degradation in triggering transcriptional adaptation, we use genetic and pharmacological approaches and find that mRNA degradation is indeed required for this process. Notably, uncapped RNAs, themselves subjected to rapid degradation, can also induce transcriptional adaptation. Next, we generate alleles that fail to transcribe the mutated gene and find that they do not show transcriptional adaptation, and exhibit more severe phenotypes than those observed in alleles displaying mutant mRNA decay. Transcriptome analysis of these different alleles reveals the upregulation of hundreds of genes with enrichment for those showing sequence similarity with the mutated gene’s mRNA, suggesting a model whereby mRNA degradation products induce the response via sequence similarity. These results expand the role of the mRNA surveillance machinery in buffering against mutations by triggering the transcriptional upregulation of related genes. Besides implications for our understanding of disease-causing mutations, our findings will help design mutant alleles with minimal transcriptional adaptation-derived compensation.

Post-transcriptional regulation is crucial to shape gene expression. During the Maternal-to-Zygotic Transition (MZT), thousands of maternal transcripts are regulated upon fertilization and genome activation. Transcript stability can be influenced bycis-elements andtrans-factors, but how these inputs are integrated to determine the overall mRNA stability is unclear. Here, we show that most transcripts are under combinatorial regulation by multiple decay pathways during zebrafish MZT. To identifycis-regulatory elements, we performed a massively parallel reporter assay for stability-influencing sequences, which revealed that 3’-UTR poly-U motifs are associated with mRNA stability. In contrast, miR-430 target sequences, UAUUUAUU AU-rich elements (ARE), CCUC and CUGC elements emerged as the main destabilizing motifs in the embryo, with miR-430 and AREs causing mRNA deadenylation in a genome activation-dependent manner. To identify thetrans-factors interacting with thesecis-elements, we comprehensively profiled RNA-protein interactions and their associated regulatory activities across the transcriptome during the MZT. We find that poly-U binding proteins are preferentially associated with 3’-UTR sequences and stabilizing motifs, and that antagonistic sequence contexts for poly-C and poly-U binding proteins shape the binding landscape and magnitude of regulation across the transcriptome. Finally, we integrate these regulatory motifs into a machine learning model that accurately predicts the stability of mRNA reportersin vivo. Our findings reveal how mechanisms of post-transcriptional regulation are coordinated to direct changes in mRNA stability within the early zebrafish embryo.

RNA folding plays a crucial role in RNA function. However, our knowledge of the global structure of the transcriptome is limited to steady-state conditions, hindering our understanding of how RNA structure dynamics influences gene function. Here, we have characterized mRNA structure dynamics during zebrafish development. We observe that on a global level, translation guides structure rather than structure guiding translation. We detect a decrease in structure in translated regions, and we identify the ribosome as a major remodeler of RNA structure in vivo. In contrast, we find that 3’-UTRs form highly folded structures in vivo, which can affect gene expression by modulating miRNA activity. Furthermore, we find that dynamic 3’-UTR structures encode RNA decay elements, including regulatory elements in nanog and cyclin A1, key maternal factors orchestrating the maternal-to-zygotic transition. These results reveal a central role of RNA structure dynamics in gene regulatory programs.

Regulation of gene expression is fundamental in establishing cellular diversity and a target of natural selection. Untranslated mRNA regions (UTRs) are key mediators of post-transcriptional regulation. Previous studies have predicted thousands of ORFs in 5'UTRs, the vast majority of which have unknown function. Here, we present a systematic analysis of the translation and function of upstream open reading frames (uORFs) across vertebrates. Using high-resolution ribosome footprinting, we find that (i)uORFs are prevalent within vertebrate transcriptomes, (ii) the majority show signatures of active translation, and (iii)uORFs act as potent regulators of translation and RNA levels, with a similar magnitude to miRNAs. Reporter experiments reveal clear repression of downstream translation by uORFs/oORFs. uORF number, intercistronic distance, overlap with the CDS, and initiation context most strongly influence translation. Evolution has targeted these features to favor uORFs amenable to regulation over constitutively repressive uORFs/oORFs. Finally, we observe that the regulatory potential of uORFs on individual genes is conserved across species. These results provide insight into the regulatory code within mRNA leader sequences and their capacity to modulate translation across vertebrates.

Autism spectrum disorders (ASDs) are a group of devastating neurodevelopmental syndromes that affect up to 1 in 68 children. Despite advances in the identification of ASD risk genes, the mechanisms underlying ASDs remain unknown. Homozygous loss-of-function mutations in Contactin Associated Protein-like 2 (CNTNAP2) are strongly linked to ASDs. Here we investigate the function of Cntnap2 and undertake pharmacological screens to identify phenotypic suppressors. We find that zebrafish cntnap2 mutants display GABAergic deficits, particularly in the forebrain, and sensitivity to drug-induced seizures. High-throughput behavioral profiling identifies nighttime hyperactivity in cntnap2 mutants, while pharmacological testing reveals dysregulation of GABAergic and glutamatergic systems. Finally, we find that estrogen receptor agonists elicit a behavioral fingerprint anti-correlative to that of cntnap2 mutants and show that the phytoestrogen biochanin A specifically reverses the mutant behavioral phenotype. These results identify estrogenic compounds as phenotypic suppressors and illuminate novel pharmacological pathways with relevance to autism.

MicroRNAs have emerged as critical regulators of gene expression. Originally shown to regulate developmental timing, microRNAs have since been implicated in a wide range of cellular functions including cell identity, migration and signaling. miRNA-430, the earliest expressed microRNA during zebrafish embryogenesis, is required to undergo morphogenesis and has previously been shown to regulate maternal mRNA clearance, Nodal signaling, and germ cell migration. The functions of miR-430 in brain morphogenesis, however, remain unclear. Herein we find that miR-430 instructs oriented cell divisons in the neural rod required for neural midline formation. Loss of miR-430 function results in mitotic spindle misorientation in the neural rod, failed neuroepithelial integration after cell division, and ectopic cell accumulation in the dorsal neural tube. We propose that miR-430, independently of canonical apicobasal and planar cell polarity (PCP) pathways, coordinates the stereotypical cell divisons that instruct neural tube morphogenesis.

CRISPR-Cas9 technology provides a powerful system for genome engineering. However, variable activity across different single guide RNAs (sgRNAs) remains a significant limitation. We analyzed the molecular features that influence sgRNA stability, activity and loading into Cas9 in vivo. We observed that guanine enrichment and adenine depletion increased sgRNA stability and activity, whereas differential sgRNA loading, nucleosome positioning and Cas9 off-target binding were not major determinants. We also identified sgRNAs truncated by one or two nucleotides and containing 5' mismatches as efficient alternatives to canonical sgRNAs. On the basis of these results, we created a predictive sgRNA-scoring algorithm, CRISPRscan, that effectively captures the sequence features affecting the activity of CRISPR-Cas9 in vivo. Finally, we show that targeting Cas9 to the germ line using a Cas9-nanos 3' UTR led to the generation of maternal-zygotic mutants, as well as increased viability and decreased somatic mutations. These results identify determinants that influence Cas9 activity and provide a framework for the design of highly efficient sgRNAs for genome targeting in vivo.

As the range and duration of human ventures into space increase, it becomes imperative that we understand the effects of the cosmic environment on astronaut health. Molecular technologies now widely used in research and medicine will need to become available in space to ensure appropriate care of astronauts. The polymerase chain reaction (PCR) is the gold standard for DNA analysis, yet its potential for use on-orbit remains under-explored. We describe DNA amplification aboard the International Space Station (ISS) through the use of a miniaturized miniPCR system. Target sequences in plasmid, zebrafish genomic DNA, and bisulfite-treated DNA were successfully amplified under a variety of conditions. Methylation-specific primers differentially amplified bisulfite-treated samples as would be expected under standard laboratory conditions. Our findings establish proof of concept for targeted detection of DNA sequences during spaceflight and lay a foundation for future uses ranging from environmental monitoring to on-orbit diagnostics.

Pre-mRNA splicing is a critical step of gene expression in eukaryotes. Transcriptome-wide splicing patterns are complex and primarily regulated by a diverse set of recognition elements and associated RNA-binding proteins. The retention and splicing (RES) complex is formed by three different proteins (Bud13p, Pml1p and Snu17p) and is involved in splicing in yeast. However, the importance of the RES complex for vertebrate splicing, the intronic features associated with its activity, and its role in development are unknown. In this study, we have generated loss-of-function mutants for the three components of the RES complex in zebrafish and showed that they are required during early development. The mutants showed a marked neural phenotype with increased cell death in the brain and a decrease in differentiated neurons. Transcriptomic analysis ofbud13,snip1(pml1) andrbmx2(snu17) mutants revealed a global defect in intron splicing, with strong mis-splicing of a subset of introns. We found these RES-dependent introns were short, rich in GC and flanked by GC depleted exons, all of which are features associated with intron definition. Using these features we developed a predictive model that classifies RES dependent introns. Altogether, our study uncovers the essential role of the RES complex during vertebrate development and provides new insights into its function during splicing.

Genome-wide RNA structure maps have recently become available through the coupling of in vivo chemical probing reagents with next-generation sequencing. Initial analyses relied on the identification of truncated reverse transcription reads to identify the chemically modified nucleotides, but recent studies have shown that mutational signatures can also be used. While these two methods have been employed interchangeably, here we show that they actually provide complementary information. Consequently, analyses using exclusively one of the two methodologies may disregard a significant portion of the structural information. We also show that the identity and sequence environment of the modified nucleotide greatly affect the odds of introducing a mismatch or causing reverse transcriptase drop-off. Finally, we identify specific mismatch signatures generated by dimethyl sulfate probing that can be exploited to remove false positives typically produced in RNA structurome analyses, and how these signatures vary depending on the reverse transcription enzyme used.

Proper functioning of an organism requires cells and tissues to behave in uniform, well-organized ways. How this optimum of phenotypes is achieved during the development of vertebrates is unclear. Here, we carried out a multi-faceted and single-cell resolution screen of zebrafish embryonic blood vessels upon mutagenesis of single and multi-gene microRNA (miRNA) families. We found that embryos lacking particular miRNA-dependent signaling pathways develop a vascular trait similar to wild-type, but with a profound increase in phenotypic heterogeneity. Aberrant trait variance in miRNA mutant embryos uniquely sensitizes their vascular system to environmental perturbations. We discovered a previously unrecognized role for specific vertebrate miRNAs to protect tissue development against phenotypic variability. This discovery marks an important advance in our comprehension of how miRNAs function in the development of higher organisms.

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.

Summary MicroRNAs (miRNAs) are typically generated as ∼22-nucleotide double-stranded RNAs via the processing of precursor hairpins by the ribonuclease III enzyme Dicer, after which they are loaded into Argonaute (Ago) proteins to form an RNA-induced silencing complex (RISC). However, the biogenesis of miR-451, an erythropoietic miRNA conserved in vertebrates, occurs independently of Dicer and instead requires cleavage of the 3′ arm of the pre-miR-451 precursor hairpin by Ago2. The 3′ end of the Ago2-cleaved pre-miR-451 intermediate is then trimmed to the mature length by an unknown nuclease. Here, using a classical chromatographic approach, we identified poly(A)-specific ribonuclease (PARN) as the enzyme responsible for the 3′–5′ exonucleolytic trimming of Ago2-cleaved pre-miR-451. Surprisingly, our data show that trimming of Ago2-cleaved precursor miRNAs is not essential for target silencing, indicating that RISC is functional with miRNAs longer than the mature length. Our findings define the maturation step in the miRNA biogenesis pathway that depends on Ago2-mediated cleavage.

This protocol describes how to generate and genotype mutants using an optimized CRISPR–Cas9 genome-editing system in zebrafish (CRISPRscan). Because single guide RNAs (sgRNAs) have variable efficiency when targeting specific loci, our protocol starts by explaining how to use the web tool CRISPRscan to design highly efficient sgRNAs. The CRISPRscan algorithm is based on the results of an integrated analysis of more than 1000 sgRNAs in zebrafish, which uncovered highly predictive factors that influence Cas9 activity. Next, we describe how to easily generate sgRNAs in vitro, which can then be injected in vivo to target specific loci. The use of highly efficient sgRNAs can lead to biallelic mutations in the injected embryos, causing lethality. We explain how targeting Cas9 to the germline increases viability by reducing somatic mutations. Finally, we combine two methods to identify F1 heterozygous fish carrying the desired mutations: (i) Mut-Seq, a method based on high-throughput sequencing to detect F0 founder fish; and (ii) a polymerase chain reaction–based fragment analysis method that identifies F1 heterozygous fish characterized by Mut-Seq. In summary, this protocol includes the steps to generate and characterize mutant zebrafish lines using the CRISPR–Cas9 genome engineering system.

Gene expression is regulated extensively at the level 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 based on their activity in vivo and uses high-throughput sequencing to provide quantitative measurement of their regulatory function with 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). The RESA platform can be broadly applicable to uncover the regulatory features shaping gene expression and cellular function.

During the maternal-to-zygotic transition, maternal mRNAs are cleared by multiple distinct but interrelated pathways. A recent study in Nature by Zhao et al. (2017) finds that YTHDF2, a reader of N6- methylation, facilitates maternal mRNA decay, introducing an additional facet of control over transcript fate and developmental reprogramming.

MicroRNAs (miRNAs) regulate gene expression by pairing with complementary sequences in the 3' untranslated regions (UTRs) of transcripts. Although the molecular mechanism underlying miRNA biogenesis and activity is becoming better understood, determining the physiological role of the interaction of an miRNA with its target remains a challenge. A number of methods have been developed to inhibit individual miRNAs, but it can be difficult to determine which specific targets are responsible for any observed phenotypes. To address this problem, we use target protector (TP) morpholinos that interfere with a single miRNA-mRNA pair by binding specifically to the miRNA target sequence in the 3' UTR. In this protocol, we detail the steps for identifying miRNA targets, validating their regulation and using TPs to interrogate their function in zebrafish. Depending on the biological context, this set of experiments can be completed in 6-8 weeks.

microRNAs function as genetic rheostats to control gene output. Based on their role as modulators, it has been postulated that microRNAs canalize development and provide genetic robustness. Here, we uncover a novel regulatory layer of chemokine signaling by microRNAs that confers genetic robustness on primordial-germ-cell (PGC) migration. In zebrafish, PGCs are guided to the gonad by the ligand Sdf1a, which is regulated by sequestration receptor Cxcr7b. We find that miR-430 regulates sdf1a- and cxcr7-mRNAs. Using Target Protectors, we demonstrate that miR-430-mediated regulation of endogenous sdf1a and cxcr7b (i) facilitates dynamic expression of sdf1a by clearing its mRNA from previous expression domains, (ii) modulates the levels of the decoy receptor Cxcr7b to avoid excessive depletion of Sdf1a and (iii) buffers against variation in gene dosage of chemokine signaling components to ensure accurate PGC migration. Our results indicate that losing microRNA-mediated regulation can expose otherwise buffered genetic lesions leading to developmental defects.

Nature 535, 294–298 (2016); doi:10.1038/nature18614 In the Author Information section of this Letter, LT571435 from the SRA was provided as the accession number for the PacBio whole-genome sequencing data. Immediately after publication, we noted that LT571435 is instead the GenBank entry for the assembled genomic sequence of the cloche locus and this information is now included in the legend of Extended Data Fig. 2.