Expulsion of Trichuris muris is associated with increased expression of angiogenin 4 in the gut and increased acidity of mucins within the goblet cell

  • Riccardo D'Elia1, 3Email author,

    Affiliated with

    • Matthew L deSchoolmeester1Email author,

      Affiliated with

      • Leo AH Zeef1,

        Affiliated with

        • Steven H Wright2,

          Affiliated with

          • Alan D Pemberton2 and

            Affiliated with

            • Kathryn J Else1

              Affiliated with

              BMC Genomics200910:492

              DOI: 10.1186/1471-2164-10-492

              Received: 13 October 2008

              Accepted: 24 October 2009

              Published: 24 October 2009

              Abstract

              Background

              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.

              Results

              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.

              Conclusion

              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.

              Background

              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 [13].

              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) [8]. 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[9], and the responses of key antigen presenting cells to parasite antigens[10]. 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 [13]. 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 [14]. The association of IFN-γ with susceptibility during a T. muris infection has also been shown via reverse transcription (RT)-PCR analyses [15] and its functional importance revealed by blocking studies [16].

              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.

              Results

              Enhanced survival of the S isolate of T. muris compared to the E isolate in C57BL/6 mice

              As previously reported [7], 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).
              http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-492/MediaObjects/12864_2008_Article_2376_Fig1_HTML.jpg
              Figure 1

              E and S isolate infection in C57BL/6 mice. C57BL/6 mice were infected with ~200 embryonated eggs from either the E or the S isolate of T. muris. At days 14, 21 and 42 p.i. the mice were sacrificed, caecae and proximal colons were removed and the number of worms counted. Each symbol represents an individual animal (5 mice per group) and the line represents the mean. **, significant difference between the E isolate and S isolate-infected mice (Mann Whitney U test, p < 0.01).

              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 [17]. 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).
              http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-492/MediaObjects/12864_2008_Article_2376_Fig2_HTML.jpg
              Figure 2

              Analysis of microarray data of Trichuris infected gut tissue with either the E isolate or S isolate shown in Venn Diagrams and scatter plots. PPLR values were calculated for E isolated-infected gut tissue compared to naïve (uninfected) gut tissue (EvN), S isolated-infected gut tissue compared to naïve gut tissue (SvN) and S isolated-infected gut tissue compared to E isolate-infected gut tissue (SvE). Distribution of probe sets with a PPLR value greater than 0.999 (A), PPLR less than 0.001 (B) and PPLR greater than 0.999 and/or less than 0.001 (C) in any one of the three comparisons (EvN, SvN or SvE) are shown in Venn diagrams. The scatter plot (D) shows the average log2-fold changes in probe set expression in gut tissue from E isolate-infected C57BL/6 mice at day 7 compared to that from uninfected C57BL/6 mice (x axis) versus the average log2-fold changes in probe set expression in gut tissue from S isolate C57BL/6 mice at day 7 compared to that from uninfected C57BL/6 mice (y axis). Gut tissue was pooled from five individual animals within a group, and three independent infection experiments were run. The average log2-fold changes presented represent values calculated from all three infection experiments.

              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).
              Table 1

              Selection of probe sets from each section of the Venn diagram with a PPLR greater than 0.999 and/or less than 0.001.

                 

              PPLR Value

              Probe Set ID

              Gene Symbol

              Gene Description

              E v N

              S v N

              S v E

              A. Probe sets of PPLR > 0.999 or < 0.001 in E v N only

                 

              1417655_a_at

              Ars2

              arsenate_resistance_protein_2

              0.999501

              0.887027

              0.005792

              1418652_at

              Cxcl9

              chemokine_(C-X-C_motif)_ligand_9

              0.999702

              0.997352

              0.083513

              1449102_at

              Ebf2

              early_B-cell_factor_2

              9.14E-04

              0.001533

              0.670498

              1435906_x_at

              Gbp2

              guanylate_nucleotide_binding_protein_2

              0.999957

              0.920883

              0.003407

              1419060_at

              Gzmb

              granzyme_B

              0.999274

              0.834789

              0.004865

              1425477_x_at

              H2-Ab1

              histocompatibility_2,_class_II_antigen_A,_beta_1

              0.999653

              0.917835

              0.009449

              1420437_at

              Indo

              indoleamine-pyrrole_2,3_dioxygenase

              0.999996

              0.998805

              0.057528

              1447927_at

              Mpa2l

              macrophage_activation_2_like

              0.999643

              0.721523

              0.002485

              B. Probe sets of PPLR > 0.999 or < 0.001 in E v N and S v N

                 

              1425519_a_at

              Cd74

              CD74_antigen_(invariant_polypeptide_of_major_histo

              compatibility_complex,_class_II_antigen-associated)

              1

              0.999649

              0.002733

              1456612_at

              Il17d

              Interleukin_17D

              4.12E-06

              4.19E-06

              0.535351

              1427615_at

              Itga4

              integrin_alpha_4

              8.56E-04

              5.59E-04

              0.422541

              1428056_at

              Klra7

              Killer_cell_lectin-like_receptor,_subfamily_A,_member_7

              5.32E-04

              6.98E-04

              0.561618

              1430132_at

              Krt28

              keratin_28

              9.00E-04

              2.22E-04

              0.244397

              1441873_at

              Prlpe

              prolactin-like_protein_E

              3.25E-05

              3.65E-05

              0.605783

              1428111_at

              Slc38a4

              solute_carrier_family_38,_member_4

              1.51E-04

              3.85E-04

              0.6009

              1429901_at

              Tcba1

              T-cell_lymphoma_breakpoint_associated_target_1

              2.14E-04

              1.20E-04

              0.438909

              C. Probe sets of PPLR > 0.999 or < 0.001 in S v N only

                 

              1419549_at

              Arg1

              arginase_1,_liver

              0.23588

              5.87E-04

              0.003669

              1425532_a_at

              Bin1

              bridging_integrator_1

              0.985409

              0.999463

              0.941873

              1449637_at

              Cdh4

              cadherin_4

              0.001517

              4.54E-04

              0.293429

              1422542_at

              Gpr34

              G_protein-coupled_receptor_34

              0.0022

              5.97E-04

              0.51351

              1427850_x_at

              Igh-VJ558

              immunoglobulin_heavy_chain_(J558_family)

              0.819554

              0.999188

              0.995084

              1422250_at

              Map3k2

              mitogen_activated_protein_kinase_kinase_kinase_2

              0.790984

              0.999799

              0.998933

              1417051_at

              Pcdh8

              protocadherin_8

              0.704933

              0.999094

              0.997809

              1416576_at

              Socs3

              suppressor_of_cytokine_signaling_3

              0.927402

              0.999427

              0.989809

              D. Probe sets of PPLR > 0.999 or < 0.001 in E v N and S v E

                 

              1440832_at

              Ang4

              angiogenin,_ribonuclease_A_family,_member_4

              0.999347

              0.119804

              3.84E-06

              1419684_at

              Ccl8

              chemokine_(C-C_motif)_ligand_8

              0.999999

              0.98578

              4.66E-04

              1447620_at

              Cno

              cappuccino

              0.999021

              0.473507

              4.52E-04

              1426731_at

              Des

              desmin

              0.999964

              0.270284

              1.93E-06

              1435290_x_at

              H2-Aa

              histocompatibility_2,_class_II_antigen_A,_alpha

              1

              0.928626

              4.61E-08

              1419043_a_at

              Iigp1

              interferon_inducible_GTPase_1

              0.999969

              0.644236

              2.52E-04

              1449490_at

              Mbd4

              methyl-CpG_binding_domain_protein_4

              0.999147

              0.449218

              1.50E-04

              1423505_at

              Tagln

              transgelin

              0.999831

              0.031314

              1.35E-08

              E. Probe sets of PPLR > 0.999 or < 0.001 in E v N and S v E and S v N

                 

              1455106_a_at

              Ckb

              creatine_kinase,_brain

              0.999225

              6.03E-04

              2.32E-10

              1418240_at

              Gbp2

              guanylate_nucleotide_binding_protein_2

              1

              0.999995

              1.49E-07

              1449009_at

              Tgtp

              T-cell_specific_GTPase

              1

              0.999989

              4.83E-11

              F. Probe sets of PPLR > 0.999 or < 0.001 in S v N and S v E

                 

              1425103_at

              Ace2

              angiotensin_I_converting_enzyme_(peptidyl-

              dipeptidase_A)_2

              0.275223

              0.999142

              0.999784

              1416330_at

              Cd81

              CD_81_antigen

              0.599285

              4.84E-06

              1.05E-06

              1457644_s_at

              Cxcl1

              chemokine_(C-X-C_motif)_ligand_1

              0.378865

              0.999467

              0.999866

              1416685_s_at

              Fbl

              fibrillarin

              0.873786

              1.01E-06

              1.11E-09

              1416221_at

              Fstl1

              follistatin-like_1

              0.716963

              1.49E-04

              1.05E-05

              1415988_at

              Hdlbp

              high_density_lipoprotein_(HDL)_binding_protein

              0.531698

              5.56E-04

              4.14E-04

              1447692_x_at

              Map3k7ip1

              mitogen-activated_protein_kinase_kinase_kinase

              _7_interacting_protein_1

              0.321227

              0.999609

              0.999863

              1450641_at

              Vim

              vimentin

              0.727751

              5.04E-05

              2.09E-06

              G. Probe sets of PPLR > 0.999 or < 0.001 in S v E only

                 

              1443341_at

              Anxa4

              Annexin_A4

              0.994522

              0.054991

              3.27E-06

              1460218_at

              Cd52

              CD52_antigen

              0.979891

              0.167083

              9.70E-04

              1418457_at

              Cxcl14

              chemokine_(C-X-C_motif)_ligand_14

              0.931634

              0.008835

              1.23E-04

              1422412_x_at

              Ear3

              eosinophil-associated,_ribonuclease_A_family,_

              member_3

              0.947348

              0.011739

              3.24E-04

              1423407_a_at

              Fbln2

              fibulin_2

              0.919939

              0.025244

              2.66E-04

              1457346_at

              Foxp1

              Forkhead_box_P1

              0.998173

              0.287033

              3.68E-04

              1427351_s_at

              Igh-6

              immunoglobulin_heavy_chain_6_(heavy_chain_of_IgM)

              0.974353

              0.029092

              6.25E-05

              1454665_at

              Irf2bp2_///_LOC672960

              interferon_regulatory_factor_2_binding_protein_2_///_

              similar_to_interferon_regulatory_factor_2_binding_

              protein_2

              0.660352

              0.003366

              8.80E-04

              1415961_at

              Itm2c

              integral_membrane_protein_2C

              0.699876

              0.003609

              4.85E-04

              1451407_at

              Jam4

              junction_adhesion_molecule_4

              0.130394

              0.98967

              0.999784

              1449989_at

              Mcpt2

              mast_cell_protease_2

              0.997228

              0.27691

              3.33E-04

              1443037_at

              Nptn

              Neuroplastin

              0.971795

              0.087496

              5.36E-04

              1445827_at

              Prkcbp1

              Protein_kinase_C_binding_protein_1

              0.952881

              0.098262

              9.82E-04

              1456262_at

              Rbm5

              RNA_binding_motif_protein_5

              0.995795

              0.026747

              1.31E-06

              1417654_at

              Sdc4

              syndecan_4

              0.78476

              0.003683

              2.83E-04

              1447034_at

              Tnpo1

              Transportin_1

              0.932695

              0.022353

              2.24E-04

              A = Probe sets of PPLR > 0.999 or < 0.001 in E isolate-infected gut tissue compared to naïve gut tissue only, B = Probe sets of PPLR > 0.999 or < 0.001 in E isolate-infected gut tissue compared to naïve gut tissue and S isolate-infected gut tissue compared to naïve gut tissue, C = of PPLR > 0.999 or < 0.001 in S isolate-infected gut tissue compared to naïve gut tissue only, D = Probe sets of PPLR > 0.999 or < 0.001 in E isolate-infected gut tissue compared to naïve gut tissue and S isolate-infected gut tissue compared to E isolate-infected gut tissue, E = Probe sets of PPLR > 0.999 or < 0.001 in E isolate-infected gut tissue compared to naïve gut tissue and S isolate-infected gut tissue compared to E isolate-infected gut tissue and S isolate-infected gut tissue compared to naïve gut tissue, F = Probe sets of PPLR > 0.999 or < 0.001 in S isolate-infected gut tissue compared to naïve gut tissue and S isolate-infected gut tissue compared to E isolate-infected gut tissue, G = Probe sets of PPLR > 0.999 or < 0.001 in S isolate-infected gut tissue compared to E isolate-infected gut tissue only.

              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).
              http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-492/MediaObjects/12864_2008_Article_2376_Fig3_HTML.jpg
              Figure 3

              Clusters, profiles and Gene Ontology groupings of the 409 probe sets that had a PPLR value greater than 0.999 and/or less than 0.001 in at least one of the three comparisons (EvN, SvN or SvE). A: Flow chart showing method of generating the clusters from the probe sets. B: The 409 genes were assigned to one of six distinct clusters using k-means clustering algorithms. Data for each cluster are represented with a z-transformed Eisen colour plot with the number of genes in each cluster shown. Red and green indicate positive and negative change from zero, respectively, with colour intensity indicating the degree of deviation. Profiles of the z-transformed data (for each probe set, the mean set to 0 and SD to 1) for each of the experimental groups are shown next to the Eisen plot. The most significantly overrepresented Gene Ontology terms are then shown for the genes within each cluster. C: Expanded version of cluster C indicating the 43 probe sets, showing probe set number, name and symbol along with Eisen plot.

              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 [14] 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).
              http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-492/MediaObjects/12864_2008_Article_2376_Fig4_HTML.jpg
              Figure 4

              Quantification of gene expression change by qPCR and Immunohistichemistry for Angiogenin 4 and indoleamine 2,3-dioxygenase. Gene expression levels for INDO in whole gut and epithelial cells (A) comparing E isolate (white bars) to S isolate (black bars) as a fold change over naïve/unstimulated. IDO+ cells were detected in proximal colon tissue sections via immunohistochemistry and data presented as mean number of IDO+ cells per 20 ccu (B). Example of IDO staining (×200 magnification) in naïve, E isolate- and S isolate-infected mice at d7 p.i, with examples of positive cells identified by black arrowheads (C). Gene expression levels for Ang4 in whole gut and epithelial cell fraction (D) comparing E isolate-infected mice (white bars) to S isolate-infected mice (black bars) as a fold change over naïve. *, significant difference between the E isolate and S isolate-infected mice (p < 0.05). Ang4+ cells were detected in proximal colon tissue sections via immunohistochemistry and data presented as mean number of Ang4+ cells per 20 ccu (E). Example of Ang4 staining (×200 magnification) in naïve, E isolate- and S isolate-infected mice at d21 p.i (F). Experiments were performed at least twice and data shown is a representative experiment. Values represent the mean ± SD for 3-4 mice per group.

              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 [21]. 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.
              http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-10-492/MediaObjects/12864_2008_Article_2376_Fig5_HTML.jpg
              Figure 5

              Changes in number and acidity of goblet cells in E isolate-infected C57BL/6 mice compared to S isolate-infected C57BL/6 mice. C57BL/6 mice were infected with either the E or S isolates of Trichuris muris. At each time point mice infected with either the E or S isolate were sacrificed and goblet cells were detected in proximal colon tissue sections via immunohistochemistry. Data are presented as mean number of goblet cells per 20 ccu (A), mean number of neutral (stain pink following alcian blue and periodic acid treatment) goblet cells per 20 ccu (B) and percentage of neutral goblet cells (C). Examples of goblet cell staining (×200 and ×400 magnification) in E isolate-infected or S isolate-infected mice at d7 p.i (D). *; p < 0.05 and **; p < 0.01 for comparison between the E isolate and S isolate-infected mice (Mann Whitney U test). Experiments were performed at least twice and data shown is a representative experiment. Values represent the mean ± SD for five mice per group.

              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.
              Table 2

              Selections of probe sets which may explain differences associated with neutral goblet cells seen in S isolate-infected mice.

                 

              PPLR

              Fold Change

              Probe Set ID

              Gene Symbol

              Gene Description

              E v N

              S v N

              S v E

              S v E

              1420903_at

              St6galnac3

              ST6 (alpha-N-acetyl-neuraminyl-2,3-beta-galactosyl-1,3)-N-acetylgalactosaminide alpha-2,6-sialyltransferase 3

              0.7554

              0.1668

              0.0473

              -3.7792

              1425668_a_at

              St3gal4

              ST3 beta-galactoside alpha-2,3-sialyltransferase 4

              0.8429

              0.4391

              0.1049

              -1.4277

              1442929_at

              St8sia5

              ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 5

              0.7053

              0.2824

              0.1351

              -3.0470

              1443740_at

              St3gal6

              ST3 beta-galactoside alpha-2,3-sialyltransferase 6

              0.1933

              0.0122

              0.0629

              -6.5657

              PPLR values for all three comparisons (EvN, SvN or SvE) and fold change in gene expression levels of probe sets from the gut of an S isolate-infected mouse compared to expression levels in an E isolate-infected gut are shown along with probe set identification (ID) number, gene symbol and gene description.

              Discussion

              The intimate interactions between hosts and parasites have been studied for many years [22]. The ability of helminths, in particular T. muris to affect the immune environment within the host has a broad associated literature [2325]. 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 [2830]. IDO also has immunoregulatory roles [14, 31, 32], indeed its ability to dampen cell proliferation [31] 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 [33]. 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 [34]. They belong to the RNase superfamily [35] which includes eosinophil secretory granules proteins which are toxic to some gastrointestinal nematodes [36]. 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 [21]. 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 [3739]. 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 [40]. 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 [42]. 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.

              Conclusion

              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.

              Methods

              Mice

              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).

              Parasites

              Both the E and S isolate of T. muris were maintained as previously described [43]. 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 [44].

              Sample collection/protocol

              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, [14]. 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.

              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.

              DNase treatment

              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).

              Microarray procedure

              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.

              Microarray analysis

              Technical quality control was performed with dChip (V2005) http://​www.​dchip.​org[45] Normalisation and expression analysis was done using multi-mgmos [46]. 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 [18]. 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/​) [19].

              MIAME

              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.

              Histology

              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.

              Statistics

              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)

              List of abbreviations

              Ang4: 

              Angiogenin 4

              d: 

              day

              EASE: 

              expression analysis systematic explorer

              E/S: 

              Excretory/Secretory

              GO: 

              gene ontology

              INDO

              Indoleamine 2,3-dioxygenase gene

              IDO: 

              Indoleamine 2,3-dioxygenase protein

              MIAME: 

              minimum information about a microarray experiment

              p.i.: 

              post infection

              PPLR: 

              probability of positive log-ratio

              qPCR: 

              quantitative real-time PCR

              Declarations

              Acknowledgements

              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.

              Authors’ Affiliations

              (1)
              Faculty of Life Sciences, University of Manchester
              (2)
              Division of Veterinary Clinical Sciences, University of Edinburgh, Easter Bush Veterinary Centre
              (3)
              DSTL, Porton Down

              References

              1. 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.View ArticlePubMed
              2. Else KJ, Grencis RK: Helper T-cell subsets in mouse trichuriasis. Parasitology today (Personal ed) 1991,7(11):313–316.View Article
              3. 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.PubMed
              4. 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.View ArticlePubMed
              5. 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.View ArticlePubMed
              6. 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.PubMed
              7. 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.View ArticlePubMed
              8. 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.View ArticlePubMed
              9. 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.View ArticlePubMed
              10. 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.View ArticlePubMed
              11. 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.View ArticlePubMed
              12. 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.View ArticlePubMed
              13. 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.PubMed
              14. 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.View ArticlePubMed
              15. Schopf LR, Hoffmann KF, Cheever AW, Urban JF Jr, Wynn TA: IL-10 is critical for host resistance and survival during gastrointestinal helminth infection. J Immunol 2002,168(5):2383–2392.PubMed
              16. 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.View ArticlePubMed
              17. 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.View ArticlePubMed
              18. 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.View Article
              19. Dennis G Jr, 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.View ArticlePubMed
              20. Hosack DA, Dennis G Jr, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome biology 2003,4(10):R70.View ArticlePubMed
              21. 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.View ArticlePubMed
              22. Gupta S: Parasite immune escape: new views into host-parasite interactions. Current opinion in microbiology 2005,8(4):428–433.View ArticlePubMed
              23. 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.View ArticlePubMed
              24. 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.View ArticlePubMed
              25. Else KJ, deSchoolmeester ML: Immunity to Trichuris muris in the laboratory mouse. Journal of helminthology 2003,77(2):95–98.View ArticlePubMed
              26. 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.View ArticlePubMed
              27. 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.View Article
              28. 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.View ArticlePubMed
              29. Grohmann U, Fallarino F, Puccetti P: Tolerance, DCs and tryptophan: much ado about IDO. Trends in immunology 2003,24(5):242–248.View ArticlePubMed
              30. 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.View ArticlePubMed
              31. 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.View ArticlePubMed
              32. 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.View ArticlePubMed
              33. 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.View ArticlePubMed
              34. Strydom DJ: The angiogenins. Cell Mol Life Sci 1998,54(8):811–824.View ArticlePubMed
              35. 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.View ArticlePubMed
              36. 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.View ArticlePubMed
              37. Else KJ, Finkelman FD: Intestinal nematode parasites, cytokines and effector mechanisms. International journal for parasitology 1998,28(8):1145–1158.View ArticlePubMed
              38. Bancroft AJ, Grencis RK: Th1 and Th2 cells and immunity to intestinal helminths. Chemical immunology 1998, 71:192–208.View ArticlePubMed
              39. 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.View ArticlePubMed
              40. 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.View ArticlePubMed
              41. 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.PubMed
              42. 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.View ArticlePubMed
              43. Wakelin D: Acquired immunity to Trichuris muris in the albino laboratory mouse. Parasitology 1967,57(3):515–524.View ArticlePubMed
              44. 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.View ArticlePubMed
              45. 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.View ArticlePubMed
              46. 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.View Article
              47. 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.View ArticlePubMed

              Copyright

              © D'Elia et al. 2009

              This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

              Advertisement