Pharmaceutical Sciences, Faculty of
Non UBC
DSpace
BMC Molecular and Cell Biology. 2020 Apr 03;21(1):24
The Author(s)
Thejer, Bashar M
Adhikary, Partho P
Kaur, Amandeep
Teakel, Sarah L
Van Oosterum, Ashleigh
Seth, Ishith
Pajic, Marina
Hannan, Katherine M
Pavy, Megan
Poh, Perlita
Jazayeri, Jalal A
Zaw, Thiri
Pascovici, Dana
Ludescher, Marina
Pawlak, Michael
Cassano, Juan C
Turnbull, Lynne
Jazayeri, Mitra
James, Alexander C
Coorey, Craig P
Roberts, Tara L
Kinder, Simon J
Hannan, Ross D
Patrick, Ellis
Molloy, Mark P
New, Elizabeth J
Fehm, Tanja N
Neubauer, Hans
Goldys, Ewa M
Weston, Leslie A
Cahill, Michael A
Background Progesterone Receptor Membrane Component 1 (PGRMC1) is expressed in many cancer cells, where it is associated with detrimental patient outcomes. It contains phosphorylated tyrosines which evolutionarily preceded deuterostome gastrulation and tissue differentiation mechanisms. Results We demonstrate that manipulating PGRMC1 phosphorylation status in MIA PaCa-2 (MP) cells imposes broad pleiotropic effects. Relative to parental cells over-expressing hemagglutinin-tagged wild-type (WT) PGRMC1-HA, cells expressing a PGRMC1-HA-S57A/S181A double mutant (DM) exhibited reduced levels of proteins involved in energy metabolism and mitochondrial function, and altered glucose metabolism suggesting modulation of the Warburg effect. This was associated with increased PI3K/AKT activity, altered cell shape, actin cytoskeleton, motility, and mitochondrial properties. An S57A/Y180F/S181A triple mutant (TM) indicated the involvement of Y180 in PI3K/AKT activation. Mutation of Y180F strongly attenuated subcutaneous xenograft tumor growth in NOD-SCID gamma mice. Elsewhere we demonstrate altered metabolism, mutation incidence, and epigenetic status in these cells. Conclusions Altogether, these results indicate that mutational manipulation of PGRMC1 phosphorylation status exerts broad pleiotropic effects relevant to cancer and other cell biology.
RESEARCH ARTICLE Open AccessPGRMC1 phosphorylation affects cell shape,motility, glycolysis, mitochondrial form andfunction, and tumor growthBashar M. Thejer1,2†, Partho P. Adhikary1,3†, Amandeep Kaur4,5, Sarah L. Teakel1, Ashleigh Van Oosterum6,Ishith Seth1, Marina Pajic7,8, Katherine M. Hannan9,10, Megan Pavy10, Perlita Poh10, Jalal A. Jazayeri1, Thiri Zaw11,Dana Pascovici11, Marina Ludescher12, Michael Pawlak13, Juan C. Cassano14, Lynne Turnbull15,16, Mitra Jazayeri17,Alexander C. James18,19,20, Craig P. Coorey18,21, Tara L. Roberts18,19,20,21, Simon J. Kinder22, Ross D. Hannan9,10,23,24,25,Ellis Patrick26, Mark P. Molloy11,27, Elizabeth J. New4, Tanja N. Fehm12, Hans Neubauer12, Ewa M. Goldys28,29,Leslie A. Weston30,31 and Michael A. Cahill1,10*AbstractBackground: Progesterone Receptor Membrane Component 1 (PGRMC1) is expressed in many cancer cells, whereit is associated with detrimental patient outcomes. It contains phosphorylated tyrosines which evolutionarilypreceded deuterostome gastrulation and tissue differentiation mechanisms.Results: We demonstrate that manipulating PGRMC1 phosphorylation status in MIA PaCa-2 (MP) cells imposesbroad pleiotropic effects. Relative to parental cells over-expressing hemagglutinin-tagged wild-type (WT) PGRMC1-HA, cells expressing a PGRMC1-HA-S57A/S181A double mutant (DM) exhibited reduced levels of proteins involvedin energy metabolism and mitochondrial function, and altered glucose metabolism suggesting modulation of theWarburg effect. This was associated with increased PI3K/AKT activity, altered cell shape, actin cytoskeleton, motility,and mitochondrial properties. An S57A/Y180F/S181A triple mutant (TM) indicated the involvement of Y180 in PI3K/AKT activation. Mutation of Y180F strongly attenuated subcutaneous xenograft tumor growth in NOD-SCID gammamice. Elsewhere we demonstrate altered metabolism, mutation incidence, and epigenetic status in these cells.Conclusions: Altogether, these results indicate that mutational manipulation of PGRMC1 phosphorylation statusexerts broad pleiotropic effects relevant to cancer and other cell biology.Keywords: Mitochondria, Migration, Invasion, Metabolism, Cytochrome P450, Mesenchymal amoeboid transition,Proteomics, Tumor biology© The Author(s). 2020 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License,which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you giveappropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate ifchanges were made. The images or other third party material in this article are included in the article's Creative Commonslicence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commonslicence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtainpermission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to thedata made available in this article, unless otherwise stated in a credit line to the data.* Correspondence: mcahill@csu.edu.au†Bashar M. Thejer and Partho P. Adhikary are Co-first author.1School of Biomedical Sciences, Charles Sturt University, Wagga Wagga, NSW2650, Australia10ACRF Department of Cancer Biology and Therapeutics, The John CurtinSchool of Medical Research, The Australian National University, ACT, Canberra2601, AustraliaFull list of author information is available at the end of the articleBMC Molecular andCell BiologyThejer et al. BMC Molecular and Cell Biology (2020) 21:24 https://doi.org/10.1186/s12860-020-00256-3Summary statementPGRMC1 phosphorylation site mutations cause pleio-tropic plasticity-related changes including altered mito-chondrial form and function, PI3K/AKT activity,migration, and glucose consumption. Efficient mousexenograft tumor formation required Y180.BackgroundProgesterone (P4) Receptor Membrane Component 1(PGRMC1) is a cytochrome b5-related heme-bindingprotein with multiple functions including interactionwith cytochrome P450 enzymes. Briefly, (the reader isreferred to previous reviews) non-comprehensivePGRMC1 functions include membrane trafficking, P4responsiveness and steroidogenesis, fertility, lipidtransport, neural axon migration, synaptic function, andanti-apoptosis. Its subcellular localization can be cyto-plasmic, nuclear/nucleolar, mitochondrial, endoplasmicreticulum, cytoplasmic vesicles, or extracellular [1–3]. Itis involved in cell cycle processes at the G1 checkpointand during mitosis [4–9], and elevated PGRMC1 expres-sion has been associated with poor prognosis in multipletypes of cancer [2, 10–15].Predicted binding site motifs for Src homology 2(SH2) and Src homology 3 (SH3) proteins in PGRMC1can potentially be negatively regulated by phosphoryl-ation at adjacent casein kinase 2 (CK2) consensus sites[2, 16, 17]. However, while CK2 knockdown leads to re-duced phosphorylation of the corresponding C-terminalCK2 site of PGRMC2, PGRMC1 phosphorylation atS181 was unaffected by CK2 knockout in C2C12 mousemyoblast cells [18]. These and Y180 can all be phos-phorylated in vivo, and constitute a potential regulatedsignaling module [19].We hypothesized that PGRMC1 is a signal hub proteinwith wide ranging effects on cell biology [2, 19–21]. Thehighly conserved motif at Y180/S181 phosphorylationarose early in animal evolution concurrently with theembryological organizer of gastrulation (e.g. Spemann-Mangold organizer), and prior to the evolution of deute-rostomes [20, 22]. Since all of the descendants of theorganism that first acquired S181 exhibit bilateral bodysymmetry, this mutation predated the rules governing cellgrowth and division leading to bilateral body plan andmore types of tissues. The mechanisms that generate andmaintain cell- and tissue-type differentiation states arevery often perturbed in cancer [23]. Therefore, PGRMC1tyrosine phosphorylation could potentially strongly in-fluence mammalian cell function and differentiation status.This present study was prompted by our discovery ofdifferential PGRMC1 phosphorylation status between es-trogen receptor-positive and -negative breast cancers[24]. PGRMC1 was induced in the hypoxic zone ofductal carcinoma in situ breast lesions at precisely thetime and place that cells require a switch to glycolyticmetabolism known as the Warburg effect, leading us topredict a Warburg-mediating role for PGRMC1.Furthermore, a PGRMC1 S57A/S181A double CK2 sitemutant (DM, Fig. 1a) enabled the survival of peroxidetreatment [24]. Sabbir [25] recently reported thatPGRMC1 induced a P4-dependent metabolic change re-sembling the Warburg effect in HEK293 cells, which wasassociated with changes in PGRMC1 stability, post-translational modifications, and subcellular locations.PGRMC1 regulation of glucose metabolism is supportedby its implicated mediation of the placental P4-dependent shift from aerobic towards anaerobic glucosemetabolism in gestational diabetes [26], association withthe insulin receptor and glucose transporters [27], andprobable involvement (based upon AG-205 sensitivity)in P4-mediated increase in neuronal glycolysis [28].We previously observed that MIA PaCa-2 pancreaticcancer (MP) cells [29–32] exhibited marked morpho-logical and metabolic changes when the DM protein wasexpressed [33]. MP cells exist in culture as a mixed ad-herent population of elongated “fibroblast-shaped”morphology, a minority population of rounded morph-ology with bleb-like protrusions, and some multicellularclumps, as well as some rounded suspension cells. Theyhave undergone epithelial-mesenchymal transition [34],and can further undergo mesenchymal-amoeboidtransition (MAT), which requires Rho Kinase- (ROCK)-dependent morphological change from “elongated”mesenchymal cells to rounded amoeboid cells [35].Here, we mutagenically examined the effects of alteredPGRMC1 phosphorylation status on MP cells to gaininsights into PGRMC1-dependent signaling in vitro, re-vealing an important role for Y180. Our objectives wereto initially characterize the effects of mutations at thetwo DM putative negative regulatory phosphorylationsites, S57 and S181, on cell biology, as well as mutatingY180 which forms the center of a putative SH2 targetmotif to test the hypothesis that the DM effect requiredY180. Importantly this investigation was extended to anin vivo system, specifically a subcutaneous mousexenograft tumorigenesis model, in order to address therequirement of Y180 for tumor development in a com-plex environment. In a companion publication [36], wedescribe differences in metabolism, genomic mutationrates, and epigenetic genomic CpG methylation levelsassociated with the PGRMC1 phosphorylation mutantsdescribed in this present study.ResultsPGRMC1 phosphorylation status alters cellularmorphologyWe stably transfected MP cells with the hemagglutinin(HA) epitope-tagged PGRMC1-HA plasmids includingThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 2 of 24the wild-type (WT) sequence [37], the S57A/S181A DM[24], or a novel S57A/Y180F/S181A triple mutant (TM),which removed the phosphate acceptor of Y180 [16, 19,20] (Fig. 1a). Three independent stable cell lines fromeach group expressed both 32 kDa 3xHA-tagged exogen-ous and a 24 kDa endogenous PGRMC1 species, whereasonly the 24 kDa species was present in MP cells (Fig. 1b-d). Both species were present at approximately equimo-lar ratios, and an anti-HA antibody detected only the 32kDa species (Fig. 1d). We reason that any consistentdifferences between biological triplicates should be dueto PGRMC1-HA mutations, rather than clonal artifacts.Fig. 1 MIA PaCa-2 pancreatic cancer cells morphology is affected by PGRMC1 phosphorylation status. a PGRMC1-HA proteins constructed for this figure.TMH: Trans-membrane helix. HA: the C-terminal 3x hemaglutinin tag. b Detection of exogenous PGRMC1 expression levels by western blot (upper panel).Equal loading is controlled by quantifying beta actin (lower panel). The results show three totally independent stably transfected cell lines per plasmidfrom (A). Open arrow: Exogenous PGRMC1-HA (Ex.). Shaded arrow: endogenous PGRMC1 (End.). Filled arrow: beta actin. The molecular weight ladder isBio-Rad 1610377 Dual Xtra Standards. c Box plots quantification of replicate gels of (B) with signals normalized to beta actin from the same respectivelanes. n= 4 lanes for MP and n= 6 for WT, DM and TM (replicates of respective lines 1–3 per condition). There were no significant differences (ns) exceptfor the exogenous band in MP (ANOVA, post-hoc Dunnet’s T3). dWestern blot quantification of HA-tagged exogenous PGRMC1, following B but detectingPGRMC1 with anti-HA antibody. The molecular weight ladder is Abcam ab116028 Prestained Protein Ladder. e PGRMC1 mutant protein expression altersMIA PaCa-2 cell morphology. PGRMC1-HA-expressing stable cells (respective lines 1 from B) or MP cells were stained with a FITC-tagged anti-HA antibody(Anti-HA) and imaged by confocal microscopy. DNA was stained with DAPI. Cells were also imaged in differential interference contrast (DIC) microscopymode. The respective left panels show merged images of all 3 channels. f The rounded phenotype of double and triple mutant (E) was reversed toelongated phenotype after 125 μM ROCKI addition, but not by addition of DMSO vehicle controlThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 3 of 24Subsequent experiments were performed using re-spective cell line triplicates 1–3 per PGRMC1-HAcondition. However, note that MP cells have notundergone hygromycin selection for stably-integratedplasmids, therefore cells in the WT/MP comparisondiffer by two variables (selection and PGRMC1-HA).We focus our analysis onto comparisons betweenWT, DM and TM cells.Like MP cells, freshly seeded WT cells exhibitedpredominantly elongated cell morphology with somerounded cells. DM and TM cells exhibited primarilyrounded morphology (Fig. 1e), which was reminiscent ofthe reported MAT of MP cells [35]. After 72 h of culturethe proportion of round cells in DM and TM cultureswas reduced, but still elevated relative to WT or MP(not shown). Transient transfections with the DM andTM plasmids (but not WT) led to similar increasedlevels of cell rounding across the entire populations ofcells by 24 h after transfection (data not shown), indi-cating that the phosphorylation status of exogenousPGRMC1-HA affects cell morphology.PGRMC1-dependent altered morphology requires rhokinaseThe ROCK pathway is required for amoeboid phenotypeand migration and its inhibition reverses MAT in MPcells [35, 38]. ROCK inhibitor (ROCKI) reversed therounded phenotype to elongated for DM and TM (Fig.1f), supporting the hypothesis that morphological transi-tion involves altered actin organization. It remainsunclear whether the process is truly MAT.PGRMC1 phosphorylation affects cell motility andinvasionTo further investigate cell plasticity imposed byPGRMC1-HA phosphorylation mutants, we examinedcell motility via a scratch assay [39]. MP cells exhibitedthe lowest migration, while DM cell migration was sub-stantially greater than other cell lines (Fig. 2a-b, File S1).DM cells migrated predominantly as rounded cells, usingextended filopodia and small pseudopodia, however, a mi-nority of flattened cells exhibited more pronouncedpseudopodia. Video imaging demonstrated that these cellshapes could rapidly interconvert (File S1C). Conversely,DM cells exhibited the lowest ability to invade through aGeltrex pseudo-basal membrane (Fig. 2c-d).PGRMC1 phosphorylation imposes broad changesA total of 1330 proteins were reliably identified by pro-teomics in at least one sample with at least 2 peptidesusing a combination of information dependent acquisi-tion (IDA) and data independent SWATH-MS (Sequen-tial Window Acquisition of all Theoretical MassSpectrometry) acquisition. Results are provided as FileS2. Approximately 50% of variation was explained bytwo principal components (PCs) in PC analysis, whichcorresponded approximately to “ribosomes and trans-lation” (PC1, separated MP&DM from WT&TM) and“mRNA splicing and processing” (PC2, separatedMP&WT from DM&TM) (Fig. S1A, File S3). Of theidentified proteins, 243 differed by 1.5 fold or more be-tween one or more comparisons with p < 0.05 (t-test),and 235 of these withstood PC multiple sample correc-tion. The heat map clustering of those 243 proteins(Fig. 3, File S6) revealed a suite of proteins whichstrongly discriminated between the different PGRMC1-HA-induced conditions. Biological replicates clusteredtightly in clades of the same cell type, with largedistances between clades. We conclude that thesedifferences are primarily specific PGRMC1-HA mutant-dependent effects.Results from those six comparisons of protein abun-dance between the four sample types [1) MP v. WT, 2)MP v. DM, 3) MP v. TM, 4) WT v. DM, 5) WT v. TM,and 6) DM v. TM] provide lists of significantly differen-tially abundant proteins in each pair-wise comparison(File S2). WebGestalt enrichment analyses were per-formed to identify pathways or features either signifi-cantly more or less abundant (respectively the “red” and“blue” lists of proteins for each comparison from col-umns B of File S4) between each of the respective sixcomparisons at the Benjamini–Hochberg adjusted p-value (adjP) < 0.001 level in at least one pairwise com-parison. WebGestalt mapped features (File S5) are plot-ted against the heat map in File S6 (for primaryWebGestalt data see File S5 and File S6). The results areschematically mapped against the 243 protein heat mapin Fig. 3, using all pathways that were detected by Web-Gestalt in any of the inter-sample comparisons at theadjP< 0.001 level, and including all proteins detected inany of those pathways across all 12 comparisons at theadjP< 0.1 level for respective red and blue protein listsfrom File S4.WT PGRMC1-HA protein induced the elevated abun-dance of many proteins involved in energy metabolism,including proteasomal components involved in proteindegradation, and pathways for amino acid, carbohydrate,and fatty acid catabolism (Fig. 3). These proteins wereannotated as both cytoplasmic and mitochondrial. Per-oxisomal and lysosomal proteins were also upregulatedin WT cells. A suite of proteins putatively involved inthe recognition of mRNA by ribosomes, tRNA aminoa-cylation, ribosomal protein translation, and chaperone-mediated protein folding, was generally less abundant inWT and TM than MP and DM cells (Fig. 3, File S6).Many of the changes in fatty acid and glucose metabol-ism enzymes resemble the effects of the insulin/glucagonsystem of metabolic regulation.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 4 of 24Many of the above proteins were inversely regulated inthe comparison of WT and DM cells (Fig. 3). The energymetabolizing suite of proteins and some proteins associ-ated with translation were less abundant in DM, whereasthe T-complex chaperone complex as well as somenuclear exportins and importins were elevated in DMcells (Fig. S1D-E). Some differentially enriched pathwayswere specific for the WT vs. DM comparison. Differentenzymes associated with heme metabolism were bothup- and down-regulated while enzymes associated withglycosaminoglycan metabolism were up-regulated inWT relative to DM (Fig. 3, File S6).Mitochondrial proteins accounted for a large percent-age of the proteins more abundant in WT than DMcells. Intriguingly, many cytoplasmic proteins were moreabundant in WT cells than DM (Fig. 3, File S6). We alsonoted higher abundance components of ATP synthase inWT and TM cells (Fig. S1C), changes of proteinsFig. 2 PGRMC1 phosphorylation status affects motility and invasion. a DM cells exhibit enhanced motility in a scratch assay. Representative resultsfor the scratch assay after 36 h of cell migration. Cells in the boxed area of the scratched void area for each image were counted. b Scratch assaycell migration results. Box plots of cell migration into the scratch void box areas depicted in (A) for multiple replicates. n = 12 (MP), or 4 replicateseach of sublines 1–3 for WT and TM (n = 3 × 4 = 12). The table shows the results of 1 way ANOVA with post-hoc Dunnet’s T3 p-values. Video filesof cell migration are available in File S1. c DM cells exhibit reduced invasion in Geltrex invasion assay. Representative images of crystal-violetstained cells in the lower surface of the transwell insert. d Boxplots of invasion assay results from (C) for replicates, as produced by SPSS software.n = 30 (MP), or 10 replicates each of sublines 1–3 for WT and TM (n = 3 × 10 = 30). The table shows the Kruskal-Wallis p-values for pairwise comparisonsThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 5 of 24involved in chaperonin and microtubule function (Fig.S1D), and a group of proteins involved in major histo-compatibility complex antigen processing and presenta-tion, and proteolysis (Fig. S1F). The latter are reduced inDM cells and may be associated with their reduced inva-siveness (Fig. 2c-d), which requires future confirmation.Consideration of the extreme (highest and lowestabundance) differential proteins for each cell type offersuseful insight into biology (Fig. S2). WT and TM cellsexhibited overlap in the subset of most abundantproteins, which included PSIP1 transcriptional coactiva-tor, TOM40 mitochondrial import channel, as well asCDIPT which catalyzes the biosynthesis of phos-phatidylinositol (circles in Fig. S2). One of the WTabundant proteins was phosphofructokinase (UniProtP08237), which catalyzes the rate limiting reaction andfirst committed step of glycolysis. The most abundantDM proteins included keratin 19, ubiquitin-associatedprotein 2-like, which is involved in stem cell mainten-ance [40], and methyl-CpG-binding protein 2, whichFig. 3 Pathways analysis of SWATH-MS proteomics results. Pathways significantly enriched at the adjP> 0.001 level between all 6 comparisons of“red” and “blue” differential proteins (red = higher abundance, blue = lower abundance, white = equal abundance). Top left: the proteomic heatmap of 243 significantly differential proteins. A color code for WebGestalt pathways is given at top right. Bottom: WebGestalt pathways mapping.This image is derived from File S6, which contains all protein and WebGestalt pathway identitiesThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 6 of 24was more abundant in WT and TM cells, suggestingPGRMC1-mediated changes in genomic methylation(see accompanying paper [36]). The least abundant pro-teins shared a surprising mixed overlap between celltypes (triangles in Fig. S2). TM and WT cells shared lowlevels of Ephrin type-A receptor 2, a tyrosine kinase re-ceptor, which was higher in DM (and MP) cells withoutbeing a top abundant protein in those cells. WT exhi-bited low levels of signal recognition particle 54 kDaprotein, suggesting altered translation of endoplasmicreticulum proteins, and AL1A1 retinal dehydrogenase.This was notable because both DM and TM exhibitedlow levels of AL1A3 NAD-dependent aldehydedehydrogenase involved in the formation of retinoicacid, suggesting alterations in retinoic acid metabolismby mutating S57/S181. DM and TM also shared lowlevels of ApoC3 and ApoA1 (Fig. S2), probably reflectingcommon lower lipoprotein synthesis by those cells.Taken together, our proteomics analysis revealed signifi-cant differences in the abundance of enzymes involvedin diverse cell processes, many of which are directlyimplicated in cancer biology. The resemblance of WTand TM differential proteomics profiles suggests that theDM mutation activates signaling processes that arelargely dependent upon Y180. (The sole difference be-tween DM and TM proteins is the phosphate acceptorFig. 4 RPPA measurements of protein and phosphorylated protein levels. a-d, f-h Average RPPA normalized fluorescent intensity (NFI) from theindicated antibodies (described in methods) is plotted from 6 replicate measurements. NFI is normalized to protein content. Statisticalcalculations for normally distributed data were made using one way ANOVA, and post-hoc Bonferroni (BF) for equal variances (all variances wereequal). For non-parametric data, Kruskall-Wallis (KW) pairwise comparisons were calculated for 24 unrelated samples. * p < 0.05; ** p < 0.01, *** p <0.001. Non-phosphorylated BAD levels could not be accurately determined because signal values were less than three times background. e Theratio of average NFI of D relative to C. Labels follow the aboveThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 7 of 24oxygen of Y180). Overall, this study indicated thatPGRMC1 phosphorylation status exerts higher ordereffects in MP cells. To assess the validity of the pathwaysanalysis, selected enriched pathways were functionallyexamined, as described below.ERR1 activity is not directly affected by PGRMC1Some mitochondrial proteins associated with energymetabolism were predicted by pathways enrichment ana-lysis to be regulated by estrogen receptor related 1(ERR1) transcription factor (Fig. S3A) in the compari-sons of DM cells with both WT (adjP = 0.004) and TM(adjP = 0.04) (File S4). Since ERR1 is a steroid receptor,we investigated any potential link between the biology ofPGRMC1 and ERR1 by attenuating ERR1 levels in WTcells by shRNA. This changed cell morphology frompredominantly elongated to rounded cells (Fig. S3B-D).SWATH-MS proteomics revealed that ERR1 indeedregulated genes differentially abundant between WT andDM cells observed in Fig. S3A and Fig. 3. However,PGRMC1 phosphorylation status affected the abundanceof only a subset of ERR1-driven proteins.PGRMC1 phosphorylation affects PI3K/AKT signalingStrikingly, proteins associated with PI3K/AKT activityfrom Fig. 3 and File S6 were revealed by File S5 to ex-hibit lower abundance in TM cells (Fig. S1B) relative toMP (adjP = 0.0063), WT (adjP = 0.0001) and DM (adjP =0.0002) cells. We assayed the phosphorylation status oftwo AKT substrates by Reverse Phase Protein Array(RPPA). BAD is phosphorylated by both PKA at S112[41], and at S136 by AKT [42]. Whereas there was nosignificant difference between BAD S112 phosphoryl-ation between DM and TM (Fig. 4a), phosphorylatedS136 levels were lower in TM cells than in all other cells(p < 0.001, Fig. 4b). It was not possible to reliably quan-tify these levels relative to BAD itself since BAD signalswere too low (not shown). One of the best attested sub-strates of AKT is glycogen synthase kinase 3 beta(GSK3β), which is phosphorylated on S9 by Akt leadingto inactivation of GSK3β [43]. Levels of phosphorylatedGSK3β S9 were elevated in WT over MP cells, elevatedonce more by removing the inhibitory CK2 sites in theDM mutation, and reduced by the further mutation ofY180 in TM mutant cells (Fig. 4c-e). These resultsstrengthen the model arising from pathways enrichmentFig. 5 PGRMC1 phosphorylation status affects glucose uptake, lactate secretion, and mitochondrial function. a Glucose uptake by cell lines usingthe Cayman “Glycolysis” kit. The boxplots represent four technical replicates each for each of independent stably transfected sub-lines 1–3 of eachcondition WT, DM and TM (lines from Fig. 1). i.e. n = 4 × 3 = 12 per condition. For MP cells, 12 replicates of the MP cell line were performed. Therewas significant difference between all means (Kruskal-Wallis p < 0.0001, pairwise two-sample Kolmogorov-Smirnov Tests p < 0.0002). b Lactatesecretion by cell lines. Details follow (A), except duplicates of each stable cell line were measured (n = 3 × 2 = 6 per condition). Inter cell-typecomparison tests revealed that the means of all pairwise comparisons were significantly different from one another (ANOVA, post-hoc Dunnet’sT3, p < 0.003), except the WT-TM comparison which was not significant (p = 0.211). c The maximal respiratory capacity of mitochondria is affectedby PGRMC1 phosphorylation status. Mitochondrial oxygen consumption of respective independent clonal stable lines C1 of each PGRMC1-HAmutant condition WT, DM and TM. Arrows indicate the time of addition of ATP synthase inhibitor oligomycin, Δψm uncoupler FCCP, andelectron transport chain inhibitors rotenone & antimycin A. OCR: oxygen consumption rate (pmol/min normalized per μg protein), n = 5;mean +/− s.d.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 8 of 24analysis, that PI3K/AKT pathway which is activated inDM cells requires phosphorylated PGRMC1 Y180, andis therefore attenuated in TM cells.PGRMC1 phosphorylation affects FAK activation and HSFlevelsPathways mapping (File S5) suggested that transcriptionfactor HSF1 activity could be involved in the differencebetween WT v. DM (adjP = 0.006). HSF1 has been linkedwith Focal Adhesion Kinase (FAK) activity [44], andFAK activity is dependent upon Rho/ROCK signalingwhich influences focal adhesion dynamics and tumor cellmigration and invasion [45]. RPPA measurementsshowed that FAK1 tyrosine phosphorylation and in-creased HSF1 levels were all significantly elevated in WTand TM (Fig. 4f-h). Notably, this profile resembled thedifferential proteomics profile of Fig. 3, rather than theROCK-dependent rounded morphology of Fig. 1d, andonce more provided validation of a pathways prediction.Notably, whereas FAK activity has been associated withincreased motility [45], its phosphorylation did not differsignificantly between WT and DM cells, which exhibitedaltered motility. Elevated FAK phosphorylation in TMcells, which exhibited lower migration than DM, couldpossibly be due to feedback regulation onto FAK. Thebasis for this merits future investigation, however pre-sumably reflects an impingement of PGRMC1 phosphor-ylation status on the migratory process.Vinculin elevated in DM is required for cell motilityProteins of the actin cytoskeleton were more abundantin DM cells (Fig. S4A), one of which was vinculin (Fig.S1B), an actin filament-binding protein associated withcell differentiation status, locomotion, and PI3K/AKT,E-cadherin, and β-catenin-regulated WNT signaling incolon carcinoma [46, 47]. We attenuated vinculin levelsin MP, WT, and DM cells via shRNA. Scrambled shRNAcontrol (shScr) DM cells exhibited elevated scratch assaymotility relative to MP cells. However in anti-vinculinshRNA (shVCL) cells, both MP and DM cell motilitywas reduced (Fig. S4). These results are consistent withelevated levels of proteins involved in the actin cytoske-leton (Fig. S4A) contributing directly to the enhancedmotility of DM cells (Fig. 1f-g). However, that hypothesisremains untested except for vinculin. Future studiesshould examine the effect of PGRMC1 signaling on actincytoskeleton properties.PGRMC1 affects glucose metabolismFigure 3 predicted altered glycolysis activity, which weinvestigated by glucose uptake and lactate production as-says. Expression of all PGRMC1-HA proteins (WT, DM,and TM) led to significantly lower levels of both mea-sures relative to MP, with DM cells exhibiting the lowestlevels (Fig. 5a-b). It is unclear whether the MP/WT dif-ference is due to PGRMC1-HA expression or to thehygromycin selection of WT cells. However, for theother cells PGRMC1 phosphorylation status regulatesboth features, consistent with recently reported regula-tion of Warburg metabolism by PGRMC1 [25], whichour data indicates can be regulated by PGRMC1 phos-phorylation status.PGRMC1 phosphorylation affects mitochondrial functionFigure 3 also implied that mitochondria may be affectedby PGRMC1 phosphorylation status. Naphthalimide-flavin redox sensor 2 (NpFR2) is a fluorophore targetedto the mitochondrial matrix. Its fluorescence is elevatedapproximately 100-fold when oxidized, providing anassay for mitochondrial matrix redox state [48]. NpFR2revealed that the matrix of WT and TM cells was moreoxidizing than MP and DM cells (Fig. S5A, B), whichcorresponded with the elevated expression of manynuclear-encoded mitochondrial proteins in Fig. 3.We then examined mitochondria using the fluorescentmarker MitoTracker, whose affinity for mitochondria isaffected by mitochondrial membrane potential (Δψm)[49]. Flow cytometry revealed the presence of two popu-lations of MitoTracker-binding cells in each cell type:low and high MitoTracker-binding (Fig. S5C-E). Thesepopulations were about equal for MP, WT, and TMcells, however, DM cells exhibited overall lower relativefluorescence level in each population Fig. S5C-D) and ahigher proportion of cells with higher MitoTracker bind-ing (Fig. S5C,E). Notably, higher levels of mitochondrialproteins in WT and to some extent TM cells apparentlydid not correspond with higher Δψm caused by activelyrespiring mitochondria.Relative to MP cells, the maximal respiratory rate wasreduced (between 2 and 3 fold in Fig. 5c) by expressionof DM or WT PGRMC1-HA, but not TM cells (Fig. 5c).The relative profiles of basal (Fig. S5F) and maximal(Fig. S5G) respiratory rates for WT and DM cells weresimilar, with TM cells exhibiting rates intermediate tothose of MP. This profile was observed on three inde-pendent comparisons WT/DM/TM. The single ex-periment including MP is shown. We conclude thatthe altered abundance of mitochondrial proteinsdue to PGRMC1 phosphorylation status detected inFig. 3 was accompanied by altered mitochondrialfunction. However, the relationship is not as simpleas lower glucose uptake being associated withhigher mitochondrial oxygen consumption, orNpFR2 redox status (a function of the rates ofsuperoxide leakage from complex IV, and antioxi-dant defense systems [50]) and may involve alter-ations in mitochondrial permeability to protons orother uncoupling mechanisms, for instance byThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 9 of 24altered cholesterol content [21], and/or altered elec-tron transport chain properties. This result clearlyreveals an important area requiring further studythat was revealed by our proteomics study.PGRMC1 phosphorylation affects mitochondrialmorphology and functionMitochondria were further explored by measuring mito-chondrial content (area per cell), size (perimeter), andFig. 6 PGRMC1 phosphorylation influences mitochondrial and cell morphology together. a Images show cells with close to average form factor (avFF) valuesfor elongate or rounded cells of each cell type. The numbers of cells scored in each class are indicated. All cell images are reproduced to the same scale.b Elongated and round cells from (a) examined according to cell type. The top left panel shows average percentage elongated cells scored for threeindependent biological replicate cell lines per cell type across all 3 replicates (n(Σ1–3)) corresponding to the given n values in (A). The remaining panels showingthe distribution of cells in each shape and type category of Area, Perimeter and FF. Kruskal-Wallis analysis revealed no significant differences at the p< 0.05 levelbetween any elongated cell type comparisons (ns), whereas all cell types exhibited significant differences for round cells (Kruskal-Wallis, P< 0.001). Theaccompanying Kruskal-Wallis post-hoc pairwise comparison p-values for round cells are given in the respective Tables. c AvFF per cell plotted including cellshape. The values where cells transition between round and elongated morphology (avFF 2.2–2.6) are indicated by dotted linesThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 10 of 24morphology, or form factor (FF). FF is a parameter de-rived from individual mitochondrial area and perimeter,where higher values correspond to a greater degree offilamentous than fragmented mitochondria [51, 52].Representative images of mitochondria are shown inFig. 6a. Numbers of mitochondria per cell varied greatly,with no significant differences detected between celltypes (not shown). Over the entire data set, elongatedcells exhibited greater mitochondrial area, larger mito-chondria, and greater average FF (avFF) (Kolmogorov-Smirnov p < 0.0001; not shown). When analyzed accord-ing to PGRMC1 status (cell type), MP and WT cellswere predominantly elongated, and DM and TM werepredominantly rounded, as expected (Fig. 6a-b). Wedetected no significant differences in average mitochon-drial area, perimeter or avFF between cell types forelongated cells, however rounded cell types exhibitedsignificant differences between the cell types for area,perimeter, and avFF (Fig. 6b). All cells with avFF < 2.2exhibited rounded cell shape, while all cells with avFF> 2.6 exhibited elongated shape (Fig. 6c). The observedavFF-associated transition from round to elongated cellshape was discrete for all cells except WT, occurring atavFF = 2.4 (MP), 2.2–2.6 (WT), 2.7 (DM) and 2.6 (TM).Notably, the single elongated TM cell also exhibited thehighest avFF value for TM (Fig. 6c). Holo-tomographictime-lapse videos [53] for each cell type show livemitochondria (File S7). PGRMC1 phosphorylation statusprobably influences mitochondrial content, size, and FF(degree of filamentation) by the same mechanisms thataffect cell shape, consistent with the proposed influence ofcytoskeleton on mitochondrial morphology and function[54]. We conclude that our proteomics pathways analysisrevealed biologically relevant functional differences inmitochondria, once more indicating a new area thatrequires further research.PGRMC1 Y180 is required for subcutaneous mousexenograft tumor growthNo significant differences in cell proliferation betweencell types were observed in culture IncuCyte imaging(Fig. 7a), or repeated MTT assays (not shown). Weestablished subcutaneous xenograft tumors in replicatemice carrying each of sub-lines 1–3 of WT, DM andTM (e.g. 4x line 1, 3x line 2, 3x line 3, n = 13 perPGRMC1-HA condition), as well as n = 5 mice with cellsexpressing the single Y180F mutant [24], or MP cells.Tumors produced by both TM and Y180F cells weresignificantly smaller than those produced by WT or DMcells (Fig. 7b-c), indicating that PGRMC1 Y180 wasrequired for optimal tumor growth, and demonstratingthat the cellular responses to altered PGRMC1 phos-phorylation strongly influences cancer biology. All WT,DM, TM and Y180F tumor tissue expressed thePGRMC1-HA proteins (not shown), and therefore arosefrom the injected cells. There were no obvious differ-ences in histology between cell types based uponhematoxylin and eosin staining (not shown).Fig. 7 PGRMC1 Y180 contributed to growth of subcutaneous mousexenograft tumors. a There are no differences in cell doubling time inculture. Replicates of stable cell lines 1–3 for each condition weremeasure to give n = 6. There were no significant differences betweencell conditions. b Typical tumors produced by each respective class ofcell. c Box plot of the distribution of tumor sizes among mice injectedwith 2 × 106 cells of each of the cell lines. For WT, DM and TM theresults depict 4x each of lines 2 and 3, and 5 mice from line 1. MP andY180F each represent 5 replicates of a single cell line. The box showspairwise post-hoc Dunnet’s T3 p-values after one way ANOVAThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 11 of 24DiscussionWe report new biology associated with the phosphoryl-ation status of PGRMC1-HA proteins from Fig. 1a,which profoundly affected cell morphology and migra-tory behavior. The morphotypic change from WT toDM resembled MAT in MP cells, being sensitive toROCKI (Fig. 1d). The DM and TM altered morphologyis dependent upon activated ROCK, which leads to stiff-ening of cortical actomyosin [55]. In human glioma cellsover-expression of CD99 is implicated in MAT, resultingin rounded morphology, increased Rho activity, and en-hanced migration [56]. These properties superficially re-semble the phenotype of our DM cells (Fig. 1e-f),however DM cell migration involved pseudopodia andcell adhesion, as evidenced in cell migration videos (FileS1). Very little else is known at the molecular level aboutthe events that promote MAT and altered cell motility[57–60]. We provide here a global proteomics expres-sion study of a possibly MAT-related process, and showthat PGRMC1 phosphorylation status dramaticallyaffects cell morphology and behavior.Proteomics pathways enrichment analysis suggestedthat PGRMC1 also influenced mitochondria, leading usto detect altered function and morphology (Fig. 6). Itwas previously shown that PGRMC1 co-evolved with asuite of mitochondrial proteins [21]. The elongated cellmorphology that predominated in MP and WT cells wasassociated with a higher index of filamentous rather thanfragmented mitochondria. Cells with rounded morph-ology and more fragmented mitochondria predominatedin DM and TM cells (Fig. 6). Such changes in mitochon-drial function are driven by altered relative rates of mito-chondrial fission and fusion, leading to mitochondrialfragmentation or elongated hypertubulation, respectively[61]. Fragmented mitochondria are associated withpathological conditions including cardiovascular andneuromuscular disorders, cancer, obesity, and theprocess of aging, associated largely with altered cell dif-ferentiation [61]. One of the proteins more abundant inWT and TM cells was Opa1 (O60313) (File S6), a proteinknown specifically to regulate mitochondrial fission/fusion,and one that has been reported to interact directly withPGRMC1 [62]. This suggests areas of future research.The strongest driver of mitochondrial morphologyappears to have been cell shape, or vice versa (Fig. 6).The cytoskeleton is thought to influence mitochondrialfunction and morphology [54], and proteomics pathwaysmapping suggested cytoskeletal changes. WT cells hadelevated levels of Tubulin 1alpha (Q71U36), and de-creased levels of Tubulin alpha-1C (Q9BQE3), Tubulinbeta-2A chain (Q13885), Tubulin beta-4B chain(P68371) (Fig. S1E) relative to DM cells. Rounded DMcells exhibited more abundant type II cytoskeletal kera-tin 8 (P05787), type I cytoskeletal keratin 18 (P05783),and type I cytoskeletal keratin 19 (P08727) relative toMP cells (which are epithelial markers [63]), with WTand TM cells exhibiting intermediate keratin levels (seeFile S6 associated with Fig. 3). These keratin abundanceprofiles paralleled the elevated levels of proteins in anenriched actin cytoskeleton pathway in DM cells (Fig.S1B). We have recently demonstrated that many pro-teins related to the actin cytoskeleton are co-immunoprecipitated with PGRMC1 in MP cells [64],and that the PGRMC1 membrane-associated progester-one receptor (MAPR) family cytochrome b5 domaincontains a region of predicted coiled-coil formation thatresembles regions from several myosins [22], consistentwith an influence of PGRMC1 on actin-related functionsand interacting with proteins that interact with myosins.DM cells also displayed elevated levels of proteins in theT-complex protein-1 ring complex (TRiC, also known asCCT) (Fig. S1E), which contributes to folding of proteinsincluding actin and microtubules, and influences deregu-lated growth control, apoptosis, and genomic instability[65, 66]. TRiC additionally contributes to obligatorygrowth/survival functions in breast [67] and liver [68] can-cers. TCP1 (P17987, Fig. S1E) expression is driven byoncogenic PI3K signaling in breast cancer [67], and weobserve both elevated PGRMC1-dependent PI3K/AKTactivity and TRiC abundance in DM cells (Fig. S1E).In summary, many of the mitochondrial differencesobserved could be attributable to altered cytoskeletalproperties. However, the differential mitochondrial func-tions of Fig. 5 and Fig. S5 did not correspond well withcell shape, indicating that PGRMC1 also changes com-plex causative processes driven by more than mitochon-drial morphology. In this context, we note the recentobservation that PGRMC2 mediates heme shuttling tothe nucleus of adipocytes, leading to altered transcrip-tion of genes encoding mitochondrial proteins [69]. It ispossible that PGRMC1 exerts related effects in thesecells, which merits future study.Results depicted in Fig. 3 (as presented in File S5 andFile S6) revealed that both the ATP synthase subunitbeta of the F1 catalytic domain as well as the F0 protonpore domain were up-regulated in WT cells relative toDM cells (Fig. S1C). It is possible that the higher Δψmof DM cells is related to low levels of F0/F1 ATPase pro-ton channels (Fig. S1C), resulting in relatively inefficientproton gradient clearance. Mitochondrial cholesterol de-creases the permeability of the inner membrane to pro-tons, increasing the efficiency of electron transport chainyield [21]. Altered metabolism could be due to activityof the PI3K/AKT pathway, suggesting direction for fur-ther work that will be required to explain the mecha-nisms underlying the observed responses.Loss of PGRMC1 affects the SREBP-1/fatty acidhomeostasis system [70], and PGRMC1 influences cellThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 12 of 24surface localization of insulin receptor and glucosetransporters [27]. Chemical proteomics showed thatPGRMC2 but not PGRMC1 promotes adipogenesis in 3T3-L1 preadipocytes following a gain of function inter-action with a novel small molecule which displacedheme [71]. It will be interesting to examine whether thattreatment mimics the effects of phosphorylation.(PGRMC2 possesses cognates to PGRMC1 Y180 andS181, as well as heme-chelating Y113 [20].)Our results indicate that PI3K/AKT signaling in DMcells required PGRMC1 Y180, which was the sole differ-ence to TM cells. PGRMC1 has long been recognized asa modulator of AKT activity, with cell type-specific ef-fects [24, 27, 72–76]. This predicted activation of signalsby removal of the putative inhibitory CK2 consensussites in the DM protein [19, 24] was dependent uponY180 because TM (which differs from DM by a singleoxygen atom) exhibited a protein expression profile thatwas more similar to WT than DM (Fig. 3). Furthermore,our in vivo studies evaluating the growth of subcutane-ous tumors confirmed the importance of Y180 in an in-tact biological system thus supporting and extending ourin vitro observations to reveal the central importance ofY180 (Fig. 7). Because of the high conservation of Y180in eumetazoan animals [22], these results are likely tobear relevance to other species, and multiple humandiseases [2]. In a methylomics study of these cells in thecompanion paper, the most significantly down-regulatedKEGG pathway in the TM/DM comparison was PI3K-AKT [36], consistent with PI3K/AKT activation requir-ing PGRMC1 Y180, which was required for tumorgrowth (Fig. 7).Interestingly, PGRMC1 knockdown in human pluripo-tent stem cells (hPSCs) led to an increase in GSK3β in-hibitory phosphorylation [77]. Examination of PGRMC1phosphorylation status in that system is merited, wherePGRMC1 suppressed the p53 and WNT pathways tomaintain hPSC pluripotency. Similarly to our results,those authors concluded “that PGRMC1 is able to sup-press broad networks necessary for multi-lineage fatespecification.” Our hypothesis suggests that PGRMC1Y180 phosphorylation and PI3K/AKT activity could beassociated with elevated GSK-3β Ser9 phosphorylation andβ-catenin signaling as observed in in some cancers [47].Our initial hypothesis related to differential phosphor-ylation of PGRMC1, being potentially spatially and tem-porally associated with regulation of the Warburg effect[24]. It is notable that the Warburg effect resembles areversion to stem-cell-like metabolism [78]. While thismanuscript was in preparation, Sabbir showed thatPGRMC1 post-translational modification status inHEK293 cells responds to P4 treatment, which was ac-companied by a PGRMC1-dependent increase in gly-colysis [25]. Our cells exhibited a large reduction inglycolysis between MP and WT cells, which may be dueto either PGRMC1-HA expression or to hygromycinstable selection of WT cells. However phosphorylationstate differences between WT, DM and TM elicitedmarked differences in both glucose consumption andlactate secretion (Fig. 5). It is important to note that en-dogenous PGRMC1 is expressed in all of these cells, andso we cannot determine to what extent these mutatedproteins may elicit effects on their own, or interfere withendogenous PGRMC1. However, we conclude thatPGRMC1 can modulate glycolysis, in line with our hy-pothesis, although the mechanism requires further re-search which should examine the possibility that PI3K/AKT activity elevated in DM cells negatively regulatesthe Warburg effect.ConclusionThe phospho-acceptor amino acid Y180 has been con-served in PGRMC1 proteins since the evolutionary ap-pearance of the differentiation-inducing gastrulation-organizer [22]. We therefore believe it likely thatthe PGRMC1 Y180-regulated modulation of metabolicand growth control that we have manipulated could re-present a major newly identified foundational axis ofanimal cell biology, whose perturbation is inconsistentwith the maintenance of differentiated states acquiredduring the subsequent evolution of complex body plans.Although we can confidently deduce the existence of aPGRMC1 signal network, as yet we have identified nei-ther immediate upstream PGRMC1 effectors nor down-stream targets. In a companion study [36], we show thatthe cells characterized in this paper differ dramatically ingenomic methylation and mutation rates. Future studiesshould urgently explore the relationship betweenPGRMC1 signaling and diseases such as cancer, diabetes,Alzheimer’s disease, and others [2].MethodsPlasmid preparationPlasmids pcDNA3.1-PGRMC1-HA (wild type: WT),pcDNA3.1-PGRMC1-HA_Y180F (Y180F), pcDNA3.1-PGRMC1-HA_S57A/S181A (double mutant, DM) andpcDNA3.1-PGRMC1-HA_Y180F (Y180F) [24] have beendescribed. The triple mutant S57A/Y180F/S181 (TM)was constructed by Genscript (Hong Kong) using theDM plasmid as the template and introducing codonTTC for F180. PGRMC1-HA open reading frames werereconfirmed by DNA sequencing at Monash MicromonDNA Sequencing Facility (Clayton, Vic., Australia) usingthe 5′ T7 and 3′ BGH sequencing primers specific forthe parental vector. PGRMC1-HA plasmids were trans-formed into Escherichia coli Top10 strain, and culturedovernight at 37 °C on 1% agar plates containing Luriabroth (LB) media (Invitrogen) and 50 μg/mL ampicillin.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 13 of 24A single colony was picked and bacteria were grown in250 ml culture by aeration overnight at 37 °C in LBmedia. Plasmid DNA was isolated by GeneJet MaxiprepKit (ThermoScientific) following the manufacturer’sprotocol. Plasmid DNA concentration was measured byusing Nanodrop (Thermo Scientific).Cell cultureMIA PaCa-2 (MP) cell identity was verified as MIA PaCa-2(ATCC CRL-1420) by the MHTP Medical Genomics Facil-ity (Monash University, Melbourne) following the ATCCStandards Development Organization document ASN-0002for cell line identification via short tandem repeat profiling.MP cells were maintained in Dulbecco’s Modified Eagle’smedium (DMEM-high glucose, Sigma-Aldrich, D5796)supplemented with 10% Foetal bovine calf serum (Sigma-Aldrich, F9423) and 1% penicillin-streptomycin (Sigma-Al-drich, P4333) (complete DMEM) at 37 °C and 5% CO2 in a150i CO2 incubator (Heracell, Lane Cove NSW). Cell doub-ling times were estimated by 3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide (MTT) assay or by Incu-Cyte, as described [79]. Mitochondrial respiratory capacitywas measured using a Seahorse Extracellular Flux analyzerXF24 (Seahorse Biosciences).Transfection and stable cell line generationUnder our culture conditions, MP cells exist in cultureas flattened adherent cells of mesenchymal shape, a mi-nority of rounded adherent cells, and a small populationof rounded suspension cells. Plating the suspension cellsregenerates a similar population distribution (notshown). On the day before transfection, 2 × 106 MP cellswere seeded onto a 6-well plate. The cells were trans-fected at 70–80% confluency. Before transfection, cellswere washed with Dulbecco’s phosphate-buffered saline(PBS, Sigma-Aldrich, D8537) and maintained inantibiotic-free Dulbecco’s Modified eagle Medium highglucose (DMEM, Sigma-Aldrich, D5796) containing 10%bovine calf serum (Sigma-Aldrich, 12133C) and 1% peni-cillin/streptomycin (Sigma-Aldrich, P4458) (completeDMEM). In separate transfections, 4 μg each of respect-ive PGRMC1-HA plasmids (WT, DM or TM) and Lipo-fectamine 2000 (Life Technologies, 11,668–019) weremixed at 1:2 ratio and incubated for 25 min at roomtemperature. The mixture was then added drop-wise tothe wells of the culture plate. After 6 h of incubation,cells were washed with PBS and cultured at 37 °C and5% CO2 in complete DMEM for 48 h, after which cellswere harvested and plated in three fold limiting dilutionin complete DMEM containing 50 μg/ml Hygromycin B(EMD Millipore, 400,052) in 96 well plates. Cells werecultured at 37 °C and 5% CO2 for 2 weeks, with regularmedia changes containing complete DMEM with Hygro-mycin B every 3 days to select for stable integrationevents. Typically 8 independent stably transfected celllines were expanded for each of PGRMC1-HA WT, DMand TM and 3 lines with similar levels of PGRMC1-HAexpression were selected by Western blot.Cells were frozen 0.5–1.0 mL at − 80 °C in Bambanker(Novachem, #306–14,684) at 2 × 106 cells/mL. Frozencells were introduced back into culture by thawing at37 °C for 20 s followed by addition to 5 ml of completemedia and low speed centrifugation at 180 x g for 3mins at 25 °C. Pelleted cells were resuspended in 6 mLfresh complete media and seeded in 25 cm2 flasks.Because of the dramatic effects observed, MP cells areincluded in our experiments as a literature referencepoint. MP differ from WT cells by not having undergonehygromycin selection, and by lack of overexpression ofPGRMC1-HA. Therefore we cannot ascribe differencesbetween MP and WT cells to PGRMC1-HA expression.The effects of the DM and TM PGRMC1 mutations areassessed relative to WT control levels.shRNA lentiviral productionLentiviral-delivered shRNAs were constructed using Mis-sion TRC2-pLKO-Puro series lentiplasmids (SHCLND,Sigma-Aldrich) targeting ERR1/ESRRA (TRCN0000330191,GAGAGGAGTATGTTCTACTAA), vinculin (shVCL)(TRCN0000290168, CGGTTGGTACTGCTAATAAAT)or non-target scrambled shRNA (shScr) (SHC202;CAACAAGATGAAGAGCACCAA). MP and DM shVCLcells were obtained at first attempt, however, despite threeattempts, WT cells could not be established. To generatevirus particles, we co-transfected HEK293 cells with theshRNA plasmids and helper plasmids using Lipofectamine2000 (Invitrogen, 11668-027). Prior to transfection, 6 wellplates were treated with 50 μg/mL D-Lysine overnight.Next day, 1 × 106 HEK293 cells were seeded per well andincubated overnight at 37 °C in complete medium. Trans-fection mixture A contained 4 μg plasmid mixture consist-ing of 2.5 μg target or shScr shRNA lentiplasmid, 0.75 μgPax, 0.3 μg Rev. and 0.45 μg VSV-G helper plasmids [80]in 250 μL antibiotic-free medium. Transfection mixture Bcontained 8 μL Lipofectamine and 242 μL antibiotic-freemedium. After 5min incubation at 25 °C, mixtures A andB were gently mixed and incubated for 25min at 25 °C.HEK293 culture medium was removed, cells were washedwith PBS, and 2mL fresh antibiotic-free medium wasadded followed by addition of combined transfection mix-ture to the cells, dropwise with gentle shaking, followed byincubation for 6 h at 37 °C. After incubation, the mediumwas replaced with complete medium overnight at 37 °C.Virus particles were harvested by collecting culturemedium followed by the addition of new medium every24 h for 72 h. Collected media for each culture werepooled, filtered through a 22 μM filter, aliquoted into1 mL fractions, and frozen at − 80 °C.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 14 of 24shRNA lentiviral transductionBriefly, 1 × 105 MP, WT, or DM cells per 24 plate wellwere seeded in 1 mL complete DMEM medium andgrown to 60% confluency. The medium was removed,and replaced by 1 mL of medium per well, containing 2-fold serially diluted virus particles in adjacent wells, plus5 μg/mL Polybrene (hexadimethrine bromide, Sigma-Aldrich 107689) to enhance viral transduction. Afterincubation for 24 h, the medium was removed and thecells were washed twice with PBS after which freshmedium was added supplemented with 1.5 μg/mL Puro-mycin, which was replaced every 48 h for 1 week. Cellsfrom wells transduced with the lowest dilutions ofrespective virus particles that survived selection wereexpanded, and stocks frozen at − 80 °C in Bambanker.Scratch migration assayMP cells or stable transfected monoclonal MP cell linesexpressing PGRMC1-HA WT, DM or TM proteins (1 ×104 cells) were seeded in a 24 well plate. The monolayerof cells at more than 90% confluency was subjected toserum starvation for 2 h. A scratch was created in the mid-dle of the monolayer by a sterile p200 tip and washedtwice with PBS to remove floating cells. Complete mediawas then added. The cell monolayer was incubated for 36h to allow cell migration into the scratched area. Photo-graphic images were taken at 0 and 36 h using an invertedphase microscope (Nikon Eclipse Ti-U). Cells in the boxedareas of Fig. 2a were manually scored from printed images.Cell treatments included 125 nM Y-27632 dihydrochloride(Abcam, ab120129, Rho Kinase inhibitor: ROCKI) or ve-hicle control (DMSO). For video files the cells were incu-bated at 37 °C and 5% CO2 for 36 h in a stage topelectrically heated chamber (Okolab H301-NIKON-NZ100/200/500-N) including transparent heated lid(H301-EC-HG-LID), with a 24-well Nunc/Greiner platebase adapter (24MW-NUNC) and a chamber riser, for aworking distance of 28mm. The chamber was regulatedby a Control Unit (Okolab H301-TC1-HMTC) withDigital CO2 controller (Okolab DGT-CO2 BX) and AirPump (Okolab, OKO-AP) and was inserted to a Nano-Z100-N Piezo stage (Mad City Labs) on a motorized XYstage (Nikon TI-S-ER). Images were taken every 10minfor 36 h with a 10× (0.45 NA) Plan Apo objective usingthe transmitted light detector (TD) on a Nikon Ti EclipseConfocal microscope controlled by NIS Elements V4.10software (Nikon).Proteomics sample preparationThree independent stable transfected lines of eachPGRMC1-HA-expressing cell type, as well as triplicatesof the MP parental cell line, were measured in technicalreplicate data-dependent and independent data acquisi-tion SWATH-MS modes on a 5600 TripleTof™ massspectrometer (ABSciex). Global proteomics analysis wascarried out at the Australian Proteome Analysis Facility(APAF). Cells were grown in Wagga Wagga to 80% con-fluency in 75 cm2 flasks. Three separate cultures of MPcells (passages 8, 9 & 11) and three lines of eachPGRMC1-HA WT, DM and TM cells were used (inde-pendent biological triplicates). Cells were harvested andfrozen cell pellets shipped on dry ice to APAF for Massspectrometric analysis. Cell pellets were lysed using200 μL of sodium deoxycholate buffer (1% in 0.03M tri-ethyl ammonium bicarbonate), and DNA digested using0.5 μg of benzonase. Direct detect assay (EMD Millipore,DDAC00010-8P) was performed on the samples and100 μg of each sample was taken for digestion. Sampleswere reduced with dithiothreitol (5 mM), alkylated withiodoacetamide (10 mM) and then digested with 4 μgtrypsin for 16 h at 37 °C. The digested sample was acid-ified and centrifuged to remove the sodium deoxycho-late. Samples were then dried and resuspended in 100 μLof loading buffer (2% acetonitrile 0.1% formic acid). Indi-vidual samples for SWATH analyses were diluted 1:4into loading buffer and transferred to a vial. Each samplewas measured in technical replicate. For IDA runs a poolwas made for each group (MP, WT, DM, and TM) bytaking equal portions from each biological replicate anddiluting 1:4 in loading buffer.Proteomic information dependent acquisitionTryptic peptides were analyzed on a 5600 TripleTof™mass spectrometer (ABSciex). Chromatographic separ-ation of peptides was performed on a NanoLC-2DplusHPLC system (Eksigent, Dublin, CA) coupled to a self-packed analytical column (Halo C18, 160 Å, 2.7 μm,75 μm× 10 cm). Peptide samples (4 μg of total peptideamount) were loaded onto a peptide trap (Opti-trap Cap0.5 mm × 1.3 mm, Optimize Technologies) for pre-concentration and desalted with 0.1% formic acid, 2%acetonitrile, at 10 μL/min for 5 min. The peptide trapwas then switched into line with the analytical columnand peptides were eluted from the column using linearsolvent gradients with steps, from 98% Buffer A (0.1%formic acid) and 2% Buffer B (99.9% acetonitrile, 0.1%formic acid) to 90% Buffer A and 10% Buffer B for 10min, then to 65% Buffer A and 35% Buffer B at 500 nL/min over a 78 min period. After peptide elution, the col-umn was cleaned with 95% Buffer B for 15 min and thenequilibrated with 98% Buffer A for 15 min before thenext sample injection. The reverse phase nanoLC eluentwas subject to positive ion nanoflow electrospray ana-lysis in an IDA mode. Sample analysis order for LC/MSwas DM1, DM2, DM2, DM3, DM3, MP2, MP2, MP3,MP3, TM2, TM2, TM3, TM3, WT1, WT1, WT2, WT2,WT3, WT3, DM1, MP1, MP1, TM1, TM1.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 15 of 24In the IDA mode a TOFMS survey scan was acquired(m/z 350–1500, 0.25 s), with the 10 most intense multi-ply charged ions (2+ − 5+; counts > 150) in the surveyscan sequentially subjected to MS/MS analysis. Theselected precursors were then added to a dynamic exclu-sion list for 20s. MS/MS spectra were accumulated for50 milliseconds in the mass range m/z 100–1500 withrolling collision energy.Proteomic data independent acquisition (SWATH)Samples were analyzed by SWATH-MS [81] proteomicsprofiling in duplicate, with the chromatographic condi-tions as for IDA analysis above. The reverse phasenanoLC eluent was subject to positive ion nanoflowelectrospray analysis in a data independent acquisitionmode (SWATH). For SWATH MS, m/z window sizeswere determined based on precursor m/z frequencies(m/z 400–1250) in previous IDA data (SWATH variablewindow acquisition, 60 windows in total). In SWATHmode, first a TOFMS survey scan was acquired (m/z350–1500, 0.05 s) then the 60 predefined m/z rangeswere sequentially subjected to MS/MS analysis. MS/MSspectra were accumulated for 96 milliseconds in themass range m/z 350–1500 with rolling collision energyoptimized for lower m/z in m/z window + 10%. Tominimize instrument condition caused bias, SWATHdata were acquired in random order for the sampleswith one blank run between every sample.SWATH library generationThe LC-MS/MS data of the IDA data were searchedusing ProteinPilot (version 4.2) (Sciex) and combinedinto a single search report file. The files were searchedagainst Human entries in the Swissprot 2014_04 data-base (released 16/04/2015, containing 545,388 entries).The search parameters were selected as follows: iodoace-tamide cysteine alkylation, trypsin digestion, Triple TOF5600 instrumentation, biological modifications, thoroughsearch and false discovery rate enabled.SWATH data processingSWATH data were extracted using PeakView (version2.1, Sciex) with the following parameters: Top 6 mostintense fragments of each peptide were extracted fromthe SWATH data sets (75 ppm mass tolerance, 10 minretention time window). Shared peptides were excluded.After data processing, peptides (max 50 peptides perprotein) with confidence ≥99% and FDR ≤1% (based onchromatographic feature after fragment extraction) wereused for quantitation. The extracted SWATH proteinpeak areas were normalized to the total peak area foreach run and subjected to t-test to compare relative pro-tein peak area between the samples. Protein t-test withp-value smaller than 0.05 and fold change larger than1.5 were highlighted as differentially expressed. The ana-lysis of four different cell types was treated as six separ-ate paired comparisons: 1) MP vs. WT, 2) MP vs. DM,3) MP vs. TM, 4) WT vs. DM, 5) WT vs. TM, and 6)DM vs. TM. Additionally, a similar method for deter-mining differential expression was run at the peptidelevel, with the peptide level fold changes thenaveraged for each protein. The peptide level analysisis more conservative, as peptide level p-values areonly generated when a protein was identified by atleast two proteins. Peptide level data were used forthis study.WebGestalt enrichment analysesGene Ontology (GO) and pathway enrichment analysiswere conducted on differentially abundant proteins fromFigure S2D using the WEB-based GEne SeT AnaLysisToolkit (WebGestalt) platform (http://bioinfo.vanderbilt.edu/webgestalt/) [82]. Figures or Supplementary Figuresdepicting or relying on WebGestalt results describe sup-plementary data files available online that contain the re-spective original WebGestalt analysis files. Attemptingto identify the most important driving contributions tothe phenotypic alterations observed, a complementaryWebGestalt analysis was also performed at the adjP<0.001% level, but employing the subsets of proteinswhich were detected to be either significantly up- ordown-regulated (“red” and “blue” lists of proteins foreach comparison from File S1). The UniProt protein IDsfrom respective proteomics comparisons were uploadedas text files which accompany the respective WebGestaltsupplementary data files. KEGG, Pathway Commons,and Wiki Pathways were analyzed at the 5% level. GOand transcription Factor analyses were also performedfor comparisons of differentially up- or down-regulatedproteins for each pair-wise comparison between celltypes at the 0.1 and 10% levels. The following WebGes-talt settings were employed: Organism: hsapiens; gene IdType: entrezgene; Reference Set for Enrichment Ana-lysis; entrezgene_protein-coding; Significance Level:(variable see individual analysis descriptions), StatisticalMethod: Hypergeometric, Multiple Test Adjustment:Benjamini-Hochberg (BH), Minimum Number of Genesfor a Category: 3.Pathway enrichment analyses can consider either a) alldifferential proteins together (including both up- anddown-regulated proteins in the one analysis), or b) canexamine the higher abundance and lower abundanceproteins in separate analyses for each comparison (“red”or “blue” analyses for each comparison). Our WebGes-talt pathway analysis strategy of Fig. 3 pursued the sec-ond of these alternatives. In the first analysis, all proteinsfound to be differential (both up- and down-regulated)between two samples in any six of the analyses wereThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 16 of 24entered as a protein list to WebGestalt, and pathway en-richment analysis was performed with multiple samplecorrection at the Benjamini and Hochberg (BH) adjustedp (adjP) significance level of 0.001 for KEGG, PathwayCommons, and Wikipathways for each of the six basiccell-type comparisons.Twelve separate WebGestalt analyses were performedfor Fig. 3, with enrichment analyses including KEGG,PC, WikiPathways, as well as Transcription Factor andGO cellular component. The schematic representationof Fig. 3 shows pathways identified as significantlyenriched for at least one of 12 comparisons (6x moreabundant “red” and 6x less abundant “blue”) between 1)MP vs. WT, 2) MP vs. DM, 3) MP vs. TM, 4) WT vs.DM, 5) WT vs. TM, and 6) DM vs. TM. The pathwaymapping for all significantly detected features betweenall comparisons detected at adjP< 0.001 is available asFile S5A. Another WebGestalt analysis was performed atadjP< 0.1 and the significances of each pathway identi-fied in File S5A at the adjP< 0.001 level were then re-corded for all comparisons at the adjP< 0.1 level in FileS5B. The analysis results are presented for reference asFile S4B. That information was used to assign statisticalsignificance in the pathways map for all features thatwere identified at the adjP< 0.001 level in any one com-parison, across all 12 comparisons at the adjP< 0.1 level.These data were then used to map all proteins from allpathways and all comparisons of File S5B to produce theoriginal image of Fig. 3 in Microsoft Excel, which isavailable as File S6 and contains all protein and pathwayidentities.Mapping pathways to the expression heat mapThe matrix of protein membership to pathway or func-tional group category is a resulting sparse matrix with 0/1 indicating that the respective protein is/is not presentin the respective category. This matrix was clusteredusing the hclust implementation in the R Base Package(www.r-project.org/), using a binary distance andcomplete linkage, to reorder the columns (pathways inthis case) according to the proportion of shared proteins.The resulting cladogram including overlapping featuresidentified by all WebGestalt analyses appears to the rightof pathways in Fig. 3, and with complete accompanyingprotein and pathway identities in File S5.Principal components analysis on proteomics resultsPrincipal component analysis was used to examinethe largest contributions to variation in the proteinmeasurements. Wilcoxon rank-sum tests were usedto identify the pathways that were positively ornegatively associated with the principal componentscores.Sample preparation for Western blotsApproximately 70% confluent cells in a T75 flask werewashed twice with chilled PBS buffer and incubated with500 μL radio immunoprecipitation assay buffer (RIPAbuffer) (Sigma-Aldrich, R0278) supplemented with pro-tease and phosphatase inhibitor cocktail (Thermoscienti-fic, 88668) following manufacturer’s recommendations.After scraping, the lysate was centrifuged at 8000 g for20 min (Hermle Centrifuge Z233 M-2) at 4 °C. Proteinconcentration was determined using the Pierce BCAprotein assay kit (ThermoFisher, 23225) following themanufacturer’s instructions. Twenty microgram cell ly-sates were each mixed with 2x Laemmli loading buffer(Sigma-Aldrich, S3401) at a 1:1 ratio to give final volume20 μL, followed by denaturation at 95 °C for 5 min in adigital dry bath heater. Lysates were loaded immediatelyto a 10% SDS-PAGE gel. For HA western we used 17well gels (Life technologies, NW04127), for PGRMC1Western blots we used 15 well gels (Bio Rad. 456–1069),for vinculin and ERR1 Western blots we used 10 wellgels (Bio Rad 456–1096). Electrophoresis was at 150 Vfor 45 min. Protein was transferred onto PVDF mem-branes (Bio-Rad, 1620174) with a Trans-Blot Turbotransfer system (Bio-Rad, Gladesville NSW) for 7 min byTrans-Blot® Turbo RTA Mini LF PVDF Transfer Kit(Bio-Rad, 1704274) or wet transferred in 25 mM Tris,192 mM glycine, 20% (v/v) methanol (pH 8.3) (1x Tow-bin buffer) at 20 V for 2.5 h on mini trans-blot cell (BioRad, 1703930) cooled on ice.Western blotsMembranes were blocked with TBS-T (0.1% Tween-20in 1× Tris-buffered saline) containing 5% WoolworthsInstant Skim Milk Powder (Woolworths, Wagga Wagga,NSW, Australia) for 1 h and incubated overnight at 4 °Cwith primary (1°) antibody. After washing 3 times withTBS-T, blots were incubated with secondary (2°) anti-body for 1 h at room temperature. Proteins were de-tected by the following methods. For chemiluminescencedetection the membranes were incubated with ClarityMax Western ECL Substrate (Bio-Rad, #1705062) for 5min for detection by enhanced chemiluminescence usinga Bio-Rad ChemiDoc MP imaging system (Bio-Rad, Gla-desville NSW). Fluorescence detection was performedon the ChemiDoc (at the indicated wavelength). Molecu-lar weight standard proteins for gels imaged for fluores-cence or chemiluminescence were detected on theChemiDoc using the IRDy680 channel. Colorimetric de-tection (Vinculin and ERR1 Western blots) was per-formed by incubation of membranes with 3 mLTetramethylbenzidine (TMB) (Sigma-Aldrich, T0565)for 5 min. These images were captured by MolecularImager Gel Doc XR+ System (Bio-Rad, GladesvilleNSW). Multi-channel ChemiDoc images were generatedThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 17 of 24with the Bio-Rad Image Lab Software. Some dual chan-nels images were manipulated in Adobe Photoshop CC2018 (Adobe Systems Inc.) by reducing intensity in ei-ther red or green channel to lower background in thepublished image. Adjustments were applied identicallyover all image pixels so as to not alter the relative inten-sities of any bands.The following primary (1°) and secondary (2°) antibodypairs were used (at the specified dilutions) with the indi-cated detection methods. For PGRMC1 Western: 1° goatanti-PGRMC1 antibody (Abcam, ab48012) (1:1000) and2° rabbit anti-goat secondary antibody (Abcam, ab6741)(1:4000) detected by chemiluminescence. After detec-tion, membranes were blocked with TBS-T overnightand then incubated with 1° mouse anti-beta actin(Sigma-Aldrich, A5541) (1:2000) and 2° goat anti-mouseIgG H&L (IRDye 800CW) (1:5000) detected by fluores-cence (IRDye 800CW). For HA epitope Western: 1°mouse anti-HA (Sigma-Aldrich, H3663) (1:2000) and 2°goat anti-mouse IgG H&L (IRDye 800CW) (dilution 1:5000) detected by fluorescence (IRDye 800CW), as well as1° rabbit anti-beta-actin (Cell Signaling, 4967) (1:2000)and 2° donkey anti-rabbit IgG (Abcam, ab16284) (1:2000)detected by chemiluminescence. For vinculin Western: 1°anti-vinculin (E1E9V) XP (Cell Signaling Technology, 13,901) (1:1000) and 2° donkey anti-rabbit IgG (Abcam,ab16284) (1:2000) detected colorimetrically. For ERR1Western: 1° rabbit anti-ERRα (E1G1J) (Cell Signaling,13826S) (1:1000) and 2° donkey anti-rabbit IgG (Abcam,ab16284) (1:2000) detected colorimetrically.Reverse Phase Protein Array analysisRPPA arrays using Zeptosens technology (Bayer AG, Le-verkusen, Germany) were used for analysis of signalingprotein expression and activity profiling as described [83–86]. For the analysis, flash frozen cell pellets were lysed byincubation with 100 μl cell lysis buffer CLB1 (Bayer,Germany) for 30min at room temperature. Total proteinconcentrations of the lysate supernatants were determinedby Bradford Assay (Coomassie Plus, Thermo Scientific).Cell lysate samples were adjusted to uniform protein con-centration in CLB1, diluted 10-fold in RPPA spotting buf-fer CSBL1 (Bayer) and subsequently printed as series offour dilutions (starting concentration at 0.3 μg/μl plus 1.6-fold dilutions) and in two replicates each. All sampleswere printed as replicate microarrays onto Zeptosenshydrophobic chips (Bayer) using a NanoPlotter 2 (GeSim,Grosserkmannsdorf, Germany) applying single droplet de-positions (0.4 nL volume per spot). After printing, the mi-croarrays were blocked with 3% w/v albumin, washedthoroughly with double distilled H2O, dried in a stream ofnitrogen and stored in the dark at 4 °C until further use.Protein expression and activity levels were measuredusing a direct two-step sequential immunoassay andsensitive, quantitative fluorescence read-out. A singlearray was probed for each protein. Highly specific and up-front validated primary antibodies were incubated at the re-spective dilution in Zeptosens assay buffer overnight (15 h)at room temperature. Arrays were washed once in assaybuffer and incubated for 45min with Alexa647-labeledanti-species secondary antibody (Invitrogen, Paisley, UK).Arrays were then washed as before and imaged using aZeptoREADER instrument (Bayer) in the red laser channel.Typically, six fluorescence images were recorded for eacharray at exposure times of between 0.5 and 16 s. Negativecontrol assays incubated in the absence of primary antibody(blank assays) were also performed to measure the non-specific signal contributions of the secondary antibody. Inaddition, one chip out of the print series was stained tomeasure the relative amount of immobilized protein perspot (protein stain assay). The following primary antibodies(provider and reagent number, dilution) were used: BAD(CST 9239, 1:200), BAD-P-Ser112 (CST 5284, 1:100),BAD-P-Ser136 (CST 4366, 1:100), FAK1-P-Tyr577 (Invi-trogen 44-614ZG, 1:100), FAK1-P-Tyr861 (Epitomics2153–1, 1:100), GSK3beta (CST 9315, 1:200), GSK3beta-P-Ser9 (CST 9336, 1:100), HSF1 (Epitomics 2043–1, 1:1000).After assay measurements (one protein per array), arrayimages and data were analyzed with the software Zepto-VIEW 3.1 (Bayer). For each array/antibody, the image takenat the longest exposure time without showing any satur-ation effect was analyzed with the spot diameters set to160 μm. Mean fluorescence signal intensity (MFI) of eachsample was calculated from referenced, background-corrected mean intensities of the single spots (eight spotsper sample) applying a linear fit and interpolating to themean of the four printed protein concentrations. Blank-corrected MFI signals of the samples were normalized forthe relative protein concentration printed on the chip toobtain normalized fluorescence intensity signals (NFI). NFIvalues were used for all subsequent statistical analyses.Glucose uptake & lactate production assayGlucose uptake and lactate production assays were per-formed by using commercially available kits from Caymanchemical (#600470, #700510) following manufacturer’sprotocols. Glucose uptake was measured with a FluostarOmega fluorescence microplate reader (BMG Labtech,Ortenberg, Germany) and lactate production wasquantified with a Molecular Devices Spectra Max 190microplate reader (Bio-Strategy P/L, Campbellfield,Vic., Australia).NpFR2 redox assayIntramitochondrial redox status was determined usingnaphthalimide flavin redox sensor 2 (NpFR2) [48]. MiaPaCa-2 and PGRMC1-HA-expressing stable cells (1 ×106) were suspended in 2 mL complete media andThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 18 of 24seeded in six well plates and cultured for 24 h at 37 °Cand 5% CO2. Cells were washed with PBS, trypsinsed,harvested, and resuspended in 1 mL of fresh media con-taining 25 μM NpFR2 in a 1.5 mL microcentrifuge tube,followed by incubation for 20 min at 37 °C. Cells werethen centrifuged in a microcentrifuge at 180 x g, the pel-let was resuspended once with 1 mL PBS followed byrecentrifugation, and the washed pellet was againresuspended in 1mL PBS. Five hundred microliter cellsuspension was loaded to a Gallios Flow Cytometer(Beckman Coulter) and fluorescence of 2 × 104 cells wasdetected using FL1 (green) channel.Immunofluorescence microscopyTo detect the expression of exogenous HA taggedPGRMC1 in Fig. 1e, cells were seeded on coverslips on asix well plate. The cells were washed with ice-cold PBS,mildly fixed with 3.7% formaldehyde for 5 min at 4 °C.The cells were then permeabilized with ice-cold 100%methanol for 10 min at − 20 °C, followed by overnightincubation with anti-HA tag antibody (Sigma, H3663).The cells were washed extensively and incubated withFITC conjugated secondary antibody (Sigma, F8521) indark for 1 h at 4 °C. Cells were washed three times withPBS and counterstained with DAPI mounting solution.Images were captured using a Nikon Ti Eclipse Confocalmicroscope (Nikon Australia Pty Ltd).Analysis of mitochondrial morphologyMitochondria were quantified for cell shape (elongated/round), mitochondrial content (sum of mitochondrialarea/cell), mitochondrial size (average perimeter/cell),and mitochondrial morphology or Formfactor (FF): ameasure where higher values correspond to a greaterlevel of filamentous mitochondria and lower values cor-respond to more highly fragmented mitochondria [52].FF (calculated as the P2/4πA) measures mitochondrialmorphology based on the perimeter and area of shape.The calculation takes in to account not only changes inlength, but also the degree of branching, making at anideal form of measurement for the quantification ofmitochondrial morphology.To measure FF, 1 × 105 cells were seeded onto Nunc176740 four well plates with a 22x22mm #1.5 glasscoverslip on the bottom. Cells were fixed and perme-abilized as above, then incubated with Abcam mouseanti-mitochondrial IgG1 antibody (Abcam ab3298) andthen with FITC-conjugated goat anti-mouse secondaryantibody (Sigma-Aldrich F4018) and DAPI, followedphalloidin red staining and imaged with 3D-StructuredIllumination Microscopy (SIM) on a DeltaVisionOMXBlaze microscope as described [87]. Images were proc-essed using Fiji/ImageJ software [88], and Area andPerimeter values were extracted to calculate form factor.Cell morphology was scored as either ‘round’ or ‘elon-gated’ by JCC as part of the mitochondrial quantificationprocess.Holo-tomographic imagingHolo-tomographic video imaging was performed on aNanoLive (Switzerland) 3D Cell Explorer fluo (AXT PtyLtd., Warriewood, NSW) equipped with a NanoLive livecell incubator (AXT Pty Ltd). 1 × 104 cells were seededinto a FluoroDish cell culture dish 35 mm, 23 mm well(World Precision Instruments, FD35) and maintained inphenol red free DMEM medium (Sigma-Aldrich, D1145)supplemented with 10% fetal bovine calf serum (Sigma-Aldrich, F9423), 2 mM glutamine (Sigma-Aldrich,G7513) and 1% penicillin-streptomycin for 48 h. Imme-diately prior to imaging the medium was removed andreplaced with 400 μL of the same medium, followed bytransfer to the live cell incubator chamber of the 3D CellExplorer. Cells were incubated at 37 °C, 5% CO2 and100% humidity for the duration of the time-lapse. Threedimensional holo-tomographic images were capturedevery 20 s for the duration of the time-lapse using theNanolive STEVE software. For File S7 the center planeof each 96 slice stack was exported after capture usingthe built in STEVE export wizard as an .avi movie file.These files were exported at 5 frames per second (100xactual speed) to visualize cellular dynamics.Subcutaneous mouse xenograft tumorsAnimal experiments and welfare conformed to the Aus-tralian National Health and Medical Research Councilguidelines and institutional ethics approval from theAnimal Ethics Committee at the Australian NationalUniversity (E2017/16) and Charles Sturt UniversityAnimal Care and Ethics Committee (A17046). NOD/Shi-SCID/IL-2Rγnull (NSG) mice were bred and sup-plied by the Australian Phenomics Facility at the Austra-lian National University (ANU), an Australian academiccenter funded by the Australian Government’s NationalCollaborative Research Infrastructure Strategy (NCRS)program, and contributions from ANU. Experimentswere performed in the John Curtin School of MedicalResearch satellite animal facility. Animals were housedin individually ventilated cages with access to food andwater ad libitum, enriched with corn cob bedding and amaximum of 5 mice per cage. Animals were injected infour cohorts randomized by weight, and the number ofmice per cohort was based on availability. Mice wereinjected with independently generated cell lines MP,WT1-WT3, DM1-DM3, TM1-TM-3, and Y180F. Cellswere expanded in culture for a maximum of 2 weeks be-fore injection. Cells were trypsinized, pelleted, washedwith PBS and stored on ice. Cell count was determinedThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 19 of 24using a hemocytometer and trypan blue. Two millioncells were resuspended in 100 μL 50:50 Matrigel:PBS forsub-cutaneous injection into the left flank of femaleNSG mice (8–12 weeks) via a 27-gauge needle. The totalnumber of mice injected per cell line was based on thereplicate coefficient of variation (CV) for solid tumourmodels evaluating disease burden = 10%, CV being themeasure of spread that describes the amount of variabil-ity relative to the mean. A t-test-based statistical analysisof simulated differences in mean tumour progressionshows that a minimal number of replicates of n = 10 isnecessary to identify a greater than or equal to 20%change between replicate groups at a pre-determinedstatistical significance level of 5%. However this was re-fined in order to reduce animal usage where possible ac-cording to the NHMRC Australian code for the care anduse of animals for scientific purposes. For cohort 1 fouranimals each were injected with MP, WT1, DM1, andTM1. For cohort 2 four animals each were injected withWT2, DM2 and TM2, and one animal with MP. For co-hort 3 four animals each were injected with WT3, DM3and TM3, one animal with MP. For cohort 4 five ani-mals were injected with Y180F cells, and one animalwith each of MP, WT1, DM1, and TM1. Injection of eachcohort commenced mid-morning. Mice were monitoreddaily and once the tumor was palpable (~ 2.5–4 weeks),tumor growth was measured 3 times a week using calipersuntil at least one tumor in the cohort reached 1 cm3 involume. Once any animal in any group in the cohortreached the maximum size of tumor all animals wereeuthanized for tissue recovery by cervical dislocation andtumors were harvested, weighed, photographed, and fixedwith formalin.Statistical analysesUnless specified otherwise, statistical analysis wasperformed using the SPSS package (IBM). Results forWT, DM and TM cells represent equal numbers ofthree independently derived cell lines for eachPGRMC1 condition (e.g. n = 6 for WT includes 2xWT1, 2xWT2, and 2xWT3). For boxplot data depic-tions, whiskers represent quartiles 1 and 4 (withmaximum and minimum values). Boxes representquartiles 2 and 3, separated by the median, as gener-ated by SPSS analyze data function. Datasets con-forming with normal distribution were analyzed byANOVA and post-hoc Bonferroni or Tukey HSD test(equal variance) or post-hoc Dunnett’s T3 Test (un-equal variance). Statistical differences between diver-gent treatments of different cell lines were calculatedusing two way ANOVA and post-hoc pairwise com-parisons. For non-parametric data sets Kruskal-Wallis or Kolmogorov-Smirnov tests were performed,as indicated in relevant figure legends.Supplementary informationSupplementary information accompanies this paper at https://doi.org/10.1186/s12860-020-00256-3.Additional file 1: Figure S1. Detailed views of selected pathwaysidentified by WebGestalt analyses. Related to Fig. 3. All panels areadapted from File S6. Heat map colors follow Fig. 3. (A) Principalcomponent (PC) analysis of SWATH-MS proteomics results showing distri-bution of PC1 and PC2. PC1 corresponded to pathways associated withribosomes and translation, while PC2 corresponded to pathways associ-ated with mRNA splicing processing (see File S3). (B) Proteins associatedwith PI3K/AKT activity (WebGestalt Database: PC, DB_ID:1648, “Class I PI3Ksignaling events mediated by AKT”) are less abundant in TM cells. (C) F1/F0 ATPase subunits elevated in WT and TM cells. (D) Abundances of pro-teins associated with protein folding and microtubule function are alteredby PGRMC1 phosphorylation status. Proteins detected in any of the fol-lowing WebGestalt pathways or functions (1–4) or a manual search (5)are mapped against their expression profiles. 1) cellular componentchaperonin-containing T-complex GO:0005832. 2) PC pathwayChaperonin-mediated protein folding DB_ID:710. 3) cellular componentmicrotubule GO:0005874. 4) PC pathway Protein folding DB_ID:712. 5)Description from the list of 243 proteins (File S6) contains keywords“tubulin” or “microtubule” (manual search) (Adapted from File S6). 2-tailedt-test p-values for all sample comparisons are available in File S4. (E) Pro-teins associated with nuclear import/export that are elevated in DM cells.(F) Antigen processing and presentation enzymes are affected byPGRMC1 phosphorylation status. Manual additions to KEGG pathwayID:04612 “Antigen processing and presentation" (no yellow shading: fromFile S6 and File S5) are indicated with yellow highlighting. Figure S2.Highest and lowest differentially abundant proteins. Panels show the sixmost (+) and least (−) abundant proteins for each cell type that were sig-nificantly differentially abundant between cell types. Related to Fig. 3.Identical colored symbols depict the same protein in different cell types,where circles represent high abundance and triangles represent lowabundance. E.g. Mitochondrial import receptor subunit TOM40 (O96008),CDP-diacylglycerol-inositol 3-phosphatidyltransferase (CDIPT, O14735)(phosphatidylinositol synthesis) and transcriptional coactivator PSIP1(O75475) are more abundant in WT and TM. Aldehyde dehydrogenase1A3 (P47895), APOC3 (P02656) and APOA1 (P02647) are less abundant inDM and TM, whereas the receptor tyrosine kinase ephrin type-A receptor2 (EPHA2, P29317) is among the lowest abundance differential proteins inWT and TM. Protein abundances (measured ion intensities) for all proteinsare available in File S2. Figure S3. Proteins associated with predicted Es-trogen Receptor Related 1 (ERR1) transcription factor activity are moreabundant in WT and TM cells. Related to Fig. 3 and File S6. Data are avail-able via ProteomeXchange with identifier PXD014789. (A) Proteins quanti-fied by SWATH-MS and predicted by WebGestalt Transcription Factortarget enrichment analysis (WebGestalt Database: Transcription Target,Name: hsa_TGACCTY_V$ERR1_Q2, ID:DB_ID:2414) to be dependent uponERR1 transcription (Adapted from File S6). Heat map colors follow Fig. 3.(B) Western blot of shRNA attenuation of ERR1 (top panel) in WT cells,compared to scramble shRNA control. Lower panel: beta actin loadingcontrol. (C) ERR1 attenuation by shRNA induces morphological changesto WT cells. (D) Ratio of rounded to other cells as scored for replicate im-ages from each of the cell lines from B. p < 0.0005, 1 tailed t-test. (E)SWATH MS quantification of proteins significantly differentially abundantbetween cells expressing scramble shRNA or anti-ERR1 shRNA (accordingto the selection criteria of File S4). Fifty-seven proteins became moreabundant after shRNA depletion of ERR1, and 19 proteins became lessabundant. In this panel green depicts lower abundance and red depictshigher abundance. Brown depicts similar abundance. (F) Proteins presentin (E) which are also present in the 243 protein list of File S6. Heat mapcolors follow Fig. 3. The left panel shows expression in the original resultof File S6. The right panel shows expression in the presence of scrambleshRNA or anti-ESSR1 shRNA. All proteins from Fig. S3A except Q9BPW8were detected in the shRNA experiment. Double headed arrows indicateproteins which differ in expression tendency between WT (left) and WT-scramble shRNA cells (right). This may be caused by the puromycin selec-tion of both sh-scr and shERR1 cells but not the parental WT cells, how-ever this requires further investigation. Proteins with Uniprot IDThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 20 of 24highlighted by asterisk (bold red) are those both originally predicted byWebGestalt to be associated with ERR1 transcription factor (A), and whichexhibited significantly altered abundance after shRNA attenuation ofERR1 protein. Figure S4. The DM migratory phenotype is dependentupon Vinculin. Related to Fig. 3. (A) Proteins associated with actin-myosincontraction (WebGestalt Database: PC, DB_ID:144, “Smooth Muscle Con-traction”) are more abundant in DM cells. Heat map colors follow Fig. 3.(B) Western blot of Vinculin (VCL) and actin levels in cells expressingscramble control (shScr) or anti Vinclulin (shVNC) shRNA. (C) Scratch as-says as indicated were performed following the methods of Fig. 1e-f. (D)Boxplot showing results of 12 replicates of (C). ANOVA gave F = 78.7, p <0.00001. The table shows post-hoc Dunnet’s T3 p-values. Figure S5.PGRMC1 phosphorylation status affects mitochondrial function. Relatedto Fig. 5. (A) Representative flow cytometry results of cells labeled withNpFR2. The percentage of the cell population to the right of the dashedreference line (interval labeled “B”, marked by the dotted line) is quanti-fied for each measurement. White arrows indicate more oxidized NpFR2fluorescence in WT and TM cells. (B) Boxplots of the percentage of cellsexhibiting > 10 fluorescent intensity units to the right of the referenceline in (A). n = 6 for each cell type, being 6 replicates of MP cells, or dupli-cate measurements of each of 3 independent lines 1–3 (n = 3 × 2 = 6) forWT, DM and TM cells. White arrows indicate the same differences as in(A). There was significant difference between the means (Kruskal-WallisTest p < 0.0001). Independent sample median tests revealed that all me-dians were significantly different from one another (p < 0.001). (C) Repre-sentative flow cytometry results of cells labeled with MitoTracker. Therespective panels depict fluorescence intensity (x axis) plotted against ei-ther side scatter (left/black plots) or cell number (right/red plots), with orwithout the addition of MitoTracker as indicated. MitoTracker affinity formitochondria is increased with higher mitochondrial membrane potential(Δψm), such that respective left MitoTracker-treated populations/peaksrepresent low Δψm populations, and respective right MitoTracker-treatedpopulations/peaks represent high ψm populations. Vertical dotted linesfor Mitotracker-treated cells depict the MFI for MP cell low and highΔψm populations. The black arrow for DM cells with Mitotracker high-lights reduced MFI of both the low and high Δψm DM cell populations.The white arrow highlights the reduced proportion of DM cells in thelow Δψm peak. The same arrows are depicted for replicates in (D) and (E)respectively. (D) Median Fluorescent Intensity (MFI) for each of the cellpopulations observed from (C). n = 6 measurements per cell type, as per(B). Error bars depict standard deviation. The low and high Δψm popula-tions correspond to the respective populations from (C). For low Δψmcells the means were significantly different between cell types (ANOVA,post-hoc Dunnet’s T3, p < 0.003). For high Δψm cells WT vs. TM (p < 0.05)and other comparisons (p < 0.004) were significantly different (Kruskall-Wallis). Tables D(L) and D(H) show pairwise comparison p-values for lowΔψm (F(L)) and high Δψm (F(H)) cells. (E) The percentage of cells in eachof the low and high Δψm populations from (C) and (D), from the samedata. Standard deviation error bars for the respective high Δψm (uppererror bar) and low Δψm (lower error bar) cell populations are given foreach cell type. Color-coding and all other terminology follows (C). Forlow Δψm cells, Shapiro-Wilk’s test showed normal distributions, andLevene’s test revealed non-homogenous variance (p = 0.04). Table Eshows p-values for ANOVA/post-hoc Dunnet’s T3 test pairwise compari-sons. (F, G) Area under curve (AUC) plots for basal and FCCP curves from(Fig. 5c) showing one way ANOVA/post-hoc Bonferroni p-values betweensamples in the respective indicated regions.Additional file 2: File S1. A zip archive containing time lapse mp4movies of migrating cells in scratch assays. Related to Fig. 1. Images weretaken at 10 min intervals over 36 h and are replayed over 65 s at 9 framesper second (2000x real time). The presented 32.4 MB MP4 files weregenerated from the original 666 MB (6000x real time) .avi files (length 21s, frame width x height 1024 × 1024, Data rate/bit rate 287,462 kbps,frame rate 100 fps) by processing in Adobe Premiere Pro 2017 using thefollowing settings: Preset “custom”; Width: 1024; Height 1024; AspectSquare Pixels (1.0); Field Order: progressive; Profile: high; Target Bitrate[Mbps] 4.08; Maximum Bitrate [Mbps] 4.39; 9 fps; No Audio; Output asstandard .mp4; Time Interpolation: Frame Sampling; Stream Compatibility:Standard; Variable bit rate, single pass; Number of frames 00:01:04:09. Thezip archive contains 4 files with filenames representing cell type and dateof measurement (yyyymmdd). (A) A_MP_20170807.mp4. (B)B_WT_20170809.mp4. (C) C_DM_20170815.mp4. (D)D_TM_20170817.mp4.Additional file 3: File S2. An Excel file showing experimental design,normalized ion intensities for 1330 proteins identified by SWATH-MS pro-teomics, and six pairwise comparisons between the 4 sample types [1)MP v. WT, 2) MP v. DM, 3) MP v. TM, 4) WT v. DM, 5) WT v. TM, and 6) DMv. TM]. Related to Fig. 3. The first tab contains a detailed descriptive le-gend. Data are available via ProteomeXchange with identifier PXD014716.Additional file 4: File S3. Principal component analysis results forpathways associated with SWATH-MS proteomics results.Additional file 5: File S4. An excel file containing proteomics resultsfor 243 proteins which fulfil stringency criteria of t-test p-value of lessthan 0.05, and a fold change greater than 1.5 by both the protein andpeptide approaches from File S2. Related to Fig. 3. Column B shows “red”(more abundant in comparative sample 1) and “blue” (less abundant insample 1) significantly differential proteins for each pair wise comparisonwhich were later used for “red” and ‘blue” WebGestalt pathways enrich-ment analysis (File S5). Comparisons follow File S2.Additional file 6: File S5. WebGestalt pathway mapping excel resultsfile for red and blue proteins from File S4. Related to Fig. 3. (A)WebGestalt Features (GO, KEGG Pathways, Wikipathways, PathwayCommons, Transcription Factors) significant at the adjP≤0.001 levelbetween any 2 comparisons. (B) Features from A, viewed at the adjP<0.05 level for each red and blue comparison, and showing adjusted p-values (adjP) for each comparison where adjP< 0.05 (or adjP< 0.1 for twocomparisons as indicated by paler coloring). Blue means that proteinsassociated with a given the feature were less abundant in that cell line,with red indicating higher abundance. Because separate WebGestaltanalyses were performed for red and blue lists of proteins from File S4,some features were significant for both red and blue. In that case thecolor and adjP for the most significant analysis is given, with the othercell being colored black.Additional file 7: File S6. An excel file with heat map protein IDs andpathways for red vs. blue pathways adjP< 0.001. Related to Fig. 3. This fileis derived from the results of File S5, and is the source file for Fig. 3.Proteins suggested by clustering by inferred models of evolution (CLIME)analysis to co-evolve with PGRMC1 with log likelihood ratio greater than12 [21] are present in the list, marked yellow for mitochondriallocalization (WebGestalt GO:0005739) or green for cytoplasmic (P00387).Additional file 8: File S7. A zip archive containing time lapse Holo-tomographic video.avi files of cells. These images are based upon differ-ences in refractive index [53], and are provided for the dynamicvisualization of mitochondria. Prominent visible features include smallwhite lipid droplets and cholesterol-rich mitochondria [21], as well as nu-clear membrane and nucleoli. The previously described MIA PaCa-2 cellbleb-like protrusions [34] are apparent as highly dynamic rearrangementsof the cytoplasmic membrane, which may contribute to intercellularcommunication. (A) MP cells. (B) WT cells. (C) DM cells. (D) TM cells.AbbreviationsadjP: Benjamini-Hochberg adjusted p-value; avFF: average mitochondrialform factor; CK2: Casein kinase 2; DM: Hemagglutinin-tagged PGRMC1 S57A/S181 double mutant; EMT: Epithelial-mesenchymal transition; ERR1: Estrogenreceptor related 1; FF: Mitochondrial form factor; hPSCs: human pluripotentstem cells; IDA: Information dependent acquisition; LDLR: Low densitylipoprotein receptor; MAT: Mesenchymal-amoeboid transition; MP: MIA PaCa-2 pancreatic cancer; MS/MS: Tandem mass spectrometry; MTT: 3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide; NpFR2: Naphthalimide-flavin redox sensor 2; P4: Progesterone; PAQR7: Progestin and adipoQreceptor 7; PGRMC1: Progesterone Receptor Membrane Component 1;PGRMC2: Progesterone Receptor Membrane Component 2; ROCK: Rhokinase; ROCKI: Rho kinase inhibitor; RPPA: Reverse Phase Protein Array;S2R: Sigma-2 receptor; SH2: Src homology 2; SH3: Src homology 3; SWATH-MS: Sequential Window Acquisition of all Theoretical Mass Spectrometry;TM: Hemagglutinin-tagged PGRMC1 S57A/Y180F/S181 triple mutant; TRiC: T-complex protein-1 ring complex; WT: Hemagglutinin-tagged PGRMC1 wildtype; Δψm: Mitochondrial membrane potentialThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 21 of 24AcknowledgementsThe publication hiatus of more than a decade on this project has beencaused by the Australian competitive grants system failing to support asingle funding application from M.A.C. after his relocation with the PGRMC1project from Germany to a regional Australian university. We thank PadraigStrappe for the gift of lentiviral helper packaging plasmids and HEK293 cells.MIA PaCa-2 (MP) pancreatic cancer cells were obtained from Dr. Patsy Soon,Kolling Institute of Medical Research, Sydney. We thank Jean Yang for assist-ance in establishing the project.Authors’ contributionsAll authors have read and approved the manuscript. Conceptualization,M.A.C.; Methodology, M.A.C, P.P.A, A.V.O., M.P., E.J.N, S.J.K., K.M.H, R.D.H., andM.P.M; Software, D.P. (proteomics); Formal Analysis, P.P.A., D.P., M.J., E.P.M.A.C., and M.J.; Investigation, B.M.T., P.P.A., S.L.T., A.V.O., A.K., I.S., T.Z., C.P.C.,D.P., L.T., M.P., P.P., and J.C.C.; Resources, M.P.M., M.P., E.J.N., H.N., and T.F.;Writing – Original Draft, M.A.C.; Writing – Review & Editing, M.P.M., M.A.C.,P.P.A., E.M.G., C.P.C., B.M.T., M.P., E.M.G, and T.L.R.; Visualization, M.A.C.;Supervision, M.A.C, J.A.J., E.M.G., E.J.N., T.L.R, and M.P.M.; ProjectAdministration, M.A.C.; Funding Acquisition, M.A.C., J.A.J., B.M.T., E.M.G., M.P.M.,T.L.R., and E.J.N.FundingThis work has received no direct Australian competitive grant support sinceM.A.C.’s relocation to the country in 2008. The present results have beencompiled largely due to the generosity of collaborating authors, and by thePhD stipends of B.M.T. and P.P.A. Research was primarily supported byCharles Sturt University internal funds to M.A.C and J.A.J, and bycollaborating labs. Open access publication fees were jointly supported byCSU’s Faculty of Science, School of Biomedical Sciences, and Research Office.B.M.T. was supported by a PhD scholarship from the Ministry of HigherEducation and Research, through the University of Wasit, Iraq. A.K.acknowledges the University of Sydney for a World Scholars Scholarship.E.M.G. acknowledges partial support of Australian Research Council awardCE140100003. E.J.N. acknowledges the support of the Australian ResearchCouncil (DE120102687) and the Ramaciotti Foundation (ES2012/0051). T.L.R issupported by a Cancer Institute New South Wales Future Research LeaderFellowship. M.P. is supported by a NHMRC RD Wright Biomedical CareerDevelopment Fellowship and the Cancer Institute NSW Career DevelopmentFellowship. The funding bodies played no role in the design of the studyand collection, analysis, and interpretation of data and in writing themanuscript.Availability of data and materialsThe mass spectrometry proteomics data have been deposited to theProteomeXchange Consortium via the PRIDE partner repository (http://www.ebi.ac.uk/pride) with the dataset identifiers PXD014716 (Fig. 3) andPXD014789 (Fig. S3). Further data are provided in supporting informationfiles available from the journal web page.Ethics approval and consent to participateMouse experiments were approved by The Australian National UniversityAnimal Experimentation Ethics Committee ethics protocol ANU A2017/16,and by Charles Sturt University Animal Care and Ethics Committee ethicsprotocol CSU A17046. ANU ethics approval and the provision agreementwith the Australian Phenomics Facility at the Australian National University(ANU) (the provider of animals), included permission to use the animals inthis study.Consent for publicationNot applicable.Competing interestsM.A.C. is scientific advisor to and minor shareholder of CognitionTherapeutics, a company developing sigma-2 receptor ligands against Alzhei-mer’s disease. This work was performed independently of and without inputfrom the company. The authors declare that they have no other competinginterests.Author details1School of Biomedical Sciences, Charles Sturt University, Wagga Wagga, NSW2650, Australia. 2Department of Biology, College of Science, University ofWasit, Wasit 00964, Iraq. 3Present address: Faculty of Pharmaceutical Sciences,The University of British Columbia, Vancouver V6T 1Z4, Canada. 4School ofChemistry, University of Sydney, Sydney, NSW 2006, Australia. 5Presentaddress: School of Medical Sciences, University of Sydney, Sydney, NSW 2006,Australia. 6Life Sciences and Health, Faculty of Science, Charles SturtUniversity, Wagga Wagga, NSW 2650, Australia. 7The Kinghorn Cancer Centre,Cancer Division, Garvan Institute of Medical Research, 384 Victoria St,Darlinghurst, Sydney, NSW 2010, Australia. 8St Vincent’s Clinical School,Faculty of Medicine, University of NSW, Darlinghurst 2010, NSW, Australia.9Department of Biochemistry and Molecular Biology, University ofMelbourne, Parkville, Victoria 3010, Australia. 10ACRF Department of CancerBiology and Therapeutics, The John Curtin School of Medical Research, TheAustralian National University, ACT, Canberra 2601, Australia. 11AustralianProteome Analysis Facility, Macquarie University, Sydney, NSW 2109, Australia.12Department of Gynecology and Obstetrics, University Women’s Hospital ofDusseldorf, 40225 Dusseldorf, Germany. 13NMI TT Pharmaservices, ProteinProfiling, 72770 Reutlingen, Germany. 14Particles-Biology InteractionsLaboratory, Department of Materials Meet Life, Swiss Federal Laboratories forMaterials Science & Technology (Empa), Lerchenfeldstrasse 5, CH-9014 StGallen, Switzerland. 15The ithree institute, University of Technology Sydney,Ultimo, NSW 2007, Australia. 16Present address: GE Healthcare Life Sciences,Issaquah, WA 98027, USA. 17Department of Mathematics and Statistics, LaTrobe University, Melbourne, VIC 3086, Australia. 18Ingham Institute forApplied Medical Research, Liverpool, NSW 2170, Australia. 19School ofMedicine, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751,Australia. 20South Western Sydney Clinical School, Faculty of Medicine, TheUniversity of New South Wales, Liverpool, NSW 2170, Australia. 21School ofMedicine and University of Queensland Centre for Clinical Research, Herston,QLD 4006, Australia. 22AXT Pty Ltd, Warriewood (Sydney), NSW 2102,Australia. 23Division of Research, Peter MacCallum Cancer Centre, Melbourne,Victoria 3000, Australia. 24Sir Peter MacCallum Department of Oncology,University of Melbourne, Parkville, Victoria 3010, Australia. 25Department ofBiochemistry and Molecular Biology, Monash University, Clayton, Victoria3168, Australia. 26School of Mathematics and Statistics, University of Sydney,Sydney, NSW 2006, Australia. 27Present address: The Kolling Institute, TheUniversity of Sydney, St Leonards (Sydney), NSW 2064, Australia. 28ARCCentre of Excellence for Nanoscale BioPhotonics, Macquarie University,Sydney, NSW 2109, Australia. 29Present address: The Graduate School ofBiomedical Engineering, University of New South Wales, Sydney, Kensington,NSW 2052, Australia. 30Graham Centre for Agricultural Innovation, CharlesSturt University, Boorooma Street, Wagga Wagga, NSW 2678, Australia.31School of Agricultural and Wine Sciences, Charles Sturt University,Boorooma Street, Wagga Wagga, NSW 2678, Australia.Received: 16 January 2020 Accepted: 4 March 2020References1. Riad A, Zeng C, Weng CC, Winters H, Xu K, Makvandi M, Metz T, Carlin S,Mach RH. Sigma-2 receptor/TMEM97 and PGRMC-1 increase the rate ofinternalization of LDL by LDL receptor through the formation of a ternarycomplex. Sci Rep. 2018;8(1):16845.2. Cahill MA, Jazayeri JA, Catalano SM, Toyokuni S, Kovacevic Z, Richardson DR.The emerging role of progesterone receptor membrane component 1(PGRMC1) in cancer biology. Biochim Biophys Acta. 2016;1866(2):339–49.3. Ryu CS, Klein K, Zanger UM. Membrane associated progesterone receptors:promiscuous proteins with pleiotropic functions - focus on interactions withcytochromes P450. Front Pharmacol. 2017;8:159.4. Peluso JJ, Griffin D, Liu X, Horne M. Progesterone receptor membranecomponent-1 (PGRMC1) and PGRMC-2 interact to suppress entry into thecell cycle in spontaneously immortalized rat granulosa cells. Biol Reprod.2014;91(5):104.5. Sueldo C, Liu X, Peluso JJ. Progestin and AdipoQ receptor 7, progesteronemembrane receptor component 1 (PGRMC1), and PGRMC2 and their role inregulating Progesterone's ability to suppress human Granulosa/luteal cellsfrom entering into the cell cycle. Biol Reprod. 2015;93(3):63.6. Luciano AM, Lodde V, Franciosi F, Ceciliani F, Peluso JJ. Progesteronereceptor membrane component 1 expression and putative function inThejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 22 of 24bovine oocyte maturation, fertilization, and early embryonic development.Reproduction (Cambridge, England). 2010;140(5):663–72.7. Luciano AM, Peluso JJ. PGRMC1 and the faithful progression throughmitosis and meiosis. Cell cycle (Georgetown, Tex). 2016;15(17):2239–40.8. Terzaghi L, Luciano AM, Dall'Acqua PC, Modina SC, Peluso JJ, Lodde V.PGRMC1 localization and putative function in the nucleolus of bovinegranulosa cells and oocytes. Darlinghurst: Reproduction (Cambridge,England). 2018;155(3):273–82.9. Terzaghi L, Tessaro I, Raucci F, Merico V, Mazzini G, Garagna S, Zuccotti M,Franciosi F, Lodde V. PGRMC1 participates in late events of bovinegranulosa cells mitosis and oocyte meiosis. Cell cycle (Georgetown, Tex).2016;15(15):2019–32.10. Rohe HJ, Ahmed IS, Twist KE, Craven RJ. PGRMC1 (progesterone receptormembrane component 1): a targetable protein with multiple functions insteroid signaling, P450 activation and drug binding. Pharmacol Ther. 2009;121(1):14–9.11. Ahmed IS, Rohe HJ, Twist KE, Mattingly MN, Craven RJ. Progesteronereceptor membrane component 1 (Pgrmc1): a heme-1 domain protein thatpromotes tumorigenesis and is inhibited by a small molecule. J PharmacolExp Ther. 2010;333(2):564–73.12. Losel RM, Besong D, Peluso JJ, Wehling M. Progesterone receptormembrane component 1--many tasks for a versatile protein. Steroids. 2008;73(9–10):929–34.13. Ruan X, Zhang Y, Mueck AO, Willibald M, Seeger H, Fehm T, Brucker S,Neubauer H. Increased expression of progesterone receptor membranecomponent 1 is associated with aggressive phenotype and poor prognosisin ER-positive and negative breast cancer. Menopause (New York, NY). 2017;24(2):203–9.14. Willibald M, Wurster I, Meisner C, Vogel U, Seeger H, Mueck AO, Fehm T,Neubauer H. High Level of Progesteron Receptor Membrane Component 1(PGRMC 1) in Tissue of Breast Cancer Patients is Associated with WorseResponse to Anthracycline-Based Neoadjuvant Therapy. Horm Metab Res.2017;49(8):595–603.15. Shih CC, Chou HC, Chen YJ, Kuo WH, Chan CH, Lin YC, Liao EC, Chang SJ,Chan HL. Role of PGRMC1 in cell physiology of cervical cancer. Life Sci.2019;231(19):116541.16. Cahill MA. Progesterone receptor membrane component 1: an integrativereview. J Steroid Biochem Mol Biol. 2007;105(1–5):16–36.17. Peluso JJ, Pappalardo A, Losel R, Wehling M. Progesterone membranereceptor component 1 expression in the immature rat ovary and its role inmediating progesterone's antiapoptotic action. Endocrinology. 2006;147(6):3133–40.18. Franchin C, Borgo C, Cesaro L, Zaramella S, Vilardell J, Salvi M, Arrigoni G,Pinna LA. Re-evaluation of protein kinase CK2 pleiotropy: new insightsprovided by a phosphoproteomics analysis of CK2 knockout cells. Cell MolLife Sci. 2018;75(11):2011–26.19. Cahill MA, Jazayeri JA, Kovacevic Z, Richardson DR. PGRMC1 regulation byphosphorylation: potential new insights in controlling biological activity.Oncotarget. 2016;7(32):50822–7.20. Cahill MA. The evolutionary appearance of signaling motifs in PGRMC1.Biosci Trends. 2017;11(2):179–92.21. Cahill MA, Medlock AE. Thoughts on interactions between PGRMC1 anddiverse attested and potential hydrophobic ligands. J Steroid Biochem MolBiol. 2017;171:11–33.22. Hehenberger E, Eitel M, Fortunato SAV, Miller DJ, Keeling PJ, Cahill MA. Earlyeukaryotic origins and metazoan elaboration of MAPR family proteins. bioRxiv.2019: bioRxiv preprint (not peer-reviewed). https://doi.org/10.1101/737684.23. Etchevers HC, Dupin E, Le Douarin NM. The diverse neural crest: fromembryology to human pathology. Development (Cambridge, England).2019;146(5):dev169821.24. Neubauer H, Clare SE, Wozny W, Schwall GP, Poznanovic S, Stegmann W,Vogel U, Sotlar K, Wallwiener D, Kurek R, et al. Breast cancer proteomicsreveals correlation between estrogen receptor status and differentialphosphorylation of PGRMC1. Breast Cancer Res. 2008;10(5):R85.25. Sabbir MG. Progesterone induced Warburg effect in HEK293 cells isassociated with post-translational modifications and proteasomaldegradation of progesterone receptor membrane component 1. J SteroidBiochem Mol Biol. 2019;191:105376.26. Gras F, Brunmair B, Quarre L, Szocs Z, Waldhausl W, Furnsinn C.Progesterone impairs cell respiration and suppresses a compensatoryincrease in glucose transport in isolated rat skeletal muscle: a non-genomicmechanism contributing to metabolic adaptation to late pregnancy?Diabetologia. 2007;50(12):2544–52.27. Hampton KK, Anderson K, Frazier H, Thibault O, Craven RJ. Insulin receptorplasma membrane levels increased by the progesterone receptormembrane component 1. Mol Pharmacol. 2018;94(1):665–73.28. Wu H, Wu ZG, Shi WJ, Gao H, Wu HH, Bian F, Jia PP, Hou YN. Effects ofprogesterone on glucose uptake in neurons of Alzheimer's disease animalsand cell models. Life Sci. 2019;238:116979.29. Yunis AA, Arimura GK, Russin DJ. Human pancreatic carcinoma (MIA PaCa-2)in continuous culture: sensitivity to asparaginase. Int J Cancer J Int DuCancer. 1977;19(1):128–35.30. Han EK, McGonigal T, Butler C, Giranda VL, Luo Y. Characterization of Aktoverexpression in MiaPaCa-2 cells: prohibitin is an Akt substrate bothin vitro and in cells. Anticancer Res. 2008;28(2A):957–63.31. Duong HQ, Hong YB, Kim JS, Lee HS, Yi YW, Kim YJ, Wang A, Zhao W, ChoCH, Seong YS, et al. Inhibition of checkpoint kinase 2 (CHK2) enhancessensitivity of pancreatic adenocarcinoma cells to gemcitabine. J Cell MolMed. 2013;17(10):1261–70.32. Iwagami Y, Eguchi H, Nagano H, Akita H, Hama N, Wada H, Kawamoto K,Kobayashi S, Tomokuni A, Tomimaru Y, et al. miR-320c regulatesgemcitabine-resistance in pancreatic cancer via SMARCC1. Br J Cancer. 2013;109(2):502–11.33. Gosnell ME, Anwer AG, Mahbub SB, Menon Perinchery S, Inglis DW,Adhikary PP, Jazayeri JA, Cahill MA, Saad S, Pollock CA, et al. Quantitativenon-invasive cell characterisation and discrimination based on multispectralautofluorescence features. Sci Rep. 2016;6:23453.34. Gradiz R, Silva HC, Carvalho L, Botelho MF, Mota-Pinto A. MIA PaCa-2 andPANC-1 - pancreas ductal adenocarcinoma cell lines with neuroendocrinedifferentiation and somatostatin receptors. Sci Rep. 2016;6:21648.35. Fujita M, Otsuka Y, Yamada S, Iwakawa M, Imai T. X-ray irradiation and rho-kinase inhibitor additively induce invasiveness of the cells of the pancreaticcancer line, MIAPaCa-2, which exhibits mesenchymal and amoeboidmotility. Cancer Sci. 2011;102(4):792–8.36. Thejer BM, Adhikary PP, Teakel SL, Fang J, Weston PA, Gurusinghe S, AnwerAG, Gosnell M, Jazayeri JA, Ludescher M, et al. PGRMC1 effects onmetabolism, genomic mutation and CpG methylation imply crucial roles inanimal biology and disease. ResearchSquare preprint (not peer-reviewed);2020. https://dx.doi.org/10.21203/rs.2.20008/v2.37. Suchanek M, Radzikowska A, Thiele C. Photo-leucine and photo-methionineallow identification of protein-protein interactions in living cells. NatMethods. 2005;2(4):261–7.38. Matsuoka T, Yashiro M. Rho/ROCK signaling in motility and metastasis ofgastric cancer. World J Gastroenterol. 2014;20(38):13756–66.39. Cha D, O'Brien P, O'Toole EA, Woodley DT, Hudson LG. Enhancedmodulation of keratinocyte motility by transforming growth factor-alpha(TGF-alpha) relative to epidermal growth factor (EGF). J Invest Dermatol.1996;106(4):590–7.40. Bordeleau ME, Aucagne R, Chagraoui J, Girard S, Mayotte N, Bonneil E,Thibault P, Pabst C, Bergeron A, Barabe F, et al. UBAP2L is a novel BMI1-interacting protein essential for hematopoietic stem cell activity. Blood.2014;124(15):2362–9.41. Harada H, Becknell B, Wilm M, Mann M, Huang LJ, Taylor SS, Scott JD,Korsmeyer SJ. Phosphorylation and inactivation of BAD by mitochondria-anchored protein kinase a. Mol Cell. 1999;3(4):413–22.42. Hayakawa J, Ohmichi M, Kurachi H, Kanda Y, Hisamoto K, Nishio Y, Adachi K,Tasaka K, Kanzaki T, Murata Y. Inhibition of BAD phosphorylation either atserine 112 via extracellular signal-regulated protein kinase cascade or atserine 136 via Akt cascade sensitizes human ovarian cancer cells to cisplatin.Cancer Res. 2000;60(21):5988–94.43. Cross DA, Alessi DR, Cohen P, Andjelkovich M, Hemmings BA. Inhibition ofglycogen synthase kinase-3 by insulin mediated by protein kinase B. Nature.1995;378(6559):785–9.44. Antonietti P, Linder B, Hehlgans S, Mildenberger IC, Burger MC, Fulda S,Steinbach JP, Gessler F, Rodel F, Mittelbronn M, et al. Interference withthe HSF1/HSP70/BAG3 pathway primes Glioma cells to matrixdetachment and BH3 mimetic-induced apoptosis. Mol Cancer Ther.2017;16(1):156–68.45. Joshi B, Strugnell SS, Goetz JG, Kojic LD, Cox ME, Griffith OL, Chan SK, JonesSJ, Leung SP, Masoudi H, et al. Phosphorylated caveolin-1 regulates rho/ROCK-dependent focal adhesion dynamics and tumor cell migration andinvasion. Cancer Res. 2008;68(20):8210–20.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 23 of 2446. Le Clainche C, Dwivedi SP, Didry D, Carlier MF. Vinculin is a dually regulatedactin filament barbed end-capping and side-binding protein. J Biol Chem.2010;285(30):23420–32.47. Pal I, Rajesh Y, Banik P, Dey G, Dey KK, Bharti R, Naskar D, Chakraborty S,Ghosh SK, Das SK, et al. Prevention of epithelial to mesenchymal transitionin colorectal carcinoma by regulation of the E-cadherin-beta-catenin-vinculin axis. Cancer Lett. 2019;452:254–63.48. Kaur A, Brigden KW, Cashman TF, Fraser ST, New EJ. Mitochondriallytargeted redox probe reveals the variations in oxidative capacity of thehaematopoietic cells. Org Biomol Chem. 2015;13(24):6686–9.49. Perry SW, Norman JP, Barbieri J, Brown EB, Gelbard HA. Mitochondrialmembrane potential probes and the proton gradient: a practical usageguide. BioTechniques. 2011;50(2):98–115.50. Turrens JF. Mitochondrial formation of reactive oxygen species. J Physiol.2003;552(Pt 2):335–44.51. Gosnell ME, Anwer AG, Cassano JC, Sue CM, Goldys EM. Functionalhyperspectral imaging captures subtle details of cell metabolism in olfactoryneurosphere cells, disease-specific models of neurodegenerative disorders.Biochim Biophys Acta. 2016;1863(1):56–63.52. Koopman WJ, Visch HJ, Smeitink JA, Willems PH. Simultaneous quantitativemeasurement and automated analysis of mitochondrial morphology, mass,potential, and motility in living human skin fibroblasts. Cytometry A. 2006;69(1):1–12.53. Ali A, Abouleila Y, Amer S, Furushima R, Emara S, Equis S, Cotte Y, MasujimaT. Quantitative live single-cell mass spectrometry with spatial evaluation bythree-dimensional holographic and tomographic laser microscopy. Anal Sci.2016;32(2):125–7.54. Anesti V, Scorrano L. The relationship between mitochondrial shape andfunction and the cytoskeleton. Biochim Biophys Acta. 2006;1757(5–6):692–9.55. Sahai E, Marshall CJ. Differing modes of tumour cell invasion have distinctrequirements for rho/ROCK signalling and extracellular proteolysis. Nat CellBiol. 2003;5(8):711–9.56. Seol HJ, Chang JH, Yamamoto J, Romagnuolo R, Suh Y, Weeks A, Agnihotri S,Smith CA, Rutka JT. Overexpression of CD99 increases the migration andinvasiveness of human malignant Glioma cells. Genes Cancer. 2012;3(9–10):535–49.57. Taddei ML, Giannoni E, Morandi A, Ippolito L, Ramazzotti M, Callari M,Gandellini P, Chiarugi P. Mesenchymal to amoeboid transition is associatedwith stem-like features of melanoma cells. Cell Commun Signal. 2014;12:24.58. Friedl P. Prespecification and plasticity: shifting mechanisms of cellmigration. Curr Opin Cell Biol. 2004;16(1):14–23.59. Parri M, Taddei ML, Bianchini F, Calorini L, Chiarugi P. EphA2 reexpressionprompts invasion of melanoma cells shifting from mesenchymal toamoeboid-like motility style. Cancer Res. 2009;69(5):2072–81.60. Friedl P, Wolf K. Plasticity of cell migration: a multiscale tuning model. J CellBiol. 2010;188(1):11–9.61. Wai T, Langer T. Mitochondrial dynamics and metabolic regulation. TrendsEndocrinol Metab. 2016;27(2):105–17.62. Piel RB 3rd, Shiferaw MT, Vashisht AA, Marcero JR, Praissman JL, Phillips JD,Wohlschlegel JA, Medlock AE. A novel role for progesterone receptormembrane component 1 (PGRMC1): a partner and regulator ofFerrochelatase. Biochemistry. 2016;55(37):5204–17.63. Kim HJ, Choi WJ, Lee CH. Phosphorylation and reorganization of keratinnetworks: implications for carcinogenesis and epithelial Mesenchymaltransition. Biomol Ther. 2015;23(4):301–12.64. Teakel SL, Ludescher M, Thejer BM, Poschmann G, Forwood JK, Neubauer H,Cahill MA. Protein complexes including PGRMC1 and actin-associatedproteins are disrupted by AG-205. Biochem Biophys Res Commun. 2020;524(1):64–69.65. Roh SH, Kasembeli M, Bakthavatsalam D, Chiu W, Tweardy DJ. Contribution of thetype II Chaperonin, TRiC/CCT, to Oncogenesis. Int J Mol Sci. 2015;16(11):26706–20.66. Boudiaf-Benmammar C, Cresteil T, Melki R. The cytosolic chaperonin CCT/TRiC and cancer cell proliferation. PLoS One. 2013;8(4):e60895.67. Guest ST, Kratche ZR, Bollig-Fischer A, Haddad R, Ethier SP. Two members of theTRiC chaperonin complex, CCT2 and TCP1 are essential for survival of breastcancer cells and are linked to driving oncogenes. Exp Cell Res. 2015;332(2):223–35.68. Zhang Y, Wang Y, Wei Y, Wu J, Zhang P, Shen S, Saiyin H, Wumaier R, YangX, Wang C, et al. Molecular chaperone CCT3 supports proper mitoticprogression and cell proliferation in hepatocellular carcinoma cells. CancerLett. 2016;372(1):101–9.69. Galmozzi A, Kok BP, Kim AS, Montenegro-Burke JR, Lee JY, Spreafico R, MosureS, Albert V, Cintron-Colon R, Godio C, et al. PGRMC2 is an intracellular haemchaperone critical for adipocyte function. Nature. 2019;576(7785):138–42.70. Lee SR, Kwon SW, Kaya P, Lee YH, Lee JG, Kim G, Lee GS, Baek IJ, Hong EJ.Loss of progesterone receptor membrane component 1 promotes hepaticsteatosis via the induced de novo lipogenesis. Sci Rep. 2018;8(1):15711.71. Parker CG, Galmozzi A, Wang Y, Correia BE, Sasaki K, Joslyn CM, Kim AS,Cavallaro CL, Lawrence RM, Johnson SR, et al. Ligand and target discoveryby fragment-based screening in human cells. Cell. 2017;168(3):527–41 e529.72. Hand RA, Craven RJ. Hpr6.6 protein mediates cell death from oxidative damage inMCF-7 human breast cancer cells. J Cell Biochem. 2003;90(3):534–47.73. He Y, Zhang P, Zhang D, Xia Z, Wang X, Deng S, Li H, Zhu H, Xu N, Liang S.Combined assessment of low PGRMC1/positive ATP1A1 levels has enhancedprognostic value for renal cell carcinoma. Oncol Rep. 2018;40(3):1467–76.74. Zhu X, Ji M, Han Y, Guo Y, Zhu W, Gao F, Yang X, Zhang C. PGRMC1-dependent autophagy by hyperoside induces apoptosis and sensitizesovarian cancer cells to cisplatin treatment. Int J Oncol. 2017;50(3):835–46.75. Zhu X, Han Y, Fang Z, Wu W, Ji M, Teng F, Zhu W, Yang X, Jia X, Zhang C.Progesterone protects ovarian cancer cells from cisplatin-induced inhibitoryeffects through progesterone receptor membrane component 1/2 as wellas AKT signaling. Oncol Rep. 2013;30(5):2488–94.76. Liu L, Wang J, Zhao L, Nilsen J, McClure K, Wong K, Brinton RD.Progesterone increases rat neural progenitor cell cycle gene expression andproliferation via extracellularly regulated kinase and progesterone receptormembrane components 1 and 2. Endocrinology. 2009;150(7):3186–96.77. Kim JY, Kim SY, Choi HS, Kim MK, Lee HM, Jang YJ, Ryu CJ. Progesterone receptormembrane component 1 suppresses the p53 and Wnt/beta-catenin pathways topromote human pluripotent stem cell self-renewal. Sci Rep. 2018;8(1):3048.78. Riester M, Xu Q, Moreira A, Zheng J, Michor F, Downey RJ. The Warburgeffect: persistence of stem-cell metabolism in cancers as a failure ofdifferentiation. Ann Oncol. 2018;29(1):264–70.79. Quin J, Chan KT, Devlin JR, Cameron DP, Diesch J, Cullinane C, Ahern J, KhotA, Hein N, George AJ, et al. Inhibition of RNA polymerase I transcriptioninitiation by CX-5461 activates non-canonical ATM/ATR signaling.Oncotarget. 2016;7(31):49800–18.80. Gurusinghe S, Young P, Michelsen J, Strappe P. Suppression ofdedifferentiation and hypertrophy in canine chondrocytes through lentiviralvector expression of Sox9 and induced pluripotency stem cell factors.Biotechnol Lett. 2015;37(7):1495–504.81. Gillet LC, Navarro P, Tate S, Rost H, Selevsek N, Reiter L, Bonner R, AebersoldR. Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: a new concept for consistent and accurateproteome analysis. Mol Cell Proteomics. 2012;11(6):O111.016717.82. Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit(WebGestalt): update 2013. Nucleic Acids Res. 2013;41(Web Server issue):W77–83.83. Pawlak M, Schick E, Bopp MA, Schneider MJ, Oroszlan P, Ehrat M. Zeptosens’protein microarrays: a novel high performance microarray platform for lowabundance protein analysis. Proteomics. 2002;2(4):383–93.84. Pirnia F, Pawlak M, Thallinger GG, Gierke B, Templin MF, Kappeler A,Betticher DC, Gloor B, Borner MM. Novel functional profiling approachcombining reverse phase protein microarrays and human 3-D ex vivo tissuecultures: expression of apoptosis-related proteins in human colon cancer.Proteomics. 2009;9(13):3535–48.85. Bader S, Zajac M, Friess T, Ruge E, Rieder N, Gierke B, Heubach Y, Thomas M,Pawlak M. Evaluation of protein profiles from treated Xenograft tumormodels identifies an antibody panel for formalin-fixed and paraffin-embedded (FFPE) tissue analysis by reverse phase protein arrays (RPPA). MolCell Proteomics. 2015;14(10):2775–85.86. Kaistha BP, Lorenz H, Schmidt H, Sipos B, Pawlak M, Gierke B, Kreider R, Lankat-Buttgereit B, Sauer M, Fiedler L, et al. PLAC8 localizes to the inner plasma membraneof pancreatic Cancer cells and regulates cell growth and disease progressionthrough critical cell-cycle regulatory pathways. Cancer Res. 2016;76(1):96–107.87. Strauss MP, Liew AT, Turnbull L, Whitchurch CB, Monahan LG, Harry EJ. 3D-SIM super resolution microscopy reveals a bead-like arrangement for FtsZand the division machinery: implications for triggering cytokinesis. PLoSBiol. 2012;10(9):e1001389.88. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T,Preibisch S, Rueden C, Saalfeld S, Schmid B, et al. Fiji: an open-sourceplatform for biological-image analysis. Nat Methods. 2012;9(7):676–82.Publisher’s NoteSpringer Nature remains neutral with regard to jurisdictional claims inpublished maps and institutional affiliations.Thejer et al. BMC Molecular and Cell Biology (2020) 21:24 Page 24 of 24
Article
eng
Reviewed
Vancouver : University of British Columbia Library
BioMed Central
Attribution 4.0 International (CC BY 4.0)
Faculty
Mitochondria
Migration
Invasion
Metabolism
Cytochrome P450
Mesenchymal amoeboid transition
Proteomics
Tumor biology
PGRMC1 phosphorylation affects cell shape, motility, glycolysis, mitochondrial form and function, and tumor growth
Text