- Open Access
Exploring the abomasal lymph node transcriptome for genes associated with resistance to the sheep nematode Teladorsagia circumcincta
Veterinary Research volume 44, Article number: 68 (2013)
This study exploited Blackface lambs that varied in their resistance to the abomasal nematode parasite, Teladorsagia circumcincta. Infection of these lambs over 3 months identified susceptible (high adult worm count, high faecal egg count and low IgA antibody) and resistant animals that had excluded all parasites. Previous work had shown that susceptibility and resistance is dependent on the differential immune response to the parasite, which occurs within the abomasal (gastric) lymph node (ALN) that drains the site of infection. The Affymetrix ovine gene array was used to interrogate the transcriptome of the ALN to identify genes and physiological pathways associated with resistance. We used a bovine RT-qPCR array of 84 genes to validate the gene array, and also report digital gene expression analysis on the same tissues, reanalysed using the Oar v3.1 sheep genome assembly. These analyses identified Humoral Immune Response, Protein Synthesis, Inflammatory Response and Hematological System Development and Function as the two top-ranked networks associated with resistance. Central genes within these networks were IL4, IL5, IL13RA2 and in particular IL13, which confirmed that differential activation of Th2 polarized responses is critical to the resistance phenotype. Furthermore, in resistant sheep there was up-regulation of genes linked to control and suppression of inflammation. The identity of differentially-expressed chemokines and receptors in the resistant and susceptible sheep also begins to explain the cellular nature of the host response to infection. This work will greatly help in the identification of candidate genes as potential selectable markers of genetic resistance.
Gastrointestinal nematode parasites are the cause of major economic losses to the sheep agricultural industry  and the major species in cool temperate regions is the abomasal strongylid Teladorsagia circumcincta[1, 2]. The control of this parasite is largely by the use of broad-spectrum anthelmintics [3, 4] but the increasing incidence of drug-resistant parasites and concern of drug residues in meat  has led to the search for alternative methods of parasite management . The animals most susceptible to T. circumcincta are weaned lambs . Most lambs eventually suppress infection  through the development of IgE and IgA anti-parasite antibodies; but this takes more than 6 weeks of persistent infection with infectious larvae [9–11]. Mucosal mast cells have also been shown to play an important role in the limitation of larval colonization and expulsion of helminths [12, 13] and these also function largely in association with parasite-specific antibodies . However, some sheep in most flocks develop only low levels of helminth-specific antibodies and fail to control larval colonization and egg production. Indeed, IgA levels and faecal egg counts (FEC) have been used as selectable markers for resistance [8, 15, 16]. Furthermore, different sheep breeds show marked diversity in resistance to helminth infection [17–19]. Consequently, one strategy for the non-pharmacological control of parasites is the exploitation of genetic variation for resistance found within and between different sheep breeds [20, 21].
Selection for resistance can be based on quantitative measurements of one or more phenotypic traits such as FEC and IgA antibody levels [16, 22] but the identification of molecular markers is potentially a more reliable approach for high resolution selection . There are three approaches for the identification of such markers, quantitative trait locus (QTL) mapping, genome-wide association studies (GWAS) and candidate gene analysis . QTL mapping is of low power and requires extensive further work to identify candidate genes . GWAS is expensive, requiring very large numbers of samples; in addition, lack of current sheep genomic resources mitigates against high resolution analysis . The alternative candidate gene approach aims to evaluate the relationship between a phenotypic trait and a variation in a gene; this gene is selected by measuring differential expression in relation to a relevant phenotype. A number of studies have used sheep microarrays to identify genes and molecular pathways associated with host responses to abomasal nematode parasites in sheep. Most have analysed the transcriptome of the Haemonchus contortus or T. circumcincta infected abomasal mucosa [27–29] or afferent lymph cells draining that mucosa . In addition, a RT-qPCR assay has been developed to analyse a limited number of immune-inflammatory genes . However, the immune response to parasites in the abomasum takes place within the abomasal (gastric) lymph node (ALN) and the events within that node determine the quality and quantity of the immune response and consequently the clinical outcome of infection.
This current study exploited parasite-naïve Blackface lambs with diversity in their predicted genetic resistance to T. circumcincta, which were trickle-infected with L3 larvae to mimic natural infection . This regime resulted in lambs with a range of resistance as assessed by adult worm counts, FEC and IgA levels. Previous studies with these sheep used digital gene expression (DGE)  and RT-qPCR  to conclude that both resistance (no FEC/high IgA) and susceptibility (high FEC/low IgA) are active responses to infection; and that the inflammatory lesions of the susceptible sheep are associated with differential activation of Th17 T cells. Consequently the aim of this project was to investigate genes and physiological pathways associated with the differential activation of the immune response linked to the different disease outcomes. These pathways are likely to contain candidate genes as potential selectable markers for resistance to T. circumcincta infection. However, in this new study we use the novel Affymetrix Ovine Gene 1.1 ST whole-genome array, based on the homologous Oar v2.0 assembly, and focus on gene and pathway identity in relation to resistance and susceptibility.
Materials and methods
Animals and experimental design
Fifty-five female Blackface lambs (10–13 weeks old), from a flock previously used for quantitative genetic and QTL analyses , were housed in worm-free conditions. Ten lambs were sham infected controls; 45 lambs were infected experimentally with ~2300 infective L3 T. circumcincta larvae three times a week for 12 weeks and sacrificed two days after the last infection. The sham-infected controls (C) were twins of lambs in the infected group. At post mortem ten infected lambs had no detectable adult worms in the total abomasal contents, while the other infected lambs had a range of adult worm counts up to 11 300. The lambs selected for analysis were chosen to maximize the power of detecting differential expression. Consequently, animals were ranked according to their infection levels . The 7 most resistant lambs (R) had no detectable abomasal adult worms or faecal egg count (FEC), high IgA antibody levels and high body weight. The 7 most susceptible lambs (S) were those with the highest adult worm count (mean 6000, maximum 11 300), high FEC (mean 414, maximum 950), low IgA antibody levels and low body weight. Details of the animals, infection protocols, trait and population genetic analyses have been described previously [32, 33]. Animal experiments were approved by University of Edinburgh Ethical Review Committee and conducted under an Animals (Scientific Procedures) Act 1986 Project Licence.
Sample collection and total RNA isolation
Abomasal (gastric) lymph nodes (ALN) were removed immediately post mortem and stored at −80 °C in RNAlater (Ambion, Huntingdon, UK). Total RNA was isolated using the Ribopure Kit (Ambion) as described previously . RNA quality and integrity were assessed using a RNA 6000 Nano LabChip on the Agilent 2100 Bioanalyzer and quantified using a NanoDrop ND-1000 spectrophotometer; all had an RNA Integrity Number of > 7.5.
Whole-transcript expression analysis and profiling
Primary transcriptome analysis was by Affymetrix Ovine Gene 1.1 ST Array. Sense-strand cDNA was generated from total RNA (500 ng) subjected to two rounds of amplification (Ambion® WT Expression Kit). The obtained cDNA was used for biotin labelling and fragmentation by Affymetrix GeneChip® WT Terminal Labelling and Hybridization kit (Affymetrix). Biotin-labelled fragments of cDNA (5.5 μg) were hybridized to Affymetrix Ovine Gene 1.1 ST Array plates using the appropriate Hyb-Wash-Scan protocol for this plate and the Gene Titan Hyb Wash Stain kit for the reagents (Affymetrix). After hybridization every array plate was washed and stained before the array plates were scanned by the Imaging Station of GeneTitan System. Image generation and the resulting CEL files for analysis were produced in Affymetrix® GeneChip® Command Console® Software (AGCC) version 3.0.1. Initial QCs were performed in Expression Console. The obtained Affymetrix .CEL files were imported into the Genomics Suite software package version 6.13.0213 (Partek, St. Louis, MO, USA). The imported data were analysed at the gene-level, with exons summarized to genes, using the mean expression of all the exons of a gene. Background correction was carried out using the robust multiarray average (RMA) algorithm, with quantile normalization, median polish probe summarization, and log2 probe transformation. Differentially expressed genes were identified by ANOVA, genes with a fold change > 1.5 or < −1.5, and a false discovery rate (FDR) > 0.05 were kept, calculated using the Benjamini–Hochberg method to adjust P-values . Gene annotation was performed based on similarity scores in BLASTN comparisons against ovine or bovine sequences in GenBank.
Inflammatory cytokines & receptors RT2 profiler™ PCR array
Quantitative real-time RT-PCR (RT-qPCR) analysis was performed using the SABiosciences Cow Inflammatory Cytokines & Receptors RT2 Profiler™ PCR Array (Qiagen, Crawley, UK), which measures the expression of 84 genes that mediate the inflammatory response (Cat. no. 330231 PABT-011ZR). Total RNA was extracted as described above and treated with RNase-free DNase I (Qiagen) and RNeasy MinElute Cleanup Kit (Qiagen) according to the manufacturer’s protocol to eliminate DNA contamination. Each 0.8 μg sample of RNA was reverse transcribed using a RT2 First Strand Kit (Qiagen) before dilution with RNase-free water according to the manufacturer’s protocol. A real time PCR was performed on each cDNA sample using the RT2 Profiler™ PCR Array with the RT2 SYBR Green ROX FAST Mastermix on a Rotor-Gene Q cycler (Qiagen). The cycling profile was performed at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Melting curve analysis of PCR products confirmed the absence of secondary product. RT2 Profiler PCR Array Data Analysis v3.5 was used for data analysis. The data analysis was based on the ΔΔCt method with gene expression normalized to the reference gene YWHAZ.
Cloning of ovine gene fragments
Amplicons from the RT2 Profiler™ PCR Array were cloned and sequenced to confirm PCR primer specificity using the TOPO® TA Cloning® Kit for Sequencing (Life Technologies) and BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems) according to the manufacturer’s instructions. Cloned amplicon sequences were used in BLASTN comparisons against ovine or bovine sequences in GenBank to confirm identity.
Illumina digital gene expression
Full details of the RNAseq methods and analysis have been described previously . These sheep sequences were originally aligned against the Btau 4.0 bovine genome. Detailed protocols, metadata and all raw data are deposited at the ArrayExpress database  accession number E-MTAB-445. For the current study the raw data were reanalysed by alignment against the sheep genome assembly Oar v3.1  using Bowtie v0.12.8 . Only tags with phred mapping quality of at least 30, with a maximum of one base-pair mismatch and mapped to less than two genome locations were retained for further analysis. Furthermore, only genes that were mapped from five samples or more were included in the final analysis. Statistical analysis was performed using “R v2.15”, and Limma  within BioConductor 2.11  was used to calculate the differential gene expression between resistant (R), susceptible (S) and control (C) groups, including fold change and q value, which is analogous to an adjusted p value or false discovery rate (FDR). Only genes with a fold change greater than 1.5 and q value ≤ 0.05 (FDR 5%) were annotated.
Molecular network and pathway analysis
The Ingenuity Pathways Analysis (IPA) Spring Release (2013) Software (Qiagen) was used to identify networks of interacting genes and other functional groups from the datasets of differentially-expressed genes. DEG was analysed, by uploading the HUGO Gene Nomenclature Committee (HGNC) gene symbols for the sheep orthologues and fold change data to IPA.
Expression analysis by Affymetrix ovine gene array
All microarray data, metadata and protocols are available in the ArrayExpress database under accession number E-MTAB-1580. The Affymetrix Ovine Gene 1.1 ST Arrays identified 43 genes (Table 1) in the ALN that showed significant difference (fold change ≥ 1.5 and adjusted p value ≤ 0.05) in the R vs. C comparison, four genes in the S vs. C comparison and only three genes in the R vs. S comparison. All differentially-expressed genes in the S vs. C and R vs. S comparisons were also significantly differentially-expressed in the R vs. C comparison except VIRP2 that was 2.04 fold increased in susceptible sheep compared to controls.
Presentation of these data by heat map (Figure 1) illustrates that the overwhelming majority of the differentially-expressed genes are increased in the infected groups. Forty-one of the 43 genes in the R vs. C comparison are increased in the R group; all four of the genes in the S vs. C comparison were increased in the S group and all three in the R vs. S comparison were increased in the R group. IL13, COL6A5 and ACTG2 were significantly increased in the R animals when compared to both the S and C groups. CHI3L2, IGHE and FCER2 were significantly increased in both the R vs. C and S vs. C comparisons but not in the R vs. S comparison. VIPR2 was the only gene significantly increased in the S animals and not in the R group when compared to uninfected controls.
RT-qPCR array validation
Validation of the Affymetrix gene arrays was performed using the SABiosciences Cow Inflammatory Cytokines & Receptors RT2 Profiler™ PCR Array which consists of 84 key genes that mediate inflammation. This assay is optimized for bovine genes and data were obtained for 83 genes when used with sheep cDNA (Additional file 1). The mean sequence identity between these 83 bovine and ovine genes is 96% (minimum 91%, IL27; maximum 100%, TNFSF13B). Single distinct melt curves obtained for all genes, except for NAMPT (no data), in all samples confirmed the specificity of the assay. Furthermore, analysis of the sequences of ten sheep amplicons selected randomly (CCL11, CXCL13, IL2RB, IL2RG, IL4, IL6R, IL13, IL16, PF4, TNFSF4) and produced using the assay primers, showed them to be 100% identical to sequences of the respective sheep homologues.
Table 2 identifies those genes showing ≥ 1.5 fold differential expression in the three comparisons. Three genes were significantly differentially expressed (p ≤ 0.05) in the R vs. C comparison, the Th2 cytokines IL4 and IL13 were increased by 11.9 and 3.96 fold respectively and MIP was repressed by 2.16 fold, in resistant sheep. IL4 (3.83 fold) and IL13 (3.37 fold) were also significantly increased, and CCL5 (2.43 fold), CCL11 (2.11 fold), CXCR1 (2.49 fold), CXCR3 (2.25 fold) and TNFSF10 (2.03 fold) were significantly repressed in the susceptible animals when compared to uninfected controls. In the comparison of the two infected groups the Th2 cytokines IL5 (3.4 fold) and IL13 (3.54 fold) were significantly increased and the chemokine CCL5 (3.16 fold) was repressed in resistant sheep.
Spearman’s rank correlation analysis of the fold changes calculated for all 83 genes in the R vs. C and S vs. C comparisons (comparing 164 values) gave a low, but highly significant positive correlation (ρ = 0.37, P < 0.0001) between the data obtained by the Affymetrix and PCR arrays. A direct comparison of quantitative expression of eight selected genes (Figure 2) shows the close relationship between the relative levels of expression measured by the two independent methods.
Ingenuity pathway analysis
IPA was used to help characterize how individual differentially-expressed genes interact, and consequently influence biological processes that affect the development of resistance or the maintenance of susceptibility to chronic T. circumcincta infection. Analysis of the data obtained using the Affymetrix gene array identified seven networks (Table 3), five with the R vs. C comparison and one each with the S vs. C and R vs. S comparisons. The top two ranked networks for the R vs. C dataset were Humoral Immune Response, Protein Synthesis, Inflammatory Response (Figure 3) and Hematological System Development and Function, with P-scores of 48 and 28 respectively. The three top Bio Functions, beyond the p ≤ 10-10 threshold, within these networks (Table 4a) were Immunological Disease with 24 genes with the highest p value of 2.73 × 10-13, Inflammatory Disease with 41 genes p = 8.34 × 10-13, and Hypersensitivity Response with 19 genes p = 7.98 × 10-11 associated with Diseases and Disorders and Hematological System Development and Function with 37 genes p = 4.79 × 10-12, Tissue Morphology with 33 genes p = 4.79 × 10-12, and Humoral Immune with 13 genes p = 2.28 × 10-11, associated with Physiological System Development and Function.
The network Cell-To-Cell Signaling and Interaction, Hematological System Development and Function, Immune Cell Trafficking with a P score of 24 was the only network identified for the S vs. C dataset (Table 4b); and Antigen Presentation, Lipid Metabolism, Small Molecule Biochemistry with a P score of 18 was identified for the R vs. S comparison. Within these networks there were no Bio Functions beyond the 10-10 threshold (Table 4c).
The most significant network was Humoral Immune Response, Protein Synthesis, Inflammatory Response (R vs. C comparison), with IL13 as the central gene (Figure 2). IL13 was the top ranked gene in the R vs. C comparison, up-regulated 3.17 fold (p = 7.44 × 10-7) in resistant sheep. This was confirmed by the RT-qPCR analysis where it was 11.92 fold (p = 0.001) increased in the resistant group (Table 2). In this analysis it was also significantly up-regulated in susceptible animals in the S vs. C comparison (3.37 fold p = 0.03) and consequently 3.54 fold increased in resistant vs. susceptible sheep (p = 0.002).
Mapping of sequencing tags to the sheep Oar v3.1 genome assembly
A previous study  also performed transcriptome analysis on the abomasal lymph nodes of these resistant and susceptible sheep using Illumina digital gene expression analysis. These data were originally analysed in relation to the Bos taurus genome assembly (Btau4.0). Here, we reanalysed the same primary data (ArrayExpress E-MTAB-445) against the most recent Ovis aries genome assembly (Oar v3.1) to obtain more accurate gene mapping. The mean total number of reads of the samples from 15 sheep (5 resistant, 5 susceptible and 5 control) was 1 473 000, of which 288 753 mapped to Btau4.0, with a maximum of 1 mismatch from an average tag length of 17 bases; in contrast 865 844 tags mapped the Oar v3.1. A ~2-4 fold increase in the number of significantly differentially expressed genes (Additional file 2) was also noted when comparing Oar v3.1 to Btau4.0; 229 in Oar v3.1 and 131 in Btau4.0 in the R vs. C comparison, 150 and 37 in the S vs. C comparison, and 146 and 83 in the R vs. S comparison.
IPA analysis of the revised digital gene expression data identified the top two networks (Additional file 3) in the R vs. C comparison as Cellular Growth and Proliferation, Cell Morphology, Cell-mediated Immune Response with a P-score of 56 and Post-Translational Modification, Hematological Disease, Cell Cycle, P-score of 52. The top network in the S vs. C comparison was Post-Translational Modification, Cell Signaling, DNA Replication, Recombination, and Repair with a P-score of 54. The top Bio Functions (Additional file 4) within these networks include Cancer, Gastrointestinal Disease, Dermatological Disease and Conditions and Immunological Disease; however, none of these Bio Functions reached the p ≤ 10-10 threshold.
This project is the logical extension of our three previous studies [24, 32, 33] concerned with the genetics and immunology of resistant and susceptible Blackface sheep persistently infected with the common abomasal nematode parasite, T. circumcincta. The first  described host-parasite interactions in a single intensively-phenotyped cohort with variable susceptibility; which were exploited to help analyse the nature of the mature host response associated with differential resistance [24, 33]; although this misses early events associated with the development of immune responses The aim of the current study was to identify genes and physiological pathways associated with the differential activation of the immune response, linked to the maintenance of resistance and susceptibility. It is part of a larger project that eventually aims to identify candidate genes that could be used as selectable markers of resistance in these Blackface sheep as well as other commercial sheep breeds.
The study exploits a new genomics resource for the analysis of the sheep transcriptome; the homologous Affymetrix Gene 1.1 ST Array based on the Ovis aries Oar v2 genome assembly. This consists of 508 538, 25mer probes for 22 047 genes, which interrogate approximately 625 bases per gene covering all exons of each transcript. Also new to this study is the use and validation of the SABiosciences Cow Inflammatory Cytokines & Receptors RT2 Profiler™ PCR Array in sheep; enabling 83 RT-qPCR assays to be used to validate the arrays rather than the more usual 8 or 10 genes. In addition we reanalysed first generation Illumina digital gene expression data  using the latest sheep genome assembly (Oar v3.1) with significantly different results; more than 3 times the number of tags mapped to Oar v3.1 as originally mapped to Btau 4.0. The only other study  that has examined the sheep ALN transcriptome used a small array of 1480 annotated probes with 5373 unannotated expressed sequence tags. This identified only one differentially expressed gene (HSPA1A) in relation to H. contortus infection.
The sheep in this study that were predicted to show variation in resistance to T. circumcincta were trickle-infected regularly for 3 months to mimic natural infection and continual exposure. They were analysed when the mature immune response of the resistant animals had controlled and/or eliminated that infection. At the same time susceptible animals did not control infection, and retained adult nematodes that produced large numbers of eggs . A previous study showed that these susceptible sheep generated an active Th17 immune and inflammatory response that failed to control infection . In this current study the results of the array analysis focused on the nature of the mature host response in the resistant animals, which controls nematode colonization despite chronic exposure to infectious larvae.
Murine models of nematode parasite immunity, resistance and susceptibility, using the nematodes Heligmosomoides polygyrus, Nippostrongylus brasiliensis and Trichuris muris are linked to strong Th2 responses, indicated by high levels of the cytokines interleukin 4 (IL4) and IL13 with parasite immunity and resistance, and high levels of interferon γ (IFNγ) with susceptibility [41–43]. It is clear from the IPA in this study that the host response linked to parasite control is strongly associated with the network “Humoral Immune Response, Protein Synthesis, Inflammatory Response” and that many of the network genes are key regulators of a Th2 cell response and are strongly up-regulated in the resistant sheep. This is in agreement with the results from the array analysis of the mucosa from challenge-treated-reinfected “immune” sheep . IL4 and IL13 are the two cytokines that principally control Th2 differentiation [44, 45]; both were significantly increased in resistant sheep, IL4 (3.96 fold in R vs. C) in the RT-qPCR and IL13 in both array (4.76 fold in R vs. C and 3.17 fold in R vs. S) and RT-qPCR (11.92 fold in R vs. C and 3.54 in S vs. C). These cytokines have related receptors and IL13RA2 is also up-regulated in resistant sheep; binding to their cognate receptors stimulates the activation of STAT6 that controls the expression of GATA3, the prime transcription factor for Th2 differentiation . IL5 is the third archetypal Th2 cytokine, and is also increased in resistant sheep (3.4 fold in R vs. S by RT-qPCR). A major function of IL5 is in stimulating eosinophil maturation and localization; these are crucial effector cells in the exclusion of T. circumcincta through binding to parasite-specific IgG, IgE and IgA antibodies and subsequent degranulation . Eosinophils are most numerous in the mucosa of the resistant sheep  evidenced by the increased expression, by array, of both FCER2 (1.72 fold in R vs. S) the low affinity receptor of IgE, and the IL5 receptor (IL5RA 3.15 fold in R vs. C) both expressed by eosinophils (Table 1).
Sheep control larval colonization, worm development and egg production through the generation of parasite-specific IgA and IgE antibodies [7, 9, 11, 20]. Levels of these two antibody classes are highly negatively correlated with worm length and fecundity and FEC  and both IgA (IGHA 2.01 fold) and IgE (IGHE 17.21 fold) are significantly increased (both by array) in the resistant sheep. The generation of Th2 responses in the gastrointestinal tract seem to be induced by the events at the mucosal epithelial surface . The intestinal epithelium is activated by parasite antigens via innate receptors and secrete the cytokines IL25 and IL33 that act on the newly identified innate lymphoid cells (ILC2) , and TSLP that acts on dendritic cells (DC) . This leads to IL13 expression by ILC2 and Th2 activation by the DCs. Major parts of the receptors for these cytokines are up-regulated (by array) in resistant sheep; IL17RB (receptor for IL25) is increased 4.73 fold, IL1RL1 (IL33 receptor and the principal marker for ILC2) in increased 2.7 fold and CRLF2 (TSLP receptor) is increased 2.03 fold, in the R vs. C comparison. The endothelia and DCs in the abomasal lymph nodes of resistant sheep express significantly increased levels of HRH1 (1.57 in R vs. C). Histamine receptor positive DCs modulate Th1/Th2 balance by inhibiting IL12R1 signalling thus promoting Th2 responses .
However, resistance and susceptibility is not just a matter of Th1/Th2 discrimination as susceptible sheep also show increased expression of IL4 (3.83 fold in S vs. C) and IL13 (3.37 fold in S vs. C) by RT-qPCR and IGHE (13.25 fold in S vs. C) by array. IGHA is only marginally increased (but non-significant) in these susceptible animals. Nevertheless the mature tissue response of the resistant and susceptible sheep is distinct . The abomasal mucosa of resistant sheep, at least 6 weeks after the last positive FEC, has evidence of only minor pathological change with small numbers of infiltrating lymphocytes and eosinophils. In contrast the mucosa of susceptible sheep was grossly inflamed.
These pathological differences are reflected in the differential expression of the chemokines, molecules which control leukocyte movement. The high expression of CCL17 (3.33 fold in R vs. C in array) and CCL26 or eotaxin-3 (15.5 fold in R vs. C in array) begins to explain the nature of mucosal infiltrate of resistant sheep, which consists largely of lymphocytes and eosinophils [52, 53]. Indeed CCL17 is chemotactic for, and activates, CCR4+ Th2 cells; their expression of IL4 and IL13 stimulates the expression of CCL26, which is chemotactic for eosinophils and basophils. IL13 also stimulates the expression of CCL5 (CCL5 3.15 fold in R vs. C in RT-qPCR) by activated T and NK cells, which is chemotactic for memory T cells and promotes Th2 responses .
Changes to these and many of the other genes indicate that a principal component of the response of resistant animals also includes a repression of acute inflammation and tissue healing. The low levels of MIF (−2.16 in R vs. C in array) and SELE (−1.8 in R vs. C in array) begins to explain the paucity of neutrophils in resistant animals . A principal function of MIF is induction of pro-inflammatory cytokines, which stimulates SELE (E-selectin) expression and hence neutrophil localization ; however the level of SELE in resistant sheep is significantly repressed, possibly explaining the lack of neutrophils in the abomasal mucosa of these animals . Three other genes which play significant roles in the regulation of inflammation are increased in resistant sheep. ALOX15 (lipoxygenase 15) is induced by IL4 and IL13 and inhibits pro-inflammatory leukotrienes and suppresses neutrophil chemotaxis; NR1H4 plays an key role in cholesterol homeostasis and inhibits IL1β induced inflammation ; CD200R is the Ox-2 receptor expressed on myeloid cells and functions to down-regulate myeloid cell activation and therefore depress inflammation . Both resistant and susceptible sheep are equally affected for the first 6 weeks of trickle infection  but only resistant animals control infection. Consequently the abomasum of those animals begins to heal; and this is seen by the significant changes (R vs. S and R vs. C) to genes associated with healing, with an 8 – 10 fold increase in type VI collagen (COL6A5) and a 3 fold increase in enteric smooth muscle actin (ACTG2).
VIPR2 (2.04 in S vs. C) is the one gene in the array that is differentially-expressed in susceptible sheep but not in the resistant group. This is the vasoactive intestinal peptide receptor; ligand-receptor interactions leads to inhibition of IL2-driven T cell proliferation and chemoattraction , consistent with the immunopathology of the susceptible abomasum. The excretory/secretory antigens of T. circumcincta have been shown to induce the expression of FoxP3 by murine T cells in vitro ; however there is no evidence for any increase in Tregs in the parasite-free resistant or highly-infected susceptible sheep.
The ultimate aim of these studies was to identify candidate genes as potential selectable markers for resistance to T. circumcincta infection. Within this study we have highlighted a small number of physiological pathways associated with differential susceptibility in Blackface sheep; in particular identified genes linked to differential T cell polarization. The Affymetrix ovine gene array incorporates probes for all the annotated exons within the sheep genome and a major advantage of using this technology is that differential transcript usage can be identified directly. We are currently analysing the array data to identify transcript variants of genes within the top-networks that are differentially-expressed in relation to resistance and susceptibility.
Brunsdon RV: The economic impact of nematode infection in sheep: implications for future research and control. The economic importance of parasites of livestock in New Zealand. Edited by: Heath ACG. 1988, New Zealand: New Zealand Society for Parasitology, 4-16.
Urquhart GM, Armour J, Duncan JL, Dunn AM, Jennings FW: Veterinary Parasitology. 1987, Avon: Longman Scientific and Technical
Coles GC, Jackson F, Pomroy WE, Prichard RK, von-Samson-Himmelstjerna G, Silvestre A, Taylor MA, Vercruysse J: The detection of anthelmintic resistance in nematodes of veterinary importance. Vet Parasitol. 2006, 136: 167-185. 10.1016/j.vetpar.2005.11.019.
Sargison ND, Scott PR, Wilson DJ, Macrae AI, Penny CD: Teladorsagia circumcincta resistance to moxidectin and multiple anthelmintic groups in ewes following use of the persistent drug before lambing. Vet Rec. 2010, 167: 523-527. 10.1136/vr.c4478.
Cooper K, Whyte M, Danaher M, Kennedy D: Emergency slaughter of casualty cattle increases the prevalence of anthelmintic drug residues in muscle. Food Addit Contam Part A Chem Anal Control Expo Risk Assess. 2012, 29: 1263-1271. 10.1080/19440049.2012.685893.
Eady SJ, Woolaston RR, Barger IA: Comparison of genetic and nongenetic strategies for control of gastrointestinal nematodes of sheep. Livest Prod Sci. 2003, 81: 11-23. 10.1016/S0301-6226(02)00197-5.
Stear MJ, Strain S, Bishop SC: Mechanisms underlying resistance to nematode infection. Int J Parasitol. 1999, 29: 51-56. 10.1016/S0020-7519(98)00179-9.
Stear MJ, Strain S, Bishop SC: How lambs control infection with Ostertagia circumcincta. Vet Immunol Immunopathol. 1999, 72: 213-218. 10.1016/S0165-2427(99)00134-8.
Smith WD, Jackson F, Jackson E, Graham R, Williams J, Willadsen SM, Fehilly CB: Transfer of immunity to Ostertagia circumcincta and IgA memory between identical sheep by lymphocytes collected from gastric lymph. Res Vet Sci. 1986, 41: 300-306.
Halliday AM, McAllister HC, Smith WD: Kinetics of the local immune response in the gastric lymph of lambs after primary and challenge infection with Teladorsagia circumcincta. Parasite Immunol. 2010, 32: 81-90. 10.1111/j.1365-3024.2009.01173.x.
Murphy L, Eckersall PD, Bishop SC, Pettit JJ, Huntley JF, Burchmore R, Stear MJ: Genetic variation among lambs in peripheral IgE activity against the larval stages of Teladorsagia circumcincta. Parasitology. 2010, 137: 1249-1260. 10.1017/S0031182010000028.
Huntley JF: Mast cells and basophils: a review of their heterogeneity and function. J Comp Pathol. 1992, 107: 349-372. 10.1016/0021-9975(92)90010-R.
Henderson NG, Stear MJ: Eosinophil and IgA responses in sheep infected with Teladorsagia circumcincta. Vet Immunol Immunopathol. 2006, 112: 62-66. 10.1016/j.vetimm.2006.03.012.
Stear MJ, Bishop SC, Doligalska M, Duncan JL, Holmes PH, Irvine J, McCririe L, McKellar QA, Sinski E, Murray M: Regulation of egg production, worm burden, worm length and worm fecundity by host responses in sheep infected with Ostertagia circumcincta. Parasite Immunol. 1995, 17: 643-652. 10.1111/j.1365-3024.1995.tb01010.x.
Shaw R, Morris C, Wheeler M, Tate M, Sutherland I: Salivary IgA: a suitable measure of immunity to gastrointestinal nematodes in sheep. Vet Parasitol. 2012, 186: 109-117. 10.1016/j.vetpar.2011.11.051.
Sayers G, Sweeney T: Gastrointestinal nematode infection in sheep–a review of the alternatives to anthelmintics in parasite control. Anim Health Res Rev. 2005, 6: 159-171. 10.1079/AHR2005108.
Mugambi J, Bain R, Wanyangu S, Ihiga M, Duncan J, Murray M, Stear M: Resistance of four sheep breeds to natural and subsequent artificial Haemonchus contortus infection. Vet Parasitol. 1997, 69: 265-273. 10.1016/S0304-4017(96)01128-4.
Terefe G, Lacroux C, Andreoletti O, Grisez C, Prevot F, Bergeaud JP, Penicaud J, Rouillon V, Gruner L, Brunel JC, Francois D, Bouix J, Dorchies P, Jacquiet P: Immune response to Haemonchus contortus infection in susceptible (INRA 401) and resistant (Barbados Black Belly) breeds of lambs. Parasite Immunol. 2007, 29: 415-424. 10.1111/j.1365-3024.2007.00958.x.
Sayers G, Good B, Hanrahan JP, O’Donovan J, Mulcahy G, Sweeney T: Breed differences in mucosal and systemic antibody response to nematode infection in sheep: an important role for IgE?. Parasitology. 2008, 135: 71-80.
Gasbarre LC, Leighton EA, Sonstegard T: Role of the bovine immune system and genome in resistance to gastrointestinal nematodes. Vet Parasitol. 2001, 98: 51-64. 10.1016/S0304-4017(01)00423-X.
Dominik S: Quantitative trait loci for internal nematode resistance in sheep: a review. Genet Sel Evol. 2005, 37 (Suppl 1): S83-S96. 10.1186/1297-9686-37-S1-S83.
Stear MJ, Abuagob O, Benothman M, Bishop SC, Innocent G, Kerr A, Mitchell S: Variation among faecal egg counts following natural nematode infection in Scottish Blackface lambs. Parasitology. 2006, 132: 275-280.
Davies G, Stear MJ, Benothman M, Abuagob O, Kerr A, Mitchell S, Bishop SC: Quantitative trait loci associated with parasitic infection in Scottish blackface sheep. Heredity. 2006, 96: 252-258. 10.1038/sj.hdy.6800788.
Pemberton JM, Beraldi D, Craig BH, Hopkins J: Digital gene expression analysis of gastrointestinal helminth resistance in Scottish blackface lambs. Mol Ecol. 2011, 20: 910-919. 10.1111/j.1365-294X.2010.04992.x.
Slate J, Gratten J, Beraldi D, Stapley J, Hale M, Pemberton JM: Gene mapping in the wild with SNPs: guidelines and future directions. Genetica. 2009, 136: 97-107. 10.1007/s10709-008-9317-z.
Brown EA, Pilkington JG, Nussey DH, Watt KA, Hayward AD, Tucker R, Graham AL, Paterson S, Beraldi D, Pemberton JM, Slate J: Detecting genes for variation in parasite burden and immunological traits in a wild population: testing the candidate gene approach. Mol Ecol. 2013, 22: 757-773. 10.1111/j.1365-294X.2012.05757.x.
Diez-Tascon C, Keane OM, Wilson T, Zadissa A, Hyndman DL, Baird DB, McEwan JC, Crawford AM: Microarray analysis of selection lines from outbred populations to identify genes involved with nematode parasite resistance in sheep. Physiol Genomics. 2005, 21: 59-69. 10.1152/physiolgenomics.00257.2004.
Rowe A, Gondro C, Emery D, Sangster N: Genomic analyses of Haemonchus contortus infection in sheep: abomasal fistulation and two Haemonchus strains do not substantially confound host gene expression in microarrays. Vet Parasitol. 2008, 154: 71-81. 10.1016/j.vetpar.2008.02.014.
Knight PA, Griffith SE, Pemberton AD, Pate JM, Guarneri L, Anderson K, Talbot RT, Smith S, Waddington D, Fell M, Archibald AL, Burgess ST, Smith DW, Miller HR, Morrison IW: Novel gene expression responses in the ovine abomasal mucosa to infection with the gastric nematode Teladorsagia circumcincta. Vet Res. 2011, 42: 78-10.1186/1297-9716-42-78.
Knight JS, Baird DB, Hein WR, Pernthaner A: The gastrointestinal nematode Trichostrongylus colubriformis down-regulates immune gene expression in migratory cells in afferent lymph. BMC Immunol. 2010, 11: 51-10.1186/1471-2172-11-51.
Ingham A, Reverter A, Windon R, Hunt P, Menzies M: Gastrointestinal nematode challenge induces some conserved gene expression changes in the gut mucosa of genetically resistant sheep. Int J Parasitol. 2008, 38: 431-442. 10.1016/j.ijpara.2007.07.012.
Beraldi D, Craig BH, Bishop SC, Hopkins J, Pemberton JM: Phenotypic analysis of host-parasite interactions in lambs infected with Teladorsagia circumcincta. Int J Parasitol. 2008, 38: 1567-1577. 10.1016/j.ijpara.2008.04.011.
Gossner AG, Venturina VM, Shaw DJ, Pemberton JM, Hopkins J: Relationship between susceptibility of Blackface sheep to Teladorsagia circumcincta infection and an inflammatory mucosal T cell response. Vet Res. 2012, 43: 26-10.1186/1297-9716-43-26.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
ArrayExpress database. 2013, [http://www.ebi.ac.uk/arrayexpress/]
Sheep genome assembly Oar v3.1. 2013, [http://www.livestockgenomics.csiro.au/cgi-bin/gbrowse/oarv3.1/]
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Smyth GK: Limma: linear models for microarray data. Bioinformatics and computational biology solutions using R and bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Andronicos N, Hunt P, Windon R: Expression of genes in gastrointestinal and lymphatic tissues during parasite infection in sheep genetically resistant or susceptible to Trichostrongylus colubriformis and Haemonchus contortus. Int J Parasitol. 2010, 40: 417-429. 10.1016/j.ijpara.2009.09.007.
Maizels RM, Yazdanbakhsh M: Immune regulation by helminth parasites: cellular and molecular mechanisms. Nat Rev Immunol. 2003, 3: 733-744. 10.1038/nri1183.
Cliffe LJ, Grencis RK: The Trichuris muris system: a paradigm of resistance and susceptibility to intestinal nematode infection. Adv Parasitol. 2004, 57: 255-307.
Anthony RM, Rutitzky LI, Urban JF, Stadecker MJ, Gause WC: Protective immune mechanisms in helminth infection. Nat Rev Immunol. 2007, 7: 975-987. 10.1038/nri2199.
Finkelman FD, Shea-Donohue T, Morris SC, Gildea L, Strait R, Madden KB, Schopf L, Urban JF: Interleukin-4- and interleukin-13-mediated host protection against intestinal nematode parasites. Immunol Rev. 2004, 201: 139-155. 10.1111/j.0105-2896.2004.00192.x.
Wynn TA: IL-13 effector functions. Annu Rev Immunol. 2003, 21: 425-456. 10.1146/annurev.immunol.21.120601.141142.
Ho IC, Tai TS, Pai SY: GATA3 and the T-cell lineage: essential functions before and after T-helper-2-cell differentiation. Nat Rev Immunol. 2009, 9: 125-135. 10.1038/nri2476.
Beraldi D, McRae AF, Gratten J, Slate J, Visscher PM, Pemberton JM: Development of a linkage map and mapping of phenotypic polymorphisms in a free-living population of Soay sheep (Ovis aries). Genetics. 2006, 173: 1521-1537. 10.1534/genetics.106.057141.
Allen JE, Maizels RM: Diversity and dialogue in immunity to helminths. Nat Rev Immunol. 2011, 11: 375-388. 10.1038/nri2992.
Spits H, Artis D, Colonna M, Diefenbach A, Di-Santo JP, Eberl G, Koyasu S, Locksley RM, McKenzie AN, Mebius RE, Powrie F, Vivier E: Innate lymphoid cells - a proposal for uniform nomenclature. Nat Rev Immunol. 2013, 13: 145-149. 10.1038/nri3365.
Taylor BC, Zaph C, Troy AE, Du Y, Guild KJ, Comeau MR, Artis D: TSLP regulates intestinal immunity and inflammation in mouse models of helminth infection and colitis. J Exp Med. 2009, 206: 655-667. 10.1084/jem.20081499.
Gutzmer R, Langer K, Lisewski M, Mommert S, Rieckborn D, Kapp A, Werfel T: Expression and function of histamine receptors 1 and 2 on human monocyte-derived dendritic cells. J Allergy Clin Immunol. 2002, 109: 524-531. 10.1067/mai.2002.121944.
Yoneyama H, Harada A, Imai T, Baba M, Yoshie O, Zhang Y, Higashi H, Murai M, Asakura H, Matsushima K: Pivotal role of TARC, a CC chemokine, in bacteria-induced fulminant hepatic failure in mice. J Clin Invest. 1998, 102: 1933-1941. 10.1172/JCI4619.
Komiya A, Nagase H, Yamada H, Sekiya T, Yamaguchi M, Sano Y, Hanai N, Furuya A, Ohta K, Matsushima K, Yoshie O, Yamamoto K, Hirai K: Concerted expression of eotaxin-1, eotaxin-2, and eotaxin-3 in human bronchial epithelial cells. Cell Immunol. 2003, 225: 91-100. 10.1016/j.cellimm.2003.10.001.
Tekkanat KK, Maassab H, Miller A, Berlin AA, Kunkel SL, Lukacs NW: RANTES (CCL5) production during primary respiratory syncytial virus infection exacerbates airway disease. Eur J Immunol. 2002, 32: 3276-3284. 10.1002/1521-4141(200211)32:11<3276::AID-IMMU3276>3.0.CO;2-5.
Calandra T, Roger T: Macrophage migration inhibitory factor: a regulator of innate immunity. Nat Rev Immunol. 2003, 3: 791-800. 10.1038/nri1200.
Huang RB, Eniola-Adefeso O: Shear stress modulation of IL-1β-induced E-selectin expression in human endothelial cells. PLoS One. 2012, 7: e31874-10.1371/journal.pone.0031874.
Kim I, Morimura K, Shah Y, Yang Q, Ward JM, Gonzalez FJ: Spontaneous hepatocarcinogenesis in farnesoid X receptor-null mice. Carcinogenesis. 2007, 28: 940-946.
Rosenblum MD, Olasz E, Woodliff JE, Johnson BD, Konkol MC, Gerber KA, Orentas RJ, Sandford G, Truitt RL: CD200 is a novel p53-target gene involved in apoptosis-associated immune tolerance. Blood. 2004, 103: 2691-2698. 10.1182/blood-2003-09-3184.
Voice JK, Dorsam G, Chan RC, Grinninger C, Kong Y, Goetzl EJ: Immunoeffector and immunoregulatory activities of vasoactive intestinal peptide. Regul Pept. 2002, 109: 199-208. 10.1016/S0167-0115(02)00182-9.
Grainger JR, Smith KA, Hewitson JP, McSorley HJ, Harcus Y, Filbey KJ, Finney CA, Greenwood EJ, Knox DP, Wilson MS, Belkaid Y, Rudensky AY, Maizels RM: Helminth secretions induce de novo T cell Foxp3 expression and regulatory function through the TGF-beta pathway. J Exp Med. 2010, 207: 2331-2341. 10.1084/jem.20101074.
This project was funded initially by the Natural Environment Research Council (NERC) through its “Post-Genomics and Proteomics” programme (grant NE/D000645/1) and subsequently by the Biotechnology and Biological Sciences Research Council (BBSRC) Institute Strategic Programme Grant to The Roslin Institute. HW is a BBSRC Industrial CASE scholar, sponsored by Sheep Improved Genetics Ltd and Eblex (grant BB/J012653/1). AJ acknowledges funding from the 3SR project (Sustainable Solutions for Small Ruminants) – 7th Framework Programme, call identifier FP7-KBBE-2009-3, Project 245140. We thank Professor Josephine Pemberton (Institute of Evolutionary Biology, University of Edinburgh) and Dr Dario Beraldi (CRUK, Li Ka Shing Centre, Cambridge) for original collaboration on the larger project and critical appraisal of the manuscript; Joan Docherty (The Marshall Building) for animal husbandry and collection of weight data and blood samples; Barbara Craig for the original parasitological measurements; staff at Genepool and ARKGenomics for help with the RNAseq and Affymetrix respectively.
The authors declare that they have no competing interests.
JH conceived the study, collected the samples and is the principal investigator. AGG and JH designed the current experiments; AGG performed the Affymetrix and RT-qPCR experiments and was responsible for the statistical, bioinformatic and pathway analysis. DB performed the original digital gene expression study, including the bioinformatics. HW and AJ reanalyzed the digital gene expression in relation to Oar v3.1. JH drafted the manuscript. All the authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: T. circumcincta infected sheep. Differentially expressed genes in the R vs. C, S vs. C and R vs. S comparisons analysed using the Ovis aries Oar v3.1 genome assembly. Mean; ΔΔCt (mean CtGOI - meanCYWHAZ). (XLSX 19 KB)
Additional file 2: T. circumcincta infected sheep. Differentially expressed genes in the R vs. C and S vs. C comparisons. Mean R, S, C; mean tag numbers per 106 tags in Resistant (R), Susceptible (S), Control (C). FC R/C, S/C, R/S; fold change Resistant vs. Control (R/C), Susceptible vs. Control (S/C, Resistant vs. Susceptible (R/S). (XLSX 61 KB)
Additional file 3: Top networks identified by Ingenuity Pathway Analysis, from the Illumina digital gene expression data. Network P-scores [−log10 (P-value)] is the probability of a network being randomly generated. (DOCX 15 KB)
Top Bio Functions identified by Ingenuity Pathway Analysis, from Illumina digital gene expression data.
Additional file 4: P-values calculated by Fisher’s exact test. (DOCX 28 KB)
About this article
Cite this article
Gossner, A., Wilkie, H., Joshi, A. et al. Exploring the abomasal lymph node transcriptome for genes associated with resistance to the sheep nematode Teladorsagia circumcincta. Vet Res 44, 68 (2013). https://doi.org/10.1186/1297-9716-44-68
- Ingenuity Pathway Analysis
- Digital Gene Expression
- Digital Gene Expression Analysis
- Digital Gene Expression Data
- Abomasal Mucosa