Association between the ENPP1 K121Q Polymorphism and Risk of Diabetic Kidney Disease: A Systematic Review and Meta-Analysis

The potential association between the K121Q (A/C, rs1044498) polymorphism in the ectonucleotide pyrophosphatase/phosphodiesterase (ENPP1) gene and risk of diabetic kidney disease (DKD) has been investigated. Nevertheless, the effect of this variant on DKD risk is still under debate, and conflicting results have been reported. To this date, no meta-analysis has evaluated the association of the K121Q polymorphism with DKD. This paper describes the first meta-analysis conducted to evaluate whether the ENPP1K121Q polymorphism is associated with DKD. A literature search was conducted to identify all case-control or cross-sectional studies that evaluated associations between the ENPP1K121Q polymorphism and DKD. Pooled odds ratios (OR) and 95% confidence intervals (95% CI) were calculated for allele contrast, additive, dominant and recessive inheritance models. Seven studies were eligible for inclusion in the meta-analysis, providing data on 3571 type 1 or type 2 diabetic patients (1606 cases with DKD and 1965 diabetic controls without this complication). No significant heterogeneity was observed among the studies included in the meta-analysis when assuming different inheritance models (I^² < 50% or P > 0.10 for the entire sample and after stratification by ethnicity). Meta-analysis results revealed significant associations between the K121Q polymorphism and risk of DKD in Asians and Europeans when assuming the different inheritance models analyzed. The most powerful association was observed for the additive model (OR = 1.74, 95% CI 1.27-2.38 for the total sample). In conclusion, the present meta-analysis detected a significant association between the ENPP1K121Q polymorphism and increased susceptibility of DKD in European and Asian populations.     

The Q121/Q121 genotype of ENPP1/PC-1 is associated with lower BMI in non-diabetic whites

This study investigated the role of the ENPP1/PC-1 gene K121Q polymorphism in predicting BMI (kg/m2) in non-diabetic individuals. Three independent samples (n = 631, n = 304, and n = 505) of adult whites were analyzed. Selection criteria were fasting plasma glucose level <126 mg/dL, absence of severe obesity (BMI > or =40 kg/m2), and lack of treatment known to modulate BMI. In Sample 1, BMI values were different in individuals carrying the K121/K121 (KK), K121/Q121 (KQ), and Q121/Q121 (QQ) genotypes (25.5 +/- 4.3, 25.3 +/- 4.1, and 22.8 +/- 2.5 kg/m2, respectively (adjusted p = 0.022); BMI values in Samples 2 and 3 also tended to be different, although the differences, after adjustment for age and sex, did not reach statistical significance. When data were pooled, BMI values were 25.8 +/- 4.4, 25.6 +/- 4.4, and 23.6 +/- 3.3 kg/m2 in KK, KQ, and QQ individuals (adjusted p = 0.029). According to a recessive model, QQ individuals had lower BMI values than KK and KQ individuals combined (23.6 +/- 3.3 kg/m2 vs. 25.7 +/- 4.4 kg/m2; adjusted p = 0.008). These data suggest that the QQ genotype of the ENPP1/PC-1 gene is associated with lower BMI. If similar results are confirmed in prospective studies, the K121Q polymorphism may help identify people at risk for obesity.     

The ENPP1 K121Q polymorphism is not associated with type 2 diabetes and related metabolic traits in an Iranian population

The K121Q polymorphism of the ectoenzyme nucleotide pyrophosphate phosphodiesterase 1 (ENPP1) gene has been variably associated with insulin resistance and type 2 diabetes (T2D) in several populations. However, this association has not been studied in Iranian subjects and we hypothesized that the K121Q variant might be associated with T2D and related metabolic traits in this population. The K121Q genotypes were determined by PCR-restriction fragment length polymorphism in 377 normoglycemic controls and 155 T2D patients. T2D patients had significantly higher values for systolic and diastolic blood pressure, BMI, glucose, cholesterol, triglyceride, LDL, apoB, insulin, and HOMA-IR, and lower levels of HDL than the normoglycemic subjects. The frequency of the Q allele did not differ between T2D and normoglycemic subjects (OR 0.96, 95% CI 0.90-2.00, P?=?0.70). The Q allele frequency was 16.5% in T2D and 15.2% in normoglycemic subjects. The ENPP1 genotype (KQ?+?QQ) was not associated with the systolic and diastolic blood pressure, glucose, triglyceride, cholesterol, LDL-C and HDL-C, apo B, BMI, HOMA-IR, and insulin levels in both normoglycemic and T2D groups. Our results suggest that the ENPP1 121Q allele might not be associated with T2D and related metabolic traits among Iranian subjects.     

The Q121 variant of ENPP1 may protect from childhood overweight/obesity in the Italian population

Ectonucleotide Pyrophosphatase Phosphodiesterase 1 (ENPP1) downregulates insulin signaling by inhibiting the insulin receptor's tyrosine-kinase. K121Q and other ENPP1 single-nucleotide polymorphisms (SNPs), IVS20delT-11 and A/G+1044TGA, have been previously associated with obesity in French children, and the risk haplotype QdelTG has also been associated with this condition in both French and German children. Our aim was to perform a case-control replication study in order to assess the possible association of childhood obesity and overweight with the above-mentioned ENPP1 SNPs, and with the QdelTG haplotype, in the Italian population. A total of 865 healthy Italian children were studied: 453 normal-weight, 243 overweight and 169 obese subjects. Genotyping was performed by Taq-Man or Light-Cycler Technology. The Q variant of K121Q showed a negative association with overweight-obesity under both additive (odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.57-0.97, P = 0.030) and recessive (OR = 0.32, 95% CI = 0.10-0.97, P = 0.035) modes of inheritance. The Z-score of BMI showed a significant decreasing trend from children K/K homozygous to K/Q heterozygous, and to Q/Q homozygous (0.45 vs. 0.28 vs. -0.19; P = 0.009), according to the additive model. The two other SNPs and the QdelTG haplotype did not exhibit any association with overweight/obesity. This is the first child-based study showing a protective role of the 121Q variant of ENPP1 against overweight/obesity.     

Further evidence for the role of ENPP1 in obesity: association with morbid obesity in Finns

The aim of this study was to investigate a series of single-nucleotide polymorphisms (SNPs) in the genes MC2R, MC3R, MC4R, MC5R, POMC, and ENPP1 for association with obesity. Twenty-five SNPs (2-7 SNPs/gene) were genotyped in 246 Finns with extreme obesity (BMI > or = 40 kg/m2) and in 481 lean subjects (BMI 20-25 kg/m2). Of the obese subjects, 23% had concomitant type 2 diabetes. SNPs and SNP haplotypes were tested for association with obesity and type 2 diabetes. Allele frequencies differed between obese and lean subjects for two SNPs in the ENPP1 gene, rs1800949 (P = 0.006) and rs943003 (P = 0.0009). These SNPs are part of a haplotype (rs1800949 C-rs943003 A), which was observed more frequently in lean subjects compared to obese subjects (P = 0.0007). Weaker associations were detected between the SNPs rs1541276 in the MC5R, rs1926065 in the MC3R genes and obesity (P = 0.04 and P = 0.03, respectively), and between SNPs rs2236700 in the MC5R, rs2118404 in the POMC, rs943003 in the ENPP1 genes and type 2 diabetes (P = 0.03, P = 0.02 and P = 0.02, respectively); these associations did not, however, remain significant after correction for multiple testing. In conclusion, a previously unexplored ENPP1 haplotype composed of SNPs rs1800949 and rs943003 showed suggestive evidence for association with adult-onset morbid obesity in Finns. In this study, we did not find association between the frequently studied ENPP1 K121Q variant, nor SNPs in the MCR or POMC genes and obesity or type 2 diabetes.     

Analysis of Ectonucleotide Pyrophosphatase/Phosphodiesterase 1 Gene K121Q Polymorphism Association with Some Risk Factors of Atherosclerosis in Patients with Acute Coronary Syndrome

The present study was performed to investigate whether common single-nucleotide polymorphism K121Q (rs1044498) of the ENPP1 gene is associated with the known risk factors of atherosclerosis (overweight, dyslipoproteinemia, hypertension, diabetes, smoking, and hypercoagulability) in persons with acute coronary syndrome. Venous blood of 118 patients was genotyped for the polymorphism by PCR and restriction fragment length polymorphism method. In patients divided into two subgroups according to their genotype (KK and KQ + QQ), the statistically significant differences were revealed only for plasma LDL-cholesterol level and fibrinolytic activity. The carriers of minor allele (KQ + QQ) had lower LDL-cholesterol concentration and the time of fibrinolysis than the major allele homozygotes (KK). The division of patients into subgroups according to presence or absence of some risk factors for atherosclerosis showed no statistically significant differences between K121Q genotype distributions for any of the comparison.     

The K121Q polymorphism of the plasma cell glycoprotein-1 gene is not associated with diabetes mellitus type 2 in German Caucasians.

The K121Q polymorphism of the human plasma cell membrane glycoprotein 1 (PC-1) gene is known to be associated with diabetes mellitus type 2 in some populations studied, with contradictory results. The purpose of the present study was to examine a possible association between the presence of diabetes and the PC-1 K121Q polymorphism in a German Caucasian population. Associations between the polymorphism and various metabolic and anthropometric parameters were also examined. The presence of the K121Q variant was investigated using polymerase chain reaction restriction fragment-length polymorphism in 402 subjects with diabetes (231 men, 171 women, age 63+/-11 yrs, body mass index 28.7+/-5.1 kg/m2) and in 432 age- and sex-matched controls (247 men, 185 women, age 64+/-7 yrs, BMI 26.5+/-3.7 kg/m2). Ninety-seven subjects were carriers of the K121Q polymorphism in the control and 110 in the diabetic group (allelic frequency 11.9% and 14.7%, respectively, P=0.25). The polymorphism had no significant influence on the presence of atherosclerotic disease, body mass index, and blood pressure, both, in diabetics and in non-diabetic controls. Our data suggest that the K121Q polymorphism of the PC-1 gene is not associated with diabetes, obesity, hypertension or atherosclerosis in a German Caucasian population.     

中国汉族人群ENPP1基因K121Q多态与2型糖尿病的关联及Meta分析

多个欧洲白人的Meta分析表明核苷酸焦磷酸酶1(Ectonucleotide pyrophosphatase/phosphodiesterase 1, ENPP1)基因K121Q多态与2型糖尿病相关联, 但在日本人、韩国人和中国台湾人的研究中没有发现相关性, 而在中国大陆人群中二者的关联研究结果不尽一致。文章调查了湖北地区539例2型糖尿病患者和404名正常人ENPP1基因K121Q多态性。基因型及等位基因频率在病例组和对照组间没有显著差异(P＞0.05), 但经性别、年龄和体重指数调整后的Logistic回归分析揭示XQ基因型与2型糖尿病显著相关(OR=1.5, 95%CI: 1.39~1.62, P<0.001)。对性别进行的亚组分析显示, 女性病例组Q等位基因和XQ基因型的频率显著高于对照组(Q: 12.4% vs. 6.1%, P=0.001; XQ: 23.7% vs. 11.7%, P=0.001)。结果表明ENPP1基因K121Q多态与湖北汉族人2型糖尿病的关联存在性别差异, 在女性中更明显。文章是对中国大陆人群进行的第一个Meta分析, 结果显示Q等位基因增加2型糖尿病的发病风险(OR=1.42, P=0.042)。     

ENPP1/PC-1 K121Q polymorphism and genetic susceptibility to type 2 diabetes in North Indians

Genetic susceptibility may be responsible for high prevalence of type 2 diabetes worldwide. A common missense single nucleotide polymorphism, K121Q in the ectoenzyme nucleotide pyrophosphate phosphodiesterase (ENPP1) gene, has recently been associated with type 2 diabetes in Italian, South Indian, and American populations. The objective of this study was to investigate the possible role of K121Q polymorphism in ENPP1 gene with type 2 diabetes in North Indians. The genotype of the ENPP1/PC-1 K121Q polymorphism was determined using polymerase chain reaction-restriction fragment length polymorphism analysis for 328 T2DM patients and 326 non-diabetic participants. Anthropometric and clinical characteristics (Body mass index (BMI), glucose, total cholesterol, triglyceride, high-density lipoprotein cholesterol (HDL), Creatinine, HbA1c, and insulin levels) were measured using standard protocols. Their Chi-square analyses were used to test the significance differences in genotypic and allelic frequencies. Association studies were undertaken using the t test and logistic regression analyses. Our results revealed there was no significant difference in the genotypic distribution between T2DM patients and control subjects. The KK and KQ genotype frequencies were similar in T2DM cases and controls (60.7 and 39.3% in T2DM and 59.8 and 40.2% in controls). No subjects with the QQ genotype were found. Binary logistic regression analysis of data did not show any association of K121Q polymorphism with type 2 diabetes (OR; 0.97, 95% CI; 0.7–1.32, P = 0.82). No significant correlation among the BMI, WHR, BP, TG, TC, HDL-C, LDL-C, Glucose, HbA1c, Creatinine, and insulin indices (HOMA-IR) was observed in the individuals carrying KK and KQ genotypes. In conclusion, our results showed that ENPP1/PC-1 K121Q polymorphism is not associated with type 2 diabetes and related quantitative metabolic traits in North Indian Punjabi population.     

Association of Genetic Variants of MTHFR,ENPP1, and ADIPOQ with Myocardial Infarction in Egyptian Patients

The study aimed to investigate the association between MTHFR C677T, ENPP1 K121Q, and ADIPOQ 45 T/G gene polymorphisms and incidence of myocardial infarction (MI) in Egyptian patients. The study included 60 unrelated patients suffering from their first MI and 60 unrelated controls. Patients were recruited from Kasr-El Eini hospital, Cairo University. The previously mentioned polymorphisms were determined in all participants by PCR–RFLP. There was no significant difference in the distribution of genotypes and alleles of MTHFR C677T between groups. In contrast, significant difference was found in the distributions of genotypes and alleles of ENPP1 K121Q and ADIPOQ 45 T/G between MI patients and controls (P = 0.01, P = 0.004, P = 0.009, P = 0.001, respectively). Univariate analysis revealed that 121Q ENPP1 and 45 G ADIPOQ alleles were associated with the increased risk of MI (OR = 3; 95 % CI = 1.45–6.2; P = 0.004 and OR = 5.8; 95 % CI = 1.92–17.54; P = 0.001, respectively). The mutant homozygous genotypes of MTHFR, ENPP1, and ADIPOQ were more prevalent in diabetic hypertensive MI patients than it was among non-diabetic normotensive MI patients. Regarding the coagulation profile, INR (P = 0.009) and PC % (P = 0.022) were significantly different among the three genotypes of MTHFR C677T. The 677 T, 121 Q, and 45G variants were associated with MI in Egyptian patients; however, more studies are needed to determine the possible protective effect for these polymorphisms in our population.     

Functional Consequences of a Novel Variant of PCSK1

Background

Common single nucleotide polymorphisms (SNPs) in proprotein convertase subtilisin/kexin type 1 with modest effects on PC1/3 in vitro have been associated with obesity in five genome-wide association studies and with diabetes in one genome-wide association study. We here present a novel SNP and compare its biosynthesis, secretion and catalytic activity to wild-type enzyme and to SNPs that have been linked to obesity.

Methodology/Principal Findings

A novel PC1/3 variant introducing an Arg to Gln amino acid substitution at residue 80 (within the secondary cleavage site of the prodomain) (rs1799904) was studied. This novel variant was selected for analysis from the 1000 Genomes sequencing project based on its predicted deleterious effect on enzyme function and its comparatively more frequent allele frequency. The actual existence of the R80Q (rs1799904) variant was verified by Sanger sequencing. The effects of this novel variant on the biosynthesis, secretion, and catalytic activity were determined; the previously-described obesity risk SNPs N221D (rs6232), Q665E/S690T (rs6234/rs6235), and the Q665E and S690T SNPs (analyzed separately) were included for comparative purposes. The novel R80Q (rs1799904) variant described in this study resulted in significantly detrimental effects on both the maturation and in vitro catalytic activity of PC1/3.

Conclusion/Significance

Our findings that this novel R80Q (rs1799904) variant both exhibits adverse effects on PC1/3 activity and is prevalent in the population suggests that further biochemical and genetic analysis to assess its contribution to the risk of metabolic disease within the general population is warranted.     

Amerindians show no association of PC-1 gene Gln121 allele and obesity: a thrifty gene population genetics

PC-1 Gln121 gene is a risk factor for type 2 diabetes, obesity and insulin resistance in European/American Caucasoids and Orientals. We have aimed to correlate for the first time this gene in Amerindians with obesity and their corresponding individuals genotypes with obesity in order to establish preventive medicine programs for this population and also studying the evolution of gene frequencies in world populations. Central obesity was diagnosed by waist circumference perimeter and food intake independent HDL-cholesterol plasma levels were measured. HLA genes were determined in order to more objectively ascertain participants Amerindians origin. 321 Amerindian blood donors who were healthy according to the blood doning parameters were studied. No association was found between PC-1 Gln121 variant and obesity. Significant HDL-cholesterol lower values were found in the PC-1 Lys121 bearing gene individuals versus PC-1 Gln121 bearing gene ones (45.1 ± 12.7 vs. 48.7 ± 15.2 mg/dl, p < 0.05). Population analyses showed a world geographical gradient in the PC-1 Gln121 allele frequency: around 9% in Orientals, 15% in European Caucasoids and 76% in Negroids. The conclusions are: (1) No association of PC-1 Gln121 gene is found with obesity in Amerindians when association is well established in Europeans. (2) PC-1 Gln121 gene is associated to higher levels of HDL-cholesterol than the alternative PC-1 Lys121 allele. This may be specific for Amerindians. (3) Amerindians have an intermediate frequency of this possible PC-1 Gln121 thrifty gene when compared with Negroid African Americans (78.5%) or Han Chinese (7.5%, p < 0.0001). Historical details of African and other groups may support the hypothesis that PC-1 Gln121 is indeed a thrifty gene.     

ENPP1 affects insulin action and secretion: evidences from in vitro studies

The aim of this study was to deeper investigate the mechanisms through which ENPP1, a negative modulator of insulin receptor (IR) activation, plays a role on insulin signaling, insulin secretion and eventually glucose metabolism. ENPP1 cDNA (carrying either K121 or Q121 variant) was transfected in HepG2 liver-, L6 skeletal muscle- and INS1E beta-cells. Insulin-induced IR-autophosphorylation (HepG2, L6, INS1E), Akt-Ser(473), ERK1/2-Thr(202)/Tyr(204) and GSK3-beta Ser(9) phosphorylation (HepG2, L6), PEPCK mRNA levels (HepG2) and 2-deoxy-D-glucose uptake (L6) was studied. GLUT 4 mRNA (L6), insulin secretion and caspase-3 activation (INS1E) were also investigated. Insulin-induced IR-autophosphorylation was decreased in HepG2-K, L6-K, INS1E-K (20%, 52% and 11% reduction vs. untransfected cells) and twice as much in HepG2-Q, L6-Q, INS1E-Q (44%, 92% and 30%). Similar data were obtained with Akt-Ser(473), ERK1/2-Thr(202)/Tyr(204) and GSK3-beta Ser(9) in HepG2 and L6. Insulin-induced reduction of PEPCK mRNA was progressively lower in untransfected, HepG2-K and HepG2-Q cells (65%, 54%, 23%). Insulin-induced glucose uptake in untransfected L6 (60% increase over basal), was totally abolished in L6-K and L6-Q cells. GLUT 4 mRNA was slightly reduced in L6-K and twice as much in L6-Q (13% and 25% reduction vs. untransfected cells). Glucose-induced insulin secretion was 60% reduced in INS1E-K and almost abolished in INS1E-Q. Serum deficiency activated caspase-3 by two, three and four folds in untransfected INS1E, INS1E-K and INS1E-Q. Glyburide-induced insulin secretion was reduced by 50% in isolated human islets from homozygous QQ donors as compared to those from KK and KQ individuals. Our data clearly indicate that ENPP1, especially when the Q121 variant is operating, affects insulin signaling and glucose metabolism in skeletal muscle- and liver-cells and both function and survival of insulin secreting beta-cells, thus representing a strong pathogenic factor predisposing to insulin resistance, defective insulin secretion and glucose metabolism abnormalities.     

Replication of association between polymorphisms of the pancreatic ATP-sensitive potassium channel and susceptibility to type 2 diabetes in two Russian urban populations

The KCNJ11 and ABCC8 genes encode components of the pancreatic ATP-sensitive potassium (KATP) channel. Previously, we reported association of the KCNJ11 E23K and ABCC8 R1273R G/A variants with type 2 diabetes (T2D) in a small Russian population sample (n=244). Here we replicated association between these genetic variants and T2D in a larger cohort (588 diabetic and 597 non-diabetic subjects). Using the ANCOVA analysis, Odds Ratios (ORs) and relationships between the carriage of a genotype and biochemical parameters of the patients were assessed and then adjusted for confounders (age, gender, HbA1c, hypertension, and obesity). The KCNJ11 K23 variant and the ABCC8 R1273R allele A showed association with higher risk of T2D (adjusted OR of 1.41 and 2.03, P<0.0001, respectively). Diabetic patients homozygous for K/K had lower 2h insulin (P_adjusted=0.044). The ABCC8 A/A variant was associated with increased 2h serum insulin in diabetic and non-diabetic subjects (P_adjusted=0.027 and 0.033, respectively). The carriage of the risk variant K/K of KCNJ11 E23K or A/A of ABCC8 G/A R1273R was associated with reduced response to nonsulfonylurea and sulfonylurea blockers of the pancreatic KATP channel. Adjusted attributable population risk was 3.0% (KCNJ11 E23K) and 4.8% (ABCC8 G/A) suggesting for the modest effects of these genetic variants on diabetes susceptibility.     

Analysis of the SIM1 Contribution to Polygenic Obesity in the French Population

SIM1 (single‐minded 1) haploinsufficiency is responsible for obesity in both humans and mice, but the contribution of frequent DNA variation to polygenic obesity is unknown. Sequencing of all exons, exon/intron boundaries, 870 base pairs (bp) of the putative promoter, and 1,095 bp of the 3′UTR of SIM1 gene in 143 obese children and 24 control adults identified 13 common variants. After analysis of the linkage disequilibrium (LD) structure, association study of eight variants was performed in 1,275 obese children and severely obese adults, in 1,395 control subjects, and in 578 obesity‐selected pedigrees. A nominal evidence of association was found for the nonsynonymous variant P352T C/A (rs3734354) (P = 0.01, OR = 0.81 (0.70–0.95)), the +2,004 TGA ?/insT SNP (rs35180395) (P = 0.02, OR = 1.21 (1.02–1.43)), the +2,215A/G TGA SNP (rs9386126) (P = 0.002, OR = 0.81 (0.71–0.93)), and pooled childhood/adult obesity. Even though transmission disequilibrium test (TDT) further supported the association of P352T and +2,004 ?/inst T with obesity, none of these nominal associations remained significant after a multiple testing Bonferroni correction. Therefore, our study excludes a major contribution of SIM1 common variants in exons, 5′ and 3′ UTR regions in polygenic obesity susceptibility in French Europeans.     

The R453Q and D151A polymorphisms of hexose-6-phosphate dehydrogenase gene (H6PD) influence the polycystic ovary syndrome (PCOS) and obesity

Hexose-6-phosphate dehydrogenase (H6PDH) influences 11β-hydroxysteroid dehydrogenase activity, a key enzyme in the peripheral metabolism of cortisol that modulates insulin sensitivity in adipose tissue. To study the associations of R453Q and D151A polymorphisms in the H6PDH gene (H6PD) with polycystic ovary syndrome (PCOS) and their influence on clinical and metabolic variables, we genotyped 237 patients with PCOS and 135 control women for the R453Q (rs6688832) and D151A (rs34603401) variants in H6PD. The R453Q genotypes were distributed differently in patients and controls (χ(2)=9.55, P=0.002). Genotypes of D151A were distributed evenly in women with PCOS and controls, but showed a different distribution in non-obese and obese women (χ(2)=3.95, P=0.047), especially within the PCOS subgroup (χ(2)=4.65, P=0.031). A backward stepwise likelihood ratio logistic regression model (Nagelkerke's R(2)=0.490; χ(2)=164; P<0.0001) retained free testosterone (OR=1.13; 95% CI: 1.10-1.17) and H6PD Q453 alleles (OR=0.46; 95% CI: 0.27-0.79) as statistically significant predictors for PCOS, whereas homeostasis model assessment of insulin resistance and the H6PD D151A variant were excluded by the model. Both H6PD variants were associated with several phenotypic variables, including fasting insulin, homeostasis model assessment of insulin resistance and androstenedione levels. In summary, the R453Q and D151A variants of the H6PD gene are associated with PCOS and obesity, respectively, and may contribute to the PCOS phenotype by influencing obesity, insulin resistance and hyperandrogenism.     

Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r² ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4⁺ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.     

Multiple Autoimmune-Associated Variants Confer Decreased IL-2R Signaling in CD4+CD25hi T Cells of Type 1 Diabetic and Multiple Sclerosis Patients

IL-2 receptor (IL-2R) signaling is essential for optimal stability and function of CD4⁺CD25^hiFOXP3⁺ regulatory T cells (Treg); a cell type that plays an integral role in maintaining tolerance. Thus, we hypothesized that decreased response to IL-2 may be a common phenotype of subjects who have autoimmune diseases associated with variants in the IL2RA locus, including T1D and MS, particularly in cells expressing the high affinity IL-2R alpha chain (IL-2RA or CD25). To examine this question we used phosphorylation of STAT5 (pSTAT5) as a downstream measure of IL-2R signaling, and found a decreased response to IL-2 in CD4⁺CD25^hi T cells of T1D and MS, but not SLE patients. Since the IL2RArs2104286 haplotype is associated with T1D and MS, we measured pSTAT5 in controls carrying the rs2104286 risk haplotype to test whether this variant contributed to reduced IL-2 responsiveness. Consistent with this, we found decreased pSTAT5 in subjects carrying the rs2104286 risk haplotype. Reduced IL-2R signaling did not result from lower CD25 expression on CD25^hi cells; instead we detected increased CD25 expression on naive Treg from controls carrying the rs2104286 risk haplotype, and subjects with T1D and MS. However the rs2104286 risk haplotype correlated with increased soluble IL-2RA levels, suggesting that shedding of the IL-2R may account in part for the reduced IL-2R signaling associated with the rs2104286 risk haplotype. In addition to risk variants in IL2RA, we found that the T1D-associated risk variant of PTPN2rs1893217 independently contributed to diminished IL-2R signaling. However, even when holding genotype constant at IL2RA and PTPN2, we still observed a significant signaling defect in T1D and MS patients. Together, these data suggest that multiple mechanisms converge in disease leading to decreased response to IL-2, a phenotype that may eventually lead to loss of tolerance and autoimmunity.     

Association between MSX1 variants and oral clefts in Han Chinese in western China

Previous animal and association studies have shown that the MSX1 gene is associated with oral clefts. Our aim was to investigate association between variants in the MSX1 gene and oral clefts in a Han Chinese population. Our study group consisted of 206 nonsyndromic oral cleft (NSOC) nuclear families (including the patients and their parents) and 224 controls. The three variants evaluated in this study were single-nucleotide polymorphisms rs3821949 and rs12532 and a missense mutation P147Q. Polymerase chain reaction-restriction fragment length polymorphism was used to genotype the three markers. Case-control and family-based association analyses were carried out. In the case-control analysis, no significant differences in genotypic or allelic frequencies were observed in any of the two single-nucleotide polymorphisms between patients and controls. Although the homozygous T allele for P147Q was not detected in any sample in this population, heterozygotes were more prevalent in NSOC (1.2%) when compared with the controls (0%). The analyses for family-based association did not suggest association between any of the three variants and NSOC. No significant association was found between NSOC and rs3821949 or rs12532 in MSX1 gene, whereas an association was observed between the P147Q variant and cleft lip with cleft palate in the case-control analysis.