Sequence variation, linkage disequilibrium and association with Crohn's disease on chromosome 5q31

Chromosome 5q31 contains a cluster of genes involved in immune response, including a 250 kb risk haplotype associated with Crohn's disease (CD) susceptibility. Recently, two functional variants in SLC22A4 and SLC22A5 (L503F and G-207C), encoding the cation transporters OCTN1 and OCTN2, were proposed as causal variants for CD, but with conflicting genetic evidence regarding their contribution. We investigated this locus by resequencing the coding regions of 10 genes in 24 CD cases and deriving a linkage disequilibrium (LD) map of the 27 single nucleotide polymorphisms (SNPs) detected. Ten SNPs representative of the LD groups observed, were tested for CD association. L503F in SLC22A4 was the only nonsynonymous SNP significantly associated with CD (P=0.003), but was not associated with disease in the absence of other markers of the 250 kb risk haplotype. Two other SNPs, rs11242115 in IRF1 and rs17166050 in RAD50, lying outside the 250 kb risk haplotype, also showed CD association (P=0.019 and P=0.0080, respectively). The RAD50 gene contains a locus control region regulating expression of the Th2 cytokine genes at this locus. Other as yet undiscovered SNPs in this region may therefore modulate gene expression and contribute to the risk of CD, and perhaps of other inflammatory phenotypes.     

Association of the organic cation transporter OCTN genes with Crohn's disease in the Spanish population

The SLC22A4 and SLC22A5 genes within the IBD5 risk locus encode the organic cation transporters OCTN1 and OCTN2. Two variants, 1672C>T in SLC22A4 and -207G>C in SLC22A5, were shown to alter these genes' functions and were identified as genetic susceptibility factors for Crohn's disease (CD). We pursued to check both putative etiologic variants in an independent population through a case-control study with 309 Spanish CD patients and 408 ethnically matched healthy subjects. Both polymorphisms were found in partial linkage disequilibrium (D'=0.86). The separate analysis of each OCTN variant evidenced no association. However, when the simultaneous presence of mutant variants in both genes was analyzed, an effect on CD susceptibility was observed (P=0.026, odds ratio (OR) (95% confidence interval (CI))=1.59 (1.03-2.45)). The previously described predisposition conferred by the 5q31-risk haplotype increased in the absence of the etiologic 1672T and -207C alleles (P=0.0006, OR (95% CI)=10.14 (1.97-98.04)). Moreover, the risk contributed by these polymorphisms was higher in the IBD5 wild-type population (P=0.003, OR (95% CI)=2.65 (1.32-5.35)), arguing against the exclusive etiological role of the OCTN variants. The haplotype pattern inferred led to the consideration of these variants as susceptibility markers only in a defined genetic context. Our data support the interpretation of the 1672C>T SLC22A4 and -207G>C SLC22A5 polymorphisms as genetic markers of susceptibility/protection haplotypes.     

Direct or indirect association in a complex disease: the role of SLC22A4 and SLC22A5 functional variants in Crohn disease

A common haplotype spanning 250 kb on chromosome 5q31 is strongly associated with Crohn disease (CD). Recently, two functional variants within the SLC22A4 and SLC22A5 genes at this locus (IBD5), L503F (c.1507C > T) and G-207C (c.-207G > C), have been proposed to contribute directly to susceptibility to CD. However, extensive linkage disequilibrium at the IBD5 locus has complicated efforts to distinguish causal variants from association of the general risk haplotype. We genotyped the SLC22A4 and SLC22A5 variants and other polymorphisms across the risk haplotype in four populations of European origin, and applied regression-based haplotype analysis to over 1,200 fully genotyped case-control pairs, modeling case/control status on the presence of one or more SNPs to test for conditional association and to identify risk haplotypes. We found highly significant association of SNPs at the IBD5 locus with Crohn disease in all populations tested. However, the frequencies of L503F and G-207C in individuals who did not carry the general IBD5 risk haplotype were not significantly different in cases and controls, with associated disease odds ratios (ORs) of 0.90 (95% CI, 0.57-1.40) and 0.90 (95% CI, 0.65-1.23), respectively. Haplotype analysis showed that addition of the SLC22A4 and SLC22A5 variants to a null model that included the background risk haplotype did not significantly improve the model fit. In addition to the common risk haplotype, several rare haplotypes had an increased frequency in cases compared to controls. This study suggests that the molecular basis for Crohn disease susceptibility at the IBD5 locus remains to be defined, and highlights the challenge of the identification of causal variants in a complex disease in regions of extensive linkage disequilibrium.     

The CREM gene is involved in genetic predisposition to inflammatory bowel disease in the Tunisian population

The identification of susceptibility genes for inflammatory bowel disease (IBD) is key to understanding pathogenic mechanisms. Recently, the results of genetic association studies have highlighted many loci that are shared among several autoimmune diseases. We aimed to study the genetic epidemiology of polymorphisms in specific genes previously associated with other autoimmune diseases, namely the CREM, STAT4, STAT5a, Stat5b, and IRF5 genes. Twelve polymorphisms in the CREM, STAT4, STAT5a, Stat5b, and IRF5 genes were genotyped in a cohort of 107 IBD patients (39 Crohn's disease [CD] and 68 ulcerative colitis [UC]) and 162 controls from southern Tunisia. One CREM single nucleotide polymorphism (SNP) displayed evidence for genetic association with IBD (p=8.7×10(-4), odds ratio [OR]=2.84 [1.58; 5.09]). One STAT4 SNP (p=0.026; OR=1.65 [1.06; 2.58]) exhibited a marginal association with UC but not with CD. No significant association was observed with the SNPs in STAT5a, IRF5, and STAT5b. These results suggest that common variants of the CREM gene are involved in the genetic component conferring general susceptibility to IBD, whereas STAT4 appears to be more specifically associated with UC. This work provides motivation for studies aiming to replicate these findings in larger populations.     

Haplotype variation at the IBD5/SLC22A4 locus (5q31) in coeliac disease in the Irish population

In addition to the well-established association of coeliac disease (CD) with HLA-DQ (6p21) and possibly CTLA4 (2q33), there is considerable evidence for a susceptibility locus on chromosome 5q, which contains many potential candidates for inflammatory disease, including a cluster of cytokine genes in 5q31. CD cases and controls were genotyped for four single-nucleotide polymorphism (SNP) markers that together characterize >90% of the haplotype variation at the IBD5 locus encoding, among others, the SLC22A4 gene. IBD5 and SLC22A4 map to 5q31 and have recently been associated with Crohn's disease and rheumatoid arthritis. Haplotype frequencies do not differ significantly between CD cases and controls in the Irish population, and therefore the chromosome 5 CD susceptibility locus most likely lies elsewhere on 5q.     

Genome-wide association of serum uric acid concentration: replication of sequence variants in an island population of the Adriatic coast of Croatia

A genome-wide association study of serum uric acid (SUA) laevels was performed in a relatively isolated population of European descent from an island of the Adriatic coast of Croatia. The study sample included 532 unrelated and 768 related individuals from 235 pedigrees. Inflation due to relatedness was controlled by using genomic control. Genetic association was assessed with 2,241,249 single nucleotide polymorphisms (SNPs) in 1300 samples after adjusting for age and gender. Our study replicated four previously reported SUA loci (SLC2A9, ABCG2, RREB1, and SLC22A12). The strongest association was found with a SNP in SLC2A9 (rs13129697, P=2.33×10(-19)), which exhibited significant gender-specific effects, 35.76 μmol/L (P=2.11×10(-19)) in females and 19.58 μmol/L (P=5.40×10(-5)) in males. Within this region of high linkage disequilibrium, we also detected a strong association with a nonsynonymous SNP, rs16890979 (P=2.24×10(-17)), a putative causal variant for SUA variation. In addition, we identified several novel loci suggestive of association with uric acid levels (SEMA5A, TMEM18, SLC28A2, and ODZ2), although the P-values (P<5×10(-6)) did not reach the threshold of genome-wide significance. Together, these findings provide further confirmation of previously reported uric-acid-related genetic variants and highlight suggestive new loci for additional investigation.     

Catalog of 238 variations among six human genes encoding solute carriers ( hSLCs) in the Japanese population

We screened DNAs of 48 Japanese individuals for single-nucleotide polymorphisms (SNPs) in six genes encoding proteins of the solute carrier (SLC) family by direct sequencing of their entire genomic regions except for repetitive-sequence elements. This approach identified 213 SNPs and 25 insertion/deletion polymorphisms among the six genes. On average, we identified 1 SNP in every 509 nucleotides. Of the 213 SNPs, 14 were identified in the SLC10A1 gene, 51 in SLC15A1, 29 in SLC22A1, 27 in SLC22A2, 54 in SLC22A4, and 38 in SLC22A5. Eight were located in 5' flanking regions, 172 in introns, 25 in exons, and 8 in 3' flanking regions. These variants should contribute to investigations of possible correlations between genotypes and phenotypes as regards disease susceptibilities or responsiveness to drug therapy.     

Evidence for the association of the SLC22A4 and SLC22A5 genes with Type 1 Diabetes: a case control study

Background  

Type 1 diabetes (T1D) is a chronic, autoimmune and multifactorial disease characterized by abnormal metabolism of carbohydrate and fat. Diminished carnitine plasma levels have been previously reported in T1D patients and carnitine increases the sensitivity of the cells to insulin. Polymorphisms in the carnitine transporters, encoded by the SLC22A4 and SLC22A5 genes, have been involved in susceptibility to two other autoimmune diseases, rheumatoid arthritis and Crohn's disease. For these reasons, we investigated for the first time the association with T1D of six single nucleotide polymorphisms (SNPs) mapping to these candidate genes: slc2F2, slc2F11, T306I, L503F, OCTN2-promoter and OCTN2-intron.     

Pathway analysis of genome-wide association studies on uric acid concentrations

Objective

The aims of this study were to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms influencing uric acid level and to generate hypotheses for the SNP to gene to pathways that influence uric acid concentrations.

Methods

Meta-analysis data of 954 SNPs with genome-wide significance in 14 genome-wide association studies (GWASs) comprising 28,141 individuals of European ancestry was subjected to ICSNPathway (Identify candidate Causal SNPs and Pathways) analysis to establish associations between pathways and uric acid concentrations.

Results

ICSNPathway analysis identified 14 candidate causal SNPs, five genes, and two candidate causal pathways, which provided two hypothetical biologic mechanisms: (1) rs2728121 (regulatory region) to polycystic kidney disease 2 (PKD2) to ion transmembrane transporter activity; (2) rs942377, rs3799346, rs3799344, rs2762353, rs13197601, rs3757131, rs1165215, rs1165196 to SLC17A1 to ion transmembrane transporter activity and secondary active transmembrane transporter activity. SLC17A1, SLC17A3, SLC17A4, SLC22A11, and SLC2A9 were involved in both pathways, and PKD2 and SLC16A9 in one pathway.

Conclusion

By applying ICSNPathway analysis to GWAS data on uric acid levels, 14 SNPs, five genes, and two pathways involving the PKD2, SLC17A1, SLC17A3, SLC17A4, and SLC2A9 genes were identified that might contribute to the condition in Europeans.     

Catalog of 238 variations among six human genes encoding solute carriers (hSLCs) in the Japanese population

We screened DNAs of 48 Japanese individuals for single-nucleotide polymorphisms (SNPs) in six genes encoding proteins of the solute carrier (SLC) family by direct sequencing of their entire genomic regions except for repetitive-sequence elements. This approach identified 213 SNPs and 25 insertion/deletion polymorphisms among the six genes. On average, we identified 1 SNP in every 509 nucleotides. Of the 213 SNPs, 14 were identified in the SLC10A1 gene, 51 in SLC15A1, 29 in SLC22A1, 27 in SLC22A2, 54 in SLC22A4, and 38 in SLC22A5. Eight were located in 5′ flanking regions, 172 in introns, 25 in exons, and 8 in 3′ flanking regions. These variants should contribute to investigations of possible correlations between genotypes and phenotypes as regards disease susceptibilities or responsiveness to drug therapy. Received: July 24, 2002 / Accepted: July 25, 2002 Correspondence to:Y. Nakamura     

DLG5 variants contribute to Crohn disease risk in a Canadian population

Variants in the gene encoding the DLG5 scaffolding protein have been reported to be associated with increased risk of inflammatory bowel disease (IBD) and particularly Crohn's disease (Crohn disease; CD). These findings have not been uniformly replicated in follow-up studies. In this study we genotyped a cohort of 402 Canadian CD and 179 ulcerative colitis (UC) patients and 537 healthy controls for three IBD/CD-associated DLG5 variants. Our data reveal that the common DLG5 haplotype (A), which was previously considered protective for IBD, is associated with modest increases in risk for IBD (P=0.02) and CD (P=0.04). The effects of haplotype copy number on risk for IBD were minor, with the odds ratio (ORs) being 1.37 for the heterozygous risk genotype and 1.7 for the homozygous risk genotype. While we were unable to replicate the proposed association between the DLG5 c.113G>A variant and IBD, an association of IBD (P=0.02) and CD (P=0.04) with the rarer c.4136C>A variant was replicated in this cohort. These associations were restricted to the non-Jewish subjects in this cohort and were not detected in the Ashkenazi Jewish population studied here. Within the non-Jewish group, no associations were detected between the DLG5 variants and specific phenotypic features, such as site of disease, and there was no evidence of epistasis between DLG5 and any of the CD-associated CARD15 or SLC22A4/A5 gene variants. Together, the results indicate a role for DLG5 variants in IBD susceptibility and suggest that further studies are warranted to evaluate this role in different IBD populations and to determine the functional pathways that couple DLG5 variants to IBD.     

Evaluation of polymorphisms in predicted target sites for micro RNAs differentially expressed in endometriosis

Previous microarray analyses identified 22 microRNAs (miRNAs) differentially expressed in paired ectopic and eutopic endometrium of women with and without endometriosis. To investigate further the role of these miRNAs in women with endometriosis, we conducted an association study aiming to explore the relationship between endometriosis risk and single-nucleotide polymorphisms (SNPs) in miRNA target sites for these differentially expressed miRNAs. A panel of 102 SNPs in the predicted miRNA binding sites were evaluated for an endometriosis association study and an ingenuity pathway analysis was performed. Fourteen rare variants were identified in this study. We found SNP rs14647 in the Wolf-Hirschhorn syndrome candidate gene1 (WHSC1) 3'UTR (untranslated region) was associated with endometriosis-related infertility presenting an odds ratio of 12.2 (95% confidence interval = 2.4-60.7, P = 9.03 × 10(-5)). SNP haplotype AGG in the solute carrier family 22, member 23 (SLC22A23) 3'UTR was associated with endometriosis-related infertility and more severe disease. With the individual genotyping data, ingenuity pathways analysis identified the tumour necrosis factor and cyclin-dependant kinase inhibitor as major factors in the molecular pathways. Significant associations between WHSC1 alleles and endometriosis-related infertility and SLC22A23 haplotypes and the disease severe stage were identified. These findings may help focus future research on subphenotypes of this disease. Replication studies in independent large sample sets to confirm and characterize the involvement of the gene variation in the pathogenesis of endometriosis are needed.     

Comprehensive evaluation of the genetic variants of interferon regulatory factor 5 (IRF5) reveals a novel 5 bp length polymorphism as strong risk factor for systemic lupus erythematosus

We analyzed a comprehensive set of single-nucleotide polymorphisms (SNPs) and length polymorphisms in the interferon regulatory factor 5 (IRF5) gene for their association with the autoimmune disease systemic lupus erythematosus (SLE) in 485 Swedish patients and 563 controls. We found 16 SNPs and two length polymorphisms that display association with SLE (P < 0.0005, OR > 1.4). Using a Bayesian model selection and averaging approach we identified parsimonious models with exactly two variants of IRF5 that are independently associated with SLE. The variants of IRF5 with the highest posterior probabilities (1.00 and 0.71, respectively) of being causal in SLE are a SNP (rs10488631) located 3' of IRF5, and a novel CGGGG insertion-deletion (indel) polymorphism located 64 bp upstream of the first untranslated exon (exon 1A) of IRF5. The CGGGG indel explains the association signal from multiple SNPs in the IRF5 gene, including rs2004640, rs10954213 and rs729302 previously considered to be causal variants in SLE. The CGGGG indel contains three or four repeats of the sequence CGGGG with the longer allele containing an additional SP1 binding site as the risk allele for SLE. Using electrophoretic mobility shift assays we show increased binding of protein to the risk allele of the CGGGG indel and using a minigene reporter assay we show increased expression of IRF5 mRNA from a promoter containing this allele. Increased expression of IRF5 protein was observed in peripheral blood mononuclear cells from SLE patients carrying the risk allele of the CGGGG indel. We have found that the same IRF5 allele also confers risk for inflammatory bowel diseases and multiple sclerosis, suggesting a general role for IRF5 in autoimmune diseases.     

CD226 Gly307Ser association with multiple autoimmune diseases

Genome-wide association studies provide insight into multigenic diseases through the identification of susceptibility genes and etiological pathways. In addition, the identification of shared variants among autoimmune disorders provides insight into common disease pathways. We previously reported an association of a nonsynonymous single nucleotide polymorphism (SNP) rs763361/Gly307Ser in the immune response gene CD226 on chromosome 18q22 with type 1 diabetes (T1D) susceptibility. Here, we report efforts toward identifying the causal variant by exonic resequencing and tag SNP mapping of the 18q22 region in both T1D and multiple sclerosis (MS). In addition to the analysis of newly available samples in T1D (2088 cases and 3289 controls) and autoimmune thyroid disease (AITD) (821 cases and 1920 controls), resulting in strong support for the Ser(307) association with T1D (P=3.46 x 10(-9)) and continued potential evidence for AITD (P=0.0345), we provide evidence for association of Gly307Ser with MS (P=4.20 x 10(-4)) and rheumatoid arthritis (RA) (P=0.017). The Ser(307) allele of rs763361 in exon 7 of CD226 predisposes to T1D, MS, and possibly AITD and RA, and based on the tag SNP analysis, could be the causal variant.     

Association analysis of <Emphasis Type="Italic">SLC22A4</Emphasis>, <Emphasis Type="Italic">SLC22A5</Emphasis> and <Emphasis Type="Italic">DLG5</Emphasis> in Japanese patients with Crohn disease

Crohn disease (CD) is an inflammatory bowel disease characterized by chronic transmural, segmental, and typically granulomatous inflammation of the gut. Recently, two novel candidate gene loci associated with CD, SLC22A4 and SLC22A5 on chromosome 5 known as IBD5 and DLG5 on chromosome 10, were identified through association analysis of Caucasian CD patients. We validated these candidate genes in Japanese patients with CD and found a weak but possible association with both SLC22A4 (P=0.028) and DLG5 (P=0.023). However, the reported genetic variants that were indicated to be causative in the Caucasian population were completely absent in or were not associated with Japanese CD patients. These findings imply significant differences in genetic background with CD susceptibility among different ethnic groups and further indicate some difficulty of population-based studies.     

The 5q31 variants associated with psoriasis and Crohn's disease are distinct

Predisposition to psoriasis is known to be affected by genetic variation in HLA-C, IL12B and IL23R, but other genetic risk factors also exist. We recently reported three psoriasis-associated single nucleotide polymorphisms (SNPs) in the 5q31 locus, a region of high linkage disequilibrium laden with inflammatory pathway genes. The aim of this study was to assess whether other variants in the 5q31 region are causal to these SNPs or make independent contributions to psoriasis risk by genotyping a comprehensive set of tagging SNPs in a 725 kb region bounded by IL3 and IL4 and testing for disease association. Ninety SNPs, capturing 86.4% of the genetic diversity, were tested in one case-control sample set (467 cases/460 controls) and significant markers (P(allelic) < 0.05) (n = 9) were then tested in two other sample sets (981 cases/925 controls). All nine SNPs were significant in a meta-analysis of the combined sample sets. Pair-wise conditional association tests showed rs1800925, an intergenic SNP located just upstream of IL13 (Mantel-Haenszel P(combined) = 1.5 x 10(-4), OR = 0.77 [0.67-0.88]), could account for observed significant association of all but one other SNP, rs11568506 in SLC22A4 [Mantel-Haenszel P(combined) = 0.043, OR = 0.68 (0.47-0.99)]. Haplotype analysis of these two SNPs showed increased significance for the two common haplotypes (rs11568506-rs1800925: GC, P(combined) = 5.67 x 10(-6), OR = 1.37; GT, P(combined) = 6.01 x 10(-5), OR = 0.75; global haplotype P = 8.93 x 10(-5)). Several 5q31-region SNPs strongly associated with Crohn's disease (CD) in the recent WTCCC study were not significant in the psoriasis sample sets tested here. These results identify the most significant 5q31 risk variants for psoriasis and suggest that distinct 5q31 variants contribute to CD and psoriasis risk.     

Investigation of NOD1/CARD4 variation in inflammatory bowel disease using a haplotype-tagging strategy

Both NOD1/CARD4 and NOD2/CARD15 are intracellular pattern-recognition receptors involved in the innate immune response. Germline NOD2/CARD15 variation has a definite effect on susceptibility to Crohn's disease (CD) and phenotype, although this contribution is weak in Scotland and Scandinavia. The NOD1/CARD4 gene lies within the putative inflammatory bowel disease (IBD) locus at 7p14.3. We have assessed, in detail, the influence of germline NOD1/CARD4 variation on IBD susceptibility and phenotype in the Scottish population. Two thousand two hundred and ninety-six subjects, including 356 children with IBD, were involved in parallel case-control and family-based association studies. Nine tagging single-nucleotide polymorphisms (SNPs) were selected based on HapMap data spanning the whole of the NOD1/CARD4 gene. Our case-control SNP analysis was powered to detect an effect size with OR 1.5 for IBD and OR 1.6 for CD. No significant associations were observed between any of nine the NOD1/CARD4 SNPs studied and IBD, CD or ulcerative colitis (UC) (P > 0.05 for all). Haplotype case-control analysis was also negative (P > 0.05 in IBD, CD and UC). Multimarker transmission disequilibrium testing analysis was negative (P > 0.05 in IBD, CD and UC). NOD2/CARD15 variant carriage had no influence on NOD1/CARD4 effect on IBD susceptibility. This study represents the first application of a gene -wide haplotype-tagging approach to assess, in detail, the contribution of NOD1/CARD4 in IBD.