Diabetes 56:3020-3026, 2007 DOI: 10.2337/db07-0338 © 2007 by the American Diabetes Association
LPIN2 Is Associated With Type 2 Diabetes, Glucose Metabolism, and Body Composition
1 Genetic Epidemiology Unit, Department of Epidemiology and Biostatistics, Clinical Genetics, and Forensic Molecular Biology, Erasmus Medical Centre, Rotterdam, the Netherlands Address correspondence and reprint requests to Yurii S. Aulchenko, Department of Epidemiology and Biostatistics, Erasmus Medical Centre Rotterdam, Postbus 2040, 3000 CA Rotterdam, Netherlands. E-mail: i.aoultchenko{at}erasmusmc.nl
Abbreviations:
DEXA, dual-energy X-ray absorptiometry; ERF, Erasmus Rucphen Family; ISI, insulin sensitivity index; OGTT, oral glucose tolerance test; PPAR, peroxisome proliferator–activated receptor; SNP, single nucleotide polymorphism; UTR, untranslated region; WHR, waist-to-hip ratio
OBJECTIVE—To identify the type 2 diabetes gene located at chromosome 18p11. RESEARCH DESIGN AND METHODS—We investigated the region in a young genetically isolated population by genotyping 34 single nucleotide polymorphisms (SNPs) in 78 case subjects and 101 control subjects. Two SNPs were selected and followed up in two cohorts. The first cohort came from a general Dutch population. In this cohort, association with type 2 diabetes was investigated using 616 type 2 diabetic case subjects and 2,890 control subjects; association with oral glucose tolerance test data was performed in 361 normoglycemic people. Association with fat distribution was studied in the second replication cohort, consisting of 836 people from the genetically isolated population. RESULTS—At the initial step, we found that the common C allele of SNP rs3745012 was associated with type 2 diabetes (odds ratio 2.01, P = 0.03). This SNP is located at the 3' untranslated region of the LPIN2 gene, which is a plausible candidate for type 2 diabetes and obesity. In the cohort from the general Dutch population, we demonstrated that rs3745012 interacts with BMI in determination of type 2 diabetes: whereas in subjects with high BMI, the common C allele is associated with type 2 diabetes, the same allele exhibits a neutral or protective effect in lean subjects (P = 0.05 overall effect, P = 0.02 interaction). Most remarkably, rs3745012 strongly affected composite insulin sensitivity index (P = 0.006 for overall effect, P = 0.004 for interaction). In the second replication cohort, we found that the allele C of rs3745012 increases trunk-to-legs fat mass ratio (P = 0.001) and may affect other fat-related measurements. CONCLUSIONS—rs3745012 SNP of the LPIN2 gene is associated with type 2 diabetes and fat distribution. The majority of type 2 diabetic patients are obese, and obesity by itself may be a cause of insulin resistance (1,2). An important determinant of insulin resistance is not only total fat mass, but also fat distribution (3,4). Given the same degree of obesity, insulin sensitivity varies largely depending on visceral adipose tissue content (5,6). Abdominal obesity measures and visceral fat content may be better predictors of insulin resistance when compared with BMI (7). Also, abnormally decreased adipose tissue mass, such as that observed in lipodystrophy, leads to severe insulin resistance in humans (8,9) and animal models (10–13).
Several known genes for type 2 diabetes are also implicated in lipodistrophy. Examples include peroxisome proliferator–activated receptor We reported a genome-wide scan using a combined association and linkage test in a sample of type 2 diabetic patients derived from a recently genetically isolated population from the Southwest of the Netherlands (25). Evidence for association with type 2 diabetes was found for one region at chromosome 18p11 (P = 0.001 at marker D18S63). The locus at 18p was originally found to be linked to type 2 diabetes in families from Finland and Southern Sweden (26) and was confirmed in a second series derived from the Southwest of the Netherlands (27). In all three studies, the effect of the 18p region was strongest in the obese subpopulation. One of the human homologs of mouse Lpin1 gene, LPIN2, is located in the 18p11 region. In this study, we examine the association between type 2 diabetes and related traits and LPIN2 gene in three cohorts.
Study population. Two series of patients and control subjects were studied. The first was derived from a young genetically isolated population located in the Southwest of the Netherlands. The characteristics of the population regarding linkage disequilibrium and drift are described elsewhere (28–30). The original sample was formed by 79 small families (117 patients) that could be traced to a common ancestor within 13 generations. The diagnosis of type 2 diabetes was based on American Diabetes Association criteria (31). The details of the recruitment procedure and analysis are described elsewhere (25). The study protocol was approved by the medical ethics committee of the Erasmus Medical Centre Rotterdam, and written consent was obtained from all participants. In the present study, we used 78 patients who were probands in our previous research and who were not or were remotely (>10 meioses) related. The unrelated spouses of the type 2 diabetic patients, together with the spouses of patients from other ongoing studies, were used as control subjects (n = 101). This schema justifies the use of the sample as independent case and control subjects. No significant differences in sex ratio or age between patients and control subjects were observed (P > 0.1). The BMI was significantly higher in patients (30.21) than in control subjects (27.53, P = 0.03). The second series was derived from the Rotterdam Study. This is a large cohort from outbred population that is aimed to study prevalence, incidence, and determinants of chronic disease in the elderly (32). The medical ethics committee of the Erasmus Medical Centre Rotterdam approved the study protocol. Baseline examinations, including a detailed interview, physical examination, and blood sampling, were conducted between 1990 and 1993. In the Rotterdam Study, 3,506 people were genotyped. Diagnosis of type 2 diabetes was made based on the use of antidiabetic medication or random glucose level >11.1 mmol/l, according to the previous work of Rietveld et al. (33). Of the investigated group, 616 were affected. The effect of sex was not significant. Type 2 diabetic people were on average older (73.46 vs. 68.09 years in control subjects, P < 0.001) and had higher BMI (26.76 vs. 26.14 kg/m2 in control subjects, P < 0.01) and waist-to-hip ratio (WHR) (0.93 vs. 0.90 in control subjects, P < 0.001). Within the Rotterdam Study, oral glucose tolerance test (OGTT) was performed for 361 normoglycemic people, thus making fasting glucose and insulin and postload glucose and insulin at 30 and 120 min available.
To study body fat distribution, we used data on 836 diabetes-free individuals from the family-based Erasmus Rucphen Family (ERF) study. The study protocol was approved by the medical ethics board of Erasmus Medical Centre Rotterdam. The ERF study consists of
Measurements.
Glucose level was measured in millimoles per liter and insulin in milliunits per liter. As indicator of insulin resistance we used HOMA-IR, as suggested by Matthews et al. (36):
Genotyping. In the initial scan, all single nucleotide polymorphisms (SNPs) were genotyped using the Pyrosequencing technology according to the manufacturer's instructions (Biotage, Uppsala, Sweden). For each assay, a standard PCR was performed in a total volume of 10 µl using 20 ng of genomic DNA, 10 pmol of each primer, one of which was biotinylated at the 5' end, and AmpiTaq Gold DNA polymerase (Applied Biosystems, Foster City, CA). A full list of primer sequences can be found in Supplementary Table 1 (available in an online appendix at http://dx.doi.org/10.2337/db07-0338). The biotinylated PCR product was bound to magnetic streptavidin-coated Dynabeads (Dynal, Oslo, Norway) by shaking at 65°C for 15 min and then denatured in 0.5 mol/l NaOH. DNA bound to the Dynabeads was transferred between solutions using the magnetic PSQ Sample Tool (Biotage). After washing to remove excess NaOH, the single-stranded template was hybridized to 15 nmol sequencing primer by incubating at 80°C for 2 min and cooling for 15 min. The pyrosequencing reaction was then carried out in a PSQ96 instrument (Biotage) using the PSQ SNP Reagent kit (Biotage). Genotypes were analyzed using the PSQ96 Analysis Software (Biotage). In the ERF and Rotterdam Study, genotyping was performed using Taqman allelic discrimination Assay-By-Design (Applied Biosystems) (38). Forward and reverse primer sequences and the minor groove binding probes are presented in Supplementary Table 2. We used the reverse-strand design for SNP9 and SNP10. The assays used 5 ng genomic DNA and 5-µl reaction volumes. The amplification and extension protocol was as follows: an initial activation step of 10 min at 95°C preceded 40 cycles of denaturation at 95°C for 15 s and annealing and extension at 50°C for 60 s. Allele-specific fluorescence was then analyzed on an ABI Prism 7900HT Sequence Detection System with SDS v. 2.1 (Applied Biosystems). Based on the analysis of 5% of blind duplicates, there was a 100% concordance in genotyping this polymorphism.
Expression analysis. To study expression of LPIN2, whole-blood samples were collected from seven participants of the Rotterdam Study. Ten milliliters of peripheral venous blood was drawn from each donor into heparin-containing Vacutainer tubes (BD Biosciences, Breda, the Netherlands). Selective lysis of erythrocytes was achieved by the addition of three volumes of hypotonic buffer (150 mmol/l ammonium chloride, 10 mmol/l potassium bicarbonate, and 1 mmol/l EDTA disodium salt) to whole blood and incubation on ice for 20 min. Nucleated blood cells were then collected with centrifugation at 500g for 10 min at 4°C and washed twice with ice-cold PBS before RNA extraction. Total RNA was extracted from leukocytes with RNABee reagent (Tel-Test, Friendswood, TX) according to the manufacturer's protocol. RNA quality was determined with an A260-to-A280 ratio and capillary electrophoresis on an Agilent 2100 Bioanalyzer automated analysis system (Agilent Technologies, Palo Alto, CA). Hybridization of the RNA to Affymetrix GeneChip Human Genome U133 Plus 2.0 array (Affymetrix, Santa Clara, CA) was performed by the microarray core facility of the Erasmus MC Center for Biomics according to Affymetrix specifications. Background subtraction, probe signal summarization, and normalization were calculated according to the RMA algorithm using the R Bioconductor (41). LPIN2 gene is represented by three probe sets on Affymetrix U133 plus2 microarray: probe sets 202459_s_at and 202460_s_at target different regions of the 3' UTR; the probe set 244799_s_at targets the exon IV of LPIN2 transcript. The position of Affymetrix target sequences for LPIN2 gene on genomic alignment is shown in Supplementary Fig. 1.
Statistical analysis.
In the family-based ERF study, we assessed association using a variance-components model, as implemented in the SOLAR software (42). The null hypothesis of no association was formulated using the polygenic model with covariates. The hypothesis of association was obtained by adding SNP (coded as –1 for AA, 0 for AB, and 1 for BB) as a covariate. To test significance of association, the likelihood ratio test was performed: under the null, twice the difference between log-likelihood of these nesting models is distributed as
Fine-mapping of 18p11 in a young genetically isolated population. With 26 SNPs, we covered the telomeric region 18p11.22–18p11.32 spanning from 0.2 to 12.3 Mb. Table 1 shows SNPs typed, their position, counts of subjects with different genotypes, and ORs with P values. In total, 78 type 2 diabetic patients and 101 control subjects were genotyped. The average heterozygosity of the SNP markers was 38% with a minimum of 6% and maximum of 50%. The average coverage was one SNP per 465 kb, with a maximal gap of 1.3 Mb (Table 1 and Supplementary Fig. 2A and B).
In our previous analysis, the strongest association was found for D18S63 (position 3,428,540 in the Ensembl database) and D18S1105 (position 1,823,900), flanking the LPIN2 gene region. In the current SNP analysis, association was found for two SNPs, SNP9 and SNP12 (Table 1). For SNP9, after adjustment for sex and age, the common allele C was associated with odds ratio (OR) 2.01 (95% CI 1.09–3.71, P = 0.03); the frequency of C was 0.84 in case subjects and 0.73 in control subjects. For SNP12, the G allele was associated with type 2 diabetes (OR = 1.66, P = 0.03, after adjustment for age and sex). The region showing association with type 2 diabetes (SNP9–SNP12) contains six hypothetical and seven known genes: EMILIN2, LPIN2, MYOM1, MRCL3, MRLC2, TGIF, and DLGAP1. From these seven, LPIN2 was studied further because of the association of its 3' UTR SNP9 with type 2 diabetes and the analogy to mouse Lpin1, the gene involved in insulin resistance in the fdl mouse (13,43). According to the HapMap data, LPIN2 is included into a large block of linkage disequilibrium; the boundaries of the block coincide well with the gene boundaries (Supplementary Fig. 3). We typed eight additional SNPs in the LPIN2 gene (Table 1; SNPs with study name starting with "LPN;" Supplementary Fig. 2B). However, no strong association with additional SNPs was found.
Role of LPIN2 in the general population.
We studied the association of SNP9 with early type 2 diabetes pathology in 361 normoglycemic subjects who underwent OGTT (descriptives shown in Supplementary Table 3). For fasting glucose, insulin, and insulin resistance (HOMA-IR), the results, although insignificant, point in the same direction as those of the analysis of type 2 diabetes: the common SNP9C allele is a risk allele at high BMI, whereas in those with low BMI, it appears to decrease fasting glucose, insulin, and insulin resistance (Fig. 2).
We found a significant effect of SNP9 on glucose (P = 0.04) and insulin (P = 0.02), when studying the area under the curve, and on the insulin sensitivity index (ISI) (P = 0.006) (Fig. 2). For all three measures, interaction with BMI was significant. Consistent with our previous findings on type 2 diabetes, the common SNP9C allele was associated with increased risk of early pathology in participants with high BMI, whereas the effect was reversed in those with low BMI (Fig. 2).
SNP9 and body composition. We further studied association between LPIN2 and body composition in 836 diabetes-free participants of the family-based cohort derived from the EFR study. Table 2 shows studied traits and results of analysis. The family structure was adjusted for within variance components framework. Similar to the findings in the Rotterdam Study, when adjusting for age and sex, the SNP9C was associated with an increased BMI, albeit only marginally significantly (effect of 0.51 kg/m2, P = 0.08). In the ERF study, SNP9C increased both waist (by 1.67 cm, P = 0.03) and hip circumference (by 1.03 cm, P = 0.06).
In the analysis of DEXA data, we found that SNP9C increased total fat mass by 1.1 kg (P = 0.05). This happened almost exclusively through increase in trunk fat mass (by 0.9 kg, P = 0.01). A strong positive effect of SNP9C was observed on various adiposity indexes: fat mass index (P = 0.02), ratio of total fat to total tissue (P = 0.01), ratio of trunk fat to total tissue (P = 0.001), and ratio of trunk to legs fat mass (P = 0.001) (Table 2). After adjustment for trunk-to-legs fat mass ratio, the association of SNP9C with waist circumstance, fat mass, trunk fat mass, fat mass index, ratio of total fat to total tissue, and ratio of trunk fat to total tissue became insignificant, suggesting that association to these variables can be explained by the trunk-to-legs fat ratio. The same effect was observed when adjusting for the trunk fat–to–total tissue ratio. These analyses indicate that trunk-to-leg fat (or trunk fat–to–total tissue) mass ratio may be the most important endophenotype associated with SNP9 and type 2 diabetes.
Potential functionality of SNP9. We studied expression of LPIN2 in seven healthy individuals from the Rotterdam Study. Four of these were homozygous for the SNP9C, and three were C/T heterozygotes. Using exact Wilcoxon's nonparametric test, we showed that the probe sets targeting the 3' UTR did not differ in expression between T carriers and noncarriers (P > 0.2). However, the probe set 244799_s_at, located in exon IV, did show marginally significant (P = 0.06) downregulation in the carriers of the T allele.
In our previous study in a young genetically isolated population, we found evidence for a type 2 diabetes locus at chromosome 18p11. This region was also previously reported by Parker et al. (26) in a Scandinavian study and by van Tilburg et al. (27) in the general Dutch population. Both studies suggested that the effect of the locus is modified by BMI. In the fine-typing of the 18p11 region with 34 SNPs, the association was found to be confined to the LPIN2 gene. In a small series of case and control subjects, nominally significant association was found between type 2 diabetes and SNP9 (rs3745012), located in the 3' UTR of the LPIN2 gene. We extended the analysis of association between LPIN2 and type 2 diabetes in an independent series of patients from a large population-based outbred cohort, the Rotterdam Study. In this outbred population, we observed interaction between SNP9 and BMI. Consistent with our findings in the genetically isolated population, in which patients had a high BMI, SNP9C allele was associated with type 2 diabetes in obese participants. In the Rotterdam Study, SNP9 was also associated with BMI and early type 2 diabetes pathology measured using OGTT. Lastly, in a family-based study conducted in the genetically isolated population, we demonstrated that SNP9C is associated with increased adiposity indexes, most notably trunk fat. Together, these studies provide evidence that the SNP9C allele may be a thrifty allele (47,48), i.e., it is common and increases obesity and risk of type 2 diabetes in obese subjects. The issue of multiple testing is worth consideration when judging the results of our study. For example, in the first stage, we tested 34 SNPs and found two to be associated to type 2 diabetes with nominal P < 0.05. When corrected for multiple testing, these results would have become insignificant. However, in our study, we did not rely on experiment-wise significance achieved at any stage. As suggested by others (49), we attempted to replicate our finding: in stage one we selected SNP9 and followed it up in two other cohorts, achieving positive results. Such replication procedure strongly reduces the chance of a false-positive finding (49). The LPIN2 gene is one of the three human homologs of the mouse Lpin1 gene, which is responsible for the lipodystrophy in the fld line (13), which, among other traits, is characterized by severe insulin resistance. The product of Lpin1 is required upstream of PPARG for normal adipocyte differentiation (24). The effects of increased Lpin1 expression on insulin resistance are opposite, depending on the tissue. Whereas Lpin1 overexpression in muscle produces obesity-associated insulin resistance, enhanced Lpin1 expression in mature adipocytes improves insulin sensitivity (43). A recent study of human homolog of mouse Lpin1 gene, LPIN1, has shown its potential involvement in glucose metabolism (23).
LPIN2, being a member of the same lipin family as the LPIN1 and bearing notable sequence similarity to both mouse Lpin1 and human LPIN1, is therefore a plausible candidate gene for type 2 diabetes and body composition. Human LPIN2 contains 20 exons and spans The most important question is how SNP9 may be related to the risk of type 2 diabetes. SNP9 is a noncoding polymorphism located in the 3' UTR of the gene at the third nucleotide downstream of the stop codon (Supplementary Fig. 4). Strong computational evidences suggested that the region containing SNP9 is a potential target site for miRNA miR-371. However, our experiments on coexpression of luciferase reporter and miRNA do not reveal detectable interaction between miR-371 and the LPIN2 3' UTR. We therefore have to conclude that, in line with recent observations (46), perfect seed matching is not sufficient to confer regulation in case of miR-371. Alternatively, our heterologous expression system may be not sensitive enough to detect possible subtle effects of miR-371 on expression of LPIN2. In our limited expression experiments, we observed increased levels of LPIN2 expression in people homozygous for the "risk" C alleles. Although these results reach only marginal significance (P = 0.06), we may hypothesize that SNP9 polymorphism (or other LPIN2 polymorphism in linkage disequilibrium with SNP9) leads to altered expression of the LPIN2 gene. This, by analogy with the mouse model, may lead to (possibly obesity-induced) type 2 diabetes. Further investigation, including larger expression and protein studies are required to establish the link between LPIN2, obesity, and type 2 diabetes.
This research was supported by the Netherlands Organisation for Scientific Research and the Centre for Medical Systems Biology. Genotyping for the isolated population was conducted in the laboratories of Oxagen; the work of C. Vince is gratefully acknowledged. We thank C. Zillikens, F. Rivadeneira, A. Isaacs, and F. Sayed-Tabatabaei for valuable discussions. Above all, we thank the participants from the Genetic Research in Isolated Populations, the ERF study, and the Rotterdam Study who made this work possible.
Published ahead of print at http://diabetes.diabetesjournals.org on 5 September 2007. DOI: 10.2337/db07-0338. I.M. is currently affiliated with the National Institute of Agricultural Botany, Cambridge, U.K. Additional information for this article can be found in an online appendix at http://dx.doi.org/10.2337/db07-0338. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Received for publication March 11, 2007 and accepted in revised form August 26, 2007
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||