Genetic analysis of the Trichuris muris-induced model of colitis reveals QTL overlap and a novel gene cluster for establishing colonic inflammation
- Scott E Levison1,
- Paul Fisher2,
- Jenny Hankinson3,
- Leo Zeef4,
- Steve Eyre5,
- William E Ollier3,
- John T McLaughlin1,
- Andy Brass6,
- Richard K Grencis†7 and
- Joanne L Pennock†1Email author
© Levison et al.; licensee BioMed Central Ltd. 2013
Received: 30 October 2012
Accepted: 14 February 2013
Published: 26 February 2013
Genetic susceptibility to colonic inflammation is poorly defined at the gene level. Although Genome Wide Association studies (GWAS) have identified loci in the human genome which confer susceptibility to Inflammatory Bowel Disease (Crohn’s and Ulcerative Colitis), it is not clear if precise loci exist which confer susceptibility to inflammation at specific locations within the gut e.g. small versus large intestine. Susceptibility loci for colitis in particular have been defined in the mouse, although specific candidate genes have not been identified to date. We have previously shown that infection with Trichuris muris (T. muris) induces chronic colitis in susceptible mouse strains with clinical, histological, and immunological homology to human colonic Crohn’s disease. We performed an integrative analysis of colitis susceptibility, using an F2 inter-cross of resistant (BALB/c) and susceptible (AKR) mice following T. muris infection. Quantitative Trait Loci (QTL), polymorphic and expression data were analysed alongside in silico workflow analyses to discover novel candidate genes central to the development and biology of chronic colitis.
7 autosomal QTL regions were associated with the establishment of chronic colitis following infection. 144 QTL genes had parental strain SNPs and significant gene expression changes in chronic colitis (expression fold-change ≥ +/-1.4). The T. muris QTL on chromosome 3 (Tm3) mapped to published QTL in 3 unrelated experimental models of colitis and contained 33 significantly transcribed polymorphic genes. Phenotypic pathway analysis, text mining and time-course qPCR replication highlighted several potential cis-QTL candidate genes in colitis susceptibility, including FcgR1, Ptpn22, RORc, and Vav3.
Genetic susceptibility to induced colonic mucosal inflammation in the mouse is conserved at Tm3 and overlays Cdcs1.1. Genes central to the maintenance of intestinal homeostasis reside within this locus, implicating several candidates in susceptibility to colonic inflammation. Combined methodology incorporating genetic, transcriptional and pathway data allowed identification of biologically relevant candidate genes, with Vav3 newly implicated as a colitis susceptibility gene of functional relevance.
KeywordsTrichuris muris Colitis Genetic susceptibility Cdcs1 Crohn’s
Many diseases result from the complex interaction of environmental and genetic factors (e.g. Crohn’s disease, diabetes mellitus) [1, 2]. Phenotypic expression is influenced by multiple genes, which individually may increase or decrease the probability of disease development. Gene variation and gene-gene interactions, additionally results in non-linear contributions to phenotypic variation. Discovering the genetic architecture of complex traits thus represents a true challenge  and requires collaborative multi-disciplinary investigation and a variety of experimental approaches [4, 5]. The exploration of new animal models of colitis with well-defined phenotypes and homology to human pathology, provide a comparative approach to refine biological discoveries for subsequent human translation.
Trichuris muris, a natural intestinal parasite of mice has been extensively studied as a model for human whipworm (Trichuris trichiura) infection. In dissecting the immune response to Trichuris infection, a paradigm of resistance and susceptibility to chronic colonic inflammation has emerged . Following the ingestion of parasite ova, acute colitis develops in all mice, but it is the genetic composition of mouse strain which dictates the presence of colitis. BALB/c mice mount immune-mediated TH2 dependent parasite expulsion (IL4, IL5 and IL13 expression) [6, 7] with full resolution within 20 days. Conversely AKR mice sustain a chronic Trichuris infection, respond with a TH1 immune response (IFNγ, IL12), and subsequent establishment of colitis . These polarized outcomes occur despite identical treatment and conditioning, and are almost certainly determined by host genetic variation. Importantly, we have recently characterised differences in colonic tissue transcription between susceptible and resistant mice and demonstrated phenotypic, immunological and biological pathway homology to human Crohn’s disease . These data present T. muris infection not as an aetiological factor in the pathogenesis of Crohn’s disease, nor solely a model of infection but as a viable and relevant colitis model to investigate and study mucosal inflammation.
The multifactorial and complex nature of Crohn’s disease remains to be fully characterised, but it is evident that disease can be initiated anywhere along the digestive tract. It is likely that precise environmental triggers determine the site of initiation, but it is also possible that host genetics play a part. A variety of experimental models have been developed to study pathogenic mechanisms responsible for the induction and perpetuation of Crohn’s disease. Phenotypic and biological factors common between colonic Crohn’s disease and chronic T. muris induced colitis, present a novel opportunity to characterise the genetic architecture central to disease susceptibility in the colon. The aim of the current study was to identify genome wide genetic elements and mechanistic pathways which underpin the development and maintenance of such chronic inflammation.
Systemic and colonic phenotyping of chronic T. muriscolitis in an F2 population of resistant and susceptible mice
Colonic histological assessment demonstrated persistent T. muris infection and large bowel inflammation. Mild-to-moderate inflammatory changes included: transmural tissue oedema and associated leukocytic infiltration (lymphocytes, macrophages, neutrophils); prominent mucosal and submucosal reactive lymphoid aggregates; colonic crypt hyperplasia and hypertrophy (Figure 1C-E). Significant correlation between histological parameters of inflammation (e.g. crypt length), immune response phenotype (Figure 1F: Spearman’s Rs = -0.54) and worm burden (Figure 1G: Spearman’s Rs = 0.84), were demonstrated. 98.5% of mice with persistent helminthosis demonstrated colonic inflammatory changes.
Whole genome Linkage analysis
Summary of Trichuris muris QTL ( Tm ) found across the genome
T. Muris QTL
Gene no. within QTL (Ensembl)
AKR vs BALB/c genes with SNPs
Serum IgG’s Worm count
Serum IgG’s Worm count
Serum IgG’s Worm count
Serum IgG’s Worm count
Serum IgG’s Worm count
Prioritization of QTL candidate genes via pathway-driven workflow analysis
The cross correlation of functional pathways containing QTL genes and genes demonstrating significant expression were identified, linking genotype and phenotype trait interactions (Figure 3.3). Finally, polymorphic genes between parental AKR and BALB/c mice were identified within each locus (Figure 3.4).
As an example, 344 Ensembl gene ID’s were detected within Tm3. Of these, 97 (28%) were designated as functionally important within molecular interaction networks, as assigned by the KEGG pathway database. Significantly expressed microarray genes were similarly assigned biological pathways. For Tm3, the cross correlation of common pathway data and the exclusion of any gene which lacked SNPs between parental strains, identified 17 Quantitative Trait genes (Figure 3.4, Column D). In comparison, 61 KEGG-assigned polymorphic genes did not demonstrate any change in transcriptional activity (Column C). Of the genes yet to be allocated a KEGG pathway, 16 of 191 genes displayed significant transcription (Column B). The same process was undertaken for all 7 QTL (Figure 3).
Chromosome 3 candidates
Significantly expressed Tm3 genes possessing strain-specific SNPs and a designated biological (KEGG) pathway
KEGG designated pathway
KEGG designated pathway
Chemokine signalling pathway
Fc gamma R-mediated phagocytosis
Fc gamma R-mediated phagocytosis
Hematopoietic cell lineage
Leukocyte transendothelial migration
T cell receptor signalling pathway
Systemic lupus erythematosus
B cell receptor signalling pathway
Regulation of actin cytoskeleton
Cell adhesion moleculaes (CAMs)
Leukocyte transendothelial migration
Synthesis and degradation of ketone
Terpenoid backbone biosynthesis
Valine, leuckine and isoleucine degradation
Antigen processing and presentation
Valine, leucine and isoleucine degradation
Drug metabolism – cytochrome P450
Metabolism of xenobiotics by cytochrome p450
Drug metabolism – cytochrome P450
Steroid hormone biosynthesis
Metabolism of xenobiotcis by cytochrome P450
Drug metabolism – cytochrome p450
Drug metabolism – other enzymes
Drug metabolism – cytochrome p450
Pantothenate and CoA biosynthesis
Systemic lupus erythematosus
Significantly expressed genes possessing strain-specific SNPs but as of yet, undesignated a biological (KEGG) pathway
Text mining (cosine vector score)
Candidate gene validation
Trichuris muris-induced colitis represents a tractable murine model for understanding the patho-biological mechanisms of chronic intestinal inflammation [9, 19]. The use of Quantitative Trait Loci (QTL) mapping based on continuous phenotypic variation has proved a useful technique in many murine polygenic traits including intestinal inflammation [12, 20, 21]. Yet, of more than 2000 QTL documented within the mouse genome database  fewer than 1% of studies have actually been characterized at a gene or molecular level, due to the small effect size of the susceptibility locus in question (<10% penetrance), or the large interval size defined . New multi-factorial approaches have been discussed in the literature  and demonstrate that understanding complex genetic traits requires an integrative analysis.
Specific steps were taken in our experimental design to consider recent reports concerning the QTL/microarray approach in the identification of QTL candidate genes . First, QTL were defined with regards to experimental phenotype (pQTL), and correlated with transcriptional expression activity in parental strains. Second, the use of high density Affymetrix exon array, which targets approximately 40 exonic probes per gene, overcame any problem of potential allelic-biased probe binding. Third, a hypothesis-free pathway analysis, backed up by additional text-mining, was employed in the secondary filtering of potential candidate genes to reduce bias. Fourth, any genes lacking polymorphisms (coding and non-coding) between parental strains were excluded from analyses, and lastly, positional overlap with a previously replicated major colitis susceptibility quantitative trait locus (Cdcs1) prioritised Tm3 for targeted analysis.
With regards to this shared locus, Cdcs1 on chromosome 3 was first noted in a QTL study of spontaneous colitis using IL10 deficient mice . This locus has been shown to contain at least 3 distinct regions (Cdcs1.1, 1.2 & 1.3) that contribute to a severe colitic phenotype [14, 15]. Interestingly, all three regions contribute to caecal and proximal colonic inflammation strongly suggesting that this locus is a colitis ‘hotspot’ for susceptibility and/or regulation. Here we show complete overlap of Tm3 with Cdcs1.1 (Figure 2). Although NF-kB1 has been suggested as a candidate gene for the Cdcs1 locus, it is clear that it is not responsible entirely for the severe pathology observed . To date, FcgR1 remains the key candidate gene described in the Cdcs1.1 region [14, 15] and is corroborated by our findings. Additional association with colitis susceptibility in Gnai2-/- mice  suggests that this locus may govern key inflammatory pathways in disease development, irrespective of trigger. QTL mapping specifically highlighted the Cdcs1.3 region in the spontaneous colitis and colorectal cancer development of TRUC mice . However, more distal colonic disease (distal third of the colon) or the potential for malignant transformation may not be represented at this sub-locus.
We have shown that at least 6 biologically significant and polymorphic candidate genes lie within the Cdcs1.1 autosomal region. Importantly, 4 of these candidate genes are key in pathways relevant in the context of human Crohn’s disease (FcgR1, Vav3, Vcam1 and Ctss), a disease with highly similar pathology to both the IL10 deficient and T. muris models of colonic inflammation [9, 24]. The remaining 2 genes are highly polymorphic and known to be important in inflammation (RORc, Ptpn22). As individual candidate genes, each demonstrates interesting biological functionality that could play a role in mucosal inflammation. For instance, FcgR1 codes for a high affinity IgG receptor, key to IgG2a-induced phagocytosis and antigen specific immune responses. In the mouse, FcgR1 has been associated with autoimmune disorders such as rheumatoid arthritis and bacterial infection . In humans the closely related FcgR2a and FcgR3 have been associated with IBD . The protein tyrosine phosphatase gene (Ptpn22) is of particular interest, as in humans a mis-sense SNP (C1858T) has already demonstrated strong correlation with rheumatoid arthritis , type-1 diabetes mellitus , and other autoimmune disease . Interestingly, the C1858T gene variant is not associated with the establishment of human Crohn’s disease  and may even represent protection . In our study, Ptpn22 demonstrated progressively increased expression within the colonic tissue of susceptible mice following the establishment of colitis.
The unbiased approach we have used to select candidate genes has also highlighted a gene whose currently assigned pathway (circadian rhythm) does not overtly relate to mucosal inflammation. The Retinoic acid-related orphan receptor-C (RORc/RORγ) gene encodes for RORγt (RORγ2), a lineage-specific transcription factor of CD4+ TH17 cell differentiation . Excessive TH17 cell activity has been implicated in both autoimmune  and inflammatory bowel diseases .
Finally, Vav3 was the primary candidate revealed by integrative pathway and SNP analysis and is of particular interest, as in six week old Vav1/2/3 triple knockout mice altered gut enterocyte differentiation and morphology has been shown, along with spontaneous colitis and ulceration in the caecum and ascending colon . Vav3 is also involved in at least 7 known biological pathways, all of which could play a role in mucosal homeostasis and regulation. Some of these pathways involve other candidate genes in this region, for instance FcgR1 (Fc-gamma receptor mediated phagocytosis), Vcam1 (leukocyte transendothelial migration, focal adhesion) and Ptpn22 (negative regulation of T-cell receptor signalling) [36, 37]. We hypothesise therefore that Cdcs1 is in fact a ‘colitis hotspot’ containing several genes which if dysregulated through genetic variation, could adversely affect gut inflammation. It is possible that the specific candidate genes for each colitis model are not the same. However, the biological interaction between genes at this locus, demonstrates the importance of Cdcs1 and why this region appears in unrelated models of gut inflammation.
Interestingly, Tm3 (Cdcs1) does not correlate with any known nematode infection susceptibility QTL, but instead appears exclusive to colonic inflammatory disease.
For instance, expulsion and resistance to the small intestinal nematode Heligmosomoides bakeri in mice has been characterized at murine chromosome 1 and 17  corresponding to Tm1 and Tm17. Similarly, a study of Trichinella spiralis infection in rats, which causes acute and transient small bowel inflammation, identified a single significant QTL region homologous to the murine chromosome 1 locus (Tm1) . Lastly, resistance to small bowel and abomasum/gastric nematode infections of sheep, have highlighted a number of suggestive QTL [40–42], homologous to Tm17, and downstream of Tm10. All studies demonstrated that resistance/susceptibility to GI nematode infection is under multi-genetic control, with MHC and non-MHC loci important in outcome . However, these studies also highlight the importance of the Cdcs1 locus with the establishment of a large bowel inflammatory phenotype, separate to precise anti-parasitic mechanisms.
In conclusion, we have corroborated three previously published studies which associate the locus Cdcs1 with colonic mucosal inflammation in the mouse. Furthermore, we have shown that in the AKR and BALB/c, genetic variation in this region has the potential to affect mucosal homeostasis through several different pathways. Most importantly, we have demonstrated that an unbiased integrative analysis can be beneficial in candidate gene identification and prioritization, particularly cis-regulated genes, even in large regions. This approach is particularly useful for hypothesis generation, and has positionally implicated Vav3 as a biologically relevant gene candidate in colitis.
Mice were housed with free access to food and water under specific pathogen free conditions. All experiments were performed under regulations of The UK Home Office Animals (Scientific Procedures) Act of 1986.
For QTL analysis, AKR/OlaHsd (susceptible, hereafter referred to AKR) and BALB/cOlaHsd (resistant, hereafter referred to BALB/c) mice (Harlan Olac, UK) were interbred. To generate an F1 population of mice, equal numbers of AKR males vs BALB/c females (F1a offspring), and AKR females vs BALB/c males (F1b offspring) were mated. At least 50 breeding-pairs of F1-mice were then interbred. All F1 vs F1 breeding was performed over the same time-period. To maintain genetic balance, F1a males were bred with F1a and F1b females; and, F1b males with F1a and F1b females. A single generation of 307 F2 mice (male and female) was created for study. All F2 mice were infected at the same time with T. muris ova at 6-8 weeks of age.
Trichuris muris parasites were harvested and ova collected and maintained as previously described . All infected mice received 300 T. muris ova in distilled water (200 μl) by oral gavage.
Phenotypic analysis was performed for all 307 F2 mice. Day 35 post-infection, serum samples and intestines were taken at autopsy. Resistance (0 worm load) and susceptibility (>0 worms) were defined. All worm counts were performed by a single investigator over 1 week from caeca frozen at autopsy. This method of storage and counting is used routinely for large experiments and does not affect quantification. Parasite-specific antibody ELISA was performed as described previously , using in-house T. muris excretory-secretory (ES) protein. T. muris specific IgG1 (TH2 specific, driven by IL4) and IgG2a (TH1 specific, driven by IFNγ) optical density (OD) was measured simultaneously for all samples. All 307 serum ELISAs were performed in one run. For histology, 0.5cm of whole colonic segments from the proximal ascending colon was snap-frozen, thawed in 4% paraformaldehyde, paraffin embedded, and 5μm transverse tissue sections stained with Haematoxylin and Eosin (H&E) simultaneously. 50 randomly assigned colonic specimens were assessed. Proximal colonic specimens were scored according to colonic crypt length (μm), immune cell infiltration and tissue inflammation by a single investigator. Colonic crypt length per individual was taken as a mean across at least 20 crypt units and 3 separate sections. Crypt units were measured using Image-J software . Spearman’s rank correlation coefficient was performed to measure the statistical dependence of (a) worm count and colonic crypt length variables, and (b) IgG1:IgG2a ratio and crypt length variables.
DNA was isolated (Promega Wizard DNA isolation kit) from tail snips digested in proteinase K digestion buffer (20 mg/ml). DNA concentration was determined by Nanodrop spectrophotometer and then stored at -80°C until analysis.
165 polymorphic murine microsatellite markers distinguishing between AKR and BALB/c were selected . Whole genome coverage was 85% and median inter-marker distance 12.3 cM. Conversion of marker positions from recombination fraction (cM) to physical position (Mb) was achieved using the Ensembl database .
Microsatellite amplification and genotype analysis
Forward polymerase chain reaction (PCR) primers were fluorescently labelled with 6-FAM, HEX or NED (MWG Biotec, Applied Biosystems). 25 ng of genomic DNA was used for each marker. Semi-automated analysis of genotypes on pooled panels of PCR products was performed using an Applied Biosystems 3100 Capillary sequencer with Genescan analysis and Genotyper software.
IgG1:IgG2a ratios were log10 transformed to achieve parametric distribution. Median, mean and kurtosis values were calculated using QStat, Windows QTL Cartographer 2.5. Normalised data were analysed using multiple interval mapping to optimize and refine QTL positions. A genome-wide permutation test (1000 repeats) determined thresholds for significance; a logarithm of odds (LOD) score of 4.0 or a p value of <5.2×10-5 was considered significant. A LOD score of 2.5 or p value of 1.6×10-3 was considered suggestive of linkage according to published guidelines . All significant LOD scores were confirmed by 1-way ANOVA with pairwise comparison, using the Bonferroni correction method. Kruskal-Wallis analysis was used for worm burden and IgG2a data, and converted to an LOD score .
Genome-wide colonic transcriptional activity of parental murine strains
Naïve and infected 6-to-8 week old male AKR and BALB/c mice (Harlan Olac, UK) were monitored through to day 35 post-infection (n = 6, 3 experimental replicates for each conditional cohort) as described previously . 3 replicate pooled samples of colonic RNA (ascending colon) were generated for each experimental group. Whole transcriptome microarray expression analysis (Affymetrix Genechip Mouse Exon 1.0 ST Array®) and bioinformatic analysis was performed. The entire genome-wide expression dataset was used for subsequent analysis [9, 50].
The in-silicoprioritization of QTL candidate genes
The use of workflows in the analysis of large-scale genomic data provides a systematic and un-biased mechanism for hypothesis generation . Previously constructed workflows were re-used for the analysis of QTL and gene expression data, to identify biological pathways which correlated with Trichuris muris infection. The identification of candidate genes underlying each QTL was carried out by firstly determining the precise co-ordinates of each genetic marker (Mbp) (Table 1). Each QTL was subsequently entered into the workflow qtl_to_pathway (Additional file 1: Figure S2 ). Genes located within each QTL were annotated with additional accession number identifiers (including UniProt ID and Entrez Gene IDs), in order to cross-reference Ensembl database identifiers to KEGG (Kyoto Encyclopaedia of Genes and Genomes)  pathway identifiers. As a result, annotated biological pathways were extracted from the KEGG database for inclusion in further analysis.
In parallel, differentially expressed genes identified from the T. muris microarray study  were analysed using the refseq_ids_to_pathways workflow (Additional file 1: Figure S3 ). This workflow required preliminary analysis of the gene expression data  (Partek Genomic Solution version 6.5, 2009, Partek, USA) and conversion of Affymetrix probe-set identification markers to their recognised NCBI RefSeq identification code (refseq ids). An identical process to that of the qtl_to_pathway workflow for gene annotation was then carried out.
The mapping of gene expression data to KEGG highlighted biological pathway activity in the pathogenesis of colonic disease. All genes with significant transcriptional differences between resistant and susceptible strains, in naïve and infected states (ANOVA, factor interaction, p <0.05), were included for analysis. To identify cis-QTL genes of biological relevance to phenotype, those genes with a higher degree of over/under expression (Fold Change ≥ +/-1.4 over naïve levels) during chronic T. muris intestinal inflammation, were used in the workflow analysis (see Figure 3).
The workflow common_pathways (Additional file 1: Figure S4 ) was used to identify candidate pathways containing differentially expressed genes within a QTL, in order to obtain an overall view of the mechanisms which may be influencing the expression of the phenotype.
Additional text mining was used to prevent potential candidate genes which lacked KEGG pathway annotation from being discarded. Transcribed QTL genes were analysed using a text mining workflow (Additional file 1: Figure S5 ). Briefly, published abstracts were identified from a PubMed search using the term “(“Colitis” AND “Inflammation”) AND (“Human” OR “Mouse”)”. All scientifically relevant keywords contained within individual abstracts were extracted, constructing a phenotype concept profile and allowing the calculation of inverse document frequency (IDF) scores ie a score relating the number of resulting documents which contained the keywords in question. In parallel, abstracts pertaining to selected genes were similarly recorded. The identification of phenotype keywords within individual gene abstracts allowed for the generation of a cosine vector score for each gene ranging from +1 to -1 (+1 = causation of phenotype; 0 = unknown association with phenotype; -1 = preventative of phenotype). Ranked by their cosine vector score, the association with phenotype of a particular gene was displayed. Similarly, individual phenotype keywords were also ranked according to the IDF scores, identifying possible correlations between each gene and the phenotype. All data regarding text mining and workflow approaches are published online .
Only QTL genes known to possess SNP variation between parental AKR and BALB/c  were subject to further analysis.
Independent replication of candidate gene expression by qPCR (Tm3)
Infected parental strains AKR and BALB/c (Harlan Olac, UK) received 300 T. muris ova by oral gavage. Mice were culled days 0 (naive), 7, 14, 21 and 35 post-infection for analysis (n = 3 for each cohort). mRNA was extracted from 0.5 cm of whole colonic tissue segments, from the ascending colon, according to manufacturer’s instruction (TRIZOL®, Invitrogen). cDNA was synthesised. A full list of gene primers (Eurofins-MWG-Operon, Germany) and their sequences are provided (Additional file 1: Table S1). Samples were quantitatively analysed using KAPA SYBR FAST qPCR Master Mix (Kapa Biosystems Inc., USA) and a Bio-Rad MyIQ™ PCR detection system (Bio-Rad IQ5 optical system software, version 2; Bio-Rad Laboratories Inc.,©). Three replicate cDNA samples were run at a 1:20, a 1:100, and a 1:500 dilutions for each time-point. Threshold cycles were calculated; gene detection within the three serially diluted samples was standardized, and then normalized against housekeeping gene beta-actin (Act-b). Relative fold change in gene quantity was calculated using naïve resistant mice as a reference.
Richard K Grencis and Joanne L Pennock co-senior author.
The authors acknowledge the Wellcome Trust (RKG, JLP) and MRC (SL Clinical Fellowship) for funding this work. Also the staff in the BSF of the University of Manchester for technical assistance and support.
- Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared control. Nature. 2007, 447: 661-678. 10.1038/nature05911.View Article
- McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ionnidis JP, Hirschhorn JN: Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008, 9: 356-369. 10.1038/nrg2344.View ArticlePubMed
- Glazier AM, Nadeau JH, Altman TJ: Finding genes that underlie complex traits. Science. 2002, 298: 2345-2349. 10.1126/science.1076641.View ArticlePubMed
- Stylianou IM, Affouritt JP, Shockley KR, Wilpan RY, Abdi FA, Bhardwaj S, Rollins J, Churchill GA, Paigen B: Applying gene expression, proteomics and single-nucleotide polymorphism analysis for complex trait gene identification. Genetics. 2008, 178: 1795-1805. 10.1534/genetics.107.081216.PubMed CentralView ArticlePubMed
- Olofsson P, Holmberg J, Tordsson J, Lu S, Akerstrom B, Holmdahl R: Positional identification of Ncf1 as a gene that regulates arthritis severity in rats. Nat Genet. 2003, 33: 25-32. 10.1038/ng1058.View ArticlePubMed
- Cliffe LJ, Grencis RK: The Trichuris muris system: a paradigm of resistance and susceptibility to intestinal nematode infection. Adv Parasitol. 2004, 57: 255-307.View ArticlePubMed
- Bancroft AJ, McKenzie AN, Grencis RK: A critical role for IL-13 in resistance to intestinal nematode infection. J Immunol. 1998, 160: 3463-3461.
- Else KJ, Finkelman FD, Malizewski CR, Grencis RK: Cytokine mediated regulation of chronic intestinal helminth infection. J Exp Med. 1994, 179: 347-351. 10.1084/jem.179.1.347.View ArticlePubMed
- Levison SE, McLaughlin JT, Zeef LA, Fisher P, Grencis RK, Pennock JL: Colonic transcriptional profiling in resistance and susceptibility to trichuriasis: phenotyping a chronic colitis and lessons for iatrogenic helminthosis. Inflamm Bowel Dis. 2010, 16: 2065-2079. 10.1002/ibd.21326.View ArticlePubMed
- Awasthi S, Bundy DA, Savioli L: Helminthic infections. BMJ. 2003, 327: 431-433.PubMed CentralView ArticlePubMed
- Borm ME, He J, Kelsall B, Pena AS, Strober W, Bouma G: A major quantitative trait locus on mouse chromosome 3 is involved in disease susceptibility in different colitis models. Gastroenterology. 2005, 128: 74-85. 10.1053/j.gastro.2004.10.044.View ArticlePubMed
- Farmer MA, Sundberg JP, Bristol IJ, Churchill GA, Elson CO, Leiter EH: A major quantitative locus on chromosome 3 controls colitis severity in IL-10-deficient mice. Proc Natl Acad Sci USA. 2001, 98: 13820-13825. 10.1073/pnas.241258698.PubMed CentralView ArticlePubMed
- Ermann J, Garrett WS, Kuchroo J, Rorida K, Glickman JN, BLeich A, Glimcher LH: Severity of innate immune-mediated colitis is controlled by the cytokine deficiency-induced colitis susceptibility-1 (Cdcs1) locus. Proc Natl Acad Sci U S A. 2011, 108: 7137-7141. 10.1073/pnas.1104234108.PubMed CentralView ArticlePubMed
- Bleich A, Buchler G, Beckwith J, Petell LM, Affourtit JP, King BL, Shaffer DJ, Roopenian DC, Hedrich HJ, Sundberg JP, Leiter EH: Cdcs1 a major colitis susceptibility locus in mice; subcongenic analysis reveals genetic complexity. Inflamm Bowel Dis. 2010, 16: 765-775. 10.1002/ibd.21146.PubMed CentralView ArticlePubMed
- Beckwith J, Cong Y, Sunderberg JP: Cdcs1, a major colitogenic locus in mice, regulates innate and adaptive immune responses to enteric bacterial antigens. Gastroenterology. 2005, 129: 1473-1484. 10.1053/j.gastro.2005.07.057.View ArticlePubMed
- Fisher P: Pathways and Gene annotations for QTL region. http://www.myexperment.org/workflows/1661.html.
- Fisher P: Pathways and Gene annotations for RefSeq ids. http://www.myexperiment.org/workflows/1662.html,
- Wakelin D: Acquired immunity to Trichuris muris in the albino laboratory mouse. Parasitology. 1967, 57: 515-524. 10.1017/S0031182000072395.View ArticlePubMed
- Artis D, Grencis RK: The intestinal epithelium: sensors to effectors in nematode infection. Mucosal Immunol. 2008, 1: 252-264. 10.1038/mi.2008.21.View ArticlePubMed
- Kowaiwa K, Sugawara K, Smith MF, Carl V, Yamschikov V, Belyea B, McEwen SB, Moskaluk CA, Pizarro TT, Cominelli F, McDuffie M: Identification of a quantitative trait locus for ileitis in a spontaneous mouse model of crohn’s disease: SAMP/YitFc. Gastroenterology. 2003, 125: 477-490. 10.1016/S0016-5085(03)00876-X.View Article
- Flint J, Valdar W, Shifman S, Mott R: Strategies for mapping and cloning quantitative trait genes in rodent. Nat Rev Genetics. 2005, 6: 271-286.View ArticlePubMed
- Mouse Phenome Database. http://www.informatics.jax.org,
- Verdugo RA, Farber CR, Warden CH, Medrano JF: Serious limitations of the QTL/Microarray approach for QTL gene discovery. BMC Biol. 2010, 8: 96-10.1186/1741-7007-8-96.PubMed CentralView ArticlePubMed
- Bristol IJ, Farmer MA, Cong Y, Zheng XX, Strom TB, Elson CO, Sundberg JP, Leiter EH: Heritable susceptibility for colitis in mice induced by IL-10 deficiency. Inflamm Bowel Dis. 2000, 6: 290-302.PubMed
- Iaon-Facsinay A, de Kimpe SJ, van Lent PL, Hofhuis FM, Van Ojik HH, Sedlik C, da Silveira SA, Gerver J, de Jong YF, Roozendaal R, Aarden LA, van den Berg WB, Saito T, Mosser D, Amigorena S, Izui S, van Ommen GJ, van Vugt M, van de Winkel JG, Verbeek JS: Fc-gamma-R1 (CD64) contributes substantially to severity of arthritis, hypersensitivity responses, and protection from bacterial infection. Immunity. 2002, 16: 391-402. 10.1016/S1074-7613(02)00294-7.View Article
- Weersma RK, Crusius JB, Roberts RL, Koeleman BP, Palomino-Morales R, Wolkamp S, Hollis-Moffatt JE, Festen EA, Meisneris S, Heijmans R, Noble CL, Gearry RB, Barclay ML, Gomez-Garcia M, Lopez-Nevot MA, Nieto A, Rodrigo L, Radstake TR, van Bodegraven AA, Wijmenga C, Merriman TR, Stokkers PC, Pena AS, Martin J, Alizadeh BZ: Association of FcgR2a, but not FcgR3a, with inflammatory bowel disease across three Caucasian populations. Inflamm Bowel Dis. 2010, 16: 2080-2089. 10.1002/ibd.21342.View ArticlePubMed
- Begovich AB, Carlton VE, Honigberg LA, Schrodi SJ, CHokkalingam AP, Alexander HC, Ardlie KG, Huang Q, Smith AM, Spoerke JM, Conn MT, Chang M, Chang SY, Siki RK, Catanese JJ, Leong DU, Garcia VE, McAllister LB, Jeffery DA, Lee AT, Batliwalla F, Remmers E, Criswell LA, Seldin MF, Kastner DL, Amos CI, Sninsky JJ, Gregersen PK: A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet. 2004, 75: 330-337. 10.1086/422827.PubMed CentralView ArticlePubMed
- Bottini N, Musumeci L, Alonso A, Rahmouni S, Nika K, Rostamkhani M, MacMurray J, Meloni GH, Lucarelli P, Pellecchia M, Eisenbarth GS, Comings D, Mustelin T: A functional variant in lymphoid tyrosine phosphatase is associated with type I diabetes. Nat Genet. 2004, 36: 337-338. 10.1038/ng1323.View ArticlePubMed
- Lee YH, Rho YH, Choi SJ, Ji JD, Song GG, Nath SK, Harley JB: The Ptpn22 C1858T functional polymorphism and autoimmune disease – a meta-analysis. Rheumatology. 2007, 46: 49-56. 10.1093/rheumatology/kel170.View ArticlePubMed
- van Oene M, Wintle RF, Lui X, Yazdanpanah M, Gu X, Newman B, Kwan A, Johnsn B, Owen J, Greer W, Mosher D, Maksymowych W, Keystone E, Rubin LA, Amos CI, Siminovitch KA: Association of the lymphoid tyrosine phosphatase R620W variant with rheumatoid arthritis, but not Crohn’s disease, in Canadian populations. Arthritis Rheum. 2005, 52: 1993-1998. 10.1002/art.21123.View ArticlePubMed
- Barrett JC, Hansoul S, Nicolae DL, Cho JH, Duerr R, Rioux JD, Brant SR, Silverberg MS, Raylor D, Barmada MM, Bitton A, Dasopoulos T, Datta LW, Green T, Riffiths AM, Kistner EO, Murtha MT, Reguieiro MD, Rotter JI, Schumm LP, Steinhart AH, Targan SR, Xavier RJ, Libioulle C, Sandor C, Lthrop M, Belaiche J, Dewit O, Gut I, NIHHKIBD Genetics Consortium: Genome-wide association defines more than 30 distinct susceptibility loci for Crohn’s disease. Nat Genet. 2008, 40: 955-962. 10.1038/ng.175.PubMed CentralView ArticlePubMed
- Ivanov I, McKenzie BS, Zhou L, Tadokoro CE, Lepelley A, Lafaille JJ, Cua DJ, Littman DR: The orphan nuclear receptor RORγt directs the differentiation program of proinflammatory IL-17+ T-helper cells. Cell. 2006, 126: 1121-1133. 10.1016/j.cell.2006.07.035.View ArticlePubMed
- Martinez GJ, Nurieva RI, Yang XO, Dong C: Regulation and function of proinflammatory TH17 cells. Ann N Y Acad Sci. 2008, 1143: 188-211. 10.1196/annals.1443.021.View ArticlePubMed
- Sarra M, Pallone F, MacDonald TT, Monteleone G: I23/IL17 axis in IBD. Inflamm Bowel Dis. 2010, 16: 1808-1013. 10.1002/ibd.21248.View ArticlePubMed
- Liu JY, Seno H, Miletic AV, Mills JC, Swat W, Stappenbeck TS: Vav proteins are necessary for correct differentiation of mouse caecal and colonic enterocytes. J Cell Sci. 2009, 122: 324-334. 10.1242/jcs.033720.PubMed CentralView ArticlePubMed
- Zeng L, Sachdev P, Yan L, Chan JL, Trengle T, McClelland M, Welsh J, Wang LH: Vav3 mediates receptor protein tyrosine kinase signalling, regulates GTPase activity, modulates cell morphology and induces cell transformation. Mol Cell Biol. 2000, 20: 9212-9224. 10.1128/MCB.20.24.9212-9224.2000.PubMed CentralView ArticlePubMed
- Tybulewicz VLJ: Vav-family proteins in T-cell signalling. Curr Opin Immunol. 2005, 17: 267-274. 10.1016/j.coi.2005.04.003.View ArticlePubMed
- Behnke JM, Iraqi FA, Mugambi JM, Clifford S, Nagda S, Wakelin D, Kemp SJ, Baker RL, Gibson JP: High resolution mapping of chromosomal regions controlling resistance to gastrointestinal nematode infections in an advanced intercross line of mice. Mamm Genome. 2006, 17: 584-597. 10.1007/s00335-005-0174-0.View ArticlePubMed
- Suzuki T, Ishih A, Kino H, Muregi FW, Takabayashi S, Nishikawa T, Takagi H, Terada M: Chromosomal mapping of host resistance to Trichinella spiralis nematode infection in rats. Immunogenetics. 2006, 58: 26-30. 10.1007/s00251-005-0079-9.View ArticlePubMed
- Beh KJ, Hulme DJ, Callaghan MJ, Leish Z, Lenane I, Windon RG, Maddox JF: A genome scan for quantitative trait loci affecting resistance to Trichostrongylus colubriformis in sheep. Anim Genet. 2002, 33: 97-106. 10.1046/j.1365-2052.2002.00829.x.View ArticlePubMed
- Dominik S: Quantitative trait loci for internal nematode resistance in sheep: a review. Genet Sel Evol. 2005, 37 (Suppl.1): S83-S96.PubMed CentralView ArticlePubMed
- Stear MJ, Boag B, Cattadori I, Murphy L: Genetic variation in resistance to mixed, predominantly Teladorsagia circumcincta nematode infection of sheep: from heritabilities to gene identification. Parasite Immunol. 2009, 31: 274-282. 10.1111/j.1365-3024.2009.01105.x.View ArticlePubMed
- Behnke JM, Menge DM, Noyes H: Heligmosomoides bakeri: a model for exploring the biology and genetics of resistance to chronic gastrointestinal nematode infections. Parasitology. 2009, 136: 1565-1580. 10.1017/S0031182009006003.View ArticlePubMed
- Else KJ, Entwistle GM, Grencis RK: Correlation between worm burden and markers of Th1 and Th2 cell subset induction in an inbred strain of mouse infected with Trichuris muris. Parasite Immunol. 1993, 15: 595-600.PubMed
- Image-J Software. http://rsbweb.nih.gov/ij,
- Broad Institute. http://www.broad.mit.edu/mouse,
- Ensembl Database: Ensembl Database. http://www.ensembl.org] Ensembl 64 assembly accessed September 2011
- Lander E, Kruglyak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995, 11: 241-247. 10.1038/ng1195-241.View ArticlePubMed
- Broman KW: Mapping quantitative trait loci in the case of a spike in the phenotype distribution. Genetics. 2003, 163: 1169-1175.PubMed CentralPubMed
- Genome-wide array data: Genome-wide array data. available at: [http://www.ebi.ac.uk/arrayexpress/experiments/E-MEXP-3098]
- Fisher P, Hedeler C, Wolstencroft K, Hulme H, Noyes H, Kemp S, Stevens R, Brass A: A systematic strategy for large-scale analysis of genotype–phenotype correlations: identification of candidate genes involved in African trypanosomiasis. Nucleic Acids Res. 2007, 35: 5625-5633. 10.1093/nar/gkm623.PubMed CentralView ArticlePubMed
- Kanehisa M, Goto S: KEGG: Kyoto encyclopaedia of genes and genome. Nucleic Acid Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMed
- Fisher P: KEGG pathways common to both QTL and microarray based investigations. http://www.myexperiment.org/workflows/1663,
- Fisher P: Pathway and Gene to PubMed; a text mining workflow. http://www.myexperiment.org/workflows/1846,
- Fisher P: http://www.myexperiment.org/packs/169,
- Sanger Mouse SNP Repository. http://www.sanger.ac.uk/cgi-bin/modelorgs/mousegenomes/snps/pl,
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.