Epithelial organs are the first barrier against microorganisms and genotoxic stress, in which the p53 family members p63 and p73 have both overlapping and distinct functions. Intriguingly, p73 displays a very specific localization to basal epithelial cells in human tissues, while p63 is expressed in both basal and differentiated cells. Investigations of p63 and p73 protein–protein interactions reveal distinct functions underlying the aforementioned distribution. The p73 and p63 cooperate in the genome stability surveillance in proliferating cells; p73 specific interactors contribute to the transcriptional repression, anaphase promoting complex and spindle assembly checkpoint, whereas p63 specific interactors play roles in the regulation of mRNA processing and splicing in both proliferating and differentiated cells suggesting diversification of the RNA and DNA specific functions within the p53 family.
The p53-family is a group of proteins, which consists of p53, p63, and p73 . p63 and p73 are paralogs of p53 and have similar structure: transactivation (TA) at the N-terminus following by DNA-binding, nuclear localization, tetramerization domains and sterile alpha motif (SAM) domain, a specific feature of the p63 and p73 at the C-terminus (Figure 1) .
p53, p63 and p73 proteins are represented by multiple isoforms that regulate each other’s transcription and functions . Full length p73 and p63 contain TA domain and alternative intragenic transcriptional start sites between p63/p73 exons 3 and 4 produce N-truncated “DN” isoforms lacking the TA domain. Full length “-alfa” isoforms contain SAM domain and alternative C-terminal splicing produces beta, gamma, etc., isoforms lacking the SAM domain. Thus, an isoform nomenclature “TAp73α” means that a p73 isoform contains N-terminal TA domain with full length SAM domain at the C-terminal. Accordingly, DNp73β and above (gamma, zeta etc.) means that the activation domain is missing and the SAM domain is truncated. Thereby, combination of N- and C-terminal variants creates a variety of isoforms (Figure 1) .
Figure 1. Schematic representation of the p73 domain structure for the main isoform expressed in the epithelial tissue. The p63 follows the same nomenclature. Altogether, 13 isoform of p63 and 13 isoforms of p73 exist. DN- isoforms begin with the stretch of 13 amino acids that replace TA domain at 62 bp. The TAp73a coordinates are for the O15350 protein and DNp73a coordinates are for the BAB87244.1 protein. TA – transactivation domain, DBD – DNA binding domain, OD - oligomerization domain, SAM – sterile alpha motif domain and TID - transcription inhibitory domain.
In the full length TAp73α and TAp63a, N-terminal TA domain acts as the transcriptional activator whereas C- terminal containing SAM and inhibitory domains repress TA activity . Generally, isoforms lacking TA domain DNp63 and DNp73 repress activity of isoforms containing TA domain . Also, p73 isoforms can act as conditional activators or repressors . For example, both TAp73β and DNp73β but neither p53, TAp73α, TAp63α, nor TAp63β induce the expression of caspase-2S . Similarly, TAp73α co-activates c-Jun whereas TAp73β can act as a repressor of c-Jun mediated transactivation .
In particular, DNA-binding domain is the most conserved structural element across the p53-family (Figure 4), leading to identical DNA recognition , underlining the functional relatedness of these proteins as transcription factors. p53 functions as a transcription factor in charge of the activation of genes responsible for the cell-cycle arrest in G1 and G2/M checkpoints and apoptosis induction. At the same time, both p73 and p63 are known for regulation of the p53-responsive genes: MDM2, IGFBP3, p21 . However, apart from the large amount of evidence that p53, p63 and p73 proteins have overlapping and interdependent functions in the apoptosis control and genome stability , p63 and p73 also have distinct specific roles in developmental processes , DNA repair  and ciliogenesis .
Specifically, p63 and p73 play distinct roles in epithelial morphogenesis. The p63 is a major player in skin biology and its role in epidermis was thoroughly investigated during the last two decades .
The p63 knockout in mice leads to E1 lethality and pups born with non-differentiated skin . Total p63 knockout has been rescued by specific p63 isoforms expression under the control of KRT5 promoter . Double TAp63a/DNp63a rescue showed regions of differentiated skin, while p63-/TAp63 did not rescue, whereas p63-/DNp63 showed only basal keratinocytes formation. Consistently, DNp63 regulated basal undifferentiated keratinocytes marker KRT14, whereas TAp63 regulated markers of differentiated keratinocytes (Ets-1, KRT1, transglutaminases, involucrin) . In addition, p63 function is required for the establishment of the cell division perpendicular to the basal membrane at E18.5 day of the embryo . Moreover, both β1-integrin and laminin are distributed around the cells and not localised to the basal membrane in the p63 knockout mice .
In contrast, the p73 knockout mice are characterised by runting phenotype, hippocampal dysgenesis and hydrocephalus, sterility, chronic inflammation and infection in the lungs, sinus, and ears . Total p73 knockout mice do not have profound skin phenotype except the somewhat slower wound-healing response . TAp73 knockout mice demonstrate sterility, hippocampal dysgenesis, hydrocephalus, premature aging, genomic instability, and increased frequency of tumors . Finally DNp73 knockout causes neurodegeneration, specifically hippocampal dysgenesis and hydrocephalus .
In the human skin, some keratinocytes fail to divide and become polyploid upon differentiation that can be induced by the DNA damage and is associated with overexpression of Cyclin E . Using fluorescence in situ hybridization for DNA, it was shown that human skin keratinocytes become multinucleated with their genomes increasing up to 4 fold .
Intriguingly, p73 displays a very specific localization to basal epithelial cells in all human tissue where it’s localization was determined by immunohistochemistry. In contrast, the p63 is expressed in both basal and differentiated cells .
The role of p63 in polyploidy has not been evaluated, however, the presence of p63 in differentiating keratinocytes suggests that p63 is not interfering with the polyploidy of these cells. At the same time, p73 suppresses polyploidy and aneuploidy in the absence of functional p53 . Also, cells from TAp73−/− mice exhibit genomic instability associated with the enhanced aneuploidy . Therefore, it is possible that the differential presence of the p63 and p73 in the differentiated and basal cells is related to the strict polyploidy control by the p73. Specifically, the absence of p73 allows cells to loosen the mitotic spindle assembly checkpoint and hence become polyploid and multinucleated.
According to the human protein and RNA expression atlas http://www.proteinatlas.org , the p53 protein could be detected at extremely low levels in the oral mucosa, oesophagus, urinary bladder, skin, and tonsils, and this expression is variable across samples and antibodies. In contrast, p63 and p73 proteins are consistently expressed. Moreover, proteins exhibit a characteristic expression pattern. In the skin, p63 levels are the highest in the basal cells, but also detected in the squamous, more differentiated cells whereas p73 is found only in the basal cells (Figure 2A,B). Similarly, in the nasopharynx (Figure 2D,E), vagina (not shown) and oral mucosa (Figure 2G,H) the p73 is restricted to the basal cells and p63 is expressed throughout the basal epithelium and differentiated cells. The p73 expression can be detected in a few cells of the basal layer of the salivary glands lumens, whereas p63 is determined in the different cells (not shown).
In the bronchial epithelium of the lungs, p63 and p73 staining patterns seem to be similar (Figure 2J–K). Whereas p63 antibody strongly stains these tissues, with the different staining patterns than described above, it also stains myoepithelial cells of the breast, urinary bladder, epididymis, seminal vesicles, prostate, and placenta. Interestingly, in these organs, p63 stained cells immediately next to the basal membranes, resembling the p73 staining pattern in the organs expressing both p63 and p73.
In vivo, mitosis mostly occurs in the basal layer, where the p73 is, but a few cells go through mitosis just above the basal layer (Figure 2C,F,I,L) marking the last asymmetric division, leading to the terminal differentiation for one cell and an ability to proliferate for another .
How do these data correlate with RNA expression levels of p63 and p73?
According to GTEx portal RNA-seq data, expression of p73 is quite tissue specific, with relatively high levels of DNp73 in epithelial tissue: skin, vagina, esophagus (similar to the protein Atlas data) and TAp73 is expressed in the brain https://www.gtexportal.org/home/gene/TP73.
Similarly to the immunohistochemistry images, the level of tp73 expression in skin is 7.8 transcripts per million (TPM), whereas the level of tp63 is 127 TPM and tp53 is 36 TPM (Supplementary Excel File). In comparison, MDM2 is 11 TPM and MDM4 is 24 TPM and CDKN1A (p21) is 225 TPM. According to the single cell RNA-seq of the human skin, the p73 RNA is under detection threshold . However, to release single cells, the skin is digested at 37 °C for 2 h, which can lead to RNA degradation, especially for the unstable transcripts .
Given the low p73 RNA levels in skin and other epithelial organs, the question arises: how reliable is p73 antibody staining? We have thoroughly analyzed the available data using different antibodies and conclude that p73 expression has the same characteristic patterns in skin and other epithelial organs, as defined by experiments using antibodies raised against distinct p73 epitopes in independent studies.
Analysis of the literature and web resources revealed that there are 4 different antibodies against different epitopes that generate exactly the same p73 staining pattern in skin. Chemicon rabbit polyclonal AB7824 , Leica Biosystems mouse monoclonal (Leica Biosystems Cat# NCL-p73, RRID:AB_563940), and two Abcam’s: EP436Y (aa 50–150) , EPR19884 (C-terminus, skin IHC image is at the Abcam website). Reactivity of the Abcam’s antibodies to p73 is validated by the knockout/knockdown and overexpression experiments either by manufacturer or by independent investigators using immunohistochemistry and immunoblotting . The main concern is that p73 antibodies can cross react with p63. There are several lines of evidence against it. Firstly, p63 immunostaining appeared earlier during skin development than EP436Y p73 staining . Secondly, during oocyte development, EP436Y p73 antibody stains granulosa cells and these staining disappear in the p73 −/− mice. In contrast, p63 stains primordial follicles in both p73+/+ and p73 −/− mice . These data were supported by immunoblotting of the wild type and p73 −/− ovaries. Third, the p73 EP436Y antibody does not react with p63 by immunohistochemistry  as it stains a subset of the p63 positive cells in human and mouse epidermis . Moreover, closer examination revealed that a few cells have brighter p73 staining than p63 . The absence of reactivity of the EP436Y p73 antibody to p63 by immunoblotting was also demonstrated . Therefore, it is possible to conclude that at least EP436Y p73 antibody is reliable.
In addition, several tissues display strong p63 staining and no AB7824 or RRID:AB_563940 p73 antibody staining (breast, prostate, placenta, seminal vesicles) was observed, whereas other epithelial organs display the same staining of the basal cells as in skin (Figure 2) http://www.proteinatlas.org) .
Thus, p73 protein is expressed in the basal layer of epithelial organs, whereas p63 appears in the basal as well as in the differentiated cells.
Differential and common functions of p63 and p73 are reflected by the proteins they interact with. In order to have the available information of their binding partners, we downloaded protein-protein interactions for the p53, p63 and p73 from the BioGrid database. We excluded interactions from the genetic screens, keeping only physical associations. Shared and unique interactions are presented as an Euler diagram (Figure 3A).
Furthermore, we annotated shared and unique p53, p63 and p73 interactions using Metascape software . The MCODE algorithm within Metascape detects densely interconnected regions in protein-protein interaction networks that may represent distinct molecular complexes (Figure 3B–D). In addition, Metascape applies enrichment analysis of pathways and processes to identified networks of protein interactions and uses the top three enriched terms for the annotation of biological roles of the putative molecular complexes. In addition, we compared protein interaction data with mRNA expression using publicly available TCGA RNA-seq datasets and single cell sequencing of human skin (Supplementary File) . Below we review what is known about the functional impact of the p53 family protein interactions focusing on the biology of the epithelial tissue.
Figure 3. Annotation of the unique p73 and p63 protein interactions reveal distinct physiological functions. (A). Euler diagram of the p53, p63, and p73 protein interactions from the BioGrid database. Coloring of the Euler diagram corresponds to the coloring of dots representing proteins in the subsequent sections of the figure (B–D) describing protein interactions. (B). Annotation of proteins uniquely interacting with p63 revealed putative complexes involved in the protein catabolism and synthesis including protein folding and transport (CCT8, HSPA5, HSP90AB, VCP, CCT8, ATP5F1A), translation—EEF1A1 and degradation (PSMD2, PSMD14). The second complex contains ribonucleoproteins involved in the mRNA splicing and Nf-kb transcription factors and the third represents ribosomal proteins involved in translation and nonsense mRNA decay (UPF2). (C). Annotation of proteins interacting with p73 revealed putative complexes involved in the regulation of the cell cycle (Rb, CDKN1A (p21), Cyclin D1, Cyclin A2, SPK2) as well as mitosis (BUB1B, MAD2L1, DSN1, STAG1). The next complex is FOXM1 mediated cell-cycle regulation. Interestingly, the adaptive immunity cluster is composed exclusively of the p73 interacting proteins and three chromatin regulators can bind both p53 and p73. (D). Annotation of proteins interacting with both p63 and p73 revealed putative complexes involved in the regulation (repression) of the p53 stability (MDM2, MDM4,RCHY1, ITCH, SUMO1, UBE2I) and transcriptional activity (HDAC2) transcriptional coactivators (EP300, KAT2B, HIPK2), DNA repair (BRCA1), DNA binding transcription factors (HIF1a, YAP, ATF3, p63, p73). Two proteins that are bound by p63 and p73 and not p53 are CDC20—a pivotal protein in the mitosis and FUS—the splicing regulating protein.
It is widely accepted that transcription is regulated by the interplay of the multiprotein complexes . Moreover, due to the presence of dominant activators or repressors co-occupying DNA regulatory sites, functions of a particular transcription factor might be limited to a subset of its binding sites . For the p53 transcription factors family, interactions with distinct proteins may skew the binding of p53, p63 and p73 isoforms towards different sites in the genome and determine their specific activities .
As discussed above, p63 and p73 colocalize in the basal layers of the epithelial organs, while p73 is not present in the differentiating cells. Indeed, the p63 and p73 can hetero-oligomerize through their tetramerization domains, and the resulting hetero-tetramer consisting of two p63 and two p73 molecules is a thermodynamically preferred structure .
Calcium differentiation of the keratinocytes induces p73 protein and promotes its interaction with p63 at days 2–4 of calcium-induced differentiation . The interaction between p63 and p73 is initiated just before cell culture stratification and might play a role in the induction of differentiation markers. However, it is important to mention that at 48 h in calcium-enriched media, cultures of human keratinocytes typically contain a mixture of proliferating and differentiated cells, thereby, it is possible that interaction between p63 and p73 occurs in the proliferating cells, consistent with in vivo p73/p63 colocalization in the basal cells.
In support of this data, it was shown that in HaCat and NHEK cells p73 delta and gamma isoforms as well as p63 were able to activate involucrin and loricrin promoters suggesting that both of these genes may play a role in the induction of differentiation .
Interactions of DNp63a and TAp73β led to reduced TAp73β transactivation of the luciferase reporter in neuroblastoma cells . However, in skin and other epithelia p73 and p63 are mostly represented by the DN- isoforms, whose transcriptional activation is determined by interaction with other transcription factors . The distinction between DN- isoforms of p63 and p73 is demonstrated by co-expression of the KLF4 with DNp73β (but not DNp63α or either KLF4 or DNp73β alone) that dramatically induced fibroblast–keratinocytes reprogramming and expression of keratinocyte specific marks (KRT5, KRT14, FLG, SPN) in both fibroblasts and mesenchymal triple-negative breast cancer cell line MDA-MB-231 .
Several signaling transcription factors/pathways cooperate with the p63/p73 to regulate gene expression. Of these, the most studied are CEBP/AP1/NfKb transcription factors that play a pivotal role in the epithelial differentiation .
It is likely that in skin cells DNp63a and DNp73α interact with either c-Jun or c-Rel and localise to the AP-1 binding sites similarly to what was observed in neck squamous cell carcinoma and osteosarcoma Saos cells . In the squamous cancer cells, proinflammatory cytokine TNF-α modulates interaction of DNp63α/TAp73 with c-REL and induces redistribution of TAp73 from the p53 to AP-1 DNA binding sites repressing TAp73 pro-apoptotic functions and inducing an oncogenic gene expression in the absence of functional p53 . Thus, TAp73 in complex with AP-1 versus DNp63a switch from anti- to pro-oncogenic activities .
Accordingly, a consensus p53 binding motif is found next to the composite CEBP/c-Jun binding site enriched in the promoters bound by CEBPb and c-Jun in mice basal keratinocytes . Notably, promoters with differentiation induced RNA polymerase 2 binding are preferentially occupied by CEBPb and c-Jun. However, the combinatorial p53/CEBP/AP1 motif was overrepresented in the CEBPb bound promoters also active in the basal keratinocytes.
In addition, it was shown that p63 interacts with CEBPb and Nf-kb in lung cancer cells upon cigarette smoke extract exposure and regulates COX-2 expression, supporting possible functional consequence of the p53/CEBP/AP1 binding sites proximity .
The localization of the p63 binding sites as suggested by the ChIP-seq data does not change upon the keratinocytes differentiation and modulation of the activity of the p63-bound enhancers is likely determined by other protein complexes that cooperate with p63 . Upon progression through the cell cycle during epidermal differentiation, high expression of the chromatin remodelers occurs in the mitotic phase of the cell cycle, concomitant with induction of the early differentiation markers KRT1 and KRT10 , possibly marking a subset of epidermal asymmetrical division, in which p63 plays a pivotal role . Accordingly, p73 interacts with chromatin remodelers and transcription factors (Jun, RB1, SP1, HDAC1) involved in the cell cycle regulation, FOXM1 pathway and mitosis (Figure 3), suggesting that the p73 specific interactions might modulate p63 activity either by promoting or repressing transcription .
Indeed, additional evidence from the chromatin immunoprecipitation sequencing (ChIP-seq) data suggests that p63 and p73 occupy a similar set of promoters in cancer cell lines . Even though genome-wide colocalization of the p63 and p73 in the primary keratinocytes or in skin has not been yet analyzed, it is assumed that p63 and p73 will be localized similarly . Thus, both hetero-tetramerization mediated and distinct DNp63 and DNp73 DNA binding are likely to occur in the epithelial organs. While a plausible hypothesis would be that these complexes differentially regulate transcription in the basal epithelial and in the differentiating cells, details of the regulation are just beginning to emerge .
Complexes of 52 proteins that interact with both p63 and p73 correspond to the shared pathways regulating the p53 family. These include regulators of the p53 stability and transcriptional activation by ubiquitination and sumoylation: MDM2, MDM4, RCHY1, ITCH, SUMO1, UBE2I , transcriptional co-regulators: HDAC2, EP300, KAT2B, HIPK2 , DNA repair gene BRCA1 that interacts with p63  while repressing TAp73 expression and chemoresistance , DNA binding transcription factors: HIF1a, YAP, ATF3, p63, p73 . Two proteins that are bound by p63 and p73 and not p53 are CDC20—a pivotal protein in the mitosis  and splicing regulating protein FUS .
It was shown that HDAC1 and HDAC2 depletion phenocopy p63 knockout in mice derepressing negatively regulated DNp63 targets CDKN1A(p21), SFN(14-3-3σ), and CDKN2A . Interestingly, classical p53 target and cell cycle inhibitor CDKN1A  is expressed strictly in the suprabasal cells in human epithelial organs, especially in the esophagus and vagina. A p53/p73/CDKN1A inhibitor MDM2 is also in this group with unrestricted expression in all skin cells .
Analysis of 89 p73 or p53 interacting proteins revealed two putative protein complexes involved in the cell-cycle progression, chromatin modification complex and adaptive immunity complex (Figure 3C).
Analysis of 36 p73-only interacting proteins revealed clusters of interactors pointing to their potential role in the mitotic prometaphase, and resolution of the sister chromatin cohesion and segregation. These include protein of the spindle assembly checkpoint complex (Figure 3C, complex 1) (CCNB1, MAD2L1, Bub1B, STAG1, DSN1) .
In turn, 53 p73 interacting proteins that are also known to bind p53 but not p63, represent histone modifiers including: histone deacetylases: HDAC1, CHD3 , histone acetylases: KAT5, CREBBP , H3K4me2 or H3K9me2 lysine demethylase KDM1A . The p73 and p53 interact with proteins involved in G1/S transition of the cell cycle (SP1, MYC, CCND1, RB1)  as well as proteins involved in the G2 to M transition Cyclin A and CDK1 . Also, p73 and p53 interact with Aurora A kinase that represses p73 mediated transcriptional activation and SAC functions .
This group includes a classical p53/p73 target p21 (CDKN1A)—a major regulator of the G1/S cell cycle progression that acts by dephosphorylation and repression of multiple proteins involved in the cell cycle progression . Accordingly its expression is restricted to the suprabasal layer of skin, vagina and other epithelial organs with no proliferation . However, p21 protein does not interact with p73 directly, only p73 binding to the p21 promoter has been described, although p21 interacts with many proteins that do interact with p73, including SP1, RB1, CyclinD1 and CyclinA2.
Analysis of 58 p53/p73 interacting proteins that are highly expressed in the skin  revealed putative protein complexes involved in the cell cycle regulation and FOXM1 pathway: CCND1, CCNA2, CCNB1, CDK1, CDK2, Rb, transcriptional regulation: MYC, HDAC1, KAT5, CREBBP, cellular response to stress: HSPA4, HSPA9, NFYB, microtubule organizing center and G2/M transition:PLK1, AURKA.
Gene annotation of 187 p63-only interacting proteins reveals significant fraction of them being involved in RNA splicing and TGFb response. Chrondate and nervous system developmental pathways are highly represented in this group as well.
The TAp73 interacting with APC/C(CDC20) repress progression to anaphase and aberrant chromosome segregation, while APC/C(CDC20) interaction with DNp63a promotes DNp63a degradation and is required for keratinocytes differentiation . Notably, p63 interacting APC/C components CDC20 and ANAPC2 are located strictly in a few cells of the basal and suprabasal layer, likely marking the last division of the cells committed to differentiation.
However, in the suprabasal layer, after keratinocytes differentiate, p63 protein is abundant suggesting that DNp63a isoform that is sensitive to the APC/C (CDC20, CDH1) mediated degradation is more restricted to the basal keratinocytes, whereas DNp63 present in the differentiated layers is likely other isoforms such as DNp63b and DNp63g .
Interestingly, one of the p63 interacting proteins is chaperonin containing TCP1, subunit 8 (theta) (CCT8) , that according to high throughput screening in turn interacts with gamma tubulin1  and component of SAC complex BUB1 , suggesting yet unexplored implication of the p63 in the spindle checkpoint that might interact with or compete with p73.
Among proteins interacting with p63 are WWP1, C/EBPb, SMAD4 and KLF5. The WWP1 is an E3 ubiquitin ligase that is required for p63 transcriptional activity and WWP1 depletion causes cell cycle arrest . Consistently, WWP1 displays an interesting staining pattern in skin, being nuclear in the basal cells and strictly peri-nuclear in the differentiated keratinocytes. Moreover, WWP1 ubiquitinates and leads to degradation and repression of KLF5  that is highly expressed in the nuclei and cytoplasm of the basal and differentiated keratinocytes and of SMAD4  that is localised to the cytoplasm in the epithelial tissue.
C/EBPb and C/EBPa are pivotal factors in the induction of genes during keratinocytes differentiation . While, interaction between C/EBP and DNp63 leads to the induction of COX-2 expression in the lung epithelial cells upon smoke extract exposure , COX-2 induction in the skin keratinocytes leads to PGE2 induction and inflammation in response to ultraviolet irradiation . Therefore, the functions of WWP1 interacting proteins is this subgroup are likely to be a switch of transcriptional regulation from basal to differentiated program.
Another protein interacting with p63 and WWP1 is KLF5. Both KLF5 and KlF4—a well-known pluripotency Yamanaka factor —are expressed in skin with distinct staining patterns. The KLF5 decorates nuclear borders just like WWP1, whereas KLF4 is highly expressed throughout the differentiated cells. KLF4 is not known to interact directly with p63, however, it was shown that p63 regulates KLF4 expression and combined p63/KLF4 expression differentiate fibroblasts into keratinocytes . In contrast, the role of p73 in reprogramming is controversial, because while one study used p73 depletion technique to demonstrate that p73 has no role in reprogramming , the other concluded that overexpression of DNp73β and KLF4 leads to reprogramming .
Common p63 and p53 interacting proteins relate to autophagy and protein folding (CHEK2, RELA, HSP90AB1) Figure 3D. Among p63 interacting proteins are several 60S ribosomal subunits including RPLP0, RPL3, RPL12, RPL13A, RPL24 as well as 40S ribosomal protein SA and UPF2 that controls nonsense-mediated mRNA decay . In addition, p63 binds several human ribonucleoproteins HNRNPK, HNRNPL HNRNPA3, HNRNPR , HNRNPD  and HNRNPAB  suggesting that p63 isoforms function in the splicing regulation (Figure 3D).
It was not possible to annotate 971 unique p53 interactions by the Metascape. When we limited our analysis to 604 proteins that are detected in skin by the single cell RNA-seq , the Metascape recognized 11 putative protein complexes with functions ranging from the DNA repair to mRNA splicing and metabolism (Supplementary Excel File) consistently with known p53 functions in the genotoxic stress responses, transcription and RNA regulation pathways .