Genome-wide association study of height-adjusted BMI in childhood identifies functional variant in ADCY3

Objective: Genome-wide association studies (GWAS) of BMI are mostly undertaken under the assumption that "kg/m2" is an index of weight fully adjusted for height, but in general this is not true. The aim here was to assess the contribution of common genetic variation to a adjusted version of that phenotype which appropriately accounts for covariation in height in children. Methods: A GWAS of height-adjusted BMI (BMI[x]=weight/heightx), calculated to be uncorrelated with height, in 5809 participants (mean age 9.9 years) from the Avon Longitudinal Study of Parents and Children (ALSPAC) was performed. Results: GWAS based on BMI[x] yielded marked differences in genomewide results profile. SNPs in ADCY3 (adenylate cyclase 3) were associated at genome-wide significance level (rs11676272 (0.28 kg/m3.1 change per allele G (0.19, 0.38), P=6 × 10-9). In contrast, they showed marginal evidence of association with conventional BMI [rs11676272 (0.25 kg/m2 (0.15, 0.35), P=6 × 10-7)]. Results were replicated in an independent sample, the Generation R study. Conclusions: Analysis of BMI[x] showed differences to that of conventional BMI. The association signal at ADCY3 appeared to be driven by a missense variant and it was strongly correlated with expression of this gene. Our work highlights the importance of well understood phenotype use (and the danger of convention) in characterising genetic contributions to complex traits

A Comparison of Three Methods of Mendelian Randomization when the Genetic Instrument, the Risk Factor and the Outcome Are All Binary

The method of instrumental variable (referred to as Mendelian randomization when the instrument is a genetic variant) has been initially developed to infer on a causal effect of a risk factor on some outcome of interest in a linear model. Adapting this method to nonlinear models, however, is known to be problematic. In this paper, we consider the simple case when the genetic instrument, the risk factor, and the outcome are all binary. We compare via simulations the usual two-stages estimate of a causal odds-ratio and its adjusted version with a recently proposed estimate in the context of a clinical trial with noncompliance. In contrast to the former two, we confirm that the latter is (under some conditions) a valid estimate of a causal odds-ratio defined in the subpopulation of compliers, and we propose its use in the context of Mendelian randomization. By analogy with a clinical trial with noncompliance, compliers are those individuals for whom the presence/absence of the risk factor X is determined by the presence/absence of the genetic variant Z (i.e., for whom we would observe X = Z whatever the alleles randomly received at conception). We also recall and illustrate the huge variability of instrumental variable estimates when the instrument is weak (i.e., with a low percentage of compliers, as is typically the case with genetic instruments for which this proportion is frequently smaller than 10%) where the inter-quartile range of our simulated estimates was up to 18 times higher compared to a conventional (e.g., intention-to-treat) approach. We thus conclude that the need to find stronger instruments is probably as important as the need to develop a methodology allowing to consistently estimate a causal odds-ratio

Adult height variants affect birth length and growth rate in children

Previous studies identified 180 single nucleotide polymorphisms (SNPs) associated with adult height, explaining ∼10% of the variance. The age at which these begin to affect growth is unclear. We modelled the effect of these SNPs on birth length and childhood growth. A total of 7768 participants in the Avon Longitudinal Study of Parents and Children had data available. Individual growth trajectories from 0 to 10 years were estimated using mixed-effects linear spline models and differences in trajectories by individual SNPs and allelic score were determined. The allelic score was associated with birth length (0.026 cm increase per ‘tall’ allele, SE = 0.003, P = 1 × 10−15, equivalent to 0.017 SD). There was little evidence of association between the allelic score and early infancy growth (0–3 months), but there was evidence of association between the allelic score and later growth. This association became stronger with each consecutive growth period, per ‘tall’ allele per month effects were 0.015 SD (3 months–1 year, SE = 0.004), 0.023 SD (1–3 years, SE = 0.003) and 0.028 SD (3–10 years, SE = 0.003). By age 10, the mean height difference between individuals with ≤170 versus ≥191 ‘tall’ alleles (the top and bottom 10%) was 4.7 cm (0.8 SD), explaining ∼5% of the variance. There was evidence of associations with specific growth periods for some SNPs (rs3791675, EFEMP1 and rs6569648, L3MBTL3) and supportive evidence for previously reported age-dependent effects of HHIP and SOCS2 SNPs. SNPs associated with adult height influence birth length and have an increasing effect on growth from late infancy through to late childhood. By age 10, they explain half the height variance (∼5%) of that explained in adults (∼10%)

FTO Is Expressed in Neurones throughout the Brain and Its Expression Is Unaltered by Fasting

Single-nucleotide polymorphisms in the first intron of the ubiquitously expressed FTO gene are associated with obesity. Although the physiological functions of FTO remain unclear, food intake is often altered when Fto expression levels are manipulated. Furthermore, deletion of FTO from neurones alone has a similar effect on food intake to deletion of FTO in all tissues. These results indicate that FTO expression in the brain is particularly important. Considerable focus has been placed on the dynamic regulation of Fto mRNA expression in the hypothalamus after short-term (16–48 hour) fasting, but results have been controversial. There are no studies that quantify FTO protein levels across the brain, and assess its alteration following short-term fasting. Using immunohistochemistry, we found that FTO protein is widely expressed in mouse brain, and present in the majority of neurones. Using quantitative Western blotting and RT-qPCR we show that FTO protein and mRNA levels in the hypothalamus, cerebellum and rostral brain are relatively uniform, and levels in the brain are higher than in skeletal muscles of the lower limbs. Fasting for 18 hours does not alter the expression pattern, or levels, of FTO protein and mRNA. We further show that the majority of POMC neurones, which are critically involved in food intake regulation, also express FTO, but that the percentage of FTO-positive POMC neurones is not altered by fasting. In summary, we find no evidence that Fto/FTO expression is regulated by short-term (18-hour) fasting. Thus, it is unlikely that the hunger and increased post-fasting food intake caused by such food deprivation is driven by alterations in Fto/FTO expression. The widespread expression of FTO in neurones also suggests that physiological studies of this protein should not be limited to the hypothalamus

A meta-analysis of the associations between common variation in the PDE8B gene and thyroid hormone parameters, including assessment of longitudinal stability of associations over time and effect of thyroid hormone replacement

Objective Common variants in PDE8B are associated with TSH but apparently without any effect on thyroid hormone levels that is difficult to explain. Furthermore, the stability of the association has not been examined in longitudinal studies or in patients on levothyroxine (l-T4). Design Totally, four cohorts were used (n=2557): the Busselton Health Study (thyroid function measured on two occasions), DEPTH, EFSOCH (selective cohorts), and WATTS (individuals on l-T4). Methods Meta-analysis to clarify associations between the rs4704397 single nucleotide polymorphism in PDE8B on TSH, tri-iodothyronine (T3), and T4 levels. Results Meta-analysis confirmed that genetic variation in PDE8B was associated with TSH (P=1.64×10−10 0.20 s.d./allele, 95% confidence interval (CI) 0.142, 0.267) and identified a possible new association with free T4 (P=0.023, −0.07 s.d./allele, 95% CI −0.137, −0.01), no association was seen with free T3 (P=0.218). The association between PDE8B and TSH was similar in 1981 (0.14 s.d./allele, 95% CI 0.04, 0.238) and 1994 (0.20 s.d./allele, 95% CI 0.102, 0.300) and even more consistent between PDE8B and free T4 in 1981 (−0.068 s.d./allele, 95% CI −0.167, 0.031) and 1994 (−0.07 s.d./allele, 95% CI −0.170, 0.030). No associations were seen between PDE8B and thyroid hormone parameters in individuals on l-T4. Conclusion Common genetic variation in PDE8B is associated with reciprocal changes in TSH and free T4 levels that are consistent over time and lost in individuals on l-T4. These findings identify a possible genetic marker reflecting variation in thyroid hormone output that will be of value in epidemiological studies and provides additional evidence that PDE8B is involved in TSH signaling in the thyroid

Thyroid Function and Body Weight: A Community-Based Longitudinal Study

OBJECTIVE: Body weight and overt thyroid dysfunction are associated. Cross-sectional population-based studies have repeatedly found that thyroid hormone levels, even within the normal reference range, might be associated with body weight. However, for longitudinal data, the association is less clear. Thus, we tested the association between serum thyrotropin (TSH) and body weight in a community-based sample of adult persons followed for 11 years. METHODS: A random sample of 4,649 persons aged 18-65 years from a general population participated in the DanThyr study in 1997-8. We included 2,102 individuals who participated at 11-year follow-up, without current or former treatment for thyroid disease and with measurements of TSH and weight at both examinations. Multiple linear regression models were used, stratified by sex and adjusted for age, smoking status, and leisure time physical activity. RESULTS: Baseline TSH concentration was not associated with change in weight (women, P = 0.17; men, P = 0.72), and baseline body mass index (BMI) was not associated with change in TSH (women, P = 0.21; men, P = 0.85). Change in serum TSH and change in weight were significantly associated in both sexes. Weight increased by 0.3 kg (95% confidence interval [CI] 0.1, 0.4, P = 0.005) in women and 0.8 kg (95% CI 0.1, 1.4, P = 0.02) in men for every one unit TSH (mU/L) increase. CONCLUSIONS: TSH levels were not a determinant of future weight changes, and BMI was not a determinant for TSH changes, but an association between weight change and TSH change was present

Maturation in serum thyroid function parameters over childhood and puberty:results of a longitudinal study

Context: Serum thyroid hormone levels differ between children and adults, but have not been studied longitudinally through childhood. Objective: To assess changes in thyroid-stimulating hormone (TSH) and thyroid hormone levels over childhood and their interrelationships. Design: Cohort study. Setting: The Avon Longitudinal Study of Parents and Children, a population-based birth cohort. Participants: A total of 4442 children who had thyroid function measured at age 7, and 1263 children who had thyroid function measured at age 15. Eight hundred eighty-four children had measurements at both ages. Main Outcome Measures: Reference ranges for TSH, free tri-iodothyronine (FT3), free thyroxine (FT4), their longitudinal stability, and interrelationships. Results: Children at age 7 years had a higher FT3 [6.17 pmol/L, standard deviation (SD) 0.62] than children at age 15 (5.83 pmol/L, SD 0.74); P , 0.0001 with 23.2% of children at age 7 having FT3 above the adult reference range. Higher FT3 levels at age 7 in boys (P = 0.0001) and girls (P = 0.04) were associated with attainment of a more advanced pubertal stage at age 13. TSH was positively associated with FT3 at age 7 and age 15 even after adjusting for confounders. In contrast, TSH was negatively associated with FT

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

A genome-wide association meta-analysis of self-reported allergy identifies shared and allergy-specific susceptibility loci

Allergic disease is very common and carries substantial public-health burdens. We conducted a meta-analysis of genome-wide associations with self-reported cat, dust-mite and pollen allergies in 53,862 individuals. We used generalized estimating equations to model shared and allergy-specific genetic effects. We identified 16 shared susceptibility loci with association P < 5 × 10-8, including 8 loci previously associated with asthma, as well as 4p14 near TLR1, TLR6 and TLR10 (rs2101521, P = 5.3 × 10 -21); 6p21.33 near HLA-C and MICA (rs9266772, P = 3.2 × 10 -12); 5p13.1 near PTGER4 (rs7720838, P = 8.2 × 10 -11); 2q33.1 in PLCL1 (rs10497813, P = 6.1 × 10-10), 3q28 in LPP (rs9860547, P = 1.2 × 10-9); 20q13.2 in NFATC2 (rs6021270, P = 6.9 × 10-9), 4q27 in ADAD1 (rs17388568, P = 3.9 × 10-8); and 14q21.1 near FOXA1 and TTC6 (rs1998359, P = 4.8 × 10-8). We identified one locus with substantial evidence of differences in effects across allergies at 6p21.32 in the class II human leukocyte antigen (HLA) region (rs17533090, P = 1.7 × 10-12), which was strongly associated with cat allergy. Our study sheds new light on the shared etiology of immune and autoimmune disease