Diabetes 54:1171-1178, 2005 © 2005 by the American Diabetes Association, Inc. Genetic Variations in the Gene Encoding ELMO1 Are Associated With Susceptibility to Diabetic Nephropathy
1 SNP Research Center, The Institute of Physical and Chemical Research, Yokohama, Kanagawa, Japan
To search for a gene(s) conferring susceptibility to diabetic nephropathy (DN), we genotyped over 80,000 gene-based single nucleotide polymorphisms (SNPs) in Japanese patients and identified that the engulfment and cell motility 1 gene (ELMO1) was a likely candidate for conferring susceptibility to DN, in view of the significant association of an SNP in this gene with the disease (intron 18+9170, GG vs. GA+AA, 2 = 19.9, P = 0.000008; odds ratio 2.67, 95% CI 1.714.16). In situ hybridization (ISH) using the kidney of normal and diabetic mice revealed that ELMO1 expression was weakly detectable mainly in tubular and glomerular epithelial cells in normal mouse kidney and was clearly elevated in the kidney of diabetic mice. Subsequent in vitro analysis revealed that ELMO1 expression was elevated in cells cultured under high glucose conditions (25 mmol/l) compared with cells cultured under normal glucose conditions (5.5 mmol/l). Furthermore, we identified that the expression of extracellular matrix protein genes, such as type 1 collagen and fibronectin, were increased in cells that overexpress ELMO1, whereas the expression of matrix metalloproteinases was decreased. These results indicate that ELMO1 is a novel candidate gene that both confers susceptibility to DN and plays an important role in the development and progression of this disease.
Diabetic nephropathy (DN) is the leading cause of end-stage renal disease in Western countries and Japan (1,2). The pathogenesis of DN appears to be multifactorial. Several genetic and environmental factors are likely to contribute to its development and progression, although the precise mechanisms for this contribution are unknown. Because cumulative epidemiological findings have provided evidence that genetic susceptibility plays an important role in the pathogenesis of this renal disease (3,4), there have been extensive efforts to identify the gene(s) involved in the development and progression of nephropathy, in both type 1 and type 2 diabetes (5,6), but no definitive results have yet emerged. Single nucleotide polymorphisms (SNPs) are the most common genetic variation and are considered to be potentially useful as markers to identify genetic variants that may confer susceptibility to common etiologically complex diseases. After developing a high-throughput system for genotyping SNPs that combines the Invader assay with multiplex PCRs (7), we carried out genome-wide association studies using SNPs to identify loci involved in susceptibility to common diseases. In the study reported here, we performed a genome-wide SNP genotyping analysis of a large panel of Japanese patients with type 2 diabetes in an effort to identify the gene(s) conferring susceptibility to DN. Our data suggest that the engulfment and cell motility 1 gene (ELMO1) is likely to contribute to genetic susceptibility to DN.
DNA preparation and SNP genotyping. DNA samples were obtained from the peripheral blood of patients with type 2 diabetes who regularly attend outpatient clinics at Shiga University of Medical Science, Tokyo Womens Medical University, Juntendo University, Kawasaki Medical School, Iwate Medical University, Toride Kyodo Hospital, Kawai Clinic, Osaka City General Hospital, or Chiba Tokusyukai Hospital. All subjects provided informed consent before enrolling in this study, and DNA extraction was performed according to standard phenol-chloroform procedures. Diabetic patients were divided into two groups according to the following diagnostic criteria: 1) cases of DN, i.e., patients with diabetic retinopathy as well as overt nephropathy, indicated by urinary albumin excretion rates (AERs) 200 µg/min or urinary albumin-to-creatinine ratios (Alb/Cr) 300 mg/gCr, or patients under chronic renal-replacement therapy; and 2) control patients, or patients with diabetic retinopathy but showing no evidence of renal dysfunction, i.e., AER <20 µg/min or Alb/Cr <30 mg/gCr. The SNPs for genotyping were randomly selected from our gene-based Japanese SNP database (available at http://snp.ims.u-tokyo.ac.jp) (8,9). The genotype of each SNP locus was analyzed with Invader assays, as previously described (7). Our first screening involved genotyping 94 DN patients (63 men and 31 women, age 57.9 ± 12.5 years, duration of diabetes 18.6 ± 9.7 years, HbA1c 7.7 ± 1.3%, mean ± SD) and 94 control patients (37 men and 57 women, age 62.7 ± 9.9 years, duration of diabetes 16.2 ± 8.4 years, HbA1c 7.4 ± 1.1%) for 81,315 SNP loci. By evaluating the statistical data using 2 x 3 or 2 x 2 contingency tables, we selected SNPs that showed significant differences in genotypic or allelic frequencies between the DN and control groups. Then we analyzed the SNPs in a larger number of subjects466 DN and 266 control patients (DN patients: 305 men and161 women, age 59.6 ± 13.5 years, duration of diabetes 17.3 ± 10.4 years, HbA1c 7.8 ± 1.1%; control patients: 125 men and 141 women, age 62.9 ± 12.0 years, duration of diabetes 14.6 ± 9.3 years, HbA1c 7.1 ± 1.2%) (10). The ethics committees of the Institute of Physical and Chemical Research and each participating institution approved the study protocol.
Identification of polymorphisms in ELMO1 and genotyping.
ISH.
ISH study. The probes were precipitated and DIG incorporation was assessed by dot blotting. After dewaxing, the tissue sections were dipped into 0.01 mol/l citrate buffer (pH6.0) at 95°C for 30 min, rinsed with 0.1 mol/l PBS (pH 7.4), and permeabilized with 2.5 µg/ml proteinase K (Roche Diagnostics) at 37°C for 10 min. Next, they were briefly washed with 0.1 mol/l PBS, rinsed in 0.1 mol/l triethanolamine (pH 8.0), and then acetylated with 0.1 mol/l triethanolamine containing 0.25% acetic anhydride for 10 min. Hybridization was conducted by the method described previously (13). Briefly, sections were prehybridized and then hybridized with 500 ng/ml DIG-labeled antisense or sense cRNA probe for 17 h at 50°C in a hybridization buffer containing 50% deionized formamide, 1 x Denhardts solution, 10% dextran sulfate, 600 mmol/l NaCl, 0.025% SDS, 5 mmol/l EDTA (pH 8.0), 0.25 mg/ml yeast tRNA, and 10 mmol/l Tris-HCl (pH 7.6). During hybridization, the slides were covered with parafilm and kept in a closed, moist chamber. After hybridization, the samples were rinsed three times with 5 x sodium chloridesodium citrate (SSC) for 10 min at 50°C. The samples were stringently washed with 2 x SSC containing 50% formamid for 20 min at 50°C and thoroughly washed twice with 0.2 x SSC for 20 min at 50°C. Then the samples were immersed in 1.5% blocking reagent dissolved in DIG buffer 1 (100 mmol/l Tris-HCl, pH 7.5, containing 150 mmol/l NaCl) for 60 min at room temperature, preincubated in normal rabbit serum at a dilution of 1:500 in DIG buffer 1 for 30 min, and subsequently incubated in anti-DIG sheep polyclonal antibodies at a dilution of 1:1,000 in DIG buffer 1 for 30 min at room temperature and rinsed with DIG buffer 1. These samples were then incubated in biotinylated anti-sheep rabbit polyclonal antibodies at a dilution of 1:1,000 in DIG buffer 1 for 30 min at room temperature and rinsed again with DIG buffer 1. The samples were next treated with avidin-biotinylated horseradish peroxidase complex solution (Vector Laboratories, Burlingame, CA) at room temperature for 60 min. After immunological incubation, the samples were extensively washed with DIG buffer 1, then processed using 0.1% 3,3'-diaminobenzidine hydrochloride substrate dissolved in 50 mmol/l Tris-HCl (pH 7.4) containing 0.05% H2O2 for 5 min at room temperature in the dark and examined for expression of the specific gene of interest. After a brownish color developed, the reaction was stopped by rinsing in Tris-EDTA buffer (10 mmol/l Tris-HCl, pH 7.6, containing 1 mmol/l EDTA), and the target mRNA signals were visualized. The sections were counterstained with periodic acid-Shiff reagent, dehydrated, and finally mounted in Entellan new (Merck).
Photomicrographic evaluation.
Real-time quantitative RT-PCR.
Statistical analysis.
Genome-wide association study. Among the 81,315 SNP loci we tested in our original screening, 1,615 SNP loci had significant P values <0.01 between DN and control patients (see supplementary Table 1 in the online appendix available at http://diabetes.diabetesjournals.org) and were analyzed further in another larger group of patients (10). Of these SNPs, one SNP locus in the 18th intron of ELMO1 on chromosome 7p14 was strongly associated with DN (GG vs. GA+AA, 2 = 16.3, P = 0.00005; odds ratio [OR] 3.33, 95% CI 1.81 6.15) (Table 1). Correction of multiple testing error was made using the following calculation: Overall P values (P1st x P2nd) x number of tests (1st + 2nd) = 0.001 x 0.00005 x (81,315 x 4 + 1,615 x 4) = 0.017. After performing this correction, we concluded that the association of this landmark SNP with DN was statistically significant.
Subsequent mapping of LD around the landmark SNP in ELMO1 revealed that LD in this region extended 100 kb upstream and 100 kb downstream of the landmark site (Fig. 1). Therefore, the critical region for susceptibility to DN seemed likely to lie within this 200-kb high-LD block. Because ELMO1 is the only gene within this LD block, it was concluded that ELMO1 was the most likely candidate for conferring susceptibility to DN.
We next screened for additional polymorphisms in ELMO1 and identified another 516 polymorphisms (448 SNPs, 49 insertion/deletion polymorphisms, 18 tandem repeat polymorphisms, and 1 other polymorphism) (see supplementary Fig. 1 in the online appendix). Genotyping of these variations using the maximum number of DN (n = 640) and control (n = 426) subjects showed that several of them were associated with DN (Tables 2 and 3; 2/T3information on those SNPs can be found in supplementary Table 2 in the online appendix); in particular, SNP in intron 18 (+9170 A/G) had the strongest association with DN (GG vs. AG+AA, 2 = 19.9, P = 0.000008; OR 2.67, 95% CI 1.714.16). We also analyzed haplotype structure using the expectation maximization algorithm and found that nine SNPs in ELMO1, with an allelic frequency of >0.15, constituted one haplotype block, and that the five common haplotypes accounted for >90% of the population (Fig. 2). Subsequent association studies with DN for each haplotype identified a significant association of haplotypes 2 and 6 with DN. However, the association of these haplotypes was not stronger than that found at the +9,170 A/G locus.
ISH. To identify a possible role of ELMO1 in the pathogenesis of DN, we examined the expression pattern of this gene in normal (db/+m) and diabetic (db/db) mouse kidney by ISH and identified that the expression of ELMO1 could be observed mainly in glomerular epithelial cells and tubular epithelial cells (Fig. 3). The signal of ELMO1 mRNA was remarkably increased in the kidney of diabetic mice (db/db 30 weeks old) (Figs. 3E and F) compared with that of controls (db/+m 30 weeks, Figs. 3B and C), whereas the expression of ELMO2 was not different between these mice (Fig. 3H, I, K, and L). We also examined ELMO1 expression by RT-PCR in various kinds of cultured cells, including glomerular mesangial cells, proximal tubular epithelial cells, COS cells, and retinal pigment epithelial cells, and could only detect definite expression in COS cells.
Real-time quantitative RT-PCR for ELMO1. We next examined the effects of glucose on the expression of ELMO1 in COS cells cultured under either normal (5.5 mmol/l) or high (25 mmol/l) concentrations of glucose. As shown in Fig. 4, the expressions of ELMO1 were significantly increased in the cells cultured under high glucose conditions compared with those in the cells under normal conditions or conditions in which osmolarity had been adjusted with mannitol (+19.5 mmol/l). These results suggested that the increase in ELMO1 expression under high glucose conditions might contribute to the development and progression of DN.
Expression of extracellular matrix genes in ELMO1-overexpressing cells. From the results described above, it was suggested that the excess expression of ELMO1 under high glucose conditions might contribute to the development and progression of DN. Although the precise mechanism of the development and progression of DN is still unknown, it has been suggested that TGF-ß1induced overaccumulation of extracellular matrix proteins is, at least in part, involved in the pathogenesis of DN (17). Therefore, we investigated the effects of excess expression of ELMO1 on the expression of extracellular matrix protein genes, as well as on that of the TGF-ß1 gene, using cells that overexpressed ELMO1. As shown in Fig. 5, the expression of these genes in cells stably transfected with pcDNA3.1-ELMO1 were remarkably increased compared with those in the cells transfected with pcDNA3.1-LacZ (TGF-ß1 1.36-fold vs. 1.5-fold; collagen type 1 [ 1] 5.5-fold vs. 5.8-fold; fibronectin 11.1-fold vs. 12.1-fold; ELMO1-line 1 and ELMO1-line 2 vs. LacZ, respectively) (Figs. 5BD), whereas the expression of MMPs were reduced (MMP2 0.66 vs. 0.27 and MMP3 0.54 vs. 0.57; ELMO1-line 1 and ELMO1-line 2 vs. LacZ, respectively) (Figs. 5E and F).
In a genome-wide case-control association study using SNPs as genetic markers, we identified ELMO1 as a candidate gene conferring susceptibility to DN. Our data also suggest that one intronic SNP (intron 18+9170) might affect susceptibility to DN. The results presented here, as well as in recent publications (10,18,19), provide evidence that a genome-wide case-control association study using gene-based SNPs as genetic markers is a powerful strategy for identifying genes associated with susceptibility to common diseases. Epidemiological findings had strongly suggested a contribution of genetic factors to the development and progression of DN (3,4), but worldwide efforts have so far failed to identify any solid evidence to indicate a genetic susceptibly to the disease. In many cases, the results were conflicting, probably because sample sizes were often inappropriate. However, the major cause of failures to obtain solid conclusions is possibly that multiple genetic factors are involved in DN in a complex manner and the influence of each individual factor is too weak to be identified. Therefore, approaches other than standard candidate-gene analysis or family-based linkage analysis are required to identify genes that are involved in susceptibility to common diseases such as DN. As we began this study by analyzing a huge number of loci (81,315 SNPs) on a genome-wide basis, using a high-throughput genotyping system developed in our institute, and because all of the SNPs contributing to this report are the gene-based SNPs from the Japanese population (8,9), we were able to screen for candidate genes for disease susceptibility more efficiently than would be possible using other SNP databases. The ELMO1 gene, on chromosome 7p14, is a known mammalian homologue of the C. elegans gene, ced-12, which is required for engulfment of dying cells and for cell migration (20). ELMO1 has also been reported to cooperate with CrkII and Dock180, which are homologues of C. elegans ced-2 and ced-5, respectively, to promote phagocytosis and cell shape changes (20,21). However, until now no evidence has been reported to suggest a role for this gene in the pathogenesis of DN. In this study, we showed that the expression of ELMO1 was increased in the kidney of diabetic mice compared with that of control mice, whereas the expression of ELMO2 was not different between diabetic and nondiabetic mice. We also identified that the increased expression of ELMO1 could be observed in the glomeruli isolated from diabetic mice by real-time quantitative PCR (data not shown). These results suggested some role of ELMO1 in the pathogenesis of DN. We next examined the effects of glucose concentrations on the expression of this gene in COS cells, and we found that the expression of ELMO1 was significantly elevated in cells under high glucose conditions compared with cells under normal glucose conditions. The effect of glucose on the expression of ELMO1 was not due to increased osmolality, because the expression was not elevated in cells under normal glucose conditions with 19.5 mmol/l mannitol. From these observations, we speculated that increases in the expression of ELMO1 under high glucose conditions might contribute to the development and progression of DN. Interestingly, the excess expression of ELMO1 resulted in increased expression of extracellular matrix genes and in the decreased expression of MMP genes. Therefore, it is suggested that persistent excess of ELMO1 leads to the overaccumulation of extracellular matrix proteins and to the development and progression of diabetic glomerulosclerosis. Because the increase in the expression of extracellular matrix genes is much greater than that of TGF-ß1, a TGF-ß1independent mechanism was thought to exist. The precise mechanism of ELMO1 influence on the expression of these genes is still unknown. Previous reports suggested that ELMO1 cooperates with Dock180 to function as an activator for Rac-1 (21). Accordingly, we next examined the effects of Rac-1 and coexpression of Dock180 with ELMO1 on the expression of both extracellular matrix genes and MMPs. However, no significant additional effect of these factors on the expression of these genes was detected (data not shown), suggesting that the effects of ELMO1 on the regulation of expression of these genes were independent of Rac-1. The mechanism whereby ELMO1 polymorphisms affect the susceptibility to DN also requires further examination. ELMO1 is >500 kb, and the polymorphisms associated with DN are located far from the transcription start site. Because we could not identify any variations within the coding region or splicing area, we thought that the identified variations might in fact affect ELMO1 gene expression under a high glucose milieu. Indeed, several modulators for transcriptional activity were reported to be located as far as hundreds of kilobasepairs from the transcription initiation site. We then searched for specific protein binding to the DNA corresponding to several candidate regions within ELMO1 and identified high glucoseinduced specific protein binding at the landmark SNP site (data not shown). Because the sequence including this landmark SNP site (intron 18+1970) was compatible with the binding site for GCR-1 (see supplementary Table 2 in the online appendix), which could control the transcriptional activity of glycolytic enzyme genes in response to the alteration of extracellular concentrations of glucose in Saccharomyces cerevisiae (22), this region might be directly involved in conferring susceptibility to DN. However, the elucidation of the precise mechanism requires further study. In summary, we identified ELMO1 as a candidate gene for conferring susceptibility to DN based on results from a genome-wide case-control association study using gene-based SNPs. The evidence presented here suggests that ELMO1 could be a useful target for new drugs to aid in the prevention and treatment of DN.
This work was supported by a grant from the Japanese Millennium Project. The authors would like to thank Shuichi Tsukada, Kana Ohashi, and the technical staff at the SNP Research Center and the Institute of Physical and Chemical Research.
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org. Address correspondence and reprint requests to Shiro Maeda, Laboratory for Diabetic Nephropathy, SNP Research Center, The Institute of PhysicalChemical Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan. E-mail: smaeda{at}src.riken.jp Received for publication November 11, 2004 and accepted in revised form December 20, 2004
Abbreviations: AER, albumin excretion rate; Alb/Cr, albumin-to-creatinine ratio; DIG, digoxigenine; DN, diabetic nephropathy; ISH, in situ hybridization; LD, linkage disequilibrium; MMP, matrix metalloproteinase; SNP, single nucleotide polymorphism; SSC, sodium chloridesodium citrate; TGF, transforming growth factor
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||