Transcriptomic analysis of the temporal host response to skin infestation with the ectoparasitic mite Psoroptes ovis

Genes up-regulated with infestation, gene expression clusters 1-4

Cluster 1 - One hour post infestation

Cluster 1 represents the immediate early (IE) group of genes with peak expression at 1 hpi. The top biological function represented was "cancer" with 34 of the 51 pathway-eligible genes in this cluster associated with this category. The enrichment for cancer related genes is related to the abundance of transcription factors represented in this cluster and to highlight this, the top 2 networks from the Ingenuity pathway analysis were combined, forming a single network highly enriched for transcriptional regulators (Figure 3). This combined network demonstrates key roles for tumour necrosis factor alpha (TNFα), NF-kB and JUN each forming distinct hubs in the network (see red circles, Figure 3). A number of the genes in this network are known to be transcriptionally regulated by NF-kB, playing an important role in the instigation of the pro-inflammatory response. Further classification of cluster 1 using the DAVID Pathway Analysis software [19, 20] identified enrichment for factors in the toll-like receptor (TLR) pathway. This was supported by IPA which identified that six genes from cluster 1 are members of the TLR signalling pathway (FOS, JUN, JUNB, JUND, ATF3 and MAP2K6). Indeed a number of genes from network 1 are classified as primary response genes for bacterial lipopolysaccharide (LPS), i.e. FOS, JUN, TNF, EGR1 and NR4A1 (Figure 3) [21, 22].

The NF-kB family members NFKB1, NFKB2 and RELB and the NF-kB inhibitor complex factors NFKBIA, NFKBIB, NFKBIE and NFKBIZ were up-regulated following infestation with sheep scab with expression levels peaking at 3 hpi before declining towards baseline levels by 24 hours. It is well recognised that NF-kB is activated following signalling through the TLRs [23] (e.g. TLR4), to investigate this we mapped the temporal gene expression changes within the TLR signalling pathway following infestation. At 3 hpi, six genes (MYD88, NFKB1, NFKB2, NFKBIA, TLR4 and TLR2) from within the TLR signalling pathway showed peak expression (cluster 2). At 6 hpi (cluster 3) a further three members of this pathway showed peak expression (CD14, EIF2AK2 and MYD88) and at 24 hpi (cluster 4) peak expression of a further two genes (IRAK4 and TLR6) was observed. This indicates that early interactions between mites and host skin may lead to activation of elements of the TLR pathway, in turn leading to NF-kB activation. These findings are supported by a recent study by Trompette et al., [24] which demonstrated that the HDM allergen Der p 2 was capable of functional mimicry of lymphocyte antigen 96 (LY96 or MD-2) an integral part of the TLR4 signalling complex involved in the cellular response to bacterial LPS [25]. In addition, the P. ovis antigen Pso o 2, which shows a high level of sequence similarity with Der p 2 [26], instigates an NF-kB mediated pro-inflammatory response in ovine keratinocytes in vitro (S.T.G Burgess, unpublished observations).

Due to the increase in IL12 production following stimulation of APCs with LPS, TLR4 signalling has been associated with a Th1-type response [27]. However, TLR4 has been shown to be important in the induction of Th2 mediated allergic responses and for an ability to induce DC2 maturation (dendritic cells that support a Th2-type response) following exposure to helminth antigen [28]. Eisenbarth et al., [29] showed that the level of LPS exposure determines the type of inflammatory response, with exposure to high levels of LPS triggering a Th1 type response and exposure to low levels of LPS leading to Th2 type responses. It is likely that in a sheep scab infection occurring on the surface of the skin, that low levels of LPS from commensal gram negative bacteria would be present on the host skin. In addition, Hogg et al., [30] using 16S rRNA analysis demonstrated the presence of a number of bacterial species associated with P. ovis, including the gram negative bacterium Serratia marcescens, as such it is possible that LPS (derived from either host skin or mites) is present at the site of infection, and in conjunction with mite-derived allergens could trigger TLR4 instigating a pro-inflammatory response. Finally, Hammad et al., [31], demonstrated that HDM extract induced asthma via a TLR4-mediated mechanism resulting in activation of the Th2 promoting cytokines thymic stromal lymphopoeitin (TSLP), CSF2 (granulocyte-macrophage colony stimulating factor (GM-CSF)), IL25 and IL33. IL25, IL33 and TSLP were not represented on the ovine microarray used in this study and therefore we were unable to determine their expression levels. However CSF2 was activated within 1 hpi (9-fold increase) and these cytokines (IL33, IL25, TSLP & CSF2) have previously been linked to allergic immunopathology [32–34].

TNFα is up-regulated in human epithelial cells following exposure to HDM allergens and has also been detected in early stage lesions in sheep scab [11, 35]. In support of these previous findings TNF was up-regulated (3.4-fold) in our dataset within 1 hpi, a response that declined by 3 hpi. TNFα, a pro-inflammatory cytokine the expression of which is controlled by NF-kB [36], is involved in cellular proliferation, differentiation and apoptosis and is also a potent pyrogen via stimulation of IL1 secretion [37]. Another factor of interest up-regulated at 1 hpi is tachykinin 1 (TAC1 or substance P, 2.3-fold). TAC1 is a neurotransmitter involved in induction of behavioural responses and as a vasodilator and secretagogue with well characterised roles in the itch response [38]. Up-regulation of TAC1 within the first hour of exposure could explain the rapid itching response observed with sheep scab [11].

In summary cluster 1 highlights the transcriptional response in the instigation of the host pro-inflammatory reaction to P. ovis infestation within the first hour. A number of key transcriptional regulators are involved in this reaction including NF-kB, AP-1, EGR1-3, TNFα and ATF3. The key signalling event, which is likely to be instigated by a mite-derived factor, closely mirrors that of the host response to LPS, potentially acting via a TLR to trigger NF-kB.

Cluster 2 - Three hours post infestation

Cluster 2 represents the early (E) group of genes showing a peak expression at 3 hpi. Pathway analysis identified the rapidly expanding host immune response with 62 of the 124 (50%) genes in this cluster associated with "inflammatory disease". The merging of the top 4 networks from the analysis of cluster 2 highlighted the up-regulation of NF-kB and pro-inflammatory genes (Figure 4). To enable the identification of the main hub factors this figure highlights only the direct interactions (Figure 4). Key roles for TLR2 (4.6-fold up), TLR4 (2.3-fold up), NF-kB and the pro-inflammatory cytokines IL1A (3-fold up), IL1B (15-fold up) and IL6 were identified (see red circles and oval, Figure 4), and are supported by enrichment for the "TLR" (p = 4.18E^-06), "NF-kB" (p = 6.69E^-09), "IL1" (p = 1.54E^-06) and "IL6" (p = 8.09E^-08) signalling pathways (Table 4). IL6 is a potent inducer of the acute phase response, with well characterised roles in a number of inflammatory diseases and is involved in a number of immune signalling pathways, i.e. TLR and IL10 signalling pathways [39]. IL6 demonstrated a 100-fold up-regulation at 1 hpi and a peak expression level of nearly 800-fold by 3 hpi before declining back to 15-fold up-regulated by 24 hpi. The implications for the sizeable and early up-regulation of IL6 following exposure to P. ovis are clear with the ability of this factor to induce acute phase protein synthesis, T-cell activation, stimulation of B cell Ig production and proliferation of keratinocytes [39]. This final point is of interest as sheep scab lesions are characterised by hyperkeratosis, evident by the presence of crusty scabs at the site of infection [7]. A similar up-regulation of IL6 has been observed in vitro in human monocyte and epithelial cells following exposure to the HDM allergen, Der p 1 [40, 41]. The expression pattern of IL6 was mirrored by that of suppressor of cytokine signalling 3 (SOCS3), a well characterised inhibitor of IL6 signalling [42]. Like IL6, SOCS3 expression peaked at 3 hpi and then declined towards baseline levels by 24 hpi. It is plausible that SOCS3 is co-activated with IL6 as part of a negative feedback mechanism, possibly acting to control the level of IL6 expression induced following infestation with P. ovis. Seki et al., [43] demonstrated that SOCS3 regulates the onset and maintenance of Th2-mediated allergic responses and that this expression correlated with the pathology of the Th2-mediated diseases atopic dermatitis and asthma. SOCS3 is a potent negative regulator of interferon gamma (IFN-γ) signalling in keratinocytes and may play a role in the polarisation of immune responses towards Th2 via suppression of IFN-γ production [44].

IL1-α and IL1-β are pro-inflammatory cytokines proteolytically processed to form potent mediators of inflammation [45]. IL1 is one of the main participants in the initiation of inflammation and is important in the pathogenesis of skin diseases such as psoriasis and atopic dermatitis [46–48]. As well as being up-regulated here, IL1α and IL1β protein levels are increased in human skin equivalents stimulated with the scabies mites S. scabiei[49]. IL1 is expressed and stored intracellularly in keratinocytes, being quickly released in response to infection or injury [45]. Once released, IL1 triggers expression of prostaglandins, i.e. cyclooxygenase-2 (COX-2 or PTGS2, also differentially regulated here, 15-fold up at 3 hpi) and a number of selectin molecules, guiding chemotaxis of immune cells to the site of infection [45]. IL1 also induces expression of the chemokine IL8 which was highly up-regulated here (temporal cluster 4 and Table 3) and is likely to result in activation of surrounding keratinocytes, release of growth factors, such as fibroblast growth factor 7 (FGF7, 2.8-fold up at 24 hpi) and CSF2 (10-fold up between 1-6 hpi), and stimulation of keratinocyte proliferation [50, 51]. Interestingly, extracts from S. scabiei have been shown to reduce IL8 production in keratinocytes, indicating the presence of differing responses to infestation with these two related mite species [52]. Suppression of IL8 by S. scabiei may be explained by the fact that unlike P. ovis, S. scabiei is a burrowing mite forming a more intimate interaction with the host than that observed with P. ovis which does not penetrate deeper than the stratum corneum [53]. Therefore S. scabiei may have greater need to evade the host immune response, unlike P. ovis which appears to rely on it as a food source. The IL1 surface receptor IL1R1 was also up-regulated here (2.2-fold, peaking at 6 hpi). This receptor binds both IL1α and IL1β leading to the activation of NF-kB, AP-1, MYC and EGR1 (all of which were up-regulated between 1 and 3 hpi) and the subsequent instigation of a pro-inflammatory response. The receptors IL1R2 (3.3-fold) and IL1RN (3.6-fold) were up-regulated at 6 hpi; these factors act as decoy receptors, dampening the response to IL1 through negative feedback mechanisms [54, 55]. IL1 activity is thought to be controlled by IL4 via its ability to induce the expression of IL1R2 [55]. These mechanisms may act to attenuate the effects of IL1 at the site of infestation and could be linked to the increasing Th2 type response through IL4 induction.

As well as peak expression for a number of genes involved in TLR and NF-kB signalling pathways, cluster 2 included genes involved in the instigation of the host inflammatory response. The expression levels of selectin-E (SELE, 7-fold up) and selectin-P (SELP, average 18-fold up) both peaked at 3 hpi whilst selectin-L (SELL, 12-fold up) expression peaked at 6 hpi. Selectins are cellular adhesion molecules with roles in initiating leukocyte rolling along the endothelial barrier prior to intercellular adhesion molecule-1 (ICAM1) mediated adhesion and extravasation. Expression of ICAM1 (9-fold) also peaked at 3 hpi, falling back towards baseline levels by 24 hpi. The neutrophil chemo-attractant molecule IL8 showed a biphasic pattern of expression with one peak of up-regulation at 3 hpi (average 300-fold up), followed by a brief hiatus at 6 hpi (average 230-fold up) before finally peaking at 24 hpi (average 450-fold up). Biphasic patterns of IL8 mRNA expression have been observed following LPS stimulation of bovine alveolar macrophages and human white blood cells with initial peaks in mRNA expression between 1-4 hours followed by second peaks between 6-24 hours [56, 57]. These studies demonstrated that the first wave of IL8 expression was due to the initial stimulation with LPS, whilst the secondary phase was due to local production of TNFα and IL-1 [56, 57]. Our results support these findings with IL1A, IL1B and TNF all upregulated between 1-3 hours post-infestation due to stimulation by mite derived factors and then a second wave of IL8 expression likely to be stimulated by this initial pro-inflammatory response. A biphasic pattern of IL8 mRNA and protein expression has also been demonstrated in afferent lymph following infection of sheep skin with the ectoparasite L. cuprinia, indicating that upregulation of IL8 mRNA may also be reflected at the protein level [58]. The rapid increase in expression of factors involved in immune cell extravasation, i.e. selectins and adhesion molecules, at the site of infestation with P. ovis has obvious implications for disease progression, especially when combined with the significant up-regulation of chemo-attractant molecules like IL8 (300-fold), IL6 (700-fold up), colony stimulating factor 3 (CSF3, 18-fold up), CSF2 (GM-CSF, 10-fold up), chemokine (C-C) ligand 2 (CCL2 or MCP-1, 36-fold up) and chemokine (C-X-C) ligand 2 (CXCL2 or GRO2, 23-fold up).

CCL2 is chemo-attractant for monocytes and basophils, but not for neutrophils and eosinophils, and has previously been shown to be involved in the pathogenesis of psoriasis and atopic dermatitis [59]. Previous studies have demonstrated similar responses to HDM antigen, including up-regulation of adhesion molecules and cytokines, such as CSF2, ICAM1, IL6 and CCL2[60]. CXCL2 was up-regulated in our dataset and has been shown to be chemotactic for neutrophils and to be involved in modulation of skin inflammatory responses [61, 62]. A recent study also demonstrated the up-regulation of the related chemokine CXCL1 in human keratinocytes following exposure to HDMs, indicating a conserved response [63]. Sheep scab is characterised by a pronounced dermal oedema dominated by infiltration of eosinophils, accompanied by macrophages, lymphocytes, plasma cells, neutrophils and mast cells [11, 12, 64]. Therefore these findings help us to further understand lesion development, elucidating the mechanisms by which a local inflammatory response in the skin results in the recruitment and activation of circulating and skin resident immune cells. In contrast to these findings S. scabiei has been shown to down regulate the expression of SELE and this combined with the down regulation of IL8 by S. scabiei has been hypothesised to contribute to inhibition of neutrophil extravasation in early scabies lesions and may be related to immune evasion [65].

Mite-derived proteases, including serine, cysteine and aspartic proteases trigger host inflammatory responses [10, 66]. A number of these act by triggering the protease activated receptors (PARs) leading to release of chemokines, such as IL8 and CCL2 and expression of the selectins SELE and SELP and ICAM1[67–69]. Probes for the PAR family members were absent on the ovine microarray in this study so we were unable to determine any changes in their expression. However evidence points to their involvement as other components of these signalling cascades were up-regulated during the time course of infestation with P. ovis and homologues of the HDM allergens involved in the triggering of PARs have also been identified in P. ovis, i.e. Der p 1, Der p 3 and Der p 9 [67–70].

Sheep scab is classically associated with a Th2 type response, supported by infiltration of CD4⁺ T-cells and γδ-T-cells, combined with a lack of CD8⁺ T-cell infiltration and an influx of eosinophils and mast cells [71]. As well as induction of a number of pro-inflammatory factors with peak expression at 3 hpi, we detected a number of pro-allergic factors, in particular those involved in the development of a Th2 type response: IL4 (9-fold up) and the IL4 receptor (IL4R, 6-fold up) showed peak expression at 3 hpi and combined with a lack of up-regulation of the Th1-inducing factors, [IFN-γ, IL12 and IL18] this supports the Th2-based nature of sheep scab. These findings are combined with the up-regulation of IL18 binding protein (IL18BP, 8-fold at 6-hpi) which sequesters IL18, further inhibiting IL18 and IFN-γ production and also with the up-regulation of the canonical Th2 marker, IL13 (20-fold at 6 hpi). IL4 is released from mast cells and basophils following exposure to Der p 1, pointing to a link between activation of the innate immune system and the subsequent activation of an adaptive immune response [72].

In summary cluster 2 highlights the gathering pace of the pro-inflammatory response to infestation with P. ovis, likely to be triggered by the pro-inflammatory transcription factors identified in cluster 1. This is supported by up-regulation of pro-inflammatory cytokines such as IL1A, IL1B, IL8 and IL6 and increased expression of selectins and adhesion molecules, all of which are proposed to be involved in migration of immune cells to the site of infestation in sheep scab. These events may further exacerbate the early pro-inflammatory response, leading to oedema and serous exudate production, providing a ready source of food for the mites. These early changes in gene expression are in accord with an early pro-inflammatory reaction, followed by a Th2 associated pro-allergic response, laying the foundations for the allergic nature of sheep scab pathogenesis.

Cluster 3 - Six hours post infestation

Cluster 3 represents the intermediate (IM) group of genes with peak expression at 6 hpi. As in cluster 2, a number of these genes (91 of the 207 pathway mapping eligible genes in this cluster) are involved in inflammatory disease, along with "inflammatory response" (80/207), "cell-cell signalling and interaction" (79/207), "immunological disease" (81/207) and "cellular growth and proliferation" (104/207). The main canonical pathway represented here was "T-helper cell differentiation" (p = 6.38E^-08), followed by the "IL10" (p = 3.72E^-06), "NF-kB" (p = 1.57E^-05), and "acute phase response" signalling pathways (p = 6.13E^-05) and "communication between innate and adaptive immune cells" (p = 2.6E^-04) (Table 4). Development of a pro-allergic Th2 type response appears to be a critical step in disease progression in sheep scab. Analysis of the genes enriched within the canonical pathway for "T-helper cell differentiation" supports this, with increased expression of CD86 (6-fold), IL2RG (6-fold), IL2RA (11-fold), IL10 (2.8-fold) and IL13 (20-fold).

Network 2 from cluster 3, showed enrichment (network score = 38) for genes involved in infectious and dermatological disease (Figure 5). This network highlighted the role of IL1 signalling at 6 hpi and the pivotal role played by NF-kB in this process is clear (see red circle, Figure 5). Increased expression of specific IFN-γ receptor genes, IFNGR1 and IFNGR2, can also be seen from this network (see red oval, Figure 5). However, in the absence of increased expression of IFN-γ, along with the skewing towards a Th2 type response it is probable that their expression is indicative of a general T-cell activation. This network also highlights the activation of pentraxin 3 (PTX3) which may be induced by an increase in local IL1β [73]. PTX3 is a soluble pattern recognition receptor (PRR) activated by inflammatory signals and also by TLR activation [74] and has roles in modulating the inflammatory response, interacting with complement component C1Q, enabling it to either activate or inhibit the complement cascade [75]. PTX3 is up-regulated in H292 cells exposed to HDM extract, indicating a conserved role in sheep scab pathogenesis [17]. PTX3 expression is down-regulated by the Th2 cytokines IL4 and IL13 and, in our dataset, PTX3 expression begins to decline by 24 hpi with an average 60-fold increase in expression compared to the average 270-fold peak at 6 hpi and coinciding with the peak up-regulation of IL13 at 6 hpi [76] and matrix metalloproteinase 1 (MMP1, average 270-fold up) which is itself regulated by IL1β [77]. The expression of MMP1 is induced by TNFα, and IL6 is released in response to cytokine signalling in skin subjected to UV-induced DNA damage [78, 79]. MMP1, also known as fibroblast collagenase, hydrolyses collagen fibrils in the extracellular matrix of the dermis, with important implications in sheep scab pathogenesis [80]. Other MMPs, notably MMP7 (3-fold) and MMP13 (2.3 fold), were up-regulated at 6 hpi. MMPs have been implicated in other allergic diseases, for example MMP9 is the predominant MMP in asthma and its expression is associated with enhanced airway inflammation [81]. Additional genes clustering with PTX3 in a sub-cluster within cluster 3 were superoxide dismutase 2 (SOD2, 30-fold up), metallothionein 1E (MT1E, 77-fold up) and metallothionein 2A (MT2A, 10-fold up). The acute phase protein SOD2 is activated in response to oxidative stress acting to destroy toxic radicals, protecting cells from toxic oxygen, its anti-oxidant functions have also been linked to its anti-inflammatory properties [82]. Metallothioneins have a cytoprotective role following the instigation of an inflammatory response, and can protect against allergic airway inflammation induced by ovalbumin (OVA) by suppression of IL1B expression [83]. This cluster of genes may be indicative of an orchestrated attempt to modulate the excessive inflammatory response induced by the presence of mite factors and may inform novel methods by which sheep scab could be controlled. These interactions can be observed in more detail in cluster 3, network 4 (see red ovals, Figure 6).

The TNF superfamily member, CD40 was up-regulated with peak expression (6.5-fold) at 6 hpi and a 5.4-fold induction at 3 hpi. The HDM allergen Der p 1 proteolytically cleaves CD40, rendering dendritic cells (DCs) less responsive to signalling through the CD40-CD40L pathway [84]. DCs matured in the presence of Der p 1 also induced production of less IFNG and IL12 and more IL4 from CD4+ T-cells, directing a Th2 type response to mite exposure [84]. CSF2 (GM-CSF), although not classified as significantly differentially expressed in our dataset did show a strong up-regulation following infestation, confirmed by qRT-PCR (9-fold up at 1 and 3 hpi, peaking at 10-fold by 6 hpi). CSF2 is over-expressed by keratinocytes in atopic dermatitis lesions and its release is up-regulated from keratinocytes stimulated with Der p 1 [85, 86]. Pastore et al., [87] showed that this up-regulation was due to the dysregulated activation of the transcription factor AP-1 (composed of JUN and FOS both up-regulated at 1 hpi). It is plausible that the up-regulation of FOS and JUN at 1 hpi with P. ovis could lead to increased expression of CSF2. CSF2 has been shown to be essential for DC development and is involved in the regulation of DC function; this could have important implications for sheep scab pathogenesis, by reducing DC function and interfering with antigen presentation by DCs [88]. Expression of IL10 peaked at 6hpi (2.8-fold). IL10 is an anti-inflammatory cytokine with roles in NF-kB inhibition, prevention of Th1 T-cell clonal expansion, activation of mast cells and Th2 T-cells [89–91]. The up-regulation of IL10 may be indicative of a concerted effort to dampen Th1 responses, further biasing the immune system towards a Th2 pro-allergic response and may influence the down-regulation of the NF-kB genes found in cluster 2 by 6 hpi.

IL2 receptor alpha (IL2RA) showed peak expression at 6 hpi (11-fold up), indicative of an influx of immune cells at the site of infestation with P. ovis. IL2RA (also known as CD25) has been shown to be enzymatically cleaved by Der p 1 and promotes development of allergic inflammatory responses [92]. Cleavage of IL2RA by Der p 1 resulted in diminished proliferation of T-cells and a reduction in IFN-γ secretion [92], presumably promoting a bias towards a Th2-type response due to the impaired proliferation of Th1-type cells. This is further supported by the observations that Der p 1 fails to impair T-cell proliferation in response to exogenous IL4 in mice, promoting a Th2-type response by modulating IL4 and IFN-γ [93, 94]. Lee et al., [95] previously characterised a P. ovis orthologue of Der p 1, termed Pso o 1 and it is possible that this may act in a similar manner by enzymatically cleaving IL2RA.

In summary, peak expression of the Th2 cytokines IL13 and IL4 in cluster 3 highlights the continuing bias towards a Th2 based immune response to infestation with P. ovis and shows the up-regulation of a number of factors involved in wound repair and dampening of the inflammatory response. Cluster 3 highlights methods by which mite-allergens may modulate the host immune response, for example through the cleavage of immune factors, i.e. CD25.

Cluster 4 - Twenty four hours post infestation

Cluster 4 represents the late (L) genes with peak up-regulation at 24 hpi. The top biological function represented was "inflammatory response" with 89 of the 235 genes associated with this category. A number of additional biological functions including, "cell-cell signalling and interaction" (86/235), "cellular growth and proliferation" (111/235) and "cellular movement" (65/235) were enriched. Enriched canonical pathways included "leukocyte extravasation signalling" (p = 2.96E^-05), "antigen presentation" (p = 3.9E^-05), "T-helper cell differentiation" (p = 1.74E^-04), "dendritic cell maturation" (p = 5.96E^-04) and "IL8 signalling" (p = 3.1E^-03) (Table 4). The largest fold change was observed for IL8 (600-fold induction) which represented its peak expression level. IL8 is a potent chemo-attractant and angiogenic factor involved in the chemotaxis of neutrophils, basophils and T-cells, and has also been shown to be involved in eosinophil chemotaxis [96]. Increased IL8 expression supports our previous findings from ovine keratinocyte cultures stimulated with P. ovis mite wash and whole mite extract [14]. However, the timing of this response differs between the in vitro and in vivo systems as the previous study had shown that IL8 expression levels began to decline by 24 hpi; whereas this time point represented the peak level in vivo[14]. However, it cannot be ruled out that this response continues to be activated beyond 24 hpi. This difference in timing may arise due to additional signals provided by immune cells present in vivo in the skin layer absent in the in vitro ovine keratinocyte cultures. Similar responses to those seen in this study with IL8 have also been observed in epithelial cells treated with HDM extract [17, 49].

Another cytokine with peak expression at 24 hpi was regakine (8-fold) which acts in synergy with either C5a or IL8 (both of which were up-regulated in our dataset) to recruit circulating leukocytes into inflamed tissue [97, 98]. This finding has implications for sheep scab which is characterised by high levels of the pro-inflammatory cytokine IL8 [14], and it is plausible that IL8 and regakine could work synergistically with other cytokines to enhance the inflammatory response to infestation with P. ovis. In addition, as described for cluster 6 (below), the HDM derived allergen Der p 3 cleaves complement components C3 and C5, producing the anaphylatoxins C3a and C5a, respectively [99]. Local production of IL8 could combine with the anaphylatoxin C5a, derived from allergen-mediated cleavage of C5, and with regakine acting synergistically to produce a more pronounced inflammatory response. Regakine has also been shown to be up-regulated in cattle skin infected with the tick Rhipicephalus (Boophilus) microplus[100]. Co-activation of C5a, regakine and IL8 by the recently identified Der p 3 homologue, Pso o 3 (S.T.G. Burgess, unpublished data) following infestation with P. ovis could potentially lead to increased levels of inflammation via enhanced migration of inflammatory cells to the site of infestation, further exacerbating disease pathogenesis [97, 98].

Figure 7 shows the canonical pathway for leukocyte extravasation overlaid with the up-regulated genes from cluster 4 and shows an increase in many of the factors involved, indicating the strong influence of this activity at 24 hpi. These findings are consistent with previous reports highlighting a rapid influx of immune cells at the site of infestation with P. ovis[7] and the expression of a number of adhesion molecules, such as selectins, E, L and P at 3-6 hpi as described here. Consistent with our findings at 3-6 hpi there was also a significant enrichment for factors involved in "T-helper cell differentiation" towards a Th2 type response (p = 1.74E^-04). An additional pathway highlighted at this stage was "dendritic cell maturation" (p = 5.96E^-04). Ten molecules (CD86, FCER1G, FCGR1A, FCGR3A, HLA-A, HLA-DMA, HLA-DMB, IRF8, PIK3CD and TRA@) from this pathway were up-regulated at 24 hpi and elements of this pathway overlap closely with network 3, cluster 4 highlighting the induction of an IgE based pro-allergic Th2-type response to infestation with P. ovis at this time point (see red ovals, Figure 8).

Cytotoxic T lymphocyte associated protein 4 (CTLA4) showed peak up-regulation (7.7-fold) at 24 hpi (Figure 8). CTLA4 is a cell surface marker, associated with T-regulatory (T-reg) cell function, which binds to CD80 and CD86 on antigen presenting cells (APCs) sequestering them and preventing their interaction with CD28 on T-cells [101–103]. T-reg cells may be involved in sheep scab pathogenesis with the infiltration of T-reg cells in lesional skin observed within 2 weeks of infestation with P. ovis[104]. Although probes for the T-reg cell marker FOXP3 were absent from the ovine microarray used here we observed a 1.8-fold increase in its expression at 3 hpi, rising to 2.4-fold by 6 hpi as determined by qRT-PCR, offering support for the potential involvement of T-reg cells in sheep scab pathogenesis. Figure 8 highlights the increased expression (3.7-fold) of the C-type lectin molecule CLEC7A (Dectin-1) at 24 hpi. Further C-type lectins were up-regulated at this time with increased expression of CLEC4E (Mincle) (26-fold), CLEC4A (DCIR) (6.4-fold), CLEC4G (6-fold), CLEC6A (Dectin-2) (5-fold), CLEC2D (4.4-fold) and the mannose receptor C type 1 (MRC1 or CD206) (3.3-fold). C-type lectins are non-toll PRRs involved in recognition of pathogen associated molecular patterns (PAMPs) [105]. They have roles in the migration of DCs and in determining the outcome of their interactions with lymphocytes [106].

A previous study demonstrated the involvement of MRC1 in the uptake of Der p 1 in human monocyte derived DCs [107]. This interaction was more efficient in HDM-allergic patients, suggesting that its involvement in pathogenesis of allergic disease [107]. Dectin-1 is a receptor for beta-glucan and this interaction is activated in respiratory epithelial cells following exposure to HDM extracts [108]. These findings are supported by Kobayashi et al., [61], who showed that activation of keratinocytes by fungal beta-glucan resulted in activation of Dectin-1 and also the non-toll PRR, NOD2 (up-regulated here, 2.2-fold at 3 hpi). Therefore activation of Dectin-1, along with additional C-type lectins and other non-toll PRRs, i.e. NOD2, by mite derived allergens, combined with activation of TLR4 by HDM Der p 2 [(cluster 2, above) and as described by Trompette et al., [24]] could trigger activation of DCs, providing a link between the innate and adaptive immune systems in sheep scab. The activation of innate immune factors, such as IL1A, IL1B, TNFα, CSF2 and IL8 (all up-regulated here within 24 hpi) activates cutaneous acquired immunity by enhancing antigen presentation on DCs [109]. In addition, Der p 1 has been shown to cleave the C-type lectins DC-SIGN and DC-SIGNR, resulting in loss of cell surface DC-SIGN and reduced ICAM3 binding [110]. ICAM3 is involved in the trafficking of DCs and in particular with polarisation towards a Th1 based immune response, thus reduced binding of ICAM3 following mite allergen-mediated cleavage of DC-SIGN and DC-SIGNR could further bias towards a Th2 based immune response to P. ovis[111, 112]. Genes encoding the major histocompatibility complex factors, HLA-A, HLA-DMA and HLA-DMB are highlighted in Figure 8, with peak expression at 24 hpi supporting the enrichment of antigen presentation. Figure 8 also shows the up-regulation of a number of genes involved in antibody dependent responses, namely the Fc fragment of IgE receptors, FCER1A (3-fold up), FCER1G (5.7-fold up) and the Fc fragment of IgG receptors, FCGR1A (CD64, 5.4-fold up) and FCGR3A (CD16, 4.7-fold up). The Fc gamma receptors, FCGR1A and FCGR3A are activated by immune complexes involving IgG and are involved in maturation of DCs via triggering of IL10 production [113, 114]. These molecules can trigger production of ICAM1, CD86 and CD40, involved in further propagation of the immune response [84, 92, 115], and were also up-regulated here. FCER1A and FCER1G are IgE receptors that play key roles in allergic diseases by coupling allergen with mast cells initiating inflammatory and immediate hypersensitivity responses [113, 116]. This results in release of inflammatory mediators and secretion of lymphokines, furthering the development of the local inflammatory reaction [113]. Although the expression of these Fc receptors peaks at 24 hpi they are up-regulated by greater than 2-fold by 3-6 hpi in our dataset. The implications of the increased expression of these receptors provides further indication of progression towards a Th2-type response at the site of infestation and highlights the speed at which these events may unfold. Therefore the up-regulation of these factors (with the exception of IL5) shown here within the first 24 hours is indicative of early activation of a Th2-mediated pro-allergic response in sheep scab.

Cluster 4 contains genes involved in terminal differentiation of keratinocytes, with the small proline rich proteins (SPRRs) SPRR2A (324-fold up) and SPRR2E (average 136-fold up) showing high levels of expression. In addition the S100 calcium binding proteins, S100A9 (90-fold) and S100A12 (66-fold) were up-regulated. These genes are components of the epidermal differentiation complex (EDC), a genomic locus encoding for proteins involved in terminal differentiation of keratinocytes and ensuring integrity of skin barrier function [117]. A similar up-regulation of these genes has been observed in psoriasis and atopic dermatitis, and as in this study these changes were accompanied by the down-regulation of other EDC genes such as loricrin and filaggrin (see cluster 8 for further details (below), Sugiura et al., [118]). SPRRs are predominantly expressed in squamous epithelium, where they contribute to the formation of the cornified envelope, providing structural integrity and reducing permeability [119]. They have previously been shown to be up-regulated in response to epidermal injury (i.e. UV radiation) and in response to pro-inflammatory molecules e.g. IL1 [120, 121]. In addition, Zimmermann et al., [122] showed that SPRR2A and SPRR2B are involved in regulation of allergic inflammation and that they are induced by IL13 in experimental allergic responses.

The up-regulation of SPRRs in scab infected skin is likely to be controlled by the inflammatory response, either through the action of IL1, IL13, or both. The up-regulation of EDC genes at this time may aid skin barrier remodelling in an attempt to reconstruct an effective barrier against exogenous factors following mite and/or mite allergen-mediated damage. In relation to this, previous studies have highlighted the ability of HDM allergens (Der p 1 and Der p 3), for which P. ovis homologues have been identified (Pso o 1 and Pso o 3, respectively) to cleave intercellular tight junctions, thus opening up the tight junction barrier and increasing epithelial permeability [123, 124]. The proteolytic activity of Der p 1 also reduces skin barrier function in nude mice [125], disrupting the critical barrier function of the skin mediated by the stratum corneum [126] increasing the risk of further allergen sensitisation [125].

The S100 factors S100A9 and S100A12 (90-fold and 66-fold up at 24 hpi) represent two of fourteen S100 factors encoded for within the EDC [127]. S100A2, an oxidative stress-regulated factor over-expressed in psoriasis [128], was also up-regulated here however its expression peaked at 6 hpi (2.6-fold). S100A9 is a stress induced factor which, together with S100A8, forms the heterodimeric complex calprotectin [129]. S100A9 is expressed at low levels in normal skin, however like S100A2 and S100A12; these factors are both highly expressed in psoriatic and atopic skin lesions [130–132]. Like S100A9, S100A8 is a marker of keratinocyte activation and has roles in wound healing and inflammatory cell chemoattraction [132], in addition, S100A9 displays anti-microbial activity and has a role in resistance to invasion by pathogenic bacteria [133]. It also acts as a pro-inflammatory mediator, up-regulating IL8 release and surface expression of ICAM1 [134]. As well as local release from keratinocytes, S100A8 and S100A9 are produced by circulating neutrophils further enhancing inflammatory responses by influencing local trafficking of leukocytes [135]. Importantly, S100A9 up-regulates transcription of genes under NF-kB control enhancing development of endotoxic shock in response to LPS [136]. Ehrchen et al., [137] demonstrated that calprotectin amplifies innate immunity by triggering TLR4, thus further promoting the inflammatory response and enhancing extravasation of leukocytes to the site of inflammation. The S100 proteins S100A8 and S100A9 play roles in the hyper-proliferation and abnormal differentiation of keratinocytes observed in psoriasis, an effect that is mirrored in the pathogenesis of sheep scab with the development of severe crusted lesions [7, 11, 138].

The increased expression of S100 genes in scab infected skin could lead to exacerbation of an already prolific inflammatory response.

In summary, cluster 4 highlights the activation and maturation of DCs, indicating that the initial pro-inflammatory response is beginning to influence an adaptive immune response to P. ovis infestation. This cluster demonstrates that this is biased towards a Th2 based pro-allergic response to infestation. Cluster 4 also shows the increased involvement of antigen presentation and leukocyte extravasation. It is clear that the early pro-inflammatory response observed in sheep scab develops into a Th2, pro-allergic response to the presence of mites which is more likely to lead to disease exacerbation than to effective parasite clearance.

Genes down-regulated with infestation, gene expression clusters 5-8

Cluster 5 - One hour post infestation

Cluster 5 represents the immediate early (IE) repressed genes with peak repressed expression at 1 hpi. There were insufficient genes in this cluster for effective pathway mapping. However, five genes were down-regulated at 1 hpi, namely apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 2 (APOBEC2), zinc activated ligand-gated ion channel (ZACN), heat shock 70kDa protein 8 (HSPA8), ADAM metallopeptidase with thrombospondin type 1 motif, 2 (ADAMTS2) and fibrinogen-like 2 (FGL2).

Cluster 6 - Three hours post infestation

Cluster 6 represents the early (E) repressed genes with peak repressed expression at 3 hpi. The biological functions enriched in cluster 6 were "infectious disease" (6 of the 19 pathway mapping eligible genes), "connective tissue disorders" (5/19), "dermatological disease and conditions" (4/19) and "antigen presentation" (2/19). This cluster showed enrichment for the complement system (p = 5.29^-12), with six complement genes (C3, C1QA, C1QB, C1QC, C1R and C1S) down-regulated at 3 hpi (Table 4). This is highlighted in Figure 9, which shows the components of the complement pathway regulated at this time.

The complement genes (C3, C1QA, C1QB, C1QC, C1R and C1S) clustered together, sharing a very similar expression profile. These genes have a repressed level of expression between 1-6 hpi. However, by 24 hpi their expression had returned to pre-infestation levels (between 1-2 fold up-regulated), indicating the potential presence of a transcriptional control mechanism to repress complement in the early stages of infestation with P. ovis. Interestingly, Vroling et al., [17] also observed a down-regulation of C3 in airway epithelial cells following exposure to HDM extract. The complement system is a well characterised target of the HDM cysteine protease Der p 1 which proteolytically cleaves complement components C3 and C5 resulting in production of the anaphylatoxins C3a and C5a and formation of the membrane attack complex [99]. C3a and C5a are biologically active peptides involved in recruitment and activation of DCs, eosinophils and mast cells [139]. C3a and C5a are pro-allergic and their generation could further influence the immune response in the direction of Th2. In addition to their pro-inflammatory properties, the production of C5a and C3a may, through exogenous proteases, result in down-regulation of complement components at the transcriptional level as part of a negative feedback mechanism, limiting the potentially damaging effects of high levels of anaphylatoxin. As discussed earlier, high levels of IL8 and regakine could act in synergy with C5a to further amplify the inflammatory response. The gene encoding the C5a receptor, C5AR, was also up-regulated 10-fold by 3 hpi, remaining high at over 7-fold up-regulated at 24 hpi. C5AR is involved in promotion of immune cell recruitment and inflammation and intriguingly cross-talk between this pathway and the TLR pathways, notably TLR4 and TLR2, has been reported [140, 141].

A family of inactivated serine protease paralogs, termed scabies mite-inactivated protease paralogs (SMIPPs), has been identified in the scabies mite S. scabiei[142]. These factors inhibit the human complement system by binding to molecules involved in separate arms of the complement cascade, C1q (classical pathway), mannose binding lectin (lectin pathway) and properdin (complement factor P (CFP), alternative pathway) [142, 143]. Inhibition of these factors is proposed to allow scabies mites to limit host complement production, thus limiting complement-mediated gut damage and creating a favourable micro-environment for mite survival [142, 143]. The genes encoding the three subunits that make up C1q, namely C1QA, C1QB and C1QC were down-regulated at 3 hpi here, indicating a potential inhibitory effect on the classical complement pathway in sheep scab infestation. However no change in expression was observed for mannose binding lectin and conversely CFP was up-regulated (6-fold) by 24 hpi, indicating a potentially positive effect on the alternative complement pathway in sheep scab. To date no P. ovis homologues of the S. scabiei SMIPPs have been identified, however EST data from the sheep scab mite remains limited.

In summary cluster 6 highlights the reduced expression of complement genes following infestation with P. ovis. The mechanism for this remains unclear however it is conceivable that mites specifically target the complement pathway to evade this potentially damaging host response or in order to further exacerbate the localised immune response.

Cluster 7 - Six hours post infestation

Cluster 7 represents the intermediate (IM) repressed genes with peak down-regulated expression at 6 hpi. The top biological function represented was "cancer", (78 of the 150 pathway mapping eligible genes) indicating a general down regulation in gene expression. Additional enriched biological functions included "cellular movement" (46/150), "gastrointestinal disease" (30/150) and "cellular growth and proliferation" (64/150). This cluster is dominated by genes encoding collagens, including COL4A1, COL5A2, COL1A2, COL1A1 and COL3A1. In a similar pattern to complement genes, these genes are repressed between 1-6 hpi with expression levels returning towards pre-infestation levels by 24 hpi. P. ovis and HDM extracts contain proteases capable of collagen degradation [144, 145]. In addition, interstitial collagenase (MMP1) which was 270-fold up-regulated at 6 hpi is involved in collagen breakdown in the extra-cellular matrix and is also accompanied by the increased expression of its inhibitory factor TIMP1 (6.5-fold up-regulated) at this time. Suppression of collagen related genes could be in response to elevated levels of inflammatory activity or as part of a general remodelling/wound healing event in response to infestation. Separation and disruption of collagen bundles has been observed in sheep skin following P. ovis infestation, indicating a role for collagen disruption in disease pathogenesis [11].

In summary cluster 7 highlights the reduced expression of a number of collagen genes at 6 hpi with P. ovis. The repressed expression of these genes may be part of a host strategy to affect localised tissue repair and wound healing, however it may also be due to the requirement for increased tissue access for immune cells extravasating to the skin.

Cluster 8 - Twenty four hours post infestation

The final cluster of genes contained the late (L) repressed genes with peak repressed expression at 24 hpi and as with cluster 4 (24 hpi up-regulated) represents the end stage of the time course. The top biological function represented was "cancer" (60 of the 109 pathway mapping eligible genes). Additional biological functions associated with this cluster included "cell death" (44/109), "dermatological diseases and conditions" (12/109) and "cellular growth and proliferation" (45/109). As observed with cluster 4, cluster 8 was enriched for genes involved in terminal differentiation of keratinocytes with a number of genes located within the EDC locus. The most repressed of these genes were loricrin (LOR, 26-fold down-regulated), filaggrin (FLG, 5.8-fold down), S100 calcium binding protein A3 (S100A3, 2-fold down) and late cornified envelope 1B (LCE1B, 2-fold down), however a further three EDC-associated genes; PDZ domain containing 1 (PDZK1); selenium binding protein 1 (SELENBP1) and endosulfine alpha (ENSA) were repressed 2-fold, at 24 hpi.

Down regulation of filaggrin expression is interesting from the perspective of mite proteases, in particular for the serine protease Der p 3 found in HDM extracts. A homologue of Der p 3 has been identified in S. scabiei (Sar s 3) and is capable of cleaving recombinant human filaggrin in vitro[146]. Using an antibody specific to human filaggrin the authors were able to demonstrate the presence of human filaggrin within the scabies mite gut, indicating active ingestion of this protein [146]. As stated earlier we recently identified a P. ovis homologue of Der p 3, termed Pso o 3, although the serine protease activity of this allergen has yet to be confirmed (S.T.G. Burgess, unpublished observations). Digestion of filaggrin by mite proteases could damage the epidermal barrier, reducing skin barrier function and rendering the epidermis susceptible to external factors, such as allergens and secondary infections, a scenario likely to be important in disease pathogenesis [7].

The mechanisms responsible for down-regulation of filaggrin and loricrin at this stage are unclear. However, expression of these factors, along with the involucrin (IVL), which was also down-regulated here (1.98-fold between control and 24 hpi) have been shown to be down-regulated through the suppressive action of Th2 cytokines, such as IL4 and IL13 [147, 148]. It is possible that this mechanism exists to prevent further cell differentiation in response to wound healing and the expanding inflammatory response. Previous studies on atopic dermatitis have shown that defective skin barrier function allows allergen penetration and initiates immunological reactions [149]. The altered expression of EDC genes and defective epidermal differentiation are implicated in the pathogenesis of other skin and allergic diseases, including atopic dermatitis, ichthyosis vulgaris, psoriasis and asthma and these changes in expression are linked to mutations and single nucleotide polymorphisms (SNPs) within EDC genes [150–154]. It is currently believed that loss-of-function mutations in EDC genes, such as filaggrin, enable increased allergen penetration and therefore heighten the risk of allergen sensitisation [155]. These findings could potentially lead to the identification of host genes responsible for the observed innate susceptibility or resistance of different breeds of sheep to P. ovis infestation. Loricrin, filaggrin and involucrin are crucial components of the epidermal barrier and their down-regulation in scab infected skin, potentially due to the increase in Th2 cytokine expression, may lead to further barrier dysfunction, increasing mite allergen exposure and further exacerbating pathogenesis. The reduced expression of filaggrin, loricrin and involucrin has also been observed in human atopic dermatitis skin lesions and, as found here, this was accompanied by up-regulation of the S100 calcium binding proteins S100A2 and S100A9 (see cluster 4 and also [118]).

Also down regulated (2-fold) at 24 hpi was claudin-3 (CLDN3); the claudins are involved in the formation of tight junctions in the epithelium and claudin-1 is enzymatically cleaved by Der p 1 resulting in tight junction barrier disruption [123]. In addition "tight junction signalling" (p = 2.35E^-03) was amongst the top canonical pathways associated with cluster 8 (Table 4), with the expression of six genes from this pathway repressed (ACTC1, CLDN3, CNKSR3, JUN, MYL1 and TIAM1). The disruption of tight junction signalling may be crucial in pathogenesis of sheep scab where allergen penetration beyond the stratum corneum would be dependent on their enzymatic breakdown.

A number of keratin encoding genes showed an approximate 2-fold down-regulation at 24 hpi. These included keratin 2 (KRT2), KRT34, KRT71, KRT83 and genes encoding keratin-associated proteins (KRTAP), namely KRTAP1-1, KRTAP12-2, KRTAP13-1, KRTAP16-1, KRTAP24-1, KRTAP6-1, KRTAP6-2, KRTAP7-1, KRTAP8-1 and KRTAP8-2. Keratin 1 (KRT1, 7-fold down) and the gene encoding keratin-associated protein (KRTAP16-3, 4.4-fold down) were more highly repressed and this repression began at an earlier stage (2-fold repressed at 1 hpi) peaking at 24 hpi. The down-regulation of these genes following infestation with P. ovis may result in further perturbation of the skin barrier and could be a response to the increasing inflammatory reaction and mechanical skin damage. Interestingly, KRT1 is a predicted target for digestion by the S. scabiei allergen, Sar s 3 (a homologue of Der p 3 (HDM) and of Pso o 3 from P. ovis), if this proves to be the case, the implications for skin barrier function following this digestion may be similar to those expected with the down-regulation of filaggrin [146].

In summary cluster 8 highlights the importance of EDC genes and combined with the up-regulation of other EDC genes within cluster 4 (also 24 hpi), indicates an important role for this gene complex in sheep scab pathogenesis. Other roles of note within this cluster are the reduced expression of genes involved in the formation of tight junctions and also for keratin and keratin-associated genes.