doi: 10.1038/s41598-022-17443-1.
A systematic review of skin ageing genes: gene pleiotropy and genes on the chromosomal band 16q24.3 may drive skin ageing
Affiliations
- PMID:35907981
- PMCID: PMC9338925
- DOI: 10.1038/s41598-022-17443-1
Item in Clipboard
A systematic review of skin ageing genes: gene pleiotropy and genes on the chromosomal band 16q24.3 may drive skin ageing
Jun Yan Ng et al. Sci Rep..
Display options
Format
Display options
Format
Abstract
Skin ageing is the result of intrinsic genetic and extrinsic lifestyle factors. However, there is no consensus on skin ageing phenotypes and ways to quantify them. In this systematic review, we first carefully identified 56 skin ageing phenotypes from multiple literature sources and sought the best photo-numeric grading scales to evaluate them. Next, we conducted a systematic review on all 44 Genome-wide Association Studies (GWAS) on skin ageing published to date and identified genetic risk factors (2349 SNPs and 366 genes) associated with skin ageing. We identified 19 promising SNPs found to be significantly (p-Value < 1E-05) associated with skin ageing phenotypes in two or more independent studies. Here we show, using enrichment analyses strategies and gene expression data, that (1) pleiotropy is a recurring theme among skin ageing genes, (2) SNPs associated with skin ageing phenotypes are mostly located in a small handful of 44 pleiotropic and hub genes (mostly on the chromosome band 16q24.3) and 32 skin colour genes. Since numerous genes on the chromosome band 16q24.3 and skin colour genes show pleiotropy, we propose that (1) genes traditionally identified to contribute to skin colour have more than just skin pigmentation roles, and (2) further progress towards understand the development of skin pigmentation requires understanding the contributions of genes on the chromosomal band 16q24.3. We anticipate our systematic review to serve as a hub to locate primary literature sources pertaining to the genetics of skin ageing and to be a starting point for more sophisticated work examining pleiotropic genes, hub genes, and skin ageing phenotypes.
© 2022. The Author(s).
Conflict of interest statement
F.T.C. reports grants from Singapore Ministry of Education Academic Research Fund, Singapore Immunology Network, National Medical Research Council (Singapore), Biomedical Research Council (Singapore), and the Agency for Science Technology and Research (Singapore), during the conduct of the study; and has received consultancy fees from Sime Darby Technology Centre, First Resources Ltd, Genting Plantation, and Olam International, outside the submitted work. The other authors declare no other competing interests.
Figures
Workflow for screening through publications to include in the review.
(a) SNPs associated with the phenotype: wrinkles, site: whole face. This is a phenotype under Category A Phenotype—skin wrinkling and sagging-related phenotypes. Pooled Odds Ratio and the 95% Confidence Interval are calculated from both studies which reported significant associations between the SNP and the phenotype. An inconsistency index (I2 index) value ≥ 50% and a Het p-Value < 0.05 are considered statistically significant for heterogeneity. Begg’s funnel plots and Egger’s test were used to assess publication bias; Begg’s funnel plots were only drawn for meta-analysis of ≥ 2 studies or cohorts and Egger’s test was only performed for meta-analysis of ≥ 3 studies or cohorts to estimate asymmetry of data; an Egger’s test p-Value < 0.05 was considered as that publication bias possibly exist. (b) SNPs associated with the phenotype: skin type—inability to tan, site: whole body. This is a phenotype under Category B Phenotype—skin pigmentation-related phenotypes. Pooled Odds Ratio and the 95% Confidence Interval are calculated from both studies which reported significant associations between the SNP and the phenotype. An inconsistency index (I2 index) value ≥ 50% and a Het p-Value < 0.05 are considered statistically significant for heterogeneity. Begg’s funnel plots and Egger’s test were used to assess publication bias; Begg’s funnel plots were only drawn for meta-analysis of ≥ 2 studies or cohorts and Egger’s test was only performed for meta-analysis of ≥ 3 studies or cohorts to estimate asymmetry of data; an Egger’s test p-Value < 0.05 was considered as that publication bias possibly exist. (c) SNPs associated with the phenotype: skin colour, site: upper inner arm. This is a phenotype under Category B Phenotype—skin pigmentation-related phenotypes. Pooled Odds Ratio and the 95% Confidence Interval are calculated from both studies which reported significant associations between the SNP and the phenotype. An inconsistency index (I2 index) value ≥ 50% and a Het p-Value < 0.05 are considered statistically significant for heterogeneity. Begg’s funnel plots and Egger’s test were used to assess publication bias; Begg’s funnel plots were only drawn for meta-analysis of ≥ 2 studies or cohorts and Egger’s test was only performed for meta-analysis of ≥ 3 studies or cohorts to estimate asymmetry of data; an Egger’s test p-Value < 0.05 was considered as that publication bias possibly exist. (d) Forest plot for all skin ageing phenotypes associated with rs258322 under the additive allele model. Pooled Odds Ratios (OR) and 95% Confidence Intervals (CI) under the random effects model were computed for each phenotype with ≥ 2 studies or cohorts. Effect alleles were streamlined as allele G; in cases where the primary literature reports the effect allele to be allele A, reciprocals of the OR and CI reported in the primary literature were taken where appropriate to ensure compatibility. (e) Begg’s funnel plots and Egger’s test p-Values for rs258322 and skin ageing.
(a) SNPs associated with the phenotype: wrinkles, site: whole face. This is a phenotype under Category A Phenotype—skin wrinkling and sagging-related phenotypes. Pooled Odds Ratio and the 95% Confidence Interval are calculated from both studies which reported significant associations between the SNP and the phenotype. An inconsistency index (I2 index) value ≥ 50% and a Het p-Value < 0.05 are considered statistically significant for heterogeneity. Begg’s funnel plots and Egger’s test were used to assess publication bias; Begg’s funnel plots were only drawn for meta-analysis of ≥ 2 studies or cohorts and Egger’s test was only performed for meta-analysis of ≥ 3 studies or cohorts to estimate asymmetry of data; an Egger’s test p-Value < 0.05 was considered as that publication bias possibly exist. (b) SNPs associated with the phenotype: skin type—inability to tan, site: whole body. This is a phenotype under Category B Phenotype—skin pigmentation-related phenotypes. Pooled Odds Ratio and the 95% Confidence Interval are calculated from both studies which reported significant associations between the SNP and the phenotype. An inconsistency index (I2 index) value ≥ 50% and a Het p-Value < 0.05 are considered statistically significant for heterogeneity. Begg’s funnel plots and Egger’s test were used to assess publication bias; Begg’s funnel plots were only drawn for meta-analysis of ≥ 2 studies or cohorts and Egger’s test was only performed for meta-analysis of ≥ 3 studies or cohorts to estimate asymmetry of data; an Egger’s test p-Value < 0.05 was considered as that publication bias possibly exist. (c) SNPs associated with the phenotype: skin colour, site: upper inner arm. This is a phenotype under Category B Phenotype—skin pigmentation-related phenotypes. Pooled Odds Ratio and the 95% Confidence Interval are calculated from both studies which reported significant associations between the SNP and the phenotype. An inconsistency index (I2 index) value ≥ 50% and a Het p-Value < 0.05 are considered statistically significant for heterogeneity. Begg’s funnel plots and Egger’s test were used to assess publication bias; Begg’s funnel plots were only drawn for meta-analysis of ≥ 2 studies or cohorts and Egger’s test was only performed for meta-analysis of ≥ 3 studies or cohorts to estimate asymmetry of data; an Egger’s test p-Value < 0.05 was considered as that publication bias possibly exist. (d) Forest plot for all skin ageing phenotypes associated with rs258322 under the additive allele model. Pooled Odds Ratios (OR) and 95% Confidence Intervals (CI) under the random effects model were computed for each phenotype with ≥ 2 studies or cohorts. Effect alleles were streamlined as allele G; in cases where the primary literature reports the effect allele to be allele A, reciprocals of the OR and CI reported in the primary literature were taken where appropriate to ensure compatibility. (e) Begg’s funnel plots and Egger’s test p-Values for rs258322 and skin ageing.
Visual representation of the pleiotropic diversity of genes on the chromosomal band 16q24.3. Each gene shown in this figure is associated with morphologically distinct phenotypes from different phenotypic categories (e.g., Phenotype Category A, B, C, and D colour coded in yellow, blue, orange, and green respectively).
Functional enrichment analysis of skin ageing genes. Each gene shown in this figure is associated with two or more morphologically distinct phenotypes from different phenotypic categories (e.g., Phenotype Category A, B, C, and D). Nodes coloured red represent genes on the chromosomal band 16q24.3 and nodes coloured green represent skin colour genes. Edges (i.e., lines) linking the node indicate evidence of enrichment between the two connected nodes.
Network enrichment of skin ageing hub proteins. Only one strong network (i.e., Network score ≥ 4) is identified for each key skin ageing phenotype. Red gene nodes are genes on the chromosomal band 16q24.3 and yellow gene nodes are skin colour genes. Edges linking the node indicate evidence of enrichment between the two connected nodes.
Functional enrichment of skin ageing genes reveals that skin colour genes are pleiotropic for morphologically diverse (Category A and B) skin ageing phenotypes.
Functional enrichment of skin ageing genes reveals that genes on the chromosomal band 16q24.3 are pleiotropic for morphologically diverse (Category A, B, C, and D) skin ageing phenotypes.
Gene expression of the pleiotropic and hub gene SPIRE2 is downregulated in non-infant skin, aged skin, and melanoma skin, when compared against expression levels in infant skin, young skin, and non-melanoma control skin respectively. Melanoma is a skin ageing phenotype. Expression ratio is the ratio of the average expression level in the cases (e.g., aged skin tissue) against the average expression level of the reference group (e.g., young skin tissue). Ratios above 1 indicate higher expression levels in the cases group while ratios below 1 indicate lower expression levels in the cases group.
Similar articles
- Principal component analysis of seven skin-ageing features identifies three main types of skin ageing.Pardo LM, Hamer MA, Liu F, Velthuis P, Kayser M, Gunn DA, Nijsten T.Pardo LM, et al.Br J Dermatol. 2020 Jun;182(6):1379-1387. doi: 10.1111/bjd.18523. Epub 2019 Dec 23.Br J Dermatol. 2020.PMID:31519034
- GWA-based pleiotropic analysis identified potential SNPs and genes related to type 2 diabetes and obesity.Zeng Y, He H, Zhang L, Zhu W, Shen H, Yan YJ, Deng HW.Zeng Y, et al.J Hum Genet. 2021 Mar;66(3):297-306. doi: 10.1038/s10038-020-00843-4. Epub 2020 Sep 18.J Hum Genet. 2021.PMID:32948839Free PMC article.
- Additional common loci associated with stroke and obesity identified using pleiotropic analytical approach.Wang L, Xu F, Brickell A, Sun N, Mao X, Zhang Q, Wang G, Zhou Q, Yang B, Li F, Yue L, Zhang W, Hao Y, Sun C.Wang L, et al.Mol Genet Genomics. 2020 Mar;295(2):439-451. doi: 10.1007/s00438-019-01630-3. Epub 2019 Dec 7.Mol Genet Genomics. 2020.PMID:31813042
- Abundant pleiotropy in human complex diseases and traits.Sivakumaran S, Agakov F, Theodoratou E, Prendergast JG, Zgaga L, Manolio T, Rudan I, McKeigue P, Wilson JF, Campbell H.Sivakumaran S, et al.Am J Hum Genet. 2011 Nov 11;89(5):607-18. doi: 10.1016/j.ajhg.2011.10.004.Am J Hum Genet. 2011.PMID:22077970Free PMC article.Review.
- Comprehensive identification of pleiotropic loci for body fat distribution using the NHGRI-EBI Catalog of published genome-wide association studies.Kaur Y, Wang DX, Liu HY, Meyre D.Kaur Y, et al.Obes Rev. 2019 Mar;20(3):385-406. doi: 10.1111/obr.12806. Epub 2018 Nov 22.Obes Rev. 2019.PMID:30565845Review.
Cited by
- Comparisons between eyebags, droopy eyelids, and eyebrow positioning identified by photo-numeric scales or identified by written descriptive scales: Insights from the Singapore/Malaysia cross-sectional genetics epidemiology study (SMCGES) cohort.Ng JY, Chew FT.Ng JY, et al.Skin Res Technol. 2024 Feb;30(2):e13620. doi: 10.1111/srt.13620.Skin Res Technol. 2024.PMID:38376131Free PMC article.
- Epigenetics insights from perceived facial aging.Vladimir K, Perišić MM, Štorga M, Mostashari A, Khanin R.Vladimir K, et al.Clin Epigenetics. 2023 Nov 3;15(1):176. doi: 10.1186/s13148-023-01590-x.Clin Epigenetics. 2023.PMID:37924108Free PMC article.
- Genetical Signature-An Example of a Personalized Skin Aging Investigation with Possible Implementation in Clinical Practice.Sepetiene R, Patamsyte V, Valiukevicius P, Gecyte E, Skipskis V, Gecys D, Stanioniene Z, Barakauskas S.Sepetiene R, et al.J Pers Med. 2023 Aug 25;13(9):1305. doi: 10.3390/jpm13091305.J Pers Med. 2023.PMID:37763073Free PMC article.
- Comparisons between wrinkles and photo-ageing detected and self-reported by the participant or identified by trained assessors reveal insights from Chinese individuals in the Singapore/Malaysia Cross-sectional Genetics Epidemiology Study (SMCGES) cohort.Ng JY, Zhou H, Li T, Chew FT.Ng JY, et al.J Physiol Anthropol. 2024 May 18;43(1):14. doi: 10.1186/s40101-024-00361-8.J Physiol Anthropol. 2024.PMID:38762735Free PMC article.
- Comparisons between Caucasian-validated and Chinese-validated photo-numeric scales for assessing facial wrinkles.Ng JY, Zhou H, Li T, Chew FT.Ng JY, et al.Sci Rep. 2024 Nov 16;14(1):28293. doi: 10.1038/s41598-024-78945-8.Sci Rep. 2024.PMID:39550380Free PMC article.
References
Publication types
MeSH terms
Related information
LinkOut - more resources
Full Text Sources
Medical