Functional genomics of AP-2α and AP-2γ in cancers: in silico study
BMC Medical Genomics volume 13, Article number: 174 (2020)
Among all causes of death, cancer is the most prevalent and is only outpaced by cardiovascular diseases. Molecular theory of carcinogenesis states that apoptosis and proliferation are regulated by groups of tumor suppressors or oncogenes. Transcription factors are example of proteins comprising representatives of both cancer-related groups. Exemplary family of transcription factors which exhibits dualism of function is Activating enhancer-binding Protein 2 (AP-2). Scientific reports concerning their function in carcinogenesis depend on particular family member and/or tumor type which proves the issue to be unsolved. Therefore, the present study examines role of the best-described AP-2 representatives, AP-2α and AP-2γ, through ontological analysis of their target genes and investigation what processes are differentially regulated in 21 cancers using samples deposited in Genomic Data Analysis Center (GDAC) Firehose.
Expression data with clinical annotation was collected from TCGA-dedicated repository GDAC Firehose. Transcription factor targets were obtained from Gene Transcription Regulation Database (GTRD), TRANScription FACtor database (TRANSFAC) and Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining (TRRUST). Monocle3 R package was used for global samples profiling while Protein ANalysis THrough Evolutionary Relationships (PANTHER) tool was used to perform gene ontology analysis.
With RNA-seq data and Monocle3 or PANTHER tools we outlined differences in many processes and signaling pathways, separating tumor from normal tissues or tumors from each other. Unexpectedly, a number of alterations in basal-like breast cancer were identified that distinguished it from other subtypes, which could bring future clinical benefits.
Our findings indicate that while the AP-2α/γ role remains ambiguous, their activity is based on processes that underlie the cancer hallmarks and their expression could have potential in diagnosis of selected tumors.
Following cardiovascular diseases, cancer is characterized by the second highest global mortality rate of any disease . Many paradigms have been developed to explain the process of cancer development [2,3,4]—among them is the molecular theory of carcinogenesis associated with tumor suppressors and oncogenes , having its origin in the twentieth century . It postulates that carcinogenesis is determined by alterations in cancer regulatory genes, of which two crucial groups are tumor suppressors and oncogenes, both responsible for apoptosis and proliferation regulation  being utmost importance in the model of cancer platform . Nonetheless, occurrence of cancer-related genes demonstrating concurrently suppressive and oncogenic characteristics is definite for instance in transcription factors (TFs) . Of these, probably the best known is Guardian of the Genome (p53 protein), whose correct activity requires that all of its subunits are not mutated and that no dominant negative effect occur . In addition, the members of the Activating enhancer-binding Protein 2 (AP-2) family also demonstrate ambiguity; of these, AP-2α exhibits dualism of function depending on its related signaling pathway [11, 12]. Another AP-2 protein, AP-2γ, was thought to have oncogenic activity since its transcriptional activity is inhibited by WW Domain Containing Oxidoreductase (WWOX) tumor suppressor after its sequestration outside the nucleus . However, the other reports indicating that AP-2γ induces expression of p21 protein suggest that it may also demonstrate anti-tumoral traits , presenting contrasting function which prompts consideration.
The entire AP-2 family fits within the basic Helix–Span–Helix (bHSH) superclass and comprises the five members (AP-2α, AP-2β, AP-2γ, AP-2δ, AP-2ε) encoded by TFAP2A-E genes . The activation domain is located in the amino terminus, whereas the sequences responsible for dimerization and binding of DNA are on the side of carboxyl terminus . All members recognize specific G/C-rich evolutionarily conserved motifs i.e. GCCN3/4GGC, GCCN3/4GGG or CCCCAGGC [17, 18] while the binding of transcription factors’ themselves is dictated by a proline-rich motif located in the activation domain (excluding AP-2δ, which lacks these critical residues) . Cytogenetically, AP-2 family members are encoded on the sixth chromosome, except AP-2γ and AP-2ε which are on the twentieth and first chromosomes, respectively . When functioning correctly, the AP-2 transcription factors regulate appropriate gene expression during early developmental processes such as face, eye or limb development . It was suggested that in the case of mutation, the loss of TF activity of AP-2 members can lead to the impairment of proliferation, differentiation and apoptosis processes , suggesting AP-2 activity may play role in development. Indeed, both AP-2α and AP-2γ have prognostic value for some tumor types. AP-2α overexpression is recognized as a prognostic indicator of shorter patient survival in papillary thyroid carcinoma and epithelial ovarian cancer but longer survival in gastric adenocarcinoma [20,21,22]. In terms of AP-2γ, its overexpression was also associated with poorer overall survival in breast cancer patients  with the other literature data to support this statement and enrich it with correlation presenting worse anti-hormone therapy response . The other case of chemoresistance (to 5-fluorouracil) was shown in colorectal cancer upregulating AP-2γ  while endometrial cancer example demonstrated that knockdown of this AP-2 member sensitizes cells to megestrol acetate via Estrogen receptor alpha (ERα) expression upregulation .
Literature data indicates that the activity of AP-2 in carcinogenesis is very much dependent on specific family members. The most recent scientific reports concern the first three transcription factors, while the least is known about AP-2δ or AP-2ε. The best described factors, AP-2α and AP-2γ, play distinct roles depending on tumor type, as noted in a recent review ; however; this summary has since been enriched with additional tumor models not available at that time. These findings confirm that AP-2α activity is influenced by tumor tissue type: while it was found to demonstrate oncogenic activity in cervical, gallbladder and ovarian cancer [27,28,29], it has been also seen to act as a suppressor in many other tumors [11, 30,31,32,33,34,35]. AP-2γ has been conclusively demonstrated to play an oncogenic role in breast cancer  (which coincides with previous evaluation) and furthermore expanded with its unfavorable characteristics in colorectal cancer . The last representative of the AP-2 family abundantly described in the literature is beta member. Latest data regarding contribution in cancer indicates that AP-2β overexpression has been found to promote tumor growth in both breast and thyroid cancer and predicted poor prognosis or tumor progression, respectively [37, 38]. Finally, although little is known of AP-2δ and AP-2ε, the former is thought to be associated with progression and genomic instability of prostate cancer while the latter acts as a tumor suppressor in neuroblastoma [39, 40].
The precise functionality of AP-2 factors in cancer clearly remains unknown. Therefore, the aim of the present study was to clarify the role of transcription factors AP-2α and AP-2γ in carcinogenesis using samples (RNA-seq) acquired from 21 cancers of The Cancer Genome Atlas (TCGA) and Monocle3 or Protein Analysis Through Evolutionary Relationships (PANTHER) bioinformatics tools. The study identifies the processes which are differentially regulated between the studied cancer types, compare the findings with those of normal tissue and identifies whether any differences exist in the expression of target genes for these factors in specific types of tumors.
Acquisition of tumor patients’ data
RNA-seq expression data with corresponding clinical annotation was collected from 21 tumors (Table 1) of TCGA-dedicated GDAC Firehose Repository (level 3 RNA-seqV2, RSEM normalized, data status of 28th Jan 2016 available at https://gdac.broadinstitute.org/). Patients that lacked expression or clinical data were discarded from the study. Available normal, paired solid tissues were additionally retrieved through R-dedicated package TCGA-Assembler .
AP-2α and AP-2γ target genes identification
Targets for both transcription factors were identified through several databases: GTRD (version 19.10 [42, 43]), TRANSFAC (version 2019.2) and TRRUST (version v2). Excluding duplicates, there were 4810 and 5175 targets for AP-2α and AP-2γ, respectively.
Global profiling of the tumors and ontological annotation
Phenotype heterogeneity between selected tumors accompanied by normal tissues was studied and visualized by applying the UMAP method preceded by principle component analysis (PCA) regarding the expression of the targets of both AP-2α and AP-2γ using the Monocle3 R package (https://cole-trapnell-lab.github.io/monocle3/) . The analyses of AP-2α and AP-2γ targets were performed as two separate entities accordingly. The PCA pre-processing step (preprocess_cds()) was done with the dimensionality of the reduced space of 100 (num_dim). Reduction of dimensions (reduce_dimension()), and clustering of the individuals (cluster_cells()) within spaces were applied with the UMAP algorithm for dimensionality reduction method upon which to base clustering (reduction_method). The clusters of individuals were compared with graph_test() function based on Moran's I spatial autocorrelation analysis with knn neighbor graph and q-value threshold of 0.05. Furthermore, the genes varying across the clusters selected at the previous step were grouped into modules through Louvain community analysis (find_gene_modules()) with parameters set to default. The modules were clustered in two ways: 1) clustering involving all of the individuals enabling to compare between tumors and tumor vs non-tumor samples of specific cancer type and 2) clustering of the modules restricted to the BRCA foci enabling to differentiate between PAM50 subtypes, separately for AP-2α and AP-2γ targets. Finally, the results were visualized with pheatmap() and clustered with the Ward D2 method. The whole pipeline was performed according to the Monocle3 tutorial (https://cole-trapnell-lab.github.io/monocle3/). To annotate the findings in the context of the biological processes in which they are involved, overrepresentation test available at PANTHER Classification System (http://www.pantherdb.org/) with Fischer's Exact Test and Bonferroni correction for multiple testing were performed. The methodology is graphically summarized in Fig. 1.
Identification of global differences using toolkit for spatial analysis
Considered cohorts were clustered according to designed variable discriminating tumors and their corresponding normal tissues (marked with an additional prefix "n_"). Their distribution across Uniform Manifold Approximation and Projection (UMAP) dimensions was presented according to the genetic targets for a given transcription factor (AP-2α or AP-2γ). The analysis based on AP-2α identified a group of cancers distinct from other separated clusters but to a lesser extent from the corresponding normal tissue: LIHC, KIRP, KIRC, PRAD, SKCM, OV, THCA, BRCA, STAD, LUAD and PAAD. In some cases, one cancer type coincided in the plot with another (forming a “cluster”) e.g. GBM and LGG, COAD and READ. Surprisingly, the remaining tumor reservoir (BLCA, CESC, UCEC, ESCA, HNSC, LUSC) formed a heterogeneous mixed cluster together with their normal tissues, albeit with deviations for ESCA and LUSC. The same tendencies as of AP-2α were observed in terms of AP-2γ targets, although with minor differences—these are collectively presented in Fig. 2.
In the next stage, distribution in the UMAP1 and UMAP2 dimensions was further analyzed with regard to differences in expression of the transcription factor itself. Global similarities were observed in the case of both AP-2α or AP-2γ expression; however, some exceptions were observed e.g. AP-2α expression was higher in SKCM while AP-2γ was higher in THCA (Fig. 3a, b, respectively). Evident differences in expression were observed between kidney tumors (KIRC, KIRP) and normal kidney tissue; this is not surprising in the case of AP-2α which has been suggested as a biomarker in renal carcinoma [45, 46] however our findings suggest that AP-2γ may have the same properties.
Furthermore, due to noticeable separation of breast cancer samples visible for both target gene lists (Fig. 2), possession of extraordinary tumor profiling signature (Prosigna Breast Cancer Prognostic Gene Signature Assay—PAM50) prompted us to perform additional analyses to assess the subtypes’ separation by TFs target genes. For both AP-2α or AP-2γ, differences were observed in the PAM50 classifiers of breast cancer indicating that the basal-like subtype was different from the others (Fig. 4a, b, respectively).
Lastly, heatmaps were generated to arrange the transcription factor target genes into modules with common expression profile, indicating how these modules differ between clusters; this approach allowed more specific differences between individual cancers, the tumor and the corresponding normal tissue or intrinsic subtypes to be identified (Figs. 5, 6). The content of modules from particular analysis type has been summarized in the Additional file 1. At first glance, signaling by the AP-2 targets in tumor tissue is congruent to signaling in the corresponding normal tissue for PAAD, LIHC, THCA and SKCM (in the case of AP-2α) or for PAAD, THCA and SKCM (in the case of AP-2γ). However, the analysis of signaling through TF targets in the remaining clusters indicated a large number of modules that could be used to distinguish the tumor from normal tissue (the latter taken as a reference), or the basal-like BRCA subtype from the luminal A/B and HER2-enriched subtypes. The proposed modules taken for further comparisons between tumor and normal tissue or between different tumor tissues are summarized in Table 2 (AP-2α) and Table 3 (AP-2γ). In the case of differences between BRCA subtypes, the modules 6, 12 and 13 were taken for subsequent consideration for the AP-2α target gene list while 5, 13, 19 modules for AP-2γ.
Another approach was to distinguish different cancer types with similar sites of origin. The following comparisons were used for gene modules generated from AP-2α and AP-2γ targets: COAD versus READ, LGG versus GBM, LUSC versus LUAD and KIRC versus KIRP. For AP-2α, COAD and READ gave similar results; however, LGG versus GBM were differentiated by modules 10 and 13, LUSC versus LUAD for modules 6, 9, 11, 15 and KIRC versus KIRP for modules 1, 3, 4, 6, 13, 16, 19. Regarding AP-2γ, COAD versus READ were differentiated by module 9, LGG versus GBM with module 5, LUSC versus LUAD with modules 7, 10, 12, 15 and KIRC versus KIRP with modules 3, 4, 12, 16.
The differences between tumors were then identified from a heterogeneous mixed cluster. Using modules 2, 11, 19 with target genes for AP-2α, the entire selection was divided into two different “communities”—LUSC, ESCA, HNSC and UCEC, BLCA, CESC. The same tendency could be explained using module 17 of AP-2γ targets with additional dissimilarities in UCEC versus LUSC or ESCA versus BLCA which can be explained by the inverse of module 15 or 6, respectively.
Gene ontology of modules distinguishing specific comparisons
Analyzing a multitude of comparisons in the form of modules diversifying groups (from the previous section) led to the extraction of genes for their subsequent ontological analysis. The greatest attention was paid to modules showing the opposite tendency in the context of expression of genes contained therein. As previously, the following distinctions were made: between the cancer and the corresponding normal tissue, between similar tumors or between PAM50 classifiers for breast cancer.
Comparisons by means of target genes for AP-2α transcription factor
Differences between tumors and normal tissues indicated changes in processes (e.g. apoptosis, proliferation, cell migration, cell cycle, cytoskeleton organization, cell adhesion or angiogenesis) and signal transduction pathways (e.g. ERBB, Wnt, MAPK and Notch) which are visualized in Fig. 7. In regards to alterations between different types of tumors, they were approximate to above mentioned however additional processes such as regulation of cell shape were also noticeable (Fig. 8). Further data such as specific genes that described a particular module during ontological analysis, as well as statistical significance, are reviewed in Additional file 2.
Comparisons by means of target genes for AP-2γ transcription factor
In the case of second AP-2 family representative, tumor versus normal tissue changes also affected processes (e.g. autophagy, adhesion, vesicle budding, extracellular matrix organization, apoptosis, cell growth and spreading) and signaling pathways (e.g. TNF, EGFR, Wnt, JAK-STAT, mTOR, NFkB and TGFβ) which is summarized in Fig. 9. Regarding variation between tumors, no additional processes other than those mentioned above were identified, even though different modules were used (i.e. module 6) (Fig. 10). The specific genes that characterizing individual modules through ontological analysis, together with their statistical significance, are recapitulated in Additional file 3.
Comparisons of breast cancer intrinsic PAM50 subtypes
Pursuing the scheme of previous subsections, gene ontology of modules distinguishing basal-like BRCA signature from other subtypes indicated differences e.g. in apoptosis, cytoskeleton organization, autophagy, but also in EGFR, TNF, Wnt or cadherin signaling pathway. The results are summarized in Fig. 11 and Additional file 4.
Validation of the Results
Validation of the findings based on independent adequate cohorts for all considered tumors was not performed due to lack of relevant expression data.
Carcinogenesis is shaped by both tumor suppressors and oncogenes where one group of crucial proteins that demonstrate such duality of function are transcription factors. Available literature indicates that the character of AP-2 factors depends on the particular family representative or considered tumor type. Assurance definitely applies to fact that AP-2α and AP-2γ have been described to the greatest extent compared to other family members.
Unlike previous papers, the present study examines the activities of the target genes of AP-2α and AP-2γ in more than 20 solid tumors from the TCGA database. The observed alterations in the expression of transcription factors and their effect on target genes identified processes that can potentially distinguish the tumor from its corresponding normal tissue or from other tumor types.
Regarding the first type of comparisons, a multitude of differences were observed; these allowed the identification of more than one module that could distinguish tumor tissue from corresponding normal tissue, with ontological descriptors concerning both biological processes and molecular pathways. A few modules from AP-2α target genes list explained the differences for several tissues at once: modules 2, 11 and 12 respectively differentiated CESC and UCEC, HNSC and KIRC, as well as COAD, ESCA and STAD from their corresponding normal tissues. Module 2 has been described as implicated in cell cycle regulation as well as the apoptotic pathway, with the cell cycle being differentially regulated between cervical cancer versus solid normal tissue samples  or uterine corpus endometrial carcinoma versus normal controls . Following, in module 11 the processes involving e.g. MAPK activity was noticed which could support previous observations that MAPK may play a role in HNSC and KIRC carcinogenesis [49, 50]. Lastly, module 12 was annotated with processes involved in cell cycle or cell adhesion. Differences in cell cycle regulation were previously found between tumor tissue and paired normal samples of colon adenocarcinoma , while cell adhesion processes were found to distinguish tumor samples from normal tissue in esophageal carcinoma  or gastric cancer .
The remaining modules selected for gene ontology analysis indicated differences in PRAD, LUSC, BLCA, GBM, READ, LUAD, BRCA and KIRP between tumors and their matched solid normal tissues. Module 18 shown alterations in apoptosis, the ERBB pathway and regulation of cytoskeleton organization between PRAD and its paired control which is consistent with previous findings for programmed cell death  and ERBB signaling . The other pathway shown was Wnt, describing lung squamous cell carcinoma versus normal tissue differences with module 15 which corresponded to literature data on cellular levels of this protein . Furthermore, “regulation of cell cycle” or “response to drug” descriptors of module 13 explained bladder cancer versus normal bladder tissue; this is consistent with current scientific reports where the drug resistance is dependent on the JUND representative of the Jun family [57, 58], which was a part of AP-2α target genes list forming module 13 along with two other Jun members. Module 10 concerns the misfolding of proteins which appeared to be crucial mechanism of lipid storage or increased cellular proliferation and varied between glioblastoma multiforme and corresponding normal tissue . Module 8 revealed differences in intrinsic apoptosis, among others, between rectum adenocarcinoma and controls as noted previously ; the same signaling pathway was annotated in module 6 explaining lung adenocarcinoma versus normal lung, although with an additional cell adhesion descriptor which might be related to changes in adhesive properties during epithelial-to-mesenchymal transition, as confirmed via AP-2α–promoted tumorigenicity of LUAD . Finally, kidney renal papillary and breast invasive carcinomas were distinguished from their paired normal tissues using module 1 and 3, respectively. Among the others, signaling pathways such as Wnt, Notch (for KIRP) or ERBB (for BRCA) were implicated in explanation which is coherent with other research [62, 63].
Since the magnitude of the results for AP-2γ is greater than in the case of AP-2α, obtained differences were confronted with the literature only in terms of signaling pathways. Differences in certain pathways, such as TNF or ERBB (for liver hepatocellular carcinoma and prostate adenocarcinoma), EGFR or Wnt (for lung squamous cell carcinoma), JAK-STAT (for kidney renal clear cell carcinoma), NFkB (for breast invasive carcinoma) and TGFβ, mTOR or NFkB (for stomach adenocarcinoma) were noted between tumor and normal tissue. Of these, TNF for LIHC, NFkB for BRCA and TGFβ for STAD have been identified previously [64,65,66].
Both AP-2α and AP-2γ demonstrated differences associated with signal transduction pathways or biological processes in tumor versus tumor comparisons. Regarding first transcription factor, many of the modules generated through Monocle3 analysis can distinguish KIRC from KIRP; in this, ontological analysis identified changes in processes e.g. cell adhesion, proliferation, angiogenesis, cell migration, cell cycle or pathways i.e. Notch, Wnt, ERBB signaling (Fig. 8). Differences in adhesion or cell cycle have been noted between two renal cell carcinomas, showing specific gene expression to be lower in KIRP . Increased expression of AP-2α targets genes from module 13 (ontologically indicating regulation of cell cycle, proliferation or development) in KIRC compared to KIRP was also noticeable in LGG versus GBM comparison although with inverse tendency (Fig. 8). Considering the cell cycle differences, the alteration frequencies indicate on average higher percentage of pathway disruption in GBM (86%) compared to LGG (46%) . The cell adhesion and intrinsic apoptotic signaling pathway (module 6) remained on trend in both KIRC versus KIRP and LUSC versus LUAD, implying that lung carcinomas differ not only in terms of adhesion or apoptosis but also regarding cell migration or Wnt signaling pathways, as demonstrated by module 15 (Fig. 8). This is consistent with a previous comparison of genes linked to adhesion and migration that were found to be differentially expressed between LUSC and LUAD . Of the two types of lung cancer, LUSC was also found within a mixed cluster (Fig. 2a) created by more similar tumors (LUSC, ESCA, HNSC, UCEC, BLCA, CESC). Nevertheless, the tumors could be distinguished using modules 2, 11 and 19 which included genes related to cell adhesion, regulation of cell cycle arrest and inactivation of MAPK activity.
Concerning differences outlined through AP-2γ targets, identical subgroups could be distinguished within the mixed tumor cluster (Fig. 2b) formed by the same six cancer types described above—distinction is explained by module 17 that post-ontologically shown differences in cell adhesion (Fig. 10), which is concordant with the analysis of the same tumors through AP-2α targets (Fig. 2a). Additionally, module 6 or 15 distinguished two subdivisions inside the cluster i.e. ESCA versus BLCA or UCEC versus LUSC, respectively. Using ontology descriptors, changes concerned regulation of mTOR signaling and extrinsic apoptosis for module 6 or EGFR and Wnt pathways for module 15 (Fig. 10). In regards to the latter, the alteration frequencies indicate on average lower percentage of Wnt pathway disruption in LUSC (18%) compared to UCEC (47%) . Other comparisons applied to COAD versus READ explained by vesicle budding (module 9) or LGG versus GBM elucidated through mTOR signaling, apoptosis and cellular adhesion (module 5) with the latter confirmed by another bioinformatics analysis . Lastly, numerous modules distinguished KIRC from KIRP or LUSC from LUAD with module 12 (response to hypoxia, regulation of autophagy, adhesion-dependent cell spreading) differentiating both. Referring the literature, renal cell carcinomas have already been differentiated in terms of hypoxia indicating specific gene expression to be higher in KIRC . Other modules for KIRC versus KIRP comparison have been associated with processes such as autophagy (module 3), hormone-mediated signaling pathway (module 4) and cell growth regulation or NFkB signaling (module 16). For LUSC versus LUAD, they differ according to genes associated with Wnt pathway (module 10), regulation of cytoskeleton organization (module 7) or EGFR signaling (module 15). The last mentioned might be supported by a threefold difference in the frequency of EGFR mutations between LUSC and LUAD during a previous pan-cancer analysis .
A separate topic was the division of breast cancer subtypes noticed during the global analysis of all cancer cohorts in relation to their corresponding normal tissues. Both target gene lists i.e. for AP-2α and AP-2γ, have separated two subgroups of BRCA samples. Additional spatial analysis and assignment of PAM50 signatures to samples revealed that the distinguished group was a basal-like subtype (Fig. 4). Due to the fact that this subtype is the most aggressive form of breast cancer , finding potential differences from the others could support clinical aspects in the future. Ontological analysis was performed on modules whose constituent genes demonstrated opposite expression profiles (basal-like vs. luminal A/B or HER2-enriched). Mutual part indicated through target genes lists of both AP-2α and AP-2γ shown that the processes differentiating the basal-like subtype concern regulation of programmed cell death or EGFR and TNF signaling pathways (Fig. 11), with the latter observed during previous enrichment analysis . Nevertheless, the cadherin and Wnt signaling pathways play a key role in subtype differentiation due to the number of genes described ontologically in module 5 (for AP-2γ targets’ comparisons), as noted previously .
On the whole, the analysis presented in this research indicate that many differentiating processes or pathways between tumor and normal tissue or between different tumor types/subtypes find their reference in current scientific reports. However, some of the obtained alterations not having regard in the literature, could serve as a novel preliminary point for the development of new anti-cancer therapies.
Collectively, our research indicates both novel and previously described differences between tumors or between a specific tumor and normal matched solid tissue. In addition, an unplanned division of samples from the breast cancer cohort was performed, differentiating the basal-like subtype from the others. These findings could have the clinical applications.
Analyses also proved that the regulation of gene expression by AP-2α and AP-2γ is very complex and involves many biological processes and signaling pathways. The role of AP-2 in carcinogenesis is ambiguous, but it is associated with the regulation of processes that underlie cancer hallmarks (proliferation, apoptosis, angiogenesis, cell growth), and the target genes of both transcription factors allow for further analysis of changes that occur in individual tumors. Finally, the differences in AP-2α and AP-2γ expression observed between tissue types (e.g. THCA and SKCM or KIRC/KIRP vs. normal kidney tissue) suggest their potential usefulness in diagnosing and treating specific cancers.
Nonetheless, we believe that our findings regarding AP-2α and AP-2γ functional genomics sheds new light on specific comparisons yet it is entirely possible that subsequent discoveries through in-depth data mining could have a further impact on the future of medical practice.
Availability of data and materials
TCGA gene expression profiles (level 3 RNA-seqV2, RSEM normalized) are publicly available in the Broad GDAC Firehose repository (https://gdac.broadinstitute.org/) wherein cohort abbreviations used in the present study are accession identifiers (data status of 28th Jan 2016). Direct access to respective cohort is possible through http://firebrowse.org/?cohort=X where “X” should be changed to cohort abbreviation (e.g. BLCA, BRCA, CESC, etc.).
Activating enhancer-binding Protein 2
WW Domain Containing Oxidoreductase
Estrogen receptor alpha
The Cancer Genome Atlas
Protein Analysis Through Evolutionary Relationships
Bladder urothelial carcinoma
Breast invasive carcinoma
Cervical and endocervical cancers
Head and neck squamous cell carcinoma
Kidney renal clear cell carcinoma
Kidney renal papillary cell carcinoma
Brain lower grade glioma
Liver hepatocellular carcinoma
Lung squamous cell carcinoma
Ovarian serous cystadenocarcinoma
Skin cutaneous melanoma
Uterine corpus endometrial carcinoma
Uniform Manifold Approximation and Projection
Prosigna Breast Cancer Prognostic Gene Signature Assay
Gene Transcription Regulation Database
TRANScription FACtor database
Transcriptional Regulatory Relationships Unraveled by Sentence-based Text mining
Principle Component Analysis
Nagai H, Kim YH. Cancer prevention from the perspective of global cancer burden patterns. J Thorac Dis. 2017;9(3):448–51.
Hyndman IJ. Review: the contribution of both nature and nurture to carcinogenesis and progression in solid tumours. Cancer Microenviron. 2016;9(1):63–9.
Jilkine A, Gutenkunst RN. Effect of dedifferentiation on time to mutation acquisition in stem cell-driven cancers. PLoS Comput Biol. 2014;10(3):e1003481.
Cortessis VK, Thomas DC, Levine AJ, Breton CV, Mack TM, Siegmund KD, et al. Environmental epigenetics: prospects for studying epigenetic mediation of exposure-response relationships. Hum Genet. 2012;131(10):1565–89.
Paduch R. Theories of cancer origin. Eur J Cancer Prev. 2015;24(1):57–67.
Sonnenschein C, Soto AM. Carcinogenesis explained within the context of a theory of organisms. Prog Biophys Mol Biol. 2016;122(1):70–6.
Zhu K, Liu Q, Zhou Y, Tao C, Zhao Z, Sun J, et al. Oncogenes and tumor suppressor genes: comparative genomics and network perspectives. BMC Genom. 2015;16(Suppl 7):S8.
Pedraza-Farina LG. Mechanisms of oncogenic cooperation in cancer initiation and metastasis. Yale J Biol Med. 2006;79(3–4):95–103.
Lou X, Zhang J, Liu S, Xu N, Liao DJ. The other side of the coin: the tumor-suppressive aspect of oncogenes and the oncogenic aspect of tumor-suppressive genes, such as those along the CCND-CDK4/6-RB axis. Cell Cycle. 2014;13(11):1677–93.
Muller PA, Vousden KH. Mutant p53 in cancer: new functions and therapeutic opportunities. Cancer Cell. 2014;25(3):304–17.
Ruiz M, Pettaway C, Song R, Stoeltzing O, Ellis L, Bar-Eli M. Activator protein 2alpha inhibits tumorigenicity and represses vascular endothelial growth factor transcription in prostate cancer cells. Cancer Res. 2004;64(2):631–8.
Shi D, Xie F, Zhang Y, Tian Y, Chen W, Fu L, et al. TFAP2A regulates nasopharyngeal carcinoma growth and survival by targeting HIF-1alpha signaling pathway. Cancer Prev Res (Phila). 2014;7(2):266–77.
Aqeilan RI, Palamarchuk A, Weigel RJ, Herrero JJ, Pekarsky Y, Croce CM. Physical and functional interactions between the Wwox tumor suppressor protein and the AP-2gamma transcription factor. Cancer Res. 2004;64(22):8256–61.
Li H, Goswami PC, Domann FE. AP-2gamma induces p21 expression, arrests cell cycle, and inhibits the tumor growth of human carcinoma cells. Neoplasia. 2006;8(7):568–77.
Kolat D, Kaluzinska Z, Bednarek AK, Pluciennik E. The biological characteristics of transcription factors AP-2alpha and AP-2gamma and their importance in various types of cancers. Biosci Rep. 2019;39(3):BSR20181928.
Eckert D, Buhl S, Weber S, Jager R, Schorle H. The AP-2 family of transcription factors. Genome Biol. 2005;6(13):246.
Mohibullah N, Donner A, Ippolito JA, Williams T. SELEX and missing phosphate contact analyses reveal flexibility within the AP-2[alpha] protein: DNA binding complex. Nucleic Acids Res. 1999;27(13):2760–9.
Mitchell PJ, Wang C, Tjian R. Positive and negative regulation of transcription in vitro: enhancer-binding protein AP-2 is inhibited by SV40 T antigen. Cell. 1987;50(6):847–61.
Hilger-Eversheim K, Moser M, Schorle H, Buettner R. Regulatory roles of AP-2 transcription factors in vertebrate development, apoptosis and cell-cycle control. Gene. 2000;260(1–2):1–12.
Wu HR, Zhang J. AP-2alpha expression in papillary thyroid carcinoma predicts tumor progression and poor prognosis. Cancer Manag Res. 2018;10:2615–25.
Anttila MA, Kellokoski JK, Moisio KI, Mitchell PJ, Saarikoski S, Syrjanen K, et al. Expression of transcription factor AP-2alpha predicts survival in epithelial ovarian cancer. Br J Cancer. 2000;82(12):1974–83.
Wang W, Lv L, Pan K, Zhang Y, Zhao JJ, Chen JG, et al. Reduced expression of transcription factor AP-2alpha is associated with gastric adenocarcinoma prognosis. PLoS ONE. 2011;6(9):e24897.
Perkins SM, Bales C, Vladislav T, Althouse S, Miller KD, Sandusky G, et al. TFAP2C expression in breast cancer: correlation with overall survival beyond 10 years of initial diagnosis. Breast Cancer Res Treat. 2015;152(3):519–31.
Gee JM, Eloranta JJ, Ibbitt JC, Robertson JF, Ellis IO, Williams T, et al. Overexpression of TFAP2C in invasive breast cancer correlates with a poorer response to anti-hormone therapy and reduced patient survival. J Pathol. 2009;217(1):32–41.
Wang X, Sun D, Tai J, Chen S, Yu M, Ren D, et al. TFAP2C promotes stemness and chemotherapeutic resistance in colorectal cancer via inactivating hippo signaling pathway. J Exp Clin Cancer Res. 2018;37(1):27.
Lin CY, Chao A, Wang TH, Lee LY, Yang LY, Tsai CL, et al. Nucleophosmin/B23 is a negative regulator of estrogen receptor alpha expression via AP2gamma in endometrial cancer cells. Oncotarget. 2016;7(37):60038–52.
Zhang P, Hou Q, Yue Q. MiR-204-5p/TFAP2A feedback loop positively regulates the proliferation, migration, invasion and EMT process in cervical cancer. Cancer Biomark. 2020;28:381–90.
Huang HX, Yang G, Yang Y, Yan J, Tang XY, Pan Q. TFAP2A is a novel regulator that modulates ferroptosis in gallbladder carcinoma cells via the Nrf2 signalling axis. Eur Rev Med Pharmacol Sci. 2020;24(9):4745–55.
Sliwa A, Kubiczak M, Szczerba A, Walkowiak G, Nowak-Markwitz E, Burczynska B, et al. Regulation of human chorionic gonadotropin beta subunit expression in ovarian cancer. BMC Cancer. 2019;19(1):746.
Pellikainen J, Kataja V, Ropponen K, Kellokoski J, Pietilainen T, Bohm J, et al. Reduced nuclear expression of transcription factor AP-2 associates with aggressive breast cancer. Clin Cancer Res. 2002;8(11):3487–95.
Su W, Xia J, Chen X, Xu M, Nie L, Chen N, et al. Ectopic expression of AP-2alpha transcription factor suppresses glioma progression. Int J Clin Exp Pathol. 2014;7(12):8666–74.
Hallberg AR, Vorrink SU, Hudachek DR, Cramer-Morales K, Milhem MM, Cornell RA, et al. Aberrant CpG methylation of the TFAP2A gene constitutes a mechanism for loss of TFAP2A expression in human metastatic melanoma. Epigenetics. 2014;9(12):1641–7.
Lian W, Zhang L, Yang L, Chen W. AP-2alpha reverses vincristine-induced multidrug resistance of SGC7901 gastric cancer cells by inhibiting the Notch pathway. Apoptosis. 2017;22(7):933–41.
Makhov PB, Golovine KV, Kutikov A, Canter DJ, Rybko VA, Roshchin DA, et al. Reversal of epigenetic silencing of AP-2alpha results in increased zinc uptake in DU-145 and LNCaP prostate cancer cells. Carcinogenesis. 2011;32(12):1773–81.
Li Q, Dashwood RH. Activator protein 2alpha associates with adenomatous polyposis coli/beta-catenin and Inhibits beta-catenin/T-cell factor transcriptional activity in colorectal cancer cells. J Biol Chem. 2004;279(44):45669–75.
Hu J, Tan SK, Lim MGL, Chang SH, Cui G, Liu S, et al. Identification of a Wells-Dawson polyoxometalate-based AP-2gamma inhibitor with pro-apoptotic activity. Biochem J. 2018;475(11):1965–77.
Li Z, Xu X, Luo M, Hao J, Zhao S, Yu W, et al. Activator protein-2beta promotes tumor growth and predicts poor prognosis in breast cancer. Cell Physiol Biochem. 2018;47(5):1925–35.
Fu X, Zhang H, Chen Z, Yang Z, Shi D, Liu T, et al. TFAP2B overexpression contributes to tumor growth and progression of thyroid cancer through the COX-2 signaling pathway. Cell Death Dis. 2019;10(6):397.
Fraune C, Harms L, Buscheck F, Hoflmayer D, Tsourlakis MC, Clauditz TS, et al. Upregulation of the transcription factor TFAP2D is associated with aggressive tumor phenotype in prostate cancer lacking the TMPRSS2:ERG fusion. Mol Med. 2020;26(1):24.
Hoshi R, Watanabe Y, Ishizuka Y, Hirano T, Nagasaki-Maeoka E, Yoshizawa S, et al. Depletion of TFAP2E attenuates adriamycin-mediated apoptosis in human neuroblastoma cells. Oncol Rep. 2017;37(4):2459–64.
Wei L, Jin Z, Yang S, Xu Y, Zhu Y, Ji Y. TCGA-assembler 2: software pipeline for retrieval and processing of TCGA/CPTAC data. Bioinformatics. 2018;34(9):1615–7.
Yevshin I, Sharipov R, Valeev T, Kel A, Kolpakov F. GTRD: a database of transcription factor binding sites identified by ChIP-seq experiments. Nucleic Acids Res. 2017;45(D1):D61–7.
Yevshin I, Sharipov R, Kolmykov S, Kondrakhin Y, Kolpakov F. GTRD: a database on gene transcription regulation-2019 update. Nucleic Acids Res. 2019;47(D1):D100–5.
Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. 2014;32(4):381–6.
Qin S, Shi X, Wang C, Jin P, Ma F. Transcription factor and miRNA interplays can manifest the survival of ccRCC patients. Cancers (Basel). 2019;11(11):1668.
Lommen K, Vaes N, Aarts MJ, van Roermund JG, Schouten LJ, Oosterwijk E, et al. Diagnostic DNA methylation biomarkers for renal cell carcinoma: a systematic review. Eur Urol Oncol. 2019. https://doi.org/10.1016/j.euo.2019.07.011.
Chen P, Zhang W, Chen Y, Zheng X, Yang D. Comprehensive analysis of aberrantly expressed long noncoding RNAs, microRNAs, and mRNAs associated with the competitive endogenous RNA network in cervical cancer. Mol Med Rep. 2020;22(1):405–15.
Shen L, Liu M, Liu W, Cui J, Li C. Bioinformatics analysis of RNA sequencing data reveals multiple key genes in uterine corpus endometrial carcinoma. Oncol Lett. 2018;15(1):205–12.
Li ZX, Zheng ZQ, Wei ZH, Zhang LL, Li F, Lin L, et al. Comprehensive characterization of the alternative splicing landscape in head and neck squamous cell carcinoma reveals novel events associated with tumorigenesis and the immune microenvironment. Theranostics. 2019;9(25):7648–65.
Huang D, Ding Y, Luo WM, Bender S, Qian CN, Kort E, et al. Inhibition of MAPK kinase signaling pathways suppressed renal cell carcinoma growth and angiogenesis in vivo. Cancer Res. 2008;68(1):81–8.
Hu M, Fu X, Si Z, Li C, Sun J, Du X, et al. Identification of differently expressed genes associated with prognosis and growth in colon adenocarcinoma based on integrated bioinformatics analysis. Front Genet. 2019;10:1245.
Zeng JH, Xiong DD, Pang YY, Zhang Y, Tang RX, Luo DZ, et al. Identification of molecular targets for esophageal carcinoma diagnosis using miRNA-seq and RNA-seq data from The Cancer Genome Atlas: a study of 187 cases. Oncotarget. 2017;8(22):35681–99.
Marimuthu A, Jacob HK, Jakharia A, Subbannayya Y, Keerthikumar S, Kashyap MK, et al. Gene expression profiling of gastric cancer. J Proteomics Bioinform. 2011;4(4):74–82.
Kulik G. Personalized prostate cancer therapy based on systems analysis of the apoptosis regulatory network. Asian J Androl. 2015;17(3):471–4.
Miller DR, Ingersoll MA, Lin MF. ErbB-2 signaling in advanced prostate cancer progression and potential therapy. Endocr Relat Cancer. 2019;26(4):R195–209.
Bartis D, Csongei V, Weich A, Kiss E, Barko S, Kovacs T, et al. Down-regulation of canonical and up-regulation of non-canonical Wnt signalling in the carcinogenic process of squamous cell lung carcinoma. PLoS ONE. 2013;8(3):e57393.
Mitra AP, Hansel DE, Cote RJ. Prognostic value of cell-cycle regulation biomarkers in bladder cancer. Semin Oncol. 2012;39(5):524–33.
Peng Y, Chen Y, Chen S, Wang J, Jiang C, Hou W, et al. JUND-dependent up-regulation of HMOX1 is associated with cisplatin resistance in muscle-invasive bladder cancer. J Biochem. 2020;168(1):73–82.
Madden E, Logue SE, Healy SJ, Manie S, Samali A. The role of the unfolded protein response in cancer progression: From oncogenesis to chemoresistance. Biol Cell. 2019;111(1):1–17.
Benard A, Zeestraten EC, Goossens-Beumer IJ, Putter H, van de Velde CJ, Hoon DS, et al. DNA methylation of apoptosis genes in rectal cancer predicts patient survival and tumor recurrence. Apoptosis. 2014;19(11):1581–93.
Yuanhua L, Pudong Q, Wei Z, Yuan W, Delin L, Yan Z, et al. TFAP2A induced KRT16 as an oncogene in lung adenocarcinoma via EMT. Int J Biol Sci. 2019;15(7):1419–28.
Fendler A, Bauer D, Busch J, Jung K, Wulf-Goldenberg A, Kunz S, et al. Inhibiting WNT and NOTCH in renal cancer stem cells and the implications for human patients. Nat Commun. 2020;11(1):929.
Zghair AN, Sinha DK, Kassim A, Alfaham M, Sharma AK. Differential Gene Expression of BRCA1, ERBB2 and TP53 biomarkers between human breast tissue and peripheral blood samples of breast cancer. Anticancer Agents Med Chem. 2016;16(4):519–25.
Tan W, Luo X, Li W, Zhong J, Cao J, Zhu S, et al. TNF-alpha is a potential therapeutic target to overcome sorafenib resistance in hepatocellular carcinoma. EBioMedicine. 2019;40:446–56.
Sau A, Cabrita MA, Pratt MAC. NF-kappaB at the crossroads of normal mammary gland biology and the pathogenesis and prevention of BRCA1-mutated breast cancer. Cancer Prev Res (Phila). 2018;11(2):69–80.
Ishimoto T, Miyake K, Nandi T, Yashiro M, Onishi N, Huang KK, et al. Activation of transforming growth factor beta 1 signaling in gastric cancer-associated fibroblasts increases their motility, via expression of rhomboid 5 homolog 2, and ability to induce invasiveness of gastric cancer cells. Gastroenterology. 2017;153(1):191.e16-204.e16.
Chen F, Zhang Y, Senbabaoglu Y, Ciriello G, Yang L, Reznik E, et al. Multilevel genomics-based taxonomy of renal cell carcinoma. Cell Rep. 2016;14(10):2476–89.
Sanchez-Vega F, Mina M, Armenia J, Chatila WK, Luna A, La KC, et al. Oncogenic signaling pathways in the cancer genome atlas. Cell. 2018;173(2):321.e10-337.e10.
Chen M, Liu X, Du J, Wang XJ, Xia L. Differentiated regulation of immune-response related genes between LUAD and LUSC subtypes of lung cancers. Oncotarget. 2017;8(1):133–44.
Xu Y, Geng R, Yuan F, Sun Q, Liu B, Chen Q. Identification of differentially expressed key genes between glioblastoma and low-grade glioma by bioinformatics analysis. PeerJ. 2019;7:e6560.
Liu H, Zhang B, Sun Z. Spectrum of EGFR aberrations and potential clinical implications: insights from integrative pan-cancer analysis. Cancer Commun (Lond). 2020;40(1):43–59.
Yang K, Gao J, Luo M. Identification of key pathways and hub genes in basal-like breast cancer using bioinformatics analysis. Onco Targets Ther. 2019;12:1319–31.
Jiang S, Zhang M, Zhang Y, Zhou W, Zhu T, Ruan Q, et al. WNT5B governs the phenotype of basal-like breast cancer by activating WNT signaling. Cell Commun Signal. 2019;17(1):109.
This research was funded by the Medical University of Lodz, grant number 503/0–078-02/503–01-001–19-00. The funding body had no role in study design, collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
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.
. Content of modules from particular Monocle3 approach.
. Detailed ontological analysis of selected modules differentiating tumor and corresponding normal tissue by means of AP-2α target genes.
. Detailed ontological analysis of selected modules differentiating tumor and corresponding normal tissue by means of AP-2γ target genes.
. Detailed ontological analysis of selected modules differentiating basal-like breast cancer from other subtypes by means of AP-2α and AP-2γ lists of target genes.
About this article
Cite this article
Kołat, D., Kałuzińska, Ż., Orzechowska, M. et al. Functional genomics of AP-2α and AP-2γ in cancers: in silico study. BMC Med Genomics 13, 174 (2020). https://doi.org/10.1186/s12920-020-00823-9
- Transcription factors
- Target genes