Functional proteomic analysis reveals the involvement of KIAA1199 in breast cancer growth, motility and invasiveness

Background KIAA1199 is a recently identified novel gene that is up-regulated in human cancer with poor survival. Our proteomic study on signaling polarity in chemotactic cells revealed KIAA1199 as a novel protein target that may be involved in cellular chemotaxis and motility. In the present study, we examined the functional significance of KIAA1199 expression in breast cancer growth, motility and invasiveness. Methods We validated the previous microarray observation by tissue microarray immunohistochemistry using a TMA slide containing 12 breast tumor tissue cores and 12 corresponding normal tissues. We performed the shRNA-mediated knockdown of KIAA1199 in MDA-MB-231 and HS578T cells to study the role of this protein in cell proliferation, migration and apoptosis in vitro. We studied the effects of KIAA1199 knockdown in vivo in two groups of mice (n = 5). We carried out the SILAC LC-MS/MS based proteomic studies on the involvement of KIAA1199 in breast cancer. Results KIAA1199 mRNA and protein was significantly overexpressed in breast tumor specimens and cell lines as compared with non-neoplastic breast tissues from large-scale microarray and studies of breast cancer cell lines and tumors. To gain deeper insights into the novel role of KIAA1199 in breast cancer, we modulated KIAA1199 expression using shRNA-mediated knockdown in two breast cancer cell lines (MDA-MB-231 and HS578T), expressing higher levels of KIAA1199. The KIAA1199 knockdown cells showed reduced motility and cell proliferation in vitro. Moreover, when the knockdown cells were injected into the mammary fat pads of female athymic nude mice, there was a significant decrease in tumor incidence and growth. In addition, quantitative proteomic analysis revealed that knockdown of KIAA1199 in breast cancer (MDA-MB-231) cells affected a broad range of cellular functions including apoptosis, metabolism and cell motility. Conclusions Our findings indicate that KIAA1199 may play an important role in breast tumor growth and invasiveness, and that it may represent a novel target for biomarker development and a novel therapeutic target for breast cancer.


Background
Breast cancer as the most commonly diagnosed and the second leading cause of cancer-related death in women, is responsible for approximately 40,000 deaths in the United States each year [1]. At the time of diagnosis, a majority of patients have metastases to regional and distant sites, which is a major cause of cancer-related mortality [2]. Chemotaxis, cellular migration driven by chemokine gradients, is a critical process involved in tumor invasion and metastasis in various types of cancers including breast cancer [2]. Cell migration is a highly polarized process characterized by protrusion of a leading pseudopodium at the front and establishment of a trailing rear compartment or tail region at the back. Our earlier, comprehensive proteomic analysis of the pseudopodium and cell body in chemotactic cells provided a rich source of information for investigating key signaling pathways and proteins involved in chemotaxis and cancer metastasis [3]. When we compared our pseudopodium proteome dataset with the breast cancer gene expression dataset [4], a protein without a defined function in breast cancer, KIAA1199, caught our attention, as only identified in pseudopodium and highly up-regulated in aggressive breast cancer tissues and cells.
The KIAA1199 gene which was first discovered to be involved in non-syndromic hearing loss is expressed in a wide range of normal human tissues, with the highest expression level in brain [5]. The KIAA1199 gene is located on 15q25, where a brain tumor suppressor gene has been mapped [6]. It is highly expressed in three basal type B breast cancer cell lines (HS578T, MDA-MB-231, and BT549) and the expression of this gene is significantly correlated with the invasive ductal carcinoma type of breast cancer [7]. Also, the high expression of KIAA1199 in gastric tumors is associated with a poor prognosis and with lymph node metastasis [8]. These findings are consistent with a recent report which showed that repression of KIAA1199 attenuates Wnt-signaling and decreases the proliferation of colon cancer cells [9]. Other studies have shown that up-regulation of the KIAA1199 gene is associated with cellular mortality [10] and that the KIAA1199 expression level is significantly elevated upon p53 activation [11]. Based on these observations, we hypothesized that KIAA1199 is a novel regulator of breast cancer growth and aggressiveness.
In this report, we demonstrated the overexpression of KIAA1199 mRNA and protein in breast tumors and invasive cell lines as compared to non-neoplastic tissue and non-invasive cells. Knockdown of KIAA1199 inhibited cell proliferation and motility in vitro and tumor incidence and growth in vivo. Our comprehensive functional proteomic study to analyze the consequences of KIAA1199 knockdown in the breast cancer cell line MDA-MB-231 demonstrate that KIAA1199 may play an important role in the pathogenesis of breast cancer and that it may represent a novel therapeutic target for breast cancer.
MDA-MB-231 and Hs578T cells (obtained from ATCC (Manassas, VA)) were cultured in DMEM containing 10% FBS, 100 U/ml penicillin and 100 μg/ml streptomycin at 37˚C in an atmosphere containing 5% CO 2 . The SILAC labeling was performed according to the manufacture's protocol. The lysine and arginine depleted DMEM medium supplemented with L-[ 12 C 6 ]arginine and L-[ 12 C 6 ] lysine was used for light condition and the depleted DMEM medium supplemented with L-[ 13 C 6 ]arginine and L-[ 13 C 6 ]lysine was used for heavy condition.

Knockdown of KIAA1199 by shRNA-mediated RNA interference
Four different sets of annealed oligonucleotides specific for the KIAA1199 gene sequence were cloned into the pGPH1/GFP/NEO shRNA expression vector obtained from Genepharma (Shanghai, China). These vector constructs (in addition to an empty vector) were transfected into MDA-MB-231 and Hs578T cells to generate the KIAA1199 knockdown cells (ShA and ShB) and control (ShNC) cells respectively. Since the shRNA plasmids contain the neomycin resistance gene and green fluorescence protein (GFP) expression cassette the transfected cells were selected using 400 μg/ml of G418 (Invitrogen, Carlsbad, MD) and monitored by fluorescent microscopy (Leica, Bannockburn, IL) and flow cytometry.

mRNA analysis
Total-RNA was extracted from 1×10 7 cells (cultured in triplicate) using Trizol reagent (Invitrogen,Carlsbad, CA) according to the manufacturer's instructions. RNA (2-5 μg) was treated with DNAse I (Promega), then reverse transcribed, using 200 U Superscript II (Invitrogen) and 250 ng random primers (Invitrogen), according to the manufacturer's instructions. The resulting cDNA diluted 1:5 in nuclease-free water and stored in aliquots at −80°C until used. The RT-PCR amplification of KIAA1199 was performed with a denaturation step at 95°C for 10 min, followed by 32 cycles of denaturation at 95°C for 1 min, primer annealing at 56°C for 30 s, and primer extension at 72°C for 30 s. The PCR conditions varied for S100A11 (35 cycles, annealing at 60°C for 30 s, and primer extension at 72°C for 45 s), WASL (28 cycles, annealing at 60°C for 45 s, and primer extension at 72°C for 90 s), PPP1R9B (30 cycles, annealing at 60°C for 30 s, and primer extension at 72°C for 60 s) and GAPDH (30 cycles, annealing at 53°C for 30 s, and primer extension at 72°C for 30 s). Upon completion of the cycling steps, a final extension at 72°C for 5 min was done for all of the reactions and then the reactions were stored at 4°C. The bands obtained after electrophoresis were quantified by densitometry and their intensities were normalized to that provided by the GAPDH (Glyceraldehyde 3-phosphate dehydrogenase) band (relative integral optical density (IOD)) as described before [13]. The average intensity value of the transcripts obtained from the negative control cells were set to 100%. A list of primers is provided in Additional file 1: Table S1.

Cell motility and migration assay
Wound healing assay was performed to determine cellular motility as described before [14]. Briefly, cells were separately seeded at a density of 5 × 10 5 cell/well in a 6-well plate (triplicate for knockdown and control cells) and grown to confluence in serum containing DMEM media. The monolayer was scratched using a pipette tip and washed with PBS to remove floating cells and refed with serum containing DMEM media. The wounds were photographed immediately after scratching and again 24 h refeeding. The inhibition in wound closure was qualitatively evaluated.
In order to quantitatively examine the effect of KIAA1199 knockdown in breast cancer cells, we performed trans-well motility assays utilizing 6.5 mm Transwell® with 8.0 μm pore polycarbonate membrane filters (Corning Corp, Corning, NY). Single cell suspensions were seeded onto the upper surface of the filters in supplemental serum free McCoy's 5A medium. The bottom chamber contained 1.0 ml serum containing media. MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) was added and cells were incubated for an additional 3 h. Cells from the top of the transwell chambers were removed using a cotton swab (residual cells). The transwell chambers (migrated cells) and cotton swab containing residual cells were plated in separate well of a 24-well plate containing 400 μl of DMSO. Following 1 h of gentle shaking, 100 μl samples were removed and absorbancy was determined at 570 nm using a microtiter plate reader. The percent migratory activity was calculated as: percent migration = [(A / B) -1 × 100], where A is the number of migrated cells and B is the number of residual cells. Percent migratory activity was compared between different groups. The assay was performed in triplicate.

Cell proliferation and apoptosis assay
MDA-MB-231 and Hs578T stable cell lines were plated at 2 × 10 3 cells/well in 96-well plates (triplicate for knockdown and control cells). Following overnight adherence, cells were incubated with serum containing media for various durations. Cell proliferation was determined by MTT (3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide, a yellow tetrazole) assay. The differences in absorbance were compared in vector control transfected cells and KIAA1199 knockdown cells. To determine the role of KIAA1199 in apoptosis, isogenic variants of MDA-MB-231 and Hs578T stable cell lines were grown in DMEM with 10% FBS. A total of 1×10 6 cells were washed with PBS (phosphate buffered saline), collected and double-stained for Propidium Iiodide (PI) and Annexin V using the Annexin-V-FLUOS Staining Kit (Roche Applied Science, Mannheim, Germany) according to the manufacturer's instructions. The frequency of apoptotic cells was analyzed using the FACSCalibur flow cytometer (BD Biosciences, San Jose, CA) with CellQuest Pro software (BD Biosciences).

Tumor growth assay
Mice were housed and handled according to protocols approved by the University of Nebraska Medical Center Institutional Animal Care and Use Committee. Two groups (n = 5) of female BALB/C nude mice (Charles River, Wilmington, MA), 6-8 weeks of age, housed under pathogen free conditions were used. MDA-MB-231-ShNC and MDA-MB-231-ShB cell monolayers were trypsinized and washed with Hank's balanced salt solution (HBSS) 3 times and counted using trypan blue (Sigma) exclusion dye. Single cell suspensions of 1x10 6 cells (>95% viability) in 100 μL were injected into the mammary fat pad. Twice a week tumor size was measured using digital calipers (Fisher Scientific, Pittsburgh, PA). Tumor volume was calculated according to the formula Volume = W 2 × L/2, where W = short diameter, and L = long diameter. Mice were euthanized and primary tumors were removed and processed by formalin fixation with subsequent embedding in paraffin for immunohistochemistry.

Immunohistochemical analysis
IHC analysis was performed as described previously [15] using the rabbit polyclonal anti-KIAA1199 (Sigma-Aldrich; 1:10 dilution), the rabbit monoclonal anticleaved caspase 3 (CASP3; Invitrogen; 1:500 dilution) and the mouse monoclonal anti-proliferating cell nuclear antigen (PCNA; Santa Cruz, CA; 1:40 dilution) as primary antibodies. Tumor sections were deparaffinized by incubation in EZ-Dewax (BioGenex Laboratories Inc, San Ramon, CA) and rinsed in distilled water to remove residual EZ-Dewax. Following nonspecific blocking for 30 min, sections were incubated with primary antibodies overnight at 4°C. Sections were then washed and subsequently incubated at room temperature with the respective biotinylated secondary antibodies (1:500 in PBS) for 45 min. Immunoreactivity was visualized by incubating the avidin-biotin complex with diaminobenzidine tetrahydrochloride substrate (Vector Laboratories, Burlingame, CA). The sections were observed microscopically (Nikon, Melville, NY) using 5 × 5 reticle grid (Klarmann Rulings, Litchfield, NH) and stained cells and vessels were identified. The slides were lightly counterstained with Harris hematoxylin and viewed under a light microscope.
The breast cancer TMA slide (catalog number A712(12) and A712(13)) was purchased from AccuMax (Seoul, Korea). A human kidney tissue was used as positive control. The slide was processed for IHC detection of KIAA1199 expression with a polyclonal anti-KIAA1199 primary antibody (1:10 dilution; SigmaAldrich). An iSan Coreo slide scanner (Ventana Medical Systems, AR) was used to scan the slide at 40× and the resulting images were analyzed by Metamorph Imaging Software (Molecular Devices, CA) to determine the intensity of immunostaining. Immunostaining index (arbitrary unit) was calculated by considering the level of immunostaining intensity and the area with KIAA1199 positivity.

Quantitative proteomic analysis
MDA-MB-231-ShNC (cultured in light medium) and MDA-MB-231-ShB Cells (cultured in heavy medium) were grown in doublet SILAC conditions and the proteomic samples were prepared as previously described [16]. Briefly, MDA-MB-231-ShNC and MDA-MB-231-ShB cells were seeded at 20-30% confluence and harvested when cell density reached 90%. After 10 passages, heavy (Arg6, Lys6) labeled MDA-MB-231-ShB and MDA-MB-231-ShNC cells (Light) were harvested separately in 7 M urea, 2 M thiourea and 50 mM ammonium bicarbonate. Equal amounts of protein were combined from each condition. Following tryptic digestion and chromatography separation via strong cation exchange (SCX), a total of 21 fractions of peptide mixtures were subjected to C18 reverse-phase liquid chromatography (Eksigent, Dublin, CA) coupled online to an LTQ-Orbitrap mass spectrometer (Thermo Scientific, Bremen, Germany). Two biological replicates were performed. The MS data were analyzed using the UNiquant software pipeline [16]. Briefly, DeconMSn (http://omics.pnl.gov/software/) was used to determine and refine the monoisotopic mass and charge state of parent ions from the LTQ-Orbitrap raw data, and to create a peak list of these ions in .mgf format. The peak list contained the fragment information such as the MS/MS spectra, refined precursor ion and charge state. DtaRefinery (http://omics.pnl.gov/software/) was used to improve mass measurement errors for parent ions of tandem MS/MS data by modeling systematic errors based on putative peptide identifications using the algorithm as described [16]. A script written in Python (programming language) was used to automate the process of generating .mgf files from raw data using DeconMSn and DtaRefinery. The resulting .mgf file was submitted to Mascot (version 2.2, Matrix Science, London, U.K.) database searching against (i) a concatenated database containing 73,928 proteins from international protein index (IPI) database (version 3.52), (ii) the commonly observed 262 contaminants (forward database), and (iii) the reversed sequences of all proteins (reverse database). Carbamidomethylation was set as the fixed modification and oxidation of methionine was set as the variable modification. The initial mass deviation tolerance of precursor ion was set to 10 ppm and fragment ion tolerance was set to 0.5 Da. A maximum of 2 missed cleavages were allowed in peptide identification. Identified peptides were sorted by a descending order of Quality of Peptide Identification (QPI) which is defined by the Mascot peptide identification score (a minimum of 10) divided by the square root of the precursor ion mass error. A cutoff of QPI was applied to ensure a total false discovery rate (FDR) for peptide identification < 0.01 evaluated by reverse database approach [16].

Statistical analysis
In vivo data analysis was performed using the Mann-Whitney U-test for significance. For the quantitative analysis of differentially expressed proteins identified by LC-MS/MS, a mixed-effects model with random effects from the two experimental runs was fit to the log2 of the protein fold changes to test whether the log2 of protein fold change was significantly different from zero. Note that a differentially expressed protein is expected to have a non-zero log2 fold change. The p-value was calculated and further corrected by the Benjamini-Hochberg (BH) procedure [17] to control the false discovery rate to be no more than 0.05. A protein with a BH corrected p-value equal-to-or-less-than 0.05 was considered to be statistically significant. For the TMA analysis immunostaining index was tested using the paired t-test to determine the significance of difference between the carcinoma and non-neoplastic cores. The TMA results were reviewed by three independent pathologists.

Ethics statement
All procedures performed in vivo tumor growth and metastasis studies were in accordance with institutional guidelines and approved by the University of Nebraska Medical Center Institutional Animal Care and Use Committee.

Expression of KIAA1199 in breast cancer specimens
In order to assess the clinical relevance of KIAA1199 in breast cancer we performed a bioinformatics study of a large database of microarray data from cancer experiments available at the Oncomine website (www.oncomine.org). We observed the overexpression of KIAA1199 mRNA in breast tumor tissues (see Discussion) as compared to nonneoplastic tissue (Table 1). We performed a tissue microarray (TMA) analysis to examine the KIAA1199 protein expression level in breast carcinoma and normal tissues ( Table 2). As shown in the Additional file 2: Figure S1 a human kidney tissue was used as positive (cells in tubules) and negative (cells in glomeruli) control for immunohistochemical staining (according to the human protein atlas at http://www.proteinatlas.org KIAA1199 has the highest expression level in renal tubules). Figure 1 illustrates the cytosolic localization of KIAA1199 and results of immunohistochemical staining of a TMA slide containing 12 breast tumor tissue cores (rows a, c and e) and 12 corresponding normal tissues (rows b, d and f). We quantified and evaluated the KIAA1199 protein expression by analyzing the intensity of immunostatining and positive areas percentage in each core image using the Metamorph software (Zeiss). We observed a 14.66 fold overexpression of KIAA1199 protein in breast tumor tissues (t-test, p = 0.025) compared to non-neoplastic breast tissues (Figure 1).

Knockdown of KIAA1199 in breast cancer cell lines
The construction of the silencing vector pGPH1/GFP/ NEO is shown in Additional file 3: Figure S2. Two different sets of annealed oligonucleotides (ShA and ShB) were used to knockdown the KIAA1199 gene in both MDA-MB-231 and Hs578T cells. We evaluated the efficiency of knockdown through both RT-PCR and Western blotting approaches in triplicate. As shown in the Additional file 3: Figure S2, we observed an average of 86% and 92% decrease in the level of KIAA1199 transcription in  Next, we examined whether KIAA1199 knockdown modulated breast cancer cell proliferation. KIAA1199 knockdown in both MDA-MB-231 and Hs578T cells (the experiment was performed in three biological replicates) significantly inhibited the cell proliferation ( Figure 2C) as compared to the vector control transfected cells (t-test, P < 0.05).
In order to study the effect of KIAA1199 knockdown on apoptosis, we performed flow cytometric analysis using AnnexinV + (early apoptosis marker) and AnnexinV + /PI + (late apoptosis) cells. We observed higher frequency of cells programmed for both early and late phases of apoptosis in KIAA1199 knockdown cells as compared to vector controls ( Figure 3A). We observed an average of 1.72 and  monoclonal antibody (Cell Signaling) which detects both pro-caspase-3 and cleaved caspase-3. As shown in Figure 3B, we observed an overrepresentation of cleaved caspase-3 in KIAA1199 knockdown cells compared to control cells. Together these data suggest that KIAA1199 knockdown inhibited cellular migration and proliferation and enhanced apoptosis. Since the MDA-MB-231-ShB seemed to be more efficiently affected during the KIAA1199 we choose to use this cell line together with MDA-MB-231-ShNC for further in vivo studies and proteomic analyses.

KIAA1199 knockdown inhibits tumor incidence/growth and cell proliferation
To determine whether KIAA1199 depletion modulates tumor growth, we implanted the MDA-MB-231-ShNC (control) and MDA-MB-231-ShB cells into the mammary fat pads of nude mice (n = 5). We observed significant reduction in tumor incidence following KIAA1199 knockdown ( Figure 4A)

Hs578T-ShNC
Hs578T-ShA Hs578T-ShB software pipeline [16]. Although the SILAC based proteomic study was limited to MDA-MB-231-ShNC and MDA-MB-231-ShB cells, the experiment was performed in two independent biological replicates to increase confidence. Total numbers of 1217 and 1404 proteins were identified in replicate 1 and 2 respectively. Among them, 91 proteins were differentially expressed in both replicate experiments (p < 0.05). Using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Uniprot Database, the differentially expressed proteins were classified into eight major categories based on their biological roles and their Gene Ontology (GO) (see Figure 6A). Figure 6B shows the results for representative peptides associated with three of the differentially expressed proteins. Our SILAC-based LC-MS/MS study showed the average upregulation to 1.85 fold for protein S100A11, downregulation to 0.10 fold for WASL and 0.25 fold for PPP1R9B. In order to validate the protein level alteration we performed the semi-quantitative RT-PCR as a standard method to evaluate the transcription level of these proteins [13]. We observed the average of 1.75 and 2.1 fold over-expression of S100A11 in mRNA level in MDA-MB-231-ShA and MDA-MB-231-ShB respectively. Also the transcription level of WASL/PPP1R9B was decreased to 0.14/0.38 and 0.46/0.43 fold in MDA-MB-231-ShA and MDA-MB-231-ShB respectively. These findings showed the accuracy of normalization method used by the UNiquant software pipeline and validated the mass spectral observations (Figure 6C and D). Further data describing the protein changes are detailed in Additional file 4: Table S2 and summarized in Table 3. As shown in Figure 6A, the functions of proteins differentially expressed between MDA-MB-231-ShNC and MDA-MB-231-ShB cells can be assigned to eight categories including Apoptosis, DNA repair and cell cycle, Gene expression and regulation, Cytoskeleton, cell adhesion and motility, Ubiquitin proteasome pathway, Metabolism, Oxidative stress and other proteins. This data suggest that KIAA1199 may affect a broad range of cellular functions.

Discussion
In order to identify new biomarkers for the improvement of new diagnosis strategies and targeted therapy, it is essential to better understand breast cancer biology and the molecular profiles that will respond to targeted treatment. Molecular markers such as progesterone receptor, estrogen receptor, and ErbB2 have been associated with the five major subtypes of breast cancer: luminal A, luminal B, ErbB2+/ER-, basal-like, and normal breast-like [18]. However, molecular pathways involved in incidence, progression and clinical outcomes remain elusive.
Several microarray based expression studies have previously shown the overexpression of KIAA1199 in breast cancer ( Table 1). The results of a recent study from The Cancer Genome Atlas (TCGA) on 593 samples shows 9.094 fold (p = 3.39E-28) overexpression in invasive breast carcinoma, 8.233 fold (p = 1.71E-36) in invasive ductal breast carcinoma and 5.527 fold (p = 7.32E-12) in invasive lobular breast carcinoma compared to corresponding normal breast tissues. Another comparison between invasive breast carcinoma and normal tissue in 158 samples by Gluck and co-workers showed a 2.926 fold (p = 2.48E-7) overexpression of KIAA1199 in invasive breast carcinoma [19]. Furthermore, Richardson and co-workers have reported a 4.125 (p = 1.06E-6) fold overexpression of KIAA1199 in ductal breast carcinoma [20]. In addition to these data, our immunohistochemical study on clinical breast cancer specimens showed 14.66 fold (p = 0.025) overexpression of this protein.
Based on these findings, we examined the role of KIAA1199 in the MDA-MB-231 and Hs578T breast cancer cell lines using two sets of shRNA-mediated knockdown cells for each cell line. We observed that knockdown of KIAA1199 enhanced apoptosis and inhibited cell proliferation and survival in both cell lines in vitro. Additionally, using immunohistochemical staining against cleaved caspase-3 (CASP3) and PCNA we respectively confirmed the apoptosis enhancement and inhibition of cell proliferation in vivo. Interestingly, our proteomic study showed that while the negative control cells expressed higher levels of the apoptosis inhibitors, several proteins involved in apoptosis were overrepresented in the knockdown cells justifying the higher apoptotic activity we observed in vitro and in vivo. For instance the apoptosis regulator BAX which promotes programmed cell death after binding to, and antagonizing the apoptosis repressor BCL2 is upregulated. BAX also accelerates the activation of CASP3, and thereby promotes apoptosis [21][22][23][24]. In addition, we observed the up-regulation of FADD (FAS-Associated Death Domain protein) which is another apoptotic adaptor molecule. FADD bridges the death receptors (e.g. Fas-receptor) to the death-inducing signaling complex (DISC) and activates caspase-8. Active caspase-8 initiates a cascade of caspases which mediate apoptosis [25]. Another example is a large increase in the expression of DIO-1 (death inducer-obliterator-1) that translocates to the nucleus and activates apoptosis in cell culture [26]. KIAA1199 knockdown also led to up-regulation of Akinase anchor protein 8 (AKAP95) that is a potential carrier protein for active caspase 3, carrying it from the cytoplasm into the nuclei in apoptotic cells and is involved in the process of apoptotic nuclear morphological change [27]. It is noteworthy that we found progesterone (P4) receptor membrane component-1 (PGRMC1) down-regulated upon KIAA1199 knockdown. This protein promotes cell survival in human cancer after chemotherapy [28]. PGRMC1 was reported to be over-expressed in breast tumors and other cancer cell lines [29].
It is known that high expression of BAX is associated with improved chemotherapy responsiveness [30] whereas PGRMC1 has a negative impact on chemotherapy by promoting the survival of treated cancer cells [28]. This knowledge plus the fact that KIAA1199 knockdown alters the expression level of these proteins, suggests that KIAA1199 depletion may potentially improve cellular response to chemotherapy.
Our wound healing and transwell cell motility assays showed lower motility in the MDA-MB-231-ShB cells. These findings can be explained by the observation of altered levels of proteins involved in cellular shape change, filopodia extension, nuclear migration and adhesion inhibition in the knockdown cells. We observed the upregulation of S100A11 protein which functions in tubulin polymerization, motility, and tumor invasion [31] and down-regulation of the transforming acidic coiled-coilcontaining protein 3 (TACC3). The latter plays a role in the microtubule-dependent coupling of the nucleus and the centrosome, and it has been demonstrated to be overexpressed in various cancer cell lines [32]. Furthermore, Figure 6 Representative data from the proteomics study on MDA-MB-231 cells after knockdown of KIAA1199. A) Differentially expressed proteins were classified based on their Gene Ontology (GO) and function. B) MS and MS2 spectra from one up-regulated protein (S100A11) and two down-regulated proteins (WASL and PPP1R9B). The left-hand panels show MS spectra, where the red *symbols show intensities of the monoisotopic peak for the light and heavy SILAC labeled peptides. The right-hand panels show the MS2 spectra corresponding to the peptide with the most intense signal of the pair. C) Validation of differential expression of S100A11, WASL and PPP1R9B genes using semi-quantitative RT-PCR analysis in triplicate. D) The relative expression of the bands is shown as bars (*p < 0.01). Intensities were normalized to GAPDH band. The average intensity values of the transcripts obtained from the MDA-MB-231-ShNC cells were set to 100%.
TACC3 depletion has been reported to strongly sensitize cells to chemotherapy [33], therefore KIAA1199 depletion can also potentially affect the cellular response to chemotherapy via TACC3.
Neural Wiskott-Aldrich syndrome protein (WASL) is dramatically down-regulated (0.10 fold) in the KIAA1199 knockdown cells. WASL activates the Arp2/3 complex required for the extension of lamellipodia and filopodia during cell movement [34]. Another down-regulated protein is Neurabin-2 (PPP1R9B) which binds along the sides of F-actin and plays a role in linking the actin cytoskeleton to the plasma membrane at the synaptic junction. PPP1R9B therefore might be involved in cell shape change and migration [35]. A member of the tenascin protein family, the glycoprotein tenascin X (TNXB) is also dramatically down-regulated in the KIAA1199 knockdown cells. As opposed to fibronectin which is adhesive, the tenascins have anti-adhesive effects. TNXB mediates interactions between cells and the extracellular matrix and may support the growth of epithelial tumors [36]. Overall, these findings suggest that KIAA1199 may be involved in determination of cellular morphology and motility.
However, unlike in MDA-MB-231-ShB cells the cell motility was not affected in Hs578Tcell after KIAA1199 knockdown. Although both of these cell lines belong to basal type B breast cancer, MDA-MB-231 cells was originated from invasive ductal carcinomas (IDC) whilst Hs578TT cells originated from a breast carcinosarcoma, and they highly differ in migration and invasion capability [37]. These data suggest discrete cell migratory mechanisms in these cell lines in which KIAA1199 may or may not participate.
In this work we studied the effects of KIAA1199 knockdown for the first time in vivo. We demonstrated the inhibition in tumor incidence and growth rate. Our findings are in concordance with the results of the proteomic study where we observed modulation of several proteins involved in cell cycle progression and division such as ANAPC10 (Anaphase-promoting complex subunit 10), PPP1CB (Serine/threonine-protein phosphatase PP1beta protein catalytic subunit) and PPP2R1A (Serine/ threonine-protein phosphatase 2A regulatory subunit) upon KIAA1199 knockdown. All of these proteins play role in cell cycle regulation and cell division. For example ANAPC10 participates in the progression through mitosis and the G1 phase of the cell cycle [38]. PPP1CB is a component of the PTW/PP1 phosphatase complex, which plays a role in the control of chromatin structure and cell cycle progression during the transition from mitosis into interphase [39] and PPP2R1A is required for proper chromosome segregation and for centromeric localization in mitosis [40]. These data suggest an important role for KIAA1199 in breast cancer incidence, growth and progression.
Mass spectrometry based proteomics holds special promise to provide better insights into biological pathways. In this study, we pursued the functional analysis of KIAA1199 in breast cancer cells as a novel target screened in our previous proteomic study [3]. Although the detailed mechanism of KIAA1199-mediated cellular responses is still obscure, our proteomic study shed light on how different biological pathways may be influenced by KIAA1199 directly or indirectly. For instance alteration of components of MAPK, NF-k-B and apoptosis pathways can potentially affect other cellular phenomena such as angiogenesis.
Furthermore, our findings suggest that KIAA1199 knockdown may also affect the cellular metabolism. It is known that tumor cells typically have much higher rates of glycolysis compared to their normal tissues of origin; consequently they secrete glucose-derived carbon mostly as lactate instead of completely oxidizing glucose. This phenomenon is known as the Warburg effect [41,42]. In this study we observed the modulation of several metabolism associated enzymes. The KIAA1199 knockdown cells have lower expression of proteins involved in glycolysis and cytosolic break down of glucose (such as PGAM4) and instead tend to the mitochondrial oxidation. Therefore, the Warburg effect which is a fundamental character of cancer cells also seems to be negatively influenced by KIAA1199 depletion.
We utilized various approaches and techniques to comprehensively evaluate the major consequences of KIAA1199 depletion in breast cancer cells in vitro and in vivo. Despite the limitations in study sizes (the number of TMA tissues, animals, cell lines etc.) we studied several aspects of cancer development and progression following KIAA1199 knockdown. Further studies on each aspect with larger sample sizes will help to uncover the mechanism of KIAA1199 function and provide more evidences. Taken together, our findings presented here suggest that KIAA1199 may represent a novel target for biomarker development and a novel therapeutic target to control breast cancer progression and metastasis.