Genetic Ancestry, Skin Reflectance and Pigmentation Genotypes in Association with Serum Vitamin D Metabolite Balance

Cross-ancestry analyses identify new genetic loci associated with 25-hydroxyvitamin D

Vitamin D status-a complex trait influenced by environmental and genetic factors-is tightly associated with skin colour and ancestry. Yet very few studies have investigated the genetic underpinnings of vitamin D levels across diverse ancestries, and the ones that have, relied on small sample sizes, resulting in inconclusive results. Here, we conduct genome-wide association studies (GWAS) of 25 hydroxyvitamin D (25OHD)-the main circulating form of vitamin D-in 442,435 individuals from four broad genetically-determined ancestry groups represented in the UK Biobank: European (N = 421,867), South Asian (N = 9,983), African (N = 8,306) and East Asian (N = 2,279). We identify a new genetic determinant of 25OHD (rs146759773) in individuals of African ancestry, which was not detected in previous analysis of much larger European cohorts due to low minor allele frequency. We show genome-wide significant evidence of dominance effects in 25OHD that protect against vitamin D deficiency. Given that key events in the synthesis of 25OHD occur in the skin and are affected by pigmentation levels, we conduct GWAS of 25OHD stratified by skin colour and identify new associations. Lastly, we test the interaction between skin colour and variants associated with variance in 25OHD levels and identify two loci (rs10832254 and rs1352846) whose association with 25OHD differs in individuals of distinct complexions. Collectively, our results provide new insights into the complex relationship between 25OHD and skin colour and highlight the importance of diversity in genomic studies. Despite the much larger rates of vitamin D deficiency that we and others report for ancestry groups with dark skin (e.g., South Asian), our study highlights the importance of considering ancestral background and/or skin colour when assessing the implications of low vitamin D.

What can Mendelian randomization contribute to biological anthropology?

Uncovering causal relationships between exposures and outcomes can be difficult in observational studies because of the potential for confounding and reverse causation to produce biased estimates. Conversely, randomized controlled trials (RCTs) provide the strongest evidence for causality but they are not always feasible. Mendelian randomization (MR) is a method that aims to strengthen causal inference using genetic variants as proxies or instrumental variables (IVs) for exposures, to overcome the above-mentioned biases. Since allele segregation occurs at random from parents to offspring, and alleles for a trait assort independently from those for other traits, MR studies have frequently been compared to "natural" RCTs. In biological anthropology (BA) relationships between variables of interest are usually evaluated using observational data, often remaining descriptive, and other approaches to causal inference have seldom been implemented. Here, we propose the use of MR to investigate cause and effect relationships in BA studies and provide examples to show how that can be done across areas of BA relevance, such as adaptation to the environment, nutrition and life history theory. While we consider MR a useful addition to the biological anthropologist's toolbox, we advocate the adoption of a wide range of methods, affected by different types of biases, in order to better answer the important causal questions for the discipline.

25-Hydroxyvitamin D Plasma Levels in Natural Populations of Pigmented and Partially Pigmented Land Iguanas from Galápagos (Conolophus spp.)

We report the first data on 25-hydroxyvitamin D plasma levels in natural populations of three species of land iguana endemic to the Galápagos Islands (Conolophus marthae, C. subcristatus, and C. pallidus). The pigment is present throughout the whole body in the skin of C. subcristatus and C. pallidus. On the contrary, pigment is not present in the skin of an extended part of the body in C. marthae. The only existing population of C. marthae is syntopic with a population of C. subcristatus, and the two species are closely related. These circumstances would suggest that, under the assumption that the species show a similar basking behavior and in the absence of compensatory mechanisms, lighter pigmentation should favor higher vitamin D levels. Thus, C. marthae, compared with C. subcristatus in Wolf Volcano, could show higher levels of 25(OH)D plasma levels, or equal, if compensatory mechanisms exist. The three species showed levels in the range of average values for healthy iguanas. However, contrary to the expectation, C. marthae consistently exhibited the lowest 25(OH)D plasma levels. We discuss possible factors affecting vitamin concentration and hypothesize that C. marthae may use the habitat to limit exposure to the high UVB irradiation at Wolf Volcano.

Ancestry-Adjusted Vitamin D Metabolite Concentrations in Association With Cytochrome P450 3A Polymorphisms

We investigated the association between genetic polymorphisms in cytochrome P450 (CYP2R1, CYP24A1, and the CYP3A family) with nonsummer plasma concentrations of vitamin D metabolites (25-hydroxyvitamin D3 (25(OH)D3) and proportion 24,25-dihydroxyvitamin D3 (24,25(OH)2D3)) among healthy individuals of sub-Saharan African and European ancestry, matched on age (within 5 years; n = 188 in each ancestral group), in central suburban Pennsylvania (2006-2009). Vitamin D metabolites were measured using high-performance liquid chromatography with tandem mass spectrometry. Paired multiple regression and adjusted least-squares mean analyses were used to test for associations between genotype and log-transformed metabolite concentrations, adjusted for age, sex, proportion of West-African genetic ancestry, body mass index, oral contraceptive (OC) use, tanning bed use, vitamin D intake, days from summer solstice, time of day of blood draw, and isoforms of the vitamin D receptor (VDR) and vitamin D binding protein. Polymorphisms in CYP2R1, CYP3A43, vitamin D binding protein, and genetic ancestry proportion remained associated with plasma 25(OH)D3 after adjustment. Only CYP3A43 and VDR polymorphisms were associated with proportion 24,25(OH)2D3. Magnitudes of association with 25(OH)D3 were similar for CYP3A43, tanning bed use, and OC use. Significant least-squares mean interactions (CYP2R1/OC use (P = 0.030) and CYP3A43/VDR (P = 0.013)) were identified. A CYP3A43 genotype, previously implicated in cancer, is strongly associated with biomarkers of vitamin D metabolism. Interactive associations should be further investigated.