Erratum: A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling.

A correction to this article has been published and is linked from the HTML version of this paper. The error has been fixed in the paper.

A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling.

We conducted a genome-wide association study (GWAS) of anorexia nervosa (AN) using a stringently defined phenotype. Analysis of phenotypic variability led to the identification of a specific genetic risk factor that approached genome-wide significance (rs929626 in EBF1 (Early B-Cell Factor 1); P = 2.04 × 10-7; OR = 0.7; 95% confidence interval (CI) = 0.61-0.8) with independent replication (P = 0.04), suggesting a variant-mediated dysregulation of leptin signaling may play a role in AN. Multiple SNPs in LD with the variant support the nominal association. This demonstrates that although the clinical and etiologic heterogeneity of AN is universally recognized, further careful sub-typing of cases may provide more precise genomic signals. In this study, through a refinement of the phenotype spectrum of AN, we present a replicable GWAS signal that is nominally associated with AN, highlighting a potentially important candidate locus for further investigation.

Association of a rare NOTCH4 coding variant with systemic sclerosis: a family-based whole exome sequencing study.

BACKGROUND: Systemic sclerosis (SSc) is a rheumatologic disease with a multifactorial etiology. Genome-wide association studies imply a polygenic, complex mode of inheritance with contributions from variation at the human leukocyte antigen locus and non-coding variation at a locus on chromosome 6p21, among other modestly impactful loci. Here we describe an 8-year-old female proband presenting with diffuse cutaneous SSc/scleroderma and a family history of SSc in a grandfather and maternal aunt. METHODS: We employed whole exome sequencing (WES) of three members of this family. We examined rare missense, nonsense, splice-altering, and coding indels matching an autosomal dominant inheritance model. We selected one missense variant for Sanger sequencing confirmation based on its predicted impact on gene function and location in a known SSc genetic locus. RESULTS: Bioinformatic analysis found eight candidate variants meeting our criteria. We identified a very rare missense variant in the regulatory NODP domain of NOTCH4 located at the 6p21 locus, c.4245G > A:p.Met1415Ile, segregating with the phenotype. This allele has a frequency of 1.83 × 10-5 by the data of the Exome Aggregation Consortium. CONCLUSION: This family suggests a novel mechanism of SSc pathogenesis in which a rare and penetrant coding variation can substantially elevate disease risk in contrast to the more modest non-coding variation typically found at this locus. These results suggest that modulation of the NOTCH4 gene might be responsible for the association signal at chromosome 6p21 in SSc.

Genomic copy number variation association study in Caucasian patients with nonsyndromic cryptorchidism.

BACKGROUND: Copy number variation (CNV) is a potential contributing factor to many genetic diseases. Here we investigated the potential association of CNV with nonsyndromic cryptorchidism, the most common male congenital genitourinary defect, in a Caucasian population. METHODS: Genome wide genotyping were performed in 559 cases and 1772 controls (Group 1) using Illumina HumanHap550 v1, HumanHap550 v3 or Human610-Quad platforms and in 353 cases and 1149 controls (Group 2) using the Illumina Human OmniExpress 12v1 or Human OmniExpress 12v1-1. Signal intensity data including log R ratio (LRR) and B allele frequency (BAF) for each single nucleotide polymorphism (SNP) were used for CNV detection using PennCNV software. After sample quality control, gene- and CNV-based association tests were performed using cleaned data from Group 1 (493 cases and 1586 controls) and Group 2 (307 cases and 1102 controls) using ParseCNV software. Meta-analysis was performed using gene-based test results as input to identify significant genes, and CNVs in or around significant genes were identified in CNV-based association test results. Called CNVs passing quality control and signal intensity visualization examination were considered for validation using TaqMan CNV assays and QuantStudio® 3D Digital PCR System. RESULTS: The meta-analysis identified 373 genome wide significant (p < 5X10-4) genes/loci including 49 genes/loci with deletions and 324 with duplications. Among them, 17 genes with deletion and 1 gene with duplication were identified in CNV-based association results in both Group 1 and Group 2. Only 2 genes (NUCB2 and UPF2) containing deletions passed CNV quality control in both groups and signal intensity visualization examination, but laboratory validation failed to verify these deletions. CONCLUSIONS: Our data do not support that structural variation is a major cause of nonsyndromic cryptorchidism.

Variants in CXCR4 associate with juvenile idiopathic arthritis susceptibility.

BACKGROUND: Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic disease among children, the etiology of which involves a strong genetic component, but much of the underlying genetic determinants still remain unknown. Our aim was to identify novel genetic variants that predispose to JIA. METHODS: We performed a genome-wide association study (GWAS) and replication in a total of 1166 JIA cases and 9500 unrelated controls of European ancestry. Correlation of SNP genotype and gene expression was investigated. Then we conducted targeted resequencing of a candidate locus, among a subset of 480 cases and 480 controls. SUM test was performed to evaluate the association of the identified rare functional variants. RESULTS: The CXCR4 locus on 2q22.1 was found to be significantly associated with JIA, peaking at SNP rs953387. However, this result is subjected to subpopulation stratification within the subjects of European ancestry. After adjusting for principal components, nominal significant association remained (p < 10(-4)). Because of its interesting known function in immune regulation, we carried out further analyses to assess its relationship with JIA. Expression of CXCR4 was correlated with CXCR4 rs953387 genotypes in lymphoblastoid cell lines (p = 0.014) and T-cells (p = 0.0054). In addition, rare non-synonymous and stop-gain sequence variants in CXCR4, putatively damaging for CXCR4 function, were significantly enriched in JIA cases (p = 0.015). CONCLUSION: Our results suggest the association of CXCR4 variants with JIA, implicating that this gene may be involved in the pathogenesis of autoimmune disease. However, because this locus is subjected to population stratification within the subjects of European ancestry, additional replication is still necessary for this locus to be considered a true risk locus for JIA. This cell-surface chemokine receptor has already been targeted in other diseases and may serve as a tractable therapeutic target for a specific subset of pediatric arthritis patients with additional replication and functional validation of the locus.

Genetic sharing and heritability of paediatric age of onset autoimmune diseases.

Autoimmune diseases (AIDs) are polygenic diseases affecting 7-10% of the population in the Western Hemisphere with few effective therapies. Here, we quantify the heritability of paediatric AIDs (pAIDs), including JIA, SLE, CEL, T1D, UC, CD, PS, SPA and CVID, attributable to common genomic variations (SNP-h(2)). SNP-h(2) estimates are most significant for T1D (0.863±s.e. 0.07) and JIA (0.727±s.e. 0.037), more modest for UC (0.386±s.e. 0.04) and CD (0.454±0.025), largely consistent with population estimates and are generally greater than that previously reported by adult GWAS. On pairwise analysis, we observed that the diseases UC-CD (0.69±s.e. 0.07) and JIA-CVID (0.343±s.e. 0.13) are the most strongly correlated. Variations across the MHC strongly contribute to SNP-h(2) in T1D and JIA, but does not significantly contribute to the pairwise rG. Together, our results partition contributions of shared versus disease-specific genomic variations to pAID heritability, identifying pAIDs with unexpected risk sharing, while recapitulating known associations between autoimmune diseases previously reported in adult cohorts.

Meta-analysis of shared genetic architecture across ten pediatric autoimmune diseases.

Genome-wide association studies (GWASs) have identified hundreds of susceptibility genes, including shared associations across clinically distinct autoimmune diseases. We performed an inverse χ(2) meta-analysis across ten pediatric-age-of-onset autoimmune diseases (pAIDs) in a case-control study including more than 6,035 cases and 10,718 shared population-based controls. We identified 27 genome-wide significant loci associated with one or more pAIDs, mapping to in silico-replicated autoimmune-associated genes (including IL2RA) and new candidate loci with established immunoregulatory functions such as ADGRL2, TENM3, ANKRD30A, ADCY7 and CD40LG. The pAID-associated single-nucleotide polymorphisms (SNPs) were functionally enriched for deoxyribonuclease (DNase)-hypersensitivity sites, expression quantitative trait loci (eQTLs), microRNA (miRNA)-binding sites and coding variants. We also identified biologically correlated, pAID-associated candidate gene sets on the basis of immune cell expression profiling and found evidence of genetic sharing. Network and protein-interaction analyses demonstrated converging roles for the signaling pathways of type 1, 2 and 17 helper T cells (TH1, TH2 and TH17), JAK-STAT, interferon and interleukin in multiple autoimmune diseases.

Copy Number Variations in CTNNA3 and RBFOX1 Associate with Pediatric Food Allergy.

Food allergy is a significant public health concern, especially among children. Previous candidate gene studies suggested a few susceptibility loci for food allergy, but no study investigated the contribution of copy number variations (CNVs) to food allergy on a genome-wide scale. To investigate the genetics of food allergy, we performed CNV assessment using high-resolution genome-wide single nucleotide polymorphism arrays. CNV calls from a total of 357 cases with confirmed food allergy and 3980 controls were analyzed within a discovery cohort, followed by a replication analysis composed of 167 cases and 1573 controls. We identified that CNVs in CTNNA3 were significantly associated with food allergy in both the discovery cohort and the replication cohort. Of particular interest, CTNNA3 CNVs hit exons or intron regions rich in histone marker H3K4Me1. CNVs in a second gene (RBFOX1) showed a significant association (p = 7.35 × 10(-5)) with food allergy at the genome-wide level in our meta-analysis of the European ancestry cohorts. The presence of these CNVs was confirmed by quantitative PCR. Furthermore, knockdown of CTNNA3 resulted in upregulation of CD63 and CD203c in mononuclear cells upon PMA stimulation, suggesting a role in sensitization to allergen. We uncovered at least two plausible genes harboring CNV loci that are enriched in pediatric patients with food allergies. The novel gene candidates discovered in this study by genome-wide CNV analysis are compelling drug and diagnostic targets for food allergy.

Pathway analysis supports association of nonsyndromic cryptorchidism with genetic loci linked to cytoskeleton-dependent functions.

STUDY QUESTION: What are the genetic loci that increase susceptibility to nonsyndromic cryptorchidism, or undescended testis? SUMMARY ANSWER: A genome-wide association study (GWAS) suggests that susceptibility to cryptorchidism is heterogeneous, with a subset of suggestive signals linked to cytoskeleton-dependent functions and syndromic forms of the disease. WHAT IS KNOWN ALREADY: Population studies suggest moderate genetic risk of cryptorchidism and possible maternal and environmental contributions to risk. Previous candidate gene analyses have failed to identify a major associated locus, although variants in insulin-like 3 (INSL3), relaxin/insulin-like family peptide receptor 2 (RXFP2) and other hormonal pathway genes may increase risk in a small percentage of patients. STUDY DESIGN, SIZE, DURATION: This is a case-control GWAS of 844 boys with nonsyndromic cryptorchidism and 2718 control subjects without syndromes or genital anomalies, all of European ancestry. PARTICIPANTS/MATERIALS, SETTING, METHODS: All boys with cryptorchidism were diagnosed and treated by a pediatric specialist. In the discovery phase, DNA was extracted from tissue or blood samples and genotyping performed using the Illumina HumanHap550 and Human610-Quad (Group 1) or OmniExpress (Group 2) platform. We imputed genotypes genome-wide, and combined single marker association results in meta-analyses for all cases and for secondary subphenotype analyses based on testis position, laterality and age, and defined genome-wide significance as P = 7 × 10(-9) to correct for multiple testing. Selected markers were genotyped in an independent replication group of European cases (n = 298) and controls (n = 324). We used several bioinformatics tools to analyze top (P < 10(-5)) and suggestive (P < 10(-3)) signals for significant enrichment of signaling pathways, cellular functions and custom gene lists after multiple testing correction. MAIN RESULTS AND THE ROLE OF CHANCE: In the full analysis, we identified 20 top loci, none reaching genome-wide significance, but one passing this threshold in a subphenotype analysis of proximal testis position (rs55867206, near SH3PXD2B, odds ratio = 2.2 (95% confidence interval 1.7, 2.9), P = 2 × 10(-9)). An additional 127 top loci emerged in at least one secondary analysis, particularly of more severe phenotypes. Cytoskeleton-dependent molecular and cellular functions were prevalent in pathway analysis of suggestive signals, and may implicate loci encoding cytoskeletal proteins that participate in androgen receptor signaling. Genes linked to human syndromic cryptorchidism, including hypogonadotropic hypogonadism, and to hormone-responsive and/or differentially expressed genes in normal and cryptorchid rat gubernaculum, were also significantly overrepresented. No tested marker showed significant replication in an independent population. The results suggest heterogeneous, multilocus and potentially multifactorial susceptibility to nonsyndromic cryptorchidism. LIMITATIONS, REASONS FOR CAUTION: The present study failed to identify genome-wide significant markers associated with cryptorchidism that could be replicated in an independent population, so further studies are required to define true positive signals among suggestive loci. WIDER IMPLICATIONS OF THE FINDINGS: As the only GWAS to date of nonsyndromic cryptorchidism, these data will provide a basis for future efforts to understand genetic susceptibility to this common reproductive anomaly and the potential for additive risk from environmental exposures. STUDY FUNDING/COMPETING INTERESTS: This work was supported by R01HD060769 (the Eunice Kennedy Shriver National Institute for Child Health and Human Development (NICHD)), P20RR20173 (the National Center for Research Resources (NCRR), currently P20GM103464 from the National Institute of General Medical Sciences (NIGMS)), an Institute Development Fund to the Center for Applied Genomics at The Children's Hospital of Philadelphia, and Nemours Biomedical Research. The authors have no competing interests to declare.

Mutations in SPECC1L, encoding sperm antigen with calponin homology and coiled-coil domains 1-like, are found in some cases of autosomal dominant Opitz G/BBB syndrome.

BACKGROUND: Opitz G/BBB syndrome is a heterogeneous disorder characterised by variable expression of midline defects including cleft lip and palate, hypertelorism, laryngealtracheoesophageal anomalies, congenital heart defects, and hypospadias. The X-linked form of the condition has been associated with mutations in the MID1 gene on Xp22. The autosomal dominant form has been linked to chromosome 22q11.2, although the causative gene has yet to be elucidated. METHODS AND RESULTS: In this study, we performed whole exome sequencing on DNA samples from a three-generation family with characteristics of Opitz G/BBB syndrome with negative MID1 sequencing. We identified a heterozygous missense mutation c.1189A>C (p.Thr397Pro) in SPECC1L, located at chromosome 22q11.23. Mutation screening of an additional 19 patients with features of autosomal dominant Opitz G/BBB syndrome identified a c.3247G>A (p.Gly1083Ser) mutation segregating with the phenotype in another three-generation family. CONCLUSIONS: Previously, SPECC1L was shown to be required for proper facial morphogenesis with disruptions identified in two patients with oblique facial clefts. Collectively, these data demonstrate that SPECC1L mutations can cause syndromic forms of facial clefting including some cases of autosomal dominant Opitz G/BBB syndrome and support the original linkage to chromosome 22q11.2.

Phenotype specific association of the TGFBR3 locus with nonsyndromic cryptorchidism.

PURPOSE: Based on a genome-wide association study of testicular dysgenesis syndrome showing a possible association with TGFBR3, we analyzed data from a larger, phenotypically restricted cryptorchidism population for potential replication of this signal. MATERIALS AND METHODS: We excluded samples based on strict quality control criteria, leaving 844 cases and 2,718 controls of European ancestry that were analyzed in 2 separate groups based on genotyping platform (ie Illumina® HumanHap550, version 1 or 3, or Human610-Quad, version 1 BeadChip in group 1 and Human OmniExpress 12, version 1 BeadChip platform in group 2). Analyses included genotype imputation at the TGFBR3 locus, association analysis of imputed data with correction for population substructure, subsequent meta-analysis of data for groups 1 and 2, and selective genotyping of independent cases (330) and controls (324) for replication. We also measured Tgfbr3 mRNA levels and performed TGFBR3/betaglycan immunostaining in rat fetal gubernaculum. RESULTS: We identified suggestive (p ≤ 1× 10(-4)) association of markers in/near TGFBR3, including rs9661103 (OR 1.40; 95% CI 1.20, 1.64; p = 2.71 × 10(-5)) and rs10782968 (OR 1.58; 95% CI 1.26, 1.98; p = 9.36 × 10(-5)) in groups 1 and 2, respectively. In subgroup analyses we observed strongest association of rs17576372 (OR 1.42; 95% CI 1.24, 1.60; p = 1.67 × 10(-4)) with proximal and rs11165059 (OR 1.32; 95% CI 1.15, 1.38; p = 9.42 × 10(-4)) with distal testis position, signals in strong linkage disequilibrium with rs9661103 and rs10782968, respectively. Association of the prior genome-wide association study signal (rs12082710) was marginal (OR 1.13; 95% CI 0.99, 1.28; p = 0.09 for group 1), and we were unable to replicate signals in our independent cohort. Tgfbr3/betaglycan was differentially expressed in wild-type and cryptorchid rat fetal gubernaculum. CONCLUSIONS: These data suggest complex or phenotype specific association of cryptorchidism with TGFBR3 and the gubernaculum as a potential target of TGFß signaling.

The missense variation landscape of FTO, MC4R, and TMEM18 in obese children of African Ancestry.

OBJECTIVE: Common variation at the loci harboring fat mass and obesity (FTO), melanocortin receptor 4 (MC4R), and transmembrane protein 18 (TMEM18) is consistently reported as being statistically most strongly associated with obesity. Investigations if these loci also harbor rarer missense variants that confer substantially higher risk of common childhood obesity in African American (AA) children were conducted. DESIGN AND METHODS: The exons of FTO, MC4R, and TMEM18 in an initial subset of our cohort were sequenced, that is, 200 obese (BMI ≥ 95 th percentile) and 200 lean AA children (BMI ≤ 5 th percentile). Any missense exonic variants that were uncovered went on to be further genotyped in a further 768 obese and 768 lean (BMI≤50th percentile) children of the same ethnicity. RESULTS: A number of exonic variants were observed from our sequencing effort: seven in FTO, of which four were non-synonymous (A163T, G182A, M400V, and A405V), thirteen in MC4R, of which six were non-synonymous (V103I, N123S, S136A, F202L, N240S, and I251L), and four in TMEM18, of which two were non-synonymous (P2S and V113L). Follow-up genotyping of these missense variants revealed only one significant difference in allele frequency between cases and controls, namely with N240S in MC4R (Fisher's exact P = 0.0001). CONCLUSION: In summary, moderately rare missense variants within the FTO, MC4R, and TMEM18 genes observed in our study did not confer risk of common childhood obesity in African Americans except for a degree of evidence for one known loss-of-function variant in MC4R.

GWAS of blood cell traits identifies novel associated loci and epistatic interactions in Caucasian and African-American children.

Hematological traits are important clinical indicators, the genetic determinants of which have not been fully investigated. Common measures of hematological traits include red blood cell (RBC) count, hemoglobin concentration (HGB), hematocrit (HCT), mean corpuscular hemoglobin (MCH), MCH concentration (MCHC), mean corpuscular volume (MCV), platelet count (PLT) and white blood cell (WBC) count. We carried out a genome-wide association study of the eight common hematological traits among 7943 African-American children and 6234 Caucasian children. In African Americans, we report five novel associations of HBE1 variants with HCT and MCHC, the alpha-globin gene cluster variants with RBC and MCHC, and a variant at the ARHGEF3 locus with PLT, as well as replication of four previously reported loci at genome-wide significance. In Caucasians, we report a novel association of variants at the COPZ1 locus with PLT as well as replication of four previously reported loci at genome-wide significance. Extended analysis of an association observed between MCH and the alpha-globin gene cluster variants demonstrated independent effects and epistatic interaction at the locus, impacting the risk of iron deficiency anemia in African Americans with specific genotype states. In summary, we extend the understanding of genetic variants underlying hematological traits based on analyses in African-American children.

Examination of genetic variants influencing lipid traits in pediatric populations.

Previous large-scale genome-wide association studies in adult populations have implicated ∽100 loci in determining high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol, or triglyceride levels. However, whether these loci also contribute to variations of lipid traits in pediatric populations remain unknown. Here we assayed a population of Philadelphia children by high-density single nucleotide polymorphism arrays, and performed association analysis on lipid traits ascertained from lipid measurements stored in electronic medical records. We examined previously reported lipid trait associations, and found that most of them show identical direction of association in our pediatric cohorts, including genome-wide significant association on cholesteryl ester transfer protein with HDL-C levels (rs3764261, P = 2.1 × 10(-8)) and other significant associations on oxysterol-binding protein-like protein 7, low-density lipoprotein receptor-related protein 4 and low-density lipoprotein receptor-related protein 1. Additionally, we identified suggestive association on low-density lipoprotein receptor-related protein 1B with HDL-C levels (rs17736712, P = 2.1 × 10(-7)), but this signal is not supported by previous meta-analysis on adult cohorts. Finally, we examined rare copy number variants and identified deletions encompassing tetratricopeptide repeat domain 39B in two children with extreme lipid measures. Our results highlight the commonalities and differences of genetic components in determining lipid traits in pediatric versus adult populations. Furthermore, our study demonstrates the unique utility of automated information retrieval from electronic medical records in facilitating the identification of genotype-phenotype associations.

Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder.

Attention deficit hyperactivity disorder (ADHD) is a common, heritable neuropsychiatric disorder of unknown etiology. We performed a whole-genome copy number variation (CNV) study on 1,013 cases with ADHD and 4,105 healthy children of European ancestry using 550,000 SNPs. We evaluated statistically significant findings in multiple independent cohorts, with a total of 2,493 cases with ADHD and 9,222 controls of European ancestry, using matched platforms. CNVs affecting metabotropic glutamate receptor genes were enriched across all cohorts (P = 2.1 × 10(-9)). We saw GRM5 (encoding glutamate receptor, metabotropic 5) deletions in ten cases and one control (P = 1.36 × 10(-6)). We saw GRM7 deletions in six cases, and we saw GRM8 deletions in eight cases and no controls. GRM1 was duplicated in eight cases. We experimentally validated the observed variants using quantitative RT-PCR. A gene network analysis showed that genes interacting with the genes in the GRM family are enriched for CNVs in ∼10% of the cases (P = 4.38 × 10(-10)) after correction for occurrence in the controls. We identified rare recurrent CNVs affecting glutamatergic neurotransmission genes that were overrepresented in multiple ADHD cohorts.

A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci.

Diabetes impacts approximately 200 million people worldwide, of whom approximately 10% are affected by type 1 diabetes (T1D). The application of genome-wide association studies (GWAS) has robustly revealed dozens of genetic contributors to the pathogenesis of T1D, with the most recent meta-analysis identifying in excess of 40 loci. To identify additional genetic loci for T1D susceptibility, we examined associations in the largest meta-analysis to date between the disease and ∼2.54 million SNPs in a combined cohort of 9,934 cases and 16,956 controls. Targeted follow-up of 53 SNPs in 1,120 affected trios uncovered three new loci associated with T1D that reached genome-wide significance. The most significantly associated SNP (rs539514, P = 5.66×10⁻¹¹) resides in an intronic region of the LMO7 (LIM domain only 7) gene on 13q22. The second most significantly associated SNP (rs478222, P = 3.50×10⁻9 resides in an intronic region of the EFR3B (protein EFR3 homolog B) gene on 2p23; however, the region of linkage disequilibrium is approximately 800 kb and harbors additional multiple genes, including NCOA1, C2orf79, CENPO, ADCY3, DNAJC27, POMC, and DNMT3A. The third most significantly associated SNP (rs924043, P = 8.06×10⁻9 lies in an intergenic region on 6q27, where the region of association is approximately 900 kb and harbors multiple genes including WDR27, C6orf120, PHF10, TCTE3, C6orf208, LOC154449, DLL1, FAM120B, PSMB1, TBP, and PCD2. These latest associated regions add to the growing repertoire of gene networks predisposing to T1D.

Role of BMI-associated loci identified in GWAS meta-analyses in the context of common childhood obesity in European Americans.

Obesity is a serious health concern for children and adolescents, particularly in Western societies, where its incidence is now considered to have reached epidemic proportions. A number of genetic determinants of adult BMI have already been established through genome wide association studies (GWAS), most recently from the GIANT meta-analysis of such datasets combined. In this current study of European Americans, we examined the 32 loci detected in that GIANT study in the context of common childhood obesity within a cohort of 1,097 cases (defined as BMI ≥95th percentile), together with 2,760 lean controls (defined as BMI <50th percentile), aged between 2 and 18 years old. Nine of these single-nucleotide polymorphims (SNPs) yielded at least nominal evidence for association with common childhood obesity, namely at the FTO, TMEM18, NRXN3, MC4R, SEC16B, GNPDA2, TNNI3K, QPCTL, and BDNF loci. However, overall 28 of the 32 loci showed directionally consistent effects to that of the adult BMI meta-analysis. We conclude that among the 32 loci that have been reported to associate with adult BMI in the largest meta-analysis of BMI to date, at least nine also contribute to the determination of common obesity in childhood in European Americans, as demonstrated by their associations in our pediatric cohort.

BMD-associated variation at the Osterix locus is correlated with childhood obesity in females.

Recent genome wide association studies (GWAS) have revealed a number of genetic variants robustly associated with bone mineral density (BMD) and/or osteoporosis. Evidence from epidemiological and clinical studies has shown an association between BMD and BMI, presumably as a consequence of bone loading. We investigated the 23 previously published BMD GWAS-derived loci in the context of childhood obesity by leveraging our existing genome-wide genotyped European American cohort of 1106 obese children (BMI ≥ 95th percentile) and 5997 controls (BMI < 95th percentile). Evidence of association was only observed at one locus, namely Osterix (SP7), with the G allele of rs2016266 being significantly over-represented among childhood obesity cases (P = 2.85 × 10(-3)). When restricting these analyses to each gender, we observed strong association between rs2016266 and childhood obesity in females (477 cases and 2867 controls; P = 3.56 × 10(-4)), which survived adjustment for all tests applied. However, no evidence of association was observed among males. Interestingly, Osterix is the only GWAS locus uncovered to date that has also been previously implicated in the determination of BMD in childhood. In conclusion, these findings indicate that a well established variant at the Osterix locus associated with increased BMD is also associated with childhood obesity primarily in females.

A genome-wide study reveals copy number variants exclusive to childhood obesity cases.

The prevalence of obesity in children and adults in the United States has increased dramatically over the past decade. Genomic copy number variations (CNVs) have been strongly implicated in subjects with extreme obesity and coexisting developmental delay. To complement these previous studies, we addressed CNVs in common childhood obesity by examining children with a BMI in the upper 5(th) percentile but excluding any subject greater than three standard deviations from the mean in order to reduce severe cases in the cohort. We performed a whole-genome CNV survey of our cohort of 1080 defined European American (EA) childhood obesity cases and 2500 lean controls (< 50(th) percentile BMI) who were genotyped with 550,000 SNP markers. Positive findings were evaluated in an independent African American (AA) cohort of 1479 childhood obesity cases and 1575 lean controls. We identified 17 CNV loci that were unique to at least three EA cases and were both previously unreported in the public domain and validated via quantitative PCR. Eight of these loci (47.1%) also replicated exclusively in AA cases (six deletions and two duplications). Replicated deletion loci consisted of EDIL3, S1PR5, FOXP2, TBCA, ABCB5, and ZPLD1, whereas replicated duplication loci consisted of KIF2B and ARL15. We also observed evidence for a deletion at the EPHA6-UNQ6114 locus when the AA cohort was investigated as a discovery set. Although these variants may be individually rare, our results indicate that CNVs contribute to the genetic susceptibility of common childhood obesity in subjects of both European and African ancestry.

The role of height-associated loci identified in genome wide association studies in the determination of pediatric stature.

BACKGROUND: Human height is considered highly heritable and correlated with certain disorders, such as type 2 diabetes and cancer. Despite environmental influences, genetic factors are known to play an important role in stature determination. A number of genetic determinants of adult height have already been established through genome wide association studies. METHODS: To examine 51 single nucleotide polymorphisms (SNPs) corresponding to the 46 previously reported genomic loci for height in 8,184 European American children with height measurements. We leveraged genotyping data from our ongoing GWA study of height variation in children in order to query the 51 SNPs in this pediatric cohort. RESULTS: Sixteen of these SNPs yielded at least nominally significant association to height, representing fifteen different loci including EFEMP1-PNPT1, GPR126, C6orf173, SPAG17, Histone class 1, HLA class III and GDF5-UQCC. Other loci revealed no evidence for association, including HMGA1 and HMGA2. For the 16 associated variants, the genotype score explained 1.64% of the total variation for height z-score. CONCLUSION: Among 46 loci that have been reported to associate with adult height to date, at least 15 also contribute to the determination of height in childhood.