Intestinal Serum amyloid A suppresses systemic neutrophil activation and bactericidal activity in response to microbiota colonization.

The intestinal microbiota influences the development and function of myeloid lineages such as neutrophils, but the underlying molecular mechanisms are unresolved. Using gnotobiotic zebrafish, we identified the immune effector Serum amyloid A (Saa) as one of the most highly induced transcripts in digestive tissues following microbiota colonization. Saa is a conserved secreted protein produced in the intestine and liver with described effects on neutrophils in vitro, however its in vivo functions remain poorly defined. We engineered saa mutant zebrafish to test requirements for Saa on innate immunity in vivo. Zebrafish mutant for saa displayed impaired neutrophil responses to wounding but augmented clearance of pathogenic bacteria. At baseline, saa mutants exhibited moderate neutrophilia and altered neutrophil tissue distribution. Molecular and functional analyses of isolated neutrophils revealed that Saa suppresses expression of pro-inflammatory markers and bactericidal activity. Saa's effects on neutrophils depended on microbiota colonization, suggesting this protein mediates the microbiota's effects on host innate immunity. To test tissue-specific roles of Saa on neutrophil function, we over-expressed saa in the intestine or liver and found that sufficient to partially complement neutrophil phenotypes observed in saa mutants. These results indicate Saa produced by the intestine in response to microbiota serves as a systemic signal to neutrophils to restrict aberrant activation, decreasing inflammatory tone and bacterial killing potential while simultaneously enhancing their ability to migrate to wounds.

Introduction

The vertebrate intestine is densely colonized with complex communities of micro-organisms, collectively referred to as the intestinal microbiota. Studies using gnotobiotic animals have demonstrated that microbiota colonization is required for the normal development of an innate immune system capable of mounting appropriate responses to diverse challenges such as infection and injury [1]. Despite being spatially confined to the intestinal lumen by physical and chemical barriers such as the intestinal epithelium and mucus, the microbiota influences both local and systemic host immune development and function [2, 3]. However, the specific molecular mechanisms by which the microbiota impacts local and systemic host immune responses remain poorly defined.

Intestinal epithelial cells (IECs) serve as the primary host interface with the intestinal microbiota and secrete a myriad of factors following microbial colonization [4]. We and others have hypothesized that these microbiota-induced IEC products may mediate the microbiota’s influences on the host immune system [5, 6]. Previous studies have identified a secreted host factor, Serum amyloid A (Saa), that is potently up-regulated in the intestine following microbial colonization in zebrafish and mice [7-11]. SAA is highly conserved amongst vertebrates, existing as a single gene in fishes and birds and as a multi-gene family in mammals [12, 13]. While basal SAA production is stimulated by the microbiota, SAA production is also markedly augmented following acute injury and infection as a part of the acute phase response, whereby circulating levels can reach 1 mg/mL [13-15]. Moreover, SAA is elevated in chronic pathological conditions where both local and circulating concentrations are positively correlated with inflammation. Accordingly, SAA is an established biomarker for chronic inflammatory diseases such as diabetes, atherosclerosis, and inflammatory bowel disease (IBD) [16-20]. Taken together, SAA’s high degree of evolutionary conservation coupled with its strong induction following inflammatory stimuli suggests important roles for SAA in health and disease.

Previous studies have reported both pro- and anti-inflammatory effects of SAA on host immune responses. The vast majority of these studies have been conducted in vitro using recombinant human SAA (rhSAA), which has been shown to directly influence granulocytes, including monocytes and neutrophils, promoting the production of inflammatory cytokines, reactive oxygen species (ROS), and directing motility [21-25]. Recent reports have also shown that mammalian SAAs can bind retinol and mediate host responses during infection [8]. Moreover, induction of SAA in the intestine following colonization with specific microorganisms such as segmented filamentous bacteria (SFB) can shape local adaptive immune cell development by promoting Th17 differentiation [26, 27]. However, a fuller assessment of SAA’s in vivo functional roles has remained elusive due to the existence of multiple SAA gene paralogs in mammals (3 in humans, 4 in mice) and the use of cell-culture based assays performed with rhSAA that behaves dissimilarly to endogenous SAA protein [28, 29].

The existence of a single SAA ortholog in fishes provides interesting opportunities to define SAA’s in vivo functional roles. We previously demonstrated that partial (~30%) knockdown of saa transcript in zebrafish influenced tissue-specific neutrophil behaviors in vivo, leading us to hypothesize that Saa regulates neutrophil activity in homeostasis [12]. However, Saa’s influence on systemic neutrophil function in homeostasis and in the relevant contexts of injury and infection remained unresolved. Neutrophils are professional phagocytic myeloid cells that play critical roles in host defense against pathogens. The most abundant immune cell in circulation and the first to be recruited to sites of injury, neutrophils eliminate microbial invaders and debris through a variety of mechanisms including phagocytosis, generation of ROS, and secretion of anti-microbial peptides [30, 31]. Neutrophils are conditioned by host and microbially derived signals, including pathogen associated molecular patterns (PAMPs) and damage associated molecular patterns (DAMPs), allowing for proper responses to inflammatory stimuli [32, 33]. Previous studies have shown that microbial colonization of the intestine promotes neutrophil differentiation, activation, and response to peripheral injury [10, 12, 34–39]. However the signaling molecules that mediate these interactions in vivo remain largely unknown.

Here, using zebrafish, we demonstrate that Saa is a host factor that signals microbial status in the intestine to extra-intestinal populations of immune cells and directs their responses to inflammatory stimuli. Zebrafish share highly conserved hematopoietic programs with other vertebrates, including specification of myeloid lineages, which can be coupled with optical transparency and genetic tractability to allow for high resolution in vivo imaging of innate immune processes [40, 41]. Moreover, the zebrafish genome encodes a single saa ortholog, allowing us to generate the first-ever saa null vertebrate model. By comparing wild-type (WT) zebrafish to those that lack saa or express it only in intestinal epithelial cells or hepatocytes, we reveal Saa’s impact on systemic neutrophil activity in homeostasis and following bacterial infection and wounding. Using tissue specific rescue, we demonstrate that liver- and intestinally-derived Saa can shape systemic neutrophil function and intestinal Saa can even restore neutrophil defects observed in germ-free zebrafish.

Results

Saa promotes neutrophil migration during injury and homeostasis

To investigate Saa’s effects on neutrophil function in vivo, we first generated saa mutant zebrafish, identifying three independent deletion alleles all resulting in frameshift mutations within saa exon 2 (S1A–S1C Fig). The largest saa deletion allele (22 bp, designated rdu60, homozygous mutants hereafter referred to as saa) resulted in 90% reduced saa mRNA (S1D–S1F Fig). saa zebrafish survived to adulthood and exhibited no gross developmental abnormalities (S1H–S1L Fig). We performed caudal fin amputations on WT and saa Tg(lyz:EGFP) larvae and quantified neutrophil recruitment to the wound margin, observing fewer neutrophils at 6 hours post-wounding in saa larvae (Fig 1A and 1B, S1N Fig). Further, in vivo imaging revealed that neutrophils in the vicinity of the wound moved with reduced mean speed in saa mutant larvae (Fig 1C). Importantly, saa mRNA was not upregulated at 6 hours post amputation in WT larvae (S1M Fig), demonstrating acute saa induction did not affect neutrophil activity. Thus, Saa is required for neutrophil mobilization to sites of injury independent of systemic induction.

Fig 1

Saa mediates neutrophil behavior in vivo.

(A-B) Imaging and quantification of lyz:EGFP+ neutrophils recruited to tail wound margin over 6 hours following caudal fin amputation (dashed red line indicates wound margin; scale bar = 250 μm) (n ≥ 24 larvae / genotype at 6 hour time point). (C) Measurement of lyz:DsRed+ neutrophil speed from time-lapse imaging in caudal fin tissue over 6 hour period following amputation (n = 4 larvae / genotype, 87–112 cells tracked / genotype). (D-E) Representative images of 6 dpf Tg(lyz:DsRed) WT and saa-/- larvae (scale bar = 500 μm). (F) Enumeration of intestine-associated lyz:EGFP+ cells in 6 dpf larvae (n = 32–40 larvae / genotype). (G-J) Quantitative analysis of lyz:EGFP+ neutrophil behavior from time-lapse imaging of distinct anatomical compartments (intestine and trunk, ROIs in panel D) in 6 dpf larval zebrafish (6 larvae / genotype, ≥ 23 cells analyzed / genotype / tissue). Data analyzed by t-test. For panel B, statistical comparisons were performed within each time point. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

Given that Saa affects neutrophil responses to injury, we tested if Saa regulates basal neutrophil behavior. Analysis of neutrophils in homeostasis revealed Saa promoted neutrophil speed in the trunk and linear migration in the intestine (Fig 1G–1J), consistent with our prior saa morpholino knockdown data [12]. Considering saa is expressed in IECs [10], we reasoned that Saa may promote neutrophil recruitment to the intestine. Indeed, we observed fewer intestine-associated neutrophils in saa larvae (Fig 1D–1F). These data demonstrate that Saa promotes neutrophil recruitment to distinct tissues such as the intestine during homeostasis as well as to peripheral tissues following injury.

Saa restricts systemic neutrophil abundance and bactericidal activity

Given that saa loss is associated with impaired neutrophil recruitment to wounds and healthy tissues, we asked whether systemic neutrophil abundance is altered in saa mutants. We enumerated systemic neutrophils by flow cytometry (Fig 2A, S2 Fig) and consistently observed elevated neutrophil abundance in saa larvae. This was corroborated by increased expression of the granulocyte marker genes lysozyme C (lyz), l-plastin (lcp), and the granulopoietic cytokine colony stimulating factor 3a (csf3a, also known as gcsf1a) in 6 days post fertilization (dpf) whole larvae (Fig 2D–2F). Morphological classification of lyz+ neutrophils into sub-populations from cytospin preparations (adapted from [42]) revealed an over-representation of immature lyz+ neutrophils in saa deficient animals compared to WT controls (Fig 2B and 2C). Together, these data demonstrate a novel role for Saa regulating neutrophil maturation in vivo.

Fig 2

Saa regulates neutrophil abundance and maturation.

(A) Flow cytometry analysis of lyz:EGFP+ neutrophil abundance from whole 6 dpf WT and saa mutant zebrafish larvae (results are combined from 3 independent experiments, ≥ 4 replicates / genotype / experiment, 60–90 larvae / replicate). (B-C) Morphological analysis of lyz:DsRed+ neutrophil cytospins stained with Wright-Giemsa from adult WT and saa mutant zebrafish kidneys (5–6 adult zebrafish kidneys / genotype / experiment, 2 independent experiments, n ≥ 199 cells analyzed / genotype / experiment) (scale bar = 10 μm). (D-F) qRT-PCR of leukocyte-associated transcripts lysozyme C (lyz), l-plastin (lcp), and colony-stimulating factor 3a (csf3a) from 6 dpf whole zebrafish larvae (n = 4 replicates / genotype / timepoint, 25–30 larvae / replicate). Data in panel A analyzed by t-test. Data in panel C analyzed by chi-squared test. Data in panels D-F analyzed by one-way ANOVA with Tukey’s multiple comparisons test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

Exposure to PAMPs or inflammatory host molecules can elicit defined neutrophil transcriptional responses, reflecting their activation state [43-47]. Gene expression analysis of FACS-isolated neutrophils revealed elevated expression of genes encoding pro-inflammatory cytokines (tnfa, il1b), antimicrobial peptides (pglyrp2, pglyrp5), and regulators of ROS production (mpx, ncf1) in saa larvae (Fig 3A). These transcriptional differences suggest Saa restricts basal neutrophil activation. As neutrophils primarily function to clear microbial infections [48, 49], we co-cultured neutrophils isolated from adult zebrafish kidneys with non-pathogenic Escherichia coli then assessed bacterial viability. Isolated neutrophils from both WT and saa fish were viable ex vivo and exhibited morphological responses to bacteria (e.g., extending cytosolic projections; Fig 3B, S3A–S3E and S3J Fig). Co-culture with E. coli induced il1b mRNA in WT neutrophils, demonstrating zebrafish neutrophils respond transcriptionally to bacteria ex vivo (Fig 3C). Moreover, il1b transcript levels were significantly increased in unstimulated saa mutant adult kidney neutrophils compared to WT, consistent with our observations from larval neutrophils (Fig 3A and 3C). However, following co-culture with E. coli, il1b in WT neutrophils reached similar levels measured in saa neutrophils (Fig 3C).

Fig 3

Saa suppresses neutrophil transcriptional activation and bactericidal activity.

(A) qRT-PCR of pro-inflammatory mRNAs from sorted neutrophils from 6 dpf WT and saa mutant zebrafish larvae (5,000–12,000 lyz:EGFP+ cells / replicate, 3–6 replicates / genotype / experiment, 60–90 larvae / replicate). (B) Microscopic analysis revealed neutrophils isolated from adult zebrafish kidneys extend protrusions in response to bacterial signals ex vivo (white arrows; scale bar = 20 μm). (C) il1b expression in un-stimulated and E. coli exposed lyz:EGFP+ neutrophils from WT and saa mutant zebrafish following 4 hours ex vivo culture (3–5 replicates / genotype / condition). (D) CFU quantification of bacterial concentration following 4 hour co-culture of isolated lyz:EGFP+ neutrophils from WT and saa mutant zebrafish with E. coli (MOI 2). (E) Quantification of intracellular ROS levels by CellROX fluorescence from neutrophils cultured ex vivo with and without E. coli (lyz:EGFP+ cells isolated from 6 zebrafish adult kidneys / genotype, quantification of ≥ 177 cells / condition). In panels C and E, a one-way ANOVA with Tukey’s multiple comparisons test was used. Data in panels A and D were analyzed with a t-test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

To assess neutrophil bactericidal activity, we enumerated CFUs following 4 hours of co-culture and found saa neutrophils killed significantly more bacteria than WT neutrophils (Fig 3D). To interrogate possible mechanisms of bacterial clearance, we labeled neutrophils with CellROX ex vivo and measured levels of intracellular ROS by confocal microscopy (S3F and S3G Fig). Bacterial exposure resulted in decreased ROS in both WT and saa neutrophils, indicating bacteria stimulated neutrophil degranulation (Fig 3E). Interestingly, neutrophils from saa mutant animals had elevated levels of ROS relative to WT controls both at baseline and after bacterial stimulation (Fig 3E), which is consistent with their augmented bacterial killing activity (Fig 3D). Collectively, these data indicate neutrophils from saa mutants are aberrantly activated as evidenced by elevated pro-inflammatory mRNA expression, augmented bactericidal activity, and elevated ROS production, and suggest Saa restricts systemic neutrophil inflammatory tone in vivo.

Intestinally-derived Saa alters neutrophil distribution and activity

Since saa is highly expressed in the larval zebrafish intestine and liver following microbiota colonization [10], we reasoned that Saa produced by specific tissues might differentially affect systemic neutrophil conditioning. Saa produced by the liver enters circulation and has systemic affects [50], however the distal influences of intestinally-derived saa remained unresolved. We generated a transgenic line to drive expression of zebrafish saa in hepatocytes using a 2.8 kb promoter fragment from the zebrafish fabp10a gene [Tg(-2.8fabp10a:saa);cmlc2:EGFP, hereafter referred to as Tg(fabp10a:saa)] [51-54]. We crossed this transgene into the saa mutant background and assessed whether hepatic saa was sufficient to rescue neutrophil defects observed in saa mutant animals.

Despite elevated levels of saa transcript in saa;Tg(fabp10a:saa) relative to saa mutant larvae (Fig 4A), we still observed reduced intestinal neutrophil recruitment compared to WT controls (Fig 4B). To test the impact of liver-derived saa on neutrophil function, we quantified neutrophil recruitment to the wound margin following caudal fin amputation. Neutrophil recruitment at 6 hours post wounding was restored to WT levels in saa;Tg(fabp10a:saa) larvae, demonstrating that hepatic saa is sufficient to restore neutrophil response to peripheral injury (Fig 4C). To further investigate the effects of hepatic-derived Saa on neutrophil function, we quantified bactericidal activity of neutrophils from WT, saa-/-, and saa;Tg(fabp10a:saa) zebrafish. Whereas E. coli is cleared quickly by the zebrafish immune system, Pseudomonas aeruginosa is capable of establishing systemic infections [55, 56]. Interestingly, transgenic hepatic saa was sufficient to restore P. aeruginosa killing activity of isolated neutrophils to WT levels in ex vivo co-culture (Fig 4D). These data establish that saa produced in the liver is capable of influencing a distinct subset of systemic neutrophil functions, promoting their ability to migrate to wounds yet restricting bacterial killing.

Fig 4

Intestinally-derived Saa regulates systemic neutrophil activity.

(A) qRT-PCR of saa from whole 6 dpf larvae of the indicated genotypes (n = 4 replicates / genotype, 25–30 larvae / replicate) (B) Enumeration of intestine-associated lyz:DsRed+ neutrophils along the anterior to posterior axis (segment 1 to segment 3) in 6 dpf larvae (n ≥ 25 larvae / genotype). (C) lyz:DsRed+ neutrophil recruitment to caudal fin wound 6 hours following amputation in 6 dpf zebrafish larvae (n ≥ 25 larvae / genotype at 6 hour time point). (D) CFU quantification of bacterial concentration following 4 hour co-culture of lyz:DsRed+ adult zebrafish neutrophils with P. aeruginosa (P.a., MOI 0.2) (8 replicates / genotype). (E) Representative stereoscope images of IEC specific mCherry expression in 6 dpf Tg(-0.349cldn15la:mCherry) larvae compared to non-transgenic (NTG) controls. White dashed line indicates the intestine (scale bar = 500 μm). (F) Representative confocal micrograph of immunostained transverse section of Tg(-0.349cldn15la:mCherry) 6 dpf larvae labeled with the absorptive cell brush border-specific antibody 4E8 (scale bar = 20 μm). (G,H) qRT-PCR of cldn15la, fabp2, and fabp10a (G) from sorted Tg(-0.349cldn15la:mCherry)+ IECs and saa (H) from cldn15la:mCherry and negative cells isolated from 6 dpf larvae of indicated genotypes (13,000 –0.349cldn15la:mCherry+ or mCherry negative cells / replicate, 4 replicates / genotype, 30 larvae / replicate). (I) qRT-PCR of saa from 6 dpf larval dissected digestive tissue of the indicated genotypes (n = 4 replicates / genotype, 25–30 larvae / replicate) (J) Enumeration of intestine-associated lyz:DsRed+ neutrophils in 6 dpf larvae (n = 30 larvae / genotype). (K) lyz:DsRed+ neutrophil recruitment to caudal fin wound 6 hours following amputation in 6 dpf zebrafish larvae (n ≥ 18 larvae / genotype at 6 hour time point). (L) CFU quantification of bacterial concentration following 4 hour co-culture of lyz:DsRed+ adult zebrafish neutrophils with P. aeruginosa (P.a., MOI 0.2) (3–6 replicates / genotype). (M) il1b qRT-PCR from lyz:DsRed+ neutrophils co-cultured with and without P.a. ex vivo for 4 hours (n ≥ 2 replicates / condition). (N) CFU quantification of in vivo P.a. bacterial burden following systemic infection of larval zebrafish at 5 days post infection (dpi) (data from 3 independent experiments, n ≥ 30 larvae / genotype). Data in panels A-D and H-N were analyzed by one-way ANOVA with Tukey’s multiple comparisons test. A Mann-Whitney test was applied to panel G. For panels C and K, statistical comparisons were performed amongst samples within the same time point. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

Given that the intestine serves as a primary interface between the host and microbiota, and that saa is transcriptionally upregulated in the intestine following microbial colonization, we sought to determine if intestinally-derived saa could also influence systemic neutrophil phenotypes. To do so, we first needed to engineer transgenic zebrafish in which saa is expressed specifically in IECs. While several zebrafish promoters have been identified with intestine-restricted activity, they only drive transgene expression in subsets of IECs (e.g., the fabp2/ifabp promoter is active in enterocytes in the anterior intestine) [10, 57]. Since zebrafish cldn15la is expressed more broadly and specifically in IECs [58, 59], we queried adult zebrafish IEC FAIRE-seq data and identified a 349 bp open chromatin region in the cldn15la promoter which contains predicted binding sites for intestine-specific transcription factors including Cdx2 (S4A and S4B Fig) [60]. To confirm intestinal specificity of this cldn15la promoter fragment, we used this element to drive expression of mCherry [Tg(-0.349cldn15la:mCherry)], and observed mCherry fluorescence restricted to IECs (Fig 4E and 4F, S4C–S4F Fig, S5A–S5E Fig) [60, 61]. Gene expression analysis of enterocyte marker fabp2 and pan-IEC marker cldn15la in FACS-isolated mCherry+ and negative cells from Tg(-0.349cldn15la:mCherry) larvae revealed strong enrichment in the mCherry+ fraction, further validating the IEC specificity of this promoter (Fig 4G, S4G and S4H Fig).

We next generated Tg(-0.349cldn15la:saa;cmlc2:EGFP) zebrafish [subsequently denoted Tg(cldn15la:saa)], which express zebrafish saa in IECs (S6A and S6B Fig). To confirm that saa expression was restricted to IECs, we crossed Tg(-0.349cldn15la:mCherry) and Tg(cldn15la:saa) adults, then sorted mCherry+ IECs and measured levels of saa expression in mCherry+ and negative fractions from WT and Tg(cldn15la:saa) 6 dpf larvae (S4G and S4H Fig). We observed significantly higher saa levels in mCherry+ IECs isolated from Tg(cldn15la:saa) larvae compared to WT IECs and both negative fractions (Fig 4H). This data indicates that saa expression in Tg(cldn15la:saa) larvae is driven specifically in IECs.

We crossed this transgene into the saa mutant background, and asked whether IEC-derived Saa could complement saa neutrophil defects. As expected, we observed fewer intestine-associated neutrophils in each intestinal segment of saa larvae (from anterior to posterior, segments 1–3). Intestinally-derived Saa was sufficient to partially complement this mutant phenotype, with an increased abundance of intestine-associated neutrophils in all segments (Fig 4I and 4J). Thus Saa produced in IECs, unlike hepatic Saa, is sufficient to promote neutrophil recruitment to the intestine, suggesting Saa is a neutrophil chemoattractant in vivo. To determine if systemic neutrophil function was altered by intestinally-derived Saa, we performed caudal fin amputations. At 6 hours post amputation, neutrophil recruitment to the wound in saa;Tg(cldn15la:saa) larvae was equivalent to the WT response (Fig 4K), demonstrating intestinally-derived saa is sufficient to restore neutrophil mobilization to the caudal fin in otherwise saa deficient larvae.

To test the effects of intestinally-derived Saa on neutrophil function, we measured bactericidal activity of neutrophils from WT, saa-/-, and saa;Tg(cldn15la:saa) zebrafish. We observed reduced survival of both P. aeruginosa and E. coli following co-culture with saa neutrophils vs WT controls. However, neutrophils from Tg(cldn15la:saa)+ saa zebrafish exhibited comparable bactericidal activity to WT neutrophils (Fig 4L, S6D Fig). Additionally, following P. aeruginosa co-culture, neutrophils from WT, saa mutant, and saa;Tg(cldn15la:saa) zebrafish exhibited comparable il1b transcript induction (Fig 4M), demonstrating that neutrophils from each genotype are responsive to bacterial stimulation ex vivo. These findings demonstrate intestinally-derived Saa is sufficient to constrain bactericidal activity of neutrophils.

Considering the profound effects of saa levels and source on neutrophil antibacterial function ex vivo, we asked if intestinal Saa influenced bacterial clearance in vivo using systemic P. aeruginosa infection. As neutrophils play an important role in P. aeruginosa clearance in vivo [62, 63], we postulated that differences in bacterial clearance between WT and saa deficient larvae would be largely driven by differences in neutrophil activation and behavior. Larval zebrafish were injected with P. aeruginosa via the caudal vein to achieve systemic infection and bacterial burden was assessed from whole zebrafish larvae. Consistent with our ex vivo results, we observed enhanced bacterial clearance in saa larvae which is restored to WT levels in saa;Tg(cldn15la:saa) larvae (Fig 4N).

Flow cytometry revealed that transgenic intestinal Saa did not significantly alter systemic neutrophil abundance in saa larvae whereas liver derived Saa returned neutrophil abundance to WT levels (S7A and S7C Fig). Moreover, qRT-PCR analysis of neutrophils isolated from both saa;Tg(cldn15la:saa) and saa;Tg(fabp10a:saa) larvae revealed persistently elevated expression of pro-inflammatory and anti-microbial effectors (pglyrp2, tnfa, il1b, ncf1) (S7B and S7D Fig). Thus, liver and intestinal expression of saa was sufficient to rescue only a subset of neutrophil defects observed in saa mutant animals. These observations highlight the potential requirement for other tissue sources, presentation, or temporal control of Saa to condition different aspects of systemic neutrophil function. Collectively, these data demonstrate that expression of saa in the intestinal epithelium and liver is sufficient to promote neutrophil recruitment to a peripheral wound and to restrict bactericidal activity in vivo and ex vivo, but is unable to dampen elevated neutrophil pro-inflammatory mRNA profiles observed in saa mutant zebrafish. These results further resolve the local and systemic roles for intestinal Saa, promoting neutrophil recruitment to the intestine as well as shaping systemic neutrophil migratory and bacterial killing activity.

Microbiota-induced Saa suppresses neutrophil pro-inflammatory mRNA production and promotes migration to a tail wound

We and others have shown that the intestinal microbiota influences a variety of neutrophil phenotypes both in homeostasis and following injury [12, 36–38]. Since our results indicate that Saa suppresses neutrophil activation (Fig 3A, 3D and 3E) and previous studies have reported Saa has direct bactericidal activity [64-66] we speculated that saa deficiency may impact microbiota composition. We used 16S rRNA gene sequencing to compare bacterial communities in the digestive tracts of co-housed saa and WT sibling zebrafish at both larval and adult stages. We found that saa genotype had no significant effects on gut bacterial community composition as measured by alpha- or beta-diversity metrics at larval (6 dpf) or adult (70 dpf) stages (S2 and S3 Tables). These findings indicate that Saa does not broadly impact gut microbiota composition in zebrafish.

Since saa is potently induced following microbial colonization (Fig 5A), we asked if Saa regulates neutrophil activation in response to the microbiota. We isolated neutrophils from 6 dpf gnotobiotic WT and saa larvae, and found expression of pro-inflammatory mRNAs was significantly elevated in WT neutrophils from conventionalized (CV) zebrafish vs germ-free (GF) controls, confirming microbiota-derived signals induce neutrophil pro-inflammatory mRNA expression (Fig 5B). Comparison of the same transcripts in neutrophils from saa zebrafish reared under CV and GF conditions revealed augmented induction of mRNAs following microbiota colonization. Since neutrophil pro-inflammatory mRNA levels were comparable between WT GF and saa GF larvae, we conclude the microbiota potentiates transcriptional activation observed in saa larvae (Fig 5B). To assess the impact of the microbiota on neutrophil function in WT and saa mutant animals, we measured neutrophil recruitment to peripheral tail wound injury in GF larvae. Unlike our previous observations in conventionally reared larvae (Fig 1B), neutrophil recruitment to the fin wound was indistinguishable in GF saa and GF WT larvae at 6 hours post wounding (S8A Fig). This indicates that neutrophil recruitment defects only manifest in saa mutants when colonized with a microbiota. These data further demonstrate that Saa functions to restrict aberrant activation of neutrophils by the microbiota at homeostasis, thus allowing neutrophils to respond appropriately to injury.

Fig 5

Microbiota-induced Saa conditions neutrophils in vivo.

(A) qRT-PCR of saa from gnotobiotic zebrafish larvae following microbiota colonization (CV) at 3 dpf versus germ free (GF) (3 replicates / condition / timepoint, n ≥ 20 larvae / replicate). (B) qRT-PCR of neutrophils isolated from 6 dpf gnotobiotic WT and saa mutant zebrafish larvae (3 replicates / genotype / condition, n ≥ 27 larvae / replicate). (C) lyz:EGFP+ neutrophil recruitment to caudal fin wound margin 6 hours after amputation in 6 dpf gnotobiotic WT and Tg(cldn15la:saa) sibling zebrafish larvae (n ≥ 18 larvae / genotype / condition at the 6 hour timepoint). (D) Working model depicting signals from the microbiota evoking production of intestinal and hepatic Saa leading to shared and distinct effects on systemic neutrophil function. Statistical comparisons of data in panel in A were performed within each time point and analyzed by t-test. Data in panel B analyzed by two-way ANOVA with p values reported in the table. Data in panel C analyzed by one-way ANOVA with Tukey’s multiple comparisons test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

Given that the microbiota induces saa expression in the intestine, and GF WT zebrafish larvae have diminished neutrophil responses to wounding compared to conventionalized WT siblings [10], we asked if transgenic saa expression in IECs was sufficient to rescue neutrophil deficiencies in GF animals. Consistent with previous findings, we observed reduced neutrophil wound recruitment to caudal fin amputation in GF WT larvae vs CV controls [12]. Transgenic intestinal saa expression in GF larvae rescued neutrophil mobilization to the wound margin 6 hours after injury (Fig 5C). These data support a working model wherein Saa produced by IECs in response to the microbiota promotes local recruitment of neutrophils to the intestine while limiting aberrant systemic neutrophil activation by the microbiota, thereby allowing neutrophils to effectively respond to subsequent challenges such as peripheral injury (Fig 5D).

Discussion

Previous reports in mice and zebrafish have demonstrated that the microbiota influences diverse aspects of neutrophil biology, including increased abundance and longevity, enhanced wound recruitment, and elevated bacterial killing [10, 12, 34–36, 38, 39, 67]. However, few studies conclusively identify host or microbial factors that mediate the microbiota’s influences on neutrophils. Here we demonstrated that Saa induced following colonization transduces information regarding intestinal microbiota, regulating diverse aspects of neutrophil biology (Fig 5D). In colonized zebrafish at homeostasis, Saa promoted neutrophil maturation and mobilization to the intestine while suppressing systemic neutrophil abundance and pro-inflammatory gene expression. During inflammatory challenge, Saa restricted neutrophil anti-bacterial activity yet promoted neutrophil recruitment to a wound. While in vitro studies suggest rhSAA (rhSAA; Peprotech, Rocky Hill, NJ, USA) promotes pro-inflammatory cytokine (TNFα, IL-1β, and IL-8) and ROS production [23, 68–72], the pro-inflammatory effects of rhSAA on granulocytes contrasts with reports using purified endogenous SAA [28, 29]. Thus, our in vivo analysis clarifies conflicting reports of Saa’s effects on neutrophils.

Our results establish Saa as a major host factor mediating microbial control of neutrophil function. Expression of saa is negligible in germ-free zebrafish larvae, but is potently induced following microbial colonization [7, 10]. We therefore reasoned that neutrophil defects in colonized saa larvae would phenocopy those exhibited in WT germ-free larvae [10, 12, 73, 74]. In accord, we observed decreases in neutrophil migratory behavior, intestinal association, and wound recruitment in GF WT and colonized saa larvae. We demonstrated intestinally-derived Saa was sufficient to partially complement neutrophil deficiencies in saa zebrafish and could restore neutrophil wound recruitment in defects GF WT larvae. These results indicate that intestinally-derived Saa conditions neutrophils in vivo following microbiota colonization. As we have previously shown that colonization with different bacterial taxa leads to varied levels of saa induction [11], this signaling axis may facilitate taxa-specific effects on host innate immune development. Further, it is possible that mammalian SAA paralogs similarly condition neutrophils, as a subset of Saa genes are induced in mouse intestinal tissues following microbial colonization [7–9, 64]. As broad-spectrum antibiotics usage can dramatically impact the microbiota [75-77] and antibiotic treatment results in reduced intestinal Saa in mice [78], our model predicts antibiotic treatment could be associated with Saa-mediated aberrations in neutrophil function. We speculate that secondary infections that can occur following antibiotic use [79] could be due in part to concomitant alterations in SAA production.

The mammalian SAA gene family is comprised of both constitutively-expressed and acutely-inducible forms. There are four described paralogs in humans and mice, with much focus placed upon the acute forms which are upregulated following inflammatory stimuli as part of the acute phase response [13]. Following injury or trauma, hepatic production of the canonical acute SAA paralogs, SAA1 and SAA2, is dramatically induced. These secreted proteins enter circulation and associate with high density lipoproteins (HDL) [80], although the functional consequences of this interaction remain unclear. Notably, these acute forms are not only expressed in the liver, as they have also been detected in tissues including the small and large intestine where IECs produce them in response to microbial and inflammatory stimuli [7, 8, 26, 64, 81–83]. In contrast, SAA4 is constitutively expressed in the liver and is not induced by the microbiota or injury [13, 81, 84]. Our previous and present analyses indicate that the single zebrafish saa homolog is inducible by microbial stimuli with negligible constitutive expression in germ-free animals, suggesting zebrafish have minimal or no constitutive Saa activity.

Whole animal Saa1Saa2 double knockout mice (which still possess Saa3 –an acute extra-hepatic SAA that is a pseudogene in humans) and whole animal Saa3 knockout mice are both sensitive to chemically-induced enterocolitis [64, 85], suggesting SAA protects against intestinal inflammation. Moreover, following colonization with segmented filamentous bacteria (SFB), induction of SAA1/2 in mouse small intestinal epithelial cells stimulates expression of effector cytokines (IL-17) in T helper cells and promotes Th17 cell expansion and mucosal defense [26, 27, 86]. Our data indicate that IEC-derived Saa is important for neutrophil recruitment to the larval zebrafish intestine, which may play an important homeostatic role in host defense, mitigating breaches in barrier by commensal or pathogenic microbes. Mouse and human SAAs have also been described as retinol binding proteins that can bind dietary vitamin A with important implications in bacterial infections [8]. Furthermore, expression of SAA in IECs and the liver of mice requires dietary vitamin A [8]. We speculate that vitamin A, the microbiota, and perhaps other environmental factors interactively regulate SAA expression and subsequently influence the development and function of the innate immune system.

In this study we uncover tissue specific influences of Saa on both local and systemic neutrophil biology. Although our transgenic lines express saa at levels above endogenous saa transcript (Fig 4A and 4I, S6A and S6B Fig), we believe these transgenic levels are still within the physiologic range as circulating levels of Saa protein can reach 1 mg/mL [13-15]. Using these tissue specific promoters allowed us to disentangle the contribution of different tissue sources of Saa on host neutrophil responses. Our data indicate that Saa from the intestine and liver can impact neutrophil activities, such as promoting neutrophil recruitment to peripheral injury and suppressing bacterial killing. However, only Saa produced from IECs was sufficient to promote local neutrophil recruitment to the intestine. We speculate that both intestinal and hepatic Saa shape systemic neutrophil responses to injury by entering circulation from their respective origins. In zebrafish, microbiota induction of Saa is NF-κB and Myd88 dependent [10], and microbial factors, including LPS and flagellar function, are sufficient to induce Saa [11, 87]. Future studies aimed at delineating tissue specific induction of Saa by microbial signals are needed to understand how specific bacterial factors promote Saa production. Differential activity of SAA from the intestine versus liver may be due in part to differences in protein complex formation and downstream receptor recognition, as mammalian SAAs are known to associate with other proteins in HDL particles [88]. Our model suggests that in cases of chronic intestinal inflammation (e.g. inflammatory bowel disease) when SAA is expressed at high levels in the intestine, there may be SAA-mediated impacts on host systemic innate immune function [16-18].

Saa is thought to directly and indirectly influence in epithelial homeostasis in the mammalian intestine through direct bactericidal activity in the intestinal lumen and by shaping host mucosal innate and adaptive immunity [27, 64, 65]. Studies in mice have demonstrated that acute SAAs are expressed in IECs, and can be secreted both into the lumen and basolaterally, potentially aiding in the clearance of pathogenic bacteria [8, 64, 81, 83]. We did not detect significant alterations in intestinal microbiota composition in our saa mutant zebrafish when co-housed with WT controls. However, it remains possible that saa mutation may impact microbiota composition or density when genotypes are housed separately.

Our discovery that Saa functions in the microbiota-neutrophil axis motivates interest in underlying mechanisms. Neutrophil functions are largely regulated through “priming”, whereby microbial products (e.g. LPS, peptidoglycan, and flagellin) and host factors (e.g. TNFα, IL-1β, IL-8, GM-CSF, and ATP) signal to neutrophils, preparing them for response to additional stimuli [32, 89]. Intriguingly, many priming factors are induced in the zebrafish and mammalian intestine by the microbiota (e.g. TNFα, IL-1β, and GM-CSF) [7, 8, 10, 60, 81]. Primed neutrophils exhibit enhanced bacterial killing, altered motility, and transcriptional changes [33]. Once believed to be transcriptionally quiescent, several studies have reported neutrophil transcriptional responses to various priming stimuli in vitro [43–45, 47] and in vivo [46, 90, 91]. We demonstrate Saa restricts neutrophil pro-inflammatory gene expression and show induction of pro-inflammatory genes in neutrophils from colonized vs GF WT zebrafish, which is augmented in saa larvae. Consistent with our gene expression results, Saa limits primed neutrophil phenotypes of ROS production and bacterial killing in vivo and ex vivo. We propose Saa is upregulated following microbiota colonization to temper aberrant neutrophil priming by microbial products and production of host inflammatory mediators, thus limiting collateral damage to host tissues [92-94]. Similarly, a recent report using Saa3 knockout mice demonstrated that SAA3 suppresses bone marrow derived dendritic cell response to LPS [95], mirroring the anti-inflammatory effects of Saa on zebrafish neutrophils demonstrated here. It will be interesting to determine if SAA’s pleiotropic effects on diverse cell types are mediated by different receptors, oligomeric state, or binding other molecules such as retinol [96].

Our data reveal that Saa promotes neutrophil maturation in adult zebrafish. A recent study showed that mature neutrophils exhibit increased motility and response to stimuli [97]. While it is possible that altered neutrophil maturation may underlie the phenotypes we observe in saa mutants, some maturation-associated phenotypes (e.g. ROS production and bacterial killing) were elevated in neutrophils from saa zebrafish. Thus, Saa may differentially affect neutrophils at different stages of development and maturation.

Collectively, our findings highlight the importance of intestinal and hepatic Saa effecting systemic neutrophil development and function, suppressing their inflammatory tone and increasing mobilization to wounds. More broadly, our findings suggest that the ontogenetic and microbial control of priming factors is important for vertebrate immunological development.

Materials and methods

Ethics statement

Zebrafish studies were approved by the Institutional Animal Care and Use Committees of Duke University Medical Center (protocol number A115-16-05) in accordance with the Public Health Service Policy on the Human Care and Use of Laboratory Animals under the Unites States of America National Institutes of Health (NIH) Office of Laboratory Animal Welfare (OLAW).

Animal husbandry

All zebrafish lines were maintained on a mixed Tübingen (Tü) / TL background on a 14:10 hour light:dark cycle in a recirculating aquaculture system. From 5 dpf to 14 dpf, larvae were fed Zeigler AP100 <50-micron larval diet (Pentair, LD50-AQ) twice daily and Skretting Gemma Micro 75 (Bio-Oregon, B5676) powder once daily. Beginning at 14 dpf, larvae were fed Artemia (Brine Shrimp Direct, BSEACASE) twice daily, supplemented with a daily feed of Skretting Gemma Micro 75. From 28 dpf, the Gemma Micro 75 diet was replaced with Gemma Micro 300 (Bio-Oregon, B2809). At the onset of breeding age or sexual maturity, adult fish were given a 50/50 mix of Skretting Gemma Micro 500 (Bio-Oregon, B1473) and Skretting Gemma Wean 0.5 (Bio-Oregon, B2818) and one feeding of Artemia daily.

Larvae were also maintained on a 14:10 hour light:dark cycle in a 28.5°C incubator, and are of indeterminate sex. Gnotobiotic zebrafish were generated following natural mating and reared as described previously [98] with the following exception: GZM with antibiotics (AB-GZM) was supplemented with 50 μg/ml gentamycin (Sigma, G1264). Conventionally raised zebrafish were maintained at a density of ≤ 1 larva / mL, and at 3 dpf groups of 30 larvae were transferred to 10 cm petri dishes containing 20 mL gnotobiotic zebrafish media (GZM), inoculated with 3 mL filtered system water (5μm filter, SLSV025LS, Millipore) and fed autoclaved ZM-000 zebrafish diet (1% w/v stock concentration (in RO H2O), 0.0025% w/v final concentration, Zebrafish Management Ltd.) [10]. TgBAC(cldn15la:EGFP) (which expresses a Cldn15la-EGFP fusion protein), Tg(lyz:DsRed) and Tg(lyz:GFP) have been previously characterized [40, 58].

Zebrafish mutagenesis

Targeted deletion of the saa gene was performed using CRISPR/Cas9 gene editing targeting the second exon of saa as described [7]. Briefly, the guide RNA sequence was designed using the “CRISPR Design Tool” (http://crispr.mit.edu/). Guide RNA oligos (S1 Table, primers P1 and P2) were ligated into pT7-gRNA plasmid (Addgene, 46759), which, following BamHI (New England Biolabs, R0136L) digest, was in vitro transcribed using MEGAshortscript T7 kit (ThermoFisher, AM1354) [99]. Cas9 mRNA was generated from XbaI (New England Biolabs, R0145S) digested pT3TS-nls-zCas9-nls plasmid (Addgene, 46757), and in vitro transcribed using mMESSAGE mMACHINE T3 kit (ThermoFisher, AM1348) [99]. A cocktail consisting of 150 ng/μL of nls-zCas9-nls and 120 ng/μL of gRNA, 0.05% phenol red, 120 mM KCl, and 20 mM HEPES (pH 7.0) was prepared, and approximately 1–2 nL was injected directly into the cell of one cell stage Tü zebrafish embryos. Mutagenesis was initially screened using Melt Doctor High Resolution Melting Assay (ThermoFisher, 4409535), and identified three independent alleles which were confirmed as deletions by Sanger sequencing of TOPO-cloned PCR products. Subsequent screening of the Δ22 (allele designation rdu60) was performed by PCR (primers P3 and P4, S1 Table) and products were resolved on 2% agarose TBE gels. Screening of the Δ2 and Δ5 alleles (allele designations rdu61 and rdu62 respectively) was performed by PCR amplification using primers P3 and P4 (S1 Table), followed by purification and Hha1 digest (New England Biolabs, R0139S). All mutant alleles (rdu60, rdu61, and rdu62) result in the loss of a single Hha1 restriction site present in the WT sequence.

Zebrafish transgenesis

For generation of transgenic zebrafish expressing saa in intestinal epithelial cells (IECs) or liver hepatocytes, the following strategy was employed utilizing Tol2 mediated transgenesis (https://onlinelibrary.wiley.com/doi/abs/10.1002/dvdy.21343). A 349 bp region of the zebrafish cldn15la promoter was PCR amplified from Tü genomic DNA using primers P5 and P6 (S1 Table), digested with FseI and AscI (New England Biolabs, R0588 and R0558), and ligated into the p5E 381 vector that had been linearized with the same enzymes to generate p5E-0.349cldn15la. The -2.8 kb fabp10a p5E vector was kindly provided by Brian Link. The full-length zebrafish saa coding sequence was PCR amplified from Tübingen whole-larvae cDNA using primers P7 and P8 (S1 Table) and subcloned into the plasmid pENTR-AleI using In Fusion (Takara Bio, 638909) to generate pME-saa. Both p5E-0.349cldn15la and pME-saa were verified by PCR and Sanger sequencing. A 4-way Gateway LR reaction was performed using LR Clonase II (ThermoFisher, 12538120) to recombine p5E-0.349cldn15la or p5E-2.8fabp10a, with pME-saa, and p3E-polyA 229 into pDEST 395 (which contains a bicistronic cmlc2:EGFP reporter), yielding the following constructs: -0.349cldn15la:saa:polyA;cmlc2:EGFP and -2.8fabp10a:saa:polyA;cmlc2:EGFP. These plasmids were co-injected with transposase mRNA into Tü embryos at the single cell stage respectively, as described elsewhere [60]. Injected F0 larvae were subsequently screened for mosaic EGFP expression in the heart, raised to adulthood, and used to establish stable lines for three independent alleles of Tg(-0.349cldn15la:saa;cmlc2:EGFP) (rdu64, rdu67, and rdu68) and Tg(-2.8fabp10a:saa;cmlc2:EGFP) (rdu66). Experiments were conducted with both larval and adult zebrafish positive for rdu64, rdu67, rdu68 and rdu66 and non-transgenic siblings used as controls.

To confirm IEC specificity of the 349 bp cldn15la promoter fragment, we also recombined p5E-0.349cldn15la, pME-mCherry 386 (cytosolic mCherry), and p3E-polyA with pDEST 394, and injected this construct into single cell Tü embryos. Mosaic F0 larvae were raised to adulthood, screened for germline transmission, and used to establish stable lines for three independent alleles of Tg(-0.349cldn15la:mCherry:polyA). All three alleles displayed mCherry expression restricted to the intestine, and allele rdu65 was maintained for further analysis.

Caudal fin amputation assay

Larval zebrafish were anaesthetized with 0.75 mM Tricaine and mounted in 3% (w/v) methylcellulose (in GZM). Caudal fins were amputated using a surgical scalpel (Surgical Specialties Sharpoint, 72–2201) and fish were revived into 10 cm dishes containing either GZM (for experiments with conventionally reared larvae) or AB-GZM without gentamycin (for experiments with gnotobiotic larvae). At 15 minutes, 3 hours, and 6 hours post wounding, animals were euthanized and fixed in 4% PFA / 1x PBS overnight at 4°C on an orbital platform. Larvae were subsequently washed with 1x PBS 3–5 times at room temperature. Imaging was performed with a Leica M205 FA equipped with a Leica DFC 365FX camera using either GFP or mCherry filters. Lyz+ cells were enumerated at the tail wound margin using the Fiji Cell Counter plugin. Fish that were wounded at the notochord or were moribund were excluded from analysis.

In vivo imaging and cell tracking analysis

Larval zebrafish were anaesthetized and mounted in 100 μL 0.75% (w/v in GZM) low melt agarose (Fisher Scientific, BP165-25) with 0.6 mM Tricaine in 96-well clear-bottom black-walled plates (Greiner Bio-One, 655090), and overlaid with 100 μl GZM containing 0.375 mM Tricaine. For homeostatic behavioral analysis, time lapse imaging was performed for 10 or 15 minutes and frames acquired at 30s intervals on a Zeiss Axio Observer with a Photometrics Evolve EMCCD camera and a 5x objective (NA 0.16, WD 18.5 mm). For live imaging following caudal fin amputation, larvae were mounted as described above and imaged for 6 hours at 2 or 5 minute intervals on an inverted Zeiss Axio Observer Z1 microscope equipped with an Xcite 120Q light source (Lumen Dynamics), an MRm camera (Zeiss), and a 20X objective lens (NA 0.4, WD 7.9 mm). Fish that were wounded at the notochord, moribund, or damaged during mounting were excluded from analysis. Automated cell tracking was performed using the MTrack2 plugin for Fiji. For cell tracking following caudal fin amputation, a region of interest (ROI) was drawn from the posterior end of the notochord to the wound margin. For neutrophil behavior analysis in homeostasis, two distinct ROIs were analyzed. An ROI (396 pixels by 52 pixels, w x h) was drawn over the intestine ending at the cloaca, and another ROI (364 pixels by 152 pixels, w x h) was positioned dorsally to the intestine in the trunk. We subsequently filtered tracking results to include cells that were tracked for ≥ 3 consecutive frames. Speed and meandering index were calculated as previously described [12].

Epifluorescence stereomicroscopy

Larval zebrafish from pooled clutches were anaesthetized in 0.75 mM Tricaine and mounted in 3% methylcellulose. Between 15 and 30 fish / genotype were imaged on a Leica M205 FA microscope equipped with a Leica DFC 365FX camera using identical magnification and exposures. Neutrophil recruitment to the intestine was quantified from 8-bit images with Fiji software using the Cell Counter plug in as described previously [12].

For quantifying mCherry fluorescence mean grey values were measured in FIJI. Equal size ROIs were drawn in the liver, trunk, or eye. To quantify intestinal fluorescence, an ROI was drawn around the entire intestine. All mean grey values were normalized by subtracting background signal.

Light sheet microscopy

For in vivo light sheet microscopy, 6 dpf zebrafish larvae were anesthetized in 0.75 mM Tricaine in GZM, mixed with 1% LMP Agarose in GZM supplemented with 0.75 mM Tricaine, and drawn up into a glass capillary tube. The agarose was allowed to polymerize for at least 10 minutes prior to imaging. Agarose-embedded larvae were extruded from the capillary into an imaging chamber filled with GZM supplemented with 0.75 mM Tricaine and heated to 28.5°C. Single Plane Illumination Microscopy (SPIM) was performed with a Zeiss Lightsheet Z.1 Microscope equipped with a 20x aqueous immersion objective (1.0 NA, 2.4 mm WD). Two channel acquisition was performed with frame switching using 488 nm and 561 nm excitation lines, and z-series were acquired with a 7.32 μm interval. Image-processing was performed using 64-bit FIJI.

Fluorescence Activated Cell Sorting (FACS)

For FACS, replicate pools of 60–90 lyz+ larvae or 30 Tg(-0.349cldn15la:mCherry)+ larvae of the indicated genotypes were euthanized with 3 mM Tricaine and washed for 5 minutes with deyolking buffer (55 mM NaCl, 1.8 mM KCl and 1.25 mM NaHCO3). Larvae were transferred to gentleMACS “C” tubes (Miltenyi Biotec, 130-096-334) containing 5 mL Buffer 1 [HBSS supplemented with 5% heat-inactivated fetal bovine serum (HI-FBS, Sigma, F2442) and 10 mM HEPES (Gibco, 15630–080)]. Larvae were dissociated using a combination of enzymatic and mechanical disruption. Following addition of Liberase (Roche, 05 401 119 001, 5 μg/mL final), DNaseI (Sigma, D4513, 2 μg/mL final), Hyaluronidase (Sigma, H3506, 6 U/mL final) and Collagenase XI (Sigma, C7657, 12.5 U/mL final), samples were dissociated using pre-set program C_01 on a gentleMACS dissociator (Miltenyi Biotec, 130-093-235), then incubated at 30°C on an orbital platform at 75 RPM for 10 minutes. The disruption-incubation process was repeated 4–6 times, after which 400 μL of ice-cold 120 mM EDTA (in 1x PBS) was added to each sample. Following addition of 10 mL Buffer 2 [HBSS supplemented with 5% HI-FBS, 10 mM HEPES and 2 mM EDTA], samples were filtered through 30 μm cell strainers (Miltenyi Biotec, 130-098-458) into 50 mL conical tubes. Filters were washed with 10 mL Buffer 2, and samples were centrifuged at 1800 rcf for 15 minutes at room temperature. The supernatant was decanted, and cell pellets were resuspended in 500 μl Buffer 2. For CellROX (Invitrogen C10491 or C10444) labeling experiments, cells were then resuspended in 500 μL Buffer 2 and transferred to individual wells of a 24-well plate. CellROX was added to a final concentration of 1 μM, and samples were incubated for 45 minutes in the dark at 28.5°C on a tilting platform. Samples were transferred to FACS tubes (Falcon, 352052), and DNaseI (5 μg/mL; Sigma, D4513) and 7-AAD (Sigma, A9400, 5 μg/mL) were added. FACS was performed with either Beckman Coulter MoFlo XDP, Beckman Coulter Astrios, Becton Dickinson Diva, or Sony SH800S cell sorters at the Duke Cancer Institute Flow Cytometry Shared Resource. Single-color control samples were used for compensation and gating. Viable neutrophils or IECs were identified as 7-AAD- lyz+ or 7-AAD- mCherry+ respectively. Data were analyzed with FloJo v10 (Treestar, CA).

Ex vivo zebrafish neutrophil culture

Kidneys were dissected from adult male and female transgenic lyz+ zebrafish of the indicated genotypes. Fish were of standard length 25.72mm±1.18mm (mean±S.D.). Single cell suspensions were generated by enzymatic treatment of dissected kidneys with DNaseI (2 μg/mL) and Liberase (5 μg/mL final) with gentle agitation on a fixed speed orbital platform (VWR, 82007–202) for 20 minutes at room temperature. Enzymes were deactivated by the addition of EDTA as described above and cell suspensions were filtered through 30 μm filters and either stained with 1 μM CellROX for 1 hour at 28°C and/or stained with 7-AAD (5 μg/mL) or Propidium Iodide (PI, Sigma, P21493, 5 μg/mL). Viable (7-AAD- or PI-) lyz+ cells were collected into poly-D-lysine (Sigma, P7280-5MG, 33 μg/mL) coated 96 well black-wall clear-bottom plates (Corning, 3603) containing 100 μl RPMI1640 (Gibco, 11835030) supplemented with 10% HI-FBS at a density of 15,000 cells per well. A total of 6–9 kidneys per genotype were pooled and used to seed 4–6 wells. DsRed-expressing Escherichia coli MG1655 (pRZT3, [87]) or GFP-expressing Pseudomonas aeruginosa PA01 (pMF230, [62]) were cultured aerobically overnight shaking at 37°C in LB supplemented with either tetracycline (10 μg/mL) or carbenicillin (100 μg/mL), respectively. Bacteria (100 μL) were subsequently sub-cultured into 5 mL selective LB media and grown at 37°C with shaking under aerobic conditions to an OD600 of 0.7–1, diluted in sterile 1x PBS (Gibco, 14190) to a concentration of 104 bacterial / μL (E. coli) or 103 bacteria / μL (P. aeruginosa), and added isolated neutrophils at an MOI of 2 or 0.2 respectively. Neutrophils were co-cultured with bacteria for 2 to 4 hours at 28°C with gentle agitation in the dark. Serial dilutions of co-culture supernatants were prepared in sterile 1x PBS and plated on selective media. Plates were incubated aerobically at 28°C for 24 hours and CFUs were enumerated.

For imaging studies, neutrophils were collected as described above with the following modification: cells were collected into poly-D-lysine coated thin-bottom 96 well plates (Greiner, 655090). Neutrophils were imaged on a Zeiss 710 inverted confocal microscope with 10x (NA 0.45) or 63x oil objectives (NA 1.40) with or without addition of bacteria.

To measure neutrophil viability, 96 well black wall clear bottom plates containing 15,000 cells / well were incubated with or without E. coli MG1655 (MOI 2, grown as described above) for 3.5 hours. Propidium Iodide (PI) was added to a final concentration of 6 μg/mL to each well, and plates were incubated for an additional 30 minutes at 28.5°C before reading fluorescence at 535/620 (Ex/Em, nm) with a BioTek Synergy2 plate reader. As a positive control, lyz+ cells were incubated with Triton X-100 (20% v/v) for 3 hours, then incubated at 65°C for 30 minutes prior to PI staining. Data are shown as % of maximum PI signal.

Cytospins and scoring of neutrophil maturation

For morphological assessment, 15,000–30,000 viable lyz+ cells isolated by FACS from adult dissected kidneys as described above were sorted into 500 μl Buffer 2. Cytospins were performed immediately following collection with a Cytospin 3 (Shandon) by centrifuging cell suspensions for 3 minutes at 800 rcf. Slides were dried overnight at room temperature, fixed with absolute methanol, then stained with Wright Giemsa (Sigma, WG16) according to the manufacturer’s instructions. Slides were imaged with a Leica DMRA2 compound microscope with a PL APO 40x air objective (NA 0.85, WD 0.11mm) and Q Imaging Micropublisher Digital color camera. Twenty individual ROIs were imaged per genotype, and cells were classified based on distinct nuclear morphology by a blinded investigator as described [42].

Gene expression analysis

Pools of 20–30 whole larvae, dissected digestive tracts, or carcasses were collected into 1 mL of TRIzol (ThermoFisher, 15596026) and stored at -80°C. For experiments measuring transcripts from digestive tract versus carcass, dissections were performed as described previously [73] on 6 dpf larvae to remove the digestive tract (including intestine, liver, and pancreas). Digestive tracts were dissected and pooled into TRIzol (15–20 per replicate). The remaining carcasses were collected and pooled into TRIzol. Tissues were homogenized by passing samples 10–15 times through a 27-gauge needle. RNA was isolated following the manufacturer’s protocol with the following modification: a second wash with 70% ethanol (prepared in DEPC-treated H2O) was performed. For gene expression analysis of sorted neutrophils, viable lyz+ (4000–6000 cells / replicate) cells were collected into 750 μL TRIzol LS (ThermoFisher 10296010). For validation of the -0.349cldn15la promoter, 13,000 Tg(-0.349cldn15la:mCherry)+ and mCherry negative cells per replicate were collected into 750 μL of TRIzol LS. RNA was isolated using a NORGEN RNA Cleanup and Concentrator Micro-Elute Kit (Norgen Biotek, 61000), and samples eluted in 10 μL from which 8 μL of RNA was treated with DNaseI (New England Biolabs, M0303L). cDNA was synthesized using the iScript kit (Bio-Rad, 1708891). Quantitative PCR was performed in duplicate 25 μl reactions using 2X SYBR Green SuperMix (PerfeCTa, Hi Rox, Quanta Biosciences, 95055) run on an ABI Step One Plus qPCR instrument using gene specific primers (S1 Table). Data were analyzed with the ΔΔCt method.

In vivo bacterial infections

Pseudomonas aeruginosa PA01 carrying a constitutively expressed GFP plasmid (pMF230) [62] was grown in LB media supplemented with 100 μg/mL carbenicillin overnight shaking at 37°C. Overnight culture was concentrated to an OD600 of 5 (approximately 2x109 bacteria / mL) then frozen in aliquots. At 1 dpf, larvae were treated with 45 μg/mL 1-phenyl-2-thiourea (PTU) to inhibit melanization. Larvae were infected at 2 dpf by a genotype-blinded investigator. To achieve systemic infection, bacteria were injected into circulation via the caudal vein with borosilicate needles along with a phenol red tracer (3% w/v) as described previously [100], and any larvae that were damaged by injections were excluded from analysis. Approximately 150–300 CFU of P. aeruginosa PA01 GFP was injected per larvae. To enumerate CFUs in the inoculum, the injection dose was plated before and after infections on LB agar plates supplemented with 100 μg/mL carbenicllin. Immediately following infection larvae were screened for even dosing by fluorescence microscopy, and significantly under- or over-infected larvae were excluded from further analysis by an investigator blinded to genotype. To quantify bacterial burden, individual larvae were homogenized in 500 μL sterile 1x PBS using a Tissue-Tearor (BioSpec Products, 985370) at 24 hour intervals, and serial dilutions were plated on selective media (LB agar supplemented with 100 μg/mL carbenicillin). P. aeruginosa CFUs were enumerated following aerobic incubation for 24 hours at 28°C.

Zebrafish immunohistochemistry

For immunostaining, 6 dpf larval zebrafish were fixed in 4% PFA/1x PBS overnight at 4°C on an orbital platform, then washed 3-5x with 1x PBS the following day. Larvae were mounted in 4% low melt agarose in cryo-molds molds (Tissue-Tek) and 200 μm transverse sections cut using a Leica VT1000S vibratome. Sections were transferred to 24-well plates, washed 3x in 1x PBS at room temperature, and permeabilized with 1x with PBS containing 0.1% (v/v) Triton x-100 (PBS-T) for 30 minutes at RT. Sections were blocked in 5% (v/v) donkey serum in PBS-T for 1 hour at room temperature, then incubated with primary antibodies diluted in blocking buffer overnight at 4°C with agitation (mouse anti 4E8, Abcam, ab73643, 1:200; mouse anti 2F11, Abcam, ab71286, 1:200; rabbit anti DsRed, Clontech, 632496, 1:200; chicken anti GFP, AVES, GFP-10x0, 1:200). Sections were washed with PBS, then incubated with species-specific secondary antibodies [(ThermoFisher, A10042, A32728, A11039, 1:200) and Hoechst 33258 (ThermoFisher, H3569, 1:1000)] diluted in PBS-T for 2–4 hours at room temperature on a tilting platform. Sections were then washed 3x with 1x PBS then mounted and coverslipped on slides using mounting media containing DAPI (Vector Laboratories, Inc, H-1200). Slides were imaged with a Zeiss LSM 780 upright confocal microscope equipped with a GaAsP array detector using a 63x oil objective (NA 1.4, WD 0.19 mm).

16S rRNA gene sequencing

Adult heterozygous saa zebrafish were bred naturally and the resulting embryos were collected into system water and pooled at 0 dpf. Fertilized embryos were sorted at equal densities into autoclaved 3L tanks filled with system water at 1 dpf (50 embryos per 3 liter tank). Zebrafish of both WT and saa mutant genotypes were co-housed for the duration of the experiment. Larvae were maintained under static conditions until 6 dpf, at which time water flow and feeding were begun. Dissected digestive tracts and environmental water samples were collected at two time points: 6 dpf and 70 dpf. At 6 dpf zebrafish were sampled prior to first feeding, and fish sampled at 70 dpf were fed using the standard facility diet (as described above in animal husbandry) beginning at 6 dpf. For each timepoint, fish were sampled from a minimum of 5 tanks. For the 6 dpf timepoint DNA was isolated from 14 WT and 15 saa mutant intestines. For the 70 dpf timepoint, DNA was isolated from 15 WT and 12 saa mutant intestines. For water samples, 50 mL of water was collected from each tank and filtered using 0.2 μm MicroFunnel Filter Units (Pall Corporation, 4803). Filters were removed from filter units with sterile forceps, transferred to Eppendorf tubes and snap frozen in a dry ice/ethanol bath. For intestinal samples, digestive tracts were dissected and flash frozen and either carcasses or fin were reserved for genotyping. Fish samples were genotyped to identify homozygous saa and saa zebrafish prior to submission for genomic DNA extraction.

The Duke Microbiome Shared Resource (MSR) extracted bacterial DNA from gut and water samples using a MagAttract PowerSoil DNA EP Kit (Qiagen, 27100-4-EP) that allows for the isolation of samples in a 96 well plate format using a Retsch MM400 plate shaker. DNA was extracted from ≥12 fish per genotype per time point, and from 5 to 8 different tanks per timepoint to control for tank effects. Sample DNA concentration was assessed using a Qubit dsDNA HS assay kit (ThermoFisher, Q32854) and a PerkinElmer Victor plate reader. Bacterial community composition in isolated DNA samples was characterized by amplification of the V4 variable region of the 16S rRNA gene by polymerase chain reaction using the forward primer 515 and reverse primer 806 following the Earth Microbiome Project protocol (http://www.earthmicrobiome.org/). These primers (515F and 806R) carry unique barcodes that allow for multiplexed sequencing. Equimolar 16S rRNA PCR products from all samples were quantified and pooled prior to sequencing. Sequencing was performed by the Duke Sequencing and Genomic Technologies shared resource on an Illumina MiSeq instrument configured for 150 base-pair paired-end sequencing runs.

Subsequent data analysis was conducted in QIIME2 (https://qiime2.org), the successor of QIIME [101]. Paired reads were demultiplexed with qiime demux emp-paired, and denoised with qiime dada2 denoise-paired [102]. Taxonomy was assigned with qiime feature-classifier classify-sklearn [103], using a naive Bayes classifier, trained against the 99% clustered 16S reference sequence set of SILVA, v. 1.19 [104]. A basic statistical diversity analysis was performed, using qiime diversity core-metrics-phylogenetic, including alpha- and beta-diversity, as well as relative taxa abundances in sample groups. The determined relative taxa abundances were further analyzed with LEfSe (Linear discriminant analysis effect size) [105], to identify differential biomarkers in sample groups.

Statistical methods

All experiments were repeated at least two times and statistical analyses were performed with GraphPad Prism v.7. Data are presented as mean ± SEM. For comparisons between 2 groups a two tailed student’s t-test or Mann-Whitney test was applied. For comparisons between 3 or more groups, a one-way ANOVA with Tukey’s multiple comparisons test was used. For experiments with 2 independent variables, a two-way ANOVA was performed. Significance was set as p < 0.05, and denoted as: * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001. Sample sizes are indicated in the figure legends.

Generation of saa mutant zebrafish.

(A) Sequence of CRISPR/Cas9 induced deletions in the coding region of zebrafish saa exon 2. (B-C) In silico translation of mutant allele rdu60 revealed a frame shift mutation in exon 2 with a predicted early stop codon, indicated by the asterisk. (D) qRT-PCR of 6 dpf whole larvae demonstrated negligible saa mRNA expression in saa homozygous mutant zebrafish (saa or saa). (E) qRT-PCR of 6 dpf whole larvae from trans-heterozygous saa in-crosses revealed reduced saa mRNA levels relative to WT controls. (F) qRT-PCR of dissected digestive tracts from 6 dpf WT and saa larvae demonstrated significantly reduced saa mRNA levels. (G) qRT-PCR of 6 dpf whole larvae from homozygous saa in-crosses demonstrated reduced saa mRNA levels relative to WT controls. (qRT-PCR shown in panels D-G included 4–8 replicates / genotype, n ≥ 20 larvae / replicate). (H) Representative genotype distributions exhibited no deviation from expected Mendelian outcomes indicating no difference in viability of saa animals (p = 0.8257) (from two independent experiments, actual number of animals from each genotype overlaid on bars). (I-J) Morphometric analysis illustrated loss of saa does not impact growth [standard length (SL), height at anterior of anal fin (HAA)] in 6 dpf larvae (n ≥ 26 larvae). (K-L) Confocal micrographs of transverse sections from transgenic WT and saa 6 dpf larvae positive for TgBAC(cldn15la:EGFP) (which labels IEC basolateral membranes with a Cldn15la-GFP fusion protein), and immunofluorescence labeling with the brush-border antibody 4E8 demonstrated intestinal architecture is qualitatively normal in mutant animals (scale bar = 50 μm). (M) qRT-PCR of 6 dpf whole larvae following a tail amputation showed no significant induction of saa. (N) lyz:GFP+ neutrophil recruitment to caudal fin wound 6 hours following amputation in 6 dpf zebrafish larvae revealed decreased neutrophil recruitment in saa compound mutants (saa) versus WT controls (n ≥ 35 larvae / genotype at 6 hour time point). In panels D, F, G, I, J, M and N a t-test was used. In panel E, a Mann-Whitney test was applied. A chi-squared test was used to test for deviation from Mendelian distribution in H. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

(TIF)

Click here for additional data file.

Isolation and characterization of neutrophils from WT and saa mutant zebrafish larvae.

(A) Gating strategy for isolation of lyz:EGFP+ neutrophils from 6 dpf zebrafish larvae. (B) The mean fluorescence intensity (MFI) of the lyz:EGFP+ neutrophil population was not significantly different between WT and saa mutant larvae. (C) qRT-PCR revealed no significant difference in lysozyme C (lyz) transcript levels in sorted lyz+ neutrophils from WT and saa mutant larvae. (D) Quantification of intracellular ROS levels as indicated by CellROX staining measured by flow cytometry in lyz:EGFP+ neutrophils from WT and saa mutant larvae showed no significant difference. (E-F) qRT-PCR analysis of sorted neutrophils revealed no differential expression of genes associated with pro-myelocyte progenitors (cpa4, cpb1, cel.2) or apoptotic markers (cyba, p53, bcl2a, baxa) between WT or saa mutants. (For panels B—F: n ≥ 4 replicates / genotype, n = 60–90 larvae / genotype). In panels B-F data was analyzed by t-test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

(TIF)

Click here for additional data file.

Isolation and culture of zebrafish neutrophils ex vivo.

(A) Gating strategy for isolation of CellROX labeled lyz:DsRed+ neutrophils from adult zebrafish kidneys. (B-E) Low magnification (10x) confocal images of lyz:EGFP+ neutrophils labeled with CellROX ex vivo (scale bar = 200 μm). (F-G) Imaging of CellROX labeled lyz:EGFP+ neutrophils isolated from WT and saa mutant zebrafish revealed cytoplasmic punctae (indicated by white arrows). Red dashed box indicates region enlarged to show cytoplasmic CellROX punctae (Scale bar = 20 μm). (H) The mean fluorescence intensity (MFI) of the lyz:DsRed+ population was unchanged between WT, saa, or saa;Tg(cldn15la:saa) neutrophils (4 replicates / genotype). (I) Quantification of intracellular ROS levels as indicated by CellROX staining measured by flow cytometry of lyz:DsRed+ neutrophils illustrated no significant difference in baseline CellROX levels between genotypes (4 replicates / genotype). (J) Measurement of lyz:EGFP+ neutrophil viability as assessed by Propidium Iodide (PI) staining demonstrated no significant differences between genotypes after 4 hours of co-culture with E. coli (relative to maximum PI signal from lysed neutrophils) (4 replicates / genotype). In panels H-I, data were analyzed by one-way ANOVA with Tukey’s multiple comparisons test. In panel J, data were analyzed by t-test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

(TIF)

Click here for additional data file.

Characterization of the zebrafish cldn15la promoter used to drive intestine-specific saa expression.

(A) UCSC genome browser view of the zebrafish cldn15la gene locus with the translational start indicated by the red arrow. Pink bar represents the cloned 349 bp promoter region upstream of the cldn15la gene used to drive intestine-specific transgene expression. Tracks for vertebrate conservation, FAIRE-seq and motifs for transcription factors important of IEC gene expression programs (identified by HOMER) are shown below the locus [60]. (B) Expression pattern of endogenous cldn15la along the length of the intestine in adult zebrafish, as measured by microarray in Wang et al., 2010 [59]. (C) Representative stereoscope images of IEC specific cytosolic mCherry expression in 5 dpf Tg(-0.349cldn15la:mCherry) larvae compared to non-transgenic (NTG) controls (scale bar = 500 μm). White dashed line indicates the intestine. (D) Quantification of mCherry fluorescence in the indicated tissues of 6 dpf Tg(-0.349cldn15la:mCherry) and non-transgenic control larvae demonstrated intestine-restricted mCherry reporter activity (n = 7 larvae / genotype). (E) Representative confocal micrographs of immunostained transverse sections of Tg(-0.349cldn15la:mCherry) 6 dpf larvae along the anterior-posterior axis labeled with the absorptive cell brush border-specific antibody 4E8 illustrated transgene expression in absorptive enterocytes (scale bar = 20 μm). (F) Representative immunofluorescence images of transverse sections from Tg(-0.349cldn15la:mCherry) 6 dpf larvae stained with secretory cell-specific antibody 2F11 demonstrated weak expression in secretory cells, including enteroendocrine cells (outlined by red dashed lines) and goblet cells (outlined by green dashed lines) (scale bar = 20 μm). (G) Gating strategy for isolation of -0.349cldn15la:mCherry+ IECs from larval zebrafish. (H) qPCR for mCherry transcript in sorted Tg(-0.349cldn15la:mCherry) IECs showed enrichment as compared to negative fraction (13,000 –0.349cldn15la:mCherry+ or mCherry negative cells / replicate, 4 replicates / genotype, 30 larvae / replicate). Data in panel D analyzed by one-way ANOVA with Tukey’s multiple comparison’s test. Data in panel H analyzed by a Mann-Whitney test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

(TIF)

Click here for additional data file.

Comparison of Tg(-0.349cldn15la:mCherry) and TgBAC(cldn15la:EGFP) expression patterns.

(A) Widefield fluorescence images of IEC-specific cytosolic mCherry expression in 6 dpf double transgenic Tg(-0.349cldn15la:mCherry); TgBAC(cldn15la:EGFP) larvae demonstrated overlap in mCherry and GFP expression domains (scale bar = 500 μm). (B) Representative confocal micrographs of immunolabeled transverse sections from anterior (upper) and posterior (lower) intestinal segments of Tg(-0.349cldn15la:mCherry);TgBAC(cldn15la:EGFP) 6 dpf larvae revealed expression of both mCherry and EGFP in IECs (scale bar = 20 μm). (C) Representative maximum intensity projections from single plane illumination microscopy (SPIM) z-stacks of Tg(-0.349cldn15la:mCherry);TgBAC(cldn15la:EGFP) 6 dpf larvae (scale bar = 500 μm). (D-E) High magnification single slice lightsheet images of the intestinal epithelium in a representative 6 dpf Tg(-0.349cldn15la:mCherry);TgBAC(cldn15la:EGFP) larva.

(TIF)

Click here for additional data file.

Intestinal expression of zebrafish saa using cldn15la promoter.

(A) qRT-PCR analysis of dissected digestive tracts and carcasses (remaining tissues following removal of the intestine) from 6 dpf larvae of indicated genotypes revealed significant increase in saa expression in transgenic larvae is restricted to the gut (4 replicates / genotype, 20 larvae / replicate). [73] (B) qRT-PCR of saa from whole 6 dpf larvae of the indicated genotypes (n = 4 replicates / genotype, 15–20 larvae / replicate). (C) Enumeration of lyz:DsRed+ wound-associated neutrophils at 6 hours post amputation revealed intestinally-derived saa does not affect neutrophil recruitment in WT larvae (n ≥ 32 larvae / genotype at 6 hour time point). (D) CFU quantification of bacterial concentration following 4 hour co-culture of lyz:DsRed+ adult zebrafish neutrophils with E. coli (MOI 2) (3–6 replicates / genotype). Data in panel A analyzed by a two-way ANOVA with p values reported in the table. Data in panel B was analyzed with a Kruskal-Wallis test. Data in panels C-D analyzed by one-way ANOVA with Tukey’s multiple comparison’s test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

(TIF)

Click here for additional data file.

Transgenic Saa is insufficient to dampen transcriptional activation of neutrophils in saa mutant zebrafish larvae.

(A) FACS revealed no significant difference in abundance of lyz:DsRed+ neutrophils in 6 dpf saa larvae as compared to saa;Tg(cldn15la:saa) larvae (4 replicates / genotype / experiment, n ≥ 60 larvae / replicate, data pooled from 2 independent experiments). (B) qRT-PCR of pro-inflammatory mRNAs from lyz:DsRed+ neutrophils isolated from both saa and saa;Tg(cldn15la:saa) larvae revealed persistent transcriptional activation as compared to WT (4 replicates / genotype, n ≥ 60 larvae / replicate). (C) FACS demonstrated increased abundance of lyz:DsRed+ neutrophils in 6 dpf saa larvae which is restored to WT levels in saa;Tg(fabp10a:saa) larvae (4 replicates / genotype / experiment, n ≥ 60 larvae / replicate, data pooled from 2 independent experiments). (D) qRT-PCR of pro-inflammatory mRNAs from lyz:DsRed+ neutrophils isolated from both saa and saa;Tg(fabp10a:saa) larvae revealed persistent transcriptional activation as compared to WT (4 replicates / genotype, n ≥ 60 larvae / replicate). Data in panels A-D were analyzed by one-way ANOVA with Tukey’s multiple comparisons test. Data are presented as mean ± SEM. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.

(TIF)

Click here for additional data file.

Neutrophil recruitment defects in saa deficient larvae are microbiota dependent.

(A) lyz:EGFP+ neutrophil recruitment to caudal fin wound margin 6 hours after amputation in 6 dpf gnotobiotic GF WT and GF saa zebrafish larvae was not significantly different (n ≥ 33 larvae / genotype / condition at the 6 hour timepoint). Data analyzed by a two-way ANOVA.

(TIF)

Click here for additional data file.

Primers used in this study (qRT-PCR, cloning and genotyping).

(PDF)

Click here for additional data file.

Statistical analysis of alpha-diversity comparisons from 16S rRNA gene sequencing (Kruskal-Wallis test).

Bacterial communities in WT and saa-/- mutant (MUT) guts were not significantly different by any of these metrics, whereas WT/MUT guts were significantly different from their housing water (H2O) as expected.

(PDF)

Click here for additional data file.

Statistical analysis of beta-diversity comparisons from 16S rRNA gene sequencing (PERMANOVA).

Bacterial communities in WT and saa-/- mutant (MUT) guts were not significantly different by any of these metrics, whereas WT/MUT guts were significantly different from their housing water (H2O) as expected.

(PDF)

Click here for additional data file.