Next Article in Journal
The Involvement of Endoplasmic Reticulum Stress during the Interaction between Calcium Oxalate Crystals and Renal Tubular Epithelial Cells
Next Article in Special Issue
Translating Lupus: Comparative Transcriptional Profiles of Preclinical Lupus Models and Their Relevance to Human Disease
Previous Article in Journal
Bioactive Peptides from Corn (Zea mays L.) with the Potential to Decrease the Risk of Developing Non-Communicable Chronic Diseases: In Silico Evaluation
Previous Article in Special Issue
Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Catching the Big Fish in Big Data: A Meta-Analysis of Zebrafish Kidney scRNA-Seq Datasets Highlights Conserved Molecular Profiles of Macrophages and Neutrophils in Vertebrates

by
Aleksandr V. Bobrovskikh
1,2,*,
Ulyana S. Zubairova
2,3,
Ludmila G. Naumenko
1,2 and
Alexey V. Doroshkov
2,4
1
Department of Physics, Novosibirsk State University, 630090 Novosibirsk, Russia
2
The Federal Research Center Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, 630090 Novosibirsk, Russia
3
Department of Information Technologies, Novosibirsk State University, 630090 Novosibirsk, Russia
4
Department of Genomics and Bioinformatics, Institute of Fundamental Biology and Biotechnology, Siberian Federal University, 660036 Krasnoyarsk, Russia
*
Author to whom correspondence should be addressed.
Submission received: 13 August 2024 / Revised: 10 September 2024 / Accepted: 26 September 2024 / Published: 27 September 2024
(This article belongs to the Special Issue Immune Response Regulation in Animals)

Abstract

:

Simple Summary

This paper presents a meta-analysis of currently available single-cell RNA sequencing (scRNAseq) datasets from zebrafish kidneys. Our work aims to identify key marker genes for various immune cell types and transcription factors (TFs) involved in fish myelopoiesis. Our results validate previous studies and expand them; newly discovered markers and TFs have experimental confirmation of their roles in complementary cell types of mammals and fish species. Reconstructed gene networks of macrophage- and neutrophil-specific genes revealed the potential interactions of identified TFs and marker genes. Obtained results could be taken into account during the design of immunogenetics experiments (e.g., creating zebrafish lines with TF knockouts) and serve as a basis for novel applications for zebrafish. The revealed cell-specific markers could help in the accurate determination of fish cell types in the future. Our methodology of data integration showed its reliability and could be further adapted to meta-analyses of various scRNAseq datasets.

Abstract

The innate immune system (IIS) is an ancient and essential defense mechanism that protects animals against a wide range of pathogens and diseases. Although extensively studied in mammals, our understanding of the IIS in other taxa remains limited. The zebrafish (Danio rerio) serves as a promising model organism for investigating IIS-related processes, yet the immunogenetics of fish are not fully elucidated. To address this gap, we conducted a meta-analysis of single-cell RNA sequencing (scRNA-seq) datasets from zebrafish kidney marrow, encompassing approximately 250,000 immune cells. Our analysis confirms the presence of key genetic pathways in zebrafish innate immune cells that are similar to those identified in mammals. Zebrafish macrophages specifically express genes encoding cathepsins, major histocompatibility complex class II proteins, integral membrane proteins, and the V-ATPase complex and demonstrate the enrichment of oxidative phosphorylation ferroptosis processes. Neutrophils are characterized by the significant expression of genes encoding actins, cytoskeleton organizing proteins, the Arp2/3 complex, and glycolysis enzymes and have demonstrated their involvement in GnRH and CLR signaling pathways, adherents, and tight junctions. Both macrophages and neutrophils highly express genes of NOD-like receptors, phagosomes, and lysosome pathways and genes involved in apoptosis. Our findings reinforce the idea about the existence of a wide spectrum of immune cell phenotypes in fish since we found only a small number of cells with clear pro- or anti-inflammatory signatures.

Graphical Abstract

1. Introduction

The innate immune system (IIS) serves as the first line of defense in multicellular animals [1]. IIS cells, such as macrophages and neutrophils, not only protect organisms from pathogens but also regulate development, homeostasis, and tissue repair [2]. Among multicellular animals, Danio rerio (zebrafish) continue to garner interest as a promising model organism for studying the properties of the immune system [3]. Teleostei, including zebrafish (at both adult and larval stages), have the potential to become key organisms for IIS-related studies, particularly with the growing interest in their single-cell transcriptomics [4]. Due to their similarity to mammals and ease of handling, zebrafish are widely used as models for human diseases but remain understudied as fish themselves [5]. The overall antiquity and evolutionary conservatism of the IIS make zebrafish a potentially powerful model for studying animal pathogenesis, injuries, and wound healing; macrophages and neutrophils are known to play a crucial role in these processes [6]. In addition, there is strong evidence highlighting the crucial role of zebrafish inflammatory pathways in developmental and aging processes [7].
The fish IIS comprises various cell types such as macrophages, neutrophils, mast cells, thrombocytes, dendritic cells, B cells, and rodlet cells. These diverse immune cells are synthesized in the fish kidney, which is considered a functional ortholog of mammalian bone marrow and serves as the main hematopoietic organ, producing almost all known immune cell types and performing both innate and trained immune responses [8,9]. A persistent challenge in immune system biology is the difficulty in defining and distinguishing different blood cell populations. There are numerous immune cell subtypes that flow and morph into one another, making them inseparable in vivo and irreproducible in vitro [10,11,12]. Therefore, single-cell transcriptomics data can serve as a powerful foundation for developing novel research strategies in fish immunology. In recent years, researchers have increasingly collected diverse transcriptomic single-cell data for fish, the combined analysis of which could help to identify key markers of different immune cells that are not dependent on the specifics of individual experiments but rather are inherent in the general population of these cell types [4]. A number of single-cell transcriptomic studies have identified key marker genes in various immune cells in zebrafish, emphasizing the important roles of macrophages and neutrophils in fish immunity [9,13,14,15]. The detailed molecular profiles of macrophages, neutrophils, and B cells in the zebrafish kidney were recently investigated [8]. Additionally, cross-organ scRNAseq experiments revealed highly conserved expression profiles of immune cells [16].
Summarizing the facts above, we already possess considerable detailed knowledge about the organization of the IIS in fish thanks to single-cell transcriptomics, although the depth of study varies. However, we are still far from having a comprehensive, systematic view of fish immunogenetics compared to the mammalian one. There are several dozen experiments available for adult and immature developmental stages of zebrafish and a few experiments for other fish species [4]. Individual scRNAseq datasets undoubtedly provide valuable insights into fish immunogenetics. However, simple summarizing of marker gene lists of various immune cells from different experiments is not an optimal strategy for analyzing such large and heterogeneous data. The proper meta-analysis of scRNAseq data could provide much better results through the integration and generalization of heterogeneous cellular samples obtained from fish with different genotypes, ages, and experimental characteristics. This approach could identify crucial genetic systems expressed in various immune cell subpopulations. There are many examples of productive applications of meta-analysis in the context of single-cell transcriptomics, including the study of leukocyte diversity in atherosclerotic mouse aortas [17], the identification of SARS-CoV-2 entry genes across tissues and demographics [18], and the study of responses of tumor-reactive CXCL13+ T cells to immune-checkpoint blockade [19]. However, for zebrafish, there is only one recent scRNAseq meta-analysis dedicated to the identification of genes involved in the heart regeneration process [20].
Therefore, conducting a meta-analysis of IIS-related scRNAseq datasets for zebrafish is an urgent task. In addition to the informative and predictive power of such meta-analysis, the identified genetic systems can be compared with those of mammals to investigate the degree of evolutionary conservatism of the IIS in Vertebrata. In such complex cases, modern approaches to the integration of transcriptome and interactome data can help in generating new knowledge and hypotheses about living systems [21]. Our work aims to bridge the current gap in understanding the core molecular genetic systems of macrophages and neutrophils in zebrafish compared to mammals. We identified novel marker genes for macrophages, neutrophils, and other cells of the IIS in zebrafish and transcription factors involved in myelopoiesis in fish. The discovered genes and pathways of fish neutrophils and macrophages require further study in the field of immunogenetics and transgenesis to evaluate their functions in pathogenesis, regeneration, and aging. Thus, our research could contribute to the development of novel applications of zebrafish as a model of human immunity. In addition, our work methodologically complements existing strategies for the meta-analysis of scRNAseq data. Our analysis strategy, which consisted of initial soft filtering of the data followed by more stringent filtering, proved to be flexible and effective, providing a diversity of cell types and high read coverage in the final dataset. We think that in future meta-analyses, researchers should select the threshold for filtering cells by the number of reads based on the quality of the original data rather than blindly using certain thresholds in data analysis pipelines. Therefore, the general steps of our meta-analysis can be used in future studies devoted to the joined analysis of scRNAseq data.

2. Materials and Methods

The overall pipeline of carried meta-analysis and post-analysis is presented in Figure 1. The main steps are described below, and in separate Section 2.1, Section 2.2, Section 2.3, Section 2.4 and Section 2.5, we explain them in detail.
1.
For our meta-analysis, we searched and downloaded the suitable scRNAseq datasets described in our previous study [4] and experiments that were added to GEO NCBI during the last year;
2.
The downloaded datasets were converted to zebrafish gene identifiers (GRCz11). For each experimental dataset, we made a separate dataframe used to create the corresponding Seurat object. Seurat objects were annotated according to the general information for the experiments (GEO identifiers, fish age, and genotype);
3.
Individual Seurat objects went through a soft initial filtering procedure and were merged into a single Seurat. Afterward, to correct batch effects, we used normalization and harmonization procedures. In the resulting space, the main types of immune and non-immune cells of the zebrafish kidney were identified;
4.
Further filtering involved the removal of non-immune cells, immune cells with insufficient coverage, and the removal of genes not specific to immune cells. One experiment was also excluded from the final dataset because of the strong individual batch effect, which cannot be completely eliminated by Harmony and may lead to incorrect data integration;
5.
The final dataset was normalized and harmonized. In the resulting clusters, we identified immune cell types based on literature data, their marker genes, and TFs;
6.
The identified marker genes of different cell types were ranked and classified into one of four groups based on their expression level and the strength of upregulation in the target cell type;
7.
Next, for the most highly specific and highly expressed genes of macrophages and neutrophils, co-expression (Step 7 in Figure 1) and protein–protein interactions (Step 8 in Figure 1) were extracted, which allowed us to reconstruct gene networks relevant to the immune cells of these two types (Step 9 in Figure 1).

2.1. Data Collection

Following the search methodology from our earlier study [4], we selected relevant datasets that contained at least 5000 cells from zebrafish kidneys to ensure comparability in the meta-analysis. Thirteen experiments from the GEO NCBI database [28] met the inclusion criteria and were analyzed further: GSE100910 [29], GSE112438 [30], GSE130487 [31], GSE150373 and GSE176036 [32], GSE151231 [33], GSE166646 [34], GSE179401 [35], GSE190794 [36], GSE191029 [37], GSE242133 [8], GSE246039 [38], and GSE252788 [9]. Detailed information about these experiments is given in Table 1 and Supplementary Table S6.
For the further meta-analysis, we downloaded and used the count matrices or Seurat objects provided by the authors in the Supplementary Materials in the corresponding GEO NCBI records (available by their GEO NCBI identifiers).

2.2. Data Meta-Analysis

2.2.1. Software and Packages Used

All data analysis steps were performed using the R programming language v. 4.4.0 and the RStudio desktop environment v. 2023.12.1. The following packages were used: Seurat v. 5 [23], Harmony [24], Genekitr [22], hdWGCNA [26], tidyverse [39], dplyr, readr, ggplot2, cowplot, and patchwork. Bar plots and heatmaps were visualized using Python v. 3.11 and the Seaborn (function seaborn.heatmap) and Matplotlib (function matplotlib.pyplot.bar) packages. A Sankey diagram was created using the SankeyMATIC web service.

2.2.2. Dataset Integration and Cell Type Identification

The initial step of our analysis involved merging individual sample count matrices from the selected experiments. Due to variations in gene identifiers across the analyzed experiments, the TransId function from the Genekitr package was used to convert all identifiers to an Ensembl format. Subsequently, individual count matrices were used to create Seurat objects via the CreateSeuratObject function from the Seurat package, applying soft filtering parameters: min.cells = 10, min.features = 200, and mitochondrial gene expression less than 10%. Multiple Seurat objects were then combined into a unified Seurat object using the merge function. Normalization and variance stabilization were performed using SCTransform v. 2. To initially estimate the number of clusters prior to harmonization, principal components were calculated (RunPCA function, npcs = 30, verbose = F), followed by the calculation of UMAP projections of cells (RunUMAP function, reduction = “pca”, dims = 1:30, verbose = FALSE). Preliminary clusters (prior to harmonization) were identified using the FindNeighbors (reduction = “pca”, dims = 1:30, verbose = F) and FindClusters (resolution = 1.0, verbose = FALSE) functions.
Subsequently, we used the Harmony algorithm to eliminate batch effects from individual datasets. The RunHarmony function was applied to the SCT assay of the unified Seurat object using the following parameters:
  • group.by.vars = "orig.ident", project.dim = TRUE, plot_convergence = TRUE, verbose = TRUE, assay.use = "SCT".
After harmonization, Harmony components were utilized for uniform manifold approximation and projection (UMAP) calculation and cluster identification: RunUMAP, reduction = “harmony”, dims = 1:23. Clusters were then identified using the FindNeighbors (reduction = “harmony”, dims = 1:23) and FindClusters (resolution = 1, verbose = TRUE) functions. The cell types of these clusters were annotated based on marker genes from [8].

2.2.3. Analysis of the Immune Cell Subset and Identification of Marker Genes

The subsequent analysis aimed to identify immune cell-specific genes that are highly expressed within immune cell types. To calculate immune-related marker genes, we used the FindMarkers function with the following parameters: only.pos = TRUE (upregulated), min.pct = 0.02 (expressed at least in 2 percent cells of the immune population), and logfc.threshold = 0.5 (upregulated by log2FC ≥ 0.5 in the immune population). Afterward, cells from non-immune clusters were excluded from the Seurat object. Further filtration was used to remove immune cells with low coverage. In particular, we selected immune cells with a high number of counts (3000–30,000 per cell), which expressed 1000–4000 genes at a non-zero level using the subset function. Further calculations and analyses were conducted on this most informative and immune cell-specific subspace (matrix consisting of selected immune cells × immune-related marker genes). Subsequently, we repeated the harmonization process on the immune subset using the RunHarmony function. The first 26 harmony dimensions were used for UMAP projection (RunUMAP) and cluster identification (FindNeighbors; FindClusters). Marker genes of individual clusters were calculated using the functions PrepSCTFindMarkers and FindAllMarkers: only.pos = TRUE (upregulated), min.pct = 0.25 (expressed at least in 25 percent cells of target clusters), and logfc.threshold = 0.5 (upregulated by log2FC ≥ 0.5 in target clusters compared to all other clusters).
To precisely investigate the cell types of the revealed immune cell subset, we compared expression within all Seurat clusters to assess their similarity. In particular, the function AverageExpression was used on an immune cell subset to calculate the expression matrix of all Seurat clusters in average counts-per-cell (ACPC) values (matrix of Seurat clusters × ACPC values of immune-related marker genes). This matrix was used to calculate the correlation matrix using the pandas.DataFrame.corr function. Based on the correlation matrix, we calculated the distance matrix (1 − abs(correlation)). For hierarchical clustering, the square form of the distance matrix was given to the function scipy.cluster.hierarchy.linkage with the following parameters: method = ‘ward’, metric = ‘Euclidean’, and optimal ordering = ‘True’. The majority of clusters were classified based on the expression of their marker genes [8] into neutrophils, macrophages, B cells, and hematopoietic stem cells. The minor cell types that we also detected were macrophage-like cells, neutrophil-like cells, immature neutrophils, monocytes, and T cells. After initialization of all cell types, cell-type-specific markers were calculated using the following functions FindMarkers: ident.1 = c(target cluster numbers), only.pos = TRUE (upregulated), min.pct = 0.25 (expressed at least in 25 percent cells of the target cell type), and logfc.threshold = 0.5 (upregulated by log2FC ≥ 0.5 in target cell type); ident.1 were different vectors corresponding to identified cell types that were compared with the rest of the immune clusters. Once we had identified markers in the cell type cluster groups, we identified transcription factors (TFs) among them using a list of transcription factors known for zebrafish in the database AnimalTFDB4 [25].

2.2.4. Classification of Cell-Type-Specific Markers

All cell-type cluster markers were classified into four groups. The first group consisted of cell-type-specific genes with log2FC ≥ 2 and average counts-per-cell (ACPC) ≥ 2. This group included the most specific markers whose expression could be reliably detected in small datasets of immune cells. The second group included genes that were highly expressed in a given cell type with ACPC ≥ 2 but were less specific, with 0.5 ≤ log2FC ≤ 2. The third group consisted of cell-type-specific genes with log2FC ≥ 2 but low expression levels with ACPC ≤ 2. These cannot be fully considered marker genes, but their expression may be inducible under specific conditions and could play a role in the immune status of the cell. The fourth group comprised the remaining differential genes with low basal expression, where ACPC was ≤2 and with 0.5 ≤ log2FC ≤ 2.

2.3. Identification of Human Orthologs of Fish Immune Genes

All gene identifiers of detected marker genes from Groups 1–4 were used to search for corresponding orthological genes in the H. sapiens genome. The Ensembl Biomart data mining tool [40] was used, and lists of zebrafish genes were provided as input sources. The following search options were used: Database: Ensembl Genes 112; Dataset: Zebrafish Genes (GRCz11); Filters: Gene: Input external references ID list: Gene stable ID(s); Attributes: Homologues: Orthologues (F–J): Human Orthologues: Human gene name, human orthology confidence.

2.4. Reconstruction of Highly Expressed Macrophage and Neutrophil Gene Networks

The identified genes from Group 1 and Group 2 for macrophages and neutrophils were used to reconstruct cell-type-specific gene networks. The STRING-db database v. 12.0. was used as a source of interactions of studied protein-encoding genes [41]. Combined lists of macrophage-specific and neutrophil-specific genes were submitted to the web server version of the database using the option Search Multiple Proteins by Names/Identifiers. For the reconstruction of networks, we used the following parameters: minimum required interaction score: 0.4 (medium); sources of active interactions: experiments and databases; network display options: show your query protein names. The resulting protein–protein interaction networks were saved as a short tabular text output.
To search for additional associations between immune genes in zebrafish, we performed co-expression analysis on the final dataset of immune cells using the hdWGCNA [26] package. The analysis was performed according to a basic tutorial vignette, and its main steps consisted of setting up a Seurat object for hdWGCNA, the construction of metacells, and the construction of a co-expression network. After reconstruction of the co-expression matrix, we saved the topological overlap matrix and used the exportNetworkToCytoscape option of this matrix to export the resulting list of all edges in a Cytoscape-usable format. From the full co-expression edge list, co-expression interactions for macrophage-specific and neutrophil-specific genes with a co-expression threshold of at least 0.1 were extracted for further analysis.
Then, once we had data on protein–protein and co-expression interactions of neutrophil and macrophage genes, these edgelists were combined and used as input into the Cytoscape environment v. 3.10.2 [27] using option File–Import–Network from File. When importing edgelists, additional attributes were specified for interactions: co-expression or protein–protein interaction. Additional meta-information, such as gene expression levels, their annotations, and belonging to a certain group, was imported into the networks using the File–Import–Table from the File option. The layout of the networks was carried out according to the following principles: transcription factors were placed at the top of the network, known and candidate marker genes and their interactors were placed below TFs, and the remaining genes were placed below marker genes. Node degrees were estimated using the built-in NetworkAnalyzer tool. The final network layouts were saved and exported from the Cytoscape environment to Figures.

2.5. Functional Annotation

For gene-set enrichment analysis of genes of interest, the ShinyGO v. 0.8.0 tool was used. The following options were used: Species: Zebrafish genes GRCz11; FDR cut-off = 0.05; pathway database: KEGG; minimal pathway size = 2; remove redundancy: yes; abbreviate pathways: yes.

3. Results and Discussion

3.1. Zebrafish Kidney Marrow scRNA-Seq Dataset Integration and IIS Cell Classification

Based on the search methodology described in our prior study [4], we selected 13 relevant datasets consisting of at least 5000 cells from zebrafish kidneys [8,9,29,30,31,32,33,34,35,36,37,38]. The main steps of our initial analysis, filtration procedures, and cellular compositions of the final subset are presented in Figure 2. After applying soft filtration to 13 datasets, we obtained a total of 382,295 cells, which were subsequently normalized, harmonized, and clustered into 65 distinct cell groups (UMAP is shown in Figure 2B). Based on known gene markers for non-immune and immune cell types in the zebrafish kidney described in [8], we determined that 38 clusters correspond to immune cells, including 254,256 cells and representing five main cell types: hematopoietic stem cells (HSCs), B cells, macrophages, T cells, and neutrophils (showed as non-black-colored clusters in Figure 2B). The remaining 27 clusters correspond to non-immune cells, including 128,039 cells, and represent six main cell types: non-immune HSCs, erythrocyte progenitors, kidney mucin cells, distal tubule cells, proximal tubule cells, and multiciliated cells (showed as black-colored clusters in Figure 2B). The expression plots of the key marker genes for both immune and non-immune cells are provided in Supplementary Figure S1.
To ensure a more precise and accurate final classification of immune cells, we filtered the final dataset in four steps. The first step involved filtering out cells from non-immune clusters. Afterward, we left in the resulting Seurat object only 3737 genes upregulated in the immune cells compared to non-immune ones (log2FC 0.5 ; see Supplementary Table S5.1). Then, due to the sufficient number of immune cells (254,256 cells in total), more stringent filtering was performed prior to the final data analysis. Specifically, cells with the best available read coverage, ranging from 3000 to 30,000 read counts per cell and expressing between 400 and 4000 unique genes per cell, were selected. The filtered dataset was composed of 96,509 immune cells. The different colors in Figure 2B) highlighted the presence of immune cells from initial clusters in this subset. During the subsequent step of harmonization and clustering of this subset, it was noted that the selected immune cells from the GSE112438 experiment revealed significant differences in embeddings compared to the cells from other experiments. This led to the exclusion of this experiment from the final analysis (see Supplementary Figure S2). Consequently, an additional 2111 cells were excluded, resulting in a final dataset comprising 94,398 cells from 12 experiments and 3737 relevant immune genes (the expression matrix of the final dataset is provided in Supplementary File S12). The cellular composition of the final dataset, resulting from the further integration of individual scRNA-seq experiments, is detailed in Figure 2C.
To sum up, our filtering approach allowed us to effectively select the highest quality cells available, which subsequently facilitated the accurate classification of the specific molecular genetic systems of individual cell types. This process reduced the representation of several cell clusters in the final dataset, marked in gray in Figure 2B, which appear to be in a transitional state among stem and mature immune cells, most likely neutrophils and T cells. Nevertheless, our final immune dataset consists of all terminal types characterized by pronounced genetic signatures in satisfactory quantities, marked in pink, or large quantities, marked in red or dark red, relative to the original dataset. Consequently, we maintained the balance between the qualitative composition and read coverage of individual cells. This filtering procedure suggests that researchers need to reasonably select thresholds when filtering scRNA-seq data, as these decisions can greatly affect the quality and results of the final analysis. Using the known markers [8], nine immune cell types were identified: hematopoietic stem cells, macrophages, neutrophils, B cells, T cells, monocytes, immature neutrophils, neutrophil-like cells, and macrophage-like cells (Figure 2C). The selected experiments are diverse in their design, including different mutant genotypes, various cell isolation and sorting protocols, and fish of different ages, as detailed in Supplementary Table S6. Therefore, it is not surprising that these experiments were quite heterogeneous in their cellular composition. Due to the inclusion of a sufficient number of experiments in our analysis, it was possible to compile a representative dataset balanced in the composition of individual cell types.
The overall information about the properties of the final dataset is presented in Figure 3. In particular, we used best-quality immune cells marked in red in Figure 3A. Classification of these immune cells revealed 44 clusters shown in Figure 3B, and corresponding immune cell types were revealed based on marker genes provided in [8]; the expression of certain important markers is shown in Figure 3D. The expression levels of the top 10 marker genes for each cell type are provided in Supplementary Figure S3. To achieve a more precise definition of the cell types, hierarchical clustering was performed; the results are depicted in Figure 3C.
Thus, we revealed immune markers of the main IIS cell types and compared them with the study [8]. Five primary cell populations were clearly identified in our dataset: neutrophils, macrophages, hematopoietic stem cells, B cells, and immature neutrophils. Additionally, based on specific known markers, UMAP clustering, and hierarchical clustering, monocytes, immature T cells, macrophage-like, and neutrophil-like cells were identified. These cell types, their corresponding clusters, the number of cells in individual clusters, and the total number of cells by cell type are presented in Table 2. Detailed information about cell type markers and their classification is provided in Supplementary Table S4. In the sections below, we provide comprehensive information on the most relevant immune cell markers identified in our analysis.

3.1.1. Hematopoietic Stem Cells

Hematopoietic stem cells are located in the top right area of the UMAP projection (Figure 3B, Table 2). HSCs are characterized by the increased expression of numerous ribosomal genes, including 17 rpl and 10 rps genes. Additionally, well-known marker genes such as foxp4, rsl1d1, csf1rb, gtpbp4, rps7, and rpl7 exhibit strong upregulation in the cells of HSC clusters (Supplementary Figure S3).

3.1.2. Macrophages and Macrophage-like Cells

Macrophages are located in the middle right area of the UMAP projection (Figure 3B, Table 2). The expression of known macrophage markers, including c1qc, c1qa, marco, c1qb, lgals3bp.1, hmox1a, fgl2a, grna, grna.1, ctss2.2, fthl27, ccl34a.4, s100a10b, and ctsd, was confirmed. Macrophage-specific genes were strongly upregulated in clusters 2 and 4. Additionally, clusters 10, 32, and 40 showed significant similarity to the main clusters and were also identified as macrophages (Figure 3C).
Macrophage-like cells were identified as a highly heterogeneous group of small clusters occupying the bottom right area of the UMAP projection (Figure 3B). These cells exhibit high expression of several known macrophage markers, including grna.1, grna.2, ccl34a.4, lgmn, s100a10b, ctss2.2, and ctsd, but not classical markers of mature macrophages such as marco, c1qa, c1qc, and c1qb.

3.1.3. Mature Neutrophils, Immature Neutrophils, and Neutrophil-like Cells

Neutrophils form a large group of cells located in the middle left area of the UMAP projection (Figure 3B). These cells are characterized by increased expression of markers such as lyz, cpa5, sms, lect2.1, npsn, mpx, nccrp1, ncf1, il1b, and scpp8. This group of cells is distinctly different from stem cells and other immune subtypes. Based on clustering (Figure 3C) and the expression of specific markers, such as abcc13, hsd3b7, and si:ch1073-429i10.1 (Supplementary Figure S3), ten clusters were classified as mature neutrophils (see Table 2). Additionally, two clusters, which are positioned closer to T cells and stem cells on the UMAP projection (Figure 3B), were identified as immature neutrophils. Among the highly expressed markers of immature neutrophils are well-known genes such as ncf1, il1b, mmp9, scpp8, lect2.1, npsn, mpx, and nccrp1. Compared to mature neutrophils, these cells additionally express genes such as rgs2 and gpr84 but do not express key markers of mature neutrophils, including lyz, cpa5, and sms. We also identified neutrophil-like cells in the bottom right of the UMAP projection that express certain neutrophil-specific genes, including lect2.1, npsn, cpa5, mpx, and sms, but at lower levels than classical neutrophils (Supplementary Figure S3).

3.1.4. B Cells, T Cells, and Monocytes

B cells occupy a compact area in the upper left of the UMAP projection (Figure 3B). They express many known markers, such as igic1s1, plaat1l, cxcr4a, bhlhe40, bcl2l10, irf8, ccr9a, and cd74a. Interestingly, these cells also exhibit high levels of several rpl and rps genes, similar to stem cells, indicating the occurrence of translation processes within this cell population.
T cells were identified by the expression of the runx3 gene and the lymphocyte-specific gene ccr9a. These cells are most likely immature T cells, as they represent a minor component of kidney immune cells and, like B cells, express numerous genes related to translation processes.
Monocytes were identified based on the primary cell marker cd34, which corresponds to the monocytic lineage of progenitors [42]. This extensive group of cells is located between stem cells and macrophages, retaining some progenitor properties and characterized by the elevated expression of six rpl genes and five rps genes. The corresponding plots are provided in Supplementary Figure S3.

3.2. Detailed Analysis of the Molecular Genetic Systems of Immune Cells

During our analysis, we revealed that macrophages and neutrophils are the two most terminal and distinctive cell types presented in our datasets, and they are enriched in the expression of many specific immune genes. We also revealed minor subpopulations of B cells and immature T cells, confirming the hypothesis that the zebrafish kidney could serve as a potential source of lymphocytes [8]. The other cell types either show similarities to macrophages and neutrophils or exhibit significant stemness and do not express specific immune genes as strongly. Hence, in the following sections, we primarily focused on the in-depth analysis of the molecular–genetic systems of macrophages and neutrophils.
For further analysis, we classified immune markers into four groups by their average expression in identified cell types and by their upregulation folds. This classification is described in Materials and Methods, Section 2.2.4. In short, genes in Group 1 (cell-type-specific highly expressed genes: log2FC ≥ 2; ACPC ≥ 2) and Group 2 (highly expressed genes: ACPC ≥ 2; 0.5 ≤ log2FC ≤ 2) should be considered the highest priority for further experimental and theoretical research. The roles of the markers from Group 3 (cell-type-specific but weakly expressed) and Group 4 (less specific and weakly expressed) require additional meta-analyses or experiments to verify their involvement in myelopoiesis processes in zebrafish. The expression heatmap of identified cell-type-specific highly expressed genes (Group 1) in all Seurat clusters is depicted in Figure 4.
We identified 62 highly specific marker genes for macrophages and macrophage-like cells (see Figure 4A), thirteen of which were previously highlighted in the list of the 20 most specific macrophage marker genes [8]. For mature neutrophils and immature neutrophils, 58 specific genes were revealed, depicted in Figure 4B; among them were thirteen of the most relevant neutrophil genes mentioned previously [8]. B-cells and HSCs can be identified by their specific twelve and five markers, respectively, as shown in Figure 4C. Seven B-cell markers were previously reported by Hu et al., as well as two markers of HSCs, which were validated in our analysis [8]. Therefore, our empirical classification is reliable since we used rather stringent thresholds for both fold-changes and basal expression levels for selecting these genes. Furthermore, our classification and selection of these genes is strongly supported by previous data [8] and expands key marker genes for neutrophils and macrophages in the zebrafish kidney.
Thus, the number of newly discovered key markers of immune cells in the zebrafish kidney does not exceed one hundred genes, most of which are expressed in macrophages and neutrophils; all marker genes and their classification are given in Supplementary Table S4. Therefore, in Section 3.2.1, Section 3.2.2 and Section 3.2.3, we discuss in detail the potential roles of our discovered genes in fish immune cells based on the available literature on different fish species, mice, and humans. In Section 3.2.4, we review and discuss, in the context of recent discoveries, our newly discovered transcription factors involved in myelopoiesis in zebrafish. Throughout Section 3.2.1, Section 3.2.2, Section 3.2.3 and Section 3.2.4, we used the generally accepted gene nomenclature spelling: human gene names are marked by upper case (NAME), mouse gene names are marked by title case (Name), and fish gene names are marked by lower case (name).

3.2.1. The Novel Marker Genes of Zebrafish Macrophages

For fish macrophages (Figure 4A), our meta-analysis verified the expression of 13 out of the 20 primary kidney marrow macrophage marker genes from the study [8]. Additionally, the expression of 41 newly identified marker genes was confirmed, including six cathepsin genes (ctsa, ctsba, ctsc, ctsh, ctsk, and ctsz), three annexin genes (anxa1a, anxa3b, and anxa5b), two cd74 genes (cd74a and cd74b), two keratin genes (krt4 and krt18b), two solute carrier family genes (slc2a6 and slc22a21), and two major histocompatibility complex class II genes (mhc2a and mhc2dab). Cathepsins are essential as a defense against intracellular bacteria. The upregulation of cathepsins CTSB, CTSC, CTSH, and CTSZ is a marker of the M1 phenotype of human macrophages, while the upregulation of CTSK is more pronounced in M2 macrophages [43]. In mice, annexin Anxa1 has been shown to induce macrophage migration to participate in muscle regeneration through the AMPK cascade [44], annexin Anxa3 upregulates the infiltrated neutrophil-lymphocyte ratio to remodel the immune microenvironment in hepatocellular carcinoma [45], and annexin Anxa5 regulates the polarization of liver macrophages by acting on its primary target, pyruvate kinase M2 [46]. Increased CD74 gene expression in atherosclerosis is an indicator of macrophage inflammatory responses [47]. A cellular model of head and neck squamous cell carcinoma revealed the high correlation between the expression levels of keratin genes (KRT4, KRT24, KRT78, and KRT13) and the infiltration levels of CD8+ T cells and macrophages [48]. One of the solute carrier family proteins in mice, Slc37a2, has been shown to suppress macrophage inflammation through the regulation of glycolysis [49]. Among the newly identified macrophage marker genes, four genes encoding receptors were found: havcr1 (hepatitis A virus cellular receptor 1), mrc1b (mannose receptor C type 1b), which is upregulated in M2 macrophages in fish [50], prosaposin psap, and pglyrp5 (peptidoglycan recognition protein 5), which has been shown to play a role in antimicrobial activity and enhancing macrophage phagocytosis in Larimichthys crocea [51]. Additionally, four transcription factors were identified: junba, which is upregulated during zebrafish heart regeneration [52], mafba, known for its role in regulating microglia colonization in zebrafish [53], tcima, recently identified as a macrophage-enriched gene in zebrafish [54], and id1, previously identified as relevant to NK cells in zebrafish [15]. Among other relevant genes, the AHNAK nucleoprotein has been observed to suppress tumor proliferation [55], while Cfp (complement factor properdin) knockout mice exhibited a shift towards M2 macrophage polarization [56] and galectin lgals2a, whose ortholog in an in vivo mouse study demonstrated an increase in macrophage numbers [57]. Additionally, we detected the strong expression of lysozyme g lygl1, whose expression in fish immune cells was confirmed in a study [58], indicating membrane-spanning 4-domains member ms4a17a.10, previously shown to be expressed in fish B cells [15], surfactant protein Bb sftpbb, which was confirmed to be upregulated in zebrafish macrophages [59], syndecan 4 (sdc4), whose ortholog knockdown in mice (Sdc4) leads to a pro-inflammatory phenotype in macrophages and the development of atherosclerosis [60], vacuole membrane protein 1 (vmp1), whose ortholog knockout (VMP1) in human cell lines triggered the release of pro-inflammatory cytokines [61], and gene zgc:174904, whose elevated expression was previously detected in fish B cells [62]. The roles of the genes discussed above are confirmed to varying degrees in macrophages and, in some cases, in other immune cells of zebrafish or only for mammals. However, for five of the newly identified genes, there is currently no information about their specific roles in macrophages: brain-subtype creatine kinase ckbb, chloride intracellular channel protein clic2, phorbol-12-myristate-13-acetate-induced protein 1 pmaip1, tetraspanin-36 (tspan36), and tubulin polymerization promoting protein (LOC100333482).
Macrophage-like cells (Figure 4A) express a set of genes typical for the classical macrophage population, including lgals2a, lygl1, cd74a, and cd74b, albeit at lower levels. In addition to the known markers grna.1, grna.2, cxcl12a, and lgmn, we identified higher expression of new potential markers in these cells, such as beta-microseminoprotein-like (LOC110437731), mtp family protein (LOC100151049), uncharacterized protein LOC101883708, and si:ch211-194m7.4. The beta-microseminoprotein-like protein has demonstrated antibacterial activity in chicken eggs [63], while the other genes remain largely uncharacterized.

3.2.2. The Novel Marker Genes of Zebrafish Neutrophils and Immature Neutrophils

The expression of genes typical for mature neutrophils is illustrated in the top part of Figure 4B. Newly identified proteins include two ATP-binding cassette proteins (abcb9 and abcc13), cathepsin S (ctss1), 15-hydroxyprostaglandin dehydrogenase (hpgd), 3′-phosphoadenosine 5′-phosphosulfate synthase 2b (papss2b), arachidonate 5-lipoxygenase-activating protein (alox5ap), cholesterol 25-hydroxylase-like protein 2 (ch25hl2), hydroxy-delta-5-steroid dehydrogenase (hsd3b7), phospholipid phosphatase 1 (plpp1a), plac8 onzin-related protein 6 (ponzr6), protein kinase C delta type (prkcda), S-adenosylmethionine decarboxylase proenzyme (amd1), serglycin (srgn), sulfotransferase (sult2st1), thy-1 cell surface antigen (thy1), very-long-chain (3R)-3-hydroxyacyl-CoA dehydratase (hacd4), vesicle-associated membrane protein 8 (vamp8), and a gene with predicted oxidoreductase activity si:ch1073-429i10.1. The expression of abcb9 has been validated in zebrafish neutrophils [64]. Studies using mouse models have shown that knockout of cathepsin Ctss reduces neutrophil recruitment [65]. Inhibition of the Hpgd gene has been found to enhance neutrophil regeneration in mice [66], indicating its potential as a marker for neutrophil maturation. The ALOX5AP gene plays a role in neutrophil ferroptosis and compromises cell wall integrity during intracranial aneurysms [67]. In humans, HSD3B7 gene expression is positively correlated with the infiltration of immune cells, including neutrophils [68]. The ponzr6 gene has been previously confirmed in zebrafish neutrophils [69], while the study [70] demonstrated the high specificity of prkcda expression in fish neutrophils. The role of amd1 in neutrophil maturation is noted in the research [9]. In tumor tissues, the suppression of SRGN gene expression correlated with lower immune infiltration of neutrophils and other immune cells [71]. In mice, cholesterol sulfate produced by the sulfotransferase Sult2b1 acts on inflammatory neutrophils and prevents excessive intestinal inflammation [72]. Thy-1 (CD90) is involved in the skin wound healing mechanism by activating various immune cell types in mammals [73]. The role of the Hacd4 gene in inflammatory processes has been confirmed in mouse models [74]. Different Vamp genes in mice are engaged in neutrophil exocytosis, suggesting a similar function for the vamp8 gene in fish [75]. The functions of the identified genes abcc13, ch25hl2, and plpp1a in zebrafish neutrophils need to be investigated in future research.
Similar to mature neutrophils, immature neutrophils (Figure 4B, middle section) highly express several established markers, such as lyz, cpa5, lect2.1, npsn, mpx, and srgn. We identified a distinct set of genes specific to immature neutrophils: disintegrin and metalloproteinase domain 8a (adam8a), activating transcription factor 3 (atf3), arginase (arg2), immune-related lectin-like receptor 4 (illr4), interleukin 6 receptor (il6r), MCL1 apoptosis regulator (mcl1b), ornithine decarboxylase 1 (odc1), purine nucleoside phosphorylase (pnp5a), a member of the RAS oncogene family rab44, scinderin-like b (scinlb), serine/threonine kinase 17B (stk17b), solute carrier family 7 member 3a (slc7a3a), sphingosine-1-phosphate receptor 4 (s1pr4), suppressor of cytokine-signaling 3a (socs3a), TIMP metallopeptidase inhibitor 2b (timp2b), transcobalamin beta b (tcnbb), tumor necrosis factor b (tnfb), v-set immunoregulatory receptor (vsir), zgc:158343, si:ch211-264f5.2, and BI1 family proteins si:ch211-284o19.8, si:ch73-248e21.5, si:ch73-343l4.8, si:dkey-195m11.11, si:dkey-27h10.2, and zgc:158343. The deletion or inhibition of the Adam8 gene leads to impaired neutrophil transmigration in mice [76]. Atf3-deficient mice show weakened antibacterial immunity due to their involvement in the modulation of macrophage killing and migration functions [77]. The ARG2 gene is a signature of N2 neutrophils in a rosacea-like model and is induced by LPS+ IFN- γ and IL-4 [78]. The illr4 gene has been identified as a chemotaxis gene for neutrophils in fish [79]. It is known that human neutrophils produce interleukin 6 and its receptor IL6R during the inflammatory response [80]. While the anti-apoptotic gene mcl1b was previously recognized for its role in neutrophil maturation, our analysis indicates a more stable expression in immature neutrophil clusters [9]. Deletion of the Odc1 gene in mice resulted in delayed neutrophil clearance and a reduction in the proresolving cytokine interleukin-10 [81]. The fibrinogen and endopeptidase activity of the pnp5a gene in fish during the inflammation and immune response has been confirmed [82]. Rab44 gene expression has been shown to decrease during neutrophil maturation in the mouse model [83], a finding that our study supports for its ortholog in zebrafish. The proteomic analysis of zebrafish kidney neutrophils showed a significant increase in the expression of several proteins, including scinlb, in response to chemically-induced inflammation [84]. The Stk17b gene in mice and STK17B in humans have been identified as signatures of activated neutrophils in metastatic tumors [85]. The slc7a3a gene regulates lipid metabolism through the mediation of arginine-dependent nitric oxide synthesis in zebrafish cells [86], and nitric oxide subsequently modulates numerous physiological properties of neutrophils in response to various pathogenic stimuli, which may be a crucial feature of immature neutrophils, reflecting their functional adaptability [87]. The S1PR4 gene in humans and the S1pr4 gene in mice are known participants in immunity modulation, the regulation of homeostasis, and neutrophil migration [88]. In mice, the Socs3 gene serves as a negative regulator of Tbk1 and the interferon pathway and is expressed in neutrophils in response to the granulocyte–macrophage colony-stimulating factor combined with Prostaglandin E2, supporting the anti-inflammatory phenotype of neutrophils [89] and potentially maintaining their immature characteristics. The expression of timp2b could indicate its critical role in the development of local inflammation, resolution, and tissue repair in zebrafish, which is consistent with the revealed role of its ortholog TIMP-2 in human leukocytes [90]. The tcnbb gene’s expression in neutrophils was confirmed in [8]. It is well-known that an injection of tnf-alpha precursor in zebrafish results in induced neutrophil recruitment [91], and tnfb may similarly stimulate immature neutrophils. Fish orthologs of the vsir gene have been confirmed in neutrophils and other immune cells in both humans and mice [92]. Huang’s dissertation [93] reports on the upregulation of the zgc:158343 gene during bacterial infection in zebrafish, and the upregulation of the si:ch211-284o19.8 gene was revealed in the fish’s response to bacteria [94].
Neutrophil-like cells (Figure 4B, lower section) exhibit high expression levels of several neutrophil-related genes, including lect2.1, npsn, mpx, srgn, ncf1, scinlb, and illr4. The genes specific to these cells are capzb, LOC103910140, si:ch211-9d9.1, and si:ch211-223l2.4. The upregulation of CAPZB gene expression has been detected in human B cells in the context of lymphoma [95]. The expression of the glycoprotein-like gene cd59 (LOC103910140) in fish is downregulated in response to alginate oligosaccharides [96], which similarly leads to reduced neutrophil infiltration and lower levels of inflammatory markers in mice [97].

3.2.3. The Novel Markers of Other Cell Types

In B cells (Figure 4C, upper section), we detected the expression of dusp2, ccl35.1, syk, grap2b, and kctd5a. The Dusp genes, including Dusp2, function as regulators and specifiers of innate immune responses in mice, with their expression being controlled via the MAPK pathway following the binding of foreign agents to pattern recognition receptors [98]. High expression of ccl35.1 has been identified in zebrafish dendritic cells [16]. The Syk gene has been shown to play an essential role in B cell activation in mice [99]. The downregulation of the Grap2 gene was observed in the transcriptome of B cells from pregnant mice [100]. In the tumor microenvironment, KCTD5 gene expression was positively correlated with macrophages and T cells while negatively correlated with plasma cells, mast cells, and memory B cells [101].
Increased expression of the cd34 marker gene was observed in some monocytic clusters and HSCs (Figure 4C, middle section). Sialomucin (CD34) is known as a marker for the early maturity stages of human monocytic precursors, suggesting that the monocytic cells detected in our dataset may possess a high degree of immaturity and stemness [102].
We found that in HSCs (Figure 4C, lower section), in addition to the well-known csf1rb and nanos1, the genes hist2h2l, myca, and si:ch211-161c3.6 are highly expressed. During Aeromonas salmonicida infection in Cyprinus carpio, the expression of hist2h2l significantly decreases, which may indicate the mobilization of stem cell differentiation in the immune system [103]. The role of the cell cycle regulator myca in the proliferation of stem cells in zebrafish was demonstrated in [104]. The transcription factor si:ch211-161c3.6 was previously identified as a marker of bone marrow stem cells in fish [36].

3.2.4. The Novel Transcription Factors Involved in Zebrafish Myelopoesis

Some of the specific genes identified as markers for immune cell types in our study, including several transcription factors (see Figure 5), are discussed in the context of myelopoiesis in fish across multiple studies. The review [105] covers various genes that regulate neutrophil homeostasis, including transcription factors cebpa, cebpb, stat1a, myca, and spi1b, which were confirmed in our analysis. In the study of the transcriptional network of the zbtb14 factor, which regulates the development of monocytes and macrophages, factors cebpa, myca, spi1b, irf8, and mafba were found to be involved [106]. The genes spi1a, mafba, mafbb, irf8, and hif1ab identified in our analysis are potential immunotherapeutic targets for myeloid cells involved in neurovascular repair following ischemic stroke [107]. Transcription factors klf6a, cebpb, cebpa, hif1ab, and spi1b are highlighted as key players in granulopoiesis in fish [108]. These eleven factors are represented in Figure 5. Additionally, we identified 27 transcription factors to be involved in fish myelopoiesis and discussed them below.
According to our data, three of the factors broadly mentioned in the literature, cebpa, cebpb, and spi1b, are common to most types of immune cells. We identified that transcription factors ybx1, zgc:162730, xbp1, and fosab are also ubiquitously expressed in all immune cell types. The YBX1 gene is known for its association with the immune system and its pro-oncogenic properties in human hepatocellular carcinoma [109]. Upregulation of xbp1 expression has been observed in response to INH-lipopolysaccharide in zebrafish [110]. The fosab gene expression is upregulated under inflammatory conditions induced by Bisphenol AP in zebrafish [111].
A total of 17 transcription factors exhibited specific upregulation in the macrophage group (Figure 5, second section from the top), including five that are known to play key roles in monocyte maturation: klf6a, spi1a, mafba, irf8, and mafbb. The other 12 factors are junba (discussed above), csde1, nr4a1, fosl2, atf4b, atf4a, nfe2l2a, id1, jdp2b, klf2a, zfp36l1b, and nr1h3. The human CSDE1 gene is implicated in the immune recognition of nascent tumorigenic cells [112]. In mice, Nr4a1 deletion supports pro-angiogenic macrophage polarization and is essential for the maturation of inflammatory monocytes [113]. FOSL2 contributes to M2 macrophage polarization and enhances the tumorigenic capacity of glioblastoma cells [114]. The transcription factor ATF4 enhances M2 macrophage infiltration through the CCL2-dependent pathway [115]. Nfe2l2 plays a role in preventing the upregulation of pro-inflammatory cytokine expression and suppressing the pro-inflammatory phenotype of macrophages in mice [116]. Macrophages expressing the ID1 transcription factor have been shown to limit T cell infiltration in colorectal cancer [117]. Neomycin exposure has been shown to lead to the sustained activation of several genes in zebrafish macrophages, including jdp2b, which collectively mediate a switch to an anti-inflammatory state [118]. The expression of miR-2188-5p triggers the degradation of klf2a, leading to macrophage differentiation in Ctenopharyngodon idellus [119], suggesting that klf2a may support macrophage stemness. The anti-inflammatory properties of mouse macrophages are supported by the transcription factor Zfp36l1, whose expression significantly increases upon LPS treatment [120]. The NR1H3 gene is associated with a pro-inflammatory phenotype in human macrophages [121].
The next group of TFs identified by our meta-analysis is predominantly expressed in neutrophils: atf3, plek, id3, lrrfip1a, and raraa (Figure 5, second section from the bottom). An increase in the expression of PLEK along with genes associated with neutrophil activation and degranulation has been observed in response to osteonecrosis of the femoral head [122], suggesting a potential role for this TF in neutrophil activation in fish as well. Although the role of id3 in neutrophils remains unclear, it is known that exogenous supplementation with Id2 protein in mice can effectively protect them from colitis and inhibit NF-kB activation in neutrophils [123]. In human cells, LRRFIP1 has been shown to inhibit the NLRP1 inflammasome and regulate hematopoiesis [124]. The elevated expression of raraa in zebrafish neutrophils has also been confirmed in [125], and the expression of the Rarα gene is associated with neutrophil differentiation in mice [126].
The lower section of Figure 5 highlights additional marker TFs expressed in various cell types, including si:ch211-9d9.1, which is expressed in neutrophil-like cells; bhlhe40, found in B cells; foxp4 and si:ch211-161c3.6, both expressed in HSCs; and mych, which is present in both HSCs and B cells. It has been shown that Bhlhe40 plays a key role in activated mouse B cells by restraining the germinal center reaction and preventing lymphomagenesis [127]. The gene Foxp4 has been associated with regulating robust recall responses in mouse T lymphocytes [128]. The expression of mych was previously identified in zebrafish neutrophils, where it plays a role in inflammation and the formation of mycobacterial granulomas [129].
Thus, among identified genes and TFs, most have been functionally validated in various immune roles in mammals, which likely accounts for their detection as immune markers in zebrafish. Additionally, some genes already known to be specific to fish immune cells from the study [8] and TFs from four other studies [105,106,107,108] have been confirmed. However, these scattered data do not yet allow us to piece together a coherent picture of myelopoiesis in fish, indicating a need for further investigation. The newly identified candidate transcription factors for macrophages and neutrophils in fish suggest complex regulation of these cell types and their diverse functional roles. Therefore, in the subsequent section of the results, we focus on the reconstruction of gene networks consisting of highly expressed marker genes of neutrophils and macrophages, as well as transcription factors identified in our study. Using available data on protein–protein and co-expression interactions, we aim to visualize their molecular–genetic systems and characterize their functional roles.

3.3. Gene Networks and Functional Annotations of Key Macrophage and Neutrophil Marker Genes Identified through Integration of Zebrafish Kidney scRNA-Seq Datasets

During our meta-analysis, we encountered two main issues:
1.
Initial low coverage: The limited coverage of our data analysis was due to the inherent limitations of scRNA-seq technology (see Section 3.1). This limitation was partially addressed through a filtering procedure that allowed us to exclude cells with a minimal number of reads.
2.
Lack of precise information: There was insufficient detailed information on the immune properties and functions of individual genes in the identified immune cell type markers (see Section 3.2). This gap is often addressed by utilizing data from specific studies that examine the function of particular orthologous genes in humans or mouse models.
Given these challenges, we believe that the most effective strategy for bioinformatics analysis and interpretation of such data is to reconstruct gene networks of two major cell types identified during our meta-analysis (macrophages and neutrophils). This approach can provide a valuable blueprint for understanding the genetic systems involved in the peculiarities and regulation of macrophages and neutrophils. Combining data on protein–protein and co-expression interactions for macrophage-specific and neutrophil-specific genes (see Supplementary Table S4) from Group 1 (cell-type-specific highly expressed genes) and Group 2 (highly expressed genes), we reconstructed gene networks for highly expressed macrophage genes (Figure 6) and neutrophil genes (Figure 7). Genes without any protein–protein or co-expression associations were excluded from these networks.

3.3.1. Macrophage-Specific Gene Network

The resulting gene network of highly expressed macrophage genes (Figure 6) comprises 113 protein-coding genes connected through 764 interactions, of which 95 are protein–protein interactions, and 669 are co-expression associations. The Cytoscape session file (.cys) with this network is provided in Supplementary File S10. The network includes 13 transcription factors (highlighted at the top of the figure as hexagons). The transcription factors spi1a, cebpb, and mafba are already known for their roles in myelopoiesis [105,106,107,108]. We predicted the involvement of ten new potential regulators of macrophage gene expression: nr4a1, atf4a, fosl2, junba, fosab, xbp1, nfe2l2a, ybx1, nme2b.1, and id1. Collectively, the transcription factors are associated with 50 potential target genes through five protein–protein and 71 co-expression interactions.
It is worth noting that the newly predicted factors fosab and fosl2 were found to have protein–protein interactions with the already known important regulator of fish myelopoiesis, mafba, which, in turn, interacts with the candidate factors atf4a and junba. Additionally, the transcription factor cebpb, noted to be involved in myleopoiesis, has a protein–protein interaction with atf4a. Taken together, such interactomic data may hint at the complex nature of the regulation of maturation and polarization of zebrafish macrophages, which are able to restructure their metabolic properties depending on their spatial context [16] and the presence of foreign agents, wounds, and other factors [15].
Among the genes presented in the macrophage-specific network, there are three small clusters highly enriched with internal protein–protein interactions:
1.
Twelve genes are related to immune functions. Six cathepsins: ctsa, ctsba, ctsc, ctsh, ctsk, and ctsz; four genes encoding major histocompatibility complex class II proteins (si:busm1-266f07.2 (mhc2a), mhc2dab, cd74a, and cd74b); peptidoglycan recognition protein 5 (pglyrp5) and macrophage migration inhibitory factor (mif);
2.
Eight genes encode integral membrane proteins. Three genes of transporters of proteins into the endoplasmic reticulum (sec61a1, sec61b, and sec61g); three genes of signal sequence receptors (ssr2, ssr3, and ssr4); transmembrane protein 258 (tmem258); a subunit of the oligosaccharyl transferase complex dad1;
3.
Six genes of proton-transporting V-type ATPase complex: atp6ap2, atp6v0e1, atp6v0ca, atp6v1g1, atp6v1e1b, and rnaseka.
The network includes 15 out of the 20 key macrophage markers identified in the study [8], which are marked by ellipses outlined in black. Additionally, the network encompasses 34 potential macrophage markers, highlighted by ellipses outlined in blue. The most highly connected nodes in the network with more than five protein–protein interactions are cd74a, ctsc, ctsh, ctsk, mhc2dab, sec61a1, sec61b, sec61g, ssr2, and ssr3, while the most co-expressed genes with more than 35 associations include atp6v0ca, c1qc, cfp, ctsa, ctsba, ctss2.2, fgl2a, grna, mafba, ms4a17a.10, psap, and si:dkey-5n18.1. Genes with highest expression in network (log2(ACPC) > 5 ) are nme2b.1, s100a10b, grna.1, fthl27, cd74a, and hspa8. Therefore, the significant enrichment of protein–protein interactions is primarily observed due to associations between protein-coding genes rather than transcription factor–gene interactions, highlighting the need for further research into the regulatory mechanisms of myelopoiesis and the specific roles of individual transcription factors within this network.
The identified co-expression patterns may provide insights into potential targets for specific transcription factors, which are detailed in Supplementary Table S7.1. Among the transcription factors, mafba exhibits the highest number of co-expression associations within the network, with 48 interactions, followed by id1 with 17 connections. Additionally, the transcription factor xbp1 engages in four protein–protein interactions with cd74a, sec61a1, sec61b, and sec61g.

3.3.2. Neutrophil-Specific Gene Network

The reconstructed gene network of highly expressed neutrophil genes (Figure 7) encompasses 124 protein-coding genes with 2436 interactions, 67 of which are protein–protein interactions, and 2369 are co-expression interactions. The Cytoscape session file (.cys) with this network is provided in Supplementary File S11. The network includes 13 transcription factors, five of which are known: cebpa, cebpb, hif1ab, spi1b, and stat1a [105,106,107,108], and nine newly predicted factors: atf3, atf4b, id3, lrrfip1a, nme2b.1, plek, xbp1, and ybx1. The factors atf4b and cebpb are linked by a protein–protein interaction, while nme2b.1 and ybx1 are co-expressed in neutrophils. The identified transcription factors are connected to 79 potential targets through four protein–protein and 201 co-expression interactions. The network includes 11 of the 20 key neutrophil markers from [8], highlighted by black-bordered ellipses. Additionally, 18 potential neutrophil markers, indicated by blue-bordered ellipses, are present in the network. The genes with more than five protein–protein interactions in this network are actb1, actb2, actr2a, arpc1b, arpc3, arpc4l, and arpc5b, while the genes with more than 85 co-expression associations are scpp8, alox5ap, si:ch73-343l4.8, ppdpfa, si:ch211-284o19.8, illr4, adam8a, mapre1a, ctss1, vamp8, hacd4, lta4h, and mpx. The most highly expressed genes in the network (log2(ACPC) > 5 ) are lyz, lect2.1, actb2, tmsb4x, pfn1, npsn, srgn, icn, actb1, and cpa5.
Neutrophils demonstrate more pronounced overall co-expression compared to macrophages, which might indirectly suggest their more deterministic transcriptional regulation compared to macrophages. On the other hand, the presence of only one protein–protein interaction between their TFs (cebpb-atf4b) may speak in favor of their weaker co-regulation, which could mediate their specific pattern of target gene regulation.
In this neutrophil-specific network, there are 21 genes with a large number of internal protein–protein interactions:
1.
Two genes encoding actin (actb1, actb2) and five genes of actin-related protein 2/3 complex subunit: actr2a, arpc1b, arpc3, arpc4l, and arpc5b;
2.
Six genes related to actin cytoskeleton organization: plastin-2 lcp1, cofilin 1-like cfl1l, profilin pfn1, thymosin beta tmsb4x, myosin heavy chain 9a myh9a, and tropomyosin 1 alpha tpm1;
3.
Four genes involved in glycolytic process: enolase 1a eno1a, enolase 3 eno3, glyceraldehyde-3-phosphate dehydrogenase 2 gapdhs, and glucose-6-phosphate isomerase gpia;
4.
Glutathione reductase gsr, myosin light chain 12 myl12.1, transketolase b tktb, and zgc:153867.
The interactions of neutrophil transcription factors and their potential targets are detailed in Supplementary Table S7.2. The transcription factors plek (47 connections), id3 (41 connections), atf3 (32 connections), lrrfip1a (28 connections), spi1b (19 connections), and xbp1 (16 connections) exhibit the highest number of co-expression associations with genes in the network.

3.3.3. Functional Annotation of Highly Expressed Genes in Zebrafish Macrophages and Neutrophils

For macrophage and neutrophil marker genes from Group 1 and Group 2, we identified 15 and 25 enriched KEGG pathway terms, respectively. Seven of these pathways are shared in both cell types, though they often involve different genes. A comprehensive list of all enriched pathways and associated genes can be found in Supplementary Table S8, with key pathways highlighted in Figure 8 below.
Macrophage cells show significant enrichment in pathways such as “Protein processing in the endoplasmic reticulum”, “Protein export”, and “RNA degradation”, which are essential for their functional roles [130,131]. These cells also display specific enrichment in metabolic pathway components, indicating high metabolic activity and suggesting activation or polarization processes similar to those observed in humans. For instance, the metabolic pathways of glycine, serine, threonine, and pyruvate, which are linked to macrophage polarization in humans [132,133], are also important in zebrafish, where serine metabolism plays a crucial role in responding to bacterial infections [134]. Another enriched pathway, oxidative phosphorylation (Figure 8A), is implicated in macrophage polarization in humans [135], marking anti-inflammatory macrophages in contrast to pro-inflammatory macrophages, which rely on glycolysis. Additionally, during zebrafish tissue regeneration, macrophages initially exhibit glucocorticoid activation, followed by IL-10 signaling, and eventually, the induction of oxidative phosphorylation via IL-4/Polyamine signaling [118]. Our analysis further reveals that regulators of ferroptosis are specifically associated with zebrafish macrophages. While data on this process in fish immune cells are limited [136], in mammals, macrophages are crucial in regulating ferroptosis within tissues [137]. The enrichment of the term “Herpes simplex virus 1 infection” in macrophages is particularly noteworthy, as zebrafish are increasingly recognized as a valuable model for studying herpesvirus infections in mammals [138]. In humans, macrophages play a central role in the immune response against this virus [139], utilizing the enzyme acid ceramidase, which is localized in vesicles [140]. Recent studies indicate that acid ceramidase is also present in zebrafish and localizes to lysosomes [141], further supporting zebrafish as an effective model for investigating the antiviral activities of the vertebrate innate immune system.
We also identified the enrichment of common pathways in zebrafish macrophages and neutrophils (Figure 8B). In particular, the NOD-like receptor signaling pathway, which plays a crucial role in the antimicrobial activity of these cells across a wide range of multicellular organisms [142], was prominently enriched. Additionally, components of the lysosome, phagosome, and apoptosis pathways, all essential for the function of these cell types [143,144,145], were also significantly represented. Both macrophages and neutrophils exhibited a marked increase in the expression of components associated with heightened metabolic and energy activity, including carbon metabolism, metabolic pathways, and glycolysis/gluconeogenesis. For example, mammalian neutrophils store glycogen, which is critical for their survival and function, with the glycogen cycle being fundamental to the metabolism of these cells [146]. Collectively, these findings confirm that the innate immune cells analyzed in this study exhibit consistent transcriptomic signatures.
We observed specific enrichment of terms that reflect the unique biology of neutrophils (Figure 8C). Notably, genes involved in vascular smooth muscle contraction were found to be specific to this cell type. Both neutrophils [147] and macrophages are known to participate in the regulation of vascular smooth muscle cells in humans, with several molecular, genetic, and biochemical mechanisms underlying their adhesive interactions with these cells [148]. Neutrophils influence smooth muscle cells through the release of extracellular traps and are considered crucial players in the development of pathological conditions such as hypertension [147]. Additionally, monocytes, macrophages, and vascular smooth muscle cells are key contributors to the progression of atherosclerosis [149]. Our findings, along with existing data on the role of innate immune cells in the regulation of vascular smooth muscle cell proliferation [147] and cardiomyocyte proliferation during zebrafish cardiac regeneration [150], suggest that the interaction of immune cells with vascular smooth muscle cells in fish may occur similarly to that in humans, making zebrafish a promising model for studying a wide range of cardiovascular pathologies. Furthermore, the enrichment of neutrophil marker genes with components of the GnRH (gonadotropin-releasing hormone) signaling pathway in fish reveals parallels with the functional characteristics of these cells in humans. In humans, GnRH influences the formation of neutrophil extracellular traps and plays a role in processes such as wound healing in diabetes [151]. Although the effect of GnRH on neutrophils in fish has yet to be studied, these parallels suggest a similar functional presence. Additionally, among the pathways enriched in fish neutrophils, the three genes odc1, pgd, and gsr are responsible for glutathione metabolism. In mouse neutrophils, reduced glutathione is essential for maintaining their oxidative burst capabilities and antimicrobial activity [152]. The C-type lectin receptor signaling pathway is uniquely enriched in zebrafish neutrophils. In mice, activated dermal neutrophils, especially in cases of epidermolysis bullosa acquisita, show upregulation of multiple C-type lectin receptors [153]. Lastly, genes associated with adherens junctions and tight junctions were also significantly enriched in fish neutrophils. In mammalian neutrophils, these junctions play a critical role during inflammation and in neutrophil–epithelial interactions [154].

3.3.4. Identified Molecular Profiles of Pro- and Anti-Inflammatory Immune Cells

In our dataset, the majority of macrophages and neutrophils are characterized by the expression of typical marker genes associated with mature immune cells of these types. A crucial functional aspect of these cells’ homeostasis is the lability of their transcriptional programs and their capacity to polarize into subtypes with contrasting properties. A simplified classification divides macrophages into pro-inflammatory M1 macrophages, which stimulate pathogen defense, and anti-inflammatory M2 macrophages, which facilitate wound healing and regeneration. However, M1-like human macrophages are involved in both tumor-promoting and anti-tumor processes, and recent studies have revealed extensive immunophenotypes and plasticity in both M1-like and M2-like cells [155]. A similar classification exists for neutrophils, dividing them into pro- and anti-inflammatory subtypes, with the main cytokines regulating their polarization being well-known [156]. In zebrafish, key genes expressed in M1 macrophages include tnfa, tnfb, il1b, and il6, while M2 macrophages are known to express tgfb1, ccr2, and cxcr4b at high levels [157]. Similarly, pro-inflammatory neutrophils in fish express tnfa, tnfb, il1b, and il6 [158], whereas anti-inflammatory neutrophils are characterized by the expression of sgk1 and hif1a [159].
Therefore, we sought to identify macrophage and neutrophil cells within our dataset that were strongly polarized toward either the pro-inflammatory or anti-inflammatory phenotype, aiming to reveal the typical molecular signatures of these distinct subpopulations. The lists of genes that differentiate these subpopulations from the broader immune cell population are provided in Supplementary Table S9. Among the entire macrophage community (15,365 cells), we identified 39 cells with an extreme pro-inflammatory phenotype, characterized by the high-level expression of all four markers (il1b, il6, tnfa, and tnfb) (marker genes for these cells are listed in Supplementary Table S9.1). However, we did not detect any macrophages with the M2 phenotype described previously, which would be indicated by the simultaneous non-zero expression of tgfb1a, tgfb1b, cxcr4b, and ccr2. Among neutrophils (26,556 cells), we identified a small subset of pro-inflammatory cells (27 in total) expressing il1b, il6, tnfa, and tnfb, with their marker genes provided in Supplementary Table S9.2. We also found a group of anti-inflammatory neutrophils (28 cells in total) that expressed sgk1 and hif1a at the highest levels (Supplementary Table S9.3). These findings suggest that only a very small percentage of immune cells in our sample can be definitively categorized into one polarization phenotype or the other. The vast majority of mature immune cells in the fish kidney retain transcriptomic lability, supporting the notion that the binary classification of the innate immune system (IIS) into pro-inflammatory and anti-inflammatory cells may be outdated. This concept is further supported by recent studies highlighting the overall immunophenotypic plasticity of these cells [155,160]. Although this classification of immune cells has its limitations, we further examined which of the macrophage-specific (Section 3.3.1) and neutrophil-specific (Section 3.3.2) markers identified in our study are associated with pro- or anti-inflammatory processes.
Macrophage-specific genes: Ten identified markers were significantly upregulated in M1 macrophages, including ctsz, vmp1, tspan36, sftpbb, ccl34a.4, lgmn, fgl2a, mfap4.1, sdc4, and lgals3bp.1. Additionally, three potential macrophage markers, atf4a, pglyrp5, and fthl27, were upregulated across all three polarized immune populations (M1 macrophages, pro-inflammatory neutrophils, and anti-inflammatory neutrophils), suggesting their critical role in the zebrafish innate immune response. The well-known macrophage-specific gene lygl1 and transcription factor junba were not upregulated during M1 macrophage polarization but were instead expressed in both pro-inflammatory and anti-inflammatory neutrophils. Furthermore, nr4a1, ckbb, and ctsd were specifically upregulated in pro-inflammatory neutrophils, while the transcription factor nfe2l2a was upregulated in anti-inflammatory neutrophils.
Neutrophil-specific genes: Six potential neutrophil marker genes were significantly upregulated in all polarized subsets, including il1b, rgs2, atf3, hif1ab, stat1a, and ncf1. This set intriguingly includes both the pro-inflammatory marker il1b and the anti-inflammatory hif1ab, which suggests that transcriptional regulation of immune cell polarization in zebrafish may not be strictly discrete but instead may involve a gradient expression of key interleukins and mediators. Such insights are crucial for constructing more accurate cell-based models in the future, particularly when relying on scRNA-seq data [161]. The transcription factors atf4b and lrrfip1a, as well as the well-known marker mpx, were upregulated in both polarized subsets of neutrophils. Additionally, the gene srgn and transcription factor id3 were upregulated in M1 macrophages. The transcription factor plek was upregulated in both anti-inflammatory neutrophils and M1 macrophages, while the genes amd1 and hsd3b7 were specifically upregulated in anti-inflammatory neutrophils.
In summary, a small proportion of the markers identified in our gene networks are specifically associated with cell polarization. Additionally, several marker genes show upregulated expression during the polarization of another cell type, suggesting their functional significance in both neutrophils and macrophages. It is remarkable that the majority of macrophage (39 out of 62) and neutrophil (25 out of 43) markers did not exhibit differential expression during polarization, indicating their universal roles in mature immune cells. The high expression of these genes appears to be less dependent on the polarization state of the cell. Therefore, we propose that future studies should shift their focus from the individual properties of specific immune cell populations and subtypes to the complex interactions and regulatory mechanisms governing entire communities and subtypes of fish immune cells across various processes. As highlighted in a recent review, there is a notable gap in research on the cooperation and interactions among zebrafish neutrophils and macrophages [6]. Our meta-analysis underscores the dominant role these immune cells play in the zebrafish IIS and reveals that their classification is far more complex than the simple dichotomy of pro- and anti-inflammatory cells.

4. Conclusions

Our study demonstrates the effectiveness of meta-analysis as an approach to identifying key properties of large heterogeneous immune cell populations. The flexibility and iterativeness of this strategy in filtering cell types, identifying cell types and marker genes, and the ability to combine such data with interactomic analysis make it a powerful tool for predicting molecular–genetic systems of certain cell populations. In total, we identified more than two hundred highly expressed genes for macrophages and neutrophils of zebrafish kidneys. For more than half of the newly identified marker genes of them, we found literature evidence of the function of its direct or indirect ortholog in the same immune types in humans or mice. These data are still fragmentary and allow only a rough statement about the general similarity of mammalian and fish immune cells, but they can become the basis for further detailed immunogenetics and evolutionary studies. Also, macrophages and neutrophils account for about half of the cells in the immune dataset, highlighting their prevalence and importance for fish immunity. Our analysis also highlights the high complexity and heterogeneity within these cell types and the weakness of their binary classification into pro- and anti-inflammatory phenotypes.
On the other hand, meta-analysis is not an optimal tool for studying rare cell types. Small cell populations may be excluded or mistakenly merged with other cell subtypes during such integrative analyses, making their identification and annotation difficult. Our analysis did not detect zebrafish eosinophils and NK cells. To study the functional roles and features of such minor cell types, separate special experiments are needed, including additional steps such as cell sorting, immunostaining, and probably the usage of specific transgenic fish lines.
One of the main limitations of scRNAseq data is low read coverage. Our analysis involved about 250,000 immune cells, less than half of which passed the medium-stringency filtering threshold. Therefore, the integrative analysis of large scRNAseq datasets cannot yet guarantee accurate prediction of more than several hundred marker genes per cell type. Our results reflect only the highest expressed genes of these cell types. In turn, markers with low counts per cell cannot be considered unambiguous, and they require additional validation. We could try to bypass this limitation by merging individual cells, considering them metacells. However, this approach is quite slippery since it introduces significant Bayes associated with the specific characteristics of the analyzed dataset and does not guarantee that low-expressed marker genes identified by this method can be identified in another experiment. However, single-cell technologies continue to rapidly develop, and this issue could be overcome in the next few years.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/biology13100773/s1, Figure S1: Expression of the key known markers in harmonized UMAP space; Figure S2: Initial analysis of selected immune cells; Figure S3. Expression of the top 10 markers for revealed major immune cell types; Table S4: Detected marker genes for revealed populations of immune cells; Table S5: Marker genes of immune cells and their expression in individual Seurat clusters; Table S6: The detailed information about selected experiments used in meta-analysis; Table S7. Potential targets of transcription factors of zebrafish macrophages and neutrophils; Table S8: Significantly enriched KEGG pathways in neutrophil- and macrophage-specific genes; Table S9: Genes upregulated in subpopulations of M1 macrophages, pro-inflammatory neutrophils, and anti-inflammatory neutrophils; File S10: Reconstructed gene network of highly expressed genes of zebrafish kidney macrophages in Cytoscape format; File S11: Reconstructed gene network of highly expressed genes of zebrafish neutrophils in Cytoscape format; File S12: The joined expression matrix consisting of 94,938 immune cells × 3737 immune-related genes in .tsv.zip format.

Author Contributions

Conceptualization, A.V.D., A.V.B. and U.S.Z.; methodology, A.V.B., A.V.D. and U.S.Z.; formal analysis, A.V.B.; investigation, A.V.B., A.V.D., U.S.Z. and L.G.N.; writing—original draft preparation, A.V.B., U.S.Z., A.V.D. and L.G.N.; writing—review and editing, A.V.B., U.S.Z. and A.V.D.; visualization, A.V.B., U.S.Z., L.G.N. and A.V.D.; preparation of supplementary materials, A.V.B.; supervision, A.V.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian Science Foundation, grant number 23-24-00536 (accessed on 28 July 2024). The APC was funded by Russian Science Foundation, grant number 23-24-00536.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this study are available from a publicly accessible repository Gene Expression Omnibus database. The identifiers of particular experiments and corresponding articles are available in Table 1. The joined expression matrix, which was used in our final analysis and consisted of 94,938 immune cells × 3737 immune-related genes, is provided in Supplementary File S12. Other additional materials and results of our calculations can be provided to those interested upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
IISInnate immune system
HSPCsHematopoietic stem and progenitor cells
scRNAseqSingle-cell RNA sequencing
UMAPUniform manifold approximation and projection
ACPCAverage gene counts per cell
TFTranscription factor

References

  1. Riera Romo, M.; Pérez-Martínez, D.; Castillo Ferrer, C. Innate immunity in vertebrates: An overview. Immunology 2016, 148, 125–139. [Google Scholar] [CrossRef] [PubMed]
  2. De Kleer, I.; Willems, F.; Lambrecht, B.; Goriely, S. Ontogeny of myeloid cells. Front. Immunol. 2014, 5, 423. [Google Scholar] [CrossRef]
  3. Campos-Sánchez, J.C.; Esteban, M.Á. Review of inflammation in fish and value of the zebrafish model. J. Fish Dis. 2021, 44, 123–139. [Google Scholar] [CrossRef]
  4. Bobrovskikh, A.V.; Zubairova, U.S.; Doroshkov, A.V. Fishing innate immune system properties through the transcriptomic single-cell data of Teleostei. Biology 2023, 12, 1516. [Google Scholar] [CrossRef]
  5. Bjørgen, H.; Koppang, E.O. Anatomy of teleost fish immune structures and organs. In Principles of Fish Immunology: From Cells and Molecules to Host Protection; Springer: Cham, Switzerland, 2022; pp. 1–30. [Google Scholar] [CrossRef]
  6. Speirs, Z.C.; Loynes, C.A.; Mathiessen, H.; Elks, P.M.; Renshaw, S.A.; von Gersdorff Jørgensen, L. What can we learn about fish neutrophil and macrophage response to immune challenge from studies in zebrafish. Fish Shellfish Immunol. 2024, 148, 109490. [Google Scholar] [CrossRef]
  7. Mastrogiovanni, M.; Martínez-Navarro, F.J.; Bowman, T.V.; Cayuela, M.L. Inflammation in development and aging: Insights from the zebrafish model. Int. J. Mol. Sci. 2024, 25, 2145. [Google Scholar] [CrossRef]
  8. Hu, C.; Zhang, N.; Hong, Y.; Tie, R.; Fan, D.; Lin, A.; Chen, Y.; Xiang, L.x.; Shao, J.z. Single-cell RNA sequencing unveils the hidden powers of zebrafish kidney for generating both hematopoiesis and adaptive antiviral immunity. Elife 2024, 13, RP92424. [Google Scholar] [CrossRef]
  9. Kirchberger, S.; Shoeb, M.R.; Lazic, D.; Wenninger-Weinzierl, A.; Fischer, K.; Shaw, L.E.; Nogueira, F.; Rifatbegovic, F.; Bozsaky, E.; Ladenstein, R.; et al. Comparative transcriptomics coupled to developmental grading via transgenic zebrafish reporter strains identifies conserved features in neutrophil maturation. Nat. Commun. 2024, 15, 1792. [Google Scholar] [CrossRef]
  10. Kiss, M.; Van Gassen, S.; Movahedi, K.; Saeys, Y.; Laoui, D. Myeloid cell heterogeneity in cancer: Not a single cell alike. Cellular Immunol. 2018, 330, 188–201. [Google Scholar] [CrossRef]
  11. Ma, R.Y.; Black, A.; Qian, B.Z. Macrophage diversity in cancer revisited in the era of single-cell omics. Trends Immunol. 2022, 43, 546–563. [Google Scholar] [CrossRef]
  12. McLaren, A.S.; Fetit, R.; Wood, C.S.; Falconer, J.; Steele, C.W. Single cell sequencing of neutrophils demonstrates phenotypic heterogeneity and functional plasticity in health, disease, and cancer. Chin. Clin. Oncol. 2023, 12, 18. [Google Scholar] [CrossRef] [PubMed]
  13. Hernández, P.P.; Strzelecka, P.M.; Athanasiadis, E.I.; Hall, D.; Robalo, A.F.; Collins, C.M.; Boudinot, P.; Levraud, J.P.; Cvejic, A. Single-cell transcriptional analysis reveals ILC-like cells in zebrafish. Sci. Immunol. 2018, 3, eaau5265. [Google Scholar] [CrossRef] [PubMed]
  14. Rougeot, J.; Torraca, V.; Zakrzewska, A.; Kanwal, Z.; Jansen, H.J.; Sommer, F.; Spaink, H.P.; Meijer, A.H. RNAseq profiling of leukocyte populations in zebrafish larvae reveals a cxcl11 chemokine gene as a marker of macrophage polarization during mycobacterial infection. Front. Immunol. 2019, 10, 832. [Google Scholar] [CrossRef]
  15. Hu, C.B.; Wang, J.; Hong, Y.; Li, H.; Fan, D.D.; Lin, A.F.; Xiang, L.X.; Shao, J.Z. Single-cell transcriptome profiling reveals diverse immune cell populations and their responses to viral infection in the spleen of zebrafish. FASEB J. 2023, 37, e22951. [Google Scholar] [CrossRef] [PubMed]
  16. Zhou, Q.; Zhao, C.; Yang, Z.; Qu, R.; Li, Y.; Fan, Y.; Tang, J.; Xie, T.; Wen, Z. Cross-organ single-cell transcriptome profiling reveals macrophage and dendritic cell heterogeneity in zebrafish. Cell Rep. 2023, 42, 112793. [Google Scholar] [CrossRef]
  17. Zernecke, A.; Winkels, H.; Cochain, C.; Williams, J.W.; Wolf, D.; Soehnlein, O.; Robbins, C.S.; Monaco, C.; Park, I.; McNamara, C.A.; et al. Meta-analysis of leukocyte diversity in atherosclerotic mouse aortas. Circ. Res. 2020, 127, 402–426. [Google Scholar] [CrossRef]
  18. Muus, C.; Luecken, M.D.; Eraslan, G.; Sikkema, L.; Waghray, A.; Heimberg, G.; Kobayashi, Y.; Vaishnav, E.D.; Subramanian, A.; Smillie, C.; et al. Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics. Nat. Med. 2021, 27, 546–559. [Google Scholar] [CrossRef]
  19. Liu, F.; Yuwono, C.; Tay, A.C.Y.; Wehrhahn, M.C.; Riordan, S.M.; Zhang, L. Analysis of global Aeromonas veronii genomes provides novel information on source of infection and virulence in human gastrointestinal diseases. BMC Genom. 2022, 23, 166. [Google Scholar] [CrossRef]
  20. Botos, M.A.; Arora, P.; Chouvardas, P.; Mercader, N. Transcriptomic data meta-analysis reveals common and injury model specific gene expression changes in the regenerating zebrafish heart. Sci. Rep. 2023, 13, 5418. [Google Scholar] [CrossRef]
  21. Bobrovskikh, A.V.; Zubairova, U.S.; Bondar, E.I.; Lavrekha, V.V.; Doroshkov, A.V. Transcriptomic data meta-analysis sheds light on high light response in Arabidopsis thaliana L. Int. J. Mol. Sci. 2022, 23, 4455. [Google Scholar] [CrossRef]
  22. Liu, Y.; Li, G. Empowering biologists to decode omics data: The Genekitr R package and web server. BMC Bioinform. 2023, 24, 214. [Google Scholar] [CrossRef] [PubMed]
  23. Satija, R.; Farrell, J.A.; Gennert, D.; Schier, A.F.; Regev, A. Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 2015, 33, 495–502. [Google Scholar] [CrossRef] [PubMed]
  24. Korsunsky, I.; Millard, N.; Fan, J.; Slowikowski, K.; Zhang, F.; Wei, K.; Baglaenko, Y.; Brenner, M.; Loh, P.r.; Raychaudhuri, S. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat. Methods 2019, 16, 1289–1296. [Google Scholar] [CrossRef] [PubMed]
  25. Shen, W.K.; Chen, S.Y.; Gan, Z.Q.; Zhang, Y.Z.; Yue, T.; Chen, M.M.; Xue, Y.; Hu, H.; Guo, A.Y. AnimalTFDB 4.0: A comprehensive animal transcription factor database updated with variation and expression annotations. Nucleic Acids Res. 2023, 51, D39–D45. [Google Scholar] [CrossRef] [PubMed]
  26. Morabito, S.; Reese, F.; Rahimzadeh, N.; Miyoshi, E.; Swarup, V. hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data. Cell Rep. Methods 2023, 3, 100498. [Google Scholar] [CrossRef] [PubMed]
  27. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  28. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for functional genomics data sets—update. Nucleic Acids Res. 2012, 41, D991–D995. [Google Scholar] [CrossRef] [PubMed]
  29. Tang, Q.; Iyer, S.; Lobbardi, R.; Moore, J.C.; Chen, H.; Lareau, C.; Hebert, C.; Shaw, M.L.; Neftel, C.; Suva, M.L.; et al. Dissecting hematopoietic and renal cell heterogeneity in adult zebrafish at single-cell resolution using RNA sequencing. J. Exp. Med. 2017, 214, 2875–2887. [Google Scholar] [CrossRef]
  30. Baron, C.S.; Barve, A.; Muraro, M.J.; van der Linden, R.; Dharmadhikari, G.; Lyubimova, A.; de Koning, E.J.; van Oudenaarden, A. Cell type purification by single-cell transcriptome-trained sorting. Cell 2019, 179, 527–542. [Google Scholar] [CrossRef]
  31. Jiang, M.; Xiao, Y.; Ma, L.; Wang, J.; Chen, H.; Gao, C.; Liao, Y.; Guo, Q.; Peng, J.; Han, X.; et al. Characterization of the zebrafish cell landscape at single-cell resolution. Front. Cell Dev. Biol. 2021, 9, 743421. [Google Scholar] [CrossRef] [PubMed]
  32. Avagyan, S.; Henninger, J.; Mannherz, W.; Mistry, M.; Yoon, J.; Yang, S.; Weber, M.; Moore, J.; Zon, L. Resistance to inflammation underlies enhanced fitness in clonal hematopoiesis. Science 2021, 374, 768–772. [Google Scholar] [CrossRef]
  33. Avagyan, S.; Weber, M.C.; Ma, S.; Prasad, M.; Mannherz, W.P.; Yang, S.; Buenrostro, J.D.; Zon, L.I. Single-cell ATAC-seq reveals GATA2-dependent priming defect in myeloid and a maturation bottleneck in lymphoid lineages. Blood Adv. 2021, 5, 2673–2686. [Google Scholar] [CrossRef] [PubMed]
  34. Amanda, S.; Tan, T.K.; Ong, J.Z.L.; Theardy, M.S.; Wong, R.W.J.; Huang, X.Z.; Ali, M.Z.; Li, Y.; Gong, Z.; Inagaki, H.; et al. IRF4 drives clonal evolution and lineage choice in a zebrafish model of T-cell lymphoma. Nat. Commun. 2022, 13, 2420. [Google Scholar] [CrossRef] [PubMed]
  35. Yan, C.; Yang, Q.; Zhang, S.; Millar, D.G.; Alpert, E.J.; Do, D.; Veloso, A.; Brunson, D.C.; Drapkin, B.J.; Stanzione, M.; et al. Single-cell imaging of T cell immunotherapy responses in vivo. J. Exp. Med. 2021, 218, e20210314. [Google Scholar] [CrossRef] [PubMed]
  36. Rubin, S.A.; Baron, C.S.; Pessoa Rodrigues, C.; Duran, M.; Corbin, A.F.; Yang, S.P.; Trapnell, C.; Zon, L.I. Single-cell analyses reveal early thymic progenitors and pre-B cells in zebrafish. J. Exp. Med. 2022, 219, e20220038. [Google Scholar] [CrossRef] [PubMed]
  37. Binder, V.; Li, W.; Faisal, M.; Oyman, K.; Calkins, D.L.; Shaffer, J.; Teets, E.M.; Sher, S.; Magnotte, A.; Belardo, A.; et al. Microenvironmental control of hematopoietic stem cell fate via CXCL8 and protein kinase C. Cell Rep. 2023, 42, 112528. [Google Scholar] [CrossRef]
  38. Li, G.; Sun, Y.; Kwok, I.; Yang, L.; Wen, W.; Huang, P.; Wu, M.; Li, J.; Huang, Z.; Liu, Z.; et al. Cebp1 and Cebpβ transcriptional axis controls eosinophilopoiesis in zebrafish. Nat. Commun. 2024, 15, 811. [Google Scholar] [CrossRef] [PubMed]
  39. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  40. Harrison, P.W.; Amode, M.R.; Austine-Orimoloye, O.; Azov, A.G.; Barba, M.; Barnes, I.; Becker, A.; Bennett, R.; Berry, A.; Bhai, J.; et al. Ensembl 2024. Nucleic Acids Res. 2024, 52, D891–D899. [Google Scholar] [CrossRef]
  41. Szklarczyk, D.; Kirsch, R.; Koutrouli, M.; Nastou, K.; Mehryary, F.; Hachilif, R.; Gable, A.L.; Fang, T.; Doncheva, N.T.; Pyysalo, S.; et al. The STRING database in 2023: Protein–protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023, 51, D638–D646. [Google Scholar] [CrossRef] [PubMed]
  42. Pozzoli, O.; Vella, P.; Iaffaldano, G.; Parente, V.; Devanna, P.; Lacovich, M.; Lamia, C.L.; Fascio, U.; Longoni, D.; Cotelli, F.; et al. Endothelial fate and angiogenic properties of human CD34+ progenitor cells in zebrafish. Arterioscler. Thromb. Vasc. Biol. 2011, 31, 1589–1597. [Google Scholar] [CrossRef] [PubMed]
  43. Szulc-Dąbrowska, L.; Bossowska-Nowicka, M.; Struzik, J.; Toka, F.N. Cathepsins in bacteria-macrophage interaction: Defenders or victims of circumstance? Front. Cell. Infect. Microbiol. 2020, 10, 601072. [Google Scholar] [CrossRef] [PubMed]
  44. McArthur, S.; Juban, G.; Gobbetti, T.; Desgeorges, T.; Theret, M.; Gondin, J.; Toller-Kawahisa, J.E.; Reutelingsperger, C.P.; Chazaud, B.; Perretti, M.; et al. Annexin A1 drives macrophage skewing to accelerate muscle regeneration through AMPK activation. J. Clin. Investig. 2020, 130, 1156–1167. [Google Scholar] [CrossRef] [PubMed]
  45. Zhu, Q.; Pan, Q.Z.; Zhong, A.L.; Hu, H.; Zhao, J.J.; Tang, Y.; Hu, W.M.; Li, M.; Weng, D.S.; Chen, M.Y.; et al. Annexin A3 upregulates the infiltrated neutrophil-lymphocyte ratio to remodel the immune microenvironment in hepatocellular carcinoma. Int. Immunopharmacol. 2020, 89, 107139. [Google Scholar] [CrossRef]
  46. Xu, F.; Guo, M.; Huang, W.; Feng, L.; Zhu, J.; Luo, K.; Gao, J.; Zheng, B.; Kong, L.; Pang, T.; et al. Annexin A5 regulates hepatic macrophage polarization via directly targeting PKM2 and ameliorates NASH. Redox Biol. 2020, 36, 101634. [Google Scholar] [CrossRef] [PubMed]
  47. Ke, K.; Wu, Z.; Lin, J.; Lin, L.; Huang, N.; Yang, W. Increased expression of CD74 in atherosclerosis associated with inflammatory responses of endothelial cells and macrophages. Biochem. Genet. 2024, 62, 294–310. [Google Scholar] [CrossRef]
  48. Guo, Q.; Lu, T.; Xu, H.; Luo, Q.; Liu, Z.; Jiang, S.; Pan, J.; Lin, S.; Lin, M.; Guo, F. Identification of immune-related genes contributing to head and neck squamous cell carcinoma development using weighted gene co-expression network analysis. Cancer Rep. 2023, 6, e1808. [Google Scholar] [CrossRef]
  49. Wang, Z.; Zhao, Q.; Nie, Y.; Yu, Y.; Misra, B.B.; Zabalawi, M.; Chou, J.W.; Key, C.C.C.; Molina, A.J.; Quinn, M.A.; et al. Solute carrier family 37 member 2 (SLC37A2) negatively regulates murine macrophage inflammation by controlling glycolysis. Iscience 2020, 23, 101125. [Google Scholar] [CrossRef]
  50. Wentzel, A.S.; Petit, J.; van Veen, W.G.; Fink, I.R.; Scheer, M.H.; Piazzon, M.C.; Forlenza, M.; Spaink, H.P.; Wiegertjes, G.F. Transcriptome sequencing supports a conservation of macrophage polarization in fish. Sci. Rep. 2020, 10, 13470. [Google Scholar] [CrossRef]
  51. Li, Q.; Cui, K.; Xu, D.; Wu, M.; Mai, K.; Ai, Q. Molecular identification of peptidoglycan recognition protein 5 and its functional characterization in innate immunity of large yellow croaker, Larimichthys crocea. Dev. Comp. Immunol. 2021, 124, 104130. [Google Scholar] [CrossRef]
  52. Beisaw, A.; Kuenne, C.; Guenther, S.; Dallmann, J.; Wu, C.C.; Bentsen, M.; Looso, M.; Stainier, D.Y. AP-1 contributes to chromatin accessibility to promote sarcomere disassembly and cardiomyocyte protrusion during zebrafish heart regeneration. Circ. Res. 2020, 126, 1760–1778. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, X.; Huang, Y.; Gao, P.; Lv, Y.; Jia, D.; Sun, K.; Han, Y.; Hu, H.; Tang, Z.; Ren, X.; et al. Knockout of mafba causes inner-ear developmental defects in zebrafish via the impairment of proliferation and differentiation of ionocyte progenitor cells. Biomedicines 2021, 9, 1699. [Google Scholar] [CrossRef]
  54. Shaw, D.K.; Saraswathy, V.M.; McAdow, A.R.; Zhou, L.; Park, D.; Mote, R.; Johnson, A.N.; Mokalled, M.H. Elevated phagocytic capacity directs innate spinal cord repair. bioRxiv 2024. [Google Scholar] [CrossRef]
  55. Zhu, J.; Wang, H.; Ma, T.; He, Y.; Shen, M.; Song, W.; Wang, J.J.; Shi, J.P.; Wu, M.Y.; Liu, C.; et al. Identification of immune-related genes as prognostic factors in bladder cancer. Sci. Rep. 2020, 10, 19695. [Google Scholar] [CrossRef] [PubMed]
  56. Mangogna, A.; Varghese, P.M.; Agostinis, C.; Alrokayan, S.H.; Khan, H.A.; Stover, C.M.; Belmonte, B.; Martorana, A.; Ricci, G.; Bulla, R.; et al. Prognostic value of complement properdin in cancer. Front. Immunol. 2021, 11, 614980. [Google Scholar] [CrossRef]
  57. Ji, P.; Gong, Y.; Jin, M.l.; Wu, H.l.; Guo, L.W.; Pei, Y.C.; Chai, W.J.; Jiang, Y.Z.; Liu, Y.; Ma, X.Y.; et al. In vivo multidimensional CRISPR screens identify Lgals2 as an immunotherapy target in triple-negative breast cancer. Sci. Adv. 2022, 8, eabl8247. [Google Scholar] [CrossRef]
  58. Zhao, L.; Wang, L.f.; Wang, Y.c.; Liu, A.; Xiao, Q.w.; Hu, M.C.; Sun, M.z.; Hao, H.y.; Gao, Q.; Zhao, X.; et al. Deficiency of P2RY11 causes narcolepsy and attenuates the recruitment of neutrophils and macrophages in the inflammatory response in zebrafish. Cell Biol. Toxicol. 2024, 40, 36. [Google Scholar] [CrossRef] [PubMed]
  59. Rovira, M.; Ferrero, G.; Miserocchi, M.; Montanari, A.; Wittamer, V. A single-cell transcriptomic atlas reveals resident dendritic-like cells in the zebrafish brain parenchyma. bioRxiv 2023. [Google Scholar] [CrossRef]
  60. Hu, J.; Zhang, Y.; Hu, L.; Chen, H.; Wu, H.; Chen, J.; Xie, J.; Xu, B.; Wei, Z. A reduction of Syndecan-4 in macrophages promotes atherosclerosis by aggravating the proinflammatory capacity of macrophages. J. Transl. Med. 2022, 20, 319. [Google Scholar] [CrossRef]
  61. Zack, S.R.; Venkatesan, M.; Nikolaienko, R.; Cook, B.; Melki, R.; Zima, A.V.; Campbell, E.M. Altered vacuole membrane protein 1 (VMP1) expression is associated with increased NLRP3 inflammasome activation and mitochondrial dysfunction. Inflamm. Res. 2024, 73, 563–580. [Google Scholar] [CrossRef]
  62. Taylor, R.S.; Ruiz Daniels, R.; Dobie, R.; Naseer, S.; Clark, T.C.; Henderson, N.C.; Boudinot, P.; Martin, S.A.; Macqueen, D.J. Single cell transcriptomics of Atlantic salmon (Salmo salar L.) liver reveals cellular heterogeneity and immunological responses to challenge by Aeromonas salmonicida. Front. Immunol. 2022, 13, 984799. [Google Scholar] [CrossRef] [PubMed]
  63. Guyot, N.; Labas, V.; Harichaux, G.; Chessé, M.; Poirier, J.C.; Nys, Y.; Réhault-Godbert, S. Proteomic analysis of egg white heparin-binding proteins: Towards the identification of natural antibacterial molecules. Sci. Rep. 2016, 6, 27974. [Google Scholar] [CrossRef]
  64. Foulkes, M.J.; Henry, K.M.; Rougeot, J.; Hooper-Greenhill, E.; Loynes, C.A.; Jeffrey, P.; Fleming, A.; Savage, C.O.; Meijer, A.H.; Jones, S.; et al. Expression and regulation of drug transporters in vertebrate neutrophils. Sci. Rep. 2017, 7, 4967. [Google Scholar] [CrossRef] [PubMed]
  65. McKelvey, M.C.; Abladey, A.A.; Small, D.M.; Doherty, D.F.; Williams, R.; Scott, A.; Spek, C.A.; Borensztajn, K.S.; Holsinger, L.; Booth, R.; et al. Cathepsin S contributes to lung inflammation in acute respiratory distress syndrome. Am. J. Respir. Crit. Care Med. 2022, 205, 769–782. [Google Scholar] [CrossRef]
  66. Smith, J.N.; Dawson, D.M.; Christo, K.F.; Jogasuria, A.P.; Cameron, M.J.; Antczak, M.I.; Ready, J.M.; Gerson, S.L.; Markowitz, S.D.; Desai, A.B. 15-PGDH inhibition activates the splenic niche to promote hematopoietic regeneration. JCI Insight 2021, 6, e143658. [Google Scholar] [CrossRef] [PubMed]
  67. Ji, H.; Han, Y.; Jie, D.; Li, Y.; Yang, H.; Sun, H.; You, C.; Xiao, A.; Liu, Y. Decoding the biology and clinical implication of neutrophils in intracranial aneurysm. Ann. Clin. Transl. Neurol. 2024, 11, 958–972. [Google Scholar] [CrossRef] [PubMed]
  68. Hu, Z.; Wang, J.; Fu, Y.; Li, Y.; Jiang, Q. HSD3B7 as a Prognostic-Related Biomarker Predicts Poor Prognostic in ccRCC; Research Square: Durham, NC, USA, 2022. [Google Scholar] [CrossRef]
  69. García-López, J.P.; Grimaldi, A.; Chen, Z.; Meneses, C.; Bravo-Tello, K.; Bresciani, E.; Banderas, A.; Burgess, S.M.; Hernández, P.P.; Feijoo, C.G. Ontogenetically distinct neutrophils differ in function and transcriptional profile in zebrafish. Nat. Commun. 2023, 14, 4942. [Google Scholar] [CrossRef]
  70. Li, W.; Faisal, M.; Shaffer, J.L.; Teets, E.M.; Blaser, B.W. Epigenetic regulation of protein kinase C expression in the hematopoietic niche. Blood 2022, 140, 8624. [Google Scholar] [CrossRef]
  71. Wang, X.; Xiong, H.; Liang, D.; Chen, Z.; Li, X.; Zhang, K. The role of SRGN in the survival and immune infiltrates of skin cutaneous melanoma (SKCM) and SKCM-metastasis patients. BMC Cancer 2020, 20, 378. [Google Scholar] [CrossRef]
  72. Morino, K.; Kunimura, K.; Sugiura, Y.; Izumi, Y.; Matsubara, K.; Akiyoshi, S.; Maeda, R.; Hirotani, K.; Sakata, D.; Mizuno, S.; et al. Cholesterol sulfate limits neutrophil recruitment and gut inflammation during mucosal injury. Front. Immunol. 2023, 14, 1131146. [Google Scholar] [CrossRef] [PubMed]
  73. Pérez, L.A.; Leyton, L.; Valdivia, A. Thy-1 (CD90), integrins and syndecan 4 are key regulators of skin wound healing. Front. Cell Dev. Biol. 2022, 10, 810474. [Google Scholar] [CrossRef]
  74. Su, T.; He, Y.; Huang, Y.; Ye, M.; Guo, Q.; Xiao, Y.; Cai, G.; Chen, L.; Li, C.; Zhou, H.; et al. Myeloid-derived grancalcin instigates obesity-induced insulin resistance and metabolic inflammation in male mice. Nat. Commun. 2024, 15, 97. [Google Scholar] [CrossRef]
  75. Catz, S.D.; McLeish, K.R. Therapeutic targeting of neutrophil exocytosis. J. Leukoc. Biol. 2020, 107, 393–408. [Google Scholar] [CrossRef]
  76. Conrad, C.; Yildiz, D.; Cleary, S.J.; Margraf, A.; Cook, L.; Schlomann, U.; Panaretou, B.; Bowser, J.L.; Karmouty-Quintana, H.; Li, J.; et al. ADAM8 signaling drives neutrophil migration and ARDS severity. JCI Insight 2022, 7, e149870. [Google Scholar] [CrossRef]
  77. Du, Y.; Ma, Z.; Zheng, J.; Huang, S.; Yang, X.; Song, Y.; Dong, D.; Shi, L.; Xu, D. ATF3 positively regulates antibacterial immunity by modulating macrophage killing and migration functions. Front. Immunol. 2022, 13, 839502. [Google Scholar] [CrossRef] [PubMed]
  78. Zhao, Z.; Liu, T.; Liang, Y.; Cui, W.; Li, D.; Zhang, G.; Deng, Z.; Chen, M.; Sha, K.; Xiao, W.; et al. N2-polarized neutrophils reduce inflammation in rosacea by regulating vascular factors and proliferation of CD4+ T cells. J. Investig. Dermatol. 2022, 142, 1835–1844. [Google Scholar] [CrossRef]
  79. Wei, K.H.; Lin, I.T.; Chowdhury, K.; Lim, K.L.; Liu, K.T.; Ko, T.M.; Chang, Y.M.; Yang, K.C.; Lai, S.L.B. Comparative single-cell profiling reveals distinct cardiac resident macrophages essential for zebrafish heart regeneration. Elife 2023, 12, e84679. [Google Scholar] [CrossRef]
  80. Wang, Y.; Gu, Y.; Alexander, J.S.; Lewis, D.F. Preeclampsia status controls interleukin-6 and soluble IL-6 receptor release from neutrophils and endothelial cells: Relevance to increased inflammatory responses. Pathophysiology 2021, 28, 202–211. [Google Scholar] [CrossRef] [PubMed]
  81. Yurdagul, A., Jr.; Kong, N.; Gerlach, B.D.; Wang, X.; Ampomah, P.; Kuriakose, G.; Tao, W.; Shi, J.; Tabas, I. ODC (ornithine decarboxylase)-dependent putrescine synthesis maintains MerTK (MER tyrosine-protein kinase) expression to drive resolution. Arterioscler. Thromb. Vasc. Biol. 2021, 41, e144–e159. [Google Scholar] [CrossRef]
  82. Huo, Y.; Hu, X.; Lü, J.; Luo, F.; Liang, J.; Lei, H.; Lv, A. Single-cell transcriptome, phagocytic activity and immunohistochemical analysis of crucian carp (Carassius auratus) in response to Rahnella aquatilis infection. Fish Shellfish Immunol. 2023, 140, 108970. [Google Scholar] [CrossRef]
  83. Tokuhisa, M.; Kadowaki, T.; Ogawa, K.; Yamaguchi, Y.; Kido, M.A.; Gao, W.; Umeda, M.; Tsukuba, T. Expression and localisation of Rab44 in immune-related cells change during cell differentiation and stimulation. Sci. Rep. 2020, 10, 10728. [Google Scholar] [CrossRef] [PubMed]
  84. Singh, S.K.; Aravamudhan, S.; Armant, O.; Krüger, M.; Grabher, C. Proteome dynamics in neutrophils of adult zebrafish upon chemically-induced inflammation. Fish Shellfish Immunol. 2014, 40, 217–224. [Google Scholar] [CrossRef] [PubMed]
  85. Fetit, R.; White, M.; Mills, M.L.; Cortes-Lavaud, X.; McLaren, A.; Falconer, J.; Gilroy, K.; Nixon, C.; Kirschner, K.; Jackstadt, R.; et al. Characterising neutrophil subtypes in cancer using human and murine single-cell RNA sequencing datasets. bioRxiv 2023. [Google Scholar] [CrossRef]
  86. Chang, C.; Li, H.; Zhang, R. Zebrafish facilitate non-alcoholic fatty liver disease research: Tools, models and applications. Liver Int. 2023, 43, 1385–1398. [Google Scholar] [CrossRef] [PubMed]
  87. Saini, R.; Singh, S. Inducible nitric oxide synthase: An asset to neutrophils. J. Leukoc. Biol. 2019, 105, 49–61. [Google Scholar] [CrossRef]
  88. Olesch, C.; Ringel, C.; Brüne, B.; Weigert, A. Beyond immune cell migration: The emerging role of the sphingosine-1-phosphate receptor S1PR4 as a modulator of innate immune cell activation. Mediat. Inflamm. 2017, 2017, 6059203. [Google Scholar] [CrossRef]
  89. Loynes, C.A.; Lee, J.A.; Robertson, A.L.; Steel, M.J.; Ellett, F.; Feng, Y.; Levy, B.D.; Whyte, M.K.; Renshaw, S.A. PGE2 production at sites of tissue injury promotes an anti-inflammatory neutrophil phenotype and determines the outcome of inflammation resolution in vivo. Sci. Adv. 2018, 4, eaar8320. [Google Scholar] [CrossRef] [PubMed]
  90. Trojanek, J.B.; Michałkiewicz, J.; Grzywa-Czuba, R.; Jańczyk, W.; Gackowska, L.; Kubiszewska, I.; Helmin-Basa, A.; Wierzbicka-Rucińska, A.; Szalecki, M.; Socha, P. Expression of matrix metalloproteinases and their tissue inhibitors in peripheral blood leukocytes and plasma of children with nonalcoholic fatty liver disease. Mediat. Inflamm. 2020, 2020, 8327945. [Google Scholar] [CrossRef]
  91. Roca, F.J.; Mulero, I.; Lopez-Munoz, A.; Sepulcre, M.P.; Renshaw, S.A.; Meseguer, J.; Mulero, V. Evolution of the inflammatory response in vertebrates: Fish TNF-α is a powerful activator of endothelial cells but hardly activates phagocytes. J. Immunol. 2008, 181, 5071–5081. [Google Scholar] [CrossRef]
  92. Liu, Y.; Zhang, J.; Wang, Z.; Zhang, X.; Dai, Z.; Wu, W.; Zhang, N.; Liu, Z.; Zhang, J.; Luo, P.; et al. Identify the prognostic and immune profile of VSIR in the tumor microenvironment: A pan-cancer analysis. Front. Cell Dev. Biol. 2022, 10, 821649. [Google Scholar] [CrossRef] [PubMed]
  93. Huang, C. Tissue Damage Signaling Is Essential for Protective Neutrophil Responses to Microbial Infection in Zebrafish (Danio rerio). Ph.D. Thesis, Weill Medical College of Cornell University, New York, NY, USA, 2018. [Google Scholar]
  94. Dirks, R.P.; Ordas, A.; Jong-Raadsen, S.; Brittijn, S.A.; Haks, M.C.; Henkel, C.V.; Oravcova, K.; Racz, P.I.; Tuinhof-Koelma, N.; Korzeniowska nee Wiweger, M.I.; et al. The human pathogen Mycobacterium tuberculosis and the fish pathogen Mycobacterium marinum trigger the same core set of late innate immune response genes in zebrafish larvae. bioRxiv 2024. [Google Scholar] [CrossRef]
  95. Zhao, H.; Jin, Z.; Li, J.; Fang, J.; Wu, W. Disulfidptosis-Mediated Immune Microenvironment Regulation Characteristics in Atherosclerosis: Novel Insights from Bioinformatics Analyses; Research Square: Durham, NC, USA, 2024. [Google Scholar] [CrossRef]
  96. Rehman, S.; Gora, A.H.; Abdelhafiz, Y.; Dias, J.; Pierre, R.; Meynen, K.; Fernandes, J.M.; Sørensen, M.; Brugman, S.; Kiron, V. Potential of algae-derived alginate oligosaccharides and β-glucan to counter inflammation in adult zebrafish intestine. Front. Immunol. 2023, 14, 1183701. [Google Scholar] [CrossRef]
  97. He, N.; Yang, Y.; Wang, H.; Liu, N.; Yang, Z.; Li, S. Unsaturated alginate oligosaccharides (UAOS) protects against dextran sulfate sodium-induced colitis associated with regulation of gut microbiota. J. Funct. Foods 2021, 83, 104536. [Google Scholar] [CrossRef]
  98. Schmidt, C.; Berger, T.; Groettrup, M.; Basler, M. Immunoproteasome inhibition impairs T and B cell activation by restraining ERK signaling and proteostasis. Front. Immunol. 2018, 9, 2386. [Google Scholar] [CrossRef] [PubMed]
  99. Weißenberg, S.Y.; Szelinski, F.; Schrezenmeier, E.; Stefanski, A.L.; Wiedemann, A.; Rincon-Arevalo, H.; Welle, A.; Jungmann, A.; Nordström, K.; Walter, J.; et al. Identification and characterization of post-activated B cells in systemic autoimmune diseases. Front. Immunol. 2019, 10, 2136. [Google Scholar] [CrossRef]
  100. Valeff, N.; Muzzio, D.O.; Matzner, F.; Dibo, M.; Golchert, J.; Homuth, G.; Abba, M.C.; Zygmunt, M.; Jensen, F. B cells acquire a unique and differential transcriptomic profile during pregnancy. Genomics 2021, 113, 2614–2622. [Google Scholar] [CrossRef] [PubMed]
  101. Shi, Y.X.; Yan, J.H.; Liu, W.; Deng, J. Identifies KCTD5 as a novel cancer biomarker associated with programmed cell death and chemotherapy drug sensitivity. BMC Cancer 2023, 23, 408. [Google Scholar] [CrossRef]
  102. Lambert, C.; Preijers, F.W.; Yanikkaya Demirel, G.; Sack, U. Monocytes and macrophages in flow: An ESCCA initiative on advanced analyses of monocyte lineage using flow cytometry. Cytom. Part B Clin. Cytom. 2017, 92, 180–188. [Google Scholar] [CrossRef]
  103. Dietrich, M.A.; Adamek, M.; Teitge, F.; Teich, L.; Jung-Schroers, V.; Malinowska, A.; Świderska, B.; Rakus, K.; Kodzik, N.; Chadzińska, M.; et al. Proteomic analysis of carp seminal plasma provides insights into the immune response to bacterial infection of the male reproductive system. Fish Shellfish Immunol. 2022, 127, 822–835. [Google Scholar] [CrossRef] [PubMed]
  104. Grainger, S.; Richter, J.; Palazón, R.E.; Pouget, C.; Lonquich, B.; Wirth, S.; Grassme, K.S.; Herzog, W.; Swift, M.R.; Weinstein, B.M.; et al. Wnt9a is required for the aortic amplification of nascent hematopoietic stem cells. Cell Rep. 2016, 17, 1595–1606. [Google Scholar] [CrossRef]
  105. Malengier-Devlies, B.; Metzemaekers, M.; Wouters, C.; Proost, P.; Matthys, P. Neutrophil homeostasis and emergency granulopoiesis: The example of systemic juvenile idiopathic arthritis. Front. Immunol. 2021, 12, 766620. [Google Scholar] [CrossRef]
  106. Deng, Y.; Wang, H.; Liu, X.; Yuan, H.; Xu, J.; de Thé, H.; Zhou, J.; Zhu, J. Zbtb14 regulates monocyte and macrophage development through inhibiting pu. 1 expression in zebrafish. Elife 2022, 11, e80760. [Google Scholar] [CrossRef] [PubMed]
  107. Zhu, Z.; Zheng, L.; Li, Y.; Huang, T.; Chao, Y.C.; Pan, L.; Zhu, H.; Zhao, Y.; Yu, W.; Li, P. Potential immunotherapeutic targets on myeloid cells for neurovascular repair after ischemic stroke. Front. Neurosci. 2019, 13, 758. [Google Scholar] [CrossRef] [PubMed]
  108. Thind, M.K.; Uhlig, H.H.; Glogauer, M.; Palaniyar, N.; Bourdon, C.; Gwela, A.; Lancioni, C.L.; Berkley, J.A.; Bandsma, R.H.; Farooqui, A. A metabolic perspective of the neutrophil life cycle: New avenues in immunometabolism. Front. Immunol. 2024, 14, 1334205. [Google Scholar] [CrossRef]
  109. Zhang, Q.; Zhu, B.; Yang, H.; Li, F.; Qu, Y.; Lu, L.; Zhang, Q. Exploration of YBX1 role in the prognostic value and immune characteristics by single-cell and bulk sequencing analysis for liver hepatocellular carcinoma. J. Gene Med. 2024, 26, e3680. [Google Scholar] [CrossRef]
  110. Zhang, Y.; Cen, J.; Jia, Z.; Hsiao, C.D.; Xia, Q.; Wang, X.; Chen, X.; Wang, R.; Jiang, Z.; Zhang, L.; et al. Hepatotoxicity induced by isoniazid-lipopolysaccharide through endoplasmic reticulum stress, autophagy, and apoptosis pathways in zebrafish. Antimicrob. Agents Chemother. 2019, 63, 10–1128. [Google Scholar] [CrossRef] [PubMed]
  111. Ding, Q.; Wang, H. Bisphenol AP Exposure-Induced Inflammatory Response in Zebrafish Embryos Revealed by Transcriptomic Analysis. Environ. Health 2024, 2, 373–380. [Google Scholar] [CrossRef]
  112. Lv, J.; Zhou, Y.; Zhou, N.; Wang, Z.; Chen, J.; Chen, H.; Wang, D.; Zhou, L.; Wei, K.; Zhang, H.; et al. Epigenetic modification of CSDE1 locus dictates immune recognition of nascent tumorigenic cells. Sci. Transl. Med. 2023, 15, eabq6024. [Google Scholar] [CrossRef]
  113. Droho, S.; Voigt, A.P.; Sterling, J.K.; Rajesh, A.; Chan, K.S.; Cuda, C.M.; Perlman, H.; Lavine, J.A. NR4A1 deletion promotes pro-angiogenic polarization of macrophages derived from classical monocytes in a mouse model of neovascular age-related macular degeneration. J. Neuroinflammation 2023, 20, 238. [Google Scholar] [CrossRef]
  114. Du, H.; Sun, J.; Wang, X.; Zhao, L.; Liu, X.; Zhang, C.; Wang, F.; Wu, J. FOSL2-mediated transcription of ISG20 induces M2 polarization of macrophages and enhances tumorigenic ability of glioblastoma cells. J. Neuro-Oncol. 2024, 169, 659–670. [Google Scholar] [CrossRef]
  115. Liu, B.; Chen, P.; Xi, D.; Zhu, H.; Gao, Y. ATF4 regulates CCL2 expression to promote endometrial cancer growth by controlling macrophage infiltration. Exp. Cell Res. 2017, 360, 105–112. [Google Scholar] [CrossRef] [PubMed]
  116. Kobayashi, E.H.; Suzuki, T.; Funayama, R.; Nagashima, T.; Hayashi, M.; Sekine, H.; Tanaka, N.; Moriguchi, T.; Motohashi, H.; Nakayama, K.; et al. Nrf2 suppresses macrophage inflammatory response by blocking proinflammatory cytokine transcription. Nat. Commun. 2016, 7, 11624. [Google Scholar] [CrossRef]
  117. Shang, S.; Yang, C.; Chen, F.; Xiang, R.S.; Zhang, H.; Dai, S.Y.; Liu, J.; Lv, X.X.; Zhang, C.; Liu, X.T.; et al. ID1 expressing macrophages support cancer cell stemness and limit CD8+ T cell infiltration in colorectal cancer. Nat. Commun. 2023, 14, 7661. [Google Scholar] [CrossRef] [PubMed]
  118. Denans, N.; Tran, N.T.; Swall, M.E.; Diaz, D.C.; Blanck, J.; Piotrowski, T. An anti-inflammatory activation sequence governs macrophage transcriptional dynamics during tissue injury in zebrafish. Nat. Commun. 2022, 13, 5356. [Google Scholar] [CrossRef]
  119. Li, Y.; Chen, L.; Li, Y.; Deng, P.; Yang, C.; Li, Y.; Liao, L.; Zhu, Z.; Wang, Y.; Huang, R. miR-2188-5p promotes GCRV replication by the targeted degradation of klf2a in Ctenopharyngodon idellus. Dev. Comp. Immunol. 2023, 138, 104516. [Google Scholar] [CrossRef] [PubMed]
  120. Cao, H. Lipopolysaccharide regulation of antiinflammatory tristetraprolin family and proinflammatory gene expression in mouse macrophages. BMC Res. Notes 2024, 17, 82. [Google Scholar] [CrossRef] [PubMed]
  121. Vegliante, M.C.; Mazzara, S.; Zaccaria, G.M.; De Summa, S.; Esposito, F.; Melle, F.; Motta, G.; Sapienza, M.R.; Opinto, G.; Volpe, G.; et al. NR1H3 (LXRα) is associated with pro-inflammatory macrophages, predicts survival and suggests potential therapeutic rationales in diffuse large b-cell lymphoma. Hematol. Oncol. 2022, 40, 864–875. [Google Scholar] [CrossRef]
  122. Dinglong, Y.; Shuai, C.; Yujing, C.; Beiyang, W.; Guohao, Z.; Zhiqiang, Z. Identification of Key Genes and Competitive Endogenous RNA Network Construction in Osteonecrosis of the Femoral Head by Integrated Bioinformatics Analysis; Research Square: Durham, NC, USA, 2021. [Google Scholar] [CrossRef]
  123. Ren, J.; Yan, D.; Wang, Y.; Zhang, J.; Li, M.; Xiong, W.; Jing, X.; Li, P.; Zhao, W.; Xiong, X.; et al. Inhibitor of differentiation-2 protein ameliorates DSS-induced ulcerative colitis by inhibiting NF-κB activation in neutrophils. Front. Immunol. 2021, 12, 760999. [Google Scholar] [CrossRef] [PubMed]
  124. Rodríguez-Ruiz, L.; Lozano-Gil, J.M.; Naranjo-Sánchez, E.; Martínez-Balsalobre, E.; Martínez-López, A.; Lachaud, C.; Blanquer, M.; Phung, T.K.; García-Moreno, D.; Cayuela, M.L.; et al. ZAKα/P38 kinase signaling pathway regulates hematopoiesis by activating the NLRP1 inflammasome. EMBO Mol. Med. 2023, 15, e18142. [Google Scholar] [CrossRef]
  125. Hou, Y.; Khatri, P.; Rindy, J.; Schultz, Z.; Gao, A.; Chen, Z.; Gibson, A.L.; Huttenlocher, A.; Dinh, H.Q. Single-cell transcriptional landscape of temporal neutrophil response to burn wound in larval zebrafish. bioRxiv 2024. [Google Scholar] [CrossRef]
  126. Brown, G. Retinoic acid receptor regulation of decision-making for cell differentiation. Front. Cell Dev. Biol. 2023, 11, 1182204. [Google Scholar] [CrossRef]
  127. Rauschmeier, R.; Reinhardt, A.; Gustafsson, C.; Glaros, V.; Artemov, A.V.; Dunst, J.; Taneja, R.; Adameyko, I.; Månsson, R.; Busslinger, M.; et al. Bhlhe40 function in activated B and TFH cells restrains the GC reaction and prevents lymphomagenesis. J. Exp. Med. 2021, 219, e20211406. [Google Scholar] [CrossRef] [PubMed]
  128. Wiehagen, K.R.; Corbo-Rodgers, E.; Li, S.; Staub, E.S.; Hunter, C.A.; Morrisey, E.E.; Maltzman, J.S. Foxp4 is dispensable for T cell development, but required for robust recall responses. PLoS ONE 2012, 7, e42273. [Google Scholar] [CrossRef] [PubMed]
  129. Meijer, A.H.; van der Sar, A.M.; Cunha, C.; Lamers, G.E.; Laplante, M.A.; Kikuta, H.; Bitter, W.; Becker, T.S.; Spaink, H.P. Identification and real-time imaging of a myc-expressing neutrophil population involved in inflammation and mycobacterial granuloma formation in zebrafish. Dev. Comp. Immunol. 2008, 32, 36–49. [Google Scholar] [CrossRef] [PubMed]
  130. Sukhorukov, V.N.; Khotina, V.A.; Bagheri Ekta, M.; Ivanova, E.A.; Sobenin, I.A.; Orekhov, A.N. Endoplasmic reticulum stress in macrophages: The vicious circle of lipid accumulation and pro-inflammatory response. Biomedicines 2020, 8, 210. [Google Scholar] [CrossRef]
  131. Lai, H.C.; James, A.; Luff, J.; De Souza, P.; Quek, H.; Ho, U.; Lavin, M.F.; Roberts, T.L. Regulation of RNA degradation pathways during the lipopolysaccharide response in macrophages. J. Leucoc. Biol. 2021, 109, 593–603. [Google Scholar] [CrossRef] [PubMed]
  132. Ran, L.; Zhang, S.; Wang, G.; Zhao, P.; Sun, J.; Zhou, J.; Gan, H.; Jeon, R.; Li, Q.; Herrmann, J.; et al. Mitochondrial pyruvate carrier-mediated metabolism is dispensable for the classical activation of macrophages. Nat. Metab. 2023, 5, 804–820. [Google Scholar] [CrossRef] [PubMed]
  133. Shan, X.; Hu, P.; Ni, L.; Shen, L.; Zhang, Y.; Ji, Z.; Cui, Y.; Guo, M.; Wang, H.; Ran, L.; et al. Serine metabolism orchestrates macrophage polarization by regulating the IGF1–p38 axis. Cell. Mol. Immunol. 2022, 19, 1263–1278. [Google Scholar] [CrossRef]
  134. Yang, D.; Yang, M.; Yin, Y.; Kou, T.; Peng, L.; Chen, Z.; Zheng, J.; Peng, B. Serine metabolism tunes immune responses to promote Oreochromis niloticus survival upon Edwardsiella tarda infection. mSystems 2021, 6, e0042621. [Google Scholar] [CrossRef]
  135. Wculek, S.K.; Heras-Murillo, I.; Mastrangelo, A.; Mañanes, D.; Galán, M.; Miguel, V.; Curtabbi, A.; Barbas, C.; Chandel, N.S.; Enríquez, J.A.; et al. Oxidative phosphorylation selectively orchestrates tissue macrophage homeostasis. Immunity 2023, 56, 516–530. [Google Scholar] [CrossRef]
  136. Lv, P.; Liu, F. Heme-deficient primitive red blood cells induce HSPC ferroptosis by altering iron homeostasis during zebrafish embryogenesis. Development 2023, 150, dev201690. [Google Scholar] [CrossRef] [PubMed]
  137. Yang, Y.; Wang, Y.; Guo, L.; Gao, W.; Tang, T.L.; Yan, M. Interaction between macrophages and ferroptosis. Cell Death Dis. 2022, 13, 355. [Google Scholar] [CrossRef]
  138. Antoine, T.E.; Jones, K.S.; Dale, R.M.; Shukla, D.; Tiwari, V. Zebrafish: Modeling for herpes simplex virus infections. Zebrafish 2014, 11, 17–25. [Google Scholar] [CrossRef]
  139. Lee, D.H.; Ghiasi, H. Roles of M1 and M2 macrophages in herpes simplex virus 1 infectivity. J. Virol. 2017, 91, 10–1128. [Google Scholar] [CrossRef] [PubMed]
  140. Lang, J.; Bohn, P.; Bhat, H.; Jastrow, H.; Walkenfort, B.; Cansiz, F.; Fink, J.; Bauer, M.; Olszewski, D.; Ramos-Nascimento, A.; et al. Acid ceramidase of macrophages traps herpes simplex virus in multivesicular bodies and protects from severe disease. Nat. Commun. 2020, 11, 1338. [Google Scholar] [CrossRef] [PubMed]
  141. Rajput, V.B.; Karthikeyan, M.; Ramasamy, S. Zebrafish acid ceramidase: Expression in Pichia pastoris GS115 and biochemical characterization. Int. J. Biol. Macromol. 2019, 122, 587–593. [Google Scholar] [CrossRef]
  142. Almeida-da Silva, C.L.C.; Savio, L.E.B.; Coutinho-Silva, R.; Ojcius, D.M. The role of NOD-like receptors in innate immunity. Front. Immunol. 2023, 14, 1122586. [Google Scholar] [CrossRef]
  143. Watts, C. Lysosomes and lysosome-related organelles in immune responses. FEBS Open Bio 2022, 12, 678–693. [Google Scholar] [CrossRef] [PubMed]
  144. Lee, H.J.; Woo, Y.; Hahn, T.W.; Jung, Y.M.; Jung, Y.J. Formation and maturation of the phagosome: A key mechanism in innate immunity against intracellular bacterial infection. Microorganisms 2020, 8, 1298. [Google Scholar] [CrossRef]
  145. Kemp, M.G. Crosstalk between apoptosis and autophagy: Environmental genotoxins, infection, and innate immunity. J. Cell Death 2017, 10, 1179670716685085. [Google Scholar] [CrossRef] [PubMed]
  146. Sadiku, P.; Willson, J.A.; Ryan, E.M.; Sammut, D.; Coelho, P.; Watts, E.R.; Grecian, R.; Young, J.M.; Bewley, M.; Arienti, S.; et al. Neutrophils fuel effective immune responses through gluconeogenesis and glycogenesis. Cell Metab. 2021, 33, 411–423. [Google Scholar] [CrossRef]
  147. Fang, X.; Ma, L.; Wang, Y.; Ren, F.; Yu, Y.; Yuan, Z.; Wei, H.; Zhang, H.; Sun, Y. Neutrophil extracellular traps accelerate vascular smooth muscle cell proliferation via Akt/CDKN1b/TK1 accompanying with the occurrence of hypertension. J. Hypertens. 2022, 40, 2045–2057. [Google Scholar] [CrossRef] [PubMed]
  148. Wirrig, C.; McKean, J.S.; Wilson, H.M.; Nixon, G.F. Sphingosylphosphorylcholine inhibits macrophage adhesion to vascular smooth muscle cells. Biochem. Pharmacol. 2016, 115, 43–50. [Google Scholar] [CrossRef]
  149. Bonetti, J.; Corti, A.; Lerouge, L.; Pompella, A.; Gaucher, C. Phenotypic Modulation of Macrophages and Vascular Smooth Muscle Cells in Atherosclerosis—Nitro-Redox Interconnections. Antioxidants 2021, 10, 516. [Google Scholar] [CrossRef]
  150. Bruton, F.A.; Kaveh, A.; Ross-Stewart, K.M.; Matrone, G.; Oremek, M.E.; Solomonidis, E.G.; Tucker, C.S.; Mullins, J.J.; Lucas, C.D.; Brittan, M.; et al. Macrophages trigger cardiomyocyte proliferation by increasing epicardial vegfaa expression during larval zebrafish heart regeneration. Dev. Cell 2022, 57, 1512–1528. [Google Scholar] [CrossRef]
  151. Lee, Y.S.; Kang, S.U.; Lee, M.H.; Kim, H.J.; Han, C.H.; Won, H.R.; Park, Y.U.; Kim, C.H. GnRH impairs diabetic wound healing through enhanced NETosis. Cell. Mol. Immunol. 2020, 17, 856–864. [Google Scholar] [CrossRef]
  152. Morrison, T.; Watts, E.R.; Sadiku, P.; Walmsley, S.R. The emerging role for metabolism in fueling neutrophilic inflammation. Immunol. Rev. 2023, 314, 427–441. [Google Scholar] [CrossRef] [PubMed]
  153. Guerrero-Juarez, C.F.; Schilf, P.; Li, J.; Zappia, M.P.; Bao, L.; Patel, P.M.; Gieseler-Tillmann, J.; Murthy, S.; Cole, C.; Sverdlov, M.; et al. C-type lectin receptor expression is a hallmark of neutrophils infiltrating the skin in epidermolysis bullosa acquisita. Front. Immunol. 2023, 14, 1266359. [Google Scholar] [CrossRef]
  154. Parkos, C.A. Neutrophil-epithelial interactions: A double-edged sword. Am. J. Pathol. 2016, 186, 1404–1416. [Google Scholar] [CrossRef]
  155. Strizova, Z.; Benesova, I.; Bartolini, R.; Novysedlak, R.; Cecrdlova, E.; Foley, L.K.; Striz, I. M1/M2 macrophages and their overlaps–myth or reality? Clin. Sci. 2023, 137, 1067–1093. [Google Scholar] [CrossRef]
  156. Jones, H.R.; Robb, C.T.; Perretti, M.; Rossi, A.G. The role of neutrophils in inflammation resolution. In Proceedings of the Seminars in Immunology; Elsevier: Amsterdam, The Netherlands, 2016; Volume 28, pp. 137–145. [Google Scholar] [CrossRef]
  157. Nguyen-Chi, M.; Laplace-Builhe, B.; Travnickova, J.; Luz-Crawford, P.; Tejedor, G.; Phan, Q.T.; Duroux-Richard, I.; Levraud, J.P.; Kissa, K.; Lutfalla, G.; et al. Identification of polarized macrophage subsets in zebrafish. Elife 2015, 4, e07288. [Google Scholar] [CrossRef]
  158. Yang, L.; Zhou, X.; Huang, W.; Fang, Q.; Hu, J.; Yu, L.; Ma, N.; Zhang, W. Protective effect of phillyrin on lethal LPS-induced neutrophil inflammation in zebrafish. Cell. Physiol. Biochem. 2018, 43, 2074–2087. [Google Scholar] [CrossRef]
  159. Xie, Y.; Meijer, A.H.; Schaaf, M.J. Modeling inflammation in zebrafish for the development of anti-inflammatory drugs. Front. Cell Dev. Biol. 2021, 8, 620984. [Google Scholar] [CrossRef] [PubMed]
  160. Kolonics, F.; Kajdácsi, E.; Farkas, V.J.; Veres, D.S.; Khamari, D.; Kittel, Á.; Merchant, M.L.; McLeish, K.R.; Lõrincz, Á.M.; Ligeti, E. Neutrophils produce proinflammatory or anti-inflammatory extracellular vesicles depending on the environmental conditions. J. Leucoc. Biol. 2021, 109, 793–806. [Google Scholar] [CrossRef] [PubMed]
  161. Bobrovskikh, A.; Doroshkov, A.; Mazzoleni, S.; Cartenì, F.; Giannino, F.; Zubairova, U. A sight on single-cell transcriptomics in plants through the prism of cell-based computational modeling approaches: Benefits and challenges for data analysis. Front. Genet. 2021, 12, 652974. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Key steps of the conducted scRNA-seq meta-analysis. The first stage of the analysis consisted of searching and selecting data from 13 experiments from the GEO NCBI database in the form of count matrices or Seurat objects. The second stage consisted of bringing the data into a uniform format, combining individual subsets of experiments into single matrices, and converting gene identifiers into a GRCz11 format using Genekitr [22]. Next, the procedures for integrating the 13 matrices, normalizing them, harmonizing them to eliminate the batch effect, and clustering were carried out using Seurat v. 5 [23] and Harmony [24]. The fourth stage consisted of filtering out non-immune cells, cells with low coverage, and genes irrelevant to immune cells. After filtering, the harmonization and clustering procedures were repeated on the final dataset. The fifth stage was to determine the cell types of immune cells in the final dataset using the main known marker genes from the recent work [8]. Transcription factors were identified according to AnimalTFDB4 [25]. Next, in stage six, the marker genes of the identified cell types were calculated and classified by expression level and fold changes in target types. Stages 7–9 consisted of reconstructing gene networks for macrophage- and neutrophil-specific genes using the STRING-db v. 12.0 database and co-expression data calculated from our final dataset using the hdWGCNA [26] package. Network layout and visualization were performed using Cytoscape v. 3.10.2 [27].
Figure 1. Key steps of the conducted scRNA-seq meta-analysis. The first stage of the analysis consisted of searching and selecting data from 13 experiments from the GEO NCBI database in the form of count matrices or Seurat objects. The second stage consisted of bringing the data into a uniform format, combining individual subsets of experiments into single matrices, and converting gene identifiers into a GRCz11 format using Genekitr [22]. Next, the procedures for integrating the 13 matrices, normalizing them, harmonizing them to eliminate the batch effect, and clustering were carried out using Seurat v. 5 [23] and Harmony [24]. The fourth stage consisted of filtering out non-immune cells, cells with low coverage, and genes irrelevant to immune cells. After filtering, the harmonization and clustering procedures were repeated on the final dataset. The fifth stage was to determine the cell types of immune cells in the final dataset using the main known marker genes from the recent work [8]. Transcription factors were identified according to AnimalTFDB4 [25]. Next, in stage six, the marker genes of the identified cell types were calculated and classified by expression level and fold changes in target types. Stages 7–9 consisted of reconstructing gene networks for macrophage- and neutrophil-specific genes using the STRING-db v. 12.0 database and co-expression data calculated from our final dataset using the hdWGCNA [26] package. Network layout and visualization were performed using Cytoscape v. 3.10.2 [27].
Biology 13 00773 g001
Figure 2. Primary steps of zebrafish kidney marrow scRNA-seq dataset integration and IIS cell classification. (A) Sankey diagram showing the total number of zebrafish kidney cells pooled from each of the 13 selected experiments. The initial analysis included 382,295 cells that passed the soft filtering threshold (nFeatureRNA greater than 200 and the percent of reads mapped to mitochondrial genes less than 10%) and were further classified into immune (a total of 254,256 cells) and non-immune (a total of 128,039 cells) categories. The final analysis included 94,398 immune cells from the selected group, characterized by high quality and satisfying the following criteria: 3000 ≤ nCountRNA ≤ 30,000 and 400 ≤ nFeatureRNA ≤ 4000. (B) Distribution of the selected 94,398 cells within the original dataset of 382,295 cells. The color indicates cluster characteristics after filtration. Cells from clusters marked in black were excluded from further analysis as they were not immune, corresponding to the bottom stream in the Sankey diagram (A). Clusters with less than 20 percent of cells remaining after filtration are marked in gray. Clusters with 20 to 40 percent of cells remaining are marked in pink. Clusters with 40 to 60 percent of cells remaining are marked in red. Clusters with more than 60 percent of cells remaining are marked in dark red. Cell types were labeled according to the key markers from [8]. Expression plots of the main marker genes are provided in Supplementary Figure S1. (C) Total number of cells selected from individual experiments and their cell-type composition. Stacked bar charts display the overall composition of the datasets by cell type. HSCs refer to hematopoietic stem cells.
Figure 2. Primary steps of zebrafish kidney marrow scRNA-seq dataset integration and IIS cell classification. (A) Sankey diagram showing the total number of zebrafish kidney cells pooled from each of the 13 selected experiments. The initial analysis included 382,295 cells that passed the soft filtering threshold (nFeatureRNA greater than 200 and the percent of reads mapped to mitochondrial genes less than 10%) and were further classified into immune (a total of 254,256 cells) and non-immune (a total of 128,039 cells) categories. The final analysis included 94,398 immune cells from the selected group, characterized by high quality and satisfying the following criteria: 3000 ≤ nCountRNA ≤ 30,000 and 400 ≤ nFeatureRNA ≤ 4000. (B) Distribution of the selected 94,398 cells within the original dataset of 382,295 cells. The color indicates cluster characteristics after filtration. Cells from clusters marked in black were excluded from further analysis as they were not immune, corresponding to the bottom stream in the Sankey diagram (A). Clusters with less than 20 percent of cells remaining after filtration are marked in gray. Clusters with 20 to 40 percent of cells remaining are marked in pink. Clusters with 40 to 60 percent of cells remaining are marked in red. Clusters with more than 60 percent of cells remaining are marked in dark red. Cell types were labeled according to the key markers from [8]. Expression plots of the main marker genes are provided in Supplementary Figure S1. (C) Total number of cells selected from individual experiments and their cell-type composition. Stacked bar charts display the overall composition of the datasets by cell type. HSCs refer to hematopoietic stem cells.
Biology 13 00773 g002
Figure 3. General information on the final dataset of the selected IIS cells. (A) Distribution of the selected immune cells, totaling 94,398, in the harmonized UMAP space containing 382,295 cells. The selected immune cells are marked in red within the original dataset space. (B) Distribution of selected zebrafish kidney immune cells in the harmonized UMAP space of immune-specific differentially expressed genes, comprising 94,398 cells. A total of 44 identified Seurat clusters are numbered in descending order based on their cell count. The areas corresponding to the identified cell types are outlined by a dashed line. HSCs refer to hematopoietic stem cells. (C) Hierarchical clustering dendrogram of immune cell clusters based on the expression of immune-specific differentially expressed genes. This clustering is based on the correlating gene expression of 3737 immune genes between Seurat clusters. The separated subgroups of the same cell types are numbered with (1)–(3). (D) Distribution of the expression of key marker genes, selected based on [8], used to annotate the main immune cell types.
Figure 3. General information on the final dataset of the selected IIS cells. (A) Distribution of the selected immune cells, totaling 94,398, in the harmonized UMAP space containing 382,295 cells. The selected immune cells are marked in red within the original dataset space. (B) Distribution of selected zebrafish kidney immune cells in the harmonized UMAP space of immune-specific differentially expressed genes, comprising 94,398 cells. A total of 44 identified Seurat clusters are numbered in descending order based on their cell count. The areas corresponding to the identified cell types are outlined by a dashed line. HSCs refer to hematopoietic stem cells. (C) Hierarchical clustering dendrogram of immune cell clusters based on the expression of immune-specific differentially expressed genes. This clustering is based on the correlating gene expression of 3737 immune genes between Seurat clusters. The separated subgroups of the same cell types are numbered with (1)–(3). (D) Distribution of the expression of key marker genes, selected based on [8], used to annotate the main immune cell types.
Biology 13 00773 g003
Figure 4. The expression heatmap for identified cell-type-specific highly expressed marker genes in different cell types. The heatmap displays highly expressed cell-type-specific genes with log2FC ≥ 2 and average counts-per-cell (ACPC) ≥ 2 in the target cell type. The expression values of genes in 44 clusters, grouped by cell types, are shown using a blue–red color scheme. A colored bar to the right of the heatmap indicates the cell type to which the gene is specific. The cell type labels are provided at the bottom of the heatmap. Red asterisks mark known marker genes for this cell type, as reported by [8]. Common names for Danio rerio genes are provided; genes with high-confidence human orthologs are marked in bold blue font. (A) Identified 62 marker genes for macrophages and macrophage-like cells. (B) Identified 62 marker genes for mature neutrophils, immature neutrophils, and neutrophil-like cells. (C) Identified 18 marker genes specific to other immune cell types (B cells, monocytes, and hematopoietic stem cells).
Figure 4. The expression heatmap for identified cell-type-specific highly expressed marker genes in different cell types. The heatmap displays highly expressed cell-type-specific genes with log2FC ≥ 2 and average counts-per-cell (ACPC) ≥ 2 in the target cell type. The expression values of genes in 44 clusters, grouped by cell types, are shown using a blue–red color scheme. A colored bar to the right of the heatmap indicates the cell type to which the gene is specific. The cell type labels are provided at the bottom of the heatmap. Red asterisks mark known marker genes for this cell type, as reported by [8]. Common names for Danio rerio genes are provided; genes with high-confidence human orthologs are marked in bold blue font. (A) Identified 62 marker genes for macrophages and macrophage-like cells. (B) Identified 62 marker genes for mature neutrophils, immature neutrophils, and neutrophil-like cells. (C) Identified 18 marker genes specific to other immune cell types (B cells, monocytes, and hematopoietic stem cells).
Biology 13 00773 g004
Figure 5. Heatmap of identified transcription factors (TFs) expressed in at least one immune cell cluster in zebrafish kidney (ACPC 1 ). Gene names are listed on the left, with those having the Homo sapiens ortholog(s) with high confidence (according to Ensembl Orthologues) highlighted in bold blue. Genes that were previously mentioned in related papers [105,106,107,108] are marked with plus signs with certain colors. The expression values of genes across 44 clusters, grouped by cell types, are depicted using a blue–red color scheme. TFs are categorized based on their expression patterns into sections from top to bottom: common (expressed in multiple identified immune cell types), macrophage-specific, neutrophil-specific, and TFs of other cell types. Cell type labels are provided at the bottom of the heatmaps.
Figure 5. Heatmap of identified transcription factors (TFs) expressed in at least one immune cell cluster in zebrafish kidney (ACPC 1 ). Gene names are listed on the left, with those having the Homo sapiens ortholog(s) with high confidence (according to Ensembl Orthologues) highlighted in bold blue. Genes that were previously mentioned in related papers [105,106,107,108] are marked with plus signs with certain colors. The expression values of genes across 44 clusters, grouped by cell types, are depicted using a blue–red color scheme. TFs are categorized based on their expression patterns into sections from top to bottom: common (expressed in multiple identified immune cell types), macrophage-specific, neutrophil-specific, and TFs of other cell types. Cell type labels are provided at the bottom of the heatmaps.
Biology 13 00773 g005
Figure 6. Reconstructed gene network of highly expressed genes in the population of zebrafish kidney macrophages. Nodes in the network represent either transcription factors (hexagons) or other protein-coding genes (ellipses). The red fill intensity of the nodes indicates the level of ACPC expression within the macrophage population. Ellipses are highlighted to denote highly specific macrophage genes (log2FC   2 ): those outlined in black represent key macrophage marker genes confirmed by [8], while those outlined in blue denote potential new macrophage markers. Transcription factors previously identified in the literature [105,106,107,108] are outlined in black, whereas those identified in our analysis are outlined in blue. The color of the connections represents the type of interaction: black lines indicate protein–protein associations based on STRING-db v. 12.0 (minimum required interaction score: 0.4 (medium); sources of active interactions: experiments and databases); gray lines denote co-expression interactions predicted by our analysis using hdWGCNA (minimum co-expression score = 0.1); red lines represent protein–protein associations between transcription factors and protein-coding genes; pink lines indicate co-expression among transcription factors and protein-coding genes.
Figure 6. Reconstructed gene network of highly expressed genes in the population of zebrafish kidney macrophages. Nodes in the network represent either transcription factors (hexagons) or other protein-coding genes (ellipses). The red fill intensity of the nodes indicates the level of ACPC expression within the macrophage population. Ellipses are highlighted to denote highly specific macrophage genes (log2FC   2 ): those outlined in black represent key macrophage marker genes confirmed by [8], while those outlined in blue denote potential new macrophage markers. Transcription factors previously identified in the literature [105,106,107,108] are outlined in black, whereas those identified in our analysis are outlined in blue. The color of the connections represents the type of interaction: black lines indicate protein–protein associations based on STRING-db v. 12.0 (minimum required interaction score: 0.4 (medium); sources of active interactions: experiments and databases); gray lines denote co-expression interactions predicted by our analysis using hdWGCNA (minimum co-expression score = 0.1); red lines represent protein–protein associations between transcription factors and protein-coding genes; pink lines indicate co-expression among transcription factors and protein-coding genes.
Biology 13 00773 g006
Figure 7. Reconstructed gene network of highly expressed genes in the population of zebrafish kidney neutrophils. Nodes in the network represent either transcription factors (hexagons) or other protein-coding genes (ellipses). The intensity of the red fill in the nodes reflects the level of ACPC expression within the neutrophil population. Ellipses are highlighted to denote highly specific neutrophil genes (log2FC   2 ): those outlined in black are key neutrophil marker genes confirmed by [8], while those outlined in blue represent potential new neutrophil markers. Transcription factors previously identified in the literature [105,106,107,108] are outlined in black, while those identified in our analysis are outlined in blue. The color of the connections indicates the type of interaction: black lines represent protein–protein associations based on STRING-db (minimum required interaction score: 0.4 (medium); sources of active interactions include experiments and databases); gray lines denote co-expression interactions predicted by our analysis using hdWGCNA (minimum co-expression score = 0.1); red lines indicate protein–protein associations between transcription factors and protein-coding genes; pink lines represent co-expression among transcription factors and protein-coding genes.
Figure 7. Reconstructed gene network of highly expressed genes in the population of zebrafish kidney neutrophils. Nodes in the network represent either transcription factors (hexagons) or other protein-coding genes (ellipses). The intensity of the red fill in the nodes reflects the level of ACPC expression within the neutrophil population. Ellipses are highlighted to denote highly specific neutrophil genes (log2FC   2 ): those outlined in black are key neutrophil marker genes confirmed by [8], while those outlined in blue represent potential new neutrophil markers. Transcription factors previously identified in the literature [105,106,107,108] are outlined in black, while those identified in our analysis are outlined in blue. The color of the connections indicates the type of interaction: black lines represent protein–protein associations based on STRING-db (minimum required interaction score: 0.4 (medium); sources of active interactions include experiments and databases); gray lines denote co-expression interactions predicted by our analysis using hdWGCNA (minimum co-expression score = 0.1); red lines indicate protein–protein associations between transcription factors and protein-coding genes; pink lines represent co-expression among transcription factors and protein-coding genes.
Biology 13 00773 g007
Figure 8. Specific KEGG pathways enriched in macrophages and neutrophils. Green bars and green genes represent macrophage-specific pathways, while orange bars and orange genes represent neutrophil-specific pathways. Significant enrichment within these pathways is indicated by white asterisks. The names of the genes included in each pathway are displayed above the corresponding bars. (A) Macrophage-specific enriched pathways; (B) commonly enriched pathways; (C) neutrophil-specific enriched pathways.
Figure 8. Specific KEGG pathways enriched in macrophages and neutrophils. Green bars and green genes represent macrophage-specific pathways, while orange bars and orange genes represent neutrophil-specific pathways. Significant enrichment within these pathways is indicated by white asterisks. The names of the genes included in each pathway are displayed above the corresponding bars. (A) Macrophage-specific enriched pathways; (B) commonly enriched pathways; (C) neutrophil-specific enriched pathways.
Biology 13 00773 g008
Table 1. Single-cell transcriptomics experiments of zebrafish kidney used in the meta-analysis.
Table 1. Single-cell transcriptomics experiments of zebrafish kidney used in the meta-analysis.
GEO NCBI ID 1AgeGenotype/StrainNumber of Used SamplesInitial Number of CellsRef.
GSE1009103–9 monthsWT; prkdc611,327[29]
GSE112438not availableAB; CD41:GFP3713,824[30]
GSE1304874–12 monthsWT120,000[31]
GSE1503738 monthsWT; runx1839,424[32]
GSE1512314 monthsWT; gata2b914,463[33]
GSE166646adultWT16422[34]
GSE1760368 monthsrunx1835,178[32]
GSE1794012 monthsWT; rag; rag il2rga947,832[35]
GSE1907944 monthsGESTALT1051,540[36]
GSE191029adultWT; prkcda; cxcl8820,695[37]
GSE2421331–1.5 monthsAB336,600[8]
GSE2460393 monthsWT; cebpb10105,194[38]
GSE2527886 monthscebpb228,534[9]
1 Detailed information about these experiments is provided in Supplementary Table S6.
Table 2. IIS cell types and their presence in Seurat clusters identified in our meta-analysis.
Table 2. IIS cell types and their presence in Seurat clusters identified in our meta-analysis.
Cell TypeTotal Cell NumberSeurat Cluster NoCell Number in Cluster
HSCs11,01307852
161736
251202
41223
Macrophages15,36525831
45232
103149
32877
40276
Macrophage-like cells10,595151769
201600
211552
231545
27992
30910
31879
34819
37408
43121
Neutrophils26,55616372
54812
64466
93506
113095
191684
221550
36486
38403
42182
Immature neutrophils747735383
142094
Neutrophil-like cells3876123035
33841
B cells541684137
28982
39297
T cells567874209
241469
Monocytes8422132334
171710
181705
261019
29960
35697
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bobrovskikh, A.V.; Zubairova, U.S.; Naumenko, L.G.; Doroshkov, A.V. Catching the Big Fish in Big Data: A Meta-Analysis of Zebrafish Kidney scRNA-Seq Datasets Highlights Conserved Molecular Profiles of Macrophages and Neutrophils in Vertebrates. Biology 2024, 13, 773. https://doi.org/10.3390/biology13100773

AMA Style

Bobrovskikh AV, Zubairova US, Naumenko LG, Doroshkov AV. Catching the Big Fish in Big Data: A Meta-Analysis of Zebrafish Kidney scRNA-Seq Datasets Highlights Conserved Molecular Profiles of Macrophages and Neutrophils in Vertebrates. Biology. 2024; 13(10):773. https://doi.org/10.3390/biology13100773

Chicago/Turabian Style

Bobrovskikh, Aleksandr V., Ulyana S. Zubairova, Ludmila G. Naumenko, and Alexey V. Doroshkov. 2024. "Catching the Big Fish in Big Data: A Meta-Analysis of Zebrafish Kidney scRNA-Seq Datasets Highlights Conserved Molecular Profiles of Macrophages and Neutrophils in Vertebrates" Biology 13, no. 10: 773. https://doi.org/10.3390/biology13100773

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop