Comparison of static immersion and intravenous injection systems for exposure of zebrafish embryos to the natural pathogen Edwardsiella tarda
© van Soest et al; licensee BioMed Central Ltd. 2011
Received: 9 May 2011
Accepted: 17 October 2011
Published: 17 October 2011
The zebrafish embryo is an important in vivo model to study the host innate immune response towards microbial infection. In most zebrafish infectious disease models, infection is achieved by micro-injection of bacteria into the embryo. Alternatively, Edwardsiella tarda, a natural fish pathogen, has been used to treat embryos by static immersion. In this study we used transcriptome profiling and quantitative RT-PCR to analyze the immune response induced by E. tarda immersion and injection.
Mortality rates after static immersion of embryos in E. tarda suspension varied between 25-75%, while intravenous injection of bacteria resulted in 100% mortality. Quantitative RT-PCR analysis on the level of single embryos showed that expression of the proinflammatory marker genes il1b and mmp9 was induced only in some embryos that were exposed to E. tarda in the immersion system, whereas intravenous injection of E. tarda led to il1b and mmp9 induction in all embryos. In addition, microarray expression profiles of embryos subjected to immersion or injection showed little overlap. E. tarda-injected embryos displayed strong induction of inflammatory and defense genes and of regulatory genes of the immune response. E. tarda-immersed embryos showed transient induction of the cytochrome P450 gene cyp1a. This gene was also induced after immersion in Escherichia coli and Pseudomonas aeruginosa suspensions, but, in contrast, was not induced upon intravenous E. tarda injection. One of the rare common responses in the immersion and injection systems was induction of irg1l, a homolog of a murine immunoresponsive gene of unknown function.
Based on the differences in mortality rates between experiments and gene expression profiles of individual embryos we conclude that zebrafish embryos cannot be reproducibly infected by exposure to E. tarda in the immersion system. Induction of il1b and mmp9 was consistently observed in embryos that had been systemically infected by intravenous injection, while the early transcriptional induction of cyp1a and irg1l in the immersion system may reflect an epithelial or other tissue response towards cell membrane or other molecules that are shed or released by bacteria. Our microarray expression data provide a useful reference for future analysis of signal transduction pathways underlying the systemic innate immune response versus those underlying responses to external bacteria and secreted virulence factors and toxins.
In the last decade the zebrafish has been firmly established as a model for infectious diseases [1–4]. The increasing popularity of the zebrafish is due to its many useful characteristics. The embryos develop fast ex utero and are transparent, making it possible to follow infection in vivo. The real-time analysis of infection processes in this model is facilitated by the development of transgenic zebrafish lines with fluorescently marked immune cell populations that can be used in combination with differential fluorescently labeled pathogens [5–8]. In addition, reverse and forward mutagenesis screens are possible, as are antisense knock-down techniques using morpholinos.
Like all jawed vertebrates the zebrafish possesses an innate and adaptive immune system. Innate immunity forms the first line of defense against invading microorganisms. Humoral components of the innate immune system, such as complement and acute phase proteins, were shown to be expressed in embryos and larvae and could be induced by lipopolysaccharide (LPS) challenge or infection [9, 10]. The major cell types required for cell-mediated innate immunity, macrophages and neutrophils, also develop during the first days of zebrafish embryogenesis [11–13]. An essential step in innate immunity is the recognition of invading microorganisms by pattern recognition receptor families, the most well studied being the Toll-like receptor (TLR) family. The TLRs activate a signaling pathway leading to a cytokine response and the activation of antimicrobial defense genes . The TLR signaling components are highly conserved between zebrafish and humans [15, 16]. In adults the innate and adaptive immune systems are tightly connected, however in the zebrafish embryo there is a temporal segregation. Whereas innate immunity is functional as early as 1 day post fertilization (dpf) [11, 17, 18], adaptive immunity does not reach full maturity until approximately 4 weeks post fertilization [13, 19, 20]. This makes the zebrafish embryo a useful in vivo model to study vertebrate innate immunity separate from adaptive immunity .
Bacterial infection models that have been developed in zebrafish differ in mode and time of infection, inoculum size, pathogenicity and host response [2–4]. The most common method of infection is injection, with the caudal vein as injection site at 1 dpf or the yolk circulation valley at 2 dpf . Salmonella typhimurium, a mammalian pathogen, was shown to be lethal to zebrafish embryos after caudal vein injection of a low dose of 25-50 bacteria . In contrast, injection of E. coli or an LPS-mutant of Salmonella typhimurium (Ra-mutant) was not lethal and the bacteria were cleared efficiently by the embryonic innate immune system . Pseudomonas aeruginosa, a broad host range pathogen, capable of infecting plants, invertebrates, and vertebrates, was lethal after injection into the yolk circulation valley at 10-100-fold higher injection inocula than used for S. typhimurium, while Burkholderia cenocepacia was recently shown to cause a lethal infection upon intravenous injection at a dose of less than 10 bacteria [23–25]. At relatively high doses, also gram-positive bacteria such as Streptococcus and Staphylococcus species were shown to be capable of causing lethality upon injection in both adults and embryos [26–29]. Injection of embryos with Mycobacterium marinum does not lead to a lethal infection, but the immune system is unable to clear this bacterium, leading to a chronic infection. This chronic infection is characterized by aggregation of macrophages into granuloma-like structures similar to the tuberculous granulomas found in human tuberculosis patients . The different infection models were useful to study bacterial virulence factors and the response of the host immune system [3, 9, 30, 31].
For experimental screening, intravenous injection of zebrafish embryos is a relatively low throughput method. For high throughput analysis, such as mutant or drug screens, it is highly desirable to have an easier method of infection like static immersion. Thus far, the only bacterial pathogens that were reported to be capable of infecting zebrafish embryos without the need of injection are Edwardsiella tarda and Flavobacterium columnare [32, 33], which are Gram-negative naturally occurring fish pathogens. E. tarda is primarily known for infecting channel catfish, Japanese eel and flounder, in which it causes edwarsiellosis, a generalized septicemia. Pressley and colleagues showed that 24 hpf zebrafish embryos immersed for five hours in a suspension of E. tarda had a cumulative mortality rate of 31% after 14 days, compared to 11% in the control embryos . In addition, the zebrafish embryos showed peaks in the expression of tnfa and il1b at 2 and 4 hours post exposure, respectively. In adults, E. tarda is capable of causing infection by static immersion in combination with dermal abrasion .
The aim of this study was to compare the robustness of immersion and injection methods for treatment of 1-day-old zebrafish embryos with E. tarda and to identify marker genes that provide a reproducible read-out for the immune response. We set out with a microarray analysis of embryos subjected to immersion in E. tarda, and used E. coli and P. aeruginosa, both non-lethal in the immersion method, for comparison. Several markers were selected for a qPCR time-course analysis of the immersion method and for comparison with caudal vein injection. Marker expression analysis at single embryo level revealed high variation between individuals in response to static immersion. In contrast, qPCR and microarray analysis of single embryos that were systemically infected by caudal vein injection showed a consistent profile of strong activation of the proinflammatory marker genes il1b and mmp9. We conclude that the injection method is best suited for studying the innate immune response towards systemic infection, while the immersion system is useful for studying epithelial or other tissue responses towards cell membrane or other molecules that are shed or released by bacteria.
Survival of zebrafish embryos after immersion in E. tarda suspension
Microarray analysis of embryos subjected to the immersion system
Surprisingly, very few of the genes up-regulated in the zebrafish embryo after exposure to E. tarda were immune related. Although transient induction of il1b and tnfa was previously observed by Pressley et al. , no induction of these genes was detected in our microarray analysis. Furthermore, expression of mmp9, one of the most strongly induced markers after Salmonella infection , was only slightly up-regulated (1.4 times). In total only 21 genes showed 2-fold or higher levels of up-regulation (P < 1.0 E-4) after E. tarda exposure (Additional file 1, Table S1). Some of these genes have a possible immune-related function. The highest induced gene after exposure to E. tarda was cyp1a (9.8-fold induction), which encodes a cytochrome P450 enzyme known to be involved in the toxic response [34, 35]. As shown in Additional file 1, Table S1, this gene is also highly induced after P. aeruginosa and E. coli exposure. The second highest induced gene was zgc:154020 (6.8-fold), which shows 62.1% identity with immunoresponsive gene 1 (irg1) from Mus musculus, a gene with homology to bacterial methylcitrate dehydratase, which is up-regulated in murine macrophages after exposure to LPS, cytokines, and mycobacteria [36–39]. Zgc:154020 will hereafter be referred to as irg1-like (irg1l). Like cyp1a, irg1l was also highly up-regulated after P. aeruginosa and E. coli exposure. A third gene with a possible immune-related function is stanniocalcin 1 (stc1), which was only induced after E. tarda exposure (2.1-fold). Stanniocalcin is involved in Ca2+ homeostasis in fish [40, 41], but in humans has also been implicated in inflammatory responses [42–44].
To compare the responses of zebrafish embryos to immersion with the different bacterial strains, we performed a gene ontology analysis on all genes showing differential expression in the microarray analysis (Additional file 2, Table S2). In embryos immersed in P. aeruginosa PAO1 and PA14, and in E. coli, but not in embryos immersed in E. tarda, genes with the GO-term "response to stimulus" were significantly enriched. The largest group of up-regulated genes with this GO-term (61 genes) was observed in the case of immersion with P. aeruginosa PA14. Further analysis into the "response to stimulus" GO category revealed that in particular genes with the GO-term "response to stress" were up-regulated (41 genes in the case of PA14), while only few genes were associated with the GO-term "immune response" (6 genes in the case of PA14). An overview of the genes with the GO-term "response to stimulus" that were up-regulated in response to the different bacteria is given in Additional file 3, Table S3. The lack of induction of many of the known immune response genes after 5 hours of exposure to E. tarda suggests that at that time, tissue infection has not yet been established.
Time course analysis of marker gene expression in the immersion system
Immune response in single embryos after static immersion in E. tarda
The variability in mortality rates in the static immersion system, led us to hypothesize that not all embryos become systemically infected with this method. At 4 days after E. tarda immersion, none of the surviving embryos, even those that were close to dying, showed clear fluorescence of the mCherry marker plasmid. Subsequently, we plated individual surviving embryos for CFU counting. From five surviving embryos, of which three showed a slow heart beat indicative of approaching death, we obtained CFU counts of 140 to 690 per individual embryo. In contrast, the egg water medium of these embryos, kept individually in well plates, contained between 80,000 and 300,000 CFUs. It cannot be ascertained from CFU plating if the surviving embryos were actually infected with low numbers of bacteria or that the low CFU counts resulted from bacteria sticking to the surface epithelium of these embryos. However, it is clear that the surviving embryos did not carry heavy infections.
Immune response in single embryos after caudal vein injection of E. tarda
Microarray analysis of embryos infected by caudal vein injection
Microarray analysis was used to further characterize the immune response in response to microinjection of E. tarda bacteria and compare this with the previous microarray results of the immersion system and with our published data of the response to Salmonella typhimurium injection . Single infected and mock-injected embryos were analyzed at 8 hpi in triplicate. In gene ontology analysis we observed significant enrichment of the GO-terms "immune system process" and "response to stimulus" (Additional file 5, Table S4), whereas these GO-terms were not enriched in results of the immersion method (Additional file 2, Table S2). In addition, functional annotation using DAVID  showed significant enrichment of the KEGG pathways for apoptosis and for Toll-like receptor, adipocytokine, NOD-like receptor, insulin, MAP kinase, RIG-I-like receptor, ErbB, and Jak-Stat signaling. Manual annotation of the induced gene group showed several representatives of the categories complement activation and acute phase response, immune-related transcription factors and signaling components, cytokines and chemokines, apoptosis, and defense response (Additional file 6, Figure S2). In addition, many genes that were not previously linked to the immune response were differentially expressed, including genes involved in signal transduction, transporting activity and metabolism (Additional file 6, Figure S2). Out of 498 significantly regulated probes at 8 hpi (Additional file 7, Table S5), only 2 down-regulated probes (for vtg6 and an unannotated transcribed locus) and 1 up-regulated probe (for an unknown gene) were also significantly changed in the immersion system at 5 hpe. The microarray comparison supports that the transcriptional signatures of embryos subjected to immersion and injection are markedly different, although it should be noted that the immersion and injection data are not directly comparable due to a few hours difference in time to respond to the bacteria and in the developmental stage of assessment. The E. tarda injection microarray data were also compared with our previous microarray data set of intravenous Salmonella typhimurium infection of embryos at 2, 5 and 8 hpi . This comparison showed an overlap of 141 probes with significantly changed expression in response to both pathogens (Additional file 7, Table S5). These probes represented among others tnfb, il1b, cxcl-c1c, mmp9, ncf1, mxc, pglyrp5, hamp1 and several signal transduction (e.g. tlr5b, irak3, nfkbiaa, pim1, socs1/3a/3b) and transcription factor genes (e.g. atf3, elf3, fos, junb, irg9/11, rel, stat1) (Additional file 7, Table S5).
Zebrafish is being established as an alternative vertebrate model to murine models for infection research. To enable large scale mutant and chemical screening the development of an easily applicable infection test system is highly desired. In this report we studied the effectiveness and variability of treatment of zebrafish embryos by static immersion in Edwardsiella tarda, a method previously described by Pressley et al. , in comparison with the caudal vein injection method.
In order to perform large scale screenings, a model test system should be optimized for a reproducible response. Our results confirmed the ability of E. tarda to cause mortality in zebrafish embryos after static immersion. However, the mortality rate was highly variable between different experiments, ranging from 25 - 75%, comparable to the mortality rate of 31% reported by Pressley et al . In order to find a more reproducible readout, we performed microarray analysis on zebrafish embryos that had been exposed to E. tarda by static immersion. Surprisingly, only a small number of genes showed differential expression. In contrast, a much larger number of genes were regulated by immersion in bacterial suspensions of E. coli and P. aeruginosa strains PAO1 and PA14 that do not cause any mortality. In addition, very few immune-related genes were induced by immersion in E. tarda and we observed no induction of il1b and tnfa, which showed transient induction patterns between 2 and 12 hpi in the study of Pressley et al. .
Interestingly, cyp1a was highly induced by all tested bacteria. In E. tarda immersion experiments the induction of this gene preceded that of il1b and mmp9 induction. Our results suggest that this gene is not induced by direct exposure to the bacteria, but by released cell membrane components or other molecules. Expression of cyp1a was most strongly induced by P. aeruginosa. Cyp1a is known to be induced by toxic chemicals in vascular endothelium, but also in the epithelium of the gills [34, 35]. P. aeruginosa PAO1 and PA14 are known to secrete large amounts of toxins and protein virulence factors [47–50]. Since cyp1a belongs to the cytochrome P450 family, its induction might be involved in a detoxification response. The observation that many of the genes regulated by P. aeruginosa are associated with the GO term "response to stress", and the lack of enrichment of genes with the GO-term "immune system process" is consistent with a response to toxins rather than an immune response to systemic infection.
A further time-course analysis by qPCR of pools of embryos subjected to E. tarda immersion showed strong induction of il1b and mmp9 after 48 hours. In addition, the irg1l gene, one of the few immune-related genes identified in the microarray study, was also induced at later time points after exposure to E. tarda. The irg1l gene is homologous to mammalian irg1, expression of which in murine macrophages is induced by cytokines, agonists of TLR signalling, and by mycobacterial infections [36–39]. Sequence similarity of irg1l and mammalian irg1 with bacterial methylcitrate dehydratases suggests an important role in metabolism, but the function in vertebrates remains unknown. When we analyzed gene expression at the level of single embryos we observed that il1b and mmp9 were highly expressed in only one out of five treated embryos. Expression of irg1l was induced in all embryos, but only at a high level in the same embryos that also showed induction of il1b and mmp9. One possible explanation for the variable results of the static immersion assay is that embryos can individually differ in their resistance towards E. tarda. To test this, we compared the immersion system with intravenous injection of bacteria. In contrast to the relatively low and highly variable mortality rates that we observed with the immersion method, injection of bacteria resulted in a reproducible rate of 100% mortality within 2 days. Strong individual differences in levels of gene expression were also observed in the injection system, but nevertheless, induction of the proinflammatory marker genes il1b and mmp9 was positive in all embryos and their induction levels correlated with the dose of live bacteria injected. Furthermore, microarray experiments with single injected embryos showed a consistent profile of strong activation of proinflammatory and defense genes and regulatory genes of the immune response. The observed gene expression profiles are concordant with those observed for intravenous Salmonella typhimurium infection of embryos at similar time periods after injection . Detailed comparisons of the responses to E. tarda and S. typhimurium infections will be part of a follow-up study that will also address the function of essential immune regulators in these models.
Since only a subset of embryos in the immersion assay showed induction of immune response markers and mortality it is conceivable that only these embryos were systemically invaded by E. tarda bacteria or that non-responsive embryos were invaded by a much lower number of bacteria. Neither fluorescence monitoring nor CFU plating indicated that embryos become heavily infected close before dying. On the contrary, bacteria were present in high abundance in the egg water medium and only few were associated with dying embryos. It therefore remains uncertain whether infection or toxic insult is the actual cause of mortality in the immersion system. It is possible that the variable immune gene inductions and mortality rates resulted from slight epithelial damage to embryos that occurred during dechorionating and washing procedures, providing sites of entry for bacteria. Instead of exposure at 1 dpf, we used the same immersion protocol on embryos of 3 dpf, which is the developmental stage when the mouth opens and the gut begins to be colonized by environmental bacteria . We followed survival until 5 dpf, which is the time-point up to which larvae do not fall under the European animal experimentation law, but did not observe mortality within that time (data not shown).
Besides being more practical for high-throughput screening, an immersion system might be preferred as a more natural route of infection compared to injection methods. However, we conclude that the E. tarda immersion method as applied here on 1-day-old zebrafish embryos is not suitable to achieve reproducible systemic infection. Therefore, unless a more virulent strain can be identified, injection remains the preferred method of infection for screening purposes. On the other hand, the immersion system is shown to be useful for studying epithelial or other tissue responses towards cell membrane or other molecules that are shed or released by bacteria. An alternative solution for high-throughput screening of systemic infection is the use of robotic yolk injection system recently developed for Mycobacterium marinum infection . However, the wild type E. tarda FL6-60 strain used here causes early lethality after yolk injection (data not shown). The use of less virulent (wild type or mutant) strains might provide a solution for this problem. In any case, our gene expression profiling data sets will be necessary for comparisons to the immune response in such alternative yolk infection methods.
Zebrafish embryos proved to be remarkably resistant to becoming systemically infected after immersion in bacterial suspensions of E. tarda, whereas they are strongly susceptible to intravenous injection of this pathogen. While the microarray expression profile of intravenously infected embryos indicates a strong inflammatory response, the transcriptional signature of embryos subjected to immersion was markedly different. Our data suggest that most of the early transcriptional responses in the immersion system may reflect an epithelial or other tissue response towards cell membrane or other molecules that are shed or released by bacteria. Therefore, our studies on the expression analysis in the bacterial immersion system will be useful for future analysis of signal transduction pathways underlying responses to external bacteria and secreted putative virulence factors and toxins. Transient induction of the cytochrome P450 gene cyp1a was specifically observed in immersion experiments but not when embryos were systemically infected by injection. In addition, our identification of the irg1l gene as a rapid response factor to externally added bacteria deserves further study of the underlying signal transduction pathway as compared to systemic tissue responses. Although irg1l is also up-regulated during systemic infection, its expression kinetics in embryos immersed in E. tarda is very different from that of well-known inflammation genes such as il1b and mmp9. Considering the important function of epithelial cells in cross talk with cells of the innate immune system, as recently underscored by studies in zebrafish , further analysis of infection modes using the identified marker genes will help to better understand the systemic response of tissues toward an infection in a whole organism context.
Zebrafish were handled in compliance with the local animal welfare regulations and maintained according to standard protocols (zfin.org). An albino strain was used for all immersion and injection experiment, except for the microarray study of injected embryos that was performed with wild type zebrafish. Embryos were grown at 28.5-30°C in egg water (60 μg/ml Instant Ocean salts). For the duration of bacterial injections embryos were kept under anesthesia in egg water containing 0.02% buffered 3-aminobezoic acid ethyl ester (tricaine; Sigma-Aldrich).
Bacterial immersion and injection experiments
Edwardsiella tarda strain FL6-60 obtained from Dr. P. Klesius (USDA, Auburn, AL) is the identical strain as used in the study of Pressley et al. . Identity of this strain was confirmed by performing nucleotide sequencing of the entire genome using Illumina technology with a 180-fold coverage (Genbank accessions CP002154 and CP002155). FL6-60 was grown over night on tryptic soy agar (Difco) at 28°C and subsequently a liquid culture in tryptic soy broth (TSB, Difco) was inoculated and grown overnight at 28°C with shaking at 150 rpm. Pseudomonas aeruginosa PAO1 and PA14 and Escherichia coli were grown over night in Luria-Bertani broth (LB)  at 37°C. For immersion experiments bacterial cultures were centrifuged in 50 ml tubes and the pellet was subsequently suspended in egg water to a final 108 CFU/ml for E. tarda and E. coli, and 109 CFU/ml for P. aeruginosa. Embryos were dechorionated at 24 hpf by a 3-5 min pronase treatment (2 mg/ml in embryo medium prewarmed to 30°C) and left to recover for one hour in egg water. Subsequently pools of 20 embryos in 6-well plates were immersed in 5 ml of the bacterial suspension and incubated for 5 hours at 28°C. After 5 hours of incubation, the embryos were either snap-frozen in liquid nitrogen or transferred to a new 6-wells plate, washed 3 times in egg water, and incubated at 28°C in 5 ml of egg-water. For CFU plating experiments, embryos were kept individually in 2.5 ml of egg water in 24-well plates.
For caudal vein injection experiments, E. tarda labeled with mCherry  was washed and subsequently suspended in PBS (phosphate-buffered saline) to a final 108 CFU/ml. Embryos were manually dechorionated at 24 hpf. Approximately 200 CFUs of E. tarda were injected into the blood island after the onset of blood flow at 28 hpf, or PBS was injected as a control. After injection, embryos were kept at 28°C and snap-frozen in liquid nitrogen at the required time points.
RNA isolation from pools of embryos
Pools of 20 - 30 embryos were snap-frozen in liquid nitrogen and subsequently stored at -80°C. Embryos were homogenized in 1 ml of TRI reagent (Ambion), and subsequently total RNA was extracted according to the manufacturer's instructions. The RNA samples were incubated for 20 min at 37°C with 10 U of DNaseI (Roche Applied Science) to remove residual genomic DNA before purification using the RNeasy MinElute Cleanup kit (Qiagen) according to the RNA clean-up protocol. The integrity of the RNA was confirmed by lab-on-chip analysis using the 2100 Bioanalyzer (Agilent Technologies). Samples used for microarray analysis had an average RNA integrity number value of 9 and a minimum RNA integrity number value of 8.
RNA isolation from single embryos
The single embryo RNA isolation procedure was performed according to de Jong et al. . Embryos were individually snap-frozen in liquid nitrogen and subsequently stored at -80°C. A frozen embryo was crushed with a chilled pestle and homogenized in 300 μl of TRI reagent (Ambion). 60 μl of chloroform was added and the mixture was transferred to a 1.5 ml reaction tube containing 50 mg phase lock gel (Eppendorf) and incubated at room temperature for 5 minutes. The mixture was centrifuged at 12000 g at 4°C for 15 minutes, after which the aqueous phase was transferred to a fresh tube. 1 volume of 70% ethanol was added and the mixture transferred to a RNeasy MinElute Cleanup kit (Qiagen) column which was centrifuged 15 seconds at 8000 g. 500 μl RPE buffer from the kit was applied to the column and centrifuged 15 seconds at 8000 g. 500 μl 80% ethanol was applied to the column and centrifuged 2 minutes at 8000 g. The collection tube was replaced and the column centrifuged 5 minutes at 14000 g. 14 μl H2O was applied to the column and centrifuged 1 minute at 14000 g. The average amount of RNA isolated from a single embryo was 500 ng.
The microarray slides were custom-designed by Agilent Technologies as previously described . The slides contained in total 43,371 probes of a 60-oligonucleotide length.
Amino-allyl-modified amplified RNA (aRNA) was synthesized in one amplification round from total RNA using the Amino Allyl MessageAmp II aRNA Amplification kit (Ambion). The amount of total RNA used was 1 μg in experiments using RNA from pooled embryos and 400 ng in experiments using RNA from single embryos. Subsequently, 6 μg of amino-allyl-modified aRNA was used for coupling of monoreactive Cy3 and Cy5 dyes (GE Healthcare) and column purified. Samples from embryos immersed in E. tarda, E. coli, or P. aeruginosa suspensions or untreated control embryos were labeled with Cy5 and hybridized against a Cy3-labeled common reference that consisted of a mixture of all samples from the immersion experiments. E. tarda and control immersions were analyzed in triplicate using pools of 20 embryos and compared with single experiments of E. coli, P. aeruginosa PAO1 and P. aeruginosa PA14 immersion. For the E. tarda injection study, infected embryos and control embryos injected with the PVP-carrier solution were labeled with Cy5 and analyzed in triplicate against a Cy3-labeled common reference. Dual-color hybridization of the microarray chips was performed at ServiceXS according to Agilent protocol G4140-90050 version 5.7 (http://www.Agilent.com) for two-color microarray-based gene expression analysis.
Microarray data were processed from raw data image files with Feature Extraction Software 9.5.3 (Agilent Technologies). Processed data were subsequently imported into Rosetta Resolver 7.0 (Rosetta Biosoftware) and subjected to default ratio error modeling. The raw data were submitted to the Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo) under accession no. GSE28486. To compare samples from treatment groups to the control samples re-ratio analyses were performed using the Rosetta built-in re-ratio with common reference application. Data were analyzed at the level of UniGene clusters (UniGene build no. 105) and at probe level. Significance cut-offs for the ratios were set at 1.5-fold change at P < 10-4 for analysis at UniGene cluster level and P < 10-5 for analysis at probe level.
Gene ontology (GO) analysis was performed using the GeneTools eGOn v2.0 web-based gene ontology analysis software (http://www.genetools.microarray.ntnu.no) . KEGG pathway analysis was performed using DAVID tools for functional annotation (http://david.abcc.ncifcrf.gov/) . In addition, genes were manually annotated based on information in the ZFIN (zfin.org) and NCBI Entrez Gene databases, and PubMed abstracts.
cDNA synthesis and quantitative reverse transcriptase PCR
For RNA samples from pooled embryos, cDNA synthesis reactions were performed in a 20 μl mixture of 500 ng of RNA, 4 μl of 5x iScript reaction mix (Bio-Rad Laboratories), and 1 μl of iScript reverse transcriptase (Bio-Rad Laboratories). For RNA samples from single embryos, cDNA synthesis reactions were performed in a 10 μl mixture of 100 ng of RNA, 2 μl of 5x iScript reaction mix (Bio-Rad Laboratories), and 0.5 μl of iScript reverse transcriptase (Bio-Rad Laboratories). The reaction mixtures were incubated at 25°C for 5 min, 42°C for 30 min, and 85°C for 5 min.
Real-time PCR was performed using the Chromo4 Real-time PCR detection system (Bio-Rad Laboratories) according to the manufacturer's instructions. Each reaction was performed in a 25-μl volume comprised of 1 μl of cDNA, 12.5 μl of 2x iQ SYBR Green Supermix (Bio-Rad Laboratories), and 10 pmol of each primer. Cycling parameters were 95°C for 3 min to activate the polymerase, followed by 40 cycles of 95°C for 15 s and 59°C for 45 s. Fluorescence measurements were taken at the end of each cycle. Melting curve analysis was performed to verify that no primer dimers were amplified. All reactions were performed as technical duplicates. For normalization, peptidylprolyl isomerase A-like (ppial), which showed no changes over the infection time course series, was taken as reference. Results were analyzed using the ΔΔC t method. Sequences of forward and reverse primers are described in Additional file 8, Table S6.
We thank Dr. Philip Klesius (USDA, Auburn, AL) for providing us with E. tarda strain FL6-60, Roel de Haan for help with qPCR, and Davy de Wit and Ulrike Nehrdich for fish care. This work was supported by the European Commission 6th framework project ZF-TOOLS (LSHG-CT-2006-037220).
- Kanther M, Rawls JF: Host-microbe interactions in the developing zebrafish. Curr Opin Immunol. 2010, 22: 10-19. 10.1016/j.coi.2010.01.006.PubMedPubMed CentralView ArticleGoogle Scholar
- Meeker ND, Trede NS: Immunology and zebrafish: Spawning new models of human disease. Dev Comp Immunol. 2008, 32: 745-757. 10.1016/j.dci.2007.11.011.PubMedView ArticleGoogle Scholar
- Meijer AH, Spaink HP: Host-pathogen interactions made transparent with the zebrafish model. Curr Drug Targets. 2011, 12: 1000-1007.PubMedPubMed CentralView ArticleGoogle Scholar
- Sullivan C, Kim CH: Zebrafish as a model for infectious disease and immune function. Fish Shellfish Immunol. 2008, 25: 341-350. 10.1016/j.fsi.2008.05.005.PubMedView ArticleGoogle Scholar
- Hall C, Flores MV, Crosier K, Crosier P: Live cell imaging of zebrafish leukocytes. Methods Mol Biol. 2009, 546: 255-271. 10.1007/978-1-60327-977-2_16.PubMedView ArticleGoogle Scholar
- Meijer AH, van der Sar AM, Cunha C, Lamers GEM, Laplante MA, Kikuta H, Bitter W, Becker TS, Spaink HP: Identification and real-time imaging of a myc-expressing neutrophil population involved in inflammation and mycobacterial granuloma formation in zebrafish. Dev Comp Immunol. 2008, 32: 36-49. 10.1016/j.dci.2007.04.003.PubMedView ArticleGoogle Scholar
- Renshaw SA, Loynes CA, Trushell DMI, Elworthy S, Ingham PW, Whyte MKB: A transgenic zebrafish model of neutrophilic inflammation. Blood. 2006, 108: 3976-3978. 10.1182/blood-2006-05-024075.PubMedView ArticleGoogle Scholar
- Ellett F, Pase L, Hayman JW, Andrianopoulos A, Lieschke GJ: mpeg1 promoter transgenes direct macrophage-lineage expression in zebrafish. Blood. 2011, 117: e49-56. 10.1182/blood-2010-10-314120.PubMedPubMed CentralView ArticleGoogle Scholar
- Stockhammer OW, Zakrzewska A, Hegedus Z, Spaink HP, Meijer AH: Transcriptome profiling and functional analyses of the zebrafish embryonic innate immune response to Salmonella infection. J Immunol. 2009, 182: 5641-5653. 10.4049/jimmunol.0900082.PubMedView ArticleGoogle Scholar
- Wang Z, Zhang S, Wang G: Response of complement expression to challenge with lipopolysaccharide in embryos/larvae of zebrafish Danio rerio: Acquisition of immunocompetent complement. Fish Shellfish Immunol. 2008, 25: 264-270. 10.1016/j.fsi.2008.05.010.PubMedView ArticleGoogle Scholar
- Herbomel P, Thisse B, Thisse C: Ontogeny and behaviour of early macrophages in the zebrafish embryo. Development. 1999, 126: 3735-3745.PubMedGoogle Scholar
- Le Guyader D, Redd MJ, Colucci-Guyon E, Murayama E, Kissa K, Briolat V, Mordelet E, Zapata A, Shinomiya H, Herbomel P: Origins and unconventional behavior of neutrophils in developing zebrafish. Blood. 2008, 111: 132-141. 10.1182/blood-2007-06-095398.PubMedView ArticleGoogle Scholar
- Willett CE, Cortes A, Zuasti A, Zapata AG: Early hematopoiesis and developing lymphoid organs in the zebrafish. Dev Dyn. 1999, 214: 323-336. 10.1002/(SICI)1097-0177(199904)214:4<323::AID-AJA5>3.0.CO;2-3.PubMedView ArticleGoogle Scholar
- Mogensen TH: Pathogen recognition and inflammatory signaling in innate immune defenses. Clin Microbiol Rev. 2009, 22: 240-273. 10.1128/CMR.00046-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Meijer AH, Gabby Krens SF, Medina Rodriguez IA, He S, Bitter W, Ewa Snaar-Jagalska B, Spaink HP: Expression analysis of the Toll-like receptor and TIR domain adaptor families of zebrafish. Mol Immunol. 2004, 40: 773-783. 10.1016/j.molimm.2003.10.003.PubMedView ArticleGoogle Scholar
- Stein C, Caccamo M, Laird G, Leptin M: Conservation and divergence of gene families encoding components of innate immune response systems in zebrafish. Genome Biol. 2007, 8: R251-10.1186/gb-2007-8-11-r251.PubMedPubMed CentralView ArticleGoogle Scholar
- Davis JM, Clay H, Lewis JL, Ghori N, Herbomel P, Ramakrishnan L: Real-time visualization of Mycobacterium-macrophage interactions leading to initiation of granuloma formation in zebrafish embryos. Immunity. 2002, 17: 693-702. 10.1016/S1074-7613(02)00475-2.PubMedView ArticleGoogle Scholar
- Herbomel P, Thisse B, Thisse C: Zebrafish early macrophages colonize cephalic mesenchyme and developing brain, retina, and epidermis through a M-CSF receptor-dependent invasive process. Dev Biol. 2001, 238: 274-88. 10.1006/dbio.2001.0393.PubMedView ArticleGoogle Scholar
- Davidson AJ, Zon LI: The 'definitive' (and 'primitive') guide to zebrafish hematopoiesis. Oncogene. 2004, 23: 7233-7246. 10.1038/sj.onc.1207943.PubMedView ArticleGoogle Scholar
- Lam SH, Chua HL, Gong Z, Lam TJ, Sin YM: Development and maturation of the immune system in zebrafish, Danio rerio: a gene expression profiling, in situ hybridization and immunological study. Dev Comp Immunol. 2004, 28: 9-28. 10.1016/S0145-305X(03)00103-4.PubMedView ArticleGoogle Scholar
- Cui C, Benard EL, Kanwal Z, Stockhammer OW, Van der Vaart M, Zakrzewska A, Spaink HP, Meijer AH: Infectious disease modeling and innate immune function in zebrafish. Methods Cell Biol. 2011, 105: 273-308.PubMedView ArticleGoogle Scholar
- van der Sar AM, Musters RJP, van Eeden FJM, Appelmelk BJ, Vandenbroucke-Grauls CMJE, Bitter W: Zebrafish embryos as a model host for the real time analysis of Salmonella typhimurium infections. Cell Microbiol. 2003, 5: 601-611. 10.1046/j.1462-5822.2003.00303.x.PubMedView ArticleGoogle Scholar
- Brannon MK, Davis JM, Mathias JR, Hall CJ, Emerson JC, Crosier PS, Huttenlocher A, Ramakrishnan L, Moskowitz SM: Pseudomonas aeruginosa Type III secretion system interacts with phagocytes to modulate systemic infection of zebrafish embryos. Cell Microbiol. 2009, 11: 755-768. 10.1111/j.1462-5822.2009.01288.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Clatworthy AE, Lee JS-W, Leibman M, Kostun Z, Davidson AJ, Hung DT: Pseudomonas aeruginosa infection of zebrafish involves both host and pathogen determinants. Infect Immun. 2009, 77: 1293-1303. 10.1128/IAI.01181-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Vergunst AC, Meijer AH, Renshaw SA, O'Callaghan D: Burkholderia cenocepacia creates an intra-macrophage replication niche in zebrafish embryos, followed by bacterial dissemination and establishment of systemic infection. Infect Immun. 2010, 78: 1495-1508. 10.1128/IAI.00743-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin A, Loughman JA, Zinselmeyer BH, Miller MJ, Caparon MG: Streptolysin S inhibits neutrophil recruitment during the early stages of Streptococcus pyogenes Infection. Infect Immun. 2009, 77: 5190-5201. 10.1128/IAI.00420-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin B, Chen S, Cao Z, Lin Y, Mo D, Zhang H, Gu J, Dong M, Liu Z, Xu A: Acute phase response in zebrafish upon Aeromonas salmonicida and Staphylococcus aureus infection: Striking similarities and obvious differences with mammals. Mol Immunol. 2007, 44: 295-301. 10.1016/j.molimm.2006.03.001.PubMedView ArticleGoogle Scholar
- Neely MN, Pfeifer JD, Caparon M: Streptococcus-zebrafish model of bacterial pathogenesis. Infect Immun. 2002, 70: 3904-3914. 10.1128/IAI.70.7.3904-3914.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Prajsnar TK, Cunliffe VT, Foster SJ, Renshaw SA: A novel vertebrate model of Staphylococcus aureus infection reveals phagocyte-dependent resistance of zebrafish to non-host specialized pathogens. Cell Microbiol. 2008, 10: 2312-2325. 10.1111/j.1462-5822.2008.01213.x.PubMedView ArticleGoogle Scholar
- Phelps HA, Neely MN: Evolution of the zebrafish model: from development to immunity and infectious disease. Zebrafish. 2005, 2: 87-103. 10.1089/zeb.2005.2.87.PubMedView ArticleGoogle Scholar
- Volkman HE, Pozos TC, Zheng J, Davis JM, Rawls JF, Ramakrishnan L: Tuberculous granuloma induction via interaction of a bacterial secreted protein with host epithelium. Science. 2010, 327: 466-469. 10.1126/science.1179663.PubMedPubMed CentralView ArticleGoogle Scholar
- Pressley ME, Phelan PE, Witten PE, Mellon MT, Kim CH: Pathogenesis and inflammatory response to Edwardsiella tarda infection in the zebrafish. Dev Comp Immunol. 2005, 29: 501-513. 10.1016/j.dci.2004.10.007.PubMedView ArticleGoogle Scholar
- Chang MX, Nie P: RNAi suppression of zebrafish peptidoglycan recognition protein 6 (zfPGRP6) mediated differentially expressed genes involved in Toll-like receptor signaling pathway and caused increased susceptibility to Flavobacterium columnare. Vet Immunol Immunopathol. 2008, 124: 295-301. 10.1016/j.vetimm.2008.04.003.PubMedView ArticleGoogle Scholar
- Guiney PD, Smolowitz RM, Peterson RE, Stegeman JJ: Correlation of 2,3,7,8-tetrachlorodibenzo-p-dioxin induction of cytochrome P4501A in vascular endothelium with toxicity in early life stages of lake trout. Toxicol Appl Pharmacol. 1997, 143: 256-273. 10.1006/taap.1996.8051.PubMedView ArticleGoogle Scholar
- Jönsson ME, Brunström B, Brandt I: The zebrafish gill model: Induction of CYP1A, EROD and PAH adduct formation. Aquat Toxicol. 2009, 91: 62-70. 10.1016/j.aquatox.2008.10.010.PubMedView ArticleGoogle Scholar
- Basler T, Jeckstadt S, Valentin-Weigand P, Goethe R: Mycobacterium paratuberculosis, Mycobacterium smegmatis, and lipopolysaccharide induce different transcriptional and post-transcriptional regulation of the IRG1 gene in murine macrophages. J Leukoc Biol. 2006, 79: 628-638.PubMedView ArticleGoogle Scholar
- Degrandi D, Hoffmann R, Beuter-Gunia C, Pfeffer K: The proinflammatory cytokine-induced IRG1 protein associates with mitochondria. J Interferon Cytokine Res. 2009, 29: 55-68. 10.1089/jir.2008.0013.PubMedView ArticleGoogle Scholar
- Lee CG, Jenkins NA, Gilbert DJ, Copeland NG, O'Brien WE: Cloning and analysis of gene regulation of a novel LPS-inducible cDNA. Immunogenetics. 1995, 41: 263-270.PubMedView ArticleGoogle Scholar
- Shi S, Nathan C, Schnappinger D, Drenkow Jr, Fuortes M, Block E, Ding A, Gingeras TR, Schoolnik G, Akira S, et al: MyD88 primes macrophages for full-scale activation by Interferon-gamma yet mediates few responses to Mycobacterium tuberculosis. J Exp Med. 2003, 198: 987-997. 10.1084/jem.20030603.PubMedPubMed CentralView ArticleGoogle Scholar
- Tseng DY, Chou MY, Tseng YC, Hsiao CD, Huang CJ, Kaneko T, Hwang PP: Effects of stanniocalcin 1 on calcium uptake in zebrafish (Danio rerio) embryo. Am J Physiol Regul Integr Comp Physiol. 2009, 296: R549-557.PubMedView ArticleGoogle Scholar
- Wagner GF, Dimattia GE, Davie JR, Copp DH, Friesen HG: Molecular cloning and cDNA sequence analysis of coho salmon stanniocalcin. Mol Cell Endocrinol. 1992, 90: 7-15. 10.1016/0303-7207(92)90095-N.PubMedView ArticleGoogle Scholar
- Chakraborty A, Brooks H, Zhang P, Smith W, McReynolds MR, Hoying JB, Bick R, Truong L, Poindexter B, Lan H, Elbjeirami W, Sheikh-Hamad D: Stanniocalcin-1 regulates endothelial gene expression and modulates transendothelial migration of leukocytes. Am J Physiol Renal Physiol. 2007, 292: F895-904.PubMedView ArticleGoogle Scholar
- Kanellis J, Bick R, Garcia G, Truong L, Tsao CC, Etemadmoghadam D, Poindexter B, Feng L, Johnson RJ, Sheikh-Hamad D: Stanniocalcin-1, an inhibitor of macrophage chemotaxis and chemokinesis. Am J Physiol Renal Physiol. 2004, 286: F356-F362. 10.1152/ajprenal.00138.2003.PubMedView ArticleGoogle Scholar
- Wang Y, Huang L, Abdelrahim M, Cai Q, Truong A, Bick R, Poindexter B, Sheikh-Hamad D: Stanniocalcin-1 suppresses superoxide generation in macrophages through induction of mitochondrial UCP2. J Leukoc Biol. 2009, 86: 981-988. 10.1189/jlb.0708454.PubMedPubMed CentralView ArticleGoogle Scholar
- de Jong M, Rauwerda H, Bruning O, Verkooijen J, Spaink H, Breit T: RNA isolation method for single embryo transcriptome analysis in zebrafish. BMC Res Notes. 2010, 3: 73-10.1186/1756-0500-3-73.PubMedPubMed CentralView ArticleGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.PubMedView ArticleGoogle Scholar
- Darby C, Cosma CL, Thomas JH, Manoil C: Lethal paralysis of Caenorhabditis elegans by Pseudomonas aeruginosa. Proc Natl Acad Sci USA. 1999, 96: 15202-15207. 10.1073/pnas.96.26.15202.PubMedPubMed CentralView ArticleGoogle Scholar
- Galloway DR: Pseudomonas aeruginosa elastase and elastolysis revisited: recent developments. Mol Microbiol. 1991, 5: 2315-2321. 10.1111/j.1365-2958.1991.tb02076.x.PubMedView ArticleGoogle Scholar
- Peters JE, Galloway DR: Purification and characterization of an active fragment of the LasA protein from Pseudomonas aeruginosa: enhancement of elastase activity. J Bacteriol. 1990, 172: 2236-2240.PubMedPubMed CentralGoogle Scholar
- Wick MJ, Frank DW, Storey DG, Iglewski BH: Structure, function, and regulation of Pseudomonas Aeruginosa exotoxin A. Ann Rev Microbiol. 1990, 44: 335-363. 10.1146/annurev.mi.44.100190.002003.View ArticleGoogle Scholar
- Rawls JF, Mahowald MA, Goodman AL, Trent CM, Gordon JI: In vivo imaging and genetic analysis link bacterial motility and symbiosis in the zebrafish gut. Proc Natl Acad Sci USA. 2007, 104: 7622-7627. 10.1073/pnas.0702386104.PubMedPubMed CentralView ArticleGoogle Scholar
- Carvalho R, de Sonneville J, Stockhammer OW, Savage NDL, Veneman WJ, Ottenhoff THM, Dirks RP, Meijer AH, Spaink HP: A high-throughput screen for tuberculosis progression. PLoS ONE. 2011, 6: e16779-10.1371/journal.pone.0016779.PubMedPubMed CentralView ArticleGoogle Scholar
- Sambrook JaDR: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor NY, USA: Cold Spring Laboratory PressGoogle Scholar
- Lagendijk EL, Validov S, Lamers GEM, Weert Sd, Bloemberg GV: Genetic tools for tagging Gram-negative bacteria with mCherry for visualization in vitro and in natural habitats, biofilm and pathogenicity studies. FEMS Microbiol Lett. 2010, 305: 81-90. 10.1111/j.1574-6968.2010.01916.x.PubMedView ArticleGoogle Scholar
- Beisvag V, Junge F, Bergum H, Jolsum L, Lydersen S, Gunther CC, Ramampiaro H, Langaas M, Sandvik A, Laegreid A: GeneTools - application for functional annotation and statistical hypothesis testing. BMC Bioinformatics. 2006, 7: 470-10.1186/1471-2105-7-470.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.