Identification of Genes Promoting Skin Youthfulness by Genome-Wide Association Study
Total Page:16
File Type:pdf, Size:1020Kb
ORIGINAL ARTICLE Identification of Genes Promoting Skin Youthfulness by Genome-Wide Association Study Anne L.S. Chang1, Gil Atzmon2,AvivBergman2, Samantha Brugmann3, Scott X. Atwood1, Howard Y. Chang1 and Nir Barzilai2 To identify genes that promote facial skin youthfulness (SY), a genome-wide association study on an Ashkenazi Jewish discovery group (n ¼ 428) was performed using Affymetrix 6.0 Single-Nucleotide Polymorphism (SNP) Array. After SNP quality controls, 901,470 SNPs remained for analysis. The eigenstrat method showed no stratification. Cases and controls were identified by global facial skin aging severity including intrinsic and extrinsic parameters. Linear regression adjusted for age and gender, with no significant differences in smoking history, body mass index, menopausal status, or personal or family history of centenarians. Six SNPs met the À 8 À 8 Bonferroni threshold with Palleleo10 ; two of these six had Pgenotypeo10 . Quantitative trait loci mapping confirmed linkage disequilibrium. The six SNPs were interrogated by MassARRAY in a replication group (n ¼ 436) with confirmation of rs6975107, an intronic region of KCND2 (potassium voltage-gated channel, Shal-related family member 2) (Pgenotype ¼ 0.023). A second replication group (n ¼ 371) confirmed rs318125, downstream of DIAPH2 (diaphanous homolog 2 (Drosophila)) (Pallele ¼ 0.010, Pgenotype ¼ 0.002) and rs7616661, downstream of EDEM1 (ER degradation enhancer, mannosidase a-like 1) (Pgenotype ¼ 0.042). DIAPH2 hasbeenassociatedwith premature ovarian insufficiency, an aging phenotype in humans. EDEM1 associates with lifespan in animal models, although not humans. KCND2 is expressed in human skin, but has not been associated with aging. These genes represent new candidate genes to study the molecular basis of healthy skin aging. Journal of Investigative Dermatology (2014) 134, 651–657; doi:10.1038/jid.2013.381; published online 17 October 2013 INTRODUCTION Prior skin aging studies have categorized skin aging Genes promoting overall youthfulness in the form of delayed parameters into intrinsic versus extrinsic characteristics, aging have been reported in the centenarian population believed to reflect genetic versus environmental factors (Perls et al., 2002; Barzilai et al., 2003; Atzmon et al., contributing to skin aging features, respectively (e.g., Fisher 2009b; Tazearslan et al., 2011), but whether these genes et al., 2002; Vierko¨tter and Krutmann, 2012). However, also promote skin youthfulness (SY) is unclear. In fact, genes data suggest that extrinsic skin aging parameters can have a that promote SY have, to date, not been identified. Multiple strong genetic basis. For instance, MC1R gene mutations sources have estimated that over US$140 billion yearly is can promote severe photoaging (Elfakir et al., 2010). spent on cosmetics worldwide (Oh, 2006), testifying to the Heritability analyses in twins have shown that extrinsic intense and broad interest in promoting SY. Anecdotally, some skin aging parameters such ‘‘sun damage’’ exhibit 60% centenarians possess facial skin that appears decades younger heritability (Gunn et al., 2009). Recently, extrinsic skin aging than their chronologic age. This dissociation between parameters of facial photoaging have been linked to a phenotypic and chronologic age may be a unique window genetic change, a single-nucleotide polymorphism (SNP) in to examine whether particular genotypes associate with SY. the intronic area of the STXBP5L gene, in a genome- wide association study of French women (Le Clerc et al., 2013). Hence, this study aims to examine skin aging 1Department of Dermatology, Stanford University School of Medicine, in a global sense, and includes both intrinsic and extrinsic 2 Redwood City, California, USA; Albert Einstein College of Medicine, Bronx, skin aging parameters, both of which have genetic New York, USA and 3Department of Pediatrics, Cincinnati Children’s Hospital and Medical Center, University of Cincinnati, Cincinnati, Ohio, USA underpinnings. Correspondence: Anne L.S. Chang, Department of Dermatology, Stanford The primary objective of this study was to assess whether University School of Medicine, 450 Broadway Street, Mail Code 5334, there may be genes that confer SY in a unique centenarian Redwood City, California 94603, USA. E-mail: [email protected] population, by identifying individuals that lack global skin Abbreviations: GWAS, genome-wide association study; LC, Langerhans cell; aging features (both intrinsic and extrinsic as both have a SAS, skin aging score; SNP, single-nucleotide polymorphism; SY, skin genetic basis) compared with individuals of similar chrono- youthfulness logic age. The secondary objective of this study was to Received 31 January 2013; revised 10 August 2013; accepted 12 August 2013; accepted article preview online 13 September 2013; published online determine whether SY genes overlap with longevity genes in 17 October 2013 this centenarian population. & 2014 The Society for Investigative Dermatology www.jidonline.org 651 ALS Chang et al. Genes Promoting Skin Youthfulness RESULTS 2012). Identity-by-descent analysis showed that all individuals Cases and controls for facial SY phenotype were identified in the discovery group were unrelated (Supplementary Figure using a global skin aging score (SAS) incorporating intrinsic S3 online). Quantile–quantile plots of the P-values are shown and extrinsic skin aging parameters (Guinot et al., 2002). The in Supplementary Figure S4 online. À 8 average difference between chronologic age and SAS (i.e., CA Six SNPs met the Bonferroni threshold with Palleleo10 ; À 8 minus SAS), or average ‘‘delta,’’ for the entire discovery group two of these six had Pgenotypeo10 (Table 2a). A Manhattan was 15.5. To avoid bias in assigning SY cases and controls for plot of the discovery group is shown in Figure 1. Quantitative the genome-wide association study (GWAS), individuals trait loci mapping confirmed linkage disequilibrium (Figure 2). whose ‘‘deltas’’ were greater than the average ‘‘delta’’ were To explore the effects of the asymmetry in percent of defined as cases (individuals with SY), and those whose centenarians between SY cases and SY controls in the ‘‘deltas’’ were less than or equal to the group average were discovery group, we corrected for this and report the effect defined as controls (individuals without SY). There were no on the Palleles for each of the six SNPs Supplementary Table S1 individuals whose ‘‘delta’’ equaled the group average. online. These P-values remain significant when the asymmetry The global SASs were determined by a board-certified is corrected for. dermatologist experienced in assessing skin aging parameters A gene network to assess coexpression of the six significant (see Materials and Methods section). Bland–Altman analysis SNPs identified in the discovery group is shown in Supple- showed the differences in skin aging ratings at two time points mentary Figure S5 online. The six SNPs form four clusters of (Supplementary Figure S1 online). The rating differences at coexpression but no physical interaction. two time points that fell outside the 95th percentile were To confirm these findings of the discovery group, the six distributed across the range of average skin ages, indicating candidate SNPs were interrogated by MassARRAY in a lack of bias in rating for any particular skin age (Supple- replication group (n ¼ 436). Clinical characteristics of the case mentary Figure S1 online). and control groups of the replication group are reported in Clinical characteristics of the case and control groups are Supplementary Table S2a online. After covariate correction reported in Table 1. Linear regression adjusted for age, gender, for age and gender, the significance of rs6975107 was con- and personal or family history of centenarians. There were no firmed. This SNP is located in the intronic region of KCND2 significant differences in smoking history, body mass index, or (potassium voltage-gated channel, Shal-related family personal history of skin cancer between cases and controls. All member 2) (Pgenotype ¼ 0.023). female subjects in the discovery group were postmenopausal. The six candidate SNPs were interrogated by MassARRAY A conventional genome-wide association study on an in a second replication group (n ¼ 371). Clinical character- Ashkenazi Jewish discovery group (n ¼ 428) was performed istics of the case and control groups of the second replication using Affymetrix 6.0 SNP Array. SNP quality controls included group are reported in Supplementary Table S2B online. After removal of 22,118 SNPs that were not in Hardy–Weinberg covariate correction for age and gender, the significance of equilibrium, with 901,470 SNPs remaining for analysis. rs318125 and rs7616661 were confirmed (Tables 2b and c). Eigenstrat method showed no stratification, with good homo- The SNP rs318125 is downstream of DIAPH2 (diaphanous geneity between cases and controls (Supplementary Figure S2 homolog 2 (Drosophila)) (Pallele ¼ 0.010, Pgenotype ¼ 0.002). online). Isolation-by-distance studies in this Jewish diaspora The SNP rs7616661 is downstream of EDEM1 (ER degradation were reported previously (Atzmon et al., 2010; Kenny et al., enhancer, mannosidase a-like 1) (Pgenotype ¼ 0.042). Table 1. Characteristics of cases (SY group) and controls (no-SY group) in the discovery group (n ¼ 428) Cases: SY Controls: no-SY Characteristic Delta4avg. delta Deltapavg. delta P-value N 218 210 Chronologic age, mean in years±SD 79.3±15.4 93.3±10.4 o0.0001 Delta SY, mean±SD 24.6±6.6