Transcriptional profiling of the murine cutaneous response during initial and subsequent infestations with Ixodes scapularis nymphs
© Heinze et al; licensee BioMed Central Ltd. 2012
Received: 2 November 2011
Accepted: 6 February 2012
Published: 6 February 2012
Ixodes scapularis ticks are hematophagous arthropods capable of transmitting many infectious agents to humans. The process of blood feeding is an extended and continuous interplay between tick and host responses. While this process has been studied extensively in vitro, no global understanding of the host response to ticks has emerged.
To address this issue, we used PCR-arrays to measure skin-specific expression of 233 discrete genes at 8 time points during primary and secondary infestations of mice with pathogen-free I. scapularis nymphs. Selected results were then validated at the mRNA and protein levels by additional real-time PCR and bioplex assay.
Primary infestation was characterized by the late induction of an innate immune response. Lectin pattern recognition receptors, cytokines, and chemokines were upregulated consistent with increased neutrophil and macrophage migration. Gene ontology and pathway analyses of downregulated genes suggested inhibition of gene transcription and Th17 immunity. During the secondary infestation, additional genes were modulated suggesting a broader involvement of immune cells including CD8 and CD4 positive T lymphocytes. The cytokine response showed a mixed Th1/Th2 profile with a potential for T regulatory cell activity. Key gene ontology clusters observed during the secondary infestation were cell migration and activation. Matrix metalloproteinases were upregulated, apoptosis-related genes were differentially modulated, and immunoreceptor signaling molecules were upregulated. In contrast, transcripts related to mitogenic, WNT, Hedgehog, and stress pathways were downregulated.
Our results support a model of tick feeding where lectin pattern recognition receptors orchestrate an innate inflammatory response during primary infestation that primes a mixed Th1/Th2 response upon secondary exposure. Tick feeding inhibits gene transcription and Th17 immunity. Salivary molecules may also inhibit upregulation of mitogenic, WNT, Hedgehog, and stress pathways and enhance the activity of T regulatory cells, production of IL-10, and suppressors of cytokine signaling molecules (SOCS). This study provides the first comprehensive transcriptional analysis of the murine host response at the I. scapularis bite site and suggests both a potential model of the host cutaneous response and candidate genes for further description and investigation.
KeywordsTick(s) Ixodes scapularis Skin Cutaneous response Tick-host interface Gene expression profiling
The process of tick feeding activates a highly complex sequence of events at the bite site that facilitate the acquisition of a blood meal and create a suitable microenvironment for pathogen transmission and establishment . These events are governed by an array of salivary molecules secreted by the tick and the responses of the host to those molecules. It is a dynamic relationship with outcomes ranging from successful tick engorgement and potential pathogen transmission to tick rejection and greatly reduced pathogen acquisition. A critical factor that controls this variability is the host response to tick feeding. Laboratory animals with prior exposure to ticks may be significantly protected from pathogen acquisition from infected ticks ; after a single feeding with Dermacentor variabilis, rabbits develop an anti-tick immunity that greatly reduces successful blood feeding during future infestations . These observations suggest the host response to infestation may yield vital insights for tick and tick-borne disease control.
During the course of blood feeding, ticks have been shown to inhibit host pain/itch responses, hemostasis, angiogenesis, complement activation, and both innate and adaptive immune responses. In vitro experiments suggest tick saliva inhibits the production of cytokines (IL-1, IL-2, IL-6, IL-12, TNF-α, and IFN-γ) and adhesion molecules (ICAM-1, VCAM-1, P-selectin, LFA-1, and VLA-4) with the notable exception of IL-4 and IL-10 . The production of IL-4 in response to tick feeding has been supported in vivo . Tick salivary molecules also inhibit the function of immune cells present at the bite site. Salp15, an I. scapularis salivary protein, inhibits CD4-mediated activation of helper T-cells  and modulates dendritic cell activation through the lectin receptor DC-SIGN . Similarly, salivary gland disintegrin-like proteins ISL 929 and ISL 1373 inhibit neutrophil function  while salivary gland extracts have been shown to inhibit dendritic cell maturation, migration, and cutaneous turnover . Despite the ability of tick saliva to suppress host responses, some animals develop successful immunity dependent in part on T-cells, antibodies, complement, mast cells, and basophils . Piper and colleagues have compared the gene expression profile in skin [10, 11] and white blood cells  of tick resistant Bos indicus and tick susceptible Bos taurus cattle after multiple artificial and natural infestations with Rhipicephalus microplus. These studies suggest T-cell mediated immunity, integrity of the dermis, and calcium signaling are important aspects of tick resistance, while innate immune responses may contribute to susceptibility. Thus our present understanding indicates host immunity to ticks is characterized by a complex interplay between host-effector responses and tick evasion strategies.
The tick-host interface is the skin, an organ increasingly recognized to have a significant role in immunity, acting as a sentinel organ that also shapes the ensuing immune response . Anatomically, the skin is divided into two compartments, the epidermis and dermis. The barrier function of the epidermis is maintained by keratinocytes, while keratinocytes, lymphocytes, and langerhans cells play a role responding to epidermal invasion . The dermal compartment is much more heterogeneous, with lymphocytes, macrophages, mast cells, natural killer cells, fibroblasts, and multiple types of dendritic cells . In addition, lymphatic and vascular channels allow the migration of many additional cell types into the dermis. Thus the skin presents a complex array of resident and circulating cells that participate in homeostasis, immunosurveillance, and immune responses. In the case of tick feeding, the cutaneous response represents both the initiation and effector functions of the host. In an effort to understand the spectrum and temporal patterns of the in vivo host response to ticks, we used a PCR-array based approach to characterize the patterns of cutaneous bite-site gene expression during the course of primary and secondary infestations of mice with I. scapularis nymphs.
Pathogen-free I. scapularis colonies were maintained in our laboratory as described . All life cycle stages were kept in sterile glass vials with mesh tops in desiccators at 22°C containing saturated salt solutions to obtain 97% relative humidity with a 16:8 hour photoperiod. For routine colony maintenance adult ticks were fed on New Zealand white rabbits and nymphs and larvae were fed on mice.
Time course infestations
To perform time course infestations, six-week-old female BALB/c mice were placed in individual restrainers and infested with 10-15 pathogen-free I. scapularis nymphs. Ticks were allowed to attach for at least one hour and unattached ticks were discarded. Mice were then removed from restraints and housed individually. Secondary infestations involved two rounds of infestation. Mice were infested with 10-15 I. scapularis nymphs that were allowed to complete their feeding cycle (4-5 days). Fourteen days after the last primary infestation tick dropped off the animals, mice were re-infested with 10-15 I. scapularis nymphs. For tissue harvesting, infested mice were euthanized by CO2 inhalation followed by cervical dislocation and 4 mm punch biopsies were taken from the feeding lesion at 12, 48, 72, and 96 hr post-infestation. Three mice were measured at each time point; controls consisted of 3 similarly housed but tick-free mice. Biopsies were stored in RNAlater (Ambion) at -20°C for RNA and snap frozen in liquid nitrogen and stored at -80°C for cytokine analysis. The Institutional Animal Care and Use Committee of the University of Texas Medical Branch approved all animal experiments (protocol #0907054).
Mouse tissues stored in RNAlater (Ambion) were used for RNA extraction by a combination of Trizol reagent (Invitrogen) and RNeasy (Qiagen) protocols that included an in-column DNase digestion step. Quality and integrity of RNA was verified by the ratio of readings at A260/A280 and A260/A230 (all samples read ≥ 1.8 by both ratios), and by denaturing (formaldehyde) agarose gel electrophoresis followed by staining with Sybr Gold stain (Invitrogen). All samples had readily visible 18S and 20S RNA bands, indicating minimal degradation. Eluted RNA samples were aliquoted and stored at -80°C until use.
Host gene expression profiling using pathway-specific PCR Array analysis
Host cutaneous gene expression was assessed at each time point using three commercially available RT2 Profiler PCR Arrays (SAbiosciences, now Qiagen). Arrays were chosen to measure biological pathways related to T-helper cell differentiation (PAMM-073), wound healing (PAMM-013), and signal transduction (PAMM-014). Each 96-well array contains 84 test and five housekeeping genes. Each array also included controls to assess genomic DNA contamination, RNA quality, and general qRT-PCR performance (for more information, see http://www.sabiosciences.com/). For each array, 1 μg total RNA purified from skin biopsies was converted into cDNA using the RT2 First strand kit (Qiagen). Template cDNAs were mixed with RT2 SYBR Green/Fluorescein qPCR Master Mix (Qiagen) and loaded onto the array using an 8-channel pipette. Arrays were run on an iCycler iQ5 real-time PCR System (Bio-Rad) under standard cycling conditions (10 min at 95°C, 15 s at 95°C, 1 min 60°C for 40 cycles, and an 80-cycle 55-95°C melt curve). The instrument's software was used to calculate the threshold cycle (C t ) values for all molecules analyzed.
Array data analysis
Fold-changes in gene expression between test and control mice were calculated using the ΔΔCt method. For each included gene, individual measurements that were below the threshold selected (C t > 34) were excluded from further analysis. This was done to reduce the impact of stochastic variations in rare transcripts on the calculated fold change and its associated p-value. Data normalization was based on correcting all C t values for the average C t values of the invariant endogenous control genes hypoxanthine guanine phosphoribosyl transferase (Hprt) and heat shock protein 90 alpha (Hsp90ab1). These genes were selected based on analysis of tested housekeeping genes in geNorm . Statistical significance was assessed using LIMMA (linear models in micro-array analysis) in HTqPCR, an R-based program designed for real-time PCR array data analysis . Statistical comparisons were generated for all time points vs. uninfested controls, between time-points (i.e., 12 hr-48 hr, 48 hr-72 hr, etc.), and between infestations (e.g. 12 hr primary vs. 12 hr secondary, etc.). Data sets were filtered with the following criteria: fold change ≥ 3 or ≤ -2 with an adjusted p-value ≤ 0.01. Array data was made publicly available through Gene Expression Omnibus accession number GSE33345.
Gene ontology analysis was conducted on the resulting lists of significantly modulated genes. All significant results from any time point during the primary infestation were divided into three lists: all modulated, upregulated, and downregulated. Similar lists were created for the secondary infestation. Each list was then submitted to the Database for Annotation, Visualization, and Integrated Discovery (DAVID) website using all genes measured as a background list [17, 18]. The functional annotation chart and functional annotation clustering tools were used to assess enriched gene ontology terms; due to the small background list, terms with p-values ≤ 0.05 were considered significant.
Validation of array data
Array results were validated by an additional experiment. Skin biopsies from tick bite-sites were collected as before from two time points during primary (48 and 96 hr p.i.) and secondary (48 and 72 hr p.i.) infestation. Four mice were used at each time point. Twenty-five genes were selected from the list of significantly modulated genes from the array experiment and assayed by additional real-time PCR (see additional file 1). Primer assays and SYBR green master mix were purchased from Qiagen and added to PCR plates to make custom-made arrays. These primer assays contain pre-optimized primer pairs but the primer sequences are proprietary information of Qiagen. Custom-made arrays measured the same 5 house-keeping genes as the original arrays, and included both no template and no first strand controls. In contrast to the arrays, each gene was measured in triplicate. These plates were run and analyzed as the PCR arrays, including the melt curve. To keep data analysis consistent with the PCR arrays, Hprt and Hsp90ab1 were used as normalization genes without additional analysis with geNorm.
The relative concentrations of interleukin-1β (IL-1β), IL-3, IL-4, IL-6, IL-10, IL-17A, interferon-γ (IFN-γ), and monocyte chemoattractant protein-1 (MCP-1 or CCL2) at the tick bite site were quantified using an 8-analyte bioplex assay and the Bioplex 200 system (Bio-Rad). Samples represented two time points (4 mice/time point) in the primary and secondary infestations (48 hr and 96 hr primary, and 48 hr and 72 hr secondary). Biopsies were removed from storage at -80°C and immediately homogenized in 1 mL protein extraction buffer containing 0.5% BSA (Fisher Scientific), 0.1% Igepal-630 (Sigma), and 1% Halt protease inhibitor (Thermo Scientific) in PBS. Homogenates were centrifuged at 20,000 × g at 4°C for 20 min. The resulting supernatants were divided into aliquots and stored at -80°C until use. Standards, blanks, and samples were analyzed in duplicate according to the manufactures' instructions. Analyte concentrations were determined from the standard curve(s) by analysis of mean fluorescent intensity values using the Bio-Plex Manager™ 6.0 software. Individual time points were compared to controls using a 2-tailed T-test in Prism (GraphPad Software, Inc.).
Gene ontology results from genes upregulated during primary infestation.
Chemokine receptor binding
Chemokine signaling pathway
Small chemokine, interleukin-8-like
Cytokine-cytokine receptor interaction
NOD-like receptor signaling pathway
Extracellular region part
Response to wounding
Gene ontology results from genes downregulated during primary infestation.
Positive regulation of biosynthetic process
Positive regulation of cellular biosynthetic process
Positive regulation of gene expression
Positive regulation of macromolecule biosynthetic process
Positive regulation of macromolecule metabolic process
Positive regulation of nitrogen compound metabolic process
Positive regulation of nucleobase, nucleoside, nucleotide and nucleic- acid metabolic process
Positive regulation of RNA metabolic process
Positive regulation of transcription
Positive regulation of transcription from RNA polymerase II- promoter
Positive regulation of transcription, DNA-dependent
Regulation of RNA metabolic process
Regulation of transcription
Regulation of transcription from RNA polymerase II promoter
Regulation of transcription, DNA-dependent
Sequence-specific DNA binding
IPR013568:SEFIR (IL-17R domain)
Gene ontology results from genes upregulated during secondary infestation.
Positive regulation of signal transduction
Positive regulation of cell communication
Four-helical cytokine, core
Regulation of cell activation
Regulation of peptidyl-tyrosine phosphorylation
Regulation of leukocyte proliferation
Regulation of T cell activation
Positive regulation of peptidyl-tyrosine phos.
Cytokine-cytokine receptor interaction
Regulation of lymphocyte activation
Hematopoietic cell lineage
Regulation of leukocyte activation
Regulation of hematopoiesis by cytokines
Regulation of hematopoiesis by cytokines
Regulation of mononuclear cell proliferation
Regulation of lymphocyte proliferation
Regulation of protein kinase cascade
Positive regulation of T cell activation
Regulation of cellular localization
Positive regulation of protein kinase cascade
Regulation of protein modification process
Regulation of cellular protein metabolic process
Positive regulation of immune system process
Regulation of chemokine production
Regulation of chemokine biosynthetic process
Monocyte and its Surface Molecules
Positive regulation of cell activation
Cell adhesion molecules (CAMs)
Positive regulation of protein modification process
External side of plasma membrane
Regulation of T cell proliferation
Regulation of peptidyl-tyrosine phosphorylation
Positive regulation of T cell activation
Positive regulation of peptidyl-tyrosine phos.
Regulation of cellular localization
Monocyte and its Surface Molecules
Regulation of alpha-beta T cell differentiation
Cell adhesion molecules (CAMs)
Positive regulation of alpha-beta T cell diff.
Neutrophil and Its Surface Molecules
Adhesion Molecules on Lymphocyte
Natural killer cell mediated cytotoxicity
In contrast with the primary infestation, upregulated transcripts during secondary exposure were enriched for GO-terms related to a wide array of categories as shown in Table 3. From this data, four major categories are evident: cytokine, chemotaxis, immune cell signaling and activation, and leukocyte adhesion. Other groups of enriched terms included T-cell regulation and cell surface molecules while the remaining terms clustered into a number of minor categories related to putative secreted, hematopoietic lineage, inflammation, protein-protein interactions (sushi domain), activation peptide, and tyrosine kinase phosphorylation. On the other hand, only GO terms 'negative regulation of cell proliferation' (GO_BP:0008285) and SEFIR (IPR013568:SEFIR) were significantly enriched from genes downregulated in the secondary exposure.
Modulation of gene expression during primary infestations
Among downregulated genes, the most notable were members of the IL-17 receptor family (Il17rb, rc, and re), which were consistently downregulated during the entire feeding process while the pro-inflammatory cytokines IL-17a, c, d, and f were not expressed. Other downregulated genes were T-cell molecules CD28, GATA3, and retinoic acid-related orphan nuclear hormone receptor C (RORC), DNA repair molecule telomerase reverse transcriptase (TERT), basement membrane/ECM structural components HAPLN1 and LAMA1, neural adhesion molecule NCAM2, mitogenic pathway member JUN, WNT pathway members JUN and LEF1, hedgehog pathway members BMP2 and 4, and retinoic acid pathway member HOXA1.
Modulation of gene expression during secondary infestations
Array result validation
During tick feeding, the cutaneous environment responds to skin injury by initiating innate defense mechanisms, shaping the ensuing adaptive immune response, and accommodating effector responses of adaptive immunity. In contrast, the feeding tick secretes an arsenal of salivary molecules that pharmacologically inhibit potentially unfavorable host responses . The late initiation of host responses during primary infestation compared to secondary infestation is a striking example of tick-induced suppression of the host response. Early events at the bite site become measurable by 48 hours p.i. and include upregulation of CLEC7a, a lectin pattern recognition receptor. I. scapularis SALP15 has been shown to modulate dendritic cell (DC) function through the lectin receptor DC-SIGN . Together these results suggest lectin pattern-recognition receptors may be important in initiation and modulation of anti-tick immunity. Ligation of CLEC7a induces the up-regulation of CXCL2 and IL-10 , molecules that were also upregulated in our study. Tick-induced expression of IL-10 has been previously reported and may represent a method of immune evasion by dampening pro-inflammatory responses [22, 23].
Other cytokines upregulated early in the host response were IL-1β and IL-6. These are both potent pro-inflammatory cytokines suggesting a balance between anti-inflammatory IL-10 and pro-inflammatory IL-1β and IL-6 during the early host response to ticks. The presence of IL-1β and IL-6 at the bite site is supported by previous studies [24, 25] and the concomitant upregulation of ICAM1, PTGS2 (COX2), CXCL1, CXCL2, CXCL5, and MMP13, molecules known to be induced by these cytokines [26–28]. The CXCL chemokines are potent chemoattractants for neutrophils, although they have also been shown to attract monocytes and mast cells . CCL2 and CCL7 were originally described as macrophage chemotactic proteins (MCP) 1 and 3, reflecting their primary role as chemoattractants for macrophages, but they are also known to recruit basophils, eosinophils, NK cells, and DCs . Recruitment of these cells into the bite site could be facilitated by the upregulation of SELL and ITGB2. These results suggest a model of immune activation during primary infestation where CLEC7a initiates neutrophil chemotaxis and anti-inflammatory cytokine production. Increased production of IL-1β and IL-6 by unknown mechanisms could play a role in promoting upregulation of chemokines specific for neutrophils and macrophages which in turn produce matrix metalloproteinases (MMP) and prostaglandins. Neutrophils are known to be present at the bite-site , but their role in anti-tick immunity is not well understood. Based on the previous identification of I. scapularis salivary proteins (ISL929 and 1373) that reduce super-oxide formation and expression of β-2 integrins in neutrophils treated with TNF-α , it is reasonable to assume they are important components of anti-tick immunity. These changes suggest reduced neutrophil ability to respond to tissue insult and destroy phagocytosed infectious agents.
Matrix metalloproteinases have a wide range of potential functions at the tick bite-site. MMP cleavage of ECM components exposes cryptic sites that have been associated with increased migration of leukocytes to the inflammatory focus; cleavage can also release bioactive molecules from the ECM. I. scapularis has been shown to possess a large family of salivary serine-protease inhibitors (serpins) that may be important in inhibiting host responses . Immunization of rabbits with a serpin from I. ricinus resulted in increased tick mortality and reduced weight and fecundity in female ticks . Since MMPs degrade and inactivate endogenous serpins , it is reasonable to hypothesize that MMPs contribute to host immunity by degrading tick-secreted serpins. MMPs also aid in angiogenesis and wound healing, processes that are inhibited by tick feeding . Gene ontology gives general support to this analysis of the primary infestation. Significant terms from genes upregulated during primary infestation clustered into host response and biomineral formation groups. The host response category was dominated by chemokine, chemotaxis, cytokine, and immune response terms, although none of these terms were specific for any cell type. GO analysis also supported the role of upregulated genes as secreted molecules acting in the extracellular space. Analysis of downregulated genes during primary infestation identified nucleotide metabolism/transcription and SEFIR domain (IPR013568: SEFIR, an IL-17R domain) as significant. These terms are of interest for CD4 T cell differentiation (discussed below), and the possibility that tick feeding suppresses transcription during primary infestation. This is a potential mechanism behind the late induction of host responses during primary infestation.
A second exposure to feeding by I. scapularis nymphs resulted in a faster and stronger host response as shown in Figures 3A and 3B. In contrast to the primary infestation, very significant gene modulation was evident by 12 hrs p.i. The genes modulated during primary infestations were also modulated during secondary infestations and were, in general, the genes with the highest fold-changes. Thus we postulate that genes upregulated during the primary infestation form a core host-response that drives anti-tick immunity even on repeated exposure.
The migration of cells into an inflammatory focus is an important aspect of host immunity . Resident cells must recognize skin damage by the feeding tick and secrete factors that enhance the recruitment of immune effectors to the bite site. Gene ontology analyses of upregulated genes during the secondary infestation strongly support the important role of chemotaxis in the anti-tick immune response. Specific GO terms suggested the migration of neutrophils, monocytes, other leukocytes, and lymphocytes into the bite site (Table 3). The upregulation of CCL1 was the only observed change in chemokine expression between primary and secondary infestation. Interestingly, this chemokine has been shown to attract Th2 and T regulatory cells . Other upregulated genes known to support cell migration included selectins, integrins, and the integrin ligand ICAM1. While a number of alpha chain integrins were upregulated, the only beta chain upregulated was β-2. In support of previous reports that I. scapularis saliva inhibited endothelial cell expression of P selectin , our study showed only minimal upregulation of SELP that was not supported by later validation.
Many additional cytokines were modulated during the secondary infestation when compared to the primary exposure (Figure 3A). These transcripts group together to form the cytokine cluster on gene ontology analysis, lending formal support to their importance in the anti-tick response. In particular, IL-4 and IL-13 were upregulated; these cytokines can be produced by Th2 cells, but also by basophils, eosinophils, and mast cells. Basophils have been shown to be indispensible for anti-tick immunity in models of infestation where acquired resistance occurs [36, 37], and their migration into the bite site was supported by the upregulation of CCL chemokines and IL-3, which are chemotactic factors for basophils . It is interesting to note the significant increase in expression of IL-4 and IL-13 during the secondary exposure despite the down regulation of IL-25, an important inducer of type-2 immunity . In contrast to these type-2 cytokines, upregulation of IFN-γ and IL-27 could be due to the presence of T cells, Th1 cells, NK cells, and antigen presenting cells [40, 41]. The upregulation of IFN-γ is surprising in light of previous reports of suppression by tick saliva , although negligible increases in expression have been previously reported in BALB/c mice infested with I. scapularis . The mechanisms behind IFN-γ and IL-4 upregulation were strong enough to overcome the downregulation of IL-18, a known inducer of both cytokines . Upregulation of colony stimulating factors 2 (GM-CSF or CSF2) and 3 (G-CSF or CSF3) and IL-3 suggests tick feeding may stimulate increased hematopoiesis and/or myelopoiesis. This possibility was supported by the gene ontology analysis (Table 3), previous reports of extramedullary erythropoiesis in tick-infested mice , and the downregulation of IL-17d, an inhibitor of hematopoietic progenitor colony formation . Finally, our study also supports previously reported repression in the expression of tumor necrosis factor family members by tick salivary molecules [6, 44]. In summary, the cytokine profile during secondary infestation presents a complex interplay between inducers and repressors of type 1 and type 2 immunity.
Th2 polarization of the cytokine response to tick feeding has been thoroughly documented by in vitro and in vivo studies . For this reason, we sought to characterize the modulation of genes associated with T-cell and helper T-cell differentiation. During primary infestation, classic T-cell markers such as CD3, CD4, and CD8 did not significantly change, suggesting early T-cell involvement is minimal. Interestingly, the expression of co-stimulatory molecule CD28 was downregulated, which could be due to a lack of CD4 T-cell activation at the bite site, or the migration of CD28-expressing cells out of the skin. Genes related to Th17 differentiation, including the transcription factor RORC, IL-17, and the IL-17 receptors were either unchanged or downregulated, despite the high levels of IL-1β and IL-6. Most genes related to Th2 development were unchanged with the exception of GATA3, which was downregulated. GATA3 is an important transcription factor in Th2 development. Transcripts related to Th1 and T reg development were unchanged. These results suggest that during primary infestation of mice with I. scapularis nymphs, the cutaneous environment is not strongly polarized (or polarizing) toward any helper T-cell sub-set.
On secondary infestation, the upregulation of T-cell markers CD2, CD3, CD4, and CD8 suggested T-cell involvement at the bite-site. However, the polarization of CD4 T-cells remained equivocal. While Th2 cytokines IL-4 and IL-13 were upregulated, GATA3 remained repressed (Figure 3A). Similarly, the Th1 cytokine IFN-γ was upregulated, but IL-12 and TBX21 remained unresponsive. Th17-related transcripts were downregulated or unchanged. Interestingly, FOXP3 and IL-10 were upregulated, supporting a possible role for T regulatory cells at the bite site. In summary, results from the secondary exposure strongly suggests Th17 involvement at the bite-site is unlikely, while the remaining data shows a mixed Th1/Th2 cytokine profile and suggests the involvement of T regs. Failure to produce a polarized CD4 T cell response was also observed when keyhole limpet haemocyanin (KLH)-specific T cells were stimulated with KLH-loaded DCs in the presence of Rhipicephalus sanguineus tick saliva . This implies that non-polarized CD4 T cell responses may be a common trait of anti-tick immunity and also supports our results at the protein/cellular level. Sialostatin L, an I. scapularis salivary protein, suppressed IL-17 production by lymph node cells during the induction of experimental autoimmune encephalomyelitis in mice . In our results, significant Th17 suppression was observed (Figure 2, 3A) even from a naïve state, supporting the possibility that tick saliva contains potent suppressors of Th17 immunity.
Another focus of the present study was to uncover novel signaling pathways activated at the tick bite-site. Surprisingly, most genes related to the signaling pathways tested were either downregulated or unresponsive. Immunoreceptor signaling was a significant exception. Gene ontology results showed the largest gene cluster was related to immune cell signaling and activation. This is consistent with the rest of our results and suggests immunoreceptor signaling as a potential major pathway induced by tick feeding. However, we were unable to show any modulation of signal transducers and activators of transcription (STAT) or NFκB pathway molecules. The lack of STAT modulation in our study was surprising since STAT molecules are important effector molecules of cytokine signaling that induce their own expression . Modulation or silencing of the NFκB pathway could be significant because of its vital role in the induction and regulation of immunity . These results paired with the increase of SOCS transcripts suggest the tick bite site is characterized by both suppression and activation of immunoreceptor signaling.
Gene ontology analysis of the downregulated genes during secondary infestation showed only two significant terms: negative regulation of cell proliferation (GO_BP: 0008285) and SEFIR (IPR013568). This suggests the genes downregulated during secondary infestation do not fit into a common theme for GO enrichment. However, many groups and pathways were qualitatively downregulated. Transcripts related to integrin binding and neural adhesion were downregulated, indicating potential impairment of cell migration or adhesion-related signaling at the tick bite site. Increased cell turnover was suggested by the downregulation of genes encoding anti-apoptotic and DNA repair molecules. The downregulation of BM/ECM structure and ECM protease inhibitor groups combined with the upregulation of ECM proteases suggest significant modulation of ECM components. In addition to these groups, genes in classical pathways such as mitogenic, WNT, hedgehog, stress, and metabolism were downregulated. The WNT signaling pathway regulates a variety of cellular processes including cell proliferation, migration, and tissue morphogenesis. In canonical signaling, WNT stabilizes β-catenin that acts as a transcriptional co-activator by interacting with Lef/T-cell transcription factors (Tcfs) to regulate WNT target gene expression. Non-canonical signaling, on the other hand, is calcium-dependent and leads to activation of c-jun N-terminal kinase which plays a role in cell proliferation, differentiation, and apoptosis . In addition to its role in developmental biology, the Hedgehog pathway has been shown to play a role in regulating regenerating cell populations . Since cell proliferation, regeneration, and morphogenesis are involved in wound healing, epithelial maintenance, and hair follicle cycling, tick feeding may influence these processes. However, it is unclear whether this is a result of tick-saliva induced repression or a consequence of the inflammatory process at the bite-site lesion. In this regard, our infestation protocol avoided the use of capsules or any device to restrain the ticks during feeding that might have influenced the inflammatory reaction. In either case, our results qualitatively suggest the tick bite-site is characterized in part by the suppression of signaling molecule transcription.
Our study supports a model of the tick-host interface where tick saliva inhibits gene transcription, Th17 immunity, and signal transduction molecule upregulation. In contrast, the host senses infestation through lectin PRRs and is primarily focused on the recruitment and subsequent activation of immune cells. During primary infestation, neutrophils and macrophages are recruited, while many additional cell types are recruited during secondary infestation. Host effector responses include a mixed Th1/Th2 CD4 T cell response, innate effector functions, a highly proteolytic environment, and increased cell turn-over. These responses are dampened by the action of T regulatory cells, SOCS, and IL-10.
To our knowledge, this is the first report of in vivo transcriptome profiling at the I. scapularis tick-host interface. Our results suggest tick feeding may activate favorable host responses such as the inhibition of gene transcription, downregulation of signaling molecules, and upregulation of inhibitors of inflammation while repressing unfavorable responses such as Th17 immunity. The mixed Th1/Th2 profile reported here is a novel finding that implies greater complexity to the host cutaneous response than previously reported. We believe this study will allow the rational design of further work probing in vivo mechanisms at the tick-host-pathogen interface.
The authors would like to acknowledge Nicole Hausser, Natacha Milhano, and Tais Saito for their assistance in rearing ticks. The project described was supported by the National Institute of Allergy and Infectious Diseases awards number T32AI007526 to DMH and RO1AI062735 to SKW. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Allergy and Infectious Diseases or the National Institutes of Health.
- Brossard M, Wikel S: Tick Immunobiology. Ticks: Biology, Disease and Control. Edited by: Bowman AS, Nuttall PA. 2008, Cambridge, UK: Cambridge University Press, 186-204.View ArticleGoogle Scholar
- Wikel SK, Ramachandra RN, Bergman DK, Burkot TR, Piesman J: Infestation with pathogen-free nymphs of the tick Ixodes scapularis induces host resistance to transmission of Borrelia burgdorferi by ticks. Infect Immun. 1997, 65: 335-338.PubMed CentralPubMedGoogle Scholar
- Trager W: Acquired immunity to ticks. Journal of Parasitology. 1939, 25: 57-81. 10.2307/3272160.View ArticleGoogle Scholar
- Muller-Doblies UU, Maxwell SS, Boppana VD, Mihalyo MA, McSorley SJ, Vella AT, Adler AJ, Wikel SK: Feeding by the tick, Ixodes scapularis, causes CD4(+) T cells responding to cognate antigen to develop the capacity to express IL-4. Parasite Immunol. 2007, 29: 485-499. 10.1111/j.1365-3024.2007.00966.x.View ArticlePubMedGoogle Scholar
- Anguita J, Ramamoorthi N, Hovius JW, Das S, Thomas V, Persinski R, Conze D, Askenase PW, Rincon M, Kantor FS, Fikrig E: Salp15, an ixodes scapularis salivary protein, inhibits CD4(+) T cell activation. Immunity. 2002, 16: 849-859. 10.1016/S1074-7613(02)00325-4.View ArticlePubMedGoogle Scholar
- Hovius JW, de Jong MA, den Dunnen J, Litjens M, Fikrig E, van der Poll T, Gringhuis SI, Geijtenbeek TB: Salp15 binding to DC-SIGN inhibits cytokine expression by impairing both nucleosome remodeling and mRNA stabilization. PLoS Pathog. 2008, 4: e31-10.1371/journal.ppat.0040031.PubMed CentralView ArticlePubMedGoogle Scholar
- Guo X, Booth CJ, Paley MA, Wang X, DePonte K, Fikrig E, Narasimhan S, Montgomery RR: Inhibition of neutrophil function by two tick salivary proteins. Infect Immun. 2009, 77: 2320-2329. 10.1128/IAI.01507-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Oliveira CJ, Cavassani KA, More DD, Garlet GP, Aliberti JC, Silva JS, Ferreira BR: Tick saliva inhibits the chemotactic function of MIP-1alpha and selectively impairs chemotaxis of immature dendritic cells by down-regulating cell-surface CCR5. Int J Parasitol. 2008, 38: 705-716. 10.1016/j.ijpara.2007.10.006.View ArticlePubMedGoogle Scholar
- Wikel SK: Host immunity to ticks. Annu Rev Entomol. 1996, 41: 1-22. 10.1146/annurev.en.41.010196.000245.View ArticlePubMedGoogle Scholar
- Kongsuwan K, Piper EK, Bagnall NH, Ryan K, Moolhuijzen P, Bellgard M, Lew A, Jackson L, Jonsson NN: Identification of genes involved with tick infestation in Bos taurus and Bos indicus. Dev Biol (Basel). 2008, 132: 77-88.Google Scholar
- Piper EK, Jackson LA, Bagnall NH, Kongsuwan KK, Lew AE, Jonsson NN: Gene expression in the skin of Bos taurus and Bos indicus cattle infested with the cattle tick, Rhipicephalus (Boophilus) microplus. Vet Immunol Immunopathol. 2008, 126: 110-119. 10.1016/j.vetimm.2008.06.011.View ArticlePubMedGoogle Scholar
- Piper EK, Jonsson NN, Gondro C, Lew-Tabor AE, Moolhuijzen P, Vance ME, Jackson LA: Immunological profiles of Bos taurus and Bos indicus cattle infested with the cattle tick, Rhipicephalus (Boophilus) microplus. Clin Vaccine Immunol. 2009, 16: 1074-1086. 10.1128/CVI.00157-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Nestle FO, Di Meglio P, Qin JZ, Nickoloff BJ: Skin immune sentinels in health and disease. Nat Rev Immunol. 2009, 9: 679-691.PubMed CentralPubMedGoogle Scholar
- Bouchard KR, Wikel SK: Care, maintenance, and experimental infestation of ticks in the laboratory setting. Biology of Disease Vectors. Edited by: Marquart WC. 2005, San Diego: Elsevier Academic Press, 705-711. SecondGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-PubMed CentralView ArticlePubMedGoogle Scholar
- Dvinge H, Bertone P: HTqPCR: high-throughput analysis and visualization of quantitative real-time PCR data in R. Bioinformatics. 2009, 25: 3325-3326. 10.1093/bioinformatics/btp578.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticlePubMedGoogle Scholar
- Maitra A, Shen F, Hanel W, Mossman K, Tocker J, Swart D, Gaffen SL: Distinct functional motifs within the IL-17 receptor regulate signal transduction and target gene expression. Proc Natl Acad Sci USA. 2007, 104: 7506-7511. 10.1073/pnas.0611589104.PubMed CentralView ArticlePubMedGoogle Scholar
- Torcia MG, Santarlasci V, Cosmi L, Clemente A, Maggi L, Mangano VD, Verra F, Bancone G, Nebie I, Sirima BS: Functional deficit of T regulatory cells in Fulani, an ethnic group with low susceptibility to Plasmodium falciparum malaria. Proc Natl Acad Sci USA. 2008, 105: 646-651. 10.1073/pnas.0709969105.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown GD: Dectin-1: a signalling non-TLR pattern-recognition receptor. Nat Rev Immunol. 2006, 6: 33-43. 10.1038/nri1745.View ArticlePubMedGoogle Scholar
- Kopecky J, Kuthejlova M, Pechova J: Salivary gland extract from Ixodes ricinus ticks inhibits production of interferon-gamma by the upregulation of interleukin-10. Parasite Immunol. 1999, 21: 351-356. 10.1046/j.1365-3024.1999.00233.x.View ArticlePubMedGoogle Scholar
- Schoeler GB, Manweiler SA, Wikel SK: Ixodes scapularis: effects of repeated infestations with pathogen-free nymphs on macrophage and T lymphocyte cytokine responses of BALB/c and C3H/HeN mice. Exp Parasitol. 1999, 92: 239-248. 10.1006/expr.1999.4426.View ArticlePubMedGoogle Scholar
- Kovar L, Kopecky J, Rihova B: Salivary gland extract from Ixodes ricinus tick polarizes the cytokine profile toward Th2 and suppresses proliferation of T lymphocytes in human PBMC culture. J Parasitol. 2001, 87: 1342-1348.View ArticlePubMedGoogle Scholar
- Schoeler GB, Manweiler SA, Wikel SK: Cytokine responses of C3H/HeN mice infested with Ixodes scapularis or Ixodes pacificus nymphs. Parasite Immunol. 2000, 22: 31-40.View ArticlePubMedGoogle Scholar
- Dinarello CA: Immunological and inflammatory functions of the interleukin-1 family. Annu Rev Immunol. 2009, 27: 519-550. 10.1146/annurev.immunol.021908.132612.View ArticlePubMedGoogle Scholar
- Leeman MF, Curran S, Murray GI: The structure, regulation, and function of human matrix metalloproteinase-13. Crit Rev Biochem Mol Biol. 2002, 37: 149-166. 10.1080/10409230290771483.View ArticlePubMedGoogle Scholar
- Strieter RM, Burdick MD, Gomperts BN, Belperio JA, Keane MP: CXC chemokines in angiogenesis. Cytokine Growth Factor Rev. 2005, 16: 593-609. 10.1016/j.cytogfr.2005.04.007.View ArticlePubMedGoogle Scholar
- Viola A, Luster AD: Chemokines and their receptors: drug targets in immunity and inflammation. Annu Rev Pharmacol Toxicol. 2008, 48: 171-197. 10.1146/annurev.pharmtox.48.121806.154841.View ArticlePubMedGoogle Scholar
- Krause PJ, Grant-Kels JM, Tahan SR, Dardick KR, Alarcon-Chaidez F, Bouchard K, Visini C, Deriso C, Foppa IM, Wikel S: Dermatologic changes induced by repeated Ixodes scapularis bites and implications for prevention of tick-borne infection. Vector Borne Zoonotic Dis. 2009, 9: 603-610. 10.1089/vbz.2008.0091.PubMed CentralView ArticlePubMedGoogle Scholar
- Mulenga A, Khumthong R, Chalaire KC: Ixodes scapularis tick serine proteinase inhibitor (serpin) gene family; annotation and transcriptional analysis. BMC Genomics. 2009, 10: 217-10.1186/1471-2164-10-217.PubMed CentralView ArticlePubMedGoogle Scholar
- Prevot PP, Couvreur B, Denis V, Brossard M, Vanhamme L, Godfroid E: Protective immunity against Ixodes ricinus induced by a salivary serpin. Vaccine. 2007, 25: 3284-3292. 10.1016/j.vaccine.2007.01.008.View ArticlePubMedGoogle Scholar
- Sternlicht MD, Werb Z: How matrix metalloproteinases regulate cell behavior. Annu Rev Cell Dev Biol. 2001, 17: 463-516. 10.1146/annurev.cellbio.17.1.463.PubMed CentralView ArticlePubMedGoogle Scholar
- Francischetti IM, Mather TN, Ribeiro JM: Tick saliva is a potent inhibitor of endothelial cell proliferation and angiogenesis. Thromb Haemost. 2005, 94: 167-174.PubMed CentralPubMedGoogle Scholar
- Maxwell SS, Stoklasek TA, Dash Y, Macaluso KR, Wikel SK: Tick modulation of the in-vitro expression of adhesion molecules by skin-derived endothelial cells. Ann Trop Med Parasitol. 2005, 99: 661-672. 10.1179/136485905X51490.View ArticlePubMedGoogle Scholar
- Brown SJ, Galli SJ, Gleich GJ, Askenase PW: Ablation of immunity to Amblyomma americanum by anti-basophil serum: cooperation between basophils and eosinophils in expression of immunity to ectoparasites (ticks) in guinea pigs. J Immunol. 1982, 129: 790-796.PubMedGoogle Scholar
- Wada T, Ishiwata K, Koseki H, Ishikura T, Ugajin T, Ohnuma N, Obata K, Ishikawa R, Mukai K, Yoshikawa S, Kawano Y, Minegishi Y, Yokozeki H, Watanabe N, Karasuyama H: Selective ablation of basophils in mice reveals their nonredundant role in acquired immunity against ticks. J Clin Invest. 2010, 120: 2867-2875. 10.1172/JCI42680.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim S, Prout M, Ramshaw H, Lopez AF, LeGros G, Min B: Cutting edge: basophils are transiently recruited into the draining lymph nodes during helminth infection via IL-3, but infection-induced Th2 immunity can develop without basophil lymph node recruitment or IL-3. J Immunol. 2010, 184: 1143-1147. 10.4049/jimmunol.0902447.PubMed CentralView ArticlePubMedGoogle Scholar
- Barlow JL, McKenzie AN: IL-25: a key requirement for the regulation of type-2 immunity. Biofactors. 2009, 35: 178-182. 10.1002/biof.24.View ArticlePubMedGoogle Scholar
- Xu M, Mizoguchi I, Morishima N, Chiba Y, Mizuguchi J, Yoshimoto T: Regulation of antitumor immune responses by the IL-12 family cytokines, IL-12, IL-23, and IL-27. Clin Dev Immunol. 2010, 2010:Google Scholar
- McLaren JE, Ramji DP: Interferon gamma: a master regulator of atherosclerosis. Cytokine Growth Factor Rev. 2009, 20: 125-135. 10.1016/j.cytogfr.2008.11.003.View ArticlePubMedGoogle Scholar
- Dash Y, Maxwell SS, Rajan TV, Wikel SK: Murine extramedullary erythropoiesis induced by tick infestation. Ann Trop Med Parasitol. 2005, 99: 518-531. 10.1179/136485905X51319.View ArticlePubMedGoogle Scholar
- Starnes T, Broxmeyer HE, Robertson MJ, Hromas R: Cutting edge: IL-17D, a novel member of the IL-17 family, stimulates cytokine production and inhibits hemopoiesis. J Immunol. 2002, 169: 642-646.View ArticlePubMedGoogle Scholar
- Ramachandra RN, Wikel SK: Modulation of host-immune responses by ticks (Acari: Ixodidae): effect of salivary gland extracts on host macrophages and lymphocyte cytokine production. J Med Entomol. 1992, 29: 818-826.View ArticlePubMedGoogle Scholar
- Sa-Nunes A, Bafica A, Antonelli LR, Choi EY, Francischetti IM, Andersen JF, Shi GP, Chavakis T, Ribeiro JM, Kotsyfakis M: The immunomodulatory action of sialostatin L on dendritic cells reveals its potential to interfere with autoimmunity. J Immunol. 2009, 182: 7422-7429. 10.4049/jimmunol.0900075.PubMed CentralView ArticlePubMedGoogle Scholar
- Li WX: Canonical and non-canonical JAK-STAT signaling. Trends Cell Biol. 2008, 18: 545-551. 10.1016/j.tcb.2008.08.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayden MS, Ghosh S: Shared principles in NF-kappaB signaling. Cell. 2008, 132: 344-362. 10.1016/j.cell.2008.01.020.View ArticlePubMedGoogle Scholar
- Bikkavilli RK, Malbon CC: Mitogen-activated protein kinases and Wnt/beta-catenin signaling: Molecular conversations among signaling pathways. Commun Integr Biol. 2009, 2: 46-49. 10.4161/cib.2.1.7503.PubMed CentralView ArticlePubMedGoogle Scholar
- Ingham PW, Nakano Y, Seger C: Mechanisms and functions of Hedgehog signalling across the metazoa. Nat Rev Genet. 2011, 12: 393-406. 10.1038/nrg2984.View ArticlePubMedGoogle 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.