EPHA2 polymorphisms and age-related cataract in India

Objective

We investigated whether previously reported single nucleotide polymorphisms (SNPs) of EPHA2 in European studies are associated with cataract in India.

Methods

We carried out a population-based genetic association study. We enumerated randomly sampled villages in two areas of north and south India to identify people aged 40 and over. Participants attended a clinical examination including lens photography and provided a blood sample for genotyping. Lens images were graded by the Lens Opacification Classification System (LOCS III). Cataract was defined as a LOCS III grade of nuclear ≥4, cortical ≥3, posterior sub-capsular (PSC) ≥2, or dense opacities or aphakia/pseudophakia in either eye. We genotyped SNPs rs3754334, rs7543472 and rs11260867 on genomic DNA extracted from peripheral blood leukocytes using TaqMan assays in an ABI 7900 real-time PCR. We used logistic regression with robust standard errors to examine the association between cataract and the EPHA2 SNPs, adjusting for age, sex and location.

Results

7418 participants had data on at least one of the SNPs investigated. Genotype frequencies of controls were in Hardy-Weinberg Equilibrium (p>0.05). There was no association of rs3754334 with cataract or type of cataract. Minor allele homozygous genotypes of rs7543472 and rs11260867 compared to the major homozygote genotype were associated with cortical cataract, Odds ratio (OR) = 1.8, 95% Confidence Interval (CI) (1.1, 3.1) p = 0.03 and 2.9 (1.2, 7.1) p = 0.01 respectively, and with PSC cataract, OR = 1.5 (1.1, 2.2) p = 0.02 and 1.8 (0.9, 3.6) p = 0.07 respectively. There was no consistent association of SNPs with nuclear cataract or a combined variable of any type of cataract including operated cataract.

Conclusions

Our results in the Indian population agree with previous studies of the association of EPHA2 variants with cortical cataracts. We report new findings for the association with PSC which is particularly prevalent in Indians.     

Effects of metabotropic glutamate receptor 3 genotype on phonetic mismatch negativity

Background

The genetic and molecular basis of glutamatergic dysfunction is one key to understand schizophrenia, with the identification of an intermediate phenotype being an essential step. Mismatch negativity (MMN) or its magnetic counterpart, magnetic mismatch field (MMF) is an index of preattentive change detection processes in the auditory cortex and is generated through glutamatergic neurotransmission. We have previously shown that MMN/MMF in response to phoneme change is markedly reduced in schizophrenia. Variations in metabotropic glutamate receptor (GRM3) may be associated with schizophrenia, and has been shown to affect cortical function. Here we investigated the effect of GRM3 genotypes on phonetic MMF in healthy men.

Methods

MMF in response to phoneme change was recorded using magnetoencephalography in 41 right-handed healthy Japanese men. Based on previous genetic association studies in schizophrenia, 4 candidate SNPs (rs6465084, rs2299225, rs1468412, rs274622) were genotyped.

Results

GRM3 rs274622 genotype variations significantly predicted MMF strengths (p = 0.009), with C carriers exhibiting significantly larger MMF strengths in both hemispheres compared to the TT subjects.

Conclusions

These results suggest that variations in GRM3 genotype modulate the auditory cortical response to phoneme change in humans. MMN/MMF, particularly those in response to speech sounds, may be a promising and sensitive intermediate phenotype for clarifying glutamatergic dysfunction in schizophrenia.     

Impact of genetic variation in SORCS1 on memory retention

Objective

We previously reported that genetic variants in SORCS1 increase the risk of AD, that over-expression of SorCS1 reduces γ-secretase activity and Aβ levels, and that SorCS1 suppression increases γ-secretase processing of APP and Aβ levels. We now explored the effect of variation in SORCS1 on memory.

Methods

We explored associations between SORCS1-SNPs and memory retention in the NIA-LOAD case control dataset (162 cases,670 controls) and a cohort of Caribbean Hispanics (549 cases,544 controls) using single marker and haplotype analyses.

Results

Three SNPs in intron 1, were associated with memory retention in the NIA-LOAD dataset or the Caribbean Hispanic dataset (rs10884402(A allele:β = −0.15,p = 0.008), rs7078098(C allele:β = 0.18,p = 0.007) and rs950809(C allele:β = 0.17,p = 0.008)) and all three SNPs were significant in a meta-analysis of both datasets (0.002ConclusionsVariation in intron 1 in SORCS1 is associated with memory changes in AD.     

Common genetic variation in the SERPINF1 locus determines overall adiposity, obesity-related insulin resistance, and circulating leptin levels

Objective

Pigment epithelium-derived factor (PEDF) belongs to the serpin family of peptidase inhibitors (serpin F1) and is among the most abundant glycoproteins secreted by adipocytes. In vitro and mouse in vivo data revealed PEDF as a candidate mediator of obesity-induced insulin resistance. Therefore, we assessed whether common genetic variation within the SERPINF1 locus contributes to adipose tissue-related prediabetic phenotypes in humans.

Subjects/Methods

A population of 1,974 White European individuals at increased risk for type 2 diabetes was characterized by an oral glucose tolerance test with glucose and insulin measurements (1,409 leptin measurements) and genotyped for five tagging SNPs covering 100% of common genetic variation (minor allele frequency ≥0.05) in the SERPINF1 locus. In addition, a subgroup of 486 subjects underwent a hyperinsulinaemic-euglycaemic clamp and a subgroup of 340 magnetic resonance imaging (MRI) and spectroscopy (MRS).

Results

After adjustment for gender and age and Bonferroni correction for the number of SNPs tested, SNP rs12603825 revealed significant association with MRI-derived total adipose tissue mass (p = 0.0094) and fasting leptin concentrations (p = 0.0035) as well as nominal associations with bioelectrical impedance-derived percentage of body fat (p = 0.0182) and clamp-derived insulin sensitivity (p = 0.0251). The association with insulin sensitivity was completely abolished by additional adjustment for body fat (p = 0.8). Moreover, the fat mass-increasing allele of SNP rs12603825 was significantly associated with elevated fasting PEDF concentrations (p = 0.0436), and the PEDF levels were robustly and positively associated with all body fat parameters measured and with fasting leptin concentrations (p<0.0001, all).

Conclusion

In humans at increased risk for type 2 diabetes, a functional common genetic variant in the gene locus encoding PEDF contributes to overall body adiposity, obesity-related insulin resistance, and circulating leptin levels.     

Variants in the FFAR1 gene are associated with beta cell function

Background

The FFAR1 receptor is expressed mainly in pancreatic beta cells and is activated by medium to long chain free fatty acids (FFAs), as well as by thiazolidinediones, resulting in elevated Ca²⁺ concentrations and promotion of insulin secretion. These properties suggest that FFAR1 could be a mediator of lipotoxicity and a potential candidate gene for Type 2 diabetes (T2D). We therefore investigated whether variations at the FFAR1 locus are associated with T2D and beta cell function.

Methodology/Principal Findings

We re-sequenced the FFAR1 region in 96 subjects (48 healthy and 48 T2D individuals) and found 13 single nucleotide polymorphisms (SNPs) 8 of which were not previously described. Two SNPs located in the upstream region of the FFAR1 gene (rs1978013 and rs1978014) were chosen and genotyped in 1929 patients with T2D and 1405 healthy control subjects. We observed an association of rs1978013 and rs1978014 with insulinogenic index in males (p = 0.024) and females (p = 0.032), respectively. After Bonferroni corrections, no association with T2D was found in the case-control material, however a haplotype consisting of the T-G alleles conferred protection against T2D (p = 0.0010).

Conclusions/Significance

Variation in the FFAR1 gene may contribute to impaired beta cell function in T2D.     

Variance of the SGK1 gene is associated with insulin secretion in different European populations: results from the TUEF, EUGENE2, and METSIM studies

Hypothesis

Serum- and Glucocorticoid-inducible Kinase 1 (SGK1) is involved in the regulation of insulin secretion and may represent a candidate gene for the development of type 2 diabetes mellitus in humans.

Methods

Three independent European populations were analyzed for the association of SGK1 gene (SGK) variations and insulin secretion traits. The German TUEF project provided the screening population (N = 725), and four tagging SNPs (rs1763527, rs1743966, rs1057293, rs9402571) were investigated. EUGENE2 (N = 827) served as a replication cohort for the detected associations. Finally, the detected associations were validated in the METSIM study, providing 3798 non-diabetic and 659 diabetic (type 2) individuals.

Results

Carriers of the minor G allele in rs9402571 had significantly higher C-peptide levels in the 2 h OGTT (+10.8%, p = 0.04; dominant model) and higher AUC_C-Peptide/AUC_Glc ratios (+7.5%, p = 0.04) compared to homozygous wild type TT carriers in the screening population. As interaction analysis for BMI×rs9402571 was significant (p = 0.04) for the endpoint insulin secretion, we stratified the TUEF cohort for BMI, using a cut off point of BMI = 25. The effect on insulin secretion only remained significant in lean TUEF participants (BMI≤25). This finding was replicated in lean EUGENE2 rs9402571 minor allele carriers, who had a significantly higher AUC_Ins/AUC_Glc (TT: 226±7, XG: 246±9; p = 0.019). Accordingly, the METSIM trial revealed a lower prevalence of type 2 diabetes (OR: 0.85; 95%CI: 0.71–1.01; p = 0.065, dominant model) in rs9402571 minor allele carriers.

Conclusions

The rs9402571 SGK genotype associates with increased insulin secretion in lean non-diabetic TUEF/EUGENE2 participants and with lower diabetes prevalence in METSIM. Our study in three independent European populations supports the conclusion that SGK variability affects diabetes risk.     

A common HLA-DPA1 variant is associated with hepatitis B virus infection but fails to distinguish active from inactive Caucasian carriers

Background and Aims

Chronic infection with the hepatitis B virus (HBV) is a major health issue worldwide. Recently, single nucleotide polymorphisms (SNPs) within the human leukocyte antigen (HLA)-DP locus were identified to be associated with HBV infection in Asian populations. Most significant associations were observed for the A alleles of HLA-DPA1 rs3077 and HLA-DPB1 rs9277535, which conferred a decreased risk for HBV infection. We assessed the implications of these variants for HBV infection in Caucasians.

Methods

Two HLA-DP gene variants (rs3077 and rs9277535) were analyzed for associations with persistent HBV infection and with different clinical outcomes, i.e., inactive HBsAg carrier status versus progressive chronic HBV (CHB) infection in Caucasian patients (n = 201) and HBsAg negative controls (n = 235).

Results

The HLA-DPA1 rs3077 C allele was significantly associated with HBV infection (odds ratio, OR = 5.1, 95% confidence interval, CI: 1.9–13.7; p = 0.00093). However, no significant association was seen for rs3077 with progressive CHB infection versus inactive HBsAg carrier status (OR = 2.7, 95% CI: 0.6–11.1; p = 0.31). In contrast, HLA-DPB1 rs9277535 was not associated with HBV infection in Caucasians (OR = 0.8, 95% CI: 0.4–1.9; p = 1).

Conclusions

A highly significant association of HLA-DPA1 rs3077 with HBV infection was observed in Caucasians. However, as a differentiation between different clinical courses of HBV infection was not possible, knowledge of the HLA-DPA1 genotype cannot be translated into personalized anti-HBV therapy approaches.     

The NOD2 single nucleotide polymorphisms rs2066843 and rs2076756 are novel and common Crohn's disease susceptibility gene variants

Background

The aims were to analyze two novel NOD2 variants (rs2066843 and rs2076756) in a large cohort of patients with inflammatory bowel disease and to elucidate phenotypic consequences.

Methodology/Principal Findings

Genomic DNA from 2700 Caucasians including 812 patients with Crohn''s disease (CD), 442 patients with ulcerative colitis (UC), and 1446 healthy controls was analyzed for the NOD2 SNPs rs2066843 and rs2076756 and the three main CD-associated NOD2 variants p.Arg702Trp (rs2066844), p.Gly908Arg (rs2066847), and p.Leu1007fsX1008 (rs2066847). Haplotype and genotype-phenotype analyses were performed. The SNPs rs2066843 (p = 3.01×10⁻⁵, OR 1.48, [95% CI 1.23-1.78]) and rs2076756 (p = 4.01×10⁻⁶; OR 1.54, [95% CI 1.28-1.86]) were significantly associated with CD but not with UC susceptibility. Haplotype analysis revealed a number of significant associations with CD susceptibility with omnibus p values <10⁻¹⁰. The SNPs rs2066843 and rs2076756 were in linkage disequilibrium with each other and with the three main CD-associated NOD2 mutations (D''>0.9). However, in CD, SNPs rs2066843 and rs2076756 were more frequently observed than the other three common NOD2 mutations (minor allele frequencies for rs2066843 and rs2076756: 0.390 and 0.380, respectively). In CD patients homozygous for these novel NOD2 variants, genotype-phenotype analysis revealed higher rates of a penetrating phenotype (rs2076756: p = 0.015) and fistulas (rs2076756: p = 0.015) and significant associations with CD-related surgery (rs2076756: p = 0.003; rs2066843: p = 0.015). However, in multivariate analysis only disease localization (p<2×10⁻¹⁶) and behaviour (p = 0.02) were significantly associated with the need for surgery.

Conclusion/Significance

The NOD2 variants rs2066843 and rs2076756 are novel and common CD susceptibility gene variants.     

Genetic determinants of cardiovascular events among women with migraine: a genome-wide association study

Background

Migraine is associated with an increased risk for cardiovascular disease (CVD). Both migraine and CVD are highly heritable. However, the genetic liability for CVD among migraineurs is unclear.

Methods

We performed a genome-wide association study for incident CVD events during 12 years of follow-up among 5,122 migraineurs participating in the population-based Women''s Genome Health Study. Migraine was self-reported and CVD events were confirmed after medical records review. We calculated odds ratios (OR) and 95% confidence intervals (CI) and considered a genome-wide p-value <5×10⁻⁸ as significant.

Results

Among the 5,122 women with migraine 164 incident CVD events occurred during follow-up. No SNP was associated with major CVD, ischemic stroke, myocardial infarction, or CVD death at the genome-wide level; however, five SNPs showed association with p<5×10⁻⁶. Among migraineurs with aura rs7698623 in MEPE (OR = 6.37; 95% CI 3.15–12.90; p = 2.7×10⁻⁷) and rs4975709 in IRX4 (OR = 5.06; 95% CI 2.66–9.62; p = 7.7×10⁻⁷) appeared to be associated with ischemic stroke, rs2143678 located close to MDF1 with major CVD (OR = 3.05; 95% CI 1.98–4.69; p = 4.3×10⁻⁷), and the intergenic rs1406961 with CVD death (OR = 12.33; 95% CI 4.62–32.87; p = 5.2×10⁻⁷). Further, rs1047964 in BACE1 appeared to be associated with CVD death among women with any migraine (OR = 4.67; 95% CI 2.53–8.62; p = 8.0×10⁻⁷).

Conclusion

Our results provide some suggestion for an association of five SNPs with CVD events among women with migraine; none of the results was genome-wide significant. Four associations appeared among migraineurs with aura, two of those with ischemic stroke. Although our population is among the largest with migraine and incident CVD information, these results must be treated with caution, given the limited number of CVD events among women with migraine and the low minor allele frequencies for three of the SNPs. Our results await independent replication and should be considered hypothesis generating for future research.     

An exploratory study on CLU, CR1 and PICALM and Parkinson disease

Background

Recent GWAS and subsequent confirmation studies reported several single-nucleotide polymorphisms (SNPs) at the CLU, CR1 and PICALM loci in association with late-onset Alzheimer''s disease (AD). Parkinson disease (PD) shares several clinical and pathologic characteristics with AD; we therefore explored whether these SNPs were also associated with PD risk.

Methodology/Principal Findings

791 non-Hispanic Whites cases and 1,580 matched controls were included in the study. Odds ratios (OR) and 95% confidence intervals (CI) were obtained from logistic regression models. rs11136000 at the CLU locus was associated with PD risk under the recessive model (comparing TT versus CC+CT: OR = 0.71, 95% CI: 0.55-0.92, p = 0.008) after adjusting for year of birth, gender, smoking, and caffeine intake. Further adjustment for family history of PD and ApoE ε4 status did not change the result. In addition, we did not find evidence for effect modification by ApoE or known PD risk factors. The association, however, appeared to be stronger for PD with dementia (OR = 0.49, 95% CI: 0.27-0.91) than for PD without dementia (OR = 0.81, 95% CI: 0.61-1.06). The two other SNPs, rs6656401 from CR1, and rs3851179 from PICALM region were not associated with PD (p>0.05).

Conclusion

Our exploratory analysis suggests an association of CLU with PD. This exploratory finding and the role of dementia in explaining this finding needs further investigation.     

A Single Nucleotide Polymorphism within DUSP9 Is Associated with Susceptibility to Type 2 Diabetes in a Japanese Population

Aims

The DUSP9 locus on chromosome X was identified as a susceptibility locus for type 2 diabetes in a meta-analysis of European genome-wide association studies (GWAS), and GWAS in South Asian populations identified 6 additional single nucleotide polymorphism (SNP) loci for type 2 diabetes. However, the association of these loci with type 2 diabetes have not been examined in the Japanese. We performed a replication study to investigate the association of these 7 susceptibility loci with type 2 diabetes in the Japanese population.

Methods

We genotyped 11,319 Japanese participants (8,318 with type 2 diabetes and 3,001 controls) for each of the 7 SNPs–rs5945326 near DUSP9, rs3923113 near GRB14, rs16861329 in ST6GAL1, rs1802295 in VPS26A, rs7178572 in HMG20A, rs2028299 near AP3S2, and rs4812829 in HNF4A–and examined the association of each of these 7 SNPs with type 2 diabetes by using logistic regression analysis.

Results

All SNPs had the same direction of effect (odds ratio [OR]>1.0) as in the original reports. One SNP, rs5945326 near DUSP9, was significantly associated with type 2 diabetes at a genome-wide significance level (p = 2.21×10⁻⁸; OR 1.39, 95% confidence interval [CI]: 1.24−1.56). The 6 SNPs derived from South Asian GWAS were not significantly associated with type 2 diabetes in the Japanese population by themselves (p≥0.007). However, a genetic risk score constructed from 6 South Asian GWAS derived SNPs was significantly associated with Japanese type 2 diabetes (p = 8.69×10⁻⁴, OR  = 1.06. 95% CI; 1.03−1.10).

Conclusions/interpretation

These results indicate that the DUSP9 locus is a common susceptibility locus for type 2 diabetes across different ethnicities, and 6 loci identified in South Asian GWAS also have significant effect on susceptibility to Japanese type 2 diabetes.     

Possible Associations of NTRK2 Polymorphisms with Antidepressant Treatment Outcome: Findings from an Extended Tag SNP Approach

Background

Data from clinical studies and results from animal models suggest an involvement of the neurotrophin system in the pathology of depression and antidepressant treatment response. Genetic variations within the genes coding for the brain-derived neurotrophic factor (BDNF) and its key receptor Trkb (NTRK2) may therefore influence the response to antidepressant treatment.

Methods

We performed a single and multi-marker association study with antidepressant treatment outcome in 398 depressed Caucasian inpatients participating in the Munich Antidepressant Response Signature (MARS) project. Two Caucasian replication samples (N = 249 and N = 247) were investigated, resulting in a total number of 894 patients. 18 tagging SNPs in the BDNF gene region and 64 tagging SNPs in the NTRK2 gene region were genotyped in the discovery sample; 16 nominally associated SNPs were tested in two replication samples.

Results

In the discovery analysis, 7 BDNF SNPs and 9 NTRK2 SNPs were nominally associated with treatment response. Three NTRK2 SNPs (rs10868223, rs1659412 and rs11140778) also showed associations in at least one replication sample and in the combined sample with the same direction of effects (P_corr = .018, P_corr = .015 and P_corr = .004, respectively). We observed an across-gene BDNF-NTRK2 SNP interaction for rs4923468 and rs1387926. No robust interaction of associated SNPs was found in an analysis of BDNF serum protein levels as a predictor for treatment outcome in a subset of 93 patients.

Conclusions/Limitations

Although not all associations in the discovery analysis could be unambiguously replicated, the findings of the present study identified single nucleotide variations in the BDNF and NTRK2 genes that might be involved in antidepressant treatment outcome and that have not been previously reported in this context. These new variants need further validation in future association studies.     

Association of new loci identified in European genome-wide association studies with susceptibility to type 2 diabetes in the Japanese

Background

Several novel susceptibility loci for type 2 diabetes have been identified through genome-wide association studies (GWAS) for type 2 diabetes or quantitative traits related to glucose metabolism in European populations. To investigate the association of the 13 new European GWAS-derived susceptibility loci with type 2 diabetes in the Japanese population, we conducted a replication study using 3 independent Japanese case-control studies.

Methodology/Principal Findings

We examined the association of single nucleotide polymorphisms (SNPs) within 13 loci (MTNR1B, GCK, IRS1, PROX1, BCL11A, ZBED3, KLF14, TP53INP1, KCNQ1, CENTD2, HMGA2, ZFAND6 and PRC1) with type 2 diabetes using 4,964 participants (2,839 cases and 2,125 controls) from 3 independent Japanese samples. The association of each SNP with type 2 diabetes was analyzed by logistic regression analysis. Further, we performed combined meta-analyses for the 3 studies and previously performed Japanese GWAS data (4,470 cases vs. 3,071 controls). The meta-analysis revealed that rs2943641 in the IRS1 locus was significantly associated with type 2 diabetes, (P = 0.0034, OR = 1.15 95% confidence interval; 1.05–1.26) and 3 SNPs, rs10930963 in the MTNR1B locus, rs972283 in the KLF14 locus, and rs231362 in the KCNQ1 locus, had nominal association with type 2 diabetes in the present Japanese samples (P<0.05).

Conclusions

These results indicate that IRS1 locus may be common locus for type 2 diabetes across different ethnicities.     

A common variant in the telomerase RNA component is associated with short telomere length

Background

Telomeres shorten as cells divide. This shortening is compensated by the enzyme telomerase. We evaluated the effect of common variants in the telomerase RNA component (TERC) gene on telomere length (TL) in the population-based Health Aging and Body Composition (Health ABC) Study and in two replication samples (the TwinsUK Study and the Amish Family Osteoporosis Study, AFOS).

Methodology

Five variants were identified in the TERC region by sequence analysis and only one SNP was common (rs2293607, G/A). The frequency of the G allele was 0.26 and 0.07 in white and black, respectively. Testing for association between TL and rs2293607 was performed using linear regression models or variance component analysis conditioning on relatedness among subjects.

Results

The adjusted mean TL was significantly shorter in 665 white carriers of the G allele compared to 887 non-carriers from the Health ABC Study (4.69±0.05 kbp vs. 4.86±0.04 kbp, measured by quantitative PCR, p = 0.005). This association was replicated in another white sample from the TwinsUK Study (6.90±0.03 kbp in 301 carriers compared to 7.06±0.03 kbp in 395 non-carriers, measured by Southern blots, p = 0.009). A similar pattern of association was observed in whites from the family-based AFOS and blacks from the Health ABC cohort, although not statistically significant, possibly due to the lower allele frequency in these populations. Combined analysis using 2,953 white subjects from 3 studies showed a significant association between TL and rs2293607 (β = −0.19±0.04 kbp, p = 0.001).

Conclusion

Our study shows a significant association between a common variant in TERC and TL in humans, suggesting that TERC may play a role in telomere homeostasis.     

The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease

Background

Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients.

Methodology/Principal Findings

Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn''s disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction.

Conclusions/Significance

Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.