Figures
Abstract
Ixodes species ticks are competent vectors of tick-borne viruses including tick-borne encephalitis and Powassan encephalitis. Tick saliva has been shown to facilitate and enhance viral infection. This likely occurs by saliva-mediated modulation of host responses into patterns favorable for viral infection and dissemination. Because of the rapid kinetics of tick-borne viral transmission, this modulation must occur as early as tick attachment and initiation of feeding. In this study, cutaneous bite-site lesions were analyzed using Affymetrix mouse genome 430A 2.0 arrays and histopathology at 1, 3, 6, and 12 hours after uninfected Ixodes scapularis nymphal tick attachment. At 1 and 3 hrs after attachment, the gene expression profile is markedly different than at later time points. Upregulated gene ontology term clusters enriched at 1 and 3 hrs were related to post-translational modification. At 6 and 12 hrs, cytoskeletal rearrangements, DNA replication/cell division, inflammation, and chemotaxis were prominent clusters. At 6 and 12 hrs, extracellular matrix, signaling, and DNA binding clusters were downregulated. Histopathological analysis shows minimal inflammation at 1 and 3 hrs but an appreciable neutrophil infiltrate at 6 and 12 hrs. In addition, putative hyperemia, localized necrosis, and increased ECM deposition were identified. Putting the gene expression and histopathology analysis together suggests early tick feeding is characterized by modulation of host responses in resident cells that merges into a nascent, neutrophil-driven immune response by 12 hrs post-attachment.
Citation: Heinze DM, Carmical JR, Aronson JF, Thangamani S (2012) Early Immunologic Events at the Tick-Host Interface. PLoS ONE 7(10): e47301. https://doi.org/10.1371/journal.pone.0047301
Editor: Bradley S. Schneider, Metabiota, United States of America
Received: July 6, 2012; Accepted: September 10, 2012; Published: October 15, 2012
Copyright: © Heinze et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: DH was supported by the National Institute of Allergy and Infectious Diseases award number T32AI007526. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Emerging and re-emerging diseases transmitted by blood feeding arthropods are significant global public health problems. Ticks transmit the greatest variety of pathogenic spirochetes, rickettsiae and viruses of any blood feeding arthropod [1]. Infectious agents transmitted by ticks are delivered to the vertebrate host together with saliva at the bite site. Tick salivary glands produce a complex repertoire of bioactive molecules that creates an immunologically privileged microenvironment facilitating blood feeding and pathogen transmission [2]. Ticks remain attached to their hosts for a few hours in the case of soft ticks or several days in the case of hard ticks. To successfully complete their feeding, ticks have evolved strategies to circumvent innate immune responses when feeding on naïve hosts and both innate and adaptive immune responses when feeding on tick-experienced animals. Skin is the interface between the tick and the host. Skin acts as a physical barrier and also contains an array of resident immune cells such as eosinophils, mast cells, dendritic cells, macrophages and keratinocytes. Tick feeding and salivary gland molecules skew host immune response away from TH1 and toward a TH2 profile [3]. Salivary gland proteins and extracts have been shown to inhibit immune cells at the bite site and dendritic cell maturation and migration [4].
During tick feeding and attachment, significant morphological changes occur in salivary glands [5]. Salivation is not a continuous process during blood feeding [6] and the repertoire of salivary proteins changes during the course of feeding [7], [8]. These temporal patterns allow the saliva proteins of the tick to “prime” the feeding site to different degrees prior to introduction of infectious agents. Not all tick-borne pathogens are transmitted at the same time during the feeding process. Tick-borne encephalitis was observed to be transmitted within the first few hours of attachment/feeding, while Borrelia burgdorferi was observed to be transmitted between 24 and 48 hours post tick attachment/feeding [9], [10], [11]. We believe that the temporal expression of immunomodulatory tick salivary proteins secreted into the bite site depends on the pathogen it transmits. In this study, we sought to characterize tick-induced changes in cutaneous gene expression at the earliest stages of attachment/feeding by I. scapularis nymphs using Mouse Genome microarrays. This will allow us to understand immunomodulation at the tick-host interface induced by tick saliva that facilitates tick-borne virus transmission.
Methods
Ethics Statement
All experiments were conducted in an arthropod containment level 2 (ACL-2) facility in accordance with an animal use protocol approved by the University of Texas Medical Branch (UTMB) Institutional Animal Care and Use Committee (IACUC).
Animals
BALB/c mice used in this study were obtained from The Jackson Laboratory (Bar Harbor, ME). Mice were cared for in accordance with guidelines of the Committee on Care and Use of Laboratory Animals (Institute of Laboratory Animal Resources National Research Council, Washington, DC).
Tick Colony Maintenance and Infestations
Ixodes scapularis ticks were maintained as described in [12] and [13]. Eight-week old BALB/c mice were anesthetized with a 150 υL intraperitoneal injection containing 10 mg/mL ketamine (Fort Dodge Animal Health, Fort Dodge, IA) and 1 mg/mL xylazine (Phoenix Pharmaceutical, St. Joseph, MO) in PBS (Gibco, Life Technologies, Carlsbad, CA). Eight Ixodes scapularis nymphs were placed on the ears of each mouse and observed until attachment. Mice were euthanized by an additional injection of 200 υL anesthetic followed by cervical dislocation at 1, 3, 6, and 12 hours after tick attachment. Skin biopsies with attached ticks were harvested using 4 mm biopsy punches (Premier Products Co., Plymouth Meeting, PA). Control samples were harvested in the same way from tick-free mice. For the micro-array experiment, 3 mice were used at each time point; for the real-time PCR validation experiment, 4 mice were used at each time point. Skin biopsies were preserved in RNALater (Ambion, Life Technologies, Carlsbad, CA) or 10% neutral buffered formalin. The Institutional Animal Care and Use Committee of the University of Texas Medical Branch approved all animal protocols.
RNA Extraction
Ticks were removed from all skin biopsies before RNA extraction. Tissue samples were homogenized individually in 1 mL Trizol (Life Technologies, Carlsbad, CA) using an Ultra-Turrax T8 (Ika, Wilmington, NC) tissue disperser. After 5 min incubation, 200υL chloroform (Fisher Scientific, Waltham, MA) was added and the samples incubated for 3 min. Samples were centrifuged at 12,000×g, 4°C for 15 min, and the aqueous phase was retained. One volume 70% ethanol (Acros Organics) was added and the samples applied to RNeasy micro kit (Qiagen) columns. The RNeasy protocol was then followed, including the in-column DNase digestion step. After extraction, RNA was quantitated spectrophotometrically using a NanoDrop ND-1000 (NanoDrop Technologies, DE). All samples were required to read greater than 1.8 on both A260/A280 and A260/A230 ratios. For subsequent microarray analysis, quality of the purified RNA was assessed by visualization of 18 S and 28 S RNA bands using an Agilent BioAnalyzer 2100 (Agilent Technologies, CA). Resulting electropherograms were used in the calculation of the 28S/18S ratio and the RNA Integrity Number, which was greater than 6.8 in all samples [14]. For subsequent real-time PCR analysis, RNA integrity was determined by denaturing (formaldehyde) agarose gel electrophoresis followed by staining with Sybr Gold stain (Invitrogen). Visualization of clear ribosomal bands indicated minimal degradation. Eluted RNA samples were aliquoted and stored at −80°C until use.
Affymetrix GeneChip Gene Expression Analysis
Total RNA (500 ng) was converted to cRNA for microarray analysis using the Ambion MessageAmp™ Premier RNA Amplification Kit (Life Technologies Corporation, CA) according to manufacturer's instructions. Total fragmented cRNA (10υg) was hybridized to the Affymetrix GeneChip Mouse Genome 430A 2.0 array according to the manufacturer’s (Affymetrix, CA) conditions. The chips were washed and stained in a GeneChip Fluidics Station 450 and fluorescence detected with an Affymetrix-7G Gene Array scanner using the Affymetrix GeneChip Command Console software (AGCC1.1). Raw data can be accessed through Gene Expression Omnibus record GSE39100.
Gene Expression Data Analysis
Gene expression changes in comparison to tick-free mice were identified using Partek Genomics Suite (Partek, MO) following the default gene expression workflow. The resulting values were then filtered for p-values ≤0.05 and a fold change ≤ −1.5 or ≥ +1.5. Lists of up and downregulated genes at each time point were individually submitted to the Database for Annotation, Visualization, and Integrated Discovery (DAVID) [15], [16] website using the Mouse Genome 430A 2.0 array as a background list. The functional annotation clustering tool was used to cluster gene ontology terms with shared genes into groups to allow an easier functional understanding of the array data. Gene expression data was also entered into ingenuity pathway analysis software (Ingenuity Systems, Redwood City, CA).
Real-time PCR Validation
A list of 28 genes (Table 1) were chosen for further real-time PCR validation based on significant fold change in the array data at 12 hpi or previously identified genes of interest in host anti-tick responses [13]. Pre-optimized primer pairs for these genes were purchased from SABiosciences (Qiagen); primer sequences are property of Qiagen. Primers were mixed with RT2 SYBR green qPCR master mix (Qiagen) and aliquoted into 20 iCycler iQ PCR plates (Bio-Rad) using an epMotion 5075 automated pipetting system (Eppendorf). Plates were sealed and stored at −20°C until use. For each real-time PCR run, the RT2 First Strand Kit (Qiagen) was used to convert 1 µg total RNA into cDNA, which was then loaded onto PCR plates using the epMotion 5075 automated pipetting system (Eppendorf). These plates were run on an iCycler iQ5 real-time PCR instrument (Bio-Rad) with the following cycling protocol: 10 min at 95°C; 15 s at 95°C, 1 min 60°C for 40 cycles, and an 80-cycle (+0.5°C/cycle) 55–95°C melt curve. Every run included hypoxanthine guanine phosphoribosyl transferase (Hprt) and heat shock protein 90 alpha (Hsp90ab1) as endogenous control genes, and ‘no template’ and ‘no first strand’ controls. HTqPCR, an R-based program for real-time PCR data analysis [17], was used to analyze data using the delta-delta Ct method for gene expression normalization and measurement, and the linear models in microarray analysis (LIMMA) package for statistical comparisons between infested and tick-free mice, as previously described [13].
Histological Analysis
Skin biopsies were fixed for a minimum of 48 hours in 10% neutral buffered formalin, treated for 2 hours with decal (Decal Chemical Corp, Tallman, NY), and embedded in paraffin with special attention to orient the tissue sample to allow a cross-section of the skin and a longitudinal section of the tick upon sectioning. Paraffin blocks were carefully sectioned and the slides containing the entry of the hypostome into the epidermis were retained for H&E and Giemsa staining.
Results and Discussion
Microarray Analysis of Host Immune Responses to Early Tick Feeding
The immune response at the tick-host interface was investigated at 1, 3, 6 and 12 hours post nymphal tick infestation (hpi) using microarrays. Significantly modulated genes increased across time (Figure 1a), reflecting the development of host responses as the infestation progressed. A higher percentage of modulated genes were shared with adjacent than distant time points (Figure 1b). The specific gene expression profiles were similar between 1 and 3 hours, but showed appreciable change between 3 and 6 hours, and again between 6 and 12 hours post-infestation (Figure 1c). These results suggest cutaneous responses undergo rapid changes in gene expression profiles in the early stages of tick feeding.
The immune response at the tick-host interface was investigated at 1, 3, 6 and 12 hours post nymphal tick infestation (hpi) using mouse Affymetrix GeneChip microarrays. A: Number of significantly up and downregulated genes measured at each time point during tick infestations of mice with I. scapularis nymphs compared to tick-free mice; B: Venn diagram showing overlap of significantly modulated genes between time points; C: Differential gene expression data was used to generate a heat map using Partek Genomics analysis suite showing temporal changes in gene expression profiles.
Significantly modulated genes were divided into up and downregulated lists at each time point and submitted to the Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics database. The functional annotation-clustering tool was used to group similar terms together. These clusters were then named according to the gene ontology terms in each cluster (Table 2). At 1 and 3 hrs post-infestation, the only significantly up-regulated cluster was “post-translational modification: ubiquitin, isopeptide;” no significantly downregulated clusters were observed. Changes in gene expression at 1 and 3 hrs post-infestation were related to signaling pathways such as NFκB and cation homeostatsis, suggesting pro-inflammatory pathways are already activated. At 6 hrs post-infestation, clusters related to cytoskeletal elements, keratinocyte migration, cell signaling, transcription, and immune responses were prominent. At 12 hrs post-infestation, cell cycle, cytoskeletal elements, and immune response clusters were observed. Immune response clusters differed between 6 and 12 hrs post-infestation. At 6 hrs, anti-microbial responses, immunoreceptor signaling, and negative regulation of myeloid differentiation were significant, while at 12 hrs, inflammation and chemotaxis were significant. These results suggest a rapid, pro-inflammatory evolution of the early host response beginning soon after attachment that progress from early inflammatory signaling and pre-programmed anti-microbial responses to the infiltration of innate immune cells by 12 hpi.
Cytoskeletal Changes
At both 6 at 12 hpi, the most significantly upregulated gene ontology clusters were related to components of the cytoskeleton such as intermediate filaments. A closer look at these genes revealed many keratin intermediate filament transcripts. Keratin intermediate filaments have been shown to protect epithelial tissues from mechanical and non-mechanical stresses, modulate apoptosis, regulate some aspects of skin pigmentation, and control keratinocyte migration during the process of wound healing [18], [19], [20], [21]. Because the initiation of the feeding lesion necessitates significant local damage to epithelial tissues, we believe these ontology terms likely reveal early epithelial attempts to close the wound. Interestingly, Krt6, a gene upregulated at both 6 and 12 hpi has already been shown to be upregulated in keratinocytes at wound margins as early as 6 hrs post injury, and to pair with Krt16/17 during wound healing [21]. The upregulation of all these genes at 12 hpi in our study supports the possibility these gene ontology clusters describe a nascent wound healing response at the tick bite site.
Transcription Factors and Cell Signaling Pathways
A number of gene ontology clusters were related to signaling pathways and downstream targets such as transcription factors. Upregulated clusters at 6 hpi were transcription/transcription factors and protein phosphatase, while downregulated clusters included G-protein signaling domain. While many transcription factors were modulated, the largest group belonged to members of the activator protein-1 (AP-1) transcription factor family. AP-1 is composed of homo- or hetero- dimers of Fos, Jun, Atf3, and Maff proteins. It can activate an extremely broad range of processes that include oncogenesis, regulation of bone homeostasis, and immunity. In the skin, loss of JunB leads to a systemic lupus erythematosis-like phenotype, while loss of both c-Jun and JunB leads to postnatal death. An inducible deletion of c-Jun and JunB in the skin leads to a psoriasis-like phenotype characterized by joint inflammation, hyperkeratosis, and epidermal infiltrates [22]. Thus AP-1 members play a significant role in regulating epidermal inflammatory responses. Our data suggests AP-1 may be activated through the mitogen-activated protein kinase (MAPK) pathway, based on the upregulation of Mak3k6 and negative regulators of the MAPK pathway such as dual-specificity phosphatases (DUSP). In addition to AP-1, members of the NFκB family such as Nfkbia, Ikbz, and Nfkbiz were also upregulated. NFκB target genes include many pro-inflammatory cytokines and chemokines, including IL-1β, IL-6, and CCL2 that are upregulated in this study.
G-protein signaling domain was the primary downregulated cluster related to signaling. Two groups of G-protein signaling modulators were represented: regulators of G-protein signaling (RGS) and G-protein coupled receptor kinases (GRK). Both of these groups act to dampen GPCR activity. RGS molecules accomplish this by enhancing the intrinsic GTPase activity of activated Ga subunits [23] while GRK proteins phosphorylate the active GPCR, making it a high-affinity target for arrestin binding which blocks G-protein binding and activation [24]. Thus the downregulation of these molecules may potentiate GPCR signaling at the tick bite site.
Immune Response
A number of clusters in the gene ontology analysis at 6 hpi were related to the immune response. These were cation homeostasis, anti-microbial response, negative regulation of myeloid cell differentiation, and B-cell, T-cell, and Toll-like receptor signaling. Within the cation cluster were transcripts for genes involved with iron, zinc, calcium, and proton transport or regulation. In particular, lactotransferritin, metallothionein 1, and metallothionein 2 have been shown to function in regulating reactive oxygen species production and scavenging [25], [26]. While some of the genes in this cluster are related calcium transport and may function in cell signaling, we suspect that regulating the oxidative status of the tissues near the bite site is the primary function of these genes.
Genes of interest in the anti-microbial cluster were beta-3 defensin (Def3b) and peptidoglycan recognition protein (Pglyrp1). Defensins are small positively charged cysteine-rich peptides with antimicrobial activity; interestingly, epithelial tissues but not neutrophils were the primary sources of mouse beta-defensins [27]. Def3b has wide spectrum anti-microbial activity against bacteria [28], fungi [29], and viruses [30]. Pglyrp1 has been shown to enhance intracellular killing of bacteria in neutrophils [31]. Thus early host responses to tick feeding include upregulation of potent anti-microbial proteins that could impact the transmission of tick-borne pathogens.
Genes within the negative regulation of myeloid cell differentiation and B-cell, T-cell and Toll-like receptor signaling clusters were transcription factors and signaling intermediates mentioned above (see Transcription factors and cell signaling pathways heading).
At 12 hpi, the only gene ontology cluster related to the immune response was inflammation and chemotaxis. The genes in this cluster were cytokines, chemokines, and related molecules. Chemokines are small peptides that are potent activators and chemoattractants for leukocytes, and play an important role at the sites of inflammation. Overall, nymphal tick feeding induced the expression of chemokines specific for neutrophil (Cxcl1 and 5) and monocyte (Ccl2, 6, 7, and 12) recruitment. In addition, Cxcl14 was upregulated, a chemokine specific for dendritic cell precursors [32] but without a defined function in the skin [33]. Increasing evidence suggests that small inflammatory mediators such as leukotrienes, prostaglandins, platelet activating factor, and complement initiate chemotaxis to sites of inflammation. This initial response is amplified by cytokine production that drives chemokine synthesis [34]. Our results support the upregulation of IL-1b, IL-6, and C1qb that may interact with the chemokine profile to maintain and amplify the chemotactic response. While the sequence of events could not be defined in this study, the gene expression profile strongly suggests the recruitment of neutrophils and monocytes to the bite site.
Ingenuity Pathways Analysis
All significantly modulated genes were submitted to Ingenuity Pathways Analysis. IPA has some advantages over DAVID in that it takes into account the fold change associated with each modulated transcript when determining biological significance. The most significant pathway in this dataset was related to acute inflammation and immune cell recruitment, supporting the DAVID analysis. The interactions between genes in this pathway were mapped over time (Figure 2), showing temporal increases in gene modulation in this pathway. In addition, all the genes modulated at any time point in the inflammatory response pathway were plotted to show temporal changes in gene expression (Figure 3). We believe this data suggests early tick feeding is characterized by an inflammatory response from the earliest time point that intensifies as feeding continues.
Gene expression data was entered into ingenuity pathway analysis software to identify potential canonical pathways modulated during tick feeding. Upregulated genes are orange/pink, downregulated genes are green, and unchanged or unsignificant genes are grey.
Gene expression data was entered into ingenuity pathway analysis software. All genes significantly modulated at any time point in the inflammatory response pathway from the IPA analysis were plotted showing temporal changes in gene expression. Genes significant at each time point are marked with an asterisk.
Validation of Microarray Data
Based on the significance of chemotaxis and inflammation in the gene ontology analyses, we chose a list of chemokines, cytokines, inflammatory molecules, and genes in the ontology term “viral reproduction” (GO:0016032) for validation by real-time PCR (Table 1). Specific validation targets were chosen by significant fold change in the array data and/or previously identified genes of interest in host anti-tick responses [13]. All significantly modulated genes at any time point in the validation experiment are shown in figure 4. While exact fold changes did not correlate, trends of up or down regulation were well preserved between the array and validation studies (Table S1). Ctse, Foxp1, IL-6, and Muc1 were modulated in the array but were not validated. Other genes such as Ifng, IL-3, Jak1, Stat6, and Vwf were not modulated in either study. Significant increases in chemokine gene expression begin at 1 hr post-infestation and intensify as time progresses, supporting the importance of chemotaxis seen in the array. Upregulation of IL-1β, IL-10, and IL-4ra along with the general inflammation group supports early inflammatory changes at the bite site. While genes in the viral reproduction group showed some significant modulation, changes between time points were not consistent, casting some doubt on the biological significance of these differences.
Validation targets were chosen based on significant fold change in the array study at 12 hpi and/or genes previously identified as important in host anti-tick responses [13]. Pre-optimized primers were purchased from Qiagen, and real-time PCR was performed as described in the methods section. All significantly modulated genes at any time point in the validation study are plotted. Significant changes in gene expression at individual time points are marked with an asterisk. Significance was assessed using the delta-delta Ct method for gene expression normalization and measurement, and the linear models in microarray analysis (LIMMA) package for statistical comparisons between infested and tick-free mice, as previously described [13].
Histopathology
Concurrently with gene expression measurements, we pursued histopathological analysis of the bite site to see if morphological changes could be correlated with array data. Resident cells such as keratinocytes, fibroblasts, dendritic cells, and mast cells initially sense cutaneous tissue injury or antigenic stimuli. These cells release rapidly synthesized or pre-formed pro-inflammatory and chemotactic molecules including histamine and eicosanoids. These molecules can stimulate cytokine production and endothelial cells to mobilize Weibel-Palade bodies containing pre-formed selectins to the cell surface. While these early events cannot be measured at the transcriptional level, the histopathological analysis clearly shows tissue damage from the insertion of the hypostome and degranulating mast cells (Figure S1) as early as 1 hr post attachment. Minor inflammatory changes consisting of a few inflammatory cells and a small amount of eosinophilic material near the tick hypostome were also observed. By 3 hrs post-infestation, inflammatory cells were readily evident, the eosinophilic material near the hypostome was much more pronounced, and the tissue architecture had the appearance of streaming toward the bite site, even in hypodermal muscle layers. This appearance suggests that ticks may initially insert the hypostome deeply and then retract it, pulling deeper tissues towards the epidermis. These changes intensify at 6 hrs post-infestation, leading to a very intense, neutrophil dominated inflammatory lesion by 12 hrs of tick feeding. Also visible at 12 hrs were potential areas of myositis, muscle necrosis, and increased congestion in blood vessels near the hypostome (Figure 5).
Skin biopsies were fixed in formaldehyde followed by decalcification prior to paraffin embedding. Five micron sections were stained with hematoxylin and eosin, as described in the methods section. The arrowhead marks a marginating neutrophil at 6 hrs PI 1000x, while the arrow marks areas of putative myositis/muscle necrosis at 12 hrs PI 100x.
Early Immunologic Response to Tick Bites
The early appearance of pro-inflammatory changes in transcription and histopathology was unexpected. Previous studies in our laboratory had suggested a minimal early host response [13], supporting many studies that have shown tick salivary components are capable of inhibiting nearly every aspect of cell recruitment. Ixodes ricinus saliva contains leukotriene B4 binding proteins that have been shown to reduce neutrophil migration [35], histamine binding proteins have been described from Rhipicephalus appendiculatus saliva [36], and numerous tick anti-complement molecules have been described [37], [38], [39]. The release of histamine, eicosanoids, and complement fragments are likely some of the earliest events in the chemotactic cascade. In addition, I. scapularis saliva has been shown to downregulate neutrophil beta-2 integrins, reduce phagocytic efficiency, and inhibit intracellular killing of Borrelia burgdorgeri [40]. The reduction in intracellular killing may be explained by salivary proteins that block super-oxide formation [41], and detoxify reactive oxygen species [42]. Tick salivary proteins have also been shown to bind human IL-8 [43] and chemokines such as Cxcl8 [44]. These studies show tick saliva can inhibit later events in the chemotactic cascade and also effector functions of neutrophils. Against this backdrop, the present study shows leukocytes such as neutrophils and pro-inflammatory gene transcription was initiated before 3 hours post-infestation. Thus despite the impressive arsenal of inhibitory tick salivary proteins, the host is able to mount a surprisingly timely immune response. Studies in mice with labeled neutrophils (enhanced GFP expression under the control of the lysozyme M promoter) have shown that neutrophils migrate into sites of sterile cutaneous injury as soon as 20 minutes post-injury. Neutrophil numbers then increased rapidly for 2 hrs when a plateau-phase was reached [45]. In a similar model using a larger wound (6 mm) and EGFP- labeled neutrophils, influx was measurable at 4 hrs and did not plateau until 2–3 days post wounding [46]. These studies suggest neutrophil chemotaxis into sites of cutaneous injury was initiated within 20 minutes, but the subsequent kinetics and final concentration of neutrophils may depend on other factors such as the size of the wound. In our study, very few neutrophils were visible at the bite site by 1 hpi. It should be noted that this is 1 hour after apparent tick attachment and hence represents the maximum length of attachment. Even so, it seems likely that the very early phase of neutrophil recruitment to nymphal tick bite sites is slower than that reported to sterile cutaneous wounds. At 3 hpi, appreciable neutrophils are present, and their numbers increase across our study, suggesting the plateau phase may not be reached during the time scale of the experiment.
Conclusions
Our gene-expression analysis reveals modulation of some pro-inflammatory genes at 1 and 3 hpi that intensifies to include genes related to antimicrobial responses, regulation of reactive oxygen species, wound healing, and signaling through AP-1, NFκB, MAPK, and G-protein pathways at 6 hpi. At 12 hpi, chemokines and cytokines consistent with a neutrophil-dominated immune response are prominent. Histopathology showed mast cell degranulation by 1 hpi and increasing neutrophil influx beginning at 3 hpi, supporting the array data. These results suggest early cutaneous host responses to tick feeding are more pro-inflammatory than expected, and highlight the importance of neutrophils and related pathways in tick engorgement and tick-borne viral transmission.
Supporting Information
Figure S1.
Histopathology of Ixodes scapularis nymphal bite sites at 1 hr PI. Skin biopsies were fixed in formaldehyde prior to decalcification and paraffin embedding. Sections were stained with Geimsa, as described in methods section. Mast cells appear as collections of dark purple granules.
https://doi.org/10.1371/journal.pone.0047301.s001
(TIF)
Table S1.
Correlation between the microarray and the real-time PCR analysis. FC, fold change; P val, p-value; hpi, hours post-infestation.
https://doi.org/10.1371/journal.pone.0047301.s002
(DOCX)
Author Contributions
Conceived and designed the experiments: DH ST. Performed the experiments: DH ST JRC. Analyzed the data: DH ST JRC JFA. Contributed reagents/materials/analysis tools: DH ST JRC JFA. Wrote the paper: DH ST.
References
- 1. Jongejan F, Uilenberg G (2004) The global importance of ticks. Parasitology 129 Suppl: S3–14
- 2. Wikel SK (1999) Tick modulation of host immunity: an important factor in pathogen transmission. International journal for parasitology 29: 851–859.
- 3. Schoeler GB, Manweiler SA, Wikel SK (1999) Ixodes scapularis: effects of repeated infestations with pathogen-free nymphs on macrophage and T lymphocyte cytokine responses of BALB/c and C3H/HeN mice. Experimental parasitology 92: 239–248.
- 4. Skallova A, Iezzi G, Ampenberger F, Kopf M, Kopecky J (2008) Tick saliva inhibits dendritic cell migration, maturation, and function while promoting development of Th2 responses. Journal of immunology 180: 6186–6192.
- 5. Binnington KC (1978) Sequential changes in salivary gland structure during attachment and feeding of the cattle tick, Boophilus microplus. International journal for parasitology 8: 97–115.
- 6. Kaufman WR (1989) Tick-host interaction: a synthesis of current concepts. Parasitology today 5: 47–56.
- 7. McSwain JL, Essenberg RC, Sauer JR (1982) Protein changes in the salivary glands of the female lone star tick, Amblyomma americanum, during feeding. The Journal of parasitology 68: 100–106.
- 8. Ribeiro JM, Alarcon-Chaidez F, Francischetti IM, Mans BJ, Mather TN, et al. (2006) An annotated catalog of salivary gland transcripts from Ixodes scapularis ticks. Insect biochemistry and molecular biology 36: 111–129.
- 9. Alekseev AN, Burenkova LA, Vasilieva IS, Dubinina HV, Chunikhin SP (1996) Preliminary studies on virus and spirochete accumulation in the cement plug of ixodid ticks. Experimental & applied acarology 20: 713–723.
- 10. des Vignes F, Piesman J, Heffernan R, Schulze TL, Stafford KC 3rd, et al (2001) Effect of tick removal on transmission of Borrelia burgdorferi and Ehrlichia phagocytophila by Ixodes scapularis nymphs. The Journal of infectious diseases 183: 773–778.
- 11. Ebel GD, Kramer LD (2004) Short report: duration of tick attachment required for transmission of powassan virus by deer ticks. The American journal of tropical medicine and hygiene 71: 268–271.
- 12.
Bouchard KR, Wikel SK (2005) Care, maintenance, and experimental infestation of ticks in the laboratory setting. In: Marquart WC, editor. Biology of Disease Vectors. San Diego: Elsevier Academic Press. 705–711.
- 13. Heinze DM, Wikel SK, Thangamani S, Alarcon-Chaidez FJ (2012) Transcriptional profiling of the murine cutaneous response during initial and subsequent infestations with Ixodes scapularis nymphs. Parasites & vectors 5: 26.
- 14. Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, et al. (2006) The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC molecular biology 7: 3.
- 15. Huang da W, Sherman BT, Lempicki RA (2009) Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic acids research 37: 1–13.
- 16. Huang da W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols 4: 44–57.
- 17. Dvinge H, Bertone P (2009) HTqPCR: high-throughput analysis and visualization of quantitative real-time PCR data in R. Bioinformatics. 25: 3325–3326.
- 18. Gu LH, Coulombe PA (2007) Keratin function in skin epithelia: a broadening palette with surprising shades. Current opinion in cell biology 19: 13–23.
- 19. Coulombe PA, Tong X, Mazzalupo S, Wang Z, Wong P (2004) Great promises yet to be fulfilled: defining keratin intermediate filament function in vivo. European journal of cell biology 83: 735–746.
- 20. Wong P, Coulombe PA (2003) Loss of keratin 6 (K6) proteins reveals a function for intermediate filaments during wound repair. The Journal of cell biology 163: 327–337.
- 21. Paladini RD, Takahashi K, Bravo NS, Coulombe PA (1996) Onset of re-epithelialization after skin injury correlates with a reorganization of keratin filaments in wound edge keratinocytes: defining a potential role for keratin 16. The Journal of cell biology 132: 381–397.
- 22. Schonthaler HB, Guinea-Viniegra J, Wagner EF (2011) Targeting inflammation by modulating the Jun/AP-1 pathway. Annals of the rheumatic diseases 70 Suppl 1i109–112.
- 23. Sierra DA, Gilbert DJ, Householder D, Grishin NV, Yu K, et al. (2002) Evolution of the regulators of G-protein signaling multigene family in mouse and human. Genomics 79: 177–185.
- 24. Gurevich EV, Tesmer JJG, Mushegian A, Gurevich VV (2012) G protein-coupled receptor kinases: More than just kinases and not only for GPCRs. Pharmacology & Therapeutics 133: 40–69.
- 25. Ward PP, Mendoza-Meneses M, Park PW, Conneely OM (2008) Stimulus-dependent impairment of the neutrophil oxidative burst response in lactoferrin-deficient mice. The American journal of pathology 172: 1019–1029.
- 26. Inoue K, Takano H, Shimada A, Satoh M (2009) Metallothionein as an anti-inflammatory mediator. Mediators of inflammation 2009: 101659.
- 27. Ganz T (2003) Defensins: antimicrobial peptides of innate immunity. Nature reviews Immunology 3: 710–720.
- 28. Bals R, Wang X, Meegalla RL, Wattler S, Weiner DJ, et al. (1999) Mouse beta-defensin 3 is an inducible antimicrobial peptide expressed in the epithelia of multiple organs. Infection and immunity 67: 3542–3547.
- 29. Jiang Y, Wang Y, Wang B, Yang D, Yu K, et al. (2010) Antifungal activity of recombinant mouse beta-defensin 3. Letters in applied microbiology 50: 468–473.
- 30. Jiang Y, Wang Y, Kuang Y, Wang B, Li W, et al. (2009) Expression of mouse beta-defensin-3 in MDCK cells and its anti-influenza-virus activity. Archives of virology 154: 639–647.
- 31. Dziarski R, Platt KA, Gelius E, Steiner H, Gupta D (2003) Defect in neutrophil killing and increased susceptibility to infection with nonpathogenic gram-positive bacteria in peptidoglycan recognition protein-S (PGRP-S)-deficient mice. Blood 102: 689–697.
- 32. Meuter S, Moser B (2008) Constitutive expression of CXCL14 in healthy human and murine epithelial tissues. Cytokine 44: 248–255.
- 33. Meuter S, Schaerli P, Roos RS, Brandau O, Bosl MR, et al. (2007) Murine CXCL14 is dispensable for dendritic cell function and localization within peripheral tissues. Molecular and cellular biology 27: 983–992.
- 34. Sadik CD, Luster AD (2012) Lipid-cytokine-chemokine cascades orchestrate leukocyte recruitment in inflammation. Journal of leukocyte biology 91: 207–215.
- 35. Beaufays J, Adam B, Menten-Dedoyart C, Fievez L, Grosjean A, et al. (2008) Ir-LBP, an ixodes ricinus tick salivary LTB4-binding lipocalin, interferes with host neutrophil function. PloS one 3: e3987.
- 36. Paesen GC, Adams PL, Harlos K, Nuttall PA, Stuart DI (1999) Tick histamine-binding proteins: isolation, cloning, and three-dimensional structure. Molecular cell 3: 661–671.
- 37. Schuijt TJ, Coumou J, Narasimhan S, Dai J, Deponte K, et al. (2011) A tick mannose-binding lectin inhibitor interferes with the vertebrate complement cascade to enhance transmission of the lyme disease agent. Cell host & microbe 10: 136–146.
- 38. Valenzuela JG, Charlab R, Mather TN, Ribeiro JM (2000) Purification, cloning, and expression of a novel salivary anticomplement protein from the tick, Ixodes scapularis. The Journal of biological chemistry 275: 18717–18723.
- 39. Tyson KR, Elkins C, de Silva AM (2008) A novel mechanism of complement inhibition unmasked by a tick salivary protein that binds to properdin. Journal of immunology 180: 3964–3968.
- 40. Montgomery RR, Lusitani D, De Boisfleury Chevance A, Malawista SE (2004) Tick saliva reduces adherence and area of human neutrophils. Infection and immunity 72: 2989–2994.
- 41. Guo X, Booth CJ, Paley MA, Wang X, DePonte K, et al. (2009) Inhibition of neutrophil function by two tick salivary proteins. Infection and immunity 77: 2320–2329.
- 42. Narasimhan S, Sukumaran B, Bozdogan U, Thomas V, Liang X, et al. (2007) A tick antioxidant facilitates the Lyme disease agent's successful migration from the mammalian host to the arthropod vector. Cell host & microbe 2: 7–18.
- 43. Hajnicka V, Kocakova P, Slavikova M, Slovak M, Gasperik J, et al. (2001) Anti-interleukin-8 activity of tick salivary gland extracts. Parasite immunology 23: 483–489.
- 44. Vancova I, Hajnicka V, Slovak M, Kocakova P, Paesen GC, et al. (2010) Evasin-3-like anti-chemokine activity in salivary gland extracts of ixodid ticks during blood-feeding: a new target for tick control. Parasite immunology 32: 460–463.
- 45. Ng LG, Qin JS, Roediger B, Wang Y, Jain R, et al. (2011) Visualizing the neutrophil response to sterile tissue injury in mouse dermis reveals a three-phase cascade of events. The Journal of investigative dermatology 131: 2058–2068.
- 46. Kim MH, Liu W, Borjesson DL, Curry FR, Miller LS, et al. (2008) Dynamics of neutrophil infiltration during cutaneous wound healing and infection using fluorescence imaging. The Journal of investigative dermatology 128: 1812–1820.