Per- and polyfluoroalkyl substances alter innate immune function: evidence and data gaps.

Per- and polyfluoroalkyl substances (PFASs) are a large class of compounds used in a variety of processes and consumer products. Their unique chemical properties make them ubiquitous and persistent environmental contaminants while also making them economically viable and socially convenient. To date, several reviews have been published to synthesize information regarding the immunotoxic effects of PFASs on the adaptive immune system. However, these reviews often do not include data on the impact of these compounds on innate immunity. Here, current literature is reviewed to identify and incorporate data regarding the effects of PFASs on innate immunity in humans, experimental models, and wildlife. Known mechanisms by which PFASs modulate innate immune function are also reviewed, including disruption of cell signaling, metabolism, and tissue-level effects. For PFASs where innate immune data are available, results are equivocal, raising additional questions about common mechanisms or pathways of toxicity, but highlighting that the innate immune system within several species can be perturbed by exposure to PFASs. Recommendations are provided for future research to inform hazard identification, risk assessment, and risk management practices for PFASs to protect the immune systems of exposed organisms as well as environmental health.

Automated, high-throughput quantification of EGFP-expressing neutrophils in zebrafish by machine learning and a highly-parallelized microscope.

Normal development of the immune system is essential for overall health and disease resistance. Bony fish, such as the zebrafish (Danio rerio), possess all the major immune cell lineages as mammals and can be employed to model human host response to immune challenge. Zebrafish neutrophils, for example, are present in the transparent larvae as early as 48 hours post fertilization and have been examined in numerous infection and immunotoxicology reports. One significant advantage of the zebrafish model is the ability to affordably generate high numbers of individual larvae that can be arrayed in multi-well plates for high throughput genetic and chemical exposure screens. However, traditional workflows for imaging individual larvae have been limited to low-throughput studies using traditional microscopes and manual analyses. Using a newly developed, parallelized microscope, the Multi-Camera Array Microscope (MCAM™), we have optimized a rapid, high-resolution algorithmic method to count fluorescently labeled cells in zebrafish larvae in vivo. Using transgenic zebrafish larvae, in which neutrophils express EGFP, we captured 18 gigapixels of images across a full 96-well plate, in 75 seconds, and processed the resulting datastream, counting individual fluorescent neutrophils in all individual larvae in 5 minutes. This automation is facilitated by a machine learning segmentation algorithm that defines the most in-focus view of each larva in each well after which pixel intensity thresholding and blob detection are employed to locate and count fluorescent cells. We validated this method by comparing algorithmic neutrophil counts to manual counts in larvae subjected to changes in neutrophil numbers, demonstrating the utility of this approach for high-throughput genetic and chemical screens where a change in neutrophil number is an endpoint metric. Using the MCAM™ we have been able to, within minutes, acquire both enough data to create an automated algorithm and execute a biological experiment with statistical significance. Finally, we present this open-source software package which allows the user to train and evaluate a custom machine learning segmentation model and use it to localize zebrafish and analyze cell counts within the segmented region of interest. This software can be modified as needed for studies involving other zebrafish cell lineages using different transgenic reporter lines and can also be adapted for studies using other amenable model species.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution.

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

Automated, high-throughput quantification of EGFP-expressing neutrophils in zebrafish by machine learning and a highly-parallelized microscope.

Normal development of the immune system is essential for overall health and disease resistance. Bony fish, such as the zebrafish (Danio rerio), possess all the major immune cell lineages as mammals and can be employed to model human host response to immune challenge. Zebrafish neutrophils, for example, are present in the transparent larvae as early as 48 hours post fertilization and have been examined in numerous infection and immunotoxicology reports. One significant advantage of the zebrafish model is the ability to affordably generate high numbers of individual larvae that can be arrayed in multi-well plates for high throughput genetic and chemical exposure screens. However, traditional workflows for imaging individual larvae have been limited to low-throughput studies using traditional microscopes and manual analyses. Using a newly developed, parallelized microscope, the Multi-Camera Array Microscope (MCAM™), we have optimized a rapid, high-resolution algorithmic method to count fluorescently labeled cells in zebrafish larvae in vivo. Using transgenic zebrafish larvae, in which neutrophils express EGFP, we captured 18 gigapixels of images across a full 96-well plate, in 75 seconds, and processed the resulting datastream, counting individual fluorescent neutrophils in all individual larvae in 5 minutes. This automation is facilitated by a machine learning segmentation algorithm that defines the most in-focus view of each larva in each well after which pixel intensity thresholding and blob detection are employed to locate and count fluorescent cells. We validated this method by comparing algorithmic neutrophil counts to manual counts in larvae subjected to changes in neutrophil numbers, demonstrating the utility of this approach for high-throughput genetic and chemical screens where a change in neutrophil number is an endpoint metric. Using the MCAM™ we have been able to, within minutes, acquire both enough data to create an automated algorithm and execute a biological experiment with statistical significance. Finally, we present this open-source software package which allows the user to train and evaluate a custom machine learning segmentation model and use it to localize zebrafish and analyze cell counts within the segmented region of interest. This software can be modified as needed for studies involving other zebrafish cell lineages using different transgenic reporter lines and can also be adapted for studies using other amenable model species.

Systemic Inflammation and Normocytic Anemia in DOCK11 Deficiency.

BACKGROUND: Increasing evidence links genetic defects affecting actin-regulatory proteins to diseases with severe autoimmunity and autoinflammation, yet the underlying molecular mechanisms are poorly understood. Dedicator of cytokinesis 11 (DOCK11) activates the small Rho guanosine triphosphatase (GTPase) cell division cycle 42 (CDC42), a central regulator of actin cytoskeleton dynamics. The role of DOCK11 in human immune-cell function and disease remains unknown. METHODS: We conducted genetic, immunologic, and molecular assays in four patients from four unrelated families who presented with infections, early-onset severe immune dysregulation, normocytic anemia of variable severity associated with anisopoikilocytosis, and developmental delay. Functional assays were performed in patient-derived cells, as well as in mouse and zebrafish models. RESULTS: We identified rare, X-linked germline mutations in DOCK11 in the patients, leading to a loss of protein expression in two patients and impaired CDC42 activation in all four patients. Patient-derived T cells did not form filopodia and showed abnormal migration. In addition, the patient-derived T cells, as well as the T cells from Dock11-knockout mice, showed overt activation and production of proinflammatory cytokines that were associated with an increased degree of nuclear translocation of nuclear factor of activated T cell 1 (NFATc1). Anemia and aberrant erythrocyte morphologic features were recapitulated in a newly generated dock11-knockout zebrafish model, and anemia was amenable to rescue on ectopic expression of constitutively active CDC42. CONCLUSIONS: Germline hemizygous loss-of-function mutations affecting the actin regulator DOCK11 were shown to cause a previously unknown inborn error of hematopoiesis and immunity characterized by severe immune dysregulation and systemic inflammation, recurrent infections, and anemia. (Funded by the European Research Council and others.).

A chromosome-level genome assembly of longnose gar, Lepisosteus osseus.

Holosteans (gars and bowfins) represent the sister lineage to teleost fishes, the latter being a clade that comprises over half of all living vertebrates and includes important models for comparative genomics and human health. A major distinction between the evolutionary history of teleosts and holosteans is that all teleosts experienced a genome duplication event in their early evolutionary history. As the teleost genome duplication occurred after teleosts diverged from holosteans, holosteans have been heralded as a means to bridge teleost models to other vertebrate genomes. However, only three species of holosteans have been genome-sequenced to date, and sequencing of more species is needed to fill sequence sampling gaps and provide a broader comparative basis for understanding holostean genome evolution. Here we report the first high quality reference genome assembly and annotation of the longnose gar (Lepisosteus osseus). Our final assembly consists of 22,709 scaffolds with a total length of 945 bp with contig N50 of 116.61 kb. Using BRAKER2, we annotated a total of 30,068 genes. Analysis of the repetitive regions of the genome reveals the genome to contain 29.12% transposable elements, and the longnose gar to be the only other known vertebrate outside of the spotted gar and bowfin to contain CR1, L2, Rex1, and Babar. These results highlight the potential utility of holostean genomes for understanding the evolution of vertebrate repetitive elements, and provide a critical reference for comparative genomic studies utilizing ray-finned fish models.

Legacy and emerging per- and polyfluoroalkyl substances suppress the neutrophil respiratory burst.

Per- and polyfluoroalkyl substances (PFASs) are used in a multitude of processes and products, including nonstick coatings, food wrappers, and fire-fighting foams. These chemicals are environmentally-persistent, ubiquitous, and can be detected in the serum of 98% of Americans. Despite evidence that PFASs alter adaptive immunity, few studies have investigated their effects on innate immunity. The report here presents results of studies that investigated the impact of nine environmentally-relevant PFASs [e.g. perfluorooctanoic acid (PFOA), perfluorooctane sulfonic acid potassium salt (PFOS-K), perfluorononanoic acid (PFNA), perfluorohexanoic acid (PFHxA), perfluorohexane sulfonic acid (PFHxS), perfluorobutane sulfonic acid (PFBS), ammonium perfluoro(2-methyl-3-oxahexanoate) (GenX), 7H-perfluoro-4-methyl-3,6-dioxa-octane sulfonic acid (Nafion byproduct 2), and perfluoromethoxyacetic acid sodium salt (PFMOAA-Na)] on one component of the innate immune response, the neutrophil respiratory burst. The respiratory burst is a key innate immune process by which microbicidal reactive oxygen species (ROS) are rapidly induced by neutrophils in response to pathogens; defects in the respiratory burst can increase susceptibility to infection. The study here utilized larval zebrafish, a human neutrophil-like cell line, and primary human neutrophils to ascertain whether PFAS exposure inhibits ROS production in the respiratory burst. It was observed that exposure to PFHxA and GenX suppresses the respiratory burst in zebrafish larvae and a human neutrophil-like cell line. GenX also suppressed the respiratory burst in primary human neutrophils. This report is the first to demonstrate that these PFASs suppress neutrophil function and support the utility of employing zebrafish larvae and a human cell line as screening tools to identify chemicals that may suppress human immune function.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors.

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Placing human gene families into their evolutionary context.

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Knockdown of Transmembrane Protein 150A (TMEM150A) Results in Increased Production of Multiple Cytokines.

Lipopolysaccharide (LPS)-induced signaling through Toll-like receptor 4 (TLR4) is mediated by the plasma membrane lipid, phosphatidylinositol (4,5)-bisphosphate [PI(4,5)P2] and its derivatives diacylglycerol and inositol trisphosphate. Levels of PI(4,5)P2 are controlled enzymatically and fluctuate in LPS-stimulated cells. Recently, transmembrane protein 150A (TMEM150A/TM6P1/damage-regulated autophagy modulator 5) has been shown to regulate PI(4,5)P2 production at the plasma membrane by modifying the composition of the phosphatidylinositol 4-kinase enzyme complex. To determine if TMEM150A function impacts TLR4 signaling, TMEM150A was knocked down in TLR4-expressing epithelial cells and cytokine expression quantified after LPS stimulation. In general, decreased expression of TMEM150A led to increased levels of LPS-induced cytokine secretion and transcript levels. Unexpectedly, knockdown of TMEM150A in a lung epithelial cell line (H292) also led to increased cytokine levels in the unstimulated conditions suggesting TMEM150A plays an important role in cellular homeostasis. Future studies will investigate if TMEM150A plays a similar role for other TLR agonists and in other cell lineages.

Therapeutic targeting of LCK tyrosine kinase and mTOR signaling in T-cell acute lymphoblastic leukemia.

Relapse and refractory T-cell acute lymphoblastic leukemia (T-ALL) has a poor prognosis, and new combination therapies are sorely needed. Here, we used an ex vivo high-throughput screening platform to identify drug combinations that kill zebrafish T-ALL and then validated top drug combinations for preclinical efficacy in human disease. This work uncovered potent drug synergies between AKT/mTORC1 (mammalian target of rapamycin complex 1) inhibitors and the general tyrosine kinase inhibitor dasatinib. Importantly, these same drug combinations effectively killed a subset of relapse and dexamethasone-resistant zebrafish T-ALL. Clinical trials are currently underway using the combination of mTORC1 inhibitor temsirolimus and dasatinib in other pediatric cancer indications, leading us to prioritize this therapy for preclinical testing. This combination effectively curbed T-ALL growth in human cell lines and primary human T-ALL and was well tolerated and effective in suppressing leukemia growth in patient-derived xenografts (PDX) grown in mice. Mechanistically, dasatinib inhibited phosphorylation and activation of the lymphocyte-specific protein tyrosine kinase (LCK) to blunt the T-cell receptor (TCR) signaling pathway, and when complexed with mTORC1 inhibition, induced potent T-ALL cell killing through reducing MCL-1 protein expression. In total, our work uncovered unexpected roles for the LCK kinase and its regulation of downstream TCR signaling in suppressing apoptosis and driving continued leukemia growth. Analysis of a wide array of primary human T-ALLs and PDXs grown in mice suggest that combination of temsirolimus and dasatinib treatment will be efficacious for a large fraction of human T-ALLs.

On the relationship between extant innate immune receptors and the evolutionary origins of jawed vertebrate adaptive immunity.

For over half a century, deciphering the origins of the genomic loci that form the jawed vertebrate adaptive immune response has been a major topic in comparative immunogenetics. Vertebrate adaptive immunity relies on an extensive and highly diverse repertoire of tandem arrays of variable (V), diversity (D), and joining (J) gene segments that recombine to produce different immunoglobulin (Ig) and T cell receptor (TCR) genes. The current consensus is that a recombination-activating gene (RAG)-like transposon invaded an exon of an ancient innate immune VJ-bearing receptor, giving rise to the extant diversity of Ig and TCR loci across jawed vertebrates. However, a model for the evolutionary relationships between extant non-recombining innate immune receptors and the V(D)J receptors of the jawed vertebrate adaptive immune system has only recently begun to come into focus. In this review, we provide an overview of non-recombining VJ genes, including CD8ß, CD79b, natural cytotoxicity receptor 3 (NCR3/NKp30), putative remnants of an antigen receptor precursor (PRARPs), and the multigene family of signal-regulatory proteins (SIRPs), that play a wide range of roles in immune function. We then focus in detail on the VJ-containing novel immune-type receptors (NITRs) from ray-finned fishes, as recent work has indicated that these genes are at least 50 million years older than originally thought. We conclude by providing a conceptual model of the evolutionary origins and phylogenetic distribution of known VJ-containing innate immune receptors, highlighting opportunities for future comparative research that are empowered by this emerging evolutionary perspective.

Transcriptome annotation reveals minimal immunogenetic diversity among Wyoming toads, Anaxyrus baxteri.

Briefly considered extinct in the wild, the future of the Wyoming toad (Anaxyrus baxteri) continues to rely on captive breeding to supplement the wild population. Given its small natural geographic range and history of rapid population decline at least partly due to fungal disease, investigation of the diversity of key receptor families involved in the host immune response represents an important conservation need. Population decline may have reduced immunogenetic diversity sufficiently to increase the vulnerability of the species to infectious diseases. Here we use comparative transcriptomics to examine the diversity of toll-like receptors and major histocompatibility complex (MHC) sequences across three individual Wyoming toads. We find reduced diversity at MHC genes compared to bufonid species with a similar history of bottleneck events. Our data provide a foundation for future studies that seek to evaluate the genetic diversity of Wyoming toads, identify biomarkers for infectious disease outcomes, and guide breeding strategies to increase genomic variability and wild release successes.

A Zebrafish Model of Metastatic Colonization Pinpoints Cellular Mechanisms of Circulating Tumor Cell Extravasation.

Metastasis is a multistep process in which cells must detach, migrate/invade local structures, intravasate, circulate, extravasate, and colonize. A full understanding of the complexity of this process has been limited by the lack of ability to study these steps in isolation with detailed molecular analyses. Leveraging a comparative oncology approach, we injected canine osteosarcoma cells into the circulation of transgenic zebrafish with fluorescent blood vessels in a biologically dynamic metastasis extravasation model. Circulating tumor cell clusters that successfully extravasated the vasculature as multicellular units were isolated under intravital imaging (n = 6). These extravasation-positive tumor cell clusters sublines were then molecularly profiled by RNA-Seq. Using a systems-level analysis, we pinpointed the downregulation of KRAS signaling, immune pathways, and extracellular matrix (ECM) organization as enriched in extravasated cells (p < 0.05). Within the extracellular matrix remodeling pathway, we identified versican (VCAN) as consistently upregulated and central to the ECM gene regulatory network (p < 0.05). Versican expression is prognostic for a poorer metastasis-free and overall survival in patients with osteosarcoma. Together, our results provide a novel experimental framework to study discrete steps in the metastatic process. Using this system, we identify the versican/ECM network dysregulation as a potential contributor to osteosarcoma circulating tumor cell metastasis.

Single-minded 2 is required for left-right asymmetric stomach morphogenesis.

The morphogenesis of left-right (LR) asymmetry is a crucial phase of organogenesis. In the digestive tract, the development of anatomical asymmetry is first evident in the leftward curvature of the stomach. To elucidate the molecular events that shape this archetypal laterality, we performed transcriptome analyses of the left versus right sides of the developing stomach in frog embryos. Besides the known LR gene pitx2, the only gene found to be expressed asymmetrically throughout all stages of curvature was single-minded 2 (sim2), a Down Syndrome-related transcription factor and homolog of a Drosophila gene (sim) required for LR asymmetric looping of the fly gut. We demonstrate that sim2 functions downstream of LR patterning cues to regulate key cellular properties and behaviors in the left stomach epithelium that drive asymmetric curvature. Our results reveal unexpected convergent cooption of single-minded genes during the evolution of LR asymmetric morphogenesis, and have implications for dose-dependent roles of laterality factors in non-laterality-related birth defects.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system.

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

The bowfin genome illuminates the developmental evolution of ray-finned fishes.

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

In vivo assessment of respiratory burst inhibition by xenobiotic exposure using larval zebrafish.

Currently, assessment of the potential immunotoxicity of a given agent involves a tiered approach for hazard identification and mechanistic studies, including observational studies, evaluation of immune function, and measurement of susceptibility to infectious and neoplastic diseases. These studies generally use costly low-throughput mammalian models. Zebrafish, however, offer an excellent alternative due to their rapid development, ease of maintenance, and homology to mammalian immune system function and development. Larval zebrafish also are a convenient model to study the innate immune system with no interference from the adaptive immune system. In this study, a respiratory burst assay (RBA) was utilized to measure reactive oxygen species (ROS) production after developmental xenobiotic exposure. Embryos were exposed to non-teratogenic doses of chemicals and at 96 h post-fertilization, the ability to produce ROS was measured. Using the RBA, 12 compounds with varying immune-suppressive properties were screened. Seven compounds neither suppressed nor enhanced the respiratory burst; five reproducibly suppressed global ROS production, but with varying potencies: benzo[a]pyrene, 17ß-estradiol, lead acetate, methoxychlor, and phenanthrene. These five compounds have all previously been reported as immunosuppressive in mammalian innate immunity assays. To evaluate whether the suppression of ROS by these compounds was a result of decreased immune cell numbers, flow cytometry with transgenic zebrafish larvae was used to count the numbers of neutrophils and macrophages after chemical exposure. With this assay, benzo[a]pyrene was found to be the only chemical that induced a change in the number of immune cells by increasing macrophage but not neutrophil numbers. Taken together, this work demonstrates the utility of zebrafish larvae as a vertebrate model for identifying compounds that impact innate immune function at non-teratogenic levels and validates measuring ROS production and phagocyte numbers as metrics for monitoring how xenobiotic exposure alters the innate immune system.