Large trans-ethnic meta-analysis identifies AKR1C4 as a novel gene associated with age at menarche.

STUDY QUESTION: Does the expansion of genome-wide association studies (GWAS) to a broader range of ancestries improve the ability to identify and generalise variants associated with age at menarche (AAM) in European populations to a wider range of world populations? SUMMARY ANSWER: By including women with diverse and predominantly non-European ancestry in a large-scale meta-analysis of AAM with half of the women being of African ancestry, we identified a new locus associated with AAM in African-ancestry participants, and generalised loci from GWAS of European ancestry individuals. WHAT IS KNOWN ALREADY: AAM is a highly polygenic puberty trait associated with various diseases later in life. Both AAM and diseases associated with puberty timing vary by race or ethnicity. The majority of GWAS of AAM have been performed in European ancestry women. STUDY DESIGN, SIZE, DURATION: We analysed a total of 38 546 women who did not have predominantly European ancestry backgrounds: 25 149 women from seven studies from the ReproGen Consortium and 13 397 women from the UK Biobank. In addition, we used an independent sample of 5148 African-ancestry women from the Southern Community Cohort Study (SCCS) for replication. PARTICIPANTS/MATERIALS, SETTING, METHODS: Each AAM GWAS was performed by study and ancestry or ethnic group using linear regression models adjusted for birth year and study-specific covariates. ReproGen and UK Biobank results were meta-analysed using an inverse variance-weighted average method. A trans-ethnic meta-analysis was also carried out to assess heterogeneity due to different ancestry. MAIN RESULTS AND THE ROLE OF CHANCE: We observed consistent direction and effect sizes between our meta-analysis and the largest GWAS conducted in European or Asian ancestry women. We validated four AAM loci (1p31, 6q16, 6q22 and 9q31) with common genetic variants at P < 5 × 10-7. We detected one new association (10p15) at P < 5 × 10-8 with a low-frequency genetic variant lying in AKR1C4, which was replicated in an independent sample. This gene belongs to a family of enzymes that regulate the metabolism of steroid hormones and have been implicated in the pathophysiology of uterine diseases. The genetic variant in the new locus is more frequent in African-ancestry participants, and has a very low frequency in Asian or European-ancestry individuals. LARGE SCALE DATA: N/A. LIMITATIONS, REASONS FOR CAUTION: Extreme AAM (<9 years or >18 years) were excluded from analysis. Women may not fully recall their AAM as most of the studies were conducted many years later. Further studies in women with diverse and predominantly non-European ancestry are needed to confirm and extend these findings, but the availability of such replication samples is limited. WIDER IMPLICATIONS OF THE FINDINGS: Expanding association studies to a broader range of ancestries or ethnicities may improve the identification of new genetic variants associated with complex diseases or traits and the generalisation of variants from European-ancestry studies to a wider range of world populations. STUDY FUNDING/COMPETING INTEREST(S): Funding was provided by CHARGE Consortium grant R01HL105756-07: Gene Discovery For CVD and Aging Phenotypes and by the NIH grant U24AG051129 awarded by the National Institute on Aging (NIA). The authors have no conflict of interest to declare.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index.

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10-5, Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10-06/Pfemales: 3.45 × 10-07/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

Dissecting heterogeneity in paediatric Type 1 diabetes: association of TCF7L2 rs7903146 TT and low-risk human leukocyte antigen (HLA) genotypes.

AIMS: To test the hypothesis that non-obese individuals with childhood-onset Type 1 diabetes and the rs7903146 TT genotype would be less likely to have high-risk human leukocyte antigen (HLA) genotypes and alleles. METHODS: We studied a cohort of 105 non-obese participants in the T1D Exchange Biobank Residual Insulin Study who had childhood-onset Type 1 diabetes [mean (sd) age at onset and recruitment, respectively, 9.9 (4.15) and 14.4 (4.13) years; 84.8% non-Hispanic white]. We analysed islet autoantibodies (glutamic acid decarboxylase 65, islet cell autoantigen 512/islet antigen-2 and zinc transporter 8), non-fasting random C-peptide levels, HLA type and TCF7L2 single nucleotide polymorphism rs7903146 in this cohort. RESULTS: None of the 13 individuals with the rs7903146 TT genotype carried the highest Type 1 diabetes risk HLA genotype, i.e. DRB1*03:01/DR4 (DRB1*0401, *04:05 or *04:02), compared with 29.4% (27/92) of those without it (P=0.023). The DRB1*03:01 allele was present in 15.4% (2/13) of individuals with the single nucleotide polymorphism, compared with 59.8% (55/92) of those without it (P=0.003). Analyses restricted to autoantibody-positive individuals (n=80) yielded similar results. The HLA DRB1*15:01 allele, which affords dominant protection against Type 1 diabetes, was found in one participant, who had multiple islet autoantibodies and carried the rs7903146 TT genotype. CONCLUSIONS: These findings further support the hypothesis that TCF7L2 gene variation contributes to diabetogenesis in a subset of young people with Type 1 diabetes, opening possible new pathways for therapy and prevention.

Targeted resequencing identifies defective variants of decoy receptor 3 in pediatric-onset inflammatory bowel disease.

Genome-wide association studies have implicated common variation at the 20q13 locus in inflammatory bowel disease, particularly for the pediatric Crohn's form. This locus harbors tumor necrosis factor receptor superfamily (TNFRSF6B), encoding a secreted protein, decoy receptor 3 (DcR3), which binds to and neutralizes pro-inflammatory cytokines of the tumor necrosis factor superfamily. We sought to further the evidence of DcR3's role in pediatric IBD by identifying missense mutations with functional significance within TNFRSF6B. We sequenced the exons of the gene in 528 Caucasian pediatric IBD cases and 549 Caucasian healthy controls to establish the frequency of such events in each population. Sequencing revealed that our IBD cohort harbored a greater number of missense variants, yielding an odds ratio of 3.9 (P-value=0.005). Using functional assays, we established that the frequency of mutants defective in secretion from cultured cells was greater in the Crohn's category than in the controls, yielding an odds ratio of 7.1 (P-value=0.004). These results suggest that rare defective variants in TNFRSF6B have a role in the pathogenesis of some cases of IBD and that interventions targeting this group of tumor necrosis factor-family members may benefit patients with IBD.

Reduced bone density and osteoporosis associated with a polymorphic Sp1 binding site in the collagen type I alpha 1 gene.

Osteoporosis is a common disease with a strong genetic component, characterized by reduced bone mass and increased fracture risk. Current evidence suggests that the inheritance of bone mass is under polygenic control but the genes responsible are poorly defined. Type I collagen is the major protein of bone encoded by the COLIA1 and COLIA2 genes. While these are strong candidates for the genetic regulation of bone mass, no abnormality of either gene has so far been defined in osteoporosis. In this study, we describe a novel G-->T polymorphism in a regulatory region of COLIA1 at a recognition site for the transcription factor Sp1(7) that is significantly related to bone mass and osteoporotic fracture. G/T heterozygotes at the polymorphic Sp1 site (Ss) had significantly lower bone mineral density (BMD) than G/G homozygotes (SS) in two populations of British women and BMD was lower still in T/T homozygotes (ss). The unfavourable Ss and ss genotypes were over-represented in patients with severe osteoporosis and vertebral fractures (54%), as compared with controls (27%), equivalent to a relative risk of 2.97 (95% confidence interval 1.63-9.56) for vertebral fracture in individuals who carry the 's' allele. While the mechanisms that underlie this association remain to be defined, the COLIA1 Sp1 polymorphism appears to be an important marker for low bone mass and vertebral fracture, raising the possibility that genotyping at this site may be of value in identifying women who are at risk of osteoporosis.

Perturbation of the insomnia WDR90 GWAS locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q.

Although genome wide association studies (GWAS) have been crucial for the identification of loci associated with sleep traits and disorders, the method itself does not directly uncover the underlying causal variants and corresponding effector genes. The overwhelming majority of such variants reside in non-coding regions and are therefore presumed to impact the activity of cis-regulatory elements, such as enhancers. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated PIG-Q as a functionally relevant gene at the insomnia 'WDR90' locus. However, importantly that effort did not characterize the corresponding underlying causal variant at this GWAS signal. Specifically, our genome-wide ATAC-seq and high-resolution promoter-focused Capture C datasets generated in this cell setting brought our attention to a shortlist of three tightly neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium in a candidate intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. The objective of this study was to investigate the influence of the proxy SNPs collectively and then individually on PIG-Q modulation and to pinpoint the causal "regulatory" variant among the three SNPs. Starting at a gross level perturbation, deletion of the entire region harboring all three SNPs in human iPSC-derived neural progenitor cells via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from more refined individual luciferase reporter assays for each of the three SNPs in iPSCs revealed that the intronic region with the rs3752495 risk allele induced a ~2.5-fold increase in luciferase expression (n=10). Importantly, rs3752495 also exhibited an allele specific effect, with the risk allele increasing the luciferase expression by ~2-fold compared to the non-risk allele. In conclusion, our variant-to-function approach and subsequent in vitro validation implicates rs3752495 as a causal insomnia risk variant embedded at the WDR90-PIG-Q locus.

Common genetic variants differentially influence the transition from clinically defined states of fasting glucose metabolism.

AIMS/HYPOTHESIS: Common genetic variants have been associated with type 2 diabetes. We hypothesised that a subset of these variants may have different effects on the transition from normal fasting glucose (NFG) to impaired fasting glucose (IFG) than on that from IFG to diabetes. METHODS: We identified 16 type 2 diabetes risk variants from the Illumina Broad Candidate-gene Association Resource (CARe) array genotyped in 26,576 CARe participants. Participants were categorised at baseline as NFG, IFG or type 2 diabetic (n = 16,465, 8,017 or 2,291, respectively). Using Cox proportional hazards and likelihood ratio tests (LRTs), we compared rates of progression by genotype for 4,909 (NFG to IFG) and 1,518 (IFG to type 2 diabetes) individuals, respectively. We then performed multinomial regression analyses at baseline, comparing the risk of assignment to the NFG, IFG or diabetes groups by genotype. RESULTS: The rate of progression from NFG to IFG was significantly greater in participants carrying the risk allele at MTNR1B (p = 1 × 10(-4)), nominally greater at GCK and SLC30A8 (p < 0.05) and nominally smaller at IGF2BP2 (p = 0.01) than the rate of progression from IFG to diabetes by the LRT. Results of the baseline, multinomial regression model were consistent with these findings. CONCLUSIONS/INTERPRETATION: Common genetic risk variants at GCK, SLC30A8, IGF2BP2 and MTNR1B influence to different extents the development of IFG and the transition from IFG to type 2 diabetes. Our findings may have implications for understanding the genetic contribution of these variants to the development of IFG and type 2 diabetes.

Transferability and fine-mapping of glucose and insulin quantitative trait loci across populations: CARe, the Candidate Gene Association Resource.

AIMS/HYPOTHESIS: Hyperglycaemia disproportionately affects African-Americans (AfAs). We tested the transferability of 18 single-nucleotide polymorphisms (SNPs) associated with glycaemic traits identified in European ancestry (EuA) populations in 5,984 non-diabetic AfAs. METHODS: We meta-analysed SNP associations with fasting glucose (FG) or insulin (FI) in AfAs from five cohorts in the Candidate Gene Association Resource. We: (1) calculated allele frequency differences, variations in linkage disequilibrium (LD), fixation indices (F(st)s) and integrated haplotype scores (iHSs); (2) tested EuA SNPs in AfAs; and (3) interrogated within ± 250 kb around each EuA SNP in AfAs. RESULTS: Allele frequency differences ranged from 0.6% to 54%. F(st) exceeded 0.15 at 6/16 loci, indicating modest population differentiation. All iHSs were <2, suggesting no recent positive selection. For 18 SNPs, all directions of effect were the same and 95% CIs of association overlapped when comparing EuA with AfA. For 17 of 18 loci, at least one SNP was nominally associated with FG in AfAs. Four loci were significantly associated with FG (GCK, p = 5.8 × 10(-8); MTNR1B, p = 8.5 × 10(-9); and FADS1, p = 2.2 × 10(-4)) or FI (GCKR, p = 5.9 × 10(-4)). At GCK and MTNR1B the EuA and AfA SNPs represented the same signal, while at FADS1, and GCKR, the EuA and best AfA SNPs were weakly correlated (r(2) <0.2), suggesting allelic heterogeneity for association with FG at these loci. CONCLUSIONS/INTERPRETATION: Few glycaemic SNPs showed strict evidence of transferability from EuA to AfAs. Four loci were significantly associated in both AfAs and those with EuA after accounting for varying LD across ancestral groups, with new signals emerging to aid fine-mapping.

Lack of relationship between mitochondrial heteroplasmy or variation and childhood obesity.

Mitochondrial electron transport has a central role in regulating energy supply within a cell. We hypothesized that mitochondrial variants or increased levels of mitochondrial heteroplasmy could be associated with common childhood obesity through their effects on mitochondrial function. To investigate this question, we queried two genome-wide genotyped childhood obesity datasets, consisting of 1080 European-American (EA) obese children (defined as BMI ≥ 95th percentile) together with 2500 EA lean controls (defined as BMI<50th percentile) and 1479 African-American (AA) obese children and 1575 AA lean controls. Association was not observed between childhood obesity and any of the assayed mitochondrial polymorphisms in either ethnicity. We also found no observable differences in heteroplasmy between each obese and non-obese group. Finally, we analyzed the quantitative mitochondrial genotype cells generated, whether they exceeded the heteroplasmy threshold or not. With this more lenient test, we found six positions with a significant difference between EA cases and controls (P<1 × 10(-4)). However, when evaluating the AA data set, no differences were noted at these sites, suggesting that our initial observations were because of chance rather than a meaningful relationship to childhood obesity. As such, it is unlikely that common mitochondrial polymorphisms or heteroplasmy have a role in childhood obesity.

Rare structural variants found in attention-deficit hyperactivity disorder are preferentially associated with neurodevelopmental genes.

Attention-deficit/hyperactivity disorder (ADHD) is a common and highly heritable disorder, but specific genetic factors underlying risk remain elusive. To assess the role of structural variation in ADHD, we identified 222 inherited copy number variations (CNVs) within 335 ADHD patients and their parents that were not detected in 2026 unrelated healthy individuals. Although no excess CNVs, either deletions or duplications, were found in the ADHD cohort relative to controls, the inherited rare CNV-associated gene set was significantly enriched for genes reported as candidates in studies of autism, schizophrenia and Tourette syndrome, including A2BP1, AUTS2, CNTNAP2 and IMMP2L. The ADHD CNV gene set was also significantly enriched for genes known to be important for psychological and neurological functions, including learning, behavior, synaptic transmission and central nervous system development. Four independent deletions were located within the protein tyrosine phosphatase gene, PTPRD, recently implicated as a candidate gene for restless legs syndrome, which frequently presents with ADHD. A deletion within the glutamate receptor gene, GRM5, was found in an affected parent and all three affected offspring whose ADHD phenotypes closely resembled those of the GRM5 null mouse. Together, these results suggest that rare inherited structural variations play an important role in ADHD development and indicate a set of putative candidate genes for further study in the etiology of ADHD.

Disease-associated loci are significantly over-represented among genes bound by transcription factor 7-like 2 (TCF7L2) in vivo.

AIMS/HYPOTHESIS: Transcription factor 7-like 2 (TCF7L2) has been strongly implicated in type 2 diabetes and cancer. Our goal was to identify the DNA sequences bound by this transcription factor in vivo. METHODS: We applied chromatin immunoprecipitation and sequencing to globally identify and map human DNA sequences bound by TCF7L2 in the colorectal carcinoma cell line, HCT116, where it is abundantly expressed. RESULTS: We identified 1,095 discrete binding sites across the genome, of which a subset were within 5 kb of 548 annotated NCBI Reference Sequence (RefSeq) genes. Despite using a cancer cell line, the most significant functions represented using pathway analysis software were related to diabetes, genetic disorders and coronary artery disease. As one of the enriched categories was related to genetic disorders, we queried our results against all published genome-wide association studies (GWAS) and found a highly significant over-representation of reported loci from among the genes bound by TCF7L2 within 5 kb (p = 7.50 × 10⁻¹5). This observation was primarily driven by excess loci revealed from GWAS of metabolic and cardiovascular traits; however, there was no or only minor enrichment of GWAS-derived loci for cancer, and inflammatory or neurological diseases. Of the specific traits, the most enriched loci were for type 2 diabetes and height. When defining the distance from genes at 50 kb or 500 kb, this enrichment pattern persisted, with some additional evidence for enrichment of cancer-related loci. CONCLUSIONS/INTERPRETATION: A highly significant proportion of genes bound by TCF7L2 are known disease-associated loci. These findings suggest that TCF7L2 is a central node in the regulation of human diabetes and other disease-associated genes.

Association of RASGRP1 with type 1 diabetes is revealed by combined follow-up of two genome-wide studies.

BACKGROUND: The two genome-wide association studies published by us and by the Wellcome Trust Case-Control Consortium (WTCCC) revealed a number of novel loci, but neither had the statistical power to elucidate all of the genetic components of type 1 diabetes risk, a task for which larger effective sample sizes are needed. METHODS: We analysed data from two sources: (1) The previously published second stage of our study, with a total sample size of the two stages consisting of 1046 Canadian case-parent trios and 538 multiplex families with 929 affected offspring from the Type 1 Diabetes Genetics Consortium (T1DGC); (2) the Rapid Response 2 (RR2) project of the T1DGC, which genotyped 4417 individuals from 1062 non-overlapping families, including 2059 affected individuals (mostly sibling pairs) for the 1536 markers with the highest statistical significance for type 1 diabetes in the WTCCC results. RESULTS: One locus, mapping to a linkage disequilibrium (LD) block at chr15q14, reached statistical significance by combining results from two markers (rs17574546 and rs7171171) in perfect LD with each other (r2 = 1). We obtained a joint p value of 1.3 x 10(-6), which exceeds by an order of magnitude the conservative threshold of 3.26 x 10(-5) obtained by correcting for the 1536 single nucleotide polymorphisms (SNPs) tested in our study. Meta-analysis with the original WTCCC genome-wide data produced a p value of 5.83 x 10(-9). CONCLUSIONS: A novel type 1 diabetes locus was discovered. It involves RASGRP1, a gene known to play a crucial role in thymocyte differentiation and T cell receptor (TCR) signalling by activating the Ras signalling pathway.

Remapping the type I diabetes association of the CTLA4 locus.

The Type I Diabetes Genetics Consortium genotyped 24 single-nucleotide polymorphisms (SNPs) in the CTLA4 locus in 2298 type I diabetes (T1D) nuclear families (11 159 individuals, 5003 affected) to evaluate the recognized T1D association. The 24 CTLA4 SNPs span approximately 43 kb from the 5' flanking to 3' flanking region of the gene in the middle of an extended region of linkage disequilibrium of more than 100 kb. The genotyping was performed using two technologies (Illumina GoldenGate and Sequenom iPlex) on the same samples. The genotype calls by both the methods were highly consistent (the majority >99%). Previously reported T1D association from both the +49G>A and the CT60 SNPs was replicated. The reported association of the -319C>T SNP was not replicated. Although associated with T1D risk, it is likely that neither SNP is causative, as the peak of T1D association was from the SNP rs231727 at 3' flanking of the CTLA4 gene. Comprehensive resequencing and fine mapping of the CTLA4 region are still needed to clarify the causal variants.

The type I diabetes association of the IL2RA locus.

To confirm and fine map previous reports of association, the Type I Diabetes (T1D) Genetics Consortium (T1DGC) assembled a large collection of DNA samples from affected sib-pair (ASP) families with T1D (5003 affected individuals) and genotyped polymorphic markers. One of these loci, involving the IL2RA gene, had been reported to be due to three independent effects. The T1DGC genotyped 69 single-nucleotide polymorphisms (SNPs) that span approximately 88 kb from the 5' flanking to 3' flanking region of the IL2RA locus. The most highly associated SNP reported earlier (ss52580101) was not included in the genotyping list; however, a 5-SNP (rs3134883, rs3118470, rs7072793, rs4749955 and rs12251307) haplotype (H5) was identified that strongly tagged its minor allele with r(2)=0.869 (95% CI, 0.850-0.885). This haplotype was significantly protective (P=3.2 x 10(-5)) in the T1D ASP families, with an odds ratio virtually identical to that reported for ss52580101. The SNP marking the second independent locus, (rs11594656) showed no association in the T1DGC set and the third (rs2104286) could not be distinguished, by conditional regression, from H5. Instead, the most significant independent effect was detected from the 5' flanking IL2RA SNP rs4749955, which remained significant after regression for H5. Thus, we confirm independent effects at the IL2RA locus.

A COL1A1 Sp1 binding site polymorphism predisposes to osteoporotic fracture by affecting bone density and quality.

Osteoporosis is a common disease with a strong genetic component. We previously described a polymorphic Sp1 binding site in the COL1A1 gene that has been associated with osteoporosis in several populations. Here we explore the molecular mechanisms underlying this association. A meta-analysis showed significant associations between COL1A1 "s" alleles and bone mineral density (BMD), body mass index (BMI), and osteoporotic fractures. The association with fracture was stronger than expected on the basis of the observed differences in BMD and BMI, suggesting an additional effect on bone strength. Gel shift assays showed increased binding affinity of the "s" allele for Sp1 protein, and primary RNA transcripts derived from the "s" allele were approximately three times more abundant than "S" allele--derived transcripts in "Ss" heterozygotes. Collagen produced from osteoblasts cultured from "Ss" heterozygotes had an increased ratio of alpha 1(I) protein relative to alpha 2(I), and this was accompanied by an increased ratio of COL1A1 mRNA relative to COL1A2. Finally, the yield strength of bone derived from "Ss" individuals was reduced when compared with bone derived from "SS" subjects. We conclude that the COL1A1 Sp1 polymorphism is a functional genetic variant that predisposes to osteoporosis by complex mechanisms involving changes in bone mass and bone quality.

FGF21 deficiency is associated with childhood obesity, insulin resistance and hypoadiponectinaemia: The BCAMS Study.

OBJECTIVE: Fibroblast growth factor 21 (FGF21) exerts beneficial effects on metabolic homoeostasis and has been reported to be regulated by adiponectin, leptin and resistin. However, while an association between increased circulating FGF21 and metabolic disorders has been reported in adults, paediatric-specific data are lacking. DESIGN AND METHODS: This study investigated the relationship between FGF21 levels and obesity, insulin resistance (IR), the metabolic syndrome (MetS) and adipokines (adiponectin, leptin and resistin) in a cohort of 3231 Chinese youngsters aged 6-18. RESULTS: There were gender- and puberty-related differences in FGF21 levels. Unexpectedly, FGF21 levels were decreased in children with obesity, and negatively correlated with insulin, HOMA-IR and leptin levels after adjusting for age, gender, puberty and lifestyle factors. Moreover, multiple regression analyses showed that serum FGF21 positively predicted adiponectin levels while resistin positively predicted FGF21 levels independent of BMI (P<0.05). Children in the lowest FGF21 quintile were more likely to have IR (OR: 1.85, 95% CI: 1.41-2.42; P=0.002) and MetS (OR: 1.62, 95% CI: 1.14-2.28; P=0.007) than those in the highest quintile. Further adjusting for BMI and/or the three adipokines modified the association of FGF21 with MetS (P>0.10) but not with IR (P<0.01). CONCLUSION: Although the associations between adiponectin, leptin, resistin and metabolic abnormalities in our paediatric population were similar to those in adults, correlations of FGF21 levels with obesity, IR and MetS were the inverse of those found in adults. Our present findings suggest that FGF21 deficiency, rather than resistance, contribute to IR and hypoadiponectinaemia independently of obesity in young people.

Influence of SNP*SNP interaction on BMI in European American adolescents: findings from the National Longitudinal Study of Adolescent Health.

BACKGROUND: Adolescent obesity is predictive of future weight gain, obesity and adult onset severe obesity (body mass index [BMI] ≥40 kg m(-2) ). Despite successful efforts to identify Single Nucleotide Polymorphisms (SNPs) influencing BMI, <5% of the 40-80% heritability of the phenotype has been explained. Identification of gene-gene (G-G) interactions between known variants can help explain this hidden heritability as well as identify potential biological mechanisms affecting weight gain during this critical developmental period. OBJECTIVE: We have recently shown distinct genetic effects on BMI across the life course, and thus it is important to examine the evidence for epistasis in adolescence. METHODS: In adolescent participants of European descent from wave II of the National Longitudinal Study of Adolescent Health (Add Health, n = 5072, ages 12-21, 52.5% female), we tested 34 established BMI-related SNPs for G-G interaction effects on BMI z-score. We used mixed-effects regression, assuming multiplicative interaction models adjusting for age, sex and geographic region, with random effects for family and school. RESULTS: For 28 G-G interactions that were nominally significant (P < 0.05), we attempted to replicate our results in an adolescent sample from the Childhood European American Cohort from Philadelphia. In the replication study, one interaction (PRKD1-FTO) was significant after correction for multiple testing. CONCLUSIONS: Our results are suggestive of epistatic effects on BMI during adolescence and point to potentially interactive effects between genes in biological pathways important in obesity.

Genes and osteoporosis.

Osteoporosis is a common disease associated with reduced bone mineral density, affecting up to 40% of women and 12% of men at some point during life. Although osteoporosis is multifactorial, genetic factors play an important role in the pathogenesis of osteoporosis, as up to 75% of the variance in bone mass-a major risk factor for osteoporotic fracture-is genetically determined. Segregation analysis of bone mass in families has suggested a model whereby bone mass is under the control of a large number of genes with small effect, rather than a few genes with large effect. Although the molecular-genetic basis by which bone mass is regulated is incompletely defined, polymorphisms of several candidate genes have been linked to bone mass in clinical studies. Polymorphisms in the gene encoding the vitamin D receptor have been extensively investigated with differing results, and interleukin-6, transforming growth factor beta, and the estrogen receptor have been associated with bone mass in isolated studies. Recent work has identified a polymorphism in the binding site for the transcription factor Sp1 in the collagen type I alpha 1 gene. This polymorphism is associated with bone mass and osteoporotic fracture in three distinct populations, suggesting that the polymorphism may act as a marker for low bone density and fracture risk. Further studies will uncover many new genes and candidate loci that relate to bone mass. This information may be of value in defining new therapeutic targets in the prevention and management of osteoporosis and in developing genetic tests to assess osteoporotic fracture risk.

An Sp1 binding site polymorphism in the COLIA1 gene predicts osteoporotic fractures in both men and women.

Genetic factors play an important role in the pathogenesis of osteoporosis, and recent studies have shown that a polymorphic Sp1 binding site in collagen type I alpha1 (COLIA1) gene is associated with bone mass and vertebral fractures in women from the U.K. Information on the predictive value of the COLIA1 Sp1 polymorphism in other populations is limited, however, and no studies have yet been performed in osteoporotic males. In view of this, we analyzed COLIA1 genotypes in relation to bone density and biochemical markers of bone turnover and the presence of osteoporotic fractures in a case-control study of Danish men and women. COLIA1 genotype was determined by polymerase chain reaction analysis of genomic DNA extracted from peripheral blood samples and related to bone mass, biochemical markers of bone turnover, and the presence of fracture in a study of 375 osteoporotic vertebral fracture patients and normal controls. There was no significant effect of COLIA1 genotype on bone mass or biochemical markers when data from the control group (n = 195) and fracture group (n = 180) were analyzed separately. However, the genotype distribution was significantly different in the fracture cases compared with age-matched controls (chi2 = 16.48, n = 249,p = 0.0003) due mainly to over-representation of the ss genotype in the fracture patients (14.3% vs. 1.4%), equivalent to an odds ratio for vertebral fracture of 11.83 (95% confidence interval 2.64-52.97) in those with the ss genotype. Similar differences in genotype distribution between osteoporotic patients and controls were observed in both men (chi2 = 11.52, n = 95, p = 0.0032, OR = 2.04) and women (chi2 = 6.90, n = 154, p = 0.032, OR = 1.37). In keeping with the above, logistic regression analysis showed that the ss genotype was an independent predictor of osteoporotic fracture (p = 0.028). This study confirms that the COLIA1 Sp1 polymorphism is significantly associated with osteoporotic vertebral fractures. The association is seen in both men and women, and the effect on fracture risk appears to be partly independent of bone mineral density. Our results raise the possibility that genotyping at the Sp1 site could be of clinical value in identifying individuals at risk of osteoporotic fractures in both genders.