|
|
||||||||
is associated with type 2 diabetes in Swedish CaucasiansDepartment of Molecular Medicine and Surgery, Karolinska Institutet, Rolf Luft Center for Diabetes Research, Karolinska University Hospital Solna, 171 76 Stockholm, Sweden
(Correspondence should be addressed to E-C Långberg; Email: ewa-carin.langberg{at}ki.se)
| Abstract |
|---|
|
|
|---|
(RPTP
) is involved in glucose homeostasis and insulin signaling in several animal models. The aim of this study was to evaluate whether polymorphisms in this gene influence the development of type 2 diabetes (T2D) in humans.
Design: We investigated how genetic variations in the RPTP
gene influence susceptibility to impaired glucose tolerance (IGT) and T2D, in Swedish men and women.
Methods: Genotyping of single nucleotide polymorphisms (SNPs) was performed by dynamic allele-specific hybridization in a total of 1057 Swedish Caucasians including 497 subjects with normal glucose tolerance (NGT), 262 subjects with IGT, and 298 patients with T2D.
Results: SNPs rs1143699, rs4807015, and rs1978237 were found to be associated with T2D. SNP rs1143699 was associated with male T2D patients when compared with NGT controls (odds ratio; OR = 1.57; P = 0.029). SNP rs4807015 showed association with T2D patients when compared with NGT controls (OR = 1.32; P = 0.025). Finally, SNP rs1978237 was associated with T2D patients when compared with NGT controls (OR = 1.59; P = 0.002). Logistic regression analysis demonstrated that for SNP rs1143699 in men, C/C homozygosity conveys an increased risk of T2D (OR = 2.19; P = 0.035), while SNP rs4807015 was associated with an increased risk of T2D in both men and women (OR = 1.74; P = 0.029). SNP rs1978237 also demonstrated a risk of T2D in men and women (OR = 1.59; P = 0.026).
Conclusions: This study provides evidence for association of SNPs in the RPTP
gene with T2D in Swedish Caucasians. SNPs rs1143699, rs4807015, and rs1978237 confer an increased risk of developing T2D.
| Introduction |
|---|
|
|
|---|
We have previously demonstrated that the receptor protein tyrosine phosphatase
(RPTP
) gene, also known as PTP
, is overexpressed by ~60% in pancreatic islets and liver of spontaneously diabetic Goto–Kakizaki (GK) rats that have impaired insulin secretion. When islet RPTP
was inhibited by antisense, improved glucose-induced insulin secretion was seen in GK islets (5). In addition, RPTP
knockout (–/–) mice have shown decreased plasma glucose and insulin levels in the fasted state when compared with wild-type controls. The mice also had increased whole-body insulin sensitivity, suggesting that RPTP
affects insulin signaling in insulin-sensitive tissues (6).
RPTP
has a function in the cellular receptor signaling, which is critical for signal transduction in both normal and pathophysiological conditions (7). It is a member of the leukocyte antigen-related (LAR) RPTP family, which has been suggested to act in key steps of neural development and also in diabetes and cancer (8). LAR, a protein closely related to RPTP
, has negative regulatory effects in the insulin signaling pathway when it is overexpressed (9). The RPTP
gene gives rise to several different distinct isoforms (10). PTPs are key regulators of the insulin receptor signal transduction pathway and the RPTP
gene has been shown to be highly expressed in insulin target tissues, such as liver, adipose tissue, skeletal muscle, and endothelial cells (11). In humans, the RPTP
gene is located on chromosome 19p13.3, a region that may influence traits underlying lipid abnormalities associated with T2D (12).
Until now, there has not been a genetic report of RPTP
in T2D patients. We have selected SNPs in the RPTP
gene to study in a Swedish cohort including T2D patients, subjects with impaired glucose tolerance (IGT), and normal glucose tolerance (NGT) control subjects. Our study hypothesis was to evaluate the genetic influence of polymorphisms of the RPTP
gene in development of T2D among Swedish men and women.
| Subjects and methods |
|---|
|
|
|---|
The study includes 497 unrelated subjects with NGT (246 men/251 women), 262 subjects with IGT (107/155), and 298 patients with T2D (241/57). All individuals were Swedish Caucasians collected from the Stockholm Diabetes Prevention Program with a baseline study in men 1992–1994 and in women 1996–1998. A follow-up study in about 70% of the baseline participants was conducted 10 years later (13–15). Patients were diagnosed according to the World Health Organization criteria (1998) (16). NGT control subjects had a mean birth body weight of 3742 g (3000–5800), body mass index
26 kg/m2, and no relatives diagnosed with T2D. All clinical data in T2D patients were used from the baseline study (those diagnosed at baseline and incident cases between baseline and follow-up studies) or at follow-up study (diagnosed at follow-up). The data were not influenced by use of anti-diabetic medication. Clinical parameters for each study group are given in Table 1
. All participants gave their informed consent to take part in the study. The procedures followed were in accordance with the declaration of Helsinki II and approved by the ethics committee of Karolinska Hospital. DNA was extracted from peripheral blood using Puregene DNA purification kit (Gentra Systems, Minneapolis, MN, USA).
|
We selected a total of 28 SNPs within the RPTP
gene from the National Center for Biotechnology Information (USA; http://www.ncbi.nlm.nih.gov/SNP/) database based on validation status, region, and function for our study. SNPs chosen were picked to cover the whole gene. The Tagger program from the International HapMap project (release no. 21) was also used to select and evaluate tagSNPs from genotype data in HapMap (17). Here, pairwise tagging was used together with an r2 cutoff of 0.8 and a minor allele frequency of 5%. The majority of the chosen SNPs were captured by tagger.
Validation of 28 chosen SNPs in 32 Swedish DNA samples including 16 T2D patients and 16 NGT controls was performed. Twelve SNPs representing at least 4% allele frequencies were used for further genotyping. Genotyping was performed by dynamic allele-specific hybridization (DASH) (18). PCR-DASH assay design and SNP genotyping protocol were used as described previously (19). Information on primer and probe sequences used in PCR-DASH assays are available, upon request. All PCRs were run in 20 µl volumes using 10–20 ng genomic DNA. Negative controls (water blanks) were included on each plate. For genotyping quality control, the subjects were distributed randomly across plates with the same number of cases and controls per plate. Successful genotype calls were
95% and plates were randomly genotyped twice for duplication accuracy, which was calculated to be 98%.
Calculations and statistical analysis
The homeostasis model assessment (HOMA) was used to assess insulin resistance. Based on fasting glucose and insulin levels, it was calculated according to the following equation: fasting plasma glucose (mmol/l) x fasting plasma insulin (mU/ml)/22.5 (20).
Normal probability plots were created and parameter distributions were transformed to natural logarithm as required to improve skewness and to obtain a normal distribution before performing statistical analysis. Homogeneity of variances was tested by Levenes test. Genotype and allele distribution were compared between cases (IGT subjects and/or T2D patients) and NGT control subjects using
2-test. OR and 95% confidence intervals (CI) were calculated to test for relative risk for association. Testing for differences in clinical parameters between genotype groups was performed using ANOVA and covariance analysis adjusting for age. P values less than 0.05 were interpreted as statistically significant. Multivariate logistic regression analysis was used to predict the susceptibility to T2D between case and control subjects. Adjustments for age, sex, BMI, and blood pressure (where appropriate) were applied. The effects of the SNPs were tested individually and/or as haplotypes. Linkage disequilibrium (LD) and haplotype frequencies were calculated using Arlequin (21). LD between different markers/SNPs was summarized using |D'|. Haplotypes prevalent at 5% were used for further haplotype analysis. All statistical analyses were performed using Statistica version 7.1 (Statsoft Inc., Tulsa, OK, USA) and/or Biomedical package (BMDP) version 1.1 (BMDP Statistical Software Inc., Los Angeles, CA, USA).
| Results |
|---|
|
|
|---|
gene and all of them were in Hardy-Weinberg equilibrium (P > 0.05). Genotype and allele frequencies of all valid SNPs in NGT and IGT subjects and T2D patients are presented in Table 2
|
SNP rs1143699 showed significant allelic association with male T2D patients when compared with NGT controls (0.89 vs 0.84), where the major allele C was present at a higher frequency in T2D (P = 0.029; OR = 1.57; 95% CI = 1.05–2.36; Table 3
). No allelic association could be noted when comparing male IGT subjects and controls, or in female, IGT and T2D, when compared with controls. Comparison of genotypes and quantitative traits within groups was performed by ANOVA and covariance analysis adjusting for age. No significant associations with quantitative traits were detected (data not shown). Multivariate logistic regression analysis of SNP rs1143699 with adjustments for age, BMI, and blood pressure indicated that in men the C/C genotype was significantly associated with T2D when compared with C/T+T/T (P = 0.035; OR = 2.19; 95% CI = 1.055–4.523; Table 4
).
|
|
A significant difference in allele frequency was found in SNP rs4807015 comparing all T2D patients with NGT controls (0.49 vs 0.56), C being the major allele present at a lower frequency in T2D (P = 0.025; OR = 1.32; 95% CI = 1.034–1.678; Table 3
). Multiple logistic regression analysis showed that this SNP was associated with increased risk of T2D in men and women carrying the C/C genotype (P = 0.029; OR = 1.74; 95% CI = 1.058–2.865) adjusted for age and blood pressure (Table 5
). When adding sex as a covariate to the model, this association was no longer seen. In IGT men and women subjects, the association was not significant (data not shown).
|
SNP rs1978237 showed significant allelic association with all T2D patients when compared with NGT controls (0.80 vs 0.72), where the major allele G was present at a higher frequency in T2D (P = 0.002; OR = 1.59; 95% CI = 1.184–2.126; Table 3
). Multiple logistic regression analysis showed that this SNP was associated with increased risk of T2D in men and women carrying the G/G genotype (P = 0.026; OR = 1.591; 95% CI = 1.056–2.398) adjusted for age, BMI, and blood pressure (Table 6
). A similar result was retrieved when adding the covariate sex to the analysis. In IGT men and women subjects, the association was not significant.
|
LD between SNPs was found to be in concordance with HapMap data (release no. 21). We found three common haplotypes for SNP rs1143699, rs4807015, and rs1978237 with a frequency above 5%. However, there was no significant difference in haplotype distribution between groups for SNP rs1143699, rs4807015, and rs1978237 (data not shown).
| Discussion |
|---|
|
|
|---|
in Swedish NGT subjects, IGT subjects, and T2D patients, and found that three polymorphisms confer susceptibility to development of T2D. SNP rs1143699 was associated with an increased risk of T2D in male patients carrying the C/C genotype. SNP rs4807015 and rs1978237 was associated with an increased risk of T2D in both male and female patients. The moderate OR for these three SNPs, which were around 1.3–1.6, support the hypothesis that there is restricted influence of a single SNP in T2D. However, when comparing sample groups and adjusting for clinical features separating them (age, BMI, and blood pressure), the relative risk to develop T2D was even more increased for SNP rs1143699 (OR 2.2), but still the same for SNP rs1978237 (OR 1.6). SNP rs4807015 showed the same trend in results. This is a truer reflection of how the different SNPs studied contribute to susceptibility of T2D, since the data have been adjusted for covariates. In all three SNPs, the allele being responsible for the increased risk was the major allele.
SNP rs1143699 is located in exon 33, which is within the PTP catalytic (PTPc) domain, an important functional part of the RPTP
gene being involved in phosphorylation of interacting proteins. SNP rs4807015 is located in intron 13 and rs1978237 in intron 34. In other genes, intronic SNPs have been suggested to contribute to gene expression, and therefore, this is a possible explanation for the role of SNP rs4807015 and rs1978237 (22–24). Finally, susceptibility is often a result of the combined action of many different variants in a gene, one example being calpain-10 in T2D (25). These considerations need to be emphasized when interpreting association data of RPTP
.
We did not detect any association between genotype groups and clinical phenotypes for the SNPs in our study. This is most probably because our male patient group contains clinical data at baseline also for patients who were diagnosed with T2D during the 10 years between baseline and follow-up studies. This subgroup of T2D men did not have diabetes at baseline, which explains the low fasting plasma glucose in the diabetes group (6.6 mmol/l). However, we chose not to include follow-up data of these patients, since anti-diabetic drugs and lifestyle changes would then be affecting the quantitative traits of the patients.
LAR is a member of the same PTP family as RPTP
and has a role in negative modulation of insulin receptor signaling (26). Increased activity of LAR or related PTPs in insulin target tissues could contribute to development of insulin resistance (27). The fact that RPTP
shows similarity to LAR, 85% at amino acid level, and is expressed to a higher level than LAR in insulin-sensitive tissues, makes this protein a possible candidate for regulating insulin signaling and for playing a role in glucose homeostasis (28). To our knowledge, the present study is the first genetic investigation of RPTP
in T2D. However, there is a genetic association study of an SNP in the promoter region of LAR in two non-diabetic populations in Italy (29), showing association with diabetes-related characteristics (29). Thus, LAR and RPTP
seem to be important genes in the development of T2D.
SNP rs1143699 was found to be associated with T2D in men but not in women, suggesting that gender differences may be an important factor in the genetic influence of RPTP
in T2D. Our group and others have previously shown that there are sex differences in genetic association with T2D. Polymorphisms in the insulin-degrading enzyme gene are associated with increased BMI and insulin resistance in men but not in women (30). Also, a polymorphism in neuropeptide Y, a gene involved in satiety and body weight, is related to IGT and T2D in men (31). Thus, sex specificity could be an important confounding factor in analyzing genetics of T2D and related disorders (30–34).
In conclusion, this study provides evidence that three polymorphisms in the RPTP
gene are associated with the development of T2D in Swedish Caucasians, which support our hypothesis that the RPTP
gene may be important for the development of T2D in humans. Further investigation in other populations and determination of the functional influence of variation in RPTP
are necessary to make conclusive implications.
| Acknowledgements |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Hofmeyer and J. E. Treisman The receptor protein tyrosine phosphatase LAR promotes R7 photoreceptor axon targeting by a phosphatase-independent signaling mechanism PNAS, November 17, 2009; 106(46): 19399 - 19404. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |