Abstract
Ambient air particulate matter (PM) induces senescence in human skin cells. However, the underlying mechanisms remain largely unknown. We investigated how epigenetic regulatory mechanisms participate in cellular senescence induced by PM with a diameter <2.5 (PM2.5) in human keratinocytes and mouse skin tissues. PM2.5-treated cells exhibited characteristics of cellular senescence. PM2.5 induced a decrease in DNA methyltransferase (DNMT) expression and an increase in DNA demethylase (ten–eleven translocation; TET) expression, leading to hypomethylation of the p16INK4A promoter region. In addition, PM2.5 led to a decrease in polycomb EZH2 histone methyltransferase expression, whereas the expression of the epigenetic transcriptional activator MLL1 increased. Furthermore, binding of DNMT1, DNMT3B, and EZH2 to the promoter region of p16INK4A decreased in PM2.5-treated keratinocytes, whereas TET1 and MLL1 binding increased, leading to decreased histone H3 lysine 27 trimethylation (H3K27Me3) and increased H3K4Me3 in the promoter of p16INK4A. PM2.5-induced senescence involved aryl hydrocarbon receptor (AhR)-induced reactive oxygen species (ROS) production. ROS scavenging dampened PM2.5-induced cellular senescence through regulation of DNA and histone methylation. Altogether, our work shows that skin senescence induced by environmental PM2.5 occurs through ROS-dependent the epigenetic modification of senescence-associated gene expression. Our findings provide information for the design of preventive and therapeutic strategies against skin senescence, particularly in light of the increasing problem of PM2.5 exposure due to air pollution.
Similar content being viewed by others
Introduction
Urban air pollution continues to be a serious health issue. Air pollution comprises a large proportion of particulate matter (PM) generated by coal combustion and diesel automobile exhaust1. Respiratory and cardiovascular diseases caused by PM are widely recognized. In addition, the skin, as the largest organ and outer layer of the human body, is constantly exposed to external stimuli and is thus a direct target of PM. Polycyclic aromatic hydrocarbons (PAHs) attached to PM represent important active agents that induce skin allergy, inflammatory dermatitis, eczema, and senescence2,3,4,5. PAHs such as benzo[a]pyrene induce toxicity in various cell types by activating aryl hydrocarbon receptor (AhR) signaling. AhR is a ligand-activated transcription factor. Upon ligand binding, AhR translocates from the cytosol to the nucleus and interacts with the AhR nuclear translocator to initiate transcription6. AhR activation can lead to oxidative stress as a result of excess generation of reactive oxygen species (ROS)7. These ROS and free radicals play a major role in skin senescence, inflammatory skin diseases, and carcinogenesis8,9,10.
Senescent cells remain metabolically active but cannot express genes required for cell proliferation. Senescent cells are typically characterized by a large and flat morphology, β-galactosidase activity, and senescence-associated heterochromatin foci (SAHF)11. The tumor suppressor protein p16INK4A is a key regulator of senescence and arrests the cell cycle by functioning as a specific inhibitor of cyclin-dependent kinases (CDKs) 4 and 612. In addition, the removal of p16INK4A-positive senescent cells has been found to delay senescence-associated disorders13.
Mammalian gene expression is tightly modulated by distinct epigenetic modifications, including DNA methylation and histone modifications. DNA methylation, which generally leads to transcriptional silencing, is one of the best-studied epigenetic modifications. DNA methylation is catalyzed by three DNA methyltransferases (DNMTs): DNMT1, DNMT3A, and DNMT3B, all of which add a methyl group to the C5 position of the cytosine ring of DNA, resulting in 5-methylcytosine (5-mC). This methylation process can be reversed by three DNA demethylases (ten–eleven translocations; TETs): TET1, TET2, and TET3. These enzymes can convert 5-mC to 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine, and 5-carboxylcytosine, ultimately yielding cytosine14.
Expression of the p16INK4a protein, encoded by the cyclin-dependent kinase inhibitor 2A gene15, is regulated by the DNA methylation of CpGs in its promoter, which reduces its expression16,17. In addition, the p16INK4A locus is silenced in actively growing cells by histone H3 lysine 27 trimethylation (H3K27Me3), which is induced by EZH2, a histone methyltransferase and transcriptional suppressor18. When cells are exposed to cellular stress, H3K27Me3 levels at the locus decrease, resulting in p16INK4A expression19. In contrast, MLL1 histone methyltransferase binds to the p16INK4A locus and induces the H3K4Me3 mark, activating p16INK4A transcription during replicative and premature senescence20.
Although many studies have led to a better understanding of the biochemical and cellular functions, as well as the epigenetic gene regulation of p16INK4A, the epigenetic regulation of p16INK4A expression during PM2.5-induced skin senescence remains poorly understood. Here, we focused on the epigenetic regulation of p16INK4A transcription and the molecular mechanisms by which they regulate p16INK4A transcription during PM2.5-induced cellular senescence in human keratinocytes and mouse skin tissues.
Materials and methods
PM2.5 preparation
Diesel particulate matter (PM2.5), which is a standard diesel PM (SRM 1650b) according to the National Institute of Standards and Technology (NIST, USA), was purchased from Sigma-Aldrich (St. Louis, MO, USA). With a mean diameter of 0.18 µm, 1650b diesel PM was composed predominantly of PAHs and nitro-PAHs. Information on the certified mass fraction values of PAHs and nitro-PAHs in SRM 1650b is available in our recently published paper21. PM2.5 was dissolved in dimethyl sulfoxide (DMSO) at a concentration of 100%, and when applied to cells, the concentration of DMSO did not exceed 0.01%.
Keratinocyte culture
The human keratinocyte cell lines HaCaT and HEK001 were obtained from the CLS Cell Lines Service (Eppelheim, Germany) and the American Type Culture Collection (Manassas, VA, USA), respectively. Normal human epidermal keratinocytes (NHEKs) derived from a foreskin specimen of a healthy male donor (18 years old) were obtained from Professor Jin Ho Chung (Seoul National University Hospital, Seoul, Republic of Korea). HaCaT cells were cultured in RPMI-1640 medium containing 10% heat-inactivated fetal calf serum at 37 °C in an incubator with a humidified atmosphere of 5% CO2. HEK001 cells were cultured in keratinocyte-SFM containing human recombinant EGF (Thermo Fisher Scientific, Inc., Waltham, MA, USA). Primary human HEKs were cultured in EpiLife serum-free medium with EpiLife growth supplement (Thermo Fisher Scientific, Inc.) and were used at the third or fourth passage. The procedures involving human subjects received prior approval from the Institutional Review Board of Seoul National University Hospital. The study was conducted according to the principles described in the Declaration of Helsinki.
Murine experiment
HR-1 hairless male mice weighing 22–24 g (OrientBio, Gyeonggi-do, Republic of Korea) were maintained under a controlled temperature of 25–28 °C and a 12 h light/dark cycle and received a standard diet and water ad libitum. All experimental procedures were conducted in accordance with the guidelines for the care and use of laboratory animals of Jeju National University (Jeju, Republic of Korea) (permit number: 2017-0026). Mice were divided randomly into groups of four animals each. PM2.5 was dissolved in 100 μg/mL propylene glycol and 100 μM NAC, spread on a nonwoven polyethylene pad with an area of 1 cm2, and then applied to the mouse dorsal region for 7 consecutive days. The skin tissue was immediately dissected for histological and molecular analyses22.
Detection of β-galactosidase activity
Cells were plated on coverslips, and SPiDER-βGal working solution (Dojindo Molecular Technologies, Inc., Rockville, MD, USA) was added to the plate for 15 min at 37 °C. The cells were mounted on microscope slides in mounting medium containing DAPI to label the nuclei. Images were collected using a Zeiss confocal microscope and Zeiss LSM 510 software.
Cell viability and colony formation
To assess cell viability, cells were seeded in a 60-mm dish at a density of 0.1 × 105 cells/mL and treated for 7 days with 50 μg/mL PM2.5 without medium change. The cells were then resuspended in a 0.4% trypan blue solution and incubated at room temperature for 15 min. Stained (dead) cells were counted under a light microscope. Next, to assess colony formation, cells were plated in 60-mm dishes at a density producing ~500 colonies per dish and treated with 50 μg/mL PM2.5 for 7 days without medium change. The resultant colonies were fixed with 75% ethanol and 25% acetic acid and stained with trypan blue. Colonies containing 50 or more cells were considered viable.
DAPI staining for SAHF detection
Cells were mounted on microscope slides in mounting medium containing DAPI to label the nuclei. Images were collected using a confocal microscope and Zeiss LSM 510 software.
Western blot analysis
Cell lysates (50 μg protein) were electrophoresed, transferred, and immunoblotted with specific antibodies. Anti-AhR, anti-TET1, and anti-TBP antibodies were purchased from Thermo Fisher Scientific, Inc.; the anti-DNMT1 antibody was from Abcam (Burlingame, CA, USA); the anti-DNMT3B antibody was from Santa Cruz Biotechnology (Santa Cruz, CA, USA); the anti-p16INK4A and anti-EZH2 antibodies were from Cell Signaling Technology, Inc. (Danvers, MA, USA), and the anti-MLL1 antibody was purchased from Active Motif (Carlsbad, CA, USA). The membranes with bound primary antibodies (1:1000) were reacted with horseradish peroxidase (HRP)-conjugated secondary antibodies (Pierce, Rockland, IL, USA), and the protein bands were assessed by a western blotting detection kit (Amersham, Little Chalfont, Buckinghamshire, UK).
Reverse transcription-PCR
The conditions for PCR using the forward primer 5′-GCAGCATGGAGCCTTCGGCT-3′ and the reverse primer 5′-TGCAGCACCACCAGCGTGTC-3′ for the amplification of p16INK4A were as follows: 5 min at 94 °C for initial denaturation; 35 cycles of 1 min at 94 °C, 1 min at 56 °C, and 1 min at 72 °C; and final elongation for 7 min at 72 °C. PCR amplification was performed in a programmable thermal cycler.
Immunofluorescence
After we fixed the cells with 4% paraformaldehyde for 30 min, they were permeabilized with PBS with 0.1% Triton X-100 for 2.5 min. The cells were treated with blocking medium (PBS containing 3% bovine serum albumin) for 1 h and incubated for 2 h with primary antibodies diluted in blocking medium. Primary antibody binding was assessed by an FITC- or Alexa594-bound secondary antibody (Santa Cruz Biotechnology) for 1 h. After the cells were mounted on microscope slides in mounting medium containing DAPI, they were imaged on a Zeiss confocal microscope using the LSM 510 program.
Chromatin immunoprecipitation and quantitative PCR
TET1, DNMT1, DNMT3B, EZH2, H3K27Me3, MLL1, and H3K4Me3 antibodies and normal rabbit IgG were used for chromatin immunoprecipitation (ChIP) assays, which were assessed with a specific ChIP kit (Cell Signaling Technology). DNA (100 ng) recovered from the immune-precipitated complexes was subjected to quantitative PCR (qPCR). The primers for the p16INK4A locus are described in Table 123.
Quantitative methylation-specific PCR and bisulfite sequencing
For quantitative methylation-specific PCR (qMSP) analysis, DNA (2 μg) was modified by bisulfite with a DNA methylation kit (Zymo Research, Orange, CA, USA), which converts nonmethylated cytosine to uracil and protects methylated cytosine. The methylation of gene promoter sites was assessed using qMSP primer pairs located close to the putative transcription start site in the 5′ CpG island. Bisulfite-treated DNA as the template and JumpStart REDTaq DNA polymerase (Sigma-Aldrich Co.) were used for amplification. For bisulfite sequencing, the PCR products were purified with a gel extraction kit (Qiagen GmbH, Hilden, Germany) and cloned using the TOPO TA vector system (Invitrogen, Carlsbad, CA, USA). Each clone was isolated and purified using a NucleoSpin plasmid isolation kit (Macherey-Nagel, Düren, Germany). Randomly selected positive clones (10–15 from each sample) were sequenced using the M13F primer, and the methylation status of each CpG dinucleotide was analyzed. For the quantification of p16INK4A methylation, qMSP amplification was performed on bisulfite-treated samples and normalized based on Alu element amplification. qPCR was performed using a CFX96TM real-time system (BioRad, Hercules, CA, USA). MSP and bisulfite sequencing primers were from published reports24,25.
Transfection of small interfering RNA
Control small interfering RNA (siRNA) (SS-1001), AhR sense (5′-UCAUGCAGCUGAUAUGCUU-3′) and antisense (5′-AAGCAUAUCAGCUGCAUGA-3′) siRNA, TET1#1 sense (5′-CAGUGUAACCAGCACAGUU-3′) and antisense (5′-AACUGUGCUGGUUACACUG-3′) siRNA, MLL1#1 sense (5′-GUCACAGUAGGUGAUCCUU-3′) and antisense (5′-AAGGAUCACCUACUGUGAC-3′) siRNA, and MLL1#2 sense (5′-CUAUUCUCGGGUCAUCAAU-3′) and antisense (5′-AUUGAUGACCCGAGAAUAG-3′) siRNA were purchased from Bioneer Corporation (Daejeon, Republic of Korea). TET1#2 siRNA (AM16708) was obtained from Invitrogen. Cells were transfected with Lipofectamine RNAiMax reagent (Invitrogen) in serum-free Opti-MEM. Then, 3 μL of lipofectamine was incubated with 1.5 mL of Opti-MEM for 5 min at room temperature. siRNA was added to the Opti-MEM-lipofectamine solution to a final concentration of 20 nM. The mixture was incubated for 5 min at room temperature. The transfection solution was incubated with the plates for 24 h.
Immunoprecipitation
Cell lysates (100 μg protein) were reacted overnight at 4 °C with MLL, TET1, or IgG antibodies. Immune complexes were collected with protein A/G PLUS beads (Santa Cruz Biotechnology) overnight at 4 °C and washed with immunoprecipitation buffer. Equal amounts of the precipitates were separated by SDS-polyacrylamide gel electrophoresis, followed by western blot analysis with antibodies specific for MLL1 and TET1.
Proximity ligation assay
The mouse/rabbit red starter Duolink kit (Sigma-Aldrich Co.) was used for this experiment. Cells were permeabilized according to the manufacturer’s instructions and incubated with goat anti-TET1 (1:100) or rabbit anti-MLL (1:1000) primary antibodies in fluorescence dilution buffer (5% fetal calf serum, 5% normal donkey serum, 2% bovine serum albumin in PBSCM, pH 7.6) for 2 h at room temperature and incubated with the proximity ligation assay (PLA) probes for 1 h at 37 °C in a humidified chamber. The ligation reaction was performed at 37 °C for 1 h in a humidified chamber. The cells were then incubated with the amplification mix for 2 h at 37 °C in a darkened humidified chamber, mounted using mounting medium, and imaged using a Zeiss confocal microscope and the LSM 510 program.
Detection of intracellular ROS
The intracellular levels of ROS were evaluated by spectrophotometric assay using an Amplex® Red kit (Invitrogen). Briefly, 0.1 × 105 cells/mL were seeded in a six-well plate. After 16 h of culture, the cells were treated with PM2.5. The cell pellet was resuspended in buffer (20 mM HEPES, 10 mM KCl, 1.5 mM MgCl2, 1 mM EDTA, 1 mM EGTA, 1 mM DTT, 0.1 mM PMSF, and protease inhibitors supplemented with 250 mM sucrose), and the cells were disrupted with a Dounce homogenizer. The reaction in the cell lysate was started with Amplex® Red solution prepared according to the manufacturer’s instructions (Amplex® Red reagent, 10-acetyl-3,7-dihydroxyphenoxazine in DMSO and 0.2 U/mL HRP). After incubation at 37 °C and 5% CO2 for 30 min, the absorbance at 560 nm was measured with a plate reader. To quantify ROS production, a standard curve for H2O2 was established in parallel. The results were expressed as H2O2 nmol released in 30 min per μg of protein.
Histological analysis
After skin sections (5 µm) were stained with hematoxylin and eosin, the thickness of the epidermis (from the stratum basale to the stratum corneum) was measured in ten randomly chosen fields taken from three representative sections per group, followed by microscopic evaluation at X100 optical magnification with a digital camera. Immunohistochemistry was performed using an ABC Elite kit (Vector Labs). The primary antibodies, including antikeratin 10 (Invitrogen), antiproliferating cell nuclear antigen (PCNA) (Santa Cruz Biotechnology), anti-DNMT3B, anti-TET1, anti-EZH2, anti-MLL1, and anti-p16INK4A antibodies, were incubated with the skin sections for 1 h. The sections were counterstained with hematoxylin before mounting.
Statistical analysis
The results are presented as the mean ± the standard error of the mean (SEM). SigmaStat software v12 (SPSS, Chicago, IL, USA) was used for statistical analysis. The data were analyzed using one-way analysis of variance (ANOVA) with Tukey’s post hoc test. Differences were considered statistically significant at p < 0.05.
Results
PM2.5 induces keratinocyte senescence
Recently, we investigated the potential role of oxidative stress induced by PM2.5 by measuring ROS generation at various concentrations (25–100 μg/mL) of PM2.5. At 50 μg/mL PM2.5, ROS levels, oxidative stress-induced cell damage and inflammatory responses were clearly observed21,26. Therefore, PM2.5 at a concentration of 50 μg/mL was used as the optimal concentration for the induction of senescence in cells.
The human keratinocyte cell lines HaCaT and HEK001 and NHEKs were treated with PM2.5 at 50 μg/mL, the optimal concentration for the induction of senescence in each of the cell types. PM2.5-treated cells were evaluated in various cellular senescence assays (Fig. 1a). PM2.5-treated cells exhibited β-galactosidase activity in the cytosol, a characteristic of cellular senescence, as evidenced by higher levels of green signal in PM2.5-treated cells than in untreated cells (Fig. 1b). Moreover, the expression of p16INK4A, a CDK inhibitor and senescence inducer, was considerably higher in PM2.5-treated cells than in the corresponding untreated cells (Fig. 1c). Furthermore, PM2.5-treated cells showed characteristics of cellular senescence, including an enlarged and flattened cell shape and irregular size (Fig. 1d), low proliferation (Fig. 1e), and decreased colony-forming ability (Fig. 1f). Chromatin in senescent cells undergoes large-scale rearrangement, forming dense nuclear domains called SAHF26. PM2.5-treated cells displayed significantly more SAHF-like chromatin foci in the nuclei than the controls (Fig. 1g).
P16INK4A is regulated epigenetically via DNA methylation during PM2.5-induced cellular senescence
The mRNA expression of p16INK4A, a senescence inducer, was upregulated in PM2.5-treated cells from 3 h onwards, continuing up to 72 h (Fig. 2a). In agreement with the RT-PCR results, western blotting also revealed strong induction of the p16INK4A protein (Fig. 2b). To investigate whether epigenetic DNA methylation could be involved in PM2.5-induced p16INK4A expression, we measured the expression of DNA methylation-related proteins. As shown in Fig. 2c, DNMT1 and DNMT3B expression decreased 3–12 h after PM2.5 treatment, whereas TET1 expression increased transiently at 3–12 h after PM2.5 treatment. The decrease in the expression of DNMT1 and DNMT3B and the increase in TET1 observed at 3 h after treatment by western blot analysis was also confirmed by immunofluorescence (Fig. 2d–f). Concomitant with the PM2.5-induced increase in TET1 and decrease in DNMT expression, we observed increased conversion of 5-mC to 5-hmC (Fig. 2g). These results suggest that these changes in DNMT1 and TET1 expression affect the transcription of senescence-associated genes in response to PM2.5. Using ChIP-qPCR analysis, we then assessed whether DNMT1, DNMT3B, and TET1 could directly bind to the p16INK4A locus in HaCaT cells and NHEKs (Fig. 3a). Binding of DNMT1 and DNMT3B to the p16INK4A promoter region occurred in the absence of PM2.5. This binding decreased in the proximal promoter regions of p16INK4A in HaCaT cells and NHEKs following PM2.5 treatment (Fig. 3b, c). The loss of DNMT1 and DNMT3B induced by PM2.5 prompted us to investigate whether this is accompanied by the appearance of TET1. We found that concomitant with the decreases in DNMT1 and DNMT3B, TET1 strongly binds to the B and C regions of the p16INK4A promoter region in PM2.5-treated HaCaT cells and NHEKs (Fig. 3d). Collectively, these findings suggest that PM2.5-mediated p16INK4A expression involves the recruitment of the DNA demethylase TET and the release of the DNMT, which participates in the induction of senescence.
We next examined the methylation status of the p16INK4A promoter using qMSP and bisulfite sequencing analysis. As shown in Fig. 3e, the DNA methylation of the p16INK4A promoter region was significantly lower in PM2.5-treated cells than in control cells. We then used bisulfite sequencing to examine DNA methylation in the p16INK4A promoter region before and after PM2.5 treatment. Interestingly, the promoter region (from −150 to +200 bp, including regions B and C for ChIP and 35 CpG sites of p16INK4A) showed lower methylation in PM2.5-treated cells than in control HaCaT cells (58% vs. 31%) and control NHEKs (55% vs. 43%) (Fig. 3f). Therefore, the change in the transcription of p16INK4A in PM2.5-treated HaCaT cells and PM2.5-treated NHEKs is likely due to epigenetic regulation through promoter methylation.
Because the increase in TET expression and hypomethylation of p16INK4A preceded the induction of p16INK4A during PM2.5-induced cellular senescence, we investigated whether the expression of p16INK4A depended on TET. TET1 siRNA reduced the expression of p16INK4A in PM2.5-treated cells, leading to a decrease in PM2.5-induced senescence (Fig. 3g, h). Thus, TET1 markedly contributes to the transcriptional activation of the p16INK4A locus, which participates in the induction of cellular senescence.
P16INK4A is regulated epigenetically via histone methylation in PM2.5-induced cellular senescence
To investigate whether epigenetic histone methylation could also be involved in PM2.5-induced p16INK4A expression, we measured the expression of the histone methyltransferases EZH2 and MLL1 in response to PM2.5 by western blot analysis. As shown in Fig. 4a, the expression of EZH2, a component of the polycomb complex with H3K27 methyltransferase activity, decreased 3–12 h after PM2.5 treatment, whereas the expression of MLL1, which is a transcriptional activator and has H3K4 methyltransferase activity, increased transiently at 3–12 h after PM2.5 treatment. The decrease in EZH2 and the increase in MLL1 expression were confirmed by immunofluorescence at 3 h after PM2.5 treatment (Fig. 4b, c). We then wondered whether these changes in EZH2 and MLL1 could participate in the induction of the p16INK4A locus in response to PM2.5. First, we assessed whether the silencer EZH2 occupies the p16INK4A locus by ChIP-qPCR analysis. We observed the binding of EZH2 and the presence of H3K27Me3 in the p16INK4A promoter region in the absence of PM2.5. The binding of EZH2 to the B and C regions of the p16INK4A promoter in HaCaT cells and to the C and D regions in NHEKs significantly decreased upon PM2.5 treatment in HaCaT cells and NHEKs (Fig. 4d) and was linked to a decrease in H3K27Me3 at the same positions (Fig. 4e). The loss of EZH2 and the H3K27Me3 suppressive histone mark induced by PM2.5 prompted us to investigate whether this is accompanied by the appearance of the MLL1 activator and H3K4Me3-positive histone mark. We found that concomitant with the decrease in EZH2 and H3K27Me3, the active MLL1 and H3K4Me3 mark strongly increased in the B and C regions of the p16INK4A promoter in PM2.5-treated HaCaT cells and NHEKs (Fig. 4f, g). Collectively, our results show that during PM2.5-mediated p16INK4A expression, the transcriptional activator MLL1 can replace the transcriptional silencer EZH2 in its proximal promoter.
Because the increase in MLL expression precedes the induction of p16INK4A during PM2.5-induced cellular senescence, we examined whether the expression of p16INK4A depended on MLL1. MLL1 siRNA reduced the expression of p16INK4A in cells treated with PM2.5 and limited the induction of cellular senescence (Fig. 4h, i). Thus, MLL significantly contributes to the transcriptional activation of the p16INK4A locus. We then evaluated whether TET1 interacts with MLL in response to PM2.5. Using both immunoprecipitation and PLA, we showed that the interaction between TET1 and MLL was greater in PM2.5-treated cells than in untreated cells (Fig. 4j, k).
PM2.5-induced cellular senescence is associated with oxidative stress via AhR-ROS signaling
PAHs interact directly with AhR, promoting AhR translocation into the nucleus and leading to ROS production27. We thus hypothesized that PM2.5-induced cellular senescence might occur via the AhR–ROS signaling pathway. Interestingly, immunoblotting data showed that the nuclear accumulation of AhR in PM2.5-treated cells occurred at 0.5 h and continued until 24 h (Fig. 5a). The nuclear accumulation of AhR was confirmed by immunofluorescence at 0.5 h (Fig. 5b). In response to PM2.5, the ROS level increased from 0.5 h, and the increase continued until 24 h in HaCaT cells and NHEKs, concomitant with AhR expression (Fig. 5c, d). To further identify the role of AhR in PM2.5-induced skin cellular senescence, we used siRNA against AhR. AhR siRNA transfection attenuated the increase in ROS at 0.5 h after PM2.5 treatment in HaCaT cells and NHEKs (Fig. 5e, f) and attenuated the increase in β-galactosidase activity, a classic marker of senescence (Fig. 5g). In addition to changes in gene expression and metabolic control during the cellular senescence process, the senescence rate is associated with the production of high levels of ROS28. PM2.5-treated HaCaT and HEK001 cells and NHEKs showed significantly higher intracellular ROS production at 24 h than control cells, and N-acetylcysteine (NAC), a well-known antioxidant, attenuated this increase (Fig. 5h). Interestingly, NAC treatment suppressed AhR expression and its translocation into the nucleus induced by PM2.5 (Fig. 5i, j). Furthermore, NAC pretreatment of HaCaT and HEK001 cells and NHEKs reduced senescence, as evidenced by a decrease in the β-galactosidase activity induced by PM2.5 treatment (Fig. 5k). Finally, unlike control mice, PM2.5-treated mice showed hyperkeratotic epidermis in their skin tissue, a skin senescence phenomenon29 (Fig. 5l). The hyperkeratotic epidermis induced by PM2.5 occurred via increases in the expression levels of keratin-10, which is a differentiation marker, and PCNA, which is a proliferation marker (Fig. 5m). However, NAC pretreatment of PM2.5-treated mice reduced this condition (Fig. 5l, m). Thus, PM2.5 induces skin cellular senescence, at least in part, by activating the AhR–ROS signaling pathway.
ROS are involved in the regulation of epigenetic modifiers by PM2.5
We investigated whether ROS induction by PM2.5 affects cellular senescence-related epigenetic regulators. As shown in Fig. 6a, NAC pretreatment of PM2.5-treated cells reversed the changes in DNMT, DNMT3B, TET1, EZH2, and MLL1 expression and limited the induction of p16INK4 expression. These changes in the expression levels of these enzymes resulted in similar changes in their respective binding to the p16INK4 promoter, as measured by ChIP-qPCR (Fig. 6b–f). Finally, changes in the different epigenetic enzymes and the induction of p16 were also observed in the epidermal keratinocytes of mice treated with PM2.5, and importantly, NAC reversed these changes (Fig. 6g).
Discussion
The skin is targeted directly by PM. PAHs present on the surface of PM cause harmful effects on the skin, including skin senescence, allergies, and inflammatory dermatitis. PM is known to induce the translocation of AhR to the nucleus, the activation of ERK and c-Jun, and the transcription of senescence-related proteins, such as CYP1A and MMP122. Recently, we reported that HaCaT keratinocytes internalized PM2.5, induced oxidative stress30, and triggered inflammatory cytokine production via the activation of NF-κB by TLR5-NOX4-ROS signaling21. However, the underlying mechanism that regulates skin cell senescence induced by PAH-mediated oxidative stress has not been analyzed in terms of epigenetic control, including DNA methylation and histone modification. In the present study, we demonstrated that PAHs on PM2.5 induced human skin epidermal keratinocyte senescence through the expression of p16INK4A, a senescence inducer, and the induction of oxidative stress.
P16INK4A is a specific CDK4/CDK6 inhibitor that plays a pivotal role in inducing cellular senescence. The p16INK4A protein is relatively stable, and its expression is primarily regulated at the transcriptional level. Aberrant hypermethylation of the p16INK4A gene is widely detected in most tumor types, leading to reduced gene expression31,32,33. In mammalian cells, DNMT1, DNMT3A, and DNMT3B participate in the maintenance of global DNA methylation and gene-specific de novo DNA methylation34. This methylation process can be reversed by DNA demethylases known as TETs. During PM2.5 exposure, TET replaces DNMT in the promoter region of p16INK4A, thereby promoting p16INK4A transcription. In addition to DNA methylation, histone modification via the histone methyltransferases EZH2 and MLL1 is also crucial for gene transcription. The H3K27 methyltransferase EZH2 binds to the p16INK4A locus, leading to the suppression of p16INK4A transcription. Conversely, the binding of the H3K4 methyltransferase MLL1 to the p16INK4A locus induces its transcription during replicative and premature senescence35. In the absence of PM2.5, EZH2 binds the proximal promoter (−0.3 to +85 bp) of p16INK4A in HaCaT cells and a broader region (from +85 to +5.6 kb) in NHEKs, probably reflecting differential regulation between these cell lines. However, upon PM2.5 exposure, MLL1 replaced EZH2 in both cell lines, thereby promoting p16INK4A transcription. Our previous studies showed that the DNA demethylase TET1 and the histone methyltransferase MLL1 interact to upregulate antioxidant transcription factors under oxidative stress36. This interaction also occurred upon PM2.5 exposure. Because MLL1- and TET1-binding sites overlap to a large extent with H3K4Me3 sites in the p16INK4A promoter, PM2.5-induced interactions between these proteins might participate in the induction of p16INK4A transcription.
ROS are known regulators of the methylation status of the p16INK4A promoter via JNK-DNMT1 pathway activation in the lung epithelium and in myelodysplastic syndrome37,38. PM2.5 induces AhR expression and translocation, leading to ROS generation. ROS contribute to the induction and maintenance of the cellular senescence process, as evidenced by the delayed onset of replicative senescence when ROS production is counteracted through antioxidants28. PM2.5-induced toxic effects in humans involve oxidative stress mechanisms. PM2.5 has been shown to be associated with DNA damage, increased ROS formation, and accelerated cellular senescence39. We also demonstrated that the antioxidant NAC attenuated the senescence induced by PM2.5-generated ROS. NAC pretreatment of PM2.5-treated cells reversed the changes in DNMT, DNMT3B, TET1, EZH2, and MLL1 expression, thereby reversing the decreased p16INK4 expression and enabling cellular senescence escape.
To the best of our knowledge, this is the first report that skin senescence induced by the PM2.5–AhR–ROS pathway involves epigenetic regulation of the expression of the senescence-associated p16INK4 gene (Fig. 6h). These results may provide insights into the development of preventive and therapeutic strategies against skin senescence, which may help address skin aging in the current scenario of ever-increasing air pollution.
References
Son, J. Y., Lee, J. T., Kim, K. H., Jung, K. & Bell, M. L. Characterization of fine particulate matter and associations between particulate chemical constituents and mortality in Seoul, Korea. Environ. Health Perspect. 120, 872–878 (2012).
Bartsch, N., Heidler, J., Vieth, B., Hutzler, C. & Luch, A. Skin permeation of polycyclic aromatic hydrocarbons: a solvent-based in vitro approach to assess dermal exposures against benzo[a]pyrene and dibenzopyrenes. J. Occup. Environ. Hyg. 13, 969–979 (2016).
Krutmann, J. et al. Pollution and skin: from epidemiological and mechanistic studies to clinical implications. J. Dermatol. Sci. 76, 163–168 (2014).
Pan, T. L. et al. The impact of urban particulate pollution on skin barrier function and the subsequent drug absorption. J. Dermatol. Sci. 78, 51–60 (2015).
Qiao, Y. et al. Airborne polycyclic aromatic hydrocarbons trigger human skin cells aging through aryl hydrocarbon receptor. Biochem. Biophys. Res. Commun. 488, 445–452 (2017).
Abiko, Y., Lin, F. Y., Lee, H., Puga, A. & Kumagai, Y. Quinone-mediated induction of cytochrome P450 1A1 in HepG2 cells through increased interaction of aryl hydrocarbon receptor with aryl hydrocarbon receptor nuclear translocator. J. Toxicol. Sci. 41, 775–781 (2016).
Costa, C. et al. Exposure of human skin to benzo[a]pyrene: role of CYP1A1 and aryl hydrocarbon receptor in oxidative stress generation. Toxicology 271, 83–86 (2010).
Birch-Machin, M. A. & Bowman, A. A. Oxidative stress and ageing. Br. J. Dermatol. 2, 26–29 (2016).
Kim, K. E., Cho, D. & Park, H. J. Air pollution and skin diseases: adverse effects of airborne particulate matter on various skin diseases. Life Sci. 152, 126–134 (2016).
Venza, M. et al. Cellular mechanisms of oxidative stress and action in melanoma. Oxid. Med. Cell. Longev. 2015, 481782 (2015).
Zhao, L. et al. JMJD3 promotes SAHF formation in senescent WI38 cells by triggering an interplay between demethylation and phosphorylation of RB protein. Cell Death Differ. 22, 1630–1640 (2015).
He, S. & Sharpless, N. E. Senescence in health and disease. Cell 169, 1000–1011 (2017).
Baker, D. J. et al. Clearance of p16Ink4a-positive senescent cells delays ageing-associated disorders. Nature 479, 232–236 (2011).
Ito, S. et al. Tet proteins can convert 5-methylcytosine to 5-formylcytosine and 5-carboxylcytosine. Science 333, 1300–1303 (2011).
Zhao, R., Choi, B. Y., Lee, M. H., Bode, A. M. & Dong, Z. Implications of genetic and epigenetic alterations of CDKN2A (p16(INK4a)) in cancer. EBioMedicine 8, 30–39 (2016).
Song, W. et al. Increased p16 DNA methylation in mouse thymic lymphoma induced by irradiation. PLoS One 9, e93850 (2014).
Zhu, B. et al. Atorvastatin treatment modulates p16 promoter methylation to regulate p16 expression. FEBS J. 284, 1868–1881 (2017).
D’Arcangelo, D., Tinaburri, L. & Dellambra, E. The role of p16INK4a pathway in human epidermal stem cell self-renewal, aging and cancer. Int. J. Mol. Sci. 18, E1591 (2017).
Serra, R. W., Fang, M., Park, S. M., Hutchinson, L. & Green, M. R. A KRAS-directed transcriptional silencing pathway that mediates the CpG island methylator phenotype. Elife 3, e02313 (2014).
Kotake, Y., Zeng, Y. & Xiong, Y. DDB1-CUL4 and MLL1 mediate oncogene-induced p16INK4a activation. Cancer Res. 69, 1809–1814 (2009).
Ryu, Y. S. et al. Particulate matter induces inflammatory cytokine production via activation of NFκB by TLR5-NOX4-ROS signaling in human skin keratinocyte and mouse skin. Redox Biol. 21, 101080 (2019).
Lee, C. W. et al. Urban particulate matter down-regulates filaggrin via COX2 expression/PGE2 production leading to skin barrier dysfunction. Sci. Rep. 6, 27995 (2016).
Kia, S. K., Gorski, M. M., Giannakopoulos, S. & Verrijzer, C. P. SWI/SNF mediates polycomb eviction and epigenetic reprogramming of the INK4b-ARF-INK4a locus. Mol. Cell. Biol. 28, 3457–3464 (2008).
Herman, J. G. et al. a novel PCR assay for methylation status of CpG islands. Proc. Natl Acad. Sci. USA 93, 9821–9826 (1996).
McGarvey, K. M., Greene, E., Fahrner, J. A., Jenuwein, T. & Baylin, S. B. DNA methylation and complete transcriptional silencing of cancer genes persist after depletion of EZH2. Cancer Res. 67, 5097–5102 (2007).
Aird, K. M. & Zhang, R. Detection of senescence-associated heterochromatin foci (SAHF). Methods Mol. Biol. 965, 185–196 (2013).
Zhou, B. et al. Mitochondrial activity and oxidative stress functions are influenced by the activation of AhR-induced CYP1A1 overexpression in cardiomyocytes. Mol. Med. Rep. 16, 174–180 (2017).
Davalli, P., Mitic, T., Caporali, A., Lauriola, A. & D’Arca, D. ROS, cell senescence, and novel molecular mechanisms in aging and age-related diseases. Oxid. Med. Cell. Longev. 2016, 3565127 (2016).
Singh, B., Schoeb, T. R., Bajpai, P., Slominski, A. & Singh, K. K. Reversing wrinkled skin and hair loss in mice by restoring mitochondrial function. Cell Death Dis. 9, 735 (2018).
Piao, M. J. et al. Particulate matter 2.5 damages skin cells by inducing oxidative stress, subcellular organelle dysfunction, and apoptosis. Arch. Toxicol. 92, 2077–2091 (2018).
Emmett, R. A., Davidson, K. L., Gould, N. J. & Arasaradnam, R. P. DNA methylation patterns in ulcerative colitis-associated cancer: a systematic review. Epigenomics 9, 1029–1042 (2017).
Han, Y. D. et al. Associations of P16INK4a promoter hypermethylation with squamous intra-epithelial lesion, cervical cancer and their clinicopathological features: a meta-analysis. Oncotarget 8, 1871–1883 (2017).
Nikolic, N. et al. High frequency of p16 and p14 promoter hypermethylation and marked telomere instability in salivary gland tumors. Arch. Oral. Biol. 60, 1662–1666 (2015).
Zhou, W. et al. Repeated PM2.5 exposure inhibits BEAS-2B cell P53 expression through ROS-Akt-DNMT3B pathway-mediated promoter hypermethylation. Oncotarget 7, 20691–20703 (2016).
Kotake, Y., Naemura, M., Murasaki, C., Inoue, Y. & Okamoto, H. Transcriptional regulation of the p16 tumor suppressor gene. Anticancer Res. 35, 4397–4401 (2015).
Kang, K. A. et al. Interaction of DNA demethylase and histone methyltransferase upregulates Nrf2 in 5-fluorouracil-resistant colon cancer cells. Oncotarget 7, 40594–40620 (2016).
Gonçalves, A. C. et al. Oxidative stress levels are correlated with P15 and P16 gene promoter methylation in myelodysplastic syndrome patients. Clin. Exp. Med. 116, 333–343 (2016).
Soberanes, S. et al. Particulate matter air pollution induces hypermethylation of the p16 promoter via a mitochondrial ROS-JNK-DNMT1 pathway. Sci. Rep. 2, 275 (2012).
Gao, Z. X. et al. Assessment of DNA damage and cell senescence in corneal epithelial cells exposed to airborne particulate matter (PM2.5) collected in Guangzhou, China. Investig. Ophthalmol. Vis. Sci. 57, 3093–3102 (2016).
Acknowledgements
This work was supported by a grant from the Basic Research Laboratory Program (NRF-2017R1A4A1014512) of the National Research Foundation of Korea (NRF) funded by the Korean government (MSIP).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ryu, Y.S., Kang, K.A., Piao, M.J. et al. Particulate matter-induced senescence of skin keratinocytes involves oxidative stress-dependent epigenetic modifications. Exp Mol Med 51, 1–14 (2019). https://doi.org/10.1038/s12276-019-0305-4
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s12276-019-0305-4
This article is cited by
-
Tubular aryl hydratocarbon receptor upregulates EZH2 to promote cellular senescence in cisplatin-induced acute kidney injury
Cell Death & Disease (2023)
-
Molecular mechanisms of environmental exposures and human disease
Nature Reviews Genetics (2023)
-
Aging of mesenchymal stem cell: machinery, markers, and strategies of fighting
Cellular & Molecular Biology Letters (2022)
-
Staphylococcus epidermidis WF2R11 Suppresses PM2.5-Mediated Activation of the Aryl Hydrocarbon Receptor in HaCaT Keratinocytes
Probiotics and Antimicrobial Proteins (2022)
-
Lnc-IL7R alleviates PM2.5-mediated cellular senescence and apoptosis through EZH2 recruitment in chronic obstructive pulmonary disease
Cell Biology and Toxicology (2022)