The role of Tcfap2c in tumorigenesis and cancer growth in an activated Neu model of mammary carcinogenesis

ANOS1 accelerates the progression of esophageal cancer identified by multi-omic approaches

To assess the role of ANOS1 in esophageal cancer (ESCA) progression, multi-omic analysis and experimental validation were employed. It was revealed that ANOS1 expression is significantly enhanced in ESCA patients and cell lines. The expression level of ANOS1 in ESCA patients can distinguish the malignancy from normal tissue with an area under curve (AUC) >0.75. Moreover, increased expression of ANOS1 is associated with advanced T stage and worse disease-free survival of ESCA patients. Therefore, a clinically applicable nomogram with ANOS1 was established with strong predictive power. Furthermore, high expression of ANOS1 in ESCA is correlated with (i) the enrichment of epithelial-mesenchymal transition by gene set enrichment analysis, (ii) the involvement in hypoxia, angiogenesis, WNT signaling pathway, and TGFβ signaling pathway by gene set variation analysis, (iii) the presence of the small insertion and deletion mutational signature ID9, associated with chromothripsis, in the single-nucleotide polymorphism analysis, (iv) the amplification of 11q13.3 in the copy number variants analysis, (v) the enrichment of cancer-associated fibroblasts and mesenchymal stromal cells in the tumor microenvironment. All the results from multi-omic analysis indicate that ANOS1 plays a pivotal role in accelerating the progression of ESCA. Results from in vivo and in vitro experiments show that the knockdown of ANOS1 hampers the proliferation of ESCA cells, further validating the oncogenic role of ANOS1 in ESCA. Additionally, potential chemotherapeutics with sensitivity were identified in the high-ANOS1 group. In conclusion, ANOS1 accelerates the progression of ESCA.

Decoding gene regulatory circuitry underlying TNBC chemoresistance reveals biomarkers for therapy response and therapeutic targets

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype characterised by extensive intratumoral heterogeneity, high rates of metastasis and chemoresistance, leading to poor clinical outcomes. Despite progress, the mechanistic basis of chemotherapy resistance in TNBC patients remains poorly understood. Here, leveraging single-cell transcriptome datasets of matched longitudinal TNBC chemoresponsive and chemoresistant patient cohorts, we unravel distinct cell subpopulations intricately associated with chemoresistance and the signature genes defining these populations. Notably, using genome-wide mapping of the H3K27ac mark, we show that the expression of these chemoresistance genes is driven via a set of TNBC super-enhancers and associated transcription factor networks across TNBC subtypes. Furthermore, genetic screens reveal that a subset of these transcription factors is essential for the survival of TNBC cells, and their loss increases sensitivity to chemotherapeutic agents. Overall, our study has revealed epigenetic and transcription factor networks underlying chemoresistance and suggests novel avenues to stratify and improve the treatment of patients with a high risk of developing resistance.

Tyrosine phosphorylation-mediated YAP1-TFAP2A interactions coordinate transcription and trastuzumab resistance in HER2+ breast cancer

Trastuzumab resistance in HER2+ breast cancer (BC) is the major reason leading to poor prognosis of BC patients. Oncogenic gene overexpression or aberrant activation of tyrosine kinase SRC is identified to be the key modulator of trastuzumab response. However, the detailed regulatory mechanisms underlying SRC activation-associated trastuzumab resistance remain poorly understood. In the present study, we discover that SRC-mediated YAP1 tyrosine phosphorylation facilitates its interaction with transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha, TFAP2A), which in turn promotes YAP1/TEAD-TFAP2A (YTT) complex-associated transcriptional outputs, thereby conferring trastuzumab resistance in HER2+ BC. Inhibition of SRC kinase activity or disruption of YTT complex sensitizes cells to trastuzumab treatment in vitro and in vivo. Additionally, we also identify YTT complex co-occupies the regulatory regions of a series of genes related to trastuzumab resistance and directly regulates their transcriptions, including EGFR, HER2, H19 and CTGF. Moreover, YTT-mediated transcriptional regulation is coordinated by SRC kinase activity. Taken together, our study reveals that SRC-mediated YTT complex formation and transcriptions are responsible for multiple mechanisms associated with trastuzumab resistance. Therefore, targeting HER2 signaling in combination with the inhibition of YTT-associated transcriptional outputs could serve as the treatment strategy to overcome trastuzumab resistance caused by SRC activation.

Structural basis for specific DNA sequence motif recognition by the TFAP2 transcription factors

The TFAP2 family regulates gene expression during differentiation, development, and organogenesis, and includes five homologs in humans. They all possess a highly conserved DNA binding domain (DBD) followed by a helix-span-helix (HSH) domain. The DBD-HSH tandem domain specifically binds to a GCC(N3)GGC consensus sequence, but the precise recognition mechanisms remain unclear. Here, we found that TFAP2 preferred binding to the GCC(N3)GGC sequence, and the pseudo-palindromic GCC and GGC motifs and the length of the central spacer between the two motifs determined their binding specificity. Structural studies revealed that the two flat amphipathic α-helical HSH domains of TFAP2A stacked with each other to form a dimer via hydrophobic interactions, while the stabilized loops from both DBD domains inserted into two neighboring major grooves of the DNA duplex to form base-specific interactions. This specific DNA binding mechanism controlled the length of the central spacer and determined the DNA sequence specificity of TFAP2. Mutations of the TFAP2 proteins are implicated in various diseases. We illustrated that reduction or disruption of the DNA binding ability of the TFAP2 proteins is the primary cause of TFAP2 mutation-associated diseases. Thus, our findings also offer valuable insights into the pathogenesis of disease-associated mutations in TFAP2 proteins.

Crucial role of the transcription factors family activator protein 2 in cancer: current clue and views

The transcription factor family activator protein 2 (TFAP2) is vital for regulating both embryonic and oncogenic development. The TFAP2 family consists of five DNA-binding proteins, including TFAP2A, TFAP2B, TFAP2C, TFAP2D and TFAP2E. The importance of TFAP2 in tumor biology is becoming more widely recognized. While TFAP2D is not well studied, here, we mainly focus on the other four TFAP2 members. As a transcription factor, TFAP2 regulates the downstream targets directly by binding to their regulatory region. In addition, the regulation of downstream targets by epigenetic modification, posttranslational regulation, and interaction with noncoding RNA have also been identified. According to the pathways in which the downstream targets are involved in, the regulatory effects of TFAP2 on tumorigenesis are generally summarized as follows: stemness and EMT, interaction between TFAP2 and tumor microenvironment, cell cycle and DNA damage repair, ER- and ERBB2-related signaling pathway, ferroptosis and therapeutic response. Moreover, the factors that affect TFAP2 expression in oncogenesis are also summarized. Here, we review and discuss the most recent studies on TFAP2 and its effects on carcinogenesis and regulatory mechanisms.

TFAP2C Knockdown Sensitizes Bladder Cancer Cells to Cisplatin Treatment via Regulation of EGFR and NF-κB

Simple Summary

Bladder cancer (BCa) is considered one of the most common neoplasms of the urology system. Cisplatin-based chemotherapy has been the primary treatment for patients with advanced or metastatic BCa. Nevertheless, cisplatin resistance often limits the treatment of bladder cancer. We expect to find approaches to improve the therapeutic efficacy of cisplatin in bladder cancer. In recent years, many studies have shown that transcription factor AP-2 gamma (TFAP2C) acts as a key player in cancer development and and its expression level is closely related to the sensitivity of tumors to cisplatin. Our study investigated whether TFAP2C affects the sensitivity of BCa cells to cisplatin and the possible mechanisms. We found that TFAP2C expression was significantly upregulated in most BCa tissues compared to adjacent normal tissues. The present study confirmed that TFAP2C knockdown enhanced the anti-tumor effects of cisplatin by decreasing cisplatin-induced activation levels of epidermal growth factor receptor (EGFR) and nuclear factor kappaB (NF-κB). Specifically, this study provides a novel approach to improve the efficacy of cisplatin.

Abstract

Cisplatin is the first-line chemotherapy for advanced or metastatic bladder cancer. Nevertheless, approximately half of patients with BCa are insensitive to cisplatin therapy or develop cisplatin resistance during the treatment process. Therefore, it is especially crucial to investigate ways to enhance the sensitivity of tumor cells to cisplatin. Transcription factor AP-2 gamma (TFAP2C) is involved in cancer development and chemotherapy sensitivity. However, its relationship with chemotherapy has not been studied in BCa. In this study, we aimed to investigate the therapeutic potential of TFAP2C in human BCa. Results based on TCGA (The Cancer Genome Atlas), GTEx (The Genotype-Tissue Expression) and GEO (Gene Expression Omnibus) data showed that TFAP2C expression was upregulated in BCa tissues and that its high expression was associated with poor prognosis. Meanwhile, we demonstrated the overexpression of TFAP2C in BCa clinical specimens. Subsequently, in vitro, we knocked down TFAP2C in BCa cells and found that TFAP2C knockdown further increased cell cycle arrest and apoptosis caused by cisplatin. In addition, the inhibitory effect of cisplatin on BCa cell migration and invasion was enhanced by TFAP2C knockdown. Our data indicated that cisplatin increased epidermal growth factor receptor (EGFR) and nuclear factor-kappaB (NF-κB) activation levels, but TFAP2C knockdown suppressed this effect. Finally, in vivo data further validated these findings. Our study showed that TFAP2C knockdown affected the activation levels of EGFR and NF-κB and enhanced the anti-tumor effects of cisplatin in vivo and in vitro. This provides a new direction to improve the efficacy of traditional cisplatin chemotherapy.

TRIM37 Augments AP-2γ Transcriptional Activity and Cellular Localization via K63-linked Ubiquitination to Drive Breast Cancer Progression

Activator Protein 2 gamma (AP-2γ) is a master transcription factor that plays a critical role in the development and progression of breast cancer. However, the underlying mechanism is still unclear. Herein, using a proteomics approach, we identified Tripartite motif-containing 37 (TRIM37) as a novel coactivator of AP-2γ-mediated transcription in breast cancer cells. We demonstrate that TRIM37 facilitates AP-2γ chromatin binding to directly regulate the AP-2γ mediated transcriptional program. We also show that TRIM37 achieves this by stimulating K63 chain-linked ubiquitination of AP-2γ, promoting protein localization from the cytoplasm to the nucleus. In clinical analyses, we find TRIM37 is upregulated in multiple breast cancer datasets, supporting our findings that the TRIM37-AP-2γ interaction is essential for breast cancer tumor growth. Overall, our work reveals that TRIM37 is an oncogenic coactivator of AP-2γ in breast cancer and provides a novel therapeutic target for treating the disease.

WWOX Loses the Ability to Regulate Oncogenic AP-2γ and Synergizes with Tumor Suppressor AP-2α in High-Grade Bladder Cancer

The cytogenic locus of the WWOX gene overlaps with the second most active fragile site, FRA16D, which is present at a higher frequency in bladder cancer (BLCA) patients with smoking habit, a known risk factor of this tumor. Recently, we demonstrated the relevance of the role of WWOX in grade 2 BLCA in collaboration with two AP-2 transcription factors whose molecular actions supported or opposed pro-cancerous events, suggesting a distinct character. As further research is needed on higher grades, the aim of the present study was to examine WWOX-AP-2 functionality in grade 3 and 4 BLCA using equivalent in vitro methodology with additional transcriptome profiling of cellular variants. WWOX and AP-2α demonstrated similar anti-cancer functionality in most biological processes with subtle differences in MMP-2/9 regulation; this contradicted that of AP-2γ, whose actions potentiated cancer progression. Simultaneous overexpression of WWOX and AP-2α/AP-2γ revealed that single discrepancies appear in WWOX-AP-2α collaboration but only at the highest BLCA grade; WWOX-AP-2α collaboration was considered anti-cancer. However, WWOX only appeared to have residual activity against oncogenic AP-2γ in grade 3 and 4: variants with either AP-2γ overexpression alone or combined WWOX and AP-2γ overexpression demonstrated similar pro-tumoral behavior. Transcriptome profiling with further gene ontology certified biological processes investigated in vitro and indicated groups of genes consisting of AP-2 targets and molecules worth investigation as biomarkers. In conclusion, tumor suppressor synergism between WWOX and AP-2α is unimpaired in high-grade BLCA compared to intermediate grade, yet the ability of WWOX to guide oncogenic AP-2γ is almost completely lost.

Interaction of transcription factor AP-2 gamma with proto-oncogene PELP1 promotes tumorigenesis by enhancing RET signaling

A significant proportion of estrogen receptor-positive (ER+) breast cancer (BC) initially responds to endocrine therapy but eventually evolves into therapy-resistant BC. Transcription factor AP-2 gamma (TFAP2C) is a known regulator of ER activity, and high expression of TFAP2C is associated with a decreased response to endocrine therapies. PELP1 is a nuclear receptor coregulator, commonly overexpressed in BC, and its levels are correlated with poorer survival. In this study, we identified PELP1 as a novel interacting protein of TFAP2C. RNA-seq analysis of PELP1 knockdown BC cells followed by transcription factor motif prediction pointed to TFAP2C being enriched in PELP1-regulated genes. Gene set enrichment analysis (GSEA) revealed that the TFAP2C-PELP1 axis induced a subset of common genes. Reporter gene assays confirmed PELP1 functions as a coactivator of TFAP2C. Mechanistic studies showed that PELP1-mediated changes in histone methylation contributed to increased expression of the TFAP2C target gene RET. Furthermore, the TFAP2C-PELP1 axis promoted the activation of the RET signaling pathway, which contributed to downstream activation of AKT and ERK pathways in ER+ BC cells. Concomitantly, knockdown of PELP1 attenuated these effects mediated by TFAP2C. Overexpression of TFAP2C contributed to increased cell proliferation and therapy resistance in ER+ BC models, while knockdown of PELP1 mitigated these effects. Utilizing ZR75-TFAP2C xenografts with or without PELP1 knockdown, we provided genetic evidence that endogenous PELP1 is essential for TFAP2C-driven BC progression in vivo. Collectively, our studies demonstrated that PELP1 plays a critical role in TFAP2C transcriptional and tumorigenic functions in BC and blocking the PELP1-TFAP2C axis could have utility for treating therapy resistance.

TFAP2C facilitates somatic cell reprogramming by inhibiting c-Myc-dependent apoptosis and promoting mesenchymal-to-epithelial transition

Transcription factors are known to mediate the conversion of somatic cells to induced pluripotent stem cells (iPSCs). Transcription factor TFAP2C plays important roles in the regulation of embryonic development and carcinogenesis; however, the roles of Tfap2c in regulating somatic cell reprogramming are not well understood. Here we demonstrate Tfap2c is induced during the generation of iPSCs from mouse fibroblasts and acts as a facilitator for iPSCs formation. Mechanistically, the c-Myc-dependent apoptosis, which is a roadblock to reprogramming, can be significantly mitigated by Tfap2c overexpression. Meanwhile, Tfap2c can greatly promote mesenchymal-to-epithelial transition (MET) at initiation stage of OSKM-induced reprogramming. Further analysis of gene expression and targets of Tfap2c during reprogramming by RNA-sequencing (RNA-seq) and ChIP-qPCR indicates that TFAP2C can promote epithelial gene expression by binding to their promoters directly. Finally, knockdown of E-cadherin (Cdh1), an important downstream target of TFAP2C and a critical regulator of MET antagonizes Tfap2c-mediated reprogramming. Taken together, we conclude that Tfap2c serves as a strong activator for somatic cell reprogramming through promoting the MET and inhibiting c-Myc-dependent apoptosis.

Transcriptomics-based screening of molecular signatures associated with patients overall survival and their key regulators in subtypes of breast cancer

Molecular subtypes of breast cancer are associated with differences in prognosis and strategies of molecular targeted therapies. Gene regulatory mechanisms as one of the reasons might modulate these differences. In the present study, we proposed a comprehensive data analysis and systems biology approach to explore molecular signatures which reduce the chance of patients overall survival and the possible mechanisms of their regulation by transcription factors (TFs) and microRNAs (miRNAs) in the main subtypes of breast tumor consist of Basal like, Her2 enriched, Luminal A and Luminal B breast cancer. In this regards, we used available microarray datasets to assess common differentially expressed genes (DEGs) in breast cancer subtypes. Using Kaplan-Meier curve analysis we identified common DEGs which are associated with decreasing in the overall survival of breast cancer patients. Furthermore, gene regulatory networks (GRNs) were depicted based on TFs and miRNAs with interest target genes. Then GRNs were analyzed and using five algorithms (Control centrality, Betweenness, Degree, Classification, and MCDS) the key regulators were identified for each subtype. In this study, we highlighted mechanisms underlying the regulation of breast cancer molecular signatures by TFs and miRNAs which their alteration reduce the chance of survival rate in each subtype of breast cancer. Our current study in a holistic insight revealed the importance of some genes and their regulators as potential prognostic markers and/or therapeutic targets in breast cancer patients.

The biological characteristics of transcription factors AP-2α and AP-2γ and their importance in various types of cancers

The Activator Protein 2 (AP-2) transcription factor (TF) family is vital for the regulation of gene expression during early development as well as carcinogenesis process. The review focusses on the AP-2α and AP-2γ proteins and their dualistic regulation of gene expression in the process of carcinogenesis. Both AP-2α and AP-2γ influence a wide range of physiological or pathological processes by regulating different pathways and interacting with diverse molecules, i.e. other proteins, long non-coding RNAs (lncRNA) or miRNAs. This review summarizes the newest information about the biology of two, AP-2α and AP-2γ, TFs in the carcinogenesis process. We emphasize that these two proteins could have either oncogenic or suppressive characteristics depending on the type of cancer tissue or their interaction with specific molecules. They have also been found to contribute to resistance and sensitivity to chemotherapy in oncological patients. A better understanding of molecular network of AP-2 factors and other molecules may clarify the atypical molecular mechanisms occurring during carcinogenesis, and may assist in the recognition of new diagnostic biomarkers.

Prospective Isolation of Poised iPSC Intermediates Reveals Principles of Cellular Reprogramming

Cellular reprogramming converts differentiated cells into induced pluripotent stem cells (iPSCs). However, this process is typically very inefficient, complicating mechanistic studies. We identified and molecularly characterized rare, early intermediates poised to reprogram with up to 95% efficiency, without perturbing additional genes or pathways, during iPSC generation from mouse embryonic fibroblasts. Analysis of these cells uncovered transcription factors (e.g., Tfap2c and Bex2) that are important for reprogramming but dispensable for pluripotency maintenance. Additionally, we observed striking patterns of chromatin hyperaccessibility at pluripotency loci, which preceded gene expression in poised intermediates. Finally, inspection of these hyperaccessible regions revealed an early wave of DNA demethylation that is uncoupled from de novo methylation of somatic regions late in reprogramming. Our study underscores the importance of investigating rare intermediates poised to produce iPSCs, provides insights into reprogramming mechanisms, and offers a valuable resource for the dissection of transcriptional and epigenetic dynamics intrinsic to cell fate change.

Identification of a Wells-Dawson polyoxometalate-based AP-2γ inhibitor with pro-apoptotic activity

AP-2 gamma (AP-2γ) is a transcription factor that plays pivotal roles in breast cancer biology. To search for small molecule inhibitors of AP-2γ, we performed a high-throughput fluorescence anisotropy screen and identified a polyoxometalate compound with Wells-Dawson structure K₆[P₂Mo₁₈O₆₂] (Dawson-POM) that blocks the DNA-binding activity of AP-2γ. We showed that this blocking activity is due to the direct binding of Dawson-POM to AP-2γ. We also provided evidence to show that Dawson-POM decreases AP-2γ-dependent transcription similar to silencing the gene. Finally, we demonstrated that Dawson-POM contains anti-proliferative and pro-apoptotic effects in breast cancer cells. In summary, we identified the first small molecule inhibitor of AP-2γ and showed Dawson-POM-mediated inhibition of AP-2γ as a potential avenue for cancer therapy.

Upregulation of microRNA-137 expression by Slug promotes tumor invasion and metastasis of non-small cell lung cancer cells through suppression of TFAP2C

The epithelial-mesenchymal transition (EMT) regulator, Slug, plays multifaceted roles in controlling lung cancer progression, but its downstream targets and mechanisms in promoting lung cancer progression have not been well defined. In particular, the miRNAs downstream of Slug in non-small cell lung cancer (NSCLC) remain undetermined. Here, we report that miR-137 is downstream of the EMT regulator, Slug, in lung cancer cells. Slug binds directly to the E-box of the miR-137 promoter and up-regulates its expression in lung cancer cells. Knockdown of miR-137 abolished Slug-induced cancer invasion and migration, whereas upregulation of miR-137 was found to trigger lung cancer cell invasion and progression by direct suppressing TFAP2C (transcription factor AP-2 gamma). Clinical data showed that lung adenocarcinoma patients with low-level expression of Slug and miR-137 but high-level expression of TFAP2C experienced significantly better survival. miR-137 is a Slug-induced miRNA that relays the pro-metastatic effects of Slug by targeting TFAP2C. Our findings add new components to the Slug-mediated regulatory network in lung cancer, and suggest that Slug, miR-137, and TFAP2C may be useful prognostic markers in lung adenocarcinoma.

Identification of transcription factors that may reprogram lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma is one of most threatening disease to human health. Although many efforts have been devoted to its genetic study, few researches have been focused on the transcription factors which regulate tumor initiation and progression by affecting multiple downstream gene transcription. It is proved that proper transcription factors may mediate the direct reprogramming of cancer cells, and reverse the tumorigenesis on the epigenetic and transcription levels.

METHODS

In this paper, a computational method is proposed to identify the core transcription factors that can regulate as many as possible lung adenocarcinoma associated genes with as little as possible redundancy. A greedy strategy is applied to find the smallest collection of transcription factors that can cover the differentially expressed genes by its downstream targets. The optimal subset which is mostly enriched in the differentially expressed genes is then selected.

RESULTS

Seven core transcription factors (MCM4, VWF, ECT2, RBMS3, LIMCH1, MYBL2 and FBXL7) are detected, and have been reported to contribute to tumorigenesis of lung adenocarcinoma. The identification of the transcription factors provides a new insight into its oncogenic role in tumor initiation and progression, and benefits the discovery of functional core set that may reverse malignant transformation and reprogram cancer cells.

Genome-wide profiling of long non-coding RNA expression patterns in the EGFR-TKI resistance of lung adenocarcinoma by microarray

Mutations in the epidermal growth factor receptor (EGFR) make lung adenocarcinoma cells sensitive to EGFR tyrosine kinase inhibitors (TKIs). Long-term cancer therapy may cause the occurrence of acquired resistance to EGFR TKIs. Long non-coding RNAs (lncRNAs) play important roles in tumor formation, tumor metastasis and the development of EGFR-TKI resistance in lung cancer. To gain insight into the molecular mechanisms of EGFR-TKI resistance, we generated an EGFR-TKI-resistant HCC827-8-1 cell line and analyzed expression patterns by lncRNA microarray and compared it with its parental HCC827 cell line. A total of 1,476 lncRNA transcripts and 1,026 mRNA transcripts were dysregulated in the HCC827‑8-1 cells. The expression levels of 7 chosen lncRNAs were validated by real-time quantitative PCR. As indicated by functional analysis, several groups of lncRNAs may be involved in the bio-pathways associated with EGFR-TKI resistance through their cis- and/or trans‑regulation of protein-coding genes. Thus, lncRNAs may be used as novel candidate biomarkers and potential targets in EGFR-TKI therapy in the future.

EGFR Is Regulated by TFAP2C in Luminal Breast Cancer and Is a Target for Vandetanib

Expression of TFAP2C in luminal breast cancer is associated with reduced survival and hormone resistance, partially explained through regulation of RET. TFAP2C also regulates EGFR in HER2 breast cancer. We sought to elucidate the regulation and functional role of EGFR in luminal breast cancer. We used gene knockdown (KD) and treatment with a tyrosine kinase inhibitor (TKI) in cell lines and primary cancer isolates to determine the role of RET and EGFR in regulation of p-ERK and tumorigenesis. KD of TFAP2C decreased expression of EGFR in a panel of luminal breast cancers, and chromatin immunoprecipitation sequencing (ChIP-seq) confirmed that TFAP2C targets the EGFR gene. Stable KD of TFAP2C significantly decreased cell proliferation and tumor growth, mediated in part through EGFR. While KD of RET or EGFR reduced proliferation (31% and 34%, P < 0.01), combined KD reduced proliferation greater than either alone (52% reduction, P < 0.01). The effect of the TKI vandetanib on proliferation and tumor growth response of MCF-7 cells was dependent upon expression of TFAP2C, and dual KD of RET and EGFR eliminated the effects of vandetanib. The response of primary luminal breast cancers to TKIs assessed by ERK activation established a correlation with expression of RET and EGFR. We conclude that TFAP2C regulates EGFR in luminal breast cancer. Response to vandetanib was mediated through the TFAP2C target genes EGFR and RET. Vandetanib may provide a therapeutic effect in luminal breast cancer, and RET and EGFR can serve as molecular markers for response.

Does mouse embryo primordial germ cell activation start before implantation as suggested by single-cell transcriptomics dynamics?

STUDY HYPOTHESIS

Does primordial germ cell (PGC) activation start before mouse embryo implantation, and does the possible regulation of the DNA (cytosine-5-)-methyltransferase 3-like (Dnmt3l) by transcription factor AP-2, gamma (TCFAP2C) have a role in this activation and in the primitive endoderm (PE)-epiblast (EPI) lineage specification?

STUDY FINDING

A burst of expression of PGC markers, such as Dppa3/Stella, Ifitm2/Fragilis, Fkbp6 and Prdm4, is observed from embryonic day (E) 3.25, and some of them, together with the late germ cell markers Zp3, Mcf2 and Morc1, become restricted to the EPI subpopulation at E4.5, while the dynamics analysis of the PE-EPI transitions in the single-cell data suggests that TCFAP2C transitorily represses Dnmt3l in EPI cells at E3.5 and such repression is withdrawn with reactivation of Dnmt3l expression in PE and EPI cells at E4.5.

WHAT IS KNOWN ALREADY

In the mouse preimplantation embryo, cells with the same phenotype take different fates based on the orchestration between topological clues (cell polarity, positional history and division orientation) and gene regulatory rules (at transcriptomics and epigenomics level), prompting the proposal of positional, stochastic and combined models explaining the specification mechanism. PGC specification starts at E6.0-6.5 post-implantation. In view of the important role of DNA methylation in developmental events, the cross-talk between some transcription factors and DNA methyltransferases is of particular relevance. TCFAP2C has a CpG DNA methylation motif that is not methylated in pluripotent cells and that could potentially bind on DNMT3L, the stimulatory DNA methyltransferase co-factor that assists in the process of de novo DNA methylation. Chromatin-immunoprecipitation analysis has demonstrated that Dnmt3l is indeed a target of TCFAP2C.

STUDY DESIGN, SAMPLES/MATERIALS, METHODS

We aimed to assess the timing of early preimplantation events and to understand better the segregation of the inner cell mass (ICM) into PE and EPI. We designed a single-cell transcriptomics dynamics computational study to identify markers of the PE-EPI bifurcation in ICM cells through searching for statistically significant (using the Student's t-test method) differently expressed genes (DEGs) between PE and EPI cells from E3.5 to E4.5. The DEGs common for E3.5 and E4.5 were used as the markers defining the steady states. We collected microarray and next-generation sequencing transcriptomics data from public databases from bulk populations and single cells from mice at E3.25, E3.5 and E4.5. The results are based on three independent single-cell transcriptomics data sets, with a fold change of 3 and P-value <0.01 for the DEG selection.

MAIN RESULTS AND THE ROLE OF CHANCE

The dynamics analysis revealed new transitory E3.5 and steady PE and EPI markers. Among the transitory E3.5 PE markers (Dnmt3l, Dusp4, Cpne8, Akap13, Dcaf12l1, Aaed1, B4galt6, BC100530, Rnpc3, Tfpi, Lgalsl, Ckap4 and Fbxl20), several (Dusp4, Akap13, Cpn8, Dcaf12l1 and Tfpi) are related to the extracellular regulated kinase pathway. We also identified new transitory E3.5 EPI markers (Sgk1, Mal, Ubxn2a, Atg16l2, Gm13102, Tcfap2c, Hexb, Slc1a1, Svip, Liph and Mier3), six new stable PE markers (Sdc4, Cpn1, Dkk1, Havcr1, F2r/Par1 and Slc7a6os) as well as three new stable EPI markers (Zp3, Mcf2 and Hexb), which are known to be late stage germ cell markers. We found that mouse PGC marker activation starts at least at E3.25 preimplantation. The transcriptomics dynamics analyses support the regulation of Dnmt3l expression by TCFAP2C.

LIMITATIONS, REASONS FOR CAUTION

Since the regulation of Dnmt3l by TCFAP2C is based on computational prediction of DNA methylation motifs, Chip-Seq and transcriptomics data, functional studies are required to validate this result.

WIDER IMPLICATIONS OF THE FINDINGS

We identified a collection of previously undescribed E3.5-specific PE and EPI markers, and new steady PE and EPI markers. Identification of these genes, many of which encode cell membrane proteins, will facilitate the isolation and characterization of early PE and EPI populations. Since it is so well established in the literature that mouse PGC specification is a post-implantation event, it was surprising for us to see activation of PGC markers as early as E3.25 preimplantation, and identify the newly found steady EPI markers as late germ cell markers. The discovery of such early activation of PGC markers has important implications in the derivation of germ cells from pluripotent cells (embryonic stem cells or induced pluripotent stem cells), since the initial stages of such derivation resemble early development. The early activation of PGC markers points out the difficulty of separating PGC cells from pluripotent populations. Collectively, our results suggest that the combining of the precision of single-cell omics data with dynamic analysis of time-series data can establish the timing of some developmental stages as earlier than previously thought.

LARGE-SCALE DATA

Not applicable.

STUDY FUNDING AND COMPETING INTERESTS

This work was supported by grants DFG15/14 and DFG15/020 from Diputación Foral de Gipuzkoa (Spain), and grant II14/00016 from I + D + I National Plan 2013-2016 (Spain) and FEDER funds. The authors declare no conflict of interest.