Novel common copy number variation for early onset extreme obesity on chromosome 11q11 identified by a genome-wide analysis.

Heritability of obesity is substantial and recent meta-analyses of genome-wide association studies (GWASs) have been successful in detecting several robustly associated genomic regions for obesity using single-nucleotide polymorphisms (SNPs). However, taken together, the SNPs explain only a small proportion of the overall heritability. Copy number variations (CNVs) might contribute to the 'missing heritability'. We searched genome-wide for association between common CNVs and early-onset extreme obesity. Four hundred and twenty-four case-parents obesity trios and an independent sample of 453 extremely obese children and adolescents and 435 normal-weight and lean adult controls were genotyped by the Affymetrix Genome-Wide Human SNP Array 6.0. We detected 20 common copy number variable regions (CNVRs) which were associated with obesity. The most promising CNVRs were followed-up in an independent sample of 365 obesity trios, confirming the association for two candidate CNVRs. We identified a common CNVR exclusively covering the three olfactory receptor genes OR4P4, OR4S2 and OR4C6 to be associated with obesity (combined P-value = 0.015 in a total of 789 families; odds ratio for the obesity effect allele = 1.19; 95% confidence interval = 1.016-1.394). We also replicated two common deletions (near NEGR1 and at chromosome 10q11.22) that have previously been reported to be associated with body weight. Additionally, we support a rare CNV on chromosome 16 that has recently been reported by two independent groups. However, rare CNVs had not been the focus of our study. We conclude that common CNVs are unlikely to contribute substantially to the genetic basis of early-onset extreme obesity.

Two new Loci for body-weight regulation identified in a joint analysis of genome-wide association studies for early-onset extreme obesity in French and german study groups.

Meta-analyses of population-based genome-wide association studies (GWAS) in adults have recently led to the detection of new genetic loci for obesity. Here we aimed to discover additional obesity loci in extremely obese children and adolescents. We also investigated if these results generalize by estimating the effects of these obesity loci in adults and in population-based samples including both children and adults. We jointly analysed two GWAS of 2,258 individuals and followed-up the best, according to lowest p-values, 44 single nucleotide polymorphisms (SNP) from 21 genomic regions in 3,141 individuals. After this DISCOVERY step, we explored if the findings derived from the extremely obese children and adolescents (10 SNPs from 5 genomic regions) generalized to (i) the population level and (ii) to adults by genotyping another 31,182 individuals (GENERALIZATION step). Apart from previously identified FTO, MC4R, and TMEM18, we detected two new loci for obesity: one in SDCCAG8 (serologically defined colon cancer antigen 8 gene; p = 1.85x10(-8) in the DISCOVERY step) and one between TNKS (tankyrase, TRF1-interacting ankyrin-related ADP-ribose polymerase gene) and MSRA (methionine sulfoxide reductase A gene; p = 4.84x10(-7)), the latter finding being limited to children and adolescents as demonstrated in the GENERALIZATION step. The odds ratios for early-onset obesity were estimated at approximately 1.10 per risk allele for both loci. Interestingly, the TNKS/MSRA locus has recently been found to be associated with adult waist circumference. In summary, we have completed a meta-analysis of two GWAS which both focus on extremely obese children and adolescents and replicated our findings in a large followed-up data set. We observed that genetic variants in or near FTO, MC4R, TMEM18, SDCCAG8, and TNKS/MSRA were robustly associated with early-onset obesity. We conclude that the currently known major common variants related to obesity overlap to a substantial degree between children and adults.

Genome-wide association study in German patients with attention deficit/hyperactivity disorder.

The heritability of attention deficit hyperactivity disorder (ADHD) is approximately 0.8. Despite several larger scale attempts, genome-wide association studies (GWAS) have not led to the identification of significant results. We performed a GWAS based on 495 German young patients with ADHD (according to DSM-IV criteria; Human660W-Quadv1; Illumina, San Diego, CA) and on 1,300 population-based adult controls (HumanHap550v3; Illumina). Some genes neighboring the single nucleotide polymorphisms (SNPs) with the lowest P-values (best P-value: 8.38 × 10(-7)) have potential relevance for ADHD (e.g., glutamate receptor, metabotropic 5 gene, GRM5). After quality control, the 30 independent SNPs with the lowest P-values (P-values ≤ 7.57 × 10(-5) ) were chosen for confirmation. Genotyping of these SNPs in up to 320 independent German families comprising at least one child with ADHD revealed directionally consistent effect-size point estimates for 19 (10 not consistent) of the SNPs. In silico analyses of the 30 SNPs in the largest meta-analysis so far (2,064 trios, 896 cases, and 2,455 controls) revealed directionally consistent effect-size point estimates for 16 SNPs (11 not consistent). None of the combined analyses revealed a genome-wide significant result. SNPs in previously described autosomal candidate genes did not show significantly lower P-values compared to SNPs within random sets of genes of the same size. We did not find genome-wide significant results in a GWAS of German children with ADHD compared to controls. The second best SNP is located in an intron of GRM5, a gene located within a recently described region with an infrequent copy number variation in patients with ADHD.

Analysis of Genes Involved in Body Weight Regulation by Targeted Re-Sequencing.

INTRODUCTION: Genes involved in body weight regulation that were previously investigated in genome-wide association studies (GWAS) and in animal models were target-enriched followed by massive parallel next generation sequencing. METHODS: We enriched and re-sequenced continuous genomic regions comprising FTO, MC4R, TMEM18, SDCCAG8, TKNS, MSRA and TBC1D1 in a screening sample of 196 extremely obese children and adolescents with age and sex specific body mass index (BMI) ≥ 99th percentile and 176 lean adults (BMI ≤ 15th percentile). 22 variants were confirmed by Sanger sequencing. Genotyping was performed in up to 705 independent obesity trios (extremely obese child and both parents), 243 extremely obese cases and 261 lean adults. RESULTS AND CONCLUSION: We detected 20 different non-synonymous variants, one frame shift and one nonsense mutation in the 7 continuous genomic regions in study groups of different weight extremes. For SNP Arg695Cys (rs58983546) in TBC1D1 we detected nominal association with obesity (pTDT = 0.03 in 705 trios). Eleven of the variants were rare, thus were only detected heterozygously in up to ten individual(s) of the complete screening sample of 372 individuals. Two of them (in FTO and MSRA) were found in lean individuals, nine in extremely obese. In silico analyses of the 11 variants did not reveal functional implications for the mutations. Concordant with our hypothesis we detected a rare variant that potentially leads to loss of FTO function in a lean individual. For TBC1D1, in contrary to our hypothesis, the loss of function variant (Arg443Stop) was found in an obese individual. Functional in vitro studies are warranted.

Fine Mapping of a GWAS-Derived Obesity Candidate Region on Chromosome 16p11.2.

INTRODUCTION: Large-scale genome-wide association studies (GWASs) have identified 97 chromosomal loci associated with increased body mass index in population-based studies on adults. One of these SNPs, rs7359397, tags a large region (approx. 1MB) with high linkage disequilibrium (r2>0.7), which comprises five genes (SH2B1, APOBR, sulfotransferases: SULT1A1 and SULT1A2, TUFM). We had previously described a rare mutation in SH2B1 solely identified in extremely obese individuals but not in lean controls. METHODS: The coding regions of the genes APOBR, SULT1A1, SULT1A2, and TUFM were screened for mutations (dHPLC, SSCP, Sanger re-sequencing) in 95 extremely obese children and adolescents. Detected non-synonymous variants were genotyped (TaqMan SNP Genotyping, MALDI TOF, PCR-RFLP) in independent large study groups (up to 3,210 extremely obese/overweight cases, 485 lean controls and 615 obesity trios). In silico tools were used for the prediction of potential functional effects of detected variants. RESULTS: Except for TUFM we detected non-synonymous variants in all screened genes. Two polymorphisms rs180743 (APOBR p.Pro428Ala) and rs3833080 (APOBR p.Gly369_Asp370del9) showed nominal association to (extreme) obesity (uncorrected p = 0.003 and p = 0.002, respectively). In silico analyses predicted a functional implication for rs180743 (APOBR p.Pro428Ala). Both APOBR variants are located in the repetitive region with unknown function. CONCLUSION: Variants in APOBR contributed as strongly as variants in SH2B1 to the association with extreme obesity in the chromosomal region chr16p11.2. In silico analyses implied no functional effect of several of the detected variants. Further in vitro or in vivo analyses on the functional implications of the obesity associated variants are warranted.

Mitochondrial DNA variants in obesity.

Heritability estimates for body mass index (BMI) variation are high. For mothers and their offspring higher BMI correlations have been described than for fathers. Variation(s) in the exclusively maternally inherited mitochondrial DNA (mtDNA) might contribute to this parental effect. Thirty-two to 40 mtDNA single nucleotide polymorphisms (SNPs) were available from genome-wide association study SNP arrays (Affymetrix 6.0). For discovery, we analyzed association in a case-control (CC) sample of 1,158 extremely obese children and adolescents and 435 lean adult controls. For independent confirmation, 7,014 population-based adults were analyzed as CC sample of n = 1,697 obese cases (BMI ≥ 30 kg/m2) and n = 2,373 normal weight and lean controls (BMI<25 kg/m2). SNPs were analyzed as single SNPs and haplogroups determined by HaploGrep. Fisher's two-sided exact test was used for association testing. Moreover, the D-loop was re-sequenced (Sanger) in 192 extremely obese children and adolescents and 192 lean adult controls. Association testing of detected variants was performed using Fisher's two-sided exact test. For discovery, nominal association with obesity was found for the frequent allele G of m.8994G/A (rs28358887, p = 0.002) located in ATP6. Haplogroup W was nominally overrepresented in the controls (p = 0.039). These findings could not be confirmed independently. For two of the 252 identified D-loop variants nominal association was detected (m.16292C/T, p = 0.007, m.16189T/C, p = 0.048). Only eight controls carried the m.16292T allele, five of whom belonged to haplogroup W that was initially enriched among these controls. m.16189T/C might create an uninterrupted poly-C tract located near a regulatory element involved in replication of mtDNA. Though follow-up of some D-loop variants still is conceivable, our hypothesis of a contribution of variation in the exclusively maternally inherited mtDNA to the observed larger correlations for BMI between mothers and their offspring could not be substantiated by the findings of the present study.

Gene set of nuclear-encoded mitochondrial regulators is enriched for common inherited variation in obesity.

There are hints of an altered mitochondrial function in obesity. Nuclear-encoded genes are relevant for mitochondrial function (3 gene sets of known relevant pathways: (1) 16 nuclear regulators of mitochondrial genes, (2) 91 genes for oxidative phosphorylation and (3) 966 nuclear-encoded mitochondrial genes). Gene set enrichment analysis (GSEA) showed no association with type 2 diabetes mellitus in these gene sets. Here we performed a GSEA for the same gene sets for obesity. Genome wide association study (GWAS) data from a case-control approach on 453 extremely obese children and adolescents and 435 lean adult controls were used for GSEA. For independent confirmation, we analyzed 705 obesity GWAS trios (extremely obese child and both biological parents) and a population-based GWAS sample (KORA F4, n = 1,743). A meta-analysis was performed on all three samples. In each sample, the distribution of significance levels between the respective gene set and those of all genes was compared using the leading-edge-fraction-comparison test (cut-offs between the 50(th) and 95(th) percentile of the set of all gene-wise corrected p-values) as implemented in the MAGENTA software. In the case-control sample, significant enrichment of associations with obesity was observed above the 50(th) percentile for the set of the 16 nuclear regulators of mitochondrial genes (p(GSEA,50) = 0.0103). This finding was not confirmed in the trios (p(GSEA,50) = 0.5991), but in KORA (p(GSEA,50) = 0.0398). The meta-analysis again indicated a trend for enrichment (p(MAGENTA,50) = 0.1052, p(MAGENTA,75) = 0.0251). The GSEA revealed that weak association signals for obesity might be enriched in the gene set of 16 nuclear regulators of mitochondrial genes.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture.

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Mutation screen in the GWAS derived obesity gene SH2B1 including functional analyses of detected variants.

BACKGROUND: The SH2B1 gene (Src-homology 2B adaptor protein 1 gene) is a solid candidate gene for obesity. Large scale GWAS studies depicted markers in the vicinity of the gene; animal models suggest a potential relevance for human body weight regulation. METHODS: We performed a mutation screen for variants in the SH2B1 coding sequence in 95 extremely obese children and adolescents. Detected variants were genotyped in independent childhood and adult study groups (up to 11,406 obese or overweight individuals and 4,568 controls). Functional implications on STAT3 mediated leptin signalling of the detected variants were analyzed in vitro. RESULTS: We identified two new rare mutations and five known SNPs (rs147094247, rs7498665, rs60604881, rs62037368 and rs62037369) in SH2B1. Mutation g.9483C/T leads to a non-synonymous, non-conservative exchange in the beta (ßThr656Ile) and gamma (γPro674Ser) splice variants of SH2B1. It was additionally detected in two of 11,206 (extremely) obese or overweight children, adolescents and adults, but not in 4,506 population-based normal-weight or lean controls. The non-coding mutation g.10182C/A at the 3' end of SH2B1 was only detected in three obese individuals. For the non-synonymous SNP rs7498665 (Thr484Ala) we observed nominal over-transmission of the previously described risk allele in 705 obesity trios (nominal p = 0.009, OR = 1.23) and an increased frequency of the same allele in 359 cases compared to 429 controls (nominal p = 0.042, OR = 1.23). The obesity risk-alleles at Thr484Ala and ßThr656Ile/γPro674Ser had no effect on STAT3 mediated leptin receptor signalling in splice variants ß and γ. CONCLUSION: The rare coding mutation ßThr656Ile/γPro674Ser (g.9483C/T) in SH2B1 was exclusively detected in overweight or obese individuals. Functional analyzes did not reveal impairments in leptin signalling for the mutated SH2B1.

A genome-wide association meta-analysis identifies new childhood obesity loci.

Multiple genetic variants have been associated with adult obesity and a few with severe obesity in childhood; however, less progress has been made in establishing genetic influences on common early-onset obesity. We performed a North American, Australian and European collaborative meta-analysis of 14 studies consisting of 5,530 cases (≥95th percentile of body mass index (BMI)) and 8,318 controls (<50th percentile of BMI) of European ancestry. Taking forward the eight newly discovered signals yielding association with P < 5 × 10(-6) in nine independent data sets (2,818 cases and 4,083 controls), we observed two loci that yielded genome-wide significant combined P values near OLFM4 at 13q14 (rs9568856; P = 1.82 × 10(-9); odds ratio (OR) = 1.22) and within HOXB5 at 17q21 (rs9299; P = 3.54 × 10(-9); OR = 1.14). Both loci continued to show association when two extreme childhood obesity cohorts were included (2,214 cases and 2,674 controls). These two loci also yielded directionally consistent associations in a previous meta-analysis of adult BMI(1).

Genetic variation of the ghrelin activator gene ghrelin O-acyltransferase (GOAT) is associated with anorexia nervosa.

The gastrointestinal peptide hormone ghrelin promotes food intake and increases body weight and adiposity through activation of the growth hormone secretagogue receptor (GHSR1a). To promote its biological action ghrelin is acylated at its serine 3 residue by the recently discovered ghrelin O-acyltransferase (GOAT, a.k.a. membrane-bound O-acyltransferase 4, MBOAT4). Plasma levels of total and acyl-ghrelin are negatively correlated with body-mass-index (BMI); as lower the BMI as higher plasma levels of total and acylated ghrelin and vice versa. Accordingly, plasma levels of total and acyl-ghrelin are elevated in patients with anorexia nervosa (AN) and decline upon weight regain. The importance of the endogenous Goat/ghrelin system in the neuroendocrine adaptation to fasting was recently highlighted by the observation that acyl-ghrelin mediated elevation of growth hormone (GH) release prevents starvation induced hypoglycemia in Goat(-/-) mice. The aim of this study was to test if genetic variation of GOAT is implicated in the etiology of AN. We therefore assessed association of 6 tagging single nucleotide polymorphisms (tagSNPs), which were predicted to cover 96% the common genetic variability of GOAT plus 50 kb of the 5' and 3' flanking region, in 543 German patients with AN and 612 German normal and underweight healthy controls. Based on a recessive mode of inheritance we observed some evidence for association of the G/G genotype at SNP rs10096097 with AN (nominal two-sided p = 0.031). Based on our results we conclude that genetic variation in GOAT might be implicated in the etiology of AN.

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus.

BACKGROUND: Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3' and 5' of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association. CONCLUSIONS/SIGNIFICANCE: A haplotype reaching from a region 5' of the MC4R to a region at least 150 kb from the 3' end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.