A high-resolution mRNA expression time course of embryonic development in zebrafish

Establishment and transcriptome analysis of single blastomere-derived cell lines from zebrafish

Maintaining chromosome euploidy in zebrafish embryonic cells is challenging because of the degradation of genomic integrity during cell passaging. In this study, we report the derivation of zebrafish cell lines from single blastomeres. These cell lines have a stable chromosome status attributed to BMP4 and exhibit continuous proliferation in vitro. Twenty zebrafish cell lines are successfully established from single blastomeres. Single-cell transcriptome sequencing analysis confirms the fidelity of gene expression profiles throughout long-term culturing of at least 45 passages. The long-term cultured cells are specialized into epithelial cells, exhibiting similar expression patterns validated by integrative transcriptomic analysis. Overall, this work provides a protocol for establishing zebrafish cell lines from single blastomeres, which can serve as valuable tools for in vitro investigations of epithelial cell dynamics in terms of life-death balance and cell fate determination during normal homeostasis.

Rhabdomyosarcoma fusion oncoprotein initially pioneers a neural signature in vivo

Fusion-positive rhabdomyosarcoma is an aggressive pediatric cancer molecularly characterized by arrested myogenesis. The defining genetic driver, PAX3::FOXO1, functions as a chimeric gain-of-function transcription factor. An incomplete understanding of PAX3::FOXO1's in vivo epigenetic mechanisms has hindered therapeutic development. Here, we establish a PAX3::FOXO1 zebrafish injection model and semi-automated ChIP-seq normalization strategy to evaluate how PAX3::FOXO1 initially interfaces with chromatin in a developmental context. We investigated PAX3::FOXO1's recognition of chromatin and subsequent transcriptional consequences. We find that PAX3::FOXO1 interacts with inaccessible chromatin through partial/homeobox motif recognition consistent with pioneering activity. However, PAX3::FOXO1-genome binding through a composite paired-box/homeobox motif alters chromatin accessibility and redistributes H3K27ac to activate neural transcriptional programs. We uncover neural signatures that are highly representative of clinical rhabdomyosarcoma gene expression programs that are enriched following chemotherapy. Overall, we identify partial/homeobox motif recognition as a new mode for PAX3::FOXO1 pioneer function and identify neural signatures as a potentially critical PAX3::FOXO1 tumor initiation event.

Mitochondrial abnormalities contribute to muscle weakness in a Dnajb6 deficient zebrafish model

Mutations in DNAJB6 are a well-established cause of limb girdle muscular dystrophy type D1 (LGMD D1). Patients with LGMD D1 develop progressive muscle weakness with histology showing fibre damage, autophagic vacuoles, and aggregates. Whilst there are many reports of LGMD D1 patients, the role of DNAJB6 in the muscle is still unclear. In this study, we developed a loss of function zebrafish model in order to investigate the role of Dnajb6. Using a double dnajb6a and dnajb6b mutant model, we show that loss of Dnajb6 leads to a late onset muscle weakness. Interestingly, we find that adult fish lacking Dnajb6 do not have autophagy or myofibril defects, however, they do show mitochondrial changes and damage. This study demonstrates that loss of Dnajb6 causes mitochondrial defects and suggests that this contributes to muscle weakness in LGMD D1. These findings expand our knowledge of the role of Dnajb6 in the muscle and provides a model to screen novel therapies for LGMD D1.

Pathogenic variants in TMEM184B cause a neurodevelopmental syndrome via alteration of metabolic signaling

Transmembrane protein 184B (TMEM184B) is an endosomal 7-pass transmembrane protein with evolutionarily conserved roles in synaptic structure and axon degeneration. We report six pediatric patients who have de novo heterozygous variants in TMEM184B. All individuals harbor rare missense or mRNA splicing changes and have neurodevelopmental deficits including intellectual disability, corpus callosum hypoplasia, seizures, and/or microcephaly. TMEM184B is predicted to contain a pore domain, wherein many human disease-associated variants cluster. Structural modeling suggests that all missense variants alter TMEM184B protein stability. To understand the contribution of TMEM184B to neural development in vivo, we suppressed the TMEM184B ortholog in zebrafish and observed microcephaly and reduced anterior commissural neurons, aligning with patient symptoms. Ectopic TMEM184B expression resulted in dominant effects for K184E and G162R. However, in vivo complementation studies demonstrate that all other variants tested result in diminished protein function and indicate a haploinsufficiency basis for disease. Expression of K184E and other variants increased apoptosis in cell lines and altered nuclear localization of transcription factor EB (TFEB), a master regulator of lysosomal biogenesis, suggesting disrupted nutrient signaling pathways. Together, our data indicate that TMEM184B variants cause cellular metabolic disruption likely through divergent molecular effects that all result in abnormal neural development.

Mapping the chromatin accessibility landscape of zebrafish embryogenesis at single-cell resolution by SPATAC-seq

Currently, the dynamic accessible elements that determine regulatory programs responsible for the unique identity and function of each cell type during zebrafish embryogenesis lack detailed study. Here we present SPATAC-seq: a split-pool ligation-based assay for transposase-accessible chromatin using sequencing. Using SPATAC-seq, we profiled chromatin accessibility in more than 800,000 individual nuclei across 20 developmental stages spanning the sphere stage to the early larval protruding mouth stage. Using this chromatin accessibility map, we identified 604 cell states and inferred their developmental relationships. We also identified 959,040 candidate cis-regulatory elements (cCREs) and delineated development-specific cCREs, as well as transcription factors defining diverse cell identities. Importantly, enhancer reporter assays confirmed that the majority of tested cCREs exhibited robust enhanced green fluorescent protein expression in restricted cell types or tissues. Finally, we explored gene regulatory programs that drive pigment and notochord cell differentiation. Our work provides a valuable open resource for exploring driver regulators of cell fate decisions in zebrafish embryogenesis.

Transcriptomic profiling of gastrointestinal tracts in dairy cattle during lactation reveals molecular adaptations for milk synthesis

During lactation, dairy cattle's digestive tract requires significant adaptations to meet the increased nutrient demands for milk production. As we attempt to improve milk-related traits through selective pressure, it is crucial to understand the biological functions of the epithelia of the rumen, small intestine, and colonic tissues in response to changes in physiological state driven by changes in nutrient demands for milk synthesis. In this study, we obtained a total of 108 transcriptome profiles from three tissues (epithelia of the colon, duodenum, and rumen) of five Holstein cows, spanning eight time points from the early, mid, late lactation periods to the dry period. On average 97.06% of reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We analyzed 27,607 gene expression patterns at multiple periods, enabling direct comparisons within and among tissues during different lactation stages, including early and peak lactation. We identified 1645, 813, and 2187 stage-specific genes in the colon, duodenum, and rumen, respectively, which were enriched for common or specific biological functions among different tissues. Time series analysis categorized the expressed genes within each tissue into four clusters. Furthermore, when the three tissues were analyzed collectively, 36 clusters of similarly expressed genes were identified. By integrating other comprehensive approaches such as gene co-expression analyses, functional enrichment, and cell type deconvolution, we gained profound insights into cattle lactation, revealing tissue-specific characteristics of the gastrointestinal tract and shedding light on the intricate molecular adaptations involved in nutrient absorption, immune regulation, and cellular processes for milk synthesis during lactation.

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.

CRIMP: a CRISPR/Cas9 insertional mutagenesis protocol and toolkit

Site-directed insertion is a powerful approach for generating mutant alleles, but low efficiency and the need for customisation for each target has limited its application. To overcome this, we developed a highly efficient targeted insertional mutagenesis system, CRIMP, and an associated plasmid toolkit, CRIMPkit, that disrupts native gene expression by inducing complete transcriptional termination, generating null mutant alleles without inducing genetic compensation. The protocol results in a high frequency of integration events and can generate very early targeted insertions, during the first cell division, producing embryos with expression in one or both halves of the body plan. Fluorescent readout of integration events facilitates selection of successfully mutagenized fish and, subsequently, visual identification of heterozygous and mutant animals. Together, these advances greatly improve the efficacy of generating and studying mutant lines. The CRIMPkit contains 24 ready-to-use plasmid vectors to allow easy and complete mutagenesis of any gene in any reading frame without requiring custom sequences, modification, or subcloning.

Bioinformatic Analysis of Roquin Family Reveals Their Potential Role in Immune System

The Roquin family is a recognized RNA-binding protein family that plays vital roles in regulating the expression of pro-inflammatory target gene mRNA during the immune process in mammals. However, the evolutionary status of the Roquin family across metazoans remains elusive, and limited studies are found in fish species. In this study, we discovered that the RC3H genes underwent a single round of gene duplication from a primitive ancestor during evolution from invertebrates to vertebrates. Furthermore, there were instances of species-specific gene loss events or teleost lineage-specific gene duplications throughout evolution. Domain/motif organization and selective pressure analysis revealed that Roquins exhibit high homology both within members of the family within the same species and across species. The three rc3h genes in zebrafish displayed similar expression patterns in early embryos and adult tissues, with rc3h1b showing the most prominent expression among them. Additionally, the promoter regions of the zebrafish rc3h genes contained numerous transcription factor binding sites similar to those of mammalian homologs. Moreover, the interaction protein network of Roquin and the potential binding motif in the 3'-UTR of putative target genes analysis both indicated that Roquins have the potential to degrade target mRNA through mechanisms similar to those of mammalian homologs. These findings shed light on the evolutionary history of Roquin among metazoans and hypothesized their role in the immune systems of zebrafish.

Reusable tutorials for using cloud-based computing environments for the analysis of bacterial gene expression data from bulk RNA sequencing

This manuscript describes the development of a resource module that is part of a learning platform named "NIGMS Sandbox for Cloud-based Learning" https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox at the beginning of this Supplement. This module delivers learning materials on RNA sequencing (RNAseq) data analysis in an interactive format that uses appropriate cloud resources for data access and analyses. Biomedical research is increasingly data-driven, and dependent upon data management and analysis methods that facilitate rigorous, robust, and reproducible research. Cloud-based computing resources provide opportunities to broaden the application of bioinformatics and data science in research. Two obstacles for researchers, particularly those at small institutions, are: (i) access to bioinformatics analysis environments tailored to their research; and (ii) training in how to use Cloud-based computing resources. We developed five reusable tutorials for bulk RNAseq data analysis to address these obstacles. Using Jupyter notebooks run on the Google Cloud Platform, the tutorials guide the user through a workflow featuring an RNAseq dataset from a study of prophage altered drug resistance in Mycobacterium chelonae. The first tutorial uses a subset of the data so users can learn analysis steps rapidly, and the second uses the entire dataset. Next, a tutorial demonstrates how to analyze the read count data to generate lists of differentially expressed genes using R/DESeq2. Additional tutorials generate read counts using the Snakemake workflow manager and Nextflow with Google Batch. All tutorials are open-source and can be used as templates for other analysis.

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

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

Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

Shedding a Light on Dark Genes: A Comparative Expression Study of PRR12 Orthologues during Zebrafish Development

Haploinsufficiency of the PRR12 gene is implicated in a human neuro-ocular syndrome. Although identified as a nuclear protein highly expressed in the embryonic mouse brain, PRR12 molecular function remains elusive. This study explores the spatio-temporal expression of zebrafish PRR12 co-orthologs, prr12a and prr12b, as a first step to elucidate their function. In silico analysis reveals high evolutionary conservation in the DNA-interacting domains for both orthologs, with significant syntenic conservation observed for the prr12b locus. In situ hybridization and RT-qPCR analyses on zebrafish embryos and larvae reveal distinct expression patterns: prr12a is expressed early in zygotic development, mainly in the central nervous system, while prr12b expression initiates during gastrulation, localizing later to dopaminergic telencephalic and diencephalic cell clusters. Both transcripts are enriched in the ganglion cell and inner neural layers of the 72 hpf retina, with prr12b widely distributed in the ciliary marginal zone. In the adult brain, prr12a and prr12b are found in the cerebellum, amygdala and ventral telencephalon, which represent the main areas affected in autistic patients. Overall, this study suggests PRR12's potential involvement in eye and brain development, laying the groundwork for further investigations into PRR12-related neurobehavioral disorders.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

Aspartame Causes Developmental Defects and Teratogenicity in Zebra Fish Embryo: Role of Impaired SIRT1/FOXO3a Axis in Neuron Cells

Aspartame, a widely used artificial sweetener, is present in many food products and beverages worldwide. It has been linked to potential neurotoxicity and developmental defects. However, its teratogenic effect on embryonic development and the underlying potential mechanisms need to be elucidated. We investigated the concentration- and time-dependent effects of aspartame on zebrafish development and teratogenicity. We focused on the role of sirtuin 1 (SIRT1) and Forkhead-box transcription factor (FOXO), two proteins that play key roles in neurodevelopment. It was found that aspartame exposure reduced the formation of larvae and the development of cartilage in zebrafish. It also delayed post-fertilization development by altering the head length and locomotor behavior of zebrafish. RNA-sequencing-based DEG analysis showed that SIRT1 and FOXO3a are involved in neurodevelopment. In silico and in vitro analyses showed that aspartame could target and reduce the expression of SIRT1 and FOXO3a proteins in neuron cells. Additionally, aspartame triggered the reduction of autophagy flux by inhibiting the nuclear translocation of SIRT1 in neuronal cells. The findings suggest that aspartame can cause developmental defects and teratogenicity in zebrafish embryos and reduce autophagy by impairing the SIRT1/FOXO3a axis in neuron cells.

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.

Biallelic loss-of-function variants in CACHD1 cause a novel neurodevelopmental syndrome with facial dysmorphism and multisystem congenital abnormalities

PURPOSE

We established the genetic etiology of a syndromic neurodevelopmental condition characterized by variable cognitive impairment, recognizable facial dysmorphism, and a constellation of extra-neurological manifestations.

METHODS

We performed phenotypic characterization of 6 participants from 4 unrelated families presenting with a neurodevelopmental syndrome and used exome sequencing to investigate the underlying genetic cause. To probe relevance to the neurodevelopmental phenotype and craniofacial dysmorphism, we established two- and three-dimensional human stem cell-derived neural models and generated a stable cachd1 zebrafish mutant on a transgenic cartilage reporter line.

RESULTS

Affected individuals showed mild cognitive impairment, dysmorphism featuring oculo-auriculo abnormalities, and developmental defects involving genitourinary and digestive tracts. Exome sequencing revealed biallelic putative loss-of-function variants in CACHD1 segregating with disease in all pedigrees. RNA sequencing in CACHD1-depleted neural progenitors revealed abnormal expression of genes with key roles in Wnt signaling, neurodevelopment, and organ morphogenesis. CACHD1 depletion in neural progenitors resulted in reduced percentages of post-mitotic neurons and enlargement of 3D neurospheres. Homozygous cachd1 mutant larvae showed mandibular patterning defects mimicking human facial dysmorphism.

CONCLUSION

Our findings support the role of loss-of-function variants in CACHD1 as the cause of a rare neurodevelopmental syndrome with facial dysmorphism and multisystem abnormalities.

Molecular functions of ANKLE2 and its implications in human disease

Ankyrin repeat and LEM domain-containing 2 (ANKLE2) is a scaffolding protein with established roles in cell division and development, the dysfunction of which is increasingly implicated in human disease. ANKLE2 regulates nuclear envelope disassembly at the onset of mitosis and its reassembly after chromosome segregation. ANKLE2 dysfunction is associated with abnormal nuclear morphology and cell division. It regulates the nuclear envelope by mediating protein-protein interactions with barrier to autointegration factor (BANF1; also known as BAF) and with the kinase and phosphatase that modulate the phosphorylation state of BAF. In brain development, ANKLE2 is crucial for proper asymmetric division of neural progenitor cells. In humans, pathogenic loss-of-function mutations in ANKLE2 are associated with primary congenital microcephaly, a condition in which the brain is not properly developed at birth. ANKLE2 is also linked to other disease pathologies, including congenital Zika syndrome, cancer and tauopathy. Here, we review the molecular roles of ANKLE2 and the recent literature on human diseases caused by its dysfunction.

Secreted novel AID/APOBEC-like deaminase 1 (SNAD1) - a new important player in fish immunology

The AID/APOBECs are a group of zinc-dependent cytidine deaminases that catalyse the deamination of bases in nucleic acids, resulting in a cytidine to uridine transition. Secreted novel AID/APOBEC-like deaminases (SNADs), characterized by the presence of a signal peptide are unique among all of intracellular classical AID/APOBECs, which are the central part of antibody diversity and antiviral defense. To date, there is no available knowledge on SNADs including protein characterization, biochemical characteristics and catalytic activity. We used various in silico approaches to define the phylogeny of SNADs, their common structural features, and their potential structural variations in fish species. Our analysis provides strong evidence of the universal presence of SNAD1 proteins/transcripts in fish, in which expression commences after hatching and is highest in anatomical organs linked to the immune system. Moreover, we searched published fish data and identified previously, "uncharacterized proteins" and transcripts as SNAD1 sequences. Our review into immunological research suggests SNAD1 role in immune response to infection or immunization, and interactions with the intestinal microbiota. We also noted SNAD1 association with temperature acclimation, environmental pollution and sex-based expression differences, with females showing higher level. To validate in silico predictions we performed expression studies of several SNAD1 gene variants in carp, which revealed distinct patterns of responses under different conditions. Dual sensitivity to environmental and pathogenic stress highlights its importance in the fish and potentially enhancing thermotolerance and immune defense. Revealing the biological roles of SNADs represents an exciting new area of research related to the role of DNA and/or RNA editing in fish biology.

Error modelled gene expression analysis (EMOGEA) provides a superior overview of time course RNA-seq measurements and low count gene expression

Temporal RNA-sequencing (RNA-seq) studies of bulk samples provide an opportunity for improved understanding of gene regulation during dynamic phenomena such as development, tumor progression or response to an incremental dose of a pharmacotherapeutic. Moreover, single-cell RNA-seq (scRNA-seq) data implicitly exhibit temporal characteristics because gene expression values recapitulate dynamic processes such as cellular transitions. Unfortunately, temporal RNA-seq data continue to be analyzed by methods that ignore this ordinal structure and yield results that are often difficult to interpret. Here, we present Error Modelled Gene Expression Analysis (EMOGEA), a framework for analyzing RNA-seq data that incorporates measurement uncertainty, while introducing a special formulation for those acquired to monitor dynamic phenomena. This method is specifically suited for RNA-seq studies in which low-count transcripts with small-fold changes lead to significant biological effects. Such transcripts include genes involved in signaling and non-coding RNAs that inherently exhibit low levels of expression. Using simulation studies, we show that this framework down-weights samples that exhibit extreme responses such as batch effects allowing them to be modeled with the rest of the samples and maintain the degrees of freedom originally envisioned for a study. Using temporal experimental data, we demonstrate the framework by extracting a cascade of gene expression waves from a well-designed RNA-seq study of zebrafish embryogenesis and an scRNA-seq study of mouse pre-implantation and provide unique biological insights into the regulation of genes in each wave. For non-ordinal measurements, we show that EMOGEA has a much higher rate of true positive calls and a vanishingly small rate of false negative discoveries compared to common approaches. Finally, we provide two packages in Python and R that are self-contained and easy to use, including test data.

Paternal-effect-genes revealed through sperm cryopreservation in Perca fluviatilis

Knowledge about paternal-effect-genes (PEGs) (genes whose expression in the progeny is influenced by paternal factors present in the sperm) in fish is very limited. To explore this issue, we used milt cryopreservation as a specific challenge test for sperm cells, thus enabling selection amidst cryo-sensitivity. We created two groups of Eurasian perch (Perca fluviatilis) as a model - eggs fertilized either with fresh (Fresh group) or cryopreserved (Cryo group) milt from the same male followed by phenotypic-transcriptomic examination of consequences of cryopreservation in obtained progeny (at larval stages). Most of the phenotypical observations were similar in both groups, except the final weight which was higher in the Cryo group. Milt cryopreservation appeared to act as a "positive selection" factor, upregulating most PEGs in the Cryo group. Transcriptomic profile of freshly hatched larvae sourced genes involved in the development of visual perception and we identified them as PEGs. Consequently, larvae from the Cryo group exhibited enhanced eyesight, potentially contributing to more efficient foraging and weight gain compared to the Fresh group. This study unveils, for the first time, the significant influence of the paternal genome on the development of the visual system in fish, highlighting pde6g, opn1lw1, and rbp4l as novel PEGs.

Transcriptome alterations in sf3b4-depleted zebrafish: Insights into cataract formation in retinitis pigmentosa model

Posterior subcapsular cataract (PSC) frequently develops as a complication in patients with retinitis pigmentosa (RP). Despite numerous scientific investigations, the intricate pathomechanisms underlying cataract formation in individuals affected by RP remain elusive. Therefore, our study aims to elucidate the potential pathogenesis of cataracts in an RP model using splicing factor subunit 3b (sf3b4) mutant zebrafish. By analyzing our previously published transcriptome dataset, we identified that, in addition to RP, cataract was listed as the second condition in our transcriptomic analysis. Furthermore, we confirmed the presence of nucleus retention in the lens fiber cells, along with abnormal cytoskeleton expression in both the lens fiber cells and lens epithelial cells in sf3b4-depleted fish. Upon closer examination, we identified 20 differentially expressed genes (DEGs) that played a role in cataract formation, with 95 % of them related to the downregulation of structural lens proteins. Additionally, we also identified that among all the DEGs, 13 % were associated with fibrotic processes. It seems that the significant upregulation of inflammatory mediators, in conjunction with TGF-β signaling, plays a central role in the cellular biology of PSC and posterior capsular opacification (PCO) in sf3b4 mutant fish. In summary, our study provides valuable insights into cataract formation in the RP model of sf3b4 mutants, highlighting its complexity driven by changes in structural lens proteins and increased cytokines/growth factors.

Anp32e protects against accumulation of H2A.Z at Sox motif containing promoters during zebrafish gastrulation

Epigenetic regulation of chromatin states is crucial for proper gene expression programs and progression during development, but precise mechanisms by which epigenetic factors influence differentiation remain poorly understood. Here we find that the histone variant H2A.Z accumulates at Sox motif-containing promoters during zebrafish gastrulation while neighboring genes become transcriptionally active. These changes coincide with reduced expression of anp32e, the H2A.Z histone removal chaperone, suggesting that loss of Anp32e may lead to increases in H2A.Z binding during differentiation. Remarkably, genetic removal of Anp32e in embryos leads to H2A.Z accumulation prior to gastrulation and developmental genes become precociously active. Accordingly, H2A.Z accumulation occurs most extensively at Sox motif-associated genes, including many which are normally activated following gastrulation. Altogether, our results provide compelling evidence for a mechanism in which Anp32e preferentially restricts H2A.Z accumulation at Sox motifs to regulate the initial phases of developmental differentiation in zebrafish.

Optimal Performance Objectives in the Highly Conserved Bone Morphogenetic Protein Signaling Pathway

Throughout development, complex networks of cell signaling pathways drive cellular decision-making across different tissues and contexts. The transforming growth factor β (TGF-β) pathways, including the BMP/Smad pathway, play crucial roles in these cellular responses. However, as the Smad pathway is used reiteratively throughout the life cycle of all animals, its systems-level behavior varies from one context to another, despite the pathway connectivity remaining nearly constant. For instance, some cellular systems require a rapid response, while others require high noise filtering. In this paper, we examine how the BMP- Smad pathway balances trade-offs among three such systems-level behaviors, or "Performance Objectives (POs)": response speed, noise amplification, and the sensitivity of pathway output to receptor input. Using a Smad pathway model fit to human cell data, we show that varying non-conserved parameters (NCPs) such as protein concentrations, the Smad pathway can be tuned to emphasize any of the three POs and that the concentration of nuclear phosphatase has the greatest effect on tuning the POs. However, due to competition among the POs, the pathway cannot simultaneously optimize all three, but at best must balance trade-offs among the POs. We applied the multi-objective optimization concept of the Pareto Front, a widely used concept in economics to identify optimal trade-offs among various requirements. We show that the BMP pathway efficiently balances competing POs across species and is largely Pareto optimal. Our findings reveal that varying the concentration of NCPs allows the Smad signaling pathway to generate a diverse range of POs. This insight identifies how signaling pathways can be optimally tuned for each context.

Metabonomics and Transcriptomics Analyses Reveal the Development Process of the Auditory System in the Embryonic Development Period of the Small Yellow Croaker under Background Noise

Underwater noise pollution has become a potential threat to aquatic animals in the natural environment. The main causes of such pollution are frequent human activities creating underwater environmental noise, including commercial shipping, offshore energy platforms, scientific exploration activities, etc. However, in aquaculture environments, underwater noise pollution has also become an unavoidable problem due to background noise created by aquaculture equipment. Some research has shown that certain fish show adaptability to noise over a period of time. This could be due to fish's special auditory organ, i.e., their "inner ear"; meanwhile, otoliths and sensory hair cells are the important components of the inner ear and are also essential for the function of the auditory system. Recently, research in respect of underwater noise pollution has mainly focused on adult fish, and there is a lack of the research on the effects of underwater noise pollution on the development process of the auditory system in the embryonic development period. Thus, in this study, we collected embryo-larval samples of the small yellow croaker (Larimichthys polyactis) in four important stages of otic vesicle development through artificial breeding. Then, we used metabonomics and transcriptomics analyses to reveal the development process of the auditory system in the embryonic development period under background noise (indoor and underwater environment sound). Finally, we identified 4026 differentially expressed genes (DEGs) and 672 differential metabolites (DMs), including 37 DEGs associated with the auditory system, and many differences mainly existed in the neurula stage (20 h of post-fertilization/20 HPF). We also inferred the regulatory mode and process of some important DEGs (Dnmt1, CPS1, and endothelin-1) in the early development of the auditory system. In conclusion, we suggest that the auditory system development of L. polyactis begins at least in the neurula stage or earlier; the other three stages (tail bud stage, caudal fin fold stage, and heart pulsation stage, 28-35 HPF) mark the rapid development period. We speculate that the effect of underwater noise pollution on the embryo-larval stage probably begins even earlier.

CRK and NCK adaptors may functionally overlap in zebrafish neurodevelopment, as indicated by common binding partners and overlapping expression patterns

CRK adaptor proteins are important for signal transduction mechanisms driving cell proliferation and positioning during vertebrate central nervous system development. Zebrafish lacking both CRK family members exhibit small, disorganized retinas with 50% penetrance. The goal of this study was to determine whether another adaptor protein might functionally compensate for the loss of CRK adaptors. Expression patterns in developing zebrafish, and bioinformatic analyses of the motifs recognized by their SH2 and SH3 domains, suggest NCK adaptors are well-positioned to compensate for loss of CRK adaptors. In support of this hypothesis, proteomic analyses found CRK and NCK adaptors share overlapping interacting partners including known regulators of cell adhesion and migration, suggesting their functional intersection in neurodevelopment.

Deficiency of Acetyltransferase nat10 in Zebrafish Causes Developmental Defects in the Visual Function

Purpose

N4-acetylcytidine (ac4C) is a post-transcriptional RNA modification catalyzed by N-acetyltransferase 10 (NAT10), a critical factor known to influence mRNA stability. However, the role of ac4C in visual development remains unexplored.

Methods

Analysis of public datasets and immunohistochemical staining were conducted to assess the expression pattern of nat10 in zebrafish. We used CRISPR/Cas9 and RNAi technologies to knockout (KO) and knockdown (KD) nat10, the zebrafish ortholog of human NAT10, and evaluated its effects on early development. To assess the impact of nat10 knockdown on visual function, we performed comprehensive histological evaluations and behavioral analyses. Transcriptome profiling and real-time (RT)-PCR were utilized to detect alterations in gene expression resulting from the nat10 knockdown. Dot-blot and RNA immunoprecipitation (RIP)-PCR analyses were conducted to verify changes in ac4C levels in both total RNA and opsin mRNA specifically. Additionally, we used the actinomycin D assay to examine the stability of opsin mRNA following the nat10 KD.

Results

Our study found that the zebrafish NAT10 protein shares similar structural properties with its human counterpart. We observed that the nat10 gene was prominently expressed in the visual system during early zebrafish development. A deficiency of nat10 in zebrafish embryos resulted in increased mortality and developmental abnormalities. Behavioral and histological assessments indicated significant vision impairment in nat10 KD zebrafish. Transcriptomic analysis and RT-PCR identified substantial downregulation of retinal transcripts related to phototransduction, light response, photoreceptors, and visual perception in the nat10 KD group. Dot-blot and RIP-PCR analyses confirmed a pronounced reduction in ac4C levels in both total RNA and specifically in opsin messenger RNA (mRNA). Additionally, by evaluating mRNA decay in zebrafish treated with actinomycin D, we observed a significant decrease in the stability of opsin mRNA in the nat10 KD group.

Conclusions

The ac4C-mediated mRNA modification plays an essential role in maintaining visual development and retinal function. The loss of NAT10-mediated ac4C modification results in significant disruptions to these processes, underlining the importance of this RNA modification in ocular development.

RNA methyltransferase SPOUT1/CENP-32 links mitotic spindle organization with the neurodevelopmental disorder SpADMiSS

SPOUT1/CENP-32 encodes a putative SPOUT RNA methyltransferase previously identified as a mitotic chromosome associated protein. SPOUT1/CENP-32 depletion leads to centrosome detachment from the spindle poles and chromosome misalignment. Aided by gene matching platforms, we identified 24 individuals with neurodevelopmental delays from 18 families with bi-allelic variants in SPOUT1/CENP-32 detected by exome/genome sequencing. Zebrafish spout1/cenp-32 mutants showed reduction in larval head size with concomitant apoptosis likely associated with altered cell cycle progression. In vivo complementation assays in zebrafish indicated that SPOUT1/CENP-32 missense variants identified in humans are pathogenic. Crystal structure analysis of SPOUT1/CENP-32 revealed that most disease-associated missense variants mapped to the catalytic domain. Additionally, SPOUT1/CENP-32 recurrent missense variants had reduced methyltransferase activity in vitro and compromised centrosome tethering to the spindle poles in human cells. Thus, SPOUT1/CENP-32 pathogenic variants cause an autosomal recessive neurodevelopmental disorder: SpADMiSS ( SPOUT1 Associated Development delay Microcephaly Seizures Short stature) underpinned by mitotic spindle organization defects and consequent chromosome segregation errors.

A clinical and biologic review of congenital melanocytic nevi

Congenital melanocytic nevi (CMN) are the result of aberrations in the mitogen-activated protein kinase signal transduction pathway caused by postzygotic somatic mutations. The estimated incidence of newborns with CMN is 1%-2%. The main complications of CMN include proliferative nodules, melanomas, and neurocutaneous melanosis, and the latter two are the most troublesome issues to address. Treatments are primarily taken into account for aesthetic purposes and the reduction of melanoma risk. Due to the much lower incidence of malignant transformation observed in recent studies than in previous data, clinical management paradigms for CMN patients have gradually shifted towards conservative observation and close monitoring. Surgery and lasers are still the main treatments, and targeted therapy may be a promising strategy to help manage complications. With the increase in awareness of mental health, increasing focus has been placed on the quality of life (QoL) and psychological issues of both CMN patients and their parents. Recent studies have revealed that families coping with CMN might endure intense pressure, a major loss in QoL, and psychological problems after diagnosis and during treatment. Here, we sought to present an overview of genetic basis, complications, treatments, and psychological issues related to CMN and hope to provide better management for patients with CMN.

Zebrafish as a model for study of disorders in pyrimidine nucleotide metabolism

Pyrimidine nucleotides are not only the building blocks of DNA and RNA but also participate in multiple cellular metabolic processes, including protein, lipid and polysaccharide biosynthesis. Pyrimidine nucleotides are synthesized by two distinct pathways-the de novo and salvage pathways. Disorders in pyrimidine nucleotide metabolism cause severe neurodegenerative disorders in human. For example, deficiency in thymidylate kinase, an essential enzyme in dTTP synthesis, causes severe microcephaly in human patients. Zebrafish mutants selected by insertion mutagenesis that results in inactive enzymes in pyrimidine metabolism showed also neurological and developmental disorders. In this work I have summarized current data on neurological and developmental disorders caused by defects in enzymes in pyrimidine nucleotide metabolism in zebrafish and compared to human. All these data suggest that zebrafish is a useful animal model to study pathogenic mechanism of neurological disorders due to defect in pyrimidine nucleotide metabolism.

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

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

Translational control of furina by an RNA regulon is important for left-right patterning, heart morphogenesis and cardiac valve function

Heart development is a complex process that requires asymmetric positioning of the heart, cardiac growth and valve morphogenesis. The mechanisms controlling heart morphogenesis and valve formation are not fully understood. The pro-convertase FurinA functions in heart development across vertebrates. How FurinA activity is regulated during heart development is unknown. Through computational analysis of the zebrafish transcriptome, we identified an RNA motif in a variant FurinA transcript harbouring a long 3' untranslated region (3'UTR). The alternative 3'UTR furina isoform is expressed prior to organ positioning. Somatic deletions in the furina 3'UTR lead to embryonic left-right patterning defects. Reporter localisation and RNA-binding assays show that the furina 3'UTR forms complexes with the conserved RNA-binding translational repressor, Ybx1. Conditional ybx1 mutant embryos show premature and increased Furin reporter expression, abnormal cardiac morphogenesis and looping defects. Mutant ybx1 hearts have an expanded atrioventricular canal, abnormal sino-atrial valves and retrograde blood flow from the ventricle to the atrium. This is similar to observations in humans with heart valve regurgitation. Thus, the furina 3'UTR element/Ybx1 regulon is important for translational repression of FurinA and regulation of heart development.

Estrogenic regulation of claudin 5 and tight junction protein 1 gene expression in zebrafish: A role on blood-brain barrier?

The blood-brain barrier (BBB) is a physical interface between the blood and the brain parenchyma, playing key roles in brain homeostasis. In mammals, the BBB is established thanks to tight junctions between cerebral endothelial cells, involving claudin, occludin, and zonula occludens proteins. Estrogens have been documented to modulate BBB permeability. Interestingly, in the brain of zebrafish, the estrogen-synthesizing activity is strong due to the high expression of Aromatase B protein, encoded by the cyp19a1b gene, in radial glial cells (neural stem cells). Given the roles of estrogens in BBB function, we investigated their impact on the expression of genes involved in BBB tight junctions. We treated zebrafish embryos and adult males with 17β-estradiol and observed an increased cerebral expression of tight junction and claudin 5 genes in adult males only. In females, treatment with the nuclear estrogen receptor antagonist (ICI182,780 ) had no impact. Interestingly, telencephalic injuries performed in males decreased tight junction gene expression that was partially reversed with 17β-estradiol. This was further confirmed by extravasation experiments of Evans blue showing that estrogenic treatment limits BBB leakage. We also highlighted the intimate links between endothelial cells and neural stem cells, suggesting that cholesterol and peripheral steroids could be taken up by endothelial cells and used as precursors for estrogen synthesis by neural stem cells. Together, our results show that zebrafish provides an alternative model to further investigate the role of steroids on the expression of genes involved in BBB integrity, both in constitutive and regenerative physiological conditions. The link we described between capillaries endothelial cells and steroidogenic neural cells encourages the use of this model in understanding the mechanisms by which peripheral steroids get into neural tissue and modulate neurogenic activity.

Evaluating the association of biallelic OGDHL variants with significant phenotypic heterogeneity

BACKGROUND

Biallelic variants in OGDHL, encoding part of the α-ketoglutarate dehydrogenase complex, have been associated with highly heterogeneous neurological and neurodevelopmental disorders. However, the validity of this association remains to be confirmed. A second OGDHL patient cohort was recruited to carefully assess the gene-disease relationship.

METHODS

Using an unbiased genotype-first approach, we screened large, multiethnic aggregated sequencing datasets worldwide for biallelic OGDHL variants. We used CRISPR/Cas9 to generate zebrafish knockouts of ogdhl, ogdh paralogs, and dhtkd1 to investigate functional relationships and impact during development. Functional complementation with patient variant transcripts was conducted to systematically assess protein functionality as a readout for pathogenicity.

RESULTS

A cohort of 14 individuals from 12 unrelated families exhibited highly variable clinical phenotypes, with the majority of them presenting at least one additional variant, potentially accounting for a blended phenotype and complicating phenotypic understanding. We also uncovered extreme clinical heterogeneity and high allele frequencies, occasionally incompatible with a fully penetrant recessive disorder. Human cDNA of previously described and new variants were tested in an ogdhl zebrafish knockout model, adding functional evidence for variant reclassification. We disclosed evidence of hypomorphic alleles as well as a loss-of-function variant without deleterious effects in zebrafish variant testing also showing discordant familial segregation, challenging the relationship of OGDHL as a conventional Mendelian gene. Going further, we uncovered evidence for a complex compensatory relationship among OGDH, OGDHL, and DHTKD1 isoenzymes that are associated with neurodevelopmental disorders and exhibit complex transcriptional compensation patterns with partial functional redundancy.

CONCLUSIONS

Based on the results of genetic, clinical, and functional studies, we formed three hypotheses in which to frame observations: biallelic OGDHL variants lead to a highly variable monogenic disorder, variants in OGDHL are following a complex pattern of inheritance, or they may not be causative at all. Our study further highlights the continuing challenges of assessing the validity of reported disease-gene associations and effects of variants identified in these genes. This is particularly more complicated in making genetic diagnoses based on identification of variants in genes presenting a highly heterogenous phenotype such as "OGDHL-related disorders".

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 regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

Highlights

In vivo MPRA systematically interrogates the regulatory potential of endogenous 5' UTRs The 5' UTR alone is sufficient to regulate the dynamics of ribosome recruitment during early embryogenesis The MPRA identifies 5' UTR cis -regulatory motifs for translation initiation control 5' UTR length, upstream AUGs and motif grammar contribute to the differential regulatory capability of 5' UTR switching isoforms.

Competition for H2A.Z underlies the developmental impacts of repetitive element de-repression

The histone variant H2A.Z is central to early embryonic development, determining transcriptional competency through chromatin regulation of gene promoters and enhancers. In addition to genic loci, we find that H2A.Z resides at a subset of evolutionarily young repetitive elements, including DNA transposons, long interspersed nuclear elements and long terminal repeats, during early zebrafish development. Moreover, increases in H2A.Z occur when repetitive elements become transcriptionally active. Acquisition of H2A.Z corresponds with a reduction in the levels of the repressive histone modification H3K9me3 and a moderate increase in chromatin accessibility. Notably, however, de-repression of repetitive elements also leads to a significant reduction in H2A.Z over non-repetitive genic loci. Genic loss of H2A.Z is accompanied by transcriptional silencing at adjacent coding sequences, but remarkably, these impacts are mitigated by augmentation of total H2A.Z protein via transgenic overexpression. Our study reveals that levels of H2A.Z protein determine embryonic sensitivity to de-repression of repetitive elements, that repetitive elements can function as a nuclear sink for epigenetic factors and that competition for H2A.Z greatly influences overall transcriptional output during development. These findings uncover general mechanisms in which counteractive biological processes underlie phenotypic outcomes.

Gene expression in notochord and nuclei pulposi: a study of gene families across the chordate phylum

The transition from notochord to vertebral column is a crucial milestone in chordate evolution and in prenatal development of all vertebrates. As ossification of the vertebral bodies proceeds, involutions of residual notochord cells into the intervertebral discs form the nuclei pulposi, shock-absorbing structures that confer flexibility to the spine. Numerous studies have outlined the developmental and evolutionary relationship between notochord and nuclei pulposi. However, the knowledge of the similarities and differences in the genetic repertoires of these two structures remains limited, also because comparative studies of notochord and nuclei pulposi across chordates are complicated by the gene/genome duplication events that led to extant vertebrates. Here we show the results of a pilot study aimed at bridging the information on these two structures. We have followed in different vertebrates the evolutionary trajectory of notochord genes identified in the invertebrate chordate Ciona, and we have evaluated the extent of conservation of their expression in notochord cells. Our results have uncovered evolutionarily conserved markers of both notochord development and aging/degeneration of the nuclei pulposi.

Cereblon influences the timing of muscle differentiation in Ciona tadpoles

Thalidomide has a dark history as a teratogen, but in recent years, its derivates have been shown to function as potent chemotherapeutic agents. These drugs bind cereblon (CRBN), the substrate receptor of an E3 ubiquitin ligase complex, and modify its degradation targets. Despite these insights, remarkably little is known about the normal function of cereblon in development. Here, we employ Ciona, a simple invertebrate chordate, to identify endogenous Crbn targets. In Ciona, Crbn is specifically expressed in developing muscles during tail elongation before they acquire contractile activity. Crbn expression is activated by Mrf, the ortholog of MYOD1, a transcription factor important for muscle differentiation. CRISPR/Cas9-mediated mutations of Crbn lead to precocious onset of muscle contractions. By contrast, overexpression of Crbn delays contractions and is associated with decreased expression of contractile protein genes such as troponin. This reduction is possibly due to reduced Mrf protein levels without altering Mrf mRNA levels. Our findings suggest that Mrf and Crbn form a negative feedback loop to control the precision of muscle differentiation during tail elongation.

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.

MCPIP1 functions as a safeguard of early embryonic development

Monocyte chemoattractant protein-induced protein 1 (MCPIP1), also called Regnase-1, is an RNase that has been described as a key negative modulator of inflammation. MCPIP1 also controls numerous tumor-related processes, such as proliferation, apoptosis and differentiation. In this study, we utilized a zebrafish model to investigate the role of Mcpip1 during embryogenic development. Our results demonstrated that during embryogenesis, the expression of the zc3h12a gene encoding Mcpip1 undergoes dynamic changes. Its transcript levels gradually increase from the 2-cell stage to the spherical stage and then decrease rapidly. We further found that ectopic overexpression of wild-type Mcpip1 but not the catalytically inactive mutant form resulted in an embryonic lethal phenotype in zebrafish embryos (24 hpf). At the molecular level, transcriptomic profiling revealed extensive changes in the expression of genes encoding proteins important in the endoplasmic reticulum stress response and in protein folding as well as involved in the formation of primary germ layer, mesendoderm and endoderm development, heart morphogenesis and cell migration. Altogether, our results demonstrate that the expression of zc3h12a must be tightly controlled during the first cell divisions of zebrafish embryos and that a rapid decrease in its mRNA expression is an important factor promoting proper embryo development.

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

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

Systematic Proteomics Study on the Embryonic Development of Danio rerio

The early development of zebrafish (Danio rerio) is a complex and dynamic physiological process involving cell division, differentiation, and movement. Currently, the genome and transcriptome techniques have been widely used to study the embryonic development of zebrafish. However, the research of proteomics based on proteins that directly execute functions is relatively vacant. In this work, we apply label-free quantitative proteomics to explore protein profiling during zebrafish's embryogenesis, and a total of 5961 proteins were identified at 10 stages of zebrafish's early development. The identified proteins were divided into 11 modules according to weighted gene coexpression network analysis (WGCNA), and the characteristics between modules were significantly different. For example, mitochondria-related functions enriched the early development of zebrafish. Primordial germ cell-related proteins were identified at the 4-cell stage, while the eye development event is dominated at 5 days post fertilization (dpf). By combining with published transcriptomics data, we discovered some proteins that may be involved in activating zygotic genes. Meanwhile, 137 novel proteins were identified. This study comprehensively analyzed the dynamic processes in the embryonic development of zebrafish from the perspective of proteomics. It provided solid data support for further understanding of the molecular mechanism of its development.

Reduced GS Domain Serine/Threonine Requirements of Fibrodysplasia Ossificans Progressiva Mutant Type I BMP Receptor ACVR1 in the Zebrafish

Fibrodysplasia ossificans progressiva (FOP) is a rare human genetic condition characterized by altered skeletal development and extraskeletal bone formation. All cases of FOP are caused by mutations in the type I bone morphogenetic protein (BMP) receptor gene ACVR1 that result in overactivation of the BMP signaling pathway. Activation of the wild-type ACVR1 kinase requires assembly of a tetrameric type I and II BMP receptor complex followed by phosphorylation of the ACVR1 GS domain by type II BMP receptors. Previous studies showed that the FOP-mutant ACVR1-R206H required type II BMP receptors and presumptive glycine/serine-rich (GS) domain phosphorylation for overactive signaling. Structural modeling of the ACVR1-R206H mutant kinase domain supports the idea that FOP mutations alter the conformation of the GS domain, but it is unclear how this leads to overactive signaling. Here we show, using a developing zebrafish embryo BMP signaling assay, that the FOP-mutant receptors ACVR1-R206H and -G328R have reduced requirements for GS domain phosphorylatable sites to signal compared to wild-type ACVR1. Further, ligand-independent and ligand-dependent signaling through the FOP-mutant ACVR1 receptors have distinct GS domain phosphorylatable site requirements. ACVR1-G328R showed increased GS domain serine/threonine requirements for ligand-independent signaling compared to ACVR1-R206H, whereas it exhibited reduced serine/threonine requirements for ligand-dependent signaling. Remarkably, while ACVR1-R206H does not require the type I BMP receptor partner, Bmpr1, to signal, a ligand-dependent GS domain mutant of ACVR1-R206H could signal independently of Bmpr1 only when Bmp7 ligand was overexpressed. Of note, unlike human ACVR1-R206H, the zebrafish paralog Acvr1l-R203H does not show increased signaling activity. However, in domain-swapping studies, the human kinase domain, but not the human GS domain, was sufficient to confer overactive signaling to the Acvr1l-R203H receptor. Together these results reflect the importance of GS domain activation and kinase domain functions in regulating ACVR1 signaling and identify mechanisms of reduced regulatory constraints conferred by FOP mutations. © 2023 American Society for Bone and Mineral Research (ASBMR).

Downregulation of Zebrafish Cytosolic Sialidase Neu3.2 Affects Skeletal Muscle Development

Sialidases remove terminal sialic acids residues from the non-reducing ends of glycoconjugates. They have been recognized as catabolic enzymes that work within different subcellular compartments and can ensure the proper turn-over of glycoconjugates. Four mammalian sialidases (NEU1-4) exist, with different subcellular localization, pH optimum and substrate specificity. In zebrafish, seven different sialidases, with high homology to mammalian counterparts, have been identified. Zebrafish Neu3.2 is similar to the human cytosolic sialidase NEU2, which is involved in skeletal muscle differentiation and exhibits a broad substrate specificity toward gangliosides and glycoproteins. In zebrafish neu3.2, mRNA is expressed during somite development, and its enzymatic activity has been detected in the skeletal muscle and heart of adult animals. In this paper, 1-4-cell-stage embryos injected with neu3.2 splice-blocking morpholino showed severe embryonic defects, mainly in somites, heart and anterior-posterior axis formation. Myog and myod1 expressions were altered in morphants, and impaired musculature formation was associated with a defective locomotor behavior. Finally, the co-injection of Neu2 mouse mRNA in morphants rescued the phenotype. These data are consistent with the involvement of cytosolic sialidase in pathologies related to muscle formation and support the validity of the model to investigate the pathogenesis of the diseases.

Betaglycan promoter activity is differentially regulated during myogenesis in zebrafish embryo somites

BACKGROUND

Betaglycan, also known as the TGFβ type III receptor (Tgfbr3), is a co-receptor that modulates TGFβ family signaling. Tgfbr3 is upregulated during C2C12 myoblast differentiation and expressed in mouse embryos myocytes.

RESULTS

To investigate tgfbr3 transcriptional regulation during zebrafish embryonic myogenesis, we cloned a 3.2 kb promoter fragment that drives reporter transcription during C2C12 myoblasts differentiation and in the Tg(tgfbr3:mCherry) transgenic zebrafish. We detect tgfbr3 protein and mCherry expression in the adaxial cells concomitantly with the onset of their radial migration to become slow-twitch muscle fibers in the Tg(tgfbr3:mCherry). Remarkably, this expression displays a measurable antero-posterior somitic gradient expression.

CONCLUSIONS

tgfbr3 is transcriptionally regulated during somitic muscle development in zebrafish with an antero-posterior gradient expression that preferentially marks the adaxial cells and their descendants.

Effect of caldesmon mutations in the development of zebrafish embryos

Cancer is a profound medical concern and better treatments are needed for cancer patients. Therefore, new cancer targets are constantly being studied. These targets need not only be relevant for cancer progression, but their modulation needs to be tolerated reasonably well by the host. Caldesmon is one of these proposed novel targets for cancer therapy. Therefore, we analyzed effects of caldesmon mutations in normal development using genetically modified zebrafish embryos. We analyzed mutations in both zebrafish caldesmon genes, cald1a and cald1b and analyzed effects of either mutation alone or as in combination in double homozygous embryos using molecular, morphological and functional analyses. The effects of caldesmon mutations were mild and the gross development of zebrafish embryos was normal. The caldesmon mutant embryos had, however, alterations in response to light-stimulus in behavioural assays. Taken together, the effects of caldesmon mutations in the development of zebrafish embryos were reasonably well tolerated and did not indicate significant concerns for caldesmon being a potential target for cancer therapy.

Transposable elements drive the evolution of metazoan zinc finger genes

Cys2-His2 zinc finger genes (ZNFs) form the largest family of transcription factors in metazoans. ZNF evolution is highly dynamic and characterized by the rapid expansion and contraction of numerous subfamilies across the animal phylogeny. The forces and mechanisms underlying rapid ZNF evolution remain poorly understood, but there is growing evidence that, in tetrapods, the targeting and repression of lineage-specific transposable elements (TEs) plays a critical role in the evolution of the Krüppel-associated box ZNF (KZNF) subfamily. Currently, it is unknown whether this function and coevolutionary relationship is unique to KZNFs or is a broader feature of metazoan ZNFs. Here, we present evidence that genomic conflict with TEs has been a central driver of the diversification of ZNFs in animals. Sampling from 3221 genome assemblies, we show that the copy number of retroelements correlates with that of ZNFs across at least 750 million years of metazoan evolution. Using computational predictions, we show that ZNFs preferentially bind TEs in diverse animal species. We further investigate the largest ZNF subfamily found in cyprinid fish, which is characterized by a conserved sequence we dubbed the fish N-terminal zinc finger-associated (FiNZ) domain. Zebrafish possess approximately 700 FiNZ-ZNFs, many of which are evolving adaptively under positive selection. Like mammalian KZNFs, most zebrafish FiNZ-ZNFs are expressed at the onset of zygotic genome activation, and blocking their translation using morpholinos during early embryogenesis results in derepression of transcriptionally active TEs. Together, these data suggest that ZNF diversification has been intimately connected to TE expansion throughout animal evolution.

Molecular evolution of Keap1 was essential for adaptation of vertebrates to terrestrial life

Reactive oxygen species (ROS) posed a risk for the transition of vertebrates from aquatic to terrestrial life. How ancestral organisms adapted to such ROS exposure has remained a mystery. Here, we show that attenuation of the activity of the ubiquitin ligase CRL3Keap1 for the transcription factor Nrf2 during evolution was key to development of an efficient response to ROS exposure. The Keap1 gene was duplicated in fish to give rise to Keap1A and the only remaining mammalian paralog Keap1B, the latter of which shows a lower affinity for Cul3 and contributes to robust Nrf2 induction in response to ROS exposure. Mutation of mammalian Keap1 to resemble zebrafish Keap1A resulted in an attenuated Nrf2 response, and most knock-in mice expressing such a Keap1 mutant died on exposure as neonates to sunlight-level ultraviolet radiation. Our results suggest that molecular evolution of Keap1 was essential for adaptation to terrestrial life.