Ivonne Jarick

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.85×10−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.84×10−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 ∼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.

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.

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.

Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder. Genetic loci have not yet been identified by genome-wide association studies. Rare copy number variations (CNVs), such as chromosomal deletions or duplications, have been implicated in ADHD and other neurodevelopmental disorders. To identify rare (frequency ⩽1%) CNVs that increase the risk of ADHD, we performed a whole-genome CNV analysis based on 489 young ADHD patients and 1285 adult population-based controls and identified one significantly associated CNV region. In tests for a global burden of large (>500 kb) rare CNVs, we observed a nonsignificant (P=0.271) 1.126-fold enriched rate of subjects carrying at least one such CNV in the group of ADHD cases. Locus-specific tests of association were used to assess if there were more rare CNVs in cases compared with controls. Detected CNVs, which were significantly enriched in the ADHD group, were validated by quantitative (q)PCR. Findings were replicated in an independent sample of 386 young patients with ADHD and 781 young population-based healthy controls. We identified rare CNVs within the parkinson protein 2 gene (PARK2) with a significantly higher prevalence in ADHD patients than in controls (P=2.8 × 10−4 after empirical correction for genome-wide testing). In total, the PARK2 locus (chr 6: 162 659 756–162 767 019) harboured three deletions and nine duplications in the ADHD patients and two deletions and two duplications in the controls. By qPCR analysis, we validated 11 of the 12 CNVs in ADHD patients (P=1.2 × 10−3 after empirical correction for genome-wide testing). In the replication sample, CNVs at the PARK2 locus were found in four additional ADHD patients and one additional control (P=4.3 × 10−2). Our results suggest that copy number variants at the PARK2 locus contribute to the genetic susceptibility of ADHD. Mutations and CNVs in PARK2 are known to be associated with Parkinson disease.

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.

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

BackgroundThe 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.MethodsWe 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.ResultsWe 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 γ.ConclusionThe 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.

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. Fishers 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 Fishers 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 50th and 95th 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 50th percentile for the set of the 16 nuclear regulators of mitochondrial genes (pGSEA,50 = 0.0103). This finding was not confirmed in the trios (pGSEA,50 = 0.5991), but in KORA (pGSEA,50 = 0.0398). The meta-analysis again indicated a trend for enrichment (pMAGENTA,50 = 0.1052, pMAGENTA,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.

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 (r²>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.