An improved zebrafish transcriptome annotation for sensitive and comprehensive detection of cell type-specific genes

Single-cell Transcriptional Landscape of Temporal Neutrophil Response to Burn Wound in Larval Zebrafish

Neutrophils accumulate early in tissue injury. However, the cellular and functional heterogeneity of neutrophils during homeostasis and in response to tissue damage remains unclear. In this study, we use larval zebrafish to understand neutrophil responses to thermal injury. Single-cell transcriptional mapping of myeloid cells during a 3-d time course in burn and control larvae revealed distinct neutrophil subsets and their cell-cell interactions with macrophages across time and conditions. The trajectory formed by three zebrafish neutrophil subsets resembles human neutrophil maturation, with varying transition patterns between conditions. Through ligand-receptor cell-cell interaction analysis, we found that neutrophils communicate more in burns in a pathway and temporal manner. Finally, we identified the correlation between zebrafish myeloid signatures and human burn severity, establishing GPR84+ neutrophils as a potential marker of early innate immune response in burns. This work builds a comparative single-cell transcriptomic framework to identify neutrophil markers of tissue damage using model organisms.

Single-cell analysis of innate spinal cord regeneration identifies intersecting modes of neuronal repair

Adult zebrafish have an innate ability to recover from severe spinal cord injury. Here, we report a comprehensive single nuclear RNA sequencing atlas that spans 6 weeks of regeneration. We identify cooperative roles for adult neurogenesis and neuronal plasticity during spinal cord repair. Neurogenesis of glutamatergic and GABAergic neurons restores the excitatory/inhibitory balance after injury. In addition, a transient population of injury-responsive neurons (iNeurons) show elevated plasticity 1 week post-injury. We found iNeurons are injury-surviving neurons that acquire a neuroblast-like gene expression signature after injury. CRISPR/Cas9 mutagenesis showed iNeurons are required for functional recovery and employ vesicular trafficking as an essential mechanism that underlies neuronal plasticity. This study provides a comprehensive resource of the cells and mechanisms that direct spinal cord regeneration and establishes zebrafish as a model of plasticity-driven neural repair.

Cohesin composition and dosage independently affect early development in zebrafish

Cohesin, a chromatin-associated protein complex with four core subunits (Smc1a, Smc3, Rad21 and either Stag1 or 2), has a central role in cell proliferation and gene expression in metazoans. Human developmental disorders termed 'cohesinopathies' are characterized by germline variants of cohesin or its regulators that do not entirely eliminate cohesin function. However, it is not clear whether mutations in individual cohesin subunits have independent developmental consequences. Here, we show that zebrafish rad21 or stag2b mutants independently influence embryonic tailbud development. Both mutants have altered mesoderm induction, but only homozygous or heterozygous rad21 mutation affects cell cycle gene expression. stag2b mutants have narrower notochords and reduced Wnt signaling in neuromesodermal progenitors as revealed by single-cell RNA sequencing. Stimulation of Wnt signaling rescues transcription and morphology in stag2b, but not rad21, mutants. Our results suggest that mutations altering the quantity versus composition of cohesin have independent developmental consequences, with implications for the understanding and management of cohesinopathies.

The axillary lymphoid organ - an external, experimentally accessible immune organ in the zebrafish

Lymph nodes and other secondary lymphoid organs play critical roles in immune surveillance and immune activation in mammals, but the deep internal locations of these organs make it challenging to image and study them in living animals. Here, we describe a previously uncharacterized external immune organ in the zebrafish ideally suited for studying immune cell dynamics in vivo , the axillary lymphoid organ (ALO). This small, translucent organ has an outer cortex teeming with immune cells, an inner medulla with a mesh-like network of fibroblastic reticular cells along which immune cells migrate, and a network of lymphatic vessels draining to a large adjacent lymph sac. Noninvasive high-resolution imaging of transgenically marked immune cells can be carried out in the lobes of living animals, and the ALO is readily accessible to external treatment. This newly discovered tissue provides a superb model for dynamic live imaging of immune cells and their interaction with pathogens and surrounding tissues, including blood and lymphatic vessels.

Teaser

A newly characterized external zebrafish lymphoid organ provides a powerful model for live imaging of immune cell dynamics.

RNaseH-based ribodepletion of total planarian RNA improves detection of longer and non-polyadenylated transcripts

The overwhelming majority of RNA species isolated from cells or tissues using organic extraction are ribosomal RNAs (rRNA), whereas a relatively small percentage are messenger RNAs (mRNA). For studies that seek to detect mRNA transcripts and measure changes in their expression, this lopsided ratio of desired transcripts to undesired transcripts creates a significant challenge to obtaining sensitive and reproducible results. One method for improving mRNA detection is to selectively amplify polyadenylated (polyA) mRNA molecules when generating RNA-seq libraries, a strategy that is generally very successful in many species. However, this strategy is less effective when starting with total RNA from some species e.g., the planarian species Schmidtea mediterranea (S.med), as it generates libraries that still contain significant and variable amounts of rRNA reads. Further, commercially available ribodepletion kits do not efficiently deplete rRNAs from these samples because their sequences are divergent from mammalian rRNAs. Here we report a customized, optimized, and economical ribodepletion strategy than allows the generation of comprehensive RNA-seq libraries with less than one percent rRNA contamination. We show that this method improves transcript detection, particularly for those without polyA tails (e.g., core histones) and those that are relatively long (e.g., microtubule motor proteins). Using this custom ribodepletion approach, we also detected many transcripts that are not represented in the most recent set of S.med gene annotations, including a subset that are likely expressed transposable elements (TEs). To facilitate future differential expression analyses of these newly identified loci, we created both an annotation file of the new loci we identified and a bioinformatic pipeline for generating additional annotations from future libraries. As significant recent research shows that TE activation is regulated and functionally important, the resources provided here will provide a starting point for investigating such mechanisms in planarians and other species with less conserved rRNA sequences.

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.

Elevated phagocytic capacity directs innate spinal cord repair

Immune cells elicit a continuum of transcriptional and functional states after spinal cord injury (SCI). In mammals, inefficient debris clearance and chronic inflammation impede recovery and overshadow pro-regenerative immune functions. We found that, unlike mammals, zebrafish SCI elicits transient immune activation and efficient debris clearance, without causing chronic inflammation. Single-cell transcriptomics and inducible genetic ablation showed zebrafish macrophages are highly phagocytic and required for regeneration. Cross-species comparisons between zebrafish and mammalian macrophages identified transcription and immune response regulator ( tcim ) as a macrophage-enriched zebrafish gene. Genetic deletion of zebrafish tcim impairs phagocytosis and regeneration, causes aberrant and chronic immune activation, and can be rescued by transplanting wild-type immune precursors into tcim mutants. Conversely, genetic expression of human TCIM accelerates debris clearance and regeneration by reprogramming myeloid precursors into activated phagocytes. This study establishes a central requirement for elevated phagocytic capacity to achieve innate spinal cord repair.

Reconstruction of macroglia and adult neurogenesis evolution through cross-species single-cell transcriptomic analyses

Macroglia fulfill essential functions in the adult vertebrate brain, producing and maintaining neurons and regulating neuronal communication. However, we still know little about their emergence and diversification. We used the zebrafish D. rerio as a distant vertebrate model with moderate glial diversity as anchor to reanalyze datasets covering over 600 million years of evolution. We identify core features of adult neurogenesis and innovations in the mammalian lineage with a potential link to the rarity of radial glia-like cells in adult humans. Our results also suggest that functions associated with astrocytes originated in a multifunctional cell type fulfilling both neural stem cell and astrocytic functions before these diverged. Finally, we identify conserved elements of macroglial cell identity and function and their time of emergence during evolution.

Distinct features of the regenerating heart uncovered through comparative single-cell profiling

Adult humans respond to heart injury by forming a permanent scar, yet other vertebrates are capable of robust and complete cardiac regeneration. Despite progress towards characterizing the mechanisms of cardiac regeneration in fish and amphibians, the large evolutionary gulf between mammals and regenerating vertebrates complicates deciphering which cellular and molecular features truly enable regeneration. To better define these features, we compared cardiac injury responses in zebrafish and medaka, two fish species that share similar heart anatomy and common teleost ancestry but differ in regenerative capability. We used single-cell transcriptional profiling to create a time-resolved comparative cell atlas of injury responses in all major cardiac cell types across both species. With this approach, we identified several key features that distinguish cardiac injury response in the non-regenerating medaka heart. By comparing immune responses to injury, we found altered cell recruitment and a distinct pro-inflammatory gene program in medaka leukocytes, and an absence of the injury-induced interferon response seen in zebrafish. In addition, we found a lack of pro-regenerative signals, including nrg1 and retinoic acid, from medaka endothelial and epicardial cells. Finally, we identified alterations in the myocardial structure in medaka, where they lack primordial layer cardiomyocytes and fail to employ a cardioprotective gene program shared by regenerating vertebrates. Our findings reveal notable variation in injury response across nearly all major cardiac cell types in zebrafish and medaka, demonstrating how evolutionary divergence influences the hidden cellular features underpinning regenerative potential in these seemingly similar vertebrates.

Single-cell transcriptional landscape of temporal neutrophil response to burn wound in larval zebrafish

Neutrophils accumulate early in tissue injury. However, the cellular and functional heterogeneity of neutrophils during homeostasis and in response to tissue damage remains unclear. Here, we use larval zebrafish to understand neutrophil responses to thermal injury. Single-cell transcriptional mapping of myeloid cells during a 3-day time course in burn and control larvae revealed distinct neutrophil subsets and their cell-cell interactions with macrophages across time and conditions. The trajectory formed by three zebrafish neutrophil subsets resembles human neutrophil maturation, with varying transition patterns between conditions. Through ligand-receptor cell-cell interaction analysis, we found neutrophils communicate more in burns in a pathway and temporal manner. Finally, we identified the correlation between zebrafish myeloid signatures and human burn severity, establishing GPR84+ neutrophils as a potential marker of early innate immune response in burns. This work builds the molecular foundation and a comparative single-cell genomic framework to identify neutrophil markers of tissue damage using model organisms.

(Zebra)fishing for nephrogenesis genes

Kidney disease is a devastating condition affecting millions of people worldwide, where over 100,000 patients in the United States alone remain waiting for a lifesaving organ transplant. Concomitant with a surge in personalized medicine, single-gene mutations, and polygenic risk alleles have been brought to the forefront as core causes of a spectrum of renal disorders. With the increasing prevalence of kidney disease, it is imperative to make substantial strides in the field of kidney genetics. Nephrons, the core functional units of the kidney, are epithelial tubules that act as gatekeepers of body homeostasis by absorbing and secreting ions, water, and small molecules to filter the blood. Each nephron contains a series of proximal and distal segments with explicit metabolic functions. The embryonic zebrafish provides an ideal platform to systematically dissect the genetic cues governing kidney development. Here, we review the use of zebrafish to discover nephrogenesis genes.

The development of brain pericytes requires expression of the transcription factor nkx3.1 in intermediate precursors

Brain pericytes are one of the critical cell types that regulate endothelial barrier function and activity, thus ensuring adequate blood flow to the brain. The genetic pathways guiding undifferentiated cells into mature pericytes are not well understood. We show here that pericyte precursor populations from both neural crest and head mesoderm of zebrafish express the transcription factor nkx3.1 develop into brain pericytes. We identify the gene signature of these precursors and show that an nkx3.1-, foxf2a-, and cxcl12b-expressing pericyte precursor population is present around the basilar artery prior to artery formation and pericyte recruitment. The precursors later spread throughout the brain and differentiate to express canonical pericyte markers. Cxcl12b-Cxcr4 signaling is required for pericyte attachment and differentiation. Further, both nkx3.1 and cxcl12b are necessary and sufficient in regulating pericyte number as loss inhibits and gain increases pericyte number. Through genetic experiments, we have defined a precursor population for brain pericytes and identified genes critical for their differentiation.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Cephalopod-omics: Emerging Fields and Technologies in Cephalopod Biology

Few animal groups can claim the level of wonder that cephalopods instill in the minds of researchers and the general public. Much of cephalopod biology, however, remains unexplored: the largest invertebrate brain, difficult husbandry conditions, and complex (meta-)genomes, among many other things, have hindered progress in addressing key questions. However, recent technological advancements in sequencing, imaging, and genetic manipulation have opened new avenues for exploring the biology of these extraordinary animals. The cephalopod molecular biology community is thus experiencing a large influx of researchers, emerging from different fields, accelerating the pace of research in this clade. In the first post-pandemic event at the Cephalopod International Advisory Council (CIAC) conference in April 2022, over 40 participants from all over the world met and discussed key challenges and perspectives for current cephalopod molecular biology and evolution. Our particular focus was on the fields of comparative and regulatory genomics, gene manipulation, single-cell transcriptomics, metagenomics, and microbial interactions. This article is a result of this joint effort, summarizing the latest insights from these emerging fields, their bottlenecks, and potential solutions. The article highlights the interdisciplinary nature of the cephalopod-omics community and provides an emphasis on continuous consolidation of efforts and collaboration in this rapidly evolving field.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Multi-omics analysis of zebrafish response to tick saliva reveals biological processes associated with alpha-Gal syndrome

The alpha-Gal syndrome (AGS) is a tick-borne allergy. A multi-omics approach was used to determine the effect of tick saliva and mammalian meat consumption on zebrafish gut transcriptome and proteome. Bioinformatics analysis using R software was focused on significant biological and metabolic pathway changes associated with AGS. Ortholog mapping identified highly concordant human ortholog genes for the detection of disease-enriched pathways. Tick saliva treatment increased zebrafish mortality, incidence of hemorrhagic type allergic reactions and changes in behavior and feeding patterns. Transcriptomics analysis showed downregulation of biological and metabolic pathways correlated with anti-alpha-Gal IgE and allergic reactions to tick saliva affecting blood circulation, cardiac and vascular smooth muscle contraction, behavior and sensory perception. Disease enrichment analysis revealed downregulated orthologous genes associated with human disorders affecting nervous, musculoskeletal, and cardiovascular systems. Proteomics analysis revealed suppression of pathways associated with immune system production of reactive oxygen species and cardiac muscle contraction. Underrepresented proteins were mainly linked to nervous and metabolic human disorders. Multi-omics data revealed inhibition of pathways associated with adrenergic signaling in cardiomyocytes, and heart and muscle contraction. Results identify tick saliva-related biological pathways supporting multisystemic organ involvement and linking α-Gal sensitization with other illnesses for the identification of potential disease biomarkers.

Warming affects routine swimming activity and novel odour response in larval zebrafish

Temperature is a primary factor affecting the physiology of ectothermic animals and global warming of water bodies may therefore impact aquatic life. Understanding the effects of near-future predicted temperature changes on the behaviour and underlying molecular mechanisms of aquatic animals is of particular importance, since behaviour mediates survival. In this study, we investigate the effects of developmental temperature on locomotory behaviour and olfactory learning in the zebrafish, Danio rerio. We exposed zebrafish from embryonic stage to either control (28 °C) or elevated temperature (30 °C) for seven days. Overall, warming reduced routine swimming activity and caused upregulation of metabolism and neuron development genes. When exposed to olfactory cues, namely catfish cue, a non-alarming but novel odour, and conspecifics alarming cue, warming differently affected the larvae response to the two cues. An increase in locomotory activity and a large transcriptional reprogramming was observed at elevated temperature in response to novel odour, with upregulation of cell signalling, neuron development and neuron functioning genes. As this response was coupled with the downregulation of genes involved in protein translation and ATP metabolism, novel odour recognition in future-predicted thermal conditions would require energetic trade-offs between expensive baseline processes and responsive functions. To evaluate their learning abilities at both temperatures, larvae were conditioned with a mixture of conspecifics alarm cue and catfish cue. Regardless of temperature, no behavioural nor gene expression changes were detected, reinforcing our findings that warming mainly affects zebrafish molecular response to novel odours. Overall, our results show that future thermal conditions will likely impact developing stages, causing trade-offs following novel olfactory detection in the environment.

Molecular analyses of zebrafish V0v spinal interneurons and identification of transcriptional regulators downstream of Evx1 and Evx2 in these cells

BACKGROUND

V0v spinal interneurons are highly conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells.

METHODS

To identify candidate members of V0v gene regulatory networks, we FAC-sorted wild-type and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

RESULTS

Our data reveal two molecularly distinct subtypes of zebrafish V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuron expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression.

CONCLUSIONS

This study identifies two molecularly distinct subsets of zebrafish V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of these cells. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

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.

Single-cell RNAseq analysis of spinal locomotor circuitry in larval zebrafish

Identification of the neuronal types that form the specialized circuits controlling distinct behaviors has benefited greatly from the simplicity offered by zebrafish. Electrophysiological studies have shown that in addition to connectivity, understanding of circuitry requires identification of functional specializations among individual circuit components, such as those that regulate levels of transmitter release and neuronal excitability. In this study, we use single-cell RNA sequencing (scRNAseq) to identify the molecular bases for functional distinctions between motoneuron types that are causal to their differential roles in swimming. The primary motoneuron, in particular, expresses high levels of a unique combination of voltage-dependent ion channel types and synaptic proteins termed functional 'cassettes.' The ion channel types are specialized for promoting high-frequency firing of action potentials and augmented transmitter release at the neuromuscular junction, both contributing to greater power generation. Our transcriptional profiling of spinal neurons further assigns expression of this cassette to specific interneuron types also involved in the central circuitry controlling high-speed swimming and escape behaviors. Our analysis highlights the utility of scRNAseq in functional characterization of neuronal circuitry, in addition to providing a gene expression resource for studying cell type diversity.

Single-cell analysis reveals distinct fibroblast plasticity during tenocyte regeneration in zebrafish

Despite their importance in tissue maintenance and repair, fibroblast diversity and plasticity remain poorly understood. Using single-cell RNA sequencing, we uncover distinct sclerotome-derived fibroblast populations in zebrafish, including progenitor-like perivascular/interstitial fibroblasts, and specialized fibroblasts such as tenocytes. To determine fibroblast plasticity in vivo, we develop a laser-induced tendon ablation and regeneration model. Lineage tracing reveals that laser-ablated tenocytes are quickly regenerated by preexisting fibroblasts. By combining single-cell clonal analysis and live imaging, we demonstrate that perivascular/interstitial fibroblasts actively migrate to the injury site, where they proliferate and give rise to new tenocytes. By contrast, perivascular fibroblast-derived pericytes or specialized fibroblasts, including tenocytes, exhibit no regenerative plasticity. Active Hedgehog (Hh) signaling is required for the proliferation of activated fibroblasts to ensure efficient tenocyte regeneration. Together, our work highlights the functional diversity of fibroblasts and establishes perivascular/interstitial fibroblasts as tenocyte progenitors that promote tendon regeneration in a Hh signaling-dependent manner.

Electrical synapse structure requires distinct isoforms of a postsynaptic scaffold

Electrical synapses are neuronal gap junction (GJ) channels associated with a macromolecular complex called the electrical synapse density (ESD), which regulates development and dynamically modifies electrical transmission. However, the proteomic makeup and molecular mechanisms utilized by the ESD that direct electrical synapse formation are not well understood. Using the Mauthner cell of zebrafish as a model, we previously found that the intracellular scaffolding protein ZO1b is a member of the ESD, localizing postsynaptically, where it is required for GJ channel localization, electrical communication, neural network function, and behavior. Here, we show that the complexity of the ESD is further diversified by the genomic structure of the ZO1b gene locus. The ZO1b gene is alternatively initiated at three transcriptional start sites resulting in isoforms with unique N-termini that we call ZO1b-Alpha, -Beta, and -Gamma. We demonstrate that ZO1b-Beta and ZO1b-Gamma are broadly expressed throughout the nervous system and localize to electrical synapses. By contrast, ZO1b-Alpha is expressed mainly non-neuronally and is not found at synapses. We generate mutants in all individual isoforms, as well as double mutant combinations in cis on individual chromosomes, and find that ZO1b-Beta is necessary and sufficient for robust GJ channel localization. ZO1b-Gamma, despite its localization to the synapse, plays an auxiliary role in channel localization. This study expands the notion of molecular complexity at the ESD, revealing that an individual genomic locus can contribute distinct isoforms to the macromolecular complex at electrical synapses. Further, independent scaffold isoforms have differential contributions to developmental assembly of the interneuronal GJ channels. We propose that ESD molecular complexity arises both from the diversity of unique genes and from distinct isoforms encoded by single genes. Overall, ESD proteomic diversity is expected to have critical impacts on the development, structure, function, and plasticity of electrical transmission.

Rtf1-dependent transcriptional pausing regulates cardiogenesis

During heart development, a well-characterized network of transcription factors initiates cardiac gene expression and defines the precise timing and location of cardiac progenitor specification. However, our understanding of the post-initiation transcriptional events that regulate cardiac gene expression is still incomplete. The PAF1C component Rtf1 is a transcription regulatory protein that modulates pausing and elongation of RNA Pol II, as well as cotranscriptional histone modifications. Here we report that Rtf1 is essential for cardiogenesis in fish and mammals, and that in the absence of Rtf1 activity, cardiac progenitors arrest in an immature state. We found that Rtf1's Plus3 domain, which confers interaction with the transcriptional pausing and elongation regulator Spt5, was necessary for cardiac progenitor formation. ChIP-seq analysis further revealed changes in the occupancy of RNA Pol II around the transcription start site (TSS) of cardiac genes in rtf1 morphants reflecting a reduction in transcriptional pausing. Intriguingly, inhibition of pause release in rtf1 morphants and mutants restored the formation of cardiac cells and improved Pol II occupancy at the TSS of key cardiac genes. Our findings highlight the crucial role that transcriptional pausing plays in promoting normal gene expression levels in a cardiac developmental context.

Single cell RNA-seq analysis of spinal locomotor circuitry in larval zebrafish

Identification of the neuronal types that form the specialized circuits controlling distinct behaviors has benefited greatly from the simplicity offered by zebrafish. Electrophysiological studies have shown that additional to connectivity, understanding of circuitry requires identification of functional specializations among individual circuit components, such as those that regulate levels of transmitter release and neuronal excitability. In this study we use single cell RNA sequencing (scRNAseq) to identify the molecular bases for functional distinctions between motoneuron types that are causal to their differential roles in swimming. The primary motoneuron (PMn) in particular, expresses high levels of a unique combination of voltage-dependent ion channel types and synaptic proteins termed functional 'cassettes'. The ion channel types are specialized for promoting high frequency firing of action potentials and augmented transmitter release at the neuromuscular junction, both contributing to greater power generation. Our transcriptional profiling of spinal neurons further assigns expression of this cassette to specific interneuron types also involved in the central circuitry controlling high speed swimming and escape behaviors. Our analysis highlights the utility of scRNAseq in functional characterization of neuronal circuitry, in addition to providing a gene expression resource for studying cell type diversity.

Single-cell analysis of Rohon-Beard neurons implicates Fgf signaling in axon maintenance and cell survival

Peripheral sensory neurons are a critical part of the nervous system that transmit a multitude of sensory stimuli to the central nervous system. During larval and juvenile stages in zebrafish, this function is mediated by Rohon-Beard somatosensory neurons (RBs). RBs are optically accessible and amenable to experimental manipulation, making them a powerful system for mechanistic investigation of sensory neurons. Previous studies provided evidence that RBs fall into multiple subclasses; however, the number and molecular make up of these potential RB subtypes have not been well defined. Using a single-cell RNA sequencing (scRNA-seq) approach, we demonstrate that larval RBs in zebrafish fall into three, largely non-overlapping classes of neurons. We also show that RBs are molecularly distinct from trigeminal neurons in zebrafish. Cross-species transcriptional analysis indicates that one RB subclass is similar to a mammalian group of A-fiber sensory neurons. Another RB subclass is predicted to sense multiple modalities, including mechanical stimulation and chemical irritants. We leveraged our scRNA-seq data to determine that the fibroblast growth factor (Fgf) pathway is active in RBs. Pharmacological and genetic inhibition of this pathway led to defects in axon maintenance and RB cell death. Moreover, this can be phenocopied by treatment with dovitinib, an FDA-approved Fgf inhibitor with a common side effect of peripheral neuropathy. Importantly, dovitinib-mediated axon loss can be suppressed by loss of Sarm1, a positive regulator of neuronal cell death and axonal injury. This offers a molecular target for future clinical intervention to fight neurotoxic effects of this drug.

Molecular Analyses of V0v Spinal Interneurons and Identification of Transcriptional Regulators Downstream of Evx1 and Evx2 in these Cells

Background

V0v spinal interneurons are highly conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells.

Methods

To identify candidate members of V0v gene regulatory networks, we FAC-sorted WT and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

Results

Our data reveal two molecularly distinct subtypes of V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuronal expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression.

Conclusions

This study identifies two molecularly distinct subsets of V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of these cells. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

Haematopoietic stem and progenitor cell heterogeneity is inherited from the embryonic endothelium

Definitive haematopoietic stem and progenitor cells (HSPCs) generate erythroid, lymphoid and myeloid lineages. HSPCs are produced in the embryo via transdifferentiation of haemogenic endothelial cells in the aorta-gonad-mesonephros (AGM). HSPCs in the AGM are heterogeneous in differentiation and proliferative output, but how these intrinsic differences are acquired remains unanswered. Here we discovered that loss of microRNA (miR)-128 in zebrafish leads to an expansion of HSPCs in the AGM with different cell cycle states and a skew towards erythroid and lymphoid progenitors. Manipulating miR-128 in differentiating haemogenic endothelial cells, before their transition to HSPCs, recapitulated the lineage skewing in both zebrafish and human pluripotent stem cells. miR-128 promotes Wnt and Notch signalling in the AGM via post-transcriptional repression of the Wnt inhibitor csnk1a1 and the Notch ligand jag1b. De-repression of cskn1a1 resulted in replicative and erythroid-biased HSPCs, whereas de-repression of jag1b resulted in G2/M and lymphoid-biased HSPCs with long-term consequence on the respective blood lineages. We propose that HSPC heterogeneity arises in the AGM endothelium and is programmed in part by Wnt and Notch signalling.

The RNA-binding protein Adad1 is necessary for germ cell maintenance and meiosis in zebrafish

The double stranded RNA binding protein Adad1 (adenosine deaminase domain containing 1) is a member of the adenosine deaminase acting on RNAs (Adar) protein family with germ cell-specific expression. In mice, Adad1 is necessary for sperm differentiation, however its function outside of mammals has not been investigated. Here, through an N-ethyl-N-nitrosourea (ENU) based forward genetic screen, we identified an adad1 mutant zebrafish line that develops as sterile males. Further histological examination revealed complete lack of germ cells in adult mutant fish, however germ cells populated the gonad, proliferated, and entered meiosis in larval and juvenile fish. Although meiosis was initiated in adad1 mutant testes, the spermatocytes failed to progress beyond the zygotene stage. Thus, Adad1 is essential for meiosis and germline maintenance in zebrafish. We tested if spermatogonial stem cells were affected using nanos2 RNA FISH and a label retaining cell (LRC) assay, and found that the mutant testes had fewer LRCs and nanos2-expressing cells compared to wild-type siblings, suggesting that failure to maintain the spermatogonial stem cells resulted in germ cell loss by adulthood. To identify potential molecular processes regulated by Adad1, we sequenced bulk mRNA from mutants and wild-type testes and found mis-regulation of genes involved in RNA stability and modification, pointing to a potential broader role in post-transcriptional regulation. Our findings suggest that the RNA regulatory protein Adad1 is required for fertility through regulation of spermatogonial stem cell maintenance in zebrafish.

3' RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

Introduction: The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Methods and results: Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. Conclusions: These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

Interplay between calcium and sarcomeres directs cardiomyocyte maturation during regeneration

Zebrafish hearts can regenerate by replacing damaged tissue with new cardiomyocytes. Although the steps leading up to the proliferation of surviving cardiomyocytes have been extensively studied, little is known about the mechanisms that control proliferation and redifferentiation to a mature state. We found that the cardiac dyad, a structure that regulates calcium handling and excitation-contraction coupling, played a key role in the redifferentiation process. A component of the cardiac dyad called leucine-rich repeat-containing 10 (Lrrc10) acted as a negative regulator of proliferation, prevented cardiomegaly, and induced redifferentiation. We found that its function was conserved in mammalian cardiomyocytes. This study highlights the importance of the underlying mechanisms required for heart regeneration and their application to the generation of fully functional cardiomyocytes.

A tessellated lymphoid network provides whole-body T cell surveillance in zebrafish

Homeostatic trafficking to lymph nodes allows T cells to efficiently survey the host for cognate antigen. Nonmammalian jawed vertebrates lack lymph nodes but maintain diverse T cell pools. Here, we exploit in vivo imaging of transparent zebrafish to investigate how T cells organize and survey for antigen in an animal devoid of lymph nodes. We find that naïve-like T cells in zebrafish organize into a previously undescribed whole-body lymphoid network that supports streaming migration and coordinated trafficking through the host. This network has the cellular hallmarks of a mammalian lymph node, including naïve T cells and CCR7-ligand expressing nonhematopoietic cells, and facilitates rapid collective migration. During infection, T cells transition to a random walk that supports antigen-presenting cell interactions and subsequent activation. Our results reveal that T cells can toggle between collective migration and individual random walks to prioritize either large-scale trafficking or antigen search in situ. This lymphoid network thus facilitates whole-body T cell trafficking and antigen surveillance in the absence of a lymph node system.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 hours post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and suggest new long-term cycling populations. Focused analyses of non-skeletal muscle and the endoderm identified transcriptional profiles of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and homologs of recently discovered human best4+ enterocytes. The transcriptional regulators of these populations remain unknown, so we reconstructed gene expression trajectories to suggest candidates. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

A single-cell transcriptome atlas of the maturing zebrafish telencephalon

The zebrafish telencephalon is composed of highly specialized subregions that regulate complex behaviors such as learning, memory, and social interactions. The transcriptional signatures of the neuronal cell types in the telencephalon and the timeline of their emergence from larva to adult remain largely undescribed. Using an integrated analysis of single-cell transcriptomes of approximately 64,000 cells obtained from 6-day-postfertilization (dpf), 15-dpf, and adult telencephalon, we delineated nine main neuronal cell types in the pallium and eight in the subpallium and nominated novel marker genes. Comparing zebrafish and mouse neuronal cell types revealed both conserved and absent types and marker genes. Mapping of cell types onto a spatial larval reference atlas created a resource for anatomical and functional studies. Using this multiage approach, we discovered that although most neuronal subtypes are established early in the 6-dpf fish, some emerge or expand in number later in development. Analyzing the samples from each age separately revealed further complexity in the data, including several cell types that expand substantially in the adult forebrain and do not form clusters at the larval stages. Together, our work provides a comprehensive transcriptional analysis of the cell types in the zebrafish telencephalon and a resource for dissecting its development and function.

Tris(4-chlorophenyl)methane and tris(4-chlorophenyl)methanol disrupt pancreatic organogenesis and gene expression in zebrafish embryos

OBJECTIVES

Tris(4-chlorophenyl) methane (TCPM) and tris(4-chlorophenyl)methanol (TCPMOH) are anthropogenic environmental contaminants believed to be manufacturing byproducts of the organochlorine pesticide dichlorodiphenyltrichloroethane (DDT) due to environmental co-occurrence. TCPM and TCPMOH are persistent, bioaccumulate in the environment, and are detected in human breast milk and adipose tissues. DDT exposures have been previously shown to disrupt insulin signaling and glucoregulation, increasing risk for diabetes. We have previously shown that embryonic exposures organochlorines such as polychlorinated biphenyls disrupted pancreatic development and early embryonic glucoregulatory networks. Here, we determined the impacts of the similar compounds TCPM and TCPMOH on zebrafish pancreatic growth and gene expression following developmental exposures.

METHODS

Zebrafish embryos were exposed to 50 nM TCPM or TCPMOH beginning at 24 hr postfertilization (hpf) and exposures were refreshed daily. At 96 hpf, pancreatic growth and islet area were directly visualized in Tg(ptf1a::GFP) and Tg(insulin::GFP) embryos, respectively, using microscopy. Gene expression was assessed at 100 hpf with RNA sequencing.

RESULTS

Islet and total pancreas area were reduced by 20.8% and 13% in embryos exposed to 50 nM TCPMOH compared to controls. TCPM did not induce significant morphological changes to the developing pancreas, indicating TCPMOH, but not TCPM, impairs pancreatic development despite similarity in molecular responses. Transcriptomic responses to TCPM and TCPMOH were correlated (R²  = .903), and pathway analysis found downregulation of processes including retinol metabolism, circadian rhythm, and steroid biosynthesis.

CONCLUSION

Overall, our data suggest that TCPM and TCPMOH may be hazardous to embryonic growth and development.

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.

Transgenic IDH2R172K and IDH2R140Q zebrafish models recapitulated features of human acute myeloid leukemia

Isocitrate dehydrogenase 2 (IDH2) mutations occur in more than 15% of cytogenetically normal acute myeloid leukemia (CN-AML) but comparative studies of their roles in leukemogenesis have been scarce. We generated zebrafish models of IDH2^R172K and IDH2^R140Q AML and reported their pathologic, functional and transcriptomic features and therapeutic responses to target therapies. Transgenic embryos co-expressing FLT3^ITD and IDH2 mutations showed accentuation of myelopoiesis. As these embryos were raised to adulthood, full-blown leukemia ensued with multi-lineage dysplasia, increase in myeloblasts and marrow cellularity and splenomegaly. The leukemia cells were transplantable into primary and secondary recipients and resulted in more aggressive disease. Tg(Runx1:FLT3^ITDIDH2^R172K) but not Tg(Runx1:FLT3^ITDIDH2^R140Q) zebrafish showed an increase in T-cell development at embryonic and adult stages. Single-cell transcriptomic analysis revealed increased myeloid skewing, differentiation blockade and enrichment of leukemia-associated gene signatures in both zebrafish models. Tg(Runx1:FLT3^ITDIDH2^R172K) but not Tg(Runx1:FLT3^ITDIDH2^R140Q) zebrafish showed an increase in interferon signals at the adult stage. Leukemic phenotypes in both zebrafish could be ameliorated by quizartinib and enasidenib. In conclusion, the zebrafish models of IDH2 mutated AML recapitulated the morphologic, clinical, functional and transcriptomic characteristics of human diseases, and provided the prototype for developing zebrafish leukemia models of other genotypes that would become a platform for high throughput drug screening.

Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3' RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

A tessellated lymphoid network provides whole-body T cell surveillance in zebrafish

Homeostatic trafficking to lymph nodes allows T cells to efficiently survey the host for cognate antigen. Non-mammalian jawed vertebrates lack lymph nodes but maintain similarly diverse T cell pools. Here, we exploit in vivo imaging of transparent zebrafish to investigate how T cells organize and survey for antigen in an animal devoid of lymph nodes. We find that naïve-like T cells in zebrafish organize into a previously undescribed whole-body lymphoid network that supports streaming migration and coordinated trafficking through the host. This network has the cellular hallmarks of a mammalian lymph node, including naïve T cells and CCR7-ligand expressing non-hematopoietic cells, and facilitates rapid collective migration. During infection, T cells transition to a random walk that supports antigen presenting cell interactions and subsequent activation. Our results reveal that T cells can toggle between collective migration and individual random walks to prioritize either large-scale trafficking or antigen search in situ . This novel lymphoid network thus facilitates whole-body T cell trafficking and antigen surveillance in the absence of a lymph node system.

Significance Statement

In mammals, lymph nodes play a critical role in the initiation of adaptive immune responses by providing a dedicated place for T cells to scan antigen-presenting cells. Birds, reptiles, amphibians, and fish all maintain diverse repertoires of T cells but lack lymph nodes, raising questions about how adaptive immunity functions in lower jawed vertebrates. Here, we describe a novel network of lymphocytes in zebrafish that supports whole-body T cell trafficking and provides a site for antigen search, mirroring the function of mammalian lymph nodes. Within this network, T cells can prioritize large-scale trafficking or antigen scanning by toggling between two distinct modes of migration. This network provides valuable insights into the evolution of adaptive immunity.

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.

Cell type diversity in a developing octopus brain

Octopuses are mollusks that have evolved intricate neural systems comparable with vertebrates in terms of cell number, complexity and size. The brain cell types that control their sophisticated behavioral repertoire are still unknown. Here, we profile the cell diversity of the paralarval Octopus vulgaris brain to build a cell type atlas that comprises mostly neural cells, but also multiple glial subtypes, endothelial cells and fibroblasts. We spatially map cell types to the vertical, subesophageal and optic lobes. Investigation of cell type conservation reveals a shared gene signature between glial cells of mouse, fly and octopus. Genes related to learning and memory are enriched in vertical lobe cells, which show molecular similarities with Kenyon cells in Drosophila. We construct a cell type taxonomy revealing transcriptionally related cell types, which tend to appear in the same brain region. Together, our data sheds light on cell type diversity and evolution in the octopus brain.

Zebrafish Danio rerio myotomal muscle structure and growth from a spatial transcriptomics perspective

Fish exhibit different muscle structures and growth characteristics compared with mammals. We used a spatial transcriptomics approach and examined myotomal muscle sections from zebrafish. Adult muscles were divided into eight regions according to spatial gene expression characteristics. Slow muscle was located in the wedge-shaped region near the lateral line and at the base of the dorsal fin, intermediate muscle was located in a ribbon-shaped region adjacent to slow muscle, and fast muscle was located in the deep region of the trunk, surrounded by intermediate muscle; the interior of fast muscle was further divided into 6 parts by their transcriptomic features. Combined analysis of adult and larval data revealed that adult muscles contain specific regions similar to larval muscles. These regions showed active myogenesis and a high expression of genes associated with muscle hyperplasia. This is the first study to apply spatial transcriptomics to fish myotomal muscle structure and growth.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Single-cell profiling of transcriptome and histone modifications with EpiDamID

Recent advances in single-cell sequencing technologies have enabled simultaneous measurement of multiple cellular modalities, but the combined detection of histone post-translational modifications and transcription at single-cell resolution has remained limited. Here, we introduce EpiDamID, an experimental approach to target a diverse set of chromatin types by leveraging the binding specificities of single-chain variable fragment antibodies, engineered chromatin reader domains, and endogenous chromatin-binding proteins. Using these, we render the DamID technology compatible with the genome-wide identification of histone post-translational modifications. Importantly, this includes the possibility to jointly measure chromatin marks and transcription at the single-cell level. We use EpiDamID to profile single-cell Polycomb occupancy in mouse embryoid bodies and provide evidence for hierarchical gene regulatory networks. In addition, we map H3K9me3 in early zebrafish embryogenesis, and detect striking heterochromatic regions specific to notochord. Overall, EpiDamID is a new addition to a vast toolbox to study chromatin states during dynamic cellular processes.

GH Action in Prostate Cancer Cells Promotes Proliferation, Limits Apoptosis, and Regulates Cancer-related Gene Expression

Previous studies investigating the effects of blocking the growth hormone (GH)/insulin-like growth factor-1 (IGF-1) axis in prostate cancer found no effects of the growth hormone receptor (GHR) antagonist, pegvisomant, on the growth of grafted human prostate cancer cells in vivo. However, human GHR is not activated by mouse GH, so direct actions of GH on prostate cancer cells were not evaluated in this context. The present study addresses the species specificity of GH-GHR activity by investigating GH actions in prostate cancer cell lines derived from a mouse Pten-deletion model. In vitro cell growth was stimulated by GH and reduced by pegvisomant. These in vitro GH effects were mediated at least in part by the activation of JAK2 and STAT5. When Pten-mutant cells were grown as xenografts in mice, pegvisomant treatment dramatically reduced xenograft size, and this was accompanied by decreased proliferation and increased apoptosis. RNA sequencing of xenografts identified 1765 genes upregulated and 953 genes downregulated in response to pegvisomant, including many genes previously implicated as cancer drivers. Further evaluation of a selected subset of these genes via quantitative reverse transcription-polymerase chain reaction determined that some genes exhibited similar regulation by pegvisomant in prostate cancer cells whether treatment was in vivo or in vitro, indicating direct regulation by GH via GHR activation in prostate cancer cells, whereas other genes responded to pegvisomant only in vivo, suggesting indirect regulation by pegvisomant effects on the host endocrine environment. Similar results were observed for a prostate cancer cell line derived from the mouse transgenic adenocarcinoma of the mouse prostate (TRAMP) model.

Vinculin controls endothelial cell junction dynamics during vascular lumen formation

Blood vessel morphogenesis is driven by coordinated endothelial cell behaviors. Active remodeling of cell-cell junctions promotes cellular plasticity while preserving vascular integrity. Here, we analyze the dynamics of endothelial adherens junctions during lumen formation in angiogenic sprouts in vivo. Live imaging in zebrafish reveals that lumen expansion is accompanied by the formation of transient finger-shaped junctions. Junctional fingers are positively regulated by blood pressure, whereas flow inhibition prevents their formation. Using fluorescent reporters, we show that junctional fingers contain the mechanotransduction protein vinculin. Furthermore, genetic deletion of vinculin prevents finger formation, a junctional defect that could be rescued by transient endothelial expression of vinculin. Our findings suggest a mechanism whereby lumen expansion leads to an increase in junctional tension, triggering recruitment of vinculin and formation of junctional fingers. We propose that endothelial cells employ force-dependent junctional remodeling to counteract external forces in order to maintain vascular integrity during sprouting angiogenesis.

Single-cell Iso-Sequencing enables rapid genome annotation for scRNAseq analysis

Single-cell RNA sequencing is a powerful technique that continues to expand across various biological applications. However, incomplete 3'-UTR annotations can impede single-cell analysis resulting in genes that are partially or completely uncounted. Performing single-cell RNA sequencing with incomplete 3'-UTR annotations can hinder the identification of cell identities and gene expression patterns and lead to erroneous biological inferences. We demonstrate that performing single-cell isoform sequencing in tandem with single-cell RNA sequencing can rapidly improve 3'-UTR annotations. Using threespine stickleback fish (Gasterosteus aculeatus), we show that gene models resulting from a minimal embryonic single-cell isoform sequencing dataset retained 26.1% greater single-cell RNA sequencing reads than gene models from Ensembl alone. Furthermore, pooling our single-cell sequencing isoforms with a previously published adult bulk Iso-Seq dataset from stickleback, and merging the annotation with the Ensembl gene models, resulted in a marginal improvement (+0.8%) over the single-cell isoform sequencing only dataset. In addition, isoforms identified by single-cell isoform sequencing included thousands of new splicing variants. The improved gene models obtained using single-cell isoform sequencing led to successful identification of cell types and increased the reads identified of many genes in our single-cell RNA sequencing stickleback dataset. Our work illuminates single-cell isoform sequencing as a cost-effective and efficient mechanism to rapidly annotate genomes for single-cell RNA sequencing.

Regeneration of Jaw Joint Cartilage in Adult Zebrafish

The poor intrinsic repair capacity of mammalian joint cartilage likely contributes to the high incidence of arthritis worldwide. Adult zebrafish can regenerate many structures that show limited or no healing capacity in mammals, including the jawbone. To test whether zebrafish can also regenerate damaged joints, we developed a surgical injury model in which the zebrafish jaw joint is destabilized via transection of the major jaw joint ligament, the interopercular-mandibular (IOM). Unilateral transection of the IOM ligament in 1-year-old fish resulted in an initial reduction of jaw joint cartilage by 14 days, with full regeneration of joint cartilage by 28 days. Joint cartilage regeneration involves the re-entry of articular chondrocytes into the cell cycle and the upregulated expression of sox10, a marker of developing chondrocytes in the embryo that becomes restricted to a subset of joint chondrocytes in adults. Genetic ablation of these sox10-expressing chondrocytes shows that they are essential for joint cartilage regeneration. To uncover the potential source of new chondrocytes during joint regeneration, we performed single-cell RNA sequencing of the uninjured adult jaw joint and identified multiple skeletal, connective tissue, and fibroblast subtypes. In particular, we uncovered a joint-specific periosteal population expressing coch and grem1a, with the jaw joint chondrocytes marked by grem1a expression during regeneration. Our findings demonstrate the capacity of zebrafish to regenerate adult joint cartilage and identify candidate cell types that can be tested for their roles in regenerative response.

Choice of pre-processing pipeline influences clustering quality of scRNA-seq datasets

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) has quickly become one of the most dominant techniques in modern transcriptome assessment. In particular, 10X Genomics' Chromium system, with its high throughput approach, turn key and thorough user guide made this cutting-edge technique accessible to many laboratories using diverse animal models. However, standard pre-processing, including the alignment and cell filtering pipelines might not be ideal for every organism or tissue. Here we applied an alternative strategy, based on the pseudoaligner kallisto, on twenty-two publicly available single cell sequencing datasets from a wide range of tissues of eight organisms and compared the results with the standard 10X Genomics' Cell Ranger pipeline.

RESULTS

In most of the tested samples, kallisto produced higher sequencing read alignment rates and total gene detection rates in comparison to Cell Ranger. Although datasets processed with Cell Ranger had higher cell counts, outside of human and mouse datasets, these additional cells were routinely of low quality, containing low gene detection rates. Thorough downstream analysis of one kallisto processed dataset, obtained from the zebrafish pineal gland, revealed clearer clustering, allowing the identification of an additional photoreceptor cell type that previously went undetected. The finding of the new cluster suggests that the photoreceptive pineal gland is essentially a bi-chromatic tissue containing both green and red cone-like photoreceptors and implies that the alignment and pre-processing pipeline can affect the discovery of biologically-relevant cell types.

CONCLUSION

While Cell Ranger favors higher cell numbers, using kallisto results in datasets with higher median gene detection per cell. We could demonstrate that cell type identification was not hampered by the lower cell count, but in fact improved as a result of the high gene detection rate and the more stringent filtering. Depending on the acquired dataset, it can be beneficial to favor high quality cells and accept a lower cell count, leading to an improved classification of cell types.