Leptin gene promoter DNA methylation in WNIN obese mutant rats
© Kalashikam et al.; licensee BioMed Central Ltd. 2014
Received: 4 October 2013
Accepted: 29 January 2014
Published: 5 February 2014
Obesity has become an epidemic in worldwide population. Leptin gene defect could be one of the causes for obesity. Two mutant obese rats WNIN/Ob and WNIN/GROb, isolated at National Centre for Laboratory Animal Sciences (NCLAS), Hyderabad, India, were found to be leptin resistant. The present study aims to understand the regulatory mechanisms underlying the resistance by promoter DNA methylation of leptin gene in these mutant obese rats.
Male obese mutant homozygous, carrier and heterozygous rats of WNIN/Ob and WNIN/GROb strain of 6 months old were studied to check the leptin gene expression (RT-PCR) and promoter DNA methylation (MassARRAY Compact system, SEQUENOM) of leptin gene by invivo and insilico approach.
Homozygous WNIN/Ob and WNIN/GROb showed significantly higher leptin gene expression compared to carrier and lean counterparts. Leptin gene promoter DNA sequence region was analyzed ranging from transcription start site (TSS) to-550 bp length and found four CpGs in this sequence among them only three CpG loci (-309, -481, -502) were methylated in these WNIN mutant rat phenotypes.
The increased percentage of methylation in WNIN mutant lean and carrier phenotypes is positively correlated with transcription levels. Thus genetic variation may have effect on methylation percentages and subsequently on the regulation of leptin gene expression which may lead to obesity in these obese mutant rat strains.
KeywordsLeptin expression Promoter DNA methylation Insilico analysis Obesity WNIN
Obesity has become an alarming epidemic to global health with increased hypertension, coronary atherosclerosis, myocardial hypertrophy, type II diabetes and cardiovascular morbidity & mortality. Although obesity is caused by both genetic and environmental factors the molecular mechanism underlying obesity is not yet well established. It is certainly a multifaceted disease with complex interactions between genetic traits and environmental factors, regulated by epigenetic mechanisms.
Earlier many experiments have been conducted to study obesity and the underlying mechanisms involved in various animal models especially in Ob/Ob and db/db mouse strains, Zucker and Koletsky rat strains. Two mutant obese rat strains were isolated from WNIN stock, designated as WNIN/Ob and WNIN/GROb in 1994, at National Institute of Nutrition, Hyderabad, India. These mutant rats are characterized by hyperphagia, hyperinsulinemia, hypertriglyceridemia, higher body fat% (47%) and are found to be leptin resistant as well. Additionally WNIN/GROb rat show impaired glucose tolerance (IGT) too. The inheritance of obesity in these rats follows a co-dominant pattern as per Mendelian laws and the phenotypic ratio of the offspring is 1:2:1 (1 lean+/+ : 2 carrier+/- : 1 Obese-/-). Both the homozygous obese male and females are infertile while the carriers are fertile and are used for breeding and maintaining the stock.
Interest on the pathophysiology of obesity has been intensified recently with the discovery of leptin. Leptin, a 167 amino acid protein and the product of the obese gene Leptin (LEP), is a key hormone primarily secreted by adipocytes. Leptin levels increase exponentially with increasing fat mass[4, 5]. Leptin inhibits food intake and regulates energy intake and expenditure by central action on the hypothalamus. Activation of central leptin receptors increases the activity of the sympathetic nervous system[6, 7] which stimulates energy expenditure in adipose tissue[8, 9]. Leptin levels significantly increase in obese humans and this endogenous hyperleptinemia does not reduce appetite or increase energy expenditure suggesting that they are leptin resistant. Recent studies demonstrate that leptin promoter methylation regulates leptin expression. Melzner et al reported that CpGs in the proximal leptin promoter were highly methylated in pre-adipocytes and gets highly demethylated as they mature to terminally differentiated adipocytes. This suggests that the methylation of CpGs inhibits leptin expression, while their demethylation activates the expression. Although WNIN/Ob rats and WNIN/GROb were found to be leptin resistant like other rodent models the regulatory mechanisms underlying this resistance is not yet known. Thus the present study aims to understand the leptin resistance by determining the promoter DNA methylation pattern of leptin gene which presumably should explain the extent of leptin expression seen in these animals.
The experimental procedure was approved by the ethics committee on animal experiments at the National Institute of Nutrition (NIN), Hyderabad, India and was performed in accordance with the principles of laboratory animal care. Male obese WNIN rats i.e.homozygous, WNIN/Ob obese (-/-), heterozygous WNIN/Ob carrier (+/-), homozygous WNIN/Ob lean (+/+) and WNIN/GROb obese (-/-), WNIN/GROb carrier (+/-), WNIN/GROb lean (+/+) of 6 months old (3 per strain) were obtained from National Centre for Laboratory Animal Sciences (NCLAS), NIN, Hyderabad, India. They were housed individually in polypropylene cages with wire mesh bottom cases maintained at 22 ± 2°C, under standard lighting conditions (12-h light/dark cycle). Animals were fed ad libitum with stock rat chow and free access to deionized water. The animals were sacrificed and their retroperitoneal adipose tissue was collected and snap frozen in liquid nitrogen and stored at -80°C until analysis.
Leptin mRNA expression in retroperitoneal adipose tissue by semi-quantitative PCR
Total RNA was isolated from 100 mg of adipose tissue (retroperitoneal) using Qiazol reagent. This was followed by the synthesis of cDNA from 2 μg of total RNA using Invitrogen kit (Invitrogen Life technologies, Carlsbad, CA). Primers were designed with the aid of primer quest software (Integrated DNA Technologies, Corolville, IOWA). Semi-quantitative PCR was conducted to analyze the expression of leptin (5′ACTTCATTCCCGGGCTTC3′; 5′GGTCTCGCAGGTTCTCCA3′), with the internal standard18SrRNA (5′CCAGAGCGAAAGCATTTGCCAAGA3′;5′AATCAACGCAAGCTTATGACCCGC3′). The amplified products were resolved on 1.2% agarose gel electrophoresis and the image was quantified with the Bio-Rad gel documentation system using Quantity One software (Bio-Rad Laboratories, Herculus, CA). Results were expressed as the ratio of the intensities of the band of the gene of interest to that of the 18s rRNA.
DNA was isolated from 100 mg of retroperitoneal adipose according to the manufacturer’s instructions of QIAGEN DNeasy tissue kit (Qiagen India Pvt. Ltd., India). The isolated gDNA was quantified using a ND-1000 spectrophotometer (Nano Drop Technologies, Inc).
Sodium bisulfite modification
Bisulfite modified gDNA was prepared using the EZ DNA Methylation-Gold kit (Zymo Research, catalog #D5005) according to the manufacturer’s instructions. The bisulfate reaction was performed on 0.5 μg gDNA and the reaction volume was adjusted to 20 μl with sterile water and 130 μl of CT conversion reagent was added. The sample tubes were placed in a thermal cycler (Bio-Rad) and the following steps were performed: 10 min at 98°C, 2 h 30 min at 64°C and stored at 4°C. The converted samples were added into Zymo-Spin ICTM Column containing 600 μl of the M-Binding Buffer and mixed by inverting the column several times. After mixing, the column was centrifuged at full speed for 30 sec and discarded the flow-through. The column was then washed by adding 200 μl of M-Desulphonation buffer and incubated at room temperature (20-30°C) for 15-20 min. After incubation, the column was centrifuged at full speed for 30 sec. The column was again washed by adding 200 μl of M-Wash buffer and spun at full speed (this step was repeated). The converted gDNA was eluted by adding 20 μl of M-Elution buffer into the column and spin. The bisulfite converted gDNA samples were stored at -20°C until methylation analysis. Pyro Q-CpG Software (Biotage) was used to check automatically the generation of dispensation order (sequence) and it also performs quality control of each sample for completion of bisulfite treatment, eliminating the need for manual inspection and estimation levels of non-converted DNA.
Quantitative methylation analysis
Quantitative methylation analysis of the leptin gene promoter was performed with the MassARRAY Compact system (Sequenom, San Diego, CA, USA). The target regions of the sodium bisulfate treated gDNA were amplified by PCR with a set of primers: leptin_10 left primer 5′TTATGGATTAGTAATGAAGATTTTAATAGA3′; leptin_10 right primer 5′TCAATAAAATACCCACCCTTAAAAA3′. Each forward primer contained a 10-nucleotide tag (5′-AGGAAGAGAG-3′) to balance the PCR, and each reverse primer contained a T7 promoter tag and sequence insert (5′-CAGTAATACGACTCACTATAGGGAGAAGGCT-3′) for in vitro transcription. The amplification protocol comprised an initial incubation at 94°C for 4 min; 45 cycles of denaturation at 95°C for 20 s, annealing at 62°C for 30 s, and extension at 72°C for 1 min; and a final incubation at 72°C for 3 min. Unincorporated deoxynucleoside triphosphates were dephosphorylated by the addition of 2 μl of Premix for in vitro transcription including 0.3 U of shrimp alkaline phosphatase (Sequenom, San Diego, CA, USA). The reaction mixture was incubated at 37°C for 40 min, after which the phosphatase was inactivated by incubation for 5 min at 85°C. A portion of the PCR products (2 μl) was then subjected to in vitro transcription, with RNase A cleavage being used for the T-reverse reaction (Sequenom). Samples were spotted onto a 384-pad Spectro-CHIP (Sequenom) with the use of a MassARRAYnanodispenser (Sequenom) and spectra were acquired with a MassARRAY compact matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometer (Sequenom). The resulting methylation calls were analyzed with EpiTyper software (Sequenom) to generate quantitative results for each CpG site or an aggregate of multiple CpG sites.
Insilico analysis of leptin promoter
The leptin gene promoter DNA sequence of rat was obtained from the UCSC genome browser (http://genome.ucsc.edu/). In this study the promoter was defined as the one which can extend from transcription starting site (TSS) to -550bp upstream. The CpG island searcher (http://cpgislands.usc.edu) was used to screen for CpG islands (D.T and P.A.J). The methylator software (http://bio.dfci.harvard.edu/Methylator) was used to identify CpG and predict whether CpG in a promoter DNA sequence are likely to be methylated or not.
All values were represented as mean ± S.E. Data was analyzed by one-way analysis of variance (ANOVA) followed by the multiple range test or least significant difference method. Wherever heterogeneity of variance was observed, differences between groups were tested using non-parametric Mann-Whitney U test. The differences were considered significant at p < 0.05.
Insilico analysis of leptin gene promoter DNA sequence
The leptin promoter DNA sequence in silico analysis - The CpG locus and predicted methylated CpG sites
Promoter length (bp)
No. of CpG
Position of CpG
No. of methylated CpG
Methylated CpG position
TSS to -100
-101 to -200
-201 to -300
-301 to -400
-401 to -500
-501 to -550
Leptin gene expression in adipose tissue
Leptin gene promoter DNA methylation
The fundamental feature of the epigenetic markers is their dynamic flux, which is the basic mechanism of environmental adaptability. Methylation of CpG sites have been widely observed in vertebrates and occasionally reported in invertebrates. Regulatory motifs on promoters can be identified by virtue of their preferential location relative to the TSS. The changes in methylation status of the CpG islands overlapping gene promoter enable epigenetic control of transcription. At the same time, established methylation pattern has to be maintained with high fidelity since aberrant methylation events (hypo- and hyper-methylation) lead to dysregulation of the downstream gene. Methylation of DNA essentially leads to a repression of transcription by interfering with the binding sequence of transcription factors and through the binding of methyl-CpG binding proteins.
A comprehensive insilico analysis of epigenetic properties of the rat leptin gene promoter DNA sequence necessitates experimental testing that affords direct measurement of the methylation state of this promoter region. In the present study, 550bp length leptin gene promoter DNA sequence was analysed and identified four CpGs which were in close proximity to the TSS and had the potentials to be methylated. The present in silico analysis is in line with a study which noted that methylation surrounding the TSS is tightly linked to transcriptional silencing while methylation of more downstream regions is unassociated with the magnitude of gene expression[16, 17]. Considering that the comparative genomics approach is widely used to identify functional sequence elements and regulatory networks this analysis was performed as a prelude to the experimental testing of the leptin gene promoter methylation in WNIN mutant rat genomic samples which may help to design the experiments, to minimize possible experimental bias, and thus improve validity of the experimental data.
Leptin, a multifunctional hormone, is predominantly delivered by mature adipocytes. The present study demonstrated that the leptin mRNA expression levels were varying with their adipose levels in WNIN mutant rat strains i.e. WNIN/Ob & WNIN/GROb and their carrier and lean counterparts. The varying adipose tissue levels indicate that the genes associated with adipogenesis regulation may be different in these mutant WNIN rat phenotypes. These observations are in agreement with few reports which showed strong correlation between the levels of leptin mRNA expression and adipocyte size supporting the importance of leptin in altering the fat mass[19, 20]. Another study on human subjects showed that the developmental increase in leptin mRNA expression in adipose tissue during childhood, reaching maximal capacity in adulthood. In addition a study also depicted an evidence linking leptin to a direct regulation of adipose tissue metabolism through inhibition of lipogenesis and stimulation of lipolysis. Thus these results support that the differentialexpression of leptin mRNA in these mutant phenotypes could be due to the genetic differences among the strains studied although these animals were fed with the same rat chow. According to our knowledge this is the first report on leptin mRNA expression in genetically mutant rat models originated from WNIN that follow the Mendelian genetic ratio’s (1:2:1). However the underlying mechanism for the close correlation between leptin mRNA level and adipocyte mass is not yet understood in these WNIN obese mutant rats.
Many researchers have suggested that the expression of genes controlling energy homeostasis could be regulated by epigenetic mechanisms, which may play a role in the development of obesity[23, 24]. Considering this, the regulation of leptin gene expression in these rat models was studied by epigenetic approach i.e. by analysing promoter DNA methylation in silico and in vivo and correlated with gene expression. It was noticed that the leptin gene promoter DNA methylation percentages were varying in mutant WNIN obese rat phenotypes. The decreased percentage of methylation in leptin gene promoter in WNIN obese than in lean and carriers are in accordance with a recent study which reported genetic polymorphisms associated with differential methylation in the H19/IGF2 locus. Another report indicates that genetic variation may have a substantial impact on local methylation patterns[26, 27], but neither the extent to which methylation is affected by genetic variation, nor the mechanisms are yet clear. Our results also demonstrated that methylated CpGs are in close proximity to the TSS in the rat promoter studied. This work is in agreement with Melzner et al who showed that the three specific CpG sites within the minimal LEP promoter region played an important role in regulating reporter gene expression in human adipose cell line.
The present study clearly demonstrated that the body weights, adipocyte mass and leptin gene expression are in line with leptin promoter DNA methylation percentages in WNIN obese mutant rats suggesting that the genetic differences can have profound influence on methylation patterns. The present in silico methylation prediction was in agreement with in vivo methylation status. The sole exception was at the locus -481 which showed methylation in in vivo, while the in silico, analysis of methylation could not predict. In silico analysis also showed the presence of CpG at the locus -309 and its predicted methylation, which was identified proximal to the TSS in the rat promoter DNA sequence. Interestingly, our in vivo studies also demonstrated clearly that methylation was more at locus -309 which regulated leptin gene expression in WNIN mutant obese rat phenotypes. These observations are in line with several studies which showed that methylation of CpG sites that are close to the TSS in promoter DNA sequence will have more effect on regulation of gene expression. Melzner et al. also reported that the CpGs in the proximal leptin promoter were highly methylated in pre-adipocytes and during maturation toward terminally differentiated adipocytes, this promoter region was found to be highly demethylated. Hence it was clear that the CpG sites which are close to TSS methylated may have profound effect on gene expression.
The present study clearly demonstrated that the leptin gene expression is in line with leptin promoter DNA methylation patterns in WNIN obese mutant rats suggesting that genetic variation seemed to have an effect on methylation percentages and subsequently on the regulation of gene expression. These results may have implications for the functional interpretation of mechanisms underlying association of genetic variants with the development of obesity.
Transcription start site
Wistar National Institute of Nutrition.
This work was supported by National Institute of Nutrition, ICMR, Government of India, New Delhi, India. We also acknowledge Meher Vani Chaduvula for her inputs and editing of the manuscript. Technical support was provided by Aceprobe Technologies (India) Pvt Ltd to analyze promoter methylation on Sequenom Epityper.
- Giridharan NVHNSM: A new rat model for the study of obesity. Scand J Lab Anim Sci. 1996, 23: 131-139.
- Zhang Y, Proenca R, Maffei M, Barone M, Leopold L, Friedman JM: Positional cloning of the mouse obese gene and its human homologue. Nature. 1994, 372: 425-432. 10.1038/372425a0View ArticlePubMed
- Mantzoros CS: The role of leptin in human obesity and disease: a review of current evidence. Ann Intern Med. 1999, 130: 671-680. 10.7326/0003-4819-130-8-199904200-00014View ArticlePubMed
- Lonnqvist F, Arner P, Nordfors L, Schalling M: Overexpression of the obese (ob) gene in adipose tissue of human obese subjects. Nat Med. 1995, 1: 950-953. 10.1038/nm0995-950View ArticlePubMed
- Considine RV, Sinha MK, Heiman ML, Kriauciunas A, Stephens TW, Nyce MR, Ohannesian JP, Marco CC, McKee LJ, Bauer TL: Serum immunoreactive-leptin concentrations in normal-weight and obese humans. N Engl J Med. 1996, 334: 292-295. 10.1056/NEJM199602013340503View ArticlePubMed
- Haynes WG, Morgan DA, Walsh SA, Mark AL, Sivitz WI: Receptor-mediated regional sympathetic nerve activation by leptin. J Clin Invest. 1997, 100: 270-278. 10.1172/JCI119532PubMed CentralView ArticlePubMed
- Scarpace PJ, Matheny M, Pollock BH, Tumer N: Leptin increases uncoupling protein expression and energy expenditure. Am J Physiol. 1997, 273: E226-E230.PubMed
- Commins SP, Marsh DJ, Thomas SA, Watson PM, Padgett MA, Palmiter R, Gettys TW: Norepinephrine is required for leptin effects on gene expression in brown and white adipose tissue. Endocrinology. 1999, 140: 4772-4778.PubMed
- Mizuno TM, Kleopoulos SP, Bergen HT, Roberts JL, Priest CA, Mobbs CV: Hypothalamic pro-opiomelanocortin mRNA is reduced by fasting and [corrected] in ob/ob and db/db mice, but is stimulated by leptin. Diabetes. 1998, 47: 294-297. 10.2337/diab.47.2.294View ArticlePubMed
- Stoger R: In vivo methylation patterns of the leptin promoter in human and mouse. Epigenetics. 2006, 1: 155-162. 10.4161/epi.1.4.3400View ArticlePubMed
- Melzner I, Scott V, Dorsch K, Fischer P, Wabitsch M, Brüderlein S, Hasel C, Möller P: Leptin gene expression in human pre-adipocytes is switched on by maturation-induced demethylation of distinct CpGs in its proximal promoter. J Biol Chem. 2002, 277: 45420-45427. 10.1074/jbc.M208511200View ArticlePubMed
- Bird AP: CpG-rich islands and the function of DNA methylation. Nature. 1986, 321: 209-213. 10.1038/321209a0View ArticlePubMed
- Cooper DN, Youssoufian H: The CpG dinucleotide and human genetic disease. Hum Genet. 1988, 78: 151-155. 10.1007/BF00278187View ArticlePubMed
- Galm O, Herman JG, Baylin SB: The fundamental role of epigenetics in hematopoietic malignancies. Blood Rev. 2006, 20: 1-13. 10.1016/j.blre.2005.01.006View ArticlePubMed
- Wade PA: Methyl CpG-binding proteins and transcriptional repression. Bioessays. 2001, 23: 1131-1137. 10.1002/bies.10008View ArticlePubMed
- Hellman A, Chess A: Gene body-specific methylation on the active X chromosome. Science. 2007, 315: 1141-1143. 10.1126/science.1136352View ArticlePubMed
- Hsieh CL: Dependence of transcriptional repression on CpG methylation density. Mol Cell Biol. 1994, 14: 5487-5494.PubMed CentralView ArticlePubMed
- Gregory SG, Barlow KF, McLay KE, Kaul R, Swarbreck D: The DNA sequence and biological annotation of human chromosome 1. Nature. 2006, 441: 315-321. 10.1038/nature04727View ArticlePubMed
- Zhang Y, Guo KY, Diaz PA, Heo M, Leibel RL: Determinants of leptin gene expression in fat depots of lean mice. Am J PhysiolRegulIntegr Comp Physiol. 2002, 282: R226-R234.View Article
- Guo KY, Halo P, Leibel RL, Zhang Y: Effects of obesity on the relationship of leptin mRNA expression and adipocyte size in anatomically distinct fat depots in mice. Am J PhysiolRegulIntegr Comp Physiol. 2004, 287: R112-R119. 10.1152/ajpregu.00028.2004.View Article
- Schoof E, Stuppy A, Harig F, Carbon R, Horbach T, Stöhr W, Rascher W, Dötsch J: Comparison of leptin gene expression in different adipose tissues in children and adults. Eur J Endocrinol. 2004, 150: 579-584. 10.1530/eje.0.1500579View ArticlePubMed
- Yang WS, Lee WJ, Funahashi T, Tanaka S, Matsuzawa Y, Chao CL, Chen CL, Tai TY, Chuang LM: Weight reduction increases plasma levels of an adipose-derived anti-inflammatory protein, adiponectin. J Clin Endocrinol Metab. 2001, 86: 3815-3819. 10.1210/jcem.86.8.7741View ArticlePubMed
- Campion J, Milagro FI, Martinez JA: Individuality and epigenetics in obesity. Obes Rev. 2009, 10: 383-392. 10.1111/j.1467-789X.2009.00595.xView ArticlePubMed
- Lomba A, Milagro FI, Garcia-Diaz DF, Marti A, Campion J, Martínez JA: Obesity induced by a pair-fed high fat sucrose diet: methylation and expression pattern of genes related to energy homeostasis. Lipids Health Dis. 2010, 9 (9): 60-PubMed CentralView ArticlePubMed
- Heijmans BT, Kremer D, Tobi EW, Boomsma DI, Slagboom PE: Heritable rather than age-related environmental and stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus. Hum Mol Genet. 2007, 16: 547-554. 10.1093/hmg/ddm010View ArticlePubMed
- Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai SL, Arepalli S, Dillman A, Rafferty IP, Troncoso J, Johnson R, Zielke HR, Ferrucci L, Longo DL, Cookson MR, Singleton AB: Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010, 6 (5): e1000952-13, 10.1371/journal.pgen.1000952PubMed CentralView ArticlePubMed
- Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, Sommer IE, Kahn RS, Ophoff RA: The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS One. 2009, 4 (8): e6767-26, 10.1371/journal.pone.0006767PubMed CentralView ArticlePubMed
- Harishankar N, Kumar PU, Sesikeran B, Giridharan N: Obesity associated pathophysiological & histological changes in WNIN obese mutant rats. Indian J Med Res. 2011, 134: 330-340.PubMed CentralPubMed
- Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, Degner JF, Gilad Y, Pritchard JK: DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol. 2011, 12: R10- 10.1186/gb-2011-12-1-r10PubMed CentralView ArticlePubMed
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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.