Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10−16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.
Introduction
Leptin is an adipocyte-derived hormone that helps maintain homeostatic control of fat tissue mass by signaling the status of body energy stores to the appetite-regulating circuits of the brain (1). Rare homozygous mutations in the leptin (LEP) gene can cause complete leptin deficiency that results in hyperphagia and severe early-onset obesity, which can be treated effectively by exogenous leptin administration (2,3). Mice and patients heterozygous for these mutations show partial leptin deficiency and increased body weight (4–6).
In the general population, leptin concentrations correlate closely with body fat mass. However, there is wide interindividual variability; ∼10–20% of obese individuals have leptin concentrations that are similar to those observed in nonobese individuals, which is in part due to genetic differences (7,8). Twin and family studies suggest that 30–50% of variation in leptin at any given level of adiposity and across different ethnic groups is explained by genetic differences (8). The implications of this variability for body weight regulation remain poorly understood.
Identification of genetic variants associated with circulating leptin may shed new light on the role of variability in leptin levels in the general population. In a recent genome-wide association study (GWAS) of leptin concentrations, we identified four loci associated with leptin concentrations independent of BMI (9). The variant most strongly associated with leptin concentrations was rs10487505, located 21 kb upstream from LEP, in a region shown to harbor a long noncoding RNA (EST EL947753) that influences the transcriptional control of leptin expression (10). The leptin-decreasing allele of rs10487505 was nominally associated with ∼0.03 kg/m2 higher BMI in adults and 1.05-fold increased risk of early-onset obesity (9). More recently, the association of the leptin-decreasing allele of rs10487505 with higher adult BMI, body fat percentage, and risk of extreme obesity was replicated in the UK Biobank (10). The most pronounced association, however, was observed for body size at 10 years of age; carriers of the leptin-decreasing allele reported being “plumper” at age 10 years (compared with peers) more frequently than carriers of the allele associated with higher leptin concentration. The association between rs10487505 and childhood body size was recently replicated in 14,521 Norwegian children, and the peak effect of rs10487505 on BMI was observed in 1.5-year-old children (11).
In the current study, we sought to elucidate the genetic basis of leptin concentrations through screening genetic variants with an exome-targeted array in up to 57,232 individuals of European, African, East Asian, or Hispanic ancestry. We confirm five previously established and identify five novel variants associated with leptin concentrations, including four missense variants in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The novel LEP variant, Val94Met (rs17151919), is associated with leptin concentrations in adults of African ancestry only. The leptin-lowering Met94 allele of the rs17151919 variant is associated with higher BMI in young children but shows a weak or no association with BMI in adulthood, suggesting that leptin regulates early adiposity.
Research Design and Methods
Study Design
We performed an exome-based association study using data from 35 cohorts comprising up to 57,232 adults (≥18 years old), of whom 50,321 were of European descent, 4,387 of African descent, 2,036 of East Asian descent, and 488 of Hispanic descent. We carried out additional analyses in men and women separately. All analyses were performed in models combining studies of all ancestries and in European ancestry cohorts only, for both additive and recessive genetic models. All participating institutions and coordinating centers approved the project. Informed consent was obtained from all study participants. We report the study-specific design, sample quality control, and descriptive statistics in Supplementary Tables 1 and 2.
Outcome Traits
The participating studies acquired residuals for leptin concentrations (in ng/mL) using linear regression, adjusting for age, genome-wide principal components, and any study-specific covariates (e.g., study center). The residuals were calculated with and without adjustment for BMI. Studies with unrelated individuals acquired the residuals in men and women separately, whereas family-based studies additionally acquired sex-combined residuals with adjustment for sex as a covariate. Case-control studies acquired the residuals in case and control subjects separately. Finally, we rank transformed the residuals using inverse normal transformation to follow a distribution with a mean of 0 and an SD of 1.
Genotyping
All participating studies performed genotyping using the Illumina HumanExome BeadChip. The genotype calling was performed using the designated manufacturer’s software, followed by zCall. Study-specific quality control measures were implemented before the association analyses to remove poorly genotyped variants (Supplementary Table 3).
Study-Level Association Analyses
Associations of the exome-wide variants with the residuals of leptin concentrations were examined using linear mixed models implemented in either RAREMETALWORKER (12) or RVTESTS (13) (Supplementary Table 3). The model accounted for potential cryptic relatedness by incorporating a kinship matrix. We performed the single variant association analyses using both additive and recessive genotypic models. We also calculated covariance matrices capturing linkage disequilibrium (LD) relationships between markers within 1 Mb for use in gene-level meta-analyses.
Quality Control of Study-Level Association Results
We applied the EasyQC package in R to association summary statistics from each participating study to identify cohort-specific quality control issues. This included 1) identifying issues with calculation of leptin residuals and transformation of the residuals, 2) identifying strand issues by comparing allele frequencies against reference alleles from the 1000 Genomes Project phase 1, and 3) identifying issues arising from population stratification.
Single-Variant Meta-analyses
The meta-analyses of summary statistics from the participating studies were carried out with use of RAREMETAL (14) by two different analysts in parallel. We excluded all variants with a call rate <98%, Hardy-Weinberg equilibrium P value <1 × 10−6, or an allele frequency that strongly deviated from the 1000 Genomes reference frequency (>0.60 for all-ancestry analyses and >0.30 for ancestry-specific analyses). To identify the leptin-associated variants, we used the array-wide Bonferroni-corrected threshold of P < 2 × 10−7 for ∼250,000 variants in the single variant analyses.
Gene-Based Meta-analyses
We performed gene-based analyses using the sequence kernel association test (15) (SKAT) and variable threshold (16) (VT) methods in RAREMETAL. The analyses were performed with two different sets of criteria (broad and strict) to select predicted damaging rare and low-frequency variants with minor allele frequency (MAF) <5% annotated using five prediction algorithms: PolyPhen-2, HumDiv, HumVar, LRT, MutationTaster, and SIFT. The broad gene-based tests included nonsense, stop-loss, splice-site, and missense variants that were annotated as damaging by at least one of the five algorithms, whereas the strict tests only included variants predicted as damaging by all of the five algorithms. The statistical significance for the gene-based tests was set at a Bonferroni-corrected threshold of P < 2.5 × 10−6 for 20,000 genes.
Age-Stratified BMI Analyses of Variants in and Near LEP
To study the influence of age on the association of the Val94Met variant in LEP and the rs10487505 variant near LEP with childhood BMI, we performed age-stratified analyses in children with African and European ancestry from the Center for Applied Genomics at The Children’s Hospital of Philadelphia (CHOP) cohort recruited from 2006 to present (17). The participants had multiple BMI measurements at different ages, and analyses were performed with measurements in 1-year age bins. The number of BMI measurements in each age bin is shown in Supplementary Tables 9 and 10. Statistical significance was defined as P < 0.05. The Val94Met and rs10487505 variants were genotyped with the Illumina Infinium II HumanHap550 and Human610-Quad BeadChip and imputed to the HRC r1.1 reference panel with use of the Sanger Imputation Server. All participants were biologically unrelated, aged between 2 and 18 years, and between −3 and +3 SDs of corrected BMI using Centers for Disease Control and Prevention z scores. The study was approved by the institutional review board of CHOP. Parental informed consent was given for each study participant.
Additionally, we used information on comparative body size at age 10 years (data field 1,687) for 452,264 individuals of European ancestry and 8,154 individuals of African ancestry from the UK Biobank. The participants were asked to choose one of the three categories of “about average,” “thinner,” or “plumper” to describe their body size compared with average when they were 10 years old.
Pathway Enrichment Analyses
We used the ExomeChip-Data-driven Expression-Prioritized Integration for Complex Traits (EC-DEPICT) (18,19) gene set enrichment analysis method to evaluate nonsynonymous index variants (strongest nonsynonymous variant within ±1 Mb boundary) with P < 5 × 10−4 for association with either 1) leptin unadjusted for BMI or 2) leptin adjusted for BMI. EC-DEPICT’s primary innovation is the use of “reconstituted” gene sets, which consist of gene sets downloaded from several databases that have been extended based on publicly available large-scale coexpression data (18). Two analyses were performed: 1) all coding variants (N = 93 loci for leptin unadjusted for BMI and N = 91 loci for leptin adjusted for BMI) and 2) coding variants with MAF <5% only (N = 77 loci for leptin unadjusted for BMI and N = 65 loci for leptin adjusted for BMI).
We also used PASCAL (pathway scoring algorithm) (20) to study the enrichment of exome-wide association results in gene sets and pathways with two estimation approaches: MAX and SUM. The MAX estimation is more powerful for single variant–driven associations, whereas the SUM estimation is more powerful when multiple variants are driving the signal (20). We used reconstituted gene sets from DEPICT and the reference data from UK10K (TwinsUK [21] and Avon Longitudinal Study of Parents and Children [ALSPAC] [22]) to estimate LD. The PASCAL analyses were performed for all exome-chip variants (Nall = 265,780 for leptin adjusted for BMI, Nall = 265,780 for leptin unadjusted) and for coding variants only (Ncoding = 176,035 for leptin adjusted for BMI, Ncoding = 180,864 for leptin unadjusted). No allele frequency or P value thresholds were used to select variants for the PASCAL analyses. The pathway scoring method used by PASCAL combines individual gene scores without the need for a tuneable threshold parameter to determine inclusion of genes in the enrichment analysis (20).
Leptin adjusted for BMI is correlated with body fat–free mass (correlation with fat-free mass index [FFMI] in the Fenland cohort = −0.39). The initial pathway analyses for leptin adjusted for BMI using EC-DEPICT and PASCAL suggested enrichment of skeletal muscle–related pathways. To make sure that the gene set enrichment results were not due to correlation between leptin adjusted for BMI and FFMI, we corrected the effect sizes using the following equation (23): βcorrected = βleptin − (βFFMI × rFFMIvs.LEPTIN), where rFFMIvs.LEPTIN = −0.39 (Pearson correlation coefficient in the Fenland Study). The βFFMI coefficients were extracted from an ongoing exome-wide association study of FFMI in ∼500,000 individuals.
Collider Bias
Given that we adjusted leptin concentrations for BMI in our exome-based analyses and leptin and BMI are strongly correlated (r = 0.5–0.8) (9), we tested all exome-based significant loci for evidence of collider bias (23–25). For each index, we extracted the association results from our BMI-unadjusted leptin analyses and from the largest published exome-wide analysis for BMI (19). We corrected BMI-adjusted associations for potential bias due to phenotypic correlation between leptin concentrations and BMI and compared the strength and significance of association with leptin concentrations unadjusted for BMI, leptin adjusted for BMI, and association with BMI (Supplementary Table 6).
Expression Quantitative Trait Locus Colocalization Analyses
The cis expression quantitative trait locus (eQTL) analyses were carried out by use of abdominal subcutaneous adipose tissue from 770 participants of the METSIM (METabolic Syndrome In Men) study, all of whom were Finnish men from Kuopio, Finland (26). The eQTL mapping in 770 METSIM individuals was performed by EPACTS implementing a linear mixed model to account for the population structure among the samples. The eQTLs were defined as cis (local) if the peak association was within 1 Mb on either side of the exon boundaries of the gene. We also identified variants most strongly associated with genes/transcripts from the index variant (“eSNP”). We used METSIM LD (based on n = 770, HRC imputation) to assess LD r2 between the index variant and the lead eSNP. If the pairwise LD was r2 > 0.80, we performed a reciprocal conditional analysis. We tested association between the lead SNP and transcript level when the lead eSNP was included in the model and vice versa.
Expression of the Potential Causal Genes in Preadipocytes and Mature Adipocytes
We compared the expression of the candidate causal genes in the novel leptin-associated loci, including ZNF800, KLF14, KLHL31, ACTL9, CNTD1, and DNAJC18 in preadipocytes and mature adipocytes, two major constituent cell types of adipose tissue. Human preadipocytes isolated from adipose tissue were induced to undergo adipocyte differentiation in vitro (27). RNA samples were obtained from preadipocytes and lipid-laden mature adipocytes at postdifferentiation day 12.
Impact of Val94Met Variant in LEP on Leptin Protein Stability
We used UCSF Chimera 1.13.1 to model the three-dimensional (3D) protein structure and valine-to-methionine substitution in the leptin protein (28). The Rotamers tool and the Dunbrack Rotamer Library were used to view and evaluate amino acid sidechain rotamers. The displayed orientation of methionine was chosen based on the clashes and contacts observed in the protein and hydrogen bonds (29). To predict protein stability, we used SDM (30), iStable (31), CUPSAT (32), and I-Mutant 2.0 (33). All analyses applied the 3D structure for leptin (identifier 1AX8) from the RSCB Protein Data Bank as the reference data set.
Effect of the Val94Met Variant in LEP on Leptin Protein Stability and Secretion Rate
We tested the effect of the Val94Met variant on leptin protein stability and secretion rate in human embryonic kidney (HEK)293 cells in vitro. Human leptin cDNA clone was obtained from Open Biosystems, Inc. (Huntsville, AL) and subcloned into pcDNA3.1 vector. The original cDNA clone encodes the Val94 variant. The 94Met variant was created using Quikchange II Site-Directed Mutagenesis Kit (Agilent Technologies, Santa Clara, CA), with the Val94 plasmid as template and the following mutagenesis primers: forward, 5′-atgccttccagaaacatgatccaaa tatccaac-3′; reverse, 5′-gttggatatttggatcatgtttctggaaggcat-3′). Plasmids carrying Val94 or 94Met cDNA (0.05 μg) were introduced into HEK293 cells (0.65 million cells/well in 12-well plate) using Lipofectamine 2000 as previously described (34). For measurement of intracellular leptin protein turnover and secretion rates, cells were treated with protein synthesis inhibitor cycloheximide (CHX) (20 μg/mL) in fresh media for 0.5 and 1 h at 72 h posttransfection. Cells incubated with fresh media for 1 h without CHX were used as untreated controls. Conditioned media were saved for leptin assay, and cell lysate were prepared using NP-40 lysis buffer (50 mmol/L Tris-HCl, pH 8.0; 150 mmol/L NaCl; 1 mmol/L DTT; 1 mmol/L EDTA; 0.5% NP-40; 10% glycerol; and 1× Roche protein inhibitor mixture). Leptin concentrations in cell lysates and the amount of leptin in conditioned media were determined with a human leptin ELISA kit (R&D Systems, Minneapolis, MN). Little or no cell debris was observed in the conditioned media after centrifugation, suggesting little or no cell breakage during the incubation. The experiments were carried out in duplicates or triplicates and repeated four times.
Data and Resource Availability
The meta-analysis summary results are available in the NHGRI-EBI GWAS Catalog.
Results
Five Novel Genetic Variants Show Association With Leptin Concentrations Independent of Adiposity
To identify genetic variants associated with leptin concentrations, we tested the associations of 246,328 single nucleotide variants, genotyped on an exome-targeted genotyping array, with leptin concentrations in up to 57,232 individuals of European (n = 50,321), African (n = 4,387), East Asian (n = 2,036), or Hispanic (n = 488) ancestry from 35 studies (Supplementary Tables 1 and 3). The exome array provides a detailed coverage of gene-coding regions and includes tags for variants identified in previously published GWAS for human complex traits. Given the strong correlation between leptin and BMI (r = 0.5–0.8) (9), we examined associations with leptin concentrations with and without adjustment for BMI. Additional analyses were performed in men (n = 23,862) and women (n = 32,940) separately. All the analyses were performed in all ancestries combined and in European-ancestry individuals only.
We confirmed five previously established (9) and identified five novel variants associated with leptin concentrations. The novel associations include four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14 (Table 1 and Supplementary Table 4). The associations at already established loci include intergenic variants near LEP and CCNL1, a missense variant in GCKR, and intronic variants in COBLL1 and FTO (Table 1 and Supplementary Table 4). To detect additional independent signals at the 10 leptin-associated loci, we performed conditional analyses, but no further signals were identified.
Leptin-associated loci identified in exome-based association analyses
SNP . | Chr . | Position . | Nearest gene . | Trait . | Most significant model . | Annotation . | EA . | OA . | EAF . | β . | SE . | P . | N . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Novel variants | |||||||||||||
rs3799260 | 6 | 53519605 | KLHL31 | LeptinAdjBMI | Additive/all ancestries/women | Missense | C | T | 0.175 | 0.055 | 0.010 | 1.0E−07 | 32,886 |
rs62621812 | 7 | 127015083 | ZNF800 | LeptinAdjBMI | Additive/all ancestries | Missense | A | G | 0.028 | −0.127 | 0.018 | 2.0E−12 | 56,708 |
rs17151919 | 7 | 127894592 | LEP | LeptinAdjBMI | Additive/all ancestries | Missense | A | G | 0.007 | −0.333 | 0.040 | 1.5E−16 | 49,034 |
rs972283 | 7 | 130466854 | KLF14 | LeptinAdjBMI | Additive/European | Intergenic | A | G | 0.479 | 0.056 | 0.006 | 3.8E−18 | 49,830 |
rs2340550 | 19 | 8808942 | ACTL9 | LeptinAdjBMI | Recessive/European/men | Missense | A | G | 0.316 | 0.071 | 0.014 | 2.0E−07 | 21,883 |
Previously identified variants | |||||||||||||
rs1260326 | 2 | 27730940 | GCKR | LeptinAdjBMI | Additive/all ancestries | Missense | T | C | 0.375 | −0.050 | 0.006 | 2.7E−15 | 56,708 |
rs13389219 | 2 | 165528876 | COBLL1 | LeptinAdjBMI | Additive/all ancestries | Intronic | T | C | 0.410 | 0.053 | 0.007 | 3.0E−15 | 50,297 |
rs900399 | 3 | 156798732 | CCNL1 | LeptinAdjBMI | Additive/all ancestries/women | Intergenic | G | A | 0.391 | −0.054 | 0.008 | 1.2E−10 | 29,510 |
rs791600 | 7 | 127865816 | LEP | LeptinAdjBMI | Additive/all ancestries | Intergenic | A | G | 0.422 | −0.066 | 0.007 | 1.1E−23 | 49,282 |
rs1121980 | 16 | 53809247 | FTO | Leptin | Additive/European | Intronic | A | G | 0.432 | 0.055 | 0.007 | 7.7E−17 | 49,909 |
SNP . | Chr . | Position . | Nearest gene . | Trait . | Most significant model . | Annotation . | EA . | OA . | EAF . | β . | SE . | P . | N . |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Novel variants | |||||||||||||
rs3799260 | 6 | 53519605 | KLHL31 | LeptinAdjBMI | Additive/all ancestries/women | Missense | C | T | 0.175 | 0.055 | 0.010 | 1.0E−07 | 32,886 |
rs62621812 | 7 | 127015083 | ZNF800 | LeptinAdjBMI | Additive/all ancestries | Missense | A | G | 0.028 | −0.127 | 0.018 | 2.0E−12 | 56,708 |
rs17151919 | 7 | 127894592 | LEP | LeptinAdjBMI | Additive/all ancestries | Missense | A | G | 0.007 | −0.333 | 0.040 | 1.5E−16 | 49,034 |
rs972283 | 7 | 130466854 | KLF14 | LeptinAdjBMI | Additive/European | Intergenic | A | G | 0.479 | 0.056 | 0.006 | 3.8E−18 | 49,830 |
rs2340550 | 19 | 8808942 | ACTL9 | LeptinAdjBMI | Recessive/European/men | Missense | A | G | 0.316 | 0.071 | 0.014 | 2.0E−07 | 21,883 |
Previously identified variants | |||||||||||||
rs1260326 | 2 | 27730940 | GCKR | LeptinAdjBMI | Additive/all ancestries | Missense | T | C | 0.375 | −0.050 | 0.006 | 2.7E−15 | 56,708 |
rs13389219 | 2 | 165528876 | COBLL1 | LeptinAdjBMI | Additive/all ancestries | Intronic | T | C | 0.410 | 0.053 | 0.007 | 3.0E−15 | 50,297 |
rs900399 | 3 | 156798732 | CCNL1 | LeptinAdjBMI | Additive/all ancestries/women | Intergenic | G | A | 0.391 | −0.054 | 0.008 | 1.2E−10 | 29,510 |
rs791600 | 7 | 127865816 | LEP | LeptinAdjBMI | Additive/all ancestries | Intergenic | A | G | 0.422 | −0.066 | 0.007 | 1.1E−23 | 49,282 |
rs1121980 | 16 | 53809247 | FTO | Leptin | Additive/European | Intronic | A | G | 0.432 | 0.055 | 0.007 | 7.7E−17 | 49,909 |
The chromosomal positions are based on hg19. Chr, chromosome; EA, effect allele; EAF, effect allele frequency; LeptinAdjBMI, leptin adjusted for BMI; OA, other allele.
The association between rs1121980 near FTO and leptin concentrations became nonsignificant after adjustment for BMI (Punadjusted = 8 × 10−17; PBMI adjusted = 0.45). The effects of all other known and novel loci were independent of BMI (Supplementary Table 5). We tested whether the adjustment for BMI, a strongly correlated covariate (23), may have introduced collider bias but found no evidence of such bias (Supplementary Table 6).
The strongest variant associated with leptin concentrations was rs791600, an intergenic variant near the LEP gene. The rs791600 variant is in LD (European [EUR] r2 = 0.70) with the rs10487505 variant identified in our previously published GWAS (9), which is not included in the exome array and was therefore not available for analyses in the current study. In the prior GWAS, the rs10487505 variant showed a more significant association with BMI-adjusted leptin concentrations (β = 0.034 per allele, P = 2.7 × 10−11, n = 29,252) than rs791600 (β = 0.029 per allele, P = 3.0 × 10−9, n = 31,800) and thus is still considered the lead variant at this locus (Supplementary Fig. 1).
Nine of the 10 identified loci showed an association with leptin concentrations in all ancestries combined and in European-ancestry-only analyses. However, the novel LEP variant Val94Met (rs17151919) only showed a significant association for all ancestries combined (P = 2 × 10−16) and not in European-ancestry individuals alone (P = 0.47) (Supplementary Table 7). In further ancestry-stratified analyses, we observed that the Met94-coding allele is common in populations of African ancestry (MAF = 8%), less common in those with Hispanic ancestry (MAF = 2%), very rare in those with European ancestry (MAF = 0.02%), and monomorphic in people with East Asian ancestry (35). In individuals of African descent, each Met94-coding allele was associated with 0.34 SDs lower leptin concentrations (P = 2 × 10−16, n = 3,901) (Supplementary Fig. 2). The direction of effect was consistent in individuals with Hispanic (−0.21 SD effect per allele, P = 0.29, n = 488) and individuals with European (−0.19 SD effect per allele, P = 0.47, n = 44,401) ancestry but did not reach statistical significance, most likely because very few carriers were available (NHispanic = 24, NEUR = 15) (Supplementary Table 7).
Gene-Based Analysis Identifies Two Novel Genes With Sex-Specific Effect on Leptin
In addition to single variant–based association tests, we performed gene-based tests using rare and low-frequency coding variants in aggregate (15,16) (research design and methods). We identified two genes associated with leptin concentrations. CNTD1 showed association with leptin concentrations unadjusted for BMI in men (P = 1 × 10−7) but not in women (P = 0.27) (Table 2 and Supplementary Table 8). The association in men was driven by five coding variants and was strongly attenuated by adjustment for BMI (P = 0.007), suggesting that the association of CNTD1 with leptin concentrations may be due to a link between CNTD1 and adiposity, although no such connection has previously been reported. The CNTD1 gene encodes cyclin N-terminal domain-containing 1, which is critical for meiotic crossover maturation and deselection of excess precrossover sites.
Leptin-associated genes identified by gene-based exome-wide association analyses
Gene . | Chr . | Position . | Trait . | Most significant model . | Method . | N . | P . | β . | SE . | N variants . |
---|---|---|---|---|---|---|---|---|---|---|
CNTD1 | 17 | 40950810–40963605 | Leptin | Additive/European/men | SKAT broad | 18,882 | 1.3E−07 | 0.898 | 0.165 | 5 |
DNAJC18 | 5 | 138743559–138780898 | LeptinAdjBMI | Additive/all ancestries/women | SKAT strict | 29,510 | 5.5E−08 | 0.757 | 0.169 | 2 |
Gene . | Chr . | Position . | Trait . | Most significant model . | Method . | N . | P . | β . | SE . | N variants . |
---|---|---|---|---|---|---|---|---|---|---|
CNTD1 | 17 | 40950810–40963605 | Leptin | Additive/European/men | SKAT broad | 18,882 | 1.3E−07 | 0.898 | 0.165 | 5 |
DNAJC18 | 5 | 138743559–138780898 | LeptinAdjBMI | Additive/all ancestries/women | SKAT strict | 29,510 | 5.5E−08 | 0.757 | 0.169 | 2 |
The chromosomal positions are based on hg19. Chr, chromosome; LeptinAdjBMI, leptin adjusted for BMI.
Another gene, DNAJC18, showed association with BMI-adjusted leptin concentrations in women (P = 6 × 10−8) but not men (P = 0.02). The association in women was driven by two coding variants (Table 2 and Supplementary Table 8). DNAJC18 is part of the DnaJ Hsp family. However, no function has yet been described to C18 subfamily.
LEP Val94Met Regulates Leptin Secretion and Early Adiposity
The Val94Met (rs17151919) LEP variant was associated with BMI-adjusted leptin concentrations in individuals of African ancestry. A previous study in 2,129 African Americans in the Coronary Artery Risk Development Study in Young Adults (CARDIA) (not included in the present meta-analyses) reported a significant association between the leptin-decreasing Met94 allele of the Val94Met (rs17151919) variant in LEP and up to 1.12 kg/m2 higher BMI in adulthood (P = 0.018) (36). However, results from two larger studies of adult BMI by the African Ancestry Anthropometry Genetics Consortium (n = 42,752, P = 0.88) (37) and the African ancestry population of the UK Biobank study (N = 7,820, P = 0.17) did not replicate the association. Nevertheless, among the African ancestry population in the UK Biobank, carriers of the leptin-decreasing Met94 allele reported more often that at age 10 years they were “plumper” compared with peers (odds ratio 1.11, P = 0.04), suggesting that the effect of this variant may be age dependent. To study the influence of age, we performed age-stratified analyses in up to 2,726 children with African ancestry from the Center for Applied Genomics at CHOP cohort (17). Comparison of the effect sizes across different age points revealed that each leptin-decreasing Met94 allele was associated with 0.12–0.20 units higher BMI z score between the ages 3 and 7 years (P < 0.05). The most pronounced effect was reached at age 6 years (Fig. 1 and Supplementary Table 9), and no association with BMI was observed after age 8 years (β = −0.04 to 0.05) (Fig. 1 and Supplementary Table 9), suggesting that the BMI-increasing effect of the Met94 allele wanes shortly before puberty. The rs10487575 variant near LEP showed a trajectory of association with childhood BMI similar to that of the Val94Met variant, but the effect sizes were much more modest (Fig. 1 and Supplementary Table 10), consistent with the fivefold smaller effect of rs10487505 on leptin concentrations compared with Val94Met in adults (Table 1).
Association of the leptin-decreasing alleles of the LEP Val94Met (rs17151919) variant (on the left) and the rs10487505 variant near LEP (on the right) with BMI SD score (BMI-SDS) in the CHOP cohort. The analyses for the Val94Met variant were performed in up to 2,726 African-ancestry participants and the analyses for the rs10487505 variant in up to 3,681 African- and European-ancestry participants of the CHOP cohort. The y-axis shows the effect of each leptin-decreasing allele on BMI at each age year. The error bars indicate 1 SEM.
Association of the leptin-decreasing alleles of the LEP Val94Met (rs17151919) variant (on the left) and the rs10487505 variant near LEP (on the right) with BMI SD score (BMI-SDS) in the CHOP cohort. The analyses for the Val94Met variant were performed in up to 2,726 African-ancestry participants and the analyses for the rs10487505 variant in up to 3,681 African- and European-ancestry participants of the CHOP cohort. The y-axis shows the effect of each leptin-decreasing allele on BMI at each age year. The error bars indicate 1 SEM.
The Val94Met variant is located at position 94 in the 167 amino acid leptin precursor protein and results in a valine-to-methionine change at position 73 of the mature protein (Fig. 2A). Position 73 is situated at the leptin protein surface and is not believed to be involved in binding of leptin to its receptor. Nevertheless, structural prediction tools (30–33) suggested that the substitution of valine with methionine at this position is likely to lead to reduced stability of the mature leptin protein (Fig. 2A, Supplementary Figs. 3 and 4, and Supplementary Table 11). This is consistent with our observation that the methionine-coding allele is associated with lower leptin concentrations.
Impact of Val94Met transversion at LEP rs17151919 on leptin secretion rate in HEK293 cells. The rs17151919 variant changes valine to methionine in position 73 of the mature leptin protein. A: The 3D illustration of leptin structure derived from RSCB Protein Data Bank (PDB) and modified with UCSF Chimera1.13.1. The prediction of protein stability is derived from the SDM server (30). B: Leptin secretion rates for Val94 and Met94 expressed as the amount of leptin secreted (ng) during a 1-h incubation (72–73 h posttransfection) (LEPs/h) normalized by the respective cellular leptin content (LEPc) in untreated control cells at the end of incubation. Individual data points from four separate experiments (each with two to three technical replicates) are plotted. The normality of data distribution was examined using D’Agostino and Pearson normality test (P = 0.65 and 0.54 for LEPV94 (Val94) and LEPM94 (Met94), respectively), and repeated-measures one-way ANOVA was performed to assess the difference in secretion rate between the genotypes. Mean ± SD and AVOVA results (F and P values) are reported in the table below the graph. C: Intracellular leptin turnover rates for Val94 and Met94 alleles, obtained by measuring the relative cellular leptin contents in the untreated control cells (defined as 1 for the respective LEP variant) and in samples treated with the protein synthesis inhibitor CHX (20 μg/mL) for 0.5 and 1.0 h. Mean ± SD at each time point from four separate experiments (each with two to three technical replicates) is plotted. Paired t test was used to assess the genotype effect on the fractions of cellular LEP that remained after 0.5 h and 1.0 h of CHX treatment (P values are reported in the table below the graph). The average hourly turnover rates for Val94 and Met94 were 61 ± 2% and 60 ± 3%, respectively, calculated by subtraction of the percent cellular LEP that remained after 1 h of CHX treatment from those of the respective untreated controls (defined as 100%). hr, hour(s); ID, identifier; M94, Met94; pseudo ΔΔG, predicted stability difference score; SNP, single nucleotide polymorphism; V94, Val94; WT, wild type.
Impact of Val94Met transversion at LEP rs17151919 on leptin secretion rate in HEK293 cells. The rs17151919 variant changes valine to methionine in position 73 of the mature leptin protein. A: The 3D illustration of leptin structure derived from RSCB Protein Data Bank (PDB) and modified with UCSF Chimera1.13.1. The prediction of protein stability is derived from the SDM server (30). B: Leptin secretion rates for Val94 and Met94 expressed as the amount of leptin secreted (ng) during a 1-h incubation (72–73 h posttransfection) (LEPs/h) normalized by the respective cellular leptin content (LEPc) in untreated control cells at the end of incubation. Individual data points from four separate experiments (each with two to three technical replicates) are plotted. The normality of data distribution was examined using D’Agostino and Pearson normality test (P = 0.65 and 0.54 for LEPV94 (Val94) and LEPM94 (Met94), respectively), and repeated-measures one-way ANOVA was performed to assess the difference in secretion rate between the genotypes. Mean ± SD and AVOVA results (F and P values) are reported in the table below the graph. C: Intracellular leptin turnover rates for Val94 and Met94 alleles, obtained by measuring the relative cellular leptin contents in the untreated control cells (defined as 1 for the respective LEP variant) and in samples treated with the protein synthesis inhibitor CHX (20 μg/mL) for 0.5 and 1.0 h. Mean ± SD at each time point from four separate experiments (each with two to three technical replicates) is plotted. Paired t test was used to assess the genotype effect on the fractions of cellular LEP that remained after 0.5 h and 1.0 h of CHX treatment (P values are reported in the table below the graph). The average hourly turnover rates for Val94 and Met94 were 61 ± 2% and 60 ± 3%, respectively, calculated by subtraction of the percent cellular LEP that remained after 1 h of CHX treatment from those of the respective untreated controls (defined as 100%). hr, hour(s); ID, identifier; M94, Met94; pseudo ΔΔG, predicted stability difference score; SNP, single nucleotide polymorphism; V94, Val94; WT, wild type.
To study the impact of the Val94Met variant on the intracellular turnover of the leptin protein and its secretion rate, we performed in vitro experiments in HEK293 cells. Leptin secretion rate—calculated as the amount of leptin secreted in 1 h normalized to the respective cellular leptin content—was 20.4% lower in Met94 than in Val94 cells (P = 0.0007 by repeated-measures one-way ANOVA) 72 h posttransfection (Fig. 2B). Leptin secretion rates between 48 and 72 h posttransfection and during a 1-h treatment with CHX were 11.8% (P = 0.0005) and 17.9% (P = 0.0002) lower, respectively, in Met94 compared with Val94 (Supplementary Fig. 5). Notably, no difference was found in the intracellular turnover rate of leptin between Val94 and Met94 cells during a 0.5- or 1.0-h incubation with CHX to impair protein synthesis (Fig. 2C). The unchanged turnover rate incorporates protein secretion and degradation, suggesting that decreased leptin secretion rate was likely associated with increased intracellular leptin degradation in Met94 cells. Overall, these in vitro experiments suggest that methionine substitution in position 73 of the mature leptin protein decreases the rate of leptin secretion from the cells, which may contribute to the association of the Met94 allele with lower leptin concentrations.
ZNF800 Locus Regulates Adipose Gene Expression and Body Composition
The Pro103Ser (rs62621812) variant in ZNF800 changes the amino acid sequence of CH2 zinc finger protein, a putative transcription factor (38). We found that the Ser103 allele (frequency 2.8%) is associated with lower BMI-adjusted leptin concentrations (P = 2.0 × 10−12). As shown before (26), Pro103Ser is the lead variant associated with expression of ZNF800 in subcutaneous adipose tissue in the Finnish METSIM study (P = 2.4 × 10−16); the Ser103 allele is associated with higher ZNF800 expression levels (Supplementary Table 12). ZNF800 is a master regulator in subcutaneous adipose tissue, as the Pro103Ser variant has also been associated with adipose tissue expression of nine other genes (26). In the eQTL data, the leptin-decreasing Ser103 allele was not significantly associated with the expression of LEP (P = 0.20), located 866 kb downstream, and the observed direction of the effect on LEP expression was opposite that observed for leptin concentrations (β = 0.14 SD/allele vs. β = −0.13 SD/allele, respectively), suggesting that the leptin-lowering effect of the Ser103 allele on leptin concentrations is unlikely to be mediated by direct transcriptional regulation of LEP.
In the UK Biobank study, we found that each leptin-decreasing Ser103 allele is associated with 0.14 kg/m2 higher BMI (P = 8.1 × 10−6). However, there was no association between Ser103 allele and body fat percentage (−0.045% per allele, P = 0.25), indicating that the variant impacts BMI primarily by increasing fat-free body mass. Indeed, the leptin-decreasing Ser103 allele was associated with a 0.33 kg higher fat-free mass (P = 4.6 × 10−20) and only 0.13 kg higher fat mass (P = 0.023). The Ser103 allele is associated with higher expression of the ZNF800 gene in the tibial nerve (GTEx V8) (P = 1.4 × 10−6, n = 532) that innervates the muscles of the leg and was previously identified for association with increased appendicular lean mass (39). There was no association between the Ser103 allele and self-reported body size at age 10 years (P = 0.75).
The KLF14 Locus Regulates Adipogenesis and Fat Distribution
The rs972283 variant (MAFEUR = 48%), associated with leptin concentrations, is located 51 kb upstream from KLF14 and 2.5 Mb downstream from LEP and is in near-perfect LD with previously reported GWAS variants for type 2 diabetes (40), insulin resistance (41), HDL cholesterol (42), and body fat distribution (43). As previously reported (26), rs972283 is associated with KLF14 expression in subcutaneous adipose tissue (Supplementary Table 12). As KLF14 is a master regulator in adipose tissue, rs972283 is also associated with the expression of multiple other genes in trans (44). No significant association was observed between rs972283 and LEP expression in the METSIM eQTL study (44), suggesting that KLF14 may not regulate leptin production at the transcriptional level, at least not in men. Lower expression of KLF14 has been implicated in impaired adipogenesis due to defective adipocyte glucose uptake in women, characterized by the presence of fewer but larger adipocytes and a shift in fat distribution from gynoid stores to abdominal tissues (44). However, while the effects of KLF14 on adipogenesis and adipose redistribution have been found to be specific to women, there was no difference in the association of rs972283 with leptin levels between men and women (Supplementary Table 4).
Interestingly, the carriers of the rs972283-G allele more frequently reported being plumper (P = 2.8 × 10−5) and shorter (P = 0.014) than average at age 10 years in the UK Biobank than noncarriers, whereas the same allele was associated with a lower BMI (P = 6.8 × 10−9) and increased height (P = 0.010) in adults, suggesting that the effect of the rs972283 variant on body size may change during life course. In previous GWAS of adults, the rs972283-G allele was identified as associated with higher risk of type 2 diabetes (40) and insulin resistance (41) and lower hip circumference (adjusted for BMI) (43) and HDL cholesterol (42). In the UK Biobank study, the rs972283-G allele was associated with lower body fat percentage in adults (P = 5.9 × 10−22).
The KLHL31 Locus Is Implicated in Adipogenesis in Adult Females
The Val156Ile (rs3799260) variant (MAFEUR = 18%) in KLHL31, associated with leptin concentrations in female-only analyses, changes the amino acid sequence of the Kelch-like family member 31 protein. KLHL31 suppresses Wnt–β-catenin signaling that is involved in promoting adipocyte differentiation and suppressing oxidative metabolism in adipocytes. The Val156Ile variant is predicted to be benign/tolerated by SIFT/PolyPhen (45,46). Previous genetic associations have identified a variant in low LD (rs7739232) (EUR r2 = 0.27) to be associated with BMI-adjusted hip circumference, also specific to women (43). The rs7739232 variant was not included in the exome array and was thus not analyzed in the current study. Our in vitro experiments showed that KLHL31 is only expressed in mature adipocytes, but not in preadipocytes (Supplementary Fig. 6), suggesting that the gene is developmentally regulated.
In the UK Biobank, similar to the variants in and near LEP, the leptin-decreasing Ile156 allele carriers more often reported being plumper than average at age 10 years (P = 5.6 × 10−6), but there was a weaker association with higher BMI (P = 0.045) in adulthood.
In Men, the ACTL9 Locus May Regulate Leptin Concentrations in a Cell-Nonautonomous Fashion
Homozygosity for the minor allele of the Ser37Phe (rs2340550) variant in ACTL9 was associated with leptin concentrations in men only in a recessive genetic model. While the Ser37Phe variant is predicted to be benign/tolerated by SIFT/PolyPhen, another missense variant, Ala51Val (rs10410943), in high LD (EUR r2 = 0.99) is predicted to be deleterious/probably damaging and could be the causal variant at the locus. The Ser37Phe variant is also in high LD (r2 > 0.8) with several nearby noncoding variants (Supplementary Fig. 7). However, none of these overlap with regulatory elements in adipocytes. The expression of ACTL9 is restricted to the testis, and it is therefore likely to act in an adipocyte-nonautonomous fashion to influence leptin concentrations. Actin proteins have cytoskeletal functions and have also been implicated in signaling and nuclear activities.
Gene Set Analyses Implicate Adipocyte-Related Pathways
We performed gene set enrichment analyses using EC-DEPICT (18,19,47) and PASCAL (20) to identify biological processes and candidate pathways enriched for loci associated with leptin unadjusted or adjusted for BMI. Among coding variants associated with BMI-unadjusted leptin concentrations, PASCAL identified significant enrichment of the gene set for “positive regulation of reproductive success” (Pempirical = 1.6 × 10−5) (Supplementary Table 13), consistent with the crucial permissive role of leptin in the integrity of the gonadal axis (48). Among coding variants associated with leptin adjusted for BMI, we found enrichment of the immune-related TRIM39 protein-protein interaction subnetwork (49,50) (Pempirical = 8.4 × 10−6) (Supplementary Table 14). No gene sets were found to be significantly enriched in PASCAL analyses where all exome-wide variants (coding and noncoding) for leptin adjusted for BMI were included or in the EC-DEPICT analyses.
Discussion
We identified 10 genetic variants associated with leptin concentrations and two gene-based associations using an exome-based genotyping array in up to 57,232 individuals with varying ancestries. The two independent variants most strongly associated with leptin concentrations were located in and near the LEP gene. The African ancestry–driven variant within LEP, Val94Met (rs17151919), was found to decrease leptin secretion in HEK293 cells, whereas rs10487505, located near LEP, overlaps a long noncoding (lnc)RNA that regulates LEP expression (10). Both variants showed significant association with increased adiposity in children, whereas only a nominal or no association was observed in adults.
Previous analyses have shown that the leptin-lowering allele of rs10487505 is only weakly associated with higher BMI in adulthood but shows a pronounced association with BMI in early childhood (10,11). Similarly, we showed that the LEP Met94 allele, associated with lower leptin concentration, is associated with early childhood BMI. Our results suggest that leptin has an impact specifically on early adiposity, encouraging further studies to uncover the molecular mechanisms that underlie this age-dependent relationship between leptin and BMI.
The Val94Met and rs10487505 variants in and near LEP are likely to influence leptin concentrations by different molecular mechanisms. The novel African ancestry–driven variant Val94Met may affect circulating levels of leptin by reducing leptin secretion. The rs10487505 variant is associated with leptin mRNA levels in adipose tissue. We hypothesize that this variant, located upstream of LEP within an lncRNA (EL947753), interacts with enhancer regions to regulate the expression of LEP. Defects in LEP regulation in mice lead to a relative hypoleptinemic form of obesity that is responsive to leptin administration (10).
We identified four new loci associated with leptin concentrations, located in or near the ZNF800, KLF14, KLHL31, and ACTL9 genes. Two additional genes, CNTD1 and DNAJC18, were identified in gene-based analyses. The ZNF800 and KLF14 genes are master trans regulators of adipose tissue gene expression (26) and located in the proximity of the LEP gene (866 kb and 2.5 Mb away, respectively). The variants in ZNF800 and near KLF14 were not associated with LEP mRNA levels, however, suggesting that they may be involved in translational or posttranslational rather than transcriptional regulation of leptin production. KLHL31 has been shown to promote adipocyte differentiation and suppress oxidative metabolism in adipocytes, whereas ACTL9 is not expressed in adipocytes and could affect circulating levels by a non–cell-autonomous mechanism. The KLHL31 and ACTL9 loci, and the CNTD1 and DNAJC18 genes, were only identified in sex-specific models and narrowly passed the array-wide significance threshold. Further validation of the association of these loci with leptin concentrations is warranted.
In summary, we identified a new genetic association of an African ancestry–specific missense variant rs17151919 in LEP with leptin concentrations and replicated the association of the rs10487505 variant near LEP. The pronounced association of these variants with BMI in early childhood implicates genetic regulation of LEP in early growth and suggests that young children may be particularly sensitive to the metabolic/behavioral effects of leptin. We also identified novel loci at ZNF800, KLF14, KLHL31, ACTL9, CNTD1, and DNAJC18 associated with leptin concentrations, providing additional insights into leptin physiology.
L.O.H. is currently affiliated with Steno Diabetes Center Copenhagen, Gentofte, Denmark. R.S.F. is currently affiliated with Vertex Pharmaceuticals Incorporated, Boston, MA. M.E.G. is currently affiliated with the National Heart, Lung, and Blood Institute (NHLBI), Bethesda, MD.
This work was conducted prior to M.E.G.’s current affiliation with the NHLBI, and, as such, the views expressed in this article do not represent the views of the NHLBI, National Institutes of Health, or other government entity.
This article contains supplementary material online at https://doi.org/10.2337/figshare.12933737.
H.Y., Y.Z., R.J.F.L., and T.O.K. contributed equally to this study.
Article Information
Acknowledgments. The authors thank all investigators, study members, and study participants for their contributions to the participating studies. The acknowledgments for all participating studies are provided in the Supplementary Material. This research has been conducted using the UK Biobank resource and carried out under UK Biobank project 9072. Details on patient and public involvement in the UK Biobank are available online at www.ukbiobank.ac.uk/about-biobank-uk and www.ukbiobank.ac.uk/wp-content/uploads/2011/07/Summary-EGF-consultation.pdf?phpMyAdmin=trmKQlYdjjnQIgJ%2CfAzikMhEnx6.
Funding and Duality of Interest. The funding information for all participating studies are provided in the Supplementary Material. D.M.-K. is a part-time clinical research consultant for Metabolon, Inc. M.A.N.’s participation is supported by a consulting contract between Data Tecnica International and the National Institute on Aging, National Institutes of Health. V.S. has served on advisory boards for Novo Nordisk and Sanofi and received honoraria from these companies. V.S. also has ongoing research collaboration with Bayer Ltd. B.M.P. serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. J.R.K. reports stock ownership in Bristol-Myers Squibb, Johnson & Johnson, Merck, and Pfizer. R.S.F. is currently an employee of Vertex Pharmaceuticals Incorporated, though all work for this article was done prior to Vertex affiliation. This work was supported by NHGRI NIH F31HG009850. E.I. is now an employee of GlaxoSmithKline, but this work was done prior to joining the company. No other potential conflicts of interest relevant to this article were reported.
The ongoing research by V.S. in collaboration with Bayer Ltd. is unrelated to the current study.
Author Contributions. H.Y., Y.Z., R.J.F.L., and T.O.K. comprised the central analysis and writing team. R.L.-G, N.G., X.G., L.-P.L., M.F.F., A.M., T.Kac., D.I.C., N.A., L.B., N.R., K.L.Y., J.B.B., J.B.-J, P.C., A.D., M.G., Y.-D.I.C., P.J., J.K., J.Lu., S.Mä., C.M.-G., D.O.M.-K., B.M.P., R.A.S., K.S., A.T., A.G.U., Y.W., P.K., C.M.v.D., P.M.R., G.H., E.I., M.A.P., T.L., L.S.A., J.I.R., J.G.W., H.G., T.H., N.J.W., C.L., C.M.L., and K.L.M. contributed to genotyping. H.Y., Y.Z., C.N.S., T.Kar., L.O.H., J.B., C.S., R.S.F., M.H.P., Z.K., L.B.L.W., Y.U., S.M., S.M.W., R.L.-G., N.G., S.W., S.Mo., A.A.S., M.A.N., L.A.L., J.Ha., X.G., L.-P.L., M.F.F., C.M.S., C.V., A.M., T.Kac., C.A.W., D.I.C., N.A., L.B., N.R., K.L.Y., P.L.A., J.B.-J., A.D., C.A.D., K.G., H.H., H.J., J.Lu., M.M., C.M.-G., A.P.M., I.P., K.S., A.T., J.Y., H.G., K.L.M., R.L.L., and T.O.K. contributed to statistical analysis within cohorts. M.F.F., C.V., N.A., K.L.Y., M.A., M.B., J.B.B., G.D.C., A.D., M.E.G., M.G., J.Ho., R.J., P.J., A.E.J., M.K., J.R.K., C.A.L., J. Li, L.L., D.A.M., S.Mä., J.F.M.C., C.M.-G., D.O.M.-K., R.d.M., M.N., A.D.P., B.M.P., O.T.R., T.S., K.D.T., Y.W., M.A.I., C.M.v.D., M.A.P., T.L., L.S.A., J.I.R., J.G.W., H.G., J.B.M., N.J.W., T.M.F., and C.L. contributed to sample collection and phenotyping. N.G., D.I.C., R.J., M.K., L.L., M.M., C.E.P., O.T.R., K.S., A.G.U., M.W., K.E.N., M.A.I., P.M.R., S.L., E.I., T.D.S., B.M., T.L., L.S.A., J.I.R., A.P.R., T.B.H., S.R., H.G., J.B.M., V.S., T.H., K.W.v.D., N.J.W., S.F.A.G., C.L., and K.L.M. were the Principal Investigators within cohorts. H.Y. and T.O.K. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.