Mutational characterization of the bile acid receptor TGR5 in primary sclerosing cholangitis

Background

TGR5, the G protein-coupled bile acid receptor 1 (GPBAR1), has been linked to inflammatory pathways as well as bile homeostasis, and could therefore be involved in primary sclerosing cholangitis (PSC) a chronic inflammatory bile duct disease. We aimed to extensively investigate TGR5 sequence variation in PSC, as well as functionally characterize detected variants.

Methodology/Principal Findings

Complete resequencing of TGR5 was performed in 267 PSC patients and 274 healthy controls. Six nonsynonymous mutations were identified in addition to 16 other novel single-nucleotide polymorphisms. To investigate the impact from the nonsynonymous variants on TGR5, we created a receptor model, and introduced mutated TGR5 constructs into human epithelial cell lines. By using confocal microscopy, flow cytometry and a cAMP-sensitive luciferase assay, five of the nonsynonymous mutations (W83R, V178M, A217P, S272G and Q296X) were found to reduce or abolish TGR5 function. Fine-mapping of the previously reported PSC and UC associated locus at chromosome 2q35 in large patient panels revealed an overall association between the TGR5 single-nucleotide polymorphism rs11554825 and PSC (odds ratio  = 1.14, 95% confidence interval: 1.03–1.26, p = 0.010) and UC (odds ratio  = 1.19, 95% confidence interval 1.11–1.27, p = 8.5×10⁻⁷), but strong linkage disequilibrium precluded demarcation of TGR5 from neighboring genes.

Conclusions/Significance

Resequencing of TGR5 along with functional investigations of novel variants provided unique insight into an important candidate gene for several inflammatory and metabolic conditions. While significant TGR5 associations were detected in both UC and PSC, further studies are needed to conclusively define the role of TGR5 variation in these diseases.

Independent and population-specific association of risk variants at the IRGM locus with Crohn's disease

DNA polymorphisms in a region on chromosome 5q33.1 which contains two genes, immunity related GTPase related family, M (IRGM) and zinc finger protein 300 (ZNF300), are associated with Crohn's disease (CD). The deleted allele of a 20 kb copy number variation (CNV) upstream of IRGM was recently shown to be in strong linkage disequilibrium (LD) with the CD-associated single nucleotide polymorphisms and is itself associated with CD (P < 0.01). The deletion was correlated with increased or reduced expression of IRGM in transformed cells in a cell line-dependent manner, and has been proposed as a likely causal variant. We report here that small insertion/deletion polymorphisms in the promoter and 5' untranslated region of IRGM are, together with the CNV, strongly associated with CD (P = 1.37 x 10(-5) to 1.40 x 10(-9)), and that the CNV and the 5'-untranslated region variant -308(GTTT)(5) contribute independently to CD susceptibility (P = 2.6 x 10(-7) and P = 2 x 10(-5), respectively). We also show that the CD risk haplotype is associated with a significant decrease in IRGM expression (P < 10(-12)) in untransformed lymphocytes from CD patients. Further analysis of these variants in a Japanese CD case-control sample and of IRGM expression in HapMap populations revealed that neither the IRGM insertion/deletion polymorphisms nor the CNV was associated with CD or with altered IRGM expression in the Asian population. This suggests that the involvement of the IRGM risk haplotype in the pathogenesis of CD requires gene-gene or gene-environment interactions which are absent in Asian populations, or that none of the variants analysed are causal, and that the true causal variants arose after the European-Asian split.

Common variants at five new loci associated with early-onset inflammatory bowel disease

The inflammatory bowel diseases (IBD) Crohn's disease and ulcerative colitis are common causes of morbidity in children and young adults in the western world. Here we report the results of a genome-wide association study in early-onset IBD involving 3,426 affected individuals and 11,963 genetically matched controls recruited through international collaborations in Europe and North America, thereby extending the results from a previous study of 1,011 individuals with early-onset IBD. We have identified five new regions associated with early-onset IBD susceptibility, including 16p11 near the cytokine gene IL27 (rs8049439, P = 2.41 x 10(-9)), 22q12 (rs2412973, P = 1.55 x 10(-9)), 10q22 (rs1250550, P = 5.63 x 10(-9)), 2q37 (rs4676410, P = 3.64 x 10(-8)) and 19q13.11 (rs10500264, P = 4.26 x 10(-10)). Our scan also detected associations at 23 of 32 loci previously implicated in adult-onset Crohn's disease and at 8 of 17 loci implicated in adult-onset ulcerative colitis, highlighting the close pathogenetic relationship between early- and adult-onset IBD.

Searching for genotype-phenotype structure: using hierarchical log-linear models in Crohn disease

There has been considerable recent success in the detection of gene-disease associations. We consider here the development of tools that facilitate the more detailed characterization of the effect of a genetic variant on disease. We replace the simplistic classification of individuals according to a single binary disease indicator with classification according to a number of subphenotypes. This more accurately reflects the underlying biological complexity of the disease process, but it poses additional analytical difficulties. Notably, the subphenotypes that make up a particular disease are typically highly associated, and it becomes difficult to distinguish which genes might be causing which subphenotypes. Such problems arise in many complex diseases. Here, we concentrate on an application to Crohn disease (CD). We consider this problem as one of model selection based upon log-linear models, fitted in a Bayesian framework via reversible-jump Metropolis-Hastings approach. We evaluate the performance of our suggested approach with a simple simulation study and then apply the method to a real data example in CD, revealing a sparse disease structure. Most notably, the associated NOD2.908G-->R mutation appears to be directly related to more severe disease behaviors, whereas the other two associated NOD2 variants, 1007L-->FS and 702R-->W, are more generally related to disease in the small bowel (ileum and jejenum). The ATG16L1.300T-->A variant appears to be directly associated with only disease of the small bowel.

Genome-wide association defines more than 30 distinct susceptibility loci for Crohn's disease

Several risk factors for Crohn's disease have been identified in recent genome-wide association studies. To advance gene discovery further, we combined data from three studies on Crohn's disease (a total of 3,230 cases and 4,829 controls) and carried out replication in 3,664 independent cases with a mixture of population-based and family-based controls. The results strongly confirm 11 previously reported loci and provide genome-wide significant evidence for 21 additional loci, including the regions containing STAT3, JAK2, ICOSLG, CDKAL1 and ITLN1. The expanded molecular understanding of the basis of this disease offers promise for informed therapeutic development.

Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn's disease

We report results of a nonsynonymous SNP scan for ulcerative colitis and identify a previously unknown susceptibility locus at ECM1. We also show that several risk loci are common to ulcerative colitis and Crohn's disease (IL23R, IL12B, HLA, NKX2-3 and MST1), whereas autophagy genes ATG16L1 and IRGM, along with NOD2 (also known as CARD15), are specific for Crohn's disease. These data provide the first detailed illustration of the genetic relationship between these common inflammatory bowel diseases.

Systematic association mapping identifies NELL1 as a novel IBD disease gene

Crohn disease (CD), a sub-entity of inflammatory bowel disease (IBD), is a complex polygenic disorder. Although recent studies have successfully identified CD-associated genetic variants, these susceptibility loci explain only a fraction of the heritability of the disease. Here, we report on a multi-stage genome-wide scan of 393 German CD cases and 399 controls. Among the 116,161 single-nucleotide polymorphisms tested, an association with the known CD susceptibility gene NOD2, the 5q31 haplotype, and the recently reported CD locus at 5p13.1 was confirmed. In addition, SNP rs1793004 in the gene encoding nel-like 1 precursor (NELL1, chromosome 11p15.1) showed a consistent disease-association in independent German population- and family-based samples (942 cases, 1082 controls, 375 trios). Subsequent fine mapping and replication in an independent sample of 454 French/Canadian CD trios supported the authenticity of the NELL1 association. Further confirmation in a large German ulcerative colitis (UC) sample indicated that NELL1 is a ubiquitous IBD susceptibility locus (combined p<10⁻⁶; OR = 1.66, 95% CI: 1.30–2.11). The novel 5p13.1 locus was also replicated in the French/Canadian sample and in an independent UK CD patient panel (453 cases, 521 controls, combined p<10⁻⁶ for SNP rs1992660). Several associations were replicated in at least one independent sample, point to an involvement of ITGB6 (upstream), GRM8 (downstream), OR5V1 (downstream), PPP3R2 (downstream), NM_152575 (upstream) and HNF4G (intron).

Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility

A genome-wide association scan in individuals with Crohn's disease by the Wellcome Trust Case Control Consortium detected strong association at four novel loci. We tested 37 SNPs from these and other loci for association in an independent case-control sample. We obtained replication for the autophagy-inducing IRGM gene on chromosome 5q33.1 (replication P = 6.6 x 10(-4), combined P = 2.1 x 10(-10)) and for nine other loci, including NKX2-3, PTPN2 and gene deserts on chromosomes 1q and 5p13.

A nonsynonymous SNP in ATG16L1 predisposes to ileal Crohn's disease and is independent of CARD15 and IBD5

BACKGROUND & AIMS

A genome-wide association scan of nonsynonymous DNA polymorphisms identified association of a threonine-to-alanine substitution (T300A) in the autophagy-related 16-like gene ATG16L1 with Crohn's disease. We investigated this association in independent U.K. cohorts of Crohn's disease and ulcerative colitis.

METHODS

The T300A variant (rs2241880) was genotyped in an independent sample of 727 Crohn's disease and 877 ulcerative colitis cases, and in 579 controls. We then performed an extension analysis combining these data with the U.K. data from the initial study to give a total of 1236 U.K. Crohn's disease cases and 1235 controls to estimate disease risk and test for interaction with the CARD15 and IBD5 risk loci and for association with disease subtypes.

RESULTS

The association of T300A was replicated in the independent sample of 727 Crohn's disease cases (P = .001), and was strongly associated in the extended analysis of 1236 Crohn's cases (P = 2.4 x 10(-6)). The 300A/A genotype conferred a 1.65-fold risk of Crohn's disease, with a 2.2-fold risk of ileal disease. Analysis of the interaction of ATG16L1 with CARD15 and IBD5 indicated that all 3 loci contribute independently to disease risk. Homozygosity for the risk allele at all 3 loci conferred a combined risk of 20.4 (95% confidence interval: 8.71, 47.7) for Crohn's disease. The ATG16L1 risk genotype showed a modest but significant association with ulcerative colitis (P = .026).

CONCLUSIONS

The association of ATG16L1 with Crohn's disease and possibly with ulcerative colitis supports a role for autophagy in the pathogenesis of inflammatory bowel disease.

IL23R variation determines susceptibility but not disease phenotype in inflammatory bowel disease

BACKGROUND & AIMS

Identification of inflammatory bowel disease (IBD) susceptibility genes is key to understanding pathogenic mechanisms. Recently, the North American IBD Genetics Consortium provided compelling evidence for an association between ileal Crohn's disease (CD) and the IL23R gene using genome-wide association scanning. External replication is a priority, both to confirm this finding in other populations and to validate this new technique. We tested for association between IL23R and IBD in a large independent UK panel to determine the size of the effect and explore subphenotype correlation and interaction with CARD15.

METHODS

Eight single nucleotide polymorphism markers in IL23R tested in the North American study were genotyped in 1902 cases of Crohn's disease (CD), 975 cases of ulcerative colitis (UC), and 1345 controls using MassARRAY. Data were analyzed using chi(2) statistics, and subgroup association was sought.

RESULTS

A highly significant association with CD was observed, with the strongest signal at coding variant Arg381Gln (allele frequency, 2.5% in CD vs 6.2% in controls [P = 1.1 x 10(-12)]; odds ratio, 0.38; 95% confidence interval, 0.29-0.50). A weaker effect was seen in UC (allele frequency, 4.6%; odds ratio, 0.73; 95% confidence interval, 0.55-0.96). Analysis accounting for Arg381Gln suggested that other loci within IL23R also influence IBD susceptibility. Within CD, there were no subphenotype associations or evidence of interaction with CARD15.

CONCLUSIONS

This study shows an association between IL23R and all subphenotypes of CD with a smaller effect on UC. This extends the findings of the North American study, providing clear evidence that genome-wide association scanning can successfully identify true complex disease genes.

Investigation of association of the DLG5 gene with phenotypes of inflammatory bowel disease in the British population

BACKGROUND AND AIMS

Mutations in the DLG5 gene are associated with an increased risk of inflammatory bowel disease (IBD) in some European populations. Initial investigation of a British IBD population showed evidence of association of one of three DLG5 variants, R30Q, in a family-based collection but not in a case-control cohort. We have now examined the association of the R30Q polymorphism in a large cohort of British IBD cases, tested for interaction between the DLG5 and CARD15 genes and assessed possible association of DLG5 with clinical features of Crohn's disease (CD) and ulcerative colitis (UC).

MATERIALS AND METHODS

DLG5 R30Q and the CARD15 polymorphisms, Arg702Trp, Gly908Arg and Leu1007fs were genotyped in 1,148 IBD cases and 749 controls. DLG5 R30Q was also analysed in cases stratified by CARD15 genotype, disease subtype and smoking history.

RESULTS/FINDINGS

No significant difference in frequencies of the R30Q variant was observed between IBD cases (9.9%) and controls (10.1%) or in cases analysed separately as CD and UC. There was also no significant difference in the frequency of R30Q between CD cases carrying risk-associated CARD15 alleles and those that did not. The frequency of R30Q was higher in CD cases with ileal disease than cases without (p=0.042) and higher in CD cases who had smoked than in nonsmokers (p=0.009).

INTERPRETATION/CONCLUSION

The R30Q variant in the DLG5 gene does not appear to be associated with an overall increase in the risk of disease in a British IBD cohort, but differences in its frequency in subgroups of CD patients warrant further investigation.

A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P <or= 0.01 with an allele-based disease association test in 380 independent Crohn disease trios, 498 Crohn disease singleton cases and 1,032 controls. Disease association of rs2241880 in the autophagy-related 16-like 1 gene (ATG16L1) was replicated in these samples (P = 4.0 x 10(-8)) and confirmed in a UK case-control sample (P = 0.0004). By haplotype and regression analysis, we found that marker rs2241880, a coding SNP (T300A), carries virtually all the disease risk exerted by the ATG16L1 locus. The ATG16L1 gene encodes a protein in the autophagosome pathway that processes intracellular bacteria. We found a statistically significant interaction with respect to Crohn disease risk between rs2241880 and the established CARD15 susceptibility variants (P = 0.039). Together with the lack of association between rs2241880 and ulcerative colitis (P > 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.

Genetic variation in myosin IXB is associated with ulcerative colitis

BACKGROUND & AIMS

Common germline genetic variation in the 3' region of myosin IXB (MYO9B) has been associated recently with susceptibility to celiac disease, with a hypothesis that MYO9B variants might influence intestinal permeability. These findings suggested the current study investigating a possible further role for MYO9B variation in inflammatory bowel disease.

METHODS

Eight single-nucleotide polymorphisms (SNPs) were selected to tag common haplotypes from the 35-kb 3' region of MYO9B. These included the strongest celiac disease-associated variants reported in a Dutch cohort. These SNPs were studied in 3 independently collected and genotyped case-control cohorts of European descent (UK, Dutch, and Canadian/Italian), comprising in total 2717 inflammatory bowel disease patients (1197 with Crohn's disease, 1520 with ulcerative colitis) and 4440 controls.

RESULTS

Common variation in MYO9B was associated with susceptibility to inflammatory bowel disease in all 3 cohorts examined (most associated SNP, rs1545620; meta-analysis P = 1.9 x 10(-6); odds ratio, 1.2), with the same alleles showing association as reported for celiac disease.

CONCLUSIONS

MYO9B genetic variants predispose to inflammatory bowel disease. Interestingly, rs1545620 is a nonsynonymous variant leading to an amino acid change (Ala1011Ser) in the third calmodulin binding IQ domain of MYO9B. Unlike previous variants (in other genes) reported to predispose to inflammatory bowel disease, the association at MYO9B was considerably stronger with ulcerative colitis, although weaker association with Crohn's disease also was observed. These data imply shared causal mechanisms underlying intestinal inflammatory diseases.

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.

Associations of allelic variants of the multidrug resistance gene (ABCB1 or MDR1) and inflammatory bowel disease and their effects on disease behavior: a case-control and meta-analysis study

BACKGROUND

Allelic variants of the ATP-binding cassette, subfamily B member 1 (ABCB1), also known as the multidrug resistance gene (MDR1) that encodes the membrane-bound efflux transporter P-glycoprotein 170 (PGP-170), have been associated with inflammatory bowel disease but with conflicting results.

METHODS

The present study examined the association of ABCB1 C3435T and G2677T/A in a large British case-control cohort of 828 Crohn's disease, 580 ulcerative colitis (UC) cases, and 285 healthy controls. The effect of these variants was further examined with respect to phenotypic and epidemiological characteristics. A meta-analysis was carried out of our results and those from 8 previously published association studies of the C3435T variant in inflammatory bowel disease.

RESULTS

The 2677T allele was significantly increased in British UC cases compared with controls (45.2% vs. 39.6%; P = 0.034). In particular, the TT genotype was significantly associated with severe UC (odds ratio [OR] 1.90; 95% CI 1.01-3.55) and the use of steroids in UC (OR 1.77; 95% CI 1.08-2.88). No significant association was seen with C3435T and UC, Crohn's disease, or any clinical subgroup. A meta-analysis of 9 association studies of C3435T showed a significant association of the 3435T allele with UC (OR 1.12; 95% CI 1.02-1.23; P = 0.013) but not with CD.

CONCLUSIONS

These results indicate that ABCB1 sequence variants are associated with a small increase in the risk of developing UC and may influence disease behavior.

Mutation, selection, and evolution of the Crohn disease susceptibility geneCARD15

Three common mutations in the CARD15 (NOD2) gene are known to be associated with susceptibility to Crohn disease (CD), and genetic data suggest a gene dosage model with an increased risk of 2-4-fold in heterozygotes and 20-40-fold in homozygotes. However, the discovery of numerous rare variants of CARD15 indicates that some heterozygotes for the common mutations have a rare mutation on the other CARD15 allele, which would support a recessive model for CD. We addressed this issue by screening CARD15 for mutations in 100 CD patients who were heterozygous for one of the three common mutations. We also developed a strategy for evaluating potential disease susceptibility alleles (DSAs) that involves assessing the degree of evolutionary conservation of involved residues, predicted effects on protein structure and function, and genotyping in a large sample of cases and controls. The evolutionary analysis was aided by sequencing the entire coding region of CARD15 in three primates (chimp, gibbon, and tamarin) and aligning the human sequence with these and orthologs from other species. We found that 11 of the 100 CD patients screened had a second potential pathogenic mutation within the exonic and periexonic sequences examined. Assuming that there are no additional pathogenic mutations in noncoding regions, our study suggests that most carriers of the common DSAs are true heterozygotes, and supports previous evidence for a gene dosage model. Four novel nonsynonymous mutations were detected, one of which would produce premature termination of translation c.2686C>T (p.Arg896X). Two potential DSAs--c.2107C>T (p.Arg703Cys) and g.2238T>A (c.74-7T>A)--were significantly associated with CD in the case control sample. Analysis of the evolution of CARD15 revealed strong conservation of the encoded protein, with identity to the human sequence ranging from 99.1% in the chimp to 44.5% in fugu. Higher primates possess an open reading frame (ORF) upstream of the putative initiation site in other species that encodes a further 27 N-terminal amino acids, while four regions of high conservation are observed outside of the known domains of CARD15, indicative of additional residues of functional importance. The strategy developed here may have general application to the assessment of mutation pathogenicity and genetic models in other complex disorders.