Estimated prevalence of potentially damaging variants in the leptin gene

Mutations in the leptin gene (LEP) can alter the secretion or interaction of leptin with its receptor, leading to extreme early-onset obesity. The purpose of this work was to estimate the prevalence of heterozygous and homozygous mutations in the leptin gene with the help of the Exome Aggregation Consortium (ExAC) database ( http://exac.broadinstitute.org/about ).

The ExAC database encompasses exome sequencing data from 60,706 individuals. We searched for listed leptin variants and identified 36 missense, 1 in-frame deletion, and 3 loss-of-function variants. The functional relevance of these variants was assessed by the in silico prediction tools PolyPhen-2, Sorting Intolerant from Tolerant (SIFT), and Loss-Of-Function Transcript Effect Estimator (LOFTEE). PolyPhen-2 predicted 7 of the missense variants to be probably damaging and 10 to be possibly damaging. SIFT predicted 7 of the missense variants to be deleterious. Three loss-of-function variants were predicted by LOFTEE. Excluding double counts, we can summarize 21 variants as potentially damaging. Considering the allele count, we identified 31 heterozygous but no homozygous subjects with at least probably damaging variants. In the ExAC population, the estimated prevalence of heterozygous carriers of these potentially damaging variants was 1:2000. The probability of homozygosity was 1:15,000,000. We furthermore tried to assess the functionality of ExAC-listed leptin variants by applying a knowledge-driven approach. By this approach, additional 6 of the ExAC-listed variants were considered potentially damaging, increasing the number of heterozygous subjects to 58, the prevalence of heterozygosity to 1:1050, and the probability of homozygosity to 1:4,400,000.

Using exome sequencing data from ExAC, in silico prediction tools and by applying a knowledge-driven approach, we identified 27 probably damaging variants in the leptin gene of 58 heterozygous subjects. With this information, we estimate the prevalence for heterozygosity at 1:1050 corresponding to an incidence of homozygosity of 1:4,400,000 in this large pluriethnic cohort.