- 1Department of Dermatology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2Department of Gastroenterology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
Vitiligo is an pigmentation disorder caused by a variety of pathogenic factors; its main pathophysiological conditions include oxidative stress, immune activation, and genetic background. Additionally, DNA methylation is often associated with the pathogenesis of vitiligo; however, the underlying mechanism remains unknown. In the present study, we used the Human Methylation 850K BeadChip platform to detect DNA methylation changes in the vitiligo melanocytes. We then integrated the results with the transcriptome data of vitiligo melanocytes and lesions to analyse the correlation between differentially methylated levels and differentially expressed genes. The results showed that there was a significant negative correlation between methylation levels and differentially expressed genes. Subsequently, we enriched GO and KEGG based on methylated differentially expressed genes (MDEGs) using R package ClusterProfiler, and the results were closely related to the pathogenesis of vitiligo. In addition, we also constructed a PPI network of MDEGs and excavated three important functional epigenetic modules, involving a total of 12 (BCL2L1, CDK1, ECT2, HELLS, HSP90AA1, KIF23, MC1R, MLANA, PBK, PTGS2, SOX10, and TYRP1) genes. These genes affect melanocyte melanogenesis, cellular oxidative stress and other important biological processes. Our comprehensive analysis results support the significant contribution of the status of DNA methylation modification to vitiligo, which will help us to better understand the molecular mechanism of vitiligo and explore new therapeutic strategies.
Introduction
Vitiligo is an acquired chronic skin pigmentation disease that affects 0.5–2% of the world's population (1). The reported epidemiological data varies according to region, and the disease impact may be related to the respective social and cultural stigmas (2). The pathogenesis of vitiligo manifests as oxidative stress, infiltration of inflammatory mediators, melanocyte loss, and autoimmune responses (3). It is caused due to the complex interaction among environmental and genetic factors that eventually leads to melanocyte dysfunction (4). However, the inherent defect in melanocytes is an early event in vitiligo; the gradual disappearance of melanocytes may involve multiple processes, including immune system attack, cell degeneration, and exfoliation among others (1).
Epigenetic factors such as DNA methylation, histone modification, and gene silencing through microRNAs play an important role in the pathogenesis and progression of autoimmune dermatoses, such as psoriasis and atopic dermatitis (5). Increasing number of studies have shown that failure to maintain the level and pattern of DNA methylation can lead to abnormal cell function and proliferative activity (6). Although epigenetics plays a key role in the pathogenesis of skin tumors, researchers have recently begun to pay more attention to the epigenetics in the pathogenesis of inflammatory skin diseases, such as psoriasis, atopic dermatitis, and other inflammatory skin diseases (7–9). Abnormal DNA methylation has been identified as a major epigenetic change leading to the development of various skin cancers (10). However, DNA methylation is reversible and is therefore considered an attractive therapeutic intervention (11, 12) to provide the possibility treatment methods for various skin diseases.
Evidences have shown that abnormal DNA methylation is involved in the development of vitiligo. For example, the DNA methylation level of peripheral blood mononuclear cells (PBMC) in vitiligo have been investigated and the emerging evidence suggests an epigenetic regulation of CD8+ T cells. Zhao M et al. reported that in vitiligo patients, DNA methylation of peripheral blood mononuclear cells affects the mRNA expression level of DNMT, MBD, IL-10, and other genes (13). Another study suggests that in autoimmune diseases such as type 1 diabetes (T1D), systemic lupus erythematosus (SLE), and vitiligo genes associated with the proliferation or activation of CD8+ T cells are affected by epigenetic modification (14). However, the important biological pathway regulated by abnormal methylation has not been elucidated yet, and the role of epigenetic factors in the vitiligo pathogenesis needs further research.
In the present study, we integrated and analyzed the differential transcriptomic data (from vitiligo melanocyte cell lines and vitiligo lesions) and differential methylation data (from vitiligo melanocyte cell lines). We used Infinium Methylation EPIC BeadChip for the methylation sequencing of vitiligo melanocyte PIG3V and normal melanocyte PIG1. The 850K microarray covered the gene promoter region, gene coding region, CpG island, and enhancer regions present in ENCODE and FANTOM5. The data of the transcriptional group were derived from vitiligo lesion transcriptional sequencing datasets GSE75819 (15 non-segmental vitiligo patients with lesion skin and peri-lesion skin biopsies, the peri-lesion skin biopsies as controls) in the GEO database and our cell line transcriptional sequencing data (PIG3V and PIG1). Subsequently, we conducted GO/KEGG pathway analysis, PPI, and functional epigenetic module analysis of MDEGs to further explore the potential mechanism driven by abnormal methylation of vitiligo. Flowchart of bioinformatics analysis is shown in Figure 1.
Figure 1. Flowchart of bioinformatics analysis. DMGs, Differentially methylated genes; DEGs, Differentially expressed genes; PIG1, human normal melanocyte; PIG3V, human vitiligo melanocyte.
Materials and Methods
Study Subjects
The human vitiligo melanocyte cell line PIG3V and the human normal melanocyte cell line PIG1 (gifted by Dr. Li Chunying, Xijing Hospital, Fourth Military Medical University, Xi'an, Shanxi Province) were cultured in human melanocyte growth medium 254 supplemented with 5% fetal bovine serum and penicillin- streptomycin at 37°C in a 5% CO2 incubator. To integrate and study the correlation between the methylation status and transcriptome dataset, we sequenced the whole genome of two cases of PIG1 cell line and two cases of PIG3V cell line, and used it as cell line methylation dataset and expression dataset. Vitiligo lesion expression GSE75819 was downloaded from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/). Fifteen vitiligo focal skin samples and fifteen vitiligo peri-focal skin samples were included in the GSE75819 to verify the expression dataset.
Infinium Human Methylation EPIC Array Using 850K BeadChip and Data Analysis
DNA was isolated from cell samples, using DNeasy Blood & Tissue Kit (Qiagen, Germany). The purity and concentration of DNA was estimated using Nanodrop 2000 (Thermo Fisher Scientific, China). Approximately 500 ng of genomic DNA from each sample was used for sodium bisulfite treatment using the EZ DNA methylation Gold Kit (Zymo Research, USA) following the manufacturer's protocol. Genome-wide DNA methylation was assessed using the Infinium Human Methylation 850K BeadChip (Illumina Inc, USA) according to the manufacturer's instructions. The array data (IDAT files) were analyzed using the ChAMP package in R for deriving the level of methylation. The methylation status of all the probes was denoted as β value, which is the ratio of the methylated probe intensity to the overall probe intensity (sum of methylated and unmethylated probe intensities plus constant α, where α = 100). A CpG site with |Δβ| ≥ 0.20 (in test vs. control) and adjusted P-value ≤ 0.05 was considered as a differentially methylated site. A CpG was considered hypermethylated if Δβ ≥ 0.20 or hypomethylated if Δβ ≤ −0.20. The average β value of promoters and CpG islands was compared between diseased and normal conditions. Promoters and CGIs with |Δβ| ≥ 0.20 and adjusted P-value ≤ 0.05 were considered for further analysis.
Copy Number Variation Analysis
We used the consumer package to analyse the human measurement epic data to identify copy number variation (CNV). CNV analysis was carried out in two steps. First, the combined signal values of “methylated” and “unmethylated” CPGs were normalized using the control group (PIG1 as control). This step required correcting the probe and sample deviation. Second, the adjacent probes were combined to generate a probe bin with the minimum size and quantity to reduce the remaining technical changes and achieve meaningful segmentation results. For CNV results and differential methylation sites, we used the Circos map for joint analysis.
Acquisition and Processing of Expression Spectrum Data
RNA-Seq strand-specific libraries were constructed using a VAHTS Total RNA-Seq (H/M/R) Library Prep Kit (Vazyme, China). The original gene expression data set GSE75819 was downloaded from the GEO public database. We used R 4.0.0 to analyse the expression matrix. The robust multiarray average (RMA) method was used to pre-process data, including background adjustment, normalization, and log value conversion. The Limma package was used to search for differential genes. The threshold of up and down genes was set as |log FC| ≥ 1 and P ≤ 0.05.
FunRich
FunRich 3.1.3 (http://www.funrich.org) is an independent software tool, used mainly for the functional enrichment and interaction network analysis of genes and proteins (15). The software was used to identify hypomethylation and hypermethylation genes for further analysis.
GO and KEGG Enrichment
GO and KEGG enrichment were carried out for hypomethylation and hypermethylation genes using R package clusterProfiler. Results with P < 0.05 were considered statistically significant. We used R package ggplot2 to visualize the important projects of each group function and pathway enrichment analysis.
Construction of PPI Network and Related Analysis
To further explain the interactions between the two sets of genes in the pathogenesis of vitiligo and the specific molecular mechanisms, PPI analysis was performed. Since the number of hypomethylation and hypermethylation genes was small, we integrated them to construct PPI networks. STRING 11.0 was used to generate a PPI network (https://string-db.org/). The cutoff value was set to 0.4 of the interaction score. The results were imported into Cytoscape 3.8.0 for subsequent analysis. CytoHubba app (http://apps.cytoscape.org/apps/cytohubba) was used to screen the hub genes. The first 3 modules were applied in Cytoscape by the MCODE app software package. We used FunRich 3.1.3 to identify and visualize the potential genes related to each module.
Statistical Analysis
The correlation between differentially expressed genes and methylation levels was analyzed by Pearson's correlation coefficient. Results with P < 0.05 were considered statistically significant. The statistical analyses were performed using R 4.0.0. GraphPad Prism 8 (GraphPad Prism Software Inc., San Diego, California) to correlate analysis and graphic display.
Results
Global Methylation Pattern of Vitiligo Melanocytes and Normal Melanocytes
The characteristics of DNA methylation were significantly different between vitiligo melanocyte cell line and normal melanocyte cell line. The analysis of differential methylation position (DMP) showed 83 451 differential methylation positions in PIG3V compared with PIG1. Among them, 53 895 were hypermethylated and 29 556 were hypomethylated. CpGs showed bimodal methylation distribution in PIG3V, and most of them showed hyper/hypo methylation levels. The differential methylation positions are shown in the cluster heat map (Figure 2A) and the volcano map (Figure 2B). The genes corresponding to the top 20 differences in methylation positions are listed at the two tips of the volcanic map. In addition, to examine the methylation changes throughout the genome, we used Circos maps for joint analysis (Figure 2C).
Figure 2. (A) The heat map, (B) the Volcano map analysis of genes with differential methylation between vitiligo melanocyte and normal melanocyte. (C) A joint analysis of CNV results and differential methylation sites performde by the Circos diagram. CNV and the outermost circle represent chromosome regions, and the second inward circle represents the CNV copy number (red represents ratio > 0, green represents ratio < 0), The inner circle represents the DMR area.
Integration of Epigenomic and the Transcriptomic Data
The intersection of transcriptomic data on vitiligo cell lines and lesions was visualized in the veen map (Figure 3A); there were 117 down-regulated and 149 up-regulated genes. The intersection of transcriptomic methylation data of vitiligo melanocyte cell line was also visualized by veen map (Figure 3B) and 126 overlapping genes were obtained. We analyzed the correlation between intersectional genes and their methylation levels, and the results showed that DNA methylation in vitiligo melanocyte cell line is negatively correlated with gene expression (Figure 3C), but there was no significant positive correlation (Figure 3D). Therefore, we selected the genes with a high expression of low methylation and low expression of high methylation for further study. There were 69 genes that met these requirements, including 14 down-regulated genes with hypermethylation and 55 up-regulated genes with low methylation. It is important to note that GEO currently does not have DNA methylation data for vitiligo, so there is no verification set for vitiligo methylation.
Figure 3. Integrative analysis of methylome and transcriptome of vitiligo. (A) Veen diagram of differentially expressed genes between vitiligo skin lesion microarray and vitiligo melanocyte microarray. (B) Venn diagram between differentially methylated genes and differentially expressed genes of vitiligo. Correlation analysis between differentially methylated genes and differentially expressed genes: (C) negative correlation analysis, (D) positive correlation analysis.
Enrichment Analysis of GO and KEGG Pathways
In PIG3V, the biological processes of MDEGs, such as oxidation–reduction oxidoreductase activity, are mainly related to oxidative stress, immune responses, and melanogenesis (Figures 4A,B). The response to oxidative stress is related to the positive regulation of T cell proliferation, whereas regulation of acute inflammatory responses is involved in inflammation, and melanocyte differentiation and melanosome transport affect melanogenesis. The specific genes involved in each GO term are shown in the Supplementary Table 1. Apart from these, some other processes such as the regulation of cell aging, response to vitamin D, autophagy, and the positive regulation of fatty acid metabolic processes are related to the pathogenesis of vitiligo.
Figure 4. GO and KEGG enrichment analysis of methylated-differentially expressed genes related with vitiligo melanocyte. Top 10 GO terms in (A) GO enrichment analysis of Hyper-MDGEs, (B) GO enrichment analysis of Hypo-MDGEs. (C) KEGG pathway analysis of Hyper-MDGEs, (D) KEGG pathway analysis of Hypo-MDGEs. The horizontal axes shows -log10 transformed P-value and p < 0.05 is considered significant.
KEGG analysis identified a variety of pathways in enrichment analysis, and we visualized the top 10 (Figures 4C,D). According to KEGG classification, the pathways rich in MDEGs are commonly involved in IL-17 signaling, Th17 cell differentiation, TNF signaling, and NOD-like receptor signaling pathway. Further, different pathways may overlap due to common elements. For the transcripts of differential regulation, some pathways are concentrated in the glucagon signaling such as insulin resistance, insulin signaling, and thyroid hormone signaling. This may be related to the crosstalk among vitiligo, diabetes, thyroid, and other autoimmune diseases, and the MDEGs involved in these pathways may be important biomarkers for the mutual crosstalk in immune diseases.
PPI Network and Functional Epigenetic Modules
The PPI network analyzed by STRING showed 69 nodes and 66 edges. Cytoscape 3.8.0 was used to visualize the PPI network (Figure 5A), CytoHubba app in Cytoscape software is used to select the hub nodes genes of the PPI network. In the results, a total of 7 genes of the top 10 hub genes detected by six ranked methods in cytoHubba were overlap hub genes (Figure 5B), including CDK1 (cyclin dependent kinase 1), HSP90AA1 (heat shock protein 90 alpha family class A member 1), AKT1 (AKT serine/threonine kinase 1), BCL2L1 (BCL2 like 1), HDAC2 (histone deacetylase 2), HELLS (helicase, lymphoid specific), and KIF23 (kinesin family member 23). Among them, BCL2L1 was hypermethylated and its expression was downregulated, whereas the other six genes are hypomethylated and their expression was upregulated. Finally, based on the PPI network, three important epigenetic function modules were identified and were established as statistically significant. The results of enrichment analysis of functional epigenetic modules showed that the MDEGs of module 1 (Figure 5C) played an important role in the regulation of cell cycle and cell division, and the MDEGs of module 2 (Figure 5D) affected melanogenesis and pigmentation. Module 3 (Figure 5E) involved the regulation of cellular response to stress, negative regulation of intrinsic apoptotic signaling pathway, and complex pathway regulation. The enrichment results are presented in detail in Table 1.
Figure 5. (A) Protein-protein interaction network of MDEGs, Disconnected nodes were hid in the network. (B) Hub genes for MDEGs ranked in cytoHubba, functional epigenetic modules of the protein-protein interaction network: (C) module1, (D) module2, (E) module3.
Table 1. GO and KEGG enrichment analysis of functional epigenetic modules related with vitiligo melanocyte.
Discussion
Vitiligo is considered to be an autoimmune disease because cytokines such as (IFN)-γ, IL-1, IL-6, IL-8, and IL-10 are overexpressed in lesions, and activated CD8 + T lymphocytes, TH17, and other immune cells are significantly aggregated in the lesion area (16–19). However, some researchers have suggested that autoimmunity may be secondary to high oxidative stress in vitiligo skin, leading to inherent defects of melanocytes and changes in its microenvironment (4). At present, the unified view is that the inherent defect of melanocytes is an early event in vitiligo, and the gradual disappearance of epidermal melanocytes leads to skin depigmentation of vitiligo lesions.
GO enrichment analyses of the oxidation-reduction process, response to oxidative stress, and oxidoreductase activity have been partially studied in vitiligo (20–23). Inflammation-related processes such as positive regulation of T cell proliferation, regulation of acute inflammatory response, and melanogenesis-related processes such as cellular pigmentation, melanocyte differentiation, melanosome transport were also widely studied in vitiligo (24–26). However, previous studied did not consider methylation modification, and our study shows that the methylation levels of key genes involved in these important regulatory processes, such as TYP1, IL17, and MC1R, are altered in vitiligo melanocyte PIG3V. Therefore, based on our findings, we propose that the change in methylation level may regulate the differential expression of functional genes in vitiligo. To the best of our knowledge, this is the first study to analyze the methylation profile of vitiligo melanocytes.
The infiltration of immune cells and the release of pro-and anti-inflammatory cytokines is key to the pathogenesis of vitiligo. KEGG results showed that MDEGs were significantly enriched in IL-17 signaling pathway, Th17 cell differentiation, TNF signaling pathway, and NOD-like receptor signaling pathways, and these signaling pathways have been proved to be involved in mediating the immune regulation of vitiligo (18, 27, 28). In addition, differential methylation of AKT1, PYGB, HDAC2 were frequently observed in the insulin signaling and thyroid hormone signaling pathways. These pathways play an important role in abnormal glucose metabolism in diabetes and thyroid diseases (29–31), while vitiligo is frequently associated with other autoimmune diseases, particularly autoimmune thyroid diseases (Hashimoto's thyroiditis and Graves' disease), adult-onset type 1 diabetes mellitus, psoriasis, and so on (32). Thus, genes in these signaling pathways are regulated by methylation modification may play an important role in crosstalk between immune diseases. Previous studies have shown that the effective treatment of psoriasis, an immune inflammatory skin disease, can reverse DNA methylation in the epidermis (33, 34), suggesting that further studies on vitiligo methylation may be potentially identify new treatment strategies for vitiligo.
Futhermore, the PPI network showed that CDK1, HSP90AA1, AKT1, BCL2L1, HDAC2, HELLS, and KIF23 are 7 central genes identified by the top 10 high hub nodes after gene overlap. Among these central genes, CDK1 is a set of Ser/Thr kinase system corresponding to cell cycle progression (35), and its differential expression was found in another vitiligo study (36); HSP90AA1 plays a key role in signal transduction, protein folding, protein degradation, and morphological evolution (37); AKT1 and BCL2L1 are involved in apoptosis (38, 39), and AKT methylation has been shown to promote skin tumors (40); HDAC2 interferes with histone deacetylase function (41); KIF23 plays an important role in the changes of cellular motor cytoskeleton (42); HELLS is considered to be a regulator of cell senescence and is necessary for DNMT1-mediated methylation maintenance and DNMT3A/DNMT3B-mediated remethylation (43). Therefore, we speculate that these hub genes may be involved in the pathological process in the early stage of melanocyte loss in vitiligo, but further functional experiments are needed to further explore.
The analysis of the epigenetic module shows that the MDEGs in module 1 is highly enriched in the processes of cell cycle regulation, whereas the MDEGs in module 2 is involved in the process of pigmentation and the pathway of regulating melanin production. Interestingly, the functions of MDEGs in these modules are relatively concentrated, which suggests that these biological processes may be collectively regulated by methylation. Whether demethylation can reverse these processes needs to be verified by further studies. Module 3 involves cellular stress, apoptosis, and regulation of complex biological processes and pathways (44, 45) key to the pathogenesis of vitiligo. These findings suggest that MDEGs in vitiligo melanocytes may have regulatory functions in these biological processes and molecular functions. However, some genes and pathways identified in this study have not been formally studied as targets of the vitiligo process, and need to be evaluated further.
Although our research attempts to bridge an important gap, it does have several limitations that can be addressed in the future studies. Firstly, the sample size was relatively small as only four microarray profiles were analyzed and there was no vitiligo methylation data available on the GEO platform. Hence, a replication with larger samples will be required to validate the findings. Secondly, the study lacked experimental verification of the effects of aberrant methylation on gene expression and functions in vitiligo melanocytes. Therefore, supplementary molecular experiments should be conducted to verify the results of our investigation. Furthermore, due to the limited sample size, the GEO verification set of differentially expressed genes selected only differential datasets rather than general high and low expression sets due to the differences between cell lines and human samples. This strict intersection may mask genes that play significant roles.
In conclusion, using a series of bioinformatics databases and tools, we found that the interactions among MDEGs with different functions and signaling pathways is related to the pathogenesis of vitiligo melanocytes. The hub genes of vitiligo melanocytes include CDK1, HSP90AA1, AKT1, BCL2L1, HDAC2, HELLS, and KIF23. The genes involved in three important epigenetic modules include BCL2L1, CDK1, ECT2, HELLS, HSP90AA1, KIF23, MC1R, MLANA, PBK, PTGS2, SOX10, and TYRP1. This study provides hypothetical and biological characteristic insight into the pathogenesis of vitiligo. However, additional molecular-level studies are needed to confirm the identified genes and pathways in vitiligo to elucidate potential mechanisms, and in vitro and in vivo functional studies are also required to find the crucial role of the identified genes in vitiligo pathogenesis.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
YP and XC prepared the figures and writing. YC and LZ collected and organized the data. JiaC cultured cells. YZ and XS prepared the tables. JinC critically revised the data. All authors contributed to the article and approved the submitted version.
Funding
This work was financially supported by the National Natural Science Foundation of China (No. 81773307) and Chongqing Science and Technology Commission (No. cstc2018jcyjAX0195), which offered the financial support.
Conflict of Interest
The 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.
Acknowledgments
We would like to thank the Sinotech Genomics Biotechnology Company for technical support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.587440/full#supplementary-material
References
1. Picardo M, Dell'Anna ML, Ezzedine K, Hamzavi I, Harris JE, Parsad D, et al. Vitiligo. Nat Rev Dis Primers. (2020) 312:461–6.
2. Alikhan A, Felsten LM, Daly M, Petronic-Rosic V. Vitiligo: a comprehensive overview Part I. Introduction, epidemiology, quality of life, diagnosis, differential diagnosis, associations, histopathology, etiology, and work-up. J Am Acad Dermatol. (2011) 65:473–91. doi: 10.1016/j.jaad.2010.11.061
3. Schallreuter KU, Bahadoran P, Picardo M, Slominski A, Elassiuty YE, Kemp EH, et al. Vitiligo pathogenesis: autoimmune disease, genetic defect, excessive reactive oxygen species, calcium imbalance, or what else? Exp Dermatol. (2008) 17:139–60. doi: 10.1111/j.1600-0625.2007.00666.x
4. Abdel-Malek ZA, Jordan C, Ho T, Upadhyay PR, Fleischer A, Hamzavi I. The enigma and challenges of vitiligo pathophysiology and treatment. Pigment Cell Melanoma Res. (2020) 33:778–87. doi: 10.1111/pcmr.12878
5. Mervis JS, McGee JS. DNA methylation and inflammatory skin diseases. Arch Dermatol Res. (2020) 312:461–6. doi: 10.1007/s00403-019-02005-9
6. Gopisetty G, Ramachandran K, Singal R. DNA methylation and apoptosis. Mol Immunol. (2006) 43:1729–40. doi: 10.1016/j.molimm.2005.11.010
7. Millington GW. Epigenetics and dermatological disease. Pharmacogenomics. (2008) 9:1835–50. doi: 10.2217/14622416.9.12.1835
8. Liang Y, Chang C, Lu Q. The genetics and epigenetics of atopic dermatitis-filaggrin and other polymorphisms. Clin Rev Allergy Immunol. (2016) 51:315–28. doi: 10.1007/s12016-015-8508-5
9. Chandra A, Senapati S, Roy S, Chatterjee G, Chatterjee R. Epigenome-wide DNA methylation regulates cardinal pathological features of psoriasis. Clin Epigenetics. (2018) 10:108. doi: 10.1186/s13148-018-0541-9
10. Sang Y, Deng Y. Current insights into the epigenetic mechanisms of skin cancer. Dermatol Ther. (2019) 32:e12964. doi: 10.1111/dth.12964
11. Zhang X, Hu M, Lyu X, Li C, Thannickal VJ, Sanders YY. DNA methylation regulated gene expression in organ fibrosis. Biochim Biophys Acta Mol Basis Dis. (2017) 1863:2389–2397. doi: 10.1016/j.bbadis.2017.05.010
12. Yung PYK, Elsasser SJ. Evolution of epigenetic chromatin states. Curr Opin Chem Biol. (2017) 41:36–42. doi: 10.1016/j.cbpa.2017.10.001
13. Zhao M, Gao F, Wu X, Tang J, Lu Q. Abnormal DNA methylation in peripheral blood mononuclear cells from patients with vitiligo. Br J Dermatol. (2010) 163:736–42. doi: 10.1111/j.1365-2133.2010.09919.x
14. Deng Q, Luo Y, Chang C, Wu H, Ding Y, Xiao R. The emerging epigenetic role of CD8+T cells in autoimmune diseases: a systematic review. Front Immunol. (2019) 10:856. doi: 10.3389/fimmu.2019.00856
15. Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CYJ, Williamson NA, et al. FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics. (2015) 15:2597–601. doi: 10.1002/pmic.201400515
16. Martins LT, Frigeri HR, de Castro CCS, Mira MT. Association study between vitiligo and autoimmune-related genes CYP27B1, REL, TNFAIP3, IL2 and IL21. Exp Dermatol. (2020) 29:535–8. doi: 10.1111/exd.14100
17. Fraczek A, Owczarczyk-Saczonek A, Placek W. The role of TRM cells in the pathogenesis of vitiligo-a review of the current state-of-the-art. Int J Mol Sci. (2020) 21:3552. doi: 10.3390/ijms21103552
18. Gholijani N, Yazdani MR, Dastgheib L. Predominant role of innate pro-inflammatory cytokines in vitiligo disease. Arch Dermatol Res. (2020) 312:123–294. doi: 10.1007/s00403-019-01996-9
19. Frisoli ML, Essien K, Harris JE. Vitiligo: mechanisms of pathogenesis and treatment. Annu Rev Immunol. (2020) 38:621–48. doi: 10.1146/annurev-immunol-100919-023531
20. Speeckaert R, Dugardin J, Lambert J, Lapeere H, Verhaeghe E, Speeckaert MM, et al. Critical appraisal of the oxidative stress pathway in vitiligo: a systematic review and meta-analysis. J Eur Acad Dermatol Venereol. (2018) 32:1089–98. doi: 10.1111/jdv.14792
21. Qiao Z, Wang X, Xiang L, Zhang C. Dysfunction of autophagy: a possible mechanism involved in the pathogenesis of vitiligo by breaking the redox balance of melanocytes. Oxid Med Cell Longev. (2016) 2016:3401570. doi: 10.1155/2016/3401570
22. Yang L, Yang F, Teng L, Katayama I. 6-Shogaol protects human melanocytes against oxidative stress through activation of the Nrf2-antioxidant response element signaling pathway. Int J Mol Sci. (2020) 21:3537. doi: 10.3390/ijms21103537
23. Passeron T, Ortonne JP. Physiopathology and genetics of vitiligo. J Autoimmun. (2005) 25:63–68. doi: 10.1016/j.jaut.2005.10.001
24. Liu XM, Zhou Q, Xu SZ, Wakamatsu K, Lei TC. Maintenance of immune hyporesponsiveness to melanosomal proteins by DHICA-mediated antioxidation: possible implications for autoimmune vitiligo. Free Radic Biol Med. (2011) 50:1177–85. doi: 10.1016/j.freeradbiomed.2011.01.017
25. Tang L, Li J, Lin X, Wu W, Kang K, Fu W. Oxidation levels differentially impact melanocytes: low versus high concentration of hydrogen peroxide promotes melanin synthesis and melanosome transfer. Dermatology. (2012) 224:145–53. doi: 10.1159/000336777
26. Tofani LB, Depieri LV, Campos PM, Riul TB, Antonietto KS, de Abreu Fantini MC, et al. In vitro TyRP-1 knockdown based on siRNA carried by liquid crystalline nanodispersions: an alternative approach for topical treatment of vitiligo. Pharm Res. (2018) 35:104. doi: 10.1007/s11095-017-2330-0
27. Pyo JJ, Ahn S, Jin SH, An S, Lee E, Choi J, et al. Keratinocyte-derived IL-36gamma plays a role in hydroquinone-induced chemical leukoderma through inhibition of melanogenesis in human epidermal melanocytes. Arch Toxicol. (2019) 93:2307–20. doi: 10.1007/s00204-019-02506-6
28. Yang CA, Chiang B-L. Inflammasomes and human autoimmunity: a comprehensive review. J Autoimmun. (2015) 61:1–8. doi: 10.1016/j.jaut.2015.05.001
29. Xuehan Wu, Weilu Huang, Minxue Quan, Yongqi Chen, Jiaxin Tu, Jialu Zhou, et al. Inhibition of brain-type glycogen phosphorylase ameliorates high glucose-induced cardiomyocyte apoptosis via Akt-HIF-1α activation. Biochem Cell Biol. (2020) 6:458–65. doi: 10.1139/bcb-2019-0247
30. Sasaki S, Lesoon-Wood LA, Dey A, Kuwata T, Weintraub BD, Humphrey G, et al. Ligand-induced recruitment of a histone deacetylase in the negative-feedback regulation of the thyrotropin beta gene. EMBO J. (1999) 18:5389–98. doi: 10.1093/emboj/18.19.5389
31. Bergman RN, Iyer MS. Indirect regulation of endogenous glucose production by insulin: the single gateway hypothesis revisited. Diabetes. (2017) 66:1742–7. doi: 10.2337/db16-1320
32. Narita T, Oiso N, Fukai K, Kabashima K, Kawada A, Suzuki T. Generalized vitiligo and associated autoimmune diseases in Japanese patients and their families. Allergol Int. (2011) 60:505–8. doi: 10.2332/allergolint.11-OA-0303
33. Ovejero-Benito MC, Reolid A, Sanchez-Jimenez P, Saiz-Rodriguez M, Munoz-Aceituno E, Llamas-Velasco M, et al. Histone modifications associated with biological drug response in moderate-to-severe psoriasis. Exp Dermatol. (2018) 27:1361–71. doi: 10.1111/exd.13790
34. Gu X, Nylander E, Coates PJ, Fahraeus R, Nylander K. Correlation between reversal of DNA methylation and clinical symptoms in psoriatic epidermis following narrow-band UVB phototherapy. J Invest Dermatol. (2015) 135:2077–83. doi: 10.1038/jid.2015.128
35. Prevo R, Pirovano G, Puliyadi R, Herbert KJ, Rodriguez-Berriguete G, O'Docherty A, et al. CDK1 inhibition sensitizes normal cells to DNA damage in a cell cycle dependent manner. Cell Cycle. (2018) 17:1513–23. doi: 10.1080/15384101.2018.1491236
36. Li S, Zeng H, Huang J, Lu J, Chen J, Zhou Y, et al. Identification the competing endogenous RNA networks in oxidative stress injury of melanocytes. DNA Cell Biol. (2021) 40:192–208. doi: 10.1089/dna.2020.5455
37. Zuehlke AD, Beebe K, Neckers L, Prince T. Regulation and function of the human HSP90AA1 gene. Gene. (2015) 570:8–16. doi: 10.1016/j.gene.2015.06.018
38. Duggal S, Jailkhani N, Midha MK, Agrawal N, Rao KVS, Kumar A. Defining the Akt1 interactome and its role in regulating the cell cycle. Sci Rep. (2018) 8:1303. doi: 10.1038/s41598-018-19689-0
39. Warren CFA, Wong-Brown MW, Bowden NA. BCL-2 family isoforms in apoptosis and cancer. Cell Death Dis. (2019) 10:177. doi: 10.1038/s41419-019-1407-6
40. Guo J, Dai X, Laurent B, Zheng N, Gan W, Zhang J, et al. AKT methylation by SETDB1 promotes AKT kinase activity and oncogenic functions. Nat Cell 323 Biol. (2019) 21:226–37. doi: 10.1038/s41556-018-0261-6
41. Guan JS, Haggarty SJ, Giacometti E, Dannenberg JH, Joseph N, Gao J, et al. HDAC2 negatively regulates memory formation and synaptic plasticity. Nature. (2009) 459:55–60. doi: 10.1038/nature07925
42. Capalbo L, Bassi ZI, Geymonat M, Todesca S, Copoiu L, Enright AJ, et al. The midbody interactome reveals unexpected roles for PP1 phosphatases in cytokinesis. Nat Commun. (2019) 10:4513. doi: 10.1038/s41467-019-12507-9
43. Jung HJ, Byun HO, Jee BA, Min S, Jeoun UW, Lee YK, et al. The ubiquitin-like with PHD and ring finger domains 1 (UHRF1)/DNA methyltransferase 1 (DNMT1) axis is a primary regulator of cell senescence. J Biol Chem. (2017) 292:3729–39. doi: 10.1074/jbc.M116.750539
44. Harris JE. Cellular stress and innate inflammation in organ-specific autoimmunity: lessons learned from vitiligo. Immunol Rev. (2016) 269:11–25. doi: 10.1111/imr.12369
Keywords: vitiligo melanocyte, 850K, MDEGs, PPI, functional epigenetic modules
Citation: Pu Y, Chen X, Chen Y, Zhang L, Chen J, Zhang Y, Shao X and Chen J (2021) Transcriptome and Differential Methylation Integration Analysis Identified Important Differential Methylation Annotation Genes and Functional Epigenetic Modules Related to Vitiligo. Front. Immunol. 12:587440. doi: 10.3389/fimmu.2021.587440
Received: 26 July 2020; Accepted: 19 February 2021;
Published: 10 March 2021.
Edited by:
Reinhart Speeckaert, Ghent University Hospital, BelgiumReviewed by:
Mitesh Dwivedi, Uka Tarsadia University, IndiaStanca Birlea, University of Colorado, United States
Copyright © 2021 Pu, Chen, Chen, Zhang, Chen, Zhang, Shao and Chen. 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: Jin Chen, Y2hlbmppbjc3MTAxMCYjeDAwMDQwOzE2My5jb20=
†These authors have contributed equally to this work and share first authorship