Common and specific downstream signaling targets controlled by Tlr2 and Tlr5 innate immune signaling in zebrafish
© Yang et al. 2015
Received: 11 May 2015
Accepted: 30 June 2015
Published: 25 July 2015
Although the responses to many pathogen associated molecular patterns (PAMPs) in cell cultures and extracted organs are well characterized, there is little known of transcriptome responses to PAMPs in whole organisms. To characterize this in detail, we have performed RNAseq analysis of responses of zebrafish embryos to injection of PAMPs in the caudal vein at one hour after exposure. We have compared two ligands that in mammals have been shown to specifically activate the TLR2 and TLR5 receptors: Pam3CSK4 and flagellin, respectively.
We identified a group of 80 common genes that respond with high stringency selection to stimulations with both PAMPs, which included several well-known immune marker genes such as il1b and tnfa. Surprisingly, we also identified sets of 48 and 42 genes that specifically respond to either Pam3CSK4 or flagellin, respectively, after a comparative filtering approach. Remarkably, in the Pam3CSK4 specific set, there was a set of transcription factors with more than 2 fold-change, as confirmed by qPCR analyses, including cebpb, fosb, nr4a1 and egr3. We also showed that the regulation of the Pam3CSK4 and flagellin specifically responding sets is inhibited by knockdown of tlr2 or tlr5, respectively.
Our studies show that Pam3CSK4 and flagellin can stimulate the Tlr2 and Tlr5 signaling pathways leading to common and specific responses in the zebrafish embryo system.
The innate immune system is referred to as the first line in host defense against invading pathogens [1, 2]. Its highly developed ability to recognize microbial patterns and host-derived danger signals relies on so-called pattern recognition receptors (PRRs), especially on the Toll-like receptors (TLRs) [3–5]. In humans, the TLR family is composed of 10 members, which are located at the cell surface with the exception of TLR3, 7, 8 and 9, which are localized on intracellular endosomal membranes [6–8]. The TLRs are involved in the recognition of a wide variety of ligands, including pathogen associated molecular patterns (PAMPs), such as bacterial cell wall components and viral RNA, as well as damage associated molecular patterns (DAMPs) [9, 10]. This leads to subsequent intracellular signal transduction, triggering the production of inflammatory cytokines and chemokines, but it can also lead to anti-inflammatory responses as has been recently shown for TLR10 acting as a heterodimer with TLR2 [11, 12].
The recognition of PAMPs and DAMPs by different TLRs is directed by structurally conserved leucine-rich repeats (LRR) motifs of the TLRs ectodomains (ECDs) [13, 14]. Cell-surface TLRs can mediate binding to PAMPs by homodimerizing, like TLR4 and TLR5 that recognize lipopolysaccharide (LPS)  and flagellin , respectively. In contrast, other TLRs form heterodimers like TLR2 that associates with TLR1, TLR6 or TLR10, in conjunction recognizing lipopeptides and lipoproteins [12, 17, 18]. The diversity of TLR2 dimer combinations is thought to be responsible for its extensive recognition ability, ranging from the diverse components of various pathogens to the host heat-shock protein 70 . TLR2 plays an important role in resistance to the infection induced by Mycobacterium tuberculosis (Mtb) [19, 20]. For example, McBride et al. demonstrated that TLR2 knockout mice show enhanced cell infiltration and inflammation in lungs upon low dose chronic infection with Mtb, and fail to stably control the bacterial burden . A series of components of Mtb can trigger the TLR2 signaling pathway upon infection, such as triacylated lipoprotein LprG, LpqH and PhoS1 [22, 23], and glycolipid lipoarabinomannan . Pam3CysSerLys4 (Pam3CSK4) is a synthetic triacylated lipopeptide that mimics the triacylated lipoprotein of mycobacteria and classical gram-positive bacteria, which can be recognized by TLR2/TLR1 heterodimers and induce the production and release of pro- and anti-inflammatory cytokines (IL-6, IL-12, TNF-α and IL-10), chemokines (IL-8) and interferon (IFN-γ) [25–28]. Most of these studies on the responses of the TLR2 signaling pathway have been performed in cell culture systems. A notable exception is the reported transcriptome response of mouse mononuclear phagocytes recruited to lungs challenged with Pam3CSK4 as measured by micro-arrays . As another example 8 day old mice were treated with Pam3CSK4 and analyzed for the expression of several inflammatory genes using qPCR .
The zebrafish embryonic model has much potential to study the ligand specificity of TLRs at the organism level . Importantly, zebrafish offers the possibility to study the innate immune system separated from the adaptive immune system in their embryonic and larval stages (up to 3-4 weeks post fertilization) [31–33]. To date, a number of TLR signaling pathway mediators have been identified and studied in zebrafish such as the adaptor proteins Myd88, Tirap (Mal), Trif and Sarm1 [34–38], and the downstream signaling intermediates Irak and Traf6 [39, 40].
In work previously published by our group, we demonstrated that the function and regulation of the zebrafish homologs of human TLR5, tlr5a and tlr5b, are conserved with their mammalian counterparts. Both tlr5a and tlr5b are strongly up-regulated in response to Salmonella typhimurium infection . Furthermore, in the same study it was shown that knockdown of these two genes prevented or weakened the activation of genes for several inflammatory mediators like mmp9, cxcl-C1c, il1b and il8 upon flagellin stimulation .
In this study, we aimed to study TLR2 function, in comparison with TLR5, in zebrafish using transcriptome analysis. Injection into the blood stream of the tlr2/tlr1 ligand Pam3CSK4, was followed by transcriptome profiling to characterize key genes involved in the early response to this PAMPs. In addition, by comparing the transcriptome response towards treatment with flagellin, we were able to discriminate non-specific immune responsive genes from a set of genes which are regulated by tlr2 but not by tlr5.
Wild-type zebrafish of the AB/TL strain were handled in compliance with the local animal welfare regulations and maintained according to standard protocols (zfin. org). Embryos were raised in egg water (60 g/ml Instant Ocean sea salts) at 28.5 °C. For the duration of bacterial injections, embryos were kept under anesthesia in egg water containing 0.02 % buffered 3-aminobenzoic acid ethyl ester (Tricaine). The breeding of adult fish was approved by the local animal welfare committee (DEC) of the University of Leiden. All protocols adhered to the international guidelines specified by the EU Animal Protection Directive 2010/63/EU.
Morpholino oligonucleotides (Gene Tools) were diluted to desired concentrations in 1× Danieu’s buffer (58 mM NaCl, 0.7 mM KCl, 0.4 mM MgSO4, 0.6 mM Ca (NO3)2, 5.0 mM HEPES (pH 7.6)) containing 1× phenol red (Sigma-Aldrich). For knockdown experiments, tlr2 ATG-morpholino (tlr2 mo, Additional file 1: Table S1) was injected with the optimal concentration at 0.5 mM and 1 nl volume per embryo at 0 ~ 2 cell stage. tlr5 translation was blocked via injecting 1 nl of the tlr5a and tlr5b ATG-morpholinos at a dose of 0.1 mM and 0.5 mM at 0 ~ 2 cell stage, as previously published by Stockhammer and coworkers . Control embryos were injected with the standard control morpholino (Sc mo, Additional file 1: Table S1).
For the ligands injection assay, purified Pam3CSK4 (InvivoGen) and flagellin from S. typhimurium (Flagellin FliC VacciGrade™, Invitrogen) were respectively diluted to 1 mg/ml and 100 μg/ml in sterile water. For their administration, 1 nl of the ligands was injected into the blood stream at 27 hpf, and sterile water was used as control. Injections were perormed using a FemtoJet microinjector (Eppendorf) and a micromanipulator with pulled microcapillary pipettes.
RNA isolation, cDNA synthesis and qPCR
RNA was extracted using TRIzol Reagent (Life Technologies) and purified by column according to the manufacturer’s instructions of RNeasy MinElute Cleanup Kit (Qiagen). The concentration and quality of RNA were detected by NANODROP 2000/2000c (Thermo Scientific). 1 μg cDNA synthesis reactions and qPCR were performed as described in the manufacturer’s instructions (iScript™ cDNA Synthesis Kit and iQ™ SYBR® Green Supermix, BioRad) and normalized against the expression of ppial as a housekeeping gene . PCR analysis was performed using the following protocol: 95 °C 3 min, 40 cycles real time of 95 °C 15 s and 60 °C 45 s, and final melting curve of 81 cycles from 95 °C 1 min to 55 °C 10 s. Results were analyzed using the ΔΔCt method. Primer sequences used can be found in Additional file 1: Table S1.
Deep sequencing and data analyzing
Triplicates of 10 ~ 20 embryos of AB/TL or tlr morphants from three injection conditions, Pam3CSK4, flagellin or water injection, were homogenized in 500 μl of Trizol reagent (Qiagen). Total RNA was extracted and column-purified according to the manufacturer’s instructions of the RNeasy MinElute Cleanup Kit (Qiagen). The subsequent sample preparation and Illumina RNA sequencing were as previously described . RNA samples were treated with DNaseI (Life Technologies) to remove residual genomic DNA. RNA integrity was analyzed by Lab-on-a-chip analysis (Agilent, Amstelveen, The Netherlands). A total of 2 μg of RNA was used to make RNAseq libraries using the Illumina TruSeq RNA Sample Preparation Kit v2 (Illumina, Inc., San Diego, CA, USA). The manufacturer’s instructions were followed with the exception of two modifications. In the adapter ligation step, 1 μl, instead of 2.5 μl, adaptor was used. In the library size-selection step, the library fragments were isolated with a double Ampure XP purification with a 0.7× beads to library ratio (Beckman Coulter, Woerden, The Netherlands). The resulting mRNAseq library was sequenced using an Illumina HiSeq2500 Instrument (Illumina, Inc.) according to the manufacturer’s instructions with a read length of 2 × 50 nucleotides. Image analysis and base-calling were done using the Illumina HCS version 2.0.12. The raw data has been submitted to the GEO database (accession number GSE64570). The total number of reads for each sample is summarized in Additional file 2: Table S4 and quality control was according to the sequencing company guidelines (ZF-sceens.com). The data was analyzed using the GeneTiles software (http://www.genetiles.com)  using a cut-off setting of 2 fold-change and a p-value <0.01. In brief, Genetiles used fastq files as input for the program Bowtie2  to align the reads to the zebrafish genome (obtained from Ensembl version Zv9). Subsequently, the programs SAMtools , DESeq and DEXSeq [46, 47] are used for data processing. The complete data processing pipeline for Genetiles, including the used parameters, is available for download at www.genetiles.com and can also be found in Veneman et al . Using these settings we have mapped the numbers of reads as shown in Additional file 2: Table S4. The triplicate data sets of Pam3CSK4, flagellin and control treatments were mapped to 27104, 26583, and 26409 ENSEMBL genes, respectively. The difference between the mapped reads of the individual samples compared with the mapped reads of triplicate samples was always lower than 12 % (Additional file 2: Table S4). GO analysis was performed using the software package DAVID available at http://david.abcc.ncifcrf.gov/home.jsp .
The immune response of zebrafish embryos to injection of PAMPs in the caudal vein
The function of Tlr5 and Tlr2 in the immune response towards flagellin and Pam3CSK4
The function of tlr2 in recognition of Pam3CSK4 was tested in the same manner. Our results showed that tlr2 morphants did not exhibit up-regulation of il1b expression upon Pam3CSK4 stimulation, while the control morphants did (Fig. 2b). In contrast, tlr2 morphants stimulated with flagellin, still showed a significant induction of il1b expression (Fig. 2c).
Identification of a common response gene set for Pam3CSK4 and flagellin stimulation
Since comparisons of the transcriptome response to PAMPs that activate different TLRs has not been described before in a whole organism we decided to perform RNA deep sequencing (RNAseq) of embryos treated with flagellin and Pam3CSK4 at 1 hpi.
Identification of gene sets that are specifically regulated by Pam3CSK4 and flagellin
This study unexpectedly revealed that there is a relatively large group of genes that are only induced or repressed by either Pam3CSK4 or flagellin (Fig. 3b). To further test whether this specifically induced group is completely unaffected by the other PAMP treatment we used a rigorous filtering approach as shown in Fig. 3d and e. For this approach, the DEGs from the Pam3CSK4 (264 genes) and flagellin (306 genes) stimulation groups were compared with the group of genes that were not affected by flagellin (22611 genes) and Pam3CSK4 (22391 genes), respectively, with a cut-off setting at <1.4-fold change and p > 0.05. By taking the overlap of these sets we thereby exclude genes that were inducible by the other ligand even at very low stringency. This resulted in a set of 48 genes (Fig. 3d) for which the response is specific to Pam3CSK4 and a set of 42 genes (Fig. 3e) for which the response is specific to flagellin. GO analysis indicated that these two groups contain different categories (Fig. 3d, e). Most notably, genes with high fold-change (>2) from the Pam3CSK4 specific group include many transcription factors involved in the TLR signaling pathway, such as junbb, cebpb, fosb, fosl1a, egr3 and nr4a1. In the flagellin specific group of genes this is not the case and an obvious enriched category could not be identified.
Function of the tlr2 and tlr5a genes in the transcriptome responses to Pam3CSK4 and flagellin
Discussion and Conclusion
The signaling pathways underlying recognition of PAMPs have been studied intensively and this has led to a broad understanding of key regulators of innate immunity based on studies of cell cultures and the use of knockout rodent mutants. The new possibilities for analysis of transcriptomes using RNA deep sequencing make it highly attractive to analyze the responses of an entire test animal model at the system biology level. In this manuscript we have chosen the zebrafish embryo model for such an approach and have included functional analysis of Tlr5 and Tlr2 in the response towards two well-known PAMPs, flagellin and Pam3CSK4.
The results show that there is a relatively limited overlap between the transcriptome responses towards flagellin and Pam3CSK4 (Fig. 3b). The overlap includes well known downstream immune mediators that were previously shown to be induced by flagellin  such as il1b, tnfa, irak3, mmp9, cxcl-c1c and il8. In contrast, il6 and il10, that are associated with an anti-inflammatory response, were induced much stronger by Pam3CSK4 than by flagellin. A relatively much larger group of genes showed a differential response to flagellin or Pam3CSK4, including a group of genes of which the transcription is specific for activation by one of the two treatments (Fig. 4).
GO terms analysis of the genes specifically regulated by these two PAMPs show that there is an enriched category of transcription factors in the Pam3CSK4-specific group, of which most of genes are up-regulated and only one is down-regulated. Additionally, a less enriched category of immune response genes is found in this group as well, which include a down-regulated CRP (C-reactive protein). For the flagellin-specific group of genes, there are only five genes under the GO-term regulation of transcription of which two are down-regulated. Many of the genes specifically induced by Pam3CSK4 have also been shown to be strongly regulated by infection in the zebrafish embryo system  and therefore we would like to further study this group of genes in more detail in future research.
For functional analysis of the transcriptome response towards flagellin and Pam3CSK4 we used morpholinos that were selected on basis of their blocking effect on downstream signaling using qPCR and subsequently confirmed by RNAseq analysis. Surprisingly for tlr5 we found that a morpholino against each of the two copies of this gene, tlr5a and tlr5b had an effect on downstream signaling, with the tlr5a morpholino giving a complete block of induction of il1b by injection with flagellin and a partial effect of the tlr5b morpholino. These data suggests that these two tlr5 copies function in concert, perhaps by forming heterodimers.
We thank Ulrike Nehrdich, Davy de Witt and Laura van Hulst for fish caretaking and Hulya Ozupek, Dr. Hans Jansen and Dr. Ron Dirks (ZF-screens B. V.) for assistance with RNAseq analysis. We also want to thank to Dr. Jan de Sonneville and Kees-Jan van der Kolk for making the Genetiles software available.
S. Y. was supported by a grant from the China Scholarship Council (CSC) and R. M. J by the European Marie-Curie Initial Training Network FishForPharma (contract number PITN-GA-2011-289209).
- Ulevitch RJ. Molecular mechanisms of innate immunity. Immunol Res. 2000;21(2-3):49–54.PubMedView ArticleGoogle Scholar
- Masters SL, De Nardo D. Innate immunity. Curr Opin Immunol. 2014;26:v–vi.PubMedView ArticleGoogle Scholar
- Wright SD. Toll, a new piece in the puzzle of innate immunity. J Exp Med. 1999;189(4):605–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Kopp EB, Medzhitov R. The Toll-receptor family and control of innate immunity. Curr Opin Immunol. 1999;11(1):13–8.PubMedView ArticleGoogle Scholar
- Tan RS, Ho B, Leung BP, Ding JL. TLR cross-talk confers specificity to innate immunity. Int Rev Immunol. 2014;33(6):443–53.PubMed CentralPubMedView ArticleGoogle Scholar
- Matsumoto M, Funami K, Tanabe M, Oshiumi H, Shingai M, Seto Y, et al. Subcellular localization of Toll-like receptor 3 in human dendritic cells. J Immunol. 2003;171(6):3154–62.PubMedView ArticleGoogle Scholar
- Heil F, Ahmad-Nejad P, Hemmi H, Hochrein H, Ampenberger F, Gellert T, et al. The Toll-like receptor 7 (TLR7)-specific stimulus loxoribine uncovers a strong relationship within the TLR7, 8 and 9 subfamily. Eur J Immunol. 2003;33(11):2987–97.PubMedView ArticleGoogle Scholar
- Ahmad-Nejad P, Hacker H, Rutz M, Bauer S, Vabulas RM, Wagner H. Bacterial CpG-DNA and lipopolysaccharides activate Toll-like receptors at distinct cellular compartments. Eur J Immunol. 2002;32(7):1958–68.PubMedView ArticleGoogle Scholar
- Yu L, Wang L, Chen S. Endogenous toll-like receptor ligands and their biological significance. J Cell Mol Med. 2010;14(11):2592–603.PubMed CentralPubMedView ArticleGoogle Scholar
- Akira S, Takeda K. Toll-like receptor signalling. Nat Rev Immunol. 2004;4(7):499–511.PubMedView ArticleGoogle Scholar
- Striz I, Brabcova E, Kolesar L, Sekerkova A. Cytokine networking of innate immunity cells: a potential target of therapy. Clin Sci. 2014;126(9):593–612.PubMedView ArticleGoogle Scholar
- Oosting M, Cheng SC, Bolscher JM, Vestering-Stenger R, Plantinga TS, Verschueren IC, et al. Human TLR10 is an anti-inflammatory pattern-recognition receptor. Proc Natl Acad Sci U S A. 2014;111(42):E4478–84.PubMed CentralPubMedView ArticleGoogle Scholar
- Bell JK, Mullen GE, Leifer CA, Mazzoni A, Davies DR, Segal DM. Leucine-rich repeats and pathogen recognition in Toll-like receptors. Trends Immunol. 2003;24(10):528–33.PubMedView ArticleGoogle Scholar
- Gong J, Wei T. Structure modeling of Toll-like receptors. Methods Mol Biol. 2014;1169:45–53.PubMedView ArticleGoogle Scholar
- Poltorak A, He X, Smirnova I, Liu MY, Van Huffel C, Du X, et al. Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: mutations in Tlr4 gene. Science. 1998;282(5396):2085–8.PubMedView ArticleGoogle Scholar
- Hayashi F, Smith KD, Ozinsky A, Hawn TR, Yi EC, Goodlett DR, et al. The innate immune response to bacterial flagellin is mediated by Toll-like receptor 5. Nature. 2001;410(6832):1099–103.PubMedView ArticleGoogle Scholar
- Takeuchi O, Kawai T, Muhlradt PF, Morr M, Radolf JD, Zychlinsky A, et al. Discrimination of bacterial lipoproteins by Toll-like receptor 6. Int Immunol. 2001;13(7):933–40.PubMedView ArticleGoogle Scholar
- Takeuchi O, Sato S, Horiuchi T, Hoshino K, Takeda K, Dong Z, et al. Cutting edge: role of Toll-like receptor 1 in mediating immune response to microbial lipoproteins. J Immunol. 2002;169(1):10–4.PubMedView ArticleGoogle Scholar
- Bafica A, Scanga CA, Feng CG, Leifer C, Cheever A, Sher A. TLR9 regulates Th1 responses and cooperates with TLR2 in mediating optimal resistance to Mycobacterium tuberculosis. J Exp Med. 2005;202(12):1715–24.PubMed CentralPubMedView ArticleGoogle Scholar
- Ding A, Yu H, Yang J, Shi S, Ehrt S. Induction of macrophage-derived SLPI by Mycobacterium tuberculosis depends on TLR2 but not MyD88. Immunology. 2005;116(3):381–9.PubMed CentralPubMedView ArticleGoogle Scholar
- McBride A, Konowich J, Salgame P. Host defense and recruitment of Foxp3(+) T regulatory cells to the lungs in chronic Mycobacterium tuberculosis infection requires toll-like receptor 2. Plos Pathog. 2013;9(6):e1003397. doi:10.1371/journal.ppat.1003397.PubMed CentralPubMedView ArticleGoogle Scholar
- Drage MG, Tsai HC, Pecora ND, Cheng TY, Arida AR, Shukla S, et al. Mycobacterium tuberculosis lipoprotein LprG (Rv1411c) binds triacylated glycolipid agonists of Toll-like receptor 2. Nat Struct Mol Biol. 2010;17(9):1088–95. doi:10.1038/nsmb.1869.PubMed CentralPubMedView ArticleGoogle Scholar
- Drage MG, Pecora ND, Hise AG, Febbraio M, Silverstein RL, Golenbock DT, et al. TLR2 and its co-receptors determine responses of macrophages and dendritic cells to lipoproteins of Mycobacterium tuberculosis. Cell Immunol. 2009;258(1):29–37. doi:10.1016/j.cellimm.2009.1003.1008.PubMed CentralPubMedView ArticleGoogle Scholar
- Means TK, Wang S, Lien E, Yoshimura A, Golenbock DT, Fenton MJ. Human toll-like receptors mediate cellular activation by Mycobacterium tuberculosis. J Immunol. 1999;163(7):3920–7.PubMedGoogle Scholar
- Kanczkowski W, Morawietz H, Ziegler CG, Funk RH, Schmitz G, Zacharowski K, et al. Pam3CSK4 and LTA-TLRs ligands associated with microdomains induce IL8 production in human adrenocortical cancer cells. Horm Metab Res. 2007;39(6):457–60.PubMedView ArticleGoogle Scholar
- Lombardi V, Van Overtvelt L, Horiot S, Moussu H, Chabre H, Louise A, et al. Toll-like receptor 2 agonist Pam3CSK4 enhances the induction of antigen-specific tolerance via the sublingual route. Clin Exp Allergy. 2008;38(11):1819–29.PubMedGoogle Scholar
- Mintz M, Mintz D, Ilia-Ezra R, Shpigel NY. Pam3CSK4/TLR2 signaling elicits neutrophil recruitment and restricts invasion of Escherichia coli P4 into mammary gland epithelial cells in a murine mastitis model. Vet Immunol Immunopathol. 2013;152(1-2):168–75.PubMedView ArticleGoogle Scholar
- Agrawal S, Gollapudi S, Gupta S, Agrawal A. Dendritic cells from the elderly display an intrinsic defect in the production of IL-10 in response to lithium chloride. Exp Gerontol. 2013;48(11):1285–92.PubMedView ArticleGoogle Scholar
- Cabanski M, Wilhelm J, Zaslona Z, Steinmuller M, Fink L, Seeger W, et al. Genome-wide transcriptional profiling of mononuclear phagocytes recruited to mouse lungs in response to alveolar challenge with the TLR2 agonist Pam3CSK4. Am J Physiol Lung Cell Mol Physiol. 2009;297(4):L608–18.PubMedView ArticleGoogle Scholar
- Stridh L, Ek CJ, Wang X, Nilsson H, Mallard C. Regulation of Toll-like receptors in the choroid plexus in the immature brain after systemic inflammatory stimuli. Translational Stroke Res. 2013;4(2):220–7.View ArticleGoogle Scholar
- Willett CE, Cortes A, Zuasti A, Zapata AG. Early hematopoiesis and developing lymphoid organs in the zebrafish. Dev Dyn. 1999;214(4):323–36.PubMedView ArticleGoogle Scholar
- Davidson AJ, Zon LI. The ’definitive’ (and ’primitive’) guide to zebrafish hematopoiesis. Oncogene. 2004;23(43):7233–46.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(1):9–28.PubMedView ArticleGoogle Scholar
- van der Vaart M, van Soest JJ, Spaink HP, Meijer AH. Functional analysis of a zebrafish myd88 mutant identifies key transcriptional components of the innate immune system. Dis Model Mech. 2013;6(3):841–54.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu Y, Li M, Fan S, Lin Y, Lin B, Luo F, et al. A unique feature of Toll/IL-1 receptor domain-containing adaptor protein is partially responsible for lipopolysaccharide insensitivity in zebrafish with a highly conserved function of MyD88. J Immunol. 2010;185(6):3391–400.PubMedView ArticleGoogle Scholar
- Meijer AH, Gabby Krens SF, Medina Rodriguez IA, He S, Bitter W, Ewa Snaar-Jagalska B, et al. Expression analysis of the Toll-like receptor and TIR domain adaptor families of zebrafish. Mol Immunol. 2004;40(11):773–83.PubMedView ArticleGoogle Scholar
- Fan S, Chen S, Liu Y, Lin Y, Liu H, Guo L, et al. Zebrafish TRIF, a Golgi-localized protein, participates in IFN induction and NF-kappaB activation. J Immunol. 2008;180(8):5373–83.PubMedView ArticleGoogle Scholar
- Jault C, Pichon L, Chluba J. Toll-like receptor gene family and TIR-domain adapters in Danio rerio. Mol Immunol. 2004;40(11):759–71.PubMedView ArticleGoogle Scholar
- Stockhammer OW, Rauwerda H, Wittink FR, Breit TM, Meijer AH, Spaink HP. Transcriptome analysis of Traf6 function in the innate immune response of zebrafish embryos. Mol Immunol. 2010;48(1-3):179–90.PubMedView ArticleGoogle Scholar
- Phelan PE, Mellon MT, Kim CH. Functional characterization of full-length TLR3, IRAK-4, and TRAF6 in zebrafish (Danio rerio). Mol Immunol. 2005;42(9):1057–71.PubMedView 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(9):5641–53.PubMedView ArticleGoogle Scholar
- van Soest JJ, Stockhammer OW, Ordas A, Bloemberg GV, Spaink HP, Meijer AH. Comparison of static immersion and intravenous injection systems for exposure of zebrafish embryos to the natural pathogen Edwardsiella tarda. BMC Immunol. 2011;12:58.PubMed CentralPubMedView ArticleGoogle Scholar
- Marin-Juez R, Jong-Raadsen S, Yang S, Spaink HP. Hyperinsulinemia induces insulin resistance and immune suppression via Ptpn6/Shp1 in zebrafish. J Endocrinol. 2014;222(2):229–41.PubMedView ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Anders S, Reyes A, Huber W. Detecting differential usage of exons from RNA-seq data. Genet Res. 2012;22(10):2008–17.View ArticleGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols. 2008;4(1):44–57.View ArticleGoogle Scholar
- Veneman WJ, de Sonneville J, van der Kolk KJ, Ordas A, Al-Ars Z, Meijer AH, et al. Analysis of RNAseq datasets from a comparative infectious disease zebrafish model using GeneTiles bioinformatics. Immunogenetics. 2014;67(3):135–47.PubMed CentralPubMedView ArticleGoogle Scholar
- Ballouz S, Verleyen W, Gillis J. Guidance for RNA-seq co-expression network construction and analysis: safety in numbers. Bioinformatics. 2015;31(13):2123–30.PubMedView ArticleGoogle Scholar
- Veneman WJ, Stockhammer OW, de Boer L, Zaat SA, Meijer AH, Spaink HP. A zebrafish high throughput screening system used for Staphylococcus epidermidis infection marker discovery. BMC Genomics. 2013;14:255.PubMed CentralPubMedView ArticleGoogle Scholar
- Choi YJ, Jung J, Chung HK, Im E, Rhee SH. PTEN regulates TLR5-induced intestinal inflammation by controlling Mal/TIRAP recruitment. Faseb J. 2013;27(1):243–54.PubMed CentralPubMedView ArticleGoogle Scholar
- Bielinski SJ, Hall JL, Pankow JS, Boerwinkle E, Matijevic-Aleksic N, He M, et al. Genetic variants in TLR2 and TLR4 are associated with markers of monocyte activation: the atherosclerosis risk in communities MRI study. Hum Genet. 2011;129(6):655–62.PubMed CentralPubMedView ArticleGoogle Scholar
- Fang J, Fang D, Silver PB, Wen F, Li B, Ren X, et al. The role of TLR2, TRL3, TRL4, and TRL9 signaling in the pathogenesis of autoimmune disease in a retinal autoimmunity model. Invest Ophthalmol Vis Sci. 2010;51(6):3092–9.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.