Anna-Liisa Hartikainen

Common genetic determinants of vitamin D insufficiency: a genome-wide association study

BACKGROUND Vitamin D is crucial for maintenance of musculoskeletal health, and might also have a role in extraskeletal tissues. Determinants of circulating 25-hydroxyvitamin D concentrations include sun exposure and diet, but high heritability suggests that genetic factors could also play a part. We aimed to identify common genetic variants affecting vitamin D concentrations and risk of insufficiency. METHODS We undertook a genome-wide association study of 25-hydroxyvitamin D concentrations in 33 996 individuals of European descent from 15 cohorts. Five epidemiological cohorts were designated as discovery cohorts (n=16 125), five as in-silico replication cohorts (n=9367), and five as de-novo replication cohorts (n=8504). 25-hydroxyvitamin D concentrations were measured by radioimmunoassay, chemiluminescent assay, ELISA, or mass spectrometry. Vitamin D insufficiency was defined as concentrations lower than 75 nmol/L or 50 nmol/L. We combined results of genome-wide analyses across cohorts using Z-score-weighted meta-analysis. Genotype scores were constructed for confirmed variants. FINDINGS Variants at three loci reached genome-wide significance in discovery cohorts for association with 25-hydroxyvitamin D concentrations, and were confirmed in replication cohorts: 4p12 (overall p=1.9x10(-109) for rs2282679, in GC); 11q12 (p=2.1x10(-27) for rs12785878, near DHCR7); and 11p15 (p=3.3x10(-20) for rs10741657, near CYP2R1). Variants at an additional locus (20q13, CYP24A1) were genome-wide significant in the pooled sample (p=6.0x10(-10) for rs6013897). Participants with a genotype score (combining the three confirmed variants) in the highest quartile were at increased risk of having 25-hydroxyvitamin D concentrations lower than 75 nmol/L (OR 2.47, 95% CI 2.20-2.78, p=2.3x10(-48)) or lower than 50 nmol/L (1.92, 1.70-2.16, p=1.0x10(-26)) compared with those in the lowest quartile. INTERPRETATION Variants near genes involved in cholesterol synthesis, hydroxylation, and vitamin D transport affect vitamin D status. Genetic variation at these loci identifies individuals who have substantially raised risk of vitamin D insufficiency. FUNDING Full funding sources listed at end of paper (see Acknowledgments).

Genome-wide association analysis of metabolic traits in a birth cohort from a founder population.

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.

Genome-wide association study for early-onset and morbid adult obesity identifies three new risk loci in European populations.

We analyzed genome-wide association data from 1,380 Europeans with early-onset and morbid adult obesity and 1,416 age-matched normal-weight controls. Thirty-eight markers showing strong association were further evaluated in 14,186 European subjects. In addition to FTO and MC4R, we detected significant association of obesity with three new risk loci in NPC1 (endosomal/lysosomal Niemann-Pick C1 gene, P = 2.9 × 10−7), near MAF (encoding the transcription factor c-MAF, P = 3.8 × 10−13) and near PTER (phosphotriesterase-related gene, P = 2.1 × 10−7).

A variant near MTNR1B is associated with increased fasting plasma glucose levels and type 2 diabetes risk

In genome-wide association (GWA) data from 2,151 nondiabetic French subjects, we identified rs1387153, near MTNR1B (which encodes the melatonin receptor 2 (MT2)), as a modulator of fasting plasma glucose (FPG; P = 1.3 × 10−7). In European populations, the rs1387153 T allele is associated with increased FPG (β = 0.06 mmol/l, P = 7.6 × 10−29, N = 16,094), type 2 diabetes (T2D) risk (odds ratio (OR) = 1.15, 95% CI = 1.08–1.22, P = 6.3 × 10−5, cases N = 6,332) and risk of developing hyperglycemia or diabetes over a 9-year period (hazard ratio (HR) = 1.20, 95% CI = 1.06–1.36, P = 0.005, incident cases N = 515). RT-PCR analyses confirm the presence of MT2 transcripts in neural tissues and show MT2 expression in human pancreatic islets and beta cells. Our data suggest a possible link between circadian rhythm regulation and glucose homeostasis through the melatonin signaling pathway.

A new highly penetrant form of obesity due to deletions on chromosome 16p11.2

Obesity has become a major worldwide challenge to public health, owing to an interaction between the Western ‘obesogenic’ environment and a strong genetic contribution. Recent extensive genome-wide association studies (GWASs) have identified numerous single nucleotide polymorphisms associated with obesity, but these loci together account for only a small fraction of the known heritable component. Thus, the ‘common disease, common variant’ hypothesis is increasingly coming under challenge. Here we report a highly penetrant form of obesity, initially observed in 31 subjects who were heterozygous for deletions of at least 593 kilobases at 16p11.2 and whose ascertainment included cognitive deficits. Nineteen similar deletions were identified from GWAS data in 16,053 individuals from eight European cohorts. These deletions were absent from healthy non-obese controls and accounted for 0.7% of our morbid obesity cases (body mass index (BMI) ≥ 40 kg m-2 or BMI standard deviation score ≥ 4; P = 6.4 × 10-8, odds ratio 43.0), demonstrating the potential importance in common disease of rare variants with strong effects. This highlights a promising strategy for identifying missing heritability in obesity and other complex traits: cohorts with extreme phenotypes are likely to be enriched for rare variants, thereby improving power for their discovery. Subsequent analysis of the loci so identified may well reveal additional rare variants that further contribute to the missing heritability, as recently reported for SIM1 (ref. 3). The most productive approach may therefore be to combine the ‘power of the extreme’ in small, well-phenotyped cohorts, with targeted follow-up in case-control and population cohorts.

Infant vitamin D supplementation and allergic conditions in adulthood: Northern Finland Birth Cohort 1966

Abstract: Allergen‐induced secretion of Th2‐type cytokines and IgE production have recently been reported to be increased in mice treated with 1,25(OH)2D, the active form of vitamin D. Our objective was to investigate whether vitamin D supplementation in infancy is associated with the risk of atopy, allergic rhinitis, and asthma. The Northern Finland Birth Cohort consists of all individuals in the two most northern provinces of Finland who were due to be born in 1966. Data on vitamin D supplementation during the first year of life was obtained in 1967. Current asthma and allergic rhinitis were reported at age 31 years (n= 7,648), and atopy determined by skin‐prick test in a sub‐sample still living in northern Finland or the Helsinki area (n= 5,007). The prevalence of atopy and allergic rhinitis at age 31 years was higher in participants who had received vitamin D supplementation regularly during the first year compared to others (OR 1.46, 95%CI 1.4‐2.0, and OR 1.66, 95%CI 1.1‐1.6, respectively). A similar association was observed for asthma (OR 1.35, 95%CI 0.99‐1.8). These associations persisted after adjustment for a wide range of behavioral and social factors (adjusted: OR 1.33 for all, P= 0.01 for atopy, P= 0.001 for allergic rhinitis, and P= 0.08 for asthma). We observed an association between vitamin D supplementation in infancy and an increased risk of atopy and allergic rhinitis later in life. Further study is required to determine whether these observations reflect long‐term effects on immune regulation or differences in unmeasured determinants of vitamin D supplementation.

A Comprehensive Evaluation of Potential Lung Function Associated Genes in the SpiroMeta General Population Sample

Rationale Lung function measures are heritable traits that predict population morbidity and mortality and are essential for the diagnosis of chronic obstructive pulmonary disease (COPD). Variations in many genes have been reported to affect these traits, but attempts at replication have provided conflicting results. Recently, we undertook a meta-analysis of Genome Wide Association Study (GWAS) results for lung function measures in 20,288 individuals from the general population (the SpiroMeta consortium). Objectives To comprehensively analyse previously reported genetic associations with lung function measures, and to investigate whether single nucleotide polymorphisms (SNPs) in these genomic regions are associated with lung function in a large population sample. Methods We analysed association for SNPs tagging 130 genes and 48 intergenic regions (+/−10 kb), after conducting a systematic review of the literature in the PubMed database for genetic association studies reporting lung function associations. Results The analysis included 16,936 genotyped and imputed SNPs. No loci showed overall significant association for FEV1 or FEV1/FVC traits using a carefully defined significance threshold of 1.3×10−5. The most significant loci associated with FEV1 include SNPs tagging MACROD2 (P = 6.81×10−5), CNTN5 (P = 4.37×10−4), and TRPV4 (P = 1.58×10−3). Among ever-smokers, SERPINA1 showed the most significant association with FEV1 (P = 8.41×10−5), followed by PDE4D (P = 1.22×10−4). The strongest association with FEV1/FVC ratio was observed with ABCC1 (P = 4.38×10−4), and ESR1 (P = 5.42×10−4) among ever-smokers. Conclusions Polymorphisms spanning previously associated lung function genes did not show strong evidence for association with lung function measures in the SpiroMeta consortium population. Common SERPINA1 polymorphisms may affect FEV1 among smokers in the general population.

Two new Loci for body-weight regulation identified in a joint analysis of Genome-Wide Association Studies for Early-Onset Extreme Obesity in French and German Study Groups

Meta-analyses of population-based genome-wide association studies (GWAS) in adults have recently led to the detection of new genetic loci for obesity. Here we aimed to discover additional obesity loci in extremely obese children and adolescents. We also investigated if these results generalize by estimating the effects of these obesity loci in adults and in population-based samples including both children and adults. We jointly analysed two GWAS of 2,258 individuals and followed-up the best, according to lowest p-values, 44 single nucleotide polymorphisms (SNP) from 21 genomic regions in 3,141 individuals. After this DISCOVERY step, we explored if the findings derived from the extremely obese children and adolescents (10 SNPs from 5 genomic regions) generalized to (i) the population level and (ii) to adults by genotyping another 31,182 individuals (GENERALIZATION step). Apart from previously identified FTO, MC4R, and TMEM18, we detected two new loci for obesity: one in SDCCAG8 (serologically defined colon cancer antigen 8 gene; p = 1.85×10−8 in the DISCOVERY step) and one between TNKS (tankyrase, TRF1-interacting ankyrin-related ADP-ribose polymerase gene) and MSRA (methionine sulfoxide reductase A gene; p = 4.84×10−7), the latter finding being limited to children and adolescents as demonstrated in the GENERALIZATION step. The odds ratios for early-onset obesity were estimated at ∼1.10 per risk allele for both loci. Interestingly, the TNKS/MSRA locus has recently been found to be associated with adult waist circumference. In summary, we have completed a meta-analysis of two GWAS which both focus on extremely obese children and adolescents and replicated our findings in a large followed-up data set. We observed that genetic variants in or near FTO, MC4R, TMEM18, SDCCAG8, and TNKS/MSRA were robustly associated with early-onset obesity. We conclude that the currently known major common variants related to obesity overlap to a substantial degree between children and adults.

Early Life Factors and Blood Pressure at Age 31 Years in the 1966 Northern Finland Birth Cohort

Data on the birth weight–blood pressure relationship are inconsistent. Although an inverse association has been suggested in several large studies, interpretation is complicated by publication and other biases. Few data are available on the relationship between other early growth measures and blood pressure. We examined the shape and size of association between determinants of fetal growth, size at birth, growth in infancy, and adult systolic and diastolic blood pressure at 31 years in the prospective northern Finnish 1966 birth cohort of 5960 participants. Birth weight, birth length, gestational age, ponderal index, and birth weight relative to gestational age showed a significant inverse association with blood pressure at age 31. Rapid growth in infancy (“change-up”) was positively associated with blood pressure. Adjusted regression coefficients for birth weight indicated systolic/diastolic blood pressure lower by −1.7 (95% confidence interval [CI], −2.5, −1.0)/−0.7 (95% CI, −1.4, −0.02) mm Hg for 1 kg higher birth weight. The significant inverse association between birth weight and systolic blood pressure persisted without adjustment for adult body mass index for males. Among females, gestational age showed a stronger association with blood pressure than birth weight: gestational age higher by 7 weeks (equivalent to an average of 1 kg higher birth weight) among singletons associated with −2.9 (95% CI, −4.7, −1.1) mm Hg lower systolic blood pressure. Our results support the concept that birth weight, other birth measures, and infant growth are important determinants of blood pressure and hence cardiovascular disease risk in later life.

A polymorphism within the G6PC2 gene is associated with fasting plasma glucose levels

Several studies have shown that healthy individuals with fasting plasma glucose (FPG) levels at the high end of the normal range have an increased risk of mortality. To identify genetic determinants that contribute to interindividual variation in FPG, we tested 392,935 single-nucleotide polymorphisms (SNPs) in 654 normoglycemic participants for association with FPG, and we replicated the most strongly associated SNP (rs560887, P = 4 × 10–7) in 9353 participants. SNP rs560887 maps to intron 3 of the G6PC2 gene, which encodes glucose-6-phosphatase catalytic subunit–related protein (also known as IGRP), a protein selectively expressed in pancreatic islets. This SNP was associated with FPG (linear regression coefficient β = –0.06 millimoles per liter per A allele, combined P = 4 × 10–23) and with pancreatic β cell function (Homa-B model, combined P = 3 × 10–13) in three populations; however, it was not associated with type 2 diabetes risk. We speculate that G6PC2 regulates FPG by modulating the set point for glucose-stimulated insulin secretion in pancreatic β cells.