Structure and sequence variation at the human leptin receptor gene in lean and obese Pima Indians
Structure and sequence variation at the human leptin receptor gene in lean and obese Pima IndiansD. Bruce Thompson*, Eric Ravussin, Peter H. Bennett and Clifton Bogardus
Phoenix Epidemiology and Clinical Research Branch, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 4212 N 16th Street, Phoenix, AZ 85016, USA
Received November 21, 1996;Revised and Accepted February 3, 1997
The cloning of human and mouse cDNAs from brain that encode high affinity leptin receptors was recently reported. We have physically localized the human leptin receptor gene (LEPR) to a region at 1p31, between the anonymous microsatellite markers D1S515 and D1S198. The genomic structure of the human leptin receptor gene, corresponding to the published human brain cDNA sequence, spans over 70 kb and includes 20 exons. Since the leptin receptor gene is a candidate gene for obesity, and because of its proximity to D1S198, a marker previously linked to insulin secretion, the LEPR gene was sequenced in 20 non-diabetic Pima Indians chosen for extremes in percent body fat and in their acute insulin response to intravenous glucose. Seven polymorphic sites were identified. Two of these polymorphisms, Lys109Arg and Gln223Arg, are amino acid substitutions in the extracellular domain of the leptin receptor, one polymorphism is a silent substitution, and four occur in non-coding regions of the leptin receptor. Four of these sites are in linkage disequilibrium with one another. Nucleotides at three non-coding polymorphic sites were found exclusively in obese Pima Indians. This demonstrates an association between variation at the leptin receptor gene and obesity in humans.
The Pima Indians of Arizona have the highest reported prevalence of non-insulin-dependent diabetes mellitus (NIDDM) in the world, and obesity is also prevalent (1 ,2 ). The Pimas farmed along the banks of the Gila and Salt Rivers in central Arizona for ~2000 years (3 ). In this arid environment, cycles of feast and famine were frequent, depending on the availability of water (3 ). As originally hypothesized by Neel (4 ), natural selection may have favored a metabolically efficient phenotype, based on a `thrifty genotype' that increased the chances of surviving periods when food was scarce. The resulting thrifty phenotype could be characterized by a lower rate of energy expenditure and/or hyperphagia. Either or both of these characteristics would result in positive energy balance and increase fat stores during `times of plenty' which would be advantageous for survival during periods of drought and limited food supplies. However, in westernized societies, where food is plentiful, and typically high fat, this `thrifty genotype' would be disadvantageous, resulting in marked obesity and increased risk of NIDDM.
Two mouse models of obesity, ob and db, are recessively inherited mutations which are characterized by increased rates of insulin secretion, low rates of energy expenditure, hyperphagia, obesity and diabetes (5 ). Parabiosis experiments between control, ob/ob and db/db animals led Coleman to hypothesize that the ob/ob mouse lacked a circulating satiety factor and that food regulatory centers in the brain of the db/db mouse were unresponsive to the same factor (6 ). These genes encoding this circulating satiety factor, now called leptin, and its receptor, were recently cloned, and mutations have been identified that result in absence of leptin production in fat cells in the ob/ob mouse and decreased expression of this receptor in the db/db mouse (7 ). The human leptin (LEP) and leptin receptor genes (LEPR) are therefore good candidate genes potentially resulting in thrifty phenotypes in the Pima Indians. Studies of the leptin gene in the Pima Indians did not reveal any amino acid substitutions (8 ), so we have now examined the leptin receptor gene.
The human leptin receptor gene is located ~1.5 Mb telomeric from D1S198 (9 ,10 ). Based on comparison of the cDNA sequence and the genomic sequence, the human brain leptin receptor gene is encoded in 20 exons and spans over 70 kb of DNA (Fig. 1 ). The first two exons are non-coding and are capable of forming several alternative secondary structures. Exon 1 has a GC content of 70% and includes two 8 bp palindromes beginning at positions 28 and 57 (11 ). Homology between the rodent and human leptin receptor begins in exon 3, which contains the initiation codon and the putative signal sequence that is also encoded in exons 3 and 4. A newly described imperfect dinucleotide repeat is located in the third intron, ~230 nucleotides (nt) downstream of exon 3 (11 ). In a sample of 88 chromosomes genotyped in the Pima, only two alleles were observed with frequencies of 0.49 and 0.51. The single putative transmembrane region lies in exon 18 and the previously described splice site for the alternate forms of the leptin receptor occurs at the exon/intron boundary of exon 19 (13 ). The intracellular domain is encoded in exons 19 and 20. Exon 20, the largest exon, spans over 900 nt and encodes the last 274 amino acids of the leptin receptor (Fig. 1 ). The intron/exon boundaries conform to the established consensus sequences and all the splice donor and acceptor sites contain GT or AG dinucleotides, respectively (14 ).
Twenty Pima Indian volunteers who are participating in longitudinal studies of NIDDM were admitted to the Clinical Diabetes and Nutrition Section, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health (26 ,27 ). Written informed consent was obtained from all subjects who were then given physical examinations. Subjects were healthy, taking no medications and were non-diabetic. DNA for sequencing was extracted from Epstein-Barr virus transformed lymphocytes, or from blood samples drawn in conjunction with other studies (28 ). Body fat was determined by hydrostatic underwater weighing and corrected for age and sex (29 ). AIR (the log10 of the mean plasma insulin concentration above basal, over 3, 4, and 5 min after the injection of 25 g of intravenous glucose) was determined in all individuals (30 ). The subjects were divided into four groups based on percent body fat and acute insulin response. Percent body fat +- SD averaged 40 +- 5 and 22 +- 8 in the obese and lean groups, respectively. The AIR averaged 2.7 +- 0.2 and 1.7 +- 0.1 [mu]U/ml in the high and low insulin secretion groups.
Oligonucleotide primers were designed from the cDNA sequence and used to sequence across the intron/exon boundaries of the 20 exons of the leptin receptor gene. YAC clones containing the leptin receptor gene (9 ) were obtained from Research Genetics (Huntsville, AL). XL PCR (Perkin Elmer, Norwalk, CA) was performed according to manufacturers' instructions. Alu-vector PCR and Inverse PCR were also used to isolate the ends of YAC arms (26 ). In some instances extra long PCR was combined with inverse PCR to isolate large segments of the leptin receptor. PCR products were amplified in parallel from Pima genomic DNA and the yeast artificial chromosomes 940H12, 930F4 and 895B12 containing the LEPR gene. All PCR products were sequenced to assure they were from the leptin receptor. Statistical analyses were performed using programs of the SAS institute (Cary, NC).
Oligonucleotide primers for sequence analysis were designed on the basis of the intron sequence, synthesized on an Applied Biosystems oligonucleotide synthesizer (model 391, Applied Biosystems, Foster, CA), and used to PCR amplify the corresponding sequence from Pima genomic DNA. DNA sequence analysis was carried out with an ABI 377 automated sequencer (Perkin Elmer, Norwalk, CA). PCR templates for sequence analysis were separated and isolated from 0.8% agarose gels and purified using QIAquick gel extraction kits (Qiagen, Chatsworth, CA). The individual exons and intron sequence flanking each exon have been submitted to GenBank under the following accession numbers: exon 1, U59246; exon 2, U59247; exon 3, U59248; exon 4, U59249; exon 5, U59250; exon 6, U59251; exon 7, U59252; exon 8, U59253; exons 9 and 10, U59254; exon 11, U59255; exon 12, U59256; exons 13 and 14, U59257; exon 15, U59258; exon 16, U59259; exon 17, U59260; exon 18, U59261; exon 19, U59262; exon 20, U59263. The following oligonucleotides were used to amplify regions of the leptin receptor:
We wish to acknowledge the technical and analytical assistance of Warren Apel, Victoria Ossowski, Jeff Sutherland, Jennifer Biesterfeildt, Larry Grahmn and Rachel Janssen. We thank the clinical staff of the Clinical Diabetes and Nutrition Section for collection of the phenotypic information and Dr Michal Prochazka for helpful comments and discussion. Finally, we wish to thank the members of the Gila River Indian Community, without whose assistance and dedication this work would not be possible.
1 Knowler W.C., Bennett P.H., Hamman, R.F. and M. Miller. (1978) Diabetes incidence and prevalence in Pima Indians: a 19-fold greater incidence than in Rochester, Minnesota. Am. J. Epidemiol.108, 497-505.
2 Knowler, W.C., Pettitt, D. J., Saad, M. F., Charles, M. A., Nelson, R. G., Howard, B. V., Bogardus, C. and P. H. Bennett. (1991) Obesity in the Pima Indians: its magnitude and relationship with diabetes. Am. J. Clin. Nutr.53, 1543S-1551S.MEDLINE Abstract
3 Fontana, B.L. (1979) Handbook of North American Indians. Smithsonian Institution. Volume 10, pages 125-135. Fink, T.M. The prehistoric Irrigation canals and reservoirs of southern Arizona and their possible impact on Hohokam health. Health and Disease in the Prehistoric Southwest. Edited by Charles F. Mervs and Rover J. Miller, Arizona State University Anthropological Research Papers, number 23, 359-379.
4 Neel. J.V. (1962) Diabetes mellitus: a `thrifty' genotype rendered detrimental by progress. Am. J. Hum. Genet.14, 353-362.
5 Coleman, D.L. (1978) Obese and diabetes: Two mutant genes causing diabetes-obesity syndromes in mice. Diabetologia14, 141-148.MEDLINE Abstract
6 Coleman, D.L. (1973) Effects of Parabiosis of obese with diabetes and normal mice. Diabetologia9, 294-298.MEDLINE Abstract
7 Chau, S.C., Chung, W.K., Wu-Peng, S., Zhang, Y., Liu, S., Tartaglia, L. and R. L. Leibel. (1996) Phenotypes of mouse diabetes and rat fatty due to mutations in the OB (leptin) receptor. Science271, 994-996.
8 Maffei, M., Stoffel, M., Barone, M., Dammerman, M., Ravussin, E., Bogardus, C., Ludwig, D.S., Flier, J.S., Talley, M., Auerbach, S. and J.M. Friedman. (1996) Absence of mutations in the human OB gene in obese/diabetic subjects. Diabetes45, 679-682.MEDLINE Abstract
9 Thompson. D. B., Sutherland, J., Appel, W. and V. Ossowski. (1997) A physical map at 1p31 encompassing the acute insulin response locus and the leptin receptor. Genomics39, 227-230.
10 Winick, J.D., Stoffel, M. and J.F. Friedman. (1996) Identification of microsatellites linked to the human leptin receptor gene on chromosome 1. Genomics36, 221-222.MEDLINE Abstract
11 Thompson, D.B., Ossowski, V., Sutherland, J., Apel, W. and Biesterfeldt, J. (1996) Human leptin receptor (LEPR) gene, exons 1-20. GenBank Accession numbers U59246-U59263.
12 Kozak, M. (1986) Influences of mRNA secondary structure on initiation by eukaryotic ribosomes. Proc. Natl. Acad. Sci. USA83, 2850-2854.MEDLINE Abstract
13 Lee, G.H., Proenca, R., Montez, J.M., Carroll, K.M., Darvishzadeh, J.G. and J.M. Friedman. (1996) Abnormal splicing of the leptin receptor in diabetic mice. Nature379, 632-635.MEDLINE Abstract
14 R.A. Padgett, P.J. Grabowski, M.M. Konarska, S. Seiler and P.A. Sharp. (1986) Splicing of messenger RNA precursors. Annu. Rev. Biochem.55, 1119-1150.
15 Yallow, R.S. and W.A. Bauman. (1990) Plasma insulin in health and disease. In Rifkin, H. and Porter, D. Jr (eds), Diabetes mellitus: Theory and Practice. New York, Elsevier p 119-142.
16 Tartaglia, L.A., Dembski, M., Weng, X, Deng, N, Culpepper, J., Devos, R., Richards, G.J., Campfield, L.A., Clark, F.T., Deeds, J., Muir, C., Sanker, S., Moriarty, A., Moore, K.J., Smutko, J.S., Mays, G.G., Woolf, E.A., Monroe, C.A. and Tepper, R.A. (1995) Identification and expression cloning of a leptin receptor, OB-R. Cell83, 1263-1271.MEDLINE Abstract
17 Cioffi, J.A., Shafer, A.W., Zupancic, T.J., Smith-Gbur, J., Mikhail, A., Platika, D. and Snodgrass, H.R. (1996) Novel B219/OB receptor isoforms: Possible role of leptin in hematopoiesis and reproduction. Nature Med.2, 585-589.
18 Chau, S.C., White, D.W., Wu-Peng, X.S., Lui, S., Okada, N., Kershaw, E.E., Chung, W.K., Power-Kehoe, L., Chau, M., Tartagila, L.A. and Leibel, R.L. (1996) Phenotype of fatty due to Gln269Pro mutation in the leptin receptor (Lepr). Diabetes45, 1141-1143.
19 Kozak, M. (1989) Circumstances and mechanisms of inhibition of translation by secondary structure in eucaryotic mRNAs. Mol. Cell. Biol.9, 5134-5142.MEDLINE Abstract
20 Kozak, M. (1987) An analysis of 5'-noncoding sequences from 699 vertebrate messenger RNAs. Nucleic Acids Res.15, 8125-8148.MEDLINE Abstract
21 Considine, R.V., Considine, E.L., Williams, C.J. and Caro, J.F. (1996) The hypothalamic leptin receptor in humans-identification of incidental sequence polymorphisms and absence of the db/db mouse and fa/fa rat mutations. Diabetes45, 992-994.MEDLINE Abstract
22 Thompson, D.B., Janssen, R.C., Ossowski, V.M., Prochazka, M., Knowler, W.C. and Bogardus, C. (1995) Evidence for linkage between a region on chromosome 1 and the acute insulin response in Pima Indians. Diabetes44, 478-481.MEDLINE Abstract
23 Hummel, K.P., Dickie, M.M. and Coleman, D. (1996) Diabetes, a new mutation in the mouse, Science153, 1127-1128.
25 Norman, R.A., Leibel, R.L., Devoto, M., Knowler, W.C., Thompson, D.B., Bogardus, C. and Ravussin, E. (1996) Linkage of obesity and energy metabolism to markers flanking homologues of rodent obesity genes is not evident in Pima Indians. Diabetes45, 1299-1232.
26 Lillioja, S., Mott, D.M., Spraul, M., Ferraro, R., Foley, J.E., Ravussin, E., Knowler, W.C., Bennett, P.H. and Bogardus, C. (1993) Insulin resistance and insulin secretory dysfunction as precursors of non-insulin-dependent diabetes mellitus. N. Engl. J. Med.329, 1988-1992.MEDLINE Abstract
27 Knowler, W.C., Pettitt, D.J., Saad, M.F. and Bennett, P.H. (1990) Diabetes mellitus in the Pima Indians: incidence, risk factors and pathogenesis. Diabetes Metab. Rev.6, 1-27.MEDLINE Abstract
28 DeFronzo, R.A., Bonadonna, R.C. and Ferrannini, E. (1992) Pathogenesis of NIDDM: a balanced overview. Diabetes Care15, 318-68.MEDLINE Abstract
29 Siri, W.E. (1961) Body composition from fluid spaces and density: analysis of methods. pp.223-244. Brozek, J and Henschel, A. In Techniques for Measuring Body Composition: Proceeding of a Conference. National Research Council. Washington DC, USA.
30 Janssen, R., Bogardus, C., Takeda, J., Knowler, W.C. and Thompson, D.B. (1994) Linkage analysis of acute insulin secretion with GLUT2 and glucokinase in Pima Indians and the identification of a missense mutation in GLUT2. Diabetes43, 558-563.MEDLINE Abstract
31 Zoghbi, H. and Chinault, C. (1994) Generation of YAC contigs by walking. In D.L. Nelson and B.H. Brownstein(eds) YAC Libraries a Users Guide. W.H. Freeman and Company, New York pp 93-112.
*To whom correspondence should be addressed. Tel: +1 602 200 5300; Fax: +1 602 200 5335; Email: bthompson@phx.niddk.nih.gov
-->
This page is maintained by OUP admin. Last updated Fri Apr 11 08:44:24 BST 1997. Part of the OUP Journals World Wide Web service.
Copyright
Oxford University Press, 1996