Article Text

Download PDFPDF

Analysis of the influence of OCTN1/2 variants within the IBD5 locus on disease susceptibility and growth indices in early onset inflammatory bowel disease
  1. R K Russell1,
  2. H E Drummond1,
  3. E R Nimmo1,
  4. N H Anderson2,
  5. C L Noble1,
  6. D C Wilson3,
  7. P M Gillett4,
  8. P McGrogan5,
  9. K Hassan5,
  10. L T Weaver6,
  11. W M Bisset7,
  12. G Mahdi7,
  13. J Satsangi1
  1. 1Gastrointestinal Unit, Molecular Medicine Centre, University of Edinburgh, Western General Hospital, Edinburgh, UK
  2. 2Public Health Sciences, University of Edinburgh
  3. 3Department of Child life and Health, University of Edinburgh and Department of Paediatric Gastroenterology and Nutrition, Royal Hospital for Sick Children, Edinburgh
  4. 4Department of Paediatric Gastroenterology and Nutrition, Royal Hospital for Sick Children, Edinburgh
  5. 5Department of Paediatric Gastroenterology, Yorkhill Hospital, Glasgow, UK
  6. 6Department of Child Health, University of Glasgow
  7. 7Department of Paediatric Gastroenterology, Royal Aberdeen Children’s Hospital, Aberdeen, UK
  1. Correspondence to:
    Dr R K Russell
    Gastrointestinal Unit, Molecular Medicine Centre, University of Edinburgh, Western General Hospital, Crewe Road, Edinburgh EH4 2XU, UK; richardkrussell71{at}hotmail.com

Abstract

Background and aims: The OCTN1 (SLC22A4 1672C→T) and OCTN2 (SLC22A5 −207G→C) variants within the IBD5 locus have been associated with susceptibility to adult onset Crohn’s disease (CD), but their contribution in children has not been examined.

Methods: These OCTN1/2 variants and IBD5 marker single nucleotide polymorphisms (SNPs) (IGR2096a_1, IGR2198a_1, and IGR2230a_1) were examined in 299 Scottish children (200 with CD, 74 with ulcerative colitis (UC), and 25 with indeterminate colitis (IC)), together with 502 parents (for transmission disequilibrium testing) and 256 controls.

Results: All SNPs were in strong linkage disequilibrium (D′ >0.94). TDT analysis showed association of the OCTN1 variant with inflammatory bowel disease (IBD) (p = 0.01) and CD (p = 0.04). Allele frequencies of the OCTN1/2 variants were significantly higher in IBD/CD cases (p<0.04). The homozygous mutant OCTN1/2 haplotype was increased in IBD (24.3% v 16.1%, p = 0.02) and UC (28.2% v 16.1%, p = 0.02) compared with controls. The OCTN1/2 variants were not independent of the background IBD5 risk haplotype in conferring disease susceptibility. Unifactorial analysis in CD patients showed that carriage of the TC haplotype was associated with lower weight, height, and BMI centile (<9th centile) at diagnosis (weight: 87.9% v 67.3% (p = 0.002), odds ratio (OR) = 3.52 (95% confidence interval, 1.51 to 8.22); height: 84.1% v 68.4% (p<0.05), OR = 2.44 (1.00 to 5.99); BMI: 79.6% v 61.1% (p = 0.02), OR = 2.49 (1.14 to 5.44)), and lower weight centile at follow up (87.5% v 64.6% (p = 0.03), OR = 3.83 (1.03 to 14.24)). Multifactorial binary logistic regression analysis confirmed association of the TC haplotype with lower weight centile at diagnosis (p = 0.02, OR = 3.41 (1.20 to 9.66)).

Conclusions: These data implicate variants within the IBD5 haplotype, as determinants of disease susceptibility and growth indices in early onset IBD. The OCTN1/2 variants remain potential positional candidate genes, but require further analysis.

  • CD, Crohn’s disease
  • IBD, inflammatory bowel disease
  • IC, indeterminate colitis
  • SNP, single nucleotide polymorphism
  • UC, ulcerative colitis
  • Crohn’s disease
  • ulcerative colitis
  • OCTN
  • IBD5
  • growth

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

The inflammatory bowel diseases (IBD) are chronic relapsing inflammatory conditions affecting the gastrointestinal tract, comprising Crohn’s disease (CD), ulcerative colitis (UC), and indeterminate colitis (IC). The aetiology of IBD is unknown, although there is strong evidence for a gene–environment interaction.1 Around 15–25% of cases will present in childhood.2 The incidence of paediatric IBD in Scotland is higher than elsewhere in the United Kingdom and has increased threefold in the past 30 years.3–5

IBD starting in childhood can have profound effects on a child’s growth, pubertal development, and education. Growth failure is a common feature of childhood IBD, particularly in Crohn’s disease.6 Several factors have been linked to growth failure in Crohn’s disease including delay to diagnosis,7 jejunal inflammation,8 disease severity,9 age at diagnosis,10 and genotype.11

Genome-wide scanning in IBD has identified various susceptibility loci, with two that are particularly noteworthy in early onset disease—IBD1 and IBD5.12 The contribution of the NOD2/CARD15 gene to genetic susceptibility to Crohn’s disease varies between populations. In Scotland and northern Europe this effect is significantly less strong than in many other paediatric and adult IBD populations.13

The IBD5 locus (5q31–33) was first identified in genome-wide scanning of North American Crohn’s disease patients.14,15 Rioux et al stratified the genetic data from the Canadian genome-wide scan by age, and showed that the highest LOD (log of odds) score was found in patients with Crohn’s disease diagnosed under 16 years of age.14

In a detailed study of the IBD5 locus, Daly and colleagues reported strong linkage disequilibrium across the region, and derived a risk haplotype for Crohn’s disease that as represented by 11 marker single nucleotide polymorphisms (SNPs) in separate haplotype blocks that spanned the whole 250 kb interval.16 Heterozygotes for the IBD5 risk haplotype had a twofold increased risk of Crohn’s disease, and homozygotes a sixfold increase, but with no increased risk of ulcerative colitis.17 Several European studies have now replicated the association of IBD5 with susceptibility to adult onset Crohn’s disease,18–21 and additionally one study has shown an association with ulcerative colitis.18 Genotype–phenotype studies in adult onset disease have shown association with perianal Crohn’s disease20 and earlier age of disease onset.19 IBD5 epistasis has been demonstrated with the IBD6 locus22 and with NOD2/CARD15, for both Crohn’s disease19 and ulcerative colitis.18,23

Two variants within the IBD5 interval have been suggested to be independently associated with Crohn’s disease; variant alleles of the OCTN1 gene (SLC22A4 C/T, missense mutation) and OCTN2 (SLC22A5 −207 G/C, promoter mutation).24 Both of these genes have been suggested to play a role in carnitine transport but critical expression and functional data in IBD patients are still awaited. In the initial publication from Peltekova and colleagues, the resulting two allele risk (TC) haplotype was independently associated with susceptibility to Crohn’s disease when Crohn’s disease patients and controls who were homozygous wild type for marker SNP IGR2078a_1 were compared.24 Moreover, several subsequent adult studies have been unable to confirm that the OCTN1/2 effect is independent of the other potential determinants within the extended IBD5 haplotype.25,26

In this study we have analysed the contribution to disease susceptibility and phenotype of three markers on the IBD5 haplotype, together with the OCTN1/2 variants and the TC haplotype within a large homogenous paediatric IBD population. We have specifically examined whether the OCTN1/2 effect is independent of other determinants within the IBD5 locus. In addition our detailed phenotypic data have allowed us to examine the effect of these markers on growth indices.

METHODS

Patients

We recruited 299 patients with IBD diagnosed at less than 16 years of age from Scottish paediatric gastroenterology centres and from the Western General Hospital, Edinburgh. Two hundred patients had an established diagnosis of Crohn’s disease, 74 ulcerative colitis, and 25 indeterminate colitis.

Parents and controls

We also enrolled 502 parents to construct family trios for transmission disequilibrium testing (TDT) (71% of the patients had complete family trios). DNA from 256 healthy adult controls was also available for case–control analysis.26

Disease phenotype

Standardised criteria were used for IBD diagnosis.27 A patient was categorised as having “indeterminate colitis” if definite evidence of chronic inflammatory bowel disease affecting the colon only was present, but the patient remained unclassifiable as either Crohn’s disease or ulcerative colitis after considering all clinical, radiological, endoscopic, and pathological findings. The location and behaviour of Crohn’s disease was assessed at the time of diagnosis and at two yearly intervals after diagnosis, categorised using the Vienna classification.28 Disease location in Crohn’s disease was also categorised and analysed according to disease location regardless of disease involvement elsewhere in the gastrointestinal tract, as previously described.13 The use of this comprehensive complementary classification system overcomes some of the recognised limitations of the Vienna classification, notably the difficulties of categorising patients with disease involving both the upper gastrointestinal tract and other sites.29 In this classification, patients were diagnosed with Crohn’s disease of the upper gastrointestinal tract (oesophagus, stomach, and duodenum) when biopsies from any of these sites confirmed the presence of epithelioid granulomas (not merely in the presence of chronic inflammation), or when features of macroscopic disease were present,30–32 or both. Perianal disease was defined by the presence of anal fissures, perianal abscesses, fistulae, and perianal ulcers, consistent with recommendations of Lennard-Jones et al.27 To allow comparison with the results published by other groups,20,33 we also analysed data using the more limited definition of perianal disease applied by these investigators, who excluded anal fissures.

Data collection

Disease phenotype was determined using a combination of patient questionnaire, interview, and retrospective case note review, as previously described.13 Information collected included weight, height, and body mass index (BMI) values, which were plotted on a centile chart (© Child Growth Foundation, 1996) and an appropriate centile band allocated. Each patient had their height measured on a wall mounted stadiometer. Anthropometric data were also defined using z scores. A z score for a growth parameter indicates how far and in what direction that item deviates from its distribution’s mean, expressed in units of its distribution’s standard deviation. The demographics details of the patients with Crohn’s disease and ulcerative colitis in the study are shown in tables 1 and 2, respectively.

Table 1

 Demographic variables at the time of diagnosis in the Crohn’s disease patients included in the study

Table 2

 Demographic variables at the time of diagnosis in the ulcerative colitis patients included in the study

Ethical approval was given for all participating centres. All families gave informed written consent before participation in the study.

Genotyping

Genomic DNA was extracted from blood using a modified salting-out technique,34 and resuspended in 1×TE (10 mM Tris (pH 8.0), 1 mEDTA (pH 8.0)) at a final concentration of 100 ng/μl. Three marker SNPs on the IBD5 haplotype were typed IGR2096a_1, IGR2198a_1, and IGR2230a_1. The rs1050152 polymorphism of the OCTN1 gene (SCL22A4 exon 9 1672C→T) and the rs26313667 (SLC22A5 promoter, −207G→C) polymorphism of the OCTN2 gene were also typed. The primers used in the study are given in table 3.

Table 3

 Primers used to examine each of the five single nucleotide polymorphisms involved in this study

The relative positions of these markers within the haplotype blocks described by Daly et al are shown in fig 1. All genotyping was carried out using the TaqMan system (©Biosystems). The three NOD2/CARD15 mutations (Leu1007InsC, G908R, and R702W) were genotyped as previously described.35

Figure 1

 The relative position within the haplotype blocks of the IBD5 locus (as described by Daly et al16) of the five single nucleotide polymorphisms (SNPs) within the haplotype examined in this study. Also shown is the marker SNP IGR 2078a_1 examined in the study by Peltekova et al.24 The corresponding rs numbers are given below each SNP examined in this study. In Daly’s original study, 11 separate haplotype blocks were described; the numbers shown here are the ones referred to in the original publication. IGR 2222a_1 refers to the OCTN2 variant (−207G→C) and IGR 3002a_1 refers to the OCTN1 variant (1672C→T).

Statistical analysis

Each of the five SNPs was analysed individually using allele frequencies and carriage rates for association with IBD, Crohn’s disease, and ulcerative colitis. As the total number of patients with indeterminate colitis was small these data were included in IBD analysis overall and were not analysed as an individual disease group. TDT analysis was carried out using Transmit (version 2.5)36 after the Pedcheck software program was used to exclude any potential cases of non-paternity or genotyping error.37 Genotype–phenotype relations were determined by χ2 or Fisher’s exact test using Minitab v.13 (Minitab Ltd, Coventry, UK). Multifactorial analysis was done using binary logistic regression analysis. Haplotypes were calculated using Haploview, v.3.2. All SNPs in patients and controls were in Hardy–Weinberg equilibrium. We calculated z scores using a program designed by Dr C Wright (Glasgow University) and values were compared by the Mann–Whitney test. Odds ratios (OR) and 95% confidence intervals (CI) are given.

RESULTS

Strong linkage disequilibrium was seen between the three IBD5 haplotype SNPs and the two OCTN variants (fig 2). The pairwise D′ scores for all five markers were 0.94 or above.

Figure 2

 The D′ scores of the five single nucleotide polymorphisms studied across the IBD5 haplotype and the common haplotype frequencies.

Transmission disequilibrium testing analysis

These results are shown in table 4. The IGR2198a_1 allele was associated with susceptibility to IBD (p = 0.02). Homozygosity for the mutant IBD5 haplotype (patients who were homozygous mutants for all three IBD5 haplotype marker variants IGR2096a_1, IGR2198a_1, and IGR2230a_1) was associated with susceptibility to Crohn’s disease (p<0.05). The OCTN1 variant was associated with susceptibility to IBD and Crohn’s disease (p = 0.01 and p = 0.04, respectively). The homozygous OCTN1/2 mutant haplotype was associated with susceptibility to IBD (p = 0.01).

Table 4

 Transmission disequilibrium testing results for IBD5 and OCTN variants in patients with inflammatory bowel disease

Case–control analysis

Inflammatory bowel disease

The results are given in table 5. The allele frequencies were significantly higher in IBD patients for all three IBD5 SNPs and for the OCTN1/2 variants than in healthy controls (p = 0.01 to 0.04). Homozygosity rates for the IBD5 variant alleles IGR2096a_1 and IGR2198a_1, as well as homozygosity rates for the mutant TC haplotype in patients with IBD, were significantly higher than in healthy controls:

Table 5

 Allele frequencies and carriage rates for variants in the IBD5 haplotype in patients with inflammatory bowel disease compared with healthy controls

IGR2096a_1: 24.4% v 15.2% (p = 0.008), odds ratio (OR) = 1.79 (1.16 to 2.78);

IGR2198a_1: 22.8% v 15.2% (p = 0.03), OR = 1.64 (1.06 to 2.54);

mutant TC haplotype: 24.3% v 16.1% (p = 0.02), OR = 1.67 (1.08 to 2.58).

Crohn’s disease

These results are shown in table 6. The frequencies of variant alleles in Crohn’s disease patients compared with healthy controls were significantly higher for the variants:

Table 6

 Allele frequencies and carriage rates for variants in the IBD5 haplotype in Crohn’s disease patients compared with healthy controls

IGR2096a_1: 49.2% v 42.0% (p = 0.04), OR = 1.34 (1.02 to 1.75);

IGR2198a_1: 47.7% v 41.0% (p = 0.04), OR = 1.31 (1.01 to 1.71);

OCTN2: 54.8% v 47.9% (p = 0.04), OR = 1.32 (1.01 to 1.72).

Homozygous carriage of mutant IGR2096a_1 was higher in Crohn’s disease patients than in the controls, at 23.6% v 15.2% (p = 0.03), OR = 1.72 (1.06 to 2.79).

Ulcerative colitis

These results are shown in table 7. The frequency of the IGR2198a_1 variant allele was higher in the patients with ulcerative colitis than in the controls (50.7% v 41.0% (p = 0.04), OR = 1.48 (1.02 to 2.15)). Homozygosity for the variant alleles was higher in ulcerative colitis patients than controls, as follows:

Table 7

 Allele frequencies and carriage rates for variants in the IBD5 haplotype in patients with ulcerative colitis compared with healthy controls

IGR2096a_1: 27.9% v 15.2% (p = 0.01), OR = 2.16 (1.15 to 4.05);

IGR2198a_1: 26.8% v 15.2% (p = 0.04), OR = 1.92 (1.03 to 3.58).

Homozygous carriage of the mutant TC haplotype and of a combined IBD5 haplotype (homozygosity for all three IBD5 marker SNP variants) was also higher in patients with ulcerative colitis than in controls:

mutant TC haplotype: 28.2% v 16.1% (p = 0.02), OR = 2.04 (1.10 to 3.78);

combined IBD5 haplotype: 26.6% v 14.7% (p = 0.02), OR = 2.10 (1.09 to 4.05).

Association between the OCTN1/2 TC haplotype and disease susceptibility

There was no independent association between the TC haplotype and susceptibility to IBD, Crohn’s disease, or ulcerative colitis in individuals lacking the extended IBD5 risk haplotype markers (IGR2198a_1, IGR2230a_1, and IGR2096a_1) (table 8). The data shown in table 8 do not show a significant differences between IBD, Crohn’s disease, or ulcerative colitis compared with controls in the carriage of the variant alleles representing the TC haplotype on a background of having no IGR2198a_1 risk alleles. Thus for Crohn’s disease only four of 47 patients carried the TC haplotype but no IGR2198a_1 variants, compared with 10 of 82 IGR2198a_1 variant negative controls (p = 0.57). The same lack of association was also seen with the TC haplotype if either of the two other variants (IGR2230a_1 and IGR2096a_1) was analysed similarly, or if all three IBD5 SNPs were combined (data available on request).

Table 8

 Comparison of patients and controls who are homozygous wild type for the IGR2198a_1 IBD5 haplotype marker single nucleotide polymorphism but have a risk allele in the TC haplotype

No independent effect was demonstrated if the OCTN variants 1 and 2 were considered individually rather than as the TC haplotype (data available on request).

Epistasis with NOD2/CARD15 mutations

There was no evidence of epistasis between IBD/CD patients with regard to the TC haplotype and patients carrying the common NOD2/CARD15 mutations (Leu1007InsC, G908R, and R702W) in NOD2/CARD15 carriers v non-carriers for IBD and Crohn’s disease (70.3% v 76.8% (p = 0.39) and 71.0% v 76.9% (p = 0.49), respectively).

Genotype–phenotype analysis: unifactorial analysis

For clarity, the results presented in the genotype–phenotype analysis (table 9) describe associations with the TC haplotype only, with the exception of growth indices. The phenotypic associations with the TC haplotype generally reflected the genotype–phenotype results for all SNPs analysed across the IBD5 haplotype (all SNPs in close linkage disequilibrium). Supplementary information on the IBD5 alleles and OCTN variant alleles individually is available on request.

Table 9

 Genotype–phenotype of carriage rates of the TC haplotype and homozygous TC haplotype in patients with Crohn’s disease at the time of diagnosis

Growth indices

A comparison was made between lowest three (<9th) centile groupings and the remaining seven (>9th) centiles for weight, height, and BMI. Carriage of the TC haplotype was more common in patients with lower weight, height, and BMI than in those in higher centiles at diagnosis of Crohn’s disease:

lower weight: 87.9% v 67.3% (p = 0.002), OR = 3.52 (1.51 to 8.22);

lower height: 84.1% v 68.4% (p<0.05), OR = 2.44 (1.00 to 5.99);

lower BMI: 79.6% v 61.1% (p = 0.02), OR = 2.49 (1.14 to 5.44).

Analysis using median z scores comparing homozygous carriers of the TC haplotype with non-carriers confirmed a significant difference between the groups in BMI (−1.32 v −0.51, p = 0.03) and for weight (−1.16 v −0.49, p = 0.05) but not for height (p = 0.89).

At two year follow up the TC haplotype was more common in lower than in higher weight centiles, at 87.5% v 64.6% (p = 0.03), OR = 3.83 (1.03 to 14.24).

Carriage of the TC haplotype showed a significance difference for weight centile (<25th centile v >25th centile): 83.0% v 65.8% (p = 0.01), OR = 2.54 (1.23 to 5.23) at diagnosis; 88.9% v 58.5% (p =  0.002), OR = 5.68 (1.75 to 18.37) at the two year follow up.

Disease location

Homozygous carriage of the TC haplotype was associated with gastric antral Crohn’s disease, when comparison was made with patients assessed and found to have no antral disease at diagnosis (36.7% v 14.2% (p = 0.006), OR = 3.51 (1.40 to 8.83)).

Asthma, eczema, hayfever, and atopy

Homozygosity for the TC haplotype was less often seen in patients with asthma than in those without (11.5% v 26.9% (p = 0.02), OR = 0.35 (0.14 to 0.90)); was less common in patients with hayfever than in those without (7.5% v 27.0% (p = 0.009), OR = 0.22 (0.06 to 0.75)); and was less common in patients who had all three atopic diseases combined (asthma, eczema, and hayfever) than in those without (0% v 25.3% (p = 0.02)).

Disease behaviour and need for surgery

At Crohn’s disease diagnosis and at the most recent follow up no association was seen between carriage of the TC haplotype and disease behaviour. During the period of analysis 26.9% of patients underwent abdominal surgery for Crohn’s disease. Carriage of the TC haplotype was no higher in those patients needing surgery than those who did not (77.6% v 73.9% (p = 0.61)).

Age at diagnosis

The median age at diagnosis of carriers of the TC haplotype was 11.6 years (v 10.9 years for non-carriers (p = 0.12)) and for TC homozygotes it was 12.1 years (v 11.3 years in non-carriers (p = 0.06)).

Other phenotypic variables

No association was seen with the TC haplotype and family history of IBD, extraintestinal manifestations, granulomas, or blood test abnormalities (full blood count, erythrocyte sedimentation rate, C reactive protein, or albumin concentration).

Ulcerative colitis

Carriers of the TC haplotype were more common in patients in lower height centiles at diagnosis (<25th centile v >25th centile: 85.7% v 55.9% (p = 0.02), OR = 4.74 (1.17 to 19.16)) and showed a trend towards significance for weight centile (83.3% v 59.5% (p = 0.07), OR = 3.40 (0.85 to 13.57)). Genotype–phenotype analysis in this group was limited by the small number of patients available for analysis (n = 69). No relations were found between the TC haplotype and age at diagnosis, extraintestinal manifestations, IBD family history, need for surgery, and extent of disease (data not shown).

Multifactorial analysis

Using multifactorial binary logistic regression considering the weight centile as either less than or more than the 9th centile, disease location and behaviour as defined by the Vienna classification, family history of IBD, erythema nodosum, and gastric antral Crohn’s disease, a low weight centile was predicted by carriage of the TC haplotype (p = 0.02; OR = 3.41 (1.20 to 9.66)) and the presence of antral Crohn’s disease (p = 0.003; OR = 4.02 (1.58 to 10.21)).

Using a logistic regression model to compare propensity across all 11 weight centile groupings, a very strong trend towards higher carriage of the TC haplotype at lower weight centiles and lower TC haplotype carriage at higher weight centiles (p = 0.05) was apparent at Crohn’s disease diagnosis. Homozygous carriers of the extended IBD5 risk haplotype (patients who were homozygous mutants for all three IBD5 haplotype markers, IGR2096a_1, IGR2198a_1, and IGR2230a_1) were more likely to lie within lower BMI centiles and less likely to be within the highest BMI centiles in a logistic regression model (p = 0.04).

DISCUSSION

Our data provide a series of novel and important insights into the contribution of the IBD5 locus. We confirm a role in determining susceptibility to early onset Crohn’s disease and report the first evidence of association in early onset ulcerative colitis. In addition, we have presented the first consistent genotype–phenotype data suggesting that these variants influence growth in both Crohn’s disease and ulcerative colitis. We have shown that the OCTN1/2 variants do not act independently of variants within this region in determining susceptibility, emphasising the complexity of the IBD5 contribution.

Perhaps the unique strengths of our study lie in the robust phenotypic data that accompany the genetic data and in the size of the early onset population studied. A strong association between the IBD5 haplotype (and TC haplotype) with low weight and body mass index at diagnosis of Crohn’s disease is evident, with a less significant association with height centile. The relation was demonstrated on unifactorial and multifactorial analysis for weight centile (with the TC haplotype) and for homozygosity for the extended IBD5 haplotype with BMI. Although the relation was most marked in Crohn’s disease, there was a strong suggestion that it was also present in ulcerative colitis. Given the relatively small numbers of patients with ulcerative colitis, the study may have been underpowered to assess this relation fully.

Our data represent perhaps the strongest to date relating genotype to abnormal growth variables in inflammatory bowel disease. Previous paediatric studies have tried to correlate patient genotype with growth indices but overall results have been not been consistent between populations. One study has suggested variant alleles of the −238 promoter polymorphism in the TNFα gene protect patients from height retardation.10 Tomer et al showed a relation between carriage of NOD2/CARD15 variant alleles and weight below the 5th centile (but not for height; BMI was not examined) and no follow up data were presented.11 The data in Tomer’s study were not subjected to multifactorial analysis; the findings have been confounded by the strong relation between NOD2/CARD15 carriers with ileal disease, which in turn is associated with lower weight at diagnosis.38 Two subsequent publications have failed to replicate these findings.9,39

In our own population, we have studied the relation of NOD2/CARD15 genotype with growth indices; no relation was found at diagnosis, while at follow up carriers of NOD2/CARD15 variant alleles were more common in a low centile band for height (2–9th) and for weight (<0.4th) only.13

There is no universally accepted definition of growth failure or good evidence of how it should be treated.40 Growth failure nevertheless remains a common problem, especially in paediatric Crohn’s disease, and can persist into adult life.41 There is also evidence from prospective paediatric studies that the growth failure is “preprogrammed”, strongly suggesting a genetic influence.42 Although our data do not establish that the genotype predicts “growth failure”—because our strongest growth associations were with weight and BMI, not height—they do show that genetic markers have the potential to be used to predict growth problems, and even perhaps determine treatment strategies.

The association of markers in the IBD5 haplotype with low anthropometric centiles signifying a more severe disease phenotype are entirely consistent with our recent findings in Scottish adults with Crohn’s disease.26 Possession of the variant alleles of the TC haplotype in 374 Scottish adults with Crohn’s disease was associated with disease progression at follow up and, on multifactorial analysis, the need for surgery—both phenotypic variables representing a severe disease phenotype. As discussed in an accompanying editorial to our paper,43 other adult datasets may also implicate the IBD5 locus as a determinant of severe disease.

We describe a novel association of the TC haplotype with gastric antral Crohn’s disease. Previous adult studies have linked possession of the TC haplotype with ileal,44 colonic,25 and perianal33 Crohn’s disease, but none of these studies considered a relation with Crohn’s disease of the proximal gastrointestinal tract. The establishment of such a relation gives some justification for using a phenotypic classification of Crohn’s disease location in genotype–phenotype studies that is more extensive than the Vienna classification—an issue that has been intensively examined by the 2005 Montreal World congress Working Party into the classification of IBD.29 The relation of antral disease with low weight centile at diagnosis on multifactorial analysis is intriguing in the light of previous prospective paediatric studies establishing a link between jejunal Crohn’s disease and low weight centile at diagnosis.8

Our data show an incidence of proximal gastrointestinal Crohn’s disease that may appear high in comparison with published adult series,26 but the figures are similar to the prospective paediatric Crohn’s disease studies from Europe and North America.8,45 The difference may in part reflect a different disease phenotype between paediatric and adult Crohn’s disease, but may also reflect the different investigation practices of paediatric gastroenterologists who would regard routine upper gastrointestinal endoscopy at diagnosis as standard practice.30

The association between the IBD5 haplotype and susceptibility to Crohn’s disease has now been clearly demonstrated in adult studies18–21,25,26 but here we present further data from an exclusive early onset population to replicate the findings in Rioux’s genome-wide scan.14 Our data represent the first paediatric study in ulcerative colitis to date and complement adult data.18 It is noteworthy that the index study from Canada may have been underpowered to assess the association between IBD5 and ulcerative colitis, as it contained only 20 families with ulcerative colitis.14 Our own dataset is considerably larger than the index study; nevertheless is relatively small for a case–control association study, and therefore the data have to be interpreted carefully because of the small number of patients compared with controls. The TC haplotype has been associated with Crohn’s disease susceptibility in three adult studies,24–26 but no previous association has been shown with ulcerative colitis. This difference may be explained by a greater influence of the TC haplotype in early onset disease, or may indeed reflect differing disease phenotypes in the paediatric and adult ulcerative colitis populations, notably the higher proportion of paediatric patients with more extensive colitis.8,46

Despite demonstrating an association between the OCTN1/2 variants and disease susceptibility, our data do not support these variants as independent predictors of disease when the IBD5 risk markers are considered. This finding is contrary to the findings of Peltekova’s study,24 but consistent with subsequent publications and presentations.25,26,47,48 These apparent inconsistencies are likely to reflect the haplotype block structure of the region as described by Daly et al (fig 1).16 The published studies to date, apart from our own studies in the Scottish IBD population, have only examined one haplotype marker, IGR2078a_1 (haplotype block 4), and explored its relation (or lack of relation) with the OCTN1/2 genes (haplotype block 7).24,25 Our study examined markers in block 4 (IGR2096a_1), block 5 (IGR2198a_1), and block 7 (IGR2230a_1), representing an extended haplotype across the region. Thus the apparently independent effect reported in Peltekova’s data may reflect recombination events between haplotype block 4 and 7 in the Canadian population. This issue would be best resolved by further studies in the Canadian dataset.

Given this uncertainty about the identity of the causal gene within the IBD5 locus it remains pertinent to examine the role of other genes within the region: interleukin 4,49–51 interleukin 13,52 and interleukins 3 and 5, together with GMCSF, are all plausible candidate genes within the linkage interval.53,54 The protective association of the homozygous TC haplotype with asthma, hayfever, and atopy would be consistent with the evidence from genome-wide scans that genes within the 5q31–33 region are not only important in IBD but also in asthma and atopy.55–58 Indeed, candidate genes within the 5q31–33 region, including IL12B gene59 and IL13,60 have been associated with asthma and atopy.

Conclusions

In summary, we have demonstrated evidence for genetic association between the OCTN1/2 variants within the IBD5 locus and susceptibility to inflammatory bowel disease, Crohn’s disease, and ulcerative colitis. Of note, our study also describes novel and consistent phenotypic associations with low anthropometric centiles. However, the effect of the OCTN1/2 variants is not independent of the extended IBD5 region and further investigation of these and other candidate genes within the IBD5 linkage interval is necessary.

Acknowledgments

We acknowledge the help of all patients and parents who participated in the study, together with the specialist nurses, dieticians, and secretaries in each of the teaching hospitals, as well as the paediatricians, practice nurses, and GPs throughout Scotland, whose support for the study was invaluable. We would also like to thank John Rioux and Mark Daly for their advice.

The University of Edinburgh Medical Faculty Fellowship funds RKR. This study is supported by a Wellcome Trust programme grant (072789/Z/03/Z) with additional support from Schering Plough and the GI/Nutrition research fund, Child Life and Health, University of Edinburgh.

These data were presented as an oral presentation at the 36th Digestive Diseases Week in Chicago, USA, May 2005.

REFERENCES

Footnotes

  • Published Online First 9 February 2006

  • Conflict of interest: None declared.