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.     

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.     

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.     

Refined genomic localization and ethnic differences observed for the IBD5 association with Crohn's disease

Although the general association of the inflammatory bowel disease (IBD) 5 region on chromosome 5q31 to Crohn's disease (CD) has been replicated repeatedly, the identity of the precise causal variant within the region remains unknown. A recent report proposed polymorphisms in solute carrier family 22, member 4 (SLC22A4) organic cation transporter 1(OCTN1) and solute carrier family 22, member 5 (SLC22A5) (OCTN2) as responsible for the IBD5 association, but definitive, large-sample comparison of those polymorphisms with others known to be in strong linkage disequilibrium was not performed. We evaluated 1879 affected offspring and parents ascertained by a North American IBD Genetics Consortium for six IBD5 tag single nucleotide polymorphisms (SNPs) to evaluate association localization and ethnic and subphenotypic specificity. We confirm association to the IBD5 region (best SNP IGR2096a_1/rs12521868, P<0.0005) and show this association to be exclusive to the non-Jewish (NJ) population (P=0.00005) (risk allele undertransmitted in Ashkenazi Jews). Using Phase II HapMap data, we demonstrate that there are a set of polymorphisms, spanning genes from prolyl 4-hydroxylase (P4HA2) through interferon regulatory factor 1 (IRF1) with equivalent statistical evidence of association to the reported SLC22A4 variant and that each, by itself, could entirely explain the IBD5 association to CD. Additionally, the previously reported SLC22A5 SNP is rejected as the potential causal variant. No specificity of association was seen with respect to disease type and location, and a modest association to ulcerative colitis is also observed. We confirm the importance of IBD5 to CD susceptibility, demonstrate that the locus may play a role in NJ individuals only, and establish that IRF1, PDLIM, and P4HA2 may be equally as likely to contain the IBD5 causal variant as the OCTN genes.     

Evidence for unique association signals in SLE at the CD28-CTLA4-ICOS locus in a family-based study

CD28, CTLA4 (cytotoxic T lymphocyte-associated protein 4) and ICOS (inducible T cell co-stimulator) are good candidate genes for systemic lupus erythematosus (SLE) because of their role in regulating T cell activation. CTLA4 inhibits CD28-mediated T cell activation. CTLA4 is expressed on CD4+ and CD8+ activated T cells, and also B cells, but CD28 and ICOS are largely restricted to T cells. An interval encompassing the CD28-CTLA4-ICOS locus on chromosome 2q33 was linked to lupus in two genome-wide linkage scans. This large family-based association study in 532 UK SLE families represents the first high-density genetic screen of 80 SNPs at this locus. There are seven haplotype blocks across the locus. In CTLA4, the strongest signal comes from two variants, located 2.1 kb downstream from the 3'-UTR. These polymorphisms, rs231726 (SNP 43) and rs231726 (SNP 44), are in complete linkage disequilibrium (LD) (r(2)=1) and are associated with SLE P=0.0008 (GH) and P=0.01 (family-based association test). There is also a signal in the distal 3' flanking region of CTLA4/ICOS promoter (P=0.003). There was no confirmation of published associations for SLE in the promoter or coding region of CTLA4. These SLE risk alleles are more distal than those identified in Graves' disease and are in LD with Graves' disease protective alleles identified in both of these regions of CTLA4 (Ueda et al. 2003). These factors suggest an SLE-specific pattern of association. The functional consequences of the associated polymorphisms are likely to influence CTLA4 expression, although it is possible that genetically modulated ICOS expression is involved in SLE susceptibility.     

Association between variants in IDE-KIF11-HHEX and plasma amyloid β levels

Genetic linkage and association studies in late-onset Alzheimer’s disease (LOAD) or its endophenotypes have pointed to several regions on chromosome 10q, among these the ∼ 250 kb linkage disequilibrium (LD) block harboring the genes IDE, KIF1, and HHEX. We explored the association between variants in the genomic region harboring the IDE-KIF11-HHEX complex with plasma Aβ40 and Aβ42 levels in a case-control cohort of Caribbean Hispanics. First, we performed single marker linear regression analysis relating the individual single nucleotide polymorphisms (SNPs) with plasma Aβ40 and Aβ42 levels. Then we performed 3-SNP sliding window haplotype analyses, correcting all analyses for multiple testing. Out of 32 SNPs in this region, 3 SNPs in IDE (rs2421943, rs12264682, rs11187060) were associated with plasma Aβ40 or Aβ42 levels in single marker and haplotype analyses after correction for multiple testing. All these SNPs lie within the same LD block, and are in LD with the previously reported haplotypes. Our findings provide support for an association in the IDE region on chromosome 10q with Aβ40 and 42 levels.     

Identification of independent risk loci for Graves' disease within the MHC in the Japanese population

To identify genetic variants that confer the risk of Graves' disease (GD) in the Japanese population, we conducted a two-stage genome-wide association study (GWAS) using 1119 Japanese individuals with GD and 2718 unrelated controls, and a subsequent replication study using independent 432 GD cases and 1157 controls. We identified 34 single nucleotide polymorphisms (SNPs) to be significantly associated with GD in the GWAS phase. Twenty-two out of 34 SNPs remained positive in the replication study. All 22 SNPs were located within the major histocompatibility complex (MHC) locus on chromosome 6p21. No strong long-range linkage disequilibrium (LD) was observed among the 22 SNPs, indicating independent involvement of multiple loci within the MHC with the risk of GD. Multivariate stepwise logistic regression analysis selected rs3893464, rs4313034, rs3132613, rs4248154, rs2273017, rs9394159 and rs4713693, as markers for independent risk loci for GD. The analysis of LD between these seven SNPs and tagging SNPs for GD-associated human leukocyte antigen (HLA) alleles in the Japanese population (HLA-DPB1(*)0501 and HLA-A(*)0206) demonstrated that all of and five of seven SNPs were not in strong LD with HLA-DPB1(*)0501 and HLA-A(*)0206, respectively. Although causal variants remain to be identified, our results demonstrate the existence of multiple GD susceptibility loci within the MHC region.     

Haplotype patterns in cancer-related genes with long-range linkage disequilibrium: no evidence of association with breast cancer or positive selection

The average length of linkage disequilibrium (LD) blocks in European populations is about 22 kb. In this study, we have selected 20 genes with LD blocks larger than 60 kb (with a median length of 88 kb) from a total of 121 cancer-related genes. We observed limited haplotype diversity, with an average of three haplotypes per gene accounting for more than 90% of the diversity, two of these being a Yin-Yang pair in 95% of the LD blocks. The mean frequency of the most common haplotype in the Spanish population was just below 50%, similar to those for the HapMap CEU and African samples, but lower than the 60% observed in Asian samples. Genes involved in the regulation of nucleobases and nucleic acid metabolism were overrepresented among these 20 genes with long LD blocks (eight genes ATM, BRCA1, BRCA2, ERCC6, MLH1, MSH3, RAD54B and XRCC4) relative to the other 101 cancer-related genes studied (P=1.23 x 10(-6)). The ancestral haplotype was observed at a frequency greater than 3 in 67% of the genes either in the Spanish or one of the HapMap sampled populations. When observed, the ancestral haplotype had an average 15% frequency in the Spanish sample, less than half that observed in Asian and African samples. The Spanish Yin-Yang haplotype pair represented over 35% of haplotypes in African samples and over 65% in non-African samples. We detected differences in SNP frequencies between populations for five genes (ALDH2, APC, PIK3CB, RB1 and XRCC4, all with Fst>0.4); however, these genes did not show evidence of positive selection. Finally, we found no evidence that the haplotypes formed by SNPs in the 20 genes are associated with breast cancer.     

Association of the TSHR gene with Graves' disease: the first disease specific locus

The development of autoimmune thyroid disease (AITD) is associated with autoantibodies directed against the thyroid stimulating hormone receptor (TSHR). Previous studies have failed to demonstrate a consistent association between the TSHR and AITD, or any of its sub-phenotypes. In the present study, we analysed the linkage disequilibrium (LD) structure encompassing the TSHR, to identify LD 'blocks' and SNPs, which capture the majority of intra-block haplotype diversity. The haplotype tagging SNPs, plus all common SNPs reported in previous studies were genotyped in 1,059 AITD Caucasian cases and 971 Caucasian controls. A haplotype, across two LD blocks, showed association (P<1 x 10(-6), OR 1.7) with Graves' disease (GD) but not autoimmune hypothyroidism (AIH). We replicated these findings by genotyping the most associated GD SNP, rs2268458, in a separate UK Caucasian cohort of 1,366 AITD cases and 1,061 controls (GD, P=2 x 10(-6), OR 1.3; AIH, P=NS). These results in two independent Caucasian data sets suggest that the TSHR is the first replicated GD-specific locus meriting further fine mapping and functional analysis to identify the aetiological variants.     

Clinical and molecular characterization of Wilson's disease in China: identification of 14 novel mutations

Background

Abdominal aortic aneurysm (AAA) is a complex disorder with multiple genetic risk factors. Using affected relative pair linkage analysis, we previously identified an AAA susceptibility locus on chromosome 19q13. This locus has been designated as the AAA1 susceptibility locus in the Online Mendelian Inheritance in Man (OMIM) database.

Methods

Nine candidate genes were selected from the AAA1 locus based on their function, as well as mRNA expression levels in the aorta. A sample of 394 cases and 419 controls was genotyped for 41 SNPs located in or around the selected nine candidate genes using the Illumina GoldenGate platform. Single marker and haplotype analyses were performed. Three genes (CEBPG, PEPD and CD22) were selected for DNA sequencing based on the association study results, and exonic regions were analyzed. Immunohistochemical staining of aortic tissue sections from AAA and control individuals was carried out for the CD22 and PEPD proteins with specific antibodies.

Results

Several SNPs were nominally associated with AAA (p < 0.05). The SNPs with most significant p-values were located near the CCAAT enhancer binding protein (CEBPG), peptidase D (PEPD), and CD22. Haplotype analysis found a nominally associated 5-SNP haplotype in the CEBPG/PEPD locus, as well as a nominally associated 2-SNP haplotype in the CD22 locus. DNA sequencing of the coding regions revealed no variation in CEBPG. Seven sequence variants were identified in PEPD, including three not present in the NCBI SNP (dbSNP) database. Sequencing of all 14 exons of CD22 identified 20 sequence variants, five of which were in the coding region and six were in the 3'-untranslated region. Five variants were not present in dbSNP. Immunohistochemical staining for CD22 revealed protein expression in lymphocytes present in the aneurysmal aortic wall only and no detectable expression in control aorta. PEPD protein was expressed in fibroblasts and myofibroblasts in the media-adventitia border in both aneurysmal and non-aneurysmal tissue samples.

Conclusions

Association testing of the functional positional candidate genes on the AAA1 locus on chromosome 19q13 demonstrated nominal association in three genes. PEPD and CD22 were considered the most promising candidate genes for altering AAA risk, based on gene function, association evidence, gene expression, and protein expression.     

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.     

Screening of genetic variations of SLC15A2, SLC22A1, SLC22A2 and SLC22A6 genes

A growing list of membrane-spanning proteins involved in the transport of a large variety of drugs has been recognized and characterized to include peptide and organic anion/cation transporters. Given such an important role of transporter genes in drug disposition process, the role of single-nucleotide polymorphisms (SNPs) in such transporters as potential determinants of interindividual variability in drug disposition and pharmacological response has been investigated. To define the distribution of transporter gene SNPs across ethnic groups, we screened 450 DNAs in cohorts of 250 Korean, 50 Han Chinese, 50 Japanese, 50 African-American and 50 European-American ancestries for 64 SNPs in four transporter genes encoding proteins of the solute carrier family (SLC15A2, SLC22A1, SLC22A2 and SLC22A6). Of the 64 SNPs, 19 were core pharmacogenetic variants and 45 were HapMap tagging SNPs. Polymorphisms were genotyped using the golden gate genotyping assay. After genetic variability, haplotype structures and ethnic diversity were analyzed, we observed that the distributions of SNPs in a Korean population were similar to other Asian groups (Chinese and Japanese), and significantly different from African-American and European-American cohorts. Findings from this study would be valuable for further researches, including pharmacogenetic studies for drug responses.     

Positive associations of polymorphisms in the metabotropic glutamate receptor type 8 gene (GRM8) with schizophrenia.

The glutamatergic dysfunction has been implicated in pathophysiology of schizophrenia. The Group III metabotropic glutamate receptor 4 (mGluR4), 6, 7, and 8 are thought to modulate glutamatergic transmission in the brain by inhibiting glutamate release at the synapse. We tested association of schizophrenia with GRM8 using 22 single nucleofide polymorphisms (SNPs) with the average intervals of 40.3 kb in the GRM8 region in 100 case-control pairs for the SNPs. Although we observed significant associations of schizophrenia with two SNPs, SNP18 (rs2237748, allele: P = 0.0279; genotype: P = 0.0124) and SNP19 (rs2299472, allele: P = 0.0302; genotype: P = 0.0127), none of two SNPs showed significant association with disease after Bonferroni correction. Both SNP18 and SNP19 were included in a large region (>330 kb) in which SNPs are in linkage disequilibrium (LD) at the 3' region of GRM8. We also tested haplotype association of schizophrenia with constructed haplotypes of the SNPs in LD. Significant associations were detected for the combinations of SNP5-SNP6 (chi(2) = 18.12, df = 3, P = 0.0004, P corr = 0.0924 with Bonferroni correction), SNP4-SNP5-SNP6 (chi(2) = 27.50, df = 7, P = 0.0075, P corr = 0.015 with Bonferroni correction), and SNP5-SNP6-SNP7 (chi(2) = 23.92, df = 7, P = 0.0011, P corr = 0.0022 with Bonferroni correction). Thus, we conclude that at least one susceptibility locus for schizophrenia is located within the GRM8 region in Japanese.     

Identification of genetic variants that influence circulating IGF1 levels: a targeted search strategy

An important class of genetic variants that affect disease susceptibility may lie within regulatory elements that influence gene expression. Regulatory sequences are difficult to identify and may be distant from the genes they regulate, but many lie within evolutionarily conserved regions (ECRs). We used comparative genomics to identify 12 ECRs up to 75 kb 5' to and within introns of IGF1. These were screened by high-resolution melting curve analysis, and 18 single-nucleotide polymorphisms (SNPs) were identified, including five novel variants. We analysed two large population-based series of healthy women to test the nine SNPs with minor allele frequency (MAF) >1% within ECRs. Three of the nine SNPs within ECRs (rs35455143, rs35765817 and rs3839984) were significantly associated with circulating IGF1 levels in a multivariate analysis (P 70 kb 5' of IGF1. Two (rs35455143 and rs35765817) are in strong LD with each other and appear to have opposite effects on circulating IGF1. Our results on a subset of other SNPs in or near IGF1 were consistent with previously reported associations with IGF1 levels, although only one (rs35767: P = 0.05) was statistically significant. We believe that this is the first systematic study of an association between a phenotype and SNPs within ECRs extending over a large region adjacent to a gene. Targeting ECRs appears to be a useful strategy for identifying a subset of potentially functional non-coding regulatory SNPs.     

Linkage disequilibrium and haplotype diversity in the genes of the renin-angiotensin system: findings from the family blood pressure program

Association studies of candidate genes with complex traits have generally used one or a few single nucleotide polymorphisms (SNPs), although variation in the extent of linkage disequilibrium (LD) within genes markedly influences the sensitivity and precision of association studies. The extent of LD and the underlying haplotype structure for most candidate genes are still unavailable. We sampled 193 blacks (African-Americans) and 160 whites (European-Americans) and estimated the intragenic LD and the haplotype structure in four genes of the renin-angiotensin system. We genotyped 25 SNPs, with all but one of the pairs spaced between 1 and 20 kb, thus providing resolution at small scale. The pattern of LD within a gene was very heterogeneous. Using a robust method to define haplotype blocks, blocks of limited haplotype diversity were identified at each locus; between these blocks, LD was lost owing to the history of recombination events. As anticipated, there was less LD among blacks, the number of haplotypes was substantially larger, and shorter haplotype segments were found, compared with whites. These findings have implications for candidate-gene association studies and indicate that variation between populations of European and African origin in haplotype diversity is characteristic of most genes.