- Open Access
Characterization of the molecular changes associated with the overexpression of a novel epithelial cadherin splice variant mRNA in a breast cancer model using proteomics and bioinformatics approaches: identification of changes in cell metabolism and an increased expression of lactate dehydrogenase B
Cancer & Metabolism volume 7, Article number: 5 (2019)
Breast cancer (BC) is the most common female cancer and the leading cause of cancer death in women worldwide. Alterations in epithelial cadherin (E-cadherin) expression and functions are associated to BC, but the underlying molecular mechanisms have not been fully elucidated. We have previously reported a novel human E-cadherin splice variant (E-cadherin variant) mRNA. Stable transfectants in MCF-7 human BC cells (MCF7Ecadvar) depicted fibroblast-like cell morphology, E-cadherin wild-type downregulation, and other molecular changes characteristic of the epithelial-to-mesenchymal transition process, reduced cell-cell adhesion, and increased cell migration and invasion. In this study, a two-dimensional differential gel electrophoresis (2D-DIGE) combined with mass spectrometry (MS) protein identification and bioinformatics analyses were done to characterize biological processes and canonical pathways affected by E-cadherin variant expression.
By 2D-DIGE and MS analysis, 50 proteins were found differentially expressed (≥ Δ1.5) in MCF7Ecadvar compared to control cells. Validation of transcript expression was done in the ten most overexpressed and underexpressed proteins. Bioinformatics analyses revealed that 39 of the 50 proteins identified had been previously associated to BC. Moreover, metabolic processes were the most affected, and glycolysis the canonical pathway most altered. The lactate dehydrogenase B (LDHB) was the highest overexpressed protein, and transcript levels were higher in MCF7Ecadvar than in control cells. In agreement with these findings, MCF7Ecadvar conditioned media had lower glucose and higher lactate levels than control cells. MCF7Ecadvar cell treatment with 5 mM of the glycolytic inhibitor 2-deoxy-glucose led to decreased cell viability, and modulation of LDHB expression in MCF7Ecadvar cells with a specific small interfering RNA resulted in decreased cell proliferation. Finally, a positive association between expression levels of the E-cadherin variant and LDHB transcripts was demonstrated in 21 human breast tumor tissues, and breast tumor samples with higher Ki67 expression showed higher LDHB mRNA levels.
Results from this investigation contributed to further characterize molecular changes associated to the novel E-cadherin splice variant expression in BC cells. They also revealed an association between expression of the novel variant and changes related to BC progression and aggressiveness, in particular those associated to cell metabolism.
Breast cancer (BC) is the most common cancer in women, representing 25% of female cancers. Moreover, it is the fifth cause of cancer death and women death leading cause by cancer . Taking into account the high incidence and mortality associated to BC, it is still a great challenge to delve into the molecular mechanisms involved in the origin and progression of this heterogeneous disease, to better understand its pathogenesis and to identify novel biomarkers and new therapeutic targets.
Epithelial cadherin (E-cadherin) is a calcium-ion dependent cell-cell adhesion molecule that plays a key role in the maintenance of normal epithelia intercellular junctions and in a variety of processes fundamental to tissue homeostasis, such as differentiation, proliferation, migration, cell death, and survival . A partial or total loss of E-cadherin expression and/or function(s) has been reported in numerous tumors of epithelial origin, among them in breast tumors, leading to destabilization of cell-cell interactions and aberrant activation of cellular signaling pathways [3, 4]. A reduced expression of E-cadherin is a hallmark of epithelial-mesenchymal transition (EMT), a process that induces changes in gene expression and in cell behavior that have been associated with tumor progression and aggressiveness . Despite its relevance, the underlying molecular alterations leading to deregulation of E-cadherin expression and function(s) have not been completely characterized. Among the most frequent mechanisms responsible of these changes are a loss of heterozygosity, mutations, epigenetic silencing, expression of transcriptional repressors, and post-translational modifications [6, 7]. In recent years, a differential regulation of alternative splicing (AS) and its association with molecular and phenotypic changes has been reported during EMT [8, 9]. Under normal conditions, AS regulates gene expression and contributes to the expansion of proteomic diversity in a spatially and temporally controlled fashion. However, during malignant transformation, AS events have been reported in tumor suppressor genes or oncogenes, which result in association with mutations at specific splice sites and/or with changes in the expression of AS regulatory factors, among other causes, changes that contribute to tumor progression and the acquisition of therapeutic resistance [10, 11].
We have recently reported, for the first time, a novel transcript of human E-cadherin that differs from the wild-type E-cadherin (E-cadherin wt) mRNA by lacking the first 34 bp of Exon 14. This novel variant, called E-cadherin variant, would result from an AS event of the immature E-cadherin mRNA. Its expression levels would be regulated, at least in part, by the nonsense-mediated mRNA decay (NMD) surveillance mechanism. The E-cadherin variant mRNA was detected in cell lines derived from different human tumors types. Moreover, stable transfection of the E-cadherin variant in the MCF-7 human BC cell line (MCF7Ecadvar) induces a significant decrease in the E-cadherin wt (transcript and protein) expression, as well as changes that resemble the EMT process. MCF7Ecadvar cells depict a fibroblast-like morphology, and show changes in the expression levels of several EMT markers (i.e. cytokeratins, vimentin, and neural cadherin), reduced cell-cell interaction, increased cell migration and cell invasive properties when compared to the control cell line .
The present study aimed to thoroughly characterize molecular changes associated to the novel E-cadherin splice variant mRNA expression in MCF-7 cells, by performing a proteomic analysis using the two-dimensional differential gel electrophoresis (2D-DIGE) technology combined with protein identification by mass spectrometry (MS). In addition, this study aimed to identify other cellular processes and cell signaling pathways affected by the expression of the novel E-cadherin splice variant mRNA using bioinformatics tools, and to compare these changes with those previously associated with BC. Finally, studies on the lactate dehydrogenase B chain (LDHB) enzyme were done; this protein was identified in the proteomic analysis with the highest increased expression in MCF7Ecadvar compared to MCF7pcDNA3 control cells.
2D-DIGE and MS proteomic analysis of MCF7Ecadvar cells
To assess the global protein profile of MCF-7 cells stably transfected with the novel E-cadherin variant transcript, an expression analysis was carried out using the 2D-DIGE technique coupled to protein identification by MS. For this purpose, MCF7Ecadvar and MCF7pcDNA3 cell cultures (three per cell line) were prepared and processed by means of total protein precipitation and quantification, as detailed in “Materials and Methods” section. Downregulation of E-cadherin wt in MCF7Ecadvar cells was confirmed by Western immunoblotting of total protein extracts (Additional file 1: Figure S1).
Total protein extracts from MCF7pcDNA3 and MCF7Ecadvar cells were then labeled with fluorescent dyes, mixed, and separated by isoelectric focusing (first-dimension) followed by denaturing polyacrylamide gel electrophoresis (second-dimension). After gel scanning, protein alignment, and analysis, 64 protein spots were identified as differentially expressed between both cell lines (1.5-fold change cutoff and an analysis of variance (ANOVA) ≤ 0.005). Of these spots, 60 were recovered and their spectra were obtained by MS. The identity of each protein spot was determined using the Mascot Database Search Program (Matrix Science, London, UK). As result of these studies, 50 proteins were identified: ten proteins were found overexpressed and 40 underexpressed in MCF7Ecadvar cells relative to MCF7pcDNA3 cells. The name and symbol of these proteins, their UniProt number, the symbol of the gene encoding each protein, a brief protein functional description, and the fold change value associated to each protein are included in Table 1.
Biological characterization of the proteomic analysis results
To analyze biological characteristics of the 50 differentially expressed proteins found in MCF7Ecadvar cells, a set of bioinformatics tools were applied. Firstly, proteins were classified using the Protein ANalysis THrough Evolutionary Relationships (PANTHER) tool, by means of their molecular function (Fig. 1a) and the biological processes (Fig. 1b) in which they were involved. As result, catalytic activity was the most represented molecular function (56.0%; 27/50 proteins). Other categories listed were binding, structural molecule, antioxidant activity, transporter, and translation regulator. The energy metabolism was identified as the most affected biological process (34.5%), followed by cellular process (32.2%).
Additional evaluations with the PANTHER tool included cellular component and protein class analysis; results are listed in Additional file 2: Figure S2. With regard to the cellular component classification, the cell part category represented a 52.5% (23/50). Among the protein class categories, 21 (SODC, PRDX2, ALDR, LDHB, 6PGL, G6PD, ECI1, MAOX, PRDX6, CATD, EF2, PSB3, PSA6, F16P1, ACPH, PMM2, TPIS, HPRT, G6PI, GDIB, and UBA1) of the 50 differentially expressed proteins in MCF7Ecadvar cells were included in the six enzymatic categories listed in the protein class classification.
The ingenuity pathway analysis (IPA) bioinformatics tool was run to evaluate the ten canonical pathways most affected in the MCF7Ecadvar cellular model; results were in agreement with the PANTHER analysis (Table 2).
In particular, the most altered pathways listed in IPA were glycolysis and gluconeogenesis, two metabolic processes composed by a set of consecutive enzymatic reactions. The pentose phosphate metabolic pathway, an alternative branch of glycolysis to produce sugars that make up DNA and RNA, was also listed in the seventh position. In third position, the protein ubiquitination pathway was found, involving 8/50 proteins identified in the MCF7Ecadvar cell model. This pathway regulates degradation of cellular proteins by the ubiquitin-proteasome system, controlling protein’s half-life and expression levels as well as key cellular processes . The p70S6 kinase (p70S6K) signaling pathway (fifth position) is an important regulator of the cell cycle, which activation requires phosphatidylinositol-4,5-bisphosphate 3-kinase (PI3K)/RAC-alpha serine/threonine-protein kinase (AKT)-dependent signal(s), a pathway found in sixth position . Furthermore, key cellular functions are regulated by the other four canonical pathways identified by the IPA tool, such as cell growth (pathways number 4 and 10), ion transport (pathways number 4 and 8), inflammation, immunity, and metabolism (pathways number 4, 9, and 10) [15,16,17,18,19].
Using IPA, four networks were generated with the 50 proteins differentially expressed in the MCF7Ecadvar cell line. The score associated to each network, the number of molecules among the 50 identified by 2D-DIGE and MS that were part of them (in bold in the second column), and the main diseases and functions associated to each of them are detailed in Table 3.
Among the diseases associated to one of the two networks described with the highest score (network number 2) was the term “Cancer.” On the other hand, networks 3 and 4 were related with cellular processes usually altered during tumorigenesis. A graphical representation of each network was also generated with IPA and shown in Fig. 2, in which the type of connection between molecules and their function is indicated.
Association between proteomic analysis results and BC disease
In order to determine whether the proteins identified as differentially expressed in MCF7Ecadvar cells have already been reported to be associated with BC, the DisGeNET bioinformatics tool was used. A comprehensive study was conducted, in which all databases available in DisGeNET (curated, predicted, and literature) were considered to search for all the genes associated with terms that include the words “Breast” or “Mammary”: “Carcinoma,” “Cancer,” “Neoplasms,” or “Tumorigenesis”. In total, 36 terms were identified (March 24, 2018) and, from the gene lists associated with each of these terms, a list of 5261 genes (Additional file 3: Figure S3) was generated using the Venny 2.1.0 interactive tool. A total of 39 from the 50 identified genes had already been reported associated with these terms. Using the DisGeNET tool, it was also found that for ten of the 11 genes not previously related to the 36 analyzed terms, an altered expression had been reported in other carcinomas (Table 4).
Finally, the PSMB3 gene was not found in the DisGeNET database, but a manual search performed in PubMed led to the identification of one scientific paper that related PSMB3 with BC , and two reports with colorectal cancer and multiple myeloma [21, 22].
Validation studies of 2-DIGE and MS proteomic analysis results
Validation of the results obtained by the 2D-DIGE and MS technologies was done at the mRNA level, by quantitative real-time PCR. Specific primers were designed for the genes encoding the top 10 up- and downregulated proteins in MCF7Ecadvar cells relative to MCF7pcDNA3 (Additional file 4: Figure S4). Among the ten upregulated proteins, the same trend to an increased expression found at the protein level was determined for LDHB, AKR1B1, VIM, PDIA3, and LMNA transcripts in MCF7Ecadvar cells compared with MCF7pcDNA3 (Fig. 3a). On the other hand, no significant differences were observed in the expression levels of the EEF2, HSP90AB1, YWHAQ, and PDIA6 mRNAs between both cell lines (Fig. 3b). For the UBA1 protein-coding transcript, a significant decreased expression (p < 0.05) was registered in MCF7Ecadvar cells relative to the control cell line, contrasting with the increased protein level found by proteomics (Fig. 3c). Among the top 10 downregulated proteins, a decreased mRNA expression was also found in CRABP2, YWHAZ, G6PD, HSP90AA1, CTSD, and ANXA4 (Fig. 4a), while TPI1, PRDX2, and S100A11 transcript levels did not show significant differences between MCF7Ecadvar and MCF7pcDNA3 cells (Fig. 4b). Finally, an increased expression of the GGCT mRNA was found in MCF7Ecadvar cells compared to the control, results that contrasted with the decreased protein expression levels registered in the proteomic analysis (Fig. 4c).
LDHB expression analysis
Among the 50 proteins identified as differentially expressed in MCF7Ecadvar cells, the LDHB enzyme showed the highest level of increase (Fold = 9.34; Table 1). LDHB is a cytoplasmic protein that catalyzes the interconversion of pyruvate and lactate, with the concomitant interconversion of nicotinamide adenine dinucleotide (reduced form; NADH) and nicotinamide adenine dinucleotide (oxidized form; NAD+). Previous reports have associated changes in LDHB expression with tumor progression. However, while decreased LDHB levels were reported in prostate, bladder, and hepatocellular carcinomas [23,24,25], an increased LDHB expression was described in lung and pancreatic adenocarcinomas, non-small-cell lung cancer and osteosarcoma [26,27,28,29]. In BC, its expression has been reported in different cell lines and in human tumors [30, 31].
Immunodetection of the LDHB protein by Western immunoblotting revealed, under the conditions evaluated, no detectable levels of the enzyme in total protein extracts of MCF7pcDNA3 cells, whereas a 35 kDa signal was observed for the MCF7Ecadvar cell line (Fig. 5a). By immunofluorescence microscopy, a specific signal localized at the cytoplasm of MCF7pcDNA3 and MCF7Ecadvar cells was detected, which was 3.45 times higher in MCF7Ecadvar than in control cells (Fig. 5b).
Different types of cancers, including BC, show high rates of glucose consumption associated with lactate production, even in the presence of normal levels of oxygen, a phenomenon known as aerobic glycolysis or Warburg effect [32, 33]. When glucose concentration (mg/dL) was measured in the conditioned media of MCF7Ecadvar and MCF7pcDNA3 cell lines, a significant decrease (p < 0.05) was observed in the former relative to the control (Fig. 5c). The assessment of the lactate concentration revealed higher levels (p < 0.05) of this metabolite in the conditioned medium of MCF7Ecadvar cells than in the control cell medium (Fig. 5d). In agreement with these findings, a significant increase (p < 0.05) in the expression of the monocarboxylate transporters 1 (MCT1) and 4 (MCT4) transcripts was determined in MCF7Ecadvar cells compared with MCF7pcDNA3 (Additional file 5: Figure S5).
In order to evaluate whether the metabolic differences observed between MCF7Ecadvar and MCF7pcDNA3 cells affected cellular viability, cultures of both cell lines were incubated with the glycolytic inhibitor 2-deoxy-d-glucose (2DG), a modified-glucose molecule (substitution of hydroxyl group for hydrogen at the second carbon atom) that blocks glycolysis by inhibiting hexokinase activity, impairs cancer cell growth, and induces cell death through apoptosis . Incubation of MCF7Ecadvar cells with 5 mM 2DG led to a decrease (p < 0.001) in cell viability but not in control cells. Contrasting, the addition of 0.5 mM 2DG did not induce significant changes in MCF7Ecadvar cell viability compared with the control, suggesting a concentration-dependent effect of the inhibitor (Fig. 5e). Finally, the association between LDHB expression and MCF7Ecadvar cell proliferation was evaluated by transfecting cells with a LDHB specific small interfering RNA (siRNA). As shown in Fig. 5f, g, a decreased expression of the transcript and protein of LDHB was found in MCF7Ecadvar cells treated with 100 pmoL/mL of LDHB siRNA relative to the MCF7Ecadvar cells transfected with a scramble siRNA (siScr) used as negative control after 72 h of treatment. This negative modulation of LDHB expression resulted in a significant reduction (p < 0.05) of MCF7Ecadvar cell viability, when compared with results in the MCF7Ecadvar siScr treated cells (Fig. 5h).
E-cadherin variant and LDHB mRNA expression analysis in tumor tissue samples from BC patients
To translate the results obtained with the MCF7Ecadvar cellular model to the clinic, 21 breast tumor tissues from women who underwent biopsy and/or surgery for BC management before receiving radiotherapy or chemotherapy were analyzed. The histological and molecular characteristics of the samples are detailed in Additional file 6: Table S1.
Tissues were processed for total RNA extraction and the E-cadherin variant and LDHB transcripts expression levels were evaluated by quantitative real-time PCR. As a result, detectable amounts of both mRNAs were registered in the 21 breast tumors analyzed (Fig. 6a and Additional file 7: Table S2), although expression levels of each transcript differed among tumors.
Samples were numbered and organized in ascending order according to the 2-ΔCt values obtained for the E-cadherin variant transcript, which ranged from 0.000078 to 0.008669. For the LDHB mRNA, 0.011010 was the lowest 2-ΔCt value registered and 0.384302 the highest. Using the median value of E-cadherin variant and LDHB transcripts levels (0.000353 and 0.1073, respectively) as cutoff values to group samples in “low expression” or “high expression” for each molecule, a 2 × 2 contingency table was made and the Fisher’s exact test was applied. This analysis revealed a significant association (odds ratio 0.06, confidence interval 0.01–0.40, p < 0.0089) between the expression levels of the E-cadherin splice variant and LDHB transcripts. As shown in Fig. 6b, 80% of the breast tumor samples with low E-cadherin variant mRNA expression showed low levels of LDHB mRNA, while 81.82% of the samples with high E-cadherin variant expression levels presented high LDHB mRNA expression.
Finally, taking into account that an increased expression of LDHB has been associated with a high cellular viability and proliferation in previous reports [28, 29], and considering that this behavior was also observed in the MCF7Ecadvar cellular model (Fig. 5h), LDHB mRNA levels were compared in the 21 breast tumor tissues grouped according to the Ki-67 labeling index (samples with Ki-67 < 20% and samples with Ki-67 ≥ 20%). The immunohistochemical analysis of Ki-67 is regularly utilized to evaluate tumor cell proliferation, and is considered a prognostic and predictive marker in BC . As result of this analysis, breast tumor tissues from the Ki-67 ≥ 20% group showed higher LDHB transcript expression levels than those from the Ki-67 < 20% group.
In the mammary gland, E-cadherin is one of the main cellular adhesion molecules required for morphogenesis and maintenance of the epithelial structure as well as for epithelial cell survival . BC is a complex disease that involves cell dedifferentiation and proliferation of the mammary epithelium, leading to tissue disorganization. Deregulation of E-cadherin expression has been associated with BC progression and aggressiveness, and the underlying mechanisms involved in E-cadherin expression changes have been extensively studied [37, 38]. Dr. Vazquez-Levin’s group identified a novel human E-cadherin splice variant transcript that differs from the E-cadherin wt mRNA in 34 bp deleted in exon 14. Stable transfection of the variant transcript in MCF-7 cells showed a detrimental effect over the E-cadherin wt expression, and additional molecular changes characteristic of the EMT process and of an aggressive cellular behavior . Based on this background information, a global protein analysis was carried out using the 2D-DIGE and MS technologies to unleash the molecules and mechanisms of the cellular alterations observed associated with the overexpression of the E-cadherin variant in MCF-7 cells.
In the 2D-DIGE technique, protein mixtures of test and control samples are obtained, labeled with different fluorescent dyes, mixed, and subjected to 2D-electrophoresis, by which protein components are resolved based on their isoelectric point and apparent molecular weight. The inclusion of an internal standard containing all the proteins of the samples under study reduces gel-to-gel variations, allows for individual experiments normalization, and confers a high sensitivity and reproducibility for differential quantitative analysis of protein expression between experimental conditions . Despite having some limitations, such as under-representation of hydrophobic, membrane, and low abundance proteins, and the difficulty to resolve proteins with extreme isoelectric points or with very high or low molecular weight, the comparative 2D-DIGE method has been applied in biological and biomedical research areas to discover novel biomarkers, improve diagnostics and prognostics, and for therapy-monitoring. The potential of this approach lies in an accurate detection and quantification of protein isoforms, giving information that can be lost when using standard immunochemical assays and genomic analysis for protein detection [39, 40]. This electrophoresis procedure is followed by protein spot digestion and MS analysis, allowing peptide identification with high specificity and sensitivity. The 2D-DIGE technique coupled to MS protein identification was previously applied by members of the team , as well as by others, for the characterization of the global proteome of diverse experimental models. In particular, this approach was used for the assessing of EMT protein profiles [42,43,44,45]. In this regard, Vergara and collaborators  carried out a stable knockdown of E-cadherin wt expression in MCF-7 cells by transfecting with an shRNA, and changes in the proteomic signature were evaluated using a 2D-gel electrophoresis and MS approach, also observing changes in the expression of proteins involved in carbohydrate metabolism, such as PGAM1, also identified in this study, but with the opposite trend.
The proteomic study done on MCF7Ecadvar and MCF7pcDNA3 cells led to the identification of 50 proteins depicting a significant differential expression level between cell lines, which are involved in several signaling pathways and participate in the regulation of a variety of important cellular processes. In particular, 29 of the 50 proteins identified were associated with cell metabolism, and the glycolysis and gluconeogenesis metabolic pathways were listed as the most altered in the MCF7Ecadvar cellular model. Other processes categorized as response to stimulus, localization and organization of cellular components, biological regulation and development, were also listed. Among the 50 proteins, there were found enzymes (LDHB, ALDR, APT, TERA, HPRT, UBA1, G6PD, ACPH, NDKB, 6PGL, F16P1, G6PI, PGAM1, ECI1, ALDOA, MAOX, GGCT, PMM2, GPDM, TPIS, SODC, PRDX2, and PRDX6), proteases (CATD, PSB3, and PSA6), binding (ANXA4 and EMAL2) and transport (RABP2 and CLIC1) molecules, chaperones (HS90A, HS90B, HS71A, HS71B, GRP75, PDIA3, PDIA6, ERP29), structural molecules (ACTB, ACTA, VIME and LMNA), adapter/regulatory molecules (1433Z, 1433 T, GRB2, GDIB and EF2), splicing (HNRH1 and HNRPF), and differentiation (S10AB) factors. Transcript expression analysis showed the same trend of change for 11 (LDHB, AKRIB1, VIM, PDIA3, LMNA, CRABP2, YWHAZ, G6PD, HSP90AA1, CTSD, and ANXA4) of the 20 molecules evaluated in MCF7Ecadvar compared with MCF7pcDNA3 cells. Interestingly, the specific knockdown of E-cadherin wt was also associated with expression changes in a significant number of metabolic proteins and the catalytic activity was one of the most altered molecular function , suggesting that some of the molecular changes triggered by E-cadherin variant mRNA overexpression in MCF-7 could be associated with its role as a negative E-cadherin wt modulator. Part of the expression changes recorded for the molecules differentially expressed in the MCF7Ecadvar cell line were reported in samples from patients with different cancer types, in particular with BC. In some cases, these changes have also been associated with an aggressive tumor phenotype and a worse patient prognosis.
In the last decade, it was demonstrated that tumor cells exhibit an altered expression of different metabolic enzymes and metabolite transporters to increase the generation of energy necessary for a high proliferation and migration, as well as for survival and invasion . Glycolysis is a critical catabolic process through which a glucose molecule is cleaved to produce two molecules of pyruvate, two molecules of adenosine triphosphate (ATP), and two of NADH. In normal cells, the fate of pyruvate depends, to a large extent, on oxygen availability. In its presence, pyruvate is oxidized to acetyl coenzyme A (acetil-CoA), which is then completely oxidized to carbon dioxide and water in the tricarboxylic acid cycle, generating 3 NADH and 1 flavin adenine dinucleotide hydroquinone form (FADH2). As a result, the reducing power produced through these processes is used to yield large amounts of energy in the form of ATP through the mitochondrial oxidative phosphorylation system. In contrast, in the absence of oxygen, glycolysis is decoupled from oxidative phosphorylation and pyruvate is transformed into lactate through the anaerobic glycolysis metabolic pathway, becoming the primary source of ATP production.
Almost 100 years ago, Warburg and collaborators reported that aerobic glycolysis and lactate generation were increased in tumor cells, allowing a faster production of energy than through the oxidative phosphorylation carried out by normal cells . Among the ten proteins upregulated in the MCF7Ecadvar cellular model, LDHB showed the highest positive fold change. Since LDHB catalyzes interconversion of pyruvate (the glycolysis final product) and lactate, quantification of glucose and lactate levels in MCF7Ecadvar and MCF7pcDNA3 conditioned media was done. These studies revealed an increased glucose consumption and lactate production by MCF7Ecadvar cells. Moreover, the decrease in the expression of enzymes involved in gluconeogenesis (FBP1, GPI, PGAM1, ALDOA, ME1, and TPI1) and in the pentose phosphate pathway (G6PD and PGLS) registered in MCF7Ecadvar cells with respect to control would contribute to the increased availability of glucose-6-phosphate for the aerobic glycolysis over the pathways responsible of the synthesis of glycogen, fatty acids, and nucleic acids. An association between E-cadherin variant and LDHB levels was also demonstrated through the analysis of both transcripts in 21 human breast tumor samples.
To achieve the reprogramming of metabolism, tumor cells depend on the action of specific transcription factors, such as hypoxia-inducible factor 1-alpha (HIF1A) and MYC [49,50,51]. The activation of the PI3K, AKT, mTOR, Ras, Raf, Src, and EGFR oncogenes [52, 53], as well as the loss of tumor suppressor genes like PTEN and TP53 [54, 55], were also demonstrated to exacerbate tumor cell growth and metabolism. Some of these molecules were involved in the list of the ten canonical pathways identified by the IPA tool as the most altered in MCF7Ecadvar compared with MCF7pcDNA3 cells (mTOR, PI3K, and AKT) or were included as connectors into the IPA networks to link some of the proteomic molecules (EGFR, HIF1A, AKT1, MYC, and TP53). Particularly, it was reported that LDHB is targeted by the mTOR complexes and plays a fundamental role in the tumorigenesis triggered by the PI3K/AKT signaling pathway , all major mediators of cell survival. These findings could be related to the diminished cellular viability found in MCF7Ecadvar cells after the specific knockdown of LDHB expression. Moreover, the analysis of LDHB mRNA levels in 21 human breast tumor tissues showed that samples with Ki-67 labeling index ≥ 20% had higher levels of the LDHB transcript than those with a Ki-67 < 20%, reinforcing that the overexpression of the LDHB enzyme is necessary for a faster tumor growth [30, 31] and a more aggressive tumoral behavior [57, 58].
Taken together, results obtained from the 2D-DIGE and MS proteomic analysis done in MCF7Ecadvar and MCF7pcDNA3 cells revealed the complexity of the molecular effects associated with the overexpression of the novel E-cadherin splice variant mRNA in MCF-7 cells. While a large proportion (39 of 50) of the proteins identified as differentially expressed between both cell lines had been previously reported in the literature associated with BC, for some of them the available information is scarce and sometimes contradictory. Thus, the stable transfection of the novel E-cadherin variant transcript in MCF-7 cells may represent a useful experimental model to deepen in the knowledge of the cellular alterations associated with the deregulation of E-cadherin wt in this complex and heterogeneous disease. Specifically, this is the first report that describes the involvement of a novel E-cadherin splice variant in triggering molecular changes related to BC progression and aggressiveness and its relationship with LDHB expression.
Chemicals were of analytical or tissue culture grade and purchased from Sigma-Aldrich (Sigma; St. Louis, MO, USA). Molecular biology and electrophoresis reagents were purchased from Invitrogen (Thermo Fisher Scientific; Waltham, MA, USA), Qiagen (Hilden, Germany), and Bio-Rad (Hercules, CA, USA), unless specified.
The following antibodies were used: anti E-cadherin (610,181, mouse, monoclonal; Becton Dickinson Biosciences [BD], San Diego, CA, USA), anti LDHB (431.1, mouse, monoclonal; Santa Cruz Biotechnology [SCBT], Santa Cruz, CA, USA), and anti β-tubulin (clone D66, mouse, monoclonal; Sigma). For immunocytochemistry, a Cy3-labeled anti-mouse (Sigma) IgG was used as secondary antibody. An anti-mouse (Vector Laboratories Inc.; Burlingame, CA, USA) IgG coupled to horseradish peroxidase was employed as secondary antibody in Western immunoblotting assays.
In silico analysis
The following bioinformatics tools were used in this study:
Comprehensive software system for classifying proteins (and their genes) depending on the family or subfamily they belong, the molecular function they have, the biological process they are involved, or the way in which they participate on it. The classifications result from human data healing and sophisticated bioinformatics algorithms [59, 60].
Discovery platform that integrates information on gene-disease associations by text-mining approaches from expert-curated databases and scientific literature. Data is organized according to their type and level of curation in curated (gene-disease associations from UNIPROT and CTD human databases), predicted (CTD mouse and rat data, RGD and MGD), literature (GAD, LHGDN, and BeFree), and all (curated, predicted, and literature). The gene-disease associations can be ranked according to a 0 to 1 DisGeNET score and are annotated with the DisGeNET gene-disease association type ontology. The score takes into account the number and type of sources (i.e., level of curation and organisms), as well as the number of publications supporting the association. It can be accessed through the web interface  or using a Cytoscape plugin  (see below).
Online tool to compare lists of genes or proteins by performing Venn diagrams. It also allows visualization of specific elements of each list, as well as common elements between them .
Web-based software application for the analysis, integration, and interpretation of data derived from omics experiments. It also transforms a list of genes or gene products into a set of networks based on the ingenuity knowledge base, which has been abstracted into a large network called global molecular network composed of thousands of genes that interact with each other. The relationship between two genes is called connection and two genes are considered to be connected if there is a path between them in the network, i.e., a set of genes connecting one to the other. If the relationship between genes of a list of interest is not direct, the IPA program relates them through other neighboring genes present in the global molecular network. In addition, IPA assigns a score value to each network and the higher the score the higher the degree of confidence for the connections identified in the network .
All procedures reported in this study were done following protocols that complied with biosafety and ethical guidelines accepted worldwide.
Twenty-one tumor tissue samples from patients diagnosed with BC of different subtypes and stages of the disease were obtained from the Hospital Italiano of Buenos Aires and the Instituto de Oncología Ángel H. Roffo, Buenos Aires, Argentina. Sample collection protocols were approved by the institutional Ethics Committees and patient informed consents were obtained in all cases. The study was conducted in accordance with the principles of the Declaration of Helsinki. Samples were from women older than 18 years old, who underwent biopsy and/or surgery because of BC diagnosis. The breast tumor tissues of patients that received radiotherapy or chemotherapy before surgery were excluded from the study.
Cell lines were grown in appropriate culture conditions as previously described . Briefly, MCF-7 human BC cells, as well as MCF7Ecadvar and MCF7pcDNA3 stably transfected cells, were cultured in sterile culture dishes (Nunc, Thermo Fisher Scientific) at 37 °C in a humid atmosphere containing 5% CO2 in air. Cells were grown in the Dulbecco’s modified Eagle’s medium Ham’s Nutrient Mixture F-12 (DMEM F-12; Gibco, Thermo Fisher Scientific) cell culture medium, supplemented with 10% fetal calf serum (FBS), 100 U/mL penicillin, 100 mg/mL streptomycin, and 2 mM glutamine. Upon reaching 80% confluence, cells were sub-cultured using a 0.25% sterile solution of trypsin and 0.025% ethylenediaminetetraacetic acid (EDTA) in phosphate-buffered saline (PBS).
Proteomic analysis was done employing 2D-DIGE technology coupled to protein identification by MS at the Proteomics Service of the Vall d’Hebron Institute of Oncology, Vall d’Hebron University Hospital, Barcelona, Spain, basically as previously reported .
For sample preparation, triplicates of MCF7Ecadvar and MCF7pcDNA3 cell cultures were grown in Petri dishes (Nunc, Thermo Fisher Scientific) as described before. When cells reached 60–70% confluence, the growth medium was replaced by medium without FBS supplementation. After 24 h, cultures were collected with 5 mM EDTA, centrifuged 5 min at 200×g and stored at − 70 °C until use.
For the proteomic analysis, cellular pellets were lysed with 400 μL of DIGE lysis buffer (pH 8.5; 30 mM Tris, 7 M urea, 2 M thiourea, 4% CHAPS) to achieve protein extraction. After sonication, protein extracts were purified using the 2D-CleanUp kit (GE Healthcare; Buckinghamshire, UK) and re-dissolved in DIGE lysis buffer. Protein concentration was determined with the Bio-Rad RC DC Protein Assay (Bio-Rad). Proteins were then labeled using the Ettan DIGE method (GE Healthcare). A pool consisting of equal amounts of each of the samples analyzed in this study was prepared to be used as an internal standard for quantitative comparisons. Triplicates from MCF7Ecadvar and MCF7pcDNA3 cell cultures were alternatively labeled with Cy3 and Cy5 cyanine dyes to avoid possible bias by labeling efficiency, and the internal standard was labeled with the Cy2 dye. The recommended marker/protein ratio of 400 picomoles of marker per 50 μg of total proteins was used.
MCF7Ecadvar and MCF7pcDNA3 total protein samples were separated in a first-dimension by the IEF technique using the IPGphor system (GE Healthcare). The IPGphor strips have a length of 24 cm and cover a range of 3 to 10 pH units. Isoelectrical focusing was carried out overnight or until reaching 8000 Vh. Total proteins separation in a second-dimension was carried out by laying the strips on 12.5% isocratic Laemmli gels (24 × 20 cm) on an Ettan DALTsix Electrophoresis System (GE Healthcare). Gels were run at 20 °C at constant power (2.5 W) per gel for 30 min, followed by 17 W per gel until the bromophenol blue tracking front reached the end of the gel. When the run was finished, gels were scanned with the Typhoon FLA 7000 Scanner (GE Healthcare; Cy2: 488/520 nm, Cy3: 532/580 nm and Cy5: 633/670 nm excitation/emission wavelengths). Image analysis and statistical quantification of relative protein abundances were performed using Progenesis Samespots (Nonlinear Dynamics, UK) and a DeCyder 6.0 version software (GE Healthcare). Protein spots that met a 1.5-fold change cutoff and an ANOVA ≤ 0.005 were considered. Once detected, the spots differentially expressed between MCF7Ecadvar and MCF7pcDNA3 samples were excised from the gels using the Ettan SpotPicker System (GE Healthcare). In-gel trypsin digestion was performed as previously described  using autolysis stabilized trypsin (Promega, Madison, WI, USA) and tryptic digests were purified with ZipTip Pipette Tips (EMD Millipore, Cork, Ireland). Spectra of the protein spots selected were obtained by the matrix-assisted laser desorption ionization-time of flight (MALDI-TOF), MALDI/TOF-TOF, and liquid chromatography-electrospray ionization-MS in tandem (LC-ESI-MS/MS) techniques. Protein identity determination in each spot recovered was done through the Mascot Database Search Program.
Transfection with LDHB siRNA
To knockdown endogenous LDHB gene expression, 250,000 MCF7Ecadvar cells per well were seeded in 6-well plates (Nunc, Thermo Fisher Scientific). After 24 h, cells were transfected with 100 pmoL/mL of a specific LDHB or a scramble (SCBT) siRNA, using Lipofectamine 2000 (Thermo Fisher Scientific) and following manufacturer’s instructions. Seventy-two hours after transfection, cells were recovered with EDTA 5 mM to perform LDHB mRNA and protein expression analysis. In addition, assessment of cellular viability was done with the CellTiter 96 AQueous Non-Radioactive Cell Proliferation Assay (Promega; Madison, WI, USA). Briefly, three replicates of 5000 cells were done from both LDHB and control siRNA transfections, and cultured for 5 days at 37 °C in 5%CO2 in air. Cells were harvested and each replicate was seeded in 5-wells of a 96-well plate. Then, 20 μL of the reagent solution of 2 mL MTS [3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium] and 0.1 mL PMS (phenazine methosulfate) were added to each well and incubated for 1 h, after which absorbance was read at 490 nm on a microtiter spectrophotometer plate reader (Thermo Fisher Scientific). Controls without cells were included for background estimation.
Total RNA extraction and real-time PCR procedures
Total RNA from MCF-7, MCF7pcDNA3, and MCF7Ecadvar cell cultures, as well as from breast tumor tissue samples, was isolated with the Trizol reagent (Thermo Fisher Scientific) following the manufacturer’s instructions and quantified (Nanodrop, Thermo Fisher Scientific). Synthesis of complementary DNA (cDNA) was performed with oligo-dT and SuperScript™ III Reverse Transcriptase (Thermo Fisher Scientific). Negative controls omitting the RNA or the reverse transcriptase enzyme were included in all cases. Standard PCR amplification protocols were carried out using the TaqI DNA polymerase (Qiagen). Reverse transcription procedures were verified by amplification of a glyceraldehyde 3-phosphate dehydrogenase (GAPDH) cDNA fragment.
Quantitative assessment of gene expression was done by real-time PCR with the Bio-Rad Real Time PCR unit using the SYBR Select Master Mix (Thermo Fisher Scientific). All samples were run in triplicates and negative controls (no template) were included in all cases. Primers used in this study were designed using the Primer-BLAST program (http://www.ncbi.nlm.nih.gov/tools/primer-blast) and synthesized by Eurofins MWG Operon (Louiswille, KY, USA).
Relative expression was calculated using GAPDH as endogenous control. Expression levels of Ecadvar mRNA were determined with primers originally designed for its identification : forward: 5′-GACCAAGTGACCACCTTAGA-3′, reverse: 5′GACCACCGCTCTCTTAGC3 ′. When indicated, the MCF7pcDNA3 cell line or the MCF7Ecadvar cells transfected with the scramble siRNA (MCF7Ecadvar siScr) were used as references as specifically indicated. The calculation used to describe these relations was 2-ΔΔCt, where ΔΔCt = [ΔCt test sample − ΔCt reference sample] and ΔCt = [Ct gene under study − Ct endogenous gene].
One-dimension protein electrophoresis and Western immunoblotting
For total protein extraction, cell cultures were washed twice at 4 °C with PBS (pH 7.4) supplemented with 2 mM CaCl2 and placed at − 20 °C until processing. Cells were then resuspended in Laemmli sample buffer, boiled for 5 min, and sonicated three times for 30 s at maximum power (Sonifier Cell Disruptor, model W-140; Heat Systems-Ultrasonics Inc., Plainview, NY, USA). Cell lysates were centrifuged 30 min at 10,000×g at 4 °C to eliminate cellular debris and stored at − 70 °C until used. Protein concentration was determined by the Bradford Protein Assay (Bio-Rad). Protein mixtures of cellular lysates were supplemented with 5% β-mercaptoethanol final concentration, boiled for 10 min, and subjected to SDS-PAGE in 8 or 10% polyacrylamide gels. Molecular weight standards (Bio-Rad) were included in all runs. After Western immunoblotting, a replica of protein patterns was obtained on nitrocellulose membranes (Hybond-ECL; GE Healthcare) using standard procedures. Membranes were placed for 1 h at room temperature with PBS containing 0.02% Tween-20 and 5% skimmed milk powder (blocking buffer), followed by overnight incubation at 4 °C with the specific primary antibody diluted in blocking buffer. Blots were washed, placed for 1 h at room temperature with the secondary antibody in blocking buffer, and developed with enhanced chemiluminescence (ECL kit; GE Healthcare), following the procedure suggested by the manufacturer. Replicates of three experiments were obtained and a densitometric analysis of the bands was performed using the ImageJ software (Wright Cell Imaging Facility, UHNR, CA, USA), when indicated. A representative image of each experiment is shown.
MCF7Ecadvar and MCF7pcDNA3 cells grown on glass coverslips until 80% confluence were fixed with 4% (v/v) paraformaldehyde in PBS for 10 min and permeabilized in PBS supplemented with 0.2% (v/v) Triton X-100. For fluorescence immunocytochemistry, non-specific binding sites were blocked with 4% (w/v) bovine serum albumin in PBS followed by overnight incubation at 4 °C with an anti LDHB antibody (2 μg/mL) or mouse IgG (control). After several washes, cells were incubated with purified secondary antibody coupled to Cy3-fluorophor for an additional hour. Nuclear staining was done by incubation with Hoechst 33342 (Sigma). Coverslips were mounted with Vectashield anti-fade solution (Vector). Images were acquired with a Nikon laser confocal microscope C1 (Nikon, Tokyo, Japan; excitation lines: 488 nm and 544 nm, emission filters: 515–530 nm and 570-LP nm) using a × 20 objective and were analyzed using the ImageJ software when specified.
Assessment of glucose and lactate concentrations
Conditioned media of MCF7Ecadvar and MCF7pcDNA3 cell cultures were collected and subjected to assessment of glucose and lactate concentrations (mg/dL). Briefly, cells were grown to 70% confluence as previously described. Cultures were washed twice with PBS and media was replaced by fresh medium supplemented with 2% fetal bovine serum. After 72 h incubation, the media were collected in sterility, centrifuged at 200×g for 10 min to remove cellular debris and used. The glucose and lactate concentrations were determined at the Dr. Enrique Rossi Diagnostic Centre, Buenos Aires, Argentina, using a Cobas 8000 equipment (Roche Diagnostics, Basel, Switzerland).
Cell treatment with 2DG
To evaluate the effect of 2DG (Sigma D6134) upon growth of MCF7Ecadvar and control cell cultures, 5000 cells were seeded in triplicates in a 96-well plate in DMEM medium freed of phenol red containing 5% FBS. After 18–20 h, culture medium was replaced by the same medium containing 2DG at 0.5 mM or 5 mM final concentrations and incubated for 96 h at 37 °C in 5% CO2 in air. Cells incubated in the absence of 2DG were also included and considered as the control condition. At the end of the incubation, cell viability was evaluated with the CellTiter 96 AQueous Non-Radioactive Cell Proliferation Assay, adding the reagent and reading the absorbance after 4-h incubation as described above.
All experiments were run in triplicates. Data are presented as the mean ± standard error of the mean (SEM). Evaluations done on MCF7Ecadvar and MCF7pcDNA3 cells were carried out using the Mann-Whitney test or the two-way ANOVA test followed by Bonferroni’s post-hoc test analysis. Studies performed on human samples were done using the Fisher’s exact test or the Mann-Whitney test, as detailed. A p < 0.05 was considered significant for each comparison. Statistical analyses and graphs were performed using the GraphPad Prism version 5.01 software (GraphPad Software, La Jolla, CA, USA).
Two-dimensional differential gel electrophoresis
N-protected amino acids
Actin, aortic smooth muscle
Actin, cytoplasmic 1
Disintegrin and metalloproteinase domain-containing protein 10
Disintegrin and metalloproteinase domain-containing protein 17
RAC-alpha serine/threonine-protein kinase
Fructose-bisphosphate aldolase A
Analysis of variance
Apoptosis regulator Bcl-2
Chloride intracellular channel protein 1
Cellular retinoic acid-binding protein 2
- E-cadherin wt:
Epithelial cadherin wild type
Enoyl-CoA delta isomerase 1, mitochondrial
Elongation factor 2
Elongation factor 2
Epidermal growth factor
Epidermal growth factor receptor
Echinoderm microtubule-associated protein-like 2
Extracellular signal-regulated kinase
Endoplasmic reticulum resident protein 29
- FADH2 :
Flavin adenine dinucleotide hydroquinone form
Fetal calf serum
Hepatocyte nuclear factor 3-alpha
Glyceraldehyde 3-phosphate dehydrogenase
Rab GDP dissociation inhibitor beta
Golgi-associated plant pathogenesis-related protein 1
Glycerol-3-phosphate dehydrogenase, mitochondrial
Growth factor receptor-bound protein 2
Stress-70 protein, mitochondrial
Hypoxia-inducible factor 1-alpha
High mobility group protein HMG-I/HMG-Y
High mobility group protein HMGI-C
Heterogeneous nuclear ribonucleoprotein H
Heterogeneous nuclear ribonucleoprotein F
Heat shock 70 kDa protein 1A
Heat shock 70 kDa protein 1B
Heat shock protein HSP 90-alpha
Heat shock protein HSP 90-beta
Heat shock protein HSP 90-alpha
Heat shock protein HSP 90-beta
Insulin-like growth factor 1 receptor
Ingenuity pathway analysis
Liquid chromatography-electrospray ionization-ms in tandem
Lactate dehydrogenase B chain
Matrix-assisted laser desorption ionization-time of flight
NADP-dependent malic enzyme
MCF-7 human breast cancer cells stably transfected with a pcDNA3 eukaryotic plasmid containing an E-cadherin variant mRNA
MCF-7 human breast cancer cells stably transfected with a pcDNA3 empty plasmid
Monocarboxylate transporter 1
Monocarboxylate transporter 4
NADP-dependent malic enzyme
Dual specificity mitogen-activated protein kinase
Mammalian target of rapamycin
- NAD+ :
Nicotinamide adenine dinucleotide, oxidized form
Nicotinamide adenine dinucleotide, reduced form
Nicotinamide adenine dinucleotide phosphate
Nucleoside diphosphate kinase B
Nuclear factor kappa-light-chain-enhancer of activated B cells
Nonsense-Mediated mRNA Decay
Protein ANalysis THrough Evolutionary Relationships
Polymerase chain reaction
Protein disulfide-isomerase A3
Protein disulfide-isomerase A6
Phosphoglycerate mutase 1
Protein kinase C
Phospholipase C gamma
Peroxisome proliferator-activated Receiver α
PR domain zinc finger protein 1
Proteasome subunit alpha type-6
Proteasome subunit beta type-3
Proteasome subunit beta type-3
Cellular retinoic acid-binding protein 2
E3 ubiquitin-protein ligase RING
Retinoid X receiver α
Standard error of the mean
small interfering RNA
Superoxide dismutase [Cu-Zn]
Signal transducer and activator of transcription
Transitional endoplasmic reticulum ATPase
Transforming growth factor
Cellular tumor antigen p53
Ubiquitin-like modifier-activating enzyme 1
14-3-3 protein theta
14-3-3 protein zeta/delta
14-3-3 protein theta
14-3-3 protein zeta/delta
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359–86.
van Roy F, Berx G. The cell-cell adhesion molecule E-cadherin. Cell Mol Life Sci. 2008;65:3756–88.
Wong SHM, Fang CM, Chuah LH, Leong CO, Ngai SC. E-cadherin: its dysregulation in carcinogenesis and clinical implications. Crit Rev Oncol Hematol. 2018;121:11–22.
Birchmeier W, Behrens J. Cadherin expression in carcinomas: role in the formation of cell junctions and the prevention of invasiveness. Biochim Biophys Acta. 1994;1198:11–26.
Prieto-García E, Díaz-García CV, García-Ruiz I, Agulló-Ortuño MT. Epithelial-to-mesenchymal transition in tumor progression. Med Oncol. 2017;34:122.
Hirohashi S. Inactivation of the E-cadherin-mediated cell adhesion system in human cancers. Am J Pathol. 1998;153:333–9.
Lee JY, Kong G. Roles and epigenetic regulation of epithelial-mesenchymal transition and its transcription factors in cancer initiation and progression. Cell Mol Life Sci. 2016;73:4643–60.
Biamonti G, Bonomi S, Gallo S, Ghigna C. Making alternative splicing decisions during epithelial-to-mesenchymal transition (EMT). Cell Mol Life Sci. 2012;69:2515–26.
Pradella D, Naro C, Sette C, Ghigna C. EMT and stemness: flexible processes tuned by alternative splicing in development and cancer progression. Mol Cancer. 2017;16:8.
Zhang J, Manley JL. Misregulation of pre-mRNA alternative splicing in cancer. Cancer Discov. 2013;3:1228–37.
Pagliarini V, Naro C, Sette C. Splicing regulation: a molecular device to enhance cancer cell adaptation. Biomed Res Int. 2015;2015:543067.
Matos ML, Lapyckyj L, Rosso M, Besso MJ, Mencucci MV, Marin-Briggiler CI, Giustina S, Furlong LI, Vazquez-Levin MH. Identification of a novel human E-cadherin splice variant and assessment of its effects upon EMT-related events. J Cell Physiol. 2017;232:1368–86.
Tu Y, Chen C, Pan J, Xu J, Zhou ZG, Wang CY. The ubiquitin proteasome pathway (UPP) in the regulation of cell cycle control and DNA damage repair and its implication in tumorigenesis. Int J Clin Exp Pathol. 2012;5:726–38.
Tavares MR, Pavan IC, Amaral CL, Meneguello L, Luchessi AD, Simabuco FM. The S6K protein family in health and disease. Life Sci. 2015;131:1–10.
Gardino AK. Yaffe MB.14-3-3 proteins as signaling integration points for cell cycle control and apoptosis. Semin Cell Dev Biol. 2011;22:688–95.
Moutsatsou P, Papavassiliou AG. The glucocorticoid receptor signaling in breast cancer. J Cell Mol Med. 2008;12:145–63.
Sacta MA, Chinenov Y, Rogatsky I. Glucocorticoid signaling: an update from a genomic perspective. Annu Rev Physiol. 2016;78:155–80.
Hermidorff MM, de Assis LV, Isoldi MC. Genomic and rapid effects of aldosterone: what we know and do not know thus far. Heart Fail Rev. 2017;22:65–89.
Pyper SR, Viswakarma N, Yu S, Reddy JK. PPARalpha: energy combustion, hypolipidemia, inflammation and cancer. Nucl Recept Signal. 2010;8:e002.
Dressman MA, Baras A, Malinowski R, Alvis LB, Kwon I, Walz TM, Polymeropoulos MH. Gene expression profiling detects gene amplification and differentiates tumor types in breastcancer. Cancer Res. 2003;63:2194–9.
Gemoll T, Habermann JK, Becker S, Szymczak S, Upender MB, Bruch HP, Hellman U, Ried T, Auer G, Jörnvall H, Roblick UJ. Chromosomal aneuploidy affects the global proteome equilibrium of colorectal cancer cells. Anal Cell Pathol (Amst). 2013;36:149–61.
Zhu YX, Tiedemann R, Shi CX, Yin H, Schmidt JE, Bruins LA, Keats JJ, Braggio E, Sereduk C, Mousses S, Stewart AK. RNAi screen of the druggable genome identifies modulators of proteasome inhibitor sensitivity in myeloma including CDK5. Blood. 2011;117:3847–57.
Leiblich A, Cross SS, Catto JW, Phillips JT, Leung HY, Hamdy FC, Rehman I. Lactate dehydrogenase-B is silenced by promoter hypermethylation in human prostate cancer. Oncogene. 2006;25:2953–60.
Liao AC, Li CF, Shen KH, Chien LH, Huang HY, Wu TF. Loss of lactate dehydrogenase B subunit expression is correlated with tumour progression and independently predicts inferior disease-specific survival in urinary bladder urothelial carcinoma. Pathology. 2011;43:707–12.
Chen R, Zhou X, Yu Z, Liu J, Huang G. Low expression of LDHB correlates with unfavorable survival in hepatocellular carcinoma: strobe-compliant article. Medicine (Baltimore). 2015;94:e1583.
Luo Y, Yang Z, Li D, Liu Z, Yang L, Zou Q, Yuan Y. LDHB and FABP4 are associated with progression and poor prognosis of pancreatic ductal adenocarcinomas. Appl Immunohistochem Mol Morphol. 2017;25:351–7.
Koh YW, Lee SJ, Park SY. Prognostic significance of lactate dehydrogenase B according to histologic type of non-small-cell lung cancer and its association with serum lactate dehydrogenase. Pathol Res Pract. 2017;213:1134–8.
Li C, Chen Y, Bai P, Wang J, Liu Z, Wang T, Cai Q. LDHB may be a significant predictor of poor prognosis in osteosarcoma. Am J Transl Res. 2016;8:4831–43.
McCleland ML, Adler AS, Deming L, Cosino E, Lee L, Blackwood EM, Solon M, Tao J, Li L, Shames D, Jackson E, Forrest WF, Firestein R. Lactate dehydrogenase B is required for the growth of KRAS-dependent lung adenocarcinomas. Clin Cancer Res. 2013;19:773–84.
McCleland ML, Adler AS, Shang Y, Hunsaker T, Truong T, Peterson D, Torres E, Li L, Haley B, Stephan JP, Belvin M, Hatzivassiliou G, Blackwood EM, et al. An integrated genomic screen identifies LDHB as an essential gene for triple-negative breast cancer. Cancer Res. 2012;72:5812–23.
Dennison JB, Molina JR, Mitra S, González-Angulo AM, Balko JM, Kuba MG, Sanders ME, Pinto JA, Gómez HL, Arteaga CL, Brown RE, Mills GB. Lactate dehydrogenase B: a metabolic marker of response to neoadjuvant chemotherapy inbreast cancer. Clin Cancer Res. 2013;19:3703–13.
Seyfried TN, Shelton LM. Cancer as a metabolic disease. Nutr Metab (Lond). 2010;7:7.
Koppenol WH, Bounds PL, Dang CV. Otto Warburg's contributions to current concepts of cancer metabolism. Nat Rev Cancer. 2011;11:325–37.
Aft RL, Zhang FW, Gius D. Evaluation of 2-deoxy-D-glucose as a chemotherapeutic agent: mechanism of cell death. Br J Cancer. 2002;87:805–12.
Cusatia MP, De la Muela MH, Hernaezb DH, Dionisioa MC, Guindoa AR, De Santiago Garcia FJ. Correlación entre la expresión de Ki67 con factores clásicos pronósticos y predictivos en el cáncer de mama precoz. Rev Senol Patol Mamar. 2014;27:163–9.
Shamir ER, Ewald AJ. Adhesion in mammary development: novel roles for E-cadherin in individual and collective cell migration. Curr Top Dev Biol. 2015;112:353–82.
Baranwal S, Alahari SK. Molecular mechanisms controlling E-cadherin expression in breast cancer. Biochem Biophys Res Commun. 2009;384:6–11.
Ashaie MA, Chowdhury EH. Cadherins: the superfamily critically involved in breast cancer. Curr Pharm Des. 2016;22:616–38.
Meleady P. Two-dimensional gel electrophoresis and 2D-DIGE. Methods Mol Biol. 2018;1664:3–14.
Arentz G, Weiland F, Oehler MK, Hoffmann P. State of the art of 2D DIGE. Proteomics Clin Appl. 2015;9(3–4):277–88.
Monge M, Doll A, Colas E, Gil-Moreno A, Castellvi J, Garcia A, Colome N, Perez-Benavente A, Pedrola N, Lopez-Lopez R, Dolcet X. Ramon y Cajal S, Xercavins J. et al. subtractive proteomic approach to the endometrial carcinoma invasion front. J Proteome Res. 2009;8:4676–84.
Mathias RA, Simpson RJ. Towards understanding epithelial-mesenchymal transition: a proteomics perspective. Biochim Biophys Acta. 2009;1794:1325–31.
Biarc J, Gonzalo P, Mikaelian I, Fattet L, Deygas M, Gillet G, Lemoine J, Rimokh R. Combination of a discovery LC-MS/MS analysis and a label-free quantification for the characterization of an epithelial-mesenchymal transition signature. J Proteome. 2014;110:183–94.
Bottoni P, Isgrò MA, Scatena R. The epithelial-mesenchymal transition in cancer: a potential critical topic for translationalproteomic research. Expert Rev Proteomics. 2016;13:115–33.
Vergara D, Simeone P, Franck J, Trerotola M, Giudetti A, Capobianco L, Tinelli A, Bellomo C, Fournier I, Gaballo A, Alberti S, Salzet M, Maffia M. Translating epithelial mesenchymal transition markers into the clinic: novel insights from proteomics. EuPA Open Proteom. 2016;10:31–41.
Vergara D, Simeone P, Latorre D, Cascione F, Leporatti S, Trerotola M, Giudetti AM, Capobianco L, Lunetti P, Rizzello A, Rinaldi R, Alberti S, Maffia M. Proteomics analysis of E-cadherin knockdown in epithelial breast cancer cells. J Biotechnol. 2015;202:3–11.
Sreedhar A, Zhao Y. Dysregulated metabolic enzymes and metabolic reprogramming in cancer cells. Biomed Rep. 2018;8:3–10.
Warburg O, Wind F, Negelein E. The metabolism of tumors in the body. J Gen Physiol. 1927;8:519–30.
Wise DR, Ward PS, Shay JE, Cross JR, Gruber JJ, Sachdeva UM, Platt JM, DeMatteo RG, Simon MC, Thompson CB. Hypoxia promotes isocitrate dehydrogenase-dependent carboxylation of α-ketoglutarate to citrate to support cell growth and viability. Proc Natl Acad Sci U S A. 2011;108:19611–6.
Terunuma A, Putluri N, Mishra P, Mathe EA, Dorsey TH, Yi M, Wallace TA, Issaq HJ, Zhou M, Killian JK, Stevenson HS, Karoly ED, Chan K, et al. MYC-driven accumulation of 2- hydroxyglutarate is associated with breast cancer prognosis. J Clin Invest. 2014;124:398–412.
Goetzman ES, Prochownik EV. The role for Myc in coordinating glycolysis, oxidative phosphorylation, glutaminolysis, and fatty acid metabolism in Normal and neoplastic tissues. Front Endocrinol (Lausanne). 2018;9:129.
Davis NM, Sokolosky M, Stadelman K, Abrams SL, Libra M, Candido S, Nicoletti F, Polesel J, Maestro R, D'Assoro A, Drobot L, Rakus D, Gizak A, et al. Deregulation of the EGFR/PI3K/PTEN/Akt/mTORC1 pathway in breast cancer: possibilities for therapeutic intervention. Oncotarget. 2014;5:4603–50.
Wee P, Wang Z. Epidermal growth factor receptor cell proliferation signaling pathways. Cancers (Basel). 2017;9:52.
Wallace JA, Li F, Leone G, Ostrowski MC. Pten in the breast tumor microenvironment: modeling tumor-stroma coevolution. Cancer Res. 2011;71:1203–7.
Moulder DE, Hatoum D, Tay E, Lin Y, McGowan EM. The roles of p53 in mitochondrial dynamics and cancer metabolism: the pendulum between survival and death in breast Cancer? Cancers (Basel). 2018;10:189.
Zha X, Hu Z, He S, Wang F, Shen H, Zhang H. TSC1/TSC2 inactivation inhibits AKT through mTORC1-dependent up-regulation of STAT3-PTEN cascade. Cancer Lett. 2011;313:211–7.
Choi J, Jung WH, Koo JS. Metabolism-related proteins are differentially expressed according to the molecular subtype of invasive breast cancer defined by surrogate immunohistochemistry. Pathobiology. 2013;80:41–52.
Long JP, Li XN, Zhang F. Targeting metabolism in breast cancer: how far we can go? World J Clin Oncol. 2016;7:122–30.
Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, Thomas PD. PANTHER version 11: expanded annotation data from gene ontology and reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2017;45:D183–9.
PANTHER Classification System, 13.1 version. http://pantherdb.org. Accessed Mar 2018.
DisGeNET, v5.0 version. http://www.disgenet.org. Accessed Mar 2018.
Piñero J, Queralt-Rosinach N, Bravo À, Deu-Pons J, Bauer-Mehren A, Baron M, Sanz F, Furlong LI. DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes. Database (Oxford). 2015;2015:bav028.
VENNY, 2.1 version. http://bioinfogp.cnb.csic.es/tools/venny. Accessed Apr 2018.
Ingenuity Pathway Analysis, Qiagen. http://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/. Accessed Oct 2015.
García-Ramírez M, Canals F, Hernández C, Colomé N, Ferrer C, Carrasco E, García-Arumí J, Simó R. Proteomic analysis of human vitreous fluid by fluorescence-based difference gel electrophoresis (DIGE): a new strategy for identifying potential candidates in the pathogenesis of proliferative diabetic retinopathy. Diabetologia. 2007;50:1294–303.
Shevchenko A, Wilm M, Vorm O, Mann M. Mass spectrometric sequencing of proteins silver-stained polyacrylamide gels. Anal Chem. 1996;68:850–8.
We sincerely thank Mariana Regueira, PhD, for providing the 2DG. In addition, we thank Dr. Liliana Maggi and Dr. Alejandro Celia, from Centro de Diagnóstico Rossi, for lactate and glucose concentrations measurements. Moreover, we thank M.Sc. María Florencia Abascal, for technical assistance in bioinformatics analysis. Finally, we thank Dr. Roberto Ortí and Dr. Alejandra Wernicke from the Hospital Italiano of Buenos Aires, and Dr. Pilar Carballo, Dr. Liliana Gimenez, and MSc. María de Nazareth Juárez Rusjan from the Instituto Roffo Tissue Bank, for providing us with the BC tissue samples evaluated in this study.
Studies were supported by grants from the Agencia Nacional de Promoción Científica y Tecnológica de Argentina [Grant PICTSU1072], Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) [Grant PIP740 and PIP887], Instituto Nacional del Cáncer [Grant INC 2014–2015 and INC 2016–2017], Fundación Roemmers [Grant 2011–2013], Fundación R. Barón (donation) and Fundación Williams (donation) to MHVL, and IRSES [International Research Staff Exchange Scheme; Ref. 269285] to JR and MHVL. The Proteomics Unit at VHIO belongs to ProteoRed, PRB2-ISCIII, supported by grant PT13/0001.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
Tissue samples from patients diagnosed with BC were obtained from the Hospital Italiano of Buenos Aires and the Instituto de Oncología Ángel H. Roffo, Buenos Aires, Argentina. Samples collection protocols were approved by the institutional Ethics Committees and patient informed consents were obtained in all cases. The study was conducted in accordance with the principles of the Declaration of Helsinki.
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.
Figure S1. E-cadherin protein expression analysis in MCF7Ecadvar and MCF7pcDNA3 cells. Western immunoblotting of the 120 kDa E-cadherin wild type (E-cadherin wt) full length form in total protein extracts obtained from 3 MCF7pcDNA3 (1 to 3, 10 μg per lane) and 3 MCF7Ecadvar (4 to 6, 20 μg per lane) cell cultures. The immunodetection of β-tubulin (55 kDa) was used as loaded control. (TIF 2449 kb)
Figure S2. Cellular components and protein classes associated with the 50 proteins identified. Results obtained with PANTHER. (A) Column graph bar in which the percentage (%) of each cellular component was determined from the number of proteins included in each category (cell part 52.5%, organelle 26.1%, macromolecular complex 10.9%, extracellular region 6.5% and membrane 4.3%). A table including the number and symbol of the proteins involved in each protein class category is also shown. (B) Column graph bar in which the percentage (%) of each protein class was determined from the number of proteins included in each category (oxidoreductase 19.6%; hydrolase 15.2%; isomerase 10.9%; chaperone 8.7%; transferase 8.7%; cytoskeletal protein 6.5%; nucleic acid binding 6.5%; ligase, enzyme modulator 4.3% and lyase, 4.3% each; calcium-binding protein, membrane-traffic protein signaling molecule, transfer/carrier protein and transporter, 2.2% each). A table including the number and symbol of the proteins involved in each protein class category is also shown. (TIF 2499 kb)
Figure S3. Genes associated with human BC. Results obtained with DisGeNET. List of the 5261 human genes that emerged from the 36 terms (“Breast Carcinoma”, “Female Breast Carcinoma”, “Stage 0 Breast Carcinoma”, “Stage IIIA Breast Carcinoma”, “Stage IIIB Breast Carcinoma”, “Invasive Ductal Breast Carcinoma”, “Invasive Lobular Breast Carcinoma”, “Secretory Breast Carcinoma”, “Inflammatory Breast Carcinoma”, “Adenoid Cystic Breast Carcinoma”, “Apocrine Breast Carcinoma”, “Invasive Apocrine Breast Carcinoma”, “Intermediate Grade Ductal Breast Carcinoma In Situ”, “Breast Carcinoma Metastatic in the Skin”, “Breast Cancer 3”, “Breast Cancer Stage II”, “Stage III Breast Cancer AJCC v6”, “Breast Cancer Recurrent”, “Bilateral Breast Cancer”, “Breast Cancer and Pregnancy”, “Breast Cancer, Familial”, “Breast Cancer (non-specific) Premenopausal”, “Contralateral Breast Cancer”, “Unilateral Breast Neoplasms”, “Malignant Neoplasm of Breast”, “Malignant Neoplasm of Female Breast”, “Malignant Neoplasm of Breast Stage I”, “Malignant Neoplasm of Breast Staging”, “Secondary Malignant Neoplasm of Female Breast”, “Triple Negative Breast Neoplasms”, “Mammary Carcinoma, Human”, “Mammary Ductal Carcinoma”, “Mammary Neoplasms”, “Mammary Neoplasms, Human”, “Mammary Neoplasms, Experimental” and “Mammary Tumorigenesis”) found in DisGeNET containing the words “Breast” or “Mammary”, and “Carcinoma”, “Cancer”, “Neoplasms” or “Tumorigenesis”. The 39 genes in common with those that code for the 50 proteins identified by 2D-DIGE and MS as differentially expressed between the MCF7Ecadvar and MCF7pcDNA3 cell lines are highlighted. (PDF 168 kb)
Figure S4. Specific primers for the top-10 upregulated and downregulated molecules among the 50 identified. List of the primers used for the evaluation of the expression levels of the mRNAs that code for the 10 most (A) upregulated and (B) downregulated proteins identified with differential expression levels between MCF7Ecadvar and MCF7pcDNA3 cells by 2D-DIGE and MS. The sequence of the forward and reverse primers, as well as the size of each amplified product, are indicated. (TIF 2455 kb)
Figure S5. Transcripts expression analysis of MCT1 and MCT4 in MCF7Ecadvar and MCF7pcDNA3 cells. Quantitative expression analysis of (A) MCT1 and (B) MCT4 lactate transporters by real time PCR in MCF7pcDNA3 and MCF7Ecadvar cells. The relative expression was calculated as described in the Materials and Methods section, using GAPDH as the endogenous gene and the MCF7pcDNA3 cell line as reference. *p < 0.05. (TIF 4993 kb)
Table S1. Histological and molecular characteristics of breast tumor tissue samples. Description table in which the histological type and grade, the presence or absence of metastasis, the expression pattern of ER, PR and HER2, the percentage (%) of Ki-67 staining, as well as the molecular subtype of 21 breast tumor tissues, are detailed. (DOCX 17 kb)
Table S2. Transcripts expression analysis of E-cadherin variant and LDHB in human breast tumor samples. Quantitative expression analysis of E-cadherin variant and LDHB by real time PCR in 21 human breast tumor tissues. Samples were numbered and ordered in an ascending form according with the 2-ΔCt values obtained for the E-cadherin variant mRNA. (DOCX 15 kb)
About this article
Cite this article
Rosso, M., Lapyckyj, L., Besso, M.J. et al. Characterization of the molecular changes associated with the overexpression of a novel epithelial cadherin splice variant mRNA in a breast cancer model using proteomics and bioinformatics approaches: identification of changes in cell metabolism and an increased expression of lactate dehydrogenase B. Cancer Metab 7, 5 (2019). https://doi.org/10.1186/s40170-019-0196-9
- Breast cancer
- Epithelial cadherin
- Alternative splicing
- Epithelial to mesenchymal transition
- Proteomic analysis
- Mass spectrometry
- Lactate dehydrogenase B