- Open Access
Exploring the abomasal lymph node transcriptome for genes associated with resistance to the sheep nematode Teladorsagia circumcincta
© Gossner et al.; licensee BioMed Central Ltd. 2013
- Received: 21 May 2013
- Accepted: 2 August 2013
- Published: 8 August 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.
- Ingenuity Pathway Analysis
- Digital Gene Expression
- Digital Gene Expression Analysis
- Digital Gene Expression Data
- Abomasal Mucosa
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.
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
Differentially-expressed genes in ALN as assessed by Affymetrix ovine gene array.
R vs. C
S vs. C
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.
Differential expression of inflammatory cytokines and receptors in ALN as assessed by RT-qPCR.
Resistant vs. Control
Susceptible vs. Control
Resistant vs. Susceptible
Ingenuity pathway analysis
Top networks identified by ingenuity pathway analysis, from the Affymetrix ovine gene array.
Resistant vs. Control
Humoral Immune Response, Protein Synthesis, Inflammatory Response
Haematological System Development and Function, Haematopoiesis, Tissue Morphology
Cardiovascular System Development and Function, Cellular Movement, Gene Expression
Hereditary Disorder, Skeletal and Muscular Disorders, Tissue Morphology
Cell Signalling, Molecular Transport, Vitamin and Mineral Metabolism
Susceptible vs. Control
Cell-To-Cell Signalling and Interaction, Haematological System Development and Function, Immune Cell Trafficking
Resistant vs. Susceptible
Antigen Presentation, Lipid Metabolism, Small Molecule Biochemistry
Top bio functions identified by ingenuity pathway analysis, from the Affymetrix ovine gene array.
a. R vs. C comparison
DISEASES AND DISORDERS
CSF3R, SELE, PTGDR2, IL1RL1, IGHE, VIPR2, CCL23, NR1H4, IL13RA2, CCL17, CCL22, IL13, ADRA1D, HRH1, IL17RB, CCR4, CCL2, NLRP12, CD1A, CCL26, CFH,HRH4, CYSLTR2, FCER2.
2.73 × 10-13 - 6.71 × 10-03
GATA1, CRLF2, NR1H4, VIPR2, IL1RL1, IL13RA2, SOCS2, P2RX1, CCL17, SLC9A4, CCL22, HRH1, NFIL, TFF2, CCL2, CD1A, HPGD, CFH, CYSLTR2, FCER2, CSF3R, ALOX15, SELE, PTGDR2, IL5RA, CCL23, IGHE, CD200R1, CD1B, ULBP1, P2RX7, IL13, IL17RB, CCR4, NT5E, NLRP12, PILRA, CCL26, IGHA1, SH2D1B, HRH4.
8.34 × 10-13 - 6.71 × 10-03
ALOX15, PTGDR2, SELE, IL5RA, GATA1, VIPR2, IL1RL1, IGHE, CD200R1, CCL17, CCL22, P2RX7, IL13, IL17RB, CCR4, CCL2, CCL26, HRH4, FCER2
7.98 × 10-11 - 4.94 × 10-03
Dermatological Diseases and Conditions
ALOX15, PTGDR2, SELE, IGHE, VIPR2, NR1H4, CCL23, CCL17, CCL22, IL13, ADRA1D, HRH1, IL17RB, CCL2, CCR4, NT5E, CD1A, HPGD, HRH4, CYSLTR2, FCER2
1.78 × 10-09 - 6.71 × 10-03
CRLF2, NR1H4, IL1RL1, VIPR2, IL13RA2, CCL17, CCL22, HRH1, TH, COL6A1, CCL2, CD1A, HPGD, CFH, CYSLTR2, FCER2, CSF3R, ALOX15, SELE, PTGDR2, IL5RA, CCL23, IGHE, CD200R1, P2RX7, IL13, ADRA1D, GLIPR2, IL17RB, CCR4, NT5E, NLRP12, CCL26.
1.78 × 10-09 - 5.7 × 10-03
PHYSIOLOGICAL SYSTEM DEVELOPMENT AND FUNCTION
Haematological System Development and Function
GATA1, CRLF2, VIPR2, IL1RL1, SOCS2, IL13RA2, P2RX1, CCL17, CCL22, HRH1, NFIL3, TFF2, CCL2, CD1A, CFH, CYSLTR2, FCER2, CSF3R, ALOX15, PTGDR2, SELE, IL5RA, IGHE, CCL23, CD200R1, CD1B, ULBP1, P2RX7, IL13, ZBTB32, IL17RB, CCR4, NLRP12, NT5E, MYOCD, CCL26, CYP4A11.
4.79 × 10-12 - 6.71 × 10-03
GATA1, CRLF2, VIPR2, IL1RL1, NR1H4, IL13RA2, SLC9A4, CCL22, TH, HRH1, NFIL3,TFF2, CCL2, CFH, SMPD3, FCER2, CSF3R, ALOX15, SELE, PTGDR2, IL5RA, IGHE, CD200R1, P2RX, IL13, OXT, IL17RB, CCR4, NLRP12, NT5E, FLNC, CNN1, CYP4A11.
4.79 × 10-12 - 6.49 × 10-03
Humoral Immune Response
HRH1,PTGDR2, NFIL3, IL5RA, IL17RB, CRLF2, CCR4, IGHE, IL1RL1, IL13RA2, IL13, HRH4, FCER2.
2.28 × 10-11 - 6.53 × 10-03
Immune Cell Trafficking
IL1RL1, VIPR2, SOCS2, CCL17, CCL22, HRH1, NFIL3, TFF2, CCL2, CD1A, CFH, CYSLTR2, FCER2, CSF3R, ALOX15, PTGDR2, SELE, IGHE, CCL23, CD200R1, ULBP1, CD1B, P2RX7, IL13, IL17RB, CCR4, NT5E, NLRP12, CCL26.
5.04 × 10-10 - 6.71 × 10-03
Cell-mediated Immune Response
PTGDR2, SELE, GATA1, CCR4, CCL2, IL1RL1, CCL23, IL13RA2, CCL17, CCL22, IL13.
1.09 × 10-08 - 5.63 × 10-03
Lymphoid Tissue Structure and Development
CSF3R, ALOX15, SELE, GATA1, CRLF2, IL1RL1, IL13RA2, CCL17, CCL22, P2RX7, IL13, CCL2, CCR4, NT5E, CYSLTR2
3.82 × 10-07 - 5.7 × 10-03
Digestive System Development and Function
TFF2, NR1H4, IL13RA2, SLC9A4, IL13
2.41 × 10-05 - 4.94 × 10-03
b. S vs. C comparison
DISEASES AND DISORDERS
FCER2, IGHE, SELE, VIPR2.
1.76 × 10-06 - 3.31 × 10-02
FCER2, IGHE, IGKC, NT5E, SELE, VIPR2
4.64 × 10-05 - 4.10 × 10-02
PHYSIOLOGICAL SYSTEM DEVELOPMENT AND FUNCTION
Hematological System Development and Function
FCER2, IGHE, IGKC, SELE, NT5E, VIPR2.
1.76 × 10-06 - 4.57 × 10-02
Immune Cell Trafficking
FCER2, IGHE, NT5E, SELE, VIPR2.
1.76 × 10-06 - 4.57 × 10-02
c. R vs. S comparison
DISEASES AND DISORDERS
9.85 × 10-06 - 2.62 × 10-02
CCL26, CD1A, IL13, FCER2, IGHE, IGKC, NT5E, SELE, VIPR2.
9.85 × 10-06 - 2.62 × 10-02
PHYSIOLOGICAL SYSTEM DEVELOPMENT AND FUNCTION
Hematological System Development and Function
CCL26, CD1A, IL13.
2.35 × 10-05 - 4.76 × 10-02
Immune Cell Trafficking
CCL26, CD1A, IL13.
2.35 × 10-05 - 4.76 × 10-02
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.
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.
- 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.Google Scholar
- Urquhart GM, Armour J, Duncan JL, Dunn AM, Jennings FW: Veterinary Parasitology. 1987, Avon: Longman Scientific and TechnicalGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Maizels RM, Yazdanbakhsh M: Immune regulation by helminth parasites: cellular and molecular mechanisms. Nat Rev Immunol. 2003, 3: 733-744. 10.1038/nri1183.View ArticlePubMedGoogle Scholar
- Cliffe LJ, Grencis RK: The Trichuris muris system: a paradigm of resistance and susceptibility to intestinal nematode infection. Adv Parasitol. 2004, 57: 255-307.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Wynn TA: IL-13 effector functions. Annu Rev Immunol. 2003, 21: 425-456. 10.1146/annurev.immunol.21.120601.141142.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Allen JE, Maizels RM: Diversity and dialogue in immunity to helminths. Nat Rev Immunol. 2011, 11: 375-388. 10.1038/nri2992.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Calandra T, Roger T: Macrophage migration inhibitory factor: a regulator of innate immunity. Nat Rev Immunol. 2003, 3: 791-800. 10.1038/nri1200.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.