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

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.

Author Summary

Genetic information is transmitted to the cell only through RNA molecules. A special class of RNAs is comprised of the small (up to 30 nucleotide) ones, known to be potent regulators of various cellular processes. At the same time, they have not been as widely studied as messenger RNAs—we do not know how much variation in their sequence and expression level occurs naturally in human populations or how this variability influences other traits. We measured small RNA levels and genetic variability in fat tissue from 131 individuals by high-throughput sequencing. We could associate the expression levels with genetic background of the individuals, as well as changes in metabolic traits. Surprisingly, we found no large scale influence of small RNA variation on mRNA levels, their main regulatory target. Overall, our study is the first to give a quantitative picture of the naturally occurring variation in these important regulatory molecules in human fat tissue.

Genome-wide association studies have identified many genetic variants associated with complex traits. However, at only a minority of loci have the molecular mechanisms mediating these associations been characterized. In parallel, whilst cis-regulatory patterns of gene expression have been extensively explored, the identification of trans-regulatory effects in humans has attracted less attention. We demonstrate that the Type 2 diabetes and HDL-cholesterol associated cis-acting eQTL of the maternally-expressed transcription factor KLF14 acts as a master trans-regulator of adipose gene expression. Expression levels of genes regulated by this trans-eQTL are highly-correlated with concurrently-measured metabolic traits, and a subset of the trans-genes harbor variants directly-associated with metabolic phenotypes. This trans-eQTL network provides a mechanistic understanding of the effect of the KLF14 locus on metabolic disease risk, providing a potential model for other complex traits.

Background

Vitamin D is crucial for maintaining musculoskeletal health. Recently, vitamin D insufficiency has been linked to a number of extraskeletal disorders, including diabetes, cancer, and cardiovascular disease. Determinants of circulating 25-hydroxyvitamin D (25-OH D) include sun exposure and dietary intake, but its high heritability suggests that genetic determinants may also play a role.

Methods

We performed a genome-wide association study of 25-OH D among ∼30,000 individuals of European descent from 15 cohorts. Five cohorts were designated as discovery cohorts (n=16,125), five as in silico replication cohorts (n=9,366), and five as de novo replication cohorts (n=8,378). Association results were combined using z-score-weighted meta-analysis. Vitamin D insufficiency was defined as 25-OH D <75 nmol/L or <50 nmol/L.

Findings

Variants at three loci reached genome-wide significance in the discovery cohorts, and were confirmed in the replication cohorts: 4p12 (overall P=1.9 × 10-109 for rs2282679, in GC); 11q12 (P=2.1 × 10-27 for rs12785878, near DHCR7); 11p15 (P=3.3 × 10-20 for rs10741657, near CYP2R1). Variants at an additional locus (20q13, CYP24A1) were genome-wide significant in the pooled sample (P=6.0 × 10-10 for rs6013897). A genotype score was constructed using the three confirmed variants. Those in the top quartile of genotype scores had 2- to 2.5-fold elevated odds of vitamin D insufficiency (P≤1 × 10-26).

Interpretation

Variants near genes involved in cholesterol synthesis (DHCR7), hydroxylation (CYP2R1, CYP24A1), and vitamin D transport (GC) influence vitamin D status. Genetic variation at these loci identifies individuals of European descent who have substantially elevated risk of vitamin D insufficiency.

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Author Summary

Regulation of gene expression is a fundamental cellular process determining a large proportion of the phenotypic variance. Previous studies have identified genetic loci influencing gene expression levels (eQTLs), but the complexity of their tissue-specific properties has not yet been well-characterized. In this study, we perform cis-eQTL analysis in a unique matched co-twin design for three human tissues derived simultaneously from the same set of individuals. The study design allows validation of the substantial discoveries we make in each tissue. We explore in depth the tissue-dependent features of regulatory variants and estimate the proportions of shared and specific effects. We use continuous measures of eQTL sharing to circumvent the statistical power limitations of comparing direct overlap of eQTLs in multiple tissues. In this framework, we demonstrate that 30% of eQTLs are shared among tissues, while 29% are exclusively tissue-specific. Furthermore, we show that the fold change in expression between eQTL genotypic classes differs between tissues. Even among shared eQTLs, we report a substantial proportion (10%–20%) of significant tissue differences in magnitude of these effects. The complexities we highlight here are essential for understanding the impact of regulatory variants on complex traits.

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

Author Summary

We measured human eye color to hue and saturation values from high-resolution, digital, full-eye photographs of several thousand Dutch Europeans. This quantitative approach, which is extremely cost-effective, portable, and time efficient, revealed that human eye color varies along more dimensions than the one represented by the blue-green-brown categories studied previously. Our work represents the first genome-wide study of quantitative human eye color. We clearly identified 3 new loci, LYST, 17q25.3, TTC3/DSCR9, in contributing to the natural and subtle eye color variation along multiple dimensions, providing new leads towards a more detailed understanding of the genetic basis of human eye color. Our quantitative prediction model explained over 50% of eye color variance, representing the highest accuracy achieved so far in genomic prediction of human complex and quantitative traits, with relevance for future forensic applications.