- Research article
- Open Access
Gene expression profiling of Naïve sheep genetically resistant and susceptible to gastrointestinal nematodes
BMC Genomicsvolume 7, Article number: 42 (2006)
Gastrointestinal nematodes constitute a major cause of morbidity and mortality in grazing ruminants. Individual animals or breeds, however, are known to differ in their resistance to infection. Gene expression profiling allows us to examine large numbers of transcripts simultaneously in order to identify those transcripts that contribute to an animal's susceptibility or resistance.
With the goal of identifying genes with a differential pattern of expression between sheep genetically resistant and susceptible to gastrointestinal nematodes, a 20,000 spot ovine cDNA microarray was constructed. This array was used to interrogate the expression of 9,238 known genes in duodenum tissue of four resistant and four susceptible female lambs. Naïve animals were used in order to look at genes that were differentially expressed in the absence of infection with gastrointestinal nematodes. Forty one unique known genes were identified that were differentially expressed between the resistant and susceptible animals. Northern blotting of a selection of the genes confirmed differential expression. The differentially expressed genes had a variety of functions, although many genes relating to the stress response and response to stimulus were more highly expressed in the susceptible animals.
We have constructed the first reported ovine microarray and used this array to examine gene expression in lambs genetically resistant and susceptible to gastrointestinal nematode infection. This study indicates that susceptible animals appear to be generating a hyper-sensitive immune response to non-nematode challenges. The gastrointestinal tract of susceptible animals is therefore under stress and compromised even in the absence of gastrointestinal nematodes. These factors may contribute to the genetic susceptibility of these animals.
Grazing ruminants are constantly exposed to natural challenge by gastrointestinal nematodes. Infection by such parasites leads to clinical disease and production losses and is a serious problem in animal husbandry. The economically important gastrointestinal parasites of sheep belong to the Order Strongylida and the family Trichostrongyloidea and include Teladorsagia (Ostertagia), Trichostrongylus, Nematodirus and Haemonchus spp. The host response to parasite challenge is complex and poorly understood due to its polygenic nature. The response varies from sheep breed to sheep breed and from animal to animal . The response can result in rapid or delayed expulsion of parasites and the host response is critical for determining subsequent parasite status . Resistance to nematodes is primarily due to acquired immunity and is largely a Th2 type response [3–8] although innate immunity is also thought to play a role . Sheep lines which have been selected to be resistant to one particular nematode species have been shown to have increased resistance to other nematode species [10–13] primarily due to a non-specific mechanism of parasite expulsion . Historically anthelmintic drenching has been used to control nematode infection, however, the efficacy of this treatment is decreasing as parasite resistance to anthelmintics increases and nematodes resistant to multiple classes of anthlemintics are now found worldwide . The use of chemicals in food production is also increasingly encountering public hostility .
An alternative method of helminth control is vaccination of host animals. This strategy has successfully combated many bacterial and viral diseases, however, despite extensive research there is currently no vaccine available against the major gastrointestinal parasitic nematode species [15–17]. This is most likely due to the number of nematode species infecting the host and the complex parasite life cycle which expresses different antigens at each stage.
A natural method of parasite control is breeding for host resistance. Resistance to internal parasites, as measured by faecal egg count (FEC), has moderate heritability (h2 = 0.23–0.41, ) and this method has been shown to make significant genetic gains in a variety of sheep breeds [18–24] due to the combination of reduced FEC and reduced pasture contamination. The current method of animal selection is cumbersome and inefficient, however, as animals are selected based on their phenotype. Genetic gain could be accelerated if animals were selected on genotype rather than phenotype.
Lines of Perendale sheep have been divergently selected for parasite resistance and susceptibility at AgResearch, New Zealand since 1986 and differ in FEC by 4.9 fold . The number of adult nematodes in the gastrointestinal tract of animals from the resistant line is also lower than that from the susceptible line . Resistant animals therefore have an increased genetic capability to respond to and subsequently reject parasites when challenged. The identification of genes involved in this process would allow the development of genetic markers, which could be used in marker assisted selection breeding programs. A previous study identified genes differentially expressed between the selection lines in response to parasite challenge, but the use of challenged animals confounded the resistance status of the animals with level of infection. This study also did not distinguish innate from acquired immunity . We therefore chose to look at those genes that were differentially expressed between the selection lines regardless of the level of infection by using naïve animals that had never been exposed to gastrointestinal nematodes. We did this using high-throughput DNA microarray technology. This allowed us to examine the relative expression of thousands of genes in a single experiment and to group genes into expression classes, providing insight into their biological function.
The microarray experiment examined gene expression in duodenum tissue from four resistant and four susceptible naïve Perendale lambs. Poly(A)+ RNA isolated from duodenum tissue was reverse transcribed into cDNA, fluorescently labelled and hybridised to ovine 20 k cDNA microarrays. The experiment design was a factorial dye swap design involving 16 slides, where every animal was compared to every animal in the opposite selection line, as described previously . One hundred ESTs showed differential expression (Figure 1) between the resistant and susceptible animals. Forty ESTs were more highly expressed in the resistant line while 60 were more highly expressed in the susceptible line. These ESTs were resequenced in order to verify the insert sequence. Confirmed ESTs were annotated using BLASTN against the human RefSeq database and a cut-off E value of 1 × 10-18. Thirty of the ESTs more highly expressed in susceptible animals could be verified and annotated in this manner giving 25 unique human RefSeqs. Sixteen of the ESTs more highly expressed in resistant animals could also be annotated. The lists of differentially expressed human RefSeqs is given in Tables 1 and 2.
Northern blotting was carried out to confirm differential expression of a number of genes (Figure 2). TFF3 encodes intestinal trefoil factor 3. This gene is a member of a family of trefoil factors that is involved in mucosal maintenance and repair and is known to be overexpressed during inflammatory processes . In accordance with the microarray results, this gene had a 20% increase in expression in the duodenum of the susceptible animals compared to the resistant (Figure 2A). While the degree of upregulation was not marked it was still significant (P = 0.05).
Differential expression of the pancreatic secretory trypsin inhibitor gene (SPINK1) was also confirmed. SPINK1 encodes a gene also known as TATI (tumour associated trypsin inhibitor) that is expressed mainly in the pancreas but is also expressed in the mucosa of the small intestine . One of the major roles of SPINK1 is thought to be the prevention of premature activation of pancreatic proteases. This has the effect of decreasing the rate of mucus digestion by luminal proteases within the stomach and colon. SPINK1 is also known to increase the proliferation of a variety of cell lines and to stimulate cell migration, implying that it may be involved in the healing response following injury . When measured by Northern blot, SPINK1 expression was increased two fold in the susceptible animals compared to the resistant animals (Figure 2B) and this increased expression was highly significant (P = 0.005). In particular two of the resistant animals had very low levels of SPINK1 expression although expression in the other two resistant animals was still lower than expression in any of the susceptible animals.
Differential expression of the immunity associated GTPase, GIMAP8, was also confirmed. This gene is a member of a novel family of GTPases conserved among higher plants and vertebrates . Human GIMAP proteins are known to be expressed most highly in the spleen and lymph nodes, but expression has also been detected in the digestive tract . GIMAP proteins are thought to be involved in the control of cell survival and response to infection and GIMAP8 has been shown to have anti-apoptotic functions. Expression of this gene is known to be decreased in the spleen of mice infected with the protozoan parasite Plasmodium chabaudi . GIMAP8 was expressed 2.6 fold more highly in the intestine of resistant animals compared to susceptible (Figure 2C) as measured by Northern blot. This differential expression is also highly significant (P = 0.0018). It is notable that the fold increase in expression of both GIMAP8 and SPINK1 is markedly higher when quantitated by Northern blot rather than when quantitated by microarray. Differential transcript expression of GIMAP8 was also observed with the four resistant animals predominantly expressing an approximately 3,300 bp transcript, although a larger transcript of approximately 4,100 bp could still be detected. In contrast three of the susceptible animals predominantly expressed the larger transcript although the smaller transcript was detectable, while one susceptible animal expressed the smaller transcript. The precise size of the GIMAP8 transcript is unknown in sheep. Transcript variants of this gene, with sizes of 4,700 bp and 4,200 bp, have been previously reported in human and transcript variants of size 4,200 bp and 2,900 bp have been reported in mouse [29, 30].
Gene Ontology (GO) terms significantly associated with the differentially expressed genes were found using EASEonline  and are listed in Table 3. A number of GO terms were significantly associated (P < 0.05) with the genes more highly expressed in the intestine of susceptible animals. These terms pertain to an organism's response to stress and stimulus and imply that despite the absence of nematodes from the environment of the susceptible animals, their gastrointestinal tract is still responding to insults or injuries, possibly by other pathogens. The only GO term significantly associated with genes more highly expressed in resistant animals was "cellular process" indicating that the intestine of these animals does not appear to be under stress.
Promoter regions of the 16 genes more highly expressed in resistant animals and the 25 genes more highly expressed in susceptible animals were analysed for common cis-regulatory motifs. The analysis was performed using the MEME motif prediction program . Only two motifs were significant and both were found in the promoter regions of the genes more highly expressed in susceptible animals. The sequence logo for these motifs is shown in Figure 3. MAST analysis indicated that despite individual motifs not being significant, the combination of motifs identified in promoter regions of both groups were unique to the group in which they were discovered as they were not significantly associated with the opposing group or with all the genes on the array by Fisher's exact test (Table 4). A separate MAST analysis indicated that the motifs identified in this study were also not significantly associated (P > 0.05) with genes previously reported to be differentially expressed in sheep exposed to gastrointestinal nematodes  (data not shown).
In order to identify which transcription factors could bind to the motifs, each motif was compared to the transcription factor binding sites in the TRANSFAC database . Each motif was shuffled 1,000 times and the comparison process repeated in order to establish its significance. Figure 3 lists the significant motifs discovered along with their best TRANSFAC hits. The full list of identified motifs and their TRANSFAC hits are available in Additional files 1 and 2. The two significant motifs found in the promoter regions of the genes more highly expressed in susceptible animals give B cell lineage specific activator protein (BSAP) and Peroxisome proliferative activated receptor gamma (PPARG) as their best hit. BSAP, also known as PAX5 is a transcription factor required for B cell development , however, it should be noted that the P value of this hit is not significant and so this motif may be bound by a transcription factor not represented in the TRANSFAC dataset. Interestingly, five of the six genes associated with the GO term "response to biotic stimulus" have this motif in their promoter regions, indicating this motif may be involved in the regulation of this response. Peroxisome proliferative activated receptor gamma (PPARG) is a ligand activated transcription factor which has an important role in adipocyte differentiation . PPARG is known to have anti-inflammatory effects and to play an important role in the maintenance of mucosal integrity in the intestine . This transcription factor may therefore play a role in co-ordinately regulating genes more highly expressed in the intestine of the susceptible animals.
This study identified a number of genes differentially expressed between lines of lambs differing genetically in their ability to become resistant to gastrointestinal nematodes. These genes were differentially expressed in the absence of nematode challenge and so represent basal expression differences between the selection lines. It is noteworthy that significant differences were observed in naïve animals as much previous work has shown that animals differing in host resistance primarily differ in the rapidity and strength of their acquired immune response rather than differing in their innate immunity [3, 5, 6]. To our knowledge this is the first report of constitutive differences in gene expression in naïve sheep differing in their genetic ability to respond to host infection.
Gene Ontology terms associated with the genes more highly expressed in susceptible animals included "response to stimulus", "response to stress", "defence response" and "response to pests, pathogens and parasites". This implies that the gastrointestinal tract of the susceptible animals is responding to stress even in the absence of nematode challenge. This is evidenced by the elevated expression of seven genes that are induced in response to stimulus, HLA-A, MSH6, GPX2, IFI35, UBD, SERPING1 and TFF3. HLA-A encodes an MHC class I heavy chain molecule that, in conjunction with B2M, presents endogenously derived peptides to CD8+ cytotoxic T cells. HLA-A has broad tissue expression and has been shown to be induced by a number of cytokines, in response to infection and in patients with the inflammatory bowel syndrome Crohn's disease [37–39]. MSH6 encodes a MutS homologue protein which, in complex with MSH2, forms the MutSα heterodimer, while MSH2 in complex with MSH3 forms the MutSβ heterodimer. Both these complexes are involved in mismatch repair and repair mutagen-induced lesions in DNA as well as errors in DNA replication . MutSα repairs both DNA mismatches and short 1–2 bp insertions or deletions (indels) while MutSβ repairs longer 2–6 bp indels . A change in expression of either MSH6 or MSH3 can subsequently lead to a change in the MutSα to MutSβ ratio. The MutSα complex is induced at the transcriptional level in response to radiation . Expression of MSH6 may be higher in the duodenum of susceptible animals in order to assist these animals to cope with toxic and mutagenic insults present in ingested food. GPX2 is another stress response gene more highly expressed in the intestine of the susceptible animals. This gene encodes the selenoprotein, gastrointestinal glutathione peroxidase 2, which catalyses the reduction of peroxides by reduced glutathione and protects the cell against oxidative damage. This gene is known to be induced in response to oxidative stress . The interferon inducible gene IFI35 is induced at the transcriptional level in response to interferons, and complexes with the N-Myc interacting protein, NMI, to form a high molecular weight cytosolic complex . The precise function of this complex, however, remains unknown. UBD encodes a small ubiquitin-like modifier protein also known as FAT10. UBD is induced by the cytokines IFNγ and TNFα  and can bind proteins and target them for degradation by the proteasome in a cytokine inducible, irreversible, ubiquitin-independent manner . SERPING1, also known as C1 inhibitor (C1-INH), encodes an IFNγ inducible  serine protease inhibitor of the complement and contact systems. This protein therefore exhibits an anti-inflammatory effect [48, 49]. SERPING1 is known to help prevent endotoxic shock  and is under investigation as a clinical treatment for a variety of diseases . TFF3 encodes intestinal trefoil factor 3. Trefoil factors are involved in mucosal protection and healing and are induced during inflammation and in response to gastrointestinal mucosa damage . The increased expression of so many stress response genes in the intestine of the susceptible animals appears to indicate that these animals are responding to gut insult and inflammation. Alternatively, these animals may constitutively express these stress response genes at inappropriate levels in the absence of challenge. Interestingly, stress response genes did not appear to be more highly expressed in the duodenum of susceptible animals after natural challenge with nematodes  indicating that these genes may be inappropriately regulated in the gastrointestinal tract of susceptible animals.
The source of stress to the intestinal tract of the susceptible animals remains unknown. Some of the stress response genes more highly expressed in susceptible animals are known to play a role in protecting the cell against oxidative damage. MSH6 is involved in repairing DNA lesions caused by oxidation  while GPX2 reduces DNA damaging agents . HMOX1 is also strongly induced by oxidative stress and is involved in heme degradation . Chronic gut inflammation is associated with enhanced production of leukocyte derived oxidants . Therefore the gastrointestinal tract of the susceptible animals may be suffering damage from reactive oxygen species. This could come from a number of sources such as aberrant cellular metabolism or phagocytic leukocytes responding to an infection. Indeed some of the stress proteins more highly expressed in the gastrointestinal tract of the susceptible animals are involved in the immune response and are inducible by pro-inflammatory cytokines. Therefore although free from gastrointestinal nematodes, the susceptible animals may be responding to infection by other viral or bacterial pathogens.
A number of pro-apoptotic genes are more highly expressed in the intestine of the susceptible animals. These genes promote apoptosis in response to stress. Apoptosis is a primary form of defence against infection, stress, damage or injury in the cell. MSH6, TRAF4, UBD and BMF can all induce a protective apoptotic response in cells [45, 57–59]. The soluable form of HLA-A has also been shown to trigger apoptosis of CD8+ T cells . Apoptosis can be induced directly by DNA damage or by cytokine signalling.
The genes more highly expressed in intestinal tissue of susceptible animals also included a number that encode proteins involved in protein degradation: these are PSMD12, which encodes a non-ATPase component of the 26S proteasome subunit, UBD which encodes a ubiquitin-like protein that targets proteins for degradation and BMSC-UbP which encodes a bone marrow stromal cell derived ubiquitin-like protein [46, 61, 62]. Another gene which can play a role in protein degradation is SEC61B, which encodes the beta subunit of the Sec61 protein translocator, which transports proteins across the endoplasmic reticulum. This complex is also known to play a role in retrograde transport of misfolded or degraded proteins back into the cytoplasm for proteasome degradation . Therefore there appears to be increased protein degradation in the intestinal tract of the susceptible animals.
In summary, a number of stress response genes appear to be induced in the intestinal tract of naïve susceptible animals. How the induction of these genes is regulated is unclear at present, however, the identification of two significant motifs in the promoter regions of these genes indicated that transcription may be co-ordinately regulated by the BSAP or PPARG transcription factors. The role of PPARG in regulating genes involved in nematode susceptibility could potentially be examined by treating animals with PPARG agonists in order to modulate PPARG activation. Motifs for BSAP or PPARG binding were not found in the promoter region of genes whose expression was not elevated in the susceptible animals, showing they are unique to these genes.
Only 16 genes were more highly expressed in the intestine of the resistant animals than the susceptible animals. These genes had a variety of biological functions and the only GO term significantly associated with these genes was "cellular process". However two of genes more highly expressed in the resistant animals are vital for maintaining a functioning and healthy immune system. RAC2 plays an important role in response to pathogens as it regulates neutrophil chemotaxis and superoxide production and deficiency of RAC2 leads to impaired host defences and neutrophilia [64, 65]. Deficiency of ITGB2, also known as CD18, can also lead to neutrophilia  and ITGB2 is also known to be important in innate immunity . Two apoptosis inducing genes are more highly expressed in the intestine of the resistant animals. These genes are DAP3 and TRADD. DAP3, a proposed nucleotide binding protein, is a major positive regulator of apoptosis and has been shown to be critical for anoikis . TRADD encodes a tumour necrosis factor receptor adaptor protein. This protein links the TNF receptor to the caspase pathway initiating apoptosis . Despite a few genes with similar functions being more highly expressed in the intestine of resistant animals, no significant motifs were identified in the promoter regions of these genes indicating that they may not be co-ordinately regulated at a transcriptional level. Alternatively, the small number of genes with significantly elevated expression provided little power to detect motifs responsible for co-ordinate regulation.
A previous study identified genes differentially expressed between the resistant and susceptible lines in response to nematode challenge . Interestingly, there were no genes that were consistently differentially expressed both pre and post-infection. However, the ubiquitin-like modifier gene, UBD, was more highly expressed in the susceptible animals prior to infection but more highly expressed in resistant animals after infection with gastrointestinal nematodes. UBD expression is induced by the cytokines IFNγ and TNFα. Expression of TNFα is known to be induced in intestinal lymph of genetically resistant Romney sheep during primary infection of naïve animals with T. colubriformis . Therefore, while UBD is expressed more highly in naïve susceptible animals, upon infection with nematodes the resistant animals may induce UBD expression to a higher level than the susceptible animals in a TNFα-dependent manner. The previous study also identified smooth muscle function and Major Histocompatibility Complex II expression as important mediators of parasite resistance . The current study did not identify smooth muscle or MHC II genes as differentially expressed between the lines in the absence of challenge and so shows that these genes are induced in response to challenge. Resistance to gastrointestinal nematodes may therefore be mediated, in part, by the ability to induce expression of these genes and this response is generated upon exposure to infection.
Despite divergent selection over many years relatively few genes were differentially expressed between the selection lines in the absence of nematode challenge. Many more genes were found to be differentially expressed between the lines in response to natural parasite challenge  however, the genes identified in the present study may contribute to an animals innate resistance or susceptibility. The genes more highly expressed in resistant animals had a variety of functions, but some were involved in maintaining a healthy immune system, while some were pro-apoptotic genes. A number of genes more highly expressed in the susceptible animals were related to cellular response to stress and infection indicating that the susceptible animals may have a compromised gastrointestinal tract, even in the absence of nematode infection, and this may contribute to their innate susceptibility. It is notable that the human orthologs of the differentially expressed genes are located on a number of different human chromosomes (Tables 1 and 2). It is not known whether the polymorphisms that give rise to the differences in expression reported here are due to cis or trans mutations. Although it is possible that a mutation in one or a small number of genes may give rise to the expression differences, the diversity of the actions and pathways of the genes differentially expressed is consistent with genetic resistance to nematodes being due to many genes with small effects rather than a mutation in a single locus. The recent release of the bovine genome , a closely related species, will greatly aid mapping and resequencing of sheep genes responsible for parasite resistance.
Selection lines and tissue collection
Resistant and susceptible lines of Perendale sheep have been selected based on faecal egg count (FEC) since 1986. These lines now differ in faecal egg count by 4.9 fold . Pregnant ewes were adjusted to concentrate feed, treated with anthelmintic, brought indoors and subsequently lambed indoors. The lambs were raised indoors to insure their immune system remained naïve with respect to gastrointestinal nematodes. During this period the ewes and lambs had access to commercially formulated sheep nuts and chopped Lucerne hay ad libitum. All animals were faecal sampled periodically to ensure no parasitic nematodes were present. Four female lambs per selection line were chosen and at an average age of 84 days (standard deviation, 6.8 days) these animals were sacrificed and duodenum tissue promptly collected. The tissue was frozen in liquid nitrogen and stored at -80°C. No adult nematodes could be detected in the abomasums or intestines of the animals confirming their naïve status. The average live weight of the lambs pre-slaughter was 23.4 kg (standard deviation, 3.8 kg). No significant differences were observed between the weights of the animals in the two lines at any stage. All procedures were approved by the AgResearch Invermay Animal Ethics Committee, formally constituted under the New Zealand Animal Welfare Act.
Total RNA was isolated from the duodenum of each animal using TRIzol (Invitrogen) and was cleaned using an RNeasy kit (Qiagen). RNA integrity was confirmed by denaturing agarose gel electrophoresis and RNA was quantitated using a NanoDrop® ND-1000 Spectrophotomoter (NanoDrop Technologies). First-strand cDNA was made from 25 μg of total RNA by anchored oligo(dT)20-primed reverse transcription incorporating amino-modified dNTPs, and was subsequently labelled indirectly by fluorescent coupling of Cy™3 and Cy™5 mono-reactive dyes (Amersham) using the SuperScript™ Indirect cDNA Labeling System (Invitrogen) according to the manufacturer's instructions.
Ovine cDNA libraries were prepared from 27 tissues by MWG Biotech (Germany). These libraries were single-pass sequenced from the 5' end generating expressed sequence tags (ESTs). The inserts, representing ovine expressed sequences, were then amplified in 96 well plates using the universal primers SP6 and T7 in 50 μl reactions (1.5 mM MgCl2, 2 mM each dNTP, 45 pmol primers and 2.5 units AB Red Hot Taq polymerase). The reactions were denatured at 94°C for 3 minutes then cycled 36 times at 94°C for 45 sec, 55°C for 45 sec and 72°C for 60 sec. Finally products were extended at 72°C for 5 min. The PCR products were verified by agarose gel electrophoresis, precipitated and resuspended in printing solution as described previously . At the University of Otago Genomics Facility an ovine microarray consisting of 19,968 spots was printed onto poly-L-lysine coated glass slides using an ESI array robot with 32 split pinheads depositing 0.6 nl with a 100 μm spot size. After printing the slides were UV irradiated to cross-link the DNA to the polylysine coating.
Slide hybridisation, scanning and normalisation
Slides were prehybridised by incubation for 45 min at 42°C in 50 ml of 5 X SSC, 0.1% SDS, 1% BSA. They were then rinsed twice in deionised water, once in isopropanol and dried. The labelled cDNA was denatured by heating at 95°C and then combined with 50 μl SlideHyb #1 (Ambion) and applied to slides. Slides were hybridised at 54°C for 16 hours in humidified chambers (Monterey). Post-hybridisation the slides were washed in the dark at 54°C for 10 min in 2 X SSC, 0.1% SDS, 5 min in 1 X SSC and 10 min in 0.1 X SSC. All buffers were filtered through 0.22 μm filters. The slides were dried and scanned in a ScanArray 5000 (Packard Biosciences) and the dual images collected in TIFF format. The combination and processing of the images was carried out using GenePix Pro (Axon Instruments). The data for each slide was normalised following the procedure of Baird et al., . For each EST on the array, the normalised data from all 16 slides was combined and a number of average statistics calculated. ESTs were excluded from further analysis if they had more than 6/16 bad spots. The remaining ESTs were sorted based on the modified T value  for the log ratio of the mean [Additional file 3]. ESTs where the tail of the probability plot was more extreme than the 95% confidence limit  were counted as differentially expressed (Figure 1). All the microarray information has been submitted into the National Centre for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) website . The accession number for the experiment series is GSE3738.
20 μg of total RNA from duodenum tissue was separated on a 1% formaldehyde agarose gel and transferred to Hybond nylon membrane (Amersham) by capillary transfer. Probes were generated by amplification from the original EST using the gene specific primers listed in Table 5. The PCR products were verified by sequencing, labelled with α32P-dCTP (Amersham) using the RadPrime DNA labelling system (Invitrogen) and unincorporated radioactivity removed using the High Pure PCR product purification kit (Roche). Membranes were prehybridised for 45 min at 42°C in 10 ml ULTRAhyb (Ambion). The probes were diluted 1:10 in 10 mM EDTA and denatured by boiling for 10 min followed by cooling on ice for 5 min. After addition of 0.5 ml of ULTRAhyb, the probe was added to the prehyb solution and incubated overnight at 42°C. The membranes were then washed twice for 5 min at room temperature in 2 X SSC, 0.1% SDS followed by 2 X 15 min in 0.1 X SSC, 0.1% SDS at 42°C. The membranes were then exposed to BioMax XAR film (Kodak). The images were scanned using an ImageScanner (Amersham) and the resulting images quantitated using ImageQuant TL (Amersham) and normalised to GAPDH mRNA levels. This gene was not differentially expressed in the microarray experiment making it a suitable choice of housekeeping gene for normalisation.
Microarray data interpretation
All ESTs in the ovine libraries, along with all ovine ESTs deposited in NCBI were assembled into contigs using CAP3  after an initial clustering step using BLAST. ESTs on the array were annotated by finding the human RefSeq (RefSeq release as at 11/4/2005) corresponding to the contig to which they belonged using BLASTN  and the following options: -e 0.01 -v 5 -b 5. Each EST was annotated with the top human RefSeq hit. In the case of the differentially expressed ESTs, a more stringent E value cutoff of 1 × 10-18 was applied. In cases where the EST matched more than one transcript variant of a gene then the top hit is listed. All analysis hereafter refers to the annotated RefSeq genes. Gene Ontology (GO) terms significantly associated with the differentially expressed genes were found using the Expression Analysis Systematic Explorer (EASE). The list of differentially expressed genes was submitted to EASEonline . The background list submitted included all human RefSeqs on the array. EASE calculates overrepresented functional gene categories compared to all the genes on the array .
A region of 1,500 bp upstream of the transcription start site of each of the RefSeq genes on the array was retrieved from the human genome browser at University of California-Santa Cruz  and sequences were masked for repeats through the retrieval process. The choice of promoter length was based on the report that approximately 75% of human core promoters lie within 1,500 bp of the transcription start site . Of the 9,238 unique RefSeqs on the array, promoter regions for 8,809 RefSeqs were retrieved. Motifs in the promoter regions of the differentially expressed RefSeq genes were detected using the MEME motif detection program . The motif lengths ranged between 6 and 12 bases, any number of repetitions of the motif was permitted, the reverse complement was allowed and the top 15 motifs were identified in each group The identified motifs were also compared to the transcription factor database, TRANSFAC release 8.1,  using the method described in Aerts et al., . Each motif was permuted 1,000 times and best matches and their level of significance are reported. The full list of motifs identified in the promoter regions of the genes more highly expressed in susceptible and resistant animals, and their top TRANSFAC hits are available in Additional files 1 and 2 respectively. In order to carry out a MAST analysis  motifs with a Pearson's correlation coefficient greater than 0.6 were identified and excluded from further analysis. Promoter regions of both groups of differentially expressed genes and all RefSeq genes on the array were subsequently screened for the presence of all unique motifs using MAST. Once again, reverse complement orientation of the motifs was allowed and the results obtained were analysed using Fisher's exact test based on the number of RefSeq promoters with a combined motif significance of P < 0.0001 in each of the groups.
Dominik S: Quantitative trait loci for internal nematode resistance in sheep: a review. Genet Sel Evol. 2005, 37 Suppl 1: S83-96. 10.1051/gse:2004027.
Balic A, Bowles VM, Meeusen EN: The immunobiology of gastrointestinal nematode infections in ruminants. Adv Parasitol. 2000, 45: 181-241.
Gill HS: Genetic control of acquired resistance to haemonchosis in Merino lambs. Parasite Immunol. 1991, 13 (6): 617-628.
Pernthaner A, Cole SA, Morrison L, Hein WR: Increased expression of interleukin-5 (IL-5), IL-13, and tumor necrosis factor alpha genes in intestinal lymph cells of sheep selected for enhanced resistance to nematodes during infection with Trichostrongylus colubriformis. Infect Immun. 2005, 73 (4): 2175-2183. 10.1128/IAI.73.4.2175-2183.2005.
Courtney CH, Parker CF, McClure KE, Herd RP: Resistance of exotic and domestic lambs to experimental infection with Haemonchus contortus. Int J Parasitol. 1985, 15 (1): 101-109. 10.1016/0020-7519(85)90107-9.
Stear MJ, Strain S, Bishop SC: Mechanisms underlying resistance to nematode infection. Int J Parasitol. 1999, 29 (1): 51-6; discussion 73-5. 10.1016/S0020-7519(98)00179-9.
McClure SJ, Emery DL, Wagland BM, Jones WO: A serial study of rejection of Trichostrongylus colubriformis by immune sheep. Int J Parasitol. 1992, 22 (2): 227-234. 10.1016/0020-7519(92)90106-U.
Bisset SA, Vlassoff A, Douch PG, Jonas WE, West CJ, Green RS: Nematode burdens and immunological responses following natural challenge in Romney lambs selectively bred for low or high faecal worm egg count. Vet Parasitol. 1996, 61 (3-4): 249-263. 10.1016/0304-4017(95)00836-5.
Aumont G, Gruner L, Hostache G: Comparison of the resistance to sympatric and allopatric isolates of Haemonchus contortus of Black Belly sheep in Guadeloupe (FWI) and of INRA 401 sheep in France. Vet Parasitol. 2003, 116 (2): 139-150. 10.1016/S0304-4017(03)00259-0.
Woolaston RR, Barger IA, Piper LR: Response to helminth infection of sheep selected for resistance to Haemonchus contortus. Int J Parasitol. 1990, 20 (8): 1015-1018. 10.1016/0020-7519(90)90043-M.
Gray GD, Barger IA, Le Jambre LF, Douch PG: Parasitological and immunological responses of genetically resistant Merino sheep on pastures contaminated with parasitic nematodes. Int J Parasitol. 1992, 22 (4): 417-425. 10.1016/0020-7519(92)90142-8.
Sreter T, Kassai T, Takacs E: The heritability and specificity of responsiveness to infection with Haemonchus contortus in sheep. Int J Parasitol. 1994, 24 (6): 871-876. 10.1016/0020-7519(94)90013-2.
Gruner L, Bouix J, Brunel JC: High genetic correlation between resistance to Haemonchus contortus and to Trichostrongylus colubriformis in INRA 401 sheep. Vet Parasitol. 2004, 119 (1): 51-58. 10.1016/j.vetpar.2003.10.014.
Kaplan RM: Drug resistance in nematodes of veterinary importance: a status report. Trends Parasitol. 2004, 20 (10): 477-481. 10.1016/j.pt.2004.08.001.
Knox DP, Redmond DL, Newlands GF, Skuce PJ, Pettit D, Smith WD: The nature and prospects for gut membrane proteins as vaccine candidates for Haemonchus contortus and other ruminant trichostrongyloids. Int J Parasitol. 2003, 33 (11): 1129-1137. 10.1016/S0020-7519(03)00167-X.
Meeusen EN, Piedrafita D: Exploiting natural immunity to helminth parasites for the development of veterinary vaccines. Int J Parasitol. 2003, 33 (11): 1285-1290. 10.1016/S0020-7519(03)00162-0.
Newton SE, Meeusen EN: Progress and new technologies for developing vaccines against gastrointestinal nematode parasites of sheep. Parasite Immunol. 2003, 25 (5): 283-296. 10.1046/j.1365-3024.2003.00631.x.
McEwan JC, Dodds KG, Watson TG, Greer GJ, Hosking BC, Douch PGC: Selection for resistance to roundworms by the New Zealand sheep breeding industry; the WormFec service. Proc Aust Assoc Anim Breed Genet. 1995, 11: 70-73.
Albers GA, Gray GD: Breeding for worm resistance: a perspective. Int J Parasitol. 1987, 17 (2): 559-566. 10.1016/0020-7519(87)90132-9.
Windon RG: Selective breeding for the control of nematodiasis in sheep. Rev Sci Tech. 1990, 9 (2): 555-576.
Bisset SA, Morris CA, McEwan JC, Vlassoff A: Breeding sheep in New Zealand that are less reliant on anthelmintics to maintain health and productivity. N Z Vet J. 2001, 49 (6): 236-246.
McEwan JC, Bisset SA, Morris CA: The selection of sheep for natural resistance to internal parasites. Sustainable Control of Internal Parasites in Ruminants. Edited by: G.K B. 1997, Canterbury, New Zealand , Lincoln University Press, 161-182.
Woolaston RRPLR: Selection of Merino sheep for resistance to Haemonchus contortus: genetic variation. Anim Sci. 1996, 62: 451-460.
Morris CA, Wheeler M, Watson TG, Hosking BC, Leathwick DM: Direct and correlated responses to selection for high or low faecal nematode egg count in Perendale sheep. N Z J Agric Res. 2005, 48: 1-10.
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 (1): 59-69. 10.1152/physiolgenomics.00257.2004.
Okada H, Kimura MT, Tan D, Fujiwara K, Igarashi J, Makuuchi M, Hui AM, Tsurumaru M, Nagase H: Frequent trefoil factor 3 (TFF3) overexpression and promoter hypomethylation in mouse and human hepatocellular carcinomas. Int J Oncol. 2005, 26 (2): 369-377.
Bohe H, Bohe M, Lundberg E, Polling A, Ohlsson K: Production and secretion of pancreatic secretory trypsin inhibitor in normal human small intestine. J Gastroenterol. 1997, 32 (5): 623-627.
Marchbank T, Freeman TC, Playford RJ: Human pancreatic secretory trypsin inhibitor. Distribution, actions and possible role in mucosal integrity and repair. Digestion. 1998, 59 (3): 167-174. 10.1159/000007485.
Krucken J, Epe M, Benten WP, Falkenroth N, Wunderlich F, Schroetel RM, Muller IU, Saidani N, Marinovski P, Stamm O: Malaria-suppressible expression of the anti-apoptotic triple GTPase mGIMAP8: Comparative analysis of the human gimap gene cluster encoding a novel GTPase family. J Cell Biochem. 2005, 341: 291-304.
Krucken J, Schroetel RM, Muller IU, Saidani N, Marinovski P, Benten WP, Stamm O, Wunderlich F: Comparative analysis of the human gimap gene cluster encoding a novel GTPase family. Gene. 2004, 341: 291-304. 10.1016/j.gene.2004.07.005.
Hosack DA, Dennis GJ, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.
Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.
Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, Kloos DU, Land S, Lewicki-Potapov B, Michael H, Munch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31 (1): 374-378. 10.1093/nar/gkg108.
Johnson K, Shapiro-Shelef M, Tunyaplin C, Calame K: Regulatory events in early and late B-cell differentiation. Mol Immunol. 2005, 42 (7): 749-761. 10.1016/j.molimm.2004.06.039.
Hong JH, Hwang ES, McManus MT, Amsterdam A, Tian Y, Kalmukova R, Mueller E, Benjamin T, Spiegelman BM, Sharp PA, Hopkins N, Yaffe MB: TAZ, a transcriptional modulator of mesenchymal stem cell differentiation. Science. 2005, 309 (5737): 1074-1078. 10.1126/science.1110955.
Dubuquoy L, Dharancy S, Nutten S, Pettersson S, Auwerx J, Desreumaux P: Role of peroxisome proliferator-activated receptor gamma and retinoid X receptor heterodimer in hepatogastroenterological diseases. Lancet. 2002, 360 (9343): 1410-1418. 10.1016/S0140-6736(02)11395-X.
Koretz K, Momburg F, Otto HF, Moller P: Sequential induction of MHC antigens on autochthonous cells of ileum affected by Crohn's disease. Am J Pathol. 1987, 129 (3): 493-502.
Arnold SJ, Osvath SR, Hall RA, King NJ, Sedger LM: Regulation of antigen processing and presentation molecules in West Nile virus-infected human skin fibroblasts. Virology. 2004, 324 (2): 286-296. 10.1016/j.virol.2004.03.036.
Johnson DR: Locus-specific constitutive and cytokine-induced HLA class I gene expression. J Immunol. 2003, 170 (4): 1894-1902.
Young LC, Thulien KJ, Campbell MR, Tron VA, Andrew SE: DNA mismatch repair proteins promote apoptosis and suppress tumorigenesis in response to UVB irradiation: an in vivo study. Carcinogenesis. 2004, 25 (10): 1821-1827. 10.1093/carcin/bgh191.
Mark SC, Sandercock LE, Luchman HA, Baross A, Edelmann W, Jirik FR: Elevated mutant frequencies and predominance of G:C to A:T transition mutations in Msh6(-/-) small intestinal epithelium. Oncogene. 2002, 21 (46): 7126-7130. 10.1038/sj.onc.1205861.
Li Z, Xia L, Lee LM, Khaletskiy A, Wang J, Wong JY, Li JJ: Effector genes altered in MCF-7 human breast cancer cells after exposure to fractionated ionizing radiation. Radiat Res. 2001, 155 (4): 543-553.
Kelner MJ, Bagnell RD, Montoya MA, Lanham KA: Structural organization of the human gastrointestinal glutathione peroxidase (GPX2) promoter and 3'-nontranscribed region: transcriptional response to exogenous redox agents. Gene. 2000, 248 (1-2): 109-116. 10.1016/S0378-1119(00)00137-2.
Chen J, Shpall RL, Meyerdierks A, Hagemeier M, Bottger EC, Naumovski L: Interferon-inducible Myc/STAT-interacting protein Nmi associates with IFP 35 into a high molecular mass complex and inhibits proteasome-mediated degradation of IFP 35. J Biol Chem. 2000, 275 (46): 36278-36284. 10.1074/jbc.M006975200.
Raasi S, Schmidtke G, Groettrup M: The ubiquitin-like protein FAT10 forms covalent conjugates and induces apoptosis. J Biol Chem. 2001, 276 (38): 35334-35343. 10.1074/jbc.M105139200.
Hipp MS, Kalveram B, Raasi S, Groettrup M, Schmidtke G: FAT10, a ubiquitin-independent signal for proteasomal degradation. Mol Cell Biol. 2005, 25 (9): 3483-3491. 10.1128/MCB.25.9.3483-3491.2005.
Liu D, Scafidi J, Prada AE, Zahedi K, Davis AE: Nuclear phosphatases and the proteasome in suppression of STAT1 activity in hepatocytes. Biochem Biophys Res Commun. 2002, 299 (4): 574-580. 10.1016/S0006-291X(02)02694-3.
Cai S, Davis AE: Complement regulatory protein C1 inhibitor binds to selectins and interferes with endothelial-leukocyte adhesion. J Immunol. 2003, 171 (9): 4786-4791.
Caliezi C, Wuillemin WA, Zeerleder S, Redondo M, Eisele B, Hack CE: C1-Esterase inhibitor: an anti-inflammatory agent and its potential use in the treatment of diseases other than hereditary angioedema. Pharmacol Rev. 2000, 52 (1): 91-112.
Liu D, Cramer CC, Scafidi J, Davis AE: N-linked glycosylation at Asn3 and the positively charged residues within the amino-terminal domain of the c1 inhibitor are required for interaction of the C1 Inhibitor with Salmonella enterica serovar typhimurium lipopolysaccharide and lipid A. Infect Immun. 2005, 73 (8): 4478-4487. 10.1128/IAI.73.8.4478-4487.2005.
Davis AE: Biological effects of C1 inhibitor. Drug News Perspect. 2004, 17 (7): 439-446. 10.1358/dnp.2004.17.7.863703.
Lafontaine PO, Arnal M, Buron N, Solary E, Bron AM, Westley BR, May FE, Bara J, Gespach C, Creuzot-Garcher C: Trefoil factor family mRNA and protein expression in pterygium. Int J Oncol. 2005, 27 (4): 997-1003.
Tornaletti S, Maeda LS, Kolodner RD, Hanawalt PC: Effect of 8-oxoguanine on transcription elongation by T7 RNA polymerase and mammalian RNA polymerase II. DNA Repair (Amst). 2004, 3 (5): 483-494. 10.1016/j.dnarep.2004.01.003.
Chu FF, Esworthy RS, Chu PG, Longmate JA, Huycke MM, Wilczynski S, Doroshow JH: Bacteria-induced intestinal cancer in mice with disrupted Gpx1 and Gpx2 genes. Cancer Res. 2004, 64 (3): 962-968. 10.1158/0008-5472.CAN-03-2272.
Lee YM, Cheng PY, Hong SF, Chen SY, Lam KK, Sheu JR, Yen MH: Oxidative stress induces vascular heme oxygenase-1 expression in ovariectomized rats. Free Radic Biol Med. 2005, 39 (1): 108-117. 10.1016/j.freeradbiomed.2005.02.033.
Kimura H, Hokari R, Miura S, Shigematsu T, Hirokawa M, Akiba Y, Kurose I, Higuchi H, Fujimori H, Tsuzuki Y, Serizawa H, Ishii H: Increased expression of an inducible isoform of nitric oxide synthase and the formation of peroxynitrite in colonic mucosa of patients with active ulcerative colitis. Gut. 1998, 42 (2): 180-187.
Christmann M, Tomicic MT, Kaina B: Phosphorylation of mismatch repair proteins MSH2 and MSH6 affecting MutSalpha mismatch-binding activity. Nucleic Acids Res. 2002, 30 (9): 1959-1966. 10.1093/nar/30.9.1959.
Puthalakath H, Villunger A, O'Reilly LA, Beaumont JG, Coultas L, Cheney RE, Huang DC, Strasser A: Bmf: a proapoptotic BH3-only protein regulated by interaction with the myosin V actin motor complex, activated by anoikis. Science. 2001, 293 (5536): 1829-1832. 10.1126/science.1062257.
Sax JK, El-Deiry WS: Identification and characterization of the cytoplasmic protein TRAF4 as a p53-regulated proapoptotic gene. J Biol Chem. 2003, 278 (38): 36435-36444. 10.1074/jbc.M303191200.
Puppo F, Contini P, Ghio M, Indiveri F: Soluble HLA class I molecules/CD8 ligation trigger apoptosis of CD8+ cells by Fas/Fas-ligand interaction. ScientificWorldJournal. 2002, 2: 421-423. 10.1100/tsw.2002.122.
Coux O, Tanaka K, Goldberg AL: Structure and functions of the 20S and 26S proteasomes. Annu Rev Biochem. 1996, 65: 801-847. 10.1146/annurev.bi.65.070196.004101.
Liu S, Yu Y, An H, Li N, Lin N, Wang W, Zhang W, Wan T, Cao X: Cloning and identification of a novel ubiquitin-like protein, BMSC-UbP, from human bone marrow stromal cells. Immunol Lett. 2003, 86 (2): 169-175. 10.1016/S0165-2478(03)00004-X.
Kalies KU, Allan S, Sergeyenko T, Kroger H, Romisch K: The protein translocation channel binds proteasomes to the endoplasmic reticulum membrane. Embo J. 2005, 24 (13): 2284-2293. 10.1038/sj.emboj.7600731.
Roberts AW, Kim C, Zhen L, Lowe JB, Kapur R, Petryniak B, Spaetti A, Pollock JD, Borneo JB, Bradford GB, Atkinson SJ, Dinauer MC, Williams DA: Deficiency of the hematopoietic cell-specific Rho family GTPase Rac2 is characterized by abnormalities in neutrophil function and host defense. Immunity. 1999, 10 (2): 183-196. 10.1016/S1074-7613(00)80019-9.
Carstanjen D, Yamauchi A, Koornneef A, Zang H, Filippi MD, Harris C, Towe J, Atkinson S, Zheng Y, Dinauer MC, Williams DA: Rac2 regulates neutrophil chemotaxis, superoxide production, and myeloid colony formation through multiple distinct effector pathways. J Immunol. 2005, 174 (8): 4613-4620.
Forlow SB, Schurr JR, Kolls JK, Bagby GJ, Schwarzenberger PO, Ley K: Increased granulopoiesis through interleukin-17 and granulocyte colony-stimulating factor in leukocyte adhesion molecule-deficient mice. Blood. 2001, 98 (12): 3309-3314. 10.1182/blood.V98.12.3309.
Lee SH, Corry DB: Homing alone? CD18 in infectious and allergic disease. Trends Mol Med. 2004, 10 (6): 258-262. 10.1016/j.molmed.2004.04.002.
Miyazaki T, Shen M, Fujikura D, Tosa N, Kim HR, Kon S, Uede T, Reed JC: Functional role of death-associated protein 3 (DAP3) in anoikis. J Biol Chem. 2004, 279 (43): 44667-44672. 10.1074/jbc.M408101200.
Dempsey PW, Doyle SE, He JQ, Cheng G: The signaling adaptors and pathways activated by TNF superfamily. Cytokine Growth Factor Rev. 2003, 14 (3-4): 193-209. 10.1016/S1359-6101(03)00021-2.
Bovine Genome Project[http://www.hgsc.bcm.tmc.edu/projects/bovine/].
Baird D, Johnstone P, Wilson T: Normalization of microarray data using a spatial mixed model analysis which includes splines. Bioinformatics. 2004, 20 (17): 3196-3205. 10.1093/bioinformatics/bth384.
Smyth GK: Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments. Stat Appl Genet Mol Biol. 2004, 3 (1): 1-25.
Michael JR: The stabilized probability plot. Biometrika. 1983, 70: 11-17.
Gene Expression Omnibus[http://www.ncbi.nlm.nih.gov/geo/].
Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1006/jmbi.1990.9999.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12 (6): 996-1006. 10.1101/gr.229102. Article published online before print in May 2002.
Ma XT, Qian MP, Tang HX: Predicting polymerase II core promoters by cooperating transcription factor binding sites in eukaryotic genes. Acta Biochim Biophys Sin (Shanghai). 2004, 36 (4): 250-258.
Aerts S, Thijs G, Coessens B, Staes M, Moreau Y, De Moor B: Toucan: deciphering the cis-regulatory logic of coregulated genes. Nucleic Acids Res. 2003, 31 (6): 1753-1764. 10.1093/nar/gkg268.
Bailey TL, Gribskov M: Combining evidence using p-values: application to sequence homology searches. Bioinformatics. 1998, 14 (1): 48-54. 10.1093/bioinformatics/14.1.48.
We thank Roger Wheeler, Wendy Bain, Chris Morris, Kevin Knowler, Mary Wheeler and Nadia McLean for assistance with maintenance of the Perendale selection lines, animal handling and slaughter. This work was funded by the Ovita Consortium Ltd.
OMK extracted all RNA, carried out the microarray experiments, Northern blotting and Gene Ontology analysis as well as preparing the manuscript. AZ carried out the MEME and MAST analysis. TW and DLH designed, constructed and validated the ovine cDNA arrays. GJG was responsible for all animal care, handling and experimentation. DBB developed the microarray normalisation procedure and carried out statistical analysis, AFM curated the EST database, assembled the contigs and provided bioinformatics support. AMC and JCM conceived the study, participated it its design, coordination, analysis and writing. All authors read and approved the final manuscript.