Unravelling Structure, Localization, and Genetic Crosstalk of KLF3 in Human Breast Cancer

Khan, Khushbukhat; Safi, Sadia; Abbas, Asma; Badshah, Yasmin; Dilshad, Erum; Rafiq, Mehak; Zahra, Kainat; Shabbir, Maria

doi:https://doi.org/10.1155/2020/1354381

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 1354381 | https://doi.org/10.1155/2020/1354381

Unravelling Structure, Localization, and Genetic Crosstalk of KLF3 in Human Breast Cancer

Khushbukhat Khan,¹Sadia Safi,¹Asma Abbas,¹Yasmin Badshah,¹Erum Dilshad,²Mehak Rafiq,³Kainat Zahra,¹and Maria Shabbir¹

Academic Editor: Claudio Acu a Castillo

Received01 Jul 2020

Revised08 Dec 2020

Accepted17 Dec 2020

Published28 Dec 2020

Abstract

Breast cancer is the most prevailing disease among women. It actually develops from breast tissue and has heterogeneous and complex nature that constitutes multiple tumor quiddities. These features are associated with different histological forms, distinctive biological characteristics, and clinical patterns. The predisposition of breast cancer has been attributed to a number of genetic factors, associated with the worst outcomes. Unfortunately, their behavior with relevance to clinical significance remained poorly understood. So, there is a need to further explore the nature of the disease at the transcriptome level. The focus of this study was to explore the influence of Krüppel-like factor 3 (KLF3), tumor protein D52 (TPD52), microRNA 124 (miR-124), and protein kinase C epsilon (PKCε) expression on breast cancer. Moreover, this study was also aimed at predicting the tertiary structure of KLF3 protein. Expression of genes was analyzed through real-time PCR using the delta cycle threshold method, and statistical significance was calculated by two-way ANOVA in Graphpad Prism. For the construction of a 3D model, various bioinformatics software programs, Swiss Model and UCSF Chimera, were employed. The expression of KLF3, miR-124, and PKCε genes was decreased (fold change: 0.076443, 0.06969, and 0.011597, respectively). However, there was 2-fold increased expression of TPD52 with value < 0.001 relative to control. Tertiary structure of KLF3 exhibited 80.72% structure conservation with its template KLF4 and was 95.06% structurally favored by a Ramachandran plot. These genes might be predictors of stage, metastasis, receptor, and treatment status and used as new biomarkers for breast cancer diagnosis. However, extensive investigations at the tissue level and in in vivo are required to further strengthen their role as a potential biomarker for prognosis of breast cancer.

1. Introduction

Breast cancer is a malignancy that originates from breast tissue. It comprises a mass of cells that along with invading surrounding tissues of the body also metastasize to distant parts such as the lungs, bone marrow, regional lymph nodes, and liver [1]. It is presumed that breast cancer results from the accumulation of several genetic aberrations, resulting in inactivation of tumor suppressor genes and activation of protooncogenes. Current evidences from literature indicate that approximately 10-20% of cases are due to a germline mutation in a gene that leads to increased predisposition for breast cancer [2, 3] while 15-20% of breast cancer incidences are familial [4]. Other than germ line genes, dysregulation in several other genes encoding tumor suppressor proteins, transcription regulators, microRNAs, kinases, and phosphatases has led to breast cancer [5, 6].

Genome-wide studies of human cancer have shown that Receptor Tyrosine Kinase (RTK) and G Protein-Coupled Receptor (GPCR) are the major receptors that regulate the effector protein [7, 8]. Phosphoinositide 3-kinases (PI3K) and Kras via channeling signal transduction from various growth factors and cytokines into intracellular signaling lead to the major activation of Akt and downstream pathways that will ultimately result in cancer [9]. Through a literature review, it is also revealed that KLF3, TPD52, PKCε, and miR-124 interact with different downstream signaling molecules of these pathways in different cancers and lead to cancer proliferation, stage progression and metastasis, and drug resistance [10–15].

KLF3 is a DNA-binding transcription repressor that belongs to the seventeen-member family of KLFs. It contains three C₂H₂ zinc finger elements at the C-terminal region. It binds to GC-rich elements and CACCC-binding motif at specific DNA through these motifs and regulates erythropoiesis, lymphopoiesis (particularly B lymphocyte), and adipogenesis [16, 17]. KLF3 engages with corepressor C-terminal binding protein (CtBP) and recruits other repressive cofactors, like histone deacetylases (HDAC), histone methyltransferases (HMT), and histone lysine-specific demethylases (HLSD) [18]. KLF3 acts as an oncogene in cervical and lung tumors [11] while in acute myeloid leukemia and pancreatic and prostate cancer, it displays tumor-suppressive behavior [19].

miR-124 is one of the widely expressed microRNAs in the nervous system, which is more abundant in differentiating neurons of the brain, spinal cord, and retina and lasts in mature neurons [20]. Previous reports had highlighted its role as oncomiRs and oncosuppressors in various cancers. Its dysregulation is associated with aberrant signal transduction in various cellular pathways in cancers such as STAT3 in glioblastoma [21], phosphoinositide 3-kinase (PI3K) and protein kinase B (PKB) (also known as Akt) in prostate cancer, and Notch signaling pathway in gastric cancer [22]. TPD52 is a 180 amino acid-long, highly charged acidic protein. Its genetic sequence is located on chromosome 8 (8q21). It is actually a peripheral membrane and cytosolic protein of the mammalian gene family that performs several biological functions including vesicle trafficking and lipid metabolism [23]. Its aberrant expression in several cancers such as prostate, ovarian, and colorectal cancer, adenocarcinoma, hepatocellular carcinoma, and renal and adrenocortical carcinoma had been widely reported [12].

PKCε is one of the isoforms of PKC which is the family of serine/threonine kinases [24]. The gene of PKCε (PRKCE) resides on chromosome 2 (2p21). It needs diacylglycerol/phobol esters and Ca²⁺ for its activity [25]. PKCε is involved actively in adhesion, migration, transportation, proliferation, differentiation, gene expression, inflammation, immune regulation, secretory processes, and various signal transduction processes [26]. PKCε involvement in the onset of various tumors such glioblastoma; renal, prostate, colorectal, and follicular thyroid cancer; and non-small-cell lung carcinoma is also reported [27]. To date, no study is conducted that investigated the coexpression profile of TPD52, KLF3, miR-124, and PKCε. So, the goal of the current study is to analyze the expression pattern of TPD52, KLF3, miR-124, and PKCε in breast cancer patients as well as to predict their possible correlation with clinicopathological characteristics of breast cancer. Furthermore, the three-dimensional structure of KLF3 is not predicted yet. In this study, we have also predicted the tertiary structure of KLF3 by Swiss modeling. Understanding the expression profile of these genes will provide further insight into complex molecular dynamics. Additionally, this data can be a step towards constructing a molecular diagnostic strategy for breast cancer.

2. Materials and Methods

2.1. In Vitro Method

2.1.1. Blood Sample Collection

After approval from the Institutional Review Board of Atta-ur-Rahman School of Applied Biosciences (Ref. no: IRB-110) and Combined Military Hospital, blood samples were collected from 50 breast cancer patients as well as 50 controls. All the procedures were performed according to the guidelines provided by ethical review board. A consent form was signed by patients who were willing to give their blood.

History forms of patients were filled after collecting all the relevant information which comprised patient’s age, tumor grade, cancer type, treatment status (pretreated and naïve-treated), and receptor subtypes. Other than pathological features, information regarding family history, breastfeeding, post/premenopausal age, and smoking/alcohol consumption was also collected from patients.

2.1.2. RNA Isolation and Synthesis of cDNA

RNA was extracted from 50 collected blood samples of breast cancer patients as well as 50 controls, i.e., healthy individuals, using the TRIzol reagent manufactured by Thermo Fisher Scientific (Waltham, MA, USA). cDNA was synthesized from the extracted mRNA using MLV reverse transcriptase (Promega). The synthesized cDNA was stored at -20°C.

2.1.3. Real-Time PCR

Expression analysis was performed by a real-time PCR (qPCR) quantification assay (Applied Biosystems 7300) using SYBR Premix ExTaq II (TaKaRa). The relative expression of TPD52, KLF3, miR-124, and PKCε was calculated using the Livak method () while the beta-actin gene was used as the housekeeping gene.

Primers used for understudied genes are as follows: KLF3 forward 5-TGTCTCAGTGTCATACCCATCT-3 and reverse 5-CCTTCTGGGGTCTGAAAGAACTT-3, miR-124 forward 5-GATACTCATAAGGCACGCGG-3 and reverse 5-GTGCAGGGTCCGAGGT-3, TPD52 forward 5-GAGGAAGGAGAAGATGTTGC-3 and reverse 5-GCCGAATTCAAGACTTCTCC-3, and PKCε forward 5-AGCCTCGTTCACGGTTCT-3 and reverse 5-TGTCCAGCCATCATCTCG-3.

2.1.4. Statistical Analysis

One-way and two-way ANOVA was performed via Graphpad Prism 6.0 for analysis of association between relative expression of TPD52, KLF3, miR-124, and PKCε and clinicopathological characteristics of breast cancer.

2.1.5. In Silico Method

The 345-amino acid sequence of KLF3 (Homo sapiens) in the FASTA format was retrieved from the National Center for Biotechnology Information (NCBI) having reference sequence id NP_057615.3. The multiple sequence alignment tool “Clustal Omega” was used for aligning homologous primary protein sequences of KLF family members by progressive alignments [28]. The alignment revealed conserved domains within the family. Various different tools were used to predict secondary structure such as Porter 5.0 [29], SOPMA [30], and I-TASSER [31]. For better understanding of evolutionary conservation, phylogenetic analysis of whole KLF family proteins was performed by using Molecular Evolutionary Genetics Analysis version 7.0 (MEGA7) software [32]. For the construction of 3D protein structure, KLF4 was selected as a template (the only predicted structure of the KLF family). The target sequence (KLF3) in the FASTA format and the PDB file of template (KLF4) were added and run in Swiss Model. The obtained 3D structure was visualized in UCSF Chimera 1.13 using template structure as the reference.

3. Results

3.1. Clinicopathological Features of Breast Cancer and Patient Distribution

After written and oral contest, blood samples of 50 breast cancer patients were collected. Patients were categorized according to clinical features of the diseases and risk factors (see Tables 1 and 2). Clinical features of the disease included receptor status, cancer stage and type, and treatment status while risk factors included were family history, age, breast feeding status, pre- or postmenopause, and smoking or alcohol consumption habits.

3.2. Subcellular Localization of KLF3

The subcellular localization of the KLF3 (NP_057615.3) protein is within “nucleus” (see Figure 1) as predicted by the server DeepLoc-1.0 that covers various organelles and also differentiates between them for localizations. The location of KLF3 is then checked and confirmed through the UniProtKB database.

3.3. Conserved Domains

The result summary from Clustal Omega has revealed that within the amino acid sequence of KLF3 protein, the areas that link three zinc fingers are zinc linker motifs with sequence: TGEKP and TGIKP. The amino acids at position 61-65 make the CtBP-binding motif while the motifs PVKQE and IKIE are involved in repression activity of KLF3 (see Figure 2).

Among all members of KLFs exist three regions that correspond to zinc finger domains with C₂H₂-type elements showing coordination with three zinc ions (see Figure 3). These all domains have regulatory function via interacting with promoters of targeted genes. The presence of these conserved domains was also confirmed through UniProtKB.

Figure 3

Sequence alignment of Krüppel-like factors depicting conserved domains obtained from Clustal Omega. Certain sequence alignments have to be deleted during formatting. Zinc figure domain 1 (labeled as green), zinc figure domain 2 (labeled as purple), and zinc figure domain 3 (labeled as blue) have been highlighted with underlined C₂H₂ elements that coordinate with zinc ions. Different residues of zinc fingers are involved in making coils predicted by I-TASSER, Porter 5.0, and SOPMA. In between the three zinc fingers, two linker motifs are present for link fingers (highlighted in red).

3.4. Phylogenetic Analysis

The results of evolutionary conservation between KLF members in the form of phylogeny had placed KLF3 in group 2 based on its transcriptional activity and domain difference from other members. This division among members was confirmed by Molecular Evolutionary Genetics Analysis version 7.0 (MEGA7) software. The software relates members according to evolution and organizes them into three groups in the phylogenetic tree (see Figure 4).

Figure 4

Phylogenetic analysis of KLF proteins by MEGA7. The evolutionary history was predicted by using the hierarchical clustering method unweighted pair group method with arithmetic mean (UPGMA). The evolutionary tree is calibrated to scale with the sum branch length of 7.54632578. In the optimal tree, the branch lengths are of the same unit that are used for evolutionary distances and computed by Poisson correction method by considering the number of amino acid substitutions per site. A total of 602 positions were present in the final amino acid sequences after removing ambiguous positions (pairwise deletion), and final evolutionary analyses indicating relationship of taxa were conducted in MEGA X [33].

3.5. 3-Dimensional Structure and Visualization of KLF3 Protein

The homology modeling between KLF3 and template KLF4 by using the Swiss Model server revealed 80.72% structural conservation. Along with sequence identity, the server also takes solvation factor, torsion angle, and structure quality estimate of KLF3 protein. The 3D structure of KLF3 is given in Figure 5. The quality estimate of the KLF3 model was observed by a Ramachandran plot and MolProbity scores which show that KLF3 protein structure is 95.06% favored by a Ramachandran plot; i.e., the predicted model was appropriate (see Figure 6). The built model of KLF3 was visualized in UCSF Chimera 1.13 for finding structural similarity between KLF4 (template) and the model (KLF3). The comparison of structure was evaluated by superimposition of the model (red) on the template (blue) as shown in Figure 7. The result of superimposition revealed that there exist 80.72% homology between both proteins after matching.

Figure 6

Ramachandran plot with -180 to +180 axes for determining secondary structures. The disfavored and allowed values of dihedral angles ( and ) are shown in a two-dimensional plot. In quadrant I (Q-I), some left-handed helices are allowed (KLF3 lacks loops), quadrant II (Q-II) contained amino acids that favor beta strand conformation, quadrant III (Q-III) favors the region where right-handed alpha helices lie, and quadrant IV (Q-IV) has two amino acids (THR and ASP) clashing with the protein structure. Overall, this plot showed 724 interactions (bonds) within residues and 970 angles with 95.06% structure favorability.

3.6. Expression of TPD52, KLF3, miR-124, and PKCε in Blood of Breast Cancer Patients

Expression level of TPD52 was found to be elevated while KLF3, miR-124, and PKCε genes were downregulated in breast cancer patients in comparison to controls. By using two-way analysis of variance (ANOVA), statistical significance was computed between both groups. The results indicated high expression of TPD52 () (Figure 8(a)) and low expression of KLF3 () (Figure 8(b)), miR-124 () (Figure 8(c)), and PKCε () (Figure 8(d)) in blood of breast cancer patients as compared to healthy controls. In general, there was 2-fold increased expression of TPD52, and 0.076443-, 0.06969-, and 0.011597-fold decrease in the expression of KLF3, miR-124, and PKCε, respectively, was observed in patients in comparison with healthy controls.

(a)

(b)

(c)

(d)

3.7. Relative Expression of TPD52, KLF3, miR-124, and PKCε in Clinical-Pathological Features of Breast Cancer

The expression of TPD52, KLF3, miR-124, and PKCε in breast cancer patients was measured (see Supplementary Table 1). The significant results of expression were found between tumor grades I/II and III/IV, metastatic vs. nonmetastatic features, pretreated vs. nontreated patients, and breast cancer subtypes (luminal A, luminal B, and triple negative). In the low-grade cancer group (grades I and II), the expression of TPD52 was found to be significantly elevated () as compared to that in the advance stage (III and IV). Moreover, the expression of TPD52 was significantly upregulated in nonmetastatic and metastatic groups (). Among ER/PR/Her2-neu receptors, the most enhanced expression of TPD52 was present in the luminal B subtype group and the outcome was statistically significant (). The patients that were naïve to any treatment (chemotherapy/radiation) carried advance expression () as compared to controls (see Figure 9).

(a)

(b)

(c)

(d)

Figure 9

qPCR analysis of TPD52 expression in clinical characteristics of breast cancer cases. Correlation of TPD52 expression with (a) tumor stage, (b) metastasis, (c) receptor status, and (d) treatment condition. There was upregulated expression of TPD52 in lower tumor stages (I and II), nonmetastasis, cancer subtype (luminal B), and naïve-treated groups. The data expressed as fold change represents the experiments performed in triplicate. Ordinary one-way ANOVA was used to establish significance ().

As indicated in Figure 10, the expression of KLF3 was downregulated in stage I-IV, metastatic vs. nonmetastatic groups, in every molecular subtype (luminal A, luminal B, and triple negative), and in patients that were either on treatment or naïve to any kind of breast cancer treatment as compared to healthy controls. KLF3 expression was lower in the luminal B group in comparison to luminal A and triple negative groups, and its expression was further downregulated in the metastasis group in comparison to the nonmetastasis group. The expression of miR-124 was also downregulated in stage I-IV, metastatic vs. nonmetastatic groups, in every molecular subtype, and in patients that were either on treatment or naïve to any kind of breast cancer treatment as compared to healthy controls as shown in Figure 11. But compared to KLF3, the expression of miR-124 was downregulated in all molecular subtypes. miR-124 expression was also downregulated in both treated and nontreated groups in comparison to the control. But its expression was upregulated in the pretreated group relative to the nontreated group.

(a)

(b)

(c)

(d)

Figure 10

qPCR analysis of KLF3 expression with clinical characteristics of breast cancer cases. Correlation of KLF3 expression with (a) tumor stage, (b) metastasis, (c) receptor status, and (d) treatment condition. There is upregulated expression of KLF3 in lower tumor grades (I and II), nonmetastasis, cancer subtype (triple negative), and naïve-treated groups. The data expressed as fold change represents the experiments performed in triplicate. Ordinary one-way ANOVA was used to establish significance ().

(a)

(b)

(c)

(d)

Figure 11

qPCR analysis of miR-124 expression with clinical characteristics of breast cancer cases. Correlation of miR-124 expression with (a) tumor grade, (b) metastasis, (c) receptor status, and (d) treatment condition. There is upregulated expression of miR-124 in low-grade tumor (I and II), nonmetastasis, cancer subtype (luminal A), and pretreated groups. The data expressed as fold change represents the experiments performed in triplicate. Ordinary one-way ANOVA was used to establish significance ().

The expression of PKCε was downregulated in all cancer stages, in molecular subtypes, in metastatic and nonmetastatic groups, and in patients that were either on treatment or naïve to any kind of breast cancer treatment as compared to healthy controls (see Figure 12).

(a)

(b)

(c)

(d)

Figure 12

qPCR analysis of PKCε expression with clinical characteristics of breast cancer cases. Correlation of PKCε expression with (a) tumor stage, (b) metastasis, (c) receptor status, and (d) treatment condition. There is upregulated expression of PKCε in high-grade tumor (III and IV), nonmetastasis, cancer subtype (luminal A), and naïve-treated groups. The data expressed as fold change represents the experiments performed in triplicate. Ordinary one-way ANOVA was used to establish significance ().

3.8. Risk Factors for Breast Cancer

Chance of breast cancer risk relative to age is shown in Figure 13. Mean age data shows that risk of breast cancer is higher in patients of . The data is significant with calculated by Student’s -test. Simple data analysis regarding other risk factors (see Table 2) has shown that 44% patients with age more than 50 years were on postmenopause and 56% patients with age less than 49 were on premenopause at the time of encountering breast cancer. Six percent of the cases were those who consume cigarette while 94% patients had never took cigarette in their entire life. Moreover, 56% breast cancer patients had history of cancer in their family and 88% patients had breastfed their children (see Figure 14).

4. Discussion

Current diagnostic analysis of breast cancer includes pathological assessments through clinical examination as well as imaging techniques. Numerous techniques are being tried for diagnostic imaging and screening. Among them, some techniques are not in operational use and the rest often times fail to detect breast carcinoma at the treatable stage. Increasing cancer incidences and inadequacy of efficient diagnostic strategy have given rise to a strong need to establish novel diagnostic and prognostic approaches for breast cancer [34]. Hence, the present study intends to find the expression pattern of KLF3, miR-124, TPD52, and PKCε in breast cancer patients which may serve as effective diagnostic biomarkers. Although the expression pattern of understudied genes have been analyzed individually and reported in several cancers following different signals, mechanisms, and factors, the coexpression of these genes in breast cancer was never reported. Further, results of the literature review and several other databases including the Protein Data Bank (RCSB PDB) revealed that the X-ray crystallography structure of KLF3 is yet to be determined. Hence, this study is conducted to establish 3D structure of KLF3 as well as the other features of KLF3 protein such as localization, conserved domains, and functional and phylogenetic analysis to build evolutionary relationship with other members of the family.

In silico results of our study have shown that KLF3 is located in the nucleus which was confirmed through UniProtKB (entry no. P57682). Evolutionary conservation was confirmed by the presence of conserved domains: three zinc fingers among members of KLFs, and it was also relatable to literature [18] and the UniProtKB database. Our study is also analogous to the Pearson review about the placement of KLF3 in group 2, and this study established the relation among family members in the form of a phylogenetic tree. The Ramachandran plot was used to assess the quality of 3D structure of KLF3, and the score of 95.06% depicted good quality of model. The proposed structure of KLF3 is 80.72% identical to the template model of KLF4. KLF4 anticancer activity in breast cancer is previously reported by several studies. It has a growth inhibitory role in ER⁺ breast cancer [35]. In vitro knockdown studies reported its role as a repressor of cell proliferation and angiogenesis [36, 37]. Considering the homology ratio of KLF3 with KLF4 and KLF4 tumor-suppressor role, it can be assumed that KLF3 might also function to inhibit carcinogenicity in breast cancer. KLF3 expression analysis revealed that its expression is downregulated with a fold change of 0.07. Present study expression analysis also suggests the tumor-suppressor role of KLF3 in breast cancer. Similar to our findings, significant reduction in expression of KLF3 is also reported in human metastatic sarcomas [19]. However, validation of KLF3 expression on larger cohort size and unraveling of its molecular role in breast cancer will further strengthen its significance as a biological marker.

Contrary to KLF3, in vitro results of the present study have shown the upregulated expression of TPD52 (2-fold) in breast cancer patients as compared to healthy controls. Our results were consistent with previous studies in which TPD52 expression was found to be significantly increased in ovarian and prostate cancer [38]. Moreover, the current study has depicted downregulation of miR-123 expression with a fold change of 0.06. Wang et al. [39] reported underexpression of miR-124 in pancreatic duct adenocarcinoma (PDAC) and showed its significant association with poor survival of patients. Their study outcomes further revealed that miR-124 targets the 3-untranslated region of Rac1 (supposed tumor promoter in PDAC) directly. Rac1 miR-124-mediated downregulation causes inactivation of the MKK4-JNK-c-Jun signaling pathway and prevents uncontrolled cell proliferation. miR-124 serves as a tumor suppressor, and its silencing promotes cell proliferation in pancreatic cancer.

Various studies reported upregulation of PKCε expression in numerous carcinomas such as prostate, brain, and lung cancer [40]. According to previous evidences on the expression of PKCε, there exists an overexpression of PKCε in breast cancer [41]. However, our study outcomes revealed a significant downregulation (0.01-fold) of PKCε in breast cancer patients which are contradictory to past studies.

All of the genes involved in our study are known to be part of major cancer signaling pathways. Studies revealed that RTKs and GPCR are chiefly involved in modulating diverse signaling cascades. The downstream effector molecules of both receptors PI3K and Kras ultimately result in the activation of Akt and its associated pathways [9]. Akt signaling is a key cellular pathway responsible for cancer development, growth, survival, stage progression, drug resistance, and invasiveness [42]. RTKs are activated in the signaling pathways, so its activation leads to the stimulation of TPD52, located on chromosome 8. When TPD52 increases its copy number on one chromosome, it results in the increment of the MAL2 gene copy number [43]. These genes ultimately induce the activation of the Akt gene. Akt further activates downstream genes and leads the cell to its survival. miR-124—a tumor suppressor miRNA—plays a role in cell apoptosis, but in unfavorable conditions, the deregulated expression of miR-124 results in the proliferation of cells by aberrant activation of Kras and Akt pathways through the negative regulation of PTPN12 [13]. miR-124 downregulation results in the activation of cell cycle kinase (CDK4) which was an effector of Kras signaling and was activated by PKCε as a result of GPCR. The activation of CDK4 signals the activation of the transcription factor E2F1 which in turn is responsible for cell proliferation [44]. In the same way, miR-124 also results in the increased expression of SLUG which is then bound to the promoter of E-cadherin and promotes cell invasion [45]. miR-124 activates the alpha catalytic subunit of PI3K, and the activated PI3K further promotes Akt activation [46]. Moreover, the deregulated expression of the transcriptional corepressor KLF3 is involved in Kras signaling [10]. Reduced expression of KLF3 in breast cancer upregulates the KLF8 gene, which is actually downstream of FAK. These sets of genes stimulate the expression of metalloproteinase 9 (MMP 9) [47] which further promotes angiogenesis and is the one source of nutrition for tumor cells (see Figure 15).

Figure 15

The involvement of the Akt and Kras pathway in breast cancer. Akt and Kras pathways were majorly involved in the progression of breast cancer. These two pathways were regulated by TPD52, miR-124, KLF3, PKCε, and several other downstream effectors. The deregulated expression of TPD52, KLF3, miR-124, and PKCε led to the overexpression of these two oncogenic pathways which in turn promoted cell survival and proliferation. Moreover, the partner genes of KLF3, TPD52, miR-124, and PKCε were extracted through STRING software [53].

These molecules TPD52, KLF3, miR-124, and PKCε have functional significance in the regulation of the Kras and Akt pathway. So, in the current study, the association between the expressions of TPD52, KLF3, miR-124, and PKCε and clinicopathological attributes of breast cancer patients is also inferred. We found that the expression of TPD52 is higher in advanced-stage and metastatic breast cancer while the expression of KLF3 and miR-124 is lower in advanced clinical stages of tumor in comparison to initial stages. This suggests that the lower expression of KLF3 and miR-124 and enhanced expression of TPD52 have association with tumor metastatic potential and stage progression. Previously, elevated TPD52 expression in prostate cancer was reported to have association with tumor progression and metastasis [48]. Similarly, lower KLF3 and miR-124 expression was indicated in aggressiveness and metastatic uterine and cervical and ovarian cancers [49, 50]. Independent expression of TPD52 was evaluated by Roslan and colleagues in hormone receptor-positive breast cancer [51]. Their study findings also reported upregulated expression of TPD52 which is in accord with the present study result. But none of these studies evaluated the combined expression of these molecules in breast cancer, and the current study provides insight into functional relationship of these molecules with tumor aggressiveness and stage progression.

Similarly, we also investigated the expression pattern of these genes in treated and treatment-naïve patients. We found an elevated expression of TPD52 in treatment-naïve patients who were either on radiation therapy, chemotherapy, or chemoradiation therapy. This can be presumed that radio-, chemo-, or combined therapy in breast cancer influences the expression of oncogene TPD52 which suggests that TPD52 can be a potential therapeutic target for the metastatic and luminal B group of breast cancers.

The expression of PKCε in the present study is found significantly downregulated in all tumor stages (I-IV), metastatic/nonmetastatic group of patients, molecular subtypes, and both naïve/treated groups. This data is in contrast to evidences provided by Pan and colleagues who reported elevated expression of PKCε in breast cancer [41]. So, our study has reported a new response of the PKCε gene in terms of its expression pattern.

The extent of genetics influence on the expression of oncogenes is boosted by a favorable environment. So, in this study, some of the risk factors that could possibly cause breast cancer were also analyzed. According to the analysis generated by our study, older age () is one of the factors that increase the likelihood of breast cancer. The facts and figures of 2017-2018 by the American Cancer Society (ACS) have also highlighted older age as a risk factor for breast cancer. According to previous studies, 74% females who were at their premenopausal age (35-49 years), before encountering breast cancer, were at definite risk of developing breast cancer [52].

The analysis of these genes provides novel insight into the coexpression pattern of these genes in different types and stages of breast cancer. Through extensive molecular studies, their role in disease progression can be understood. Further, the study outcomes also highlighted the potential application of these genes as biomarkers for cancer diagnosis and prognosis. In terms of therapeutic mediation in breast cancer, PI3K and Kras pathways are the very attractive targets because the components of these pathways are well suited for pharmacological intervention.

5. Conclusions

The RNA samples isolated from the blood of breast cancer patients were analyzed by qPCR for the expression of TPD52, KLF3, miR-124, and PKCε. According to our results, the dysregulated expression pattern of all these genes is significantly coupled with disease progression. The fold change observed in samples has spotted the upregulation of TPD52 and downregulation of KLF3, miR-124, and PKCε in breast cancer. Likewise, enhanced TPD52 expression and lower KLF3 expression in the metastatic group of patients in comparison to the control and nonmetastatic group can be considerably utilized for diagnosis purpose of breast cancer at the molecular level. As these genes might have application as a predictor of the stage, metastasis, receptor, and treatment status, they can be used as a new set of biomarkers for breast cancer diagnosis. Further, the in vitro knockdown studies and transcriptome analysis of these genes will provide deep insight into the functioning of these molecules in breast cancer.

Data Availability

The underlying raw data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This study is financially supported by the National Research Support Programme (NRPU 10067; HEC) and National University of Sciences and Technology, Islamabad, Pakistan.

Supplementary Materials

Expression of TPD52, KLF3, miR-124, and PKCε in different clinical attributes (such as cancer stage, type, and treatment status) of breast cancer is available in Supplementary Table 1. (Supplementary Materials)

References

R. Chetty and M. R. Kalan, “Malignant granular cell tumor of the breast,” Journal of Surgical Oncology, vol. 49, no. 2, pp. 135–137, 1992.
View at: Publisher Site | Google Scholar
N. Tung, N. U. Lin, J. Kidd et al., “Frequency of germline mutations in 25 cancer susceptibility genes in a sequential series of patients with breast cancer,” Journal of Clinical Oncology, vol. 34, no. 13, pp. 1460–1468, 2016.
View at: Publisher Site | Google Scholar
C. Hu, E. C. Polley, S. Yadav et al., “The contribution of germline predisposition gene mutations to clinical subtypes of invasive breast cancer from a clinical genetic testing cohort,” JNCI: Journal of the National Cancer Institute, vol. 112, no. 12, pp. 1231–1241, 2020.
View at: Publisher Site | Google Scholar
C. Wendt and S. Margolin, “Identifying breast cancer susceptibility genes–a review of the genetic background in familial breast cancer,” Acta Oncologica, vol. 58, no. 2, pp. 135–146, 2019.
View at: Publisher Site | Google Scholar
Z. Javed, K. Khan, M. Iqbal et al., “Long non-coding RNA regulation of TRAIL in breast cancer: a tangle of non-coding threads,” Oncology Letters, vol. 20, no. 4, 2020.
View at: Publisher Site | Google Scholar
B. Alday-Parejo, F. Richard, J. Wörthmüller et al., “MAGI1, a new potential tumor suppressor gene in estrogen receptor positive breast cancer,” Cancers, vol. 12, no. 1, p. 223, 2020.
View at: Publisher Site | Google Scholar
K. Sriram, C. Salmerón, S. Z. Wiley, and P. A. Insel, “GPCRs in pancreatic adenocarcinoma: contributors to tumour biology and novel therapeutic targets,” British Journal of Pharmacology, vol. 177, no. 11, pp. 2434–2455, 2020.
View at: Publisher Site | Google Scholar
E. Pupo, D. Avanzato, E. Middonti, F. Bussolino, and L. Lanzetti, “KRAS-driven metabolic rewiring reveals novel actionable targets in cancer,” Frontiers in Oncology, vol. 9, p. 848, 2019.
View at: Publisher Site | Google Scholar
F. G. Alanazi, B. Q. Alenazi, A. D. Al-Faim et al., “Ras signaling pathway, historical view,” Journal of Cell Science & Therapy, vol. 5, no. 4, p. 1, 2014.
View at: Publisher Site | Google Scholar
M. E. Fernandez-Zapico, G. A. Lomberk, S. Tsuji et al., “A functional family-wide screening of SP/KLF proteins identifies a subset of suppressors of KRAS-mediated cell growth,” Biochemical Journal, vol. 435, no. 2, pp. 529–537, 2011.
View at: Publisher Site | Google Scholar
R. Wang, J. Xu, J. Xu et al., “miR-326/Sp1/KLF3: a novel regulatory axis in lung cancer progression,” Cell Proliferation, vol. 52, no. 2, article e12551, 2019.
View at: Publisher Site | Google Scholar
Z. Zhao, H. Liu, J. Hou et al., “Tumor protein D52 (TPD52) inhibits growth and metastasis in renal cell carcinoma cells through the PI3K/Akt signaling pathway,” Oncology Research Featuring Preclinical and Clinical Cancer Therapeutics, vol. 25, no. 5, pp. 773–779, 2017.
View at: Publisher Site | Google Scholar
A. K. Mehta, K. Hua, W. Whipple et al., “Regulation of autophagy, NF-κB signaling, and cell viability by miR-124 in KRAS mutant mesenchymal-like NSCLC cells,” Science Signaling, vol. 10, no. 496, article eaam6291, 2017.
View at: Publisher Site | Google Scholar
R. Garg, M. Cooke, F. Benavides et al., “PKCε is required for KRAS-driven lung tumorigenesis,” Cancer Research, vol. 80, no. 23, pp. 5166–5173, 2020.
View at: Publisher Site | Google Scholar
L. Chauvin, C. Goupille, C. Blanc et al., “Long chain n-3 polyunsaturated fatty acids increase the efficacy of docetaxel in mammary cancer cells by downregulating Akt and PKCε/δ-induced ERK pathways,” Biochimica et Biophysica Acta (BBA) - Molecular and Cell Biology of Lipids, vol. 1861, no. 4, pp. 380–390, 2016.
View at: Publisher Site | Google Scholar
M. D. Ilsley, K. R. Gillinder, G. W. Magor et al., “Krüppel-like factors compete for promoters and enhancers to fine-tune transcription,” Nucleic Acids Research, vol. 45, no. 11, pp. 6572–6588, 2017.
View at: Publisher Site | Google Scholar
R. Pearson, J. Fleetwood, S. Eaton, M. Crossley, and S. Bao, “Krüppel-like transcription factors: a functional family,” The International Journal of Biochemistry & Cell Biology, vol. 40, no. 10, pp. 1996–2001, 2008.
View at: Publisher Site | Google Scholar
R. C. Pearson, A. P. Funnell, and M. Crossley, “The mammalian zinc finger transcription factor Krüppel-like factor 3 (KLF3/BKLF),” IUBMB Life, vol. 63, no. 2, pp. 86–93, 2011.
View at: Publisher Site | Google Scholar
M. Sachdeva, R. D. Dodd, Z. Huang et al., “Epigenetic silencing of Kruppel like factor-3 increases expression of pro-metastatic miR-182,” Cancer Letters, vol. 369, no. 1, pp. 202–211, 2015.
View at: Publisher Site | Google Scholar
L. Smirnova, A. Gräfe, A. Seiler, S. Schumacher, R. Nitsch, and F. G. Wulczyn, “Regulation of miRNA expression during neural cell specification,” European Journal of Neuroscience, vol. 21, no. 6, pp. 1469–1477, 2005.
View at: Publisher Site | Google Scholar
J. Wei, F. Wang, L.-Y. Kong et al., “miR-124 inhibits STAT3 signaling to enhance T cell–mediated immune clearance of glioma,” Cancer Research, vol. 73, no. 13, pp. 3913–3926, 2013.
View at: Publisher Site | Google Scholar
Y. Fang, J. Qiu, Z. B. Jiang, S. R. Xu, Z. H. Zhou, and R. L. He, “Increased serum levels of miR-214 in patients with PCa with bone metastasis may serve as a potential biomarker by targeting PTEN,” Oncology Letters, vol. 17, no. 1, pp. 398–405, 2019.
View at: Publisher Site | Google Scholar
Y. Chen, S. Frost, and J. A. Byrne, “Dropping in on the lipid droplet-tumor protein D52 (TPD52) as a new regulator and resident protein,” Adipocytes, vol. 5, no. 3, pp. 326–332, 2016.
View at: Publisher Site | Google Scholar
M. A. Gorin and Q. Pan, “Protein kinase Cε: an oncogene and emerging tumor biomarker,” Molecular Cancer, vol. 8, no. 1, p. 9, 2009.
View at: Publisher Site | Google Scholar
X.-N. Chen, J. A. Knauf, R. Gonsky et al., “From amplification to gene in thyroid cancer: a high-resolution mapped bacterial-artificial-chromosome resource for cancer chromosome aberrations guides gene discovery after comparative genome hybridization,” The American Journal of Human Genetics, vol. 63, no. 2, pp. 625–637, 1998.
View at: Publisher Site | Google Scholar
Y. Akita, “Protein kinase Cε: multiple roles in the function of, and signaling mediated by, the cytoskeleton,” The FEBS Journal, vol. 275, no. 16, pp. 3995–4004, 2008.
View at: Publisher Site | Google Scholar
B. Huang, S. J. Fu, W. Z. Fan et al., “PKCε inhibits isolation and stemness of side population cells via the suppression of ABCB1 transporter and PI3K/Akt, MAPK/ERK signaling in renal cell carcinoma cell line 769P,” Cancer Letters, vol. 376, no. 1, pp. 148–154, 2016.
View at: Publisher Site | Google Scholar
F. Sievers and D. G. Higgins, “Clustal Omega for making accurate alignments of many protein sequences,” Protein Science, vol. 27, no. 1, pp. 135–145, 2018.
View at: Publisher Site | Google Scholar
G. Pollastri and A. Mclysaght, “Porter: a new, accurate server for protein secondary structure prediction,” Bioinformatics, vol. 21, no. 8, pp. 1719-1720, 2005.
View at: Publisher Site | Google Scholar
C. Geourjon and G. Deleage, “SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments,” Bioinformatics, vol. 11, no. 6, pp. 681–684, 1995.
View at: Publisher Site | Google Scholar
J. Yang, R. Yan, A. Roy, D. Xu, J. Poisson, and Y. Zhang, “The I-TASSER Suite: protein structure and function prediction,” Nature Methods, vol. 12, no. 1, p. 7, 2015.
View at: Publisher Site | Google Scholar
S. Kumar, G. Stecher, and K. Tamura, “MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets,” Molecular Biology and Evolution, vol. 33, no. 7, pp. 1870–1874, 2016.
View at: Publisher Site | Google Scholar
S. Kumar, G. Stecher, M. Li, C. Knyaz, and K. Tamura, “MEGA X: molecular evolutionary genetics analysis across computing platforms,” Molecular Biology and Evolution, vol. 35, no. 6, pp. 1547–1549, 2018.
View at: Publisher Site | Google Scholar
S. H. Jafari, Z. Saadatpour, A. Salmaninejad et al., “Breast cancer diagnosis: imaging techniques and biochemical markers,” Journal of Cellular Physiology, vol. 233, no. 7, pp. 5200–5213, 2018.
View at: Publisher Site | Google Scholar
K. Akaogi, Y. Nakajima, I. Ito et al., “KLF4 suppresses estrogen-dependent breast cancer growth by inhibiting the transcriptional activity of ERα,” Oncogene, vol. 28, no. 32, pp. 2894–2902, 2009.
View at: Publisher Site | Google Scholar
A. Ray, M. Alalem, and B. K. Ray, “Loss of epigenetic Kruppel-like factor 4 histone deacetylase (KLF-4-HDAC)-mediated transcriptional suppression is crucial in increasing vascular endothelial growth factor (VEGF) expression in breast cancer,” Journal of Biological Chemistry, vol. 288, no. 38, pp. 27232–27242, 2013.
View at: Publisher Site | Google Scholar
J. Ferralli, R. Chiquet-Ehrismann, and M. Degen, “KLF4α stimulates breast cancer cell proliferation by acting as a KLF4 antagonist,” Oncotarget, vol. 7, no. 29, pp. 45608–45621, 2016.
View at: Publisher Site | Google Scholar
J. A. Byrne, S. Maleki, J. R. Hardy et al., “MAL2 and tumor protein D52 (TPD52) are frequently overexpressed in ovarian carcinoma, but differentially associated with histological subtype and patient outcome,” BMC Cancer, vol. 10, no. 1, p. 497, 2010.
View at: Publisher Site | Google Scholar
P. Wang, L. Chen, J. Zhang et al., “Methylation-mediated silencing of the miR-124 genes facilitates pancreatic cancer progression and metastasis by targeting Rac1,” Oncogene, vol. 33, no. 4, pp. 514–524, 2014.
View at: Publisher Site | Google Scholar
F. Zheng, Y.-J. Liao, M.-Y. Cai et al., “The putative tumour suppressor microRNA-124 modulates hepatocellular carcinoma cell aggressiveness by repressing ROCK2 and EZH2,” Gut, vol. 61, no. 2, pp. 278–289, 2011.
View at: Publisher Site | Google Scholar
Q. Pan, L. W. Bao, C. G. Kleer et al., “Protein kinase Cε is a predictive biomarker of aggressive breast cancer and a validated target for RNA interference anticancer therapy,” Cancer Research, vol. 65, no. 18, pp. 8366–8371, 2005.
View at: Publisher Site | Google Scholar
G. Hoxhaj and B. D. Manning, “The PI3K–AKT network at the interface of oncogenic signalling and cancer metabolism,” Nature Reviews Cancer, vol. 20, no. 2, pp. 74–88, 2020.
View at: Publisher Site | Google Scholar
H. Chen, G. Pimienta, Y. Gu et al., “Proteomic characterization of Her2/neu-overexpressing breast cancer cells,” Proteomics, vol. 10, no. 21, pp. 3800–3810, 2010.
View at: Publisher Site | Google Scholar
T. Feng, F. Shao, Q. Wu et al., “miR-124 downregulation leads to breast cancer progression via LncRNA-MALAT1 regulation and CDK4/E2F1 signal activation,” Oncotarget, vol. 7, no. 13, pp. 16205–16216, 2016.
View at: Publisher Site | Google Scholar
Y.-J. Liang, Q.-Y. Wang, C.-X. Zhou et al., “miR-124 targets Slug to regulate epithelial–mesenchymal transition and metastasis of breast cancer,” Carcinogenesis, vol. 34, no. 3, pp. 713–722, 2013.
View at: Publisher Site | Google Scholar
Q. Lang and C. Ling, “miR-124 suppresses cell proliferation in hepatocellular carcinoma by targeting PIK3CA,” Biochemical and Biophysical Research Communications, vol. 426, no. 2, pp. 247–252, 2012.
View at: Publisher Site | Google Scholar
R. Limame, K. Op de Beeck, F. Lardon, O. De Wever, and P. Pauwels, “Krüppel-like factors in cancer progression: three fingers on the steering wheel,” Oncotarget, vol. 5, no. 1, pp. 29–48, 2014.
View at: Publisher Site | Google Scholar
M. A. Rubin, S. Varambally, R. Beroukhim et al., “Overexpression, amplification, and androgen regulation of TPD52 in prostate cancer,” Cancer Research, vol. 64, no. 11, pp. 3814–3822, 2004.
View at: Publisher Site | Google Scholar
H. Lyng, R. S. Brøvig, D. H. Svendsrud et al., “Gene expressions and copy numbers associated with metastatic phenotypes of uterine cervical cancer,” BMC Genomics, vol. 7, no. 1, 2006.
View at: Publisher Site | Google Scholar
H. Zhang, Q. Wang, Q. Zhao, and W. Di, “miR-124 inhibits the migration and invasion of ovarian cancer cells by targeting SphK1,” Journal of Ovarian Research, vol. 6, no. 1, p. 84, 2013.
View at: Publisher Site | Google Scholar
N. Roslan, I. Bièche, R. K. Bright, R. Lidereau, Y. Chen, and J. A. Byrne, “TPD52 represents a survival factor in ERBB2-amplified breast cancer cells,” Molecular Carcinogenesis, vol. 53, no. 10, pp. 807–819, 2014.
View at: Publisher Site | Google Scholar
D. Apter, M. Reinilä, and R. Vihko, “Some endocrine characteristics of early menarche, a risk factor for breast cancer, are preserved into adulthood,” International Journal of Cancer, vol. 44, no. 5, pp. 783–787, 1989.
View at: Publisher Site | Google Scholar
G. C. K. W. Koh, P. Porras, B. Aranda, H. Hermjakob, and S. E. Orchard, “Analyzing protein–protein interaction networks,” Journal of Proteome Research, vol. 11, no. 4, pp. 2014–2031, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Khushbukhat Khan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

961

Downloads

1177

Citations