Characterization of the Zebrafish Cell Landscape at Single-Cell Resolution

CUT&Tag applied to zebrafish adult tail fins reveals a return of embryonic H3K4me3 patterns during regeneration

Regenerative potential is governed by a complex process of transcriptional reprogramming, involving chromatin reorganization and dynamics in transcription factor binding patterns throughout the genome. The degree to which chromatin and epigenetic changes contribute to this process remains only partially understood. Here we provide a modified CUT&Tag protocol suitable for improved characterization and interrogation of changes in chromatin modifications during adult fin regeneration in zebrafish. Our protocol generates data that recapitulates results from previously published ChIP-Seq methods, requires far fewer cells as input, and significantly improves signal to noise ratios. We deliver high-resolution enrichment maps for H3K4me3 of uninjured and regenerating fin tissues. During regeneration, we find that H3K4me3 levels increase over gene promoters which become transcriptionally active and genes which lose H3K4me3 become silenced. Interestingly, these reprogramming events recapitulate the H3K4me3 patterns observed in developing fin folds of 24-h old zebrafish embryos. Our results indicate that changes in genomic H3K4me3 patterns during fin regeneration occur in a manner consistent with reactivation of developmental programs, demonstrating CUT&Tag to be an effective tool for profiling chromatin landscapes in regenerating tissues.

scMultiome analysis identifies a single caudal hindbrain compartment in the developing zebrafish nervous system

BACKGROUND

A key step in nervous system development involves the coordinated control of neural progenitor specification and positioning. A long-standing model for the vertebrate CNS postulates that transient anatomical compartments - known as neuromeres - function to position neural progenitors along the embryonic anteroposterior neuraxis. Such neuromeres are apparent in the embryonic hindbrain - that contains six rhombomeres with morphologically apparent boundaries - but other neuromeres lack clear morphological boundaries and have instead been defined by different criteria, such as differences in gene expression patterns and the outcomes of transplantation experiments. Accordingly, the caudal hindbrain (CHB) posterior to rhombomere (r) 6 has been variably proposed to contain from two to five 'pseudo-rhombomeres', but the lack of comprehensive molecular data has precluded a detailed definition of such structures.

METHODS

We used single-cell Multiome analysis, which allows simultaneous characterization of gene expression and chromatin state of individual cell nuclei, to identify and characterize CHB progenitors in the developing zebrafish CNS.

RESULTS

We identified CHB progenitors as a transcriptionally distinct population, that also possesses a unique profile of accessible transcription factor binding motifs, relative to both r6 and the spinal cord. This CHB population can be subdivided along its dorsoventral axis based on molecular characteristics, but we do not find any molecular evidence that it contains multiple pseudo-rhombomeres. We further observe that the CHB is closely related to r6 at the earliest embryonic stages, but becomes more divergent over time, and that it is defined by a unique gene regulatory network.

CONCLUSIONS

We conclude that the early CHB represents a single neuromere compartment that cannot be molecularly subdivided into pseudo-rhombomeres and that it may share an embryonic origin with r6.

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.

Vitamin E alleviates chlorpyrifos induced glutathione depletion, lipid peroxidation and iron accumulation to inhibit ferroptosis in hepatocytes and mitigate toxicity in zebrafish

Organophosphates, a widely used group of pesticides, can cause severe toxicity in human beings and other non-target organisms. Liver, being the primary site for xenobiotic metabolism, is extremely vulnerable to xenobiotic-induced toxicity. Considering the numerous vital functions performed by the liver, including xenobiotic detoxification, protecting this organ from the ubiquitous pesticides in our food and environment is essential for maintaining homeostasis. In this study, we have investigated the impact of the organophosphate pesticide, Chlorpyrifos (CPF), on zebrafish liver at a concentration (300 μg/L) which is environmentally realistic. We have also demonstrated the role of dietary supplementation of α-tocopherol or Vitamin E (Vit E) (500 mg/kg feed) in mitigating pesticide-induced liver toxicity. Mechanistically, we showed that Vit E resulted in significant elevation of the Nrf2 and its downstream antioxidant enzyme activities and gene expressions, especially that of GST and GPx, resulting in reduction of CPF-induced intracellular lipid ROS and hepatic LPO. Further interrogation, such as analysis of GSH: GSSG ratio, intracellular iron concentration, iron metabolizing genes, mitochondrial dysfunction etc. revealed that CPF induces ferroptosis which can be reversed by Vit E supplementation. Ultimately, reduced concentration of CPF in zebrafish serum and flesh highlighted the role of Vit E in ameliorating CPF toxicity.

scRNA-seq reveals the diversity of the developing cardiac cell lineage and molecular players in heart rhythm regulation

We utilized scRNA-seq to delineate the diversity of cell types in the zebrafish heart. Transcriptome profiling of over 50,000 cells at 48 and 72 hpf defined at least 18 discrete cell lineages of the developing heart. Utilizing well-established gene signatures, we identified a population of cells likely to be the primary pacemaker and characterized the transcriptome profile defining this critical cell type. Two previously uncharacterized genes, atp1b3b and colec10, were found to be enriched in the sinoatrial cardiomyocytes. CRISPR/Cas9-mediated knockout of these two genes significantly reduced heart rate, implicating their role in cardiac development and conduction. Additionally, we describe other cardiac cell lineages, including the endothelial and neural cells, providing their expression profiles as a resource. Our results established a detailed atlas of the developing heart, providing valuable insights into cellular and molecular mechanisms, and pinpointed potential new players in heart rhythm regulation.

Status of single-cell RNA sequencing for reproductive toxicology in zebrafish and the transcriptomic trade-off

The utilization of transcriptomic studies identifying profiles of gene expression, especially in toxicogenomics, has catapulted next-generation sequencing to the forefront of reproductive toxicology. An innovative yet underutilized RNA sequencing technique emerging into this field is single-cell RNA sequencing (scRNA-seq), which provides sequencing at the individual cellular level of gonad tissue. ScRNA-seq provides a novel and unique perspective for identifying distinct cellular profiles, including identification of rare cell subtypes. The specificity of scRNA-seq is a powerful tool for reproductive toxicity research, especially for translational animal models including zebrafish. Studies to date not only have focused on 'tissue atlassing' or characterizing what cell types make up different tissues but have also begun to include toxicant exposure as a factor that this review aims to explore. Future scRNA-seq studies will contribute to understanding exposure-induced outcomes; however, the trade-offs with traditional methods need to be considered.

Single cell RNA analysis uncovers the cell differentiation and functionalization for air breathing of frog lung

The evolution and development of vertebrate lungs have been widely studied due to their significance in terrestrial adaptation. Amphibians possess the most primitive lungs among tetrapods, underscoring their evolutionary importance in bridging the transition from aquatic to terrestrial life. However, the intricate process of cell differentiation during amphibian lung development remains poorly understood. Using single-cell RNA sequencing, we identify 13 cell types in the developing lungs of a land-dwelling frog (Microhyla fissipes). We elucidate the differentiation trajectories and mechanisms of mesenchymal cells, identifying five cell fates and their respective driver genes. Using temporal dynamics analyses, we reveal the gene expression switches of epithelial cells, which facilitate air breathing during metamorphosis. Furthermore, by integrating the published data from another amphibian and two terrestrial mammals, we illuminate both conserved and divergent cellular repertoires during the evolution of tetrapod lungs. These findings uncover the frog lung cell differentiation trajectories and functionalization for breathing in air and provide valuable insights into the cell-type evolution of vertebrate lungs.

Single cell sequencing of zebrafish kidney marrows reveals AHR2-dependent endogenous regulation of hematopoiesis

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor that mediates a wide range of biological and toxicological responses. While largely studied in ligand-activated toxicant responses, AHR also plays important roles in endogenous physiological processes. We leveraged single cell sequencing and an AHR2 knockout zebrafish line to investigate the role of AHR2 in regulating hematopoiesis (production and differentiation of red and white blood cells from hematopoietic stem cells). Our objectives were to determine if absence of AHR2- 1) alters proportions of immune cell populations and/or 2) impacts gene expression within individual immune cell types. We dissected kidney marrow (organ of hematopoiesis in zebrafish) from adult wildtype and AHR2 knockout zebrafish (N=3/genotype), isolated single cells and sequenced ~ 5000 cells/sample (10X Genomics). We identified 14 cell clusters representing the expected major blood (erythrocytes, thrombocytes), immune (B cells, macrophages, lymphoid cells, granulocytes, etc), progenitors and kidney cell populations. We focused our analyses only on the progenitor and mature immune cell types. While there were no genotype-specific differences in proportion of individual cell types, gene expression differences were observed within several cell types. For known genes, such as rrm2, changes were up to 2000-fold, signifying their importance in AHR2-hematopoesis interaction. Several of the known genes are also identified as markers of carcinoma cells for an array of cancer types. However, many of the dysregulated genes are poorly annotated, limiting our ability to examine biological processes and pathways dysregulated on AHR2 mutation. Nevertheless, our study indicates that AHR2 plays an important endogenous role in hematopoiesis. Future work will focus on better characterizing anatomy of dysregulated genes and their functions in hematopoiesis.

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.

CUT&Tag Applied to Zebrafish Adult Tail Fins Reveals a Return of Embryonic H3K4me3 Patterns During Regeneration

Regenerative potential is governed by a complex process of transcriptional reprogramming, involving chromatin reorganization and dynamics in transcription factor binding patterns throughout the genome. The degree to which chromatin and epigenetic changes contribute to this process remains partially understood. Here we provide a modified CUT&Tag protocol suitable for improved characterization and interrogation of epigenetic changes during adult fin regeneration in zebrafish. Our protocol generates data that recapitulates results from previously published ChIP-Seq methods, requires far fewer cells as input, and significantly improves signal to noise ratios. We deliver high-resolution enrichment maps for H3K4me3 of uninjured and regenerating fin tissues. During regeneration, we find that H3K4me3 levels increase over gene promoters which become transcriptionally active and genes which lose H3K4me3 become silenced. Interestingly, these epigenetic reprogramming events recapitulate the H3K4me3 patterns observed in developing fin folds of 24-hour old zebrafish embryos. Our results indicate that changes in genomic H3K4me3 patterns during fin regeneration occur in a manner consistent with reactivation of developmental programs, demonstrating CUT&Tag to be an effective tool for profiling chromatin landscapes in regenerating tissues.

Slc43a2+ T cell metastasis from spleen to brain in RGNNV infected teleost

The origin of T cells in the teleost's brain is unclear. While viewing the central nervous system (CNS) as immune privileged has been widely accepted, previous studies suggest that T cells residing in the thymus but not in the spleen of the teleost play an essential role in communicating with the peripheral organs. Here, we identified nine T cell subpopulations in the thymus and spleen of orange-spotted grouper (Epinephelus coioices) through single-cell RNA-sequencing analysis. After viral CNS infection with red-spotted grouper nervous necrosis virus (RGNNV), the number of slc43a2+ T cells synchronously increased in the spleen and brain. During the infection tests in asplenic zebrafish (tlx1▲ zebrafish model), no increase in the number of slc43a2+ T cells was observed in the brain. Single-cell transcriptomic analysis indicated that slc43a2+ T cells mature and functionally differentiate within the spleen and then migrate into the brain to trigger an immune response. This study suggests a novel route for T cell migration from the spleen to the brain during viral infection in fish.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

Paternal starvation affects metabolic gene expression during zebrafish offspring development and lifelong fitness

Dietary restriction in the form of fasting is a putative key to a healthier and longer life, but these benefits may come at a trade-off with reproductive fitness and may affect the following generation(s). The potential inter- and transgenerational effects of long-term fasting and starvation are particularly poorly understood in vertebrates when they originate from the paternal line. We utilised the externally fertilising zebrafish amenable to a split-egg clutch design to explore the male-specific effects of fasting/starvation on fertility and fitness of offspring independently of maternal contribution. Eighteen days of fasting resulted in reduced fertility in exposed males. While average offspring survival was not affected, we detected increased larval growth rate in F1 offspring from starved males and more malformed embryos at 24 h post-fertilisation in F2 offspring produced by F1 offspring from starved males. Comparing the transcriptomes of F1 embryos sired by starved and fed fathers revealed robust and reproducible increased expression of muscle composition genes but lower expression of lipid metabolism and lysosome genes in embryos from starved fathers. A large proportion of these genes showed enrichment in the yolk syncytial layer suggesting gene regulatory responses associated with metabolism of nutrients through paternal effects on extra-embryonic tissues which are loaded with maternal factors. We compared the embryo transcriptomes to published adult transcriptome datasets and found comparable repressive effects of starvation on metabolism-associated genes. These similarities suggest a physiologically relevant, directed and potentially adaptive response transmitted by the father, independently from the offspring's nutritional state, which was defined by the mother.

miR-210 promotes immune- and suppresses oocyte meiosis-related genes in the zebrafish ovarian cells

microRNA-210 (miRNA), a well-documented miRNA, has been implicated in a myriad of biological processes, including responses to hypoxia, angiogenesis, cell proliferation, and male infertility in humans. However, a comprehensive understanding of its functions in fish requires further investigation. This study pursued to elucidate the downstream effect of dre-miR-210-5p on primary ovarian cell culture in zebrafish (Danio rerio), an animal model. A protocol was settled down by incubations with either an miR-210 mimic or a scrambled miRNA in the isolated ovaries. RNA-sequencing analysis identified ∼6000 differentially expressed target genes revealing that downregulated genes were associated with reproduction-related pathways while immune-related pathways displayed an upregulated pattern. To identify molecular markers, predicted target genes were classified into reproduction and immune cell types. These findings underscore the existence of a profound interplay between the reproductive and immune systems, with miR-210 emerging as a pivotal player in orchestrating transcriptomic alterations within fish ovaries.

Compartmentalization and synergy of osteoblasts drive bone formation in the regenerating fin

Zebrafish regenerate their fins which involves a component of cell plasticity. It is currently unclear how regenerate cells divide labor to allow for appropriate growth and patterning. Here, we studied lineage relationships of fluorescence-activated cell sorting-enriched epidermal, bone-forming (osteoblast), and (non-osteoblast) blastemal fin regenerate cells by single-cell RNA sequencing, lineage tracing, targeted osteoblast ablation, and electron microscopy. Most osteoblasts in the outgrowing regenerate derive from osterix+ osteoblasts, while mmp9+ cells reside at segment joints. Distal blastema cells contribute to distal osteoblast progenitors, suggesting compartmentalization of the regenerating appendage. Ablation of osterix+ osteoblasts impairs segment joint and bone matrix formation and decreases regenerate length which is partially compensated for by distal regenerate cells. Our study characterizes expression patterns and lineage relationships of rare fin regenerate cell populations, indicates inherent detection and compensation of impaired regeneration, suggests variable dependence on growth factor signaling, and demonstrates zonation of the elongating fin regenerate.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Convergent differentiation of multiciliated cells

Multiciliated cells (MCCs) are epithelial cells that control body fluid flow and contribute to the clearance of pathogenic microbes and other particles from the airways, egg transport in oviducts, and circulation of cerebrospinal fluid in the central nervous system. Although MCCs have shared functions to control fluid flow via coordinated motility of multiple ciliary structures, they are found in multiple mammalian tissues originating from distinct germ layers and differentiate via distinct developmental pathways. To understand the similarities and differences of MCCs in multiple tissues, we investigated single-cell transcriptome data of nasal epithelial cells, bronchial tubes, fallopian tubes, and ependymal cells in the subventricular zone from humans and mice by cross-species data integration. Expression of cilia-associated genes was indistinguishable between these MCCs, although cell populations had unique properties by the species and tissue, demonstrating that they share the same final differentiation status for ciliary functions. We further analyzed the final differentiation step of MCCs from their distinctive progenitors and confirmed their convergent gene set expression for ciliogenesis at the final step. These results may provide new insight into understanding ciliogenesis during the developmental process.

Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei

The innate immune system is the first line of defense in multicellular organisms. Danio rerio is widely considered a promising model for IIS-related research, with the most amount of scRNAseq data available among Teleostei. We summarized the scRNAseq and spatial transcriptomics experiments related to the IIS for zebrafish and other Teleostei from the GEO NCBI and the Single-Cell Expression Atlas. We found a considerable number of scRNAseq experiments at different stages of zebrafish development in organs such as the kidney, liver, stomach, heart, and brain. These datasets could be further used to conduct large-scale meta-analyses and to compare the IIS of zebrafish with the mammalian one. However, only a small number of scRNAseq datasets are available for other fish (turbot, salmon, cavefish, and dark sleeper). Since fish biology is very diverse, it would be a major mistake to use zebrafish alone in fish immunology studies. In particular, there is a special need for new scRNAseq experiments involving nonmodel Teleostei, e.g., long-lived species, cancer-resistant fish, and various fish ecotypes.

scMultiome analysis identifies embryonic hindbrain progenitors with mixed rhombomere identities

Rhombomeres serve to position neural progenitors in the embryonic hindbrain, thereby ensuring appropriate neural circuit formation, but the molecular identities of individual rhombomeres and the mechanism whereby they form has not been fully established. Here, we apply scMultiome analysis in zebrafish to molecularly resolve all rhombomeres for the first time. We find that rhombomeres become molecularly distinct between 10hpf (end of gastrulation) and 13hpf (early segmentation). While the embryonic hindbrain transiently contains alternating odd- versus even-type rhombomeres, our scMultiome analyses do not detect extensive odd versus even molecular characteristics in the early hindbrain. Instead, we find that each rhombomere displays a unique gene expression and chromatin profile. Prior to the appearance of distinct rhombomeres, we detect three hindbrain progenitor clusters (PHPDs) that correlate with the earliest visually observed segments in the hindbrain primordium that represent prospective rhombomere r2/r3 (possibly including r1), r4, and r5/r6, respectively. We further find that the PHPDs form in response to Fgf and RA morphogens and that individual PHPD cells co-express markers of multiple mature rhombomeres. We propose that the PHPDs contain mixed-identity progenitors and that their subdivision into individual rhombomeres requires the resolution of mixed transcription and chromatin states.

Widespread sex dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism for aging research. Here, we describe a multitissue, single-cell gene expression atlas of female and male blood, kidney, liver, and spleen. We annotate 22 cell types, define marker genes, and infer differentiation trajectories. We find pervasive sex-dimorphic gene expression across cell types. Sex-dimorphic genes tend to be linked to lipid metabolism, consistent with clear differences in lipid storage in female vs. male turquoise killifish livers. We use machine learning to predict sex using single-cell gene expression and identify potential markers for molecular sex identity. As a proof of principle, we show that our atlas can be used to deconvolute existing bulk RNA sequencing (RNA-seq) data to obtain accurate estimates of cell type proportions. This atlas can be a resource to the community that could be leveraged to develop cell-type-specific expression in transgenic animals.

Benchmarking strategies for cross-species integration of single-cell RNA sequencing data

The growing number of available single-cell gene expression datasets from different species creates opportunities to explore evolutionary relationships between cell types across species. Cross-species integration of single-cell RNA-sequencing data has been particularly informative in this context. However, in order to do so robustly it is essential to have rigorous benchmarking and appropriate guidelines to ensure that integration results truly reflect biology. Here, we benchmark 28 combinations of gene homology mapping methods and data integration algorithms in a variety of biological settings. We examine the capability of each strategy to perform species-mixing of known homologous cell types and to preserve biological heterogeneity using 9 established metrics. We also develop a new biology conservation metric to address the maintenance of cell type distinguishability. Overall, scANVI, scVI and SeuratV4 methods achieve a balance between species-mixing and biology conservation. For evolutionarily distant species, including in-paralogs is beneficial. SAMap outperforms when integrating whole-body atlases between species with challenging gene homology annotation. We provide our freely available cross-species integration and assessment pipeline to help analyse new data and develop new algorithms.

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.

Zebrafish (Danio rerio) as a translational model for neuro-immune interactions in the enteric nervous system in autism spectrum disorders

Autism spectrum disorders (ASD) affect about 1% of the population and are strongly associated with gastrointestinal diseases creating shortcomings in quality of life. Multiple factors contribute to the development of ASD and although neurodevelopmental deficits are central, the pathogenesis of the condition is complex and the high prevalence of intestinal disorders is poorly understood. In agreement with the prominent research establishing clear bidirectional interactions between the gut and the brain, several studies have made it evident that such a relation also exists in ASD. Thus, dysregulation of the gut microbiota and gut barrier integrity may play an important role in ASD. However, only limited research has investigated how the enteric nervous system (ENS) and intestinal mucosal immune factors may impact on the development of ASD-related intestinal disorders. This review focuses on the mechanistic studies that elucidate the regulation and interactions between enteric immune cells, residing gut microbiota and the ENS in models of ASD. Especially the multifaceted properties and applicability of zebrafish (Danio rerio) for the study of ASD pathogenesis are assessed in comparison to studies conducted in rodent models and humans. Advances in molecular techniques and in vivo imaging, combined with genetic manipulation and generation of germ-free animals in a controlled environment, appear to make zebrafish an underestimated model of choice for the study of ASD. Finally, we establish the research gaps that remain to be explored to further our understanding of the complexity of ASD pathogenesis and associated mechanisms that may lead to intestinal disorders.

Sex-biased gene expression at single-cell resolution: cause and consequence of sexual dimorphism

Gene expression differences between males and females are thought to be key for the evolution of sexual dimorphism, and sex-biased genes are often used to study the molecular footprint of sex-specific selection. However, gene expression is often measured from complex aggregations of diverse cell types, making it difficult to distinguish between sex differences in expression that are due to regulatory rewiring within similar cell types and those that are simply a consequence of developmental differences in cell-type abundance. To determine the role of regulatory versus developmental differences underlying sex-biased gene expression, we use single-cell transcriptomic data from multiple somatic and reproductive tissues of male and female guppies, a species that exhibits extensive phenotypic sexual dimorphism. Our analysis of gene expression at single-cell resolution demonstrates that nonisometric scaling between the cell populations within each tissue and heterogeneity in cell-type abundance between the sexes can influence inferred patterns of sex-biased gene expression by increasing both the false-positive and false-negative rates. Moreover, we show that, at the bulk level, the subset of sex-biased genes that are the product of sex differences in cell-type abundance can significantly confound patterns of coding-sequence evolution. Taken together, our results offer a unique insight into the effects of allometry and cellular heterogeneity on perceived patterns of sex-biased gene expression and highlight the power of single-cell RNA-sequencing in distinguishing between sex-biased genes that are the result of regulatory change and those that stem from sex differences in cell-type abundance, and hence are a consequence rather than a cause of sexual dimorphism.

Widespread sex-dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism to study vertebrate aging. Here we describe the first multi-tissue, single-cell gene expression atlas of female and male turquoise killifish tissues comprising immune and metabolic cells from the blood, kidney, liver, and spleen. We were able to annotate 22 distinct cell types, define associated marker genes, and infer differentiation trajectories. Using this dataset, we found pervasive sex-dimorphic gene expression across cell types, especially in the liver. Sex-dimorphic genes tended to be involved in processes related to lipid metabolism, and indeed, we observed clear differences in lipid storage in female vs. male turquoise killifish livers. Importantly, we use machine-learning to predict sex using single-cell gene expression in our atlas and identify potential transcriptional markers for molecular sex identity in this species. As proof-of-principle, we show that our atlas can be used to deconvolute existing liver bulk RNA-seq data in this species to obtain accurate estimates of cell type proportions across biological conditions. We believe that this single-cell atlas can be a resource to the community that could notably be leveraged to identify cell type-specific genes for cell type-specific expression in transgenic animals.

Zebrafish Galectin 3 binding protein is the target antigen of the microglial 4C4 monoclonal antibody

BACKGROUND

Two decades ago, the fish-specific monoclonal antibody 4C4 was found to be highly reactive to zebrafish microglia, the macrophages of the central nervous system. This has resulted in 4C4 being widely used, in combination with available fluorescent transgenic reporters to identify and isolate microglia. However, the target protein of 4C4 remains unidentified, which represents a major caveat. In addition, whether the 4C4 expression pattern is strictly restricted to microglial cells in zebrafish has never been investigated.

RESULTS

Having demonstrated that 4C4 is able to capture its native antigen from adult brain lysates, we used immunoprecipitation/mass-spectrometry, coupled to recombinant expression analyses, to identify its target. The cognate antigen was found to be a paralog of Galectin 3 binding protein (Lgals3bpb), known as MAC2-binding protein in mammals. Notably, 4C4 did not recognize other paralogs, demonstrating specificity. Moreover, our data show that Lgals3bpb expression, while ubiquitous in microglia, also identifies leukocytes in the periphery, including populations of gut and liver macrophages.

CONCLUSIONS

The 4C4 monoclonal antibody recognizes Lgals3bpb, a predicted highly glycosylated protein whose function in the microglial lineage is currently unknown. Identification of Lgals3bpb as a new pan-microglia marker will be fundamental in forthcoming studies using the zebrafish model.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

scMultiome analysis identifies embryonic hindbrain progenitors with mixed rhombomere identities

Rhombomeres serve to position neural progenitors in the embryonic hindbrain, thereby ensuring appropriate neural circuit formation, but the molecular identities of individual rhombomeres and the mechanism whereby they form have not been fully established. Here we apply scMultiome analysis in zebrafish to molecularly resolve all rhombomeres for the first time. We find that rhombomeres become molecularly distinct between 10hpf (end of gastrulation) and 13hpf (early segmentation). While the mature hindbrain consists of alternating odd- versus even-type rhombomeres, our scMultiome analyses do not detect extensive odd versus even characteristics in the early hindbrain. Instead, we find that each rhombomere displays a unique gene expression and chromatin profile. Prior to the appearance of distinct rhombomeres, we detect three hindbrain progenitor clusters (PHPDs) that correlate with the earliest visually observed segments in the hindbrain primordium and that represent prospective rhombomere r2/r3 (possibly including r1), r4 and r5/r6, respectively. We further find that the PHPDs form in response to Fgf and RA morphogens and that individual PHPD cells co-express markers of multiple mature rhombomeres. We propose that the PHPDs contain mixed-identity progenitors and that their subdivision into individual mature rhombomeres requires resolution of mixed transcription and chromatin states.

Construction of a cross-species cell landscape at single-cell level

Individual cells are basic units of life. Despite extensive efforts to characterize the cellular heterogeneity of different organisms, cross-species comparisons of landscape dynamics have not been achieved. Here, we applied single-cell RNA sequencing (scRNA-seq) to map organism-level cell landscapes at multiple life stages for mice, zebrafish and Drosophila. By integrating the comprehensive dataset of > 2.6 million single cells, we constructed a cross-species cell landscape and identified signatures and common pathways that changed throughout the life span. We identified structural inflammation and mitochondrial dysfunction as the most common hallmarks of organism aging, and found that pharmacological activation of mitochondrial metabolism alleviated aging phenotypes in mice. The cross-species cell landscape with other published datasets were stored in an integrated online portal-Cell Landscape. Our work provides a valuable resource for studying lineage development, maturation and aging.

Alterations in immune cell heterogeneities in the brain of aged zebrafish using single-cell resolution

Immunocytes, including the microglia, are crucial in the neurodegenerative process in old people. However, the understanding of regarding microglia heterogeneity and other involved immunocytes remains elusive. We analyzed 26,456 immunocytes from 12-and 26-month-old zebrafish brains at single-cell resolution. Microglia and T lymphocytes were detected in the brain at both time points. Two types of microglia were annotated, namely, ac⁺ microglia and xr⁺ microglia, which were clustered into subsets 1, 2, 3, 4, 5, and subsets 6, 7, 8, 9, respectively. Diversified microglia predominated the adult brains and cooperated with T cells to perform the functions of immune response and neuronal nutrition. We validated the specific microglia markers. The novel transgenic lines, Tg(lgals3bpb:eGFP) and Tg(apoc1:eGFP), were created, which faithfully labeled ac⁺ microglia and served as valuable labeling tools. However, the microglia population reduced while T cells of six subtypes intriguingly increased to serve as the primary immune cells in aged brains. Unlike in 12-month-old brains, T cells, together with microglia, exhibited a coordinated signature of inflammation in the 26-month-old brains. Our findings revealed the immunocytes atlas in aged zebrafish brains. It implied the involvement of microglia and T cells in the progression of neurodegeneration in aging.

Bulk RNA-seq and scRNA-seq analysis reveal an activation of immune response and compromise of secretory function in major salivary glands of obese mice

Obesity affects the function of multiple organs/tissues including the exocrine organ salivary glands. However, the effects of obesity on transcriptomes and cell compositions in the salivary glands have yet been studied by bulk RNA-sequencing and single-cell RNA-sequencing. Besides, the cell types in the sublingual gland, one of the three major salivary glands, have yet been characterized by the approach of single-cell RNA-sequencing. In this report, we find that the histological structure of the three major salivary glands are not obviously affected in the obese mice. Bulk RNA-sequencing analysis shows that the most prominent changes observed in the three major salivary glands of the obese mice are the mobilization of transcriptomes related to the immune response and down-regulation of genes related to the secretory function of the salivary glands. Based on single-cell RNA-sequencing analysis, we identify and annotate 17 cell clusters in the sublingual gland for the first time, and find that obesity alters the relative compositions of immune cells and secretory cells in the major glands of obese mice. Integrative analysis of the bulk RNA-sequencing and single-cell RNA-sequencing data confirms the activation of immune response genes and compromise of secretory function in the three major salivary glands of obese mice. Consequently, the secretion of extracellular matrix proteins is significantly reduced in the three major salivary glands of obese mice. These results provide new molecular insights into understanding the effect of obesity on salivary glands.

Deep learning of cross-species single-cell landscapes identifies conserved regulatory programs underlying cell types

Despite extensive efforts to generate and analyze reference genomes, genetic models to predict gene regulation and cell fate decisions are lacking for most species. Here, we generated whole-body single-cell transcriptomic landscapes of zebrafish, Drosophila and earthworm. We then integrated cell landscapes from eight representative metazoan species to study gene regulation across evolution. Using these uniformly constructed cross-species landscapes, we developed a deep-learning-based strategy, Nvwa, to predict gene expression and identify regulatory sequences at the single-cell level. We systematically compared cell-type-specific transcription factors to reveal conserved genetic regulation in vertebrates and invertebrates. Our work provides a valuable resource and offers a new strategy for studying regulatory grammar in diverse biological systems.

Unraveling Differential Transcriptomes and Cell Types in Zebrafish Larvae Intestine and Liver

The zebrafish intestine and liver, as in other vertebrates, are derived from the endoderm. Great effort has been devoted to deciphering the molecular mechanisms controlling the specification and development of the zebrafish intestine and liver; however, genome-wide comparison of the transcriptomes between these two organs at the larval stage remains unexplored. There is a lack of extensive identification of feature genes marking specific cell types in the zebrafish intestine and liver at 5 days post-fertilization, when the larval fish starts food intake. In this report, through RNA sequencing and single-cell RNA sequencing of intestines and livers separately dissected from wild-type zebrafish larvae at 5 days post-fertilization, together with the experimental validation of 47 genes through RNA whole-mount in situ hybridization, we identified not only distinctive transcriptomes for the larval intestine and liver, but also a considerable number of feature genes for marking the intestinal bulb, mid-intestine and hindgut, and for marking hepatocytes and cholangiocytes. Meanwhile, we identified 135 intestine- and 97 liver-enriched transcription factor genes in zebrafish larvae at 5 days post-fertilization. Our findings provide rich molecular and cellular resources for studying cell patterning and specification during the early development of the zebrafish intestine and liver.

Single-cell transcriptomic analysis of vascular endothelial cells in zebrafish embryos

Vascular endothelial cells exhibit substantial phenotypic and transcriptional heterogeneity which is established during early embryogenesis. However, the molecular mechanisms involved in establishing endothelial cell diversity are still not well understood. Zebrafish has emerged as an advantageous model to study vascular development. Despite its importance, the single-cell transcriptomic profile of vascular endothelial cells during zebrafish development is still missing. To address this, we applied single-cell RNA-sequencing (scRNA-seq) of vascular endothelial cells isolated from zebrafish embryos at the 24 hpf stage. Six distinct clusters or subclusters related to vascular endothelial cells were identified which include arterial, two venous, cranial, endocardial and endothelial progenitor cell subtypes. Furthermore, we validated our findings by characterizing novel markers for arterial, venous, and endocardial cells. We experimentally confirmed the presence of two transcriptionally different venous cell subtypes, demonstrating heterogeneity among venous endothelial cells at this early developmental stage. This dataset will be a valuable resource for future functional characterization of vascular endothelial cells and interrogation of molecular mechanisms involved in the establishment of their heterogeneity and cell-fate decisions.

Cell landscape of larval and adult Xenopus laevis at single-cell resolution

The rapid development of high-throughput single-cell RNA sequencing technology offers a good opportunity to dissect cell heterogeneity of animals. A large number of organism-wide single-cell atlases have been constructed for vertebrates such as Homo sapiens, Macaca fascicularis, Mus musculus and Danio rerio. However, an intermediate taxon that links mammals to vertebrates of more ancient origin is still lacking. Here, we construct the first Xenopus cell landscape to date, including larval and adult organs. Common cell lineage-specific transcription factors have been identified in vertebrates, including fish, amphibians and mammals. The comparison of larval and adult erythrocytes identifies stage-specific hemoglobin subtypes, as well as a common type of cluster containing both larval and adult hemoglobin, mainly at NF59. In addition, cell lineages originating from all three layers exhibits both antigen processing and presentation during metamorphosis, indicating a common regulatory mechanism during metamorphosis. Overall, our study provides a large-scale resource for research on Xenopus metamorphosis and adult organs.

Insight into 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced disruption of zebrafish spermatogenesis via single cell RNA-seq

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a potent and environmentally persistent endocrine disrupting chemical. Our previous work demonstrated the latent reproductive maladies of early-life TCDD exposure in zebrafish. Zebrafish acutely exposed to low, environmentally relevant levels of TCDD (50 pg/mL) during two windows of sexual differentiation in development (1 hour of exposure at 3 and 7 weeks postfertilization) were later infertile, showed a reduction in sperm, and exhibited gene expression consistent with an altered microenvironment, even months after exposure. Due to the highly heterogeneous cell- type and -stage landscape of the testes, we hypothesized various cell types contribute markedly different profiles toward the pathology of TCDD exposure. To investigate the contributions of the diverse cell types in the adult zebrafish testes to TCDD-induced pathology, we utilized single-cell RNA-seq and the 10x Genomics platform. The method successfully captured every stage of testicular germ cell development. Testes of adult fish exposed during sexual differentiation to TCDD contained sharply decreased populations of late spermatocytes, spermatids, and spermatozoa. Spermatogonia and early spermatocyte populations were, in contrast, enriched following exposure. Pathway analysis of differentially expressed genes supported previous findings that TCDD exposure resulted in male infertility, and suggested this outcome is due to apoptosis of spermatids and spermatozoa, even years after exposure cessation. Increased germ cell apoptosis was confirmed histologically. These results provide support for an environmental exposure explanation of idiopathic male infertility.

RNA-Seq of Single Fish Cells – Seeking Out the Leukocytes Mediating Immunity in Teleost Fishes

The immune system is a complex and sophisticated biological system, spanning multiple levels of complexity, from the molecular level to that of tissue. Our current understanding of its function and complexity, of the heterogeneity of leukocytes, is a result of decades of concentrated efforts to delineate cellular markers using conventional methods of antibody screening and antigen identification. In mammalian models, this led to in-depth understanding of individual leukocyte subsets, their phenotypes, and their roles in health and disease. The field was further propelled forward by the development of single-cell (sc) RNA-seq technologies, offering an even broader and more integrated view of how cells work together to generate a particular response. Consequently, the adoption of scRNA-seq revealed the unexpected plasticity and heterogeneity of leukocyte populations and shifted several long-standing paradigms of immunology. This review article highlights the unprecedented opportunities offered by scRNA-seq technology to unveil the individual contributions of leukocyte subsets and their crosstalk in generating the overall immune responses in bony fishes. Single-cell transcriptomics allow identifying unseen relationships, and formulating novel hypotheses tailored for teleost species, without the need to rely on the limited number of fish-specific antibodies and pre-selected markers. Several recent studies on single-cell transcriptomes of fish have already identified previously unnoticed expression signatures and provided astonishing insights into the diversity of teleost leukocytes and the evolution of vertebrate immunity. Without a doubt, scRNA-seq in tandem with bioinformatics tools and state-of-the-art methods, will facilitate studying the teleost immune system by not only defining key markers, but also teaching us about lymphoid tissue organization, development/differentiation, cell-cell interactions, antigen receptor repertoires, states of health and disease, all across time and space in fishes. These advances will invite more researchers to develop the tools necessary to explore the immunology of fishes, which remain non-conventional animal models from which we have much to learn.