Looking young for one's age has been a desire since time immemorial. This desire is attributable to the belief that appearance reflects health and fecundity. Indeed, perceived age predicts survival [1] and associates with molecular markers of aging such as telomere length [2]. Understanding the underlying molecular biology of perceived age is vital for identifying new aging therapies among other purposes, but studies are lacking thus far. As a first attempt, we performed genome-wide association studies (GWASs) of perceived facial age and wrinkling estimated from digital facial images by analyzing over eight million SNPs in 2,693 elderly Dutch Europeans from the Rotterdam Study. The strongest genetic associations with perceived facial age were found for multiple SNPs in the MC1R gene (p < 1 × 10(-7)). This effect was enhanced for a compound heterozygosity marker constructed from four pre-selected functional MC1R SNPs (p = 2.69 × 10(-12)), which was replicated in 599 Dutch Europeans from the Leiden Longevity Study (p = 0.042) and in 1,173 Europeans of the TwinsUK Study (p = 3 × 10(-3)). Individuals carrying the homozygote MC1R risk haplotype looked on average up to 2 years older than non-carriers. This association was independent of age, sex, skin color, and sun damage (wrinkling, pigmented spots) and persisted through different sun-exposure levels. Hence, a role for MC1R in youthful looks independent of its known melanin synthesis function is suggested. Our study uncovers the first genetic evidence explaining why some people look older for their age and provides new leads for further investigating the biological basis of how old or young people look.

Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. By performing whole-exome sequence association analyses of hematologic quantitative traits in 15,459 community-dwelling individuals, followed by in silico replication in up to 52,024 independent samples, we identified two previously undescribed coding variants associated with lower platelet count: a common missense variant in CPS1 (rs1047891, MAF = 0.33, discovery + replication p = 6.38 × 10(-10)) and a rare synonymous variant in GFI1B (rs150813342, MAF = 0.009, discovery + replication p = 1.79 × 10(-27)). By performing CRISPR/Cas9 genome editing in hematopoietic cell lines and follow-up targeted knockdown experiments in primary human hematopoietic stem and progenitor cells, we demonstrate an alternative splicing mechanism by which the GFI1B rs150813342 variant suppresses formation of a GFI1B isoform that preferentially promotes megakaryocyte differentiation and platelet production. These results demonstrate how unbiased studies of natural variation in blood cell traits can provide insight into the regulation of human hematopoiesis.

Folate and its synthetic form folic acid function as donor of one-carbon units and have been, together with other B-vitamins, implicated in programming of epigenetic processes such as DNA methylation during early development. To what extent regulation of DNA methylation can be altered via B-vitamins later in life, and how this relates to health and disease, is not exactly known. The aim of this study was to identify effects of long-term supplementation with folic acid and vitamin B12 on genome-wide DNA methylation in elderly subjects. This project was part of a randomized, placebo-controlled trial on effects of supplemental intake of folic acid and vitamin B12 on bone fracture incidence (B-vitamins for the PRevention Of Osteoporotic Fractures (B-PROOF) study). Participants with mildly elevated homocysteine levels, aged 65-75 years, were randomly assigned to take 400 μg folic acid and 500 μg vitamin B12 per day or a placebo during an intervention period of 2 years. DNA was isolated from buffy coats, collected before and after intervention, and genome-wide DNA methylation was determined in 87 participants (n = 44 folic acid/vitamin B12, n = 43 placebo) using the Infinium HumanMethylation450 BeadChip.

Tobacco smoking, a risk factor for diabetes, is an established modifier of DNA methylation. We hypothesised that tobacco smoking modifies DNA methylation of genes previously identified for diabetes.

Decline in muscle strength with aging is an important predictor of health trajectory in the elderly. Several factors, including genetics, are proposed contributors to variability in muscle strength. To identify genetic contributors to muscle strength, a meta-analysis of genomewide association studies of handgrip was conducted. Grip strength was measured using a handheld dynamometer in 27 581 individuals of European descent over 65 years of age from 14 cohort studies. Genomewide association analysis was conducted on ~2.7 million imputed and genotyped variants (SNPs). Replication of the most significant findings was conducted using data from 6393 individuals from three cohorts. GWAS of lower body strength was also characterized in a subset of cohorts. Two genomewide significant (P-value< 5 × 10(-8) ) and 39 suggestive (P-value< 5 × 10(-5) ) associations were observed from meta-analysis of the discovery cohorts. After meta-analysis with replication cohorts, genomewide significant association was observed for rs752045 on chromosome 8 (β = 0.47, SE = 0.08, P-value = 5.20 × 10(-10) ). This SNP is mapped to an intergenic region and is located within an accessible chromatin region (DNase hypersensitivity site) in skeletal muscle myotubes differentiated from the human skeletal muscle myoblasts cell line. This locus alters a binding motif of the CCAAT/enhancer-binding protein-β (CEBPB) that is implicated in muscle repair mechanisms. GWAS of lower body strength did not yield significant results. A common genetic variant in a chromosomal region that regulates myotube differentiation and muscle repair may contribute to variability in grip strength in the elderly. Further studies are needed to uncover the mechanisms that link this genetic variant with muscle strength.

The bone-sparing effect of estrogens is mediated primarily via estrogen receptor (ER)α, which stimulates gene transcription through activation function (AF)-1 and AF-2. The role of ERαAF-1 for the estradiol (E2) effects is tissue specific. The selective ER modulators (SERMs) raloxifene (Ral), lasofoxifene (Las), and bazedoxifene (Bza) can be used to treat postmenopausal osteoporosis. They all reduce the risk for vertebral fractures, whereas Las and partly Bza, but not Ral, reduce the risk for nonvertebral fractures. Here, we have compared the tissue specificity of Ral, Las, and Bza and evaluated the role of ERαAF-1 for the effects of these SERMs, with an emphasis on bone parameters. We treated ovariectomized (OVX) wild-type (WT) mice and OVX mice lacking ERαAF-1 (ERαAF-1(0)) with E2, Ral, Las, or Bza. All three SERMs increased trabecular bone mass in the axial skeleton. In the appendicular skeleton, only Las increased the trabecular bone volume/tissue volume and trabecular number, whereas both Ral and Las increased the cortical bone thickness and strength. However, Ral also increased cortical porosity. The three SERMs had only a minor effect on uterine weight. Notably, all evaluated effects of these SERMs were absent in ovx ERαAF-1(0) mice. In conclusion, all SERMs had similar effects on axial bone mass. However, the SERMs had slightly different effects on the appendicular skeleton since only Las increased the trabecular bone mass and only Ral increased the cortical porosity. Importantly, all SERM effects require a functional ERαAF-1 in female mice. These results could lead to development of more specific treatments for osteoporosis.

Obesity is highly heritable. Genetic variants showing robust associations with obesity traits have been identified through genome-wide association studies. We investigated whether a composite score representing healthy diet modifies associations of these variants with obesity traits. Totally, 32 body mass index (BMI)- and 14 waist-hip ratio (WHR)-associated single nucleotide polymorphisms were genotyped, and genetic risk scores (GRS) were calculated in 18 cohorts of European ancestry (n = 68 317). Diet score was calculated based on self-reported intakes of whole grains, fish, fruits, vegetables, nuts/seeds (favorable) and red/processed meats, sweets, sugar-sweetened beverages and fried potatoes (unfavorable). Multivariable adjusted, linear regression within each cohort followed by inverse variance-weighted, fixed-effects meta-analysis was used to characterize: (a) associations of each GRS with BMI and BMI-adjusted WHR and (b) diet score modification of genetic associations with BMI and BMI-adjusted WHR. Nominally significant interactions (P = 0.006-0.04) were observed between the diet score and WHR-GRS (but not BMI-GRS), two WHR loci (GRB14 rs10195252; LYPLAL1 rs4846567) and two BMI loci (LRRN6C rs10968576; MTIF3 rs4771122), for the respective BMI-adjusted WHR or BMI outcomes. Although the magnitudes of these select interactions were small, our data indicated that associations between genetic predisposition and obesity traits were stronger with a healthier diet. Our findings generate interesting hypotheses; however, experimental and functional studies are needed to determine their clinical relevance.

Screening gene function in vivo is a powerful approach to discover novel drug targets. We present high-throughput screening (HTS) data for 3 762 distinct global gene knockout (KO) mouse lines with viable adult homozygous mice generated using either gene-trap or homologous recombination technologies. Bone mass was determined from DEXA scans of male and female mice at 14 weeks of age and by microCT analyses of bones from male mice at 16 weeks of age. Wild-type (WT) cagemates/littermates were examined for each gene KO. Lethality was observed in an additional 850 KO lines. Since primary HTS are susceptible to false positive findings, additional cohorts of mice from KO lines with intriguing HTS bone data were examined. Aging, ovariectomy, histomorphometry and bone strength studies were performed and possible non-skeletal phenotypes were explored. Together, these screens identified multiple genes affecting bone mass: 23 previously reported genes (Calcr, Cebpb, Crtap, Dcstamp, Dkk1, Duoxa2, Enpp1, Fgf23, Kiss1/Kiss1r, Kl (Klotho), Lrp5, Mstn, Neo1, Npr2, Ostm1, Postn, Sfrp4, Slc30a5, Slc39a13, Sost, Sumf1, Src, Wnt10b), five novel genes extensively characterized (Cldn18, Fam20c, Lrrk1, Sgpl1, Wnt16), five novel genes with preliminary characterization (Agpat2, Rassf5, Slc10a7, Slc26a7, Slc30a10) and three novel undisclosed genes coding for potential osteoporosis drug targets.

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common, complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified ten new risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with new secondary signals at four of these loci). Notably, the new loci include candidate genes with roles in the regulation of innate host defenses and T cell function, underscoring the important contribution of (auto)immune mechanisms to atopic dermatitis pathogenesis.

The third generation selective estrogen receptor modulators lasofoxifene (las) and bazedoxifene (bza) are indicated for treatment of postmenopausal osteoporosis. 17β-Estradiol (E2) and the second generation SERM raloxifene (ral) have major effects on the immune system, particularly on B cells. Treatment with E2 or ral inhibits B lymphopoiesis and treatment with E2, but not ral, stimulates antibody production. The effects of las and bza on the immune system have not been studied. Therefore, the aim of this study was to investigate their role in B cell development, maturation, and function. C57BL/6 mice were sham-operated or ovariectomized (ovx) and treated with vehicle, E2, ral, las, or bza. All substances increased total bone mineral density in ovx mice, as measured by peripheral quantitative computed tomography. In uterus, bza alone lacked agonistic effect in ovx mice and even acted as an antagonist in sham mice. As expected, E2 decreased B cell numbers at all developmental stages from pre-BI cells (in bone marrow) to transitional 1 (T1) B cells (in spleen) and increased marginal zone (MZ) B cells as determined by flow cytometry. However, treatment with las or bza only decreased the last stages of bone marrow B cell development and splenic T1 B cells, but had no effect MZ B cells. E2 increased antibody-producing cells quantified by ELISPOT, but las or bza did not. In conclusion, las and bza differ from E2 by retaining normal number of cells at most B cell stages during B lymphopoiesis and maturation and by not increasing antibody-producing cells.

Treatment with anti-inflammatory glucocorticoids is associated with osteoporosis. Many of the treated patients are postmenopausal women, who even without treatment have an increased risk of osteoporosis. Lymphocytes have been shown to play a role in postmenopausal and arthritis-induced osteoporosis, and they are targeted by glucocorticoids. The aim of this study was to investigate the mechanisms behind effects of glucocorticoids on bone during health and menopause, focusing on lymphocytes. Female C57BL/6 or SCID mice were therefore sham-operated or ovariectomized and 2 weeks later treatment with dexamethasone (dex), the nonsteroidal anti-inflammatory drug carprofen, or vehicle was started and continued for 2.5 weeks. At the termination of experiments, femurs were phenotyped using peripheral quantitative computed tomography and high-resolution micro-computed tomography, and markers of bone turnover were analyzed in serum. T and B lymphocyte populations in bone marrow and spleen were analyzed by flow cytometry. Dex-treated C57BL/6 mice had increased trabecular bone mineral density, but lower cortical content and thickness compared with vehicle-treated mice. The dex-treated mice also had lower levels of bone turnover markers and markedly decreased numbers of spleen T and B lymphocytes. In contrast, these effects could not be repeated when mice were treated with the nonsteroidal anti-inflammatory drug carprofen. In addition, dex did not increase trabecular bone in ovariectomized SCID mice lacking functional T and B lymphocytes. In contrast to most literature, the results from this study indicate that treatment with dex increased trabecular bone density, which may indicate that this effect is associated with corticosteroid-induced alterations of the lymphocyte populations.

In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations.

Many genetic loci affect circulating lipid levels, but it remains unknown whether lifestyle factors, such as physical activity, modify these genetic effects. To identify lipid loci interacting with physical activity, we performed genome-wide analyses of circulating HDL cholesterol, LDL cholesterol, and triglyceride levels in up to 120,979 individuals of European, African, Asian, Hispanic, and Brazilian ancestry, with follow-up of suggestive associations in an additional 131,012 individuals. We find four loci, in/near CLASP1, LHX1, SNTA1, and CNTNAP2, that are associated with circulating lipid levels through interaction with physical activity; higher levels of physical activity enhance the HDL cholesterol-increasing effects of the CLASP1, LHX1, and SNTA1 loci and attenuate the LDL cholesterol-increasing effect of the CNTNAP2 locus. The CLASP1, LHX1, and SNTA1 regions harbor genes linked to muscle function and lipid metabolism. Our results elucidate the role of physical activity interactions in the genetic contribution to blood lipid levels.

Cognition in adults shows variation due to developmental and degenerative components. A recent genome-wide association study identified genetic variants for general cognitive function in 148 independent loci. Here, we aimed to elucidate possible developmental and neurodegenerative pathways underlying these genetic variants by relating them to functional, clinical and neuroimaging outcomes. This study was conducted within the population-based Rotterdam Study (N=11,496, mean age 65.3±9.9 years, 58.0% female). We used lead variants for general cognitive function to construct a polygenic score (PGS), and additionally excluded developmental variants at multiple significance thresholds. A higher PGS was related to more years of education (β=0.29, p=4.3x10-7) and a larger intracranial volume (β=0.05, p=7.5x10-4). To a smaller extent, the PGS was associated with less cognitive decline (βΔG-factor=0.03, p=1.3x10-3), which became non-significant after adjusting for education (p=1.6x10-2). No associations were found with daily functioning, dementia, parkinsonism, stroke or microstructural white matter integrity. Excluding developmental variants attenuated nearly all associations. In conclusion, this study suggests that the genetic variants identified for general cognitive function are acting mainly through the developmental pathway of cognition. Therefore, cognition, assessed cross-sectionally, seems to have limited value as a biomarker for neurodegeneration.

Both fibroblast growth factors (FGFs), by binding to FGF receptors (FGFRs), and activation of the gravitostat, by artificial loading, decrease the body weight (BW). Previous studies demonstrate that both the FGF system and loading have the capacity to regulate BW independently of leptin. The aim of the current study was to determine the possible interactions between the effect of increased loading and the FGF system for the regulation of BW. We observed that the BW-reducing effect of increased loading was abolished in mice treated with a monoclonal antibody directed against FGFR1c, suggesting interactions between the two systems. As serum levels of endocrine FGF21 and hepatic FGF21 mRNA were increased in the loaded mice compared with the control mice, we first evaluated the loading response in FGF21 over expressing mice with constant high FGF21 levels. Leptin treatment, but not increased loading, decreased the BW in the FGF21-overexpressing mice, demonstrating that specifically the loading effect is attenuated in the presence of high activity in the FGF system. However, as FGF21 knockout mice displayed a normal loading response on BW, FGF21 is neither mediating nor essential for the loading response. In conclusion, the BW-reducing effect of increased loading but not of leptin treatment is blocked by high activity in the FGF system. We propose that both the gravitostat and the FGF system regulate BW independently of leptin and that pharmacologically enhanced activity in the FGF system reduces the sensitivity of the gravitostat.

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Age at adiposity rebound (AR) is associated with obesity and Type 2 Diabetes in adults. The aim of the present study was to investigate the role of age at AR in adult fat mass, fat distribution and pubertal timing for a Swedish cohort.

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes--PRDM16, PAX3, TP63, C5orf50, and COL17A1--in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

To assess candidate genes for association with osteoarthritis (OA) and identify promising genetic factors and, secondarily, to assess the candidate gene approach in OA.