A mendelian randomization study based on data from multiple cohorts conducted by Karani Santhanakrishnan Vimaleswaran and colleagues re-examines the causal nature of the relationship between vitamin D levels and obesity.

Background

Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis.

Methods and Findings

We used information from 21 adult cohorts (up to 42,024 participants) with 12 BMI-related SNPs (combined in an allelic score) to produce an instrument for BMI and four SNPs associated with 25(OH)D (combined in two allelic scores, separately for genes encoding its synthesis or metabolism) as an instrument for vitamin D. Regression estimates for the IVs (allele scores) were generated within-study and pooled by meta-analysis to generate summary effects.

Associations between vitamin D scores and BMI were confirmed in the Genetic Investigation of Anthropometric Traits (GIANT) consortium (n = 123,864). Each 1 kg/m2 higher BMI was associated with 1.15% lower 25(OH)D (p = 6.52×10−27). The BMI allele score was associated both with BMI (p = 6.30×10−62) and 25(OH)D (−0.06% [95% CI −0.10 to −0.02], p = 0.004) in the cohorts that underwent meta-analysis. The two vitamin D allele scores were strongly associated with 25(OH)D (p≤8.07×10−57 for both scores) but not with BMI (synthesis score, p = 0.88; metabolism score, p = 0.08) in the meta-analysis. A 10% higher genetically instrumented BMI was associated with 4.2% lower 25(OH)D concentrations (IV ratio: −4.2 [95% CI −7.1 to −1.3], p = 0.005). No association was seen for genetically instrumented 25(OH)D with BMI, a finding that was confirmed using data from the GIANT consortium (p≥0.57 for both vitamin D scores).

Conclusions

On the basis of a bi-directional genetic approach that limits confounding, our study suggests that a higher BMI leads to lower 25(OH)D, while any effects of lower 25(OH)D increasing BMI are likely to be small. Population level interventions to reduce BMI are expected to decrease the prevalence of vitamin D deficiency.

Please see later in the article for the Editors' Summary

Editors' Summary

Background

Obesity—having an unhealthy amount of body fat—is increasing worldwide. In the US, for example, a third of the adult population is now obese. Obesity is defined as having a body mass index (BMI, an indicator of body fat calculated by dividing a person's weight in kilograms by their height in meters squared) of more than 30.0 kg/m2. Although there is a genetic contribution to obesity, people generally become obese by consuming food and drink that contains more energy than they need for their daily activities. Thus, obesity can be prevented by having a healthy diet and exercising regularly. Compared to people with a healthy weight, obese individuals have an increased risk of developing diabetes, heart disease and stroke, and tend to die younger. They also have a higher risk of vitamin D deficiency, another increasingly common public health concern. Vitamin D, which is essential for healthy bones as well as other functions, is made in the skin after exposure to sunlight but can also be obtained through the diet and through supplements.

Why Was This Study Done?

Observational studies cannot prove that obesity causes vitamin D deficiency because obese individuals may share other characteristics that reduce their circulating 25-hydroxy vitamin D [25(OH)D] levels (referred to as confounding). Moreover, observational studies cannot indicate whether the larger vitamin D storage capacity of obese individuals (vitamin D is stored in fatty tissues) lowers their 25(OH)D levels or whether 25(OH)D levels influence fat accumulation (reverse causation). If obesity causes vitamin D deficiency, monitoring and treating vitamin D deficiency might alleviate some of the adverse health effects of obesity. Conversely, if low vitamin D levels cause obesity, encouraging people to take vitamin D supplements might help to control the obesity epidemic. Here, the researchers use bi-directional “Mendelian randomization” to examine the direction and causality of the relationship between BMI and 25(OH)D. In Mendelian randomization, causality is inferred from associations between genetic variants that mimic the influence of a modifiable environmental exposure and the outcome of interest. Because gene variants do not change over time and are inherited randomly, they are not prone to confounding and are free from reverse causation. Thus, if a lower vitamin D status leads to obesity, genetic variants associated with lower 25(OH)D concentrations should be associated with higher BMI, and if obesity leads to a lower vitamin D status, then genetic variants associated with higher BMI should be associated with lower 25(OH)D concentrations.

What Did the Researchers Do and Find?

The researchers created a “BMI allele score” based on 12 BMI-related gene variants and two “25(OH)D allele scores,” which are based on gene variants that affect either 25(OH)D synthesis or breakdown. Using information on up to 42,024 participants from 21 studies, the researchers showed that the BMI allele score was associated with both BMI and with 25(OH)D levels among the study participants. Based on this information, they calculated that each 10% increase in BMI will lead to a 4.2% decrease in 25(OH)D concentrations. By contrast, although both 25(OH)D allele scores were strongly associated with 25(OH)D levels, neither score was associated with BMI. This lack of an association between 25(OH)D allele scores and obesity was confirmed using data from more than 100,000 individuals involved in 46 studies that has been collected by the GIANT (Genetic Investigation of Anthropometric Traits) consortium.

What Do These Findings Mean?

These findings suggest that a higher BMI leads to a lower vitamin D status whereas any effects of low vitamin D status on BMI are likely to be small. That is, these findings provide evidence for obesity as a causal factor in the development of vitamin D deficiency but not for vitamin D deficiency as a causal factor in the development of obesity. These findings suggest that population-level interventions to reduce obesity should lead to a reduction in the prevalence of vitamin D deficiency and highlight the importance of monitoring and treating vitamin D deficiency as a means of alleviating the adverse influences of obesity on health.

Additional Information

Please access these Web sites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.1001383.

The US Centers for Disease Control and Prevention provides information on all aspects of overweight and obesity (in English and Spanish); a data brief provides information about the vitamin D status of the US population

The World Health Organization provides information on obesity (in several languages)

The UK National Health Service Choices website provides detailed information about obesity and a link to a personal story about losing weight; it also provides information about vitamin D

The International Obesity Taskforce provides information about the global obesity epidemic

The US Department of Agriculture's ChooseMyPlate.gov website provides a personal healthy eating plan; the Weight-control Information Network is an information service provided for the general public and health professionals by the US National Institute of Diabetes and Digestive and Kidney Diseases (in English and Spanish)

The US Office of Dietary Supplements provides information about vitamin D (in English and Spanish)

MedlinePlus has links to further information about obesity and about vitamin D (in English and Spanish)

Wikipedia has a page on Mendelian randomization (note: Wikipedia is a free online encyclopedia that anyone can edit; available in several languages)

Overview and details of the collaborative large-scale genetic association study (D-CarDia) provide information about vitamin D and the risk of cardiovascular disease, diabetes and related traits

To identify previously unknown genetic loci associated with fasting glucose concentrations, we examined the leading association signals in ten genome-wide association scans involving a total of 36,610 individuals of European descent. Variants in the gene encoding melatonin receptor 1B (MTNR1B) were consistently associated with fasting glucose across all ten studies. The strongest signal was observed at rs10830963, where each G allele (frequency 0.30 in HapMap CEU) was associated with an increase of 0.07 (95% CI = 0.06-0.08) mmol/l in fasting glucose levels (P = 3.2 = × 10−50) and reduced beta-cell function as measured by homeostasis model assessment (HOMA-B, P = 1.1 × 10−15). The same allele was associated with an increased risk of type 2 diabetes (odds ratio = 1.09 (1.05-1.12), per G allele P = 3.3 × 10−7) in a meta-analysis of 13 case-control studies totaling 18,236 cases and 64,453 controls. Our analyses also confirm previous associations of fasting glucose with variants at the G6PC2 (rs560887, P = 1.1 × 10−57) and GCK (rs4607517, P = 1.0 × 10−25) loci.

Objectives

Prevention of obesity should start as early as possible after birth. We aimed to build clinically useful equations estimating the risk of later obesity in newborns, as a first step towards focused early prevention against the global obesity epidemic.

Methods

We analyzed the lifetime Northern Finland Birth Cohort 1986 (NFBC1986) (N = 4,032) to draw predictive equations for childhood and adolescent obesity from traditional risk factors (parental BMI, birth weight, maternal gestational weight gain, behaviour and social indicators), and a genetic score built from 39 BMI/obesity-associated polymorphisms. We performed validation analyses in a retrospective cohort of 1,503 Italian children and in a prospective cohort of 1,032 U.S. children.

Results

In the NFBC1986, the cumulative accuracy of traditional risk factors predicting childhood obesity, adolescent obesity, and childhood obesity persistent into adolescence was good: AUROC = 0·78[0·74–0.82], 0·75[0·71–0·79] and 0·85[0·80–0·90] respectively (all p<0·001). Adding the genetic score produced discrimination improvements ≤1%. The NFBC1986 equation for childhood obesity remained acceptably accurate when applied to the Italian and the U.S. cohort (AUROC = 0·70[0·63–0·77] and 0·73[0·67–0·80] respectively) and the two additional equations for childhood obesity newly drawn from the Italian and the U.S. datasets showed good accuracy in respective cohorts (AUROC = 0·74[0·69–0·79] and 0·79[0·73–0·84]) (all p<0·001). The three equations for childhood obesity were converted into simple Excel risk calculators for potential clinical use.

Conclusion

This study provides the first example of handy tools for predicting childhood obesity in newborns by means of easily recorded information, while it shows that currently known genetic variants have very little usefulness for such prediction.

Background

The CHRNA5-CHRNA3-CHRNB4 gene cluster on 15q25 has consistently been associated with smoking quantity, nicotine dependence and lung cancer. Recent research also points towards its involvement in cardiovascular homeostasis, but studies in large human samples are lacking, especially on the role of the gene cluster in blood pressure regulation.

Methodology/Principal Findings

We studied the associations between 18 single nucleotide polymorphisms (SNPs) in CHRNA5-CHRNA3-CHRNB4 and systolic blood pressure (SBP), diastolic blood pressure (DBP), and body mass index (BMI) in 5402 young adults from the Northern Finland Birth Cohort 1966. We observed some evidence for associations between two SNPs and SBP and between six SNPs and BMI; the evidence for associations with DBP was weaker. The associations with the three phenotypes were driven by different loci with low linkage disequilibrium with each other. The associations appeared more pronounced in smokers, such that the smoking-increasing alleles would predict lower SBP and BMI. Each additional copy of the rs1948 G-allele and the rs950776 A-allele reduced SBP on average by −1.21 (95% CI −2.01, −0.40) mmHg in smokers. The variants associated with BMI included rs2036534, rs6495309, rs1996371, rs6495314, rs4887077 and rs11638372 and had an average effect size of −0.38 (−0.68, −0.08) kg/m2 per an additional copy of the risk allele in smokers. Formal assessments of interactions provided weaker support for these findings, especially after adjustment for multiple testing.

Conclusions

Variation at 15q25 appears to interact with smoking status in influencing SBP and BMI. The genetic loci associated with SBP were in low linkage disequilibrium with those associated with BMI suggesting that the gene cluster might regulate SBP through biological mechanisms that partly differ from those regulating BMI. Further studies in larger samples are needed for more precise evaluation of the possible interactions, and to understand the mechanisms behind.

Maternal smoking during pregnancy is associated with low birth weight. Common variation at rs1051730 is robustly associated with smoking quantity and was recently shown to influence smoking cessation during pregnancy, but its influence on birth weight is not clear. We aimed to investigate the association between this variant and birth weight of term, singleton offspring in a well-powered meta-analysis. We stratified 26 241 European origin study participants by smoking status (women who smoked during pregnancy versus women who did not smoke during pregnancy) and, in each stratum, analysed the association between maternal rs1051730 genotype and offspring birth weight. There was evidence of interaction between genotype and smoking (P = 0.007). In women who smoked during pregnancy, each additional smoking-related T-allele was associated with a 20 g [95% confidence interval (95% CI): 4–36 g] lower birth weight (P = 0.014). However, in women who did not smoke during pregnancy, the effect size estimate was 5 g per T-allele (95% CI: −4 to 14 g; P = 0.268). To conclude, smoking status during pregnancy modifies the association between maternal rs1051730 genotype and offspring birth weight. This strengthens the evidence that smoking during pregnancy is causally related to lower offspring birth weight and suggests that population interventions that effectively reduce smoking in pregnant women would result in a reduced prevalence of low birth weight.

Onset of alcohol use at an early age increases the risk for later alcohol dependence. We investigated the role of the glucocorticoid receptor (GR) gene (NR3C1) in onset of alcohol use and abuse in 14-year-old adolescents (n = 4534). Several NR3C1 polymorphisms were associated with onset of alcohol drinking or drunkenness at this age. Strongest associations were observed in females, with one marker (rs244465) remaining significant after correction for multiple testing (Padj = 0.0067; odds ratio = 1.7, for drunkenness). Our data provide the first evidence that GR modulates initiation of alcohol abuse and reveal a polymorphism that might contribute to susceptibility to addiction.

Atopic dermatitis (AD) is a common chronic skin disease with high heritability. Apart from filaggrin (FLG), the genes influencing AD are largely unknown. We conducted a genome-wide association meta-analysis of 5,606 cases and 20,565 controls from 16 population-based cohorts and followed up the ten most strongly associated novel markers in a further 5,419 cases and 19,833 controls from 14 studies. Three SNPs met genome-wide significance in the discovery and replication cohorts combined: rs479844 upstream of OVOL1 (OR=0.88, p=1.1×10−13) and rs2164983 near ACTL9 (OR=1.16, p=7.1×10−9), genes which have been implicated in epidermal proliferation and differentiation, as well as rs2897442 in KIF3A within the cytokine cluster on 5q31.1 (OR=1.11, p=3.8×10−8). We also replicated the FLG locus and two recently identified association signals at 11q13.5 (rs7927894, p=0.008) and 20q13.3 (rs6010620, p=0.002). Our results underline the importance of both epidermal barrier function and immune dysregulation in AD pathogenesis.

Stature is a classical and highly heritable complex trait, with 80%–90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ2 = 83.89, df = 1; p = 5.2×10−20). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

Author Summary

Studies investigating the extent to which genetics influences human characteristics such as height have concentrated mainly on common variants of genes, where having one or two copies of a given variant influences the trait or risk of disease. This study explores whether a different type of genetic variant might also be important. We investigate the role of recessive genetic variants, where two identical copies of a variant are required to have an effect. By measuring genome-wide homozygosity—the phenomenon of inheriting two identical copies at a given point of the genome—in 35,000 individuals from 21 European populations, and by comparing this to individual height, we found that the more homozygous the genome, the shorter the individual. The offspring of first cousins (who have increased homozygosity) were predicted to be up to 3 cm shorter on average than the offspring of unrelated parents. Height is influenced by the combined effect of many recessive variants dispersed across the genome. This may also be true for other human characteristics and diseases, opening up a new way to understand how genetic variation influences our health.

Background

We previously identified via a genome wide association study variants near LEKR and CCNL1 and in the ADCY5 genes lead to lower birthweight. Here, we study the impact of these variants and social stress during pregnancy, defined as social adversity and neighborhood disparity, on infant birth size. We aimed to determine whether the addition of genetic variance magnified the observed associations.

Methodology/Principal Findings

We analyzed data from the Northern Finland Birth Cohort 1986 (n = 5369). Social adversity was defined by young maternal age (<20 years), low maternal education (<11 years), and/or single marital status. Neighborhood social disparity was assessed by discrepancy between neighborhoods relative to personal socio-economic status. These variables are indicative of social and socioeconomic stress, but also of biological risk. The adjusted multiple regression analysis showed smaller birth size in both infants of mothers who experienced social adversity (birthweight by −40.4 g, 95%CI −61.4, −19.5; birth length −0.14 cm, 95%CI −0.23, −0.05; head circumference −0.09 cm 95%CI −0.15, −0.02) and neighborhood disparity (birthweight −28.8 g, 95%CI −47.7, −10.0; birth length −0.12 cm, 95%CI −0.20, −0.05). The birthweight-lowering risk allele (SNP rs900400 near LEKR and CCNL1) magnified this association in an additive manner. However, likely due to sample size restriction, this association was not significant for the SNP rs9883204 in ADCY5. Birth size difference due to social stress was greater in the presence of birthweight-lowering alleles.

Conclusions/Significance

Social adversity, neighborhood disparity, and genetic variants have independent associations with infant birth size in the mutually adjusted analyses. If the newborn carried a risk allele rs900400 near LEKR/CCNL1, the impact of stress on birth size was stronger. These observations give support to the hypothesis that individuals with genetic or other biological risk are more vulnerable to environmental influences. Our study indicates the need for further research to understand the mechanisms by which genes impact individual vulnerability to environmental insults.

Background

CHRNA5-CHRNA3-CHRNB4 and TTC12-ANKK1-DRD2 gene-clusters influence smoking behavior. Our aim was to test developmental changes in their effects as well as the interplays between them and with non-genetic factors.

Methods

Participants included 4,762 subjects from a general population based prospective Northern Finland 1966 Birth Cohort (NFBC 1966). Smoking behavior was collected at age 14 and 31 years(y). Information on maternal smoking, socio-economic status, and novelty seeking were also collected. Structural equation modeling was used to construct an integrative etiological model including genetic and non-genetic factors.

Results

Several SNPs in both gene-clusters were significantly associated with smoking. The most significant were in CHRNA3 (rs1051730, P=1.1×10−5) and in TTC12 (rs10502172, P=9.1×10−6). CHRNA3-rs1051730[A] was more common among heavy/regular smokers than non-smokers with similar effect-sizes at age 14y [OR(95%CI):1.27(1.06–1.52)] and 31y [1.28(1.13–1.44)]. TTC12-rs10502172[G] was more common among smokers than non-smokers with stronger association at 14y [1.33(1.11–1.60)] than 31y [1.14(1.02–1.28)]. In adolescence, carriers of three-four risk alleles at either CHRNA3-rs1051730 or TTC12-rs10502172 had almost 3-fold odds of smoking regularly than subjects with no risk alleles. TTC12-rs10502172 effect on smoking in adulthood was mediated by its effect on smoking in adolescence and via novelty seeking. Effect of CHRNA3-rs1051730 on smoking in adulthood was direct.

Conclusions

TTC12-ANKK1-DRD2s seemed to influence smoking behavior mainly in adolescence and its effect is partially mediated by personality characteristics promoting drug-seeking behavior. In contrast, CHRNA5-CHRNA3-CHRNB4 is involved in the transition towards heavy smoking in mid-adulthood and in smoking persistence. Factors related to familial and social disadvantages were strong independent predictors of smoking.

Ruth Loos and colleagues report findings from a meta-analysis of multiple studies examining the extent to which physical activity attenuates effects of a specific gene variant, FTO, on obesity in adults and children. They report a fairly substantial attenuation by physical activity on the effects of this genetic variant on the risk of obesity in adults.

Background

The FTO gene harbors the strongest known susceptibility locus for obesity. While many individual studies have suggested that physical activity (PA) may attenuate the effect of FTO on obesity risk, other studies have not been able to confirm this interaction. To confirm or refute unambiguously whether PA attenuates the association of FTO with obesity risk, we meta-analyzed data from 45 studies of adults (n = 218,166) and nine studies of children and adolescents (n = 19,268).

Methods and Findings

All studies identified to have data on the FTO rs9939609 variant (or any proxy [r2>0.8]) and PA were invited to participate, regardless of ethnicity or age of the participants. PA was standardized by categorizing it into a dichotomous variable (physically inactive versus active) in each study. Overall, 25% of adults and 13% of children were categorized as inactive. Interaction analyses were performed within each study by including the FTO×PA interaction term in an additive model, adjusting for age and sex. Subsequently, random effects meta-analysis was used to pool the interaction terms. In adults, the minor (A−) allele of rs9939609 increased the odds of obesity by 1.23-fold/allele (95% CI 1.20–1.26), but PA attenuated this effect (pinteraction  = 0.001). More specifically, the minor allele of rs9939609 increased the odds of obesity less in the physically active group (odds ratio  = 1.22/allele, 95% CI 1.19–1.25) than in the inactive group (odds ratio  = 1.30/allele, 95% CI 1.24–1.36). No such interaction was found in children and adolescents.

Conclusions

The association of the FTO risk allele with the odds of obesity is attenuated by 27% in physically active adults, highlighting the importance of PA in particular in those genetically predisposed to obesity.

Please see later in the article for the Editors' Summary

Editors’ Summary

Background

Two in three Americans are overweight, of whom half are obese, and the trend towards increasing obesity is now seen across developed and developing countries. There has long been interest in understanding the impact of genes and environment when it comes to apportioning responsibility for obesity. Carrying a change in the FTO gene is common (found in three-quarters of Europeans and North Americans) and is associated with a 20%–30% increased risk of obesity. Some overweight or obese individuals may feel that the dice are loaded and there is little point in fighting the fat; it has been reported that those made aware of their genetic susceptibility to obesity may still choose a poor diet. A similar fatalism may occur when overweight and obese people consider physical activity. But disentangling the influence of physical activity on those genetically susceptible to obesity from other factors that might impact weight is not straightforward, as it requires large sample sizes, could be subject to publication bias, and may rely on less than ideal self-reporting methods.

Why Was This Study Done?

The public health ramifications of understanding the interaction between genetic susceptibility to obesity and physical activity are considerable. Tackling the rising prevalence of obesity will inevitably include interventions principally aimed at changing dietary intake and/or increasing physical activity, but the evidence for these with regards to those genetically susceptible has been lacking to date. The authors of this paper set out to explore the interaction between the commonest genetic susceptibility trait and physical activity using a rigorous meta-analysis of a large number of studies.

What Did the Researchers Do and Find?

The authors were concerned that a meta-analysis of published studies would be limited both by the data available to them and by possible bias. Instead of this more widely used approach, they took the literature search as their starting point, identified other studies through their collaborators’ network, and then undertook a meta-analysis of all available studies using a new and standardized analysis plan. This entailed an extremely large number of authors mining their data afresh to extract the relevant data points to enable such a meta-analysis. Physical activity was identified in the original studies in many different ways, including by self-report or by using an external measure of activity or heart rate. In order to perform the meta-analysis, participants were labeled as physically active or inactive in each study. For studies that had used a continuous scale, the authors decided that the bottom 20% of the participants were inactive (10% for children and adolescents). Using data from over 218,000 adults, the authors found that carrying a copy of the susceptibility gene increased the odds of obesity by 1.23-fold. But the size of this influence was 27% less in the genetically susceptible adults who were physically active (1.22-fold) compared to those who were physically inactive (1.30-fold). In a smaller study of about 19,000 children, no such effect of physical activity was seen.

What Do these Findings Mean?

This study demonstrates that people who carry the susceptibility gene for obesity can benefit from physical activity. This should inform health care professionals and the wider public that the view of genetically determined obesity not being amenable to exercise is incorrect and should be challenged. Dissemination, implementation, and ensuring uptake of effective physical activity programs remains a challenge and deserves further consideration. That the researchers treated “physically active” as a yes/no category, and how they categorized individuals, could be criticized, but this was done for pragmatic reasons, as a variety of means of assessing physical activity were used across the studies. It is unlikely that the findings would have changed if the authors had used a different method of defining physically active. Most of the studies included in the meta-analysis looked at one time point only; information about the influence of physical activity on weight changes over time in genetically susceptible individuals is only beginning to emerge.

Additional Information

Please access these websites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.1001116.

This study is further discussed in a PLoS Medicine Perspective by Lennert Veerman

The US Centers for Disease Control and Prevention provides obesity-related statistics, details of prevention programs, and an overview on public health strategy in the United States

A more worldwide view is given by the World Health Organization

The UK National Health Service website gives information on physical activity guidelines for different age groups, while similar information can also be found from US sources

An age-dependent association between variation at the FTO locus and BMI in children has been suggested. We meta-analyzed associations between the FTO locus (rs9939609) and BMI in samples, aged from early infancy to 13 years, from 8 cohorts of European ancestry. We found a positive association between additional minor (A) alleles and BMI from 5.5 years onwards, but an inverse association below age 2.5 years. Modelling median BMI curves for each genotype using the LMS method, we found that carriers of minor alleles showed lower BMI in infancy, earlier adiposity rebound (AR), and higher BMI later in childhood. Differences by allele were consistent with two independent processes: earlier AR equivalent to accelerating developmental age by 2.37% (95% CI 1.87, 2.87, p = 10−20) per A allele and a positive age by genotype interaction such that BMI increased faster with age (p = 10−23). We also fitted a linear mixed effects model to relate genotype to the BMI curve inflection points adiposity peak (AP) in infancy and AR. Carriage of two minor alleles at rs9939609 was associated with lower BMI at AP (−0.40% (95% CI: −0.74, −0.06), p = 0.02), higher BMI at AR (0.93% (95% CI: 0.22, 1.64), p = 0.01), and earlier AR (−4.72% (−5.81, −3.63), p = 10−17), supporting cross-sectional results. Overall, we confirm the expected association between variation at rs9939609 and BMI in childhood, but only after an inverse association between the same variant and BMI in infancy. Patterns are consistent with a shift on the developmental scale, which is reflected in association with the timing of AR rather than just a global increase in BMI. Results provide important information about longitudinal gene effects and about the role of FTO in adiposity. The associated shifts in developmental timing have clinical importance with respect to known relationships between AR and both later-life BMI and metabolic disease risk.

Author Summary

Variation at the FTO locus is reliably associated with BMI and adiposity-related traits, but little is still known about the effects of variation at this gene, particularly in children. We have examined a large collection of samples for which both genotypes at rs9939609 and multiple measurements of BMI are available. We observe a positive association between the minor allele (A) at rs9939609 and BMI emerging in childhood that has the characteristics of a shift in the age scale leading simultaneously to lower BMI during infancy and higher BMI in childhood. Assessed in cross section and longitudinally, we find evidence of variation at rs9939609 being associated with the timing of AR and the concert of events expected with such a change to the BMI curve. Importantly, the apparently negative association between the minor allele (A) and BMI in early life, which is then followed by an earlier AR and greater BMI in childhood, is a pattern known to be associated with both the risk of adult BMI and metabolic disorders such as type 2 diabetes (T2D). These findings are important in our understanding of the contribution of FTO to adiposity, but also in light of efforts to appreciate genetic effects in a lifecourse context.

OBJECTIVE

A1C is widely considered the gold standard for monitoring effective blood glucose levels. Recently, a genome-wide association study reported an association between A1C and rs7072268 within HK1 (encoding hexokinase 1), which catalyzes the first step of glycolysis. HK1 deficiency in erythrocytes (red blood cells [RBCs]) causes severe nonspherocytic hemolytic anemia in both humans and mice.

RESEARCH DESIGN AND METHODS

The contribution of rs7072268 to A1C and the RBC-related traits was assessed in 6,953 nondiabetic European participants. We additionally analyzed the association with hematologic traits in 5,229 nondiabetic European individuals (in whom A1C was not measured) and 1,924 diabetic patients. Glucose control–related markers other than A1C were analyzed in 18,694 nondiabetic European individuals. A type 2 diabetes case-control study included 7,447 French diabetic patients.

RESULTS

Our study confirms a strong association between the rs7072268–T allele and increased A1C (β = 0.029%; P = 2.22 × 10−7). Surprisingly, despite adequate study power, rs7072268 showed no association with any other markers of glucose control (fasting- and 2-h post-OGTT–related parameters, n = 18,694). In contrast, rs7072268–T allele decreases hemoglobin levels (n = 13,416; β = −0.054 g/dl; P = 3.74 × 10−6) and hematocrit (n = 11,492; β = −0.13%; P = 2.26 × 10−4), suggesting a proanemic effect. The T allele also increases risk for anemia (836 cases; odds ratio 1.13; P = 0.018).

CONCLUSIONS

HK1 variation, although strongly associated with A1C, does not seem to be involved in blood glucose control. Since HK1 rs7072268 is associated with reduced hemoglobin levels and favors anemia, we propose that HK1 may influence A1C levels through its anemic effect or its effect on glucose metabolism in RBCs. These findings may have implications for type 2 diabetes diagnosis and clinical management because anemia is a frequent complication of the diabetes state.

INTRODUCTORY PARAGRAPH

To identify genetic variants associated with birth weight, we meta-analyzed six genome-wide association (GWA) studies (N=10,623 Europeans from pregnancy/birth cohorts) and followed up two lead signals in thirteen replication studies (N=27,591). Rs900400 near LEKR1 and CCNL1 (P=2×10−35), and rs9883204 in ADCY5 (P=7×10−15) were robustly associated with birth weight. Correlated SNPs in ADCY5 were recently implicated in regulation of glucose levels and type 2 diabetes susceptibility,1 providing evidence that the well described association between lower birth weight and subsequent type 2 diabetes2,3 has a genetic component, distinct from the proposed role of programming by maternal nutrition. Using data from both SNPs, the 9% of Europeans with 4 birth weight-lowering alleles were, on average, 113g (95%CI 89-137g) lighter at birth than the 24% with 0 or 1 allele (Ptrend=7×10−30). The impact on birth weight is similar to that of a mother smoking 4-5 cigarettes per day in the third trimester of pregnancy.4

The association between variation in the fat mass and obesity-associated (FTO) gene and adulthood body mass index (BMI; weight (kg)/height (m)2) is well-replicated. More thorough analyses utilizing phenotypic data over the life course may deepen our understanding of the development of BMI and thus help in the prevention of obesity. The authors used a structural equation modeling approach to explore the network of variables associated with BMI from the prenatal period to age 31 years (1965–1997) in 4,435 subjects from the Northern Finland Birth Cohort 1966. The use of structural equation modeling permitted the easy inclusion of variables with missing values in the analyses without separate imputation steps, as well as differentiation between direct and indirect effects. There was an association between the FTO single nucleotide polymorphism rs9939609 and BMI at age 31 years that persisted after controlling for several relevant factors during the life course. The total effect of the FTO variant on adult BMI was mostly composed of the direct effect, but a notable part was also arising indirectly via its effects on earlier BMI development. In addition to well-established genetic determinants, many life-course factors such as physical activity, in spite of not showing mediation or interaction, had a strong independent effect on BMI.

Tooth development is a highly heritable process which relates to other growth and developmental processes, and which interacts with the development of the entire craniofacial complex. Abnormalities of tooth development are common, with tooth agenesis being the most common developmental anomaly in humans. We performed a genome-wide association study of time to first tooth eruption and number of teeth at one year in 4,564 individuals from the 1966 Northern Finland Birth Cohort (NFBC1966) and 1,518 individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC). We identified 5 loci at P<5×10−8, and 5 with suggestive association (P<5×10−6). The loci included several genes with links to tooth and other organ development (KCNJ2, EDA, HOXB2, RAD51L1, IGF2BP1, HMGA2, MSRB3). Genes at four of the identified loci are implicated in the development of cancer. A variant within the HOXB gene cluster associated with occlusion defects requiring orthodontic treatment by age 31 years.

Author Summary

Genome-wide association studies have been used to identify genetic variants conferring susceptibility to diseases, intermediate phenotypes, and physiological traits such as height, hair color, and age at menarche. Here we analyze the NFBC1966 and ALSPAC birth cohorts to investigate the genetic determinants of a key developmental process: primary tooth development. The prospective nature of our studies allows us to exploit accurate measurements of age at first tooth eruption and number of teeth at one year, and also provides the opportunity to assess whether genetic variants affecting these traits are associated with dental problems later in the life course. Of the genes that we find to be associated with primary tooth development, several have established roles in tooth development and growth, and almost half have proposed links with the development of cancer. We find that one of the variants is also associated with occlusion defects requiring orthodontic treatment later in life. Our findings should provide a strong foundation for the study of the genetic architecture of tooth development, which as well as its relevance to medicine and dentistry, may have implications in evolutionary biology since teeth represent important markers of evolution.

To identify novel genetic loci associated with fasting glucose concentrations, we examined the leading association signals in 10 genome-wide association scans involving a total of 36,610 individuals of European descent. Variants in the gene encoding the melatonin receptor 1B (MTNR1B) were consistently associated with fasting glucose across all ten studies. The strongest signal was observed at rs10830963, where each G-allele (frequency 0.30 in HapMap CEU) was associated with an increase of 0.07 (95%CI 0.06–0.08) mmol/L in fasting glucose levels (P=3.2×10−50) and reduced beta-cell function as measured by homeostasis model assessment (HOMA-B, P=1.1×10−15). The same allele was associated with an increased risk of type 2 diabetes (odds ratio = 1.09 (1.05–1.12), per G allele P=3.3×10−7) in a meta-analysis of thirteen case-control studies totalling 18,236 cases and 64,453 controls. Our analyses also confirm previous associations of fasting glucose with variants at the G6PC2 (rs560887, P=1.1×10−57) and GCK (rs4607517, P=1.0×10−25) loci.

The INSIG2 rs7566605 polymorphism was identified for obesity (BMI≥30 kg/m2) in one of the first genome-wide association studies, but replications were inconsistent. We collected statistics from 34 studies (n = 74,345), including general population (GP) studies, population-based studies with subjects selected for conditions related to a better health status (‘healthy population’, HP), and obesity studies (OB). We tested five hypotheses to explore potential sources of heterogeneity. The meta-analysis of 27 studies on Caucasian adults (n = 66,213) combining the different study designs did not support overall association of the CC-genotype with obesity, yielding an odds ratio (OR) of 1.05 (p-value = 0.27). The I2 measure of 41% (p-value = 0.015) indicated between-study heterogeneity. Restricting to GP studies resulted in a declined I2 measure of 11% (p-value = 0.33) and an OR of 1.10 (p-value = 0.015). Regarding the five hypotheses, our data showed (a) some difference between GP and HP studies (p-value = 0.012) and (b) an association in extreme comparisons (BMI≥32.5, 35.0, 37.5, 40.0 kg/m2 versus BMI<25 kg/m2) yielding ORs of 1.16, 1.18, 1.22, or 1.27 (p-values 0.001 to 0.003), which was also underscored by significantly increased CC-genotype frequencies across BMI categories (10.4% to 12.5%, p-value for trend = 0.0002). We did not find evidence for differential ORs (c) among studies with higher than average obesity prevalence compared to lower, (d) among studies with BMI assessment after the year 2000 compared to those before, or (e) among studies from older populations compared to younger. Analysis of non-Caucasian adults (n = 4889) or children (n = 3243) yielded ORs of 1.01 (p-value = 0.94) or 1.15 (p-value = 0.22), respectively. There was no evidence for overall association of the rs7566605 polymorphism with obesity. Our data suggested an association with extreme degrees of obesity, and consequently heterogeneous effects from different study designs may mask an underlying association when unaccounted for. The importance of study design might be under-recognized in gene discovery and association replication so far.

Author Summary

A polymorphism of the INSIG2 gene was identified as being associated with obesity in one of the first genome-wide association studies. However, this association has since then been highly debated upon inconsistent subsequent reports. We collected association information from 34 studies including a total of 74,000 participants. In a meta-analysis of the 27 studies including 66,000 Caucasian adults, we found no overall association of this polymorphism rs7566605 with obesity, comparing subjects with a body-mass-index (BMI)≥30 kg/m2 with normal BMI subjects (BMI<30 kg/m2). Our data suggested an association of this polymorphism with extreme obesity (e.g., BMI≥37.5 kg/m2) compared to normal controls. Such an association with extreme obesity might induce heterogeneous effects from different study designs depending on the proportion of extreme obesity included by the design. However, further studies would be required to substantiate this finding. The importance of study design might be under-recognized in gene discovery and association replication so far.

Genome-wide association studies (GWAS) of longitudinal birth cohorts enable joint investigation of environmental and genetic influences on complex traits. We report GWAS results for nine quantitative metabolic traits (triglycerides, high-density lipoprotein, low-density lipoprotein, glucose, insulin, C-reactive protein, body mass index, and systolic and diastolic blood pressure) in the Northern Finland Birth Cohort 1966 (NFBC1966), drawn from the most genetically isolated Finnish regions. We replicate most previously reported associations for these traits and identify nine new associations, several of which highlight genes with metabolic functions: high-density lipoprotein with NR1H3 (LXRA), low-density lipoprotein with AR and FADS1-FADS2, glucose with MTNR1B, and insulin with PANK1. Two of these new associations emerged after adjustment of results for body mass index. Gene-environment interaction analyses suggested additional associations, which will require validation in larger samples. The currently identified loci, together with quantified environmental exposures, explain little of the trait variation in NFBC1966. The association observed between low-density lipoprotein and an infrequent variant in AR suggests the potential of such a cohort for identifying associations with both common, low-impact and rarer, high-impact quantitative trait loci.

Recent genome-wide association (GWA) studies have identified dozens of common variants associated with adult height. However, it is unknown how these variants influence height growth during childhood. We derived peak height velocity in infancy (PHV1) and puberty (PHV2) and timing of pubertal height growth spurt from parametric growth curves fitted to longitudinal height growth data to test their association with known height variants. The study consisted of N = 3,538 singletons from the prospective Northern Finland Birth Cohort 1966 with genotype data and frequent height measurements (on average 20 measurements per person) from 0–20 years. Twenty-six of the 48 variants tested associated with adult height (p<0.05, adjusted for sex and principal components) in this sample, all in the same direction as in previous GWA scans. Seven SNPs in or near the genes HHIP, DLEU7, UQCC, SF3B4/SV2A, LCORL, and HIST1H1D associated with PHV1 and five SNPs in or near SOCS2, SF3B4/SV2A, C17orf67, CABLES1, and DOT1L with PHV2 (p<0.05). We formally tested variants for interaction with age (infancy versus puberty) and found biologically meaningful evidence for an age-dependent effect for the SNP in SOCS2 (p = 0.0030) and for the SNP in HHIP (p = 0.045). We did not have similar prior evidence for the association between height variants and timing of pubertal height growth spurt as we had for PHVs, and none of the associations were statistically significant after correction for multiple testing. The fact that in this sample, less than half of the variants associated with adult height had a measurable effect on PHV1 or PHV2 is likely to reflect limited power to detect these associations in this dataset. Our study is the first genetic association analysis on longitudinal height growth in a prospective cohort from birth to adulthood and gives grounding for future research on the genetic regulation of human height during different periods of growth.

Author Summary

Family studies have shown that adult height is largely genetically determined. Identification of common genetic factors has been expedited with recent advances in genotyping techniques. However, factors regulating childhood height growth remain unclear. We investigated genetic variants of adult height for associations with peak height velocity in infancy (PHV1) and puberty (PHV2) and timing of pubertal growth spurt in a population based sample of 3,538 Finns born in 1966. Most variants studied associated with adult height in this sample. Of the 48 genetic variants tested, seven of them associated with PHV1 and five with PHV2. However, only one of these associated with both, and we found suggestive evidence for differential effects at different stages of growth for some of the variants. In this sample, less than half of the variants associated with adult height had a measurable effect on PHV1 or PHV2. However, these differences may reflect lower statistical power to detect associations with height velocities compared to adult height. This study provides a foundation for further biological investigation into the genes acting at each stage of height growth.