Expulsion of Trichuris muris is associated with increased expression of angiogenin 4 in the gut and increased acidity of mucins within the goblet cell
BMC Genomics volume 10, Article number: 492 (2009)
Trichuris muris in the mouse is an invaluable model for infection of man with the gastrointestinal nematode Trichuris trichiura. Three T. muris isolates have been studied, the Edinburgh (E), the Japan (J) and the Sobreda (S) isolates. The S isolate survives to chronicity within the C57BL/6 host whereas E and J are expelled prior to reaching fecundity. How the S isolate survives so successfully in its host is unclear.
Microarray analysis was used as a tool to identify genes whose expression could determine the differences in expulsion kinetics between the E and S T. muris isolates. Clear differences in gene expression profiles were evident as early as day 7 post-infection (p.i.). 43 probe sets associated with immune and defence responses were up-regulated in gut tissue from an E isolate-infected C57BL/6 mouse compared to tissue from an S isolate infection, including the message for the anti-microbial protein, angiogenin 4 (Ang4). This led to the identification of distinct differences in the goblet cell phenotype post-infection with the two isolates.
Differences in gene expression levels identified between the S and E-infected mice early during infection have furthered our knowledge of how the S isolate persists for longer than the E isolate in the C57BL/6 mouse. Potential new targets for manipulation in order to aid expulsion have been identified. Further we provide evidence for a potential new marker involving the acidity of the mucins within the goblet cell which may predict outcome of infection within days of parasite exposure.
Studies of Trichuris muris focus on one particular isolate, the Edinburgh (E) isolate. Infections of resistant mice, such as BALB/c, with the E isolate results in a protective Th2 response. Susceptibility to infection is associated with the host mounting a IFN-γ-dominated response inappropriate for worm expulsion and this is seen in mouse strains such as AKR [1–3].
However, other laboratory isolates of T. muris exist, the Japan (J) isolate, sub-cultured from the E isolate and the Sobreda (S) isolate discovered in Portugal. Interestingly, these isolates provoke different immune responses within the same host [4, 5], such as the C57BL/6 mouse strain. In the C57BL/6 mouse the S isolate is able to survive to chronicity, whereas the other two isolate are expelled prior day 21 p.i [6, 7]. This therefore, gives us a rare opportunity to study innate and adaptive immune responses to T. muris in the context of a resistant or susceptible outcome within one mouse strain without altering worm burden levels. The only other such model available involves manipulating egg dose to generate high or low dose infections and thus resistance (Th2) or susceptibility (Th1) . It has been previously reported that the S isolate survival, within the C57BL/6 mouse, is associated with a dampened effector Th2 response and an increased Th1 responses [6, 7].
Little is known however, about the underlying mechanisms evolved by the S isolate to enhance its survival within the host. Data from our laboratory suggests that the S isolate has evolved methods of manipulating the hosts T regulatory cell arm of the immune response, and the responses of key antigen presenting cells to parasite antigens. However differences in gene expression locally in gut tissue p.i. have not been analysed, despite the fact that they may underlie subsequent infection outcome.
Microarray analysis is a useful tool to look at global gene expression changes and indeed has been utilised to usefully inform research in many infections and diseases including Helicobacter pylori infection [11, 12] and inflammatory bowel disease . Interestingly, data already published from our laboratory has used oligonucleotide microarrays to determine gene expression changes in either resistant (BALB/c) or susceptible (AKR) mice infected with the E isolate of T. muris at 19 or 60 days post infection (p.i). Results indicated that AKR mice had a Th1- dominated mucosa, with up-regulated expression of genes associated with IFN-γ and BALB/c mice up-regulated the expression of genes coding for potential anti-parasitic proteins including intelectin and angiogenins . The association of IFN-γ with susceptibility during a T. muris infection has also been shown via reverse transcription (RT)-PCR analyses  and its functional importance revealed by blocking studies .
Here we analyse gut tissue, from C57BL/6 mice infected with either the E isolate or the S isolate. A time point of 7 days p.i was chosen, as it is a time point where the host will not have expelled either of the isolates (D'Elia et al. unpublished data), yet early gene expression changes may be occurring which later determine whether the parasite is expelled or not. Data presented here highlights the possible importance of two genes in particular, indoleamine 2,3-dioxygenase (INDO) and angiogenin 4 (Ang4). INDO is a gene whose expression is up-regulated by both isolates and may aid early survival within the host. Ang4 is a gene whose expression is differentially regulated by the two isolates and may determine outcome of infection within the C57BL/6 host. Further, histological analysis of goblet cells demonstrated differences in their mucin content p.i with the E or S isolates of T. muris. Changes in Ang4 expression and/or other genes described here may explain the alterations of mucin acidity within the goblet cell and these changes may contribute to S isolate survival within the C57BL/6 host.
Enhanced survival of the S isolate of T. muris compared to the E isolate in C57BL/6 mice
As previously reported , the S isolate of T. muris is able to persist for longer and indeed reach fecundity in a C57BL/6 mouse unlike the E isolate of T. muris. As shown in Fig. 1, worm expulsion is clearly evident at d21 p.i in an E isolate infection with only two mice harbouring high number of worms. In contrast, all five mice infected with the S isolate have high numbers of worms at d21 p.i and still have worms present at d42 p.i where as, at this time point, all five C57BL/6 mice have expelled the E isolate (p = 0.0079, Mann-Whitney test).
Commonality and differences in RNA expression profiles following infection with the E or S isolates of T. muris
Mice infected with either the E or S isolates of T. muris were sacrificed at d7 p.i. Microarray analysis was carried out at this time point, prior to expulsion of either isolates to identify changes in gene expression which might reflect the early immune interactions which underlie the subsequent differences in survival. Microarray data from RNA extracted from three gut samples (naïve C57BL/6 mice, E isolate-infected C57BL/6 mice d7 p.i and S isolate-infected C57BL/6 mice d7 p.i), done in triplicate, was analysed by the PUMA Bayesian statistical method . The probability of positive log-ratio (PPLR) statistical scores were generated for each probe set by two-way comparisons between the three sample groups. PPLR values range from 0 to 1, with values closest to 0 representing the most significantly down-regulated probe sets and values closest to 1, representing the most significantly up-regulated probe sets. Values of 0.5 represent no significant change [17, 18]. Venn diagrams were generated to show the distribution of the probe sets. Probe sets with a probability of positive log-ratio (PPLR) value greater than 0.999 (significantly up-regulated) in at least one of the comparisons (probe set expression of gut tissue from an E isolate infection versus that of naïve gut tissue [EvN], probe set expression of gut tissue from an S isolate infection versus that of naïve gut tissue [SvN] and probe set expression of gut tissue from an S isolate infection versus that of gut tissue from an E isolate infection [SvE]) indicate that 35 probe sets are uniquely up-regulated when comparing the gene expression profile of gut tissue from a naïve mouse to the gene expression profile of gut tissue from an E isolate-infected mouse (EvN). When comparing SvN, 14 probe sets are unique, interestingly however, 3 probe sets are up-regulated in both an E isolate infection and an S isolate infection. Importantly there were 17 probe sets that were uniquely up-regulated when comparing the gene expression profiles of gut tissue taken from mice infected with the two isolates - SvE (Fig. 2A). The number of probe sets associated with a significantly down-regulated PPLR value (0.001) was higher than the number for up-regulated probe sets. 168 probe sets alone were significantly down-regulated in an S isolate infection compared to an E isolate infection- SvE (Fig. 2B).
Combining probe sets that were either significantly up-regulated (0.999) or significantly down-regulated (0.001) in at least one of the comparisons (EvN, SvN and SvE) generated a list of 409 probe sets. 164 of these probe sets were unique to SvE (Fig. 2C and Table 1G). A total of 29 genes were significantly up-regulated or down-regulated in gut tissue p.i with the more commonly studied E isolate compared to an uninfected (Naïve) mouse (Fig. 2C and Table 1A). Only 3 genes, however, were common in all three comparisons (EvN, SvN and SvE) (Fig. 2C and Table 1E). Examples of the 409 probe sets are shown in Table 1, with the table split into sections representing the segments shown in Fig. 2E (full list of probe sets available in additional data file 1).
A scatter plot was generated to show the general profile of the microarray data set. Log fold change is shown for S isolate infection over naïve (y axis) vs. log fold change of E isolate infection over naïve (x axis). The majority of probe sets fit within the centre quadrant, representing probe sets with small changes between the isolates. The four quadrants on the extremities of the scatter plot identify the probe sets that may be of interest, i.e. up or down in both the E isolate-infected gut and the S isolate-infected gut compared to naive or up in one isolate infection and down in the other isolate compared to naïve (Fig 2D).
Microarray probe sets associated with immune response respond similarly in the gut of an S isolate-infected mouse to an un-infected (naïve) mouse
The 409 probe sets that had a PPLR value greater than 0.999 or less than 0.001 in at least one of the comparisons (EvN, SvN and SvE) underwent clustering analysis (Fig. 3A). Six distinct clusters were generated and probe sets from each cluster were subjected to expression analysis systematic explorer (EASE) online tool and DAVID [19, 20]. For each cluster, overrepresented gene ontology (GO) groups were identified (Fig 3B).
Cluster A contained 59 probe sets associated with mating behaviour but respond in a similar pattern for both an E isolate infection and an S isolate infection. In cluster B, the S isolate down-regulates 122 probe sets associated with cell cycle and cytoskeletal organisation compared to naïve and E isolate probe sets. In cluster C, containing 43 probe sets, the S isolates expression levels in the gut are similar to naïve levels, whereas the E isolate over express these probe sets. Interestingly, these 43 probe sets are associated with GO groups for immune response, defence response, antigen presentation and response to stress. In cluster D, the S isolate infection causes an increase in gene expression for probe sets associated with growth regulation but down-regulates genes associated with muscle development in cluster E. In both clusters D and E, the E isolate and naïve samples are responding similarly. The final cluster, F, shows no change in naïve levels, a down-regulation with an E isolate infection and an up regulation with the S isolate. These 39 probe sets were associated with potassium ion transport (Fig. 3B). Cluster C is expanded and shown in detail in Fig. 3C. Genes include Angiogenin 4 (Ang4), indoleamine-pyrrole 2,3 dioxygenase (INDO), chemokine ligand 8 (CCL8) and interferon inducible GTPase 1 (Iigp1).
Confirmation of genes from cluster C; angiogenin 4 and indoleamine-pyrrole 2,3 dioxygenase
Given that we had previously demonstrated a correlation between Ang4 expression and resistance and INDO expression and susceptibility in a study using different strains of mouse infected with the E isolate  we selected these two genes to corroborate our microarray data by qPCR and immunohistochemistry.
INDO was a highly expressed gene for both E and S isolate compared to naïve (PPLR value = 0.99999, 0.998805 respectively). However when comparing SvE the PPLR value implied that the E isolate transcribed more INDO than the S isolate (PPLR = 0.057528). These findings were confirmed by qPCR on RNA extracted from the gut and the epithelial cell fraction (Fig. 4A) (Whole gut - E isolate = 6.68 ± 1.88, S isolate = 3.65 ± 3.4; pooled epithelial cells - E isolate = 20.03, S isolate 10.14).
Although qPCR revealed a non-significant elevation in INDO in the gut and epithelial cell fraction of mice infected with the E isolate compared to the S isolate, analysis of IDO at the protein level revealed significantly higher numbers of positive staining cells in an E isolate infection. Thus, infection with either isolate resulted in significantly increased levels of IDO+ cells compared to naïve (E, p < 0.01; S, p < 0.05), with the E isolate infection provoking higher numbers compared to infection with the S isolate (p < 0.05) (Fig. 4B and 4C).
By microarray, gene expression of Ang4 was higher in the gut of mice infected with the E isolate compared to the levels seen in the gut of S isolate-infected C57BL/6 mice (PPLR value = 3.84E-06, Table 1). Ang4 expression was also elevated in whole gut and epithelial cell fractions derived from E-infected mice compared to S-infected mice as identified by qPCR (Fig. 4D). Further, immunohistochemical staining revealed significantly more Ang4+ cells at day 21 p.i. in E isolate-infected mice (Fig. 4E and 4F).
Increased numbers and acidity of goblet cells in C57BL/6 mice infected with the E isolate compared to the S isolate of T. muris
Quantification of Ang4 expression has previously focused on small intestinal Paneth cells as the cellular source . Fig. 4D shows higher expression of Ang4 in the epithelial cell fraction and in Fig. 4F this expression can clearly be seen to localise to the goblet cell suggesting that as no Paneth cells are present in the large intestine, the goblet cell is the most likely source.
Therefore, goblet cell staining was carried out on proximal colon sections from naïve, E isolate-infected and S isolate-infected C57BL/6 mice. Total numbers of goblet cells per 20 colonic crypt units (ccu) is shown throughout an infection time course (Fig. 5A). At day 21 p.i an E isolate-infected mice had significantly higher number of goblet cells compared to mice infected with the S isolate (p < 0.05 Mann Whitney U test). Following infection with either the E or S isolate total numbers of goblet cells increased steadily over the full time course from naïve levels.
An interesting phenotypical change in the acidity of mucins was within the goblet cells was evident between the infections with the different isolates. At d7 p.i mice infected with the S isolate of T. muris had significantly higher numbers of goblet cells containing neutral (pink) mucins compared to mice infected with the E isolate (Fig. 5B) (p < 0.01 Mann Whitney U test). When these values are converted to percentages of total goblet cell numbers, the pattern is exaggerated at d7 and d13 p.i (Fig. 5C). Examples of typical sections from mice harbouring either the E isolate or S isolate, stained for goblet cells, taken on d7 p.i are shown in Fig. 5D at two magnifications (×200 and ×400) clearly highlight the acidic mucin differences in the goblet cells.
Acidity and/or neutrality of goblet cells are determined by the addition or removal of sugars on mucins by transferases. Increased transferase function would cause increased acidity of mucins within the goblet cells. Table 2 indicates possible transferases that could be causing the phenotypical changes of the goblet cell mucins seen between the two isolates at d7 p.i from the whole gut microarray analyses. Thus, analyses of gene expression levels in whole gut tissue from S and E infected mice reveals a down regulation of all four transferases in tissue from mice infected with the S isolate compared to the E isolate.
The intimate interactions between hosts and parasites have been studied for many years . The ability of helminths, in particular T. muris to affect the immune environment within the host has a broad associated literature [23–25]. Genome wide approaches to understand the changes occurring in the host at a gene expression level have been utilised in many infectious models and diseases [26, 27]. Indeed microarray and RT-PCR analysis have been used to explore the responses of mouse strains infected with T. muris [14, 15]. However, these experiments use different mouse strains in order to compare the gene expression profiles associated with resistance or susceptibility. Here we use one mouse strain, C57BL/6, and different T. muris isolates to achieve susceptibility or resistant outcomes. In addition we look early in infection to identify changes in gene expression, that could later determine susceptibility or resistance of the host.
We demonstrate that the S isolate of T. muris, which is able to persist chronically in a C57BL/6 host, provokes a different gene expression profile in gut tissue compared to a C57BL/6 mouse infected with the E isolate, which is ultimately expelled. These changes are already apparent by d7 p.i. prior to E isolate expulsion. In particular, differences in 43 probe sets associated with GO groups - immune response, defence response, antigen presentation and response to stress were identified.
Interestingly the expression of genes from an S isolate infected gut resemble those seen in a naïve (uninfected) gut, whereas there was an elevation in the expression of the genes represented by the 43 probe sets during an E isolate infection. This data would imply that even at d7 p.i., the S isolate is dampening down expression of genes that would normally be associated with expulsion of the parasite. Microarray analyses on gut tissue from different strains of mice infected with the E isolate had previously identified differences in the expression of INDO and Ang4. As these two genes were also identified as being differentially expressed in this current study we investigated their expression further.
Indoleamine 2,3-dioxygenase plays an important part in the kynurenine pathway where it metabolises tryptophan which has been associated with controlling parasite growth in a number of infections [28–30]. IDO also has immunoregulatory roles [14, 31, 32], indeed its ability to dampen cell proliferation  may extend to the control of gut epithelial cell turnover and thus delayed or absent worm expulsion given the importance of epithelial cell turnover in the expulsion of T. muris . As the expression of IDO was elevated in gut tissue from both E and S isolate-infected mice its presence may reflect a less efficient ability, rather than a complete inability, to expel the parasite.
In contrast, Ang4 expression did differentiate between worm survival and expulsion with expression levels only elevated in E-infected mice. The angiogenins are a family of closely related proteins described as Paneth cell-derived and encoded for, in the mouse, by a gene cluster on chromosome 1 . They belong to the RNase superfamily  which includes eosinophil secretory granules proteins which are toxic to some gastrointestinal nematodes . Originally implicated in the growth of tumours, the normal physiological role of the angiogenins was unclear until the demonstration of a role for Paneth cell derived Ang4 as an endogenous anti-microbial protein central in epithelial host defence against gut-dwelling bacteria . Immunohistochemical staining of gut tissue in this study localised Ang4 production to the goblet cell, identifying the cell type as a novel cellular source in the large intestine. This interesting discovery expands on the already extensive literature on the importance of goblet cells in the expulsion of nematode infection, where increased numbers of goblet cells correlate with resistance [37, 38].
A goblet cell hyperplasia is reported in the context of many nematode infections including T. muris, T. spiralis, H. polygyrus and N. brasiliensis and this is thought to be under the control of a Th2 response [37–39]. Several goblet cell factors have been identified to play an important role in nematode infection including Relm-β and Muc2. Relm-β expression is linked with the production of Th2 cytokines and worm expulsion in T. muris, T. spiralis and N. brasiliensis . Thus effective worm expulsion may be achieved through combined effects of antimicrobial proteins including the angiogenins, Relm-β and indeed the intelectins [14, 41].
Goblet cell numbers did not differ quantitatively between E and S isolate infection until d21 p.i. However, the goblet cells were qualitatively strikingly different at d7 p.i. Thus the number of neutral goblet cells present in the gut of an S isolate infected C57BL/6 mouse are significantly higher than that seen in an E isolate infected mouse. By d21 the number of neutral goblet cells is similar for both isolates, if not higher in the E isolate infection. At this time point the majority of E isolate worms have been expelled.
The acidity of goblet cells is determined by the addition and removal of sugars to the mucins within the goblet cell . Transferases are required to do this and Table 2 highlights candidate transferases that may explain the differences seen in the quality of the goblet cells. Increased gene expression of the transferases is seen in the microarray data in the gut of an E isolate-infected mouse compared to an S isolate-infected gut. The increases in gene expression of the transferases listed could explain the increased acidity of the goblet cell mucins seen in an E isolate infection and may even underlie subsequent expulsion.
The C57BL/6 mouse strain differs in its ability to expel the E and S isolates of T. muris. Microarray analysis of gut tissue taken from infected mice just one week p.i. revealed significant differences in gene expression profiles. Thus, under similar levels of parasite exposure in a single mouse strain changes which ultimately correlate with resistance or susceptibility are apparent.
We report an association between the expression of the anti-microbial protein Ang4 and resistance to infection and further identify the goblet cell as a novel cellular source. Qualitative differences in goblet cell mucins early p.i. are also reported and may act as an early predictor of future resistance or susceptibility.
C57BL6 male mice were obtained from Harlan (Bicester, UK) at 6-8 weeks of age. Mice were specific pathogen free and maintained in sterile conditions in individually ventilated cages by the Biological Services Faculty (BSF), University of Manchester, UK. All work was performed under the regulations of the Home Office Scientific Procedures Act (1986).
Both the E and S isolate of T. muris were maintained as previously described . Excretory/secretory (E/S) antigen was collected by incubating adult worms in RPMI media for 4 hrs. Mice were infected orally with around 150-200 embryonated eggs. Worm burden analysis was carried out as previously described .
C57BL/6 mice were infected with ~200 embryonated eggs from the E isolate or the S isolate of T. muris. 7 days p.i. with the E or S isolate of T. muris, groups of 5 mice were killed and gut samples were taken from each mouse, including naïve (uninfected) mice. Samples were placed in trizol and snap frozen in liquid nitrogen for RNA extraction or placed in Neutral buffer formalin (NBF) for immunohistochemistry.
Two groups of five mice infected with either the E isolate or the S isolate of T. muris were allowed to progress to day 14 p.i to determine worm burden. Three independent infection experiments were run.
Isolation of intestinal epithelial cells
As previously described, . Briefly, the caecum and approximately 5 cm of colon were removed. The tissue was then slit longitudinally and rinsed in calcium- and magnesium-free Hanks balanced salt solution containing 2% foetal calf serum (FCS) (CMF2%), cut into 1-cm pieces, and placed into ice-cold CMF2%. Samples were washed until supernatant was clear. The tissue was then placed into calcium- and magnesium-free Hanks balanced salt solution containing 10% FCS, 1 mM EDTA, 1 mM dithiothreitol, 100 units/ml penicillin, and 100 μg/ml streptomycin (CMF10%) and incubated at 37°C for 20 min. The supernatant was then passed through a 100-μm cell strainer (Becton Dickinson, Oxford, United Kingdom) and centrifuged at 200 × g for 10 min. The cells were resuspended in ice-cold RPMI 1640 (Invitrogen, Carlsbad, CA), 5% Foetal calf serum (FCS), 100 U/ml penicillin and 100 μg/ml streptomycin, 1% L-glutamine, 0.1% MTG (Sigma Aldrich, St Louis, MO). The cell concentration determined using a CASY-1 Coulter Counter. The final cell concentration was adjusted to give 5 × 106 cells/ml. An aliquot was then taken for RNA extraction.
Gut Samples were removed from -80°C freezer and allowed to thaw on ice. Each sample was placed into a sterile falcon tube and homogenised using an Ultra -Turrax (T25 Basic) homogeniser. The homogeniser was cleaned with chloroform, ethanol and 3× DEPC water (Promega) between each sample. The samples were then placed in RNAse-free 1.5 ml eppendorfs. Following tissue homogenisation, cell debris was pelleted by spinning at 13000 rpm for 10 minutes at 4°C. The supernatant was then recovered into a new tube. The samples were incubated at room temperature for 5 minutes. 0.2 ml of chloroform (Sigma Aldrich) was added to the supernatant and samples were shook vigorously for 15 seconds. Samples were then centrifuged for 15 minutes at 4°C at 13000 rpm. The aqueous phase was then transferred to a tube containing 0.5 ml isopropyl alcohol (Sigma Aldrich) and incubated at RT for 10 min. Samples were then centrifuged at 13000 rpm for 10 min at 4°C. The resulting supernatant was removed and 1 ml of 75% ethanol was added. Samples were mixed by vortexing and centrifuged at 12000 rpm for 5 min at 4°C. This was repeated once more. The supernatant was removed and the pellet was air dried for 15-20 minutes. Pellets were resuspended in 50 μl of nuclease free water (Promega). Samples were then incubated at 60°C for 10 minutes. RNA was pooled from five individual animals within a group and each group run as an individual microarray chip. Three independent infection experiments were run and therefore a total of 9 microarray chips were run and analysed.
A turbo DNase kit (Ambion Turbo DNA-free) was used to remove all pieces of contaminating DNA. To each 50 μl of RNA, 5 μl of Turbo DNase buffer and 1 μl of DNase enzyme was added. Samples were then incubated at 37°C for 30 minutes. Following incubation, 5 μl of turbo DNase inactivation buffer was added and left for 2 minutes. Samples were then spun down for 1 minute at 13000 g. The supernatant was then placed into a clean 200 μl eppendorf and stored at -80°C. The concentration of RNA was determined via Nanodrop spectrometry (Nanodrop technologies). The quality of RNA was checked on a 2% agarose gel in 0.5× TBE buffer (Invitrogen) with 0.01% Ethidium Bromide. Samples were loaded with Orange G loading buffer (4% sucrose in dH20 + 0.01% Orange G).
RNA quality was checked using the RNA 6000 Nano Assay, and analyzed on an Agilent 2100 Bioanalyser (Agilent Technologies). RNA was quantified using a Nanodrop ultra-low-volume spectrophotometer (Nanodrop Technologies).
Hybridization cocktail was hybridised to HG-U133 PLUS2 oligonucleotide arrays (Affymetrix) according to manufacturer's instructions. Arrays were read using Agilent GeneArray scanner 3000 7G using Affymetrix GCOS (V1.4) software.
Technical quality control was performed with dChip (V2005) http://www.dchip.org Normalisation and expression analysis was done using multi-mgmos . Differential expression between the sample groups (E, S, N) was assessed with a Bayesian method which includes probe-level measurement error when assessing statistical significance . Analysis was performed with the PUMA package in R [17, 47].
A gene list of differentially expressed genes (409 probe sets) was created by filtering for probe sets with a probability of positive log-ratio (PPLR) value less than 0.001 (or greater than 0.999) in any of the three-way comparisons between E, S and N samples. This data set was segregated into 6 clusters based on similarity of expression profile across the dataset using a k-means clustering algorithm. Clustering was performed on the means of each sample group (log 2) that had been z-transformed (for each probe set the mean set to zero, standard deviation to 1). K-means clustering was done on the basis of similarity of profiles across the dataset using the "Super Grouper" plug-in of maxdView software (available from http://bioinf.man.ac.uk/microarray/maxd/)). Functional annotation of the genes was performed using DAVID version 2. http://david.abcc.ncifcrf.gov/) .
Microarray data has been submitted in a MIAME compliant standard to the Array Express database (Experiment E-MEXP-1795, http://www.ebi.ac.uk/microarray-as/ae/)
Quantitative real-time polymerase chain reaction (qPCR)
1.0 μg of total RNA was reverse transcribed using Bioscript (Bioline, London, U.K.) in a final volume of 40 μl according to the manufacturer's instructions and stored at -20° until used. Quantitative PCR was performed using SensiMix plus SYBR (Bioline) on an OPTICON DNA engine with OPTICON MONITOR software version 2·03 (Real-Time systems; MJ Research, Hemel Hempstead, UK). Amplification of mRNA encoding 18S was performed to control for the starting amount of cDNA. Expression levels of genes of interest are shown as fold change over that seen in naïve animals after normalization to housekeeping gene levels using the ΔΔCt method. Primers sequences were AGTCCCTGCCTTTGTACACA and GATCCGAGGGCCTCACTAAC for 18S, CTGCACGACATAGCTACCAGTCTG and ACATTTGAGGGCTCTTCCGACTTG for IDO and CTCTGGCTCAGAATGTAAGGTACGA and GAAATCTTTAAAGGCTCGGTACCC for Ang4. All sequences are 5'-3' with the sense primer given first.
Histological sections were prepared from proximal colon tissue fixed in 10% buffered formalin and embedded in paraffin. 6-μm sections were cut using a microtome and added to gelatin-coated microscope slides. Sections were dewaxed with citroclear and taken to water through decreasing concentrations of ethanol.
For IDO and Ang4 staining slides were washed in PBS and endogenous peroxidase activity was quenched using 0.064 mg/ml sodium azide, 1.5 U/ml glucose oxidase, and 1.8 mg/ml D-glucose (Sigma-Aldrich) in PBS for 20 min at 37°C. Endogenous avidin and biotin binding sites were blocked using a commercial kit according to the manufacturer's instructions (Vector Laboratories). For IDO staining this was followed by a mouse on mouse kit, M.O.M (Vector Laboratories). Briefly, slides were incubated for 1 hr in mouse Ig blocking reagent, 5 minutes in working solution and 30 minutes in primary antibody (MAB5412, MS × Indoleamine 2,3-dioxygenase). For Ang4 staining slides were incubated consecutively with 1.5% donkey serum (Jackson Immunoresearch), sheep anti-mouse Ang4 and biotinylated F(ab')2 donkey anti-sheep IgG (Jackson Immunoresearch). Colour development in both cases was by incubation with ABC (Avidin Biotin Complex), DAB (3,3'-diaminobenzidine) and haematoxylin (all from Vector Laboratories).
For goblet cell staining, mucins in goblet cells were stained with 1% alcian blue (Sigma) in 3% acetic acid, washed and treated with 1% periodic acid (Sigma) followed by counterstaining in Mayer's haematoxylin (Sigma). Slides were dehydrated and mounted in aquamount (BDH Laboratory Supplies, Poole, UK). For enumeration of immunohistochemistry and goblet cell staining, the average number of cells from 20 crypts was taken from three different sections per mouse.
Differences between groups were tester by non-parametric methods, the Mann-Whitney U test (two factor comparisons) and Kruskal Wallis test (for more than two factors) with Dunns post test for multiple parameter comparisons. All tests were performed using GraphPad Prism software (San Diego, CA)
expression analysis systematic explorer
- INDO :
Indoleamine 2,3-dioxygenase gene
Indoleamine 2,3-dioxygenase protein
minimum information about a microarray experiment
probability of positive log-ratio
quantitative real-time PCR.
Else K, Wakelin D: The effects of H-2 and non-H-2 genes on the expulsion of the nematode Trichuris muris from inbred and congenic mice. Parasitology. 1988, 96 (Pt 3): 543-550. 10.1017/S0031182000080173.
Else KJ, Grencis RK: Helper T-cell subsets in mouse trichuriasis. Parasitology today (Personal ed). 1991, 7 (11): 313-316. 10.1016/0169-4758(91)90268-S.
Else KJ, Hultner L, Grencis RK: Cellular immune responses to the murine nematode parasite Trichuris muris. II. Differential induction of TH-cell subsets in resistant versus susceptible mice. Immunology. 1992, 75 (2): 232-237.
Bellaby T, Robinson K, Wakelin D, Behnke JM: Isolates of Trichuris muris vary in their ability to elicit protective immune responses to infection in mice. Parasitology. 1995, 111 (Pt 3): 353-357. 10.1017/S0031182000081907.
Koyama K, Ito Y: Comparative studies on immune responses to infection in susceptible B10.BR mice infected with different strains of the murine nematode parasite Trichuris muris. Parasite immunology. 1996, 18 (5): 257-263. 10.1046/j.1365-3024.1996.d01-92.x.
Bellaby T, Robinson K, Wakelin D: Induction of differential T-helper-cell responses in mice infected with variants of the parasitic nematode Trichuris muris. Infection and immunity. 1996, 64 (3): 791-795.
Johnston CE, Bradley JE, Behnke JM, Matthews KR, Else KJ: Isolates of Trichuris muris elicit different adaptive immune responses in their murine host. Parasite immunology. 2005, 27 (3): 69-78. 10.1111/j.1365-3024.2005.00746.x.
Bancroft AJ, Else KJ, Grencis RK: Low-level infection with Trichuris muris significantly affects the polarization of the CD4 response. European journal of immunology. 1994, 24 (12): 3113-3118. 10.1002/eji.1830241230.
D'Elia R, Behnke JM, Bradley JE, Else KJ: Regulatory T cells: a role in the control of helminth-driven intestinal pathology and worm survival. J Immunol. 2009, 182 (4): 2340-2348. 10.4049/jimmunol.0802767.
D'Elia R, Else KJ: In vitro antigen presenting cell-derived IL-10 and IL-6 correlate with Trichuris muris isolate-specific survival. Parasite immunology. 2009, 31 (3): 123-131. 10.1111/j.1365-3024.2008.01088.x.
Mueller A, O'Rourke J, Chu P, Kim CC, Sutton P, Lee A, Falkow S: Protective immunity against Helicobacter is characterized by a unique transcriptional signature. Proc Natl Acad Sci USA. 2003, 100 (21): 12289-12294. 10.1073/pnas.1635231100.
Mueller A, O'Rourke J, Grimm J, Guillemin K, Dixon MF, Lee A, Falkow S: Distinct gene expression profiles characterize the histopathological stages of disease in Helicobacter-induced mucosa-associated lymphoid tissue lymphoma. Proc Natl Acad Sci USA. 2003, 100 (3): 1292-1297. 10.1073/pnas.242741699.
Dieckgraefe BK, Stenson WF, Korzenik JR, Swanson PE, Harrington CA: Analysis of mucosal gene expression in inflammatory bowel disease by parallel oligonucleotide arrays. Physiological genomics. 2000, 4 (1): 1-11.
Datta R, deSchoolmeester ML, Hedeler C, Paton NW, Brass AM, Else KJ: Identification of novel genes in intestinal tissue that are regulated after infection with an intestinal nematode parasite. Infection and immunity. 2005, 73 (7): 4025-4033. 10.1128/IAI.73.7.4025-4033.2005.
Schopf LR, Hoffmann KF, Cheever AW, Urban JF, Wynn TA: IL-10 is critical for host resistance and survival during gastrointestinal helminth infection. J Immunol. 2002, 168 (5): 2383-2392.
Else KJ, Finkelman FD, Maliszewski CR, Grencis RK: Cytokine-mediated regulation of chronic intestinal helminth infection. The Journal of experimental medicine. 1994, 179 (1): 347-351. 10.1084/jem.179.1.347.
Pearson RD, Liu X, Sanguinetti G, Milo M, Lawrence ND, Rattray M: puma: a Bioconductor package for propagating uncertainty in microarray analysis. BMC bioinformatics. 2009, 10: 211-10.1186/1471-2105-10-211.
Liu X, Milo M, Lawrence ND, Rattray M: A tractable probabilistic model for Affymetrix probe-level analysis across multiple chips. Bioinformatics (Oxford, England). 2005, 21 (18): 3637-3644. 10.1093/bioinformatics/bti583.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome biology. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome biology. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.
Hooper LV, Stappenbeck TS, Hong CV, Gordon JI: Angiogenins: a new class of microbicidal proteins involved in innate immunity. Nature immunology. 2003, 4 (3): 269-273. 10.1038/ni888.
Gupta S: Parasite immune escape: new views into host-parasite interactions. Current opinion in microbiology. 2005, 8 (4): 428-433. 10.1016/j.mib.2005.06.011.
Deschoolmeester ML, Else KJ: Cytokine and chemokine responses underlying acute and chronic Trichuris muris infection. International reviews of immunology. 2002, 21 (4-5): 439-467. 10.1080/08830180213278.
Cliffe LJ, Grencis RK: The Trichuris muris system: a paradigm of resistance and susceptibility to intestinal nematode infection. Advances in parasitology. 2004, 57: 255-307. full_text.
Else KJ, deSchoolmeester ML: Immunity to Trichuris muris in the laboratory mouse. Journal of helminthology. 2003, 77 (2): 95-98. 10.1079/JOH2002162.
Ng HH, Frantz CE, Rausch L, Fairchild DC, Shimon J, Riccio E, Smith S, Mirsalis JC: Gene expression profiling of mouse host response to Listeria monocytogenes infection. Genomics. 2005, 86 (6): 657-667. 10.1016/j.ygeno.2005.07.005.
Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, De La Vega P, Holder AA, Batalov S, Carucci DJ, et al: Discovery of gene function by expression profiling of the malaria parasite life cycle. Science (New York, NY). 2003, 301 (5639): 1503-1508.
Fujigaki S, Saito K, Takemura M, Maekawa N, Yamada Y, Wada H, Seishima M: L-tryptophan-L-kynurenine pathway metabolism accelerated by Toxoplasma gondii infection is abolished in gamma interferon-gene-deficient mice: cross-regulation between inducible nitric oxide synthase and indoleamine-2,3-dioxygenase. Infection and immunity. 2002, 70 (2): 779-786. 10.1128/IAI.70.2.779-786.2002.
Grohmann U, Fallarino F, Puccetti P: Tolerance, DCs and tryptophan: much ado about IDO. Trends in immunology. 2003, 24 (5): 242-248. 10.1016/S1471-4906(03)00072-3.
Silva NM, Rodrigues CV, Santoro MM, Reis LF, Alvarez-Leite JI, Gazzinelli RT: Expression of indoleamine 2,3-dioxygenase, tryptophan degradation, and kynurenine formation during in vivo infection with Toxoplasma gondii: induction by endogenous gamma interferon and requirement of interferon regulatory factor 1. Infection and immunity. 2002, 70 (2): 859-868. 10.1128/IAI.70.2.859-868.2002.
Munn DH, Shafizadeh E, Attwood JT, Bondarev I, Pashine A, Mellor AL: Inhibition of T cell proliferation by macrophage tryptophan catabolism. The Journal of experimental medicine. 1999, 189 (9): 1363-1372. 10.1084/jem.189.9.1363.
Fallarino F, Grohmann U, Hwang KW, Orabona C, Vacca C, Bianchi R, Belladonna ML, Fioretti MC, Alegre ML, Puccetti P: Modulation of tryptophan catabolism by regulatory T cells. Nature immunology. 2003, 4 (12): 1206-1212. 10.1038/ni1003.
Cliffe LJ, Humphreys NE, Lane TE, Potten CS, Booth C, Grencis RK: Accelerated intestinal epithelial cell turnover: A new mechanism of parasite expulsion. Science. 2005, 308 (5727): 1463-1465. 10.1126/science.1108661.
Strydom DJ: The angiogenins. Cell Mol Life Sci. 1998, 54 (8): 811-824. 10.1007/s000180050210.
Holloway DE, Hares MC, Shapiro R, Subramanian V, Acharya KR: High-Level Expression of Three Members of the Murine Angiogenin Family in Escherichia coli and Purification of the Recombinant Proteins. Protein Expression and Purification. 2001, 22 (2): 307-317. 10.1006/prep.2001.1434.
Hamann KJ, Barker RL, Loegering DA, Gleich GJ: Comparative toxicity of purified human eosinophil granule proteins for newborn larvae of Trichinella spiralis. The Journal of parasitology. 1987, 73 (3): 523-529. 10.2307/3282130.
Else KJ, Finkelman FD: Intestinal nematode parasites, cytokines and effector mechanisms. International journal for parasitology. 1998, 28 (8): 1145-1158. 10.1016/S0020-7519(98)00087-3.
Bancroft AJ, Grencis RK: Th1 and Th2 cells and immunity to intestinal helminths. Chemical immunology. 1998, 71: 192-208. full_text.
Yamauchi J, Kawai Y, Yamada M, Uchikawa R, Tegoshi T, Arizono N: Altered expression of goblet cell- and mucin glycosylation-related genes in the intestinal epithelium during infection with the nematode Nippostrongylus brasiliensis in rat. Apmis. 2006, 114 (4): 270-278. 10.1111/j.1600-0463.2006.apm_353.x.
Artis D, Wang ML, Keilbaugh SA, He W, Brenes M, Swain GP, Knight PA, Donaldson DD, Lazar MA, Miller HR, et al: RELMbeta/FIZZ2 is a goblet cell-specific immune-effector molecule in the gastrointestinal tract. Proc Natl Acad Sci USA. 2004, 101 (37): 13596-13600. 10.1073/pnas.0404034101.
Pemberton AD, Knight PA, Gamble J, Colledge WH, Lee JK, Pierce M, Miller HR: Innate BALB/c enteric epithelial responses to Trichinella spiralis: inducible expression of a novel goblet cell lectin, intelectin-2, and its natural deletion in C57BL/10 mice. J Immunol. 2004, 173 (3): 1894-1901.
Ishikawa N, Wakelin D, Mahida YR: Role of T helper 2 cells in intestinal goblet cell hyperplasia in mice infected with Trichinella spiralis. Gastroenterology. 1997, 113 (2): 542-549. 10.1053/gast.1997.v113.pm9247474.
Wakelin D: Acquired immunity to Trichuris muris in the albino laboratory mouse. Parasitology. 1967, 57 (3): 515-524. 10.1017/S0031182000072395.
Else KJ, Wakelin D: Genetically-determined influences on the ability of poor responder mice to respond to immunization against Trichuris muris. Parasitology. 1990, 100 (Pt 3): 479-489. 10.1017/S0031182000078793.
Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98 (1): 31-36. 10.1073/pnas.011404098.
Liu X, Milo M, Lawrence ND, Rattray M: Probe-level measurement error improves accuracy in detecting differential gene expression. Bioinformatics (Oxford, England). 2006, 22 (17): 2107-2113. 10.1093/bioinformatics/btl361.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome biology. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
RD is supported by a BBSRC PhD studentship. This work was also supported by BBSRC research grant number BB/E012647/1. We thank L. Wardleworth from the University of Manchester Core Microarray Facility. We also thank Dr. Lora Hooper for providing the Angiogenin 4 plasmid and Professors Jerzy Behnke and Jan Bradley for providing the S isolate of T. muris to KJE. The authors declare no conflict of interest.
RD was involved in study design, performed experiments and data analysis, and was involved in manuscript preparation. MLD performed experiments and data analysis and was involved manuscript preparation. LZ was involved in microarray analysis and manuscript preparation. SHW and AP produced the sheep anti-mouse Ang4 antibody. KJE conceived of the study and participated in its design and coordination and helped to draft the manuscript. Both RD and KJE have had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Riccardo D'Elia, Matthew L deSchoolmeester contributed equally to this work.
Electronic supplementary material
Additional file 1: Significantly altered gene probe set expression in the gut following infection with the E or S isolate of T. muris. The fold-change values for all 409 probe sets identified with a PPLR value greater than 0.999 or less than 0.001 in at least one of the comparisons (EvN, SvN and SvE) are presented. (XLS 119 KB)
About this article
Cite this article
D'Elia, R., deSchoolmeester, M.L., Zeef, L.A. et al. Expulsion of Trichuris muris is associated with increased expression of angiogenin 4 in the gut and increased acidity of mucins within the goblet cell. BMC Genomics 10, 492 (2009). https://doi.org/10.1186/1471-2164-10-492