Skip to main content

Copy number variation at leptin receptor gene locus associated with metabolic traits and the risk of type 2 diabetes mellitus

Abstract

Background

Recent efforts have been made to link complex human traits and disease susceptibility to DNA copy numbers. The leptin receptor (LEPR) has been implicated in obesity and diabetes. Mutations and genetic variations of LEPR gene have been discovered in rodents and humans. However, the association of DNA copy number variations at the LEPR gene locus with human complex diseases has not been reported. In an attempt to study DNA copy number variations associated with metabolic traits and type 2 diabetes mellitus (T2DM), we targeted the LEPR gene locus in DNA copy number analyses.

Results

We identified DNA copy number variations at the LEPR gene locus among a Korean population using genome-wide SNP chip data, and then quantified copy numbers of the E2 DNA sequence in the first two exons overlapped between LEPR and LEPROT genes by the quantitative multiplex PCR of short fluorescent fragment (QMPSF) method. Among the non-diabetic subjects (n = 1,067), lower E2 DNA copy numbers were associated with higher fasting glucose levels in men (p = 1.24 × 10-7) and women (p = 9.45 × 10-5), as well as higher total cholesterol levels in men (p = 9.96 × 10-7). In addition, the significant association between lower E2 DNA copy numbers and lower level of postprandial 2hr insulin was evident only in non-diabetic women, whereas some obesity-related phenotypes and total cholesterol level exhibited significant associations only in non-diabetic men. Logistic regression analysis indicated that lower E2 DNA copy numbers were associated with T2DM (odds ratio, 1.92; 95% CI, 1.26~2.96; p < 0.003) in our nested case-control study. Interestingly, the E2 DNA copy number exhibited a negative correlation with LEPR gene expression, but a positive correlation with LEPROT gene expression.

Conclusions

This work suggests that a structural variation at the LEPR gene locus is functionally associated with complex metabolic traits and the risk of T2DM.

Background

Leptin and the leptin receptor (LEPR) are involved in satiety and energy expenditure via central and peripheral mechanisms. The primary site of leptin action is the hypothalamus where the leptin receptor interacts with the adipocyte-derived leptin signal to regulate appetite, energy balance, and metabolism. LEPRs also regulate energy homeostasis in peripheral tissues including skeletal muscle, liver, pancreas, and adipose tissue. Leptin prevents obesity via LEPRs by stimulating glucose uptake and fatty acid oxidation in skeletal muscle and liver [1–3], and inhibits insulin secretion of pancreatic β-cells [4].

Mutations and genetic variations of the LEPR gene have been discovered in rodents and humans. Sequencing-based detection of LEPR mutations demonstrated that frame-shift or missense mutations of the LEPR gene caused obesity, pituitary dysfunction, hyperphagia, or hypogonadism in humans [5, 6]. In rodents, leptin receptor gene mutations resulted in obesity, hyperglycemia, hyperinsulinemia, and reduced fertility [7, 8]. Common genetic variants (e.g., SNPs) at the LEPR gene locus have been associated with obesity, hyperinsulinemia, type 2 diabetes mellitus (T2DM), and variations in leptin levels in different populations. For example, three non-synonymous SNPs (Arg109Lys, Arg223Gln, and Lys656Asn) have been evaluated for association studies [9–13]. Therefore, most association studies of LEPR polymorphisms have been performed on SNPs or deletion/insertion polymorphisms (DIPs). However, the association of DNA copy number variations (CNVs) at the LEPR gene locus with human diseases has not yet been reported.

In recent years, comprehensive human CNV maps were assembled using various experimental platforms [14–16]. While the growing number of genome-wide association data sets can be utilized to detect clinically-relevant CNVs, absolute copy number information cannot be easily determined by current quantitative assays, with the exception of Fiber-FISH. Nonetheless, recent studies have shown that CNVs are implicated in human diseases including glomerulonephritis (FCGR3B) [17], HIV-1/AIDS (CCL3L1) [18], bipolar disorder and schizophrenia (GLUR7, CACNG2 and AKAP5) [19], neoplasia (14q12) [20], psoriasis (DEFB) [21], and myocardial infarction (C4B) [22].

Results

On employing a candidate gene approach to the association study of CNVs with diabetes and metabolic traits, we targeted the genomic locus of the LEPR gene encompassing approximately 200 kb of chromosome 1. DNA copy number analysis of the Affymetrix 50 K SNP array data showed that the Korean population exhibited CNVs at the LEPR gene locus (Additional file 1). Our QMPSF experiments generated continuous variables of DNA copy number at the LEPR locus for individuals (Additional file 2). Next, the experimental copy number information was transformed to dichotomized DNA copy numbers (low or high) for CNV association studies. An appropriate copy number value (low or high) was then assigned to each individual using either male or female median copy number value among non-diabetic subjects (n = 1,067) as the cutoff. The same cutoff value was used to dichotomize copy numbers among diabetes subjects in the nested case-control study. DNA copy numbers were determined using only a short PCR fragment (herein referred to the "E2 DNA", located near exon 2 of the LEPR gene, see Figure 1) of the entire LEPR gene sequence. To address the question as to whether this short sequence was representative of the whole LEPR gene, we investigated the CNV boundaries using PCR-based (QMPSF) and array-based (Affymetrix SNP array 5.0) analyses. Relative DNA copy numbers were obtained by QMPSF using primers targeting several regions along the LEPR gene sequence. QMPSF results showed that the copy number of the E2 DNA sequence was correlated to that of the promoter-containing exon region at the LEPR locus, but not to the other downstream exonic regions (Figure 1). On the other hand, CNV analysis of the Affymetrix SNP array 5.0 data revealed that the whole LEPR gene locus region was highly CNV-affected, with a varied length of CNV boundaries over the LEPR gene sequence (Additional file 1). This result indicated that the E2 DNA copy number did not represent DNA copy numbers of the whole LEPR gene, but a partial DNA sequence containing the LEPROT (leptin receptor overlapping transcript) gene at the LEPR locus.

Figure 1
figure 1

Schematic diagram of CNVs at the LEPR gene locus. Genomic structures of LEPR and LEPROT genes, including exons, are shown. E2 DNA content relative to the Factor VIII gene from 20 subjects is displayed in the lower panel. Primer target sites (P1 to P7) for QMPSF are indicated along the LEPR sequence. P2* indicates the primer target site for detection of the initial copy number at the LEPR locus. P1 and P2 primer target sites were highly variable for copy number among the 19 individuals tested, but the other primer sites were not.

Next, we tested the statistical significance of the association of the E2 DNA copy number with metabolic traits. In male non-diabetic subjects, significant differences for some demographic characteristics were observed between low (1X) and high (2X) copy number individuals (Table 1). Fasting glucose and total cholesterol levels were higher in individuals with low copy numbers of the E2 DNA at LEPR gene locus in three different age groups of male individuals (Figure 2). The association of low copy numbers with higher levels of obesity-related phenotypes and total cholesterol levels remained significant in linear regression analysis with age-and BMI-adjustments followed by Bonferroni correction for multiple testing among male individuals (Table 1). Among female subjects, E2 DNA copy numbers were also associated with levels of fasting glucose and postprandial 2hr insulin. On the other hand, beneficial phenotypes for obesity-related traits and total cholesterol levels were associated with higher DNA copy numbers (2X) of the LEPR gene locus in male subjects.

Table 1 Association of LEPR CNVs with metabolic traits in nondiabetic subjects.
Figure 2
figure 2

Total cholesterol and fasting glucose levels for individuals with low or high copy numbers of the E2 DNA sequence. Non-T2DM subjects were divided into different age groups (40 s, 50 s, and 60s) and their metabolic phenotypes were compared. (A) Total cholesterol levels. ***: P = 2.13 × 10-6, **: P = 4.33 × 10-3, *: P < 0.05. (B) Fasting glucose levels. ***: P = 1.72 × 10-6, **: P = 8.39 × 10-4, *: P < 0.03.

The copy numbers of the E2 DNA sequence at the LEPR gene locus were used in the nested case-control study of T2DM, in which subjects (n = 137) were matched for age, sex, and BMI at a ratio of approximately 1:2 to control subjects (n = 258). Our case groups showed typical demographic characteristics to diabetes subjects (Additional file 3). For example, diabetic case subjects had higher triglyceride, higher total cholesterol, and lower HDL-cholesterol than those levels in non-diabetic subjects. Logistic regression analysis showed that low E2 DNA copy numbers were associated with T2DM (odds ratio, 1.92; 95% CI, 1.26~2.96), with a significance of P < 0.003 (Table 2).

Table 2 The nested case-control study of T2DM

To study biological roles of the E2 DNA copy numbers, we compared mRNA expression levels of both LEPR and LEPROT genes with E2 DNA copy numbers in lymphoblastoid cell lines. The E2 DNA sequence exists immediately downstream of the overlapped exon 2 of the LEPR and LEPROT genes. Real-time PCR analysis showed a significant negative correlation of the E2 DNA copy number with the LEPR transcript level (R2 = 0.439, p = 1.52 × 10-4), but a weak positive correlation with the LEPROT transcript level (R2 = 0.135, p < 0.05) (Figure 3). Moreover, the LEPR expression level was inversely correlated with the LEPROT expression level (R2 = 0.126, p < 0.05). This result suggests that higher copy numbers of the E2 DNA may contribute to higher transcriptional activity of the LEPROT gene by gene dosage effect, which may be accordingly responsible for downregulation of LEPR gene.

Figure 3
figure 3

Correlation of the E2 DNA copy number with LEPR and LEPROT transcription levels. Expression levels of LEPR and LEPROT relative to GAPDH were compared to the relative E2 DNA content in 32 strains of LCLs. Relative E2 DNA content was calculated from QMPSF data in relation to DNA content at the Factor VIII gene locus in LCLs.

Discussion

Our candidate gene approach to CNV association studies identified the E2 DNA copy number variation at the LEPR gene locus, which was associated with metabolic traits and the risk of T2DM. Lower copy numbers of the E2 DNA were associated with higher fasting glucose in both male and female non-diabetic subjects. However, obesity-related traits and total cholesterol levels exhibited significant association with the E2 DNA copy numbers only in male non-T2DM subjects whereas postprandial 2hr insulin level was significantly higher only in female non-T2DM subjects. These results suggest that the E2 DNA copy number variations at the LEPR gene locus may contribute to glucose homeostasis and energy metabolism in a different manner between male and female subjects.

In our CNV association study, we observed that the beneficial high E2 DNA copy number was associated with lower fasting glucose levels (hypoglycemia) in both genders which are beneficial in terms of diabetes. Similar results were reported for mouse models with LEPR signaling defects [23, 24]. For example, hypothalamus and pancreas-specific LEPR signaling disruption led to fasting hypoglycemia, hyperinsulinemia, and hypertriglyceridemia, highlighting the important role of beta cells in LEPR mediated glucose homeostasis [24]. Hypothalamic LEPR KOVMH mice fed a low fat diet also exhibited hypertriglyceridemia without concomitant hypercholesterolemia [25]. Neuronal LEPR deletion induces diabetic phenotypes such as body weight increase, adiposity, and hyperglycemia in proportion to hypothalamic LEPR deficiency [25]. In addition, antisense-RNA mediated adipocyte-specific reduction of LEPR leads to decrease in LEPR signaling, resulting in increased adipocyte, hypertriglyceridemia, and insulin resistance [26]. As presented in our association study, the high E2 copy number individuals shared similar phenotypes with LEPR deficient animal models. Moreover, human LEPR signaling mutations also caused obesity as well as impairments of pubertal development and growth hormone secretion [5]. Thus, it is conceivable that higher E2 DNA copy numbers may result in a decrease in LEPR signaling, which functions in a complex manner in different tissues.

Indeed, our expression analysis of LEPR gene in LCLs showed that the high E2 DNA copy number was significantly correlated with downregulation of the LEPR gene and concomitant upregulation of the LEPROT gene. Hypothalamic LEPR mRNA expression is indicative of functional LEP-LEPR signaling [27]. The LEPR and LEPROT transcripts share the first two exons, which are not translated in the LEPR gene [28]. Among the LEPR isoforms, the long form (LEPRb) is primarily expressed in specific nuclei of the hypothalamus, whereas the short isoform (LEPRa) is expressed in most tissues. These two isoforms of LEPR have identical extracellular domains but differ in the length of their intracellular domains, which contributes to their different capacity for signal transduction [29]. On the other hand, our CNV analyses of the Affymetrix SNP array 5.0 and QMPSF revealed a varied length of the CNV region at the LEPR gene locus. Also, we found that the E2 DNA copy numbers did not represent DNA copy numbers of the whole LEPR gene but only a partial sequence containing LEPROT gene. Thus, this finding suggests that E2 DNA copy numbers may negatively regulate LEPR gene expression through a gene dosage effect of LEPROT.

However, LEPROT overexpression decreases only cell surface expression of LEPR without modifying the total amount of LEPR [30]. In our study, the observed weak correlation of E2 DNA copy numbers with LEPR mRNA levels in LCLs might be ascribed to inappropriate selection of target tissues for LEPR and LEPROT expression studies. Thus, the correlation of high E2 DNA copy number with LEPROT gene expression may be more significant in an appropriate LEPR target tissues such as the hypothalamus, pancreatic β-cells, or adipocytes. Elucidating tissue-specific regulation of LEPR gene expression could explain the complex pattern of association of phenotypes with the E2 DNA copy number variation. It remains to be investigated whether the E2 copy number confers a gene dosage effect of LEPROT on LEPR signaling using an appropriate target tissue, and whether this occurs via either transcriptional regulation or intracellular transport of LEPR. Alternatively, the E2 DNA sequence may merely play a role as a direct negative transcriptional regulatory sequence for the LEPR gene expression.

Conclusion

A prior CNV analysis of Affymetrix 50 K SNP and 5.0 array chip data indicated the potential presence of CNVs at the LEPR gene locus as well as a varied length of the CNV region. DNA copy number measurements at the LEPR gene locus and subsequent association studies revealed that lower copy numbers of the E2 DNA sequence around LEPROT gene were associated with detrimental phenotypes of metabolic traits (e.g., higher fasting glucose and total cholesterol levels) as well as the risk of T2DM. In addition, the E2 DNA copy number exhibited a negative correlation with LEPR gene expression but a positive correlation with LEPROT gene expression. More research may be warranted to determine if LEPROT expression impacts LEPR expression.

Methods

Study samples

Subjects described in this study include 1,204 individuals from the Korean Genome and Epidemiology Study (KoGES) for CNV detection and subsequent association studies. The KoGES is a community-based prospective cohort study of chronic diseases (e.g., diabetes and hypertension) and their relationship with potential risk factors or lifestyles, which was started in 2001. The study design, sampling, concept and consent have been previously described [31, 32]. Briefly, a total of 10,038 individuals participated in baseline examinations. Approximately 86.6% (n = 8,693) and 67.8% (n = 6,085) of the baseline study subjects participated in the first and the second 2-year follow-up studies, respectively.

For the study of CNV association with metabolic traits, we randomly selected 1,067 non-diabetic subjects who did not develop T2DM during either baseline or two rounds of 2-year follow-up period from the KoGES. WHO criteria were used to define T2DM using an oral glucose tolerance test. Diabetes cases were defined as subjects who were determined to have T2DM at baseline, and then diabetes was diagnosed at least one time during two rounds of 2-year follow-ups. For the nested case-control study, T2DM cases (n = 137) were matched for sex, age±1, and BMI±0.1, with control subjects (n = 258) who did not have diabetes in either baseline or during the 2-year follow-ups. This study was approved by the institutional review boards of Korea Centers for Disease Control and Prevention.

DNA copy number measurement

DNA samples used for DNA copy number measurement were provided by the National Biobank of Korea in Korea Centers for Disease Control and Prevention. We extracted CNV information for candidate genes (diabetes-and obesity-related) from the Affymetrix 50 K SNP chip data in order to estimate whether CNVs were detectable within or near these genes in the Korean population. The Affymetrix 50 K SNP chip data was obtained from the Korean HapMap samples [33]. CNV analysis using CNAT3 (Affymetrix) indicated that the LEPR gene locus contained potential CNV regions (CNVRs) in this Korean population (Additional file 1) [34]. In addition, we selected a copy number invariant region (CN2-2) in chromosome 2 (2q36.1) that exhibited the least copy number variation of the CN2-2 locus among the 90 individuals tested (Additional file 4, 5).

A potential CNV-affected region near rs2025804, upstream of the long isoform of LEPR (LEPRb, NM_002303), was chosen for experimental detection of DNA copy numbers using the QMPSF method [35]. The QMPSF reaction included the Factor VIII gene and CN2-2 as references [34, 36, 37]. The relative copy numbers of a target locus were calculated in relation to the copy numbers of the Factor VIII gene or CN2-2. These two normalizing genes yielded almost same results when used for calculating relative copy numbers. All samples were subjected to three rounds of QMPSF reactions. All reactions were performed in duplicate.

For the justification of further CNV calling, DNA copy numbers of Factor VIII gene were compared between male and female subjects (Additional file 4). Based on the relative copy numbers of Factor VIII gene to CN2-2 region, individuals of lower or higher copy number groups were assigned to be 1X or 2X of DNA copy number for Factor VIII gene. The copy numbers (1X or 2X) of Factor VIII gene was perfectly matched to individuals' sex.

Fluorescent primers were designed to target a specific region of the LEPR gene locus, encompassing over 216 kb of chromosome 1p31, as well as the Factor VIII gene and CN2-2 loci used for normalization of the QMPSF reactions. Following optimization of the PCR reaction, multiplex PCR was performed with the appropriate primers (Additional file 6). The amplicon size for the LEPR locus was 242bp or 245bp depending on the presence of a trinucleotide deletion (-/AGG, rs60086513, at chromosomal position 65603354-65603358) and 199bp for the Factor VIII gene. The DNA copy number of the LEPR locus in relative to the Factor VIII reference was calculated according to the peak height measurements of target and reference loci. The histogram of the LEPR copy numbers indicates two subpopulations among the non-diabetes subjects who had either low or high copy numbers at the LEPR locus (Additional file 2). Thus, median values of the DNA copy number relative to Factor VIII in non-diabetic male and female groups were used as cutoff values to call DNA copy number values in male and female groups, respectively. Finally, either a 'low (1X)' or 'high (2X)' DNA copy number value at the LEPR locus was assigned to all individuals.

Quantitative real-time PCR

Total RNA was isolated from EBV-transformed lymphoblastoid cell lines (LCLs) that were generated from peripheral blood mononuclear cells of participants. cDNA was synthesized as previously described [38] and then used in SYBR green PCR with the appropriate primers (Additional file 6). The difference in threshold cycles (Ct) between GAPDH and the test genes was determined from three experiments, each performed in triplicate. The expression levels of the target genes were determined relative to the expression level of GAPDH.

Statistical analyses

All data are presented as mean±SD. Student's t test was used to compare means of the baseline measurements between two groups. For statistical analysis, DNA copy numbers were dichotomized as '1' for low copy number and '2' for high copy number at the LEPR locus. Multiple regressions were performed for association analyses of metabolic phenotypes in non-diabetic subjects while adjusting for age and sex, which was then corrected by the Bonferroni method for multiple testing. Logistic regression analysis was used for calculating P-values, controlling for age, sex, and BMI as covariates. All analyses were performed with the SPSS 12.0 Statistical Software package.

Abbreviations

CNV:

copy number variation

QMPSF:

quantitative multiplex PCR of short fluorescent fragment

LCL:

lymphoblastoid cell line

LEPR:

leptin receptor

T2DM:

type 2 diabetes mellitus.

References

  1. Aiston S, Agius L: Leptin enhances glycogen storage in hepatocytes by inhibition of phosphorylase and exerts an additive effect with insulin. Diabetes. 1999, 48: 15-20. 10.2337/diabetes.48.1.15.

    Article  CAS  PubMed  Google Scholar 

  2. Wauters M, Considine RV, Van Gaal LF: Human leptin: from an adipocyte hormone to an endocrine mediator. Eur J Endocrinol. 2000, 143: 293-311. 10.1530/eje.0.1430293.

    Article  CAS  PubMed  Google Scholar 

  3. Minokoshi Y, Kim YB, Peroni OD, Fryer LG, Müller C, Carling D, Kahn BB: Leptin stimulates fatty-acid oxidation by activating AMP-activated protein kinase. Nature. 2002, 415: 339-343. 10.1038/415339a.

    Article  CAS  PubMed  Google Scholar 

  4. Seufert J: Leptin effects on pancreatic beta-cell gene expression and function. Diabetes. 2004, 53 (Suppl 1): S152-8. 10.2337/diabetes.53.2007.S152.

    Article  CAS  PubMed  Google Scholar 

  5. Clément K, Vaisse C, Lahlou N, Cabrol S, Pelloux V, Cassuto D, Gourmelen M, Dina C, Chambaz J, Lacorte JM, Basdevant A, Bougnères P, Lebouc Y, Froguel P, Guy-Grand BA: A mutation in the human leptin receptor gene causes obesity and pituitary dysfunction. Nature. 1998, 392: 398-401. 10.1038/32911.

    Article  PubMed  Google Scholar 

  6. Farooqi IS, Wangensteen T, Collins S, Kimber W, Matarese G, Keogh JM, Lank E, Bottomley B, Lopez-Fernandez J, Ferraz-Amaro I, Dattani MT, Ercan O, Myhre AG, Retterstol L, Stanhope R, Edge JA, McKenzie S, Lessan N, Ghodsi M, De Rosa V, Perna F, Fontana S, Barroso I, Undlien DE, O'Rahilly S: Clinical and molecular genetic spectrum of congenital deficiency of the leptin receptor. N Engl J Med. 2007, 356: 237-247. 10.1056/NEJMoa063988.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Chen H, Charlat O, Tartaglia LA, Woolf EA, Weng X, Ellis SJ, Lakey ND, Culpepper J, Moore KJ, Breitbart RE, Duyk GM, Tepper RI, Morgenstern JP: Evidence that the diabetes gene encodes the leptin receptor: identification of a mutation in the leptin receptor gene in db/db mice. Cell. 1996, 84: 491-495. 10.1016/S0092-8674(00)81294-5.

    Article  CAS  PubMed  Google Scholar 

  8. Lee GH, Proenca R, Montez JM, Carroll KM, Darvishzadeh JG, Lee JI, Friedman JM: Abnormal splicing of the leptin receptor in diabetic mice. Nature. 1996, 379: 632-635. 10.1038/379632a0.

    Article  CAS  PubMed  Google Scholar 

  9. Rosmond R, Chagnon YC, Holm G, Chagnon M, Pérusse L, Lindell K, Carlsson B, Bouchard C, Björntorp P: Hypertension in obesity and the leptin receptor gene locus. J Clin Endocrinol Metab. 2000, 85: 3126-3131. 10.1210/jc.85.9.3126.

    CAS  PubMed  Google Scholar 

  10. Takahashi-Yasuno A, Masuzaki H, Miyawaki T, Ogawa Y, Matsuoka N, Hayashi T, Hosoda K, Inoue G, Yoshimasa Y, Nakao K: Leptin receptor polymorphism is associated with serum lipid levels and impairment of cholesterol lowering effect by simvastatin in Japanese men. Diabetes Res Clin Pract. 2003, 62: 169-175. 10.1016/S0168-8227(03)00163-3.

    Article  CAS  PubMed  Google Scholar 

  11. van Rossum CT, Hoebee B, van Baak MA, Mars M, Saris WH, Seidell JC: Genetic variation in the leptin receptor gene, leptin, and weight gain in young Dutch adults. Obes Res. 2003, 11: 377-386. 10.1038/oby.2003.51.

    Article  CAS  PubMed  Google Scholar 

  12. Loos RJ, Rankinen T, Chagnon Y, Tremblay A, Pérusse L, Bouchard C: Polymorphisms in the leptin and leptin receptor genes in relation to resting metabolic rate and respiratory quotient in the Québec Family Study. Int J Obes. 2006, 30: 183-190. 10.1038/sj.ijo.0803127.

    Article  CAS  Google Scholar 

  13. Park KS, Shin HD, Park BL, Cheong HS, Cho YM, Lee HK, Lee JY, Lee JK, Oh B, Kimm K: Polymorphisms in the leptin receptor (LEPR)--putative association with obesity and T2DM. J Hum Genet. 2006, 51: 85-91. 10.1007/s10038-005-0327-8.

    Article  CAS  PubMed  Google Scholar 

  14. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36: 949-951. 10.1038/ng1416.

    Article  CAS  PubMed  Google Scholar 

  15. Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Månér S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science. 2004, 305: 525-528. 10.1126/science.1098918.

    Article  CAS  PubMed  Google Scholar 

  16. Human Genome Structural Variation Working Group, Eichler EE, Nickerson DA, Altshuler D, Bowcock AM, Brooks LD, Carter NP, Church DM, Felsenfeld A, Guyer M, Lee C, Lupski JR, Mullikin JC, Pritchard JK, Sebat J, Sherry ST, Smith D, Valle D, Waterston RH: Completing the map of human genetic variation. Nature. 2007, 447: 161-165. 10.1038/447161a.

    Article  Google Scholar 

  17. Aitman TJ, Dong R, Vyse TJ, Norsworthy PJ, Johnson MD, Smith J, Mangion J, Roberton-Lowe C, Marshall AJ, Petretto E, Hodges MD, Bhangal G, Patel SG, Sheehan-Rooney K, Duda M, Cook PR, Evans DJ, Domin J, Flint J, Boyle JJ, Pusey CD, Cook HT: Copy number polymorphism in Fcgr3 predisposes to glomerulonephritis in rats and humans. Nature. 2006, 439: 851-815. 10.1038/nature04489.

    Article  CAS  PubMed  Google Scholar 

  18. Gonzalez E, Kulkarni H, Bolivar H, Mangano A, Sanchez R, Catano G, Nibbs RJ, Freedman BI, Quinones MP, Bamshad MJ, Murthy KK, Rovin BH, Bradley W, Clark RA, Anderson SA, O'connell RJ, Agan BK, Ahuja SS, Bologna R, Sen L, Dolan MJ, Ahuja SK: The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility. Science. 2005, 307: 1434-1440. 10.1126/science.1101160.

    Article  CAS  PubMed  Google Scholar 

  19. Wilson GM, Flibotte S, Chopra V, Melnyk BL, Honer WG, Holt RA: DNA copy number analysis in bipolar disorder and schizophrenia reveals aberrations in genes involved in glutamate signaling. Hum Mol Genet. 2006, 15: 743-749. 10.1093/hmg/ddi489.

    Article  CAS  PubMed  Google Scholar 

  20. Braude I, Vukovic B, Prasad M, Marrano P, Turley S, Barber D, Zielenska M, Squire JA: Large scale copy number variation (CNV) at 14q12 is associated with the presence of genomic abnormalities in neoplasia. BMC Genomics. 2006, 7: 138-10.1186/1471-2164-7-138.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Hollox EJ, Huffmeier U, Zeeuwen PL, Palla R, Lascorz J, Rodijk-Olthuis D, van de Kerkhof PC, Traupe H, de Jongh G, den Heijer M, Reis A, Armour JA, Schalkwijk J: Psoriasis is associated with increased beta-defensin genomic copy number. Nat Genet. 2008, 40: 23-25. 10.1038/ng.2007.48.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Blaskó B, Kolka R, Thorbjornsdottir P, Sigurdarson ST, Sigurdsson G, Rónai Z, Sasvári-Székely M, Bödvarsson S, Thorgeirsson G, Prohászka Z, Kovács M, Füst G, Arason GJ: Low complement C4B gene copy number predicts short-term mortality after acute myocardial infarction. Int Immunol. 2007, 20: 31-37. 10.1093/intimm/dxm117.

    Article  PubMed  Google Scholar 

  23. Covey SD, Wideman RD, McDonald C, Unniappan S, Huynh F, Asadi A, Speck M, Webber T, Chua SC, Kieffer TJ: The pancreatic beta cell is a key site for mediating the effects of leptin on glucose homeostasis. Cell Metab. 2006, 4: 291-302. 10.1016/j.cmet.2006.09.005.

    Article  CAS  PubMed  Google Scholar 

  24. Bingham NC, Anderson KK, Reuter AL, Stallings NR, Parker KL: Selective loss of leptin receptors in the ventromedial hypothalamic nucleus results in increased adiposity and a metabolic syndrome. Endocrinology. 2008, 149: 2138-2148. 10.1210/en.2007-1200.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. McMinn JE, Liu SM, Liu H, Dragatsis I, Dietrich P, Ludwig T, Boozer CN, Chua SC: Neuronal deletion of Lepr elicits diabesity in mice without affecting cold tolerance or fertility. Am J Physiol Endocrinol Metab. 2005, 289: E403-411. 10.1152/ajpendo.00535.2004.

    Article  CAS  PubMed  Google Scholar 

  26. Huan JN, Li J, Han Y, Chen K, Wu N, Zhao AZ: Adipocyte-selective reduction of the leptin receptors induced by antisense RNA leads to increased adiposity, dyslipidemia, and insulin resistance. J Biol Chem. 2003, 278: 45638-45650. 10.1074/jbc.M304165200.

    Article  CAS  PubMed  Google Scholar 

  27. Baskin DG, Seeley RJ, Kuijper JL, Lok S, Weigle DS, Erickson JC, Palmiter RD, Schwartz MW: Increased expression of mRNA for the long form of the leptin receptor in the hypothalamus is associated with leptin hypersensitivity and fasting. Diabetes. 1998, 47: 538-543. 10.2337/diabetes.47.4.538.

    Article  CAS  PubMed  Google Scholar 

  28. Mercer JG, Moar KM, Hoggard N, Strosberg AD, Froguel P, Bailleul B: B219/OB-R 5'-UTR and leptin receptor gene-related protein gene expression in mouse brain and placenta: tissue-specific leptin receptor promoter activity. J Neuroendocrinol. 2000, 12: 649-655. 10.1046/j.1365-2826.2000.00501.x.

    Article  CAS  PubMed  Google Scholar 

  29. Ghilardi N, Ziegler S, Wiestner A, Stoffel R, Heim MH, Skoda RC: Defective STAT signaling by the leptin receptor in diabetic mice. Proc Nat Acad Sci. 1996, 93: 6231-6235. 10.1073/pnas.93.13.6231.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  30. Couturier C, Sarkis C, Séron K, Belouzard S, Chen P, Lenain A, Corset L, Dam J, Vauthier V, Dubart A, Mallet J, Froguel P, Rouillé Y, Jockers R: Silencing of OB-RGRP in mouse hypothalamic arcuate nucleus increases leptin receptor signaling and prevents diet-induced obesity. Pro Natl Acad Sci. 2007, 104: 19476-19481. 10.1073/pnas.0706671104.

    Article  CAS  Google Scholar 

  31. Shin C, Abbott RD, Lee H, Kim J, Kimm K: Prevalence and correlates of orthostatic hypotension in middle-aged men and women in Korea: the Korean Health and Genome Study. J Hum Hypertens. 2004, 18: 717-723. 10.1038/sj.jhh.1001732.

    Article  CAS  PubMed  Google Scholar 

  32. Lim S, Jang HC, Lee HK, Kimm KC, Park C, Cho NH: The relationship between body fat and C-reactive protein in middle-aged Korean population. Atherosclerosis. 2006, 184: 171-177. 10.1016/j.atherosclerosis.2005.04.003.

    Article  CAS  PubMed  Google Scholar 

  33. Kim KJ, Lee HJ, Park MH, Cha SH, Kim KS, Kim HT, Kimm K, Oh B, Lee JY: SNP identification, linkage disequilibrium, and haplotype analysis for a 200-kb genomic region in a Korean population. Genomics. 2006, 88: 535-540. 10.1016/j.ygeno.2006.03.003.

    Article  CAS  PubMed  Google Scholar 

  34. Jeon JP, Shim SM, Jung J, Nam HY, Oh B, Kimm K, Kim HL, Han BG: A Comprehensive Profile of DNA Copy Number Variations in a Korean Population: Identification of Copy Number Invariant Regions among Koreans. Exp Mol Med. 2009, 41 (9): 618-628. 10.3858/emm.2009.41.9.068.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Vaurs-Barriere C, Bonnet-Dupeyron MN, Combes P, Gauthier-Barichard F, Reveles XT, Schiffmann R, Bertini E, Rodriguez D, Vago P, Armour JA, Saugier-Veber P, Frebourg T, Leach RJ, Boespflug-Tanguy O: Golli-MBP copy number analysis by FISH, QMPSF and MAPH in 195 patients with hypomyelinating leukodystrophies. Ann Hum Genet. 2006, 70: 66-77. 10.1111/j.1529-8817.2005.00208.x.

    Article  CAS  PubMed  Google Scholar 

  36. Levine RL, Wadleigh M, Cools J, Ebert BL, Wernig G, Huntly BJ, Boggon TJ, Wlodarska I, Clark JJ, Moore S, Adelsperger J, Koo S, Lee JC, Gabriel S, Mercher T, D'Andrea A, Frohling S, Dohner K, Marynen P, Vandenberghe P, Mesa RA, Tefferi A, Griffin JD, Eck MJ, Sellers WR, Meyerson M, Golub TR, Lee SJ, Gilliland DG: Activating mutation in the tyrosine kinase JAK2 in polycythemia vera, essential thrombocythemia, and myeloid metaplasia with myelofibrosis. Cancer Cell. 2005, 7: 387-397. 10.1016/j.ccr.2005.03.023.

    Article  CAS  PubMed  Google Scholar 

  37. Jeon JP, Shim SM, Nam HY, Baik SY, Kim JW, Han BG: Copy number increase of 1p36.33 and mitochondrial genome amplification in Epstein-Barr virus-transformed lymphoblastoid cell lines. Cancer Genet Cytogenet. 2007, 173: 122-130. 10.1016/j.cancergencyto.2006.10.010.

    Article  CAS  PubMed  Google Scholar 

  38. Baik SY, Yun HS, Lee HJ, Lee MH, Jung SE, Kim JW, Jeon JP, Shin YK, Rhee HS, Kimm KC, Han BG: Identification of stathmin 1 expression induced by Epstein-Barr virus in human B lymphocytes. Cell Prolif. 2007, 40: 268-281. 10.1111/j.1365-2184.2007.00429.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Dr. Joo-Young Lee for helpful advice during the statistical analysis. DNA samples for this study were provided by National Biobank of Korea. Epidemiological data used in this study were provided from the Korean Genome and Epidemiology Study, KoGES (2001-6111-221, 2002-6111-221). This work was supported by an intramural grant (2007-N-00353-00/2910-212-207) from Korea National Institute of Health, Korea Centers for Disease Control and Prevention.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Jae-Pil Jeon or Bok-Ghee Han.

Additional information

Authors' contributions

JPJ developed the concept and wrote the paper. JPJ and BGH designed the study. SMS generated and analyzed experimental data, and drafted the manuscript. HYN and EJH participated in the QMPSF experiment. GMR extracted and analyzed CNV information from the chip data. HLK and BGH helped with interpretation of the results. BGH supervised the project. All authors read and approved the final manuscript.

Jae-Pil Jeon, Sung-Mi Shim contributed equally to this work.

Electronic supplementary material

12864_2009_3020_MOESM1_ESM.PPT

Additional file 1: Copy number variations at the LEPR locus in a Korean population. LEPR CNVs were analyzed using the Affymetrix 50 K SNP (A, B) and 5.0 SNP arrays (C) for 90 Korean individuals. (A) Affymetrix 50 K SNP array data was obtained from 90 individuals from the Korean HapMap project and then used to extract copy number information using the CNAT3 in reference to the Korean reference genome assembly of 90 individuals. For each SNP probe, standard deviation (STD) of the copy number (CN) values from 90 individuals was calculated and plotted along with the physical positions of SNP probes near the LEPR locus on chromosome 1. The STD of 0.25 was chosen as a cutoff value for CNV calling. A dashed box indicates copy number variations of the LEPR gene. (B) The Affymetrix 50 K SNP array data were analyzed to define CNV regions in reference to the Affymetrix reference genome assembly (provided by Affymetrix Inc.). The P value of < 0.01 was chosen to call CNVs. The CNV regions at the LEPR locus for 90 individuals were schematically depicted near the LEPR locus. A dashed box indicates copy number variations at the LEPR gene. Red bars: copy number gains, blue bars: copy number losses. (C) The whole genomes of 90 individuals were analyzed for CNV detection using the Affymetrix SNP array 5.0, and then compared with one specific reference genome (NA07357). The DNA-Chip Analyzer (dChip) was used for CNV detection http://www.dchip.org. dChip is a Windows software package for probe-level and high-level analysis of gene expression microarrays and SNP microarrays. dChip performed invariant set normalization on the 91 chip data (90 Koreans and one reference). Next it calculated model-based expression values. Copy number variations were inferred by applying median smoothing method. Only CNV regions covering LEPR gene were shown. Assuming that CN states are 2N in the reference genome, CN states indicate the copy numbers of corresponding genomic locus for each individual. (PPT 761 KB)

12864_2009_3020_MOESM2_ESM.PPT

Additional file 2: Histogram of the E2 DNA content among non-T2DM subjects. E2 DNA content was calculated relative to the DNA content at the Factor VIII gene locus from QMPSF data and then plotted to display the distribution pattern of E2 DNA copy numbers among the non-diabetes subjects. The median copy numbers in the men's and women's groups were used to divide them into two subpopulations of lower and higher E2 DNA copy number individuals. The E2 DNA content was then dichotomized into 'lower' or 'higher' copy number groups. (PPT 126 KB)

12864_2009_3020_MOESM3_ESM.XLS

Additional file 3: Characteristics of case-control subjects. *: Student's t-test for difference between low and high CN groups within either cases or controls. **: Student's t-test for difference between T2DM cases and non-T2DM controls. (XLS 34 KB)

12864_2009_3020_MOESM4_ESM.PPT

Additional file 4: Justification of the QMPSF-based CNV calling using Factor VIII and CN2-2 references. A) Experimental validation of the CN2-2 copy number invariant region using QMPSF. CNV analysis of the Affymetrix 50 k SNP array data for 90 individuals identified copy number invariant regions with the least standard deviation of copy number (CN) values in the CNAT3 analysis. Relative copy numbers of the CN2-2 region were calculated in relative to the Factor VIII gene reference from 35 individuals. Error bars indicate standard deviations from three independent QMPSF reactions. Each reaction was performed in duplicate. Relative copy numbers of CN2-2 region to Factor VIII gene exhibited 7.3% of coefficient of variation (CV), suggesting that CN2-2 can be used as a good internal control of further QMPSF reactions. B) Copy number distribution of Factor VIII gene in male and female subjects. Total individuals (n = 1,202) were clearly divided into two subgroups (male and female groups) by relative DNA copy numbers of Factor VIII gene to CN2-2 region. Individuals of lower or higher copy number groups were assigned to be 1X or 2X for Factor VIII gene. The DNA copy numbers (1X or 2X) of Factor VIII gene were perfectly matched to individuals' sex. (PPT 194 KB)

12864_2009_3020_MOESM5_ESM.XLS

Additional file 5: Copy number invariant regions. a: Copy number values (CN) of each SNP probe were obtained from CNV analysis of 90 individuals using the CNAT3. b: Average of copy number values of 90 individuals for each SNP probe. c: Standard deviation of copy number values of 90 individuals for each SNP probe. d: Average copy number values of all SNP probes for 90 individuals within the CNV region. e: Standard deviations of all probes for 90 individuals within the CNV region. (XLS 20 KB)

Additional file 6: Primers used in QMPSF and real-time PCR. (XLS 23 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Jeon, JP., Shim, SM., Nam, HY. et al. Copy number variation at leptin receptor gene locus associated with metabolic traits and the risk of type 2 diabetes mellitus. BMC Genomics 11, 426 (2010). https://doi.org/10.1186/1471-2164-11-426

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-426

Keywords