The role of NPY2R/NFATc1/DYRK1A regulatory axis in sebaceous glands for sebum synthesis
Cellular & Molecular Biology Letters volume 28, Article number: 60 (2023)
Sebaceous glands (SGs) synthesize and secret sebum to protect and moisturize the dermal system via the complicated endocrine modulation. Dysfunction of SG are usually implicated in a number of dermal and inflammatory diseases. However, the molecular mechanism behind the differentiation, development and proliferation of SGs is far away to fully understand.
Herein, the rat volar and mammary tissues with abundant SGs from female SD rats with (post-natal day (PND)-35) and without puberty onset (PND-25) were arrested, and conducted RNA sequencing. The protein complex of Neuropeptide Y receptor Y2 (NPY2R)/NPY5R/Nuclear factor of activated T cells 1 (NFATc1) was performed by immunoprecipitation, mass spectrum and gel filtration. Genome-wide occupancy of NFATc1 was measured by chromatin immunoprecipitation sequencing. Target proteins’ expression and localization was detected by western blot and immunofluorescence.
NPY2R gene was significantly up-regulated in volar and mammary SGs of PND-25. A special protein complex of NPY2R/NPY5R/NFATc1 in PND-25. NFATc1 was dephosphorylated and activated, then localized into nucleus to exert as a transcription factor in volar SGs of PND-35. NFATc1 was especially binding at enhancer regions to facilitate the distal SG and sebum related genes’ transcription. Dual specificity tyrosine phosphorylation regulated kinase 1A (DYRK1A) contributed to NFATc1 phosphorylation in PND-25, and inactivated of DYRK1A resulted in NFATc1 dephosphorylation and nuclear localization in PND-35.
Our findings unmask the new role of NPY2R/NFATc1/DYRK1A in pubertal SG, and are of benefit to advanced understanding the molecular mechanism of SGs’ function after puberty, and provide some theoretical basis for the treatment of acne vulgaris from the perspective of hormone regulation.
As the largest organ, skin serves as the body’s natural protective barrier, which is composed of the epidermal cells as well as the epidermal surface lipid film and natural moisturizing factors that populate them. Epidermal surface lipid film is a kind of transparent lipid “membrane” with weak acidity mixed by the emulsification of sebum and sweat . Sebum is synthesized by sebaceous glands (SGs), epidermal keratinocytes and sweat glands via a unique holocrine manner, in which sebocytes undergo the differentiation, maturation, transformation into sebum vacuoles, programmed death and disintegration followed by the release of intracellular substances . Different from lipid substances (triglycerides, free fatty acid and cholesterol) derived from the other sources, sebum produced by SGs contains special squalene, wax ester, galactose, vitamins and some antimicrobial peptides, but extremely low cholesterol. Sebum also serves as a barrier to retard the release of skin odour signals produced by sweat glands and SGs [3, 4]. In a nutshell, sebum plays a crucial role in protecting and moisturizing the dermal system.
Due to the major distribution in scalp and central face (400–900 glands/cm2) compared to the average density 100/cm2 throughout the entire skin system [5, 6], dysfunction of SGs and irregular sebum secretion are usually implicated in chronic inflammatory diseases, including particular facial acne, but also atopic dermatitis, seborrhea, psoriasis and rosacea. Previous observations in clinical and experimental dermatology have highlighted the importance of nervous and neuroendocrine system majorly three neuroendocrine axes including hypothalamic-pituitary-thyroid (HPT) axis, hypothalamic-pituitary-adrenocortical (HPA) axis and hypothalamic-pituitary-gonad (HPG) axis that are implicated in modulation of cutaneous physiology and pathology. Numerous hormones, peptides and transmitters from pituitary and hypothalamus communicate with skin and its appendages through a large number of receptors [7, 8]. Beyond that, various neuropeptides and sterol hormones and cytokines from peripheral nerve fiber, capillaries, lymphatic capillaries and immune cells are also involved in the epidermal microenvironment as well as affect the proliferation and differentiation of sebocytes in burn or wound healing and tissue regeneration [9,10,11]. Nevertheless, we still have a long way to understand the molecular mechanism behind the differentiation, development and proliferation of sebaceous glands systematically and in depth.
We are more concerned about the abnormal function of SGs in adolescent acne. During pubertal onset, a number of hormones and neuropeptides are changing [12, 13]. In addition to the known distal modulations of androgen, growth factor, prolactin, adrenocorticotropin and melanotropin , we hope to seek more novel proximal regulators unmasking the substantial effect on SGs. In present research, we isolated the sebaceous gland tissue from rats with (post-natal day (PND)-35) and without puberty onset (PND-25) as previously described [14, 15] and conducted transcriptomic sequencing, and focused on neuropeptide Y receptor Y2 (NPY2R) gene. Then we further studied the subsequent intracellular and endonuclear signaling pathway for final gene expression regulation in SGs. Our findings are of benefit to advanced understanding the molecular mechanism of SGs’ function after puberty, and provide some theoretical basis for the treatment of acne vulgaris from the perspective of hormone regulation.
Total twenty female Sprague Dawley (SD) rats with post-natal day (PND)-25 (n = 10) and PND-35 (n = 10) were housed at 25 °C with enough food and water. Ten mg/kg NPY2R antagonist SF-11 (sc-311535, Santa Cruz Biotechnology)  was subcutaneously administrated to PND-35 (n = 5) in designated position for 24 h. One mg/kg DYRK1-IN-1 (HY-132308, MedChemExpress)  was subcutaneously administrated with PND-25 (n = 5) in designated position for 24 h. After sacrificing by cervical dislocation, skin tissues of volar, nipple and abdominal region were disinfected with alcohol (abdominal needs to be removed before disinfection). The whole skin was separated from the underlying tissue, and cut into 1 cm2, and stored in − 80 °C for sequencing and Oil red O staining, or in 4% paraformaldehyde for morphological detection. Three rats were used for all morphological experiments such as H&E, Oil red O, immunofluorescence assay. Two rats were used for molecular study including western blot, gel filtration, immunoprecipitation and sequencing. All animal experiments were performed in accordance with the Declaration of Helsinki and “Guide for The Care and Use of Laboratory Animals” 8th Edition, and approved by the Committee on the Ethics of Suzhou Ninth People’s Hospital (Assigned number: SZJY-2022-08).
Tissues were processed paraffin embedded, and sectioned 4 μm on a microtome. The sections were dewaxed through xylene—gradient alcohol, and stained using Hematoxylin and Eosin Kit (Beyotime Biotechnology). After that, sections were dehydrated through gradient alcohol–xylene.
Oil red O staining
Frozen skin tissues were embedded by OCT (SAKURA), and sectioned 4 μm directly on a freezing microtome. Sections were washed by 50% alcohol and added improved oil Red O dye solution (Beyotime Biotechnology) to cover the whole tissue for 20 min, then washed by washing buffer for 30 s and PBS for 20 s. After that, sections were further processed by Hematoxylin and Eosin Kit without dehydration.
The steps prior to dewaxing were consistent with H&E staining. Sections were processed by 3% hydrogen peroxide for 3 min, and rinsed with running water. After 1 min antigen repair using Tris–EDTA Antigen Repair Solution (PR30002, Proteintech), sections were washed by PBS for three times, and blocked by 5% BSA in room temperature for 30 min, then the primary antibodies including KRT5 (1: 200, #66727-1-lg, Proteintech), KRT10 (1: 200, #18343-1-AP, Proteintech), NPY1R (1: 200, ab91262, Abcam), NPY2R (1: 200,SAB4502029, MilliporeSigma), NPY5R (1: 200, ab133757, Abcam), p-NFATc1 (Ser294) (1: 200, AF8012, Affinity Biosciences), NFATc1 (1: 200, #66963-1-lg, Proteintech) or DYRK1A (1: 200, A0595, Abclonal) were added to incubate at 4 °C overnight. Appropriate secondary antibodies including FITC-conjugated affinipure goat anti-mouse IgG (H + L) (1: 10000, SA00003-1, Proteintech), FITC-conjugated affinipure goat anti-rabbit IgG (H + L) (1: 10000, SA00003-2, Proteintech), CoraLite 594-conjugated Donkey anti-mouse IgG (H + L) (1: 20000, SA00013-7, Proteintech) and CoraLite 594-conjugated Donkey anti-rabbit IgG (H + L) (1: 20000, SA00013-8, Proteintech) were incubated at room temperature avoiding light for 1 h. After PBS washing three times and drying in air, sections were added by mounting medium with DAPI-aqueous, fluoroshield (ab104139, Abcam) and covered by cover glass.
The main procedures of this experiment included RNA extraction, library construction, quality control, sequencing and analysis as described in our previous studies [15, 18]. The clean reads were aligned to rat genome assembly Rnor_6.0.
General western blotting was described previously [19, 20]. The information, concentration and treatment condition of antibodies were listed as follows (antibodies were sorted according to the order in which they appeared), NPY1R (1: 2000, ab91262, Abcam), NPY2R (1: 2000, SAB4502029, MilliporeSigma), NPY5R (1: 2000, ab133757, Abcam), GAPDH (1: 5000, AC001, Abclonal), NFATc1 (1: 2000, #66963-1-lg, Proteintech), p-NFATc1 (Ser294) (1: 1000, AF8012, Affinity Biosciences), CDC2 (1: 2000, #9116, CST), p-CDC2 (Tyr15) (1: 1000, #9111, CST), PIM-1 (1: 2000, #54523, CST), p-PIM-1 (Tyr218) (1: 500, orb6723, Biorbyt), mTORC1 (1: 2000, A2445, Abclonal), p-mTORC1 (Ser2448) (1: 1000, AP0094, Abclonal), IKKε (1: 2000, A16470, Abclonal), p-IKKε (Ser172) (1: 1000, #8766, CST), MAP3K8 (1: 2000, A15623, Abclonal), p-MAP3K8 (Ser400) (1: 1000, #4491, CST), DYRK1A (1: 2000, A0595, Abclonal), p-DYRK1A (Tyr271/319) (1: 1000, AF3507, Affinity Biosciences), HRP-conjugated affinipure goat anti-rabbit IgG (H + L) (1: 5000, SA00001-2, Proteintech) and HRP-conjugated affinipure goat anti-mouse IgG (H + L) (1: 5000, SA00001-1, Proteintech).
The gel filtration assay was performed as previously described . In brief, hypotonic buffer with NP-40 was used to separate the cytoplasm and nuclei of tissues containing 1 × 108 cells and the cytoplasmic proteins were isolated using low salt extraction buffer [20 mM HEPES (pH 7.9), 50 mM NaCl, 25% glycerol, 1.5 mM MgCl2, 0.2 mM EDTA, 0.5 mM dithiothreitol, and protease inhibitors]. The cytoplasmic extracts (4 mg) were directly applied to a Sepharose 6B column (MilliporeSigma) equilibrated with column running buffer containing 20 mM HEPES (pH 7.9), 200 mM NaCl, 1 mM dithiothreitol, 0.1 mM phenylmethylsulfonyl fluoride, and 10% glycerol. Fractions of 1 mL each were collected and NPY1R, NPY2R, NPY5R and NFATc1 were detected by western blot.
Frozen tissues were pulverized using ceramic grinding rod in liquid nitrogen, and fixed by 1 mL 4% paraformaldehyde for 10 min, then added glycine to 0.2 M for quenching fixation. After 14,000 rpm centrifugation for 10 min as well as supernatant removal, cell precipitation was added 1 mL cold standard RIPA buffer (#20-188, Millipore) and incubated on ice for 30 min. After 14,000 rpm centrifugation again, supernatant was transferred into another new eppendorf tube. Ultrasonication optional for protein Immunoprecipitation, but necessary for chromatin immunoprecipitation was set as high power 5 s, break 30 s for 35 cycles by CHIP-300 T System (Jiangsu Gaozhou Technology). 20 μL was stored as input, and could be subsequently extracted the total protein or genomic DNA. The rest supernatant was incubated with 1 μg appropriated primary antibodies including NPY5R (ab133757, Abcam) or NFATc1 (MA3-024, Thermo Fisher Scientific) at 4 °C overnight. In next morning, 20 μL Magna ChIP Protein A + G Magnetic Beads (#16-663, Millipore) were added and incubated for another 2 h. Beads were washed by modified washing buffer I [RIPA buffer plus 500 mM LiCl for ChIP, plus 350 mM NaCl for protein IP] for three times, then washed by washing buffer II [RIPA buffer plus 150 mM NaCl] for three times, and finally added 20 μL RIPA buffer. For ChIP, another 480 μL RIPA buffer was added, and added 500 μL Tris-phenol for DNA isolation by general phenol chloroform. Library construction, quality control, sequencing and analysis were described in our previous studies . For protein IP, 5 μL 5 × SDS-PAGE Protein Loading Buffer (Sangon Biotech) was added, and run 10% SDS-PAGE gel. Coomassie brilliant blue staining for interested protein band by mass spectrum (OE Biotech), or western blot for verification of protein–protein interaction.
Data are presented as means ± standard deviations. Student’s t-test was used to analyze the intergroup differences. The p value less than 0.05 was considered as statistical significance.
Sebaceous glands became vibrant after puberty onset
Initially, skin tissue of volar, nipple and abdominal regions from juvenile (PND-25) and pubertal (PND-35) female SD rats were obtained. In general, SGs were mainly enriched in volar and mammary skin. We observed that SGs were gigantic and numerous especially in volar skin in PND-35 compared to PND-25 (Fig. 1A). Moreover, robust expression of KRT5 (biomarker of progenitor cells) and KRT10 (biomarker of ductal and sebocytic cells)  in PND-35 compared to PND-25 (Fig. 1B). Consistently, Oil red O staining was conducted to investigate the status of sebum synthesis and secretion, and sebum levels in tissues with volar and mammary skin were higher in PND-35 compared to PND-25 (Fig. 1C). The results above determined that differentiation of SGs and sebum synthesis and secretion initiated were substantially activated in pubertal rats.
NPY2R connected with the suppressive effects on sebaceous glands
Next, we studied the transcriptome changes of these skin tissues by RNA sequencing. Differential expressed genes (DEGs) compared between PND-25 and PND-35 were obtained (Fig. 2A, Additional file 1: Table S1). Due to the primary concern on DEGs for SG, abdominal skin samples without SGs was excluded, and the remaining subset of 362 genes (127 up-regulated and 235 down-regulated in PND-25 vs PND-35) were considered as the putative special DEGs for SGs (Fig. 2B). “Metascape” showed that GO and KEGG analysis on these 362 DEGs showed that hormone regulation, gland development, cell junction and adhesion as well as epithelial differentiation were affected along with puberty onset (Fig. 2C). Here, we observed that neuropeptide Y receptor Y2 (Npy2r) was significantly down-regulated in PND-35 (PND-25 vs PND-35: log2 FC = 1.62, p = 0.018 in volar skin; log2 FC = 1.73, p = 0.005 in mammary skin) (Additional file 1: Table S1), implying that the corresponding substrate NPY might exert as an important ligand for the function of SG. Consistently, only the protein expression of NPY2R in NPY receptor family (NPY1R, NPY2R and NPY5R) was reduced both in volar and mammary skin tissues of PND-35 compared to PND-25 (Fig. 2D, E). Interestingly, NPY2R antagonist SF-11 was subcutaneously administrated to PND-35 rat for 24 h, and resulted in the reduced expression of KRT5 and KRT10, as well as compromised the sebum level (Fig. 1B, C). All these data above suggested that NPY2R was supposed to be attributed to sebum synthesis.
Heterodimer of NPY2R/NPY5R interacted with phosphorylated NFATc1
Due to the dimer form of NPY receptor on cytomembrane , we studied the potential heterodimer of NPY2R with NPY1R or NPY5R using gel filtration assay. Our findings indicated that NPY1R was located in different protein fraction from NPY2R, and was likely to form multiple complex with other unknown proteins in vola of PND-25 because only multiple peaks superimposed would show the sustained high expression in a series of continuous fraction (Fig. 3A). In turn, NPY2R and NPY5R clearly paired up around Fraction 23 and 35 (Fig. 3A). By contrast, due to the reduced expression of NPY2R in PND-35, NPY5R still appeared exactly in the same Fraction 35 as PND-25, but completely lost in Fraction 25, whereas the distribution of NPY1R was not affected (Fig. 3B). The discovery above suggested that one complex containing heterodimer of NPY2R/NPY5R was deteriorated in PND-35 owing to the reduction of NPY2R. To figure out this variable complex and the underlying mechanism behind, NPY5R was captured via immunoprecipitation by antibody. We noticed an interested band around 140 kDa, which obviously weakened in PND-35 compared to PND-25 (Fig. 3C), and identified it as nuclear factor of activated T cells 1 (Nfatc1) by mass spectrum (Fig. 3D). The interaction between NFATc1 and NPY2R/NPY5R was validated by gel filtration (Fig. 3A) and immunoprecipitation (Fig. 3E). In PND-35, we observed that NFATc1 was dispersed and the total protein was substantially reduced, suggesting that NFATc1 was more likely to exist in the free form (Fig. 3B). As previously reported, the function of NFATc1 was closely related to its phosphorylation [24, 25]. We further observed that p-NFATc1 at Ser294 interacting with NPY2R was abundant in PND-25 but almost deficient in PND-35 (Fig. 3A, B, E). Immunofluorescence also uncovered a remarkable reduction of phosphorylated NFATc1 in PND-35 (Fig. 3F). Taken together, we concluded that NFATc1 could substantially interact with NPY2R in sebocytes in PND-25.
Genome-wide occupancy of NFATc1 enhanced transcription activities of sebocytes
NFATc1 enters the nucleus and acts as a transcription factor after dephosphorylation. Therefore, we further studied the genome-wide occupancy of NFATc1 in sebocytes. Our observation indicated that the abundance of NFATc1 on genome was more robust in PND-35 than PND-25 (Fig. 4A). Moreover, we noticed that NFATc1 binding at the intergenic regions was increased in PND-35 compared PND-25 (Fig. 4B). Compared with known enhancer elements of rat species using EnhancerAtlas Database (http://www.enhanceratlas.org/) , it was notable that there were 154 enhancer regions overlapping with these significantly increased binding peaks of NFATc1 compared between PND-35 and PND-25 (log2FC > 1, p < 0.05). Simultaneously, 126 of these 154 enhancers contained at least three DEGs within 10 kb range (Fig. 4C). IGV showed two canonical enhancers (Chr1: 168825516–168927272, chr20: 5798616–5800552, Rnor5) followed by a series of genes involved in sebaceous gland synthesis Ucp2 , and Pou5f1  (Fig. 4D, E). On the other hand, no significant correlation or weak correlation was found between differential gene expression and differential NFATc1 binding between PND-25 and PND-35 on other regions except enhancer (Fig. 4F). These data above suggested that the differentially binding of NFATc1 on enhancers could indeed govern the target genes’ transcription in the system of sebaceous glands.
Inhibition of DYRK1A promoted the nuclear translocation and activity of NFATc1
Next, we focused on which phosphatase or phosphokinase was attributed to NFATc1 dephosphorylation in sebaceous glands of PND-35. Kinases such as CDC2 , PIM-1 , mTORC1  and IκB Kinase ε , MAP3K8  as well as DYRK1A [33, 34] have been reported to affect the phosphorylation and nuclear translocation of NFATc1 in other tissues or diseases. Expression of these kinases were detected by western blot. We observed that all these proteins appeared to be equal between PND-25 and PND-35 (Fig. 5A). Strikingly, the phosphorylation of DYRK1A (Tyr271/319) was significantly reduced, while mTOR (Ser2448) was slightly elevated in PND-35 compared to PND-25 (Fig. 5A). The activation of mTOR apparently had the opposite effect on nuclear translocation of NFATc1. Therefore, we only examined the effect of DYRK1A using DYRK1A inhibitor to the subcutaneous vola of PND-25, and observed that the phosphorylation of NFATc1 was also compromised (Fig. 5B). Meanwhile, the staining of NFATc1 was abundant in nucleus accompanied with reduced localization of DYRK1A in PND-25 induced by DYRK1A inhibitor or PND-35 (Fig. 5C). Taken together, our findings suggested that DYRK1A was responsible for the nuclear translocation and activity of NFATc1 in sebocytes.
Relevant epidemiological studies have figured out that hormone, inflammation, allergies, stress, irritation, cosmetics, blood lipids, dairy food intake, smoking, hyperpigmentation, ethnicity, and genetics are implicated in acne . Nevertheless, hormone imbalances may render acne unresponsive to conventional therapy, suggesting that endocrine regulation plays a vital role in acnegenesis. Moderate-to-severe acne affects around 20% of young people and severity correlates with pubertal maturity . Exception of the androgen-stimulated sebum production, other hormones, including estrogens, growth hormone, insulin, insulin-like growth factor-1, glucocorticoids, adrenocorticotropic hormone, and melanocortins all participate in the pathophysiology of acne . Corresponding to them, sebum cells express a variety of receptors, including peptide hormone receptors, neurotransmitter receptors, and steroid and thyroid hormone receptors. Various hormones and mediators acting through sebum cell receptors finally determine the choreography of acne pathologies.
NPY has been shown to be the most widely distributed neuropeptide throughout multiple tissues in a variety of organisms, and has pleiotropic roles throughout the body in the central biological processes . NPY is elevated in the affected skin and/or circulation of humans suffering from various skin pathologies, including atopic dermatitis, cutaneous melanoma, psoriasis, and vitiligo . The biological behavior of NPY exerts through the interaction with five, G protein-coupled, membrane-bound, NPY Y receptors (Y1R, Y2R, Y4R, Y5R) . However, the molecular role of NPY and receptors in sebaceous regulation and acne is still poorly understood to date. NPY4R is not expressed in all skin samples by RNA sequencing, therefore we majorly focus on the other three receptors in current study.
Now our observations provide novel molecular evidence of NPY and NPY2R in SGs. RNA sequencing has indicated that the transcript of Npy2r is reduced in volar and mammary tissue with abundant SGs after pubertal initiation (Additional file 1: Table S1). Furthermore, NPY2R inhibition and compensation can indeed affect the growth of sebaceous progenitors and sebum production (Fig. 1B, C). NPY signaling occurrs via the NPY2 receptor (NPY2R), stimulating PI3K, MAPK, mTORC1 and NFAT activation . NPY-NPY2R signaling is benefit to vasodilation, and cell proliferation, but not good for cell differentiation and migration [42,43,44]. Therefore, the functional activation of SGs as well as the reduced Npy2r after puberty onset are logically self-consistent.
Moreover, the biological significance of the differential expression of Npy2r lies in the NPYR heterodimer. Different combinations of NPYR heterodimer showing the tissue and disease specificity are usually subject to the subcellular localization of NPYR , NPY concentrations and target proteins coordinated with the receptor-ligand interactions . Our findings suggest that dimerization of NPY2R and NPY5R is specific in system of SG, and loss of NPY2R substantially deteriorate the interaction with NPY5R. Beyond that, we have unexpectedly discovered that NPY5R appear to participate in two kinds of complex, the larger of which independent of NPY2R (Fraction 35) is likely a homodimer, and is worthy of further investigation into its function on NPY signal conduction.
Furthermore, subsequent observations of IP and mass spectrum indicate that NFATc1 is a time-specific regulator interacting with NPY2R/NPY5R. NFATc1 is calcium–calcineurin dependent regulator that contributes to osteoclastogenesis , cardiac valves formation , and regulatory T-cells suppression . In classical osteoclastogenesis activated by RANKL, PLCγ binds to the IP3 receptor (IP3R), induces calcium release from the endoplasmic reticulum, and then increases calcium levels to activate calcineurin, which is followed by the dephosphorylation and nuclear translocation of NFATc1 . Although the issue about calcium–calcineurin is not covered in this paper, we still believe that NFATc1 is dependent on their regulation whatever it is inhibited in SGs in the early puberty or activated after puberty onset. The relationship between calcium and the action of DYRK1A or other phosphokinases is worthy of further exploration in future studies. Moreover, as a master transcription factor, NFATc1 can bind to many enhancers for remote transcriptional regulation of terrestrial sebaceous gland related genes in our system. But this does not mean that NFATc1 has similar genome-wide occupancy in other disease states and different tissue types. For example, we compared the genomic binding data of NFATc1 in one article on pancreatic cancer , and the results were completely different (Data now shown). We have reason to believe that NFATc1 has a tissue-specific regulatory mode in the nucleus, showing more unknown epigenetic regulatory factors interacting with it, which is a new research direction developed from NPY2R.
In summary, we figure out a novel molecular mechanism of NPY2R/NFATc1/DYRK1A regulatory axis in sebum synthesis. Our findings establish a new regulatory relationship in the field of sebaceous gland secretion, and have a certain theoretical value for the advancement of the research on skin diseases, skin repair, medical cosmetology, pharmaceuticals and cosmetics.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional file.
Neuropeptide Y receptor
Differential expressed gene
Nuclear factor of activated T cells 1
Cell division cycle 2
Pim-1 proto-oncogene, serine/threonine kinase
Mammalian target of rapamycin C1
Dual specificity tyrosine phosphorylation regulated kinase 1A
Integrative genomics viewer
Wang T, Zhou Z, Luo E, Zhong J, Zhao D, Dong H, Yao B. Comprehensive RNA sequencing in primary murine keratinocytes and fibroblasts identifies novel biomarkers and provides potential therapeutic targets for skin-related diseases. Cell Mol Biol Lett. 2021;26:42.
Zouboulis CC, Picardo M, Ju Q, Kurokawa I, Torocsik D, Biro T, Schneider MR. Beyond acne: current aspects of sebaceous gland biology and function. Rev Endocr Metab Disord. 2016;17:319–34.
Mollebjerg A, Palmen LG, Gori K, Meyer RL. The bacterial life cycle in textiles is governed by fiber hydrophobicity. Microbiol Spectr. 2021;9: e0118521.
Kippenberger S, Havlicek J, Bernd A, Thaci D, Kaufmann R, Meissner M. “Nosing Around” the human skin: what information is concealed in skin odour? Exp Dermatol. 2012;21:655–9.
Ludovici M, Kozul N, Materazzi S, Risoluti R, Picardo M, Camera E. Influence of the sebaceous gland density on the stratum corneum lipidome. Sci Rep. 2018;8:11500.
Abramovits W, Gonzalez-Serva A. Sebum, cosmetics, and skin care. Dermatol Clin. 2000;18:617–20.
Pondeljak N, Lugovic-Mihic L. Stress-induced interaction of skin immune cells, hormones, and neurotransmitters. Clin Ther. 2020;42:757–70.
Clayton RW, Langan EA, Ansell DM, de Vos I, Gobel K, Schneider MR, Picardo M, Lim X, van Steensel MAM, Paus R. Neuroendocrinology and neurobiology of sebaceous glands. Biol Rev Camb Philos Soc. 2020;95:592–624.
Yannas IV, Tzeranis DS, So PTC. Regeneration of injured skin and peripheral nerves requires control of wound contraction, not scar formation. Wound Repair Regen. 2017;25:177–91.
Liu B, Chen R, Wang J, Li Y, Yin C, Tai Y, Nie H, Zeng D, Fang J, Du J, et al. Exploring neuronal mechanisms involved in the scratching behavior of a mouse model of allergic contact dermatitis by transcriptomics. Cell Mol Biol Lett. 2022;27:16.
Nowak NC, Menichella DM, Miller R, Paller AS. Cutaneous innervation in impaired diabetic wound healing. Transl Res. 2021;236:87–108.
Herbison AE. Control of puberty onset and fertility by gonadotropin-releasing hormone neurons. Nat Rev Endocrinol. 2016;12:452–66.
Fudvoye J, Lopez-Rodriguez D, Franssen D, Parent AS. Endocrine disrupters and possible contribution to pubertal changes. Best Pract Res Clin Endocrinol Metab. 2019;33: 101300.
Zhou S, Zang S, Hu Y, Shen Y, Li H, Chen W, Li P, Shen Y. Transcriptome-scale spatial gene expression in rat arcuate nucleus during puberty. Cell Biosci. 2022;12:8.
Zhou S, Shen Y, Zang S, Yin X, Li P. The epigenetic role of HTR1A antagonist in facilitaing GnRH expression for pubertal initiation control. Mol Ther Nucleic Acids. 2021;25:198–206.
Domin H, Piergies N, Pieta E, Wyska E, Pochwat B, Wlaz P, Smialowska M, Paluszkiewicz C, Szewczyk B. Characterization of the brain penetrant neuropeptide Y Y2 receptor antagonist SF-11. ACS Chem Neurosci. 2019;10:3454–63.
Henderson SH, Sorrell F, Bennett J, Fedorov O, Hanley MT, Godoi PH, Ruela de Sousa R, Robinson S, Ashall-Kelly A, Hopkins Navratilova I, et al. Discovery and characterization of selective and ligand-efficient DYRK inhibitors. J Med Chem. 2021;64:11709–28.
Shen Y, Zhou S, Zhao X, Li H, Sun J. Characterization of genome-wide DNA methylation and hydroxymethylation in mouse arcuate nucleus of hypothalamus during puberty process. Front Genet. 2020;11: 626536.
Hu Y, Shen Y, Zhao Y, Tang Y, Liu C, Gu Y, Yang T, Shen Y. Genomic distribution of signal transducer and activator of transcription (STAT) family in colorectal cancer. Hum Cell. 2023;36:286–95.
Shen Y, Zhao H, Zhang L, Hu Y, Cai L, Li J, Zhou S. The roles of DNA methylation and hydroxymethylation at short interspersed nuclear elements in the hypothalamic arcuate nucleus during puberty. Mol Ther Nucleic Acids. 2021;26:242–52.
Han B, Park HK, Ching T, Panneerselvam J, Wang H, Shen Y, Zhang J, Li L, Che R, Garmire L, Fei P. Human DBR1 modulates the recycling of snRNPs to affect alternative RNA splicing and contributes to the suppression of cancer development. Oncogene. 2017;36:5382–91.
Boecker W, Reusch M, Mielke V, Reusch U, Loening T, Tiemann M, Buchwalow I. Spatial analysis of p63, K5 and K7 defines two groups of progenitor cells that differentially contribute to the maintenance of normal sebaceous glands, extraocular sebaceous carcinoma and benign sebaceous tumors. J Dermatol. 2019;46:249–58.
Czarnecka M, Lu C, Pons J, Maheswaran I, Ciborowski P, Zhang L, Cheema A, Kitlinska J. Neuropeptide Y receptor interactions regulate its mitogenic activity. Neuropeptides. 2019;73:11–24.
Eerola SK, Santio NM, Rinne S, Kouvonen P, Corthals GL, Scaravilli M, Scala G, Serra A, Greco D, Ruusuvuori P, et al. Phosphorylation of NFATC1 at PIM1 target sites is essential for its ability to promote prostate cancer cell migration and invasion. Cell Commun Signal. 2019;17:148.
Liu Z, Li H, He L, Xiang Y, Tian C, Li C, Tan P, Jing J, Tian Y, Du L, et al. Discovery of small-molecule inhibitors of the HSP90-Calcineurin-NFAT pathway against glioblastoma. Cell Chem Biol. 2019;26(352–365): e357.
Gao T, He B, Liu S, Zhu H, Tan K, Qian J. EnhancerAtlas: a resource for enhancer annotation and analysis in 105 human cell/tissue types. Bioinformatics. 2016;32:3543–51.
Mori S, Yoshizuka N, Takizawa M, Takema Y, Murase T, Tokimitsu I, Saito M. Expression of uncoupling proteins in human skin and skin-derived cells. J Invest Dermatol. 2008;128:1894–900.
Takechi T, Hirota T, Fujii K, Nakahara T, Sakai T, Maeda N, Furue M, Ieiri I. Effect of genetic polymorphisms of human SLC22A3 in the 5’-flanking region on OCT3 expression and sebum levels in human skin. J Dermatol Sci. 2021;101:4–13.
Kim HM, He L, Lee S, Park C, Kim DH, Han HJ, Han J, Hwang J, Cha-Molstad H, Lee KH, et al. Inhibition of osteoclasts differentiation by CDC2-induced NFATc1 phosphorylation. Bone. 2020;131: 115153.
Huynh H, Wan Y. mTORC1 impedes osteoclast differentiation via calcineurin and NFATc1. Commun Biol. 2018;1:29.
Zhang J, Feng H, Zhao J, Feldman ER, Chen SY, Yuan W, Huang C, Akbari O, Tibbetts SA, Feng P. IkappaB kinase epsilon Is an NFATc1 kinase that inhibits t cell immune response. Cell Rep. 2016;16:405–18.
Kuroda Y, Hisatsune C, Mizutani A, Ogawa N, Matsuo K, Mikoshiba K. Cot kinase promotes Ca2+ oscillation/calcineurin-independent osteoclastogenesis by stabilizing NFATc1 protein. Mol Cell Biol. 2012;32:2954–63.
Liu H, Wang K, Chen S, Sun Q, Zhang Y, Chen L, Sun X. NFATc1 phosphorylation by DYRK1A increases its protein stability. PLoS ONE. 2017;12: e0172985.
Lee Y, Ha J, Kim HJ, Kim YS, Chang EJ, Song WJ, Kim HH. Negative feedback Inhibition of NFATc1 by DYRK1A regulates bone homeostasis. J Biol Chem. 2009;284:33343–51.
Habeshian KA, Cohen BA. Current issues in the treatment of acne vulgaris. Pediatrics. 2020;145:S225–30.
Heng AHS, Chew FT. Systematic review of the epidemiology of acne vulgaris. Sci Rep. 2020;10:5754.
Thiboutot D. Regulation of human sebaceous glands. J Invest Dermatol. 2004;123:1–12.
Groneberg DA, Folkerts G, Peiser C, Chung KF, Fischer A. Neuropeptide Y (NPY). Pulm Pharmacol Ther. 2004;17:173–80.
Anderson ZT, Dawson AD, Slominski AT, Harris ML. Current insights into the role of neuropeptide Y in skin physiology and pathology. Front Endocrinol. 2022;13: 838434.
Lindner D, Stichel J, Beck-Sickinger AG. Molecular recognition of the NPY hormone family by their receptors. Nutrition. 2008;24:907–17.
Lay AC, Barrington AF, Hurcombe JA, Ramnath RD, Graham M, Lewis PA, Wilson MC, Heesom KJ, Butler MJ, Perrett RM, et al. A role for NPY-NPY2R signaling in albuminuric kidney disease. Proc Natl Acad Sci USA. 2020;117:15862–73.
Stephens DP, Saad AR, Bennett LA, Kosiba WA, Johnson JM. Neuropeptide Y antagonism reduces reflex cutaneous vasoconstriction in humans. Am J Physiol Heart Circ Physiol. 2004;287:H1404-1409.
Hodges GJ, Jackson DN, Mattar L, Johnson JM, Shoemaker JK. Neuropeptide Y and neurovascular control in skeletal muscle and skin. Am J Physiol Regul Integr Comp Physiol. 2009;297:R546-555.
Kang X, Ma X, Li H, Jin X, Gao X, Feng D, Wu S. Neuropeptide Y promotes mTORC1 to regulate chondrocyte proliferation and hypertrophy. Endocrinology. 2023;164: bqac213.
Kobayashi Y. New perspectives on GPCRs: GPCR heterodimer formation (melanocortin receptor) and GPCR on primary cilia (melanin concentrating hormone receptor). Gen Comp Endocrinol. 2020;293: 113474.
Gehlert DR, Schober DA, Morin M, Berglund MM. Co-expression of neuropeptide Y Y1 and Y5 receptors results in heterodimerization and altered functional properties. Biochem Pharmacol. 2007;74:1652–64.
Kang JY, Kang N, Yang YM, Hong JH, Shin DM. The role of Ca(2+)-NFATc1 signaling and its modulation on osteoclastogenesis. Int J Mol Sci. 2020;21:3646.
Wu B, Baldwin HS, Zhou B. Nfatc1 directs the endocardial progenitor cells to make heart valve primordium. Trends Cardiovasc Med. 2013;23:294–300.
Giri PS, Bharti AH, Begum R, Dwivedi M. Calcium controlled NFATc1 activation enhances suppressive capacity of regulatory T cells isolated from generalized vitiligo patients. Immunology. 2022;167:314–27.
Asagiri M, Takayanagi H. The molecular understanding of osteoclast differentiation. Bone. 2007;40:251–64.
Patil S, Forster T, Reutlinger K, Kopp W, Versemann L, Spitalieri J, Gaedcke J, Strobel P, Singh SK, Ellenrieder V, et al. Chromatin-independent interplay of NFATc1 and EZH2 in pancreatic cancer. Cells. 2021;10:3463.
We are also grateful to Shanghai GeneFund Co., Ltd. for sequencing data analysis.
This project is supported by Suzhou Medical Health Science and Technology Innovation Project (SKYD2022023, SKJY2021016), Science and Technology Planning Fund of Affiliated Hospital of Xuzhou Medical University (XYFY202219, XYFM202217), Gusu Health Talent Research Project (GSWS2022111), and Start-up Funding of Suzhou Ninth People’s Hospital, College Project of Suzhou Ninth People’s Hospital (YK202302, YK202305).
Ethics approval and consent to participate
All animal experiments were approved by the Committee on the Ethics of Suzhou Ninth People’s Hospital (Approval number: SZJY-2022-08) on February 22nd 2022. All animal experiments were performed in accordance with the Declaration of Helsinki and “Guide for The Care and Use of Laboratory Animals” 8th Edition.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Gene information of DEGs compared between PND-25 and PND-35. DEGs are defined as log2 fold change > 1 or < − 1, and p < 0.05. MS: mammary skin; VS: volar skin; AS: abdominal skin.
About this article
Cite this article
Yang, T., Hei, R., Li, X. et al. The role of NPY2R/NFATc1/DYRK1A regulatory axis in sebaceous glands for sebum synthesis. Cell Mol Biol Lett 28, 60 (2023). https://doi.org/10.1186/s11658-023-00467-4