Identification of the zebrafish maternal and paternal transcriptomes

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Role of ZFHX4 in orofacial clefting based on human genetic data and zebrafish models

Orofacial clefting (OFC) is a frequent congenital anomaly and can occur either in the context of underlying syndromes or in isolation (nonsyndromic). The two common OFC phenotypes are cleft lip with/without cleft palate (CL/P) and cleft palate only (CPO). In this study, we searched for penetrant CL/P genes, by evaluating de novo copy number variants (CNV) from an exome sequencing dataset of 50 nonsyndromic patient-parent trios. We detected a heterozygous 86 kb de novo deletion affecting exons 4-11 of ZFHX4, a gene previously associated with OFC. Genetic and phenotypic data from our in-house and the AGORA cohort (710 and 229 individuals with nonsyndromic CL/P) together with literature and database reviews demonstrate that ZFHX4 variants can lead to both nonsyndromic and syndromic forms not only of CL/P but also CPO. Expression analysis in published single-cell RNA-sequencing data (mouse embryo, zebrafish larva) at relevant time-points support an important role of Zfhx4/zfhx4 in craniofacial development. To characterize the role of zfhx4 in zebrafish craniofacial development, we knocked out/down the zebrafish orthologue. Cartilage staining of the zfhx4 CRISPR F0 knockout and morpholino knockdown at 4 days post-fertilization showed an underdeveloped and abnormally shaped ethmoid plate and cartilaginous jaw (resembling micrognathia). While there is evidence for the dominant inheritance of ZFHX4 variants in OFC, we here present a patient with a possible recessive inheritance. In conclusion, ZFHX4 has a highly heterogeneous phenotypic spectrum and variable mode of inheritance. Our data highlight that ZFHX4 should be considered in genetic testing in patients with nonsyndromic clefting.

Ribosomal RNA Depletion for Poly(A)-Tail-Independent Quantification of Genome Activation

High-throughput RNA sequencing (RNA-seq) is commonly used to quantify gene expression transcriptome-wide. While usually paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA) to the exclusion of highly abundant ribosomal RNA (rRNA) in the cell, this strategy will under-quantify mRNA with short or absent poly(A) tails and can conflate changes in poly(A) tail length with changes in RNA level. This is notably an issue during early development, when cytoplasmic polyadenylation of maternal mRNA over time can be mistaken for genome activation in poly(A) + RNA-seq time courses. Here, we present a method to perform total RNA-seq using a streamlined rRNA depletion strategy customizable to any taxon. Antisense DNA oligos are designed with the aid of our Oligo-ASST web tool to sparsely tile the length of the rRNA, which are used with thermostable RNaseH to digest rRNA from a total RNA sample. After column cleanup, the mRNA-enriched sample is ready for sequencing library construction.

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

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

Protein profiling of zebrafish embryos unmasks regulatory layers during early embryogenesis

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

Screening of functional maternal-specific chromatin regulators in early embryonic development of zebrafish

The early stages of embryonic development rely on maternal products for proper regulation. However, screening for functional maternal-specific factors is challenging due to the time- and labor-intensive nature of traditional approaches. Here, we combine a computational pipeline and F0 null mutant technology to screen for functional maternal-specific chromatin regulators in zebrafish embryogenesis and identify Mcm3l, Mcm6l, and Npm2a as playing essential roles in DNA replication and cell division. Our results contribute to understanding the molecular mechanisms underlying early embryo development and highlight the importance of maternal-specific chromatin regulators in this critical stage.

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

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

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

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

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

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

Removal of pomt1 in zebrafish leads to loss of α-dystroglycan glycosylation and dystroglycanopathy phenotypes

Biallelic mutations in Protein O-mannosyltransferase 1 (POMT1) are among the most common causes of a severe group of congenital muscular dystrophies (CMDs) known as dystroglycanopathies. POMT1 is a glycosyltransferase responsible for the attachment of a functional glycan mediating interactions between the transmembrane glycoprotein dystroglycan and its binding partners in the extracellular matrix (ECM). Disruptions in these cell-ECM interactions lead to multiple developmental defects causing brain and eye malformations in addition to CMD. Removing Pomt1 in the mouse leads to early embryonic death due to the essential role of dystroglycan during placental formation in rodents. Here, we characterized and validated a model of pomt1 loss of function in the zebrafish showing that developmental defects found in individuals affected by dystroglycanopathies can be recapitulated in the fish. We also discovered that pomt1 mRNA provided by the mother in the oocyte supports dystroglycan glycosylation during the first few weeks of development. Muscle disease, retinal synapse formation deficits, and axon guidance defects can only be uncovered during the first week post fertilization by generating knock-out embryos from knock-out mothers. Conversely, maternal pomt1 from heterozygous mothers was sufficient to sustain muscle, eye, and brain development only leading to loss of photoreceptor synapses at 30 days post fertilization. Our findings show that it is important to define the contribution of maternal mRNA while developing zebrafish models of dystroglycanopathies and that offspring generated from heterozygous and knock-out mothers can be used to differentiate the role of dystroglycan glycosylation in tissue formation and maintenance.

miR-430 regulates zygotic mRNA during zebrafish embryogenesis

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Zebrafish raptor mutation inhibits the activity of mTORC1, inducing craniofacial defects due to autophagy-induced neural crest cell death

The mechanistic target of rapamycin (mTOR) coordinates metabolism and cell growth with environmental inputs. mTOR forms two functional complexes: mTORC1 and mTORC2. Proper development requires both complexes but mTORC1 has unique roles in numerous cellular processes, including cell growth, survival and autophagy. Here, we investigate the function of mTORC1 in craniofacial development. We created a zebrafish raptor mutant via CRISPR/Cas9, to specifically disrupt mTORC1. The entire craniofacial skeleton and eyes were reduced in size in mutants; however, overall body length and developmental timing were not affected. The craniofacial phenotype associates with decreased chondrocyte size and increased neural crest cell death. We found that autophagy is elevated in raptor mutants. Chemical inhibition of autophagy reduced cell death and improved craniofacial phenotypes in raptor mutants. Genetic inhibition of autophagy, via mutation of the autophagy gene atg7, improved facial phenotypes in atg7;raptor double mutants, relative to raptor single mutants. We conclude that finely regulated levels of autophagy, via mTORC1, are crucial for craniofacial development.

Maternal vgll4a regulates zebrafish epiboly through Yap1 activity

Gastrulation in zebrafish embryos commences with the morphogenetic rearrangement of blastodermal cells, which undergo a coordinated spreading from the animal pole to wrap around the egg at the vegetal pole. This rearrangement, known as epiboly, relies on the orchestrated activity of maternal transcripts present in the egg, compensating for the gradual activation of the zygotic genome. Epiboly involves the mechano-transducer activity of yap1 but what are the regulators of yap1 activity and whether these are maternally or zygotically derived remain elusive. Our study reveals the crucial role of maternal vgll4a, a proposed Yap1 competitor, during zebrafish epiboly. In embryos lacking maternal/zygotic vgll4a (MZvgll4a), the progression of epiboly and blastopore closure is delayed. This delay is associated with the ruffled appearance of the sliding epithelial cells, decreased expression of yap1-downstream targets and transient impairment of the actomyosin ring at the syncytial layer. Our study also shows that, rather than competing with yap1, vgll4a modulates the levels of the E-cadherin/β-catenin adhesion complex at the blastomeres' plasma membrane and hence their actin cortex distribution. Taking these results together, we propose that maternal vgll4a acts at epiboly initiation upstream of yap1 and the E-cadherin/β-catenin adhesion complex, contributing to a proper balance between tissue tension/cohesion and contractility, thereby promoting a timely epiboly progression.

The zebrafish as a potential model for vaccine and adjuvant development

INTRODUCTION

Zebrafishes represent a proven model for human diseases and systems biology, exhibiting physiological and genetic similarities and having innate and adaptive immune systems. However, they are underexplored for human vaccinology, vaccine development, and testing. Here we summarize gaps and challenges.

AREAS COVERED

Zebrafish models have four potential applications: 1) Vaccine safety: The past successes in using zebrafishes to test xenobiotics could extend to vaccine and adjuvant formulations for general safety or target organs due to the zebrafish embryos' optical transparency. 2) Innate immunity: The zebrafish offers refined ways to examine vaccine effects through signaling via Toll-like or NOD-like receptors in zebrafish myeloid cells. 3) Adaptive immunity: Zebrafishes produce IgM, IgD,and two IgZ immunoglobulins, but these are understudied, due to a lack of immunological reagents for challenge studies. 4) Systems vaccinology: Due to the availability of a well-referenced zebrafish genome, transcriptome, proteome, and epigenome, this model offers potential here.

EXPERT OPINION

It remains unproven whether zebrafishes can be employed for testing and developing human vaccines. We are still at the hypothesis-generating stage, although it is possible to begin outlining experiments for this purpose. Through transgenic manipulation, zebrafish models could offer new paths for shaping animal models and systems vaccinology.

Mature mRNA processing that deletes 3' end sequences directs translational activation and embryonic development

Eggs accumulate thousands of translationally repressed mRNAs that are translated into proteins after fertilization to direct diverse developmental processes. However, molecular mechanisms underlying the translation of stored mRNAs after fertilization remain unclear. Here, we report a previously unknown RNA processing of 3' end sequences of mature mRNAs that activates the translation of stored mRNAs. Specifically, 9 to 72 nucleotides at the 3' ends of zebrafish pou5f3 and mouse Pou5f1 mRNAs were deleted in the early stages of development. Reporter assays illustrated the effective translation of the truncated forms of mRNAs. Moreover, promotion and inhibition of the shortening of 3' ends accelerated and attenuated Pou5f3 accumulation, respectively, resulting in defective development. Identification of proteins binding to unprocessed and/or processed mRNAs revealed that mRNA shortening acts as molecular switches. Comprehensive analysis revealed that >250 mRNAs underwent this processing. Therefore, our results provide a molecular principle that triggers the translational activation and directs development.

The translation initiation factor homolog eif4e1c regulates cardiomyocyte metabolism and proliferation during heart regeneration

The eIF4E family of translation initiation factors bind 5' methylated caps and act as the limiting step for mRNA translation. The canonical eIF4E1A is required for cell viability, yet other related eIF4E families exist and are utilized in specific contexts or tissues. Here, we describe a family called Eif4e1c, for which we find roles during heart development and regeneration in zebrafish. The Eif4e1c family is present in all aquatic vertebrates but is lost in all terrestrial species. A core group of amino acids shared over 500 million years of evolution forms an interface along the protein surface, suggesting that Eif4e1c functions in a novel pathway. Deletion of eif4e1c in zebrafish caused growth deficits and impaired survival in juveniles. Mutants surviving to adulthood had fewer cardiomyocytes and reduced proliferative responses to cardiac injury. Ribosome profiling of mutant hearts demonstrated changes in translation efficiency of mRNA for genes known to regulate cardiomyocyte proliferation. Although eif4e1c is broadly expressed, its disruption had most notable impact on the heart and at juvenile stages. Our findings reveal context-dependent requirements for translation initiation regulators during heart regeneration.

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.

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

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

Maternal patterns of inheritance alter transcript expression in eggs

BACKGROUND

Modifications to early development can lead to evolutionary diversification. The early stages of development are under maternal control, as mothers produce eggs loaded with nutrients, proteins and mRNAs that direct early embryogenesis. Maternally provided mRNAs are the only expressed genes in initial stages of development and are tightly regulated. Differences in maternal mRNA provisioning could lead to phenotypic changes in embryogenesis and ultimately evolutionary changes in development. However, the extent that maternal mRNA expression in eggs can vary is unknown for most developmental models. Here, we use a species with dimorphic development- where females make eggs and larvae of different sizes and life-history modes-to investigate the extent of variation in maternal mRNA provisioning to the egg.

RESULTS

We find that there is significant variation in gene expression across eggs of different development modes, and that there are both qualitative and quantitative differences in mRNA expression. We separate parental effects from allelic effects, and find that both mechanisms contribute to mRNA expression differences. We also find that offspring of intraspecific crosses differentially provision their eggs based on the parental cross direction (a parental effect), which has not been previously demonstrated in reproductive traits like oogenesis.

CONCLUSION

We find that maternally controlled initiation of development is functionally distinct between eggs of different sizes and maternal genotypes. Both allele-specific effects and parent-of-origin effects contribute to gene expression differences in eggs. The latter indicates an intergenerational effect where a parent's genotype can affect gene expression in an egg made by the next generation.

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

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

Nanog organizes transcription bodies

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

Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata)

Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.

Sexually Dimorphic Gene Expression in X and Y Sperms Instructs Sexual Dimorphism of Embryonic Genome Activation in Yellow Catfish (Pelteobagrus fulvidraco)

Paternal factors play an important role in embryonic morphogenesis and contribute to sexual dimorphism in development. To assess the effect of paternal DNA on sexual dimorphism of embryonic genome activation, we compared X and Y sperm and different sexes of embryos before sex determination. Through transcriptome sequencing (RNA-seq) and whole-genome bisulfite sequencing (WGBS) of X and Y sperm, we found a big proportion of upregulated genes in Y sperm, supported by the observation that genome-wide DNA methylation level is slightly lower than in X sperm. Cytokine-cytokine receptor interaction, TGF-beta, and toll-like receptor pathways play important roles in spermatogenesis. Through whole-genome re-sequencing (WGRS) of parental fish and RNA-seq of five early embryonic stages, we found the low-blastocyst time point is a key to maternal transcriptome degradation and zygotic genome activation. Generally, sexual differences emerged from the bud stage. Moreover, through integrated analysis of paternal SNPs and gene expression, we evaluated the influence of paternal inheritance on sexual dimorphism of genome activation. Besides, we screened out gata6 and ddx5 as potential instructors for early sex determination and gonad development in yellow catfish. This work is meaningful for revealing the molecular mechanisms of sex determination and sexual dimorphism of fish species.

Dynamically regulated focal adhesions coordinate endothelial cell remodelling in developing vasculature

The assembly of a mature vascular network involves coordinated endothelial cell (EC) shape changes, including the process of EC elongation. How EC elongation is dynamically regulated in vivo is not fully understood. Here, we have generated a zebrafish mutant that is deficient for the integrin adaptor protein Talin 1 (Tln1). Using a new focal adhesion (FA) marker line expressing endothelial Vinculinb-eGFP, we demonstrate that EC FAs function dynamically and are lost in our tln1 mutants, allowing us to uncouple the primary roles of FAs in EC morphogenesis from the secondary effects that occur due to systemic vessel failure or loss of blood flow. Tln1 loss led to compromised F-actin rearrangements, perturbed EC elongation and disrupted cell-cell junction linearisation in vessel remodelling. Finally, chemical induction of actin polymerisation restored actin dynamics and EC elongation during vascular morphogenesis. Together, we identify that FAs are essential for EC elongation and junction linearisation in flow-pressured vessels and that they influence actin polymerisation in cellular morphogenesis. These observations can explain the severely compromised vessel beds and vascular leakage observed in mutant models that lack integrin signalling. This article has an associated 'The people behind the papers' interview.

Patterns and Crucial Regulation of Alternative Splicing During Early Development in Zebrafish

Many vertebrate genes generate multiple transcript variants that may encode functionally distinct protein isoforms, but the transcriptomes of various developmental stages are poorly defined. Identifying the transcriptome and its regulation during the normal developmental process is the key to deciphering the developmental stage-specific functions of genes. Here we presented a systematic assessment of the temporal alternative splicing (AS) events during the critical development stages to capture the dynamic gene expression changes and AS in zebrafish. An unexpected transcriptome complexity generated by AS was observed during zebrafish development. The patterns of AS events varied substantially among developmental stages despite the similarities in the total proportion of AS genes. We further found that AS afforded substantial functional diversification of genes through the generation of stage-specific AS events from broadly protein-coding genes as an essential developmental regulatory mechanism. Skipped exon (SE) showed the strongest signals among developmental AS (devAS), suggesting that devAS events generated by SE may be necessary for the normal development of zebrafish. Most developmental genes regulated by AS mechanisms were not modulated in terms of their overall expression levels, indicating that AS shaped the transcriptome independently from transcriptional regulation during development. 128-cell stage was a critical stage for gene transcription during embryonic development. Splicing factors as an essential developmental regulator underwent AS in the potential autoregulatory feedback loop and expressed multiple isoforms. Thus, zebrafish development was shaped by an interplay of programs controlling gene expression levels and AS. Overall, we provided a global view of developmental patterns of AS during zebrafish development and revealed that AS transitions were the crucial regulatory component of zebrafish embryonic development.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Translatome and transcriptome co-profiling reveals a role of TPRXs in human zygotic genome activation

Translational regulation plays a critical role during the oocyte-to-embryo transition (OET) and zygotic genome activation (ZGA). Here, we integrated ultra-low-input Ribo-seq with mRNA-seq to co-profile the translatome and transcriptome in human oocytes and early embryos. Comparison with mouse counterparts identified widespread differentially translated genes functioning in epigenetic reprogramming, transposon defense, and small RNA biogenesis, in part driven by species-specific regulatory elements in 3' untranslated regions. Moreover, PRD-like homeobox transcription factors, including TPRXL, TPRX1, and TPRX2, are highly translated around ZGA. TPRX1/2/L knockdown leads to defective ZGA and preimplantation development. Ectopically expressed TPRXs bind and activate key ZGA genes in human embryonic stem cells. These data reveal the conservation and divergence of translation landscapes during OET and identify critical regulators of human ZGA.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

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

Identification of embryonic RNA granules that act as sites of mRNA translation after changing their physical properties

Fertilized eggs begin to translate mRNAs at appropriate times and placements to control development, but how the translation is regulated remains unclear. Here, we found that pou5f3 mRNA encoding a transcriptional factor essential for development formed granules in a dormant state in zebrafish oocytes. Although the number of pou5f3 granules remained constant, Pou5f3 protein accumulated after fertilization. Intriguingly, signals of newly synthesized peptides and a ribosomal protein became colocalized with pou5f3 granules after fertilization and, moreover, nascent Pou5f3 was shown to be synthesized in the granules. This functional change was accompanied by changes in the state and internal structure of granules. Dissolution of the granules reduced the rate of protein synthesis. Similarly, nanog and sox19b mRNAs in zebrafish and Pou5f1/Oct4 mRNA in mouse assembled into granules. Our results reveal that subcellular compartments, termed embryonic RNA granules, function as activation sites of translation after changing physical properties for directing vertebrate development.

Pluripotency factors determine gene expression repertoire at zygotic genome activation

Awakening of zygotic transcription in animal embryos relies on maternal pioneer transcription factors. The interplay of global and specific functions of these proteins remains poorly understood. Here, we analyze chromatin accessibility and time-resolved transcription in single and double mutant zebrafish embryos lacking pluripotency factors Pou5f3 and Sox19b. We show that two factors modify chromatin in a largely independent manner. We distinguish four types of direct enhancers by differential requirements for Pou5f3 or Sox19b. We demonstrate that changes in chromatin accessibility of enhancers underlie the changes in zygotic expression repertoire in the double mutants. Pou5f3 or Sox19b promote chromatin accessibility of enhancers linked to the genes involved in gastrulation and ventral fate specification. The genes regulating mesendodermal and dorsal fates are primed for activation independently of Pou5f3 and Sox19b. Strikingly, simultaneous loss of Pou5f3 and Sox19b leads to premature expression of genes, involved in regulation of organogenesis and differentiation.

Zygotic genome activation in zebrafish relies on pluripotency transcription factors Pou5f3 and Sox19b. Here the authors investigate how these factors interact in vivo by analyzing the changes in chromatin state and time-resolved transcription in Pou5f3 and Sox19b single and double mutant embryos.

brca2-mutant zebrafish exhibit context- and tissue-dependent alterations in cell phenotypes and response to injury

Cancer cells frequently co-opt molecular programs that are normally activated in specific contexts, such as embryonic development and the response to injury. Determining the impact of cancer-associated mutations on cellular phenotypes within these discrete contexts can provide new insight into how such mutations lead to dysregulated cell behaviors and subsequent cancer onset. Here we assess the impact of heritable BRCA2 mutation on embryonic development and the injury response using a zebrafish model (Danio rerio). Unlike most mouse models for BRCA2 mutation, brca2-mutant zebrafish are fully viable and thus provide a unique tool for assessing both embryonic and adult phenotypes. We find that maternally provided brca2 is critical for normal oocyte development and embryonic survival in zebrafish, suggesting that embryonic lethality associated with BRCA2 mutation is likely to reflect defects in both meiotic and embryonic developmental programs. On the other hand, we find that adult brca2-mutant zebrafish exhibit aberrant proliferation of several cell types under basal conditions and in response to injury in tissues at high risk for cancer development. These divergent effects exemplify the often-paradoxical outcomes that occur in embryos (embryonic lethality) versus adult animals (cancer predisposition) with mutations in cancer susceptibility genes such as BRCA2. The altered cell behaviors identified in brca2-mutant embryonic and adult tissues, particularly in adult tissues at high risk for cancer, indicate that the effects of BRCA2 mutation on cellular phenotypes are both context- and tissue-dependent.

Developmental series of gene expression clarifies maternal mRNA provisioning and maternal-to-zygotic transition in a reef-building coral

BACKGROUND

Maternal mRNA provisioning of oocytes regulates early embryogenesis. Maternal transcripts are degraded as zygotic genome activation (ZGA) intensifies, a phenomenon known as the maternal-to-zygotic transition (MZT). Here, we examine gene expression over nine developmental stages in the Pacific rice coral, Montipora capitata, from eggs and embryos at 1, 4, 9, 14, 22, and 36 h-post-fertilization (hpf), as well as swimming larvae (9d), and adult colonies.

RESULTS

Weighted Gene Coexpression Network Analysis revealed four expression peaks, identifying the maternal complement, two waves of the MZT, and adult expression. Gene ontology enrichment revealed maternal mRNAs are dominated by cell division, methylation, biosynthesis, metabolism, and protein/RNA processing and transport functions. The first MZT wave occurs from ~4-14 hpf and is enriched in terms related to biosynthesis, methylation, cell division, and transcription. In contrast, functional enrichment in the second MZT wave, or ZGA, from 22 hpf-9dpf, includes ion/peptide transport and cell signaling. Finally, adult expression is enriched for functions related to signaling, metabolism, and ion/peptide transport. Our proposed MZT timing is further supported by expression of enzymes involved in zygotic transcriptional repression (Kaiso) and activation (Sox2), which peak at 14 hpf and 22 hpf, respectively. Further, DNA methylation writing (DNMT3a) and removing (TET1) enzymes peak and remain stable past ~4 hpf, suggesting that methylome programming occurs before 4 hpf.

CONCLUSIONS

Our high-resolution insight into the coral maternal mRNA and MZT provides essential baseline information to understand parental carryover effects and the sensitivity of developmental success under increasing environmental stress.

Core pathway proteins and the molecular basis of planar polarity in the zebrafish gastrula

The planar polarization of cells and subcellular structures is critical for embryonic development. Coordination of this polarity can provide cells a sense of direction in relation to the anterior-posterior and dorsal-ventral body axes. Fly epithelia use a core pathway comprised of transmembrane (Van Gogh/Strabismus, Frizzled, and Flamingo/Starry night) and cytoplasmic (Prickle or Spiny-legs, Dishevelled, and Diego) proteins to communicate directional information between cells and thereby promote the uniform orientation of structures such as hairs. In the zebrafish gastrula, planar polarity underlies complex cellular processes, including directed migration and intercalation, that are required to shape the embryo body. Like other vertebrates, the zebrafish genome encodes homologs of each core protein, and it is well-established that polarized gastrula cell behaviors are regulated by some of them. However, it is unknown whether a conserved six-member core protein pathway regulates planar polarity during zebrafish gastrulation. Here, we review our current understanding of core protein function as it relates to two specific examples of planar polarity, the dorsal convergence of lateral gastrula cells and the mediolateral intercalation of midline cells. We consider the hallmarks of fly planar polarity and discuss data regarding asymmetric protein localization and function, and the intercellular communication of polarity information.

Conserved germ plasm characteristics across the Danio and Devario lineages

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

RNA-induced inflammation and migration of precursor neurons initiates neuronal circuit regeneration in zebrafish

Tissue regeneration and functional restoration after injury are considered as stem- and progenitor-cell-driven processes. In the central nervous system, stem cell-driven repair is slow and problematic because function needs to be restored rapidly for vital tasks. In highly regenerative vertebrates, such as zebrafish, functional recovery is rapid, suggesting a capability for fast cell production and functional integration. Surprisingly, we found that migration of dormant "precursor neurons" to the injury site pioneers functional circuit regeneration after spinal cord injury and controls the subsequent stem-cell-driven repair response. Thus, the precursor neurons make do before the stem cells make new. Furthermore, RNA released from the dying or damaged cells at the site of injury acts as a signal to attract precursor neurons for repair. Taken together, our data demonstrate an unanticipated role of neuronal migration and RNA as drivers of neural repair.

Loss of zebrafish atp6v1e1b, encoding a subunit of vacuolar ATPase, recapitulates human ARCL type 2C syndrome and identifies multiple pathobiological signatures

The inability to maintain a strictly regulated endo(lyso)somal acidic pH through the proton-pumping action of the vacuolar-ATPases (v-ATPases) has been associated with various human diseases including heritable connective tissue disorders. Autosomal recessive (AR) cutis laxa (CL) type 2C syndrome is associated with genetic defects in the ATP6V1E1 gene and is characterized by skin wrinkles or loose redundant skin folds with pleiotropic systemic manifestations. The underlying pathological mechanisms leading to the clinical presentations remain largely unknown. Here, we show that loss of atp6v1e1b in zebrafish leads to early mortality, associated with craniofacial dysmorphisms, vascular anomalies, cardiac dysfunction, N-glycosylation defects, hypotonia, and epidermal structural defects. These features are reminiscent of the phenotypic manifestations in ARCL type 2C patients. Our data demonstrates that loss of atp6v1e1b alters endo(lyso)somal protein levels, and interferes with non-canonical v-ATPase pathways in vivo. In order to gain further insights into the processes affected by loss of atp6v1e1b, we performed an untargeted analysis of the transcriptome, metabolome, and lipidome in early atp6v1e1b-deficient larvae. We report multiple affected pathways including but not limited to oxidative phosphorylation, sphingolipid, fatty acid, and energy metabolism together with profound defects on mitochondrial respiration. Taken together, our results identify complex pathobiological effects due to loss of atp6v1e1b in vivo.

Cutis laxa syndromes are pleiotropic disorders of the connective tissue, characterized by skin redundancy and variable systemic manifestations. Cutis laxa syndromes are caused by pathogenic variants in genes encoding structural and regulatory components of the extracellular matrix or in genes encoding components of cellular trafficking, metabolism, and mitochondrial function. Pathogenic variants in genes coding for vacuolar-ATPases, a multisubunit complex responsible for the acidification of multiple intracellular vesicles, cause type 2 cutis laxa syndromes, a group of cutis laxa subtypes further characterized by neurological, skeletal, and rarely cardiopulmonary manifestations. To investigate the pathomechanisms of vacuolar-ATPase dysfunction, we generated zebrafish models that lack a crucial subunit of the vacuolar-ATPases. The mutant zebrafish models show morphological and functional features reminiscent of the phenotypic manifestations in cutis laxa patients carrying pathogenic variants in ATP6V1E1. In-depth analysis at multiple -omic levels identified biological signatures that indicate impairment of signaling pathways, lipid metabolism, and mitochondrial respiration. We anticipate that these data will contribute to a better understanding of the pathogenesis of cutis laxa syndromes and other disorders involving defective v-ATPase function, which may eventually improve patient treatment and management.

Nanopore long-read RNA-seq and absolute quantification delineate transcription dynamics in early embryo development of an insect pest

The olive fruit fly, Bactrocera oleae, is the most important pest for the olive fruit but lacks adequate transcriptomic characterization that could aid in molecular control approaches. We apply nanopore long-read RNA-seq with internal RNA standards allowing absolute transcript quantification to analyze transcription dynamics during early embryo development for the first time in this organism. Sequencing on the MinION platform generated over 31 million reads. Over 50% of the expressed genes had at least one read covering its entire length validating our full-length approach. We generated a de novo transcriptome assembly and identified 1768 new genes and a total of 79,810 isoforms; a fourfold increase in transcriptome diversity compared to the current NCBI predicted transcriptome. Absolute transcript quantification per embryo allowed an insight into the dramatic re-organization of maternal transcripts. We further identified Zelda as a possible regulator of early zygotic genome activation in B. oleae and provide further insights into the maternal-to-zygotic transition. These data show the utility of long-read RNA in improving characterization of non-model organisms that lack a fully annotated genome, provide potential targets for sterile insect technic approaches, and provide the first insight into the transcriptome landscape of the developing olive fruit fly embryo.

The adaptor protein Grb2b is an essential modulator for lympho-venous sprout formation in the zebrafish trunk

Vegfc/Vegfr3 signaling is critical for lymphangiogenesis, the sprouting of lymphatic vessels. In zebrafish, cells sprouting from the posterior cardinal vein can either form lymphatic precursor cells or contribute to intersegmental vein formation. Both, the Vegfc-dependent differential induction of Prox1a in sprouting cells as well as a Notch-mediated pre-pattern within intersegmental vessels have been associated with the regulation of secondary sprout behavior. However, how exactly a differential lymphatic versus venous sprout cell behavior is achieved is not fully understood. Here, we characterize a zebrafish mutant in the adaptor protein Grb2b, and demonstrate through genetic interaction studies that Grb2b acts within the Vegfr3 pathway. Mutant embryos exhibit phenotypes that are consistent with reduced Vegfr3 signaling outputs prior to the sprouting of endothelial cells from the vein. During secondary sprouting stages, loss of grb2b leads to defective cell behaviors resulting in a loss of parachordal lymphangioblasts, while only partially affecting the number of intersegmental veins. A second GRB2 zebrafish ortholog, grb2a, contributes to the development of lymphatic structures in the meninges and in the head, but not in the trunk. Our results illustrate an essential role of Grb2b in vivo for cell migration to the horizontal myoseptum and for the correct formation of the lymphatic vasculature, while being less critically required in intersegmental vein formation. Thus, there appear to be higher requirements for Grb2b and therefore Vegfr3 downstream signaling levels in lymphatic versus vein precursor-generating sprouts.

Germ cell differentiation requires Tdrd7-dependent chromatin and transcriptome reprogramming marked by germ plasm relocalization

In many animal models, primordial germ cell (PGC) development depends on maternally deposited germ plasm, which prevents somatic cell fate. Here, we show that PGCs respond to regulatory information from the germ plasm in two distinct phases using two distinct mechanisms in zebrafish. We demonstrate that PGCs commence zygotic genome activation together with the somatic blastocysts with no demonstrable differences in transcriptional and chromatin opening. Unexpectedly, both PGC and somatic blastocysts activate germ-cell-specific genes, which are only stabilized in PGCs by cytoplasmic germ plasm determinants. Disaggregated perinuclear relocalization of germ plasm during PGC migration is regulated by the germ plasm determinant Tdrd7 and is coupled to dramatic divergence between PGC and somatic transcriptomes. This transcriptional divergence relies on PGC-specific cis-regulatory elements characterized by promoter-proximal distribution. We show that Tdrd7-dependent reconfiguration of chromatin accessibility is required for elaboration of PGC fate but not for PGC migration.

•

No evidence for transcriptional activation delay in zebrafish PGCs

•

Germ-plasm-associated post-transcriptional divergence during ZGA

•

Epigenetic reprogramming marks onset of PGC migration

•

Epigenetic reprogramming in PGCs relies on Tdrd7, coupled to germ plasm relocalization

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.

The X-linked acrogigantism-associated gene gpr101 is a regulator of early embryonic development and growth in zebrafish

We recently described X-linked acrogigantism (X-LAG), a condition of early childhood-onset pituitary gigantism associated with microduplications of the GPR101 receptor. The expression of GPR101 in hyperplastic pituitary regions and tumors in X-LAG patients, and GPR101's normally transient pituitary expression during fetal development, suggest a role in the regulation of growth. Nevertheless, little is still known about GPR101's physiological functions, especially during development. By using zebrafish models, we investigated the role of gpr101 during embryonic development and somatic growth. Transient ectopic gpr101 expression perturbed the embryonic body plan but did not affect growth. Loss of gpr101 led to a significant reduction in body size that was even more pronounced in the absence of maternal transcripts, as well as subfertility. These changes were accompanied by gastrulation and hypothalamic defects. In conclusion, both gpr101 loss- and gain-of-function affect, in different ways, fertility, embryonic patterning, growth and brain development.

Decoding Dishevelled-Mediated Wnt Signaling in Vertebrate Early Development

Dishevelled proteins are key players of Wnt signaling pathways. They transduce Wnt signals and perform cellular functions through distinct conserved domains. Due to the presence of multiple paralogs, the abundant accumulation of maternal transcripts, and the activation of distinct Wnt pathways, their regulatory roles during vertebrate early development and the mechanism by which they dictate the pathway specificity have been enigmatic and attracted much attention in the past decades. Extensive studies in different animal models have provided significant insights into the structure-function relationship of conserved Dishevelled domains in Wnt signaling and the implications of Dishevelled isoforms in early developmental processes. Notably, intra- and inter-molecular interactions and Dishevelled dosage may be important in modulating the specificity of Wnt signaling. There are also distinct and redundant functions among Dishevelled isoforms in development and disease, which may result from differential spatiotemporal expression patterns and biochemical properties and post-translational modifications. This review presents the advances and perspectives in understanding Dishevelled-mediated Wnt signaling during gastrulation and neurulation in vertebrate early embryos.

The Zebrafish Perivitelline Fluid Provides Maternally-Inherited Defensive Immunity

In the teleost egg, the embryo is immersed in an extraembryonic fluid that fills the space between the embryo and the chorion and partially isolates it from the external environment, called the perivitelline fluid (PVF). The exact composition of the PVF remains unknown in vertebrate animals. The PVF allows the embryo to avoid dehydration, to maintain a safe osmotic balance and provides mechanical protection; however, its potential defensive properties against bacterial pathogens has not been reported. In this work, we determined the global proteomic profile of PVF in zebrafish eggs and embryos, and the maternal or zygotic origin of the identified proteins was studied. In silico analysis of PVF protein composition revealed an enrichment of protein classes associated with non-specific humoral innate immunity. We found lectins, protease inhibitors, transferrin, and glucosidases present from early embryogenesis until hatching. Finally, in vitro and in vivo experiments done with this fluid demonstrated that the PVF possessed a strong agglutinating capacity on bacterial cells and protected the embryos when challenged with the pathogenic bacteria Edwardsiella tarda. Our results suggest that the PVF is a primitive inherited immune extraembryonic system that protects the embryos from external biological threats prior to hatching.

Maternal contributions to gastrulation in zebrafish

Gastrulation is a critical early morphogenetic process of animal development, during which the three germ layers; mesoderm, endoderm and ectoderm, are rearranged by internalization movements. Concurrent epiboly movements spread and thin the germ layers while convergence and extension movements shape them into an anteroposteriorly elongated body with head, trunk, tail and organ rudiments. In zebrafish, gastrulation follows the proliferative and inductive events that establish the embryonic and extraembryonic tissues and the embryonic axis. Specification of these tissues and embryonic axes are controlled by the maternal gene products deposited in the egg. These early maternally controlled processes need to generate sufficient cell numbers and establish the embryonic polarity to ensure normal gastrulation. Subsequently, after activation of the zygotic genome, the zygotic gene products govern mesoderm and endoderm induction and germ layer patterning. Gastrulation is initiated during the maternal-to-zygotic transition, a process that entails both activation of the zygotic genome and downregulation of the maternal transcripts. Genomic studies indicate that gastrulation is largely controlled by the zygotic genome. Nonetheless, genetic studies that investigate the relative contributions of maternal and zygotic gene function by comparing zygotic, maternal and maternal zygotic mutant phenotypes, reveal significant contribution of maternal gene products, transcripts and/or proteins, that persist through gastrulation, to the control of gastrulation movements. Therefore, in zebrafish, the maternally expressed gene products not only set the stage for, but they also actively participate in gastrulation morphogenesis.

Transcriptome Analysis of Maternal Gene Transcripts in Unfertilized Eggs of Misgurnus anguillicaudatus and Identification of Immune-Related Maternal Genes

Maternal genes are important in directing early development and determining egg quality in fish. We here report the de novo transcriptome from four tissue libraries of the cyprinid loach, Misgurnus anguillicaudatus, and for the first time identified maternal gene transcripts in unfertilized eggs and suggest their immune system involvement. Expression profiles and functional enrichment revealed a total 24,116 transcripts were expressed as maternal transcripts in unfertilized eggs, which were involved in a wide range of biological functions and pathways. Comparison expression profiles and analysis of tissue specificity revealed that the large numbers of maternal transcripts were stored in unfertilized eggs near the late phase of ovarian maturation and before ovulation. Functional classification showed a total of 279 maternal immune-related transcripts classified with immune system process GO term and immune system KEGG pathway. qPCR analysis showed that transcript levels of identified maternal immune-related candidate genes were dynamically modulated during development and early ontogeny of M. anguillicaudatus. Taken together, this study could not only provide knowledge on the protective roles of maternal immune-related genes during early life stage of M. anguillicaudatus but could also be a valuable transcriptomic/genomic resource for further analysis of maternally provisioned genes in M. anguillicaudatus and other related teleost fishes.

Tissue-Specific Requirement for the GINS Complex During Zebrafish Development

Efficient and accurate DNA replication is particularly critical in stem and progenitor cells for successful proliferation and survival. The replisome, an amalgam of protein complexes, is responsible for binding potential origins of replication, unwinding the double helix, and then synthesizing complimentary strands of DNA. According to current models, the initial steps of DNA unwinding and opening are facilitated by the CMG complex, which is composed of a GINS heterotetramer that connects Cdc45 with the mini-chromosome maintenance (Mcm) helicase. In this work, we provide evidence that in the absence of GINS function DNA replication is cell autonomously impaired, and we also show that gins1 and gins2 mutants exhibit elevated levels of apoptosis restricted to actively proliferating regions of the central nervous system (CNS). Intriguingly, our results also suggest that the rapid cell cycles during early embryonic development in zebrafish may not require the function of the canonical GINS complex as neither zygotic Gins1 nor Gins2 isoforms seem to be present during these stages.

RNA structural dynamics regulate early embryogenesis through controlling transcriptome fate and function

BACKGROUND

Vertebrate early embryogenesis is initially directed by a set of maternal RNAs and proteins, yet the mechanisms controlling this program remain largely unknown. Recent transcriptome-wide studies on RNA structure have revealed its pervasive and crucial roles in RNA processing and functions, but whether and how RNA structure regulates the fate of the maternal transcriptome have yet to be determined.

RESULTS

Here we establish the global map of four nucleotide-based mRNA structures by icSHAPE during zebrafish early embryogenesis. Strikingly, we observe that RNA structurally variable regions are enriched in the 3' UTR and contain cis-regulatory elements important for maternal-to-zygotic transition (MZT). We find that the RNA-binding protein Elavl1a stabilizes maternal mRNAs by binding to the cis-elements. Conversely, RNA structure formation suppresses Elavl1a's binding leading to the decay of its maternal targets.

CONCLUSIONS

Our study finds that RNA structurally variable regions are enriched in mRNA 3' UTRs and contain cis-regulatory elements during zebrafish early embryogenesis. We reveal that Elavl1a regulates maternal RNA stability in an RNA structure-dependent fashion. Overall, our findings reveal a broad and fundamental role of RNA structure-based regulation in vertebrate early embryogenesis.

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

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

Functional Characterization and Expression Analyses Show Differential Roles of Maternal and Zygotic Dgcr8 in Early Embryonic Development

Dgcr8 is involved in the biogenesis of canonical miRNAs to process pri-miRNA into pre-miRNA. Previous studies have provided evidence that Dgcr8 plays an essential role in different biological processes. However, the function of maternal and zygotic Dgcr8 in early embryonic development remains largely unknown. Recently, we have reported a novel approach for generating germline-specific deletions in zebrafish. This germline knockout model offers an opportunity to investigate into the differential roles of maternal or zygotic Dgcr8. Although germline specific dgcr8 deletion has no influence on gonad development, maternal or zygotic dgcr8 is essential for embryonic development in the offspring. Both maternal dgcr8 (Mdgcr8) and maternal zygotic dgcr8 (MZdgcr8) mutants display multiple developmental defects and die within 1 week. Moreover, MZdcgr8 mutant displays more severe morphogenesis defects. However, when a miR-430 duplex (the most abundantly expressed miRNA in early embryonic stage) is used to rescue the maternal mutant phenotype, the Mdgcr8 embryos could be rescued successfully and grow into adulthood and achieve sexual maturation, whereas the MZdgcr8 embryos are only partially rescued and they all die within 1 week. The differential phenotypes between the Mdgcr8 and MZdgcr8 embryos provide us with an opportunity to study the roles of individual miRNAs during early development.

Globins and nitric oxide homeostasis in fish embryonic development

Since the discovery of new members of the globin superfamily such as Cytoglobin, Neuroglobin and Globin X, in addition to the most well-known members, Hemoglobin and Myoglobin, different hypotheses have been suggested about their function in vertebrates. Globins are ubiquitously found in living organisms and can carry out different functions based on their ability to bind ligands such as O2, and nitric oxide (NO) and to catalyze reactions scavenging NO or generating NO by reducing nitrite. NO is a highly diffusible molecule with a central role in signaling important for egg maturation, fertilization and early embryonic development. The globins ability to scavenge or generate NO makes these proteins ideal candidates in regulating NO homeostasis depending on the micro environment and tissue NO demands. Different amounts of various globins have been found in zebrafish eggs and developing embryos where it's unlikely that they function as respiratory proteins and instead could play a role in maintaining embryonic NO homeostasis. Here we summarize the current knowledge concerning the role of NO in adult fish in comparison to mammals and we discuss NO function during embryonic development with possible implications for globins in maintaining embryonic NO homeostasis.