<title>Background</title> <p>Identification of the mechanisms of autoimmune inflammatory disease development in the gastrointestinal tract is an emerging topic of research. Inflammatory bowel disease (IBD) is thought to result from immune-mediated tissue injury, primed by the enteric microflora, in genetically predisposed subjects. IBD phenotypically occurs in the form of Crohn's disease (CD) or ulcerative colitis (UC). In these disorders, chronic mucosal inflammation results from inappropriate and overreactive mucosal response to intestinal bacteria followed by activation of inflammatory cells and production of inflammatory mediators [##REF##12167685##1##]. Primary sclerosing cholangitis (PSC) is a chronic inflammatory disease of the hepatic bile ducts that likely develops as a result of an inappropriate immune mediated process [##REF##17100974##2##]. Nearly 80% of northern European PSC patients also have concomitant IBD, mostly in the form of UC [##UREF##0##3##]. Primary biliary cirrhosis (PBC) is another autoimmune chronic cholestatic liver disorder that may also develop as a result of an abnormal immune response to stimulating environmental or infectious agents [##REF##16126958##4##].</p> <p>As in other autoimmune diseases, the specific alleles of the human leukocyte antigen (HLA) complex represent genes associated with susceptibility to PSC [##REF##11584356##5##, ####REF##12235090##6##, ##REF##17907284##7####17907284##7##], PBC [##REF##16034472##8##,##REF##15884119##9##] and IBD [##REF##17628615##10##]. It has been reported that the genetic predisposition to CD is caused by alterations in several genes, including <italic>NOD2/CARD15, OCTN1</italic>, <italic>OCTN2</italic>, <italic>DLG</italic>, <italic>ATG16L1 </italic>and <italic>IL23R </italic>[##REF##11385577##11##, ####REF##17068223##12##, ##REF##17786191##13##, ##REF##17206667##14####17206667##14##]. Several other non-HLA genes are also involved in PCS and PBC disease susceptibility [##REF##17907284##7##,##REF##17074031##15##].</p> <p>According to our current understanding of CD, PSC and PBC are complex diseases in which the genetic determinants contributing to disease susceptibility interact with environmental factors. However, the mechanisms underlying these interactions remain unclear. Considering high frequency at which IBD and PSC are diagnosed together [##UREF##0##3##], close spatial proximity of affected organs, and the fact that common immune mediated processes are likely involved in development of these diseases [##REF##15557349##16##], searching for common risk alleles for IBD and autoimmune cholestatic liver disorders seems to be relevant</p> <p>Most of the susceptibility genes are considered to not be directly involved in disease pathogenesis, but rather act as environmental response modifiers. Given their subtle individual impact on disease susceptibility, independently they may have little to no effect on disease progression and likely act in collaboration with other factors to promote disease development. Moreover, disease genetic markers vary among the populations studied, despite the phenotypic similarity of disease. Thus, analyzing the genetic background of disease is an important challenge, and several studies have been performed with the purpose of finding similarities in the molecular mechanism of development of different inflammatory disorders of the alimentary tract. The aim of this study was to investigate whether polymorphisms of CD susceptibility genes are predisposing factors to the development of PSC and PBC in a group of Polish patients.</p>
<title>Methods</title> <title>Patients</title> <p>The study included 60 patients with CD (28 men, range of 22 – 78 years old; median 39), 77 patients with PSC (45 men, range of 20 – 68 years old; median 33), of which 61 exhibited IBD as confirmed by endoscopy and histology (40 UC, 8 CD, and 13 indeterminate colitis), and 144 patients with PBC (8 men, range of 40 – 72 years old; median 58). All patients were diagnosed at the Department of Gastroenterology and Hepatology, Medical Center for Postgraduate Education and Cancer Center, Warsaw, Poland. Diagnosis of CD was based on generally accepted clinical, endoscopical and histological criteria. The diagnosis of PSC was based on characteristic biochemical and radiological features (irregularity of the intrahepatic and extrahepatic bile ducts in endoscopic retrograde or magnetic resonance cholangiopancreatography). PBC was diagnosed based on generally accepted clinical, biochemical, serological (antimitochondrial antibody positivity by indirect immunofluorescence) and histological findings [##REF##8900092##17##].</p> <p>The control group consisted of 139 healthy volunteers (74 men, range of 21 – 66 years old, median 32) recruited from hospital staff, medical students, and healthy subjects selected for screening colonoscopy. All patients and controls were Polish Caucasians. The study was approved by the local ethics committee (Medical Center for Postgraduate Education and Cancer Center, Warsaw, Poland), and all the participants provided appropriate consent. The study protocol conforms to the ethical guidelines of the 1975 Declaration of Helsinki.</p> <title>Genotyping</title> <p>Genomic DNA was extracted from whole blood treated with EDTA using the NucleoSpin Blood Quick Pure kit (Macherey-Nagel, Germany), following the manufacturer's protocol.</p> <p>The selected polymorphisms of genes with previously confirmed evidence for association with CD were as follows: <italic>NOD2/CARD15</italic>: rs2066842 (Pro268Ser), rs2066844 (Arg702Trp), rs2066845 (Gly908Arg), rs5743293 (1007fs); in the <italic>SLC22A4/OCTN1</italic>/<italic>SLC22A5/OCTN2 </italic>genes: rs1050152 (Leu503Phe)/rs2631367 (-207G&gt;C); in the <italic>DLG5 </italic>gene: rs1248696 (Arg30Gln); in the <italic>ATG16L1 </italic>gene: rs2241880 (Thr300Ala); and in the <italic>IL23R </italic>gene: rs11209026 (Arg381Gln), rs1884444 (His3Gln), rs10889677 (exon-3'UTR).</p> <p>Genotyping was carried out using TaqMan SNP Genotyping Assays (Applied Biosystems; Foster City, CA). The PCR reactions were performed in 96-well plates on the ABI PRISM 7000 Sequence Detection System (SDS) (Applied Biosystems; Foster City, CA). Each well contained 20 ng genomic DNA, 1.25 μL TaqMan SNP Genotyping Assay (probes and primers mix), 12.5 μL TaqMan Universal PCR Master Mix, No AmpErase UNG (Applied Biosystems; Foster City, CA), and 9.25 μL water. Two non-template-control wells were included on each plate. After DNA amplification (95°C for 10 min, followed by 40 cycles of 92°C for 15 sec and 60°C for 1 min), fluorescence was acquired and analyzed for allelic discrimination using the ABI Prism 7000 SDS Software.</p> <title>Statistical analyses</title> <p>The frequency distribution of alleles, genotypes, and haplotypes was compared using standard or Yates corrected χ<sup>2 </sup>test and Fisher's exact test when appropriate. Statistical significance threshold level was Bonferroni corrected for multiple hypothesis testing, according to the number of genes tested. The p-value significance threshold level was defined as 0.008. For the Fisher's exact tests, the STATISTICA software was used. ORs with 95% confidence intervals (95% CI) were calculated using a calculator available at <ext-link ext-link-type="uri" xlink:href="http://www.hutchon.net/ConfidOR.htm"/>. LD analysis and calculation of the Hardy-Weinberg equilibrium were performed using the Haploview v3.2 software, available at <ext-link ext-link-type="uri" xlink:href="http://www.broad.mit.edu/mpg/haploview/"/>. Statistical power analyses were done using G*power's (ver. 3.0.10) [##REF##17695343##18##] post-hoc procedure for the Fisher's exact test.</p>
<title>Results and Discussion</title> <p>Allelic distributions for all polymorphisms studied in the patient and control groups are shown in Table ##TAB##0##1##, while genotype counts for each of the studied groups are presented in Additional file ##SUPPL##0##1## [see Additional file ##SUPPL##0##1##]. In the control group, the genotype distributions for all polymorphisms were in Hardy-Weinberg equilibrium (p &gt; 0.05).</p> <p>The <italic>Caspase Recruitment Activation Domain 15 </italic>(<italic>NOD2/CARD15</italic>) gene is located in the proximal region of chromosome 16 (16q12) and encodes the nucleotide-binding oligomerization domain protein 2 (NOD2), which is a cytosolic receptor for a bacterial peptidoglycan response pathway [##REF##17379562##19##]. Mutations in NOD2 result in immune system dysfunctions [##REF##11385577##11##,##REF##11385576##20##] and are considered causative genetic factors in the development of CD. <italic>NOD2/CARD15 </italic>mutations are also associated with susceptibility to other granulomatous inflammatory disorders, such as early-onset sarcoidosis and Blau syndrome [##REF##11528384##21##,##REF##15459013##22##].</p> <p>Three of the four <italic>NOD2/CARD15 </italic>variants tested (rs2066842, rs2066844, rs5743293) were strongly associated with CD (Table ##TAB##1##2##). In particular, a strong association was observed with Pro268Ser and 1007fs SNPs with odds ratio (OR) values over 10. Increased values of OR for homozygous risk genotypes confirm previously described codominance of these variants. The observed allele distribution was similar to those reported in other Caucasian CD patient groups [##UREF##1##23##, ####REF##12055616##24##, ##REF##12105838##25####12105838##25##].</p> <p><italic>The autophagy-related 16-like 1 </italic>(<italic>ATG16L1</italic>) gene encodes the protein component of the autophagosome pathway of intracellular bacteria processing. The T300A variant (rs2241880) was previously considered to be a risk factor for CD development [##REF##17200669##26##,##REF##17435756##27##]. In our CD patients, however, this variant demonstrated only a weak associaf2tion with CD, as the p-value (= 0.022) did not pass the Bonferroni corrected significance threshold level of 0.008.</p> <p>The organic cation transporters from the family of solute carrier protein (<italic>SLC22A4/OCTN1 </italic>and <italic>SLC22A5/OCTN2</italic>) genes are located on chromosome 5 (5q31) within the IBD5 region [##REF##16344054##28##]. The 1672C/T (Leu503Phe) missense variant of <italic>SLC22A4/OCTN1 </italic>(rs1050152) potentially alters the carnitine adsorption and the exchange of other positively charged compounds between the cell and extracellular matrix [##REF##16835882##29##]. Also included in this study is the G-to-C transversion in the <italic>SLC22A5/OCTN2 </italic>promoter (rs2631367), which is thought to affect a heat shock transcription factor binding element [##REF##18005709##30##]. The variants of both genes were previously found in linkage disequilibrium, enabling the selection of a two-allele CD risk haplotype [##REF##16344054##28##,##REF##15107849##31##]. This finding, however, has not been confirmed by other studies [##REF##16835882##29##].</p> <p>As shown in Table ##TAB##2##3##, there was a weak association between the protective <italic>OCTN1/OCTN2 </italic>CC haplotype and CD (OR = 0.28; CI = 0.08 – 0.94; p = 0.0298). A strong linkage disequilibrium was found between the <italic>NOD2/CARD15 </italic>Pro268Ser loci and three other studied loci of <italic>NOD2/CARD15</italic>. Analyses of reconstructed haplotypes for each patient group showed that the odds of TC haplotype of Pro268Ser and 1007fs loci is approximately ten times higher (OR = 10.26; CI = 4.52 – 23.29; p = 0.0000000000671) among patients with CD compared to healthy controls.</p> <p>The <italic>interleukine-23 receptor </italic>(<italic>IL23R</italic>) gene is located on chromosome region 1p31, and encodes a subunit of the IL23 receptor. IL23 and IL6 activate the STAT3 transcription factor which in turn leads to differentiation of CD4+ helper T cells into Th17 cells that function in driving autoimmune inflammation by producing pro-inflammatory cytokines, such as IL-17A, IL-17F and IL-22 [##REF##17786191##13##,##REF##16670770##32##, ####REF##17030949##33##, ##REF##17277312##34##, ##REF##18039580##35####18039580##35##]. Genetic variants of <italic>IL23R </italic>were identified as disease-associated factors in chronic inflammatory disorders, such as IBD [##REF##17068223##12##], psoriasis, T-cell-mediated inflammatory disease of the skin [##REF##17236132##36##], and ankylosing spondylitis [##REF##17952073##37##], indicating the involvement of the IL23/IL23R pathway in the etiology of various autoimmune-related disorders. Several genetic non-coding variants of this gene demonstrated independent risk association with CD, while the rare variant in the coding region of <italic>IL23R </italic>(1142G/A; Arg381Gln; rs11209026) has been regarded as a CD protective variant [##REF##17068223##12##,##REF##17786191##13##].</p> <p>We did not find any association of CD not only with Arg381Gln, but also two other more frequent variants of <italic>IL23R</italic>. The frequency of the G/A genotype of rs11209026 SNP was 3.3% in CD (2 of 60) compared to 5.8% (8 of 139) in controls, and only one of the patients and controls was homozygous for the AA genotype (not significant). Thus, although this rare variant was found at a similar proportion to that described from the genome-wide scan-based population studies (2.5–3.0% <italic>versus </italic>6.2–6.8%; [##REF##17786191##13##,##REF##17484863##38##]), our small group of CD patients (N = 60) was insufficient to reach statistical significance for this disproportion.</p> <p>The <italic>Drosophila </italic>disc large homolog 5 (<italic>DLG5</italic>) gene, mapped to chromosome 10 (10q23), is a member of the membrane-associated guanylate kinase gene family [##REF##17628615##10##,##REF##16773680##39##] and encodes a protein involved in maintaining the correct shape, polarity and growth of cells. The functional variant (113G/A, Arg30Gln, rs1248696) has been proposed to potentially lead to serious abnormalities in both the structure and function of cells. Although the first study of Stoll <italic>et al. </italic>[##REF##15107852##40##] reported an association between genetic variants in <italic>DLG5 </italic>and IBD, the more recent meta-analysis of published data on Arg30Gln question this association [##REF##17628615##10##]. Our study did not show any association between the Arg30Gln variant of <italic>DLG5 </italic>and CD.</p> <p>In contrast to CD patients, none of the polymorphisms showed association with PSC and PBC, two disorders of autoimmune etiology. Furthermore, no significant associations were detected in PSC patients with concurrent IBD. These negative results are consistent with the recently published studies in Scandinavian PSC patients, the largest PSC population in which genotyping has been performed [##REF##17100974##2##], and in a group of Hungarian and Polish patients with PBC [##UREF##2##41##].</p> <p>The present study confirms a strong association of <italic>NOD2/CARD15 </italic>gene variants, and to a lesser extent the coding <italic>ATG16L1 </italic>variant and the <italic>OCTN1/OCTN2 </italic>haplotype, with CD in a relatively small Polish Caucasian patient group. There was no evidence for significant association between variants of <italic>IL23R </italic>or <italic>DLG5 </italic>and CD. However, because small sample size significantly limits estimation of genetic variant disproportion, especially with a low prevalence of allelic frequency, whether the result represents a true negative association or a false negative finding is somewhat disputable, mainly due to insufficient power of statistical tests [##REF##17907284##7##]. In fact, calculation of the statistical testing power reached decent values only for the three main <italic>NOD2/CARD15 </italic>variants (Table ##TAB##3##4##). Therefore, considering SNP frequencies, small patient groups and measured effect size, our studies might not be able to confirm existing weak and even moderate disease associations.</p>
<title>Results and Discussion</title> <p>Allelic distributions for all polymorphisms studied in the patient and control groups are shown in Table ##TAB##0##1##, while genotype counts for each of the studied groups are presented in Additional file ##SUPPL##0##1## [see Additional file ##SUPPL##0##1##]. In the control group, the genotype distributions for all polymorphisms were in Hardy-Weinberg equilibrium (p &gt; 0.05).</p> <p>The <italic>Caspase Recruitment Activation Domain 15 </italic>(<italic>NOD2/CARD15</italic>) gene is located in the proximal region of chromosome 16 (16q12) and encodes the nucleotide-binding oligomerization domain protein 2 (NOD2), which is a cytosolic receptor for a bacterial peptidoglycan response pathway [##REF##17379562##19##]. Mutations in NOD2 result in immune system dysfunctions [##REF##11385577##11##,##REF##11385576##20##] and are considered causative genetic factors in the development of CD. <italic>NOD2/CARD15 </italic>mutations are also associated with susceptibility to other granulomatous inflammatory disorders, such as early-onset sarcoidosis and Blau syndrome [##REF##11528384##21##,##REF##15459013##22##].</p> <p>Three of the four <italic>NOD2/CARD15 </italic>variants tested (rs2066842, rs2066844, rs5743293) were strongly associated with CD (Table ##TAB##1##2##). In particular, a strong association was observed with Pro268Ser and 1007fs SNPs with odds ratio (OR) values over 10. Increased values of OR for homozygous risk genotypes confirm previously described codominance of these variants. The observed allele distribution was similar to those reported in other Caucasian CD patient groups [##UREF##1##23##, ####REF##12055616##24##, ##REF##12105838##25####12105838##25##].</p> <p><italic>The autophagy-related 16-like 1 </italic>(<italic>ATG16L1</italic>) gene encodes the protein component of the autophagosome pathway of intracellular bacteria processing. The T300A variant (rs2241880) was previously considered to be a risk factor for CD development [##REF##17200669##26##,##REF##17435756##27##]. In our CD patients, however, this variant demonstrated only a weak associaf2tion with CD, as the p-value (= 0.022) did not pass the Bonferroni corrected significance threshold level of 0.008.</p> <p>The organic cation transporters from the family of solute carrier protein (<italic>SLC22A4/OCTN1 </italic>and <italic>SLC22A5/OCTN2</italic>) genes are located on chromosome 5 (5q31) within the IBD5 region [##REF##16344054##28##]. The 1672C/T (Leu503Phe) missense variant of <italic>SLC22A4/OCTN1 </italic>(rs1050152) potentially alters the carnitine adsorption and the exchange of other positively charged compounds between the cell and extracellular matrix [##REF##16835882##29##]. Also included in this study is the G-to-C transversion in the <italic>SLC22A5/OCTN2 </italic>promoter (rs2631367), which is thought to affect a heat shock transcription factor binding element [##REF##18005709##30##]. The variants of both genes were previously found in linkage disequilibrium, enabling the selection of a two-allele CD risk haplotype [##REF##16344054##28##,##REF##15107849##31##]. This finding, however, has not been confirmed by other studies [##REF##16835882##29##].</p> <p>As shown in Table ##TAB##2##3##, there was a weak association between the protective <italic>OCTN1/OCTN2 </italic>CC haplotype and CD (OR = 0.28; CI = 0.08 – 0.94; p = 0.0298). A strong linkage disequilibrium was found between the <italic>NOD2/CARD15 </italic>Pro268Ser loci and three other studied loci of <italic>NOD2/CARD15</italic>. Analyses of reconstructed haplotypes for each patient group showed that the odds of TC haplotype of Pro268Ser and 1007fs loci is approximately ten times higher (OR = 10.26; CI = 4.52 – 23.29; p = 0.0000000000671) among patients with CD compared to healthy controls.</p> <p>The <italic>interleukine-23 receptor </italic>(<italic>IL23R</italic>) gene is located on chromosome region 1p31, and encodes a subunit of the IL23 receptor. IL23 and IL6 activate the STAT3 transcription factor which in turn leads to differentiation of CD4+ helper T cells into Th17 cells that function in driving autoimmune inflammation by producing pro-inflammatory cytokines, such as IL-17A, IL-17F and IL-22 [##REF##17786191##13##,##REF##16670770##32##, ####REF##17030949##33##, ##REF##17277312##34##, ##REF##18039580##35####18039580##35##]. Genetic variants of <italic>IL23R </italic>were identified as disease-associated factors in chronic inflammatory disorders, such as IBD [##REF##17068223##12##], psoriasis, T-cell-mediated inflammatory disease of the skin [##REF##17236132##36##], and ankylosing spondylitis [##REF##17952073##37##], indicating the involvement of the IL23/IL23R pathway in the etiology of various autoimmune-related disorders. Several genetic non-coding variants of this gene demonstrated independent risk association with CD, while the rare variant in the coding region of <italic>IL23R </italic>(1142G/A; Arg381Gln; rs11209026) has been regarded as a CD protective variant [##REF##17068223##12##,##REF##17786191##13##].</p> <p>We did not find any association of CD not only with Arg381Gln, but also two other more frequent variants of <italic>IL23R</italic>. The frequency of the G/A genotype of rs11209026 SNP was 3.3% in CD (2 of 60) compared to 5.8% (8 of 139) in controls, and only one of the patients and controls was homozygous for the AA genotype (not significant). Thus, although this rare variant was found at a similar proportion to that described from the genome-wide scan-based population studies (2.5–3.0% <italic>versus </italic>6.2–6.8%; [##REF##17786191##13##,##REF##17484863##38##]), our small group of CD patients (N = 60) was insufficient to reach statistical significance for this disproportion.</p> <p>The <italic>Drosophila </italic>disc large homolog 5 (<italic>DLG5</italic>) gene, mapped to chromosome 10 (10q23), is a member of the membrane-associated guanylate kinase gene family [##REF##17628615##10##,##REF##16773680##39##] and encodes a protein involved in maintaining the correct shape, polarity and growth of cells. The functional variant (113G/A, Arg30Gln, rs1248696) has been proposed to potentially lead to serious abnormalities in both the structure and function of cells. Although the first study of Stoll <italic>et al. </italic>[##REF##15107852##40##] reported an association between genetic variants in <italic>DLG5 </italic>and IBD, the more recent meta-analysis of published data on Arg30Gln question this association [##REF##17628615##10##]. Our study did not show any association between the Arg30Gln variant of <italic>DLG5 </italic>and CD.</p> <p>In contrast to CD patients, none of the polymorphisms showed association with PSC and PBC, two disorders of autoimmune etiology. Furthermore, no significant associations were detected in PSC patients with concurrent IBD. These negative results are consistent with the recently published studies in Scandinavian PSC patients, the largest PSC population in which genotyping has been performed [##REF##17100974##2##], and in a group of Hungarian and Polish patients with PBC [##UREF##2##41##].</p> <p>The present study confirms a strong association of <italic>NOD2/CARD15 </italic>gene variants, and to a lesser extent the coding <italic>ATG16L1 </italic>variant and the <italic>OCTN1/OCTN2 </italic>haplotype, with CD in a relatively small Polish Caucasian patient group. There was no evidence for significant association between variants of <italic>IL23R </italic>or <italic>DLG5 </italic>and CD. However, because small sample size significantly limits estimation of genetic variant disproportion, especially with a low prevalence of allelic frequency, whether the result represents a true negative association or a false negative finding is somewhat disputable, mainly due to insufficient power of statistical tests [##REF##17907284##7##]. In fact, calculation of the statistical testing power reached decent values only for the three main <italic>NOD2/CARD15 </italic>variants (Table ##TAB##3##4##). Therefore, considering SNP frequencies, small patient groups and measured effect size, our studies might not be able to confirm existing weak and even moderate disease associations.</p>
<title>Conclusion</title> <p>Although inappropriate immune mediated processes are also related with PSC and PBC, this study and previously published [##REF##17100974##2##,##UREF##2##41##] studies have not identified common allelic variants. In IBD, risk variants of genes encoding proteins that are involved in signalling pathways activated by bacterial products have been described. Consequently, IBD is thought to result from dysfunction of the mucosal immune response triggered by intestinal bacteria, and defects in the early immune response may play a role in the pathogenesis of CD [##REF##12167685##1##,##REF##16819502##42##,##REF##16920636##43##]. However, this hypothesis is based on statistical considerations rather than functional studies. It is noteworthy that only some CD patients carry mutated <italic>NOD2/CARD15 </italic>and/or other risk gene variants, and therefore, mechanisms for CD development likely are not simply related to mutations within described risk genes. Furthermore, the genetic background of CD possibly is based on more causative genetic factors than studied thus far, and each of them independently may have a relatively weak impact on disease development.</p> <p>Such high complexity of disease background requires mapping of genetic predisposition using genome-wide and large-scale studies in well-defined populations. Despite the fact that these kind of studies are still relatively expensive, they are likely the most cost effective giving real capacity to provide comprehensive insight into disease development mechanisms.</p> <p>Even though familial aggregation is observed in many complex diseases, including IBD, PSC and PBC, genetic factors in complex diseases should be considered as "risk-factor genes" rather than genes responsible for the development of a particular disease. Their presence is not synonymous with the development of the disease. Thus, functional studies should be initiated in order to clarify the contribution of genetic background in development of these diseases.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Numerous papers have addressed the association of mutations and polymorphisms of susceptibility genes with autoimmune inflammatory disorders. We investigated whether polymorphisms that confer susceptibility to Crohn's disease could be classified also as predisposing factors for the development of primary sclerosing cholangitis and primary biliary cirrhosis in Polish patients.</p> <title>Methods</title> <p>The study included 60 patients with CD, 77 patients with PSC, of which 61 exhibited IBD (40 UC, 8 CD, and 13 indeterminate colitis), and 144 patients with PBC. All the patients were screened against Crohn's disease associating genetic polymorphisms.</p> <p>The polymorphisms were chosen according to previously confirmed evidence for association with Crohn's disease, including Pro268Ser, Arg702Trp, Gly908Arg and 1007fs in <italic>NOD2/CARD15</italic>, Leu503Phe/-207G&gt;C in <italic>SLC22A4/OCTN1</italic>/<italic>SLC22A5/OCTN2</italic>, Arg30Gln in <italic>DLG5</italic>, Thr300Ala in <italic>ATG16L1</italic>, and Arg381Gln, His3Gln and exon-3'UTR in <italic>IL23R</italic>. Genotyping was carried out using TaqMan SNP genotyping assays.</p> <title>Results</title> <p>We confirmed a strong association between three <italic>NOD2/CARD15 </italic>gene variants (Pro268Ser, OR = 2.52, 95% CI = 1.34 – 4.75); (Arg702Trp, OR = 6.65, 95% CI = 1.99 – 22.17); (1007fs, OR = 9.59, 95% CI = 3.94 – 23.29), and a weak association between both the protective <italic>OCTN1/OCTN2 </italic>CC haplotype (OR = 0.28, 95% CI = 0.08 – 0.94), and a variant of <italic>ATG16L1 </italic>gene (Thr300Ala, OR = 0.468, 95% CI = 0.24 – 0.90) with Crohn's disease. In contrast, none of the polymorphisms exhibited association with susceptibility to primary sclerosing cholangitis and primary biliary cirrhosis, including a group of primary sclerosing cholangitis patients with concurrent IBD.</p> <title>Conclusion</title> <p>Although the clinical data indicate non-random co-occurrence of inflammatory bowel disease and primary sclerosing cholangitis, consistently with the previously published studies, no genetic association was found between the genetic variants predisposing to Crohn's disease and hepatobiliary autoimmune disorders. However, since estimation of genetic variant disproportion is limited by sample size, these negative results may also indicate that eventually shared genetic predispositions are too little to be captured by small patient groups.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>PG participated in preparation of the study design, carried out the molecular assays as well as performed statistical analyses and helped to draft the manuscript. AH provided diagnoses for the patients and enrolled patients eligible for the study. MM participated in preparation of the study design. JO conceived of the study, and participated in its design and coordination and drafted the manuscript. All authors read and approved the final manuscript.</p> <title>Availability &amp; requirements</title> <p><ext-link ext-link-type="uri" xlink:href="http://www.hutchon.net/ConfidOR.htm"/></p> <p><ext-link ext-link-type="uri" xlink:href="http://www.broad.mit.edu/mpg/haploview/"/></p> <title>Pre-publication history</title> <p>The pre-publication history for this paper can be accessed here:</p> <p><ext-link ext-link-type="uri" xlink:href="http://www.biomedcentral.com/1471-2350/9/81/prepub"/></p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We would like to thank MSc. Karolina Hanusek for cooperation in performing molecular assays.</p> <p>This work was supported by the CMKP grant 501-1-1-09-17/05.</p>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Allelic distribution in patients and healthy controls for all polymorphisms</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Polymorphism</td><td align="left">Genotype</td><td align="left">CD</td><td align="left">PSC</td><td align="left">PBC</td><td align="left">Controls</td></tr></thead><tbody><tr><td align="left">NOD2/CARD15 Pro268Ser</td><td align="left">C/C</td><td align="left">0.33</td><td align="left">0.48</td><td align="left">0.58</td><td align="left">0.56</td></tr><tr><td/><td align="left">C/T</td><td align="left">0.40</td><td align="left">0.45</td><td align="left">0.36</td><td align="left">0,41</td></tr><tr><td/><td align="left">T/T</td><td align="left">0.27</td><td align="left">0.06</td><td align="left">0.06</td><td align="left">0.03</td></tr><tr><td align="left">NOD2/CARD15 Arg702Trp</td><td align="left">C/C</td><td align="left">0.83</td><td align="left">0.92</td><td align="left">0.93</td><td align="left">0.97</td></tr><tr><td/><td align="left">C/T</td><td align="left">0.17</td><td align="left">0.08</td><td align="left">0.07</td><td align="left">0.03</td></tr><tr><td/><td align="left">T/T</td><td align="left">0.00</td><td align="left">0.00</td><td align="left">0.00</td><td align="left">0.00</td></tr><tr><td align="left">NOD2/CARD15 Gly908Arg</td><td align="left">G/G</td><td align="left">0.95</td><td align="left">0.96</td><td align="left">0.97</td><td align="left">0.94</td></tr><tr><td/><td align="left">G/C</td><td align="left">0.05</td><td align="left">0.04</td><td align="left">0.03</td><td align="left">0.06</td></tr><tr><td/><td align="left">C/C</td><td align="left">0.00</td><td align="left">0.00</td><td align="left">0.00</td><td align="left">0.00</td></tr><tr><td align="left">NOD2/CARD15 1007fs</td><td align="left">-/-</td><td align="left">0.63</td><td align="left">0.91</td><td align="left">0.92</td><td align="left">0.94</td></tr><tr><td/><td align="left">-/C</td><td align="left">0.27</td><td align="left">0.08</td><td align="left">0.08</td><td align="left">0.06</td></tr><tr><td/><td align="left">C/C</td><td align="left">0.10</td><td align="left">0.01</td><td align="left">0.00</td><td align="left">0.00</td></tr><tr><td align="left">OCTN1 Leu503Phe</td><td align="left">C/C</td><td align="left">0.27</td><td align="left">0.35</td><td align="left">0.37</td><td align="left">0.37</td></tr><tr><td/><td align="left">C/T</td><td align="left">0.50</td><td align="left">0.52</td><td align="left">0.43</td><td align="left">0.46</td></tr><tr><td/><td align="left">T/T</td><td align="left">0.23</td><td align="left">0.13</td><td align="left">0.20</td><td align="left">0.17</td></tr><tr><td align="left">OCTN2-207G&gt;C</td><td align="left">G/G</td><td align="left">0.25</td><td align="left">0.28</td><td align="left">0.30</td><td align="left">0.27</td></tr><tr><td/><td align="left">G/C</td><td align="left">0.47</td><td align="left">0.54</td><td align="left">0.43</td><td align="left">0.46</td></tr><tr><td/><td align="left">C/C</td><td align="left">0.27</td><td align="left">0.18</td><td align="left">0.27</td><td align="left">0.27</td></tr><tr><td align="left">DLG5 Arg30Gln</td><td align="left">C/C</td><td align="left">0.80</td><td align="left">0.77</td><td align="left">0.82</td><td align="left">0.80</td></tr><tr><td/><td align="left">C/T</td><td align="left">0.18</td><td align="left">0.21</td><td align="left">0.17</td><td align="left">0.19</td></tr><tr><td/><td align="left">T/T</td><td align="left">0.02</td><td align="left">0.03</td><td align="left">0.01</td><td align="left">0.01</td></tr><tr><td align="left">IL23R Arg381Gln</td><td align="left">G/G</td><td align="left">0.97</td><td align="left">0.90</td><td align="left">0.92</td><td align="left">0.95</td></tr><tr><td/><td align="left">G/A</td><td align="left">0.03</td><td align="left">0.09</td><td align="left">0.08</td><td align="left">0.05</td></tr><tr><td/><td align="left">A/A</td><td align="left">0.00</td><td align="left">0.01</td><td align="left">0.00</td><td align="left">0.00</td></tr><tr><td align="left">IL23R His3Gln</td><td align="left">C/C</td><td align="left">0.19</td><td align="left">0.26</td><td align="left">0.23</td><td align="left">0.30</td></tr><tr><td/><td align="left">C/A</td><td align="left">0.59</td><td align="left">0.57</td><td align="left">0.54</td><td align="left">0.49</td></tr><tr><td/><td align="left">A/A</td><td align="left">0,22</td><td align="left">0.17</td><td align="left">0.23</td><td align="left">0.22</td></tr><tr><td align="left">IL23R exon-3'UTR</td><td align="left">C/C</td><td align="left">0.40</td><td align="left">0.53</td><td align="left">0.51</td><td align="left">0.49</td></tr><tr><td/><td align="left">C/A</td><td align="left">0.52</td><td align="left">0.41</td><td align="left">0.40</td><td align="left">0.42</td></tr><tr><td/><td align="left">A/A</td><td align="left">0.08</td><td align="left">0.07</td><td align="left">0.08</td><td align="left">0.09</td></tr><tr><td align="left">ATG16L1 Thr300Ala</td><td align="left">C/C</td><td align="left">0.39</td><td align="left">0.29</td><td align="left">0.25</td><td align="left">0.23</td></tr><tr><td/><td align="left">C/T</td><td align="left">0.42</td><td align="left">0.51</td><td align="left">0.52</td><td align="left">0.50</td></tr><tr><td/><td align="left">T/T</td><td align="left">0.19</td><td align="left">0.20</td><td align="left">0.23</td><td align="left">0.27</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Statistical analyses for significant genotype associations</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Patient group</td><td align="center" colspan="2">SNP</td><td align="center">Genotype</td><td align="center">OR</td><td align="center">95% CI</td><td align="center">p-value</td><td align="center">Stat.<break/> Power</td></tr></thead><tbody><tr><td align="center">Crohn's<break/> disease</td><td align="center" colspan="1">NOD2/CARD15</td><td align="center">Pro268Ser</td><td align="center">C/T+T/T <break/>T/T</td><td align="center">2.52 <break/>12.18</td><td align="center">1.34 – 4.75 <break/>3.86 – 38.37</td><td align="center">0.005 <break/>0.0000003</td><td align="center">0.59 <break/>0.98</td></tr><tr><td/><td align="center">NOD2/CARD15</td><td align="center">Arg702Trp</td><td align="center">C/T+T/T <break/>T/T</td><td align="center">6.65 <break/>n/a</td><td align="center">1.99 – 22.17<break/> n/a</td><td align="center">0.0013 <break/>n/a</td><td align="center">0.71 <break/>n/a</td></tr><tr><td/><td align="center">NOD2/CARD15</td><td align="center">1007fs</td><td align="center">-/C+C/C <break/>C/C</td><td align="center">9.59<break/> &gt;16</td><td align="center">3.94 – 23.29<break/> ---</td><td align="center">0.00000002 <break/>0.0006</td><td align="center">0.99<break/> 0.86</td></tr><tr><td/><td align="center">ATG16L1</td><td align="center">Thr300Ala</td><td align="center">T/C+T/T <break/>T/T</td><td align="center">0.468 <break/>0.632</td><td align="center">0.24 – 0.90 <break/>0.30 – 1.34</td><td align="center">0.022 <break/>0.231</td><td align="center">0.32 <break/>0.05</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Statistical analyses using Haploview software (ver. 3.2) for significant haplotype associations.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Patient group</td><td align="center">SNP</td><td align="center">Haplotype</td><td align="center">OR</td><td align="center">95% CI</td><td align="center">p-value</td></tr></thead><tbody><tr><td align="center">Crohn's disease</td><td align="center">NOD2/CARD15 Pro268Ser/Arg702Trp</td><td align="center">TC</td><td align="center">2.25</td><td align="center">1.41 – 3.58</td><td align="center">0.0005</td></tr><tr><td/><td align="center">NOD2/CARD15 Pro268Ser/Gly908Arg</td><td align="center">TG</td><td align="center">3.01</td><td align="center">1.89 – 4.78</td><td align="center">0.000002</td></tr><tr><td/><td align="center">NOD2/CARD15 Pro268Ser/1007fs</td><td align="center">TC</td><td align="center">10.26</td><td align="center">4.52 – 23.29</td><td align="center">0.0000000000671</td></tr><tr><td/><td align="center">OCTN1/OCTN2</td><td align="center">CC</td><td align="center">0.28</td><td align="center">0.08 – 0.94</td><td align="center">0.0298</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Power of genotype association statistics given as 1-β error probability. The β represents the probability of falsely accepting H0 hypothesis (lack of association) when in fact H1 hypothesis (association) is true.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Polymorphism</td><td align="center">Genet. model</td><td align="center">CD</td><td align="center">PSC</td><td align="center">PBC</td></tr></thead><tbody><tr><td align="left">NOD2/CARD15 Pro268Ser</td><td align="center">dominant</td><td align="center">0.59</td><td align="center">0.04</td><td align="center">&lt;0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">0.98</td><td align="center">0.03</td><td align="center">0.04</td></tr><tr><td align="left">NOD2/CARD15 Arg702Trp</td><td align="center">dominant</td><td align="center">0.71</td><td align="center">0.11</td><td align="center">0.08</td></tr><tr><td/><td align="center">recessive</td><td align="center">n/a</td><td align="center">&lt;0.01</td><td align="center">n/a</td></tr><tr><td align="left">NOD2/CARD15 Gly908Arg</td><td align="center">dominant</td><td align="center">&lt;0.01</td><td align="center">&lt;0.01</td><td align="center">0.05</td></tr><tr><td/><td align="center">recessive</td><td align="center">n/a</td><td align="center">n/a</td><td align="center">n/a</td></tr><tr><td align="left">NOD2/CARD15 1007fs</td><td align="center">dominant</td><td align="center">0.99</td><td align="center">0.02</td><td align="center">0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">0.86</td><td align="center">&lt;0.01</td><td align="center">n/a</td></tr><tr><td align="left">OCTN1 Leu503Phe</td><td align="center">dominant</td><td align="center">0.08</td><td align="center">&lt;0.01</td><td align="center">&lt;0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">0.04</td><td align="center">0.02</td><td align="center">0.02</td></tr><tr><td align="left">OCTN2-207G&gt;C</td><td align="center">dominant</td><td align="center">&lt;0.01</td><td align="center">&lt;0.01</td><td align="center">0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">&lt;0.01</td><td align="center">0.09</td><td align="center">&lt;0.01</td></tr><tr><td align="left">DLG5 Arg30Gln</td><td align="center">dominant</td><td align="center">&lt;0.01</td><td align="center">0.02</td><td align="center">0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">&lt;0.01</td><td align="center">0.02</td><td align="center">&lt;0.01</td></tr><tr><td align="left">IL23R Arg381Gln</td><td align="center">dominant</td><td align="center">&lt;0.01</td><td align="center">0.08</td><td align="center">0.03</td></tr><tr><td/><td align="center">recessive</td><td align="center">&lt;0.01</td><td align="center">&lt;0.01</td><td align="center">n/a</td></tr><tr><td align="left">IL23R His3Gln</td><td align="center">dominant</td><td align="center">0.09</td><td align="center">0.01</td><td align="center">0.06</td></tr><tr><td/><td align="center">recessive</td><td align="center">&lt;0.01</td><td align="center">0.02</td><td align="center">&lt;0.01</td></tr><tr><td align="left">IL23R exon-3'UTR</td><td align="center">dominant</td><td align="center">0.05</td><td align="center">0.01</td><td align="center">0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">&lt;0.01</td><td align="center">0.01</td><td align="center">&lt;0.01</td></tr><tr><td align="left">ATG16L1 Thr300Ala</td><td align="center">dominant</td><td align="center">0.32</td><td align="center">0.04</td><td align="center">0.01</td></tr><tr><td/><td align="center">recessive</td><td align="center">0.05</td><td align="center">0.04</td><td align="center">0.02</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>genotype counts. genotype counts for each of the studied groups</p></caption></supplementary-material>
<media xlink:href="1471-2350-9-81-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Schrumpf", "Boberg"], "given-names": ["E", "KM"], "article-title": ["Epidemiology of primary sclerosing cholangitis"], "source": ["Best practice & research"], "year": ["2001"], "volume": ["15"], "fpage": ["553"], "lpage": ["562"]}, {"collab": ["Dobrowolska-Zachwieja"], "article-title": ["The sequence variant of NOD2/CARD15 in a Polish family on the background of Polish patients with Crohn\u2019s disease"], "source": ["Gastroenterologia Polska"], "year": ["2004"], "volume": ["11"], "fpage": ["325"], "lpage": ["331"]}, {"collab": ["Lakatos"], "article-title": ["NOD2/CARD15 SNP8, 12 and 13 and other exon4 mutations and primary biliary cirrhosis (PBC) in Hungarian and Polish patients."], "source": ["Z Gastroenterol"], "year": ["2004"], "volume": ["42"]}]
BMC Med Genet. 2008 Aug 21; 9:81