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.

Variants of DENND1B associated with asthma in children.

BACKGROUND: Asthma is a complex disease that has genetic and environmental causes. The genetic factors associated with susceptibility to asthma remain largely unknown. METHODS: We carried out a genomewide association study involving children with asthma. The sample included 793 North American children of European ancestry with persistent asthma who required daily inhaled glucocorticoid therapy and 1988 matched controls (the discovery set). We also tested for genomewide association in an independent cohort of 917 persons of European ancestry who had asthma and 1546 matched controls (the replication set). Finally, we tested for an association between 20 single-nucleotide polymorphisms (SNPs) at chromosome 1q31 and asthma in 1667 North American children of African ancestry who had asthma and 2045 ancestrally matched controls. RESULTS: In our meta-analysis of all samples from persons of European ancestry, we observed an association, with genomewide significance, between asthma and SNPs at the previously reported locus on 17q21 and an additional eight SNPs at a novel locus on 1q31. The SNP most strongly associated with asthma was rs2786098 (P=8.55x10(-9)). We observed replication of the association of asthma with SNP rs2786098 in the independent series of persons of European ancestry (combined P=9.3x10(-11)). The alternative allele of each of the eight SNPs on chromosome 1q31 was strongly associated with asthma in the children of African ancestry (P=1.6x10(-13) for the comparison across all samples). The 1q31 locus contains the 1q31 locus contains DENND1B, a gene expressed by natural killer cells and dendritic cells. DENND1B protein is predicted to interact with the tumor necrosis factor α receptor [corrected]. CONCLUSIONS: We have identified a locus containing DENND1B on chromosome 1q31.3 that is associated with susceptibility to asthma.

Rare recurrent copy number variations in metabotropic glutamate receptor interacting genes in children with neurodevelopmental disorders.

BACKGROUND: Neurodevelopmental disorders (NDDs), such as attention deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD), are examples of complex and partially overlapping phenotypes that often lack definitive corroborating genetic information. ADHD and ASD have complex genetic associations implicated by rare recurrent copy number variations (CNVs). Both of these NDDs have been shown to share similar biological etiologies as well as genetic pleiotropy. METHODS: Platforms aimed at investigating genetic-based associations, such as high-density microarray technologies, have been groundbreaking techniques in the field of complex diseases, aimed at elucidating the underlying disease biology. Previous studies have uncovered CNVs associated with genes within shared candidate genomic networks, including glutamate receptor genes, across multiple different NDDs. To examine shared biological pathways across two of the most common NDDs, we investigated CNVs across 15,689 individuals with ADHD (n = 7920), ASD (n = 4318), or both (n = 3,416), as well as 19,993 controls. Cases and controls were matched by genotype array (i.e., Illumina array versions). Three case-control association studies each calculated and compared the observed vs. expected frequency of CNVs across individual genes, loci, pathways, and gene networks. Quality control measures of confidence in CNV-calling, prior to association analyses, included visual inspection of genotype and hybridization intensity. RESULTS: Here, we report results from CNV analysis in search for individual genes, loci, pathways, and gene networks. To extend our previous observations implicating a key role of the metabotropic glutamate receptor (mGluR) network in both ADHD and autism, we exhaustively queried patients with ASD and/or ADHD for CNVs associated with the 273 genomic regions of interest within the mGluR gene network (genes with one or two degrees protein-protein interaction with mGluR 1-8 genes). Among CNVs in mGluR network genes, we uncovered CNTN4 deletions enriched in NDD cases (P = 3.22E - 26, OR = 2.49). Additionally, we uncovered PRLHR deletions in 40 ADHD cases and 12 controls (P = 5.26E - 13, OR = 8.45) as well as clinically diagnostic relevant 22q11.2 duplications and 16p11.2 duplications in 23 ADHD + ASD cases and 9 controls (P = 4.08E - 13, OR = 15.05) and 22q11.2 duplications in 34 ADHD + ASD cases and 51 controls (P = 9.21E - 9, OR = 3.93); those control samples were not with previous 22qDS diagnosis in their EHR records. CONCLUSION: Together, these results suggest that disruption in neuronal cell-adhesion pathways confers significant risk to NDDs and showcase that rare recurrent CNVs in CNTN4, 22q11.2, and 16p11.2 are overrepresented in NDDs that constitute patients predominantly suffering from ADHD and ASD. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT02286817 First Posted: 10 November 14, ClinicalTrials.gov Identifier: NCT02777931 first posted: 19 May 2016, ClinicalTrials.gov Identifier: NCT03006367 first posted: 30 December 2016, ClinicalTrials.gov Identifier: NCT02895906 first posted: 12 September 2016.

High-resolution mapping and analysis of copy number variations in the human genome: a data resource for clinical and research applications.

We present a database of copy number variations (CNVs) detected in 2026 disease-free individuals, using high-density, SNP-based oligonucleotide microarrays. This large cohort, comprised mainly of Caucasians (65.2%) and African-Americans (34.2%), was analyzed for CNVs in a single study using a uniform array platform and computational process. We have catalogued and characterized 54,462 individual CNVs, 77.8% of which were identified in multiple unrelated individuals. These nonunique CNVs mapped to 3272 distinct regions of genomic variation spanning 5.9% of the genome; 51.5% of these were previously unreported, and >85% are rare. Our annotation and analysis confirmed and extended previously reported correlations between CNVs and several genomic features such as repetitive DNA elements, segmental duplications, and genes. We demonstrate the utility of this data set in distinguishing CNVs with pathologic significance from normal variants. Together, this analysis and annotation provides a useful resource to assist with the assessment of CNVs in the contexts of human variation, disease susceptibility, and clinical molecular diagnostics.

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.

A genome-wide association study identifies a locus for nonsyndromic cleft lip with or without cleft palate on 8q24.

OBJECTIVE: To identify, in a non-hypothesis manner, novel genetic factors associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P). STUDY DESIGN: We performed a genome-wide association study in a pediatric cohort of European decent consisting of 111 NSCL/P cases and 5951 control subjects. All subjects were consecutively recruited from the Greater Philadelphia area from 2006 to 2009. High throughput genome-wide single nucleotide polymorphism genotyping was carried out with the Illumina Infinium II HumanHap550 BeadChip technology. RESULTS: We observed association at the genome-wide significance level with SNP rs987525 at a locus on 8q24, which harbors no characterized genes to date (P = 9.18 x 10(-8); odds ratio = 2.09, 95% confidence interval = 1.59 to 2.76). While searching for a replication cohort, the same genetic determinant was established through a genome-wide association study of NSCL/P in Germany, so this previous report acts as a de novo replication for our independent observation outlined here. CONCLUSIONS: These results strongly suggest that a locus on 8q24 is involved in the pathogenesis of NSCL/P.

Association of variants of the interleukin-23 receptor gene with susceptibility to pediatric Crohn's disease.

BACKGROUND & AIMS: Recently an association was shown between the single nucleotide polymorphism (SNP), rs11209026, within the interleukin-23 receptor (IL23R) locus and Crohn's disease (CD) as a consequence of a genome-wide association study of this disease in adults. We examined the effects of this and other previously reported SNPs at this locus with respect to CD in children. METHODS: By using data from our ongoing genome-wide association study in our cohort of 142 pediatric CD patients and 281 matched controls, we investigated the association of the previously reported SNPs at the IL23R locus with the childhood form of this disease. RESULTS: By using the Fisher exact test, the minor allele frequency of rs11209026 in the patients was 1.75%, whereas it was 6.61% in the controls, yielding a protective odds ratio (OR) of 0.25 (95% confidence interval, 0.10-0.65; 1-sided P = 9.2 x 10(-4)). Furthermore, of all the SNPs previously reported, rs11209026 was associated the most strongly. A subsequent family based association test (which is more resistant to population stratification) with 65 sets of trios derived from our initial patient cohort yielded significant association with rs11209026 in a transmission disequilibrium test (1-sided P = .0017). In contrast, no association was detected to the caspase-recruitment domain 15 gene for the inflammatory bowel disease phenotype. CONCLUSIONS: The OR of the IL23R variant in our pediatric study is highly comparable with that reported previously in a non-Jewish adult inflammatory bowel disease case-control cohort (OR, 0.26). As such, variants in the IL23R gene confer a similar magnitude of risk of CD to children as for their adult counterparts.

Examination of all type 2 diabetes GWAS loci reveals HHEX-IDE as a locus influencing pediatric BMI.

OBJECTIVE: A number of studies have found that BMI in early life influences the risk of developing type 2 diabetes later in life. Our goal was to investigate if any type 2 diabetes variants uncovered through genome-wide association studies (GWAS) impact BMI in childhood. RESEARCH DESIGN AND METHODS: Using data from an ongoing GWAS of pediatric BMI in our cohort, we investigated the association of pediatric BMI with 20 single nucleotide polymorphisms at 18 type 2 diabetes loci uncovered through GWAS, consisting of ADAMTS9, CDC123-CAMK1D, CDKAL1, CDKN2A/B, EXT2, FTO, HHEX-IDE, IGF2BP2, the intragenic region on 11p12, JAZF1, KCNQ1, LOC387761, MTNR1B, NOTCH2, SLC30A8, TCF7L2, THADA, and TSPAN8-LGR5. We randomly partitioned our cohort exactly in half in order to have a discovery cohort (n = 3,592) and a replication cohort (n = 3,592). RESULTS: Our data show that the major type 2 diabetes risk-conferring G allele of rs7923837 at the HHEX-IDE locus was associated with higher pediatric BMI in both the discovery (P = 0.0013 and survived correction for 20 tests) and replication (P = 0.023) sets (combined P = 1.01 x 10(-4)). Association was not detected with any other known type 2 diabetes loci uncovered to date through GWAS except for the well-established FTO. CONCLUSIONS: Our data show that the same genetic HHEX-IDE variant, which is associated with type 2 diabetes from previous studies, also influences pediatric BMI.

Investigation of the locus near MC4R with childhood obesity in Americans of European and African ancestry.

Recently a modest, but consistently, replicated association was demonstrated between obesity and the single-nucleotide polymorphism (SNP), rs17782313, 3' of the MC4R locus as a consequence of a meta-analysis of genome-wide association (GWA) studies of the disease in white populations. We investigated the association in the context of the childhood form of the disease utilizing data from our ongoing GWA study in a cohort of 728 European-American (EA) obese children (BMI > or =95th percentile) and 3,960 EA controls (BMI <95th percentile), as well as 1,008 African-American (AA) obese children and 2,715 AA controls. rs571312, rs10871777, and rs476828 (perfect surrogates for rs17782313) yielded odds ratios in the EA cohort of 1.142 (P = 0.045), 1.137 (P = 0.054), and 1.145 (P = 0.042); however, there was no significant association with these SNPs in the AA cohort. When investigating all 30 SNPs present on the Illumina BeadChip at this locus, again there was no evidence for association in AA cases when correcting for the number of tests employed. As such, variants 3' to the MC4R locus present on the genotyping platform utilized confer a similar magnitude of risk of obesity in white children as to their adult white counterparts but this observation did not extend to AAs.

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

The prevalence of obesity in children and adults in the United States has increased dramatically over the past decade. Besides environmental factors, genetic factors are known to play an important role in the pathogenesis of obesity. A number of genetic determinants of adult BMI have already been established through genome-wide association (GWA) studies. In this study, we examined 25 single-nucleotide polymorphisms (SNPs) corresponding to 13 previously reported genomic loci in 6,078 children with measures of BMI. Fifteen of these SNPs yielded at least nominally significant association to BMI, representing nine different loci including INSIG2, FTO, MC4R, TMEM18, GNPDA2, NEGR1, BDNF, KCTD15, and 1q25. Other loci revealed no evidence for association, namely at MTCH2, SH2B1, 12q13, and 3q27. For the 15 associated variants, the genotype score explained 1.12% of the total variation for BMI z-score. We conclude that among 13 loci that have been reported to associate with adult BMI, at least nine also contribute to the determination of BMI in childhood as demonstrated by their associations in our pediatric cohort.

Association of the BANK 1 R61H variant with systemic lupus erythematosus in Americans of European and African ancestry.

Recently an association was demonstrated between the single nucleotide polymorphism (SNP), rs10516487, within the B-cell gene BANK1 and systemic lupus erythematosus (SLE) as a consequence of a genome wide association study of this disease in European and Argentinean populations. In a bid for replication, we examined the effects of the R61H non-synonymous variant with respect to SLE in our genotyped American cohorts of European and African ancestry. Utilizing data from our ongoing genome-wide association study in our cohort of 178 Caucasian SLE cases and 1808 Caucasian population-based controls plus 148 African American (AA) SLE cases and 1894 AA population-based controls we investigated the association of the previously described non-synonymous SNP at the BANK1 locus with the disease in the two ethnicities separately. Using a Fisher's exact test, the minor allele frequency (MAF) of rs10516487 in the Caucasian cases was 22.6% while it was 31.2% in Caucasian controls, yielding a protective odds ratio (OR) of 0.64 (95% CI 0.49-0.85; one-sided p = 7.07 × 10(-4)). Furthermore, the MAF of rs10516487 in the AA cases was 18.7% while it was 23.3% in AA controls, yielding a protective OR of 0.75 (95% CI 0.55-1.034; one-sided p = 0.039). The OR of the BANK1 variant in our study cohorts is highly comparable with that reported previously in a South American/European SLE case-control cohort (OR = 0.72). As such, R61H in the BANK1 gene confers a similar magnitude of SLE protection, not only in European Americans, but also in African Americans.

Examination of type 2 diabetes loci implicates CDKAL1 as a birth weight gene.

OBJECTIVE: A number of studies have found that reduced birth weight is associated with type 2 diabetes later in life; however, the underlying mechanism for this correlation remains unresolved. Recently, association has been demonstrated between low birth weight and single nucleotide polymorphisms (SNPs) at the CDKAL1 and HHEX-IDE loci, regions that were previously implicated in the pathogenesis of type 2 diabetes. In order to investigate whether type 2 diabetes risk-conferring alleles associate with low birth weight in our Caucasian childhood cohort, we examined the effects of 20 such loci on this trait. RESEARCH DESIGN AND METHODS: Using data from an ongoing genome-wide association study in our cohort of 5,465 Caucasian children with recorded birth weights, we investigated the association of the previously reported type 2 diabetes-associated variation at 20 loci including TCF7L2, HHEX-IDE, PPARG, KCNJ11, SLC30A8, IGF2BP2, CDKAL1, CDKN2A/2B, and JAZF1 with birth weight. RESULTS: Our data show that the minor allele of rs7756992 (P = 8 x 10(-5)) at the CDKAL1 locus is strongly associated with lower birth weight, whereas a perfect surrogate for variation previously implicated for the trait at the same locus only yielded nominally significant association (P = 0.01; r(2) rs7756992 = 0.677). However, association was not detected with any of the other type 2 diabetes loci studied. CONCLUSIONS: We observe association between lower birth weight and type 2 diabetes risk-conferring alleles at the CDKAL1 locus. Our data show that the same genetic locus that has been identified as a marker for type 2 diabetes in previous studies also influences birth weight.

Follow-up analysis of genome-wide association data identifies novel loci for type 1 diabetes.

OBJECTIVE: Two recent genome-wide association (GWA) studies have revealed novel loci for type 1 diabetes, a common multifactorial disease with a strong genetic component. To fully utilize the GWA data that we had obtained by genotyping 563 type 1 diabetes probands and 1,146 control subjects, as well as 483 case subject-parent trios, using the Illumina HumanHap550 BeadChip, we designed a full stage 2 study to capture other possible association signals. RESEARCH DESIGN AND METHODS: From our existing datasets, we selected 982 markers with P < 0.05 in both GWA cohorts. Genotyping these in an independent set of 636 nuclear families with 974 affected offspring revealed 75 markers that also had P < 0.05 in this third cohort. Among these, six single nucleotide polymorphisms in five novel loci also had P < 0.05 in the Wellcome Trust Case-Control Consortium dataset and were further tested in 1,303 type 1 diabetes probands from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) plus 1,673 control subjects. RESULTS: Two markers (rs9976767 and rs3757247) remained significant after adjusting for the number of tests in this last cohort; they reside in UBASH3A (OR 1.16; combined P = 2.33 x 10(-8)) and BACH2 (1.13; combined P = 1.25 x 10(-6)). CONCLUSIONS: Evaluation of a large number of statistical GWA candidates in several independent cohorts has revealed additional loci that are associated with type 1 diabetes. The two genes at these respective loci, UBASH3A and BACH2, are both biologically relevant to autoimmunity.

Association analysis of the FTO gene with obesity in children of Caucasian and African ancestry reveals a common tagging SNP.

Recently an association was demonstrated between the single nucleotide polymorphism (SNP), rs9939609, within the FTO locus and obesity as a consequence of a genome wide association (GWA) study of type 2 diabetes in adults. We examined the effects of two perfect surrogates for this SNP plus 11 other SNPs at this locus with respect to our childhood obesity cohort, consisting of both Caucasians and African Americans (AA). Utilizing data from our ongoing GWA study in our cohort of 418 Caucasian obese children (BMI>or=95th percentile), 2,270 Caucasian controls (BMI<95th percentile), 578 AA obese children and 1,424 AA controls, we investigated the association of the previously reported variation at the FTO locus with the childhood form of this disease in both ethnicities. The minor allele frequencies (MAF) of rs8050136 and rs3751812 (perfect surrogates for rs9939609 i.e. both r(2) = 1) in the Caucasian cases were 0.448 and 0.443 respectively while they were 0.391 and 0.386 in Caucasian controls respectively, yielding for both an odds ratio (OR) of 1.27 (95% CI 1.08-1.47; P = 0.0022). Furthermore, the MAFs of rs8050136 and rs3751812 in the AA cases were 0.449 and 0.115 respectively while they were 0.436 and 0.090 in AA controls respectively, yielding an OR of 1.05 (95% CI 0.91-1.21; P = 0.49) and of 1.31 (95% CI 1.050-1.643; P = 0.017) respectively. Investigating all 13 SNPs present on the Illumina HumanHap550 BeadChip in this region of linkage disequilibrium, rs3751812 was the only SNP conferring significant risk in AA. We have therefore replicated and refined the association in an AA cohort and distilled a tag-SNP, rs3751812, which captures the ancestral origin of the actual mutation. As such, variants in the FTO gene confer a similar magnitude of risk of obesity to children as to their adult counterparts and appear to have a global impact.

A novel susceptibility locus for type 1 diabetes on Chr12q13 identified by a genome-wide association study.

OBJECTIVE: In stage 1 of our genome-wide association (GWA) study for type 1 diabetes, one locus at 16p13 was detected (P = 1.03 x 10(-10)) and confirmed in two additional cohorts. Here we describe the results of testing, in these additional cohorts, 23 loci that were next in rank of statistical significance. RESEARCH DESIGN AND METHODS: Two independent cohorts were studied. The Type 1 Diabetes Genetics Consortium replication cohort consisted of 549 families with at least one child diagnosed with diabetes (946 total affected) and DNA from both parents. The Canadian replication cohort consisted of 364 nuclear family trios with one type 1 diabetes-affected offspring and two parents (1,092 individuals). RESULTS: One locus at 12q13, with the highest statistical significance among the 23, was confirmed. It involves type 1 diabetes association with the minor allele of rs1701704 (P = 9.13 x 10(-10), OR 1.25 [95% CI 1.12-1.40]). CONCLUSIONS: We have discovered a type 1 diabetes locus at 12q13 that is replicated in an independent cohort of type 1 diabetic patients and confers a type 1 diabetes risk comparable with that of the 16p13 locus we recently reported. These two loci are identical to two loci identified by the whole-genome association study of the Wellcome Trust Case-Control Consortium, a parallel independent discovery that adds further support to the validity of the GWA approach.

Loci on 20q13 and 21q22 are associated with pediatric-onset inflammatory bowel disease.

Inflammatory bowel disease (IBD) is a common inflammatory disorder with complex etiology that involves both genetic and environmental triggers, including but not limited to defects in bacterial clearance, defective mucosal barrier and persistent dysregulation of the immune response to commensal intestinal bacteria. IBD is characterized by two distinct phenotypes: Crohn's disease (CD) and ulcerative colitis (UC). Previously reported GWA studies have identified genetic variation accounting for a small portion of the overall genetic susceptibility to CD and an even smaller contribution to UC pathogenesis. We hypothesized that stratification of IBD by age of onset might identify additional genes associated with IBD. To that end, we carried out a GWA analysis in a cohort of 1,011 individuals with pediatric-onset IBD and 4,250 matched controls. We identified and replicated significantly associated, previously unreported loci on chromosomes 20q13 (rs2315008[T] and rs4809330[A]; P = 6.30 x 10(-8) and 6.95 x 10(-8), respectively; odds ratio (OR) = 0.74 for both) and 21q22 (rs2836878[A]; P = 6.01 x 10(-8); OR = 0.73), located close to the TNFRSF6B and PSMG1 genes, respectively.