From the Research Division, Joslin Diabetes Center, and the Department of Medicine, Harvard Medical School, Boston, Massachusetts
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
A polymorphism in the ecto-nucleotide pyrophosphatase/phosphodiesterase 1 gene (ENPP1), previously known as PC-1, has been demonstrated to be associated with insulin resistance (5). ENPP1 is one of five cell membrane proteins with an extracellular active site catalyzing the release of nucleoside 5'-monophosphatase from nucleotides and their derivatives. These proteins consist of a short NH2-terminal intracellular domain, a single transmembrane domain, two somatomedin-Blike domains, a catalytic domain, and a COOH-terminal nuclease-like domain. Three of these proteins (ENPP13) have been characterized and implicated in various processes, including bone mineralization, signaling by insulin and by nucleotides, and the differentiation and motility of cells (6). However, the precise mechanisms through which the proteins affect these processes are not known.
ENPP1 is by far the best-characterized ENPP. It was first discovered as a surface marker of antibody-secreting B-cells, hence the name plasma-cell differentiation antigen-1 or PC-1 (6). Now, however, it is known that ENPP1 is expressed in various tissues, including muscle, fat, liver, and kidney (7). The physiologic functions of ENPP1 in these tissues have not been described. However, it has been found that overexpression of ENPP1 inhibits insulin receptor tyrosine kinase activity in various cells (810) and causes insulin resistance (8,11,12). ENPP1 seems to inhibit insulin signal transduction by interacting directly with the insulin receptor subunit (12).
The gene encoding for ENPP1 has 25 exons and is located on human chromosome 6q2223 (5,13). Recently, Pizzuti et al. (5) described a DNA polymorphism in exon 4 that causes an amino acid change from lysine to glutamine at codon 121 (K121Q). This amino acid change is located in the second somatomedin-Blike domain of ENPP1 and may interfere with protein-protein interactions (6). Studies in vitro have shown that the Q variant of ENPP1 interacts more strongly with the insulin receptor than the K variant and reduces insulin receptor autophosphorylation (5). The Q variant also reduces insulin-induced phosphorylation of insulin receptor substrate (IRS)-1, phosphatidylinositol-3 kinase activity, glycogen synthesis, and cell proliferation (14). Furthermore, human studies demonstrated that Sicilian carriers of the Q variant had lower insulin sensitivity than noncarriers. Similar results were obtained in studies of Finns and Swedes (15) but not of Danes (16).
De Cosmo et al. (17) described an effect of the ENPP1 polymorphism on the rate of loss of renal function in Caucasians with type 1 diabetes and proteinuria. During 6.5 years of follow-up, the glomerular filtration rate (GFR) declined more rapidly in carriers of the Q variant (QQ/KQ, n = 22) than in noncarriers (KK, n = 55). The median rates of decline were 7.2 and 3.7 ml · min1 · year 1, respectively. With such rapid loss of renal function, diabetic carriers of the Q variant would progress through the stage of proteinuria to end-stage renal disease (ESRD) in fewer years than noncarriers. Therefore, the effect of this variant would be most evident when ESRD develops early in the course of diabetes.
The objective of this study was to investigate the association of advanced diabetic nephropathy with the K121Q polymorphism in ENPP1. Specifically, we sought to test whether the Q variant confers an increased risk of the development of ESRD early in the course of type 1 diabetes. Both case-control and family-based (transmission disequilibrium test [TDT]) study designs were used.
![]() |
Research Design and Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
As unrelated control subjects for genetic studies on diabetic nephropathy, we selected Caucasian patients with normoalbuminuria and type 1 diabetes duration 15 years. Among the 1,080 normoalbuminuric participants in the Natural History of Microalbuminuria Study, 420 met these additional criteria, and we examined 307 (73%) of them.
As unrelated cases for the genetic studies on diabetic nephropathy, we selected Caucasian patients with persistent proteinuria or ESRD. We were able to enroll and examine 152 of the 206 patients with advanced nephropathy identified in the Natural History Study described above. In addition, we enrolled and examined 201 (86%) of 235 Caucasian patients with type 1 diabetes and proteinuria or ESRD attending the clinic between 1994 and 2000 who had not been included in the 50% random sample screened in 19911993 for the Natural History Study. The two groups of patients did not differ with regard to sex, age at examination, duration of diabetes, or values of HbA1c, total serum cholesterol, or serum triglycerides. Therefore, we pooled the two groups, bringing the total number of unrelated subjects with advanced diabetic nephropathy to 352.
For the family-based association study, both parents of 174 patients and both parents of 78 control subjects were examined.
Examination of study participants.
After consenting to participate in the study, each subject had a standardized physical examination and provided a diabetes history (diagnosis, treatment, and complications). Each individual provided a blood sample for biochemical measurements and DNA extraction. For the family-based study, parents provided blood samples for DNA extraction. The Committee on Human Subjects of the Joslin Diabetes Center approved the protocol and informed consent procedures.
Diagnosis of diabetic nephropathy.
The diabetic nephropathy status of each patient was determined on the basis of questionnaires, medical records of the Joslin Clinic (supplemented with records of other physicians if necessary), and measurements of the ACR. Methods for measuring the ACR have been described previously (18). Patients were classified as control subjects (n = 307) if they had diabetes duration 15 years and ACR <17 mg/g (men) or <25 mg/g (women) in at least two of the three most recent urine specimens. Patients with microalbuminuria or intermittent proteinuria were excluded from the study. Patients were considered case subjects (n = 352) if they had persistent proteinuria or if they had ESRD due to diabetic nephropathy. Persistent proteinuria was defined as two of three successive urinalyses positive by reagent strip (
2+ on Multistix; Bayer, Elkhart, IN) or as ACR
250 mg/g (men) or
355 mg/g (women). At the time of the first examination (19922000), subjects were divided into those with persistent proteinuria (n = 260) and those with ESRD (n = 92). Patients with proteinuria were followed until the end of July 2001, by which time ESRD had developed in 60 of them, i.e., renal dialysis had been initiated. These patients were considered incident case subjects of ESRD. The individuals who had begun dialysis or received a renal transplant before the first examination were considered prevalent case subjects of ESRD. For prevalent or incident cases of ESRD, the year when dialysis began was used to calculate the duration of diabetes. In total, we had 152 subjects with ESRD (92 prevalent and 60 incident).
Genotyping protocol.
The DNA polymorphism in exon 4 of ENPP1 has been described previously (5). A transversion from A to C at the first position of codon 121, resulting in an amino acid change from lysine to glutamine (K121Q), was genotyped with a previously described protocol (5). PCR was performed using the described conditions and primers: forward 5'-CTGTGTTCACTTTGGACATGTTG-3' and reverse 5' -GACGTTGGAAGATACCAGGTTG-3'. PCR products were digested by the restriction enzyme AvaII and separated on agarose gel. The allele encoding the K variant appeared as a single fragment of 238 bp, and the allele encoding the Q variant appeared as two fragments of 148 and 90 bp each.
Q carriers (KQ and QQ) were analyzed together, given that both genotypes were previously described to be associated with faster decline in GFR (17).
Analytic methods.
The data from this study were analyzed first as a case-control study. The characteristics of the case and control subjects were compared by 2 test or Fishers exact test for frequencies and one-way ANOVA for continuous variables. Odds ratios (ORs) and 95% CIs for persistent proteinuria and ESRD were estimated for carriers of the Q variant (19). Logistic regression analysis and Breslow-Day test were performed in SAS (version 8.02 for Windows; SAS Institute, Cary, NC).
Transmission of the risk variant from heterozygous parents to offspring was assessed by TDT (20). Based on an a priori hypothesis that the Q variant confers risk for diabetic nephropathy, a one-sided exact test based on binomial distribution was used to test the transmission of the risk allele. We ascertained the total number of heterozygous parents (i.e., the number of Bernoulli trials) and the total number of these parents transmitting the Q allele by direct counting and calculated statistical significance using the standard Excel spreadsheet function (BINOMDIST).
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
To examine the hypothesis that the Q variant is associated with the early development of advanced diabetic nephropathy in type 1 diabetes, we divided each of the study groups at the median duration of diabetes, 24 years (Table 2). Among subjects with diabetes duration <24 years, carriers of the Q variant had a significantly higher overall risk of advanced diabetic nephropathy (proteinuria and ESRD combined) than noncarriers: OR 1.9 (95% CI 1.13.3). Moreover, in agreement with our hypothesis, the OR for ESRD was 2.3 (95% CI 1.24.4) and only 1.6 for persistent proteinuria, the latter being nonsignificant (95% CI 0.92.9). Among subjects with diabetes duration 24 years, carriers of the Q variant had a higher overall risk of advanced diabetic nephropathy than noncarriers, but the OR was not significantly different from 1: OR 1.6 (95% CI 0.92.7). However, their risk of persistent proteinuria was significantly increased (OR 1.9, 95% CI 1.13.3), whereas their risk of ESRD was not (OR 1.3, 95% CI 0.72.4). With the sample sizes available, the OR for ESRD in subjects with diabetes duration <24 years was not significantly different from that in subjects with duration
24 years (Breslow-Day test). Similarly, when we used a multiple logistic model with duration of diabetes as a continuous variable instead of as dichotomous, the interaction between the presence of the Q variant and diabetes duration was still not statistically significant (data not shown).
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Because case-control findings may be biased by unrecognized population stratification, a family-based approach, the TDT, can be used to examine the association of the Q variant with diabetic nephropathy. For this test, DNA from both parents is examined, and only those families with at least one heterozygous parent are used in the analysis. Excess transmission of a variant beyond the expected 50% confirms its role as a risk allele. In agreement with the results of the case-control study, transmission of the allele that encodes the Q variant to index subjects with early-onset ESRD exceeds the expected 50%. In contrast, transmission of the Q variant from heterozygous parents to subjects with persistent proteinuria is consistent with the expected 50%. Thus, if the Q variant does have an impact on the risk of proteinuria, the effect is too weak to be detected in our relatively small number of TDT families.
Our findings can be interpreted as supporting two alternative disease models. The first invokes two effects of the Q variant: 1) carriers of Q have a moderately increased risk of developing persistent proteinuria, regardless of the duration of diabetes, and 2) the Q variant interacts in some people with a factor to accelerate progression of proteinuria to ESRD. Curiously, that factor is relatively isolated to individuals with short-duration diabetes. This model attributes the negative TDT analysis in late-onset proteinuria to a lack of power (type 2 error). The second model is simpler and, therefore, considered more likely: the Q variant accelerates the natural history of diabetic nephropathy (earlier onset and faster progression through each stage), so the association is strongest in patients with the earliest onset of ESRD. This model attributes the association of the Q variant with late-onset proteinuria in the case-control study to type 1 error and accepts the negative TDT analysis as truth rather than type 2 error. With the available evidence, distinction between these two models is not possible. Much larger groups for case-control comparisons and a much larger set of families, with both parents examined for TDT analysis, will be required.
The molecular mechanisms accounting for the association between the Q variant of ENPP1 and the development of advanced stages of diabetic nephropathy can only be hypothesized. ENPP1 is expressed in various tissues, including muscle, fat, liver, and kidney (7). Recently, we demonstrated that this gene is also expressed in human mesangial and endothelial cells (L.H.C., D.P.K.N., A.S.K., unpublished data). This pattern of expression is compatible with the localization of diabetic nephropathy disease processes. The DNA sequence difference in exon 4 of ENPP1 results in an amino acid change from lysine to glutamine at codon 121 (K121Q). In a recent review, Bollen et al. (6), using a different transcriptional start site, suggested that the actual position of the polymorphism is at codon 173 (K173Q). This amino acid change is located in the second somatomedin-Blike domain of ENPP1, and it seems to interfere with the homodimerization of ENPP1 and interactions between ENPP1 and other proteins (6).
For example, studies in vitro show that the Q variant of ENPP1 interacts more strongly with the insulin receptor than the common K variant. Because of this, cells with the Q variant have reduced insulin receptor autophosphorylation (5). The Q variant also reduces insulin-induced phosphorylation of the insulin receptor substrate (IRS)-1, phosphatidylinositol 3-kinase activity, glycogen synthesis, and cellular proliferation (22). In humans, carriers of the Q variant are more insulin resistant and hyperinsulinemic than noncarriers (15). Hyperinsulinemia may stimulate renal sodium reabsorption, leading to volume expansion, enhanced sympathetic adrenergic activity, and upregulation of the angiotensin II type 1 receptor, causing impaired peripheral vasodilation (23). Volume overload and impaired peripheral vasodilation may predispose to blood pressure elevation and decrease nocturnal dipping of the blood pressure (23,24). These are known phenotypic characteristics of ensuing diabetic nephropathy. However, in agreement with De Cosmo et al. (17), carriers and noncarriers of the Q variant in our study did not differ with regard to phenotypic characteristics of the insulin resistance syndrome, such as BMI, blood pressure, lipid profile, or HbA1c values. However, because we did not measure insulin sensitivity directly, we cannot exclude the possibility that the early development of ESRD in carriers of the Q variant might be accounted for by factors associated with insulin resistance.
Similar to the insulin receptor, ENPP1 may be a binding partner of integrins, which may be implicated in the pathologic processes involved in progression of diabetic nephropathy. Integrins are extracellular-matrix proteins involved in cell adhesion (25). Hyperglycemia is reported to influence expression of integrins in mesangial cells (26), and alterations in expression of these molecules has been hypothesized to be involved in the pathogenesis of diabetic nephropathy (27). The Q variant of ENPP1 may interfere with this putative process, leading in some unknown way to increased risk of diabetic nephropathy.
Also noteworthy is the cysteine richness of somatomedin-Blike domains, which may be involved in the dimerization of ENPPs via disulfide bonds (6). However, this has not been studied, and the impact of the Q variant on structure and function of dimerized ENPP1 is unknown. ENPP1 does regulate the concentration of pharmacologically active extracellular compounds such as adenosine and its derivatives, and they affect cell motility (28). Other possible functions of ENPP1 are associated with sulfation of glycosaminoglycans and regulation of cellular proliferation through modulation of local concentrations of ATP. Any of these functions may be relevant to the development of diabetic nephropathy in patients with type 1 diabetes.
Finally, we cannot exclude the possibility that ENPP1 polymorphism is in linkage disequilibrium with flanking DNA sequence differences that may be responsible for susceptibility to diabetic nephropathy. Screening of 25 exons of this gene has not yielded any DNA sequence differences beyond the polymorphism in exon 4 that was investigated in this report (17). A relatively infrequent haplotype (6.8% in the general population) in the 3' end of the gene is associated with insulin resistance (22). This sequence difference was not examined by us, but confounding of our findings by this haplotype is unlikely because it is in linkage equilibrium with the polymorphism in exon 4 (22).
To investigate the possibility that our findings could be explained by linkage disequilibrium with DNA sequence difference in flanking genes, 50-kb sequences 5' and 3' from the ENPP1 locus were searched in silico (http://www.ncbi.nlm.nih.gov/cgi-bin/Entrez/map_search). No known genes were found in this region. Although three human UniGene EST clusters were identified (Hs.25557, Hs.44230, and Hs.183276), BLAST searches (http://www.ncbi.nlm.nih.gov/BLAST) did not reveal any similarity between these clusters and any known gene or protein. However, one gene from the nucleotide pyrophosphatase/phosphodiesterase family, ENPP3, is located 60 kb centromeric to ENPP1. ENPP3 has 25 exons and spans
60 kb (29). No systematic search for DNA polymorphisms in this gene has been conducted. Further work is needed to exclude the possibility that DNA sequence differences in this gene could be in linkage disequilibrium with the Q variant of ENPP1.
Because patients with diabetic nephropathy have increased mortality (30,31), the possibility of spurious findings due to survival bias has to be evaluated, especially when a majority of patients are prevalent case subjects. In our study, this would imply that the excess of Q carriers in subjects with early-onset ESRD was a spurious result due to better survival of Q carriers than KK homozygotes. While definitive exclusion of this possibility is feasible only in a long-term follow-up study, several features of our investigation provide strong arguments against it. First, we studied relatively young individuals with ESRD, in whom mortality is still low. Second, the frequency of Q variant carriers was not higher in prevalence than among incident cases of ESRD, despite their greater possibility to be affected by selective survival. Finally, a survival advantage associated with the Q variant would be most evident in prevalent case subjects of ESRD and long duration of diabetes, where the cumulative effect of mortality would be highest.
![]() |
ACKNOWLEDGMENTS |
---|
We are indebted to Lucjan Wyrwicz for his assistance with in silico work.
![]() |
FOOTNOTES |
---|
Received for publication 31 July 2001 and accepted in revised form 4 January 2002.
ACR, albumin/creatinine ratio; ESRD, end-stage renal disease; GFR, glomerular filtration rate; IRS, insulin receptor substrate; TDT, transmission disequilibrium test.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|