Abstract
One of the earliest events in epithelial carcinogenesis is the dissolution of tight junctions and cell polarity signals that are essential for normal epithelial barrier function. Here, we report that EFA6B, a guanine nucleotide exchange factor for the Ras superfamily protein Arf6 that helps assemble and stabilize tight junction, is required to maintain apico-basal cell polarity and mesenchymal phenotypes in mammary epithelial cells. In organotypic three-dimensional cell cultures, endogenous levels of EFA6B were critical to determine epithelial–mesenchymal status. EFA6B downregulation correlated with a mesenchymal phenotype and ectopic expression of EFA6B hampered TGFβ-induced epithelial-to-mesenchymal transition (EMT). Transcriptomic and immunohistochemical analyses of human breast tumors revealed that the reduced expression of EFA6B was associated with loss of tight junction components and with increased signatures of EMT, cancer stemness, and poor prognosis. Accordingly, tumors with low levels of EFA6B were enriched in the aggressive triple-negative and claudin-low breast cancer subtypes. Our results identify EFA6B as a novel antagonist in breast cancer and they point to its regulatory and signaling pathways as rational therapeutic targets in aggressive forms of this disease. Cancer Res; 74(19); 5493–506. ©2014 AACR.
Introduction
Gene expression profiling has helped define five intrinsic tumor subtypes; essentially based on the expression of the two hormonal receptors, estrogen receptor (ER) and progesterone receptor (PR), and HER2. The triple-negative breast cancer (TNBC) subtype does not express any of these receptors, rendering it unsuitable for the hormonal and anti-HER2 treatments (1). TNBC is associated with low overall survival and high recurrence. Transcriptomic analyses have identified a new breast cancer subtype with poor prognosis named claudin-low (Cld-low), defined by the loss of expression of the proteins constitutive of the tight junction (2, 3). This subtype, which represents 12% of total breast cancer, is 75% TNBC, characterized by a core epithelial-to-mesenchymal transition (EMT) signature, and enriched in cancer stem cells (CSC), which are associated with relative resistance to conventional anticancer treatments.
We have reported that exchange factor for Arf6 (EFA6) contributes to the assembly and maintenance of the tight junction in the canine kidney cell line MDCK (4, 5). Four isoforms of EFA6 (A, B, C, and D) are encoded by separate genes (PSD, PSD4, PSD2, and PSD3, respectively). EFA6A, B, and D have ubiquitous expression while EFA6C is restricted to neuronal cells. We found that EFA6B levels increase at newly formed cell–cell contacts during epithelial polarity development to help promote tight junction assembly (4, 5). Once tight junctions are assembled, the levels of EFA6B decreases, yet at steady-state EFA6B still contributes to stabilize the tight junctions. Knockdown of EFA6B is sufficient to impair the assembly and stability of the tight junctions (5). Hence, the regulation of the levels of EFA6B proteins is key to tight junction homeostasis.
In this study, we have asked whether the EFA6B protein could act as an antagonist to malignant progression. To answer this question, we combined in vivo expression analyses in human breast tumor samples with in vitro approaches. On the basis of our results, we propose that EFA6B acts as an antagonist at the early stages of breast cancer development by hampering tight junction disassembly and loss of epithelial polarity.
Materials and Methods
Cells
MCF7, T47D, SKBR3, and MCF10A cells were obtained from the ATCC, authenticated by short-tandem repeat profiling by the vendor, and passaged for fewer than 6 months before experiments. Cells were grown in DMEM containing 10% FCS, except MCF10A cells grown in DMEM-F12, 5% horse serum, insulin, hydrocortisone, EGF, and cholera toxin (6). Transient transfection was performed by nucleofection (AMAXA) and stable transfection by lipofection using JETPEI (Polyplus Transfection). For stable expression of EFA6Bvsvg, cells were selected with Geneticin (Invitrogen) 2 days after transfection and low expressor clones were isolated. They all gave similar results in two-dimensional (2D) and three-dimensional (3D) culture assays. One clone was used in all the assays, and two other clones when randomly used in the different assays displayed similar results. When indicated, cells were grown on permeable filters (Costar) at 3 × 105 cells/12 mm filter. For 3D-cell culture, 2 × 103 cells were mixed with 25 μL of Matrigel (BD Biosciences) deposited as a drop on a 12-mm glass coverslip. MDCK-EFA6Bvsvg cells were described elsewhere (5). Before TGFβ treatment, the cells were serum-starved overnight and 10 μg/mL TGFβ (Sigma-Aldrich) was added for the indicated period of times.
DNA, siRNAs, and shRNAs
pEGFP-Arf6T27N and pcDNA3-EFA6Bvsvg were described elsewhere (7, 8). EFA6B was cloned into pEGFPC3 by PCR. The silencing mutations to generate a construct insensitive to the siRNAs were introduced using the QuickChange Mutagenesis Kit (Agilent Technologies). Extinction of EFA6B expression was obtained in MCF7 cells using a pool of four siRNA or in MCF10A, T47D, and SKNBR3 cells using a pool of three human EFA6B-specific shRNA cloned into the lentiviral vector pLB, independently expressing GFP (Addgene; #11619). Sequences of all the primers, siRNAs, and shRNAs are described in Supplementary Methods. The pWZL-Blast-mouse-E-cadherin lentiviral vector was obtained from Addgene (#18804; ref. 9).
Antibodies
Rabbit polyclonal sera specific for ZO-1, claudin-1, occludin (Invitrogen), EFA6B (Sigma-Aldrich), Smad2, and phospho-Smad2 (Ser465/467; Cell Signaling Technology), mouse monoclonals specific for E-cadherin, β-catenin, GM130 (BD Biosciences), RhoA (Santa Cruz Biotechnology), actin, and vsvg tag (Sigma-Aldrich), rabbit monoclonal specific for phospho-Ezrin/Radixin/Moesin (Cell Signaling Technology) were used. The mouse monoclonal anti-Arf6 was described elsewhere (10).
Transepithelial electrical resistance
The transepithelial electrical resistance (TER) was measured with an EVOM (WPI) using triplicates or quadruplicates for each measurement as previously described (11).
Immunofluorescence
Cells were fixed in paraformaldehyde, the samples processed as previously described (4) and imaged on a confocal microscope (Leica TCS-SP5). Tight junction length and fluorescence intensity quantifications were performed on the basis of the occludin fluorescence staining using a home-made ImageJ program described in details in Supplementary Methods (12).
Immunoblot
Cells were solubilized in SDS lysis buffer and the lysates were prepared for immunoblot as previously described (5). The proteins were revealed by chemiluminescence (Amersham) and the membranes were analyzed with the luminescent analyzer LAS-3000 (Fujifilm).
EFA6B/PSD4 mRNA expression analysis in breast cancer samples
To determine EFA6B/PSD4 mRNA expression in breast cancer samples and search for correlations with histoclinical and molecular features, we analyzed DNA microarray–based transcriptomic data. Expression in mammary cell lines was analyzed in our database of 43 cell lines profiled using our Affymetrix platform (U133 Plus 2.0; ref. 13). Expression in clinical tissue samples was analyzed in our series of 352 patients with invasive adenocarcinoma and four pools of normal breast tissue samples (11 healthy women), and in 32 public datasets collected from the National Center for Biotechnology Information/Genbank Gene Expression Omnibus (GEO) database, the European Bioinformatics Institute ArrayExpress database, or at the authors' websites (Supplementary Table S1). This resulted in a total of 5,252 nonredundant primary invasive breast cancers with both EFA6B/PSD4 mRNA expression and histoclinical data available for analysis. For the Agilent-based datasets, we applied quantile normalization to available processed data. Regarding the Affymetrix-based datasets, we used robust multichip average (RMA) with the nonparametric quantile algorithm as normalization parameter. Comparison of mean EFA6B/PSD4 mRNA expression level according to classical histoclinical factors was done using Student t test. To compare the distribution according to categorical variables, we used the Fisher exact test. All statistical tests were two-sided at the 5% level of significance. Statistical analysis was done using the survival package (version 2.30) in the R software. We followed the reporting REcommendations for tumor MARKer prognostic studies (REMARK criteria). A fully detailed description of the methods and statistical procedures is provided in Supplementary Methods.
Immunohistochemical expression of EFA6B and claudin-3 in breast cancer samples
A total of 412 human primary invasive breast adenocarcinoma were analyzed by immunohistochemistry (IHC) using tissue microarray (TMA). Using a TMA instrument (Beecher Instruments), three representative areas of each tumor were punched from the archival paraffin-embedded formalin-fixed tissue blocks. IHC was performed on serial 4-μm deparaffinized TMA sections using an automated single-staining procedure. The extent of EFA6B and claudin-3 immunostaining signals was expressed semi-quantitatively as the sum of scores representing the proportion and staining intensity of negative and positive tumor cells. The score of the immunoreactivity for the two proteins was tested for correlation with traditional clinicopathologic parameters by using the χ2 test for categorical variables, and the Mann–Whitney U test for continuous variables. The coefficient of correlation (ρ) between variables was calculated using the Spearman rank test. Association between protein levels and survival was calculated using the Kaplan–Meier method and curves were compared with the log-rank test. The collection of breast cancer samples, detailed methods, and statistical procedures are described in Supplementary Methods.
Results
Exogenous expression of EFA6B in MCF7 cells restores a normal polarized epithelial monolayer with functional tight junction
To address the role of EFA6B on breast cancer development, we first asked whether the exogenous expression of EFA6B would enable tumoral mammary cells to assemble a functional tight junction. We chose to study the MCF7 cells because it is the only human mammary tumoral cell line known to assemble, although very partially, a tight junction. Owing to their weakly tumoral phenotype, they adopt an incomplete polarized morphology when grown on filters and an unpolarized compact cell cluster organization under 3D-culture conditions (14, 15). After stable transfection of EFA6Bvsvg, we selected subclones with no more than five times the level of the endogenous protein (Fig. 1B). At this level of expression, EFA6Bvsvg remained localized at the apical pole of the cells consistent with our previous observations for the endogenous protein in MDCK cells (Fig. 1A; ref. 4). When grown to confluence on filters, about 10% of MCF7 cells organized in small patches of 5 to 10 cells that assembled a tight junction. In contrast, MCF7 cells expressing EFA6Bvsvg (MCF7-EFA6Bvsvg) demonstrated consistent assembly of a continuous tight junction throughout the entire monolayer. The adherens junctions correctly assembled and were not affected by EFA6Bvsvg expression (Fig. 1A and Supplementary Fig. S1A and S1B). The expression of EFA6Bvsvg did not alter the levels of expression of the components of the tight junction, nor those of the adherens junction, nor the substrate of EFA6B, Arf6 (Fig. 1B). An automated quantification based on occludin staining showed that tight junction size and mean intensity increased at least three and four times, respectively (Fig. 1C). The assembled tight junctions were functional as MCF7-EFA6Bvsvg gained TER faster than control cells with a 3.2-fold increase at day 3 (Fig. 1D). More strikingly, the MCF7-EFA6Bvsvg cells adopted characteristics of a well-differentiated apico-basal epithelial morphology: (i) regular cuboidal cells organized as a monolayer, (ii) apical punctate actin staining typical of the microvilli of polarized cells, (iii) nuclei regularly aligned basally, (iv) adherens junction proteins restricted to the lateral cell–cell contacts, (v) apical tight junction proteins (Fig. 1A and Supplementary Fig. S1A and S1B). In contrast, MCF7 cells displayed a heterogeneous and rather round morphology, tended to form multilayers, displayed a continuous plasma membrane staining of the adherens junctions all around the cell, and formed little tight junctions and no microvilli (Fig. 1A and Supplementary Fig. S1A and S1B). Thus, expression of EFA6Bvsvg is sufficient to elicit an epithelial transition of the MCF7 cells from their intermediate epithelio-mesenchymal phenotype.
Exogenous expression of EFA6B in MCF7 cells induces a phenotypic reversion in 3D culture
Although 2D culture on porous filters is best suited to analyze the formation and barrier function of the tight junction, 3D-organotypic cell culture where cells organize as acini is more appropriate to assess the luminogenic and polarizing capacities of epithelial cells. As previously described, in 3D culture, MCF7 cells formed aggregates with no radial symmetry and random nuclei distribution; a phenotype named “mass” representative of rather weakly tumoral mammary cell lines (14, 15). The F-actin staining was stronger in peripheral membrane protrusions (Fig. 1E), and the rare lumens were virtually always between just two cells. In contrast, EFA6Bvsvg expression stimulated lumen formation and MCF7-EFA6Bvsvg cells organized as nearly normal mammary acini-like structures (Fig. 1E and Supplementary Fig. S1C and S1D) comparable with those obtained with the normal mammary epithelial MCF10A cells (data not shown; refs. 14, 15). The quantification of the number of aggregates containing at least one lumen showed a 3-fold increase in MCF7 cells when expressing EFA6Bvsvg; few lumens were formed just in between two cells (Fig. 1F). Most lumens were connected to multiple cells to form either one central lumen (9.4%), or two large lumens (62.5%) and sometimes three or more lumens (28.1%). Regardless of the number of lumens, within the cysts each individual cell was connected to a lumen and perfectly polarized. Their apical surface facing the lumen displayed a strong actin staining colabeled for the apical phosphorylated proteins of the ERM family (Fig. 1E and Supplementary Fig. S1C and S1D). The adherens junction proteins were restricted to the lateral surface and the tight junction markers confined to the apex of the lateral membrane delineating the apical domain (Fig. 1E and Supplementary Fig. S1C and S1D). Finally, the Golgi apparatus was positioned above the nucleus, consistent with the correct orientation of the apico-basal axis (Supplementary Fig. S1D). Thus, EFA6Bvsvg expression is not only capable of facilitating the assembly of the tight junction, it also promotes the reversion to an epithelial phenotype, resulting in cells organizing collectively into normal mammary acini-like structures. This property was totally blocked upon coexpression of the dominant-negative Arf6T27N-GFP (Fig. 1G), indicating that EFA6B effects are dependent on Arf6 activation.
EFA6B contributes to the formation and fusion of VACs to generate lumens in polarized MCF7 acini
Although, MCF7 cells do not form acinis, they do form aggregates, indicating that the limiting step to make normal acinis is likely the ability to form lumens. Thus, we investigated how MCF7-EFABvsvg cell formed their lumens. When grown on coverslips and stained for F-actin, we observed that a higher proportion of cells overexpressing EFA6Bvsvg contained large intracellular vacuoles reminiscent of the vacuolar apical compartments (VAC) as well as small lumens formed in between two and four cells. Several VACs could be observed in certain cells, and VACs were already visible at the two-cell stage (Fig. 2A). VACs share characteristics with the apical plasma membrane displaying microvilli and containing apically, but not basolaterally, targeted proteins. VACs are proposed to fuse to cell–cell contacts to form the apical membrane (16–18). Quantification of cell aggregates grown in 3D culture showed that MCF7-EFA6Bvsvg contained three times more VACs and five times more lumens than control MCF7 cells (Fig. 2F). The lumens formed in MCF7-EFA6Bvsvg cells displayed the typical characteristics of an apical plasma membrane with actin-rich microvilli, selective staining for apical markers such as the P-ERM proteins (Fig. 2B), absence of the basolateral proteins β-catenin and E-cadherin (Fig. 2B and C), and surrounded by a tight junction ring (Fig. 2C). To assess EFA6B localization during luminogenesis, we generated a stable MCF7 cell line expressing EFA6BGFP. Interestingly, at early stages, EFA6BGFP was found at the cell surface, then on the VACs and at the apical plasma membrane (Fig. 2D). To confirm that the lumens were formed by fusion of the VACs with the plasma membrane, we mixed MCF7-EFA6BGFP and MCF7 control cells. We observed that whole lumens stained for F-actin and P-ERM were made of one half labeled with EFA6BGFP and the other half unlabeled. Thus, upon fusion, the VACs provide one half of the lumen, the other half being formed on the other side by the facing cell (Fig. 2E). In addition, upon fusion of the VACs, the new lumens were rapidly limited to the cell–cell contact by a ring of tight junction (Fig. 2E).
Wild-type MCF7 cells form typical VACs and small lumens albeit essentially seen in between two cells (Fig. 1F and Supplementary Fig. S1E). We asked whether the endogenous EFA6B was also necessary to generate these VACs and small lumens. We used a pool of four siRNAs that specifically downregulated EFA6B but not the other EFA6 isoforms expressed in MCF7 cells, EFA6A and EFA6D (Supplementary Fig. S2A and S2B). We observed that siRNA knockdown of EFA6B led to half fewer VACs and three times less lumens (Fig. 2G).
Collectively, our results showed that EF6AB is involved in the formation and fusion of the VACs to form lumens within MCF7 cell aggregates, providing a mechanistic explanation as to how EFA6B contributes to the epithelial phenotypic reversion of MCF7 cells in 3D culture.
EFA6B repression in MCF7 cells abrogates the assembly of tight junctions and leads to highly unorganized aggregates
Next, we analyzed the effects of EFA6B downregulation by siRNA. When grown on filters, EFA6B-depleted cells displayed a less organized phenotype with a higher propensity to form multilayers. No tight junction formation was observed, although the adherens junction was unaffected as well as the total amount of tight junction proteins (Fig. 3A–C), indicating that EFA6B is regulating the tight junction assembly per se as previously observed in MDCK cells (4, 5). When grown in Matrigel, the aggregates became much less compact and appeared multilobular reminiscent of the “grape-like” morphotype typically found in advanced tumoral cell lines (14, 15). In these unorganized aggregates, adherens junction staining was still visible while the tight junction and the apical membrane domain were absent (Fig. 3D). The proportion of grape-like aggregates was reduced upon EFA6Bvsvg expression and increased in EFA6B knockdown cells. As a consequence, the percentage of aggregates with lumen(s) was proportionally reduced in EFA6B knockdown cells (Fig. 3E and F). The specificity of EFA6B knockdown was verified by a rescue experiment in which the effects of the depletion of the endogenous EFA6B protein by siRNA were reverted by the coexpression of a siRNA-insensitive EFA6B construct (Supplementary Fig. S2C and S2D). In conclusion, while overexpression of EFA6B supports the epithelial phenotype, downregulation of endogenous EFA6B promotes the mesenchymal phenotype.
EFA6B repression disturbs the epithelial phenotype of MCF10A, T47D, and SKBR3/E-cadherin cells
We then assessed the effects of EFA6B depletion on the collective organization of other human mammary cell lines grown in Matrigel. Normal mammary epithelial cells, MCF10A, were analyzed at early-stage aggregate formation and after luminogenesis. At early stages, wild-type and pLB-transduced cells formed round structures composed of rather cuboidal cells organized as a regular peripheral layer surrounding the inner cells of the aggregates. In contrast, EFA6B knocked down cells formed irregularly shaped structures of the “mass” type made of round and heterogeneous cells, similar to MCF7 aggregates (Fig. 4A–C). At later stages, EFA6B depletion did not prevent the cells from eventually forming a single central lumen, most likely because MCF10A cells do not form their lumen by exocytosis of VACs, but by cavitation, a process dependent on apoptosis (19), which we found not to be affected by EFA6B levels of expression (data not shown). Next, we studied the T47D cell line that is similar to MCF7 (ER/PR+, HER2-low) with an intermediate epithelial–mesenchymal phenotype although they do not assemble any tight junction. As previously reported (14, 15), T47D wild-type as well as pLB-transduced cells adopted the “mass” phenotype. Similarly to MCF7, upon depletion of EFA6B, T47D cells no longer formed round and compact aggregates but rather unorganized loose structures (Fig. 4D–F). Subsequently, we studied the more tumoral cell line SKBR3 (ER/PR−, HER2-high) that does not express E-cadherin and organizes as archetypical grape-like structure (14, 15). However, upon exogenous expression of E-cadherin, SKBR3 cells were shown to organize in the “mass” phenotype (Fig. 4G–I; ref. 20). Although knockdown of EFA6B in SKBR3 cells did not affect their grape-like phenotype (data not shown), EFA6B knocked down in SKBR3/E-cadherin cells reverted the “mass” phenotype back to a grape-like phenotype. This result supports our previous observations that EFA6B acts downstream of E-cadherin to regulate the tight junction assembly and epithelial cell polarization (4, 5). Altogether, our results demonstrate that EFA6B expression is required for the proper collective organization of mammary epithelial cells.
Exogenous expression of EFA6Bvsvg abrogates the TGFβ-induced EMT
We showed that increasing EFA6Bvsvg expression supports the epithelial phenotype, whereas reciprocally repression of EFA6B promotes the mesenchymal phenotype, suggesting that EFA6B might regulate the EMT. TGFβ is a major inducer of EMT during mammary gland morphogenesis and breast cancer progression (21–23). TGFβ induces early on the disassembly of the tight junction to which the TGFβ receptor (TGFβR) is physically associated through the polarity protein Par6 and occludin (24–27). When MCF7 cells were exposed to TGFβ, they downregulated the expression of E-cadherin and occludin, without affecting the endogenous EFA6B (Fig. 5A and B). This was accompanied by a major change in cell morphology with the apparition of stress fibers and an elongated cell shape (Fig. 5C). In contrast, MCF7-EFA6Bvsvg cells did not respond to TGFβ treatment (Fig. 5A and C). The expression of EFA6Bvsvg did not alter the levels of expression of TGFβR, thus the inhibitory effect of EFA6Bvsvg was not due to a decrease in the expression of TGFβR (data not shown). The inhibitory effect of EFA6B was totally blocked upon coexpression of the dominant-negative Arf6T27N-GFP (Fig. 5D). Note that in those cells, as in MCF7 cells expressing Arf6T27N-GFP (data not shown), stress fibers were no longer observed even upon TGFβ treatment. Nevertheless, the change in morphology and drop of E-cadherin expression were clearly apparent, indicating that Arf6 activation is required to hamper the TGFβ-induced EMT. To assess whether this inhibitory effect was cell-specific, we analyzed the non-mammary epithelial MDCK cell line. Like in MCF7 cells, in MDCK cells overexpressing EFA6Bvsvg the morphologic changes and the decrease in the levels of E-cadherin and occludin in response to TGFβ were totally blocked (Supplementary Fig. S4).
Induction of EMT by TGFβ is believed to take place through the activation of the canonical SMAD pathway that transduces gene expression program, and noncanonical pathways responsible for the rapid morphologic changes including the Par6, PI3K, and ERK pathways (23, 28). We observed a steady increase in Smad2 phosphorylation in both MCF7 and MCF7-EFA6Bvsvg cells, indicating that EFA6Bvsvg did not affect the SMAD pathway (Fig. 5E). TGFβ induces the disassembly of the tight junction by downregulating the levels of RhoA through Par6 and the ubiquitin-ligase Smurf1 (26). Similarly to previous reports (29, 30), we found a transient but reproducible decrease in the levels of RhoA after 30-minute exposure to TGFβ, which was comparable in both cell lines, indicating that EFA6Bvsvg did not act through RhoA downregulation (Fig. 5E). MCF7 cells carry an activating mutation of the PIK3CA, leading to a constitutive phosphorylation of its downstream effector Akt (31, 32). This constitutive phosphorylation was unaffected by EFA6Bvsvg expression and we could not detect further phosphorylation upon TGFβ treatment in either cell lines. Similarly, ERK1/2 phosphorylation was not increased by TGFβ in MCF7 and MCF7EFA7Bvsvg cell lines (data not shown). Collectively, our data indicate that EFA6Bvsvg inhibits the TGFβ-induced EMT through a yet unidentified mechanism.
EFA6B/PSD4 mRNA expression in breast cancer
The fact that EFA6B regulates the tight junction assembly, apico-basal polarity, and inhibits the TGFβ-induced EMT led us to determine whether EFA6B might act as an antagonist to mammary tumoral progression in vivo. A recent transcriptomic analysis revealed the existence of a new subtype of breast cancer called Cld-low characterized by the loss of expression of the tight junction components (2). Thus, we postulated that loss of EFA6B expression could contribute to the downregulation of the tight junction during tumoral progression. Analysis of EFA6B mutations and copy number alterations from the cBioPortal website (33) showed a very low alteration rate (0.2%) in breast cancer (Supplementary Fig. S5). Then, we analyzed the mRNA expression level of EFA6B/PSD4 in our collection of 43 human mammary cell lines previously profiled using DNA microarrays and of which 14 are characterized as Cld-low. EFA6B/PSD4 expression was lower in Cld-low cell lines than in non–Cld-low cell lines (P < 0.001, Student t test; Fig. 6A). Subsequently, we measured the expression levels of EFA6B/PSD4 in our large dataset of 5,252 invasive breast carcinomas (Supplementary Table S1). From this cohort, a total of 500 tumors showed a 2-fold or greater upregulation and 281 showed a 2-fold or greater downregulation of EFA6B/PSD4, using normal breast tissue as the standard. Correlation analysis between EFA6B/PSD4 expression and histoclinical features of tumors (Table 1) showed that EFA6B/PSD4 downregulation was associated (Fisher exact test) with higher grade (P < 0.001), higher frequency of ER− status (P < 0.001), PR− status (P < 0.001), triple-negative status (P < 0.001), and more frequently basal (P < 0.001) and Cld-low subtypes (P < 0.001). Correlation also existed with disease-free survival (DFS) within the 2,930 non–stage IV patients with follow-up available. A total of 1,056 patients experienced a disease event, and 1,874 remained disease-free after a median follow-up of 82 months. The 5-year metastasis-free survival (MFS) was 67% [95% confidence interval (CI), 65%–69%] for the whole population, and 54%, 67%, and 76% in cases of downregulation, no deregulation, and upregulation, respectively (Fig. 6B). Importantly, for all of the histoclinical correlations described above, a continuum existed from downregulated to neutral to upregulated levels of EFA6B/PSD4 (Table 1).
In agreement with the survival data, we found that tumors with EFA6B/PSD4 downregulation were more frequently classified “high risk” of relapse by two major prognostic gene expression signatures, the Amsterdam 70-gene signature and the Oncotype Recurrence Score. Here too, a continuum was present between the three groups defined by EFA6B/PSD4 expression (downregulation, neutral, and upregulation) and the above gene signatures (Table 1).
We next analyzed correlations between expression levels of EFA6B/PSD4 and those of 20 genes involved in tight junction or EMT. As shown by the correlation matrix (Fig. 6C), EFA6B/PSD4 expression positively correlated with that of genes involved in tight junction, including claudin genes, and negatively correlated with that of genes involved in EMT. Of note, the three other EFA6/PSD genes did not correlate with the “tight junction” and “EMT” gene clusters. We also found correlations between EFA6B/PSD4 expression and gene expression signatures related to epithelial and mesenchymal phenotypes (Table 1). Tumors with EFA6B/PSD4 downregulation had a higher expression score of the Stroma signature, the Stromal differentiation score, and the EMT score than tumors without no deregulation or tumors with upregulation. They also showed a differentiation score closer to mammary stem cells than mature luminal cells. Together with the low expression of tight junction genes and high expression of EMT genes, these findings are consistent with an EMT phenotype in the EFA6B/PSD4–downregulated tumors.
EFA6B and claudin-3 immunostaining expression in breast cancer
To extend the transcriptome analyses, we carried out IHC studies on human mammary tumor samples and assessed the protein levels of EFA6B in relation to the subtype and prognosis. We analyzed a collection 412 human primary invasive breast adenocarcinoma spotted onto a TMA. The Cld-low subtype has only been defined by gene profiling; nevertheless, claudin-3 appears to be one of the most significantly affected genes in this subset (2, 34), so we quantified the levels of claudin-3 protein by IHC as a marker for the Cld-low subtype. We found a positive correlation between the presence of a triple-negative IHC for ER, PR, and HER2 and low expression of both EFA6B and claudin-3 proteins in primary breast tumors (Fig. 6D). Moreover, in tumors with a low expression of EFA6B, a low claudin-3 expression was concomitantly noted (Fig. 6D). Finally, patients with low levels of EFA6B expression correlated significantly with a lower DFS expressed as the cumulative incidence of relapse (Fig. 6E). Together with our transcriptomic analyses, these data showed that EFA6B is lost in aggressive mammary tumors downregulated in tight junction.
Discussion
We have provided evidence that the tight junction regulator EFA6B acts as an antagonist to human breast cancer development. In vitro, we have identified two biologic processes regulated by EFA6B through which it might hamper or even revert breast cancer progression: EFA6B inhibits the TGFβ-induced EMT, and EFA6B increases the formation and fusion of VACs to promote epithelial apico-basal polarity of mammary acinis. In vivo, breast cancer analyses show that loss of expression of EFA6B is a marker of poor prognosis and associated with the TNBC and Cld-low subtypes.
On the basis of these results and previous observations from us and others, we propose that PSD4/EFA6B is a new potential tumor suppressor and that an early loss of expression of EFA6B contributes to the malignant progression of breast cancer and might participate in the orientation toward the TNBC and Cld-low subtypes. Three lines of evidence corroborate this model: (i) in vitro manipulation of EFA6B levels affects the tight junction, the apico-basal polarity, and the EMT status of various tumoral and nontumoral human mammary cell lines: all properties that are shown to correlate, in vivo, with tumor progression and found to be necessary for tumor cell migration and invasion (35–38); (ii) likewise, in vivo, the gradual loss of expression of EFA6B correlates with a progressive increase in the appearance of breast cancer tumors characterized by a downregulation of the tight junction along with the emergence of an EMT signature; (iii) finally, loss of expression of EFA6B (mRNA and protein) in human mammary tumors is correlated with the formation of tumors of the Cld-low subtype primarily defined by a loss of tight junction and an upregulation of the EMT signature.
To further substantiate the antitumoral role of EFA6B, we set out to decipher the molecular mechanism(s) by which it reverted the MCF7 cell aggregates into normal-like mammary acinis in 3D culture. Our results demonstrate that EFA6B regulates both the formation of VACs, and lumens. The EFA6–Arf6 complex is known to control cell surface transport, therefore it is possible that EFA6B overexpression stimulates the formation of endocytic vesicles from the basal surface of the plasma membrane to form VACs, onto which EFA6B is localized, and then helps their subsequent fusion at cell–cell contacts to generate the lumens. Because the newly formed lumens are delineated by tight junctions, EF6AB might also help form the lumen by strengthening its surrounding tight junction, which have independently been shown to promote lumen formation (39).
In addition, the levels of expression of EFA6B determined the epithelial–mesenchymal morphologic status of MCF7 cell aggregates, which indicated that EFA6B might control EMT. This is of primary importance because EMT has been largely implicated in carcinogenesis. TGFβ is a potent EMT-inducer of human mammary cancer (21, 37, 40). One of the earliest consequences of TGFβ binding to its receptor is the dissolution of the tight junction (26). Interestingly, exogenous expression of EFA6B blocked the TGFβ-induced EMT in MCF7 and MDCK cells. Preliminary examination of the major downstream effectors of TGFβ implicated in tight junction remodeling did not allow us to identify the signaling pathway(s) inhibited by EFA6B. In the past, we have shown that EFA6B and Arf6 control the tight junction by regulating the apical ring of acto-myosin (4, 41). Thus, the complex EFA6B–Arf6 may act directly on the apical tight junction actin ring, possibly counteracting or bypassing the effects of RhoA or other actin regulators.
Collectively, our results indicate that downregulation of EFA6B may contribute to breast cancer malignant progression. We tested this hypothesis by analyzing the expression of EFA6B in human breast tumor samples. Not only did we observe a reduced expression of EFA6B in human breast cancer, but in agreement with our in vitro studies, it was correlated with a decrease in tight junction components. This was associated with an increase in the signatures for EMT and CSC, and a lower DFS. The transcriptomic analysis showed that tumors with low EFA6B expression were more frequently of the Cld-low subtype. This breast cancer subtype is also defined by transcriptomic signatures associated with EMT and CSC, and a poor prognosis. Thus, collectively our transcriptomic, IHC, and in vitro data support a model whereby the loss of EFA6B might participate in the initiation or maintenance of the tight junction–downregulated Cld-low subtype. In any case, independent of a possible role of EFA6B on the subtype orientation, it is striking that the same biologic changes associated with EFA6B downregulation in human breast cancer and the ones induced by a decrease in EFA6B levels in cell culture—namely loss of tight junction and apico-basal polarity, increased EMT, stromal and CSC signatures—are all found in breast cancer tumoral progression.
In summary, EFA6B appears as a potential tumor antagonist whose downregulation may trigger the progression toward the Cld-low subtype. Inhibition of the mechanisms that downregulate EFA6B levels may offer a much needed therapeutic option.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: F. Luton
Development of methodology: J. Zangari, F. Bertucci, O. Cabaud, M. Lopez, P. Hofman, F. Luton
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J. Zangari, M. Partisani, F. Bertucci, C. Berruyer-Pouyet, P. Finetti, E. Long, O. Cabaud, D. Birnbaum, M. Lopez, P. Hofman
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J. Zangari, F. Bertucci, G. Bidaut, P. Finetti, F. Brau, B. Chetaille, D. Birnbaum, M. Lopez, P. Hofman, M. Franco, F. Luton
Writing, review, and/or revision of the manuscript: F. Bertucci, P. Finetti, F. Luton
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J. Zangari, J. Milanini, M. Franco
Study supervision: F. Luton
Acknowledgments
The authors thank Drs. E. Van Obberghen-Schilling and J. Mazella for the generous gift of various reagents.
Grant Support
This work was supported by the Centre National de la Rechereche de Scientifique (CNRS) and the National Research Agency through the “Investments for the Future” LABEX SIGNALIFE (ANR-11-LABX-0028-01).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.