- 1Beiersdorf AG, Research and Development, Hamburg, Germany
- 2Institute for Bioinformatics, University Medicine Greifswald, Greifswald, Germany
- 3Division of Epigenetics, DKFZ-ZMBH Alliance, German Cancer Research Center, Heidelberg, Germany
- 4Institute for Community Medicine, University Medicine Greifswald, Greifswald, Germany
Changes in DNA methylation patterning have been reported to be a key hallmark of aged human skin. The altered DNA methylation patterns are correlated with deregulated gene expression and impaired tissue functionality, leading to the well-known skin aging phenotype. Searching for small molecules, which correct the aged methylation pattern therefore represents a novel and attractive strategy for the identification of anti-aging compounds. DNMT1 maintains epigenetic information by copying methylation patterns from the parental (methylated) strand to the newly synthesized strand after DNA replication. We hypothesized that a modest inhibition of this process promotes the restoration of the ground-state epigenetic pattern, thereby inducing rejuvenating effects. In this study, we screened a library of 1800 natural substances and 640 FDA-approved drugs and identified the well-known antioxidant and anti-inflammatory molecule dihydromyricetin (DHM) as an inhibitor of the DNA methyltransferase DNMT1. DHM is the active ingredient of several plants with medicinal use and showed robust inhibition of DNMT1 in biochemical assays. We also analyzed the effect of DHM in cultivated keratinocytes by array-based methylation profiling and observed a moderate, but significant global hypomethylation effect upon treatment. To further characterize DHM-induced methylation changes, we used published DNA methylation clocks and newly established age predictors to demonstrate that the DHM-induced methylation change is associated with a reduction in the biological age of the cells. Further studies also revealed re-activation of age-dependently hypermethylated and silenced genes in vivo and a reduction in age-dependent epidermal thinning in a 3-dimensional skin model. Our findings thus establish DHM as an epigenetic inhibitor with rejuvenating effects for aged human skin.
Introduction
Aging is a time-dependent general loss of functionality and fitness across lifespan, and the parallel accumulation of molecular damages, including altered epigenomic regulation is considered one of the key hallmarks of aging (Booth and Brunet, 2016). DNA methylation is a key epigenetic modification that has been associated with many aspects of human health and disease (Bergman and Cedar, 2013). DNA methylation is catalyzed by a conserved class of enzymes that are known as DNA methyltransferases (DNMTs) and that cooperate to establish and maintain DNA methylation patterns during embryonic development and tissue homeostasis (Lyko, 2018). Broadly, DNA methylation patterns are established during early embryogenesis by DNMT3A and DNMT3B, and then maintained by DNMT1. DNMT1 is an essential enzyme that functions to copy methylation patterns from the parental (methylated) strand to the newly synthesized strand after DNA replication (Jurkowska and Jeltsch, 2016).
The combined activity of DNMTs results in widespread methylation of the human genome. However, there are also regions that are characterized by low or absent levels of methylation. For example, CpG islands, which are CpG-dense regions that are often associated with the promoters of housekeeping genes, are usually devoid of methylation (Jones, 2012). Furthermore, gene regulatory elements, such as transcription factor binding sites, often show dynamic and context-dependent methylation levels. Indeed, DNA methylation has been shown to modulate transcription factor binding (Yin et al., 2017), which provides an explanation for how DNA methylation modulates gene expression (Schübeler, 2015).
While the role of DNA methylation in mammalian cell fate-specification is now widely accepted (Smith and Meissner, 2013), its role in environmental adaption has been more difficult to substantiate. In this context, the analysis of human skin can be considered particularly relevant, as skin represents the most environmentally exposed human organ and its phenotype is intricately linked to epigenetic regulatory mechanisms (Köhler and Rodríguez-Paredes, 2020). For example, we have shown that skin aging is associated with defined DNA methylation changes, especially hypermethylation (Gronniger et al., 2010; Bormann et al., 2016). A detailed comparison of single-base resolution epidermis methylomes also revealed age-related DNA methylation changes in gene regulatory regions (Raddatz et al., 2013). This provided a mechanistic explanation for how altered DNA methylation in old skin could affect skin phenotypes. Finally, we also showed that age-related methylation changes in the human epidermis can be used to predict the chronological age of sample donors and used these findings to establish an epidermis DNA methylation clock (Bormann et al., 2016). Together, these studies comprehensively illustrate the role of altered DNA methylation patterns in human skin aging.
It has been suggested that healthy lifestyle, which is based on the avoidance of environmental stressors such as UV light and smoking, but also intake of dietary antioxidants, may improve skin appearance through epigenetic modulation (Alegría-Torres et al., 2011). Improved epigenetic plasticity through a moderate inhibition of the methylation pattern could plausibly support this process. As DNA methylation needs to be maintained through an active mechanism, it can be modified by inhibiting the corresponding enzymes, most notably DNMT1. We hypothesize that a moderate reduction of DNMT1 activity might result in a reduction of age-dependent hypermethylation and promote the restoration of the epigenetic pattern to the ground-state. This provides an attractive strategy for identifying new anti-aging compounds.
We have now screened a library of natural substances to identify active compounds that inhibit DNMT1 in a biochemical assay. This identified dihydromyricetin (DHM), a compound from plants that are popular in traditional medicines that is already known for its beneficial anti-cancer, antioxidant, and anti-inflammatory properties (Kang et al., 2011; Zhang et al., 2011). Further analyses showed that this nutritional flavonoid inhibited DNA methylation in cultured primary human keratinocytes and induced a reduction in DNA methylation age, which is an epigenetic predictor of the skin aging phenotype. Additionally, the beneficial modulation of the methylation pattern by DHM is reflected by a reactivation of the expression of marker genes for age-related epigenetic gene silencing in vivo. Finally, our results show that DHM has rejuvenating effects in human skin models, which indicates considerable potential for cosmetical applications.
Results
Identification of DHM as a natural small-molecule DNMT inhibitor
We used full-length recombinant human DNMT1 to establish a homogenous 384-well scintillation proximity assay (SPA) with a biotinylated 42mer hemi-methylated oligonucleotide that was captured with yttrium silicate (YSi) streptavidin coated SPA beads. This assay showed clear dose-response curves for S-adenosyl-homocysteine, a known DNMT1 product inhibitor, and doxorubicin, a DNA intercalating compound (Figure 1A). We therefore used this assay to screen a library of 1800 natural products and 640 FDA-approved drugs. Two replicate screens were carried out and showed high reproducibility (r = 0.78, Figure 1B). Screening hits contained a number of substances that are known DNA intercalators, such as doxorubicin, as well as known natural compound DNMT inhibitors, such as epigallocatechinol-3-gallate (Fang et al., 2003), and dietary polyphenols like quercetin, which is known for its anti-cancer effect (Kedhari Sundaram et al., 2019). Interestingly, the screen also identified myricetin (Figure 1C) with a mean IC50 of 43.37 µM (Supplementary Figure S1A), a common flavonoid found in grapes, berries, fruits, vegetables, and herbs, among the most potent hits. We then obtained three highly (>95%) pure extracts containing the colorless myricetin derivative DHM (3,3′,4′,5,5’,7-hexahydroxy-2,3-dihydroflavanonol, Figure 1C). DHM can be extracted from plants that are being used in traditional medicines (Zhang et al., 2018), e.g., the Vine tea Ampelopsis grossedentata (DHM2 and DHM3) and the bark, leaves and flowers of the Wax Myrtle tree (Myrica cerifera, DHM1). Testing of the three extracts in an independent biochemical DNMT assay (see Materials and Methods for details) indicated significant (Padj < 0.01, one-way ANOVA with Bonferroni correction) DNMT1 inhibitory activity for DHM in the similar concentration range like myricetin (Figure 1D; Supplementary Figure S1B) and no significant differences in efficacy between the three tested extracts (one-way ANOVA). Due to its availability, favorable material properties and known tolerability in humans (Zhang et al., 2018), we prioritized DHM for further analyses. Taken together, these finding provide first evidence for DHM being a DNMT1 inhibitor.
FIGURE 1. Identification of dihydromyricetin (DHM) as a DNMT inhibitor. (A) Establishment of a biochemical assay with human recombinant DNMT1. The plot shows dose-response curves for S-adenosyl-homocysteine (SAH) and doxorubicin (Doxo). (B) Scatter plot showing results from two independent screens. (C) Chemical structures of myricetin (Myr) and dihydromyricetin (DHM). (D) Testing of three independent DHM extracts in an independent biochemical DNMT assay (n = 10). *Padj < 0.01 (one-way ANOVA with Bonferroni correction).
DHM causes DNA hypomethylation in primary human keratinocytes
To test the ability of DHM to modulate DNA methylation in human cells, we cultivated 13 independent primary human keratinocyte lines in the presence or absence of the highest tolerable concentration of DHM (Supplementary Figure S2). Subsequently, we isolated genomic DNA from the cultured cells and analyzed genomic DNA methylation patterns using Infinium 850 k arrays. Further data analysis that considered the heterogeneity between individual cell lines identified 101,067 differentially methylated probes (see Materials and Methods for details) that were used for a more detailed analysis. Principal Component Analysis based on the differentially methylated probes showed a detectable separation of DHM-treated and control (DMSO) cells (Figure 2A) which is in agreement with the overall accuracy obtained with random-forest based classification (Supplementary Figure S3), indicating directional effects across all individual cell lines. Further analysis of the differentially methylated probes in all primary human keratinocytes from 13 different donors also showed a moderate, but robust and highly significant (p = 1.4 × 10−33, t-test) reduction of DNA methylation in the DHM-treated cells (Figure 2B). Interestingly, this effect was evident and statistically significant in all epigenomic features, except CpG islands and the surrounding CpG island shores, which represent the features with the lowest baseline methylation levels in the human genome (Figure 2C). Together, these results demonstrate that DHM causes moderate changes of DNA methylation in primary keratinocytes and provide important proof-of mechanism for DHM inhuman cells.
FIGURE 2. DHM inhibits DNA methylation in primary human keratinocytes. (A) Principal component analysis based on 101,067 differentially methylated probes in the complete dataset obtained with primary human keratinocytes from different donors (n = 13). (B) DHM-dependent hypomethylation in primary human keratinocytes. Box plots indicate the distribution of beta values of the 101,067 differentially methylated probes in the complete dataset. *p = 1.4 × 10−33, t-test. (C) DHM-dependent hypomethylation in epigenomic substructures. *Padj < 6.7 × 10−4, t-test.
DHM treatment does not affect general markers of epigenetic and genetic integrity in primary human keratinocytes
Profound losses of DNA methylation have been shown to affect epigenetic and genetic integrity, resulting in apoptosis (Yokochi and Robertson, 2004; Laranjeira et al., 2023). We therefore analyzed whether DHM induced global loss of methylation in keratinocytes, using LINE-1 retroelement methylation as an established marker (Zhang et al., 2020). The results showed no reduction of LINE-1 methylation in DHM-treated primary keratinocytes (Figure 3A and Supplementary Figure S4). In addition, we analyzed the development of double-strand breaks after DHM treatment by immunofluorescence staining of gamma-H2AX, an established marker of DNA damage. Our results showed no DHM-dependent increase in gamma-H2AX staining (Figure 3B) in cultured primary keratinocytes from six independent donors. Finally, we investigated whether DHM activates Caspase 3/7 in keratinocytes, which is a widely used marker for apoptosis induction. The results showed no increase in Caspase 3/7 activity after treatment (Figure 3C). Together, these findings suggest that DHM does not affect general markers of epigenetic and genetic integrity in keratinocytes.
FIGURE 3. No evidence for critical epigenetic changes upon DHM treatment in keratinocytes. (A). Direct pyrosequencing after bisulfite conversion of the LINE-1 sequence covering methylation of 4 representative CpGs. Bar plot represents mean methylation difference of the analyzed LINE-1 sequence between 20 µM DHM and DMSO treatment in % (n = 4). (B) DNA damage assay. Bar plots depicting number of DNA double strand breaks per cell in primary human keratinocytes upon DHM treatment (blue) and positive control H2O2 (yellow) compared to the DMSO control (gray). *p < 0.05 (n = 6; 2-way ANOVA). (C) Apoptosis assay. Bar plots depicting caspase 3/7 activity in primary human keratinocytes upon DHM treatment (blue) and positive control staurosporine (orange) compared to the DMSO control (gray). *p < 0.01 (n = 10; one-way ANOVA, Dunnett’s multiple comparison test).
DHM affects epigenetic predictors of the skin phenotype
In subsequent analyses, we investigated the relationship between DNA methylation and skin aging phenotypes. In initial analyses, we used a previously established epidermis DNA methylation clock (Bormann et al., 2016) to calculate the epigenetic age of 461 female volunteers from a cohort study (see Material and Methods for details). For the same group of volunteers, we then quantitatively determined their skin aging phenotypes by assessing the visual skin appearance (visual age) and by measuring skin elasticity, respectively. Partial correlation analysis identified a significant (p < 0.05, Pearson correlation) positive relationship between epigenetic age and visual age after controlling the influence of the chronological age on both variables (Figure 4A). Similarly, we found a significant (p < 0.05, Pearson correlation) inverse association between methylation age and skin elasticity adjusting for chronological age (Figure 4B). These results identify a correlation between epigenetic age and skin aging phenotypes independent of chronological age and establish a suitable framework for evaluating the effect of DHM on epigenetic predictors of the skin phenotype.
FIGURE 4. DHM reduces DNA methylation age and predicted wrinkle grade, epigenetic predictors of the skin aging phenotype. Association between DNA methylation age and (A) visual age and (B) skin elasticity, adjusting for the influence of the chronological age using a skin specific clock (Bormann et al., 2016). Pearson correlation coefficient (r) and significance (p) are depicted. DHM reduces the DNA methylation age of cultured primary human keratinocytes determined (C) by the Skin & Blood clock (Horvath et al., 2018) and (D) the newly established human epidermis clock. (E) Induced reduction of the predicted wrinkle grade upon DHM treatment determined with the recently established wrinkle predictor (see accompanying manuscript by Bienkowska et al.). Every symbol represents an individual cell line from an independent donor. *p < 0.05 (Wilcoxon rank-sum test).
Subsequent analyses then addressed the effects of DHM on DNA methylation clocks. To this end, we used 1) the established Skin & Blood Clock (Horvath et al., 2018), 2) a human epidermis clock, which we generated based on data from the 461 female volunteers (Supplementary Figure S5) and 3) a wrinkle predictor, which we recently established by training a model that predicts the wrinkle grade based on DNA methylation profiles of epidermal samples (see accompanying manuscript by Bienkowska et al.). These models were used to analyze methylation patterns from DHM-treated primary human keratinocytes, which reflect the main cell type of the epidermis. Correlating the donor age to the predicted biological age of the cell populations revealed a clear trend, which justifies applying the ex vivo developed age clocks to cell culture samples, particularly to identify relative changes after treatment (Supplementary Figure S6). The results showed a consistent and statistically significant (p < 0.05, Wilcoxon rank-sum test) reduction in the DNA biological age of around 2 years with both the Skin & Blood clock and with the human epidermis clock (Figures 4C, D). A similar effect was also observed for the predicted wrinkle grade (Figure 4E). Based on the observed correlation between chronological age and observed wrinkle grade (Supplementary Figure S7) this finding indicates a cell rejuvenation by 3.7 years of wrinkle accumulation in female volunteers. Taken together, the applied epigenetic predictors indicate a consistent DHM-dependent cell rejuvenation in the range of 2–3.7 years in keratinocytes.
In vivo and functional effects of DHM on skin aging phenotypes
Before we investigated the effect of DHM in vivo, an evaluation of the toxicological profile of DHM was performed to ensure the safe use of the intended cosmetic product application. From our studies, we observed no discernible cytotoxicity at in-use concentrations in toxicological in vitro assays with DHM, such as the phototoxicity test (Supplementary Figure S8), which we performed equivalent to the OECD Guideline 498. The safety assessment (see Material and Methods for details) concluded that DHM has a low potential for acute and repeated dose toxicity. Additionally, limits on the concentration of DHM in the formulation have been applied to ensure consumer safety.
To assess the bioavailability of DHM in human skin, we verified that the molecular properties fulfill the criteria that serve as an indicator for a potential skin barrier penetration (Lee et al., 1994; Bos and Meinardi, 2000; Cross et al., 2003; Hadgraft, 2004; Chen et al., 2018): With a molecular weight of 320.25 g/mol and a logP of 1.23 (ChemSpider) dermal absorption into the living skin cells is expected.
To characterize the effect of DHM in vivo, 30 female volunteers were recruited for a vehicle-controlled study where DHM was topically applied to the skin on their lower backs for 2 weeks before receiving a repetitive low-dose UV exposure as an aging stimulus. Subsequently, epidermis samples were obtained from the treated area and their methylation profiles were analyzed. When these profiles were intersected with epidermis methylation profiles from young and old donors (n = 204, see Materials and Methods for details), this identified a set of 12,444 age-related DNA methylation probes that were specifically differentially methylated by the DHM treatment. Notably, when compared to vehicle controls, DHM changed the methylation values towards the pattern that is representative of the young epidermis for more than 95% of these probes (Figure 5A). These findings provide an important confirmation for our observations in cultured keratinocytes and suggest an in vivo rejuvenating effect on the epigenetic profile of the skin.
FIGURE 5. In vivo and functional effects of DHM on skin aging phenotypes. (A) DHM induces the DNA methylation pattern of the young epidermis in vivo. Volcano plot of hypomethylated CpGs upon topically DHM treatment overlapping with hypomethylated CpGs in young (pink) and old (shaded pink) epidermis, respectively (left). Overlap of hypermethylated CpGs after DHM treatment with hypermethylated CpGs in young and old epidermis, respectively (right). (B) Topical DHM treatment induces the reactivation of selected marker genes (blue bars) that are subject to age-related epigenetic silencing in vivo (shaded pink bars). *p < 0.05 (Wilcoxon rank-sum test). # Significant negative correlation (Padj < 0.05, Pearson correlation with Holm correction) of gene expression with in vivo wrinkle grade and ## being among the top 10% of genes negatively correlated with in vivo wrinkle grade (see Material and Methods, and Supplementary Figure S9 for details). (C) DHM increases epidermis thickness in an organoid skin aging model. Microscopic images of 3D skin models cultivated with primary fibroblasts either form old donors (top left) or from young donors (top middle) and the resulting epidermal thickness measurements (top right). The effect of DHM was analyzed in the same experimental setup and is shown in the lower panels. Vertical turquoise bars: epidermal thickness measurements. Horizontal black bars: 50 μm *p < 0.005 (Wilcoxon rank-sum test). (D) DHM increases the proliferation of fibroblasts incubated with DHM-conditioned medium obtained from keratinocytes. Schematic representation of the crosstalk experiment is depicted on the left panel. The effect of DHM on number of counted fibroblasts is shown in the right panel (n = 6). *p < 0.05 (Wilcoxon’s signed rank test).
To further analyze the effect of DHM on in vivo gene expression, 19 female volunteers (50–65 years) were recruited to topically treat their inner forearm with a DHM containing formulation (or with the corresponding vehicle) for 8 weeks. Subsequently, epidermis samples were obtained by suction blistering, RNA was isolated and the expression level of 34 marker genes for age-related epigenetic silencing, of which 23 were associated with wrinkle development in vivo, were determined by qPCR (see Materials and Methods and Supplementary Figure S9 for details). The results showed an increased mRNA expression level for all 34 genes in DHM-treated epidermis when compared to vehicle controls (Figure 5B), indicating a substantial capacity for DHM to ameliorate age-related epigenetic gene silencing in the human epidermis.
Finally, to functionally investigate the effect of DHM on skin aging, we established a 3-dimensional skin aging model. To this end, we compared skin models that were grown with primary keratinocytes from a single donor (42 years) and primary fibroblasts of donors from two different age groups (18–33 years and 58–67 years). After cultivating the skin equivalents for 12 weeks, comparison of epidermal thickness showed a significantly (p < 0.01, Wilcoxon rank-sum test) thinner epidermis in skin models with fibroblasts from old donors (Figure 5C), which corresponds to a well-known skin aging phenotype in vivo (Giangreco et al., 2010). Treatment of old skin models (fibroblasts: 63–68 years, keratinocytes: 23–59 years) with DHM resulted in a significantly (p < 0.05, Wilcoxon rank-sum test) increased epidermal thickness (Figure 5C). To further investigate the functional effect of DHM on dermal fibroblasts, we used an experimental setup to analyze epithelial/mesenchymal interactions (Harrison et al., 2005). The results showed that DHM-conditioned supernatant from keratinocytes significantly (p < 0.05, Wilcoxon rank-sum test) stimulated fibroblast proliferation (Figure 5D). As decreased proliferation and a reduced cell number of fibroblasts is a marker for skin aging (Varani et al., 2000) these results further support a rejuvenating effect of DHM in cell culture.
Taken together, these in vivo and in vitro studies provide evidence for a skin rejuvenating effect of DHM on a molecular and functional level.
Discussion
As natural compounds are well known for their beneficial effects on human physiology, this class of substances represents an attractive resource for the screening of new anti-aging active ingredients. And because altered epigenomic regulation is considered a key hallmark of aging (Booth and Brunet, 2016), the identification of compounds that modulate epigenetic mechanisms is considered a particularly impactful approach to slow down or reverse the process. We have focused our small-molecule screen on natural products and FDA-approved drugs, two classes of compounds with established safety profiles, that do not require further animal testing and tested them for their inhibitory effect on the key epigenetic modulator DNMT1. This approach identified DHM, a flavonoid with known beneficial effects like anti-cancer, anti-oxidative and anti-inflammatory properties. DHM has also been discussed in the context of anti-aging therapy, as it has been shown that DHM protects against cognitive impairment and ameliorates brain aging in rats, as it modulates apoptosis and dysfunctional autophagy of hippocampal neurons, and astrogliosis (Kou et al., 2016; Qian et al., 2021). Additionally, DHM has been shown to activate FOXO and AOP to modulate longevity in Drosophila (Fan et al., 2020). Nevertheless, a deep understanding of the mechanisms of action is still missing. Our study could serve as a starting point to investigate how DHM could counteract or prevent aging on a cellular level, even though, our study focused on a specific part of the body, the skin. Beyond the scope of our study, DHM has also been proposed by others to serve as a promising compound to treat human aging conditions (Martínez-Coria et al., 2019). However, there are only few studies that demonstrate the health-promoting effects of DHM in a placebo-controlled manner in humans. As an example, Chen et al. observed after 12 weeks of consuming DHM a decrease in serum TNF-α (Chen et al., 2015), a blood factor which increases with aging (Paolisso et al., 1998; Bruunsgaard et al., 1999). Investigating the potential of DHM to treat aging or age-related diseases systemically would be an exciting approach for future evaluation.
DHM is considered the active ingredient of Vine tea, a health-promoting herbal tea with increasing global popularity and consumption. As DHM and its derivatives are considered to make up 1/5th of dry Vine tea leaves depending on extraction method (Ye et al., 2015; Carneiro et al., 2021), the amount of DHM in a prepared cup of Vine tea can reach mmol concentrations, thus illustrating its excellent tolerability. For topical application of DHM a safety assessment following next-generation risk assessment approaches for new cosmetic ingredients concluded that the substance has only a low potential for acute and repeated dose toxicity, which supports a safe usage up to 0.15% on human skin (see Materials and Methods for details and Supplementary Figure S8). This distinguishes DHM from the most widely used DNA demethylating agents 5-azacytidine and 2′-deoxy-5-azacytidine, which have been developed as clinical anti-cancer drugs targeting aberrant DNA methylation (Jones et al., 2016). However, the broad and indirect mode of action of these hypomethylating agents is associated with noticeable toxicity (Stresemann and Lyko, 2008). More selective small-molecule inhibitors, such as RG-108 (Brueckner et al., 2005) and GSK3685032 (Pappalardi et al., 2021), should provide important alternatives for therapeutic applications. And although a strong reduction in DNA methylation was shown to cause genetic instability and apoptosis (Jackson-Grusby et al., 2001; Eden et al., 2003; Chen et al., 2007), partial inhibition by the highly selective small-molecule DNMT1 inhibitor GSK3685032 showed good in vivo tolerability (Pappalardi et al., 2021). In agreement with this notion, mice that carry a heterozygous null allele for Dnmt1 are viable and fertile, without any detectable phenotypic aberrations (Li et al., 1992). While additional work will be required to precisely define the therapeutic window for nutraceutical or cosmetical applications of DNA methylation inhibitors, exposure to DHM resulted in substantially more moderate changes in DNA methylation patterning compared to some of the established DNA methylation inhibitors. These results indicate a considerable epigenetic safety of DHM. This, and the good tolerability of DHM represent distinctive features that support the use of the compound in nutraceutical and cosmetical applications.
Our results also show that DHM reduces the DNA methylation age in cultured primary human keratinocytes. We confirmed the DHM-dependent reduction of DNA methylation age with two independent clocks, with the newly trained human epidermis clock showing a slightly stronger effect on the DNA methylation age than the well-established Skin & Blood clock (Horvath et al., 2018). This difference may be because we trained our algorithm only on epidermis samples to achieve improved prediction accuracy. Importantly, this finding adds DHM to a growing list of active substances that have a measurable impact on the biological age of human cells and tissues (Fahy et al., 2019). Furthermore, our analysis revealed an association between DNA methylation age and two key phenotypic features of human skin aging, visual appearance, and skin elasticity, independent of the chronological age. These findings are consistent with the notion that a reduction in the epigenetic age is indicative of an improved functionality of human tissues and organs (Horvath and Raj, 2018). This notion is further supported by the results obtained with our recently developed wrinkle predictor, which links the DNA methylation pattern to the arguably most relevant skin aging phenotype (see accompanying manuscript by Bienkowska et al.).
Interestingly, the moderate DNMT1 inhibition by DHM was associated with a measurable rejuvenating effect on the level of DNA methylation and gene expression in vivo. The vast majority of differentially methylated CpGs revealed a methylation pattern that was similar to the young epidermis. Although DHM was originally identified as a DNMT1 inhibitor and subsequently shown to induce moderate global hypomethylation in vitro, hypo-as well as hypermethylated CpGs were observed after the 2-weeks of topical DHM treatment in vivo. This can be explained by secondary epigenetic effects leading to re-methylation of distinct CpGs and/or through a reset of the genetically encoded methylation pattern. For example, our data indicate that the expression of TET2 was upregulated in keratinocytes after DHM treatment (Supplementary Figure S10). TET2 is an epigenetic enzyme which antagonizes DNMTs to shape methylation patterns (Lyko, 2018). In addition to TET2, the expression of two other epigenetic modifiers, KDM5C and PRDM2, was also upregulated upon DHM treatment in keratinocytes (Supplementary Figure S10). While KDM5C demethylates histones H3K4 and H3K9, PRDM2 methylates histone H3K4. The upregulation of both enzymes may lead to downstream changes in H3K4 and H3K9 methylation patterns, which can influence the binding of DNMTs, and consequently DNA methylation (Li et al., 2021). Further research will be required to determine the mechanisms of DHM-associated CpG hypermethylation in detail.
In agreement with the rejuvenating effects observed at the level of DNA methylation, topical DHM treatment also had a measurable effect on age-related silencing of gene expression in vivo. This effect was observed for genes from various pathways, indicating a broad impact across various skin aging-related processes and several of the analyzed marker genes were related to the skin aging phenotype of wrinkle formation. For example, LYPD3 acts on cell-cell adhesion and has been described to play a role in epidermis differentiation and stratification (Jacobsen et al., 2017). In addition, our studies showed an association between reduced expression of LYPD3 with increased wrinkle development in human skin. Also, HAS3 expression is related to wrinkle formation and beside HAS1 it is suggested to play a crucial role in the regulation of hyaluronic acid synthesis in the epidermis (Sayo et al., 2002), thus being important for keratinocyte differentiation, as hyaluronic acid is involved in the formation, organization, and stabilization of the epidermal extracellular matrix (Malaisse et al., 2014). SPON2 is modulating cell communication between keratinocytes and fibroblasts, as it has been shown to be secreted by keratinocytes and induces fibrogenesis (Evrard et al., 2021). IRS2 is important for adequate wound healing (Lima et al., 2012; Manzano-Nunez et al., 2019), which is delayed in elderly skin (Amini-Nik et al., 2022) and loss of EDA is associated with dry skin and eczema (Kowalczyk-Quintas et al., 2015), which are more frequently observed in elderly people (White-Chu and Reddy, 2011). Loss of SLC25A25 is associated with a reduced basal mitochondrial respiration and ATP content in mice (Anunciado-Koza et al., 2011), which represents another key hallmark of human skin aging (Wei et al., 2009; Schniertshauer et al., 2018) and our studies indicate an association with increased wrinkle development in human skin.
Restoring the regulatory landscape of an aged cell could thus have a more comprehensive impact on the multidimensional process of aging than targeting single factors or pathways. This hypothesis was further supported by our studies using an organotypic skin culture. 3-dimensional skin models capture the complex process of skin aging better than 2-dimensional cell culture models. The robust dermal compartment, created by fibroblasts grown within an organized extracellular matrix, supports the development of a differentiated, stratified, and keratinized epidermis, with spatiotemporal cellular communication and barrier properties (Costello et al., 2022). By integrating cell populations from older donors intrinsic aging is simulated and skin equivalents exhibited many aspects of the ageing phenotype including epidermal atrophy and reduced organization of the epidermal strata (Ahlers et al., 2022). Hence, the 3D model is suitable to investigate holistic anti-aging effects, specifically for the human skin. Interestingly, using such a 3-dimensional aging model, we observed a DHM-dependent transformation of the epidermal structure towards a more juvenile, which further establishes DHM as an interesting natural compound that links epigenetic modulation to skin rejuvenation. Future studies will be required to investigate whether DHM also improves healthy skin aging in vivo, either as a standalone application or in combination with healthy lifestyles.
Materials and Methods
Ethic approval
Ethical approval was obtained in consideration of the Declaration of Helsinki and the guideline of the International Conference on Harmonization Good Clinical Practice (ICH GCP) by the International Medical & Dental Ethics Commission in Hamburg (Std. no. 67686) and by the Independent Ethics Committee Freiburg (feki code 08/2610).
Production of DNMT1 protein
Full-length human DNMT1 protein was produced and purified as described before (Castellano et al., 2008). Briefly, the protein was expressed in insect cells and purified by affinity chromatography and gel filtration. The protein concentration of purified DNMT was determined by Bradford assay and verified by using Coomassie stained SDS/polyacrylamide gels and suitable molecular mass markers of known concentration.
Compound screening
A homogenous 384-well scintillation proximity assay (SPA) was developed and optimized for screening of DNMT-1 inhibitors. In this assay 3H-SAM (PerkinElmer) is used as a methyl-group donor in the methylation of a biotinylated 42mer hemi-methylated substrate comprising the following two oligonucleotides (IBA): 5′ (bio-LC)GAT CCG ACG ACG ACG ACG AXG ACG ACG ACG ACG ACG ACG ATC (X = 5 methylcytosine) and 5′GAT CGT CGT CGT CGT CGT CGT CGT CGT CGT CGT CGT CGG ATC. The methylated is captured with yttrium silicate (YSi) streptavidin coated SPA beads (PerkinElmer). The library containing the subsets: FDA approved drugs (Enzo Life Sciences), Naturstoffe (Enzo Life Sciences), Natural products (Cfm Tropitzsch), Flavonoid derivatives (TimTec) and a collection of additional compounds (Beiersdorf AG) was screened on two separate experimental days against recombinant human DNMT1 at the DKFZ/EMBL Chemical Biology Core Facility. Compounds/substances were diluted 1:25 in H2O and dispensed into 384-well Optiplates. 800 nM 3H-SAM was added to 100 nM 42mer oligo and 30 nM DNMT1 to start the methylation reaction (180 min at room temperature). Subsequently, YSi beads were added with 372 mM NaCl to stop the methylation reaction and the plates were read on a MicroBeta LumiJet.
Biochemical DNMT assay
Myriceline™ SPE90 extract from Myrica cerifera (>95% DHM, Provital, “DHM1”) and two dihydromyricetin extracts from Ampelopsis grossedentata (>95% DHM, Naturalin “DHM2”, BOC Sciences “DHM3”) were dissolved in DMSO. DNMT inhibition was tested using the DNMT Activity Quantification Kit (Abcam) following the manufacturer’s instructions. 0.1% DMSO was used for controls.
DNA methylation analysis of primary human keratinocytes
Primary human keratinocyte lines from 13 independent donors were cultivated for 3 days in KGM Gold media (Lonza) supplemented with 20 µM DHM (Provital) in DMSO (final concentration: 0.1%). Genomic DNA was isolated with the QIAamp DNA Investigator Kit (Qiagen) and processed for Methylation EPIC arrays (Illumina). Methylation data analysis was carried out using the R Bioconductor package Minfi (Aryee et al., 2014) and limma (Ritchie et al., 2015). Specifically, raw. idat files were read and preprocessed. Methylation loci (probes) were filtered for high detection p-value [p > 0.01, as provided by Minfi (Aryee et al., 2014)], location on sex chromosomes, ability to self-hybridize, and potential SNP contamination. Array normalization was performed using the preprocessFunnorm function, available in Minfi (Aryee et al., 2014). Quality control was performed after every preprocessing step. Subsequently, differentially methylated probes were identified by fitting a linear model with a paired sample design followed by statistical analysis using an empirical Bayes method to moderate standard errors. Lastly, differentially methylated probes were filtered by significance threshold (Padj < 0.05, t-test, after correction for multiple testing using the Benjamini–Hochberg method).
Additionally, we trained a classifier to discriminate between DHM-treated and DMSO control cells using a machine learning approach. This analysis was performed using the randomForest function provided by the randomForest R package (Liaw and Wiener, 2002). Training was performed using the differentially methylated probes from the previous analysis, with the treatment status of the samples converted into a factor variable. The random forest classifier was built with 500 trees, and 10-fold cross-validation was applied.
Primary human keratinocyte lines from 2 independent donors and 2 keratinocyte cell lines were cultivated for 3 days in KGM Gold media (Lonza) supplemented with 20 µM DHM (Provital) in DMSO (final concentration: 0.1%). Genomic DNA was isolated with the QIAamp DNA Investigator Kit (Qiagen) and LINE-1 promotor methylation (sequence to analyze: TTYGTGGTGYGTYGTTTTTTAAGTYGGTTTGAAAAG) was analyzed by EpigenDX using direct pyrosequencing after bisulfite conversion.
Training of DNA methylation clock and methylation analysis
We made use of a methylation dataset (Methylation EPIC arrays, Illumina) obtained from epidermis samples of the population-based Study of Health in Pomerania (SHIP-TREND-1) (Volzke et al., 2022) comprising 461 female participants. The data were processed as described previously (see section analysis of primary human keratinocytes) using quantile normalization (preprocessQuantile). Before training, one probe (cg26614073) was excluded from the dataset to avoid conflict with Horvath’s patent claims (EP 3 049 535 B1). After splitting 378 females of the dataset (80% for training and 20% for testing), the training of the age clock was conducted with the beta-values of the training dataset applying the cv. glmnet function (alpha = 0, lambda = 1408.78) of the glmnet R package (Friedman et al., 2010) in a 10-fold cross-validation mode. The trained clock was validated with the test set and an independent publicly available dataset (Holzscheck et al., 2020a) (Supplementary Figure S5A) and outperformed the Skin & Blood clock established by (Horvath et al., 2018) (Supplementary Figure S5B). The methylation levels of the independent dataset were obtained and pre-processed in the same way as the dataset used for the training. Finally, the obtained DNA methylation clock was used to determine the methylation age of the primary cultured keratinocytes. Additionally, the methylation age of the primary cultured keratinocytes was determined using the published Skin & Blood clock (Horvath et al., 2018). Finally, the predicted wrinkle grade of the primary cultured keratinocytes was determined using our recently developed wrinkle predictor (see accompanying manuscript by Bienkowska et al.).
Correlation of skin aging phenotype and DNA methylation age
We used the methylation and phenotypic data from epidermal samples of SHIP-TREND-1 (Volzke et al., 2022). We calculated the DNA methylation age via a published skin specific DNA methylation clock (Bormann et al., 2016) to make use of the whole dataset, as the published clock was trained on data from another study. For the same study participant their visual age was determined as the average rated age appearance based on the volunteer’s portrait pictures and assessed by more than 30 experts. Additionally, before epidermal samples were taken, the skin elasticity on the sample area was measured for each study participant using a cutometer® (Courage & Khazaka Electronic GmbH, Cologne, Germany) as previously described (Dobrev, 2002). The elasticity was defined as the average of 3 consecutive measurements of the R7 value. Partial correlation analysis was performed as follows: First, the residuals of the linear regression between the parameters of interest and the chronological age were calculated (y ∼ chronological age; with y = methylation age, visual age, and skin elasticity, respectively). Subsequently, the association between the residuals was analyzed by Pearson correlation analysis to determine the significance of the correlation.
In vivo DNA methylation analysis
The skin of the lower back from 30 female volunteers was treated for 2 weeks by topical application with a formulation based on oil-in-water emulsions (o/w) containing 0.15% DHM (w/w) or vehicle, followed by repetitive UV exposure at a sub-erythemal dose (Holzscheck et al., 2020b) as an aging stimulus. 24 h after the last UV exposure, suction blister samples were obtained (Südel et al., 2003), DNA was isolated and DNA methylation patterns were analyzed using Infinium Methylation EPIC arrays (Illumina). Data were processed using Minfi (Aryee et al., 2014) and quantile normalization was applied. For the subsequent analyses M values were used to describe CpG methylation levels. CpGs were considered significantly differentially methylated, if the adjusted p-value was <0.05 using the Benjamini–Hochberg procedure. Data analysis revealed 16,389 CpGs, which were specifically differentially methylated in a DHM-dependent manner on the irradiated skin area. In parallel, we used the methylation data from the big cohort study (see section training of DNA methylation clock and methylation age analysis) analyzing the methylation patterns from 102 younger (28–40 years) and 102 older volunteers (65–88 years) and we found 270,642 probes to be differentially methylated in aged samples. The intersection of the 16,389 DHM-dependent probes with the 270,642 age-dependent probes defined the set of 12,444 DHM-dependent and age-dependent probes.
In vivo gene expression analysis
The skin of the inner forearm from 19 female volunteers (50–65 years) was treated for 8 weeks by topical application with a formulation based on oil-in-water emulsions (o/w) containing 0.15% DHM (w/w) or vehicle. After completion of the treatment, suction blisters were obtained for RNA isolation. In parallel, we analyzed gene expression datasets (RNA-Seq) obtained from epidermis samples of SHIP-TREND-1 (Volzke et al., 2022), which were analyzed with the R Bioconductor package DESeq2 (Love et al., 2014). This identified 2,820 genes which were downregulated (logFC <0, Padj ≤ 0.05) in the epidermis of older volunteers (≥70 years) compared to the younger group (≤40 years) and which were associated with hypermethylated CpGs. The expression of 34 of these genes (selected based on keratinocyte-relevant pathway association) was analyzed using qRT-PCR. To this end, total RNA of suction blister roofs was isolated using the RNeasy Fibrous Tissue Mini Kit (Qiagen) according to the manufacturer’s protocol. After reverse transcription using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems) gene expression was detected using the Real-Time TaqMan PCR (Applied Biosystems). Ct values were normalized to endogenous housekeeping gene 18sRNA and gene expression fold changes were calculated via the 2−ΔΔCt method.
3D skin model and histology
Three-dimensional skin models were generated as previously described (Boehnke et al., 2007). Hyalograft-3D was replaced by the non-woven fibrous scaffold M3-II Bemcot (Asahi Kasei, Tokyo, Japan). After 6 weeks of culture, DHM at a final concentration of 20 µM dissolved in DMSO (0.1%) was added to the medium for additional 6 weeks. Sections of 3D skin equivalents were stained using AEMA fast staining Kit (LT Lehmann, Berlin, Germany) according to manufacturer’s protocol and epidermal thickness was determined using Fiji software (Schindelin et al., 2012).
Crosstalk experiment
Fibroblasts were seeded in 96-well plates at a density of 5 × 103 cells/well in DMEM (Gibco) with 10% FCS (Biowest) and 1% Penicillin/Streptomycin (Sigma-Aldrich) and keratinocytes were seeded in six-well-plates at a density of 5 × 103 cells/well in KGM Gold (Lonza). Keratinocytes were treated with 39 µM DHM for 24 h. Keratinocyte-conditioned media were filtered (Sartorius Vivaspin Turbo4 10,000 MW Cutoff) and retentates were mixed with fresh KGM and transferred to fibroblasts. After 24 h fibroblasts were PFA-fixed, stained (EarlyTox DeadGreen, Molecular Devices) and imaged (SpectraMax, Molecular Devices).
DNA damage assay
Primary human keratinocytes were seeded in 96-well plates at a density of 10 × 103 cells/well in KGM Gold (Lonza) and preincubated overnight at 37°C/5% CO2. Cells were then treated with 10 μM and 20 µM DHM for 72 h at 37°C/5% CO2. DNA double strand breaks were induced in positive control wells with 250 µM hydrogen peroxide (AppliChem) for 15 min. Cells were then fixed with formaldehyde followed by immunostaining with Anti-Phospho-H2A.X-AlexaFluor(™)555 antibody (Merck Millipore) and counterstaining with Hoechst 33,342 (ThermoFisher). Image acquisition and analysis was performed on the semi-automated high-content imaging system scanR (Olympus).
Apoptosis assay
Primary human keratinocytes were seeded in 96-well plates at a density of 10 × 103 cells/well in KGM Gold (Lonza) and preincubated overnight at 37°C/5% CO2. Cells were treated with 10 μM and 20 µM DHM for 72 h at 37°C/5% CO2. Caspase 3/7 activity was assessed using the Caspase 3/7-Glo luminescence assay (Promega).
Safety assessment
Safety assessment followed next-generation risk assessment approaches for new cosmetic ingredients without animal testing. These combine any existing and relevant data from in vivo animal studies with data generated with in vitro assays, in silico predictions, read-across analyses and other biological or computational models (Alexander-White et al., 2022; Westmoreland et al., 2022) to cover toxicological endpoints such as acute toxicity, skin irritation, skin sensitization, genotoxicity and systemic toxicity.
Phototoxicity assay
To determine the phototoxicity of UV absorbing compounds the Phototoxicity Test according to the OECD Guideline 498 was performed using the human reconstructed epidermis model EpiDerm™ (MatTek; InVitro Life Science Laboratories, s. r.o.) according to the previously described protocol (Kandarova and Liebsch, 2017). Briefly, DHM was applied topically to epidermal surface, 24 h later the skin equivalents were exposed to UV irradiation (Hönle SOL 500; applied dosis: 1.7 mW/cm2 für 60 min = 6 J/cm2). After additional 24 h cytotoxicity has been evaluated by reduction of mitochondrial conversion of MTT to formazan (MatTek; InVitro Life Science Laboratories, s. r.o.).
Cell viability assay
To determine the most tolerable in vitro concentration of DHM, we measured endogenous esterase activity. Briefly, keratinocytes of three different donors were cultured in 96-well plates at a density of 10 × 103 cells/well in KGM Gold (Lonza) for 72 h in the presence of 0.125µM–160 µM DHM. Subsequently, cells were washed with 1× Dulbecco’s Phosphate Buffered Saline (DPBS) and incubated for 15 min at 37°C in 100μL/well fluorescein diacetate (FDA) (Sigma, Taufkirchen, Germany) solution (15 μg/mL FDA in 1 × DPBS). Fluorescence was then determined at 517 nm in the 96-well plate reader infinite M200 (Tecan, Crailsheim, Germany).
In vitro gene expression analysis
To determine if DHM stimulates gene expression of epigenetic regulators in skin cells, human primary epidermal keratinocytes were isolated from skin biopsies obtained from six different donors (55–72 years) as described previously (Reichert et al., 2017). After isolation, keratinocytes were seeded in six-well plates (7.5 × 104 cells/well) and cultured in KGM-Gold (Lonza) for 24 h. In the following, cells were treated with 20 µM DHM for 72 h. Total RNA was then extracted using the RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. After hybridization using target-specific Reporter and Capture probes, reactions were purified and imaged on the nCounter system (NanoString). Raw mRNA counts were normalized using five endogenous housekeeping genes (GAPDH, GUSB, OAZ1, PUM1, UBC) and six NanoString positive controls that were spiked into every sample. Fold Change values were calculated as the ration of normalized mRNA counts (DHM vs. 0.1% DMSO treatment).
Correlation of wrinkle grade and gene expression
To further demonstrate the relevance of the selected age-related epigenetically silenced genes for the skin aging phenotype, we analyzed the relation between their expression and the grade of in vivo observed wrinkle formation. Like for the correlation analysis of the skin aging phenotype with the DNA methylation age, we made use of the data from the SHIP-TREND-1 study (Volzke et al., 2022). To ensure generalizability of the results, we used four different approaches to represent the gene expression for the correlation analysis: i) raw read counts, ii) transcripts per million (TPMs), iii) normalized read counts and iv) normalized TPMs. Normalization was performed in analogy to RT-PCR approaches by dividing the gene expression of each gene by the geometric mean of the housekeeping genes (GAPDH, GUSB, OAZ1, PUM1 and UBC). The Pearson correlation between gene expression and wrinkle grade was determined with the cor. test function in R and the calculated p-values were then adjusted for multiple testing with the R function p. adjust. Finally, only the expression of genes which had an adjusted p-value <0.05 in all four approaches were considered as significantly correlating with the wrinkle grade.
Data availability statement
The data of the cohort study are under special EU privacy legislation due to their sensitive nature and hence, only available from the SHIP consortium upon official request. The datasets analyzed for this study can be found in the Array Express database https://www.ebi.ac.uk/biostudies/arrayexpress via E-MTAB-8993.
Ethics statement
The studies involving humans were approved by International Medical & Dental Ethics Commission in Hamburg (Std. no. 67686) and by the Independent Ethics Committee Freiburg (feki code 08/2610). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
CF: Formal Analysis, Supervision, Visualization, Writing–original draft, Writing–review and editing. AB: Formal Analysis, Visualization, Writing–review and editing. JS: Investigation, Writing–review and editing. KW: Investigation, Writing–review and editing. GR: Formal Analysis, Visualization, Writing–review and editing. BK: Data curation, Formal Analysis, Visualization, Writing–review and editing. DK: Investigation, Writing–review and editing. RS: Investigation, Writing–review and editing. RK: Investigation, Writing–review and editing. JK: Investigation, Writing–review and editing. SB: Investigation, Writing–review and editing. DL: Data curation, Supervision, Writing–review and editing. AS: Resources, Writing–review and editing. HV: Resources, Writing–review and editing. LK: Supervision, Writing–review and editing. MW: Resources, Writing–review and editing. FL: Resources, Writing–original draft, Writing–review and editing. EG: Conceptualization, Project administration, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. SHIP is part of the Community Medicine Research Network for the University Medicine Greifswald and is supported by the German Federal State of Mecklenburg–West Pomerania.
Acknowledgments
We like to thank Ahmed Kamal for his support with the supply of the compounds. Additionally, we are grateful to the SHIP consortium for granting us the access to the phenotypic, DNA methylation and gene expression data of the epidermis samples. Furthermore, we thank the DKFZ Chemical Biology Core Facility for their support in the compound screening. And we also thank Steffen Erler for supporting the safety evaluation of the ingredient.
Conflict of interest
Authors CF, AB, JS, KW, BK, RS, RK, JK, SB, DL, AS, MW, and EG were employed by the company Beiersdorf AG. LK and FL received consultation fees from Beiersdorf AG.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fragi.2023.1258184/full#supplementary-material
References
Ahlers, J. M. D., Falckenhayn, C., Holzscheck, N., Solé-Boldo, L., Schütz, S., Wenck, H., et al. (2022). Single-cell RNA profiling of human skin reveals age-related loss of dermal sheath cells and their contribution to a juvenile phenotype. Front. Genet. 12, 797747. doi:10.3389/fgene.2021.797747
Alegría-Torres, J. A., Baccarelli, A., and Bollati, V. (2011). Epigenetics and lifestyle. Epigenomics 3, 267–277. doi:10.2217/epi.11.22
Alexander-White, C., Bury, D., Cronin, M., Dent, M., Hack, E., Hewitt, N. J., et al. (2022). A 10-step framework for use of read-across (RAX) in next generation risk assessment (NGRA) for cosmetics safety assessment. Regul. Toxicol. Pharmacol. 129, 105094. doi:10.1016/j.yrtph.2021.105094
Amini-Nik, S., Abdullahi, A., Vinaik, R., Yao, R. J. R., Yu, N., Datu, A., et al. (2022). Aging impairs the cellular interplay between myeloid cells and mesenchymal cells during skin healing in mice. Aging Dis. 13, 540–551. doi:10.14336/AD.2021.1008
Anunciado-Koza, R. P., Zhang, J., Ukropec, J., Bajpeyi, S., Koza, R. A., Rogers, R. C., et al. (2011). Inactivation of the mitochondrial carrier SLC25A25 (ATP-Mg2+/Pi transporter) reduces physical endurance and metabolic efficiency in mice. J. Biol. Chem. 286, 11659–11671. doi:10.1074/jbc.M110.203000
Aryee, M. J., Jaffe, A. E., Corrada-Bravo, H., Ladd-Acosta, C., Feinberg, A. P., Hansen, K. D., et al. (2014). Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 30, 1363–1369. doi:10.1093/bioinformatics/btu049
Bergman, Y., and Cedar, H. (2013). DNA methylation dynamics in health and disease. Nat. Struct. Mol. Biol. 20, 274–281. doi:10.1038/nsmb.2518
Boehnke, K., Mirancea, N., Pavesio, A., Fusenig, N. E., Boukamp, P., and Stark, H. J. (2007). Effects of fibroblasts and microenvironment on epidermal regeneration and tissue function in long-term skin equivalents. Eur. J. Cell Biol. 86, 731–746. doi:10.1016/j.ejcb.2006.12.005
Booth, L. N., and Brunet, A. (2016). The aging epigenome. Mol. Cell 62, 728–744. doi:10.1016/j.molcel.2016.05.013
Bormann, F., Rodriguez-Paredes, M., Hagemann, S., Manchanda, H., Kristof, B., Gutekunst, J., et al. (2016). Reduced DNA methylation patterning and transcriptional connectivity define human skin aging. Aging Cell 15, 563–571. doi:10.1111/acel.12470
Bos, J. D., and Meinardi, M. M. (2000). The 500 Dalton rule for the skin penetration of chemical compounds and drugs. Exp. Dermatol 9, 165–169. doi:10.1034/j.1600-0625.2000.009003165.x
Brueckner, B., Garcia Boy, R., Siedlecki, P., Musch, T., Kliem, H. C., Zielenkiewicz, P., et al. (2005). Epigenetic reactivation of tumor suppressor genes by a novel small-molecule inhibitor of human DNA methyltransferases. Cancer Res. 65, 6305–6311. doi:10.1158/0008-5472.CAN-04-2957
Bruunsgaard, H., Andersen-Ranberg, K., Jeune, B., Pedersen, A. N., Skinhøj, P., and Pedersen, B. K. (1999). A high plasma concentration of TNF-alpha is associated with dementia in centenarians. J. Gerontol. A Biol. Sci. Med. Sci. 54, M357–M364. doi:10.1093/gerona/54.7.m357
Carneiro, R. C. V., Ye, L., Baek, N., Teixeira, G. H. A., and O'Keefe, S. F. (2021). Vine tea (Ampelopsis grossedentata): a review of chemical composition, functional properties, and potential food applications. J. Funct. Foods 76, 104317. doi:10.1016/j.jff.2020.104317
Castellano, S., Kuck, D., Sala, M., Novellino, E., Lyko, F., and Sbardella, G. (2008). Constrained analogues of procaine as novel small molecule inhibitors of DNA methyltransferase-1. J. Med. Chem. 51, 2321–2325. doi:10.1021/jm7015705
Chen, C. P., Chen, C. C., Huang, C. W., and Chang, Y. C. (2018). Evaluating molecular properties involved in transport of small molecules in stratum corneum: a quantitative structure-activity relationship for skin permeability. Molecules 23, 911. doi:10.3390/molecules23040911
Chen, S., Zhao, X., Wan, J., Ran, L., Qin, Y., Wang, X., et al. (2015). Dihydromyricetin improves glucose and lipid metabolism and exerts anti-inflammatory effects in nonalcoholic fatty liver disease: a randomized controlled trial. Pharmacol. Res. 99, 74–81. doi:10.1016/j.phrs.2015.05.009
Chen, T., Hevi, S., Gay, F., Tsujimoto, N., He, T., Zhang, B., et al. (2007). Complete inactivation of DNMT1 leads to mitotic catastrophe in human cancer cells. Nat. Genet. 39, 391–396. doi:10.1038/ng1982
Costello, L., Dicolandrea, T., Tasseff, R., Isfort, R., Bascom, C., Von Zglinicki, T., et al. (2022). Tissue engineering strategies to bioengineer the ageing skin phenotype in vitro. Aging Cell 21, e13550. doi:10.1111/acel.13550
Cross, S. E., Magnusson, B. M., Winckle, G., Anissimov, Y., and Roberts, M. S. (2003). Determination of the effect of lipophilicity on the in vitro permeability and tissue reservoir characteristics of topically applied solutes in human skin layers. J. Invest. Dermatol 120, 759–764. doi:10.1046/j.1523-1747.2003.12131.x
Dobrev, H. P. (2002). A study of human skin mechanical properties by means of Cutometer. Folia Med. Plovdiv. 44, 5–10.
Eden, A., Gaudet, F., Waghmare, A., and Jaenisch, R. (2003). Chromosomal instability and tumors promoted by DNA hypomethylation. Science 300, 455. doi:10.1126/science.1083557
Evrard, C., Lambert De Rouvroit, C., and Poumay, Y. (2021). Epidermal hyaluronan in barrier alteration-related disease. Cells 10, 3096. doi:10.3390/cells10113096
Fahy, G. M., Brooke, R. T., Watson, J. P., Good, Z., Vasanawala, S. S., Maecker, H., et al. (2019). Reversal of epigenetic aging and immunosenescent trends in humans. Aging Cell 18, e13028. doi:10.1111/acel.13028
Fang, M. Z., Wang, Y., Ai, N., Hou, Z., Sun, Y., Lu, H., et al. (2003). Tea polyphenol (-)-epigallocatechin-3-gallate inhibits DNA methyltransferase and reactivates methylation-silenced genes in cancer cell lines. Cancer Res. 63, 7563–7570.
Fan, X., Zeng, Y., Fan, Z., Cui, L., Song, W., Wu, Q., et al. (2020). Dihydromyricetin promotes longevity and activates the transcription factors FOXO and AOP in Drosophila. Aging (Albany NY) 13, 460–476. doi:10.18632/aging.202156
Friedman, J. H., Hastie, T., and Tibshirani, R. (2010). Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22. doi:10.18637/jss.v033.i01
Giangreco, A., Goldie, S. J., Failla, V., Saintigny, G., and Watt, F. M. (2010). Human skin aging is associated with reduced expression of the stem cell markers beta1 integrin and MCSP. J. Invest. Dermatol 130, 604–608. doi:10.1038/jid.2009.297
Gronniger, E., Weber, B., Heil, O., Peters, N., Stab, F., Wenck, H., et al. (2010). Aging and chronic sun exposure cause distinct epigenetic changes in human skin. PLoS Genet. 6, e1000971. doi:10.1371/journal.pgen.1000971
Hadgraft, J. (2004). Skin deep. Eur. J. Pharm. Biopharm. 58, 291–299. doi:10.1016/j.ejpb.2004.03.002
Harrison, C. A., Dalley, A. J., and Mac Neil, S. (2005). A simple in vitro model for investigating epithelial/mesenchymal interactions: keratinocyte inhibition of fibroblast proliferation and fibronectin synthesis. Wound Repair Regen. 13, 543–550. doi:10.1111/j.1524-475X.2005.00076.x
Holzscheck, N., Söhle, J., Kristof, B., Grönniger, E., Gallinat, S., Wenck, H., et al. (2020a). Multi-omics network analysis reveals distinct stages in the human aging progression in epidermal tissue. AGING 12, 12393–12409. doi:10.18632/aging.103499
Holzscheck, N., Sohle, J., Schlager, T., Falckenhayn, C., Gronniger, E., Kolbe, L., et al. (2020b). Concomitant DNA methylation and transcriptome signatures define epidermal responses to acute solar UV radiation. Sci. Rep. 10, 12918. doi:10.1038/s41598-020-69683-8
Horvath, S., Oshima, J., Martin, G. M., Lu, A. T., Quach, A., Cohen, H., et al. (2018). Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria Syndrome and ex vivo studies. AGING 10, 1758–1775. doi:10.18632/aging.101508
Horvath, S., and Raj, K. (2018). DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat. Rev. Genet. 19, 371–384. doi:10.1038/s41576-018-0004-3
Jackson-Grusby, L., Beard, C., Possemato, R., Tudor, M., Fambrough, D., Csankovszki, G., et al. (2001). Loss of genomic methylation causes p53-dependent apoptosis and epigenetic deregulation. Nat. Genet. 27, 31–39. doi:10.1038/83730
Jacobsen, B., Larouche, D., Rochette-Drouin, O., Ploug, M., and Germain, L. (2017). Expression of C4.4A in an in vitro human tissue-engineered skin model. Biomed. Res. Int. 2017, 2403072. doi:10.1155/2017/2403072
Jones, P. A. (2012). Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 13, 484–492. doi:10.1038/nrg3230
Jones, P. A., Issa, J. P., and Baylin, S. (2016). Targeting the cancer epigenome for therapy. Nat. Rev. Genet. 17, 630–641. doi:10.1038/nrg.2016.93
Jurkowska, R. Z., and Jeltsch, A. (2016). Enzymology of mammalian DNA methyltransferases. Adv. Exp. Med. Biol. 945, 87–122. doi:10.1007/978-3-319-43624-1_5
Kandarova, H., and Liebsch, M. (2017). “The EpiDerm™ phototoxicity test (EpiDerm™ H3D-PT),” in Alternatives for dermal toxicity testing. Editors C. Eskes, E. Van Vliet, and H. I. Maibach (Cham: Springer International Publishing).
Kang, N. J., Jung, S. K., Lee, K. W., and Lee, H. J. (2011). Myricetin is a potent chemopreventive phytochemical in skin carcinogenesis. Ann. N. Y. Acad. Sci. 1229, 124–132. doi:10.1111/j.1749-6632.2011.06122.x
Kedhari Sundaram, M., Hussain, A., Haque, S., Raina, R., and Afroze, N. (2019). Quercetin modifies 5'CpG promoter methylation and reactivates various tumor suppressor genes by modulating epigenetic marks in human cervical cancer cells. J. Cell Biochem. 120, 18357–18369. doi:10.1002/jcb.29147
Köhler, F., and Rodríguez-Paredes, M. (2020). DNA methylation in epidermal differentiation, aging, and cancer. J. Invest. Dermatol 140, 38–47. doi:10.1016/j.jid.2019.05.011
Kou, X., Liu, X., Chen, X., Li, J., Yang, X., Fan, J., et al. (2016). Ampelopsin attenuates brain aging of D-gal-induced rats through miR-34a-mediated SIRT1/mTOR signal pathway. Oncotarget 7, 74484–74495. doi:10.18632/oncotarget.12811
Kowalczyk-Quintas, C., Schuepbach-Mallepell, S., Willen, L., Smith, T. K., Huttner, K., Kirby, N., et al. (2015). Pharmacological stimulation of Edar signaling in the adult enhances sebaceous gland size and function. J. Invest. Dermatol 135, 359–368. doi:10.1038/jid.2014.382
Laranjeira, A. B. A., Hollingshead, M. G., Nguyen, D., Kinders, R. J., Doroshow, J. H., and Yang, S. X. (2023). DNA damage, demethylation and anticancer activity of DNA methyltransferase (DNMT) inhibitors. Sci. Rep. 13, 5964. doi:10.1038/s41598-023-32509-4
Lee, C. K., Uchida, T., Kitagawa, K., Yagi, A., Kim, N.-S., and Goto, S. (1994). Skin permeability of various drugs with different lipophilicity. J. Pharm. Sci. 83, 562–565. doi:10.1002/jps.2600830424
Li, E., Bestor, T. H., and Jaenisch, R. (1992). Targeted mutation of the DNA methyltransferase gene results in embryonic lethality. Cell 69, 915–926. doi:10.1016/0092-8674(92)90611-f
Lima, M. H., Caricilli, A. M., De Abreu, L. L., Araújo, E. P., Pelegrinelli, F. F., Thirone, A. C., et al. (2012). Topical insulin accelerates wound healing in diabetes by enhancing the AKT and ERK pathways: a double-blind placebo-controlled clinical trial. PLoS One 7, e36974. doi:10.1371/journal.pone.0036974
Li, Y., Chen, X., and Lu, C. (2021). The interplay between DNA and histone methylation: molecular mechanisms and disease implications. EMBO Rep. 22, e51803. doi:10.15252/embr.202051803
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi:10.1186/s13059-014-0550-8
Lyko, F. (2018). The DNA methyltransferase family: a versatile toolkit for epigenetic regulation. Nat. Rev. Genet. 19, 81–92. doi:10.1038/nrg.2017.80
Malaisse, J., Bourguignon, V., De Vuyst, E., Lambert De Rouvroit, C., Nikkels, A. F., Flamion, B., et al. (2014). Hyaluronan metabolism in human keratinocytes and atopic dermatitis skin is driven by a balance of hyaluronan synthases 1 and 3. J. Invest. Dermatol 134, 2174–2182. doi:10.1038/jid.2014.147
Manzano-Nunez, F., Arambul-Anthony, M. J., Galan Albinana, A., Leal Tassias, A., Acosta Umanzor, C., Borreda Gasco, I., et al. (2019). Insulin resistance disrupts epithelial repair and niche-progenitor Fgf signaling during chronic liver injury. PLoS Biol. 17, e2006972. doi:10.1371/journal.pbio.2006972
Martínez-Coria, H., Mendoza-Rojas, M. X., Arrieta-Cruz, I., and López-Valdés, H. E. (2019). Preclinical research of dihydromyricetin for brain aging and neurodegenerative diseases. Front. Pharmacol. 10, 1334. doi:10.3389/fphar.2019.01334
Paolisso, G., Rizzo, M. R., Mazziotti, G., Tagliamonte, M. R., Gambardella, A., Rotondi, M., et al. (1998). Advancing age and insulin resistance: role of plasma tumor necrosis factor-alpha. Am. J. Physiol. 275, E294–E299. doi:10.1152/ajpendo.1998.275.2.E294
Pappalardi, M. B., Keenan, K., Cockerill, M., Kellner, W. A., Stowell, A., Sherk, C., et al. (2021). Discovery of a first-in-class reversible DNMT1-selective inhibitor with improved tolerability and efficacy in acute myeloid leukemia. Nat. cancer 2, 1002–1017. doi:10.1038/s43018-021-00249-x
Qian, J., Wang, X., Cao, J., Zhang, W., Lu, C., and Chen, X. (2021). Dihydromyricetin attenuates D-galactose-induced brain aging of mice via inhibiting oxidative stress and neuroinflammation. Neurosci. Lett. 756, 135963. doi:10.1016/j.neulet.2021.135963
Raddatz, G., Hagemann, S., Aran, D., Söhle, J., Kulkarni, P. P., Kaderali, L., et al. (2013). Aging is associated with highly defined epigenetic changes in the human epidermis. Epigenetics Chromatin 6, 36. doi:10.1186/1756-8935-6-36
Reichert, O., Fleming, T., Neufang, G., Schmelz, M., Genth, H., Kaever, V., et al. (2017). Impaired glyoxalase activity is associated with reduced expression of neurotrophic factors and pro-inflammatory processes in diabetic skin cells. Exp. Dermatol 26, 44–50. doi:10.1111/exd.13118
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, e47. doi:10.1093/nar/gkv007
Sayo, T., Sugiyama, Y., Takahashi, Y., Ozawa, N., Sakai, S., Inoue, S., et al. (2002). Hyaluronan synthase 3 regulates hyaluronan synthesis in cultured human keratinocytes. J. Investigative Dermatology 118, 43–48. doi:10.1046/j.0022-202x.2001.01613.x
Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682. doi:10.1038/nmeth.2019
Schniertshauer, D., Gebhard, D., and Bergemann, J. (2018). Age-dependent loss of mitochondrial function in epithelial tissue can Be reversed by coenzyme Q(10). J. Aging Res. 2018, 6354680. doi:10.1155/2018/6354680
Schübeler, D. (2015). Function and information content of DNA methylation. Nature 517, 321–326. doi:10.1038/nature14192
Smith, Z. D., and Meissner, A. (2013). DNA methylation: roles in mammalian development. Nat. Rev. Genet. 14, 204–220. doi:10.1038/nrg3354
Stresemann, C., and Lyko, F. (2008). Modes of action of the DNA methyltransferase inhibitors azacytidine and decitabine. Int. J. Cancer 123, 8–13. doi:10.1002/ijc.23607
Südel, K. M., Venzke, K., Knussmann-Hartig, E., Moll, I., Stäb, F., Wenck, H., et al. (2003). Tight control of matrix metalloproteinase-1 activity in human skin. Photochem Photobiol. 78, 355–360. doi:10.1562/0031-8655(2003)078<0355:tcomma>2.0.co;2
Varani, J., Warner, R. L., Gharaee-Kermani, M., Phan, S. H., Kang, S., Chung, J. H., et al. (2000). Vitamin A antagonizes decreased cell growth and elevated collagen-degrading matrix metalloproteinases and stimulates collagen accumulation in naturally aged human skin. J. Invest. Dermatol 114, 480–486. doi:10.1046/j.1523-1747.2000.00902.x
Volzke, H., Schossow, J., Schmidt, C. O., Jurgens, C., Richter, A., Werner, A., et al. (2022). Cohort profile update: the study of health in Pomerania (SHIP). Int. J. Epidemiol. 51, e372–e383. doi:10.1093/ije/dyac034
Wei, Y. H., Wu, S. B., Ma, Y. S., and Lee, H. C. (2009). Respiratory function decline and DNA mutation in mitochondria, oxidative stress and altered gene expression during aging. Chang. Gung Med. J. 32, 113–132.
Westmoreland, C., Bender, H. J., Doe, J. E., Jacobs, M. N., Kass, G. E. N., Madia, F., et al. (2022). Use of new approach methodologies (NAMs) in regulatory decisions for chemical safety: report from an EPAA deep dive workshop. Regul. Toxicol. Pharmacol. 135, 105261. doi:10.1016/j.yrtph.2022.105261
White-Chu, E. F., and Reddy, M. (2011). Dry skin in the elderly: complexities of a common problem. Clin. Dermatol 29, 37–42. doi:10.1016/j.clindermatol.2010.07.005
Ye, L., Wang, H., Duncan, S. E., Eigel, W. N., and O'Keefe, S. F. (2015). Antioxidant activities of Vine Tea (Ampelopsis grossedentata) extract and its major component dihydromyricetin in soybean oil and cooked ground beef. Food Chem. 172, 416–422. doi:10.1016/j.foodchem.2014.09.090
Yin, Y., Morgunova, E., Jolma, A., Kaasinen, E., Sahu, B., Khund-Sayeed, S., et al. (2017). Impact of cytosine methylation on DNA binding specificities of human transcription factors. Science 356, eaaj2239. doi:10.1126/science.aaj2239
Yokochi, T., and Robertson, K. D. (2004). Doxorubicin inhibits DNMT1, resulting in conditional apoptosis. Mol. Pharmacol. 66, 1415–1420. doi:10.1124/mol.104.002634
Zhang, D., Xie, L., Jia, G., Cai, S., Ji, B., Liu, Y., et al. (2011). Comparative study on antioxidant capacity of flavonoids and their inhibitory effects on oleic acid-induced hepatic steatosis in vitro. Eur. J. Med. Chem. 46, 4548–4558. doi:10.1016/j.ejmech.2011.07.031
Zhang, J., Chen, Y., Luo, H., Sun, L., Xu, M., Yu, J., et al. (2018). Recent update on the pharmacological effects and mechanisms of dihydromyricetin. Front. Pharmacol. 9, 1204. doi:10.3389/fphar.2018.01204
Keywords: dihydromyricetin, DNMT1 inhibition, DNA methylation, DNAm age clock, skin, rejuvenation
Citation: Falckenhayn C, Bienkowska A, Söhle J, Wegner K, Raddatz G, Kristof B, Kuck D, Siegner R, Kaufmann R, Korn J, Baumann S, Lange D, Schepky A, Völzke H, Kaderali L, Winnefeld M, Lyko F and Grönniger E (2024) Identification of dihydromyricetin as a natural DNA methylation inhibitor with rejuvenating activity in human skin. Front. Aging 4:1258184. doi: 10.3389/fragi.2023.1258184
Received: 13 July 2023; Accepted: 12 December 2023;
Published: 04 March 2024.
Edited by:
Sigrid Le Clerc, Conservatoire National des Arts et Métiers (CNAM), FranceReviewed by:
Jing Liang, University of Southern California, United StatesLei Chen, Guangdong Ocean University, China
Copyright © 2024 Falckenhayn, Bienkowska, Söhle, Wegner, Raddatz, Kristof, Kuck, Siegner, Kaufmann, Korn, Baumann, Lange, Schepky, Völzke, Kaderali, Winnefeld, Lyko and Grönniger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elke Grönniger, elke.groenniger@beiersdorf.com